Metagenomic views of microbial dynamics influenced by hydrocarbon seepage in sediments of the Gulf of Mexico

Active microbial communities facilitate carbon turnover in brine pools found in the deep Southeastern Mediterranean Sea

Discharge of gas-rich brines fuels productive chemosynthetic ecosystems in the deep sea. In these salty, methanic and sulfidic brines, microbial communities adapt to specific niches along the physicochemical gradients. However, the molecular mechanisms that underpin these adaptations are not fully known. Using metagenomics, we investigated the dense (∼106 cell ml-1) microbial communities that occupy small deep-sea brine pools found in the Southeastern Mediterranean Sea (1150 m water depth, ∼22 °C, ∼60 PSU salinity, sulfide, methane, ammonia reaching millimolar levels, and oxygen usually depleted), reaching high productivity rates of 685 μg C L-1 d-1 ex-situ. We curated 266 metagenome-assembled genomes of bacteria and archaea from the several pools and adjacent sediment-water interface, highlighting the dominance of a single Sulfurimonas, which likely fuels its autotrophy using sulfide oxidation or inorganic sulfur disproportionation. This lineage may be dominant in its niche due to genome streamlining, limiting its metabolic repertoire, particularly by using a single variant of sulfide: quinone oxidoreductase. These primary producers co-exist with ANME-2c archaea that catalyze the anaerobic oxidation of methane. Other lineages can degrade the necromass aerobically (Halomonas and Alcanivorax), or anaerobically through fermentation of macromolecules (e.g., Caldatribacteriota, Bipolaricaulia, Chloroflexota, etc). These low-abundance organisms likely support the autotrophs, providing energy-rich H2, and vital organics such as vitamin B12.

Microbial community response to hydrocarbon exposure in iron oxide mats: an environmental study

Hydrocarbon pollution is a widespread issue in both groundwater and surface-water systems; however, research on remediation at the interface of these two systems is limited. This interface is the oxic-anoxic boundary, where hydrocarbon pollutant from contaminated groundwaters flows into surface waters and iron mats are formed by microaerophilic iron-oxidizing bacteria. Iron mats are highly chemically adsorptive and host a diverse community of microbes. To elucidate the effect of hydrocarbon exposure on iron mat geochemistry and microbial community structure and function, we sampled iron mats both upstream and downstream from a leaking underground storage tank. Hydrocarbon-exposed iron mats had significantly higher concentrations of oxidized iron and significantly lower dissolved organic carbon and total dissolved phosphate than unexposed iron mats. A strong negative correlation between dissolved phosphate and benzene was observed in the hydrocarbon-exposed iron mats and water samples. There were positive correlations between iron and other hydrocarbons with benzene in the hydrocarbon-exposed iron mats, which was unique from water samples. The hydrocarbon-exposed iron mats represented two types, flocculent and seep, which had significantly different concentrations of iron, hydrocarbons, and phosphate, indicating that iron mat is also an important context in studies of freshwater mats. Using constrained ordination, we found the best predictors for community structure to be dissolved oxygen, pH, and benzene. Alpha diversity and evenness were significantly lower in hydrocarbon-exposed iron mats than unexposed mats. Using 16S rDNA amplicon sequences, we found evidence of three putative nitrate-reducing iron-oxidizing taxa in microaerophile-dominated iron mats (Azospira, Paracoccus, and Thermomonas). 16S rDNA amplicons also indicated the presence of taxa that are associated with hydrocarbon degradation. Benzene remediation-associated genes were found using metagenomic analysis both in exposed and unexposed iron mats. Furthermore, the results indicated that season (summer vs. spring) exacerbates the negative effect of hydrocarbon exposure on community diversity and evenness and led to the increased abundance of numerous OTUs. This study represents the first of its kind to attempt to understand how contaminant exposure, specifically hydrocarbons, influences the geochemistry and microbial community of freshwater iron mats and further develops our understanding of hydrocarbon remediation at the land-water interface.

A vast repertoire of secondary metabolites potentially influences community dynamics and biogeochemical processes in cold seeps

