# Dynamic RNA profiling in Plasmodium falciparum synchronized blood stages exposed to lethal doses of artesunate.

Onguma Natalang, Emmanuel Bischoff, Guillaume Deplaine, Caroline Proux, Marie-Agnès Dillies, Odile Sismeiro, Ghislaine Guigon, Serge Bonnefoy, Jintana Patarapotikul, Odile Mercereau-Puijalon, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Onguma Natalang, Emmanuel Bischoff, Guillaume Deplaine, Caroline Proux, Marie-Agnès Dillies, et al. Dynamic RNA profiling in Plasmodium falciparum synchronized blood stages exposed to lethal doses of artesunate. BMC Genomics, 2008, 9 (1), pp.388. 10.1186/1471-2164-9-388. pasteur-01472184

HAL Id: pasteur-01472184

https://pasteur.hal.science/pasteur-01472184

Submitted on 20 Feb 2017

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
Research article

Open Access

# Dynamic RNA profiling in *Plasmodium falciparum* synchronized blood stages exposed to lethal doses of artesunate

Onguma Natalang<sup>1,3</sup>, Emmanuel Bischoff<sup>2</sup>, Guillaume Deplaine<sup>1</sup>, Caroline Proux<sup>2</sup>, Marie-Agnès Dillies<sup>2</sup>, Odile Sismeiro<sup>2</sup>, Ghislaine Guigon<sup>2</sup>, Serge Bonnefoy<sup>1</sup>, Jintana Patarapotikul<sup>3</sup>, Odile Mercereau-Puijalon<sup>1</sup>, Jean-Yves Coppée<sup>2</sup> and Peter H David*<sup>1</sup>

Address: <sup>1</sup>Institut Pasteur, Unité d'Immunologie Moléculaire des Parasites; CNRS URA 2581, 28 Rue du Docteur Roux, F-75724, Paris, Cedex 15, France, <sup>2</sup>Institut Pasteur, Plate-Forme 2 - Puces à ADN, 28 Rue du Docteur Roux, F-75724, Paris, Cedex 15, France and <sup>3</sup>Faculty of Tropical Medicine, Mahidol University, 420/6 Rajvithi Road, Bangkok, 10400, Thailand

Email: Onguma Natalang - onguma382@yahoo.com; Emmanuel Bischoff - bischoff@pasteur.fr; Guillaume Deplaine - gdep@pasteur.fr; Caroline Proux - karo@pasteur.fr; Marie-Agnès Dillies - mdillies@pasteur.fr; Odile Sismeiro - osisme@pasteur.fr; Ghislaine Guigon - gguigon@pasteur.fr; Serge Bonnefoy - sbf@pasteur.fr; Jintana Patarapotikul - tmjpt@mahidol.ac.th; Odile Mercereau-Puijalon - omp@pasteur.fr; Jean-Yves Coppée - jcoppee@pasteur.fr; Peter H David* - pdavid@pasteur.fr

Published: 18 August 2008

Received: 14 February 2008

BMC Genomics 2008, 9:388 doi:10.1186/1471-2164-9-388

Accepted: 18 August 2008

This article is available from: http://www.biomedcentral.com/1471-2164/9/388

© 2008 Natalang et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## Abstract

**Background:** Translation of the genome sequence of *Plasmodium* sp. into biologically relevant information relies on high through-put genomics technology which includes transcriptome analysis. However, few studies to date have used this powerful approach to explore transcriptome alterations of *P. falciparum* parasites exposed to antimalarial drugs.

**Results:** The rapid action of artesunate allowed us to study dynamic changes of the parasite transcriptome in synchronous parasite cultures exposed to the drug for 90 minutes and 3 hours. Developmentally regulated genes were filtered out, leaving 398 genes which presented altered transcript levels reflecting drug-exposure. Few genes related to metabolic pathways, most encoded chaperones, transporters, kinases, Zn-finger proteins, transcription activating proteins, proteins involved in proteasome degradation, in oxidative stress and in cell cycle regulation. A positive bias was observed for over-expressed genes presenting a subtelomeric location, allelic polymorphism and encoding proteins with potential export sequences, which often belonged to subtelomeric multi-gene families. This pointed to the mobilization of processes shaping the interface between the parasite and its environment. In parallel, pathways were engaged which could lead to parasite death, such as interference with purine/pyrimidine metabolism, the mitochondrial electron transport chain, proteasome-dependent protein degradation or the integrity of the food vacuole.

**Conclusion:** The high proportion of over-expressed genes encoding proteins exported from the parasite highlight the importance of extra-parasitic compartments as fields for exploration in drug research which, to date, has mostly focused on the parasite itself rather than on its intra and extra erythrocytic environment. Further work is needed to clarify which transcriptome alterations observed reflect a specific response to overcome artesunate toxicity or more general perturbations on the path to cellular death.

<footnote>* Corresponding author †Equal contributors</footnote>
## Background

Efforts at extracting biologically meaningful information from the genome sequence of *Plasmodium* sp. are fueled by the necessity to find new methods of malaria control. *Plasmodium falciparum* resistance to most common antimalarials such as chloroquine, sulfadoxine or pyrimethamine, is now-days widespread, and resistance to other antimalarials is increasing alarmingly. There is recent evidence for emerging in vitro resistance to one of the major groups of drugs, artemisinins [1] and treatment failures to recently introduced drug combinations including artemisinin derivatives [2,3]. Research priorities include the discovery of new drugs, the understanding of drugs' mode of action and resistance mechanisms. Whatever the drug's mode of action, a better appraisal of the biological processes taking place in a parasite dying within an erythrocyte is needed to uncover additional drug targets and potentiate the effect of existing weapons. Towards reaching such goals, much hope rests on high through-put technologies such as DNA microarrays aimed at the study of gene expression at the genome scale and proteomics. Both approaches are complementary, and each has its own limitations. Proteomic approaches are complicated by solubility issues, detection of low abundance proteins and post-translational alterations. Both proteome and transcriptome approaches investigate steady state levels, a net sum of synthesis and decay of the gene product. Transcriptome studies are limited by the uncertain correlation of steady state RNA with protein levels. They nevertheless deliver a comprehensive and sensitive, genome-wide exploration of expression profiling. This proved invaluable in deciphering the developmental cycle gene expression patterns in malaria parasites [4,5].

Few studies have shown that the transcriptome of *P. falciparum* can be altered by exposure of the parasite to a drug. Analysis of the parasite transcriptome under doxycyclin showed that apicoplast gene expression was deeply altered [6]. This study is so far the only published example of a specific transcriptome response to a drug related to the drug's known target-pathway. No significant changes of the mRNA levels for enzymes involved in the lipid biosynthesis pathway could be evidenced upon exposure for 24-36 hours to the antimalarial choline analog T4, which targets the inhibition of phosphatidyl choline biosynthesis [7]. In a recent study of the transcriptome of asynchrónous *P. falciparum* cultures under chloroquine pressure [8], around 600 genes were differentially expressed in presence of the drug, but only 38 of these were observed in two different experiments. Such poor reproducibility of results along with the low expression ratios observed for these 38 genes was interpreted as reflecting limited reactivity of the parasite transcriptome to environmental stimuli and the possible importance of post-transcriptional gene expression regulation in *Plasmodium*. However, the

observed discrepancies may relate to different proportions of developmental stages in the 2 populations explored. Indeed, parasite maturation stage is an essential variable to control, as a majority of *P. falciparum* genes are differentially expressed during the 48 hour erythrocytic cycle [4]. Such developmentally regulated control of gene expression hinders the interpretation of results if drug treatment exerts a retardation effect on the parasite erythrocytic cycle.

We have used here a novel strategy and analysed the alterations of the parasite transcriptome in synchronous parasite cultures shortly after exposure to a lethal dose of artesunate, a rapidly acting drug. Artemisinin derivatives are potent antimalarials that are the cornerstone of currently recommended drug combination to treat *P. falciparum* malaria [9]. Artemisinins, which are endoperoxide-containing sesquiterpene lactones are structurally distinct from all other anti-malarials. Importantly they are the most rapidly acting antimalarials known to date, which moreover are active on young as well mature erythrocytic stages [10]. These unique properties allowed the kinetic analysis of the events occurring in synchronous cultures of parasites exposed to a lethal drug concentration at different time points along the intraerythrocytic developmental cycle. We explored the dynamics of transcriptome alteration during the process leading to parasite death in different parasite stages. This identified 398 genes with dynamic transcriptome alterations after exposure to artesunate.

## Results

### Choice of experimental design

In order to work with homogeneous parasite populations, tightly synchronized parasites were used in all experiments, along with an artesunate concentration which irreversibly damages 100% of the parasites within 3 hours. The lethal dose of artesunate, defined here as the minimal concentration needed for a 3 hour exposure to completely prevent reinvasion (see "Materials and methods" section), was 780 nM. Importantly, it corresponds to 100-1000 times the IC₅₀, this dose remains physiological, being in the range of peak plasma concentration reached in patients administered artesunate (1300-8400 nM) [11].

Each experiment consisted in an artesunate-free culture (control) and an artesunate treated culture. For each time point explored, cDNA from parasites with and without the drug were labeled with different cyanines, mixed and hybridized. To analyse dynamic transcriptome alterations, 2 pilot experiments were performed in which RNA was harvested at 30, 90, 180 minutes and 10 hours of incubation with the drug. No significant changes in gene expression could be detected at 30 minutes. In contrast, after 10 hour incubation time, major transcriptome alter-
ations occurred under artesunate, but their analysis was un-interpretable due to drug-induced massive slowing-down of parasite development (data not shown). We thus decided to focus on 2 drug-exposure time, namely 90 minutes and 3 hours after drug addition.

Since artesunate is active on ring stages as well as on older trophozoite stages, we decided to explore different time points along the erythrocytic cycle and search for genes affected at each time point. Therefore, 5 drug treatment experiments, staggered between 20 hours and 30 hours of parasite development, were performed for the comparison of drug versus no drug at 90 minutes and 3 hours.

### Transcript changes under artesunate exposure

We looked for genes displaying similar changes across all experiments, with the hypothesis that responses to artesunate would not be stage-dependent within the explored 10 hour window of the developmental cycle. The results of 5 experiments performed at time points staggered between 20-30 hours post-invasion, were analyzed using ANOVA with a Holm-Bonferroni method of p-value adjustment, p < 0.05. This led to a list of genes with similar ratios of expression (i. e. low variance). From this list, only genes differentially expressed at 3 hours under artesunate and not after 3 hours in the absence of the drug were retained (p = 1). It is important to stress that what we refer to, for convenience, as over- or under-expression in fact refers to higher or lower steady state RNA levels.

Due to the high number of microarrays analyzed (32 at 3 hours), power of the statistical analysis was such that expression log ratios as low as +/- 0.3 could be found for genes that displayed significant changes by ANOVA [for distribution of log ratios, see Additional file 1]. This appeared questionable from a biological standpoint and we chose to use a cut-off log-ratio of +/- 0.8. This value is lower than what is usually (and often arbitrarily) chosen in most transcriptome studies, i.e. +/- 1, but we felt it was justified given the large amount of data analyzed here and the low amplitude of fold-changes reported in previous *P. falciparum* transcriptome studies [8,12].

