Gene Expression Changes and Community Turnover Differentially Shape the Global Ocean Metatranscriptome

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

Diversity in the utilization of different molecular classes of dissolved organic matter by heterotrophic marine bacteria

Heterotrophic marine bacteria utilize and recycle dissolved organic matter (DOM), impacting biogeochemical cycles. It is currently unclear to what extent distinct DOM components can be used by different heterotrophic clades. Here, we ask how a natural microbial community from the Eastern Mediterranean Sea (EMS) responds to different molecular classes of DOM (peptides, amino acids, amino sugars, disaccharides, monosaccharides, and organic acids) comprising much of the biomass of living organisms. Bulk bacterial activity increased after 24 h for all treatments relative to the control, while glucose and ATP uptake decreased or remained unchanged. Moreover, while the per-cell uptake rate of glucose and ATP decreased, that of Leucin significantly increased for amino acids, reflecting their importance as common metabolic currencies in the marine environment. Pseudoalteromonadaceae dominated the peptides treatment, while different Vibrionaceae strains became dominant in response to amino acids and amino sugars. Marinomonadaceae grew well on organic acids, and Alteromonadaseae on disaccharides. A comparison with a recent laboratory-based study reveals similar peptide preferences for Pseudoalteromonadaceae, while Alteromonadaceae, for example, grew well in the lab on many substrates but dominated in seawater samples only when disaccharides were added. We further demonstrate a potential correlation between the genetic capacity for degrading amino sugars and the dominance of specific clades in these treatments. These results highlight the diversity in DOM utilization among heterotrophic bacteria and complexities in the response of natural communities.

IMPORTANCE

A major goal of microbial ecology is to predict the dynamics of natural communities based on the identity of the organisms, their physiological traits, and their genomes. Our results show that several clades of heterotrophic bacteria each grow in response to one or more specific classes of organic matter. For some clades, but not others, growth in a complex community is similar to that of isolated strains in laboratory monoculture. Additionally, by measuring how the entire community responds to various classes of organic matter, we show that these results are ecologically relevant, and propose that some of these resources are utilized through common uptake pathways. Tracing the path between different resources to the specific microbes that utilize them, and identifying commonalities and differences between different natural communities and between them and lab cultures, is an important step toward understanding microbial community dynamics and predicting how communities will respond to perturbations.

Microorganisms Involved in Methylmercury Demethylation and Mercury Reduction are Widely Distributed and Active in the Bathypelagic Deep Ocean Waters

The ocean's mercury (Hg) content has tripled due to anthropogenic activities, and although the dark ocean (>200 m) has become an important Hg reservoir, concentrations of the toxic and bioaccumulative methylmercury (MeHg) are low and therefore very difficult to measure. As a consequence, the current understanding of the Hg cycle in the deep ocean is severely data-limited, and the factors controlling MeHg, as well as its transformation rates, remain largely unknown. By analyzing 52 globally distributed bathypelagic deep-ocean metagenomes and 26 new metatranscriptomes from the Malaspina Expedition, our study reveals the widespread distribution and expression of bacterial-coding genes merA and merB in the global bathypelagic ocean (∼4000 m depth). These genes, associated with HgII reduction and MeHg demethylation, respectively, are particularly prevalent within the particle-attached fraction. Moreover, our results indicate that water mass age and the organic matter composition shaped the structure of the communities harboring merA and merB genes living in different particle size fractions, their abundance, and their expression levels. Members of the orders Corynebacteriales, Rhodobacterales, Alteromonadales, Oceanospirillales, Moraxellales, and Flavobacteriales were the main taxonomic players containing merA and merB genes in the deep ocean. These findings, together with our previous results of pure culture isolates of the deep bathypelagic ocean possessing the metabolic capacity to degrade MeHg, indicated that both methylmercury demethylation and HgII reduction likely occur in the global dark ocean, the largest biome in the biosphere.

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.

Refocusing the microbial rare biosphere concept through a functional lens

The influential concept of the rare biosphere in microbial ecology has underscored the importance of taxa occurring at low abundances yet potentially playing key roles in communities and ecosystems. Here, we refocus the concept of rare biosphere through a functional trait-based lens and provide a framework to characterize microbial functional rarity, a combination of numerical scarcity across space or time and trait distinctiveness. We demonstrate how this novel interpretation of the rare biosphere, rooted in microbial functions, can enhance our mechanistic understanding of microbial community structure. It also sheds light on functionally distinct microbes, directing conservation efforts towards taxa harboring rare yet ecologically crucial functions.

Hunting for pigments in bacterial settlers of the Great Pacific Garbage Patch

The Great Pacific Garbage Patch, a significant collection of plastic introduced by human activities, provides an ideal environment to study bacterial lifestyles on plastic substrates. We proposed that bacteria colonizing the floating plastic debris would develop strategies to deal with the ultraviolet-exposed substrate, such as the production of antioxidant pigments. We observed a variety of pigmentation in 67 strains that were directly cultivated from plastic pieces sampled from the Garbage Patch. The genomic analysis of four representative strains, each distinct in taxonomy, revealed multiple pathways for carotenoid production. These pathways include those that produce less common carotenoids and a cluster of photosynthetic genes. This cluster appears to originate from a potentially new species of the Rhodobacteraceae family. This represents the first report of an aerobic anoxygenic photoheterotrophic bacterium from plastic biofilms. Spectral analysis showed that the bacteria actively produce carotenoids, such as beta-carotene and beta-cryptoxanthin, and bacteriochlorophyll a. Furthermore, we discovered that the genetic ability to synthesize carotenoids is more common in plastic biofilms than in the surrounding water communities. Our findings suggest that plastic biofilms could be an overlooked source of bacteria-produced carotenoids, including rare forms. It also suggests that photoreactive molecules might play a crucial role in bacterial biofilm communities in surface water.

Challenges in computational discovery of bioactive peptides in 'omics data

