Re-assembly, quality evaluation, and annotation of 678 microbial eukaryotic reference transcriptomes

Long-read DNA sequencing reveals the organization of the mitochondrial genome in the early-branching dinoflagellate Oxyrrhis marina

The mitochondrial genomes of dinoflagellate protists are remarkable for their highly fragmented and heterogeneous organization. Early attempts to determine their structure without 'next-generation' DNA sequencing failed to recover a defined genome. Still, it coincided in showing that the proteins coding genes, three in total, and parts of the ribosomal RNA genes were spread across a diffuse assortment of small linear fragments. In contrast, a recent study employed Illumina sequencing to assemble a 326 kbp long single-molecule, circular mitochondrial genome in the symbiotic dinoflagellate Breviolum minutum. Here, we used a combination of short- and long-read massively-parallel DNA sequencing to analyze further the mitochondrial DNA of the early-branching dinoflagellate Oxyrrhis marina. We found that the mitochondrial genome of O. marina consists of 3 linear chromosomes sized 15.9, 33.8 and 40.6 kbp for a total of 90.3 kbp. It contains the cox1, cox3 and cob genes, the same three proteins encoded in the mitochondrion of all myzozoans (Apicomplexa and Dinophyceae), some fragments of ribosomal RNA genes as well as many non-functional gene fragments and extensive noncoding DNA. Our analysis unveiled segments syntenic patterns and rearrangements encompassing coding and non-coding regions, suggesting that recombination is a pervasive process driving the evolution of these genomes.

Missing microbial eukaryotes and misleading meta-omic conclusions

Meta-omics is commonly used for large-scale analyses of microbial eukaryotes, including species or taxonomic group distribution mapping, gene catalog construction, and inference on the functional roles and activities of microbial eukaryotes in situ. Here, we explore the potential pitfalls of common approaches to taxonomic annotation of protistan meta-omic datasets. We re-analyze three environmental datasets at three levels of taxonomic hierarchy in order to illustrate the crucial importance of database completeness and curation in enabling accurate environmental interpretation. We show that taxonomic membership of sequence clusters estimates community composition more accurately than returning exact sequence labels, and overlap between clusters can address database shortcomings. Clustering approaches can be applied to diverse environments while continuing to exploit the wealth of annotation data collated in databases, and selecting and evaluating these databases is a critical part of correctly annotating protistan taxonomy in environmental datasets. We argue that ongoing curation of genetic resources is crucial in accurately annotating protists in in situ meta-omic datasets. Moreover, we propose that precise taxonomic annotation of meta-omic data is a clustering problem rather than a feasible alignment problem.

Microeukaryote metabolism across the western North Atlantic Ocean revealed through autonomous underwater profiling

Microeukaryotes are key contributors to marine carbon cycling. Their physiology, ecology, and interactions with the chemical environment are poorly understood in offshore ecosystems, and especially in the deep ocean. Using the Autonomous Underwater Vehicle Clio, microbial communities along a 1050 km transect in the western North Atlantic Ocean were surveyed at 10-200 m vertical depth increments to capture metabolic signatures spanning oligotrophic, continental margin, and productive coastal ecosystems. Microeukaryotes were examined using a paired metatranscriptomic and metaproteomic approach. Here we show a diverse surface assemblage consisting of stramenopiles, dinoflagellates and ciliates represented in both the transcript and protein fractions, with foraminifera, radiolaria, picozoa, and discoba proteins enriched at >200 m, and fungal proteins emerging in waters >3000 m. In the broad microeukaryote community, nitrogen stress biomarkers were found at coastal sites, with phosphorus stress biomarkers offshore. This multi-omics dataset broadens our understanding of how microeukaryotic taxa and their functional processes are structured along environmental gradients of temperature, light, and nutrients.

Seasonal metabolic dynamics of microeukaryotic plankton: a year-long metatranscriptomic study in a temperate sea

Seasonal fluctuations profoundly affect marine microeukaryotic plankton composition and metabolism, but accurately tracking these changes has been a long-standing challenge. In this study, we present a year-long metatranscriptomic data set from the Southern Bight of the North Sea, shedding light on the seasonal dynamics in temperate plankton ecosystems. We observe distinct shifts in active plankton species and their metabolic processes in response to seasonal changes. We characterized the metabolic signatures of different seasonal phases in detail, thereby revealing the metabolic versatility of dinoflagellates, the heterotrophic dietary strategy of Phaeocystis during its late-stage blooms, and stark variations in summer and fall diatom abundance and metabolic activity across nearby sampling stations. Our data illuminate the varied contributions of microeukaryotic taxa to biomass production and nutrient cycling at different times of the year and allow delineation of their ecological niches.

IMPORTANCE

Ecosystem composition and metabolic functions of temperate marine microeukaryote plankton are strongly influenced by seasonal dynamics. Although monitoring of species composition of microeukaryotes has expanded recently, few methods also contain seasonally resolved information on ecosystem functioning. We generated a year-long spatially resolved metatranscriptomic data set to assess seasonal dynamics of microeukaryote species and their associated metabolic functions in the Southern Bight of the North Sea. Our study underscores the potential of metatranscriptomics as a powerful tool for advancing our understanding of marine ecosystem functionality and resilience in response to environmental changes, emphasizing its potential in continuous marine ecosystem monitoring to enhance our ecological understanding of the ocean's eukaryotic microbiome.

First regional reference database of northern Adriatic diatom transcriptomes

Marine microbial communities form the basis for the functioning of marine ecosystems and the conservation of biodiversity. With the application of metagenomics and metatranscriptomics in marine environmental studies, significant progress has been made in analysing the functioning of microbial communities as a whole. These molecular techniques are highly dependent on reliable, well-characterised, comprehensive and taxonomically diverse sequenced reference transcriptomes of microbial organisms. Here we present a set of 12 individual transcriptome assemblies derived from 6 representative diatom species from the northern Adriatic Sea grown under 2 environmentally relevant growth conditions (phosphate replete vs. phosphate deprived). After filtering the reads and assembly, an average number of 64,932 transcripts per assembly was obtained, of which an average of 8856 were assigned to functionally known proteins. Of all assigned transcripts, an average of 6483 proteins were taxonomically assigned to diatoms (Bacillariophyta). On average, a higher number of assigned proteins was detected in the transcriptome assemblies of diatoms grown under replete media condition. On average, 50% of the mapped proteins were shared between the two growth conditions. All recorded proteins in the dataset were classified into 24 COG categories, with approximately 25% belonging to the unknown function and the remaining 75% belonging to all other categories. The resulting diatom reference database for the northern Adriatic, focussing on the response to nutrient limitation as characteristic for the region and predicted for the future world oceans, provides a valuable resource for analysing environmental metatranscriptome and metagenome data. Each northern Adriatic transcriptome can also be used by itself as a reference database for the (meta)transcriptomes and gene expression studies of the associated species that will be generated in the future.

DeepMicroClass sorts metagenomic contigs into prokaryotes, eukaryotes and viruses

Sequence classification facilitates a fundamental understanding of the structure of microbial communities. Binary metagenomic sequence classifiers are insufficient because environmental metagenomes are typically derived from multiple sequence sources. Here we introduce a deep-learning based sequence classifier, DeepMicroClass, that classifies metagenomic contigs into five sequence classes, i.e. viruses infecting prokaryotic or eukaryotic hosts, eukaryotic or prokaryotic chromosomes, and prokaryotic plasmids. DeepMicroClass achieved high performance for all sequence classes at various tested sequence lengths ranging from 500 bp to 100 kbps. By benchmarking on a synthetic dataset with variable sequence class composition, we showed that DeepMicroClass obtained better performance for eukaryotic, plasmid and viral contig classification than other state-of-the-art predictors. DeepMicroClass achieved comparable performance on viral sequence classification with geNomad and VirSorter2 when benchmarked on the CAMI II marine dataset. Using a coastal daily time-series metagenomic dataset as a case study, we showed that microbial eukaryotes and prokaryotic viruses are integral to microbial communities. By analyzing monthly metagenomes collected at HOT and BATS, we found relatively higher viral read proportions in the subsurface layer in late summer, consistent with the seasonal viral infection patterns prevalent in these areas. We expect DeepMicroClass will promote metagenomic studies of under-appreciated sequence types.

