On the reversibility of parasitism: adaptation to a free-living lifestyle via gene acquisitions in the diplomonad Trepomonas sp. PC1
© Xu et al. 2016
Received: 20 May 2016
Accepted: 13 July 2016
Published: 1 August 2016
It is generally thought that the evolutionary transition to parasitism is irreversible because it is associated with the loss of functions needed for a free-living lifestyle. Nevertheless, free-living taxa are sometimes nested within parasite clades in phylogenetic trees, which could indicate that they are secondarily free-living. Herein, we test this hypothesis by studying the genomic basis for evolutionary transitions between lifestyles in diplomonads, a group of anaerobic eukaryotes. Most described diplomonads are intestinal parasites or commensals of various animals, but there are also free-living diplomonads found in oxygen-poor environments such as marine and freshwater sediments. All these nest well within groups of parasitic diplomonads in phylogenetic trees, suggesting that they could be secondarily free-living.
We present a transcriptome study of Trepomonas sp. PC1, a diplomonad isolated from marine sediment. Analysis of the metabolic genes revealed a number of proteins involved in degradation of the bacterial membrane and cell wall, as well as an extended set of enzymes involved in carbohydrate degradation and nucleotide metabolism. Phylogenetic analyses showed that most of the differences in metabolic capacity between free-living Trepomonas and the parasitic diplomonads are due to recent acquisitions of bacterial genes via gene transfer. Interestingly, one of the acquired genes encodes a ribonucleotide reductase, which frees Trepomonas from the need to scavenge deoxyribonucleosides. The transcriptome included a gene encoding squalene-tetrahymanol cyclase. This enzyme synthesizes the sterol substitute tetrahymanol in the absence of oxygen, potentially allowing Trepomonas to thrive under anaerobic conditions as a free-living bacterivore, without depending on sterols from other eukaryotes.
Our findings are consistent with the phylogenetic evidence that the last common ancestor of diplomonads was dependent on a host and that Trepomonas has adapted secondarily to a free-living lifestyle. We believe that similar studies of other groups where free-living taxa are nested within parasites could reveal more examples of secondarily free-living eukaryotes.
KeywordsFree-living Parasite Diplomonad Dollo’s law Reversibility Trepomonas Horizontal gene transfer Ribonucleotide reductase
The word parasite originates from Greek parasitos meaning “a person who eats at the table of another”. In biology, the word is used for a relationship where an organism (the parasite) uses resources of another organism (the host), and lives on or inside that organism. The historical view of parasites is that they are simplified versions of free-living organisms. This view is, however, outdated, as it has become increasingly clear that parasites are organisms highly adapted to their specific niches . The transition from a free-living to a parasitic lifestyle is an evolutionary process that includes the loss of some existing functions as well as the gain of new functions needed to survive on or within the host, transmit between hosts and exploit the resources from the host . It is often argued that this evolutionary transition from a free-living state to parasitism is irreversible. The rationale is that parasites take advantage of resources from the host, leading to specialized and reductive evolution including, often, a simplified metabolism [2–6]. Once such dependence has evolved, it would seem to be nearly impossible to revert to a more complex metabolism as was found in free-living ancestors. This has sometimes been taken as an example of Dollo’s law, which states that a complex trait cannot re-evolve in the same form [7, 8].
The idea of irreversibility of parasitism is widespread; in an overview of 15 parasitology books, only four mentioned reversals to a free-living state as a possible, but unlikely, evolutionary path . However, this paradigm in biology has been questioned . Free-living house dust mites and certain nematodes have been proposed to have evolved from parasitic ancestors [9, 10]. There is strong phylogenetic support for house dust mites being secondarily free-living , but the genetics behind this lifestyle transition remains unknown. Diplomonads are another group in which free-living members may have evolved from host-associated ancestors, based on phylogenetic analyses [1, 8, 11]. Herein, we examine the hypothesis of a parasitic ancestry for free-living diplomonads using a transcriptome sequencing approach, with the aim of revealing the genomic basis and evolutionary origins of the lifestyle differences within the group.
G. intestinalis and S. salmonicida share adaptations to a parasitic lifestyle such as a large family of cysteine-rich proteins that likely provide antigenic variation to escape the host immune system and an encystation pathway that enables them to transmit between hosts [19, 20]. They are dependent on scavenging of metabolites from the host, exemplified by the absence of a ribonucleotide reductase (RNR) needed for synthesis of deoxyribonucleotides . These similarities indicate that they could share a common parasitic ancestor. Comparison between G. intestinalis and S. salmonicida also revealed differences. For example, the latter has an extended metabolic repertoire with more potential for gene regulation. This difference is likely a result of their different habitats, as Giardia lives in a stable gut environment, while S. salmonicida is adapted to a fluctuating environment .
The diplomonads studied with molecular methods can be tentatively divided into four monophyletic groups based on molecular data (Fig. 1) [22, 23]. Three of the groups (I–III) only contain taxa described as parasites or commensals (Fig. 1). Group I includes G. intestinalis, a parasite adapted to the environment of the intestine of humans and other mammals [19, 24, 25]. Groups II and III contain members of the genus Spironucleus, which are mostly dependent on fish for their survival as parasites or commensals . S. vortens is the causative agent of hole-in-the-head disease in ornamental fish [27, 28], S. barkhanus, S. salmonicida and S. salmonis are associated with wild and farmed salmonid fish [16, 18, 29–31], and S. torosa is associated with gadidae . Group IV contains both diplomonads found in association with host species, such as S. muris , and free-living diplomonads, for example (most representatives of) the genera Trepomonas and Hexamita . The fact that free-living diplomonads are exclusively found in group IV, nested within host-associated lineages with strong statistical support (Fig. 1), has invited two possible explanations: (1) that diplomonads have adapted to life within a host several times independently, or (2) that the free-living lifestyle of diplomonads such as Trepomonas and Hexamita is a secondary adaptation from a host-associated ancestor [11, 34]. If diplomonads were ancestrally free-living, we expect genomic features associated with this lifestyle, such as enzymes for degradation of bacterial prey, to be shared with free-living eukaryotes. If, on the other hand, Trepomonas recently adapted to life outside an animal host, these features should have evolved after the divergence from the lineages leading to G. intestinalis and S. salmonicida. We could annotate almost 8000 gene fragments in the Trepomonas transcriptome and among these found hundreds of genes recently acquired from various sources. Many of the laterally transferred genes are associated with degradation of phagocytosed bacteria, an important trait for a free-living heterotrophic protist.
Identification of Trepomonas sp. PC1 genes in a transcriptome of a mixed culture
We assembled 41 million single reads into 18,527 transcripts with length ≥ 200 nt using Inchworm ; 9980 genes were called from the transcripts and annotated based on the S. salmonicida genome , supplemented with information from protein and domain databases. The annotation algorithm is described in detail in the Method section. The Trepomonas culture was monoeukaryotic, but included a mixture of bacteria. Thus, mRNA extracted for sequencing was expected to contain bacterial contamination, even though the RNA was polyA-purified before library construction.
The fact that hexamitine diplomonads (groups II–IV in Fig. 1) use a non-canonical genetic code [36, 37] greatly aided the detection of contamination in the dataset. Trepomonas utilizes the codons TAA and TAG to encode glutamine instead of termination, leaving TGA as the only stop codon. This alternative genetic code has been identified in a few other eukaryotic lineages (e.g., some ciliates ), but not, to our knowledge, in any prokaryote. We excluded from further analyses a total of 1995 fragments of genes as possible contaminants because they had higher similarity to non-eukaryotic sequences than to any eukaryotic sequence and contained less than two in-frame TAA/TAG codons.
Details of the annotation
With BLAST hit
With only Pfam hit
# In-frame TAA/TAG in hit ≥ 2
# In-frame TAA/TAG ≥ 2
We used BUSCO  to estimate how complete the gene content is with the Trepomonas transcriptome; 143 of the 429 conserved eukaryotic proteins provided by BUSCO were identified within the Trepomonas transcriptome, compared to 138 and 173 for the genomes of S. salmonicida and G. intestinalis [19, 20], respectively. This suggests that the Trepomonas transcriptome dataset includes the majority of the protein-coding genes in the genome.
Trepomonas sp. PC1 has an extended coding capacity compared to parasitic diplomonads
Functional pathway differences identified in the KAAS analysis
S. salmonicida a
G. intestinalis a
Amino acid metabolism
Metabolism of other amino acids
Glycan biosynthesis and metabolism
Metabolism of cofactors and vitamins
Metabolism of terpenoids and polyketides
Biosynthesis of other secondary metabolites
Xenobiotics biodegradation and metabolism
Transport and catabolism
Gene transfer of prokaryotic genes into eukaryotic lineages is a common mechanism for adaptation that is acting on different evolutionary timescales [40–42]. For example, analyses of the genome of the moss Physcomitrella patens suggested that gene transfer was important for the plant colonization of land , and a recent study showed that transfer of genes have occurred in historical times between fungi used in cheese making . Gene transfer has been identified as an evolutionary mechanism affecting diplomonad genomes [19, 45]. Reported cases include ancient events shared with the parabasalid Trichomonas vaginalis , as well as an example of recent acquisition where an isolate of G. intestinalis encodes a functional bacterial gene flanked by two pseudogenes of bacterial origin .
We performed a phylogenomic analysis of all Trepomonas sp. PC1 transcripts to detect genes that have been gained relatively recently from non-eukaryotic sources. Genes unique to Trepomonas among sampled diplomonads would indicate probable cases of functions gained by the free-living diplomonad. We used a combination of the programs Phylogenie  and Darkhorse  to identify transfer candidates. RAxML trees were constructed for the candidates predicted by either of the programs, and were then inspected manually. We identified 423 transfer candidates that corresponded to 271 transfer events (Additional file 2: Table S1; Additional file 3). Among them there are 40 transfer events, corresponding to 61 genes, for which we could infer a putative donor lineage with bootstrap support ≥ 70 %. No dominating lineage could be identified among the putative donors, suggesting that these genes have been acquired in a large number of individual events. The majority of the transfers were from Bacteria, with only eight from Archaea and one from a virus. Among Bacteria, Proteobacteria and Firmicutes were the most common donor groups, with 12 and 10 cases, respectively. Functions associated with degradation of bacterial prey were abundant among the transfer candidates.
General lysosomal degradation
Phagocytized prey needs to be digested to recover components for anabolic processes. Eukaryotes fuse their phagocytic vacuoles with lysosomes that deliver a cocktail of hydrolytic enzymes that degrade captured prey. We found many Trepomonas proteins that are potentially involved in this process, for example, lysosomal proteases (cathepsins), glycosidases (GBA, HexA/B), palmitoyl-thioesterases, lipases (LYPLA3), nucleases (DNaseII), saposins, hydrolases, and hemolysins. Several of these proteins have been transferred from bacterial sources where they perform these functions in the absence of a lysosome (Additional file 2: Table S1).
Bactericidal permeability-increasing proteins (BPI) have antimicrobial activity in humans . These proteins bind to the lipid A moiety of lipopolysaccharides of the outer membrane of Gram-negative bacteria and penetrate into the inner membrane, thereby causing cell death. However, members of this protein family have diverse functional roles and very little is known about their function outside animals . We detect 40 BPI-like proteins in Trepomonas and they are also present in multiple copies in the S. salmonicida and Giardia genomes (Additional file 4: Figure S2). The high sequence diversity and large number of BPI-like proteins in Trepomonas sp. PC1 suggest that the family is functionally diverse. Some of these proteins may have bactericidal functions, although biochemical studies are needed to understand the role of BPI-like proteins in various diplomonads.
Transglutaminases are enzymes that catalyze formation of isopeptide bonds, although some members of this family perform the reverse reaction . The formation of isopeptide bonds may serve to entrap or clot bacteria to prevent escape from the phagosome or to establish contact at the interface between the eukaryotic microbe and a bacterium . Trepomonas encodes a set of proteins with transglutaminase domains that have been laterally transferred from bacteria and which may serve these functions (Additional file 2: Table S1). Functional studies are needed to test whether these proteins are functional within a lysosome in Trepomonas.
Degradation of bacterial cell walls
Trepomonas also encodes four bacterial N-acetylmuramoyl-L-alanine amidases (peptidoglycan amidases) that are predicted to degrade peptidoglycan cross-links at the stem peptide of MurNAc and L-alanine (Fig. 2). Peptidoglycan amidases are very common in bacteria, where they assist in the maintenance, elongation and division of the peptidoglycan sacculus, but are rarely found in eukaryotes [55, 56]. Mammals and some other animals encode homologs of peptidoglycan amidase (peptidoglycan recognition proteins) that bind to the cell wall and kill the bacterium by activating a protein-sensing two-component system . The four Trepomonas peptidoglycan amidases are found in a single cluster in the phylogenetic tree, nested among bacterial and phage proteins (Fig. 3c), indicating that they have been acquired via gene transfer.
Lysozyme-family proteins cleave β-1,4-glycosidic bonds between the MurNAc and GlcNAc groups in peptidoglycan (or between the GlcNAc residues of chitin) and thereby cause lysis of bacterial cells [55, 56]. The five Trepomonas lysozyme-family proteins are of at least two different origins (Fig. 3d, e). The first group has four members and multiple close homologs are found in the parabasalid T. vaginalis (Fig. 3d). These Trepomonas lysozyme-family proteins are likely of one or two bacterial origins because they are absent from other diplomonad genomes, although a common origin with the Trichomonas genes cannot be excluded. The second type is distantly related and found in a group of bacteria and a few sequences from fungi and water fleas (Fig. 3e). This lysozyme appears to have been acquired from bacteria by an ancestor of S. vortens and Trepomonas (Fig. 3e). A further adaptation to digest peptidoglycan constituents is suggested by the presence of a GlcNAc kinase of bacterial origin (Additional file 5: Figure S3) that might allow S. vortens and Trepomonas to activate GlcNAc for further utilization. Interestingly, activated GlcNAc may shuttle into the pathway to synthesize N-acetylgalactosamine for building a cyst wall or into glucose metabolism.
In summary, the phylogenetic analyses support the hypothesis that Trepomonas recently acquired enhanced capabilities to degrade bacterial prey through gene transfers from bacteria (Figs. 2 and 3). There are no indications from the analyses that any of these genes have a eukaryotic ancestry, except that two of the genes are present in the parabasalid T. vaginalis. However, these genes are absent in all other sampled diplomonads. Thus, shared ancestry of these genes in Trichomonas and Trepomonas would imply three independent losses of each gene within diplomonads. Therefore, independent origins in Trepomonas and Trichomonas via gene transfer appear as a more likely scenario explaining the distribution of these two genes.
Carbohydrate and amino acid metabolism
The carbohydrate and amino acid metabolism of Trepomonas sp. PC1 is similar to the predicted pathways of S. salmonicida, although Trepomonas has additional enzymes that have been acquired via lateral gene transfer (Table 2 and Additional file 2: Table S1). Trepomonas is predicted to be able to utilize glycerol as a carbon source due to possessing a laterally transferred glycerol dehydrogenase. It also seems to be able to utilize more carbon sources than the fish parasite S. salmonicida due to an extended set of 18 glycosyltransferases and glucosidases of uncertain specificity but of apparent bacterial origin. Furthermore, lateral transfers of α-amylase, α- and β-galactosidase as well as glucan endo-1,6-β-glucosidase contribute to an increased capacity to catabolize oligo- and polysaccharides. Together, this indicates an expansion of metabolic capacity to break down larger biomolecules to metabolites that feed into glycolysis.
Extended nucleotide metabolism in Trepomonas sp. PC1
All life is dependent on a supply of deoxyribonucleotides as building blocks for DNA. Almost all organisms encode at least one RNR to catalyze the reduction of ribonucleotides to deoxyribonucleotides. The very few organisms that lack RNRs are all parasites or endosymbionts, including the amoebozoan Entamoeba histolytica and the diplomonads G. intestinalis and S. salmonicida [19–21]. This suggests that these parasitic diplomonads rely on deoxyribonucleotide scavenging and nucleoside kinases in the absence of a RNR [20, 24, 58]. S. barkhanus and Spironucleus vortens encode a class III anaerobic RNR of bacterial origin (Fig. 4b), indicating that they have re-gained the capability to reduce ribonucleotides. Three Trepomonas sp. PC1 transcripts together make a putative class III anaerobic RNR of the NrdD class, fused with the RNR-activating protein NrdG (Fig. 4b). Two considerations indicate that this fused gene has a bacterial origin distinct from Spironucleus RNR. First, the Trepomonas RNR is nested among bacterial sequences in the phylogenetic analysis, with a weak affinity to Clostridium and Clostridium phage sequences (Fig. 4b). Second, the nrdD and nrdG genes are frequently found in an operon in phages and bacteria . The presence of RNR is predicted to make Trepomonas independent of scavenging of three of the four deoxynucleosides, however, the absence of thymidylate synthetase transcripts suggests that it may need a source of deoxythymidine (Fig. 4a).
Trepomonas encodes a squalene-tetrahymanol cyclase (STC)
We have used a transcriptomics approach to examine the hypothesis that the free-living lifestyle of Trepomonas is a secondary adaptation. Our analyses show that this diplomonad encodes a large number of enzymes potentially involved in degradation of prey, including several enzyme families that degrade various parts of the bacterial cell walls and that are missing from the studied parasitic diplomonads (Fig. 2). Phylogenetic analyses showed that most prey-degradation enzymes were introduced recently into the Trepomonas genome via gene transfer, as expected if the organism is secondarily free-living (Fig. 3). The increased capacity of degradation of prey is complemented with a general increase of the metabolic capacity of Trepomonas compared to parasitic diplomonads, making it capable of utilizing more metabolites (Table 2 and Additional file 2: Table S1). This evolutionary path is indeed the reverse of what has been observed for the transition from a free-living to a parasitic lifestyle; the largest differences in the coding potential of the genome of the free-living Bodo saltans compared to parasitic kinetoplastids were shown to be within macromolecular degradation and catabolism of various metabolites from bacterial pray .
The increased metabolic capability in Trepomonas compared to parasitic diplomonads is especially striking for nucleotide metabolism, where the presence of nine additional enzymes, including an RNR, indicates that Trepomonas has adapted to a life less dependent on a host association than its ancestors (Fig. 4). In addition, we detected a homolog of the gene encoding STC, the enzyme required for the synthesis of tetrahymanol, a sterol surrogate . This hallmark protein of free-living anaerobic phagotrophs possibly originated via eukaryote-to-eukaryote gene transfer (Fig. 5). We conclude that the transcriptome data favor the hypothesis that Trepomonas has adapted secondarily to a free-living lifestyle, over the alternative that the ancestral diplomonad was a free-living organism. This is supported by the organismal phylogeny where Trepomonas free-living diplomonads are nested within host-associated species (Fig. 1). However, there could exist free-living diplomonads within groups I–III, and a diplomonad normally found within a host could, in principle, have a cryptic free-living stage in its life cycle. To our knowledge, no such diplomonad has been observed. Only 10 diplomonad ribosomal RNA sequences from culture-independent surveys of environmental samples are present in the Silva ribosomal RNA gene database . All of these belong to group IV (Additional file 6: Figure S4). This is in agreement with the view that group I–III diplomonads are parasites in the sense that they are taking resources from another organism, although they do not necessarily cause their host any negative effect.
The results challenge the assumption that parasitism is irreversible. Our data suggest that the adaptation to a free-living lifestyle occurred at least partly by acquisition of bacterial genes coding for prey-degrading enzymes needed by a free-living phagotroph. Such genes could have been acquired stepwise, and may have had a selective advantage also in the intestine of a host, if the diplomonad was able to ingest bacteria. Once the diplomonad has acquired some of the genes associated with a free-living lifestyle, it may be able to grow outside the host as trophozoites for short periods. The order of the events leading to adaptation to a free-living lifestyle in the ancestors of Trepomonas sp. PC1 could be deciphered by studying additional free-living diplomonads together with their closest host-associated relatives.
Adaptation by acquisition of prokaryotic genes is common in parasitic diplomonads [19, 45], suggesting that the ancestor of Trepomonas was exposed to bacterial genes. Such frequent lateral gene transfer probably was a precondition for the evolution of secondary free-living diplomonads. Lateral gene transfer has been proposed to be important in other parasites [66–69], which hints that evolution of secondary free-living taxa by gene acquisition could be a general phenomenon. There could, indeed, be examples of secondary free-living lineages in protist groups that include important human parasites such as Entamoeba (i.e., Entamoebidae) and Trichomonas (i.e., Parabasalia). Free-living lineages are nested within parasites in the phylogenetic trees of these groups [70–72], and lateral gene transfer has been shaping the metabolism of the parasites in the groups . Interestingly, RNR, an essential enzyme for life independent of a host, has been lost in the human parasite Entamoeba histolytica , whereas we identified homologs of RNR of bacterial origins in three divergent Entamoeba species (Fig. 4b), including E. terrapinae, which is considered to be free-living . This lineage might have adapted to a free-living lifestyle secondarily, similar to Trepomonas. If so, E. terrapinae is expected to harbor more recently acquired genes associated with a free-living lifestyle. This prediction could be tested by comparative studies of Entamoeba genomes.
Transitions from parasitic to free-living lifestyles might not be restricted to protist parasites. There are strong phylogenetic indications that it has happened at least once in the evolution of dust mites , and it has been suggested that the nematode genus Rhabditophanes is secondarily free-living . Nematodes include both free-living taxa and members that are parasitic on animals or plants. The different lifestyles are admixed in the nematode phylogeny and it is estimated that parasitism has arisen at least 15 times independently . Lateral gene transfer has contributed to parasitism in nematodes [67, 74], and it is reasonable to assume that it has also contributed to adaptation in free-living nematodes. Genomic studies are needed to understand the genetic basis and evolutionary history of the different lifestyles in this animal group.
The argument against reversibility of parasitism is that it is improbable that an organism regains the same traits that were lost during the evolution to parasitism. Our study has shown that adaptation to a free-living lifestyle can occur via introduction of ‘foreign’ genes that the pre-parasitism free-living ancestor probably never encoded, thereby resolving the paradox. We believe that diplomonads and their closest relatives are a suitable group of organisms for future studies of evolutionary transitions between parasitic and free-living lifestyles.
Materials and sequencing
Trepomonas sp. PC1 was isolated from marine sediment near Peggy’s Cove, Nova Scotia, Canada, and grown in the lab on a mixed culture of bacteria. Total RNA was collected from several independently grown cultures in 50 mL Falcon tubes. Messenger RNA was purified from total RNA using the Poly(A)Purist™ MAG system (Ambion). Directed and size-fractionated cDNA libraries were made using the CloneMiner™ cDNA Library Construction Kit (Invitrogen), which were then sequenced with Sanger technology. The rest of the mRNAs were sequenced with Illumina Genome Analyzer IIx instrument, producing 41 million 100 bp long single reads. Raw RNA sequence reads were deposited at NCBI Sequence Read Archive (SRA) under accession number SRR2079337. The isolate is no longer in culture.
The adapter sequence and the low quality bases at the ends of RNA-Seq reads affected the quality of the assembly, and were trimmed by cutadapt v2.6  and prinseq v0.20.3 , respectively. Approximately 2 % of the low quality reads were removed and the remaining reads of high quality were assembled using Inchworm, the first part of Trinity . Inchworm reconstructs full-length transcripts from RNA-Seq data without creating alternatively spliced isoforms; evidence from other sequenced diplomonads indicates that alternative splicing is a rare event . Inchworm assembled 18,523 transcripts with size ≥ 200 nt. All EST sequences were found in the RNA-Seq assembly, except nine singleton EST reads that contribute no additional gene information. The annotations and analysis were thus only based on RNA-Seq data in this paper. Kmer abundance calculated by Inchworm was used to represent transcript expression.
Transcripts assembled from transcriptome data are often partial and the translating frame is not always the frame containing the longest open reading frame (ORF). BLASTX results could better capture the potential gene information in all six reading frames, however, BLASTX matches can become too fragmented to be significant. BLASTP using the longest ORF is a good complementary strategy to BLASTX when the longest ORF is the correct gene. In order to take advantage of both BLASTP and BLASTX results, we annotated the transcripts by combining the BLASTP results of the longest ORFs and the BLASTX results. We favor annotation from S. salmonicida because it is the closest sequenced relative with a manually annotated genome sequence . The UniProtKB 20130905 database  and a database containing only S. salmonicida proteins were used for BLAST separately.
Genes that are positioned very close to each other or overlap each other can be assembled into a single transcript, since we do not have strand-specific or paired-end reads. Some transcript fusions could also result from incorrect assembly. To handle both cases, we first gathered all significant annotations belonging to the same transcript, allowing annotations to be partially or non-overlapping, and then systematically analyzed the reliability of an alternative translation. Manual efforts were also employed to better understand certain transcript fusions.
We annotated 7905 genes using this approach, with 42 transcripts coding for two genes. The longest ORFs were then searched against Pfam 25.0  using HMMER3 , and 196 extra genes were added. For the rest of the transcripts that lack similarity to other genes or domains, we kept the longest ORFs with size ≥ 300 aa, which adds another 1879 genes, all of which contain in-frame TAA/TAG codons. In total, we had 9980 genes after the annotation step.
There were obvious cases of contamination among the annotated genes, which is not surprising as Trepomonas was grown with mixed bacteria. Contamination was first removed on the nucleotide sequence level. All the assembled transcripts were searched against the human genome and all the 2680 complete bacterial genomes available from NCBI (20131003) using BLASTN; 594 transcripts were determined to be from bacteria and 304 from human, using criteria that the e-value should be < 1 × 10–10 and > 50 % of the query sequence should be matched. Additional screening for contamination was then performed at the gene level. Since Trepomonas has TGA as its only stop codon, genes with in-frame TAA/TAG codons were likely to be true Trepomonas genes. A gene was excluded if it had no in-frame TAA/TAG codons, and its best BLASTP hit was to a bacterial sequence. This resulted in exclusion of 1605 likely bacterial genes. The two steps combined removed 1709 bacterial genes and 259 human genes, leaving 8012 Trepomonas genes.
All called genes with only a single in-frame TAA/TAG codon in the BLAST alignment and with best matches to prokaryotes were checked manually. The called gene was retained if the TAA/TAG coding for glutamine was in-frame and well within the BLAST alignment. If the TAA/TAG codon was in the beginning or the end of the BLAST alignment and there was no convincing alignment up- or downstream of the codon, the gene was removed from further analyses. In total, 61 genes were manually checked and 27 of them were removed. The numbers of in-frame TAA/TAG codons are reported for the genes detected as putative lateral gene transfers (Additional file 2: Table S1). In total, 7985 genes we retained for downstream analyses.
BUSCO  compares the input to its lineage-specific profile to quantitatively measure the completeness of the input in terms of the expected gene content. The lineage-specific profile “Eukaryota” containing 429 conserved eukaryotic proteins was downloaded from the BUSCO homepage, and setting “–m OGS” was used as gene set was used as the input.
This Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GDID00000000. The version described in this paper is the first version, GDID01000000. Since TSA does not accept annotation on the complementary strand or an annotation that is internally partial, the transcript contigs were reverse-complemented for annotation on the complementary strand, were split into two if they coded for two genes, and were trimmed if they were internally partial.
Protein families and orthologous groups
OrthoMCL v2.0.2  was used to identify protein families in Trepomonas, as well as homologous groups shared by Trepomonas and S. salmonicida. E-value and match cutoff were set to be 1 × 10–10 and 40 %, respectively. OrthoMCL orthologous groups with only one member from each species were included in the protein identity analysis. Protein identities were extracted from the reciprocal BLAST results between Trepomonas and S. salmonicida.
Analysis of pathways
The KEGG Automatic Annotation Server (KAAS) v1.69x using the bi-directional best hit (BBH) assignment method  was employed to predict pathways for each of Trepomonas, S. salmonicida and G. intestinalis, and the results were compared. The pathways reported in the main text were manually examined.
Screening for lateral gene transfer candidates
PhyloGenie  and Darkhorse v1.5  were combined to predict potential lateral gene transfer candidates. Both programs used the BLASTP results (against the nr database) of the predicted Trepomonas protein sequences. PhyloGenie constructs phylogenetic trees to search for gene transfer candidates. Darkhorse computes a lineage probability index to predict the potential gene transfer donors as well as recipients. This lineage probability index is supposed to be proportional to the phylogenetic distance between the database match sequence and the query sequence.
PhyloGenie uses HMMer v2.3.1  to build a Hidden Markov Model (HMM) profile from the alignment in the BLAST results, and the full-length BLAST hits were used to search against the HMM profile to further select a subset to build the alignment. Neighbor-joining trees were first generated by mtrees provided in the package, and fed into Phylome Analysis Tool (PHAT), part of PhyloGenie, to pre-select the lateral gene transfer candidates.
RAxML v8.1.15  was used to construct phylogenetic trees based on the alignments with the setting “-m PROTGAMMALG4X -f a -n 100” and PHAT was applied again to select trees of interest. Those trees were then manually inspected.
We selected trees with nodes that contain prokaryotes, diplomonads and up to four other eukaryotes, which selects the transfers from prokaryotes into diplomonads only, or into diplomonads, but shared with up to four other eukaryotes.
Phylogenie was run with the default parameter settings, except that we allowed sequences with maximum 70 % identity at the genus level (taxlevel = 1, maxsim = 0.7); and a maximum 90 % pairwise sequence similarity regardless of the species (maxhmmsim = 0.9). Further, we allowed up to 200 sequences in the alignment and 200 sequences to build the HMM profile.
During the manual inspection, donors and other eukaryotes in the same cluster were noted. Donor information at the genus level is also noted if there was a consistent genus and the bootstrap value was ≥ 70 %.
Individual gene trees
The phylogenetic trees shown in the figures were generated as follows. Taxa were manually selected from the RAxML trees and the BLAST hits. Our draft genome sequences of Giardia muris, Spironucleus barkhanus and Spironucleus vortens were searched using tBLASTn, and identified homologs were added to the datasets. MAFFT v7.215  with the L-INS-i option was used to align the sequences; BMGE v1.12  with BLOSUM30 option was used to select sites to include in the phylogenetic analyses; RAxML was used with the same settings as described earlier to reconstruct the trees, except that the number of bootstrap replicates was increased to 500.
We thank Katarina Roxström-Lindquist for assistance with the polyA selection and cDNA library construction.
This work was supported by a grant from The Swedish Research Council Formas (www.formas.se; 2010-899). Illumina data was sequenced at SNP SEQ Technology Platform in Uppsala, which is supported by Uppsala University (www.uu.se), Uppsala University Hospital (www.akademiska.se), Science for Life Laboratory (www.scilifelab.se) and the Swedish Research Council (www.vr.se). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
Raw RNA sequence reads were deposited at NCBI Sequence Read Archive (SRA) under accession number SRR2079337. This Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GDID00000000. The version described in this paper is the first version, GDID01000000.
FX, JJH, MK, SGS, and JOA conceived and designed the experiments. MK, AGBS and AJR isolated the organism and prepared nucleic acids. JJH prepared mRNA for sequencing. FX performed the bioinformatics analyses. FX, JJH, SGS, and JOA analyzed the data. FX, JJH, AGBS, SGS, and JOA wrote the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Poulin R. Evolutionary ecology of parasites. Princeton: Princeton University Press; 2007.Google Scholar
- Poulin R, Randhawa HS. Evolution of parasitism along convergent lines: from ecology to genomics. Parasitology. 2015;142:S6–15.View ArticlePubMedGoogle Scholar
- Heinz E, Williams TA, Nakjang S, Noel CJ, Swan DC, Goldberg AV, et al. The genome of the obligate intracellular parasite Trachipleistophora hominis: new insights into microsporidian genome dynamics and reductive evolution. PLoS Pathog. 2012;8:e1002979.View ArticlePubMedPubMed CentralGoogle Scholar
- Rohmer L, Hocquet D, Miller SI. Are pathogenic bacteria just looking for food? Metabolism and microbial pathogenesis. Trends Microbiol. 2011;19:341–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Jackson A, Otto T, Aslett M, Armstrong S, Bringaud F, Schlacht A, et al. Kinetoplastid phylogenomics reveals the evolutionary innovations associated with the origins of parasitism. Curr Biol. 2016;26:161–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Woo YH, Ansari H, Otto TD, Klinger CM, Kolisko M, Michálek J, et al. Chromerid genomes reveal the evolutionary path from photosynthetic algae to obligate intracellular parasites. eLife. 2015;4:e06974.View ArticlePubMedPubMed CentralGoogle Scholar
- Gould SJ. Dollo on Dollo’s law: irreversibility and the status of evolutionary laws. J Hist Biol. 1970;3:189–212.View ArticlePubMedGoogle Scholar
- Cruickshank RH, Paterson AM. The great escape: do parasites break Dollo’s law? Trends Parasitol. 2006;22:509–15.View ArticlePubMedGoogle Scholar
- Dorris M, Viney ME, Blaxter ML. Molecular phylogenetic analysis of the genus Strongyloides and related nematodes. Int J Parasitol. 2002;32:1507–17.View ArticlePubMedGoogle Scholar
- Klimov PB, OConnor B. Is permanent parasitism reversible? – critical evidence from early evolution of house dust mites. Syst Biol. 2013;62:411–23.View ArticlePubMedGoogle Scholar
- Siddall ME, Brooks DR, Desser SS. Phylogeny and the reversibility of parasitism. Evolution. 1993;47:308–13.View ArticleGoogle Scholar
- Brugerolle G, Lee JJ. Order Diplomonadida. In: Lee JJ, Leedale GF, Bradbury P, editors. An Illustrated Guide to the Protozoa. 2nd ed. Lawrence: Society of Protozoologists; 2002. p. 1125–35.
- Adl SM, Simpson AG, Lane CE, Lukes J, Bass D, Bowser SS, et al. The revised classification of eukaryotes. J Eukaryot Microbiol. 2012;59:429–514.View ArticlePubMedPubMed CentralGoogle Scholar
- Monis PT, Caccio SM, Thompson RC. Variation in Giardia: towards a taxonomic revision of the genus. Trends Parasitol. 2009;25:93–100.View ArticlePubMedGoogle Scholar
- Caccio SM, Ryan U. Molecular epidemiology of giardiasis. Mol Biochem Parasitol. 2008;160:75–80.View ArticlePubMedGoogle Scholar
- Jørgensen A, Sterud E. The marine pathogenic genotype of Spironucleus barkhanus from farmed salmonids redescribed as Spironucleus salmonicida n. sp. J Eukaryot Microbiol. 2006;53:531–41.View ArticlePubMedGoogle Scholar
- Sterud E, Mo TA, Poppe TT. Systemic spironucleosis in sea-farmed Atlantic Salmon Salmo salar, caused by Spironucleus barkhanus transmitted from feral Arctic Char Salvelinus alpinus? Dis Aquat Organ. 1998;33:63–6.View ArticlePubMedGoogle Scholar
- Sterud E, Poppe TT, Bornø G. Intracellular infection with Spironucleus barkhanus (Diplomonadida, Hexamitidae) in farmed Arctic char Salvelinus alpinus. Dis Aquat Org. 2003;56:155–61.View ArticlePubMedGoogle Scholar
- Morrison HG, McArthur AG, Gillin FD, Aley SB, Adam RD, Olsen GJ, et al. Genomic minimalism in the early diverging intestinal parasite Giardia lamblia. Science. 2007;317:1921–6.View ArticlePubMedGoogle Scholar
- Xu F, Jerlström-Hultqvist J, Einarsson E, Astvaldsson A, Svärd SG, Andersson JO. The genome of Spironucleus salmonicida highlights a fish pathogen adapted to fluctuating environments. PLoS Genet. 2014;10:e1004053.View ArticlePubMedPubMed CentralGoogle Scholar
- Lundin D, Torrents E, Poole AM, Sjoberg BM. RNRdb, a curated database of the universal enzyme family ribonucleotide reductase, reveals a high level of misannotation in sequences deposited to Genbank. BMC Genomics. 2009;10:589.View ArticlePubMedPubMed CentralGoogle Scholar
- Kolisko M, Silberman JD, Cepicka I, Yubuki N, Takishita K, Yabuki A, et al. A wide diversity of previously undetected free-living relatives of diplomonads isolated from marine/saline habitats. Environ Microbiol. 2010;12:2700–10.PubMedGoogle Scholar
- Takishita K, Kolisko M, Komatsuzaki H, Yabuki A, Inagaki Y, Cepicka I, et al. Multigene phylogenies of diverse Carpediemonas-like organisms identify the closest relatives of ‘amitochondriate’ diplomonads and retortamonads. Protist. 2012;163:344–55.View ArticlePubMedGoogle Scholar
- Adam RD. Biology of Giardia lamblia. Clin Microbiol Rev. 2001;14:447–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Ankarklev J, Jerlström-Hultqvist J, Ringqvist E, Troell K, Svärd SG. Behind the smile: cell biology and disease mechanisms of Giardia species. Nat Rev Microbiol. 2010;8:413–22.PubMedGoogle Scholar
- Williams CF, Lloyd D, Poynton SL, Jorgensen A, Millet COM, Cable J. Spironucleus species: economically-important fish pathogens and enigmatic single-celled eukaryotes. J Aquac Res Devel. 2011;S2:002. http://www.omicsonline.org/spironucleus-species-economically-important-fish-pathogens-and-enigmatic-single-celled-eukaryotes-2155-9546.S2-002.php?aid=2762.
- Poynton SL, Fraser W, Francis-Floyd R, Rutledge P, Reed P, Nerad TA. Spironucleus vortens n. sp. from fresh-water angel fish Pterophyllum scalare. Morphology and culture. J Eukaryot Microbiol. 1995;42:731–42.View ArticleGoogle Scholar
- Paull GC, Matthews RA. Spironucleus vortens, a possible cause of hole-in-the-head disease in cichlids. Dis Aquat Organ. 2001;45:197–202.View ArticlePubMedGoogle Scholar
- Sterud E, Mo TA, Poppe TT. Ultrastructure of Spironucleus barkhanus n. sp. (Diplomonadida: Hexamitidae) from grayling Thymallus thymallus (L.) (Salmonidae) and Atlantic salmon Salmo salar L (Salmonidae). J Eukaryot Microbiol. 1997;44:399–407.View ArticleGoogle Scholar
- Poynton SL, Fard MRS, Jenkins J, Ferguson HW. Ultrastructure of Spironucleus salmonis n. comb. (formerly Octomitus salmonis sensu Moore 1922, Davis 1926, and Hexamita salmonis sensu Ferguson 1979), with a guide to Spironucleus species. Dis Aquat Organ. 2004;60:49–64.View ArticlePubMedGoogle Scholar
- Fard MRS, Jorgensen A, Sterud E, Bleiss W, Poynton SL. Ultrastructure and molecular diagnosis of Spironucleus salmonis (Diplomonadida) from rainbow trout Oncorhynchus mykiss in Germany. Dis Aquat Organ. 2007;75:37–50.View ArticlePubMedGoogle Scholar
- Poynton SL, Morrison CM. Morphology of diplomonad flagellates: Spironucleus torosa n. sp. from Atlantic cod Gadus morhua L., and haddock Melanogrammus aeglefinus (L.) and Hexamita salmonis Moore from brook trout Salvelinus fontinalis (Mitchill). J Protozool. 1990;37:369–83.View ArticlePubMedGoogle Scholar
- Brett SJ, Cox FEG. Immunological aspects of Giardia muris and Spironucleus muris infections in inbred and outbred strains of laboratory mice: a comparative study. Parasitology. 1982;85:85–99.View ArticlePubMedGoogle Scholar
- Siddall ME, Hong H, Desser SS. Phylogenetic analysis of the Diplomonadida (Wenyon, 1926) Brugerolle, 1975: evidence for heterochrony in protozoa and against Giardia lamblia as a “missing link”. J Protozool. 1992;39:361–7.View ArticlePubMedGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Keeling PJ, Doolittle WF. A non-canonical genetic code in an early diverging eukaryotic lineage. EMBO J. 1996;15:2285–90.PubMedPubMed CentralGoogle Scholar
- Kolisko M, Cepicka I, Hampl V, Leigh J, Roger AJ, Kulda J, et al. Molecular phylogeny of diplomonads and enteromonads based on SSU rRNA, alpha-tubulin and HSP90 genes: implications for the evolutionary history of the double karyomastigont of diplomonads. BMC Evol Biol. 2008;8:205.View ArticlePubMedPubMed CentralGoogle Scholar
- Lozupone CA, Knight RD, Landweber LF. The molecular basis of nuclear genetic code change in ciliates. Curr Biol. 2001;11:65–74.View ArticlePubMedGoogle Scholar
- Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31:3210–2.View ArticlePubMedGoogle Scholar
- Andersson JO. Gene transfer and diversification of microbial eukaryotes. Annu Rev Microbiol. 2009;63:177–93.View ArticlePubMedGoogle Scholar
- Hirt RP, Alsmark C, Embley TM. Lateral gene transfers and the origins of the eukaryote proteome: a view from microbial parasites. Curr Opin Microbiol. 2014;23C:155–62.Google Scholar
- Soucy SM, Huang J, Gogarten JP. Horizontal gene transfer: building the web of life. Nat Rev Genet. 2015;16:472–82.View ArticlePubMedGoogle Scholar
- Yue J, Hu X, Sun H, Yang Y, Huang J. Widespread impact of horizontal gene transfer on plant colonization of land. Nat Commun. 2012;3:1152.View ArticlePubMedPubMed CentralGoogle Scholar
- Ropars J, Rodríguez dela Vega Ricardo C, López-Villavicencio M, Gouzy J, Sallet E, Dumas É, et al. Adaptive horizontal gene transfers between multiple cheese-associated fungi. Curr Biol. 2015;25:2562–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Andersson JO, Sjögren ÅM, Horner DS, Murphy CA, Dyal PL, Svärd SG, et al. A genomic survey of the fish parasite Spironucleus salmonicida indicates genomic plasticity among diplomonads and significant lateral gene transfer in eukaryote genome evolution. BMC Genomics. 2007;8:51.View ArticlePubMedPubMed CentralGoogle Scholar
- Franzén O, Jerlström-Hultqvist J, Castro E, Sherwood E, Ankarklev J, Reiner D, et al. Draft genome sequencing of Giardia intestinalis assemblage B isolate GS: are human giardiasis caused by two different species? PLoS Pathog. 2009;5(8):e1000560.View ArticlePubMedPubMed CentralGoogle Scholar
- Frickey T, Lupas AN. PhyloGenie: automated phylome generation and analysis. Nucleic Acids Res. 2004;32:5231–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Podell S, Gaasterland T. DarkHorse: a method for genome-wide prediction of horizontal gene transfer. Genome Biol. 2007;8:R16.View ArticlePubMedPubMed CentralGoogle Scholar
- Elsbach P, Weiss J. Role of the bactericidal/permeability-increasing protein in host defence. Curr Opin Immunol. 1998;10:45–9.View ArticlePubMedGoogle Scholar
- Balakrishnan A, Marathe SA, Joglekar M, Chakravortty D. Bactericidal/permeability increasing protein: A multifaceted protein with functions beyond LPS neutralization. Innate Immun. 2013;19:339–47.View ArticlePubMedGoogle Scholar
- Makarova KS, Aravind L, Koonin EV. A superfamily of archaeal, bacterial, and eukaryotic proteins homologous to animal transglutaminases. Protein Sci. 1999;8:1714–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Z, Wilhelmsson C, Hyrsl P, Loof TG, Dobes P, Klupp M, et al. Pathogen entrapment by transglutaminase—a conserved early innate immune mechanism. PLoS Pathog. 2010;6:e1000763.View ArticlePubMedPubMed CentralGoogle Scholar
- Anantharaman V, Aravind L. Evolutionary history, structural features and biochemical diversity of the NlpC/P60 superfamily of enzymes. Genome Biol. 2003;4:R11.View ArticlePubMedPubMed CentralGoogle Scholar
- Firczuk M, Bochtler M. Folds and activities of peptidoglycan amidases. FEMS Microbiol Rev. 2007;31:676–91.View ArticlePubMedGoogle Scholar
- Vollmer W, Joris B, Charlier P, Foster S. Bacterial peptidoglycan (murein) hydrolases. FEMS Microbiol Rev. 2008;32:259–86.View ArticlePubMedGoogle Scholar
- Sobhanifar S, King DT, Strynadka NC. Fortifying the wall: synthesis, regulation and degradation of bacterial peptidoglycan. Curr Opin Struct Biol. 2013;23:695–703.View ArticlePubMedGoogle Scholar
- Kashyap DR, Wang M, Liu L-H, Boons G-J, Gupta D, Dziarski R. Peptidoglycan recognition proteins kill bacteria by activating protein-sensing two-component systems. Nat Med. 2011;17:676–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Baum KF, Berens RL, Marr JJ, Harrington JA, Spector T. Purine deoxynucleoside salvage in Giardia lamblia. J Biol Chem. 1989;264:21087–90.PubMedGoogle Scholar
- Dwivedi B, Xue B, Lundin D, Edwards R, Breitbart M. A bioinformatic analysis of ribonucleotide reductase genes in phage genomes and metagenomes. BMC Evol Biol. 2013;13:33.View ArticlePubMedPubMed CentralGoogle Scholar
- Sillo A, Bloomfield G, Balest A, Balbo A, Pergolizzi B, Peracino B, et al. Genome-wide transcriptional changes induced by phagocytosis or growth on bacteria in Dictyostelium. BMC Genomics. 2008;9:291.View ArticlePubMedPubMed CentralGoogle Scholar
- Castoreno AB, Wang Y, Stockinger W, Jarzylo LA, Du H, Pagnon JC, et al. Transcriptional regulation of phagocytosis-induced membrane biogenesis by sterol regulatory element binding proteins. Proc Natl Acad Sci U S A. 2005;102:13129–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Summons RE, Bradley AS, Jahnke LL, Waldbauer JR. Steroids, triterpenoids and molecular oxygen. Philos Trans R Soc B. 2006;361:951–68.View ArticleGoogle Scholar
- Mallory FB, Gordon JT, Conner RL. The isolation of a pentacyclic triterpenoid alcohol from a protozoan. J Am Chem Soc. 1963;85:1362–3.View ArticleGoogle Scholar
- Takishita K, Chikaraishi Y, Leger MM, Kim E, Yabuki A, Ohkouchi N, et al. Lateral transfer of tetrahymanol-synthesizing genes has allowed multiple diverse eukaryote lineages to independently adapt to environments without oxygen. Biol Direct. 2012;7:5.View ArticlePubMedPubMed CentralGoogle Scholar
- Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.View ArticlePubMedGoogle Scholar
- Alsmark C, Foster PG, Sicheritz-Ponten T, Nakjang S, Martin Embley T, Hirt RP. Patterns of prokaryotic lateral gene transfers affecting parasitic microbial eukaryotes. Genome Biol. 2013;14:R19.View ArticlePubMedPubMed CentralGoogle Scholar
- Paganini J, Campan-Fournier A, Da Rocha M, Gouret P, Pontarotti P, Wajnberg E, et al. Contribution of lateral gene transfers to the genome composition and parasitic ability of root-knot nematodes. PLoS One. 2012;7:e50875.View ArticlePubMedPubMed CentralGoogle Scholar
- Boto L. Horizontal gene transfer in the acquisition of novel traits by metazoans. Proc R Soc B. 2014;281:20132450.View ArticlePubMedPubMed CentralGoogle Scholar
- Wijayawardena BK, Minchella DJ, DeWoody JA. Hosts, parasites, and horizontal gene transfer. Trends Parasitol. 2013;29:329–38.View ArticlePubMedGoogle Scholar
- Stensvold CR, Lebbad M, Victory EL, Verweij JJ, Tannich E, Alfellani M, et al. Increased sampling reveals novel lineages of Entamoeba: consequences of genetic diversity and host specificity for taxonomy and molecular detection. Protist. 2011;162:525–41.View ArticlePubMedGoogle Scholar
- Clark CG, Diamond LS. Intraspecific variation and phylogenetic relationships in the genus Entamoeba as revealed by riboprinting. J Eukaryot Microbiol. 1997;44:142–54.View ArticlePubMedGoogle Scholar
- Yubuki N, Ceza V, Cepicka I, Yabuki A, Inagaki Y, Nakayama T, et al. Cryptic diversity of free-living parabasalids, Pseudotrichomonas keilini and Lacusteria cypriaca n. g., n. sp., as inferred from small subunit rDNA sequences. J Eukaryot Microbiol. 2010;57:554–61.View ArticlePubMedGoogle Scholar
- Blaxter M, Koutsovoulos G. The evolution of parasitism in Nematoda. Parasitology. 2015;142:S26–39.View ArticlePubMedGoogle Scholar
- Wu B, Novelli J, Jiang D, Dailey HA, Landmann F, Ford L, et al. Interdomain lateral gene transfer of an essential ferrochelatase gene in human parasitic nematodes. Proc Natl Acad Sci U S A. 2013;110:7748–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998;8:175–85.View ArticlePubMedGoogle Scholar
- Gordon D, Abajian C, Green P. Consed: a graphical tool for sequence finishing. Genome Res. 1998;8:195–202.View ArticlePubMedGoogle Scholar
- Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2.View ArticleGoogle Scholar
- Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27:863–4.View ArticlePubMedPubMed CentralGoogle Scholar
- UniProt C. Reorganizing the protein space at the Universal Protein Resource (UniProt). Nucleic Acids Res. 2012;40:D71–5.View ArticleGoogle Scholar
- Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2013;42:D222–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Eddy SR. Accelerated Profile HMM Searches. PLoS Comput Biol. 2011;7:e1002195.View ArticlePubMedPubMed CentralGoogle Scholar
- Li L, Stoeckert Jr CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35:W182–5.View ArticlePubMedPubMed CentralGoogle Scholar
- HMMER: biosequence analysis using profile hidden Markov models. 2016. http://hmmer.org/. Accessed 3 Feb 2016.
- Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Criscuolo A, Gribaldo S. BMGE (Block Mapping and Gathering with Entropy): a new software for selection of phylogenetic informative regions from multiple sequence alignments. BMC Evol Biol. 2010;10:210.View ArticlePubMedPubMed CentralGoogle Scholar