Euglenozoan kleptoplasty illuminates the early evolution of photoendosymbiosis.

Kleptoplasts (kP) are distinct among photosynthetic organelles in eukaryotes (i.e., plastids) because they are routinely sequestered from prey algal cells and function only temporarily in the new host cell. Therefore, the hosts of kleptoplasts benefit from photosynthesis without constitutive photoendosymbiosis. Here, we report that the euglenozoan Rapaza viridis has only kleptoplasts derived from a specific strain of green alga, Tetraselmis sp., but no canonical plastids like those found in its sister group, the Euglenophyceae. R. viridis showed a dynamic change in the accumulation of cytosolic polysaccharides in response to light-dark cycles, and 13C isotopic labeling of ambient bicarbonate demonstrated that these polysaccharides originate in situ via photosynthesis; these data indicate that the kleptoplasts of R. viridis are functionally active. We also identified 276 sequences encoding putative plastid-targeting proteins and 35 sequences of presumed kleptoplast transporters in the transcriptome of R. viridis. These genes originated in a wide range of algae other than Tetraselmis sp., the source of the kleptoplasts, suggesting a long history of repeated horizontal gene transfer events from different algal prey cells. Many of the kleptoplast proteins, as well as the protein-targeting system, in R. viridis were shared with members of the Euglenophyceae, providing evidence that the early evolutionary stages in the green alga-derived secondary plastids of euglenophytes also involved kleptoplasty.

Genomics of Preaxostyla Flagellates Illuminates the Path Towards the Loss of Mitochondria.

The notion that mitochondria cannot be lost was shattered with the report of an oxymonad Monocercomonoides exilis, the first eukaryote arguably without any mitochondrion. Yet, questions remain about whether this extends beyond the single species and how this transition took place. The Oxymonadida is a group of gut endobionts taxonomically housed in the Preaxostyla which also contains free-living flagellates of the genera Trimastix and Paratrimastix. The latter two taxa harbour conspicuous mitochondrion-related organelles (MROs). Here we report high-quality genome and transcriptome assemblies of two Preaxostyla representatives, the free-living Paratrimastix pyriformis and the oxymonad Blattamonas nauphoetae. We performed thorough comparisons among all available genomic and transcriptomic data of Preaxostyla to further decipher the evolutionary changes towards amitochondriality, endobiosis, and unstacked Golgi. Our results provide insights into the metabolic and endomembrane evolution, but most strikingly the data confirm the complete loss of mitochondria for all three oxymonad species investigated (M. exilis, B. nauphoetae, and Streblomastix strix), suggesting the amitochondriate status is common to a large part if not the whole group of Oxymonadida. This observation moves this unique loss to 100 MYA when oxymonad lineage diversified.

Contrasting outcomes of genome reduction in mikrocytids and microsporidians.

BACKGROUND: Intracellular symbionts often undergo genome reduction, losing both coding and non-coding DNA in a process that ultimately produces small, gene-dense genomes with few genes. Among eukaryotes, an extreme example is found in microsporidians, which are anaerobic, obligate intracellular parasites related to fungi that have the smallest nuclear genomes known (except for the relic nucleomorphs of some secondary plastids). Mikrocytids are superficially similar to microsporidians: they are also small, reduced, obligate parasites; however, as they belong to a very different branch of the tree of eukaryotes, the rhizarians, such similarities must have evolved in parallel. Since little genomic data are available from mikrocytids, we assembled a draft genome of the type species, Mikrocytos mackini, and compared the genomic architecture and content of microsporidians and mikrocytids to identify common characteristics of reduction and possible convergent evolution. RESULTS: At the coarsest level, the genome of M. mackini does not exhibit signs of extreme genome reduction; at 49.7 Mbp with 14,372 genes, the assembly is much larger and gene-rich than those of microsporidians. However, much of the genomic sequence and most (8075) of the protein-coding genes code for transposons, and may not contribute much of functional relevance to the parasite. Indeed, the energy and carbon metabolism of M. mackini share several similarities with those of microsporidians. Overall, the predicted proteome involved in cellular functions is quite reduced and gene sequences are extremely divergent. Microsporidians and mikrocytids also share highly reduced spliceosomes that have retained a strikingly similar subset of proteins despite having reduced independently. In contrast, the spliceosomal introns in mikrocytids are very different from those of microsporidians in that they are numerous, conserved in sequence, and constrained to an exceptionally narrow size range (all 16 or 17 nucleotides long) at the shortest extreme of known intron lengths. CONCLUSIONS: Nuclear genome reduction has taken place many times and has proceeded along different routes in different lineages. Mikrocytids show a mix of similarities and differences with other extreme cases, including uncoupling the actual size of a genome with its functional reduction.

