Next Article in Journal
Decoding the Genomic Variability among Members of the Bifidobacterium dentium Species
Next Article in Special Issue
Methane-Oxidizing Communities in Lichen-Dominated Forested Tundra Are Composed Exclusively of High-Affinity USCα Methanotrophs
Previous Article in Journal
New Insights in Laboratory Testing for COVID-19 Patients: Looking for the Role and Predictive Value of Human epididymis secretory protein 4 (HE4) and the Innate Immunity of the Oral Cavity and Respiratory Tract
Previous Article in Special Issue
Metagenomic Insight into Environmentally Challenged Methane-Fed Microbial Communities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic and Physiological Properties of a Facultative Methane-Oxidizing Bacterial Strain of Methylocystis sp. from a Wetland

1
Geologic Environment Research Division, Korea Institute of Geoscience and Mineral Resources, Daejeon 34132, Korea
2
Department of Microbiology, Chungbuk National University, Cheongju 28644, Korea
3
Department of Environmental Engineering, Chungnam National University, Daejeon 34134, Korea
*
Author to whom correspondence should be addressed.
Microorganisms 2020, 8(11), 1719; https://doi.org/10.3390/microorganisms8111719
Submission received: 2 September 2020 / Revised: 23 October 2020 / Accepted: 30 October 2020 / Published: 2 November 2020
(This article belongs to the Special Issue Biology, Diversity, and Ecology of Methanotrophic Bacteria)

Abstract

:
Methane-oxidizing bacteria are crucial players in controlling methane emissions. This study aimed to isolate and characterize a novel wetland methanotroph to reveal its role in the wetland environment based on genomic information. Based on phylogenomic analysis, the isolated strain, designated as B8, is a novel species in the genus Methylocystis. Strain B8 grew in a temperature range of 15 °C to 37 °C (optimum 30–35 °C) and a pH range of 6.5 to 10 (optimum 8.5–9). Methane, methanol, and acetate were used as carbon sources. Hydrogen was produced under oxygen-limited conditions. The assembled genome comprised of 3.39 Mbp and 59.9 mol% G + C content. The genome contained two types of particulate methane monooxygenases (pMMO) for low-affinity methane oxidation (pMMO1) and high-affinity methane oxidation (pMMO2). It was revealed that strain B8 might survive atmospheric methane concentration. Furthermore, the genome had various genes for hydrogenase, nitrogen fixation, polyhydroxybutyrate synthesis, and heavy metal resistance. This metabolic versatility of strain B8 might enable its survival in wetland environments.

1. Introduction

Wetlands, including bogs, marshes, and swamps, are typical areas where the methane cycle occurs actively and play a role in the carbon cycle [1,2]. Owing to methane production by methanogens, eutrophic wetlands are the main channel of methane fluxes [3,4]. Methane diffuses throughout wetlands [5,6] and creates a suitable environment for methanotrophs. Methanotrophic metabolism contributes to carbon cycling in wetland ecosystems [7,8]. The emission of methane is a major contributor to the greenhouse effect [9]. The importance of methanotrophic microorganisms in wetland carbon cycling and methane emissions has spurred investigations to identify and characterize aerobic methanotrophs.
Type I (including a subdivision of type X) methanotrophs belong to the class Gammaproteobacteria and Type II methanotrophs belong to the class Alphaproteobacteria [10,11]. Methanotrophs belonging to Verrucomicrobia were discovered in acidic thermal environments. Methane oxidation by methane monooxygenase is the first and most critical step in methane metabolism [12,13,14]. Methane is oxidized by particulate methane monooxygenase (pMMO) and soluble methane monooxygenase (sMMO) [15]. Whereas pMMO is bound to the cell membrane, sMMO exists in the cytoplasm [15]. Both pMMO and sMMO require distinct metal elements as a prosthetic group (copper ions and iron ions, respectively) [16,17,18]. Methanol is oxidized to formaldehyde by methanol dehydrogenase [19]. Formaldehyde is oxidized by several enzymes, including the formaldehyde-activating enzyme, through the methylene tetrahydromethanopterin (H4MPT) pathway [20,21]. Formate is oxidized by formate dehydrogenase to supply electrons to NAD [22]. Formaldehyde and formate are used as carbon sources depending on the methanotroph. In the RuMP pathway of type I methanotrophs, formaldehyde is used for cellular carbon production. Formaldehyde and ribulose 5-phosphate are converted to hexulose 6-phosphate through the pentose phosphate pathway, which is used to synthesize cellular carbohydrates. In the serine pathway of type II methanotrophs, formate is converted to methylenetetrahydrofolate (THF) and used for carbon metabolism [22].
Methanotrophs are not only able to reduce methane but can also produce value-added substances [23]. Recently, methanotrophs have become a promising option for mitigating methane emission [24]. Methylocystis parvus OBBP was known as one of the candidates for the cell factory production of Polyhydroxybutyrate (PHB) [25]. Other studies have been conducted on the synthesis of PHB by methane-oxidizing bacteria [26]. The production of organic acids (formate, acetate, lactate, and succinate) and hydrogen in Methylomicrobium alcaliphilum 20Z was considered a potential biocatalysis based on methane [27].
The genus Methylocystis includes type II methanotrophs belonging to the class Alphaproteobacteria and family Methylocystaceae [28]. The genus Methylocystis was first identified in 1970 [29] and was later validated by Bowman et al. [28]. The genus Methylocystis includes six species: M. rosea [30], M. hirsuta [31], M. bryophila [32], M. echinoides, M. parvus [28], and M. heyeri [33]. These species have been isolated from various environments, including groundwater aquifers, acidic peat-bog lakes, and wetland soils. Several Methylocystis genomes have been reported [34,35,36,37,38].
Methanotrophs were initially known as microorganisms which used only methane and methanol as carbon and energy sources. However, starting with Methylocella palustris, among the verified strains, it has been reported that facultative methanotrophs could utilize compounds with carbon-carbon bonds. The genus Methylocystis included some facultative methanotrophic species. The common source of multi-carbon was acetate. A study assumed that acetate utilization might be important for survival in environments where methane availability is variable or limited [39]. Such cases included wetlands and upland soils. Methylocystis spp. were predominant in these environments [40,41].
In this study, a novel methanotroph was isolated from a wetland and characterized to reveal its role in methane cycling in this environment. Genome properties were revealed using comparative genomics. The physiological and genomic data provide insight into the adaptations of Methylocystis methanotrophs in wetland environments.

2. Materials and Methods

2.1. Enrichment Culture and Isolation

A water sample was obtained from an artificial wetland in Ok-Cheon, Republic of Korea (36.19 N, 127.33 E) in June 2017. Temperature, electrical conductivity (EC), oxidation-reduction potential (ORP), pH, and dissolved oxygen (DO) of sampled water were checked by a Multimeter (HQ40d, Hach, Loveland, CO, USA) onsite. To examine the water composition, the water was filtered with a syringe filter (0.45 μm pore size) and analyzed with ion chromatography (IC, Metrohm, Herisau, Switzerland) and inductively coupled plasma-optical emission spectrometry (ICP-OES, Horiba Jobin Yvon, Longjumeau, France). The 50 mL sample was filtered through a 0.1 μm pore polycarbonate (PC) filter membrane to collect cells. The collected cells were incubated in 10 mL of modified nitrate mineral salt (NMS) medium in a Balch-type tube (CLS-4209, Chemglass, Vineland, NJ, USA) at 25 °C. Methane (filter sterilized, 0.22 μm pore size) was added to the headspace air (15%, v/v) through a syringe. One liter of the medium contained 0.0985 g MgSO4·7H2O, 0.0147 g CaCl2·2H2O, 0.0348 g K2SO4, 0.017 g NaNO3, 0.0136 g KH2PO4, 0.252021 g NaHCO3, 1 mL SL-10 trace element solution [42], 1 mL vitamin solution [43], and 0.028 g cerium [44]. The culture was serially diluted from 10−1 to 10−10 several times. Furthermore, the dilution-to-extinction cultivation method was conducted in 96-well plates. The emerging turbidity among the wells was used for isolation using the floating membrane method with the polycarbonate membrane (0.1 μm pore size) on the NMS medium at various methanol concentrations (0.1%, 0.5%, and 1%) at 25 °C. After 2 weeks, small colonies were observed on the floating membrane. The isolated cell was named B8. Cells were routinely cultured in a modified NMS medium. The purity of culture was checked by monitoring the T-RFLP of 16S rRNA gene PCR product [45] from cultures grown with either methane or acetate after incubation. In addition, no growth was confirmed in the modified NMS medium with 0.1% yeast extract as the sole energy source. Furthermore, contamination by other heterotrophic bacteria was routinely checked on nutrient agar and R2A agar plates.

2.2. Physiologic Analysis

