The origins of apicomplexan sequence innovation

Plasmodium, the Apicomplexa Outlier When It Comes to Protein Synthesis

Plasmodium is an obligate intracellular parasite that has numerous interactions with different hosts during its elaborate life cycle. This is also the case for the other parasites belonging to the same phylum Apicomplexa. In this study, we bioinformatically identified the components of the multi-synthetase complexes (MSCs) of several Apicomplexa parasites and modelled their assembly using AlphaFold2. It appears that none of these MSCs resemble the two MSCs that we have identified and characterized in Plasmodium. Indeed, tRip, the central protein involved in the association of the two Plasmodium MSCs is different from its homologues, suggesting also that the tRip-dependent import of exogenous tRNAs is not conserved in other apicomplexan parasites. Based on this observation, we searched for obvious differences that could explain the singularity of Plasmodium protein synthesis by comparing tRNA genes and amino acid usage in the different genomes. We noted a contradiction between the large number of asparagine residues used in Plasmodium proteomes and the single gene encoding the tRNA that inserts them into proteins. This observation remains true for all the Plasmodia strains studied, even those that do not contain long asparagine homorepeats.

Massive invasion of organellar DNA drives nuclear genome evolution in Toxoplasma

Toxoplasma gondii is a zoonotic protist pathogen that infects up to one third of the human population. This apicomplexan parasite contains three genome sequences: nuclear (65 Mb); plastid organellar, ptDNA (35 kb); and mitochondrial organellar, mtDNA (5.9 kb of non-repetitive sequence). We find that the nuclear genome contains a significant amount of NUMTs (nuclear integrants of mitochondrial DNA) and NUPTs (nuclear integrants of plastid DNA) that are continuously acquired and represent a significant source of intraspecific genetic variation. NUOT (nuclear DNA of organellar origin) accretion has generated 1.6% of the extant T. gondii ME49 nuclear genome-the highest fraction ever reported in any organism. NUOTs are primarily found in organisms that retain the non-homologous end-joining repair pathway. Significant movement of organellar DNA was experimentally captured via amplicon sequencing of a CRISPR-induced double-strand break in non-homologous end-joining repair competent, but not ku80 mutant, Toxoplasma parasites. Comparisons with Neospora caninum, a species that diverged from Toxoplasma ~28 mya, revealed that the movement and fixation of five NUMTs predates the split of the two genera. This unexpected level of NUMT conservation suggests evolutionary constraint for cellular function. Most NUMT insertions reside within (60%) or nearby genes (23% within 1.5 kb), and reporter assays indicate that some NUMTs have the ability to function as cis-regulatory elements modulating gene expression. Together, these findings portray a role for organellar sequence insertion in dynamically shaping the genomic architecture and likely contributing to adaptation and phenotypic changes in this important human pathogen.

Massive invasion of organellar DNA drives nuclear genome evolution in Toxoplasma

Toxoplasma gondii is a zoonotic protist pathogen that infects up to 1/3 of the human population. This apicomplexan parasite contains three genome sequences: nuclear (63 Mb); plastid organellar, ptDNA (35 kb); and mitochondrial organellar, mtDNA (5.9 kb of non-repetitive sequence). We find that the nuclear genome contains a significant amount of NUMTs (nuclear DNA of mitochondrial origin) and NUPTs (nuclear DNA of plastid origin) that are continuously acquired and represent a significant source of intraspecific genetic variation. NUOT (nuclear DNA of organellar origin) accretion has generated 1.6% of the extant T. gondii ME49 nuclear genome; the highest fraction ever reported in any organism. NUOTs are primarily found in organisms that retain the non-homologous end-joining repair pathway. Significant movement of organellar DNA was experimentally captured via amplicon sequencing of a CRISPR-induced double-strand break in non-homologous end-joining repair competent, but not ku80 mutant, Toxoplasma parasites. Comparisons with Neospora caninum, a species that diverged from Toxoplasma ~28 MY ago, revealed that the movement and fixation of 5 NUMTs predates the split of the two genera. This unexpected level of NUMT conservation suggests evolutionary constraint for cellular function. Most NUMT insertions reside within (60%) or nearby genes (23% within 1.5 kb) and reporter assays indicate that some NUMTs have the ability to function as cis-regulatory elements modulating gene expression. Together these findings portray a role for organellar sequence insertion in dynamically shaping the genomic architecture and likely contributing to adaptation and phenotypic changes in this important human pathogen.

Mobilome of Apicomplexa Parasites

Transposable elements (TEs) are mobile genetic elements found in the majority of eukaryotic genomes. Genomic studies of protozoan parasites from the phylum Apicomplexa have only reported a handful of TEs in some species and a complete absence in others. Here, we studied sixty-four Apicomplexa genomes available in public databases, using a ‘de novo’ approach to build candidate TE models and multiple strategies from known TE sequence databases, pattern recognition of TEs, and protein domain databases, to identify possible TEs. We offer an insight into the distribution and the type of TEs that are present in these genomes, aiming to shed some light on the process of gains and losses of TEs in this phylum. We found that TEs comprise a very small portion in these genomes compared to other organisms, and in many cases, there are no apparent traces of TEs. We were able to build and classify 151 models from the TE consensus sequences obtained with RepeatModeler, 96 LTR TEs with LTRpred, and 44 LINE TEs with MGEScan. We found LTR Gypsy-like TEs in Eimeria, Gregarines, Haemoproteus, and Plasmodium genera. Additionally, we described LINE-like TEs in some species from the genera Babesia and Theileria. Finally, we confirmed the absence of TEs in the genus Cryptosporidium. Interestingly, Apicomplexa seem to be devoid of Class II transposons.

Structures of plasmepsin X from Plasmodium falciparum reveal a novel inactivation mechanism of the zymogen and molecular basis for binding of inhibitors in mature enzyme

Plasmodium falciparum plasmepsin X (PfPMX), involved in the invasion and egress of this deadliest malarial parasite, is essential for its survival and hence considered as an important drug target. We report the first crystal structure of PfPMX zymogen containing a novel fold of its prosegment. A unique twisted loop from the prosegment and arginine 244 from the mature enzyme is involved in zymogen inactivation; such mechanism, not previously reported, might be common for apicomplexan proteases similar to PfPMX. The maturation of PfPMX zymogen occurs through cleavage of its prosegment at multiple sites. Our data provide thorough insights into the mode of binding of a substrate and a potent inhibitor 49c to PfPMX. We present molecular details of inactivation, maturation, and inhibition of PfPMX that should aid in the development of potent inhibitors against pepsin-like aspartic proteases from apicomplexan parasites.

Next-Generation Technologies and Systems Biology for the Design of Novel Vaccines Against Apicomplexan Parasites

Parasites of the phylum Apicomplexa are the causative agents of important diseases such as malaria, toxoplasmosis or cryptosporidiosis in humans, and babesiosis and coccidiosis in animals. Whereas the first human recombinant vaccine against malaria has been approved and recently recommended for wide administration by the WHO, most other zoonotic parasitic diseases lack of appropriate immunoprophylaxis. Sequencing technologies, bioinformatics, and statistics, have opened the "omics" era into apicomplexan parasites, which has led to the development of systems biology, a recent field that can significantly contribute to more rational design for new vaccines. The discovery of novel antigens by classical approaches is slow and limited to very few antigens identified and analyzed by each study. High throughput approaches based on the expansion of the "omics", mainly genomics and transcriptomics have facilitated the functional annotation of the genome for many of these parasites, improving significantly the understanding of the parasite biology, interactions with the host, as well as virulence and host immune response. Developments in genetic manipulation in apicomplexan parasites have also contributed to the discovery of new potential vaccine targets. The present minireview does a comprehensive summary of advances in "omics", CRISPR/Cas9 technologies, and in systems biology approaches applied to apicomplexan parasites of economic and zoonotic importance, highlighting their potential of the holistic view in vaccine development.

In silico identification of novel open reading frames in Plasmodium falciparum oocyte and salivary gland sporozoites using proteogenomics framework

