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

Despite substantial work, the phylogeny of malaria parasites remains debated. The matter is complicated by concerns about patterns of evolution in potentially strongly selected genes as well as the extreme AT bias of some Plasmodium genomes. Particularly contentious has been the position of the most virulent human parasite Plasmodium falciparum, whether grouped with avian parasites or within a larger clade of mammalian parasites. Here, we study 3 classes of rare genomic changes, as well as the sequences of mitochondrial ribosomal RNA (rRNA) genes. We report 3 lines of support for a clade of mammalian parasites: 1) we find no instances of spliceosomal intron loss in a hypothetical ancestor of P. falciparum and the avian parasite Plasmodium gallinaceum, suggesting against a close relationship between those species; 2) we find 4 genomic mitochondrial indels supporting a mammalian clade, but none grouping P. falciparum with avian parasites; and 3) slowly evolving mitochondrial rRNA sequences support a mammalian parasite clade with 100% posterior probability. We further report a large deletion in the mitochondrial large subunit rRNA gene, which suggests a subclade including both African and Asian parasites within the clade of closely related primate malarias. This contrasts with previous studies that provided strong support for separate Asian and African clades, and reduces certainty about the historical and geographic origins of Plasmodium vivax. Finally, we find a lack of synapomorphic gene losses, suggesting a low rate of ancestral gene loss in Plasmodium.     

Remarkable interkingdom conservation of intron positions and massive,lineage-specific intron loss and gain in eukaryotic evolution

Sequencing of eukaryotic genomes allows one to address major evolutionary problems, such as the evolution of gene structure. We compared the intron positions in 684 orthologous gene sets from 8 complete genomes of animals, plants, fungi, and protists and constructed parsimonious scenarios of evolution of the exon-intron structure for the respective genes. Approximately one-third of the introns in the malaria parasite Plasmodium falciparum are shared with at least one crown group eukaryote; this number indicates that these introns have been conserved through >1.5 billion years of evolution that separate Plasmodium from the crown group. Paradoxically, humans share many more introns with the plant Arabidopsis thaliana than with the fly or nematode. The inferred evolutionary scenario holds that the common ancestor of Plasmodium and the crown group and, especially, the common ancestor of animals, plants, and fungi had numerous introns. Most of these ancestral introns, which are retained in the genomes of vertebrates and plants, have been lost in fungi, nematodes, arthropods, and probably Plasmodium. In addition, numerous introns have been inserted into vertebrate and plant genes, whereas, in other lineages, intron gain was much less prominent.     

The genome of model malaria parasites, and comparative genomics

The field of comparative genomics of malaria parasites has recently come of age with the completion of the whole genome sequences of the human malaria parasite Plasmodium falciparum and a rodent malaria model, Plasmodium yoelii yoelii. With several other genome sequencing projects of different model and human malaria parasite species underway, comparing genomes from multiple species has necessitated the development of improved informatics tools and analyses. Results from initial comparative analyses reveal striking conservation of gene synteny between malaria species within conserved chromosome cores, in contrast to reduced homology within subtelomeric regions, in line with previous findings on a smaller scale. Genes that elicit a host immune response are frequently found to be species-specific, although a large variant multigene family is common to many rodent malaria species and Plasmodium vivax. Sequence alignment of syntenic regions from multiple species has revealed the similarity between species in coding regions to be high relative to non-coding regions, and phylogenetic footprinting studies promise to reveal conserved motifs in the latter. Comparison of non-synonymous substitution rates between orthologous genes is proving a powerful technique for identifying genes under selection pressure, and may be useful for vaccine design. This is a stimulating time for comparative genomics of model and human malaria parasites, which promises to produce useful results for the development of antimalarial drugs and vaccines.     

Plasmodium yoelii: experimental evidences for the conserved epitopes between mouse and human malaria parasite, Plasmodium falciparum

Bioinformatic analyses of gene homologues have revealed functionally conserved epitopes between human and rodent malaria parasites. Here, we present experimental evidence for the presence of functionally and antigenically conserved domains between Plasmodium falciparum and Plasmodium yoelii asexual blood-stages. Merozoite released soluble proteins (MRSPs) from both P. falciparum and P. yoelii bound to heterologous mouse or human red blood cells, respectively. The presence of conserved antigenic epitopes between the two species of parasites was evident by the inhibitory effect of antibodies, developed against P. yoelii in convalescent mice, on P. falciparum growth and merozoite reinvasion in vitro. Furthermore, mice immunized with P. falciparum MRSPs were protected from infection by a P. yoelii challenge. These data indicate that different species of Plasmodium contain antigenically conserved interspecies domains, which are immunogenic and, thus constitute a potential novel antigen source for vaccine development and testing using a mouse model.     

Comparative kinomics of Plasmodium organisms: unity in diversity

