Hydrophobicity-Driven Increases in Editing in Mitochondrial mRNAs during the Evolution of Kinetoplastids.

Kinetoplastids are a diverse group of flagellates which exhibit editing by insertion/deletion of Us in the mitochondrial mRNAs. Some mRNAs require editing to build most of their coding sequences, a process known as pan-editing. Evidence suggests that pan-editing is an ancestral feature in kinetoplastids. Here, we investigate how the transition from nonedited to pan-edited states occurred. The mitochondrial mRNAs and protein sequences from nine kinetoplastids and related groups (diplonemids, euglenids, and jakobids) were analyzed. RNA editing increased protein hydrophobicity to extreme values by introducing Us in the second codon position, despite the absence of editing preferences related to codon position. In addition, hydrophobicity was maintained by purifying selection in species that lost editing by retroposition of the fully edited mRNA. Only a few hydrophobic to hydrophilic amino acid changes were inferred for such species. In the protein secondary structure, these changes occurred spatially close to other hydrophilic residues. The analysis of coevolving sites showed that multiple changes are required together for hydrophobicity to be lost, which suggest the proteins are locked into extended hydrophobicity. Finally, an analysis of the NAD7 protein-protein interactions showed they can also influence hydrophobicity increase in the protein and where editing can occur in the mRNA. In conclusion, our results suggest that protein hydrophobicity has influenced editing site selection and how editing expanded in mRNAs. In effect, the hydrophobicity increase was entrenched by a neutral ratchet moved by a mutational pressure to introduce Us, thus helping to explain both RNA editing increase and, possibly, persistence.

Phylogenomics of novel ploeotid taxa contribute to the backbone of the euglenid tree.

Euglenids are a diverse group of flagellates that inhabit most environments and exhibit many different nutritional modes. The most prominent euglenids are phototrophs, but phagotrophs constitute the majority of phylogenetic diversity of euglenids. They are pivotal to our understanding of euglenid evolution, yet we are only starting to understand relationships amongst phagotrophs, with the backbone of the tree being most elusive. Ploeotids make up most of this backbone diversity-yet despite their morphological similarities, SSU rDNA analyses and multigene analyses show that they are non-monophyletic. As more ploeotid diversity is sampled, known taxa have coalesced into some subgroups (e.g. Alistosa), but the relationships amongst these are not always supported and some taxa remain unsampled for multigene phylogenetics. Here, we used light microscopy and single-cell transcriptomics to characterize five ploeotid euglenids and place them into a multigene phylogenetic framework. Our analyses place Decastava in Alistosa; while Hemiolia branches with Liburna, establishing the novel clade Karavia. We describe Hemiolia limna, a freshwater-dwelling species in an otherwise marine clade. Intriguingly, two undescribed ploeotids are found to occupy pivotal positions in the tree: Chelandium granulatum nov. gen. nov. sp. branches as sister to Olkasia, and Gaulosia striata nov. gen. nov. sp. remains an orphan taxon.

Molecular phylogenetics of sessile Dolium sedentarium, a petalomonad euglenid.

The euglenids are a species-rich group of flagellates with varying modes of nutrition that can be found in diverse habitats. Phagotrophic members of this group gave rise to phototrophs and hold the key to understanding the evolution of euglenids as a whole, including the evolution of complex morphological characters like the euglenid pellicle. Yet to understand the evolution of these characters, a comprehensive sampling of molecular data is needed to correlate morphological and molecular data, and to estimate a basic phylogenetic backbone of the group. While the availability of SSU rDNA and, more recently, multigene data from phagotrophic euglenids has improved, several "orphan" taxa remain without any molecular data whatsoever. Dolium sedentarium is one such taxon: It is a rarely-observed phagotrophic euglenid that inhabits tropical benthic environments and is one of few known sessile euglenids. Based on morphological characters, it has been thought of as part of the earliest branch of euglenids, the Petalomonadida. We report the first molecular sequencing data for Dolium using single-cell transcriptomics, adding another small piece in the puzzle of euglenid evolution. Both SSU rDNA and multigene phylogenies confirm it as a solitary branch within Petalomonadida.

