A global ocean atlas of eukaryotic genes

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.

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.

Self-assembling viral histones are evolutionary intermediates between archaeal and eukaryotic nucleosomes

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

Sensitive and error-tolerant annotation of protein-coding DNA with BATH

Summary

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based translated sequence annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long-read sequencing data and in the context of pseudogenes.

Availability and implementation

The software is available at https://github.com/TravisWheelerLab/BATH.

A refined picture of the native amine dehydrogenase family revealed by extensive biodiversity screening

Native amine dehydrogenases offer sustainable access to chiral amines, so the search for scaffolds capable of converting more diverse carbonyl compounds is required to reach the full potential of this alternative to conventional synthetic reductive aminations. Here we report a multidisciplinary strategy combining bioinformatics, chemoinformatics and biocatalysis to extensively screen billions of sequences in silico and to efficiently find native amine dehydrogenases features using computational approaches. In this way, we achieve a comprehensive overview of the initial native amine dehydrogenase family, extending it from 2,011 to 17,959 sequences, and identify native amine dehydrogenases with non-reported substrate spectra, including hindered carbonyls and ethyl ketones, and accepting methylamine and cyclopropylamine as amine donor. We also present preliminary model-based structural information to inform the design of potential (R)-selective amine dehydrogenases, as native amine dehydrogenases are mostly (S)-selective. This integrated strategy paves the way for expanding the resource of other enzyme families and in highlighting enzymes with original features.

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.

A distinct, high-affinity, alkaline phosphatase facilitates occupation of P-depleted environments by marine picocyanobacteria

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus, the two most abundant phototrophs on Earth, thrive in oligotrophic oceanic regions. While it is well known that specific lineages are exquisitely adapted to prevailing in situ light and temperature regimes, much less is known of the molecular machinery required to facilitate occupancy of these low-nutrient environments. Here, we describe a hitherto unknown alkaline phosphatase, Psip1, that has a substantially higher affinity for phosphomonoesters than other well-known phosphatases like PhoA, PhoX, or PhoD and is restricted to clade III Synechococcus and a subset of high light I-adapted Prochlorococcus strains, suggesting niche specificity. We demonstrate that Psip1 has undergone convergent evolution with PhoX, requiring both iron and calcium for activity and likely possessing identical key residues around the active site, despite generally very low sequence homology. Interrogation of metagenomes and transcriptomes from TARA oceans and an Atlantic Meridional transect shows that psip1 is abundant and highly expressed in picocyanobacterial populations from the Mediterranean Sea and north Atlantic gyre, regions well recognized to be phosphorus (P)-deplete. Together, this identifies psip1 as an important oligotrophy-specific gene for P recycling in these organisms. Furthermore, psip1 is not restricted to picocyanobacteria and is abundant and highly transcribed in some α-proteobacteria and eukaryotic algae, suggesting that such a high-affinity phosphatase is important across the microbial taxonomic world to occupy low-P environments.

Decoding populations in the ocean microbiome

Understanding the characteristics and structure of populations is fundamental to comprehending ecosystem processes and evolutionary adaptations. While the study of animal and plant populations has spanned a few centuries, microbial populations have been under scientific scrutiny for a considerably shorter period. In the ocean, analyzing the genetic composition of microbial populations and their adaptations to multiple niches can yield important insights into ecosystem function and the microbiome's response to global change. However, microbial populations have remained elusive to the scientific community due to the challenges associated with isolating microorganisms in the laboratory. Today, advancements in large-scale metagenomics and metatranscriptomics facilitate the investigation of populations from many uncultured microbial species directly from their habitats. The knowledge acquired thus far reveals substantial genetic diversity among various microbial species, showcasing distinct patterns of population differentiation and adaptations, and highlighting the significant role of selection in structuring populations. In the coming years, population genomics is expected to significantly increase our understanding of the architecture and functioning of the ocean microbiome, providing insights into its vulnerability or resilience in the face of ongoing global change. Video Abstract.

In situ community transcriptomics illuminates CO2-fixation potentials and supporting roles of phagotrophy and proton pump in plankton in a subtropical marginal sea

Lineage-wise physiological activities of plankton communities in the ocean are important but challenging to characterize. Here, we conducted whole-assemblage metatranscriptomic profiling at continental shelf and slope sites in the South China Sea to investigate carbon fixation potential in different lineages. RuBisCO expression, the proxy of Calvin carbon fixation (CCF) potential, was mainly contributed by Bacillariophyta, Chlorophyta, Cyanobacteria, and Haptophyta, which was differentially affected by environmental factors among lineages. CCF potential exhibited positive or negative correlations with phagotrophy gene expression, suggesting phagotrophy possibly enhances or complements CCF. Our data also reveal significant non-Calvin carbon fixation (NCF) potential, as indicated by the active expression of genes in all five currently recognized NCF pathways, mainly contributed by Flavobacteriales, Alteromonadales, and Oceanospirillales. Furthermore, in Flavobacteriales, Alteromonadales, Pelagibacterales, and Rhodobacterales, NCF potential was positively correlated with proton-pump rhodopsin (PPR) expression, suggesting that NCF might be energetically supported by PPR. The novel insights into the lineage-differential potential of carbon fixation, widespread mixotrophy, and PPR as an energy source for NCF lay a methodological and informational foundation for further research to understand carbon fixation and the trophic landscape in the ocean.IMPORTANCEMarine plankton plays an important role in global carbon cycling and climate regulation. Phytoplankton and cyanobacteria fix CO2 to produce organic compounds using solar energy and mainly by the Calvin cycle, whereas autotrophic bacteria and archaea may fix CO2 by non-Calvin cycle carbon fixation pathways. How active individual lineages are in carbon fixation and mixotrophy, and what energy source bacteria may employ in non-Calvin carbon fixation, in a natural plankton assemblage are poorly understood and underexplored. Using metatranscriptomics, we studied carbon fixation in marine plankton with lineage resolution in tropical marginal shelf and slope areas. Based on the sequencing results, we characterized the carbon fixation potential of different lineages and assessed Calvin- and non-Calvin- carbon fixation activities and energy sources. Data revealed a high number of unigenes (4.4 million), lineage-dependent differential potentials of Calvin carbon fixation and responses to environmental conditions, major contributors of non-Calvin carbon fixation, and their potential energy source.

Flexible B12 ecophysiology of Phaeocystis antarctica due to a fusion B12-independent methionine synthase with widespread homologues

