Plasmodium parasites mount an arrest response to dihydroartemisinin, as revealed by whole transcriptome shotgun sequencing (RNA-seq) and microarray study

Understanding the Antiplasmodial Action of Resistance-Refractory Xanthoquinodin A1

Our previous work identified a series of 12 xanthoquinodin analogues and 2 emodin-dianthrones with broad-spectrum activities against Trichomonas vaginalis, Mycoplasma genitalium, Cryptosporidium parvum, and Plasmodium falciparum. Analyses conducted in this study revealed that the most active analogue, xanthoquinodin A1, also inhibits Toxoplasma gondii tachyzoites and the liver stage of Plasmodium berghei, with no cross-resistance to the known antimalarial targets PfACS, PfCARL, PfPI4K, or DHODH. In Plasmodium, inhibition occurs prior to multinucleation and induces parasite death following 12 h of compound exposure. This moderately fast activity has impeded resistance line generation, with xanthoquinodin A1 demonstrating an irresistible phenotype in both T. gondii and P. falciparum.

Molecular characterization and functional analysis of Eimeria tenella ankyrin repeat-containing protein

Chicken coccidiosis causes disastrous losses to the poultry industry all over the world. Eimeria tenella is the most prevalent of these disease-causing species. Our former RNA-seq indicated that E. tenella ankyrin repeat-containing protein (EtANK) was expressed differently between drug-sensitive (DS) and drug-resistant strains. In this study, we cloned EtANK and analyzed its translational and transcriptional levels using quantitative real-time PCR (qPCR) and western blotting. The data showed that EtANK was significantly upregulated in diclazuril-resistant (DZR) strain and maduramicin-resistant (MRR) strain compared with the drug-sensitive (DS) strain. In addition, the transcription levels in the DZR strains isolated from the field were higher than in the DS strain. The translation levels of EtANK were higher in unsporulated oocysts (UO) than in sporozoites (SZ), sporulated oocysts (SO), or second-generation merozoites (SM), and the protein levels in SM were significantly higher than in UO, SO, and SZ. The results of the indirect immunofluorescence localization showed that the protein was distributed mainly at the anterior region of SZ and on the surface and in the cytoplasm of SM. The fluorescence intensity increased further with its development in vitro. An anti-rEtANK polyclonal antibody inhibited the invasive ability of E. tenella in DF-1 cells. These results showed that EtANK may be related to host cell invasion, required for the parasite's growth in the host, and may be involved in the development of E. tenella resistance to some drugs.

Molecular characterization of glyceraldehyde-3-phosphate dehydrogenase from Eimeria tenella

Chicken coccidiosis is an extremely common and lethally epidemic disease caused by Eimeria spp. The control measures of coccidiosis depend mainly on drugs. However, the ensuing drug resistance problem has brought considerable economic loss to the poultry industry. In our previous study, comparative transcriptome analyses of a drug-sensitive (DS) strain and two drug-resistant strains (diclazuril-resistant (DZR) and maduramicin-resistant (MRR) strains) of Eimeria tenella were carried out by transcriptome sequencing. The expression of glyceraldehyde-3-phosphate dehydrogenase of E. tenella (EtGAPDH) was upregulated in the two resistant strains. In this study, we cloned and characterized EtGAPDH. Indirect immunofluorescence localization was used to observe the distribution of EtGAPDH in E. tenella. The results showed that the protein was distributed mainly on the surface of sporozoites and merozoites, and in the cytoplasm of merozoites. qPCR was performed to detect the transcription level of EtGAPDH in the different developmental stages of the E. tenella DS strain. The transcription level of EtGAPDH was significantly higher in second-generation merozoites than in the other three stages. The transcription level of EtGAPDH in the different drug-resistant strains and DS strain of E. tenella was also analyzed by qPCR. The results showed that the transcription level was significantly higher in the two drug-resistant strains (MRR and DZR) than in the DS strain. As the concentration of diclazuril and maduramicin increased, the transcription levels also increased. Western blot results showed that EtGAPDH protein was upregulated in the DZR and MRR strains. Enzyme activity showed that the enzyme activity of EtGAPDH was higher in the two resistant strains than in the DS strain. These results showed that EtGAPDH possess several roles that separate and distinct from its glycolytic function and maybe involved in the development of E. tenella resistance to anticoccidial drugs.

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.

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.

