# An Origin-of-Life Reactor to Simulate Alkaline Hydrothermal Vents

Barry Herschy · Alexandra Whicher · Eloi Camprubi · Cameron Watson · Lewis Dartnell · John Ward · Julian R. G. Evans · Nick Lane

Received: 30 October 2014 / Accepted: 21 November 2014 / Published online: 27 November 2014
© The Author(s) 2014. This article is published with open access at Springerlink.com

**Abstract** Chemiosmotic coupling is universal: practically all cells harness electrochemical proton gradients across membranes to drive ATP synthesis, powering biochemistry. Autotrophic cells, including phototrophs and chemolithotrophs, also use proton gradients to power carbon fixation directly. The universality of chemiosmotic coupling suggests that it arose very early in evolution, but its origins are obscure. Alkaline hydrothermal systems sustain natural proton gradients across the thin inorganic barriers of interconnected micropores within deep-sea vents. In Hadean oceans, these inorganic barriers should have contained catalytic Fe(Ni)S minerals similar in structure to cofactors in modern metabolic enzymes, suggesting a possible abiotic origin of chemiosmotic coupling. The continuous supply of H₂ and CO₂ from vent fluids and early oceans, respectively, offers further parallels with the biochemistry of ancient autotrophic cells, notably the acetyl CoA pathway in archaea and bacteria. However, the precise mechanisms by which natural proton gradients, H₂,

CO₂ and metal sulphides could have driven organic synthesis are uncertain, and theoretical ideas lack empirical support. We have built a simple electrochemical reactor to simulate conditions in alkaline hydrothermal vents, allowing investigation of the possibility that abiotic vent chemistry could prefigure the origins of biochemistry. We discuss the construction and testing of the reactor, describing the precipitation of thin-walled, inorganic structures containing nickel-doped mackinawite, a catalytic Fe(Ni)S mineral, under prebiotic ocean conditions. These simulated vent structures appear to generate low yields of simple organics. Synthetic microporous matrices can concentrate organics by thermophoresis over several orders of magnitude under continuous open-flow vent conditions.

**Keywords** Serpentinization · Alkaline vent · Hydrothermal · CO₂ reduction · Thermophoresis

## The Origins of Biochemistry

Does the biochemistry of modern cells offer any insight into its own origins? There are good reasons to seek congruence between the origins of biochemistry and particular early Earth environments. de Duve (2005) has asked how biological catalysts, whether enzymes or ribozymes, might have first arisen, to which his “only scientifically plausible explanation” was “by selection”. If so, the first biological catalysts must have been selected in the context of protometabolism, meaning that they enhanced a process that occurred spontaneously, driven by natural disequilibria and catalysed by inorganic catalysts. Selection therefore imposes a link between protometabolism and metabolism, providing good grounds for seeking just such a link. This argument is at least parsimonious, and advocates that we

<footnote>B. Herschy · A. Whicher · E. Camprubi · C. Watson · N. Lane
(✉)
Department of Genetics, Evolution and Environment, University College London, London, UK
e-mail: nick.lane@ucl.ac.uk</footnote>

<footnote>L. Dartnell
Space Research Centre, Department of Physics and Astronomy,
University of Leicester, Leicester, UK</footnote>

<footnote>J. Ward
Department of Biochemical Engineering, University College
London, London, UK</footnote>

<footnote>J. R. G. Evans
Department of Chemistry, University College London, London,
UK</footnote>
should look to life itself for clues to how life arose (Martin et al. 2014). One factor in particular distinguishes living cells from conventional synthetic chemistry as practised by humans: biochemistry is fundamentally vectorial. It has structure and direction in space, as pointed out by Peter Mitchell from the late 1950s onwards (Mitchell 1959, 1961, 1966). Most importantly, in all known autotrophic bacteria and archaea, carbon and energy metabolism is driven by electrochemical ion (generally proton) gradients across membranes, Mitchell's chemiosmotic coupling (Maden 1995; Stetter 2006; Lane et al. 2010). Geochemical systems with analogous vectorial chemistry could therefore offer valuable insights into the origin of life.

Searching for congruence between geochemistry and biochemistry may be reasonable, but which aspects of biochemistry are genuinely ancient, and which are more recent adaptations to specific (potentially ancient) environments? It is nearly impossible to build a phylogenetic tree of life to reconstruct the biochemistry of the Last Universal Common Ancestor (LUCA), as frequent lateral gene transfer (LGT) in bacteria and archaea has founded the order of deep branches (Doolittle 1999; Doolittle and Bapteste 2007; Dagan and Martin 2007; Sousa et al. 2013; Nelson-Sathi et al. 2014). However, deep biochemical differences between the domains of life do give some insights. Because eukaryotes are secondarily derived from an endosymbiosis between an archaeal host cell and bacterial endosymbionts (Cox et al. 2008; Williams et al. 2013), LUCA was the common ancestor of bacteria and archaea (Dagan et al. 2010). While bacteria and archaea share a number of fundamental traits including transcription, ribosomal translation, aspects of amino acid biosynthesis and chemiosmotic coupling using a membrane-integral ATP synthase, several other fundamental traits are shockingly different (Martin and Russell 2003; Sousa et al. 2013). For example, most enzymes involved in DNA replication are not homologous in bacteria and archaea (Edgell and Doolittle 1997; Leipe et al. 1999).

Most importantly, the cell membrane and cell wall are radically different in their chemistry and stereochemistry (Koga et al. 1998; Peretó et al. 2004), as are the enzymes for lipid biosynthesis (Martin and Russell 2003; Peretó et al. 2004; Lombard et al. 2012). Biochemical pathways such as glycolysis (Say and Fuchs 2010), and heme and quinone synthesis (Sousa et al. 2013) are also distinct, albeit with some confounding LGT. Such considerations imply that fermentation and respiration, or more specifically the heme and quinone-dependent proteins used for respiratory ion pumping, arose independently in the bacteria and archaea, but see Ducluzeau et al. (2014) for a contrary view.

These paradoxical properties could be resolved if LUCA depended on natural (geochemically sustained) proton

gradients to drive carbon and energy metabolism, a life-style demanding membranes extremely leaky to protons (Lane and Martin 2012; Sojo et al. 2014). This requirement for leaky membranes could explain the divergence in other traits that might have coevolved later with membranes (see Sojo et al. 2014) notably DNA replication (in which the replicon is attached to the membrane in most bacteria), and cell wall synthesis, which requires membrane-integral export machinery.

