Next Article in Journal
Gut Microbial Adaptation to Varied Altitudes and Temperatures in Tibetan Plateau Yaks
Previous Article in Journal
Structure and Diversity of Endophytic Bacteria in Maize Seeds and Germinating Roots
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity, Methane Oxidation Activity, and Metabolic Potential of Microbial Communities in Terrestrial Mud Volcanos of the Taman Peninsula

by
Alexander I. Slobodkin
1,*,
Igor I. Rusanov
1,
Galina B. Slobodkina
1,
Aleksandra R. Stroeva
2,
Nikolay A. Chernyh
1,
Nikolai V. Pimenov
1 and
Alexander Y. Merkel
1
1
Winogradsky Institute of Microbiology, Research Center of Biotechnology, Russian Academy of Sciences, Leninskiy Prospect, 33, Bld. 2, 119071 Moscow, Russia
2
Lomonosov Moscow State University, 119234 Moscow, Russia
*
Author to whom correspondence should be addressed.
Microorganisms 2024, 12(7), 1349; https://doi.org/10.3390/microorganisms12071349
Submission received: 29 May 2024 / Revised: 21 June 2024 / Accepted: 26 June 2024 / Published: 1 July 2024
(This article belongs to the Section Environmental Microbiology)

Abstract

:
Microbial communities of terrestrial mud volcanoes are involved in aerobic and anaerobic methane oxidation, but the biological mechanisms of these processes are still understudied. We have investigated the taxonomic composition, rates of methane oxidation, and metabolic potential of microbial communities in five mud volcanoes of the Taman Peninsula, Russia. Methane oxidation rates measured by the radiotracer technique varied from 2.0 to 460 nmol CH4 cm−3 day−1 in different mud samples. This is the first measurement of high activity of microbial methane oxidation in terrestrial mud volcanos. 16S rRNA gene amplicon sequencing has shown that Bacteria accounted for 65–99% of prokaryotic diversity in all samples. The most abundant phyla were Pseudomonadota, Desulfobacterota, and Halobacterota. A total of 32 prokaryotic genera, which include methanotrophs, sulfur or iron reducers, and facultative anaerobes with broad metabolic capabilities, were detected in relative abundance >5%. The most highly represented genus of aerobic methanotrophs was Methyloprofundus reaching 36%. The most numerous group of anaerobic methanotrophs was ANME-2a-b (Ca. Methanocomedenaceae), identified in 60% of the samples and attaining relative abundance of 54%. The analysis of the metagenome-assembled genomes of a community with high methane oxidation rate indicates the importance of CO2 fixation, Fe(III) and nitrate reduction, and sulfide oxidation. This study expands current knowledge on the occurrence, distribution, and activity of microorganisms associated with methane cycle in terrestrial mud volcanoes.

1. Introduction

Mud volcanism is a geological phenomenon of great importance for hydrocarbon exploration, seismicity, and the atmospheric budget of methane [1]. Terrestrial mud volcanoes (TMVs) are the structures located on the Earth’s surface but geologically connected to deep hydrocarbon-bearing sedimentary basins and, thus, to the subsurface biosphere [2]. Morphological structures of TMVs are diverse and include mud pools, salsa lakes, and gryphons. The breccia, liquid, and gas emissions from TMVs contain different inorganic and organic compounds, which can be used as electron donors and acceptors in microbial metabolism. Unlike marine mud volcanoes, where the main electron acceptor is sulfate, the main terminal electron acceptors in TMVs are probably iron minerals [3,4].
Microbial communities inhabiting TMVs have been studied by cultivation and molecular approaches. Aerobic and anaerobic microorganisms were isolated from TMVs located in different geographical regions [5,6,7,8,9,10,11]. TMVs in Azerbaijan, China, Iran, Italy, Japan, Romania, Russia, Taiwan, and Trinidad were inspected in relation to microbial community composition by 16S rRNA gene amplicon sequencing [5,12,13,14,15,16,17,18,19,20,21,22,23,24]. Metagenomic analyses of the genomes and functional genes were also applied to explore the diversity and metabolic potential of prokaryotes thriving in TMVs [3,25,26]. Overall, the results obtained point to an active methane, sulfur, and iron cycling in the studied ecosystems. However, most of these studies focus on a single environmental site, and comparisons of multiple TMVs located in a given geographic region are rare.
Methane oxidation (MO) in mud volcanoes is performed by aerobic methanotrophic bacteria and anaerobic methanotrophic archaea (ANME). Molecular surveys show the presence in TMVs of aerobic methylotrophs of the genera Methylobacter, Methylomicrobium, Methylomonas, Methylophaga, Methyloprofundus, and Methylosoma [4,12,23,26]. Various groups of ANME (ANME-1, ANME-2a-b, Ca. Methanoperedenaceae, ANME-3) have been detected in TMVs worldwide [23,24,26].
Anaerobic methane oxidation (AOM) is a globally important process that prevents the release of a significant portion of methane into the atmosphere [27]. The biological mechanisms of AOM are still insufficiently studied, but it is generally assumed that syntrophic interactions between ANME and partner bacteria play a pivotal role. ANME can couple methane oxidation with the reduction of terminal electron acceptors either alone, or in concert with a bacterial partner via extracellular electron transfer. AOM dependent on the reduction of sulfate, nitrate, nitrite, Fe(III), and Mn(IV) compounds has been described [28,29,30,31,32,33,34]. The metabolism of bacteria, which participate in AOM, is poorly characterized; most of these microorganisms are still not culturable. In sulfate-reducing AOM consortia, ANME are usually associated with bacteria belonging to different clades of the phylum Desulfobacterota [35,36]. Genomes of these bacteria always contain numerous multiheme c-type cytochromes [37,38,39]. The involvement of Desulfobacterota in Fe(III)-dependent AOM in TMVs has also been proposed [4].
The rate of AOM is a significant factor that determines methane emissions into the atmosphere. To quantify AOM, environmental samples are incubated with radioactive tracers, such as 14CH4. Values of AOM rates vary from 0.0001 nmol cm−3 day−1 for lake sediments to 10,000 nmol cm−3 day−1 for Black Sea microbial reefs. Most marine mud volcanoes exhibit high rates of microbial methane oxidation (10–1500 nmol cm−3 day−1) compared to other methane-containing environments [27,40]. However, TMVs are markedly different from marine mud volcanoes due to the absence of sulfate as a potential electron acceptor for AOM. In situ MO activity in TMVs has never been measured.
The Kerch–Taman mud volcanic province is one of the most active areas of mud volcanism in Eurasia. More than 50 TMVs are located on the Taman peninsula (the Krasnodar Krai, Russia). The total thickness of the sedimentary strata of the peninsula is 10 km; the thickness of the Oligocene Neogene deposits is 7 km, with the Maikop series playing a dominant role. Hydrochemical (Li–Na and Li–Mg) geothermometers suggest that TMV waters in this area formed at depths ranging from 1 to 4.5 km [41].
For a complete understanding of the functioning of a complex microbial community, the most productive approach is a combination of two approaches: a cultivation approach based on the isolation of pure cultures and a cultivation-independent approach based on the use of molecular methods. Several new taxa of anaerobic bacteria were isolated and described from Taman TMVs [42,43,44,45,46]. Physiological studies of these microorganisms suggested the functioning of sulfur and iron biogeochemical cycles, but none of these isolates were highly represented in the community based on 16S rRNA gene profiling or metagenomic data. On the other hand, microbiological studies of the Taman TMVs via molecular methods have only been carried out for two microbial communities [4,26], and it is unclear how these limited data reflect the common features of the TMVs of peninsula.
The main objective of this study was to obtain information on biological mechanisms of methane oxidation and determine the dominant taxonomic groups of methanotrophic bacteria and archaea in Taman TMVs. In this paper, we report the taxonomic composition and methane oxidation rates of 17 microbial communities from five TMVs of the Taman peninsula. We compared prokaryotic diversity and revealed the common taxa of the different taxonomic ranks present in all TMVs and identified the major methanotrophic players of aerobic and anaerobic methane oxidation at the genus level. We obtained 29 high-quality metagenome-assembled genomes (MAGs) from the community with high AOM rate and assessed their metabolic capabilities, allowing us to obtain conclusions on possible physiological groups of microorganisms and electron donors and acceptors involved in methane oxidation. We also for the first time established the ranges and maximal values of methane oxidation activities in terrestrial mud volcanoes by the radiotracer technique.

2. Materials and Methods

2.1. Sampling Site Description, Sample Collection, and Chemical Analysis

Samples were collected in June 2022, from the fields of five mud volcanoes located in the Taman peninsula, Krasnodar Krai, Russia. Samples were taken at a depth of 15–20 cm below the surface, placed into sterile 50 mL bottles, and transported to the laboratory where they were prepared for molecular and geochemical studies. A description of the samples is given in the Section 3.
Temperature and pH of the samples were measured at the point of collection using a digital pH meter equipped with a thermocouple (AZ Instrument, Moscow, Russia). The concentrations of chloride, sulfate, and nitrate were analyzed in water extracts by HPLC with a Stayer ion chromatograph (Aquilon, Moscow, Russia) with an IonPack AS4-ASC column (Dionex, Sunnyvale, CA, USA) and conductivity detector; the eluent was bicarbonate (1.36 mM)/carbonate (1.44 mM), and the flow rate was 1.5 mL min−1.