Molecular characterization and analysis of the ATPase ASNA1 homolog gene of Eimeria tenella in a drug sensitive strain and drug resistant strains

The widespread use of drugs has exacerbated the resistance of Eimeria tenalla to anti-coccidial drugs. Using RNA-seq, we previously found the ATPase ASNA1 homolog of E. tenella (EtASNA1) was differentially expressed in resistant strains and drug sensitive (DS) strain. In our study, we used western blotting and quantitative real-time PCR (qRT-PCR) to analyze the translational and transcriptional levels of EtASNA1 in a diclazuril-resistant (DZR) strain, maduramicin-resistant (MRR) strain, salinomycin-resistant (SMR) strain, and DS strain and found EtASNA1 was highly expressed in three drug-resistant strains. The qRT-PCR and western blotting results also showed that the expression levels of EtASNA1 increased with increasing drug concentration, and the transcription levels of the DZR strains isolated from the field were higher than those of the DS strain. In addition, we used in vivo and in vitro tests to analyze the changes of EtASNA1 expression after DZR, MRR, and DS strain infections in chickens, and in vitro inoculation of DF-1 cells in the presence of drugs. The addition of drugs caused expression to be upregulated. The results of qRT-PCR and western blotting also showed that the expression levels of EtASNA1 in second-generation merozoites (SM) and unsporulated oocysts (UO) were significantly higher than those in the other two developmental stages. The immunofluorescence localization of EtASNA1 indicated that the protein was distributed throughout the sporozoites (SZ) and SM, except for the refractile bodies of SZ. In vitro inhibition experiments showed that anti-EtASNA1 antibody incubation significantly inhibited SZ invasion of DF-1 cells. The above results showed that EtASNA1 may be related to host cell invasion of E. tenella and may be involved in the development of E. tenella resistance to some drugs.

Molecular characterization and functional analysis of Eimeria tenella citrate synthase

Chicken coccidiosis, caused by an obligate intracellular protozoan parasite of the genus Eimeria, is a major parasitic disease in the intensively reared poultry industry. Due to the widespread use of anticoccidial drugs, resistance has become an inevitable problem. In our previous study, Eimeria tenella citrate synthase (EtCS) was found to be up-expressed in two drug-resistant strains (diclazuril-resistant and maduramycin-resistant strains) compared to drug-sensitive strain by RNA sequence. In this study, we cloned and expressed EtCS and obtain its polyclonal antibodies. Quantitative real-time polymerase chain (qPCR) reactions and Western blots were used to analyze the transcription and translation levels of EtCS in sensitive and three drug-resistant strains. Compared with the sensitive strain, the transcription of EtCS was both significantly upregulated in diclazuril-resistant and maduramycin-resistant strains, but was not significantly different in salinomycin-resistant strain. No significant difference was seen in translation level in the three drug-resistant strains. Indirect immunofluorescence indicated that EtCS was mainly located in the cytoplasm of sporozoites except for posterior refractile bodies and in the cytoplasm and surface of merozoites. Anti-rEtCS antibody has inhibitory effects on E. tenella sporozoite invasion of DF-1 cells and the inhibition rate is more than 83%. Binding of the protein to chicken macrophage (HD11) cells was confirmed by immunofluorescence assays. When macrophages were treated with rEtCS, secretion of nitric oxide and cell proliferation of the macrophages were substantially reduced. These results showed that EtCS may be related to host cell invasion of E. tenella and involve in the development of E.tenella resistance to some drugs.

Artemisinin-resistant K13 mutations rewire Plasmodium falciparum's intra-erythrocytic metabolic program to enhance survival

The emergence and spread of artemisinin resistance, driven by mutations in Plasmodium falciparum K13, has compromised antimalarial efficacy and threatens the global malaria elimination campaign. By applying systems-based quantitative transcriptomics, proteomics, and metabolomics to a panel of isogenic K13 mutant or wild-type P. falciparum lines, we provide evidence that K13 mutations alter multiple aspects of the parasite's intra-erythrocytic developmental program. These changes impact cell-cycle periodicity, the unfolded protein response, protein degradation, vesicular trafficking, and mitochondrial metabolism. K13-mediated artemisinin resistance in the Cambodian Cam3.II line was reversed by atovaquone, a mitochondrial electron transport chain inhibitor. These results suggest that mitochondrial processes including damage sensing and anti-oxidant properties might augment the ability of mutant K13 to protect P. falciparum against artemisinin action by helping these parasites undergo temporary quiescence and accelerated growth recovery post drug elimination.

