More protist genomes needed

Heterotrophic bacteria trigger transcriptome remodelling in the photosynthetic picoeukaryote Micromonas commoda

Marine biogeochemical cycles are built on interactions between surface ocean microbes, particularly those connecting phytoplankton primary producers to heterotrophic bacteria. Details of these associations are not well understood, especially in the case of direct influences of bacteria on phytoplankton physiology. Here we catalogue how the presence of three marine bacteria (Ruegeria pomeroyi DSS-3, Stenotrophomonas sp. SKA14 and Polaribacter dokdonensis MED152) individually and uniquely impact gene expression of the picoeukaryotic alga Micromonas commoda RCC 299. We find a dramatic transcriptomic remodelling by M. commoda after 8 h in co-culture, followed by an increase in cell numbers by 56 h compared with the axenic cultures. Some aspects of the algal transcriptomic response are conserved across all three bacterial co-cultures, including an unexpected reduction in relative expression of photosynthesis and carbon fixation pathways. Expression differences restricted to a single bacterium are also observed, with the Flavobacteriia P. dokdonensis uniquely eliciting changes in relative expression of algal genes involved in biotin biosynthesis and the acquisition and assimilation of nitrogen. This study reveals that M. commoda has rapid and extensive responses to heterotrophic bacteria in ways that are generalizable, as well as in a taxon specific manner, with implications for the diversity of phytoplankton-bacteria interactions ongoing in the surface ocean.

Towards high-throughput parallel imaging and single-cell transcriptomics of microbial eukaryotic plankton

Single-cell transcriptomics has the potential to provide novel insights into poorly studied microbial eukaryotes. Although several such technologies are available and benchmarked on mammalian cells, few have been tested on protists. Here, we applied a microarray single-cell sequencing (MASC-seq) technology, that generates microscope images of cells in parallel with capturing their transcriptomes, on three species representing important plankton groups with different cell structures; the ciliate Tetrahymena thermophila, the diatom Phaeodactylum tricornutum, and the dinoflagellate Heterocapsa sp. Both the cell fixation and permeabilization steps were adjusted. For the ciliate and dinoflagellate, the number of transcripts of microarray spots with single cells were significantly higher than for background spots, and the overall expression patterns were correlated with that of bulk RNA, while for the much smaller diatom cells, it was not possible to separate single-cell transcripts from background. The MASC-seq method holds promise for investigating "microbial dark matter", although further optimizations are necessary to increase the signal-to-noise ratio.

Characterization of Rosculus vilicus sp. nov., a rhizarian amoeba interacting with Mycobacterium avium subsp. paratuberculosis

Free-living amoebae are described as potential reservoirs for pathogenic bacteria in the environment. It has been hypothesized that this might be the case for Mycobacterium avium subsp. paratuberculosis, the bacterium responsible for paratuberculosis. In a previous work, we isolated an amoeba from a water sample in the environment of infected cattle and showed that this amoeba was associated with Mycobacterium avium subsp. paratuberculosis. While a partial 18S rRNA gene has allowed us to suggest that this amoeba was Rosculus-like, at that time we were not able to sub-cultivate it. In the present study, we succeeded in cultivating this strain at 20-25°C. This amoeba is among the smallest (5-7 μm) described. The sequencing of the whole genome allowed us to extract the full 18S rRNA gene and propose this strain as a new species of the Rosculus genus, i.e., R. vilicus. Of note, the mitochondrial genome is particularly large (184,954 bp). Finally, we showed that this amoeba was able to phagocyte Mycobacterium avium subsp. paratuberculosis and that the bacterium was still observed within amoebae after at least 3 days. In conclusion, we characterized a new environmental amoeba species at the cellular and genome level that was able to interact with Mycobacterium avium subsp. paratuberculosis. As a result, R. vilicus is a potential candidate as environmental reservoir for Mycobacterium avium subsp. paratuberculosis but further experiments are needed to test this hypothesis.

Extreme environments offer an unprecedented opportunity to understand microbial eukaryotic ecology, evolution, and genome biology

Research in extreme environments has substantially expanded our understanding of the ecology and evolution of life on Earth, but a major group of organisms has been largely overlooked: microbial eukaryotes (i.e., protists). In this Perspective, we summarize data from over 80 studies of protists in extreme environments and identify focal lineages that are of significant interest for further study, including clades within Echinamoebida, Heterolobosea, Radiolaria, Haptophyta, Oomycota, and Cryptophyta. We argue that extreme environments are prime sampling targets to fill gaps in the eukaryotic tree of life and to increase our understanding of the ecology, metabolism, genome architecture, and evolution of eukaryotic life.

Microeukaryotic predators shape the wastewater microbiome

The physicochemical parameters that shape the prokaryotic community composition in wastewater have been extensively studied. In contrast, it is poorly understood whether and how biotic interactions affect the prokaryotic community composition in wastewater. We used metatranscriptomics data from a bioreactor sampled weekly over 14 months to investigate the wastewater microbiome, including often neglected microeukaryotes. Our analysis revealed that while prokaryotes are unaffected by seasonal changes in water temperature, they are impacted by a seasonal, temperature-induced change in the microeukaryotic community. Our findings suggest that selective predation pressure exerted by microeukaryotes is a significant factor shaping the prokaryotic community in wastewater. This study underscores the importance of investigating the entire wastewater microbiome to develop a comprehensive understanding of wastewater treatment.

Genomic insights into the coupling of a Chlorella-like microeukaryote and sulfur bacteria in the chemocline of permanently stratified Lake Cadagno

Meromictic Lake Cadagno is a permanently stratified system with a persistent microbial bloom within the oxic-anoxic boundary called the chemocline. The association between oxygenic and anoxygenic photosynthesis within the chemocline has been known for at least two decades. Although anoxygenic purple and green sulfur bacteria have been well studied, reports on oxygenic phytoplankton have remained sparse since their discovery in the 1920s. Nearly a century later, this study presents the first near-complete genome of a photosynthetic microbial eukaryote from the chemocline of Lake Cadagno, provisionally named Chlorella-like MAG. The 18.9 Mbp nuclear genome displays a high GC content (71.5%), and the phylogenetic placement suggests that it is a novel species of the genus Chlorella of Chlorophytes. Functional annotation of the Chlorella-like metagenome-assembled genome predicted 10,732 protein-coding genes, with an approximate 0.6% proportion potentially involved in carbon, sulfur, and nitrogen (C, N, and S) metabolism. In addition to C4 photosynthesis, this study detected genes for heat shock proteins (HSPs) in the Chlorella-like algae, consistent with the other Chlorella species. Altogether, the genomic insights in this study suggest the cooperation of photosynthetic algae with phototrophic sulfur bacteria via C, N, and S metabolism, which may aid their collective persistence in the Lake Cadagno chemocline. Furthermore, this work additionally presents the chloroplast genome of Cryptomonas-like species, which was likely to be presumed as cyanobacteria in previous studies because of the presence of phycobilisomes.

