Gleaning Euglenozoa-specific DNA polymerases in public single-cell transcriptome data

Multiple genes encoding family A DNA polymerases (famA DNAPs), which are evolutionary relatives of DNA polymerase I (PolI) in bacteria and phages, have been found in eukaryotic genomes, and many of these proteins are used mainly in organelles. Among members of the phylum Euglenozoa, distinct types of famA DNAP, PolIA, PolIBCD+, POP, and eugPolA, have been found. It is intriguing how the suite of famA DNAPs had been established during the evolution of Euglenozoa, but the DNAP data have not been sampled from the taxa that sufficiently represent the diversity of this phylum. In particular, little sequence data were available for basal branching species in Euglenozoa until recently. Thanks to the single-cell transcriptome data from symbiontids and phagotrophic euglenids, we have an opportunity to cover the "hole" in the repertory of famA DNAPs in the deep branches in Euglenozoa. The current study identified 16 new famA DNAP sequences in the transcriptome data from 33 phagotrophic euglenids and two symbiontids, respectively. Based on the new famA DNAP sequences, the updated diversity and evolution of famA DNAPs in Euglenozoa are discussed.

DESCRIPTION, BIOLOGY, AND MEDICAL SIGNIFICANCE OF LEISHMANIA (MUNDINIA) CHANCEI N. SP. (KINETOPLASTEA: TRYPANOSOMATIDAE) FROM GHANA AND LEISHMANIA (MUNDINIA) PROCAVIENSIS N. SP. (KINETOPLASTEA: TRYPANOSOMATIDAE) FROM NAMIBIA

Genetic and phylogenetic analysis was performed on 2 isolates of Leishmania using DNA sequence data from the RNA polymerase II large subunit gene and the ribosomal protein L23a intergenic sequence. This showed the isolates to represent 2 new species within the subgenus Leishmania (Mundinia). The addition of Leishmania (Mundinia) chancei and Leishmania (Mundinia) procaviensis creates a total of 6 named species to date within this recently described subgenus of parasitic protozoa, containing both human pathogens and nonpathogens. Their widespread geographical distribution, basal phylogenetic position within the genus Leishmania, and probable non-sand fly vectors make these L. (Mundinia) species of significant medical and biological interest.

The Glycosylphosphatidylinositol Anchor: A Linchpin for Cell Surface Versatility of Trypanosomatids

The use of glycosylphosphatidylinositol (GPI) to anchor proteins to the cell surface is widespread among eukaryotes. The GPI-anchor is covalently attached to the C-terminus of a protein and mediates the protein’s attachment to the outer leaflet of the lipid bilayer. GPI-anchored proteins have a wide range of functions, including acting as receptors, transporters, and adhesion molecules. In unicellular eukaryotic parasites, abundantly expressed GPI-anchored proteins are major virulence factors, which support infection and survival within distinct host environments. While, for example, the variant surface glycoprotein (VSG) is the major component of the cell surface of the bloodstream form of African trypanosomes, procyclin is the most abundant protein of the procyclic form which is found in the invertebrate host, the tsetse fly vector. Trypanosoma cruzi, on the other hand, expresses a variety of GPI-anchored molecules on their cell surface, such as mucins, that interact with their hosts. The latter is also true for Leishmania, which use GPI anchors to display, amongst others, lipophosphoglycans on their surface. Clearly, GPI-anchoring is a common feature in trypanosomatids and the fact that it has been maintained throughout eukaryote evolution indicates its adaptive value. Here, we explore and discuss GPI anchors as universal evolutionary building blocks that support the great variety of surface molecules of trypanosomatids.

Development of Trypanosoma everetti in Culicoides biting midges

Trypanosoma species (Trypanosomatida, Kinetoplastea) are almost exclusively heteroxenous flagellated parasites, which have been extensively studied as the causative agents of severe trypanosomiasis in humans and domestic animals. However, the biology of avian trypanosomes remains insufficiently known, particularly in wildlife, despite information that some species might be pathogenic and affect the fitness of intensively infected individuals. Avian trypanosomes are cosmopolitans. Due to regular bird seasonal migrations, this host-parasite system might provide new insight for better understanding mechanisms of transcontinental dispersal of pathogens, their ecological plasticity, specificity and speciation. Trypanosoma everetti parasitizes numerous bird species globally, but data on its biology are scarce and its vectors remain unknown. This study aimed to test experimentally whether widespread Culicoides (Diptera: Ceratopogonidae) biting midges are susceptible to infection with this parasite. Two common house martins Delichon urbicum and two sedge warblers Acrocephalus schoenobaenus naturally infected with T. everetti were caught in the wild after arrival from African wintering grounds. Laboratory reared Culicoides nubeculosus and wild-caught Culicoides impunctatus biting midges were exposed by allowing them to take infected blood meals. The experimentally infected and control insects were maintained in the laboratory and dissected at intervals to follow the development of the parasite. Infections were determined using microscopic examination and PCR-based testing. Four closely related haplotypes of T. everetti were found, and each was present in different individual parasite-donor birds. These parasites readily developed and produced metacyclic trypomastigotes in C. nubeculosus and C. impunctatus biting midges. Molecular characterisation of T. everetti was developed. According to Bayesian phylogenetic analysis using a DNA fragment encoding 18S rRNA, the five species of small avian trypanosomes were closely related. Wild caught Culicoides biting midges were also collected and screened for the presence of natural infections. In all, 6.8% of wild-caught biting midges belonging to five Culicoides species were PCR-positive for kinetoplastids, including Trypanosoma species. Culicoides biting midges are readily susceptible and likely naturally transmit avian trypanosomes and thus, should be targeted in epidemiology research of avian trypanosomiasis.

