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

Revealing the intricate temporal dynamics and adaptive responses of prokaryotic and eukaryotic microbes in the coastal South China Sea

This comprehensive two-year investigation in the coastal South China Sea has advanced our understanding of marine microbes at both community and genomic levels. By combining metagenomics and metatranscriptomics, we have revealed the intricate temporal dynamics and remarkable adaptability of microbial communities and phytoplankton metagenome-assembled genomes (MAGs) in response to environmental fluctuations. We observed distinct seasonal shifts in microbial community composition and function: cyanobacteria were predominant during warmer months, whereas photosynthetic protists were more abundant during colder seasons. Notably, metabolic marker KOs of photosynthesis were consistently active throughout the year, underscoring the persistent role of these processes irrespective of seasonal changes. Our analysis reveals that environmental parameters such as temperature, salinity, and nitrate concentrations profoundly influence microbial community composition, while temperature and silicate have emerged as crucial factors shaping their functional traits. Through the recovery and analysis of 37 phytoplankton MAGs, encompassing nine prokaryotic cyanobacteria and 28 eukaryotic protists from diverse phyla, we have gained insights into their genetic diversity and metabolic capabilities. Distinct profiles of photosynthesis-related pathways including carbon fixation, carotenoid biosynthesis, photosynthesis-antenna proteins, and photosynthesis among the MAGs indicated their genetic adaptations to changing environmental conditions. This study not only enhances our understanding of microbial dynamics in coastal marine ecosystems but also sheds light on the ecological roles and adaptive responses of different microbial groups to environmental changes.

Presence of vitamin B12 metabolism in the last common ancestor of land plants

Vitamin B12, also known as cobalamin, is an essential organic cofactor for methionine synthase (METH), and is only synthesized by a subset of bacteria. Plants and fungi have an alternative methionine synthase (METE) that does not need B12 and are typically considered not to utilize it. Some algae facultatively utilize B12 because they encode both METE and METH, while other algae are dependent on B12 as they encode METH only. We performed phylogenomic analyses of METE, METH and 11 further proteins involved in B12 metabolism across more than 1600 plant and algal genomes and transcriptomes (e.g. from OneKp), demonstrating the presence of B12-associated metabolism deep into the streptophytes. METH and five further accessory proteins (MTRR, CblB, CblC, CblD and CblJ) were detected in the hornworts (Anthocerotophyta), and two (CblB and CblJ) were identified in liverworts (Marchantiophyta) in the bryophytes, suggesting a retention of B12-metabolism in the last common land plant ancestor. Our data further show more limited distributions for other B12-related proteins (MCM and RNR-II) and B12 dependency in several algal orders. Finally, considering the collection sites of algae that have lost B12 metabolism, we propose freshwater-to-land transitions and symbiotic associations to have been constraining factors for B12 availability in early plant evolution. This article is part of the theme issue 'The evolution of plant metabolism'.

How marine are Marine Stramenopiles (MAST)? A cross-system evaluation

Marine Stramenopiles (MAST) were first described two decades ago through ribosomal RNA gene (rRNA gene) sequences from marine surveys of microbial eukaryotes. MAST comprise several independent lineages at the base of the Stramenopiles. Despite their prevalence in the ocean, the majority of MAST diversity remains uncultured. Previous studies, mainly in marine environments, have explored MAST's cell morphology, distribution, trophic strategies, and genomics using culturing-independent methods. In comparison, less is known about their presence outside marine habitats. Here, we analyse the extensive EukBank dataset to assess the extent to which MAST can be considered marine protists. Additionally, by incorporating newly available rRNA gene sequences, we update Stramenopiles phylogeny, identifying three novel MAST lineages. Our results indicate that MAST are primarily marine with notable exceptions within MAST-2 and MAST-12, where certain subclades are prevalent in freshwater and soil habitats. In the marine water column, only a few MAST species, particularly within clades -1, -3, -4, and -7, dominate and exhibit clear latitudinal distribution patterns. Overall, the massive sequencing dataset analysed in our study confirms and partially expands the previously described diversity of MASTs groups and underscores the predominantly marine nature of most of these uncultured lineages.

Eukaryotic viruses encode the ribosomal protein eL40

Viruses in the phylum Nucleocytoviricota are large, complex and have an exceptionally diverse metabolic repertoire. Some encode hundreds of products involved in the translation of mRNA into protein, but none was known to encode any of the proteins in ribosomes, the central engines of translation. With the discovery of the eL40 gene in FloV-SA2, we report the first example of a eukaryotic virus encoding a ribosomal protein and show that this gene is also present and expressed in other uncultivated marine giant viruses. FloV-SA2 also encodes a "group II" viral rhodopsin, a viral light-activated protein of unknown function previously only reported in metagenomes. FloV-SA2 is thus a valuable model system for investigating new mechanisms by which viruses manipulate eukaryotic cell metabolism.

Diatom pyrenoids are encased in a protein shell that enables efficient CO2 fixation