Marine ecosystem shifts with deglacial sea-ice loss inferred from ancient DNA shotgun sequencing

Sea ice is a key factor for the functioning and services provided by polar marine ecosystems. However, ecosystem responses to sea-ice loss are largely unknown because time-series data are lacking. Here, we use shotgun metagenomics of marine sedimentary ancient DNA off Kamchatka (Western Bering Sea) covering the last ~20,000 years. We traced shifts from a sea ice-adapted late-glacial ecosystem, characterized by diatoms, copepods, and codfish to an ice-free Holocene characterized by cyanobacteria, salmon, and herring. By providing information about marine ecosystem dynamics across a broad taxonomic spectrum, our data show that ancient DNA will be an important new tool in identifying long-term ecosystem responses to climate transitions for improvements of ocean and cryosphere risk assessments. We conclude that continuing sea-ice decline on the northern Bering Sea shelf might impact on carbon export and disrupt benthic food supply and could allow for a northward expansion of salmon and Pacific herring.

Disentangling the Functional Role of Fungi in Cold Seep Sediment

Cold seeps are biological oases of the deep sea fueled by methane, sulfates, nitrates, and other inorganic sources of energy. Chemolithoautotrophic bacteria and archaea dominate seep sediment, and their diversity and biogeochemical functions are well established. Fungi are likewise diverse, metabolically versatile, and known for their ability to capture and oxidize methane. Still, no study has ever explored the functional role of the mycobiota in the cold seep biome. To assess the complex role of fungi and fill in the gaps, we performed network analysis on 147 samples to disentangle fungal-prokaryotic interactions (fungal 18S and prokaryotic 16S) in the Haima cold seep region. We demonstrated that fungi are central species with high connectivity at the epicenter of prokaryotic networks, reduce their random-attack vulnerability by 60%, and enhance information transfer efficiency by 15%. We then scavenged a global metagenomic and metatranscriptomic data set from 10 cold seep regions for fungal genes of interest (hydrophobins, cytochrome P450s, and ligninolytic family of enzymes); this is the first study to report active transcription of 2,500+ fungal genes in the cold seep sediment. The genera Fusarium and Moniliella were of notable importance and directly correlated with high methane abundance in the sulfate-methane transition zone (SMTZ), likely due to their ability to degrade and solubilize methane and oils. Overall, our results highlight the essential yet overlooked contribution of fungi to cold seep biological networks and the role of fungi in regulating cold seep biogeochemistry. IMPORTANCE The challenges we face when analyzing eukaryotic metagenomic and metatranscriptomic data sets have hindered our understanding of cold seep fungi and microbial eukaryotes. This fact does not make the mycobiota any less critical in mediating cold seep biogeochemistry. On the contrary, many fungal genera can oxidize and solubilize methane, produce methane, and play a unique role in nutrient recycling via saprotrophic enzymatic activity. In this study, we used network analysis to uncover key fungal-prokaryotic interactions that can mediate methane biogeochemistry and metagenomics and metatranscriptomics to report that fungi are transcriptionally active in the cold seep sediment. With concerns over rising methane levels and cold seeps being a pivotal source of global methane input, our holistic understanding of methane biogeochemistry with all domains of life is essential. We ultimately encourage scientists to utilize state-of-the-art tools and multifaceted approaches to uncover the role of microeukaryotic organisms in understudied systems.

Beyond the limits of the unassigned protist microbiome: inferring large-scale spatio-temporal patterns of Syndiniales marine parasites

Marine protists are major components of the oceanic microbiome that remain largely unrepresented in culture collections and genomic reference databases. The exploration of this uncharted protist diversity in oceanic communities relies essentially on studying genetic markers from the environment as taxonomic barcodes. Here we report that across 6 large scale spatio-temporal planktonic surveys, half of the genetic barcodes remain taxonomically unassigned at the genus level, preventing a fine ecological understanding for numerous protist lineages. Among them, parasitic Syndiniales (Dinoflagellata) appear as the least described protist group. We have developed a computational workflow, integrating diverse 18S rDNA gene metabarcoding datasets, in order to infer large-scale ecological patterns at 100% similarity of the genetic marker, overcoming the limitation of taxonomic assignment. From a spatial perspective, we identified 2171 unassigned clusters, i.e., Syndiniales sequences with 100% similarity, exclusively shared between the Tropical/Subtropical Ocean and the Mediterranean Sea among all Syndiniales orders and 25 ubiquitous clusters shared within all the studied marine regions. From a temporal perspective, over 3 time-series, we highlighted 39 unassigned clusters that follow rhythmic patterns of recurrence and are the best indicators of parasite community's variation. These clusters withhold potential as ecosystem change indicators, mirroring their associated host community responses. Our results underline the importance of Syndiniales in structuring planktonic communities through space and time, raising questions regarding host-parasite association specificity and the trophic mode of persistent Syndiniales, while providing an innovative framework for prioritizing unassigned protist taxa for further description.

The clubroot pathogen Plasmodiophora brassicae: A profile update

BACKGROUND

Plasmodiophora brassicae is the causal agent of clubroot disease of cruciferous plants and one of the biggest threats to the rapeseed (Brassica napus) and brassica vegetable industry worldwide.

DISEASE SYMPTOMS

In the advanced stages of clubroot disease wilting, stunting, yellowing, and redness are visible in the shoots. However, the typical symptoms of the disease are the presence of club-shaped galls in the roots of susceptible hosts that block the absorption of water and nutrients.

HOST RANGE

Members of the family Brassicaceae are the primary host of the pathogen, although some members of the family, such as Bunias orientalis, Coronopus squamatus, and Raphanus sativus, have been identified as being consistently resistant to P. brassicae isolates with variable virulence profile.

TAXONOMY

Class: Phytomyxea; Order: Plasmodiophorales; Family: Plasmodiophoraceae; Genus: Plasmodiophora; Species: Plasmodiophora brassicae (Woronin, 1877).

DISTRIBUTION

Clubroot disease is spread worldwide, with reports from all continents except Antarctica. To date, clubroot disease has been reported in more than 80 countries.

PATHOTYPING

Based on its virulence on different hosts, P. brassicae is classified into pathotypes or races. Five main pathotyping systems have been developed to understand the relationship between P. brassicae and its hosts. Nowadays, the Canadian clubroot differential is extensively used in Canada and has so far identified 36 different pathotypes based on the response of a set of 13 hosts.

EFFECTORS AND RESISTANCE

After the identification and characterization of the clubroot pathogen SABATH-type methyltransferase PbBSMT, several other effectors have been characterized. However, no avirulence gene is known, hindering the functional characterization of the five intercellular nucleotide-binding (NB) site leucine-rich-repeat (LRR) receptors (NLRs) clubroot resistance genes validated to date.

IMPORTANT LINK

Canola Council of Canada is constantly updating information about clubroot and P. brassicae as part of their Canola Encyclopedia: https://www.canolacouncil.org/canola-encyclopedia/diseases/clubroot/.

