Morphological Identification and Single-Cell Genomics of Marine Diplonemids

Description of new diplonemids (Diplonemea, Euglenozoa) and their endosymbionts: Charting the morphological diversity of these poorly known heterotrophic flagellates

Diplonemids are a hyperdiverse group of flagellated protists, but with less than two dozen formally described representatives. Here, we describe four new species of cultured diplonemids, identified on the basis of their 18S rRNA sequences, light-, fluorescence-, scanning- and transmission electron microscopy. Three new species belong to the genus Rhynchopus (R. asiaticus sp.n., R. granulatus sp.n., and R. valaseki sp.n.), while the fourth species is an unusual representative of the genus Lacrimia (L. aflagellata sp.n.). The latter organism is the first diplonemid outside the genus Rhynchopus (as defined previously) to show a gliding trophic stage with flagellar stubs concealed inside the flagellar pocket and a highly motile dispersive swimming stage. Since this character is thus no longer a genus-specific apomorphy, we provide a taxonomic revision of the genus Rhynchopus with separation of the new genus Natarhynchopus gen. n. We also identify bacterial endosymbionts of L. aflagellata and R. asiaticus as Ca. Syngnamydia medusae (Chlamydiales, Simkaniaceae) and Ca. Cytomitobacter rhynchopi sp. n. (Alphaproteobacteria, Holosporaceae), respectively, and discuss their potential functions. This is the first report of a chlamydial symbiont within a diplonemid host. We also propose that diplonemids may serve as vectors for chlamydial pathogens of marine fish.

Comprehensive analysis of the Kinetoplastea intron landscape reveals a novel intron-containing gene and the first exclusively trans-splicing eukaryote

BACKGROUND

In trypanosomatids, a group of unicellular eukaryotes that includes numerous important human parasites, cis-splicing has been previously reported for only two genes: a poly(A) polymerase and an RNA helicase. Conversely, trans-splicing, which involves the attachment of a spliced leader sequence, is observed for nearly every protein-coding transcript. So far, our understanding of splicing in this protistan group has stemmed from the analysis of only a few medically relevant species. In this study, we used an extensive dataset encompassing all described trypanosomatid genera to investigate the distribution of intron-containing genes and the evolution of splice sites.

RESULTS

We identified a new conserved intron-containing gene encoding an RNA-binding protein that is universally present in Kinetoplastea. We show that Perkinsela sp., a kinetoplastid endosymbiont of Amoebozoa, represents the first eukaryote completely devoid of cis-splicing, yet still preserving trans-splicing. We also provided evidence for reverse transcriptase-mediated intron loss in Kinetoplastea, extensive conservation of 5' splice sites, and the presence of non-coding RNAs within a subset of retained trypanosomatid introns.

CONCLUSIONS

All three intron-containing genes identified in Kinetoplastea encode RNA-interacting proteins, with a potential to fine-tune the expression of multiple genes, thus challenging the perception of cis-splicing in these protists as a mere evolutionary relic. We suggest that there is a selective pressure to retain cis-splicing in trypanosomatids and that this is likely associated with overall control of mRNA processing. Our study provides new insights into the evolution of introns and, consequently, the regulation of gene expression in eukaryotes.

Phylogenomic position of eupelagonemids, abundant, and diverse deep-ocean heterotrophs

Eupelagonemids, formerly known as Deep Sea Pelagic Diplonemids I (DSPD I), are among the most abundant and diverse heterotrophic protists in the deep ocean, but little else is known about their ecology, evolution, or biology in general. Originally recognized solely as a large clade of environmental ribosomal subunit RNA gene sequences (SSU rRNA), branching with a smaller sister group DSPD II, they were postulated to be diplonemids, a poorly studied branch of Euglenozoa. Although new diplonemids have been cultivated and studied in depth in recent years, the lack of cultured eupelagonemids has limited data to a handful of light micrographs, partial SSU rRNA gene sequences, a small number of genes from single amplified genomes, and only a single formal described species, Eupelagonema oceanica. To determine exactly where this clade goes in the tree of eukaryotes and begin to address the overall absence of biological information about this apparently ecologically important group, we conducted single-cell transcriptomics from two eupelagonemid cells. A SSU rRNA gene phylogeny shows that these two cells represent distinct subclades within eupelagonemids, each different from E. oceanica. Phylogenomic analysis based on a 125-gene matrix contrasts with the findings based on ecological survey data and shows eupelagonemids branch sister to the diplonemid subgroup Hemistasiidae.

Ultrastructure and 3D reconstruction of a diplonemid protist (Diplonemea) and its novel membranous organelle

IMPORTANCE

The knowledge of cell biology of a eukaryotic group is essential for correct interpretation of ecological and molecular data. Although diplonemid protists are one of the most species-rich lineages of marine eukaryotes, only very fragmentary information is available about the cellular architecture of this taxonomically diverse group. Here, a large serial block-face scanning electron microscopy data set complemented with light and fluorescence microscopy allowed the first detailed three-dimensional reconstruction of a diplonemid species. We describe numerous previously unknown peculiarities of the cellular architecture and cell division characteristic for diplonemid flagellates, and illustrate the obtained results with multiple three-dimensional models, comprehensible for non-specialists in protist ultrastructure.

Intragenomic diversity of the V9 hypervariable domain in eukaryotes has little effect on metabarcoding

Metabarcoding revolutionized our understanding of diversity and ecology of microorganisms in different habitats. However, it is also associated with several inherent biases, one of which is associated with intragenomic diversity of a molecular barcode. Here, we compare intragenomic variability of the V9 region of the 18S rRNA gene in 19 eukaryotic phyla abundant in marine plankton. The level of intragenomic variability is comparable across all the phyla, and in most genomes and transcriptomes one V9 sequence and one OTU is predominant. However, most of the variability observed at the barcode level is probably caused by sequencing errors and is mitigated by using a denoising tool, DADA2. The SWARM algorithm commonly used in metabarcoding studies is not optimal for collapsing genuine and erroneous sequences into a single OTU, leading to an overestimation of diversity in metabarcoding data. For an unknown reason, SWARM inflates diversity of eupelagonemids more than that of other eukaryotes.

Gene expression dynamics of natural assemblages of heterotrophic flagellates during bacterivory

BACKGROUND

Marine heterotrophic flagellates (HF) are dominant bacterivores in the ocean, where they represent the trophic link between bacteria and higher trophic levels and participate in the recycling of inorganic nutrients for regenerated primary production. Studying their activity and function in the ecosystem is challenging since most of the HFs in the ocean are still uncultured. In the present work, we investigated gene expression of natural HF communities during bacterivory in four unamended seawater incubations.

RESULTS

The most abundant species growing in our incubations belonged to the taxonomic groups MAST-4, MAST-7, Chrysophyceae, and Telonemia. Gene expression dynamics were similar between incubations and could be divided into three states based on microbial counts, each state displaying distinct expression patterns. The analysis of samples where HF growth was highest revealed some highly expressed genes that could be related to bacterivory. Using available genomic and transcriptomic references, we identified 25 species growing in our incubations and used those to compare the expression levels of these specific genes. Video Abstract CONCLUSIONS: Our results indicate that several peptidases, together with some glycoside hydrolases and glycosyltransferases, are more expressed in phagotrophic than in phototrophic species, and thus could be used to infer the process of bacterivory in natural assemblages.