Noncanonical usage of stop codons in ciliates expands proteins with structurally flexible Q-rich motifs

Serine(S)/threonine(T)-glutamine(Q) cluster domains (SCDs), polyglutamine (polyQ) tracts and polyglutamine/asparagine (polyQ/N) tracts are Q-rich motifs found in many proteins. SCDs often are intrinsically disordered regions that mediate protein phosphorylation and protein-protein interactions. PolyQ and polyQ/N tracts are structurally flexible sequences that trigger protein aggregation. We report that due to their high percentages of STQ or STQN amino acid content, four SCDs and three prion-causing Q/N-rich motifs of yeast proteins possess autonomous protein expression-enhancing activities. Since these Q-rich motifs can endow proteins with structural and functional plasticity, we suggest that they represent useful toolkits for evolutionary novelty. Comparative Gene Ontology (GO) analyses of the near-complete proteomes of 26 representative model eukaryotes reveal that Q-rich motifs prevail in proteins involved in specialized biological processes, including Saccharomyces cerevisiae RNA-mediated transposition and pseudohyphal growth, Candida albicans filamentous growth, ciliate peptidyl-glutamic acid modification and microtubule-based movement, Tetrahymena thermophila xylan catabolism and meiosis, Dictyostelium discoideum development and sexual cycles, Plasmodium falciparum infection, and the nervous systems of Drosophila melanogaster, Mus musculus and Homo sapiens. We also show that Q-rich-motif proteins are expanded massively in 10 ciliates with reassigned TAAQ and TAGQ codons. Notably, the usage frequency of CAGQ is much lower in ciliates with reassigned TAAQ and TAGQ codons than in organisms with expanded and unstable Q runs (e.g. D. melanogaster and H. sapiens), indicating that the use of noncanonical stop codons in ciliates may have coevolved with codon usage biases to avoid triplet repeat disorders mediated by CAG/GTC replication slippage.

Diurnal Rhythms in the Red Seaweed Gracilariopsis chorda are Characterized by Unique Regulatory Networks of Carbon Metabolism

Cellular and physiological cycles are driven by endogenous pacemakers, the diurnal and circadian rhythms. Key functions such as cell cycle progression and cellular metabolism are under rhythmic regulation, thereby maintaining physiological homeostasis. The photoreceptors phytochrome and cryptochrome, in response to light cues, are central input pathways for physiological cycles in most photosynthetic organisms. However, among Archaeplastida, red algae are the only taxa that lack phytochromes. Current knowledge about oscillatory rhythms is primarily derived from model species such as Arabidopsis thaliana and Chlamydomonas reinhardtii in the Viridiplantae, whereas little is known about these processes in other clades of the Archaeplastida, such as the red algae (Rhodophyta). We used genome-wide expression profiling of the red seaweed Gracilariopsis chorda and identified 3,098 rhythmic genes. Here, we characterized possible cryptochrome-based regulation and photosynthetic/cytosolic carbon metabolism in this species. We found a large family of cryptochrome genes in G. chorda that display rhythmic expression over the diurnal cycle and may compensate for the lack of phytochromes in this species. The input pathway gates regulatory networks of carbon metabolism which results in a compact and efficient energy metabolism during daylight hours. The system in G. chorda is distinct from energy metabolism in most plants, which activates in the dark. The green lineage, in particular, land plants, balance water loss and CO2 capture in terrestrial environments. In contrast, red seaweeds maintain a reduced set of photoreceptors and a compact cytosolic carbon metabolism to thrive in the harsh abiotic conditions typical of intertidal zones.

Massive genome reduction predates the divergence of Symbiodiniaceae dinoflagellates

Dinoflagellates in the family Symbiodiniaceae are taxonomically diverse, predominantly symbiotic lineages that are well-known for their association with corals. The ancestor of these taxa is believed to have been free-living. The establishment of symbiosis (i.e. symbiogenesis) is hypothesized to have occurred multiple times during Symbiodiniaceae evolution, but its impact on genome evolution of these taxa is largely unknown. Among Symbiodiniaceae, the genus Effrenium is a free-living lineage that is phylogenetically positioned between two robustly supported groups of genera within which symbiotic taxa have emerged. The apparent lack of symbiogenesis in Effrenium suggests that the ancestral features of Symbiodiniaceae may have been retained in this lineage. Here, we present de novo assembled genomes (1.2-1.9 Gbp in size) and transcriptome data from three isolates of Effrenium voratum and conduct a comparative analysis that includes 16 Symbiodiniaceae taxa and the other dinoflagellates. Surprisingly, we find that genome reduction, which is often associated with a symbiotic lifestyle, predates the origin of Symbiodiniaceae. The free-living lifestyle distinguishes Effrenium from symbiotic Symbiodiniaceae vis-à-vis their longer introns, more-extensive mRNA editing, fewer (~30%) lineage-specific gene sets, and lower (~10%) level of pseudogenization. These results demonstrate how genome reduction and the adaptation to distinct lifestyles intersect to drive diversification and genome evolution of Symbiodiniaceae.

Functional diversity of microbial eukaryotes in a meromictic lake: Coupling between metatranscriptomic and a trait-based approach

The advent of high-throughput sequencing has led to the discovery of a considerable diversity of microbial eukaryotes in aquatic ecosystems, nevertheless, their function and contribution to the trophic food web functioning remain poorly characterized especially in freshwater ecosystems. Based on metabarcoding data obtained from a meromictic lake ecosystem (Pavin, France), we performed a morpho-physio-phenological traits-based approach to infer functional groups of microbial eukaryotes. Metatranscriptomic data were also analysed to assess the metabolic potential of these groups across the diel cycle, size fraction, sampling depth, and periods. Our analysis highlights a huge microbial eukaryotic diversity in the monimolimnion characterized by numerous saprotrophs expressing transcripts related to sulfur and nitrate metabolism as well as dissolved and particulate organic matter degradation. We also describe strong seasonal variations of microbial eukaryotes in the mixolimnion, especially for parasites and mixoplankton. It appears that the water mixing (occurring during spring and autumn) which benefits photosynthetic host communities also promotes parasitic fungi dissemination and over-expression of genes involved in the zoospore phototaxis and stage transition in the parasitic cycle. Mixoplanktonic haptophytes over-expressing photosynthesis-, endocytosis- and phagosome-linked genes under nutrient limitation also suggest that phagotrophy may provide them an advantage over non-phagotrophic phytoplankton.

Multi-omics analysis reveals the molecular response to heat stress in a "red tide" dinoflagellate

BACKGROUND

"Red tides" are harmful algal blooms caused by dinoflagellate microalgae that accumulate toxins lethal to other organisms, including humans via consumption of contaminated seafood. These algal blooms are driven by a combination of environmental factors including nutrient enrichment, particularly in warm waters, and are increasingly frequent. The molecular, regulatory, and evolutionary mechanisms that underlie the heat stress response in these harmful bloom-forming algal species remain little understood, due in part to the limited genomic resources from dinoflagellates, complicated by the large sizes of genomes, exhibiting features atypical of eukaryotes.

RESULTS