Challenging the Importance of Plastid Genome Structure Conservation: New Insights From Euglenophytes.

Plastids, similar to mitochondria, are organelles of endosymbiotic origin, which retained their vestigial genomes (ptDNA). Their unique architecture, commonly referred to as the quadripartite (four-part) structure, is considered to be strictly conserved; however, the bulk of our knowledge on their variability and evolutionary transformations comes from studies of the primary plastids of green algae and land plants. To broaden our perspective, we obtained seven new ptDNA sequences from freshwater species of photosynthetic euglenids-a group that obtained secondary plastids, known to have dynamically evolving genome structure, via endosymbiosis with a green alga. Our analyses have demonstrated that the evolutionary history of euglenid plastid genome structure is exceptionally convoluted, with a patchy distribution of inverted ribosomal operon (rDNA) repeats, as well as several independent acquisitions of tandemly repeated rDNA copies. Moreover, we have shown that inverted repeats in euglenid ptDNA do not share their genome-stabilizing property documented in chlorophytes. We hypothesize that the degeneration of the quadripartite structure of euglenid plastid genomes is connected to the group II intron expansion. These findings challenge the current global paradigms of plastid genome architecture evolution and underscore the often-underestimated divergence between the functionality of shared traits in primary and complex plastid organelles.

Tiara: deep learning-based classification system for eukaryotic sequences.

MOTIVATION: With a large number of metagenomic datasets becoming available, eukaryotic metagenomics emerged as a new challenge. The proper classification of eukaryotic nuclear and organellar genomes is an essential step toward a better understanding of eukaryotic diversity. RESULTS: We developed Tiara, a deep-learning-based approach for the identification of eukaryotic sequences in the metagenomic datasets. Its two-step classification process enables the classification of nuclear and organellar eukaryotic fractions and subsequently divides organellar sequences into plastidial and mitochondrial. Using the test dataset, we have shown that Tiara performed similarly to EukRep for prokaryotes classification and outperformed it for eukaryotes classification with lower calculation time. In the tests on the real data, Tiara performed better than EukRep in analyzing the small dataset representing eukaryotic cell microbiome and large dataset from the pelagic zone of oceans. Tiara is also the only available tool correctly classifying organellar sequences, which was confirmed by the recovery of nearly complete plastid and mitochondrial genomes from the test data and real metagenomic data. AVAILABILITY AND IMPLEMENTATION: Tiara is implemented in python 3.8, available at https://github.com/ibe-uw/tiara and tested on Unix-based systems. It is released under an open-source MIT license and documentation is available at https://ibe-uw.github.io/tiara. Version 1.0.1 of Tiara has been used for all benchmarks. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

New Parabasalia symbionts Snyderella spp. and Daimonympha gen. nov. from South American Rugitermes termites and the parallel evolution of a cell with a rotating "head".

Most Parabasalia are symbionts in the hindgut of "lower" (non-Termitidae) termites, where they widely vary in morphology and degree of morphological complexity. Large and complex cells in the class Cristamonadea evolved by replicating a fundamental unit, the karyomastigont, in various ways. We describe here four new species of Calonymphidae (Cristamonadea) from Rugitermes hosts, assigned to the genus Snyderella based on diagnostic features (including the karyomastigont pattern) and molecular phylogeny. We also report a new genus of Calonymphidae, Daimonympha, from Rugitermes laticollis. Daimonympha's morphology does not match that of any known Parabasalia, and its SSU rRNA gene sequence corroborates this distinction. Daimonympha does however share a puzzling feature with a few previously described, but distantly related, Cristamonadea: a rapid, smooth, and continuous rotation of the anterior end of the cell, including the many karyomastigont nuclei. The function of this rotatory movement, the cellular mechanisms enabling it, and the way the cell deals with the consequent cell membrane shear, are all unknown. "Rotating wheel" structures are famously rare in biology, with prokaryotic flagella being the main exception; these mysterious spinning cells found only among Parabasalia are another, far less understood, example.

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.

Monophyly of diverse Bigyromonadea and their impact on phylogenomic relationships within stramenopiles.

