Methodology | Open | Published:
New Agilent platform DNA microarrays for transcriptome analysis of Plasmodium falciparum and Plasmodium berghei for the malaria research community
Malaria Journalvolume 11, Article number: 187 (2012)
DNA microarrays have been a valuable tool in malaria research for over a decade but remain in limited use in part due their relatively high cost, poor availability, and technical difficulty. With the aim of alleviating some of these factors next-generation DNA microarrays for genome-wide transcriptome analysis for both Plasmodium falciparum and Plasmodium berghei using the Agilent 8x15K platform were designed.
Probe design was adapted from previously published methods and based on the most current transcript predictions available at the time for P. falciparum or P. berghei. Array performance and transcriptome analysis was determined using dye-coupled, aminoallyl-labelled cDNA and streamlined methods for hybridization, washing, and array analysis were developed.
The new array design marks a notable improvement in the number of transcripts covered and average number of probes per transcript. Array performance was excellent across a wide range of transcript abundance, with low inter-array and inter-probe variability for relative abundance measurements and it recapitulated previously observed transcriptional patterns. Additionally, improvements in sensitivity permitted a 20-fold reduction in necessary starting RNA amounts, further reducing experimental costs and widening the range of application.
DNA microarrays utilizing the Agilent 8x15K platform for genome-wide transcript analysis in P. falciparum and P. berghei mark an improvement in coverage and sensitivity, increased availability to the research community, and simplification of the experimental methods.
Malaria remains one of the most deadly infections to humans worldwide. It is caused by several species in the genus Plasmodium, with Plasmodium falciparum being the most virulent. The malaria parasite’s complex lifecycle is initiated by an infectious bite from a female anopheline mosquito, injecting sporozoites into the bloodstream and leading to the invasion of liver cells. The subsequent growth of asexual parasites within red blood cells (RBCs) is responsible for pathology. A small percentage of infected RBCs transforms into male and female gametocytes, which can be transmitted to the mosquito vector. After fusion of gametes within the midgut, development proceeds at various sites in the mosquito ending with the invasion of the salivary glands. Upon the vector’s next bloodmeal the life cycle is completed with the injection of salivary gland sporozoites.
While malaria is among the longest studied afflictions of humans, progress towards new therapeutics and vaccines has generally been slow. The past decade has seen significant advances in the fundamental understanding of the parasite’s biology, which in turn has opened new and promising avenues for novel anti-malarial development. This resurgence has been brought about in large part by technological advances that have enhanced the ability to genetically manipulate the parasite [1, 2] as well as through the insights provided by the whole genome sequencing of many species and strains of Plasmodium[3–6]. A major benefit of genome sequencing has been the ability to assay gene transcription on a genome-wide basis through the use of DNA microarrays, which have allowed probing questions to be asked regarding the transcriptional status at distinct life stages, differences in gene expression between strains, in response to changes in environmental conditions or drug perturbations. Microarrays can also be used to find genetic alterations such as copy number variations or single nucleotide polymorphisms. Despite these possible uses, the DNA microarray has not been broadly utilized for several reasons including the large amount of biological material required, technological impracticality for many research settings, and relatively high cost.
The first DNA microarrays for P. falciparum were generated using either sheared genomic DNA  or cDNAs  and were used to compare differences across asexual blood stages. In the absence of the genome sequence, follow up of potential gene expression differences was pursued on a candidate basis by sequencing of the material spotted on the array. With the completion of the P. falciparum strain 3D7 reference genome , oligonucleotide-based arrays quickly followed and provided an unprecedented in-depth view of transcriptional changes during asexual development in red blood cells [9, 10]. Subsequent studies determined that other strains (HB3, Dd2) exhibited similar transcriptional programs in vitro and ultimately these general observations were extended to transcriptional profiles of patient isolates by either profiling parasite mRNA abundance directly from blood [12, 13] or after short-term ex vivo culturing . Such transcriptional signatures from non-culture adapted parasites have identified subtle and important differences that are the source of on-going research. Other in vitro studies have struggled to identify variations in transcriptional profiles under environmental perturbations [15–17], although larger scale efforts examining growth perturbations under dozens of conditions have yielded associations between genes based on transcriptional covariation [18–21]. More detailed studies have examined transcriptional changes associated with differences in cellular adhesion and antigenic variation [15, 22]. Ultimately, one of the real powers of DNA microarray analysis will be to characterize transcriptional responses to targeted genetic alterations [23, 24].
Since these early studies, other stages of P. falciparum development have also been explored using DNA microarrays including sporozoites  and various stages of gametocyte maturation [23, 26, 27]. Furthermore, several efforts have utilized DNA microarrays for genome-wide transcriptional analysis in the rodent models of malaria, Plasmodium berghei and Plasmodium yoelii. From these efforts, there are reports characterizing ookinete development in P. berghei[28, 29] and Mikolajczek et al have examined oocyst vs. salivary gland sporozoites in P. yoelii. The liver stage of human malaria parasites has been difficult to analyse transcriptionally. However, microarray studies of P. yoelii, or P. berghei liver stage infections have characterized both the parasite  and host transcriptional programs during this developmental stage .
In light of the many insights already gained from this relatively small number of genome-wide experiments, improving the performance and increasing the availability of Plasmodium DNA microarray to the research community remains a worthwhile effort.
Annotated transcripts where downloaded for P. falciparum strain 3D7 (PlasmoDB release 7.1)  and P. berghei ANKA (Sanger Institute, March 2011). These sequences were supplemented with exogenous genes commonly used in Plasmodium molecular biology (hdhfr, bsd, gfp, fkbp, etc.) and published P. falciparum non-protein coding RNAs not included in PlasmoDB at the time [34–38]. Transcripts were first masked for simple repeats using RepeatMasker  and then split in half until all fragments were less than 1040nt in size to obtain approximately 15,000 sub-sequences for which a single probe would be designed to achieve near-maximal coverage given the 8x15K array format. These sequences were analysed using OligoBLAST  to design 60mer oligonucleotide probes followed by probe selection using OligoRankPick  with the following parameters: WBLAST= 6,7,8,9 / WGC = 6,7,8,9 / WSW:1,2 / WLZ = 1 / length = 60 / GC% = 31.4 (see  for details). Probes were then BLASTed against the Crick strand of the target transcripts using NCBI BLAST  and probes with multiple hits of E-value < 1e-15 were marked as non-unique (Additional files 1 and 2). Probes were then randomly distributed across the array (to reduce spatial bias) using the 8x15K format after submission to Agilent Technologies eArray . For this design, no linker sequence was used and any empty features were filled with random duplicate probes. Slides with printed arrays were ordered directly from Agilent Technologies (Santa Clara, CA, USA).
Parasite growth and culturing
Plasmodium falciparum strain 3D7 was cultivated in vitro using standard techniques  with the following modifications. Parasites were grown at 4–5 % haematocrit at 5 % CO2/ 6 % O2 in RPMI media supplemented with 0.25 % Albumax I (Gibco, San Diego, CA, USA), 0.1 mM hypoxanthine, 2 g/L sodium bicarbonate, 25 mM HEPES pH7.4, and 50 μg/L gentamycin. For assessing transcripts throughout the intra-erythrocytic development, parasites were double synchronized by L-alanine treatment  12 hours apart during the preceding cycle and harvested every six hours throughout the 48-hour cycle, starting approximately three hours post-invasion.
Cells were harvested from 50 ml of culture suspension by centrifugation at 1,500 rpm for 5 min. Total RNA was extracted and purified using TriZol reagent (Invitrogen, Grand Island, NY, USA) at a volume of 5 ml reagent to 1 ml of packed parasitized erythrocytes as previously described . Quality and quantity of total RNA extracted was assessed by agarose gel electrophoresis and using a ND-1000 (NanoDrop Technologies, Thermo Scientific, Wilmington, DE, USA) then stored at −80 °C for later use.
cDNA generation and dye coupling
Detailed step-by-step protocols can be found in the latest edition of Methods in Malaria Research Chapter 2.6 (in preparation) and on the Llinás laboratory website . Single-strand aminoallyl-containing cDNA synthesis and Amersham CyDye-coupling (GE Healthcare, Piscataway, NJ, USA) was carried out as previously described . To eliminate Cy5 degradation by ozone , all steps starting with dye resuspension were carried out in an ozone-free environment.
Array hybridization and washing
Final cDNA concentration and dye-incorporation was assessed on a NanoDrop ND-1000 spectrophotometer. These measurements were utilized to calculate the frequency of incorporation and samples that met the manufacturer’s suggested target of 10–20 fmol of dye/ng cDNA were further used for array hybridization. Equal amounts between 50 and 1000 ng of Cy3 and Cy5-labelled cDNA were hybridized on each array for 16 h in a rotating hybridization oven (10 rpm) at 65 °C. Prior to scanning, arrays were washed in 6X and 0.06X SSPE (both containing 0.005 % N-lauryl-sarcosine (Sigma-Aldrich, St. Louis, MO, USA), followed by an acetonitrile rinse. For the 48-hour time-course, cDNA from each time-point was labelled with Cy5 and hybridized to an equal amount of Cy3-labelled cDNA reference pool generated from equal amounts of ring, trophozoite, and schizont stage mRNA.
Scanning, data acquisition and analysis
Arrays were scanned on an Agilent G2505B Microarray Scanner (Agilent Technologies, Santa Clara, CA, USA) with 5 μm resolution at wavelengths of 532 nm (Cy3) and 633 nm (Cy5) using the extended dynamic range (10–100 %) setting. Normalized intensities were extracted using Agilent Feature Extractor Software Version 9.5 employing the GE2-v5_95_Feb07_no_spikein extraction protocol and uploaded to the Princeton University Microarray Database (PUMA.princeton.edu) for analysis.
Results & discussion
Three main goals provided the impetus for switching DNA microarray platforms from an in-house spotted 70mer array to the Agilent Technologies 60mer SurePrint platform; (1) to improve array performance in terms of the number of features, sensitivity, and inter-array variability, (2) to have the ability to easily and affordably update microarray design as genome annotation improves, and (3) to make Plasmodium transcriptome analysis widely available to the research community at-large.
Switching to the Agilent Technologies platform enables researchers with access to a two-color microarray scanner to perform genome-wide transcriptome analysis for P. falciparum or P. berghei without the need for collaboration or the considerable cost of synthesizing large number of oligonucleotides for spotting and the even greater cost of microarray spotting equipment. Due to Agilent’s manufacturing process any number of slides, from one to hundreds, can be printed and ordered at identical cost per slide. Time from order to delivery varied from three to six weeks.
Unlike the photolithographic manufacturing of Affymetrix arrays, which require expensive masks to print, or spotted oligo-arrays, which required the one-time bulk purchase of synthesized oligonucleotides for a large number of arrays with each redesign ( and ), Agilent arrays can be redesigned without any additional cost, allowing researchers to alter or update array designs quickly and affordably as the need arises.
Agilent Technologies offers a variety of slide designs with one or more arrays per slide ranging from 15,000 to 1 million features per array. Given the good performance of the previous spotted 10,000-feature array and to keep cost at a minimum, the 8x15K slide format was chosen, which features eight arrays with 15,000+ features on each slide. This provides the added advantage of being able to hybridize multiple samples on the same slide, minimizing variation of hybridization conditions thereby reducing array-to-array variability. In addition, the chosen format is optimal for a full 48-hour intraerythrocytic development cycle (IDC) to be assayed on a single slide when collecting samples every six hours.
For this array design, the annotated transcripts from P. falciparum strain 3D7 (PlasmoDB  release 7.1) and P. berghei strain ANKA (Welcome Trust Sanger Centre, March 2010 Annotation) were used. These sequences were supplemented with exogenous genes commonly used in Plasmodium molecular biology (hdhfr, bsd, gfp, fkbp, etc.) (Additional files 1 and 2) and for P. falciparum published non-protein coding RNAs not included in the PlasmoDB release [34–38]. These transcripts were masked for simple repeats using RepeatMasker  and split depending on sequence length to yield ~15,000 sub-sequences, for each of which a single probe would be designed. Transcripts of fewer than 1040 nt were not split, 1041–2080 nt were split into two sub-sequences, 2081–3120 nt were split into three and so forth, with the goal of just under 15,000 probes and 1+ probes for every ~1 kb of transcript length.
These sub-sequences were run through OligoBLAST  to generate all possible non-masked 60mers (maximum length for Agilent arrays) and score their uniqueness. From this set the best quality probe was chosen for each sub-sequence using OligoRankPick . The resulting probes were BLASTed against the Crick strand of the initial target transcripts and any probes with hits of E-value < 1e-15 were marked as non-unique to indicate the potential for cross-hybridization. Furthermore, several additional handpicked probes for targets of particular interest to this laboratory were also included (e.g. DNA-binding domains of ApiAP2 proteins). This final set of probes was randomly distributed across the 8x15K array (one copy per array) format with the remaining available features filled by random repeats of existing probes. The 8x15K design includes an Agilent Control Grid , which serves various functions including manufacturing controls, alignment during feature extraction, and features for spike-in experiments using the Agilent RNA Spike-In Kit. Additionally, a small number of transcripts were repeated 15+ times each across the array for use in further spatial normalization during feature extraction.
The final designs of the new 8x15K P. falciparum and P. berghei arrays are summarized in Table 1 and details on the 8x15K arrays can be found in Additional files 1 and 2. Briefly, 14,000+ probes were designed for each array, with an average of 2.5+ probes per transcript (Table 1, Additional file 3). Probe specificity was estimated at greater than 97 % for both arrays and over 97 % of transcripts have one or more unique probes. Less than 1 % of transcripts are not represented on these arrays either because they were too short (< 60 nt) or failed to contain a single stretch of 60 nt of non-repeat sequence. Additionally, roughly 2 % of transcripts did not have any unique probes for one of the following reasons: very short (ncRNAs), recent gene duplication (e.g. falcipain-2a/b (PF11_0161/PF11_0164) or elongation factor 1a (PF13_0304/PF13_0305)), or belonging to large gene families (vars, rifins, BIRs, etc.). It should be noted that for most genes belonging to such a gene family this does not indicate an absence of unique probes. For example, the P. falciparum array contains unique probes for 58 of the 63 var genes (var gene fragments or var pseudogenes were not included in this number but were included when analysing probes for specificity). In total there are unique probes for 5,752 transcripts including 930 transcripts not previously represented on the spotted array design.
Performance analysis was limited to the P. falciparum array as large quantities of P. falciparum total RNA were easily accessible. The P. berghei 8x15K array is not expected to differ in terms of array-wide performance measurements as probes were designed using the same methodology. In fact, collaborators have successfully used these arrays for transcript analysis (Andrew Waters, personal communication) validating their utility using similar methods.
A cDNA pool derived from an even mixture of ring, trophozoite, and schizont-stage RNA was self-hybridized onto three P. falciparum 8x15K arrays to assess inter-probe and inter-array variability. For further in-depth analysis, six genes with coverage by at least 15 probes each were chosen as they represent a wide range of transcript abundance encompassing the full dynamic range of the measurements. Inter-array consistency was good overall for both the mean intensity of each gene (σ ≤ 15 %, Figure 1A top) or each probe individually (σ ≤ 23 %, Figure 1B top). While probe-to-probe intensity can vary significantly across a given gene (Figure 1B top), but has little effect on log-ratio measurements when averaged across probes (Figure 1A bottom) or even when considering individual probes (Figure 1B bottom). When extending this analysis to the entire array, technical reproducibility was good between replicates (Pearson’s r > 0.80) at both the gene and probe level for the log2(Cy5/Cy3) ratio as well as for the Cy5 and Cy3 intensities across all probes (Additional file 4). For this reason a common cDNA reference pool should be used when doing multi-sample comparisons of transcript abundance. Furthermore, even for a large transcript such as MAL8P1.139 (17.8 kb), no 3’ bias (Pearson’s r = − 0.36, Figure 1A) was found during cDNA generation.
This paper describes two new long-oligonucleotide arrays for transcriptome analysis in P. falciparum and P. berghei using the Agilent Technologies SurePrint 60mer 8x15K platform. The advantages of these arrays include 50 % more probes compared to our previous in-house spotted-array, the ability to easily update the array design as genome annotation evolves at no cost, and wide availability to the research community. The new 8x15K P. falciparum array was demonstrated to have excellent reproducibility and sensitivity, which allows for transcript analysis with considerably lower amounts of materials compared to previously used methods, thus minimizing the need for high volume samples or RNA/cDNA amplification.
Limenitakis J, Soldati-Favre D: Functional genetics in Apicomplexa: Potentials and limits. FEBS Lett. 2011, 585: 1579-1588. 10.1016/j.febslet.2011.05.002.
Carvalho TG, Menard R: Manipulating the Plasmodium genome. Curr Issues Mol Biol. 2005, 7: 39-55.
Gardner MJ, Hall N, Fung E, White O, Berriman M, Hyman RW, Carlton JM, Pain A, Nelson KE, Bowman S, Paulsen IT, James K, Eisen JA, Rutherford K, Salzberg SL, Craig A, Kyes S, Chan MS, Nene V, Shallom SJ, Suh B, Peterson J, Angiuoli S, Pertea M, Allen J, Selengut J, Haft D, Mather MW, Vaidya AB, Martin DM, Fairlamb AH, Fraunholz MJ, Roos DS, Ralph SA, McFadden GI, Cummings LM, Subramanian GM, Mungall C, Venter JC, Carucci DJ, Hoffman SL, Newbold C, Davis RW, Fraser CM, Barrell B: Genome sequence of the human malaria parasitePlasmodium falciparum. Nature. 2002, 419: 498-511. 10.1038/nature01097.
Carlton JM, Angiuoli SV, Suh BB, Kooij TW, Pertea M, Silva JC, Ermolaeva MD, Allen JE, Selengut JD, Koo HL, Peterson JD, Pop M, Kosack DS, Shumway MF, Bidwell SL, Shallom SJ, van Aken SE, Riedmuller SB, Feldblyum TV, Cho JK, Quackenbush J, Sedegah M, Shoaibi A, Cummings LM, Florens L, Yates JR, Raine JD, Sinden RE, Harris MA, Cunningham DA, Preiser PR, Bergman LW, Vaidya AB, van Lin LH, Janse CJ, Waters AP, Smith HO, White OR, Salzberg SL, Venter JC, Fraser CM, Hoffman SL, Gardner MJ, Carucci DJ: Genome sequence and comparative analysis of the model rodent malaria parasitePlasmodium yoelii yoelii. Nature. 2002, 419: 512-519. 10.1038/nature01099.
Carlton J: ThePlasmodium vivaxgenome sequencing project. Trends Parasitol. 2003, 19: 227-231. 10.1016/S1471-4922(03)00066-7.
Pain A, Bohme U, Berry AE, Mungall K, Finn RD, Jackson AP, Mourier T, Mistry J, Pasini EM, Aslett MA, Balasubrammaniam S, Borgwardt K, Brooks K, Carret C, Carver TJ, Cherevach I, Chillingworth T, Clark TG, Galinski MR, Hall N, Harper D, Harris D, Hauser H, Ivens A, Janssen CS, Keane T, Larke N, Lapp S, Marti M, Moule S, Meyer IM, Ormond D, Peters N, Sanders M, Sanders S, Sargeant TJ, Simmonds M, Smith F, Squares R, Thurston S, Tivey AR, Walker D, White B, Zuiderwijk E, Churcher C, Quail MA, Cowman AF, Turner CM, Rajandream MA, Kocken CH, Thomas AW, Newbold CI, Barrell BG, Berriman M: The genome of the simian and human malaria parasitePlasmodium knowlesi. Nature. 2008, 455: 799-803. 10.1038/nature07306.
Hayward RE, Derisi JL, Alfadhli S, Kaslow DC, Brown PO, Rathod PK: Shotgun DNA microarrays and stage-specific gene expression inPlasmodium falciparummalaria. Mol Microbiol. 2000, 35: 6-14. 10.1046/j.1365-2958.2000.01730.x.
Ben Mamoun C, Gluzman IY, Hott C, MacMillan SK, Amarakone AS, Anderson DL, Carlton JM, Dame JB, Chakrabarti D, Martin RK, Brownstein BH, Goldberg DE: Co-ordinated programme of gene expression during asexual intraerythrocytic development of the human malaria parasitePlasmodium falciparumrevealed by microarray analysis. Mol Microbiol. 2001, 39: 26-36. 10.1046/j.1365-2958.2001.02222.x.
Bozdech Z, Llinas M, Pulliam BL, Wong ED, Zhu J, DeRisi JL: The transcriptome of the intraerythrocytic developmental cycle ofPlasmodium falciparum. PLoS Biol. 2003, 1: E5-
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.
Llinas M, Bozdech Z, Wong ED, Adai AT, DeRisi JL: Comparative whole genome transcriptome analysis of threePlasmodium falciparumstrains. Nucleic Acids Res. 2006, 34: 1166-1173. 10.1093/nar/gkj517.
Siau A, Toure FS, Ouwe-Missi-Oukem-Boyer O, Ciceron L, Mahmoudi N, Vaquero C, Froissard P, Bisvigou U, Bisser S, Coppee JY, Bischoff E, David PH, Mazier D: Whole-transcriptome analysis ofPlasmodium falciparumfield isolates: identification of new pathogenicity factors. J Inf Dis. 2007, 196: 1603-1612. 10.1086/522012.
Daily JP, Scanfeld D, Pochet N, Le Roch K, Plouffe D, Kamal M, Sarr O, Mboup S, Ndir O, Wypij D, Levasseur K, Thomas E, Tamayo P, Dong C, Zhou Y, Lander ES, Ndiaye D, Wirth D, Winzeler EA, Mesirov JP, Regev A: Distinct physiological states ofPlasmodium falciparumin malaria-infected patients. Nature. 2007, 450: 1091-1095. 10.1038/nature06311.
Mackinnon MJ, Li J, Mok S, Kortok MM, Marsh K, Preiser PR, Bozdech Z: Comparative transcriptional and genomic analysis ofPlasmodium falciparumfield isolates. PLoS Pathog. 2009, 5: e1000644-10.1371/journal.ppat.1000644.
Ralph SA, Bischoff E, Mattei D, Sismeiro O, Dillies MA, Guigon G, Coppee JY, David PH, Scherf A: Transcriptome analysis of antigenic variation in Plasmodium falciparum–var silencing is not dependent on antisense RNA. Genome Biol. 2005, 6: R93-10.1186/gb-2005-6-11-r93.
Le Roch KG, Johnson JR, Ahiboh H, Chung DW, Prudhomme J, Plouffe D, Henson K, Zhou Y, Witola W, Yates JR, Mamoun CB, Winzeler EA, Vial H: A systematic approach to understand the mechanism of action of the bisthiazolium compound T4 on the human malaria parasite. Plasmodium falciparum. BMC Genomics. 2008, 9: 513-10.1186/1471-2164-9-513.
Ganesan K, Ponmee N, Jiang L, Fowble JW, White J, Kamchonwongpaisan S, Yuthavong Y, Wilairat P, Rathod PK: A genetically hard-wired metabolic transcriptome inPlasmodium falciparumfails to mount protective responses to lethal antifolates. PLoS Pathog. 2008, 4: e1000214-10.1371/journal.ppat.1000214.
Hu G, Cabrera A, Kono M, Mok S, Chaal BK, Haase S, Engelberg K, Cheemadan S, Spielmann T, Preiser PR, Gilberger TW, Bozdech Z: Transcriptional profiling of growth perturbations of the human malaria parasitePlasmodium falciparum. Nat Biotechnol. 2010, 28: 91-98. 10.1038/nbt.1597.
Chaal BK, Gupta AP, Wastuwidyaningtyas BD, Luah YH, Bozdech Z: Histone deacetylases play a major role in the transcriptional regulation of thePlasmodium falciparumlife cycle. PLoS Pathog. 2010, 6: e1000737-10.1371/journal.ppat.1000737.
Mok S, Imwong M, Mackinnon M, Sim J, Ramadoss R, Yi P, Mayxay M, Chotivanich K, Liong K-Y, Russell B, Socheat D, Newton P, Day N, White N, Preiser P, Nosten F, Dondorp A, Bozdech Z: Artemisinin resistance inPlasmodium falciparumis associated with an altered temporal pattern of transcription. BMC Genomics. 2011, 12: 391-10.1186/1471-2164-12-391.
van Brummlen T, Becker J, Mancama D, Hoppe H: Investigating transcriptional regulation ofPlasmodium falciparumupon drug perturbation. Malar J. 2010, 9: O32-10.1186/1475-2875-9-S2-O32.
Mok BW, Ribacke U, Winter G, Yip BH, Tan CS, Fernandez V, Chen Q, Nilsson P, Wahlgren M: Comparative transcriptomal analysis of isogenicPlasmodium falciparumclones of distinct antigenic and adhesive phenotypes. Mol Biochem Parasitol. 2007, 151: 184-192. 10.1016/j.molbiopara.2006.11.006.
Silvestrini F, Bozdech Z, Lanfrancotti A, Di Giulio E, Bultrini E, Picci L, Derisi JL, Pizzi E, Alano P: Genome-wide identification of genes upregulated at the onset of gametocytogenesis inPlasmodium falciparum. Mol Biochem Parasitol. 2005, 143: 100-110. 10.1016/j.molbiopara.2005.04.015.
Yuda M, Iwanaga S, Shigenobu S, Mair GR, Janse CJ, Waters AP, Kato T, Kaneko I: Identification of a transcription factor in the mosquito-invasive stage of malaria parasites. Mol Microbiol. 2009, 71: 1402-1414. 10.1111/j.1365-2958.2009.06609.x.
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 thePlasmodium falciparumlife cycle. Genome Res. 2004, 14: 2308-2318. 10.1101/gr.2523904.
Gissot M, Refour P, Briquet S, Boschet C, Coupe S, Mazier D, Vaquero C: Transcriptome of 3D7 and its gametocyte-less derivative F12Plasmodium falciparumclones during erythrocytic development using a gene-specific microarray assigned to gene regulation, cell cycle and transcription factors. Gene. 2004, 341: 267-277.
Young JA, Fivelman QL, Blair PL, de la Vega P, Le Roch KG, Zhou Y, Carucci DJ, Baker DA, Winzeler EA: ThePlasmodium falciparumsexual development transcriptome: a microarray analysis using ontology-based pattern identification. Mol Biochem Parasitol. 2005, 143: 67-79. 10.1016/j.molbiopara.2005.05.007.
Vontas J, Siden-Kiamos I, Papagiannakis G, Karras M, Waters AP, Louis C: Gene expression inPlasmodium bergheiookinetes and early oocysts in a co-culture system with mosquito cells. Mol Biochem Parasitol. 2005, 139: 1-13. 10.1016/j.molbiopara.2004.03.003.
Yuda M, Iwanaga S, Shigenobu S, Kato T, Kaneko I: Transcription Factor AP2-Sp and its Target Genes in Malarial Sporozoites. Mol Microbiol. 2010, 75: 854-863. 10.1111/j.1365-2958.2009.07005.x.
Mikolajczak SA, Silva-Rivera H, Peng X, Tarun AS, Camargo N, Jacobs-Lorena V, Daly TM, Bergman LW, de la Vega P, Williams J, Aly AS, Kappe SH: Distinct malaria parasite sporozoites reveal transcriptional changes that cause differential tissue infection competence in the mosquito vector and mammalian host. Mol Cell Biol. 2008, 28: 6196-6207. 10.1128/MCB.00553-08.
Tarun AS, Peng X, Dumpit RF, Ogata Y, Silva-Rivera H, Camargo N, Daly TM, Bergman LW, Kappe SH: A combined transcriptome and proteome survey of malaria parasite liver stages. Proc Natl Acad Sci U S A. 2008, 105: 305-310. 10.1073/pnas.0710780104.
Albuquerque S, Carret C, Grosso A, Tarun A, Peng X, Kappe S, Prudencio M, Mota M: Host cell transcriptional profiling during malaria liver stage infection reveals a coordinated and sequential set of biological events. BMC Genomics. 2009, 10: 270-10.1186/1471-2164-10-270.
Aurrecoechea C, Brestelli J, Brunk BP, Dommer J, Fischer S, Gajria B, Gao X, Gingle A, Grant G, Harb OS, Heiges M, Innamorato F, Iodice J, Kissinger JC, Kraemer E, Li W, Miller JA, Nayak V, Pennington C, Pinney DF, Roos DS, Ross C, Stoeckert CJ, Treatman C: Wang H: PlasmoDB: a functional genomic database for malaria parasites. Nucleic Acids Res. 2009, 37: D539-543. 10.1093/nar/gkn814.
Li F, Sonbuchner L, Kyes SA, Epp C, Deitsch KW: Nuclear non-coding RNAs are transcribed from the centromeres ofPlasmodium falciparumand are associated with centromeric chromatin. J Biol Chem. 2008, 283: 5692-5698.
Otto TD, Wilinski D, Assefa S, Keane TM, Sarry LR, Böhme U, Lemieux J, Barrell B, Pain A, Berriman M, Newbold C, Llinás M: New insights into the blood-stage transcriptome ofPlasmodium falciparumusing RNA-Seq. Mol Microbiol. 2010, 76: 12-24. 10.1111/j.1365-2958.2009.07026.x.
Upadhyay R, Bawankar P, Malhotra D, Patankar S: A screen for conserved sequences with biased base composition identifies noncoding RNAs in the A–T rich genome ofPlasmodium falciparum. Mol Biochem Parasitol. 2005, 144: 149-158. 10.1016/j.molbiopara.2005.08.012.
Mourier T, Carret C, Kyes S, Christodoulou Z, Gardner PP, Jeffares DC, Pinches R, Barrell B, Berriman M, Griffiths-Jones S, Ivens A, Newbold C, Pain A: Genome-wide discovery and verification of novel structured RNAs inPlasmodium falciparum. Genome Res. 2008, 18: 281-292. 10.1101/gr.6836108.
Raabe CA, Sanchez CP, Randau G, Robeck T, Skryabin BV, Chinni SV, Kube M, Reinhardt R, Ng GH, Manickam R, Kuryshev VY, Lanzer M, Brosius J, Tang TH, Rozhdestvensky TS: A global view of the nonprotein-coding transcriptome inPlasmodium falciparum. Nucleic Acids Res. 2010, 38: 608-617. 10.1093/nar/gkp895.
RepeatMasker. , ,http://www.repeatmasker.org,
Hu G, Llinas M, Li J, Preiser P, Bozdech Z: Selection of long oligonucleotides for gene expression microarrays using weighted rank-sum strategy. BMC Bioinforma. 2007, 8: 350-10.1186/1471-2105-8-350.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
eArray. , ,http://earray.chem.agilent.com,
Trager W, Jensen JB: Human malaria parasites in continuous culture. Science. 1976, 193: 673-675. 10.1126/science.781840.
Haynes JD, Moch JK: Automated synchronization ofPlasmodium falciparumparasites by culture in a temperature-cycling incubator. Methods Mol Med. 2002, 72: 489-497.
Llinás Manuel: Laboratory Protocols. http://www.molbio1.princeton.edu/labs/llinas/protocols.html
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.
Branham W, Melvin C, Han T, Desai V, Moland C, Scully A, Fuscoe J: Elimination of laboratory ozone leads to a dramatic improvement in the reproducibility of microarray gene expression measurements. BMC Biotechnol. 2007, 7: 8-10.1186/1472-6750-7-8.
Claessens A, Ghumra A, Gupta AP, Mok S, Bozdech Z, Rowe JA: Design of a variant surface antigen-supplemented microarray chip for whole transcriptome analysis of multiplePlasmodium falciparumcytoadherent strains, and identification of strain-transcendent rif and stevor genes. Malar J. 2011, 10: 180-10.1186/1475-2875-10-180.
Bozdech Z, Llinas M, Pulliam BL, Wong ED, Zhu J, DeRisi JL: The transcriptome of the intraerythrocytic developmental cycle ofPlasmodium falciparum. PLoS Biol. 2003, 1: E5-
We thank Lindsey Altenhofen, Donna Storton, Rohan Bansal and Laurie Kramer for technical support, Dr Ulrike Boehme for the March 2011 Plasmodium berghei transcript annotations, and Drs Andy Waters, Chris Janse, and Oliver Billker for input on designing the Plasmodium berghei array. We apologize to the numerous authors whose work we were unable to cite due to space limitations. BK is supported by an HHMI Fellowship of the Damon Runyon Cancer Research Foundation. ML is funded by the NIH R01 AI076276, NIH P50 GM0715078, and NIH Director’s New Innovators Award 1 DP2 OD 001315.
The authors declare they have no conflicts of interest.
BK designed the arrays, analysed the data, and contributed equally to the generation of figures and tables. HJP developed the experimental protocols, performed the experiments and contributed equally to the generation of figures and tables. ML contributed to the design of the microarrays and analysis. All authors contributed to experimental design and to writing the manuscript. All authors read and approved the final manuscript.
Björn F C Kafsack, Heather J Painter contributed equally to this work.