The origins of malaria artemisinin resistance defined by a genetic and transcriptomic background

Plasmodium falciparum artemisinin resistance: something gained in translation

Small-Saunders et al. uncovered a new facet of artemisinin resistance in Plasmodium in which parasites use a previously underexplored arm of stress response mechanisms. Through altered epitranscriptomic modifications on tRNA, changed translation patterns adapt resistant cells to facilitate entry into a quiescent-like state which provides the parasite an escape from many drugs.

Molecular markers of artemisinin resistance during falciparum malaria elimination in Eastern Myanmar

BACKGROUND

Artemisinin resistance in Plasmodium falciparum threatens global malaria elimination efforts. To contain and then eliminate artemisinin resistance in Eastern Myanmar a network of community-based malaria posts was instituted and targeted mass drug administration (MDA) with dihydroartemisinin-piperaquine (three rounds at monthly intervals) was conducted. The prevalence of artemisinin resistance during the elimination campaign (2013-2019) was characterized.

METHODS

Throughout the six-year campaign Plasmodium falciparum positive blood samples from symptomatic patients and from cross-sectional surveys were genotyped for mutations in kelch-13-a molecular marker of artemisinin resistance.

RESULT

The program resulted in near elimination of falciparum malaria. Of 5162 P. falciparum positive blood samples genotyped, 3281 (63.6%) had K13 mutations. The prevalence of K13 mutations was 73.9% in 2013 and 64.4% in 2019. Overall, there was a small but significant decline in the proportion of K13 mutants (p < 0.001). In the MDA villages there was no significant change in the K13 proportions before and after MDA. The distribution of different K13 mutations changed substantially; F446I and P441L mutations increased in both MDA and non-MDA villages, while most other K13 mutations decreased. The proportion of C580Y mutations fell from 9.2% (43/467) before MDA to 2.3% (19/813) after MDA (p < 0.001). Similar changes occurred in the 487 villages where MDA was not conducted.

CONCLUSION

The malaria elimination program in Kayin state, eastern Myanmar, led to a substantial reduction in falciparum malaria. Despite the intense use of artemisinin-based combination therapies, both in treatment and MDA, this did not select for artemisinin resistance.

What exactly does the PfK13 C580Y mutation in Plasmodium falciparum influence?

BACKGROUND

The emergence and spread of artemisinin resistance threaten global malaria control and elimination goals, and encourage research on the mechanisms of drug resistance in malaria parasites. Mutations in Plasmodium falciparum Kelch 13 (PfK13) protein are associated with artemisinin resistance, but the unique or common mechanism which results in this resistance is unclear.

METHODS

We analyzed the effects of the PfK13 mutation on the transcriptome and proteome of P. falciparum at different developmental stages. Additionally, the number of merozoites, hemozoin amount, and growth of P. falciparum 3D7C580Y and P. falciparum 3D7WT were compared. The impact of iron supplementation on the number of merozoites of P. falciparum 3D7C580Y was also examined.

RESULTS

We found that the PfK13 mutation did not significantly change glycolysis, TCA, pentose phosphate pathway, or oxidative phosphorylation, but did reduce the expression of reproduction- and DNA synthesis-related genes. The reduced number of merozoites, decreased level of hemozoin, and slowed growth of P. falciparum 3D7C580Y were consistent with these changes. Furthermore, adding iron supply could increase the number of the merozoites of P. falciparum 3D7C580Y.

CONCLUSIONS

These results revealed that the PfK13 mutation reduced hemoglobin ingestion, leading to artemisinin resistance, likely by decreasing the parasites' requirement for haem and iron. This study helps elucidate the mechanism of artemisinin resistance due to PfK13 mutations.

Genetic characteristics of P. falciparum parasites collected from 2012 to 2016 and anti-malaria resistance along the China-Myanmar border

BACKGROUNDS

The therapeutic efficacy studies of DHA-PIP for uncomplicated Plasmodium falciparum patients were implemented from 2012 to 2016 along China (Yunnan province)-Myanmar border, which verified the high efficacy of DHA-PIP. With the samples collected in these studies, the genetic characteristics of P. falciparum parasites based on in vivo parasite clearance time (PCT) was investigated to explore if these parasites had developed resistance to DHA and PIP at molecular level.