Coastal Antarctic marine ecosystems are significant in carbon cycling because of their intense seasonal phytoplankton blooms. Southern Ocean algae are primarily limited by light and iron (Fe) and can be co-limited by cobalamin (vitamin B12). Micronutrient limitation controls productivity and shapes the composition of blooms which are typically dominated by either diatoms or the haptophyte Phaeocystis antarctica. However, the vitamin requirements and ecophysiology of the keystone species P. antarctica remain poorly characterized. Using cultures, physiological analysis, and comparative omics, we examined the response of P. antarctica to a matrix of Fe-B12 conditions. We show that P. antarctica is not auxotrophic for B12, as previously suggested, and identify mechanisms underlying its B12 response in cultures of predominantly solitary and colonial cells. A combination of proteomics and proteogenomics reveals a B12-independent methionine synthase fusion protein (MetE-fusion) that is expressed under vitamin limitation and interreplaced with the B12-dependent isoform under replete conditions. Database searches return homologues of the MetE-fusion protein in multiple Phaeocystis species and in a wide range of marine microbes, including other photosynthetic eukaryotes with polymorphic life cycles as well as bacterioplankton. Furthermore, we find MetE-fusion homologues expressed in metaproteomic and metatranscriptomic field samples in polar and more geographically widespread regions. As climate change impacts micronutrient availability in the coastal Southern Ocean, our finding that P. antarctica has a flexible B12 metabolism has implications for its relative fitness compared to B12-auxotrophic diatoms and for the detection of B12-stress in a more diverse set of marine microbes.

Evolution: No end in sight for novel incredible (heterotrophic) protists

The microbial eukaryotes known as protists are of immense importance for our understanding of eukaryotic biology. Although it is often difficult to convince funding bodies to sponsor research projects aimed at finding new protist lineages, such discoveries usually provide new and fundamental insights into cell and evolutionary biology, and ecology.

SLC24A-mediated calcium exchange as an indispensable component of the diatom cell density-driven signaling pathway

Diatom bloom is characterized by a rapid increase of population density. Perception of population density and physiological responses can significantly influence their survival strategies, subsequently impacting bloom fate. The population density itself can serve as a signal, which is perceived through chemical signals or chlorophyll fluorescence signals triggered by high cell density, and their intracellular signaling mechanisms remain to be elucidated. In this study, we focused on the model diatom, Phaeodactylum tricornutum, and designed an orthogonal experiment involving varying cell densities and light conditions, to stimulate the release of chemical signals and light-induced chlorophyll fluorescence signals. Utilizing RNA-Seq and Weighted Gene Co-expression Network Analysis, we identified four gene clusters displaying density-dependent expression patterns. Within these, a potential hub gene, PtSLC24A, encoding a Na+/Ca2+ exchanger, was identified. Based on molecular genetics, cellular physiology, computational structural biology, and in situ oceanic data, we propose a potential intracellular signaling mechanism related to cell density in marine diatoms using Ca2+: upon sensing population density signals mediated by chemical cues, the membrane-bound PtSLC24A facilitates the efflux of Ca2+ to maintain specific intracellular calcium levels, allowing the transduction of intracellular density signals, subsequently regulating physiological responses, including cell apoptosis, ultimately affecting algal blooms fate. These findings shed light on the calcium-mediated intracellular signaling mechanism of marine diatoms to changing population densities, and enhances our understanding of diatom bloom dynamics and their ecological implications.

Sensitive and error-tolerant annotation of protein-coding DNA with BATH

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long read sequencing data and in the context of pseudogenes.

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

The diversity, community dynamics, and interactions of the microbiome in the world's deepest blue hole: insights into extreme environmental response patterns and tolerance of marine microorganisms

IMPORTANCE

This study comprehensively examined the community dynamics, functional profiles, and interactions of the microbiome in the world's deepest blue hole. The findings revealed a positive correlation between the α-diversities of Symbiodiniaceae and archaea, indicating the potential reliance of Symbiodiniaceae on archaea in an extreme environment resulting from a partial niche overlap. The negative association between the α-diversity and β-diversity of the bacterial community suggested that the change rule of the bacterial community was consistent with the Anna Karenina effects. The core microbiome comprised nine microbial taxa, highlighting their remarkable tolerance and adaptability to sharp environmental gradient variations. Bacteria and archaea played significant roles in carbon, nitrogen, and sulfur cycles, while fungi contributed to carbon metabolism. This study advanced our understanding of the community dynamics, response patterns, and resilience of microorganisms populating the world's deepest blue hole, thereby facilitating further ecological and evolutional exploration of microbiomes in diverse extreme environments.

DMSOP-cleaving enzymes are diverse and widely distributed in marine microorganisms

Dimethylsulfoxonium propionate (DMSOP) is a recently identified and abundant marine organosulfur compound with roles in oxidative stress protection, global carbon and sulfur cycling and, as shown here, potentially in osmotolerance. Microbial DMSOP cleavage yields dimethyl sulfoxide, a ubiquitous marine metabolite, and acrylate, but the enzymes responsible, and their environmental importance, were unknown. Here we report DMSOP cleavage mechanisms in diverse heterotrophic bacteria, fungi and phototrophic algae not previously known to have this activity, and highlight the unappreciated importance of this process in marine sediment environments. These diverse organisms, including Roseobacter, SAR11 bacteria and Emiliania huxleyi, utilized their dimethylsulfoniopropionate lyase 'Ddd' or 'Alma' enzymes to cleave DMSOP via similar catalytic mechanisms to those for dimethylsulfoniopropionate. Given the annual teragram predictions for DMSOP production and its prevalence in marine sediments, our results highlight that DMSOP cleavage is likely a globally significant process influencing carbon and sulfur fluxes and ecological interactions.

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.

Eukaryotic genomes from a global metagenomic data set illuminate trophic modes and biogeography of ocean plankton