Despite these stark disparities, some clues do exist to early biochemistry. Strikingly, there are only six known pathways of carbon fixation across all life (Fuchs 2011), but just one of these, the acetyl CoA pathway, is found in both archaea (methanogens) and bacteria (acetogens), albeit with some striking differences between them (Maden 2000; Martin 2012; Sousa and Martin 2014). Neglecting these differences, several other factors testify to the antiquity of the acetyl CoA pathway. It is the only exergonic pathway of carbon fixation, drawing on just H₂ and CO₂ as substrates to drive both carbon and energy metabolism (Fuchs and Stupperich 1985; Ragsdale and Pierce 2008; Ljungdahl 2009); what Everett Shock has called “a free lunch you’re paid to eat” (Shock et al. 1998). It is short and linear, with just a few steps leading from H₂ and CO₂ to acetyl CoA and pyruvate, the gateway to intermediary metabolism (Fuchs 2011; Morowitz et al. 2000), thereby avoiding the problem of sequentially declining yields with non-enzymic cycles that might have precluded an abiotic reductive TCA cycle (Orgel 2008). It is replete with Fe(Ni)S proteins, in which the inorganic cofactors that actually catalyse the key reactions in proteins such as ferredoxin and hydrogenases have structures essentially identical to FeS minerals such as greigite and mackinawite (Eck and Dayhoff 1966; Russell and Martin 2004; Ragsdale and Kumar 1996; Cody 2004; Major et al. 2004; Baymann et al. 2003; Kim et al. 2013; Harel et al. 2014). Finally, despite being exergonic overall, the acetyl CoA pathway is strictly dependent on chemiosmotic coupling in both methanogens and acetogens (Thauer et al. 2007; Poehlein et al. 2012). In the case of methanogens, the H⁺ (or Na⁺) gradient drives ferredoxin reduction via the Fe(Ni)S membrane protein, Ech (Kaster et al. 2011, Buckel and Thauer 2013).

All these factors—H₂, CO₂, Fe(Ni)S catalysis and electrochemical proton gradients—point to one very specific environment on the early Earth as the cradle of life: alkaline hydrothermal vents (Martin et al. 2014).

## Alkaline Hydrothermal Vents

### Early Vents as Electrochemical Reactors

Russell and colleagues (Russell et al. 1988, 1989, 1993, 1994; Russell and Hall 1997) predicted the existence and
properties of deep-ocean alkaline hydrothermal systems more than a decade before their discovery, pointing out their suitability as natural electrochemical reactors capable of driving the origin of life. While fossil vent systems had been reported in Ireland (Boyce et al. 1983), the discovery of the first active submarine system, Lost City hydrothermal field (Kelley et al. 2001, 2005), was remarkable in that its properties corresponded almost exactly to those postulated by Russell et al. (1993). Lost City is powered by a process called serpentinization, the exothermic reaction of ultramafic minerals from the upper mantle, in particular olivine, with

