Functional Analysis of Methylomonas sp. DH-1 Genome as a Promising Biocatalyst for Bioconversion of Methane to Valuable Chemicals

: Methylomonas sp. DH-1, newly isolated from the activated sludge of a brewery plant, has been used as a promising biocatalytic platform for the conversion of methane to value-added chemicals. Methylomonas sp. DH-1 can efﬁciently convert methane and propane into methanol and acetone with a speciﬁc productivity of 4.31 and 0.14 mmol/g cell/h, the highest values ever reported, respectively. Here, we present the complete genome sequence of Methylomonas sp. DH-1 which consists of a 4.86 Mb chromosome and a 278 kb plasmid. The existence of a set of genes related to one-carbon metabolism and various secondary metabolite biosynthetic pathways including carotenoid pathways were identiﬁed. Interestingly, Methylomonas sp. DH-1 possesses not only the genes of the ribulose monophosphate cycle for type I methanotrophs but also the genes of the serine cycle for type II. Methylomonas sp. DH-1 accumulated 80 mM succinate from methane under aerobic conditions, because DH-1 has 2-oxoglutarate dehydrogenase activity and the ability to operate the full TCA cycle. Availability of the complete genome sequence of Methylomonas sp. DH-1 enables further investigations on the metabolic engineering of this strain for the production of value-added chemicals from methane. pyrophosphate; IspD: synthase;


Introduction
Methane is the principal component of natural/shale gas and biogas, and recently, has attracted much attention as a chemical feedstock [1]. The chemical conversions of methane to other chemicals generally require the input of high amounts of energy because of the high stability of the carbon-hydrogen bond (C-H bond), while the biological conversion of methane to chemicals using methanotrophs can be conducted in ambient conditions [2]. In addition, the bioconversion of methane showed higher conversion of up to 75% [1,3].
In an effort to develop a biocatalytic platform strain, we have isolated a type I methanotroph, Methylomonas sp. DH-1, from the activated sludge of a brewery plant [4]. Methylomonas sp. DH-1 was reported as a highly efficient biocatalyst for the bioconversion of methane to methanol, which can be directly used as alternative fuels, antifreeze and as a precursor to other compounds, with a specific productivity of 4.31 mmol/g cell/h [4]. Furthermore, Methylomonas sp. DH-1 has high potential for productivity of 4.31 mmol/g cell/h [4]. Furthermore, Methylomonas sp. DH-1 has high potential for methanol production due to its high tolerance to methanol of up to 7% (v/v) [4]. Methylomonas sp. DH-1 has also been evaluated for its catalytic capability to convertpropane to acetone, which is used as an industrial solvent for polymers, in acetylene storage and in the pharmaceutical industry [5]. Moreover, the accumulation of acetone in the absence of chemical inhibitors is advantageous for biocatalytic gas-to-liquid conversion technology. Additionally, Methylomonas sp. DH-1 can produce yellow-to-red pigments which are expected to be numerous carotenoids (unpublished report). Thus, Methylomonas sp. DH-1 can be an important biocatalyst for methane bioconversion to chemicals/fuels. In this study, we sequenced, assembled and annotated the whole genome sequence of Methylomonas sp. DH-1 as the first step for the development of a methanotrophic platform strain. Methylomonas sp. DH-1 was also used for the production of succinate as a model compound from methane.

Genome Statistics and General Features
The complete genome of Methylomonas sp. DH-1 consists of a circular chromosome of 4,849,532 bp (56.5% G + C) and a plasmid of 277,875 bp (51.7% G + C) ( Figure 1).