PHYTOSANITARY CATEGORIZATION

PLADBR: EPPO A2 list; Annex designation 9E.

Long-Read Sequencing Improves Recovery of Picoeukaryotic Genomes and Zooplankton Marker Genes from Marine Metagenomes

Long-read sequencing offers the potential to improve metagenome assemblies and provide more robust assessments of microbial community composition and function than short-read sequencing. We applied Pacific Biosciences (PacBio) CCS (circular consensus sequencing) HiFi shotgun sequencing to 14 marine water column samples and compared the results with those for short-read metagenomes from the corresponding environmental DNA samples. We found that long-read metagenomes varied widely in quality and biological information. The community compositions of the corresponding long- and short-read metagenomes were frequently dissimilar, suggesting higher stochasticity and/or bias associated with PacBio sequencing. Long reads provided few improvements to the assembly qualities, gene annotations, and prokaryotic metagenome-assembled genome (MAG) binning results. However, only long reads produced high-quality eukaryotic MAGs and contigs containing complete zooplankton marker gene sequences. These results suggest that high-quality long-read metagenomes can improve marine community composition analyses and provide important insight into eukaryotic phyto- and zooplankton genetics, but the benefits may be outweighed by the inconsistent data quality. IMPORTANCE Ocean microbes provide critical ecosystem services, but most remain uncultivated. Their communities can be studied through shotgun metagenomic sequencing and bioinformatic analyses, including binning draft microbial genomes. However, most sequencing to date has been done using short-read technology, which rarely yields genome sequences of key microbes like SAR11. Long-read sequencing can improve metagenome assemblies but is hampered by technological shortcomings and high costs. In this study, we compared long- and short-read sequencing of marine metagenomes. We found a wide range of long-read metagenome qualities and minimal improvements to microbiome analyses. However, long reads generated draft genomes of eukaryotic algal species and provided full-length marker gene sequences of zooplankton species, including krill and copepods. These results suggest that long-read sequencing can provide greater genetic insight into the wide diversity of eukaryotic phyto- and zooplankton that interact as part of and with the marine microbiome.

Genomic differentiation of three pico-phytoplankton species in the Mediterranean Sea

For more than a decade, high-throughput sequencing has transformed the study of marine planktonic communities and has highlighted the extent of protist diversity in these ecosystems. Nevertheless, little is known relative to their genomic diversity at the species-scale as well as their major speciation mechanisms. An increasing number of data obtained from global scale sampling campaigns is becoming publicly available, and we postulate that metagenomic data could contribute to deciphering the processes shaping protist genomic differentiation in the marine realm. As a proof of concept, we developed a findable, accessible, interoperable and reusable (FAIR) pipeline and focused on the Mediterranean Sea to study three a priori abundant protist species: Bathycoccus prasinos, Pelagomonas calceolata and Phaeocystis cordata. We compared the genomic differentiation of each species in light of geographic, environmental and oceanographic distances. We highlighted that isolation-by-environment shapes the genomic differentiation of B. prasinos, whereas P. cordata is impacted by geographic distance (i.e. isolation-by-distance). At present time, the use of metagenomics to accurately estimate the genomic differentiation of protists remains challenging since coverages are lower compared to traditional population surveys. However, our approach sheds light on ecological and evolutionary processes occurring within natural marine populations and paves the way for future protist population metagenomic studies.

Discrepancies between prokaryotes and eukaryotes need to be considered in soil DNA-based studies

Metabarcoding approaches are exponentially increasing our understanding of soil biodiversity, with a major focus on the bacterial part of the microbiome. Part of the soil diversity are also eukaryotes that include fungi, algae, protists and Metazoa. Nowadays, soil eukaryotes are targeted with the same approaches developed for bacteria and archaea (prokaryotes). However, fundamental differences exist between domains. After providing a short historical overview of the developments of metabarcoding applied to environmental microbiology, we compile the most important differences between domains that prevent direct method transfers between prokaryotic and eukaryotic soil metabarcoding approaches, currently dominated by short-read sequencing. These include the existence of divergent diversity concepts and the variations in eukaryotic morphology that affect sampling and DNA extraction. Furthermore, eukaryotes experienced much more variable evolutionary rates than prokaryotes, which prevent capturing the entire eukaryotic diversity in a soil with a single amplification protocol fit for short-read sequencing. In the final part we focus on future potentials for optimization of eukaryotic metabarcoding that include superior possibility of functionally characterizing eukaryotes and to extend the current information obtained, such as by adding a real quantitative component. This review should optimize future metabarcoding approaches targeting soil eukaryotes and kickstart this promising research direction.

Phylogenomics Supports the Monophyly of Aphelids and Fungi and Identifies New Molecular Synapomorphies

The supergroup Holomycota, composed of Fungi and several related lineages of unicellular organisms (Nucleariida, Rozellida, Microsporidia, and Aphelida), represents one of the major branches in the phylogeny of eukaryotes. Nevertheless, except for the well-established position of Nucleariida as the first holomycotan branch to diverge, the relationships among the other lineages have so far remained unresolved largely owing to the lack of molecular data for some groups. This was notably the case aphelids, a poorly known group of endobiotic phagotrophic protists that feed on algae with cellulose walls. The first molecular phylogenies including aphelids supported their sister relationship with Rozellida and Microsporidia which, collectively, formed a new group called Opisthosporidia (the 'Opisthosporidia hypothesis'). However, recent phylogenomic analyses including massive sequence data from two aphelid genera, Paraphelidium and Amoeboaphelidium, suggested that the aphelids are sister to fungi (the 'Aphelida+Fungi hypothesis'). Should this position be confirmed, aphelids would be key to understanding the early evolution of Holomycota and the origin of Fungi. Here, we carry out phylogenomic analyses with an expanded taxonomic sampling for aphelids after sequencing the transcriptomes of two species of the genus Aphelidium (A. insulamus and A. tribonematis) in order to test these competing hypotheses. Our new phylogenomic analyses including species from the three known aphelid genera strongly rejected the Opisthosporidia hypothesis. Furthermore, comparative genomic analyses further supported the Aphelida+Fungi hypothesis via the identification of 19 orthologous genes exclusively shared by these two lineages. Seven of them originated from ancient horizontal gene transfer events predating the aphelid-fungal split and the remaining 12 likely evolved de novo, constituting additional molecular synapomorphies for this clade. Ancestral trait reconstruction based on our well-resolved phylogeny of Holomycota suggests that the progenitor of both fungi and rozellids, was aphelid-like, having an amoeboflagellate state and likely preying endobiotically on cellulose-containing, cell-walled organisms. Two lineages, which we propose to call Phytophagea and Opisthophagea, evolved from this ancestor. Phytophagea, grouping aphelids and classical fungi, mainly specialized in endobiotic predation of algal cells. Fungi emerged from this lineage after losing phagotrophy in favour of osmotrophy. Opisthophagea, grouping rozellids and Microsporidia, became parasites, mostly of chitin-containing hosts. This lineage entered a progressive reductive process that resulted in a unique lifestyle, especially in the highly derived Microsporidia.