Stramenopiles are a diverse but relatively well-studied eukaryotic supergroup with considerable genomic information available (Sibbald and Archibald, 2017). Nevertheless, the relationships between major stramenopile subgroups remain unresolved, in part due to a lack of data from small nanoflagellates that make up a lot of the genetic diversity of the group. This is most obvious in Bigyromonadea, which is one of four major stramenopile subgroups but represented by a single transcriptome. To examine the diversity of Bigyromonadea and how the lack of data affects the tree, we generated transcriptomes from seven novel bigyromonada species described in this study: Develocauda condao n. gen. n. sp., Develocanicus komovi n. gen. n. sp., Develocanicus vyazemskyi n. sp., Cubaremonas variflagellatum n. gen. n. sp., Pirsonia chemainus nom. prov., Feodosia pseudopoda nom. prov., and Koktebelia satura nom. prov. Both maximum likelihood and Bayesian phylogenomic trees based on a 247 gene-matrix recovered a monophyletic Bigyromonadea that includes two diverse subgroups, Developea and Pirsoniales, that were not previously related based on single gene trees. Maximum likelihood analyses show Bigyromonadea related to oomycetes, whereas Bayesian analyses and topology testing were inconclusive. We observed similarities between the novel bigyromonad species and motile zoospores of oomycetes in morphology and the ability to self-aggregate. Rare formation of pseudopods and fused cells were also observed, traits that are also found in members of labyrinthulomycetes, another osmotrophic stramenopiles. Furthermore, we report the first case of eukaryovory in the flagellated stages of Pirsoniales. These analyses reveal new diversity of Bigyromonadea, and altogether suggest their monophyly with oomycetes, collectively known as Pseudofungi, is the most likely topology of the stramenopile tree.

The Oxymonad Genome Displays Canonical Eukaryotic Complexity in the Absence of a Mitochondrion.

The discovery that the protist Monocercomonoides exilis completely lacks mitochondria demonstrates that these organelles are not absolutely essential to eukaryotic cells. However, the degree to which the metabolism and cellular systems of this organism have adapted to the loss of mitochondria is unknown. Here, we report an extensive analysis of the M. exilis genome to address this question. Unexpectedly, we find that M. exilis genome structure and content is similar in complexity to other eukaryotes and less "reduced" than genomes of some other protists from the Metamonada group to which it belongs. Furthermore, the predicted cytoskeletal systems, the organization of endomembrane systems, and biosynthetic pathways also display canonical eukaryotic complexity. The only apparent preadaptation that permitted the loss of mitochondria was the acquisition of the SUF system for Fe-S cluster assembly and the loss of glycine cleavage system. Changes in other systems, including in amino acid metabolism and oxidative stress response, were coincident with the loss of mitochondria but are likely adaptations to the microaerophilic and endobiotic niche rather than the mitochondrial loss per se. Apart from the lack of mitochondria and peroxisomes, we show that M. exilis is a fully elaborated eukaryotic cell that is a promising model system in which eukaryotic cell biology can be investigated in the absence of mitochondria.

Revisions to the Classification, Nomenclature, and Diversity of Eukaryotes.

This revision of the classification of eukaryotes follows that of Adl et al., 2012 [J. Euk. Microbiol. 59(5)] and retains an emphasis on protists. Changes since have improved the resolution of many nodes in phylogenetic analyses. For some clades even families are being clearly resolved. As we had predicted, environmental sampling in the intervening years has massively increased the genetic information at hand. Consequently, we have discovered novel clades, exciting new genera and uncovered a massive species level diversity beyond the morphological species descriptions. Several clades known from environmental samples only have now found their home. Sampling soils, deeper marine waters and the deep sea will continue to fill us with surprises. The main changes in this revision are the confirmation that eukaryotes form at least two domains, the loss of monophyly in the Excavata, robust support for the Haptista and Cryptista. We provide suggested primer sets for DNA sequences from environmental samples that are effective for each clade. We have provided a guide to trophic functional guilds in an appendix, to facilitate the interpretation of environmental samples, and a standardized taxonomic guide for East Asian users.

Dictyochophyceae Plastid Genomes Reveal Unusual Variability in Their Organization.