We present the de novo assembled genome (~ 4.75 Gbp with 85,849 protein-coding genes), transcriptome, proteome, and metabolome from Prorocentrum cordatum, a globally abundant, bloom-forming dinoflagellate. Using axenic algal cultures, we study the molecular mechanisms that underpin the algal response to heat stress, which is relevant to current ocean warming trends. We present the first evidence of a complementary interplay between RNA editing and exon usage that regulates the expression and functional diversity of biomolecules, reflected by reduction in photosynthesis, central metabolism, and protein synthesis. These results reveal genomic signatures and post-transcriptional regulation for the first time in a pelagic dinoflagellate.

CONCLUSIONS

Our multi-omics analyses uncover the molecular response to heat stress in an important bloom-forming algal species, which is driven by complex gene structures in a large, high-G+C genome, combined with multi-level transcriptional regulation. The dynamics and interplay of molecular regulatory mechanisms may explain in part how dinoflagellates diversified to become some of the most ecologically successful organisms on Earth.

Multiple parallel origins of parasitic Marine Alveolates

Microbial eukaryotes are important components of marine ecosystems, and the Marine Alveolates (MALVs) are consistently both abundant and diverse in global environmental sequencing surveys. MALVs are dinoflagellates that are thought to be parasites of other protists and animals, but the lack of data beyond ribosomal RNA gene sequences from all but a few described species means much of their biology and evolution remain unknown. Using single-cell transcriptomes from several MALVs and their free-living relatives, we show that MALVs evolved independently from two distinct, free-living ancestors and that their parasitism evolved in parallel. Phylogenomics shows one subgroup (MALV-II and -IV, or Syndiniales) is related to a novel lineage of free-living, eukaryovorous predators, the eleftherids, while the other (MALV-I, or Ichthyodinida) is related to the free-living predator Oxyrrhis and retains proteins targeted to a non-photosynthetic plastid. Reconstructing the evolution of photosynthesis, plastids, and parasitism in early-diverging dinoflagellates shows a number of parallels with the evolution of their apicomplexan sisters. In both groups, similar forms of parasitism evolved multiple times and photosynthesis was lost many times. By contrast, complete loss of the plastid organelle is infrequent and, when this does happen, leaves no residual genes.

Increased Sampling and Intracomplex Homologies Favor Vertical Over Horizontal Inheritance of the Dam1 Complex

Kinetochores connect chromosomes to spindle microtubules to ensure their correct segregation during cell division. Kinetochores of human and yeasts are largely homologous, their ability to track depolymerizing microtubules, however, is carried out by the nonhomologous complexes Ska1-C and Dam1-C, respectively. We previously reported the unique anti-correlating phylogenetic profiles of Dam1-C and Ska-C found among a wide variety of eukaryotes. Based on these profiles and the limited presence of Dam1-C, we speculated that horizontal gene transfer could have played a role in the evolutionary history of Dam1-C. Here, we present an expanded analysis of Dam1-C evolution, using additional genome as well as transcriptome sequences and recently published 3D structures. This analysis revealed a wider and more complete presence of Dam1-C in Cryptista, Rhizaria, Ichthyosporea, CRuMs, and Colponemidia. The fungal Dam1-C cryo-EM structure supports earlier hypothesized intracomplex homologies, which enables the reconstruction of rooted and unrooted phylogenies. The rooted tree of concatenated Dam1-C subunits is statistically consistent with the species tree of eukaryotes, suggesting that Dam1-C is ancient, and that the present-day phylogenetic distribution is best explained by multiple, independent losses and no horizontal gene transfer was involved. Furthermore, we investigated the ancient origin of Dam1-C via profile-versus-profile searches. Homology among 8 out of the 10 Dam1-C subunits suggests that the complex largely evolved from a single multimerizing subunit that diversified into a hetero-octameric core via stepwise subunit duplication and subfunctionalization of the subunits before the origin of the last eukaryotic common ancestor.

Reverse engineering environmental metatranscriptomes clarifies best practices for eukaryotic assembly

BACKGROUND

Diverse communities of microbial eukaryotes in the global ocean provide a variety of essential ecosystem services, from primary production and carbon flow through trophic transfer to cooperation via symbioses. Increasingly, these communities are being understood through the lens of omics tools, which enable high-throughput processing of diverse communities. Metatranscriptomics offers an understanding of near real-time gene expression in microbial eukaryotic communities, providing a window into community metabolic activity.

RESULTS

Here we present a workflow for eukaryotic metatranscriptome assembly, and validate the ability of the pipeline to recapitulate real and manufactured eukaryotic community-level expression data. We also include an open-source tool for simulating environmental metatranscriptomes for testing and validation purposes. We reanalyze previously published metatranscriptomic datasets using our metatranscriptome analysis approach.

CONCLUSION

We determined that a multi-assembler approach improves eukaryotic metatranscriptome assembly based on recapitulated taxonomic and functional annotations from an in-silico mock community. The systematic validation of metatranscriptome assembly and annotation methods provided here is a necessary step to assess the fidelity of our community composition measurements and functional content assignments from eukaryotic metatranscriptomes.

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.

Expression of nuclear-encoded, haptophyte-derived ftsH genes support extremely rapid PSII repair and high-light photoacclimation in Karenia brevis (Dinophyceae)

Karenia brevis, a neurotoxic dinoflagellate that produces brevetoxins, is endemic to the Gulf of Mexico and can grow at high irradiances typical of surface waters found there. To build upon a growing number of studies addressing high-light tolerance in K. brevis, specific photobiology and molecular mechanisms underlying this capacity were evaluated in culture. Since photosystem II (PSII) repair cycle activity can be crucial to high light tolerance in plants and algae, the present study assessed this capacity in K. brevis and characterized the ftsH-like genes which are fundamental to this process. Compared with cultures grown in low-light, cultures grown in high-light showed a 65-fold increase in PSII photoinactivation, a ∼50-fold increase in PSII repair, enhanced nonphotochemical quenching (NPQ), and depressed F_v/F_m. Repair rates were among the fastest reported in phytoplankton. Publicly available K. brevis transcriptomes (MMETSP) were queried for ftsH-like sequences and refined with additional sequencing from two K. brevis strains. The genes were phylogenetically related to haptophyte orthologs, implicating acquisition during tertiary endosymbiosis. RT-qPCR of three of the four ftsH-like homologs revealed that poly-A tails predominated in all homologs, and that the most highly expressed homolog had a 5' splice leader and amino-acid motifs characteristic of chloroplast targeting, indicating nuclear encoding for this plastid-targeted gene. High-light cultures showed a ∼1.5-fold upregulation in mRNA expression of the thylakoid-associated genes. Overall, in conjunction with NPQ mechanisms, rapid PSII repair mediated by a haptophyte-derived ftsH prevents chronic photoinhibition in K. brevis. Our findings continue to build the case that high-light photobiology-supported by the acquisition and maintenance of tertiary endosymbiotic genes-is critical to the success of K. brevis in the Gulf of Mexico.

Proton-transporting heliorhodopsins from marine giant viruses

Rhodopsins convert light into signals and energy in animals and microbes. Heliorhodopsins (HeRs), a recently discovered new rhodopsin family, are widely present in archaea, bacteria, unicellular eukaryotes, and giant viruses, but their function remains unknown. Here, we report that a viral HeR from Emiliania huxleyi virus 202 (V2HeR3) is a light-activated proton transporter. V2HeR3 absorbs blue-green light, and the active intermediate contains the deprotonated retinal Schiff base. Site-directed mutagenesis study revealed that E191 in TM6 constitutes the gate together with the retinal Schiff base. E205 and E215 form a PAG of the Schiff base, and mutations at these positions converted the protein into an outward proton pump. Three environmental viral HeRs from the same group as well as a more distantly related HeR exhibited similar proton-transport activity, indicating that HeR functions might be diverse similarly to type-1 microbial rhodopsins. Some strains of E. huxleyi contain one HeR that is related to the viral HeRs, while its viruses EhV-201 and EhV-202 contain two and three HeRs, respectively. Except for V2HeR3 from EhV-202, none of these proteins exhibit ion transport activity. Thus, when expressed in the E. huxleyi cell membranes, only V2HeR3 has the potential to depolarize the host cells by light, possibly to overcome the host defense mechanisms or to prevent superinfection. The neuronal activity generated by V2HeR3 suggests that it can potentially be used as an optogenetic tool, similarly to type-1 microbial rhodopsins.