water (Bach et al. 2006; Sleep et al. 2004; Martin et al. 2008; Russell et al. 2010). This reaction produces large volumes of H₂ (the presence of the mineral awaruite in some serpentinizing systems indicating as much as 200 mM (McCollom and Bach 2009) dissolved in warm (45–90 °C) alkaline (pH 9–11) fluids containing magnesium hydroxides (Kelley et al. 2001, 2005). Alkaline vents do not form chimneys, as in black smokers (and indeed do not normally ‘smoke’ at all) but rather are labyrinthine networks of interconnected micropores bounded by thin inorganic walls, through which hydrothermal fluids (and ocean waters) percolate.

**Fig. 1** *a* Standard reduction potentials of H₂ and CO₂ at pH 7. Transfer of electrons from H₂ to CO₂ is unfavourable as the reduction potential for CO₂ at this pH is lower (more negative) than H₂. *b* With H₂ dissolved in waters at pH 10 and dissolved CO₂ in waters at pH 6 however, the reduction potential for CO₂ becomes higher (more positive) than that of H₂ making the reduction of CO₂ favourable. This would theoretically allow for the reduction of CO₂ to form organic compounds such as formate, formaldehyde, methanol and methane. *c* How acid and alkaline fluids could interact inside hydrothermal vents across thin semi-conducting Fe(Ni)S walls, leading to the reduction of CO₂ to formaldehyde via formate

<figure><img src="image_1.png" /></figure>
Such vents should have been more common on the early Earth, as the mantle was less differentiated from the crust, hence ultramafic minerals could have been found across much of the ocean floor (Fyfe 1994; Jaffrés et al. 2007; Shields and Kasting 2007). In contrast, ultramafic minerals are mostly exposed close to the mid-ocean spreading centres today (Schrenk et al. 2013). Alkaline vents are highly stable geological systems; Lost City is estimated to be about 100,000-years old (Ludwig et al. 2005), which as noted by Russell, is $10^{17}$ microseconds, a time unit more consonant with chemistry. That gives plenty of time for abiotic chemistry to develop, especially if early vents were indeed contiguous across the ocean floor (Sleep 2010; Shields and Kasting 2007). Moreover, the fact that olivine and water are both abundant in space (de Leeuw et al. 2010), and so presumably on all wet, rocky Earth-like planets, implies that equivalent conditions could be projected to occur on as many as 40 billion exoplanets in the Milky Way alone (Lane 2015).

Lost City is composed of carbonate minerals, mostly aragonite, and magnesium hydroxide, brucite, (Kelley et al. 2001, 2005) but this is unlikely to represent the composition of ancient vents. That difference is critical and relates not to serpentinization as a process (which should have been the same), but to ocean chemistry in the Hadean and Archaean, around 4 billion years ago (Pinti 2005). There were two critical differences: oxygen was absent (Bekker et al. 2004; Kasting 2013); and the CO₂ concentration in the oceans was substantially higher (although there is little consensus on how much higher; see Russell and Arndt 2005; Sleep 2010; Arndt and Nisbet 2012). Anoxia is necessary for both thermodynamic and kinetic reasons. Thermodynamic, because the reaction between H₂ and CO₂ is only favoured under anoxic conditions (Amend et al. 2013); and kinetic, because the solubility of catalytic transition metals, notably Fe²⁺ and Ni²⁺ is much greater when oceans are anoxic (Russell and Arndt 2005; Arndt and Nisbet 2012). That the Hadean oceans were indeed replete in Fe²⁺ and Ni²⁺ (derived from volcanic systems such as black smokers) are indicated by the precipitation of vast banded-iron formations throughout the Archaean (Anbar and Holland 1992; Zahnle et al. 2007). The great availability of transition metals (along with bisulphide ions within alkaline vents; Nitschke and Russell 2009) must have resulted in the precipitation of catalytic Fe(Ni)S minerals such as mackinawite and greigite in the walls of the vents themselves; but equivalent catalytic Fe(Ni)S minerals are not found in modern vents. In early vents then, H₂-rich hydrothermal fluids must have percolated through labyrinths of micropores bounded by thin inorganic walls containing catalytic Fe(Ni)S minerals (Nitschke and Russell 2009; Lane and Martin 2012).

The higher CO₂ concentration in Hadean oceans should have increased carbon availability (modern alkaline hydrothermal vents are often carbon limited, from carbonate precipitation and removal by living cells; Proskurowski et al. 2008; Bradley et al. 2009) and lowered the pH of the oceans, probably to around pH 5–7 (Arndt and Nisbet 2012). That could have produced pH gradients of 5 or 6 pH units between the alkaline hydrothermal fluids and acidic oceans. While mixing could prevent such steep gradients being juxtaposed across single barriers, laminar flow in elongated hydrothermal pores does make it feasible for sharp gradients of several pH units to exist across distances of a few micrometres.

## The Driving Force for Organic Synthesis

Steep natural proton gradients across thin catalytic Fe(Ni)S barriers could theoretically promote organic synthesis by lowering the energetic barrier to CO₂ reduction (Lane 2014; Yamaguchi et al. 2014). Amend and McCollom (2009) calculated that anoxic alkaline hydrothermal conditions (between 25 and 125 °C) are thermodynamically conducive to the synthesis of total cell biomass (i.e. amino acids, fatty acids, carbohydrates, nucleotides) from H₂ and CO₂. Nonetheless, experimental attempts to drive the reaction of H₂ and CO₂ using Fe(Ni)S catalysts have proved unsuccessful, even at high pressures (Shock and Canovas 2010), as the reduction potential of the H₂/2H⁺ couple is not sufficiently low to reduce CO₂ to CO, formate (HCOO⁻), formaldehyde (HCHO) or similar organics with equivalent reduction potentials (Lane 2014; Lane and Martin 2012).

<figure><img src="image_2.png" /><figcaption>**Fig. 2** Photograph of the bench-top reactor containing ceramic foam within the reactor vessel. The reaction chamber is open-flow allowing for heated alkaline fluids and cool acidic fluids to be pumped into the main chamber with an outflow from the top into a collection vessel. There are several ports on the side of the reactor, which allow for addition of fluids or sampling while the reactor is in operation</figcaption></figure>
A clue might lie in the strict dependence of methanogens and acetogens on proton gradients to drive CO₂ reduction (Buckel and Thauer 2013). In methanogens, the membrane-integral energy converting hydrogenase (Ech) uses the proton-motive force to reduce ferredoxin directly, which in turn reduces CO₂, ultimately to a methyl group (Buckel and Thauer 2013). Ech could conceivably utilise proton gradients to modulate pH within the active site of the enzyme, thereby altering the reduction potential locally. Whenever protons are involved in a reduction, the reduction potential depends on pH, falling by ~59 mV per pH unit rise, according to the Nernst equation (Nicholls and Ferguson 2013). Such pH dependence is true of both H₂ and CO₂, hence at any particular pH, the reduction remains equally difficult (Fig. 1a). However, in alkaline vents, H₂ is

dissolved in hydrothermal fluids at pH 10, whereas CO₂ is dissolved in ocean waters at pH 6. This sharp difference should modulate both reduction potentials sufficiently to drive the reduction of CO₂ with H₂. If fluids of pH 6 and 10 are juxtaposed across a thin semi-conducting Fe(Ni)S barrier, it should be possible in principle to reduce CO₂ to CO, HCOO⁻ and even HCHO (Fig. 1b).

Once the energetic barrier to CO₂ reduction has been overcome, the ensuing steps of the acetyl CoA pathway are exergonic, and in methanogens and acetogens drive carbon and energy metabolism via acetyl CoA and ATP, respectively (Fuchs 2011). An abiotic equivalent of this pathway could arguably generate the reactive thioester methyl thioacetate, a simple analogue of acetyl CoA, which has been synthesised from CO and CH₃SH by Huber and Wäctershäuser (1997) using Fe(Ni)S catalysts. In modern cells, acetyl CoA can be phosphorylated without an enzyme to form acetyl phosphate, a reactive acyl phosphate that could act as an abiotic equivalent to ATP, with a higher phosphorylating potential (AcP ΔG'₀ = -43 kJ mol⁻¹, ATP ΔG'₀ADP = -31 kJ mol⁻¹), providing a source of metabolic energy for phosphorylation and condensation to form polymers such as polypeptides and RNA (de Duve 1988, 1995; Martin and Russell 2007; Lane and Martin 2012). Overall, substrate-level phosphorylation produces acetyl phosphate, as argued by Ferry and House (2006), but in this case the whole process is driven by natural proton gradients.

<table><caption>**Table 1** Table giving fluid composition for both alkaline and acidic reactor fluids</caption><thead><tr><th>Acid</th><th>Conc. (mM)</th><th>Alkali</th><th>Conc. (mM)</th></tr></thead><tbody><tr><td>FeCl<sub>2</sub></td><td>50</td><td>Na<sub>2</sub>Si<sub>3</sub>O<sub>7</sub></td><td>10</td></tr><tr><td>NaHCO<sub>3</sub></td><td>10</td><td>K<sub>2</sub>HPO<sub>4</sub></td><td>10</td></tr><tr><td>NiCl<sub>2</sub></td><td>5</td><td>Na<sub>2</sub>S</td><td>10</td></tr></tbody></table>

The fluids are made up in 2 L batches in millimolar concentrations as shown under strict anaerobic conditions. All water is degassed by bubbling nitrogen for 30 min prior to preparation

**Fig. 3** A series of photographs taken of the precipitates formed inside the reactor vessel over 4 h. **a** Initial formation of the precipitate structures. **b** After 20 min, the structures continue to form with the only disruption as they hit the surface of the fluid in the reactor. **c** After 1 h, precipitates of good structure are still forming. **d** After 4 h the precipitates become thicker around the base, probably inhibiting reduction across the barrier

<figure><img src="image_3.png" /></figure>
<figure><img src="image_4.png" /></figure>

<table><thead><tr><th>Element</th><th>Weight%</th><th>Atomic%</th></tr></thead><tbody><tr><td>O K</td><td>41.87</td><td>65.13</td></tr><tr><td>Si K</td><td>8.05</td><td>7.14</td></tr><tr><td>P K</td><td>13.66</td><td>10.98</td></tr><tr><td>Cl K</td><td>1.72</td><td>1.21</td></tr><tr><td>K K</td><td>0.45</td><td>0.29</td></tr><tr><td>Fe K</td><td>34.24</td><td>15.26</td></tr><tr><td>Totals</td><td>100.00</td><td></td></tr></tbody></table>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>
► **Fig. 4** a Powder X-ray diffraction trace showing the precipitate is amorphous in character to X-rays: there are no peaks to indicate diffraction from crystal planes. **b** Results of elemental analysis of bulk precipitate conducted by EDX analysis. **c** SEM image of the precipitate collected from the reactor at ×330 magnification. **d** SEM image of the precipitate at ×7,000 magnification. **e** TEM image of the crystalline fractions of precipitate showing the presence of long, thin tetragonal crystals. **f** ×4 magnification of previous image showing the tetragonal crystals. **g–i** TEM lattice imaging of individual crystals showing visible atomic planes in the crystals. This planar difference was measured using a Gatan Digital Micrograph. The traces show light intensity at a specific cross-section of an individual micrograph indicating the spacing between the atomic planes. Average spacings measured were **g** 0.3 nm, **h** 0.5 nm and **i** 0.5 nm

Further phosphorylation and condensation reactions are only favoured if the concentration of monomers is high. That is possible, despite the anticipated low yields of most of these reactions, because alkaline hydrothermal vents should provide a dynamic concentration mechanism known as thermophoresis (Braun and Libchaber 2002; Baaske et al. 2007). Convection currents and thermal diffusion across the interconnected microporous matrix of alkaline vents produce thermal gradients that can concentrate organic molecules in the cooler regions. In closed experimental systems, even small thermal gradients (2.3–4.4 K) concentrate large molecules, notably DNA (Reineck et al. 2010) and RNA (Mast and Braun 2010; Mast et al. 2013), while fatty acids can be concentrated sufficiently to precipitate into vesicles (Budin et al. 2009). Thermophoresis is predicted to concentrate organics in open systems such as alkaline vents, but this has not previously been tested.

At a later stage, some form of compartmentalization is also crucial for selection to act on groups of replicators (e.g. RNAs) encoding functions such as metabolism and cooperation, rather than replication speed alone, which invariably leads to the formation of ‘Spiegelman’s monsters’ (Mills et al. 1967; Branciamore et al. 2009). The natural inorganic compartments in alkaline vents could facilitate not only the concentration of organics by thermophoresis, but also the beginnings of selection for metabolism (Branciamore et al. 2009; Koonin and Martin 2005). The two processes combined could potentially drive the replication of simple organic vesicles composed of mixed amphiphiles enclosing primitive replicators within vent pores (Budin et al. 2009; Mauer and Monndard 2011). Such vesicles are capable of growth and division, while retaining RNA (Hanczyc et al. 2003; Mansy et al. 2008) and are *en route* to the known end-point, modern cells with lipid membranes.

In sum, alkaline hydrothermal vents have the potential to drive the origins of biochemistry from H₂ and CO₂ using natural proton gradients and Fe(Ni)S minerals, in a

manner remarkably analogous to the acetyl CoA pathway in methanogens and acetogens. Modern vent systems cannot replicate this chemistry, as modern oceans are aerobic, so extant vent systems lack Fe(Ni)S catalysts; they are also depleted in CO₂, starving them of carbon and diminishing natural proton gradients; and any abiotic carbon chemistry is complicated by the presence of living cells. We have therefore built a simple bench-top reactor, which operates in an anaerobic hood, to simulate pertinent conditions in alkaline hydrothermal vents and test whether such conditions could drive the origins of biochemistry.

## An Electrochemical Reactor to Simulate Alkaline Hydrothermal Vents

We report the construction and preliminary testing of a continuous, open-flow, bench-top reactor to investigate the potential of alkaline hydrothermal systems to drive the origins of biochemistry. We explore

(i) The potential of alkaline hydrothermal vents to form simple organic molecules, most importantly formaldehyde (HCHO), by reducing CO₂ with H₂ using natural proton gradients across thin, semi-conducting, inorganic barriers. Initial work reported here characterises the Fe(Ni)S precipitates and establishes sampling and detection methodologies for small organics including formate and formaldehyde.

(ii) The formation of key biochemical intermediates such as amino acids, fatty acids and sugars from methyl sulfide, CO, NH₃ and formaldehyde. Initial work reported here characterises the synthesis of sugars including ribose from HCHO via the formose reaction, which generates sugars that could be used for RNA synthesis under alkaline hydrothermal conditions.

(iii) The concentration of organic molecules within a microporous matrix, via thermophoresis under open, continuous flow conditions. Initial work reported here demonstrates substantial temperature gradients of ~50 °C across a microporous ceramic foam (diameter 9 cm), which enable the concentration of fluorescein via thermophoresis by ~5,000-fold.

These studies are preliminary, but show the scope for more sophisticated future experiments in the reactor, and reinforce the potential of alkaline hydrothermal vents as promising far-from-equilibrium electrochemical reactors for the origin of life.
# The Reactor

## Design and Fabrication

The simple bench-top reactor simulates a continuous, open-flow, alkaline hydrothermal vent (Fig. 2). The main vessel is borosilicate glass with an internal diameter 100 mm, height 100 mm and wall thickness 5 mm. Eight side ports provide for infusion of fluids or sampling. A grade-5 titanium plate is held tightly against each end of the reaction vessel, and each plate is fitted with an inlet/outlet for the reaction vessel. A drilled titanium disc flow distributor is fitted to the inlet to distribute inflow within the reactor vessel. The external titanium inlet tube is fitted with two heating elements, and thermocouples attached to this inlet tube allow temperature control of the ‘hydrothermal’ fluids. Fluids feed into the reactor through Viton™ tubing using variable speed peristaltic pumps. Flow rates are generally set between 10 and 120 mL/h, depending on requirements. The reactor is housed in an anaerobic hood under a controlled atmosphere of 98 % N₂/2 % H₂ to ensure anoxia. All solutions were prepared within the anaerobic hood using HPLC-grade water that had been deoxygenated for 24 h before use.

## Precipitation of Thin, Inorganic Barriers

To precipitate dynamic, thin-walled, semi-conducting Fe(Ni)S barriers, alkaline fluids (~pH 11, 70 °C, see Table 1 for fluid compositions) were infused into an acidic ‘ocean’ (~pH 5, 20 °C) inside the reaction vessel. Mixing of the two fluids results in the dynamic precipitation of vertical hollow tube structures composed mostly of ferrous silicates and phosphates, as well as lesser amounts of Fe(Ni)S minerals. The characteristics of these dynamic structures vary depending on the flow rates and chemical composition of the acid and alkaline fluids. To reduce CO₂ successfully, we hypothesise that hollow thin-walled structures are needed to act as semi-conducting barriers to harness the natural proton gradient from the alkaline interior to the acidic ‘ocean’ (Fig. 3). The composition of fluids and flow rate into the reactor was varied to achieve stable, thin-walled structures. The ideal alkaline flow rate required to produce suitable structures is ~50 mL/h, which allows dynamic structures to form over several hours before thickening into more solid chimney-like structures (Fig. 3).

Scanning electron microscopy (SEM) with Energy Dispersive X-ray (EDX) analysis of precipitates showed structures consisting largely of iron, phosphorus, silicon and oxygen (Fig. 4a) which are not homogeneous in composition. The precipitates are amorphous to X-rays on XRD analysis, probably because of low particle diameter. Transmission electron microscopic (TEM) inspection of some

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /><figcaption>**Fig. 5** **a** GC-MS trace showing the analysis for formate. The formate peak at 2.45 min is the propyl-ester of formate. Estimated concentration is 50 µM based on extrapolation from calibration data. **b** GC-MS trace showing analysis for formaldehyde. The formaldehyde peak at 3.8 min is the PFBOA adduct. Estimated concentration is 100 nM based on extrapolation from calibration data. **c** Graph of formaldehyde concentration over time during the course of an experimental run. After an initial increase, the concentration remains relatively constant, though repeatability of sampling, total volume and dynamic reaction environment all impact on the repeatability and consistency of results</figcaption></figure>

precipitates identified separate small crystalline structures (10 nm × 30 nm) from which lattice images were taken (Fig. 4g-i). The spacing of these planes was measured on the image using the profile tool in Digital Micrograph™. Of twelve crystals, spacings of 0.3, 0.5 and (occasionally) 0.7 nm were recorded. Generally, the low index planes are likely to be imaged. Jeong and coauthors (2008) reported the unit cell of mackinawite ([Fe,Ni]₁+x S, where x = 0 to 0.11), a tetragonal crystal structure, as a = b = 0.3 nm and c = 0.5 nm consistent with the presence of fine mackinawite crystals associated with some precipitates.
**Fig. 6** a GC–MS trace showing sugars, with the peaks for glyceraldehyde, erythrose and ribose labelled. The internal standard (IS) used was myoinositol. All enantiomers of ribose (arabinose, lyxose and xylose) were identified using known standards. b) The concentration over time of glyceraldehyde, erythrose and ribose. The area showing 0 to 5 h is with the reaction heated to 60 °C while the shaded area shows the reaction mixture *left* at ambient temperature (~20 °C) between 5 and 120 h

