Integration of population and functional genomics to understand mechanisms of artemisinin resistance in Plasmodium falciparum

Molecular characterization of the Plasmodium falciparum k13 gene helix domain in samples from native communities of Condorcanqui, Amazonas, Perú

Introduction. Resistance of Plasmodium falciparum to different antimalarial drugs is an obstacle to disease elimination. The artemisinin-resistant genotype of P. falciparum can be assessed by examining polymorphisms in the helix domain of the Pfk13 gene. The World Health Organization recommends these mutations as molecular markers to detect artemisinin-resistant in countries where P. falciparum malaria is endemic. Objective. To identify artemisinin resistance-related mutations present in the helix domain of the P. falciparum k13 gene. Materials and methods. We collected a total of 51 samples through passive case detection, positive for Plasmodium by microscopy, from six communities in the district of Río Santiago in Condorcanqui, Amazonas. Molecular species confirmation was performed by real-time PCR and Pfk13 helix domain was amplified and sequenced by capillary electrophoresis. The obtained sequences were compared with the wild type 3D7 reference strain. Results. A total of 51 positive samples were confirmed for P. falciparum from the communities of Ayambis, Chapiza, Palometa, Muchinguis, Alianza Progreso and Caterpiza. DNA sequences alignment showed the absence of resistance-associated mutations in the k13 gene of the collected samples. Discussion. The obtained results are consistent with similar studies conducted in other South American countries, including Perú, so these data provide a baseline for artemisinin-resistance molecular surveillance in the Amazon region and reinforce the efficacy of artemisinin-based combination therapy in this area.

A novel Modulator of Ring Stage Translation (MRST) gene alters artemisinin sensitivity in Plasmodium falciparum

The implementation of artemisinin (ART) combination therapies (ACTs) has greatly decreased deaths caused by Plasmodium falciparum malaria, but increasing ACT resistance in Southeast Asia and Africa could reverse this progress. Parasite population genetic studies have identified numerous genes, single-nucleotide polymorphisms (SNPs), and transcriptional signatures associated with altered artemisinin activity with SNPs in the Kelch13 (K13) gene being the most well-characterized artemisinin resistance marker. However, there is an increasing evidence that resistance to artemisinin in P. falciparum is not related only to K13 SNPs, prompting the need to characterize other novel genes that can alter ART responses in P. falciparum. In our previous analyses of P. falciparum piggyBac mutants, several genes of unknown function exhibited increased sensitivity to artemisinin that was similar to a mutant of K13. Further analysis of these genes and their gene co-expression networks indicated that the ART sensitivity cluster was functionally linked to DNA replication and repair, stress responses, and maintenance of homeostatic nuclear activity. In this study, we have characterized PF3D7_1136600, another member of the ART sensitivity cluster. Previously annotated as a conserved Plasmodium gene of unknown function, we now provide putative annotation of this gene as a Modulator of Ring Stage Translation (MRST). Our findings reveal that the mutagenesis of MRST affects gene expression of multiple translation-associated pathways during the early ring stage of asexual development via putative ribosome assembly and maturation activity, suggesting an essential role of MRST in protein biosynthesis and another novel mechanism of altering the parasite's ART drug response.IMPORTANCEPlasmodium falciparum malaria killed more than 600,000 people in 2021, though ACTs have been critical in reducing malaria mortality as a first-line treatment for infection. However, ACT resistance in Southeast Asia and emerging resistance in Africa are detrimental to this progress. Mutations to Kelch13 (K13) have been identified to confer increased artemisinin tolerance in field isolates, however, genes other than K13 are implicated in altering how the parasite responds to artemisinin prompts additional analysis. Therefore, in this study we have characterized a P. falciparum mutant clone with altered sensitivity to artemisinin and identified a novel gene (PF3D7_1136600) that is associated with alterations to parasite translational metabolism during critical timepoints for artemisinin drug response. Many genes of the P. falciparum genome remain unannotated, posing a challenge for drug-gene characterizations in the parasite. Therefore, through this study, we have putatively annotated PF3D7_1136600 as a novel MRST gene and have identified a potential link between MRST and parasite stress response mechanisms.

Chemogenomic Profiling of a Plasmodium falciparum Transposon Mutant Library Reveals Shared Effects of Dihydroartemisinin and Bortezomib on Lipid Metabolism and Exported Proteins

