Tara Oceans: towards global ocean ecosystems biology

Mixotrophic growth of a ubiquitous marine diatom

Diatoms are major players in the global carbon cycle, and their metabolism is affected by ocean conditions. Understanding the impact of changing inorganic nutrients in the oceans on diatoms is crucial, given the changes in global carbon dioxide levels. Here, we present a genome-scale metabolic model (iMK1961) for Cylindrotheca closterium, an in silico resource to understand uncharacterized metabolic functions in this ubiquitous diatom. iMK1961 represents the largest diatom metabolic model to date, comprising 1961 open reading frames and 6718 reactions. With iMK1961, we identified the metabolic response signature to cope with drastic changes in growth conditions. Comparing model predictions with Tara Oceans transcriptomics data unraveled C. closterium's metabolism in situ. Unexpectedly, the diatom only grows photoautotrophically in 21% of the sunlit ocean samples, while the majority of the samples indicate a mixotrophic (71%) or, in some cases, even a heterotrophic (8%) lifestyle in the light. Our findings highlight C. closterium's metabolic flexibility and its potential role in global carbon cycling.

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.

Chemosynthesis: a neglected foundation of marine ecology and biogeochemistry

Chemosynthesis is a metabolic process that transfers carbon to the biosphere using reduced compounds. It is well recognised that chemosynthesis occurs in much of the ocean, but it is often thought to be a negligible process compared to photosynthesis. Here we propose that chemosynthesis is the underlying process governing primary production in much of the ocean and suggest that it extends to a much wider range of compounds, microorganisms, and ecosystems than previously thought. In turn, this process has had a central role in controlling marine biogeochemistry, ecology, and carbon budgets across the vast realms of the ocean, from the dawn of life to contemporary times.

Parity in bacterial communities and resistomes: Microplastic and natural organic particles in the Tyrrhenian Sea

Petroleum-based microplastic particles (MPs) are carriers of antimicrobial resistance genes (ARGs) in aquatic environments, influencing the selection and spread of antimicrobial resistance. This research characterized MP and natural organic particle (NOP) bacterial communities and resistomes in the Tyrrhenian Sea, a region impacted by plastic pollution and climate change. MP and NOP bacterial communities were similar but different from the free-living planktonic communities. Likewise, MP and NOP ARG abundances were similar but different (higher) from the planktonic communities. MP and NOP metagenome-assembled genomes contained ARGs associated with mobile genetic elements and exhibited co-occurrence with metal resistance genes. Overall, these findings show that MPs and NOPs harbor potential pathogenic and antimicrobial resistant bacteria, which can aid in the spread of antimicrobial resistance. Further, petroleum-based MPs do not represent novel ecological niches for allochthonous bacteria; rather, they synergize with NOPs, collectively facilitating the spread of antimicrobial resistance in marine ecosystems.

Niche partitioning and plastisphere core microbiomes in the two most plastic polluted zones of the world ocean

Plastics are offering a new niche for microorganisms colonizing their surface, the so-called "plastisphere," in which diversity and community structure remain to be characterized and compared across ocean pelagic regions. Here, we compared the bacterial diversity of microorganisms living on plastic marine debris (PMD) and the surrounding free-living (FL) and organic particle-attached (PA) lifestyles sampled during the Tara expeditions in two of the most plastic polluted zones in the world ocean, i.e., the North Pacific gyre and the Mediterranean Sea. The 16S rRNA gene sequencing analysis confirmed that PMD are a new anthropogenic ocean habitat for marine microbes at the ocean-basin-scale, with clear niche partitioning compared to FL and PA lifestyles. At an ocean-basin-scale, the composition of the plastisphere communities was mainly driven by environmental selection, rather than polymer types or dispersal effect. A plastisphere "core microbiome" could be identified, mainly dominated by Rhodobacteraceae and Cyanobacteria. Predicted functions indicated the dominance of carbon, nitrogen and sulfur metabolisms on PMD that open new questions on the role of the plastisphere in a large number of important ecological processes in the marine ecosystem.

Marine Fungal Diversity and Dynamics in the Gulf of Trieste (Northern Adriatic Sea)

Fungi contribute to different important ecological processes, including decomposition of organic matter and nutrient cycling, but in the marine environment the main factors influencing their diversity and dynamics at the spatial and temporal levels are still largely unclear. In this study, we performed DNA metabarcoding on seawater sampled monthly over a year and a half in the Gulf of Trieste (northern Adriatic Sea), targeting the internal transcribed spacer (ITS) and the 18S rRNA gene regions. The fungal communities were diverse, very dynamic, and belonged predominantly to marine taxa. Samples could be clustered in two groups, mainly based on the high (> 30%) or low relative proportion of the ascomycetes Parengyodontium album, which emerged as a key taxon in this area. Dissolved and particulate organic C:N ratio played important roles in shaping the mycoplankton assemblages, suggesting that differently bioavailable organic matter pools may be utilized by different consortia. The proportion of fungal over total reads was 31% for ITS and 0.7% for 18S. ITS had the highest taxonomic resolution but low power to detect early divergent fungal lineages. Our results on composition, distribution, and environmental drivers extended our knowledge of the structure and function of the mycobiome of coastal waters.

Uncoupled: investigating the lack of correlation between the transcription of putative plastic-degrading genes in the global ocean microbiome and marine plastic pollution

BACKGROUND

Plastic pollution is a severe threat to marine ecosystems. While some microbial enzymes can degrade certain plastics, the ability of the global ocean microbiome to break down diverse environmental plastics remains limited. We employed metatranscriptomic data from an international ocean survey to explore global and regional patterns in microbial plastic degradation potential.

RESULTS

On a global oceanic scale, we found no significant correlation between levels of plastic pollution and the expression of genes encoding enzymes putatively identified as capable of plastic degradation. Even when looking at different regional scales, ocean depth layers, or plastic types, we found no strong or even moderate correlation between plastic pollution and relative abundances of transcripts for enzymes with presumed plastic biodegradation potential. Our data, however, indicate that microorganisms in the Southern Ocean show a higher potential for plastic degradation, making them more appealing candidates for bioprospecting novel plastic-degrading enzymes.

CONCLUSION

Our research contributes to understanding the complex global relationship between plastic pollution and microbial plastic degradation potential. We reveal that the transcription of putative plastic-degrading genes in the global ocean microbiome does not correlate to marine plastic pollution, highlighting the ongoing danger that plastic poses to marine environments threatened by plastic pollution.

Insights into the phage community structure and potential function in silage fermentation