Peptides have a plethora of activities in biological systems that can potentially be exploited biotechnologically. Several peptides are used clinically, as well as in industry and agriculture. The increase in available 'omics data has recently provided a large opportunity for mining novel enzymes, biosynthetic gene clusters, and molecules. While these data primarily consist of DNA sequences, other types of data provide important complementary information. Due to their size, the approaches proven successful at discovering novel proteins of canonical size cannot be naïvely applied to the discovery of peptides. Peptides can be encoded directly in the genome as short open reading frames (smORFs), or they can be derived from larger proteins by proteolysis. Both of these peptide classes pose challenges as simple methods for their prediction result in large numbers of false positives. Similarly, functional annotation of larger proteins, traditionally based on sequence similarity to infer orthology and then transferring functions between characterized proteins and uncharacterized ones, cannot be applied for short sequences. The use of these techniques is much more limited and alternative approaches based on machine learning are used instead. Here, we review the limitations of traditional methods as well as the alternative methods that have recently been developed for discovering novel bioactive peptides with a focus on prokaryotic genomes and metagenomes.

Bacterioplankton taxa compete for iron along the early spring-summer transition in the Arctic Ocean

Microbial assemblages under the sea ice of the Dease Strait, Canadian Arctic, were sequenced for metagenomes of a small size fraction (0.2-3 μm). The community from early March was typical for this season, with Alpha- and Gammaproteobacteria as the dominant taxa, followed by Thaumarchaeota and Bacteroidetes. Toward summer, Bacteroidetes, and particularly the genus Polaribacter, became increasingly dominant, followed by the Gammaproteobacteria. Analysis of genes responsible for microbial acquisition of iron showed an abundance of ABC transporters for divalent cations and ferrous iron. The most abundant transporters, however, were the outer membrane TonB-dependent transporters of iron-siderophore complexes. The abundance of iron acquisition genes suggested this element was essential for the microbial assemblage. Interestingly, Gammaproteobacteria were responsible for most of the siderophore synthesis genes. On the contrary, Bacteroidetes did not synthesize siderophores but accounted for most of the transporters, suggesting a role as cheaters in the competition for siderophores as public goods. This cheating ability of the Bacteroidetes may have contributed to their dominance in the summer.

Global biogeography of microbes driving ocean ecological status under climate change

Microbial communities play a crucial role in ocean ecology and global biogeochemical processes. However, understanding the intricate interactions among diversity, taxonomical composition, functional traits, and how these factors respond to climate change remains a significant challenge. Here, we propose seven distinct ecological statuses by systematically considering the diversity, structure, and biogeochemical potential of the ocean microbiome to delineate their biogeography. Anthropogenic climate change is expected to alter the ecological status of the surface ocean by influencing environmental conditions, particularly nutrient and oxygen contents. Our predictive model, which utilizes machine learning, indicates that the ecological status of approximately 32.44% of the surface ocean may undergo changes from the present to the end of this century, assuming no policy interventions. These changes mainly include poleward shifts in the main taxa, increases in photosynthetic carbon fixation and decreases in nutrient metabolism. However, this proportion can decrease significantly with effective control of greenhouse gas emissions. Our study underscores the urgent necessity for implementing policies to mitigate climate change, particularly from an ecological perspective.

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.

Substrate uptake patterns shape niche separation in marine prokaryotic microbiome

Marine heterotrophic prokaryotes primarily take up ambient substrates using transporters. The patterns of transporters targeting particular substrates shape the ecological role of heterotrophic prokaryotes in marine organic matter cycles. Here, we report a size-fractionated pattern in the expression of prokaryotic transporters throughout the oceanic water column due to taxonomic variations, revealed by a multi-"omics" approach targeting ATP-binding cassette (ABC) transporters and TonB-dependent transporters (TBDTs). Substrate specificity analyses showed that marine SAR11, Rhodobacterales, and Oceanospirillales use ABC transporters to take up organic nitrogenous compounds in the free-living fraction, while Alteromonadales, Bacteroidetes, and Sphingomonadales use TBDTs for carbon-rich organic matter and metal chelates on particles. The expression of transporter proteins also supports distinct lifestyles of deep-sea prokaryotes. Our results suggest that transporter divergency in organic matter assimilation reflects a pronounced niche separation in the prokaryote-mediated organic matter cycles.

Identification and validation of microbial biomarkers from cross-cohort datasets using xMarkerFinder

Microbial signatures have emerged as promising biomarkers for disease diagnostics and prognostics, yet their variability across different studies calls for a standardized approach to biomarker research. Therefore, we introduce xMarkerFinder, a four-stage computational framework for microbial biomarker identification with comprehensive validations from cross-cohort datasets, including differential signature identification, model construction, model validation and biomarker interpretation. xMarkerFinder enables the identification and validation of reproducible biomarkers for cross-cohort studies, along with the establishment of classification models and potential microbiome-induced mechanisms. Originally developed for gut microbiome research, xMarkerFinder's adaptable design makes it applicable to various microbial habitats and data types. Distinct from existing biomarker research tools that typically concentrate on a singular aspect, xMarkerFinder uniquely incorporates a sophisticated feature selection process, specifically designed to address the heterogeneity between different cohorts, extensive internal and external validations, and detailed specificity assessments. Execution time varies depending on the sample size, selected algorithm and computational resource. Accessible via GitHub ( https://github.com/tjcadd2020/xMarkerFinder ), xMarkerFinder supports users with diverse expertise levels through different execution options, including step-to-step scripts with detailed tutorials and frequently asked questions, a single-command execution script, a ready-to-use Docker image and a user-friendly web server ( https://www.biosino.org/xmarkerfinder ).

Decoupling between the genetic potential and the metabolic regulation and expression in microbial organic matter cleavage across microbiomes

Metagenomics, metatranscriptomics, and metaproteomics are used to explore the microbial capability of enzyme secretion, but the links between protein-encoding genes and corresponding transcripts/proteins across ecosystems are underexplored. By conducting a multi-omics comparison focusing on key enzymes (carbohydrate-active enzymes [CAZymes] and peptidases) cleaving the main biomolecules across distinct microbiomes living in the ocean, soil, and human gut, we show that the community structure, functional diversity, and secretion mechanisms of microbial secretory CAZymes and peptidases vary drastically between microbiomes at metagenomic, metatranscriptomic, and metaproteomic levels. Such variations lead to decoupled relationships between CAZymes and peptidases from genetic potentials to protein expressions due to the different responses of key players toward organic matter sources and concentrations. Our results highlight the need for systematic analysis of the factors shaping patterns of microbial cleavage on organic matter to better link omics data to ecosystem processes.