Gram staining was performed using a commercial kit (BD, Franklin Lakes, NJ, USA). Cell morphology and cell size were observed by differential interference contrast microscopy (Axio Scope.A1; Carl Zeiss, Jena, Germany). Catalase and oxidase tests were performed with 3% hydrogen peroxide (H2O2) and oxidase strips (Sigma-Aldrich. St. Louis, MO, USA), respectively. Strain B8 was grown at 4 °C, 10 °C, 15 °C, 20 °C, 25 °C, 30 °C, 35 °C, 37 °C, 40 °C, and 45 °C. NaCl tolerance was assessed by growth in the presence of 0% to 5% NaCl (in 0.5% increments). The pH was tested from 4 to 12 in increments of 0.5. To identify nitrogen sources used by strain B8, the sample was incubated in a modified NMS medium with various nitrogen sources (cysteine, yeast extract, casamino acid, betaine, carsitone, and ammonium sulfate) by replacing 0.05% (w/v) sodium nitrate. Carbon source utilization was assessed using 0.05% ethanol, glucose, acetate, or pyruvate. Tests were performed in duplicates. Colony formation was checked on the solid medium containing 1.5% agarose, 1% phytagel, and 1% gellan gum. PHB formation was estimated using Nile Red staining followed by examination using an Axio Scope A1 fluorescence microscope (Carl Zeiss, Oberkochen, Germany). To calculate the growth rate, strain B8 was cultured in a modified NMS medium with added methane (15%, v/v) at 30 °C and 150 rpm. The growth rate was monitored in triplicates. The methane concentration in headspace and cell growth was measured for 7 days. Cell growth (optical density at 600 nm, OD600) was determined using spectrophotometers (Optizen POP; Mecasys, Daejeon, Korea and Spectra Max 190 microplate reader; Molecular Devices, CA, USA). Methane was measured by gas chromatography (GC). The detailed condition was described in Section 2.3. Also, the quantity of 16S rRNA and pmoA gene was measured by Q-PCR (CFX96 Real-time System, Bio-rad, Hercules, CA, USA) with primers Bac518F/Bac786R [46] and A189F/A682R [47], respectively.

2.3. Hydrogen Utilization Test

Hydrogen utilization was evaluated in a modified NMS medium with a 2% H2 amended or H2 nonamended headspace and various methane and oxygen concentration based on Hakobyan et al. [48]. The conditions were as follows. Set1: 20% methane and 15% oxygen; Set2: 20% methane and 3% oxygen; Set3: 6% methane and 15% oxygen; set4: 6% methane and 3% oxygen; and set5: 20% methane and 9% oxygen. Before gas injection, headspace was replaced with Ar gas. Then, all gases were injected with a syringe through a syringe filter (0.22 μm pore size) for sterilization. The hydrogen and methane concentration and OD600 were measured during the incubation time.
Methane and hydrogen concentrations in the headspace were monitored by a GC-flame ionization detector/thermal conductivity detector (GC-FID/TCD, Agilent 6890N, Santa Clara, CA, USA) using a 12390-U column. The temperature of the injector was at 200 °C isothermally. The oven temperature was maintained at 65 °C for 1 min, then ramped at a rate of 25 °C/min to 225 °C and held for 2 min. The temperatures of FID and TCD were 250 °C and 230 °C, respectively. The methane and hydrogen were determined with 200 μL of gaseous sample of the headspace of the bottle.

2.4. Phylogenetic Analysis

The genome was extracted using a DNA extraction kit (QIAGEN, Valencia, CA, USA). The 16S ribosomal RNA gene was amplified using 27F/1492R primers [49,50], and the methane monooxygenase gene was amplified using A189F/A682R [47], mmox1/mmox2 [51], and mmox206F/mmox886R [52] primers. PCR conditions were as previously described [47,51,52]. Phylogenetic trees based on 16S rRNA genes and pmoA genes aligned using Bioedit were constructed using the neighbor-joining method [53] in MEGA 7 [54] based on the Kimura 2-parameter model and Poisson model using 1000 replicates, respectively. The 16S rRNA gene similarity was calculated using EzBioCloud [55].

2.5. Genomic Analysis

