Functional annotation of orthologs in metagenomes: a case study of genes for the transformation of oceanic dimethylsulfoniopropionate

Genome-resolved metagenomics of Venice Lagoon surface sediment bacteria reveals high biosynthetic potential and metabolic plasticity as successful strategies in an impacted environment

Bacteria living in sediments play essential roles in marine ecosystems and deeper insights into the ecology and biogeochemistry of these largely unexplored organisms can be obtained from 'omics' approaches. Here, we characterized metagenome-assembled-genomes (MAGs) from the surface sediment microbes of the Venice Lagoon (northern Adriatic Sea) in distinct sub-basins exposed to various natural and anthropogenic pressures. MAGs were explored for biodiversity, major marine metabolic processes, anthropogenic activity-related functions, adaptations at the microscale, and biosynthetic gene clusters. Starting from 126 MAGs, a non-redundant dataset of 58 was compiled, the majority of which (35) belonged to (Alpha- and Gamma-) Proteobacteria. Within the broad microbial metabolic repertoire (including C, N, and S metabolisms) the potential to live without oxygen emerged as one of the most important features. Mixotrophy was also found as a successful lifestyle. Cluster analysis showed that different MAGs encoded the same metabolic patterns (e.g., C fixation, sulfate oxidation) thus suggesting metabolic redundancy. Antibiotic and toxic compounds resistance genes were coupled, a condition that could promote the spreading of these genetic traits. MAGs showed a high biosynthetic potential related to antimicrobial and biotechnological classes and to organism defense and interactions as well as adaptive strategies for micronutrient uptake and cellular detoxification. Our results highlighted that bacteria living in an impacted environment, such as the surface sediments of the Venice Lagoon, may benefit from metabolic plasticity as well as from the synthesis of a wide array of secondary metabolites, promoting ecosystem resilience and stability toward environmental pressures.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00192-z.

AnnoPRO: a strategy for protein function annotation based on multi-scale protein representation and a hybrid deep learning of dual-path encoding

Protein function annotation has been one of the longstanding issues in biological sciences, and various computational methods have been developed. However, the existing methods suffer from a serious long-tail problem, with a large number of GO families containing few annotated proteins. Herein, an innovative strategy named AnnoPRO was therefore constructed by enabling sequence-based multi-scale protein representation, dual-path protein encoding using pre-training, and function annotation by long short-term memory-based decoding. A variety of case studies based on different benchmarks were conducted, which confirmed the superior performance of AnnoPRO among available methods. Source code and models have been made freely available at: https://github.com/idrblab/AnnoPRO and https://zenodo.org/records/10012272.

Complete genome sequence of one novel marine Pseudomonas sp. BSw22131 growing with dimethylsulfoniopropionate (DMSP) as the sole carbon source

Members of the genus Pseudomonas have been frequently isolated from the marine environment, indicating their ecological role in native habitats. One bacterial strain, Pseudomonas sp. BSw22131, was isolated from seawater in Kongsfjorden, Svalbard. The bacterium can grow with algae-derived dimethylsulfoniopropionate (DMSP) as the sole carbon source. Here, we sequenced the complete genome of strain BSw22131, which contained a single circular chromosome of 5,739,290 (G + C content of 58.23 mol%) without any plasmids. A total of 5362 protein-coding genes, 65 tRNA genes, and 16 rRNA genes were obtained. Genome sequence analysis revealed that strain BSw22131 was not only a potential novel species of the genus Pseudomonas but also different from Pseudomonas sp. DMSP-1 that was isolated from the same habitat and also utilized DMSP as the sole carbon source for growth. The results can be helpful for understanding the catabolism of the genus Pseudomonas in sulfur cycling in the Arctic fjord ecosystem.

Oceanospirillales containing the DMSP lyase DddD are key utilisers of carbon from DMSP in coastal seawater

BACKGROUND

Ubiquitous and diverse marine microorganisms utilise the abundant organosulfur molecule dimethylsulfoniopropionate (DMSP), the main precursor of the climate-active gas dimethylsulfide (DMS), as a source of carbon, sulfur and/or signalling molecules. However, it is currently difficult to discern which microbes actively catabolise DMSP in the environment, why they do so and the pathways used.

RESULTS