Phosphorylation by protein kinases is a very common and crucial process in many signal transduction pathways in eukaryotes. This review describes comparative protein kinase analysis of two apicomplexa Plasmodium falciparum (3D7 strain) and Plasmodium yoelii yoelii (17XNL strain) which are causative agents of malaria in human and African rat respectively. Sensitive bioinformatics techniques enable identification of 82 and 60 putative protein kinases in P. falciparum and P. yoelii yoelii respectively and these sequences could be classified into known subfamilies of protein kinases. The most populated kinase subfamilies in both the plasmodium species correspond to CAMK and CMGC groups. Analysis of domain architectures enables detection of uncommon domain organization in kinases of both the organisms such as kinase domain tethered to EF hands as well as PH domain. Components of MAPK signaling pathway is not well conserved in plasmodium organisms. Such observations suggest that plasmodium protein kinases are highly divergent from other eukaryotes. A transmembrane kinase with 6 membrane spanning segments in P. falciparum seems to have no orthologue in P. yoelii yoelii. 19 P. falciparum kinases have been found to cluster separately from P. yoelii yoelii kinases and hence these kinases are unique to P. falciparum genome. Only 28 orthologous pairs of kinases seem to be present between these two plasmodium organisms. Comparative kinome analysis of two plasmodium species has thus provided clues to the function of many protein kinases based upon their classification and domain organization and also implicate marked differences even between two plasmodium organisms.     

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.     

Sequence analysis of the Rhop-3 gene of Plasmodium yoelii

The 110 kDa/Rhop-3 rhoptry protein of Plasmodium falciparum is non-covalently associated with two other proteins, the 140 kDa Rhop-1 and the 130 kDa Rhop-2. cDNAs encoding Rhop-3 from Plasmodium yoelii were isolated using rhoptry-specific antisera from Plasmodium falciparum, P. yoelii, and Plasmodium chabaudi. The cDNAs encoded peptides with partial homology to the C-terminal region (residues 541-861) of P. falciparum Rhop-3. Core regions of homology to the P. falciparum gene will be useful in determining the biological role of Rhop-3 and its potential as a vaccine candidate for malaria.     

Chimeric epitopes delivered by polymeric synthetic linear peptides induce protective immunity to malaria

Polymeric linear peptide chimeras (LPCs) that incorporate Plasmodium vivax promiscuous T cell epitopes and the P. falciparum circumsporozoite protein B cell epitope have been shown to induce a high level of immunogenicity and overcome genetic restriction when tested as vaccine immunogens in BALB/c mice. The present study evaluates the biological relevance of several LPCs using a well characterized rodent malaria model. Polymeric peptide constructs based on P. berghei and P. yoelii sequences, and orthologous to the human malaria sequences included in the original LPCs, were designed and tested for immunogenicity in mice of different H-2 haplotypes. We demonstrate that robust immune responses are induced and that peptides containing the orthologous rodent Plasmodium sequences exhibited similar immunogenic capabilities. Unique to this report, we show that LPCs can also prime MHC class I-restricted cytotoxic T lymphocytes (CTLs) and, most relevantly, that a peptide construct prototype incorporating single B, T and CTL epitopes induced protection against an experimental challenge with P. berghei or P. yoelii sporozoites. Collectively, these results suggest that polymeric polypeptide chimeras can be used as a platform to deliver subunit vaccines.     

Comparative histopathology of mice infected with the 17XL and 17XNL strains of Plasmodium yoelii

Plasmodium yoelii 17XL was used to investigate the mechanism of Plasmodium falciparum-caused cerebral malaria, although its histological effect on other mouse organs is still unclear. Here, histological examination was performed on mice infected with P. yoelii 17XL; the effect of P. yoelii 17XL infection on anemia and body weight loss, as well as its lesions in the brain, liver, kidney, lung, and spleen, also was investigated. Plasmodium yoelii 17XL-infected red blood cells were sequestered in the microcirculation of the brain and in the kidney. Compared with the nonlethal P. yoelii 17XNL strain, infection by P. yoelii 17XL caused substantial pulmonary edema, severe anemia, and significant body weight loss. Although P. yoelii 17XNL and 17XL produced a similar focal necrosis in the mouse liver, infection of P. yoelii 17XL induced coalescing of red and white pulp. Mortality caused by P. yoelii 17XL may be due to cerebral malaria, as well as respiratory distress syndrome and severe anemia. Plasmodium yoelii 17XL-infected rodent malaria seems to be a useful model for investigating severe malaria caused by P. falciparum.     

A Plasmodium whole-genome synteny map: indels and synteny breakpoints as foci for species-specific genes