From Genes to Transcripts, a Tightly Regulated Journey in Plasmodium

Over the past decade, we have witnessed significant progresses in understanding gene regulation in Apicomplexa including the human malaria parasite, Plasmodium falciparum. This parasite possesses the ability to convert in multiple stages in various hosts, cell types, and environments. Recent findings indicate that P. falciparum is talented at using efficient and complementary molecular mechanisms to ensure a tight control of gene expression at each stage of its life cycle. Here, we review the current understanding on the contribution of the epigenome, atypical transcription factors, and chromatin organization to regulate stage conversion in P. falciparum. The adjustment of these regulatory mechanisms occurring during the progression of the life cycle will be extensively discussed.

Describing the current status of Plasmodium falciparum population structure and drug resistance within mainland Tanzania using molecular inversion probes

High-throughput Plasmodium genomic data is increasingly useful in assessing prevalence of clinically important mutations and malaria transmission patterns. Understanding parasite diversity is important for identification of specific human or parasite populations that can be targeted by control programmes, and to monitor the spread of mutations associated with drug resistance. An up-to-date understanding of regional parasite population dynamics is also critical to monitor the impact of control efforts. However, this data is largely absent from high-burden nations in Africa, and to date, no such analysis has been conducted for malaria parasites in Tanzania countrywide. To this end, over 1,000 P. falciparum clinical isolates were collected in 2017 from 13 sites in seven administrative regions across Tanzania, and parasites were genotyped at 1,800 variable positions genome-wide using molecular inversion probes. Population structure was detectable among Tanzanian P. falciparum parasites, approximately separating parasites from the northern and southern districts and identifying genetically admixed populations in the north. Isolates from nearby districts were more likely to be genetically related compared to parasites sampled from more distant districts. Known drug resistance mutations were seen at increased frequency in northern districts (including two infections carrying pfk13-R561H), and additional variants with undetermined significance for antimalarial resistance also varied by geography. Malaria Indicator Survey (2017) data corresponded with genetic findings, including average region-level complexity-of-infection and malaria prevalence estimates. The parasite populations identified here provide important information on extant spatial patterns of genetic diversity of Tanzanian parasites, to which future surveys of genetic relatedness can be compared.

Identification of compounds active against quiescent artemisinin-resistant Plasmodium falciparum parasites via the quiescent-stage survival assay (QSA)

BACKGROUND

Quiescence is an unconventional mechanism of Plasmodium survival, mediating artemisinin resistance. This phenomenon increases the risk of clinical failures following artemisinin-based combination therapies (ACTs) by slowing parasite clearance and allowing the selection of parasites resistant to partner drugs.

OBJECTIVES

To thwart this multiresistance, the quiescent state of artemisinin-resistant parasites must be taken into consideration from the very early stages of the drug discovery process.

METHODS

We designed a novel phenotypic assay we have named the quiescent-stage survival assay (QSA) to assess the antiplasmodial activity of drugs on quiescent parasites. This assay was first validated on quiescent forms from different artemisinin-resistant parasite lines (laboratory strain and field isolates), using two reference drugs with different mechanisms of action: chloroquine and atovaquone. Furthermore, the efficacies of different partner drugs of artemisinins used in ACTs were investigated against both laboratory strains and field isolates from Cambodia.

RESULTS

Our results highlight that because of the mechanism of quiescence and the respective pharmacological targets of drugs, drug efficacies on artemisinin-resistant parasites may be different between quiescent parasites and their proliferating forms.

CONCLUSIONS

These data confirm the high relevance of adding the chemosensitivity evaluation of quiescent parasites by the specific in vitro QSA to the antiplasmodial drug development process in the current worrisome context of artemisinin resistance.

A Biosynthetic Platform for Antimalarial Drug Discovery