In deep-sea cold seeps, microbial communities thrive on the geological seepage of hydrocarbons and inorganic compounds, differing from photosynthetically driven ecosystems. However, their biosynthetic capabilities remain largely unexplored. Here, we analyzed 81 metagenomes, 33 metatranscriptomes, and 7 metabolomes derived from nine different cold seep areas to investigate their secondary metabolites. Cold seep microbiomes encode diverse and abundant biosynthetic gene clusters (BGCs). Most BGCs are affiliated with understudied bacteria and archaea, including key mediators of methane and sulfur cycling. The BGCs encode diverse antimicrobial compounds that potentially shape community dynamics and various metabolites predicted to influence biogeochemical cycling. BGCs from key players are widely distributed and highly expressed, with their abundance and expression levels varying with sediment depth. Sediment metabolomics reveals unique natural products, highlighting uncharted chemical potential and confirming BGC activity in these sediments. Overall, these results demonstrate that cold seep sediments serve as a reservoir of hidden natural products and sheds light on microbial adaptation in chemosynthetically driven ecosystems.

Metagenomic analysis of carbohydrate-active enzymes and their contribution to marine sediment biodiversity

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 sites. 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 that extracellular CAZyme modules targeting the degradation of plant and algal detritus, necromass, and host glycans were abundant 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 metagenome assembled 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 discovered in our metagenomes were prevalent in the Bacteroidia, Gammaproteobacteria, and Alphaproteobacteria classes. We compared those MAGs to organisms from the same taxonomic classes found in soil, and we found that they were similar in its CAZyme repertoire, 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.

Metagenome-assembled genomes of deep-sea sediments: changes in microbial functional potential lag behind redox transitions

Hadal sediments are hotspots of microbial activity in the deep sea and exhibit strong biogeochemical gradients. But although these gradients are widely assumed to exert selective forces on hadal microbial communities, the actual relationship between biogeochemistry, functional traits, and microbial community structure remains poorly understood. We tested whether the biogeochemical conditions in hadal sediments select for microbes based on their genomic capacity for respiration and carbohydrate utilization via a metagenomic analysis of over 153 samples from the Atacama Trench region (max. depth = 8085 m). The obtained 1357 non-redundant microbial genomes were affiliated with about one-third of all known microbial phyla, with more than half belonging to unknown genera. This indicated that the capability to withstand extreme hydrostatic pressure is a phylogenetically widespread trait and that hadal sediments are inhabited by diverse microbial lineages. Although community composition changed gradually over sediment depth, these changes were not driven by selection for respiratory or carbohydrate degradation capability in the oxic and nitrogenous zones, except in the case of anammox bacteria and nitrifying archaea. However, selection based on respiration and carbohydrate degradation capacity did structure the communities of the ferruginous zone, where aerobic and nitrogen respiring microbes declined exponentially (half-life = 125-419 years) and were replaced by subsurface communities. These results highlight a delayed response of microbial community composition to selective pressure imposed by redox zonation and indicated that gradual changes in microbial composition are shaped by the high-resilience and slow growth of microbes in the seafloor.

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.

Evolutionary ecology of microbial populations inhabiting deep sea sediments associated with cold seeps

Deep sea cold seep sediments host abundant and diverse microbial populations that significantly influence biogeochemical cycles. While numerous studies have revealed their community structure and functional capabilities, little is known about genetic heterogeneity within species. Here, we examine intraspecies diversity patterns of 39 abundant species identified in sediment layers down to 430 cm below the sea floor across six cold seep sites. These populations are grouped as aerobic methane-oxidizing bacteria, anaerobic methanotrophic archaea and sulfate-reducing bacteria. Different evolutionary trajectories are observed at the genomic level among these physiologically and phylogenetically diverse populations, with generally low rates of homologous recombination and strong purifying selection. Functional genes related to methane (pmoA and mcrA) and sulfate (dsrA) metabolisms are under strong purifying selection in most species investigated. These genes differ in evolutionary trajectories across phylogenetic clades but are functionally conserved across sites. Intrapopulation diversification of genomes and their mcrA and dsrA genes is depth-dependent and subject to different selection pressure throughout the sediment column redox zones at different sites. These results highlight the interplay between ecological processes and the evolution of key bacteria and archaea in deep sea cold seep extreme environments, shedding light on microbial adaptation in the subseafloor biosphere.