IMPORTANCE

Omics tools are used to explore adaptive mechanism of microbes in diverse systems, but the advantages and limitations of different omics tools remain skeptical. Here, we reported distinct profiles in microbial secretory enzyme composition revealed by different omics methods. In general, the predicted function from metagenomic analysis decoupled from the expression of corresponding transcripts/proteins. Linking omics results to taxonomic origin, functional capability, substrate specificity, secretion preference, and enzymatic activity measurement suggested the substrate's source, concentration and stoichiometry impose strong filtering on the expression of extracellular enzymes, which may overwrite the genetic potentials. Our results present an integrated perspective on the need for multi-dimensional characterization of microbial adaptation in a changing environment.

Timing the evolution of phosphorus-cycling enzymes through geological time using phylogenomics

Phosphorus plays a crucial role in controlling biological productivity, but geological estimates of phosphate concentrations in the Precambrian ocean, during life's origin and early evolution, vary over several orders of magnitude. While reduced phosphorus species may have served as alternative substrates to phosphate, their bioavailability on the early Earth remains unknown. Here, we reconstruct the phylogenomic record of life on Earth and find that phosphate transporting genes (pnas) evolved in the Paleoarchean (ca. 3.6-3.2 Ga) and are consistent with phosphate concentrations above modern levels ( > 3 µM). The first gene optimized for low phosphate levels (pstS; <1 µM) appeared around the same time or in the Mesoarchean depending on the reconstruction method. Most enzymatic pathways for metabolising reduced phosphorus emerged and expanded across the tree of life later. This includes phosphonate-catabolising CP-lyases, phosphite-oxidising pathways and hypophosphite-oxidising pathways. CP-lyases are particularly abundant in dissolved phosphate concentrations below 0.1 µM. Our results thus indicate at least local regions of declining phosphate levels through the Archean, possibly linked to phosphate-scavenging Fe(III), which may have limited productivity. However, reduced phosphorus species did not become widely used until after the Paleoproterozoic Great Oxidation Event (2.3 Ga), possibly linked to expansion of the biosphere at that time.

Sulfur oxidation and reduction are coupled to nitrogen fixation in the roots of the salt marsh foundation plant Spartina alterniflora

Heterotrophic activity, primarily driven by sulfate-reducing prokaryotes, has traditionally been linked to nitrogen fixation in the root zone of coastal marine plants, leaving the role of chemolithoautotrophy in this process unexplored. Here, we show that sulfur oxidation coupled to nitrogen fixation is a previously overlooked process providing nitrogen to coastal marine macrophytes. In this study, we recovered 239 metagenome-assembled genomes from a salt marsh dominated by the foundation plant Spartina alterniflora, including diazotrophic sulfate-reducing and sulfur-oxidizing bacteria. Abundant sulfur-oxidizing bacteria encode and highly express genes for carbon fixation (RuBisCO), nitrogen fixation (nifHDK) and sulfur oxidation (oxidative-dsrAB), especially in roots stressed by sulfidic and reduced sediment conditions. Stressed roots exhibited the highest rates of nitrogen fixation and expression level of sulfur oxidation and sulfate reduction genes. Close relatives of marine symbionts from the Candidatus Thiodiazotropha genus contributed ~30% and ~20% of all sulfur-oxidizing dsrA and nitrogen-fixing nifK transcripts in stressed roots, respectively. Based on these findings, we propose that the symbiosis between S. alterniflora and sulfur-oxidizing bacteria is key to ecosystem functioning of coastal salt marshes.

Genomic and transcriptomic characterization of methylmercury detoxification in a deep ocean Alteromonas mediterranea ISS312

Methylmercury (MeHg) is one of the most worrisome pollutants in marine systems. MeHg detoxification is mediated by merB and merA genes, responsible for the demethylation of MeHg and the reduction of inorganic mercury, respectively. Little is known about the biological capacity to detoxify this compound in marine environments, and even less the bacterial transcriptional changes during MeHg detoxification. This study provides the genomic and transcriptomic characterization of the deep ocean bacteria Alteromonas mediterranea ISS312 with capacity for MeHg degradation. Its genome sequence revealed four mer operons containing three merA gene and two merB gene copies, that could be horizontally transferred among distant related genomes by mobile genetic elements. The transcriptomic profiling in the presence of 5 μM MeHg showed that merA and merB genes are within the most expressed genes, although not all mer genes were equally transcribed. Besides, we aimed to identify functional orthologous genes that displayed expression profiles highly similar or identical to those genes within the mer operons, which could indicate they are under the same regulatory controls. We found contrasting expression profiles for each mer operon that were positively correlated with a wide array of functions mostly related to amino acid metabolism, but also to flagellar assembly or two component systems. Also, this study highlights that all merAB genes of the four operons were globally distributed across oceans layers with higher transcriptional activity in the mesopelagic deeper waters. Our study provides new insights about the transcriptional patterns related to the capacity of marine bacteria to detoxify MeHg, with important implications for the understanding of this process in marine ecosystems.

Decoding populations in the ocean microbiome

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

Bacteria contribute exopolysaccharides to an algal-bacterial joint extracellular matrix

Marine ecosystems are influenced by phytoplankton aggregation, which affects processes like marine snow formation and harmful events such as marine mucilage outbreaks. Phytoplankton secrete exopolymers, creating an extracellular matrix (ECM) that promotes particle aggregation. This ECM attracts heterotrophic bacteria, providing a nutrient-rich and protective environment. In terrestrial environments, bacterial colonization near primary producers relies on attachment and the formation of multidimensional structures like biofilms. Bacteria were observed attaching and aggregating within algal-derived exopolymers, but it is unclear if bacteria produce an ECM that contributes to this colonization. This study, using Emiliania huxleyi algae and Phaeobacter inhibens bacteria in an environmentally relevant model system, reveals a shared algal-bacterial ECM scaffold that promotes algal-bacterial aggregation. Algal exudates play a pivotal role in promoting bacterial colonization, stimulating bacterial exopolysaccharide (EPS) production, and facilitating a joint ECM formation. A bacterial biosynthetic pathway responsible for producing a specific EPS contributing to bacterial ECM formation is identified. Genes from this pathway show increased expression in algal-rich environments. These findings highlight the underestimated role of bacteria in aggregate-mediated processes in marine environments, offering insights into algal-bacterial interactions and ECM formation, with implications for understanding and managing natural and perturbed aggregation events.

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.

