Metagenomic Analysis of Carbohydrate-Active Enzymes and Their Contribution to Marine Sediment Biodiversity

doi:10.21203/rs.3.rs-3129032/v1

Download PDF

Research Article

Metagenomic Analysis of Carbohydrate-Active Enzymes and Their Contribution to Marine Sediment Biodiversity

https://doi.org/10.21203/rs.3.rs-3129032/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 13 Feb, 2024

Read the published version in World Journal of Microbiology and Biotechnology →

You are reading this latest preprint version

Marine sediments constitute the world's most substantial long-term carbon repository. The microorganisms dwelling in these sediments mediate the transformation of fixed oceanic carbon, but their contribution to the carbon cycle is not fully understood. Previous culture-independent investigations into sedimentary microorganisms have underscored the significance of carbohydrates in the carbon cycle. In this study, we employ a metagenomic methodology to investigate the distribution and abundance of carbohydrate-active enzymes (CAZymes) in 37 marine sediments. These sediments exhibit varying oxygen availability and were isolated in diverse regions worldwide. Our comparative analysis is based on the metabolic potential for oxygen utilisation, derived from genes present in both oxic and anoxic environments. We found extracellular CAZyme modules targeting the degradation of plant and algal detritus, necromass, and host glycans were more prevalent across all metagenomic samples. The analysis of these results indicates that the oxic / anoxic conditions not only influence the taxonomic composition of the microbial communities, but also affect the occurrence of CAZyme modules involved in the transformation of necromass, algae and plant detritus. To gain insight into the sediment microbial taxa, we reconstructed metagenomic assembly genomes (MAG) and examined the presence of primary extracellular carbohydrate active enzyme (CAZyme) modules. Our findings reveal that the primary CAZyme modules and the CAZyme gene clusters (CGC) discovered in our metagenomes were prevalent in the Bacteroidia, Gammaproteobacteria and Alphaproteobacteria classes. We compared those MAG to MAG from the same classes found in soil, and we found that they were similar, but the soil MAG contained a more abundant and diverse CAZyme content. Furthermore, the data indicate that abundant classes in our metagenomic samples, namely Alphaproteobacteria, Bacteroidia and Gammaproteobacteria, play a pivotal role in carbohydrate transformation within the initial few metres of the sediments.

Anoxic

bioinformatics

CAZymes

oxic

marine sediments

metagenomics.

The ocean floor is the recipient of all the organic matter coming from the water column and is considered the major carbon repository on the planet. Therefore, microorganisms that live in marine sediments control the storage of massive amounts of carbon (Orcutt et al.2011). These microbes that live on and below the sea floor represent more than 10²⁹ cells, a number roughly equal to the number of microorganisms in seawater and soil (Kallmeyer et al., 2012).

Marine sediments can be classified, depending on the availability of electron acceptors such as oxygen and sulphur, into oxic or anoxic. In the oxic subseafloor, the penetration of O₂ and the resulting limitation of the electron donor result in a unique community structure, compared to anoxic sediments (Orsi, 2018). In regions of the seabed with subseafloor anoxia, oxygen is typically consumed in the upper centimetres of the sediment below (D’Hont et al., 2004; Froelich et al., 1979).

In both types of marine sediments (oxic and anoxic), microbial communities process both organic and inorganic carbon and contribute to the cycling of nutrients such as sulphur, nitrogen, and iron (Parkes et al.,2014). Despite the global importance of these organisms, marine sediments are among the least understood environments. This is in part due to the difficulty of sampling, especially in the deep sea, and to the complexity of their inhabiting communities. However, recent examination of prokaryote genes, transcripts, and metagenomes has highlighted the importance of polysaccharides and their transformations for carbon metabolism in the ocean (Teeling et al., 2012; Teeling et al., 2016). Therefore, a closer examination of the marine polysaccharide cycle and the communities driving their degradation is necessary. Although polysaccharides constitute a large fraction of phytoplankton and macroalgae bodies (Biersmith & Benner, 1998) as well as dissolved and particulate organic matter (DOM and POM, respectively) (Lee et al., 2001), little is known about their biogeochemical processing compared to other major compound classes, such as proteins, lipids, and nucleic acids. Carbohydrate-active enzymes (CAZymes) are proteins with known activities involved in the synthesis and degradation of glycoconjugates, oligo- and polysaccharides. They typically correspond to 1–5% of the genes of a living organism. These enzymes play essential roles in life not only as structure and energy reserve components but also in many intracellular and intercellular recognition events. CAZymes are often involved in immune and host-pathogen interactions and are involved in human and agricultural-related diseases. CAZymes have been classified and annotated in the CAZy database since 1998. This is a specialist database dedicated to the display and analysis of genomic, structural, and biochemical information on carbohydrate-activated enzymes (CAZymes) (Lombard et al., 2013).

Here, we present a comparative study of carbohydrate active enzymes (CAZymes) from sediment metagenomes from different locations in the world to better understand their role in the storage or degradation of carbohydrates and derivatives.

Selection of metagenomic data.

Upon identification of suitable BioProjects, the metagenome shotgun sequences were downloaded from the NCBI database. The raw data recovered from 37 metagenome samples from 12 BioProjects representing marine samples included valuable metadata associated with each sample, including latitude/longitude coordinates, metres below the sea floor, and metres below sea level. This additional information provided an important context for understanding the spatial distribution and environmental characteristics of the marine ecosystems sampled. Sediment samples were taken from all over the world with a depth range of 0-7942 metres below sea level (mbsl) and 0 to 2.23 metres below the sediment floor (mbsf) (Fig. 1). To reduce bias from sequencing, samples not sequenced by Illumina were discarded (Supplementary Table 1).

Quality Control and Pre-processing.

Quality control procedures were executed using widely used tools and software, such as Trimmomatic (Bolger et al., 2014) and FastQC (Andrews, 2010). All reads from samples that did not pass the QC filters (read quality > = Q20) were discarded.

Taxonomic analysis of metagenomic reads.

For taxonomic analysis of the reads, we used the Kraken2 specific database based on k-mer spectra from complete RefSeq genomes and the NCBI nt database (downloaded 7/09/2021) (Wood et al., 2019). The annotation tables were formatted for the R ggplot2 library to generate stacked bar plots at different taxonomic levels. The integrated matrices obtained for the 37 samples were written using R, bash, Perl, and Python and are available at https://github.com/jenniferlu717/KrakenTools.

Metagenomic Sequence Read Assembly and Functional Analysis of the Metagenome.