<figure><img src="image_19.png" /></figure>

## Formation of Organics

We now consider the reduction of CO₂ under these simulated vent conditions. Fluid samples were collected close to the precipitates within the reactor every 20 min for periods of up to 4 h, and then analysed using various different methods and techniques, including high performance liquid chromatography (HPLC) and gas chromatography–mass spectrometry (GC–MS), to identify simple organics such as formate and formaldehyde. We established the presence of formate in low micromolar quantities using a head-space GC–MS analysis by derivatization of an acidified sample with 1-propanol to form the propyl-ester (Fig. 5a). We also detected formaldehyde in samples, again by head-space GC–MS analysis, in low nanomolar quantities by derivatization of the sample with o-(2,3,4,5,6-pentafluorobenzyl)-hydroxylamine (Fig. 5b). While these preliminary findings are encouraging, the formation of formaldehyde, especially, is variable and inconsistent between runs, and these methods are still being optimised.

## Formation of Ribose Via the Formose Reaction Under Alkaline Hydrothermal Conditions

The production of formaldehyde in the reactor is a proof of principle. Introducing formaldehyde at higher concentration should drive the synthesis of various sugars under simulated vent conditions via the formose reaction. Using the method of Kopetzki and Antonietti (2011), which is consistent with internal alkaline vent conditions, formaldehyde (0.5 M) was heated in alkaline fluids (initially pH 12) at 60 °C for 5 h. A variety of sugars up to C7 (septoses) were identified within the reaction fluids by GC–MS analysis (Fig. 6) after reduction with sodium borohydride and derivatization with acetic anhydride and acetic acid. Sugars identified by this process include ribose and deoxyribose, albeit at low yield (~0.06 % of formaldehyde was converted to ribose). Ribose is stable under these alkaline conditions for at least several hours, with any breakdown balanced by new synthesis. This is significant as it allows ribose to be formed under alkaline conditions, at low yields, and then potentially concentrated
**Fig. 7** a SEM image of the internal structure of the ceramic foam at ×70 magnification. The foam has a microporous, highly permeable structure with interconnected cavities ~100 µm in diameter. b SEM image of the foam structure at ×3,500 magnification showing ~10 µm cavities within the foam structure and also sub-micron holes inside the ceramic struts of the foam. c SEM image of a foam strut at ×6,500 magnification showing in greater detail the sub-micron cavities in a foam strut