In addition, as genes of potential interest could be differentially expressed in the controls but more-so or differently under artesunate, we needed to identify these and filter out the developmentally regulated gene expression profiles. Indeed, the extent to which the drug slowed-down parasite development is unknown. Thus, different expression levels in artesunate-exposed vs. control (unexposed) cultures could reflect the direct effect of artesunate together with a possible difference in growth rate/developmental stage. To identify the developmentally regulated genes in the artesunate-free control culture, the 3 hours RNA was hybridized against the time 0 RNA (3 con-

trol experiments). Even within such a short time window, 81 genes showed differential expression in the control culture. Filtering-out this potential confounding "slowing-down effect" established a second list of genes differentially expressed in the artesunate experiments and in the control experiments, but with significantly different expression-ratios under these two conditions (p < 0.05) [see Additional file 2 for the list of differentially expressed genes].

This resulted in the following: out of the 4703 genes analyzed, 398 (approx 8.5%) were differentially expressed after 3 hours in presence of artesunate, 244 over-expressed (approx 5%) and 154 (approx 3.5%) under-expressed. Of these 398 genes differentially expressed at 3 hours, 42 showed significant expression changes at 90 minutes, presumably reflecting an early response to the drug.

Our technically demanding and relatively stringent approach may have generated false negatives, but we wanted to be as confident as possible about the list of genes selected as being differentially expressed. This confidence was supported by the results of qRT-PCR performed with 12 over-expressed and 4 under-expressed genes, which confirmed the microarray data (Spearman correlation coefficient: 0.81) [see Additional file 3 for the comparison between microarray and qRT-PCR results].

### Classification of differentially expressed genes

Interesting features were observed in the group of differentially expressed genes. For genes expressed differentially at 3 hours, there was a strong positive bias in the sub-set of over-expressed genes towards polymorphic genes located in sub-telomeric position and carrying PEXEL/ HCT motifs [13] used for export beyond the parasitophorous vacuole into the erythrocyte cytoplasm, or to the red blood cell membrane/surface (Figure 1). This was determined after excluding from the analysis the strain-specific *var*, *stevor* and *rifin* multigene families – which are not adequately represented in our case since the parasite strain used differs from 3D7. Concomitant with this positive bias, there was a negative bias towards conserved genes in the over-expressed genes. In contrast, the group of under-expressed genes displayed a positive bias towards conserved genes and a negative bias towards polymorphic genes.

Functions (in most cases putative) were assigned to 213 of the 398 differentially expressed genes [12,14-21] [see Additional file 2]. Apart from lipid and purine/pyrimidine metabolism, few genes belonged to the metabolic pathways of the parasite (as defined by Ginsburg [16]). Most genes with altered expression under artesunate pressure were related to chaperones, transporters, cell cycle, kinases, Zn finger proteins, transcription activating pro-
A

<figure><img src="image_3.png" /></figure>

B

<figure><img src="image_4.png" /></figure>

C

<figure><img src="image_5.png" /></figure>

**Figure 1**
Association plots of gene expression with gene chromosomal localization, predicted export sequences and polymorphism. A: Bias in gene expression versus gene chromosomal localization (subtel = subtelomeric position of gene defined as <150 kb from telomere). B: Bias in gene expression versus presence of predicted export signals [13]. C: Bias in gene expression versus gene polymorphism (as defined by presence or absence of non synonymous SNPs, from analysis of 3D7, Dd2, HB3, Ghana1 and IT genomes). Bar width is proportional to the number of genes per category and bar height to the Pearson residuals for an independence model. Blue and red colors indicate significant bias.

teins, proteins involved in proteasome degradation, oxidative stress and in cell cycle regulation. A positive bias towards over-expression could be observed for genes related to oxidative stress, kinases, transcription associated proteins and chaperones, and towards under-expression for genes related to proteasome degradation and transporters (Figure 2). Slowing down of parasite development may be linked to the differential expression of three genes involved in cell cycle regulation [see Additional file 4]. Two genes were over-expressed, namely PFF1125 encoding a putative RNA binding protein mei2 homologue and PFL1855 encoding a putative cell cycle control protein, while PFA0345w, encoding a putative centrin, was under-expressed.

**Oxidative stress**

Oxidative stress is likely to occur in parasites submitted to artesunate, whether linked to the direct action of the drug or to the reduced capacity of the damaged parasite to inactivate free radicals. The antioxidant defense in *P. falciparum* involves both the glutathione and the thioredoxin

<figure><img src="image_6.png" /></figure>

**Figure 2**
Differential gene expression under artesunate exposure: transcript distribution among different gene functional categories. "Total" pie-chart: genes present in the 8 different color-coded functional categories as distributed in the genome represented on the DNA array. "over-exp." and "under-exp." pie-charts: distribution of over-expressed and under-expressed genes respectively, among the 8 different functional categories. p-value: p-value of chi-square test of over-expressed gene and under-expressed gene compared to the total genes in the genome represented on the DNA array.
system, the relative contribution of each remaining unknown [22]. Four genes involved in the parasite antioxidant defense [22,23], showed altered transcription profiles (Table 1) : gamma glutamyl cystein synthetase and glutathione synthetase (both involved in the glutathione system) were over-expressed. In contrast, the genes encoding protein disulfide isomerase related protein (potentially involved in oxidative protein folding, catalyzing both the oxidation and isomerization of disulfides on nascent polypeptides [24]) and glutathione peroxidase (shown to rather encode a thioredoxin peroxidase [25] and thus involved in the thioredoxin system) were under-expressed. The consecutive action of gamma glutamyl synthetase and glutathione synthetase leads to the production of reduced glutathione (GSH), which plays a pivotal role in the antioxidant defense through maintenance of the red-ox state of protein -SH moieties, the reduction of the noxious hydrogen and lipid peroxides and the extrusion of toxic compounds (including drugs) [22]. The predicted outcome of the over-expression of these 2 genes would be an increased supply of GSH in response to the oxidative stress induced by artesunate.

### Chaperones, protein trafficking and related genes

Genome mining on parasite chaperone analogues has led to the conclusion that approximately 2% of the *P. falciparum* genes encoded potential or confirmed chaperones [19]. Chaperone proteins have diverse functions, such as serving as transcription factors, regulating the cell cycle, being involved in protein degradation, binding damaged proteins, protein transport across membranes and protein folding. As expected in a parasite under severe injury, 13 genes encoding chaperone or chaperone-related proteins (Table 2) were found up-regulated at 3 hours of artesunate pressure, among which PfHsp70, PfHsp90 and 6 genes belonging to the Hsp40 chaperone machinery of *P. falciparum* [20]. The first four were already over-expressed at 90 minutes, along with PFD0080c (a member of the resa-like DnaJ family) and thus can be considered part of the early response to the drug. Interestingly, most (9/13), chaperone-encoding genes up-regulated under artesunate have potential export signals, suggesting possible export

of their corresponding proteins from the parasite to the erythrocyte cytoplasm or membrane. Chaperone involvement in trafficking of membrane or exported proteins has been highlighted in a recent analysis of the parasite chaperone network [21]. In particular, analysis of the interac-tome suggested that a member of the Hsp40 family (PF14_0700), which interacts with PfHsp70 (PF08_0054) and PfHsp90 (PF07_0029) could function as a co-chaperone in the Hsp90 complex of the parasite to transport erythrocyte membrane exported proteins such as Antigen 332 (PF11_0507). It is remarkable to note that PF14_0700, PF08_0054, PF07_0029 and PF11_0507 were all over-expressed under artesunate pressure [see Additional file 2]. This suggests artesunate-triggered altered intracellular trafficking, contributing to remodeling of the parasitized erythrocyte membrane, a conclusion further substantiated by the observed increased steady state levels of genes encoding proteins interacting with the erythrocyte cytoskeleton such as MESA (PFE0040c), several members of the RESA family and number of genes carrying PEXEL/HT motifs. Alterations of the intracellular trafficking is also suggested by altered transcription profiles of genes such as signal peptidase (MAL13P1.167) as well as transporters, including organellar importers (see below).

### Transporters

Eleven genes involved in transport were under-expressed (Table 3). The putative UDP galactose antiporter (PF11_0141) down-regulated at 3 hours, was the only transporter found down-regulated at 90 minutes, suggesting that alteration of its transcript levels is an early event. Four subunits of the vacuolar ATP synthase were under-expressed. This constitutes one of the few cases of committant regulation of individual components of a specific machinery or pathway in this study, suggesting that indeed the activity of the vacuolar ATP synthase was severely impaired. This enzyme contributes to acidification of vacuolar and organelle contents.

Only 4 genes encoding transporters were over-expressed, amongst which *pfmdr1* (PFE1150w), the P-glycoprotein

<table><caption>Table 1: Differentially expressed genes involved in antioxidant defence in artesunate-treated cells</caption><thead><tr><th>Gene ID</th><th>Description</th><th>ART 90 min</th><th>ART 3 h</th><th>HS</th><th>CQ</th></tr></thead><tbody><tr><td>PF10925w</td><td>gamma-glutamylcysteine synthetase</td><td></td><td>1,41</td><td>+</td><td></td></tr><tr><td>PFE0605c</td><td>glutathione synthetase</td><td></td><td>0,94</td><td></td><td></td></tr><tr><td>PFL0595c</td><td>glutathione peroxidase</td><td></td><td>-0,95</td><td></td><td></td></tr><tr><td>PF11_0352</td><td>protein disulfide isomerase related protein</td><td></td><td>-0,89</td><td>-</td><td></td></tr></tbody></table>

Genes coding for proteins involved in antioxidant defence defined in [23].

Gene ID: identifier as found in PlasmoDB.

Art columns: log ratios of gene expression under artesunate exposure at 90 minutes and 3 hours.

HS column: over (+) or under (-) expression.

Results were compared with transcriptome modifications induced by heat-shock (HS) [12] and chloroquine (CQ) [8].
<table><caption>Table 2: Differentially expressed chaperone and chaperone-related genes in artesunate-treated cells</caption><thead><tr><th>Gene ID</th><th>Description</th><th>Pexel</th><th>family</th><th>ref</th><th>ART 90 mn</th><th>ART 3 h</th><th>HS</th><th>CQ</th></tr></thead><tbody><tr><td>PFE1605w</td><td>DNAJ protein</td><td>yes</td><td>resa like/dnaJ</td><td>1, 3</td><td></td><td>1,01</td><td></td><td></td></tr><tr><td>PFD0080c</td><td>hypothetical protein, conserved in P.falciparum</td><td>yes</td><td>resa like</td><td>3</td><td>1,31</td><td>1,82</td><td>+</td><td></td></tr><tr><td>PF10_0378</td><td>DnaJ protein, putative</td><td>yes</td><td>Hsp40, type III</td><td>1, 2, 4</td><td>1,33</td><td>1,53</td><td>+</td><td></td></tr><tr><td>PFE0040c</td><td>mature parasite-infected erythrocyte surface antigen (MESA) or PfEMP2</td><td>yes</td><td>Hsp40, type IV</td><td>2, 4</td><td>0,92</td><td>1,87</td><td></td><td></td></tr><tr><td>PFA0660w</td><td>protein with DNAJ domain, dnj1/sis1 family</td><td>yes</td><td>Hsp40, type II</td><td>1, 2, 4</td><td>1,01</td><td>1,34</td><td>+</td><td></td></tr><tr><td>MAL7P1.7</td><td>RESA-like protein</td><td>yes</td><td>resa like</td><td>3</td><td></td><td>0,83</td><td></td><td></td></tr><tr><td>PFL0050c</td><td>hypothetical protein</td><td>yes</td><td>resa like</td><td>3</td><td></td><td>0,86</td><td></td><td></td></tr><tr><td>PF07_0029</td><td>heat shock protein 86</td><td>no</td><td>hsp90</td><td>1</td><td></td><td>0,9</td><td>+</td><td></td></tr><tr><td>PF07_0030</td><td>heat shock protein 86 family protein</td><td>no</td><td>Hsp86 family</td><td>5</td><td></td><td>1,63</td><td>+</td><td></td></tr><tr><td>PF08_0054</td><td>heat shock 70 kDa protein</td><td>no</td><td>hsp70</td><td>1</td><td></td><td>0,88</td><td>+</td><td></td></tr><tr><td>PF11_0034</td><td>hypothetical protein</td><td>yes</td><td>Hsp40, type IV</td><td>1, 2, 4</td><td></td><td>0,81</td><td>+</td><td></td></tr><tr><td>PFA0675w</td><td>P. falciparum RESA-like protein with DnaJ domain</td><td>yes</td><td>Hsp40, type IV</td><td>1, 2, 4</td><td></td><td>1,96</td><td>+</td><td></td></tr><tr><td>PF14_0700</td><td>hypothetical protein, conserved</td><td>no</td><td>Hsp40, type III</td><td>1, 2, 4</td><td></td><td>0,81</td><td></td><td></td></tr><tr><td>PF11_0509</td><td>ring-infected erythrocyte surface antigen putative</td><td>yes</td><td>Hsp40, type IV</td><td>1, 2, 4</td><td></td><td>-1,14</td><td></td><td></td></tr><tr><td>PF13_0102</td><td>DNAJ-like Sec63 homologue</td><td>no</td><td>Hsp40, type III</td><td>1, 2, 4</td><td></td><td>-1,04</td><td></td><td></td></tr><tr><td>PF11_0216</td><td>hypothetical protein</td><td>no</td><td>Heat shock factor binding protein I</td><td>6</td><td></td><td>-0,93</td><td>-</td><td></td></tr><tr><td>PFC0975c</td><td>PFCYP19, cyclophilin, peptidyl-prolyl cis-trans isomerase</td><td>no</td><td>cyclophilin</td><td>1</td><td></td><td>-0,97</td><td></td><td></td></tr></tbody></table>