METHODS

The genetic characteristics were investigated based on K13 genotypes, copy numbers of genes pfpm2 and pfmdr1, and nine microsatellite loci (Short Tandem Repeats, STR) flanking the K13 gene on chromosome 13. The PCT 50s were compared based on different K13 genotypes, sites, periods and copy numbers.

RESULTS

In the NW (North-West Yunnan province bordering with Myanmar) region, F446I was the main K13 genotype. No significant differences for PCT 50s presented among three K13 genotypes. In SW (South-West Yunnan province bordering with Myanmar) region, only wild K13 genotype was detected in all parasite isolates whose PCT 50s was significantly longer than those in NW region. For the copy numbers of genes, parasite isolates containing multiple copies of pfmdr1 gene were found in both regions, but only single copy of pfpm2 gene was detected. Though the prevalence of parasite isolates with multiple copies of pfmdr1 gene in SW region was higher than that in NW region, no difference in PCT 50s were presented between isolates with single and multiple copies of pfmdr1 gene. The median He values of F446I group and Others (Non-F446I K13 mutation) group were 0.08 and 0.41 respectively. The mean He values of ML group (Menglian County in SW) and W (wild K13 genotype in NW) group were 0 and 0.69 respectively. The mean Fst values between ML and W groups were significantly higher than the other two K13 groups.

CONCLUSIONS

P. falciparum isolates in NW and SW regions had very different genetic characteristics. The F446I was hypothesized to have independently appeared and spread in NW region from 2012 and 2016. The high susceptibility of PIP had ensured the efficacy of DHA-PIP in vivo. Multiple copy numbers of pfmdr1 gene might be a potential cause of prolonged clearance time of ACTs drugs along China-Myanmar border.

TRIAL REGISTRATION

Trial registration: ISRCTN, ISRCTN 11775446. Registered 17 April 2020-Retrospectively registered, the registered name was Investigating resistance to DHA-PIP for the treatment of Plasmodium falciparum malaria and chloroquine for the treatment of Plasmodium vivax malaria in Yunnan, China. http://www.isrctn.com/ISRCTN11775446.

Sampling for malaria molecular surveillance

Strategic use of Plasmodium falciparum genetic variation has great potential to inform public health actions for malaria control and elimination. Malaria molecular surveillance (MMS) begins with a strategy to identify and collect parasite samples, guided by public-health priorities. In this review we discuss sampling design practices for MMS and point out epidemiological, biological, and statistical factors that need to be considered. We present examples for different use cases, including detecting emergence and spread of rare variants, establishing transmission sources and inferring changes in malaria transmission intensity. This review will potentially guide the collection of samples and data, serve as a starting point for further methodological innovation, and enhance utilization of MMS to support malaria elimination.

MalariaSED: a deep learning framework to decipher the regulatory contributions of noncoding variants in malaria parasites

Malaria remains one of the deadliest infectious diseases. Transcriptional regulation effects of noncoding variants in this unusual genome of malaria parasites remain elusive. We developed a sequence-based, ab initio deep learning framework, MalariaSED, for predicting chromatin profiles in malaria parasites. The MalariaSED performance was validated by published ChIP-qPCR and TF motifs results. Applying MalariaSED to ~ 1.3 million variants shows that geographically differentiated noncoding variants are associated with parasite invasion and drug resistance. Further analysis reveals chromatin accessibility changes at Plasmodium falciparum rings are partly associated with artemisinin resistance. MalariaSED illuminates the potential functional roles of noncoding variants in malaria parasites.

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.

Heat-shock responses: systemic and essential ways of malaria parasite survival