Hypometabolism to survive the long polar night and subsequent successful return to light in the diatom Fragilariopsis cylindrus

Diatoms, the main eukaryotic phytoplankton of the polar marine regions, are essential for the maintenance of food chains specific to Arctic and Antarctic ecosystems, and are experiencing major disturbances under current climate change. As such, it is fundamental to understand the physiological mechanisms and associated molecular basis of their endurance during the long polar night. Here, using the polar diatom Fragilariopsis cylindrus, we report an integrative analysis combining transcriptomic, microscopic and biochemical approaches to shed light on the strategies used to survive the polar night. We reveal that in prolonged darkness, diatom cells enter a state of quiescence with reduced metabolic and transcriptional activity, during which no cell division occurs. We propose that minimal energy is provided by respiration and degradation of protein, carbohydrate and lipid stores and that homeostasis is maintained by autophagy in prolonged darkness. We also report internal structural changes that manifest the morphological acclimation of cells to darkness, including the appearance of a large vacuole. Our results further show that immediately following a return to light, diatom cells are able to use photoprotective mechanisms and rapidly resume photosynthesis, demonstrating the remarkable robustness of polar diatoms to prolonged darkness at low temperature.

Marine picoplankton metagenomes and MAGs from eleven vertical profiles obtained by the Malaspina Expedition

The Ocean microbiome has a crucial role in Earth's biogeochemical cycles. During the last decade, global cruises such as Tara Oceans and the Malaspina Expedition have expanded our understanding of the diversity and genetic repertoire of marine microbes. Nevertheless, there are still knowledge gaps regarding their diversity patterns throughout depth gradients ranging from the surface to the deep ocean. Here we present a dataset of 76 microbial metagenomes (MProfile) of the picoplankton size fraction (0.2-3.0 µm) collected in 11 vertical profiles covering contrasting ocean regions sampled during the Malaspina Expedition circumnavigation (7 depths, from surface to 4,000 m deep). The MProfile dataset produced 1.66 Tbp of raw DNA sequences from which we derived: 17.4 million genes clustered at 95% sequence similarity (M-GeneDB-VP), 2,672 metagenome-assembled genomes (MAGs) of Archaea and Bacteria (Malaspina-VP-MAGs), and over 100,000 viral genomic sequences. This dataset will be a valuable resource for exploring the functional and taxonomic connectivity between the photic and bathypelagic tropical and sub-tropical ocean, while increasing our general knowledge of the Ocean microbiome.

Functional and evolutionary significance of unknown genes from uncultivated taxa

Many of the Earth's microbes remain uncultured and understudied, limiting our understanding of the functional and evolutionary aspects of their genetic material, which remain largely overlooked in most metagenomic studies1. Here we analysed 149,842 environmental genomes from multiple habitats2-6 and compiled a curated catalogue of 404,085 functionally and evolutionarily significant novel (FESNov) gene families exclusive to uncultivated prokaryotic taxa. All FESNov families span multiple species, exhibit strong signals of purifying selection and qualify as new orthologous groups, thus nearly tripling the number of bacterial and archaeal gene families described to date. The FESNov catalogue is enriched in clade-specific traits, including 1,034 novel families that can distinguish entire uncultivated phyla, classes and orders, probably representing synapomorphies that facilitated their evolutionary divergence. Using genomic context analysis and structural alignments we predicted functional associations for 32.4% of FESNov families, including 4,349 high-confidence associations with important biological processes. These predictions provide a valuable hypothesis-driven framework that we used for experimental validatation of a new gene family involved in cell motility and a novel set of antimicrobial peptides. We also demonstrate that the relative abundance profiles of novel families can discriminate between environments and clinical conditions, leading to the discovery of potentially new biomarkers associated with colorectal cancer. We expect this work to enhance future metagenomics studies and expand our knowledge of the genetic repertory of uncultivated organisms.

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.

Disentangling microbial networks across pelagic zones in the tropical and subtropical global ocean

Microbial interactions are vital in maintaining ocean ecosystem function, yet their dynamic nature and complexity remain largely unexplored. Here, we use association networks to investigate possible ecological interactions in the marine microbiome among archaea, bacteria, and picoeukaryotes throughout different depths and geographical regions of the tropical and subtropical global ocean. Our findings reveal that potential microbial interactions change with depth and geographical scale, exhibiting highly heterogeneous distributions. A few potential interactions were global, meaning they occurred across regions at the same depth, while 11-36% were regional within specific depths. The bathypelagic zone had the lowest proportion of global associations, and regional associations increased with depth. Moreover, we observed that most surface water associations do not persist in deeper ocean layers despite microbial vertical dispersal. Our work contributes to a deeper understanding of the tropical and subtropical global ocean interactome, which is essential for addressing the challenges posed by global change.

Microbial regulation of feedbacks to ecosystem change

Microbes are key biodiversity components of all ecosystems and control vital ecosystem functions. Although we have just begun to unravel the scales and factors that regulate microbial communities, their role in mediating ecosystem stability in response to disturbances remains underexplored. Here, we review evidence of how, when, and where microbes regulate or drive disturbance feedbacks. Negative feedbacks dampen the impacts of disturbance, which maintain ecosystem stability, whereas positive feedbacks instead erode stability by amplifying the disturbance. Here we describe the processes underlying the responses to disturbance using a hierarchy of functional traits, and we exemplify how these may drive biogeochemical feedbacks. We suggest that the feedback potential of functional traits at different hierarchical levels is contingent on the complexity and heterogeneity of the environment.

New chemical and microbial perspectives on vitamin B1 and vitamer dynamics of a coastal system

