Skip to main content

Global response of Plasmodium falciparum to hyperoxia: a combined transcriptomic and proteomic approach

Abstract

Background

Over its life cycle, the Plasmodium falciparum parasite is exposed to different environmental conditions, particularly to variations in O2 pressure. For example, the parasite circulates in human venous blood at 5% O2 pressure and in arterial blood, particularly in the lungs, at 13% O2 pressure. Moreover, the parasite is exposed to 21% O2 levels in the salivary glands of mosquitoes.

Methods

To study the metabolic adaptation of P. falciparum to different oxygen pressures during the intraerythrocytic cycle, a combined approach using transcriptomic and proteomic techniques was undertaken.

Results

Even though hyperoxia lengthens the parasitic cycle, significant transcriptional changes were detected in hyperoxic conditions in the late-ring stage. Using PS 6.0â„¢ software (Ariadne Genomics) for microarray analysis, this study demonstrate up-expression of genes involved in antioxidant systems and down-expression of genes involved in the digestive vacuole metabolism and the glycolysis in favour of mitochondrial respiration. Proteomic analysis revealed increased levels of heat shock proteins, and decreased levels of glycolytic enzymes. Some of this regulation reflected post-transcriptional modifications during the hyperoxia response.

Conclusions

These results seem to indicate that hyperoxia activates antioxidant defence systems in parasites to preserve the integrity of its cellular structures. Moreover, environmental constraints seem to induce an energetic metabolism adaptation of P. falciparum. This study provides a better understanding of the adaptive capabilities of P. falciparum to environmental changes and may lead to the development of novel therapeutic targets.

Background

Plasmodium falciparum is a protozoan parasite responsible for the most severe form of human malaria. This infection causes between 708,000 and 1,003,000 human deaths each year, most of them occurring in African children under the age of five years [1]. Several anti-malarial agents are used for malaria treatment and prophylaxis in endemic regions. However, the expansion of drug-resistance remains a serious problem. To develop new anti-malarial drugs, a better understanding of P. falciparum biology is required [2]. Some unique properties of the P. falciparum mitochondrion indicate that its respiratory metabolism could be exploited to generate chemotherapeutic targets [3]. Indeed, atovaquone, a mitochondrial cytochrome bc1 complex inhibitor, is currently used in combination with proguanil (Malaroneâ„¢, GlaxoSmithKline) for malaria treatment and prophylaxis [4, 5].

In vitro, P. falciparum growth is maximal in limited oxygen content (0.5%-5.0% O2), so the parasite is considered as a microaerophilic organism. In vivo, the parasite life cycle involves two hosts: the mosquito and the human, during which the parasite undergoes different cellular morphological changes and experiences oxygen pressure variations. The passage from one host to another implies metabolic adaptation and changes in the ultrastructural and physiological organization of mitochondria [6–8]. In humans, the parasite is exposed to varying oxygen pressures, which can reach up to 13% O2 in lung capillaries [9]. In mosquitoes, the parasite is exposed to 21% O2 levels in salivary glands. To adapt to these environmental constraints, the parasite has developed metabolic adaptations essential for it survival. During the intraerythrocytic cycle, two metabolic pathways are the major sources of ROS (superoxide anions, hydrogen peroxide and hydroxyl radicals) on P. falciparum - the mitochondrial respiratory chain and haemoglobin digestion [10]. In the presence of oxygen, P. falciparum can produce ATP by aerobic respiration and through glycolysis. Accordingly, microaerophilic metabolism may be a metabolic adaptation to prevent oxidative stress generation [10, 11]. The parasite also consumes haemoglobin in its digestive vacuole for protein biosynthesis. This metabolic pathway is a source of superoxide anions and ferriprotoporphyrin IX (FIX) accumulation [12], and thus ROS produced in the mitochondria could interact with these products of haemoglobin digestion and increase oxidative damage to the parasite cells. Therefore, P. falciparum has developed a preventive defence system to reduce cellular damage.

To identify metabolic pathways involved in the hyperoxia response, the effect of oxygen on P. falciparum was studied using high-throughput transcriptomic and proteomic analyses in the late-ring stage. These approaches were designed to minimize non-specific responses [13], and they revealed that a stress response occurs following parasite exposure to hyperoxia and that P. falciparum modifies the metabolism of two organelles (the mitochondrion and the digestive vacuole) as a metabolic adaptation to this environmental challenge.

Methods

Plasmodium falciparum in vitro culture

The 3D7 P. falciparum strain (the reference strain used for the genome sequencing project) was obtained from the Malaria Research and Reference Reagent source centre (MR4, managed by the American type culture collection). Parasitized human red blood cells (RBC type A+) were maintained in culture in RPMI 1640 medium (Invitrogen, Paisley, United Kingdom) supplemented with 10% human serum and buffered with 25 mM HEPES (Sigma-Aldrich, St Louis, MI, USA) and 25 mM NaCO3 (Sigma) and in an atmosphere of 5% O2, 5% CO2, and 90% N2 [14]. The haematocrit was maintained at 6% and the parasitaemia at 3-5%. Culture medium was changed every day. Strain clonality was verified every month using PCR genotyping of polymorphic genetic markers (msp1, msp2 and microsatellite loci) [15, 16].

Plasmodium falciparum culture synchronization

To obtain tightly synchronized parasite cultures, several synchronization steps were successively employed. First, parasitized erythrocytes were treated with D-sorbitol (ICN Biomedicals, Inc., CA, United States of America) as previously described [17]. This step enriched cultures in the ring parasite stage. Secondly, in the next parasitic cycle, schizonts were selected using CS columns on a VarioMACS unit (Miltenyi Biotec, Germany) according to standard procedures [18]. Five hours later, the culture was treated with 5% D-sorbitol to eliminate parasites in mature stages (schizonts). These successive parasite synchronization steps allowed us to obtain parasites tightly synchronized in the ring stage (the window from 0 to 5 h after parasite invasion of the erythrocyte). The synchronized parasites were maintained in standard culture conditions before treatment.

Hyperoxia exposure of P. falciparum cultures

Tightly synchronized cultures (ring stage aged between 4-9 hours) were split and subjected to two different conditions, either a normoxic atmosphere (5% O2, 5% CO2, 90% N2 gas mixture) or a hyperoxic atmosphere (21% O2, 5% CO2, 74% N2 gas mixture) at 37°C in two series II incubators (Model 3131, Forma Scientific, Inc.). Culture conditions were maintained over eight hours for transcript analysis (parasites aged between 12-17 hours) and for 12 hours for protein analysis (parasites aged between 16-21 hours). For each culture condition, four biological replicates were performed. Viability, parasitaemia and erythrocytic cycle stage proportions were monitored daily by examining blood smears stained with RAL® 555 (RAL, Martillac, France). Blood smears were taken at 0, 24, 32, 36, 48, and 78 hours.

RNA extraction

After incubation under normoxic or hyperoxic atmospheric conditions, total RNA from parasitized erythrocytes was extracted with TRIzol reagent (Invitrogen) according to the manufacturer's instructions. RNA extracts were treated with 1 U of RNase-free DNase I (Applied Biosystems, CA, United States of America) and quantified using a NanoDrop ND-1000 (Nanodrop Technologies, Wilmington, United States of America). The integrity of the RNA was controlled with an RNA nano chip (2100 Bioanalyzer, Agilent Biotechnologies, Wilmington, DE). Samples were immediately used or stored at - 80°C.

Microarray experiments and analysis