Recent expansion of metabolic versatility in Diplonema papillatum, the model species of a highly speciose group of marine eukaryotes

BACKGROUND

Diplonemid flagellates are among the most abundant and species-rich of known marine microeukaryotes, colonizing all habitats, depths, and geographic regions of the world ocean. However, little is known about their genomes, biology, and ecological role.

RESULTS

We present the first nuclear genome sequence from a diplonemid, the type species Diplonema papillatum. The ~ 280-Mb genome assembly contains about 32,000 protein-coding genes, likely co-transcribed in groups of up to 100. Gene clusters are separated by long repetitive regions that include numerous transposable elements, which also reside within introns. Analysis of gene-family evolution reveals that the last common diplonemid ancestor underwent considerable metabolic expansion. D. papillatum-specific gains of carbohydrate-degradation capability were apparently acquired via horizontal gene transfer. The predicted breakdown of polysaccharides including pectin and xylan is at odds with reports of peptides being the predominant carbon source of this organism. Secretome analysis together with feeding experiments suggest that D. papillatum is predatory, able to degrade cell walls of live microeukaryotes, macroalgae, and water plants, not only for protoplast feeding but also for metabolizing cell-wall carbohydrates as an energy source. The analysis of environmental barcode samples shows that D. papillatum is confined to temperate coastal waters, presumably acting in bioremediation of eutrophication.

CONCLUSIONS

Nuclear genome information will allow systematic functional and cell-biology studies in D. papillatum. It will also serve as a reference for the highly diverse diplonemids and provide a point of comparison for studying gene complement evolution in the sister group of Kinetoplastida, including human-pathogenic taxa.

Massive Accumulation of Strontium and Barium in Diplonemid Protists

Barium and strontium are often used as proxies of marine productivity in palaeoceanographic reconstructions of global climate. However, long-searched biological drivers for such correlations remain unknown. Here, we report that taxa within one of the most abundant groups of marine planktonic protists, diplonemids (Euglenozoa), are potent accumulators of intracellular barite (BaSO₄), celestite (SrSO₄), and strontiobarite (Ba,Sr)SO₄. In culture, Namystinia karyoxenos accumulates Ba²⁺ and Sr²⁺ 42,000 and 10,000 times higher than the surrounding medium, forming barite and celestite representing 90% of the dry weight, the greatest concentration in biomass known to date. As heterotrophs, diplonemids are not restricted to the photic zone, and they are widespread in the oceans in astonishing abundance and diversity, as their distribution correlates with environmental particulate barite and celestite, prevailing in the mesopelagic zone. We found diplonemid predators, the filter-feeding zooplankton that produces fecal pellets containing the undigested celestite from diplonemids, facilitating its deposition on the seafloor. To the best of our knowledge, evidence for diplonemid biomineralization presents the strongest explanation for the occurrence of particulate barite and celestite in the marine environment. Both structures of the crystals and their variable chemical compositions found in diplonemids fit the properties of environmentally sampled particulate barite and celestite. Finally, we propose that diplonemids, which emerged during the Neoproterozoic era, qualify as impactful players in Ba²⁺/Sr²⁺ cycling in the ocean that has possibly contributed to sedimentary rock formation over long geological periods. IMPORTANCE We have identified that diplonemids, an abundant group of marine planktonic protists, accumulate conspicuous amounts of Sr²⁺ and Ba²⁺ in the form of intracellular barite and celestite crystals, in concentrations that greatly exceed those of the most efficient Ba/Sr-accumulating organisms known to date. We propose that diplonemids are potential players in Ba²⁺/Sr²⁺ cycling in the ocean and have possibly contributed to sedimentary rock formation over long geological periods. These organisms emerged during the Neoproterozoic era (590 to 900 million years ago), prior to known coccolithophore carbonate biomineralization (~200 million years ago). Based on reported data, the distribution of diplonemids in the oceans is correlated with the occurrence of particulate barite and celestite. Finally, diplonemids may provide new insights into the long-questioned biogenic origin of particulate barite and celestite and bring more understanding of the observed spatial-temporal correlation of the minerals with marine productivity used in reconstructions of past global climate.

Eukaryotic evolution: Deep phylogeny does not imply morphological novelty

Eukaryotic diversity is often depicted as a molecular phylogenetic tree consisting of a few supergroups that originated over a billion years ago. A new study reveals an ancient group of tiny phagotrophic flagellates that reinforces inferences about early evolutionary history.

To Die or Not to Die-Regulated Cell Death and Survival in Cyanobacteria

Regulated cell death (RCD) is central to the development, integrity, and functionality of multicellular organisms. In the last decade, evidence has accumulated that RCD is a universal phenomenon in all life domains. Cyanobacteria are of specific interest due to their importance in aquatic and terrestrial habitats and their role as primary producers in global nutrient cycling. Current knowledge on cyanobacterial RCD is based mainly on biochemical and morphological observations, often by methods directly transferred from vertebrate research and with limited understanding of the molecular genetic basis. However, the metabolism of different cyanobacteria groups relies on photosynthesis and nitrogen fixation, whereas mitochondria are the central executioner of cell death in vertebrates. Moreover, cyanobacteria chosen as biological models in RCD studies are mainly colonial or filamentous multicellular organisms. On the other hand, unicellular cyanobacteria have regulated programs of cellular survival (RCS) such as chlorosis and post-chlorosis resuscitation. The co-existence of different genetically regulated programs in cyanobacterial populations may have been a top engine in life diversification. Development of cyanobacteria-specific methods for identification and characterization of RCD and wider use of single-cell analysis combined with intelligent image-based cell sorting and metagenomics would shed more light on the underlying molecular mechanisms and help us to address the complex colonial interactions during these events. In this review, we focus on the functional implications of RCD in cyanobacterial communities.

Typical structure of rRNA coding genes in diplonemids points to two independent origins of the bizarre rDNA structures of euglenozoans

Background

Members of Euglenozoa (Discoba) are known for unorthodox rDNA organization. In Euglenida rDNA is located on extrachromosomal circular DNA. In Kinetoplastea and Euglenida the core of the large ribosomal subunit, typically formed by the 28S rRNA, consists of several smaller rRNAs. They are the result of the presence of additional internal transcribed spacers (ITSs) in the rDNA. Diplonemea is the third of the main groups of Euglenozoa and its members are known to be among the most abundant and diverse protists in the oceans. Despite that, the rRNA of only one diplonemid species, Diplonema papillatum, has been examined so far and found to exhibit continuous 28S rRNA. Currently, the rDNA organization has not been researched for any diplonemid. Herein we investigate the structure of rRNA genes in classical (Diplonemidae) and deep-sea diplonemids (Eupelagonemidae), representing the majority of known diplonemid diversity. The results fill the gap in knowledge about diplonemid rDNA and allow better understanding of the evolution of the fragmented structure of the rDNA in Euglenozoa.

Results

