Scholarly article on topic 'Reversing methanogenesis to capture methane for liquid biofuel precursors'

Reversing methanogenesis to capture methane for liquid biofuel precursors Academic research paper on "Biological sciences"

0
0
Share paper
Academic journal
Microb Cell Fact
OECD Field of science
Keywords
{""}

Academic research paper on topic "Reversing methanogenesis to capture methane for liquid biofuel precursors"

Soo et al. Microb Cell Fact (2016) 15:11 DOI 10.1186/s12934-015-0397-z

Microbial Cell Factories

Reversing methanogenesis to capture methane for liquid biofuel precursors

Valerie W. C. Soo11, Michael J. McAnulty11, Arti Tripathi1, Fayin Zhu1, Limin Zhang2,6, Emmanuel Hatzakis3, Philip B. Smith5, Saumya Agrawal7, Hadi Nazem-Bokaee1, Saratram Gopalakrishnan1, Howard M. Salis1, James G. Ferry4, Costas D. Maranas1, Andrew D. Patterson2 and Thomas K. Wood1,4* ©

Abstract

Background: Energy from remote methane reserves is transformative; however, unintended release of this potent greenhouse gas makes it imperative to convert methane efficiently into more readily transported biofuels. No pure microbial culture that grows on methane anaerobically has been isolated, despite that methane capture through anaerobic processes is more efficient than aerobic ones.

Results: Here we engineered the archaeal methanogen Methanosarcina acetivorans to grow anaerobically on methane as a pure culture and to convert methane into the biofuel precursor acetate. To capture methane, we cloned the enzyme methyl-coenzyme M reductase (Mcr) from an unculturable organism, anaerobic methanotrophic archaeal population 1 (ANME-1) from a Black Sea mat, into M. acetivorans to effectively run methanogenesis in reverse. Starting with low-density inocula, M. acetivorans cells producing ANME-1 Mcr consumed up to 9 ± 1 % of methane (corresponding to 109 ± 12 |mol of methane) after 6 weeks of anaerobic growth on methane and utilized 10 mM FeCl 3 as an electron acceptor. Accordingly, increases in cell density and total protein were observed as cells grew on methane in a biofilm on solid FeCl3. When incubated on methane for 5 days, high-densities of ANME-1 Mcr-producing M. acetivorans cells consumed 15 ± 2 % methane (corresponding to 143 ± 16 |mol of methane), and produced 10.3 ± 0.8 mM acetate (corresponding to 52 ± 4 |mol of acetate). We further confirmed the growth on methane and acetate production using 13C isotopic labeling of methane and bicarbonate coupled with nuclear magnetic resonance and gas chromatography/mass spectroscopy, as well as RNA sequencing.

Conclusions: We anticipate that our metabolically-engineered strain will provide insights into how methane is cycled in the environment by Archaea as well as will possibly be utilized to convert remote sources of methane into more easily transported biofuels via acetate.

Keywords: Reverse methanogenesis, Anaerobic oxidation of methane, Methyl-coenzyme M reductase

Background

Advances in horizontal drilling since the late 1970s have made methane resources economical through hydraulic fracturing [1]. All continents, except Antarctica, are known to possess recoverable shale reserves, with a global amount of shale gas totaling up to 7300 trillion cubic feet [2]. The major component of shale gas is

'Correspondence: tuw14@psu.edu

+Valerie W. C. Soo and Michael J. McAnulty contributed equally

1 Department of Chemical Engineering, The Pennsylvania State University,

University Park, PA 16802-4400, USA

Full list of author information is available at the end of the article

methane, a potent greenhouse gas that is 25 times more damaging than carbon dioxide (CO2) over a span of 100 years [3]. To exploit the abundance of these energy sources and to reduce methane emissions, state-of-the-art chemical plants employ Fischer-Tropsch processes to convert methane into liquid fuels (e.g., Shell Pearl Facility in Qatar with an output of 260,000 barrels of oil per day [4]); however, these plants rely on complex technology that demands large-scale investment [5]. An alternative approach is to biologically convert methane to liquid fuels, a method that may be more suitable for the thousands of remote fracking sites. Biological conversion of

O© 2016 Soo et al. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// Centfcll creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/ zero/1.0/) applies to the data made available in this article, unless otherwise stated.

methane is a more economically and environmentally sustainable technology, since it requires a smaller footprint and is less technologically complex [6]. Hence, we explored methane activation using an anaerobe as our biological catalyst and focused on anaerobic processes, as they confer higher energy and carbon yield efficiencies with lower CO2 emissions than aerobic ones for converting methane into fuels [6]. Furthermore, creating biofuels from methane should have a high impact for second and third generation biofuels since they will compete on yield and feedstock cost rather than on fermentation time.

Anaerobic oxidation of methane (AOM) is a key regulator of global fluxes of methane and the carbon cycle, capturing up to 300 Tg of methane per year [7]. However, AOM is a little-understood biological process driven by natural consortia consisting of ANME and syntrophic bacteria. Despite decades of effort, these natural consortia have never been successfully isolated, most probably due to their long lag phase (~60 years) [8] and doubling time (~7 months) [9]. Further, separate cultivation of ANME and their syntrophic bacteria has not been reported.

The first step to anaerobically activate methane in a process known as reverse methanogenesis is likely to be catalyzed by Mcr, as suggested by the prevalence of mcr genes in ANME populations [10]. Trace AOM has also been observed in anaerobic methanogens such as Metha-nothermobacter marburgensis [11] and M. acetivorans [12, 13]; however, none of these strains utilize methane as the major energy and carbon source for growth.

Here we engineered an anaerobic archaeal strain that uses methane and bicarbonate as the major energy and carbon sources to support its growth. Our ultimate goal is to produce liquid fuels efficiently via anaerobic processes, and for this reason, we set out to create a system that could capture methane using simple technology on a laboratory scale. We cloned the genes encoding Mcr from the metagenome of an unculturable ANME-1 population originally identified in Black Sea mats [14, 15], where AOM rates are the highest in all studied aquatic environments [7]. M. acetivorans was chosen as the host, as this strain is metabolically diverse since it has the largest archaeal genome [16], is genetically tractable [17], and encodes a native Mcr for producing methane during methanogenesis [CH3-SCoM + HSCoB ^ CH4 + CoBS-SCoM, where CH3-SCoM is methyl-coenzyme M (2-(methylthio)ethanesulfonate), HSCoB is coenzyme B (7-mercaptoheptanoylthreonine phosphate), and CoBS-SCoM is the heterodisulfide of coenzyme B and coenzyme M]. Therefore, we reasoned this host has the necessary cofactors and cellular machinery to produce active ANME-1 Mcr.

Results and discussion

Growth on methane via ANME-1 Mcr

In contrast to previous studies where trace AOM was detected in M. acetivorans [12, 13], we forced the engineered cells, M. acetivorans with cloned ANME-1 mcr genes from the Black Sea metagenome, to utilize methane for growth by not providing an additional carbon source (e.g., methanol, acetate, trimethylamine). Instead, we added various electron acceptors (Fe3+, NO3-, NO2-, SO42-, and Mn4+) that could remove the electrons generated by growth on methane and to mimic the syntrophic role of sulfate-reducing bacteria associated with ANME-1 organisms. We then characterized the ability of the engineered strain to catabolize methane using growth-based measurements, 13C isotopic labeling of methane and bicarbonate, and RNA sequencing.