<figure><img src="image_20.png" /></figure>

via thermophoresis in cooler, more neutral regions of the vent, discussed below.

## Vent Matrix

As alkaline vents are chemically varied (Lost City is mainly calcium carbonate, whereas Strýtan is constructed of saponite clay), it is difficult to establish which ions or minerals are important within the vent for promotion of prebiotic reactions, especially as early mineral structures would have been very different in the absence of oxygen and in mildly acidic oceans. We are therefore using life itself as our guide, as discussed above. A porous, ceramic foam matrix allows full control of the chemical environment. The microporous foam is constructed from alumina, which is chemically inert and has minimal catalytic properties. Depending on pH, alumina allows adsorption onto its surface and can be doped with ions or minerals relevant to early biochemistry, notably Fe(Ni)S, Mg²⁺, Ca²⁺ and Mo⁴⁺, as required to investigate their effects on vent chemistry.

Foams were constructed to fit the reactor vessel, with diameter of 100 mm and height of 100 mm. The foam matrix has a ~4 % density with ~81.5 % continuous or ‘open’ porosity, meaning it has a large internal volume of interconnected pores, as in Lost City, although the morphology of our synthetic foam differs from the elongated vertical channels characteristic of the labyrinthine structure of Lost City. The ceramic foam is very hydrophilic, with water immediately and rapidly drawn in by capillary action. Analysis by SEM shows that the foam consists of a highly interconnected, porous network (100–500-µm diameters) with thin walls (1–6 µm in diameter) which themselves are highly porous with <1 µm pores in the alumina substrate (see Fig. 7 and inset).

## Temperature Gradients

Ceramic foam filling the full reactor vessel was submerged in ‘ocean’ fluid at 15 °C. ‘Hydrothermal’ fluid heated to ~75 °C was pumped into the foam at ~30 mL/h via a flow distributor. Cooling fluid was pumped into the foam in
<figure><img src="image_21.png" /><figcaption>**Fig. 8** Temperature profile inside a foam exposed to vent-like conditions. The temperature graphs present axial profiles based on the division of the foam as shown (1–5, 2–6, 3–7, 4–8) in the circle diagram (inset). Fluid at 70 °C was flowing into the foam at a rate of 15 mL/h with a coolant fluid at 20 °C being pumped into the reactor</figcaption></figure>

at a rate of 120 mL/h. The reactor was left for a period of 1 h prior to temperature readings being taken. The warmer temperatures (*red*) are observed in the lower central regions with cooler temperatures (*blue*) in the upper outer regions of the foam. The photo (inset) shows the size of the foam before the temperature profile was taken

the reactor vessel via 8 side ports, at ~180 mL/h to help maintain the maximum temperature gradient, while overflow was pumped out from the top. The reactor was allowed to equilibrate for 45 min before readings were taken. Temperature readings were taken from top to bottom, outside to centre along 8 axes shown in figure inset, at 10-mm depths. The temperature readings and locations are shown in Fig. 8.

The temperature profiles across the foam are as anticipated. The warmest region is at the point of inflow (centre, bottom) with cooler regions to the outside edges (where cooling fluids enter via the side ports) and the top of the foam. Fluid convection is also observed within the foam around the midpoint based on analysis of the axial temperature profiles. Temperature gradients of ≤0.5 °C mm⁻¹ (i.e. up to 50 °C across the whole foam) are stable for several hours.

## Organic Concentration

We investigated whether thermophoresis can operate within the ceramic foam under open-flow conditions, using

fluorescein and quinine. These compounds are UV fluorescent, so their distribution can be observed under a UV lamp, while allowing for semi-quantitative analysis of concentration. Initial testing of foam pieces soaked in standard concentrations of each indicator showed that fluorescein gives the strongest fluorescent response in the foam, but saturates at ~30 µM concentration, while quinine saturates at ~50 µM (Fig. 7). However, fluorescein is visible at concentrations >100 µM without excitation. We therefore investigated the potential thermophoretic concentration of 1 µM fluorescein.

Using standard temperature/flow conditions for thermal gradients (15 °C cooling fluids, 180 mL/h flow and 75 °C fluorescein solution, 30 mL/h flow), fluorescein (1 µM solution) was infused from the base over 7 h. Visible inspection of foams showed large regions towards the exterior (5–10-mm diameter) with intense concentration (ca. 200 and 400 µM, respectively, Fig. 9b). Thin slices of damp foam show large numbers (50–100) of 3–5-mm diameter spots in cooler regions (ca. 20–50 µM concentration). The background concentration observed is much lower than the infusion concentration. The infusion enters
<figure><img src="image_22.png" /><figcaption>**Fig. 9** a Fluorescent micrograph of foam soaked in a 50 µM fluorescein solution. Fluorescein enters the struts of the foam via the sub-micron pores in the structure and remains within the foam struts, not in the cavities as originally assumed. b Sections of the foam exposed to UV light. The foam has been infused with 0.1 µM fluorescein solution for a period of 4 h under vent conditions. The *bright blue areas* in the photos are areas of fluorescein concentration, estimated to be between ×2,500 and ×5,000 the concentration of inflow fluids (0.1 µM)</figcaption></figure>