Inventory and Evolution of Mitochondrion-localized Family A DNA Polymerases in Euglenozoa

The order Trypanosomatida has been well studied due to its pathogenicity and the unique biology of the mitochondrion. In Trypanosoma brucei, four DNA polymerases, namely PolIA, PolIB, PolIC, and PolID, related to bacterial DNA polymerase I (PolI), were shown to be localized in mitochondria experimentally. These mitochondrion-localized DNA polymerases are phylogenetically distinct from other family A DNA polymerases, such as bacterial PolI, DNA polymerase gamma (Polγ) in human and yeasts, “plant and protist organellar DNA polymerase (POP)” in diverse eukaryotes. However, the diversity of mitochondrion-localized DNA polymerases in Euglenozoa other than Trypanosomatida is poorly understood. In this study, we discovered putative mitochondrion-localized DNA polymerases in broad members of three major classes of Euglenozoa—Kinetoplastea, Diplonemea, and Euglenida—to explore the origin and evolution of trypanosomatid PolIA-D. We unveiled distinct inventories of mitochondrion-localized DNA polymerases in the three classes: (1) PolIA is ubiquitous across the three euglenozoan classes, (2) PolIB, C, and D are restricted in kinetoplastids, (3) new types of mitochondrion-localized DNA polymerases were identified in a prokinetoplastid and diplonemids, and (4) evolutionarily distinct types of POP were found in euglenids. We finally propose scenarios to explain the inventories of mitochondrion-localized DNA polymerases in Kinetoplastea, Diplonemea, and Euglenida.

Evolution of metabolic capabilities and molecular features of diplonemids, kinetoplastids, and euglenids

Background

The Euglenozoa are a protist group with an especially rich history of evolutionary diversity. They include diplonemids, representing arguably the most species-rich clade of marine planktonic eukaryotes; trypanosomatids, which are notorious parasites of medical and veterinary importance; and free-living euglenids. These different lifestyles, and particularly the transition from free-living to parasitic, likely require different metabolic capabilities. We carried out a comparative genomic analysis across euglenozoan diversity to see how changing repertoires of enzymes and structural features correspond to major changes in lifestyles.

Results

We find a gradual loss of genes encoding enzymes in the evolution of kinetoplastids, rather than a sudden decrease in metabolic capabilities corresponding to the origin of parasitism, while diplonemids and euglenids maintain more metabolic versatility. Distinctive characteristics of molecular machines such as kinetochores and the pre-replication complex that were previously considered specific to parasitic kinetoplastids were also identified in their free-living relatives. Therefore, we argue that they represent an ancestral rather than a derived state, as thought until the present. We also found evidence of ancient redundancy in systems such as NADPH-dependent thiol-redox. Only the genus Euglena possesses the combination of trypanothione-, glutathione-, and thioredoxin-based systems supposedly present in the euglenozoan common ancestor, while other representatives of the phylum have lost one or two of these systems. Lastly, we identified convergent losses of specific metabolic capabilities between free-living kinetoplastids and ciliates. Although this observation requires further examination, it suggests that certain eukaryotic lineages are predisposed to such convergent losses of key enzymes or whole pathways.

Conclusions

The loss of metabolic capabilities might not be associated with the switch to parasitic lifestyle in kinetoplastids, and the presence of a highly divergent (or unconventional) kinetochore machinery might not be restricted to this protist group. The data derived from the transcriptomes of free-living early branching prokinetoplastids suggests that the pre-replication complex of Trypanosomatidae is a highly divergent version of the conventional machinery. Our findings shed light on trends in the evolution of metabolism in protists in general and open multiple avenues for future research.

Diversity and Evolution of Paramoeba spp. and their Kinetoplastid Endosymbionts

