Transcriptome-module phenotype association study implicates extracellular vesicles biogenesis in Plasmodium falciparum artemisinin resistance

Plasmodium falciparum malaria is still an important disease in sub-Saharan Africa (sSA). Great strides have been made in its control spear-headed by artemisinin (ART)-based combination therapies (ACTs). However, concerns about the imminent spread of ART-resistant (ARTr) malaria parasites to sSA threaten gains already made. Attempts to mitigate this risk have highlighted the need to discover novel P. falciparum drug targets. Therefore, studies to deepen our understanding of the biology of P. falciparum are needed. The role of extracellular vesicles (EVs) in the biology of malaria parasites is not fully understood. Recently, the ART resistance-associated transcriptional profile has been reported to involve several biological processes connected to vesicular trafficking, proteotoxic stress, erythrocyte remodelling, and mitochondrial metabolism. We explored a role for EVs in developing the P. falciparum ARTr phenotype using bulk RNA sequencing of unsynchronized parasite cultures under untreated, 0.1% dimethyl sulfoxide and 700nM dihydroartemisinin treated conditions for six hours. As pathway and gene ontology analysis is limited in its curated knowledge repertoire on EVs biogenesis in P. falciparum, we used a modular (gene set) analysis approach to explore whether an EVs biogenesis module is associated with the ARTr phenotype in P. falciparum. We first generated well-defined EVs modules of interest and used statistical tools to determine differences in their expression among the parasite and treatment conditions. Then we used gene set enrichment analysis to determine the strength of the association between each EVs module of interest and the ARTr phenotype. This transcriptome-module phenotype association study (TMPAS) represents a well-powered approach to making meaningful discoveries out of bulk gene expression data. We identified four EVs module of interest and report that one module representing gene sets with correlated expression to PF3D7_1441800 – involved with EVs biogenesis in P. falciparum - is associated with the ARTr phenotype (R539T_DHA_treated versus R539T_untreated: normalized enrichment score (NES) = 1.1830174, FDR q-value < 0.25; C580R_DHA_treated versus C580R_untreated: NES = 1.2457103, FDR q-value < 0.25). PF3D7_1441800 has been reported to reduce EVs production when knocked out in P. falciparum. Altogether, our findings suggest a role for EVs in developing ART resistance and warrant further studies interrogating this association.