- Research article
- Open Access
Mitochondrial genes support a common origin of rodent malaria parasites and Plasmodium falciparum's relatives infecting great apes
BMC Evolutionary Biology volume 11, Article number: 70 (2011)
Plasmodium falciparum is responsible for the most acute form of human malaria. Most recent studies demonstrate that it belongs to a monophyletic lineage specialized in the infection of great ape hosts. Several other Plasmodium species cause human malaria. They all belong to another distinct lineage of parasites which infect a wider range of primate species. All known mammalian malaria parasites appear to be monophyletic. Their clade includes the two previous distinct lineages of parasites of primates and great apes, one lineage of rodent parasites, and presumably Hepatocystis species. Plasmodium falciparum and great ape parasites are commonly thought to be the sister-group of all other mammal-infecting malaria parasites. However, some studies supported contradictory origins and found parasites of great apes to be closer to those of rodents, or to those of other primates.
To distinguish between these mutually exclusive hypotheses on the origin of Plasmodium falciparum and its great ape infecting relatives, we performed a comprehensive phylogenetic analysis based on a data set of three mitochondrial genes from 33 to 84 malaria parasites. We showed that malarial mitochondrial genes have evolved slowly and are compositionally homogeneous. We estimated their phylogenetic relationships using Bayesian and maximum-likelihood methods. Inferred trees were checked for their robustness to the (i) site selection, (ii) assumptions of various probabilistic models, and (iii) taxon sampling. Our results robustly support a common ancestry of rodent parasites and Plasmodium falciparum's relatives infecting great apes.
Our results refute the most common view of the origin of great ape malaria parasites, and instead demonstrate the robustness of a less well-established phylogenetic hypothesis, under which Plasmodium falciparum and its relatives infecting great apes are closely related to rodent parasites. This study sheds light on the evolutionary history of Plasmodium falciparum, a major issue for human health.
Malaria is an overwhelming public health problem all over the world. It kills one to three million people annually and infects 200 to 500 million others . Human malaria is induced by infections caused by a range of eukaryotic protists belonging to the phylum Apicomplexa. These organisms possess an endosymbiont of red algal origin  derived into an apical organelle, the apicoplast. This organelle is specialized in host cell invasion . Within Apicomplexa, malaria parasites, also called Haemosporidia, are characterized by their infection of vertebrate hosts, haemoglobin digestion, and a complex life cycle involving dipteran vectors feeding from their vertebrate hosts' blood .
Haemosporidia include the genera Leucocytozoon (bird parasites), Haemoproteus and Parahaemoproteus (Sauria, i.e. bird and reptile parasites), Plasmodium (saurian and mammalian parasites) and Hepatocystis (mammalian parasites) [4–6]. These five genera have long been defined by their morphological differences (e.g. storage of products of haemoglobin degradation in the case of Haemoproteus, Parahaemoproteus, Plasmodium and Hepatocystis), variations in their life cycle (e.g. asexual replication stage in erythrocytes for Haemoproteus, Parahaemoproteus and Plasmodium) and host and vector specificity [4, 5, 7, 8]. However, studies which attempted to link these phenotypic and life history traits with molecular data concluded that the latter provide deeper insight into the evolutionary history of Haemosporidia, and allow the identification of cryptic species which cannot be distinguished using only microscopy observations [9, 10]. The increasing availability of molecular data has enabled numerous studies of Haemosporidia phylogeny, improving our understanding of the evolutionary history of malaria parasites.
To date, five Plasmodium species have been shown to cause human malaria: P. falciparum, P. vivax, P. malariae, P. ovale and P. knowlesi [11–14]. P. falciparum has stimulated the interest of the scientific community, largely because it is the most virulent, but also because of the very high A+T contents of its genome. Indeed, it has an average content of 80% A+T over its nuclear genome, with intergenic regions frequently displaying more than 90% A+T . Furthermore, the question of the origin of P. falciparum is intriguing: its 18 S rRNA gene diverges strikingly from that of the other Plasmodium species which infect primates, of which four cause human malaria [16–19]. Another parasite of interest is P. reichenowi, which infects chimpanzees. P. reichenowi has long been the only known close relative of P. falciparum [6, 18–24]. However, since 2009, it has been recognized that both P. falciparum and P. reichenowi belong to a wider monophyletic lineage of parasites specialized in the infection of great ape hosts: gorilla, chimpanzee, bonobo and human [25–28].
Despite intense efforts to obtain accurate molecular phylogenies, the position of P. falciparum and of great ape malaria parasites within Haemosporidia remains unclear. In recent decades, many molecular phylogenies have been produced which clustered P. falciparum (alone, or with P. reichenowi) closer to avian than to mammalian parasites (Figure 1A). This led to the hypothesis of a recent switch from avian to human (and great ape) hosts, rather than a shared ancestry of P. falciparum and other mammalian malaria parasites [16–19, 21, 29–35]. However, it has recently been recognized that, in these early works, the clustering of P. falciparum with avian malaria parasites may have resulted mainly from stochastic noise due to the insufficient number of sequences analyzed . More recent studies relied on larger data sets and demonstrated that all malaria parasites known to exclusively infect mammals, including P. falciparum, are monophyletic [6, 20–28, 36]. Hence, this monophyletic clade of mammalian malaria parasites includes (i) a lineage infecting great apes [25–28], (ii) a distantly related and less specialized lineage infecting primates , (iii) a third lineage of rodent parasites [6, 37], and (iv) species of the Hepatocystis genus infecting bats and primates [5, 6, 21]. Most studies inferred great ape parasites to be a sister-group to all other mammal parasites [6, 20–22, 24, 25, 27, 28, 36] (Figure 1B). However, a few studies yielded contradictory results. Some suggested that great ape parasites could be closer to other primate parasites  (Figure 1C), while others suggested that great ape parasites could be related to rodent parasites [23, 27] (Figure 1D).
In the present study, we distinguished between these mutually exclusive hypotheses to identify the origin of Plasmodium falciparum and other great ape parasites. We analyzed three genes from 33 publicly available complete Haemosporidian mitochondrial genomes, as well as 51 additional sets of one to three mitochondrial genes from Hepatocystis and Plasmodium species infecting various mammalian and saurian hosts. Concatenated nucleotide and amino-acid alignments were analyzed using various probabilistic models of sequence evolution, applying both maximum-likelihood (ML) and Bayesian inference (BI) methods. Statistical measurements of fit and posterior predictive experiments enabled the adequacy of model assumptions to the data to be checked. We also evaluated the robustness of the obtained phylogeny to the taxonomic sampling by removing and adding taxa. Our phylogenetic analyses provide robust support for a close relationship between great ape and rodent parasites (Figure 1D).
Results and Discussion
In this study, we analyzed coding genes from the 33 complete Haemosporidian mitochondrial genomes available in 2009 (see Methods and Additional file 1, Table S1). Each genome includes only three coding genes. In this section, we describe the general features of the resulting alignments.
In the three concatenated genes and 33 taxa data set, 54% of the 1099 amino-acid sites, and 51% of the 3308 nucleotide sites were constant. The nucleotide and amino-acid alignments are thus highly conserved. Genes displayed the typical pattern in which substitution rates are higher at the third and first codon positions, and lower at the second codon positions. The first and second codon positions are highly conserved (55% and 77% of constant positions, respectively) and the saturation plot indicates that they are slightly saturated (slopes of 0.19 and 0.89 respectively, Additional file 2, Figure S1). Although third codon positions evolved more rapidly and are thus more saturated (20% of constant sites, slope of 0.09, Additional file 2, Figure S1), another standard test for saturation (PAUP 4.0 , partition-homogeneity test by codon position) indicated that they are not significantly more saturated than the first and second codon positions (p = 1).
The nucleotide data set had a 74% A+T content and was compositionally homogeneous (p = 0.83, PAUP χ2 test of compositional homogeneity across taxa). However, the amount of A+T and the compositional homogeneity differed strikingly depending on the codon position. First and second codon positions displayed relatively low A+T contents (68% and 64% respectively) and were homogeneous (p = 1, PAUP χ2 test of compositional homogeneity across taxa). In contrast, third codon positions had a high A+T content (89%) and were compositionally heterogeneous (p = 0, PAUP χ2 test of compositional homogeneity across taxa). This difference in composition between codon positions suggests that the typical drift of Haemosporidian species toward A+T richness was negatively selected at first and second codon positions and less constrained at third positions, most likely due to constraints at the protein level.
Standard phylogenetic analyses
Phylogenetic reconstructions in this section were obtained with the most widely used models, which have both maximum-likelihood (ML) and Bayesian implementations. MrAIC  estimated that, among 56 models of nucleotide substitution, the best AIC score  was achieved by the most general homogeneous and reversible model (GTR, "general time reversible" ), combined with four discrete gamma categories of "rates across sites" , plus an additional rate category for invariant sites (model GTR nt + Γ4 + I, where subscript nt stands for nucleotides). This model was also selected by MrAIC with individual codon positions. We used the PhyloBayes 3.0 software  for Bayesian analyses, which does not implement the invariant rate category. Bayesian analyses of nucleotide alignments were thus performed under the GTR nt + Γ4 model.
ML and Bayesian analyses of the nucleotide data set under GTR nt + Γ4 + I and GTR nt + Γ4 models, respectively, strongly supported a clade containing the 20 mammal infecting Plasmodium species (Figure 2) [6, 20–25, 27, 28, 36]. The mammal malaria parasite clade comprises three strongly supported main lineages. One lineage is specialized in infecting great ape hosts and includes P. falciparum, P. reichenowi and P. gaboni [25–28]. The second lineage is characterized by African and Asian primate hosts and comprises 14 Plasmodium species . The third lineage includes P. berghei, P. yoelii and P. chabaudi, the African rodent malaria parasites . Most interestingly, the two lineages of great ape and rodent parasites clustered together with strong posterior probability (PP = 0.99, Figure 2). ML analysis also yielded significant support for this clade (Shimodaira-Hasegawa-like support SH = 0.93, bootstrap support BS = 0.82, Table 1).
First and second codon positions, and third codon positions were analyzed separately, using the same models of nucleotide evolution as previously. Both data sets supported, although weakly, a clade containing great ape and rodent malaria parasites, hereafter denoted "monophyly of great ape and rodent parasites" (Additional file 3, Table S2, column "Rodent").
Finally, ProtTest  determined that the best fit to the amino-acid alignment was provided by the JTT + Γ4 + I + F model . We used the JTT + Γ4 model for the Bayesian analysis of the amino-acid alignment. The monophyly of great ape and rodent parasites received weak SH and bootstrap support (SH = 0.46, BS = 0.41), but relatively high posterior probability (PP = 0.92, Table 1).
SH and BS supports are more conservative than posterior probabilities, which are generally expected to be higher [46, 47]. Moreover, given the high level of conservation of protein and individual codon position alignments (see above), the complete nucleotide alignment is expected to provide more phylogenetic signal and higher supports. Thus, all previous results are congruent, with differing but explainable levels of confidence, and support a monophyly of great ape and rodent parasites.
Assessment of model violation and robustness to the model choice
To evaluate the influence of potential model violations, we used improved Bayesian models implemented in PhyloBayes 3.0 software . The fits of these additional models to the data were measured using cross-validated likelihood (see Methods and Additional file 4, Table S3). Moreover, we applied posterior predictive tests which measure the model ability to accurately reproduce observed features of the data (see Methods and Additional file 5, Table S4). In this section, we only comment on experiments where posterior predictive tests were not rejected (i.e. there was significant violation of model assumptions). In this case, we conclude that a specific feature of the data is correctly anticipated by the model assumptions.
We wanted to evaluate the potential effects of site saturation on our estimates. Site-heterogeneous mixture models such as "CAT"  (see Methods) efficiently deal with violations caused by high saturation levels [48, 49]. However, this model might lack resolution power, especially in the case of small data sets . In addition, we wanted to measure the potential effect of compositional biases on our estimates. Consequently, we analyzed data sets under the time-heterogeneous model "BP" [50, 51], which is designed to deal with compositional heterogeneity across taxa (see Methods).
Saturation in the nucleotide alignment
Among the three additional models considered, CAT + GTR nt + Γ4  (see Methods) yielded the best fit to the nucleotide data set, and it outperformed GTR nt + Γ4 by 103 points of cross-validated likelihood (Additional file 4, Table S3). Moreover, posterior predictive tests showed that this model correctly anticipated the level of saturation of the nucleotide data set (p > 0.28, Additional file 5, Table S4). This suggests that the CAT + GTR nt + Γ4 model is not misled by site saturation. CAT + GTR nt + Γ4 strongly supported the monophyly of great ape and rodent parasites, considering either all codon positions (PP = 0.99, Table 1), first and second (see Additional file 6, Figure S2), or third codon positions (PP = 0.92 and PP = 0.65 respectively, Additional file 3, Table S2).
Compositional heterogeneity in the nucleotide alignment
Third codon positions were compositionally heterogeneous (see above χ2 tests), and they carried 54% of variable sites. Hence, potential convergence of sequence compositions could have misled the previously used time-homogeneous models [50–52]. Interestingly, posterior predictive tests showed that the compositional heterogeneity across taxa of first and second codon positions was correctly anticipated under the time-homogeneous models GTR nt + Γ4 and CAT + GTR nt + Γ4 (p > 0.44, Additional file 5, Table S4), suggesting that these models are relatively robust to compositional changes in this data set.
The time-heterogeneous model GTR nt + BP + Γ4  (see Methods) explicitly accounts for variations in composition across taxa. It correctly anticipated the observed compositional heterogeneity, considering either all, first and second, and third codon positions (p > 0.11, Additional file 5, Table S4), suggesting that this model is unlikely to be misled by compositional heterogeneity across taxa. The monophyly of great ape and rodent parasites was strongly supported under the GTR nt + BP + Γ4 model, considering either all codon positions (PP = 0.95, Table 1), first and second, or third codon positions (PP = 0.99 and PP = 0.72 respectively, Additional file 3, Table S2).
Model violations in the amino-acid alignment
According to the analysis of the amino-acid data set, methodological bias may arise from the use of the universal replacement model JTT. The peculiar A+T rich composition of Haemosporidian genes could lead to slightly different estimations for the exchange rate parameters and hence alter the probability of clustering great ape and rodent parasites together. Accordingly, we used the GTR aa model (where subscript aa stands for amino acids) which does not rely on a pre-estimated replacement matrix like JTT. The GTR aa + Γ4 model strongly supported the monophyly of great ape and rodent parasites (PP = 0.98, Table 1).
Cross-validations indicated that, among the 12 alternative models, the site-heterogeneous model CAT + JTT + Γ4  (see Methods) provided the best fit to the amino-acid data set (Additional file 4, Table S3). According to posterior predictive tests, this model correctly anticipated the saturation level observed in the data (p > 0.07, Additional file 5, Table S4). Moreover, the site- and time- heterogeneous model CAT + BP + Γ4  (see Methods) correctly anticipated the level of saturation of the amino-acid data set (p = 0.41). However, the posterior predictive test for compositional heterogeneity across taxa was rejected under CAT + BP + Γ4 (p = 0.02), although as expected, this model anticipated compositional heterogeneity better than do time-homogeneous models (p = 0.001, Additional file 5, Table S4). Both the two last models moderately supported the monophyly of rodent and great ape parasites (PP = 0.80 and PP = 0.75 under CAT + JTT + Γ4 and CAT + BP + Γ4 respectively, Table 1). In both cases, few variable amino-acid positions are interpreted under highly parameter rich models, and moderate support is therefore to be expected.
Most importantly, all results in this section are congruent with our initial estimate. This suggests that measured model violations do not significantly bias the relationship of great ape parasites among mammal parasites. In other words, the monophyly of great ape and rodent parasites appears to be robust to the choice of the model, as well as to its assumptions and dimensionality.
Robustness to taxon sampling
To avoid possible biases resulting from an ad hoc set of sequences, it is important to assess the robustness of this phylogenetic association with respect to taxon selection, which was achieved by analyses of 30 different taxonomic samples. Among previously considered data sets, the complete nucleotide data set had the highest number of variable positions to interpret (1627), and our experiments showed that it is unlikely to induce strong model violation, whatever the phylogenetic model considered. Hence, in the following sections, we present phylogenetic analyses of complete nucleotide alignments performed under models GTR nt + Γ4 + I (ML) and GTR nt + Γ4 (Bayesian inference).
Robustness to taxon removal
First, we checked the influence on phylogenetic reconstructions of the selected great ape parasites. According to the initial taxon selection, this lineage includes P. falciparum, P. reichenowi and P. gaboni. Six additional combinations of these three taxa were devised (see Methods). The weakest (but still relatively high) support for the monophyly of great ape and rodent parasites was obtained when P. reichenowi was considered as the only representative of its lineage (PP = 0.93, SH = 0.78, BS = 0.68, Additional file 7, Table S5). Second, we devised six combinations of the three rodent parasites, P. berghei, P. yoelii and P. chabaudi. The data set with P. berghei as the only representative of its lineage yielded the weakest support for the monophyly of great ape and rodent parasites (PP = 0.56, BS = 0.52, Additional file 7, Table S5), and the ML tree weakly supported the alternative hypothesis of a monophyly of primate and rodent parasites (i.e. Figure 1B, SH = 0.05). Third, six combinations of primate parasites were considered. Only the combination with African primate parasites (P. gonderi and P. DAJ-2004) as the only representatives of their lineage supported the alternative hypothesis of a monophyly of primate and rodent parasites (i.e. Figure 1B, PP = 0.85, SH = 0.26, BS = 0.48, Additional file 7, Table S5). The five other combinations of primate parasites supported the monophyly of great ape and rodent parasites (PP > 0.96, SH > 0.60, BS > 0.71, Additional file 7, Table S5). Finally, we investigated the robustness to the taxon composition of the outgroup (i.e. mammal and saurian parasites were considered as ingroup and outgroup, respectively). All six devised outgroups yielded high support for the monophyly of great ape and rodent parasites (PP > 0.98, SH > 0.64, BS > 0.74, Additional file 7, Table S5).
Hence, with the exception of three taxonomic samples, in which (i) P. reichenowi, (ii) P. berghei and (iii) P. gonderi and P. DAJ-2004 were considered as the only representatives of their respective lineages, all other 21 combinations of taxa provided good support for the association of great ape and rodent parasites (PP > 0.96, SH > 0.60, BS > 0.71, Additional file 7, Table S5).
Robustness to taxon addition
First, (i) 8 CytB genes from great ape parasites, (ii) 10 pairs of CytB and Cox1 genes from rodent parasites and, (iii) 27 pairs of CytB and Cox1 genes from Plasmodium species infecting a wide range of sauria hosts (Additional file 8, Table S6) were added in turn to the initial 33-taxon data set. The association of great ape and rodent parasites was still strongly supported (PP > 0.99, SH > 0.77, BS > 0.78, Additional file 7, Table S5), and the lineages of great ape, rodent, and mammal parasites were each still shown to be monophyletic (PP = 1).
Second, all previous taxa were analyzed together (33 taxa + 8 great ape + 10 rodent + 27 saurian parasites), yielding a 78-taxon tree. The monophyly of great ape and rodent parasites was still strongly supported (PP = 0.99, SH = 0.86, BS = 0.76, Additional file 7, Table S5).
Third, we added six Hepatocystis species to the initial 33-taxon nucleotide data set. These six parasites were monophyletic (PP = 1). They clustered within the clade of mammal parasites, which was then composed of four monophyletic main lineages. The monophyly of great ape and rodent parasites was weakly supported (PP = 0.40, SH = 0.05, BS = 0.38), but this low support was entirely due to high uncertainty with respect to the position of Hepatocystis within mammalian malaria parasites. Indeed, Hepatocystis were located, with weak support, in five positions on trees in which great ape and rodent parasites were located close together (e.g. Hepatocystis as a sister-group to great ape parasites: PP = 0.28, BS = 0.28, or to rodent parasites: PP = 0.31, BS = 0.29). We evaluated posterior and bootstrap support for the great ape being close to, but not necessarily monophyletic with, rodent parasites (i.e. the great ape plus rodent parasite lineage could also include Hepatocystis). When support was summed over the three possible positions of Hepatocystis relative to the association of great ape and rodent parasites, then parasites of great apes and of rodents were located close together with strong support (PP = 0.99 and BS = 0.87, Additional file 7, Table S5).
Fourth, the six Hepatocystis species were analyzed simultaneously with all the previous 78 taxa, yielding an 84-taxon tree (Figure 3). A strict monophyly of great ape and rodent parasites was weakly supported (PP = 0.28, SH = 0, BS = 0.31), due to high uncertainty with respect to the position of Hepatocystis (e.g. Hepatocystis as a sister-group to great ape parasites: PP = 0.36, BS = 0.25, or to rodent parasites: PP = 0.35, BS = 0.25). However, disregarding the exact position of Hepatocystis, great ape and rodent parasites were located close together with high support (PP = 0.99, BS = 0.75, Additional file 7, Table S5).
Thus, all six previous additions of mitochondrial genes did not alter the result indicating a likely close phylogenetic relationship of great ape and rodent parasites. Overall, statistical supports averaged over the 30 taxonomic samples considered in this section showed a robust relationship of great ape and rodent parasites . This suggests that this relationship does not depend solely on the selection of the taxa we considered here. However, the uncertainty concerning the exact position of Hepatocystis species within mammal parasites challenges the above mentioned monophyly of a clade only comprising parasites of great apes and of rodents, because it would be possible for Hepatocystis to cluster within that clade. Nonetheless, whatever the true position of Hepatocystis may be, it does not contradict our main result indicating that great ape parasites are unlikely to be a sister-group to all other mammal parasites [6, 20–22, 24, 25, 27, 28, 36] (Figure 1B), but instead, probably share a more recent common ancestor with rodent parasites [23, 27] (Figure 1D).
Comparison with previous studies
In previous studies, the hypothesis of great ape parasites being a sister-group to all other mammal parasites was defended by analyses of mitochondrial genes [6, 21, 22, 24, 27, 35] or complete mitochondrial genomes [24, 25, 28], of nuclear coding genes  and ribosomal RNA , or by combining genes from nuclear and mitochondrial genomes with genes from the apicoplast genome . Moreover, considering different rooting assumptions breaking the monophyly of mammal parasites, this result was also obtained through analyses of nuclear 18 S rRNA [11, 18], nuclear genes [29, 32, 34] and mitochondrial cytochrome b genes .
Most of these studies of the Haemosporidia phylogeny relied on a single gene, and only a few taxa data sets, which might lack phylogenetic signal . In contrast, two recent studies analyzed larger data sets. The first study analyzed a large number of taxa (40), but few concatenated genes (4) , whereas the second focused on a large number of genes (104), but considered very few taxa (8) . Next, we suggest possible reasons for the disagreement between the results of these two studies and ours.
Comparison with a taxa-wide phylogenetic analysis
As a general guideline, wider taxon sampling usually helps to resolve phylogenies more accurately, provided enough genes are available to overcome stochastic noise, and are also sufficiently conserved to avoid systematic errors . In line with this idea, Martinsen et al.  analyzed four concatenated genes for a relatively wide sample of 40 taxa. Among the previous works, their experimental conditions are thus the closest to ours. But, intriguingly, our results do not confirm theirs. We suggest that the disagreement between the two studies is due to several factors, the first being the differences in the phylogenetic markers analyzed.
Both Martinsen et al.  and our study considered CytB and Cox1 mitochondrial genes. However, Martinsen et al.  additionally analyzed adenylosuccinate lyase (ASL) and caseinolytic protease C (ClpC) genes, whereas we analyzed the third mitochondrial gene Cox3. In order to compare global rates of evolution between these genes, we measured the total lengths of gene trees , for a common sub set of eight taxa (P. falciparum, P. reichenowi, P. vivax, P. knowlesi, P. berghei, P. chabaudi, P yoelii and P. gallinaceum). Our values indicate that the ASL genes evolved 3 to 5 times faster than the slowest evolving genes: ClpC, Cox1 and Cox3 (Figure 4). The signal to noise ratio is expected to be higher for slowly evolving phylogenetic markers, and fast rates of evolution generally reduce the accuracy of inferred phylogenetic trees [49, 56]. In addition, we considered ASL and ClpC genes of 18 and 27 taxa respectively, for taxonomic samples as close as possible to our originally selected 33 taxa (Additional file 9, Table S7). The rapidly evolving gene, ASL, did not support a monophyly of mammal parasites, suggesting strong systematic errors (Additional file 10, Figure S3). In contrast, the slow evolving gene, ClpC, supported this monophyly (PP = 0.99, SH = 0.94, BS = 0.51), but did not support any particular position of P. falciparum within mammal parasites (Additional file 11, Figure S4). However, a recent study of 14 ClpC genes supported a common ancestry of great ape and rodent parasites .
Moreover, most CytB, ClpC and ASL sequences analyzed by Martinsen et al.  were partial CDS. Thus, even considering four genes, their alignment covered 2334 nucleotide sites, representing 70% of the 3308 nucleotide sites considered in the present study. Finally, they considered fewer mammalian malaria parasites than we did (11 taxa in the study of Martinsen et al. , versus 20 to 44 taxa in our work), and they considered P. falciparum as the only representative of the lineage of great ape parasites, together with P. vivax and P. knowlesi as the only representatives of primate parasites. All these differences in experimental conditions (i.e. saturation of ASL gene, fewer sites, fewer mammal parasites close to P. falciparum) could together contribute to the difference between the results of Martinsen et al.  and ours.
Comparison with a genome-wide phylogenetic analysis
Dávalos and Perkins  extracted a set of 104 putative orthologous nuclear genes from the eight complete genomes of Plasmodium species sequenced to date (P. falciparum, P. reichenowi, P. vivax, P. knowlesi, P. berghei, P. chabaudi, P yoelii and P. gallinaceum). Their phylogenetic analyses of individual genes displayed discrepancies with respect to the inferred trees. Approximately half the 104 genes supported a monophyly of primate and rodent parasites (Figure 1B) [6, 20–22, 24, 25, 27, 28, 36]. Alternative hypotheses of a monophyly of great ape and primate parasites (Figure 1C) , and of great ape and rodent parasites (Figure 1D) [23, 27], were each supported by nearly a quarter of these 104 genes. Interestingly, for the same reduced sample of eight taxa, our phylogenetic analyses of mitochondrial genes, taken separately, displayed comparable discrepancies (Table 2). Moreover, the supports obtained showed a near-random resolution of the internal branch of the 8-taxon individual gene trees. This suggests that, rather than showing a global preference of individual genes for the monophyly of primate and rodent parasites, the analysis of individual genes by Dávalos and Perkins  might have been strongly influenced by stochastic noise.
Increasing the amount of signal by concatenating genes helps to alleviate the effects of stochastic noise. The analysis of the three concatenated mitochondrial genes, for the eight taxa, supported the monophyly of great ape and rodent parasites (PP = 0.97, SH = 0.58, BS = 0.65; Table 2). In contrast, the concatenation of the 104 nuclear genes yielded strong support for the monophyly of primate and rodent parasites  (Figure 1B). We estimated the total tree lengths of each of these 104 genes, as well as of their concatenation. First, mitochondrial genes evolved as slowly as the fraction of the 104 nuclear genes which displayed the slowest rate of evolution (Figure 4). Second, total tree length estimated for the 104 gene concatenation indicated a fast average rate of evolution, about two times faster than that of the three mitochondrial genes. These observations confirm the fact that mitochondrial genes are well conserved, and corroborate the conclusion of Dávalos and Perkins  indicating that most of their 104 genes are highly saturated and evolved relatively fast. Thus, the monophyly of primate and rodent parasites (Figure 1B), obtained by the latter authors from a large concatenation of 104 genes, most likely results from systematic errors that may be due to the high saturation level of most genes [53, 58], but also presumably to the small sample of only eight taxa .
Comparison with two corroborated studies
In contrast, our results support a common origin of great ape and rodent parasites (Figure 1D). This corroborates results of a recent study published by Perkins  who, to the best of our knowledge, was the first to mention this hypothesis. This author sequenced seven new mitochondrial genomes and reconstructed the phylogeny of a sample of 24 taxa. A similar phylogeny was also obtained from 38 partial CytB sequences . However, the two previous studies did not discuss the robustness of this result, but instead suggested it should be considered with caution [23, 27], as most previous studies of mitochondrial genes supported a monophyly of primate and rodent parasites (Figure 1B) [6, 21, 22, 24, 24, 25, 27, 28, 35].
However, in our results, this hypothesis never obtained significant statistical support: at most, it reached PP = 0.35 and BS = 0.57 with the 33-taxon amino-acid alignment (Additional file 3, Table S2, column "Primate+Rodent"), and it obtained averaged supports of and over the 30 additional taxon samples considered (Additional file 7, Table S5, column "Primate+Rodent").
Interestingly, each time the monophyly of great ape and rodent parasites (Figure 1D) was not significantly supported, the only alternative hypothesis which could not be statistically rejected was the monophyly of primate and rodent parasites (i.e. support for "Primate+Rodent" in Additional files 3 and 7, Tables S2 and S5, where PP > 0.05). Moreover, nine additional samples of 19 taxa were drawn so that trees would display long branches, and we obtained five trees in which great ape parasites were shown to be a sister-group to primate and rodent parasites (Additional file 12, Table S8). This suggests that the monophyly of primate and rodent parasites (Figure 1B), rather than the monophyly of great ape and rodent parasites (Figure 1D), might result from the effect of long branch attraction. Thus, the slight tendency of mitochondrial genes to weakly support a monophyly of primate and rodent parasites, along with differences in taxon sampling and gene selection, could explain the disagreements between our results and most previous studies.
With special focus on the still unclear phylogenetic position of great ape parasites, which include Plasmodium falciparum [25–28], in this study, we explored the phylogeny of Haemosporidian species by analyzing their mitochondrial genes. We showed that these genes have evolved relatively slowly and are mostly compositionally homogeneous, which characterizes them as potentially accurate phylogenetic markers. Corroborating many results obtained over the past few years, we obtained a monophyly of mammalian malaria parasites. Within that clade, we observed four main host-specialized lineages of parasites: Plasmodium species infecting (i) primate, (ii) rodent and (iii) great ape hosts, and (iv) Hepatocystis species infecting bats and primates. The inferred relationships within host-specialized lineages of Plasmodium parasites are congruent with the literature [11, 26, 37]. Hepatocystis species have received little attention to date. According to our results, these parasites may have diverged within mammalian malaria parasites, but their exact origin remains unclear. Our results support a common ancestry of great ape and rodent parasites (Figure 1D). We showed that this phylogenetic relationship is robust to various experimental conditions, demonstrating that it is unlikely to arise from an artefact of tree reconstruction.
Our study focused on mitochondrial genes. Nevertheless, it is still unclear if mitochondrial genomes match the Haemosporidian species tree or not . This question could be answered by comparing the respective phylogenies of the three genomes hosted by Haemosporidian species (i.e. phylogenies of the mitochondrial, nuclear and apicoplast genomes). However, appropriate data sets for phylogenetic reconstruction of nuclear and apicoplast genomes are not yet available. A successful strategy for resolving a gene- and taxa- wide phylogeny would involve targeted sequencing of identified slowly evolving genes from the apicoplast and nuclear genomes. Careful phylogenetic analyses of such new and accurate phylogenetic markers will likely help to definitely resolve the phylogenetic origins of Plasmodium falciparum and other great ape parasites.
Mitochondrial gene and protein alignments
Complete mitochondrial genomes of 33 Haemosporidian parasites [22, 23, 25, 59–66] were downloaded from the NCBI website (Additional file 1, Table S1). Available taxa included four Leucocytozoon species infecting birds [23, 66]. These parasites are the closest relatives of other already identified Haemosporidia , and were used as outgroups. We collected mitochondrial genomes of four Haemoproteus and Parahaemoproteus species [23, 62], and of five Plasmodium species [22, 23, 62, 65] infecting saurian hosts (birds and reptiles). We collected a set of 20 mitochondrial genomes of mammal malaria parasites, including those of three rodent , three great ape [25, 60, 61] and 14 primate parasites [22, 63, 64].
Mitochondrial genomes of Haemosporidian species are vestigial and have a typical length of 6, 000 base pairs [67, 68]. They form linear concatemers, each repeated unit encoding fragments of ribosomal genes together with three coding genes involved in the electron transport chain: cytochrome b and cytochrome oxidase subunits I and III (denoted as CytB, Cox1 and Cox3, respectively). Both nucleotide coding sequences and their translations into amino acids were retrieved according to annotated CDS. No annotation of the mitochondrial genes was available for Plasmodium relictum jb5, Parahaemoproteus jb1.JA27 and Parahaemoproteus jb2.SEW5141. In this case, gene sequences were extracted manually and unambiguously, given the high conservation level of start and end positions shared by all 30 other annotated genes.
Nucleotide sequences were aligned using MACSE (Ranwez V, Harispe S, Delsuc F, and Douzery EJP, "MACSE: Multiple Alignment of Coding Sequences accounting for frameshifts and stop codons", manuscript in preparation). This method computes the alignment of coding nucleotide sequences with respect to their possible translations. It attempts to minimize the occurrence of frameshifts and stop codons. We applied the relevant codon table, as indicated on the NCBI website . MACSE identified and corrected three long frameshifts in the Cox3 gene of P. berghei, and in CytB genes of Ha. columbae and L. caulleryi. These frameshifts resulted in erroneous translations of the corresponding publically available amino-acid sequences. Consequently, we used the translated alignment computed by MACSE, rather than the alignment of the official NCBI translations. Individual gene and protein alignments were filtered with Gblocks 0.91 , and allowing a maximum of half gap states per site (option -b5 = h). Filtered alignments of nucleotide and amino-acid sequences were finally concatenated. This yielded two concatenations, one of 3308 nucleotide sites (number of nucleotide sites for CytB: 1125, Cox1: 1434, Cox3: 749), and one of 1099 amino-acid positions. The assembled data sets are available at the following URL: http://www.lirmm.fr/mab/blanquart
Description of phylogenetic models
Data sets were analyzed under various probabilistic models of molecular evolution. We applied the JTT  replacement model to the amino-acid data set. We applied the most general time reversible model GTR nt , where subscript nt denotes nucleotides, to the nucleotide data sets . These substitution models were run using both maximum-likelihood (ML) and Bayesian inference (BI) methods. This allowed for the use of different statistical supports with different meanings, and comparison of the phylogenies estimated with the two approaches.
ML phylogenetic reconstructions were performed using PhyML 3.0 . Irrespective of the substitution model (GTR nt , JTT) used for the analysis, the phylogenetic model additionally involved four discrete categories of gamma distributed rates across sites (denoted + Γ4, ), plus an invariant site category (denoted +I). The proportion of invariant sites and the shape parameter of the gamma distribution were estimated from the data. When the nucleotide data sets were analyzed, all eight free parameters of the GTR nt substitution models were estimated from the data (GTR nt + Γ4 + I, 10 degrees of freedom). In the case of the amino-acid alignment analyzed under JTT, stationary probabilities were set to empirical frequencies of amino acids measured over the whole data set (JTT + Γ4 + I + F, 21 degrees of freedom). Note that these models were identified as the available ML models that best fit the sequence alignments, according to the AIC criterion .
Bayesian phylogenetic reconstructions were performed using PhyloBayes 3.0 . For all Bayesian experiments performed in this study, two independent MCMC chains - each starting from a random point - were run for up to 100, 000 cycles. One MCMC sample was saved every 10 cycles, and the first 500 samples were discarded as "burnin". The eight free parameters of GTR nt and the amino-acid frequencies of JTT (19 free parameters), were estimated from the data.
We also applied more general and parameter rich models of evolution, implemented in a Bayesian framework. We used GTR aa , where subscript aa indicates a GTR model dedicated to amino-acid sequences. This model directly estimates the exchange rate parameters from the data (208 degrees of freedom). Models JTT and GTR aa homogeneously apply a single substitution model to the whole data set. However, in some cases, this parameterization is prone to violations by the data, resulting in wrong phylogenetic inferences . Consequently, we applied the site-heterogeneous mixture model CAT to the amino-acid alignment, which implements a mixture of stationary probability vectors across sites . The CAT model was combined with free (+GTR nt ) or empirical (+JTT) relative exchange rates, applied to the nucleotide and amino-acid alignments, respectively. Both CAT + GTR nt and CAT + JTT models were combined with discretized gamma rates across sites (+ Γ4). Finally, we analyzed both nucleotide and amino-acid alignments under time-heterogeneous models of evolution. The BP model component allows for changes over time of the substitution model stationary probabilities and hence, estimates the compositional drift of the sequences . We applied the GTR nt + BP + Γ4  and the CAT + BP + Γ4  models to the nucleotide and the amino-acid alignments, respectively.
Additional models with Bayesian implementations were compared using cross-validation (see below). In addition to the + Γ4 model of rate variation across sites, we applied a covarion model (+COV) which enabled us to estimate site specific rate variations (i.e. heterotachy) . In addition to JTT and GTR aa , we considered the MtREV  empirical rate matrix. Finally, in addition to the mixture model CAT, we considered the empirical mixture models UL2 and UL3 . These components allowed 13 and 3 models of evolution to be derived and applied to amino-acid and nucleotide alignments, respectively.
The fit of the models implemented in a Bayesian framework was estimated by cross-validation, as implemented in PhyloBayes 3.0 . Ten replicate data sets were randomly drawn. The learning part of each replicate data set comprised 90% of the sites of the whole alignment. The 10% of remaining sites were used to compute the cross-validated likelihood. The tree topology was considered as a free parameter. Note that, for computational reasons, fits of time-heterogeneous models GTR nt + BP + Γ4 and CAT + BP + Γ4 were not evaluated.
Posterior predictive experiments
Model violations were measured by posterior predictive experiments, as implemented in PhyloBayes 3.0 . We applied a test statistic measuring the compositional heterogeneity across taxa. The test statistic "composition" was defined as the maximum of the χ2 distances separating each sequence composition from the composition of the whole data set . We applied two test statistics to measure site saturation. The test statistic "site diversity" measures the mean state diversity across sites  (e.g. a constant site has a diversity of 1). The test statistic "homoplasy" considers the averaged number of convergence and reversion events per site, as displayed by inferred stochastic mapping . A posterior predictive test compares the value V O of a test statistic measured given the observed data, to the distribution of that test statistic measured over simulated replicate data sets. Each replicate data set was simulated given an a posteriori drawn sample of parameters. The p-value indicates the probability of observing a test statistic as extreme as V O , under the null hypothesis stating that the model assumptions are true. Failure to reject a posterior predictive test indicates that the model assumptions allow to realistically reproduce the observation V O based on real data.
Saturation of the phylogenetic signal of each codon position was illustrated by a saturation plot . For each pair of taxa in an alignment, we plotted their "pairwise similarity distance" (i.e. y-coordinates: number of sites displaying different states, normalized by the alignment length), versus the distance separating these two taxa along the tree branches (i.e. x-coordinates: the sum of branch lengths from the two taxa to their common ancestor). We used a fixed tree topology estimated from all codon positions (Figure 2). Branch lengths and other model parameters were evaluated according to separate codon positions.
Estimation of tree node supports
Under ML analysis, statistical support of tree branches was estimated from 1000 bootstrap replicates and, in addition, using the Shimodaira-Hasegawa-like test (SH) implemented in PhyML 3.0 [47, 75]. Bayesian analysis classically provides a collection of samples drawn from the a posteriori distribution. The posterior probability of observing a given phylogenetic association between two lineages is then approximated by its frequency among sampled trees. Given a monophyletic target lineage A (e.g. great ape parasites), we extracted from a tree collection the list of all its N different sister-groups B n (e.g. rodent parasites). We then computed the frequency : the posterior support of clade A + B n (e.g. great ape plus rodent parasites). The same approach was used for the bootstrap support, but not for the SH support which applies only to clades that belong to the ML tree.
Let A and B each be a monophyletic lineage within a phylogenetic tree, and let them form a well-supported monophyletic clade A + B according to an initial sample of taxa. To ensure that the target relationship between A and B does not result from stochastic noise (lack of signal) or systematic error (model violation), we checked its robustness to taxon sampling (e.g. ). Given a lineage C (possibly equal to A or B) including k taxa, we checked that every combination of 1 to k - 1 taxa of C yielded a congruent phylogeny with respect to the phylogenetic relationship A + B. If lineage C was composed of too many taxa, we selected only a few relevant taxon combinations among all those available.
We focused on the robustness of the association of lineage A: great ape parasites, with lineage B: rodent parasites. Each of these two lineages was considered in turn as a sampled lineage C. According to our initial selection of 33 taxa, both these lineages were composed of 3 taxa and, 6 combinations of single or pairs of representatives were considered. We also considered in turn the group of primate parasites (14 taxa) and the saurian parasite outgroup (13 taxa) as sampled lineages C. For each of these sampled lineages, only 6 combinations of single or pairs of sub-groups were considered. For the 14 primate parasites, the 3 sub-groups were: (a) P. malariae and P. ovale (infecting humans), (b) P. gonderi and P. DAJ-2004 (African primate parasites), and (c) 10 Plasmodium species infecting Asian primates (Additional file 1, Table S1). For the 13 saurian parasites, the 3 sub-groups were (a) the 5 Plasmodium species, (b) the 4 Haemoproteus and Parahaemoproteus species, and (c) the 4 Leucocytozoon species (Additional file 1, Table S1). All these 24 sub data sets were obtained from the nucleotide alignment by simply discarding the relevant sequences without renewed aligning.
Finally, up to 51 malaria parasites were added to the initial selection of 33 taxa. We aligned all 84 taxa following the alignment procedure described above. Note that for CytB, 38 of the 51 additional sequences were partial CDS (see Additional file 8, Table S6). In order to retain more sites, complete and partial CytB genes were filtered separately with Gblocks, and then manually reassembled into an 84-gene alignment. All 6 additions of taxa to the 33-taxon data set were obtained from the whole concatenation of genes of the 84 taxa by discarding the relevant sequences.
Volkman SK, Barry AE, Lyons EJ, Nielsen KM, Thomas SM, Choi M, Thakore SS, Day KP, Wirth DF, Hartl DL: Recent origin of Plasmodium falciparum from a single progenitor. Science. 2001, 293 (5529): 482-484. 10.1126/science.1059878.
Foth BJ, McFadden GI: The apicoplast: a plastid in Plasmodium falciparum and other Apicomplexan parasites. Int Rev Cytol. 2003, 224: 57-110. full_text.
Gardner MJ, Williamson DH, Wilson RJ: A circular DNA in malaria parasites encodes an RNA polymerase like that of prokaryotes and chloroplasts. Mol Biochem Parasitol. 1991, 44: 115-123. 10.1016/0166-6851(91)90227-W.
Valkiunas G: Avian malaria parasites and other haemosporidia. 2005, New York: CRC Press, Boca Raton
Levine ND: The protozoan phylum Apicomplexa. 1988, Florida: CRC Press, Boca Raton
Perkins SL, Schall JJ: A molecular phylogeny of malarial parasites recovered from cytochrome b gene sequences. J Parasitol. 2002, 88 (5): 972-978.
Manwell RD: The problem of species, with special reference to the malaria parasites. Ann Trop Med. 1936, 30: 435-439.
Garnham PCC: Malaria parasites and other haemosporidia. 1966, Oxford: Blackwell Scientific
Perkins SL: Species concepts and malaria parasites: detecting a cryptic species of Plasmodium. Proc Biol Sci. 2000, 267 (1459): 2345-2350. 10.1098/rspb.2000.1290.
Martinsen ES, Paperna I, Schall JJ: Morphological versus molecular identification of avian Haemosporidia: an exploration of three species concepts. Parasitology. 2006, 133 (3): 279-288. 10.1017/S0031182006000424.
Leclerc MC, Hugot JP, Durand P, Renaud F: Evolutionary relationships between 15 Plasmodium species from new and old world primates (including humans): an 18 S rDNA cladistic analysis. Parasitology. 2004, 129 (6): 677-684. 10.1017/S0031182004006146.
Singh B, Kim-Sung L, Matusop A, Radhakrishnan A, Shamsul SS, Cox-Singh J, Thomas A, Conway DJ: A large focus of naturally acquired Plasmodium knowlesi infections in human beings. Lancet. 2004, 363 (9414): 1017-1024. 10.1016/S0140-6736(04)15836-4.
Cox-Singh J, Singh B: Knowlesi malaria: newly emergent and of public health importance?. Trends Parasitol. 2008, 24 (9): 406-410. 10.1016/j.pt.2008.06.001.
Garamszegi LZ: Patterns of co-speciation and host switching in primate malaria parasites. Malar J. 2009, 8 (110): 1-15.
Aravind L, Iyer LM, Wellems TE, Miller LH: Plasmodium biology: genomic gleanings. Cell. 2003, 115 (7): 771-785. 10.1016/S0092-8674(03)01023-7.
Waters AP, Higgins DG, McCutchan TF: Plasmodium falciparum appears to have arisen as a result of lateral transfer between avian and human hosts. Proc Natl Acad Sci USA. 1991, 88 (8): 3140-3144. 10.1073/pnas.88.8.3140.
Waters AP, Higgins DG, McCutchan TF: Evolutionary relatedness of some primate models of Plasmodium. Mol Biol Evol. 1993, 10 (4): 914-923.
Escalante AA, Ayala FJ: Phylogeny of the malarial genus Plasmodium, derived from rRNA gene sequences. Proc Natl Acad Sci USA. 1994, 91 (24): 11373-11377. 10.1073/pnas.91.24.11373.
Escalante AA, Ayala FJ: Evolutionary origin of Plasmodium and other Apicomplexa based on rRNA genes. Proc Natl Acad Sci USA. 1995, 92 (13): 5793-5797. 10.1073/pnas.92.13.5793.
Qari SH, Shi YP, Pieniazek NJ, Collins WE, Lal AA: Phylogenetic relationship among the malaria parasites based on small subunit rRNA gene sequences: monophyletic nature of the human malaria parasite, Plasmodium falciparum. Mol Phylogenet Evol. 1996, 6: 157-165. 10.1006/mpev.1996.0068.
Escalante AA, Freeland DE, Collins WE, Lal AA: The evolution of primate malaria parasites based on the gene encoding cytochrome b from the linear mitochondrial genome. Proc Natl Acad Sci USA. 1998, 95 (14): 8124-8129. 10.1073/pnas.95.14.8124.
Hayakawa T, Culleton R, Otani H, Horii T, Tanabe K: Big bang in the evolution of extant malaria parasites. Mol Biol Evol. 2008, 25 (10): 2233-2239. 10.1093/molbev/msn171.
Perkins SL: Molecular systematics of the three mitochondrial protein-coding genes of malaria parasites: Corroborative and new evidence for the origins of human malaria. Mitochondrial DNA. 2008, 19 (6): 471-478.
Roy S, 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.
Ollomo B, Durand P, Prugnolle F, Douzery EJP, Arnathau C, Nkoghe D, Leroy E, Renaud F: A new malaria agent in African hominids. PLoS Pathog. 2009, 5 (5): e1000446-10.1371/journal.ppat.1000446.
Rich SM, Leendertz FH, Xu G, Lebreton M, Djoko CF, Aminake MN, Takang EE, Diffo JL, Pike BL, Rosenthal BM, Formenty P, Boesch C, Ayala FJ, Wolfe ND: The origin of malignant malaria. Proc Natl Acad Sci USA. 2009, 106 (35): 14902-14907. 10.1073/pnas.0907740106.
Prugnolle F, Durand P, Neel C, Ollomo B, Ayala FJ, Arnathau C, Etienne L, Mpoudi-Ngole E, Nkoghe D, Leroy E, Delaporte E, Peeters M, Renaud F: African great apes are natural hosts of multiple related malaria species, including Plasmodium falciparum. Proc Natl Acad Sci USA. 2010, 107 (4): 1458-1463. 10.1073/pnas.0914440107.
Krief S, Escalante AA, Pacheco MA, Mugisha L, André C, Halbwax M, Fischer A, Krief JM, Kasenene JM, Crandfield M, Cornejo OE, Chavatte JM, Lin C, Letourneur F, Grüner AC, McCutchan TF, Rénia L, Snounou G: On the diversity of malaria parasites in African apes and the origin of Plasmodium falciparum from Bonobos. PLoS Pathog. 2010, 6 (2): e1000765-10.1371/journal.ppat.1000765.
McCutchan TF, Kissinger JC, Touray MG, Rogers MJ, Li J, Sullivan M, Braga EM, Krettli AU, Miller LH: Comparison of circumsporozoite proteins from avian and mammalian malarias: biological and phylogenetic implications. Proc Natl Acad Sci USA. 1996, 93 (21): 11889-11894. 10.1073/pnas.93.21.11889.
Escalante AA, Goldman IF, De-Rijk P, De-Wachter R, Collins WE, Qari SH, Lal AA: Phylogenetic study of the genus Plasmodium based on the secondary structure-based alignment of the small subunit ribosomal RNA. Mol Biochem Parasitol. 1997, 90: 317-321. 10.1016/S0166-6851(97)00121-7.
Rathore D, Wahl AM, Sullivan M, McCutchan TF: A phylogenetic comparison of gene trees constructed from plastid, mitochondrial and genomic DNA of Plasmodium species. Mol Biochem Parasitol. 2001, 114: 89-94. 10.1016/S0166-6851(01)00241-9.
Kedzierski L, Escalante AA, Isea R, Black CG, Barnwell JW, Coppel RL: Phylogenetic analysis of the genus Plasmodium based on the gene encoding adenylosuccinate lyase. Infect Genet Evol. 2002, 1 (4): 297-301. 10.1016/S1567-1348(02)00031-X.
Kissinger JC, Souza PC, Soares CO, Paul R, Wahl AM, Rathore D, McCutchan TF, Krettli AU: Molecular phylogenetic analysis of the avian malarial parasite Plasmodium Novyella juxtanucleare. J Parasitol. 2002, 88 (4): 769-773.
Vargas-Serrato E, Corredor V, Galinski MR: Phylogenetic analysis of CSP and MSP-9 gene sequences demonstrates the close relationship of Plasmodium coatneyi and Plasmodium knowlesi. Infect Genet Evol. 2003, 3: 67-73. 10.1016/S1567-1348(03)00007-8.
Hagner SC, Misof B, Maier WA, Kampen H: Bayesian analysis of new and old malaria parasite DNA sequence data demonstrates the need for more phylogenetic signal to clarify the descent of Plasmodium falciparum. Parasitol Res. 2007, 101 (3): 493-503. 10.1007/s00436-007-0499-6.
Martinsen ES, Perkins SL, Schall JJ: A three-genome phylogeny of malaria parasites (Plasmodium and closely related genera): Evolution of life-history traits and host switches. Mol Phylogenet Evol. 2008, 47: 261-273. 10.1016/j.ympev.2007.11.012.
Perkins SL, Sarkar IN, Carter R: The phylogeny of rodent malaria parasites: simultaneous analysis across three genomes. Infect Genet Evol. 2007, 7: 74-83. 10.1016/j.meegid.2006.04.005.
Swofford DL: PAUP* 4.0 beta 5: Phylogenetic Analysis Using Parsimony and Other Methods. 2001, Sinauer Associates Sunderland, MA
MrAIC v. 1.4.4. [http://www.abc.se/~nylander/]
Akaike H: Information theory and an extension of the maximum likelihood principle. Proceedings 2nd International Symposium on Information Theory. Edited by: Petrov C. 1973, Budapest, Akademia Kiado, 267-281.
Lanave C, Preparata G, Saccone C, Serio G: A new method for calculating evolutionary substitution rates. J Mol Evol. 1984, 20: 86-93. 10.1007/BF02101990.
Yang Z: Maximum-likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods. J Mol Evol. 1994, 39 (3): 306-314. 10.1007/BF00160154.
Lartillot N, Lepage T, Blanquart S: PhyloBayes 3. A Bayesian software package for phylogenetic reconstruction and molecular dating. Bioinformatics. 2009, 25 (17): 2286-2288. 10.1093/bioinformatics/btp368.
Abascal F, Zardoya R, Posada D: ProtTest: Selection of bestfit models of protein evolution. Bioinformatics. 2005, 21 (9): 2104-2105. 10.1093/bioinformatics/bti263.
Jones DT, Taylor WR, Thomton JM: The rapid generation of mutation data matrices from protein sequences. Computational Applied Bioscience. 1992, 8 (3): 275-282.
Douady CJ, Delsuc F, Boucher Y, Doolittle WF, Douzery EJP: Comparison of Bayesian and maximum likelihood bootstrap measures of phylogenetic reliability. Mol Biol Evol. 2003, 20 (2): 248-254. 10.1093/molbev/msg042.
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate Maximum-Likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst Biol. 2010, 59 (3): 307-321. 10.1093/sysbio/syq010.
Lartillot N, Philippe H: A Bayesian mixture model for across-site heterogeneities in the amino-acid replacement process. Mol Biol Evol. 2004, 21 (6): 1095-1109. 10.1093/molbev/msh112.
Rodriguez-Ezpeleta N, Brinkmann H, Roure B, Lartillot N, Lang BF, Philippe H: Detecting and overcoming systematic errors in genome-scale phylogenies. Syst Biol. 2007, 56 (3): 389-399. 10.1080/10635150701397643.
Blanquart S, Lartillot N: A Bayesian compound stochastic process for modeling nonstationary and nonhomogeneous sequence evolution. Mol Biol Evol. 2006, 23 (11): 2058-2071. 10.1093/molbev/msl091.
Blanquart S, Lartillot N: A site- and time-heterogeneous model of amino-acid replacement. Mol Biol Evol. 2008, 25 (5): 842-858. 10.1093/molbev/msn018.
Foster PG: Modeling Compositional Heterogeneity. Systematic Biology. 2004, 53 (3): 485-495. 10.1080/10635150490445779.
Dávalos LM, Perkins SL: Saturation and base composition bias explain phylogenomic conflict in Plasmodium. Genomics. 2008, 91 (5): 433-442.
Heath TA, Hedtke SM, Hillis DM: Taxon sampling and the accuracy of phylogenetic analyses. J Syst Evol. 2008, 46 (3): 239-257.
Anisimova M, Bielawski JP, Yang Z: Accuracy and power of the likelihood ratio test in detecting adaptive molecular evolution. Mol Biol Evol. 2001, 18 (8): 1585-1592.
Felsenstein J: Cases in which parsimony or compatibility method will be positively misleading. Syst Zool. 1978, 27 (4): 401-410. 10.2307/2412923.
Hughes AL, Verra F: Malaria Parasite Sequences from Chimpanzee Support the Co-Speciation Hypothesis for the Origin of Virulent Human Malaria (Plasmodium falciparum). Mol Phylogenet Evol. 2010, 57: 135-143. 10.1016/j.ympev.2010.06.004.
Delsuc F, Brinkmann H, Philippe H: Phylogenomics and the reconstruction of the tree of life. Nat Rev Genet. 2005, 6 (5): 361-375. 10.1038/nrg1603.
Vaidya AB, Akella R, Suplick K: Sequences similar to genes for two mitochondrial proteins and portions of ribosomal RNA in tandemly arrayed 6-kilobase-pair DNA of a malarial parasite. Mol Biochem Parasitol. 1989, 35 (2): 97-107. 10.1016/0166-6851(89)90112-6.
Conway DJ, Fanello C, Lloyd JM, Al-Joubori BM, Baloch AH, Somanath SD, Roper C, Oduola AM, Mulder B, Povoa MM, Singh B, Thomas AW: Origin of Plasmodium falciparum malaria is traced by mitochondrial DNA. Mol Biochem Parasitol. 2000, 111: 163-171. 10.1016/S0166-6851(00)00313-3.
Joy DA, Feng X, Mu J, Furuya T, Chotivanich K, Krettli AU, Ho M, Wang A, White NJ, Suh E, Beerli P, Su XZ: Early origin and recent expansion of Plasmodium falciparum. Science. 2003, 300 (5617): 318-321. 10.1126/science.1081449.
Beadell JS, Fleischer R: A restriction enzyme-based assay to distinguish between avian Haemosporidians. J Parasitol. 2005, 91 (3): 683-685. 10.1645/GE-3412RN.
Jongwutiwes S, Putaporntip C, Iwasaki T, Ferreira MU, Kanbara H, Hughes AL: Mitochondrial genome sequences support ancient population expansion in Plasmodium vivax. Mol Biol Evol. 2005, 22 (8): 1733-1739. 10.1093/molbev/msi168.
Mu J, Joy DA, Duan J, Huang Y, Carlton J, Walker J, Barnwell J, Beerli P, Charleston MA, Pybus OG, Su XZ: Host switch leads to emergence of Plasmodium vivax malaria in humans. Mol Biol Evol. 2005, 22 (8): 1686-1693. 10.1093/molbev/msi160.
Omori S, Sato Y, Isobe T, Yukawa M, Murata K: Complete nucleotide sequences of the mitochondrial genomes of two avian malaria protozoa, Plasmodium gallinaceum and Plasmodium juxtanucleare. Parasitol Res. 2007, 100 (3): 661-664. 10.1007/s00436-006-0333-6.
Omori S, Sato Y, Hirakawa S, Isobe T, Yukawa M, Murata K: Two extra chromosomal genomes of Leucocytozoon caulleryi; complete nucleotide sequences of the mitochondrial genome and existence of the apicoplast genome. Parasitol Res. 2008, 103 (4): 953-957. 10.1007/s00436-008-1083-4.
Feagin JE: The extrachromosomal DNAs of apicomplexan parasites. Annu Rev Microbiol. 1994, 48: 81-104. 10.1146/annurev.mi.48.100194.000501.
Wilson RJ, Williamson DH: Extrachromosomal DNA in the Apicomplexa. Microbiol Mol Biol Rev. 1997, 61: 1-16.
Aldritt SM, Joseph JT, Wirth DF: Sequence identification of cytochrome b in Plasmodium gallinaceum. Mol Cell Biol. 1989, 9 (9): 3614-3620.
Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17 (4): 540-552.
Tuffley C, Steel M: Modeling the covarion hypothesis of nucleotide substitution. Mathematical Biosciences. 1998, 147: 63-91. 10.1016/S0025-5564(97)00081-3.
Adachi J, Hasegawa M: Model of amino acid substitution in proteins encoded by mitochondrial DNA. J Mol Evol. 1996, 42 (4): 459-468. 10.1007/BF02498640.
Le SQ, Lartillot N, Gascuel O: Phylogenetic mixture models for proteins. Philos Trans R Soc Lond B Biol Sci. 2008, 363 (1512): 3965-76. 10.1098/rstb.2008.0180.
Jeffroy O, Brinkmann H, Delsuc F, Philippe H: Phylogenomics: the beginning of incongruence?. Trends Genet. 2006, 22 (4): 225-231. 10.1016/j.tig.2006.02.003.
Anisimova M, Gascuel O: Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst Biol. 2006, 55 (4): 539-552. 10.1080/10635150600755453.
Cheesman S, Tanabe K, Sawai H, O'Mahony E, Carter R: Strain-specific immunity may drive adaptive polymorphism in the merozoite surface protein 1 of the rodent malaria parasite Plasmodium chabaudi. Infect Genet Evol. 2009, 9 (2): 248-255. 10.1016/j.meegid.2008.12.002.
Hall N, Karras M, Raine JD, Carlton JM, Kooij TW, Berriman M, Florens L, Janssen CS, Pain A, Christophides GK, James K, Rutherford K, Harris B, Harris D, Churcher C, Quail MA, Ormond D, Doggett J, Trueman HE, Mendoza J, Bidwell SL, Rajandream MA, Carucci DJ, Yates JR, Kafatos FC, Janse CJ, Barrell B, Turner CM, Waters AP, Sinden RE: A comprehensive survey of the Plasmodium life cycle by genomic, transcriptomic, and proteomic analyses. Science. 2005, 307 (5706): 82-86. 10.1126/science.1103717.
Seethamchai S, Putaporntip C, Malaivijitnond S, Cui L, Jongwutiwes S: Malaria and Hepatocystis species in wild macaques, southern Thailand. Am J Trop Med Hyg. 2008, 78 (4): 646-653.
Martinsen ES, Waite JL, Schall JJ: Morphologically defined subgenera of Plasmodium from avian hosts: test of monophyly by phylogenetic analysis of two mitochondrial genes. Parasitology. 2007, 134 (4): 483-490. 10.1017/S0031182006001922.
Dutheil J, Gaillard S, Bazin E, Glffémin S, Ranwez V, Galtier N, Belkhir K: Bio++: a set of C++ libraries for sequence analysis, phylogenetics, molecular evolution and population genetics. BMC Bioinformatics. 2006, 7 (188): 1-6.
We wish to thank Emeric Sevin, Martin Taylor, Frédéric Delsuc and Nicolas Lartillot for helpful comments on this manuscript, Julien Dutheil and Vincent Ranwez for sharing code, Liliana Dávalos for sharing data, and François Renaud for helpful discussions about Haemosporidia evolution. We wish to thank the Associate Editor, Hervé Philippe, who requested the correction of sequence frameshifts during the processing of this manuscript. We thank Nick Goldman and the European Molecular Biology Laboratory (EMBL), as well as Nicolas Galtier and the European Research Council grant 232971 (ERC PopPhyl), who supported the continuation of this work during the post-doctoral positions they offered to SB. We thank the EMBL-European Bioinformatics Institute (EBI) for the computational facilities we used in this study. This work was supported by the French "Agence Nationale de la Recherche", ANR PlasmoExplore.
SB designed and conducted all the experiments. OG provided guidance throughout the study. Both authors contributed to the writing of the paper.
Electronic supplementary material
Additional file 1: . Accession numbers of 33 complete mitochondrial genomes of Haemosporidian parasites, parasite names, and host names retrieved from NCBI annotations ( b : host names complemented from Leclerc et al. 2004 ). "P.": Plasmodium species, "Ha.": Haemoproteus species, "Pa.": Parahaemoproteus species, "L.": Leucocytozoon species. (PDF 64 KB)
Additional file 2: . Saturation plot of codon positions of the 33 taxa and 3 concatenated genes data set, computed with a Bio++ script . Each dot represents the comparison of the similarity distance (y coordinate) versus the tree distance (x coordinate), for a pair of taxa. Tree branch lengths were estimated under the GTR nt + Γ4 model (PhyloBayes 3.0, ), for 3 data sets corresponding to each codon position, and using the tree topology estimated from the whole nucleotide data set (Figure 2). (PDF 55 KB)
Additional file 3: . Dependency of clade support on codon positions ("Cod. pos."), amino-acid translation and on assumptions of various probabilistic models of substitution (GTR nt , GTR aa and JTT: single matrix model, CAT: site heterogeneous mixture model, BP: time heterogeneous model). Rates across site model components are defined as + Γ4 + I under maximum-likelihood (ML) and as + Γ4 under Bayesian (BI) methods. Cells display support as [PP, SH, BS], with PP: posterior probability (BI), SH: Shimodaira-Hasegawa-like support (ML), and BS: bootstrap support (ML). "*": not applicable. Main lineages of mammal parasites are defined according to their host preference: "Rodent", "Primate" and "Great Ape" (see Additional file 1, Table S1). (PDF 104 KB)
Additional file 4: . Cross-validation estimations of the fit of Bayesian models to the 33 taxa and the three concatenated gene data sets. Models applied to the nucleotide and amino-acid data set are compared to the best fitting ML models, GTR nt + Γ4 and JTT + Γ4, respectively. Models are defined according to their components. Substitution model: GTR nt , GTR aa , MtREV and JTT, exchange rate parameters; CAT, UL2 and UL3, site heterogeneous mixture models. Rates across sites models: + Γ4, discretized gamma rates (Yang 1994, ); +COV, covarion model (Tuffley and Steel 1998, ). (PDF 68 KB)
Supplementary Table S4, p-values of posterior predictive tests performed on the 33 taxa and 3 mitochondrial gene data sets
Additional file 5: . Data sets were analyzed under various probabilistic models of substitution (GTR nt , GTR aa and JTT: single matrix model, CAT: site heterogeneous mixture model, BP: time heterogeneous model, + Γ4: Rates across site model component). Posterior predictive test "Composition" measures compositional heterogeneity across taxa, "Site Diversity" and "Homoplasy" measure the level of saturation of the phylogenetic signal. "Cod. pos.": codon positions. "*": not applicable. (PDF 99 KB)
Additional file 6: CAT + GTR + Γ 4 . Bayesian phylogenetic reconstruction using PhyloBayes 3.0 . The CAT + GTR nt + Γ4 substitution model was applied to first and second codon positions of the 33 taxa data set. P. falciparum and 2 of its relatives infecting great ape hosts, P. reichenowi and P. gaboni, formed a monophyletic clade with 3 rodent parasites, P. yoelii, P. berghei and P. chabaudi (posterior probability PP = 0.92). Posterior probabilities equal to 1 were removed. (PDF 22 KB)
Additional file 7: . All codon positions were analyzed under GTR + Γ4 and GTR + Γ4 + I models, for Bayesian and ML methods, respectively. Addition or removal of taxa to the complete nucleotide data set comprising 33 taxa and 3 concatenated genes, 3308 sites. Phylogenetic analyses were performed under models GTR nt + Γ4 + I and GTR nt + Γ4, for maximum-likelihood (ML) and Bayesian (BI) methods, respectively. Cells display support as [PP, SH, BS], with PP: posterior probability (BI), SH: Shimodaira-Hasegawa-like support (ML), and BS: bootstrap support (ML). "*": not applicable. ( a ): PP and BS are summed over various positions of Hepatocystis species. Main lineages of mammal parasites are defined according to their host preference: "Rodent", "Primate" and "Great Ape" (see Additional file 1 and 8, Tables S1 and S6). "-" removal of species. "+" addition of species. "P. fal.": P. falciparum; "P. rei.": P. reichenowi; "P. gab.": P. gaboni; "P. yoe.": P. yoelii; "P. ber.": P. berghei; "P. cha.": P. chabaudi; "Hum": human primate parasites P. malariae and P. ovale; "Afr.": African primate parasites P. gonderi and P. DAJ-2004; "Asi": 10 Asian primate parasites; "Pla.": Plasmodium species infecting saurian hosts; "Hae.": Haemoproteus and Parahaemoproteus species; "Leu.": Leucocytozoon species; "Haemo." Haemosporidian species. (PDF 79 KB)
Additional file 8: . Accession numbers of 51, 41 and 1 additional CytB, Cox1 and Cox3 genes, respectively, and parasite and host names (1 partial CytB genes). References: ( a ) Perkins and Schall (2002) ; ( b ) Perkins et al. (2007) ; ( c ) Cheesman et al. (2009) ; ( d ) Hall et al. (2005) ; ( e ) Escalante et al. (1998) ; ( f ) Seethamchai et al. (2008) ; ( g ) Martinsen et al. (2008) ; ( h ) Martinsen et al. (2007) ; ( i ) CytB + Cox1 + Cox3, Perkins (2008) ; ( j ) Rich et al. (2009) . Abbreviation: "P.": Plasmodium species, "He.": Hepatocystis species. (PDF 119 KB)
Additional file 9: . Accession numbers of 27 ClpC and 18 ASL genes. ( a ): Taxa used for the estimation of the length of the 7-taxa tree. (PDF 60 KB)
Additional file 10: . Bayesian reconstruction under the GTR nt + Γ4 model. Edges with PP < 0.9 were collapsed, and PP = 1 are not shown. The ASL phylogeny is not congruent with a monophyly of mammal malaria parasites. (PDF 16 KB)
Additional file 11: . Bayesian reconstruction under the GTR nt + Γ4 model. Edges with PP < 0.9 were collapsed, and PP = 1 are not shown. The ClpC phylogeny supports the monophyly of mammalian malaria parasites (PP = 0.99). (PDF 16 KB)
Additional file 12: . Each of the 9 taxon samples comprised the three most distantly related Leucocytozoon species, the 14 primate parasites, and two single representatives of parasites of great apes and of rodents, respectively. In each case, 14 taxa were removed from the complete nucleotide data-set comprising 33 taxa and 3 concatenated genes, 3308 sites. Phylogenetic analyses were performed under models GTR nt + Γ4 + I and GTR nt + Γ4, for maximum likelihood (ML) and Bayesian (BI) methods, respectively. Cells display support as follows: [PP, SH, BS], with PP: posterior probability (BI), SH: Shimodaira-Hasegawa-like support ("*": not applicable, ML), and BS: bootstrap support (ML). Main lineages of mammal parasites are defined according to their host preference: "Rodent", "Primate" and "Great Ape" (see Additional file 1, Table S1). "P. fal.": P. falciparum; "P. rei.": P. reichenowi; "P. gab.": P. gaboni; "P. yoe.": P. yoelii; "P. ber.": P. berghei; "P. cha.": P. chabaudi. (PDF 75 KB)
About this article
Cite this article
Blanquart, S., Gascuel, O. Mitochondrial genes support a common origin of rodent malaria parasites and Plasmodium falciparum's relatives infecting great apes. BMC Evol Biol 11, 70 (2011). https://doi.org/10.1186/1471-2148-11-70
- Codon Position
- Plasmodium Species
- Compositional Heterogeneity
- Rodent Malaria Parasite
- Model Violation