Vitamin B1 (thiamin, B1) is an essential micronutrient for cells, yet intriguingly in aquatic systems most bacterioplankton are unable to synthesize it de novo (auxotrophy), requiring an exogenous source. Cycling of this valuable metabolite in aquatic systems has not been fully investigated and vitamers (B1-related compounds) have only begun to be measured and incorporated into the B1 cycle. Here, we identify potential key producers and consumers of B1 and gain new insights into the dynamics of B1 cycling through measurements of B1 and vitamers (HMP: 4-amino-5-hydroxymethyl-2-methylpyrimidine, HET: 4-methyl-5-thiazoleethanol, FAMP: N-formyl-4-amino-5-aminomethyl-2-methylpyrimidine) in the particulate and dissolved pool in a temperate coastal system. Dissolved B1 was not the primary limiting nutrient for bacterial production and was relatively stable across seasons with concentrations ranging from 74-117 pM, indicating a balance of supply and demand. However, vitamer concentration changed markedly with season as did transcripts related to vitamer salvage and transport suggesting use of vitamers by certain bacterioplankton, e.g. Pelagibacterales. Genomic and transcriptomic analyses showed that up to 78% of the bacterioplankton taxa were B1 auxotrophs. Notably, de novo B1 production was restricted to a few abundant bacterioplankton (e.g. Vulcanococcus, BACL14 (Burkholderiales), Verrucomicrobiales) across seasons. In summer, abundant picocyanobacteria were important putative B1 sources, based on transcriptional activity, leading to an increase in the B1 pool. Our results provide a new dynamic view of the players and processes involved in B1 cycling over time in coastal waters, and identify specific priority populations and processes for future study.

Carbohydrates and carbohydrate degradation gene abundance and transcription in Atlantic waters of the Arctic

Carbohydrates are chemically and structurally diverse, represent a substantial fraction of marine organic matter and are key substrates for heterotrophic microbes. Studies on carbohydrate utilisation by marine microbes have been centred on phytoplankton blooms in temperate regions, while far less is known from high-latitude waters and during later seasonal stages. Here, we combine glycan microarrays and analytical chromatography with metagenomics and metatranscriptomics to show the spatial heterogeneity in glycan distribution and potential carbohydrate utilisation by microbes in Atlantic waters of the Arctic. The composition and abundance of monomers and glycan structures in POM varied with location and depth. Complex fucose-containing sulfated polysaccharides, known to accumulate in the ocean, were consistently detected, while the more labile β-1,3-glucan exhibited a patchy distribution. Through 'omics analysis, we identify variations in the abundance and transcription of carbohydrate degradation-related genes across samples at the community and population level. The populations contributing the most to transcription were taxonomically related to those known as primary responders and key carbohydrate degraders in temperate ecosystems, such as NS4 Marine Group and Formosa. The unique transcription profiles for these populations suggest distinct substrate utilisation potentials, with predicted glycan targets corresponding to those structurally identified in POM from the same sampling sites. By combining cutting-edge technologies and protocols, we provide insights into the carbohydrate component of the carbon cycle in the Arctic during late summer and present a high-quality dataset that will be of great value for future comparative analyses.

Ecological function and interaction of different bacterial groups during alginate processing in coastal seawater community

The degradation of high molecular weight organic matter (HMWOM) is a core process of oceanic carbon cycle, which is determined by the activity of microbial communities harboring hundreds of different species. Illustrating the active microbes and their interactions during HMWOM processing can provide key information for revealing the relationship between community composition and its ecological functions. In this study, the genomic and transcriptional responses of microbial communities to the availability of alginate, an abundant HMWOM in coastal ecosystem, were elucidated. The main degraders transcribing alginate lyase (Aly) genes came from genera Alteromonas, Psychrosphaera and Colwellia. Meanwhile, some strains, mainly from the Rhodobacteraceae family, did not transcribe Aly gene but could utilize monosaccharides to grow. The co-culture experiment showed that the activity of Aly-producing strain could promote the growth of Aly-non-producing strain when alginate was the sole carbon source. Interestingly, this interaction did not reduce the alginate degradation rate, possibly due to the easily degradable nature of alginate. This study can improve our understanding of the relationship between microbial community activity and alginate metabolism function as well as further manipulation of microbial community structure for alginate processing.

Functional redundancy of seasonal vitamin B12 biosynthesis pathways in coastal marine microbial communities

Vitamin B12 (cobalamin) is a major cofactor required by most marine microbes, but only produced by a few prokaryotes in the ocean, which is globally B12 -depleted. Despite the ecological importance of B12 , the seasonality of B12 metabolisms and the organisms involved in its synthesis in the ocean remain poorly known. Here we use metagenomics to assess the monthly dynamics of B12 -related pathways and the functional diversity of associated microbial communities in the coastal NW Mediterranean Sea over 7 years. We show that genes related to potential B12 metabolisms were characterized by an annual succession of different organisms carrying distinct production pathways. During the most productive winter months, archaea (Nitrosopumilus and Nitrosopelagicus) were the main contributors to B12 synthesis potential through the anaerobic pathway (cbi genes). In turn, Alphaproteobacteria (HIMB11, UBA8309, Puniceispirillum) contributed to B12 synthesis potential in spring and summer through the aerobic pathway (cob genes). Cyanobacteria could produce pseudo-cobalamin from spring to autumn. Finally, we show that during years with environmental perturbations, the organisms usually carrying B12 synthesis genes were replaced by others having the same gene, thus maintaining the potential for B12 production. Such ecological insurance could contribute to the long-term functional resilience of marine microbial communities exposed to contrasting inter-annual environmental conditions.

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.

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

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

Novel RNA viruses from the Atlantic Ocean: Ecogenomics, biogeography, and total virioplankton mass contribution from surface to the deep ocean

Marine viruses play a major role in the energy and nutrient cycle and affect the evolution of their hosts. Despite their importance, there is still little knowledge about RNA viruses. Here, we have explored the Atlantic Ocean, from surface to deep (4.296 m), and used viromics and quantitative methods to unveil the genomics, biogeography, and the mass contribution of RNA viruses to the total viroplankton. A total of 2481 putative RNA viral contigs (>500 bp) and 107 larger bona fide RNA viral genomes (>2.5 kb) were identified; 88 of them representing novel viruses belonging mostly to two clades: Yangshan assemblage (sister clade to the class Alsuviricetes) and Nodaviridae. These viruses were highly endemic and locally abundant, with little or no presence in other oceans since only ≈15% of them were found in at least one of the Tara sampling metatranscriptomes. Quantitative data indicated that the abundance of RNA viruses in the surface and deep chlorophyll maximum zone was within ≈106  VLP/mL representing a potential contribution of 5.2%-24.4% to the total viroplankton community (DNA and RNA viruses), with DNA viruses being the predominant members (≈107  VLP/mL). However, for the deep sample, the observed trend was the opposite, although as further discussed, several biases should be considered. Together these results contribute to our understanding of the diversity, abundance, and distribution of RNA viruses in the oceans and provide a basis for further investigation into their ecological roles and biogeography.