A de novo assembly was made for each sample using MEGAHIT v1.1.1-2 with the parameter 'metasensitive' recommended for diverse samples (Li et al., 2015). The ORFs of each sample were predicted using the Prodigal-v2.6.3 tool (Hyatt et al., 2010). Samples with less than one million genes were discarded from this study. We used HMMER-3.2.1 (Eddy et al., 2011) (hmmscan cutoffs: E-value < 1e-15, coverage > 0.35) to annotate CAZymes against the HMM database V9 of dbCAN2 (Zhang et al.,2018). The substrate specificities of CAZymes were inferred by manual inspection of CAZy (Lombard et al., 2013; Lombard et al., 2014). Extracellular CAZymes were annotated using SignalP V-5.0 (Almagro et al., 2019). Heme-copper oxygen reductases (HCO) and nitric oxide reductases (NOR) were analysed using Diamond (parameters “ultra-sensitive”) against the HCO database (Sousa et al., 2011). Normalisation of gene counts between CAZymes samples and Heme-copper oxygen reductase genes was carried out using the equation: (Number of genes annotated / Number of total genes in sample) X 10⁶.

Diversity and statistical analysis.

Statistical analyses were performed using R-v. 4.2.3 (R Core Team 2023). Using the Bray-Curtis dissimilarity index to calculate distance matrices relative to the taxa abundance group at the class level of taxonomy and CAZyme composition, a Principal Coordinate Analysis (PCoA) was performed. To look for correlations between the metadata from the samples and the taxonomic diversity, a Mantel test of the abundance matrix vs. metres below sea level (mbsl), latitude and longitude metadata of samples, and metres below sea floor (mbsf) were calculated. Data visualisation of PCOAs was performed using the vegan, pragma (Oksanen et al.,2020) and geosphere (Hijmans et al., 2020) R packages. Linear discriminant analysis Effect Size Analysis (LEfSe) of the taxonomic matrices of archaea and bacteria and were made in the Hutlab’s Galaxy tool (Segata et al., 2011) (LEfSe cutoff: Kruskal-Wallis Alfa value Kruskal-Wallis = 0.05, Alfa value Wilcoxon test = 0.05, LDA score > 3.0). A similar analysis was performed on the extracellular CAZyme matrix with an LDA score > 3.5.

Construction of Metagenomic Assembled Genome (MAG) and Functional Analysis.

