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.

Stochastic expression of invasion genes in Plasmodium falciparum schizonts.

Genetically identical cells are known to exhibit differential phenotypes in the same environmental conditions. These phenotypic variants are linked to transcriptional stochasticity and have been shown to contribute towards adaptive flexibility of a wide range of unicellular organisms. Here, we investigate transcriptional heterogeneity and stochastic gene expression in Plasmodium falciparum by performing the quasilinear multiple annealing and looping based amplification cycles (MALBAC) based amplification and single cell RNA sequencing of blood stage schizonts. Our data reveals significant transcriptional variations in the schizont stage with a distinct group of highly variable invasion gene transcripts being identified. Moreover, the data reflects several diversification processes including putative developmental "checkpoint"; transcriptomically distinct parasite sub-populations and transcriptional switches in variable gene families (var, rifin, phist). Most of these features of transcriptional variability are preserved in isogenic parasite cell populations (albeit with a lesser amplitude) suggesting a role of epigenetic factors in cell-to-cell transcriptional variations in human malaria parasites. Lastly, we apply quantitative RT-PCR and RNA-FISH approach and confirm stochastic expression of key invasion genes, such as, msp1, msp3, msp7, eba181 and ama1 which represent prime candidates for invasion-blocking vaccines.

Artemisinin resistance in the malaria parasite, Plasmodium falciparum, originates from its initial transcriptional response.

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.

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.

A Novel Chemically Differentiated Mouse Embryonic Stem Cell-Based Model to Study Liver Stages of Plasmodium berghei.

Asymptomatic and obligatory liver stage (LS) infection of Plasmodium parasites presents an attractive target for antimalarial vaccine and drug development. Lack of robust cellular models to study LS infection has hindered the discovery and validation of host genes essential for intrahepatic parasite development. Here, we present a chemically differentiated mouse embryonic stem cell (ESC)-based LS model, which supports complete development of Plasmodium berghei exoerythrocytic forms (EEFs) and can be used to define new host-parasite interactions. Using our model, we established that host Pnpla2, coding for adipose triglyceride lipase, is dispensable for P. berghei EEF development. In addition, we also evaluated in-vitro-differentiated human hepatocyte-like cells (iHLCs) to study LS of P. berghei and found it to be a sub-optimal infection model. Overall, our results present a new mouse ESC-based P. berghei LS infection model that can be utilized to study the impact of host genetic variation on parasite development.

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

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.

Histone 4 lysine 8 acetylation regulates proliferation and host-pathogen interaction in Plasmodium falciparum.

BACKGROUND: The 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. RESULTS: Using 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. CONCLUSIONS: Overall, 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.

Plasmodium falciparum origin recognition complex subunit 1 (PfOrc1) functionally complements Δsir3 mutant of Saccharomyces cerevisiae.

Telomere position effect efficiently controls silencing of subtelomeric var genes, which are involved in antigenic variation in human malaria parasite Plasmodium falciparum. Although, PfOrc1 has been found to be associated with PfSir2 in the silencing complex, its function in telomere silencing remained uncertain especially due to an apparent lack of BAH domain at its amino-terminal region. Here we report that PfOrc1 possesses a Sir3/Orc1 like silencing activity. Using yeast as a surrogate organism we have shown that PfOrc1 could complement yeast Sir3 activity during telomere silencing in a Sir2 dependent manner. By constructing a series of chimera between PfOrc1 and ScSir3 we have observed that the amino-terminal domain of PfOrc1 harbors silencing activity similar to that present in the amino-terminal domain of ScSir3. We further generated several amino-terminal deletion mutants to dissect out such silencing activity and found that the first seventy amino acids at the amino-terminal domain are dispensable for its activity. Thus our results strongly supports that PfOrc1 may have a role in telomere silencing in this parasite. This finding will help to decipher the mechanism of telomere position effect in P. falciparum.