Microalgae-Based PUFAs for Food and Feed: Current Applications, Future Possibilities, and Constraints

Microalgae are currently considered an attractive source of highly valuable compounds for human and animal consumption, including polyunsaturated fatty acids (PUFAs). Several microalgae-derived compounds, such as ω-3 fatty acids, pigments, and whole dried biomasses are available on the market and are mainly produced by culturing microalgae in open ponds, which can be achieved with low setup and maintenance costs with respect to enclosed systems. However, open tanks are more susceptible to bacterial and other environmental contamination, do not guarantee a high reproducibility of algal biochemical profiles and productivities, and constrain massive cultivation to a limited number of species. Genetic engineering techniques have substantially improved over the last decade, and several model microalgae have been successfully modified to promote the accumulation of specific value-added compounds. However, transgenic strains should be cultured in closed photobioreactors (PBRs) to minimize risks of contamination of aquatic environments with allochthonous species; in addition, faster growth rates and higher yields of compounds of interest can be achieved in PBRs compared to open ponds. In this review, we present information collected about the major microalgae-derived commodities (with a special focus on PUFAs) produced at industrial scale, as well genetically-engineered microalgae to increase PUFA production. We also critically analyzed the main bottlenecks that make large-scale production of algal commodities difficult, as well as possible solutions to overcome the main problems and render the processes economically and environmentally safe.

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.

Reduced H+ channel activity disrupts pH homeostasis and calcification in coccolithophores at low ocean pH

Coccolithophore calcification is a major ocean biogeochemical process. While this process is likely to be sensitive to acidification-driven changes in ocean carbonate chemistry, incomplete understanding of the underlying mechanisms and constraints is a major bottleneck in predicting ocean acidification effects on calcification. We report severe disruption of pH homeostasis linked to a loss of H⁺ channel function in the coccolithophore Coccolithus braarudii acclimated to seawater pH values that are likely to be encountered currently in localized regions and more widely in future oceans. This disruption leads to specific defects in coccolith morphology. These findings provide mechanistic insight into how calcification in different coccolithophores is affected by changes in seawater carbonate chemistry.

Coccolithophores are major producers of ocean biogenic calcite, but this process is predicted to be negatively affected by future ocean acidification scenarios. Since coccolithophores calcify intracellularly, the mechanisms through which changes in seawater carbonate chemistry affect calcification remain unclear. Here we show that voltage-gated H⁺ channels in the plasma membrane of Coccolithus braarudii serve to regulate pH and maintain calcification under normal conditions but have greatly reduced activity in cells acclimated to low pH. This disrupts intracellular pH homeostasis and impairs the ability of C. braarudii to remove H⁺ generated by the calcification process, leading to specific coccolith malformations. These coccolith malformations can be reproduced by pharmacological inhibition of H⁺ channels. Heavily calcified coccolithophore species such as C. braarudii, which make the major contribution to carbonate export to the deep ocean, have a large intracellular H⁺ load and are likely to be most vulnerable to future decreases in ocean pH.

Genome evolution of a nonparasitic secondary heterotroph, the diatom Nitzschia putrida

Secondary loss of photosynthesis is observed across almost all plastid-bearing branches of the eukaryotic tree of life. However, genome-based insights into the transition from a phototroph into a secondary heterotroph have so far only been revealed for parasitic species. Free-living organisms can yield unique insights into the evolutionary consequence of the loss of photosynthesis, as the parasitic lifestyle requires specific adaptations to host environments. Here, we report on the diploid genome of the free-living diatom Nitzschia putrida (35 Mbp), a nonphotosynthetic osmotroph whose photosynthetic relatives contribute ca. 40% of net oceanic primary production. Comparative analyses with photosynthetic diatoms and heterotrophic algae with parasitic lifestyle revealed that a combination of gene loss, the accumulation of genes involved in organic carbon degradation, a unique secretome, and the rapid divergence of conserved gene families involved in cell wall and extracellular metabolism appear to have facilitated the lifestyle of a free-living secondary heterotroph.

Bioinformatic prediction of putative metallothioneins in non-ciliate protists

Intracellular ligands that bind heavy metals (HMs) and thereby minimize their detrimental effects to cellular metabolism are attracting great interest for a number of applications including bioremediation and development of HM-biosensors. Metallothioneins (MTs) are short, cysteine-rich, genetically encoded proteins involved in intracellular metal-binding and play a key role in detoxification of HMs. We searched approximately 700 genomes and transcriptomes of non-ciliate protists for novel putative MTs by similarity and structural analyses and found 21 unique proteins playing a potential role as MTs. Most putative MTs derive from heterokonts and dinoflagellates and share common features such as (i) a putative metal-binding domain in proximity of the N-terminus, (ii) two putative MT-specific domains near the C-terminus and (iii) one to three CTCGXXCXCGXXCXCXXC patterns. Although the biological function of these proteins has not been experimentally proven, knowledge of their genetic sequences adds useful information on proteins that are potentially involved in HM-binding and can contribute to the design of future biomolecular assays on HM-microbe interactions and MT-based biosensors.

The dynamic trophic architecture of open-ocean protist communities revealed through machine-guided metatranscriptomics

Mixotrophy is a ubiquitous nutritional strategy in marine ecosystems. Although our understanding of the distribution and abundance of mixotrophic plankton has improved significantly, the functional roles of mixotrophs are difficult to pinpoint, as mixotroph nutritional strategies are flexible and form a continuum between heterotrophy and phototrophy. We developed a machine learning–based method to assess the nutritional strategies of in situ planktonic populations based on metatranscriptomic profiles. We demonstrate that mixotrophic populations play varying functional roles along physicochemical gradients in the North Pacific Ocean, revealing a degree of physiological plasticity unique to aquatic mixotrophs. Our results highlight mechanisms that may dictate the flow of biogeochemical elements and ecology of the North Pacific Ocean, one of Earth's largest biogeographical provinces.

Intricate networks of single-celled eukaryotes (protists) dominate carbon flow in the ocean. Their growth, demise, and interactions with other microorganisms drive the fluxes of biogeochemical elements through marine ecosystems. Mixotrophic protists are capable of both photosynthesis and ingestion of prey and are dominant components of open-ocean planktonic communities. Yet the role of mixotrophs in elemental cycling is obscured by their capacity to act as primary producers or heterotrophic consumers depending on factors that remain largely uncharacterized. Here, we develop and apply a machine learning model that predicts the in situ trophic mode of aquatic protists based on their patterns of gene expression. This approach leverages a public collection of protist transcriptomes as a training set to identify a subset of gene families whose transcriptional profiles predict trophic mode. We applied our model to nearly 100 metatranscriptomes obtained during two oceanographic cruises in the North Pacific and found community-level and population-specific evidence that abundant open-ocean mixotrophic populations shift their predominant mode of nutrient and carbon acquisition in response to natural gradients in nutrient supply and sea surface temperature. Metatranscriptomic data from ship-board incubation experiments revealed that abundant mixotrophic prymnesiophytes from the oligotrophic North Pacific subtropical gyre rapidly remodeled their transcriptome to enhance photosynthesis when supplied with limiting nutrients. Coupling this approach with experiments designed to reveal the mechanisms driving mixotroph physiology provides an avenue toward understanding the ecology of mixotrophy in the natural environment.