Background

Plasmodium falciparum causes the deadliest form of malaria, which remains one of the most prevalent infectious diseases. Unfortunately, the only licensed vaccine showed limited protection and resistance to anti-malarial drug is increasing, which can be largely attributed to the biological complexity of the parasite’s life cycle. The progression from one developmental stage to another in P. falciparum involves drastic changes in gene expressions, where its infectivity to human hosts varies greatly depending on the stage. Approaches to identify candidate genes that are responsible for the development of infectivity to human hosts typically involve differential gene expression analysis between stages. However, the detection may be limited to annotated proteins and open reading frames (ORFs) predicted using restrictive criteria.

Methods

The above problem is particularly relevant for P. falciparum; whose genome annotation is relatively incomplete given its clinical significance. In this work, systems proteogenomics approach was used to address this challenge, as it allows computational detection of unannotated, novel Open Reading Frames (nORFs), which are neglected by conventional analyses. Two pairs of transcriptome/proteome were obtained from a previous study where one was collected in the mosquito-infectious oocyst sporozoite stage, and the other in the salivary gland sporozoite stage with human infectivity. They were then re-analysed using the proteogenomics framework to identify nORFs in each stage.

Results

Translational products of nORFs that map to antisense, intergenic, intronic, 3′ UTR and 5′ UTR regions, as well as alternative reading frames of canonical proteins were detected. Some of these nORFs also showed differential expression between the two life cycle stages studied. Their regulatory roles were explored through further bioinformatics analyses including the expression regulation on the parent reference genes, in silico structure prediction, and gene ontology term enrichment analysis.

Conclusion

The identification of nORFs in P. falciparum sporozoites highlights the biological complexity of the parasite. Although the analyses are solely computational, these results provide a starting point for further experimental validation of the existence and functional roles of these nORFs,

Proteome size reduction in Apicomplexans is linked with loss of DNA repair and host redundant pathways

Apicomplexans are alveolate parasites which include Plasmodium falciparum, the main cause of malaria, one of the world's biggest killers from infectious disease. Apicomplexans are characterized by a reduction in proteome size, which appears to result from metabolic and functional simplification, commensurate with their parasitic lifestyle. However, other factors may also help to explain gene loss such as population bottlenecks experienced during transmission, and the effect of reducing the overall genomic information content. The latter constitutes an 'informational constraint', which is proposed to exert a selective pressure to evolve and maintain genes involved in informational fidelity and error correction, proportional to the quantity of information in the genome (which approximates to proteome size). The dynamics of gene loss was examined in 41 Apicomplexan genomes using orthogroup analysis. We show that loss of genes involved in amino acid metabolism and steroid biosynthesis can be explained by metabolic redundancy with the host. We also show that there is a marked tendency to lose DNA repair genes as proteome size is reduced. This may be explained by a reduction in size of the informational constraint and can help to explain elevated mutation rates in pathogens with reduced genome size. Multiple Sequentially Markovian Coalescent (MSMC) analysis indicates a recent bottleneck, consistent with predictions generated using allele-based population genetics approaches, implying that relaxed selection pressure due to reduced population size might have contributed to gene loss. However, the non-randomness of pathways that are lost challenges this scenario. Lastly, we identify unique orthogroups in malaria-causing Plasmodium species that infect humans, with a high proportion of membrane associated proteins. Thus, orthogroup analysis appears useful for identifying novel candidate pathogenic factors in parasites, when there is a wide sample of genomes available.

Species diversity and genome evolution of the pathogenic protozoan parasite, Neospora caninum

Neospora caninum is a cyst-forming coccidian parasite of veterinary and economical significance, affecting dairy and beef cattle industries on a global scale. Comparative studies suggest that N. caninum consists of a globally dispersed, diverse population of lineages, distinguished by their geographical origin, broad host range, and phenotypic features. This viewpoint is however changing. While intraspecies diversity, and more specifically pathogenic variability, has been experimentally demonstrated in a myriad of studies, the underlying contributors and sources responsible for such diversity have remained nebulous. However, recent large-scale sequence and bioinformatics studies have aided in revealing intrinsic genetic differences distinguishing isolates of this species, that await further characterisation as causative links to virulence and pathogenicity. Furthermore, progress on N. caninum research as a non-model organism is hindered by a lack of robust, annotated genomic, transcriptomic, and proteomic data for the species, especially compared to other thoroughly studied Apicomplexa such as Toxoplasma gondii and Plasmodium species. This review explores the current body of knowledge on intra-species diversity within N. caninum. This includes the contribution of sequence variants in both coding and non-coding regions, the presence of genome polymorphic hotspots, and the identification of non-synonymous mutations. The implications of such diversity on important parasite phenotypes such as pathogenicity and population structure are also discussed. Lastly, the identification of potential virulence factors from both in-silico and next generation sequencing studies is examined, offering new insights into potential avenues for future research on neosporosis.

An immunoproteomic approach to identifying immunoreactive proteins in Leishmania infantum amastigotes using sera of dogs infected with canine visceral leishmaniasis

Visceral leishmaniasis (VL), the most severe form of leishmaniasis, is caused by Leishmania donovani and Leishmania infantum. The infected dogs with canine visceral leishmaniasis (CVL) are important reservoirs for VL in humans, so the diagnosis, treatment and vaccination of the infected dogs will ultimately decrease the rate of human VL. Proteomics and immunoproteomics techniques have facilitated the introduction of novel drug, vaccine and diagnostic targets. Our immunoproteomic study was conducted to identify new immunoreactive proteins in amastigote form of L. infantum. The strain of L. infantum (MCAN/IR/07/Moheb-gh) was obtained from CVL-infected dogs. J774 macrophage cells were infected with the L. infantum promastigotes. The infected macrophages were ruptured, and pure amastigotes were extracted from the macrophages. After protein extraction, two-dimensional gel electrophoresis was employed for protein separation followed by Western blotting. Western blotting was performed, using symptomatic and asymptomatic sera of the infected dogs with CVL. Thirteen repeatable immunoreactive spots were identified by Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Some, including prohibitin, ornithine aminotransferase, annexin A4, and apolipoprotein A-I, have been critically involved in metabolic pathways, survival, and pathogenicity of Leishmania parasites. Further investigations are required to confirm our identified immunoreactive proteins as a biomarker for CVL.

Distribution and Evolution of Peroxisomes in Alveolates (Apicomplexa, Dinoflagellates, Ciliates)

The peroxisome was the last organelle to be discovered and five decades later it is still the Cinderella of eukaryotic compartments. Peroxisomes have a crucial role in the detoxification of reactive oxygen species, the beta-oxidation of fatty acids, and the biosynthesis of etherphospholipids, and they are assumed to be present in virtually all aerobic eukaryotes. Apicomplexan parasites including the malaria and toxoplasmosis agents were described as the first group of mitochondriate protists devoid of peroxisomes. This study was initiated to reassess the distribution and evolution of peroxisomes in the superensemble Alveolata (apicomplexans, dinoflagellates, ciliates). We established transcriptome data from two chromerid algae (Chromera velia, Vitrella brassicaformis), and two dinoflagellates (Prorocentrum minimum, Perkinsus olseni) and identified the complete set of essential peroxins in all four reference species. Our comparative genome analysis provides unequivocal evidence for the presence of peroxisomes in Toxoplasma gondii and related genera. Our working hypothesis of a common peroxisomal origin of all alveolates is supported by phylogenetic analyses of essential markers such as the import receptor Pex5. Vitrella harbors the most comprehensive set of peroxisomal proteins including the catalase and the glyoxylate cycle and it is thus a promising model organism to investigate the functional role of this organelle in Apicomplexa.

Moving towards improved vaccines for Toxoplasma gondii

INTRODUCTION

Toxoplasma gondii is an intracellular parasitic protozoan that infects almost all warm-blooded animals and humans, resulting in threats to public health and economic losses. Despite continuous research efforts, there are still very few effective strategies against toxoplasmosis. In the past few years, numerous vaccination experiments have been performed to control T. gondii infection.

