Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AAC
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • AAC Podcast
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Antimicrobial Agents and Chemotherapy
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AAC
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • AAC Podcast
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Epidemiology and Surveillance

Evolution and Genetic Diversity of the k13 Gene Associated with Artemisinin Delayed Parasite Clearance in Plasmodium falciparum

M. Andreína Pacheco, Esha R. Kadakia, Zainab Chaudhary, Douglas J. Perkins, Julia Kelley, Shashidhar Ravishankar, Michael Cranfield, Eldin Talundzic, Venkatachalam Udhayakumar, Ananias A. Escalante
M. Andreína Pacheco
aDepartment of Biology/Institute for Genomics and Evolutionary Medicine (iGEM), Temple University, Philadelphia, Pennsylvania, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Esha R. Kadakia
aDepartment of Biology/Institute for Genomics and Evolutionary Medicine (iGEM), Temple University, Philadelphia, Pennsylvania, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zainab Chaudhary
aDepartment of Biology/Institute for Genomics and Evolutionary Medicine (iGEM), Temple University, Philadelphia, Pennsylvania, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Douglas J. Perkins
bCenter for Global Health, University of New Mexico Health Sciences Center, Albuquerque, New Mexico, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Julia Kelley
cAtlanta Research and Education Foundation, VAMC, Atlanta, Georgia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Shashidhar Ravishankar
dSchool of Biological Sciences, Georgia Institute of Technology, Atlanta, Georgia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael Cranfield
eGorilla Doctors, Wildlife Health Department, University of California, Davis, Davis, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eldin Talundzic
fMalaria Branch, Division of Parasitic Diseases and Malaria, Center for Global Health, Centers for Disease Control and Prevention, Atlanta, Georgia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Venkatachalam Udhayakumar
fMalaria Branch, Division of Parasitic Diseases and Malaria, Center for Global Health, Centers for Disease Control and Prevention, Atlanta, Georgia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ananias A. Escalante
aDepartment of Biology/Institute for Genomics and Evolutionary Medicine (iGEM), Temple University, Philadelphia, Pennsylvania, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/AAC.02550-18
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Mutations in the Plasmodium falciparum k13 (Pfk13) gene are linked to delayed parasite clearance in response to artemisinin-based combination therapies (ACTs) in Southeast Asia. To explore the evolutionary rate and constraints acting on this gene, k13 orthologs from species sharing a recent common ancestor with P. falciparum and Plasmodium vivax were analyzed. These comparative studies were followed by genetic polymorphism analyses within P. falciparum using 982 complete Pfk13 sequences from public databases and new data obtained by next-generation sequencing from African and Haitian isolates. Although k13 orthologs evolve at heterogeneous rates, the gene was conserved across the genus, with only synonymous substitutions being found at residues where mutations linked to the delayed parasite clearance phenotype have been reported. This suggests that those residues were under constraint from undergoing nonsynonymous changes during evolution of the genus. No fixed nonsynonymous differences were found between Pfk13 and its orthologs in closely related species found in African apes. This indicates that all nonsynonymous substitutions currently found in Pfk13 are younger than the time of divergence between P. falciparum and its closely related species. At the population level, no mutations linked to delayed parasite clearance were found in our samples from Africa and Haiti. However, there is a high number of single Pfk13 mutations segregating in P. falciparum populations, and two predominant alleles are distributed worldwide. This pattern is discussed in terms of how changes in the efficacy of natural selection, affected by population expansion, may have allowed for the emergence of mutations tolerant to ACTs.

INTRODUCTION

Antimalarial resistance is one of the greatest threats to malaria control and its eventual elimination (1, 2). The emergence of a delayed parasite clearance phenotype following artemisinin-based therapies in Plasmodium falciparum was first reported in 2007 and 2008 in western Cambodia and was later reported in Myanmar, Thailand, Vietnam, Lao People’s Democratic Republic (Lao PDR), and South China (2–7). The occurrence of this phenotype has raised concerns that it may spread to Africa or lead to resistance, following the previous path of chloroquine (CQ) and sulfadoxine-pyrimethamine (SP) resistance (8–9). However, there is currently no evidence that this phenotype has spread outside of Southeast Asia (SEA), nor is there evidence of the emergence of resistant parasites (2).

Artemisinin delayed parasite clearance was first defined and confirmed using in vitro ring-stage survival assays (RSAs) (1, 10, 11). A series of genome-wide association studies have shown that this phenotype is associated with nonsynonymous mutations in the propeller domain of a P. falciparum kelch gene (codons 441 to 726) found on chromosome 13 (Pfk13) (12, 13). To date, 14 independent Pfk13 mutations have been associated with clinical delayed parasite clearance (2), and of these, only 5 have been validated by in vivo and in vitro experiments (2): N458Y, Y493H, R539T, I543T, and C580Y. Furthermore, there is evidence of the independent emergence of the C580Y lineage across SEA (14, 15). The most frequent and widely detected Pfk13 mutations found in Cambodia, Vietnam, and Lao PDR are Y493H, R539T, I543T, P553L, and C580Y (5, 7). The mutations F446L, P553L, N458Y, R561H, P574L, and C580Y are the most frequently observed in western Thailand, Myanmar, and China (7, 16–18). Only two mutations, P553L and C580Y, are found across SEA, with the latter mutation being more widespread in eastern Thailand (16–18).

In contrast, a high number of single nonsynonymous mutations, including some associated with delayed parasite clearance, have been found at a low frequency in Africa (19). However, many of the mutations reported in Africa have not expanded to the local parasite populations (2, 19–23). In the case of South America, only Pfk13 C580Y mutant parasites have been found in Guyana, with the flanking microsatellite profiles being different from those observed in SEA, suggesting an independent origin (24). Studies in other countries in South America have searched for polymorphisms in the Pfk13 gene, but no evidence of these mutations has been found to date outside Guyana (25, 26).

Although our understanding of which specific mutations within the PfK13 domain are associated with the delayed parasite clearance phenotype is still incomplete (especially in the context of the genetic background of parasites outside of SEA), the identification of polymorphisms in the Pfk13 gene has been used as a valuable genetic marker for the rapid detection, tracking, and monitoring of this phenotype. Here, to better understand the evolution of the k13 gene, its divergence across the genus Plasmodium and its genetic variation in P. falciparum populations were studied. In particular, we explored the mode and rate of evolution of the k13 gene by analyzing its orthologs across the Plasmodium genus and 982 Pfk13 sequences, including those in public databases, along with new data obtained by next-generation sequencing of African and Haitian isolates.

RESULTS

Phylogenetic-based analyses of the k13 gene.This investigation reports 10 new k13 orthologous gene sequences, obtained by nested PCR, from macaque and gibbon parasite strains and chimpanzees infected with Plasmodium reichenowi, Plasmodium billcollinsi, and Plasmodium gaboni. The sequences obtained here were aligned with those available in public databases (see Materials and Methods and Data Set S1 in the supplemental material) (27). The alignment constructed with the Plasmodium species orthologous k13 genes showed that the encoded protein is remarkably conserved across the genus. Estimates of the evolutionary divergence among orthologous k13 gene sequences were low at both the nucleotide and amino acid levels compared to that of other genes (28). The k13 gene was almost identical among species belonging to the subgenus Laverania. For example, the genetic distance between P. falciparum and Plasmodium praefalciparum from gorillas was 0.00, and that between P. falciparum and P. reichenowi from chimpanzees was only 0.01 ± 0.00 (Data Sets S2 and S3). Furthermore, all k13 orthologous genes displayed a high degree of conservation in the kelch propeller domain among Plasmodium spp. (Fig. 1). In particular, the number of base substitutions per site between the two major malarial parasite k13 orthologs (Pfk13 and Plasmodium vivax k13 [Pvk13]) was low (0.09 ± 0.01) (Data Set S2) and showed synonymous substitutions between P. falciparum and P. vivax, with the exception of V568I versus V568G in Pfk13 and Pvk13, respectively (Fig. 1; Data Sets S1 and S4).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Alignment of the Plasmodium K13 kelch propeller domain. The Pfk13 mutations that have been associated with clinical delayed parasite clearance (2) are shown. The five mutations validated in vivo and in vitro are indicated with an asterisk. The Haemoproteus tartakovskyi (an avian haemosporidian parasite) k13 gene was included as an outgroup.

To compare all Plasmodium species k13 orthologous genes considered in this study (including different strains of some parasites), a Bayesian phylogeny was estimated using the nucleotide sequence alignment, given the high conservation at the protein level (Fig. 2). As expected, the phylogeny was consistent with earlier reports using putatively neutral genes, genes encoding merozoite antigens, and/or mitochondrial genomes (27–29). Overall, the Laverania subgenus, which includes P. falciparum, formed a monophyletic group, as did P. vivax and its closely related nonhuman primate malaria parasites from SEA and all rodent malaria parasites. All mammal Plasmodium species shared a common ancestor with avian Plasmodium species (Fig. 2).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Bayesian phylogenetic hypothesis of Plasmodium parasites based on the k13 gene (42 sequences; 1,985 bp, excluding gaps). The values at the nodes are posterior probabilities (see Materials and Methods). Haemoproteus tartakovskyi (an avian haemosporidian parasite) k13 gene was included as an outgroup.