Metagenomics is a powerful method for interpreting the ecological roles and physiological capabilities of mixed microbial communities. Yet, many tools for processing metagenomic data are neither designed to consider eukaryotes nor are they built for an increasing amount of sequence data. EukHeist is an automated pipeline to retrieve eukaryotic and prokaryotic metagenome-assembled genomes (MAGs) from large-scale metagenomic sequence data sets. We developed the EukHeist workflow to specifically process large amounts of both metagenomic and/or metatranscriptomic sequence data in an automated and reproducible fashion. Here, we applied EukHeist to the large-size fraction data (0.8-2,000 µm) from Tara Oceans to recover both eukaryotic and prokaryotic MAGs, which we refer to as TOPAZ (Tara Oceans Particle-Associated MAGs). The TOPAZ MAGs consisted of >900 environmentally relevant eukaryotic MAGs and >4,000 bacterial and archaeal MAGs. The bacterial and archaeal TOPAZ MAGs expand upon the phylogenetic diversity of likely particle- and host-associated taxa. We use these MAGs to demonstrate an approach to infer the putative trophic mode of the recovered eukaryotic MAGs. We also identify ecological cohorts of co-occurring MAGs, which are driven by specific environmental factors and putative host-microbe associations. These data together add to a number of growing resources of environmentally relevant eukaryotic genomic information. Complementary and expanded databases of MAGs, such as those provided through scalable pipelines like EukHeist, stand to advance our understanding of eukaryotic diversity through increased coverage of genomic representatives across the tree of life.IMPORTANCESingle-celled eukaryotes play ecologically significant roles in the marine environment, yet fundamental questions about their biodiversity, ecological function, and interactions remain. Environmental sequencing enables researchers to document naturally occurring protistan communities, without culturing bias, yet metagenomic and metatranscriptomic sequencing approaches cannot separate individual species from communities. To more completely capture the genomic content of mixed protistan populations, we can create bins of sequences that represent the same organism (metagenome-assembled genomes [MAGs]). We developed the EukHeist pipeline, which automates the binning of population-level eukaryotic and prokaryotic genomes from metagenomic reads. We show exciting insight into what protistan communities are present and their trophic roles in the ocean. Scalable computational tools, like EukHeist, may accelerate the identification of meaningful genetic signatures from large data sets and complement researchers' efforts to leverage MAG databases for addressing ecological questions, resolving evolutionary relationships, and discovering potentially novel biodiversity.

Non-cyanobacterial diazotrophs: global diversity, distribution, ecophysiology, and activity in marine waters

Biological dinitrogen (N2) fixation supplies nitrogen to the oceans, supporting primary productivity, and is carried out by some bacteria and archaea referred to as diazotrophs. Cyanobacteria are conventionally considered to be the major contributors to marine N2 fixation, but non-cyanobacterial diazotrophs (NCDs) have been shown to be distributed throughout ocean ecosystems. However, the biogeochemical significance of marine NCDs has not been demonstrated. This review synthesizes multiple datasets, drawing from cultivation-independent molecular techniques and data from extensive oceanic expeditions, to provide a comprehensive view into the diversity, biogeography, ecophysiology, and activity of marine NCDs. A NCD nifH gene catalog was compiled containing sequences from both PCR-based and PCR-free methods, identifying taxa for future studies. NCD abundances from a novel database of NCD nifH-based abundances were colocalized with environmental data, unveiling distinct distributions and environmental drivers of individual taxa. Mechanisms that NCDs may use to fuel and regulate N2 fixation in response to oxygen and fixed nitrogen availability are discussed, based on a metabolic analysis of recently available Tara Oceans expedition data. The integration of multiple datasets provides a new perspective that enhances understanding of the biology, ecology, and biogeography of marine NCDs and provides tools and directions for future research.

Genomic adaptation of giant viruses in polar oceans

Despite being perennially frigid, polar oceans form an ecosystem hosting high and unique biodiversity. Various organisms show different adaptive strategies in this habitat, but how viruses adapt to this environment is largely unknown. Viruses of phyla Nucleocytoviricota and Mirusviricota are groups of eukaryote-infecting large and giant DNA viruses with genomes encoding a variety of functions. Here, by leveraging the Global Ocean Eukaryotic Viral database, we investigate the biogeography and functional repertoire of these viruses at a global scale. We first confirm the existence of an ecological barrier that clearly separates polar and nonpolar viral communities, and then demonstrate that temperature drives dramatic changes in the virus-host network at the polar-nonpolar boundary. Ancestral niche reconstruction suggests that adaptation of these viruses to polar conditions has occurred repeatedly over the course of evolution, with polar-adapted viruses in the modern ocean being scattered across their phylogeny. Numerous viral genes are specifically associated with polar adaptation, although most of their homologues are not identified as polar-adaptive genes in eukaryotes. These results suggest that giant viruses adapt to cold environments by changing their functional repertoire, and this viral evolutionary strategy is distinct from the polar adaptation strategy of their hosts.

A chlorophyll c synthase widely co-opted by phytoplankton

Marine and terrestrial photosynthesis exhibit a schism in the accessory chlorophyll (Chl) that complements the function of Chl a: Chl b for green plants versus Chl c for most eukaryotic phytoplankton. The enzymes that mediate Chl c biosynthesis have long remained elusive. In this work, we identified the CHLC dioxygenase (Phatr3_J43737) from the marine diatom Phaeodactylum tricornutum as the Chl c synthase. The chlc mutants lacked Chl c, instead accumulating its precursors, and exhibited growth defects. In vitro, recombinant CHLC protein converted these precursors into Chl c, thereby confirming its identity. Phylogenetic evidence demonstrates conserved use of CHLC across phyla but also the existence of distinct Chl c synthases in different algal groups. Our study addresses a long-outstanding question with implications for both contemporary and ancient marine photosynthesis.

Linking extreme seasonality and gene expression in Arctic marine protists

At high latitudes, strong seasonal differences in light availability affect marine organisms and regulate the timing of ecosystem processes. Marine protists are key players in Arctic aquatic ecosystems, yet little is known about their ecological roles over yearly cycles. This is especially true for the dark polar night period, which up until recently was assumed to be devoid of biological activity. A 12 million transcripts catalogue was built from 0.45 to 10 μm protist assemblages sampled over 13 months in a time series station in an Arctic fjord in Svalbard. Community gene expression was correlated with seasonality, with light as the main driving factor. Transcript diversity and evenness were higher during polar night compared to polar day. Light-dependent functions had higher relative expression during polar day, except phototransduction. 64% of the most expressed genes could not be functionally annotated, yet up to 78% were identified in Arctic samples from Tara Oceans, suggesting that Arctic marine assemblages are distinct from those from other oceans. Our study increases understanding of the links between extreme seasonality and biological processes in pico- and nanoplanktonic protists. Our results set the ground for future monitoring studies investigating the seasonal impact of climate change on the communities of microbial eukaryotes in the High Arctic.

Holistic view of the seascape dynamics and environment impact on macro-scale genetic connectivity of marine plankton populations

BACKGROUND

Plankton seascape genomics studies have revealed different trends from large-scale weak differentiation to microscale structures. Previous studies have underlined the influence of the environment and seascape on species differentiation and adaptation. However, these studies have generally focused on a few single species, sparse molecular markers, or local scales. Here, we investigated the genomic differentiation of plankton at the macro-scale in a holistic approach using Tara Oceans metagenomic data together with a reference-free computational method.

