PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of malarjBioMed CentralBiomed Central Web Sitesearchsubmit a manuscriptregisterthis articleMalaria JournalJournal Front Page
 
Malar J. 2012; 11: 187.
Published online Jun 8, 2012. doi:  10.1186/1475-2875-11-187
PMCID: PMC3411454
New Agilent platform DNA microarrays for transcriptome analysis of Plasmodium falciparum and Plasmodium berghei for the malaria research community
Björn F C Kafsack,#1 Heather J Painter,#1 and Manuel Llináscorresponding author1
1Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, 08544, USA
corresponding authorCorresponding author.
#Contributed equally.
Björn F C Kafsack: bkafsack/at/princeton.edu; Heather J Painter: hpainter/at/princeton.edu; Manuel Llinás: manuel/at/genomics.princeton.edu
Received March 8, 2012; Accepted June 8, 2012.
Background
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.
Methods
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.
Results
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.
Conclusions
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 [7] or cDNAs [8] 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 [3], 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[11] 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 [14]. 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 [25] 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[30]. 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 [31] and host transcriptional programs during this developmental stage [32].
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.
Microarray design
Annotated transcripts where downloaded for P. falciparum strain 3D7 (PlasmoDB release 7.1) [33] 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 [39] 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 [40] to design 60mer oligonucleotide probes followed by probe selection using OligoRankPick [40] 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 [40] for details). Probes were then BLASTed against the Crick strand of the target transcripts using NCBI BLAST [41] 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 [42]. 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 [43] 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 [44] 12 hours apart during the preceding cycle and harvested every six hours throughout the 48-hour cycle, starting approximately three hours post-invasion.
RNA extraction
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 [9]. 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 [45]. Single-strand aminoallyl-containing cDNA synthesis and Amersham CyDye-coupling (GE Healthcare, Piscataway, NJ, USA) was carried out as previously described [46]. To eliminate Cy5 degradation by ozone [47], 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.
Array design
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 ([40] and [48]), 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 [33] 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 [39] 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 [40] 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 [40]. 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 [42], 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.
Table 1
Table 1
Overview of array features
Array performance
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.
Figure 1
Figure 1
High overall consistency of log-ratio measurements across a wide range of transcript abundance. Normalized Cy3 (white), Cy5 (black), and log2(Cy5/Cy3) (grey) intensities from triplicate self-hybridization experiments for (A) Six genes with average intensities (more ...)
Figure 2
Figure 2
High consistency of log-ratio measurements across a 20-fold range of hybridized material. Heat map of mean-centered log2(Cy5/Cy3) ratios for six genes representing a 32-fold intensity range (see Figure 1A) across a 20-fold range of hybridized (more ...)
Figure 3
Figure 3
Periodic gene expression. Relative mRNA abundance for 47 periodically expressed reference genes at eight time-points covering the intraerythrocytic development cycle of P. falciparum strain 3D7. Pearson’s coefficient of correlation (r) is indicated (more ...)
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.
Competing interests
The authors declare they have no conflicts of interest.
Authors' contributions
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.
Additional file 1
Excel file of P. falciparum v7.1 8x15K microarray probes.
click here for file(1003K, xlsx)
Additional file 2
Excel file of P. berghei Mar2011 8x15K microarray probes.
Additional file 3
Figure illustrating coverage for three arrays.
Additional file 4
Figures and table describing array-wide technical reproducibility of hybridizations across three replicates.
Additional file 5
Table of titration series of total RNA starting amounts, resultant cDNA generation, dye-coupling, and array hybridization.
Additional file 6
Table of transcripts well-above background by amount hybridized and table of signal intensity correlation across a self-hybridization dilution series.
Additional file 7
List of highly expressed periodic genes used in Figure 3.
Acknowledgements
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.
  • Limenitakis J, Soldati-Favre D. Functional genetics in Apicomplexa: Potentials and limits. FEBS Lett. 2011;585:1579–1588. doi: 10.1016/j.febslet.2011.05.002. [PubMed] [Cross Ref]
  • Carvalho TG, Menard R. Manipulating the Plasmodium genome. Curr Issues Mol Biol. 2005;7:39–55. [PubMed]
  • 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. doi: 10.1038/nature01097. [PubMed] [Cross Ref]
  • 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. doi: 10.1038/nature01099. [PubMed] [Cross Ref]
  • Carlton J. ThePlasmodium vivaxgenome sequencing project. Trends Parasitol. 2003;19:227–231. doi: 10.1016/S1471-4922(03)00066-7. [PubMed] [Cross Ref]
  • 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. doi: 10.1038/nature07306. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1046/j.1365-2958.2000.01730.x. [PubMed] [Cross Ref]
  • 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. doi: 10.1046/j.1365-2958.2001.02222.x. [PubMed] [Cross Ref]
  • 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. [PMC free article] [PubMed]
  • 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. doi: 10.1126/science.1087025. [PubMed] [Cross Ref]
  • 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. doi: 10.1093/nar/gkj517. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1086/522012. [PubMed] [Cross Ref]
  • 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. doi: 10.1038/nature06311. [PubMed] [Cross Ref]
  • 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. doi: 10.1371/journal.ppat.1000644. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1186/gb-2005-6-11-r93. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1186/1471-2164-9-513. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1371/journal.ppat.1000214. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1038/nbt.1597. [PubMed] [Cross Ref]
  • 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. doi: 10.1371/journal.ppat.1000737. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1186/1471-2164-12-391. [PMC free article] [PubMed] [Cross Ref]
  • van Brummlen T, Becker J, Mancama D, Hoppe H. Investigating transcriptional regulation ofPlasmodium falciparumupon drug perturbation. Malar J. 2010;9:O32. doi: 10.1186/1475-2875-9-S2-O32. [Cross Ref]
  • 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. doi: 10.1016/j.molbiopara.2006.11.006. [PubMed] [Cross Ref]
  • 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. doi: 10.1016/j.molbiopara.2005.04.015. [PubMed] [Cross Ref]
  • 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. doi: 10.1111/j.1365-2958.2009.06609.x. [PubMed] [Cross Ref]
  • 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. doi: 10.1101/gr.2523904. [PubMed] [Cross Ref]
  • 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. [PubMed]
  • 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. doi: 10.1016/j.molbiopara.2005.05.007. [PubMed] [Cross Ref]
  • 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. doi: 10.1016/j.molbiopara.2004.03.003. [PubMed] [Cross Ref]
  • 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. doi: 10.1111/j.1365-2958.2009.07005.x. [PubMed] [Cross Ref]
  • 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. doi: 10.1128/MCB.00553-08. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1073/pnas.0710780104. [PubMed] [Cross Ref]
  • 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. doi: 10.1186/1471-2164-10-270. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1093/nar/gkn814. [PMC free article] [PubMed] [Cross Ref]
  • 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. [PubMed]
  • 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. doi: 10.1111/j.1365-2958.2009.07026.x. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1016/j.molbiopara.2005.08.012. [PubMed] [Cross Ref]
  • 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. doi: 10.1101/gr.6836108. [PubMed] [Cross Ref]
  • 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. doi: 10.1093/nar/gkp895. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1186/1471-2105-8-350. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1093/nar/25.17.3389. [PMC free article] [PubMed] [Cross Ref]
  • eArray. , ; http://earray.chem.agilent.com.
  • Trager W, Jensen JB. Human malaria parasites in continuous culture. Science. 1976;193:673–675. doi: 10.1126/science.781840. [PubMed] [Cross Ref]
  • Haynes JD, Moch JK. Automated synchronization ofPlasmodium falciparumparasites by culture in a temperature-cycling incubator. Methods Mol Med. 2002;72:489–497. [PubMed]
  • 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. doi: 10.1186/gb-2003-4-2-r9. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1186/1472-6750-7-8. [PMC free article] [PubMed] [Cross Ref]
  • 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. doi: 10.1186/1475-2875-10-180. [PMC free article] [PubMed] [Cross Ref]
  • 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. [PMC free article] [PubMed]
Articles from Malaria Journal are provided here courtesy of
BioMed Central