2.2. Methane Content and Rates of Microbial Methane Oxidation

The activity of processes of microbial oxidation of methane in the TMVs was determined by the radioisotope method. Methane concentrations were measured by the phase-equilibrium degassing method on a gas chromatograph, Kristall-2000-M (Chromatec, Yoshkar-Ola, Russia), equipped with a flame ionization detector. The measurement error did not exceed ±5%. The rate of methane oxidation (MO) was determined by the radioisotope method with 14C-methane ([14C]-methane, specific activity 1.16 GBq/mmol, JSC Isotope, Russia) dissolved in degassed distilled water. Undisturbed sediments in amounts of 2.0 cm3 were taken with plastic syringes (total volume 5 cm3) with a rubber piston and a cut edge and hermetically sealed with a butyl rubber stopper. For the variant with the addition of air, 2.0 cm3 of the gas phase was left. Then, 0.2 mL of labeled 14C-methane (1 µCi per sample) was added with a tuberculin syringe, piercing the rubber stopper with a needle in the center and evenly distributing the substrate along the length of the syringe with the sediment. Samples fixed with 1 mL of 0.5 N KOH solution before adding the labeled substrate served as controls. After adding labeled methane, sediment samples were incubated for 1 day at a temperature close to the in situ temperature. After incubation, the samples were fixed and transported to a stationary laboratory for further determination of the products of microbial oxidation and transformation of C-methane into carbon dioxide (CO2), into microbial biomass (BM), and into dissolved organic carbon (DOC) in dissolved organic matter (DOM). Sample processing and calculations of the intensity of methane oxidation were carried out according to a previously described method [47]. Radioactivity (14C) of the products of the microbial process of MO was measured on a TRI-CarbTR 2400 liquid scintillation counter (Packard, Conroe, TX, USA).

2.3. DNA Extraction and 16S rRNA Gene Amplicon and Metagenome Library Preparation, Sequencing, and Analysis

DNA from mud samples was isolated using the FastDNA Spin Kit for Soil according to the manufacturer’s protocol (MP Biomedicals, Santa Ana, CA, USA). Preparation, sequencing, and analysis of the V3-V4 region of the 16S rRNA amplicon libraries were performed as previously described [26]. qPCR-based total prokaryote quantification was performed as previously described [48].
Preparation and sequencing of a shotgun metagenome library were performed in Laboratory “Genomed” Ltd., Moscow, Russia, using the MGIEasy Fast PCR-FREE FS Library Prep Set (MGI, China) according to the manufacturer’s protocol and the DNBSEQ-G400 system (MGI, China) with the reagent kit, which can read 150 nucleotides from each end. Raw read processing, contig assembling, binning, bin refinement, coverage calculation, bin completeness and contamination evaluation were performed as previously described [4]. Taxonomic positions were assigned to each bin according to the Genome Taxonomy Database (GTDB) classification as described earlier [4].
Metagenome-assembled genomes (MAGs) were submitted for gene calling and annotations through the RAST [49]. Genes encoding metabolic functions were queried in RASTtk using Blast search. A list of functional genes queried in this study can be found in Supplementary Table S1. All the sequencing data are deposited in NCBI BioProject PRJNA940400. The MAGs Accession Numbers are JARGFH000000000—JARGGJ000000000.

3. Results and Discussion

3.1. Geochemical Characteristics of the Studied Sites

Seventeen samples of mud volcanic fluids and sediments were taken from mud pools, gryphons, and salsa lakes located on the fields of five mud volcanoes: Gnilaya Gora (45.251 N, 37.436 E), Gladkovsky (45.005 N, 37.724 E), Kuchugursky (45.432 N, 36.922 E), Semigorsky (44.901 N, 37.597 E), and Shugo (45.070 N, 37.611 E). The main geochemical parameters that control the functioning of microbial communities were determined: temperature, pH, Eh, salinity, methane content, concentrations of chloride, sulfate, nitrate, and phosphate (Table 1). All samples had a temperature close to ambient temperature (about 23 °C) and a neutral or alkaline pH (from 6.72 to 9.10). The salinity of the samples varied from 9 to 30 g L−1, Eh from +15 to −348 mV, methane content from 3 to 917 mM, chloride concentration from 5.7 to 412 mM, and sulfate concentration from 0 to 63 mM. None of the samples had significant concentrations of nitrate and phosphate (>1 mM).

3.2. Microbial Methane Oxidation

The potential in situ rates of microbial methane oxidation (MO) were measured directly in the original samples by the radiotracer technique (Figure 1). In order to evaluate the effect of oxygen on methane oxidation, the samples were incubated in parallel in the presence and absence of atmospheric air in the gas phase. Air oxygen is supposed to inhibit the anaerobic oxidation of methane and stimulate the aerobic oxidation. The experimental results showed high MO values, both in aerobic (maximum of 213 nmol CH4 cm−3 day−1, with an average value of 50 nmol CH4 cm−3 day−1) and anaerobic conditions, reaching values of 400–460 nmol CH4 cm−3 day−1 (average value of 150 nmol CH4 cm−3 day−1). These values are several times higher than those of MO activity in the surface sediments of the underwater mud volcano Haakon Mosby (72 degrees N) and in the gas seepage fields of the Vestnesa Ridge, where the maximum values of microbial methane oxidation reached 70 nmol CH4 cm−3 day−1 [50]. The highest MO rates were observed in the samples from the TMV Gladkovsky. Moreover, the MO rates under anaerobic conditions significantly (in most cases by several times) exceeded the MO rates in the presence of air. Under aerobic conditions, methane consumption ranged from 2 to 39% of the total methane injected into the sample, whereas under anaerobic conditions, 10 to 98% of the total methane injected into the sample was consumed. These results indicate the activity of methanotrophic anaerobic microorganisms in the studied TMVs.
During methane oxidation, the carbon from methane is transformed not only into carbon dioxide but is also included in biomass as well as in various low-molecular-weight extracellular dissolved organic compounds (DOC), which can serve as a carbon source for heterotrophic microorganisms. In our experiments, the average proportion of the carbon from methane transformed into carbon dioxide was about 64% for all investigated TMVs (Figure 1). The carbon from methane in DOC was on average 35%, and only 1–2% of the methane carbon was assimilated into the biomass. These data indicate the important role of exometabolites in methane oxidation.

3.3. Microbial Community Abundance and Composition