Here, a novel DNA-stable isotope probing (SIP) approach, where only the propionate and not the DMS moiety of DMSP was 13C-labelled, was strategically applied to identify key microorganisms actively using DMSP and also likely DMS as a carbon source, and their catabolic enzymes, in North Sea water. Metagenomic analysis of natural seawater suggested that Rhodobacterales (Roseobacter group) and SAR11 bacteria were the major microorganisms degrading DMSP via demethylation and, to a lesser extent, DddP-driven DMSP lysis pathways. However, neither Rhodobacterales and SAR11 bacteria nor their DMSP catabolic genes were prominently labelled in DNA-SIP experiments, suggesting they use DMSP as a sulfur source and/or in signalling pathways, and not primarily for carbon requirements. Instead, DNA-SIP identified gammaproteobacterial Oceanospirillales, e.g. Amphritea, and their DMSP lyase DddD as the dominant microorganisms/enzymes using DMSP as a carbon source. Supporting this, most gammaproteobacterial (with DddD) but few alphaproteobacterial seawater isolates grew on DMSP as sole carbon source and produced DMS. Furthermore, our DNA-SIP strategy also identified Methylophaga and other Piscirickettsiaceae as key bacteria likely using the DMS, generated from DMSP lysis, as a carbon source.

CONCLUSIONS

This is the first study to use DNA-SIP with 13C-labelled DMSP and, in a novel way, it identifies the dominant microbes utilising DMSP and DMS as carbon sources. It highlights that whilst metagenomic analyses of marine environments can predict microorganisms/genes that degrade DMSP and DMS based on their abundance, it cannot disentangle those using these important organosulfur compounds for their carbon requirements. Note, the most abundant DMSP degraders, e.g. Rhodobacterales with DmdA, are not always the key microorganisms using DMSP for carbon and releasing DMS, which in this coastal system were Oceanospirillales containing DddD. Video abstract.

Recent insights into oceanic dimethylsulfoniopropionate biosynthesis and catabolism

Dimethylsulfoniopropionate (DMSP), a globally important organosulfur compound is produced in prodigious amounts (2.0 Pg sulfur) annually in the marine environment by phytoplankton, macroalgae, heterotrophic bacteria, some corals and certain higher plants. It is an important marine osmolyte and a major precursor molecule for the production of climate-active volatile gas dimethyl sulfide (DMS). DMSP synthesis take place via three pathways: a transamination 'pathway-' in some marine bacteria and algae, a Met-methylation 'pathway-' in angiosperms and bacteria and a decarboxylation 'pathway-' in the dinoflagellate, Crypthecodinium. The enzymes DSYB and TpMMT are involved in the DMSP biosynthesis in eukaryotes while marine heterotrophic bacteria engage key enzymes such as DsyB and MmtN. Several marine bacterial communities import DMSP and degrade it via cleavage or demethylation pathways or oxidation pathway, thereby generating DMS, methanethiol, and dimethylsulfoxonium propionate, respectively. DMSP is cleaved through diverse DMSP lyase enzymes in bacteria and via Alma1 enzyme in phytoplankton. The demethylation pathway involves four different enzymes, namely DmdA, DmdB, DmdC and DmdD/AcuH. However, enzymes involved in the oxidation pathway have not been yet identified. We reviewed the recent advances on the synthesis and catabolism of DMSP and enzymes that are involved in these processes.

Phylogenomics of SAR116 Clade Reveals Two Subclades with Different Evolutionary Trajectories and an Important Role in the Ocean Sulfur Cycle

The SAR116 clade within the class Alphaproteobacteria represents one of the most abundant groups of heterotrophic bacteria inhabiting the surface of the ocean. The small number of cultured representatives of SAR116 (only two to date) is a major bottleneck that has prevented an in-depth study at the genomic level to understand the relationship between genome diversity and its role in the marine environment. In this study, we use all publicly available genomes to provide a genomic overview of the phylogeny, metabolism, and biogeography within the SAR116 clade. This increased genomic diversity has led to the discovery of two subclades that, despite coexisting in the same environment, display different properties in their genomic makeup. One represents a novel subclade for which no pure cultures have been isolated and is composed mainly of single-amplified genomes (SAGs). Genomes within this subclade showed convergent evolutionary trajectories with more streamlined features, such as low GC content (ca. 30%), short intergenic spacers (<22 bp), and strong purifying selection (low ratio of nonsynonymous to synonymous polymorphisms [dN/dS]). Besides, they were more abundant in metagenomic databases recruiting at the deep chlorophyll maximum. Less abundant and restricted to the upper photic layers of the global ocean, the other subclade of SAR116, enriched in metagenome-assembled genomes (MAGs), included the only two pure cultures. Genomic analysis suggested that both clades have a significant role in the sulfur cycle with differences in the way both clades can metabolize dimethylsulfoniopropionate (DMSP). IMPORTANCE The SAR116 clade of Alphaproteobacteria is a ubiquitous group of heterotrophic bacteria inhabiting the surface of the ocean, but the information about their ecology and population genomic diversity is scarce due to the difficulty of getting pure culture isolates. The combination of single-cell genomics and metagenomics has become an alternative approach to study these kinds of microbes. Our results expand the understanding of the genomic diversity, distribution, and lifestyles within this clade and provide evidence of different evolutionary trajectories in the genomic makeup of the two subclades that could serve to illustrate how evolutionary pressure can drive different adaptations to the same environment. Therefore, the SAR116 clade represents an ideal model organism for the study of the evolutionary streamlining of genomes in microbes that have relatively close relatedness to each other.