Metagenome-assembled genomes reveal greatly expanded taxonomic and functional diversification of the abundant marine Roseobacter RCA cluster

BACKGROUND

The RCA (Roseobacter clade affiliated) cluster belongs to the family Roseobacteracea and represents a major Roseobacter lineage in temperate to polar oceans. Despite its prevalence and abundance, only a few genomes and one described species, Planktomarina temperata, exist. To gain more insights into our limited understanding of this cluster and its taxonomic and functional diversity and biogeography, we screened metagenomic datasets from the global oceans and reconstructed metagenome-assembled genomes (MAG) affiliated to this cluster.

RESULTS

The total of 82 MAGs, plus five genomes of isolates, reveal an unexpected diversity and novel insights into the genomic features, the functional diversity, and greatly refined biogeographic patterns of the RCA cluster. This cluster is subdivided into three genera: Planktomarina, Pseudoplanktomarina, and the most deeply branching Candidatus Paraplanktomarina. Six of the eight Planktomarina species have larger genome sizes (2.44-3.12 Mbp) and higher G + C contents (46.36-53.70%) than the four Pseudoplanktomarina species (2.26-2.72 Mbp, 42.22-43.72 G + C%). Cand. Paraplanktomarina is represented only by one species with a genome size of 2.40 Mbp and a G + C content of 45.85%. Three novel species of the genera Planktomarina and Pseudoplanktomarina are validly described according to the SeqCode nomenclature for prokaryotic genomes. Aerobic anoxygenic photosynthesis (AAP) is encoded in three Planktomarina species. Unexpectedly, proteorhodopsin (PR) is encoded in the other Planktomarina and all Pseudoplanktomarina species, suggesting that this light-driven proton pump is the most important mode of acquiring complementary energy of the RCA cluster. The Pseudoplanktomarina species exhibit differences in functional traits compared to Planktomarina species and adaptations to more resource-limited conditions. An assessment of the global biogeography of the different species greatly expands the range of occurrence and shows that the different species exhibit distinct biogeographic patterns. They partially reflect the genomic features of the species.

CONCLUSIONS

Our detailed MAG-based analyses shed new light on the diversification, environmental adaptation, and global biogeography of a major lineage of pelagic bacteria. The taxonomic delineation and validation by the SeqCode nomenclature of prominent genera and species of the RCA cluster may be a promising way for a refined taxonomic identification of major prokaryotic lineages and sublineages in marine and other prokaryotic communities assessed by metagenomics approaches. Video Abstract.

Marine nitrogen-fixers in the Canadian Arctic Gateway are dominated by biogeographically distinct noncyanobacterial communities

We describe diazotrophs present during a 2015 GEOTRACES expedition through the Canadian Arctic Gateway (CAG) using nifH metabarcoding. In the less studied Labrador Sea, Bradyrhizobium sp. and Vitreoscilla sp. nifH variants were dominant, while in Baffin Bay, a Stutzerimonas stutzeri variant was dominant. In comparison, the Canadian Arctic Archipelago (CAA) was characterized by a broader set of dominant variants belonging to Desulfobulbaceae, Desulfuromonadales, Arcobacter sp., Vibrio spp., and Sulfuriferula sp. Although dominant diazotrophs fell within known nifH clusters I and III, only a few of these variants were frequently recovered in a 5-year weekly nifH times series in the coastal NW Atlantic presented herein, notably S. stutzeri and variants belonging to Desulfobacterales and Desulfuromonadales. In addition, the majority of dominant Arctic nifH variants shared low similarity (< 92% nucleotide identities) to sequences in a global noncyanobacterial diazotroph catalog recently compiled by others. We further detected UCYN-A throughout the CAG at low-levels using quantitative-PCR assays. Temperature, depth, salinity, oxygen, and nitrate were most strongly correlated to the Arctic diazotroph diversity observed, and we found a stark division between diazotroph communities of the Labrador Sea versus Baffin Bay and the CAA, hence establishing that a previously unknown biogeographic community division can occur for diazotrophs in the CAG.

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

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

Species' functional traits and interactions drive nitrate-mediated sulfur-oxidizing community structure and functioning

IMPORTANCE

Understanding the processes and mechanisms governing microbial community assembly and their linkages to ecosystem functioning has long been a core issue in microbial ecology. An in-depth insight still requires combining with analyses of species' functional traits and microbial interactions. Our study showed how species' functional traits and interactions determined microbial community structure and functions by a well-controlled laboratory experiment with nitrate-mediated sulfur oxidation systems using high-throughput sequencing and culture-dependent technologies. The results provided solid evidences that species' functional traits and interactions were the intrinsic factors determining community structure and function. More importantly, our study established quantitative links between community structure and function based on species' functional traits and interactions, which would have important implications for the design and synthesis of microbiomes with expected functions.

Hydrodynamic regimes modulate nitrogen fixation and the mode of diazotrophy in Lake Tanganyika

The factors that govern the geographical distribution of nitrogen fixation are fundamental to providing accurate nitrogen budgets in aquatic environments. Model-based insights have demonstrated that regional hydrodynamics strongly impact nitrogen fixation. However, the mechanisms establishing this physical-biological coupling have yet to be constrained in field surveys. Here, we examine the distribution of nitrogen fixation in Lake Tanganyika - a model system with well-defined hydrodynamic regimes. We report that nitrogen fixation is five times higher under stratified than under upwelling conditions. Under stratified conditions, the limited resupply of inorganic nitrogen to surface waters, combined with greater light penetration, promotes the activity of bloom-forming photoautotrophic diazotrophs. In contrast, upwelling conditions support predominantly heterotrophic diazotrophs, which are uniquely suited to chemotactic foraging in a more dynamic nutrient landscape. We suggest that these hydrodynamic regimes (stratification versus mixing) play an important role in governing both the rates and the mode of nitrogen fixation.

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.

