Document Detail

Global response of Plasmodium falciparum to hyperoxia: a combined transcriptomic and proteomic approach.
Jump to Full Text
MedLine Citation:
PMID:  21223545     Owner:  NLM     Status:  MEDLINE    
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.
Marylin Torrentino-Madamet; Lionel Alméras; Jérôme Desplans; Yannick Le Priol; Maya Belghazi; Matthieu Pophillat; Patrick Fourquet; Yves Jammes; Daniel Parzy
Related Documents :
21529015 - A high-pressure atomic force microscope for imaging in supercritical carbon dioxide.
21276995 - Dogs mean arterial pressure and heart rate responses during high propofol plasma concen...
8943195 - Polymorphs of alumina predicted by first principles: putting pressure on the ruby press...
21501075 - Evaluation of rebound tonometer for measuring intraocular pressure at deviated angle an...
11411245 - Measurement of intercompartmental fluid shifts during haemodialysis in children.
309265 - Phentolamine infusion in cor pulmonale due to chronic obstructive pulmonary disease.
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't     Date:  2011-01-11
Journal Detail:
Title:  Malaria journal     Volume:  10     ISSN:  1475-2875     ISO Abbreviation:  Malar. J.     Publication Date:  2011  
Date Detail:
Created Date:  2011-01-31     Completed Date:  2011-04-27     Revised Date:  2013-07-02    
Medline Journal Info:
Nlm Unique ID:  101139802     Medline TA:  Malar J     Country:  England    
Other Details:
Languages:  eng     Pagination:  4     Citation Subset:  IM    
UMR-MD3 (Université de la Méditerranée), Antenne IRBA de Marseille (IMTSSA, Le Pharo), Allée du Médecin Colonel Eugène Jamot, BP 60109, 13262 Marseille cedex 07, France.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Antioxidants / metabolism
Gene Expression Profiling*
Oxidative Stress
Oxygen / metabolism*
Plasmodium falciparum / drug effects*
Proteome / analysis*
Protozoan Proteins / biosynthesis
Stress, Physiological
Reg. No./Substance:
0/Antioxidants; 0/Proteome; 0/Protozoan Proteins; 7782-44-7/Oxygen

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Malar J
ISSN: 1475-2875
Publisher: BioMed Central
Article Information
Download PDF
Copyright ©2011 Torrentino-Madamet et al; licensee BioMed Central Ltd.
Received Day: 29 Month: 9 Year: 2010
Accepted Day: 11 Month: 1 Year: 2011
collection publication date: Year: 2011
Electronic publication date: Day: 11 Month: 1 Year: 2011
Volume: 10First Page: 4 Last Page: 4
ID: 3030542
Publisher Id: 1475-2875-10-4
PubMed Id: 21223545
DOI: 10.1186/1475-2875-10-4

Global response of Plasmodium falciparum to hyperoxia: a combined transcriptomic and proteomic approach
Marylin Torrentino-Madamet1 Email:
Lionel Alméras2 Email:
Jérôme Desplans1 Email:
Yannick Le Priol1 Email:
Maya Belghazi3 Email:
Matthieu Pophillat4 Email:
Patrick Fourquet4 Email:
Yves Jammes5 Email:
Daniel Parzy1 Email:
1UMR-MD3 (Université de la Méditerranée), Antenne IRBA de Marseille (IMTSSA, Le Pharo), Allée du Médecin Colonel Eugène Jamot, BP 60109, 13262 Marseille cedex 07, France
2Unité de Recherche en Biologie et Epidémiologie Parasitaires (URBEP), Antenne IRBA de Marseille (IMTSSA, Le Pharo), Allée du Médecin Colonel Eugène Jamot, BP 60109, 13262 Marseille cedex 07, France
3Centre d'Analyse Protéomique de Marseille (CAPM), Institut Fédératif de Recherche Jean Roche, Faculté de Médecine Nord, Bd. Pierre Dramard, 13916 Marseille cedex 20, France
4Centre d'Immunologie de Marseille Luminy (CIML), Institut National de la Santé et de la Recherche Médicale, Centre National de la Recherche Scientifique, Université de la Méditerranée, Parc Scientifique de Luminy, 13288 Marseille Cedex 09, Marseille, France
5UMR-MD2, Physiologie et Physiopathologie en Conditions d'Oxygénations Extrêmes, Institut Fédératif de Recherche Jean Roche, Faculté de Médecine Nord, Bd. Pierre Dramard, 13916 Marseille cedex 20, France


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.

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 = (Cttarget - Ct18S rRNA)H- (Cttarget - Ct18S 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).

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.

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].

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).

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).

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.

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).


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).

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).


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].