RESULTS

We reconstructed the FST-based genomic differentiation of 113 marine planktonic taxa occurring in the North and South Atlantic Oceans, Southern Ocean, and Mediterranean Sea. These taxa belong to various taxonomic clades spanning Metazoa, Chromista, Chlorophyta, Bacteria, and viruses. Globally, population genetic connectivity was significantly higher within oceanic basins and lower in bacteria and unicellular eukaryotes than in zooplankton. Using mixed linear models, we tested six abiotic factors influencing connectivity, including Lagrangian travel time, as proxies of oceanic current effects. We found that oceanic currents were the main population genetic connectivity drivers, together with temperature and salinity. Finally, we classified the 113 taxa into parameter-driven groups and showed that plankton taxa belonging to the same taxonomic rank such as phylum, class or order presented genomic differentiation driven by different environmental factors.

CONCLUSION

Our results validate the isolation-by-current hypothesis for a non-negligible proportion of taxa and highlight the role of other physicochemical parameters in large-scale plankton genetic connectivity. The reference-free approach used in this study offers a new systematic framework to analyse the population genomics of non-model and undocumented marine organisms from a large-scale and holistic point of view.

The pennate diatom Pseudo-nitzschia multistriata as a model for diatom life cycles, from the laboratory to the sea

Phytoplankton dynamics are regulated by external cues, such as light and nutrients, as well as by biotic interactions and endogenous controls linked to life cycle characteristics. The planktonic pennate diatom Pseudo-nitzschia multistriata, with a heterothallic mating system with two opposite mating types (MTs), represents a model for the study of diatom life cycles. P. multistriata is a toxic species, able to produce the neurotoxin domoic acid. First described in Japan in 1993, it was detected at the long-term monitoring station MareChiara (Gulf of Naples, Italy) in 1995. Since then, P. multistriata has been reported from several worldwide coastal sites. A large body of knowledge has been produced on its ecology, genetic diversity, and life cycle characteristics. The availability of these data, the ecological relevance of the Pseudo-nitzschia genus, and its controllable life cycle with a short generation time made it an ideal species to develop a genetic model system for diatoms. To enable functional studies, a 59 Mb genome sequence and several transcriptomic data were produced, and genetic transformation was optimized. These tools allowed the discovery of the first mating-type determining gene for diatoms. Gene expression studies and metabolomics analyses defined genes and molecules underpinning different phases of the process of sexual reproduction. This model system, developed to explore the genetics of diatom life cycles, offers the opportunity to parallel experimental observations in the laboratory using in situ meta-omics analyses along space and time, empowering knowledge on the biology and ecology of the genus.

Phylogeny and biogeography of the algal DMS-releasing enzyme in the global ocean

Phytoplankton produce the volatile dimethyl sulfide (DMS), an important infochemical mediating microbial interactions, which is also emitted to the atmosphere and affecting the global climate. Albeit the enzymatic source for DMS in eukaryotes was elucidated, namely a DMSP lyase (DL) called Alma1, we still lack basic knowledge regarding its taxonomic distribution. We defined unique sequence motifs which enable the identification of DL homologs (DLHs) in model systems and environmental populations. We used these motifs to predict DLHs in diverse algae by analyzing hundreds of genomic and transcriptomic sequences from model systems under stress conditions and from environmental samples. Our findings show that the DL enzyme is more taxonomically widespread than previously thought, as it is encoded by known algal taxa as haptophytes and dinoflagellates, but also by chlorophytes, pelagophytes and diatoms, which were conventionally considered to lack the DL enzyme. By exploring the Tara Oceans database, we showed that DLHs are widespread across the oceans and are predominantly expressed by dinoflagellates. Certain dinoflagellate DLHs were differentially expressed between the euphotic and mesopelagic zones, suggesting a functional specialization and an involvement in the metabolic plasticity of mixotrophic dinoflagellates. In specific regions as the Southern Ocean, DLH expression by haptophytes and diatoms was correlated with environmental drivers such as nutrient availability. The expanded repertoire of putative DL enzymes from diverse microbial origins and geographic niches suggests new potential players in the marine sulfur cycle and provides a foundation to study the cellular function of the DL enzyme in marine microbes.

Genome analysis of Parmales, the sister group of diatoms, reveals the evolutionary specialization of diatoms from phago-mixotrophs to photoautotrophs

The order Parmales (class Bolidophyceae) is a minor group of pico-sized eukaryotic marine phytoplankton that contains species with cells surrounded by silica plates. Previous studies revealed that Parmales is a member of ochrophytes and sister to diatoms (phylum Bacillariophyta), the most successful phytoplankton group in the modern ocean. Therefore, parmalean genomes can serve as a reference to elucidate both the evolutionary events that differentiated these two lineages and the genomic basis for the ecological success of diatoms vs. the more cryptic lifestyle of parmaleans. Here, we compare the genomes of eight parmaleans and five diatoms to explore their physiological and evolutionary differences. Parmaleans are predicted to be phago-mixotrophs. By contrast, diatoms have lost genes related to phagocytosis, indicating the ecological specialization from phago-mixotrophy to photoautotrophy in their early evolution. Furthermore, diatoms show significant enrichment in gene sets involved in nutrient uptake and metabolism, including iron and silica, in comparison with parmaleans. Overall, our results suggest a strong evolutionary link between the loss of phago-mixotrophy and specialization to a silicified photoautotrophic life stage early in diatom evolution after diverging from the Parmales lineage.

Gene expression dynamics of natural assemblages of heterotrophic flagellates during bacterivory

BACKGROUND

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

RESULTS

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

High Growth Rate of Diatoms Explained by Reduced Carbon Requirement and Low Energy Cost of Silica Deposition