Codon Usage Bias in Phytoplankton

Non-random usage of synonymous codons, known as “codon bias”, has been described in many organisms, from bacteria to Drosophila, but little is known about it in phytoplankton. This phenomenon is thought to be driven by selection for translational efficiency. As the efficacy of selection is proportional to the effective population size, species with large population sizes, such as phytoplankton, are expected to have strong codon bias. To test this, we measured codon bias in 215 strains from Haptophyta, Chlorophyta, Ochrophyta (except diatoms that were studied previously), Dinophyta, Cryptophyta, Ciliophora, unicellular Rhodophyta and Chlorarachniophyta. Codon bias is modest in most groups, despite the astronomically large population sizes of marine phytoplankton. The strength of the codon bias, measured with the effective number of codons, is the strongest in Haptophyta and the weakest in Chlorarachniophyta. The optimal codons are GC-ending in most cases, but several shifts to AT-ending codons were observed (mainly in Ochrophyta and Ciliophora). As it takes a long time to reach a new equilibrium after such shifts, species having AT-ending codons show a lower frequency of optimal codons compared to other species. Genetic diversity, calculated for species with more than three strains sequenced, is modest, indicating that the effective population sizes are many orders of magnitude lower than the astronomically large census population sizes, which helps to explain the modest codon bias in marine phytoplankton. This study represents the first comparative analysis of codon bias across multiple major phytoplankton groups.

Evolutionary Genomics of Sex-Related Chromosomes at the Base of the Green Lineage

Although sex is now accepted as a ubiquitous and ancestral feature of eukaryotes, direct observation of sex is still lacking in most unicellular eukaryotic lineages. Evidence of sex is frequently indirect and inferred from the identification of genes involved in meiosis from whole genome data and/or the detection of recombination signatures from genetic diversity in natural populations. In haploid unicellular eukaryotes, sex-related chromosomes are named mating-type (MTs) chromosomes and generally carry large genomic regions where recombination is suppressed. These regions have been characterized in Fungi and Chlorophyta and determine gamete compatibility and fusion. Two candidate MT+ and MT− alleles, spanning 450–650 kb, have recently been described in Ostreococcus tauri, a marine phytoplanktonic alga from the Mamiellophyceae class, an early diverging branch in the green lineage. Here, we investigate the architecture and evolution of these candidate MT+ and MT− alleles. We analyzed the phylogenetic profile and GC content of MT gene families in eight different genomes whose divergence has been previously estimated at up to 640 Myr, and found evidence that the divergence of the two MT alleles predates speciation in the Ostreococcus genus. Phylogenetic profiles of MT trans-specific polymorphisms in gametologs disclosed candidate MTs in two additional species, and possibly a third. These Mamiellales MT candidates are likely to be the oldest mating-type loci described to date, which makes them fascinating models to investigate the evolutionary mechanisms of haploid sex determination in eukaryotes.

Revealing RNA virus diversity and evolution in unicellular algae transcriptomes

Remarkably little is known about the diversity and evolution of RNA viruses in unicellular eukaryotes. We screened a total of 570 transcriptomes from the Marine Microbial Eukaryote Transcriptome Sequencing Project that encompasses a wide diversity of microbial eukaryotes, including most major photosynthetic lineages (i.e. the microalgae). From this, we identified thirty new and divergent RNA virus species, occupying a range of phylogenetic positions within the overall diversity of RNA viruses. Approximately one-third of the newly described viruses comprised single-stranded positive-sense RNA viruses from the order Lenarviricota associated with fungi, plants, and protists, while another third were related to the order Ghabrivirales, including members of the protist and fungi-associated Totiviridae. Other viral species showed sequence similarity to positive-sense RNA viruses from the algae-associated Marnaviridae, the double-stranded RNA (ds-RNA) Partitiviridae, as well as tentative evidence for one negative-sense RNA virus related to the Qinviridae. Importantly, we were able to identify divergent RNA viruses from distant host taxa, revealing the ancestry of these viral families and greatly extending our knowledge of the RNA viromes of microalgal cultures. Both the limited number of viruses detected per sample and the low sequence identity to known RNA viruses imply that additional microalgal viruses exist that could not be detected at the current sequencing depth or were too divergent to be identified using sequence similarity. Together, these results highlight the need for further investigation of algal-associated RNA viruses as well as the development of new tools to identify RNA viruses that exhibit very high levels of sequence divergence.

Decontamination, pooling and dereplication of the 678 samples of the Marine Microbial Eukaryote Transcriptome Sequencing Project

OBJECTIVES

Complex algae are photosynthetic organisms resulting from eukaryote-to-eukaryote endosymbiotic-like interactions. Yet the specific lineages and mechanisms are still under debate. That is why large scale phylogenomic studies are needed. Whereas available proteomes provide a limited diversity of complex algae, MMETSP (Marine Microbial Eukaryote Transcriptome Sequencing Project) transcriptomes represent a valuable resource for phylogenomic analyses, owing to their broad and rich taxonomic sampling, especially of photosynthetic species. Unfortunately, this sampling is unbalanced and sometimes highly redundant. Moreover, we observed contaminated sequences in some samples. In such a context, tree inference and readability are impaired. Consequently, the aim of the data processing reported here is to release a unique set of clean and non-redundant transcriptomes produced through an original protocol featuring decontamination, pooling and dereplication steps.

DATA DESCRIPTION

We submitted 678 MMETSP re-assembly samples to our parallel consolidation pipeline. Hence, we combined 423 samples into 110 consolidated transcriptomes, after the systematic removal of the most contaminated samples (186). This approach resulted in a total of 224 high-quality transcriptomes, easy to use and suitable to compute less contaminated, less redundant and more balanced phylogenies.

Olisthodiscus represents a new class of Ochrophyta

The phylogenetic diversity of Ochrophyta, a diverse and ecologically important radiation of algae, is still incompletely understood even at the level of the principal lineages. One taxon that has eluded simple classification is the marine flagellate genus Olisthodiscus. We investigated Olisthodiscus luteus K-0444 and documented its morphological and genetic differences from the NIES-15 strain, which we described as Olisthodiscus tomasii sp. nov. Phylogenetic analyses of combined 18S and 28S rRNA sequences confirmed that Olisthodiscus constitutes a separate, deep, ochrophyte lineage, but its position could not be resolved. To overcome this problem, we sequenced the plastid genome of O. luteus K-0444 and used the new data in multigene phylogenetic analyses, which suggested that Olisthodiscus is a sister lineage of the class Pinguiophyceae within a broader clade additionally including Chrysophyceae, Synchromophyceae, and Eustigmatophyceae. Surprisingly, the Olisthodiscus plastid genome contained three genes, ycf80, cysT, and cysW, inherited from the rhodophyte ancestor of the ochrophyte plastid yet lost from all other ochrophyte groups studied so far. Combined with nuclear genes for CysA and Sbp proteins, Olisthodiscus is the only known ochrophyte possessing a plastidial sulfate transporter SulT. In addition, the finding of a cemA gene in the Olisthodiscus plastid genome and an updated phylogenetic analysis ruled out the previously proposed hypothesis invoking horizontal cemA transfer from a green algal plastid into Synurales. Altogether, Olisthodiscus clearly represents a novel phylogenetically distinct ochrophyte lineage, which we have proposed as a new class, Olisthodiscophyceae.

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.

Experimental identification and in silico prediction of bacterivory in green algae