18S rDNA sequence-structure phylogeny of the Euglenophyceae (Euglenozoa, Euglenida).

The phylogeny of Euglenophyceae (Euglenozoa, Euglenida) has been discussed for decades with new genera being described in the last few years. In this study, we reconstruct a phylogeny using 18S rDNA sequence and structural data simultaneously. Using homology modeling, individual secondary structures were predicted. Sequence-structure data are encoded and automatically aligned. Here, we present a sequence-structure neighbor-joining tree of more than 300 taxa classified as Euglenophyceae. Profile neighbor-joining was used to resolve the basal branching pattern. Neighbor-joining, maximum parsimony, and maximum likelihood analyses were performed using sequence-structure information for manually chosen subsets. All analyses supported the monophyly of Eutreptiella, Discoplastis, Lepocinclis, Strombomonas, Cryptoglena, Monomorphina, Euglenaria, and Colacium. Well-supported topologies were generally consistent with previous studies using a combined dataset of genetic markers. Our study supports the simultaneous use of sequence and structural data to reconstruct more accurate and robust trees. The average bootstrap value is significantly higher than the average bootstrap value obtained from sequence-only analyses, which is promising for resolving relationships between more closely related taxa.

Maturyoshka: A maturase inside a maturase, and other peculiarities of the novel chloroplast genomes of marine euglenophytes.

Organellar genomes often carry group II introns, which occasionally encode proteins called maturases that are important for splicing. The number of introns varies substantially among various organellar genomes, and bursts of introns have been observed in multiple eukaryotic lineages, including euglenophytes, with more than 100 introns in their plastid genomes. To examine the evolutionary diversity and history of maturases, an essential gene family among euglenophytes, we searched for their homologs in newly sequenced and published plastid genomes representing all major euglenophyte lineages. We found that maturase content in plastid genomes has a patchy distribution, with a maximum of eight of them present in Eutreptiella eupharyngea. The most basal lineages of euglenophytes, Eutreptiales, share the highest number of maturases, but the lowest number of introns. We also identified a peculiar convoluted structure of a gene located in an intron, in a gene within an intron, within yet another gene, present in some Eutreptiales. Further investigation of functional domains of identified maturases show that most of them lost at least one of the functional domains, which implies that the patchy maturase distribution is due to frequent inactivation and eventual loss over time. Finally, we identified the diversified evolutionary origin of analysed maturases, which were acquired along with the green algal plastid or horizontally transferred. These findings indicate that euglenophytes' plastid maturases have experienced a surprisingly dynamic history due to gains from diversified donors, their retention, and loss.

Toward the robust resolution of taxonomic ambiguity within Lepocinclis (Euglenida) based on DNA sequencing and morphology.

DNA sequences were analyzed for three groups of species from the Lepocinclis genus (L. acus-like, L. oxyuris-like, and L. tripteris-like) along with cellular morphology. Phylogenetic analyses were based on nuclear SSU rDNA, LSU rDNA, and plastid-encoded LSU rDNA. DNA sequences were obtained from species available in culture collections (L. acus SAG 1224-1a and UTEX 1316) and those isolated directly from the environment in Poland (48 isolates), resulting in 79 new sequences. The obtained phylogenetic tree of Lepocinclis included 27 taxa, five of which are presented for the first time (L. convoluta, L. gracillimoides, L. longissima, L. pseudospiroides, and L. torta) and nine taxonomically verified and described. Based on morphology, literature data, and phylogenetic analyses, the following species were distinguished: in the L. acus-like group, L. longissima and L. acus; in the L. tripteris-like group, L. pseudospiroides, L. torta, and L. tripteris; in the L. oxyuris-like group, L. gracillimoides, L. oxyuris var. oxyuris, and L. oxyuris var. helicoidea. For all verified species, diagnostic descriptions were emended, nomenclatural adjustments were made, and epitypes were designated.