Dictyochophyceae (silicoflagellates) are unicellular freshwater and marine algae (Heterokontophyta, stramenopiles). Despite their abundance in global oceans and potential ecological significance, discovered in recent years, neither nuclear nor organellar genomes of representatives of this group were sequenced until now. Here, we present the first complete plastid genome sequences of Dictyochophyceae, obtained from four species: Dictyocha speculum, Rhizochromulina marina, Florenciella parvula and Pseudopedinella elastica. Despite their comparable size and genetic content, these four plastid genomes exhibit variability in their organization: plastid genomes of F. parvula and P. elastica possess conventional quadripartite structure with a pair of inverted repeats, R. marina instead possesses two direct repeats with the same orientation and D. speculum possesses no repeats at all. We also observed a number of unusual traits in the plastid genome of D. speculum, including expansion of the intergenic regions, presence of an intron in the otherwise non-intron-bearing psaA gene, and an additional copy of the large subunit of RuBisCO gene (rbcL), the last of which has never been observed in any plastid genome. We conclude that despite noticeable gene content similarities between the plastid genomes of Dictyochophyceae and their relatives (pelagophytes, diatoms), the number of distinctive features observed in this lineage strongly suggests that additional taxa require further investigation.

cytb as a New Genetic Marker for Differentiation of Prototheca Species.

Achlorophyllous unicellular microalgae of the genus Prototheca (Trebouxiophyceae, Chlorophyta) are the only known plants that cause infections in both humans and animals, collectively referred to as protothecosis. Human protothecosis, most commonly manifested as cutaneous, articular, and disseminated disease, is primarily caused by Protothecawickerhamii, followed by Protothecazopfii and, sporadically, by Protothecacutis and Protothecamiyajii In veterinary medicine, however, P. zopfii is a major pathogen responsible for bovine mastitis, which is a predominant form of protothecal disease in animals. Historically, identification of Prototheca spp. has relied upon phenotypic criteria; these were later replaced by molecular typing schemes, including DNA sequencing. However, the molecular markers interrogated so far, mostly located in the ribosomal DNA (rDNA) cluster, do not provide sufficient discriminatory power to distinguish among all Prototheca spp. currently recognized. Our study is the first attempt to develop a fast, reliable, and specific molecular method allowing identification of all Prototheca spp. We propose the mitochondrial cytb gene as a new and robust marker for diagnostics and phylogenetic studies of the Prototheca algae. The cytb gene displayed important advantages over the rDNA markers. Not only did the cytb gene have the highest discriminatory capacity for resolving all Prototheca species, but it also performed best in terms of technical feasibility, understood as ease of amplification, sequencing, and multiple alignment analysis. Based on the species-specific polymorphisms in the partial cytb gene, we developed a fast and straightforward PCR-restriction fragment length polymorphism (RFLP) assay for identification and differentiation of all Prototheca species described so far. The newly proposed method is advocated to be a new gold standard in diagnostics of protothecal infections in human and animal populations.

Molecular characterization and phylogeny of four new species of the genus Trichonympha (Parabasalia, Trichonymphea) from lower termite hindguts.

Members of the genus Trichonympha are among the most well-known, recognizable and widely distributed parabasalian symbionts of lower termites and the wood-eating cockroach species of the genus Cryptocercus. Nevertheless, the species diversity of this genus is largely unknown. Molecular data have shown that the superficial morphological similarities traditionally used to identify species are inadequate, and have challenged the view that the same species of the genus Trichonympha can occur in many different host species. Ambiguities in the literature, uncertainty in identification of both symbiont and host, and incomplete samplings are limiting our understanding of the systematics, ecology and evolution of this taxon. Here we describe four closely related novel species of the genus Trichonympha collected from South American and Australian lower termites: Trichonympha hueyi sp. nov. from Rugitermes laticollis, Trichonympha deweyi sp. nov. from Glyptotermes brevicornis, Trichonympha louiei sp. nov. from Calcaritermes temnocephalus and Trichonympha webbyae sp. nov. from Rugitermes bicolor. We provide molecular barcodes to identify both the symbionts and their hosts, and infer the phylogeny of the genus Trichonympha based on small subunit rRNA gene sequences. The analysis confirms the considerable divergence of symbionts of members of the genus Cryptocercus, and shows that the two clades of the genus Trichonympha harboured by termites reflect only in part the phylogeny of their hosts.

Evolutionary Origin of Euglena.