Then, using the Bayesian tree topology, the relative evolutionary rates of the k13 gene were estimated in a maximum likelihood framework with the RelTime method (30). This revealed a pattern in which there was a faster k13 evolutionary rate in malaria parasites infecting Southeast Asian macaques and their relatives, including the human malaria parasite P. vivax (Fig. 3). The k13 relative evolutionary rates were approximately 12 times higher in malaria parasites infecting Southeast Asian macaques and their relatives than in parasites from other mammalian hosts and African ape parasites. Differences in rates have also been observed in other genes in this clade, which includes the human parasite P. vivax (31, 32). Interestingly, Pfk13 also showed a k13 relative rate of evolution faster than and similar to that in SEA parasites (Fig. 3). The acceleration in the P. falciparum k13 evolutionary rate, however, did not translate into the gene, being under diversifying natural selection when P. falciparum diverged from its closest related species (see below).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Relative evolutionary rates of the Plasmodium k13 gene. Branches are colored from high (green) to low (purple) according to their relative rates with respect to the rate at the root, which was set to 1, as estimated by RelTime without calibration constraints.

The relative evolutionary rate differences in the k13 gene between clades were further examined by estimating synonymous and nonsynonymous substitution rates per codon per million years (My) using a codon model implemented in the CodonRates (v1.0) program (33). In this analysis, the k13 nonsynonymous substitution rate (0.022 ± 0.012 per codon per My) was slightly higher (on average) than the synonymous substitution rate (0.012 ± 0.006 per codon per My) across the genus. The concordance statistic (S) showed that synonymous and nonsynonymous substitution rate changes over time were correlated (P < 0.5), suggesting that there was no pattern indicative of positive selection acting on this gene. The branch leading to P. falciparum (27, 28) from the common ancestor with its closest sister taxon (Plasmodium praefalciparum) showed a synonymous substitution rate of 0.005 per codon per My, while the nonsynonymous substitution rate was 0.001 per codon per My.

To further assess the role that selection may have had during the divergence of the k13 gene, we utilized methods implemented in the HyPhy package (34) (see Materials and Methods). A test that examines each individual site independently for evidence of episodic diversification (MEME [35]) found no evidence of episodic positive/diversifying selection acting on any site of the gene. In addition, a branch site-unrestricted statistical test for episodic diversification (BUSTED) (36), found no evidence (likelihood ratio test, P = 0.997, which is ≥0.05) of gene-wide episodic diversifying selection in any of the branches of the k13 phylogeny. Finally, RELAX (37), a test for identifying whether the strength of natural selection (positive or negative) has been relaxed or intensified along a specified set of test branches (e.g., P. falciparum), was not significant (K = 0.74, P = 0.665, likelihood ratio = 0.19). These results showed no evidence of positive diversifying selection acting on the k13 gene either during the evolution of the sampled Plasmodium species or in the lineages of P. falciparum or P. vivax.

Plasmodium falciparum population analyses.Using an alignment with 982 complete k13 gene sequences from around the globe (519 from public databases and 463 reported in this study; see Materials and Methods and Data Set S5), low genetic diversity was found in the Pfk13 gene, as measured by π (0.0005 ± 0.0; Table 1). In contrast, the polymorphism measured using the number of polymorphic (segregating) sites had an S value of 146. The number of singleton variable sites, considering two variants, was 108. It is worth noting that no mutations associated with the delayed parasite clearance phenotype were found in our samples from Africa and Haiti. In the case of P. vivax, all 56 Pvk13 sequences available from the databases were identical.

View this table:
  • View inline
  • View popup
TABLE 1

Polymorphism found in k13 gene of Plasmodium falciparum

Evidence of natural selection acting on the Pfk13 polymorphism was explored by estimating the average pairwise number of synonymous substitutions per synonymous site (dS) and the average pairwise number of nonsynonymous substitutions per nonsynonymous site (dN). This revealed a significant (P < 0.05) excess of nonsynonymous over synonymous substitutions only in Asia, the region where mutations associated with the artemisinin delayed parasite clearance phenotype were first reported (Table 1).

A total of 147 Pfk13 haplotypes were found with a haplotype (allele) diversity (Hd) of 0.742 (standard deviation [SD] = 0.012). The uneven sampling effort did not allow us to estimate the haplotype relative frequencies per country. Overall, the median-joining network for the Pfk13 gene, estimated using Network software (Fluxus Technologies, 2011), showed two star-like shapes with at least two most predominant and common haplotypes (haplotype 1 [H1] and H2; Fig. 4A) shared by the three continents (Asia, Africa, and South America) with a broad distribution. Interestingly, one of the haplotypes (H1, n = 442, 45%) is mostly distributed in Asia, and the other (H2, n = 223, 22.7%) is mostly distributed in Africa. A third haplotype (H3, n = 50, 5%), in terms of frequency, is only found in Thailand (Fig. 4A). A noteworthy finding was that haplotype H3 contains the C580Y mutation found in the SEA region. In addition, a less frequent haplotype was found only on the China-Myanmar border (H4, n = 28, 2.9%) and had the F446I mutation, one of the most common mutations found in China and on the China-Myanmar border (Fig. 4B). However, most of the haplotypes were restricted to single countries (Fig. 4A).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Median-joining network of the P. falciparum k13 haplotypes. (A) Median-joining network of the P. falciparum k13 haplotypes using 982 cosmopolitan complete sequences. All the k13 gene sequences obtained here for P. falciparum population samples (n = 463), as well as those available in the PlasmoDB and NCBI databases (n = 519), were included. Branch lengths are proportional to divergence; node sizes are proportional to the total haplotype frequencies. The network shows 147 haplotypes found in 982 sequences with a haplotype diversity (Hd) of 0.742 (SD = 0.012). Every color corresponds to a different geographic origin. Lines separating haplotypes represent mutational steps. The four most frequent haplotypes are indicated. (B) Frequencies of candidate/validated mutations in the propeller domain of the Pfk13 gene that have been associated with clinical delayed parasite clearance by country (Table 2). Frequencies were estimated using an alignment constructed with 3,916 partial (n = 2,934) and complete (n = 982) k13 gene sequences. C-M, China-Myanmar.

Finally, using 3,912 partial and complete Pfk13 sequences (Data Set S5), we estimated the frequencies (by country) of candidate/validated mutations associated with a clinical delay in parasite clearance (2). Those mutations, found in the propeller domain of a Pfk13 gene, are shown in Fig. 4B (see also Table 2). As previously reported, mutations F446I and C580Y were the most frequent mutations in this data set and were mostly found in Asia (Fig. 4B; Table 2; Data Set S6).

View this table:
  • View inline
  • View popup
TABLE 2

Number of Pfk13 sequences with mutations (candidates/validated) associated with clinical delayed parasite clearance (per codon and per country)a

DISCUSSION

Unlike other genes with mutations associated with antimalarial drug resistance, such as P. falciparum dhfr (Pfdhfr) (38, 39), the Pfk13 gene shows several single mutations at a low frequency, and many of these have no clear associations with the delayed parasite clearance phenotype (2). However, similar to genes such as Pfdhfr (28, 38), the Pfk13 mutations associated with this phenotype appeared solely in the P. falciparum lineage, indicating that they are all of recent origin (40). Although they evolve at heterogeneous rates, as indicated by the relative rate analysis, it is evident that orthologous k13 genes are highly conserved across Plasmodium species (in several species of human, nonhuman primate, other mammalian, and avian parasites). Evidence specific to the P. falciparum lineage can be found in the low nonsynonymous substitution rate estimated in this study. This finding suggests that the k13 gene was under a strong constraint (negative selection) from undergoing natural variation during its evolution in the genus Plasmodium.

There was no evidence of fixed nonsynonymous differences between P. falciparum and its closest related Laverania species (P. praefalciparum and P. reichenowi). Together, this evidence indicates that all Pfk13 mutations, including those associated with delayed parasite clearance, are of recent occurrence. The presence of several Pfk13 mutations found at a low frequency translates into the extraordinary haplotype diversity found in this study. A naive interpretation of this pattern could find it to be contradictory to the strong conservation observed among Plasmodium species at the protein level (indicative of negative selection). However, this finding is consistent with the proposed historical population expansion of P. falciparum following expansion of the human population (41). Population expansions can increase the number of segregating mutations (even semideleterious mutations) by making negative selection less effective (42, 43). Nonetheless, as parasite populations expand, negative selection becomes more efficient in eliminating those semideleterious mutations and maintaining them at a low frequency. The spatial differential distribution of those Pfk13 mutations at a low frequency is also expected because P. falciparum is not a single population, so there are local population expansions where different mutations occur. This pattern of population expansions is summarized in Tajima’s D results, where the low value of π and a high number of segregating sites (mostly singletons) lead to a significant negative estimate for the statistic D (−2.683, P < 0.001). Thus, considering the parasite demography, the Pfk13 gene is expected to have several recent geographically distinct mutations at a low frequency, with each individual parasite lineage carrying one or very few of those segregating sites (42, 43). Although in P. vivax there is also evidence of a recent population expansion (44, 45), no segregating mutations were found in Pvk13. Given the limitation of our sample (n = 53), we cannot discern whether this is a real pattern or a sampling artifact.