AREAS COVERED

In this review, the authors summarize the development of T. gondii vaccines with proper adjuvants, ranging from live or live-attenuated vaccines to protein vaccines, DNA vaccines, epitope vaccines and novel vaccines. They also highlight the challenges involved in the development of T. gondii vaccines, including specific impediments and shortcomings.

EXPERT OPINION

Moving towards the development of effective vaccines against T. gondii is not only a tedious mission but also a difficult challenge. Future studies should consider new approaches and strategies for vaccine development, particularly novel vaccines and genetic adjuvants, as well as optimizing immunization protocols and evaluation criteria.

Extreme mutation bias and high AT content in Plasmodium falciparum

For reasons that remain unknown, the Plasmodium falciparum genome has an exceptionally high AT content compared to other Plasmodium species and eukaryotes in general - nearly 80% in coding regions and approaching 90% in non-coding regions. Here, we examine how this phenomenon relates to genome-wide patterns of de novo mutation. Mutation accumulation experiments were performed by sequential cloning of six P. falciparum isolates growing in human erythrocytes in vitro for 4 years, with 279 clones sampled for whole genome sequencing at different time points. Genome sequence analysis of these samples revealed a significant excess of G:C to A:T transitions compared to other types of nucleotide substitution, which would naturally cause AT content to equilibrate close to the level seen across the P. falciparum reference genome (80.6% AT). These data also uncover an extremely high rate of small indel mutation relative to other species, primarily associated with repetitive AT-rich sequences, in addition to larger-scale structural rearrangements focused in antigen-coding var genes. In conclusion, high AT content in P. falciparum is driven by a systematic mutational bias and ultimately leads to an unusual level of microstructural plasticity, raising the question of whether this contributes to adaptive evolution.

Genomics of apicomplexan parasites

The increasing prevalence of infections involving intracellular apicomplexan parasites such as Plasmodium, Toxoplasma, and Cryptosporidium (the causative agents of malaria, toxoplasmosis, and cryptosporidiosis, respectively) represent a significant global healthcare burden. Despite their significance, few treatments are available; a situation that is likely to deteriorate with the emergence of new resistant strains of parasites. To lay the foundation for programs of drug discovery and vaccine development, genome sequences for many of these organisms have been generated, together with large-scale expression and proteomic datasets. Comparative analyses of these datasets are beginning to identify the molecular innovations supporting both conserved processes mediating fundamental roles in parasite survival and persistence, as well as lineage-specific adaptations associated with divergent life-cycle strategies. The challenge is how best to exploit these data to derive insights into parasite virulence and identify those genes representing the most amenable targets. In this review, we outline genomic datasets currently available for apicomplexans and discuss biological insights that have emerged as a consequence of their analysis. Of particular interest are systems-based resources, focusing on areas of metabolism and host invasion that are opening up opportunities for discovering new therapeutic targets.

Malaria parasites possess a telomere repeat-binding protein that shares ancestry with transcription factor IIIA

Telomere repeat-binding factors (TRFs) are essential components of the molecular machinery that regulates telomere function. TRFs are widely conserved across eukaryotes and bind duplex telomere repeats via a characteristic MYB-type domain. Here, we identified the telomere repeat-binding protein PfTRZ in the malaria parasite Plasmodium falciparum, a member of the Alveolate phylum for which TRFs have not been described so far. PfTRZ lacks an MYB domain and binds telomere repeats via a C₂H₂-type zinc finger domain instead. In vivo, PfTRZ binds with high specificity to the telomeric tract and to interstitial telomere repeats upstream of subtelomeric virulence genes. Conditional depletion experiments revealed that PfTRZ regulates telomere length homeostasis and is required for efficient cell cycle progression. Intriguingly, we found that PfTRZ also binds to and regulates the expression of 5S rDNA genes. Combined with detailed phylogenetic analyses, our findings identified PfTRZ as a remote functional homologue of the basic transcription factor TFIIIA, which acquired a new function in telomere maintenance early in the apicomplexan lineage. Our work sheds unexpected new light on the evolution of telomere repeat-binding proteins and paves the way for dissecting the presumably divergent mechanisms regulating telomere functionality in one of the most deadly human pathogens.

rROP2 from Toxoplasma gondii as a potential vaccine against oocyst shedding in domestic cats

Abstract The aim of the present study was to evaluate oocyst shedding in cats immunized by nasal route with T. gondii proteins ROP2. Twelve short hair cats (Felis catus) were divided in three groups G1, G2 and G3 (n=4). Animals from G1 received 100 μg of rROP2 proteins plus 20 μg of Quil-A, G2 received 100 μg of BSA plus 20 μg of Quil-A, and the G3 only saline solution (control group). All treatments were done by intranasal route at days 0, 21, 42, and 63. The challenge was performed in all groups on day 70 with ≅ 800 tissue cysts of ME-49 strain by oral route. Animals from G1 shed less oocysts (86.7%) than control groups. ELISA was used to detect anti-rROP2 IgG and IgA, however, there were no correlation between number of oocyst shedding by either IgG or IgA antibody levels. In the present work, in spite of lesser oocysts production in immunized group than control groups, it was not possible to associate the use of rROP2 via nostrils with protection against oocyst shedding. For the future, the use of either other recombinant proteins or DNA vaccine, in combination with rROP2 could be tested to try improving the efficacy of this kind of vaccine.

Evolution of the merozoite surface protein 7 (msp7) family in Plasmodium vivax and P. falciparum: A comparative approach

Malaria parasites (genus Plasmodium) are a diverse group found in many species of vertebrate hosts. These parasites invade red blood cells in a complex process comprising several proteins, many encoded by multigene families, one of which is merozoite surface protein 7 (msp7). In the case of Plasmodium vivax, the most geographically widespread human-infecting species, differences in the number of paralogs within multigene families have been previously explained, at least in part, as potential adaptations to the human host. To explore this in msp7, we studied its orthologs in closely related nonhuman primate parasites; investigating both paralog evolutionary history and genetic polymorphism. The emerging patterns were then compared with the human parasite Plasmodium falciparum. We found that the evolution of the msp7 family is consistent with a birth-and-death model, where duplications, pseudogenizations, and gene loss events are common. However, all paralogs in P. vivax and P. falciparum had orthologs in their closely related species in non-human primates indicating that the ancestors of those paralogs precede the events leading to their origins as human parasites. Thus, the number of paralogs cannot be explained as an adaptation to human hosts. Although there is no functional information for msp7 in P. vivax, we found evidence for purifying selection in the genetic polymorphism of some of its paralogs as well as their orthologs in closely related non-human primate parasites. We also found evidence indicating that a few of P. vivax's paralogs may have diverged from their orthologs in non-human primates by episodic positive selection. Hence, they may had been under selection when the lineage leading to P. vivax diverged from the Asian non-human primates and switched into Homininae. All these lines of evidence suggest that msp7 is functionally important in P. vivax.

Genome microsatellite diversity within the Apicomplexa phylum

The Apicomplexa phylum groups include unicellular and obligate intracellular protozoan parasites with an apical complex used for attachment and invasion to host cells. In this study, we analyze single sequence repeats (SSRs) in the whole genome of 20 apicomplexan organisms that represent four different lineages within the phylum. Only perfect SSRs with at least 12 nucleotides and composed of 2-6 mers were included. To better understand the association of SSR types with the genomic regions, the SSRs were classified accordingly with the genomic location into exon, intron and intergenic categories. Our results showed heterogeneous SSRs density within the studied genomes. However, the most frequent SSRs types were di- and tri-nucleotide repeats. The former was associated with intergenic regions, while the latter was associated with exon regions.

Origin of a major infectious disease in vertebrates: The timing of Cryptosporidium evolution and its hosts