Comparative transcriptomics reveals the molecular toolkit used by an algivorous protist for cell wall perforation

Microbial eukaryotes display a stunning diversity of feeding strategies, ranging from generalist predators to highly specialized parasites. The unicellular "protoplast feeders" represent a fascinating mechanistic intermediate, as they penetrate other eukaryotic cells (algae and fungi) like some parasites but then devour their cell contents by phagocytosis.¹ Besides prey recognition and attachment, this complex behavior involves the local, pre-phagocytotic dissolution of the prey cell wall, which results in well-defined perforations of species-specific size and structure.² Yet the molecular processes that enable protoplast feeders to overcome cell walls of diverse biochemical composition remain unknown. We used the flagellate Orciraptor agilis (Viridiraptoridae, Rhizaria) as a model protoplast feeder and applied differential gene expression analysis to examine its penetration of green algal cell walls. Besides distinct expression changes that reflect major cellular processes (e.g., locomotion and cell division), we found lytic carbohydrate-active enzymes that are highly expressed and upregulated during the attack on the alga. A putative endocellulase (family GH5_5) with a secretion signal is most prominent, and a potential key factor for cell wall dissolution. Other candidate enzymes (e.g., lytic polysaccharide monooxygenases) belong to families that are largely uncharacterized, emphasizing the potential of non-fungal microeukaryotes for enzyme exploration. Unexpectedly, we discovered various chitin-related factors that point to an unknown chitin metabolism in Orciraptor agilis, potentially also involved in the feeding process. Our findings provide first molecular insights into an important microbial feeding behavior and new directions for cell biology research on non-model eukaryotes.

What Can We Learn from -Omics Approaches to Understand Clubroot Disease?

Clubroot is one of the most economically significant diseases worldwide. As a result, many investigations focus on both curing the disease and in-depth molecular studies. Although the first transcriptome dataset for the clubroot disease describing the clubroot disease was published in 2006, many different pathogen-host plant combinations have only recently been investigated and published. Articles presenting -omics data and the clubroot pathogen Plasmodiophora brassicae as well as different host plants were analyzed to summarize the findings in the richness of these datasets. Although genome data for the protist have only recently become available, many effector candidates have been identified, but their functional characterization is incomplete. A better understanding of the life cycle is clearly required to comprehend its function. While only a few proteome studies and metabolome analyses were performed, the majority of studies used microarrays and RNAseq approaches to study transcriptomes. Metabolites, comprising chemical groups like hormones were generally studied in a more targeted manner. Furthermore, functional approaches based on such datasets have been carried out employing mutants, transgenic lines, or ecotypes/cultivars of either Arabidopsis thaliana or other economically important host plants of the Brassica family. This has led to new discoveries of potential genes involved in disease development or in (partial) resistance or tolerance to P. brassicae. The overall contribution of individual experimental setups to a larger picture will be discussed in this review.

Metagenome-assembled genomes of phytoplankton microbiomes from the Arctic and Atlantic Oceans

BACKGROUND

Phytoplankton communities significantly contribute to global biogeochemical cycles of elements and underpin marine food webs. Although their uncultured genomic diversity has been estimated by planetary-scale metagenome sequencing and subsequent reconstruction of metagenome-assembled genomes (MAGs), this approach has yet to be applied for complex phytoplankton microbiomes from polar and non-polar oceans consisting of microbial eukaryotes and their associated prokaryotes.

RESULTS

Here, we have assembled MAGs from chlorophyll a maximum layers in the surface of the Arctic and Atlantic Oceans enriched for species associations (microbiomes) with a focus on pico- and nanophytoplankton and their associated heterotrophic prokaryotes. From 679 Gbp and estimated 50 million genes in total, we recovered 143 MAGs of medium to high quality. Although there was a strict demarcation between Arctic and Atlantic MAGs, adjacent sampling stations in each ocean had 51-88% MAGs in common with most species associations between Prasinophytes and Proteobacteria. Phylogenetic placement revealed eukaryotic MAGs to be more diverse in the Arctic whereas prokaryotic MAGs were more diverse in the Atlantic Ocean. Approximately 70% of protein families were shared between Arctic and Atlantic MAGs for both prokaryotes and eukaryotes. However, eukaryotic MAGs had more protein families unique to the Arctic whereas prokaryotic MAGs had more families unique to the Atlantic.

CONCLUSION

Our study provides a genomic context to complex phytoplankton microbiomes to reveal that their community structure was likely driven by significant differences in environmental conditions between the polar Arctic and warm surface waters of the tropical and subtropical Atlantic Ocean. Video Abstract.

Functional repertoire convergence of distantly related eukaryotic plankton lineages abundant in the sunlit ocean

Marine planktonic eukaryotes play critical roles in global biogeochemical cycles and climate. However, their poor representation in culture collections limits our understanding of the evolutionary history and genomic underpinnings of planktonic ecosystems. Here, we used 280 billion Tara Oceans metagenomic reads from polar, temperate, and tropical sunlit oceans to reconstruct and manually curate more than 700 abundant and widespread eukaryotic environmental genomes ranging from 10 Mbp to 1.3 Gbp. This genomic resource covers a wide range of poorly characterized eukaryotic lineages that complement long-standing contributions from culture collections while better representing plankton in the upper layer of the oceans. We performed the first, to our knowledge, comprehensive genome-wide functional classification of abundant unicellular eukaryotic plankton, revealing four major groups connecting distantly related lineages. Neither trophic modes of plankton nor its vertical evolutionary history could completely explain the functional repertoire convergence of major eukaryotic lineages that coexisted within oceanic currents for millions of years.

The Diverse World of Protists-an Ideal Community with which to Introduce Microscopy in the Microbiology Teaching Laboratory

Protists and other eukaryotes are present in diverse terrestrial and aquatic environments. They can easily be collected from local ponds and streams. Although they are typically larger in size than prokaryotes, making them easy to study with even basic microscopes, microbiology laboratory courses often do not discuss them in detail and may only dedicate a short time to observing preserved samples. This laboratory exercise allows students to develop microscope skills, experience real world application of collecting and processing field samples, and delve deeper into the diverse world of protists and other microbial eukaryotes. Students may also engage by sharing their findings on the website iNaturalist to contribute to the scientific knowledge collected around the world. This laboratory module was initially designed for in-person learning but has been successfully adapted to remote learning and can also be applied to a complete online learning environment.

Evolution of Phytoplankton as Estimated from Genetic Diversity