The ANME-1 mcrBGA genes encoding for ANME-1 Mcr were cloned into M. acetivorans using vector pES1 [18] to form pES1-MATmcr3 which expresses mcrBGA using its native promoter Pmcr ANME-1. This native promoter was found to be superior to that of the CO dehy-drogenase/acetyl-CoA synthase promoter (Pcdh) from Methanosarcina thermophila [18] and the mcr promoter from M. acetivorans (Pmcr_M acetivorans) (Additional file 1: Figure S1). To compare the promoters, we conducted short-duration (5 days) experiments that utilized high densities of cells grown first on methanol (109 CFU/mL), which were subsequently added to minimal high-salt (HS) medium (Additional file 1: Table S1) with methane and bicarbonate and tested for the amount of methane utilization. In cells exposed to methane, ANME-1 Mcr was produced in M. acetivorans by Pcdh and Pmcr_ANME-!, as shown by the appearance of ANME-1 McrA-FLAG (the a subunit of the ANME-1 Mcr complex that was tagged with a FLAG epitope) in a Western blot (Additional file 1: Figure S2). Corroborating the production of ANME-1 Mcr by Pmcr_ANME-1, methane consumption using whole cells was the highest when Pmcr_ANME-1 was used to express ANME-1 mcr (15 ± 1.7 %, corresponding to 143 ± 16 |imol of methane, Fig. 1a). Significant differences between ANME-1 Mcr-producing cultures and those harboring an empty plasmid were observed throughout 5 days (Fig. 1b); therefore, the native ANME-1 mcr promoter was used for the remainder of the experiments. Note that absolute amounts of reacting species are indicated in |mol to provide a more direct comparison between liquid products and gaseous substrates since the volume of the headspace and the volume of the liquid in each culture tube of different experiments were not the same.

Of the tested electron acceptors, the best growth on methane was seen with 10 mM Fe3+. After 6 weeks of

Empty Mcr3 Mcr2 Merl No cell

Fig. 1 Methane consumption by ANME-1 Mcr-producing M. acetivorans during high-density growth experiments. All values are presented as mean ± sem of three cultures. **P < 0.05, *P < 0.15. a Methane consumption of cells harboring pES1 -MATmcr1 (Mcr1), pES1 -MATmcr2 (Mcr2), pES1-MATmcr3 (Mcr3), and pES1(Pmat) (Empty) after 5 days showing that pES1-MATmcr3 is the best construct in performing reverse methanogenesis in M. acetivorans.The absence of cells (No cells) in the same growth medium (HS-methane medium and 0.1 mM FeCl3 and puromycin) showing negligible methane reduction confirmed minimal gas leakage from the tubes. b Increase in methane consumption by high cell-density cultures of M. acetivorans/pES1 -MATmcr3 (Mcr3) in comparison to pES1(Pmat)-harboring cells (Empty) for 5 days showing the importance of the cloned mcr. c Higher production of extracellular acetate was observed for high cell-density cultures of M. acetivorans/pES1-MATmcr3 (Mcr3) than for M. ace-tivorans/pES1(Pmat) (Empty), M. acetivorans/pES1 -MATmcr2 (Mcr2), and M. acetivorans/pES1 -MATmcr1 (Mcr1) after 5 days on methane and 0.1 mM FeCl3. The absence of cells (Nocells) in the same growth medium (HS-methane medium with 0.1 mM FeCl3 and puromycin) showing negligible extracellular acetate confirmed minimal abiotic formation of acetate

incubation using low-density inocula, growth of M. ace-tivorans/pES1-MATmcr3 was observed in conjunction with 9 ± l % of methane consumption (corresponding to 109 ± 12 |imol of methane). Critically, cell growth was evident from steadily increasing total protein levels (Fig. 2a) and from the increase in the number of cells associated with the iron precipitates (Figs. 2b, 3a). Note that the cells were differentiated from solid precipitates via SYTO9 staining, and the morphology of cells grown on methane and 0.1 mM Fe3+ (Fig. 3b) is comparable with the unengineered strain M. acetivorans grown on methanol (Fig. 3c). Up to 97 % of cells were closely associated with salt precipitates (instead of being in suspension), which suggests that these cells are growing as a

biofilm on the precipitates in order to transfer electrons to the oxidized form of iron (Fig. 3a). We observed -fivefold increase in M. acetivorans/pESl-MATmcr3 cells on precipitates over 6 weeks of incubation (Figs. 2b, 3a); however, there was only a small two-fold increase in planktonic cell numbers. No significant changes were observed for total protein concentrations and cell counts for either cells harboring an empty plasmid or cells with pES1-MATmcr3 incubated with nitrogen (changes for cell numbers on precipitates when under a nitrogen headspace were negligible after 70 days of incubation).

The identity of M. acetivorans during growth on methane was confirmed by 16S ribosomal DNA (rDNA) sequencing (via PCR amplification using archaeal 16S

a 14-1

14 21 28

Time (day)

14 21 28 35

Time (day)

Fig. 2 Long-term growth of ANME-1 Mcr-producing M. acetivorans on methane and 10 mM FeCl3 during low-density growth experiments. a Total protein concentration of a representative culture of M. acetivorans/pES1-MATmcr3 grown on methane and 10 mM FeCl3 after 14, 21, and 28 days of incubation. Values are presented as means ± standard deviations. b Number of M. acetivorans/pES1-MATmcr3 (Mcr3) and M. acet/vorans/pES1(Pmat) (Empty) cells per area (mm2) of precipitates after 14, 21, 28, and 35 days of incubation on methane and 10 mM FeCl3. Values are presented as mean ± sem of three or more cultures

reaction for a natural consortia from marine methane-seep sediment in the Eel River Basin in California [19] presented as CH4 + 8 Fe(OH)3 + 15 H+ ^ HCO3- + 8 Fe2+ + 21 H2O, which is thermodynamically favorable (AG^n = -270 kJ/mol).

Confirmation of growth on methane using 13CH4

Using a high inoculum size, acetate was detected extracellularly in methane-grown cultures expressing ANME-1 mcr under Pmcr ANME-1 after 5 days. Acetate concentrations were the highest for cultures producing ANME-1 Mcr from Pmcr_ANME-i [10.3 ± 0.8 mM corresponding to 52 ± 4 |imol, compared to 6.1 ± 0.8 mM produced by cultures with empty vector pES1(Pmat)] among all tested strains (Fig. 1c). The use of higher inoculum size also revealed low amounts of extracellular formate and pyruvate in cells with Pmcr ANME-1; however, these amounts were not significantly different from the other tested strains. Therefore, we concluded that acetate was probably the main product of anaerobic activation of methane in our engineered strain.

To demonstrate absolutely that the methane was utilized for growth and acetate formation, 13CH4 was used as the main carbon source and the acetate product was identified three ways, by one-dimensional and multidimensional nuclear magnetic resonance (NMR) and by gas chromatography coupled with mass spectros-copy (GC/MS). 13CH4 incorporation was confirmed by a threefold higher 13C/12C ratio of acetate in ANME-1 Mcr-producing cells (a 13C/12C ratio of 0.33) in comparison to cells without ANME-1 Mcr (a 13C/12C ratio of 0.11) (Fig. 4a, b), and the structure of acetate confirmed by NMR (Fig. 4c, d). Our results are consistent with previous results based on trace utilization of methane [13], in which 13C-labeled methane is converted to acetate in M. acetivorans.

gene-specific primers). No PCR product was obtained when bacterial 16S gene-specific primers were used, providing further evidence for homogeneity of the culture. Furthermore, genome sequencing of the engineered strain that grew on methane confirmed the strain. Also, the plasmid pES1-MATmcr1 was not lost from the population after 30 days of growth on methane, as we could PCR-amplify ANME-1 mcrA from the plasmid of the methane-grown culture (Additional file 1: Figure S3).

These results indicate that we have metabolically engineered M. acetivorans to grow in the presence of methane and an electron acceptor, and that the ANME-1 Mcr is necessary for this growth. Growth on methane with ferric iron is reasonable, with one example

Confirmation of growth on methane using 13C-labeled bicarbonate

We found that 19 ± 9 |mol of bicarbonate from the growth medium (Additional file 1: Table S1, 5 mL) was consumed in high-density cultures with ANME-1 Mcr produced from Pmcr ANME-1 showing a net methane consumption of 81 ± 16 |mol (Fig. 1a) after normalization with cells harboring an empty plasmid. Therefore, we hypothesized that bicarbonate is utilized along with methane to form acetate as has been seen for trace utilization of methane [13]. Corroborating this hypothesis, when high-density cultures of M. acetivorans/pES1-MATmcr3 were incubated with 13C-labeled bicarbonate and 10 mM FeCl3 (with methane in the headspace), we found 13C incorporation into acetate using NMR and

3 Phase contrast Epifluorescence Overlay Time

Fig. 3 Biofilm growth of ANME-1 Mcr-producing M. acetivorans on methane. a Representative hemocytometer views (x400 magnification) of M. acetivorans/pES1-MATmcr3 cells growing as biofilms on the precipitates after incubation on methane and 10 mM FeCl3. All cultures were stained with SYTO9 dye for microscopic examination. Under epifluorescence settings, cells appear green but abiotic precipitates do not fluoresce. Scale bars 20 |jm. b Electron micrograph of irregular cocci M. acetivorans/pES1-MATmcr1 growing on methane and 0.1 mM FeCl3. Scale bar 0.2 |jm. c Electron micrograph of irregular-cocci M. acetivorans grown on methanol. Scale bar 0.5 jm. Red arrows show the mono-layered cell wall

Fig. 4 Extracellular production of acetate from 13CH4 in ANME-1 Mcr-producing M. acetivorans via GC/MS and NMR. a and b show GC/MS spectra of culture supernatants used to identify acetate from 13CH4. Unlabeled acetate (12C Acetate, black arrows) has a molecular mass of 60, whereas 13C-labeled acetate (13C Acetate, red arrows) has a molecular mass of 61. a M. acetivorans harboring empty plasmid pES1(Pmat) grown on 13C-labeled methane (13CH4) for 10 days. b M. acetivorans harboring pES1-MATmcr3 grown on 13CH4 for 10 days. Note that the y-axes differ for a and b. c NMR HSQC spectrum identifying acetate by showing the one bond correlation between the methyl protons and the methyl carbon of acetate (13CH3-C). d NMR HMBC spectrum identifying acetate in an independent manner by showing the two bonds correlation between the methyl protons and the carboxylic carbon of acetate (13CH3-COO). The x-axes denote 1H NMR plots, whereas the y-axes denote 13C NMR plots

GC/MS (Additional file 1: Figure S4), which was consistent with our expectation.

Metabolic flux analysis (MFA) using 13C-labeled bicarbonate as a tracer was performed to quantify metabolite flows through the methanogenesis pathways resulting in the observed labeling distribution of acetate (Additional file 1: Figure S5). MFA revealed that the higher abundance of the m/z 60 ion of acetate was due to the presence of a largely unlabeled CO2 pool arising from the complete oxidation of unlabeled methane to CO2 via the

methylotrophic pathway. The reduction of CO2 by Cdh incorporates the labeled bicarbonate into the carbonyl carbon of acetate, which is reflected in the m/z 61 ion. The abundance of m/z 62 ion of acetate can only arise from the incorporation of labeled bicarbonate into the methyl group of acetate (the 62 ion peak is small but present at higher levels than the background noise). MFA revealed a non-zero reverse flux through the methylotrophic pathway, thereby allowing the reduction of 13C-bicarbonate to the methyl group of acetate. In addition to this, MFA

estimated the in vivo co-utilization ratio of bicarbonate/ methane to be 0.06 mol/mol based on the extent of labeling of acetate; hence, complete oxidation of methane to CO2 via the methylotrophic pathway (v3, Additional file 1: Figure S5) produces an unlabeled CO2 pool resulting in a large unlabeled fraction of acetate, as seen in the 13C-labeled bicarbonate experiment for growth on methane (Additional file 1: Figure S4).

Although cysteine could be another carbon source (initial concentration of 3.2 mM or an amount of 16 ^mol, Additional file 1: Table S1), cysteine consumption in all tested strains ranged from 32 to 49 %, with no significant differences between ANME-1 Mcr-producing strains and the strain harboring an empty plasmid. Methyl sulfides were proposed as products of AOM [20], but there is not enough sulfur in the system for methyl sulfides to be generated and thermodynamically account for significant amounts of methane consumption. We do not consider trace vitamins in the growth medium as major carbon sources, as their concentration is insignificant (less than 0.005 mM, Additional file 1: Table S1).

Confirmation of growth on methane via iron reduction

Growth seen at 10 mM Fe3+ (Figs. 2, 3a) suggests the following thermodynamically possible (AGJjn = -325 kJ/ mol) stoichiometry with Fe3+ as the terminal electron acceptor shown in Eq. (1).

4 CH4 + 2 HCO- + 8 Fe:

^ 3 CH3COO" + 8 Fe2+ + 9 H+

Our results using 13C-labeled methane and bicarbonate (Fig. 4; Additional file 1: Figure S4) demonstrate that methane and bicarbonate are converted to acetate as indicated in equation [1]; therefore, we investigated whether iron is reduced as well in long-term growth experiments. In M. acetivorans/pES1-MATmcr3 cultures grown on methane and 10 mM FeCl3, 19 ± 1 ^mol of Fe3+ were reduced into Fe2+ after 35 days of incubation. No other terminal electron acceptors exist in the liquid phase in significant quantities for the observed methane consumption. Hence, growth on methane and bicarbonate proceeds via iron reduction, which is consistent with cells growing as a biofilm on the FeCl3 precipitates (Figs. 2b, 3a). 3

Genome-scale metabolic modeling

In further support of the overall reaction shown in equation [1], a recently updated genome-scale reconstruction of M. acetivorans (manuscript in preparation) indicates a thermodynamically feasible phenotype of co-metabolizing methane and bicarbonate using Fe3+ as the terminal electron acceptor. One of the possible metabolic

outcomes describing the production of only acetate follows the overall stoichiometry shown in equation [1]. Alternative products supported by the metabolic model include pyruvate and formate. The thermodynamic feasibility constraints implied by the rate of production of products (acetate, pyruvate, and formate) and the rate of Fe3+ oxidation is described in Fig. 5. For the genome-scale reconstruction analysis, we assumed that there is no cellular growth and ATP production is only used for maintenance requirements. We found that acetate and pyruvate compete with Fe3+ for electrons by reductive carboxylation of methanogenesis intermediates. Pyruvate also competes with acetate production, by lowering ATP production through substrate-level phosphorylation. Therefore, the engineered cells grow on methane and bicarbonate as indicated in equation [1]. It is interesting to note that the thermodynamic infeasibility of direct reversal of the aceticlastic pathway is overcome by coupling acetate production with oxidation of methane to CO2 via the methylotrophic pathway.

The likely biochemical pathway for the overall reaction shown in equation [1] and the basis for the genome-scale model are shown in Fig. 6. The oxidation of CH4 yielding CH3-SCoM is catalyzed by ANME-1 Mcr. The methyl group is transferred to H4SPT (tetrahydrosarcinapterin) by reversal of the reaction catalyzed by cytoplasmic methyltransferase (Cmt) [21] or membrane-bound meth-yltransferase (Mtr) [22] of methane-producing pathways. In the proposed pathway, Mtr is driven by a sodium gradient supported by oxidation of CH3-H4SPT (methyl-H4SPT) to CO2 and reduction of Fe3+ that generates a proton gradient exchanged for sodium (not shown). The CH3COSCoA (acetyl-CoA) is a product of the reversible Cdh complex [23, 24]. The carbonyl group of CH3COS-CoA is derived from the ferredoxin-dependent [25] reduction of exogenous and endogenous bicarbonate derived from oxidation of the methyl of CH3-H4SPT consistent with the 13C labeling pattern (red colored carbon atoms). Two molecules of reduced ferredoxin derive from oxidation of CHO-MF (formyl methano-furan) by formyl methanofuran dehydrogenase [26]. The remaining reduced ferredoxins derive from oxidation of HSCoB and HSCoM (coenzyme M) catalyzed by a hypothesized electron transport system bypassing reversal of the membrane-bound system in the forward pathway dependent on a sodium gradient [27]. This system is hypothesized to involve an electron bifurcating complex (EBC), which transfers electrons from the oxidation of HSCoM and HSCoB to a membrane-bound electron carrier reducing Fe3+. The EBC is postulated to couple the thermodynamically unfavorable reduction of ferre-doxin (Em = -420 mV) to the favorable reduction of Fe3+ (Em = +770 mV). Importantly, the system is postulated

Ferric Uptake (mol/mol methane)

Fig. 5 Thermodynamically feasible envelopes for the production of a acetate, b pyruvate, and c formate as a function of consumed Fe3+ as the external electron acceptor. Horizontal lines denote the measured amount of acetate and formate, respectively. Both lines intersect with the thermodynamically feasible envelope suggested by the model. Acetyl-CoA competes with iron for electrons, thus resulting in a reduction of maximum acetate production with increasing Fe3+ consumption. At lower Fe3+ consumption, the electron transport chain is not sufficient to meet the required ATP demands, thus making acetate production essential. Increasing Fe3+ uptake increases the flux through the methylotrophic, and consequently, the electron transport chain. At a Fe3+ uptake of 0.28 mol/mol methane, thermodynamic limitations restrict pyruvate production, thereby resulting in a reduced sensitivity of minimum acetate to Fe3+ uptake. A Fe3+ uptake of 0.6 mol/mol methane is sufficient to meet the ATP demands using the electron transport chain alone, thus decreasing the minimum required acetate production to zero. Since CO2 was not observed as a secreted product, the end product of the methylotrophic pathway was hypothesized to be formate produced by abiotic deformylation of formylmethano-furan, which competes with acetate and pyruvate for carbons causing the maximum pyruvate production to decline when Fe3+ uptake surpasses 0.6 mol/mol methane. However, the minimum formate production still remains zero due to the fact that formate can be oxidized to produce CO2 for acetate and pyruvate production. Beyond a Fe3+ uptake of 2.3 mol/mol methane, the thermodynamics of the system do not favor the production of acetate, thus making formate and pyruvate production mandatory

to drive the unfavorable oxidation of CH4 and transfer of the methyl group to H4SPT. The CH3COSCoA is converted to acetate by the reversible phosphotransacetylase and acetate kinase that produces ATP, supporting growth with CH4 as the energy source.

We also used RNA sequencing to determine differential gene expression in methane-grown cells producing

ANME-1 Mcr by Pmcr ANME-1 in comparison to the same cells grown on methanol. Two genes were induced by growth on methane (Additional file 1: Table S2): MA1997 (hypothetical protein; absent in methanol-grown cells) and MA0463 (ferredoxin; tenfold induction via methane). Critically, 27 genes were repressed by growth on methane, including genes specific for methanogenesis

Fig. 6 Proposed reversal of the aceticlastic pathway. The pathway shown imposes Fe3+ at a concentration producing the maximum amount of acetate (see Fig. 5). Red coded atoms indicate 13C-labeling. EBC electron bifurcating complex, Fdr reduced ferredoxin, Fdo oxidized ferredoxin, CH3-SCoM methyl-coenzyme M, HSCoM coenzyme M, HSCoB coenzyme B, CoMS-SCoB heterodisulfide of coenzyme M and coenzyme B, H4SPT5,6,7,8-tet-rahydrosarcinapterin, CH3-H4SPTmethyl-H4SPT, HSCoA coenzyme A, CHO-MFformyl methanofuran, CH3COSCoA acetyl-CoA, HCO- bicarbonate ion, HPO4~ hydrogen phosphate ion, CH3COPO4- acetyl phosphate, ATP adenosine triphosphate, ADP adenosine diphosphate, and CH3COOH acetic acid

(6- to 19-fold repression) such as mcrA, mcrB, mtaB 1, and mtaB 2 (Additional file 1: Figure S6). Altogether, since the methane-generating pathway from carbon dioxide was repressed (e.g., host mcrAB, mtaB 1B2, cdhB, fmdD) (Additional file 1: Figure S6), growth on methane appears to proceed through reverse methanogenesis.

Conclusions

By producing Mcr from ANME-1 in a methanogenic host, we are able to reverse the rest of the methanogenic pathway. This work thus provides evidence indicating the importance of AOM in the methane cycle as AOM converts methane into acetate as one byproduct. Results from these laboratory-scale experiments show that reversing methanogenesis via ANME-1 Mcr production

is a promising anaerobic technology to oxidize methane for biofuel production as well. Depending on the type of ANME populations and their syntrophic partners, natural AOM can proceed via reduction of iron, manganese, or sulfate [7, 19]. However, with our engineered M. ace-tivorans/pES1-MATmcr3 cells, growth on methane was possible only via iron reduction, which highlights the distinctiveness of our engineered population of M. ace-tivorans/pES1-MATmcr3 cells.

M. acetivorans is an excellent host for biofuel production due to its inherent tolerance to alcohol (i.e., growth on 150 mM methanol). Furthermore, acetate from our engineered strain can be used to produce ethanol and other types of biofuels; for example, fatty alcohols, alkanes, and acyl esters through acetyl-CoA and farnesane, farnesanol,

and isopentanol through the isoprenoid pathway. The cul-turability of our engineered strain also renders it amenable to further engineering efforts to biologically capture methane and potentially convert it to usable resources. Overall, our engineered strain opens up possibilities for efficient production of acetate that can be scaled up for industrial uses and provides a host that allows active Mcr from ANME-1 to be studied for the first time.

Methods

Growth conditions of M. acetivorans

All M. acetivorans strains (Additional file 1: Table S3) were routinely grown anaerobically as pre-cultures at 37 °C in an 80 % N2/19 % CO2/1 % H2 atmosphere with mild shaking in 10 mL HS medium [28] with 150 mM methanol as the carbon source. Cell growth was measured spectrophotometrically and direct cell counts were confirmed by staining cell cultures (often containing precipitates) with SYTO9 dye (Life Technologies, Carlsbad, CA, USA) and viewed microscopically using a bright-line hemocytometer (Hausser Scientific, Horsham, PA, USA) under phase-contrast and epifluorescence settings (Zeiss Axio Scope.A1, Germany). All 28-mL culture tubes (18 x 150 mm, Bellco Glass, Vinelanad, NJ, USA) were sealed by aluminum crimp seals. Plasmids were maintained with 2 |g/mL puromycin.

For long-term (ca., 40 days) growth on methane with low-density inocula (1 %), the M. acetivorans strains were grown in 5 mL HS medium in 28-mL culture tubes with additional electron acceptors at 37 °C with mild shaking; the electron acceptors tested were FeCl3, FeSO4, KNO3, NaNO3, NaNO2, and ZnSO4 (0.05-100 mM) as well as Fe2O3 (1 mM) and MnO2 (2 mM). As a frame of reference, a concentration of 0.1 mM would be 0.5 |mol in a 5-mL culture. The preferred carbon sources of M. ace-tivorans, methanol, trimethylamine, and acetate, were not present in the medium. The headspace of the tube was filled with methane (99.999 % purity, catalog no. ME5.0RS, Praxair, Danbury, CT), and crimped with aluminum seals. PCR was used to verify the presence of pES1-MATmcr1 in methane-grown cells after 30 days of incubation. Cells were used as genomic templates for PCR amplification of ANME-1 mcrA using primers B4-f and pES1-r (Additional file 1: Table S4). Cell morphology was examined on a transmission electron microscope (FEI Tecnai G2 Spirit BioTwin, Hillsboro, OR, USA) using uranyl acetate-stained cells.

To document cell growth and methane consumption as a function of time, 107 CFU/mL of M. acetivorans/ pES1-MATmcr3 and M. acetivorans/pES 1(Pmat) cells were incubated at 37 °C for 6 weeks in 8 mL HS medium and 10 mM FeCl3 in 40-mL bottles stoppered with butyl rubber stoppers and crimped with aluminum seals. Each

culture and its methane in the headspace were sampled every 1-2 weeks. After each sampling, petroleum jelly was applied to the surface of the stoppers to prevent leaking from needle punctures. All culture bottles were inverted during incubation to prevent methane from escaping from the vessels, and were shaken to ensure homogenous mixing of liquid, precipitates, and cells. As a control, cultures grown on methanol in parallel reached a saturated density of (4 ± 0.3) x 108 CFU/mL.

For short-duration (ca., 5 days) growth on methane in which high cell-density inocula were used, 2 mL of each strain was pre-grown in 200 mL of HS medium with 150 mM methanol (and 2 |g/mL puromycin when plas-mids were present) at 37 °C for 5 days (OD600 ~ 1.0). Cells were collected by centrifugation (5000 rpm for 20 min), and were washed three times with HS medium and puro-mycin alone to remove residual methanol. The final cell pellet was resuspended using 5 mL of HS medium supplemented with 0.1 or 10 mM FeCl3 and 2 |g/mL puromycin when appropriate, to yield a density of 4 x 1010 CFU/mL. After filling the headspace of each tube with methane, the tubes were incubated at 37 °C with mild shaking for 5-10 days.

Cloning ANME-1 mcrBGA

All oligonucleotides are listed in Additional file 1: Table S4. The ANME-1 mcrBGA genes (3.9 kb, locus tag fos0113c9_0022-0024, Genbank accession FP565147.1) encoding the ANME-1 Mcr whose 3D structure has been determined [15], was assembled from six DNA fragments of 600-700 bp (Integrated DNA Technologies) using the Gibson assembly method [29]. Unlike the mcr locus from M. acetivorans [16] and other methanogenic archaea [30], mcrC and mcrD are not present. ANME-1 mcrBGA were cloned downstream of promoter Pcdh using the Xbal and Bmtl sites of pES1 to form pES1-MATmcr1. After electroporating the plasmid into E. coli DH5a-\pir, the complete mcrBGA locus and promoter region were sequenced (via primers veri-p-f, pES1-f, MATmcrB2-f, MATmcrB3-f, MATmcrB4-r, and pES1-r) to confirm no errors were introduced during cloning. For pES1-MATmcr2, in which ANME-1 mcrBGA genes are under the transcription of mcr promoter from M. acetivorans (Pmcr_M. acetivorans), a 419-bp DNA fragment that corresponds to Pmcr_M acetivorans was amplified from the genomic DNA (catalog #35395D-5, American Type Culture Collection, Manassas, VA, USA) using primers Pmcr-f2 and Pmcr-r2. Pmcr_M. acetivorans is further fused to ANME-1 mcrBGA using overlap PCR via primers Pmcr-f2 and B6-r1. To place ANME-1 mcrBGA under control its native promoter (which we named Pmcr_ANME-1), a 237-bp DNA fragment upstream of the ANME-1 mcrBGA genes was synthesized and assembled with mcrBGA to

create Pmcr_ANME-1::mcrANME-1 using overlap PCR (via primers Pmat-f, Pmat-r, and B6-r1). The resulting plas-mid is pES1-MATmcr3. The empty plasmid harboring Pmcr ANME-1 was created by linearizing the vector backbone of pES1-MATmcr3 using primers pES1(Matprom)-f and pES1(Matprom)-r, followed by self-ligation of the backbone to form pES1(Pmat). All plasmids were transformed into M. acetivorans using liposome-mediated transformation [31].

Genome sequencing

Genomic DNA was isolated using the Ultraclean® Microbial DNA Isolation Kit (MO BIO Laboratories, Carlsbad, CA, USA). After shearing the genomic DNA using a Covaris ultrasonicator, DNA fragments were barcoded using TruSeq DNA Nano (Illumina, San Diego, CA, USA). The pooled, barcoded DNA library was sequenced on a MiSeq sequencing platform (Illumina) to generate 2,597,170 paired-end reads of 150 bp for M. acetivorans (ancestral strain) grown on methanol, and 2,699,253 paired-end reads of 150 bp for M. acetivorans/pES1-MATmcr1 grown on methane and 0.1 mM FeCl3.

To identify single-nucleotide polymorphisms, insertions, and deletions between the ancestral strain and the methane-grown strain, sequencing reads of each strain were mapped to the reference genome of M. acetivorans (Genbank accession NC_003552.1) using the Burrows-Wheeler Alignment tool [32]. Aligned reads were sorted based on mapped position in the reference genome using SortSam from the Picard tools (http://broadinstitute. github.io/picard). Misalignments caused by insertions and deletions were corrected locally using IndelRea-ligner from the Genome Analysis ToolKit package [33]. Duplicated reads were marked using MarkDuplicates.jar from the Picard tools to remove sequencing bias. Unified Genotyper [33] was then used to call the variants of each strain after removing variants that are present in less than 40 % of the population. The differences between the ancestral strain and the methane-grown strain were identified using vcftools [34].

16S rDNA amplification and sequencing

To verify archaeal strains, primers ARCH109-F and ARCH934-R [35] were used to amplify an 0.8-kb PCR product of 16S rDNA genes. To verify the absence of bacteria, primers 27F and 1492-R [36] were used to amplify an 1.5-kb PCR product of 16S rDNA genes.

Total protein, cysteine, bicarbonate, and iron reduction assays

Cultures (120 |L) were centrifuged briefly to remove the supernatant, and the cell pellets were resuspended with 10-24 |L of sterile water to lyse the cells. The total

protein concentration of these cell suspensions was determined using the Bradford assay (Bio-Rad Laboratories, Hercules, CA, USA).

The concentration of total cysteine (reduced and oxidized forms) in HS media (in which cells showed methane consumption) was determined spectrophoto-metrically using ninhydrin which reacts specifically with L-cysteine, even in the presence of other thiols, to form a pink-colored product (smax 560 nm) under acidic conditions [37]. The concentration of bicarbonate of HS media (in which cells showed methane consumption) was measured spectrophotometrically using the MaxDiscovery Carbon Dioxide Enzymatic Assay Kit (Bioo Scientific, Austin, TX, USA). Filtered supernatant of each culture was used as a sample solution.

The reduction of iron was measured by employing the ferrozine method for quantification of only reduced iron in the form of Fe2+, as described previously [38]. A 100-|L of culture was immediately mixed with 33 |L of 2 N HCl. Then, 10 |L of the acidified mixture was mixed with 190 |L of 1 mg/mL ferrozine in 100 mM HEPES (pH 7.0), and absorbance at 562 nm was measured. Actual Fe2+ concentrations were calculated by comparing results with those from standard solutions of Fe2+.

Gas chromatography (GC) and high-performance liquid chromatography (HPLC)

GC analyses were conducted for quantifying methane in the culture headspace. Aliquots of 50 or 100 |L volumes were passed through a 6890 N Agilent gas chromatograph equipped with a 60/80 Carboxen-1000 column (4600 x 2.1 mm, Supelco catalog no. 12390-U) and a thermal conductivity detector. The injector, column, and detector were maintained at 150, 180, and 240 °C, respectively. Carrier gas flow (nitrogen) was kept at 20 mL/min, and reference gas flow (also nitrogen) for the detector at 20 mL/min as well. Gases were identified according to their retention times and their concentrations were determined according to comparisons with standards.

HPLC analyses were conducted for the detection and quantification of three organic acids under investigation (acetic acid, formic acid, and pyruvic acid). All samples were filtered through a 0.22 |im polyvinylidene fluoride membrane before diluting 1:6 in running buffer (0.0025 M sulfuric acid in water), then 60 |L of the 1:6 dilution was fractionated by HPLC (Waters 717 autosam-pler with a model 515 pump, and a 2996 photodiode array detector) with a reversed-phase column [Phenomenex Rezex ROA-Organic Acid H + (8 %) (300 x 7.8 mm)]. Separations were conducted using an isocratic flow rate of 0.4 mL/min 0.0025 M sulfuric acid in water. Absorb-ance at 210 nm was used to detect all compounds. Chemicals used as standards for comparisons are glacial acetic

acid (EMD Millipore, catalog no. AX0073-6), sodium formate (catalog no. BP356-100, Fisher Scientific, Hampton, NH, USA), and sodium pyruvate (catalog no. S648-500, Fisher Scientific). Peaks corresponding to those of pyruvic acid, formic acid, and acetic acid were confirmed by retention time, co-elution with standards, and by comparing absorbance spectra with those from the standards. Total quantities of the compounds were calculated by comparing peak areas with standard curves made by running chemical standards.

Western blot

To demonstrate production of ANME-1 Mcr, a FLAG epitope tag was introduced into the carboxy terminus of ANME-1 McrA encoded by mcrA in pES1-MATmcr1, pES1-MATmcr2, and pES1-MATmcr3. The DNA encoding the FLAG tag was incorporated into primer B6-r-flag, which was used along with the respective forward primers to create ANME-1 mcrA-flag. After transformation, M. acetivorans harboring pES1-MATmcr1-flag, pES1-MATmcr2-flag, and pES1-MATmcr3-flag were grown on 200 mL HS-methanol for 5 days, and used for short-duration growth experiments. Methane consumption was measured after 5 days, and cells were harvested by centrifugation. Each cell pellet was resuspended in 2 mL Lysis Buffer [20 mM Tris-HCl, 0.1 mM EDTA, 500 mM s-aminocaproic acid, 10 % glycerol, 1 |L protease inhibitor cocktail (Sigma)]. Cells were sonicated on ice at a power level of 10 for 150 s (30 cycles of 5 s each, 60 Sonic Dismembrator, Fisher Scientific). Total proteins were resolved via 12 % Tris-glycine-SDS gels. Western blots were performed with monoclonal horseradish peroxidase-conjugated antibodies raised against a FLAG epitope tag (Thermo Scientific, Waltham, MA, USA). Blotted proteins were detected using the chemilumines-cence reagents from the SuperSignal West-Pico Chemilu-minescence kit (Thermo Scientific).

13C-labeled methane-grown cultures, 13C-labeled bicarbonate-grown cultures, 13C NMR, and GC/MS

Starter cultures (200 mL) of M. acetivorans/pES1(Pmat) and M. acetivorans/pES1-MATmcr3 were used for short-duration growth experiments. For cultures incubated with 13C-labeled methane, the headspace was filled with 13C-labeled methane (99 % 13C atom, Sigma). HS medium for cultures incubated with 13C-labeled bicarbonate was prepared using 13C-labeled sodium bicarbonate (99 % 13C atom, Cambridge Isotope Laboratories, Tewksbury, MA, USA). All cultures were incubated at 37 °C for 10 days. Acetate was measured using an Agilent 7890A/5975C GC/MSD using a Nukol (Supelco) capillary column (30 x 0.32, 0.25 |im phase thickness comprised of a bonded polyethylene glycol). Incorporation

was determined by integrating the peaks areas corresponding to 12C and 13C acetate. NMR experiments were conducted on a Brüker Avance III HD spectrometer operating at 500.20 MHz and 125.78 MHz for *H and 13C nuclei, respectively, using H2O:D2O as a solvent. 13C and DEPT-135 spectra were recorded with a spectral width of 220 ppm, using 64 K data points, a 90° excitation pulse (11 |is) and relaxation delay of 5 s. 1 k scans were collected and spectra zero-filled to 128 K. For all FIDs, line broadening of 1 Hz was applied prior to Fourier transform. Chemical shifts are reported in ppm from DSS (S = 0). The gradient-selected :H-13C heteronuclear multiple bond correlation (gHMBC) experiment was performed using a low-pass J-filter (3.4 ms) and delays of 65 and 36 ms to observe long-range C-H couplings with 256 increments and 64 transients of 2048 data points. The relaxation delay was 2.0 s. Zero-filling to a 2 K x 2 K matrix and n/2-shifted sine square bell multiplication was performed prior to Fourier transform. Heteronu-clear Single Quantum Coherence (gHSQC) spectra were recorded with 256 increments in F1 and 32 scans per increment, using the standard hsqcetgpsisp.2 Brüker pulse sequence. Relaxation delay of 2 s and 2 K data points was used for spectral width of 10 ppm in the proton dimension, whereas the spectral width in the carbon dimension was 180 ppm.

RNA sequencing

Differential gene analysis of two growth conditions (three biological replicates each) was performed: (1) M. ace-tivorans/pES1-MATmcr3 contacted with methane and (2) M. acetivorans/pES1-MATmcr3 grown on methanol. All starter cultures (200 mL) were grown on methanol for 5 days, and harvested by centrifugation. Cell pellets were washed three times with HS medium, and resuspended using 5 mL HS medium, 2 | g/mL puromycin, and 0.1 mM FeCl3. For condition (1), methane was filled into the headspace of the cultures. For condition (2), 150 mM methanol was added. All cultures were incubated at 37 °C for 5 days, followed by rapid centrifugation in the presence of 50 |L RNAlater solution (Ambion, Austin, TX, USA) per mL of culture. Total RNA isolated using the RNeasy Mini kit (Qiagen, Valencia, CA, USA) was digested with terminator 5'-phosphate-dependent exonu-clease (Epicentre, Madison, WI, USA) to partially remove ribosomal RNA. Digested RNA was cleaned using AgenCourt RNAClean XP beads (AgenCourt Bioscience, Beverly, MA, USA) and used for cDNA library construction using the TruSeq Stranded mRNA Library kit (Illumina). The pooled and barcoded cDNA library was sequenced on a HiSeq sequencing platform (Illumina). Obtained reads were mapped to the reference genome of M. ace-tivorans (Genbank accession NC_003552.1) and plasmid

pES1-MATmcr3 using STAR [39]. The mapped reads were assembled using Cufflink v2.2.1 [40] to identify potential novel transcripts. Assembled, unannotated novel transcripts for all the strains were combined with the list of known genes. Differential expression of genes and potential novel transcripts were determined using Cuffdiff [40] at a significance cutoff at q < 0.07 with a false discovery rate of 0.05. Expression levels of gene transcripts are expressed as fragments per kilobase of transcript per million mapped fragments (FPKM) [41], and expression changes are determined by the ratio of FPKM of culture replicates grown on methane to FPKM of culture replicates grown on methanol. Gene expression data have been deposited in the Gene Expression Omnibus under accession code GSE66445.

Additional file

Additional file 1. This file consists of four supplemental tables and six supplemental figures. Table S1 lists the components in the HS medium used to grow ANME-1 Mcr-producing M. acetivorans on methane and 0.1 mM or 10 mM FeCl3. Table S2 shows the fold changes of differentially expressed genes in M. acetivorans/pES1-MATmcr3 grown on methane and 0.1 mM FeCl3, in comparison to the same strain grown on methanol. Table S3 lists the strains and plasmids, and Table S4 lists the oligonucleotides, used in this study. Figure S1 shows the three promoters used to express ANME-1 mcrBGA genes. Figure S2 shows the detected McrA-FLAG in M. acetivorans/pES1-MATmcr3-flag grown on methane after five days. Figure S3 shows the detection of ANME-1 mcrA after 30 days of growth on methane. Figure S4 shows the GC/MS spectra of culture supernatants used to identify acetate from H13CO. Figure S5 shows the flux through the various reactions in the methanogenesis pathway of M. acetivorans estimated by 13C-metabolic flux analysis using 13C-labeled bicarbonate as the input tracer. Figure S6 shows a simplified methanogenesis pathway from CO2 and CH3OH of M. acetivorans.

Abbreviations

ANME-1: anaerobic methanotrophic archaeal population 1; AOM: anaerobic oxidation of methane; CH3-SCoM: methyl-coenzyme M (2-(methylthio) ethanesulfonate); ATP: adenosine triphosphate; CoBS-SCoM: heterodisulfide of coenzyme B and coenzyme M; Em: resting potential; Fe3+: ferric ion; Fe2+: ferrous ion; FeCl3: ferric chloride; AG^n: Gibbs free energy change upon reaction; HS: high-salt; HSCoB: coenzyme B. CH4 (methane); Mcr: methyl-coenzyme M reductase; Mn4+: manganese ion; NO3-: nitrate ion; NO2-: nitrite ion; Pmcr_

: the mcr promoter from ANME-1; Pcdh: the CO dehydrogenase/acetyl-

CoA synthase promoter from Methanosarcina thermophila; Pr

mcr_M. acetivorans•

Authors' contributions

TKW, JGF, and CDM conceived the initial idea of cloning ANME-1 mcrBGA into M. acetivorans. VWCS, AT, MJM, and TKW designed experiments and interpreted the data. VWCS, AT, MJM, and FZ performed experiments. LZ, EH, and PBS performed 13C NMR and GC/MS and analyzed these data with ADP. SA performed bioinformatics analyses. HN-B and SG performed in silico genome-scale modeling and MFA under the supervision of CDM. HMS designed the ANME-1 gBLOCK DNA sequences. JGF provided facilities for cultivation of M. acetivorans and insights for archaeal physiology. VWCS and TKW wrote the manuscript. All authors read and approved the final manuscript.

Author details

1 Department of Chemical Engineering, The Pennsylvania State University, University Park, PA 16802-4400, USA. 2 Department of Veterinary and Biomedical Sciences, The Pennsylvania State University, University Park, PA 16802-4400,

USA. 3 Department of Chemistry, The Pennsylvania State University, University Park, PA 16802-4400, USA. 4 Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, PA 16802-4400, USA. 5 The Huck Institutes of the Life Sciences, The Pennsylvania State University, University Park, PA 16802-4400, USA. 6 Key Laboratory of Magnetic Resonance in Biological Systems, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China. 7 Institute of Natural and Mathematical Sciences, Massey University, Auckland 0632, New Zealand.

Acknowledgements

This work was supported by the Department of Energy Advanced Research Projects Agency—Energy. We thank Dr. Craig Praul at the Pennsylvania State University Genomics Core Facility for his advice with genome and RNA sequencing. TEM was performed at the Pennsylvania State Microscopy and Cytometry Facility with help from Dr. Brian Kwan. We are also grateful for the comments from Dr. Christopher House. T.K.W. is the Biotechnology Endowed Chair at Pennsylvania State University.

Competing interests

The authors declare that they have no competing interests.

Received: 9 October 2015 Accepted: 13 December 2015 Published online: 14 January 2016

References

1. National Energy Technology Laboratory. DOE's unconventional gas research programs 1976-1995. An archive of important results. Washington, D.C.: US Department of Energy. 2007.

2. US Energy Information Administration. Technically recoverable shale oil and shale gas resources: an assessment of 137 shale formations in 41 countries outside the United States. Washington, D.C.: US Department of Energy; 2013.

3. Solomon S, Qin D, Manning M, Alley R, Berntsen T, Bindoff N, Chen Z, Chidthaison A, Gregory J, Hegerl G, et al. Technical summary. In: Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt K, Tignor M, Miller H, editors. Climate change 2007: the physical science basis contribution of working group I to the fourth assessment report for the intergovernmental panel on climate change. Cambridge, New York: Cambridge University Press; 2007.

4. Conrado RJ, Gonzalez R. Envisioning the bioconversion of methane to liquid fuels. Science. 2014;343:621-3.

5. Wood DA, Nwaoha C, Towler BF. Gas-to-liquids (GTL): a review of an industry offering several routes for monetizing natural gas. J Nat Gas Sci Eng. 2012;9:196-208.

6. Haynes CA, Gonzalez R. Rethinking biological activation of methane and conversion to liquid fuels. Nat Chem Biol. 2014;10:331-9.

7. Knittel K, Boetius A. Anaerobic oxidation of methane: progress with an unknown process. Annu Rev Microbiol. 2009;63:311-34.

8. Dale AW, Van Cappellen P, Aguilera DR, Regnier P. Methane efflux from marine sediments in passive and active margins: estimations from bioenergetic reaction-transport simulations. Earth Planet Sci Lett. 2008;265:329-44.

9. Nauhaus K, Albrecht M, Elvert M, Boetius A, Widdel F. In vitro cell growth of marine archaeal-bacterial consortia during anaerobic oxidation of methane with sulfate. Environ Microbiol. 2007;9:187-96.

10. Hallam SJ, Putnam N, Preston CM, Detter JC, Rokhsar D, Richardson PM, DeLong EF. Reverse methanogenesis: testing the hypothesis with environmental genomics. Science. 2004;305:1457-62.

11. Scheller S, Goenrich M, Boecher R, Thauer RK, Jaun B. The key nickel enzyme of methanogenesis catalyses the anaerobic oxidation of methane. Nature. 2010;465:606-U697.

12. Moran JJ, House CH, Freeman KH, Ferry JG. Trace methane oxidation studied in several Euryarchaeota under diverse conditions. Archaea. 2005;1:303-9.

13. Moran JJ, House CH, Thomas B, Freeman KH. Products of trace methane oxidation during nonmethyltrophic growth by Methanosarcina. J Geo-phys Res. 2007;112:G02011.

ANME-1

mcr promoter from M. acetivorans; SO42 : sulfate ion.

26. 27.

Meyerdierks A, Kube M, Kostadinov I, Teeling H, Glockner FO, Reinhardt R, Amann R. Metagenome and mRNA expression analyses of anaerobic methanotrophic archaea of the ANME-1 group. Environ Microbiol. 2010;12:422-39.

Shima S, Krueger M, Weinert T, Demmer U, Kahnt J, Thauer RK, Ermler U. Structure of a methyl-coenzyme M reductase from Black Sea mats that oxidize methane anaerobically. Nature. 2012;481:98-101. Galagan JE, Nusbaum C, Roy A, Endrizzi MG, Macdonald P, FitzHugh W, Calvo S, Engels R, Smirnov S, Atnoor D, et al. The genome of M. acetivorans reveals extensive metabolic and physiological diversity. Genome Res. 2002;12:532-42.

Kohler PR, Metcalf WW. Genetic manipulation of Methanosarcina spp. Front Microbiol. 2012;3:259.

Macauley SR, Zimmerman SA, Apolinario EE, Evilia C, Hou Y-M, Ferry JG, Sowers KR The archetype y-class carbonic anhydrase (Cam) contains iron when synthesized in vivo. Biochemistry. 2009;48:817-9. Beal EJ, House CH, Orphan VJ. Manganese- and iron-dependent marine methane oxidation. Science. 2009;325:184-7.

Moran JJ, Beal EJ, Vrentas JM, Orphan VJ, Freeman KH, House CH. Methyl sulfides as intermediates in the anaerobic oxidation of methane. Environ Microbiol. 2008;10:162-73.

Vepachedu VR, Ferry JG. Role of the fused corrinoid/methyl transfer protein CmtA during CO-dependent growth of Methanosarcina acetivorans. J Bacteriol. 2012;194:4161-8.

Lienard T, Becher B, Marschall M, Bowien S, Gottschalk G. Sodium ion translocation by N5-methyltetrahydromethanopterin: coenzyme M meth-yltransferase from Methanosarcina mazei Gö1 reconstituted in ether lipid liposomes. Eur J Biochem. 1996;239:857-64.

Raybuck SA, Ramer SE, Abbanat DR, Peters JW, Ormejohnson WH, Ferry JG, Walsh CT. Demonstration of carbon-carbon bond cleavage of acetyl coenzyme A by using isotopic exchange catalyzed by the CO dehydrogenase complex from acetate-grown Methanosarcina thermophila. J Bacteriol. 1991;173:929-32.

Abbanat DR, Ferry JG. Synthesis of acetyl coenzyme A by carbon monoxide dehydrogenase complex from acetate-grown Methanosarcina thermophila. J Bacteriol. 1990;172:7145-50.

Terlesky KC, Ferry JG. Ferredoxin requirement for electron transport from the carbon monoxide dehydrogenase complex to a membrane-bound hydrogenase in acetate-grown Methanosarcina thermophila. J Biol Chem. 1988;263:4075-9.

Matschiavelli N, Rother M. Role of a putative tungsten-dependent formylmethanofuran dehydrogenase in Methanosarcina acetivorans. Arch Microbiol. 2015;197:379-88.

Wang M, Tomb JF, Ferry JG. Electron transport in acetate-grown Methanosarcina acetivorans. BMC Microbiol. 2011;11:165.

Metcalf WW, Zhang J-K, Shi X, Wolfe RS. Molecular, genetic, and biochemical characterization of the serC gene of Methanosarcina barkeri Fusaro. J Bacteriol. 1996;178:5797-802.

Gibson DG, Young L, Chuang RY, Venter JC, Hutchison CA III, Smith HO. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat Methods. 2009;6:343-5.

Weil CF, Cram DS, Sherf BA, Reeve JN. Structure and comparative analysis of the genes encoding component C of methyl coenzyme M reductase in the extremely thermophilic archaebacterium Methanothermus fervidus. J Bacteriol. 1988;170:4718-26.

Metcalf WW, Zhang JK, Apolinario E, Sowers KR, Wolfe RS. A genetic system for Archaea of the genus Methanosarcina: liposome-mediated transformation and construction of shuttle vectors. Proc Natl Acad Sci USA. 1997;94:2626-31.

Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754-60. DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philip-pakis AA, del Angel G, Rivas MA, Hanna M, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43:491-8.

Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Hand-saker RE, Lunter G, Marth GT, Sherry ST, et al. The variant call format and VCFtools. Bioinformatics. 2011;27:2156-8.

Angel R, Claus P, Conrad R Methanogenic archaea are globally ubiquitous in aerated soils and become active under wet anoxic conditions. ISME J. 2012;6:847-62.

Lane DJ. In: Stackebrandt E, Goodfellow M, editors. Nucleic acid techniques in bacterial systematics. Chichester: Wiley. 1991. pp. 115-147. Ohtsu I, Wiriyathanawudhiwong N, Morigasaki S, Nakatani T, Kado-kura H, Takagi H. The L-cysteine/L-cystine shuttle system provides reducing equivalents to the periplasm in Escherichia coli. J Biol Chem. 2010;285:17479-87.

Ruebush SS, Brantley SL, Tien M. Reduction of soluble and insoluble iron forms by membrane fractions of Shewanella oneidensis grown under aerobic and anaerobic conditions. Appl Environ Microbiol. 2006;72:2925-35. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15-21.

Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562-78.

Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511-5.

Submit your next manuscript to BioMed Central and we will help you at every step:

• We accept pre-submission inquiries

• Our selector tool helps you to find the most relevant journal

• We provide round the clock customer support

• Convenient online submission

• Thorough peer review

• Inclusion in PubMed and all major indexing services

• Maximum visibility for your research

Submit your manuscript at ■v

www.biomedcentral.com/submit Central