In vivo transcriptional profiling of Plasmodium falciparum
© Daily et al; licensee BioMed Central Ltd. 2004
Received: 01 June 2004
Accepted: 05 August 2004
Published: 05 August 2004
Both host and pathogen factors contribute to disease outcome in Plasmodium falciparum infection. The feasibility of studying the P. falciparum in vivo transcriptome to understand parasite transcriptional response while it resides in the human host is presented.
A custom made oligonucleotide array with probes based on the P. falciparum 3D7 laboratory strain chromosome 2 sequence was used to detect in vivo P. falciparum transcripts. This study analyzed transcripts from total RNA derived from small blood samples of P. falciparum infected patients and compared the in vivo expression profile to the in vitro cultivated 3D7 strain transcriptome.
The data demonstrated that in vivo transcription can be studied from a small blood sample, despite the abundance of human RNA. The in vivo transcriptome is similar to the 3D7 ring stage transcriptome, but there are significant differences in genes encoding a sexual stage antigen and surface proteins.
Whole genome transcription analysis of P. falciparum can be carried out successfully and further studies in selected patient cohorts may provide insight into parasite in vivo biology and defense against host immunity.
Plasmodium falciparum infection remains a major health problem worldwide. Its complex life cycle has hampered standard methods for the study of pathogenesis. New approaches to elucidate parasite biology using whole genomic methods have provided insight into gene function, transcriptional regulation and stage specific biology [1–4]. Characterization of the in vivo biology of this pathogen, through adaptation of a whole genome approach, would provide insight into the host-parasite relationship, parasite virulence factors and inform new strategies for intervention. Genomic scale transcriptional profiling of P. falciparum during a natural infection is presented. Small amounts of parasite RNA, isolated from a few milliliters of a blood sample are found to be sufficient for whole genome transcriptional analysis. This data show that several genes are differentially expressed in vivo, indicating differences between the transcriptional program of 3D7 laboratory strain parasites growing in culture and naturally occurring infections in the human host.
Whole genome expression has been used in studies of bacterial pathogenesis to identify genes that are specifically transcribed under in vivo conditions [5–7]. For example, genes involved in amino acid transport and metabolism are upregulated in Pasteurella multocida in vivo as compared to in vitro conditions . Similarly, analysis of P. falciparum gene expression patterns, particularly the subset of genes that are specifically expressed in the in vivo state may identify unique parasite biology when it resides in the host environment. Processes involving parasite metabolism, immune evasion and transmission may be altered in the highly specialized environment of the human host as compared to in vitro conditions. In addition, approximately 12% of P. falciparum's predicted genes have not been found to be expressed in any of the life cycle stages previously studied . Whole genomic analysis of the parasite in vivo may reveal the unique expression of such genes in vivo, providing additional targets for intervention.
This study was conducted as part of an ongoing P. falciparum chloroquine resistance study in Senegal . Patients with mild P. falciparum malaria gave consent for the study and were enrolled at an outpatient health clinic. Patients underwent venipuncture using K3 EDTA coated Vacutainers (Beckton Dickinson) and from this sample, 1.6–2.5 ml of whole blood was collected and passed through a white cell depletion filter using a 20 ml syringe. The filtered sample was centrifuged for 5 minutes at 3,200 rpm in a clinical centrifuge and placed in Tri-Reagent BD (Molecular Research Center). The samples were vortexed and stored at minus 70°C. Samples were thawed in a room temperature bath one month later and isolation of RNA was completed per Molecular Research Centers protocol. Three samples obtained in Senegal that had the highest parasitemia and largest blood volume are presented. A 14 ml blood sample from a P. falciparum infected traveler from Nigeria was similarly processed in Boston.
Oligonucleotide array analysis
Labeling and hybridization of total RNA was performed as described . Expression levels were calculated using the Match Only Integral Distribution Algorithm (MOID) . The presence or absence of gene expression was determined using methods previously described . The design of the probes to human ESTs was based on UniGene version 116.
Real time PCR
To confirm array data, a subset of genes (PFB0120, PFB0100, PFB0270 PFB0355, and PFB0065) that vary from high to low level abundance by array were quantified using real time PCR from cDNA. PFB0120: forward primer 5'-CAG CCC TCT TAG CTC TCA ACT TC-3', reverse primer 5'-AGC AAC AGC AGA GGC TAT AGA ACT-3', PFB0100: forward primer 5'-CAC CAA ATG GCT ATG CTT ATG GA-3', reverse primer 5'-TTC CAG GAG CAC CAT TAA ATC CT-3', PFB0270: forward primer 5'-ACA CTT ACT GGT ATT TCG GAA TTT-3', reverse primer 5'-TAA TTG TCC ATA TTC TTC AAT ATA T-3', PFB0355: forward primer 5'-ATT GTA AGA AAT AGT TGG GGT-3', reverse primer 5'-TAT ATC ATG CTC CTT CTT ATC A-3', PFB0065: forward primer 5'-CGT TGG TAG TGC GTT CCT TAC AA-3', reverse primer 5'-GTT CCT GCT ATA TCA GGA GCA CCA-3'. Sequence analysis confirmed the identity of the amplification products. 3D7 strain parasites were cultivated under standard conditions and synchronized with 5% sorbitol to obtain ring stage parasites for extraction of total RNA [12, 13]. cDNA was synthesized from total RNA from the Nigerian in vivo sample and 3D7 ring stage total RNA using Super Script 1st Strand synthesis system (Invitrogen). Duplicate reactions using real time PCR were performed with 1 μl cDNA with gene specific primers in 50 μl reaction volume using fluorescent dye SYBR Green (SYBR Green PCR Master Mix, Applied Biosystems). The reactions were carried out on an ABI PRISM model 7700-sequence detector and all PCR reactions amplified a single product as determined by dissociation curve analysis (Dissociation Curve Software, Applied Biosystems).
Variation between samples was assessed using Kruskal-Wallis method (non-parametric ANOVA) to test the null hypothesis. To normalize samples, the mean gene expression level was calculated for all Plasmodium genes between the 10th and 90th percentile with at least six probes. Analysis based on rank was a second method used; in each experiment the probe intensity was ranked and this resulted in equivalent quantile distribution for all probes between two experiments. This method is more conservative and will define relative rank changes between experiments and is independent of potential normalization artifacts.
Patient blood samples were collected after informed consent was obtained. The study was approved by the institutional review boards at Harvard School of Public Health, Brigham and Women's Hospital and Cheikh Anta Diop University.
Human transcripts detected in P. falciparum infected patient blood samples and in vitro cultivated 3D7 samples Comparison of expression level of the most abundant human transcripts detected in in vivo blood and in 3D7 in vitro cultivated samples. Average expression level of in vivo isolated RNA derived from the average expression level of three samples from Senegal and one from Nigeria compared with the average expression level from three 3D7 ring stage in vitro isolated RNA samples. Transcripts are listed in order of highest abundance in vivo as detected by array. Expression levels are reported as expression units (EU), using the MOID algorithm.
in vitro 3D7
haemoglobin, delta EST, similar to tctp human translationally controlled tumor protein H. sapiens
EST, Weakly similar to KIAA0902 protein H. sapiens
haemoglobin, theta 1
Bcl-2 binding component 3
ESTs, similar to B24264 proline-rich protein MP3 -M. musculus
ESTs, similar to AF1 51 859_1 CGI-101 protein H. sapiens
haemoglobin, epsilon 1
hypothetical protein FLJ11137
peroxisomal short-chain alcohol dehydrogenase
ESTs, similar to spliceosomal protein SAP 155 H. sapiens Homo sapiens mRNA; cDNA DKFZp434H0820 (from clone DKFZp434H0820); partial cds
ESTs, similar to A48018 mucin 7 precursor, salivary – H. sapiens
hypothetical protein FLJ20124 ATPase, H+ transporting, lysosomal (vacuolar proton pump), member D
retinol dehydrogenase 5 (11-cisand 9-cis)
P. falciparum genes expressed in vivo encoded by chromosome 2. In vivo transcripts from the Nigerian sample were defined as present as compared to uninfected blood control hybridization. Asterisk (*) denotes transcripts that were also detected in a Senegal derived blood sample. Genes in bold are uniquely expressed in vivo and were not found to be expressed in the previously reported 3D7 ring stage transcriptome. Gene locus is from PlasmoDB 4.1 http://plasmodb.org/.
prt of the MAK16 family
transmission blocking target antigen
prt with 5'-3' exonucl. domain
adenylosuccinate lyase (OO)
phosphatase (acid phosphatase family)
ribosome releasing factor (OO, TP)
phospholipase A2-like a/b fold hydrolase
elF-4A-like DEAD family RNA helicase
SERA antigen/papain-like protease
GAF domain prt
SERA antigen/papain-like protease
asparaginyl-tRNA synthetase (OO, TP)
SERA antigen/papain-like protease
Leu/Phe-tRNA prt transferase
SERA antigen/papain-like protease
prt with DnaJ domain, DNJ1/SIS1 family
SERA antigen/papain-like protease
membrane protein, PfEMP3
ATP synthase alpha chain
calcium-dept. prt kinase
Pf EMP1 fragment
chromatin-binding prt (SKI/SNW family)
PfEMP1 (var gene)
prt with DnaJ domain (RESA-like)
knob-associated His-rich prt
merozoite surface antigen MSP-2
predicted multiple-TM membrane prt
predicted multiple-TM membrane prt
predicted membrane associated prt
predicted integral membrane protein
predicted secreted protein
predicted secreted protein
This data confirms that in vivo whole genomic expression can be performed despite the potential technical challenges of scarce RNA contained in a small blood volume sample and presence of abundant human RNA. Tri-Reagent BD was used to stabilize RNA and the samples were stored at -80°C before transport to the US. The denaturing gel suggests that the RNA remains intact using this reagent. Notably Kyes et al. reported that minus 80°C rather than 4°C or minus 20°C is the optimal temperature to store field sample RNA for detection of long transcripts such as var . Surprisingly, abundant human RNA was detected in the denaturing gel despite buffy coat depletion of the samples. In addition, the degree of hybridization to human probes on the oligonucleotide array used here was not seen in the previously studied in vitro sample . Haemoglobin is found to be the most abundantly expressed human transcript (Table 1) with considerably higher levels noted in the in vivo samples as compared to the 3D7 in vitro sample. Although human red cells are used for culturing in the in vitro system, reticulocytes which contribute to the haemoglobin detected are not abundant in in vitro samples. This is most likely due to the observation that reticulocyte levels display a 75% loss at 48 hours, when placed at 37°C, which is the condition of in vitro culture . In addition, other human RNA such as ribosomal RNA is more abundant in the in vivo samples and may be secondary to white cell contamination. Cross hybridization of human RNA to parasite probes may occur. However, Zhou et al have shown that human genomic DNA does not highly cross hybridize to the parasite probes on this custom array . The high specificity of this array is due to the nature of the highly AT rich parasite genome as compared to the human genome and the careful selection of parasite unique 25 mers probes [2, 18, 19]. Due to this high specificity, it is likely that buffy coat depletion is not necessary for analysis of in vivo parasite transcripts when using these probes.
The amount of blood volume necessary for comprehensive detection of transcription depends on level of parasitemia and method of microarray analysis. This analysis utilized the Affymetrix system which requires very little starting RNA as compared to other methods . Due to the presence of human RNA in the in vivo samples it was not possible to determine how much parasite RNA is required for whole genome analysis. The samples from Senegal were of low parasitemia and small volumes, whereas the 14 ml blood sample with 0.4% parasitemia was sufficient to detect a greater number of chromosome 2 transcripts. Four to five mls of packed blood in a patient with a greater than 2% parasitemia should provide sufficient material for whole genome analysis using these methods. The data demonstrates the reproducibility of the method by independent hybridizations and that maximal senstivity can be achieved with up to 60 μg of cRNA, using this array.
PFB0120w is the most abundant in vivo transcript in all samples encoded on chromosome 2. This gene is a member of a recently described gene family, etramps, expressed at early ring stage encoding a protein thought to be involved in erythrocyte remodeling; this was also the most highly expressed transcript in in vitro ring stage cultures [2, 21, 22]. The in vivo expressed genes from the Nigerian sample was compared to that of the in vitro 3D7 ring stage chromosome 2 transcriptome . As expected, there was a good correlation between the in vivo ring stage transcriptome and the 3D7 ring stage transcriptome. Most of the genes that are expressed in vivo are also expressed in vitro particularly those involved in cellular function and genes encoding hypothetical proteins (Table 2) . A number of differentially expressed genes involved in transmission and antigenic variation were identified. There is high transcription level of transmission blocking target antigen (PFB0405w) in the in vivo samples. Previously this gene has been demonstrated to be expressed and transcribed only in the sexual transmission stage [3, 9]. Since no transmissible forms were identified by microscopy this suggests that the in vivo samples may have more parasites that are undergoing or are committed to sexual development than is detectable by microscopy. Several genes encoding membrane proteins, including SERA and rifin genes, were also found to be differentially expressed. These genes are members of multigene family that encode surface proteins which are thought to be involved in immune evasion [23, 24]. The observed increase in transcription for these genes in vivo could be due either to genetic or transcriptional modulation of the parasite's defense repertoire or geographic variation. Differences between the in vivo and in vitro expression of SERA antigens may be due to differences in the in vivo stage of development as this gene family is transcribed at later stages in the in vitro life cycle [2, 4, 25]. There was overall higher hybridization intensity of the in vitro samples due to higher parasite counts and a subset of genes were found to be overexpressed in vitro after normalization. This analysis focussed on genes overexpressed in vivo as these results would not be influenced by normalization algorithms. Overall, this data demonstrates that P. falciparum transcripts can be detected from in vivo samples, and that there are potentially important differences between trancription of in vivo samples and that of the 3D7 in vitro trancription profile.
In summary, this study provides evidence that whole genome gene expression in P. falciparum can be studied in vivo from a small blood sample of an infected patient. The in vivo sample however contains human RNA, whose quantity may vary from sample to sample and therefore differenes in parasite transcript level between samples must be reported relative to a reference transcript. Despite the abundance of human RNA the genomes are sufficiently different with resultant probes specificity. Predictably, there was a high correlation of in vivo expression with the in vitro ring stage 3D7 transcriptome . Importantly these data also suggest differences between in vivo and in vitro expression levels in genes typically found in transmissable forms and encoding variant surface proteins. Evaluation of trancription of genes specific for gametocyte development in specific patient populations may uncover the in vivo conditions that favor development of transmissable forms. Similarily, a whole genome analysis can comprehensively characterize expression of multigene families that encode variant surface proteins under in vivo conditions. Further exploration of the in vivo biology of P. falciparum using specific probes to all annotated genes will be undertaken to confirm and explore other important biological differences. This new approach will further the understanding of the host-pathogen interaction and may result in the development of new strategies to combat this disease.
- Ben Mamoun C, Gluzman IY, Hott C, MacMillan SK, Amarakone AS, Anderson DL, Carlton JM, Dame JB, Chakrabarti D, Martin Brownstein B, Goldberg D: Co-ordinated programme of gene expression during asexual intraerythrocytic development of the human malaria parasite Plasmodium falciparum revealed by microarray analysis. Mol Microbiol. 2001, 39: 26-36. 10.1046/j.1365-2958.2001.02222.x.View ArticlePubMedGoogle Scholar
- Le Roch KG, Zhou Y, Batalov S, Winzeler EA: Monitoring the chromosome 2 intraerythrocytic transcriptome of Plasmodium falciparum using oligonucleotide arrays. Am J Trop Med Hyg. 2002, 67: 233-243.PubMedGoogle Scholar
- Florens L, Washburn MP, Raine JD, Anthony RM, Grainger M, Haynes JD, Moch JK, Muster N, Sacci JB, Tabb DL, Witney AA, Wolters D, Wu Y, Gardner MJ, Holder AA, Sinden RE, Yates JR, Carucci DJ: A proteomic view of the Plasmodium falciparum life cycle. Nature. 2002, 419: 520-526. 10.1038/nature01107.View ArticlePubMedGoogle Scholar
- 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: R9-10.1186/gb-2003-4-2-r9.PubMed CentralView ArticlePubMedGoogle Scholar
- Heithoff DM, Conner CP, Hentschel U, Govantes F, Hanna PC, Mahan MJ: Coordinate intracellular expression of Salmonella genes induced during infection. J Bacteriol. 1999, 181: 799-807.PubMed CentralPubMedGoogle Scholar
- Chiang SL, Mekalanos JJ, Holden DW: In vivo genetic analysis of bacterial virulence. Annu Rev Microbiol. 1999, 53: 129-154. 10.1146/annurev.micro.53.1.129.View ArticlePubMedGoogle Scholar
- Basso H, Rharbaoui F, Staendner LH, Medina E, Garcia-Del Portillo F, Guzman CA: Characterization of a novel intracellularly activated gene from Salmonella enterica serovar typhi. Infect Immun. 2002, 70: 5404-5411. 10.1128/IAI.70.10.5404-5411.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Boyce JD, Wilkie I, Harper M, Paustian ML, Kapur V, Adler B: Genomic scale analysis of Pasteurella multocida gene expression during growth within the natural chicken host. Infect Immun. 2002, 70: 6871-6879. 10.1128/IAI.70.12.6871-6879.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Le Roch KG, Zhou Y, Blair PL, Grainger M, Moch JK, Haynes JD, De La Vega P, Holder AA, Batalov S, Carucci DJ, Winzeler EA: Discovery of gene function by expression profiling of the malaria parasite life cycle. Science. 2003, 301: 1503-1508. 10.1126/science.1087025.View ArticlePubMedGoogle Scholar
- Thomas SM, Ndir O, Dieng T, Mboup S, Wypij D, Maguire JH, Wirth DF: In vitro chloroquine susceptibility and PCR analysis of pfcrt and pfmdr1 polymorphisms in Plasmodium falciparum isolates from Senegal. Am J Trop Med Hyg. 2002, 66: 474-480.PubMedGoogle Scholar
- Zhou Y, Abagyan R: Match-Only Integral Distribution (MOID) Algorithm for high-density oligonucleotide array analysis. BMC Bioinformatics. 2002, 3: 3-10.1186/1471-2105-3-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Trager W, Jensen JB: Human malaria parasites in continuous culture. Science. 1976, 193: 673-675.View ArticlePubMedGoogle Scholar
- Lambros C, Vanderberg JP: Synchronization of Plasmodium falciparum erythrocytic stages in culture. J Parasitol. 1979, 65: 418-420.View ArticlePubMedGoogle Scholar
- Snounou G, Zhu X, Siripoon N, Jarra W, Thaithong S, Brown KN, Viriyakosol S: Biased distribution of msp1 and msp2 allelic variants in Plasmodium falciparum populations in Thailand. Trans R Soc Trop Med Hyg. 1999, 93: 369-374. 10.1016/S0035-9203(99)90120-7.View ArticlePubMedGoogle Scholar
- Kyes S, Pinches R, Newbold C: A simple RNA analysis method shows var and rif multigene family expression patterns in Plasmodium falciparum. Mol Biochem Parasitol. 2000, 105: 311-315. 10.1016/S0166-6851(99)00193-0.View ArticlePubMedGoogle Scholar
- Hansson GK, Andersson M, Jarl H, Stemme S: Flow cytometric analysis of reticulocytes using an RNA-binding fluorochrome. Scand J Clin Lab Invest. 1992, 52: 35-41.View ArticlePubMedGoogle Scholar
- Zhou Y, Abagyan R: Algorithms for high-density oligonucleotide array. Curr Opin Drug Discov Devel. 2003, 6: 339-345.PubMedGoogle Scholar
- Gardner MJ, Tettelin H, Carucci DJ, Cummings LM, Aravind L, Koonin EV, Shallom S, Mason T, Yu K, Fujii C, Pederson J, Shen K, Jing J, Aston C, Lai Z, Schwartz DC, Pertea M, Salzberg S, Zhou L, Sutton GG, Clayton R, White O, Smith HO, Fraser CM, Hoffman SL: Chromosome 2 sequence of the human malaria parasite Plasmodium falciparum. Science. 1998, 282: 1126-1132. 10.1126/science.282.5391.1126.View ArticlePubMedGoogle Scholar
- Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, International Human Genome Sequencing Consortium: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.View ArticlePubMedGoogle Scholar
- Tefferi A, Bolander ME, Ansell SM, Wieben ED, Spelsberg TC: Primer on medical genomics. Part III: Microarray experiments and data analysis. Mayo Clin Proc. 2002, 77: 927-940.View ArticlePubMedGoogle Scholar
- Spielmann T, Fergusen DJ, Beck HP: Etramps, a new Plasmodium falciparum gene family coding for developmentally regulated and highly charged membrane proteins located at the parasite-host cell interface. Mol Biol Cell. 2003, 14 (4): 1529-1544. 10.1091/mbc.E02-04-0240.PubMed CentralView ArticlePubMedGoogle Scholar
- Blisnick T, Morales Betoulle ME, Barale JC, Uzureau P, Berry L, Desroses S, Fujioka H, Mattei D, Braun Breton C: Pfsbp1, a Maurer's cleft Plasmodium falciparum protein, is associated with the erythrocyte skeleton. Mol Biochem Parasitol. 2000, 111 (1): 107-121. 10.1016/S0166-6851(00)00301-7.View ArticlePubMedGoogle Scholar
- Miller SK, Good RT, Drew DR, Delorenzi M, Sanders PR, Hodder AN, Speed TP, Cowman AF, de Koning-Ward TF, Crabb BS: A subset of Plasmodium falciparum SERA genes are expressed and appear to play an important role in the erythrocytic cycle. J Biol Chem. 2002, 277: 47524-47532. 10.1074/jbc.M206974200.View ArticlePubMedGoogle Scholar
- Abdel-Latif MS, Khattab A, Lindenthal C, Kremsner PG, Klinkert MQ: Recognition of variant Rifin antigens by human antibodies induced during natural Plasmodium falciparum infections. Infect Immun. 2002, 70: 7013-7021. 10.1128/IAI.70.12.7013-7021.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Aoki S, Li J, Itagaki S, Okech BA, Egwang TG, Matsuoka H, Palacpac NM, Mitamura T, Horii T: Serine repeat antigen (SERA5) is predominantly expressed among the SERA multigene family of Plasmodium falciparum, and the acquired antibody titers correlate with serum inhibition of the parasite growth. J Biol Chem. 2002, 277: 47533-47540. 10.1074/jbc.M207145200.View ArticlePubMedGoogle Scholar
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.