Intramolecular isotopic evidence for bacterial oxidation of propane in subsurface natural gas reservoirs

Nitrate-driven anaerobic oxidation of ethane and butane by bacteria

The short-chain gaseous alkanes (ethane, propane, and butane; SCGAs) are important components of natural gas, yet their fate in environmental systems is poorly understood. Microbially mediated anaerobic oxidation of SCGAs coupled to nitrate reduction has been demonstrated for propane, but is yet to be shown for ethane or butane-despite being energetically feasible. Here we report two independent bacterial enrichments performing anaerobic ethane and butane oxidation, respectively, coupled to nitrate reduction to dinitrogen gas and ammonium. Isotopic 13C- and 15N-labelling experiments, mass and electron balance tests, and metabolite and meta-omics analyses collectively reveal that the recently described propane-oxidizing "Candidatus Alkanivorans nitratireducens" was also responsible for nitrate-dependent anaerobic oxidation of the SCGAs in both these enrichments. The complete genome of this species encodes alkylsuccinate synthase genes for the activation of ethane/butane via fumarate addition. Further substrate range tests confirm that "Ca. A. nitratireducens" is metabolically versatile, being able to degrade ethane, propane, and butane under anoxic conditions. Moreover, our study proves nitrate as an additional electron sink for ethane and butane in anaerobic environments, and for the first time demonstrates the use of the fumarate addition pathway in anaerobic ethane oxidation. These findings contribute to our understanding of microbial metabolism of SCGAs in anaerobic environments.

Carbon and hydrogen stable isotope fractionation due to monooxygenation of short-chain alkanes by butane monooxygenase of Thauera butanivorans Bu-B1211

Multi element compound-specific stable isotope analysis (ME-CSIA) is a tool to assess (bio)chemical reactions of molecules in the environment based on their isotopic fingerprints. To that effect, ME-CSIA concepts are initially developed with laboratory model experiments to determine the isotope fractionation factors specific for distinct (bio)chemical reactions. Here, we determined for the first time the carbon and hydrogen isotope fractionation factors for the monooxygenation of the short-chain alkanes ethane, propane, and butane. As model organism we used Thauera butanivorans strain Bu-B1211 which employs a non-haem iron monooxygenase (butane monooxygenase) to activate alkanes. Monooxygenation of alkanes was associated with strong carbon and hydrogen isotope effects: εbulkC = -2.95 ± 0.5 ‰ for ethane, -2.68 ± 0.1 ‰ for propane, -1.19 ± 0.18 ‰ for butane; εbulkH = -56.3 ± 15 ‰ for ethane, -40.5 ± 2.3 ‰ for propane, -14.6 ± 3.6 ‰ for butane. This resulted in lambda (Λ ≈ εHbulk/εCbulk) values of 16.2 ± 3.7 for ethane, 13.2 ± 0.7 for propane, and 11.4 ± 2.8 for butane. The results show that ME-CSIA can be used to track the occurrence and impact of monooxygenase-dependent aerobic processes converting short-chain alkanes in natural settings like marine and terrestrial seeps, gas reservoirs, and other geological formations impacted by natural gas.

Establishing the accuracy of position-specific carbon isotope analysis of propane by GC-pyrolysis-GC-IRMS

RATIONALE

Position-specific (PS) δ¹³ C values of propane have proven their ability to provide valuable information on the evolution history of natural gases. Two major approaches to measure PS δ¹³ C values of propane are isotopic ¹³ C nuclear magnetic resonance (NMR) and gas chromatography-pyrolysis-gas chromatography-isotope ratio mass spectrometry (GC-Py-GC-IRMS). Measurement accuracy of the isotopic ¹³ C NMR has been verified, but the requirements of large sample size and long experimental time limit its applications. GC-Py-GC-IRMS is a more versatile method with a small sample size, but its accuracy has not been demonstrated.

METHODS

We measured the PS δ¹³ C values of propane from nine natural gases using both ¹³ C NMR and GC-Py-GC-IRMS, then evaluated the accuracy of the GC-Py-GC-IRMS method.

RESULTS

The results show that large carbon isotope fractionations occurred for both terminal and central carbons within propane during pyrolysis. The isotope fractionations during the pyrolysis are reproducible at optimum conditions, but vary between the two GC-Py-GC-IRMS systems tested, affected by experimental conditions (e.g., pyrolysis temperature, flow rate, and reactor conditions).