List of 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.

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.

Supplementary Material Additional file 1

Primers sequence using real-time qRT-PCR.

Click here for additional data file (1475-2875-10-4-S1.DOC)

Additional file 2

Single-Peptide-Based Protein Identifications.

Click here for additional data file (1475-2875-10-4-S2.PPT)

Additional file 3

Raw data microarray at 8 hours time point under hyperoxia versus normoxia conditions on synchronized parasites.

Click here for additional data file (1475-2875-10-4-S3.XLS)

Additional file 4

List of altered genes following hyperoxia treatment using GSEA data and PS 6.0 software.

Click here for additional data file (1475-2875-10-4-S4.XLS)

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.

Click here for additional data file (1475-2875-10-4-S5.DOC)

Additional file 6

Raw data microarray at 4 hours time point under hyperoxia versus normoxia conditions on synchronized parasites.

Click here for additional data file (1475-2875-10-4-S6.XLS)


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.

WHOWord Malaria ReportWorld Health Organization, Word Malaria ReportYear: 2008
Greenwood BM,Fidock DA,Kyle DE,Kappe SH,Alonso PL,Collins FH,Duffy PE,Malaria: progress, perils, and prospects for eradicationJ Clin InvestYear: 20081181266127610.1172/JCI3399618382739
Torrentino-Madamet M,Desplans J,Travaillé C,Jammes Y,Parzy D,Microaerophilic respiratory metabolism of Plasmodium falciparum mitochondrion as a drug targetCurrent Molecular MedicineYear: 201010294610.2174/15665241079106539020205678
Srivastava IK,Rottenberg H,Vaidya AB,Atovaquone, a broad spectrum antiparasitic drug, collapses mitochondrial membrane potential in a malarial parasiteJ Biol ChemYear: 19972723961396610.1074/jbc.272.52.333609020100
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 transportMol Biochem ParasitolYear: 2008159646810.1016/j.molbiopara.2008.01.00218308406
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 falciparumSoutheast Asian J Trop Med Public HealthYear: 19993063664210928353
Krungkrai J,Prapunwattana P,Krungkrai SR,Ultrastructure and function of mitochondria in gametocytic stage of Plasmodium falciparumParasiteYear: 20007192610743643
Krungkrai J,The multiple roles of the mitochondrion of the malarial parasiteParasitologyYear: 200412951152410.1017/S003118200400588815552397
Tsai AG,Johnson PC,Intaglietta M,Oxygen gradients in the microcirculationPhysiol RevYear: 20038393396312843412
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 LettYear: 20025253610.1016/S0014-5793(02)03026-012163151
Turrens JF,Oxidative stress and antioxidant defenses: a target for the treatment of diseases caused by parasitic protozoaMol Aspects MedYear: 20042521122010.1016/j.mam.2004.02.02115051329
Lew VL,Macdonald L,Ginsburg H,Krugliak M,Tiffert T,Excess haemoglobin digestion by malaria parasites: a strategy to prevent premature host cell lysisBlood Cells Mol DisYear: 20043235335910.1016/j.bcmd.2004.01.00615121091
Wastling JM,Xia D,Sohal A,Chaussepied M,Pain A,Langsley G,Proteomes and transcriptomes of Apicomplexa - Where's the message?Int J ParasitolYear: 20083921354310.1016/j.ijpara.2008.10.00318996390
Trager W,Jensen JB,Human malaria parasites in continuous cultureJ ParasitolYear: 19769148448610.1645/0022-3395(2005)091[0484:HMPICC]2.0.CO;2
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 areasAm J Trop Med HygYear: 20067495395916760503
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 isolatesAm J Trop Med HygYear: 20067514615116837722
Lambros C,Vanderberg JP,Synchronization of Plasmodium falciparum erythrocytic stages in cultureJ ParasitolYear: 19796541842010.2307/3280287383936
Uhlemann A-C,Trine S,Klinkerna M-Q,Hviidb L,Analysis of Plasmodium falciparum-infected red blood cellsMACS&moreYear: 2000478
Weng L,Dai H,Zhan Y,He Y,Stepaniants SB,Bassett DE,Rosetta error model for gene expression analysisBioinformaticsYear: 2006221111112110.1093/bioinformatics/btl04516522673
Livak KJ,Schmittgen TD,Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) MethodMethodsYear: 20012540240810.1006/meth.2001.126211846609
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 infectionMol Cell ProteomicsYear: 200981623163710.1074/mcp.M800565-MCP20019395707
Shevchenko A,Wilm M,Vorm O,Mann M,Mass spectrometric sequencing of proteins silver-stained polyacrylamide gelsAnal ChemYear: 19966885085810.1021/ac950914h8779443
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 assayMalar JYear: 200764410.1186/1475-2875-6-4417437625
Gene Expression Omnibus
Gene Set Enrichment Analysis
Yis U,Kurul SH,Kumral A,Cilaker S,Tugyan K,Genc S,Yilmaz O,Hyperoxic exposure leads to cell death in the developing brainBrain DevYear: 20083055656210.1016/j.braindev.2008.01.01018329209
Kappes B,Suetterlin BW,Hofer-Warbinek R,Humar R,Franklin RM,Two major phosphoproteins of Plasmodium falciparum are heat shock proteinsMol Biochem ParasitolYear: 199359839410.1016/0166-6851(93)90009-M8515785
Becker K,Tilley L,Vennerstrom JL,Roberts D,Rogerson S,Ginsburg H,Oxidative stress in malaria parasite-infected erythrocytes: host-parasite interactionsInt J ParasitolYear: 20043416318910.1016/j.ijpara.2003.09.01115037104
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 expressionConference' 04Year: 2004 Month 1-2, 2004.
Muller S,Redox and antioxidant systems of the malaria parasite Plasmodium falciparumMol MicrobiolYear: 2004531291130510.1111/j.1365-2958.2004.04257.x15387810
Nirmalan N,Wang P,Sims PF,Hyde JE,Transcriptional analysis of genes encoding enzymes of the folate pathway in the human malaria parasite Plasmodium falciparumMol MicrobiolYear: 20024617919010.1046/j.1365-2958.2002.03148.x12366841
Learngaramkul P,Petmitr S,Krungkrai SR,Prapunwattana P,Krungkrai J,Molecular characterization of mitochondria in asexual and sexual blood stages of Plasmodium falciparumMol Cell Biol Res CommunYear: 19992152010.1006/mcbr.1999.014510527885
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, IndonesiaJ Infect DisYear: 200511966197410.1086/497697
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 malariaPLoS ONEYear: 20083e312610.1371/journal.pone.000312618769544
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 treatmentExp ParasitolYear: 2007117879210.1016/j.exppara.2007.03.00117475254
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 parasitesInfect ImmunYear: 2007752012202510.1128/IAI.01236-0617283083
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 artesunateBMC GenomicsYear: 2008938810.1186/1471-2164-9-38818706115
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 falciparumBMC GenomicsYear: 2008951310.1186/1471-2164-9-51318973684
Radfar A,Diez A,Bautista JM,Chloroquine mediates specific proteome oxidative damage across the erythrocytic cycle of resistant Plasmodium falciparumFree Radic Biol MedYear: 2008442034204210.1016/j.freeradbiomed.2008.03.01018397762
Prieto JH,Koncarevic S,Park SK,Yates J,Becker K,Large-scale differential proteome analysis in Plasmodium falciparum under drug treatmentPLoS ONEYear: 20083409810.1371/journal.pone.0004098
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 microarrayGenome BiolYear: 20034910.1186/gb-2003-4-2-r9
Bozdech Z,Llinas M,Pulliam BL,Wong ED,Zhu J,DeRisi JL,The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparumPLoS BiolYear: 20031510.1371/journal.pbio.0000005
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 parasitesGenome BiolYear: 2008917710.1186/gb-2008-9-12-r177
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 falciparumNatureYear: 200543810310710.1038/nature0410416267556
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 cycleGenome ResYear: 2004142308231810.1101/gr.252390415520293
Bozdech Z,Ginsburg H,Antioxidant defense in Plasmodium falciparum--data mining of the transcriptomeMalar JYear: 200432310.1186/1475-2875-3-2315245577
Horrocks P,Wong E,Russell K,Emes RD,Control of gene expression in Plasmodium falciparum - ten years onMol Biochem ParasitolYear: 200916492510.1016/j.molbiopara.2008.11.01019110008
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 cycleGenome BiolYear: 2007813410.1186/gb-2007-8-7-r134
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 remodelingPLoS ONEYear: 20072118710.1371/journal.pone.0001187
Andreyev AY,Kushnareva YE,Starkov AA,Mitochondrial metabolism of reactive oxygen speciesBiochemistry (Mosc)Year: 20057020021410.1007/s10541-005-0102-715807660
Bhattacharyya MK,Kumar N,Identification and molecular characterisation of DNA damaging agent induced expression of Plasmodium falciparum recombination protein PfRad51Int J ParasitolYear: 2003331385139210.1016/S0020-7519(03)00212-114527521
Seow F,Sato S,Janssen CS,Riehle MO,Mukhopadhyay A,Phillips RS,Wilson RJ,Barrett MP,The plastidic DNA replication enzyme complex of Plasmodium falciparumMol Biochem ParasitolYear: 200514114515310.1016/j.molbiopara.2005.02.00215926203
Rider SD Jr,Cai X,Sullivan WJ Jr,Smith AT,Radke J,White M,Zhu G,The protozoan parasite Cryptosporidium parvum possesses two functionally and evolutionarily divergent replication protein A large subunitsJ Biol ChemYear: 2005280314603146910.1074/jbc.M50446620016014411
Sugiyama T,Kantake N,Dynamic regulatory interactions of rad51, rad52, and replication protein-a in recombination intermediatesJ Mol BiolYear: 2009390455510.1016/j.jmb.2009.05.00919445949
Botha M,Pesce ER,Blatch GL,The Hsp40 proteins of Plasmodium falciparum and other apicomplexa: regulating chaperone power in the parasite and the hostInt J Biochem Cell BiolYear: 2007391781180310.1016/j.biocel.2007.02.01117428722
Acharya P,Kumar R,Tatu U,Chaperoning a cellular upheaval in malaria: heat shock proteins in Plasmodium falciparumMol Biochem ParasitolYear: 2007153859410.1016/j.molbiopara.2007.01.00917307260
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 cellsMalar JYear: 2009811310.1186/1475-2875-8-11319480682
Shonhai A,Boshoff A,Blatch GL,The structural and functional diversity of Hsp70 proteins from Plasmodium falciparumProtein SciYear: 2007161803181810.1110/ps.07291810717766381
McCarty JS,Walker GC,DnaK as a thermometer: threonine-199 is site of autophosphorylation and is critical for ATPase activityProc Natl Acad Sci USAYear: 1991889513951710.1073/pnas.88.21.95131835085
Krungkrai SR,Suraveratum N,Rochanakij S,Krungkrai J,Characterisation of carbonic anhydrase in Plasmodium falciparumInt J ParasitolYear: 20013166166810.1016/S0020-7519(01)00172-211336746
Krungkrai J,Purification, characterization and localization of mitochondrial dihydroorotate dehydrogenase in Plasmodium falciparum, human malaria parasiteBiochim Biophys ActaYear: 199512433513607727509
Vaidya AB,Mather MW,A post-genomic view of the mitochondrion in malaria parasitesCurr Top Microbiol ImmunolYear: 2005295233250 full_text. 16265893
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 yoeliiJ BiochemYear: 200914522923710.1093/jb/mvn16119060309
Mogi T,Kita K,Identification of mitochondrial Complex II subunits SDH3 and SDH4 and ATP synthase subunits a and b in Plasmodium sppMitochondrionYear: 2009944345310.1016/j.mito.2009.08.00419682605
Mather MW,Morrisey JM,Vaidya AB,Hemozoin-free Plasmodium falciparum mitochondria for physiological and drug susceptibility studiesMol Biochem ParasitolYear: 201017415015310.1016/j.molbiopara.2010.07.00620674615
Beyenbach KW,Wieczorek H,The V-type H+ ATPase: molecular structure and function, physiological roles and regulationJ Exp BiolYear: 200620957758910.1242/jeb.0201416449553
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+-pyrophosphataseJ Biol ChemYear: 20032785605561210.1074/jbc.M20864820012427765
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 erythrocytesJ Biol ChemYear: 2000275343533435810.1074/jbc.M00332320010915784
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 detoxificationProc Natl Acad Sci USAYear: 2009106133231332810.1073/pnas.090538710619666612
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 falciparumJ Biol ChemYear: 2003278273542736110.1074/jbc.M30363420012748176
Perkowski S,Sun J,Singhal S,Santiago J,Leikauf GD,Albelda SM,Gene expression profiling of the early pulmonary response to hyperoxia in miceAm J Respir Cell Mol BiolYear: 20032868269610.1165/rcmb.469212760966
Ginsburg Pathway