Advances in synthetic biology have enabled the production of a variety of compounds using bacteria as a vehicle for complex compound biosynthesis. Violacein, a naturally occurring indole pigment with antibiotic properties, can be biosynthetically engineered in Escherichia coli expressing its nonnative synthesis pathway. To explore whether this synthetic biosynthesis platform could be used for drug discovery, here we have screened bacterially derived violacein against the main causative agent of human malaria, Plasmodium falciparum We show the antiparasitic activity of bacterially derived violacein against the P. falciparum 3D7 laboratory reference strain as well as drug-sensitive and -resistant patient isolates, confirming the potential utility of this drug as an antimalarial agent. We then screen a biosynthetic series of violacein derivatives against P. falciparum growth. The varied activity of each derivative against asexual parasite growth points to the need to further develop violacein as an antimalarial. Towards defining its mode of action, we show that biosynthetic violacein affects the parasite actin cytoskeleton, resulting in an accumulation of actin signal that is independent of actin polymerization. This activity points to a target that modulates actin behavior in the cell either in terms of its regulation or its folding. More broadly, our data show that bacterial synthetic biosynthesis could become a suitable platform for antimalarial drug discovery, with potential applications in future high-throughput drug screening with otherwise chemically intractable natural products.

Insights into the intracellular localization, protein associations and artemisinin resistance properties of Plasmodium falciparum K13

The emergence of artemisinin (ART) resistance in Plasmodium falciparum intra-erythrocytic parasites has led to increasing treatment failure rates with first-line ART-based combination therapies in Southeast Asia. Decreased parasite susceptibility is caused by K13 mutations, which are associated clinically with delayed parasite clearance in patients and in vitro with an enhanced ability of ring-stage parasites to survive brief exposure to the active ART metabolite dihydroartemisinin. Herein, we describe a panel of K13-specific monoclonal antibodies and gene-edited parasite lines co-expressing epitope-tagged versions of K13 in trans. By applying an analytical quantitative imaging pipeline, we localize K13 to the parasite endoplasmic reticulum, Rab-positive vesicles, and sites adjacent to cytostomes. These latter structures form at the parasite plasma membrane and traffic hemoglobin to the digestive vacuole wherein artemisinin-activating heme moieties are released. We also provide evidence of K13 partially localizing near the parasite mitochondria upon treatment with dihydroartemisinin. Immunoprecipitation data generated with K13-specific monoclonal antibodies identify multiple putative K13-associated proteins, including endoplasmic reticulum-resident molecules, mitochondrial proteins, and Rab GTPases, in both K13 mutant and wild-type isogenic lines. We also find that mutant K13-mediated resistance is reversed upon co-expression of wild-type or mutant K13. These data help define the biological properties of K13 and its role in mediating P. falciparum resistance to ART treatment.

Validation of Plasmodium falciparum deoxyhypusine synthase as an antimalarial target

Background

Hypusination is an essential post-translational modification in eukaryotes. The two enzymes required for this modification, namely deoxyhypusine synthase (DHS) and deoxyhypusine hydrolase are also conserved. Plasmodium falciparum human malaria parasites possess genes for both hypusination enzymes, which are hypothesized to be targets of antimalarial drugs.

Methods

Transgenic P. falciparum parasites with modification of the PF3D7_1412600 gene encoding PfDHS enzyme were created by insertion of the glmS riboswitch or the M9 inactive variant. The PfDHS protein was studied in transgenic parasites by confocal microscopy and Western immunoblotting. The biochemical function of PfDHS enzyme in parasites was assessed by hypusination and nascent protein synthesis assays. Gene essentiality was assessed by competitive growth assays and chemogenomic profiling.

Results

Clonal transgenic parasites with integration of glmS riboswitch downstream of the PfDHS gene were established. PfDHS protein was present in the cytoplasm of transgenic parasites in asexual stages. The PfDHS protein could be attenuated fivefold in transgenic parasites with an active riboswitch, whereas PfDHS protein expression was unaffected in control transgenic parasites with insertion of the riboswitch-inactive sequence. Attenuation of PfDHS expression for 72 h led to a significant reduction of hypusinated protein; however, global protein synthesis was unaffected. Parasites with attenuated PfDHS expression showed a significant growth defect, although their decline was not as rapid as parasites with attenuated dihydrofolate reductase-thymidylate synthase (PfDHFR-TS) expression. PfDHS-attenuated parasites showed increased sensitivity to N ¹-guanyl-1,7-diaminoheptane, a structural analog of spermidine, and a known inhibitor of DHS enzymes.