We used available genomic (culture and single-cell) sequencing data to assemble complete or almost complete rRNA operons for three classical and six deep-sea diplonemids. The rDNA sequences acquired for several euglenids and kinetoplastids were used to provide the background for the analysis. In all nine diplonemids, 28S rRNA seems to be contiguous, with no additional ITSs detected. Similarly, no additional ITSs were detected in basal prokinetoplastids. However, we identified five additional ITSs in the 28S rRNA of all analysed metakinetoplastids, and up to twelve in euglenids. Only three of these share positions, and they cannot be traced back to their common ancestor.

Conclusions

Presented results indicate that independent origin of additional ITSs in euglenids and kinetoplastids seems to be the most likely. The reason for such unmatched fragmentation remains unknown, but for some reason euglenozoan ribosomes appear to be prone to 28S rRNA fragmentation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02014-9.

Trophic flexibility of marine diplonemids - switching from osmotrophy to bacterivory

Diplonemids are one of the most abundant groups of heterotrophic planktonic microeukaryotes in the world ocean and, thus, are likely to play an essential role in marine ecosystems. So far, only few species have been introduced into a culture, allowing basic studies of diplonemid genetics, morphology, ultrastructure, metabolism, as well as endosymbionts. However, it remains unclear whether these heterotrophic flagellates are parasitic or free-living and what are their predominant dietary patterns and preferred food items. Here we show that cultured diplonemids, maintained in an organic-rich medium as osmotrophs, can gradually switch to bacterivory as a sole food resource, supporting positive growth of their population, even when fed with a low biovolume of bacteria. We further observed remarkable differences in species-specific feeding patterns, size-selective grazing preferences, and distinct feeding strategies. Diplonemids can discriminate between low-quality food items and inedible particles, such as latex beads, even after their ingestion, by discharging them in the form of large waste vacuoles. We also detected digestion-related endogenous autofluorescence emitted by lysosomes and the activity of a melanin-like material. We present the first evidence that these omnipresent protists possess an opportunistic lifestyle that provides a considerable advantage in the generally food resource-limited marine environments.

Highly flexible metabolism of the marine euglenozoan protist Diplonema papillatum

BACKGROUND

The phylum Euglenozoa is a group of flagellated protists comprising the diplonemids, euglenids, symbiontids, and kinetoplastids. The diplonemids are highly abundant and speciose, and recent tools have rendered the best studied representative, Diplonema papillatum, genetically tractable. However, despite the high diversity of diplonemids, their lifestyles, ecological functions, and even primary energy source are mostly unknown.

RESULTS

We designed a metabolic map of D. papillatum cellular bioenergetic pathways based on the alterations of transcriptomic, proteomic, and metabolomic profiles obtained from cells grown under different conditions. Comparative analysis in the nutrient-rich and nutrient-poor media, as well as the absence and presence of oxygen, revealed its capacity for extensive metabolic reprogramming that occurs predominantly on the proteomic rather than the transcriptomic level. D. papillatum is equipped with fundamental metabolic routes such as glycolysis, gluconeogenesis, TCA cycle, pentose phosphate pathway, respiratory complexes, β-oxidation, and synthesis of fatty acids. Gluconeogenesis is uniquely dominant over glycolysis under all surveyed conditions, while the TCA cycle represents an eclectic combination of standard and unusual enzymes.

CONCLUSIONS

The identification of conventional anaerobic enzymes reflects the ability of this protist to survive in low-oxygen environments. Furthermore, its metabolism quickly reacts to restricted carbon availability, suggesting a high metabolic flexibility of diplonemids, which is further reflected in cell morphology and motility, correlating well with their extreme ecological valence.

A cell-cell atlas approach for understanding symbiotic interactions between microbes

Natural environments are composed of a huge diversity of microorganisms interacting with each other to form complex functional networks. Our understanding of the operative nature of host-symbiont associations is limited because propagating such associations in a laboratory is challenging. The advent of single-cell technologies applied to, for example, animal cells and apicomplexan parasites has revolutionized our understanding of development and disease. Such cell atlas approaches generate maps of cell-specific processes and variations within cellular populations. These methods can now be combined with cellular-imaging so that interaction stage versus transcriptome state can be quantized for microbe-microbe interactions. We predict that the combination of these methods applied to the study of symbioses will transform our understanding of many ecological interactions, including those sampled directly from natural environments.

Single-Cell Transcriptomics of Abedinium Reveals a New Early-Branching Dinoflagellate Lineage

Dinoflagellates possess many cellular characteristics with unresolved evolutionary histories. These include nuclei with greatly expanded genomes and chromatin packaged using histone-like proteins and dinoflagellate-viral nucleoproteins instead of histones, highly reduced mitochondrial genomes with extensive RNA editing, a mix of photosynthetic and cryptic secondary plastids, and tertiary plastids. Resolving the evolutionary origin of these traits requires understanding their ancestral states and early intermediates. Several early-branching dinoflagellate lineages are good candidates for such reconstruction, however these cells tend to be delicate and environmentally sparse, complicating such analyses. Here, we employ transcriptome sequencing from manually isolated and microscopically documented cells to resolve the placement of two cells of one such genus, Abedinium, collected by remotely operated vehicle in deep waters off the coast of Monterey Bay, CA. One cell corresponds to the only described species, Abedinium dasypus, whereas the second cell is distinct and formally described as Abedinium folium, sp. nov. Abedinium has classically been assigned to the early-branching dinoflagellate subgroup Noctilucales, which is weakly supported by phylogenetic analyses of small subunit ribosomal RNA, the single characterized gene from any member of the order. However, an analysis based on 221 proteins from the transcriptome places Abedinium as a distinct lineage, separate from and basal to Noctilucales and the rest of the core dinoflagellates. The transcriptome also contains evidence of a cryptic plastid functioning in the biosynthesis of isoprenoids, iron-sulfur clusters, and heme, a mitochondrial genome with all three expected protein-coding genes (cob, cox1, and cox3), and the presence of some but not all dinoflagellate-specific chromatin packaging proteins.

Heterotrophic euglenid Rhabdomonas costata resembles its phototrophic relatives in many aspects of molecular and cell biology

Euglenids represent a group of protists with diverse modes of feeding. To date, only a partial genomic sequence of Euglena gracilis and transcriptomes of several phototrophic and secondarily osmotrophic species are available, while primarily heterotrophic euglenids are seriously undersampled. In this work, we begin to fill this gap by presenting genomic and transcriptomic drafts of a primary osmotroph, Rhabdomonas costata. The current genomic assembly length of 100 Mbp is 14× smaller than that of E. gracilis. Despite being too fragmented for comprehensive gene prediction it provided fragments of the mitochondrial genome and comparison of the transcriptomic and genomic data revealed features of its introns, including several candidates for nonconventional types. A set of 39,456 putative R. costata proteins was predicted from the transcriptome. Annotation of the mitochondrial core metabolism provides the first data on the facultatively anaerobic mitochondrion of R. costata, which in most respects resembles the mitochondrion of E. gracilis with a certain level of streamlining. R. costata can synthetise thiamine by enzymes of heterogenous provenances and haem by a mitochondrial-cytoplasmic C4 pathway with enzymes orthologous to those found in E. gracilis. The low percentage of green algae-affiliated genes supports the ancestrally osmotrophic status of this species.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01035-y.