Pyrenoids are subcompartments of algal chloroplasts that increase the efficiency of Rubisco-driven CO2 fixation. Diatoms fix up to 20% of global CO2, but their pyrenoids remain poorly characterized. Here, we used in vivo photo-crosslinking to identify pyrenoid shell (PyShell) proteins, which we localized to the pyrenoid periphery of model pennate and centric diatoms, Phaeodactylum tricornutum and Thalassiosira pseudonana. In situ cryo-electron tomography revealed that pyrenoids of both diatom species are encased in a lattice-like protein sheath. Single-particle cryo-EM yielded a 2.4-Å-resolution structure of an in vitro TpPyShell1 lattice, which showed how protein subunits interlock. T. pseudonana TpPyShell1/2 knockout mutants had no PyShell sheath, altered pyrenoid morphology, and a high-CO2 requiring phenotype, with reduced photosynthetic efficiency and impaired growth under standard atmospheric conditions. The structure and function of the diatom PyShell provide a molecular view of how CO2 is assimilated in the ocean, a critical ecosystem undergoing rapid change.

Biodiversity of microorganisms in the Baltic Sea: the power of novel methods in the identification of marine microbes

Until recently, the data on the diversity of the entire microbial community from the Baltic Sea were relatively rare and very scarce. However, modern molecular methods have provided new insights into this field with interesting results. They can be summarized as follows. (i) Although low salinity causes a reduction in the biodiversity of multicellular species relative to the populations of the North-East Atlantic, no such reduction occurs in bacterial diversity. (ii) Among cyanobacteria, the picocyanobacterial group dominates when considering gene abundance, while filamentous cyanobacteria dominate in means of biomass. (iii) The diversity of diatoms and dinoflagellates is significantly larger than described a few decades ago; however, molecular studies on these groups are still scarce. (iv) Knowledge gaps in other protistan communities are evident. (v) Salinity is the main limiting parameter of pelagic fungal community composition, while the benthic fungal diversity is shaped by water depth, salinity, and sediment C and N availability. (vi) Bacteriophages are the predominant group of viruses, while among viruses infecting eukaryotic hosts, Phycodnaviridae are the most abundant; the Baltic Sea virome is contaminated with viruses originating from urban and/or industrial habitats. These features make the Baltic Sea microbiome specific and unique among other marine environments.

From nets to networks: tools for deciphering phytoplankton metabolic interactions within communities and their global significance

Our oceans are populated with a wide diversity of planktonic organisms that form complex dynamic communities at the base of marine trophic networks. Within such communities are phytoplankton, unicellular photosynthetic taxa that provide an estimated half of global primary production and support biogeochemical cycles, along with other essential ecosystem services. One of the major challenges for microbial ecologists has been to try to make sense of this complexity. While phytoplankton distributions can be well explained by abiotic factors such as temperature and nutrient availability, there is increasing evidence that their ecological roles are tightly linked to their metabolic interactions with other plankton members through complex mechanisms (e.g. competition and symbiosis). Therefore, unravelling phytoplankton metabolic interactions is the key for inferring their dependency on, or antagonism with, other taxa and better integrating them into the context of carbon and nutrient fluxes in marine trophic networks. In this review, we attempt to summarize the current knowledge brought by ecophysiology, organismal imaging, in silico predictions and co-occurrence networks using 'omics data, highlighting successful combinations of approaches that may be helpful for future investigations of phytoplankton metabolic interactions within their complex communities.This article is part of the theme issue 'Connected interactions: enriching food web research by spatial and social interactions'.

Complementary environmental analysis and functional characterization of lower glycolysis-gluconeogenesis in the diatom plastid

Organic carbon fixed in chloroplasts through the Calvin-Benson-Bassham Cycle can be diverted toward different metabolic fates, including cytoplasmic and mitochondrial respiration, gluconeogenesis, and synthesis of diverse plastid metabolites via the pyruvate hub. In plants, pyruvate is principally produced via cytoplasmic glycolysis, although a plastid-targeted lower glycolytic pathway is known to exist in non-photosynthetic tissue. Here, we characterized a lower plastid glycolysis-gluconeogenesis pathway enabling the direct interconversion of glyceraldehyde-3-phosphate and phospho-enol-pyruvate in diatoms, ecologically important marine algae distantly related to plants. We show that two reversible enzymes required to complete diatom plastid glycolysis-gluconeogenesis, Enolase and bis-phosphoglycerate mutase (PGAM), originated through duplications of mitochondria-targeted respiratory isoforms. Through CRISPR-Cas9 mutagenesis, integrative 'omic analyses, and measured kinetics of expressed enzymes in the diatom Phaeodactylum tricornutum, we present evidence that this pathway diverts plastid glyceraldehyde-3-phosphate into the pyruvate hub, and may also function in the gluconeogenic direction. Considering experimental data, we show that this pathway has different roles dependent in particular on day length and environmental temperature, and show that the cpEnolase and cpPGAM genes are expressed at elevated levels in high-latitude oceans where diatoms are abundant. Our data provide evolutionary, meta-genomic, and functional insights into a poorly understood yet evolutionarily recurrent plastid metabolic pathway.

Recovery of 240 metagenome-assembled genomes from coastal mariculture environments in South Korea

The mariculture industry has seen a rapid expansion in recent years due to the increasing global demand for seafood. However, the industry faces challenges from climate change and increased pathogen pressure. Additionally, the chemicals used to enhance mariculture productivity are changing ocean ecosystems. This study analyzed 36 surface-water metagenomes from South Korean mussel, oyster, scallop, and shrimp farms to expand our understanding of aquaculture microbial genetic resources and the potential impacts of these anthropogenic inputs. We recovered 240 non-redundant species-level metagenome-assembled genomes (MAGs), comprising 224 bacteria, 13 archaea, and three eukaryotes. Most MAGs were assigned to Proteobacteria, Bacteroidota, and Actinobacteriota, with 40.7% remaining unclassified at the species level. Among the three eukaryotic MAGs, one was identified as a novel lineage of green algae, highlighting the uncharacterized genetic diversity in mariculture environments. Additionally, 22 prokaryotic MAGs harbored 26 antibiotic and metal resistance genes, with MAGs carrying beta-lactamases being particularly prevalent in most farms. The obtained microbiome data from mariculture environments can be utilized in future studies to foster healthy, sustainable mariculture practices.