Euglenids (Excavata, Discoba, Euglenozoa, Euglenida) is a group of free-living, single-celled flagellates living in the aquatic environments. The uniting and unique morphological feature of euglenids is the presence of a cell covering called the pellicle. The morphology and organization of the pellicle correlate well with the mode of nutrition and cell movement. Euglenids exhibit diverse modes of nutrition, including phagotrophy and photosynthesis. Photosynthetic species (Euglenophyceae) constitute a single subclade within euglenids. Their plastids embedded by three membranes arose as the result of a secondary endosymbiosis between phagotrophic eukaryovorous euglenid and the Pyramimonas-related green alga. Within photosynthetic euglenids three evolutionary lineages can be distinguished. The most basal lineage is formed by one mixotrophic species, Rapaza viridis. Other photosynthetic euglenids are split into two groups: predominantly marine Eutreptiales and freshwater Euglenales. Euglenales are divided into two families: Phacaceae, comprising three monophyletic genera (Discoplastis, Lepocinclis, Phacus) and Euglenaceae with seven monophyletic genera (Euglenaformis, Euglenaria, Colacium, Cryptoglena, Strombomonas, Trachelomonas, Monomorphina) and polyphyletic genus Euglena. For 150 years researchers have been studying Euglena based solely on morphological features what resulted in hundreds of descriptions of new taxa and many artificial intra-generic classification systems. In spite of the progress towards defining Euglena, it still remains polyphyletic and morphologically almost undistinguishable from members of the recently described genus Euglenaria; members of both genera have cells undergoing metaboly (dynamic changes in cell shape), large chloroplasts with pyrenoids and monomorphic paramylon grains. Model organisms Euglena gracilis Klebs, the species of choice for addressing fundamental questions in eukaryotic biochemistry, cell and molecular biology, is a representative of the genus Euglena.

Intermediate introns in nuclear genes of euglenids - are they a distinct type?

BACKGROUND: Nuclear genes of euglenids contain two major types of introns: conventional spliceosomal and nonconventional introns. The latter are characterized by variable non-canonical borders, RNA secondary structure that brings intron ends together, and an unknown mechanism of removal. Some researchers also distinguish intermediate introns, which combine features of both types. They form a stable RNA secondary structure and are classified into two subtypes depending on whether they contain one (intermediate/nonconventional subtype) or both (conventional/intermediate subtype) canonical spliceosomal borders. However, it has been also postulated that most introns classified as intermediate could simply be special cases of conventional or nonconventional introns. RESULTS: Sequences of tubB, hsp90 and gapC genes from six strains of Euglena agilis were obtained. They contain four, six, and two or three introns, respectively (the third intron in the gapC gene is unique for just one strain). Conventional introns were present at three positions: two in the tubB gene (at one position conventional/intermediate introns were also found) and one in the gapC gene. Nonconventional introns are present at ten positions: two in the tubB gene (at one position intermediate/nonconventional introns were also found), six in hsp90 (at four positions intermediate/nonconventional introns were also found), and two in the gapC gene. CONCLUSIONS: Sequence and RNA secondary structure analyses of nonconventional introns confirmed that their most strongly conserved elements are base pairing nucleotides at positions +4, +5 and +6/ -8, -7 and -6 (in most introns CAG/CTG nucleotides were observed). It was also confirmed that the presence of the 5' GT/C end in intermediate/nonconventional introns is not the result of kinship with conventional introns, but is due to evolutionary pressure to preserve the purine at the 5' end. However, an example of a nonconventional intron with GC-AG ends was shown, suggesting the possibility of intron type conversion between nonconventional and conventional. Furthermore, an analysis of conventional introns revealed that the ability to form a stable RNA secondary structure by some introns is probably not a result of their relationship with nonconventional introns. It was also shown that acquisition of new nonconventional introns is an ongoing process and can be observed at the level of a single species. In the recently acquired intron in the gapC gene an extended direct repeats at the intron-exon junctions are present, suggesting that double-strand break repair process could be the source of new nonconventional introns.

Arginine deiminase pathway enzymes: evolutionary history in metamonads and other eukaryotes.

BACKGROUND: Multiple prokaryotic lineages use the arginine deiminase (ADI) pathway for anaerobic energy production by arginine degradation. The distribution of this pathway among eukaryotes has been thought to be very limited, with only two specialized groups living in low oxygen environments (Parabasalia and Diplomonadida) known to possess the complete set of all three enzymes. We have performed an extensive survey of available sequence data in order to map the distribution of these enzymes among eukaryotes and to reconstruct their phylogenies. RESULTS: We have found genes for the complete pathway in almost all examined representatives of Metamonada, the anaerobic protist group that includes parabasalids and diplomonads. Phylogenetic analyses indicate the presence of the complete pathway in the last common ancestor of metamonads and heterologous transformation experiments suggest its cytosolic localization in the metamonad ancestor. Outside Metamonada, the complete pathway occurs rarely, nevertheless, it was found in representatives of most major eukaryotic clades. CONCLUSIONS: Phylogenetic relationships of complete pathways are consistent with the presence of the Archaea-derived ADI pathway in the last common ancestor of all eukaryotes, although other evolutionary scenarios remain possible. The presence of the incomplete set of enzymes is relatively common among eukaryotes and it may be related to the fact that these enzymes are involved in other cellular processes, such as the ornithine-urea cycle. Single protein phylogenies suggest that the evolutionary history of all three enzymes has been shaped by frequent gene losses and horizontal transfers, which may sometimes be connected with their diverse roles in cellular metabolism.