The virus that infects bacteria known as phage, plays a crucial role in the biogeochemical cycling of nutrients. However, the community structure and potential functions of phages in silage fermentation remain largely unexplored. In this study, we utilized viral metagenomics (viromics) to investigate the types, lifestyles, functions, and nutrient utilization patterns of phages in silage. Our findings indicated a high prevalence of annotated phages belonging to Caudovirales and Geplafuvirales, as well as unclassified phages in silage. The predominant host types for these phages were Campylobacterales and Enterobacterales. Virulent phages dominated the silage environment due to their broader range of hosts and enhanced survival capabilities. All identified phages present in silage were found to be non-pathogenic. Although temperate and virulent phages carried distinct genes associated with nutrient cycling processes, the shared genes (prsA) involved in carbon metabolism underscore the potential significance of phages in regulating carbon metabolism in silage. Overall, our findings provide a valuable foundation for further exploring the complex interactions between phages and microorganisms in regulating silage fermentation quality.

Datascape: exploring heterogeneous dataspace

Data science is a powerful field for gaining insights, comparing, and predicting behaviors from datasets. However, the diversity of methods and hypotheses needed to abstract a dataset exhibits a lack of genericity. Moreover, the shape of a dataset, which structures its contained information and uncertainties, is rarely considered. Inspired by state-of-the-art manifold learning and hull estimations algorithms, we propose a novel framework, the datascape, that leverages topology and graph theory to abstract heterogeneous datasets. Built upon the combination of a nearest neighbor graph, a set of convex hulls, and a metric distance that respects the shape of the data, the datascape allows exploration of the dataset's underlying space. We show that the datascape can uncover underlying functions from simulated datasets, build predictive algorithms with performance close to state-of-the-art algorithms, and reveal insightful geodesic paths between points. It demonstrates versatility through ecological, medical, and simulated data use cases.

Biogeographical and biodiversity patterns of planktonic microeukaryotes along the tropical western to eastern Pacific Ocean transect revealed by metabarcoding

Microeukaryotic plankton (0.2-200 µm), which are morphologically and genetically highly diverse, play a crucial role in ocean productivity and carbon consumption. The Pacific Ocean (PO), one of the world's largest oligotrophic regions, remains largely unexplored in terms of the biogeography and biodiversity of microeukaryotes based on large-scale sampling. We investigated the horizontal distribution of microeukaryotes along a 16,000 km transect from the west to the east of the PO. The alpha diversity indices showed a distinct decreasing trend from west to east, which was highly correlated with water temperature. The microeukaryotic community, which was clustered into the western, central, and eastern PO groups, displayed a significant distance-decay relationship. Syndiniales, a lineage of parasitic dinoflagellates, was ubiquitously distributed along the transect and dominated the community in terms of both sequence and zero-radius operational taxonomic unit (ZOTU) proportions. The prevailing dominance of Syndiniales-affiliated ZOTUs and their close associations with dinoflagellates, diatoms, and radiolarians, as revealed by SparCC correlation analysis, suggested that parasitism may be an important trophic strategy in the surface waters of the PO. Geographical distance and temperature were the most important environmental factors that significantly correlated with community structure. Overall, our study sheds more light on the distribution pattern of both alpha and beta diversities of microeukaryotic communities and highlighted the importance of parasitisms by Syndiniales across the tropical PO.IMPORTANCEUnderstanding the biogeographical and biodiversity patterns of microeukaryotic communities is essential to comprehending their roles in biogeochemical cycling. In this study, planktonic microeukaryotes were collected along a west-to-east Pacific Ocean transect (ca. 16,000 km). Our study revealed that the alpha diversity indices were highly correlated with water temperature, and the microeukaryotic communities displayed a distinct geographical distance-driven pattern. The predominance of the parasitic dinoflagellate lineage Syndiniales and their close relationship with other microeukaryotic groups suggest that parasitism may be a crucial survival strategy for microeukaryotes in the surface waters of the Pacific Ocean. Our findings expand our understanding of the biodiversity and biogeographical pattern of microeukaryotes and highlight the significance of parasitic Syndiniales in the surface ocean.

Extensive intragenomic variations of the 18S rDNA V4 region in the toxigenic diatom species Pseudo-nitzschia multistriata revealed through high-throughput sequencing

Metabarcoding analysis is an effective technique for monitoring the domoic acid-producing Pseudo-nitzschia species in marine environments, uncovering high-levels of molecular diversity. However, such efforts may result in the overinterpretation of Pseudo-nitzschia species diversity, as molecular diversity not only encompasses interspecies and intraspecies diversities but also exhibits extensive intragenomic variations (IGVs). In this study, we analyzed the V4 region of the 18S rDNA of 30 strains of Pseudo-nitzschia multistriata collected from the coasts of China. The results showed that each P. multistriata strain harbored about a hundred of unique 18S rDNA V4 sequence varieties, of which each represented by a unique amplicon sequence variant (ASV). This study demonstrated the extensive degree of IGVs in P. multistriata strains, suggesting that IGVs may also present in other Pseudo-nitzschia species and other phytoplankton species. Understanding the scope and levels of IGVs is crucial for accurately interpreting the results of metabarcoding analysis.

Succession of diversity, assembly mechanisms, and activities of the microeukaryotic community throughout Scrippsiella acuminata (Dinophyceae) bloom phases

Harmful algal bloom (HAB) is a rapidly expanding marine ecological hazard. Although numerous studies have been carried out about the ecological impact and the ecological mechanism of HAB outbreaks, few studies have comprehensively addressed the shifts of species composition, metabolic activity level, driving factors and community assembly mechanisms of microeukaryotic plankton in the course of the bloom event. To fill the gap of research, we conducted 18S ribosomal DNA and RNA sequencing during the initiation, development, sustenance and decline stages of a Scrippsiella acuminata (S. acuminata) bloom at the coastal sea of Fujian Province, China. We found that the bloom event caused a decrease in microeukaryotic plankton species diversity and increase in community homogeneity. Our results revealed that the RNA- and DNA-inferred communities were similar, but α-diversity was more dynamic in RNA- than in DNA-inferred communities. The main taxa with high projected metabolic activity (with RNA:DNA ratio as the proxy) during the bloom included dinoflagellates, Cercozoa, Chlorophyta, Protalveolata, and diatoms. The role of deterministic processes in microeukaryotic plankton community assembly increased during the bloom, but stochastic processes were always the dominant assembly mechanism throughout the bloom process. Our findings improve the understanding of temporal patterns, driving factors and assembly mechanisms underlying the microeukarytic plankton community in a dinoflagellate bloom.

Genome-scale community modelling reveals conserved metabolic cross-feedings in epipelagic bacterioplankton communities

Marine microorganisms form complex communities of interacting organisms that influence central ecosystem functions in the ocean such as primary production and nutrient cycling. Identifying the mechanisms controlling their assembly and activities is a major challenge in microbial ecology. Here, we integrated Tara Oceans meta-omics data to predict genome-scale community interactions within prokaryotic assemblages in the euphotic ocean. A global genome-resolved co-activity network revealed a significant number of inter-lineage associations across diverse phylogenetic distances. Identified co-active communities include species displaying smaller genomes but encoding a higher potential for quorum sensing, biofilm formation, and secondary metabolism. Community metabolic modelling reveals a higher potential for interaction within co-active communities and points towards conserved metabolic cross-feedings, in particular of specific amino acids and group B vitamins. Our integrated ecological and metabolic modelling approach suggests that genome streamlining and metabolic auxotrophies may act as joint mechanisms shaping bacterioplankton community assembly in the global ocean surface.