The genomic potential of photosynthesis in piconanoplankton is functionally redundant but taxonomically structured at a global scale

Carbon fixation is a key metabolic function shaping marine life, but the underlying taxonomic and functional diversity involved is only partially understood. Using metagenomic resources targeted at marine piconanoplankton, we provide a reproducible machine learning framework to derive the potential biogeography of genomic functions through the multi-output regression of gene read counts on environmental climatologies. Leveraging the Marine Atlas of Tara Oceans Unigenes, we investigate the genomic potential of primary production in the global ocean. The latter is performed by ribulose-1,5-bisphosphate carboxylase/oxygenase (RUBISCO) and is often associated with carbon concentration mechanisms in piconanoplankton, major marine unicellular photosynthetic organisms. We show that the genomic potential supporting C4 enzymes and RUBISCO exhibits strong functional redundancy and important affinity toward tropical oligotrophic waters. This redundancy is taxonomically structured by the dominance of Mamiellophyceae and Prymnesiophyceae in mid and high latitudes. These findings enhance our understanding of the relationship between functional and taxonomic diversity of microorganisms and environmental drivers of key biogeochemical cycles.

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.

Unveiling the microbial realm with VEBA 2.0: a modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic and viral multi-omics from either short- or long-read sequencing

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible software suite that bridges the gap between genomics and biotechnological solutions.

Gene Transfer Among Viruses Substantially Contributes to Gene Gain of Giant Viruses

The phylum Nucleocytoviricota comprises a diverse group of double-stranded DNA viruses that display a wide range of gene repertoires. Although these gene repertoires determine the characteristics of individual viruses, the evolutionary processes that have shaped the gene repertoires of extant viruses since their common ancestor are poorly characterized. In this study, we aimed to address this gap in knowledge by using amalgamated likelihood estimation, a probabilistic tree reconciliation method that infers evolutionary scenarios by distinguishing origination, gene duplications, virus-to-virus horizontal gene transfer (vHGT), and gene losses. We analyzed over 4,700 gene families from 195 genomes spanning all known viral orders. The evolutionary reconstruction suggests a history of extensive gene gains and losses during the evolution of these viruses, notably with vHGT contributing to gene gains at a comparable level to duplications and originations. The vHGT frequently occurred between phylogenetically closely related viruses, as well as between distantly related viruses with an overlapping host range. We observed a pattern of massive gene duplications that followed vHGTs for gene families that was potentially related to host range control and virus-host arms race. These results suggest that vHGT represents a previously overlooked, yet important, evolutionary force that integrates the evolutionary paths of multiple viruses and affects shaping of Nucleocytoviricota virus gene repertoires.

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.

Diverse patterns of correspondence between protist metabarcodes and protist metagenome-assembled genomes

Two common approaches to study the composition of environmental protist communities are metabarcoding and metagenomics. Raw metabarcoding data are usually processed into Operational Taxonomic Units (OTUs) or amplicon sequence variants (ASVs) through clustering or denoising approaches, respectively. Analogous approaches are used to assemble metagenomic reads into metagenome-assembled genomes (MAGs). Understanding the correspondence between the data produced by these two approaches can help to integrate information between the datasets and to explain how metabarcoding OTUs and MAGs are related with the underlying biological entities they are hypothesised to represent. MAGs do not contain the commonly used barcoding loci, therefore sequence homology approaches cannot be used to match OTUs and MAGs. We made an attempt to match V9 metabarcoding OTUs from the 18S rRNA gene (V9 OTUs) and MAGs from the Tara Oceans expedition based on the correspondence of their relative abundances across the same set of samples. We evaluated several metrics for detecting correspondence between features in these two datasets and developed controls to filter artefacts of data structure and processing. After selecting the best-performing metrics, ranking the V9 OTU/MAG matches by their proportionality/correlation coefficients and applying a set of selection criteria, we identified candidate matches between V9 OTUs and MAGs. In some cases, V9 OTUs and MAGs could be matched with a one-to-one correspondence, implying that they likely represent the same underlying biological entity. More generally, matches we observed could be classified into 4 scenarios: one V9 OTU matches many MAGs; many V9 OTUs match many MAGs; many V9 OTUs match one MAG; one V9 OTU matches one MAG. Notably, we found some instances in which different OTU-MAG matches from the same taxonomic group were not classified in the same scenario, with all four scenarios possible even within the same taxonomic group, illustrating that factors beyond taxonomic lineage influence the relationship between OTUs and MAGs. Overall, each scenario produces a different interpretation of V9 OTUs, MAGs and how they compare in terms of the genomic and ecological diversity they represent.

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.

Exploring the roles of ribosomal peptides in prokaryote-phage interactions through deep learning-enabled metagenome mining

BACKGROUND

Microbial secondary metabolites play a crucial role in the intricate interactions within the natural environment. Among these metabolites, ribosomally synthesized and post-translationally modified peptides (RiPPs) are becoming a promising source of therapeutic agents due to their structural diversity and functional versatility. However, their biosynthetic capacity and ecological functions remain largely underexplored.