Nitrate determines the bacterial habitat specialization and impacts microbial functions in a subsurface karst cave

Karst caves are usually considered as natural laboratories to study pristine microbiomes in subsurface biosphere. However, effects of the increasingly detected nitrate in underground karst ecosystem due to the acid rain impact on microbiota and their functions in subsurface karst caves have remained largely unknown. In this study, samples of weathered rocks and sediments were collected from the Chang Cave, Hubei province and subjected to high-throughput sequencing of 16S rRNA genes. The results showed that nitrate significantly impacted bacterial compositions, interactions, and functions in different habitats. Bacterial communities clustered according to their habitats with distinguished indicator groups identified for each individual habitat. Nitrate shaped the overall bacterial communities across two habitats with a contribution of 27.2%, whereas the pH and TOC, respectively, structured bacterial communities in weathered rocks and sediments. Alpha and beta diversities of bacterial communities increased with nitrate concentration in both habitats, with nitrate directly affecting alpha diversity in sediments, but indirectly on weathered rocks by lowering pH. Nitrate impacted more on bacterial communities in weathered rocks at the genus level than in sediments because more genera significantly correlated with nitrate concentration in weathered rocks. Diverse keystone taxa involved in nitrogen cycling were identified in the co-occurrence networks such as nitrate reducers, ammonium-oxidizers, and N₂-fixers. Tax4Fun2 analysis further confirmed the dominance of genes involved in nitrogen cycling. Genes of methane metabolism and carbon fixation were also dominant. The dominance of dissimilatory and assimilatory nitrate reduction in nitrogen cycling substantiated nitrate impact on bacterial functions. Our results for the first time revealed the impact of nitrate on subsurface karst ecosystem in terms of bacterial compositions, interactions, and functions, providing an important reference for further deciphering the disturbance of human activities on the subsurface biosphere.

Biogeochemical Activity of Methane-Related Microbial Communities in Bottom Sediments of Cold Seeps of the Laptev Sea

Bottom sediments at methane discharge sites of the Laptev Sea shelf were investigated. The rates of microbial methanogenesis and methane oxidation were measured, and the communities responsible for these processes were analyzed. Methane content in the sediments varied from 0.9 to 37 µmol CH₄ dm^-3. Methane carbon isotopic composition (δ¹³C-CH₄) varied from -98.9 to -77.6‱, indicating its biogenic origin. The rates of hydrogenotrophic methanogenesis were low (0.4-5.0 nmol dm^-3 day^-1). Methane oxidation rates varied from 0.4 to 1.2 µmol dm^-3 day^-1 at the seep stations. Four lineages of anaerobic methanotrophic archaea (ANME) (1, 2a-2b, 2c, and 3) were found in the deeper sediments at the seep stations along with sulfate-reducing Desulfobacteriota. The ANME-2a-2b clade was predominant among ANME. Aerobic ammonium-oxidizing Crenarchaeota (family Nitrosopumilaceae) predominated in the upper sediments along with heterotrophic Actinobacteriota and Bacteroidota, and mehtanotrophs of the classes Alphaproteobacteria (Methyloceanibacter) and Gammaproteobacteria (families Methylophilaceae and Methylomonadaceae). Members of the genera Sulfurovum and Sulfurimonas occurred in the sediments of the seep stations. Mehtanotrophs of the classes Alphaproteobacteria (Methyloceanibacter) and Gammaproteobacteria (families Methylophilaceae and Methylomonadaceae) occurred in the sediments of all stations. The microbial community composition was similar to that of methane seep sediments from geographically remote areas of the global ocean.

Phylogenetically and catabolically diverse diazotrophs reside in deep-sea cold seep sediments