at 0.1 µM, not 1 µM, due to dilution by bulk ‘ocean’ fluid. It then accumulates slightly within the foam, giving a background concentration of 0.45–0.55 µM, before dramatically concentrating in the cooler regions.

These patterns suggest thermophoretic concentration. As expected, the regions of the foam with the highest

concentrations are in the cooler zones (outer edges), while the warmest region close to the central base shows the weakest concentration, despite the direct entry of fluorescein into the foam at this point. Control runs without coolant or without heating failed to show any significant thermal gradients, or high concentrations of fluorescein. Fast coolant flow through the side ports precludes the external flow of fluorescein solution around the foam (which is in any case very hydrophilic, as noted). Thermophoresis of fluorescein followed by adsorption onto alumina struts in cooler regions (Fig. 9) shows that fluorescein can be concentrated up to 5,000-fold. Such high concentration may also be achievable with more prebiotically relevant molecules such as amino acids and nucleotides.

## Conclusions

We have constructed a simple bench-top reactor to investigate the possible origins of biochemistry in alkaline hydrothermal vent systems. Within the reactor vessel, thin-walled, inorganic structures precipitated, containing catalytic Fe(Ni)S microcrystals. When transected by pH gradients, these precipitates appear to be capable of reducing CO₂ to form low yields of simple 1C organics including formate and formaldehyde. Starting with formaldehyde under mild alkaline hydrothermal conditions (60 °C, pH 11–12), a range of sugars, including ribose and deoxyribose, could be formed via the formose reaction. The concentration of these sugars is stable over several hours. While produced at low yield, organics can theoretically be concentrated by thermophoresis in a microporous matrix. We produced a simulated vent matrix from alumina foam with a highly interconnected microporous network similar to those seen in the Lost City vent system. Temperature gradients of 50 °C across this matrix can concentrate fluorescein up to ~5,000-fold under open-flow vent conditions. The studies presented here are a preliminary proof of concept, and still require optimising in many respects; they do, however, provide empirical evidence that simple organics can be generated and concentrated under mild alkaline hydrothermal conditions from H₂ and CO₂ using Fe(Ni)S catalysts transected by natural proton gradients in microporous matrices.

**Acknowledgments** We would like to acknowledge and thank the Leverhulme Trust (Grant Number: RPG-425) for their funding of this research. We thank Prof. Mike Russell, Prof. Bill Martin, Dr. Terry Kee and Dr. Laurie Barge for their helpful discussions and insight on prebiotic chemistry and geochemistry. We also thank Prof. Don Braben and Prof. David Price (VP Research) for support through the UCL Provost’s Venture Research Fellowship, and the UCL Research Frontiers programme. Thanks also to Dr. Vernon Skinner and Dr.
Kersti Karu for their assistance and advice with the analytical methods reported in this article.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.

## References

Amend JP, McCollom TM (2009) Energetics of biomolecule synthesis on early Earth. In: Zaikowski L (ed) Chemical evolution II: From the origins of life to modern society. American Chemical Society, Washington, DC, pp 63–94

Amend JP, LaRowe DE, McCollom TM, Shock EL (2013) The energetics of organic synthesis inside and outside the cell. Philos Trans R Soc B 368:20120255. doi:10.1098/rstb.2012.0255

Anbar AD, Holland HD (1992) The photochemistry of manganese and the origin of banded iron formations. Geochim Cosmochim Acta 56:2595–2603

Arndt N, Nisbet E (2012) Processes on the young earth and the habitats of early life. Annu Rev Earth Planet Sci 40:521–549

Baaske P, Weinert FM, Duhr S et al (2007) Extreme accumulation of nucleotides in simulated hydrothermal pore systems. Proc Natl Acad Sci USA 104:9346–9351

Bach W, Paulick H, Garrido CJ et al (2006) Unraveling the sequence of serpentinization reactions: petrography, mineral chemistry, and petrophysics of serpentinites from MAR 15 N (ODP Leg 209, Site 1274). Geophy Res Lett 33:L13306

Baymann F, Lebrun E, Brugna M et al (2003) The redox protein construction kit: pre-last universal common ancestor evolution of energy conserving enzymes. Philos Trans R Soc Lond B 358:267–274

Bekker A, Holland HD, Wang PL, Rumble D, Stein HJ, Hannah JL, Coetzee LL, Beukes NJ (2004) Dating the rise of atmospheric oxygen. Nature 427:117–120

Boyce A, Coleman M, Russell M (1983) Formation of fossil hydrothermal chimneys and mounds from Silvermines, Ireland. Nature 306:545–550

Bradley AS, Hayes JM, Summons RE (2009) Extraordinary C-13 enrichment of diether lipids at the Lost City Hydrothermal Field indicates a carbon-limited ecosystem. Geochim Cosmochim Acta 73:102–118

Branciamore S, Gallori E, Szathmáry E, Czárán T (2009) The origin of life: chemical evolution of a metabolic system in a mineral honeycomb? J Mol Evol 69:458–469

Braun D, Libchaber A (2002) Trapping of DNA by thermophoretic depletion and convection. Phys Rev Lett 89:188103

Buckel W, Thauer RK (2013) Energy conservation via electron bifurcating ferredoxin reduction and proton/Na + translocating ferredoxin oxidation. Biochim Biophys Acta 1827:94–113

Budin I, Bruckner RJ, Szostak JW (2009) Formation of protocell-like vesicles in a thermal diffusion column. J Am Chem Soc 131:9628–9629

Cody GD (2004) Transition metal sulfides and the origin of metabolism. Ann Rev Earth Planet Sci 32:569–599

Cox CJ, Foster PG, Hirt RP et al (2008) The archaebacterial origin of eukaryotes. Proc Natl Acad Sci USA 105:20356–20361

Dagan T, Martin W (2007) Ancestral genome sizes specify the minimum rate of lateral gene transfer during prokaryote evolution. Proc Natl Acad Sci USA 104:870–875

Dagan T, Roettger M, Bryant D, Martin W (2010) Genome networks root the tree of life between prokaryotic domains. Genome Biol Evol 2:379–392

de Duve C (1988) Did God make RNA? Nature 336:209–210

de Duve C (1995) Vital dust: life as a cosmic imperative. Basic Books, New York

de Duve C (2005) Singularities: landmarks on the pathways of life. Oxford University Press, Oxford

De Leeuw NH, Catlow CR, King HE, Putnis A, Muralidharan K, Deymier P, Stimpfl M, Drake MJ (2010) Where on Earth has our water come from? Chem Commun 46:8923–8925

Doolittle WF (1999) Phylogenetic classification and the universal tree. Science 284:2124–2128