Phytoplankton are photosynthetic, single-celled organisms producing almost half of all oxygen on Earth and play a central role as prey for higher organisms, making them irreplaceable in the marine food web. As Global Change proceeds, imposing rapidly intensifying selection pressures, phytoplankton are forced to undergo evolution, local extinction, or redistribution, with potentially cascading effects throughout the marine ecosystem. Recent results from the field of population genetics display high levels of standing genetic diversity in natural phytoplankton populations, providing ample ‘evolutionary options’ and implying high adaptive potential to changing conditions. This potential for adaptive evolution is realized in several studies of experimental evolution, even though most of these studies investigate the evolution of only single strains. This, however, shows that phytoplankton not only evolve from standing genetic diversity, but also rely on de novo mutations. Recent global sampling campaigns show that the immense intraspecific diversity of phytoplankton in the marine ecosystem has been significantly underestimated, meaning we are only studying a minor portion of the relevant variability in the context of Global Change and evolution. An increased understanding of genomic diversity is primarily hampered by the low number of ecologically representative reference genomes of eukaryotic phytoplankton and the functional annotation of these. However, emerging technologies relying on metagenome and transcriptome data may offer a more realistic understanding of phytoplankton diversity.

Evolving perspective on the origin and diversification of cellular life and the virosphere

The tree of life (TOL) is a powerful framework to depict the evolutionary history of cellular organisms through time, from our microbial origins to the diversification of multicellular eukaryotes that shape the visible biosphere today. During the past decades, our perception of the TOL has fundamentally changed, in part, due to profound methodological advances, which allowed a more objective approach to studying organismal and viral diversity and led to the discovery of major new branches in the TOL as well as viral lineages. Phylogenetic and comparative genomics analyses of these data have, among others, revolutionized our understanding of the deep roots and diversity of microbial life, the origin of the eukaryotic cell, eukaryotic diversity, as well as the origin, and diversification of viruses. In this review, we provide an overview of some of the recent discoveries on the evolutionary history of cellular organisms and their viruses and discuss a variety of complementary techniques that we consider crucial for making further progress in our understanding of the TOL and its interconnection with the virosphere.

Contrasting protist communities (Cercozoa: Rhizaria) in pristine and earthworm-invaded North American deciduous forests

Earthworms are considered ecosystem engineers due to their fundamental impact on soil structure, soil processes and on other soil biota. An invasion of non-native earthworm species has altered soils of North America since European settlement, a process currently expanding into still earthworm-free forest ecosystems due to continuous spread and increasing soil temperatures owing to climate change. Although earthworms are known to modify soil microbial diversity and activity, it is as yet unclear how eukaryote consumers in soil microbial food webs will be affected. Here, we investigated how earthworm invasion affects the diversity of Cercozoa, one of the most dominant protist taxa in soils. Although the composition of the native cercozoan community clearly shifted in response to earthworm invasion, the communities of the different forests showed distinct responses. We identified 39 operational taxonomic units (OTUs) exclusively indicating earthworm invasion, hinting at an earthworm-associated community of Cercozoa. In particular, Woronina pythii, a hyper-parasite of plant-parasitic Oomycota in American forests, increased strongly in the presence of invasive earthworms, indicating an influence of invasive earthworms on oomycete communities and potentially on forest health, which requires further research.

Metagenomics approach for Polymyxa betae genome assembly enables comparative analysis towards deciphering the intracellular parasitic lifestyle of the plasmodiophorids

Genomic knowledge of the tree of life is biased to specific groups of organisms. For example, only six full genomes are currently available in the rhizaria clade. Here, we have applied metagenomic techniques enabling the assembly of the genome of Polymyxa betae (Rhizaria, Plasmodiophorida) RES F41 isolate from unpurified zoospore holobiont and comparison with the A26-41 isolate. Furthermore, the first P. betae mitochondrial genome was assembled. The two P. betae nuclear genomes were highly similar, each with just ~10.2 k predicted protein coding genes, ~3% of which were unique to each isolate. Extending genomic comparisons revealed a greater overlap with Spongospora subterranea than with Plasmodiophora brassicae, including orthologs of the mammalian cation channel sperm-associated proteins, raising some intriguing questions about zoospore physiology. This work validates our metagenomics pipeline for eukaryote genome assembly from unpurified samples and enriches plasmodiophorid genomics; providing the first full annotation of the P. betae genome.

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.

The Genome of the Haptophyte Diacronema lutheri (Pavlova lutheri, Pavlovales): A Model for Lipid Biosynthesis in Eukaryotic Algae

Haptophytes are biogeochemically and industrially important protists with underexplored genomic diversity. We present a nuclear genome assembly for the class Pavlovales, which was assembled with PacBio long-read data into highly contiguous sequences. We sequenced strain Diacronema lutheri NIVA-4/92, formerly known as Pavlova lutheri, because it has established roles in aquaculture and has been a key organism for studying microalgal lipid biosynthesis. Our data show that D. lutheri has the smallest and most streamlined haptophycean genome assembled to date, with an assembly size of 43.503 Mb and 14,446 protein-coding genes. Together with its high nuclear GC content, Diacronema is an important genus for investigating selective pressures on haptophyte genome evolution, contrasting with the much larger and more repetitive genome of the coccolithophore Emiliania huxleyi. The D. lutheri genome will be a valuable resource for resolving the genetic basis of algal lipid biosynthesis and metabolic remodeling that takes place during adaptation and stress response in natural and engineered environments.

First Genome of Labyrinthula sp., an Opportunistic Seagrass Pathogen, Reveals Novel Insight into Marine Protist Phylogeny, Ecology and CAZyme Cell-Wall Degradation

Labyrinthula spp. are saprobic, marine protists that also act as opportunistic pathogens and are the causative agents of seagrass wasting disease (SWD). Despite the threat of local- and large-scale SWD outbreaks, there are currently gaps in our understanding of the drivers of SWD, particularly surrounding Labyrinthula spp. virulence and ecology. Given these uncertainties, we investigated the Labyrinthula genus from a novel genomic perspective by presenting the first draft genome and predicted proteome of a pathogenic isolate Labyrinthula SR_Ha_C, generated from a hybrid assembly of Nanopore and Illumina sequences. Phylogenetic and cross-phyla comparisons revealed insights into the evolutionary history of Stramenopiles. Genome annotation showed evidence of glideosome-type machinery and an apicoplast protein typically found in protist pathogens and parasites. Proteins involved in Labyrinthula SR_Ha_C's actin-myosin mode of transport, as well as carbohydrate degradation were also prevalent. Further, CAZyme functional predictions revealed a repertoire of enzymes involved in breakdown of cell-wall and carbohydrate storage compounds common to seagrasses. The relatively low number of CAZymes annotated from the genome of Labyrinthula SR_Ha_C compared to other Labyrinthulea species may reflect the conservative annotation parameters, a specialized substrate affinity and the scarcity of characterized protist enzymes. Inherently, there is high probability for finding both unique and novel enzymes from Labyrinthula spp. This study provides resources for further exploration of Labyrinthula spp. ecology and evolution, and will hopefully be the catalyst for new hypothesis-driven SWD research revealing more details of molecular interactions between the Labyrinthula genus and its host substrate.