RESULTS

Here, we aim to explore the biosynthetic profile of RiPPs and their potential roles in the interactions between microbes and viruses in the ocean, which encompasses a vast diversity of unique biomes that are rich in interactions and remains chemically underexplored. We first developed TrRiPP to identify RiPPs from ocean metagenomes, a deep learning method that detects RiPP precursors in a hallmark gene-independent manner to overcome the limitations of classic methods in processing highly fragmented metagenomic data. Applying this method to metagenomes from the global ocean microbiome, we uncover a diverse array of previously uncharacterized putative RiPP families with great novelty and diversity. Through correlation analysis based on metatranscriptomic data, we observed a high prevalence of antiphage defense-related and phage-related protein families that were co-expressed with RiPP families. Based on this putative association between RiPPs and phage infection, we constructed an Ocean Virus Database (OVD) and established a RiPP-involving host-phage interaction network through host prediction and co-expression analysis, revealing complex connectivities linking RiPP-encoding prokaryotes, RiPP families, viral protein families, and phages. These findings highlight the potential of RiPP families involved in prokaryote-phage interactions and coevolution, providing insights into their ecological functions in the ocean microbiome.

CONCLUSIONS

This study provides a systematic investigation of the biosynthetic potential of RiPPs from the ocean microbiome at a global scale, shedding light on the essential insights into the ecological functions of RiPPs in prokaryote-phage interactions through the integration of deep learning approaches, metatranscriptomic data, and host-phage connectivity. This study serves as a valuable example of exploring the ecological functions of bacterial secondary metabolites, particularly their associations with unexplored microbial interactions. Video Abstract.

Broad-spectrum hydrocarbon-degrading microbes in the global ocean metagenomes

Understanding the diversity and functions of hydrocarbon-degrading microorganisms in marine environments is crucial for both advancing knowledge of biogeochemical processes and improving bioremediation methods. In this study, we leveraged nearly 20,000 metagenome-assembled genomes (MAGs), recovered from a wide array of marine samples across the global oceans, to map the diversity of aerobic hydrocarbon-degrading microorganisms. A broad bacterial diversity was uncovered, with a notable preference for degrading aliphatic hydrocarbons over aromatic ones, primarily within Proteobacteria and Actinobacteriota. Three types of broad-spectrum hydrocarbon-degrading bacteria were identified for their ability to degrade various hydrocarbons and possession of multiple copies of hydrocarbon biodegradation genes. These bacteria demonstrate extensive metabolic versatility, aiding their survival and adaptability in diverse environmental conditions. Evidence of gene duplication and horizontal gene transfer in these microbes suggested a potential enhancement in the diversity of hydrocarbon-degrading bacteria. Positive correlations were observed between the abundances of hydrocarbon-degrading genes and environmental parameters such as temperature (-5 to 35 °C) and salinity (20 to 42 PSU). Overall, our findings offer valuable insights into marine hydrocarbon-degrading microorganisms and suggest considerations for selecting microbial strains for oil pollution remediation.

Genomic characterisation and ecological distribution of Mantoniella tinhauana: a novel Mamiellophycean green alga from the Western Pacific

Mamiellophyceae are dominant marine algae in much of the ocean, the most prevalent genera belonging to the order Mamiellales: Micromonas, Ostreococcus and Bathycoccus, whose genetics and global distributions have been extensively studied. Conversely, the genus Mantoniella, despite its potential ecological importance, remains relatively under-characterised. In this study, we isolated and characterised a novel species of Mamiellophyceae, Mantoniella tinhauana, from subtropical coastal waters in the South China Sea. Morphologically, it resembles other Mantoniella species; however, a comparative analysis of the 18S and ITS2 marker genes revealed its genetic distinctiveness. Furthermore, we sequenced and assembled the first genome of Mantoniella tinhauana, uncovering significant differences from previously studied Mamiellophyceae species. Notably, the genome lacked any detectable outlier chromosomes and exhibited numerous unique orthogroups. We explored gene groups associated with meiosis, scale and flagella formation, shedding light on species divergence, yet further investigation is warranted. To elucidate the biogeography of Mantoniella tinhauana, we conducted a comprehensive analysis using global metagenomic read mapping to the newly sequenced genome. Our findings indicate this species exhibits a cosmopolitan distribution with a low-level prevalence worldwide. Understanding the intricate dynamics between Mamiellophyceae and the environment is crucial for comprehending their impact on the ocean ecosystem and accurately predicting their response to forthcoming environmental changes.

Petabase-Scale Homology Search for Structure Prediction

The recent CASP15 competition highlighted the critical role of multiple sequence alignments (MSAs) in protein structure prediction, as demonstrated by the success of the top AlphaFold2-based prediction methods. To push the boundaries of MSA utilization, we conducted a petabase-scale search of the Sequence Read Archive (SRA), resulting in gigabytes of aligned homologs for CASP15 targets. These were merged with default MSAs produced by ColabFold-search and provided to ColabFold-predict. By using SRA data, we achieved highly accurate predictions (GDT_TS > 70) for 66% of the non-easy targets, whereas using ColabFold-search default MSAs scored highly in only 52%. Next, we tested the effect of deep homology search and ColabFold's advanced features, such as more recycles, on prediction accuracy. While SRA homologs were most significant for improving ColabFold's CASP15 ranking from 11th to 3rd place, other strategies contributed too. We analyze these in the context of existing strategies to improve prediction.