A single-cell genomics pipeline for environmental microbial eukaryotes

Single-cell sequencing of environmental microorganisms is an essential component of the microbial ecology toolkit. However, large-scale targeted single-cell sequencing for the whole-genome recovery of uncultivated eukaryotes is lagging. The key challenges are low abundance in environmental communities, large complex genomes, and cell walls that are difficult to break. We describe a pipeline composed of state-of-the art single-cell genomics tools and protocols optimized for poorly studied and uncultivated eukaryotic microorganisms that are found at low abundance. This pipeline consists of seven distinct steps, beginning with sample collection and ending with genome annotation, each equipped with quality review steps to ensure high genome quality at low cost. We tested and evaluated each step on environmental samples and cultures of early-diverging lineages of fungi and Chromista/SAR. We show that genomes produced using this pipeline are almost as good as complete reference genomes for functional and comparative genomics for environmental microbial eukaryotes.

Euglenozoa: taxonomy, diversity and ecology, symbioses and viruses

Euglenozoa is a species-rich group of protists, which have extremely diverse lifestyles and a range of features that distinguish them from other eukaryotes. They are composed of free-living and parasitic kinetoplastids, mostly free-living diplonemids, heterotrophic and photosynthetic euglenids, as well as deep-sea symbiontids. Although they form a well-supported monophyletic group, these morphologically rather distinct groups are almost never treated together in a comparative manner, as attempted here. We present an updated taxonomy, complemented by photos of representative species, with notes on diversity, distribution and biology of euglenozoans. For kinetoplastids, we propose a significantly modified taxonomy that reflects the latest findings. Finally, we summarize what is known about viruses infecting euglenozoans, as well as their relationships with ecto- and endosymbiotic bacteria.

First finding of free-living representatives of Prokinetoplastina and their nuclear and mitochondrial genomes

Kinetoplastids are heterotrophic flagellated protists, including important parasites of humans and animals (trypanosomatids), and ecologically important free-living bacterial consumers (bodonids). Phylogenies have shown that the earliest-branching kinetoplastids are all parasites or obligate endosymbionts, whose highly-derived state makes reconstructing the ancestral state of the group challenging. We have isolated new strains of unusual free-living flagellates that molecular phylogeny shows to be most closely related to endosymbiotic and parasitic Perkinsela and Ichthyobodo species that, together with unidentified environmental sequences, form the clade at the base of kinetoplastids. These strains are therefore the first described free-living prokinetoplastids, and potentially very informative in understanding the evolution and ancestral states of morphological and molecular characteristics described in other kinetoplastids. Overall, we find that these organisms morphologically and ultrastructurally resemble some free-living bodonids and diplonemids, and possess nuclear genomes with few introns, polycistronic mRNA expression, high coding density, and derived traits shared with other kinetoplastids. Their genetic repertoires are more diverse than the best-studied free-living kinetoplastids, which is likely a reflection of their higher metabolic potential. Mitochondrial RNAs of these new species undergo the most extensive U insertion/deletion editing reported so far, and limited deaminative C-to-U and A-to-I editing, but we find no evidence for mitochondrial trans-splicing.

Structural Comparison of Diplonemid Communities around the Izu Peninsula, Japan

Diplonemea (diplonemids) is one of the most abundant and species-rich protist groups in marine environments; however, their community structures among local and seasonal samples have not yet been compared. In the present study, we analyzed four diplonemid community structures around the Izu Peninsula, Japan using barcode sequences amplified from environmental DNA. These sequences and the results of statistical analyses indicated that communities at the same site were more similar to each other than those in the same season. Environmental variables were also measured, and their influence on diplonemid community structures was examined. Salinity, electrical conductivity, and temperature, and their correlated variables, appeared to influence the structures of diplonemid communities, which was consistent with previous findings; however, since the results obtained did not reach statistical significance, further studies are required. A comparison of each diplonemid community indicated that some lineages were unique to specific samples, while others were consistently detected in all samples. Members of the latter type are cosmopolitan candidates and may be better adapted to the environments of the studied area. Future studies that focus on the more adaptive members will provide a more detailed understanding of the mechanisms by which diplonemids are widely distributed in marine environments and will facilitate their utilization as indicator organisms to monitor environmental changes.

A New Type of Circular RNA derived from Nonconventional Introns in Nuclear Genes of Euglenids

Nuclear protein-coding genes of euglenids (Discoba, Euglenozoa, Euglenida) contain conventional (spliceosomal) and nonconventional introns. The latter have been found only in euglenozoans. A unique feature of nonconventional introns is the ability to form a stable and slightly conserved RNA secondary structure bringing together intron ends and placing adjacent exons in proximity. To date, little is known about the mechanism of their excision (e.g. whether it involves the spliceosome or not). The tubA gene of Euglena gracilis harbors three conventional and three nonconventional introns. While the conventional introns are excised as lariats, nonconventional introns are present in the cell solely as circular RNAs with full-length ends. Based on this discovery as well as on previous observations indicating that nonconventional introns are observed frequently at unique positions of genes, we suggest that this new type of intronic circRNA might play a role in intron mobility.

A freshwater radiation of diplonemids

Diplonemids are considered marine protists and have been reported among the most abundant and diverse eukaryotes in the world oceans. Recently we detected the presence of freshwater diplonemids in Japanese deep freshwater lakes. However, their distribution and abundances in freshwater ecosystems remain unknown. We assessed abundance and diversity of diplonemids from several geographically distant deep freshwater lakes of the world by amplicon-sequencing, shotgun metagenomics and catalysed reporter deposition-fluorescent in situ hybridization (CARD-FISH). We found diplonemids in all the studied lakes, albeit with low abundances and diversity. We assembled long 18S rRNA sequences from freshwater diplonemids and showed that they form a new lineage distinct from the diverse marine clades. Freshwater diplonemids are a sister-group to a marine clade, which are mainly isolates from coastal and bay areas, suggesting a recent habitat transition from marine to freshwater habitats. Images of CARD-FISH targeted freshwater diplonemids suggest they feed on bacteria. Our analyses of 18S rRNA sequences retrieved from single-cell genomes of marine diplonemids show they encode multiple rRNA copies that may be very divergent from each other, suggesting that marine diplonemid abundance and diversity both have been overestimated. These results have wider implications on assessing eukaryotic abundances in natural habitats by using amplicon-sequencing alone.

Environmental determinants of the distribution of planktonic diplonemids and kinetoplastids in the oceans

We analysed a widely used barcode, the V9 region of the 18S rRNA gene, to study the effect of environmental conditions on the distribution of two related heterotrophic protistan lineages in marine plankton, kinetoplastids and diplonemids. We relied on a major published dataset (Tara Oceans) where samples from the mesopelagic zone were available from just 32 of 123 locations, and both groups are most abundant in this zone. To close sampling gaps and obtain more information from the deeper ocean, we collected 57 new samples targeting especially the mesopelagic zone. We sampled in three geographic regions: the Arctic, two depth transects in the Adriatic Sea, and the anoxic Cariaco Basin. In agreement with previous studies, both protist groups are most abundant and diverse in the mesopelagic zone. In addition to that, we found that their abundance, richness, and community structure also depend on geography, oxygen concentration, salinity, temperature, and other environmental variables reflecting the abundance of algae and nutrients. Both groups studied here demonstrated similar patterns, although some differences were also observed. Kinetoplastids and diplonemids prefer tropical regions and nutrient-rich conditions and avoid high oxygen concentration, high salinity, and high density of algae.