A 2X11 k custom P. falciparum whole genome microarray was designed and manufactured using SurePrint Inkjet technology® (Agilent Technologies). In brief, the microarray was composed of 10,128 sixty-mer oligonucleotides representing 5,364 coding sequences located in the chromosomal, apicoplastic and mitochondrial genomes. Additional probes were added to control for quality. Labelling and hybridization was carried out following the manufacturer's protocol (Two-Color Microarray-Based Gene Expression Analysis, Agilent Technologies). Briefly, starting with 300 ng of total RNA, fluorescent cRNA (antisense) was generated using the Low RNA Input Fluorescent Amplification Kit (Agilent Technologies) and either cyanine 3-labeled CTP (Cy-3) or cyanine 5-labeled CTP (Cy-5) fluorescent dyes (PerkinElmer Life Sciences, Boston, MA). Dye swap hybridization was performed for 17 h at 60°C using the In situ Hybridization Kit Plus (Agilent Technologies). All processing steps were performed in an ozone-controlled environment ([O3] < 2 ppb) to avoid ozone-induced degradation of cyanine dyes on microarray slides. Slides were scanned at 5 μm resolution with a G2505B DNA microarray scanner (Agilent Technologies). Image analysis and intra-array signal correction was performed using the Agilent Feature Extractor Software A.9.1.3. Data processing, analysis and visualization were performed using the Resolver software 7.1 (Rosetta Inpharmatics). An error model-based transformation pipeline was used to map replicate reporters to genes, perform inter-array normalization and calculate fold changes (FC) as described elsewhere [19]. Using these FC values, a gene set enrichment analysis (GSEA) was performed using the Mann-Whitney-U-test enrichment algorithm in the PathwayStudio software 6.0 (PS 6.0™, Ariadne Genomics). The GSEA procedure determines whether the behaviour of an a priori set of genes shows significant concordance across two different biological states. This GSEA analysis focuses on groups of genes that share common biological function in revealing differential levels of each transcript. A gene network was generated based on information extracted from the literature using Medscan™ and the P. falciparum-specific database PS 6.0™ and using the "physical or regulatory connections" parameter between genes. Fold change values of microarray data were imported into PS 6.0™ and used to interpret the pathway with gene regulation networks.

Real-time quantitative RT-PCR

cDNA was synthesized from total RNA (DNA-free) with random hexamers using the High-Capacity cDNA Archive Kit (Applied Biosystems). Primers with a melting temperature between 55 to 60°C were designed to yield a 94- to 146-bp product preferentially containing an exon/intron boundary. Specific primer sequences used for the qRT-PCR are summarized in Additional data (Additional file 1). Real-time PCR was performed using the 7900 HT Fast Real-Time PCR System (Applied Biosystems) in a 25-μL reaction volume with the Power SYBR Green® PCR Master Mix Kit (Applied Biosystems). Each sample was assayed in triplicate and analysed with the ABI PRISM Sequence Detection System software Version SDS 2.2.1 (Applied Biosystems). Amplification of the 18S rRNA sequence served as the internal control for normalization. At the end of each reaction, cycle threshold (Ct) was manually set to the level that reflected the best kinetic PCR parameters, and melting curves were acquired and analysed. Relative quantification analysis was performed using the 2-ΔΔCt method where ΔΔCt = (Ct target - Ct 18S rRNA)H- (Ct target - Ct 18S rRNA )N., and the data reflect changes in target gene expression between two experimental conditions (N: normoxic and H: hyperoxic groups) [20].

Protein sample preparation and CyDye labelling

After incubation under normoxic or hyperoxic atmospheric conditions, parasitized erythrocytes were washed three times in PBS medium (Invitrogen) and lysed in cold H2O-saponin (0.1%, Sigma) for 10 min. The lysate was then centrifuged at 1500 g for 5 min. The supernatant was discarded and the pellet containing free parasites was recovered by washing in PBS medium followed by a centrifugation step (1500 g for 5 min). The free parasites were washed until the supernatant became colourless. The pellet was then suspended in 4% (w/v) CHAPS (Sigma) and disrupted by ultrasonication (Vibracell 72412, Bioblock Scientific, Illkirch, France) five times for 60 seconds on ice at maximum amplitude. The lysate was then centrifuged at 16000 g for 15 min. The supernatant was further precipitated with acetone 100% (Sigma). The protein concentration for each sample was estimated using the BioRad Lowry-based DC assay (Biorad, Hercules, CA, USA) according to the manufacturer's instructions. Total proteins were suspended in standard cell lysis buffer (7 M urea, 2 M thiourea, 4% CHAPS, 30 mM Tris base, pH 8.5 (Sigma)) to obtain a protein concentration adjusted to 2.5 μg/μL. Protein samples were minimally labelled with CyDye according to the manufacturer's protocols (GE Healthcare, Piscataway, NJ) [21]. The mixture of labelled proteins was then separated by two-dimensional (2D) electrophoresis.

Two-dimensional electrophoresis, image analysis and in-gel digestion

Isoelectric focusing (IEF) was performed on 18-cm 3-10 linear IPG strips (GE Healthcare). Destreak buffer containing 1% (v/v) IPG buffer 3-10 was used for overnight rehydration of IPG strips. The samples were applied at the acidic end of the IPG strip using a cup-loading technique. IEF was carried out on a Ettan IPGphor II (GE Healthcare) electrophoresis unit at 20°C for a total of 45 kVh (ramp to 300 V in 3 hrs, ramp to 1000 V in 6 hrs, ramp to 8000 V in 3 hrs, hold at 8000 V for 4 hours). IPG strips were equilibrated in equilibration buffer containing 50 mM Tris-HCl, pH 8.6, 6 M urea, and 2% SDS and 30% glycerol supplemented with 1% (w/v) DTT for 15 min at room temperature, followed by protein alkylation (carbamidomethylation) in the same equilibration buffer containing 2.5% (w/v) iodoacetamide instead of DTT for 15 min at room temperature. IPG strips were then placed on the top of 10% uniform polyacrylamide gels. Strips were overlaid with 0.5% agarose in 1× running buffer containing bromophenol blue and the proteins were further separated by SDS-PAGE (10 W per gel) at 20°C in an Ettan DALT Six electrophoresis system (GE Healthcare). After electrophoresis, the gels with Cydye-labelled proteins were directly imaged using a Typhoon™ Trio Image scanner (GE Healthcare UK). The intensity was adjusted to ensure that the maximum volume of each image was within 60,000 - 80,000 U. Analysis of 2-D DIGE was performed with DeCyder 6.5 software (GE Healthcare) using the differential in-gel analysis (DIA) and the biological variation analysis (BVA) modules. Protein spots that were expressed differentially between two experimental conditions (|ratio|≥1.5, p ≤ 0.05 t-Test) were marked with master gel numbers. Based on DeCyder v6.5 analysis, spots of interest from gels stained with Imperial Blue Stain (Pierce) were excised and digested using a Shimadzu Xcise automated gel processing platform (Shimadzu Biotech, Kyoto, Japan) as described previously [22] and stored at -20°C.

Mass spectrometry analysis

The samples were analysed by nanoscale capillary liquid chromatography-tandem mass spectrometry (nano LC-MS/MS). Purification and analysis were performed on a C18 capillary column using a CapLC system (Waters, Milford, MA) coupled to a hybrid quadrupole orthogonal acceleration time-of-flight tandem mass spectrometer (Q-TOF Ultima, Waters, MA). Chromatographic separation was conducted on a reversed-phased capillary column (Atlantis™ dC18, 3 μm, 75 μm × 150 mm Nano Ease™, Waters, MA) with a 180-200 nl min-1 flow. The gradient profile consisted of a linear gradient from 95% A (H2O, 0.1% HCOOH) to 60% B (80% ACN, 0.1% HCOOH) in 60 min followed by a linear gradient to 95% B in 10 min. Mass data acquisitions were piloted by MassLynx 4.0 software using automatic switching between MS and MS/MS modes. The internal parameters of Q-TOF were set as follows. The electro-spray capillary voltage was set to 3.2 kV, the cone voltage was set to 30 V, and the source temperature was set to 80°C. The MS survey scan was m/z 400-1300 with a scan time of 1 s and an interscan time of 0.1 s. When the intensity of a peak rose above a threshold of 15 counts, tandem mass spectra were acquired. Normalized collision energies for peptide fragmentation were set using the charge-state recognition files for +2 and +3 peptide ions. The scan range for MS/MS acquisition was from m/z 50 to 1500 with a scan time of 1 s and an interscan time of 0.1 s. Fragmentation was performed using argon as the collision gas and with the collision energy profile optimized for various mass ranges and charges of precursor ions. Mass data collected during a nano LC-MS/MS analysis were processed using ProteinLynx Global Server 2.2 software (Waters) with the following parameters: no background subtraction, smooth 3/2 Savitzky Golay and no deisotoping to generate peak lists in the micromass pkl format. Pkl files were then fed into the local search engine Mascot Daemon v2.2.2 (Matrix Science, London, UK). The data were searched against the Homo sapiens (218356 sequences) and P. falciparum (13110 sequences) National Center for Biotechnology Information non-redundant (NCBInr) protein database (March, 2010). Search parameters allowed for one missed tryptic cleavage site, the carbamidomethylation of cysteine, and the possible oxidation of methionine; precursor and product ion mass error tolerance was < 0.2 Da. All identified proteins had a Mascot score greater than 29 and 38 respectively for P. falciparum and Homo sapiens, corresponding to statistically significant (p ≤ 0.05 t-Test) identification. Identifications were considered valid when they contained at least two peptide sequences per protein. If a single peptide sequence was identified per one protein, the mascot score and sequence coverage were taken into account (Additional file 2).