Increasing transposase abundance with ocean depth correlates with a particle-associated lifestyle

Transposases are mobile genetic elements that move within and between genomes, promoting genomic plasticity in microorganisms. In marine microbial communities, the abundance of transposases increases with depth, but the reasons behind this trend remain unclear. Our analysis of metagenomes from the Tara Oceans and Malaspina Expeditions suggests that a particle-associated lifestyle is the main covariate for the high occurrence of transposases in the deep ocean, and this trend holds true for individual genomes as well as in a community-wide sense. We observed a strong and depth-independent correlation between transposase abundance and the presence of biofilm-associated genes, as well as the prevalence of secretory enzymes. This suggests that mobile genetic elements readily propagate among microbial communities within crowded biofilms. Furthermore, we show that particle association positively correlates with larger genome size, which is in turn associated with higher transposase abundance. Cassette sequences associated with transposons are enriched with genes related to defense mechanisms, which are more highly expressed in the deep sea. Thus, while transposons spread at the expense of their microbial hosts, they also introduce novel genes and potentially benefit the hosts in helping to compete for limited resources. Overall, our results suggest a new understanding of deep ocean particles as highways for gene sharing among defensively oriented microbial genomes.IMPORTANCEGenes can move within and between microbial genomes via mobile genetic elements, which include transposases and transposons. In the oceans, there is a puzzling increase in transposase abundance in microbial genomes as depth increases. To gain insight into this trend, we conducted an extensive analysis of marine microbial metagenomes and metatranscriptomes. We found a significant correlation between transposase abundance and a particle-associated lifestyle among marine microbes at both the metagenome and genome-resolved levels. We also observed a link between transposase abundance and genes related to defense mechanisms. These results suggest that as microbes become densely packed into crowded particles, mobile genes are more likely to spread and carry genetic material that provides a competitive advantage in crowded habitats. This may enable deep sea microbes to effectively compete in such environments.

Forty-nine metagenomic-assembled genomes from an aquatic virome expand Caudoviricetes by 45 potential new families and the newly uncovered Gossevirus of Bamfordvirae

Twenty complete genomes (29-63 kb) and 29 genomes with an estimated completeness of over 90 % (30-90 kb) were identified for novel dsDNA viruses in the Yangshan Harbor metavirome. These newly discovered viruses contribute to the expansion of viral taxonomy by introducing 46 potential new families. Except for one virus, all others belong to the class Caudoviricetes. The exception is a novel member of the recently characterized viral group known as Gossevirus. Fifteen viruses were predicted to be temperate. The predicted hosts for the viruses appear to be involved in various aspects of the nitrogen cycle, including nitrogen fixation, oxidation and denitrification. Two viruses were identified to have a host of Flavobacterium and Tepidimonas fonticaldi, respectively, by matching CRISPR spacers with viral protospacers. Our findings provide an overview for characterizing and identifying specific viruses from Yangshan Harbor. The Gossevirus-like virus uncovered emphasizes the need for further comprehensive isolation and investigation of polinton-like viruses.

A genome and gene catalog of the aquatic microbiomes of the Tibetan Plateau

The Tibetan Plateau supplies water to nearly 2 billion people in Asia, but climate change poses threats to its aquatic microbial resources. Here, we construct the Tibetan Plateau Microbial Catalog by sequencing 498 metagenomes from six water ecosystems (saline lakes, freshwater lakes, rivers, hot springs, wetlands and glaciers). Our catalog expands knowledge of regional genomic diversity by presenting 32,355 metagenome-assembled genomes that de-replicated into 10,723 representative genome-based species, of which 88% were unannotated. The catalog contains nearly 300 million non-redundant gene clusters, of which 15% novel, and 73,864 biosynthetic gene clusters, of which 50% novel, thus expanding known functional diversity. Using these data, we investigate the Tibetan Plateau aquatic microbiome's biogeography along a distance of 2,500 km and >5 km in altitude. Microbial compositional similarity and the shared gene count with the Tibetan Plateau microbiome decline along with distance and altitude difference, suggesting a dispersal pattern. The Tibetan Plateau Microbial Catalog stands as a substantial repository for high-altitude aquatic microbiome resources, providing potential for discovering novel lineages and functions, and bridging knowledge gaps in microbiome biogeography.

Indexing and real-time user-friendly queries in terabyte-sized complex genomic datasets with kmindex and ORA

Public sequencing databases contain vast amounts of biological information, yet they are largely underutilized as it is challenging to efficiently search them for any sequence(s) of interest. We present kmindex, an approach that can index thousands of metagenomes and perform sequence searches in a fraction of a second. The index construction is an order of magnitude faster than previous methods, while search times are two orders of magnitude faster. With negligible false positive rates below 0.01%, kmindex outperforms the precision of existing approaches by four orders of magnitude. Here we demonstrate the scalability of kmindex by successfully indexing 1,393 marine seawater metagenome samples from the Tara Oceans project. Additionally, we introduce the publicly accessible web server Ocean Read Atlas, which enables real-time queries on the Tara Oceans dataset.

The role of microbial biofilms in range shifts of marine habitat-forming organisms

Marine species, such as corals and kelp, are responding to climate change by altering their distributions. Microbial biofilms underpin key processes that affect the establishment, maintenance, and function of these dominant habitat-formers. Climate-mediated changes to microbial biofilms can therefore strongly influence species' range shifts. Here, we review emerging research on the interactions between benthic biofilms and habitat-formers and identify two key areas of interaction where climate change can impact this dynamic: (i) via direct effects on biofilm composition, and (ii) via impacts on the complex feedback loops which exist between the biofilm microbes and habitat-forming organisms. We propose that these key interactions will be fundamental in driving the speed and extent of tropicalisation of coastal ecosystems under climate change.

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.

Molecular mechanisms and evolutionary robustness of a color switch in proteorhodopsins

Proteorhodopsins are widely distributed photoreceptors from marine bacteria. Their discovery revealed a high degree of evolutionary adaptation to ambient light, resulting in blue- and green-absorbing variants that correlate with a conserved glutamine/leucine at position 105. On the basis of an integrated approach combining sensitivity-enhanced solid-state nuclear magnetic resonance (ssNMR) spectroscopy and linear-scaling quantum mechanics/molecular mechanics (QM/MM) methods, this single residue is shown to be responsible for a variety of synergistically coupled structural and electrostatic changes along the retinal polyene chain, ionone ring, and within the binding pocket. They collectively explain the observed color shift. Furthermore, analysis of the differences in chemical shift between nuclei within the same residues in green and blue proteorhodopsins also reveals a correlation with the respective degree of conservation. Our data show that the highly conserved color change mainly affects other highly conserved residues, illustrating a high degree of robustness of the color phenotype to sequence variation.