TRAPID 2.0: a web application for taxonomic and functional analysis of de novo transcriptomes

Advances in high-throughput sequencing have resulted in a massive increase of RNA-Seq transcriptome data. However, the promise of rapid gene expression profiling in a specific tissue, condition, unicellular organism or microbial community comes with new computational challenges. Owing to the limited availability of well-resolved reference genomes, de novo assembled (meta)transcriptomes have emerged as popular tools for investigating the gene repertoire of previously uncharacterized organisms. Yet, despite their potential, these datasets often contain fragmented or contaminant sequences, and their analysis remains difficult. To alleviate some of these challenges, we developed TRAPID 2.0, a web application for the fast and efficient processing of assembled transcriptome data. The initial processing phase performs a global characterization of the input data, providing each transcript with several layers of annotation, comprising structural, functional, and taxonomic information. The exploratory phase enables downstream analyses from the web application. Available analyses include the assessment of gene space completeness, the functional analysis and comparison of transcript subsets, and the study of transcripts in an evolutionary context. A comparison with similar tools highlights TRAPID’s unique features. Finally, analyses performed within TRAPID 2.0 are complemented by interactive data visualizations, facilitating the extraction of new biological insights, as demonstrated with diatom community metatranscriptomes.

Genome sequencing and de novo assembly of the giant unicellular alga Acetabularia acetabulum using droplet MDA

The macroscopic single-celled green alga Acetabularia acetabulum has been a model system in cell biology for more than a century. However, no genomic information is available from this species. Since the alga has a long life cycle, is difficult to grow in dense cultures, and has an estimated diploid genome size of almost 2 Gb, obtaining sufficient genomic material for genome sequencing is challenging. Here, we have attempted to overcome these challenges by amplifying genomic DNA using multiple displacement amplification (MDA) combined with microfluidics technology to distribute the amplification reactions across thousands of microscopic droplets. By amplifying and sequencing DNA from five single cells we were able to recover an estimated ~ 7–11% of the total genome, providing the first draft of the A. acetabulum genome. We highlight challenges associated with genome recovery and assembly of MDA data due to biases arising during genome amplification, and hope that our study can serve as a reference for future attempts on sequencing the genome from non-model eukaryotes.

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.

MOSGA: Modular Open-Source Genome Annotator

MOTIVATION

The generation of high-quality assemblies, even for large eukaryotic genomes, has become a routine task for many biologists thanks to recent advances in sequencing technologies. However, the annotation of these assemblies-a crucial step toward unlocking the biology of the organism of interest-has remained a complex challenge that often requires advanced bioinformatics expertise.

RESULTS

Here, we present MOSGA (Modular Open-Source Genome Annotator), a genome annotation framework for eukaryotic genomes with a user-friendly web-interface that generates and integrates annotations from various tools. The aggregated results can be analyzed with a fully integrated genome browser and are provided in a format ready for submission to NCBI. MOSGA is built on a portable, customizable and easily extendible Snakemake backend, and thus, can be tailored to a wide range of users and projects.

AVAILABILITY AND IMPLEMENTATION

We provide MOSGA as a web service at https://mosga.mathematik.uni-marburg.de and as a docker container at registry.gitlab.com/mosga/mosga: latest. Source code can be found at https://gitlab.com/mosga/mosga.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

•

We optimized single-cell methodology using a broad sample range, for EME

•

We combined bioinformatic and bench protocols into a concise workflow

•

We benchmarked the pipeline and used it on environmental samples

•

We selected a set of QC criteria for best genome quality prediction

Integrating morphology and metagenomics to understand taxonomic variability of Amphisorus (Foraminifera, Miliolida) from Western Australia and Indonesia

Foraminifera are a group of mostly marine protists with high taxonomic diversity. Species identification is often complex, as both morphological and molecular approaches can be challenging due to a lack of unique characters and reference sequences. An integrative approach combining state of the art morphological and molecular tools is therefore promising. In this study, we analysed large benthic Foraminifera of the genus Amphisorus from Western Australia and Indonesia. Based on previous findings on high morphological variability observed in the Soritidae and the discontinuous distribution of Amphisorus along the coast of western Australia, we expected to find multiple morphologically and genetically unique Amphisorus types. In order to gain detailed insights into the diversity of Amphisorus, we applied micro CT scanning and shotgun metagenomic sequencing. We identified four distinct morphotypes of Amphisorus, two each in Australia and Indonesia, and showed that each morphotype is a distinct genotype. Furthermore, metagenomics revealed the presence of three dinoflagellate symbiont clades. The most common symbiont was Fugacium Fr5, and we could show that its genotypes were mostly specific to Amphisorus morphotypes. Finally, we assembled the microbial taxa associated with the two Western Australian morphotypes, and analysed their microbial community composition. Even though each Amphisorus morphotype harboured distinct bacterial communities, sampling location had a stronger influence on bacterial community composition, and we infer that the prokaryotic community is primarily shaped by the microhabitat rather than host identity. The integrated approach combining analyses of host morphology and genetics, dinoflagellate symbionts, and associated microbes leads to the conclusion that we identified distinct, yet undescribed taxa of Amphisorus. We argue that the combination of morphological and molecular methods provides unprecedented insights into the diversity of foraminifera, which paves the way for a deeper understanding of their biodiversity, and facilitates future taxonomic and ecological work.

Prey capture in protists utilizing microtubule filled processes and surface motility

Surface motility, which can be visualized by the movement of live prey organisms, polystyrene microspheres or other inert particles, has been shown to occur in a wide variety of microtubule-filled extensions of the protistan cell surface, although the associated functions remain enigmatic. This article integrates an extensive but poorly known body of literature showing that surface motility, associated with microtubule-filled cell extensions such as flagella, axopodia, actinopodia, reticulopodia, and haptonema, plays a crucial role in protistan prey capture. Surface motility has been most extensively studied in Chlamydomonas where it is responsible for flagella-dependent whole cell gliding motility. The force transduction machinery for gliding motility in Chlamydomonas is intraflagellar transport. Other than in Chlamydomonas, this field has not moved far beyond the descriptive to the mechanistic because of technical challenges associated with many of the protistan organisms that utilize surface motility for prey capture. The purpose of this article is to rekindle interest in the protistan systems that utilize surface motility for prey capture at a time when newly emerging molecular tools for working with protists are poised to reinvigorate a field that has been quiescent too long.

Diversity of electron transport chains in anaerobic protists

Eukaryotic microbes (protists) that occupy low-oxygen environments often have drastically different mitochondrial metabolism compared to their aerobic relatives. A common theme among many anaerobic protists is the serial loss of components of the electron transport chain (ETC). Here, we discuss the diversity of the ETC across the tree of eukaryotes and review hypotheses for how ETCs are modified, and ultimately lost, in protists. We find that while protists have converged to some of the same metabolism as anaerobic animals, there are clear protist-specific strategies to thrive without oxygen.