Protozoan parasites of the genus Cryptosporidium infect all vertebrate groups and display some host specificity in their infections. It is therefore possible to assume that Cryptosporidium parasites evolved intimately aside with vertebrate lineages. Here we propose a scenario of Cryptosporidium–Vertebrata coevolution testing the hypothesis that the origin of Cryptosporidium parasites follows that of the origin of modern vertebrates. We use calibrated molecular clocks and cophylogeny analyses to provide and compare age estimates and patterns of association between these clades. Our study provides strong support for the evolution of parasitism of Cryptosporidium with the rise of the vertebrates about 600 million years ago (Mya). Interestingly, periods of increased diversification in Cryptosporidium coincides with diversification of crown mammalian and avian orders after the Cretaceous-Palaeogene (K-Pg) boundary, suggesting that adaptive radiation to new mammalian and avian hosts triggered the diversification of this parasite lineage. Despite evidence for ongoing host shifts we also found significant correlation between protozoan parasites and vertebrate hosts trees in the cophylogenetic analysis. These results help us to understand the underlying macroevolutionary mechanisms driving evolution in Cryptosporidium and may have important implications for the ecology, dynamics and epidemiology of cryptosporidiosis disease in humans and other animals.

Host Mitochondrial Association Evolved in the Human Parasite Toxoplasma gondii via Neofunctionalization of a Gene Duplicate

In Toxoplasma gondii, an intracellular parasite of humans and other animals, host mitochondrial association (HMA) is driven by a gene family that encodes multiple mitochondrial association factor 1 (MAF1) proteins. However, the importance of MAF1 gene duplication in the evolution of HMA is not understood, nor is the impact of HMA on parasite biology. Here we used within- and between-species comparative analysis to determine that the MAF1 locus is duplicated in T. gondii and its nearest extant relative Hammondia hammondi, but not another close relative, Neospora caninum. Using cross-species complementation, we determined that the MAF1 locus harbors multiple distinct paralogs that differ in their ability to mediate HMA, and that only T. gondii and H. hammondi harbor HMA⁺ paralogs. Additionally, we found that exogenous expression of an HMA⁺ paralog in T. gondii strains that do not normally exhibit HMA provides a competitive advantage over their wild-type counterparts during a mouse infection. These data indicate that HMA likely evolved by neofunctionalization of a duplicate MAF1 copy in the common ancestor of T. gondii and H. hammondi, and that the neofunctionalized gene duplicate is selectively advantageous.

Peroxisomes in parasitic protists

Representatives of all major lineages of eukaryotes contain peroxisomes with similar morphology and mode of biogenesis, indicating a monophyletic origin of the organelles within the common ancestor of all eukaryotes. Peroxisomes originated from the endoplasmic reticulum, but despite a common origin and shared morphological features, peroxisomes from different organisms show a remarkable diversity of enzyme content and the metabolic processes present can vary dependent on nutritional or developmental conditions. A common characteristic and probable evolutionary driver for the origin of the organelle is an involvement in lipid metabolism, notably H₂O₂-dependent fatty-acid oxidation. Subsequent evolution of the organelle in different lineages involved multiple acquisitions of metabolic processes-often involving retargeting enzymes from other cell compartments-and losses. Information about peroxisomes in protists is still scarce, but available evidence, including new bioinformatics data reported here, indicate striking diversity amongst free-living and parasitic protists from different phylogenetic supergroups. Peroxisomes in only some protists show major involvement in H₂O₂-dependent metabolism, as in peroxisomes of mammalian, plant and fungal cells. Compartmentalization of glycolytic and gluconeogenic enzymes inside peroxisomes is characteristic of kinetoplastids and diplonemids, where the organelles are hence called glycosomes, whereas several other excavate parasites (Giardia, Trichomonas) have lost peroxisomes. Amongst alveolates and amoebozoans patterns of peroxisome loss are more complicated. Often, a link is apparent between the niches occupied by the parasitic protists, nutrient availability, and the absence of the organelles or their presence with a specific enzymatic content. In trypanosomatids, essentiality of peroxisomes may be considered for use in anti-parasite drug discovery.

Pellicle formation in the malaria parasite

The intraerythrocytic developmental cycle of Plasmodium falciparum is completed with the release of up to 32 invasive daughter cells, the merozoites, into the blood stream. Before release, the final step of merozoite development is the assembly of the cortical pellicle, a multi-layered membrane structure. This unique apicomplexan feature includes the inner membrane complex (IMC) and the parasite's plasma membrane. A dynamic ring structure, referred to as the basal complex, is part of the IMC and helps to divide organelles and abscises in the maturing daughter cells. Here, we analyze the dynamics of the basal complex of P. falciparum. We report on a novel transmembrane protein of the basal complex termed BTP1, which is specific to the genus Plasmodium. It colocalizes with the known basal complex marker protein MORN1 and shows distinct dynamics as well as localization when compared to other IMC proteins during schizogony. Using a parasite plasma membrane marker cell line, we correlate dynamics of the basal complex with the acquisition of the maternal membrane. We show that plasma membrane invagination and IMC propagation are interlinked during the final steps of cell division.

Genomic changes during evolution of animal parasitism in eukaryotes

Understanding how pathogens have evolved to survive in close association with their hosts is an important step in unraveling the biology of host-pathogen interactions. Comparative genomics is a powerful tool to approach this problem as an increasing number of genomes of multiple pathogen species and strains become available. The ever-growing catalog of genome sequences makes comparison of organisms easier, but it also allows us to reconstitute the evolutionary processes occurring at the genomic level that may have led to the acquisition of pathogenic or parasitic mechanisms.

A lactate and formate transporter in the intraerythrocytic malaria parasite, Plasmodium falciparum

The intraerythrocytic malaria parasite relies primarily on glycolysis to fuel its rapid growth and reproduction. The major byproduct of this metabolism, lactic acid, is extruded into the external medium. In this study, we show that the human malaria parasite Plasmodium falciparum expresses at its surface a member of the microbial formate-nitrite transporter family (PfFNT), which, when expressed in Xenopus laevis oocytes, transports both formate and lactate. The transport characteristics of PfFNT in oocytes (pH-dependence, inhibitor-sensitivity and kinetics) are similar to those of the transport of lactate and formate across the plasma membrane of mature asexual-stage P. falciparum trophozoites, consistent with PfFNT playing a major role in the efflux of lactate and hence in the energy metabolism of the intraerythrocytic parasite.

Factors mediating plastid dependency and the origins of parasitism in apicomplexans and their close relatives

Apicomplexans are a major lineage of parasites, including causative agents of malaria and toxoplasmosis. How such highly adapted parasites evolved from free-living ancestors is poorly understood, particularly because they contain nonphotosynthetic plastids with which they have a complex metabolic dependency. Here, we examine the origin of apicomplexan parasitism by resolving the evolutionary distribution of several key characteristics in their closest free-living relatives, photosynthetic chromerids and predatory colpodellids. Using environmental sequence data, we describe the diversity of these apicomplexan-related lineages and select five species that represent this diversity for transcriptome sequencing. Phylogenomic analysis recovered a monophyletic lineage of chromerids and colpodellids as the sister group to apicomplexans, and a complex distribution of retention versus loss for photosynthesis, plastid genomes, and plastid organelles. Reconstructing the evolution of all plastid and cytosolic metabolic pathways related to apicomplexan plastid function revealed an ancient dependency on plastid isoprenoid biosynthesis, predating the divergence of apicomplexan and dinoflagellates. Similarly, plastid genome retention is strongly linked to the retention of two genes in the plastid genome, sufB and clpC, altogether suggesting a relatively simple model for plastid retention and loss. Lastly, we examine the broader distribution of a suite of molecular characteristics previously linked to the origins of apicomplexan parasitism and find that virtually all are present in their free-living relatives. The emergence of parasitism may not be driven by acquisition of novel components, but rather by loss and modification of the existing, conserved traits.

Systems-based analysis of the Sarcocystis neurona genome identifies pathways that contribute to a heteroxenous life cycle

