Phylogenomic evidence supports past endosymbiosis, intracellular and horizontal gene transfer in Cryptosporidium parvum
© Huang et al.; licensee BioMed Central Ltd. 2004
Received: 19 April 2004
Accepted: 10 September 2004
Published: 19 October 2004
The apicomplexan parasite Cryptosporidium parvum is an emerging pathogen capable of causing illness in humans and other animals and death in immunocompromised individuals. No effective treatment is available and the genome sequence has recently been completed. This parasite differs from other apicomplexans in its lack of a plastid organelle, the apicoplast. Gene transfer, either intracellular from an endosymbiont/donor organelle or horizontal from another organism, can provide evidence of a previous endosymbiotic relationship and/or alter the genetic repertoire of the host organism. Given the importance of gene transfers in eukaryotic evolution and the potential implications for chemotherapy, it is important to identify the complement of transferred genes in Cryptosporidium.
We have identified 31 genes of likely plastid/endosymbiont (n = 7) or prokaryotic (n = 24) origin using a phylogenomic approach. The findings support the hypothesis that Cryptosporidium evolved from a plastid-containing lineage and subsequently lost its apicoplast during evolution. Expression analyses of candidate genes of algal and eubacterial origin show that these genes are expressed and developmentally regulated during the life cycle of C. parvum.
Cryptosporidium is the recipient of a large number of transferred genes, many of which are not shared by other apicomplexan parasites. Genes transferred from distant phylogenetic sources, such as eubacteria, may be potential targets for therapeutic drugs owing to their phylogenetic distance or the lack of homologs in the host. The successful integration and expression of the transferred genes in this genome has changed the genetic and metabolic repertoire of the parasite.
Cryptosporidium is a member of the Apicomplexa, a eukaryotic phylum that includes several important parasitic pathogens such as Plasmodium, Toxoplasma, Eimeria and Theileria. As an emerging pathogen in humans and other animals, Cryptosporidium often causes fever, diarrhea, anorexia and other complications. Although cryptosporidial infection is often self-limiting, it can be persistent and fatal for immunocompromised individuals. So far, no effective treatment is available . Furthermore, because of its resistance to standard chlorine disinfection of water, Cryptosporidium continues to be a security concern as a potential water-borne bioterrorism agent .
Cryptosporidium is phylogenetically quite distant from the hemosporidian and coccidian apicomplexans  and, depending on the molecule and method used, is either basal to all Apicomplexa examined thus far, or is the sister group to the gregarines [4, 5]. It is unusual in several respects, notably for the lack of the apicoplast organelle which is characteristic of all other apicomplexans that have been examined [6, 7]. The apicoplast is a relict plastid hypothesized to have been acquired by an ancient secondary endosymbiosis of a pre-alveolate eukaryotic cell with an algal cell . All that remains of the endosymbiont in Coccidia and Haemosporidia is a plastid organelle surrounded by four membranes . The apicoplast retains its own genome, but this is much reduced (27-35 kilobases (kb)), and contains genes primarily involved in the replication of the plastid genome [10, 11]. In apicomplexans that have a plastid, many of the original plastid genes appear to have been lost (for example, photosynthesis genes) and some genes have been transferred to the host nuclear genome; their proteins are reimported into the apicoplast where they function . Plastids acquired by secondary endosymbiosis are scattered among eukaryotic lineages, including cryptomonads, haptophytes, alveolates, euglenids and chlorarachnions [13–17]. Among the alveolates, plastids are found in dinoflagellates and most examined apicomplexans but not in ciliates. Recent studies on the nuclear-encoded, plastid-targeted glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene suggest a common origin of the secondary plastids in apicomplexans, some dinoflagellates, heterokonts, haptophytes and cryptomonads [8, 18]. If true, this would indicate that the lineage that gave rise to Cryptosporidium contained a plastid, even though many of its descendants (for example, the ciliates) appear to lack a plastid. Although indirect evidence has been noted for the past existence of an apicoplast in C. parvum [19, 20], no rigorous phylogenomic survey for nuclear-encoded genes of plastid or algal origin has been reported.
Gene transfers, either intracellular (IGT) from an endosymbiont or organelle to the host nucleus or horizontal (HGT) between species, can dramatically alter the biochemical repertoire of host organisms and potentially create structural or functional novelties [21–23]. In parasites, genes transferred from prokaryotes or other sources are potential targets for chemotherapy due to their phylogenetic distance or lack of a homolog in the host [24, 25]. The detection of transferred genes in Cryptosporidium is thus of evolutionary and practical importance.
In this study, we use a phylogenomic approach to mine the recently sequenced genome of C. parvum (IOWA isolate; 9.1 megabases (Mb))  for evidence of the past existence of an endosymbiont or apicoplast organelle and of other independent HGTs into this genome. We have detected genes of cyanobacterial/algal origin and genes acquired from other prokaryotic lineages in C. parvum. The fate of several of these transferred genes in C. parvum is explored by expression analyses. The significance of our findings and their impact on the genetic makeup of the parasite are discussed.
Distribution of best non-apicomplexan BLAST hits in searches of the GenBank non-redundant protein database
E < 10-3
E < 10-7
Genes of algal or eubacterial origin in C. parvum
Putative gene name
Kinase; nucleotide metabolism
Hypothetical protein A†
Inosine 5' monophosphate dehydrogenase
Purine nucleotide biosynthesis
Tryptophan synthetase β chain
Amino acid biosynthesis
1,4-α-glucan branching enzyme
1,4-α-glucan branching enzyme
RNA processing and modification
Phosphoric ester hydrolase
ATPase of the AAA class
Energy production and conversion
Amino acid biosynthesis
Conserved hypothetical protein B
Amino acid biosynthesis
Asparaginyl tRNA synthetase†
Biopteridine transporter (BT-1)
Hypothetical protein C† (possible Zn-dependent metalloprotease)
Uridine kinase/uracil phosphoribosyltransferase†
Nucleotide salvage metabolism
Calcium-dependent protein kinases* †
Kinase; cell signal transduction
Genes of cyanobacterial/algal origin
Extant Cryptosporidium species do not contain an apicoplast genome or any physical structure thought to represent an algal endosymbiont or the plastid organelle it contained [6, 7]. The only possible remaining evidence of the past association of an endosymbiont or its cyanobacterially derived plastid organelle might be genes transferred from these genetic sources to the host genome prior to the physical loss of the endosymbiont or organelle itself. Several such genes were identified.
Another C. parvum nuclear-encoded gene of putative cyanobacterial origin is a protein of unknown function belonging to the biopterine transporter family (BT-1) (Table 2). Similarity searches with this protein revealed significant hits to other apicomplexans (for example, P. falciparum, Theileria annulata, T. gondii), plants (Arabidopsis, Oryza), cyanobacteria (Trichodesmium, Nostoc and Synechocystis), a ciliate (Tetrahymena) and the kinetoplastids (Leishmania and Trypanosoma). Arabidopsis thaliana apparently contains at least two copies of this gene; the protein of one (accession number NP_565734) is predicted by ChloroP  to be chloroplast-targeted, suggestive of its plastid derivation. The taxonomic distribution and sequence similarity of this protein with cyanobacterial and chloroplast homologs are also indicative of its affinity to plastids.
Only one gene of algal nuclear origin, glucose-6-phosphate isomerase (G6PI), was identified by the screen described here. Several other algal-like genes are probable, but their support was weaker (Table 2). A 'plant-like' G6PI has been described in other apicomplexan species (P. falciparum, T. gondii ) and a 'cyanobacterial-like' G6PI has been described in the diplomonads Giardia intestinalis and Spironucleus and the parabasalid Trichomonas vaginalis . Figure 4b illustrates these observations nicely. At the base of the tree, the eukaryotic organisms Giardia, Spironucleus and Trichomonas group with the cyanobacterium Nostoc, as previously published. In the midsection of the tree, the G6PI of apicomplexans and ciliates forms a well-supported monophyletic group with the plants and the heterokont Phytophthora. The multiple protein sequence alignment of G6PI identifies several conserved positions shared exclusively by apicomplexans, Tetrahymena, plants and Phytophthora. This gene does not contain a signal or transit peptide and is not predicted to be targeted to the apicoplast in P. falciparum. The remainder of the tree shows a weakly supported branch including eubacteria, fungi and several eukaryotes. The eukaryotes are interrupted by the inclusion of G6PI from the eubacterial organisms Escherichia coli and Cytophaga. This relationship of E. coli G6PI and eukaryotic G6PI has been observed before and may represent yet another gene transfer .
Genes of eubacterial (non-cyanobacterial) origin
Our study identified HGTs from several distinct sources, involving a variety of biochemical activities and metabolic pathways (Table 2). Notably, the nucleotide biosynthesis pathway contains at least two previously published, independently transferred genes from eubacteria. Inosine 5' monophosphate dehydrogenase (IMPDH), an enzyme for purine salvage, was transferred from ε-proteobacteria . Another enzyme involved in pyrimidine salvage, thymidine kinase (TK), is of α or γ-proteobacterial ancestry .
Another gene of eubacterial origin identified in C. parvum is tryptophan synthetase β subunit (trpB). This gene has been identified in both C. parvum and C. hominis, but not in other apicomplexans. The relationship of C. parvum trpB to proteobacterial sequences is well-supported as a monophyletic group by two of the three methods used in our analyses (Figure 4c).
Mode of acquisition
We examined the transferred genes for evidence of non-independent acquisition, for example, blocks of transferred genes or evidence that genes were acquired together from the same source. Examination of the chromosomal location of the genes listed in Table 2 demonstrates that the genes are currently located on different chromosomes and in most cases do not appear to have been transferred or retained in large blocks. There are two exceptions. The trpB gene and the gene for aspartate ammonia ligase are located 4,881 base-pairs (bp) apart on the same strand of a contig for chromosome V; there is no annotated gene between these two genes. Both genes are of eubacterial origin and are not found in other apicomplexan organisms. While it is possible that they have been acquired independently with this positioning, or later came to have this positioning via genome rearrangements, it is interesting to speculate that these genes were acquired together. The origin of trpB is proteobacterial. The origin of aspartate ammonia ligase is eubacterial, but not definitively of any particular lineage. In the absence of genome sequences for all organisms, throughout all of time, exact donors are extremely difficult to assess and inferences must be drawn from sequences that appear to be closely related to the actual donor.
In the second case, C. parvum encodes two genes for 1,4-α-glucan branching enzymes. Both are eubacterial in origin and both are located on chromosome VI, although not close together. They are approximately 110 kb apart and many intervening genes are present. The evidence that these genes were acquired together comes from the phylogenetic analysis presented in Figure 5. The duplication that gave rise to the two 1,4-α-glucan branching enzymes is old, and is well supported by the tree shown in Figure 5. A number of eubacteria (11), including cyanobacteria, contain this duplication. The 1,4-α-glucan branching enzymes of C. parvum and T. gondii represent one copy each of this ancient duplication. This suggests that the ancestor of C. parvum and T. gondii acquired the genes after they had duplicated and diverged in eubacteria.
Expression of transferred genes
A small published collection of 567 EST sequences for C. parvum is also available. These ESTs were searched with each of the 31 candidate genes surviving the phylogenomic screen. Three genes - aspartate ammonia ligase, BT-1 and lactate dehydrogenase - are expressed, as confirmed by the presence of an EST (Table 2).
A genome-wide search for intracellular and horizontal gene transfers in C. parvum was carried out. We systematically determined the evolutionary origins of genes in the genome using phylogenetic approaches, and further confirmed the existence and expression of putatively transferred genes with laboratory experiments. The methodology adopted in this study provides a broad picture of the extent and the importance of gene transfer in apicomplexan evolution.
The identification of gene transfers is often subject to errors introduced by methodology, data quality and taxonomic sampling. The phylogenetic approach adopted in this study is preferable to similarity searches [33, 34] but several factors, including long-branch attraction, mutational saturation, lineage-specific gene loss and acquisition, and incorrect identification of orthologs, can distort the topology of a gene tree [35, 36]. Incompleteness in the taxonomic record may also lead to false positives for IGT and HGT identification. In our study, we have attempted to alleviate these factors, as best as is possible, by sampling the GenBank nonredundant protein database, dbEST and organism-specific databases and by using several phylogenetic methods. Still, these issues remain a concern for this study as the taxonomic diversity of unicellular eukaryotes is vastly undersampled and studies are almost entirely skewed towards parasitic organisms.
The published analysis of the C. parvum genome sequence identified 14 bacteria-like and 15 plant-like genes based on similarity searches . Six of these bacterial-like and three plant-like genes were also identified as probable transferred genes in the phylogenomic analyses presented here. We have examined the fate of genes identified by one analysis and not the other to uncover the origin of the discrepancy. First, methodology is the single largest contributing factor. Genes with bacterial-like or plant-like BLAST similarities which, from the phylogenetic analyses, do not appear to be transfers were caused by the fact that PyPhy was unable to generate trees due to an insufficient number of significant hits in the database, or because of the stringent coverage length and similarity requirements adopted in this analysis. Only seven of the previously identified 15 plant-like and 11 of 14 eubacterial-like genes survived the predefined criteria for tree construction. Second, subsequent phylogenetic analyses including additional sequences from non-GenBank databases failed to provide sufficient evidence or significant support for either plant or eubacterial ancestry. Third, searches of dbEST and other organism-specific databases yielded other non-plant or non-eubacterial organisms as nearest neighbors, thus removing the possibility of a transfer.
The limitations of similarity searches and incomplete taxonomic sampling are well evidenced in our phylogenomic analyses. From similarity searches, C. parvum, like P. falciparum , is more similar to the plants Arabidopsis and Oryza than to any other single organism. Almost 800 predicted genes have best non-apicomplexan BLAST hits of at least 10-7 to plants and eubacteria (Table 1). Yet only 31 can be inferred to be transferred genes at this time with the datasets and methodology available (Table 2). In many cases (for example, phosphoglucomutase) the C. parvum gene groups phylogenetically with plant and bacterial homologs, but with only modest support. In other cases, such as pyruvate kinase and the bi-functional dehydrogenase enzyme (AdhE), gene trees obtained from automated PyPhy analyses indicate a strong monophyletic grouping of the C. parvum gene with plant or eubacterial homologs, but this topology disappears when sequences from other unicellular eukaryotes, such as Dictyostelium, Entamoeba and Trichomonas are included in the analysis (data not shown).
The list of genes in Table 2 should be considered a current best estimate of the IGTs and HGTs in C. parvum instead of a definitive list. As genomic data are obtained from a greater diversity of unicellular eukaryotes and eubacteria, phylogenetic analyses of nearest neighbors are likely to change.
Did Cryptosporidiumcontain an endosymbiont or plastid organelle?
The C. parvum sequences of cyanobacterial and algal origin reported here had to enter the genome at some point during its evolution. Formal possibilities include vertical inheritance from a plastid-containing chromalveolate ancestor, HGT from the cyanobacterial and algal sources (or from a secondary source such as a plastid-containing apicomplexan), or IGT from an endosymbiont/plastid organelle during evolution, followed by loss of the source. Cryptosporidium does not harbor an apicoplast organelle or any trace of a plastid genome ; thus an IGT scenario would necessitate loss of the organelle in Cryptosporidium or the lineage giving rise to it. The exact position of C. parvum on the tree of life has been debated, with developmental and morphological considerations placing it within the Apicomplexa, and molecular analyses locating it in various positions, both within and outside the Apicomplexa , but primarily within. If we assume that C. parvum is an apicomplexan, and if the secondary endosymbiosis which is believed to have given rise to the apicoplast occurred before the formation of the Apicomplexa, as has been suggested , C. parvum would have evolved from a plastid-containing lineage and would be expected to harbor traces of this relationship in its nuclear genome. Genes of likely cyanobacterial and algal/plant origin are detected in the nuclear genome of C. parvum (Table 2) and thus IGT followed by organelle loss cannot be ruled out.
What about other interpretations? While it is formally possible that these genes were acquired independently via HGT in C. parvum, their shared presence in other alveolates (including the non-plastidic ciliate Tetrahymena) provides the best evidence against this scenario as multiple independent transfers would be required and so far there is no evidence for intra-alveolate gene transfer. Vertical inheritance is more difficult to address as it involves distinguishing between genes acquired via IGT from a primary endosymbiotic event versus a secondary endosymbioic event. Our data, especially the analysis of G6PI and BT-1 are consistent with both primary and secondary endosymbioses, provided that the secondary endosymbiosis is pre-alveolate in origin. As more genome data become available and flanking genes can be examined for each gene in a larger context, positional information will be informative in distinguishing among the alternatives.
The plastidic nature of some genes is particularly apparent. There is a shared indel among leucine aminopeptidase protein sequences in apicomplexans, cyanobacteria and plant chloroplast precursors (Figure 3). The C. parvum leucine aminopeptidase does contain an amino-terminal extension of approximately 85-65 amino acids (depending on the alignment) relative to bacterial homologs, but this extension does not contain a signal sequence. The extension in P. falciparum is 85 amino acids and the protein is believed to be targeted to the apicoplast [26, 37]. No similarity is detected between the C. parvum and P. falciparum amino-terminal extensions (data not shown).
Other genes were less informative in this analysis. Among these, aldolase was reported in both P. falciparum  and the kinetoplastid parasite Trypanosoma  as a plant-like gene. The protein sequences of aldolase are similar in C. parvum and P. falciparum, with an identity of 60%. In our phylogenetic analyses, C. parvum clearly forms a monophyletic group with Plasmodium, Toxoplasma and Eimeria. This branch groups with Dictyostelium, Kinetoplastida and cyanobacterial lineages, but bootstrap support is not significant. The sister group to the above organisms are the plants and additional cyanobacteria, but again with no bootstrap support (see Additional data file 1 for phylogenetic tree). Another gene, enolase, contains two indels shared between land plants and apicomplexans (including C. parvum) and was suggested to be a plant-like gene , but alternative explanations exist .
The biochemical activity of the polyamine biosynthetic enzyme arginine decarboxylase (ADC), which is typically found in plants and bacteria, was previously reported in C. parvum . However, we were unable to confirm its presence by similarity searches of the two Cryptosporidium genome sequences deposited in CryptoDB using plant (Cucumis sativa, GenBank accession number AAP36992), cyanobacterial (Nostoc sp., NP-487441; Synechocystis sp., NP-439907) and other bacterial (Yersinia pestis, NP-404547) homologs.
A plethora of prokaryotic genes
Several HGTs from bacteria have been reported previously in C. parvum [25, 32, 40]. We detected many more in our screen of the completed C. parvum genome sequence (Table 2). In most cases, the exact donors of these transferred genes were difficult to determine. However, for those genes whose donors could be more reliably inferred (Table 2), several appear to be from different sources and hence represent independent transfer events. In one compelling case, both the trpB and aspartate ammonia ligase genes are located 4,881 bp apart on the same strand of a contig for chromosome V and there is no gene separating them. Both genes are of eubacterial origin and neither gene is detected in other apicomplexans. In addition, the aspartate ammonia ligase gene is expressed, as evidenced by an EST. In another case, copies of a 1,4-α-glucan branching enzyme gene duplication pair that is present in many eubacteria, were detected on the same chromosome in C. parvum. C. parvum also contains many transferred genes from distinct eubacterial sources that are not present in other apicomplexans (for example, IMPDH, TK (thymidine kinase), trpB and the gene for aspartate ammonia ligase).
The endosymbiotic event that gave rise to the mitochondrion occurred very early in eukaryotic evolution and is associated with significant IGT. However, most of these transfer events happened long before the evolutionary time window we explored in this study . Many IGTs from the mitochondrial genome that have been retained are almost universally present in eukaryotes (including C. parvum which does not contain a typical mitochondrion [7, 42–44]) and thus would not be detected in a PyPhy screen since the 'nearest phylogenetic neighbor' on the tree would be taxonomically correct and not appear as a relationship indicative of a gene transfer.
The impact of gene transfers on host evolution
Gene transfer is an important evolutionary force [21, 22, 45, 46]. Several of the transferred genes identified in C. parvum are known to be expressed. IMPDH has been shown to be essential in C. parvum purine metabolism  and TK has been shown to be functional in pyrimidine salvage . It is not yet clear whether these genes were acquired independently in this lineage, or have been lost from the rest of the apicomplexan lineage, or whether both these have happened. However, it is clear that their presence has facilitated the remodeling of nucleotide biosynthesis. C. parvum no longer possesses the ability to synthesize nucleotides; instead it relies entirely on salvage.
Many apicoplast and algal nuclear genes have been transferred to the host nuclear genome, where they were subsequently translated in the cytosol and their proteins targeted to the apicoplast organelle. However, as there is no apicoplast in C. parvum, acquired plastidic proteins are theoretically destined to go elsewhere. In the absence of an apicoplast, it is tempting to suspect that plastid-targeted proteins would have been lost, or would be detected as pseudogenes. No identifiable pseudogenes were detected and at least one gene is still viable. The C. parvum leucine aminopeptidase, which still contains an amino-terminal extension (without a signal peptide), is intact and is expressed, as shown in Figure 6. None of the cyanobacterial/algal genes identified in our study contains a canonical presequence for apicoplast targeting. One exception to this is phosphoglucomutase, a gene not present in Table 2 because of its poorly supported relationships in phylogenetic analyses. This gene exists in two copies as a tandem duplication in the C. parvum genome. One copy has a long amino-terminal extension (97 amino acids) beginning with a signal peptide. The extension does not contain characteristics of a transit peptide. Expression of a fluorescent reporter construct containing this extension in a related parasite, T. gondii, did not reveal apicoplast targeting but instead secretion via dense granules (see Additional data file 1). Exactly how and where intracellularly transferred genes (especially those that normally target the apicoplast) have become incorporated into other metabolic processes remains a fertile area for exploration.
Cryptosporidium is the recipient of a large number (31) of transferred genes, many of which are not shared by other apicomplexan parasites. The genes have been acquired from several different sources including α-, β-, and ε-proteobacteria, cyanobacteria, algae/plants and possibly the Archaea. We have described two cases of two genes that appear to have been acquired together from a eubacterial source: trpB and the aspartate ammonia ligase gene are located within 5 kb of each other, while the two copies of 1,4-α-glucan branching enzyme represent copies of an ancient gene duplication also observed in cyanobacteria.
Once thought to be a relatively rare event, reports of gene transfers in eukaryotes are increasingly common. The abundance of available eukaryotic genome sequence is providing the material for analyses that were not possible only a few years ago. Analysis of the Arabidopsis genome  has revealed potentially thousands of genes that were transferred intracellularly. HGTs are still a relatively rare class of genes among multicellular eukaryotes, most probably because of the segregation of the germ line. By definition, unicellular eukaryotes do not have a separate germ line and are obligated to tolerate the acquisition of foreign genes if they are to survive. Among unicellular eukaryotes, there are now many reports of HGTs: Giardia [48, 49], Trypanosoma , Entamoeba [21, 49], Euglena , Cryptosporidium [25, 32, 40] and other apicomplexans .
As discussed earlier, genes transferred from distant phylogenetic sources such as eubacteria could be potential therapeutic targets. In apicomplexans, transferred genes are already some of the most promising targets of anti-parasitic drugs and vaccines [7, 25, 52]. We have shown that several transferred genes are differentially expressed in the C. parvum genome, and in two cases (IMPDH and TK), the transferred genes have been shown to be functional [25, 32]. The successful integration, expression and survival of transferred genes in the Cryptosporidium genome has changed the genetic and metabolic repertoire of the parasite.
Materials and methods
Cryptosporidium sequence sources
Genomic sequences for C. parvum and C. hominis were downloaded from CryptoDB . Genes were predicted for the completed C. parvum (IOWA) sequence as previously described using the Glimmer program  trained on Cryptosporidium coding sequences . A few predicted genes that demonstrated apparent sequence incompleteness were reconstructed from genomic sequence by comparison with apicomplexan orthologs. The predicted protein encoding data set contained 5,519 sequences. A comparison of this gene set to the published annotation revealed that the Glimmer-predicted gene set contained all but 40 of the 3,396 annotated protein-encoding sequences deposited in GenBank. These 40 were added to our dataset and analyzed. Glimmer does not predict introns and some introns are present in the genome [7, 20]; thus our gene count is artificially inflated. Likewise, the official C. parvum annotation did not consider ORFs of less than 100 amino acids that did not have significant BLAST hits and thus may be a slight underestimate .
An internal database (ApiDB) containing all available apicomplexan sequence data was created . A second BLAST-searchable database, PyPhynr, was constructed that included SwissProt, TrEMBL and TrEMBL_new, as released in August 2003, predicted genes from C. parvum, ORFs of more than 120 amino acids from Theileria annulata, and more than 75 amino acids from consensus ESTs for several apicomplexan organisms. Genomic sequences for T. gondii (8x coverage) and clustered ESTs were downloaded from ToxoDB [55, 56]. Genomic data were provided by The Institute for Genomic Research (TIGR), and by the Sanger Institute. EST sequences were generated by Washington University. In addition, this study used sequence data from several general and species-specific databases. Specifically, the NCBI GenBank nr and dbEST were downloaded  and extensively searched. To provide taxonomic completeness, additional genes were obtained via searches of additional databases including: Entamoeba histolytica , D. discoideum , the kinetoplastids Leishmania major , T. brucei , T. cruzi , and a ciliate Tetrahymena thermophila . Sequence data for T. annulata, E. histolytica, D. discoideum, L. major and T. brucei were produced by the Pathogen Sequencing Unit of the Sanger Institute and can be obtained from . Preliminary sequence data for T. thermophila was obtained from TIGR and can be accessed at .
Phylogenomic analyses and similarity searches
The source code of the phylogenomic software PyPhy  was kindly provided by Thomas Sicheritz-Ponten and modified to include analyses of eukaryotic groups, and changes to improve functionality . For initial phylogenomic analyses, a BLAST cutoff of 60% sequence length coverage and 50% sequence similarity was adopted and the neighbor-joining program of PAUP 4.0b10 for Unix  was used. A detailed description of our phylogenomic pipeline and PyPhy implementation are described  and outlined in Figure 1.
Output gene trees with phylogenetic connections (that is, the nearest non-self neighbors at a distinct taxonomic rank)  to prokaryotes and algae-related groups were manually inspected. As the trees are unrooted, several factors were considered in the screen for candidate transferred genes. If the C. parvum gene does not form a monophyletic group with prokaryotic or plant-related taxa regardless of rooting, the subject gene was eliminated from further consideration. If the topology of the gene tree is consistent with a phylogenetic anomaly caused by gene transfer, but may also be interpreted differently if the tree is rooted otherwise, it was removed from consideration at this time. If the top hits of both nr and dbEST database searches are predominantly non-plant eukaryotes, and the topology of the tree was poor, the subject gene was considered an unlikely candidate. Finally, all 551 protein sequences predicted to be NEAPs in the malarial parasite P. falciparum  were used to search the C. parvum genome and the results were screened using a BLAST cutoff E-value of 10-5 and a length coverage of 50%. Sequences identified by these searches were added to the candidate list (if not already present) for manual phylogenetic analyses to verify their likely origins. It should be noted that all trees were screened for the existence of a particular phylogenetic relationship. In some cases the proteins utilized to generate a particular tree are capable of resolving relationships among many branches of the tree of life, and in others they are not. Despite these differences in resolving power, the proteins which survive our phylogenetic screen and subsequent detailed analyses described below exhibit significant support for the branches of the tree in which we are interested. Similar procedures were used to characterize the complement of nuclear-encoded genes of plastid origin in the Arabidopsis genome . BLAST searches were performed on GenBank releases 138-140 .
Detailed phylogenetic analyses of candidate genes identified by phylogenomic screening: candidate genes surviving the PyPhy phylogenomic screen were reanalyzed with careful attention to taxonomic completeness, including representative species from major prokaryotic and eukaryotic lineages when necessary and possible. New multiple sequence alignments were created with ClustalX , followed by manual refinement. Only unambiguously aligned sequence segments were used for subsequent analyses (see Additional data file 1). Phylogenetic analyses were performed with a maximum likelihood method using TREE-PUZZLE version 5.1 for Unix , a distance method using the program neighbor of PHYLIP version 3.6a package , and a maximum parsimony method with random stepwise addition using PAUP* 4.0b10 . Bootstrap support was estimated using 1,000 replicates for both parsimony and distance analyses and quartet puzzling values were obtained using 10,000 puzzling steps for maximum likelihood analyses. Distance calculation used the Jones-Taylor-Thornton (JTT) substitution matrix , and site-substitution variation was modeled with a gamma-distribution whose shape parameter was estimated from the data. For maximum likelihood analyses, a mixed model of eight gamma-distributed rates and one invariable rate was used to calculate the pairwise maximum likelihood distances. The unrooted trees presented in Figures 4 and 5 were drawn by supplying TREE-PUZZLE with the maximum parsimony tree and using TREE-PUZZLE distances as described above to calculate the branch lengths. The trees were visualized and prepared for publication with TreeView X Version 0.4.1 .
Genomic Southern analysis
C. parvum (IOWA) oocysts (108) were obtained from the Sterling Parasitology Laboratory at the University of Arizona and were lysed using a freeze/thaw method. Genomic DNA was purified using the DNeasy Tissue Kit (Qiagen). Genomic DNA (5 μg) was restricted with BamH1 and EcoR1 respectively and electrophoresed on a 0.8% gel in 1x TAE buffer, transferred to a positively charged nylon membrane (Bio-Rad), and fixed using a UVP crosslinker set at 125 mJ as described in . C. parvum genomic DNA for the probes (700-1,500 bp) was amplified by PCR (see Additional data file 1).
Semi-quantitative reverse transcription-PCR
Sterilized C. parvum (IOWA isolate) oocysts were used to infect confluent human adenocarcinoma cell monolayers at a concentration of one oocyst per cell as previously described . Total RNA was prepared from mock-infected and C. parvum-infected HCT-8 cultures at 2, 6, 12, 24, 36, 48 and 72 h post-inoculation by directly lysing the cells with 4 ml TRIzol reagent (GIBCO-BRL/Life Technologies). Purified RNA was resuspended in RNAse-free water and the integrity of the samples was confirmed by gel electrophoresis.
Primers specific for several transferred genes identified in the study were designed (see Additional data file 1) and a semi-quantitative RT-PCR analysis was carried out as previously described . Primers specific for C. parvum 18S rRNA were used to normalize the amount of cDNA product of the candidate gene to that of C. parvum rRNA in the same sample. PCR products were separated on a 4% non-denaturing polyacrylamide gel and signals from specific products were captured and quantified using a phosphorimaging system (Molecular Dynamics). The expression level of each gene at each time point was calculated as the ratio of its RT-PCR product signal to that of the C. parvum 18S rRNA. Six independent time-course experiments were used in the analysis.
Additional data files
Additional data is provided with the online version of this paper, consisting of a PDF file (Additional data file 1) containing: materials and methods for genomic Southern analysis; the amino-acid sequences of genes listed in Table 2; accession numbers for sequences used in Figure 4; accession numbers for sequences used in Figure 5; expression of C. parvum phosphoglucomutase in T. gondii; table of primers used for RT-PCR experiments; phylogenetic tree of aldolase; alignment files for phylogenetic analyses in Figure 4; and the alignment of 1,4-α-glucan branching enzyme sequences used in Figure 5.
We thank G. Buck (Virginia Commonwealth University), G. Widmer and S. Tzipori (Tufts University) for access to C. hominis genotype I genomic sequence data. Genomic sequence for Toxoplasma gondii (8X coverage) and clustered ESTs were downloaded from ToxoDB . Genomic data were provided by The Institute for Genomic Research (supported by the NIH grant #AI05093), and by the Sanger Institute (Wellcome Trust). Apicomplexan EST sequences were generated by Washington University (NIH grant #1R01AI045806-01A1). Occasionally, genes were obtained via searches of several databases containing Entamoeba histolytica , Dictyostelium discoideum , kinetoplastids Leishmania major , Trypanosoma brucei , T. cruzi , and ciliate Tetrahymena thermophila . Sequence data for T. annulata, E. histolytica, D. discoideum, Leishmania major and T. brucei were produced by the Pathogen Sequencing Unit of the Sanger Institute and can be obtained from . Preliminary sequence data for Tetrahymena thermophila was obtained from The Institute for Genomic Research and can be accessed at . We thank Fallon Hampton for her work on the phosphoglucomutase expression constructs. She was supported by a summer undergraduate fellowship in genetics (SUNFIG) award. This study was funded by a research grant from the University of Georgia Research Foundation to J.C.K. and NIH grant U01 AI 46397 to M.S.A. J.H. is supported by a postdoctoral fellowship from the American Heart Association. We thank Boris Striepen, Marc-Jan Gubbels and three anonymous reviewers for comments that greatly increased the clarity and precision of the analyses in the manuscript.
- Mead JR: Cryptosporidiosis and the challenges of chemotherapy. Drug Resistance Updates. 2002, 5: 47-57. 10.1016/S1368-7646(02)00011-0.PubMedView ArticleGoogle Scholar
- CDC: bioterrorism agents/diseases. [http://www.bt.cdc.gov/agent/agentlist.asp]
- Zhu G, Keithly JS, Philippe H: What is the phylogenetic position of Cryptosporidium?. Int J Syst Evol Microbiol. 2000, 50: 1673-1681.PubMedView ArticleGoogle Scholar
- Leander BS, Clopton RE, Keeling PJ: Phylogeny of gregarines (Apicomplexa) as inferred from small-subunit rDNA and beta-tubulin. Int J Syst Evol Microbiol. 2003, 53: 345-354. 10.1099/ijs.0.02284-0.PubMedView ArticleGoogle Scholar
- Carreno RA, Martin DS, Barta JR: Cryptosporidium is more closely related to the gregarines than to coccidia as shown by phylogenetic analysis of apicomplexan parasites inferred using small-subunit ribosomal RNA gene sequences. Parasitol Res. 1999, 85: 899-904. 10.1007/s004360050655.PubMedView ArticleGoogle Scholar
- Zhu G, Marchewka MJ, Keithly JS: Cryptosporidium parvum appears to lack a plastid genome. Microbiology. 2000, 146: 315-321.PubMedView ArticleGoogle Scholar
- Abrahamsen MS, Templeton TJ, Enomoto S, Abrahante JE, Zhu G, Lancto CA, Deng M, Liu C, Widmer G, Tzipori S, et al: Complete genome sequence of the apicomplexan, Cryptosporidium parvum. Science. 2004, 304: 441-445. 10.1126/science.1094786.PubMedView ArticleGoogle Scholar
- Fast NM, Kissinger JC, Roos DS, Keeling PJ: Nuclear-encoded, plastid-targeted genes suggest a single common origin for apicomplexan and dinoflagellate plastids. Mol Biol Evol. 2001, 18: 418-426.PubMedView ArticleGoogle Scholar
- Kohler S, Delwiche CF, Denny PW, Tilney LG, Webster P, Wilson RJ, Palmer JD, Roos DS: A plastid of probable green algal origin in apicomplexan parasites. Science. 1997, 275: 1485-1489. 10.1126/science.275.5305.1485.PubMedView ArticleGoogle Scholar
- Wilson RJ, Denny PW, Preiser PR, Rangachari K, Roberts K, Roy A, Whyte A, Strath M, Moore DJ, Moore PW, et al: Complete gene map of the plastid-like DNA of the malaria parasite Plasmodium falciparum. J Mol Biol. 1996, 261: 155-172. 10.1006/jmbi.1996.0449.PubMedView ArticleGoogle Scholar
- Denny PW, Preiser PR, Williams DC, Wilson I: Evidence for a single origin of the 35 kb plastid DNA in apicomplexans. Protist . 1998, 149: 51-59.PubMedView ArticleGoogle Scholar
- Roos DS, Crawford MJ, Donald RG, Fraunholz M, Harb OS, He CY, Kissinger JC, Shaw MK, Striepen B: Mining the Plasmodium genome database to define organellar function: what does the apicoplast do?. Philos Trans R Soc Lond B Biol Sci. 2002, 357: 35-46. 10.1098/rstb.2001.1047.PubMedPubMed CentralView ArticleGoogle Scholar
- Cavalier-Smith T: Membrane heredity and early chloroplast evolution. Trends Plant Sci. 2000, 5: 174-182. 10.1016/S1360-1385(00)01598-3.PubMedView ArticleGoogle Scholar
- Delwiche CF: Tracing the thread of plastid diversity through the tapestry of life. Am Nat. 1999, 154: S164-S177. 10.1086/303291.PubMedView ArticleGoogle Scholar
- Maier UG, Douglas SE, Cavalier-Smith T: The nucleomorph genomes of cryptophytes and chlorarachniophytes. Protist. 2000, 151: 103-109.PubMedView ArticleGoogle Scholar
- Yoon HS, Hackett JD, Pinto G, Bhattacharya D: The single, ancient origin of chromist plastids. Proc Natl Acad Sci USA. 2002, 99: 15507-15512. 10.1073/pnas.242379899.PubMedPubMed CentralView ArticleGoogle Scholar
- Bhattacharya D, Yoon HS, Hackett JD: Photosynthetic eukaryotes unite: endosymbiosis connects the dots. BioEssays. 2004, 26: 50-60. 10.1002/bies.10376.PubMedView ArticleGoogle Scholar
- Harper JT, Keeling PJ: Nucleus-encoded, plastid-targeted glyceraldehyde-3-phosphate dehydrogenase (GAPDH) indicates a single origin for chromalveolate plastids. Mol Biol Evol. 2003, 20: 1730-1735. 10.1093/molbev/msg195.PubMedView ArticleGoogle Scholar
- Keithly JS, Zhu G, Upton SJ, Woods KM, Martinez MP, Yarlett N: Polyamine biosynthesis in Cryptosporidium parvum and its implications for chemotherapy. Mol Biochem Parasitol. 1997, 88: 35-42. 10.1016/S0166-6851(97)00063-7.PubMedView ArticleGoogle Scholar
- Bankier AT, Spriggs HF, Fartmann B, Konfortov BA, Madera M, Vogel C, Teichmann SA, Ivens A, Dear PH: Integrated mapping, chromosomal sequencing and sequence analysis of Cryptosporidium parvum. Genome Res. 2003, 13: 1787-1799.PubMedPubMed CentralGoogle Scholar
- Van Der Giezen M, Cox S, Tovar J: The iron-sulfur cluster assembly genes iscS and iscU of Entamoeba histolytica were acquired by horizontal gene transfer. BMC Evol Biol. 2004, 4: 7-10.1186/1471-2148-4-7.PubMedPubMed CentralView ArticleGoogle Scholar
- Gojkovic Z, Knecht W, Zameitat E, Warneboldt J, Coutelis JB, Pynyaha Y, Neuveglise C, Moller K, Loffler M, Piskur J: Horizontal gene transfer promoted evolution of the ability to propagate under anaerobic conditions in yeasts. Mol Genet Genomics. 2004, 271: 387-393. 10.1007/s00438-004-0995-7.PubMedView ArticleGoogle Scholar
- Andersson JO, Doolittle WF, Nesbo CL: Genomics. Are there bugs in our genome?. Science. 2001, 292: 1848-1850. 10.1126/science.1062241.PubMedView ArticleGoogle Scholar
- Kissinger J, Crawford MJ, Roos D, Ajioka JW: Toxoplasma gondii: A model for evolutionary genomics and chemotherapy. In Pathogen Genomics: Impact on Human Health. Edited by: Shaw KJ. 2001, Totowa NJ: Humana Press, 255-279.Google Scholar
- Striepen B, Pruijssers AJP, Huang J, Li C, Gubbels MJ, Umejiego NN, Hedstrom L, Kissinger JC: Gene transfer in the evolution of parasite nucleotide biosynthesis. Proc Natl Acad Sci USA. 2004, 101: 3154-3159. 10.1073/pnas.0304686101.PubMedPubMed CentralView ArticleGoogle Scholar
- Gardner MJ, Hall N, Fung E, White O, Berriman M, Hyman RW, Carlton JM, Pain A, Nelson KE, Bowman S, et al: Genome sequence of the human malaria parasite Plasmodium falciparum. Nature. 2002, 419: 498-511. 10.1038/nature01097.PubMedView ArticleGoogle Scholar
- Sicheritz-Ponten T, Andersson SG: A phylogenomic approach to microbial evolution. Nucleic Acids Res. 2001, 29: 545-552. 10.1093/nar/29.2.545.PubMedPubMed CentralView ArticleGoogle Scholar
- Emanuelsson O, Nielsen H, von Heijne G: ChloroP, a neural network-based method for predicting chloroplast transit peptides and their cleavage sites. Protein Sci. 1999, 8: 978-984.PubMedPubMed CentralView ArticleGoogle Scholar
- Dzierszinski F, Popescu O, Toursel C, Slomianny C, Yahiaoui B, Tomavo S: The protozoan parasite Toxoplasma gondii expresses two functional plant-like glycolytic enzymes. Implications for evolutionary origin of apicomplexans. J Biol Chem. 1999, 274: 24888-24895. 10.1074/jbc.274.35.24888.PubMedView ArticleGoogle Scholar
- Henze K, Horner DS, Suguri S, Moore DV, Sanchez LB, Muller M, Embley TM: Unique phylogenetic relationships of glucokinase and glucosephosphate isomerase of the amitochondriate eukaryotes Giardia intestinalis, Spironucleus barkhanus and Trichomonas vaginalis. Gene. 2001, 281: 123-131. 10.1016/S0378-1119(01)00773-9.PubMedView ArticleGoogle Scholar
- Rudolph B, Hansen T, Schonheit P: Glucose-6-phosphate isomerase from the hyperthermophilic archaeon Methanococcus jannaschii: characterization of the first archaeal member of the phosphoglucose isomerase superfamily. Arch Microbiol. 2004, 181: 82-87. 10.1007/s00203-003-0626-4.PubMedView ArticleGoogle Scholar
- Striepen B, White MW, Li C, Guerini MN, Malik SB, Logsdon JM, Liu C, Abrahamsen MS: Genetic complementation in apicomplexan parasites. Proc Natl Acad Sci USA. 2002, 99: 6304-6309. 10.1073/pnas.092525699.PubMedPubMed CentralView ArticleGoogle Scholar
- Genereux DP, Logsdon JM: Much ado about bacteria-to-vertebrate lateral gene transfer. Trends Genet. 2003, 19: 191-195. 10.1016/S0168-9525(03)00055-6.PubMedView ArticleGoogle Scholar
- Stanhope MJ, Lupas A, Italia MJ, Koretke KK, Volker C, Brown JR: Phylogenetic analyses do not support horizontal gene transfers from bacteria to vertebrates. Nature. 2001, 411: 940-944. 10.1038/35082058.PubMedView ArticleGoogle Scholar
- Philippe H, Laurent J: How good are deep phylogenetic trees?. Curr Opin Genet Dev. 1998, 8: 616-623. 10.1016/S0959-437X(98)80028-2.PubMedView ArticleGoogle Scholar
- Kurland CG, Canback B, Berg OG: Horizontal gene transfer: a critical view. Proc Natl Acad Sci USA. 2003, 100: 9658-9662. 10.1073/pnas.1632870100.PubMedPubMed CentralView ArticleGoogle Scholar
- Bahl A, Brunk B, Crabtree J, Fraunholz MJ, Gajria B, Grant GR, Ginsburg H, Gupta D, Kissinger JC, Labo P, et al: PlasmoDB: the Plasmodium genome resource. A database integrating experimental and computational data. Nucleic Acids Res . 2003, 31: 212-215. 10.1093/nar/gkg081.PubMedPubMed CentralView ArticleGoogle Scholar
- Hannaert V, Saavedra E, Duffieux F, Szikora JP, Rigden DJ, Michels PA, Opperdoes FR: Plant-like traits associated with metabolism of Trypanosoma parasites. Proc Natl Acad Sci USA. 2003, 100: 1067-1071. 10.1073/pnas.0335769100.PubMedPubMed CentralView ArticleGoogle Scholar
- Keeling PJ, Palmer JD: Lateral transfer at the gene and subgenic levels in the evolution of eukaryotic enolase. Proc Natl Acad Sci USA. 2001, 98: 10745-10750. 10.1073/pnas.191337098.PubMedPubMed CentralView ArticleGoogle Scholar
- Madern D, Cai X, Abrahamsen MS, Zhu G: Evolution of Cryptosporidium parvum lactate dehydrogenase from malate dehydrogenase by a very recent event of gene duplication. Mol Biol Evol. 2004, 21: 489-497. 10.1093/molbev/msh042.PubMedView ArticleGoogle Scholar
- Baldauf SL: The deep roots of eukaryotes. Science. 2003, 300: 1703-1706. 10.1126/science.1085544.PubMedView ArticleGoogle Scholar
- LaGier MJ, Tachezy J, Stejskal F, Kutisova K, Keithly JS: Mitochondrial-type iron-sulfur cluster biosynthesis genes (IscS and IscU) in the apicomplexan Cryptosporidium parvum. Microbiology. 2003, 149: 3519-3530. 10.1099/mic.0.26365-0.PubMedView ArticleGoogle Scholar
- Riordan CE, Langreth SG, Sanchez LB, Kayser O, Keithly JS: Preliminary evidence for a mitochondrion in Cryptosporidium parvum : phylogenetic and therapeutic implications. J Eukaryot Microbiol. 1999, 46: 52S-55S.PubMedGoogle Scholar
- Riordan CE, Ault JG, Langreth SG, Keithly JS: Cryptosporidium parvum Cpn60 targets a relict organelle. Curr Genet. 2003, 44: 138-147. 10.1007/s00294-003-0432-1.PubMedView ArticleGoogle Scholar
- Martin W: Gene transfer from organelles to the nucleus: frequent and in big chunks. Proc Natl Acad Sci USA. 2003, 100: 8612-8614. 10.1073/pnas.1633606100.PubMedPubMed CentralView ArticleGoogle Scholar
- Ochman H, Lawrence JG, Groisman EA: Lateral gene transfer and the nature of bacterial innovation. Nature. 2000, 405: 299-304. 10.1038/35012500.PubMedView ArticleGoogle Scholar
- Martin W, Rujan T, Richly E, Hansen A, Cornelsen S, Lins T, Leister D, Stoebe B, Hasegawa M, Penny D: Evolutionary analysis of Arabidopsis, cyanobacterial, and chloroplast genomes reveals plastid phylogeny and thousands of cyanobacterial genes in the nucleus. Proc Natl Acad Sci USA. 2002, 99: 12246-12251. 10.1073/pnas.182432999.PubMedPubMed CentralView ArticleGoogle Scholar
- Andersson JO, Sjogren AM, Davis LA, Embley TM, Roger AJ: Phylogenetic analyses of diplomonad genes reveal frequent lateral gene transfers affecting eukaryotes. Curr Biol. 2003, 13: 94-104. 10.1016/S0960-9822(03)00003-4.PubMedView ArticleGoogle Scholar
- Nixon JE, Wang A, Field J, Morrison HG, McArthur AG, Sogin ML, Loftus BJ, Samuelson J: Evidence for lateral transfer of genes encoding ferredoxins, nitroreductases, NADH oxidase, and alcohol dehydrogenase 3 from anaerobic prokaryotes to Giardia lamblia and Entamoeba histolytica. Eukaryot Cell . 2002, 1: 181-190. 10.1128/EC.1.2.181-190.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Henze K, Badr A, Wettern M, Cerff R, Martin W: A nuclear gene of eubacterial origin in Euglena gracilis reflects cryptic endosymbioses during protist evolution. Proc Natl Acad Sci USA. 1995, 92: 9122-9126.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang J, Mullapudi N, Sicheritz-Ponten T, Kissinger JC: A first glimpse into the pattern and scale of gene transfer in Apicomplexa. Int J Parasitol. 2004, 34: 265-274. 10.1016/j.ijpara.2003.11.025.PubMedView ArticleGoogle Scholar
- Jomaa H, Wiesner J, Sanderbrand S, Altincicek B, Weidemeyer C, Hintz M, Turbachova I, Eberl M, Zeidler J, Lichtenthaler HK, et al: Inhibitors of the nonmevalonate pathway of isoprenoid biosynthesis as antimalarial drugs. Science. 1999, 285: 1573-1576. 10.1126/science.285.5433.1573.PubMedView ArticleGoogle Scholar
- Puiu D, Enomoto S, Buck GA, Abrahamsen M, Kissinger JC: CryptoDB: The Cryptosporidium genome resource. Nucleic Acids Res. 2004, 32: D329-D331. 10.1093/nar/gkh050.PubMedPubMed CentralView ArticleGoogle Scholar
- Salzberg SL, Pertea M, Delcher AL, Gardner MJ, Tettelin H: Interpolated Markov models for eukaryotic gene finding. Genomics. 1999, 59: 24-31. 10.1006/geno.1999.5854.PubMedView ArticleGoogle Scholar
- Kissinger JC, Gajria B, Li L, Paulsen IT, Roos DS: ToxoDB: accessing the Toxoplasma gondii genome. Nucleic Acids Res. 2003, 31: 234-236. 10.1093/nar/gkg072.PubMedPubMed CentralView ArticleGoogle Scholar
- ToxoDB: the Toxoplasma genome resource. [http://toxodb.org]
- GenBank FTP site. [ftp://ftp.ncbi.nih.gov]
- The Sanger Institute: Entamoeba histolytica. [http://www.sanger.ac.uk/Projects/E_histolytica]
- GeneDB. [http://www.geneDB.org]
- TcruziDB. [http://TcruziDB.org]
- TIGR Tetrahymena thermophila genome project. [http://www.tigr.org/tdb/e2k1/ttg]
- Sanger Institute Pathogen Sequencing Unit. [http://www.sanger.ac.uk/Projects/Pathogens]
- Tetrahymena thermophila sequence source. [ftp://ftp.tigr.org/pub/data/Eukaryotic_Projects/t_thermophila]
- Swofford DL: PAUP* Phylogenetic Analysis Using Parsimony (*and other methods). 1998, Sunderland: Sinauer, 4.0b10Google Scholar
- Rujan T, Martin W: How many genes in Arabidopsis come from cyanobacteria? An estimate from 386 protein phylogenies. Trends Genet. 2001, 17: 113-120. 10.1016/S0168-9525(00)02209-5.PubMedView ArticleGoogle Scholar
- Jeanmougin F, Thompson JD, Gouy M, Higgins DG, Gibson TJ: Multiple sequence alignment with ClustalX. Trends Biochem Sci. 1998, 23: 403-405. 10.1016/S0968-0004(98)01285-7.PubMedView ArticleGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.PubMedView ArticleGoogle Scholar
- Felsenstein J: PHYLIP: Phylogenetic Inference Package. 2002, Seattle: Department of Genetics, University of Washington, 3.6aGoogle Scholar
- Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci . 1992, 8: 275-282.PubMedGoogle Scholar
- Treeview X. [http://darwin.zoology.gla.ac.uk/~rpage/treeviewx/]
- Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. 1989, New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- Schroeder AA, Brown AM, Abrahamsen MS: Identification and cloning of a developmentally regulated Cryptosporidium parvum gene by differential mRNA display PCR. Gene. 1998, 216: 327-334. 10.1016/S0378-1119(98)00360-6.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.