Fever is a part of the human innate immune response that contributes to limiting microbial growth and development in many infectious diseases. For the parasite Plasmodium falciparum, survival of febrile temperatures is crucial for its successful propagation in human populations as well as a fundamental aspect of malaria pathogenesis. This review discusses recent insights into the biological complexity of the malaria parasite's heat-shock response, which involves many cellular compartments and essential metabolic processes to alleviate oxidative stress and accumulation of damaged and unfolded proteins. We highlight the overlap between heat-shock and artemisinin resistance responses, while also explaining how the malaria parasite adapts its fever response to fight artemisinin treatment. Additionally, we discuss how this systemic and essential fight for survival can also contribute to parasite transmission to mosquitoes.

Naturally Acquired Kelch13 Mutations in Plasmodium falciparum Strains Modulate In Vitro Ring-Stage Artemisinin-Based Drug Tolerance and Parasite Survival in Response to Hyperoxia

The ring-stage survival assay was utilized to assess the impact of physiological hyperoxic stress on dihydroartemisinin (DHA) tolerance for a panel of Plasmodium falciparum strains with and without Kelch13 mutations. Strains without naturally acquired Kelch13 mutations or the postulated genetic background associated with delayed parasite clearance time demonstrated reduced proliferation under hyperoxic conditions in the subsequent proliferation cycle. Dihydroartemisinin tolerance in three isolates with naturally acquired Kelch13 mutations but not two genetically manipulated laboratory strains was modulated by in vitro hyperoxic stress exposure of early-ring-stage parasites in the cycle before drug exposure. Reduced parasite tolerance to additional derivatives, including artemisinin, artesunate, and OZ277, was observed within the second proliferation cycle. OZ439 and epoxomicin completely prevented parasite survival under both hyperoxia and normoxic in vitro culture conditions, highlighting the unique relationship between DHA tolerance and Kelch13 mutation-associated genetic background. IMPORTANCE Artemisinin-based combination therapy (ACT) for treating malaria is under intense scrutiny following treatment failures in the Greater Mekong subregion of Asia. This is further compounded by the potential for extensive loss of life if treatment failures extend to the African continent. Although Plasmodium falciparum has become resistant to all antimalarial drugs, artemisinin "resistance" does not present in the same way as resistance to other antimalarial drugs. Instead, a partial resistance or tolerance is demonstrated, associated with the parasite's genetic profile and linked to a molecular marker referred to as K13. It is suggested that parasites may have adapted to drug treatment, as well as the presence of underlying population health issues such as hemoglobinopathies, and/or environmental pressures, resulting in parasite tolerance to ACT. Understanding parasite evolution and control of artemisinin tolerance will provide innovative approaches to mitigate the development of artemisinin tolerance and thereby artemisinin-based drug treatment failure and loss of life globally to malaria infections.

Cyclical regression covariates remove the major confounding effect of cyclical developmental gene expression with strain-specific drug response in the malaria parasite Plasmodium falciparum

Background

The cyclical nature of gene expression in the intraerythrocytic development cycle (IDC) of the malaria parasite, Plasmodium falciparum, confounds the accurate detection of specific transcriptional differences, e.g. as provoked by the development of drug resistance. In lab-based studies, P. falciparum cultures are synchronized to remove this confounding factor, but the rapid detection of emerging resistance to artemisinin therapies requires rapid analysis of transcriptomes extracted directly from clinical samples. Here we propose the use of cyclical regression covariates (CRC) to eliminate the major confounding effect of developmentally driven transcriptional changes in clinical samples. We show that elimination of this confounding factor reduces both Type I and Type II errors and demonstrate the effectiveness of this approach using a published dataset of 1043 transcriptomes extracted directly from patient blood samples with different patient clearance times after treatment with artemisinin.

Results

We apply this method to two publicly available datasets and demonstrate its ability to reduce the confounding of differences in transcript levels due to misaligned intraerythrocytic development time. Adjusting the clinical 1043 transcriptomes dataset with CRC results in detection of fewer functional categories than previously reported from the same data set adjusted using other methods. We also detect mostly the same functional categories, but observe fewer genes within these categories. Finally, the CRC method identifies genes in a functional category that was absent from the results when the dataset was adjusted using other methods. Analysis of differential gene expression in the clinical data samples that vary broadly for developmental stage resulted in the detection of far fewer transcripts in fewer functional categories while, at the same time, identifying genes in two functional categories not present in the unadjusted data analysis. These differences are consistent with the expectation that CRC reduces both false positives and false negatives with the largest effect on datasets from samples with greater variance in developmental stage.