Gene ID: identifier as found in PlasmoDB.

Art columns: log ratios of gene expression under artesunate exposure at 90 minutes and 3 hours.

HS column: over (+) or under (-) expression.

Results were compared with transcriptome modifications induced by heat-shock (HS) [12] and chloroquine (CQ) [8].

Chaperone- and chaperone-related genes were defined based on: 1) Acharya et al. [19], 2) Botha et al. [20], 3) resa-like annotation from GeneDB, 4) genes fitting the HMM profile of DNAJ (PF00226) from Pfam, 5) The Hsp86 family annotation from GeneDB, 6) genes fitting the HMM profile of heat shock factor binding protein I (PF06825) from Pfam.

homologue of *P. falciparum*, which modulates susceptibility to antimalarial drugs such as quinine, mefloquine and artemisinins [26-28]. An association has been shown in

the field between sensitivity to arylaminoalcohols or endoperoxides, *pfmdr1* allelic status and gene copy number [29], suggesting that an over-expressed functional

<table><caption>Table 3: Differentially expressed transport associated genes in artesunate-treated cells</caption><thead><tr><th>Gene ID</th><th>Description</th><th>ART 90 min</th><th>ART 3 h</th><th>HS</th><th>CQ</th></tr></thead><tbody><tr><td>PF07_0121</td><td>NMD3 protein, putative</td><td></td><td>0,85</td><td></td><td></td></tr><tr><td>PFB0435c</td><td>amino acid transporter, putative</td><td></td><td>0,89</td><td></td><td></td></tr><tr><td>PFE1455w</td><td>sugar transporter, putative</td><td></td><td>0,92</td><td></td><td>+</td></tr><tr><td>PFE1150w</td><td>multidrug resistance protein, Pfmdr1</td><td></td><td>1,10</td><td></td><td></td></tr><tr><td>PF11_0141</td><td>UDP-galactose transporter, putative</td><td>-0,86</td><td>-0,82</td><td>-</td><td></td></tr><tr><td>PFC0725c</td><td>formate-nitrate transporter, putative</td><td></td><td>-1,16</td><td>-</td><td></td></tr><tr><td>MAL8P1.13</td><td>integral membrane protein, conserved/folate-Biopterin transporter</td><td></td><td>-1,13</td><td></td><td></td></tr><tr><td>MAL13P1.206</td><td>Na+ -dependent Pi transporter, sodium-dependent phosphate transporter PfPiT</td><td></td><td>-0,90</td><td></td><td></td></tr><tr><td>PF13_0252</td><td>nucleoside transporter I PfNT1/PfENT1</td><td></td><td>-1,32</td><td></td><td></td></tr><tr><td>PFL0170w</td><td>Transporter, major facilitator superfamily</td><td></td><td>-0,86</td><td>-</td><td></td></tr><tr><td>PFE0410w</td><td>triose or hexose phosphate/phosphate translocator, putative</td><td></td><td>-0,88</td><td>-</td><td></td></tr><tr><td>PF13_0227</td><td>vacuolar ATP synthase subunit D, putative</td><td></td><td>-0,92</td><td></td><td></td></tr><tr><td>PF11_0412</td><td>Vacuolar ATP synthase subunit F, putative</td><td></td><td>-0,81</td><td>-</td><td></td></tr><tr><td>PF13_0130</td><td>vacuolar ATP synthase subunit G, putative</td><td></td><td>-0,82</td><td>-</td><td></td></tr><tr><td>PFE0965c</td><td>vacuolar ATP synthetase, subunit C putative</td><td></td><td>-0,85</td><td></td><td></td></tr></tbody></table>

Gene ID: identifier as found in PlasmoDB.

Art columns: log ratios of gene expression under artesunate exposure at 90 minutes and 3 hours.

HS and CQ columns: over (+) or under (-) expression.

Results were compared with transcriptome modifications induced by heat-shock (HS) [12] and chloroquine (CQ) [8].

Genes involved in transport were selected from the pathway maps of Ginsburg [16].
pfmdr1 conferred a multidrug resistance like phenotype [30,31]. Over-expression of *pfmdr1* we show to occur under artesunate pressure may have an impact on susceptibility to the drug(s) associated to artemisinins upon ACT administration. It also suggests that not only gene copy number but also gene transcription rates and mRNA stability should be monitored when exploring mechanisms of field parasite drug resistance.

## Lipid metabolism and the Apicoplast

Eight genes involved in lipid metabolism were differentially expressed (Table 4).

P. falciparum contains an unusually high number of Acyl CoA synthases (ACS) and binding proteins that might play a role in fatty acid salvage and transport from the host cell [32,33]. Amongst the 13 different ACS genes, a non subtelomeric acyl CoA synthetase gene PFB0685c has expanded into a family of 9 duplicated genes mainly located in the subtelomeric regions of the genome and transcribed during the erythrocytic stage of the parasite. At least two ACS proteins PfACS1 (PF14_0761) and PfACS3 (PFL2570w) were shown to interact with ankyrin through their lysine-rich C-terminal domain [34]. Similar domains are only present in the genes of this family, possibly targeting these enzymes to the erythrocyte cytoplasm and/or membrane, an interesting possibility as all the ACS proteins have a signal peptide but are devoid of the PEXEL/HT motif. This would contribute to the capacity of P. falciparum to activate fatty acid scavenged from the plasma or the erythrocyte membrane, a critical metabolic function for the survival of the parasite. All 4 ACS genes over-expressed under artesunate (PFB0685c, PFB695c, MAL13P1_485 and PF14_0761) belong to this family.

Two acyl-CoA binding protein genes, PF08_0099 and PF10_0016, were under-expressed. PF14_0664, an Acetyl CoA carboxylase which generates Malonyl-CoA, likely to be located in the apicoplast [32], was already found over-

expressed at 90 minutes. Interestingly, the triose phosphate transporter located in the outermost membrane of the apicoplast, PFE0410w, was under-expressed under artesunate. This transporter, also called PFoTPT, fuels the apicoplast by providing carbon, reducing power and ATP [35], but also quite importantly, dihydroxyacetone phosphate (DHAP), the precursor needed for isoprenoid and for phospholipids biosynthesis [36]. Altogether these data suggest an alteration of fatty acid metabolism both intracellularly and in the apicoplast.

## Mitochondrial genes and genes encoding proteins targeted to the mitochondrion

The main metabolic function allowed by the active mitochondrial electron transport chain maintained by *Plasmodium* is regeneration of ubiquinone, which is required as the electron acceptor for dihydroorotate dehydrogenase, an essential enzyme for pyrimidine biosynthesis [37]. Indeed, 3 important genes involved in mitochondrial electron transport are under-expressed in presence of artesunate: 2 subunits of the cytochrome c oxidase (downstream of the cytochrome bc1 complex inhibited by atovaquone) coxI and col and ubiquinol-cytochrome c reductase hinge protein, PF14_0248 (Table 5). This interference with the electron transport chain could in turn affect pyrimidine biosynthesis.

PFC0975c encoding a peptidyl-prolyl cis-trans isomerase cyclophilin-type belonging to the chaperone network of the mitochondrial matrix [16] is under-expressed in response to artesunate. Cyclophilin D/peptidylprolyl isomerase activity is a regulator of mitochondrial permeability transition, a non selective inner-membrane permeabilization occurring in response to increased calcium load and redox stress. These 2 conditions are induced by artesunate. Such a phenomenon can lead to necrosis through activation of phospholipases, proteases and nucleases [38].

<table><caption>Table 4: Differentially expressed genes involved in lipid metabolism in artesunate-treated cells</caption><thead><tr><th>Gene ID</th><th>Description</th><th>ART 90 min</th><th>ART 3 h</th><th>HS</th><th>CQ</th></tr></thead><tbody><tr><td>MAL13P1.485</td><td>acyl-coa ligase antigen</td><td></td><td>0,94</td><td></td><td></td></tr><tr><td>PFB0685c</td><td>acyl-CoA synthetase, PfACS9</td><td></td><td>0,83</td><td></td><td></td></tr><tr><td>PFB0695c</td><td>acyl-CoA synthetase</td><td>0,88</td><td>1,02</td><td>+</td><td></td></tr><tr><td>PF14_0761</td><td>fatty acyl CoA synthetase I PfACS1</td><td></td><td>0,99</td><td></td><td></td></tr><tr><td>PF14_0664</td><td>biotin carboxylase subunit of acetyl CoA carboxylase, putative</td><td>0,97</td><td>1,31</td><td>+</td><td></td></tr><tr><td>PF08_0099</td><td>acyl CoA binding protein, putative</td><td></td><td>-0,89</td><td></td><td></td></tr><tr><td>PF10_0016</td><td>acyl CoA binding protein, putative</td><td></td><td>-1,04</td><td>-</td><td></td></tr><tr><td>PFE0410w</td><td>triose or hexose phosphate/phosphate translocator, putative</td><td></td><td>-0,88</td><td>-</td><td></td></tr></tbody></table>

Gene ID: identifier as found in PlasmoDB.

Art columns: log ratios of gene expression under artesunate exposure at 90 minutes and 3 hours.

HS column: over (+) or under (-) expression.

Results were compared with transcriptome modifications induced by heat-shock (HS) [12] and chloroquine (CQ) [8].

Genes were selected based on [16].
## Purine/pyrimidine metabolism

In addition to the possibly impaired respiratory chain interfering with pyrimidine synthesis, expression of other genes related to purine/pyrimidine biosynthesis was also modified by artesunate pressure (Table 6): cytidine and deoxycytidylate deaminase (PF13_0259) hypoxanthine phophorybosyl transferase (PF10_0121), deoxyuridine 5'triphosphate nucleotidohydrolase (PF11_0282) and uridine phosphorylase (PFE0660c) are all under-expressed, while dihydroorotase (PF14_0697) is over-expressed. Interference with such key metabolic processes may readily cause parasite death.