Microbial communities of the studied sites differed significantly in the prokaryote total number estimated by 16S rRNA gene-based qPCR. The highest values of total prokaryote abundance were found in site C20 of the Semigorsky volcano (4.6 × 109 16S rRNA gene copies (g.c.) calculated per 1 mL of mud pool sediment, standard deviation (SD) 7.3%), site C05 of the Gnilaya Gora volcano (2.9 × 109 16S rRNA g.c. mL−1 of mud pool sediment, SD 12.3%), and site C14 of the Gladkovsky volcano (2.9 × 109 16S rRNA g.c. mL−1 of mud pool sediment, SD 9.9%). The lowest density of microbial communities was found in site C02 of the Gnilaya Gora volcano (1.2 × 105 16S rRNA g.c. mL−1 of gryphon mud, SD 2.4%), site C10 of the Gladkovsky volcano (3.6 × 105 16S rRNA g.c. mL−1 of mud pool sediment, SD 5.3%), and site C24 of the Shugo volcano (3.8 × 105 16S rRNA g.c. mL−1 of gryphon mud, SD 8.1%). All values are shown in Supplementary Table S2.
Using high-throughput sequencing of 16S rRNA gene amplicons, the phylogenetic composition of 17 microbial communities from five TMVs was determined. A total of 137,673 sequences of the 16S rRNA gene V3-V4 region were taken for analysis after all the initial steps of data processing. These sequences were merged into 967 amplicon sequence variants (ASVs). Alpha diversity metrics such as Chao1 richness (the predicted total number of phylotypes in a community) and the Shannon index (index of the complexity of the internal structure of a community) were calculated after normalization of the data via rarefying. The highest and the lowest values of these indexes were found in sites of the Gladkovsky mud volcano. Sites C12 and C13 had Chao1 richnesses of 121 and 142 and Shannon indices of 4.15 and 4.28, respectively, whereas sites C10 and C14 had Chao1 richnesses of 41 and 60 and Shannon indices of 2.25 and 2.59, respectively. All values are shown in Supplementary Table S3.
Bacteria dominated almost all the studied samples, accounting for 65–99%, with the exception of one sample (C11), where archaea of the phylum Halobacterota constituted 56% of all the prokaryotic sequences obtained. Twenty-one phyla (Silva and GTDB taxonomy) of Bacteria were found in relative abundance of more than 1% in at least one sample (Figure 2). Eight phyla made up a significant proportion of the community in most samples: Pseudomonadota (11–82%), Halobacterota (1–56%) Desulfobacterota (2–51%), Cyanobacteriota (1–30%), Actinomycetota (1–21%), Chloroflexota (1–20%), Bacillota (1–15%), and Bacteroidota (1–7%). In general, these data are consistent with those obtained for other TMV microbial communities on the Eurasian continent, suggesting that, due to the lack of fluid exchange, TMVs represent a patchy habitat populated by microorganisms highly adapted to local environmental conditions [20]. Taman TMVs also contain a high percentage of Halobacterota, Desulfobacterota, and Pseudomonadota, which is typical for seafloor methane seeps [51].
Thirty-one prokaryotic genera were detected in relative abundance of more than 5% in at least one sample (Figure 3). Based on the metabolic capabilities of the cultivated representative of these genera, the identified microorganisms can be divided into four physiological groups: (1) methanotrophs; (2) sulfur or iron reducers; (3) aerobes or facultative anaerobes; (4) bacteria with an unknown type of metabolism.
Aerobic methanotrophic bacteria were present in all the samples. Methyloprofundus was the most highly represented genus of methanotrophs found in all studied TMVs, with relative abundance reaching 36%. Methyloprofundus spp. have been detected by molecular methods in many sites all around the world, mainly in marine sediments, while the cultivated members of the genus are limited to a single species [52,53]. Methanotrophic aerobes of the genus Methylomicrobium also constituted a significant part of some studied microbial communities.
Anaerobic methanotrophic archaea (ANME) belonging to different phylogenetic groups [54] were detected in 15 out of 17 samples (Figure 4). The relative abundance of ANME varies from 0 to 56.2% in all prokaryotic sequences. The most widely represented and numerous group was ANME-2a-b (Ca. Methanocomedenaceae), identified in 10 samples. Other detected ANME groups were present in smaller proportions and included the following: ANME-3 (Ca. Methanovorans), ANME-2d (Ca. Methanoperedenaceae), ANME-1a (Ca. Methanophagaceae), and ANME-2c (Ca. Methanogasteraceae). The highest relative abundance of ANME was found in the mud pools of the TMV Gladkovsky (up to 54% ANME-2a-b).
The ability to dissimilate sulfur or iron compounds is characteristic for microorganisms of the phylum Desulfobacterota (Figure 3). Microbial processes performed by these bacteria can include sulfate reduction (Desulfocapsa, unclassified Desulfobulbaceae), elemental sulfur reduction (Desulfuromusa, Desulfurivibrio, Desulfocapsa, unclassified Desulfuromonadia, unclassified Dethiobacteraceae), sulfur disproportionation (Desulfurivibrio, Desulfocapsa). The majority of sulfur-reducing and sulfur-disproportionating genera were detected in the samples from the TMVs Gnilaya Gora and Shugo. Bacteria capable of Fe(III) reduction included Geothermobacter, unclassified Desulfobulbaceae, and Desulfuromusa. Recently, the capacity for iron reduction has been demonstrated for OPB41 bacteria, which were described as a new order, Anaerosomatales [55]. Also, the uncultured Sva1033 (Desulfuromonadales) bacterial group could be involved in Fe(III) reduction [56]. Possibly, iron and sulfur reducers are syntrophic partners of ANME and complete electron transfer to the terminal acceptor.
Several identified cosmopolitan genera of Gammaproteobacteria with high relative abundance are known as microaerophiles or facultative anaerobes with broad metabolic capabilities. For example, except organoheterotrophic lifestyle, Halomonas spp. have the abilities for denitrification and heterotrophic sulfur oxidation, and some species of Roseovarius can grow lithoautotrophically utilizing sulfur compounds, molecular hydrogen, and nitrate [45,57]. The genus Mariprofundus is characterized by obligate lithoautotrophic microaerobic oxidation of Fe(II) [58].

3.4. MAGs’ General Characteristics and Phylogenetic Identification

Based on a high microbial population density (2.9 × 109 16S rRNA g.c. mL−1 of mud pool sediment), high ANME-2a-2b (Ca. Methanocomedenaceae) content, and high rate of methane oxidation, site C14 of the Gladkovsky volcano was chosen for further metagenomic study. Shotgun sequencing resulted in a ~12 Gb that was assembled in 30,910 contigs with a total length of 147 Mb (N50—9401b).
Twenty-nine MAGs with >70% completeness and <4.0% contamination were recovered from the metagenome. The MAGs statistics, including the number of contigs, genome size, the presence of 16S rRNA gene, and relative abundance calculated on the basis on coverage estimation, are shown in Table 2.
According to a phylogenomic analysis based on GTDB 08-RS214 [59], 25 of 29 MAGs were assigned to the domain Bacteria and 4 were assigned to the domain Archaea. Almost all the MAGs belonged to currently uncultivated microorganisms of different taxonomic ranks. MAG C14-19, the second most abundant in metagenome, and MAG C14-15 represent novel phyla in the domain Bacteria. Three MAGs (C14-28, C14-29, and C14-12) belong to new bacterial classes: three to new orders; nine to new families; and seven to new genera. Only five MAGs could be assigned to the described genera (Methyloprofundus, Brevefilum, Ca. Methanovorans, Lutibacter, and Sedimenticola). A complete (>1400 bp) 16S rRNA gene was recovered in 7 of 29 MAGs. The most abundant taxa present in the metagenome were members of uncultivated c_Acidimicrobiia (19%), unclassified d_Bacteria (18%), two populations of anaerobic methanotrophic archaea of f_Ca. Methanocomedenaceae (14 and 12%), and aerobic methanotrophs of the genus Methyloprofundus (9%). The other 24 MAGs had relatively low abundance (0.3–3%).
Anaerobic methanotrophic archaea were represented by three MAGs. C14-17 and C14-08, belonging to Ca. Methanocomedenaceae (ANME-2a-2b), together accounted for 26% of the community. This is quite close to the assessment obtained when analyzing the V3-V4 16S rRNA gene profiles (27.9%). Phylogenetically, they appeared to be very close to the previously described MAGs from the methanotrophic enrichment culture from the Gladkovsky mud volcano [4]: ANI between C14-17 and B4-04 (GCA_023544575.1)—99.4% and between C14-08 and B4-03 (GCA_023544605.1)—96.2% (Supplementary Figure S1). The third MAG (C14-18) belonged to the Ca. Methanovorans (ANME-3), accounted for 0.4% of the community, and was also very close to the minor MAG from the same enrichment culture: ANI between C14-18 and B4-05 (GCA_023544565.1)—99.5% (Supplementary Figure S1).

3.5. Metabolic Capabilities in MAGs