Conclusions

Cyclical regression covariates have immediate application to parasite transcriptome sequencing directly from clinical blood samples and to cost-constrained in vitro experiments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08281-y.

Machine learning for artemisinin resistance in malaria treatment across in vivo-in vitro platforms

Drug resistance has been rapidly evolving with regard to the first-line malaria treatment, artemisinin-based combination therapies. It has been an open question whether predictive models for this drug resistance status can be generalized across in vivo-in vitro transcriptomic measurements. In this study, we present a model that predicts artemisinin treatment resistance developed with transcriptomic information of Plasmodium falciparum. We demonstrated the robustness of this model across in vivo clearance rate and in vitro IC50 measurement and based on different microarray and data processing modalities. The validity of the algorithm is further supported by its first placement in the DREAM Malaria challenge. We identified transcription biomarkers to artemisinin treatment resistance that can predict artemisinin resistance and are conserved in their expression modules. This is a critical step in the research of malaria treatment, as it demonstrated the potential of a platform-robust, personalized model for artemisinin resistance using molecular biomarkers.

•

Artemisinin resistance can be predicted from transcriptomes by machine learning

•

Our model can be transferred between in vivo and in vitro and different platforms

•

We identified top transcription biomarkers of artemisinin resistance

The parasitophorous vacuole nutrient channel is critical for drug access in malaria parasites and modulates the artemisinin resistance fitness cost

Intraerythrocytic malaria parasites proliferate bounded by a parasitophorous vacuolar membrane (PVM). The PVM contains nutrient permeable channels (NPCs) conductive to small molecules, but their relevance for parasite growth for individual metabolites is largely untested. Here we show that growth-relevant levels of major carbon and energy sources pass through the NPCs. Moreover, we find that NPCs are a gate for several antimalarial drugs, highlighting their permeability properties as a critical factor for drug design. Looking into NPC-dependent amino acid transport, we find that amino acid shortage is a reason for the fitness cost in artemisinin-resistant (ART^R) parasites and provide evidence that NPC upregulation to increase amino acids acquisition is a mechanism of ART^R parasites in vitro and in human infections to compensate this fitness cost. Hence, the NPCs are important for nutrient and drug access and reveal amino acid deprivation as a critical constraint in ART^R parasites.

Artemisinin and multidrug-resistant Plasmodium falciparum - a threat for malaria control and elimination

PURPOSE OF REVIEW

Artemisinin-based combination therapies (ACTs) are globally the first-line treatment for uncomplicated falciparum malaria and new compounds will not be available within the next few years. Artemisinin-resistant Plasmodium falciparum emerged over a decade ago in the Greater Mekong Subregion (GMS) and, compounded by ACT partner drug resistance, has caused significant ACT treatment failure. This review provides an update on the epidemiology, and mechanisms of artemisinin resistance and approaches to counter multidrug-resistant falciparum malaria.

RECENT FINDINGS

An aggressive malaria elimination programme in the GMS has helped prevent the spread of drug resistance to neighbouring countries. However, parasites carrying artemisinin resistance-associated mutations in the P. falciparum Kelch13 gene (pfk13) have now emerged independently in multiple locations elsewhere in Asia, Africa and South America. Notably, artemisinin-resistant infections with parasites carrying the pfk13 R561H mutation have emerged and spread in Rwanda.

SUMMARY

Enhancing the geographic coverage of surveillance for resistance will be key to ensure prompt detection of emerging resistance in order to implement effective countermeasures without delay. Treatment strategies designed to prevent the emergence and spread of multidrug resistance must be considered, including deployment of triple drug combination therapies and multiple first-line therapies.

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