## Signaling/kinases

The kinome of *P. falciparum* is atypical, with an atypical MAP kinase family and a particular R45-FIKK kinase family [39].

Fourteen genes encoding kinases, some of which potentially related to signal transduction, were differentially expressed in presence of artesunate (Table 7). Four of the 11 over-expressed genes were already differentially expressed at 90 minutes, among which the calcium calmodulin dependent protein kinase 2, PfPK2 (PFL1885c). PfPK2 belongs to the CamK group of *P. falciparum* kinases, a family of kinases which can be activated in response to increased calcium levels, possibly resulting from specific inhibition of the SERCA PfATPase6 by the drug [40]. The other 3 genes with increased expression levels at 90 minutes were members of the R45-FIKK multigene family [18,41]. This Apicomplexa-specific kinase gene family is composed of 20 genes all presenting sub-telomeric positions, scattered over 11 chromosomes. It has recently been shown that these kinases localize in different compartments of the infected erythrocyte, some being associated with the erythrocyte membrane [42]. Under artesunate pressure, 3 members of the R45-FIKK family were over-expressed at 90 minutes (PFI0095c, MAL7P1.144 and PFL0040c) and remained at elevated levels subsequently, while two additional R45-FIKK Kinases displayed increased transcript levels at 3 hours, namely PFD1165w

and PF11_0510. All 5 present export motifs. PFD1165w and PFL0040c are associated with Maurer's clefts, and PFL0040c with the infected erythrocyte membrane as well [42].

## Transcription associated proteins; zinc finger genes

Transcriptome studies of parasites during the erythrocytic stages of development have shown that in many cases, transcription is soon followed by translation, when the parasite "needs" the considered proteins to further develop [4]. In contrast, comparison between transcriptionic and proteomic studies have shown the existence of discrepancies in RNA and protein expression for about 50% of the genes expressed during the erythrocytic development [43,44]. This has been clearly confirmed for certain genes, such as *var* genes transcribed at the early ring stage and translated several hours later [45] or a series of sexual stage proteins the corresponding genes of which are transcribed long before translation occurs [46]. Genome mining revealed that relatively few genes (156, 2/3 less than expected) encode transcription associated proteins [17]. Six non-Zn finger transcription associated protein genes, as defined by Coulson [17], were over-expressed under artesunate, 3 were under-expressed (Table 8). This relatively low number suggests that transcriptional regulation may not play a major role in gene-expression modulation under drug-induced stress. Altered decay rates may modulate gene expression in drug-exposed parasites and contribute to the dynamic changes of mRNA levels observed here. In this regard, it is worth noting the recently outlined decreasing RNA decay rate as the parasite matures. A list of genes encoding putative decay components and possibly contributing to such a mechanism of post-transcriptional regulation in *P. falciparum* has been proposed [47]. However, there was no observed effect of artesunate on expression of these genes.

The low number of TAP identified by Coulson [17] contrasted with an over-representation of proteins with Zn finger motifs [48] potentially involved in regulating RNA stability. In the present study, 12 genes identified as puta-

<table><caption>Table 5: Differentially expressed genes related to the mitochondrion in artesunate treated cells</caption><thead><tr><th>Gene ID</th><th>Description</th><th>ART 90 min</th><th>ART 3 h</th><th>HS</th><th>CQ</th></tr></thead><tbody><tr><td>PF08_0054</td><td>heat shock 70 kDa protein</td><td></td><td>0,88</td><td>+</td><td></td></tr><tr><td>cox1</td><td>putative cytochrome oxidase I</td><td></td><td>-1,07</td><td></td><td></td></tr><tr><td>col</td><td>putative cytochrome oxidase I</td><td></td><td>-1,15</td><td></td><td></td></tr><tr><td>PF14_0248</td><td>ubiquinol-cytochrome c reductase hinge protein, putative</td><td></td><td>-0,99</td><td>-</td><td></td></tr><tr><td>PFC0975c</td><td>PFCYP19, cyclophilin, peptidyl-prolyl cis-trans isomerase</td><td></td><td>-0,97</td><td></td><td></td></tr></tbody></table>

Gene ID: identifier as found in PlasmoDB.

Art columns: log ratios of gene expression under artesunate exposure at 90 minutes and 3 hours.

HS column: over (+) or under (-) expression.

Results were compared with transcriptome modifications induced by heat-shock (HS) [12] and chloroquine (CQ) [8].

Genes were selected based on [16].
<table><caption>Table 6: Differentially expressed genes related to purine/pyrimidine metabolism in artesunate-treated cells</caption><thead><tr><th>Gene ID</th><th>Description</th><th>ART 90 min</th><th>ART 3 h</th><th>HS</th><th>CQ</th></tr></thead><tbody><tr><td>PF14_0697</td><td>dihydroorotase, putative</td><td>0,96</td><td>0,93</td><td></td><td></td></tr><tr><td>PF13_0259</td><td>cytidine and deoxycytidylate deaminase family, putative</td><td></td><td>-1,22</td><td></td><td></td></tr><tr><td>PF10_0121</td><td>hypoxanthine phosphoribosyltransferase</td><td></td><td>-1,24</td><td></td><td></td></tr><tr><td>PF11_0282</td><td>deoxyuridine 5'-triphosphate nucleotidohydrolase, putative</td><td>-0,88</td><td>-0,95</td><td></td><td></td></tr><tr><td>PFE0660c</td><td>purine nucleotide phosphorylase, putative</td><td>-1,07</td><td>-1,53</td><td></td><td></td></tr></tbody></table>

Gene ID: identifier as found in PlasmoDB.

Art columns: log ratios of gene expression under artesunate exposure at 90 minutes and 3 hours.

Results were compared with transcriptome modifications induced by heat-shock (HS) [12] and chloroquine (CQ) [8].

Genes were selected based on [16].

Zn finger bearing proteins were differentially expressed in the presence of artesunate; 10/12 were over-expressed, two of which (PFI0470w and PFE1245w) already at 90 minutes (Table 9).

Discussion

The methodology adopted here allowed us to find reproducible alterations of the transcriptome in *P. falciparum* parasites under artesunate pressure. To overcome the major hurdle of developmentally regulated genes, we combined the use of a fast acting drug at high dose with limited incubation time at 5 successive but synchronized developmental stages. Most observed alterations were of low amplitude, possibly reflecting particular mechanisms of gene regulation in malaria parasites. This resulted in a relatively broad perturbation, that nevertheless was highly significant and consistent across experiments. We did not detect altered expression of a single specific metabolic pathway or of the putative drug target molecules but several functionally important groups of genes exhibited

interesting dynamic changes. Within the intrinsic limitations of all transcriptome studies (in which steady-state RNA levels are studied and not protein quantities, post-translational modifications or biological activities), the overall picture that emerges from our current analysis of the transcriptome changes suggests i) perturbation of the intracellular trafficking and organisation, including changes in chaperones, transporters, remodeling the erythrocyte space beyond the parasitophorous vacuolar membrane up to the membrane of the infected erythrocyte; ii) probable down-phasing of metabolism as deduced from down-regulation of key enzymes/transporters implicated in purine, pyrimidine and isoprenoid synthesis; iii) an altered mitochondrion; iv) an altered redox homeostasis and possibly altered protein turnover. In addition, the transcriptome analysis highlighted numerous genes with unknown function, annotated as coding for hypothetical proteins [see Additional file 2]. Some display consistent high and early up-regulation and clearly deserve further investigation.

<table><caption>Table 7: Differentially expressed genes encoding kinases in artesunate-treated cells</caption><thead><tr><th>Gene ID</th><th>Description</th><th>Kinase type</th><th>Name</th><th>ART 90 min</th><th>ART 3 h</th><th>HS</th><th>CQ</th></tr></thead><tbody><tr><td>PFL0040c</td><td>protein kinase, R45/FIKK family</td><td>R45/FIKK</td><td>FIKK12</td><td>0,88</td><td>1,16</td><td></td><td></td></tr><tr><td>MAL7P1.144</td><td>protein kinase, R45/FIKK family</td><td>R45/FIKK</td><td>FIKK7.1</td><td>0,90</td><td>1,82</td><td></td><td></td></tr><tr><td>PF10095c</td><td>protein kinase, R45/FIKK family</td><td>R45/FIKK</td><td>FIKK9.1</td><td>0,83</td><td>1,28</td><td></td><td></td></tr><tr><td>PF11_0510</td><td>protein kinase, R45/FIKK family</td><td>R45/FIKK</td><td>FIKK11</td><td></td><td>0,95</td><td></td><td></td></tr><tr><td>PFD1165w</td><td>protein kinase, R45/FIKK family</td><td>R45/FIKK</td><td>FIKK4.1</td><td></td><td>1,05</td><td></td><td></td></tr><tr><td>PFL1885c</td><td>calcium/calmodulin-dependent protein kinase 2, putative</td><td>kinase</td><td>PfPK2</td><td>1,35</td><td>1,21</td><td></td><td></td></tr><tr><td>PFL2280w</td><td>cyclin g-associated kinase, putative</td><td>kinase</td><td></td><td></td><td>1,10</td><td></td><td></td></tr><tr><td>PF14_0346</td><td>cGMP-dependent protein kinase 1, beta isozyme, putative</td><td>kinase</td><td></td><td></td><td>1,03</td><td></td><td></td></tr><tr><td>PF11_0242</td><td>protein kinase</td><td>kinase</td><td></td><td></td><td>0,95</td><td></td><td></td></tr><tr><td>PF14_0294</td><td>mitogen-activated protein kinase 1, PfMAP1</td><td>kinase</td><td>PfMAP1</td><td></td><td>1,05</td><td></td><td></td></tr><tr><td>PFD0975w</td><td>ROI kinase-like protein</td><td>RIO1</td><td></td><td></td><td>0,90</td><td>+</td><td></td></tr><tr><td>PF11_0377</td><td>casein kinase 1, PfCK1</td><td>kinase</td><td></td><td></td><td>-0,91</td><td>+</td><td></td></tr><tr><td>PF11_0227</td><td>serine/threonine protein kinase, putative</td><td>kinase</td><td></td><td></td><td>-0,98</td><td></td><td></td></tr><tr><td>PFB0150c</td><td>protein kinase, putative</td><td>kinase</td><td></td><td></td><td>-1,05</td><td></td><td></td></tr></tbody></table>

Gene ID: identifier as found in PlasmoDB.

Art columns: log ratios of gene expression under artesunate exposure at 90 minutes and 3 hours.

CQ column: over (+) or under (-) expression.

Results were compared with transcriptome modifications induced by heat-shock (HS) [12] and chloroquine (CQ) [8].