While algal phago-mixotrophs play a major role in aquatic microbial food webs, their diversity remains poorly understood. Recent studies have indicated several species of prasinophytes, early diverging green algae, to be able to consume bacteria for nutrition. To further explore the occurrence of phago-mixotrophy in green algae, we conducted feeding experiments with live fluorescently labeled bacteria stained with CellTracker Green CMFDA, heat-killed bacteria stained with 5-(4,6-dichlorotriazin-2-yl) aminofluorescein (DTAF), and magnetic beads. Feeding was detected via microscopy and/or flow cytometry in five strains of prasinophytes when provided with live bacteria: Pterosperma cristatum NIES626, Pyramimonas parkeae CCMP726, Pyramimonas parkeae NIES254, Nephroselmis pyriformis RCC618, and Dolichomastix tenuilepis CCMP3274. No feeding was detected when heat-killed bacteria or magnetic beads were provided, suggesting a strong preference for live prey in the strains tested. In parallel to experimental assays, green algal bacterivory was investigated using a gene-based prediction model. The predictions agreed with the experimental results and suggested bacterivory potential in additional green algae. Our observations underline the likelihood of widespread occurrence of phago-mixotrophy among green algae, while additionally highlighting potential biases introduced when using prey proxy to evaluate bacterial ingestion by algal cells.

Iron metabolism strategies in diatoms

Diatoms are one of the most successful group of photosynthetic eukaryotes in the contemporary ocean. They are ubiquitously distributed and are the most abundant primary producers in polar waters. Equally remarkable is their ability to tolerate iron deprivation and respond to periodic iron fertilization. Despite their relatively large cell sizes, diatoms tolerate iron limitation and frequently dominate iron-stimulated phytoplankton blooms, both natural and artificial. Here, we review the main iron use strategies of diatoms, including their ability to assimilate and store a range of iron sources, and the adaptations of their photosynthetic machinery and architecture to iron deprivation. Our synthesis relies on published literature and is complemented by a search of 82 diatom transcriptomes, including information collected from seven representatives of the most abundant diatom genera in the world's oceans.

Dinoflagellates alter their carbon and nutrient metabolic strategies across environmental gradients in the central Pacific Ocean

Marine microeukaryotes play a fundamental role in biogeochemical cycling through the transfer of energy to higher trophic levels and vertical carbon transport. Despite their global importance, microeukaryote physiology, nutrient metabolism and contributions to carbon cycling across offshore ecosystems are poorly characterized. Here, we observed the prevalence of dinoflagellates along a 4,600-km meridional transect extending across the central Pacific Ocean, where oligotrophic gyres meet equatorial upwelling waters rich in macronutrients yet low in dissolved iron. A combined multi-omics and geochemical analysis provided a window into dinoflagellate metabolism across the transect, indicating a continuous taxonomic dinoflagellate community that shifted its functional transcriptome and proteome as it extended from the euphotic to the mesopelagic zone. In euphotic waters, multi-omics data suggested that a combination of trophic modes were utilized, while mesopelagic metabolism was marked by cytoskeletal investments and nutrient recycling. Rearrangement in nutrient metabolism was evident in response to variable nitrogen and iron regimes across the gradient, with no associated change in community assemblage. Total dinoflagellate proteins scaled with particulate carbon export, with both elevated in equatorial waters, suggesting a link between dinoflagellate abundance and total carbon flux. Dinoflagellates employ numerous metabolic strategies that enable broad occupation of central Pacific ecosystems and play a dual role in carbon transformation through both photosynthetic fixation in the euphotic zone and remineralization in the mesopelagic zone.

PhycoCosm, a comparative algal genomics resource