Single-cell genomics unveils a canonical origin of the diverse mitochondrial genomes of euglenozoans.

BACKGROUND: The supergroup Euglenozoa unites heterotrophic flagellates from three major clades, kinetoplastids, diplonemids, and euglenids, each of which exhibits extremely divergent mitochondrial characteristics. Mitochondrial genomes (mtDNAs) of euglenids comprise multiple linear chromosomes carrying single genes, whereas mitochondrial chromosomes are circular non-catenated in diplonemids, but circular and catenated in kinetoplastids. In diplonemids and kinetoplastids, mitochondrial mRNAs require extensive and diverse editing and/or trans-splicing to produce mature transcripts. All known euglenozoan mtDNAs exhibit extremely short mitochondrial small (rns) and large (rnl) subunit rRNA genes, and absence of tRNA genes. How these features evolved from an ancestral bacteria-like circular mitochondrial genome remains unanswered. RESULTS: We sequenced and assembled 20 euglenozoan single-cell amplified genomes (SAGs). In our phylogenetic and phylogenomic analyses, three SAGs were placed within kinetoplastids, 14 within diplonemids, one (EU2) within euglenids, and two SAGs with nearly identical small subunit rRNA gene (18S) sequences (EU17/18) branched as either a basal lineage of euglenids, or as a sister to all euglenozoans. Near-complete mitochondrial genomes were identified in EU2 and EU17/18. Surprisingly, both EU2 and EU17/18 mitochondrial contigs contained multiple genes and one tRNA gene. Furthermore, EU17/18 mtDNA possessed several features unique among euglenozoans including full-length rns and rnl genes, six mitoribosomal genes, and nad11, all likely on a single chromosome. CONCLUSIONS: Our data strongly suggest that EU17/18 is an early-branching euglenozoan with numerous ancestral mitochondrial features. Collectively these data contribute to untangling the early evolution of euglenozoan mitochondria.

Evaluation of V2 18S rDNA barcode marker and assessment of sample collection and DNA extraction methods for metabarcoding of autotrophic euglenids.

Even though the interest in metabarcoding in environmental studies is growing, euglenids are still underrepresented in both sea and freshwater bodies researches. The reason for this situation could be the unsuitability of universal eukaryotic DNA barcodes and primers as well as the lack of a verified protocol, suitable to assess euglenid diversity. In this study, using specific primers for the V2 hypervariable region of 18S rDNA for metabarcoding resulted in obtaining a high fraction (85%) of euglenid reads and species-level identification of almost 90% of them. Fifty species were detected by the metabarcoding method, including almost all species observed using a light microscope. We investigated three biomass harvesting methods (filtering, centrifugation and scraping the side of a collection vessel) and determined that centrifugation and filtration outperformed scrapes, but the choice between them is not crucial for the reliability of the analysis. In addition, eight DNA extraction methods were evaluated. We compared five commercially available DNA isolation kits, two CTAB-based protocols and a chelating resin. For this purpose, the efficiency of extraction, quality of obtained DNA, preparation time and generated costs were taken into consideration. After examination of the aforementioned criteria, we chose the GeneMATRIX Soil DNA Purification Kit as the most suitable for DNA isolation.

Multigene phylogenetics of euglenids based on single-cell transcriptomics of diverse phagotrophs.