The rapid growth of diatoms makes them one of the most pervasive and productive types of plankton in the world's ocean, but the physiological basis for their high growth rates remains poorly understood. Here, we evaluate the factors that elevate diatom growth rates, relative to other plankton, using a steady-state metabolic flux model that computes the photosynthetic C source from intracellular light attenuation and the carbon cost of growth from empirical cell C quotas, across a wide range of cell sizes. For both diatoms and other phytoplankton, growth rates decline with increased cell volume, consistent with observations, because the C cost of division increases with size faster than photosynthesis. However, the model predicts overall higher growth rates for diatoms due to reduced C requirements and the low energetic cost of Si deposition. The C savings from the silica frustule are supported by metatranscriptomic data from Tara Oceans, which show that the abundance of transcripts for cytoskeleton components in diatoms is lower than in other phytoplankton. Our results highlight the importance of understanding the origins of phylogenetic differences in cellular C quotas, and suggest that the evolution of silica frustules may play a critical role in the global dominance of marine diatoms. IMPORTANCE This study addresses a longstanding issue regarding diatoms, namely, their fast growth. Diatoms, which broadly are phytoplankton with silica frustules, are the world's most productive microorganisms and dominate in polar and upwelling regions. Their dominance is largely supported by their high growth rate, but the physiological reasoning behind that characteristic has been obscure. In this study, we combine a quantitative model and metatranscriptomic approaches and show that diatoms' low carbon requirements and low energy costs for silica frustule production are the key factors supporting their fast growth. Our study suggests that the effective use of energy-efficient silica as a cellular structure, instead of carbon, enables diatoms to be the most productive organisms in the global ocean.

In Situ Molecular Ecological Analyses Illuminate Distinct Factors Regulating Formation and Demise of a Harmful Dinoflagellate Bloom

The development and demise of a harmful algal bloom (HAB) are generally regulated by multiple processes; identifying specific critical drivers for a specific bloom is important yet challenging. Here, we conducted a whole-assemblage molecular ecological study on a dinoflagellate bloom to address the hypothesis that energy and nutrient acquisition, defense against grazing and microbial attacks, and sexual reproduction are critical to the rise and demise of the bloom. Microscopic and molecular analyses identified the bloom-causing species as Karenia longicanalis and showed that the ciliate Strombidinopsis sp. was dominant in a nonbloom plankton community, whereas the diatom Chaetoceros sp. dominated the after-bloom community, along with remarkable shifts in the community structure for both eukaryotes and prokaryotes. Metatranscriptomic analysis indicated that heightened energy and nutrient acquisition in K. longicanalis significantly contributed to bloom development. In contrast, active grazing by the ciliate Strombidinopsis sp. and attacks by algicidal bacteria (Rhodobacteracea, Cryomorphaceae, and Rhodobacteracea) and viruses prevented (at nonbloom stage) or collapsed the bloom (in after-bloom stage). Additionally, nutrition competition by the Chaetoceros diatoms plausibly contributed to bloom demise. The findings suggest the importance of energy and nutrients in promoting this K. longicanalis bloom and the failure of antimicrobial defense and competition of diatoms as the major bloom suppressor and terminator. This study provides novel insights into bloom-regulating mechanisms and the first transcriptomic data set of K. longicanalis, which will be a valuable resource and essential foundation for further elucidation of bloom regulators of this and related species of Kareniaceae in the future. IMPORTANCE HABs have increasingly occurred and impacted human health, aquatic ecosystems, and coastal economies. Despite great efforts, the factors that drive the development and termination of a bloom are poorly understood, largely due to inadequate in situ data about the physiology and metabolism of the causal species and the community. Using an integrative molecular ecological approach, we determined that heightened energy and nutrient acquisition promoted the bloom, while resource allocation in defense and failure to defend against grazing and microbial attacks likely prevented or terminated the bloom. Our findings reveal the differential roles of multiple abiotic and biotic environmental factors in driving the formation or demise of a toxic dinoflagellate bloom, suggesting the importance of a balanced biodiverse ecosystem in preventing a dinoflagellate bloom. The study also demonstrates the power of whole-assemblage metatranscriptomics coupled to DNA barcoding in illuminating plankton ecological processes and the underlying species and functional diversities.

Molecular discoveries in microbial DMSP synthesis

Dimethylsulfoniopropionate (DMSP) is one of the Earth's most abundant organosulfur compounds because many marine algae, bacteria, corals and some plants produce it to high mM intracellular concentrations. In these organisms, DMSP acts an anti-stress molecule with purported roles to protect against salinity, temperature, oxidative stress and hydrostatic pressure, amongst many other reported functions. However, DMSP is best known for being a major precursor of the climate-active gases and signalling molecules dimethylsulfide (DMS), methanethiol (MeSH) and, potentially, methane, through microbial DMSP catabolism. DMSP catabolism has been extensively studied and the microbes, pathways and enzymes involved have largely been elucidated through the application of molecular research over the last 17 years. In contrast, the molecular biology of DMSP synthesis is a much newer field, with the first DMSP synthesis enzymes only being identified in the last 5 years. In this review, we discuss how the elucidation of key DMSP synthesis enzymes has greatly expanded our knowledge of the diversity of DMSP-producing organisms, the pathways used, and what environmental factors regulate production, as well as to inform on the physiological roles of DMSP. Importantly, the identification of key DMSP synthesis enzymes in the major groups of DMSP producers has allowed scientists to study the distribution and predict the importance of different DMSP-producing organisms to global DMSP production in diverse marine and sediment environments. Finally, we highlight key challenges for future molecular research into DMSP synthesis that need addressing to better understand the cycling of this important marine organosulfur compound, and its magnitude in the environment.

On the Unknown Proteins of Eukaryotic Proteomes

To study unknown proteins on a large scale, a reference system has been set up for the three better studied eukaryotic kingdoms, built with 36 proteomes as taxonomically diverse as possible. Proteins from 362 other eukaryotic proteomes with no known homologue in this set were then analyzed, focusing noteworthy on singletons, that is, on such proteins with no known homologue in their own proteome. Consistently, for a given species, no more than 12% of the singletons thus found are known at the protein level, according to Uniprot. In addition, since they rely on the information found in the alignment of homologous sequences, predictions of AlphaFold2 for their tridimensional structure are poor. In the case of metazoan species, the number of singletons rarely exceeds 1000 for the species the closest to the reference system (divergence times below 75 Myr). Interestingly, in the cases of viridiplantae and fungi, larger amounts of singletons are found for such species, as if the timescale on which singletons are added to proteomes were different in metazoa and in other eukaryotic kingdoms. In order to confirm this phenomenon, further studies of proteomes closer to those of the reference system are, however, needed.

Marine phytoplankton downregulate core photosynthesis and carbon storage genes upon rapid mixed layer shallowing

