Enrichment, Isolation, and Identification
Different colored colonies appeared on solid media plates after 3–4 weeks of incubation. A total of four methanotroph strains were cultivated and named with the designation WWC (wetland water creek/canal). Of the four, three: WWC1, WWC3, and WWC5, showed 99% 16S rRNA similarity with Methylomonas koyamae and were classified as Methylomonas koyamae after doing 16S rRNA gene sequencing and blast. Strain WWC4 showed only 97.81% similarity to Methylomonas koyamae Fw12E-YT and hence was further subjected to polyphasic characterization being a putatively novel species of Methylomonas.
Microscopic Identification
Strain WWC4 formed pink to slightly orange-colored colonies on the NMS media plate with 2% agarose as a solidifying agent (Fig. 1). The diameter of the colonies was 2–4 mm. The live cells were thick short motile rods 2.5-3 µm long and 0.8–1.2 µm (Fig. 1). This strain formed a thin pinkish-orange turbidity in liquid NMS media with methane and air (20:80) as headspace gas (Fig. 1). The cells showed Gram Negative character.
Phylogenetic and Phylogenomic Affiliation
The 16S rRNA gene of Methylomonas sp. WWC4 (MH764454.1) showed that Methylomonas sp. WWC4 showed the closest affiliation to Methylomonas koyamae Fw12E-YT with 97.61% similarity, followed by Methylomonas rhizoryzae GJ1T 96.02% similarity, respectively, after nucleotide BLAST search. The maximum likelihood tree of the 16S rRNA gene of strain WWC4 showed a similar picture with three of these species as the closest neighbors (Fig. 2). The 16S rRNA gene of strain WWC4 and Methylomonas strain R-45383 shared 99.57% similarity, though this strain is not yet classified as a member of a new species. Methylomonas R-45383 was isolated from wetland (Ghent, Belgium) and has not been described and validated yet though the genome is described and available for comparison [25]. The pmoA gene of Methylomonas sp. WWC4 (MH806338) showed the closest affiliation to type strains; Methylomonas koyamae Fw12E-YT and Methylomonas rhizoryzae GJ1T with 89.43% and 88.78% similarity, respectively, and this can be seen by its close phylogenetic grouping with these two species within the genus Methylomonas (Fig. 3). The genomic comparison of the genome of strain WWC4 with its closest members of the genus Methylomonas was carried out using the PATRIC webserver. The phylogenomic tree of strain WWC4 showed that it grouped with Methylomonas rhizoryzae GJ1T and Methylomonas koyamae Fw12E-YT (Fig. 4).
Draft Genome Features:
The genome of strain WWC4 was sequenced using Illumina sequencing technology, which produced a total of 23,037,313 paired-end reads with average read length distribution of 141 bp and quality score distribution of 36. The 21,606,919 high-quality paired-end reads were retained after low-quality bases and adapter sequence filtration using Trimmomatic (v 0.35) and cutadapt (v 1.18). A total of 73 high-quality reads were assembled using the Soapdenovo assembler. A total of 56 scaffolds of > 500 bp were constructed, with an N50 of 223 kb. The largest scaffold assembled measured 691.7 kb. The G + C content of the draft genome was 55.9%, and the genome contains 4,672 CDS (Table 1). A total of 4,559 genes were annotated using BlastX (Table 1). These ~ 56 contigs, which belong to the genus Methylomonas were uploaded in RAST (http://rast.nmpdr.org/) and NCBI-PGAP for genome annotation (Table 1).
The nearest valid member of Methylomonas sp. WWC4 is Methylomonas koyamae Fw12E-YT with 74.45%, 75.72%, and 21.5 [19.3–23.9%] of AAI, ANIb, and dDDH values, respectively (Table 2).
Salient metabolic characteristics predicted from the draft genome
The complete details of the enzymes related to the carbon metabolism of strain WWC4 have been provided in Supplementary Table 1. Strain WWC4 contained two sets of particulate methane monooxygenase and one set of soluble methane monooxygenase for the oxidation of methane to methanol. The oxidation of methanol into formaldehyde was pursued by the presence of the methanol dehydrogenase enzyme present in the genome. Three copies of the PQQ-dependent dehydrogenase enzyme were present in the genome, out of which one of the enzymes was novel to methanotrophs and showed the closest similarity to Gemmatimonadetes bacterium. One copy of each NAD(P)-, zinc, and iron-dependent methanol dehydrogenase was also found in the genome of stain WWC4. All sets of enzymes for formaldehyde to formate and formate to carbon-dioxide conversion were present in the genome (Supplementary Table 1). The delta subunit of formate dehydrogenase was novel to methanotrophs and showed the closest similarity to Bradyrhizobium algeriense RST91.
The details of nitrogen metabolism pathway genes present in the WWC4 genome has provided in Supplementary Table 1. The genes required for environmental nitrogen fixation were found in strain WWC4. The nitrogen fixation proteins NifT, X, Q, Z, W, and NifM, were found in the genome. The presence of classical nitrogenase molybdenum-iron proteins (NifH,D, K) was also detected in the genome. Within the nitrate reduction enzymes, nitrite reductase, nitrate reductase, and nitric oxide reductase were present in the genome (Supplementary Table 1).
Other important proteins involved in hemerythrin metabolism, and these hemerythrin genes are essential in oxygen scavenging conditions and helpful in transporting oxygen to particulate methane monooxygenase, Supplementary Table 1 [26–28].
Physiology and chemotaxonomy
Methane and methanol are the only substrates used by strain WWC4 as carbon sources. This strain cannot utilize multicarbon (0.1%) such as glucose, fructose, sucrose, maltose, xylose, arabinose, and raffinose. Neither formate nor formaldehyde was utilized by strain WWC4. It can tolerate up to 0.5% of methanol present in the growth media. The concentration of 0.05% of ammonium chloride, glutamate, peptone, yeast extract, and lysine can be used as a nitrogen source. It also can fix environmental nitrogen when growing optimally in nitrogen-free media. Strain WWC4 was a mesophilic methanotroph to grow in the 15–30ºC temperature range with an optimum 25ºC at a pH range of 3–9 under a buffered medium. It can tolerate up to 1% of NaCl-salinity in the medium. The comparisons of some of the significant characteristics of strain WWC4 with other members of the Methylomonas genus are enlisted in Table 3 [13, 29].
The cell wall fatty acids profile of strain WWC4 enclosed maximum amounts of 15:1 ω8c (21.09%), 16:0 3 OH (15.7%), and 16:1ω5c (12.6%). The detail of the cell wall fatty acids profile of strain WWC4 and type strains of species of the Methylomonas genus has provided in Table 4 [13, 29]. The FAME profile of WWC4 is unique as compared to the other Methylomonas species (Table 4). Strain WWC4 did not withstand cryopreservation and could not be deposited in two international culture collections. The culture is part of our institutional, WDCM-approved culture collection as MCMB-1474, maintained live, since its isolation (from 2018- till date). It can be stored in the fridge at 4–8°C for 3–4 months in liquid and has to be sub-cultured after intervals. Similar amount of survival is seen when maintained on agarose plates in methane-air environment.
Carotenoid pigment genes:
Genus Methylomonas has been reported to produce carotenoid pigment [30–32]. Strain WWC4 showed distinct pink color indicative of carotenoid pigment and therefore, the genome of strains was analyzed for the carotenoid production pathway. The production of carotenoid in Methylomonas strains is done by the following pathway and the enzymes were detected in the draft genome (Supplementary Fig. 1and 2). In this case, the strain produces carotenoids by utilizing two molecules of C15 farnesylpyrophosphate (FPP), resulting in C30 carbon skeleton formation (ref). Terpenoid is synthesised through the MEP pathway. The critical enzyme 1-deoxy-D-xylulose 5-phosphate synthase (dxs) genes of MEP/DOXP pathways was present in strain WWC4. MEP/DOXP pathways synthesise isopentenyl diphosphate (IPP) that were converted by farnesyl diphosphate synthase (FPPS/IspA) into farnesylpyrophosphate (FPP) and geranylpyrophosphate (GGPP). Geranaylgeranyl pyrophosphate synthase would convert IPP and FPP into GGPP, which will be sequentially converted into phytoene by the enzymatic activity of squalene/ phytoene synthase (crtB).
Plant growth promotion related genes:
As methanotrophs are associated with rice plant rhizospheres they can be applied as agents to promote plant growth (Rahalkar, Monali, unpublished data). The genome of WWC4 was explored for all different types of genes including the Plant growth promotion (PGP) genes: genes for IAA production, nitrogen fixation, phosphate solubilization. were analyzed. PGP-properties in a bacterium or facilitates the plant nutrient uptake from the surrounding environment [33]. The most essential nutrient for plants is nitrogen for its development [34]. Methanotrophs are present in large numbers near rice rhizospheres [35–37] and nitrogen fixation genes from methanotrophs are detected in the transcribed genes [38] The 2nd most important nutrient for plants after nitrogen is phosphorous, which is in an insoluble form [34]. Therefore, it becomes a major growth-limiting factor in agriculture [34]. Methanotrophs contain phosphate solubilising enzymes which make the soluble phosphates available to the plants for their development. The important PGP genes were also found in the selected methanotrophic strains, such as glycoside hydrolase, alkaline phosphatase, phytase, and trehalose enzymes, to improve soil fertility and plant development (Supplementary Table 1).
Classical nitrogen fixation pathway genes (NifDKH operon) were detected in the genomes for environmental N2-fixation. From the denitrification pathway, nitrite reductase nirBD genes and nitric oxide reductase norB and norC genes were present (Supplementary Fig. 3).
Description of “ Candidatus Methylomonas sedimenticola” strain WWC4
A putative novel methanotroph strain WWC4 was isolated from a mud sample of a shallow freshwater part of a wetland of freshwater creek region, Alibag, near Nagaon beach, Maharashtra, India (18o26'30"N 72o54'20"E). Strain WWC4 showed 97.81% 16S rRNA and 89.46% pmoA gene similarity with Methylomonas koyamae Fw12E-YT. The genomic comparison showed 74.45%, 75.72%, and 21.5 [19.3–23.9%] of AAI, ANIb, and DDH values, respectively, with its closest member Methylomonas koyamae Fw12E-YT. The nearest strain is Methylomonas R45383 (genome accession number LUU01).
The genomic evaluation and phylogenetic analysis concluded the WWC4 strain represents a putative novel methanotroph species. This putative novel species was isolated from a muddy-water sample and was phylogenetically related to the Methylomonas genus; therefore, the name is tentatively proposed as a “CandidatusMethylomonas sedimenticola”sp. nov. strain WWC4. It shows 97.81% AAI and 97.01% ANIb and 77.3% DDH value with the not yet validly described strain: Methylomonas R45383, hence based on these closeness, the strain R45383 also could be a second member of the new Candidatus species.
Like almost all methanotrophs, this strain WWC4 can only utilize methane and methanol as carbon sources. The RuMP and serine pathway can be used for methane oxidation, and the presence of genes for the same was found in its genomes. While the nitrogen source, ammonium chloride, glutamate, peptone, yeast extract, and lysine can be used to grow the strain, the environmental nitrogen can be fixed by strain WWC4, which has been checked by laboratory experiments as well as the presence of genes in the genome.
Candidatus Methylomonas sedimenticola” strain WWC4 is a mesophilic methanotroph and grows in a 15–30ºC temperature range at a 3–10 pH range. The optimum temperature and pH were 25ºC and 6.8, respectively. The major fatty acids were 15:1 ω8c (21.09%), 16:0 3OH (15.7%), and 16:1ω5c (12.6%). The whole-genome shotgun project is deposited at DDBJ/ENA/GenBank, and the accession number is JAATWI0.1. The “CandidatusMethylomonas sedimenticola” strain WWC4 is maintained in our in-house WDCM-approved culture collection, (MACS collection of microorganisms), as a live culture, and designated as MCMB-1474.
Table 1: Draft genome details of strain WWC4
|
Draft genome traits
|
WWC4
|
NCBI Accession number
|
JAATWI0.1
|
Bio-project number
|
PRJNA520977
|
Bio-sample
|
SAMN10613771
|
Assembly
|
ASM1188220v1
|
Genome size
|
5.19 Mbp
|
G+C content
|
55.9%
|
N50
|
161,891 bp
|
Genes
|
4,560
|
Proteins
|
4,407
|
tRNA, nc RNA
|
39, 4
|
rRNA
|
2
|
Contigs
|
56
|
Table 2: AAI, ANIb, & DDH calculation of strain WWC4 with its closely related members
|
Closely related members
|
AAI value
|
ANIb value
|
dDDH value
|
Methylomonas strain WWC4
|
-
|
-
|
-
|
Methylomonas sp. R-45383
|
97.18%
|
97.01%
|
77.3 [74.3 - 80%]
|
Methylomonas koyamae Fw12E-YT
|
74.45%
|
75.72%
|
21.5 [19.3-23.9%]
|
Methylomonas methanica S1T
|
74.37%
|
74.77%
|
20.6 [18.4-23.1]
|
Methylomonas rhizoryzae GJ1T
|
69.06%
|
72.58%
|
21.3 [19-23.7%]
|
Table 3: Comparative Table of Methylomonas species with strain WWC4. (Included strains: Methylomonas koyamae Type strain: Fw12E-Y, JCM 16701, NBRC 105905, NCIMB 14606; Methylomonas lenta Type strain: R-45377, LMG 26260, JCM 19378; Methylomonas methanica Type strain: ACM 3307, ATCC 35067, IMET 10543, NCIMB 11130, VKM B-2110).
Characteristics
|
Strain WWC4
|
Methylomonas methanica
|
Methylomonas koyamae
|
Methylomonas rhizoryzae GJ1
|
Methylomonas paludis
|
Methylomonas lenta
|
Cell Morphology
|
Rods
|
Rods
|
Rods
|
Rods
|
Rods
|
Rods
|
Cell size
|
2.5-3 µm X
0.8-1.2 µm
|
0.5-3 µm X 0.5-1 µm
|
1.2-2.5 µm X 0.8-1.1 µm
|
0.6–1.0 µm X
0.7–2.3 µm
|
1-4 µm X
1-1.5 µm
|
1.3-2 µm X 0.6-0.9 µm
|
Motility
|
Yes
|
Yes
|
Yes
|
Yes
|
No
|
Yes
|
Surface pellicle
|
Yes
|
Yes
|
No
|
Yes
|
No
|
Not detected
|
Pigmentation
|
Pink to orange
|
Pink to orange
|
Pink to orange
|
Pink
|
Pale pink
|
White to pink sheen
|
Optimum growth temperature (ºC)
|
25 (range 15-30)
|
25-30 (range 10-35)
|
30 (range 10-40)
|
25-33 (range 16-37)
|
20-25 (range 8-30)
|
20-25 (range 15-28)
|
Optimum pH
|
6.8 (range 3-10)
|
7 (range 5.5-9)
|
6.5 (range 5.5-7)
|
6-8 (range 5.5-8.5)
|
5-8-6.4 (range 3.8-7.3)
|
6.8-7.3 (range 6.3-7.8)
|
G+C Content
|
55.9%
|
51-54%
|
57%
|
53.87%
|
48-49%
|
47%
|
Major PLFAs
|
C15:1 ω8c, 16:0 3OH
|
C14:0,
C16:1 ω8c
|
C14:0,
C16:1 ω8c
|
C16:1, C14:0
|
C16:1 ω5t, C16:1 ω8c
|
C16:1 ω8c, C16:1 ω5c
|
Table 4: Comparison of fatty acid methyl esters (FAME) profiles in cell walls of various Methylomonas species with strain WWC4.
Fatty acids
|
Strain WWC4
|
Methylomonas methanica
Methylomonas aurantiaca
Methylomonas fodinarum
|
Methylomonas koyamae
|
Methylomonas rhizoryzae GJ1
|
Methylomonas lenta
|
Methylomonas paludis
|
9:0
|
0.81
|
-
|
-
|
-
|
-
|
-
|
10:0
|
1.24
|
-
|
-
|
-
|
-
|
-
|
11:0
|
1.52
|
-
|
-
|
-
|
-
|
-
|
12:0
|
3.20
|
-
|
-
|
0.53
|
0.6-2.5
|
-
|
12:0 aldehyde
|
1.91
|
-
|
-
|
-
|
-
|
-
|
13:0
|
2.12
|
-
|
-
|
0.05
|
0.6-0.9
|
-
|
14:0
|
9.99
|
18.9-24.6
|
23
|
20.2
|
6.4-9.8
|
11.8
|
15:0
|
-
|
0-1.2
|
1.2
|
-
|
5.3-5.8
|
0.5
|
15:0 iso
|
-
|
0-2.5
|
-
|
-
|
-
|
-
|
15:0 anteiso
|
-
|
0-2.4
|
-
|
-
|
-
|
-
|
14:1 ω5c
|
4.23
|
-
|
-
|
-
|
-
|
-
|
15:1 ω8c
|
21.09
|
-
|
-
|
-
|
1.0-2.3
|
-
|
15:1 ω6c
|
3.04
|
-
|
-
|
-
|
0.5-0.6
|
-
|
15:1 ω5c
|
7.34
|
-
|
-
|
-
|
1.0-1.1
|
-
|
16:0
|
5.22
|
4.3-8.7
|
7.7
|
6.37
|
5.0
|
5.6
|
16:1ω11c
|
-
|
-
|
-
|
-
|
-
|
-
|
16:1ω8c
|
-
|
18.7-41.3
|
39.4
|
-
|
40.8-42.4
|
22.1
|
16:1ω7c
|
-
|
7.7-15.3
|
4.35
|
-
|
9.1-10.5
|
13.9
|
16:1ω6c
|
-
|
4.5-13.3
|
-
|
-
|
-
|
5.0
|
16:1ω7c/ω6c (sum)
|
8.55
|
-
|
-
|
54.42
|
-
|
-
|
16:1ω5c
|
12.60
|
1.9-6.3
|
16.7
|
16.33
|
11.7-18.3
|
1.8
|
16:1ω5t
|
-
|
7.9-16.6
|
-
|
-
|
-
|
34.8
|
17:0 cyclo
|
-
|
0-2.1
|
-
|
-
|
-
|
-
|
17:1 ω8c
|
-
|
-
|
-
|
-
|
0-0.8
|
-
|
17:1 ω7c
|
-
|
0-0.7
|
-
|
-
|
-
|
-
|
17:1 ω7t
|
-
|
0-0.3
|
-
|
-
|
-
|
-
|
18:0
|
-
|
0-0.1
|
0.72
|
-
|
-
|
1.2
|
18:1 ω7c
|
-
|
0.2-2.5
|
-
|
-
|
-
|
-
|
18:1 ω6c
|
-
|
-
|
-
|
-
|
-
|
-
|
18:1 ω5c
|
-
|
0-0.2
|
1.7
|
-
|
-
|
-
|
19:0 cyclo
|
-
|
0.2-0.4
|
-
|
-
|
-
|
-
|
19:1 branched
|
-
|
0-0.5
|
-
|
-
|
-
|
-
|
10:0 3 OH
|
1.44
|
-
|
-
|
-
|
-
|
-
|
15:0 2 OH
|
-
|
-
|
1.43
|
-
|
-
|
-
|
16:0 3 OH
|
15.70
|
-
|
3.8
|
1.68
|
4.1-4.2
|
-
|