Members of the genus Paramoeba (including Neoparamoeba) (Amoebozoa) are single-celled eukaryotes of economic and ecological importance because of their association with disease in a variety of marine animals including fish, sea urchins, and lobster. Interestingly, they harbor a eukaryotic endosymbiont of kinetoplastid ancestry, Perkinsela sp. To investigate the complex relationship between Paramoeba spp. and Perkinsela sp., as well as the relationships between different Paramoeba species, molecular data was obtained for four novel isolates. We also acquired new data from the urchin pathogen P. invadens. Comprehensive molecular phylogenetic analyses were carried out using 33 newly obtained 18S rDNA sequences from the host amoebae and 16 new 18S rDNA sequences from their corresponding Perkinsela sp., together with all publicly available 18S molecular data. Intra-isolate 18S rDNA nucleotide diversity was found to be surprisingly high within the various species of Paramoeba, but relatively low within their Perkinsela sp. endosymbionts. 18S rDNA phylogenies and ParaFit co-evolution analysis revealed a high degree of congruence between the Paramoeba and Perkinsela sp. tree topologies, strongly suggesting that a single endosymbiotic event occurred in the common ancestor of known Paramoeba species, and that the endosymbionts have been inherited vertically ever since.

The Evolutionary Loss of RNAi Key Determinants in Kinetoplastids as a Multiple Sporadic Phenomenon

We screened the genomes of a broad panel of kinetoplastid protists for genes encoding proteins associated with the RNA interference (RNAi) system using probes from the Argonaute (AGO1), Dicer1 (DCL1), and Dicer2 (DCL2) genes of Leishmania brasiliensis and Crithidia fasciculata. We identified homologs for all the three of these genes in the genomes of a subset of these organisms. However, several of these organisms lacked evidence for any of these genes, while others lacked only DCL2. The open reading frames encoding these putative proteins were structurally analyzed in silico. The alignments indicated that the genes are homologous with a high degree of confidence, and three-dimensional structural models strongly supported a functional relationship to previously characterized AGO1, DCL1, and DCL2 proteins. Phylogenetic analysis of these putative proteins showed that these genes, when present, evolved in parallel with other nuclear genes, arguing that the RNAi system genes share a common progenitor, likely across all Kinetoplastea. In addition, the genome segments bearing these genes are highly conserved and syntenic, even among those taxa in which they are absent. However, taxa in which these genes are apparently absent represent several widely divergent branches of kinetoplastids, arguing that these genes were independently lost at least six times in the evolutionary history of these organisms. The mechanisms responsible for the apparent coordinate loss of these RNAi system genes independently in several lineages of kinetoplastids, while being maintained in other related lineages, are currently unknown.

A Comparative In Silico Study of the Antioxidant Defense Gene Repertoire of Distinct Lifestyle Trypanosomatid Species

Kinetoplastids are an ancestral group of protists that contains free-living species and parasites with distinct mechanisms in response to stress. Here, we compared genes involved in antioxidant defense (AD), proposing an evolution model among trypanosomatids. All genes were identified in Bodo saltans, suggesting that AD mechanisms have evolved prior to adaptation for parasitic lifestyles. While most of the monoxenous and dixenous parasites revealed minor differences from B. saltans, the endosymbiont-bearing species have an increased number of genes. The absence of these genes was mainly observed in the extracellular parasites of the genera Phytomonas and Trypanosoma. In trypanosomes, a distinction was observed between stercorarian and salivarian parasites, except for Trypanosoma rangeli. Our analyses indicate that the variability of AD among trypanosomatids at the genomic level is not solely due to the geographical isolation, being mainly related to specific adaptations of their distinct biological cycles within insect vectors and to a parasitism of a wide range of hosts.

Heme pathway evolution in kinetoplastid protists

Background

Kinetoplastea is a diverse protist lineage composed of several of the most successful parasites on Earth, organisms whose metabolisms have coevolved with those of the organisms they infect. Parasitic kinetoplastids have emerged from free-living, non-pathogenic ancestors on multiple occasions during the evolutionary history of the group. Interestingly, in both parasitic and free-living kinetoplastids, the heme pathway—a core metabolic pathway in a wide range of organisms—is incomplete or entirely absent. Indeed, Kinetoplastea investigated thus far seem to bypass the need for heme biosynthesis by acquiring heme or intermediate metabolites directly from their environment.

Results

Here we report the existence of a near-complete heme biosynthetic pathway in Perkinsela spp., kinetoplastids that live as obligate endosymbionts inside amoebozoans belonging to the genus Paramoeba/Neoparamoeba. We also use phylogenetic analysis to infer the evolution of the heme pathway in Kinetoplastea.

Conclusion

We show that Perkinsela spp. is a deep-branching kinetoplastid lineage, and that lateral gene transfer has played a role in the evolution of heme biosynthesis in Perkinsela spp. and other Kinetoplastea. We also discuss the significance of the presence of seven of eight heme pathway genes in the Perkinsela genome as it relates to its endosymbiotic relationship with Paramoeba.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0664-6) contains supplementary material, which is available to authorized users.