A photothermal MoS2 decorated biomass carbon-based aerogel with a directionally aligned porous structure for mitigating heavy metal stress under seawater acidification

Marine animals and human are threatened by seawater acidification and metal contamination. Especially, the toxicity of copper (Cu) is expected to be boosted with seawater acidification. However, studies on the removal of Cu under seawater acidification are limited for practical applications, owing to obstacles such as instability, secondary contamination, and low adsorption efficiency. In this work, coconut shells were utilized for the synthesis of biomass carbon, which was then decorated with MoS2. A novel porous MoS2/carbon-based aerogel (MCA) with the synergistic effect of photothermal conversion and adsorption was constructed via directional freeze-drying technology. The adsorption properties of MCA were a precise match with Freundlich isotherm and pseudo-second-order kinetic models with a high correlation coefficient (R2) of more than 0.995. Under solar illumination, the surface temperature of MCA reached up to 36.3 °C and the adsorption capacity of MCA increased to 833.8 mg g-1, indicating that the remarkable thermal properties of MCA contributed to achieving high adsorption capacity. The adsorption mechanisms of MCA involved in the removal of Cu(ii) ions were dominated by chemisorption rather than surface physical adsorption. Owing to its outstanding photothermal conversion performance and directionally aligned porous structure, MCA was able to remove Cu(ii) species from seawater, and the adsorption ability of MCA reached 247.1 mg g-1 after ten adsorption cycles. MCA exhibited excellent stability to resist the complex natural environment and was easy to reuse. Overall, MCA with a series of merits, including high adsorption efficiency, excellent photothermal conversion property, and outstanding cycling stability, was confirmed to contribute to addressing heavy metal stress under seawater acidification.

The road forward to incorporate seawater microbes in predictive reef monitoring

Marine bacterioplankton underpin the health and function of coral reefs and respond in a rapid and sensitive manner to environmental changes that affect reef ecosystem stability. Numerous meta-omics surveys over recent years have documented persistent associations of opportunistic seawater microbial taxa, and their associated functions, with metrics of environmental stress and poor reef health (e.g. elevated temperature, nutrient loads and macroalgae cover). Through positive feedback mechanisms, disturbance-triggered heterotrophic activity of seawater microbes is hypothesised to drive keystone benthic organisms towards the limit of their resilience and translate into shifts in biogeochemical cycles which influence marine food webs, ultimately affecting entire reef ecosystems. However, despite nearly two decades of work in this space, a major limitation to using seawater microbes in reef monitoring is a lack of a unified and focused approach that would move beyond the indicator discovery phase and towards the development of rapid microbial indicator assays for (near) real-time reef management and decision-making. By reviewing the current state of knowledge, we provide a comprehensive framework (defined as five phases of research and innovation) to catalyse a shift from fundamental to applied research, allowing us to move from descriptive to predictive reef monitoring, and from reactive to proactive reef management.

Inferring strain-level mutational drivers of phage-bacteria interaction phenotypes

The enormous diversity of bacteriophages and their bacterial hosts presents a significant challenge to predict which phages infect a focal set of bacteria. Infection is largely determined by complementary -and largely uncharacterized- genetics of adsorption, injection, and cell take-over. Here we present a machine learning (ML) approach to predict phage-bacteria interactions trained on genome sequences of and phenotypic interactions amongst 51 Escherichia coli strains and 45 phage λ strains that coevolved in laboratory conditions for 37 days. Leveraging multiple inference strategies and without a priori knowledge of driver mutations, this framework predicts both who infects whom and the quantitative levels of infections across a suite of 2,295 potential interactions. The most effective ML approach inferred interaction phenotypes from independent contributions from phage and bacteria mutations, predicting phage host range with 86% mean classification accuracy while reducing the relative error in the estimated strength of the infection phenotype by 40%. Further, transparent feature selection in the predictive model revealed 18 of 176 phage λ and 6 of 18 E. coli mutations that have a significant influence on the outcome of phage-bacteria interactions, corroborating sites previously known to affect phage λ infections, as well as identifying mutations in genes of unknown function not previously shown to influence bacterial resistance. While the genetic variation studied was limited to a focal, coevolved phage-bacteria system, the method's success at recapitulating strain-level infection outcomes provides a path forward towards developing strategies for inferring interactions in non-model systems, including those of therapeutic significance.

Diversities and interactions of phages and bacteria in deep-sea sediments as revealed by metagenomics

Phages are found virtually everywhere, even in extreme environments, and are extremely diverse both in their virion structures and in their genomic content. They are thought to shape the taxonomic and functional composition of microbial communities as well as their stability. A number of studies on laboratory culture and viral metagenomic research provide deeper insights into the abundance, diversity, distribution, and interaction with hosts of phages across a wide range of ecosystems. Although most of these studies focus on easily accessible samples, such as soils, lakes, and shallow oceans, little is known about bathypelagic phages. In this study, through analyzing the 16S rRNA sequencing and viral metagenomic sequencing data of 25 samples collected from five different bathypelagic ecosystems, we detected a high diversity of bacteria and phages, particularly in the cold seep and hydrothermal vent ecosystems, which have stable chemical energy. The relative abundance of phages in these ecosystems was higher than in other three abyssal ecosystems. The low phage/host ratios obtained from host prediction were different from shallow ecosystems and indicated the prevalence of prophages, suggesting the complexity of phage-bacteria interactions in abyssal ecosystems. In the correlation analysis, we revealed several phages-bacteria interaction networks of potential ecological relevance. Our study contributes to a better understanding of the interactions between bathypelagic bacteria and their phages.

Vibrio cyclitrophicus phage encoding gene transfer agent fragment, representing a novel viral family

Vibrio is a prevalent bacterial genus in aquatic environments and exhibits diverse metabolic capabilities, playing a vital role in marine biogeochemical cycles. This study isolated a novel virus infecting Vibrio cyclitrophicus, vB_VviC_ZQ26, from coastal waters near Qingdao, China. The vB_VviC_ZQ26 comprises a linear double-stranded DNA genome with a length of 42,982 bp and a G + C content of 43.21 %, encoding 72 putative open reading frames (ORFs). Transmission electron microscope characterization indicates a siphoviral-morphology of vB_VviC_ZQ26. Nucleic-acids-wide analysis indicates a tetranucleotide frequency deviation for genomic segments encoding putative gene transfer agent protein (GTA) and coil-containing protein, implying divergent origins occurred in different parts of viral genomes. Phylogenetic and genome-content-based analysis suggest that vB_VviC_ZQ26 represents a novel vibriophage-specific family designated as Coheviridae. From the result of biogeographic analysis, Coheviridae is mainly colonized in the temperate and tropical epipelagic zones. This study describes a novel vibriophage infecting V. cyclitrophicus, shedding light on the evolutionary divergence of different parts of the viral genome and its ecological footprint in marine environments.