Resistance to antimalarial drugs, and in particular to the artemisinin derivatives and their partner drugs, threatens recent progress toward regional malaria elimination and eventual global malaria eradication. Population-level studies utilizing whole-genome sequencing approaches have facilitated the identification of regions of the parasite genome associated with both clinical and in vitro drug-resistance phenotypes. However, the biological relevance of genes identified in these analyses and the establishment of a causal relationship between genotype and phenotype requires functional characterization. Here we examined data from population genomic and transcriptomic studies in the context of data generated from recent functional studies, using a new population genetic approach designed to identify potential favored mutations within the region of a selective sweep (iSAFE). We identified several genes functioning in pathways now known to be associated with artemisinin resistance that were supported in early population genomic studies, as well as potential new drug targets/pathways for further validation and consideration for treatment of artemisinin-resistant Plasmodium falciparum. In addition, we establish the utility of iSAFE in identifying positively-selected mutations in population genomic studies, potentially accelerating the time to functional validation of candidate genes.

Molecular Mechanisms of Drug Resistance in Plasmodium falciparum Malaria

Understanding and controlling the spread of antimalarial resistance, particularly to artemisinin and its partner drugs, is a top priority. Plasmodium falciparum parasites resistant to chloroquine, amodiaquine, or piperaquine harbor mutations in the P. falciparum chloroquine resistance transporter (PfCRT), a transporter resident on the digestive vacuole membrane that in its variant forms can transport these weak-base 4-aminoquinoline drugs out of this acidic organelle, thus preventing these drugs from binding heme and inhibiting its detoxification. The structure of PfCRT, solved by cryogenic electron microscopy, shows mutations surrounding an electronegative central drug-binding cavity where they presumably interact with drugs and natural substrates to control transport. P. falciparum susceptibility to heme-binding antimalarials is also modulated by overexpression or mutations in the digestive vacuole membrane-bound ABC transporter PfMDR1 (P. falciparum multidrug resistance 1 transporter). Artemisinin resistance is primarily mediated by mutations in P. falciparum Kelch13 protein (K13), a protein involved in multiple intracellular processes including endocytosis of hemoglobin, which is required for parasite growth and artemisinin activation. Combating drug-resistant malaria urgently requires the development of new antimalarial drugs with novel modes of action.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

EPISPOT: An epigenome-driven approach for detecting and interpreting hotspots in molecular QTL studies

We present EPISPOT, a fully joint framework which exploits large panels of epigenetic annotations as variant-level information to enhance molecular quantitative trait locus (QTL) mapping. Thanks to a purpose-built Bayesian inferential algorithm, EPISPOT accommodates functional information for both cis and trans actions, including QTL hotspot effects. It effectively couples simultaneous QTL analysis of thousands of genetic variants and molecular traits with hypothesis-free selection of biologically interpretable annotations which directly contribute to the QTL effects. This unified, epigenome-aided learning boosts statistical power and sheds light on the regulatory basis of the uncovered hits; EPISPOT therefore marks an essential step toward improving the challenging detection and functional interpretation of trans-acting genetic variants and hotspots. We illustrate the advantages of EPISPOT in simulations emulating real-data conditions and in a monocyte expression QTL study, which confirms known hotspots and finds other signals, as well as plausible mechanisms of action. In particular, by highlighting the role of monocyte DNase-I sensitivity sites from >150 epigenetic annotations, we clarify the mediation effects and cell-type specificity of major hotspots close to the lysozyme gene. Our approach forgoes the daunting and underpowered task of one-annotation-at-a-time enrichment analyses for prioritizing cis and trans QTL hits and is tailored to any transcriptomic, proteomic, or metabolomic QTL problem. By enabling principled epigenome-driven QTL mapping transcriptome-wide, EPISPOT helps progress toward a better functional understanding of genetic regulation.

K13-Mediated Reduced Susceptibility to Artemisinin in Plasmodium falciparum Is Overlaid on a Trait of Enhanced DNA Damage Repair

Southeast Asia has been the hotbed for the development of drug-resistant malaria parasites, including those with resistance to artemisinin combination therapy. While mutations in the kelch propeller domain (K13 mutations) are associated with artemisinin resistance, a range of evidence suggests that other factors are critical for the establishment and subsequent transmission of resistance in the field. Here, we perform a quantitative analysis of DNA damage and repair in the malaria parasite Plasmodium falciparum and find a strong link between enhanced DNA damage repair and artemisinin resistance. This experimental observation is further supported when variations in seven known DNA repair genes are found in resistant parasites, with six of these mutations being associated with K13 mutations. Our data provide important insights on confounding factors that are important for the establishment and spread of artemisinin resistance and may explain why resistance has not yet arisen in Africa.