Kinases were defined as genes fitting the HMM profile kinase PF00069 (update of Ward et al. [39]) and ROI Kinases as genes fitting the HMM profile PF01163. FIKK kinases were defined as described in Schneider and Puijalon [41].
<table><caption>Table 8: Differentially expressed genes encoding transcription-associated proteins in artesunate-treated cells</caption><thead><tr><th>Gene ID</th><th>Description</th><th>Pfam</th><th>annotation</th><th>ART 90 min</th><th>ART 3 h</th><th>HS</th><th>CQ</th></tr></thead><tbody><tr><td>PF11_0053</td><td>PfSNF2L</td><td>TAP clustering</td><td>Nucleosome remodeling::SNF2L</td><td></td><td>0,82</td><td></td><td></td></tr><tr><td>PF11_0264</td><td>DNA-dependent RNA polymerase</td><td>TAP clustering</td><td>MITOCHONDRIAL RNA POLYMERASE</td><td></td><td>1,02</td><td></td><td></td></tr><tr><td>PFL0560c</td><td>minichromosome maintenance protein, putative</td><td>TAP clustering</td><td>Minichromosome maintenance protein</td><td></td><td>0,85</td><td></td><td></td></tr><tr><td>PFE0090w</td><td>hypothetical protein, conserved</td><td>TAP clustering</td><td>Histone transcription regulator</td><td></td><td>0,82</td><td></td><td></td></tr><tr><td>PF11_0297</td><td>hypothetical protein</td><td>TAP clustering</td><td>CCR4-NOT complex::NOT2</td><td></td><td>0,86</td><td>+</td><td></td></tr><tr><td>PFE1245w</td><td>hypothetical protein, conserved</td><td>PF00642</td><td>Zinc finger C-x8-C-x5-C-x3-H type (and similar)(PF00642)</td><td>1,28</td><td>0,91</td><td>+</td><td>+</td></tr><tr><td>PF14_0236</td><td>hypothetical protein</td><td>PF00642</td><td>Zinc finger C-x8-C-x5-C-x3-H type (and similar)(PF00642)</td><td></td><td>0,89</td><td>+</td><td></td></tr><tr><td>PFC0680w</td><td>hypothetical protein, conserved</td><td>PF00642</td><td>Zinc finger C-x8-C-x5-C-x3-H type (and similar)(PF00642)</td><td></td><td>1,55</td><td>+</td><td></td></tr><tr><td>PF10_0186</td><td>hypothetical protein</td><td>PF00642</td><td>Zinc finger C-x8-C-x5-C-x3-H type (and similar)(PF00642)</td><td></td><td>0,98</td><td></td><td></td></tr><tr><td>MAL13P1.122</td><td>hypothetical protein, conserved</td><td>PF00856</td><td>SET domain (PF00856)</td><td></td><td>0,81</td><td></td><td></td></tr><tr><td>PF10470w</td><td>FHA domain protein, putative</td><td>PF00097</td><td>C3HC4 type (RING finger) (PF00097)</td><td>0,89</td><td>1,05</td><td></td><td></td></tr><tr><td>PF10_0046</td><td>hypothetical protein</td><td>PF00097</td><td>C3HC4 type (RING finger) (PF00097)</td><td></td><td>1,17</td><td></td><td></td></tr><tr><td>PFF0165c</td><td>hypothetical protein, conserved</td><td>PF00097</td><td>C3HC4 type (RING finger) (PF00097)</td><td></td><td>0,94</td><td></td><td></td></tr><tr><td>PFL0440c</td><td>hypothetical protein, conserved</td><td>PF00097</td><td>C3HC4 type (RING finger) (PF00097)</td><td></td><td>0,95</td><td></td><td></td></tr><tr><td>PF11_0315</td><td>hypothetical protein</td><td>TAP clustering</td><td>APICOPLAST RNA beta, beta, &amp; beta subunits</td><td></td><td>1,64</td><td></td><td></td></tr><tr><td>PFL0145c</td><td>high mobility group protein</td><td>PF00505</td><td>HMG (high mobility group) box (PF00505)</td><td></td><td>-0,97</td><td></td><td></td></tr><tr><td>MAL13P1.76</td><td>TFIIH basal transcription factor subunit</td><td>TAP clustering</td><td>TFIIH::p44</td><td></td><td>-0,83</td><td></td><td></td></tr><tr><td>PF14_0718</td><td>hypothetical protein, conserved</td><td>TAP clustering</td><td>RNA polymerase II-associated factor SOH1</td><td></td><td>-0,86</td><td></td><td></td></tr></tbody></table>

Gene ID: identifier as found in PlasmoDB.

Art columns: log ratios of gene expression under artesunate exposure at 90 minutes and 3 hours.

HS and CQ columns: over (+) expression.

Results were compared with transcriptome modifications induced by heat-shock (HS) [12] and chloroquine (CQ) [8].

Genes encoding transcription-associated proteins (TAPs) were selected from the TAP clustering of Coulson et al. [17] and completed with an update of genes fitting the following HMM profiles defined by Coulson and colleagues [16]: PF00097, PF00642, PF00096, PF00856, PF00249, PF00439, PF00505, PF00628, PF00098, PF00808, PF01753, PF00313, PF00569, PF00643, PF08523, PF02146, PF03366.

The physiological significance of the observed perturbations need to be assessed in future studies. Transcriptome analysis is the first of a multi-step process, to be followed by scrutiny of the role of the individual genes identified here. In particular, it will be essential to establish whether or not the genes in question are implicated in a rate-limiting process and to analyse the dynamics of the corresponding protein products, including their turnover, post-translational modifications and cellular localization. An important question is to clarify which transcriptome alterations are specific to the response to artemisinins and which are related to a parasite lethally injured, as this has numerous implications for drug development.

We will discuss the transcriptome data in the frame of existing knowledge on mechanisms of action of artemisinin derivatives and mechanisms resulting in reduced susceptibility/resistance.

We will next discuss possible common changes in steady state RNA levels after different stresses inflicted to the parasite such as exposure to 41°C or treatment with chloroquine, that can cause parasite death.

The mechanism of action of artemisinins is debated. Artemisinin are endoperoxide-containing sesquiterpene lactones. Fe++-dependent activation of the endoperoxide bridge is required for the drug to be active [49]. Cleavage of the endoperoxide moiety forms highly reactive oxyl radicals that rearrange to more stable carbon-centered radical intermediates. These in turn form covalent adducts with parasite products thought to be responsible for a pleiotropic effect of the drug. Such pleiotropicity may explain the variety of genes we show here to be affected in their expression. Recent evidence suggest more specific
<table><caption>Table 9: Differentially expressed genes encoding proteins with Zn-finger motifs in artesunate-treated cells</caption><thead><tr><th>Gene ID</th><th>Description</th><th>Zn-finger</th><th>ART 90 min</th><th>ART 3 h</th><th>HS</th><th>CQ</th></tr></thead><tbody><tr><td>PFE1245w</td><td>hypothetical protein, conserved</td><td>Zinc finger C-x8-C-x5-C-x3-H type (PF00642)</td><td rowspan="5">1,28</td><td>0,91</td><td rowspan="5">+</td><td rowspan="5">+</td></tr><tr><td>PF14_0236</td><td>hypothetical protein</td><td>Zinc finger C-x8-C-x5-C-x3-H type (PF00642)</td><td>0,89</td></tr><tr><td>PFC0680w</td><td>hypothetical protein, conserved</td><td>Zinc finger C-x8-C-x5-C-x3-H type (PF00642)</td><td>1,55</td></tr><tr><td>PF10_0186</td><td>hypothetical protein</td><td>Zinc finger C-x8-C-x5-C-x3-H type (PF00642)</td><td>0,98</td></tr><tr><td>MAL13P1.122</td><td>hypothetical protein, conserved</td><td>PHD-finger (IPR001965)</td><td>0,81</td></tr><tr><td>PF14_0197</td><td>hypothetical protein</td><td>DNL zinc finger (PF05180)</td><td rowspan="6">0,89</td><td>0,87</td><td></td><td></td></tr><tr><td>PF10470w</td><td>FHA domain protein, putative</td><td>C3HC4 type (RING finger) (PF00097)</td><td>1,05</td><td></td><td></td></tr><tr><td>PF10_0046</td><td>hypothetical protein</td><td>C3HC4 type (RING finger) (PF00097)</td><td>1,17</td><td></td><td></td></tr><tr><td>PFF0165c</td><td>hypothetical protein, conserved</td><td>C3HC4 type (RING finger) (PF00097)</td><td>0,94</td><td></td><td></td></tr><tr><td>PFL0440c</td><td>hypothetical protein, conserved</td><td>C3HC4 type (RING finger) (PF00097)</td><td>0,95</td><td></td><td></td></tr><tr><td>PFB0140w</td><td>hypothetical protein</td><td>DHHC zinc finger domain (PF01529)</td><td>-0,91</td><td></td><td></td></tr><tr><td>MAL13P1.76</td><td>TFIIH basal transcription factor subunit</td><td>TFIIH C1-like domain (PF07975)</td><td>-0,83</td><td></td><td></td><td></td></tr></tbody></table>

Gene ID: identifier as found in PlasmoDB.

Art columns: log ratios of gene expression under artesunate exposure at 90 minutes and 3 hours.

HS and CQ columns: over (+) expression.

Results were compared with transcriptome modifications induced by heat-shock (HS) [12] and chloroquine (CQ) [8].

Genes containing Zn-motifs were defined as genes fitting the following Zn-finger HMM profiles: B-box zinc finger (PF00643), C3HC zinc finger-like (PF07967), C3HC4 type (RING finger) (PF00097), CHY zinc finger (PF05495), CSL zinc finger (PF05207), CW-type Zinc Finger (PF07496), DHHC zinc finger domain (PF01529), DNL zinc finger (PF05180), FYVE zinc finger (PF01363), HIT zinc finger (PF04438), MIZ/SP-RING zinc finger (PF02891), MYND finger (PF01753), PHD-finger (PF00628), Putative zinc finger motif, C2HC5-type (PF06221), Sec23/Sec24 zinc finger (PF04810), SWIM zinc finger (PF04434), TFIIH C1-like domain (PF07975), Tim10/DDP family zinc finger (PF02953), U1 zinc finger (PF06220), Zinc finger C-x8-C-x5-C-x3-H type (and similar)(PF00642), Zinc finger found in FPG and IleRS (PF06827), Zinc finger, C2H2 type (PF00096), Zinc finger, ZZ type (PF00569), ZN-finger, Zn-finger in Ran binding protein and others (PF00641), Zn-finger in ubiquitin-hydrolases and other protein (PF02148), ZPR1 zinc-finger domain (PF03367), AN1-like Zinc finger (PF01428), C2H2 and C2HC zinc fingers (SSF57667), PHD-finger (IPR001965).

mechanisms, such as inhibition of specific targets, with specific inhibition of the *P. falciparum* SERCA-type Ca++ pump (PfATPase6) [50]. Importantly, field isolates with markedly reduced *in vitro* susceptibility to artemether presented a mutant SERCA-type PfATPase6 [1,51]. Artemisinin has been shown to inhibit the endoplasmic reticulum - located SERCA-type in *Toxoplasma gondii*, a related Apicomplexan parasite as well [52]. We did not observe altered steady state levels of *PfATPase6* mRNA at any time point investigated, indicating that over-expression of this gene is not part of the parasite response to artesunate.

Ultrastructural alterations of the morphology and their kinetics are debated as well. Early-stage alterations of the mitochondrion, the endothelial reticulum and the digestive vacuole were detected in some studies [53] but not others [54]. Some of the transcriptome alterations may account for a disorganisation of the digestive vacuole and the mitochondrion. We observed increased levels of *Pfmdr1* expression concomitant with decreased mRNA levels of four of the subunits of the H⁺ vacuolar ATPase, which is implicated in regulating calcium intracellular stores of acidic compartments. Downproduction of the H⁺ vacuolar ATPase and/or alteration of its subunit ratio is predicted to negatively impact on intracellular calcium homeostasis. This is most probably on the critical path, since disruption of calcium homeostasis is central to cellular death, necrosis endophagy and apoptosis. Our transcriptome analysis also shows evidence that artesunate interferes with expression of genes related to the mito-

chondrion, in most cases by down-regulating their steady state RNA levels. This calls for additional studies on the mitochondrial activity after artesunate exposure.