CONCLUSIONS

It is necessary to evaluate and calibrate each GC-Py-GC-IRMS system using propane gases with accurately determined PS δ¹³ C values. This study also highlights a need for PS isotope standards for propane and other molecules (e.g., butane and acetic acid).

A global perspective on bacterial diversity in the terrestrial deep subsurface

While recent efforts to catalogue Earth's microbial diversity have focused upon surface and marine habitats, 12-20 % of Earth's biomass is suggested to exist in the terrestrial deep subsurface, compared to ~1.8 % in the deep subseafloor. Metagenomic studies of the terrestrial deep subsurface have yielded a trove of divergent and functionally important microbiomes from a range of localities. However, a wider perspective of microbial diversity and its relationship to environmental conditions within the terrestrial deep subsurface is still required. Our meta-analysis reveals that terrestrial deep subsurface microbiota are dominated by Betaproteobacteria, Gammaproteobacteria and Firmicutes, probably as a function of the diverse metabolic strategies of these taxa. Evidence was also found for a common small consortium of prevalent Betaproteobacteria and Gammaproteobacteria operational taxonomic units across the localities. This implies a core terrestrial deep subsurface community, irrespective of aquifer lithology, depth and other variables, that may play an important role in colonizing and sustaining microbial habitats in the deep terrestrial subsurface. An in silico contamination-aware approach to analysing this dataset underscores the importance of downstream methods for assuring that robust conclusions can be reached from deep subsurface-derived sequencing data. Understanding the global panorama of microbial diversity and ecological dynamics in the deep terrestrial subsurface provides a first step towards understanding the role of microbes in global subsurface element and nutrient cycling.

Deep Subseafloor Biogeochemical Processes and Microbial Populations Potentially Associated with the 2011 Tohoku-oki Earthquake at the Japan Trench Accretionary Wedge (IODP Expedition 343)

Post-mega-earthquake geochemical and microbiological properties in subseafloor sediments of the Japan Trench accretionary wedge were investigated using core samples from Hole C0019E, which was drilled down to 851‍ ‍m below seafloor (mbsf) at a water depth of 6,890 m. Methane was abundant throughout accretionary prism sediments; however, its concentration decreased close to the plate boundary decollement. Methane isotope systematics indicated a biogenic origin. The content of mole-cular hydrogen (H₂) was low throughout core samples, but markedly increased at specific depths that were close to potential faults predicted by logging-while-drilling ana-lyses. Based on isotopic systematics, H₂ appeared to have been abundantly produced via a low-temperature interaction between pore water and the fresh surface of crushed rock induced by earthquakes. Subseafloor microbial cell density remained constant at approximately 10⁵‍ ‍cells‍ ‍mL^-1. Amplicon sequences revealed that predominant members at the phylum level were common throughout the units tested, which also included members frequently found in anoxic subseafloor sediments. Metabolic potential assays using radioactive isotopes as tracers revealed homoacetogenic activity in H₂-enriched core samples collected near the fault. Furthermore, homoacetogenic bacteria, including Acetobacterium carbinolicum, were isolated from similar samples. Therefore, post-earthquake subseafloor microbial communities in the Japan Trench accretionary prism appear to be episodically dominated by homoacetogenic populations and potentially function due to the earthquake-induced low-temperature generation of H₂. These post-earthquake microbial communities may eventually return to the steady-state communities dominated by oligotrophic heterotrophs and hydrogenotrophic and methylotrophic methanogens that are dependent on refractory organic matter in the sediment.

Anaerobic oxidation of propane coupled to nitrate reduction by a lineage within the class Symbiobacteriia

Anaerobic microorganisms are thought to play a critical role in regulating the flux of short-chain gaseous alkanes (SCGAs; including ethane, propane and butane) from terrestrial and aquatic ecosystems to the atmosphere. Sulfate has been confirmed to act as electron acceptor supporting microbial anaerobic oxidation of SCGAs, yet several other energetically more favourable acceptors co-exist with these gases in anaerobic environments. Here, we show that a bioreactor seeded with biomass from a wastewater treatment facility can perform anaerobic propane oxidation coupled to nitrate reduction to dinitrogen gas and ammonium. The bioreactor was operated for more than 1000 days, and we used ¹³C- and ¹⁵N-labelling experiments, metagenomic, metatranscriptomic, metaproteomic and metabolite analyses to characterize the microbial community and the metabolic processes. The data collectively suggest that a species representing a novel order within the bacterial class Symbiobacteriia is responsible for the observed nitrate-dependent propane oxidation. The closed genome of this organism, which we designate as 'Candidatus Alkanivorans nitratireducens', encodes pathways for oxidation of propane to CO₂ via fumarate addition, and for nitrate reduction, with all the key genes expressed during nitrate-dependent propane oxidation. Our results suggest that nitrate is a relevant electron sink for SCGA oxidation in anaerobic environments, constituting a new microbially-mediated link between the carbon and nitrogen cycles.

