The variant antigen, Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1), expressed on the surface of P. falciparum infected Red Blood Cells (iRBCs) is a critical virulence factor for malaria1. Each parasite encodes 60 antigenically distinct var genes encoding PfEMP1s, but during infection the clonal parasite population expresses only one gene at a time before switching to the expression of a new variant antigen as an immune evasion mechanism to avoid the host’s antibody responses2,3. The mechanism by which 59 of the 60 var genes are silenced remains largely unknown4–7. Here we show that knocking out the P. falciparum variant-silencing SET gene (PfSETvs), which encodes an ortholog of Drosophila melanogaster ASH1 and controls histone H3 lysine 36 trimethylation (H3K36me3) on var genes, results in the transcription of virtually all var genes in the single parasite nuclei and their expression as proteins on the surface of individual iRBCs. PfSETvs-dependent H3K36me3 is present along the entire gene body including the transcription start site (TSS) to silence var genes. With low occupancy of PfSETvs at both the TSS of var genes and the intronic promoter, expression of var genes coincides with transcription of their corresponding antisense long non-coding RNA (lncRNA). These results uncover a novel role of the PfSETvs-dependent H3K36me3 in silencing var genes in P. falciparum that might provide a general mechanism by which orthologs of PfSETvs repress gene expression in other eukaryotes. PfSETvs knockout parasites expressing all PfEMP1s may also be applied to the development of a malaria vaccine.

Malaria, an infectious disease caused by Plasmodium parasites, still accounts for amounts of deaths annually in last decades. Despite the significance of Plasmodium falciparum as a model organism of malaria parasites, our understanding of gene expression of this parasite remains largely elusive since lots of progress on its genome and transcriptome are based on assembly with short sequencing reads. Herein, we report the new version of transcriptome dataset containing all full-length transcripts over the whole asexual blood stages by adopting a full-length sequencing approach with optimized experimental conditions of cDNA library preparation. We have identified a total of 393 alternative splicing (AS) events, 3,623 long non-coding RNAs (lncRNAs), 1,555 alternative polyadenylation (APA) events, 57 transcription factors (TF), 1,721 fusion transcripts in P. falciparum. Furthermore, the shotgun proteome was performed to validate the full-length transcriptome of P. falciparum. More importantly, integration of full-length transcriptomic and proteomic data identified 160 novel small proteins in lncRNA regions. Collectively, this full-length transcriptome dataset with high quality and accuracy and the shotgun proteome analyses shed light on the complex gene expression in malaria parasites and provide a valuable resource for related functional and mechanistic researches on P. falciparum genes.

The heterochromatin environment plays a central role in silencing genes associated with the malaria parasite’s development, survival in the host, and transmission to the mosquito vector. However, the underlying mechanism regulating the dynamic chromatin structure is not understood yet. Here, we have uncovered that Plasmodium falciparum Rrp6, an orthologue of eukaryotic RNA exosome-associated RNase, controls the silencing of heterochromatic genes. PfRrp6 knockdown disrupted the singular expression of the GC-rich ncRNA RUF6 family, a known critical regulator of virulence gene expression, through the stabilization of the nascent transcripts. Mechanistic investigation showed that the accumulation of the multiple RUF6 ncRNAs triggered local chromatin remodeling in situ, which activated their adjacent var genes. Strikingly, chromatin isolation by RNA purification analysis (ChIRP-seq) revealed that a remarkable RUF6 ncRNA had interacted with distal heterochromatin regions directly and stimulated a global derepression effect on heterochromatic genes, including all variant gene families and the sexual commitment-associated regulator ap2-g gene. Collectively, Rrp6 appears to conduct the epigenetic surveillance of heterochromatic gene expression through controlling RUF6 levels, thereby securing antigenic variation and sexual commitment of malaria parasites during the infection of the host. IMPORTANCE Malaria remains a major public health and economic burden. The heterochromatin environment controls the silencing of genes associated with the fate of malaria parasites. Previous studies have demonstrated that a group of GC-rich ncRNAs (RUF6) is associated with the mutually exclusive expression of var genes, but the underlying mechanisms remain elusive. Here, through a series of genetic manipulation and genome-wide multiomics analysis, we have identified the plasmodial orthologue of RNA exosome-associated Rrp6 as an upstream regulator of RUF6 expression and revealed that the dysregulation of RUF6 upon Rrp6 knockdown triggered local chromatin alteration, thereby activating most heterochromatic genes via direct interaction of RUF6 and distal gene loci. This finding not only uncovered the in-depth mechanism of RUF6-mediated regulation of heterochromatic genes but also identified Rrp6 as a novel regulator of gene expression in human malaria parasites, which provides a new target for developing intervention strategies against malaria.

Antisense transcription emerges as a key regulator of important biological processes in the human malaria parasite Plasmodium falciparum. RNA-processing factors, however, remain poorly characterized in this pathogen. Here, we purified the multiprotein RNA exosome complex of malaria parasites by affinity chromatography, using HA-tagged PfRrp4 and PfDis3 as the ligands. Seven distinct core exosome subunits (PfRrp41, PfMtr3, PfRrp42, PfRrp45, PfRrp4, PfRrp40, PfCsl4) and two exoribonuclease proteins PfRrp6 and PfDis3 are identified by mass spectrometry. Western blot analysis detects Dis3 and Rrp4 predominantly in the cytoplasmic fraction during asexual blood stage development. An inducible gene knock out of the PfDis3 subunit reveals the upregulation of structural and coding RNA, but the vast majority belongs to antisense RNA. Furthermore, we detect numerous types of cryptic unstable transcripts (CUTs) linked to virulence gene families including antisense RNA in the rif gene family. Our work highlights the limitations of steady-state RNA analysis to predict transcriptional activity and link the RNA surveillance machinery directly with post-transcriptional control and gene expression in malaria parasites.

BackgroundYunnan and Hainan provinces are the two major endemic regions for Plasmodiumfalciparum malaria in China. However, few studies have investigated the characteristics of this parasite. Therefore, this study aimed to evaluate the genetic diversity and population structure of P. falciparum to predict the geographic origin of falciparum malaria.MethodsThirteen highly polymorphic microsatellite loci were studied to estimate the genetic diversity and population structure of 425 P. falciparum isolates obtained from blood samples collected from Yunnan and Hainan provinces of South China. The isolates were analysed for genetic diversity, linkage disequilibrium, and population structure. The parasite populations were clustered into two subgroups (i.e., Yunnan and Hainan) and a classification algorithm was used to identify molecular markers for classifying the P. falciparum populations.ResultsAll 13 microsatellite loci were highly polymorphic, with the number of alleles per locus varying from 5 to 20. The mean expected heterozygosity (He) in Yunnan and Hainan was 0.766 ± 0.036 and 0.677 ± 0.039, respectively, revealing a moderate high level of genetic diversity. Significant linkage disequilibrium was found for some regions of Yunnan (Lazan county and Xishuangbanna region) and Hainan (Dongfang city and Sanya city) province. According to the classification algorithm, a combination of three microsatellites could be used as a discriminatory marker to identify the origin of P. falciparum isolates.ConclusionsThe results on the genetic structure of P. falciparum populations from South China provide a basis for developing a genetic marker-based tool to trace the source of the parasite infections and consequently improve malaria control and elimination strategies.Electronic supplementary materialThe online version of this article (doi:10.1186/s12936-015-0786-0) contains supplementary material, which is available to authorized users.