Global Distribution and Diversity of Marine Parmales

Parmales (Bolidophyceae) is a minor eukaryotic phytoplankton group, sister to diatoms, which exists as two distinct forms of unicellular organisms: silicified cells and naked flagellates. Since their discovery, many field studies on Parmales have been performed; however, their global distribution has not yet been examined in detail. We herein compiled more than 3,000 marine DNA metabarcoding datasets targeting the V4 region of the 18S rRNA gene from the EukBank database. By linking this large dataset with the latest morphological and genetic information, we provide updated estimates on the diversity and distribution of Parmales in the global ocean at a fine taxonomic resolution. Parmalean amplicon sequence variants (ASVs) were detected in nearly 90% of the samples analyzed. However, the relative abundance of parmaleans in the eukaryotic community was less than 0.2% on average, and the estimated true richness of parmalean ASVs was approximately 316 ASVs, confirming their low abundance and diversity. A phylogenetic ana-lysis divided these algae into four clades, and three known morphotypes of silicified cells were classified into three different clades. The abundance of Parmales is generally high in the poles and decreases towards the tropics, and individual clades/subclades show further distinctions in their distribution. Collectively, the present results suggest clade/subclade-specific adaptation to different ecological niches.

A Practical Approach to Using the Genomic Standards Consortium MIxS Reporting Standard for Comparative Genomics and Metagenomics

Comparative analysis of (meta)genomes necessitates aggregation, integration, and synthesis of well-annotated data using standards. The Genomic Standards Consortium (GSC) collaborates with the research community to develop and maintain the Minimum Information about any (x) Sequence (MIxS) reporting standard for genomic data. To facilitate the use of the GSC's MIxS reporting standard, we provide a description of the structure and terminology, how to navigate ontologies for required terms in MIxS, and demonstrate practical usage through a soil metagenome example.

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

Background

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

Results

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes were more abundant and species rich than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

Conclusions

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits.

An evolution-based framework for describing human gut bacteria

The human gut microbiome contains many bacterial strains of the same species ('strain-level variants'). Describing strains in a biologically meaningful way rather than purely taxonomically is an important goal but challenging due to the genetic complexity of strain-level variation. Here, we measured patterns of co-evolution across >7,000 strains spanning the bacterial tree-of-life. Using these patterns as a prior for studying hundreds of gut commensal strains that we isolated, sequenced, and metabolically profiled revealed widespread structure beneath the phylogenetic level of species. Defining strains by their co-evolutionary signatures enabled predicting their metabolic phenotypes and engineering consortia from strain genome content alone. Our findings demonstrate a biologically relevant organization to strain-level variation and motivate a new schema for describing bacterial strains based on their evolutionary history.

The diversity and ecological significance of microbial traits potentially involved in B12 biosynthesis in the global ocean

Cobalamin (B12), an essential nutrient and growth cofactor for many living organisms on Earth, can be fully synthesized only by selected prokaryotes in nature. Therefore, microbial communities related to B12 biosynthesis could serve as an example subsystem to disentangle the underlying ecological mechanisms balancing the function and taxonomic make-up of complex functional assemblages. By anchoring microbial traits potentially involved in B12 biosynthesis, we depict the biogeographic patterns of B12 biosynthesis genes and the taxa harboring them in the global ocean, despite the limitations of detecting de novo B12 synthesizers via metagenomes alone. Both the taxonomic and functional composition of B12 biosynthesis genes were strongly shaped by depth, differentiating the epipelagic zones from the mesopelagic layers. Functional genes related to B12 biosynthesis were relatively stably distributed across different oceans, but the taxa harboring them varied considerably, showing clear functional redundancy among microbial systems. Microbial taxa carrying B12 biosynthesis genes in the surface water were influenced by environmental factors such as temperature, oxygen, and nitrate. However, the composition of functional genes was only weakly associated with these environmental factors. Null model analyses demonstrated that determinism governed the variations in B12 biosynthesis genes, whereas a higher degree of stochasticity was associated with taxonomic variations. Significant associations were observed between the chlorophyll a concentration and B12 biosynthesis, confirming its importance in primary production in the global ocean. The results of this study reveal an essential ecological mechanism governing the assembly of microbes in nature: the environment selects for function rather than taxonomy; functional redundancy underlies stochastic community assembly.

Global change ecology: Science to heal a damaged planet

Humanity has drastically altered the biophysical systems that sustain life on Earth. We summarize progress and chart future directions in the emerging field of global change ecology, which studies interactions between organisms and their changing environment.

Surface-Enhanced Raman Scattering-Based Surface Chemotaxonomy: Combining Bacteria Extracellular Matrices and Machine Learning for Rapid and Universal Species Identification

Rapid, universal, and accurate identification of bacteria in their natural states is necessary for on-site environmental monitoring and fundamental microbial research. Surface-enhanced Raman scattering (SERS) spectroscopy emerges as an attractive tool due to its molecule-specific spectral fingerprinting and multiplexing capabilities, as well as portability and speed of readout. Here, we develop a SERS-based surface chemotaxonomy that uses bacterial extracellular matrices (ECMs) as proxy biosignatures to hierarchically classify bacteria based on their shared surface biochemical characteristics to eventually identify six distinct bacterial species at >98% classification accuracy. Corroborating with in silico simulations, we establish a three-way inter-relation between the bacteria identity, their ECM surface characteristics, and their SERS spectral fingerprints. The SERS spectra effectively capture multitiered surface biochemical insights including ensemble surface characteristics, e.g., charge and biochemical profiles, and molecular-level information, e.g., types and numbers of functional groups. Our surface chemotaxonomy thus offers an orthogonal taxonomic definition to traditional classification methods and is achieved without gene amplification, biochemical testing, or specific biomarker recognition, which holds great promise for point-of-need applications and microbial research.

Viral Diversity in Benthic Abyssal Ecosystems: Ecological and Methodological Considerations

Viruses are the most abundant 'biological entities' in the world's oceans. However, technical and methodological constraints limit our understanding of their diversity, particularly in benthic abyssal ecosystems (>4000 m depth). To verify advantages and limitations of analyzing virome DNA subjected either to random amplification or unamplified, we applied shotgun sequencing-by-synthesis to two sample pairs obtained from benthic abyssal sites located in the North-eastern Atlantic Ocean at ca. 4700 m depth. One amplified DNA sample was also subjected to single-molecule long-read sequencing for comparative purposes. Overall, we identified 24,828 viral Operational Taxonomic Units (vOTUs), belonging to 22 viral families. Viral reads were more abundant in the amplified DNA samples (38.5-49.9%) compared to the unamplified ones (4.4-5.8%), with the latter showing a greater viral diversity and 11-16% of dsDNA viruses almost undetectable in the amplified samples. From a procedural point of view, the viromes obtained by direct sequencing (without amplification step) provided a broader overview of both ss and dsDNA viral diversity. Nevertheless, our results suggest that the contextual use of random amplification of the same sample and long-read technology can improve the assessment of viral assemblages by reducing off-target reads.