[Figure ID: F1]
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).

[Figure ID: F2]
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.

[Figure ID: F3]
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.

[Figure ID: F4]
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.

[Figure ID: F5]
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.

[Figure ID: F6]
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.

[Figure ID: F7]
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.

[TableWrap ID: T1] Table 1 

Biological functions perturbed following hyperoxia exposure on P. falciparum

Functional groupa Number of entitiesb p-valuec
GPI anchor biosynthesis 29 0.0006
DNA repair 15 0.0072
Vacuolar acidification 11 0.0255
Actin filament organization 12 0.0296
Nucleosome assembly 12 0.0490
Regulation of cell shape 6 0.0519
Leading strand elongation 10 0.0533
Lysosomal H+ import 26 0.0589
Response to oxidative stress 17 0.0799

aGene Set Enrichment Analysis (GSEA) using the Mann-Whitney-U-test enrichment algorithm in PathwayStudio software indicated functional groups significantly altered (p < 0.08). For each functional group, the bnumber of genes included and corresponding cp-values are listed.

[TableWrap ID: T2] Table 2 

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

gi numberc Gene ID Protein name MW (kDa) pI Master spot number Significance (Mascot score) Average ratio t-Test
P. falciparum


gi|124512420 MAL8P1.69 14-3-3 protein homologue 29.86 4.96 3289 425 -1.54 0.0078