Results

Effect of hyperoxia on P. falciparum and experimental design

To determine the effects of hyperoxia on asexual blood stage P. falciparum parasites, tightly synchronized cultures of the 3D7 strain were exposed to normoxic (5% O2) and hyperoxic (21% O2) conditions for two life cycles. The proportion of parasitaemia and erythrocytic cycle stages were monitored by blood smears in triplicate at different times: 0, 24, 32, 36, 48, and 78 hours (Figure 1). In normoxic condition, the 3D7 P. falciparum strain had a life cycle of 45 hours with entry into schizogony at approximately 32 hours. Parasite exposure to 21% O2 increased the length of parasitic cycle and decreased the parasitaemia, but it did not alter the parasites morphology (Figure 1). At 48 hours, the second parasitic cycle began for parasites exposed to 5% O2 (nearly all of them were at ring stage (100%)), while the majority of parasites exposed to 21% O2 remained in the schizont stage (90% schizonts and 10% ring). Thus, hyperoxia induces a delay of P. falciparum cell cycle of four hours as previously described [23]. After reinvasion during the following cycle (third cycle), hyperoxia exposure did not change parasitaemia and parasites had a normal life cycle without excess lethality (unpublished data). These results seem to indicate a biological adaptation of the parasite to hyperoxia. To study the effects of hyperoxia on P. falciparum, complementary high-throughput transcriptomic and proteomic approaches were used. Transcriptome and proteome profiles from parasitized RBCs exposed to normoxic (5% O2) or hyperoxic (21% O2) atmospheric conditions were compared. The results were controlled by the experimental design on two levels: (i) the percentage of atmospheric oxygen and (ii) the timing of the parasitic stage. First, hyperoxic exposure was chosen at 21% O2 in light of the above data [23]. Second, transcriptome and proteome experiments were performed at the late-ring stage, after RBC reinvasion with a synchronization window of four hours. This stage allowed us to avoid the effects of cycle delay.

Figure 1
figure 1

Lengthening of P. falciparum cell cycle following hyperoxia exposure. Phenotypic changes in P. falciparum asexual blood stages were observed during two cycles under normoxic (A) and hyperoxic (B) conditions. The parasitaemia and parasitic stages were evaluated by blood smears at different times: 0, 24, 32, 36, 48, and 78 hours. The different intraerythrocytic stages and their corresponding percentages are indicated as follows: ring (R), trophozoite (T) and schizont (S).

Plasmodium falciparum response to hyperoxia treatment: microarray analysis