Discordant patterns between nitrogen-cycling functional traits and taxa in distant coastal sediments reveal important community assembly mechanisms

A central question in microbial ecology is how immense microbes are assembled in changing natural environments while executing critical ecosystem functions. Over the past decade, effort has been made to unravel the contribution of stochasticity and determinism to the compositional of microbial communities. However, most studies focus on microbial taxa, ignoring the importance of functional traits. By employing shotgun metagenomic sequencing and state-of-the-art bioinformatics approaches, this study comprehensively investigated the microbially mediated nitrogen (N) cycling processes in two geographically distant coastal locations. Both shotgun and 16S rRNA gene amplicon sequencing demonstrated significantly differed taxonomic compositions between the two sites. The relative abundance of major microbial phyla, such as Pseudomonadota, Thaumarchaeota, and Bacteroidota, significantly differed. In contrast, high homogeneity was observed for N-cycling functional traits. Statistical analyses suggested that N-cycling taxonomic groups were more related to geographic distance, whereas microbial functional traits were more influenced by environmental factors. Multiple community assembly models demonstrated that determinism strongly governed the microbial N-cycling functional traits, whereas their carrying taxonomic groups were highly stochastic. Such discordant patterns between N-cycling functional traits and taxa demonstrated an important mechanism in microbial ecology in which essential ecosystem functions are stably maintained despite geographic distance and stochastic community assembly.

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.

A tuned ocean color algorithm for the Arctic Ocean: a solution for waters with high CDM content

The Arctic Ocean (AO) is the most river-influenced ocean. Located at the land-sea interface wherein phytoplankton blooms are common, Arctic coastal waterbodies are among the most affected regions by climate change. Given phytoplankton are critical for energy transfer supporting marine food webs, accurate estimation of chlorophyll a concentration (Chl), which is frequently used as a proxy of phytoplankton biomass, is critical for improving our knowledge of the Arctic marine ecosystem and its response to the ongoing climate change. Due to the unique and complex bio-optical properties of the AO, efforts are still needed to obtain more accurate Chl estimates, especially for coastal waters with high colored detrital material (CDM) content. In this study, we optimized the the Garver-Siegel-Maritorena (GSM) algorithm, using an Arctic bio-optical dataset comprised of seven wavelengths (the original GSM wavelengths plus 625 nm). Results suggested that our tuned algorithm, denoted GSMA, outperformed an alternative AO GSM algorithm denoted AO.GSM, but the accuracy of Chl estimates was only improved by 8%. In addition, GSMA showed appreciable robustness when assessed using a satellite image and two non-Arctic coastal datasets.

Challenges and opportunities in sharing microbiome data and analyses

Microbiome data, metadata and analytical workflows have become 'big' in terms of volume and complexity. Although the infrastructure and technologies to share data have been established, the interdisciplinary and multi-omic nature of the field can make resources difficult to identify and use. Following best practices for data deposition requires substantial effort, with sometimes little obvious reward. Gaps remain where microbiome-specific resources for data sharing or reproducibility do not yet exist. We outline available best practices, challenges to their adoption and opportunities in data sharing in microbiome research. We showcase examples of best practices and advocate for their enforcement and incentivization for data sharing. This includes recognition of data curation and sharing endeavours by individuals, institutions, journals and funders. Opportunities for progress include enabling microbiome-specific databases to incorporate future methods for data analysis, integration and reuse.

MicroEXPERT: Microbiome profiling platform with cross-study metagenome-wide association analysis functionality

The framework of the MicroEXPERT platform. Our Platform was composed of five modules. Data management module: Users upload raw data and metadata to the system using a guided workflow. Data processing module: Uploaded data is processed to generate taxonomical distribution and functional composition results. Metagenome-wide association studies module (MWAS): Various methods, including biomarker analysis, PCA, co-occurrence networks, and sample classification, are employed using metadata. Data search module: Users can query nucleotide sequences to retrieve information in the MicroEXPERT database. Data visualization module: Visualization tools are used to illustrate the metagenome analysis results.

Mapping harmful microalgal species by eDNA monitoring: A large-scale survey across the southwestern South China Sea

A large-scale sampling was undertaken during a research cruise across the South China Sea in August 2016, covering an area of about 100,000 km2 to investigate the molecular diversity and distributions of micro-eukaryotic protists, with a focus on the potentially harmful microalgal (HAB) species along the east coast of Peninsular Malaysia. Environmental DNAs from 30 stations were extracted and DNA metabarcoding targeting the V4 and V9 markers in the 18S rDNA was performed. Many protistan molecular units, including previously unreported HAB taxa, were discovered for the first time in the water. Our findings also revealed interesting spatial distribution patterns, with a marked signal of compositional turnover between latitudinal regimes of water masses, where dinophytes and diatom compositions were among the most strongly enhanced at the fronts, leading to distinct niches. Our results further confirmed the widespread distribution of HAB species, such as the toxigenic Alexandrium tamiyavaichii and Pseudo-nitzschia species, and the fish-killing Margalefidinium polykrikoides and Karlodinium veneficum. The molecular information obtained from this study provides an updated HAB species inventory and a toolset that could facilitate existing HAB monitoring schemes in the region to better inform management decisions.

A Pipeline for Constructing Reference Genomes for Large Cohort-Specific Metagenome Compression

Metagenomic data compression is very important as metagenomic projects are facing the challenges of larger data volumes per sample and more samples nowadays. Reference-based compression is a promising method to obtain a high compression ratio. However, existing microbial reference genome databases are not suitable to be directly used as references for compression due to their large size and redundancy, and different metagenomic cohorts often have various microbial compositions. We present a novel pipeline that generated simplified and tailored reference genomes for large metagenomic cohorts, enabling the reference-based compression of metagenomic data. We constructed customized reference genomes, ranging from 2.4 to 3.9 GB, for 29 real metagenomic datasets and evaluated their compression performance. Reference-based compression achieved an impressive compression ratio of over 20 for human whole-genome data and up to 33.8 for all samples, demonstrating a remarkable 4.5 times improvement than the standard Gzip compression. Our method provides new insights into reference-based metagenomic data compression and has a broad application potential for faster and cheaper data transfer, storage, and analysis.

Globally distributed Myxococcota with photosynthesis gene clusters illuminate the origin and evolution of a potentially chimeric lifestyle