Perspectives from Ten Years of Protist Studies by High-Throughput Metabarcoding

During the last decade, high-throughput metabarcoding became routine for analyzing protistan diversity and distributions in nature. Amid a multitude of exciting findings, scientists have also identified and addressed technical and biological limitations, although problems still exist for inference of meaningful taxonomic and ecological knowledge based on short DNA sequences. Given the extensive use of this approach, it is critical to settle our understanding on its strengths and weaknesses and to synthesize up-to-date methodological and conceptual trends. This article summarizes key scientific and technical findings, and identifies current and future directions in protist research that uses metabarcoding.

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.

Novel phylogenetic methods are needed for understanding gene function in the era of mega-scale genome sequencing

Ongoing large-scale genome sequencing projects are forecasting a data deluge that will almost certainly overwhelm current analytical capabilities of evolutionary genomics. In contrast to population genomics, there are no standardized methods in evolutionary genomics for extracting evolutionary and functional (e.g. gene-trait association) signal from genomic data. Here, we examine how current practices of multi-species comparative genomics perform in this aspect and point out that many genomic datasets are under-utilized due to the lack of powerful methodologies. As a result, many current analyses emphasize gene families for which some functional data is already available, resulting in a growing gap between functionally well-characterized genes/organisms and the universe of unknowns. This leaves unknown genes on the 'dark side' of genomes, a problem that will not be mitigated by sequencing more and more genomes, unless we develop tools to infer functional hypotheses for unknown genes in a systematic manner. We provide an inventory of recently developed methods capable of predicting gene-gene and gene-trait associations based on comparative data, then argue that realizing the full potential of whole genome datasets requires the integration of phylogenetic comparative methods into genomics, a rich but underutilized toolbox for looking into the past.

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.

Why the -omic future of Apicomplexa should include gregarines

Gregarines, a polyphyletic group of apicomplexan parasites infecting mostly non-vertebrates hosts, remains poorly known at taxonomic, phylogenetic and genomic levels. However, it represents an essential group for understanding evolutionary history and adaptive capacities of apicomplexan parasites to the remarkable diversity of their hosts. Because they have a mostly extracellular lifestyle, gregarines have developed other cellular developmental forms and host-parasite interactions, compared with their much better studied apicomplexan cousins, intracellular parasites of vertebrates (Hemosporidia, Coccidia, Cryptosporidia). This review highlights the promises offered by the molecular exploration of gregarines, that have been until now left on the side of the road of the comparative -omic exploration of apicomplexan parasites. Elucidating molecular bases for both their ultrastructural, functional and behavioural similarities and differences, compared with those of the typical apicomplexan models, is expected to provide entirely novel clues on the adaptive capacities developed by Apicomplexa over evolution. A challenge remains to identify which gregarines should be explored in priority, as recent metadata from open and host-associated environments have confirmed how underestimated is our current view on true gregarine biodiversity. It is now time to turn to gregarines to widen the currently highly skewed view we have of adaptive mechanisms developed by Apicomplexa.

Transcriptome reconstruction and functional analysis of eukaryotic marine plankton communities via high-throughput metagenomics and metatranscriptomics

Large-scale metagenomic and metatranscriptomic data analyses are often restricted by their gene-centric approach, limiting the ability to understand organismal and community biology. De novo assembly of large and mosaic eukaryotic genomes from complex meta-omics data remains a challenging task, especially in comparison with more straightforward bacterial and archaeal systems. Here, we use a transcriptome reconstruction method based on clustering co-abundant genes across a series of metagenomic samples. We investigated the co-abundance patterns of ∼37 million eukaryotic unigenes across 365 metagenomic samples collected during the Tara Oceans expeditions to assess the diversity and functional profiles of marine plankton. We identified ∼12,000 co-abundant gene groups (CAGs), encompassing ∼7 million unigenes, including 924 metagenomics-based transcriptomes (MGTs, CAGs larger than 500 unigenes). We demonstrated the biological validity of the MGT collection by comparing individual MGTs with available references. We identified several key eukaryotic organisms involved in dimethylsulfoniopropionate (DMSP) biosynthesis and catabolism in different oceanic provinces, thus demonstrating the potential of the MGT collection to provide functional insights on eukaryotic plankton. We established the ability of the MGT approach to capture interspecies associations through the analysis of a nitrogen-fixing haptophyte-cyanobacterial symbiotic association. This MGT collection provides a valuable resource for analyses of eukaryotic plankton in the open ocean by giving access to the genomic content and functional potential of many ecologically relevant eukaryotic species.

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.

All Eukaryotes Are Sexual, unless Proven Otherwise: Many So-Called Asexuals Present Meiotic Machinery and Might Be Able to Have Sex

Here a wide distribution of meiotic machinery is shown, indicating the occurrence of sexual processes in all major eukaryotic groups, without exceptions, including the putative "asexuals." Meiotic machinery has evolved from archaeal DNA repair machinery by means of ancestral gene duplications. Sex is very conserved and widespread in eukaryotes, even though its evolutionary importance is still a matter of debate. The main processes in sex are plasmogamy, followed by karyogamy and meiosis. Meiosis is fundamentally a chromosomal process, which implies recombination and ploidy reduction. Several eukaryotic lineages are proposed to be asexual because their sexual processes are never observed, but presumed asexuality correlates with lack of study. The authors stress the complete lack of meiotic proteins in nucleomorphs and their almost complete loss in the fungus Malassezia. Inversely, complete sets of meiotic proteins are present in fungal groups Glomeromycotina, Trichophyton, and Cryptococcus. Endosymbiont Perkinsela and endoparasitic Microsporidia also present meiotic proteins.

Genome-wide identification of genes encoding putative secreted E3 ubiquitin ligases and functional characterization of PbRING1 in the biotrophic protist Plasmodiophora brassicae

The E3 ubiquitin ligases are key regulators of protein ubiquitination, which have been shown to be involved in a variety of cellular responses to both biotic and abiotic stresses in eukaryotes. However, the E3 ubiquitin ligase homologues in the soil-borne plant pathogen Plasmodiophora brassicae, the causal agent of clubroot disease of crucifer crops worldwide, remain largely unknown. In this study, we characterized secreted E3 ubiquitin ligases, a group of proteins known to be involved in virulence in many pathogens, in a plasmodiophorid P. brassicae. Genome-wide search in the P. brassicae genome retrieved 139 putative E3 ubiquitin ligases, comprising of 115 RING, 15 HECT, 1 HECT-like, and 8 U-box E3 ubiquitin ligases. Among these E3 ubiquitin ligases, 11 RING, 1 U-box, and 3 HECT were found to harbor signal peptide. Based on published RNA-seq data (Schwelm et al. in Sci Rep 5:11153, 2015), we found that these genes were differentially expressed in distinct life stages including germinating spores, maturing spores, and plasmodia. We characterized one potential secreted E3 ubiquitin ligase, PbRING1 (PBRA_000499). Yeast invertase assay showed that PbRING1 harbors a functional N-terminal signal peptide. PbRING1 also harbors a really interested new gene (RING) domain at its C terminus, which was found to display the E3 ligase activity in vitro. Collectively, this study provides a comprehensive insight into the reservoir of putative secreted E3 ligases in P. brassicae.