Euglenids are a well-known group of single-celled eukaryotes, with phototrophic, osmotrophic and phagotrophic members. Phagotrophs represent most of the phylogenetic diversity of euglenids, and gave rise to the phototrophs and osmotrophs, but their evolutionary relationships are poorly understood. Symbiontids, in contrast, are anaerobes that are alternatively inferred to be derived euglenids, or a separate euglenozoan group. Most phylogenetic studies of euglenids have examined the SSU rDNA only, which is often highly divergent. Also, many phagotrophic euglenids (and symbiontids) are uncultured, restricting collection of other molecular data. We generated transcriptome data for 28 taxa, mostly using a single-cell approach, and conducted the first multigene phylogenetic analyses of euglenids to include phagotrophs and symbiontids. Euglenids are recovered as monophyletic, with symbiontids forming an independent branch within Euglenozoa. Spirocuta, the clade of flexible euglenids that contains both the phototrophs (Euglenophyceae) and osmotrophs (Aphagea), is robustly resolved, with the ploeotid Olkasia as its sister group, forming the new taxon Olkaspira. Ploeotids are paraphyletic, although Ploeotiidae (represented by Ploeotia spp.), Lentomonas, and Keelungia form a robust clade (new taxon Alistosa). Petalomonadida branches robustly as sister to other euglenids in outgroup-rooted analyses. Within Spirocuta, Euglenophyceae is a robust clade that includes Rapaza, and Anisonemia is a well-supported monophyletic group containing Anisonemidae (Anisonema and Dinema spp.), 'Heteronema II' (represented by H. vittatum), and a clade of Neometanema plus Aphagea. Among 'peranemid' phagotrophs, Chasmostoma branches with included Urceolus, and Peranema with the undescribed 'Jenningsia II', while other relationships are weakly supported and consequently the closest sister group to Euglenophyceae remains unresolved. Our results are inconsistent with recent inferences that Entosiphon is the evolutionarily pivotal sister either to other euglenids, or to Spirocuta. At least three transitions between posterior and anterior flagellar gliding occurred in euglenids, with the phylogenetic positions and directions of those transitions remaining ambiguous.

Phylogenetic and Morphological Evolution of Green Euglenophytes Based on 18S rRNA.

Green euglenophytes are a group of eukaryotes with ancient origin. In order to understand the evolution of the group, it is interesting to know which characteristics are more primitive. Here, a phylogenetic tree of green euglenophytes based on the 18S rRNA gene was constructed, and ancestral states were reconstructed based on eight morphological characters. This research clarifies the phylogenetic relationships of green euglenophytes and provides a basis for the study of the origin of these plants. The phylogenetic tree, which was constructed by Bayesian inference, revealed that: Eutreptia and Eutreptiella were sister groups and that Lepocinclis, Phacus, and Discoplastis were close relatives; Euglena, Cryptoglena, Monomorphina, and Colacium were closely related in addition to Trachelomonas and Strombomonas; and Euglena was not monophyletic. An ancestral reconstruction based on morphological characters revealed seven primitive character states: ductile surface, spirally striated, slightly narrowing or sharp elongated cauda, absence of a lorica, chloroplast lamellar, shield or large discoid, pyrenoid with sheath, and with many small paramylon grains. However, the ancestral state of the length of the flagellum could not be inferred. Euglena and Euglenaria, which both possessed all of the ancestral character states, might represent the most ancient lineages of green euglenophytes.

Description of Flexiglena gen. nov. and new members of Discoplastis and Euglenaformis (Euglenida).

Environmental sampling in Poland and the United States and phylogenetic analyses based on 567 sequences of four genes (155 sequences of nuclear SSU rDNA, 139 of nuclear LSU rDNA, 135 of plastid-encoded SSU rDNA, and 138 of plastid-encoded LSU rDNA) resulted in description of the new genus Flexiglena, which has been erected by accommodating Euglena variabilis, and enriching the Discoplastis and Euglenaformis genera with five new species. Four of them have joined the Discoplastis genus, currently consisting of six representatives: D. adunca, D. angusta (=Euglena angusta), D. constricta (=Lepocinclis constricta), D. excavata (=E. excavata), D. gasterosteus (=E. gasterosteus), and D. spathirhyncha. One of them has enriched the Euglenaformis genus, currently represented by two species: Euf. chlorophoenicea (= E. chlorophoenicea) and Euf. proxima. For most studied species, the diagnostic descriptions have been emended and epitypes were designated. Furthermore, the emending of Discoplastis and Euglenaformis diagnoses was performed.

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.