The presence of the functional genes related to energy metabolism and carbon and nitrogen fixation in the MAGs is shown in Figure 5.
Primary production of organic matter in this microbial community may be partly due to methane oxidation. The genomes of methanotrophic archaea (Ca. Methanocomedenaceae, Ca. Methanovorans), which account for 26% of the prokaryotic sequences, contain the gene of methyl coenzyme M reductase (mcrA), responsible for the anaerobic oxidation of methane [60]. Genes of particulate methane monooxygenase (pMMO) as well as soluble methane monooxygenase (sMMO), performing the initial step of aerobic methane oxidation, are present in the genome of Methyloprofundus, which is also abundant in the studied community. Methanol, likely formed as an intermediate of methane oxidation abiotically or by methanotrophs [61,62], may be used aerobically by Methylophagaceae or anaerobically by Spirochaetaceae and Anaerolineae, whose genomes contain genes of methanol dehydrogenase (mdh) or methanol:corrinoid methyltransferase (mtaB). To date, the ability to utilize methanol has been found in only one bacterium isolated from a TMV: Pelomicrobium methylotrophicum [44].
The genes of the key enzymes involved in autotrophic CO2 fixation were identified in almost all MAGs, excluding two MAGs belonging to Paceibacteria and Patescibacteria. This does not mean that the majority of microorganisms in this community are autotrophs, but it does indicate the possible distribution of the pathways of carbon fixation in the studied ecosystem [63]. The genes of ribulose bisphosphate carboxylase form I (RuBisCO, rbc1), the main enzyme of the Calvin−Benson−Bassham cycle, were found in the MAGs of Methylophagaceae and Sedimenticola thiotaurini. The members of the genus Sedimenticola are facultative lithoautotrophs, and recently a new species—S. hydrogenitrophicus, capable of lithoautotrophic growth with molecular hydrogen or thiosulfate as electron donors and nitrate, nitrous oxide, or oxygen as electron acceptors—was isolated from the studied TMV Gladkovsky [46]. CO dehydrogenase/acetyl-CoA synthase (acs/cooS), the key catalytic complex of the Wood−Lyngdahl pathway (WL) was encoded in 8 out of 29 MAGs including all Archaea and the majority of Desulfobacterota. Similarly, WL predominates over other carbon assimilation pathways in AOM microbial communities of the Karabetova Gora mud volcano, also located on the Taman peninsula [26]. Marker genes of other CO2 fixation pathways were less common. ATP-citrate lyase (acl) and citryl-CoA lyase (ccl) (two variants of the reductive tricarboxylic acid cycle) were encoded in unknown Bacteria (MAG C14-19), Desulfobacterales, and Zixibacteria. The genes of 4-hydroxybutyryl-CoA dehydratase (4-hbd) (3-hydroxypropionate/4-hydroxybutyrate and dicarboxylate/4-hydroxybutyrate cycles) were present in Acidimicrobiia, Desulfobacterales, and Thorarchaeaceae. The genes of malonyl-CoA reductase (mlr, 3-hydroxypropionate bi-cycle) were not found in any MAG. Two carbon fixation pathways—the reversed oxidative tricarboxylic acid cycle and reductive glycine pathway—cannot be identified based on genomic data [64,65,66]. However, the genes of citrate synthase (cs) and glycine dehydrogenase (dgh) required for these assimilation processes were present in most of the MAGs.
Half of the MAGs contained genes of the oxygen reduction enzymatic systems: cytochrome c oxidase of aa3, cbb3, or ba3 types or cytochrome d ubiquinol oxidase (coxA, ccoN, ba3, and cydA genes, respectively). These enzymes can be used for either respiration or detoxification [67], but it is not possible to distinguish between these processes based on genomic information alone. A similar situation was observed in another Taman TMV, Karabetova Gora, where 35% of the microbial community members were potential aerobes [26]. Most probably, oxygen penetrated from the atmosphere is constantly present in the upper layers of the mud. However, the MAGs of all detected Archaea as well as of unknown Bacteria, Paceibacteria, Brevefilum, Patescibacteria, Anaerosomatales, Spirochaetaceae, and Anaerolineae did not encode any cytochrome oxidases indicating the strictly anaerobic lifestyle of these microorganisms.
Genes of uptake hydrogenases (hyaB, hybC) were present in 37% of the MAGs, indicating the importance of hydrogen metabolism in the community. One of Methanocomedenaceae, two unknown Bacteria (C14-19 and C14-15), Methylophagaceae, Bacteroidales, Sedimenticola, and Lutibacter have the potential to utilize molecular hydrogen. It should be noted that the detection of genes encoding [NiFe]-hydrogenase in ANME-2 (Methanocomedenaceae) suggests the ability for hydrogenotrophy, which is still unknown for this group of methanotrophic archaea. Previously, uptake hydrogenases genes were found in ANME-2 only in one study: in MAGs recovered from metagenome also associated with TMVs located on the Taman peninsula [4].
Genes of periplasmic Fe-Fe hydrogenases (fe-hyd) or GH1 betaglucosidase (gh1) indicative of fermentation were present in 30% of the MAGs, including the most abundant Acidimicrobia. A gene of the large chain of aerobic carbon monoxide dehydrogenase (coxL) was found in the MAGs of Acidimicrobia and Anaerolineaceae. Genes encoding the catalytic subunit of anaerobic CO-dehydrogenase (cooS) were present in five MAGs.
Anaerobic methane oxidation could be coupled to nitrate, sulfate, or Fe(III) reduction. Genes of the catalytic subunit of the respiratory nitrate reductase of the Nar type (narG) were identified in Acidimicrobiia and Methyloprofundus which are present in the community in high abundance. At present, the only known species of Methyloprofundus is M sedimenti. It was originally described as an aerobic methanotroph utilizing nitrate as a nitrogen source, and anaerobic growth with nitrate reduction was not tested [52]. Our genome analysis revealed that in addition to nasAB involved in nitrate assimilation, the genome of M sedimenti possesses narGHI, nirK, nirS, and norBC genes and does not contain nosZ. The class Acidimicrobiia has more than 4000 members according to GTDB. The order to which MAG C14-01 was assigned includes only uncultivated bacteria; so, their physiology and metabolism are unknown. Periplasmic nitrate reductase of the Nap type (napA) was encoded in Desulfurivibrionaceae and Sedimenticola. The capacity to grow with nitrate reduction is known for the representatives of these taxa [68,69]. The terminal enzyme of the complete denitrification, nitrous oxide reductase (nosZ), was found in five MAGs. Only MAG C14-22 (Sedimenticola) had both nitrate and nitrous oxide reduction genes.
The oxidative branch of the sulfur cycle was prominent in the MAGs. The key enzyme for sulfide oxidation, sulfide:quinone oxidoreductase (SQR), appears to be a ubiquitous protein present in both chemotrophic and phototrophic prokaryotes as well as in some mitochondria [70]. The presence of sqr was found in 8 of 29 MAGs (Methyloprofundus, Bacteroidales, Methylophagaceae, Anaerolineales, Trueperaceae, Bacteroidales, Sedimenticola, and Lutibacter). Additionally, the genes of sulfite dehydrogenase (soeA) and proteins of the Sox system (soxB, soxY), involved in sulfite and thiosulfate oxidation, were identified in the MAGs of Methylophagaceae and Sedimenticola. The family Methylophagaceae consists of the single genus Methylophaga, some species of which are capable of chemolithoheterotrophic growth with sulfide or thiosulfate oxidation [71,72]. All known Sedimenticola species are facultative-sulfur-compound-oxidizing lithoautotrophs [46]. The marker genes of sulfate reduction coding for alpha and beta subunits of dissimilatory sulfite reductase (dsrAB) and adenylylsulfate reductase (aprA) were detected in three MAGs with low abundance (Desulfobacterales, Desulfurivibrionaceae, and Sedimenticola). Thiosulfate/polysulfide reductase (psrA) was encoded only in Desulfobacterales.
The biochemical mechanisms of dissimilatory Fe(III) reduction are diverse, and currently, there are no universal genomic determinants of this process. However, it is well established that c-type cytochromes play a key role in electron transfer to the extracellular electron acceptor in many microorganisms [73]. The search for multiheme c-type cytochromes in the C14 metagenome revealed that 6 out of 29 MAGs contain genes of the proteins homologous to decaheme cytochrome c MtrA (mtrA), involved in iron reduction in Shewanella oneidensis [74]. Hence, it can be hypothesized that unknown Bacteria (MAG C14-19 and C14-15), Desulfuromonadales, Myxococcia, Desulfurivibrionaceae, and Zixibacteria participate in the reduction of Fe(III)-containing minerals in the studied ecosystem.
Homologs of respiratory periplasmic arsenate reductase (arrAB) were encoded in the MAG of Ca. Methanocomedenaceae C14-17. It is the second piece of evidence of the presence of Arr in ANME-2; previously, this enzyme was detected in methanotrophic Fe(III)-reducing enrichment culture obtained from Taman TMVs [4]. As(V) reduction has been also experimentally shown for a number of cultivated bacteria isolated from Taman TMVs [9,42,43,75], but other MAGs of C14 metagenome did not contain arrAB. The presence of nitrogenase genes (nifHD) in Methylophagaceae and Bacteroidales indicated the ability to fix dinitrogen.
In summary, five the most abundant members of the community (70% of the total sequences) contained genes indicative of aerobic and anaerobic methanotrophy, CO2 fixation, fermentation, Fe(III) and nitrate reduction, and sulfide oxidation. It can be assumed that two populations of Ca. Methanocomedenaceae couple AOM to Fe(III) reduction either directly or via syntrophic bacterial partner, which could be the unknown Bacteria C14-19; Methyloprofundus oxidizes methane aerobically in oxygenated zones of TMVs or links methane oxidation to nitrate reduction; facultative anaerobic Acidimicrobiia decomposes the produced organic matter.

4. Conclusions

Microbial communities of five mud volcanoes located at the Taman peninsula were investigated by molecular and radiotracer methods for the first time. Methane oxidation rates varied from 2.0 to 460 nmol CH4 cm−3 day−1 in different mud samples. MO rates under anaerobic conditions significantly exceeded the MO rates in the presence of air, indicating the activity of methanotrophic anaerobic microorganisms. Bacteria accounted for 65–99% of prokaryotic diversity in all the samples. The most abundant bacterial phyla were Pseudomonadota, Desulfobacterota, Bacteroidota, Actinomycetota, Bacillota, and Chloroflexota. Thirty-two prokaryotic genera, which include methanotrophs, sulfur or iron reducers, and facultative anaerobes with broad metabolic capabilities, were detected in relative abundance of more than 5% in at least one sample. The most highly represented genus of aerobic methanotrophs was Methyloprofundus. Its numbers were as high as 36% in the TMV Shugo. The most widely represented and numerous group of anaerobic methanotrophs was ANME-2a-b (Ca. Methanocomedenaceae), identified in 60% of samples and reaching 54% in the TMV Gladkovsky. Other detected ANME groups were present in smaller proportions and included the following: ANME-3 (Ca. Methanovorans), ANME-2d (Ca. Methanoperedenaceae), ANME-1a (Ca. Methanophagaceae), and ANME-2c (Ca. Methanogasteraceae). Twenty-nine high-quality MAGs were recovered from the metagenome of the sample with high rate of methane oxidation. Almost all the MAGs belonged to uncultivated microorganisms of different taxonomic ranks. The distribution of the functional genes between taxa points to the importance of methanotrophy, CO2 fixation, Fe(III) and nitrate reduction, and sulfide oxidation. Sulfate-reducing bacteria were not abundant, suggesting that in Taman TMVs, the biological mechanism of methane oxidation is not dependent on sulfate reduction but rather coupled to Fe(III) reduction. Overall, our study expands current knowledge on the occurrence, distribution, and activity of microorganisms associated with methane cycle in terrestrial mud volcanoes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms12071349/s1, Table S1: Functional genes; Table S2: qPCR; Table S3: Alfa diversity; Figure S1: Phylogenetic analysis of the ANME MAGs. Refs. [76,77,78,79] can be found in Supplementary Materials.