Marine phytoplankton are a diverse group of photoautotrophic organisms and key mediators in the global carbon cycle. Phytoplankton physiology and biomass accumulation are closely tied to mixed layer depth, but the intracellular metabolic pathways activated in response to changes in mixed layer depth remain less explored. Here, metatranscriptomics was used to characterize the phytoplankton community response to a mixed layer shallowing (from 233 to 5 m) over the course of two days during the late spring in the Northwest Atlantic. Most phytoplankton genera downregulated core photosynthesis, carbon storage, and carbon fixation genes as the system transitioned from a deep to a shallow mixed layer and shifted towards catabolism of stored carbon supportive of rapid cell growth. In contrast, phytoplankton genera exhibited divergent transcriptional patterns for photosystem light harvesting complex genes during this transition. Active virus infection, taken as the ratio of virus to host transcripts, increased in the Bacillariophyta (diatom) phylum and decreased in the Chlorophyta (green algae) phylum upon mixed layer shallowing. A conceptual model is proposed to provide ecophysiological context for our findings, in which integrated light limitation and lower division rates during transient deep mixing are hypothesized to disrupt resource-driven, oscillating transcript levels related to photosynthesis, carbon fixation, and carbon storage. Our findings highlight shared and unique transcriptional response strategies within phytoplankton communities acclimating to the dynamic light environment associated with transient deep mixing and shallowing events during the annual North Atlantic bloom.

Virus infection of phytoplankton increases average molar mass and reduces hygroscopicity of aerosolized organic matter

Viral infection of phytoplankton is a pervasive mechanism of cell death and bloom termination, which leads to the production of dissolved and colloidal organic matter that can be aerosolized into the atmosphere. Earth-observing satellites can track the growth and death of phytoplankton blooms on weekly time scales but the impact of viral infection on the cloud forming potential of associated aerosols is largely unknown. Here, we determine the influence of viral-derived organic matter, purified viruses, and marine hydrogels on the cloud condensation nuclei activity of their aerosolized solutions, compared to organic exudates from healthy phytoplankton. Dissolved organic material derived from exponentially growing and infected cells of well-characterized eukaryotic phytoplankton host-virus systems, including viruses from diatoms, coccolithophores and chlorophytes, was concentrated, desalted, and nebulized to form aerosol particles composed of primarily of organic matter. Aerosols from infected phytoplankton cultures resulted in an increase in critical activation diameter and average molar mass in three out of five combinations evaluated, along with a decrease in organic kappa (hygroscopicity) compared to healthy cultures and seawater controls. The infected samples also displayed evidence of increased surface tension depression at realistic cloud water vapor supersaturations. Amending the samples with xanthan gum to simulate marine hydrogels increased variability in organic kappa and surface tension in aerosols with high organic to salt ratios. Our findings suggest that the pulses of increased dissolved organic matter associated with viral infection in surface waters may increase the molar mass of dissolved organic compounds relative to surface waters occupied by healthy phytoplankton or low phytoplankton biomass.

Quantifying Shared and Unique Gene Content across 17 Microbial Ecosystems

Measuring microbial diversity is traditionally based on microbe taxonomy. Here, in contrast, we aimed to quantify heterogeneity in microbial gene content across 14,183 metagenomic samples spanning 17 ecologies, including 6 human associated, 7 nonhuman host associated, and 4 in other nonhuman host environments. In total, we identified 117,629,181 nonredundant genes. The vast majority of genes (66%) occurred in only one sample (i.e., "singletons"). In contrast, we found 1,864 sequences present in every metagenome, but not necessarily every bacterial genome. Additionally, we report data sets of other ecology-associated genes (e.g., abundant in only gut ecosystems) and simultaneously demonstrated that prior microbiome gene catalogs are both incomplete and inaccurately cluster microbial genetic life (e.g., at gene sequence identities that are too restrictive). We provide our results and the sets of environmentally differentiating genes described above at http://www.microbial-genes.bio. IMPORTANCE The amount of shared genetic elements has not been quantified between the human microbiome and other host- and non-host-associated microbiomes. Here, we made a gene catalog of 17 different microbial ecosystems and compared them. We show that most species shared between environment and human gut microbiomes are pathogens and that prior gene catalogs described as "nearly complete" are far from it. Additionally, over two-thirds of all genes only appear in a single sample, and only 1,864 genes (0.001%) are found in all types of metagenomes. These results highlight the large diversity between metagenomes and reveal a new, rare class of genes, those found in every type of metagenome, but not every microbial genome.

Mirusviruses link herpesviruses to giant viruses

DNA viruses have a major influence on the ecology and evolution of cellular organisms1-4, but their overall diversity and evolutionary trajectories remain elusive5. Here we carried out a phylogeny-guided genome-resolved metagenomic survey of the sunlit oceans and discovered plankton-infecting relatives of herpesviruses that form a putative new phylum dubbed Mirusviricota. The virion morphogenesis module of this large monophyletic clade is typical of viruses from the realm Duplodnaviria6, with multiple components strongly indicating a common ancestry with animal-infecting Herpesvirales. Yet, a substantial fraction of mirusvirus genes, including hallmark transcription machinery genes missing in herpesviruses, are closely related homologues of giant eukaryotic DNA viruses from another viral realm, Varidnaviria. These remarkable chimaeric attributes connecting Mirusviricota to herpesviruses and giant eukaryotic viruses are supported by more than 100 environmental mirusvirus genomes, including a near-complete contiguous genome of 432 kilobases. Moreover, mirusviruses are among the most abundant and active eukaryotic viruses characterized in the sunlit oceans, encoding a diverse array of functions used during the infection of microbial eukaryotes from pole to pole. The prevalence, functional activity, diversification and atypical chimaeric attributes of mirusviruses point to a lasting role of Mirusviricota in the ecology of marine ecosystems and in the evolution of eukaryotic DNA viruses.

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.

Convergent evolution and horizontal gene transfer in Arctic Ocean microalgae

Microbial communities in the world ocean are affected strongly by oceanic circulation, creating characteristic marine biomes. The high connectivity of most of the ocean makes it difficult to disentangle selective retention of colonizing genotypes (with traits suited to biome specific conditions) from evolutionary selection, which would act on founder genotypes over time. The Arctic Ocean is exceptional with limited exchange with other oceans and ice covered since the last ice age. To test whether Arctic microalgal lineages evolved apart from algae in the global ocean, we sequenced four lineages of microalgae isolated from Arctic waters and sea ice. Here we show convergent evolution and highlight geographically limited HGT as an ecological adaptive force in the form of PFAM complements and horizontal acquisition of key adaptive genes. Notably, ice-binding proteins were acquired and horizontally transferred among Arctic strains. A comparison with Tara Oceans metagenomes and metatranscriptomes confirmed mostly Arctic distributions of these IBPs. The phylogeny of Arctic-specific genes indicated that these events were independent of bacterial-sourced HGTs in Antarctic Southern Ocean microalgae.

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.