The influence of environmental factors on protistan microorganisms in grassland soils along a land-use gradient

In this study, we investigated the effect of land use intensity, soil parameters and vegetation on protistan communities in grassland soils. We performed qualitative (T-RFLP) and quantitative (qPCR) analyses using primers specifically targeting the 18S rRNA gene for all Eukarya and for two common flagellate groups, i.e. the Chrysophyceae and the Kinetoplastea. Both approaches were applied to extracted soil DNA and RNA, in order to distinguish between the potentially active protists (i.e. RNA pool) and the total protistan communities, including potentially inactive and encysted cells (i.e. DNA pool). Several environmental determinants such as site, soil parameters and vegetation had an impact on the T-RFLP community profiles and the abundance of the quantified 18S rRNA genes. Correlating factors often differed between quantitative (qPCR) and qualitative (T-RFLP) approaches. For instance the Chrysophyceae/Eukarya 18S rDNA ratio as determined by qPCR correlated with the C/N ratio, whereas the community composition based on T-RLFP analysis was not affected indicating that both methods taken together provide a more complete picture of the parameters driving protist diversity. Moreover, distinct T-RFs were obtained, which could serve as potential indicators for either active organisms or environmental conditions like water content. While site was the main determinant across all investigated exploratories, land use seemed to be of minor importance for structuring protist communities. The impact of other parameters differed between the target groups, e.g. Kinetoplastea reacted on changes to water content on all sites, whereas Chrysophyceae were only affected in the Schorfheide. Finally, in most cases different responses were observed on RNA- and DNA-level, respectively. Vegetation for instance influenced the two flagellate groups only at the DNA-level across all sites. Future studies should thus include different protistan groups and also distinguish between active and inactive cells, in order to reveal causal shifts in community composition and abundance in agriculturally used systems.

Characterization of Two Species of Trypanosomatidae from the Honey Bee Apis mellifera: Crithidia mellificae Langridge and McGhee, and Lotmaria passim n. gen., n. sp

Trypanosomatids are increasingly recognized as prevalent in European honey bees (Apis mellifera) and by default are attributed to one recognized species, Crithidia mellificae Langridge and McGhee, 1967. We provide reference genetic and ultrastructural data for type isolates of C. mellificae (ATCC 30254 and 30862) in comparison with two recent isolates from A. mellifera (BRL and SF). Phylogenetics unambiguously identify strains BRL/SF as a novel taxonomic unit distinct from C. mellificae strains 30254/30862 and assign all four strains as lineages of a novel clade within the subfamily Leishmaniinae. In vivo analyses show strains BRL/SF preferably colonize the hindgut, lining the lumen as adherent spheroids in a manner identical to previous descriptions from C. mellificae. Microscopy images show motile forms of C. mellificae are distinct from strains BRL/SF. We propose the binomial Lotmaria passim n. gen., n. sp. for this previously undescribed taxon. Analyses of new and previously accessioned genetic data show C. mellificae is still extant in bee populations, however, L. passim n. gen., n. sp. is currently the predominant trypanosomatid in A. mellifera globally. Our findings require that previous reports of C. mellificae be reconsidered and that subsequent trypanosomatid species designations from Hymenoptera provide genetic support.

Primer pairs for the specific environmental detection and T-RFLP analysis of the ubiquitous flagellate taxa Chrysophyceae and Kinetoplastea

Bacterivorous protists play a key role in microbial soil food webs, however due to the lack of specific PCR protocols targeting selected protist taxa, knowledge on the diversity and dynamics of these groups is scarce. We developed specific PCR primers in combination with a T-RFLP protocol for the cultivation-independent analysis of two important taxa of bacterivorous flagellates, the Chrysophyceae and Kinetoplastea, in soil samples. Sequence analysis of clone libraries originating from two soils in temperate regions demonstrated the specificity of the respective primer pairs. Clone sequences affiliating to the Chrysophyceae mainly clustered within the clade C2, which has been known so far for its presence mainly in cold climatic regions, whereas Kinetoplastea sequences were mainly related to the Neobodonid clade. Based on an in silico restriction analysis of database sequence entries, suitable restriction enzymes for a T-RFLP approach were selected. This in silico approach revealed the necessity to use a combination of two restriction enzymes for T-RFLP analysis of the Chrysophyceae. Soil T-RFLP profiles reflected all T-RFs of the clone library sequences obtained from the same soils and allowed to distinguish flagellate communities from different sites. We propose to use these primer pairs for PCR detection and rapid fingerprint screening in environmental samples and envisage their use also for quantitative PCR or next generation sequencing approaches.