Discussion

Loss of PfDHS function leads to reduced hypusination, which may be important for synthesis of some essential proteins. The growth defect in parasites with attenuated Pf DHS expression suggests that this gene is essential. However, the slower decline of PfDHS mutants compared with PfDHFR-TS mutants in competitive growth assays suggests that PfDHS is less vulnerable as an antimalarial target. Nevertheless, the data validate PfDHS as an antimalarial target which can be inhibited by spermidine-like compounds.

A fast and cost-effective microsampling protocol incorporating reduced animal usage for time-series transcriptomics in rodent malaria parasites

BACKGROUND

The transcriptional regulation that occurs in malaria parasites during the erythrocytic stages of infection can be studied in vivo with rodent malaria parasites propagated in mice. Time-series transcriptome profiling commonly involves the euthanasia of groups of mice at specific time points followed by the extraction of parasite RNA from whole blood samples. Current methodologies for parasite RNA extraction involve several steps and when multiple time points are profiled, these protocols are laborious, time-consuming, and require the euthanization of large cohorts of mice.

RESULTS

A simplified protocol has been designed for parasite RNA extraction from blood volumes as low as 20 μL (microsamples), serially bled from mice via tail snips and directly lysed with TRIzol reagent. Gene expression data derived from microsampling using RNA-seq were closely matched to those derived from larger volumes of leucocyte-depleted and saponin-treated blood obtained from euthanized mice with high reproducibility between biological replicates. Transcriptome profiling of microsamples taken at different time points during the intra-erythrocytic developmental cycle of the rodent malaria parasite Plasmodium vinckei revealed the transcriptional cascade commonly observed in malaria parasites.

CONCLUSIONS

Microsampling is a quick, robust and cost-efficient approach to sample collection for in vivo time-series transcriptomic studies in rodent malaria parasites.

Molecular and cellular characterization of apoptosis in flat oyster a key mechanisms at the heart of host-parasite interactions

Bonamia ostreae has been associated with the decline of flat oyster Ostrea edulis populations in some European countries. This obligatory intracellular parasite persists and multiplies into hemocytes. Previous in vitro experiments showed that apoptosis is activated in hemocytes between 1 h and 4 h of contact with the parasite. The flat oyster uses the apoptosis pathway to defend against B. ostreae. However, the parasite might be also able to modulate this response in order to survive in its host. In order to investigate this hypothesis the apoptotic response of the host was evaluated using flow cytometry, transmission electron microscopy and by measuring the response of genes involved in the apoptotic pathway after 4 h. In parallel, the parasite response was investigated by measuring the expression of B. ostreae genes involved in different biological functions including cell cycle and cell death. Obtained results allow describing molecular apoptotic pathways in O. edulis and confirm that apoptosis is early activated in hemocytes after a contact with B. ostreae. Interestingly, at cellular and molecular levels this process appeared downregulated after 44 h of contact. Concurrently, parasite gene expression appeared reduced suggesting that the parasite could inhibit its own metabolism to escape the immune response.

Antimalarial Drug Resistance: A Threat to Malaria Elimination

Increasing antimalarial drug resistance once again threatens effective antimalarial drug treatment, malaria control, and elimination. Artemisinin combination therapies (ACTs) are first-line treatment for uncomplicated falciparum malaria in all endemic countries, yet partial resistance to artemisinins has emerged in the Greater Mekong Subregion. Concomitant emergence of partner drug resistance is now causing high ACT treatment failure rates in several areas. Genetic markers for artemisinin resistance and several of the partner drugs have been established, greatly facilitating surveillance. Single point mutations in the gene coding for the Kelch propeller domain of the K13 protein strongly correlate with artemisinin resistance. Novel regimens and strategies using existing antimalarial drugs will be needed until novel compounds can be deployed. Elimination of artemisinin resistance will imply elimination of all falciparum malaria from the same areas. In vivax malaria, chloroquine resistance is an increasing problem.

Knockout of the peroxiredoxin 5 homologue PFAOP does not affect the artemisinin susceptibility of Plasmodium falciparum