•

High-throughput MalariaCometChip to measure DNA damage level in P. falciparum

•

Subpopulation of Cambodian isolates possess enhanced DNA damage repair

•

Important link between enhanced DNA damage repair and artemisinin resistance

From Metabolite to Metabolome: Metabolomics Applications in Plasmodium Research

Advances in research over the past few decades have greatly improved metabolomics-based approaches in studying parasite biology and disease etiology. This improves the investigation of varied metabolic requirements during life stages or when following transmission to their hosts, and fulfills the demand for improved diagnostics and precise therapeutics. Therefore, this review highlights the progress of metabolomics in malaria research, including metabolic mapping of Plasmodium vertebrate life cycle stages to investigate antimalarials mode of actions and underlying complex host-parasite interactions. Also, we discuss current limitations as well as make several practical suggestions for methodological improvements which could drive metabolomics progress for malaria from a comprehensive perspective.

Molecular approaches to Malaria 2020

Twenty years ago the Molecular Approaches to Malaria conference was conceived as a forum to present the very latest advances in malaria research and to consolidate and forge new collaborative links between international researchers. The 6th MAM conference, held in February 2020 in Australia, provided 5 days of stimulating scientific exchange and highlighted the incredible malaria research conducted globally that is providing the critical knowledge and cutting-edge technological tools needed to control and ultimately eliminate malaria.

The human malaria parasite genome is configured into thousands of coexpressed linear regulatory units

The human malaria parasite Plasmodium falciparum thrives in radically different host environments in mosquitoes and humans, with only a limited set of transcription factors. The nature of regulatory elements or their target genes in the P. falciparum genome remains elusive. Here, we found that this eukaryotic parasite uses an efficient way to maximally use genetic and epigenetic regulation to form regulatory units (RUs) during blood infections. Genes located in the same RU tend to have the same pattern of expression over time and are associated with open chromatin along regulatory elements. To precisely define and quantify these RUs, a novel hidden Markov model was developed to capture the regulatory structure in a genome-wide fashion by integrating expression and epigenetic evidence. We successfully identified thousands of RUs and cross-validated with previous findings. We found more genes involved in red blood cell (RBC) invasion located in the same RU as the PfAP2-I (AP2-I) transcription factor, demonstrating that AP2-I is responsible for regulating RBC invasion. Our study has provided a regulatory mechanism for a compact eukaryotic genome and offers new insights into the in vivo transcriptional regulation of the P. falciparum intraerythrocytic stage.

A comprehensive RNA handling and transcriptomics guide for high-throughput processing of Plasmodium blood-stage samples

Background

Sequencing technology advancements opened new opportunities to use transcriptomics for studying malaria pathology and epidemiology. Even though in recent years the study of whole parasite transcriptome proved to be essential in understanding parasite biology there is no compiled up-to-date reference protocol for the efficient generation of transcriptome data from growing number of samples. Here, a comprehensive methodology on how to preserve, extract, amplify, and sequence full-length mRNA transcripts from Plasmodium-infected blood samples is presented that can be fully streamlined for high-throughput studies.

Results

The utility of various commercially available RNA-preserving reagents in a range of storage conditions was evaluated. Similarly, several RNA extraction protocols were compared and the one most suitable method for the extraction of high-quality total RNA from low-parasitaemia and low-volume blood samples was established. Furthermore, the criteria needed to evaluate the quality and integrity of Plasmodium RNA in the presence of human RNA was updated. Optimization of SMART-seq2 amplification method to better suit AT-rich Plasmodium falciparum RNA samples allowed us to generate high-quality transcriptomes from as little as 10 ng of total RNA and a lower parasitaemia limit of 0.05%. Finally, a modified method for depletion of unwanted human haemoglobin transcripts using in vitro CRISPR-Cas9 treatment was designed, thus improving parasite transcriptome coverage in low parasitaemia samples. To prove the functionality of the pipeline for both laboratory and field strains, the highest  2-hour resolution RNA-seq transcriptome for P. falciparum 3D7 intraerythrocytic life cycle available to  date was generated, and the entire protocol was applied to create the largest transcriptome data from Southeast Asian field isolates.