Whole-genome comparisons are highly informative regarding genome evolution and can reveal the conservation of genome organization and gene content, gene regulatory elements, and presence of species-specific genes. Initial comparative genome analyses of the human malaria parasite Plasmodium falciparum and rodent malaria parasites (RMPs) revealed a core set of 4,500 Plasmodium orthologs located in the highly syntenic central regions of the chromosomes that sharply defined the boundaries of the variable subtelomeric regions. We used composite RMP contigs, based on partial DNA sequences of three RMPs, to generate a whole-genome synteny map of P. falciparum and the RMPs. The core regions of the 14 chromosomes of P. falciparum and the RMPs are organized in 36 synteny blocks, representing groups of genes that have been stably inherited since these malaria species diverged, but whose relative organization has altered as a result of a predicted minimum of 15 recombination events. P. falciparum-specific genes and gene families are found in the variable subtelomeric regions (575 genes), at synteny breakpoints (42 genes), and as intrasyntenic indels (126 genes). Of the 168 non-subtelomeric P. falciparum genes, including two newly discovered gene families, 68% are predicted to be exported to the surface of the blood stage parasite or infected erythrocyte. Chromosomal rearrangements are implicated in the generation and dispersal of P. falciparum-specific gene families, including one encoding receptor-associated protein kinases. The data show that both synteny breakpoints and intrasyntenic indels can be foci for species-specific genes with a predicted role in host-parasite interactions and suggest that, besides rearrangements in the subtelomeric regions, chromosomal rearrangements may also be involved in the generation of species-specific gene families. A majority of these genes are expressed in blood stages, suggesting that the vertebrate host exerts a greater selective pressure than the mosquito vector, resulting in the acquisition of diversity.     

Expression of human CD81 differently affects host cell susceptibility to malaria sporozoites depending on the Plasmodium species

Plasmodium sporozoites can enter host cells by two distinct pathways, either through disruption of the plasma membrane followed by parasite transmigration through cells, or by formation of a parasitophorous vacuole (PV) where the parasite further differentiates into a replicative exo-erythrocytic form (EEF). We now provide evidence that following invasion without PV formation, transmigrating Plasmodium falciparum and Plasmodium yoelii sporozoites can partially develop into EEFs inside hepatocarcinoma cell nuclei. We also found that rodent P. yoelii sporozoites can infect both mouse and human hepatocytes, while human P. falciparum sporozoites infect human but not mouse hepatocytes. We have previously reported that the host tetraspanin CD81 is required for PV formation by P. falciparum and P. yoelii sporozoites. Here we show that expression of human CD81 in CD81-knockout mouse hepatocytes is sufficient to confer susceptibility to P. yoelii but not P. falciparum sporozoite infection, showing that the narrow P. falciparum host tropism does not rely on CD81 only. Also, expression of CD81 in a human hepatocarcinoma cell line is sufficient to promote the formation of a PV by P. yoelii but not P. falciparum sporozoites. These results highlight critical differences between P. yoelii and P. falciparum sporozoite infection, and suggest that in addition to CD81, other molecules are specifically required for PV formation during infection by the human malaria parasite.     

Widespread intron loss suggests retrotransposon activity in ancient apicomplexans

Several facets of spliceosomal intron in apicomplexans remain mysterious. First, intron numbers vary across species by 2 orders of magnitude, indicating massive intron loss and/or gain. Second, previous studies have shown very different evolutionary patterns over different timescales, with 100-fold higher rates of intron loss/gain between genera than within genera. Third, the timing and dynamics of nearly complete intron loss in Cryptosporidium species, as well as reasons for retention of the few remaining introns, remain unknown. We compared intron positions in 785 orthologous genes between 3 moderate to intron-rich apicomplexan species. We estimate that the Theileria-Plasmodium ancestor had 4.5 times as many introns as modern Plasmodium species and 38% more than modern Theileria species, and that subsequent intron losses have outnumbered intron gains by 5.8 to 1 in Theileria and by some 56 to 1 in Plasmodium. Several patterns suggest that these intron losses occurred by recombination with reverse-transcribed mRNAs. Intriguingly, this finding suggests significant retrotransposon activity in the lineages leading to both Theileria and Plasmodium, in contrast to the dearth of known retrotransposons and intron loss within modern species from both genera. We also compared genomes from Cryptosporidium parvum and C. hominis and found no evidence of ongoing intron loss, nor of intron gain. By contrast, Cryptosporidium introns are less evolutionary conserved with Toxoplasma than are introns from other apicomplexans; thus the few remaining introns are not simply indispensable ancestral introns.     

Patterns of intron loss and gain in plants: intron loss-dominated evolution and genome-wide comparison of O. sativa and A. thaliana