Algae are a diverse, polyphyletic group of photosynthetic eukaryotes spanning nearly all eukaryotic lineages of life and collectively responsible for ∼50% of photosynthesis on Earth. Sequenced algal genomes, critical to understanding their complex biology, are growing in number and require efficient tools for analysis. PhycoCosm (https://phycocosm.jgi.doe.gov) is an algal multi-omics portal, developed by the US Department of Energy Joint Genome Institute to support analysis and distribution of algal genome sequences and other ‘omics’ data. PhycoCosm provides integration of genome sequence and annotation for >100 algal genomes with available multi-omics data and interactive web-based tools to enable algal research in bioenergy and the environment, encouraging community engagement and data exchange, and fostering new sequencing projects that will further these research goals.

Shift and Metabolic Potentials of Microbial Eukaryotic Communities Across the Full Depths of the Mariana Trench

Microbial eukaryotes are widespread and play important roles in marine ecosystems. However, their ecological characteristics in the deep sea (>1,000 m), especially hadal trenches, were largely unknown. Here, we investigated the diversity and metabolic potentials of microbial eukaryotes along the whole water column of the Mariana Trench by metagenomics. Our results showed clear depth-related distribution of microbial eukaryotic community and associated metabolic potentials. Surface seawater was dominated by phototrophic/mixotrophic groups (e.g., Dinoflagellata) and genes involved in biosynthesis (photosynthesis and fatty acid biosynthesis), while deep (bathypelagic and/or hadal) seawaters were enriched with heterotrophic groups (e.g., Bicoecea) and genes related to digestion (lysosomal enzymes and V-type ATPase) and carbohydrate metabolism. Co-occurrence analysis revealed high intra-domain connectivity, indicating that microbial eukaryotic composition was more influenced by microbial eukaryotes themselves than bacteria. Increased abundance of genes associated with unsaturated fatty acid biosynthesis likely plays a role in resisting high hydrostatic pressure. Top1 and hupB genes, responsible for the formation and stabilization of DNA structure, were unique and abundant in the hadal zone and thus may be helpful to stabilize DNA structure in the deep sea. Overall, our results provide insights into the distribution and potential adaptability of microbial eukaryotes in the hadal zone.

Streamlining data-intensive biology with workflow systems

As the scale of biological data generation has increased, the bottleneck of research has shifted from data generation to analysis. Researchers commonly need to build computational workflows that include multiple analytic tools and require incremental development as experimental insights demand tool and parameter modifications. These workflows can produce hundreds to thousands of intermediate files and results that must be integrated for biological insight. Data-centric workflow systems that internally manage computational resources, software, and conditional execution of analysis steps are reshaping the landscape of biological data analysis and empowering researchers to conduct reproducible analyses at scale. Adoption of these tools can facilitate and expedite robust data analysis, but knowledge of these techniques is still lacking. Here, we provide a series of strategies for leveraging workflow systems with structured project, data, and resource management to streamline large-scale biological analysis. We present these practices in the context of high-throughput sequencing data analysis, but the principles are broadly applicable to biologists working beyond this field.

Feeding diverse prey as an excellent strategy of mixotrophic dinoflagellates for global dominance

Microalgae fuel food webs and biogeochemical cycles of key elements in the ocean. What determines microalgal dominance in the ocean is a long-standing question. Red tide distribution data (spanning 1990 to 2019) show that mixotrophic dinoflagellates, capable of photosynthesis and predation together, were responsible for ~40% of the species forming red tides globally. Counterintuitively, the species with low or moderate growth rates but diverse prey including diatoms caused red tides globally. The ability of these dinoflagellates to trade off growth for prey diversity is another genetic factor critical to formation of red tides across diverse ocean conditions. This finding has profound implications for explaining the global dominance of particular microalgae, their key eco-evolutionary strategy, and prediction of harmful red tide outbreaks.

Overexpression of Key Sterol Pathway Enzymes in Two Model Marine Diatoms Alters Sterol Profiles in Phaeodactylum tricornutum

Sterols are a class of triterpenoid molecules with diverse functional roles in eukaryotic cells, including intracellular signaling and regulation of cell membrane fluidity. Diatoms are a dominant eukaryotic phytoplankton group that produce a wide diversity of sterol compounds. The enzymes 3-hydroxy-3-methyl glutaryl CoA reductase (HMGR) and squalene epoxidase (SQE) have been reported to be rate-limiting steps in sterol biosynthesis in other model eukaryotes; however, the extent to which these enzymes regulate triterpenoid production in diatoms is not known. To probe the role of these two metabolic nodes in the regulation of sterol metabolic flux in diatoms, we independently over-expressed two versions of the native HMGR and a conventional, heterologous SQE gene in the diatoms Thalassiosira pseudonana and Phaeodactylum tricornutum. Overexpression of these key enzymes resulted in significant differential accumulation of downstream sterol pathway intermediates in P. tricornutum. HMGR-mVenus overexpression resulted in the accumulation of squalene, cycloartenol, and obtusifoliol, while cycloartenol and obtusifoliol accumulated in response to heterologous NoSQE-mVenus overexpression. In addition, accumulation of the end-point sterol 24-methylenecholesta-5,24(24’)-dien-3β-ol was observed in all P. tricornutum overexpression lines, and campesterol increased three-fold in P. tricornutum lines expressing NoSQE-mVenus. Minor differences in end-point sterol composition were also found in T. pseudonana, but no accumulation of sterol pathway intermediates was observed. Despite the successful manipulation of pathway intermediates and individual sterols in P. tricornutum, total sterol levels did not change significantly in transformed lines, suggesting the existence of tight pathway regulation to maintain total sterol content.

Comparative Plastid Genomics of Cryptomonas Species Reveals Fine-Scale Genomic Responses to Loss of Photosynthesis

Loss of photosynthesis is a recurring theme in eukaryotic evolution. In organisms that have lost the ability to photosynthesize, nonphotosynthetic plastids are retained because they play essential roles in processes other than photosynthesis. The unicellular algal genus Cryptomonas contains both photosynthetic and nonphotosynthetic members, the latter having lost the ability to photosynthesize on at least three separate occasions. To elucidate the evolutionary processes underlying the loss of photosynthesis, we sequenced the plastid genomes of two nonphotosynthetic strains, Cryptomonas sp. CCAC1634B and SAG977-2f, as well as the genome of the phototroph Cryptomonas curvata CCAP979/52. These three genome sequences were compared with the previously sequenced plastid genome of the nonphotosynthetic species Cryptomonas paramecium CCAP977/2a as well as photosynthetic members of the Cryptomonadales, including C. curvata FBCC300012D. Intraspecies comparison between the two C. curvata strains showed that although their genome structures are stable, the substitution rates of their genes are relatively high. Although most photosynthesis-related genes, such as the psa and psb gene families, were found to have disappeared from the nonphotosynthetic strains, at least ten pseudogenes are retained in SAG977-2f. Although gene order is roughly shared among the plastid genomes of photosynthetic Cryptomonadales, genome rearrangements are seen more frequently in the smaller genomes of the nonphotosynthetic strains. Intriguingly, the light-independent protochlorophyllide reductase comprising chlB, L, and N is retained in nonphotosynthetic SAG977-2f and CCAC1634B. On the other hand, whereas CCAP977/2a retains ribulose-1,5-bisphosphate carboxylase/oxygenase-related genes, including rbcL, rbcS, and cbbX, the plastid genomes of the other two nonphotosynthetic strains have lost the ribulose-1,5-bisphosphate carboxylase/oxygenase protein-coding genes.

A Novel Single-Domain Na+-Selective Voltage-Gated Channel in Photosynthetic Eukaryotes

The evolution of Na⁺-selective four-domain voltage-gated channels (4D-Na_vs) in animals allowed rapid Na⁺-dependent electrical excitability, and enabled the development of sophisticated systems for rapid and long-range signaling. While bacteria encode single-domain Na⁺-selective voltage-gated channels (BacNa_v), they typically exhibit much slower kinetics than 4D-Na_vs, and are not thought to have crossed the prokaryote-eukaryote boundary. As such, the capacity for rapid Na⁺-selective signaling is considered to be confined to certain animal taxa, and absent from photosynthetic eukaryotes. Certainly, in land plants, such as the Venus flytrap (Dionaea muscipula) where fast electrical excitability has been described, this is most likely based on fast anion channels. Here, we report a unique class of eukaryotic Na⁺-selective, single-domain channels (EukCatBs) that are present primarily in haptophyte algae, including the ecologically important calcifying coccolithophores, Emiliania huxleyi and Scyphosphaera apsteinii The EukCatB channels exhibit very rapid voltage-dependent activation and inactivation kinetics, and isoform-specific sensitivity to the highly selective 4D-Na_v blocker tetrodotoxin. The results demonstrate that the capacity for rapid Na⁺-based signaling in eukaryotes is not restricted to animals or to the presence of 4D-Na_vs. The EukCatB channels therefore represent an independent evolution of fast Na⁺-based electrical signaling in eukaryotes that likely contribute to sophisticated cellular control mechanisms operating on very short time scales in unicellular algae.

MetaEuk-sensitive, high-throughput gene discovery, and annotation for large-scale eukaryotic metagenomics

BACKGROUND

Metagenomics is revolutionizing the study of microorganisms and their involvement in biological, biomedical, and geochemical processes, allowing us to investigate by direct sequencing a tremendous diversity of organisms without the need for prior cultivation. Unicellular eukaryotes play essential roles in most microbial communities as chief predators, decomposers, phototrophs, bacterial hosts, symbionts, and parasites to plants and animals. Investigating their roles is therefore of great interest to ecology, biotechnology, human health, and evolution. However, the generally lower sequencing coverage, their more complex gene and genome architectures, and a lack of eukaryote-specific experimental and computational procedures have kept them on the sidelines of metagenomics.

RESULTS

MetaEuk is a toolkit for high-throughput, reference-based discovery, and annotation of protein-coding genes in eukaryotic metagenomic contigs. It performs fast searches with 6-frame-translated fragments covering all possible exons and optimally combines matches into multi-exon proteins. We used a benchmark of seven diverse, annotated genomes to show that MetaEuk is highly sensitive even under conditions of low sequence similarity to the reference database. To demonstrate MetaEuk's power to discover novel eukaryotic proteins in large-scale metagenomic data, we assembled contigs from 912 samples of the Tara Oceans project. MetaEuk predicted >12,000,000 protein-coding genes in 8 days on ten 16-core servers. Most of the discovered proteins are highly diverged from known proteins and originate from very sparsely sampled eukaryotic supergroups.

CONCLUSION

The open-source (GPLv3) MetaEuk software (https://github.com/soedinglab/metaeuk) enables large-scale eukaryotic metagenomics through reference-based, sensitive taxonomic and functional annotation. Video abstract.

Phytochrome evolution in 3D: deletion, duplication, and diversification

Canonical plant phytochromes are master regulators of photomorphogenesis and the shade avoidance response. They are also part of a widespread superfamily of photoreceptors with diverse spectral and biochemical properties. Plant phytochromes belong to a clade including other phytochromes from glaucophyte, prasinophyte, and streptophyte algae (all members of the Archaeplastida) and those from cryptophyte algae. This is consistent with recent analyses supporting the existence of an AC (Archaeplastida + Cryptista) clade. AC phytochromes have been proposed to arise from ancestral cyanobacterial genes via endosymbiotic gene transfer (EGT), but most recent studies instead support multiple horizontal gene transfer (HGT) events to generate extant eukaryotic phytochromes. In principle, this scenario would be compared to the emerging understanding of early events in eukaryotic evolution to generate a coherent picture. Unfortunately, there is currently a major discrepancy between the evolution of phytochromes and the evolution of eukaryotes; phytochrome evolution is thus not a solved problem. We therefore examine phytochrome evolution in a broader context. Within this context, we can identify three important themes in phytochrome evolution: deletion, duplication, and diversification. These themes drive phytochrome evolution as organisms evolve in response to environmental challenges.

Dinoflagellates with relic endosymbiont nuclei as models for elucidating organellogenesis

Nucleomorphs are relic endosymbiont nuclei so far found only in two algal groups, cryptophytes and chlorarachniophytes, which have been studied to model the evolutionary process of integrating an endosymbiont alga into a host-governed plastid (organellogenesis). However, past studies suggest that DNA transfer from the endosymbiont to host nuclei had already ceased in both cryptophytes and chlorarachniophytes, implying that the organellogenesis at the genetic level has been completed in the two systems. Moreover, we have yet to pinpoint the closest free-living relative of the endosymbiotic alga engulfed by the ancestral chlorarachniophyte or cryptophyte, making it difficult to infer how organellogenesis altered the endosymbiont genome. To counter the above issues, we need novel nucleomorph-bearing algae, in which endosymbiont-to-host DNA transfer is on-going and for which endosymbiont/plastid origins can be inferred at a fine taxonomic scale. Here, we report two previously undescribed dinoflagellates, strains MGD and TGD, with green algal endosymbionts enclosing plastids as well as relic nuclei (nucleomorphs). We provide evidence for the presence of DNA in the two nucleomorphs and the transfer of endosymbiont genes to the host (dinoflagellate) genomes. Furthermore, DNA transfer between the host and endosymbiont nuclei was found to be in progress in both the MGD and TGD systems. Phylogenetic analyses successfully resolved the origins of the endosymbionts at the genus level. With the combined evidence, we conclude that the host-endosymbiont integration in MGD/TGD is less advanced than that in cryptophytes/chrorarachniophytes, and propose the two dinoflagellates as models for elucidating organellogenesis.

SAGER: a database of Symbiodiniaceae and Algal Genomic Resource

Symbiodiniaceae dinoflagellates are essential endosymbionts of reef building corals and some other invertebrates. Information of their genome structure and function is critical for understanding coral symbiosis and bleaching. With the rapid development of sequencing technology, genome draft assemblies of several Symbiodiniaceae species and diverse marine algal genomes have become publicly available but spread in multiple separate locations. Here, we present a Symbiodiniaceae and Algal Genomic Resource Database (SAGER), a user-friendly online repository for integrating existing genomic data of Symbiodiniaceae species and diverse marine algal gene sets from MMETSP and PhyloDB databases. Relevant algal data are included to facilitate comparative analyses. The database is freely accessible at http://sampgr.org.cn. It provides comprehensive tools for studying gene function, expression and comparative genomics, including search tools to identify gene information from Symbiodiniaceae species, and BLAST tool to find orthologs from marine algae and protists. Moreover, SAGER integrates transcriptome datasets derived from diverse culture conditions of corresponding Symbiodiniaceae species. SAGER was developed with the capacity to incorporate future Symbiodiniaceae and algal genome and transcriptome data, and will serve as an open-access and sustained platform providing genomic and molecular tools that can be conveniently used to study Symbiodiniaceae and other marine algae. Database URL: http://sampgr.org.cn.

Expansion of phycobilisome linker gene families in mesophilic red algae

The common ancestor of red algae (Rhodophyta) has undergone massive genome reduction, whereby 25% of the gene inventory has been lost, followed by its split into the species-poor extremophilic Cyanidiophytina and the broadly distributed mesophilic red algae. Success of the mesophile radiation is surprising given their highly reduced gene inventory. To address this latter issue, we combine an improved genome assembly from the unicellular red alga Porphyridium purpureum with a diverse collection of other algal genomes to reconstruct ancient endosymbiotic gene transfers (EGTs) and gene duplications. We find EGTs associated with the core photosynthetic machinery that may have played important roles in plastid establishment. More significant are the extensive duplications and diversification of nuclear gene families encoding phycobilisome linker proteins that stabilize light-harvesting functions. We speculate that the origin of these complex families in mesophilic red algae may have contributed to their adaptation to a diversity of light environments.

Widely distributed red algae have experienced massive genome reduction during evolution. Here, using an improved genome assembly of Porphyridium purpureum, Lee et al. show the role of endosymbiotic gene transfer in plastid evolution and the correlation between phycobilisome linker diversification and the red algal radiation.

Extensive Reduction of the Nuclear Pore Complex in Nucleomorphs

The nuclear pore complex (NPC) is a large macromolecular assembly situated within the pores of the nuclear envelope. Through interactions between its subcomplexes and import proteins, the NPC mediates the transport of molecules into and out of the nucleus and facilitates dynamic chromatin regulation and gene expression. Accordingly, the NPC constitutes a highly integrated nuclear component that is ubiquitous and conserved among eukaryotes. Potential exceptions to this are nucleomorphs: Highly reduced, relict nuclei that were derived from green and red algae following their endosymbiotic integration into two lineages, the chlorarachniophytes and the cryptophyceans. A previous investigation failed to identify NPC genes in nucleomorph genomes suggesting that these genes have either been relocated to the host nucleus or lost. Here, we sought to investigate the composition of the NPC in nucleomorphs by using genomic and transcriptomic data to identify and phylogenetically classify NPC proteins in nucleomorph-containing algae. Although we found NPC proteins in all examined lineages, most of those found in chlorarachniophytes and cryptophyceans were single copy, host-related proteins that lacked signal peptides. Two exceptions were Nup98 and Rae1, which had clear nucleomorph-derived homologs. However, these proteins alone are likely insufficient to structure a canonical NPC and previous reports revealed that Nup98 and Rae1 have other nuclear functions. Ultimately, these data indicate that nucleomorphs represent eukaryotic nuclei without a canonical NPC, raising fundamental questions about their structure and function.

Re-assembly, quality evaluation, and annotation of 678 microbial eukaryotic reference transcriptomes

BACKGROUND

De novo transcriptome assemblies are required prior to analyzing RNA sequencing data from a species without an existing reference genome or transcriptome. Despite the prevalence of transcriptomic studies, the effects of using different workflows, or "pipelines," on the resulting assemblies are poorly understood. Here, a pipeline was programmatically automated and used to assemble and annotate raw transcriptomic short-read data collected as part of the Marine Microbial Eukaryotic Transcriptome Sequencing Project. The resulting transcriptome assemblies were evaluated and compared against assemblies that were previously generated with a different pipeline developed by the National Center for Genome Research.

RESULTS

New transcriptome assemblies contained the majority of previous contigs as well as new content. On average, 7.8% of the annotated contigs in the new assemblies were novel gene names not found in the previous assemblies. Taxonomic trends were observed in the assembly metrics. Assemblies from the Dinoflagellata showed a higher number of contigs and unique k-mers than transcriptomes from other phyla, while assemblies from Ciliophora had a lower percentage of open reading frames compared to other phyla.

CONCLUSIONS

Given current bioinformatics approaches, there is no single "best" reference transcriptome for a particular set of raw data. As the optimum transcriptome is a moving target, improving (or not) with new tools and approaches, automated and programmable pipelines are invaluable for managing the computationally intensive tasks required for re-processing large sets of samples with revised pipelines and ensuring a common evaluation workflow is applied to all samples. Thus, re-assembling existing data with new tools using automated and programmable pipelines may yield more accurate identification of taxon-specific trends across samples in addition to novel and useful products for the community.

Keeping it light: (re)analyzing community-wide datasets without major infrastructure

DNA sequencing technology has revolutionized the field of biology, shifting biology from a data-limited to data-rich state. Central to the interpretation of sequencing data are the computational tools and approaches that convert raw data into biologically meaningful information. Both the tools and the generation of data are actively evolving, yet the practice of re-analysis of previously generated data with new tools is not commonplace. Re-analysis of existing data provides an affordable means of generating new information and will likely become more routine within biology, yet necessitates a new set of considerations for best practices and resource development. Here, we discuss several practices that we believe to be broadly applicable when re-analyzing data, especially when done by small research groups.