Targeted integration by homologous recombination enables in situ tagging and replacement of genes in the marine microeukaryote Diplonema papillatum

Diplonemids are a group of highly diverse and abundant marine microeukaryotes that belong to the phylum Euglenozoa and form a sister clade to the well-studied, mostly parasitic kinetoplastids. Very little is known about the biology of diplonemids, as few species have been formally described and just one, Diplonema papillatum, has been studied to a decent extent at the molecular level. Following up on our previous results showing stable but random integration of delivered extraneous DNA, we demonstrate here homologous recombination in D. papillatum. Targeting various constructs to the intended position in the nuclear genome was successful when 5' and 3' homologous regions longer than 1 kbp were used, achieving N-terminal tagging with mCherry and gene replacement of α- and β-tubulins. For more convenient genetic manipulation, we designed a modular plasmid, pDP002, which bears a protein-A tag and used it to generate and express a C-terminally tagged mitoribosomal protein. Lastly, we developed an improved transformation protocol for broader applicability across laboratories. Our robust methodology allows the replacement, integration as well as endogenous tagging of D. papillatum genes, thus opening the door to functional studies in this species and establishing a basic toolkit for reverse genetics of diplonemids in general.

An efficient single-cell transcriptomics workflow for microbial eukaryotes benchmarked on Giardia intestinalis cells

Background

Most diversity in the eukaryotic tree of life is represented by microbial eukaryotes, which is a polyphyletic group also referred to as protists. Among the protists, currently sequenced genomes and transcriptomes give a biased view of the actual diversity. This biased view is partly caused by the scientific community, which has prioritized certain microbes of biomedical and agricultural importance. Additionally, some protists remain difficult to maintain in cultures, which further influences what has been studied. It is now possible to bypass the time-consuming process of cultivation and directly analyze the gene content of single protist cells. Single-cell genomics was used in the first experiments where individual protists cells were genomically explored. Unfortunately, single-cell genomics for protists is often associated with low genome recovery and the assembly process can be complicated because of repetitive intergenic regions. Sequencing repetitive sequences can be avoided if single-cell transcriptomics is used, which only targets the part of the genome that is transcribed.

Results

In this study we test different modifications of Smart-seq2, a single-cell RNA sequencing protocol originally developed for mammalian cells, to establish a robust and more cost-efficient workflow for protists. The diplomonad Giardia intestinalis was used in all experiments and the available genome for this species allowed us to benchmark our results. We could observe increased transcript recovery when freeze-thaw cycles were added as an extra step to the Smart-seq2 protocol. Further we reduced the reaction volume and purified the amplified cDNA with alternative beads to test different cost-reducing changes of Smart-seq2. Neither improved the procedure, and reducing the volumes by half led to significantly fewer genes detected. We also added a 5′ biotin modification to our primers and reduced the concentration of oligo-dT, to potentially reduce generation of artifacts. Except adding freeze-thaw cycles and reducing the volume, no other modifications lead to a significant change in gene detection. Therefore, we suggest adding freeze-thaw cycles to Smart-seq2 when working with protists and further consider our other modification described to improve cost and time-efficiency.

Conclusions

The presented single-cell RNA sequencing workflow represents an efficient method to explore the diversity and cell biology of individual protist cells.

Protist diversity and function in the dark ocean - Challenging the paradigms of deep-sea ecology with special emphasis on foraminiferans and naked protists

The dark ocean and the underlying deep seafloor together represent the largest environment on this planet, comprising about 80% of the oceanic volume and covering more than two-thirds of the Earth's surface, as well as hosting a major part of the total biosphere. Emerging evidence suggests that these vast pelagic and benthic habitats play a major role in ocean biogeochemistry and represent an "untapped reservoir" of high genetic and metabolic microbial diversity. Due to its huge volume, the water column of the dark ocean is the largest reservoir of organic carbon in the biosphere and likely plays a major role in the global carbon budget. The dark ocean and the seafloor beneath it are also home to a largely enigmatic food web comprising little-known and sometimes spectacular organisms, mainly prokaryotes and protists. This review considers the globally important role of pelagic and benthic protists across all protistan size classes in the deep-sea realm, with a focus on their taxonomy, diversity, and physiological properties, including their role in deep microbial food webs. We argue that, given the important contribution that protists must make to deep-sea biodiversity and ecosystem processes, they should not be overlooked in biological studies of the deep ocean.

Tara Oceans: towards global ocean ecosystems biology

A planetary-scale understanding of the ocean ecosystem, particularly in light of climate change, is crucial. Here, we review the work of Tara Oceans, an international, multidisciplinary project to assess the complexity of ocean life across comprehensive taxonomic and spatial scales. Using a modified sailing boat, the team sampled plankton at 210 globally distributed sites at depths down to 1,000 m. We describe publicly available resources of molecular, morphological and environmental data, and discuss how an ecosystems biology approach has expanded our understanding of plankton diversity and ecology in the ocean as a planetary, interconnected ecosystem. These efforts illustrate how global-scale concepts and data can help to integrate biological complexity into models and serve as a baseline for assessing ecosystem changes and the future habitability of our planet in the Anthropocene epoch.

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.

SCGid: a consensus approach to contig filtering and genome prediction from single-cell sequencing libraries of uncultured eukaryotes

MOTIVATION

Whole-genome sequencing of uncultured eukaryotic genomes is complicated by difficulties in acquiring sufficient amounts of tissue. Single-cell genomics (SCG) by multiple displacement amplification provides a technical workaround, yielding whole-genome libraries which can be assembled de novo. Downsides of multiple displacement amplification include coverage biases and exacerbation of contamination. These factors affect assembly continuity and fidelity, complicating discrimination of genomes from contamination and noise by available tools. Uncultured eukaryotes and their relatives are often underrepresented in large sequence data repositories, further impairing identification and separation.

RESULTS

We compare the ability of filtering approaches to remove contamination and resolve eukaryotic draft genomes from SCG metagenomes, finding significant variation in outcomes. To address these inconsistencies, we introduce a consensus approach that is codified in the SCGid software package. SCGid parallelly filters assemblies using different approaches, yielding three intermediate drafts from which consensus is drawn. Using genuine and mock SCG metagenomes, we show that our approach corrects for variation among draft genomes predicted by individual approaches and outperforms them in recapitulating published drafts in a fast and repeatable way, providing a useful alternative to available methods and manual curation.

AVAILABILITY AND IMPLEMENTATION

The SCGid package is implemented in python and R. Source code is available at http://www.github.com/amsesk/SCGid under the GNU GPL 3.0 license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Gene fragmentation and RNA editing without borders: eccentric mitochondrial genomes of diplonemids