DNA barcoding in autotrophic euglenids: evaluation of COI and 18s rDNA.

Autotrophic euglenids (Euglenophyceae) are a common and abundant group of microbial eukaryotes in freshwater habitats. They have a limited number of features, which can be observed using light microscopy, thus species identification is often problematic. Establishing a barcode for this group is therefore an important step toward the molecular identification of autotrophic euglenids. Based on the literature, we selected verified species and used a plethora of available methods to validate two molecular markers: COI and 18S rDNA (the whole sequence and three fragments separately) as potential DNA barcodes. Analyses of the COI gene were performed based on the data set of 43 sequences (42 obtained in this study) representing 24 species and the COI gene was discarded as a DNA barcode mainly due to a lack of universal primer sites. For 18S rDNA analyses we used a data set containing 263 sequences belonging to 86 taxonomically verified species. We demonstrated that the whole 18S rDNA is too long to be a useful marker, but from the three shorter analyzed variable regions we recommend variable regions V2V3 and V4 of 18S rDNA as autotrophic euglenid barcodes due to their high efficiency (above 95% and 90%, respectively).

Distribution of conventional and nonconventional introns in tubulin (α and ß) genes of euglenids.

The nuclear genomes of euglenids contain three types of introns: conventional spliceosomal introns, nonconventional introns for which a splicing mechanism is unknown (variable noncanonical borders, RNA secondary structure bringing together intron ends), and so-called intermediate introns, which combine features of conventional and nonconventional introns. Analysis of two genes, tubA and tubB, from 20 species of euglenids reveals contrasting distribution patterns of conventional and nonconventional introns--positions of conventional introns are conserved, whereas those of the nonconventional ones are unique to individual species or small groups of closely related taxa. Moreover, in the group of phototrophic euglenids, 11 events of conventional intron loss versus 15 events of nonconventional intron gain were identified. A comparison of all nonconventional intron sequences highlighted the most conserved elements in their sequence and secondary structure. Our results led us to put forward two hypotheses. 1) The first one posits that mutational changes in intron sequence could lead to a change in their excision mechanism--intermediate introns would then be a transitional form between the conventional and nonconventional introns. 2) The second hypothesis concerns the origin of nonconventional introns--because of the presence of inverted repeats near their ends, insertion of MITE-like transposon elements is proposed as a possible source of new introns.

Phylogenetic Relationships and Morphological Character Evolution of Photosynthetic Euglenids (Excavata) Inferred from Taxon-rich Analyses of Five Genes.

Photosynthetic euglenids acquired chloroplasts by secondary endosymbiosis, which resulted in changes to their mode of nutrition and affected the evolution of their morphological characters. Mapping morphological characters onto a reliable molecular tree could elucidate major trends of those changes. We analyzed nucleotide sequence data from regions of three nuclear-encoded genes (nSSU, nLSU, hsp90), one chloroplast-encoded gene (cpSSU) and one nuclear-encoded chloroplast gene (psbO) to estimate phylogenetic relationships among 59 photosynthetic euglenid species. Our results were consistent with previous works; most genera were monophyletic, except for the polyphyletic genus Euglena, and the paraphyletic genus Phacus. We also analyzed character evolution in photosynthetic euglenids using our phylogenetic tree and eight morphological traits commonly used for generic and species diagnoses, including: characters corresponding to well-defined clades, apomorphies like presence of lorica and mucilaginous stalks, and homoplastic characters like rigid cells and presence of large paramylon grains. This research indicated that pyrenoids were lost twice during the evolution of phototrophic euglenids, and that mucocysts, which only occur in the genus Euglena, evolved independently at least twice. In contrast, the evolution of cell shape and chloroplast morphology was difficult to elucidate, and could not be unambiguously reconstructed in our analyses.