Unveiling the Microbial Realm with VEBA 2.0: A modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic, and viral multi-omics from either short- or long-read sequencing

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible platform that bridges the gap between genomics and biotechnological solutions.

FicD genes in invertebrates: A tale of transposons, pathogenic and integrated viruses

Many gene families are shared across the tree of life between distantly related species because of horizontal gene transfers (HGTs). However, the frequency of HGTs varies strongly between gene families and biotic realms suggesting differential selection pressures and functional bias. One gene family with a wide distribution are FIC-domain containing enzymes (FicDs). FicDs catalyze AMPylation, a post-translational protein modification consisting in the addition of adenosine monophosphate to accessible residues of target proteins. Beside the well-known conservation of FicDs in deuterostomes, we report the presence of a conserved FicD gene ortholog in a large number of protostomes and microbial eukaryotes. We also reported additional FicD gene copies in the genomes of some rotifers, parasitic worms and bivalves. A few dsDNA viruses of these invertebrates, including White spot syndrome virus, Cherax quadricarinatus iridovirus, Ostreid herpesvirus-1 and the beetle nudivirus, carry copies of FicDs, with phylogenetic analysis suggesting a common origin of these FicD copies and the duplicated FicDs of their invertebrate hosts. HGTs and gene duplications possibly mediated by endogenous viruses or genetic mobile elements seem to have contributed to the transfer of AMPylation ability from bacteria and eukaryotes to pathogenic viruses, where this pathway could have been hijacked to promote viral infection.

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

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

Integrated overview of stramenopile ecology, taxonomy, and heterotrophic origin

Stramenopiles represent a significant proportion of aquatic and terrestrial biota. Most biologists can name a few, but these are limited to the phototrophic (e.g. diatoms and kelp) or parasitic species (e.g. oomycetes, Blastocystis), with free-living heterotrophs largely overlooked. Though our attention is slowly turning towards heterotrophs, we have only a limited understanding of their biology due to a lack of cultured models. Recent metagenomic and single-cell investigations have revealed the species richness and ecological importance of stramenopiles-especially heterotrophs. However, our lack of knowledge of the cell biology and behaviour of these organisms leads to our inability to match species to their particular ecological functions. Because photosynthetic stramenopiles are studied independently of their heterotrophic relatives, they are often treated separately in the literature. Here, we present stramenopiles as a unified group with shared synapomorphies and evolutionary history. We introduce the main lineages, describe their important biological and ecological traits, and provide a concise update on the origin of the ochrophyte plastid. We highlight the crucial role of heterotrophs and mixotrophs in our understanding of stramenopiles with the goal of inspiring future investigations in taxonomy and life history. To understand each of the many diversifications within stramenopiles-towards autotrophy, osmotrophy, or parasitism-we must understand the ancestral heterotrophic flagellate from which they each evolved. We hope the following will serve as a primer for new stramenopile researchers or as an integrative refresher to those already in the field.

NMPFamsDB: a database of novel protein families from microbial metagenomes and metatranscriptomes

The Novel Metagenome Protein Families Database (NMPFamsDB) is a database of metagenome- and metatranscriptome-derived protein families, whose members have no hits to proteins of reference genomes or Pfam domains. Each protein family is accompanied by multiple sequence alignments, Hidden Markov Models, taxonomic information, ecosystem and geolocation metadata, sequence and structure predictions, as well as 3D structure models predicted with AlphaFold2. In its current version, NMPFamsDB hosts over 100 000 protein families, each with at least 100 members. The reported protein families significantly expand (more than double) the number of known protein sequence clusters from reference genomes and reveal new insights into their habitat distribution, origins, functions and taxonomy. We expect NMPFamsDB to be a valuable resource for microbial proteome-wide analyses and for further discovery and characterization of novel functions. NMPFamsDB is publicly available in http://www.nmpfamsdb.org/ or https://bib.fleming.gr/NMPFamsDB.

DiGAlign: Versatile and Interactive Visualization of Sequence Alignment for Comparative Genomics