Conclusions

Overall, the presented methodology is an inclusive pipeline for generation of good quality transcriptomic data from a diverse range of Plasmodium-infected blood samples with varying parasitaemia and RNA inputs. The flexibility of this pipeline to be adapted to robotic handling will facilitate both small and large-scale future transcriptomic studies in the field of malaria.

Arbuscular Mycorrhizal Fungus Rhizophagus irregularis Influences Artemisia annua Plant Parameters and Artemisinin Content under Different Soil Types and Cultivation Methods

Artemisinin extracted from Artemisia annua has been used efficiently in malaria treatment since 2005. In this study, the variations in plant parameters (plant biomass, glandular trichome density, essential oil total chemical content, artemisinin production, and polyphenol oxidase (PPO) activity) were tested under different soil types (Luvisol, Gleysol, Anthrosol and sterile peat) and cultivation conditions (potted plants in semi-open field, and open field experiments) for plants inoculated with arbuscular mycorrhizal fungus (AMF) Rizophagus irregularis. Under semi-open field conditions, the AMF colonization of A. annua plant roots varied, and presented the highest percentage in Luvisol and sterile peat. The increase in the root colonization rate positively influenced some plant parameters (biomass, glandular trichome density, artemisinin concentration, essential oil quantity and composition), but no effects on PPO enzyme activity were detected. AMF fungus R. irregularis significantly increased the artemisinin content and essential oil yield of plants cultivated in Luvisol, Gleysol, Anthrosol and in peat. These soil types can offer appropriate conditions for A. annua cultivation and artemisinin production even on a smaller scale. Under open field conditions, low (about 5%) AMF colonization was observed. No differences in artemisin contents were detected, but essential oil yield significantly increased compared to control plants. AMF treatment increased beta-farnesene and germacrene D concentrations in Artemisia plants in the open field experiment.

Transcriptional variation in malaria parasites: why and how

Transcriptional differences enable the generation of alternative phenotypes from the same genome. In malaria parasites, transcriptional plasticity plays a major role in the process of adaptation to fluctuations in the environment. Multiple studies with culture-adapted parasites and field isolates are starting to unravel the different transcriptional alternatives available to Plasmodium falciparum and the underlying molecular mechanisms. Here we discuss how epigenetic variation, directed transcriptional responses and also genetic changes that affect transcript levels can all contribute to transcriptional variation and, ultimately, parasite survival. Some transcriptional changes are driven by stochastic events. These changes can occur spontaneously, resulting in heterogeneity within parasite populations that provides the grounds for adaptation by dynamic natural selection. However, transcriptional changes can also occur in response to external cues. A better understanding of the mechanisms that the parasite has evolved to alter its transcriptome may ultimately contribute to the design of strategies to combat malaria to which the parasite cannot adapt.

Role of Plasmodium falciparum Kelch 13 Protein Mutations in P. falciparum Populations from Northeastern Myanmar in Mediating Artemisinin Resistance

Artemisinin resistance has emerged in Southeast Asia, endangering the substantial progress in malaria elimination worldwide. It is associated with mutations in the PfK13 protein, but how PfK13 mediates artemisinin resistance is not completely understood. Here we used a new antibody against PfK13 to show that the PfK13 protein is expressed in all stages of the asexual intraerythrocytic cycle as well as in gametocytes and is partially localized in the endoplasmic reticulum. By introducing four PfK13 mutations into the 3D7 strain and reverting these mutations in field parasite isolates, we determined the impacts of these mutations identified in the parasite populations from northern Myanmar on the ring stage using the in vitro ring survival assay. The introduction of the N458Y mutation into the 3D7 background significantly increased the survival rates of the ring-stage parasites but at the cost of the reduced fitness of the parasites. Introduction of the F446I mutation, the most prevalent PfK13 mutation in northern Myanmar, did not result in a significant increase in ring-stage survival after exposure to dihydroartemisinin (DHA), but these parasites showed extended ring-stage development. Further, parasites with the F446I mutation showed only a marginal loss of fitness, partially explaining its high frequency in northern Myanmar. Conversely, reverting all these mutations, except for the C469Y mutation, back to their respective wild types reduced the ring-stage survival of these isolates in response to in vitro DHA treatment.