Author Contributions

Conceptualization, A.I.S. and N.V.P.; investigation, I.I.R., G.B.S., N.A.C., A.R.S., N.V.P. and A.Y.M.; data curation, A.Y.M.; writing—original draft preparation, A.I.S., I.I.R. and A.Y.M.; writing—review and editing, A.I.S.; funding acquisition, A.I.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation, project no. 22-14-00011 and by the Ministry of Science and Higher Education of the Russian Federation (MAGs annotation and radioisotope experiments).

Data Availability Statement

Data are contained within the article and supplementary materials. All the sequencing data are deposited in NCBI BioProject PRJNA940400. The MAGs Ac-cession Numbers are JARGFH000000000—JARGGJ000000000.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Mazzini, A.; Etiope, G. Mud volcanism: An updated review. Earth Sci. Rev. 2017, 168, 81–112. [Google Scholar] [CrossRef]
  2. Kopf, A.J. Significance of mud volcanism. Rev. Geophys. 2002, 40, 2-1–2-52. [Google Scholar] [CrossRef]
  3. Tu, T.-H.; Wu, L.-W.; Lin, Y.-S.; Imachi, H.; Lin, L.-H.; Wang, P.-L. Microbial community composition and functional capacity in a terrestrial ferruginous, sulfate-depleted mud volcano. Front. Microbiol. 2017, 8, 2137. [Google Scholar] [CrossRef] [PubMed]
  4. Slobodkin, A.I.; Ratnikova, N.M.; Slobodkina, G.B.; Klyukina, A.A.; Chernyh, N.A.; Merkel, A.Y. Composition and metabolic potential of Fe(III)-reducing enrichment cultures of methanotrophic ANME-2a archaea and associated bacteria. Microorganisms 2023, 11, 555. [Google Scholar] [CrossRef]
  5. Alain, K.; Holler, T.; Musat, F.; Elvert, M.; Treude, T.; Krüger, M. Microbiological investigation of methane- and hydrocarbon-discharging mud volcanoes in the Carpathian Mountains, Romania. Environ. Microbiol. 2006, 8, 574–590. [Google Scholar] [CrossRef] [PubMed]
  6. Arun, A.B.; Young, C.C.; Chen, W.M.; Hung, M.H.; Lai, W.A.; Chou, J.H.; Rekha, P.D.; Shen, F.T.; Su, S.P. Belliella pelovolcani sp. nov., isolated from a mud-volcano in Taiwan. Int. J. Syst. Evol. Microbiol. 2009, 59, 2534–2537. [Google Scholar] [PubMed]
  7. Ghiasian, M.; Akhavan Sepahy, A.; Amoozegar, M.; Saadatmand, S.; Shavandi, M. Bacterial diversity determination using culture-dependent and culture-independent methods. Glob. J. Environ. Sci. Manag. 2017, 3, 153–164. [Google Scholar] [CrossRef]
  8. Sepanian, E.; Sepahy, A.A.; Hosseini, F. Isolation and characterization of bacterial species from Ain mud volcano, Iran. Microbiology 2018, 87, 282–289. [Google Scholar] [CrossRef]
  9. Frolova, A.A.; Merkel, A.Y.; Kuchierskaya, A.A.; Bonch-Osmolovskaya, E.A.; Slobodkin, A.I. Pseudodesulfovibrio alkaliphilus, sp. nov., an alkaliphilic sulfate-reducing bacterium isolated from a terrestrial mud volcano. Antonie Leeuwenhoek 2021, 114, 1387–1397. [Google Scholar] [CrossRef]
  10. Ratnikova, N.M.; Slobodkin, A.I.; Merkel, A.Y.; Kopitsyn, D.S.; Kevbrin, V.V.; Bonch-Osmolovskaya, E.A.; Slobodkina, G.B. Sulfurimonas crateris sp. nov., a facultative anaerobic sulfur-oxidizing chemolithoautotrophic bacterium isolated from a terrestrial mud volcano. Int. J. Syst. Evol. Microbiol. 2020, 70, 487–492. [Google Scholar]
  11. Manna, S.K.; Das, B.K.; Mohanty, B.P.; Bandopadhyay, C.; Das, N.; Baitha, R.; Das, A.K. Exploration of heterotrophic bacterial diversity in sediments of the mud volcano in the Andaman and Nicobar Islands, India. Environ. Nanotechnol. Monit. Manag. 2021, 16, 100465. [Google Scholar]
  12. Wrede, C.; Brady, S.; Rockstroh, S.; Dreier, A.; Kokoschka, S.; Heinzelmann, S.M.; Heller, C.; Reitner, J.; Taviani, M.; Daniel, R.; et al. Aerobic and anaerobic methane oxidation in terrestrial mud volcanoes in the Northern Apennines. Sediment. Geol. 2012, 263, 210–219. [Google Scholar] [CrossRef]
  13. Green-Saxena, A.; Feyzullayev, A.; Hubert, R.J.; Kallmeyer, J.; Krueger, M.; Sauer, P.; Schulz, H.-M.; Orphan, V.J. Active sulfur cycling by diverse mesophilic and thermophilic microorganisms in terrestrial mud volcanoes of Azerbaijan. Environ. Microbiol. 2012, 14, 3271–3286. [Google Scholar] [CrossRef] [PubMed]
  14. Cheng, T.W.; Chang, Y.H.; Tang, S.L.; Tseng, C.H.; Chiang, P.W.; Chang, K.T.; Sun, C.H.; Chen, Y.G.; Kuo, H.C.; Wang, C.H.; et al. Metabolic stratification driven by surface and subsurface interactions in a terrestrial mud volcano. ISME J. 2012, 6, 2280–2290. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, P.-L.; Chiu, Y.-P.; Cheng, T.-W.; Chang, Y.-H.; Tu, W.-X.; Lin, L.-H. Spatial variations of community structures and methane cycling across a transect of Lei-Gong-Hou mud volcanoes in eastern Taiwan. Front. Microbiol. 2014, 5, 121. [Google Scholar] [CrossRef]
  16. Yang, H.-M.; Lou, K.; Sun, J.; Zhang, T.; Ma, X.-L. Prokaryotic diversity of an active mud volcano in the Usu City of Xinjiang, China. J. Basic Microbiol. 2012, 52, 79–85. [Google Scholar] [CrossRef]
  17. Lin, Y.T.; Tu, T.H.; Wei, C.L.; Rumble, D.; Lin, L.H.; Wang, P.L. Steep redox gradient and biogeochemical cycling driven by deeply sourced fluids and gases in a terrestrial mud volcano. FEMS Microbiol. Ecol. 2018, 94, fiy171. [Google Scholar] [CrossRef]
  18. Ren, G.; Ma, A.; Zhang, Y.; Deng, Y.; Zheng, G.; Zhuang, X.; Zhuang, G.; Fortin, D. Electron acceptors for anaerobic oxidation of methane drive microbial community structure and diversity in mud volcanoes. Environ. Microbiol. 2018, 20, 2370–2385. [Google Scholar] [CrossRef]
  19. Schulze-Makuch, D.; Haque, S.; Beckles, D.; Schmitt-Kopplin, P.; Harir, M.; Schneider, B.; Stumpp, C.; Wagner, D.A. A chemical and microbial characterization of selected mud volcanoes in Trinidad reveals pathogens introduced by surface water and rain water. Sci. Total Environ. 2020, 707, 136087. [Google Scholar] [CrossRef] [PubMed]
  20. Tu, T.-H.; Chen, L.-L.; Chiu, Y.-P.; Lin, L.-H.; Wu, L.-W.; Italiano, F.; Shyu, J.B.H.; Raisossadat, S.N.; Wang, P.-L. The biogeographic pattern of microbial communities inhabiting terrestrial mud volcanoes across the Eurasian continent. Biogeosciences 2022, 19, 831–843. [Google Scholar] [CrossRef]
  21. Megyes, M.; Móga, J.; Strat, D.; Borsodi, A.K. Bacterial and archaeal taxonomic diversity of mud volcanoes (Beciu, Romania) via metagenomic approach. Geomicrobiol. J. 2021, 38, 532–539. [Google Scholar] [CrossRef]
  22. Hsu, H.-C.; Chen, J.-S.; Nagarajan, V.; Hussain, B.; Huang, S.-W.; Rathod, J.; Hsu, B.-M. Assessment of temporal effects of a mud volcanic eruption on the bacterial community and their predicted metabolic functions in the mud volcanic sites of Niaosong, Southern Taiwan. Microorganisms 2021, 9, 2315. [Google Scholar] [CrossRef] [PubMed]
  23. Nagarajan, V.; Tsai, H.C.; Chen, J.S.; Koner, S.; Kumar, R.S.; Chao, H.C.; Hsu, B.M. Systematic assessment of mineral distribution and diversity of microbial communities and its interactions in the Taiwan subduction zone of mud volcanoes. Environ. Res. 2023, 216, 114536. [Google Scholar] [CrossRef] [PubMed]
  24. Miyake, N.; Ishimaru, R.; Komatsu, G.; Matsui, T. Characterization of archaeal and bacterial communities thriving in methane-seeping on-land mud volcanoes, Niigata, Japan. Int. Microbiol. 2023, 26, 191–204. [Google Scholar] [CrossRef] [PubMed]
  25. Mardanov, A.V.; Kadnikov, V.V.; Beletsky, A.V.; Ravin, N.V. Sulfur and methane-oxidizing microbial community in a terrestrial mud volcano revealed by metagenomics. Microorganisms 2020, 8, 1333. [Google Scholar] [CrossRef] [PubMed]
  26. Merkel, A.Y.; Chernyh, N.A.; Pimenov, N.V.; Bonch-Osmolovskaya, E.A.; Slobodkin, A.I. Diversity and metabolic potential of the terrestrial mud volcano microbial community with a high abundance of archaea mediating the anaerobic oxidation of methane. Life 2021, 11, 953. [Google Scholar] [CrossRef] [PubMed]
  27. Knittel, K.; Boetius, A. Anaerobic oxidation of methane: Progress with an unknown process. Annu. Rev. Microbiol. 2009, 63, 311–334. [Google Scholar] [CrossRef]
  28. Haroon, M.F.; Hu, S.; Shi, Y.; Imelfort, M.; Keller, J.; Hugenholtz, P.; Yuan, Z.; Tyson, G.W. Anaerobic oxidation of methane coupled to nitrate reduction in a novel archaeal lineage. Nature 2013, 500, 567–570. [Google Scholar] [CrossRef]
  29. Bhattarai, S.; Cassarini, C.; Lens, P.N.L. Physiology and distribution of archaeal methanotrophs that couple anaerobic oxidation of methane with sulfate reduction. Microbiol. Mol. Biol. Rev. 2019, 83, e00074-18. [Google Scholar] [CrossRef]
  30. Raghoebarsing, A.A.; Pol, A.; Van de Pas-Schoonen, K.T.; Smolders, A.J.; Ettwig, K.F.; Rijpstra, W.I.; Schouten, S.; Damsté, J.S.; den Camp, H.J.; Jetten, M.S.; et al. A microbial consortium couples anaerobic methane oxidation to denitrification. Nature 2006, 440, 918–921. [Google Scholar] [CrossRef]
  31. Beal, E.J.; House, C.H.; Orphan, V.J. Manganese-and iron-dependent marine methane oxidation. Science 2009, 325, 184–187. [Google Scholar] [CrossRef] [PubMed]
  32. Cai, C.; Leu, A.O.; Xie, G.J.; Guo, J.; Feng, Y.; Zhao, J.X.; Tyson, G.W.; Yuan, Z.; Hu, S. A methanotrophic archaeon couples anaerobic oxidation of methane to Fe (III) reduction. ISME J. 2018, 12, 1929–1939. [Google Scholar] [CrossRef]
  33. Leu, A.O.; Cai, C.; McIlroy, S.J.; Southam, G.; Orphan, V.J.; Yuan, Z.; Hu, S.; Tyson, G.W. Anaerobic methane oxidation coupled to manganese reduction by members of the Methanoperedenaceae. ISME J. 2020, 14, 1030–1041. [Google Scholar] [CrossRef]
  34. Zhang, X.; Joyce, G.H.; Leu, A.O.; Zhao, J.; Rabiee, H.; Virdis, B.; Tyson, G.W.; Yuan, Z.; McIlroy, S.J.; Hu, S. Multi-heme cytochrome-mediated extracellular electron transfer by the anaerobic methanotroph ‘Candidatus Methanoperedens nitroreducens’. Nat. Commun. 2023, 14, 6118. [Google Scholar] [CrossRef] [PubMed]
  35. Schreiber, L.; Holler, T.; Knittel, K.; Meyerdierks, A.; Amann, R. Identification of the dominant sulfate-reducing bacterial partner of anaerobic methanotrophs of the ANME-2 clade. Environ. Microbiol. 2010, 12, 2327–2340. [Google Scholar] [CrossRef]
  36. Murali, R.; Yu, H.; Speth, D.R.; Wu, F.; Metcalfe, K.S.; Crémière, A.; Laso-Pèrez, R.; Malmstrom, R.R.; Goudeau, D.; Woyke, T.; et al. Physiological potential and evolutionary trajectories of syntrophic sulfate-reducing bacterial partners of anaerobic methanotrophic archaea. PLoS Biol. 2023, 21, e3002292. [Google Scholar] [CrossRef]
  37. Wegener, G.; Krukenberg, V.; Ruff, S.E.; Kellermann, M.Y.; Knittel, K. Metabolic capabilities of microorganisms involved in and associated with the anaerobic oxidation of methane. Front. Microbiol. 2016, 7, 869. [Google Scholar] [CrossRef]
  38. Skennerton, C.T.; Chourey, K.; Iyer, R.; Hettich, R.L.; Tyson, G.W.; Orphan, V.J. Methane-fueled syntrophy through extracellular electron transfer: Uncovering the genomic traits conserved within diverse bacterial partners of anaerobic methanotrophic archaea. mBio 2017, 8, e00530-17. [Google Scholar] [CrossRef] [PubMed]
  39. Krukenberg, V.; Harding, K.; Richter, M.; Glöckner, F.O.; Gruber-Vodicka, H.R.; Adam, B.; Berg, J.S.; Knittel, K.; Tegetmeyer, H.E.; Boetius, A.; et al. Candidatus Desulfofervidus auxilii, a hydrogenotrophic sulfate-reducing bacterium involved in the thermophilic anaerobic oxidation of methane. Environ. Microbiol. 2016, 8, 3073–3091. [Google Scholar] [CrossRef]
  40. Gao, Y.; Wang, Y.; Lee, H.-S.; Jin, P. Significance of anaerobic oxidation of methane (AOM) in mitigating methane emission from major natural and anthropogenic sources: A review of AOM rates in recent publications. Environ. Sci. Adv. 2022, 1, 401–425. [Google Scholar] [CrossRef]
  41. Lavrushin, V.Y.; Kopf, A.; Dehle, A.; Stepanets, M.I. Formation of mudvolcanic fluids in Taman (Russia) and Kakhetia (Georgia): Evidence from boron isotopes. Lithol. Miner. Resour. 2003, 38, 120–153. [Google Scholar] [CrossRef]
  42. Khomyakova, M.A.; Merkel, A.Y.; Kopitsyn, D.S.; Slobodkin, A.I. Pelovirga terrestris gen. nov., sp. nov., anaerobic, alkaliphilic, fumarate-, arsenate-, Fe(III)- and sulfur-reducing bacterium isolated from a terrestrial mud volcano. System. Appl. Microbiol. 2022, 45, 126304. [Google Scholar] [CrossRef] [PubMed]
  43. Khomyakova, M.A.; Merkel, A.Y.; Segliuk, V.S.; Slobodkin, A.I. Desulfatitalea alkaliphila sp. nov., an alkalipilic sulfate- and arsenate- reducing bacterium isolated from a terrestrial mud volcano. Extremophiles 2023, 27, 12. [Google Scholar] [CrossRef] [PubMed]
  44. Slobodkina, G.B.; Merkel, A.Y.; Novikov, A.A.; Bonch-Osmolovskaya, E.A.; Slobodkin, A.I. Pelomicrobium methylotrophicum gen. nov., sp. nov. a moderately thermophilic, facultatively anaerobic, lithoautotrophic and methylotrophic bacterium isolated from a terrestrial mud volcano. Extremophiles 2020, 24, 177–185. [Google Scholar] [CrossRef] [PubMed]
  45. Slobodkina, G.; Ratnikova, N.; Merkel, A.; Kevbrin, V.; Kuchierskaya, A.; Slobodkin, A. Lithoautotrophic lifestyle of the widespread genus Roseovarius revealed by physiological and genomic characterization of Roseovarius autotrophicus sp. nov. FEMS Microbiol. Ecol. 2022, 98, fiac113. [Google Scholar] [CrossRef]
  46. Slobodkina, G.; Merkel, A.; Ratnikova, N.; Kuchierskaya, A.; Slobodkin, A. Sedimenticola hydrogenitrophicus sp. nov. a chemolithoautotrophic bacterium isolated from a terrestrial mud volcano, and proposal of Sedimenticolaceae fam. nov. in the order. Chromatiales. System. Appl. Microbiol. 2023, 46, 126451. [Google Scholar] [CrossRef]
  47. Rusanov, I.I.; Savvichev, A.S.; Yusupov, S.K.; Pimenov, N.V.; Ivanov, M.V. Production of exometabolites in the microbial oxidation of methane in ma-rine ecosystems. Microbiology 1998, 67, 590–596. [Google Scholar]
  48. Stroeva, A.R.; Klyukina, A.A.; Vidishcheva, O.N.; Poludetkina, E.N.; Solovyeva, M.A.; Pyrkin, V.O.; Gavirova, L.A.; Birkeland, N.-K.; Akhmanov, G.G.; Bonch-Osmolovskaya, E.A.; et al. Structure of benthic microbial communities in the Northeastern part of the Barents Sea. Microorganisms 2024, 12, 387. [Google Scholar] [CrossRef] [PubMed]
  49. Brettin, T.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Olsen, G.J.; Olson, R.; Overbeek, R.; Parrello, B.; Pusch, G.D.; et al. RASTtk: A modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci. Rep. 2015, 5, 8365. [Google Scholar] [CrossRef]
  50. Pimenov, N.V.; Savvichev, A.S.; Rusanov, I.I.; Lein, A.Y.; Ivanov, M.V. Microbiological processes of the carbon and sulfur cycle on cold methane seeps of the North Atlantic. Microbiology 2000, 69, 831–843. [Google Scholar]
  51. Ruff, S.E.; Biddle, J.F.; Teske, A.P.; Knittel, K.; Boetius, A.; Ramette, A. Global dispersion and local diversification of the methane seep microbiome. Proc. Natl. Acad. Sci. USA 2015, 112, 4015–4020. [Google Scholar] [CrossRef]
  52. Tavormina, P.L.; Hatzenpichler, R.; McGlynn, S.; Chadwick, G.; Dawson, K.S.; Connon, S.A.; Orphan, V.J. Methyloprofundus sedimenti gen. nov., sp. nov., an obligate methanotroph from ocean sediment belonging to the ‘deep sea-1’ clade of marine methanotrophs. Int. J. Syst. Evol. Microbiol. 2015, 65, 251–259. [Google Scholar] [CrossRef]
  53. Hirayama, H.; Takaki, Y.; Abe, M.; Imachi, H.; Ikuta, T.; Miyazaki, J.; Tasumi, E.; Uematsu, K.; Tame, A.; Tsuda, M.; et al. Multispecies populations of methanotrophic Methyloprofundus and cultivation of a likely dominant species from the Iheya North deep-sea hydrothermal field. Appl. Environ. Microbiol. 2022, 88, e0075821. [Google Scholar] [CrossRef] [PubMed]
  54. Chadwick, G.L.; Skennerton, C.T.; Laso-Perez, R.; Leu, A.O.; Speth, D.R.; Yu, H.; Morgan-Lang, C.; Hatzenpichler, R.; Goudeau, D.; Malmstrom, R.; et al. Comparative genomics reveals electron transfer and syntrophic mechanisms differentiating methanotrophic and methanogenic archaea. PLoS Biol. 2022, 20, e3001508. [Google Scholar] [CrossRef] [PubMed]
  55. Khomyakova, M.A.; Zavarzina, D.G.; Merkel, A.Y.; Klyukina, A.A.; Pikhtereva, V.A.; Gavrilov, S.N.; Slobodkin, A.I. The first cultivated representatives of the actinobacterial lineage OPB41 isolated from subsurface environments constitute a novel order Anaerosomatales. Front. Microbiol. 2022, 13, 1047580. [Google Scholar] [CrossRef]
  56. Wunder, L.C.; Aromokeye, D.A.; Yin, X.; Richter-Heitmann, T.; Willis-Poratti, G.; Schnakenberg, A.; Otersen, C.; Dohrmann, I.; Römer, M.; Bohrmann, G.; et al. Iron and sulfate reduction structure microbial communities in (sub-)Antarctic sediments. ISME J. 2021, 12, 3587–3604. [Google Scholar] [CrossRef]
  57. Wang, L.; Shao, Z. Aerobic denitrification and heterotrophic sulfur oxidation in the genus Halomonas revealed by six novel species characterizations and genome-based analysis. Front. Microbiol. 2021, 12, 652766. [Google Scholar] [CrossRef] [PubMed]
  58. Emerson, D.; Rentz, J.A.; Lilburn, T.G.; Davis, R.E.; Aldrich, H.; Chan, C.; Moyer, C.L. A novel lineage of Proteobacteria involved in formation of marine Fe-oxidizing microbial mat communities. PLoS ONE 2007, 2, e667. [Google Scholar] [CrossRef] [PubMed]
  59. Parks, D.H.; Chuvochina, M.; Waite, D.W.; Rinke, C.; Skarshewski, A.; Chaumeil, P.A.; Hugenholtz, P. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat. Biotechnol. 2018, 36, 996–1004. [Google Scholar] [CrossRef]
  60. Timmers, P.H.A.; Welte, C.U.; Koehorst, J.J.; Plugge, C.M.; Jetten, M.S.M.; Stams, A.J.M. Reverse methanogenesis and respiration in methanotrophic archaea. Archaea 2017, 2017, 1654237. [Google Scholar] [CrossRef]
  61. Mehta, P.; Mishra, S.; Ghose, T. Methanol accumulation by resting cells of Methylosinus trichosporium (I). J. Gen. Appl. Microbiol. 1987, 33, 221–229. [Google Scholar] [CrossRef]
  62. Farr, O.; Gaudu, N.; Danger, G.; Russell, M.J.; Ferry, D.; Nitschke, W.; Duval, S. Methanol on the rocks: Green rust transformation promotes the oxidation of methane. J. R. Soc. Interface 2023, 206, 20230386. [Google Scholar] [CrossRef] [PubMed]
  63. Berg, I.A. Ecological aspects of the distribution of different autotrophic CO2 fixation pathways. Appl. Environ. Microbiol. 2011, 77, 1925–1936. [Google Scholar] [CrossRef] [PubMed]
  64. Nunoura, T.; Chikaraishi, Y.; Izaki, R.; Suwa, T.; Sato, T.; Harada, T.; Mori, K.; Kato, Y.; Miyazaki, M.; Shimamura, S.; et al. A primordial and reversible TCA cycle in a facultatively chemolithoautotrophic thermophile. Science 2018, 359, 559–563. [Google Scholar] [CrossRef] [PubMed]
  65. Mall, A.; Sobotta, J.; Huber, C.; Tschirner, C.; Kowarschik, S.; Bačnik, K.; Mergelsberg, M.; Boll, M.; Hügler, M.; Eisenreich, W.; et al. Reversibility of citrate synthase allows autotrophic growth of a thermophilic bacterium. Science 2018, 359, 563–567. [Google Scholar] [CrossRef] [PubMed]
  66. Sánchez-Andrea, I.; Guedes, I.A.; Hornung, B.; Boeren, S.; Lawson, C.E.; Sousa, D.Z.; Bar-Even, A.; Claassens, N.J.; Stams, A.J.M. The reductive glycine pathway allows autotrophic growth of Desulfovibrio desulfuricans. Nat. Commun. 2020, 11, 5090. [Google Scholar] [CrossRef] [PubMed]
  67. Pitcher, R.S.; Watmough, N.J. The bacterial cytochrome cbb3 oxidases. Biochim. Biophys. Acta 2004, 1655, 388–399. [Google Scholar] [CrossRef] [PubMed]
  68. Sorokin, D.Y.; Tourova, T.P.; Mussmann, M.; Muyzer, G. Dethiobacter alkaliphilus gen. nov. sp. nov., and Desulfurivibrio alkaliphilus gen. nov. sp. nov.: Two novel representatives of reductive sulfur cycle from soda lakes. Extremophiles 2008, 12, 431–439. [Google Scholar]
  69. Flood, B.E.; Jones, D.S.; Bailey, J.V. Sedimenticola thiotaurini sp. nov., a sulfur-oxidizing bacterium isolated from salt marsh sediments, and emended descriptions of the genus Sedimenticola and Sedimenticola selenatireducens. Int. J. Syst. Evol. Microbiol. 2015, 65, 2522–2530. [Google Scholar] [CrossRef]
  70. Frigaard, N.U.; Dahl, C. Sulfur metabolism in phototrophic sulfur bacteria. Adv. Microb. Physiol. 2009, 54, 103–200. [Google Scholar] [CrossRef] [PubMed]
  71. de Zwart, J.M.M.; Nelisse, P.N.; Kuenen, J.G. Isolation and characterization of Methylophaga sulfidovorans sp. nov.: An obligately methylotrophic, aerobic, dimethylsulfide oxidizing bacterium from microbial mat. FEMS Microbiol. Ecol. 1996, 20, 261–270. [Google Scholar] [CrossRef]
  72. Boden, R.; Kelly, D.P.; Murrell, J.C.; Schafer, H. Oxidation of dimethylsulfide to tetrathionate by Methylophaga thiooxidans sp. nov.: A new link in the sulfur cycle. Environ. Microbiol. 2010, 12, 2688–2699. [Google Scholar] [CrossRef] [PubMed]
  73. Shi, L.; Dong, H.; Reguera, G.; Beyenal, H.; Lu, A.; Liu, J.; Yu, H.Q.; Fredrickson, J.K. Extracellular electron transfer mechanisms between microorganisms and minerals. Nat. Rev. Microbiol. 2016, 14, 651–662. [Google Scholar] [CrossRef] [PubMed]
  74. Schuetz, B.; Schicklberger, M.; Kuermann, J.; Spormann, A.M.; Gescher, J. Periplasmic electron transfer via the c-type cytochromes MtrA and FccA of Shewanella oneidensis MR-1. Appl. Environ. Microbiol. 2009, 75, 7789–7796. [Google Scholar] [CrossRef] [PubMed]
  75. Frolova, A.A.; Merkel, A.Y.; Kevbrin, V.V.; Kopitsyn, D.S.; Slobodkin, A.I. Sulfurospirillum tamanensis sp. nov., a facultatively anaerobic alkaliphilic bacterium from a terrestrial mud volcano. Microbiology 2023, 92, 21–29. [Google Scholar] [CrossRef]
  76. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; von Haeseler, A.; Lanfear, R. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef] [PubMed]
  77. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.; von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef]
  78. Minh, B.Q.; Nguyen, M.A.; von Haeseler, A. Ultrafast approximation for phylogenetic bootstrap. Mol. Biol. Evol. 2013, 30, 1188–1195. [Google Scholar] [CrossRef]
  79. Anisimova, M.; Gascuel, O. Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst. Biol. 2006, 55, 539–552. [Google Scholar] [CrossRef]