Characterization of a Marine Diatom Chitin Synthase Using a Combination of Meta-Omics, Genomics, and Heterologous Expression Approaches

β-Chitin has important ecological and physiological roles and potential for widespread applications, but the characterization of chitin-related enzymes from β-chitin producers was rarely reported. Querying against the Tara Oceans Gene Atlas, 4,939 chitin-related unique sequences from 12 Pfam accessions were found in Bacillariophyta metatranscriptomes. Putative chitin synthase (CHS) sequences are decreasingly present in Crustacea (39%), Stramenopiles (16%) and Insecta (14%) from the Marine Atlas of Tara Oceans Unigenes version 1 Metatranscriptomes (MATOUv1+T) database. A CHS gene from the model diatom Thalassiosira pseudonana (Thaps3_J4413, designated TpCHS1) was identified. Homology analysis of TpCHS1 in Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP), PhycoCosm, and the PLAZA diatom omics data set showed that Mediophyceae and Thalassionemales species were potential new β-chitin producers besides Thalassiosirales. TpCHS1 was overexpressed in Saccharomyces cerevisiae and Phaeodactylum tricornutum. In transgenic P. tricornutum lines, TpCHS1-eGFP localizes to the Golgi apparatus and plasma membrane and predominantly accumulates in the cleavage furrow during cell division. Enhanced TpCHS1 expression could induce abnormal cell morphology and reduce growth rates in P. tricornutum, which might be ascribed to the inhibition of the G2/M phase. S. cerevisiae was proved to be a better system for expressing large amounts of active TpCHS1, which effectively incorporates UDP-N-acetylglucosamine in radiometric in vitro assays. Our study expands the knowledge on chitin synthase taxonomic distribution in marine eukaryotic microbes, and is the first to collectively characterize an active marine diatom CHS which may play an important role during cell division. IMPORTANCE As the most abundant biopolymer in the oceans, the significance of chitin and its biosynthesis is rarely demonstrated in diatoms, which are the main contributors to the primary productivity of the oceans, ascribed to their huge biomass and efficient photosynthesis. We retrieved genes involved in chitin-based metabolism against the Tara Oceans Gene Atlas to expand our knowledge about their diversity and distribution in the marine environment. Potential new producers of chitin were found from the analysis of various algal transcriptome and genome databases. Heterologous expression confirms that Thalassiosira pseudonana contains an active chitin synthase (CHS) which may play an important role in the cell division process of diatoms. This study provides new insight into CHS geographic and taxonomic distribution in marine eukaryotic microbes, as well as into a new CHS functioning in the biosynthesis of β-chitin in diatoms.

An ancient enzyme finds a new home: Prevalence and neofunctionalization of trypsin in marine phytoplankton

Trypsin is an ancient protease best known as a digestive enzyme in animals, and traditionally believed to be absent in plants and protists. However, our recent studies have revealed its wide presence and important roles in marine phytoplankton. Here, to gain a better understanding on the importance of trypsin in phytoplankton, we further surveyed the distribution, diversity, evolution and potential ecological roles of trypsin in global ocean phytoplankton. Our analysis indicated that trypsin is widely distributed both taxonomically and geographically in marine phytoplankton. Furthermore, by systematic comparative analyses we found that algal trypsin could be classified into two subfamilies (trypsin I and trypsin II) and exhibited highly duplicated and diversified during evolution. We also observed markedly different domain sequences and organizations between and within the subfamilies, suggesting potential neofunctionalization. Diatoms contain both subfamilies of trypsin, with higher numbers of genes and more environment-responsive expression of trypsin than other lineages. The duplication and subsequent neofunctionalization of the trypsin family may be important in diatoms for adapting to dynamical environmental conditions, contributing to diatoms' dominance in the coastal oceans. This work advances our knowledge on the distribution and neofunctionalization of this ancient enzyme and creates a new window of research on phytoplankton biology.

Structural Foundations of Potassium Selectivity in Channelrhodopsins

Potassium-selective channelrhodopsins (KCRs) are light-gated K⁺ channels recently found in the stramenopile protist Hyphochytrium catenoides. When expressed in neurons, KCRs enable high-precision optical inhibition of spiking (optogenetic silencing). KCRs are capable of discriminating K⁺ from Na⁺ without the conventional K⁺ selectivity filter found in classical K⁺ channels. The genome of H. catenoides also encodes a third paralog that is more permeable for Na⁺ than for K⁺. To identify structural motifs responsible for the unusual K⁺ selectivity of KCRs, we systematically analyzed a series of chimeras and mutants of this protein. We found that mutations of three critical residues in the paralog convert its Na⁺-selective channel into a K⁺-selective one. Our characterization of homologous proteins from other protists (Colponema vietnamica, Cafeteria burkhardae, and Chromera velia) and metagenomic samples confirmed the importance of these residues for K⁺ selectivity. We also show that Trp102 and Asp116, conserved in all three H. catenoides paralogs, are necessary, although not sufficient, for K⁺ selectivity. Our results provide the foundation for further engineering of KCRs for optogenetic needs. IMPORTANCE Recently discovered microbial light-gated ion channels (channelrhodopsins) with a higher permeability for K⁺ than for Na⁺ (potassium-selective channelrhodopsins [kalium channelrhodopsins, or KCRs]) demonstrate an alternative K⁺ selectivity mechanism, unrelated to well-characterized "selectivity filters" of voltage- and ligand-gated K⁺ channels. KCRs can be used for optogenetic inhibition of neuronal firing and potentially for the development of gene therapies to treat neurological and cardiovascular disorders. In this study, we identified structural motifs that determine the K⁺ selectivity of KCRs that provide the foundation for their further improvement as optogenetic tools.

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.

Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

High-quality genome of the basidiomycete yeast Dioszegia hungarica PDD-24b-2 isolated from cloud water

The genome of the basidiomycete yeast Dioszegia hungarica strain PDD-24b-2 isolated from cloud water at the summit of puy de Dôme (France) was sequenced using a hybrid PacBio and Illumina sequencing strategy. The obtained assembled genome of 20.98 Mb and a GC content of 57% is structured in 16 large-scale contigs ranging from 90 kb to 5.56 Mb, and another 27.2 kb contig representing the complete circular mitochondrial genome. In total, 8,234 proteins were predicted from the genome sequence. The mitochondrial genome shows 16.2% cgu codon usage for arginine but has no canonical cognate tRNA to translate this codon. Detected transposable element (TE)-related sequences account for about 0.63% of the assembled genome. A dataset of 2,068 hand-picked public environmental metagenomes, representing over 20 Tbp of raw reads, was probed for D. hungarica related ITS sequences, and revealed worldwide distribution of this species, particularly in aerial habitats. Growth experiments suggested a psychrophilic phenotype and the ability to disperse by producing ballistospores. The high-quality assembled genome obtained for this D. hungarica strain will help investigate the behavior and ecological functions of this species in the environment.