A paired-end DNA library was constructed using a Truseq Nano DNA kit (Illumina, San Diego, CA, USA). Sequencing was performed by an Illumina HiseqX device at Macrogen (Seoul, Korea). The quality of the raw reads was checked by FastQC. SPAdes (ver. 3.13) was used to assemble the raw reads [56]. The assembled methanotroph scaffolds were annotated using the NCBI Prokaryotic Genome Annotation Pipeline [57]. Gene function was analyzed using BLASTp [58] against the nonredundant (NR) database and the Kyoto Encyclopedia of Genes and Genomes [59]. A domain search was performed using PfamScan [60]. Clusters of orthologous groups (COG) classification was performed according to bac-genomics-scripts (https://github.com/aleimba/bac-genomics-scripts). The average nucleotide identities (ANI) of B8 and other Methylocystis species were analyzed using orthologous average nucleotide identity tool (OAT) software [61]. Average amino acid identity (AAI) was calculated using CompareM (https://github.com/dparks1134/CompareM). Pan-genome analysis was performed using PGAP (ver. 1.2.1) [62]. The phylogenomic tree was reconstructed using Anvio-6.1 (default parameter and hmm source: Bacteria_71) [63]. The antibiotic resistance genes were analyzed using a comprehensive antibiotic resistance database (CARD) [64].

2.6. Data Availability

The 16S rRNA gene and draft genome were deposited in the NCBI GenBank as accession numbers MN527243 and VBTZ00000000, respectively.

3. Results and Discussion

3.1. Information on Sampling Site

The artificial wetland in Ok-Cheon, Republic of Korea (36.19 N, 127.33 E) was constructed to remove nonpoint source pollution from the Daecheong lake, a source of drinking and irrigation water [65]. The water composition of the artificial wetland at the sampling time was summarized in Table 1. The temperature was over 15 °C, except for the winter season [66].

3.2. Physiological Characterization

Cells of strain B8 were 0.7 μm in width and 1.5 μm in length. Cells were gram-negative, oxidase- and catalase-positive, coccobacilli or short rod-shaped, and appeared as single cells or in aggregates, but not as rosettes. PHB formation was observed by fluorescence microscopy of Nile Red stained sections (Supplementary Figure S1). Colonies did not form on NMS agar, gellan, or phytagel plates. The harvested cell pellet was creamy white. Strain B8 grew between 15 °C and 37 °C, with optimal growth at 30–35 °C. The bacteria were not tolerant to high salinity and grew between 0% to 0.5% NaCl concentration. Growth occurred at pH values ranging from pH 6.5 to 10. The optimal pH was 8.5 to 9, indicating that strain B8 is slightly alkaliphilic. Nitrate ammonia, yeast extract, and cysteine were used as nitrogen sources, while carsitone, casamino acid, and betaine were not. Growth rate was 0.11 ± 0.01 h−1 in the modified NMS medium with methane (20%, v/v) at 30 °C and 150 rpm. Strain B8 could use methanol at levels ranging from 0.05% to 1%. The optimum level was 0.1%. Ethanol, glucose, and pyruvate were not used as carbon sources. Methane and methanol were the carbon and energy sources for strain B8. In addition, growth was observed using acetate as a carbon and energy source after 1 month of incubation. This result was also reported from other facultative methanotrophs, such as M. bryophila H2s, M. heyeri H2, M. echinoides IMET10491 [39], and Methylocystis sp. SB2 [67]. It might serve as a survival strategy in environments with limited methane concentration, such as peat [39]. Table 2 summarizes the findings and compares them with close relatives identified using the phylogenetic analysis described. The physiological properties of strain B8 indicate that this strain is slightly alkaliphilic compared to other members of Methylocystis. This might be advantageous in the development of alkaline conditions during photosynthesis in wetlands [68]. Strain B8 was deposited in the Korean Collection for Type Cultures (KCTC) as KCTC82145.

3.3. Hydrogen Utilization and Production

The hydrogen utilization test was conducted using various methane and oxygen conditions. Recently, it was reported that Methylocystis sp. SC2 uses hydrogen as an energy source [48]. The growth yield increased significantly under 6% CH4 and 3% O2 conditions when 2% H2 was added to headspace culture. However, the effect was not observed using strain B8. Instead, hydrogen was produced when oxygen decreased to an oxygen limitation level without gas exchange (Supplementary Figure S2). A similar result was found for Methylomonas sp. DH-1 [71] and Methylomicrobium alcaliphilum 20Z [27], which are type I methanotrophs. GC results showed that this result was related to the type of hydrogenase in Methylocystis sp. B8. For additional information, see Section 3.9.

3.4. Genomic Features

The genomic information of strain B8 is shown in Table 3. The G + C content was 59.9% and the lowest in the genomes of the Methylocystis species. Genes of the genome were classified according to function using COG analysis (Supplementary Table S1). Genes with ‘no functional prediction’ (S) were most prevalent (13%), followed by those with ‘general function prediction only’ (R, 11%) and related to ‘energy production and conversion’ (C, 8%). Contigs 10 (116,387 bp) and 12 (92,240 bp) displayed a cluster of plasmid replication genes (repA, repB, and repC), respectively. Genes in the clusters FEV16_14560-FEV16_14570 and FEV16_15250-FEV16_15260 of strain B8 were homologous with P1_12-P1_14 and P2_74-76 of Methylocystis sp. SC2. These findings indicate that the two contigs belong to a plasmid, further indicating that strain B8 might harbor two large plasmids. Contig 10 also contained mercury reductase genes similar to those of strain SC2 [72].

3.5. Phylogeny

Phylogenetic analysis of the 16S rRNA gene sequence showed that strain B8 clustered with the genus Methylocystis in Alphaproteobacteria. Strain B8 was closely related to M. hirsuta CSC1T, M. rosea SV97T, and M. echinoides IMET10491T with 98.8%, 98.4%, and 97.9% 16S rRNA gene similarity, respectively (Supplementary Figure S3). The genome closest to that of strain B8 was Methylocystis sp. SC2 (HE956757), with the highest ANI value of 84.76%. The lowest ANI value (73.06) was observed in the genome of M. bryophila. The highest and lowest AAI values were obtained with Methylocystis rosea BRCS1 (86.9%) and M. bryophila S285 (68.5%), respectively (Supplementary Table S2). The ANI and AAI results, combined with the phylogenetic analysis results based on 16S rRNA gene sequences, indicate that strain B8 is a novel species in the genus Methylocystis based on the threshold values suggested by Konstantidis and Tiedje [73] and Konstantidis et al. [74]. The AAI result showed that genomes of 16 species in Methylocystis have been reported or deposited (in NCBI) so far. Furthermore, the phylogenomic tree was reconstructed with genomes of Methylocystis, Methylosinus, and other genera in the Methylocystaceae family (Figure 1). It supported that strain B8 is a new species in the genus Methylocystis. Furthermore, we suggest that Methylocystis sp. LW5 should be reclassified as a member of the genus Methylosinus based on the phylogenomic tree.

3.6. Pangenome Analysis

Eleven strains of Methylocystis (Table 3) were compared based on the pan-genome analysis. The Methylocystis pan-genome consisted of 13,554 clusters with 1537 core gene clusters in our study. A total of 328 clusters (335 genes) were identified as unique genes in strain B8 (Supplementary Table S3). The COG classification of the unique genes showed that major COGs were assigned to R (general functional prediction), S (function unknown), and M (cell wall, membrane, and envelope biogenesis). As COG M predominantly includes genes for lipopolysaccharide production, these genes might be related to niche adaptation by modifying the surface structures of cells [75]. Recently, another study intensively analyzed the pangenome analysis in the genus Methylocystis/Methylosinus [76]. Briefly, the genera Methylocystis and Methylosinus had encoding genes related to methane oxidation and nitrogen fixation. However, their types were different between species. One of the major differences between Methylocystis and Methylosinus was motility. Furthermore, some species in the genus Methylocystis had metabolic flexibility, such as growth on acetate and ethanol as carbon and energy sources, as well as having photosynthetic-related genes.

3.7. Methanotrophic Pathway

Two gene clusters coding isozymes of pMMO were found in the genome of strain B8. Genes coding for sMMO were not found. The pMMO types were classified using a phylogenic tree (Figure 2). The two pMMO operons, pmoCAB1 and pmoCAB2, coded low-affinity methane oxidation (LAMO) and high-affinity methane oxidation (HAMO), respectively. Additionally, three singleton pmoC paralogs were found in the genome of strain B8. Genes related to C1 compound metabolism are shown in Supplementary Table S4. Previous studies have reported that additional gene copies of pmoC are essential for methanotrophic growth [72,77]. Furthermore, Matsen et al. [78] suggested that the role of PmoC paralogs might be related to the regulation of gene expression or methane sensing, rather than catalytic activity.
The expressions of the HAMO and LAMO genes are regulated by environmental methane concentration [80]. Methanotrophs harboring both HAMO and LAMO enzymes can modify their behavior depending on methane concentration and adapt to various situations. Both types of pMMO were detected in Methylocystis sp. SC2, M. heyeri, M. parvus, and M. bryophila S285 (Supplementary Table S5). M. bryophila, M. heyeri, and M. hirsuta also had sMMO. The presence of HAMO and LAMO in the genome of strain B8 might indicate the potential of the strain to survive various methane conditions. In addition, several studies have reported that methanotrophs (especially the Methylocystis/Methylosinus group) with HAMO activity contribute to the oxidation of atmospheric methane concentrations in soil [79,81]. The utilization of atmospheric methane by strain B8 needs to be demonstrated in a future study.
Methanol dehydrogenase types can be classified depending on whether calcium (Mxa) or lanthanides (Xox) are required for activity [82,83,84,85]. Both Mxa, and Xox were detected in strain B8 (Supplementary Table S5). All the evaluated strains, except for KS32, had both Mxa and Xox (Supplementary Table S5). XoxF reportedly oxidizes methanol to formate rather than formaldehyde [85], although disparate results were recently published [86]. One study reported that the gene expression of Mxa and Xox in Methylobacter tundripaludum may change depending on the partner microbes [87]. The function of XoxF in Methylocystis requires further study. The XoxF and MxaFI methanol dehydrogenases require pyrroloquinoline quinone (PQQ) as a catalytic cofactor [88,89]. The genome of strain B8 contained the genes pqqBCE and pqqA, which are required for PQQ biosynthesis. The following steps are described in the Supplementary Information.

3.8. Nitrogen Metabolism

Most methanotrophs possess genes related to nitrogen fixation [90,91,92]. Molybdenum-iron (Mo-Fe) nitrogenase was observed in all the genomes of Methylocystis. Vanadium-iron (V-Fe) nitrogenases were only found in other methanotrophic bacteria from wetland ecosystems (i.e., M. heyeri H2T, M. bryophila S285, and M. parvus BRSC2) [76]. Although strain B8 was obtained from wetlands, it only had genes related to molybdenum-iron (Mo-Fe) nitrogenase. These organisms can assimilate various nitrogen sources, including nitrate, nitrite, and ammonia [34,44,93]. In this study, genes related to nitrogen utilization were found in the genome of strain B8 (Supplementary Table S6). These included assimilatory nitrate/nitrite reduction-related genes, nitrate reductase (FEV16_06420) and nitrite reductase (FEV16_12950) genes, and genes for ammonium transport. These findings suggest that strain B8 could utilize N2, nitrate, nitrite, and ammonia as nitrogen sources. Strain B8 also harbored genes coding for hydroxylamine oxidoreductase (FEV16_03625 and FEV16_12730), which oxidizes hydroxylamine to nitrite, and hydroxylamine reductase (FEV16_04040), which catalyzes the reduction of hydroxylamine to ammonia, as observed in Methylocystis sp. SC2 [72]. These enzymes have roles in hydroxylamine detoxification as well as nitrogen acquisition.

3.9. Hydrogenase

Genes coding for hydrogenase and accessory proteins were observed in all Methylocystis genomes except Methylocystis sp. KS32. The phylogenetic tree, constructed using the hydrogenase gene sequences (Figure 3), revealed a variation in hydrogenase types within this genus. Methylocystis sp. SC2, which is known to use hydrogen as an energy source, had five types of hydrogenase: Ia, Id, 2b, 3b, and 1h/5 group hydrogenase. Among these genes, 1d and 2b group hydrogenases were related to hydrogen utilization by a proteomic study [48]. However, strain B8 had only two groups of hydrogenases, FEV16_01595 and FEV16_11060, which belonged to the 1h/5 and 3b groups, respectively. Nickel (Ni) and iron (Fe) ions are present in both active sites of hydrogenases [94,95,96]. Hydrogen utilization was not observed in strain B8 in our experiments, but hydrogen was produced. The production is hydrogen may have been related to the 3b group hydrogenase. The 3b group hydrogenases reversibly reduce NAD+ coupled to H2 oxidation [97]. The membrane-bound 1h/5 group hydrogenases are used for aerobic respiration under starvation conditions [97]. It has been hypothesized that hydrogen might be used under starvation conditions. In this study, we demonstrated the hydrogen production potential of B8, a type II methanotroph. The utilization and production of H2 in strain B8 need further research in combination with omics studies. The produced hydrogen by methanotrophs in oxygen-limited environments might be an energy source for non-methanotrophs [27].

3.10. Polyhydroxybutyrate (PHB) Synthesis

Some methanotrophs can synthesize biopolymers, such as PHB, to store carbon and as energy sources [98,99]. Strain B8 was expected to be able to produce and accumulate PHBs because of the presence of genes related to PHB synthesis and degradation. These genes included two poly β-hydroxybutyrate polymerases (FEV16_00425 and 06360), two poly β-hydroxybutyrate depolymerases (FEV16_05615 and 5645), acetyl-CoA acetyltransferase (FEV16_02735), acetoacetyl-CoA reductases (FEV16_02730), polyhydroxyalkanoate synthesis repressors (FEV16_02740), and three phasins (FEV16_07845, 13905, and 15530) (Supplementary Table S7). Storing PHB could lead to a reducing power when cells are starved [99,100]. Biodegradable and biocompatible artificial biopolymers have been developed to replace plastics derived from petroleum [99,101]. The production of PHB-based bioplastics from methane could be advantageous in terms of preventing environmental pollution.

3.11. Resistance Genes

Genes for glutaredoxin-dependent arsenate reductase (FEV16_04625, FEV16_04855, and FEV16_04965) and the arsenite efflux pump (FEV16_04970) involved in neutralization or exporting arsenic were demonstrated in strain B8. Arsenic is an especially toxic heavy metal because of its characteristic resemblance to phosphorus [102,103]. Further, arsenic is a redox-sensitive element with a geogenic origin [104]. The detoxification mechanism of arsenic might be necessary for survival in wetlands that feature large redox fluctuations.
Genes for cyanide detoxification were detected in strain B8. Cyanide is lethal to cells because it interrupts the electron transport systems and halts cellular energy production. Many organisms have defense mechanisms to neutralize cyanide [105,106]. The genes in strain B8 include those coding 3-mercaptopyruvate sulfurtransferase (FEV16_02590 and FEV16_16540) and rhodanese (FEV16_09915, FEV16_10965, and FEV16_13300), which catalyze the conversion of cyanide to thiocyanate.
Mercury is a toxic heavy metal, and many organisms have a mechanism to neutralize mercury toxicity [107]. Genes related to mercury detoxification, including mercury reductase (FEV16_05105, FEV16_14585, FEV16_14590, FEV16_14595, FEV16_14600, FEV16_15400, and FEV16_15405), were detected in the genome of strain B8. Mercury (II) reductase catalyzes the reduction of mercury (Hg2+) to the less toxic elemental mercury (Hg0). In addition, methanobactin (Mbn) was included in this genome (See Supplementary Information and Supplementary Figure S4). Mbn is needed under copper-limited conditions and can also be used for mercury detoxification [108]. In the genus Methylocystis, genes for mercury reduction have also been demonstrated in M. rosea, Methylocystis sp. KS32, M. hirsuta, Methylocystis sp. ATCC 49,242, and Methylocystis sp. MitZ-2018 (Table 4). The potential for mercury reduction as a means of detoxification in Methylocystis sp. HL18 was suggested based on metagenomic analysis [109]. Methylococcus capsulatus Bath (type X methanotroph) showed mercury reduction activity [110]. Methylosinus trichosporium OB3b can demethylate methylmercury but does not reduce mercury [111]. Further studies are needed to confirm the mercury reduction activity (including the demethylation of methylmercury) in Methylocystis.
The genome had genes related to tetracycline and quinolone antibiotic resistance gene (FEV16_04740) based on CARD.

4. Conclusions

The physiological and metabolic properties and genomic information of a facultative methanotrophic bacterial strain isolated from a wetland were determined. Phylogenetic/phylogenomic analysis and ANI and AAI values indicated that strain B8 is a novel species in the genus Methylocystis. The physiological and genomic analysis demonstrated the versatility of strain B8 in the adaptation to diverse environments. The presence of genes coding HAMO and LAMO indicated that B8 could survive under fluctuating methane concentrations. Strain B8 also harbors genes coding for hydrogenase to produce hydrogen under oxygen-limited conditions. The presence of genes for heavy metal resistance and PHB formation might be useful for survival under stress conditions. These collective features suggest the ubiquitous distribution of strain B8-like methanotrophs and their active roles in methane cycles in wetlands.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/8/11/1719/s1. Figure S1: The cells observed with (a) DIC microscopy and (b) fluorescence microscopy stained by Nile red. Figure S2: Monitoring of methane and hydrogen concentrations in headspace of bottles. Methane concentrations in (a) hydrogen non-amended sets and (b) hydrogen-amended sets, hydrogen concentrations in (c) hydrogen non-amended sets and (b) hydrogen-amended sets, and optical density(OD600) of (e) non-amended hydrogen sets and (f) hydrogen-amended sets with 200 μL sample using Spectra Max 190 microplate reader. Figure S3: The phylogenetic tree of Methylocystis sp. B8 and other strains in the family Methylocystaceae based on the 16S ribosomal RNA gene. The phylogenetic tree was reconstructed using the neighbor-joining method. The outgroup were Rhizobium leguminosarum and Sphingomonas aestuarii in the family Rhizobiaceae and Sphingomonadaceae respectively. Bootstrap values below 70% are not shown. Bar, 0.01 substitutions per nucleotide site. Figure S4: Methanobactin gene cluster of strain B8 and other Methylocystis species. Table S1. COG classification of the genome of strain B8. Table S2: AAI results of the genera Methylocystis and Methylosinus (excel). Table S3: COG classification based on the pan-genome and core-genome of the genus Methylocystis and the unique clusters of strain B8. Table S4: Genes related to C1 compound metabolism in strain B8. Table S5: Comparative analysis of genes related to methane and methanol utilization in Methylocystis. Table S6: Genes related to nitrogen metabolism in strain B8. Table S7: Genes related to PHB formation in strain B8.

Author Contributions

G.-Y.J.: investigation (lead); visualization (equal); writing—original draft preparation (lead); S.-K.R.: supervision (lead); writing-review and editing (equal); Y.-S.H.: investigation (equal); S.-J.K.: conceptualization (equal); data curation (lead); visualization (equal); writing-review and editing (equal). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Basic Research Project of the Korea Institute of Geoscience and Mineral Resources (20-3412-1), Republic of Korea and the National Research Council of Science & Technology (NST) grant by the Korean government (MSIP) (Project code: CAP-19-05-KIGAM).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Whiting, G.J.; Chanton, J.P. Primary production control of methane emission from wetlands. Nat. Cell Biol. 1993, 364, 794–795. [Google Scholar] [CrossRef]
  2. Mitsch, W.J.; Bernal, B.; Nahlik, A.M.; Mander, Ü.; Zhang, L.; Anderson, C.J.; Jorgensen, S.E.; Brix, H. Wetlands, carbon, and climate change. Landsc. Ecol. 2012, 28, 583–597. [Google Scholar] [CrossRef]
  3. Shao, X.; Sheng, X.; Wu, M.; Wu, H.; Ning, X. Methane production potential and emission at different water levels in the restored reed wetland of Hangzhou Bay. PLoS ONE 2017, 12, e0185709. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Kirschke, S.; Bousquet, P.; Ciais, P.; Saunois, M.; Canadell, J.G.; Dlugokencky, E.; Bergamaschi, P.; Bergmann, D.; Blake, D.R.; Bruhwiler, L.M.P.; et al. Three decades of global methane sources and sinks. Nat. Geosci. 2013, 6, 813–823. [Google Scholar] [CrossRef]
  5. Bastviken, D.; Cole, J.; Pace, M.; Tranvik, L. Methane emissions from lakes: Dependence of lake characteristics, two regional assessments, and a global estimate. Glob. Biogeochem. Cycles 2004, 18, 18. [Google Scholar] [CrossRef]
  6. Conrad, R. The global methane cycle: Recent advances in understanding the microbial processes involved. Environ. Microbiol. Rep. 2009, 1, 285–292. [Google Scholar] [CrossRef] [PubMed]
  7. King, G.M. Ecological aspects of methane oxidation, a key determinant of global methane dynamics. In Advances in Microbial Ecology; Springer: Berlin, Germany, 1992; pp. 431–468. [Google Scholar]
  8. Chowdhury, T.R.; Dick, R.P. Ecology of aerobic methanotrophs in controlling methane fluxes from wetlands. Appl. Soil Ecol. 2013, 65, 8–22. [Google Scholar] [CrossRef]
  9. Stocker, T.F.; Qin, D.; Plattner, G.-K.; Tignor, M.; Allen, S.K.; Boschung, J.; Nauels, A.; Xia, Y.; Bex, V.; Midgley, P.M. Climate change 2013: The physical science basis. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change. In Evaluation of Climate Models; Cambridge University Press: Cambridge, UK, 2013; p. 1535. [Google Scholar]
  10. Dedysh, S.N.; Eknief, C. Diversity and phylogeny of described aerobic methanotrophs. In Methane Biocatalysis: Paving the Way to Sustainability; Springer: Berlin, Germany, 2018; pp. 17–42. [Google Scholar]
  11. Semrau, J.D.; DiSpirito, A.A.; Yoon, S. Methanotrophs and copper. FEMS Microbiol. Rev. 2010, 34, 496–531. [Google Scholar] [CrossRef] [Green Version]
  12. Hanson, R.S. The obligate methanotrophic bacteria Merhvlococcus. Methylomonas, and Methylosinus. In The Procaryotes; Springer: Berlin, Germany, 1992; pp. 2350–2364. [Google Scholar]
  13. Oremland, R.S.; Culbertson, C.W. Importance of methane-oxidizing bacteria in the methane budget as revealed by the use of a specific inhibitor. Nat. Cell Biol. 1992, 356, 421–423. [Google Scholar] [CrossRef]
  14. Dalton, H. Methane Oxidation by Methanotrophs. In Methane and Methanol Utilizers; Springer: Berlin, Germany, 1992; pp. 85–114. [Google Scholar]
  15. Murrell, J.C.; Gilbert, B.; McDonald, I.R. Molecular biology and regulation of methane monooxygenase. Arch. Microbiol. 2000, 173, 325–332. [Google Scholar] [CrossRef]
  16. Fox, B.G.; Afroland, W.; Jollie, D.R.; Lipscomb, J.D. Methane monooxygenase from Methylosinus trichosporium OB3b. Methods Enzymol. 1990, 188, 191–202. [Google Scholar] [CrossRef] [PubMed]
  17. DeWitt, J.G.; Bentsen, J.G.; Rosenzweig, A.C.; Hedman, B.; Green, J.; Pilkington, S.; Papaefthymiou, G.C.; Dalton, H.; Hodgson, K.O.; Lippard, S.J. X-ray absorption, Moessbauer, and EPR studies of the dinuclear iron center in the hydroxylase component of methane monooxygenase. J. Am. Chem. Soc. 1991, 113, 9219–9235. [Google Scholar] [CrossRef]
  18. Nguyen, H.H.; Shiemke, A.K.; Jacobs, S.J.; Hales, B.J.; Lidstrom, M.E.; Chan, S.I. The nature of the copper ions in the membranes containing the particulate methane monooxygenase from Methylococcus capsulatus (Bath). J. Biol. Chem. 1994, 269, 14995–15005. [Google Scholar]
  19. Anthony, C.; Zatman, L.J. The microbial oxidation of methanol. The prosthetic group of the alcohol dehydrogenase of Pseudomonas sp. M27: A new oxidoreductase prosthetic group. Biochem. J. 1967, 104, 960–969. [Google Scholar] [CrossRef] [Green Version]
  20. Vorholt, J.A.; Chistoserdova, L.; Lidstrom, M.E.; Thauer, R.K. The NADP-dependent methylene tetrahydromethanopterin dehydrogenase in Methylobacterium extorquens AM1. J. Bacteriol. 1998, 180, 5351–5356. [Google Scholar] [CrossRef] [Green Version]
  21. Marx, C.J.; Chistoserdova, L.; Lidstrom, M.E. Formaldehyde-detoxifying role of the tetrahydromethanopterin-linked pathway in Methylobacterium extorquens AM. J. Bacteriol. 2003, 185, 7160–7168. [Google Scholar] [CrossRef] [Green Version]
  22. Chistoserdova, L.; VanWiggeren, G.D.; Roy, R. C1 transfer enzymes and coenzymes linking methylotrophic bacteria and methanogenic archaea. Science 1998, 281, 99–102. [Google Scholar] [CrossRef]
  23. Strong, P.J.; Xie, S.; Clarke, W.P. Methane as a resource: Can the methanotrophs add value? Environ. Sci. Technol. 2015, 49, 4001–4018. [Google Scholar] [CrossRef]
  24. Haynes, C.A.; Gonzalez, R. Rethinking biological activation of methane and conversion to liquid fuels. Nat. Chem. Biol. 2014, 10, 331–339. [Google Scholar] [CrossRef]
  25. Bordel, S.; Rojas, A.; Muñoz, R. Reconstruction of a genome scale metabolic model of the polyhydroxybutyrate producing methanotroph Methylocystis parvus OBBP. Microb. Cell Factories 2019, 18, 1–11. [Google Scholar] [CrossRef] [PubMed]
  26. Strong, P.J.; Laycock, B.; Mahamud, S.N.S.; Jensen, P.D.; Lant, P.A.; Tyson, G.W.; Pratt, S. The opportunity for high-performance biomaterials from methane. Microorganisms 2016, 4, 11. [Google Scholar] [CrossRef]
  27. Kalyuzhnaya, M.G.; Yang, S.; Rozova, O.N.; Smalley, N.E.; Clubb, J.; Lamb, A.; Gowda, G.A.N.; Raftery, D.; Fu, Y.; Bringel, F.; et al. Highly efficient methane biocatalysis revealed in a methanotrophic bacterium. Nat. Commun. 2013, 4, 1–7. [Google Scholar] [CrossRef]
  28. Bowman, J.P.; Sly, L.I.; Nichols, P.D.; Hayward, A.C. 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. 1993, 43, 735–753. [Google Scholar] [CrossRef]
  29. Whittenbury, R.; Davies, S.L.; Davey, J.F. Exospores and cysts formed by methane-utilizing bacteria. J. Gen. Microbiol. 1970, 61, 219–226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Wartiainen, I.; Hestnes, A.G.; McDonald, I.R.; Svenning, M.M. Methylocystis rosea sp. nov., a novel methanotrophic bacterium from Arctic wetland soil, Svalbard, Norway (78° N). Int. J. Syst. Evol. Microbiol. 2006, 56, 541–547. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Lindner, A.S.; Pacheco, A.; Aldrich, H.C.; Staniec, A.C.; Uz, I.; Hodson, D.J. Methylocystis hirsuta sp. nov., a novel methanotroph isolated from a groundwater aquifer. Int. J. Syst. Evol. Microbiol. 2007, 57, 1891–1900. [Google Scholar] [CrossRef] [PubMed]
  32. Belova, S.E.; Kulichevskaya, I.S.; Bodelier, P.L.E.; Dedysh, S.N. Methylocystis bryophila sp. nov., a facultatively methanotrophic bacterium from acidic Sphagnum peat, and emended description of the genus Methylocystis (ex Whittenbury et al. 1970) Bowman et al. Int. J. Syst. Evol. Microbiol. 2013, 63, 1096–1104. [Google Scholar] [CrossRef] [Green Version]
  33. Dedysh, S.N.; Belova, S.E.; Bodelier, P.L.E.; Smirnova, K.V.; Khmelenina, V.N.; Chidthaisong, A.; Trotsenko, Y.A.; Liesack, W.; Dunfield, P.F. Methylocystis heyeri sp. nov., a novel type II methanotrophic bacterium possessing ‘signature’ fatty acids of type I methanotrophs. Int. J. Syst. Evol. Microbiol. 2007, 57, 472–479. [Google Scholar] [CrossRef] [Green Version]
  34. Dam, B.; Dam, S.; Blom, J.; Liesack, W. Genome analysis coupled with physiological studies reveals a diverse nitrogen metabolism in Methylocystis sp. strain SC2. PLoS ONE 2013, 8, e74767. [Google Scholar] [CrossRef] [Green Version]
  35. Stein, L.Y.; Bringel, F.; DiSpirito, A.A.; Han, S.; Jetten, M.S.M.; Kalyuzhnaya, M.G.; Kits, K.D.; Klotz, M.G.; Camp, H.J.M.O.D.; Semrau, J.D.; et al. Genome sequence of the methanotrophic alphaproteobacterium Methylocystis sp. Strain Rockwell (ATCC 49242). J. Bacteriol. 2011, 193, 2668–2669. [Google Scholar] [CrossRef] [Green Version]
  36. Han, D.; Dedysh, S.N.; Liesack, W. Unusual genomic traits suggest Methylocystis bryophila S285 to be well adapted for life in peatlands. Genome Biol. Evol. 2018, 10, 623–628. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Nguyen, N.-L.; Yu, W.-J.; Gwak, J.-H.; Kim, S.-J.; Park, S.-J.; Herbold, C.W.; Kim, J.-G.; Jung, M.-Y.; Rhee, S.-K. Genomic insights into the acid adaptation of novel methanotrophs enriched from acidic forest soils. Front. Microbiol. 2018, 9, 1982. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Vorobev, A.; Jagadevan, S.; Jain, S.; Anantharaman, K.; Dick, G.J.; Vuilleumier, S.; Semrau, J.D. Genomic and Transcriptomic analyses of the facultative methanotroph Methylocystis sp. Strain SB2 grown on methane or ethanol. Appl. Environ. Microbiol. 2014, 80, 3044–3052. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Belova, S.E.; Baani, M.; Suzina, N.E.; Bodelier, P.L.; Liesack, W.; Dedysh, S.N. Acetate utilization as a survival strategy of peat-inhabiting Methylocystis spp. Environ. Microbiol. Rep. 2011, 3, 36–46. [Google Scholar] [CrossRef]
  40. Knief, C.; Dunfield, P.F. Response and adaptation of different methanotrophic bacteria to low methane mixing ratios. Environ. Microbiol. 2005, 7, 1307–1317. [Google Scholar] [CrossRef]
  41. Dunfield, P.F. The soil methane sink. In Greenhouse Gas Sinks; CABI: Wallingford, CT, USA, 2009; pp. 152–170. [Google Scholar] [CrossRef]
  42. Smith, E.J.; Davison, W.; Hamilton-Taylor, J. Methods for preparing synthetic freshwaters. Water Res. 2002, 36, 1286–1296. [Google Scholar] [CrossRef]
  43. Widdel, F.; Bak, F. Gram-Negative Mesophilic sulfate-reducing bacteria. In The Prokaryotes; Springer: Berlin, Germany, 1992; pp. 3352–3378. [Google Scholar]
  44. Whittenbury, R.; Phillips, K.C.; Wilkinson, J.F. Enrichment, isolation and some properties of methane-utilizing bacteria. J. Gen. Microbiol. 1970, 61, 205–218. [Google Scholar] [CrossRef] [Green Version]
  45. Kim, S.-J.; Park, S.-J.; Cha, I.-T.; Min, D.; Kim, J.-S.; Chung, W.-H.; Chae, J.-C.; Jeon, C.O.; Rhee, S.-K. Metabolic versatility of toluene-degrading, iron-reducing bacteria in tidal flat sediment, characterized by stable isotope probing-based metagenomic analysis. Environ. Microbiol. 2013, 16, 189–204. [Google Scholar] [CrossRef]
  46. Jung, M.-Y.; Park, S.-J.; Min, D.; Kim, J.-S.; Rijpstra, W.I.C.; Damsté, J.S.S.; Kim, G.-J.; Madsen, E.L.; Rhee, S.-K. Enrichment and characterization of an autotrophic ammonia-oxidizing archaeon of mesophilic crenarchaeal group I.1a from an agricultural soil. Appl. Environ. Microbiol. 2011, 77, 8635–8647. [Google Scholar] [CrossRef] [Green Version]
  47. Holmes, A.J.; Costello, A.; Lidstrom, M.E.; Murrell, J.C. Evidence that participate methane monooxygenase and ammonia monooxygenase may be evolutionarily related. FEMS Microbiol. Lett. 1995, 132, 203–208. [Google Scholar] [CrossRef]
  48. Hakobyan, A.; Zhu, J.; Glatter, T.; Paczia, N.; Liesack, W. Hydrogen utilization by Methylocystis sp. strain SC2 expands the known metabolic versatility of type IIa methanotrophs. Metab. Eng. 2020, 61, 181–196. [Google Scholar] [CrossRef] [PubMed]
  49. Lane, D. 16S/23S rRNA sequencing. In Nucleic Acid Techniques in Bacterial Systematics; Wiley: Hoboken, NJ, USA, 1991; pp. 115–175. [Google Scholar]
  50. Weisburg, W.G.; Barns, S.M.; Pelletier, D.A.; Lane, D.J. 16S ribosomal DNA amplification for phylogenetic study. J. Bacteriol. 1991, 173, 697–703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Miguez, C.; Bourque, D.; Sealy, J.; Greer, C.; Groleau, D. Detection and isolation of methanotrophic bacteria possessing soluble methane monooxygenase (sMMO) genes using the polymerase chain reaction (PCR). Microb. Ecol. 1997, 33, 21–31. [Google Scholar] [CrossRef] [PubMed]
  52. Hutchens, E.; Radajewski, S.; Dumont, M.G.; McDonald, I.R.; Murrell, J.C. Analysis of methanotrophic bacteria in Movile Cave by stable isotope probing. Environ. Microbiol. 2003, 6, 111–120. [Google Scholar] [CrossRef]
  53. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [CrossRef]
  54. Lewis, P.O.; Kumar, S.; Tamura, K.; Nei, M. MEGA: Molecular evolutionary genetics analysis, Version 1.02. Syst. Biol. 1995, 44, 576. [Google Scholar] [CrossRef]
  55. Yoon, S.-H.; Ha, S.-M.; Kwon, S.; Lim, J.; Kim, Y.; Seo, H.; Chun, J. Introducing EzBioCloud: A taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 2017, 67, 1613–1617. [Google Scholar] [CrossRef]
  56. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  57. Tatusova, T.; DiCuccio, M.; Badretdin, A.; Chetvernin, V.; Nawrocki, E.P.; Zaslavsky, L.; Lomsadze, A.; Pruitt, K.D.; Borodovsky, M.; Ostell, J. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 2016, 44, 6614–6624. [Google Scholar] [CrossRef]
  58. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.S.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [Green Version]
  59. Kanehisa, M.; Furumichi, M.; Tanabe, M.; Sato, Y.; Morishima, K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2016, 45, D353–D361. [Google Scholar] [CrossRef] [Green Version]
  60. Mistry, J.; Bateman, A.; Finn, R.D. Predicting active site residue annotations in the Pfam database. BMC Bioinform. 2007, 8, 1–14. [Google Scholar] [CrossRef] [Green Version]
  61. Lee, I.; Kim, Y.O.; Park, S.-C.; Chun, J. OrthoANI: An improved algorithm and software for calculating average nucleotide identity. Int. J. Syst. Evol. Microbiol. 2016, 66, 1100–1103. [Google Scholar] [CrossRef] [PubMed]
  62. Zhao, Y.; Wu, J.; Yang, J.; Sun, S.; Xiao, J.; Yu, J. PGAP: Pan-genomes analysis pipeline. Bioinformatics 2011, 28, 416–418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Eren, A.M.; Esen, Ö.C.; Quince, C.; Vineis, J.H.; Morrison, H.G.; Sogin, M.L.; Delmont, T.O. Anvi’o: An advanced analysis and visualization platform for omics data. PeerJ 2015, 3, e1319. [Google Scholar] [CrossRef] [PubMed]
  64. Alcock, B.P.; Raphenya, A.R.; Lau, T.T.Y.; Tsang, K.K.; Bouchard, M.; Edalatmand, A.; Huynh, W.; Nguyen, A.-L.V.; Cheng, A.A.; Liu, S.; et al. CARD 2020: Antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 2020, 48, D517–D525. [Google Scholar] [CrossRef]
  65. Lee, S.-P.; Park, J.-C. Management of organic matters by constructed treatment wetlands during rainfall events. J. Environ. Sci. Int. 2017, 26, 401–410. [Google Scholar] [CrossRef]
  66. Kwon, J.-H.; Han, Y.-S.; Cho, Y.-C.; Ahn, J.-S.; Yim, G.-J. Water Quality and Methane emission characteristics of aerobic wetlands constructed in coal mine area. J. Korean Soc. Miner. Energy Resour. Eng. 2018, 55, 371–382. [Google Scholar] [CrossRef]
  67. Im, J.; Lee, S.-W.; Yoon, S.; DiSpirito, A.A.; Semrau, J.D. Characterization of a novel facultative Methylocystis species capable of growth on methane, acetate and ethanol. Environ. Microbiol. Rep. 2010, 3, 174–181. [Google Scholar] [CrossRef]
  68. Reid, R.; Mosley, L.M. Comparative contributions of solution geochemistry, microbial metabolism and aquatic photosynthesis to the development of high pH in ephemeral wetlands in South East Australia. Sci. Total. Environ. 2016, 542, 334–343. [Google Scholar] [CrossRef]
  69. Dianou, D.; Adachi, K.; Dianou, D. Characterization of methanotrophic bacteria isolated from a subtropical paddy field. FEMS Microbiol. Lett. 1999, 173, 163–173. [Google Scholar] [CrossRef]
  70. Hou, C.T.; Laskin, A.I.; Patel, R.N. Growth and polysaccharide production by Methylocystis parvus OBBP on methanol. Appl. Environ. Microbiol. 1979, 37, 800–804. [Google Scholar] [CrossRef] [Green Version]
  71. Jo, S.Y.; Na Rhie, M.; Jung, S.M.; Sohn, Y.J.; Yeon, Y.J.; Kim, M.-S.; Park, C.; Lee, J.; Park, S.J.; Na, J.-G. Hydrogen production from methane by Methylomonas sp. DH-1 under micro-aerobic conditions. Biotechnol. Bioprocess Eng. 2020, 25, 71–77. [Google Scholar] [CrossRef]
  72. Dam, B.; Kube, M.; Dam, S.; Reinhardt, R.; Liesack, W. Complete sequence analysis of two methanotroph-specific repABC-containing plasmids from Methylocystis sp. strain SC2. Appl. Environ. Microbiol. 2012, 78, 4373–4379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Konstantinidis, K.T.; Tiedje, J.M. Towards a genome-based taxonomy for prokaryotes. J. Bacteriol. 2005, 187, 6258–6264. [Google Scholar] [CrossRef] [Green Version]
  74. Konstantinidis, K.T.; Ramette, A.N.; Tiedje, J.M. The bacterial species definition in the genomic era. Philos. Trans. R. Soc. B Biol. Sci. 2006, 361, 1929–1940. [Google Scholar] [CrossRef]
  75. Tatusov, R.L.; Koonin, E.V.; Lipman, D.J. A Genomic perspective on protein families. Science 1997, 278, 631–637. [Google Scholar] [CrossRef] [Green Version]
  76. Oshkin, I.Y.; Miroshnikov, K.K.; Grouzdev, D.S.; Dedysh, S.N. Pan-genome-based analysis as a framework for demarcating two closely related methanotroph genera Methylocystis and Methylosinus. Microorganisms 2020, 8, 768. [Google Scholar] [CrossRef]
  77. Stolyar, S.; Costello, A.M.; Peeples, T.L.; Lidstrom, M.E. Role of multiple gene copies in particulate methane monooxygenase activity in the methane-oxidizing bacterium Methylococcus capsulatus Bath. Microbiology 1999, 145, 1235–1244. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Matsen, J.B.; Yang, S.; Stein, L.Y.; Beck, D.A.; Kalyuzhanaya, M.G. Global molecular analyses of methane metabolism in methanotrophic alphaproteobacterium, Methylosinus trichosporium OB3b. Part I: Transcriptomic study. Front. Microbiol. 2013, 4, 40. [Google Scholar] [CrossRef] [Green Version]
  79. Cai, Y.; Zheng, Y.; Bodelier, P.L.E.; Conrad, R.; Jia, Z. Conventional methanotrophs are responsible for atmospheric methane oxidation in paddy soils. Nat. Commun. 2016, 7, 11728. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Baani, M.; Liesack, W. Two isozymes of particulate methane monooxygenase with different methane oxidation kinetics are found in Methylocystis sp. strain SC2. Proc. Natl. Acad. Sci. USA 2008, 105, 10203–10208. [Google Scholar] [CrossRef] [Green Version]
  81. Dunfield, P.F.; Liesack, W.; Henckel, T.; Knowles, R.; Conrad, R. High-affinity methane oxidation by a soil enrichment culture containing a type II methanotroph. Appl. Environ. Microbiol. 1999, 65, 1009–1014. [Google Scholar] [CrossRef] [Green Version]
  82. Anthony, C.; Ghosh, M.; Blake, C.C.F. The structure and function of methanol dehydrogenase and related quinoproteins containing pyrrolo-quinoline quinone. Biochem. J. 1994, 304, 665–674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Anthony, C.; Williams, P. The structure and mechanism of methanol dehydrogenase. Biochim. Biophys. Acta (BBA) Proteins Proteom. 2003, 1647, 18–23. [Google Scholar] [CrossRef]
  84. Nakagawa, T.; Mitsui, R.; Tani, A.; Sasa, K.; Tashiro, S.; Iwama, T.; Hayakawa, T.; Kawai, K. A Catalytic Role of XoxF1 as La3+-dependent methanol dehydrogenase in Methylobacterium extorquens strain AM1. PLoS ONE 2012, 7, e50480. [Google Scholar] [CrossRef] [Green Version]
  85. Keltjens, J.T.; Pol, A.; Reimann, J.; Camp, H.J.M.O.D. PQQ-dependent methanol dehydrogenases: Rare-earth elements make a difference. Appl. Microbiol. Biotechnol. 2014, 98, 6163–6183. [Google Scholar] [CrossRef]
  86. Picone, N.; Camp, H.J.M.O.D. Role of rare earth elements in methanol oxidation. Curr. Opin. Chem. Biol. 2019, 49, 39–44. [Google Scholar] [CrossRef]
  87. Krause, S.M.B.; Johnson, T.; Karunaratne, Y.S.; Fu, Y.; Beck, D.A.C.; Chistoserdova, L.; Lidstrom, M.E. Lanthanide-dependent cross-feeding of methane-derived carbon is linked by microbial community interactions. Proc. Natl. Acad. Sci. USA 2016, 114, 358–363. [Google Scholar] [CrossRef] [Green Version]
  88. Salisbury, S.A.; Forrest, H.S.; Cruse, W.B.T.; Kennard, O. A novel coenzyme from bacterial primary alcohol dehydrogenases. Nat. Cell Biol. 1979, 280, 843–844. [Google Scholar] [CrossRef]
  89. Anthony, C. The quinoprotein dehydrogenases for methanol and glucose. Arch. Biochem. Biophys. 2004, 428, 2–9. [Google Scholar] [CrossRef] [PubMed]
  90. Murrell, J.C.; Dalton, H. Nitrogen fixation in obligate methanotrophs. Microbiology 1983, 129, 3481–3486. [Google Scholar] [CrossRef] [Green Version]
  91. Auman, A.J.; Speake, C.C.; Lidstrom, M.E. nifH Sequences and nitrogen fixation in type I and type II methanotrophs. Appl. Environ. Microbiol. 2001, 67, 4009–4016. [Google Scholar] [CrossRef] [Green Version]
  92. Dedysh, S.N.; Ricke, P.; Liesack, W. NifH and NifD phylogenies: An evolutionary basis for understanding nitrogen fixation capabilities of methanotrophic bacteria. Microbiology 2004, 150, 1301–1313. [Google Scholar] [CrossRef] [Green Version]
  93. Bédard, C.; Knowles, R. Physiology, biochemistry, and specific inhibitors of CH4, NH4+, and CO oxidation by methanotrophs and nitrifiers. Microbiol. Rev. 1989, 53, 68–84. [Google Scholar] [CrossRef] [PubMed]
  94. Graf, E.-G.; Thauer, R.K. Hydrogenase from methanobacterium thermoautotrophicum, a nickel-containing enzyme. FEBS Lett. 1981, 136, 165–169. [Google Scholar] [CrossRef] [Green Version]
  95. Thauer, R.K.; Klein, A.A.R.; Hartmann, G.C. Reactions with molecular hydrogen in microorganisms: Evidence for a purely organic hydrogenation catalyst. Chem. Rev. 1996, 96, 3031–3042. [Google Scholar] [CrossRef]
  96. Pavlov, M.; Siegbahn, P.E.M.; Blomberg, M.R.A.; Crabtree, R.H. Mechanism of H−H activation by nickel−iron hydrogenase. J. Am. Chem. Soc. 1998, 120, 548–555. [Google Scholar] [CrossRef]
  97. Greening, C.; Biswas, A.; Carere, C.R.; Jackson, C.J.; Taylor, M.C.; Stott, M.B.; Cook, G.M.; Morales, S.E. Genomic and metagenomic surveys of hydrogenase distribution indicate H2 is a widely utilised energy source for microbial growth and survival. ISME J. 2015, 10, 761–777. [Google Scholar] [CrossRef] [Green Version]
  98. Cal, A.J.; Sikkema, W.D.; Ponce, M.I.; Franqui-Villanueva, D.; Riiff, T.J.; Orts, W.J.; Pieja, A.J.; Lee, C.C. Methanotrophic production of polyhydroxybutyrate-co-hydroxyvalerate with high hydroxyvalerate content. Int. J. Biol. Macromol. 2016, 87, 302–307. [Google Scholar] [CrossRef] [Green Version]
  99. Khosravi-Darani, K.; Mokhtari, Z.-B.; Amai, T.; Tanaka, K. Microbial production of poly(hydroxybutyrate) from C1 carbon sources. Appl. Microbiol. Biotechnol. 2013, 97, 1407–1424. [Google Scholar] [CrossRef] [PubMed]
  100. Byrom, D.P. Plastics from Microbes: Microbial Synthesis of Polymers and Polymer Precursors; Hanser: Munich, Germany, 1994; p. 5. [Google Scholar]
  101. Braunegg, G.; Sonnleitner, B.; Lafferty, R.M. A rapid gas chromatographic method for the determination of poly-?-hydroxybutyric acid in microbial biomass. Appl. Microbiol. Biotechnol. 1978, 6, 29–37. [Google Scholar] [CrossRef]
  102. Gunes, A.; Pilbeam, D.J.; Inal, A. Effect of arsenic–phosphorus interaction on arsenic-induced oxidative stress in chickpea plants. Plant Soil 2008, 314, 211–220. [Google Scholar] [CrossRef]
  103. Slyemi, D.; Bonnefoy, V. How prokaryotes deal with arsenic†. Environ. Microbiol. Rep. 2011, 4, 571–586. [Google Scholar] [CrossRef]
  104. Herath, I.; Vithanage, M.; Bundschuh, J.; Maity, J.P.; Bhattacharya, P. Natural arsenic in global groundwaters: Distribution and geochemical triggers for mobilization. Curr. Pollut. Rep. 2016, 2, 68–89. [Google Scholar] [CrossRef] [Green Version]
  105. Patil, Y.B.; Paknikar, K. Development of a process for biodetoxification of metal cyanides from waste waters. Process. Biochem. 2000, 35, 1139–1151. [Google Scholar] [CrossRef]
  106. Gupta, N.; Balomajumder, C.; Agarwal, V.K. Enzymatic mechanism and biochemistry for cyanide degradation: A review. J. Hazard. Mater. 2010, 176, 1–13. [Google Scholar] [CrossRef]
  107. Dash, H.R.; Das, S. Bioremediation of mercury and the importance of bacterial mer genes. Int. Biodeterior. Biodegrad. 2012, 75, 207–213. [Google Scholar] [CrossRef]
  108. Vorobev, A.; Jagadevan, S.; Baral, B.S.; DiSpirito, A.A.; Freemeier, B.C.; Bergman, B.H.; Bandow, N.L.; Semrau, J.D. Detoxification of mercury by methanobactin from Methylosinus trichosporium OB3b. Appl. Environ. Microbiol. 2013, 79, 5918–5926. [Google Scholar] [CrossRef] [Green Version]
  109. Shi, L.-D.; Chen, Y.-S.; Du, J.-J.; Hu, Y.-Q.; Shapleigh, J.P.; Zhao, H.-P. Metagenomic evidence for a Methylocystis species capable of bioremediation of diverse heavy metals. Front. Microbiol. 2019, 9, 3297. [Google Scholar] [CrossRef]
  110. Boden, R.; Murrell, J.C. Response to mercury (II) ions in Methylococcus capsulatus (Bath). FEMS Microbiol. Lett. 2011, 324, 106–110. [Google Scholar] [CrossRef] [Green Version]
  111. Lu, X.; Gu, W.; Zhao, L.; Haque, M.F.U.; DiSpirito, A.A.; Semrau, J.D.; Gu, B. Methylmercury uptake and degradation by methanotrophs. Sci. Adv. 2017, 3, e1700041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The phylogenetic tree of Methylocystis sp. B8 and other strains in the family Methylocystaceae based on concatenated conserved proteins. The outgroup was Rhizobium leguminosarum bv. trifolii in the family Rhizobiaceae. Closed circles indicate the bootstrap values ≥70%.
Figure 1. The phylogenetic tree of Methylocystis sp. B8 and other strains in the family Methylocystaceae based on concatenated conserved proteins. The outgroup was Rhizobium leguminosarum bv. trifolii in the family Rhizobiaceae. Closed circles indicate the bootstrap values ≥70%.
Microorganisms 08 01719 g001
Figure 2. Phylogenetic tree of PmoA of strain B8 and related strains. The phylogenetic tree was reconstructed using the neighbor-joining method. The outgroup was AmoA of Nitrosolobus multiformis Mm24C. The clades (right side of the tree) of PmoA was classified according to Cai et al. [79]. Bootstrap values below 70% are not shown. Bar: 0.05 substitutions per amino acid site.
Figure 2. Phylogenetic tree of PmoA of strain B8 and related strains. The phylogenetic tree was reconstructed using the neighbor-joining method. The outgroup was AmoA of Nitrosolobus multiformis Mm24C. The clades (right side of the tree) of PmoA was classified according to Cai et al. [79]. Bootstrap values below 70% are not shown. Bar: 0.05 substitutions per amino acid site.
Microorganisms 08 01719 g002
Figure 3. Phylogenetic tree of hydrogenase based on the amino acid sequences of Methylocystis sp. B8 and other strains. The tree was constructed with MEGA7 using the neighbor-joining method. These enzymes were classified into eight groups (3b, 3d, 1h/5, 1d, 2c, 2b, 1a, and 1i). Bootstrap values below 70% are not shown. Bar: 0.2 substitutions per amino acid site.
Figure 3. Phylogenetic tree of hydrogenase based on the amino acid sequences of Methylocystis sp. B8 and other strains. The tree was constructed with MEGA7 using the neighbor-joining method. These enzymes were classified into eight groups (3b, 3d, 1h/5, 1d, 2c, 2b, 1a, and 1i). Bootstrap values below 70% are not shown. Bar: 0.2 substitutions per amino acid site.
Microorganisms 08 01719 g003
Table 1. The properties of sampling water.
Table 1. The properties of sampling water.
PropertiesSampling Water
Sampling date29 Jun 2017
Temperature (°C)26.1
pH6.46
EC (μS/cm)670
DO (mg/L)6.82
ORP (mV)84.7
Ca2+ (mg/L)23.58
K+ (mg/L)16.43
Mg2+ (mg/L)6.49
Na+ (mg/L)83.71
Fluoride (mg/L)3.19
Chloride (mg/L)151.56
Nitrite (mg/L)1.7
Nitrate (mg/L)26.48
Sulfate (mg/L)33.22
Table 2. Characteristics of Methylocystis sp. B8 and related strains in Methylocystis.
Table 2. Characteristics of Methylocystis sp. B8 and related strains in Methylocystis.
1234567
Cell size (µm)0.7 × 1.50.3–0.6 × 0.7–1.00.8–1.1 × 1.1–2.50.5 × 0.5–2.00.5 × 0.5–1.00.9–1.4 × 1.8–3.40.8–1.2 × 1.4–4.0
Cell shapeRodDumbbellRodCoccobacillus, rodCoccobacillus, rodSmall, curved coccoids, shortstraight or regularly curved rod, ovoid
Colony colorCreamCreamPink-redWhite, buffDiffusion brown, pale pinkLight creamWhite
Temperature
(°C)
15–37<375–3720–28~20–378–375–30
Optimum
Temperature (°C)
30–35302725–2830–3725–3025
pH6.5–10.04.0–9.05.5–9.05.0–95.5–8.84.2–7.64.4–7.5
Optimum pH8.5–9.07.0NR6.5–77.56.0–6.55.8–6.2
NaCl (%)0–0.5NR0=<0–1.00–2.00–0.10–0.5
Optimum
NaCl (%)
0–0.5NRNRNRNRNRNR
Oxidase+++++++
Catalase+++++++
Methanol (%)0.05–1+-0.25≤0.15≤1
Optimum
Methanol (%)
0.1NR-NRNRNR0.1
Carbon Source
Ethanol-NR-NRNR±NR
Acetate+w *-w *-+w *
Pyruvate-NRNRNRNR±NR
Nitrogen Source
Cysteine+NRNR--+-
Yeast extract+NRNR++++
Casamino acid-NRNR++NRNR
Betaine-NRNR--NRNR
Carsitone-NRNRNRNRNRNR
(NH4)2SO4+NR+++++
Nitrate+++++++
Strains: 1, Methylocystis sp. B8 (this study); 2, Methylocystis hirsuta CSC1T [31]; 3, Methylocystis rosea SV97T [30]; 4, Methylocystis echinoides IMET 10491T [28,69]; 5, Methylocystis parvus OBBPT [28,70]; 6, Methylocystis bryophila H2sT [32,39]; 7, Methylocystis heyeri H2T [33]. All strains were gram-negative and non-motile. The characteristics of B8 and other species, Methylocystis, are provided. +, Positive reaction; -, negative reaction; w, weakly positive; ±, trace growth; NR, not reported. * Data from Belova et al. [39].
Table 3. Genomic features of Methylocystis sp. B8 and other strains of Methylocystis.
Table 3. Genomic features of Methylocystis sp. B8 and other strains of Methylocystis.
1234567891011
Genome size (bp)3,409,1644,726,0344,213,0433,773,4443,912,0504,475,9123,359,9404,707,9714,690,5664,364,3633,644,633
GC content59.962.862.463.462.563.463.163.16362.562.7
N50 (bp)386,5673,131,8073,776,0273,773,4441,614,04095,607917.8054,532,9503,287,239136,57154,643
Coverage (×)132.630252.053unknown33280.034.753.0602.0
Contig number287414108921255158
Protein coding genes31764307403635833709416035764148398740433433
rRNA
operons
363333NR65316
tRNA
genes
475349475046NR47545045
Genbankaccession numberVBTZ00000000AEVM00000000QWDD00000000HE956757ARCT00000000AJTV00000000PHSQ00000000CP019948SOPH00000000PYDU00000000AYNA00000000
Strains: 1, Methylocystis sp. B8; 2, Methylocystis sp. ATCC 49242; 3, M. hirsuta CSC1T; 4, Methylocystis sp. SC2; 5, M. rosea SV97T; 6, M. parvus OBBPT; 7, Methylocystis sp. KS32; 8, M. bryophila H2sT; 9, M. heyeri H2T; 10, Methylocystis sp. MitZ-2018; 11, Methylocystis sp. SB2.
Table 4. Genes related to heavy metal detoxification in Methylocystis sp. B8 and related species.
Table 4. Genes related to heavy metal detoxification in Methylocystis sp. B8 and related species.
Heavy Metal Related Genes123456789101112
Mercury reductase (merA)+++++-+---+-
Arsenate reductase (arsC)++++++++++++
Arsenite oxidase---+--------
Arsenite methyltransferase------+-----
Arsenite efflux transporter+---+-------
Cobalt-zinc-cadmium resistance protein-+---++--+-+
Tellurite resistance protein+-+-+---+---
Copper resistance protein+++++++++-++
Strains: 1, Methylocystis sp. B8; 2, Methylocystis sp. ATCC 49242; 3, M. hirsuta CSC1T; 4, Methylocystis sp. SC2; 5, M. rosea SV97T; 6, M. parvus OBBPT; 7, Methylocystis sp. KS32; 8, M. bryophila H2sT; 9, M. heyeri H2T; 10, Methylocystis sp. LW5; 11, Methylocystis sp. MitZ-2018; 12, Methylocystis sp. SB2. +, positive; -, negative.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jung, G.-Y.; Rhee, S.-K.; Han, Y.-S.; Kim, S.-J. Genomic and Physiological Properties of a Facultative Methane-Oxidizing Bacterial Strain of Methylocystis sp. from a Wetland. Microorganisms 2020, 8, 1719. https://doi.org/10.3390/microorganisms8111719

AMA Style

Jung G-Y, Rhee S-K, Han Y-S, Kim S-J. Genomic and Physiological Properties of a Facultative Methane-Oxidizing Bacterial Strain of Methylocystis sp. from a Wetland. Microorganisms. 2020; 8(11):1719. https://doi.org/10.3390/microorganisms8111719

Chicago/Turabian Style

Jung, Gi-Yong, Sung-Keun Rhee, Young-Soo Han, and So-Jeong Kim. 2020. "Genomic and Physiological Properties of a Facultative Methane-Oxidizing Bacterial Strain of Methylocystis sp. from a Wetland" Microorganisms 8, no. 11: 1719. https://doi.org/10.3390/microorganisms8111719

APA Style

Jung, G. -Y., Rhee, S. -K., Han, Y. -S., & Kim, S. -J. (2020). Genomic and Physiological Properties of a Facultative Methane-Oxidizing Bacterial Strain of Methylocystis sp. from a Wetland. Microorganisms, 8(11), 1719. https://doi.org/10.3390/microorganisms8111719

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