Additional factors are to be considered in light of the transcriptome data, in particular the intracellular partitioning of the drug. Several studies have demonstrated the selective uptake of artemisinin derivatives. It has been shown that artemisinin and dihydroartemisinin, two sesquiterpene endoperoxide drugs closely related to artesunate, are transported by the tubulovesicular network (TVN) formed within the erythrocyte cytoplasm for nutrient import to the parasite, and that these drugs disrupted the protein organization of the TVN [55]. Numerous genes implicated in intracellular space modelling, trafficking, chaperones and transport display altered transcription profiles in artesunate-treated parasites. The markedly bias for such genes was indeed the most salient observation of our transcriptome study. It is tempting to speculate that this results in remodelling the intracellular space, including remodelling of the infected red blood cell membrane and TVN. This could either be a response aimed at correcting drug-inflicted damage, preventing further intake of drug or reflect the dysfunction resulting from such damage. Protein trafficking is all the more important for parasite growth and survival that the parasite develops inside a parasitophorous vacuole within the erythrocyte cytoplasm. Interactions between the parasite and the erythrocyte are multiple, with an increasing number of parasite proteins recognized as being involved (reviewed in [56]).
Parasite-induced modifications of the red blood cell through interactions of parasite proteins exported to the host cell membrane or cytoskeleton play a major role in parasite survival and virulence, through induction of infected red cell cytoadhesive properties or decreased deformability. Altered transcript levels of genes involved in host cell remodelling and intracellular trafficking have been reported in parasites lethally damaged by exposure to febrile temperatures [12]. This suggests that perturbing intracellular trafficking/remodeling is either on the path to parasite death or an attempt to overcome a lethal injury. This also suggests that alterations of intracellular trafficking by artesunate may be one of the mechanisms through which the drug can be active on a wide range of parasite developmental stages.

Analysis of the proteome of artemether-treated *P. falciparum* parasites showed that 30 of the 101 proteins identified using a MALDI-TOF-MS-based analysis [57], displayed a greater than three-fold expression level compared to untreated control cultures. Two of the overexpressed proteins, namely hsp90 (PF07_0029), PF14_0425 (fructose-biphosphate aldolase) also presented an increased transcript level in our study. Direct comparison of both studies is precluded due to limitations such as different sensitivities of the transcriptome and proteome approaches, more comprehensive coverage of the expressed products in the transcriptome approach, different drug exposure protocols and/or transcription/translation uncoupling, which has been described for approx 50% of the proteins [43,44].

To put the alterations observed under artesunate in perspective with changes induced by unrelated forms of stress, we analyzed the data in the light of the transcriptome alterations reported after chloroquine [8] treatment or incubation at 41°C [12], both treatments reported as having led to parasite death. The transcriptome studies are not directly comparable, as they were performed under different experimental conditions (synchronous/asynchronous parasites, glass slide microarrays/Affymetrix chips) and corresponded to different natures of stress (rapid acting vs.slow acting). However, identification of a few genes with altered expression under such different experimental conditions provides interesting insight. Such joint analysis is illustrated in the Venn diagram (Figure 3 see also tables 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 for different functional groups). There were only 4 genes identified as differentially expressed under all 3 types of stress: PFB0095c encoding PfEMP3, PFE1245w, encoding a putative zinc finger protein, and PF14_0151 encoding a putative RNA-binding protein and MAL7P1.171 encoding a hypothetical protein. PfEMP3 and MAL7P1.171 have an export motif. PfEMP3 is exported to the cytoplasmic face of the erythrocyte membrane in the mature stages. It asso-

ciates with the erythrocyte cytoskeleton, which it destabilizes [58]. Both PfEMP3 and MAL7P1.171 have been shown to contribute to trafficking of PfEMP1 [59,60]. Their role in the response to lethal stresses deserves to be investigated. Characterization of the products encoded by the remaining 2 common genes certainly warrants particular attention, as these genes may play a crucial role in the general response to stress.

The transcriptome of chloroquine-treated and artesunate-treated parasites shared three additional altered gene profiles, namely PFE1455w encoding a putative sugar transporter, PFI0095c and PFL0040c both encoding putative kinases belonging to the R45-FIKK family. Interestingly again, the three predicted proteins have an export motif. However, since there were large inter-experiment variations in the transcriptome analysis of chloroquine-treated parasites, commonalities are probably largely underestimated.

Ninety four genes were differentially regulated both under heat shock and artesunate, 63 over-expressed and 31 under-expressed under both conditions, amongst which members of the R45-FIKK kinase family, chaperones, molecules involved in lipid metabolism and number of hypothetical proteins -including as discussed above

<figure><img src="image_7.png" /><figcaption>**Figure 3**
Venn diagram of genes differentially expressed under artesunate treatment chloroquine treatment, and heat-shock. Purple circle: 398 genes differentially expressed under artesunate (ART). Orange circle: 30 genes differentially expressed under chloroquine (CQ) [8]. Green circle: 334 genes differentially expressed under heat-shock (HS) [12]. ID of genes differentially expressed in ART and CQ, HS and CQ, ART, CQ and HS are shown.</figcaption></figure>
numerous genes coding for proteins with an export motif. Two of the 5 R45-FIKK kinase genes over-expressed under artesunate were also over-expressed under chloroquine pressure [8]. Of the 3 R45/FIKK kinase genes found over-expressed under heat shock, one was also differentially expressed under chloroquine. These data suggest that expression of genes belonging to the R45-FIKK kinase family reflects the response of the parasite to different types of environmental modifications and that the parasite response to different natures of stress may include common pathways. The same could be said about genes belonging to 2 other families. Eight of the 13 chaperone/chaperone-related encoding genes over-expressed under artesunate were also over-expressed under heat-shock. Of the 30 Hsp40 type *P. falciparum* chaperones, 6 were over-expressed under artesunate, 4 in common with heat-shock.

## Conclusion

The rapid parasiticidal activity of artesunate at high doses allowed us to study the effects of the drug on the transcriptome of synchronized parasite cultures over a short 3 hour time window. A positive bias was observed in favor of subtelomeric localization, polymorphism and presence of potential export sequences for the subset of upregulated genes. The corresponding genes, related to protein trafficking, kinases or membrane remodeling, often belong to multi-gene families [61]. Low amplitude transcript level alterations were observed, but the large number of affected genes likely induces substantial perturbations of the interface between the parasite and its host cellular environment. Whether these constitute parasite survival mechanisms or in contrast metabolic cytotoxic alterations remains to be established. The transcriptome analysis identified what may represent pathways leading to parasite death, such as inhibition of purine/pyrimidine metabolism, interference with the mitochondrial electron transport chain, with protein turn-over or with the integrity of the food vacuole and calcium stores.

The high proportion of over-expressed genes encoding proteins exported from the parasite highlight the importance of extra-parasitic compartments as fields for exploration in drug research which, to date, has mostly focused on the parasite within its plasma membrane rather than within its intra- and extra-erythrocytic environment.

## Methods

### Parasite cultures

The FCR3 strain (FUP/CB line) [62] of *P. falciparum* was used throughout the experiments. Parasites were cultured using the method described by Trager and Jensen with the modifications described in Ralph et al. [63]. Cultures were tightly synchronized by lysing maturing forms through treatment with 0.3 M Alanine at two successive cycles, starting 3 hours after rings appeared in the culture [63].

Artesunate treatment was performed the following cycle (see below).

### Drug treatment

Artesunate powder (Arenco, Belgium) was dissolved in PBS, and the 1 µg/µL aliquoted stock solution was kept at -20°C. To detect possible loss of activity during storage of stock solution, the IC₅₀ of the drug was tested periodically. The IC₅₀ remained stable at around 3.4 nM (2.2-4.6 nM, SD = 1.2). This also confirmed the artesunate sensitivity of the FCR3 strain.

Determination of minimal drug concentration irreversibly damaging 100% of the parasites: The lethal artesunate dose to be used throughout this investigation was determined using synchronized cultures incubated with increasing concentrations of artesunate (30-100-300-500-3000 ng/mL) for 3 hours, washed twice with RPMI and reincubated at 37°C for 24-30 hours, until next reinvasion. Reinvasion was assessed by examining Giemsa stained blood smears. The lethal dose was the minimal artesunate concentration needed to prevent 100% of reinvasion.

Parasite treatment for studies of the transcriptome: artesunate was added at five different time points during the cycle following parasite synchronization, staggered between 20 hours and 30 hours after appearance of the first rings in the culture (at approximately 20, 22, 26, 28 and 30 hours), while untreated cultures were cultured in parallel in the absence of drug. Parasites were extracted in Trizol after 90 minutes and 3 hours of incubation.

### Total RNA preparation

Parasite cultures were harvested by centrifugation and lysed in 5 pellet volumes of Trizol (Gibco) before freezing at -80°C. Total RNA was prepared from thawed samples following the manufacturer's instructions. RNA Quality was assessed with an Agilent 2100 Bioanalyser.

### Microarrays and hybridizations

The microarrays used here have been described elsewhere [63]. Briefly, glass slides were spotted with 7392 70-mer oligonucleotides originating from the Malaria Oligo Set (Qiagen-Operon) and custom design oligonucleotides, covering most of all *P. falciparum* genes. Of the 5542 identified *P. falciparum* genes, the array covers 4741 genes (85%). RNA labelling and hybridization were performed as described [63]. For the comparison of drug versus no drug at each time point, dye swaps with two technical replicates were performed to compensate dye effect and to assess technical reproducibility, leading to 4 hybridized slides per experiment.

In order to determine the genes differentially expressed between time 0 and 3 hours due to parasite development,
three series of hybridizations were performed with cDNAs prepared at times 20/23 hours, 23/26 hours and 25/28 hours from the control cultures left without artesunate. Microarray data were deposited in ArrayExpress [ArrayExpress:E-MEXP-1435].

### Statistical analysis

All the slides were analyzed pooled together using R [64] and Bioconductor [65] with the limma package [66,67]. After logarithm transformation of ratio of median of the intensities (without background subtraction) in the two channels, data were weighted to penalize spots of poor quality, flagged, or spots for which intensities in the two channels were close to background (median of intensity less than median of background plus 2 standard deviations for both Cy3 and Cy5 signals). In a second step, an intensity-dependent normalization was applied to each slide (loess), followed by quantile normalization [68] applied across the arrays to ensure that the log ratios had the same empirical distribution across arrays and across channels. After fitting a linear model to the expression data for each probe using the least square method, comparisons of interest were extracted with a contrast matrix. A moderated t-statistic and a log-odds of differential expression were computed for each contrast for each gene. Finally, the raw p-values were adjusted using the Holm procedure and a type I error rate of 0.05 was applied. When for a given gene significant log ratios were obtained for several oligos, the retained log ratio was one with the maximum absolute value.

### Reverse transcription and qPCR

RNA samples were treated with DNase I (Invitrogen) for 15 minutes at 37°C in presence of RNaseOUT™ inhibitor (Invitrogen). RNA was reverse transcribed with Supercript II RT (Invitrogen) using random hexamers with pre-extension step of 10 minutes at 25°C followed by extension for 50 minutes at 42°C. Then, cDNA was treated with RNase H for 20 minutes at 37°C.

Quantitative real time PCR was performed on cDNA using an ABI Prism 7900 HT thermal cycler system (Applied Biosystems) for 40 cycles (95°C for 15 s, 55°C for 15 s and 60°C for 45 s). Reactions were done in 20 µL volumes using SYBR Green PCR master Mix (Applied Biosystems) and 0.5 µM primers in triplicate with 6 concentrations of cDNA (2500 pg/µL, 625 pg/µL, 156 pg/µL, 39 pg/µL and 10 pg/µL) to assess primer amplification efficiency. Absence of DNA contamination was checked using RNA sample treated as cDNA without RT (2500 pg/µL).