gi|124513850 PF13_0304 elongation factor 1 alpha 49.16 9.12 2047 130 1.60 0.0090

gi|8918238 PF14_0486 elongation factor 2 85.03 6.30 1030 67 1.72 0.0070

gi|124810293 PF14_0655 RNA helicase-1 45.62 5.48 2268 105 1.53 0.040

Parasitophorous vacuolar membrane Transporter

gi|124810348 PF14_0678 exported protein 2 33.62 5.10 3061 88 -1.92 0.00093

3062 180 -2.06 0.0046


gi|124809201 PF14_0341 glucose-6-phosphate isomerase 67.61 6.78 1579 36 -1.63 0.00092

gi|124810131 PF14_0598 glyceraldehyde-3-phosphate dehydrogenase 37.08 7.59 2867 695 -1.69 0.033

Chaperone-assisted protein folding

gi|124512406 PF08_0054 heat shock protein 70 74.39 5.51 1314 257 1.65 0.03

1326 137 1.58 0.026

gi|505340 PF07_0029 heat shock protein 86 86.77 4.91 1028 412 1.50 0.0052

Amino acids metabolism

gi|86170756 PFF0435w ornithine aminotransferase 47 6.47 2217 96 -2.14 0.0081

gi|124513590 MAL13P1.214 phosphoethanolamine 31.31 5.43 3347 145 -2.05 0.0055