Diplonemids are highly abundant heterotrophic marine protists. Previous studies showed that their strikingly bloated mitochondrial genome is unique because of systematic gene fragmentation and manifold RNA editing. Here we report a comparative study of mitochondrial genome architecture, gene structure and RNA editing of six recently isolated, phylogenetically diverse diplonemid species. Mitochondrial gene fragmentation and modes of RNA editing, which include cytidine-to-uridine (C-to-U) and adenosine-to-inosine (A-to-I) substitutions and 3' uridine additions (U-appendage), are conserved across diplonemids. Yet as we show here, all these features have been pushed to their extremes in the Hemistasiidae lineage. For example, Namystynia karyoxenos has its genes fragmented into more than twice as many modules than other diplonemids, with modules as short as four nucleotides. Furthermore, we detected in this group multiple A-appendage and guanosine-to-adenosine (G-to-A) substitution editing events not observed before in diplonemids and found very rarely elsewhere. With >1,000 sites, C-to-U and A-to-I editing in Namystynia is nearly 10 times more frequent than in other diplonemids. The editing density of 12% in coding regions makes Namystynia's the most extensively edited transcriptome described so far. Diplonemid mitochondrial genome architecture, gene structure and post-transcriptional processes display such high complexity that they challenge all other currently known systems.

A metagenomic assessment of microbial eukaryotic diversity in the global ocean

Surveying microbial diversity and function is accomplished by combining complementary molecular tools. Among them, metagenomics is a PCR free approach that contains all genetic information from microbial assemblages and is today performed at a relatively large scale and reasonable cost, mostly based on very short reads. Here, we investigated the potential of metagenomics to provide taxonomic reports of marine microbial eukaryotes. We prepared a curated database with reference sequences of the V4 region of 18S rDNA clustered at 97% similarity and used this database to extract and classify metagenomic reads. More than half of them were unambiguously affiliated to a unique reference whilst the rest could be assigned to a given taxonomic group. The overall diversity reported by metagenomics was similar to that obtained by amplicon sequencing of the V4 and V9 regions of the 18S rRNA gene, although either one or both of these amplicon surveys performed poorly for groups like Excavata, Amoebozoa, Fungi and Haptophyta. We then studied the diversity of picoeukaryotes and nanoeukaryotes using 91 metagenomes from surface down to bathypelagic layers in different oceans, unveiling a clear taxonomic separation between size fractions and depth layers. Finally, we retrieved long rDNA sequences from assembled metagenomes that improved phylogenetic reconstructions of particular groups. Overall, this study shows metagenomics as an excellent resource for taxonomic exploration of marine microbial eukaryotes.

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.

Highly Reduced Genomes of Protist Endosymbionts Show Evolutionary Convergence

Genome evolution in bacterial endosymbionts is notoriously extreme: the combined effects of strong genetic drift and unique selective pressures result in highly reduced genomes with distinctive adaptations to hosts [1-4]. These processes are mostly known from animal endosymbionts, where nutritional endosymbioses represent the best-studied systems. However, eukaryotic microbes, or protists, also harbor diverse bacterial endosymbionts, but their genome reduction and functional relationships with their hosts are largely unexplored [5-7]. We sequenced the genomes of four bacterial endosymbionts from three species of diplonemids, poorly studied but abundant and diverse heterotrophic protists [8-12]. The endosymbionts come from two bacterial families, Rickettsiaceae and Holosporaceae, that have invaded two families of diplonemids, and their genomes have converged on an extremely small size (605-632 kilobase pairs [kbp]), similar gene content (e.g., metabolite transporters and secretion systems), and reduced metabolic potential (e.g., loss of energy metabolism). These characteristics are generally found in both families, but the diplonemid endosymbionts have evolved greater extremes in parallel. They possess modified type VI secretion systems that could function in manipulating host metabolism or other intracellular interactions. Finally, modified cellular machinery like the ATP synthase without oxidative phosphorylation, and the reduced flagellar apparatus present in some diplonemid endosymbionts and nutritional animal endosymbionts, indicates that intracellular mechanisms have converged in bacterial endosymbionts with various functions and from different eukaryotic hosts across the tree of life.

EukRef-excavates: seven curated SSU ribosomal RNA gene databases

The small subunit ribosomal RNA (SSU rRNA) gene is a widely used molecular marker to study the diversity of life. Sequencing of SSU rRNA gene amplicons has become a standard approach for the investigation of the ecology and diversity of microbes. However, a well-curated database is necessary for correct classification of these data. While available for many groups of Bacteria and Archaea, such reference databases are absent for most eukaryotes. The primary goal of the EukRef project (eukref.org) is to close this gap and generate well-curated reference databases for major groups of eukaryotes, especially protists. Here we present a set of EukRef-curated databases for the excavate protists—a large assemblage that includes numerous taxa with divergent SSU rRNA gene sequences, which are prone to misclassification. We identified 6121 sequences, 625 of which were obtained from cultures, 3053 from cell isolations or enrichments and 2419 from environmental samples. We have corrected the classification for the majority of these curated sequences. The resulting publicly available databases will provide phylogenetically based standards for the improved identification of excavates in ecological and microbiome studies, as well as resources to classify new discoveries in excavate diversity.

The need to account for cell biology in characterizing predatory mixotrophs in aquatic environments

Photosynthesis in eukaryotes first arose through phagocytotic processes wherein an engulfed cyanobacterium was not digested, but instead became a permanent organelle. Other photosynthetic lineages then arose when eukaryotic cells engulfed other already photosynthetic eukaryotic cells. Some of the resulting lineages subsequently lost their ability for phagocytosis, while many others maintained the ability to do both processes. These mixotrophic taxa have more complicated ecological roles, in that they are both primary producers and consumers that can shift more towards producing the organic matter that forms the base of aquatic food chains, or towards respiring and releasing CO2. We still have much to learn about which taxa are predatory mixotrophs as well as about the physiological consequences of this lifestyle, in part, because much of the diversity of unicellular eukaryotes in aquatic ecosystems remains uncultured. Here, we discuss existing methods for studying predatory mixotrophs, their individual biases, and how single-cell approaches can enhance knowledge of these important taxa. The question remains what the gold standard should be for assigning a mixotrophic status to ill-characterized or uncultured taxa-a status that dictates how organisms are incorporated into carbon cycle models and how their ecosystem roles may shift in future lakes and oceans. This article is part of a discussion meeting issue 'Single cell ecology'.

Unexpected mitochondrial genome diversity revealed by targeted single-cell genomics of heterotrophic flagellated protists

Most eukaryotic microbial diversity is uncultivated, under-studied and lacks nuclear genome data. Mitochondrial genome sampling is more comprehensive, but many phylogenetically important groups remain unsampled. Here, using a single-cell sorting approach combining tubulin-specific labelling with photopigment exclusion, we sorted flagellated heterotrophic unicellular eukaryotes from Pacific Ocean samples. We recovered 206 single amplified genomes, predominantly from underrepresented branches on the tree of life. Seventy single amplified genomes contained unique mitochondrial contigs, including 21 complete or near-complete mitochondrial genomes from formerly under-sampled phylogenetic branches, including telonemids, katablepharids, cercozoans and marine stramenopiles, effectively doubling the number of available samples of heterotrophic flagellate mitochondrial genomes. Collectively, these data identify a dynamic history of mitochondrial genome evolution including intron gain and loss, extensive patterns of genetic code variation and complex patterns of gene loss. Surprisingly, we found that stramenopile mitochondrial content is highly plastic, resembling patterns of variation previously observed only in plants.