Photosynthesis is a fundamental biogeochemical process, thought to be restricted to a few bacterial and eukaryotic phyla. However, understanding the origin and evolution of phototrophic organisms can be impeded and biased by the difficulties of cultivation. Here, we analyzed metagenomic datasets and found potential photosynthetic abilities encoded in the genomes of uncultivated bacteria within the phylum Myxococcota. A putative photosynthesis gene cluster encoding a type-II reaction center appears in at least six Myxococcota families from three classes, suggesting vertical inheritance of these genes from an early common ancestor, with multiple independent losses in other lineages. Analysis of metatranscriptomic datasets indicate that the putative myxococcotal photosynthesis genes are actively expressed in various natural environments. Furthermore, heterologous expression of myxococcotal pigment biosynthesis genes in a purple bacterium supports that the genes can drive photosynthetic processes. Given that predatory abilities are thought to be widespread across Myxococcota, our results suggest the intriguing possibility of a chimeric lifestyle (combining predatory and photosynthetic abilities) in members of this phylum.

Principles and challenges of modeling temporal and spatial omics data

Studies with temporal or spatial resolution are crucial to understand the molecular dynamics and spatial dependencies underlying a biological process or system. With advances in high-throughput omic technologies, time- and space-resolved molecular measurements at scale are increasingly accessible, providing new opportunities to study the role of timing or structure in a wide range of biological questions. At the same time, analyses of the data being generated in the context of spatiotemporal studies entail new challenges that need to be considered, including the need to account for temporal and spatial dependencies and compare them across different scales, biological samples or conditions. In this Review, we provide an overview of common principles and challenges in the analysis of temporal and spatial omics data. We discuss statistical concepts to model temporal and spatial dependencies and highlight opportunities for adapting existing analysis methods to data with temporal and spatial dimensions.

The Microbial Degradation of Natural and Anthropogenic Phosphonates

Phosphonates are compounds containing a direct carbon-phosphorus (C-P) bond, which is particularly resistant to chemical and enzymatic degradation. They are environmentally ubiquitous: some of them are produced by microorganisms and invertebrates, whereas others derive from anthropogenic activities. Because of their chemical stability and potential toxicity, man-made phosphonates pose pollution problems, and many studies have tried to identify biocompatible systems for their elimination. On the other hand, phosphonates are a resource for microorganisms living in environments where the availability of phosphate is limited; thus, bacteria in particular have evolved systems to uptake and catabolize phosphonates. Such systems can be either selective for a narrow subset of compounds or show a broader specificity. The role, distribution, and evolution of microbial genes and enzymes dedicated to phosphonate degradation, as well as their regulation, have been the subjects of substantial studies. At least three enzyme systems have been identified so far, schematically distinguished based on the mechanism by which the C-P bond is ultimately cleaved-i.e., through either a hydrolytic, radical, or oxidative reaction. This review summarizes our current understanding of the molecular systems and pathways that serve to catabolize phosphonates, as well as the regulatory mechanisms that govern their activity.

Top abundant deep ocean heterotrophic bacteria can be retrieved by cultivation

Traditional culture techniques usually retrieve a small fraction of the marine microbial diversity, which mainly belong to the so-called rare biosphere. However, this paradigm has not been fully tested at a broad scale, especially in the deep ocean. Here, we examined the fraction of heterotrophic bacterial communities in photic and deep ocean layers that could be recovered by culture-dependent techniques at a large scale. We compared 16S rRNA gene sequences from a collection of 2003 cultured heterotrophic marine bacteria with global 16S rRNA metabarcoding datasets (16S TAGs) covering surface, mesopelagic and bathypelagic ocean samples that included 16 of the 23 samples used for isolation. These global datasets represent 60 322 unique 16S amplicon sequence variants (ASVs). Our results reveal a significantly higher proportion of isolates identical to ASVs in deeper ocean layers reaching up to 28% of the 16S TAGs of the bathypelagic microbial communities, which included the isolation of 3 of the top 10 most abundant 16S ASVs in the global bathypelagic ocean, related to the genera Sulfitobacter, Halomonas and Erythrobacter. These isolates contributed differently to the prokaryotic communities across different plankton size fractions, recruiting between 38% in the free-living fraction (0.2-0.8 µm) and up to 45% in the largest particles (20-200 µm) in the bathypelagic ocean. Our findings support the hypothesis that sinking particles in the bathypelagic act as resource-rich habitats, suitable for the growth of heterotrophic bacteria with a copiotroph lifestyle that can be cultured, and that these cultivable bacteria can also thrive as free-living bacteria.

Inter-comparison of marine microbiome sampling protocols

