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.

The predisposition of parasites acquiring artemisinin resistance still remains unclear beyond the mutations in Pfk13 gene and modulation of the unfolded protein response pathway. To explore the chain of casualty underlying artemisinin resistance, we reanalyze 773 P. falciparum isolates from TRACI-study integrating TWAS, GWAS, and eQTL analyses. We find the majority of P. falciparum parasites are transcriptomically converged within each geographic site with two broader physiological profiles across the Greater Mekong Subregion (GMS). We report 8720 SNP-expression linkages in the eastern GMS parasites and 4537 in the western. The minimal overlap between them suggests differential gene regulatory networks facilitating parasite adaptations to their unique host environments. Finally, we identify two genetic and physiological backgrounds associating with artemisinin resistance in the GMS, together with a farnesyltransferase protein and a thioredoxin-like protein which may act as vital intermediators linking the Pfk13 C580Y mutation to the prolonged parasite clearance time.

The emergence and spread of artemisinin-resistant Plasmodium falciparum, first in the Greater Mekong Subregion (GMS), and now in East Africa, is a major threat to global malaria elimination ambitions. To investigate the artemisinin resistance mechanism, transcriptome analysis was conducted of 577 P. falciparum isolates collected in the GMS between 2016–2018. A specific artemisinin resistance-associated transcriptional profile was identified that involves a broad but discrete set of biological functions related to proteotoxic stress, host cytoplasm remodelling, and REDOX metabolism. The artemisinin resistance-associated transcriptional profile evolved from initial transcriptional responses of susceptible parasites to artemisinin. The genetic basis for this adapted response is likely to be complex.

BackgroundThe dynamics of histone modifications in Plasmodium falciparum indicates the existence of unique mechanisms that link epigenetic factors with transcription. Here, we studied the impact of acetylated histone code on transcriptional regulation during the intraerythrocytic developmental cycle (IDC) of P. falciparum.ResultsUsing a dominant-negative transgenic approach, we showed that acetylations of histone H4 play a direct role in transcription. Specifically, these histone modifications mediate an inverse transcriptional relationship between the factors of cell proliferation and host–parasite interaction. Out of the four H4 acetylations, H4K8ac is likely the rate-limiting, regulatory step, which modulates the overall dynamics of H4 posttranslational modifications. H4K8ac exhibits maximum responsiveness to HDAC inhibitors and has a highly dynamic distribution pattern along the genome of P. falciparum during the IDC. Moreover, H4K8ac functions mainly in the euchromatin where its occupancy shifts from intergenic regions located upstream of 5′ end of open reading frame into the protein coding regions. This shift is directly or indirectly associated with transcriptional activities at the corresponding genes. H4K8ac is also active in the heterochromatin where it stimulates expression of the main antigenic gene family (var) by its presence in the promoter region.ConclusionsOverall, we demonstrate that H4K8ac is a potential major regulator of chromatin-linked transcriptional changes during P. falciparum life cycle which is associated not only with euchromatin but also with heterochromatin environment. This is potentially a highly significant finding that suggests a regulatory connection between growth and parasite–host interaction both of which play a major role in malaria parasite virulence.Electronic supplementary materialThe online version of this article (doi:10.1186/s13072-017-0147-z) contains supplementary material, which is available to authorized users.

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.