The antimalarial activity of the frontline drug artemisinin involves generation of reactive oxygen species (ROS) leading to oxidative damage of parasite proteins. To achieve homeostasis and maintain protein quality control in the overwhelmed parasite, the ubiquitin-proteasome system kicks in. Even though molecular markers for artemisinin resistance like pfkelch13 have been identified, the intricate network of mechanisms driving resistance remains to be elucidated. Here, we report a forward genetic screening strategy that enables a broader identification of genetic factors responsible for altering sensitivity to dihydroartemisinin (DHA) and a proteasome inhibitor, bortezomib (BTZ). Using a library of isogenic piggyBac mutants in P. falciparum, we defined phenotype-genotype associations influencing drug responses and highlighted shared mechanisms between the two processes, which mainly included proteasome-mediated degradation and the lipid metabolism genes. Additional transcriptomic analysis of a DHA/BTZ-sensitive piggyBac mutant showed it is possible to find differences between the two response mechanisms on the specific components for regulation of the exportome. Our results provide further insight into the molecular mechanisms of antimalarial drug resistance. IMPORTANCE Malaria control is seriously threatened by the emergence and spread of Plasmodium falciparum resistance to the leading antimalarial, artemisinin. The potent killing activity of artemisinin results from oxidative damage unleashed by free heme activation released by hemoglobin digestion. Although the ubiquitin-proteasome system is considered critical for parasite survival of this toxicity, the diverse genetic changes linked to artemisinin resistance are complex and, so far, have not included the ubiquitin-proteasome system. In this study, we use a systematic forward genetic approach by screening a library of P. falciparum random piggyBac mutants to decipher the genetic factors driving malaria parasite responses to the oxidative stress caused by antimalarial drugs. This study compares phenotype-genotype associations influencing dihydroartemisinin responses with the proteasome inhibitor bortezomib to delineate the role of ubiquitin-proteasome system. Our study highlights shared and unique pathways from the complex array of molecular processes critical for P. falciparum survival resulting from the oxidative damage of artemisinin.

Malaria Molecular Surveillance in the Peruvian Amazon with a Novel Highly Multiplexed Plasmodium falciparum AmpliSeq Assay

Molecular surveillance for malaria has great potential to support national malaria control programs (NMCPs). To bridge the gap between research and implementation, several applications (use cases) have been identified to align research, technology development, and public health efforts. For implementation at NMCPs, there is an urgent need for feasible and cost-effective tools. We designed a new highly multiplexed deep sequencing assay (Pf AmpliSeq), which is compatible with benchtop sequencers, that allows high-accuracy sequencing with higher coverage and lower cost than whole-genome sequencing (WGS), targeting genomic regions of interest. The novelty of the assay is its high number of targets multiplexed into one easy workflow, combining population genetic markers with 13 nearly full-length resistance genes, which is applicable for many different use cases. We provide the first proof of principle for hrp2 and hrp3 deletion detection using amplicon sequencing. Initial sequence data processing can be performed automatically, and subsequent variant analysis requires minimal bioinformatic skills using any tabulated data analysis program. The assay was validated using a retrospective sample collection (n = 254) from the Peruvian Amazon between 2003 and 2018. By combining phenotypic markers and a within-country 28-single-nucleotide-polymorphism (SNP) barcode, we were able to distinguish different lineages with multiple resistance haplotypes (in dhfr, dhps, crt and mdr1) and hrp2 and hrp3 deletions, which have been increasing in recent years. We found no evidence to suggest the emergence of artemisinin (ART) resistance in Peru. These findings indicate a parasite population that is under drug pressure but is susceptible to current antimalarials and demonstrate the added value of a highly multiplexed molecular tool to inform malaria strategies and surveillance systems. IMPORTANCE While the power of next-generation sequencing technologies to inform and guide malaria control programs has become broadly recognized, the integration of genomic data for operational incorporation into malaria surveillance remains a challenge in most countries where malaria is endemic. The main obstacles include limited infrastructure, limited access to high-throughput sequencing facilities, and the need for local capacity to run an in-country analysis of genomes at a large-enough scale to be informative for surveillance. In addition, there is a lack of standardized laboratory protocols and automated analysis pipelines to generate reproducible and timely results useful for relevant stakeholders. With our standardized laboratory and bioinformatic workflow, malaria genetic surveillance data can be readily generated by surveillance researchers and malaria control programs in countries of endemicity, increasing ownership and ensuring timely results for informed decision- and policy-making.

Short tandem repeat polymorphism in the promoter region of cyclophilin 19B drives its transcriptional upregulation and contributes to drug resistance in the malaria parasite Plasmodium falciparum