Ubiquitin fusion proteins in algae: implications for cell biology and the spread of photosynthesis

Background

The process of gene fusion involves the formation of a single chimeric gene from multiple complete or partial gene sequences. Gene fusion is recognized as an important mechanism by which genes and their protein products can evolve new functions. The presence-absence of gene fusions can also be useful characters for inferring evolutionary relationships between organisms.

Results

Here we show that the nuclear genomes of two unrelated single-celled algae, the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans, possess an unexpected diversity of genes for ubiquitin fusion proteins, including novel arrangements in which ubiquitin occupies amino-terminal, carboxyl-terminal, and internal positions relative to its fusion partners. We explore the evolution of the ubiquitin multigene family in both genomes, and show that both algae possess a gene encoding an ubiquitin-nickel superoxide dismutase fusion protein (Ubiq-NiSOD) that is widely but patchily distributed across the eukaryotic tree of life – almost exclusively in phototrophs.

Conclusion

Our results suggest that ubiquitin fusion proteins are more common than currently appreciated; because of its small size, the ubiquitin coding region can go undetected when gene predictions are carried out in an automated fashion. The punctate distribution of the Ubiq-NiSOD fusion across the eukaryotic tree could serve as a beacon for the spread of plastids from eukaryote to eukaryote by secondary and/or tertiary endosymbiosis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5412-4) contains supplementary material, which is available to authorized users.

Genomes OnLine database (GOLD) v.7: updates and new features

The Genomes Online Database (GOLD) (https://gold.jgi.doe.gov) is an open online resource, which maintains an up-to-date catalog of genome and metagenome projects in the context of a comprehensive list of associated metadata. Information in GOLD is organized into four levels: Study, Biosample/Organism, Sequencing Project and Analysis Project. Currently GOLD hosts information on 33 415 Studies, 49 826 Biosamples, 313 324 Organisms, 215 881 Sequencing Projects and 174 454 Analysis Projects with a total of 541 metadata fields, of which 80 are based on controlled vocabulary (CV) terms. GOLD provides a user-friendly web interface to browse sequencing projects and launch advanced search tools across four classification levels. Users submit metadata on a wide range of Sequencing and Analysis Projects in GOLD before depositing sequence data to the Integrated Microbial Genomes (IMG) system for analysis. GOLD conforms with and supports the rules set by the Genomic Standards Consortium (GSC) Minimum Information standards. The current version of GOLD (v.7) has seen the number of projects and associated metadata increase exponentially over the years. This paper provides an update on the current status of GOLD and highlights the new features added over the last two years.

The evolution of ependymin-related proteins

Background

Ependymins were originally defined as fish-specific secreted glycoproteins involved in central nervous system plasticity and memory formation. Subsequent research revealed that these proteins represent a fish-specific lineage of a larger ependymin-related protein family (EPDRs). EPDRs have now been identified in a number of bilaterian animals and have been implicated in diverse non-neural functions. The recent discoveries of putative EPDRs in unicellular holozoans and an expanded EPDR family with potential roles in conspecific communication in crown-of-thorns starfish suggest that the distribution and diversity of EPDRs is significantly broader than currently understood.

Results

We undertook a systematic survey to determine the distribution and evolution of EPDRs in eukaryotes. In addition to Bilateria, EPDR genes were identified in Cnidaria, Placozoa, Porifera, Choanoflagellatea, Filasterea, Apusozoa, Amoebozoa, Charophyta and Percolozoa, and tentatively in Cercozoa and the orphan group Malawimonadidae. EPDRs appear to be absent from prokaryotes and many eukaryote groups including ecdysozoans, fungi, stramenopiles, alveolates, haptistans and cryptistans. The EPDR family can be divided into two major clades and has undergone lineage-specific expansions in a number of metazoan lineages, including in poriferans, molluscs and cephalochordates. Variation in a core set of conserved residues in EPDRs reveals the presence of three distinct protein types; however, 3D modelling predicts overall protein structures to be similar.

Conclusions

Our results reveal an early eukaryotic origin of the EPDR gene family and a dynamic pattern of gene duplication and gene loss in animals. This research provides a phylogenetic framework for the analysis of the functional evolution of this gene family.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1306-y) contains supplementary material, which is available to authorized users.

Nuclear genome sequence of the plastid-lacking cryptomonad Goniomonas avonlea provides insights into the evolution of secondary plastids

Background

The evolution of photosynthesis has been a major driver in eukaryotic diversification. Eukaryotes have acquired plastids (chloroplasts) either directly via the engulfment and integration of a photosynthetic cyanobacterium (primary endosymbiosis) or indirectly by engulfing a photosynthetic eukaryote (secondary or tertiary endosymbiosis). The timing and frequency of secondary endosymbiosis during eukaryotic evolution is currently unclear but may be resolved in part by studying cryptomonads, a group of single-celled eukaryotes comprised of both photosynthetic and non-photosynthetic species. While cryptomonads such as Guillardia theta harbor a red algal-derived plastid of secondary endosymbiotic origin, members of the sister group Goniomonadea lack plastids. Here, we present the genome of Goniomonas avonlea—the first for any goniomonad—to address whether Goniomonadea are ancestrally non-photosynthetic or whether they lost a plastid secondarily.

Results

We sequenced the nuclear and mitochondrial genomes of Goniomonas avonlea and carried out a comparative analysis of Go. avonlea, Gu. theta, and other cryptomonads. The Go. avonlea genome assembly is ~ 92 Mbp in size, with 33,470 predicted protein-coding genes. Interestingly, some metabolic pathways (e.g., fatty acid biosynthesis) predicted to occur in the plastid and periplastidal compartment of Gu. theta appear to operate in the cytoplasm of Go. avonlea, suggesting that metabolic redundancies were generated during the course of secondary plastid integration. Other cytosolic pathways found in Go. avonlea are not found in Gu. theta, suggesting secondary loss in Gu. theta and other plastid-bearing cryptomonads. Phylogenetic analyses revealed no evidence for algal endosymbiont-derived genes in the Go. avonlea genome. Phylogenomic analyses point to a specific relationship between Cryptista (to which cryptomonads belong) and Archaeplastida.

Conclusion

We found no convincing genomic or phylogenomic evidence that Go. avonlea evolved from a secondary red algal plastid-bearing ancestor, consistent with goniomonads being ancestrally non-photosynthetic eukaryotes. The Go. avonlea genome sheds light on the physiology of heterotrophic cryptomonads and serves as an important reference point for studying the metabolic “rewiring” that took place during secondary plastid integration in the ancestor of modern-day Cryptophyceae.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0593-5) contains supplementary material, which is available to authorized users.