Using single-cell transcriptomics to understand functional states and interactions in microbial eukaryotes

Understanding the diversity and evolution of eukaryotic microorganisms remains one of the major challenges of modern biology. In recent years, we have advanced in the discovery and phylogenetic placement of new eukaryotic species and lineages, which in turn completely transformed our view on the eukaryotic tree of life. But we remain ignorant of the life cycles, physiology and cellular states of most of these microbial eukaryotes, as well as of their interactions with other organisms. Here, we discuss how high-throughput genome-wide gene expression analysis of eukaryotic single cells can shed light on protist biology. First, we review different single-cell transcriptomics methodologies with particular focus on microbial eukaryote applications. Then, we discuss single-cell gene expression analysis of protists in culture and what can be learnt from these approaches. Finally, we envision the application of single-cell transcriptomics to protist communities to interrogate not only community components, but also the gene expression signatures of distinct cellular and physiological states, as well as the transcriptional dynamics of interspecific interactions. Overall, we argue that single-cell transcriptomics can significantly contribute to our understanding of the biology of microbial eukaryotes. This article is part of a discussion meeting issue 'Single cell ecology'.

The New Tree of Eukaryotes

For 15 years, the eukaryote Tree of Life (eToL) has been divided into five to eight major groupings, known as 'supergroups'. However, the tree has been profoundly rearranged during this time. The new eToL results from the widespread application of phylogenomics and numerous discoveries of major lineages of eukaryotes, mostly free-living heterotrophic protists. The evidence that supports the tree has transitioned from a synthesis of molecular phylogenetics and biological characters to purely molecular phylogenetics. Most current supergroups lack defining morphological or cell-biological characteristics, making the supergroup label even more arbitrary than before. Going forward, the combination of traditional culturing with maturing culture-free approaches and phylogenomics should accelerate the process of completing and resolving the eToL at its deepest levels.

A single-cell genome reveals diplonemid-like ancestry of kinetoplastid mitochondrial gene structure

Euglenozoa comprises euglenids, kinetoplastids, and diplonemids, with each group exhibiting different and highly unusual mitochondrial genome organizations. Although they are sister groups, kinetoplastids and diplonemids have very distinct mitochondrial genome architectures, requiring widespread insertion/deletion RNA editing and extensive trans-splicing, respectively, in order to generate functional transcripts. The evolutionary history by which these differing processes arose remains unclear. Using single-cell genomics, followed by small sub unit ribosomal DNA and multigene phylogenies, we identified an isolated marine cell that branches on phylogenetic trees as a sister to known kinetoplastids. Analysis of single-cell amplified genomic material identified multiple mitochondrial genome contigs. These revealed a gene architecture resembling that of diplonemid mitochondria, with small fragments of genes encoded out of order and or on different contigs, indicating that these genes require extensive trans-splicing. Conversely, no requirement for kinetoplastid-like insertion/deletion RNA-editing was detected. Additionally, while we identified some proteins so far only found in kinetoplastids, we could not unequivocally identify mitochondrial RNA editing proteins. These data invite the hypothesis that extensive genome fragmentation and trans-splicing were the ancestral states for the kinetoplastid-diplonemid clade but were lost during the kinetoplastid radiation. This study demonstrates that single-cell approaches can successfully retrieve lineages that represent important new branches on the tree of life, and thus can illuminate major evolutionary and functional transitions in eukaryotes.

This article is part of a discussion meeting issue ‘Single cell ecology’.

A single-cell genome perspective on studying intracellular associations in unicellular eukaryotes

Single-cell genomics (SCG) methods provide a unique opportunity to analyse whole genome information at the resolution of an individual cell. While SCG has been extensively used to investigate bacterial and archaeal genomes, the technique has been rarely used to access the genetic makeup of uncultivated microbial eukaryotes. In this regard, the use of SCG can provide a wealth of information; not only do the methods allow exploration of the genome, they can also help elucidate the relationship between the cell and intracellular entities extant in nearly all eukaryotes. SCG enables the study of total eukaryotic cellular DNA, which in turn allows us to better understand the evolutionary history and diversity of life, and the physiological interactions that define complex organisms.

This article is part of a discussion meeting issue ‘Single cell ecology’.

Single cell ecology

Cells are the building blocks of life, from single-celled microbes through to multi-cellular organisms. To understand a multitude of biological processes we need to understand how cells behave, how they interact with each other and how they respond to their environment. The use of new methodologies is changing the way we study cells allowing us to study them on minute scales and in unprecedented detail. These same methods are allowing researchers to begin to sample the vast diversity of microbes that dominate natural environments. The aim of this special issue is to bring together research and perspectives on the application of new approaches to understand the biological properties of cells, including how they interact with other biological entities. This article is part of a discussion meeting issue 'Single cell ecology'.

Combining morphology, behaviour and genomics to understand the evolution and ecology of microbial eukaryotes

Microbial eukaryotes (protists) are structurally, developmentally and behaviourally more complex than their prokaryotic cousins. This complexity makes it more difficult to translate genomic and metagenomic data into accurate functional inferences about systems ranging all the way from molecular and cellular levels to global ecological networks. This problem can be traced back to the advent of the cytoskeleton and endomembrane systems at the origin of eukaryotes, which endowed them with a range of complex structures and behaviours that still largely dominate how they evolve and interact within microbial communities. But unlike the diverse metabolic properties that evolved within prokaryotes, the structural and behavioural characteristics that strongly define how protists function in the environment cannot readily be inferred from genomic data, since there is generally no simple correlation between a gene and a discrete activity or function. A deeper understanding of protists at both cellular and ecological levels, therefore, requires not only high-throughput genomics but also linking such data to direct observations of natural history and cell biology. This is challenging since these observations typically require cultivation, which is lacking for most protists. Potential remedies with current technology include developing a more phylogenetically diverse range of model systems to better represent the diversity, as well as combining high-throughput, single-cell genomics with microscopic documentation of the subject cells to link sequence with structure and behaviour. This article is part of a discussion meeting issue 'Single cell ecology'.

There Is Treasure Everywhere: Reductive Plastid Evolution in Apicomplexa in Light of Their Close Relatives

The phylum Apicomplexa (Alveolates) comprises a group of host-associated protists, predominately intracellular parasites, including devastating parasites like Plasmodium falciparum, the causative agent of malaria. One of the more fascinating characteristics of Apicomplexa is their highly reduced (and occasionally lost) remnant plastid, termed the apicoplast. Four core metabolic pathways are retained in the apicoplast: heme synthesis, iron–sulfur cluster synthesis, isoprenoid synthesis, and fatty acid synthesis. It has been suggested that one or more of these pathways are essential for plastid and plastid genome retention. The past decade has witnessed the discovery of several apicomplexan relatives, and next-generation sequencing efforts are revealing that they retain variable plastid metabolic capacities. These data are providing clues about the core genes and pathways of reduced plastids, while at the same time further confounding our view on the evolutionary history of the apicoplast. Here, we examine the evolutionary history of the apicoplast, explore plastid metabolism in Apicomplexa and their close relatives, and propose that the differences among reduced plastids result from a game of endosymbiotic roulette. Continued exploration of the Apicomplexa and their relatives is sure to provide new insights into the evolution of the apicoplast and apicomplexans as a whole.