Sarcocystis neurona is a member of the coccidia, a clade of single-celled parasites of medical and veterinary importance including Eimeria, Sarcocystis, Neospora, and Toxoplasma. Unlike Eimeria, a single-host enteric pathogen, Sarcocystis, Neospora, and Toxoplasma are two-host parasites that infect and produce infectious tissue cysts in a wide range of intermediate hosts. As a genus, Sarcocystis is one of the most successful protozoan parasites; all vertebrates, including birds, reptiles, fish, and mammals are hosts to at least one Sarcocystis species. Here we sequenced Sarcocystis neurona, the causal agent of fatal equine protozoal myeloencephalitis. The S. neurona genome is 127 Mbp, more than twice the size of other sequenced coccidian genomes. Comparative analyses identified conservation of the invasion machinery among the coccidia. However, many dense-granule and rhoptry kinase genes, responsible for altering host effector pathways in Toxoplasma and Neospora, are absent from S. neurona. Further, S. neurona has a divergent repertoire of SRS proteins, previously implicated in tissue cyst formation in Toxoplasma. Systems-based analyses identified a series of metabolic innovations, including the ability to exploit alternative sources of energy. Finally, we present an S. neurona model detailing conserved molecular innovations that promote the transition from a purely enteric lifestyle (Eimeria) to a heteroxenous parasite capable of infecting a wide range of intermediate hosts.

Sarcocystis neurona is a member of the coccidia, a clade of single-celled apicomplexan parasites responsible for major economic and health care burdens worldwide. A cousin of Plasmodium, Cryptosporidium, Theileria, and Eimeria, Sarcocystis is one of the most successful parasite genera; it is capable of infecting all vertebrates (fish, reptiles, birds, and mammals—including humans). The past decade has witnessed an increasing number of human outbreaks of clinical significance associated with acute sarcocystosis. Among Sarcocystis species, S. neurona has a wide host range and causes fatal encephalitis in horses, marine mammals, and several other mammals. To provide insights into the transition from a purely enteric parasite (e.g., Eimeria) to one that forms tissue cysts (Toxoplasma), we present the first genome sequence of S. neurona. Comparisons with other coccidian genomes highlight the molecular innovations that drive its distinct life cycle strategies.

The role of palmitoylation for protein recruitment to the inner membrane complex of the malaria parasite

To survive and persist within its human host, the malaria parasite Plasmodium falciparum utilizes a battery of lineage-specific innovations to invade and multiply in human erythrocytes. With central roles in invasion and cytokinesis, the inner membrane complex, a Golgi-derived double membrane structure underlying the plasma membrane of the parasite, represents a unique and unifying structure characteristic to all organisms belonging to a large phylogenetic group called Alveolata. More than 30 structurally and phylogenetically distinct proteins are embedded in the IMC, where a portion of these proteins displays N-terminal acylation motifs. Although N-terminal myristoylation is catalyzed co-translationally within the cytoplasm of the parasite, palmitoylation takes place at membranes and is mediated by palmitoyl acyltransferases (PATs). Here, we identify a PAT (PfDHHC1) that is exclusively localized to the IMC. Systematic phylogenetic analysis of the alveolate PAT family reveals PfDHHC1 to be a member of a highly conserved, apicomplexan-specific clade of PATs. We show that during schizogony this enzyme has an identical distribution like two dual-acylated, IMC-localized proteins (PfISP1 and PfISP3). We used these proteins to probe into specific sequence requirements for IMC-specific membrane recruitment and their interaction with differentially localized PATs of the parasite.

The genome of Eimeria falciformis--reduction and specialization in a single host apicomplexan parasite

BACKGROUND

The phylum Apicomplexa comprises important unicellular human parasites such as Toxoplasma and Plasmodium. Eimeria is the largest and most diverse genus of apicomplexan parasites and some species of the genus are the causative agent of coccidiosis, a disease economically devastating in poultry. We report a complete genome sequence of the mouse parasite Eimeria falciformis. We assembled and annotated the genome sequence to study host-parasite interactions in this understudied genus in a model organism host.

RESULTS

The genome of E. falciformis is 44 Mb in size and contains 5,879 predicted protein coding genes. Comparative analysis of E. falciformis with Toxoplasma gondii shows an emergence and diversification of gene families associated with motility and invasion mainly at the level of the Coccidia. Many rhoptry kinases, among them important virulence factors in T. gondii, are absent from the E. falciformis genome. Surface antigens are divergent between Eimeria species. Comparisons with T. gondii showed differences between genes involved in metabolism, N-glycan and GPI-anchor synthesis. E. falciformis possesses a reduced set of transmembrane transporters and we suggest an altered mode of iron uptake in the genus Eimeria.

CONCLUSIONS

Reduced diversity of genes required for host-parasite interaction and transmembrane transport allow hypotheses on host adaptation and specialization of a single host parasite. The E. falciformis genome sequence sheds light on the evolution of the Coccidia and helps to identify determinants of host-parasite interaction critical for drug and vaccine development.

A survey of innovation through duplication in the reduced genomes of twelve parasites

We characterize the prevalence, distribution, divergence, and putative functions of detectable two-copy paralogs and segmental duplications in the Apicomplexa, a phylum of parasitic protists. Apicomplexans are mostly obligate intracellular parasites responsible for human and animal diseases (e.g. malaria and toxoplasmosis). Gene loss is a major force in the phylum. Genomes are small and protein-encoding gene repertoires are reduced. Despite this genomic streamlining, duplications and gene family amplifications are present. The potential for innovation introduced by duplications is of particular interest. We compared genomes of twelve apicomplexans across four lineages and used orthology and genome cartography to map distributions of duplications against genome architectures. Segmental duplications appear limited to five species. Where present, they correspond to regions enriched for multi-copy and species-specific genes, pointing toward roles in adaptation and innovation. We found a phylum-wide association of duplications with dynamic chromosome regions and syntenic breakpoints. Trends in the distribution of duplicated genes indicate that recent, species-specific duplicates are often tandem while most others have been dispersed by genome rearrangements. These trends show a relationship between genome architecture and gene duplication. Functional analysis reveals: proteases, which are vital to a parasitic lifecycle, to be prominent in putative recent duplications; a pair of paralogous genes in Toxoplasma gondii previously shown to produce the rate-limiting step in dopamine synthesis in mammalian cells, a possible link to the modification of host behavior; and phylum-wide differences in expression and subcellular localization, indicative of modes of divergence. We have uncovered trends in multiple modes of duplicate divergence including sequence, intron content, expression, subcellular localization, and functions of putative recent duplicates that highlight the role of duplications in the continuum of forces that have shaped these genomes.

The origin and diversification of the merozoite surface protein 3 (msp3) multi-gene family in Plasmodium vivax and related parasites

The genus Plasmodium is a diversified group of parasites with more than 200 known species that includes those causing malaria in humans. These parasites use numerous proteins in a complex process that allows them to invade the red blood cells of their vertebrate hosts. Many of those proteins are part of multi-gene families; one of which is the merozoite surface protein-3 (msp3) family. The msp3 multi-gene family is considered important in the two main human parasites, Plasmodium vivax and Plasmodium falciparum, as its paralogs are simultaneously expressed in the blood stage (merozoite) and are immunogenic. There are large differences among Plasmodium species in the number of paralogs in this family. Such differences have been previously explained, in part, as adaptations that allow the different Plasmodium species to invade their hosts. To investigate this, we characterized the array containing msp3 genes among several Plasmodium species, including P. falciparum and P. vivax. We first found no evidence indicating that the msp3 family of P. falciparum was homologous to that of P. vivax. Subsequently, by focusing on the diverse clade of nonhuman primate parasites to which P. vivax is closely related, where homology was evident, we found no evidence indicating that the interspecies variation in the number of paralogs was an adaptation related to changes in host range or host switches. Overall, we hypothesize that the evolution of the msp3 family in P. vivax is consistent with a model of multi-allelic diversifying selection where the paralogs may have functionally redundant roles in terms of increasing antigenic diversity. Thus, we suggest that the expressed MSP3 proteins could serve as "decoys", via antigenic diversity, during the critical process of invading the host red blood cells.