Numerous previous studies have elucidated 2 surprising patterns of spliceosomal intron evolution in diverse eukaryotes over the past roughly 100 Myr. First, rates of recent intron gain in a wide variety of eukaryotic lineages have been surprisingly low, far too low to explain modern intron densities. Second, intron losses have outnumbered intron gains over a variety of lineages. For several reasons, land plants might be expected to have comparatively high rates of intron gain and thus to represent a possible exception to this pattern. However, we report several studies that indicate low rates of intron gain and an excess of intron losses over intron gains in a variety of plant lineages. We estimate that intron losses have outnumbered intron gains in recent evolution in Arabidopsis thaliana (roughly 12.6 times more losses than gains), Oryza sativa (9.8 times), the green alga Chlamydomonas reinhardtii (5.1 times), and the Bigelowiella natans nucleomorph, an enslaved green algal nucleus (2.8 times). We estimate that during recent evolution, A. thaliana and O. sativa have experienced very low rates of intron gain of around one gain per gene per 2.6-8.0 billion years. In addition, we compared 8,258 pairs of putatively orthologous A. thaliana-O. sativa genes. We found that 5.3% of introns in conserved coding regions are species-specific. Observed species-specific A. thaliana and O. sativa introns tend to be exact and to lie adjacent to each other along the gene, in a pattern suggesting mRNA-mediated intron loss. Our results underscore that low intron gain rates and intron number reduction are common features of recent eukaryotic evolution. This pattern implies that rates of intron creation were higher during earlier periods of evolution and further focuses attention on the causes of initial intron proliferation.     

The evolution of spliceosomal introns in alveolates

Many issues concerning the evolution of spliceosomal introns remain poorly understood. In this respect, the reconstruction of the evolution of introns in deep branching species such as alveolates is of special significance. In this study, we inferred the intron evolution in alveolates using 3,368 intron positions in 162 orthologs from 10 species (9 alveolates and 1 outgroup, Homo sapiens). We found that although very few intron gains and losses have occurred in Theileria and Plasmodium recently, many intron gains and losses have occurred in the evolution of alveolates. Thus, the rates of intron gain and loss in alveolates have varied greatly across time and lineage. Our results seem to support the notion that massive intron gains and losses have occurred during short episodes, perhaps coinciding with major evolutionary events.     

Identifying and characterising the Plasmodium falciparum RhopH3 Plasmodium vivax homologue

Four Plasmodium species cause malaria in humans, Plasmodium falciparum being the most widely studied to date. All Plasmodium species have paired club-shaped organelles towards their apical extreme named rhoptries that contain many lipids and proteins which are released during target cell invasion. P. falciparum RhopH3 is a rhoptry protein triggering important immune responses in patients from endemic regions. It has also been shown that anti-RhopH3 antibodies inhibit in vitro invasion of erythrocytes. Recent immunisation studies in mice with the Plasmodium yoelii and Plasmodium berghei RhopH3 P. falciparum homologue proteins found that they are able to induce protection in murine models. This study described identifying and characterising RhopH3 protein in Plasmodium vivax; it is encoded by a seven exon gene and expressed during the parasite's asexual stage. PvRhopH3 has similar processing to its homologue in P. falciparum and presents a cellular immunolocalisation pattern characteristic of rhoptry proteins.     

Genome comparison of human and non-human malaria parasites reveals species subset-specific genes potentially linked to human disease

Genes underlying important phenotypic differences between Plasmodium species, the causative agents of malaria, are frequently found in only a subset of species and cluster at dynamically evolving subtelomeric regions of chromosomes. We hypothesized that chromosome-internal regions of Plasmodium genomes harbour additional species subset-specific genes that underlie differences in human pathogenicity, human-to-human transmissibility, and human virulence. We combined sequence similarity searches with synteny block analyses to identify species subset-specific genes in chromosome-internal regions of six published Plasmodium genomes, including Plasmodium falciparum, Plasmodium vivax, Plasmodium knowlesi, Plasmodium yoelii, Plasmodium berghei, and Plasmodium chabaudi. To improve comparative analysis, we first revised incorrectly annotated gene models using homology-based gene finders and examined putative subset-specific genes within syntenic contexts. Confirmed subset-specific genes were then analyzed for their role in biological pathways and examined for molecular functions using publicly available databases. We identified 16 genes that are well conserved in the three primate parasites but not found in rodent parasites, including three key enzymes of the thiamine (vitamin B1) biosynthesis pathway. Thirteen genes were found to be present in both human parasites but absent in the monkey parasite P. knowlesi, including genes specifically upregulated in sporozoites or gametocytes that could be linked to parasite transmission success between humans. Furthermore, we propose 15 chromosome-internal P. falciparum-specific genes as new candidate genes underlying increased human virulence and detected a currently uncharacterized cluster of P. vivax-specific genes on chromosome 6 likely involved in erythrocyte invasion. In conclusion, Plasmodium species harbour many chromosome-internal differences in the form of protein-coding genes, some of which are potentially linked to human disease and thus promising leads for future laboratory research.