Biogeographic traits of dimethyl sulfide and dimethylsulfoniopropionate cycling in polar oceans

BACKGROUND

Dimethyl sulfide (DMS) is the dominant volatile organic sulfur in global oceans. The predominant source of oceanic DMS is the cleavage of dimethylsulfoniopropionate (DMSP), which can be produced by marine bacteria and phytoplankton. Polar oceans, which represent about one fifth of Earth's surface, contribute significantly to the global oceanic DMS sea-air flux. However, a global overview of DMS and DMSP cycling in polar oceans is still lacking and the key genes and the microbial assemblages involved in DMSP/DMS transformation remain to be fully unveiled.

RESULTS

Here, we systematically investigated the biogeographic traits of 16 key microbial enzymes involved in DMS/DMSP cycling in 60 metagenomic samples from polar waters, together with 174 metagenome and 151 metatranscriptomes from non-polar Tara Ocean dataset. Our analyses suggest that intense DMS/DMSP cycling occurs in the polar oceans. DMSP demethylase (DmdA), DMSP lyases (DddD, DddP, and DddK), and trimethylamine monooxygenase (Tmm, which oxidizes DMS to dimethylsulfoxide) were the most prevalent bacterial genes involved in global DMS/DMSP cycling. Alphaproteobacteria (Pelagibacterales) and Gammaproteobacteria appear to play prominent roles in DMS/DMSP cycling in polar oceans. The phenomenon that multiple DMS/DMSP cycling genes co-occurred in the same bacterial genome was also observed in metagenome assembled genomes (MAGs) from polar oceans. The microbial assemblages from the polar oceans were significantly correlated with water depth rather than geographic distance, suggesting the differences of habitats between surface and deep waters rather than dispersal limitation are the key factors shaping microbial assemblages involved in DMS/DMSP cycling in polar oceans.

CONCLUSIONS

Overall, this study provides a global overview of the biogeographic traits of known bacterial genes involved in DMS/DMSP cycling from the Arctic and Antarctic oceans, laying a solid foundation for further studies of DMS/DMSP cycling in polar ocean microbiome at the enzymatic, metabolic, and processual levels. Video Abstract.

Evolutionary history of dimethylsulfoniopropionate (DMSP) demethylation enzyme DmdA in marine bacteria

Dimethylsulfoniopropionate (DMSP), an osmolyte produced by oceanic phytoplankton and bacteria, is primarily degraded by bacteria belonging to the Roseobacter lineage and other marine Alphaproteobacteria via DMSP-dependent demethylase A protein (DmdA). To date, the evolutionary history of DmdA gene family is unclear. Some studies indicate a common ancestry between DmdA and GcvT gene families and a co-evolution between Roseobacter and the DMSP-producing-phytoplankton around 250 million years ago (Mya). In this work, we analyzed the evolution of DmdA under three possible evolutionary scenarios: (1) a recent common ancestor of DmdA and GcvT, (2) a coevolution between Roseobacter and the DMSP-producing-phytoplankton, and (3) an enzymatic adaptation for utilizing DMSP in marine bacteria prior to Roseobacter origin. Our analyses indicate that DmdA is a new gene family originated from GcvT genes by duplication and functional divergence driven by positive selection before a coevolution between Roseobacter and phytoplankton. Our data suggest that Roseobacter acquired dmdA by horizontal gene transfer prior to an environment with higher DMSP. Here, we propose that the ancestor that carried the DMSP demethylation pathway genes evolved in the Archean, and was exposed to a higher concentration of DMSP in a sulfur-rich atmosphere and anoxic ocean, compared to recent Roseobacter eco-orthologs (orthologs performing the same function under different conditions), which should be adapted to lower concentrations of DMSP.

Bacteria are important dimethylsulfoniopropionate producers in marine aphotic and high-pressure environments

Dimethylsulfoniopropionate (DMSP) is an important marine osmolyte. Aphotic environments are only recently being considered as potential contributors to global DMSP production. Here, our Mariana Trench study reveals a typical seawater DMSP/dimethylsulfide (DMS) profile, with highest concentrations in the euphotic zone and decreased but consistent levels below. The genetic potential for bacterial DMSP synthesis via the dsyB gene and its transcription is greater in the deep ocean, and is highest in the sediment.s DMSP catabolic potential is present throughout the trench waters, but is less prominent below 8000 m, perhaps indicating a preference to store DMSP in the deep for stress protection. Deep ocean bacterial isolates show enhanced DMSP production under increased hydrostatic pressure. Furthermore, bacterial dsyB mutants are less tolerant of deep ocean pressures than wild-type strains. Thus, we propose a physiological function for DMSP in hydrostatic pressure protection, and that bacteria are key DMSP producers in deep seawater and sediment.

Dimethylsulfoniopropionate (DMSP) is an osmolyte produced by marine microbes that plays an important role in nutrient cycling and atmospheric chemistry. Here the authors go to the Mariana Trench—the deepest point in the ocean—and find bacteria are key DMSP producers, and that DMSP has a role in protection against high pressure.

Metagenomics-Guided Discovery of Potential Bacterial Metallothionein Genes from the Soil Microbiome That Confer Cu and/or Cd Resistance

Metallothionein (MT) genes are valuable genetic materials for developing metal bioremediation tools. Currently, a limited number of prokaryotic MTs have been experimentally identified, which necessitates the expansion of bacterial MT diversity. In this study, we conducted a metagenomics-guided analysis for the discovery of potential bacterial MT genes from the soil microbiome. More specifically, we combined resistance gene enrichment through diversity loss, metagenomic mining with a dedicated MT database, evolutionary trace analysis, DNA chemical synthesis, and functional genomic validation to identify novel MTs. Results showed that Cu stress induced a compositional change in the soil microbiome, with an enrichment of metal-resistant bacteria in soils with higher Cu concentrations. Shotgun metagenomic sequencing was performed to obtain the gene pool of environmental DNA (eDNA), which was subjected to a local BLAST search against an MT database for detecting putative MT genes. Evolutional trace analysis led to the identification of 27 potential MTs with conserved cysteine/histidine motifs different from those of known prokaryotic MTs. Following chemical synthesis of these 27 potential MT genes and heterologous expression in Escherichia coli, six of them were found to improve the hosts' growth substantially and enhanced the hosts' sorption of Cu, Cd, and Zn, among which MT5 led to a 13.7-fold increase in Cd accumulation. Furthermore, four of them restored Cu and/or Cd resistance in two metal-sensitive E. coli strains.IMPORTANCE The metagenomics-guided procedure developed here bypasses the difficulties encountered in classic PCR-based approaches and led to the discovery of novel MT genes, which may be useful in developing bioremediation tools. The procedure used here expands our knowledge on the diversity of bacterial MTs in the environment and may also be applicable to identify other functional genes from eDNA.

Dimethylsulfoniopropionate Sulfur and Methyl Carbon Assimilation in Ruegeria Species

Dimethylsulfoniopropionate (DMSP) is abundant in marine environments and an important source of reduced carbon and sulfur for marine bacteria. While both Ruegeria pomeroyi and Ruegeria lacuscaerulensis possessed genes encoding the DMSP demethylation and cleavage pathways, their responses to DMSP differed. A glucose-fed, chemostat culture of R. pomeroyi consumed 99% of the DMSP even when fed a high concentration of 5 mM. At the same time, cultures released 19% and 7.1% of the DMSP as dimethylsulfide (DMS) and methanethiol, respectively. Under the same conditions, R. lacuscaerulensis consumed only 28% of the DMSP and formed one-third of the amount of gases. To examine the pathways of sulfur and methyl C assimilation, glucose-fed chemostats of both species were fed 100 μM mixtures of unlabeled and doubly labeled [dimethyl-¹³C, ³⁴S]DMSP. Both species derived nearly all of their sulfur from DMSP despite high sulfate availability. In addition, only 33% and 50% of the methionine was biosynthesized from the direct capture of methanethiol in R. pomeroyi and R. lacuscaerulensis, respectively. The remaining methionine was biosynthesized by the random assembly of free sulfide and methyl-tetrahydrofolate derived from DMSP. Thus, although the two species possessed similar genes encoding DMSP metabolism, their growth responses were very different.IMPORTANCE Dimethylsulfoniopropionate (DMSP) is abundant in marine environments and an important source of reduced carbon and sulfur for marine bacteria. DMSP is the precursor for the majority of atmospheric dimethylsulfide (DMS), a climatically active gas that connects the marine and terrestrial sulfur cycles. Although research into the assimilation of DMSP has been conducted for over 20 years, the fate of DMSP in microbial biomass is not well understood. In particular, the biosynthesis of methionine from DMSP has been a focal point, and it has been widely believed that most methionine was synthesized via the direct capture of methanethiol. Using an isotopic labeling strategy, we have demonstrated that the direct capture of methanethiol is not the primary pathway used for methionine biosynthesis in two Ruegeria species, a genus comprised primarily of globally abundant marine bacteria. Furthermore, although the catabolism of DMSP by these species varied greatly, the anabolic pathways were highly conserved.

Metagenomic Insights Into the Cycling of Dimethylsulfoniopropionate and Related Molecules in the Eastern China Marginal Seas

The microbial cycling of dimethylsulfoniopropionate (DMSP) and its gaseous catabolites dimethylsulfide (DMS) and methanethiol (MeSH) are important processes in the global sulfur cycle, marine microbial food webs, signaling pathways, atmospheric chemistry, and potentially climate regulation. Many functional genes have been identified and used to study the genetic potential of microbes to produce and catabolize these organosulfur compounds in different marine environments. Here, we sampled seawater, marine sediment and hydrothermal sediment, and polymetallic sulfide in the eastern Chinese marginal seas and analyzed their microbial communities for the genetic potential to cycle DMSP, DMS, and MeSH using metagenomics. DMSP was abundant in all sediment samples, but was fivefold less prominent in those from hydrothermal samples. Indeed, Yellow Sea (YS) sediment samples had DMSP concentrations two orders of magnitude higher than in surface water samples. Bacterial genetic potential to synthesize DMSP (mainly in Rhodobacteraceae bacteria) was far higher than for phytoplankton in all samples, but particularly in the sediment where no algal DMSP synthesis genes were detected. Thus, we propose bacteria as important DMSP producers in these marine sediments. DMSP catabolic pathways mediated by the DMSP lyase DddP (prominent in Pseudomonas and Mesorhizobium bacteria) and DMSP demethylase DmdA enzymes (prominent in Rhodobacteraceae bacteria) and MddA-mediated MeSH S-methylation were very abundant in Bohai Sea and Yellow Sea sediments (BYSS) samples. In contrast, the genetic potential for DMSP degradation was very low in the hydrothermal sediment samples-dddP was the only catabolic gene detected and in only one sample. However, the potential for DMS production from MeSH (mddA) and DMS oxidation (dmoA and ddhA) was relatively abundant. This metagenomics study does not provide conclusive evidence for DMSP cycling; however, it does highlight the potential importance of bacteria in the synthesis and catabolism of DMSP and related compounds in diverse sediment environments.

A subfamily roadmap of the evolutionarily diverse glycoside hydrolase family 16 (GH16)

Glycoside hydrolase family (GH) 16 comprises a large and taxonomically diverse family of glycosidases and transglycosidases that adopt a common β-jelly-roll fold and are active on a range of terrestrial and marine polysaccharides. Presently, broadly insightful sequence–function correlations in GH16 are hindered by a lack of a systematic subfamily structure. To fill this gap, we have used a highly scalable protein sequence similarity network analysis to delineate nearly 23,000 GH16 sequences into 23 robust subfamilies, which are strongly supported by hidden Markov model and maximum likelihood molecular phylogenetic analyses. Subsequent evaluation of over 40 experimental three-dimensional structures has highlighted key tertiary structural differences, predominantly manifested in active-site loops, that dictate substrate specificity across the GH16 evolutionary landscape. As for other large GH families (i.e. GH5, GH13, and GH43), this new subfamily classification provides a roadmap for functional glycogenomics that will guide future bioinformatics and experimental structure–function analyses. The GH16 subfamily classification is publicly available in the CAZy database. The sequence similarity network workflow used here, SSNpipe, is freely available from GitHub.