The data are consistent with both selection and demography driving the dynamic of Pfk13 mutations conferring artemisinin-based combination therapy (ACT) tolerance (46). The parasite population expansion (demography) relaxes purifying (negative) selection, allowing for semideleterious mutations to segregate (41, 43, 46, 47), and the introduction of ACTs selects for mutations conferring delayed parasite clearance/tolerance out of those several variants present at a low frequency, even if such variants have a fitness cost in the absence of drug pressure (8, 46–48). Under such dynamics, as the data in Africa seem to indicate, an asymptomatic untreated reservoir will favor ACT-sensitive parasites, keeping a negative selective pressure on those semideleterious mutations conferring tolerance to the drugs (47–50). Consequently, the Pfk13 mutations present today are unlikely to increase in frequency in regions such as Africa, where their relative advantage is still overcome by their fitness cost in the absence of treatment (e.g., subclinical or asymptomatic infections at a high frequency [49, 50]), and parasites with those mutatiosn are still susceptible to aggressive/monitored drug treatment (8, 46). The interplay between the putative fitness cost of these mutations and the absence of drug pressure in asymptomatic individuals could be facilitated by the higher population sizes of P. falciparum in Africa, making differences in selection regimens prevail over genetic drift in the dynamics of these mutations (46, 48). Under the current scenario, interventions directed to reduce the asymptomatic reservoir, such as targeted mass drug administration (TMDA) with ACTs, will remain effective because there are not yet Pfk13 mutations conferring complete resistance (46, 51).

Although the observed patterns are robust, it is worth noting that many studies in P. falciparum directly report the frequency of mutations related to ACT-delayed parasite clearance. This practice not only may lead to underreporting of mutations occurring at a low frequency but also limits our capacity to track geographic changes in haplotype frequencies or infer other patterns. As an example, the frequency of the C580Y mutation is high in western Cambodia (52, 53), but the deposited sequence data are scarce. Whether there is an expansion of a particular haplotype in the parasite population cannot be addressed across studies by using mutation frequency data.

The proposed dynamics explain why the actions deployed thus far have limited the expansion of these mutations from Asia, and there has been only a single haplotype of C580Y identified in Guyana. Whether the spread will remain limited (17, 24) is a matter that requires continued monitoring. Importantly, reporting novel mutations is critical since they may change the parasite fitness under drug pressure. Thus, considering that all observed mutations in Pfk13 are of recent origin, understanding their evolutionary dynamics at the early stages of drug resistance could further inform containment strategies wherever putative resistance mutations start to increase in frequency.

MATERIALS AND METHODS

Samples/sequences from the genus Plasmodium.Orthologous k13 genes were sequenced from macaque and gibbon parasite laboratory strains provided by the Centers for Disease Control and Prevention (CDC): Plasmodium inui (OS strain), Plasmodium hylobati, Plasmodium fieldi (N-3 strain), Plasmodium simiovale, and Plasmodium gonderi. Information on these species and strains can be found elsewhere (54). Archived remainders of blood samples from chimpanzees infected with Plasmodium reichenowi (Shegue and 20 isolates [28]), Plasmodium billcollinsi (the Ithaito isolate [28]), and Plasmodium gaboni (Shegue and 39 isolates [28]) were also used in this study (see Data Set S1 in the supplementary material). The chimpanzees were housed at the Jane Goodall Institute (JGI) Tchimpounga Chimpanzee Rehabilitation Center in the Democratic Republic of Congo (DRC), and their blood samples were collected throughout 2009 and 2010 by the veterinary staff as part of the chimpanzees’ routine health examinations, following standards approved by the Pan Africa Sanctuary Alliance. More information on these species and ethical guidelines can be found in reference 28.

In addition, all orthologous k13 gene sequences from different nonhuman Plasmodium species and Plasmodium vivax with publicly available genomes were retrieved from the PlasmoDB (v38) (55) and MalAvi (56) databases (Data Sets S1 and S4).

Plasmodium falciparum samples/sequences.Plasmodium falciparum samples from Africa (n = 244) and Haiti (n = 103) and samples from U.S. malaria cases without travel histories (n = 116) reported to the CDC malaria surveillance network from 2014 to 2017 (n = 463) were used (Data Set S5). Genomic DNA isolated from P. falciparum strains 7G8 and HB3 was used as a positive control, and DNA isolated from whole blood from individuals without malaria was used as a negative control. The use of U.S. surveillance P. falciparum samples was approved by the Office of the Associate Director of Science, Center for Global Health, Centers for Disease Control and Prevention, as research not involving human subjects (protocol 2017-192). In addition, the methodology employed here for the amplification of the k13 gene was also approved by the same office as a surveillance activity (protocol 2017-309).

The data were expanded by retrieving complete and partial k13 gene sequences for P. falciparum from publicly available genomes/sequences found in the PlasmoDB (55) and NCBI (57) databases. All sequences were analyzed together with the new data obtained here from African and Haitian isolates (Data Set S5).

DNA extraction and Sanger sequencing for nonhuman primate parasites.Genomic DNA was extracted from whole blood of the nonhuman primate Plasmodium species using a DNeasy blood and tissue kit (Qiagen, GmbH, Hilden, Germany). Orthologous k13 genes from macaque and gibbon parasites were amplified by nested PCR using external primers forward 5′-C ATT TCC AAC TTC TCC GTC-3′ (58) and reverse 5′-TAT ATT TGC TAT TAR NAC NGA RTG-3′ and internal primers forward 5′-A TGA CGT ATG ADA GRG ART CNG-3′ and reverse 5′-AA TCT GGG AAC TAA TAR DGR DGG-3′. The primary PCR amplifications were carried out in a 50-μl reaction volume using 5 μl of total genomic DNA, 2.5 mM MgCl2, 1× PCR buffer, 1.25 mM each deoxynucleoside triphosphate, 0.4 mM each primer, and 0.03 U/μl AmpliTaq polymerase (Applied Biosystems, Roche, USA). The primary PCR conditions were a partial denaturation at 94°C for 4 min and 36 cycles with 1 min at 94°C, 1 min at 54°C, and a 2-min extension at 72°C, and a final extension of 10 min at 72°C was added in the last cycle. Then, the nested PCR mixtures were also made in 50-μl reaction volume using 1 μl of the primary PCR mixture, 2.5 mM MgCl2, 1× PCR buffer, 1.25 mM each deoxynucleoside triphosphate, 0.4 mM each primer, and 0.03 U/μl AmpliTaq polymerase. The nested PCR conditions were a partial denaturation at 94°C for 4 min and 25 cycles with 1 min at 94°C, 1 min at 56°C, and a 2-min extension at 72°C, and a final extension of 10 min at 72°C was added in the last cycle. Then, orthologous k13 genes from chimpanzee parasites, which share a recent common ancestor with P. falciparum, were amplified by using the same protocol described above, but with external primers forward 5′-GAA TTT TTC TAT NAC RTA YGA DAG-3′ and reverse 5′-ATT TGC TAT TAR NAC NGA RTG NCC-3′ and internal primers forward 5′-A TGA CGT ATG ADA GRG ART CNG-3′ and reverse 5′-AA TCT GGG AAC TAA TAR DGR DGG-3′. Then, two independent PCR products (50 μl) were excised from the gel (bands of approximately 2 kbp), purified using a QIAquick gel extraction kit (Qiagen, GmbH, Hilden, Germany), and cloned into the pGEM-T Easy vector system (Promega, USA). To detect mixed infections, at least two independent PCR products were cloned, and four clones were sequenced from each sample. Both strands for all the k13 orthologous genes were sequenced using an Applied Biosystems 3730 capillary sequencer. All the sequences obtained were identified as the k13 gene using analysis with the BLAST program (59). The sequences were submitted to GenBank and may be found under accession numbers MK103884 to MK103893.

DNA extraction from human samples and next-generation sequencing.DNA from human samples was extracted using the whole-blood protocol for the QIAamp DNA blood kit (Qiagen, GmbH, Hilden, Germany). Then, an end-to-end Illumina targeted amplicon deep sequencing (TADS) protocol and bioinformatics pipeline for molecular surveillance of drug resistance in P. falciparum, called malaria resistance surveillance (MaRS), was used to amplify the Pfk13 gene (60). First, the gene was amplified using a New England BioLabs (NEB) High Fidelity PCR kit (catalog no. 51104; New England BioLabs, USA) according to the manufacturer’s instructions with a 50.0-μl master mix preparation using 5× GC buffer. Next, 10.0 μl of each PCR product for each sample was used to generate the next-generation sequencing library. Pooled library fragment sizes were checked using an Agilent D5000 ScreenTape station (catalog numbers G2940C, 5067-4626, and G2940CA; Agilent Technologies), and the pooled library concentration was measured using a Qubit (v3.0) fluorometer (catalog numbers Q33216 and Q32853; Life Technologies Corp.). Lastly, a MiSeq desktop sequencer (Illumina Inc.) was used to obtain the sequences, and a consensus sequence for each sample was assembled using a pipeline. Additional details about the MaRS protocol are provided in Text S1 in the supplemental material of reference 60, including primer sequences, thermocycling conditions, and library amplification, purification, quantification, and normalization. The Pfk13 sequences obtained here were submitted to the GenBank database and may be found under accession numbers MK103421 to MK103883.

Phylogenetic-based method analyses of the k13 gene.First, a nucleotide sequence alignment was performed by using 42 sequences (1,985 bp, excluding gaps), including the entire orthologous k13 gene sequences obtained here, the sequences of the different haplotypes for some well-identified parasite species (e.g., P. reichenowi and P. gaboni), as well as sequences available in the PlasmoDB and MalAvi databases (Fig. 1 and 2; Data Set S1). Then, a second alignment based on 28 orthologous k13 gene sequences selected from different species of the first alignment was performed (duplicate haplotypes were not included; Fig. 3). All the alignments were generated using the ClustalX (v2.0.12) and Muscle programs implemented in the SeaView (v4.3.5) program (61) with manual editing. In both alignments, the Haemoproteus tartakovskyi (an avian haemosporidian parasite) k13 gene was included as an outgroup (56).