Environmental transcriptomics under heat stress: Can environmental RNA reveal changes in gene expression of aquatic organisms?

To safeguard biodiversity in a changing climate, taxonomic information about species turnover and insights into the health of organisms are required. Environmental DNA approaches are increasingly used for species identification, but cannot provide functional insights. Transcriptomic methods reveal the physiological states of macroorganisms, but are currently species-specific and require tissue sampling or animal sacrifice, making community-wide assessments challenging. Here, we test whether broad functional information (expression level of the transcribed genes) can be harnessed from environmental RNA (eRNA), which includes extra-organismal RNA from macroorganisms along with whole microorganisms. We exposed Daphnia pulex as well as phytoplankton prey and microorganism colonizers to control (20°C) and heat stress (28°C) conditions for 7 days. We sequenced eRNA from tank water (after complete removal of Daphnia) as well as RNA from Daphnia tissue, enabling comparisons of extra-organismal and organismal RNA-based gene expression profiles. Both RNA types detected similar heat stress responses of Daphnia. Using eRNA, we identified 32 Daphnia genes to be differentially expressed following heat stress. Of these, 17 were also differentially expressed and exhibited similar levels of relative expression in organismal RNA. In addition to the extra-organismal Daphnia response, eRNA detected community-wide heat stress responses consisting of distinct functional profiles and 121 differentially expressed genes across eight taxa. Our study demonstrates that environmental transcriptomics based on extra-organismal eRNA can noninvasively reveal gene expression responses of macroorganisms following environmental changes, with broad potential implications for the biomonitoring of health across the trophic chain.

Ecological mechanisms and current systems shape the modular structure of the global oceans' prokaryotic seascape

Major biogeographic features of the microbial seascape in the oceans have been established and their underlying ecological mechanisms in the (sub)tropical oceans and the Pacific Ocean identified. However, we still lack a unifying understanding of how prokaryotic communities and biogeographic patterns are affected by large-scale current systems in distinct ocean basins and how they are globally shaped in line with ecological mechanisms. Here we show that prokaryotic communities in the epipelagic Pacific and Atlantic Ocean, in the southern Indian Ocean, and the Mediterranean Sea are composed of modules of co-occurring taxa with similar environmental preferences. The relative partitioning of these modules varies along latitudinal and longitudinal gradients and are related to different hydrographic and biotic conditions. Homogeneous selection and dispersal limitation were identified as the major ecological mechanisms shaping these communities and their free-living (FL) and particle-associated (PA) fractions. Large-scale current systems govern the dispersal of prokaryotic modules leading to the highest diversity near subtropical fronts.

Exploring environmental intra-species diversity through non-redundant pangenome assemblies

At the genome level, microorganisms are highly adaptable both in terms of allele and gene composition. Such heritable traits emerge in response to different environmental niches and can have a profound influence on microbial community dynamics. As a consequence, any individual genome or population will contain merely a fraction of the total genetic diversity of any operationally defined "species", whose ecological potential can thus be only fully understood by studying all of their genomes and the genes therein. This concept, known as the pangenome, is valuable for studying microbial ecology and evolution, as it partitions genomes into core (present in all the genomes from a species, and responsible for housekeeping and species-level niche adaptation among others) and accessory regions (present only in some, and responsible for intra-species differentiation). Here we present SuperPang, an algorithm producing pangenome assemblies from a set of input genomes of varying quality, including metagenome-assembled genomes (MAGs). SuperPang runs in linear time and its results are complete, non-redundant, preserve gene ordering and contain both coding and non-coding regions. Our approach provides a modular view of the pangenome, identifying operons and genomic islands, and allowing to track their prevalence in different populations. We illustrate this by analysing intra-species diversity in Polynucleobacter, a bacterial genus ubiquitous in freshwater ecosystems, characterized by their streamlined genomes and their ecological versatility. We show how SuperPang facilitates the simultaneous analysis of allelic and gene content variation under different environmental pressures, allowing us to study the drivers of microbial diversification at unprecedented resolution.

A global atlas of marine antibiotic resistance genes and their expression

Oceans serve as global reservoirs of antibiotic-resistant bacteria and antibiotic resistance genes (ARGs). However, little is known about the traits and expression of ARGs in response to environmental factors. We analyzed 347 metagenomes and 182 metatranscriptomes to determine the distribution, hosts, and expression of ARGs in oceans. Our study found that the diversity and abundance of ARGs varied with latitude and depth. The core marine resistome mainly conferred glycopeptide and multidrug resistance. The hosts of this resistome were mainly limited to the core marine microbiome, with phylogenetic barriers to the horizontal transfer of ARGs, transfers being more frequent within species than between species. Sixty-five percent of the marine ARGs identified were expressed. More than 90% of high-risk ARGs were more likely to be expressed. Anthropogenic activity might affect the expression of ARGs by altering nitrate and phosphate concentrations and ocean temperature. Machine-learning models predict >97% of marine ARGs will change expression by 2100. High-risk ARGs will shift to low latitudes and regions with high anthropogenic activity, such as the Pacific and Atlantic Oceans. Certain ARGs serve a dual role in antibiotic resistance and potentially participate in element cycling, along with other unknown functions. Determining whether changes in ARG expression are beneficial to ecosystems and human health is challenging without comprehensive understanding of their functions. Our study identified a core resistome in the oceans and quantified the expression of ARGs for the development of future control strategies under global change.

Clostridium scindens secretome suppresses virulence gene expression of Clostridioides difficile in a bile acid-independent manner