Figure 1.
A circular representation of the Methylomonas sp. DH-1 chromosome and plasmid. The rings represent (from inner to outer) the nucleotide position ruler, GC skew, %GC, coding sequences (CDSs) transcribed in the counterclockwise direction, and the clockwise direction, respectively. The CDSs are colored according to the assigned clusters of orthologous genes (COG) classes. The image was rendered using the CGView server (http://stothard.afns.ualberta.ca/cgview_server/). Methylomonas sp. DH-1 was shown to be phylogenetically closely related to Methylomonas koyamae Fw12E-Y T based on 16S sequence similarity [4]. Electronic DNA-DNA hybridization (DDH) estimate between DH-1 and Fw12E-Y T (=JCM 16701 T ), calculated by the Genome-to-Genome Distance Calculator (http://ggdc.dsmz.de/distcalc2.php), was 73.9%, which suggests that DH-1 belongs to the Methylomonas koyamae species, while average nucleotide identity (ANI) between these two strains was calculated to be 97.76% using JSpecies [6]. Moreover, MUMMER whole-genome alignment [7] between DH-1 and Fw12E-Y T showed the close similarity of these two strains, aligning 298 out of 382 scaffolds of the Fw12E-Y T genome assembly (96.68% of the total length) on the DH-1 reference genome sequence (Figure 2A).

Figure 1.
A circular representation of the Methylomonas sp. DH-1 chromosome and plasmid. The rings represent (from inner to outer) the nucleotide position ruler, GC skew, %GC, coding sequences (CDSs) transcribed in the counterclockwise direction, and the clockwise direction, respectively. The CDSs are colored according to the assigned clusters of orthologous genes (COG) classes. The image was rendered using the CGView server (http://stothard.afns.ualberta.ca/cgview_server/). Methylomonas sp. DH-1 was shown to be phylogenetically closely related to Methylomonas koyamae Fw12E-Y T based on 16S sequence similarity [4]. Electronic DNA-DNA hybridization (DDH) estimate between DH-1 and Fw12E-Y T (=JCM 16701 T ), calculated by the Genome-to-Genome Distance Calculator (http://ggdc.dsmz.de/distcalc2.php), was 73.9%, which suggests that DH-1 belongs to the Methylomonas koyamae species, while average nucleotide identity (ANI) between these two strains was calculated to be 97.76% using JSpecies [6]. Moreover, MUMMER whole-genome alignment [7] between DH-1 and Fw12E-Y T showed the close similarity of these two strains, aligning 298 out of 382 scaffolds of the Fw12E-Y T genome assembly (96.68% of the total length) on the DH-1 reference genome sequence ( Figure 2A). The ANI-based genome analysis clustered 17 strains in eleven species, the M. koyamae group being the largest one, accommodating strains Fw12E-Y T , DH-1, LM6, and R-49807 ( Figure 2B,C). Multiple non-type strains originally labeled as M. koyamae or M. methanica were classified into separate groups, which implies that repositioning is required for these strains. It was found that two strain types, M. methanica and M. denitrificans, appear to form a conspecific group (100.0% ANI; 92.3% DDH estimate), while 16S rRNA sequence similarity is at 98.9%. The complete genome sequences of DH-1 and LM6 could be aligned with each other collinearly without any indication of gross rearrangement of insertion/deletion both in chromosome and in plasmid (data not shown).
The genome annotation predicted 4669 protein-coding genes, 47 tRNA and 9 rRNA (Table 1). Furthermore, there were 3638 genes assigned to different function categories based on the clusters of orthologous genes (COG) designation (Table 2) [10]. The most abundant COG category was "General function prediction only" (381 CDSs), followed by "Signal transduction mechanisms" (349 CDSs). A single gigantic gene (AYM39_10365, 32.6 kb) was found to encode a hypothetical transmembrane The ANI-based genome analysis clustered 17 strains in eleven species, the M. koyamae group being the largest one, accommodating strains Fw12E-Y T , DH-1, LM6, and R-49807 ( Figure 2B,C). Multiple non-type strains originally labeled as M. koyamae or M. methanica were classified into separate groups, which implies that repositioning is required for these strains. It was found that two strain types, M. methanica and M. denitrificans, appear to form a conspecific group (100.0% ANI; 92.3% DDH estimate), while 16S rRNA sequence similarity is at 98.9%. The complete genome sequences of DH-1 and LM6 could be aligned with each other collinearly without any indication of gross rearrangement of insertion/deletion both in chromosome and in plasmid (data not shown).
A single gigantic gene (AYM39_10365, 32.6 kb) was found to encode a hypothetical transmembrane protein with repetitive domains which have Ca 2+ and carbohydrate binding property. Eight contigs from Fw12E-Y T were aligned consecutively with this sequence and complete gene sequences were found from strains LM6, R-49807, and R-45378, while no homologous sequence could be found from other genomes, which suggests that this is a common characteristic of the Methylomonas koyamae species that could increase bacterial fitness under specific environmental niches [11]. All genes required for a type I methanotrophic lifestyle were identified. One functional operon encoding particulate methane monooxygenase (pMMO, pmoCAB), and the pxm operon (pxmABC), encoding the copper membrane monooxygenase [12] was determined in the DH-1 genome. All genes for carbon fixation via the ribulose monophosphate pathway were predicted. Genes encoding PQQ-dependent methanol dehydrogenases (mxaFJGIRSACKLDEK) along with the PQQ biosynthesis gene cluster (pqqBCDE) for methanol oxidation were detected. The tetrahydrofolate (H4F)-and tetrahydromethanopterin (H4MPT)-mediated formaldehyde oxidation pathways and formate dehydrogenase were encoded. Notably, the genome of DH-1 possesses two types of the gene cluster encoding 3-hexulose-6-phosphate synthase (hps) and 6-phospho-3-hexuloisomerase (phi) including a hps-phi operon and another hpsi gene encoding an hps-phi fused protein [13].

Functional Analysis of the Complete Genome Sequence of Methylomonas sp. DH-1 and the Production of Succinate from Methane
The existence of the complete Embden-Meyerhof-Parnas (EMP) pathway, the pentose-phosphate pathways (PPPs), and the Entner-Doudoroff pathway (EDD) along with a complete TCA cycle were confirmed. Interestingly, a complete set of genes for the serine cycle was identified together with the gene encoding phosphoenolpyruvate carboxylase (ppc) which plays a key role in the serine cycle by converting phosphoenolpyruvate (PEP) to oxaloacetate with the addition of CO 2 . Some type I methanotrophs have been predicted to have a partial serine cycle due to the absence of ppc [8,14]. The existence of PEP carboxylase together with pyruvate carboxylase and acetyl-CoA carboxylase indicates that DH-1 possesses more potential for CO 2 fixation compared to other type I methanotrophs. The ability to convert PEP to oxaloacetate, a key intermediate in the TCA cycle, is also advantageous in the production of TCA-derived products such as succinic acid. Additionally, most type I methanotrophs have no remarkable accumulation of succinate in aerobic conditions because 2-oxoglutarate cannot be converted to succinyl-CoA due to poor activity of 2-oxoglutarate dehydrogenase. It forms an incomplete "horseshoe" shaped TCA cycle [15], and consequently succinate could not be accumulated in aerobic conditions. Unusually, Methylomicoribum buryatense, Type I methanotroph, accumulated a detectable amount of succinate in aerobic conditions, because it has three different genes for succinate synthesis [16]. Even though DH-1 has only a TCA cycle as a succinate generation pathway, a large amount of succinate of up to 80 mM was successfully accumulated under aerobic growth conditions, indicating that DH-1 has 2-oxoglutarate dehydrogenase activity and the ability to operate a full TCA cycle (Figure 3). encoding 3-hexulose-6-phosphate synthase (hps) and 6-phospho-3-hexuloisomerase (phi) including a hps-phi operon and another hpsi gene encoding an hps-phi fused protein [13].

Functional Analysis of the Complete Genome Sequence of Methylomonas sp. DH-1 and the Production of Succinate from Methane
The existence of the complete Embden-Meyerhof-Parnas (EMP) pathway, the pentosephosphate pathways (PPPs), and the Entner-Doudoroff pathway (EDD) along with a complete TCA cycle were confirmed. Interestingly, a complete set of genes for the serine cycle was identified together with the gene encoding phosphoenolpyruvate carboxylase (ppc) which plays a key role in the serine cycle by converting phosphoenolpyruvate (PEP) to oxaloacetate with the addition of CO2. Some type I methanotrophs have been predicted to have a partial serine cycle due to the absence of ppc [8,14]. The existence of PEP carboxylase together with pyruvate carboxylase and acetyl-CoA carboxylase indicates that DH-1 possesses more potential for CO2 fixation compared to other type I methanotrophs. The ability to convert PEP to oxaloacetate, a key intermediate in the TCA cycle, is also advantageous in the production of TCA-derived products such as succinic acid. Additionally, most type I methanotrophs have no remarkable accumulation of succinate in aerobic conditions because 2-oxoglutarate cannot be converted to succinyl-CoA due to poor activity of 2-oxoglutarate dehydrogenase. It forms an incomplete "horseshoe" shaped TCA cycle [15], and consequently succinate could not be accumulated in aerobic conditions. Unusually, Methylomicoribum buryatense, Type I methanotroph, accumulated a detectable amount of succinate in aerobic conditions, because it has three different genes for succinate synthesis [16]. Even though DH-1 has only a TCA cycle as a succinate generation pathway, a large amount of succinate of up to 80 mM was successfully accumulated under aerobic growth conditions, indicating that DH-1 has 2-oxoglutarate dehydrogenase activity and the ability to operate a full TCA cycle (Figure 3). The existence of a set of genes related to various secondary metabolite biosynthesis pathways via the methylerythritol 4-phosphate (MEP) pathway, including isoprenoid and carotenoid pathways, was confirmed ( Figure 4). Notably, the DH-1 genome contains two genes encoding 1-Deoxy-D-xylulose 5-phosphate synthase catalyzing the first step of the MEP pathway. However, the carotenoids biosynthesis pathway in DH-1 has not been fully discovered. From the genome mining analysis, squalene/phytoene synthase (sqs) which is committed in the carotenoid synthesis pathway and the gene cluster related to 4,4′-diapolycopene biosynthesis including diapolycopene oxygenase (crtP), phytoene desaturase (crtI) and aldehyde dehydrogenase (ald) were identified in DH-1. Other Figure 3. The growth curve and succinic acid production of Methylomonas sp. DH-1. The cells were grown in nitrate mineral salts (NMS) media with a bioreactor where a gas mixture of 30% CH 4 , 55% N 2 , 15% O 2 , was continuously fed at the speed of 40 mL min −1 . The methane concentration in the reactor off-gas was indicated in the off-set y-axis.
The existence of a set of genes related to various secondary metabolite biosynthesis pathways via the methylerythritol 4-phosphate (MEP) pathway, including isoprenoid and carotenoid pathways, was confirmed ( Figure 4). Notably, the DH-1 genome contains two genes encoding 1-Deoxy-D-xylulose 5-phosphate synthase catalyzing the first step of the MEP pathway. However, the carotenoids biosynthesis pathway in DH-1 has not been fully discovered. From the genome mining analysis, squalene/phytoene synthase (sqs) which is committed in the carotenoid synthesis pathway and the gene cluster related to 4,4 -diapolycopene biosynthesis including diapolycopene oxygenase (crtP), phytoene desaturase (crtI) and aldehyde dehydrogenase (ald) were identified in DH-1. Other potential secondary metabolites that can be synthesized by Methylomonas sp. DH-1 were identified using antiSMASH [17]. The results indicated that eight possible gene clusters encoding secondary metabolites were identified in Methylomonas sp. DH-1 including aryl polyene, bacteriocins, terpene, hserlactone, and T1pks-Nrps. Among them, aryl polyene can protect the bacterium from reactive oxygen species, similar to the functionality of carotenoids [18]. potential secondary metabolites that can be synthesized by Methylomonas sp. DH-1 were identified using antiSMASH [17]. The results indicated that eight possible gene clusters encoding secondary metabolites were identified in Methylomonas sp. DH-1 including aryl polyene, bacteriocins, terpene, hserlactone, and T1pks-Nrps. Among them, aryl polyene can protect the bacterium from reactive oxygen species, similar to the functionality of carotenoids [18].

Nucleotide Sequence Accession Number
The completed genome sequence of Methylomonas sp. DH-1 was deposited at GenBank under accession number CP014360 and CP014361. In addition, the strain was deposited at the Korean Collection for Type Culture under the KCTC number 13004BP.

Conclusions
We sequenced and analyzed the whole genome of a newly isolated type I methanotroph, Methylomonas sp. DH-1 consisting of a 4.86 Mb chromosome and a 278 kb plasmid. Methylomonas sp. DH-1 accumulated a large amount of succinate (up to 80 mM) under aerobic conditions most probably due to 2-oxoglutarate dehydrogenase activity, showing its biocatalytic potential for methane bioconversion. The existence of PEP carboxylase, pyruvate carboxylase and acetyl-CoA carboxylase can enable the DH-1 strain to fix CO2 more efficiently compared to other type I methanotrophs. A set of genes related to various secondary metabolite biosynthesis pathways via the MEP pathway was also identified. The availability of a complete genome sequence of Methylomonas sp. DH-1 contributes to a system-level understanding of methanotrophic metabolism which provides valuable resources for metabolic engineering of this strain for overproduction of value-added chemicals from methane.

Nucleotide Sequence Accession Number
The completed genome sequence of Methylomonas sp. DH-1 was deposited at GenBank under accession number CP014360 and CP014361. In addition, the strain was deposited at the Korean Collection for Type Culture under the KCTC number 13004BP.

Conclusions
We sequenced and analyzed the whole genome of a newly isolated type I methanotroph, Methylomonas sp. DH-1 consisting of a 4.86 Mb chromosome and a 278 kb plasmid. Methylomonas sp. DH-1 accumulated a large amount of succinate (up to 80 mM) under aerobic conditions most probably due to 2-oxoglutarate dehydrogenase activity, showing its biocatalytic potential for methane bioconversion. The existence of PEP carboxylase, pyruvate carboxylase and acetyl-CoA carboxylase can enable the DH-1 strain to fix CO 2 more efficiently compared to other type I methanotrophs. A set of genes related to various secondary metabolite biosynthesis pathways via the MEP pathway was also identified. The availability of a complete genome sequence of Methylomonas sp. DH-1 contributes to a system-level understanding of methanotrophic metabolism which provides valuable resources for metabolic engineering of this strain for overproduction of value-added chemicals from methane.

Bacterial Growth, DNA Isolation, Genome Assembly and Annotation
Methylomonas sp. DH-1 was isolated from the activated sludge of a brewery plant based in a nitrate mineral salts (NMS) medium using enrichment culture with methane as a sole carbon source as described by Hur et al. [4]. Liquid pre-cultures were grown in a 180 mL baffled-flask with a 10 mL NMS medium containing 10 µM CuSO 4 with a supplement of 30% methane (v/v) as a sole carbon source at 30 • C and 230 rpm, sealed with a screw cap. The pre-cultures were then inoculated into 50 mL of fresh medium in a 500 mL baffled-flask for large-scale cultivation.
The genomic DNA was extracted using a Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA). The library construction and sequencing were carried out at NICEM, Seoul, Republic of Korea. The general genome properties were first obtained using Illumina HiSeq 2500 platform-based draft genome sequencing (2 × 151 nt; 1.6 Gb) and then, the PacBio RS II platform was used to obtain the complete genome sequence, which was polished and verified using the previously generated Illumina reads. The 807.6 Mb filtered polymerase reads, produced from the PacBio RS II sequencing using P6-C4 chemistry with 119-fold average coverage was assembled into two contigs using the hierarchical genome assembly process RS_HGAP.3 [19]. The identification of sequence overlap at both ends and the alignments with the Illumina assemblies revealed their circular structures. The genome annotation was performed by integrating results from Prokaryotic Genome Annotation Pipeline (PGAP) (http://www.ncbi.nlm.nih.gov/genome/annotation_prok/), Integrated Microbial Genomes (IMGs) (http://jgi.doe.gov/data-and-tools/img/), Rapid Annotation using Subsystem Technology (RAST) (http://rast.nmpdr.org/) and PROKKA [20] on the basis of stop codons. The priority for choosing functional annotation was in the order of PGAP, IMG, PROKKA, and RAST (Tables S1 and S2). The genome sequence information for comparative and phylogenomic analyses was downloaded from the RefSeq database. The ANI-based genome comparison and clustering were done using DREP [21]. The universal prokaryotic marker gene sequences, identified using the PHYLOSIFT [22], were concatenated into one and an approximately maximum-likelihood tree was constructed using the FASTTREE 2 [23].

Analytical Methods
The supernatant of cultures was separated by centrifugation. The succinate was quantified using a HPLC equipped with an Aminex HPX-87 column (Bio-Rad, Hercules, CA, USA) and a refractive index detector. Sulfuric acid to the amount of 0.005 M was used as the mobile phase with a flow rate of 0.7 mL/min at 60 • C. The bioreactor off-gas was connected to the GC (Agilent 7890A, Santa Clara, CA, USA) and the methane composition was analyzed by a GC equipped with a TCD detector.