- Research article
- Open Access
The genome of Eimeria falciformis - reduction and specialization in a single host apicomplexan parasite
BMC Genomics volume 15, Article number: 696 (2014)
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.
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.
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.
The taxon Eimeria is the largest genus of the phylum Apicomplexa with more than 1,800 species . The Apicomplexa are obligate intracellular parasites and the phylum contains many well-known pathogens of humans and livestock. Plasmodium species causing malaria can be regarded as the most threatening eukaryotes to human . The causative agent of toxoplasmosis, Toxoplasma gondii, infects about a third of the worldwide human population and causes pathology mostly in immunodeficient individuals and neonates .
Eimeriids dwell within cells of the intestine and intestine-related tissues of vertebrates and invertebrates. The family comprises human pathogenic parasites like Isospora belli and Cyclospora caytanensis. The genus Eimeria does not infect humans, but is best known for species infecting domestic animals. Eimeria species cause > 2 billion US$ damage p. a. alone in the poultry industry . Eimeriid parasites have single host life cycles and typically display a strict host and tissue specificity. For example, the chicken is parasitized by 7 species of Eimeria, each restricted to a special part of the intestine and causing severe, self-limiting infections .
The development of E. falciformis Eimer, 1870  is restricted to crypt epithelial cells of the cecum and proximal colon of the mouse [6, 7]. E. falciformis combines the typical life cycle elements of parasites of the subclass Coccidia (Figure 1): Infection occurs by ingestion of oocysts, containing eight haploid sporozoites. Sporozoites infect epithelial cells and undergo several rounds of asexual replication (“schizogony”), leading to very high numbers of progeny called merozoites. These differentiate into gametes (“gamogony”) that fuse to form diploid zygotes. After leaving the host via the feces in environmentally resistant oocysts they undergo a reduction division and mitotic divisions (“sporogony”) to yield the infective sporozoites. As its host, the mouse, is among the best-studied systems in biological research, the E. falciformis infection lends itself as a model to dissect facets of Eimeria-host interaction [8–11].
Plasmodium is a Haemosporidian only distantly related to Eimeria and the Coccidia. T. gondii, the only member of the genus Toxoplasma, belongs to the same subclass as the Eimeriidae, the Coccidia. Its life cycle shows similarities to that of Eimeria in its feline definitive hosts, but in contrast to the specialist, it can infect >180 species of warm blooded animals as intermediate hosts, including humans . Due to this flexibility, this relatively easy to handle parasite has become a well-studied model organism that has greatly furthered the understanding of apicomplexan biology. It has been argued that T. gondii needs its relatively large genome to adapt to its wide variety of intermediate hosts, requiring for example flexibility with regard to host metabolism or immune reposes. Interestingly, T. gondii has diverged into strains with different pathogenicity, supposedly reflecting adaptation to the immune responses of its intermediate hosts . The factors, which led to the divergence of Eimeriids into specialist species and to the coherence of T. gondii as one generalist species, are still not understood.
Previously, high-quality genome sequences have been reported for T. gondii and its very close relatives in the family Sarcocystidae Neospora caninum and Hammondia hammondii. In the family Eimeriidae, first studies of E. tenella chromosome I  and the description of a database for the genome assembly of Eimeria maxima were complemented by a large scale investigation of all 7 species of Eimeria infecting chickens recently . Genome data has not been reported for any Eimeria species infecting a host other than chicken.
We determined a draft genome sequence of the mouse parasite E. falciformis using a next generation sequencing approach, to promote studies on the biology of Eimeriids. We performed comparative genomic analyses and focused on the congener E. tenella and two more distant relatives, T. gondii and N. caninum, to infer the evolutionary history of genes playing roles in core processes in these parasites.
One such core process, the substrate dependent locomotion of coccidian and haemosporidian parasites is called gliding motility. An actin-dependent motor drives this process . This movement is the basis for the invasion process, during which invasive stages of apicomplexan and especially coccidian parasites establish themselves in a parasitophorous vacuole. The secretion of surface molecules from specialized apical organelles called micronemes is triggered by cyclic nucleotide and calcium signalling and transduced by calcium dependent kinases. Other apical organelles, the rhoptries secrete kinases (Rhoptry kinases; RopKs) and pseudokinases interacting with the host immune system .
Comparative genomics of the complete E. falciformis genome allow a comprehensive view of the shared gene repertoire among apicomplexan and coccidian parasites. Our data for a rodent Eimeriid species allows identifying genus specific processes and molecules. The presented genome sequence thus provides a better understanding of the biology of Eimeriids and a basis to establish E. falciformis as a suitable model-parasite.
Results and discussion
Genome sequence, size and completeness
We determined the draft genome sequence of the mouse parasite Eimeria falciformis using a next generation sequencing approach comprising one standard paired-end and one mate pair library with ~ 2 Kilobase (Kb) insert size sequenced with Illumina technology. A total nominal sequencing coverage of more than 500× allowed the construction of a largely contiguous whole genome sequence. We assembled 753 high quality contigs covering a sequence space of 43.7 Mb with a median weighted contig size (N50) of ~250 kb and thus conclude a genome size of roughly 44 Mb. The five longest contigs of the assembly cover a combined length of 4.6 Mb and could very well represent some of the typically 14 chromosomes found in Eimeria species and other Coccidia . With this size the genome of E. falciformis is smaller than those of other Eimeria species . It ranges between the size of the strongly reduced genomes of Cryptosporidium parvum and the relatively large one of T. gondii (Table 1). We have submitted the E. falciformis data to the Short Read Archive as project SRP034650, the genome assembly and relevant annotations will also be accessible through ToxoDB .
Most apicomplexan parasites harbor a mitochondrion and a rudimentary plastid organelle, the apicoplast, both of which contain extra-nuclear DNA. We estimate the size of both organellar genomes to 6.2 kb and ~33 kb, respectively. Both genomes are present in higher copy numbers relative to the nuclear genome as estimated from sequencing coverage (see Additional file 1). We predict ~180 copies per cell for the mitochondrial genome and 18–19 copies per cell for the apicoplast genome.
The mitochondrial genome of E. falciformis is in good agreement with GC content (34.5%) and size (6.2 kb) estimates for mitochondrial genomes of avian Eimeria species . For the protein encoding genes of the mitochondrial genome cytochrome c oxidase subunits 1 (cox1), 3 (cox3) and cytochrome oxidase b (cytb) expression evidence was found in RNA-Seq data.
Our attempts to assemble a single contiguous apicoplast genome sequence did not succeed with the available data. By comparison to the E. tenella plastid genome  we inferred the repeat structure of the ribosomal genes (two inverted repeats) as the likely cause of the assembly problems. Hence we reconstructed the E. falciformis apicoplast genome by assuming conserved synteny with the E. tenella genome (Additional file 2). The 33 kb apicoplast genome is AT rich with a GC content of just ~23%. Apicoplast encoded genes were not represented in our RNA-Seq data. It remains to be tested whether this is due to a lack of transcript polyadenylation in this organelle. In dinoflagellates polyadenylation has been used as an indication for nuclear origin (instead of plastid origin) before .
The overall GC content of the nuclear genome (52.9%) is similar to the E. tenella genome (51.3%). We found trimer-repeats of CAG/GTC (repeated 9–16 times) and heptamer-repeats of AAACCCT/AGGGTTT (repeated 6–14 times) as the most abundant simple sequence repeat (SSRs) instances. Another considerable proportion of the genome (9.6 Mb; 22%) consists of more complex non-SSR elements. Both SSR and complex repeats are spread in a pattern consistent with the segmental organization of the E. tenella and other chicken Eimeria genomes [16, 18], i.e. they occur together and are not confined to few telomeric regions.
Trimer repeats (and other k-mer repeats, where k is a multiple of 3) are often (58% of cases) found in protein coding sequences, heptamer and other “frame-incompatible” repeats are found outside of coding sequence in 93% of their occurrences.
The specific base combination CAG (repeated at least 7 times) is found in 1815 genes, leading to conceptual translations into homopolymeric amino acid repeats (HAARs) of alanine (A), glutamine (Q) or serine (S) in 1781 (30%) proteins. This is a smaller proportion compared to the 57% of genes containing HAARs in E. tenella. Nevertheless the identical conservation of nucleotide sequence coding for HAARs but also non-coding heptamer repeats outside of protein coding genes is remarkable. It suggests a non-protein-mediated functional role e.g. in nuclear organization or as a second possibility the conservation of the generating mechanism and a near neutrality of the SSRs and the resulting HAARs.
The nuclear genome contains 151 functional tRNA genes for 46 codons of 20 standard proteinogenic amino acids. These tRNAs genes were found on 84 different contigs, with no contig containing more than 10 tRNAs genes. Clusters of multiple rRNA genes for the nuclear large and small ribosomal subunits (LSU and SSU) were found on 3 contigs. Additional single copies of LSU and SSU were found on 6 and 2 contigs, respectively. SSU displayed polymorphism (both SNPs and indels) in alignable regions of different rRNA genes.
Combining expression evidence from RNA-seq data with ab initio predictions we inferred 5,879 protein-coding genes producing 6,586 transcripts and predicted conceptual translations for proteins. We estimate that 28.2 Mb (~64.5%) of the E. falciformis genome is contained within protein coding gene loci with 14.6 Mb (~33%) coding sequence. Coding sequence exons typically show a higher GC (55.9%) content than introns (47.1%) (Additional file 3).
Protein coding genes of E. falciformis have a median span of 3.3 kb (Table 1). As a result of both a higher number per gene and longer introns (Additional file 4) Coccidia have longer genes then e.g. the Haemosporidia. The genomes of Plasmodium species contain mainly small sized exons, but a distinct population of larger exons is also found. These are comparable in size to the single exon genes found in Cryptosporidium species. Partial or complete intron loss, in these genomes, resulting in merged exons has been ascribed to bursts of transposon activity . No similar signs of intron loss as relics of past or present transposon activity are found in protein coding genes of E. falciformis.
The location of regions with similarity to protein coding genes of transposons (TransposonPSI), however, was found to correlate with repetitive regions. This might indicate a role of these (retro-)transposon elements in the acquisition and spread of repeats. We thus suggest that transposons have not had a large role in shaping the protein coding genes of Coccidia and especially E. falciformis but potentially had a role in the emergence of the non-coding repeat repertoire. Detailed organization, categorization and origin of these repeats constitutes an interesting area of research for comparative genomics within the genus Eimeria as more complete genome sequences become available.
As observed by Reid et al.  repeats are not found within the functional domains of protein coding genes. We obtained domain annotations using InterProScan  for 3,587 (54.4%) of our predicted proteins and found 23,059 domain occurrences, representing 2,940 distinct InterPro accessions (Additional file 5). Based on these annotations Gene Ontology (GO)  terms could be obtained for 2,841 genes (Additional file 6). For a direct comparison of proteomes, we also predicted 26,584 domains in 4,498 genes covering 3,321 distinct accessions in T. gondii (Additional file 7). Comparison of the Interpro domain annotations of E. falciformis with those obtained for the T. gondii proteome and with protein domain annotations from all vertebrate, bacterial and archaeal proteomes as listed in SwissProt (Figure 2) yielded species-specific sets of domains.
Comparative functional analyses of coccidian proteomes
It is conceivable that the diversity of encoded proteins should vary with the number and complexity of the niches a pathogen meets in its host(s). If this held true, Eimeria parasites dwelling within relatively specialized cells like enterocytes or endothelial cells of a single host can be expected to have a less complex protein repertoire, as compared to parasites exploiting several very different niches in two hosts, like T. gondii. To elaborate specific biological functions of Eimeriids, we clustered 118,629 proteins from apicomplexan genomes plus proteins from Arabidopsis thaliana, Saccharomyces cerevisiae and the diatom Thalassiosira pseudonana as outgroup. We obtained 76,740 ortholog clusters, containing 4,466 E. falciformis genes in 4,205 clusters of which 4,186 contained additional genes from other species (Table 2). This analysis also demonstrates that our E. falciformis gene-set is complete for conserved genes of apicomplexan parasites. A comparison of the size of gene sets of E. falciformis and E. tenella (5,879 vs. 9,262 predicted genes) hence suggests that the smaller gene-set might be less redundant.
To obtain information on the occurrence of biological specializations, we assigned ortholog clusters into categories based on their phylogenetic profile (Figure 3). Our analysis identified a high number of genus Eimeria and Coccidia-specific gene clusters, which is in agreement with previous findings . These novel or diverged gene-families are expected to contain many of the genes responsible for the diverse biological specializations related to differences in modes of parasitism (e.g. host-usage) and life cycle (incl. reproduction). We thus profiled the clusters extensively for enrichment of particular functions (Additional file 8).
To further infer the evolutionary history of these gene clusters, we constructed phylogenetic trees for all ortholog clusters and compared these with the species phylogeny and identified gene families that underwent noticeable expansion events. We detected 22 gene clusters expanded at the species level in E. falciformis. In 135 ortholog clusters expansions at the level of the genus Eimeria was inferred (with 235 E. falciformis genes). Moreover, 135 ortholog clusters (199 E. falciformis genes) showed Coccidian specific expansions and 456 ortholog clusters (553 E. falciformis genes) expanded at the “crown group Apicomplexa” level. Within the latter ortholog clusters, copies of E. falciformis genes were obviously often lost. Additional file 9 highlights functions overabundant in clusters expanded in particular clades of the Apicomplexa.
Some ortholog clusters showed continued expansion throughout their evolutionary history: A cluster of myosin head domain containing proteins for example started expanding early in the evolution of the Apicomplexa and continued to expand to the genus level in Eimeria. This occurrence of a large cluster of genes with motor activity and myosin heads expanding throughout the phylogeny is accompanied by diversification into novel or highly diverged ortholog clusters with the same annotation especially at the Coccidia level. This can be interpreted as a continued evolutionary invention of slightly altered functions for these motility associated genes.
The phylogenetic distance between E. falciformis paralogs in ortholog clusters expanded at the class or genus level is larger compared to the most recent expansions at the species level. These most recent E. falciformis paralogs are likely less diverged simply due to a shorter time allowed for divergence. More intriguingly paralogs with a (sub-)class or genus origin have also a higher average distance than those arising basally in the Apicomplexa (Additional file 10). The function of these older paralogs gives a hint on possible reasons. Genes in clusters expanded basally in the Apicomplexa are involved in conserved processes like transcription, translation and metabolism. The acquisition of related core functions of novel paralogs could explain their low divergence. In contrast neofunctionalization of novel paralogs during the cladogenesis of the Coccidia and the genus Eimeria, accompanied by the acquisition of new niches, could have produced more divergent proteins.
Interdependent gene family expansions in Eimeriidae and Sarcocystidae
Reconciliation of gene trees with the species tree allowed us to highlight cases of independent gene family expansions in both branches of the Coccidia, namely Eimeriidae and Sarcocystidae. It can be hypothesized that in these cases functional niches similar in the Sarcocystidae and Eimeriidae could be filled by the new paralogs. These novel functions would in turn be highly interesting because of their repeated evolutionary invention.
We noticed such independently expanded gene families in total in 16 ortholog clusters, containing 26 E. falciformis genes. In ten of these families the duplicated paralog was only present in E. tenella and the remaining six clusters comprised 16 E. falciformis genes.
Two E. falciformis paralogs in one of these clusters were annotated as fructose-bisphosphate aldolase (FBA; Additional file 11). FBA typically functions as a metabolic enzyme in glycolysis and energy production. However, one of the T. gondii FBA paralogs is also essential for host cell invasion linking a micronemal protein to the cytoskeleton. This mechanism potentially allows the spatially focused supply of energy during invasion . The independently arisen expansion of Eimeria paralogs poses the question whether Eimeria has developed new functions relating glycolysis to cell invasion in a parallel evolutionary process.
As a second illustrative example, membrane alanyl aminopeptidase showed independent expansions into multiple paralogs in Eimeria and Sarcocystidae (Additional file 12). In this case the complete set of paralog copies was only retained in the genomes of E. falciformis and T. gondii (3 and 4 respectively). In Plasmodium, the single orthologous protein is involved in hemoglobin digestion and is a possible drug target due to its indispensability to the parasite . The function of these peptidases and their independent paralogous expansions and retention (and hence assumed neofunctionalisation) in some Coccidia awaits further investigation.
Surface proteins are of particular interest for studies on pathogen host interactions as they are exposed on the outer membranes of parasites and are thus potential targets of protective immune reactions. Ortholog clusters with a phylogenetic distribution limited to the genus Eimeria were enriched for traits characteristic of surface proteins, namely sequences with secretion signal (SignalP) likely to be exported by the secretory pathway. Such signal-peptide carrying E. falciformis proteins were enriched for “serine-type endopeptidase activity” (GO:0004252). Coccidia-specific clusters on the other hand were enriched for sequences with other traits of surface proteins: signal anchors (SignalP) and transmembrane helices (TMHMM). Such transmembrane helices were also overrepresented in E. falciformis species specific clusters.
The domain “Pollen allergen Poa pIX/Phl pV” (IPR001778; found in seven genes) was found restricted to E. falciformis. The same domains are found in Trypanosoma cruzi mucin associated surface proteins (MASPs). MASPs are important antigenic peptides and coordinated expression of the MASPs repertoire is likely important for the evasion of host immune responses .
Our analysis also revealed 514 domains only present in T. gondii genes, but not in the E. falciformis gene set. The most prominent example among these was the SRS (SAG related; IPR007226) domain in 107 T. gondii genes. The SRS proteins are surface proteins putatively involved in binding to host cell receptors and attachment of the asexual replicative forms present in intermediate hosts .
By overlaying our protein domain annotation and OrthoMCL clustering results, we identified two surface antigen (SAG) clusters containing “Sporozoite TA4 surface antigen” domains in large paralog clusters of E. falciformis and E. tenella. The eponymic protein SAG1 induces antibody production in chickens and has been characterized and cloned from E. tenella sporozoites as a prominent vaccine candidate . Other genes containing the same domain were later shown to be glycosylphosphatidylinositol(GPI)-anchored variant surface proteins and to be expressed rather in merozoite stages .
Merging both clusters and computing a single gene tree confirmed two distinct clades of SAGs for rodent and avian Eimeria (Additional file 13). Only one E. falciformis gene was found to be conserved across species and clustered with basal E. tenella surface antigens which are also conserved across the species infecting chicken (SAGa according to ). SAGb and SAGc versions of the domain were not found in E. falciformis.
This conserved SAGa version of the domain could very well have a conserved function in the attachment to cells  and also be a possible cause of cross-protective immune responses induced by some experimental anti-Eimeria vaccines . The function of the SAGa proteins diverged between mammalian and avian Eimeria species awaits further investigation.
The largest species-specific cluster of paralogs in E. falciformis comprised 14 paralogs without any annotation. A single Eimeria-only cluster with 103 hypothetical proteins from E. tenella had notably no member from E. falciformis. While we cannot infer detailed function of these expanded paralog families, similar gene clusters from other species suggest role in host-parasite interaction.
The restriction of surface antigens to a particular species reflects the high specificity of host-parasite interaction within the genus Eimeria and in the Apicomplexa in general. These surface antigens are prime candidates for determinants of the strict preference of the parasites for particular sites within the intestinal tract.
Micronemes, organelles of the apical complex, secrete host interacting molecules early during the invasion of a new host cell. Some of the characteristics (i.e. domain content) of these proteins both in T. gondii and E. tenella are known (reviewed in ). Functional interpretation of our phylogenetic profiling revealed genes likely involved in invasion and egress:
Ortholog clusters restricted to the Coccidia were found heavily enriched for “EF-hand-like” “Epidermal growth factor-like” and “Peptidase S8, subtilisin-related” domains. The latter are found in micronemal proteins of T. gondii cleaving other proteins as part of the secretory pathway and during invasion . The GO terms “cyclic-nucleotide phosphodiesterase activity” (GO:0004112) and “calcium ion binding” (GO:0005509) also point towards an involvement in invasion and motility of these proteins. Additionally, ortholog cluster with E. falciformis genes annotated as “pyrophosphatase” emerged in the Coccidia. Pyrophosphatase localizes to acidocalcisomes, special apical organelles storing calcium in T. gondii. Cyclic nucleotide signaling and calcium-dependent signal transduction are considered a major control element altering cytoskeleton organisation and motility of extracellular stages in order to induce host cell invasion and egress .
At the level of the genus Eimeria, DNA and protein binding were among the most enriched functions in novel genes. The AP2/ERF domains associated with predicted nucleic acid binding activity for example are known from transcription factors in T. gondiiand Plasmodium. According to our analysis we have a lack of orthologous AP2 containing genes but rather the emergence of the majority of these transcription factors at the genus level in Eimeria.
Other prominent functions divergent at the genus level were kinase and transferase activity and metal ion binding. Domains like Ankyrin, Armadillo and zinc fingers suggest further diverged control of gene-expression, signaling and protein modification.
We can thus conclude that the calcium signaling components of the invasion mechanism and the majority of the known downstream molecules are conserved throughout the Coccidia. Other aspects of signalling and control mediated by protein-nucleic acid and protein-protein binding, however, are diverged between Sarcocystidae and Eimeriidae and show conservation only within the genus Eimeria.
In T. gondii, an important role in reprogramming of host cells is ascribed to a largely expanded family of >40 kinases and pseudokinases (RopKs). Many RopKs in T. gondii are transported beyond the parasite plasma membrane in order to interact with the host’s immune system . Analyzing kinomes of T. gondii and E. tenella, Talevich et al.  found 41 and 24 RopKs, respectively. They later  refined a set of 44 Hidden Markov Models to search and classify RopKs across the Coccidia.
We found orthologs categorized as RopK for 8 genes in the genome of E. falciformis (Table 3). A search using the mentioned HMMs yielded 2 additional candidates and grouped all 10 (8 + 2) RopKs in HMM-defined families (Figure 4 and Table 3). Interestingly, the predicted RopKs were substantially enriched in potentially secreted molecules, as four contain a signal peptide for possible secretion. Only three of the other 85 general protein kinases (PF00069; score > 100) contained a secretion signal.
Close paralogs of RopK-family genes are often found in tandem repeated clusters within the genome, as shown for T. gondii Rop5 (Behnke et al., 2011). We can however rule out that the low number of E. falciformis RopKs is an artifact of sequence assembly merging highly similar paralogs into single gene models: E. falciformis RopKs are neither significantly closer to contig borders than other genes nor did we observe elevated sequencing coverage for the corresponding genomic regions. The extremely low number of RopKs in E. falciformis thus might indicate a lower number of these proteins present in rodent than in avian Eimeria parasites.
Rop18 is a highly polymorphic pathogenicity factor, whose allelic types determine the virulence of parasite strains . Rop18 of the virulent Type I T. gondii strain associates with the parasitophorous vacuole membrane of tachyzoites, where it inactivates host defence proteins, the Immunity Related GTPases (IRGs) by phosphorylation . We did not find orthologs of Rop18 in the genome of E. falciformis and searches with the respective HMM did not yield significant results. One of our predicted RopKs, however, showed some similarity to Rop18 that did not lead to inclusion in an ortholog cluster. The Rop18 similar gene was not in an ortholog cluster with an E. tenella gene and also the similarity to T. gondii Rop18 was higher than to any E. tenella gene. Recent studies show that Irgb6 is not localized to the parasitophorous vacuole of intracellular E. falciformis sporozoites . It remains to be tested whether E. falciformis needs a Rop18-like protein to release its parasitophorous vacuole from the defence mechanism executed by the IRG-system. Future studies will address if other functionally related proteins inactivate murine components in infections with E. falciformis by phosphorylation.
Neither orthologs, HMM hits or basic similar sequences were found for Rop16, another pathogenicity factor secreted by T. gondii to down regulate immune responses . It is tempting to speculate that the comparatively low number of RopKs in Eimeria is due to the fact that these parasites are specialized to a narrow spectrum of gut epithelial cells of a given host species, requiring a less diverse set of molecules to interact with the host immune system.
Additionally, Eimeria might not need to manipulate some aspects of immune responses, as these parasites typically exploit their hosts for a limited time frame, whereas T. gondii establishes life-long infections. In contrast, the large set of RopKs in T. gondii would be needed to allow infection and modulation of a large variety of host cells over a long period, as this parasite infects virtually every nucleated cell and intensively modulates host immune responses.
Intracellular parasites like E. falciformis are expected to possess a rewired metabolism with multiple heterotrophies, as they thrive on nutrients of host origin. To identify such metabolic peculiarities, we assessed the metabolic potential of E. falciformis with a combination of methods employing simple comparisons of domain and metabolic enzyme annotations but also automated genome scale metabolic reconstruction based on MetaCyc  and KEGG  pathways.
387 distinct EC numbers (in full detail to the 4th digit) for metabolic enzyme activities were obtained for 670 genes: for 301 through domain associations in InterProScan, 220 by ortholog assignment to T. gondii LAMP annotations  and 193 by SwissProt sequence similarity searches while controlling for domain agreement (Additional file 14).
We created two different versions of EfaCyc, our MetaCyc- derived model of the E. falciformis metabolism: One more accurate but incomplete and one less accurate but more comprehensive model. The comprehensive EfaCyc metabolic reconstruction comprised 271 pathways and 688 enzymes. The smaller reconstruction pruned for pathways unlikely to be present based on taxonomy listed only 125 pathways with 671 enzymes. KEGG analysis predicted 90 pathways with 325 different enzymes encoded by 365 genes.
Transporters constitute an interface between host and parasite metabolic networks. We therefore investigated especially transport annotations and also predicted transport reactions in our metabolic reconstruction . Our comprehensive EfaCyc reconstruction predicted 44 transporters for 29 transport reactions; the taxonomically pruned reconstruction recovered only 5 transport reactions and 24 transporters. KEGG analysis found 6 different transmembrane transport reactions and 8 distinct transporters.
Genes with an annotated function as transporters (GO:0005215 “transporter activity” and descendants) are more prevalent in T. gondii than in E. falciformis and slightly more prevalent in P. falciparum than in E. falciformis (Figure 5b). Especially enriched in T. gondii relative to E. falciformis were the terms amino acid-, organic acid-, organic anion-, carboxylic acid-, amine- and copper ion- transmembrane transporter activity (GO:0015171, GO:0005342, GO:0008514, GO:0046943, GO:0005275 and GO:0005375). The corresponding T. gondii genes had twelve transport-associated domains not present in any E. falciformis gene: amino acid, Ctr-copper, ABC family E, UDP-galactose, and magnesium transmembrane-, but also intraflagellar transporter domains. These suggest a richer repertoire of transmembrane transporters in T. gondii than in E. falciformis. The reduced host-range of the Eimeriidae seems a plausible reason for the loss of transporters, which would be linking host and parasite metabolic systems.
“Drug/metabolite transporter” (IPR000620), “Organic anion transporter polypeptide OATP” (IPR004156) and “Cytochrome b-561/ferric reductase transmembrane” (IPR006593) were the only transport associated domains found in E. falciformis that are absent from T. gondii, as an analysis species specific domains demonstrated (Figure 2). The drug/metabolite transporter domain absent from T. gondii is found in other Apicomplexa, the anion transporter polypeptide OATP domain is - apart from E. falciformis - only found in Trypanosomes and otherwise restricted to the Metazoa, the ferrireductase is only found outside the Apicomplexa and in E. falciformis. Based on the latter domain EfaCyc predicted a transmembrane ascorbate ferrireducatase. The proteins (2 isoforms) contain 4 and 6 transmembrane helices and are predicted to transport electrons through single membrane acting as diheme cytochrome. This is the only transporter prediction in EfaCyc which does not have an ortholog in T. gondii.
Iron uptake via lactoferrin is especially relevant in mucosal tissue, where lactoferrin is the canonical iron chelator. The transmembrane ferrireductase could be part of an alternative route for iron uptake, as part of which iron is reduced to the soluble ferrous form. A similar protein is required in Leishmania amazonensis for iron uptake from lactoferrin and replication within macrophages . We found a E. falciformis protein similar to the Leishmania iron permease Lit1 as a candidate for the ferrous iron transporter required to take up the produced ferrous iron.
We were able to identify three diversified Coccidia-specific ABC-transporters, as the KEGG module “Membrane transport” was overrepresented in genes both novel and expanded in the Coccidia. While eukaryotic ABC transporters translocate a variety of endogenous metabolites across extra- and intracellular membranes, they are also prominent for their role in xenobiotic detoxification and drug resistance .
Pathways in the KEGG module “Glycan biosynthesis and metabolism” were enriched in both ortholog clusters novel in the genus Eimeria and in the Coccidia (Figure 5a, Additional file 8). We inferred that different aspects of glycan biosynthesis show different conservation patterns.
N-Glycan biosynthesis represented by the genes beta-1,4-mannosyltransferase (ALG1; EC:22.214.171.124) and alpha-1,2-mannosyltransferase (ALG11; EC:126.96.36.199) is found in diverged ortholog groups specific to the genus Eimeria. In addition, “FAD-linked oxidase, C-terminal” (IPR004113) domains were found in three E. falciformis genes but were absent from T. gondii and all other Apicomplexa outside the genus Eimeria. The phylogenetic profile of the latter two domains and the corresponding gene-families suggests a lateral gene transfer. KEGG analysis predicted the genes containing FAD-linked oxidase domains to be part of peptidoglycan synthesis acting as UDP-N-acetylmuramate dehydrogenases (EC 188.8.131.52).
On the other hand enzymes of the pathways “Glycosylphosphatidylinositol(GPI)-anchor biosynthesis” and “Mucin type O-glycan biosynthesis” show less restricted conservation. Two subunits of the GPI transamidase complex, Phosphatidylinositol glycan, class T (PIGT) and phosphatidylinositol glycan, class U (PIGU), and N-acetylgalactosaminyltransferase (GALNT; EC:184.108.40.206) are diversified in the Coccidia but conserved throughout both Eimeriidae and Sarcocystidae. Given the between species divergence of the surface antigens it is reassuring that at least enzyme subunits involved in the synthesis of their membrane anchors are conserved within the Coccidia.
However there were also some differences in this respect between Eimeria and other Coccidia. We profiled ortholog clusters with otherwise complete apicomplexan distribution absent from all Eimeria species and scrutinized this list with basic similarity searches. Phosphomannomutase (PMM; EC:220.127.116.11) and dolichol-phosphate-mannose synthase (DPMS; EC:18.104.22.168) are prominent among the genes lost, due to their proximity in the pathways of mannose (Man) metabolism and N-glycan synthesis. Both enzymes are needed to produce activated Dol-P-Man, which serves as a substrate for N-glycan biosynthesis and the formation of glycoconjugates, among them GPI-anchored proteins.
Our results suggest a rewiring of the glycosylation pathways in the Eimeriidae as different enzymes are present, absent or show sequence divergence compared to T. gondii or other Apicomplexa. It will be interesting whether these differences are more pronounced in the synthesis of N-glycans than of GPI-anchors, as suggested by our analysis.
We have used the presented reconstructions of the E. falciformis metabolism with the clear scope to identify missing or additional metabolic modules and focused especially on transporters. A refinement of the constructed metabolic network and a combination with gene expression data from host and parasite, will allow further analysis of heterotrophies related to the uptake of host metabolites. First candidate pathways with missing intermediate reactions, divergent enzymes and enzymes acquired by lateral transfer include malate metabolism, glycosylation pathways and iron uptake.
We presented the complete genome sequence of Eimeria falciformis, a single-host parasite of the mouse. We reported on the largely contiguous reconstruction of the first Eimeria genome from a host other than chicken and predicted a reduced and remodeled set of protein coding genes in comparison to the model parasite Toxoplasma gondii (5879 vs. 8115).
Specifically, E. falciformis possesses fewer rhoptry kinases (important virulence factors) and a reduced number and diversity of transmembrane transporters in comparison to T. gondii. The difference in genome size and protein-coding gene set complexity could be a cause or consequence of the differences in host range (specialist vs. generalist) and tissue tropism (narrow vs. broad).
We observed that gene families that are involved in intracellular signaling, control of gene expression and protein modification are conserved at the level of the genus Eimeria. However, the same families diverged within the Coccidia across the two large families, the Eimeriidae and the Sarcocystidae. Mediators of invasion and motility, as well as the calcium signalling component controlling them, are however conserved and expanded in the Coccidia.
We are certain that E. falciformis with its accessible and traceable laboratory host has great potential to become a unique model for comparative host-parasite interaction analysis in Coccidia and beyond. Especially, research on avian Eimeria species with their economic importance for the poultry industry will greatly benefit from this newly established satellite model system.
We have contributed the molecular blueprint for future research on the life cycle in Eimeria and dynamics of infection in the laboratory mouse.
Sample preparation for genomic sequencing
Female NMRI mice were infected with 200 sporulated oocysts as described in Schmid et al. . At days 8 to 10 post infection faeces were collected, washed and treated with 2.5% potassium bichromate. Oocysts were allowed to sporulate for one week. For purification of sporulated oocysts, faeces were sterilized and recovered with sodium hydrochloride as described by Hoffmann and Raether (1990) and Hosek et al. (1988). Sporocysts were isolated according to the method of Kowalik and Zahner (1999) with slight modifications. Briefly, not more than 5 million sporulated oocsyts were resuspended in 0.4% pepsin-solution (Applichem), pH 3 and incubated at 37°C for 1 hour. Subsequently, sporocysts were isolated by mechanical shearing using glass beads (diameter 0.5 mm), recovered from glass beads by intense washing and seperated from oocyst cell wall components by centrifugation at 1800 g for 10 min. Sporocysts were suspended in lysis buffer and were shock frozen 3 times in liquid nitrogen. Proteins were digested with proteinase K for 45 min at 65°C. Nucleic acids were recovered by phenol/chloroform extraction and RNA was digesetd using RNase A. Genomic DNA was recovered by phenol/chloroform and chloroform extraction, followed by ethanol precipitation. The DNA pellet was resuspended in an appropriate volume of VE-water. DNA quality was assessed using an Agilent Bioanalyzer (Agilent, Santa Clara, United States). 5 μg of high quality DNA were used to prepare a paired-end DNA library according to the manufacturer’s protocol (Illumina). Additionally, 50 μg high quality DNA were subjected to a large insert mate pair preparation protocol (Illumina).
Sample preparation for RNA-Sequencing
Sporozoites were isolated from sporocysts by excystment. For this, sporocysts were incubated at 37°C in DMEM containing 0.04% tauroglycocholate (MP Biomedicals) and 0.25% trypsin (Applichem) for 30 min. Sporozoites were purified by the method of Schmatz et al. . Unsporulated oocysts were purified as described above, except for collecting the faeces every 12 hours to prevent oocysts from sporulating. For each library of merozoite/gametocyte stages the cecum of at least 3 NMRI mice was isolated. Epithelial cells were isolated as described in Schmid et al. .
Total RNA was isolated from infected epithelial cells, sporozoites and unsporulated oocysts using Trizol according to the manufacturer’s protocol (Invitrogen). High quality RNA was used to produce an mRNA library using the Illumina’s TruSeq RNA Sample Preparation guide.
Assembly and estimation of genome completeness
For genomic DNA a conventional paired-end sequencing library (fragment size 370 nt, read length 2 × 100 nt) and a large-insert mate pair library (insert size 2 kb, read length 2 × 50 nt) were sequenced. Sequencing adapters were clipped from all short reads using Flexbar . Both read sets were simultaneously assembled with Velvet 1.2.03  on a large memory machine. Apicoplast and mitochondrial genomes were reassembled from reads of first-pass contigs that aligned to organellar genomes of other Eimeria species.
We evaluated the completeness of the genome assembly by testing for the presence of a highly conserved eukaryotic core gene set (n = 248) using CEGMA  (Additional file 15). Additionally we evaluated the number of pan-apicomlexan genes based on orthologous assignment (see below) missing from our assembly. As a benchmark we performed the same analyses with publicly available datasets from EupathDB [13, 62] (version 7.3; Additional file 16) and an Eimeria maxima draft genome available in EmaxDB .
Annotation of genomic regions
We annotated transfer RNAs (tRNAs) with tRNA-Scan  and ribosomal RNAs (rRNAs) using BLAST against the Silva database . Simple Sequence Repeat (SSR) regions were identified using the MIcroSAtellite identification tool MISA (http://pgrc.ipk-gatersleben.de/misa/) and non-SSR Regions using the RepeatScout  and RepeatMasker  pipeline. Additionally, we identified transposable elements using TransposonPSI (http://transposonpsi.sourceforge.net) against profiles for proteins from (retro-) transposon families.
For annotation of protein coding genes Illumina TruSeq mRNA libraries of multiple life cycle stages were sequenced in the 2 × 100 nt format. All short reads were quality-controlled and splice-aligned to the genome sequence using TopHat 1.4.0 . Transcript structures were inferred with the Cufflinks software  and used to train the gene-finding software Augustus . Based on this training-set and additionally incorporating expression evidence from RNA-Seq, we used Augustus to predict gene-structures, coding regions and untranslated regions (UTRs).
InterProScan (version 5 RC4)  was used to predict protein domains for E. falciformis and to allow an unbiased comparison regarding annotation software also for T. gondii. InterPro domain accessions were compared for the two proteomes and enrichment in one species over the other was analyzed using Fisher’s exact test.
Annotations with Gene Ontology (GO)  terms and Enzyme Comission (EC) numbers were obtained from domain architecture if possible. Additionally, EC numbers for T. gondii genes were obtained from the hand-curated Library of Apicomplexan Metabolic Pathways (LAMP) database . Based on orthologous annotations (see OrthoMCL clustering below) these enzyme annotations were transferred to the E. falciformis proteome.
As third source of enzyme annotations we transferred functions from a BLAST search vs. SwissProt (e-value threshold 1e-5). We limited this similarity-based annotation to cases for which at least three quarters of the InterPro accessions for the Swissprot hit and the InterProScan of our E. falciformis proteins agreed.
Orthologous clusters and comparative analyses
OrthoMCL  was used with BLAST e-value cut-off of 1e-5 and inflation value of 1.5 to cluster E. falciformis proteins with the proteomes of other Apicomplexa with fully sequenced genomes (the same species used for genome comparisons (Additional file 16). OrthoMCL tries to separate out-paralogs from in-paralogs: paralogous pairs that emerged after a speciation event are contained in the same cluster whereas paralogous pairs that were present prior to a speciation event will be separated into two clusters. This usually produces a rather fine-grained clustering.
Based on these clusters proteins were phylogentically stratified and gene-set enrichment was performed for GO-terms using the R-package TopGO and InterPro domains using Fisher’s exact test.
Additionally, we constructed protein alignments for all OrthoMCL clusters using Muscle (version 3.8.31) . Alignments were trimmed using trimAl (version 1.2rev59 with trimming parameters set to “automatic”)  and phylogenetic trees were reconstructed using PhyML (version 20120412)  using the WAG model for amino-acid replacement and BIONJ as starting tree.
We excluded all C. parvum and E. maxima genes from our tree-reconciliation approach because for the first the species phylogeny in relation to the other species is not fully cleared and for the latter many alignments seemed truncated and incomplete. The remaining gene trees were reconciliated with the species tree using ranger-DTL (version ranger-dtl-U.linux 1.0)  with transfer cost set to 10, duplication cost to 2 and the cost of loss to 1. Clusters were classified as expanded if duplication events were found at a certain node. TA4 surface antigen domain containing protein alignments were built using hmmalign of HMMER3 .
For SAGs, fructose-bisphosphate aldolase and membrane alanyl aminopeptidase peptidases phylogenetic trees were built using MrBayes v3.2.1  with the mixed models setting for amino acid replacement set in priors. MCMC samples were obtained every 100 generations from 4 chains starting with random trees and running for 1,000,000 generations generations. Burnin was set to 2,500 samples, discarding the first 25% of trees. Bayesian and maximum likelihood phylogenetic trees were compared using SumTree .
EC numbers and protein coordinates were used as input to the PathwayTools Pathologic suite V17.0  and pathways were predicted based on the MetaCyc database . The transporter prediction module of PathwayTools  was used to predict transmembrane transporters and the presence of orthologs of these transporters was checked overlaying these predicitions with OrthoMCL clustering data.
Additionally we used the KEGG automatic annotation server (KAAS)  to obtain Kegg onthology (KO) identifiers and associated Kyoto Encyclopedia of Genes and Genomes (KEGG)  pathways. We tested overrepresentation of metabolic pathways in different gene-sets using fisher’s exact tests.
Production of E. falciformis and collection of infected mouse tissues are approved by the Landesamt für Gesundheit und Soziales Berlin (H0098/04).
Emanuel Heitlinger and Simone Spork as joint first authors
Richard Lucius and Christoph Dieterich as joint last authors.
Duszynski DW: Eimeria. eLS. 2011, Chichester, United Kingdom: John Wiley & Sons, Ltd
Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI: The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005, 434: 214-217. 10.1038/nature03342.
Hill D, Dubey JP: Toxoplasma gondii: transmission, diagnosis and prevention. Clin Microbiol Infect. 2002, 8: 634-640. 10.1046/j.1469-0691.2002.00485.x.
Chapman HD, Barta JR, Blake D, Gruber A, Jenkins M, Smith NC, Suo X, Tomley FM: A selective review of advances in coccidiosis research. Adv Parasitol. 2013, 83: 93-
Eimer T: Ueber Die Ei-Oder Kugelförmigen Sogenannten Psorospermien Der Wirbelthiere: Ein Beitrag Zur Entwicklungsgeschichte Der Gregarinen Und Zur Kenntniss Dieser Parasiten Als Krankheitsursache. 1870, Würzburg: Stuber, 1-58.
Haberkorn A: Die Entwicklung vonEimeria falciformis (Eimer 1870) in der weißen Maus (Mus musculus). Z Für Parasitenkd. 1970, 34: 49-67. 10.1007/BF00629179.
Owen D: Eimeria falciformis (Eimer, 1870) in specific pathogen free and gnotobiotic mice. Parasitology. 1975, 71: 293-303. 10.1017/S0031182000046734.
Pogonka T, Schelzke K, Stange J, Papadakis K, Steinfelder S, Liesenfeld O, Lucius R: CD8+ cells protect mice against reinfection with the intestinal parasite Eimeria falciformis. Microbes Infect Inst Pasteur. 2010, 12: 218-226. 10.1016/j.micinf.2009.12.005.
Stange J, Hepworth MR, Rausch S, Zajic L, Kühl AA, Uyttenhove C, Renauld J-C, Hartmann S, Lucius R: IL-22 mediates host defense against an intestinal intracellular parasite in the absence of IFN-γ at the cost of Th17-driven immunopathology. J Immunol Baltim Md 1950. 2012, 188: 2410-2418.
Schmid M, Lehmann MJ, Lucius R, Gupta N: Apicomplexan parasite, Eimeria falciformis, co-opts host tryptophan catabolism for life cycle progression in mouse. J Biol Chem. 2012, 287: 20197-20207. 10.1074/jbc.M112.351999.
Schmid M, Heitlinger E, Spork S, Mollenkopf H-J, Lucius R, Gupta N: Eimeria falciformis infection of the mouse caecum identifies opposing roles of IFNγ-regulated host pathways for the parasite development. Mucosal Immunol. 2014, 7 (4): 969-82.
Hunter CA, Sibley LD: Modulation of innate immunity by Toxoplasma gondii virulence effectors. Nat Rev Microbiol. 2012, 10: 766-778. 10.1038/nrmicro2858.
Kissinger JC, Gajria B, Li L, Paulsen IT, Roos DS: ToxoDB: accessing the Toxoplasma gondii genome. Nucleic Acids Res. 2003, 31: 234-236. 10.1093/nar/gkg072.
Reid AJ, Vermont SJ, Cotton JA, Harris D, Hill-Cawthorne GA, Könen-Waisman S, Latham SM, Mourier T, Norton R, Quail MA, Sanders M, Shanmugam D, Sohal A, Wasmuth JD, Brunk B, Grigg ME, Howard JC, Parkinson J, Roos DS, Trees AJ, Berriman M, Pain A, Wastling JM: Comparative Genomics of the Apicomplexan Parasites Toxoplasma gondii and Neospora caninum: Coccidia Differing in Host Range and Transmission Strategy. PLoS Pathog. 2012, 8: e1002567-10.1371/journal.ppat.1002567.
Walzer KA, Adomako-Ankomah Y, Dam RA, Herrmann DC, Schares G, Dubey JP, Boyle JP: Hammondia hammondi, an avirulent relative of Toxoplasma gondii, has functional orthologs of known T. gondii virulence genes. Proc Natl Acad Sci. 2013, 110: 7446-7451. 10.1073/pnas.1304322110.
Ling K-H, Rajandream M-A, Rivailler P, Ivens A, Yap S-J, Madeira AMBN, Mungall K, Billington K, Yee W-Y, Bankier AT, Carroll F, Durham AM, Peters N, Loo S-S, Isa MNM, Novaes J, Quail M, Rosli R, Nor Shamsudin M, Sobreira TJP, Tivey AR, Wai S-F, White S, Wu X, Kerhornou A, Blake D, Mohamed R, Shirley M, Gruber A, Berriman M, et al: Sequencing and analysis of chromosome 1 of Eimeria tenella reveals a unique segmental organization. Genome Res. 2007, 17: 311-319. 10.1101/gr.5823007.
Blake DP, Alias H, Billington KJ, Clark EL, Mat-Isa M-N, Mohamad A-F-H, Mohd-Amin M-R, Tay Y-L, Smith AL, Tomley FM, Wan K-L: EmaxDB: Availability of a first draft genome sequence for the apicomplexan Eimeria maxima. Mol Biochem Parasitol. 2012, 184: 48-51. 10.1016/j.molbiopara.2012.03.004.
Reid AJ, Blake DP, Ansari HR, Billington K, Browne HP, Dunn M, Hung SS, Kawahara F, Miranda-Saavedra D, Malas TB, Mourier T, Naghra H, Nair M, Otto TD, Rawlings ND, Rivailler P, Sanchez-Flores A, Sanders M, Subramaniam C, Tay Y-L, Woo Y, Wu X, Barrell B, Dear PH, Doerig C, Gruber A, Ivens AC, Parkinson J, Rajandream M-A, Shirley MW, et al: Genomic analysis of the causative agents of coccidiosis in domestic chickens. Genome Res. 2014, Epub ahead of print
Soldati D, Meissner M: Toxoplasma as a novel system for motility. Curr Opin Cell Biol. 2004, 16: 32-40. 10.1016/j.ceb.2003.11.013.
Saeij JPJ, Boyle JP, Coller S, Taylor S, Sibley LD, Brooke-Powell ET, Ajioka JW, Boothroyd JC: Polymorphic secreted kinases are key virulence factors in toxoplasmosis. Science. 2006, 314: 1780-1783. 10.1126/science.1133690.
Del Cacho E, Pages M, Gallego M, Monteagudo L, Sánchez-Acedo C: Synaptonemal complex karyotype of Eimeria tenella. Int J Parasitol. 2005, 35: 1445-1451. 10.1016/j.ijpara.2005.06.009.
Lin R-Q, Qiu L-L, Liu G-H, Wu X-Y, Weng Y-B, Xie W-Q, Hou J, Pan H, Yuan Z-G, Zou F-C, Hu M, Zhu X-Q: Characterization of the complete mitochondrial genomes of five Eimeria species from domestic chickens. Gene. 2011, 480: 28-33. 10.1016/j.gene.2011.03.004.
Cai X, Fuller AL, McDougald LR, Zhu G: Apicoplast genome of the coccidian Eimeria tenella. Gene. 2003, 321: 39-46.
Janouškovec J, Horák A, Oborník M, Lukeš J, Keeling PJ: A common red algal origin of the apicomplexan, dinoflagellate, and heterokont plastids. Proc Natl Acad Sci. 2010, 107: 10949-10954. 10.1073/pnas.1003335107.
Roy SW, Penny D: Widespread Intron Loss Suggests Retrotransposon Activity in Ancient Apicomplexans. Mol Biol Evol. 2007, 24: 1926-1933. 10.1093/molbev/msm102.
Zdobnov EM, Apweiler R: InterProScan – an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17: 847-848. 10.1093/bioinformatics/17.9.847.
Ashburner M, Ball C, Blake J, Botstein D, Butler H, Cherry J, Davis A, Dolinski K, Dwight S, Eppig J, Harris M, Hill D, Issel-Tarver L, Kasarskis A, Lewis S, Matese J, Richardson J, Ringwald M, Rubin G, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
Arabidopsis GI: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408: 796-10.1038/35048692.
Goffeau A, Barrell BG, Bussey H, Davis RW, Dujon B, Feldmann H, Galibert F, Hoheisel JD, Jacq C, Johnston M, Louis EJ, Mewes HW, Murakami Y, Philippsen P, Tettelin H, Oliver SG: Life with 6000 genes. Science. 1996, 274 (546): 563-567.
Armbrust EV, Berges JA, Bowler C, Green BR, Martinez D, Putnam NH, Zhou S, Allen AE, Apt KE, Bechner M: The genome of the diatom Thalassiosira pseudonana: ecology, evolution, and metabolism. Science. 2004, 306: 79-86. 10.1126/science.1101156.
Wasmuth J, Daub J, Peregrin-Alvarez JM, Finney CAM, Parkinson J: The origins of apicomplexan sequence innovation. Genome Res. 2009, 19: 1202-1213. 10.1101/gr.083386.108.
Li L, Stoeckert CJ, Roos DS: OrthoMCL: Identification of Ortholog Groups for Eukaryotic Genomes. Genome Res. 2003, 13: 2178-2189. 10.1101/gr.1224503.
Starnes GL, Coincon M, Sygusch J, Sibley LD: Aldolase Is Essential for Energy Production and Bridging Adhesin-Actin Cytoskeletal Interactions during Parasite Invasion of Host Cells. Cell Host Microbe. 2009, 5: 353-364. 10.1016/j.chom.2009.03.005.
Skinner-Adams TS, Stack CM, Trenholme KR, Brown CL, Grembecka J, Lowther J, Mucha A, Drag M, Kafarski P, McGowan S, Whisstock JC, Gardiner DL, Dalton JP: Plasmodium falciparum neutral aminopeptidases: new targets for anti-malarials. Trends Biochem Sci. 2010, 35: 53-61. 10.1016/j.tibs.2009.08.004.
Dos Santos SL, Freitas LM, Lobo FP, Rodrigues-Luiz GF, de Mendes TA O, Oliveira ACS, Andrade LO, Chiari É, Gazzinelli RT, Teixeira SMR, Fujiwara RT, Bartholomeu DC: The MASP Family of Trypanosoma cruzi: Changes in Gene Expression and Antigenic Profile during the Acute Phase of Experimental Infection. PLoS Negl Trop Dis. 2012, 6: e1779-10.1371/journal.pntd.0001779.
Wasmuth JD, Pszenny V, Haile S, Jansen EM, Gast AT, Sher A, Boyle JP, Boulanger MJ, Parkinson J, Grigg ME: Integrated Bioinformatic and Targeted Deletion Analyses of the SRS Gene Superfamily Identify SRS29C as a Negative Regulator of Toxoplasma Virulence. mBio. 2012, 3: e00321-12.
Brothers VM, Kuhn I, Paul LS, Gabe JD, Andrews WH, Sias SR, McCaman MT, Dragon EA, Files JG: Characterization of a surface antigen of Eimeria tenella sporozoites and synthesis from a cloned cDNA in Escherichia coli. Mol Biochem Parasitol. 1988, 28: 235-247. 10.1016/0166-6851(88)90008-4.
Tabarés E, Ferguson D, Clark J, Soon P-E, Wan K-L, Tomley F: Eimeria tenella sporozoites and merozoites differentially express glycosylphosphatidylinositol-anchored variant surface proteins. Mol Biochem Parasitol. 2004, 135: 123-132. 10.1016/j.molbiopara.2004.01.013.
Jahn D, Matros A, Bakulina AY, Tiedemann J, Schubert U, Giersberg M, Haehnel S, Zoufal K, Mock H-P, Kipriyanov SM: Model structure of the immunodominant surface antigen of Eimeria tenella identified as a target for sporozoite-neutralizing monoclonal antibody. Parasitol Res. 2009, 105: 655-668. 10.1007/s00436-009-1437-6.
Shehu K, Nowell F: Cross-reactions between Eimeria falciformis and Eimeria pragensis in mice induced by trickle infections. Parasitology. 1998, 117: 457-465. 10.1017/S0031182098003308.
Tomley FM, Soldati DS: Mix and match modules: structure and function of microneme proteins in apicomplexan parasites. Trends Parasitol. 2001, 17: 81-88. 10.1016/S1471-4922(00)01761-X.
Miller SA, Binder EM, Blackman MJ, Carruthers VB, Kim K: A conserved subtilisin-like protein TgSUB1 in microneme organelles of Toxoplasma gondii. J Biol Chem. 2001, 276: 45341-45348. 10.1074/jbc.M106665200.
Drozdowicz YM, Shaw M, Nishi M, Striepen B, Liwinski HA, Roos DS, Rea PA: Isolation and Characterization of TgVP1, a Type I Vacuolar H + −translocating Pyrophosphatase fromToxoplasma gondii - The dynamics of its subcellular localization and the cellular effects of a diphosphonate inhibitor. J Biol Chem. 2003, 278: 1075-1085. 10.1074/jbc.M209436200.
Billker O, Lourido S, Sibley LD: Calcium-Dependent Signaling and Kinases in Apicomplexan Parasites. Cell Host Microbe. 2009, 5: 612-622. 10.1016/j.chom.2009.05.017.
Walker R, Gissot M, Huot L, Alayi TD, Hot D, Marot G, Schaeffer-Reiss C, Van Dorsselaer A, Kim K, Tomavo S: Toxoplasma transcription factor TgAP2XI-5 regulates the expression of genes involved in parasite virulence and host invasion. J Biol Chem. 2013, 288: 31127-31138. 10.1074/jbc.M113.486589. published online September 10, 2013
Balaji S: Discovery of the principal specific transcription factors of Apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domains. Nucleic Acids Res. 2005, 33: 3994-4006. 10.1093/nar/gki709.
Melo MB, Jensen KDC, Saeij JPJ: Toxoplasma gondii effectors are master regulators of the inflammatory response. Trends Parasitol. 2011, 27: 487-495. 10.1016/j.pt.2011.08.001.
Talevich E, Mirza A, Kannan N: Structural and evolutionary divergence of eukaryotic protein kinases in Apicomplexa. BMC Evol Biol. 2011, 11: 321-10.1186/1471-2148-11-321.
Talevich E, Kannan N: Structural and evolutionary adaptation of rhoptry kinases and pseudokinases, a family of coccidian virulence factors. BMC Evol Biol. 2013, 13: 117-10.1186/1471-2148-13-117.
Fentress SJ, Behnke MS, Dunay IR, Mashayekhi M, Rommereim LM, Fox BA, Bzik DJ, Taylor GA, Turk BE, Lichti CF, Townsend RR, Qiu W, Hui R, Beatty WL, Sibley LD: Phosphorylation of immunity-related GTPases by a Toxoplasma gondii-secreted kinase promotes macrophage survival and virulence. Cell Host Microbe. 2010, 8: 484-495. 10.1016/j.chom.2010.11.005.
Ong Y-C, Reese ML, Boothroyd JC: Toxoplasma rhoptry protein 16 (ROP16) subverts host function by direct tyrosine phosphorylation of STAT6. J Biol Chem. 2010, 285: 28731-28740. 10.1074/jbc.M110.112359.
Caspi R, Foerster H, Fulcher CA, Kaipa P, Krummenacker M, Latendresse M, Paley S, Rhee SY, Shearer AG, Tissier C, Walk TC, Zhang P, Karp PD: The MetaCyc Database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases. Nucleic Acids Res. 2008, 36 (Database issue): D623-D631.
Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.
Shanmugasundram A, Gonzalez-Galarza FF, Wastling JM, Vasieva O, Jones AR: Library of Apicomplexan Metabolic Pathways: a manually curated database for metabolic pathways of apicomplexan parasites. Nucleic Acids Res. 2012, 41: D706-D713.
Lee TJ, Paulsen I, Karp P: Annotation-based inference of transporter function. Bioinformatics. 2008, 24: i259-i267. 10.1093/bioinformatics/btn180.
Wilson ME, Lewis TS, Miller MA, McCormick ML, Britigan BE: Leishmania chagasi: uptake of iron bound to lactoferrin or transferrin requires an iron reductase. Exp Parasitol. 2002, 100: 196-207. 10.1016/S0014-4894(02)00018-8.
Leprohon P, Légaré D, Ouellette M: ABC transporters involved in drug resistance in human parasites. Essays Biochem. 2011, 50: 121-144. 10.1042/bse0500121.
Schmatz DM, Crane MSJ, Murray PK: Purification of Eimeria Sporozoites by DE-52 Anion Exchange Chromatography. J Protozool. 1984, 31: 181-183. 10.1111/j.1550-7408.1984.tb04314.x.
Dodt M, Roehr JT, Ahmed R, Dieterich C: FLEXBAR—Flexible Barcode and Adapter Processing for Next-Generation Sequencing Platforms. Biology. 2012, 1: 895-905. 10.3390/biology1030895.
Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829. 10.1101/gr.074492.107.
Parra G, Bradnam K, Korf I: CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007, 23: 1061-1067. 10.1093/bioinformatics/btm071.
Aurrecoechea C, Heiges M, Wang H, Wang Z, Fischer S, Rhodes P, Miller J, Kraemer E, Stoeckert CJ, Roos DS, Kissinger JC: ApiDB: integrated resources for the apicomplexan bioinformatics resource center. Nucleic Acids Res. 2007, 35 (Database issue): D427-D430.
Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 0955-0964. 10.1093/nar/25.5.0955.
Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glöckner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007, 35: 7188-7196. 10.1093/nar/gkm864.
Price AL, Jones NC, Pevzner PA: De novo identification of repeat families in large genomes. Bioinformatics. 2005, 21 (suppl 1): i351-i358. 10.1093/bioinformatics/bti1018.
Tarailo-Graovac M, Chen N: Using RepeatMasker to Identify Repetitive Elements in Genomic Sequences. Curr Protoc Bioinforma. 2002, Chichester, United Kingdom: John Wiley & Sons, Ltd
Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25: 1105-1111. 10.1093/bioinformatics/btp120.
Roberts A, Pimentel H, Trapnell C, Pachter L: Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011, 27: 2325-2329. 10.1093/bioinformatics/btr355.
Stanke M, Steinkamp R, Waack S, Morgenstern B: AUGUSTUS: a web server for gene finding in eukaryotes. Nucleic Acids Res. 2004, 32 (suppl 2): W309-W312.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T: trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009, 25: 1972-1973. 10.1093/bioinformatics/btp348.
Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O: New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst Biol. 2010, 59: 307-321. 10.1093/sysbio/syq010.
Bansal MS, Alm EJ, Kellis M: Efficient algorithms for the reconciliation problem with gene duplication, horizontal transfer and loss. Bioinformatics. 2012, 28: i283-i291. 10.1093/bioinformatics/bts225.
Eddy SR: A new generation of homology search tools based on probabilistic inference. In Genome Inf Vol. 2009, 23: 205-211.
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.
Sukumaran J, Holder MT: DendroPy: a Python library for phylogenetic computing. Bioinformatics. 2010, 26: 1569-1571. 10.1093/bioinformatics/btq228.
Karp PD, Paley S, Romero P: The pathway tools software. Bioinformatics. 2002, 18 (suppl 1): S225-S232. 10.1093/bioinformatics/18.suppl_1.S225.
Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M: KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35 (suppl 2): W182-W185.
Parra G, Bradnam K, Ning Z, Keane T, Korf I: Assessing the gene space in draft genomes. Nucleic Acids Res. 2009, 37: 289-297. 10.1093/nar/gkn916.
Gardner MJ, Hall N, Fung E, White O, Berriman M, Hyman RW, Carlton JM, Pain A, Nelson KE, Bowman S, Paulsen IT, James K, Eisen JA, Rutherford K, Salzberg SL, Craig A, Kyes S, Chan M-S, Nene V, Shallom SJ, Suh B, Peterson J, Angiuoli S, Pertea M, Allen J, Selengut J, Haft D, Mather MW, Vaidya AB, Martin DMA, et al: Genome sequence of the human malaria parasite Plasmodium falciparum. Nature. 2002, 419: 498-511. 10.1038/nature01097.
We thank the sequencing service at the MDC Berlin Buch for providing access to sequencing facilities, laboratory assistant Kirsi Blank for assistance with animal handling and Nishith Gupta for comments on the manuscript. This work was supported by the Deutsche Forschungs Gemeinschaft [Graduiertenkolleg 1121, Grant to Richard Lucius].
The authors declare that they have no competing interests.
SS performed animal experiments and prepared nucleic acid samples. CD assembled the genome and predicted protein-coding genes. EH and CD performed bioinformatic analyses and interpreted the data. EH drafted the manuscript with contributions from all other authors. CD and RL conceived and coordinated the study. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Coverage from two different sequencing libraries for individual contigs of the final E. falciformis genome assembly - The read depth (coverage) estimated by mapping of raw sequencing reads back to the final assembly is given for the 2 × 100 bp short insert library on the y-axis and the 2 × 50 bp long insert (mate pair) library on the x axis. Mitochondrial and apicoplast reconstructions are highlighted in green and red color respectively, their coverage is clearly elevated. (PDF 56 KB)
Additional file 2: E. tenella (figure). Apicoplast derived reads were reassembled using velvet and the resulting contigs were compared to the E. tenella complete apicoplast genome . Sequence similarity and annotated features are visualized in the E. tenella apicoplast genome. The inferred synteny was used to reconstruct the E. falciformis apicoplast genome partitioning or reversing sequence and adding undetermined bases where necessary. (PDF 292 KB)
Additional file 3: For genome features annotated in the context of protein-coding genes using Augustus , the distribution of the per-feature GC content is given. While exons and coding sequence exons show a slightly elevated GC, introns have a lower GC. (PDF 25 KB)
Additional file 4: Protein coding features of apicomplexan genomes obtained from EupathDB (see Additional file 16 for exact versions) were compared with gene-predictions for E. falciformis obtained combining expression evidence (RNAseq) with ab initio gene-finding using Augustus . Panels show the distribution of a) the overall length of genes, b) the size of individual exons and c) introns and d) the number of coding sequence exons per gene. (PDF 144 KB)
Additional file 5: E. falciformis genes (csv table). E. falciformis genes were annotated with InterproScan (version 4 RC5) . The table is simplified and lists only distinct Interpro identifiers (IPR) for all annotated genes. (CSV 658 KB)
Additional file 6: E. falciformis genes (csv table). GO terms were obtained for E. falciformis genes from domain based annotations. The table lists all GO terms obtained for the annotated genes. (CSV 385 KB)
Additional file 7: T. gondii (csv table). T. gondii genes were annotated with InterProScan (version 4 RC5) to obtain annotations directly comparable to E. falciformis. The table is simplified and lists only distinct InterPro identifiers (IPR) for all annotated genes. (CSV 312 KB)
Enrichment analyses for ortholog clusters novel in certain clades of the apicomplexan phylogeny (figure containing tables).
Additional file 8: Gene sets identified as novel in certain clades of the apicomplexan phylogeny (Figure 3), were searched for enriched function using F-tests. The tables list enrichment results for GO molecular function, biological process, InterProScan domains and KEGG pathways. Table columns indicate IDs and names of the annotation, how often the annotation was found in the all ortholog clusterd E. falciformis genes, how many genes would be expected by chance based on this and the size of a gene-set and the p-value for the enrichment. (PDF 48 KB)
Enrichment analyses for ortholog clusters expanded in certain clades of the apicomplexan phylogeny (figure containing tables).
Additional file 9: Genes sets contained in ortholog clusters expanded at certain clades of the apicomplexan phylogeny (Figure 3), were searched for enriched function using F-tests. The tables list enrichment results for GO molecular function, biological process, InterProScan domains and KEGG pathways. Table columns indicate IDs and names of the annotation, how often the annotation was found in the all gene/species-tree reconciliated E. falciformis genes, how many genes would be expected by chance based on this and the size of a gene-set and the p-value for the enrichment. (PDF 66 KB)
Pairwise phylogenetic distance genes expanded at different nodes of the apicomplexan phylogeny (figure).
Additional file 10: Boxplots (overlayed with single datapoints) are given for the the maximal phylogenetic distance of two genes in an expanded ortholog clusters as estimated by the WAG model in PhyML . Ortholog clusters are grouped by expansions at different nodes of the apicomplexan phylogeny. (PDF 38 KB)
Additional file 11: The ortholog cluster containing E. falciformis genes annotated as fructose bisphosphate aldolase (FBA) was found to be independently expanded in both Eimeriidae and Sarcocystidae. Sequences for the ortholog group were aligned using muscle, trimmed with trimAl and a phylogenetic tree was inferred with Mister Bayes and PhyML. Labels below branches give the number of bootstrap replicates supporting the clade (out of 100), labels above branches indicate its bayesian posterior probability. For branches with 100% bootstrap support or a posterior probability of 1, labels are omitted. After gene duplication in the Sarcocystidae (FBA Sarco), both paralogous copies were retained in T. gondii and N. caninum. The same fate was experienced by the independently duplicated genes in Eimeria (FBA Eimeria), E. tenella, E. falciformis and E. maxima each retained their paralogs. (PDF 32 KB)
Additional file 12: The ortholog cluster containing E. falciformis genes annotated as M1-amylopeptidase (M1 AP) was found to be independently expanded in both Eimeriidae and Sarcocystidae. Sequences for the ortholog group were aligned using muscle, trimmed with trimAl and a phylogenetic tree was inferred with Mister Bayes and PhyML. Labels below branches give the number of bootstrap replicates supporting the clade (out of 100), labels above branches indicate its bayesian posterior probability. For branches with 100% bootstrap support or a posterior probability of 1 labels are omitted. A first duplication common to the Coccidia, preceding the split of the Sarcocystidae and Eimeriidae created two copies of the gene. In one of the resulting subclades (M1AP -A) N. caninum and T. gondii genes were conserved in one copy, the E. tenella gene was lost and the E. falciformis gene duplicated. In the other subclade (M1AP - B) the gene was further duplicated after the split of the Sarcocystidae and Eimeriidae leaving two T. gondii paralogs after the loss of one ortholog in N. caninum. Similarly, in the Eimeriidae after initial duplication one of the paralogs was lost in E. tenella, and the E. falciformis paralog expanded in an additional duplication. (PDF 33 KB)
Additional file 13: Surface antigen domain (TA4) containing genes are enriched in both E. falciformis and E. tenella expanded gene families. Phylogenetic trees were constructed based a HMM guided alignment using MrBayes and PhyML and merged. The nodes separating avian and rodent SAG clades are highlighted with stars, within these supporting bootstrap replicates (from 100) and bayesian posterior probabilities are given. This tree confirms the independent expansion of two clades in E. falciformis and E. tenella and further highlights one E. falciformis gene with higher cross-species similarity to avian than rodent Eimeria. (PDF 35 KB)
Additional file 14: The table lists annotations with EC numbers for E. falciformis genes, the source for these annotations (see Methods for details) and the level of detail (to which digit in the EC hierarchy) of the corresponding annotation. (CSV 40 KB)
Additional file 15: Apicomplexan genomes were searched for HMMs of Core Eukaryote Genes (CEGs) using CEGMA . Parra et al.  divided these CEGs into four groups according to the degree of conservation observed in the pairwise alignments (4 being the most conserved group, 1 the least). Apicomplexan genome is general show a high divergence of CEGs visible in the low recovery of less conserved groups. CEGs for which proteins are longer than 70% of the corresponding HMM alignment are recognized as fully covered, others as partially covered. Group 4, expected to be fully present even in the reduced genomes of Apicomplexa is with 97% (63 genes) represented nearly completely in our assembly. The representation of less conserved groups 1–3 compares well with the high quality genomes of Plasmodium falciparum and T. gondii. (PDF 18 KB)
Additional file 16: Genome sources used (xls table). The species and the associated genome information used in various comparative genomics analyses. Columns indicate: The whole genome sequence (in nucleotide fasta format) as it was used to estimate genome completeness (CEGMA). The annotation files (in gff format) used to compare basic genomic features and the predicted proteins (in amino acid fasta format) used for ortholog clustering (OrthoMCL) are followed by the name of the database these data was obtained from. (XLS 9 KB)
About this article
Cite this article
Heitlinger, E., Spork, S., Lucius, R. et al. The genome of Eimeria falciformis - reduction and specialization in a single host apicomplexan parasite. BMC Genomics 15, 696 (2014). https://doi.org/10.1186/1471-2164-15-696
- Gene Ontology
- Protein Code Gene
- Parasitophorous Vacuole
- Apicomplexan Parasite
- Eimeria Species