Defining species specific genome differences in malaria parasites
© Liew et al; licensee BioMed Central Ltd. 2010
Received: 10 July 2009
Accepted: 23 February 2010
Published: 23 February 2010
In recent years a number of genome sequences for different plasmodium species have become available. This has allowed the identification of numerous conserved genes across the different species and has significantly enhanced our understanding of parasite biology. In contrast little is known about species specific differences between the different genomes partly due to the lower sequence coverage and therefore relatively poor annotation of some of the draft genomes particularly the rodent malarias parasite species.
To improve the current annotation and gene identification status of the draft genomes of P. berghei, P. chabaudi and P. yoelii, we performed genome-wide comparisons between these three species. Through analyses via comparative genome hybridizations using a newly designed pan-rodent array as well as in depth bioinformatics analysis, we were able to improve on the coverage of the draft rodent parasite genomes by detecting orthologous genes between these related rodent parasite species. More than 1,000 orthologs for P. yoelii were now newly associated with a P. falciparum gene. In addition to extending the current core gene set for all plasmodium species this analysis also for the first time identifies a relatively small number of genes that are unique to the primate malaria parasites while a larger gene set is uniquely conserved amongst the rodent malaria parasites.
These findings allow a more thorough investigation of the genes that are important for host specificity in malaria.
Malaria is a disease caused by the parasitic protozoa from the genus Plasmodium. While the disease is restricted to the tropical and sub-tropical regions of the world due mainly to the natural habitat of the mosquito vector, these regions are densely populated with almost 2.2 billion people living in endemic areas and 515 million cases were expected per annum . Although well-established culture and molecular techniques have been established for Plasmodium falciparum, the use of rodent malaria parasites as in vivo models in the study of the host-parasite interactions is still as relevant today because the rodent parasites are very similar to the human and primate parasites in terms of life cycle, physiology and structure .
Since the release of the genome sequences of P. falciparum and the rodent malaria species, common features of these haploid genomes include a genome size of 22-26 Mb that are arranged in 14 chromosomes ranging from 0.5-3.0 Mb. In addition the current genomic data show a high degree of conservation between different Plasmodium species with the exception of genes located in the telomeric and subtelomeric regions that are extremely variable due to their role in antigenic variation and immune evasion [3, 4]. This suggests that genes located close to the centromeres of the chromosomes would be highly conserved. As a proof of concept, comparisons of these centrally located genes between P. falciparum and the most completely annotated rodent parasite species P. yoelii was shown to have a high degree of synteny . Thus, conservation of these common 'core' genes even amongst divergent species was demonstrated and differences were mainly due to chromosomal re-arrangements [3, 4].
In contrast, the chromosomal regions responsible for antigenic variation and host immune evasion show the most divergence. The genome of P. falciparum contains species-specific subtelomeric genes involved in host cell invasion, adhesion and antigenic variation that are not found in the P. yoelii genome. For example, the P. falciparum genes that are located in the sub-telomeric regions include the var, stevors and rifins that are responsible for antigenic variation and hence immune evasion [6, 7]. In contrast, the P. yoelii yir gene family seems to be the largest multigene family found in the sub-telomeric regions and is absent in P. falciparum. Interestingly, this multigene family is common to all the rodent species and to P. vivax . A more recent study employing probabilistic modelling in conjunction with genomic organization and protein structure analysis bridges this gap and places the P. falciparum rjf/stevor multigene family together into a conserved multigene superfamily of malaria parasites known as the Plasmodium interspersed repeats (PIRs) . This discovery further suggests that although genes are more conserved than previously thought; there is evidence of species-specific divergence that is dependent on each species' interaction with the host or host immune system. A recent phylogenetic survey of rodent malaria parasites based on DNA sequences from multiple loci in the nuclear, mitochondrial and plastid genomes placed P. berghei and P. yoelii as sister species forming a distinct clade while P. chabaudi and another rodent parasite P. vinckei forming another group  thus suggesting that P. berghei and P. yoelii seemed to be more evolutionarily related to each other while P. chabaudi appeared to be a more distinct line. While these studies address gene conservation among different species as well as confirm species-specific sequence polymorphisms a global quantitative comparisons with regards to gain or loss of gene function have not been attempted due to the incomplete genome sequences of the malaria parasite species, barring the high resolution map of P. falciparum . As a consequence, differences between the rodent parasites especially at the genomic level have not been fully elucidated as yet and thus linking genomic differences to phenotypic traits of these parasites have been difficult. In order to address this issue, genome-wide comparisons including comparative genomic hybridization (CGH) and detailed bioinformatics analysis can be employed as a valuable tool to identify similarities and differences between related species. This technique has been utilized to compare genomic differences between similar species in both prokaryotes and eukaryotes. For example, CGH has been used to genotype related yeast species to determine the presence and absence/polymorphic sequences . The usefulness of CGH in gene discovery was shown with the discovery of three thousand novel genes in Klebsiella pneuomoniae 342 using microarray hybridization with Escherichia coli K-12 open reading frames with confirmation of the presence of genes coding for conserved metabolic functions and demonstrating specificity where genes obtained via lateral transfer in K-12 were absent in 342 . In addition, CGH is also able to differentiate genome-wide variation in various virulent and avirulent Burkholderia species . Besides analyzing sequence variation, the gain/loss of DNA can also be used as a tool for elucidating evolutionary divergence between species in their natural environment.
The means for such analysis lies in the use of DNA microarray technology based on the 'spotting' of long oligonucleotides onto glass slides [15–17]. Although the techniques employed for performing DNA microarray experiments have been well established, improvements in the algorithms used for designing the array probes have been critical for advancing the reproducibility and accuracy of detecting the respective gene targets. Recently, a novel robust program called 'OligoRankPick' was able to optimize oligonucleotide design by optimizing target specificity and GC% variation. This algorithm is based on a weighted rank-sum strategy to optimize oligonucleotide selection even along genomes of many diverse organisms . The resultant P. falciparum dataset generated from this strategy was shown to be highly reproducible and also led to increased coverage of the P. falciparum genome as compared to earlier P. falciparum microarray designs. Although the incompleteness of these rodent parasite genomes could present difficulties in their comparative analysis, there is some evidence that they are highly conserved . We have thus emulated the 'OligoRankPick' strategy to design a pan-rodent cross-genome oligonucleotide microarray for the rodent malaria species P. berghei, P. chabaudi and P. yoelii. Due to the incompleteness of the genome sequence available for all three rodent genomes the microarray design approach resulted in oligonucleotide sequences predicted to be complementary to all three rodent species, two rodent malaria species or to be unique to a single species. Comparative genomic hybridization using genomic DNA from each of the three species not only validated the oligonucleotides but at the same time provided new information that allowed the closure of sequence gaps in the genomes of one or more of the rodent species based on complementary hybridization data, thereby significantly improving our overall understanding of gene content for each species. While the CGH approach proofed exceptionally powerful it was still likely that due to sequence polymorphisms found in the different species at the oligonucleotide probe region genes that are actually present would be missed using this strategy alone. For this reason the current available genome sequences were reanalyzed using bioinformatics tools to detect any additional genes missed by the microarray approach alone. Finally, a random selection of genes predicted from the microarray and bioinformatics data to now be present in all rodent genome species was validated using PCR and sequencing.
Overall this multipronged approach allowed us to significantly improve the gene predictions for all three rodent malaria species. In addition the new information obtained in this study allows the assembly of an improved core set of genes present in all plasmodium species while at the same time also identifying a rodent malaria specific gene set.
Key design features of the 'Pan-rodent' array chip
Comparative Genomic Hybridization using the Pan-rodent microarray
Since some of the oligonucleotides are designed to be capable of cross hybridizing with more than one species, we have re-sorted all the genes from the three species and matched them with their respective oligonucleotides so that an accurate normalized value can be obtained for analysis. Using this approach, we performed comparative genomic hybridization experiments and first looked at the performance of the species-specific oligonucleotides that were designed to hybridize with their target genome. After normalizing and filtering the data the percentage of successful oligonucleotides for P. yoelii- specific, P. berghei-specific and P. chabaudi-specific oligonucleotides were 91% (7,347), 90% (6,314) and 84% (6,356) respectively. It is not surprising that a low proportion of oligonucleotides were unable to hybridize with their intended targets especially if they were designed in regions with low confidence in sequence quality. A likely reason is the relatively poor genome coverage of the rodent parasite genomes as compared to the more complete P. falciparum genome. Also, there are errors in the current sequence drafts of the rodent malaria parasites due to the propensity of random sequence rearrangements of an AT-rich genome in a sequencing vector.
Using the validated set of oligonucleotides, we identified all the oligonucleotides that hybridized to the DNA of parasites they had not originally been designed against (Figure 2). A positive hybridization signal provides strong evidence that the sequences and therefore the genes represented by these oligonucleotides are also present in the genomes of the other rodent parasite species. This strongly implies that the respective genes that are represented by these oligonucleotides are actually present but are not found in the current database.
Filling the gaps in the genomes of the rodent malaria parasite species
Since it has been established that the malaria parasite genomes are well conserved , it is conceivable that genes that are missing in the current draft of either of the RMP genome have a well conserved ortholog in one or two of the other species. Using the pan-rodent microarray we wish to investigate this possibility by inspecting comparative genome hybridization (CGH) signals on "cross-species" oligonucleotide microarray elements. For example high signal from P. yoelii gDNA for an oligonucleotide that is not designed to hybridize to P. yoelii implies that this gene (sequence) is present in P. yoelii. In summary, a species where the draft genome currently does not possess a particular gene but now gives a signal on the array suggests the presence of this orthologous gene. Hence, this gene can thus be detected based on homology.
Missing genes in the draft genome could arise due to two scenarios: either the sequence information is missing, or that the sequence is present in the genome but missed by current gene prediction algorithms. Since the oligonucleotides were designed based on predicted open reading frames, a CGH signal constitutes direct experimental evidence for the presence of an orthologous genomic sequence and thus potentially the gene in the RMP genome in which this gene is missing. Based on this approach, we detect 179, 306 and 215 genes missing in the current draft sequences of the P. yoelii, the P. berghei and the P. chabaudi genome, respectively. The majority of these genes code for hypothetical proteins that lack any functional assignment based on their amino acid sequence. For those whose function has been described, genes involved in biosynthesis, protein modifications, kinases and also invasion-related proteins in the case of P. chabaudi and P. berghei have been discovered (Additional file 1: Supplemental Table S1).
Cross-species gene identification using bioinformatics
Defining the core Plasmodium genome
With the more complete coverage of the rodent malaria parasite genomes obtained via the analysis here, we now thought to define a common 'core' Plasmodium genome where orthologous genes present in all sequenced malaria parasites was defined. Plasmodium falciparum, being the most completely annotated and researched malaria parasite genome, is used as an index reference species in an attempt to consolidate genes common to as many species as possible. For each P. falciparum gene the orthologs present in the rodent parasite species were matched, and the resulting dataset was further appended with the respective orthologs from another human malaria parasite P. vivax and the simian parasite P. knowlesi as these genomes had been recently sequenced to 10× and 8× [21, 22] coverage respectively. This approach segregated P. falciparum genes into two groups: firstly, the 'core' genes containing orthologs to the rodent parasite as well as the P. vivax and P. knowlesi sequences (Additional file 2: Supplemental Table S2) and secondly, genes that have no significant alignment with rodent sequences (Additional file 3: Supplemental Table S3). While orthologous genes (with respect to P. falciparum) are highly conserved in regions where housekeeping genes predominate, the telomeric and sub-telomeric regions (i.e. non-syntenic) contain genes involved in antigenic variation that are more divergent  and this is consistent with the observation that the group representing the non-orthologous genes is dominated by var, stevor and rifin gene families that are involved in antigenic variation. What is immediately apparent is the high similarity of all six parasite species in the regions that contain housekeeping genes. In total, 4,188 genes were common between P. falciparum and the rodent parasites. Of these, 73 genes, mainly hypothetical genes, were absent in both P. vivax and P. knowlesi. In addition, 50 genes were present in P. vivax but not in P. knowlesi and 79 genes were present in P. knowlesi but not in P. vivax. Without these species-specific genes, we found 3,986 genes that are common to all six species.
Syntenic Pf genes with corresponding RMP chromosome.
No. of genes
Corresponding Pv genes
Corresponding Pk genes
Non-syntenic Pf genes without corresponding RMP chromosome.
No. of genes
Corresponding Pv genes
Corresponding Pk genes
Defining the core rodent Plasmodium genome
Overall this approach now identified 1,238 genes that are common to all rodent malaria species most of which (1,013) represent hypothetical genes. Importantly only about 8% of these predicted genes have ORF that are shorter than 100 nucleotides making it less likely that these are annotation or prediction errors but indeed represent functional genes. Taking together these gene numbers along with the 4,188 genes found in the 'core' set, the total number of conserved rodent malaria genes is 5,426. In addition it appears that there are 210 genes specific to P. yoelii and P. chabaudi (Additional file 5: Supplemental Table S5), 247 genes specific to P. yoelii and P. berghei (Additional file 6: Supplemental Table S6) and 453 unique to P. yoelii (Additional file 7: Supplemental Table S7) with most of them (> 90%) again representing hypothetical genes. Of the 453 unique P. yoelii genes 45% are less then 100 nucleotide in length, compared to 24% in the P. yoelii and P. berghei group and 17% in the P. yoelii and P. chabaudi group indicating that a number of these genes represent false gene models.
Our microarray specific to the rodent malaria species P. berghei, P. chabaudi and P. yoelii has provided a golden opportunity to verify the genomes of these rodent malaria parasite (RMP) species. The RMP draft genomes were annotated similarly where GlimmerMExon was first trained on P. falciparum data and automated assignments were based on hidden Markov model association [3, 25]. Since these genomes were annotated in a similar manner, we do not foresee any annotation bias in the three RMP draft genomes. Our oligonucleotide selection process aims to include every predicted gene model and the robustness of the design algorithm has been shown in the ability to differentiate between members of large multigene families . Evidence of significant homology via hybridization and bioinformatics show that a high proportion of these gene models are conserved amongst the three RMP species, therefore suggesting that a high proportion represent true gene models that are refractory to random mutations as compared to intronic sequences.
The definition of the core gene set of Plasmodia genes using the most completely annotated genome of P. falciparum as a reference index has enabled the use of direct hybridization and bioinformatics tools to expand the repertoire of common genes. The validity of utilizing CGH and bioinformatics to improve genome annotation is well established [12–14, 26]. Novel gene models discovered via hybridization have been validated with PCR (Figure 3) while genes discovered via both hybridization and bioinformatics were also validated (Figures 4 and 5). In fact, PCR screens of a subset of genes without at least one RMP ortholog suggest that there is high confidence that any P. falciparum gene that contains at least an ortholog from one of the RMP species can potentially contain orthologs from all three RMP species. This suggests that the high stringency thresholds in both hybridization and bioinformatics are in fact conservative and there are potentially more genes that have yet to be discovered. Hence, the results indicate that the genetic repertoires of the RMP species are indeed more common and that these genes should be present in the poorly assembled and annotated rodent parasite genomes.
The construction of a rodent specific orthology map using P. yoelii as the reference index was undertaken to study and survey rodent parasite specific genes that are distinct from the common 'core' set. Although this process would filter off genes specific to the other two RMP species, the P. yoelii genome was chosen as it is the most completely annotated rodent parasite genome while the gene set in P. chabaudi is over-predicted thus making comparisons across species difficult. This associative table clearly demonstrates certain species-specific metabolic differences and reveals gene duplications and expansions from both the core set of Plasmodia genes and the rodent parasite specific genes. Similarly, hybridization data reveals more orthologous genes than the current known set.
Due to the high AT-bias in the genome of the Plasmodia species, sequencing the entire genome, joining overlapping sequences to form contigs and gene predictions have been difficult. While the majority of the rodent malaria parasite genes have been sequenced and annotated, many gaps still remain and the data from the pan-rodent cross-genome oligonucleotide microarray provides direct experimental evidence for this case. While the genomes from P. berghei and P. chabaudi continues to be refined at the Sanger Institute, improvements on the sequence coverage of the P. yoelii genome has not been taken up by any group. Therefore, any improvements in genome coverage for P. yoelii would be extremely valuable.
Orthologs between P. falciparum and three rodent malaria parasite species.
Pb vs. Py
Pb vs. Pc
Pc vs. Py
Py vs. Pf
Pb vs. Pf
Pc vs. Pf
<3 RMP orthologs
Potential 'Core' set
Our approach here provides a significant improvement of gene coverage using existing information and well-established experimental and analysis techniques. The additional information provided here will be of particular use for future efforts using next-generation DNA sequencing technologies.
Unlike most previous studies that have focused on finding similarities between the genomes of different plasmodium species [3, 5, 21, 22, 25] the work here for the first time attempts to define the difference in gene content between the different parasites. Currently, our understanding on the genetic basis for species specific restriction of malaria parasites is limited. This is mainly due to the incomplete status of the RMP genomes that do not allow reasonable assessment of their nuclear encoded proteomes. Conducting the cross species CGH and the reciprocal homology searches allowed us to generate substantially improved dataset that enables us for the first time to investigate the subsets of genes that may be involved in Plasmodium speciation. First, we identify 117 genes that are found in only the human and primate malaria species should be considered key candidates for separating rodent from primate malarias. Similarly the 1,238 RMP specific genes are likely to provide the biological framework that separates the rodent from primate malaria. How these mainly hypothetical genes function in the parasite biology now needs to be further investigated.
In contrast to the relatively large numbers of genes that are unique to individual species (e.g. 927 gene unique to P. falciparum) or species subgoups (e.g. 1,238 specific to RMP), there is only a small number of genes that are lost from the generic core gene set 27 to 96 genes lost in one or two Plasmodium sp.). With the rodent parasites forming a distinct phylogenetic clade and genetic evidence that P. vivax and P. knowlesi are more similar to each other [27, 28] this suggests that species specific gene expansion and diversification is the main driving force of speciation amongst Plasmodium species. This is likely to apply to genes responsible for invasion and antigenic variation, but also genes involved in other species specific phenotypes such as the asynchronous character of the intra-erythrocytic development of P. yoelii, the formation of hypnozoites in P. vivax and the propensity of sequestration.
A modification of the 'OligoRankPick' program  was employed to design 60-mer probes to the three rodent parasite species P. berghei, P. chabaudi and P. yoelii (Additional file 8: Supplemental Table S8; Additional file 9 - Supplementary Table S9) using sequences deposited in PlasmoDB http://www.plasmodb.org. Cross-species oligonucleotides were selected based on the criteria as described in the results section, i.e. at least 90% homology to the target sequence, less that 37.5% to non-target sequences and capped with a 5% difference tolerance in GC content. All possible oligonucleotides for a particular gene was then ranked and sorted based on target (>90%) and non-target hit (<37.5%), GC content (±5%) and the best Smith-Waterman score (Figure 1). Hence, all possible oligonucleotides were evaluated based on every available gene model and the best sequence was then selected. Oligonucleotides were spotted onto poly-L-lysine-coated microscopic glass slides .
Rodent parasite DNA preparation for genomic DNA microarray
Balb/c mice were used as vertebrate hosts for the propagation of P. berghei ANKA, P. chabaudi AS and P. yoelii 17 × 1.1. All procedures were in accordance with the guidelines for the use of experimental animals established by the Institutional Animal Care and Use Committee (IACUC) of the Nanyang Technological University of Singapore. Leukocytes were filtered from whole blood  and gDNA extracted using the Easy-DNA™ kit (Invitrogen) according to the manufacturer's protocol. For each parasite line, 3 μg of DNA was mixed with 2.5 μg of random nanomer primer to a volume of 15.25 μl. DNA was denatured at 100°C for 5 min and snap cooled on ice for 5 min. Labeled DNA was generated by adding dNTPs to a final concentration of 1 mM dATP and 500 fM each: dCTP, dGTP, dTTP and 5-(3-aminoallyl)¬2'-deoxyuridine-5'-triophosphate, (aa-dUTP) (Biotium), with 20 units of exo-Klenow Fragment (New England Biolabs) in a total volume of 50 μl and incubated at room temperature for 10 min and then overnight at 37°C. Cy-dye coupling, cleanup, hybridization, washing and slide scanning were performed as described by Bozdech and co-workers . Experiments were performed with at least duplicates for each species Subsequently, the data were normalized using the NOMAD microarray database http://ucsf-nomad.sourceforge.net/. For the comparative genome analysis, low quality features and features with a signal level less than two-fold of the background plus two-fold of the background standard deviation was filtered from the initial raw data set.
Polymerase Chain Reaction (PCR)
Verification of genes was performed via PCR thermocycling (Eppendorf) using the following conditions: 1 cycle of 95°C for 2 min; followed by 35 cycles of 95°C for 1 min, 50°C for 1 min and 68°C for 30 s. A final extension was set for 10 min at 68°C and then the tubes were kept at 4°C. For each 20 μl reaction, 2 ng of DNA was used as a template with 10 pmol of the forward and reverse primers (Additional file 10: Supplemental Table S10), 200 μM of dNTPs with 1.5 units of Taq DNA polymerase (Kapa Biosystems).
The sequences of P. falciparum, P. knowlesi and P. vivax were obtained from PlasmoDB http://www.plasmodb.org. For any gene missing an orthologous partner in another species even after complementation with the array data, known orthologs were appended using the resource available at PlasmoDB (version 5.5). Unannotated genes were identified using the tBLASTn search where amino acid sequences from a well-annotated species (PlasmoDB version 5.5) are used to query the genome of another species. A threshold of 10-15 was used to ensure the stringency of locating matching genes.
Availability of microarray data
Microarray data are publically available at the Centre for Information Biology Gene Expression Database (CIBEX; http://cibex.nig.ac.jp). The Accession number for these data is CBX114.
This work was funded by a grant of the National Medical Research Council of Singapore.
- Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI: The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005, 434 (7030): 214-217. 10.1038/nature03342.PubMed CentralPubMedView ArticleGoogle Scholar
- Carter R, Diggs C: Plasmodia of rodents. Parasitic Protozoa. Edited by: Kreier J. 1977, New York: Academic Press, 3: 359-465.Google Scholar
- Carlton JM, Angiuoli SV, Suh BB, Kooij TW, Pertea M, Silva JC, Ermolaeva MD, Allen JE, Selengut JD, Koo HL: Genome sequence and comparative analysis of the model rodent malaria parasite Plasmodium yoelii yoelii. Nature. 2002, 419 (6906): 512-519. 10.1038/nature01099.PubMedView ArticleGoogle Scholar
- Waters AP: Orthology between the genomes of Plasmodium falciparum and rodent malaria parasites: possible practical applications. Philos Trans R Soc Lond B Biol Sci. 2002, 357 (1417): 55-63. 10.1098/rstb.2001.1011.PubMed CentralPubMedView ArticleGoogle Scholar
- Kooij TW, Carlton JM, Bidwell SL, Hall N, Ramesar J, Janse CJ, Waters AP: A Plasmodium whole-genome synteny map: indels and synteny breakpoints as foci for species-specific genes. PLoS Pathog. 2005, 1 (4): e44-10.1371/journal.ppat.0010044.PubMed CentralPubMedView ArticleGoogle Scholar
- Blythe JE, Surentheran T, Preiser PR: STEVOR--a multifunctional protein?. Molecular and biochemical parasitology. 2004, 134 (1): 11-15. 10.1016/j.molbiopara.2003.09.011.PubMedView ArticleGoogle Scholar
- Kyes S, Horrocks P, Newbold C: Antigenic variation at the infected red cell surface in malaria. Annu Rev Microbiol. 2001, 55: 673-707. 10.1146/annurev.micro.55.1.673.PubMedView ArticleGoogle Scholar
- Janssen CS, Barrett MP, Turner CM, Phillips RS: A large gene family for putative variant antigens shared by human and rodent malaria parasites. Proc R Soc Lond B Biol Sci. 2002, 269 (1489): 431-436. 10.1098/rspb.2001.1903.View ArticleGoogle Scholar
- Janssen CS, Phillips RS, Turner CM, Barrett MP: Plasmodium interspersed repeats: the major multigene superfamily of malaria parasites. Nucleic Acids Res. 2004, 32 (19): 5712-5720. 10.1093/nar/gkh907.PubMed CentralPubMedView ArticleGoogle Scholar
- Perkins SL, Sarkar IN, Carter R: The phylogeny of rodent malaria parasites: simultaneous analysis across three genomes. Infect Genet Evol. 2007, 7 (1): 74-83. 10.1016/j.meegid.2006.04.005.PubMedView ArticleGoogle Scholar
- Gardner MJ, Hall N, Fung E, White O, Berriman M, Hyman RW, Carlton JM, Pain A, Nelson KE, Bowman S: Genome sequence of the human malaria parasite Plasmodium falciparum. Nature. 2002, 419 (6906): 498-511. 10.1038/nature01097.PubMedView ArticleGoogle Scholar
- Daran-Lapujade P, Daran JM, Kotter P, Petit T, Piper MD, Pronk JT: Comparative genotyping of the Saccharomyces cerevisiae laboratory strains S288C and CEN.PK113-7D using oligonucleotide microarrays. FEMS Yeast Res. 2003, 4 (3): 259-269. 10.1016/S1567-1356(03)00156-9.PubMedView ArticleGoogle Scholar
- Dong Y, Glasner JD, Blattner FR, Triplett EW: Genomic interspecies microarray hybridization: rapid discovery of three thousand genes in the maize endophyte, Klebsiella pneumoniae 342, by microarray hybridization with Escherichia coli K-12 open reading frames. Appl Environ Microbiol. 2001, 67 (4): 1911-1921. 10.1128/AEM.67.4.1911-1921.2001.PubMed CentralPubMedView ArticleGoogle Scholar
- Ong C, Ooi CH, Wang D, Chong H, Ng KC, Rodrigues F, Lee MA, Tan P: Patterns of large-scale genomic variation in virulent and avirulent Burkholderia species. Genome Res. 2004, 14 (11): 2295-2307. 10.1101/gr.1608904.PubMed CentralPubMedView ArticleGoogle 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 (2): R9-10.1186/gb-2003-4-2-r9.PubMed CentralPubMedView ArticleGoogle Scholar
- Hughes TR, Mao M, Jones AR, Burchard J, Marton MJ, Shannon KW, Lefkowitz SM, Ziman M, Schelter JM, Meyer MR: Expression profiling using microarrays fabricated by an ink-jet oligonucleotide synthesizer. Nat Biotechnol. 2001, 19 (4): 342-347. 10.1038/86730.PubMedView ArticleGoogle Scholar
- Kane MD, Jatkoe TA, Stumpf CR, Lu J, Thomas JD, Madore SJ: Assessment of the sensitivity and specificity of oligonucleotide (50 mer) microarrays. Nucleic Acids Res. 2000, 28 (22): 4552-4557. 10.1093/nar/28.22.4552.PubMed CentralPubMedView ArticleGoogle Scholar
- Hu G, Llinas M, Li J, Preiser PR, Bozdech Z: Selection of long oligonucleotides for gene expression microarrays using weighted rank-sum strategy. BMC bioinformatics. 2007, 8: 350-10.1186/1471-2105-8-350.PubMed CentralPubMedView ArticleGoogle Scholar
- Bozdech Z, Llinas M, Pulliam BL, Wong ED, Zhu J, DeRisi JL: The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparum. PLoS Biol. 2003, 1 (1): E5-10.1371/journal.pbio.0000005.PubMed CentralPubMedView ArticleGoogle Scholar
- Bozdech Z, Mok S, Hu G, Imwong M, Jaidee A, Russell B, Ginsburg H, Nosten F, Day NP, White NJ: The transcriptome of Plasmodium vivax reveals divergence and diversity of transcriptional regulation in malaria parasites. Proceedings of the National Academy of Sciences of the United States of America. 2008, 105 (42): 16290-16295. 10.1073/pnas.0807404105.PubMed CentralPubMedView ArticleGoogle Scholar
- Carlton JM, Adams JH, Silva JC, Bidwell SL, Lorenzi H, Caler E, Crabtree J, Angiuoli SV, Merino EF, Amedeo P: Comparative genomics of the neglected human malaria parasite Plasmodium vivax. Nature. 2008, 455 (7214): 757-763. 10.1038/nature07327.PubMed CentralPubMedView ArticleGoogle Scholar
- Pain A, Bohme U, Berry AE, Mungall K, Finn RD, Jackson AP, Mourier T, Mistry J, Pasini EM, Aslett MA: The genome of the simian and human malaria parasite Plasmodium knowlesi. Nature. 2008, 455 (7214): 799-803. 10.1038/nature07306.PubMed CentralPubMedView ArticleGoogle Scholar
- Cunningham D, Lawton J, Jarra W, Preiser P, Langhorne J: The pir multigene family of Plasmodium: antigenic variation and beyond. Molecular and biochemical parasitology. 2010, 170 (2): 65-73. 10.1016/j.molbiopara.2009.12.010.PubMedView ArticleGoogle Scholar
- Cunningham DA, Jarra W, Koernig S, Fonager J, Fernandez-Reyes D, Blythe JE, Waller C, Preiser PR, Langhorne J: Host immunity modulates transcriptional changes in a multigene family (yir) of rodent malaria. Mol Microbiol. 2005, 58 (3): 636-647. 10.1111/j.1365-2958.2005.04840.x.PubMedView ArticleGoogle Scholar
- Hall N, Karras M, Raine JD, Carlton JM, Kooij TW, Berriman M, Florens L, Janssen CS, Pain A, Christophides GK: A comprehensive survey of the Plasmodium life cycle by genomic, transcriptomic, and proteomic analyses. Science. 2005, 307 (5706): 82-86. 10.1126/science.1103717.PubMedView ArticleGoogle Scholar
- Vaughan A, Chiu SY, Ramasamy G, Li L, Gardner MJ, Tarun AS, Kappe SH, Peng X: Assessment and improvement of the Plasmodium yoelii yoelii genome annotation through comparative analysis. Bioinformatics. 2008, 24 (13): i383-389. 10.1093/bioinformatics/btn140.PubMed CentralPubMedView ArticleGoogle Scholar
- Escalante AA, Ayala FJ: Phylogeny of the malarial genus Plasmodium, derived from rRNA gene sequences. Proceedings of the National Academy of Sciences of the United States of America. 1994, 91 (24): 11373-11377. 10.1073/pnas.91.24.11373.PubMed CentralPubMedView ArticleGoogle Scholar
- Roy SW, Irimia M: Origins of human malaria: rare genomic changes and full mitochondrial genomes confirm the relationship of Plasmodium falciparum to other mammalian parasites but complicate the origins of Plasmodium vivax. Mol Biol Evol. 2008, 25 (6): 1192-1198. 10.1093/molbev/msn069.PubMed CentralPubMedView ArticleGoogle Scholar
- Andrews L, Andersen RF, Webster D, Dunachie S, Walther RM, Bejon P, Hunt-Cooke A, Bergson G, Sanderson F, Hill AV: Quantitative real-time polymerase chain reaction for malaria diagnosis and its use in malaria vaccine clinical trials. Am J Trop Med Hyg. 2005, 73 (1): 191-198.PubMedGoogle Scholar