Research on marine microbial communities is growing, but studies are hard to compare because of variation in seawater sampling protocols. To help researchers in the inter-comparison of studies that use different seawater sampling methodologies, as well as to help them design future sampling campaigns, we developed the EuroMarine Open Science Exploration initiative (EMOSE). Within the EMOSE framework, we sampled thousands of liters of seawater from a single station in the NW Mediterranean Sea (Service d'Observation du Laboratoire Arago [SOLA], Banyuls-sur-Mer), during one single day. The resulting dataset includes multiple seawater processing approaches, encompassing different material-type kinds of filters (cartridge membrane and flat membrane), three different size fractionations (>0.22 µm, 0.22-3 µm, 3-20 µm and >20 µm), and a number of different seawater volumes ranging from 1 L up to 1000 L. We show that the volume of seawater that is filtered does not have a significant effect on prokaryotic and protist diversity, independently of the sequencing strategy. However, there was a clear difference in alpha and beta diversity between size fractions and between these and "whole water" (with no pre-fractionation). Overall, we recommend care when merging data from datasets that use filters of different pore size, but we consider that the type of filter and volume should not act as confounding variables for the tested sequencing strategies. To the best of our knowledge, this is the first time a publicly available dataset effectively allows for the clarification of the impact of marine microbiome methodological options across a wide range of protocols, including large-scale variations in sampled volume.

Ocean-wide comparisons of mesopelagic planktonic community structures

For decades, marine plankton have been investigated for their capacity to modulate biogeochemical cycles and provide fishery resources. Between the sunlit (epipelagic) layer and the deep dark waters, lies a vast and heterogeneous part of the ocean: the mesopelagic zone. How plankton composition is shaped by environment has been well-explored in the epipelagic but much less in the mesopelagic ocean. Here, we conducted comparative analyses of trans-kingdom community assemblages thriving in the mesopelagic oxygen minimum zone (OMZ), mesopelagic oxic, and their epipelagic counterparts. We identified nine distinct types of intermediate water masses that correlate with variation in mesopelagic community composition. Furthermore, oxygen, NO3- and particle flux together appeared as the main drivers governing these communities. Novel taxonomic signatures emerged from OMZ while a global co-occurrence network analysis showed that about 70% of the abundance of mesopelagic plankton groups is organized into three community modules. One module gathers prokaryotes, pico-eukaryotes and Nucleo-Cytoplasmic Large DNA Viruses (NCLDV) from oxic regions, and the two other modules are enriched in OMZ prokaryotes and OMZ pico-eukaryotes, respectively. We hypothesize that OMZ conditions led to a diversification of ecological niches, and thus communities, due to selective pressure from limited resources. Our study further clarifies the interplay between environmental factors in the mesopelagic oxic and OMZ, and the compositional features of communities.

Terabase-Scale Coassembly of a Tropical Soil Microbiome

Petabases of environmental metagenomic data are publicly available, presenting an opportunity to characterize complex environments and discover novel lineages of life. Metagenome coassembly, in which many metagenomic samples from an environment are simultaneously analyzed to infer the underlying genomes' sequences, is an essential tool for achieving this goal. We applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 terabases (Tbp) of metagenome data from a tropical soil in the Luquillo Experimental Forest (LEF), Puerto Rico. The resulting coassembly yielded 39 high-quality (>90% complete, <5% contaminated, with predicted 23S, 16S, and 5S rRNA genes and ≥18 tRNAs) metagenome-assembled genomes (MAGs), including two from the candidate phylum Eremiobacterota. Another 268 medium-quality (≥50% complete, <10% contaminated) MAGs were extracted, including the candidate phyla Dependentiae, Dormibacterota, and Methylomirabilota. In total, 307 medium- or higher-quality MAGs were assigned to 23 phyla, compared to 294 MAGs assigned to nine phyla in the same samples individually assembled. The low-quality (<50% complete, <10% contaminated) MAGs from the coassembly revealed a 49% complete rare biosphere microbe from the candidate phylum FCPU426 among other low-abundance microbes, an 81% complete fungal genome from the phylum Ascomycota, and 30 partial eukaryotic MAGs with ≥10% completeness, possibly representing protist lineages. A total of 22,254 viruses, many of them low abundance, were identified. Estimation of metagenome coverage and diversity indicates that we may have characterized ≥87.5% of the sequence diversity in this humid tropical soil and indicates the value of future terabase-scale sequencing and coassembly of complex environments. IMPORTANCE Petabases of reads are being produced by environmental metagenome sequencing. An essential step in analyzing these data is metagenome assembly, the computational reconstruction of genome sequences from microbial communities. "Coassembly" of metagenomic sequence data, in which multiple samples are assembled together, enables more complete detection of microbial genomes in an environment than "multiassembly," in which samples are assembled individually. To demonstrate the potential for coassembling terabases of metagenome data to drive biological discovery, we applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 Tbp of reads from a humid tropical soil environment. The resulting coassembly, its functional annotation, and analysis are presented here. The coassembly yielded more, and phylogenetically more diverse, microbial, eukaryotic, and viral genomes than the multiassembly of the same data. Our resource may facilitate the discovery of novel microbial biology in tropical soils and demonstrates the value of terabase-scale metagenome sequencing.

Roles of sulfate-reducing bacteria in sustaining the diversity and stability of marine bacterial community

Microbes play central roles in ocean food webs and global biogeochemical processes. Yet, the information available regarding the highly diverse bacterial communities in these systems is not comprehensive. Here we investigated the diversity, assembly process, and species coexistence frequency of bacterial communities in seawater and sediment across ∼600 km of the eastern Chinese marginal seas using 16S rRNA gene amplicon sequencing. Our analyses showed that compared with seawater, bacterial communities in sediment possessed higher diversity and experienced tight phylogenetic distribution. Neutral model analysis showed that the relative contribution of stochastic processes to the assembly process of bacterial communities in sediment was lower than that in seawater. Functional prediction results showed that sulfate-reducing bacteria (SRB) were enriched in the core bacterial sub-communities. The bacterial diversities of both sediment and seawater were positively associated with the relative abundance of SRB. Co-occurrence analysis showed that bacteria in seawater exhibited a more complex interaction network and closer co-occurrence relationships than those in sediment. The SRB of seawater were centrally located in the network and played an essential role in sustaining the complex network. In addition, further analysis indicated that the SRB of seawater helped maintain the high stability of the bacterial network. Overall, this study provided further comprehensive information regarding the characteristics of bacterial communities in the ocean, and provides new insights into keystone taxa and their roles in sustaining microbial diversity and stability in ocean.

VIRify: An integrated detection, annotation and taxonomic classification pipeline using virus-specific protein profile hidden Markov models

The study of viral communities has revealed the enormous diversity and impact these biological entities have on various ecosystems. These observations have sparked widespread interest in developing computational strategies that support the comprehensive characterisation of viral communities based on sequencing data. Here we introduce VIRify, a new computational pipeline designed to provide a user-friendly and accurate functional and taxonomic characterisation of viral communities. VIRify identifies viral contigs and prophages from metagenomic assemblies and annotates them using a collection of viral profile hidden Markov models (HMMs). These include our manually-curated profile HMMs, which serve as specific taxonomic markers for a wide range of prokaryotic and eukaryotic viral taxa and are thus used to reliably classify viral contigs. We tested VIRify on assemblies from two microbial mock communities, a large metagenomics study, and a collection of publicly available viral genomic sequences from the human gut. The results showed that VIRify could identify sequences from both prokaryotic and eukaryotic viruses, and provided taxonomic classifications from the genus to the family rank with an average accuracy of 86.6%. In addition, VIRify allowed the detection and taxonomic classification of a range of prokaryotic and eukaryotic viruses present in 243 marine metagenomic assemblies. Finally, the use of VIRify led to a large expansion in the number of taxonomically classified human gut viral sequences and the improvement of outdated and shallow taxonomic classifications. Overall, we demonstrate that VIRify is a novel and powerful resource that offers an enhanced capability to detect a broad range of viral contigs and taxonomically classify them.

High wax ester and triacylglycerol biosynthesis potential in coastal sediments of Antarctic and Subantarctic environments

The wax ester (WE) and triacylglycerol (TAG) biosynthetic potential of marine microorganisms is poorly understood at the microbial community level. The goal of this work was to uncover the prevalence and diversity of bacteria with the potential to synthesize these neutral lipids in coastal sediments of two high latitude environments, and to characterize the gene clusters related to this process. Homolog sequences of the key enzyme, the wax ester synthase/acyl-CoA:diacylglycerol acyltransferase (WS/DGAT) were retrieved from 13 metagenomes, including subtidal and intertidal sediments of a Subantarctic environment (Ushuaia Bay, Argentina), and subtidal sediments of an Antarctic environment (Potter Cove, Antarctica). The abundance of WS/DGAT homolog sequences in the sediment metagenomes was 1.23 ± 0.42 times the abundance of 12 single-copy genes encoding ribosomal proteins, higher than in seawater (0.13 ± 0.31 times in 338 metagenomes). Homolog sequences were highly diverse, and were assigned to the Pseudomonadota, Actinomycetota, Bacteroidota and Acidobacteriota phyla. The genomic context of WS/DGAT homologs included sequences related to WE and TAG biosynthesis pathways, as well as to other related pathways such as fatty-acid metabolism, suggesting carbon recycling might drive the flux to neutral lipid synthesis. These results indicate the presence of abundant and taxonomically diverse bacterial populations with the potential to synthesize lipid storage compounds in marine sediments, relating this metabolic process to bacterial survival.