Mutations in the Plasmodium falciparum Kelch 13 (PfK13) protein are associated with artemisinin resistance. PfK13 is essential for asexual erythrocytic development, but its function is not known. We tagged the PfK13 protein with green fluorescent protein in P. falciparum to study its expression and localization in asexual and sexual stages. We used a new antibody against PfK13 to show that the PfK13 protein is expressed ubiquitously in both asexual erythrocytic stages and gametocytes and is localized in punctate structures, partially overlapping an endoplasmic reticulum marker. We introduced into the 3D7 strain four PfK13 mutations (F446I, N458Y, C469Y, and F495L) identified in parasites from the China-Myanmar border area and characterized the in vitro artemisinin response phenotypes of the mutants. We found that all the parasites with the introduced PfK13 mutations showed higher survival rates in the ring-stage survival assay (RSA) than the wild-type (WT) control, but only parasites with N458Y displayed a significantly higher RSA value (26.3%) than the WT control. After these PfK13 mutations were reverted back to the WT in field parasite isolates, all revertant parasites except those with the C469Y mutation showed significantly lower RSA values than their respective parental isolates. Although the 3D7 parasites with introduced F446I, the predominant PfK13 mutation in northern Myanmar, did not show significantly higher RSA values than the WT, they had prolonged ring-stage development and showed very little fitness cost in in vitro culture competition assays. In comparison, parasites with the N458Y mutations also had a prolonged ring stage and showed upregulated resistance pathways in response to artemisinin, but this mutation produced a significant fitness cost, potentially leading to their lower prevalence in the Greater Mekong subregion.

Transcriptome profiling reveals functional variation in Plasmodium falciparum parasites from controlled human malaria infection studies

BACKGROUND

The transcriptome of Plasmodium falciparum clinical isolates varies according to strain, mosquito bites, disease severity and clinical history. Therefore, it remains a challenge to directly interpret the parasite's transcriptomic information into a more general biological signature in a natural human malaria infection. These confounding variations can be potentially overcome with parasites derived from controlled-human malaria infection (CHMI) studies.

METHODS

We performed CHMI studies in healthy and immunologically naïve volunteers receiving the same P. falciparum strain ((Sanaria® PfSPZ Challenge (NF54)), but with different sporozoite dosage and route of infection. Parasites isolated from these volunteers at the day of patency were subjected to in vitro culture for several generations and synchronized ring-stage parasites were subjected to transcriptome profiling.

FINDINGS

We observed clear deviations between CHMI-derived parasites from volunteer groups receiving different PfSPZ dose and route. CHMI-derived parasites and the pre-mosquito strain used for PfSPZ generation showed significant transcriptional variability for gene clusters associated with malaria pathogenesis, immune evasion and transmission. These transcriptional variation signature clusters were also observed in the transcriptome of P. falciparum isolates from acute clinical infections.

INTERPRETATION

Our work identifies a previously unrecognized transcriptional pattern in malaria infections in a non-immune background. Significant transcriptome heterogeneity exits between parasites derived from human infections and the pre-mosquito strain, implying that the malaria parasites undergo a change in functional state to adapt to its host environment. Our work also highlights the potential use of transcriptomics data from CHMI study advance our understanding of malaria parasite adaptation and transmission in humans. FUND: This work is supported by German Israeli Foundation, German ministry for education and research, MOE Tier 1 from the Singapore Ministry of Education Academic Research Fund, Singapore Ministry of Health's National Medical Research Council, National Institute of Allergy and Infectious Diseases, National Institutes of Health, USA and the German Centre for Infection Research (Deutsches Zentrum für Infektionsforschung-DZIF).