ABSTRACT: Transcriptional profiling of P. falciparum cultures treated with cyclohexylamine over time (18, 25 and 30 hours post invasion) A control experiment was also set up in which P. falciparum was not treated with cyclohexylamine, and samples were taken at 18, 25 and 30 h post invasion). Background Plasmodium falciparum, the causative agent of severe human malaria, has evolved to become resistant to previously successful antimalarial chemotherapies, most notably chloroquine and the antifolates. The prevalence of resistant strains has necessitated the discovery and development of new chemical entities with novel modes-of-action. Although much effort has been invested in the creation of analogues based on existing drugs and the screening of chemical and natural compound libraries, a crucial shortcoming in current Plasmodial drug discovery efforts remains the lack of an extensive set of novel, validated drug targets. A requirement of these targets (or the pathways in which they function) is that they prove essential for parasite survival. The polyamine biosynthetic pathway, responsible for the metabolism of highly abundant amines crucial for parasite growth, proliferation and differentiation, is currently under investigation as an antimalarial target. Chemotherapeutic strategies targeting this pathway have been successfully utilized for the treatment of Trypanosomes causing West African sleeping sickness. In order to further evaluate polyamine depletion as possible antimalarial intervention, the consequences of inhibiting P. falciparum spermidine synthase (PfSpdSyn) were examined on a morphological, transcriptomic, proteomic and metabolic level. Results Morphological analysis of P. falciparum 3D7 following application of the PfSpdSyn inhibitor cyclohexylamine confirmed that parasite development was completely arrested at the early trophozoite stage. This is in contrast to untreated parasites which progressed to late trophozoites at comparable time points. Global gene expression analyses confirmed a transcriptional arrest in the parasite. Several of the differentially expressed genes mapped to the polyamine biosynthetic and associated metabolic pathways. Differential expression of corresponding parasite proteins involved in polyamine biosynthesis was also observed. Most notably, uridine phosphorylase, adenosine deaminase, lysine decarboxylase (LDC) and S-adenosylmethionine synthetase were differentially expressed at the transcript and/or protein level. Several genes in associated metabolic pathways (purine metabolism and various methyltransferases) were also affected. The specific nature of the perturbation was additionally reflected by changes in polyamine metabolite levels. Conclusions This study details the malaria parasite’s response to PfSpdSyn inhibition on the transcriptomic, proteomic and metabolic levels. The results corroborate and significantly expand previous functional genomics studies relating to polyamine depletion in this parasite. Moreover, they confirm the role of transcriptional regulation in P. falciparum, particularly in this pathway. The findings promote this essential pathway as a target for antimalarial chemotherapeutic intervention strategies. Keywords: Time course experiment in response to a drug treatment Reference design. All timepoints compared with t = 18 hours (untreated) post invasion. Two biological replicates and one technical sample per treatment and per time point. A reference design was employed for array hybridisation, utilising the URR pool described previously. All solvent-control and drug-treated sampled were hybridised to Operon slides, along with the URR. For each time point and each untreated/treated sample, three microarray slides were processed, such that a total of eighteen slides were processed in the study. Two independent cDNA samples (biological replicates) were prepared for each untreated and drug-treated sample at each time point. One of the biological replicate cDNA samples were additionally hybridised to a third slide (representing the technical replicate). GenePix results (gpr) files were generated using GenePix 6.0 (Molecular Devices) software, without normalization. For clustering analyses, results files were normalized with DNMAD (Diagnosis and Normalization for MicroArray Data) using print-tip loess. The normalized values were subsequently downloaded and analyzed with the Multiexperiment Viewer (MeV) in the TM4 software suite. Hierarchical Clustering (HCL, average linkage) was performed to estimate technical and biological variation between samples and at which point cytostasis most likely occurred for comparative purposes in downstream analyses. Intensity data for individual slides were imported into LIMMA (linear models for microarray data) in the R computing environment. Pre- and post-normalization diagnostic plots were performed using MARRAY. Data from each microarray slide was normalized using print-tip loess. Data between microarrays was normalized using rquantile normalisation. Pearson correlations were computed in ExCel to estimate variation between technical and biological replicates. Spots excluded from slide correlations and normalisation were those weighted by the limma script or flagged in the genepix results file (gpr). Additionally, spots termed Alien, Empty, Null and Operon Use Only were excluded from the correlation analyses. These spots were similarly excluded for correlations between untreated and treated samples at each time point following normalisation. Results from biological and slide replicates within each of the time points were collated, and linear models were computed to contrast gene expression between time points. A two-fold change in gene expression was used as cut-off, in conjunction with correction for false discovery (false discovery rate (FDR) = 5%). Normalised data was deposited in the Gene Expression Omnibus (GEO) database, number GSE18075. Analysis of differentially expressed genes was performed in MADIBA (Micro Array Data Interface for Biological Annotation).