With the explosion of available genomic information, comparative genomics has become a central approach to understanding microbial ecology and evolution. We developed DiGAlign (https://www.genome.jp/digalign/), a web server that provides versatile functionality for comparative genomics with an intuitive interface. It allows the user to perform the highly customizable visualization of a synteny map by simply uploading nucleotide sequences of interest, ranging from a specific region to the whole genome landscape of microorganisms and viruses. DiGAlign will serve a wide range of biological researchers, particularly experimental biologists, with multifaceted features that allow the rapid characterization of genomic sequences of interest and the generation of a publication-ready figure.

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

IMPORTANCE

Single-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.

Multi-omics for studying and understanding polar life

Polar ecosystems are experiencing amongst the most rapid rates of regional warming on Earth. Here, we discuss 'omics' approaches to investigate polar biodiversity, including the current state of the art, future perspectives and recommendations. We propose a community road map to generate and more fully exploit multi-omics data from polar organisms. These data are needed for the comprehensive evaluation of polar biodiversity and to reveal how life evolved and adapted to permanently cold environments with extreme seasonality. We argue that concerted action is required to mitigate the impact of warming on polar ecosystems via conservation efforts, to sustainably manage these unique habitats and their ecosystem services, and for the sustainable bioprospecting of novel genes and compounds for societal gain.

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.

Fast and robust metagenomic sequence comparison through sparse chaining with skani

Sequence comparison tools for metagenome-assembled genomes (MAGs) struggle with high-volume or low-quality data. We present skani ( https://github.com/bluenote-1577/skani ), a method for determining average nucleotide identity (ANI) via sparse approximate alignments. skani outperforms FastANI in accuracy and speed (>20× faster) for fragmented, incomplete MAGs. skani can query genomes against >65,000 prokaryotic genomes in seconds and 6 GB memory. skani unlocks higher-resolution insights for extensive, noisy metagenomic datasets.

Alignment-based Protein Mutational Landscape Prediction: Doing More with Less

The wealth of genomic data has boosted the development of computational methods predicting the phenotypic outcomes of missense variants. The most accurate ones exploit multiple sequence alignments, which can be costly to generate. Recent efforts for democratizing protein structure prediction have overcome this bottleneck by leveraging the fast homology search of MMseqs2. Here, we show the usefulness of this strategy for mutational outcome prediction through a large-scale assessment of 1.5M missense variants across 72 protein families. Our study demonstrates the feasibility of producing alignment-based mutational landscape predictions that are both high-quality and compute-efficient for entire proteomes. We provide the community with the whole human proteome mutational landscape and simplified access to our predictive pipeline.

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.

Developing model systems for dinoflagellates in the post-genomic era

Dinoflagellates are a diverse group of eukaryotic microbes that are ubiquitous in aquatic environments. Largely photosynthetic, they encompass symbiotic, parasitic, and free-living lineages with a broad spectrum of trophism. Many free-living taxa can produce bioactive secondary metabolites such as biotoxins, some of which cause harmful algal blooms. In contrast, most symbiotic species are crucial for sustaining coral reef health. The year 2023 marked a decade since the first genome data of dinoflagellates became available. The growing genome-scale resources for these taxa are highlighting their remarkable evolutionary and genomic complexities. Here, we discuss the prospect of developing dinoflagellate models using the criteria of accessibility, tractability, resources, research support, and promise. Moving forward in the post-genomic era, we argue for the development of fit-to-purpose models that tailor to specific biological contexts, and that a one-size-fits-all model is inadequate for encapsulating the complex biology, ecology, and evolutionary history of dinoflagellates.

Unraveling the functional dark matter through global metagenomics

Metagenomes encode an enormous diversity of proteins, reflecting a multiplicity of functions and activities1,2. Exploration of this vast sequence space has been limited to a comparative analysis against reference microbial genomes and protein families derived from those genomes. Here, to examine the scale of yet untapped functional diversity beyond what is currently possible through the lens of reference genomes, we develop a computational approach to generate reference-free protein families from the sequence space in metagenomes. We analyse 26,931 metagenomes and identify 1.17 billion protein sequences longer than 35 amino acids with no similarity to any sequences from 102,491 reference genomes or the Pfam database3. Using massively parallel graph-based clustering, we group these proteins into 106,198 novel sequence clusters with more than 100 members, doubling the number of protein families obtained from the reference genomes clustered using the same approach. We annotate these families on the basis of their taxonomic, habitat, geographical and gene neighbourhood distributions and, where sufficient sequence diversity is available, predict protein three-dimensional models, revealing novel structures. Overall, our results uncover an enormously diverse functional space, highlighting the importance of further exploring the microbial functional dark matter.

ACR: metagenome-assembled prokaryotic and eukaryotic genome refinement tool

Microbial genome recovery from metagenomes can further explain microbial ecosystem structures, functions and dynamics. Thus, this study developed the Additional Clustering Refiner (ACR) to enhance high-purity prokaryotic and eukaryotic metagenome-assembled genome (MAGs) recovery. ACR refines low-quality MAGs by subjecting them to iterative k-means clustering predicated on contig abundance and increasing bin purity through validated universal marker genes. Synthetic and real-world metagenomic datasets, including short- and long-read sequences, evaluated ACR's effectiveness. The results demonstrated improved MAG purity and a significant increase in high- and medium-quality MAG recovery rates. In addition, ACR seamlessly integrates with various binning algorithms, augmenting their strengths without modifying core features. Furthermore, its multiple sequencing technology compatibilities expand its applicability. By efficiently recovering high-quality prokaryotic and eukaryotic genomes, ACR is a promising tool for deepening our understanding of microbial communities through genome-centric metagenomics.

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.

Long-read assembled metagenomic approaches improve our understanding on metabolic potentials of microbial community in mangrove sediments

BACKGROUND

Mangrove wetlands are coastal ecosystems with important ecological features and provide habitats for diverse microorganisms with key roles in nutrient and biogeochemical cycling. However, the overall metabolic potentials and ecological roles of microbial community in mangrove sediment are remained unanswered. In current study, the microbial and metabolic profiles of prokaryotic and fungal communities in mangrove sediments were investigated using metagenomic analysis based on PacBio single-molecule real time (SMRT) and Illumina sequencing techniques.

RESULTS

Comparing to Illumina short reads, the incorporation of PacBio long reads significantly contributed to more contiguous assemblies, yielded more than doubled high-quality metagenome-assembled genomes (MAGs), and improved the novelty of the MAGs. Further metabolic reconstruction for recovered MAGs showed that prokaryotes potentially played an essential role in carbon cycling in mangrove sediment, displaying versatile metabolic potential for degrading organic carbons, fermentation, autotrophy, and carbon fixation. Mangrove fungi also functioned as a player in carbon cycling, potentially involved in the degradation of various carbohydrate and peptide substrates. Notably, a new candidate bacterial phylum named as Candidatus Cosmopoliota with a ubiquitous distribution is proposed. Genomic analysis revealed that this new phylum is capable of utilizing various types of organic substrates, anaerobic fermentation, and carbon fixation with the Wood-Ljungdahl (WL) pathway and the reverse tricarboxylic acid (rTCA) cycle.

CONCLUSIONS

The study not only highlights the advantages of HiSeq-PacBio Hybrid assembly for a more complete profiling of environmental microbiomes but also expands our understanding of the microbial diversity and potential roles of distinct microbial groups in biogeochemical cycling in mangrove sediment. Video Abstract.

MGnify Genomes: A Resource for Biome-specific Microbial Genome Catalogues

An increasingly common output arising from the analysis of shotgun metagenomic datasets is the generation of metagenome-assembled genomes (MAGs), with tens of thousands of MAGs now described in the literature. However, the discovery and comparison of these MAG collections is hampered by the lack of uniformity in their generation, annotation and storage. To address this, we have developed MGnify Genomes, a growing collection of biome-specific non-redundant microbial genome catalogues generated using MAGs and publicly available isolate genomes. Genomes within a biome-specific catalogue are organised into species clusters. For species that contain multiple conspecific genomes, the highest quality genome is selected as the representative, always prioritising an isolate genome over a MAG. The species representative sequences and annotations can be visualised on the MGnify website and the full catalogue and associated analysis outputs can be downloaded from MGnify servers. A suite of online search tools is provided allowing users to compare their own sequences, ranging from a gene to sets of genomes, against the catalogues. Seven biomes are available currently, comprising over 300,000 genomes that represent 11,048 non-redundant species, and include 36 taxonomic classes not currently represented by cultured genomes. MGnify Genomes is available at https://www.ebi.ac.uk/metagenomics/browse/genomes/.

Petascale Homology Search for Structure Prediction

The recent CASP15 competition highlighted the critical role of multiple sequence alignments (MSAs) in protein structure prediction, as demonstrated by the success of the top AlphaFold2-based prediction methods. To push the boundaries of MSA utilization, we conducted a petabase-scale search of the Sequence Read Archive (SRA), resulting in gigabytes of aligned homologs for CASP15 targets. These were merged with default MSAs produced by ColabFold-search and provided to ColabFold-predict. By using SRA data, we achieved highly accurate predictions (GDT_TS > 70) for 66% of the non-easy targets, whereas using ColabFold-search default MSAs scored highly in only 52%. Next, we tested the effect of deep homology search and ColabFold's advanced features, such as more recycles, on prediction accuracy. While SRA homologs were most significant for improving ColabFold's CASP15 ranking from 11th to 3rd place, other strategies contributed too. We analyze these in the context of existing strategies to improve prediction.

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.

Insight into the phylogeny and metabolic divergence of Monascus species (M. pilosus, M. ruber, and M. purpureus) at the genome level

Background

Species of the genus Monascus are economically important and widely used in the production of food colorants and monacolin K. However, they have also been known to produce the mycotoxin citrinin. Currently, taxonomic knowledge of this species at the genome level is insufficient.

Methods

This study presents genomic similarity analyses through the analysis of the average nucleic acid identity of the genomic sequence and the whole genome alignment. Subsequently, the study constructed a pangenome of Monascus by reannotating all the genomes and identifying a total of 9,539 orthologous gene families. Two phylogenetic trees were constructed based on 4,589 single copy orthologous protein sequences and all the 5,565 orthologous proteins, respectively. In addition, carbohydrate active enzymes, secretome, allergic proteins, as well as secondary metabolite gene clusters were compared among the included 15 Monascus strains.

Results

The results clearly revealed a high homology between M. pilosus and M. ruber, and their distant relationship with M. purpureus. Accordingly, all the included 15 Monascus strains should be classified into two distinctly evolutionary clades, namely the M. purpureus clade and the M. pilosus-M. ruber clade. Moreover, gene ontology enrichment showed that the M. pilosus-M. ruber clade had more orthologous genes involved with environmental adaptation than the M. purpureus clade. Compared to Aspergillus oryzae, all the Monascus species had a substantial gene loss of carbohydrate active enzymes. Potential allergenic and fungal virulence factor proteins were also found in the secretome of Monascus. Furthermore, this study identified the pigment synthesis gene clusters present in all included genomes, but with multiple nonessential genes inserted in the gene cluster of M. pilosus and M. ruber compared to M. purpureus. The citrinin gene cluster was found to be intact and highly conserved only among M. purpureus genomes. The monacolin K gene cluster was found only in the genomes of M. pilosus and M. ruber, but the sequence was more conserved in M. ruber.

Conclusion

This study provides a paradigm for phylogenetic analysis of the genus Monascus, and it is believed that this report will lead to a better understanding of these food microorganisms in terms of classification, metabolic differentiation, and safety.

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

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

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.

Identifying Eukaryotes and Factors Influencing Their Biogeography in Drinking Water Metagenomes

The biogeography of eukaryotes in drinking water systems is poorly understood relative to that of prokaryotes or viruses, limiting the understanding of their role and management. A challenge with studying complex eukaryotic communities is that metagenomic analysis workflows are currently not as mature as those that focus on prokaryotes or viruses. In this study, we benchmarked different strategies to recover eukaryotic sequences and genomes from metagenomic data and applied the best-performing workflow to explore the factors affecting the relative abundance and diversity of eukaryotic communities in drinking water distribution systems (DWDSs). We developed an ensemble approach exploiting k-mer- and reference-based strategies to improve eukaryotic sequence identification and identified MetaBAT2 as the best-performing binning approach for their clustering. Applying this workflow to the DWDS metagenomes showed that eukaryotic sequences typically constituted small proportions (i.e., <1%) of the overall metagenomic data with higher relative abundances in surface water-fed or chlorinated systems with high residuals. The α and β diversities of eukaryotes were correlated with those of prokaryotic and viral communities, highlighting the common role of environmental/management factors. Finally, a co-occurrence analysis highlighted clusters of eukaryotes whose members' presence and abundance in DWDSs were affected by disinfection strategies, climate conditions, and source water types.

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.

Plankton biogeography in the 21st century and impacts of climate change: advances through genomics

This article summarizes recent advances in our knowledge of plankton biogeography obtained by genomic approaches and the impacts of global warming on it. Large-scale comparison of the genomic content of samples of different plankton size fractions revealed a partitioning of the oceans into genomic provinces and the impact of major oceanic currents on them. By defining ecological niches, these provinces are extrapolated to all oceans, with the exception of the Arctic Ocean. By the end of the 21st century, a major restructuring of these provinces is projected in response to a high emission greenhouse gas scenario over 50% of the surface of the studied oceans. Such a restructuring could lead to a decrease in export production by 4%. Finally, obtaining assembled sequences of a large number of plankton genomes defining this biogeography has allowed to better characterize the genomic content of the provinces and to identify the species structuring them. These genomes similarly enabled a better description of potential future changes of plankton communities under climate change.

High-resolution metagenomic reconstruction of the freshwater spring bloom

BACKGROUND

The phytoplankton spring bloom in freshwater habitats is a complex, recurring, and dynamic ecological spectacle that unfolds at multiple biological scales. Although enormous taxonomic shifts in microbial assemblages during and after the bloom have been reported, genomic information on the microbial community of the spring bloom remains scarce.

RESULTS

We performed a high-resolution spatio-temporal sampling of the spring bloom in a freshwater reservoir and describe a multitude of previously unknown taxa using metagenome-assembled genomes of eukaryotes, prokaryotes, and viruses in combination with a broad array of methodologies. The recovered genomes reveal multiple distributional dynamics for several bacterial groups with progressively increasing stratification. Analyses of abundances of metagenome-assembled genomes in concert with CARD-FISH revealed remarkably similar in situ doubling time estimates for dominant genome-streamlined microbial lineages. Discordance between quantitations of cryptophytes arising from sequence data and microscopic identification suggested the presence of hidden, yet extremely abundant aplastidic cryptophytes that were confirmed by CARD-FISH analyses. Aplastidic cryptophytes are prevalent throughout the water column but have never been considered in prior models of plankton dynamics. We also recovered the first metagenomic-assembled genomes of freshwater protists (a diatom and a haptophyte) along with thousands of giant viral genomic contigs, some of which appeared similar to viruses infecting haptophytes but owing to lack of known representatives, most remained without any indication of their hosts. The contrasting distribution of giant viruses that are present in the entire water column to that of parasitic perkinsids residing largely in deeper waters allows us to propose giant viruses as the biological agents of top-down control and bloom collapse, likely in combination with bottom-up factors like a nutrient limitation.

CONCLUSION

We reconstructed thousands of genomes of microbes and viruses from a freshwater spring bloom and show that such large-scale genome recovery allows tracking of planktonic succession in great detail. However, integration of metagenomic information with other methodologies (e.g., microscopy, CARD-FISH) remains critical to reveal diverse phenomena (e.g., distributional patterns, in situ doubling times) and novel participants (e.g., aplastidic cryptophytes) and to further refine existing ecological models (e.g., factors affecting bloom collapse). This work provides a genomic foundation for future approaches towards a fine-scale characterization of the organisms in relation to the rapidly changing environment during the course of the freshwater spring bloom. Video Abstract.

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

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

A 2-million-year-old ecosystem in Greenland uncovered by environmental DNA

Late Pliocene and Early Pleistocene epochs 3.6 to 0.8 million years ago1 had climates resembling those forecasted under future warming2. Palaeoclimatic records show strong polar amplification with mean annual temperatures of 11-19 °C above contemporary values3,4. The biological communities inhabiting the Arctic during this time remain poorly known because fossils are rare5. Here we report an ancient environmental DNA6 (eDNA) record describing the rich plant and animal assemblages of the Kap København Formation in North Greenland, dated to around two million years ago. The record shows an open boreal forest ecosystem with mixed vegetation of poplar, birch and thuja trees, as well as a variety of Arctic and boreal shrubs and herbs, many of which had not previously been detected at the site from macrofossil and pollen records. The DNA record confirms the presence of hare and mitochondrial DNA from animals including mastodons, reindeer, rodents and geese, all ancestral to their present-day and late Pleistocene relatives. The presence of marine species including horseshoe crab and green algae support a warmer climate than today. The reconstructed ecosystem has no modern analogue. The survival of such ancient eDNA probably relates to its binding to mineral surfaces. Our findings open new areas of genetic research, demonstrating that it is possible to track the ecology and evolution of biological communities from two million years ago using ancient eDNA.