Doolittle WF, Bapteste E (2007) Pattern pluralism and the tree of life hypothesis. Proc Natl Acad Sci USA 104:2043–2049

Ducluzeau A-L, Schoepp-Cothenet B, Baymann F, Russell MJ, Nitschke W (2014) Free energy conversion in the LUCA: Quo vadis? Biochim Biophys Acta Bioenergetics 1837:982–988

Eck RV, Dayhoff MO (1966) Evolution of the structure of ferredoxin based on living relics of primitive amino acids sequences. Science 152:363–366

Edgell DR, Doolittle WF (1997) Archaea and the origin(s) of DNA replication proteins. Cell 89:995–998

Ferry JG, House CH (2006) The step-wise evolution of early life driven by energy conservation. Mol Biol Evol 23:1286–1292

Fuchs G (2011) Alternative pathways of carbon dioxide fixation: insights into the early evolution of life? Annu Rev Microbiol 65:631–658

Fuchs G, Stupperich E (1985) Evolution of autotrophic CO₂ fixation. In: Schleifer KH, Stackebrandt E (ed) Evolution of Prokaryotes, FEMS Symposium No. 29. Academic Press, London. pp 235–251

Fyfe WS (1994) The water inventory of the Earth: fluids and tectonics. Geol Soc Lond (Spec Publ) 78:1–7

Hanczyc M, Fujikawa S, Szostak J (2003) Experimental models of primitive cellular compartments: encapsulation, growth, and division. Science 302:618–622

Harel A, Bromberg Y, Falkowski PG, Bhattacharya D (2014) Evolutionary history of redox metal-binding domains across the tree of life. Proc Natl Acad Sci USA 111:7042–7047. doi:10.1073/pnas.1403676111

Huber C, Wäctershäuser G (1997) Activated acetic acid by carbon fixation on (Fe, Ni)S under primordial conditions. Science 276:245–247

Jaffrés JBD, Shields GA, Wallmann K (2007) The oxygen isotope evolution of seawater: a critical review of a long-standing controversy and an improved geological water cycle model for the past 3.4 billion years. Earth Sci Rev 83:83–122

Kaster A-K, Moll J, Parey K, Thauer RK (2011) Coupling of ferredoxin and heterodisulfide reduction via electron bifurcation in hydrogenotrophic methanogenic Archaea. Proc Natl Acad Sci USA 108:2981–2986

Kasting JF (2013) What caused the rise in atmospheric O₂? Chem Geol 362:13–25

Kelley DS, Karson JA, Blackman DK, Früh-Green GL, Butterfield DA, Lilley MD, Olson EJ, Schrenk MO, Roe KK, Lebon GT, Rivizzigno P (2001) An off-axis hydrothermal vent field near the Mid-Atlantic Ridge at 30 degrees N. Nature 412:145–149

Kelley DS, Karson JA, Früh-Green GL et al (2005) A serpentinite-hosted submarine ecosystem: the Lost City Hydrothermal Field. Science 307:1428–1434

Kim JD, Senn S, Harel A, Jelen BI, Falkowski PG (2013) Discovering the electronic circuit diagram of life: structural relationships among transition metal binding sites in oxidoreductases. Philos Trans R Soc Lond B 368:20120257

Koga Y, Kyuragi T, Nishihara M, Sone N (1998) Did archaeal and bacterial cells arise independently from noncellular precursors? A hypothesis stating that the advent of membrane phospholipid with enantiomeric glycerophosphate backbones caused the separation of the two lines of descent. J Mol Evol 46:54–63
Koonin EV, Martin W (2005) On the origin of genomes and cells within inorganic compartments. *Trends Genet* 21:647–654

Kopetzki D, Antonietti M (2011) Hydrothermal formose reaction. *New J Chem* 35:1787–1794

Kopetzki D, Jeong HY, Lee JH, Hayes KF (2008) Characterization of synthetic nanocrystalline mackinawite: crystal structure, particle size, and specific surface area. *Geochim Cosmochim Acta* 72:493–505

Lane N (2014) Bioenergetic constraints on the evolution of complex life. *Cold Spring Harb Perspect Biol.* doi:10.1101/cshperspect. a015982

Lane N (2015) The vital question: energy, evolution and the origin of complex life. *WW Norton, New York*

Lane N, Martin WF (2012) The origin of membrane bioenergetics. *Cell* 151:1406–1416

Lane N, Allen JF, Martin W (2010) How did LUCA make a living? *Chemiosmosis in the origin of life. BioEssays* 32:271–280

Leipe DD, Aravind L, Koonin EV (1999) Did DNA replication evolve twice independently? *Nucleic Acids Res* 27:3389–3401

Ljungdahl LG (2009) A life with acetogens, thermophiles, and cellulolytic anaerobes. *Annu Rev Microbiol* 63:1–25

Lombard J, López-García P, Moreira D (2012) The early evolution of lipid membranes and the three domains of life. *Nat Rev Microbiol* 10:507–515

Ludwig KA, Kelley DS, Shen C, Cheng H, Edwards RL (2005) U/Th Geochronology of carbonate chimneys at the Lost City Hydrothermal Field. *Am Geophys Union, Fall Meeting abstract #V51B-1487*

Maden BEH (1995) No soup for starters? Autotrophy and the origins of metabolism. *Trends Biochem Sci* 20:337–341

Maden BEH (2000) Tetrahydrofolate and tetrahydromethanopterin compared: functionally distinct carriers in C1 metabolism. *Biochem J* 350:609–629

Major TA, Burd H, Whitman WB (2004) Abundance of 4Fe-4S motifs in the genomes of methanogens and other prokaryotes. *FEMS Microbiol Lett* 239:117–123

Mansy SS, Schrum JP, Krishnamurthy M, Tobé S, Treco DA et al (2008) Template-directed synthesis of a genetic polymer in a model protocell. *Nature* 454:122–125

Martin WF (2012) Hydrogen, metals, bifurcating electrons, and proton gradients: the early evolution of biological energy conservation. *FEBS Lett* 586:485–493

Martin W, Russell MJ (2003) On the origins of cells: a hypothesis for the evolutionary transitions from abiotic geochemistry to chemoautotrophic prokaryotes, and from prokaryotes to nucleated cells. *Philos Trans R Soc Lond B Biol Sci* 358:59–83

Martin W, Russell MJ (2007) On the origin of biochemistry at an alkaline hydrothermal vent. *Philos Trans R Soc Lond B* 367:1887–1925

Martin W, Baross J, Kelley D, Russell MJ (2008) Hydrothermal vents and the origin of life. *Nat Rev Microbiol* 6:805–814

Martin WF, Sousa FL, Lane N (2014) Energy at life's origin. *Science* 344:1092–1093

Mast CB, Braun D (2010) Thermal trap for DNA replication. *Phys Rev Lett* 104(188102):1–4

Mast CB, Schink S, Gerland U, Braun D (2013) Escalation of polymerization in a thermal gradient. *Proc Natl Acad Sci USA* 110:8030–8035

Mauer SE, Monndard PA (2011) Primitive membrane formation, characteristics and roles in the emergent properties of a protocell. *Entropy* 13:466–484