Molecular and Morphological Delimitation of Species in the Group of Lepocinclis Ovum-like taxa (Euglenida).

Although Lepocinclis ovum is recognized as a cosmopolitan and common species, and Lepocinclis globulus is the type species of the genus Lepocinclis, their correct identification is nearly impossible. The reason is that over 30 morphologically similar taxa appear in the literature, but no good diagnostic features exist to distinguish amongst them. Using environmental sampling and nuclear SSU rDNA sequencing, we delimited species within the group of Lepocinclis ovum-like taxa. Morphological and molecular features were analyzed for taxa isolated from Poland and six cultured strains from algal collections. In the case of environmental sampling, DNA was obtained from a small number of cells (20-400) isolated with a micropipette without setting up laboratory cultures (52 isolates), and phylogenetic analyses were based on the variation in nSSU rDNA. Apart from L. ovum and L. globulus, 13 other species were distinguished and four taxa (Lepocinclis conica comb. nov., L. fominii comb. nov., L. gracilicauda comb. nov., and L. pseudofominii nom. nov.) had their taxonomic ranks changed. For all verified species, diagnostic descriptions were emended and epitypes designated. The only exception was L. ovum, for which the epitype was questioned and thus, a new candidate for the epitype was suggested for future adoption.

Taxon-rich phylogeny and taxonomy of the genus Phacus (Euglenida) based on morphological and molecular data.

Morphological and molecular features were analyzed for a species of Phacus to better understand the phylogenetic relationships among them and establish the taxonomy. Phylogenetic analyses were based on nSSU rDNA and the research resulted in 55 new sequences. The study included species available in algal collections and those isolated directly from the environment in Poland and the Czech Republic. As a result, the obtained phylogenetic tree of Phacus includes 50 species, out of which 7 are represented on a tree for the first time (Phacus anacoelus, P. anomalus, P. curvicauda, P. elegans, P. lismorensis, P. minutus and P. stokesii) and many have been taxonomically verified. For all verified species, diagnostic descriptions were amended, the naming was reordered and epitypes were designated.

Secondary Plastids of Euglenids and Chlorarachniophytes Function with a Mix of Genes of Red and Green Algal Ancestry.

Endosymbiosis has been common all along eukaryotic evolution, providing opportunities for genomic and organellar innovation. Plastids are a prominent example. After the primary endosymbiosis of the cyanobacterial plastid ancestor, photosynthesis spread in many eukaryotic lineages via secondary endosymbioses involving red or green algal endosymbionts and diverse heterotrophic hosts. However, the number of secondary endosymbioses and how they occurred remain poorly understood. In particular, contrasting patterns of endosymbiotic gene transfer have been detected and subjected to various interpretations. In this context, accurate detection of endosymbiotic gene transfers is essential to avoid wrong evolutionary conclusions. We have assembled a strictly selected set of markers that provides robust phylogenomic evidence suggesting that nuclear genes involved in the function and maintenance of green secondary plastids in chlorarachniophytes and euglenids have unexpected mixed red and green algal origins. This mixed ancestry contrasts with the clear red algal origin of most nuclear genes carrying similar functions in secondary algae with red plastids.

Response of Anaerobic Protozoa to Oxygen Tension in Anaerobic System.

The effect of oxygen on anaerobic protozoa was studied in anaerobic batch reactors inoculated with sludge and protozoa cultures. Among the protozoa genera, Metopus, Brachonella, Plagiopyla, Trepomonas, and Vanella were more sensitive to oxygen compared to other genera. Protozoa genera Menoidium, Rhynchomonas, Cyclidium, Spathidium, and Amoeba were found to survive under aerobic conditions, and the growth rate was slightly higher or similar to anaerobic condition. O2 tension resulted in the loss of free and endosymbiotic methanogens in anaerobic system, while methanogens were observed inside the protozoan cysts. Survival of anaerobic protozoa declined considerably when the O2 tension exceeded 1% atm. sat. and showed chemosensory behavior in response to O2 exposure. Superoxide dismutase activity was detected in survived protozoa cells under O2 tension. Facultative anaerobic protozoa with SOD activity can provide a mechanism to overcome possible occurrence of oxygen toxicity in the treatment of wastewater in anaerobic reactor.