Clostridioides difficile infection (CDI) is a major health concern and one of the leading causes of hospital-acquired diarrhea in many countries. C. difficile infection is challenging to treat as C. difficile is resistant to multiple antibiotics. Alternative solutions are needed as conventional treatment with broad-spectrum antibiotics often leads to recurrent CDI. Recent studies have shown that specific microbiota-based therapeutics such as bile acids (BAs) are promising approaches to treat CDI. Clostridium scindens encodes the bile acid-induced (bai) operon that carries out 7-alpha-dehydroxylation of liver-derived primary BAs to secondary BAs. This biotransformation is thought to increase the antibacterial effects of BAs on C. difficile. Here, we used an automated multistage fermentor to study the antibacterial actions of C. scindens and BAs on C. difficile in the presence/absence of a gut microbial community derived from healthy human donor fecal microbiota. We observed that C. scindens inhibited C. difficile growth when the medium was supplemented with primary BAs. Transcriptomic analysis indicated upregulation of C. scindens bai operon and suppressed expression of C. difficile exotoxins that mediate CDI. We also observed BA-independent antibacterial activity of the secretome from C. scindens cultured overnight in a medium without supplementary primary BAs, which suppressed growth and exotoxin expression in C. difficile mono-culture. Further investigation of the molecular basis of our observation could lead to a more specific treatment for CDI than current approaches. IMPORTANCE There is an urgent need for new approaches to replace the available treatment options against Clostridioides difficile infection (CDI). Our novel work reports a bile acid-independent reduction of C. difficile growth and virulence gene expression by the secretome of Clostridium scindens. This potential treatment combined with other antimicrobial strategies could facilitate the development of alternative therapies in anticipation of CDI and in turn reduce the risk of antimicrobial resistance.

Consistency across multi-omics layers in a drug-perturbed gut microbial community

Multi-omics analyses are used in microbiome studies to understand molecular changes in microbial communities exposed to different conditions. However, it is not always clear how much each omics data type contributes to our understanding and whether they are concordant with each other. Here, we map the molecular response of a synthetic community of 32 human gut bacteria to three non-antibiotic drugs by using five omics layers (16S rRNA gene profiling, metagenomics, metatranscriptomics, metaproteomics and metabolomics). We find that all the omics methods with species resolution are highly consistent in estimating relative species abundances. Furthermore, different omics methods complement each other for capturing functional changes. For example, while nearly all the omics data types captured that the antipsychotic drug chlorpromazine selectively inhibits Bacteroidota representatives in the community, the metatranscriptome and metaproteome suggested that the drug induces stress responses related to protein quality control. Metabolomics revealed a decrease in oligosaccharide uptake, likely caused by Bacteroidota depletion. Our study highlights how multi-omics datasets can be utilized to reveal complex molecular responses to external perturbations in microbial communities.

A comprehensive genomic catalog from global cold seeps

Cold seeps harbor abundant and diverse microbes with tremendous potential for biological applications and that have a significant influence on biogeochemical cycles. Although recent metagenomic studies have expanded our understanding of the community and function of seep microorganisms, knowledge of the diversity and genetic repertoire of global seep microbes is lacking. Here, we collected a compilation of 165 metagenomic datasets from 16 cold seep sites across the globe to construct a comprehensive gene and genome catalog. The non-redundant gene catalog comprised 147 million genes, and 36% of them could not be assigned to a function with the currently available databases. A total of 3,164 species-level representative metagenome-assembled genomes (MAGs) were obtained, most of which (94%) belonged to novel species. Of them, 81 ANME species were identified that cover all subclades except ANME-2d, and 23 syntrophic SRB species spanned the Seep-SRB1a, Seep-SRB1g, and Seep-SRB2 clades. The non-redundant gene and MAG catalog is a valuable resource that will aid in deepening our understanding of the functions of cold seep microbiomes.

Human-induced salinity changes impact marine organisms and ecosystems

Climate change is fundamentally altering marine and coastal ecosystems on a global scale. While the effects of ocean warming and acidification on ecology and ecosystem functions and services are being comprehensively researched, less attention is directed toward understanding the impacts of human-driven ocean salinity changes. The global water cycle operates through water fluxes expressed as precipitation, evaporation, and freshwater runoff from land. Changes to these in turn modulate ocean salinity and shape the marine and coastal environment by affecting ocean currents, stratification, oxygen saturation, and sea level rise. Besides the direct impact on ocean physical processes, salinity changes impact ocean biological functions with the ecophysiological consequences are being poorly understood. This is surprising as salinity changes may impact diversity, ecosystem and habitat structure loss, and community shifts including trophic cascades. Climate model future projections (of end of the century salinity changes) indicate magnitudes that lead to modification of open ocean plankton community structure and habitat suitability of coral reef communities. Such salinity changes are also capable of affecting the diversity and metabolic capacity of coastal microorganisms and impairing the photosynthetic capacity of (coastal and open ocean) phytoplankton, macroalgae, and seagrass, with downstream ramifications on global biogeochemical cycling. The scarcity of comprehensive salinity data in dynamic coastal regions warrants additional attention. Such datasets are crucial to quantify salinity-based ecosystem function relationships and project such changes that ultimately link into carbon sequestration and freshwater as well as food availability to human populations around the globe. It is critical to integrate vigorous high-quality salinity data with interacting key environmental parameters (e.g., temperature, nutrients, oxygen) for a comprehensive understanding of anthropogenically induced marine changes and its impact on human health and the global economy.

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.

Glacial meltwater and seasonality influence community composition of diazotrophs in Arctic coastal and open waters

The Arctic Ocean is particularly affected by climate change with unknown consequences for primary productivity. Diazotrophs-prokaryotes capable of converting atmospheric nitrogen to ammonia-have been detected in the often nitrogen-limited Arctic Ocean but distribution and community composition dynamics are largely unknown. We performed amplicon sequencing of the diazotroph marker gene nifH from glacial rivers, coastal, and open ocean regions and identified regionally distinct Arctic communities. Proteobacterial diazotrophs dominated all seasons, epi- to mesopelagic depths and rivers to open waters and, surprisingly, Cyanobacteria were only sporadically identified in coastal and freshwaters. The upstream environment of glacial rivers influenced diazotroph diversity, and in marine samples putative anaerobic sulphate-reducers showed seasonal succession with highest prevalence in summer to polar night. Betaproteobacteria (Burkholderiales, Nitrosomonadales, and Rhodocyclales) were typically found in rivers and freshwater-influenced waters, and Delta- (Desulfuromonadales, Desulfobacterales, and Desulfovibrionales) and Gammaproteobacteria in marine waters. The identified community composition dynamics, likely driven by runoff, inorganic nutrients, particulate organic carbon, and seasonality, imply diazotrophy a phenotype of ecological relevance with expected responsiveness to ongoing climate change. Our study largely expands baseline knowledge of Arctic diazotrophs-a prerequisite to understand underpinning of nitrogen fixation-and supports nitrogen fixation as a contributor of new nitrogen in the rapidly changing Arctic Ocean.

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.