Artemisinin resistance in Plasmodium falciparum threatens global efforts to control and eliminate malaria. Polymorphisms in the kelch domain–carrying protein K13 are associated with artemisinin resistance, but the underlying molecular mechanisms are unknown. We analyzed the in vivo transcriptomes of 1043 P. falciparum isolates from patients with acute malaria and found that artemisinin resistance is associated with increased expression of unfolded protein response (UPR) pathways involving the major PROSC and TRiC chaperone complexes. Artemisinin-resistant parasites also exhibit decelerated progression through the first part of the asexual intraerythrocytic development cycle. These findings suggest that artemisinin-resistant parasites remain in a state of decelerated development at the young ring stage, whereas their up-regulated UPR pathways mitigate protein damage caused by artemisinin. The expression profiles of UPR-related genes also associate with the geographical origin of parasite isolates, further suggesting their role in emerging artemisinin resistance in the Greater Mekong Subregion.

Clonally variant expression of surface antigens allows the malaria parasite Plasmodium falciparum to evade immune recognition during blood stage infection and secure malaria transmission. We demonstrate that heterochromatin protein 1 (HP1), an evolutionary conserved regulator of heritable gene silencing, controls expression of numerous P. falciparum virulence genes as well as differentiation into the sexual forms that transmit to mosquitoes. Conditional depletion of P. falciparum HP1 (PfHP1) prevents mitotic proliferation of blood stage parasites and disrupts mutually exclusive expression and antigenic variation of the major virulence factor PfEMP1. Additionally, PfHP1-dependent regulation of PfAP2-G, a transcription factor required for gametocyte conversion, controls the switch from asexual proliferation to sexual differentiation, providing insight into the epigenetic mechanisms underlying gametocyte commitment. These findings show that PfHP1 is centrally involved in clonally variant gene expression and sexual differentiation in P. falciparum and have major implications for developing antidisease and transmission-blocking interventions against malaria.

Due to their remarkable parasitocidal activity, artemisinins represent the key components of first-line therapies against Plasmodium falciparum malaria. However, the decline in efficacy of artemisinin-based drugs jeopardizes global efforts to control and ultimately eradicate the disease. To better understand the resistance phenotype, artemisinin-resistant parasite lines were derived from two clones of the 3D7 strain of P. falciparum using a selection regimen that mimics how parasites interact with the drug within patients. This long term in vitro selection induced profound stage-specific resistance to artemisinin and its relative compounds. Chemosensitivity and transcriptional profiling of artemisinin-resistant parasites indicate that enhanced adaptive responses against oxidative stress and protein damage are associated with decreased artemisinin susceptibility. This corroborates our previous findings implicating these cellular functions in artemisinin resistance in natural infections. Genomic characterization of the two derived parasite lines revealed a spectrum of sequence and copy number polymorphisms that could play a role in regulating artemisinin response, but did not include mutations in pfk13, the main marker of artemisinin resistance in Southeast Asia. Taken together, here we present a functional in vitro model of artemisinin resistance that is underlined by a new set of genetic polymorphisms as potential genetic markers.

Epigenetic mechanisms are emerging as one of the major factors of the dynamics of gene expression in the human malaria parasite, Plasmodium falciparum. To elucidate the role of chromatin remodeling in transcriptional regulation associated with the progression of the P. falciparum intraerythrocytic development cycle (IDC), we mapped the temporal pattern of chromosomal association with histone H3 and H4 modifications using ChIP-on-chip. Here, we have generated a broad integrative epigenomic map of twelve histone modifications during the P. falciparum IDC including H4K5ac, H4K8ac, H4K12ac, H4K16ac, H3K9ac, H3K14ac, H3K56ac, H4K20me1, H4K20me3, H3K4me3, H3K79me3 and H4R3me2. While some modifications were found to be associated with the vast majority of the genome and their occupancy was constant, others showed more specific and highly dynamic distribution. Importantly, eight modifications displaying tight correlations with transcript levels showed differential affinity to distinct genomic regions with H4K8ac occupying predominantly promoter regions while others occurred at the 5′ ends of coding sequences. The promoter occupancy of H4K8ac remained unchanged when ectopically inserted at a different locus, indicating the presence of specific DNA elements that recruit histone modifying enzymes regardless of their broad chromatin environment. In addition, we showed the presence of multivalent domains on the genome carrying more than one histone mark, highlighting the importance of combinatorial effects on transcription. Overall, our work portrays a substantial association between chromosomal locations of various epigenetic markers, transcriptional activity and global stage-specific transitions in the epigenome.

Historically seen as a benign disease, it is now becoming clear that Plasmodium vivax can cause significant morbidity. Effective control strategies targeting P. vivax malaria is hindered by our limited understanding of vivax biology. Here we established the P. vivax transcriptome of the Intraerythrocytic Developmental Cycle (IDC) of two clinical isolates in high resolution by Illumina HiSeq platform. The detailed map of transcriptome generates new insights into regulatory mechanisms of individual genes and reveals their intimate relationship with specific biological functions. A transcriptional hotspot of vir genes observed on chromosome 2 suggests a potential active site modulating immune evasion of the Plasmodium parasite across patients. Compared to other eukaryotes, P. vivax genes tend to have unusually long 5′ untranslated regions and also present multiple transcription start sites. In contrast, alternative splicing is rare in P. vivax but its association with the late schizont stage suggests some of its significance for gene function. The newly identified transcripts, including up to 179 vir like genes and 3018 noncoding RNAs suggest an important role of these gene/transcript classes in strain specific transcriptional regulation.