Protist Biodiversity and Biogeography in Lakes From Four Brazilian River-Floodplain Systems.

The biodiversity and biogeography of protists inhabiting many ecosystems have been intensely studied using different sequencing approaches, but tropical ecosystems are relatively under-studied. Here, we sampled planktonic waters from 32 lakes associated with four different river-floodplains systems in Brazil, and sequenced the DNA using a metabarcoding approach with general eukaryotic primers. The lakes were dominated by the largely free-living Discoba (mostly the Euglenida), Ciliophora, and Ochrophyta. There was low community similarity between lakes even within the same river-floodplain. The protists inhabiting these floodplain systems comprise part of the large and relatively undiscovered diversity in the tropics.

Spatial and temporal variations in composition of algae assemblages with environmental variables in an urban stream (Ankara, Turkey).

Phytoplankton and epipelon assemblages form the main constituents, and they are producers in aquatic ecosystems, such as streams and rivers. This study was carried out between May 2008 and April 2009 to determine the impacts of polluted water on species variations, compositions, and community metrics in phytoplankton and epipelon at six stations on Ankara Stream. A total of 231 taxa were recorded during the study period, with 131 Bacillariophyta, 3 Charophyta, 41 Chlorophyta, 30 Cyanobacteria, 25 Euglenophyta, and 1 Ochrophyta. Heterogeneity of the stream stations was determined by the use of hierarchical cluster analysis (HCA). Community metrics were compared by using non-parametric tests, while canonical correspondence analysis (CCA) was used for the relationships between environmental variables and species. Variations in water quality and species composition along the stream flow revealed a significant spatial heterogeneity (p < 0.05). However, the upper stations of the stream were represented by unpolluted water quality with low nutrients and conductivity, and the mid- and downstream stations were characterized by high concentrations of ammonia (up to 60 mg L-1) and o-phosphate (up to 25 mg/L), with low concentrations of dissolved oxygen (< 1 mg L-1). The results, clearly supported by indicator taxa, showed that various domestic and industrial discharges affected the increase in pollution and the spatial heterogeneity. The findings obtained in this study will contribute to future improvements in Ankara Stream watershed studies.

Intrageneric Variability Between the Chloroplast Genomes of Trachelomonas grandis and Trachelomonas volvocina and Phylogenomic Analysis of Phototrophic Euglenoids.

The latest studies of chloroplast genomes of phototrophic euglenoids yielded different results according to intrageneric variability such as cluster arrangement or diversity of introns. Although the genera Euglena and Monomorphina in those studies show high syntenic arrangements at the intrageneric level, the two investigated Eutreptiella species comprise low synteny. Furthermore Trachelomonas volvocina show low synteny to the chloroplast genomes of the sister genera Monomorphina aenigmatica, M. parapyrum, Cryptoglena skujae, Euglenaria anabaena, Strombomonas acuminata, all of which were highly syntenic. Consequently, this study aims at the analysis of the cpGenome of Trachelomonas grandis and a comparative examination of T. volvocina to investigate whether the cpGenomes are of such resemblance as could be expected for a genus within the Euglenaceae. Although these analyses resulted in almost identical gene content to other Euglenaceae, the chloroplast genome showed significant novelties: In the rRNA operon, we detected group II introns, not yet found in any other cpGenome of Euglenaceae and a substantially heterogeneous cluster arrangement in the genus Trachelomonas. The phylogenomic analysis with 84 genes of 19 phototrophic euglenoids and 18 cpGenome sequences from Chlorophyta and Streptophyta resulted in a well-supported cpGenome phylogeny, which is in accordance to former phylogenetic analyses.