Using both alignments, phylogenetic relationships were inferred using the Bayesian method implemented in MrBayes (v3.2.6) software with the default priors (62). A general time-reversible model with gamma-distributed substitution rates and a proportion of invariant sites (GTR + Γ + I), which was the model with the lowest Bayesian information criterion (BIC) scores, as estimated by MEGA (v7.0.26) software, was used (63). Bayesian support for the nodes was inferred in MrBayes by sampling every 1,000 generations from two independent chains lasting 4 × 106 Markov chain Monte Carlo (MCMC) steps. The chains were assumed to have converged once the average standard deviation of the posterior probability was below 0.01 and the value of the potential scale reduction factor (PSRF) was between 1.00 and 1.02 (62). As a burn-in, 25% of the sample was then discarded once convergence was reached. The phylogenetic trees were visualized and edited using the FigTree (v1.4.3) program (64).

To detect changes in the rate of evolution, relative evolutionary rates of the k13 gene were estimated on the second alignment using a maximum likelihood framework with the RelTime method (30), as implemented in the command line version implemented in MEGA (vX) software (65). The substitution model was the same as that used for Bayesian analyses (GTR + Γ + I). Then, the tree was visualized in the FigTree (v1.4.3) program and colored using the relative evolutionary rate estimates of the k13 gene (Fig. 3). To quantitate the rate of evolution of this gene, absolute rates of synonymous and nonsynonymous substitutions per codon per My for the k13 gene were obtained using the CodonRates (v1.0) program in a Bayesian framework (33). In addition, by disentangling synonymous from nonsynonymous substitution rates, this analysis provides a first assessment of how natural selection may have affected the evolution of the gene. Because this method requires a tree topology and time calibration priors, including one for the in-group root node (the Plasmodium clade in this case), the Bayesian methods program BEAST (v2.4.1) (66) was used to estimate the origin of the Plasmodium species included in this investigation.