Unique apicomplexan IMC sub-compartment proteins are early markers for apical polarity in the malaria parasite

The phylum Apicomplexa comprises over 5000 intracellular protozoan parasites, including Plasmodium and Toxoplasma, that are clinically important pathogens affecting humans and livestock. Malaria parasites belonging to the genus Plasmodium possess a pellicle comprised of a plasmalemma and inner membrane complex (IMC), which is implicated in parasite motility and invasion. Using live cell imaging and reverse genetics in the rodent malaria model P. berghei, we localise two unique IMC sub-compartment proteins (ISPs) and examine their role in defining apical polarity during zygote (ookinete) development. We show that these proteins localise to the anterior apical end of the parasite where IMC organisation is initiated, and are expressed at all developmental stages, especially those that are invasive. Both ISP proteins are N-myristoylated, phosphorylated and membrane-bound. Gene disruption studies suggest that ISP1 is likely essential for parasite development, whereas ISP3 is not. However, an absence of ISP3 alters the apical localisation of ISP1 in all invasive stages including ookinetes and sporozoites, suggesting a coordinated function for these proteins in the organisation of apical polarity in the parasite.

Variant surface antigens of malaria parasites: functional and evolutionary insights from comparative gene family classification and analysis

Background

Plasmodium parasites, the causative agents of malaria, express many variant antigens on cell surfaces. Variant surface antigens (VSAs) are typically organized into large subtelomeric gene families that play critical roles in virulence and immune evasion. Many important aspects of VSA function and evolution remain obscure, impeding our understanding of virulence mechanisms and vaccine development. To gain further insights into VSA function and evolution, we comparatively classified and examined known VSA gene families across seven Plasmodium species.

Results

We identified a set of ultra-conserved orthologs within the largest Plasmodium gene family pir, which should be considered as high-priority targets for experimental functional characterization and vaccine development. Furthermore, we predict a lipid-binding domain in erythrocyte surface-expressed PYST-A proteins, suggesting a role of this second largest rodent parasite gene family in host cholesterol salvage. Additionally, it was found that PfMC-2TM proteins carry a novel and putative functional domain named MC-TYR, which is conserved in other P. falciparum gene families and rodent parasites. Finally, we present new conclusive evidence that the major Plasmodium VSAs PfEMP1, SICAvar, and SURFIN are evolutionarily linked through a modular and structurally conserved intracellular domain.

Conclusion

Our comparative analysis of Plasmodium VSA gene families revealed important functional and evolutionary insights, which can now serve as starting points for further experimental studies.

Patterns of prokaryotic lateral gene transfers affecting parasitic microbial eukaryotes

Background

The influence of lateral gene transfer on gene origins and biology in eukaryotes is poorly understood compared with those of prokaryotes. A number of independent investigations focusing on specific genes, individual genomes, or specific functional categories from various eukaryotes have indicated that lateral gene transfer does indeed affect eukaryotic genomes. However, the lack of common methodology and criteria in these studies makes it difficult to assess the general importance and influence of lateral gene transfer on eukaryotic genome evolution.

Results

We used a phylogenomic approach to systematically investigate lateral gene transfer affecting the proteomes of thirteen, mainly parasitic, microbial eukaryotes, representing four of the six eukaryotic super-groups. All of the genomes investigated have been significantly affected by prokaryote-to-eukaryote lateral gene transfers, dramatically affecting the enzymes of core pathways, particularly amino acid and sugar metabolism, but also providing new genes of potential adaptive significance in the life of parasites. A broad range of prokaryotic donors is involved in such transfers, but there is clear and significant enrichment for bacterial groups that share the same habitats, including the human microbiota, as the parasites investigated.

Conclusions

Our data show that ecology and lifestyle strongly influence gene origins and opportunities for gene transfer and reveal that, although the outlines of the core eukaryotic metabolism are conserved among lineages, the genes making up those pathways can have very different origins in different eukaryotes. Thus, from the perspective of the effects of lateral gene transfer on individual gene ancestries in different lineages, eukaryotic metabolism appears to be chimeric.

Extensively variable surface antigens of Sarcocystis spp. infecting Brazilian marsupials in the genus Didelphis occur in myriad allelic combinations, suggesting sexual recombination has aided their diversification

Sarcocystis neurona and Sarcocystis falcatula are very similar species of Apicomplexan protozoa that use marsupials of the genus Didelphis as definitive hosts. These mammals can serve as definitive hosts not only for these two parasites, but for other Sarcocystis such as Sarcocystis speeri and Sarcocystis lindsayi. Sarcocystis shed by opossums (with the exception of S. neurona) can cause disease in a great variety of birds, being commonly associated with acute pulmonary sarcocystosis in zoos. S. neurona is the most commonly associated parasite with the equine protozoal myeloencephalitis in horses. Herein we assessed the variability of Sarcocystis spp. isolated from opossums of the state of Rio Grande do Sul, Brazil, by sequencing fragments of genes coding for glycosylphosphatidylinositol-anchored surface antigens (termed surface antigen or SAG), SAG2, SAG3 and SAG4. Two genetic groups were identified, one of them related to S. falcatula and the other related to S. neurona. Various allelic combinations of SAG2, SAG3 and SAG4 occur among S. falcatula related isolates and strong evidences suggest that such isolates may exchange high divergent alleles in possible sexual recombination processes. Regarding the group S. neurona-like (isolates G37 and G38), none of the individuals in this group share alleles with individuals of the other group. Comparing G37 and G38 strains and North American strains of S. neurona, four polymorphisms were identified at SAG-3, five at SAG-2 and three at SAG-4. Gene sequences of locus SAG-3 from isolates G37 and G38 differed from the other sequences by an insertion 81bp long. This insertion contains several dinucleotide repeats of AT, resembling a microsatellite locus and has already been detected in SAG3 sequences of S. neurona from North America. When aligned against North American strains of S. neurona, G37 and G38 isolates have a deletion of 8 nucleotides within this intron which indicate that S. neurona strains of South America are divergent from that of North America. From the results obtained so far, we have shown extensive variability in surface antigens coding sequences among Sarcocystis eliminated by mammals of the genus Didelphis spp. In addition, such divergent alleles may be exchanged in possible sexual recombination processes between different isolates of S. falcatula related isolate. The evolutionary relationships within S. falcatula related isolates will be best clarified after markers less subjected to selection pressures are analyzed in conjunction with surface antigen genes. These results may have a striking impact on the knowledge of the Sarcocystis species that infect opossums in Brazil and also in the epidemiology of the infections caused by these protozoans.

The parasite specific substitution matrices improve the annotation of apicomplexan proteins

Background

A number of apicomplexan genomes have been sequenced successfully in recent years and this would help in understanding the biology of apicomplexan parasites. The members of the phylum Apicomplexa are important protozoan parasites (Plasmodium, Toxoplasma and Cryptosporidium etc) that cause some of the deadly diseases in humans and animals. In our earlier studies, we have shown that the standard BLOSUM matrices are not suitable for compositionally biased apicomplexan proteins. So we developed a novel series (SMAT and PfFSmat60) of substitution matrices which performed better in comparison to standard BLOSUM matrices and developed ApicoAlign, a sequence search and alignment tool for apicomplexan proteins. In this study, we demonstrate the higher specificity of these matrices and make an attempt to improve the annotation of apicomplexan kinases and proteases.

Results

The ROC curves proved that SMAT80 performs best for apicomplexan proteins followed by compositionally adjusted BLOSUM62 (PSI-BLAST searches), BLOSUM90 and BLOSUM62 matrices in terms of detecting true positives. The poor E-values and/or bit scores given by SMAT80 matrix for the experimentally identified coccidia-specific oocyst wall proteins against hematozoan (non-coccidian) parasites further supported the higher specificity of the same. SMAT80 uniquely detected (missed by BLOSUM) orthologs for 1374 apicomplexan hypothetical proteins against SwissProt database and predicted 70 kinases and 17 proteases. Further analysis confirmed the conservation of functional residues of kinase domain in one of the SMAT80 detected kinases. Similarly, one of the SMAT80 detected proteases was predicted to be a rhomboid protease.