Primers were designed using eprimer3 [69] [See Additional file 5 for the primer sequences]. Transcript abundance was compared using mean of ΔΔCt values calculated for all cDNA dilutions using PFI0425w (putative transporter) as endogenous normalizer and the con-

trol condition as reference. This gene was chosen as it remained non differentially expressed in all our experiments and showed very moderate levels of variation in whole-cycle transcriptome analysis [4].

### Abbreviations

* GSH: reduced glutathione; ACS: Acyl CoA synthetase; TAP: transcription associated protein; CQ: chloroquine; HS: heat-shock; ART: artesunate; SNP: single nucleotide polymorphism; HMM: hidden Markov model

### Authors' contributions

ON establishment of experimental protocols, parasite cultures, drug exposure, RNA extraction, q-PCR, genome mining, interpretation of results, manuscript preparation. EB establishment of experimental protocols, microarray spotting, q-PCR, genome mining, statistical analysis, interpretation of results, manuscript preparation. GD parasite cultures, drug exposure, RNA extractions, manuscript preparation. CP microarray spotting, hybridizations. MAD quality control of arrays, statistical analysis. OS microarray spotting, quality control of arrays. GG quality control of arrays, statistical analysis. SB interpretation of results, genome mining. JP interpretation of results. OMP interpretation of results, manuscript preparation. JYC spotting of microarrays, statistical analysis, interpretation of results. PHD establishment of experimental protocols, parasite cultures, RNA extraction, genome mining, interpretation of results, manuscript preparation. All authors have read and approved the final manuscript.

### Additional material

#### Additional File 1

Distribution of log ratios for statistically significant genes. Distribution of log ratios for genes differentially expressed upon 3 hour incubation with artesunate. Grey bars: genes selected (log ratio cut-off: +/- 0.8). Click here for file [http://www.biomedcentral.com/content/supplementary/1471-2164-9-388-S1.pdf]

#### Additional File 2

Mean log ratio of gene expression in presence/absence of artesunate at 90 minutes and 3 hours. List is made-up of genes displaying significant changes at 3 hour incubation time, across 5 experiments performed at time points staggered between 20–30 h of parasite development. Values shown express the mean log ratios calculated from all experiments; when for a given gene significant log ratios were obtained for several oligos, the retained log ratio was one with the maximum absolute value. In the 3 hour column, only values significant at 3 hours as defined by ANOVA are shown, with a mean log-ratio cut-off of +/- 0.8. Values in bold: log ratios of controls (0 hours vs 3 hours in absence of artesunate) have been subtracted from the log ratios obtained after 3 hours in artesunate. Red cell-background: over-expression. Green cell-background: under-expression. Gene ID and description: from PlasmoDB. Click here for file [http://www.biomedcentral.com/content/supplementary/1471-2164-9-388-S2.xls]
### Additional File 3

Confirmation by real-time quantitative rt-PCR of differential gene expression induced by 3 hour incubation with artesunate. Transcripts abundances were compared by ΔΔCt values calculated using PFI0425w (putative transporter) as endogenous control. Values for microarray and qPCR are expression fold changes between parasites in presence and in absence of artesunate (Pearson correlation coefficient between the two techniques: 0.81).

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-388-S3.xls]

### Additional File 4

Genes related to cell cycle regulation. Cell cycle associated genes were defined according to the GO biological process GO:0007049 (cell cycle). Red cell-background: over-expression. Green cell-background: under-expression.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-388-S4.xls]

### Additional File 5

Sequences of primers used for q RT-PCR validation.

Click here for file
[http://www.biomedcentral.com/content/supplementary/1471-2164-9-388-S5.xls]

## Acknowledgements

The work received financial support from the Délégation Générale pour l'Armement (DGA n°22120/DSP/SREAF et n°04 34 025), the Programme PAL+/Fonds National pour la Science, the Institut Pasteur, the Fonds dédié "Combattre les Maladies parasitaires" (Sanofi-Aventis/Ministère de l'Enseignement supérieur et de la Recherche) and the Programme Génopole. ON was supported by the Thailand Research Fund through the Royal Golden Jubilee PhD programme (Grant No. PHD/0157/2542) to JP and EB by the DGA. The authors are grateful to Dr Zbynek Bozdech for invaluable advice during the early stages of the study and to Dr Genevieve Milon for helpful comments and critical review of the manuscript.

## References

1. Jambou R, Legrand E, Niang M, Khim N, Lim P, Volney B, Ekala MT, Bouchier C, Esterre P, Fandeur T, Mercereau-Puijalon O: **Resistance of Plasmodium falciparum field isolates to in-vitro artemether and point mutations of the SERCA-type PfATPase6**. *Lancet* 2005, **366**(9501):1960-1963.

2. Zongo I, Dorsey G, Rouamba N, Tinto H, Dokomajilar C, Guiguemde RT, Rosenthal PJ, Ouedraogo JB: **Artemether-lumefantrine versus amodiaquine plus sulfadoxine-pyrimethamine for uncomplicated falciparum malaria in Burkina Faso: a randomized non-inferiority trial**. *Lancet* 2007, **369**(9560):491-498.

3. Meremikwu M, Alaribe A, Ejemot R, Oyo-lta A, Ekenjoku J, Nwachukwu C, Ordu D, Ezedinachi E: **Artemether-lumefantrine versus artesunate plus amodiaquine for treating uncomplicated childhood malaria in Nigeria: randomized controlled trial**. *Malar J* 2006, **5**(1):43.

4. Bozdech Z, Llinas M, Pulliam BL, Wong ED, Zhu J, DeRisi JL: **The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparum**. *PLoS Biol* 2003, **1**(1):e5.

5. Le Roch KG, Zhou Y, Blair PL, Grainger M, Moch JK, Haynes JD, De La Vega P, Holder AA, Batalov S, Carucci DJ, Winzeler EA: **Discovery of gene function by expression profiling of the malaria parasite life cycle**. *Science* 2003, **301**(5639):1503-1508.

6. Dahl EL, Shock JL, Shenai BR, Gut J, DeRisi JL, Rosenthal PJ: **Tetra-cyclines specifically target the apicoplast of the malaria parasite Plasmodium falciparum**. *Antimicrob Agents Chemother* 2006, **50**(9):3124-3131.

7. Deitsch K, Duraisingh M, Dzikowski R, Gunasekera A, Khan S, Le Roch K, Llinas M, Mair G, McGovern V, Roos D, Shock J, Sims J, Wiegand R, Winzeler E: **Mechanisms of gene regulation in Plasmodium**. *Am J Trop Med Hyg* 2007, **77**(2):201-208.

8. Gunasekera AM, Myrick A, Roch KL, Winzeler E, Wirth DF: **Plasmodium falciparum: Genome wide perturbations in transcript profiles among mixed stage cultures after chloroquine treatment**. *Exp Parasitol* 2007.

9. WHO: **Antimalarial drug combination therapy**. In WHO/CDS/RBM Geneva, World Health Organization; 2001.

10. Haynes RK, Krishna S: **Artemisinins: activities and actions**. *Microbes Infect* 2004, **6**(14):1339-1346.

11. Hien TT, Davis TM, Chuong LV, Illett KF, Sinh DX, Phu NH, Agus C, Chiswell GM, White NJ, Farrar J: **Comparative pharmacokinetics of intramuscular artesunate and artemether in patients with severe falciparum malaria**. *Antimicrob Agents Chemother* 2004, **48**(11):4234-4239.

12. Oakley MS, Kumar S, Anantharaman V, Zheng H, Mahajan B, Haynes JD, Moch JK, Fairhurst R, McCutchan TF, Aravind L: **Molecular factors and biochemical pathways induced by febrile temperature in intraerythrocytic Plasmodium falciparum parasites**. *Infect Immun* 2007, **75**(4):2012-2025.

13. Lingelbach K, Przyborski JM: **The long and winding road: protein trafficking mechanisms in the Plasmodium falciparum infected erythrocyte**. *Mol Biochem Parasitol* 2006, **147**(1):1-8.

14. Bahl A, Brunk B, Crabtree J, Fraunholz MJ, Gajria B, Grant GR, Ginsburg H, Gupta D, Kissinger JC, Labo P, Li L, Mailman MD, Milgram AJ, Pearson DS, Roos DS, Schug J, Stoeckert CJ Jr., Whetzel P: **PlasmoDB: the Plasmodium genome resource. A database integrating experimental and computational data**. *Nucleic Acids Res* 2003, **31**(1):212-215.

15. Hertz-Fowler C, Peacock CS, Wood V, Aslett M, Kerhornou A, Mooney P, Tivey A, Berriman M, Hall N, Rutherford K, Parkhill J, Ivens AC, Rajandream MA, Barrell B: **GeneDB: a resource for prokaryotic and eukaryotic organisms**. *Nucleic Acids Res* 2004, **32**(Data-base issue):D339-43.