N-methyltransferase 3372 252 -1.96 0.00089

3379 253 -2.17 0.00058

3503 52 -1.62 0.00056

3385 452 -2.71 0.00026

Proteasome-mediated proteolysis

gi|124512686 MAL8P1.142 proteasome beta-subunit 31.08 6.00 3216 58 1.78 0.015

gi|124513790 MAL13P1.270 proteasome subunit 27.50 6.17 3358 142 1.55 0.049

Homo sapiens

Oxygen transporter

gi|183817 Beta-globin 19.21 6.28 3451 165 7.71 0.020

3576 102 1.91 0.011

Antioxidant metabolism

gi|4502517 carbonic anhydrase I 28.91 6.59 3431 173 4.18 0.0095

3444 155 3.47 0.0019

3455 224 4.02 0.00038

gi|4557014 catalase 59.95 6.9 1648 200 6.47 0.012

gi|16306550 selenium binding protein 1 52.93 5.93 2451 73 -1.64 0.041

gi|168985379 flotillin 1 39.81 6.03 2156 157 1.64 0.040


gi|31645 glyceraldehyde-3- 36.20 8.26 2804 199 -1.94 0.0042

phosphate dehydrogenase 2806 181 -1.99 0.012

2807 132 -1.84 0.0069

2853 112 -2.17 0.0027

2854 122 -2.25 0.0032

2856 96 -2.20 0.0019

The proteins were identified by mass spectrometry following in gel trypsin digestion. The spot numbers correspond to the same numbers in figure 4. The Mascot gi number of the spots, their gene ID (PlasmoDB), their name, the theoretical MW and pI values, as well as the corresponding Mascot score are listed for MS/MS analysis (scores greater than 29 for P. falciparum and 38 for Homo sapiens are considered significant (p 0.05)). Paired average volume ratio (hyperoxic versus normoxic conditions) and p-values (t-Test) were obtained using Decyder software. MW: molecular weight.

Article Categories:
  • Research

Previous Document:  Differential endothelial cell gene expression by African Americans versus Caucasian Americans: a pos...
Next Document:  Health workforce skill mix and task shifting in low income countries: a review of recent evidence.