Conclusions

The parasite specific substitution matrices have higher specificity for apicomplexan proteins and are helpful in detecting the orthologs missed by BLOSUM matrices and thereby improve the annotation of apicomplexan proteins which are hypothetical or with unknown function.

A transcriptomic analysis of Echinococcus granulosus larval stages: implications for parasite biology and host adaptation

Background

The cestode Echinococcus granulosus - the agent of cystic echinococcosis, a zoonosis affecting humans and domestic animals worldwide - is an excellent model for the study of host-parasite cross-talk that interfaces with two mammalian hosts. To develop the molecular analysis of these interactions, we carried out an EST survey of E. granulosus larval stages. We report the salient features of this study with a focus on genes reflecting physiological adaptations of different parasite stages.

Methodology/Principal Findings

We generated ∼10,000 ESTs from two sets of full-length enriched libraries (derived from oligo-capped and trans-spliced cDNAs) prepared with three parasite materials: hydatid cyst wall, larval worms (protoscoleces), and pepsin/H⁺-activated protoscoleces. The ESTs were clustered into 2700 distinct gene products. In the context of the biology of E. granulosus, our analyses reveal: (i) a diverse group of abundant long non-protein coding transcripts showing homology to a middle repetitive element (EgBRep) that could either be active molecular species or represent precursors of small RNAs (like piRNAs); (ii) an up-regulation of fermentative pathways in the tissue of the cyst wall; (iii) highly expressed thiol- and selenol-dependent antioxidant enzyme targets of thioredoxin glutathione reductase, the functional hub of redox metabolism in parasitic flatworms; (iv) candidate apomucins for the external layer of the tissue-dwelling hydatid cyst, a mucin-rich structure that is critical for survival in the intermediate host; (v) a set of tetraspanins, a protein family that appears to have expanded in the cestode lineage; and (vi) a set of platyhelminth-specific gene products that may offer targets for novel pan-platyhelminth drug development.

Conclusions/Significance

This survey has greatly increased the quality and the quantity of the molecular information on E. granulosus and constitutes a valuable resource for gene prediction on the parasite genome and for further genomic and proteomic analyses focused on cestodes and platyhelminths.

Cestodes are a neglected group of platyhelminth parasites, despite causing chronic infections to humans and domestic animals worldwide. We used Echinococcus granulosus as a model to study the molecular basis of the host-parasite cross-talk during cestode infections. For this purpose, we carried out a survey of the genes expressed by parasite larval stages interfacing with definitive and intermediate hosts. Sequencing from several high quality cDNA libraries provided numerous insights into the expression of genes involved in important aspects of E. granulosus biology, e.g. its metabolism (energy production and antioxidant defences) and the synthesis of key parasite structures (notably, the one exposed to humans and livestock intermediate hosts). Our results also uncovered the existence of an intriguing set of abundant repeat-associated non-protein coding transcripts that may participate in the regulation of gene expression in all surveyed stages. The dataset now generated constitutes a valuable resource for gene prediction on the parasite genome and for further genomic and proteomic studies focused on cestodes and platyhelminths. In particular, the detailed characterization of a range of newly discovered genes will contribute to a better understanding of the biology of cestode infections and, therefore, to the development of products allowing their efficient control.

Integrated bioinformatic and targeted deletion analyses of the SRS gene superfamily identify SRS29C as a negative regulator of Toxoplasma virulence

The Toxoplasma gondii SRS gene superfamily is structurally related to SRS29B (formerly SAG1), a surface adhesin that binds host cells and stimulates host immunity. Comparative genomic analyses of three Toxoplasma strains identified 182 SRS genes distributed across 14 chromosomes at 57 genomic loci. Eight distinct SRS subfamilies were resolved. A core 69 functional gene orthologs were identified, and strain-specific expansions and pseudogenization were common. Gene expression profiling demonstrated differential expression of SRS genes in a developmental-stage- and strain-specific fashion and identified nine SRS genes as priority targets for gene deletion among the tissue-encysting coccidia. A Δsag1 sag2A mutant was significantly attenuated in murine acute virulence and showed upregulated SRS29C (formerly SRS2) expression. Transgenic overexpression of SRS29C in the virulent RH parent was similarly attenuated. Together, these findings reveal SRS29C to be an important regulator of acute virulence in mice and demonstrate the power of integrated genomic analysis to guide experimental investigations.

IMPORTANCE

Parasitic species employ large gene families to subvert host immunity to enable pathogen colonization and cause disease. Toxoplasma gondii contains a large surface coat gene superfamily that encodes adhesins and virulence factors that facilitate infection in susceptible hosts. We generated an integrated bioinformatic resource to predict which genes from within this 182-gene superfamily of adhesin-encoding genes play an essential role in the host-pathogen interaction. Targeted gene deletion experiments with predicted candidate surface antigens identified SRS29C as an important negative regulator of acute virulence in murine models of Toxoplasma infection. Our integrated computational and experimental approach provides a comprehensive framework, or road map, for the assembly and discovery of additional key pathogenesis genes contained within other large surface coat gene superfamilies from a broad array of eukaryotic pathogens.

Comparative genomics of the apicomplexan parasites Toxoplasma gondii and Neospora caninum: Coccidia differing in host range and transmission strategy

Toxoplasma gondii is a zoonotic protozoan parasite which infects nearly one third of the human population and is found in an extraordinary range of vertebrate hosts. Its epidemiology depends heavily on horizontal transmission, especially between rodents and its definitive host, the cat. Neospora caninum is a recently discovered close relative of Toxoplasma, whose definitive host is the dog. Both species are tissue-dwelling Coccidia and members of the phylum Apicomplexa; they share many common features, but Neospora neither infects humans nor shares the same wide host range as Toxoplasma, rather it shows a striking preference for highly efficient vertical transmission in cattle. These species therefore provide a remarkable opportunity to investigate mechanisms of host restriction, transmission strategies, virulence and zoonotic potential. We sequenced the genome of N. caninum and transcriptomes of the invasive stage of both species, undertaking an extensive comparative genomics and transcriptomics analysis. We estimate that these organisms diverged from their common ancestor around 28 million years ago and find that both genomes and gene expression are remarkably conserved. However, in N. caninum we identified an unexpected expansion of surface antigen gene families and the divergence of secreted virulence factors, including rhoptry kinases. Specifically we show that the rhoptry kinase ROP18 is pseudogenised in N. caninum and that, as a possible consequence, Neospora is unable to phosphorylate host immunity-related GTPases, as Toxoplasma does. This defense strategy is thought to be key to virulence in Toxoplasma. We conclude that the ecological niches occupied by these species are influenced by a relatively small number of gene products which operate at the host-parasite interface and that the dominance of vertical transmission in N. caninum may be associated with the evolution of reduced virulence in this species.

Evolution and architecture of the inner membrane complex in asexual and sexual stages of the malaria parasite

The inner membrane complex (IMC) is a unifying morphological feature of all alveolate organisms. It consists of flattened vesicles underlying the plasma membrane and is interconnected with the cytoskeleton. Depending on the ecological niche of the organisms, the function of the IMC ranges from a fundamental role as reinforcement system to more specialized roles in motility and cytokinesis. In this article, we present a comprehensive evolutionary analysis of IMC components, which exemplifies the adaptive nature of the IMCs' protein composition. Focusing on eight structurally distinct proteins in the most prominent "genus" of the Alveolata-the malaria parasite Plasmodium-we demonstrate that the level of conservation is reflected in phenotypic characteristics, accentuated in differential spatial-temporal patterns of these proteins in the motile stages of the parasite's life cycle. Colocalization studies with the centromere and the spindle apparatus reveal their discriminative biogenesis. We also reveal that the IMC is an essential structural compartment for the development of the sexual stages of Plasmodium, as it seems to drive the morphological changes of the parasite during the long and multistaged process of sexual differentiation. We further found a Plasmodium-specific IMC membrane matrix protein that highlights transversal structures in gametocytes, which could represent a genus-specific structural innovation required by Plasmodium. We conclude that the IMC has an additional role during sexual development supporting morphogenesis of the cell, which in addition to its functions in the asexual stages highlights the multifunctional nature of the IMC in the Plasmodium life cycle.