In BEAST, the relaxed clock (67) was applied with the GTR + Γ + I evolutionary model with a Yule model as the tree prior. Two independent chains lasting 108 MCMC steps were run until convergence. One scenario, using a combination of two calibration constraints that consider information on vertebrate hosts’ fossils (see the Paleobiology Database at http://www.paleodb.org/; last accessed 13 November 2017) was explored, and the results were compared with those of earlier studies based on the mitochondrial genomes (29). In all cases, uniform priors were used to calibrate divergences. Therefore, no particular time point within a given time interval was favored, and no punctual calibrations were used. The first calibration constraint considers that African parasites found in Mandrillus spp. (Plasmodium gonderi) diverged from those Plasmodium species found in Southeast Asia macaques (a monophyletic group) when Macaca branched from Papio (29, 68, 69). This event assumed that the divergence at that node took place at any time between 6 and 14.2 My ago (Ma) using a uniform prior. This time interval includes molecular estimates for the same Papio-Macaca divergence event (29, 69), and it is also consistent with the age of older fossils reported for Macaca spp. (see the Paleobiology Database at http://www.paleodb.org/; last accessed 13 November 2017). This calibration constraint considers the fact that P. gonderi has also been seen in Chlorocebus (a member of the Cercopithecini). The second fossil-based calibration constraint is the minimum of 23.5 Ma for the human-Macaca split (70), which is assumed to be the minimum time when the monophyletic group that includes P. malariae (a parasite found in humans) and all the Asian parasites originated. This event is used to inform a calibration prior with a maximum of 65 Ma that allows this parasite clade to be as old as the origin of primates (71). More details about this method can be found elsewhere (29, 69).

Given that the time estimated here for the avian Plasmodium-primate/rodent Plasmodium split (45.36 Ma; 95% credibility intervals, 33.23 to 61.10 Ma) was comparable to previous estimates also obtained using BEAST (45.31 Ma; credibility intervals, 38.95 to 51.96 [29]), the absolute rates of synonymous and nonsynonymous substitutions per codon per My were then estimated by using CodonRates (33) with a calibration prior of 45.4 ± 10 Ma for the in-group root (genus Plasmodium), as estimated from BEAST. This analysis was performed on a tree topology estimated with MrBayes that used Haemoproteus tartakovskyi as an outgroup.

Finally, to elucidate selective pressures during the evolution of the k13 gene, several methods were performed with the HyPhy package (34) implemented in the DataMonkey web server (http://www.datamonkey.org/; last accessed August 2018). Evidence of both episodic and pervasive positive selection at the level of an individual site was identified using a mixed-effects model of evolution, MEME (35), which detects sites evolving under positive selection under a proportion of branches. To investigate whether the k13 gene has experienced positive selection at least at one site on at least one branch, a branch site-unrestricted statistical test for episodic diversification (BUSTED) was used. This provides a gene-wide (not site-specific) test for positive selection (36). Sites positively selected were identified at a significance level (P value) of <0.05. In addition, to test whether the strength of natural selection has been relaxed or intensified along a specified set of test branches, the RELAX method (37) was performed. RELAX is most useful for identifying trends and/or shifts in the stringency of natural selection on a given gene.

Plasmodium falciparum population analyses.To study the genetic relationships among worldwide haplotypes of the k13 gene, a third alignment was performed using 982 cosmopolitan complete sequences, including all the k13 gene sequences obtained here for P. falciparum population samples (n = 463), as well as those available in the PlasmoDB and NCBI databases (n = 519). Evidence of natural selection was explored by estimating the average number of synonymous substitutions per synonymous site (dS) and the average number of nonsynonymous substitutions per nonsynonymous site (dN) between a pair of sequences under the Nei-Gojobori method with the Jukes and Cantor corrections, as implemented in the MEGA (v7) program (63). The difference between dS and dN with its standard error (estimated using bootstrap analysis with 1,000 pseudoreplications) was calculated, and significance was tested by a two-tailed codon-based Z-test on the difference between dS and dN, as described by Nei and Kumar (72). The null hypothesis was that the observed polymorphism is not under selection (dS = dN).

The alignments were used to estimate median-joining networks using Network (v4.6.1.0) software (Fluxus Technologies, 2011). Transversions were set equal to transitions, and the epsilon parameter was set equal to 0 with only one round of star contraction, which collapses star-like structures in the network into single nodes. The total number of sites included in this analysis, excluding gaps or missing data, was 2,137 bp out of 2,181 bp using the P. falciparum 3D7 strain as a reference. Estimates of nucleotide diversity, haplotype diversity, and Tajima’s D test parameter were obtained using DnaSP (v5.50) software (73). Lastly, the frequency of candidate/validated mutations located in the propeller domain of a P. falciparum kelch (Pfk13) gene that have been associated with clinical delayed parasite clearance was estimated by country, as well as by continent, using a fourth alignment constructed with 3,916 partial (n = 2,934) and complete (n = 982) k13 gene sequences (Table 2; Data Sets S5 and S6).

Accession number(s).The Pfk13 sequences obtained here were deposited in GenBank under accession numbers MK103421 to MK103893.

ACKNOWLEDGMENTS

We thank Michael Li for the silhouettes design and the DNA laboratory at the School of Life Sciences (Arizona State University) for technical support.

This study was supported in part by grants from the U.S. National Institutes of Health (2U19 AI089681-08 and R01AI130473). Computer resources were supported by a grant from the Pennsylvania Department of Health (TU-420721).

We declare that no conflicts of interest exist.

FOOTNOTES

    • Received 5 December 2018.
    • Returned for modification 24 February 2019.
    • Accepted 28 April 2019.
    • Accepted manuscript posted online 13 May 2019.
  • Supplemental material for this article may be found at https://doi.org/10.1128/AAC.02550-18.

  • Copyright © 2019 American Society for Microbiology.

All Rights Reserved.

REFERENCES

  1. 1.↵
    1. Dondorp AM,
    2. Smithuis FM,
    3. Woodrow C,
    4. Seidlein L
    . 2017. How to contain artemisinin- and multidrug-resistant falciparum malaria. Trends Parasitol 33:353–363. doi:10.1016/j.pt.2017.01.004.
    OpenUrlCrossRef
  2. 2.↵
    WHO. 2017. Artemisinin and artemisinin-based combination therapy resistance: status report. World Health Organization, Geneva, Switzerland.
  3. 3.↵
    1. Noedl H,
    2. Se Y,
    3. Schaecher K,
    4. Smith BL,
    5. Socheat D,
    6. Fukuda MM
    , Artemisinin Resistance in Cambodia 1 (ARC1) Study Consortium. 2008. Evidence of artemisinin-resistant malaria in western Cambodia. N Engl J Med 359:2619–2620. doi:10.1056/NEJMc0805011.
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    1. Dondorp AM,
    2. Nosten F,
    3. Yi P,
    4. Das D,
    5. Phyo AP,
    6. Tarning J,
    7. Lwin KM,
    8. Ariey F,
    9. Hanpithakpong W,
    10. Lee SJ,
    11. Ringwald P,
    12. Silamut K,
    13. Imwong M,
    14. Chotivanich K,
    15. Lim P,
    16. Herdman T,
    17. An SS,
    18. Yeung S,
    19. Singhasivanon P,
    20. Day NP,
    21. Lindegardh N,
    22. Socheat D,
    23. White NJ
    . 2009. Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med 361:455–467. doi:10.1056/NEJMoa0808859.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. Ashley EA,
    2. Dhorda M,
    3. Fairhurst RM,
    4. Amaratunga C,
    5. Lim P,
    6. Suon S,
    7. Sreng S,
    8. Anderson JM,
    9. Mao S,
    10. Sam B,
    11. Sopha C,
    12. Chuor CM,
    13. Nguon C,
    14. Sovannaroth S,
    15. Pukrittayakamee S,
    16. Jittamala P,
    17. Chotivanich K,
    18. Chutasmit K,
    19. Suchatsoonthorn C,
    20. Runcharoen R,
    21. Hien TT,
    22. Thuy-Nhien NT,
    23. Thanh NV,
    24. Phu NH,
    25. Htut Y,
    26. Han KT,
    27. Aye KH,
    28. Mokuolu OA,
    29. Olaosebikan RR,
    30. Folaranmi OO,
    31. Mayxay M,
    32. Khanthavong M,
    33. Hongvanthong B,
    34. Newton PN,
    35. Onyamboko MA,
    36. Fanello CI,
    37. Tshefu AK,
    38. Mishra N,
    39. Valecha N,
    40. Phyo AP,
    41. Nosten F,
    42. Yi P,
    43. Tripura R,
    44. Borrmann S,
    45. Bashraheil M,
    46. Peshu J,
    47. Faiz MA,
    48. Ghose A,
    49. Hossain MA,
    50. Samad R
    , et al. 2014. Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med 371:411–423. doi:10.1056/NEJMoa1314981.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Saunders DL,
    2. Vanachayangkul P,
    3. Lon C
    , U.S. Army Military Malaria Research Program, National Center for Parasitology Entomology and Malaria Control. Royal Cambodian Armed Forces. 2014. Dihydroartemisinin-piperaquine failure in Cambodia. N Engl J Med 371:484–485. doi:10.1056/NEJMc1403007.
    OpenUrlCrossRef
  7. 7.↵
    1. Imwong M,
    2. Suwannasin K,
    3. Kunasol C,
    4. Sutawong K,
    5. Mayxay M,
    6. Rekol H,
    7. Smithuis FM,
    8. Hlaing TM,
    9. Tun KM,
    10. van der Pluijm RW,
    11. Tripura R,
    12. Miotto O,
    13. Menard D,
    14. Dhorda M,
    15. Day NPJ,
    16. White NJ,
    17. Dondorp AM
    . 2017. The spread of artemisinin-resistant Plasmodium falciparum in the Greater Mekong subregion: a molecular epidemiology observational study. Lancet Infect Dis 17:491–497. doi:10.1016/S1473-3099(17)30048-8.
    OpenUrlCrossRef
  8. 8.↵
    1. Escalante AA,
    2. Smith DL,
    3. Kim Y
    . 2009. The dynamics of mutations associated with anti-malarial drug resistance in Plasmodium falciparum. Trends Parasitol 25:557–563. doi:10.1016/j.pt.2009.09.008.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Prosser C,
    2. Meyer W,
    3. Ellis J,
    4. Lee R
    . 2018. Evolutionary ARMS race: antimalarial resistance molecular surveillance. Trends Parasitol 34:322–334. doi:10.1016/j.pt.2018.01.001.
    OpenUrlCrossRef
  10. 10.↵
    1. Witkowski B,
    2. Amaratunga C,
    3. Khim N,
    4. Sreng S,
    5. Chim P,
    6. Kim S,
    7. Lim P,
    8. Mao S,
    9. Sopha C,
    10. Sam B,
    11. Anderson JM,
    12. Duong S,
    13. Chuor CM,
    14. Taylor WR,
    15. Suon S,
    16. Mercereau-Puijalon O,
    17. Fairhurst RM,
    18. Menard D
    . 2013. Novel phenotypic assays for the detection of artemisinin-resistant Plasmodium falciparum malaria in Cambodia: in-vitro and ex-vivo drug-response studies. Lancet Infect Dis 13:1043–1049. doi:10.1016/S1473-3099(13)70252-4.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Witkowski B,
    2. Khim N,
    3. Chim P,
    4. Kim S,
    5. Ke S,
    6. Kloeung N,
    7. Chy S,
    8. Duong S,
    9. Leang R,
    10. Ringwald P,
    11. Dondorp AM,
    12. Tripura R,
    13. Benoit-Vical F,
    14. Berry A,
    15. Gorgette O,
    16. Ariey F,
    17. Barale JC,
    18. Mercereau-Puijalon O,
    19. Menard D
    . 2013. Reduced artemisinin susceptibility of Plasmodium falciparum ring stages in western Cambodia. Antimicrob Agents Chemother 57:914–923. doi:10.1128/AAC.01868-12.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Ariey F,
    2. Witkowski B,
    3. Amaratunga C,
    4. Beghain J,
    5. Langlois A-C,
    6. Khim N,
    7. Kim S,
    8. Duru V,
    9. Bouchier C,
    10. Ma L,
    11. Lim P,
    12. Leang R,
    13. Duong S,
    14. Sreng S,
    15. Suon S,
    16. Chuor CM,
    17. Bout DM,
    18. Ménard S,
    19. Rogers WO,
    20. Genton B,
    21. Fandeur T,
    22. Miotto O,
    23. Ringwald P,
    24. Le Bras J,
    25. Berry A,
    26. Barale J-C,
    27. Fairhurst RM,
    28. Benoit-Vical F,
    29. Mercereau-Puijalon O,
    30. Ménard D
    . 2014. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature 505:50–55. doi:10.1038/nature12876.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    1. Straimer J,
    2. Gnädig NF,
    3. Witkowski B,
    4. Amaratunga C,
    5. Duru V,
    6. Ramadani AP,
    7. Dacheux M,
    8. Khim N,
    9. Zhang L,
    10. Lam S,
    11. Gregory PD,
    12. Urnov FD,
    13. Mercereau-Puijalon O,
    14. Benoit-Vical F,
    15. Fairhurst RM,
    16. Ménard D,
    17. Fidock DA
    . 2015. Drug resistance. K13-propeller mutations confer artemisinin resistance in Plasmodium falciparum clinical isolates. Science 347:428–431. doi:10.1126/science.1260867.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Miotto O,
    2. Amato R,
    3. Ashley EA,
    4. MacInnis B,
    5. Almagro-Garcia J,
    6. Amaratunga C,
    7. Lim P,
    8. Mead D,
    9. Oyola SO,
    10. Dhorda M,
    11. Imwong M,
    12. Woodrow C,
    13. Manske M,
    14. Stalker J,
    15. Drury E,
    16. Campino S,
    17. Amenga-Etego L,
    18. Thanh TN,
    19. Tran HT,
    20. Ringwald P,
    21. Bethell D,
    22. Nosten F,
    23. Phyo AP,
    24. Pukrittayakamee S,
    25. Chotivanich K,
    26. Chuor CM,
    27. Nguon C,
    28. Suon S,
    29. Sreng S,
    30. Newton PN,
    31. Mayxay M,
    32. Khanthavong M,
    33. Hongvanthong B,
    34. Htut Y,
    35. Han KT,
    36. Kyaw MP,
    37. Faiz MA,
    38. Fanello CI,
    39. Onyamboko M,
    40. Mokuolu OA,
    41. Jacob CG,
    42. Takala-Harrison S,
    43. Plowe CV,
    44. Day NP,
    45. Dondorp AM,
    46. Spencer CC,
    47. McVean G,
    48. Fairhurst RM,
    49. White NJ,
    50. Kwiatkowski DP
    . 2015. Genetic architecture of artemisinin-resistant Plasmodium falciparum. Nat Genet 47:226–234. doi:10.1038/ng.3189.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Takala-Harrison S,
    2. Jacob CG,
    3. Arze C,
    4. Cummings MP,
    5. Silva JC,
    6. Dondorp AM,
    7. Fukuda MM,
    8. Hien TT,
    9. Mayxay M,
    10. Noedl H,
    11. Nosten F,
    12. Kyaw MP,
    13. Nhien NT,
    14. Imwong M,
    15. Bethell D,
    16. Se Y,
    17. Lon C,
    18. Tyner SD,
    19. Saunders DL,
    20. Ariey F,
    21. Mercereau-Puijalon O,
    22. Menard D,
    23. Newton PN,
    24. Khanthavong M,
    25. Hongvanthong B,
    26. Starzengruber P,
    27. Fuehrer HP,
    28. Swoboda P,
    29. Khan WA,
    30. Phyo AP,
    31. Nyunt MM,
    32. Nyunt MH,
    33. Brown TS,
    34. Adams M,
    35. Pepin CS,
    36. Bailey J,
    37. Tan JC,
    38. Ferdig MT,
    39. Clark TG,
    40. Miotto O,
    41. MacInnis B,
    42. Kwiatkowski DP,
    43. White NJ,
    44. Ringwald P,
    45. Plowe CV
    . 2015. Independent emergence of artemisinin resistance mutations among Plasmodium falciparum in Southeast Asia. J Infect Dis 211:670–679. doi:10.1093/infdis/jiu491.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Talundzic E,
    2. Okoth SA,
    3. Congpuong K,
    4. Plucinski MM,
    5. Morton L,
    6. Goldman IF,
    7. Kachur PS,
    8. Wongsrichanalai C,
    9. Satimai W,
    10. Barnwell JW,
    11. Udhayakumar V
    . 2015. Selection and spread of artemisinin-resistant alleles in Thailand prior to the global artemisinin resistance containment campaign. PLoS Pathog 11:e1004789. doi:10.1371/journal.ppat.1004789.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Ménard D,
    2. Khim N,
    3. Beghain J,
    4. Adegnika AA,
    5. Shafiul-Alam M,
    6. Amodu O,
    7. Rahim-Awab G,
    8. Barnadas C,
    9. Berry A,
    10. Boum Y,
    11. Bustos MD,
    12. Cao J,
    13. Chen JH,
    14. Collet L,
    15. Cui L,
    16. Thakur GD,
    17. Dieye A,
    18. Djallé D,
    19. Dorkenoo MA,
    20. Eboumbou-Moukoko CE,
    21. Espino FE,
    22. Fandeur T,
    23. Ferreira-da-Cruz MF,
    24. Fola AA,
    25. Fuehrer HP,
    26. Hassan AM,
    27. Herrera S,
    28. Hongvanthong B,
    29. Houzé S,
    30. Ibrahim ML,
    31. Jahirul-Karim M,
    32. Jiang L,
    33. Kano S,
    34. Ali-Khan W,
    35. Khanthavong M,
    36. Kremsner PG,
    37. Lacerda M,
    38. Leang R,
    39. Leelawong M,
    40. Li M,
    41. Lin K,
    42. Mazarati JB,
    43. Ménard S,
    44. Morlais I,
    45. Muhindo-Mavoko H,
    46. Musset L,
    47. Na-Bangchang K,
    48. Nambozi M,
    49. Niaré K,
    50. Noedl H
    , et al. 2016. A worldwide map of Plasmodium falciparum K13-propeller polymorphisms. N Engl J Med 374:2453–2464. doi:10.1056/NEJMoa1513137.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Kobasa T,
    2. Talundzic E,
    3. Sug-Aram R,
    4. Boondat P,
    5. Goldman IF,
    6. Lucchi NW,
    7. Dharmarak P,
    8. Sintasath D,
    9. Fukuda M,
    10. Whistler T,
    11. MacArthur J,
    12. Udhayakumar V,
    13. Prempree P,
    14. Chinanonwait N
    . 2018. Emergence and spread of kelch13 mutations associated with artemisinin resistance in Plasmodium falciparum parasites in 12 Thai provinces from 2007 to 2016. Antimicrob Agents Chemother 62:e02141-17. doi:10.1128/AAC.02141-17.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    WWARN K13 Genotype-Phenotype Study Group. 2019. Association of mutations in the Plasmodium falciparum Kelch13 gene (Pf3D7_1343700) with parasite clearance rates after artemisinin-based treatments-a WWARN individual patient data meta-analysis. BMC Med 17:1. doi:10.1186/s12916-018-1207-3.
    OpenUrlCrossRef
  20. 20.↵
    1. Plucinski MM,
    2. Dimbu PR,
    3. Macaia AP,
    4. Ferreira CM,
    5. Samutondo C,
    6. Quivinja J,
    7. Afonso M,
    8. Kiniffo R,
    9. Mbounga E,
    10. Kelley JS,
    11. Patel DS,
    12. He Y,
    13. Talundzic E,
    14. Garrett DO,
    15. Halsey ES,
    16. Udhayakumar V,
    17. Ringwald P,
    18. Fortes F
    . 2017. Efficacy of artemether-lumefantrine, artesunate-amodiaquine, and dihydroartemisinin-piperaquine for treatment of uncomplicated Plasmodium falciparum malaria in Angola, 2015. Malar J 16:62. doi:10.1186/s12936-017-1712-4.
    OpenUrlCrossRef
  21. 21.↵
    1. Madamet M,
    2. Kounta MB,
    3. Wade KA,
    4. Lo G,
    5. Diawara S,
    6. Fall M,
    7. Bercion R,
    8. Nakoulima A,
    9. Fall KB,
    10. Benoit N,
    11. Gueye MW,
    12. Fall B,
    13. Diatta B,
    14. Pradines B
    . 2017. Absence of association between polymorphisms in the K13 gene and the presence of Plasmodium falciparum parasites at day 3 after treatment with artemisinin derivatives in Senegal. Int J Antimicrob Agents 49:754–756. doi:10.1016/j.ijantimicag.2017.01.032.
    OpenUrlCrossRef
  22. 22.↵
    1. Gupta H,
    2. Macete E,
    3. Bulo H,
    4. Salvador C,
    5. Warsame M,
    6. Carvalho E,
    7. Ménard D,
    8. Ringwald P,
    9. Bassat Q,
    10. Enosse S,
    11. Mayor A
    . 2018. Drug-resistant polymorphisms and copy numbers in Plasmodium falciparum, Mozambique, 2015. Emerg Infect Dis 24:40–48. doi:10.3201/eid2401.170864.
    OpenUrlCrossRef
  23. 23.↵
    1. de Laurent ZR,
    2. Chebon LJ,
    3. Ingasia LA,
    4. Akala HM,
    5. Andagalu B,
    6. Ochola-Oyier LI,
    7. Kamau E
    . 2018. Polymorphisms in the K13 gene in Plasmodium falciparum from different malaria transmission areas of Kenya. Am J Trop Med Hyg 98:1360–1366. doi:10.4269/ajtmh.17-0505.
    OpenUrlCrossRef
  24. 24.↵
    1. Chenet SM,
    2. Akinyi Okoth S,
    3. Huber CS,
    4. Chandrabose J,
    5. Lucchi NW,
    6. Talundzic E,
    7. Krishnalall K,
    8. Ceron N,
    9. Musset L,
    10. Macedo de Oliveira A,
    11. Venkatesan M,
    12. Rahman R,
    13. Barnwell JW,
    14. Udhayakumar V
    . 2016. Independent emergence of the Plasmodium falciparum kelch propeller domain mutant allele C580Y in Guyana. J Infect Dis 213:1472–1475. doi:10.1093/infdis/jiv752.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Mita T,
    2. Culleton R,
    3. Takahashi N,
    4. Nakamura M,
    5. Tsukahara T,
    6. Hunja CW,
    7. Win ZZ,
    8. Htike WW,
    9. Marma AS,
    10. Dysoley L,
    11. Ndounga M,
    12. Dzodzomenyo M,
    13. Akhwale WS,
    14. Kobayashi J,
    15. Uemura H,
    16. Kaneko A,
    17. Hombhanje F,
    18. Ferreira MU,
    19. Björkman A,
    20. Endo H,
    21. Ohashi J
    . 2016. Little polymorphism at the K13 propeller locus in worldwide Plasmodium falciparum populations prior to the introduction of artemisinin combination therapies. Antimicrob Agents Chemother 60:3340–3347. doi:10.1128/AAC.02370-15.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Gomes LR,
    2. Lavigne A,
    3. Peterka CL,
    4. Brasil P,
    5. Ménard D,
    6. Daniel-Ribeiro CT,
    7. Ferreira-da-Cruz MF
    . 2018. Absence of K13 polymorphism in Plasmodium falciparum from Brazilian areas where the parasite is endemic. Antimicrob Agents Chemother 62:e00354-18. doi:10.1128/AAC.00354-18.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Otto TD,
    2. Gilabert A,
    3. Crellen T,
    4. Böhme U,
    5. Arnathau C,
    6. Sanders M,
    7. Oyola SO,
    8. Okouga AP,
    9. Boundenga L,
    10. Willaume E,
    11. Ngoubangoye B,
    12. Moukodoum ND,
    13. Paupy C,
    14. Durand P,
    15. Rougeron V,
    16. Ollomo B,
    17. Renaud F,
    18. Newbold C,
    19. Berriman M,
    20. Prugnolle F
    . 2018. Genomes of all known members of a Plasmodium subgenus reveal paths to virulent human malaria. Nat Microbiol 3:687–697. doi:10.1038/s41564-018-0162-2.
    OpenUrlCrossRef
  28. 28.↵
    1. Pacheco MA,
    2. Cranfield M,
    3. Cameron K,
    4. Escalante AA
    . 2013. Malarial parasite diversity in chimpanzees: the value of comparative approaches to ascertain the evolution of Plasmodium falciparum antigens. Malar J 12:328. doi:10.1186/1475-2875-12-328.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Pacheco MA,
    2. Matta NE,
    3. Valkiunas G,
    4. Parker PG,
    5. Mello B,
    6. Stanley CE, Jr,
    7. Lentino M,
    8. Garcia-Amado MA,
    9. Cranfield M,
    10. Kosakovsky Pond SL,
    11. Escalante AA
    . 2018. Mode and rate of evolution of haemosporidian mitochondrial genomes: timing the radiation of avian parasites. Mol Biol Evol 35:383–403. doi:10.1093/molbev/msx285.
    OpenUrlCrossRef
  30. 30.↵
    1. Tamura K,
    2. Battistuzzi FU,
    3. Billing-Ross P,
    4. Murillo O,
    5. Filipski A,
    6. Kumar S
    . 2012. Estimating divergence times in large molecular phylogenies. Proc Natl Acad Sci U S A 109:19333–19338. doi:10.1073/pnas.1213199109.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Hayakawa T,
    2. Culleton R,
    3. Otani H,
    4. Horii T,
    5. Tanabe K
    . 2008. Big bang in the evolution of extant malaria parasites. Mol Biol Evol 25:2233–2239. doi:10.1093/molbev/msn171.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    1. Muehlenbein MP,
    2. Pacheco MA,
    3. Taylor JE,
    4. Prall SP,
    5. Ambu L,
    6. Nathan S,
    7. Alsisto S,
    8. Ramirez D,
    9. Escalante AA
    . 2015. Accelerated diversification of nonhuman primate malarias in Southeast Asia: adaptive radiation or geographic speciation? Mol Biol Evol 32:422–439. doi:10.1093/molbev/msu310.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Seo TK,
    2. Kishino H,
    3. Thorne JL
    . 2004. Estimating absolute rates of synonymous and non-synonymous nucleotide substitution in order to characterize natural selection and date species divergences. Mol Biol Evol 217:1201–1213. doi:10.1093/molbev/msh088.
    OpenUrlCrossRef
  34. 34.↵
    1. Pond SLK,
    2. Frost SDW,
    3. Muse SV
    . 2005. HyPhy: hypothesis testing using phylogenies. Bioinformatics 215:676–679. doi:10.1093/bioinformatics/bti079.
    OpenUrlCrossRef
  35. 35.↵
    1. Murrell B,
    2. Wertheim JO,
    3. Moola S,
    4. Weighill T,
    5. Scheffler K,
    6. Kosakovsky Pond SL
    . 2012. Detecting individual sites subject to episodic diversifying selection. PLoS Genet 87:e1002764. doi:10.1371/journal.pgen.1002764.
    OpenUrlCrossRef
  36. 36.↵
    1. Murrell B,
    2. Weaver S,
    3. Smith MD,
    4. Wertheim JO,
    5. Murrell S,
    6. Aylward A,
    7. Eren K,
    8. Pollner T,
    9. Martin DP,
    10. Smith DM,
    11. Scheffler K,
    12. Kosakovsky Pond SL
    . 2015. Gene-wide identification of episodic selection. Mol Biol Evol 32:1365–1371. doi:10.1093/molbev/msv035.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Wertheim JO,
    2. Murrell B,
    3. Smith MD,
    4. Kosakovsky Pond SL,
    5. Scheffler K
    . 2015. RELAX: detecting relaxed selection in a phylogenetic framework. Mol Biol Evol 323:820–832. doi:10.1093/molbev/msu400.
    OpenUrlCrossRef
  38. 38.↵
    1. McCollum AM,
    2. Poe AC,
    3. Hamel M,
    4. Huber C,
    5. Zhou Z,
    6. Shi YP,
    7. Ouma P,
    8. Vulule J,
    9. Bloland P,
    10. Slutsker L,
    11. Barnwell JW,
    12. Udhayakumar V,
    13. Escalante AA
    . 2006. Antifolate resistance in Plasmodium falciparum: multiple origins and identification of novel dhfr alleles. J Infect Dis 194:189–197. doi:10.1086/504687.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    1. McCollum AM,
    2. Schneider KA,
    3. Griffing SM,
    4. Zhou Z,
    5. Kariuki S,
    6. Ter-Kuile F,
    7. Shi YP,
    8. Slutsker L,
    9. Lal AA,
    10. Udhayakumar V,
    11. Escalante AA
    . 2012. Differences in selective pressure on dhps and dhfr drug resistant mutations in western Kenya. Malar J 11:77. doi:10.1186/1475-2875-11-77.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Anderson TJ,
    2. Nair S,
    3. McDew-White M,
    4. Cheeseman IH,
    5. Nkhoma S,
    6. Bilgic F,
    7. McGready R,
    8. Ashley E,
    9. Pyae Phyo A,
    10. White NJ,
    11. Nosten F
    . 2017. Population parameters underlying an ongoing soft sweep in Southeast Asian malaria parasites. Mol Biol Evol 34:131–144. doi:10.1093/molbev/msw228.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Tanabe K,
    2. Mita T,
    3. Jombart T,
    4. Eriksson A,
    5. Horibe S,
    6. Palacpac N,
    7. Ranford-Cartwright L,
    8. Sawai H,
    9. Sakihama N,
    10. Ohmae H,
    11. Nakamura M,
    12. Ferreira MU,
    13. Escalante AA,
    14. Prugnolle F,
    15. Björkman A,
    16. Färnert A,
    17. Kaneko A,
    18. Horii T,
    19. Manica A,
    20. Kishino H,
    21. Balloux F
    . 2010. Plasmodium falciparum accompanied the human expansion out of Africa. Curr Biol 20:1283–1289. doi:10.1016/j.cub.2010.05.053.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    1. Gazave E,
    2. Chang D,
    3. Clark AG,
    4. Keinan A
    . 2013. Population growth inflates the per-individual number of deleterious mutations and reduces their mean effect. Genetics 195:969–978. doi:10.1534/genetics.113.153973.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. Gravel S
    . 2016. When is selection effective? Genetics 203:451–462. doi:10.1534/genetics.115.184630.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Taylor JE,
    2. Pacheco MA,
    3. Bacon DJ,
    4. Beg MA,
    5. Machado RL,
    6. Fairhurst RM,
    7. Herrera S,
    8. Kim JY,
    9. Menard D,
    10. Póvoa MM,
    11. Villegas L,
    12. Mulyanto Snounou G,
    13. Cui L,
    14. Zeyrek FY,
    15. Escalante AA
    . 2013. The evolutionary history of Plasmodium vivax as inferred from mitochondrial genomes: parasite genetic diversity in the Americas. Mol Biol Evol 30:2050–2064. doi:10.1093/molbev/mst104.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    1. Hupalo DN,
    2. Luo Z,
    3. Melnikov A,
    4. Sutton PL,
    5. Rogov P,
    6. Escalante A,
    7. Vallejo AF,
    8. Herrera S,
    9. Arévalo-Herrera M,
    10. Fan Q,
    11. Wang Y,
    12. Cui L,
    13. Lucas CM,
    14. Durand S,
    15. Sanchez JF,
    16. Baldeviano GC,
    17. Lescano AG,
    18. Laman M,
    19. Barnadas C,
    20. Barry A,
    21. Mueller I,
    22. Kazura JW,
    23. Eapen A,
    24. Kanagaraj D,
    25. Valecha N,
    26. Ferreira MU,
    27. Roobsoong W,
    28. Nguitragool W,
    29. Sattabonkot J,
    30. Gamboa D,
    31. Kosek M,
    32. Vinetz JM,
    33. González-Cerón L,
    34. Birren BW,
    35. Neafsey DE,
    36. Carlton JM
    . 2016. Population genomics studies identify signatures of global dispersal and drug resistance in Plasmodium vivax. Nat Genet 48:953–958. doi:10.1038/ng.3588.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Kim Y,
    2. Escalante AA,
    3. Schneider KA
    . 2014. A population genetic model for the initial spread of partially resistant malaria parasites under anti-malarial combination therapy and weak intrahost competition. PLoS One 9:e101601. doi:10.1371/journal.pone.0101601.
    OpenUrlCrossRef
  47. 47.↵
    1. Straimer J,
    2. Gnädig NF,
    3. Stokes BH,
    4. Ehrenberger M,
    5. Crane AA,
    6. Fidock DA
    . 2017. Plasmodium falciparum K13 mutations differentially impact ozonide susceptibility and parasite fitness in vitro. mBio 8:e00172-17. doi:10.1128/mBio.00172-17.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Nair S,
    2. Li X,
    3. Arya GA,
    4. McDew-White M,
    5. Ferrari M,
    6. Nosten F,
    7. Anderson TJC
    . 2018. Fitness costs and the rapid spread of kelch13-C580Y substitutions conferring artemisinin resistance. Antimicrob Agents Chemother 62:e00605-18. doi:10.1128/AAC.00605-18.
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    1. Ataide R,
    2. Ashley EA,
    3. Powell R,
    4. Chan JA,
    5. Malloy MJ,
    6. O’Flaherty K,
    7. Takashima E,
    8. Langer C,
    9. Tsuboi T,
    10. Dondorp AM,
    11. Day NP,
    12. Dhorda M,
    13. Fairhurst RM,
    14. Lim P,
    15. Amaratunga C,
    16. Pukrittayakamee S,
    17. Hien TT,
    18. Htut Y,
    19. Mayxay M,
    20. Faiz MA,
    21. Beeson JG,
    22. Nosten F,
    23. Simpson JA,
    24. White NJ,
    25. Fowkes FJ
    . 2017. Host immunity to Plasmodium falciparum and the assessment of emerging artemisinin resistance in a multinational cohort. Proc Natl Acad Sci U S A 114:3515–3520. doi:10.1073/pnas.1615875114.
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    1. Klein EY,
    2. Smith DL,
    3. Boni MF,
    4. Laxminarayan R
    . 2008. Clinically immune hosts as a refuge for drug-sensitive malaria parasites. Malar J 7:67. doi:10.1186/1475-2875-7-67.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Ghani AC
    . 2018. Can improving access to care help to eliminate malaria? Lancet 391:1870–1871. doi:10.1016/S0140-6736(18)30910-3.
    OpenUrlCrossRef
  52. 52.↵
    1. Dwivedi A,
    2. Khim N,
    3. Reynes C,
    4. Ravel P,
    5. Ma L,
    6. Tichit M,
    7. Bourchier C,
    8. Kim S,
    9. Dourng D,
    10. Khean C,
    11. Chim P,
    12. Siv S,
    13. Frutos R,
    14. Lek D,
    15. Mercereau-Puijalon O,
    16. Ariey F,
    17. Menard D,
    18. Cornillot E
    . 2016. Plasmodium falciparum parasite population structure and gene flow associated to anti-malarial drugs resistance in Cambodia. Malar J 15:319. doi:10.1186/s12936-016-1370-y.
    OpenUrlCrossRef
  53. 53.↵
    1. Kheang ST,
    2. Sovannaroth S,
    3. Ek S,
    4. Chy S,
    5. Chhun P,
    6. Mao S,
    7. Nguon S,
    8. Lek DS,
    9. Menard D,
    10. Kak N
    . 2017. Prevalence of K13 mutation and day-3 positive parasitaemia in artemisinin-resistant malaria endemic area of Cambodia: a cross-sectional study. Malar J 16:372. doi:10.1186/s12936-017-2024-4.
    OpenUrlCrossRef
  54. 54.↵
    1. Coatney RG,
    2. Collins WE,
    3. Warren M,
    4. Contacos PG
    . 1971. The primate malarias. U.S. Government Printing Office, Washington, DC.
  55. 55.↵
    1. Aurrecoechea C,
    2. Brestelli J,
    3. Brunk BP,
    4. Dommer J,
    5. Fischer S,
    6. Gajria B,
    7. Gao X,
    8. Gingle A,
    9. Grant G,
    10. Harb OS,
    11. Heiges M,
    12. Innamorato F,
    13. Iodice J,
    14. Kissinger JC,
    15. Kraemer E,
    16. Li W,
    17. Miller JA,
    18. Nayak V,
    19. Pennington C,
    20. Pinney DF,
    21. Roos DS,
    22. Ross C,
    23. Stoeckert CJ, Jr,
    24. Treatman C,
    25. Wang H
    . 2009. PlasmoDB: a functional genomic database for malaria parasites. Nucleic Acids Res 37:D539–D543. doi:10.1093/nar/gkn814.
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.↵
    1. Bensch S,
    2. Hellgren O,
    3. Pérez-Tris J
    . 2009. MalAvi: a public database of malaria parasites and related haemosporidians in avian hosts based on mitochondrial cytochrome b lineages. Mol Ecol Resour 95:1353–1358. doi:10.1111/j.1755-0998.2009.02692.x.
    OpenUrlCrossRef
  57. 57.↵
    1. Benson DA,
    2. Cavanaugh M,
    3. Clark K,
    4. Karsch-Mizrachi I,
    5. Lipman DJ,
    6. Ostell J,
    7. Sayers EW
    . 2012. 2013. GenBank. Nucleic Acids Res 41:D36–D42. doi:10.1093/nar/gks1195.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Talundzic E,
    2. Chenet SM,
    3. Goldman IF,
    4. Patel DS,
    5. Nelson JA,
    6. Plucinski MM,
    7. Barnwell JW,
    8. Udhayakumar V
    . 2015. Genetic analysis and species specific amplification of the artemisinin resistance-associated kelch propeller domain in P. falciparum and P. vivax. PLoS One 10:e0136099. doi:10.1371/journal.pone.0136099.
    OpenUrlCrossRef
  59. 59.↵
    1. Altschul SF,
    2. Madden TL,
    3. Schäffer AA,
    4. Zhang J,
    5. Zhang Z,
    6. Miller W,
    7. Lipman DJ
    . 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25:3389–3402. doi:10.1093/nar/25.17.3389.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Talundzic E,
    2. Ravishankar S,
    3. Kelley J,
    4. Patel D,
    5. Plucinski M,
    6. Schmedes S,
    7. Ljolje D,
    8. Clemons B,
    9. Madison-Antenucci S,
    10. Arguin PM,
    11. Lucchi NW,
    12. Vannberg F,
    13. Udhayakumar V
    . 2018. Next-generation sequencing and bioinformatics protocol for malaria drug resistance marker surveillance. Antimicrob Agents Chemother 62:e02474-17. doi:10.1128/AAC.02474-17.
    OpenUrlAbstract/FREE Full Text
  61. 61.↵
    1. Gouy M,
    2. Guindon S,
    3. Gascuel O
    . 2010. SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol 272:221–224. doi:10.1093/molbev/msp259.
    OpenUrlCrossRef
  62. 62.↵
    1. Ronquist F,
    2. Huelsenbeck JP
    . 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 1912:1572–1574. doi:10.1093/bioinformatics/btg180.
    OpenUrlCrossRef
  63. 63.↵
    1. Kumar S,
    2. Stecher G,
    3. Tamura K
    . 2016. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol 337:1870–1874. doi:10.1093/molbev/msw054.
    OpenUrlCrossRefPubMed
  64. 64.↵
    1. Rambaut A
    . 2006. FigTree tree figure drawing tool version 1.3.1. Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom. http://tree.bio.ed.ac.uk/software/figtree/. Accessed 31 January 2018.
  65. 65.↵
    1. Kumar S,
    2. Stecher G,
    3. Li M,
    4. Knyaz C,
    5. Tamura K
    . 2018. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol 35:1547–1549. doi:10.1093/molbev/msy096.
    OpenUrlCrossRefPubMed
  66. 66.↵
    1. Bouckaert R,
    2. Heled J,
    3. Kühnert D,
    4. Vaughan T,
    5. Wu CH,
    6. Xie D,
    7. Suchard MA,
    8. Rambaut A,
    9. Drummond AJ
    . 2014. BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput Biol 10:e1003537. doi:10.1371/journal.pcbi.1003537.
    OpenUrlCrossRefPubMed
  67. 67.↵
    1. Drummond AJ,
    2. Ho SY,
    3. Phillips MJ,
    4. Rambaut A
    . 2006. Relaxed phylogenetics and dating with confidence. PLoS Biol 4:e88. doi:10.1371/journal.pbio.0040088.
    OpenUrlCrossRefPubMed
  68. 68.↵
    1. Mu J,
    2. Joy DA,
    3. Duan J,
    4. Huang Y,
    5. Carlton J,
    6. Walker J,
    7. Barnwell J,
    8. Beerli P,
    9. Charleston MA,
    10. Pybus OG,
    11. Su XZ
    . 2005. Host switch leads to emergence of Plasmodium vivax malaria in humans. Mol Biol Evol 22:1686–1693. doi:10.1093/molbev/msi160.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    1. Pacheco MA,
    2. Battistuzzi FU,
    3. Junge RE,
    4. Cornejo OE,
    5. Williams CV,
    6. Landau I,
    7. Rabetafika L,
    8. Snounou G,
    9. Jones-Engel L,
    10. Escalante AA
    . 2011. Timing the origin of human malarias: the lemur puzzle. BMC Evol Biol 11:299. doi:10.1186/1471-2148-11-299.
    OpenUrlCrossRefPubMed
  70. 70.↵
    1. Benton MJ,
    2. Donoghue PC
    . 2007. Paleontological evidence to date the tree of life. Mol Biol Evol 241:26–53. doi:10.1093/molbev/msl150.
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    1. Perelman P,
    2. Johnson WE,
    3. Roos C,
    4. Seuánez HN,
    5. Horvath JE,
    6. Moreira MA,
    7. Kessing B,
    8. Pontius J,
    9. Roelke M,
    10. Rumpler Y,
    11. Schneider MP,
    12. Silva A,
    13. O’Brien SJ,
    14. Pecon-Slattery J
    . 2011. A molecular phylogeny of living primates. PLoS Genet 73:e1001342. doi:10.1371/journal.pgen.1001342.
    OpenUrlCrossRef
  72. 72.↵
    1. Nei M,
    2. Kumar S
    . 2000. Molecular evolution and phylogenetics. Oxford University Press, New York, NY.
  73. 73.↵
    1. Librado P,
    2. Rozas J
    . 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25:1451–1452. doi:10.1093/bioinformatics/btp187.
    OpenUrlCrossRefPubMedWeb of Science
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
Evolution and Genetic Diversity of the k13 Gene Associated with Artemisinin Delayed Parasite Clearance in Plasmodium falciparum
M. Andreína Pacheco, Esha R. Kadakia, Zainab Chaudhary, Douglas J. Perkins, Julia Kelley, Shashidhar Ravishankar, Michael Cranfield, Eldin Talundzic, Venkatachalam Udhayakumar, Ananias A. Escalante
Antimicrobial Agents and Chemotherapy Jul 2019, 63 (8) e02550-18; DOI: 10.1128/AAC.02550-18

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Antimicrobial Agents and Chemotherapy article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Evolution and Genetic Diversity of the k13 Gene Associated with Artemisinin Delayed Parasite Clearance in Plasmodium falciparum
(Your Name) has forwarded a page to you from Antimicrobial Agents and Chemotherapy
(Your Name) thought you would be interested in this article in Antimicrobial Agents and Chemotherapy.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Evolution and Genetic Diversity of the k13 Gene Associated with Artemisinin Delayed Parasite Clearance in Plasmodium falciparum
M. Andreína Pacheco, Esha R. Kadakia, Zainab Chaudhary, Douglas J. Perkins, Julia Kelley, Shashidhar Ravishankar, Michael Cranfield, Eldin Talundzic, Venkatachalam Udhayakumar, Ananias A. Escalante
Antimicrobial Agents and Chemotherapy Jul 2019, 63 (8) e02550-18; DOI: 10.1128/AAC.02550-18
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Plasmodium falciparum
Plasmodium vivax
artemisinin resistance
k13 gene
kelch propeller domain
rates of evolution

Related Articles

Cited By...

About

  • About AAC
  • Editor in Chief
  • Editorial Board
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • AAC Podcast
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Ethics
  • Contact Us

Follow #AACJournal

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Print ISSN: 0066-4804; Online ISSN: 1098-6596