Figure 1. Rates of microbial methane oxidation (nmol CH4 cm−3 day−1) and distribution of carbon from methane into fractions.
Figure 1. Rates of microbial methane oxidation (nmol CH4 cm−3 day−1) and distribution of carbon from methane into fractions.
Microorganisms 12 01349 g001
Figure 2. The relative abundance of prokaryotic phyla (Silva and GTDB taxonomy), which is more than 1% in at least one of the samples.
Figure 2. The relative abundance of prokaryotic phyla (Silva and GTDB taxonomy), which is more than 1% in at least one of the samples.
Microorganisms 12 01349 g002
Figure 3. Relative abundance of prokaryotic genera (Silva and GTDB taxonomy), which is more than 5% in at least one of the samples. Taxa are grouped based on the metabolic properties of phylogenetic relatives.
Figure 3. Relative abundance of prokaryotic genera (Silva and GTDB taxonomy), which is more than 5% in at least one of the samples. Taxa are grouped based on the metabolic properties of phylogenetic relatives.
Microorganisms 12 01349 g003
Figure 4. The relative abundance of different ANME groups in the studied samples.
Figure 4. The relative abundance of different ANME groups in the studied samples.
Microorganisms 12 01349 g004
Figure 5. Metabolic capabilities encoded in MAGs of the sample C14 from the TMV Gladkovsky. Genes present in MAGs are shown in the bottom row. Gene abbreviations are explained in the text and listed in Supplementary Table S1.
Figure 5. Metabolic capabilities encoded in MAGs of the sample C14 from the TMV Gladkovsky. Genes present in MAGs are shown in the bottom row. Gene abbreviations are explained in the text and listed in Supplementary Table S1.
Microorganisms 12 01349 g005
Table 1. Characteristics of the studied sites.
Table 1. Characteristics of the studied sites.
Sample #VolcanoSampling SitepHSalinity, g L−1Eh,
mV
CH4, mMCl, mMPO4, mMSO42−, mMNO3, mM
C01Gnilaya GoraGryphon8.4720−295101.1241.10.005.790.00
C02Gnilaya GoraGryphon8.6614−175556.7158.80.071.070.00
C03Gnilaya GoraSalsa9.1020−320201.1412.00.000.000.00
C04Gnilaya GoraMud pool8.6412−348536.3162.80.020.090.92
C05Gnilaya GoraMud pool8.9420−30247.7216.10.001.950.00
C10GladkovskyMud pool6.7520−124668.0307.50.000.060.00
C11GladkovskyMud pool6.7220−110661.5300.80.000.080.00
C12GladkovskyMud pool7.0730153.3402.20.000.000.00
C13GladkovskyMud pool7.0830−255105.9393.40.000.000.00
C14GladkovskyMud pool7.0220−217439.1291.20.000.000.00
C16KuchugurskyMain crater7.4315−212640.5153.90.008.540.00
C17KuchugurskyMud pool6.9220ndnd147.40.0063.520.00
C19SemigorskyMud pool8.2010ndnd5.70.000.080.00
C20SemigorskyMud pool8.359ndnd37.20.000.020.00
C21ShugoGryphon7.2320−196456.8199.20.000.200.00
C23ShugoGryphon6.9412ndnd95.40.000.210.00
C24ShugoGryphon7.4120−157917.3193.70.000.210.00
nd—not determined.
Table 2. Overview of all genome bins >50% complete and with <5% contamination.
Table 2. Overview of all genome bins >50% complete and with <5% contamination.
Bin IDDomainTaxon *Abundance %Completeness %Contamination %# ContigsGenome Size, Mbp16S rRNA Gene, bp
C14-01Bc_Acidimicrobiia19.04100.002.141813.911171
C14-19Bd_Bacteria18.5794.780.841573.191558
C14-17Af_Ca. Methanocomedenaceae14.2395.420.001312.06-
C14-08Af_Ca. Methanocomedenaceae11.5198.690.00741.82-
C14-05Bg_Methyloprofundus9.4999.311.551893.28-
C14-10Bo_Bacteroidales3.1397.851.90503.481533
C14-03Bf_Methylophagaceae2.5699.820.65962.44-
C14-15Bd_Bacteria2.3896.400.002552.72-
C14-28Bp_Paceibacteria2.2876.851.12681.171474
C14-27Bo_Anaerolineales2.2180.182.73992.691517
C14-06Bo_Bacteroidales1.8199.281.90743.81-
C14-21Bg_Brevefilum1.2790.002.731042.661522
C14-02Bo_Desulfuromonadales1.10100.000.221182.72-
C14-14Bc_Myxococcia1.0496.552.60756.11768
C14-09Bo_Bacteroidales1.0098.571.191194.61-
C14-29Bp_Patescibacteria0.9073.760.99491.461512
C14-13Bo_Anaerosomatales0.7996.663.61521.98432
C14-04Bo_Desulfobacterales0.7899.352.262694.94556
C14-16Bo_Deferrisomatales0.7095.702.694705.78-
C14-07Bf_Desulfurivibrionaceae0.6398.802.352872.52-
C14-12Bp_Zixibacteria0.4696.701.202423.06-
C14-18Ag_Ca. Methanovorans0.4495.010.982631.88-
C14-11Bf_Spirochaetaceae0.4397.202.003443.68-
C14-24Bf_Trueperaceae0.4084.952.976663.26696
C14-20Af_Thorarchaeaceae0.3990.183.277034.031556
C14-25Bo_Bacteroidales0.3683.633.8312356.00-
C14-23Bg_Lutibacter0.3485.013.246742.49-
C14-22Bg_Sedimenticola 0.3485.631.396633.03-
C14-26Bc_Anaerolineae0.3380.490.4810214.22-
* The lowest taxon with a published name according the GTDB classification.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Slobodkin, A.I.; Rusanov, I.I.; Slobodkina, G.B.; Stroeva, A.R.; Chernyh, N.A.; Pimenov, N.V.; Merkel, A.Y. Diversity, Methane Oxidation Activity, and Metabolic Potential of Microbial Communities in Terrestrial Mud Volcanos of the Taman Peninsula. Microorganisms 2024, 12, 1349. https://doi.org/10.3390/microorganisms12071349