Resistance of the human malaria parasites, Plasmodium falciparum, to artemisinins is now fully established in Southeast Asia and is gradually emerging in Sub-Saharan Africa. Although nonsynonymous SNPs in the pfk13 Kelch-repeat propeller (KREP) domain are clearly associated with artemisinin resistance, their functional relevance requires cooperation with other genetic factors/alterations of the P. falciparum genome, collectively referred to as genetic background. Here we provide experimental evidence that P. falciparum cyclophilin 19B (PfCYP19B) may represent one putative factor in this genetic background, contributing to artemisinin resistance via its increased expression. We show that overexpression of PfCYP19B in vitro drives limited but significant resistance to not only artemisinin but also piperaquine, an important partner drug in artemisinin-based combination therapies. We showed that PfCYP19B acts as a negative regulator of the integrated stress response (ISR) pathway by modulating levels of phosphorylated eIF2α (eIF2α-P). Curiously, artemisinin and piperaquine affect eIF2α-P in an inverse direction that in both cases can be modulated by PfCYP19B towards resistance. Here we also provide evidence that the upregulation of PfCYP19B in the drug-resistant parasites appears to be maintained by a short tandem repeat (SRT) sequence polymorphism in the gene's promoter region. These results support a model that artemisinin (and other drugs) resistance mechanisms are complex genetic traits being contributed to by altered expression of multiple genes driven by genetic polymorphism at their promoter regions.

Plasmodium falciparum rosetting protects schizonts against artemisinin

Background

Artemisinin (ART) resistance in Plasmodium falciparum is thought to occur during the early stage of the parasite's erythrocytic cycle. Here, we identify a novel factor associated with the late stage parasite development that contributes to ART resistance.

Methods

Rosetting rates of clinical isolates pre- and post- brief (one hour) exposure to artesunate (AS, an ART derivative) were evaluated. The effects of AS-mediated rosetting on the post-AS-exposed parasite's replication and survival, as well as the extent of protection by AS-mediated rosetting on different parasite stages were investigated. The rosetting ligands, mechanisms, and gene mutations involved were studied.

Findings

Brief AS exposure stimulated rosetting, with AS-resistant isolates forming more rosettes in a more rapid manner. AS-mediated rosetting enabled infected erythrocytes (IRBC) to withstand AS exposure for several hours and protected the IRBC from phagocytosis. When their rosetting ability was blocked experimentally, the post-AS exposure survival advantage by the AS-resistant parasites was abrogated. Deletions in two genes coding for PfEMP1 exon 2 (PF3D7_0200300 and PF3D7_0223300) were found to be associated with AS-mediated rosetting, and these mutations were significantly selected through time in the parasite population under study, along with the K13 mutations, a molecular marker of ART-resistance.

Interpretation

Rapid ART parasite clearance is driven by the direct oxidative damages on IRBC by ART and the phagocytic destruction of the damaged IRBC. Rosetting serves as a rapid ‘buying time’ strategy that allows more parasites to complete schizont maturation, reinvasion and subsequent development into the intrinsically less ART-susceptible ring stage.

Funding

A*STAR, NMRC-OF-YIRG, HRC e-ASIA, Wellcome.

The apicoplast link to fever-survival and artemisinin-resistance in the malaria parasite

The emergence and spread of Plasmodium falciparum parasites resistant to front-line antimalarial artemisinin-combination therapies (ACT) threatens to erase the considerable gains against the disease of the last decade. Here, we develop a large-scale phenotypic screening pipeline and use it to carry out a large-scale forward-genetic phenotype screen in P. falciparum to identify genes allowing parasites to survive febrile temperatures. Screening identifies more than 200 P. falciparum mutants with differential responses to increased temperature. These mutants are more likely to be sensitive to artemisinin derivatives as well as to heightened oxidative stress. Major processes critical for P. falciparum tolerance to febrile temperatures and artemisinin include highly essential, conserved pathways associated with protein-folding, heat shock and proteasome-mediated degradation, and unexpectedly, isoprenoid biosynthesis, which originated from the ancestral genome of the parasite’s algal endosymbiont-derived plastid, the apicoplast. Apicoplast-targeted genes in general are upregulated in response to heat shock, as are other Plasmodium genes with orthologs in plant and algal genomes. Plasmodium falciparum parasites appear to exploit their innate febrile-response mechanisms to mediate resistance to artemisinin. Both responses depend on endosymbiont-derived genes in the parasite’s genome, suggesting a link to the evolutionary origins of Plasmodium parasites in free-living ancestors.

Repeating fever is a hallmark of malaria. Here, a large-scale forward genetic screen in malaria-causing Plasmodium falciparum identifies genes associated with parasite tolerance to host fever, including apicoplast targeted isoprenoid biosynthesis—sharing features with artemisinin resistance.