McCollom TM, Bach W (2009) Thermodynamic constraints on hydrogen generation during serpentinization of ultramafic rocks. *Geochim Cosmochim Acta* 73:856–875

Mills DR, Peterson RL, Spiegelman S (1967) An extracellular Darwinian experiment with a self-duplicating nucleic acid molecule. *Proc Natl Acad Sci USA* 58:217–224

Mitchell P (1959) The origin of life and the formation and organizing functions of natural membranes. In: Oparin AI et al (eds) *Proceedings of the first international symposium on the origin of life on the earth*

Mitchell P (1961) Coupling of phosphorylation to electron and hydrogen transfer by a chemi-osmotic type of mechanism. *Nature* 191:144–148

Mitchell P (1966) Chemiosmotic coupling in oxidative and photo-synthetic phosphorylation. *Biol Rev* 41:445–501

Morowitz HJ, Kostelnik JD, Yang J, Cody GD (2000) The origin of intermediary metabolism. *Proc Natl Acad Sci USA* 97:7704–7708

Nelson-Sathi S, Sousa FL, Roettger M, Lozada-Chávez N, Thiergart T, Janssen A, Bryant D, Landan G, Schönheit P, Siebers B, McInerney JO, Martin WF (2014) Origins of major archaeal clades correspond to gene acquisitions from bacteria. *Nature*. doi:10.1038/nature13805.

Nicholls DG, Ferguson SJ (2013) Bioenergetics, 4th edn. Academic Press, London

Nitschke W, Russell MJ (2009) Hydrothermal focusing of chemical and chemiosmotic energy, supported by delivery of catalytic Fe, Ni, Mo Co, S and Se forced life to emerge. *J Mol Evol* 69:481–496

Orgel LE (2008) The implausibility of metabolic cycles on the prebiotic earth. *PLoS Biol* 6(1):e18. doi:10.1371/journal.pbio.0060018

Peretó J, López-García P, Moreira D (2004) Ancestral lipid biosynthesis and early membrane evolution. *Trends Biochem Sci* 29:469–477

Pinti D (2005) The origin and evolution of the oceans. *Lect Astrobiol* 1:83–112

Poehlein A, Schmidt S, Kaster A-K, Goenrich M, Vollmers J, Thürmer A, Bertsch J, Schuchmann K, Voigt B, Hecker M, Daniel D, Thauer RK, Gottschalk G, Müller V (2012) An ancient pathway combining carbon dioxide fixation with the generation and utilization of a sodium ion gradient for ATP synthesis. *PLoS ONE* 7:e33439

Proskurowski G, Lilley MD, Seewald JS et al (2008) Abiogenic hydrocarbon production at Lost City Hydrothermal Field. *Science* 319:604–607

Ragsdale SW, Kumar M (1996) Nickel-containing carbon monoxide dehydrogenase/acetyl-CoA synthase. *Chem Rev* 96:2515–2539

Ragsdale SW, Pierce E (2008) Acetogenesis and the Wood-Ljungdahl pathway of CO₂ fixation. *Biochim Biophys Acta* 1784:1873–1898

Reineck P, Wienken CJ, Braun D (2010) Thermophoresis of single stranded DNA. *Electrophoresis* 31:279–286

Russell MJ, Arndt NT (2005) Geodynamic and metabolic cycles in the Hadean. *Biogeosciences* 2:97–111

Russell MJ, Hall AJ (1997) The emergence of life from iron monosulphide bubbles at a submarine hydrothermal redox and pH front. *J Geol Soc Lond* 154:377–402

Russell MJ, Martin W (2004) The rocky roots of the acetyl-CoA pathway. *Trends Biochem Sci* 29:358–363

Russell MJ, Hall AJ, Cairns-Smith AG, Braterman PS (1988) Submarine hot springs and the origin of life. *Nature* 336:117

Russell MJ, Hall AJ, Turner D (1989) In-vitro growth of ironsulphide chimneys: possible culture chambers for origin-of-life experiments. *Terra Nova* 1:238–241

Russell MJ, Daniel RM, Hall A (1993) On the emergence of life via catalytic iron-sulphide membranes. *Terra Nova* 5:343–347

Russell MJ, Daniel RM, Hall AJ, Sherringham J (1994) A hydrothermally precipitated catalytic iron sulphide membrane as a first step toward life. *J Mol Evol* 39:231–243

Russell MJ, Hall AJ, Martin W (2010) Serpentinization as a source of energy at the origin of life. *Geobiology* 8:355–371
Say RF, Fuchs G (2010) Fructose 1,6-bisphosphate aldolase/phosphatase may be an ancestral gluconeogenic enzyme. Nature 464:1077–1081

Schrenk MO, Brazelton WJ, Lang SQ (2013) Serpentinization, carbon and deep life. *Rev Mineral Geochem* 75:575–608

Shields GA, Kasting JF (2007) Evidence for hot early oceans? Nature 447:E1–E2

Shock E, Canovas P (2010) The potential for abiotic organic synthesis and biosynthesis at seafloor hydrothermal systems. *Geofluids* 10:161–192

Shock EL, McCollom TM, Schulte MD (1998) The emergence of metabolism from within hydrothermal systems. In: Wiegel J, Adams MWW (eds) Thermophiles: the keys to molecular evolution and the origin of life. Taylor & Francis, London, pp 59–76

Sleep NH (2010) The Hadean-Archaean environment. Cold Spring Harb Perspect Biol 2(6):a002527

Sleep NH, Meibom A, Fridriksson T et al (2004) H₂-rich fluids from serpentization: geochemical and biotic implications. Proc Natl Acad Sci USA 101:12818–12823

Sojo V, Pomiankowski A, Lane N (2014) A bioenergetic basis for membrane divergence in archaea and bacteria. PLoS Biol 12(8):e1001926

Sousa FL, Martin WF (2014) Biochemical fossils of the ancient transition from geoenergetics to bioenergetics in prokaryotic one carbon compound metabolism. *Biochim Biophys Acta Bioenergetics* 1837:964–981

Sousa FL, Thiergart T, Landan G, Nelson-Sathi S, Pereira IAC et al (2013) Early bioenergetic evolution. *Philos Trans R Soc Lond B* 368:1–30

Stetter KO (2006) Hyperthermophiles in the history of life. *Philos Trans R Soc B* 361:1837–1843

Thauer RK, Kaster A-K, Seedorf H, Buckel W, Hedderich R (2007) Methanogenic archaea: ecologically relevant differences in energy conservation. *Nat Rev Microbiol* 6:579–591

Williams TA, Foster PG, Cox CJ, Embley TM (2013) An archaeal origin of eukaryotes supports only two primary domains of life. *Nature* 504:231–236

Yamaguchi A, Yamamoto M, Takai K, Ishii T, Hashimoto K, Nakamura R (2014) Electrochemical CO₂ reduction by Ni-containing iron sulfides: how is CO₂ electrochemically reduced at bisulfide-bearing deep sea hydrothermal precipitates? *Electrochim Acta* 141:311–318

Zahnle K, Arndt N, Cockell C, Halliday A, Nisbet E, Selsis F, Sleep NH (2007) Emergence of a habitable planet. *Planet Space Sci* Rev 129:35–78