Mechanisms and evolution of virulence in oomycetes

Many destructive diseases of plants and animals are caused by oomycetes, a group of eukaryotic pathogens important to agricultural, ornamental, and natural ecosystems. Understanding the mechanisms underlying oomycete virulence and the genomic processes by which those mechanisms rapidly evolve is essential to developing effective long-term control measures for oomycete diseases. Several common mechanisms underlying oomycete virulence, including protein toxins and cell-entering effectors, have emerged from comparing oomycetes with different genome characteristics, parasitic lifestyles, and host ranges. Oomycete genomes display a strongly bipartite organization in which conserved housekeeping genes are concentrated in syntenic gene-rich blocks, whereas virulence genes are dispersed into highly dynamic, repeat-rich regions. There is also evidence that key virulence genes have been acquired by horizontal transfer from other eukaryotic and prokaryotic species.

Multiple independent fusions of glucose-6-phosphate dehydrogenase with enzymes in the pentose phosphate pathway

Fusions of the first two enzymes in the pentose phosphate pathway, glucose-6-phosphate dehydrogenase (G6PD) and 6-phosphogluconolactonase (6PGL), have been previously described in two distant clades, chordates and species of the malarial parasite Plasmodium. We have analyzed genome and expressed sequence data from a variety of organisms to identify the origins of these gene fusion events. Based on the orientation of the domains and range of species in which homologs can be found, the fusions appear to have occurred independently, near the base of the metazoan and apicomplexan lineages. Only one of the two metazoan paralogs of G6PD is fused, showing that the fusion occurred after a duplication event, which we have traced back to an ancestor of choanoflagellates and metazoans. The Plasmodium genes are known to contain a functionally important insertion that is not seen in the other apicomplexan fusions, highlighting this as a unique characteristic of this group. Surprisingly, our search revealed two additional fusion events, one that combined 6PGL and G6PD in an ancestor of the protozoan parasites Trichomonas and Giardia, and another fusing G6PD with phosphogluconate dehydrogenase (6PGD) in a species of diatoms. This study extends the range of species known to contain fusions in the pentose phosphate pathway to many new medically and economically important organisms.

BEL/Pao retrotransposons in metazoan genomes

BACKGROUND

Long terminal repeat (LTR) retrotransposons are a widespread kind of transposable element present in eukaryotic genomes. They are a major factor in genome evolution due to their ability to create large scale mutations and genome rearrangements. Compared to other transposable elements, little attention has been paid to elements belonging to the metazoan BEL/Pao subclass of LTR retrotransposons. No comprehensive characterization of these elements is available so far. The aim of this study was to describe all BEL/Pao elements in a set of 62 sequenced metazoan genomes, and to analyze their phylogenetic relationship.

RESULTS

We identified a total of 7,861 BEL/Pao elements in 53 of our 62 study genomes. We identified BEL/Pao elements in 20 genomes where such elements had not been found so far. Our analysis shows that BEL/Pao elements are the second-most abundant class of LTR retrotransposons in the genomes we study, more abundant than Ty1/Copia elements, and second only to Ty3/Gypsy elements. They occur in multiple phyla, including basal metazoan phyla, suggesting that BEL/Pao elements arose early in animal evolution. We confirm findings from previous studies that BEL/Pao elements do not occur in mammals. The elements we found can be grouped into more than 1725 families, 1623 of which are new, previously unknown families. These families fall into seven superfamilies, only five of which have been characterized so far. One new superfamily is a major subdivision of the Pao superfamily which we propose to call Dan, because it is restricted to the genome of the zebrafish Danio rerio. The other new superfamily comprises 83 elements and is restricted to lower aquatic eumetazoans. We propose to call this superfamily Flow. BEL/Pao elements do not show any signs of recent horizontal gene transfer between distantly related species.

CONCLUSIONS

In sum, our analysis identifies thousands of new BEL/Pao elements and provides new insights into their distribution, abundance, and evolution.

Parasite annexins--new molecules with potential for drug and vaccine development

In the last few years, annexins have been discovered in several nematodes and other parasites, and distinct differences between the parasite annexins and those of the hosts make them potentially attractive targets for anti-parasite therapeutics. Annexins are ubiquitous proteins found in almost all organisms across all kingdoms.Here, we present an overview of novel annexins from parasitic organisms, and summarize their phylogenetic and biochemical properties, with a view to using them as drug or vaccine targets. Building on structural and biological information that has been accumulated for mammalian and plant annexins, we describe a predicted additional secondary structure element found in many parasite annexins that may confer unique functional properties, and present a specific antigenic epitope for use as a vaccine.

The evolutionary landscape of the chromatin modification machinery reveals lineage specific gains, expansions, and losses

Model organisms such as yeast, fly, and worm have played a defining role in the study of many biological systems. A significant challenge remains in translating this information to humans. Of critical importance is the ability to differentiate those components where knowledge of function and interactions may be reliably inferred from those that represent lineage-specific innovations. To address this challenge, we use chromatin modification (CM) as a model system for exploring the evolutionary properties of their components in the context of their known functions and interactions. Collating previously identified components of CM from yeast, worm, fly, and human, we identified a "core" set of 50 CM genes displaying consistent orthologous relationships that likely retain their interactions and functions across taxa. In addition, we catalog many components that demonstrate lineage specific expansions and losses, highlighting much duplication within vertebrates that may reflect an expanded repertoire of regulatory mechanisms. Placed in the context of a high-quality protein-protein interaction network, we find, contrary to existing views of evolutionary modularity, that CM complex components display a mosaic of evolutionary histories: a core set of highly conserved genes, together with sets displaying lineage specific innovations. Although focused on CM, this study provides a template for differentiating those genes which are likely to retain their functions and interactions across species. As such, in addition to informing on the evolution of CM as a system, this study provides a set of comparative genomic approaches that can be generally applied to any biological systems.

Armadillo-repeat protein functions: questions for little creatures

Armadillo (ARM)-repeat proteins form a large family with diverse and fundamental functions in many eukaryotes. ARM-repeat proteins have largely been characterised in multicellular organisms and much is known about how a subset of these proteins function. The structure of ARM-repeats allows proteins containing them to be functionally very versatile. Are the ARM-repeat proteins in 'little creatures' as multifunctional as their better-studied relatives? The time is now right to start analysing ARM-repeat proteins in these new systems to better understand their cell biology. Here, we review recent advances in understanding the many cellular roles of both well-known and novel ARM-repeat proteins.

An infectious topic in reticulate evolution: introgression and hybridization in animal parasites

Little attention has been given to the role that introgression and hybridization have played in the evolution of parasites. Most studies are host-centric and ask if the hybrid of a free-living species is more or less susceptible to parasite infection. Here we focus on what is known about how introgression and hybridization have influenced the evolution of protozoan and helminth parasites of animals. There are reports of genome or gene introgression from distantly related taxa into apicomplexans and filarial nematodes. Most common are genetic based reports of potential hybridization among congeneric taxa, but in several cases, more work is needed to definitively conclude current hybridization. In the medically important Trypanosoma it is clear that some clonal lineages are the product of past hybridization events. Similarly, strong evidence exists for current hybridization in human helminths such as Schistosoma and Ascaris. There remain topics that warrant further examination such as the potential hybrid origin of polyploid platyhelminths. Furthermore, little work has investigated the phenotype or fitness, and even less the epidemiological significance of hybrid parasites.