AMA Style

Slobodkin AI, Rusanov II, Slobodkina GB, Stroeva AR, Chernyh NA, Pimenov NV, Merkel AY. Diversity, Methane Oxidation Activity, and Metabolic Potential of Microbial Communities in Terrestrial Mud Volcanos of the Taman Peninsula. Microorganisms. 2024; 12(7):1349. https://doi.org/10.3390/microorganisms12071349

Chicago/Turabian Style

Slobodkin, Alexander I., Igor I. Rusanov, Galina B. Slobodkina, Aleksandra R. Stroeva, Nikolay A. Chernyh, Nikolai V. Pimenov, and Alexander Y. Merkel. 2024. "Diversity, Methane Oxidation Activity, and Metabolic Potential of Microbial Communities in Terrestrial Mud Volcanos of the Taman Peninsula" Microorganisms 12, no. 7: 1349. https://doi.org/10.3390/microorganisms12071349

APA Style

Slobodkin, A. I., Rusanov, I. I., Slobodkina, G. B., Stroeva, A. R., Chernyh, N. A., Pimenov, N. V., & Merkel, A. Y. (2024). Diversity, Methane Oxidation Activity, and Metabolic Potential of Microbial Communities in Terrestrial Mud Volcanos of the Taman Peninsula. Microorganisms, 12(7), 1349. https://doi.org/10.3390/microorganisms12071349

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop