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 FeCl3 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

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 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 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 Methanothermobacter 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 dehydrogenase/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-1, as shown by the appearance of ANME-1 McrA-FLAG (the α 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 µmol 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.

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. acetivorans/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 (No cells) 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

Of the tested electron acceptors, the best growth on methane was seen with 10 mM Fe3+. After 6 weeks of incubation using low-density inocula, growth of M. acetivorans/pES1-MATmcr3 was observed in conjunction with 9 ± 1 % of methane consumption (corresponding to 109 ± 12 µmol 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/pES1-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).

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. acetivorans/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

The identity of M. acetivorans during growth on methane was confirmed by 16S ribosomal DNA (rDNA) sequencing (via PCR amplification using archaeal 16S 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 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 (\( \Delta {\text{G}}_{\text{rxn}}^{ \circ } \) = −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-1 [10.3 ± 0.8 mM corresponding to 52 ± 4 µmol, 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 multi-dimensional nuclear magnetic resonance (NMR) and by gas chromatography coupled with mass spectroscopy (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.

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. aM. acetivorans harboring empty plasmid pES1(Pmat) grown on 13C-labeled methane (13CH4) for 10 days. bM. 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

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 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).

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).

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.

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 formylmethanofuran, 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

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 methyltransferase (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 CH3COSCoA 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 methanofuran) 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 ferredoxin (Em = −420 mV) to the favorable reduction of Fe3+ (Em = +770 mV). Importantly, the system is postulated 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.

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. acetivorans/pES1-MATmcr3 cells, growth on methane was possible only via iron reduction, which highlights the distinctiveness of our engineered population of M. acetivorans/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 culturability 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.

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. acetivorans, 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/pES1(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) × 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 plasmids 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 puromycin 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 × 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 XbaI and BmtI sites of pES1 to form pES1-MATmcr1. After electroporating the plasmid into E. coli DH5α-λ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 plasmid 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].

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 IndelRealigner 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 spectrophotometrically using ninhydrin which reacts specifically with l-cysteine, even in the presence of other thiols, to form a pink-colored product (εmax 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+.

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 × 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 µm 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 autosampler with a model 515 pump, and a 2996 photodiode array detector) with a reversed-phase column [Phenomenex Rezex ROA-Organic Acid H + (8 %) (300 × 7.8 mm)]. Separations were conducted using an isocratic flow rate of 0.4 mL/min 0.0025 M sulfuric acid in water. Absorbance 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 ε-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 chemiluminescence reagents from the SuperSignal West-Pico Chemiluminescence kit (Thermo Scientific).

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 × 0.32, 0.25 μm 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 1H 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 µs) 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 (δ = 0). The gradient-selected 1H-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 × 2 K matrix and π/2-shifted sine square bell multiplication was performed prior to Fourier transform. Heteronuclear 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. acetivorans/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 exonuclease (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. acetivorans (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.

Declarations

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.

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.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://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.

Additional files

12934_2015_397_MOESM1_ESM.pdfAdditional 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.

Authors’ Affiliations

(1)

Department of Chemical Engineering, The Pennsylvania State University, University Park, USA

(2)

Department of Veterinary and Biomedical Sciences, The Pennsylvania State University, University Park, USA

(3)

Department of Chemistry, The Pennsylvania State University, University Park, USA

(4)

Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, USA

(5)

The Huck Institutes of the Life Sciences, The Pennsylvania State University, University Park, USA