Low 13C-13C abundances in abiotic ethane

Distinguishing biotic compounds from abiotic ones is important in resource geology, biogeochemistry, and the search for life in the universe. Stable isotopes have traditionally been used to discriminate the origins of organic materials, with particular focus on hydrocarbons. However, despite extensive efforts, unequivocal distinction of abiotic hydrocarbons remains challenging. Recent development of clumped-isotope analysis provides more robust information because it is independent of the stable isotopic composition of the starting material. Here, we report data from a ¹³C-¹³C clumped-isotope analysis of ethane and demonstrate that the abiotically-synthesized ethane shows distinctively low ¹³C-¹³C abundances compared to thermogenic ethane. A collision frequency model predicts the observed low ¹³C-¹³C abundances (anti-clumping) in ethane produced from methyl radical recombination. In contrast, thermogenic ethane presumably exhibits near stochastic ¹³C-¹³C distribution inherited from the biological precursor, which undergoes C-C bond cleavage/recombination during metabolism. Further, we find an exceptionally high ¹³C-¹³C signature in ethane remaining after microbial oxidation. In summary, the approach distinguishes between thermogenic, microbially altered, and abiotic hydrocarbons. The ¹³C-¹³C signature can provide an important step forward for discrimination of the origin of organic molecules on Earth and in extra-terrestrial environments.

Hydrocarbon Cycling in the Tokamachi Mud Volcano (Japan): Insights from Isotopologue and Metataxonomic Analyses

Understanding hydrocarbon cycling in the subsurface is important in various disciplines including climate science, energy resources and astrobiology. Mud volcanoes provide insights into biogeochemical processes occurring in the subsurface. They are usually associated with natural gas reservoirs consisting mainly of methane and other hydrocarbons as well as CO₂. Stable isotopes have been used to decipher the sources and sinks of hydrocarbons in the subsurface, although the interpretation can be ambiguous due to the numerous processes involved. Here we report new data for hydrocarbon isotope analysis, including position-specific isotope composition of propane, for samples from the Tokamachi mud volcano area, Japan. The data suggest that C₂₊ hydrocarbons are being biodegraded, with indirect production of methane (“secondary methanogenesis”). Data from chemical and isotopic composition are discussed with regard to 16S rRNA analysis, which exhibits the presence of hydrogenotrophic and acetoclastic methoanogens. Overall, the combination of isotopologue analysis with 16S rRNA gene data allows refining of our understanding of hydrocarbon cycling in subsurface environments.

Genome and proteome analyses show the gaseous alkane degrader Desulfosarcina sp. strain BuS5 as an extreme metabolic specialist

The metabolic potential of the sulfate-reducing bacterium Desulfosarcina sp. strain BuS5, currently the only pure culture able to oxidize the volatile alkanes propane and butane without oxygen, was investigated via genomics, proteomics and physiology assays. Complete genome sequencing revealed that strain BuS5 encodes a single alkyl-succinate synthase, an enzyme which apparently initiates oxidation of both propane and butane. The formed alkyl-succinates are oxidized to CO₂ via beta oxidation and the oxidative Wood-Ljungdahl pathways as shown by proteogenomics analyses. Strain BuS5 conserves energy via the canonical sulfate reduction pathway and electron bifurcation. An ability to utilize long-chain fatty acids, mannose and oligopeptides, suggested by automated annotation pipelines, was not supported by physiology assays and in-depth analyses of the corresponding genetic systems. Consistently, comparative genomics revealed a streamlined BuS5 genome with a remarkable paucity of catabolic modules. These results establish strain BuS5 as an exceptional metabolic specialist, able to grow only with propane and butane, for which we propose the name Desulfosarcina aeriophaga BuS5. This highly restrictive lifestyle, most likely the result of habitat-driven evolutionary gene loss, may provide D. aeriophaga BuS5 a competitive edge in sediments impacted by natural gas seeps. Etymology: Desulfosarcina aeriophaga, aério (Greek): gas; phágos (Greek): eater; D. aeriophaga: a gas eating or gas feeding Desulfosarcina.