Microbially mediated nitrogen cycling in carbon-dominated cold seep environments remains poorly understood. So far anaerobic methanotrophic archaea (ANME-2) and their sulfate-reducing bacterial partners (SEEP-SRB1 clade) have been identified as diazotrophs in deep sea cold seep sediments. However, it is unclear whether other microbial groups can perform nitrogen fixation in such ecosystems. To fill this gap, we analyzed 61 metagenomes, 1428 metagenome-assembled genomes, and six metatranscriptomes derived from 11 globally distributed cold seeps. These sediments contain phylogenetically diverse nitrogenase genes corresponding to an expanded diversity of diazotrophic lineages. Diverse catabolic pathways were predicted to provide ATP for nitrogen fixation, suggesting diazotrophy in cold seeps is not necessarily associated with sulfate-dependent anaerobic oxidation of methane. Nitrogen fixation genes among various diazotrophic groups in cold seeps were inferred to be genetically mobile and subject to purifying selection. Our findings extend the capacity for diazotrophy to five candidate phyla (Altarchaeia, Omnitrophota, FCPU426, Caldatribacteriota and UBA6262), and suggest that cold seep diazotrophs might contribute substantially to the global nitrogen balance.

Microbial nitrogen fixation could be important in the deep sea. Here the authors investigate metagenomes and metatranscriptomes of diazotrophs from deep sea cold seep sediments, reveal greater phylogenetic and functional diversity than hitherto known.

Metagenomic mining of Indian river confluence reveal functional microbial community with lignocelluloytic potential

Microbial carbohydrate-active enzymes (CAZyme) can be harnessed for valorization of Lignocellulosic biomass (LCB) to value-added chemicals/products. The two Indian Rivers Ganges and the Yamuna having different origins and flow, face accumulation of carbon-rich substrates due to the discharge of wastewater from adjoining paper and pulp industries, which could potentially contribute to the natural enrichment of LCB utilizing genes, especially at their confluence. We analyzed CAZyme diversity in metagenomic datasets across the sacred confluence of the Rivers Ganges and Yamuna. Functional annotation using CAZyme database identified a total of 77,815 putative genes with functional domains involved in the catalysis of carbohydrate degradation or synthesis of glycosidic bonds. The metagenomic analysis detected ~ 41% CAZymes catalyzing the hydrolysis of lignocellulosic biomass polymers- cellulose, hemicellulose, lignin, and pectin. The Beta diversity analysis suggested higher CAZyme diversity at downstream region of the river confluence, which could be useful niche for culture-based studies. Taxonomic origin for CAZymes revealed the predominance of bacteria (97%), followed by archaea (1.67%), Eukaryota (0.63%), and viruses (0.7%). Metagenome guided CAZyme diversity of the microflora spanning across the confluence of Ganges-Yamuna River, could be harnessed for biomass and bioenergy applications.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03190-7.

Microbial Nitrogen Transformation Potential in Sediments of Two Contrasting Lakes Is Spatially Structured but Seasonally Stable

The nitrogen (N) cycle is of global importance, as N is an essential element and a limiting nutrient in terrestrial and aquatic ecosystems. Excessive anthropogenic N fertilizer usage threatens sensitive downstream aquatic ecosystems. Although freshwater lake sediments remove N through various microbially mediated processes, few studies have investigated the microbial communities involved. In an integrated biogeochemical and microbiological study on a eutrophic and oligotrophic lake, we estimated N removal rates from pore water concentration gradients in sediments. Simultaneously, the abundance of different microbial N transformation genes was investigated using metagenomics on a seasonal and spatial scale. We observed that contrasting nutrient concentrations in sediments were associated with distinct microbial community compositions and significant differences in abundances of various N transformation genes. For both characteristics, we observed a more pronounced spatial than seasonal variability within each lake. The eutrophic Lake Baldegg showed a higher denitrification potential with higher nosZ gene (N₂O reductase) abundances and higher nirS:nirK (nitrite reductase) ratios, indicating a greater capacity for complete denitrification. Correspondingly, this lake had a higher N removal efficiency. The oligotrophic Lake Sarnen, in contrast, had a higher potential for nitrification. Specifically, it harbored a high abundance of Nitrospira, including some with the potential for comammox. Our results demonstrate that knowledge of the genomic N transformation potential is important for interpreting N process rates and understanding how the lacustrine sedimentary N cycle responds to variations in trophic conditions.