16. Ginsburg H: **Malaria Parasite Metabolic Pathways**. [http://sites.huji.ac.il/malaria/].

17. Coulson RM, Hall N, Ouzounis CA: **Comparative genomics of transcriptional control in the human malaria parasite Plasmodium falciparum**. *Genome Res* 2004, **14**(8):1548-1554.

18. Doerig C: **Protein kinases as targets for anti-parasitic chemotherapy**. *Biochim Biophys Acta* 2004, **1697**(1-2):155-168.

19. Acharya P, Kumar R, Tatu U: **Chaperoning a cellular upheaval in malaria: heat shock proteins in Plasmodium falciparum**. *Mol Biochem Parasitol* 2007, **153**(2):85-94.

20. Botha M, Pesce ER, Blatch GL: **The Hsp40 proteins of Plasmodium falciparum and other apicomplexa: regulating chaperone power in the parasite and the host**. *Int J Biochem Cell Biol* 2007, **39**(10):1781-1803.

21. Pavithra SR, Kumar R, Tatu U: **Systems analysis of chaperone networks in the malarial parasite Plasmodium falciparum**. *PLoS Comput Biol* 2007, **3**(9):1701-1715.

22. Becker K, Tilley L, Vennerstrom JL, Roberts D, Rogerson S, Ginsburg H: **Oxidative stress in malaria parasite-infected erythrocytes: host-parasite interactions**. *Int J Parasitol* 2004, **34**(2):163-189.

23. Bozdech Z, Ginsburg H: **Antioxidant defense in Plasmodium falciparum--data mining of the transcriptome**. *Malar J* 2004, **3**:23.

24. Noiva R: **Protein disulfide isomerase: the multifunctional redox chaperone of the endoplasmic reticulum**. *Semin Cell Dev Biol* 1999, **10**(5):481-493.

25. Sztajer H, Gamain B, Aumann KD, Slomianny C, Becker K, Brigelius-Flohe R, Flohe L: **The putative glutathione peroxidase gene of Plasmodium falciparum codes for a thioredoxin peroxidase**. *J Biol Chem* 2001, **276**(10):7397-7403.

26. Reed MB, Saliba KJ, Caruana SR, Kirk K, Cowman AF: **Pgh1 modulates sensitivity and resistance to multiple antimalarials in Plasmodium falciparum**. *Nature* 2000, **403**(6772):906-909.

27. Sidhu AB, Uhlemann AC, Valderramos SG, Valderramos JC, Krishna S, Fidock DA: **Decreasing pfmdr1 copy number in Plasmodium falciparum malaria heightens susceptibility to mefloquine,**
lumefantrine, halofantrine, quinine, and artemisinin. *J Infect Dis* 2006, **194**(4):528-535.

28. Rohrbach P, Sanchez CP, Hayton K, Friedrich O, Patel J, Sidhu AB, Ferding MT, Fidock DA, Lanzer M: Genetic linkage of *pfmdr1* with food vacuolar solute import in *Plasmodium falciparum*. *Embo J* 2006, **25**(13):3000-3011.

29. Price RN, Cassar C, Brockman A, Duraisingh M, van Vugt M, White NJ, Nosten F, Krishna S: The *pfmdr1* gene is associated with a multidrug-resistant phenotype in *Plasmodium falciparum* from the western border of Thailand. *Antimicrob Agents Chemother* 1999, **43**(12):2943-2949.

30. Price RN, Uhlemann AC, Brockman A, McGready R, Ashley E, Phaipun L, Patel R, Laing K, Looareesuwan S, White NJ, Nosten F, Krishna S: Mefloquine resistance in *Plasmodium falciparum* and increased *pfmdr1* gene copy number. *Lancet* 2004, **364**(9432):438-447.

31. Duraisingh MT, Cowman AF: Contribution of the *pfmdr1* gene to antimalarial drug-resistance. *Acta Trop* 2005, **94**(3):181-190.

32. Mazumdar J, Sriepen B: Make it or take it: fatty acid metabolism of apicomplexan parasites. *Eukaryot Cell* 2007, **6**(10):1727-1735.

33. Bethke LL, Zilversmit M, Nielsen K, Daily J, Volkman SK, Ndiaye D, Lozovsky ER, Hart DL, Wirth DF: Duplication, gene conversion, and genetic diversity in the species-specific acyl-CoA synthetase gene family of *Plasmodium falciparum*. *Mol Biochem Parasitol* 2006, **150**(1):10-24.

34. Tellez M, Matesanz F, Alcina A: The C-terminal domain of the *Plasmodium falciparum* acyl-CoA synthetases PfACS1 and PfACS3 functions as ligand for ankyrin. *Mol Biochem Parasitol* 2003, **129**(2):191-198.

35. Mullin KA, Lim L, Ralph SA, Spurck TP, Handman E, McFadden GI: Membrane transporters in the relict plastid of malaria parasites. *Proc Natl Acad Sci U S A* 2006, **103**(25):9572-9577.

36. Ralph SA, van Dooren GG, Waller RF, Crawford MJ, Fraunholz MJ, Foth BJ, Tonkin CJ, Roos DS, McFadden GI: Tropical infectious diseases: metabolic maps and functions of the *Plasmodium falciparum* apicoplast. *Nat Rev Med Hyg* 1993, **4**(9):485-491.

37. Painter HJ, Morrisey JM, Mather MW, Vaidya AB: Specific role of mitochondrial electron transport in blood-stage *Plasmodium falciparum*. *Nature* 2007, **446**(7131):88-91.

38. Armstrong JS: The role of the mitochondrial permeability transition in cell death. *Mitochondrion* 2006, **6**(5):225-234.

39. Ward P, Equinet L, Packer J, Doerig C: Protein kinases of the human malaria parasite *Plasmodium falciparum*: the kinome of a divergent eukaryote. *BMC Genomics* 2004, **5**(1):79.

40. Krishna S, Eckstein-Ludwig U, Joet T, Uhlemann AC, Morin C, Webb R, Woodrow C, Kun JF, Kremsner PG: Transport processes in *Plasmodium falciparum*-infected erythrocytes: potential as new drug targets. *Int J Parasitol* 2002, **32**(13):1567-1573.

41. Schneider AG, Mercereau-Puijalon O: A new Apicomplexa-specific protein kinase family: multiple members in *Plasmodium falciparum*, all with an export signature. *BMC Genomics* 2005, **6**(1):30.

42. Nunes MC, Goldring JP, Doerig C, Scherf A: A novel protein kinase family in *Plasmodium falciparum* is differentially transcribed and secreted to various cellular compartments of the host cell. *Mol Microbiol* 2007, **63**(2):391-403.

43. Le Roch KG, Johnson JR, Florens L, Zhou Y, Santrosyan A, Grainger M, Yan SF, Williamson KC, Holder AA, Carucci DJ, Yates JR 3rd, Winzeler EA: Global analysis of transcript and protein levels across the *Plasmodium falciparum* life cycle. *Genome Res* 2004, **14**(11):2308-2318.

44. Hall N, Karras M, Raine JD, Carlton JM, Kooij TW, Berriman M, Florens L, Janssen CS, Pain A, Christophides GK, James K, Rutherford K, Harris B, Harris D, Churcher C, Quail MA, Ormond D, Doggett J, Trueman HE, Mendoza J, Bidwell SL, Rajandream MA, Carucci DJ, Yates JR 3rd, Kafatos FC, Janse CJ, Barrell B, Turner CM, Waters AP, Sinden RE: A comprehensive survey of the *Plasmodium* life cycle by genomic, transcriptomic, and proteomic analyses. *Science* 2005, **307**(5706):82-86.

45. Dzikowski R, Templeton TJ, Deitsch K: Variant antigen gene expression in malaria. *Cell Microbiol* 2006, **8**(9):1371-1381.

46. Dechering KJ, Thompson J, Dodemont HJ, Eling W, Konings RN: Developmentally regulated expression of *pfs16*, a marker for sexual differentiation of the human malaria parasite *Plasmodium falciparum*. *Mol Biochem Parasitol* 1997, **89**(2):235-244.

47. Shock JL, Fischer KF, Derisi JL: Whole genome analysis of mRNA decay in *Plasmodium falciparum* reveals a global lengthening of mRNA half-life during the intraerythrocytic development cycle. *Genome Biol* 2007, **8**(7):R134.

48. Gamsjaeger R, Liew CK, Loughlin FE, Crossley M, Mackay JP: Sticky fingers: zinc-fingers as protein-recognition motifs. *Trends Biochem Sci* 2007, **32**(2):63-70.

49. Stocks PA, Bray PG, Barton VE, Al-Helal M, Jones M, Araujo NC, Gibbons P, Ward SA, Hughes RH, Biagini GA, Davies J, Amewu R, Mercer AE, Ellis G, O'Neill PM: Evidence for a common non-heme chelatable-iron-dependent activation mechanism for semisynthetic and synthetic endoperoxide antimalarial drugs. *Angew Chem Int Ed Engl* 2007, **46**(33):6278-6283.

50. Eckstein-Ludwig U, Webb RJ, Van Goethem ID, East JM, Lee AG, Kimura M, O'Neill PM, Bray PG, Ward SA, Krishna S: Artemisinins target the SERCA of *Plasmodium falciparum*. *Nature* 2003, **424**(6951):957-961.

51. Legrand E, Volney B, Meynard JB, Mercereau-Puijalon O, Esterre P: In Vitro Monitoring of *Plasmodium falciparum* Drug Resistance in French Guiana: a Synopsis of Continuous Assessment from 1994 to 2005. *Antimicrob Agents Chemother* 2008, **52**(1):288-298.

52. Nagamune K, Beatty WL, Sibley LD: Artemisinin induces calcium-dependent protein secretion in the protozoan parasite *Toxoplasma gondii*. *Eukaryot Cell* 2007, **6**(11):2147-2156.

53. Maeno Y, Toyoshima T, Fujioka H, Ito Y, Meshnick SR, Benakis A, Milhous WK, Aikawa M: Morphologic effects of artemisinin in *Plasmodium falciparum*. *Am J Trop Med Hyg* 1993, **49**(4):485-491.

54. Del Pilar Crespo M, Avery TD, Hanssen E, Fox E, Robinson TV, Valente P, Taylor DK, Tilley L: Artemisinin and a series of novel endoperoxide antimalarials exert early effects on digestive vacuole morphology. *Antimicrob Agents Chemother* 2008, **52**(1):98-109.

55. Akompong T, VanWye J, Ghori N, Haldar K: Artemisinin and its derivatives are transported by a vacuolar-network of *Plasmodium falciparum* and their anti-malarial activities are additive with toxic sphingolipid analogues that block the network. *Mol Biochem Parasitol* 1999, **101**(1-2):71-79.

56. van Ooij C, Haldar K: Protein export from Plasmodium parasites. *Cell Microbiol* 2007, **9**(3):573-582.

57. Makanga M, Bray PG, Horrocks P, Ward SA: Towards a proteomic definition of CoArtem in *Plasmodium falciparum* malaria. *Proteomics* 2005, **5**(7):1849-1858.

58. Pei X, Guo X, Coppel R, Mohandas N, An X: *Plasmodium falciparum* erythrocyte membrane protein 3 (*PfEMP3*) destabilizes erythrocyte membrane skeleton. *J Biol Chem* 2007, **282**(37):26754-26758.

59. Knuepfer E, Rug M, Klonis N, Tilley L, Cowman AF: Trafficking determinants for *PfEMP3* export and assembly under the *Plasmodium falciparum*-infected red blood cell membrane. *Mol Microbiol* 2005, **58**(4):1039-1053.

60. Maier AG, Rug M, O'Neill MT, Brown M, Chakravorty S, Szestak T, Chesson J, Wu Y, Hughes K, Coppel RL, Newbold C, Beeson JG, Craig A, Crabb BS, Cowman AF: Exported proteins required for virulence and rigidity of *Plasmodium falciparum*-infected human erythrocytes. *Cell* 2008, **134**(1):48-61.

61. Anantharaman V, Iyer LM, Balaji S, Aravind L: Adhesion molecules and other secreted host-interaction determinants in Apicomplexa: insights from comparative genomics. *Int Rev Cytol* 2007, **262**:1-74.

62. Fandeur T, Bonnefoy S, Mercereau-Puijalon O: In vivo and in vitro derived Palo Alto lines of *Plasmodium falciparum* are genetically unrelated. *Mol Biochem Parasitol* 1991, **47**(2):167-178.

63. Ralph SA, Bischoff E, Mattei D, Sismeiro O, Dillies MA, Guigon G, Coppee JY, David PH, Scherf A: Transcriptome analysis of antigenic variation in *Plasmodium falciparum*--var silencing is not dependent on antisense RNA. *Genome Biol* 2005, **6**(11):R93.

64. R Development Core Team: R: A language and environment for statistical computing. [http://www.R-project.org].

65. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. *Genome Biol* 2004, **5**(10):R80.
66. Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. *Stat Appl Genet Mol Biol* 2004, **3**:Article3.

67. Smyth GK: **Limma: linear models for microarray data.** In *Bioinformatics and Computational Biology Solutions using R and Bioconductor*
Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W.
New York, Springer; 2005:397–420.

68. Yang YH, Thorne NP: **Normalization for two-color cDNA microarray data.** *Science and Statistics: A Festschrift for Terry Speed, IMS Lecture Notes - Monograph Series* 2003, **40**:403-418.

69. Steve R, Helen JS: **Primer3.** [http://primer3.sourceforge.net/].

Publish with **BioMed Central** and every
scientist can read your work free of charge

"BioMed Central will be the most significant development for disseminating the results of biomedical research in our lifetime."

Sir Paul Nurse, Cancer Research UK

Your research papers will be:

* available free of charge to the entire biomedical community

• peer reviewed and published immediately upon acceptance

* cited in PubMed and archived on PubMed Central

* yours — you keep the copyright

Submit your manuscript here:
http://www.biomedcentral.com/info/publishing_adv.asp

<figure><img src="image_8.png" /></figure>

BioMedcentral