Trypanosomatid mitochondrial RNA editing: dramatically complex transcript repertoires revealed with a dedicated mapping tool

RNA editing by targeted insertion and deletion of uridine is crucial to generate translatable mRNAs from the cryptogenes of the mitochondrial genome of kinetoplastids. This type of editing consists of a stepwise cascade of reactions generally proceeding from 3′ to 5′ on a transcript, resulting in a population of partially edited as well as pre-edited and completely edited molecules for each mitochondrial cryptogene of these protozoans. Often, the number of uridines inserted and deleted exceed the number of nucleotides that are genome-encoded. Thus, analysis of kinetoplastid mitochondrial transcriptomes has proven frustratingly complex. Here we present our analysis of Leptomonas pyrrhocoris mitochondrial cDNA deep sequencing reads using T-Aligner, our new tool which allows comprehensive characterization of RNA editing, not relying on targeted transcript amplification and on prior knowledge of final edited products. T-Aligner implements a pipeline of read mapping, visualization of all editing states and their coverage, and assembly of canonical and alternative translatable mRNAs. We also assess T-Aligner functionality on a more challenging deep sequencing read input from Trypanosoma cruzi. The analysis reveals that transcripts of cryptogenes of both species undergo very complex editing that includes the formation of alternative open reading frames and whole categories of truncated editing products.

Neobodonids are dominant kinetoplastids in the global ocean

Kinetoplastid flagellates comprise basal mostly free-living bodonids and derived obligatory parasitic trypanosomatids, which belong to the best-studied protists. Due to their omnipresence in aquatic environments and soil, the bodonids are of ecological significance. Here, we present the first global survey of marine kinetoplastids and compare it with the strikingly different patterns of abundance and diversity in their sister clade, the diplonemids. Based on analysis of 18S rDNA V9 ribotypes obtained from 124 sites sampled during the Tara Oceans expedition, our results show generally low to moderate abundance and diversity of planktonic kinetoplastids. Although we have identified all major kinetoplastid lineages, 98% of kinetoplastid reads are represented by neobodonids, namely specimens of the Neobodo and Rhynchomonas genera, which make up 59% and 18% of all reads, respectively. Most kinetoplastids have small cell size (0.8-5 µm) and tend to be more abundant in the mesopelagic as compared to the euphotic zone. Some of the most abundant operational taxonomic units have distinct geographical distributions, and three novel putatively parasitic neobodonids were identified, along with their potential hosts.

Comparative molecular cell biology of phototrophic euglenids and parasitic trypanosomatids sheds light on the ancestor of Euglenozoa

Parasitic trypanosomatids and phototrophic euglenids are among the most extensively studied euglenozoans. The phototrophic euglenid lineage arose relatively recently through secondary endosymbiosis between a phagotrophic euglenid and a prasinophyte green alga that evolved into the euglenid secondary chloroplast. The parasitic trypanosomatids (i.e. Trypanosoma spp. and Leishmania spp.) and the freshwater phototrophic euglenids (i.e. Euglena gracilis) are the most evolutionary distant lineages in the Euglenozoa phylogenetic tree. The molecular and cell biological traits they share can thus be considered as ancestral traits originating in the common euglenozoan ancestor. These euglenozoan ancestral traits include common mitochondrial presequence motifs, respiratory chain complexes containing various unique subunits, a unique ATP synthase structure, the absence of mitochondria-encoded transfer RNAs (tRNAs), a nucleus with a centrally positioned nucleolus, closed mitosis without dissolution of the nuclear membrane and nucleoli, a nuclear genome containing the unusual 'J' base (β-D-glucosyl-hydroxymethyluracil), processing of nucleus-encoded precursor messenger RNAs (pre-mRNAs) via spliced-leader RNA (SL-RNA) trans-splicing, post-transcriptional gene silencing by the RNA interference (RNAi) pathway and the absence of transcriptional regulation of nuclear gene expression. Mitochondrial uridine insertion/deletion RNA editing directed by guide RNAs (gRNAs) evolved in the ancestor of the kinetoplastid lineage. The evolutionary origin of other molecular features known to be present only in either kinetoplastids (i.e. polycistronic transcripts, compaction of nuclear genomes) or euglenids (i.e. monocistronic transcripts, huge genomes, many nuclear cis-spliced introns, polyproteins) is unclear.

Morphological, Ultrastructural, Motility and Evolutionary Characterization of Two New Hemistasiidae Species

Until now, Hemistasia phaeocysticola was the only representative of the monogeneric family Hemistasiidae available in culture. Here we describe two new axenized hemistasiids isolated from Tokyo Bay, Japan. Like in other diplonemids, cellular organization of these heterotrophic protists is characterized by a distinct apical papilla, a tubular cytopharynx contiguous with a deep flagellar pocket, and a highly branched mitochondrion with lamellar cristae. Both hemistasiids also bear a prominent digestive vacuole, peripheral lacunae, and paraflagellar rods, are highly motile and exhibit diverse morphologies in culture. We argue that significant differences in molecular phylogenetics and ultrastructure between these new species and H. phaeocysticola are on the generic level. Therefore, we have established two new genera within Hemistasiidae - Artemidia gen. n. and Namystynia gen. n. to accommodate Artemidia motanka, sp. n. and Namystynia karyoxenos, sp. n., respectively. A. motanka permanently carries tubular extrusomes, while in N. karyoxenos, they are present only in starving cells. An additional remarkable feature of the latter species is the presence, in both the cytoplasm and the nucleus, of the endosymbiotic rickettsiid Candidatus Sneabacter namystus.

Ploeotids Represent Much of the Phylogenetic Diversity of Euglenids

Ploeotids are an assemblage of rigid phagotrophic euglenids that have 10-12 pellicular strips and glide on their posterior flagellum. Molecular phylogenies place them as a poorly resolved, likely paraphyletic assemblage outside the Spirocuta clade of flexible euglenids, which includes the well-known phototrophs and primary osmotrophs. Here, we report SSU rRNA gene sequences from 38 ploeotids, using both single-cell and culture-based methods. Several contain group I or non-canonical introns. Our phylogenetic analyses place ploeotids in 8 distinct clades: Olkasia n. gen., Hemiolia n. gen., Liburna n. gen., Lentomonas, Decastava, Keelungia, Ploeotiidae, and Entosiphon. Ploeotia vitrea, the type of Ploeotia, is closely related to P. oblonga and Serpenomonas costata, but not to Lentomonas. Ploeotia cf. vitrea sensu Lax and Simpson 2013 is not related to P. vitrea and has a different pellicle strip architecture (as imaged by scanning electron microscopy): it instead represents a novel genus and species, Olkasia polycarbonata. We also describe new genera, Hemiolia and Liburna, for the morphospecies Anisonema trepidum and A. glaciale. A recent system proposing 13 suprafamilial taxa that include ploeotids is not supported by our phylogenies. The exact relationships between ploeotid groups remain unresolved and multigene phylogenetics or phylogenomics are needed to address this uncertainty.