Karen M. Houghton
Carlo R. Carere
Matthew B. Stott
Ian R. McDonald- 1Te Pū Ao | GNS Science, Wairakei Research Centre, Taupō, New Zealand
- 2Te Aka Mātuatua | School of Science, Te Whare Wānanga o Waikato | University of Waikato, Hamilton, New Zealand
- 3Te Tari Pūhanga Tukanga Matū | Department of Chemical and Process Engineering, Te Whare Wānanga o Waitaha | University of Canterbury, Christchurch, New Zealand
- 4Te Kura Pūtaiao Koiora | School of Biological Sciences, Te Whare Wānanga o Waitaha | University of Canterbury, Christchurch, New Zealand
Geothermal areas represent substantial point sources for greenhouse gas emissions such as methane. While it is known that methanotrophic microorganisms act as a biofilter, decreasing the efflux of methane in most soils to the atmosphere, the diversity and the extent to which methane is consumed by thermophilic microorganisms in geothermal ecosystems has not been widely explored. To determine the extent of biologically mediated methane oxidation at elevated temperatures, we set up 57 microcosms using soils from 14 Aotearoa-New Zealand geothermal fields and show that moderately thermophilic (>40°C) and thermophilic (>60°C) methane oxidation is common across the region. Methane oxidation was detected in 54% (n = 31) of the geothermal soil microcosms tested at temperatures up to 75°C (pH 1.5–8.1), with oxidation rates ranging from 0.5 to 17.4 μmol g−1 d−1 wet weight. The abundance of known aerobic methanotrophs (up to 60.7% Methylacidiphilum and 11.2% Methylothermus) and putative anaerobic methanotrophs (up to 76.7% Bathyarchaeota) provides some explanation for the rapid rates of methane oxidation observed in microcosms. However, not all methane oxidation was attributable to known taxa; in some methane-consuming microcosms we detected methanotroph taxa in conditions outside of their known temperature range for growth, and in other examples, we observed methane oxidation in the absence of known methanotrophs through 16S rRNA gene sequencing. Both of these observations suggest unidentified methane oxidizing microorganisms or undescribed methanotrophic syntrophic associations may also be present. Subsequent enrichment cultures from microcosms yielded communities not predicted by the original diversity studies and showed rates inconsistent with microcosms (≤24.5 μmol d−1), highlighting difficulties in culturing representative thermophilic methanotrophs. Finally, to determine the active methane oxidation processes, we attempted to elucidate metabolic pathways from two enrichment cultures actively oxidizing methane using metatranscriptomics. The most highly expressed genes in both enrichments (methane monooxygenases, methanol dehydrogenases and PqqA precursor peptides) were related to methanotrophs from Methylococcaceae, Methylocystaceae and Methylothermaceae. This is the first example of using metatranscriptomics to investigate methanotrophs from geothermal environments and gives insight into the metabolic pathways involved in thermophilic methanotrophy.
1. Introduction
The greenhouse gas methane (CH4) is responsible for a large proportion of global climate change, being 28 times more effective than carbon dioxide (CO2) at absorbing infra-red radiation (Forster et al., 2021). The global methane budget is ~730 Tg yr−1 from both natural and anthropogenic sources; with ~5% of this methane consumed by methanotrophic bacteria before it reaches the atmosphere (Canadell et al., 2021). It is therefore important to understand the distribution and diversity of methanotrophic bacteria to gain insight into the ecological processes that modulate methane emissions and its oxidation.
Aerobic methanotrophic bacteria have been detected and/or isolated from nearly all known habitats associated with methane emissions; including wetlands, rice paddy fields, landfills, and geothermal areas (Tsubota et al., 2005; Geymonat et al., 2011; Danilova et al., 2013; Henneberger et al., 2014). While knowledge concerning methanotroph ecology within many of these ecosystems is extensive (Hanson and Hanson, 1996; Jiang et al., 2010), the understanding of the distribution and diversity of methane oxidizing bacteria in geothermal fields remains limited (Houghton et al., 2019). Geothermal methane emissions from magma degassing amounts to ≤71 Tg yr−1 (Canadell et al., 2021), but geothermal and volcanic areas are inhospitable for most known species of methanotrophs due to acidic pH, hot temperatures and low oxygen concentrations (Op den Camp et al., 2009). Evidence of methanotrophic activity in geothermal areas was first reported in 2005 (Castaldi and Tedesco, 2005), and a number of thermophilic aerobic methanotrophs have been isolated from geothermal soils, including the description of three genera (Methylacidiphilum, Methylacidimicrobium, Candidatus Methylacidithermus) within the phylum Verrucomicrobiota (van Teeseling et al., 2014; Picone et al., 2021a). Thermophilic methanotrophs from these areas demonstrate metabolic flexibility, using or consuming a wide range of energy sources in addition to methane, including hydrogen (Carere et al., 2017), hydrogen sulfide (Schmitz et al., 2023), methanethiol (Schmitz et al., 2022) and short-chain alkanes (Awala et al., 2021). However, little is known about the phylogenetic diversity of thermophilic methanotrophs from geothermal areas, perhaps due to the difficulty of culturing methane-oxidisers, from a failure to replicate environmental conditions in vitro, or because there are biological restrictions on methane oxidation at hotter temperatures, such as methane solubility or enzyme instability (Houghton et al., 2019).
Molecular techniques have frequently been used to detect methanotrophs within environmental samples. These methods commonly include DNA amplicon sequencing of the 16S rRNA gene (Bodrossy et al., 1997), or functional genes associated with methanotrophy such as methane monooxygenase (McDonald et al., 1995; Sharp et al., 2012; Gagliano et al., 2014), and metagenomics (Håvelsrud et al., 2011; Nguyen et al., 2018; Taubert et al., 2019). However, DNA sequencing techniques provide little information about the activity of microbial community members, gene expression or responses to environmental conditions (Wang et al., 2012). Considering up to 80% of soil microorganisms are believed to exist in a state of dormancy (Lennon and Jones, 2011), DNA sequencing of these microorganisms can obscure assessments of metabolic activity within an environment. Previous studies targeting methanotrophs have used transcriptomics within hydrothermal vents (Lesniewski et al., 2012; Olins et al., 2017), peatlands (Liebner and Svenning, 2013; Esson et al., 2016) and landfill soil (Henneberger et al., 2014), but not geothermal soil or hot spring ecosystems.
In this study, we combined the use of 16S rRNA gene amplicon sequencing on geothermal samples from throughout the Taupō Volcanic Zone (TVZ; Aotearoa-New Zealand) with the quantification of methane oxidation rates to detect the presence and activity of thermophilic and thermotolerant methanotrophs. The TVZ spans 6,000 km2 of the central North Island of Aotearoa-New Zealand and extends from the active volcanoes of Whakaari/White Island to Mount Ruapehu (Figure 1). The area contains more than 20 geothermal fields that are characterized by an abundance of hot springs, geysers, and mud pools (Giggenbach, 1995). Prior analysis of gases emitted from geothermal fields within the TVZ has found that most systems release 0.1 to 5.0% CH4 (v/v), with some sites, such as Golden Springs, emitting up to 19.3% CH4 (v/v) (Giggenbach, 1995).
Figure 1. Map of the Taupō Volcanic Zone (TVZ), Aotearoa-New Zealand, and a schematic of the experimental study for this research. A total of 57 samples were collected (red circles) from 14 geothermal fields. (1) Microcosms were prepared using a single soil (n = 43), sediment (n = 11) or water/biofilm (n = 3) sample with a headspace of (v/v) 10% CH4, 1.3% CO2 in air. (2, 3) The microbial community composition of the 31 microcosms that showed positive CH4 oxidation were determined by 16S rRNA gene amplicon sequencing. (4) Enrichment cultures using a subsample of the microcosm substrate were inoculated into NMS media and other media informed by methanotrophic taxa detected in the microcosm microbial communities. (5) Two enrichments, TOK1 (Tokaanu) and GDS1 (Golden Springs), showing strong CH4 consumption were further subcultured into ANMS (46°C) and NMS (37°C) media, respectively, for metatranscriptomic analyses. Individual microcosm and enrichment cultivations were undertaken as close as possible to their measured in situ temperatures (see Table 1). Map adapted from Power et al. (2018), and schematic created in www.biorender.com/.
To determine thermophilic methanotrophic activity in these geothermal areas, we emended soil and sediment microcosms with various media designed to promote methanotrophic growth. Methane oxidation in each microcosm was quantified along with microbial community composition. Two enrichments were then selected which showed both elevated rates of methane oxidation and abundant methanotroph communities. Metatranscriptomic analyses of the two enrichments were used to assess potential methanotrophic metabolic activity.
2. Materials and methods
An overview of the experimental workflow is shown in Figure 1.
2.1. Microcosms
2.1.1. Sampling location overview and characterization
Microcosm samples were collected from geothermally heated soils (with surface temperatures >55°C), hot spring waters, sediments or associated biofilms from 14 geothermal fields within the TVZ (Figure 1; Table 1). Site selection was partially informed from the 1,000 Springs Project dataset (Power et al., 2018) by screening for sites with either high in situ concentrations of methane, or with a high proportion of 16S rRNA gene sequences assigned as putative methanotrophs. Individual physicochemical observations for these sites and associated microcosms are listed in Supplementary Table 1.
Table 1. Physicochemical characteristics, number of microcosms, and location of the 14 geothermal systems from the Taupō Volcanic Zone (TVZ), Aotearoa-New Zealand sampled for this study.
Ground temperatures were measured using a 51 II thermal probe (Fluke) inserted into the soil or sediment at 20 cm depth. Soil gas samples were collected, where possible, using a custom-built stainless steel sampling tube inserted into the soil to the same depth. Before collecting the gas sample, and to avoid atmospheric contamination, 700 ml of gas (>3× tube volume) was drawn through the sampling tube and discarded. 25 ml of gas was then slowly extracted from the soil and injected through a gas-tight three-way valve into a pre-evacuated ‘Air and Gas Sampling Bag’ (Calibrated Instruments Inc). Gas samples were processed within 24 h on a Peak Performer 1 Gas Analyzer (Peak Laboratories LIC) equipped with a flame ionization detector (FID) and a Unibeads 60/80 column. Soil/sediment samples (~2 g) were collected using a trowel and spatula previously sterilized with ethanol and were immediately transferred to pre-autoclaved serum vials (Wheaton Industries) which were then capped with sterile butyl rubber stoppers and sealed with aluminum crimps. A separate duplicate soil/sediment sample was taken for pH measurement by suspending ~1 g of soil in 10 mL RO H2O at room temperature (Table 1; Supplementary Table 1).
2.1.2. Methane oxidation by microcosms
A gas mixture (v/v), consisting of 80% CH4, 15% CO2 and 5% O2 (14.3 mL), was injected into the serum vial headspace through a sterile filter resulting in final gas headspace concentrations of 10% CH4, 1.3% CO2 and 21.6% O2 of each microcosm (Supplementary Table 1). Control vials with no soil were set up with the same gas headspace. Vials were incubated at temperatures (46, 50, 60, 68, 70, 75 or 80°C) approximating the closest in situ sample temperature (controls incubated at 60°C). Microcosms from Golden Springs (GDS) were incubated at 50°C to target thermophilic methanotrophs, although the in situ temperature was 38°C. At 3–4 day intervals, vials were moved to a 20°C water bath and equilibrated to room temperature (20°C) for 5 h, to ensure accurate comparisons of pressure and moles of methane consumed Headspace gas compositions were then quantified by GC-FID by removing 500 μL with a gastight syringe (SGE Analytical Science) and diluting to 5 mL with air.
Methane oxidation rates (μmol g−1 d−1 wet weight) were calculated using linear regression. If microcosms showed a greater loss of methane than controls (following sampling without replacement) for five consecutive readings, and if the coefficient of determination of linear regression (R2) was >0.5, microcosms were designated as positive for methane oxidation. Supplementary Figure 1 shows methane oxidation rate graphs for selected soil microcosms. After positive oxidation was determined, or following at least 4 weeks of headspace measurements for microcosms that did not oxidize methane, DNA was extracted from all microcosms.
2.1.3. Microcosm DNA extraction and sequencing
DNA was extracted from soil microcosms using a modified protocol for the NucleoSpin Soil kit (Macherey-Nagel) (Houghton and Stewart, 2019). DNA was amplified using universal primers for the V4 region (515F, 806R) of the 16S rRNA gene (Caporaso et al., 2011). PCR was carried out in 50 μL reaction volumes containing 100 μM dNTPs, 0.5 μM primers, 1 U i-Taq (iNtRON Biotechnology) and 7 μL of an enhancer solution (2.7 M betaine, 0.2 M trehalose, 6.7 mM DTT, 0.06 mg ml−1 BSA and 0.07% DMSO). The final concentration of MgCl2 was 1.5 mM. DNA templates from microcosms were used at final concentrations of 10–50 ng reaction−1. Three PCR amplicons (~300 bp) for each sample were pooled. The amplicons were then purified using the NucleoSpin Gel and PCR Clean-up kit (Macherey-Nagel) and Agencourt AmPure XP (Beckman Coulter). Amplicon libraries using the PCR products were prepared and sequenced by Macrogen Inc. Sequencing data was deposited in the NCBI BioProject database (PRJNA766707 and PRJNA546003).
2.1.4. Community 16S rRNA gene sequence processing and diversity metric assessment
The quality of raw read data was assessed using FastQC (RRID:SCR_014583) (Andrews, 2010). Paired-end sequence reads were merged and filtered using USEARCH v7.0, with a maximum expected error of 1 (Edgar, 2010). Remaining sequences either >500 bp (to remove poor quality sequences at the end of long reads) or <200 bp (the minimum required for taxonomic classification) were removed using mothur (RRID:SCR_011947) v1.35.1 (Schloss et al., 2009). A de novo database of ≥97% similar sequence Operational Taxonomic Units (OTUs) was created in USEARCH (Edgar, 2010). Raw sequences were mapped against this de novo database to generate counts of sequences matching OTUs (i.e., taxa) for each sample. Using QIIME (RRID:SCR_008249) v1.9.1 (Caporaso et al., 2010), taxonomy was assigned to each OTU by using the RDP classifier v2.2 (Wang et al., 2007) with a confidence threshold of 0.5 and trained on the SILVA (RRID:SCR_006423) 16S rRNA gene database (version 123) (Quast et al., 2013). Chloroplast and mitochondrial OTUs were removed and all samples were rarefied to the lowest sample read count (n = 149,400).
All OTUs identified as belonging to the methanotrophic families Beijerinckiaceae or Methylocystaceae (Alphaproteobacteria); the order Methylococcales (Gammaproteobacteria); or as Methylacidiphilum (Verrucomicrobiota) were manually checked against the NCBI Nucleotide (RRID:SCR_004860) database using a discontiguous megablast (Altschul et al., 1990). OTUs that were > 90% related to a described methanotrophic species were identified as putative methanotrophs. No OTUs were identified from the NC10/Methylomirabilis phylum were identified.
The core diversity workflow within QIIME (Caporaso et al., 2010) was used to analyze sequencing data, by creating multiple rarefactions of the data, calculating and comparing alpha and beta diversity of the samples, and summarizing taxa across all samples and as a function of methane oxidation status. Chao1 [a non-parametric estimate of species richness (Chao, 1984)] and the Shannon (Shannon and Weaver, 1963) and Simpson (Simpson, 1949) indices (quantitative measures of diversity and abundance) were used as alpha diversity metrics. A non-parametric two-sample t-test on the Shannon and Simson indices was performed using QIIME. To assess significant differences between genera, OTUs from all microcosms were filtered to remove unclassified sequences, and Kruskal-Wallis tests (McDonald, 2014) performed on genus clusters.
2.2. Methanotroph enrichments
All microcosms displaying methane oxidation were selected for enrichment of methanotrophs using a variety of media (Supplementary Table 2; Supplementary Material 1). Each microcosm was inoculated into a modified Nitrate Mineral Salts medium (mNMS), which was designed to be a non-specific methanotroph medium. In addition to mNMS, soil from each microcosm with positive methane consumption was also inoculated into at least one other growth medium. The additional medium was selected after consideration of the soil pH and putative methanotrophs detected in the microcosm communities (see Supplementary Table 2; Supplementary Material 1 for further details). All media were adjusted to reflect the pH of the starting microcosm source material.
Between 0.3–0.5 g of the soil or sediment sample from selected microcosms (Supplementary Table 2) were aseptically transferred to between two to four serum bottles each containing 40 mL of a different sterile medium for enrichment. The final gas headspace composition (v/v) was ~10% CH4, 1% CO2 and 22% O2 (balance N2). Enrichment serum vials were incubated at temperatures reflecting in situ sample temperature (37, 46, 50, 60, 70, or 75°C), with shaking at 150 r.p.m. Every 3–4 days, gas compositions were measured as previously above. Rates of methane oxidation (μmol d−1) were calculated using linear regression. Supplementary Figure 2 shows methane oxidation rate graphs for selected enrichments. For those enrichments that did not oxidize all methane within the headspace, a coefficient of determination (R2) greater than 0.5 was designated as positive for methane oxidation.
2.3. Metatranscriptome enrichments
2.3.1. Culture conditions for RNA extraction
Two methane-oxidizing geothermal enrichment cultures (GDS1 and TOK7) were selected for metatranscriptome sequencing Supplementary Material 1 Supplementary methods, Supplementary Table 3) based upon the relative abundances of methanotroph-associated 16S rRNA gene sequences and measured methane consumption rates. New enrichment cultures were prepared by inoculating the initial GDS1 and TOK7 enrichment cultures into fresh media in triplicate. Cultures were prepared in 114 mL serum vials with 30 mL of the relevant medium. Each sealed vial contained an air headspace supplemented with (v/v) ~ 10% CH4 and ~1% CO2. GDS1 enrichments were incubated at 37°C, the optimum growth temperature for Methylococcus spp. (Bowman et al., 1993), while TOK7 enrichments were incubated at 46°C to promote growth of Methylothermaceae (Houghton et al., 2019).
The extent of methane oxidation within each metatranscriptome enrichment was assessed by measuring the headspace gas composition over a 24-h period using GC-FID as described above. RNA was stabilized by the addition of 20 mL of RNAlater (Invitrogen) following 24 h of incubation and verification of methane consumption (> 5,000 μmoles CH4 d−1). Cells were collected via centrifugation at 16,000× g for 15 min at 10°C, with the cell pellet resuspended in 10 mL of RNAlater.
2.3.2. RNA extraction and sequencing
RNA extraction and sequencing were performed by Novogene as per (Liu et al., 2015). The cDNA library was sequenced using the NovaSeq SP platform (see Supplementary Material 1 Supplementary methods).
2.3.3. Bioinformatics
cDNA sequence analyses were performed within the Galaxy (RRID:SCR_006281) web platform (Afgan et al., 2016) with all settings for the individual pipeline wrappers using default settings unless stated otherwise. Initial raw read data was assessed using FastQC (v0.69) (Andrews, 2010), with reads trimmed using Trimmomatic (RRID:SCR_011848) (v0.36.2) (Bolger et al., 2014). The average quality required was set to a Phred score of 25 for GDS samples, and 20 for TOK samples, based on their FastQC reports. Reads with a length < 120 bp after trimming were removed. Forward and reverse reads for each sample were concatenated into two files using cat (v.0.1.0) (Gruening, 2014) and then assembled into contigs using Trinity (RRID:SCR_013048) (v2.2.0) (Grabherr et al., 2011).
Predicted genes within the contigs were identified using MetaGeneMark (v3.25) (Besemer and Borodovsky, 1999; Zhu et al., 2010). Identical sequences were removed using the standalone Java applet, DuplicatesFinder (Khan, 2009). The remaining sequences were clustered using the program CD-HIT (RRID:SCR_007105) (Huang et al., 2010) with sequences clustered at 90% nucleotide sequence identity for the GDS1 sample and at 85% nucleotide sequence identity for the TOK7 sample. One representative sequence from each cluster identified by cd-hit-est was retained for further analysis. The predicted genes were analyzed using BLASTX (RRID:SCR_001653) (Altschul et al., 1997) on a downloaded BLAST+ executable (v2.7.1) (Camacho et al., 2009) and a local copy of the non-redundant (nr) protein sequences (NCBI - 19/02/2018). Predicted proteins from the GDS1 and TOK7 samples were manually searched for modules within the “Methane Metabolism,” and “Nitrogen Metabolism” pathways in the KEGG (RRID:SCR_012773) database (Kanehisa et al., 2016). Finally, transcripts from the trimmed, non-concatenated datasets were quantified against the predicted genes from the Trinity assemblies using Salmon (RRID:SCR_017036) (0.8.2) (Patro et al., 2017), and the most highly expressed transcripts, based on Transcripts Per Million (TPM) across all three biological replicates, were identified for each sample.
3. Results and discussion
3.1. Rapid methane oxidation observed within thermophilic microcosms
Microcosm experiments illustrate that methanotrophic activity in high-temperature geothermal ecosystems correlates poorly with community abundance of known methanotrophs. Following incubation within a methane headspace, 31 of the 57 geothermal microcosms (54.3%) were identified as positive for methane oxidation (Figure 2A). Methane oxidation rates ranged from 0.5 μmol g−1 d−1 wet weight (LPR16, 70°C) to a maximum of 17.4 μmol g−1 d−1 wet weight (OKO2, 70°C; Supplementary Figure 1; Supplementary Table 1). These rates are comparable with previous studies of New Zealand geothermal areas, with maximum recorded oxidation rates of 7.0 μmol g−1 d−1 wet weight (55°C) (Sharp et al., 2012), and 20.4 μmol g−1 d−1 wet weight (37°C) or 12.7 μmol g−1 d−1 wet weight (65°C) (Sharp et al., 2014b). Other surveys of geothermal soils have reported maximum methane oxidation rates between 0.8 μmol g−1 d−1 wet weight (Italy, 50°C) (Gagliano et al., 2014) and 99 μmol g−1 d−1 wet weight (Canada, 45°C) (Sharp et al., 2014a). These rates are an order of magnitude greater than those recorded for forest soils (~0.3 nmol g−1 d−1 wet weight) (Zeng et al., 2019) and lake sediments (~300 nmol g−1 d−1 wet weight) (Bornemann et al., 2016), and comparable with methane oxidation in wetlands (17.3 μmol g−1 d−1 wet weight) (Esson et al., 2016).
Figure 2. Microcosm methane oxidation rates and microbial diversity. (A) Bubble charts illustrating observed methane oxidation rates of soil microcosms and as a function of temperature and pH. Green circles represent samples that oxidized methane and are labeled with an identifier. Bubble size is proportional to the rate of observed methane oxidation. Black dots denote samples that did not oxidize methane. (B) Bar chart of phyla identified in 16S rRNA gene sequencing of microcosms, with microcosm methane oxidation rates (μmol CH4 g−1 d−1). Microcosm communities suggest unknown methanotrophs may be responsible for observed rates of methane oxidation. The relative abundance of taxa associated with methanotrophy [‘Known methanotroph taxa’, verrucomicrobial methanotrophs (Verruco Mtroph), gammaproteobacterial methanotrophs (Gamma Mtroph), alphaproteobacterial methanotrophs (Alpha Mtroph)] is bifurcated from those taxa not associated with methanotrophy (‘Non-methanotroph taxa’). Error bars represent the standard error of the slope of regression obtained from observed methane oxidation within microcosms. COM, Craters of the Moon; GDS, Golden Springs; LPR, Loop Road; NGM, Ngatamariki; OKO, Orakei Korako; RTK, Rotokawa; TKA, Te Kopia; TKT, Tikitere; TOK, Tokaanu; WAM, Waimangu; WAP, Waipahihi; WHV, Whakarewarewa Village; WKT, Waikite Valley; WTV, Wairakei Thermal Valley.
Observed methane oxidation rates suggest methanotrophs are active at temperature and pH values that exceed reported values for methanotroph isolates. In this study, the fastest rates of methane oxidation were observed in microcosms incubated between 60 and 70°C (Figure 2A). Methane oxidation was also observed in six microcosms incubated at 75°C, with LPR14 being notable for the highest oxidation rate at this temperature of 3.0 μmol g−1 d−1 wet weight. To date, Methylothermus strain HB is reported as having the hottest growth temperature for a methanotroph (72°C), but is no longer extant (Bodrossy et al., 1999). Of extant strains, Methylothermus thermalis grows at 67°C (Tsubota et al., 2005), and Methylothermus subterraneous (Hirayama et al., 2011) and Methylacidiphilum fumariolicum (Pol et al., 2007) both have maximal reported growth temperatures of 65°C.
Microcosm methane oxidation across a broad pH range further suggests undescribed methanotroph activity. Methane oxidation within microcosms was observed at pH values as acidic as pH 1.5 (LPR16, 0.5 μmol g−1 d−1 wet weight) and as alkaline as pH 8.7 (WHV13, 0.5 μmol g−1 d−1 wet weight). Thermoacidophilic Methylacidiphilum species, previously isolated in New Zealand (Dunfield et al., 2007; Carere et al., 2017) and elsewhere (Pol et al., 2007; Islam et al., 2008), are capable of growth between pH 0.8 and pH 6.0. However, these species appear to be constrained to growth between 37 and 65°C (Houghton et al., 2019; Schmitz et al., 2021a) with the closely related acidophilic Methylacidimicrobium species growing to a temperature of 55°C (Picone et al., 2021b). Three acidic pH microcosms displayed substantial rates of methane oxidation at temperatures >65°C: LPR14 (75°C, pH 2.8), TKT67 (70°C, pH 2.6), and LPR16 (70°C, pH 1.5). Proteobacterial methanotrophs generally prefer neutral pH (Houghton et al., 2019), and only Methyloferula stellata has been shown to grow at < pH 4.0 (at pH 3.5–7.2, and 4–33°C) (Vorobev et al., 2011). The high-temperature methane-oxidizing microcosms in this study also included activity at pH 8.3 (70°C) and pH 8.7 (75°C), exceeding reported values for methanotroph isolates. M. subterraneous, M. fumariolicum and M. thermalis have pH maxima of 7.5, 5.8 and 7.5, respectively (Tsubota et al., 2005; Pol et al., 2007; Hirayama et al., 2011). Methylocaldum species are known to grow at ≤ pH 8.5 (Eshinimaev et al., 2004), but not at temperatures exceeding 62°C, while Methylotuvimicrobium species can grow at ≤ pH 11.0, but only to ≤45°C (Kalyuzhnaya et al., 2001). Nevertheless, methane oxidation (without growth) may be observed at temperatures exceeding those recorded for growth (Houghton, 2018), as previously demonstrated with methanogenesis (Miller et al., 1988). In these situations, growth may be limited by critical heat labile enzymes (Wiegel, 1990), or cells may use the derived energy for processes other than growth such as maintenance activities or the production of intracellular storage polymers (Linton and Cripps, 1978; Pirt, 1982).
3.2. Microcosm communities show similar diversity, but are enriched with known methanotrophs where high rates of methane oxidation are observed
A total of 8,605 unique OTUs were detected across all microcosms (mean of microcosm OTUs 921; s.d. = 570) (Supplementary Table 1). An assessment of community composition indicated that many of the microcosm communities were highly diverse (Supplementary Table 1; Figure 2B) and dominated by Archaea (Supplementary Table 1). A large number of OTUs across multiple microcosms were affiliated with taxa known only through DNA sequencing; in WHV18 up to 99% of OTUs were from uncultured Bacteria or Archaea, with more than 50% uncultured in 47 microcosms. This abundance of OTUs with no described or isolated representatives is typical of geothermal systems and highlights the need to characterize these systems (Inskeep et al., 2010; Rinke et al., 2013; Yasir et al., 2019).
Seven archaeal and bacterial genera were found in >90% of the geothermal microcosms: Methylacidiphilum (Verrucomicrobiota), Chthonomonas (Armatimonadota), Thermus (Deinococcota), Sulfolobus (Thermoproteota), Thermoplasma (Euryarchaeota), Caldimicrobium (Thermodesulfobacteriota), and Alicyclobacillus (Bacillota). However, these ubiquitous genera comprised on average < 5% of the total reads in the microcosms and, with the exception of Methylacidiphilum, are likely to represent cosmopolitan taxa present across NZ geothermal environments rather than contributing directly to methane consumption (Huber and Stetter, 2001; Reysenbach, 2001; da Costa et al., 2009; Albuquerque and da Costa, 2014; Kojima et al., 2016; Lee et al., 2016). No significant differences in alpha diversity were observed between ‘methane-oxidizing’ and ‘non-oxidizing’ microcosms (p = 0.25 for Shannon and p = 0.20 for Simpson). Despite the similar alpha diversity metrics, the community composition differed significantly, most notably two methanotrophic genera (Methylacidiphilum and Methylothermus), four other genera (Thermoplasma, Alicyclobacillus, Sulfurimonas and several species of pseudomonads) not known for a methanotrophy phenotype were significantly enriched in ‘methane-oxidizing’ communities (Table 2). We speculate that their increased abundance could be a result of positive syntrophic interactions such as increased access to C1 by-products by methylotrophic pseudomonads (Riis et al., 2003), or other metabolic by products from the heterotrophic scavengers Thermoplasma and Alicyclobacillus (Reysenbach, 2001; Karavaiko et al., 2005). Conversely, the increased abundance may be merely a co-occurrence of conducive growth conditions, e.g., the elevated atmospheric CO2 amended to the microcosms are likely to promote the growth of autotrophs such as Sulfurimonas (Han and Perner, 2015). Sulfide-oxidizing Sulfurimonas strains may also benefit from the production of hydrogen sulfide from methanethiol by Methylacidiphilum species (Lahme et al., 2020; Schmitz et al., 2022).
Table 2. Genera enriched in methane-oxidizing microcosms.
3.3. Microcosms support both aerobic and anaerobic methanotrophs
Putative aerobic and anaerobic methanotroph OTUs were detected in all geothermal microcosms, regardless of whether they oxidized methane or not, across a wide range of both temperature and pH (Supplementary Figure 3; Supplementary Table 1). All of the major clades of aerobic methanotrophs Alphaproteobacteria, Gammaproteobacteria and Verrucomicrobiales were detected across the microcosm experiment.
Within the alphaproteobacterial clade, Methylocystis OTUs were detected in low abundance within 19 microcosms, ranging from pH 4.2 to pH 7.7 (Supplementary Figure 3) and at in situ temperatures between 38.0 and 77.3°C. Currently, all cultivated Methylocystis strains have reported maximum growth temperatures of ≤53°C (Tsypenzhapova et al., 2007). Molecular-based surveys have detected pmoA sequences similar to Methylocystis and Methylocapsa in a 51°C hot spring (Zelenkina et al., 2009) and in geothermal sediments up to 76°C (Sharp et al., 2014b), while Kizilova et al. (2014) demonstrated 14CH4 bio-assimilation in hot spring samples containing Alphaproteobacteria OTUs incubated at 75°C (Kizilova et al., 2014), further supporting the possibility of thermophilic alphaproteobacterial methanotrophs within these microcosms.
In comparison, OTUs related to methanotrophs from Gammaproteobacteria (e.g., Methylothermus, Methylococcus, Methylocaldum, Methylomonas, Crenothrix, Methylobacter) were identified in 49 microcosms and were substantially more abundant than those from Alphaproteobacteria, with an average abundance of 0.7% of all reads, and a maximum of 11.2% of normalized reads from TOK7. These microcosms ranged from pH 1.5 to pH 9.2, and from 35.0°C to 85.0°C (Supplementary Figure 3). The majority of gammaproteobacterial methanotroph reads affiliated most closely to thermophilic Methylothermus strains (Tsubota et al., 2005; Hirayama et al., 2011), but the microcosms GDS1 (38.9°C, pH 7.2) and GDS2 (38.0°C, pH 7.0) also contained a small number of OTUs most similar to thermotolerant (Methylococcus or Methylocaldum) and mesophilic (Methylomonas, Crenothrix and Methylobacter) strains (Figure 3).
Figure 3. Unrooted phylogenetic tree of the putative methanotroph OTUs from Pseudomonadota. The evolutionary history was inferred using Neighbor-Joining (Saitou and Nei, 1987), and evolutionary distances were computed using the Maximum Composite Likelihood method (Tamura et al., 2004). Bootstrap values >50% are shown at the branch nodes (1000 replicates). The scale bar represents 0.02 nucleotide substitutions per site, and analyses were conducted using MEGA7 (Kumar et al., 2016).
Methylacidiphilum strains from Verrucomicrobiota were the most ubiquitous and abundant methanotrophs detected throughout this survey; detected in every microcosm with an average abundance of 3.7% (up to a maximum of 60.7% of reads, NGM89). The majority of these OTUs were most closely related to the thermoacidophilic M. infernorum (Dunfield et al., 2007). Microcosms spanned temperature (35.0–85.0°C) and pH (pH 1.5–9.2) values that exceed the currently known methane-dependent growth ranges of Methylacidiphilum and Methylacidimicrobium isolates (Islam et al., 2008; Op den Camp et al., 2009; van Teeseling et al., 2014). While some of these Methylacidiphilum OTUs are likely relic DNA, it is also plausible that these taxa are growing/persisting at elevated temperatures on non-methane substrates, e.g., simple organic acids, alcohols, alkanes, H2, H2S or methanethiol (Carere et al., 2017, 2021; Picone et al., 2021b; Schmitz et al., 2022, 2023) and/or oxidizing methane via a persistence strategy (Houghton, 2018).
Microcosm OKO2 was collected from a steam-influenced soil and showed the greatest rate of methane oxidation (17.4 μmol g−1 d−1 wet weight) observed in this study. The soil was moderately acidic (pH 4.3) with an in situ temperature of 71.6°C. Analysis of the community composition revealed very few known methanotroph OTUs. In fact, the only reads associated with a known methanotrophic phenotype were from the genus Methylacidiphilum, but these reads were at such low abundance (13 reads total) that they cannot be responsible for the observed rates of methane oxidation. In contrast, the OKO2 microcosm was dominated (65.4% of total normalized reads) by Bathyarchaeota sequences. Recently published genomes of this phylum suggest that some strains may be capable of anaerobic methanotrophy (Evans et al., 2015). Despite this microcosm being incubated aerobically, it is plausible these candidate taxa may be responsible for the observed rates of methane oxidation. Other known anaerobic methanotrophs, Methanoperedens and ANME-1, were also detected along with several known aerobic methanotrophs at low abundance within a single microcosm (GDS2).
A minority of microcosms displayed rapid rates of methane oxidation despite a low abundance of known/putative methanotrophs (Figure 2B). TKT67 oxidized 9.4 μmole methane g−1 d−1 wet weight (the second fastest rate measured) and TKT68 oxidized 7.1 μmole methane g−1 d−1 wet weight. Both microcosms were dominated by Archaea (Thermoproteota and Nitrososphaerota, respectively) and both had ~0.6% of sequences classified to Methylacidiphilum and ~ 0.1% to Methylothermus (Figure 2B; Supplementary Table 2). In addition, TOK17 oxidized 5.5 μmole methane g−1 d−1 wet weight, but only 4.1% of sequences were classified as Methylothermus and 0.2% to Methylacidiphilum (Supplementary Table 1). Although community abundance does not necessarily correlate to methanotroph activity (Wigley et al., 2022), collectively these data suggest hitherto undescribed methanotrophs may be responsible for some of the methane-oxidizing activity observed.
3.4. Metatranscriptomic analysis of enrichments reveals active community members
To investigate the ability for thermophilic microcosms to propagate within growth media, a series of enrichment experiments were performed. Thirty-one microcosms that oxidized methane (>0.5 μmole g−1 d−1) were inoculated into different growth media for methanotroph enrichment experiments (n = 71 enrichments, Supplementary Figure 4; Supplementary Table 2). Ten microcosms which displayed methane oxidation (>0.5 μmole g−1 d−1) were unable to support this activity as enrichments (Figure 2; Supplementary Table 4). A similar loss of activity at elevated temperatures was seen in enrichment microcosms from landfill cover soil (Reddy et al., 2019). In this study, at 50°C methylotrophic communities shifted from primarily Methylocaldum species to a combination of Methylobacter, Methylocystis and non-methane-oxidizing Methylophilaceae, although these were presumably dormant organisms (Reddy et al., 2019). The loss of methane oxidation capacity is likely due to poor mass-transfer of methane (e.g., diffusivity, solubility) within a liquid medium at high temperatures (Castro et al., 1995; Duan and Mao, 2006), or an inability to provide conducive growth conditions in the growth medium, e.g., rare earth elements (lanthanides) (Keltjens et al., 2014). The lanthanides are essential as a cofactor in XoxF methanol dehydrogenases, which are commonly present in the genomes of both methanotrophic and methylotrophic bacteria, often in addition to MxaF methanol dehydrogenases (Chistoserdova, 2016). XoxF is essential for growth in verrucomicrobial Methylacidiphilum (Pol et al., 2013) and Methylacidimicrobium strains (Schmitz et al., 2021b), but only a few proteobacterial Methylosinus (Kato et al., 2020) and Methylocella (Crombie, 2022) strains. In this study, lanthanides were only added to V4 medium for enrichment of Methylacidiphilum strains (Supplementary Material 1).
To investigate the activity of the methane-oxidizing microbial communities, two enrichments (GDS1 and TOK7) were selected for metatranscriptomic analysis on the basis of observed methane-oxidizing activity, the prevalence of methanotroph taxa within initial microcosms, and continued methane oxidation in the enrichments. Sample GDS1 was collected from a site at Golden Springs (pH 7.2, 38.9°C), with the microcosm initially incubated at 50°C to promote any thermophilic methanotrophic activity. The GDS1 enrichment oxidized methane at a rate of 24.5 μmoles d−1 in NMS medium (Supplementary Table 4). The 16S rRNA gene sequences were dominated by Pseudomonadota (78.3% of total reads), with 5.5% of reads classified to Methylococcaceae (Supplementary Table 3). A total of 30.0 Gb trimmed and quality checked reads from three biological replicates of GDS1 (10.0 ± 0.7 Gb/sample) were obtained. More details are in Supplementary Table 5.
For GDS1, 18 of the 20 most abundant transcripts were assigned to Methylococcus, Methylomicrobium or Methylocystis species (Figure 4A; Supplementary Table 5). The remaining transcripts predominantly encoded for proteins assigned to Alphaproteobacteria (Rhodopseudomonas, 19.6%; Aurantimonas, 9.7%; Novosphingobium, 7.5%), and Actinomycetota (Streptomyces, 6.4%). Interestingly, these transcripts were not representative of the initial 16S rRNA gene sequences obtained from the original microcosm, which were predominately Pseudomonas spp. (46.9% of total normalized reads), and with no OTU reads affiliated to Rhodopseudomonas, Aurantimonas or Novosphingobium. This may be due to dormancy as up to 80% of soil microorganisms identified through DNA sequencing are thought to be in a low state of metabolic activity at any given time (Lennon and Jones, 2011).
Figure 4. Most abundant transcripts identified within (A) GDS1 and (B) TOK7 enrichments. Transcripts with an average transcripts per million (TPM) > 2000 across three biological replicates are shown, error bars represent 1 standard deviation. Identity is based on a BLASTx search using the NCBI non-redundant (nr) database (Supplementary Tables 5, 6).
Particulate methane monooxygenase proteins comprised between 18.5 to 22.9% of all transcripts in GDS1 replicates, and multiple methanol dehydrogenase proteins comprised 4.7 to 8.4% of all transcripts. On the basis of observed transcripts, enrichments may also be capable of methylamine oxidation to formaldehyde (K15229), as observed in multiple methanotrophic strains, including Methylomonas, Methylococcus and Methylocella species (Bowman et al., 1993; Dedysh et al., 2005) and methylotrophic Methylophilaceae spp. (Chistoserdova, 2011). The genomes of several Methylacidiphilum strains also contain methylamine dehydrogenase genes (Dunfield et al., 2007; Hou et al., 2008; Op den Camp et al., 2009). However, while trimethylamine could be oxidized by trimethylamine dehydrogenase (K00317), no transcripts encoding enzymes capable of oxidizing the resulting dimethylamine were detected. Transcripts for multiple potential mechanisms for formaldehyde oxidation were expressed, including glutathione-dependent formaldehyde dehydrogenases, and tetrahydrofolate- and tetrahydromethanopterin-linked pathways (Chistoserdova, 2011). However, transcripts for a membrane-associated quinoprotein formaldehyde dehydrogenase, previously reported to be the major formaldehyde-oxidizing enzyme in Methylococcus capsulatus during pMMO expression (Zahn et al., 2001), were not detected in this study. All genes encoding enzymes required for the RuMP and CBB cycles for carbon assimilation were identified, but a complete serine cycle pathway was not detected. Two essential modules were not detected in transcripts; K00830, which represents both serine to hydroxypyruvate and glycoxylate to pyruvate; and K08692, which represents malate-CoA ligase.
Consistent with other methane-oxidizing communities at the oxic/anoxic interface (Wigley et al., 2022), GDS1 enrichments were putatively capable of assimilatory and (partial) dissimilatory nitrate reduction (NMS medium contains 0.6 g/L−1 NO3). Multiple pathways were identified in GDS1 transcriptomes for both dissimilatory and assimilatory nitrate reduction to ammonia using genes encoding the NarGHJI (K00370), NapAB (K02567) and NasAB (K00372) forms of nitrate reductases, although these were not identified as from methanotrophs. Nitrite could be converted to ammonia by various forms of nitrite reductase, using either NADH (NirBD, K00362), cytochrome-c552 (NrfAH, K03385), NAD(P)H (NIT-6, K17877) or ferredoxin (NirA, K00366) as co-factors. However, genes encoding NO-forming nitrite reductase (NirK, K00368, or NirS, K15864) were not detected. One possible explanation could be that methanotrophy is mediated via a similar nitrate-dependent mechanism to the mesophilic genera Methylomirabilis (phylum NC10) (Guerrero-Cruz et al., 2019) via yet-to-be-described taxa when oxygen at the oxic-anoxic interface becomes limited. GDS1 also displayed evidence for nitrogen reduction to ammonia via the activity of a molybdenum-iron nitrogenase (NifDKH, K02586).
Sample TOK7 was collected from a clay soil site at Tokaanu (pH 6.8, 63.7°C, Supplementary Table 1). The TOK7 enrichment suspended in ANMS medium oxidized CH4 at a rate of 3.5 μmoles d−1 (Supplementary Table 4) while incubated at 60°C. The archaeal phylum Bathyarchaeota, which potentially includes anaerobic methanotrophs (Evans et al., 2015), was the most abundant phylum detected via 16S rRNA gene sequences (27.9%). Also noteworthy was an abundance (11.2%) of sequences classified into the gammaproteobacterial family of thermophilic methanotrophs Methylothermaceae (Houghton et al., 2019; Supplementary Table 3).
A total of 31.7 Gb trimmed and quality checked reads from three biological replicates of TOK7 (10.6 ± 1.3 Gb/sample) were obtained. More details are in Supplementary Table 6. The most highly expressed genes were related to known methanotroph taxa (Figure 4B). In TOK7 samples, 13 of the 28 most abundant transcripts were assigned to strains from Methylothermaceae, and another six to Methylococcaceae (Supplementary Table 6). However, there were only 109 transcripts in total from methanotrophs (0.2%), including 16 from Methylocystis or Methylosinus. No transcripts from Bathyarchaeota were identified. The remaining transcripts included genes encoding for proteins classified as Pseudomonas (Gammaproteobacteria, 7.8%), Anaeromyxobacter (Deltaproteobacteria, 5.2%), Streptomyces (Actinomycetota, 5.2%), Mycobacterium (Actinomycetota, 5.0%), and Rhodopseudomonas (Alphaproteobacteria, 5.0%). Similar to GDS1 samples, this did not reflect the original community 16S rRNA gene sequences, which contained only six Pseudomonas OTUs from 149,400 normalized reads, and none of the other aforementioned genera. Although transcriptome analysis specifically identifies ‘active’ members in a mixed community, taxonomic assignments based on community protein transcripts should also be assessed with caution, as there is known bias within the NCBI database for Pseudomonadota, Actinomycetota, Bacteroidota, and Bacillota (Rinke et al., 2013). Furthermore, high conservation of protein domains across taxa (or variability within families) may further complicate taxonomic placement of transcripts (Mount, 2004).
Consistent with the rapid rates of methane oxidation observed within TOK7 enrichments (Supplementary Figure 5), genes involved in one-carbon metabolism were highly expressed. The most highly expressed gene (5.8% of transcripts) encoded a methylotrophy-associated pyrroloquinoline quinone (PQQ) precursor peptide, PqqA (Chistoserdova et al., 2003). PQQ-dependent dehydrogenases (methanol/ethanol family, Methylohalobius) were also abundant (0.4–0.8% of transcripts) (Supplementary Table 6). Genes encoding the particulate form of methane monooxygenase (Methylomicrobium/Methylococcus/Methylocystis) comprised between 0.9 and 1.8% of transcripts (Supplementary Table 6). These did not match with the methanotrophs identified during 16S rRNA sequencing (primarily Methylothermus, with nine Methylocystis reads). Formaldehyde could be assimilated and detoxified via S-(hydroxymethyl) glutathione synthase (K03396), a formaldehyde activating enzyme (K10713) for the tetrahydromethanopterin pathway, or by 3-hexulose-6-phosphate synthase (K08093) for the RuMP pathway. Carbon assimilation within the TOK7 community may also be achieved via the CBB cycle, as all genes for this pathway were identified from the assembly. In contrast, the serine cycle for carbon assimilation was not complete, as genes encoding malate-CoA ligase (K08692), malyl-CoA lyase (K08691) and an alanine/serine transaminase (K00830) were all absent from the transcripts. No genes encoding enzymes for the oxidation of mono-, di- or trimethylamine were found in the TOK7 transcripts, and the only evidence for methanogenesis were genes encoding two subunits (of six) of heterodisulfide reductase.
TOK7 enrichments were putatively capable of assimilatory and (partial) dissimilatory nitrate reduction, although these were not associated with methanotrophs (ANMS medium contains 0.15 g/L−1 NO3). The TOK7 community expressed narGHI genes for dissimilatory nitrate reduction to nitrite and could potentially convert this nitrite to ammonia via NirBD (K00362) or NfrAH (K03385) nitrite reductases. The community was also capable of assimilatory nitrate reduction using NAD(P)H-dependent nitrate reductase (K10534) or the NasAB (K00372) nitrate reductase, and finally converting the resultant nitrite to ammonia using nitrite reductases NIT-6 (K17877) or NirA (K00366), which use NAD(P)H or ferredoxin, as co-factors, respectively (Exley et al., 1993; Takahashi et al., 2001). No genes for the fixation of nitrogen via nitrogenases nifDKH or anfHDGK were identified in TOK7.
4. Conclusion
The data presented in this study broadly supports the premise that thermophilic CH4 oxidation is widespread in New Zealand geothermal fields, with more than 50% of the methanotroph-targeting microcosms exhibiting methane oxidation at temperatures greater than 38°C and as hot as 75°C. Methane oxidation rates up to 17.4 μmol g−1 d−1 wet weight indicate that this is an important sink for CH4 emitted from geological sources. Microbial communities within microcosms were highly diverse, with few ubiquitous genera, but both known aerobic methanotrophs and putative anaerobic methanotrophs were identified via 16S rRNA gene sequencing of microcosms. To investigate this further, we attempted to enrich methanotrophic consortia from 32 thermophilic methane-oxidizing microcosms (< 75°C) but were unable to maintain many of the cultures. For enrichments that did oxidize methane, observed rates correlated poorly to the original microcosm, or to the presence of known methanotroph OTUs (as predicted via 16S rRNA gene amplicon sequencing). The cessation of CH4 oxidation in these enrichments was possibly due to the poor solubility of CH4/O2 at high temperatures. Consequently, poor mass transfer of these gases into the growth medium may be hampering attempts to cultivate and isolate novel thermophilic methane-oxidizing bacteria using traditional batch cultivation techniques. While we present evidence for the widespread occurrence of methane oxidation in New Zealand geothermal environments, there is also evidence of the presence of methanotrophs in high temperature environments globally (Houghton et al., 2019; Wang et al., 2022; Zhou et al., 2022).
Finally, metatranscriptomic analysis of enrichments actively oxidizing methane showed genes involved in carbon metabolism were highly transcribed. At mesophilic temperatures (GDS1, 37°C), these genes were most similar to Methylococcaceae species, while at hotter temperatures (TOK7, 46°C), the majority of transcripts were associated with the moderately thermophilic Methylothermaceae. The important role of methanotrophs in nitrogen metabolism under oxic conditions (Hoefman et al., 2014) was also highlighted; both GDS1 and TOK7 communities were carrying out assimilatory and dissimilatory nitrate reduction, although a denitrification pathway was not complete in either set of transcripts. Collectively, this study expands knowledge of thermophilic methanotrophy in geothermal areas and suggests that hitherto unidentified methanotrophs may be responsible for some of this activity. A closer examination of the mechanisms involved in control of expression of methane and nitrogen metabolism pathways will augment our understanding of how methanotrophs survive and flourish in these ecosystems.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA766707, https://www.ncbi.nlm.nih.gov/, PRJNA546003, https://doi.org/10.5281/zenodo.5535083, 5535083.
Author contributions
KH: Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. CC: Conceptualization, Funding acquisition, Methodology, Visualization, Writing – review & editing. MS: Conceptualization, Funding acquisition, Methodology, Supervision, Visualization, Writing – review & editing. IM: Conceptualization, Methodology, Supervision, Writing – review & editing.
Funding
This work was supported by the GNS Science New Zealand’s Geothermal Future Research Programme, the Royal Society of New Zealand (Marsden Grant GNS1601 to CC) and a Freemasons New Zealand Postgraduate Scholarship to KH.
Acknowledgments
The authors thank Jean Power and Dave Evans for assistance with sample collection. The authors wish to acknowledge all landowners and Māori collaborators for access and support of this research. Mana whenua (customary rights) is acknowledged for all data generated arising from geothermal features within rohe of iwi (tribal territories).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1253773/full#supplementary-material
References
Afgan, E., Baker, D., van den Beek, M., Blankenberg, D., Bouvier, D., Čech, M., et al. (2016). The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update. Nucleic Acids Res. 44, W3–W10. doi: 10.1093/nar/gkw343
Albuquerque, L., and da Costa, M. S. (2014). “The family Thermaceae” in The prokaryotes – other major lineages of bacteria and the Archaea. eds. E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson. 4th ed (Berlin: Springer-Verlag)
Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. doi: 10.1016/S0022-2836(05)80360-2
Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389
Andrews, S. (2010). FastQC: a quality control tool for high throughput sequence data. Cambridge Babraham Institute.
Awala, S. I., Gwak, J. H., Kim, Y. M., Kim, S. J., Strazzulli, A., Dunfield, P. F., et al. (2021). Verrucomicrobial methanotrophs grow on diverse C3 compounds and use a homolog of particulate methane monooxygenase to oxidize acetone. ISME J. 15, 3636–3647. doi: 10.1038/s41396-021-01037-2
Besemer, J., and Borodovsky, M. (1999). Heuristic approach to deriving models for gene finding. Nucleic Acids Res. 27, 3911–3920. doi: 10.1093/nar/27.19.3911
Bodrossy, L., Holmes, E. M., Holmes, A. J., Kovács, K. L., and Murrell, J. C. (1997). Analysis of 16S rRNA and methane monooxygenase gene sequences reveals a novel group of thermotolerant and thermophilic methanotrophs, Methylocaldum gen. nov. Arch. Microbiol. 168, 493–503. doi: 10.1007/s002030050527
Bodrossy, L., Kovács, K. L., McDonald, I. R., and Murrell, J. C. (1999). A novel thermophilic methane-oxidising γ-Proteobacterium. FEMS Microbiol. Lett. 170, 335–341.
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Bornemann, M., Bussmann, I., Tichy, L., Deutzmann, J., Schink, B., and Pester, M. (2016). Methane release from sediment seeps to the atmosphere is counteracted by highly active Methylococcaceae in the water column of deep oligotrophic Lake Constance. FEMS Microbiol. Ecol. 92:fiw123. doi: 10.1093/femsec/fiw123
Bowman, J. P., Sly, L. I., Nichols, P. D., and Hayward, A. C. (1993). Revised taxonomy of the methanotrophs: description of Methylobacter gen. nov., emendation of Methylococcus, validation of Methylosinus and Methylocystis species, and a proposal that the family Methylococcaceae includes only the group I methanotrophs. Int. J. Syst. Bacteriol. 44:375.
Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. BMC Bioinformat. 10:421. doi: 10.1186/1471-2105-10-421
Canadell, J. G., Monteiro, P. M. S., Costa, M. H., Cotrim da Cunha, L., Cox, P. M., Eliseev, A. V., et al. (2021). “Global carbon and other biogeochemical cycles and feedbacks” in Climate change 2021: the physical science basis. Contribution of working group I to the sixth assessment report of the intergovernmental panel on climate change. eds. V. Masson-Delmotte, P. Zhai, A. Pirani, S. L. Connors, C. Péan, and S. Berger, et al. (Cambridge: IPCC)
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Lozupone, C. A., Turnbaugh, P. J., et al. (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. U. S. A. 108, 4516–4522. doi: 10.1073/pnas.1000080107
Carere, C. R., Hards, K., Houghton, K. M., Power, J. F., McDonald, B., Collet, C., et al. (2017). Mixotrophy drives niche expansion of verrucomicrobial methanotrophs. ISME J. 11, 2599–2610. doi: 10.1038/ismej.2017.112
Carere, C. R., Hards, K., Wigley, K., Carman, L., Houghton, K. M., Cook, G. M., et al. (2021). Growth on formic acid is dependent on intracellular pH homeostasis for the thermoacidophilic methanotroph Methylacidiphilum sp. RTK17.1. Front. Microbiol. 12:651744. doi: 10.3389/fmicb.2021.651744
Castaldi, S., and Tedesco, D. (2005). Methane production and consumption in an active volcanic environment of southern Italy. Chemosphere 58, 131–139. doi: 10.1016/j.chemosphere.2004.08.023
Castro, M. S., Steudler, P. A., Melillo, J. M., Aber, J. D., and Bowden, R. D. (1995). Factors controlling atmospheric methane consumption by temperate forest soils. Glob. Biogeochem. Cycles 9, 1–10. doi: 10.1029/94GB02651
Chao, A. (1984). Nonparametric estimation of the number of classes in a population. Scand. J. Stat. 11, 265–270.
Chistoserdova, L. (2011). Modularity of methylotrophy, revisited. Environ. Microbiol. 13, 2603–2622. doi: 10.1111/j.1462-2920.2011.02464.x
Chistoserdova, L. (2016). Lanthanides: new life metals? World J. Microbiol. Biotechnol. 32:138. doi: 10.1007/s11274-016-2088-2
Chistoserdova, L., Chen, S. W., Lapidus, A., and Lidstrom, M. E. (2003). Methylotrophy in Methylobacterium extorquens AM1 from a genomic point of view. J. Bacteriol. 185, 2980–2987. doi: 10.1128/JB.185.10.2980-2987.2003
Crombie, A. T. (2022). The effect of lanthanum on growth and gene expression in a facultative methanotroph. Environ. Microbiol. 24, 596–613. doi: 10.1111/1462-2920.15685
da Costa, M. S., Rainey, F. A., and Albuquerque, L. (2009). “Alicyclobacillus” in Bergey's manual of systematic bacteriology. eds. P. De Vos, G. M. Garrity, D. T. Jones, N. R. Kreig, W. Ludwig, and F. A. Rainey, et al. 2nd ed (London: Springer)
Danilova, O. V., Kulichevskaya, I. S., Rozova, O. N., Detkova, E. N., Bodelier, P. L. E., Trotsenko, Y. A., et al. (2013). Methylomonas paludis sp. nov., the first acid-tolerant member of the genus Methylomonas, from an acidic wetland. Int. J. Syst. Evol. Microbiol. 63, 2282–2289. doi: 10.1099/ijs.0.045658-0
Dedysh, S. N., Knief, C., and Dunfield, P. F. (2005). Methylocella species are facultatively methanotrophic. J. Bacteriol. 187, 4665–4670. doi: 10.1128/JB.187.13.4665-4670.2005
Duan, Z., and Mao, S. (2006). A thermodynamic model for calculating methane solubility, density and gas phase composition of methane-bearing aqueous fluids from 273 to 523 K and from 1 to 2000 bar. Geochim. Cosmochim. Acta 70, 3369–3386. doi: 10.1016/j.gca.2006.03.018
Dunfield, P. F., Yuryev, A., Senin, P., Smirnova, A. V., Stott, M. B., Hou, S., et al. (2007). Methane oxidation by an extremely acidophilic bacterium of the phylum Verrucomicrobia. Nature 450, 879–882. doi: 10.1038/nature06411
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Eshinimaev, B. T., Medvedkova, K. A., Khmelenina, V. N., Suzina, N. E., Osipov, G. A., Lysenko, A. M., et al. (2004). New thermophilic methanotrophs of the genus Methylocaldum. Microbiology 73, 448–456. doi: 10.1023/B:MICI.0000036991.31677.13
Esson, K. C., Lin, X., Kumaresan, D., Chanton, J. P., Murrell, J. C., and Kostka, J. E. (2016). Alpha-and gammaproteobacterial methanotrophs codominate the active methane-oxidizing communities in an acidic boreal peat bog. Appl. Environ. Microbiol. 82, 2363–2371. doi: 10.1128/AEM.03640-15
Evans, P. N., Parks, D. H., Chadwick, G. L., Robbins, S. J., Orphan, V. J., Golding, S. D., et al. (2015). Methane metabolism in the archaeal phylum Bathyarchaeota revealed by genome-centric metagenomics. Science 350, 434–438. doi: 10.1126/science.aac7745
Exley, G. E., Colandene, J. D., and Garrett, R. H. (1993). Molecular cloning, characterization, and nucleotide sequence of nit-6, the structural gene for nitrite reductase in Neurospora crassa. J. Bacteriol. 175, 2379–2392. doi: 10.1128/jb.175.8.2379-2392.1993
Forster, P., Storelvmo, T., Armour, K., Collins, W., Dufresne, J. L., Frame, D., et al. (2021). “The earth’s energy budget, climate feedbacks, and climate sensitivity” in Climate change 2021: The physical science basis. Contribution of working group I to the sixth assessment report of the intergovernmental panel on climate change. eds. V. Masson-Delmotte, P. Zhai, A. Pirani, S. L. Connors, C. Péan, and S. Berger, et al. (Cambridge: IPCC)
Gagliano, A. L., D'Alessandro, W., Tagliavia, M., Parello, F., and Quatrini, P. (2014). Methanotrophic activity and diversity of methanotrophs in volcanic geothermal soils at Pantelleria (Italy). Biogeosciences 11, 5865–5875. doi: 10.5194/bg-11-5865-2014
Geymonat, E., Ferrando, L., and Tarlera, S. E. (2011). Methylogaea oryzae gen. nov., sp. nov., a mesophilic methanotroph isolated from a rice paddy field. Int. J. Syst. Evol. Microbiol. 61, 2568–2572. doi: 10.1099/ijs.0.028274-0
Giggenbach, W. F. (1995). Variations in the chemical and isotopic composition of fluids discharged from the Taupo Volcanic Zone, New Zealand. J. Volcanol. Geotherm. Res. 68, 89–116. doi: 10.1016/0377-0273(95)00009-J
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Guerrero-Cruz, S., Stultiens, K., van Kessel, M. A. H. J., Versantvoort, W., Jetten, M. S., den Camp, H. J. M. O., et al. (2019). Key physiology of a nitrite-dependent methane-oxidizing enrichment culture. Appl. Environ. Microbiol. 85:e00124-19. doi: 10.1128/AEM.00124-19
Han, Y., and Perner, M. (2015). The globally widespread genus Sulfurimonas: versatile energy metabolisms and adaptations to redox clines. Front. Microbiol. 6:989. doi: 10.3389/fmicb.2015.00989
Hanson, R. S., and Hanson, T. E. (1996). Methanotrophic bacteria. Microbiol. Rev. 60, 439–471. doi: 10.1128/mr.60.2.439-471.1996
Håvelsrud, O. E., Haverkamp, T. H. A., Kristensen, T., Jakobsen, K. S., and Rike, A. G. (2011). A metagenomic study of methanotrophic microorganisms in Coal Oil Point seep sediments. BMC Microbiol. 11. doi: 10.1186/1471-2180-11-221
Henneberger, R., Chiri, E., Bodelier, P. E. L., Frenzel, P., Lüke, C., and Schroth, M. H. (2014). Field-scale tracking of active methane-oxidizing communities in a landfill-cover soil reveals spatial and seasonal variability. Environ. Microbiol. 17, 1721–1737. doi: 10.1111/1462-2920.12617
Hirayama, H., Suzuki, Y., Abe, M., Miyazaki, M., Makita, H., Inagaki, F., et al. (2011). Methylothermus subterraneus sp. nov., a moderately thermophilic methanotroph isolated from a terrestrial subsurface hot aquifer. Int. J. Syst. Evol. Microbiol. 61, 2646–2653. doi: 10.1099/ijs.0.028092-0
Hoefman, S., van der Ha, D., Boon, N., Vandamme, P., De Vos, P., and Heylen, K. (2014). Niche differentiation in nitrogen metabolism among methanotrophs within an operational taxonomic unit. BMC Microbiol. 14:83. doi: 10.1186/1471-2180-14-83
Hou, S., Makarova, K. S., Saw, J. H. W., Senin, P., Ly, B. V., Zhou, Z., et al. (2008). Complete genome sequence of the extremely acidophilic methanotroph isolate V4, Methylacidiphilum infernorum, a representative of the bacterial phylum Verrucomicrobia. Biol. Direct 3. doi: 10.1186/1745-6150-3-26
Houghton, K. M. (2018). Thermophilic Methanotrophy in the Taupō volcanic zone. PhD, University of Waikato, New Zealand
Houghton, K. M., Carere, C. R., Stott, M. B., and McDonald, I. R. (2019). Thermophilic methanotrophs: in hot pursuit. FEMS Microbiol. Ecol. 95. doi: 10.1093/femsec/fiz125
Houghton, K. M., and Stewart, L. C. (2019). Temperature-gradient incubation isolates multiple competitive species from a single environmental sample. Access Microbiol. 2, –acmi000081. doi: 10.1099/acmi.0.000081
Huang, Y., Niu, B., Gao, Y., Fu, L., and Li, W. (2010). CD-HIT suite: a web server for clustering and comparing biological sequences. Bioinformatics 26, 680–682. doi: 10.1093/bioinformatics/btq003
Huber, H., and Stetter, K. O. (2001). “The order Sulfolobales” in Bergey's manual of systematic bacteriology. eds. D. R. Boone and R. W. Castenholz. 2nd ed (New York: Springer)
Inskeep, W. P., Rusch, D. B., Jay, Z. J., Herrgard, M. J., Kozubal, M. A., Richardson, T. H., et al. (2010). Metagenomes from high-temperature chemotrophic systems reveal geochemical controls on microbial community structure and function. PLoS One 5:e9773. doi: 10.1371/journal.pone.0009773
Islam, T., Jensen, S., Reigstad, L. J., Larsen, Ø., and Birkeland, N. K. (2008). Methane oxidation at 55°C and pH 2 by a thermoacidophilic bacterium belonging to the Verrucomicrobia phylum. Proc. Natl. Acad. Sci. U. S. A. 105, 300–304. doi: 10.1073/pnas.0704162105
Jiang, H., Chen, Y., Jiang, P., Zhang, C., Smith, T. J., Murrell, J. C., et al. (2010). Methanotrophs: multifunctional bacteria with promising applications in environmental bioengineering. Biochem. Eng. J. 49, 277–288. doi: 10.1016/j.bej.2010.01.003
Kalyuzhnaya, M. G., Khmelenina, V. N., Eshinimaev, B. T., Suzina, N. E., Nikitin, D., Solonin, A., et al. (2001). Taxonomic characterization of new alkaliphilic and alkalitolerant methanotrophs from soda lakes of the Southeastern Transbaikal region and description of Methylomicrobium buryatense sp. nov. Syst. Appl. Microbiol. 24, 166–176. doi: 10.1078/0723-2020-00028
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2016). KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44, D457–D462. doi: 10.1093/nar/gkv1070
Karavaiko, G. I., Bogdanova, T. I., Tourova, T. P., Kondrateva, T. F., Tsaplina, I. A., Egorova, M. A., et al. (2005). Reclassification of ‘Sulfobacillus thermosulfidooxidans subsp. thermotolerans’ strain K1 as Alicyclobacillus tolerans sp. nov. and Sulfobacillus disulfidooxidans Dufresne et al. 1996 as Alicyclobacillus disulfidooxidans comb. nov., and emended description of the genus Alicyclobacillus). Int. J. Syst. Evol. Microbiol. 55, 941–947. doi: 10.1099/ijs.0.63300-0
Kato, S., Takashino, M., Igarashi, K., and Kitagawa, W. (2020). Isolation and genomic characterization of a proteobacterial methanotroph requiring lanthanides. Microbes Environ. 35:n/a. doi: 10.1264/jsme2.ME19128
Keltjens, J. T., Pol, A., Reimann, J., and Op den Camp, H. J. M. (2014). PQQ-dependent methanol dehydrogenases: rare-earth elements make a difference. Appl. Microbiol. Biotechnol. 98, 6163–6183. doi: 10.1007/s00253-014-5766-8
Kizilova, A. K., Sukhacheva, M. V., Pimenov, N. V., Yurkov, A. M., and Kravchenko, I. K. (2014). Methane oxidation activity and diversity of aerobic methanotrophs in pH-neutral and semi-neutral thermal springs of the Kunashir Island, Russian Far East. Extremophiles 18, 207–218. doi: 10.1007/s00792-013-0603-z
Kojima, H., Umezawa, K., and Fukui, M. (2016). Caldimicrobium thiodismutans sp. nov., a sulfur-disproportionating bacterium isolated from a hot spring, and emended description of the genus Caldimicrobium. Int. J. Syst. Evol. Microbiol. 66, 1828–1831. doi: 10.1099/ijsem.0.000947
Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054
Lahme, S., Callbeck, C. M., Eland, L. E., Wipat, A., Enning, D., Head, I. M., et al. (2020). Comparison of sulfide-oxidizing Sulfurimonas strains reveals a new mode of thiosulfate formation in subsurface environments. Environ. Microbiol. 22, 1784–1800. doi: 10.1111/1462-2920.14894
Lee, K. C., Stott, M. B., Dunfield, P. F., Huttenhower, C., McDonald, I. R., and Morgan, X. C. (2016). The Chthonomonas calidirosea genome is highly conserved across geographic locations and distinct chemical and microbial environments in New Zealand's Taupō Volcanic Zone. Appl. Environ. Microbiol. 82, 3572–3581. doi: 10.1128/AEM.00139-16
Lennon, J. T., and Jones, S. E. (2011). Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat. Rev. Microbiol. 9, 119–130. doi: 10.1038/nrmicro2504
Lesniewski, R. A., Jain, S., Anantharaman, K., Schloss, P. D., and Dick, G. J. (2012). The metatranscriptome of a deep-sea hydrothermal plume is dominated by water column methanotrophs and lithotrophs. ISME J. 6, 2257–2268. doi: 10.1038/ismej.2012.63
Liebner, S., and Svenning, M. M. (2013). Environmental transcription of mmoX by methane-oxidizing Proteobacteria in a subarctic palsa peatland. Appl. Environ. Microbiol. 79, 701–706. doi: 10.1128/AEM.02292-12
Linton, J. D., and Cripps, R. E. (1978). The occurrence and identification of intracellular polyglucose storage granules in Methylococcus NCIB 11083 grown in chemostat culture on methane. Arch. Microbiol. 117, 41–48. doi: 10.1007/BF00689349
Liu, J., Wang, S., Qin, T., Li, N., Niu, Y., Li, D., et al. (2015). Whole transcriptome analysis of Penicillium digitatum strains treatmented with prochloraz reveals their drug-resistant mechanisms. BMC Genomics 16:855. doi: 10.1186/s12864-015-2043-x
McDonald, A. I. (2014). Handbook of biological statistics, Baltimore, Maryland, Sparky House Publishing.
McDonald, I. R., Kenna, E. M., and Murrell, J. C. (1995). Detection of methanotrophic bacteria in environmental samples with the PCR. Appl. Environ. Microbiol. 61, 116–121. doi: 10.1128/aem.61.1.116-121.1995
Miller, J. F., Shah, N. N., Nelson, C. M., Ludlow, J. M., and Clark, D. S. (1988). Pressure and temperature effects on growth and methane production of the extreme thermophile Methanococcus jannaschii. Appl. Environ. Microbiol. 54, 3039–3042. doi: 10.1128/aem.54.12.3039-3042.1988
Mount, D. W. (2004). Bioinformatics: Sequence and genome analysis, New York, Cold Spring Harbour Laboratory Press.
Nguyen, N.-L., Yu, W.-J., Gwak, J.-H., Kim, S.-J., Park, S.-J., Herbold, C. W., et al. (2018). Genomic insights into the acid adaptation of novel methanotrophs enriched from acidic forest soils. Front. Microbiol. 9:1982. doi: 10.3389/fmicb.2018.01982
Olins, H. C., Rogers, D. R., Preston, C., Ussler, W., Pargett, D., Jensen, S., et al. (2017). Co-registered geochemistry and metatranscriptomics reveal unexpected distributions of microbial activity within a hydrothermal vent field. Front. Microbiol. 8:1042. doi: 10.3389/fmicb.2017.01042
Op den Camp, H. J. M., Islam, T., Stott, M. B., Harhangi, H. R., Hynes, A., Schouten, S., et al. (2009). Environmental, genomic and taxonomic perspectives on methanotrophic Verrucomicrobia. Environ. Microbiol. Rep. 1, 293–306. doi: 10.1111/j.1758-2229.2009.00022.x
Patro, R., Duggal, G., Love, M. I., Irizarry, R. A., and Kingsford, C. (2017). Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419. doi: 10.1038/nmeth.4197
Picone, N., Blom, P., Hogendoorn, C., Frank, J., van Alen, T., Pol, A., et al. (2021a). Metagenome assembled genome of a novel verrucomicrobial methanotroph from Pantelleria Island. Front. Microbiol. 12:666929. doi: 10.3389/fmicb.2021.666929
Picone, N., Blom, P., Wallenius, A. J., Hogendoorn, C., Mesman, R., Cremers, G., et al. (2021b). Methylacidimicrobium thermophilum AP8, a novel methane- and hydrogen-oxidizing bacterium isolated from volcanic soil on Pantelleria Island, Italy. Front. Microbiol. 12:637762. doi: 10.3389/fmicb.2021.637762
Pirt, S. J. (1982). Maintenance energy: a general model for energy-limited and energy-sufficient growth. Arch. Microbiol. 133, 300–302. doi: 10.1007/BF00521294
Pol, A., Barends, T. R. M., Dietl, A., Khadem, A. F., Eygensteyn, J., Jetten, M. S. M., et al. (2013). Rare earth metals are essential for methanotrophic life in volcanic mudpots. Environ. Microbiol. 16, 255–264. doi: 10.1111/1462-2920.12249
Pol, A., Heijmans, K., Harhangi, H. R., Tedesco, D., Jetten, M. S. M., and Op Den Camp, H. J. M. (2007). Methanotrophy below pH 1 by a new Verrucomicrobia species. Nature 450, 874–878. doi: 10.1038/nature06222
Power, J. F., Carere, C. R., Lee, C. K., Wakerley, G. L. J., Evans, D. W., Button, M., et al. (2018). Microbial biogeography of 925 geothermal springs in New Zealand. Nat. Commun. 9:2876. doi: 10.1038/s41467-018-05020-y
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Reddy, K. R., Rai, R. K., Green, S. J., and Chetri, J. K. (2019). Effect of temperature on methane oxidation and community composition in landfill cover soil. J. Ind. Microbiol. Biotechnol. 46, 1283–1295. doi: 10.1007/s10295-019-02217-y
Reysenbach, A. L. (2001). “Thermoplasmata class. nov” in Bergey's manual of systematic bacteriology. eds. D. R. Boone and R. W. Castenholz. 2nd ed (New York: Springer)
Riis, V., Miethe, D., and Babel, W. (2003). Formate-stimulated oxidation of methanol by Pseudomonas putida 9816. Biosci. Biotechnol. Biochem. 67, 684–690. doi: 10.1271/bbb.67.684
Rinke, C., Schwientek, P., Sczyrba, A., Ivanova, N. N., Anderson, I. J., Cheng, J.-F., et al. (2013). Insights into the phylogeny and coding potential of microbial dark matter. Nature 499, 431–437. doi: 10.1038/nature12352
Saitou, N., and Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425.
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Schmitz, R. A., Mohammadi, S. S., van Erven, T., Berben, T., Jetten, M. S. M., Pol, A., et al. (2022). Methanethiol consumption and hydrogen sulfide production by the thermoacidophilic methanotroph Methylacidiphilum fumariolicum SolV. Front. Microbiol. 13:857442. doi: 10.3389/fmicb.2022.857442
Schmitz, R. A., Peeters, S. H., Mohammadi, S. S., Berben, T., van Erven, T., Iosif, C. A., et al. (2023). Simultaneous sulfide and methane oxidation by an extremophile. Nat. Commun. 14:2974. doi: 10.1038/s41467-023-38699-9
Schmitz, R. A., Peeters, S. H., Versantvoort, W., Picone, N., Pol, A., Jetten, M. S. M., et al. (2021a). Verrucomicrobial methanotrophs: ecophysiology of metabolically versatile acidophiles. FEMS Microbiol. Rev. 45:fuab007. doi: 10.1093/femsre/fuab007
Schmitz, R. A., Picone, N., Singer, H., Dietl, A., Seifert, K. A., Pol, A., et al. (2021b). Neodymium as metal cofactor for biological methanol oxidation: structure and kinetics of an XoxF1-type methanol dehydrogenase. MBio 12:e0170821. doi: 10.1128/mBio.01708-21
Shannon, C. E., and Weaver, W. (1963). The mathematical theory of communication Champaign University of Illinois Press.
Sharp, C. E., Martinez-Lorenzo, A., Brady, A. L., Grasby, S. E., and Dunfield, P. F. (2014a). Methanotrophic bacteria in warm geothermal spring sediments identified using stable isotope probing. FEMS Microbiol. Ecol. 90, 92–102. doi: 10.1111/1574-6941.12375
Sharp, C. E., Smirnova, A. V., Graham, J. M., Stott, M. B., Khadka, R., Moore, T. R., et al. (2014b). Distribution and diversity of Verrucomicrobia methanotrophs in geothermal and acidic environments. Environ. Microbiol. 16, 1867–1878. doi: 10.1111/1462-2920.12454
Sharp, C. E., Stott, M. B., and Dunfield, P. F. (2012). Detection of autotrophic verrucomicrobial methanotrophs in a geothermal environment using stable isotope probing. Front. Microbiol. 3:303. doi: 10.3389/fmicb.2012.00303
Takahashi, M., Sasaki, Y., Ida, S., and Morikawa, H. (2001). Nitrite reductase gene enrichment improves assimilation of NO2 in Arabidopsis. Plant Physiol. 126, 731–741. doi: 10.1104/pp.126.2.731
Tamura, K., Nei, M., and Kumar, S. (2004). Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc. Natl. Acad. Sci. U. S. A. 101, 11030–11035. doi: 10.1073/pnas.0404206101
Taubert, M., Grob, C., Crombie, A., Howat, A. M., Burns, O. J., Weber, M., et al. (2019). Communal metabolism by Methylococcaceae and Methylophilaceae is driving rapid aerobic methane oxidation in sediments of a shallow seep near Elba, Italy. Environ. Microbiol. 21, 3780–3795. doi: 10.1111/1462-2920.14728
Tsubota, J., Eshinimaev, B. T., Khmelenina, V. N., and Trotsenko, Y. A. (2005). Methylothermus thermalis gen. nov., sp. nov., a novel moderately thermophilic obligate methanotroph from a hot spring in Japan. Int. J. Syst. Evol. Microbiol. 55, 1877–1884. doi: 10.1099/ijs.0.63691-0
Tsypenzhapova, I. S., Eshinimaev, B. T., Khmelenina, V. N., Osipov, G. A., and Trotsenko, I. A. (2007). A new thermotolerant aerobic methanotroph from a thermal spring in Buryatia. Microbiology 76, 118–121. doi: 10.1134/S0026261707010171
van Teeseling, M. C. F., Pol, A., Harhangi, H. R., van der Zwart, S., Jetten, M. S. M., Op den Camp, H. J. M., et al. (2014). Expanding the verrucomicrobial methanotrophic world: description of three novel species of Methylacidimicrobium gen. nov. Appl. Environ. Microbiol. 80, 6782–6791. doi: 10.1128/AEM.01838-14
Vorobev, A. V., Baani, M., Doronina, N. V., Brady, A. L., Liesack, W., Dunfield, P. F., et al. (2011). Methyloferula stellata gen. nov., sp. nov., an acidophilic, obligately methanotrophic bacterium that possesses only a soluble methane monooxygenase. Int. J. Syst. Evol. Microbiol. 61, 2456–2463. doi: 10.1099/ijs.0.028118-0
Wang, Y., Bi, H. Y., Chen, H. G., Zheng, P. F., Zhou, Y. L., and Li, J. T. (2022). Metagenomics reveals dominant unusual sulfur oxidizers inhabiting active hydrothermal chimneys from the Southwest Indian Ridge. Front. Microbiol. 13:861795. doi: 10.3389/fmicb.2022.861795
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, Y., Hayatsu, M., and Fujii, T. (2012). Extraction of bacterial RNA from soil: challenges and solutions. Microbes Environ. 27, 111–121. doi: 10.1264/jsme2.ME11304
Wiegel, J. (1990). Temperature spans for growth – hypothesis and discussion. FEMS Microbiol. Lett. 75, 155–169. doi: 10.1111/j.1574-6968.1990.tb04092.x
Wigley, K., Egbadon, E., Carere, C. R., Weaver, L., Baronian, K., Burbery, L., et al. (2022). RNA stable isotope probing and high-throughput sequencing to identify active microbial community members in a methane-driven denitrifying biofilm. J. Appl. Microbiol. 132, 1526–1542. doi: 10.1111/jam.15264
Yasir, M., Qureshi, A. K., Khan, I., Bibi, F., Rehan, M., Khan, S. B., et al. (2019). Culturomics-based taxonomic diversity of bacterial communities in the hot springs of Saudi Arabia. OMICS J. Integr. Biol. 23, 17–27. doi: 10.1089/omi.2018.0176
Zahn, J. A., Bergmann, D. J., Boyd, J. M., Kunz, R. C., and DiSpirito, A. A. (2001). Membrane-associated quinoprotein formaldehyde dehydrogenase from Methylococcus capsulatus Bath. J. Bacteriol. 183, 6832–6840. doi: 10.1128/JB.183.23.6832-6840.2001
Zelenkina, T. S., Eshinimayev, B. T., Dagurova, O. P., Suzina, N. E., Namsarayev, B. B., and Trotsenko, Y. A. (2009). Aerobic methanotrophs from the coastal thermal springs of Lake Baikal. Microbiology 78, 492–497. doi: 10.1134/S0026261709040134
Zeng, L., Tian, J., Chen, H., Wu, N., Yan, Z., Du, L., et al. (2019). Changes in methane oxidation ability and methanotrophic community composition across different climatic zones. J. Soils Sediments 19, 533–543. doi: 10.1007/s11368-018-2069-1
Zhou, Z., St. John, E., Anantharaman, K., and Reysenbach, A. L. (2022). Global patterns of diversity and metabolism of microbial communities in deep-sea hydrothermal vent deposits. Microbiome 10:241. doi: 10.1186/s40168-022-01424-7
Keywords: Aotearoa-New Zealand, geothermal, methane, methanotroph, oxidation, thermophile, transcriptomics
Citation: Houghton KM, Carere CR, Stott MB and McDonald IR (2023) Thermophilic methane oxidation is widespread in Aotearoa-New Zealand geothermal fields. Front. Microbiol. 14:1253773. doi: 10.3389/fmicb.2023.1253773
Edited by:
Mark Alexander Lever, The University of Texas at Austin, United StatesReviewed by:
Marina G. Kalyuzhanaya, San Diego State University, United StatesRobert Alexander Schmitz, ETH Zürich, Switzerland
Copyright © 2023 Houghton, Carere, Stott and McDonald. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Karen M. Houghton, ay5ob3VnaHRvbkBnbnMuY3JpLm56
†ORCID: Karen M. Houghton https://orcid.org/0000-0002-5618-6321
Carlo R. Carere https://orcid.org/0000-0001-6275-1203
Matthew B. Stott https://orcid.org/0000-0002-2082-9460