MAG were annotated, reconstructed, and refined using the Squeeze-Meta with pipeline v.1.4.0 (Tamames & Puente-Sanchez et al., 2019) (parameters: mode = sequential, assembly = extassembly, doublepass, lowmem). Genomic bins with low completeness (< 75%) and high contamination were removed (> 10%). Bins were refined with the remove_duplicate_markers.pl program of the SqueezeMeta pipeline. The taxonomic classification of these bins was performed by GTDB-Tk v2.1.0 (parameters: classify_wf) against the GTDB database v-207 (Chaumeil et al., 2020). CAZyme modules and Cazyme gene clusters (CGC) were annotated using dbCAN2 (Zhang et al., 2018); (hmmscan cutoffs: E-value < 1e-15, coverage > 0.35, DIAMOND cut-offs: E-value < 1e-102, Hotpep (Frequency > 2.6, Hits > 6), CGCFinder (Distance < = 2, signature genes = CAZyme + TC).

Marker genes (MG) were annotated with FetchMG v-1.2 (Kultima et al., 2012).

Normalisation of CAZyme counts between MAG was carried out using the equation: [(Number of the CAZyme module in the MAG / Number of the CAZyme module in the metagenome sample)/ Median (MGs in metagenome sample)] X 10⁶. The list of MGs can be downloaded from the mOTU website https://motu-tool.org/fetchMG.html.

Soil MAG taxonomically assigned to Alphaproteobacteria, Gammaproteobacteria, and Bacteroidia classes (Nayfach et al., 2021) were randomly selected with the same criteria as our MAG (Completeness > 75% and Contamination < 10%). CAZyme modules and CAZyme gene clusters (CGC) were annotated using dbCAN2 (Zhang et al., 2018); (hmmscan cut-offs: E-value < 1e-15, coverage > 0.35, DIAMOND cut-offs: E-value < 1e-102, Hotpep (Frequency > 2.6, Hits > 6). The integrated matrices were written using R, bash, Perl, and Python and are available at http://github.com/Ales-ibt/Metagenomic-benchmark.

Correlations of metagenome samples based on available metadata.

We analysed 37 metagenomes from all over the world. The physicochemical variables of most of our sediment samples were not available for comparison. However, they all come from shallow sediments at the interface with the water column, for which metadata such as geographic parameters (latitude and longitude, depth in metres below the seal level (mbsl), and depth in metres below the seafloor (mbsl) are known.

Twenty samples were taken below 1000 mbsl and 17 above (Fig. 1 and Supplementary Table 1). This allowed us to test the correlation between environmental variables and the abundance diversity matrix at the class level.

We calculated a Mantel test to test whether the taxonomic community structure was correlated with geographical and spatial parameters (Supplementary Table 2). Our results showed a significant and positive correlation between depth (mbsl) and taxonomic diversity (Bray-Curtis dissimilarity matrix), while geographical distances and sediment depth (mbsf) were not significant (Supplementary Table 2). To further analyse this positive correlation between depth of the water column and taxonomic diversity, we made a linear regression of both dissimilatory matrixes (Supplementary Fig. S1). A low R-squared value (0.0404) suggests that depth below sea level does not explain much of the variation in taxonomic dissimilarity.

For that, we decided to make a metagenomic profile of the samples based on the taxonomical diversity of the community against its metabolic potential.

To this end, we assign categories to our samples (oxic/anoxic environments) based on their gene content of heme-copper oxygen reductases (HCO) and nitric oxide reductases (NOR). HCOs and NORs are enzymes found in the last complexes of many respiratory chains in microorganisms (Sousa et al., 2011). As reference, we used four sediments found in Loki's Castle labelled as anoxic and one from the South Pacific Gyre labelled as oxic (Supplementary Table 1). Everything above the oxic control was considered oxic, and everything below was considered anoxic. Our results show 18 metagenome samples that can be considered oxic and 13 anoxic. Some of the samples assigned to the oxic label were shallow samples (under 1000 mbsl) and, although there is a correlation between a greater depth of the water column and thin sediments, because the amount of organic matter is depleted and oxygen penetration is found throughout the sediment, many of our samples were taken in the first centimetres (from 0 to 2.23 mbsf) where the community utilizes oxygen. (D’Hont et al., 2015). This could be the reason why these samples below 1,000 metres below sea level are above the oxic control. (Fig. 2, Supplementary Table 1).

Once we established the abundance of HCO and NOR as a condition in the samples, a principal coordinate analysis (PCoA) based on the relative taxonomic abundance at the class level (Bray-Curtis dissimilarity matrix) showed a clear separation of the samples labelled oxic and anoxic (62.18% of the variance explained in CoA1 and CoA2) (Fig. 3a, Supplementary Table 3). Samples were clustered into two groups; In the oxic group, samples from the deep Gulf of Mexico (Godoy-Lozano et al., 2018; Zhao et al., 2020) are reported without hydrocarbon or methane seeps (Zhao et al., 2020). The South Pacific Gyre is the only sample with an oxic level and an oligotrophic layer (Tully et al., 2016). Samples from Korea and Antarctica present anthropogenic disturbances; the Korea metagenomes are beach samples, the Davis Station are shallow samples rich in nutrients, and oxygen is consumed in the first centimetres of the sediment (Leeming et al., 2015). In the anoxic group, samples from the Gulf of Mexico (Delaware University), the Basin and Loki’s Castle, the Hydrate Ridge of the Pacific, and the Santa Monica Mounds were clustered together. These have been reported to have seepages of hydrocarbons or related compounds (Zhao et al., 2020), hydrothermal vents with anaerobic metabolism (Jaeschke et al., 2012; Kauffman et al., 2018; Bäckström et al., 2019;), and mud volcanoes (Kauffman et al., 2018; Bäckström et al., 2019) (Fig. 3a).

Once we saw a clear separation between labels, we explored differences in taxonomic composition between the oxic and anoxic samples through a LEfSe analysis (Segata et al., 2011) based on bacteria and archaea abundance matrices at the class and phylum levels, respectively (Fig. 3b, Supplementary Fig. S2 and S3).

The oxic samples showed a significant difference in Alphaproteobacteria. However, anoxic samples had significant differences in several bacterial classes: Epsilonbacteria, Deltaproteobacteria, Bacilli, Clostridia, Fusobacteriia, Dehalococcoidia, Bacteroidia, Sphingobacteriia, Cytophagia, and Thermodesulfobacteria. Among Archaea phyla, Thaumarchaeota are significant indicators of oxic samples, while Candidatus Bathyarchaeota, Euryarchaeota, and Candidatus Lokiarchaeota are indicators of anoxic samples. This is consistent with the literature where it is known that anoxic sediments are enriched with strictly anaerobic groups such as sulphate-reducing bacteria of the Chloroflexota phylum and Deltaproteobacteria and methanogenic archaea, such as Euryarchaeota, while in oxic sediments there is prevalence of the Alphaproteobacteria class in bacteria and Thaumarchaeota phylum in archaea (Biddle et al., 2008; Orsi, 2018; Hoshino et al., 2020). Our results found that the classes Dehalococcoidia and Deltaproteobacteria of the Chloroflexota phylum along with other anaerobic classes such as Clostridia, Thermodesulfobacteria, Fusobacteriia bacteria and Euryarchaeota archaea indicative of an anoxic environment, while the Alphaproteobacteria class of bacteria and Thaumarchaeota archaea (Tully et al., 2016; Hoshino et al., 2020) were indicative of oxic samples (Fig. 3b).

In summary, both groups exhibited significant differences in the classes of bacteria and the archaea diversity that appear to match the anoxic / oxic conditions of the microorganisms reported in marine sediments, as well as the genes reported (Fig. 3b).

CAZyme profile of marine sediments.

We examined the distribution of CAZyme (Carbohydrate-Active Enzymes) content within the metagenomes. To accomplish this, we performed a Principal Coordinate Analysis (PCoA) using normalised counts of all CAZyme modules identified within each metagenome sample. Like our findings on beta diversity, our samples showed a clear separation between oxic and anoxic conditions (59.18% of the variance explained in CoA1 and CoA2) (Fig. 4).

Given the assumption that carbon turnover in marine sediments is carried out by microbial organisms that use secreted enzymes to store carbon over time (Orsi et al., 2018), we decided to search for extracellular CAZymes. We performed a functional annotation of CAZyme modules that had a peptide signal against the CAZyme database (Lombard et al., 2013). We categorized sequences into the six classes of the CAZy database, which are implicated in the creation, breakdown, and identification of carbohydrates. These classes are Glycoside Transferases (GTs), Glycoside Hydrolases (GHs), Carbohydrate Esterases (CEs), Carbohydrate Binding Modules (CBMs), Polysaccharide Lyases (PLs), and Auxiliary Activities (AAs). Eighteen extracellular CAZyme modules were found in more than 1% of all total annotations (accounting for 55.94% of all CAZymes annotated in our metagenome samples). Of these modules, GH109, GH23, and CE1 were the most abundant (Fig. 5a). Their abundance was particularly high in the following metagenomes: Guaymas Basin (GBGOC), Davis Station from Antarctica (DSANT), Korean beaches (KOR), South Pacific Hydrate Ridge (HRSPAC47), Loki’s Castle (LOKART) from the Arctic, Santa Monica Mounds (SMMPAC), and the Gulf of Mexico (CIGOMD18 and KJGOM6) (Fig. 5b).

The metagenomes had an extracellular inventory of CAZyme, primarily targeting algal and necromass detritus (see Fig. 5b). Among the prevalent modules engaged in the breakdown of algal debris were glycoside hydrolase modules GH2, GH3, and GH16_3, as well as carbohydrate esterase CE1. The binding modules included CBM9, CBM44, and CBM67.These modules are composed of enzyme families with β-galactosidases, β-glucuronidases, β-mannosidases, exo-β-glucosaminidases activities in the case of GH2 and GH3, where glycoside hydrolases and phosphorylases perform a wide range of functions that involve biomass degradation and remodelling of plant and bacterial cell walls. GH16_3 breaks laminarase, a carbohydrate found in brown algae (Qin et al., 2017) while CE1 has acetylxylan esterases (EC 3.1.1.72), feruloyl esterases (EC 3.1.1.73) activities, and many other esterases such as PHB depolymerases. CM9 and CBM44 are modules targeting cellulose binding domains mainly xylan and other carbohydrates cellulose binding domains and CB67 targets binding to L-rhamnose, a carbohydrate produced by microalgae (0-13.3 of algal composition%) (Brown, 1990) (Fig. 5b) (Lombard et al.,2014).

For necromass degradation, the GH23 and GH103 modules contain families of peptidoglycan lytic transglycosylases. GH23 has also been found to have chitinase activity. Furthermore, known activities of the CE4 and CE14 families include enzymes such as acetylxylan esterases, chitin deacetylases, chitooligosaccharide deacetylases, and peptidoglycan deacetylases (CE4) and diacetylchitobiose deacetylase (EC 3.5.1.-) chitin disaccharide deacetylases (CE14). (Lombard et al., 2014). Finally, for host glycan degradation, the GH29 module contains α-L-fucosidases, and the GH109 modules conform to -N-acetylgalactosaminidase, α-N-acetylgalactosaminidase, and β-N-acetylhexosaminidase. GH33 sialidase or neuraminidase (EC 3.2.1.18) targets the sialic acid of the host glycan (Fig. 5b).

It is documented that bacterial communities dominate shallow sediments, which are primarily composed of clay, cellular envelopes of planktonic organisms, and organic matter (Bienhold et al., 2016). Genes related to the degradation of recalcitrant carbon, including cellulose, chitin, or peptidoglycan, are expected to play an important role in marine sediments (Tully et al., 2016; Bradley et al., 2018; Orsi et al., 2018). Necromass contributes significantly to meeting the energy demand of up to 13% of the microbial community in shallow sediments when it is oxidised under oxic or anoxic conditions. The oxidation of one cell per year can provide sufficient energy to support the demand of thousands of cells in sediments with low energy resources, potentially positioning necromass oxidation as a primary carbon source for microorganisms unable to survive in energy-poor environments (Bradley et al., 2018). The fact that mineralization and adsorption of biopolymers in sediment particles could reduce the accessibility of other carbohydrates (Orsi et al., 2018) could make cell envelopes, such as peptidoglycan, a preferred choice for secreted CAZyme modules found to be the most abundant (Fig. 5a). Most of these CAZyme modules are found across a broad spectrum of life forms but are concentrated in bacteria (Lombard et al., 2014).

Some of the most abundant modules differed between the oxic and anoxic samples. The CAZyme modules GH23, CBM9, GH16_3, GT51, CE4, and CE14 were significantly more abundant in oxic samples. On the other hand, CBM44 and GT83 were found to be different in relation to anoxic samples. Interestingly, despite quite opposite distributions, both CBM44 and CBM9 can bind cellulose (Fig. 5c).

Reconstruction of MAG and their potential to degrade carbohydrates found in marine sediments.

To better understand the community involved in carbohydrate turnover in marine sediments, we recovered MAG from each metagenome sample. We reconstructed 494 high quality MAG (Completeness 75%, Contamination < 10%) that were assigned taxonomically. The classes Alphaproteobacteria, Bacteroidia, and Gammaproteobacteria were among the most represented in our MAG (Supplementary Table 4).

To see whether our recovered MAG are involved in carbohydrate turnover in marine sediments, the CAZymes of the 18 most abundant modules found in our annotations of the metagenome samples were searched and annotated in the MAG. We focused on secreted CAZymes modules and CAZymes modules corresponding to CAZyme Gene Clusters (CGC) (Fig. 6 Supplementary Table 5 and Supplementary Table 6).

The GH23 module was the most abundant in sediment MAG and no CBM44 modules were found in any. MAG from the classes Alphaproteobacteria, Bacteroidia and Gammaproteobacteria had more than one module of extracellular CAZymes (Fig. 6). Alphaproteobacteria MAG had GH103 and GH23 modules. The Alphaproteobacteria MAG belong to the Rhodobacteraceae and Methyloligellaceae families with species found in marine environments such as the genera Pseudorhodobacter, Sulfitobacter, Roseicyclus and Hyphomicrobium (Uchino et al., 2002; Rathgeber et al., 2005; Yoon et al., 2007; Vuilleumier et al., 2011) and other species of the genus Methyloceanibacter who had been previously reported in North Sea sediments (Vekeman et al., 2016) (Supplementary Table 4, Supplementary Table 5).

The Bacteroidia MAG contained CAZyme modules in at least one MAG except for module GT83; CAZyme composition of the main modules found in our metagenome samples (GH109, GH23 and CBM9) was higher in the family Flavobacteriaceae were MAG assigned to the genus Prevotella (DSANT95_maxbin.044), Maribacter (DSANT06_maxbin.002, DSANT06_maxbin.016 and DSANT95_maxbin.030) Pricia (DSANT95_maxbin.051, DSANT95_maxbin.024 and DSANT95_maxbin.016) Eudoraea (DSANT11_maxbin.017), Aureibaculum (DSANT06_maxbin.006 and DSANT08_maxbin.008) along with other abundant modules (Fig. 6, Supplementary Table 4, Supplementary Table 5).

The species Prevotella, Maribacter, and Aureibaculum had been recovered from marine sediments from the Pacific Ocean and Yellow Sea (Reed et al., 2002; Nedashkovskaya et al., 2004; Zhao et al., 2019). The Pricia genus had previously been isolated from a sample of sandy intertidal sediment collected from the Antarctic coast (Yu et al., 2012) which is consistent with the place it was recovered (Davis Station). Eudoraea species were isolated from coastal waters of the Adriatic Sea (Alain et al., 2008).

Finally, MAG without GH23 modules such as the classes of Phycisphaerae, UBA2214, Planctomycetia and Bacteroidia contained GH109, GH2, GH29 and CBM67. UBA2214 was also enriched with GH3 modules. MAG from UBA2214, Phycisphaerae, and Planctomycetia MAG were assigned to the Zgenome-0027, Anaerohalosphaeraceae and Thermoguttaceae families, respectively. Species from these families are found in marine sediments and low oxygen aquatic environments (Dedysh et al., 2020; Pradel et al., 2020; Chiciudean et al., 2022). Furthermore, four Bacteroidia MAG assigned to the Bacteroidales order showed a similar CAZyme inventory (Fig. 6, Supplementary Table 5).

As CAZymes are also known to work in conjunction with other CAZymes and proteins forming CGCs, we decided to search for clusters involving the main CAZyme modules found in our metagenomes.

In general, CGCs targeting the GH23 module often came attached to a CBM50 module; GH3 module often came attached to a CBM6 module and GH2 modules often came with CBM67 modules (Supplementary Table 6). CBM50, a module for the recognition of chitin or peptidoglycan (Ohnuma et al., 2008), has already been abundantly reported in marine sediments (Orsi et al., 2018). CBM6 is known for the recognition of xylanases, lichenases, β-agarases, laminarinases and deacetylases, and CBM67 is known for the recognition of rhamnose, both carbohydrates are found in algal content (Lombard et al., 2014).

MAG belonging to Bacteroidia and Gammaproteobacteria have the highest number of CGCs. Bacteroidia MAG classified as Prevotella (DSANT95_maxbin.044) and Gammaproteobacteria GCA-001735895 sp009937625 (KOR58_maxbin 012.fasta. contigs.refined) had the highest number of CGCs of all, with five including GH3, GH23 and GH2 modules in the case of Prevotella and targeting GH23, GH103, GT51 and CE4 modules. (Supplementary Table 4, Supplementary Table 5, Supplementary Table 6). The Alphaproteobacteria MAG contained CGCs targeting GH23, GH103, and GH3 modules. Gammaproteobacteria CGCs were found targeting CE4, GH103, GT51, CBM9, GH23 and GH3 modules.

Even though assembled MAG cannot cover all sediment diversity, we did find a group of MAG annotated to classes that were abundant in our samples, such as Bacteroidia, Alphaproteobacteria, and Gammaproteobacteria. We did find the CAZyme inventory and CGCs that contained the most abundant modules found in our metagenomes in these classes of bacteria. Furthermore, the MAG we found having important CAZyme modules belong to genera or families found or isolated in marine environments, making these classes some of the main drivers for carbohydrate transformation in marine sediments. It is well known that the Bacteroidota phylum is considered the primary phylum for carbohydrate degradation (Lapébie, et al., 2019). All our MAG from this phylum belonged to Bacteroidia. Phylum Proteobacteria was the most prevalent one in sediment samples. Most of the taxa we found belong to Gamma and Alpha Proteobacteria (47.75–13.88% and 33.83–6.91% of relative abundance, respectively) (Supplementary Table 7, Supplementary Fig. S2).

We successfully identified and analysed the MAG from metagenome samples, shedding light on the key players in carbohydrate turnover in marine sediments. These classes showed the presence of the most abundant CAZyme modules and CAZyme Gene Clusters (CGCs) that correspond to carbohydrate degradation in marine environments. The presence of these CAZymes and CGCs in marine-derived MAG indicates their critical role in carbohydrate transformation in marine sediments.

This highlights the importance of the Bacteroidota phylum in carbohydrate degradation, particularly the Bacteroidia class, and the significant contributions of both Gamma and Alpha Proteobacteria to the observed taxa in marine sediment samples.

CAZyme profile of marine sediment taxa vs. soil sediments.

Since Alphaproteobacteria, Gammaproteobacteria and Bacteroidia had such a rich inventory of CAZyme for carbohydrates found in marine sediments, we decided to explore how different were the CAZyme inventories of our MAG to those of MAG of Alphaproteobacteria, Bacteroidia, and Gammaproteobacteria selected from soil samples published by Nayfach et al., (2021) using the same selection criteria (Completeness > 75% Contamination < 10%).

Marine sediments and soil are rich ecosystems of microorganisms and are crucial components of the Earth's surface, as they can sequestrate carbon and play a role in carbon recycling. (Arndt et al., 2013; Bargett et al., 2014).

MAG from these classes were mainly from different families compared to the sediment MAG we recovered (Supplementary Table 8). PCoA analysis of the counts of all the CAZyme modules found in each MAG showed that the composition of CAZyme appeared to be similar between the phyla where Alfa and Gamma Proteobacteria clustered together and as did Bacteroidia MAG (29.74% of the variance explained in CoA1 and CoA2) (Fig. 7a).

The main difference between the classes recovered from the MAG of sediments compared to the MAG of soil was the number of CAZyme modules found between them and the diversity of the CAZyme modules: all classes of soil MAG had a total number of modules greater (Fig. 7b) and more diverse (Fig. 7c) compared to those we recovered from marine sediments where the Bacteroidia class was the one that had more counts and more diverse CAZyme modules. This is consistent with studies of MAG in environments where CAZyme modules are phylogenetically conserved, among microbial phyla, but some specificity toward habitat is present where soil is an ecosystem where richness in and diversity in CAZyme modules has been found in contrast to marine environments such as marine sediments (López-Mondéjar et al., 2022). Furthermore, the Bacteroidetes phylum to which the Bacteroidia class belongs has been reported as the main class for carbohydrate transformation, as it uses a large inventory of CAZyme (Lapébie et al., 2019).

This comparison between the MAG of marine sediment and soil metagenome-assembled genomes of Alphaproteobacteria, Gammaproteobacteria, and Bacteroidia reveals interesting differences that highlight the contrasting ecological roles and environmental pressures these bacteria experience in their respective habitats. The higher number and diversity of CAZyme modules found in soil MAG compared to marine sediment MAG support the idea that soil microbial communities are exposed to a wider variety of organic substrates, including plant biomass, animal detritus, and complex soil organic matter. This diversity of substrates likely drives the need for a broader suite of enzymatic capabilities in soil microorganisms, as reflected in their CAZyme repertoire.

On the contrary, marine sediment environments may be more homogeneous in terms of organic substrate availability, possibly due to the predominance of marine-derived organic matter, such as phytoplankton and other marine organisms. This could explain why marine sediment MAG possess a less diverse CAZyme profile compared to soil MAG.

Another possible explanation could be that the marine sediment environment is more energy-limited compared to the soil, leading to selective pressure for organisms that can efficiently degrade available organic matter with a smaller set of enzymes. This could potentially lead to a more streamlined CAZyme profile in marine sediment bacteria.

Despite these differences, the fact that Alphaproteobacteria, Gammaproteobacteria, and Bacteroidia from soil and marine sediments cluster together in the PCoA analysis suggests a core set of CAZyme modules that are conserved within these taxonomic groups, likely reflecting shared evolutionary histories and core metabolic functions.

This study underscores the importance of considering the ecological context when studying the functional capabilities of microbial communities. The stark differences in the CAZyme profiles between soil and marine sediment bacteria underscore how environmental factors can shape the functional potential of microbial communities. Therefore, it is essential to take these factors into account when studying the ecology and function of microorganisms in different environments.

In this study, we classified 37 metagenomes from around the world with few physicochemical metadata by comparing the community's potential to use oxygen as the last electron acceptor as a marker to classify them as oxic or anoxic. We find a clear difference between our sediment samples in terms of taxonomy and CAZyme content in the context of this classification. We established a profile of the most abundant extracellular CAZymes in our samples, where 18 CAZyme modules were abundant in all and were found to primarily target carbohydrates from necromass degradation and algae detritus, which is consistent with the environmental conditions found in sediments. We find significantly different modules targeting the same substrate depending on oxic and anoxic conditions. Most of the main abundant CAZyme modules that we found were of bacterial origin.

Finally, we recovered MAG from the samples, which were assigned to the classes Alphaproteobacteria, Gammaproteobacteria, and Bacteroidia. The MAG contained extracellular modules of the main CAZymes that were also annotated in our metagenomes, as well as CGCs that had those modules as part of the CAZyme machinery. Module GH23, which targets peptidoglycan and chitin substrates, was found in almost all our MAG. The MAG that did not contain the GH23 module had other different main modules that target host glycan and plant detritus. These taxa are the bacteria that mainly drive carbohydrate transformation in marine sediments, although further studies are needed to fully confirm this. Our findings provide valuable information on the community structure and function of carbohydrate turnover in marine sediments, highlighting the key roles that specific bacterial classes play and their associated CAZyme inventories and CGCs.

It is important to note that many of the MAG we reconstructed belonged to taxa already found in marine environments and that the classes Bacteroidia, Alpha, and Gammaproteobacteria are found in our metagenomes in abundance. When we compare the MAG from sediments with other MAG from the same taxa assembled from soil samples, we find a similar profile, as expected from taxonomy, but with fewer total and diverse CAZyme modules in sediment MAG. This is a response to the oligotrophic conditions in marine sediments, in contrast to soil conditions. Although the MAGs from our samples give us a glimpse of the microbial community, the number of recovered MAGs is not sufficient to encapsulate the great diversity of the microbial community of marine sediments. Furthermore, it is necessary to study deep subsurface sediments to better understand the CAZyme inventory compared to shallow marine sediments.

AA Auxiliary Activities

CE Carbohydrate Esterase

CBM Carbohydrate Binding Module

GH Glycoside Hydrolase

GT GlycosylTransferase

CAZymes Carbohydrate-active enzymes

CGC CAZyme Gene Cluster

MAG metagenomic assembly genomes

Mbsl Meters below sea level

Mbsf Meters below sea floor

PL Polysaccharide Lyase

Author Contributions.

Rafael López-Sánchez and Lorenzo Segovia conceived the study, selected samples, performed bioinformatic analyses, and wrote the manuscript. Eria A. Rebollar, Rosa María Gutiérrez-Ríos, Alejandro Garcíarrubio and Katy Juárez helped design the details of the study and contributed throughout. All authors read and approved the final version of the manuscript.

Funding Information.

This work was funded by UNAM grants from the Dirección General de Asuntos del Personal Académico (DGAPA) IN209921 and IV200322. Support was also provided by an infrastructure grant from the “Programa de Apoyos al Fortalecimiento y Desarrollo de la Infraestructura Científica y Tecnológica de CONACyT” “Fortalecimiento de la Infraestructura del grupo de Biocatálisis en el IBt/UNAM” Grant number 269435. RLS was a recipient of a Consejo Nacional de Ciencia y Tecnología scholarship grant.

Data Availability.

The authors confirm that all supporting data, code, and protocols have been provided within the article or through supplementary data files. Metagenomes from Illumina sequenced shotgun sediment samples were obtained from NCBI (https://www.ncbi.nlm.nih.gov/genome) and from the Instituto de Biotecnología de la UNAM. All details and accession numbers from the samples can be found in the Supplementary Material and GitHub repository: https://github.com/RafaelLopez-Sanchez/marine_sediments

Conflicts of interest.

The authors declare that there are no conflicts of interest.

Alain, K., Intertaglia, L., Catala, P., and Lebaron, P. (2008). Eudoraea adriatica gen. nov., sp. nov., a novel marine bacterium of the family Flavobacteriaceae. Int. J. Syst. Evol. Microbiol. 58, 2275–2281. doi:10.1099/ijs.0.65446-0.
Almagro Armenteros, J. J., Tsirigos, K. D., Sønderby, C. K., Petersen, T. N., Winther, O., Brunak, S., et al. (2019). SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat. Biotechnol. 37, 420–423. doi:10.1038/s41587-019-0036-z.
Andrews, S. FastQC: a quality control tool for high throughput sequence data. (2010). Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc.
Arndt, S., Jørgensen, B. B., LaRowe, D. E., Middelburg, J. J., Pancost, R. D., and Regnier, P. (2013). Quantifying the degradation of organic matter in marine sediments: A review and synthesis. Earth-Science Rev. 123, 53–86. doi:10.1016/j.earscirev.2013.02.008.
Bäckström D et al., "Virus Genomes from Deep Sea Sediments Expand the Ocean Megavirome and Support Independent Origins of Viral Gigantism.", mBio, 2019 Mar 5;10(2).
Bardgett, R. D., and Van Der Putten, W. H. (2014). Belowground biodiversity and ecosystem functioning. Nature 515, 505–511. doi:10.1038/nature13855.
Biddle, J. F., Fitz-gibbon, S., Schuster, S. C., Brenchley, J. E., and House, C. H. (2008). Metagenomic signatures of the Peru Margin subseafloor biosphere show a genetically distinct environment. 105, 10583–10588. doi:10.1073/pnas.0709942105.
Bienhold, C., Zinger, L., Boetius, A., and Ramette, A. (2016). Diversity and biogeography of bathyal and abyssal seafloor bacteria. PLoS One 11, 1–20. doi:10.1371/journal.pone.0148016.
Biersmith, A., and Benner, R. (1998). Carbohydrates in phytoplankton and freshly produced dissolved organic matter. Mar. Chem. 63, 131–144. doi:10.1016/S0304-4203(98)00057-7.
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi:10.1093/bioinformatics/btu170.
Bradley, J. A., Amend, J. P., and Larowe, D. E. (2018). for Microorganisms in Marine Sediments. 577–590.
Chiciudean, I., Russo, G., Bogdan, D. F., Levei, E. A., Faur, L., Hillebrand-Voiculescu, A., et al. (2022). Competition-cooperation in the chemoautotrophic ecosystem of Movile Cave: first metagenomic approach on sediments. Environ. Microbiomes 17, 1–18. doi:10.1186/s40793-022-00438-w.
Chaumeil, P.-A., Mussig, A. J., Hugenholtz, P., and Parks, D. H. (2020). GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics 36, 1925–1927. doi:10.1093/bioinformatics/btz848.
D’Hondt, S. D., Jørgensen, B. B., Miller, D. J., Batzke, A., Blake, R., Cragg, B. A., et al. (2004). Esearch rticle. 2216–2222.
D’hondt, S., Inagaki, F., Zarikian, C. A., Abrams, L. J., Dubois, N., Engelhardt, T., et al. (2015). Presence of oxygen and aerobic communities from sea floor to basement in deep-sea sediments. Nat. Geosci. 8, 299–304. doi:10.1038/ngeo2387.
Dedysh, S. N., Kulichevskaya, I. S., Beletsky, A. V., Ivanova, A. A., Rijpstra, W. I. C., Damsté, J. S. S., et al. (2020). Lacipirellula parvula gen. nov., sp. nov., representing a lineage of planctomycetes widespread in low-oxygen habitats, description of the family Lacipirellulaceae fam. nov. and proposal of the orders Pirellulales ord. nov., Gemmatales ord. nov. and Isosphaerales ord. nov. Syst. Appl. Microbiol. 43, 126050. doi:10.1016/j.syapm.2019.126050.
Eddy, S. R. (2011). Accelerated Profile HMM Searches. PLOS Computational Biology, 7(10), e1002195. doi:10.1371/journal.pcbi.1002195
Froelich, P. N., Klinkhammer, G. P., Bender, M. L., Luedtke, N. A., Heath, G. R., Cullen, D., et al. (1979). Early oxidation of organic matter in pelagic sediments of the eastern equatorial Atlantic: suboxic diagenesis. Geochim. Cosmochim. Acta 43, 1075–1090. doi:10.1016/0016-7037(79)90095-4.
Godoy-Lozano, E. E. et al. Bacterial diversity and the geochemical landscape in the southwestern Gulf of Mexico. Front. Microbiol. 9, 1–15 (2018).
Hijmans, R. J. (2020). Geosphere: Spherical Trigonometry [R package]. Available at: https://CRAN.R-project.org/package=geosphere.
Hyatt D, Chen GL, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11. doi:10.1186/1471-2105-11-119.
Hoshino, T., Doi, H., Uramoto, G. I., Wörmer, L., Adhikari, R. R., Xiao, N., et al. (2020). Global diversity of microbial communities in marine sediment. Proc. Natl. Acad. Sci. U. S. A. 117, 27587–27597. doi:10.1073/pnas.1919139117.
Jaeschke, A., Jørgensen, S. L., Bernasconi, S. M., Pedersen, R. B., Thorseth, I. H., and Früh-Green, G. L. (2012). Microbial diversity of Loki’s Castle black smokers at the Arctic Mid-Ocean Ridge. Geobiology 10, 548–561. doi:10.1111/gbi.12009.
Jiao, N., Herndl, G. J., Hansell, D. A., Benner, R., Kattner, G., Wilhelm, S. W., et al. (2010). Microbial production of recalcitrant dissolved organic matter: Long-term carbon storage in the global ocean. Nat. Rev. Microbiol. 8, 593–599. doi:10.1038/nrmicro2386.Kallmeyer J, Pockalny R, Adhikari RR, Smith DC,
D’Hondt S. Global distribution of microbial abundance and biomass in subseafloor sediment. Proc Natl Acad Sci U S A. 2012;109(40):16213-16216. doi:10.1073/pnas.1203849109.
Kauffman, K. M., Hussain, F. A., Yang, J., Arevalo, P., Brown, J. M., Chang, W. K., et al. (2018). A major lineage of non-tailed dsDNA viruses as unrecognized killers of marine bacteria. Nature 554, 118–122. doi:10.1038/nature25474.
Kultima, J. R., Sunagawa, S., Li, J., Chen, W., Chen, H., Mende, D. R., et al. (2012). MOCAT: A Metagenomics Assembly and Gene Prediction Toolkit. PLoS One 7, 1–6. doi:10.1371/journal.pone.0047656.
Lapébie, P., Lombard, V., Drula, E., Terrapon, N., and Henrissat, B. (2019). Bacteroidetes use thousands of enzyme combinations to break down glycans. Nat. Commun. 10, 2043. doi:10.1038/s41467-019-10068-5.
Lee, C., Hedges, J. I., Baldock, J. A., Ge, Y., Gelinas, Y., Peterson, M., et al. (2001). Evidence for non-selective preservation of organic matter in sinking marine particles. Nature 409, 801–804.
Leeming, R., Stark, J., & Smith, J. (2015). Novel use of faecal sterols to assess human faecal contamination in Antarctica: A likelihood assessment matrix for environmental monitoring. Antarctic Science, 27(1), 31-43. doi:10.1017/S0954102014000273.
Li, D., Liu, C-M., Luo, R., Sadakane, K., and Lam, T-W., (2015) MEGAHIT: An ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics, doi: 10.1093/bioinformatics/btv033.
Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42(D1):490-495. doi:10.1093/nar/gkt1178. CAZy Database. (2023). Carbohydrate-Active enZYmes Database. Retrieved from http://www.cazy.org.
Lombard, V., Golaconda Ramulu, H., Drula, E., Coutinho, P. M., & Henrissat, B. (2014). The carbohydrate-active enzymes database (CAZy) in 2013. In Cazypedia, the Carbohydrate-Active enZYmes Encyclopedia(2022). Retrieved from http://www.cazypedia.org
López-Mondéjar, R., Tláskal, V., da Rocha, U. N., and Baldrian, P. (2022). Global Distribution of Carbohydrate Utilization Potential in the Prokaryotic Tree of Life. mSystems 7. doi:10.1128/msystems.00829-22.
National Center for Biotechnology Information (2019). PubMed. In National Center for Biotechnology Information. Retrieved from https://www.ncbi.nlm.nih.gov/pubmed/
Nayfach, S., Roux, S., Seshadri, R., Udwary, D., Varghese, N., Schulz, F., et al. (2021). A genomic catalog of Earth’s microbiomes. Nat. Biotechnol. 39, 499–509. doi:10.1038/s41587-020-0718-6.
Nedashkovskaya, O. I., Kim, S. B., Han, S. K., Lysenko, A. M., Rohde, M., Rhee, M. S., et al. (2004). Maribacter gen. nov., a new member of the family Flavobacteriaceae, isolated from marine habitats, containing the species Maribacter sedimenticola sp. nov., Maribacter aquivirus sp. nov., Maribacter orientalis sp. nov. and Maribacter ulvicola sp. nov. Int. J. Syst. Evol. Microbiol. 54, 1017–1023. doi:10.1099/ijs.0.02849-0.
Ohnuma, T., Onaga, S., Murata, K., Taira, T., and Katoh, E. (2008). LysM domains from Pteris ryukyuensis chitinase-A: A stability study and characterization of the chitin-binding site. J. Biol. Chem. 283, 5178–5187. doi:10.1074/jbc.M707156200.
Orsi, W. D. (2018). Ecology and evolution of seafloor and subseafloor microbial communities. Nat. Rev. Microbiol. 16, 671–683. doi:10.1038/s41579-018-0046-8.
Orsi, W. D., Richards, T. A., and Francis, W. R. (2018). Predicted microbial secretomes and their target substrates in marine sediment. Nat. Microbiol. 3, 32–37. doi:10.1038/s41564-017-0047-9.
Oksanen, F.J., et al. (2017) Vegan: Community Ecology Package. R package Version 2.4-3. https://CRAN.R-project.org/package=vegan
Orcutt, B. N., Sylvan, J. B., Knab, N. J. & Edwards, K. J.Microbial ecology of the dark ocean above, at, and below the seafloor. Microbiol. Mol. Biol. Rev. 75, 361–422 (2011).
Pante, E., and Simon-Bouhet, B. (2013). marmap: A Package for Importing, Plotting and Analyzing Bathymetric and Topographic Data in R. PLoS One 8, e73051. Available at: https://doi.org/10.1371/journal.pone.0073051.
Parkes, R. J., Cragg, B., Roussel, E., Webster, G., Weightman, A., and Sass, H. (2014). A review of prokaryotic populations and processes in sub-seafloor sediments, including biosphere: Geosphere interactions. Mar. Geol. 352, 409–425. doi:10.1016/j.margeo.2014.02.009.
Pradel, N., Fardeau, M. L., Tindall, B. J., and Spring, S. (2020). Anaerohalosphaera lusitana gen. nov., sp. nov., and limihaloglobus sulfuriphilus gen. nov., sp. nov., isolated from solar saltern sediments, and proposal of anaerohalosphaeraceae fam. nov. within the order sedimentisphaerales. Int. J. Syst. Evol. Microbiol. 70, 1321–1330. doi:10.1099/ijsem.0.003919.
Qin, H. M., Miyakawa, T., Inoue, A., Nakamura, A., Nishiyama, R., Ojima, T., et al. (2017). Laminarinase from Flavobacterium sp. reveals the structural basis of thermostability and substrate specificity. Sci. Rep. 7, 1–9. doi:10.1038/s41598-017-11542-0.
R Core Team (2023). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/.
Rathgeber, C., Yurkova, N., Stackebrandt, E., Schumann, P., Beatty, J. T., and Yurkov, V. (2005). Roseicyclus mahoneyensis gen. nov., sp. nov., an aerobic phototrophic bacterium isolated from a meromictic lake. Int. J. Syst. Evol. Microbiol. 55, 1597–1603. doi:10.1099/ijs.0.63195-0.
Reed, D. W., Fujita, Y., Delwiche, M. E., Blackwelder, D. B., Sheridan, P. P., Uchida, T., et al. (2002). Microbial communities from methane hydrate-bearing deep marine sediments in a forearc basin. Appl. Environ. Microbiol. 68, 3759–3770. doi:10.1128/AEM.68.8.3759-3770.2002.
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., & Huttenhower, C. (2011). Metagenomic biomarker discovery and explanation. Genome Biology, 12(6), R60. Doi: 10.1186/gb-2011-12-6-r60.
Sousa, F. L., Alves, R. J., Pereira-Leal, J. B., Teixeira, M., and Pereira, M. M. (2011). A bioinformatics classifier and database for Heme-Copper oxygen reductases. PLoS One 6, 1–9. doi:10.1371/journal.pone.0019117.
Tamames, J., Puente-Sánchez, F., SqueezeMeta Development Team, et al. (2019). SqueezeMeta: a highly portable, fully automatic metagenomic analysis pipeline. PeerJ, 7, e7559. Doi: 10.7717/peerj.7559.
Teeling, H., Fuchs, B. M., Becher, D., Klockow, C., Gardebrecht, A., Bennke, C. M., et al. (2012). Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science (80-.). 336, 608–611. Doi:10.1126/science.1218344.
Teeling, H., Fuchs, B. M., Bennke, C. M., Krüger, K., Chafee, M., Kappelmann, L., et al. (2016). Recurring patterns in bacterioplankton dynamics during coastal spring algae blooms. Elife 5, 1–31. doi:10.7554/eLife.11888.
Tully, B. J., and Heidelberg, J. F. (2016). Potential mechanisms for microbial energy acquisition in oxic deep- sea sediments. Appl. Environ. Microbiol. 82, 4232–4243. doi:10.1128/AEM.01023-16.
Uchino, Y., Hamada, T., and Yokota, A. (2002). Proposal of Pseudorhodobacter ferrugineus gen. nov., comb. nov., for a non-photosynthetic marine bacterium, Agrobacterium ferrugineum, related to the genus Rhodobacter. J. Gen. Appl. Microbiol. 48, 309–319. doi:10.2323/jgam.48.309.
Vekeman, B., Kerckhof, F. M., Cremers, G., de Vos, P., Vandamme, P., Boon, N., et al. (2016). New Methyloceanibacter diversity from North Sea sediments includes methanotroph containing solely the soluble methane monooxygenase. Environ. Microbiol. 18, 4523–4536. doi:10.1111/1462-2920.13485.
Vuilleumier, S., Nadalig, T., Ul Haque, M. F., Magdelenat, G., Lajus, A., Roselli, S., et al. (2011). Complete genome sequence of the chloromethane-degrading Hyphomicrobium sp. Strain MC1. J. Bacteriol. 193, 5035–5036. doi:10.1128/JB.05627-11.
Wood, D. E., Lu, J., & Langmead, B. (2019). Improved metagenomic analysis with Kraken 2. Genome Biology, 20(1), 257. doi:10.1186/s13059-019-1891-0 Yu, Y., Li, H. R., Zeng, Y. X., Sun, K., and Chen, B. (2012). Pricia antarctica gen. nov., sp. nov., a member of the family flavobacteriaceae, isolated from antarctic intertidal sediment. Int. J. Syst. Evol. Microbiol. 62, 2218–2223. doi:10.1099/ijs.0.037515-0.
Yoon, J. H., Kang, S. J., Lee, M. H., and Oh, T. K. (2007). Description of Sulfitobacter donghicola sp. nov., isolated from seawater of the East Sea in Korea, transfer of Staleya guttiformis Labrenz et al. 2000 to the genus Sulfitobacter as Sulfitobacter guttiformis comb. nov. and emended description of the genus Sulfitobacter. Int. J. Syst. Evol. Microbiol. 57, 1788–1792. doi:10.1099/ijs.0.65071-0.
Zhang, H., Yohe, T., Huang, L., Entwistle, S., Wu, P., Yang, Z., et al. (2018). DbCAN2: A meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 46, W95–W101. doi:10.1093/nar/gky418.
Zhao, H., Wu, Y., Zhang, C., Feng, J., Xu, Z., Ding, Y., et al. (2019). Aureibaculum marinum gen. nov., sp. nov., a Novel Bacterium of the Family Flavobacteriaceae Isolated from the Bohai Gulf. Curr. Microbiol. 76, 975–981. doi:10.1007/s00284-019-01691-y.
Zhao, R., Summers, Z.M., Christman, G.D. et al. Metagenomic views of microbial dynamics influenced by hydrocarbon seepage in sediments of the Gulf of Mexico. Sci Rep 10, 5772 (2020). https://doi.org/10.1038/s41598-020-62840-z

No competing interests reported.

Download PDF

Journal Publication

published 13 Feb, 2024

Read the published version in World Journal of Microbiology and Biotechnology →

Editorial decision: Major revision
24 Oct, 2023
Reviews received at journal
23 Oct, 2023
Reviews received at journal
23 Oct, 2023
Reviewers agreed at journal
20 Oct, 2023
Reviewers agreed at journal
25 Sep, 2023
Reviewers invited by journal
03 Jul, 2023
Editor assigned by journal
01 Jul, 2023
Submission checks completed at journal
30 Jun, 2023
First submitted to journal
30 Jun, 2023

You are reading this latest preprint version

Metagenomic Analysis of Carbohydrate-Active Enzymes and Their Contribution to Marine Sediment Biodiversity

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Methods

Results and Discussion

Conclusion

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1