Artemisinins are the current mainstay of malaria chemotherapy. Their exact mode of action is an ongoing matter of debate, and several factors have recently been reported to affect an early stage of artemisinin resistance of the most important human malaria parasite Plasmodium falciparum. Here, we identified a locus on chromosome 7 that affects the artemisinin susceptibility of P. falciparum in a quantitative trait locus analysis of a genetic cross between strains 7G8 and GB4. This locus includes the peroxiredoxin gene PFAOP. However, steady-state kinetic data with recombinant PfAOP do not support a direct interaction between this peroxidase and the endoperoxide artemisinin. Furthermore, neither the overexpression nor the deletion of the encoding gene affected the IC₅₀ values for artemisinin or the oxidants diamide and tert-butyl hydroperoxide. Thus, PfAOP is dispensable for blood stage parasite survival, and the correlation between the artemisinin susceptibility and chromosome 7 is probably based on another gene within the identified locus.

Correlation between Cyclin Dependent Kinases and Artemisinin-Induced Dormancy in Plasmodium falciparum In Vitro

Background

Artemisinin-induced dormancy provides a plausible explanation for recrudescence following artemisinin monotherapy. This phenomenon shares similarities with cell cycle arrest where cyclin dependent kinases (CDKs) and cyclins play an important role.

Methods

Transcription profiles of Plasmodium falciparum CDKs and cyclins before and after dihydroartemisinin (DHA) treatment in three parasite lines, and the effect of CDK inhibitors on parasite recovery from DHA-induced dormancy were investigated.

Results

After DHA treatment, parasites enter a dormancy phase followed by a recovery phase. During the dormancy phase parasites up-regulate pfcrk1, pfcrk4, pfcyc2 and pfcyc4, and down-regulate pfmrk, pfpk5, pfpk6, pfcrk3, pfcyc1 and pfcyc3. When entering the recovery phase parasites immediately up-regulate all CDK and cyclin genes. Three CDK inhibitors, olomoucine, WR636638 and roscovitine, produced distinct effects on different phases of DHA-induced dormancy, blocking parasites recovery.

Conclusions

The up-regulation of PfCRK1 and PfCRK4, and down regulation of other CDKs and cyclins correlate with parasite survival in the dormant state. Changes in CDK expression are likely to negatively regulate parasite progression from G₁ to S phase. These findings provide new insights into the mechanism of artemisinin-induced dormancy and cell cycle regulation of P. falciparum, opening new opportunities for preventing recrudescence following artemisinin treatment.

Plasmodium parasites mount an arrest response to dihydroartemisinin, as revealed by whole transcriptome shotgun sequencing (RNA-seq) and microarray study

Background

Control of malaria is threatened by emerging parasite resistance to artemisinin and derivative drug (ART) therapies. The molecular detail of how Plasmodium malaria parasites respond to ART and how this could contribute to resistance are not well understood. To address this question, we performed a transcriptomic study of dihydroartemisinin (DHA) response in P. falciparum K1 strain and in P. berghei ANKA strain using microarray and RNA-seq technology.

Results

Microarray data from DHA-treated P. falciparum trophozoite stage parasites revealed a response pattern that is overall less trophozoite-like and more like the other stages of asexual development. A meta-analysis of these data with previously published data from other ART treatments revealed a set of common differentially expressed genes. Notably, ribosomal protein genes are down-regulated in response to ART. A similar pattern of trophozoite transcriptomic change was observed from RNA-seq data. RNA-seq data from DHA-treated P. falciparum rings reveal a more muted response, although there is considerable overlap of differentially expressed genes with DHA-treated trophozoites. No genes are differentially expressed in DHA-treated P. falciparum schizonts. The transcriptional response of P. berghei to DHA treatment in vivo in infected mice is similar to the P. falciparum in vitro culture ring and trophozoite responses, in which ribosomal protein genes are notably down-regulated.

Conclusions

Ring and trophozoite stage Plasmodium respond to ART by arresting metabolic processes such as protein synthesis and glycolysis. This response can be protective in rings, as shown by the phenomenon of dormancy. In contrast, this response is not as protective in trophozoites owing to their commitment to a highly active and vulnerable metabolic state. The lower metabolic demands of schizonts could explain why they are less sensitive and unresponsive to ART. The ART response pattern is revealed clearly from RNA-seq data, suggesting that this technology is of great utility for studying drug response in Plasmodium.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2040-0) contains supplementary material, which is available to authorized users.