IMPORTANCE Anthropogenic nitrogen (N) inputs can lead to eutrophication in surface waters, especially in N-limited coastal ecosystems. Lakes effectively remove reactive N by transforming it to N₂ through microbial denitrification or anammox. The rates and distributions of these microbial processes are affected by factors such as the amount and quality of settling organic material and nitrate concentrations. However, the microbial communities mediating these N transformation processes in freshwater lake sediments remain largely unknown. We provide the first seasonally and spatially resolved metagenomic analysis of the N cycle in sediments of two lakes with different trophic states. We show that lakes with different trophic states select for distinct communities of N-cycling microorganisms with contrasting functional potentials for N transformation.

Undervalued Pseudo-nifH Sequences in Public Databases Distort Metagenomic Insights into Biological Nitrogen Fixers

Nitrogen fixation, a distinct process incorporating the inactive atmospheric nitrogen into the active biological processes, has been a major topic in biological and geochemical studies. Currently, insights into diversity and distribution of nitrogen-fixing microbes are dependent upon homology-based analyses of nitrogenase genes, especially the nifH gene, which are broadly conserved in nitrogen-fixing microbes. Here, we report the pitfall of using nifH as a marker of microbial nitrogen fixation. We exhaustively analyzed genomes in RefSeq (231,908 genomes) and KEGG (6,509 genomes) and cooccurrence and gene order patterns of nitrogenase genes (including nifH) therein. Up to 20% of nifH-harboring genomes lacked nifD and nifK, which encode essential subunits of nitrogenase, within 10 coding sequences upstream or downstream of nifH or on the same genome. According to a phenotypic database of prokaryotes, no species and strains harboring only nifH possess nitrogen-fixing activities, which shows that these nifH genes are "pseudo"-nifH genes. Pseudo-nifH sequences mainly belong to anaerobic microbes, including members of the class Clostridia and methanogens. We also detected many pseudo-nifH reads from metagenomic sequences of anaerobic environments such as animal guts, wastewater, paddy soils, and sediments. In some samples, pseudo-nifH overwhelmed the number of "true" nifH reads by 50% or 10 times. Because of the high sequence similarity between pseudo- and true-nifH, pronounced amounts of nifH-like reads were not confidently classified. Overall, our results encourage reconsideration of the conventional use of nifH for detecting nitrogen-fixing microbes, while suggesting that nifD or nifK would be a more reliable marker. IMPORTANCE Nitrogen-fixing microbes affect biogeochemical cycling, agricultural productivity, and microbial ecosystems, and their distributions have been investigated intensively using genomic and metagenomic sequencing. Currently, insights into nitrogen fixers in the environment have been acquired by homology searches against nitrogenase genes, particularly the nifH gene, in public databases. Here, we report that public databases include a significant amount of incorrectly annotated nifH sequences (pseudo-nifH). We exhaustively investigated the genomic structures of nifH-harboring genomes and found hundreds of pseudo-nifH sequences in RefSeq and KEGG. Over half of these pseudo-nifH sequences belonged to members of the class Clostridia, which is supposed to be a prominent nitrogen-fixing clade. We also found that the abundance of nitrogen fixers in metagenomes could be overestimated by 1.5 to >10 times due to pseudo-nifH recorded in public databases. Our results encourage reconsideration of the prevalent use of nifH as a marker of nitrogen-fixing microbes.

Microbial ecology of sulfur cycling near the sulfate-methane transition of deep-sea cold seep sediments