Discovering marine biodiversity in the 21st century

We review the current knowledge of the biodiversity of the ocean as well as the levels of decline and threat for species and habitats. The lack of understanding of the distribution of life in the ocean is identified as a significant barrier to restoring its biodiversity and health. We explore why the science of taxonomy has failed to deliver knowledge of what species are present in the ocean, how they are distributed and how they are responding to global and regional to local anthropogenic pressures. This failure prevents nations from meeting their international commitments to conserve marine biodiversity with the results that investment in taxonomy has declined in many countries. We explore a range of new technologies and approaches for discovery of marine species and their detection and monitoring. These include: imaging methods, molecular approaches, active and passive acoustics, the use of interconnected databases and citizen science. Whilst no one method is suitable for discovering or detecting all groups of organisms many are complementary and have been combined to give a more complete picture of biodiversity in marine ecosystems. We conclude that integrated approaches represent the best way forwards for accelerating species discovery, description and biodiversity assessment. Examples of integrated taxonomic approaches are identified from terrestrial ecosystems. Such integrated taxonomic approaches require the adoption of cybertaxonomy approaches and will be boosted by new autonomous sampling platforms and development of machine-speed exchange of digital information between databases.

Expansion of the global RNA virome reveals diverse clades of bacteriophages

High-throughput RNA sequencing offers broad opportunities to explore the Earth RNA virome. Mining 5,150 diverse metatranscriptomes uncovered >2.5 million RNA virus contigs. Analysis of >330,000 RNA-dependent RNA polymerases (RdRPs) shows that this expansion corresponds to a 5-fold increase of the known RNA virus diversity. Gene content analysis revealed multiple protein domains previously not found in RNA viruses and implicated in virus-host interactions. Extended RdRP phylogeny supports the monophyly of the five established phyla and reveals two putative additional bacteriophage phyla and numerous putative additional classes and orders. The dramatically expanded phylum Lenarviricota, consisting of bacterial and related eukaryotic viruses, now accounts for a third of the RNA virome. Identification of CRISPR spacer matches and bacteriolytic proteins suggests that subsets of picobirnaviruses and partitiviruses, previously associated with eukaryotes, infect prokaryotic hosts.

A metabolic, phylogenomic and environmental atlas of diatom plastid transporters from the model species Phaeodactylum

Diatoms are an important group of algae, contributing nearly 40% of total marine photosynthetic activity. However, the specific molecular agents and transporters underpinning the metabolic efficiency of the diatom plastid remain to be revealed. We performed in silico analyses of 70 predicted plastid transporters identified by genome-wide searches of Phaeodactylum tricornutum. We considered similarity with Arabidopsis thaliana plastid transporters, transcriptional co-regulation with genes encoding core plastid metabolic pathways and with genes encoded in the mitochondrial genomes, inferred evolutionary histories using single-gene phylogeny, and environmental expression trends using Tara Oceans meta-transcriptomics and meta-genomes data. Our data reveal diatoms conserve some of the ion, nucleotide and sugar plastid transporters associated with plants, such as non-specific triose phosphate transporters implicated in the transport of phosphorylated sugars, NTP/NDP and cation exchange transporters. However, our data also highlight the presence of diatom-specific transporter functions, such as carbon and amino acid transporters implicated in intricate plastid-mitochondria crosstalk events. These confirm previous observations that substrate non-specific triose phosphate transporters (TPT) may exist as principal transporters of phosphorylated sugars into and out of the diatom plastid, alongside suggesting probable agents of NTP exchange. Carbon and amino acid transport may be related to intricate metabolic plastid-mitochondria crosstalk. We additionally provide evidence from environmental meta-transcriptomic/meta- genomic data that plastid transporters may underpin diatom sensitivity to ocean warming, and identify a diatom plastid transporter (J43171) whose expression may be positively correlated with temperature.

Seasonal microbial dynamics in the ocean inferred from assembled and unassembled data: a view on the unknown biosphere

In environmental metagenomic experiments, a very high proportion of the microbial sequencing data (> 70%) remains largely unexploited because rare and closely related genomes are missed in short-read assemblies. The identity and the potential metabolisms of a large fraction of natural microbial communities thus remain inaccessible to researchers. The purpose of this study was to explore the genomic content of unassembled metagenomic data and test their level of novelty. We used data from a three-year microbial metagenomic time series of the NW Mediterranean Sea, and conducted reference-free and database-guided analysis. The results revealed a significant genomic difference between the assembled and unassembled reads. The unassembled reads had a lower mean identity against public databases, and fewer metabolic pathways could be reconstructed. In addition, the unassembled fraction presented a clear temporal pattern, unlike the assembled ones, and a specific community composition that was similar to the rare communities defined by metabarcoding using the 16S rRNA gene. The rare gene pool was characterised by keystone bacterial taxa, and the presence of viruses, suggesting that viral lysis could maintain some taxa in a state of rarity. Our study demonstrates that unassembled metagenomic data can provide important information on the structure and functioning of microbial communities.

Genomic adaptation of the picoeukaryote Pelagomonas calceolata to iron-poor oceans revealed by a chromosome-scale genome sequence

The smallest phytoplankton species are key actors in oceans biogeochemical cycling and their abundance and distribution are affected with global environmental changes. Among them, algae of the Pelagophyceae class encompass coastal species causative of harmful algal blooms while others are cosmopolitan and abundant. The lack of genomic reference in this lineage is a main limitation to study its ecological importance. Here, we analysed Pelagomonas calceolata relative abundance, ecological niche and potential for the adaptation in all oceans using a complete chromosome-scale assembled genome sequence. Our results show that P. calceolata is one of the most abundant eukaryotic species in the oceans with a relative abundance favoured by high temperature, low-light and iron-poor conditions. Climate change projections based on its relative abundance suggest an extension of the P. calceolata habitat toward the poles at the end of this century. Finally, we observed a specific gene repertoire and expression level variations potentially explaining its ecological success in low-iron and low-nitrate environments. Collectively, these findings reveal the ecological importance of P. calceolata and lay the foundation for a global scale analysis of the adaptation and acclimation strategies of this small phytoplankton in a changing environment.

Genomic inference reveals potential climate change-driven range expansion of the phytoplankton species Pelagomonas calceolata.