To investigate the response of P. falciparum to hyperoxia, the 3D7 strain was cultured in vitro under normoxic (5% O2) and hyperoxic (21% O2) atmospheric conditions. Three biological replicates in each group were performed and comparisons were made with a dye-swap experimental scheme. The raw microarray data are available in the Gene Expression Omnibus database [platform GPL9482 and samples from GSM466802 to GSM466807, 2010 [24]]. Among 5,364 coding sequences represented on the microarray, 219 genes were significantly altered following hyperoxia exposure (p ≤ 0.01, Student's t-Test, |FC|≤ 1.5, Additional file 3), among which 114 were up-expressed and 105 were down-expressed. Based on the selected 219 genes, a GSEA (FDR ≤ 0.05) was performed [25], and this allowed us to define 9 functional groups that were significantly altered (following hyperoxia treatment) (Additional file 4 and Table 1). Among these functions, "DNA repair," "Vacuolar acidification" and "Response to oxidative stress" were previously reported to be involved in the hyperoxia response [10, 26–30].

Table 1 Biological functions perturbed following hyperoxia exposure on P. falciparum

GSEA data were first integrated to create a gene network based on information extracted from the literature using Medscanâ„¢ and PS 6.0â„¢ (Ariadne Genomics). Next, the expression levels of genes included in the gene network were assigned using PS 6.0â„¢ and microarray expression data. This programme gives a dynamic view of metabolism during the hyperoxia response (Figures 2A and 2B). Thus, PS 6.0 analysis clustered 28 modulated-expression genes in five metabolic groups labelled "Energetic metabolism," "Protein folding," "Signal Transduction," "DNA repair" and "Translation" (Figure 2A). This transcriptomic analysis allowed us to identify up-regulated genes involved in DNA repair and protein folding and down-regulated genes linked to PKA-dependent signal transduction and glycolysis (Figure 2A). Additionally, PS 6.0 software revealed an alteration of an ATP-dependent sub-network: specifically, up-regulation of the mitochondrial ATP synthase complex and down-regulation of the V-type ATPase complex (Figure 2B).

Figure 2
figure 2

Gene networks involved in the metabolic adaptation to hyperoxia in tightly synchronized P. falciparum cultures. Using PS 6.0â„¢ software and microarray expression data, gene networks were built and genes modulated in response to hyperoxia are represented. (A) Representation of metabolic interrelations related to adaptative hyperoxia exposure. Energetic metabolism: aspartate carbamoyltransferase (P13_0240), carbamoyl-phosphate synthetase (PF13_0044), glutamine-fructose-6-phosphate transaminase (PF10_0245), glucose-6-phosphate dehydrogenase (PF14_0511), acetyl-CoA synthetase (PFF1350c), gapdh (PF14_0598) - Signal transduction: regulatory sub-unit of cAMP-dependent protein kinase (PFL1110c), catalytic sub-unit of cAMP-dependent protein kinase (PFI1685w), ser thr protein kinase (PF13_0085) - Translation: citrate synthase (PF10_0218), translation elongation factor 1 alpha 1 (PF13_0305), proteasome 26S subunit (MAL13P1.343), dihydrolipoamide dehydrogenase (PFL1550w), translation elongation factor 2 (PF14_0486), adaptor-related protein complex 1 (PF13_0062), polymerase RNA I (PF11_0358) - DNA repair: DNA primase (PF14_0366 and PFI0530c), rpa1 (PFI0235w) - Protein folding: ferredoxin reductase (PF11_0407), Hsp10 (PFL0740c), Hsp60 (PF10_0153), Hsp70 (PF11_0351 and PF08_0054), Hsp90 (PFL1070c and PF07_0029), prohibitin (PF08_0006), DnaJ (PFF1415c). (B) Representation of ATP-dependent gene sub-networks altered in hyperoxic conditions. V-type ATPase: V-type ATPase putative (MAL13P1.271), vacuolar ATP synthase subunit h putative (PF13_0034), vacuolar ATP synthetase putative (PFE0965c), vacuolar ATP synthase subunit F putative (PF11_0412), vacuolar ATP synthase subunit D putative (PF13_0227), vacuolar ATP synthase catalytic subunit a (PF13_0065) - Mitochondrial ATP synthase F1: ATP synthase subunit putative (PF14_0615), mitochondrial ATP synthase F1 epsilon subunit (MAL7P1.75), mitochondrial ATP synthase F1 alpha subunit putative (PFB0795w). Red and blue colors correspond respectively to up- and down-regulated genes compared between hyperoxic to normoxic conditions.

To confirm the GSEA data, eight genes presenting significant variations in expression were quantified using real-time qRT-PCR (Additional file 1). These genes, which were involved in glycolysis [PlasmoDB: PF14_0598, PF10_0245], antioxidant metabolism [PlasmoDB: PF14_0187, PF11_0087], signal transduction [PlasmoDB: PFL1110c, PFI1685w] and ATP synthase activity [PlasmoDB: PF10_0059, MAL7P1.13] were chosen according to their essential functions in response to hyperoxia. qRT-PCR was performed on total RNA extracted from the same three samples used for microarray analysis. Correlation coefficients for each specific-gene standard curve were always > 0.99 (unpublished data). Gene amplification was normalized by 18S rRNA [PlasmoDB: MAL7_18Sa] levels as previously described [31]. The two analyses (qRT-PCR and microarray) yielded consistent results for all the genes evaluated (Figure 3).

Figure 3
figure 3

Validation of microarray data by qRT-PCR. Data from qRT-PCR and a microarray of eight selected genes were compared between hyperoxic and normoxic conditions. Adjacent bars correspond to the mean log2 (fold change) and respective standard deviation and present qRT-PCR and microarray results in gray scale for the respective gene. The abbreviations and their corresponding gene IDs (PlasmoDB accession numbers) are indicated below the graphic. Crude values of mean fold changes are presented in the table. The same samples were used for the qRT-PCR and microarray experiments. gapdh: glyceraldehyde-3-phosphate dehydrogenase - gf6p: glutamine-fructose-6-phosphate transaminase - gst: glutathione S-transferase - pka: protein kinase A - atp6: mitochondrial ATP synthase F0 a subunit - atp4: mitochondrial ATP synthase F0 b subunit.

Plasmodium falciparum response to hyperoxia treatment: proteomic analysis

To identify P. falciparum proteins involved in the hyperoxia response, 2D-DIGE experiments coupled to MS were performed. Four independent cultures of P. falciparum cultivated under normoxic and hyperoxic conditions were included in this analysis. After protein separation on 2-DE, each gel was individually imaged and all gel images were analysed using the DeCyder 6.5 software. Among 1840 protein spots matched, 33 spots were significantly modulated (|FC|≥1.5, p ≤ 0.05 t-Test) following hyperoxia treatment (14 and 19 spots were up- and down-modulated, respectively; Figure 4). All spots were successfully identified by MS and corresponded to 14 Homo sapiens and 19 P. falciparum proteins (Table 2 and Additional file 5). However, some proteins were detected in more than one spot, indicating different isoforms. Indeed, only six proteins were identified for Homo sapiens and 13 for P. falciparum (Table 2). These results indicated that hyperoxia induced protein modulations at two levels: namely, protein expression and post-translational modification. Among the four spots detected (spot numbers 1301, 1314, 1326, and 1331) as Pf-Hsp70 protein [PlasmoDB: PF08_0054, GenBank: gi|124512406], only two isoforms were significantly up-regulated in hyperoxic conditions (Table 2 and Figure 5). These results indicated that the hyperoxia response could induce post-translational regulation of several parasite proteins.

Figure 4
figure 4

Alterations of the P. falciparum proteome under hyperoxic exposure. The proteins from P. falciparum parasites cultivated under normoxic (A) or hyperoxic (B) conditions were labelled with Cy3 and Cy5, respectively, and separated by 2-DE using a 10% homogeneous SDS polyacrylamide gel with a pH range from 3 to 10. As determined by DeCyder software, protein spots that were down- (A) or up- (B) regulated following hyperoxic exposure (|FC|≥1.5, p ≤ 0.05) are marked with master numbers (Table 2 and Additional file 5). Bold and italic numbers correspond, respectively, to proteins identified from P. falciparum and Homo sapiens. Areas of gels containing a high density of spots down- and up-modulated are enlarged.

Table 2 Proteins identified by differential 2-D DIGE analysis following hyperoxia exposure of P. falciparum
Figure 5
figure 5

Selective post-translational modification of Pf-Hsp70 following hyperoxia exposure. (A) Enlarged 2D-DIGE gel images and their corresponding three-dimensional profiles are shown for a series of four protein spots identified as Pf-Hsp70. The amount of protein is proportional to the volume peak. Numbers correspond to master gel and significant deregulated spots (i.e., 1314 and 1326) are reported in the table 2. (B) A graphic quantification of the four spots corresponding to Pf-Hsp 70 under normoxic (light-gray bars) and hyperoxic conditions (dark-gray bars). Spot numbers are specified at the bottom. Adjacent bars correspond to the mean standard abundance and respective standard deviation. Fold change and p-values are indicated at the bottom for each spot. A.U.: arbitrary units, FC: fold change.

To determine the metabolic pathways perturbed following hyperoxia, the identified proteins were classified using the NCBI COG database. The proportion of modulated proteins involved in each functional category was determined as follows. For P. falciparum, six functional categories were found to be altered, among which were chaperone-assisted protein folding, translation, antioxidant metabolism and glycolysis, which were already identified in transcriptomic analysis. For Homo sapiens, the identified proteins were classified into three functional categories: antioxidant metabolism, glycolysis and O2 transporter (Table 2).

The 2D-DIGE analyses also indicated the accumulation of some proteins involved in digestive vacuole metabolism such as human catalase [GenBank: gi|4557014] and beta-globin [GenBank: gi|183817]. The accumulation of beta-globin suggests proteases inhibition of the beta-globin degradation pathway. To explore this hypothesis, the transcripts of four genes [PlasmoDB: PF14_0077, PF11_0161, PF11_0165, PF11_0162] were quantified using real-time qRT-PCR as described above (Additional file 1). All genes involved in digestive vacuole metabolism were found to be down-expressed (Figure 6).

Figure 6
figure 6

Analysis of the beta-globin digestion pathway. qRT-PCR and microarray data of transcripts involved in digestive vacuole metabolism were compared between hyperoxic and normoxic conditions. Adjacent bars correspond to the mean log2 (fold change) and respective standard deviation and present qRT-PCR and microarray results in gray scale for the respective gene. The abbreviated names and their corresponding gene ID (PlasmoDB accession numbers) are indicated below the graphic. Crude values of mean fold changes are presented in the table. The same samples were used for the qRT-PCR and microarray experiments.

Discussion

In vivo, the P. falciparum parasite is subjected to varying oxygen levels throughout its life cycle (i.e., from 5% O2 in human venous blood to 13% O2 in the human lungs and 21% O2 in mosquito salivary glands). In the mosquito, the metabolic adaptation of parasite to oxygen-rich environment involved mitochondrial and physiological differences [32]. These oxygen variations imply that metabolic adaptation of P. falciparum is crucial for it survival. During malaria complications such as acute respiratory distress syndrome, late-ring stage parasites are susceptible to sequestration in pulmonary capillaries and are thus exposed to hyperoxic conditions [33]. These sequestered young parasites could be exposed to higher oxygen levels than the physiologically relevant O2 tension. Additionally, Blanco et al reported that hyperbaric oxygen therapy (HBO, 100% O2) has a beneficial effect on malaria syndrome evolution [34]. Indeed, a better understanding of the metabolic adaptation of the malaria parasite to hyperoxia could help to develop new anti-malarial drug treatments that could be used in association with HBO treatment. To study the global response of P. falciparum to hyperoxia, a dual high-throughput approach combining microarray and 2D-DIGE analysis was used on parasite cultures under hyperoxic conditions (e.g., 21% O2). Accordingly Hsp90, Hsp70, GAPDH, and elongation factor 1 and 2 were found to be altered at the transcript and protein levels under hyperoxia.

Since the development of high-throughput technologies, few studies have been published regarding the transcriptome and proteome of P. falciparum in response to environmental constraints or drug treatments [35–40]. Since hyperoxia induces a P. falciparum cycle delay, the sample collections were performed before the phenotypic effect. Although the majority of genes had a periodic expression profile [41], RNA transcription is maximal between 18 and 24 hours of the parasitic cycle [42]. In the present study, sample collections were performed at the late-ring parasite stage. Moreover, as a delay exists between mRNA and protein accumulation [43–45], the time of exposure was also taken into account and a four hours delay was chosen between mRNA and protein sampling. It has been suggested that there may be a discrepancy between P. falciparum transcriptomic and proteomic responses [35, 39]. Preliminary microarray experiments were performed with transcripts from synchronized P. falciparum in the same experimental scheme, but the samples were collected after four hours of hyperoxia treatment. These microarray analyses indicated that 176 transcripts were significantly deregulated (|FC|≤ 1.5, p ≤ 0.01, Student's t-Test), and some genes significantly deregulated were involved in the early antioxidant response, such as 1-cys peroxidoxin [PlasmoDB: PF08_0131], Fe-superoxide dismutase [PlasmoDB: PF08_0071] and thioredoxin peroxidase [PlasmoDB: MAL7P1.159] (Additional file 6). This early stress response is generally observed in stress condition does not reflect a specific hyperoxia adaptation [30, 46]. Here, despite significant transcript variations (p ≤ 0.01, Student's t-Test), gene fold-changes observed were low (|FC|≤ 1.5) under hyperoxia. The low-level changes observed in the P. falciparum transcriptome could be explained by tight gene regulation [47, 48] or by post-transcriptional regulation of most P. falciparum genes [43–45]. Consequently, analysis of the parasite's adaptive response to hyperoxia requires the use of extremely successful bioinformatic tools for microarray data interpretation such as PS 6.0 software [44, 49], and this analysis was completed using a highly sensitive proteomic approach such as 2D-DIGE.

It is generally accepted that one of the first effects of hyperoxia is ROS overproduction (superoxide anions (O-2), hydrogen peroxide (H2O2) and hydroxyl radicals (OH-)), which is generated by metabolism, and particularly, by respiratory metabolism [50]. DNA, lipid, and protein alterations by ROS may be lethal to malaria parasites. Thus, to fight oxidative stress, P. falciparum has developed an adaptive defence response including repair mechanisms for nucleic acids and proteins [29, 51]. Despite high ROS defence system expression under normoxia, transcriptomic analysis suggests an up-regulation of ROS defence systems [29]. Five genes involved in DNA repair were found to be up-regulated during hyperoxia. Among them, two sub-unities of DNA primase [PlasmoDB: PF14_0366 and PFI0530c] and replication protein A1 [PlasmoDB: rpa1, PFI0235w] have been described to be involved in chromosomal replication [52, 53]. Rpa1 was reported also to interact with rad51 [PlasmoDB: PF11_0087] in nucleosomes during replication to correct DNA mismatches [51, 54]. And Hsp40 [PlasmoDB: DNAJ homologue, PFF1415c] is associated with DNA repair and the replication machinery [55]. These observations suggest that DNA repair enzymes maintain the integrity of the parasitic genome under high oxygen pressure.

Protein oxidation caused by ROS is circumvented by diverse functions such as regulation of the redox state and modulation of protein stability [56]. Chaperone proteins, Hsps known as stress response proteins, further assure this protection. Several proteins involved in chaperone activity including Hsp40 [PlasmoDB: DNAJ homologue; PFF1415c], Hsp60 [PlasmoDB: PF10_0153], Hsp70 [PlasmoDB: PF11_0351 and PF08_0054], Hsp90 [PlasmoDB: PFL1070c and PF07_0029], and protein 14-3-3 [PlasmoDB: MAL8P1.69] were found up-regulated under hyperoxia. As described by Akide-Ndunge et al [57], Hsp60 is up-regulated under oxidative stress like hyperoxia and its expression is coordinated with antioxidant enzymes in a stage-dependent manner, suggesting thus that Hsp up-regulation is implicated in ROS removal. Elsewhere, Pf-Hsp70 forms a functional network in the mitochondrial matrix with DNAJ, Hsp60 and prohibitin to be involved in post-translational modification of proteins [56]. Isoforms of Hsp70 were detected following hyperoxia exposure, which may also correspond to post-translational modification as previously predicted [58, 59]. These Hsps, which act as sensors of environmental conditions, are involved in adaptation mechanism by post-translational modification [56]. However, the role of these post-translational modifications on regulation of protein expression in P. falciparum is little known [38, 39]. Nevertheless, the chaperone activity of Hsps seems regulate during the hyperoxia response and facilitate P. falciparum adaptation to hyperoxic environments.

Under hyperoxia, down-regulation of glycolytic enzymes (glucose-6-phosphate isomerase [PlasmoDB: PF14_0341] and glyceraldehyde-3-phosphate dehydrogenase [PlasmoDB: GAPDH, PF14_0598]) was detected in this study. Additionally, three enzyme involved in de novo pyrimidine biosynthesis were found up-regulated such as carbonic anhydrase [GenBank: gi|4502517], carbamoyl-phosphate synthetase [PlasmoDB: PF13_0044] and aspartate carbamoyltransferase [PlasmoDB: PF13_0240] [60]. This last, up-expressed in microarray analysis, produces dihydroorotate oxidase (DHO), an essential substrate of mitochondrial respiratory chain complex II [61]. This pyrimidine pathway is essential for nucleic acid synthesis to repair DNA lesions caused by the oxidative stress. Microarray data indicated that alpha and epsilon ATP synthase F1 subunits [PlasmoDB: PFB0795w and MAL7P1.75] were up-regulated, indicating that parasites seem to produce ATP through mitochondrial respiration. Mitochondrial F0F1 ATP synthase is composed of two subcomplexes, F0 and F1. F1 is composed of five subunits, and these have been reported in the P. falciparum genome [62]. Recently, Mogi and Kita have identified four F0 subunits of P. falciparum ATP synthase [63, 64] and Mather et al. support the existence of ATP synthase activity [65]. Two transcripts of ATP synthase F0 subcomplexe, a subunit (ATP6) and b subunit (ATP4), were up-regulated under hyperoxia exposure, which is in favour of a mitochondrial respiration. To adapt to hyperoxia, P. falciparum seems switch from anaerobic glycolysis to aerobic respiratory metabolism (Figure 7).

Figure 7
figure 7

A schematic representation of P. falciparum metabolic adaptation to hyperoxia exposure. Metabolic pathways based on the Ginsburg Pathway [72] for glycolysis (cytoplasmic), respiratory chain (mitochondria, blue and yellow) and haemoglobin digestion (digestive vacuolar, orange) are shown. The gene IDs (PlasmoDB) are reported for each enzymatic reaction. Human protein names are underlined. Up- and down-regulated genes and reactive oxygen species (ROS) production are indicated in red and blue characters, respectively.

In addition to mitochondrial metabolism, the digestive vacuole is another ROS source organelle in P. falciparum during haemoglobin digestion [10]. Haemoglobin digestion is optimal at an acidic pH, which is also necessary for protease activity. To provide the acidic environment in the digestive vacuole, haemoglobin digestion stimulates ATP consumption by the V-type H+-ATPase pump [66–68]. In microarray analysis, the V-type ATPase complex [PlasmoDB: MAL13P1.271, PF13_0034, PFE0965c, PF11_0412, PF13_0227, PF13_0065], a membrane transporter, was found to be down-regulated under hyperoxia. Therefore, P. falciparum exposed to hyperoxia could generate a pH change in the digestive vacuole responsible for vacuolar protease activity. Additionally, beta-globin accumulation occurred in hyperoxic conditions. This beta-globin accumulation could result from a decline in protease activity. To test these two hypotheses, transcripts of the proteases involved in beta-globin degradation were quantified. Reduced expression of four genes involved in haemoglobin degradation into AAs (plasmepsin 2 [PlasmoDB: PF14_0077], falcipain 2 putative [PlasmoDB: PF11_0161], falcipain 2 precursor [PlasmoDB: PF11_0165] and falcipain 3 [PlasmoDB: PF11_0162]) was validated using qRT-PCR. In 2002, Oliveira et al hypothesized that blood-feeding parasites reduced their mitochondrial function to compensate for ROS generation from the digestive vacuole [10].

Moreover, catalase [GenBank: gi|4557014], a human protein, was found concentrated in the digestive vacuole [69]. As the P. falciparum genome does not contain a catalase gene, the parasite may import human catalase to detoxify H2O2 generated by oxidation of haem under stress conditions. Collectively, metabolism in the P. falciparum digestive vacuole would be perturbed in response to hyperoxia, and ROS production would be slowed (Figure 7).

Conclusions

Two complementary analytic approaches were used to investigate the response of P. falciparum to hyperoxia; (i) a transcriptomic study allowed us to detect whole parasite transcripts, and (ii) a proteomic study identified proteins significantly altered via post-translational modifications and accumulated host proteins. Based on all these results and according to published data mining [29, 46], a schematic representation of the adaptive response of P. falciparum following hyperoxia exposure was proposed (Figure 7). In order to prove this representation, further biochemical approaches would be required.

Hyperoxia exposure induces metabolic adaptations in P. falciparum. These adaptations seem to involve, at least, two parasite organelles, the digestive vacuole and the mitochondrion, both sources of ROS production. To preserve parasite integrity from oxidative stress, all these data suggest that the glycolysis pathway is suppressed in favour of respiratory metabolism and that digestive vacuole metabolism is slowed. Campanale et al demonstrated that stress caused by haemoglobin digestion modulates the glycolytic pathway [70]. Highly active mitochondria release H2O2, which interacts with pro-oxidant products (free iron and haem) in the digestive vacuole. These two ROS sources could be potentially synergistic. The equilibrium of oxidative stress is vital for the parasite; indeed, Hsps could be regulated to facilitate adaptation of parasite to environmental stress as observed in many organisms [56].

The knowledge of the metabolic pathways involved in stress responses to environmental conditions is fundamental to understanding the mechanisms of parasite adaptation. This study provides a starting point for investigations into new anti-malarial treatments, particularly drugs associated with hyperbaric oxygen therapy [34], which has been successfully used to treat other infections [71].

Abbreviations

2D-DIGE:

2D-differential-gel-electrophoresis

2-DE:

2-Dimensional electrophoresis

AA:

amino acid

ACN:

acetonitrile

BVA:

biological variance analysis

Ct:

cycle threshold

CQ:

chloroquine

Cy-2:

cyanine 2

Cy-3:

cyanine 3

Cy-5:

cyanine 5

DHA:

dihydroartemisinin

DHO:

dihydroorotate

DHOD:

dihydroorotate dehydrogenase

DIA:

differential in-gel analysis

FC:

fold change

FIX:

ferriprotoporphyrin IX

gapdh :

glyceraldehyde-3-phosphate dehydrogenase

GSEA:

Gene Set Enrichment Analysis

gst :

glutathione S-transferase

H2O2:

hydrogen peroxide

HBO:

hyperbaric oxygen therapy

hsp :

heat shock protein

IEF:

isoelectric focusing

MS:

mass spectrometry

Nano LC-MS/MS:

nanoscale capillary liquid chromatography-tandem mass spectrometry

O-2:

superoxide anions

OH-:

hydroxyl radicals

O3:

ozone

P. falciparum :

Plasmodium falciparum

PMT:

photomultiplier tube

NPPs:

new permeation pathways

qRT-PCR:

quantitative real time PCR

Q-TOF:

quadrupole orthogonal acceleration time-of-flight

RBC:

red blood cells

ROS:

reactive oxygen species, rpa1: replication protein A1.

References

  1. WHO: Word Malaria Report. World Health Organization, Word Malaria Report. 2008

    Google Scholar 

  2. Greenwood BM, Fidock DA, Kyle DE, Kappe SH, Alonso PL, Collins FH, Duffy PE: Malaria: progress, perils, and prospects for eradication. J Clin Invest. 2008, 118: 1266-1276. 10.1172/JCI33996.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Torrentino-Madamet M, Desplans J, Travaillé C, Jammes Y, Parzy D: Microaerophilic respiratory metabolism of Plasmodium falciparum mitochondrion as a drug target. Current Molecular Medicine. 2010, 10: 29-46. 10.2174/156652410791065390.

    Article  CAS  PubMed  Google Scholar 

  4. Srivastava IK, Rottenberg H, Vaidya AB: Atovaquone, a broad spectrum antiparasitic drug, collapses mitochondrial membrane potential in a malarial parasite. J Biol Chem. 1997, 272: 3961-3966. 10.1074/jbc.272.52.33360.

    Article  CAS  PubMed  Google Scholar 

  5. Smilkstein MJ, Forquer I, Kanazawa A, Kelly JX, Winter RW, Hinrichs DJ, Kramer DM, Riscoe MK: A drug-selected Plasmodium falciparum lacking the need for conventional electron transport. Mol Biochem Parasitol. 2008, 159: 64-68. 10.1016/j.molbiopara.2008.01.002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Krungkrai J, Burat D, Kudan S, Krungkrai S, Prapunwattana P: Mitochondrial oxygen consumption in asexual and sexual blood stages of the human malarial parasite, Plasmodium falciparum. Southeast Asian J Trop Med Public Health. 1999, 30: 636-642.

    CAS  PubMed  Google Scholar 

  7. Krungkrai J, Prapunwattana P, Krungkrai SR: Ultrastructure and function of mitochondria in gametocytic stage of Plasmodium falciparum. Parasite. 2000, 7: 19-26.

    Article  CAS  PubMed  Google Scholar 

  8. Krungkrai J: The multiple roles of the mitochondrion of the malarial parasite. Parasitology. 2004, 129: 511-524. 10.1017/S0031182004005888.

    Article  CAS  PubMed  Google Scholar 

  9. Tsai AG, Johnson PC, Intaglietta M: Oxygen gradients in the microcirculation. Physiol Rev. 2003, 83: 933-963.

    Article  CAS  PubMed  Google Scholar 

  10. Oliveira PL, Oliveira MF: Vampires, Pasteur and reactive oxygen species. Is the switch from aerobic to anaerobic metabolism a preventive antioxidant defence in blood-feeding parasites?. FEBS Lett. 2002, 525: 3-6. 10.1016/S0014-5793(02)03026-0.

    Article  CAS  PubMed  Google Scholar 

  11. Turrens JF: Oxidative stress and antioxidant defenses: a target for the treatment of diseases caused by parasitic protozoa. Mol Aspects Med. 2004, 25: 211-220. 10.1016/j.mam.2004.02.021.

    Article  CAS  PubMed  Google Scholar 

  12. Lew VL, Macdonald L, Ginsburg H, Krugliak M, Tiffert T: Excess haemoglobin digestion by malaria parasites: a strategy to prevent premature host cell lysis. Blood Cells Mol Dis. 2004, 32: 353-359. 10.1016/j.bcmd.2004.01.006.

    Article  CAS  PubMed  Google Scholar 

  13. Wastling JM, Xia D, Sohal A, Chaussepied M, Pain A, Langsley G: Proteomes and transcriptomes of Apicomplexa - Where's the message?. Int J Parasitol. 2008, 39 (2): 135-43. 10.1016/j.ijpara.2008.10.003.

    Article  PubMed  Google Scholar 

  14. Trager W, Jensen JB: Human malaria parasites in continuous culture. J Parasitol. 1976, 91: 484-486. 10.1645/0022-3395(2005)091[0484:HMPICC]2.0.CO;2.

    Article  Google Scholar 

  15. Bogreau H, Renaud F, Bouchiba H, Durand P, Assi SB, Henry MC, Garnotel E, Pradines B, Fusai T, Wade B, Adehossi E, Parola P, Kamil MA, Puijalon O, Rogier C: Genetic diversity and structure of African Plasmodium falciparum populations in urban and rural areas. Am J Trop Med Hyg. 2006, 74: 953-959.

    CAS  PubMed  Google Scholar 

  16. Henry M, Diallo I, Bordes J, Ka S, Pradines B, Diatta B, M'Baye PS, Sane M, Thiam M, Gueye PM, Wade B, Touze JE, Debonne JM, Rogier C, Fusai T: Urban malaria in Dakar, Senegal: chemosusceptibility and genetic diversity of Plasmodium falciparum isolates. Am J Trop Med Hyg. 2006, 75: 146-151.

    CAS  PubMed  Google Scholar 

  17. Lambros C, Vanderberg JP: Synchronization of Plasmodium falciparum erythrocytic stages in culture. J Parasitol. 1979, 65: 418-420. 10.2307/3280287.

    Article  CAS  PubMed  Google Scholar 

  18. Uhlemann A-C, Trine S, Klinkerna M-Q, Hviidb L: Analysis of Plasmodium falciparum-infected red blood cells. MACS&more. 2000, 4: 7-8.

    Google Scholar 

  19. Weng L, Dai H, Zhan Y, He Y, Stepaniants SB, Bassett DE: Rosetta error model for gene expression analysis. Bioinformatics. 2006, 22: 1111-1121. 10.1093/bioinformatics/btl045.

    Article  CAS  PubMed  Google Scholar 

  20. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

  21. Pastorino B, Boucomont-Chapeaublanc E, Peyrefitte CN, Belghazi M, Fusai T, Rogier C, Tolou HJ, Almeras L: Identification of cellular proteome modifications in response to West Nile virus infection. Mol Cell Proteomics. 2009, 8: 1623-1637. 10.1074/mcp.M800565-MCP200.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Shevchenko A, Wilm M, Vorm O, Mann M: Mass spectrometric sequencing of proteins silver-stained polyacrylamide gels. Anal Chem. 1996, 68: 850-858. 10.1021/ac950914h.

    Article  CAS  PubMed  Google Scholar 

  23. Briolant S, Parola P, Fusai T, Madamet-Torrentino M, Baret E, Mosnier J, Delmont JP, Parzy D, Minodier P, Rogier C, Pradines B: Influence of oxygen on asexual blood cycle and susceptibility of Plasmodium falciparum to chloroquine: requirement of a standardized in vitro assay. Malar J. 2007, 6: 44-10.1186/1475-2875-6-44.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Gene Expression Omnibus. [http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/projects/geo/]

  25. Gene Set Enrichment Analysis. [http://www.broad.mit.edu/gsea]

  26. Yis U, Kurul SH, Kumral A, Cilaker S, Tugyan K, Genc S, Yilmaz O: Hyperoxic exposure leads to cell death in the developing brain. Brain Dev. 2008, 30: 556-562. 10.1016/j.braindev.2008.01.010.

    Article  PubMed  Google Scholar 

  27. Kappes B, Suetterlin BW, Hofer-Warbinek R, Humar R, Franklin RM: Two major phosphoproteins of Plasmodium falciparum are heat shock proteins. Mol Biochem Parasitol. 1993, 59: 83-94. 10.1016/0166-6851(93)90009-M.

    Article  CAS  PubMed  Google Scholar 

  28. 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: 163-189. 10.1016/j.ijpara.2003.09.011.

    Article  CAS  PubMed  Google Scholar 

  29. Guerra R, Shaw C, Christian B, Fox G, Noyola-Martinez J, Stevens M, Garg N, Gustin M: Oxidative stress genes in Plasmodium falciparum as indicated by temporal gene expression. Conference' 04. 2004, Month 1-2, 2004

    Google Scholar 

  30. Muller S: Redox and antioxidant systems of the malaria parasite Plasmodium falciparum. Mol Microbiol. 2004, 53: 1291-1305. 10.1111/j.1365-2958.2004.04257.x.

    Article  PubMed  Google Scholar 

  31. Nirmalan N, Wang P, Sims PF, Hyde JE: Transcriptional analysis of genes encoding enzymes of the folate pathway in the human malaria parasite Plasmodium falciparum. Mol Microbiol. 2002, 46: 179-190. 10.1046/j.1365-2958.2002.03148.x.

    Article  CAS  PubMed  Google Scholar 

  32. Learngaramkul P, Petmitr S, Krungkrai SR, Prapunwattana P, Krungkrai J: Molecular characterization of mitochondria in asexual and sexual blood stages of Plasmodium falciparum. Mol Cell Biol Res Commun. 1999, 2: 15-20. 10.1006/mcbr.1999.0145.

    Article  CAS  PubMed  Google Scholar 

  33. Maguire GP, Tjandra H, Michael CFP, Enny K, Ric N, Price , Emiliana T, Nicholas MA: Lung injury in uncomplicated and severe falciparum malaria: a longitudinal study in Papua, Indonesia. J Infect Dis. 2005, 1: 1966-1974. 10.1086/497697.

    Article  Google Scholar 

  34. Blanco YC, Farias AS, Goelnitz U, Lopes SC, Arrais-Silva WW, Carvalho BO, Amino R, Wunderlich G, Santos LM, Giorgio S, Costa FT: Hyperbaric oxygen prevents early death caused by experimental cerebral malaria. PLoS ONE. 2008, 3: e3126-10.1371/journal.pone.0003126.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Gunasekera AM, Myrick A, Le Roch K, Winzeler E, Wirth DF: Plasmodium falciparum: genome wide perturbations in transcript profiles among mixed stage cultures after chloroquine treatment. Exp Parasitol. 2007, 117: 87-92. 10.1016/j.exppara.2007.03.001.

    Article  CAS  PubMed  Google Scholar 

  36. 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: 2012-2025. 10.1128/IAI.01236-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Natalang O, Bischoff E, Deplaine G, Proux C, Dillies MA, Sismeiro O, Guigon G, Bonnefoy S, Patarapotikul J, Mercereau-Puijalon O, Coppee JY, David PH: Dynamic RNA profiling in Plasmodium falciparum synchronized blood stages exposed to lethal doses of artesunate. BMC Genomics. 2008, 9: 388-10.1186/1471-2164-9-388.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Le Roch KG, Johnson JR, Ahiboh H, Chung DW, Prudhomme J, Plouffe D, Henson K, Zhou Y, Witola W, Yates JR, Mamoun CB, Winzeler EA, Vial H: A systematic approach to understand the mechanism of action of the bisthiazolium compound T4 on the human malaria parasite, Plasmodium falciparum. BMC Genomics. 2008, 9: 513-10.1186/1471-2164-9-513.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Radfar A, Diez A, Bautista JM: Chloroquine mediates specific proteome oxidative damage across the erythrocytic cycle of resistant Plasmodium falciparum. Free Radic Biol Med. 2008, 44: 2034-2042. 10.1016/j.freeradbiomed.2008.03.010.

    Article  CAS  PubMed  Google Scholar 

  40. Prieto JH, Koncarevic S, Park SK, Yates J, Becker K: Large-scale differential proteome analysis in Plasmodium falciparum under drug treatment. PLoS ONE. 2008, 3: 4098-10.1371/journal.pone.0004098.

    Article  Google Scholar 

  41. Bozdech Z, Zhu J, Joachimiak MP, Cohen FE, Pulliam B, DeRisi JL: Expression profiling of the schizont and trophozoite stages of Plasmodium falciparum with a long-oligonucleotide microarray. Genome Biol. 2003, 4: 9-10.1186/gb-2003-4-2-r9.

    Article  Google Scholar 

  42. 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: 5-10.1371/journal.pbio.0000005.

    Article  Google Scholar 

  43. Foth BJ, Zhang N, Mok S, Preiser PR, Bozdech Z: Quantitative protein expression profiling reveals extensive post-transcriptional regulation and post-translational modifications in schizont-stage malaria parasites. Genome Biol. 2008, 9: 177-10.1186/gb-2008-9-12-r177.

    Article  Google Scholar 

  44. LaCount DJ, Vignali M, Chettier R, Phansalkar A, Bell R, Hesselberth JR, Schoenfeld LW, Ota I, Sahasrabudhe S, Kurschner C, Fields S, Hughes RE: A protein interaction network of the malaria parasite Plasmodium falciparum. Nature. 2005, 438: 103-107. 10.1038/nature04104.

    Article  CAS  PubMed  Google Scholar 

  45. Le Roch KG, Johnson JR, Florens L, Zhou Y, Santrosyan A, Grainger M, Yan SF, Williamson KC, Holder AA, Carucci DJ, Yates JR, Winzeler EA: Global analysis of transcript and protein levels across the Plasmodium falciparum life cycle. Genome Res. 2004, 14: 2308-2318. 10.1101/gr.2523904.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Bozdech Z, Ginsburg H: Antioxidant defense in Plasmodium falciparum--data mining of the transcriptome. Malar J. 2004, 3: 23-10.1186/1475-2875-3-23.

    Article  PubMed Central  PubMed  Google Scholar 

  47. Horrocks P, Wong E, Russell K, Emes RD: Control of gene expression in Plasmodium falciparum - ten years on. Mol Biochem Parasitol. 2009, 164: 9-25. 10.1016/j.molbiopara.2008.11.010.

    Article  CAS  PubMed  Google Scholar 

  48. 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 intra-erythrocytic development cycle. Genome Biol. 2007, 8: 134-10.1186/gb-2007-8-7-r134.

    Article  Google Scholar 

  49. Mash DC, ffrench-Mullen J, Adi N, Qin Y, Buck A, Pablo J: Gene expression in human hippocampus from cocaine abusers identifies genes which regulate extracellular matrix remodeling. PLoS ONE. 2007, 2: 1187-10.1371/journal.pone.0001187.

    Article  Google Scholar 

  50. Andreyev AY, Kushnareva YE, Starkov AA: Mitochondrial metabolism of reactive oxygen species. Biochemistry (Mosc). 2005, 70: 200-214. 10.1007/s10541-005-0102-7.

    Article  CAS  Google Scholar 

  51. Bhattacharyya MK, Kumar N: Identification and molecular characterisation of DNA damaging agent induced expression of Plasmodium falciparum recombination protein PfRad51. Int J Parasitol. 2003, 33: 1385-1392. 10.1016/S0020-7519(03)00212-1.

    Article  CAS  PubMed  Google Scholar 

  52. Seow F, Sato S, Janssen CS, Riehle MO, Mukhopadhyay A, Phillips RS, Wilson RJ, Barrett MP: The plastidic DNA replication enzyme complex of Plasmodium falciparum. Mol Biochem Parasitol. 2005, 141: 145-153. 10.1016/j.molbiopara.2005.02.002.

    Article  CAS  PubMed  Google Scholar 

  53. Rider SD, Cai X, Sullivan WJ, Smith AT, Radke J, White M, Zhu G: The protozoan parasite Cryptosporidium parvum possesses two functionally and evolutionarily divergent replication protein A large subunits. J Biol Chem. 2005, 280: 31460-31469. 10.1074/jbc.M504466200.

    Article  CAS  PubMed  Google Scholar 

  54. Sugiyama T, Kantake N: Dynamic regulatory interactions of rad51, rad52, and replication protein-a in recombination intermediates. J Mol Biol. 2009, 390: 45-55. 10.1016/j.jmb.2009.05.009.

    Article  CAS  PubMed  Google Scholar 

  55. 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: 1781-1803. 10.1016/j.biocel.2007.02.011.

    Article  CAS  PubMed  Google Scholar 

  56. Acharya P, Kumar R, Tatu U: Chaperoning a cellular upheaval in malaria: heat shock proteins in Plasmodium falciparum. Mol Biochem Parasitol. 2007, 153: 85-94. 10.1016/j.molbiopara.2007.01.009.

    Article  CAS  PubMed  Google Scholar 

  57. Akide-Ndunge OB, Tambini E, Giribaldi G, McMillan PJ, Muller S, Arese P, Turrini F: Co-ordinated stage-dependent enhancement of Plasmodium falciparum antioxidant enzymes and heat shock protein expression in parasites growing in oxidatively stressed or G6PD-deficient red blood cells. Malar J. 2009, 8: 113-10.1186/1475-2875-8-113.

    Article  PubMed Central  PubMed  Google Scholar 

  58. Shonhai A, Boshoff A, Blatch GL: The structural and functional diversity of Hsp70 proteins from Plasmodium falciparum. Protein Sci. 2007, 16: 1803-1818. 10.1110/ps.072918107.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  59. McCarty JS, Walker GC: DnaK as a thermometer: threonine-199 is site of autophosphorylation and is critical for ATPase activity. Proc Natl Acad Sci USA. 1991, 88: 9513-9517. 10.1073/pnas.88.21.9513.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Krungkrai SR, Suraveratum N, Rochanakij S, Krungkrai J: Characterisation of carbonic anhydrase in Plasmodium falciparum. Int J Parasitol. 2001, 31: 661-668. 10.1016/S0020-7519(01)00172-2.

    Article  CAS  PubMed  Google Scholar 

  61. Krungkrai J: Purification, characterization and localization of mitochondrial dihydroorotate dehydrogenase in Plasmodium falciparum, human malaria parasite. Biochim Biophys Acta. 1995, 1243: 351-360.

    Article  PubMed  Google Scholar 

  62. Vaidya AB, Mather MW: A post-genomic view of the mitochondrion in malaria parasites. Curr Top Microbiol Immunol. 2005, 295: 233-250. full_text.

    CAS  PubMed  Google Scholar 

  63. Kawahara K, Mogi T, Tanaka TQ, Hata M, Miyoshi H, Kita K: Mitochondrial Dehydrogenases in the aerobic respiratory chain of the rodent malaria parasite Plasmodium yoelii yoelii. J Biochem. 2009, 145: 229-237. 10.1093/jb/mvn161.

    Article  CAS  PubMed  Google Scholar 

  64. Mogi T, Kita K: Identification of mitochondrial Complex II subunits SDH3 and SDH4 and ATP synthase subunits a and b in Plasmodium spp. Mitochondrion. 2009, 9: 443-453. 10.1016/j.mito.2009.08.004.

    Article  CAS  PubMed  Google Scholar 

  65. Mather MW, Morrisey JM, Vaidya AB: Hemozoin-free Plasmodium falciparum mitochondria for physiological and drug susceptibility studies. Mol Biochem Parasitol. 2010, 174: 150-153. 10.1016/j.molbiopara.2010.07.006.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. Beyenbach KW, Wieczorek H: The V-type H+ ATPase: molecular structure and function, physiological roles and regulation. J Exp Biol. 2006, 209: 577-589. 10.1242/jeb.02014.

    Article  CAS  PubMed  Google Scholar 

  67. Saliba KJ, Allen RJ, Zissis S, Bray PG, Ward SA, Kirk K: Acidification of the malaria parasite's digestive vacuole by a H+-ATPase and a H+-pyrophosphatase. J Biol Chem. 2003, 278: 5605-5612. 10.1074/jbc.M208648200.

    Article  CAS  PubMed  Google Scholar 

  68. Hayashi M, Yamada H, Mitamura T, Horii T, Yamamoto A, Moriyama Y: Vacuolar H(+)-ATPase localized in plasma membranes of malaria parasite cells, Plasmodium falciparum, is involved in regional acidification of parasitized erythrocytes. J Biol Chem. 2000, 275: 34353-34358. 10.1074/jbc.M003323200.

    Article  CAS  PubMed  Google Scholar 

  69. Koncarevic S, Rohrbach P, Deponte M, Krohne G, Prieto JH, Yates J, Rahlfs S, Becker K: The malarial parasite Plasmodium falciparum imports the human protein peroxiredoxin 2 for peroxide detoxification. Proc Natl Acad Sci USA. 2009, 106: 13323-13328. 10.1073/pnas.0905387106.

    Article  PubMed Central  PubMed  Google Scholar 

  70. Campanale N, Nickel C, Daubenberger CA, Wehlan DA, Gorman JJ, Klonis N, Becker K, Tilley L: Identification and characterization of heme-interacting proteins in the malaria parasite, Plasmodium falciparum. J Biol Chem. 2003, 278: 27354-27361. 10.1074/jbc.M303634200.

    Article  CAS  PubMed  Google Scholar 

  71. Perkowski S, Sun J, Singhal S, Santiago J, Leikauf GD, Albelda SM: Gene expression profiling of the early pulmonary response to hyperoxia in mice. Am J Respir Cell Mol Biol. 2003, 28: 682-696. 10.1165/rcmb.4692.

    Article  CAS  PubMed  Google Scholar 

  72. Ginsburg Pathway. [http://sites.huji.ac.il/malaria]

Download references

Acknowledgements

We would like to thank N. Benoit (UMR-MD3, IMTSSA, France) and C. Chapus (UMR-MD3, IMTSSA, France) for technical help. We thank the technical support team from Applied Biosystems for qRT-PCR optimization and the technical support team from Ariadne Genomics (PS 6.0™) for their help in the interpretation of microarray analysis. We acknowledge the financial support of the DGA (Délégation Générale pour l'Armement) and the Direction Centrale du Service de Santé des Armées (grant no. 06CO008). This manuscript has been reviewed and corrected by a journal expert.

All authors declare that they have no conflicts of interest.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Marylin Torrentino-Madamet.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MTM conceived the study and the design, carried out microarray and proteomic studies, participated in bioinformatics analyses and wrote the manuscript. LA carried out proteomic studies, conducted proteomic statistical analyses and revised the manuscript. JD helped in the design of molecular studies and revised the manuscript. YL participated in the microarray design and conducted statistical and bioinformatics analyses. MB, MP and PF carried out mass spectrometry identifications. YJ conceived the study and edited the manuscript. DP initiated the project, designed the method, participated in the analyses and revised the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Primers sequence using real-time qRT-PCR. (DOC 42 KB)

Additional file 2: Single-Peptide-Based Protein Identifications. (PPT 166 KB)

12936_2010_1587_MOESM3_ESM.XLS

Additional file 3: Raw data microarray at 8 hours time point under hyperoxia versus normoxia conditions on synchronized parasites. (XLS 272 KB)

Additional file 4: List of altered genes following hyperoxia treatment using GSEA data and PS 6.0 software. (XLS 24 KB)

12936_2010_1587_MOESM5_ESM.DOC

Additional file 5: MS/MS peptide sequences, respective gi number, gene ID and master spot number of proteins identified from the differential 2-D DIGE analysis following hyperoxia exposure of P. falciparum. (DOC 62 KB)

12936_2010_1587_MOESM6_ESM.XLS

Additional file 6: Raw data microarray at 4 hours time point under hyperoxia versus normoxia conditions on synchronized parasites. (XLS 60 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to 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.

Reprints and permissions

About this article

Cite this article

Torrentino-Madamet, M., Alméras, L., Desplans, J. et al. Global response of Plasmodium falciparum to hyperoxia: a combined transcriptomic and proteomic approach. Malar J 10, 4 (2011). https://0-doi-org.brum.beds.ac.uk/10.1186/1475-2875-10-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1475-2875-10-4

Keywords