Microbial sulfate reduction is largely associated with anaerobic methane oxidation and alkane degradation in sulfate-methane transition zone (SMTZ) of deep-sea cold seeps. How the sulfur cycling is mediated by microbes near SMTZ has not been fully understood. In this study, we detected a shallow SMTZ in three of eight sediment cores sampled from two cold seep areas in the South China Sea. One hundred ten genomes representing sulfur-oxidizing bacteria (SOB) and sulfur-reducing bacteria (SRB) strains were identified from three SMTZ-bearing cores. In the layers above SMTZ, SOB were mostly constituted by Campylobacterota, Gammaproteobacteria and Alphaproteobacteria that probably depended on nitrogen oxides and/or oxygen for oxidation of sulfide and thiosulfate in near-surface sediment layers. In the layers below the SMTZ, the deltaproteobacterial SRB genomes and metatranscriptomes revealed CO₂ fixation by Wood-Ljungdahl pathway, sulfate reduction and nitrogen fixation for syntrophic or fermentative lifestyle. A total of 68% of the metagenome assembled genomes were not adjacent to known species in a phylogenomic tree, indicating a high diversity of bacteria involved in sulfur cycling. With the large number of genomes for SOB and SRB, our study uncovers the microbial populations that potentially mediate sulfur metabolism and associated carbon and nitrogen cycles, which sheds light on complex biogeochemical processes in deep-sea environments.

Deep sea sediments associated with cold seeps are a subsurface reservoir of viral diversity

In marine ecosystems, viruses exert control on the composition and metabolism of microbial communities, influencing overall biogeochemical cycling. Deep sea sediments associated with cold seeps are known to host taxonomically diverse microbial communities, but little is known about viruses infecting these microorganisms. Here, we probed metagenomes from seven geographically diverse cold seeps across global oceans to assess viral diversity, virus-host interaction, and virus-encoded auxiliary metabolic genes (AMGs). Gene-sharing network comparisons with viruses inhabiting other ecosystems reveal that cold seep sediments harbour considerable unexplored viral diversity. Most cold seep viruses display high degrees of endemism with seep fluid flux being one of the main drivers of viral community composition. In silico predictions linked 14.2% of the viruses to microbial host populations with many belonging to poorly understood candidate bacterial and archaeal phyla. Lysis was predicted to be a predominant viral lifestyle based on lineage-specific virus/host abundance ratios. Metabolic predictions of prokaryotic host genomes and viral AMGs suggest that viruses influence microbial hydrocarbon biodegradation at cold seeps, as well as other carbon, sulfur and nitrogen cycling via virus-induced mortality and/or metabolic augmentation. Overall, these findings reveal the global diversity and biogeography of cold seep viruses and indicate how viruses may manipulate seep microbial ecology and biogeochemistry.

Diversity of microbial community and its metabolic potential for nitrogen and sulfur cycling in sediments of Phu Quoc island, Gulf of Thailand

Although Phu Quoc island, Gulf of Thailand possesses diverse marine and coastal ecosystems, biodiversity and metabolic capability of microbial communities remain poorly investigated. The aim of our study was to evaluate the biodiversity and metabolic potential of sediment microbial communities in Phu Quoc island. The marine sediments were collected from three different areas and analyzed by using 16S rRNA gene-based amplicon approach. A total of 1,143,939 reads were clustered at a 97% sequence similarity into 8,331 unique operational taxonomic units, representing 52 phyla. Bacteria and archaea occupied averagely around 86% and 14%, respectively, of the total prokaryotic community. Proteobacteria, Planctomycetes, Chloroflexi, and Thaumarchaeota were the dominant phyla in all sediments, which were involved in nitrogen and sulfur metabolism. Sediments harboring of higher nitrogen sources were found to coincide with increased abundance of archaeal phylum Thaumarchaeota. Predictive functional analysis showed high abundance prokaryotic genes associated with nitrogen cycling including nifA-Z, amoABC, nirA, narBIJ, napA, nxrAB, nrfA-K, nirBD, nirS, nirK, norB-Z, nlnA, ald, and ureA-J, based on taxonomic groups detected by 16S rRNA sequencing. Although the key genes involved in sulfur cycling were found to be at low to undetectable levels, the other genes encoding for sulfur-related biological processes were present, suggesting that alternative pathways may be involved in sulfur cycling at our study site. In conclusion, our study for the first time shed light on diversity of microbial communities in Phu Quoc island.