Standardization for 13 C-13 C clumped isotope analysis by the fluorination method

RATIONALE

The ¹³ C-¹³ C isotopologues of C₂ molecules have recently been measured using a fluorination method. The C₂ compound is first fluorinated into hexafluoroethane (C₂ F₆ ), and its ¹³ C-isotopologues are subsequently measured using a conventional isotope ratio mass spectrometer. Here, we present an approach for standardizing the fluorination method on an absolute reference scale by using isotopically enriched C₂ F₆ .

METHODS

We prepared physical mixtures of ¹³ C-¹³ C-labeled ethanol and natural ethanol. The enriched ethanol samples were measured using the recently developed fluorination method. Based on the difference between the calculated and measured ∆¹³ C¹³ C values, we quantified the extent to which isotopologues were scrambled during dehydration, fluorination, and ionization in the ion source.

RESULTS

The measured ∆¹³ C¹³ C value was approximately 20% lower than that expected from the amount of ¹³ C-¹³ C ethanol. The potential scrambling in the ion source was estimated to be 0.5-2%, which is lower than the observed isotopic reordering. Therefore, isotopic reordering may have occurred during either dehydration or fluorination.

CONCLUSIONS

For typical analysis of natural samples, scrambling in the ion source can only change the ∆¹³ C¹³ C value by less than 0.04‰, which is lower than the current analytical precision (±0.07‰). Therefore, the observed isotopic reordering may have occurred during the fluorination of ethene through the scrambling of isotopologues of ethene but not in the ion source of the mass spectrometer or during the dehydration of ethanol, given the small amount of C₁ and C₃₊ molecules. Thus, we obtained the empirical transfer function ∆¹³ C¹³ C_CSC  = λ × ∆¹³ C¹³ C with a λ value of 1.25 ± 0.01 for ethanol/ethene and 1.00 for ethane. Using the empirical transfer function, the developed fluorination method can provide actual differences in ∆ values.

The Grayness of the Origin of Life

In the search for life beyond Earth, distinguishing the living from the non-living is paramount. However, this distinction is often elusive, as the origin of life is likely a stepwise evolutionary process, not a singular event. Regardless of the favored origin of life model, an inherent "grayness" blurs the theorized threshold defining life. Here, we explore the ambiguities between the biotic and the abiotic at the origin of life. The role of grayness extends into later transitions as well. By recognizing the limitations posed by grayness, life detection researchers will be better able to develop methods sensitive to prebiotic chemical systems and life with alternative biochemistries.

A fluorination method for measuring the 13 C-13 C isotopologue of C2 molecules

RATIONALE

Doubly substituted isotope species ("clumped" isotopes) can provide insights into the biogeochemical history of a molecule, including its temperature of formation and/or its (bio)synthetic pathway. Here, we propose a new fluorination method for the measurement of ¹³ C-¹³ C species in C₂ molecules using a conventional isotope ratio mass spectrometer. Target molecules include ethane, ethene and ethanol.

METHODS

¹³ C-¹³ C isotope species in C₂ molecules were measured as C₂ F₆ using a conventional isotope ratio mass spectrometer. Ethane and ethene are directly fluorinated to C₂ F₆ . Ethanol is measured after dehydration to ethene and subsequent fluorination of the latter. The method enables the measurement of the Δ¹³ C¹³ C values normalized against a reference working standard.

RESULTS

The reproducibility of the whole protocol, including chemical modification steps and measurement of C₂ F₆ isotopologues, is better than ±0.14‰ for all the compounds. Ethane from natural gas samples and biologically derived ethanol show a narrow range of Δ¹³ C¹³ C values, varying from 0.72‰ to 0.90‰. In contrast, synthetic ethanol as well as putative abiotic ethane show Δ¹³ C¹³ C values significantly different from this range with values of 1.14‰ and 0.25‰, respectively.

CONCLUSIONS

The method presented here provides alternative means of measuring ¹³ C-¹³ C species to that using high-resolution mass spectrometry. Preliminary data from natural and synthetic molecules re-emphasizes the potential of ¹³ C clumped isotope species as a (bio)marker.