The origin of carbon isotope vital effects in coccolith calcite.

Calcite microfossils are widely used to study climate and oceanography in Earth's geological past. Coccoliths, readily preserved calcite plates produced by a group of single-celled surface-ocean dwelling algae called coccolithophores, have formed a significant fraction of marine sediments since the Late Triassic. However, unlike the shells of foraminifera, their zooplankton counterparts, coccoliths remain underused in palaeo-reconstructions. Precipitated in an intracellular chemical and isotopic microenvironment, coccolith calcite exhibits large and enigmatic departures from the isotopic composition of abiogenic calcite, known as vital effects. Here we show that the calcification to carbon fixation ratio determines whether coccolith calcite is isotopically heavier or lighter than abiogenic calcite, and that the size of the deviation is determined by the degree of carbon utilization. We discuss the theoretical potential for, and current limitations of, coccolith-based CO2 paleobarometry, that may eventually facilitate use of the ubiquitous and geologically extensive sedimentary archive.

The taxonomic delimitation of coccoliths into the isotopically light and heavy groups, where the calcite is respectively depleted and enriched in the heavy isotope of carbon compared with abiogenic calcite, is well established12345. However, there have been relatively few studies of carbon isotopes in coccolith calcite6, and fewer still where experiments have been carried out in vivo under closely monitored conditions6789. Although so-called vital effects are known to vary with growth rate4, and CO2 availability7, the reasons for these isotopic departures from equilibrium remain poorly understood. The relatively sparse evidence from experimental, cruise, sediment trap10 and down-core12311 studies show that members of the light group such as Coccolithus pelagicus and Calcidiscus leptoporus, tend to be large, heavily calcifying, and have a relatively slow division rate. Contrastingly, members of the heavy group such as Emiliania huxleyi are generally small, lightly calcifying and with a rapid division rate467. The co-variance of parameters across species means that their influence upon coccolith calcite vital effects cannot be decoupled without an explicit mechanistic understanding and a quantitative model.

In coccolithophores, dissolved inorganic carbon (DIC) exists in mutually exchanging intracellular reservoirs. The intracellular process responsible for the largest isotopic fractionation in these organisms is the conversion of inorganic carbon dioxide (CO2) to organic matter (carbon fixation), which has isotopic implications for carbon pools elsewhere within the cell, including, in coccolithophores, the site of calcification. The rate-limiting step of carbon fixation, catalysed by the enzyme ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) proceeds far more rapidly for CO2 containing the light isotope of carbon. Organic matter is therefore significantly depleted in the heavy isotope of carbon relative to the CO2 source121314. In addition to the kinetic isotopic discrimination effective at the enzymatic level , the overall carbon isotopic offset between the extracellular carbon substrate and the organic matter , is a function of a number of parameters that affect the isotopic composition of CO2 immediately in the vicinity of the RuBisCO enzyme15. These parameters include growth rate (μ), cell radius (r) and ambient CO2 concentration ([CO2]). So far, experiments and models attempting to understand these processes have dominantly focused on relating to [CO2] in phytoplankton (Table 1)1617181920212223242526, and the empirically established relationships have been extensively applied to ancient sediments to generate reconstructions of pCO2atm (refs 27, 28, 29, 30, 31). To quantitatively explain the origin of carbon isotopic vital effects in coccolith calcite, a holistic understanding of carbon isotopes within a coccolithophore cell is required, which considers the bidirectional isotopic effect of intracellular calcification and carbon fixation, and is consistent with the available biological evidence.

Table 1

Assumptions and conclusions of significant recent works modelling carbon fluxes in single-celled phytoplankton.

Study	Model set-up and assumptions	Conclusions
Holtz et al.4748 (E. huxleyi)	4 compartments (PY=pyrenoid, CP=chloroplast, CV=coccolith vesicle, CY=cytosol)	is used at low [CO2]
	Carbonate chemistry consists of , , CO2 and H+.	Intracellular pH gradients allow concentration of CO2 around RuBisCO without up-gradient movement of carbon.
	Hypothesized Ca2+/ CY-to-CV transporter coupled to an H+-ATPase, with no leakage of from CV	pHs: PY=5.0, CY=7.0, CP=8.0, CV=8.3–8.6
	Hypothesized upregulation of down-gradient flux into cell with decreased [CO2] in PY.	A net efflux of CO2 is not necessary to remove δ13C from cell
	Passive CO2 and fluxes. No flux from CP to PY.
	CA assumed in CP and PY but not in CY and CV
	Isotope model consists of 2 compartments (CY and CP) and does not consider isotopes of calcite.
	Membrane permeabilities to CO2 and assumed, but highly heterogeneous; different for all compartments.
Bolton & Stoll,2 (Coccolithophores)	3 compartments (CP, CV and CY)	active transport to CP increases at low [CO2], at the expense of transport to CV
	Carbonate chemistry consists of and CO2.	This effect is greatest in large cells.
	fluxes are all active and independent of []	Difference in vital effects (δ13C calcite—δ13Cmedium) between small and large cells greatest at low [CO2]
	Passive CO2 fluxes.
	Membrane permeabilities and CA activities assumed from Hopkinson et al.46.
	assumed (−27‰).
Hopkinson et al.46 (Diatoms)	1, 2 & 3 compartments (PY, CP, CY)	Membranes are highly permeable to CO2 (1.5 × 10−4–5.6 × 10−4 m s−1)
	Carbonate chemistry consists of and CO2.	Membranes are highly impermeable to (2.5 × 10−8–2.9 × 10−7 m s−1)
	Used 18O labelled DIC to track temporal evolution of carbonate system.	δ13C org is a function of passive diffusion of CO2, active movement of , kinetic fractionation factors associated with CA-catalysed hydration and dehydration, and the kinetic isotopic fractionation associated with RuBisCO.
	Passive CO2 fluxes.
	Passive and active fluxes.
	assumed (−29‰).
Schulz et al.,26 (E. huxleyi)	2 compartments (CP and CY)	Carbon concentrating mechanism relies upon active (ATP driven) uptake of CO2 and
	Carbonate chemistry consists of and CO2.	Reduction in with increased uptake into CP
	Active uptake of and CO2 independent of concentration	is larger when there is a greater degree of intracellular carbon leakage from the chloroplast.
	Passive CO2 fluxes (membrane permeability to CO2=1.8 × 10−5 from76- green algae)
	No efflux of HCO3 −
	assumed (−29‰).
Cassar et al.,25 (Diatom—Phaeodactylum tricornutum)	2 compartments	is not a unique function of
	Active and diffusive uptake of CO2	depends on leakiness of CP
	No uptake
	assumed (−27‰).
	Inferred fluxes based on an energy minimization approach.
Keller & Morel,24 (General phytoplankton)	1 compartment	Downward curvature of against is consistent with active uptake contrary to32
	No diffusion or efflux	is a poor indicator of carbon substrate
	Active uptake scales with growth rate
	inferred from model.

For earlier work and the evolving appreciation of the importance of cell size, shape and growth rate see Laws et al.15 and references therein.

To understand the origin of vital effects in coccolith calcite, we derived a cellular model of carbon isotopic fluxes, whose mechanistic underpinning is based on a consideration of the available literature. The model presented here is calibrated with new data from in vivo experiments. We found that the so-called isotopically heavy and light group coccoliths, characterized by respectively positive and negative vital effects, are produced by coccolithophores with respectively low and high calcification to photosynthesis ratios. Vital effects are dominantly the result of competing Rayleigh-type fractionation processes and fluxes between intracellular compartments that alter the isotopic composition of carbon at the site of calcification. At the equilibrium limit, and in non-calcifying cells, our model collapses to consistency with the classic literature describing carbon isotopes in phytoplankton1617182123242532. Our conclusions are compatible with observed trends in oxygen isotopes. Lastly, we discuss the potential for coccolith-based CO2 paleobarometry.

Results

Model description

We model the coccolithophore cell as consisting of three physical compartments: a chloroplast (where carbon fixation occurs) and a coccolith vesicle (a golgi-body derived compartment where calcification occurs) each contained within a third compartment, the cytosol (Fig. 1). Carbon within each compartment exists as CO2, and . As is assumed to be the substrate for calcification (justified at the rates considered here), and as membranes are assumed to be impermeable to , this latter species does not feature in flux balance equations at steady state. The dynamics of the system are described entirely by CO2 and . In every compartment, CO2 and are constantly being interconverted (Supplementary Fig. 1). Compartmental shapes and sizes relative to the cell are approximated from transition electron microscopy and assumed isomorphic across all strains and experiments. Values for the H+ concentration in each compartment are taken from Anning et al.33 The uncatalysed reacto-diffusive supply rate of CO2 alone to the surface of algal cells larger than a few microns is not high enough to account for the observed rates of photosynthesis, given the low catalytic turnover and poor substrate affinity of RuBisCO14343536. The transmembrane passive diffusive supply of CO2 is therefore supplemented by bicarbonate (; the ionically charged, isotopically heavier and vastly more abundant species of DIC)83738.

Figure 1

Cellular compartmental configuration and fluxes as modelled.

Chemical and isotopic equilibrium is assumed in the external medium, but neither is assumed inside the cell. Labelling of fluxes and compartments follows the convention of Hopkinson et al.46 and Bolton and Stoll (ref. 2). Nomenclature is as follows: C, B refer to Carbon dioxide (CO2) and Bicarbonate () respectively. Subscripts e, i, x and v refer to compartments: external, cytosol, choroplast and coccolith vesicle respectively (NB: subscript cell refers to outer cell membrane, whilst subscript i refers to cytosol. Vi=Vcell−Vx−Vv). FΘab represents the flux of carbon species Θ, from compartment a to compartment b in units of mols−1. FΘΦa represents the rate of conversion of carbon species Θ, to carbon species Φ, in compartment a, in units of mols−1. Θa represents the concentration of carbon species, Θ, in compartment a in units of molm−3. PΘa represents the permeability of compartment a membrane to carbon species Θ in units of ms−1. Va, SAa represent the volume and surface area of compartment, a in units of m3 and m2 respectively. fa, sa represents the scale and shape factor of compartment a, for inferring Va and SAa ([Dimensionless]). See ‘Methods' section for full model derivation and parameter values.

Membrane permeabilities to CO2 and to have not been measured for coccolithophores, and the isotopic fractionation of carbon by RuBisCO is known only for a single species of coccolithophore39. We therefore assume these parameters to be described by a set of functions whose coefficients are constant across species and treatments, and which can be constrained by empirical data. The model permits the bidirectional movement of and CO2 across three membranes. Sinks of carbon are calcification, carbon fixation and leakage from the cell. For a full derivation of the model and complete description of the assumptions on which is it based see ‘Methods' section, and Fig. 1.

In coccolithophores, bicarbonate is thought to enter the cell via a co-transport mechanism facilitated by proteins belonging to the SLC4 family of bicarbonate transporters4041. These proteins, which are upregulated at low [CO2]42, couple the transmembrane movement of bicarbonate to an anion moving in the opposite direction (antiport), or a cation moving in the same direction (symport)—this is secondary active transport of bicarbonate as the cotransported ion gradient must ultimately be maintained with ATP43. As this type of bicarbonate transport is electroneutral, the force driving the combined transport of and its paired ion is dependent only on the transmembrane gradients in both ions; it is independent of the cell membrane potential. Consistent with recent experimental evidence42, transport proteins are assumed to be upregulated when cellular carbon utilization increases, which increases the density of transport proteins in the membrane, and effectively increases the bidirectional permeability of any given membrane to .

The model uses a flexible description of cross-membrane transport that can account for—but does not prescribe—facilitation by cotransporter proteins4142 (including counter-gradient transport) and upregulation of transport protein density in response to CO2 limitation. Based on the results of Bach et al.42, who showed that the transcript abundance of a putative transport protein increases at low DIC concentration, we here assume the effective membrane permeability to to be a linear function of utilization (Supplementary Fig. 2), which we define as the ratio of carbon use to passive diffusive supply of CO2. Effective membrane permeability to includes two constants; a background value, and the utilization coefficient. Each of these constants is explicitly allowed to be zero—that is, representing zero background contribution of to carbon supply, and no upregulation at high utilization. As transport can be coupled to Na+ symport or Cl− antiport via different cotransporter proteins in the SLC4 family, the transmembrane transport of is potentially dependent on the concentrations of these cotransported species. By describing the flux as a product of mass action (featuring the transport protein, , and any cotransported ion species), we can capture this dependency with an asymmetric permeability. For simplicity, in the case of cotransporter-facilitated transport, we assume that electroneutral Na+-coupled symport dominates transport, driven by the large extra-cellular to intra-cellular gradient of Na+, and that intracellular transmembrane gradients of Na+ are small. We note that the electroneutrality of this process, and thus its independence from the membrane potential, is an assumption, because some Na+-coupled SLC4 transporters are known to have Na+: stoichiometries other than 1:1, and these transporters have not yet been physiologically characterized in coccolithophores. We include a third constant representing the effective extra-cellular to intra-cellular concentration ratio of the cotransported ion species, which allows for, but does not prescribe, up-gradient bicarbonate transport across the cell membrane. Whilst being able to account for usage as described, the model can fall back completely to passive diffusion of CO2 if warranted by observations used to constrain the universal constants.

Through isotopic discrimination at the enzymatic level, carbon fixation in the chloroplast preferentially removes the light isotope of carbon from the CO2 pool, which itself is isotopically lighter than the pool. These effects combine to drive both chloroplast carbon pools towards isotopically heavy values. Calcification in the coccolith vesicle meanwhile preferentially removes the heavy isotope of carbon from the pool in the coccolith vesicle, which itself is isotopically heavier than the CO2 pool. Thus calcification drives the coccolith vesicle pool towards light values. In this carbon isotope flux model, steady state is assumed whereby fluxes in and out of any given pool sum to zero, but DIC species within all compartments are allowed to be in chemical and isotopic disequilibrium. The extent of disequilibrium is a function of the rates of calcification and carbon fixation relative to those of supply and interconversion of DIC species in a given compartment. Carbonate system reactions in any given system are quantifiable functions of DIC-species-specific concentrations, H+ concentration (pH), and the concentration of carbonic anhydrase (CA), the metalloenzyme responsible for catalysing the otherwise sluggish interconversion of CO2 and (refs 44, 45, 46) (see Supplementary Fig. 1).

Modelling and experimental results

The model is driven with empirical data for cell size (r), division rate (μ), the calcification to net carbon fixation ratio (reflected in the particulate inorganic carbon to particulate organic carbon ratio; PIC:POC), external CO2 concentration ([CO2]) and pH, and is calibrated against the carbon isotopic deviation from DIC of the resultant organic matter (O) and calcite (C). To this end we conducted in vivo carbonate chemistry manipulation experiments, growing two strains each of two species of coccolithophore across four DIC levels. Here we present the first experimental data set that includes all of the parameters required to constrain the model.

The data and model show that coccolith calcite carbon vital effects, are positive (heavier than abiogenic calcite) when the PIC:POC ratio is low, and are negative (lighter than abiogenic calcite) when the PIC:POC ratio is high (Fig. 2). The value of PIC:POC where this divergence occurs is ∼1 (Fig. 3). Vital effects are largest at high cellular carbon utilization—that is, when the cellular carbon requirement is an appreciable fraction of the carbon supply to the cell—and disappear when utilization is low. The carbon isotopic composition of organic matter meanwhile plateaus at a light value of ∼−23 to −24‰ (relative to DIC) at low utilizaton, and is driven towards heavier values at low utilization.

Figure 2

Model output.

Model output across the experimental range of CO2, with empirical data superimposed. The black horizontal dashed line shows the expected composition of abiogenic calcite. Universal constants used to generate model output are constrained by the points outlined in bold, here re-sampled 1,000 times from distributions defined by the associated analytical uncertainties. The data from this study also include additional G. oceanica data (calcite only; green filled circles) and organic carbon isotopic compositions for C. pelagicus from samples from a previous study9. See Supplementary Data 1 for values. Additional literature data are from dilute batch experiments manipulated by varying DIC9 (calcite carbon isotopic compositions only: C. leptoporus—dark blue upwards triangles; C. pelagicus—blue squares; E. huxleyi—yellow downwards triangles; and P. placolithoides—orange diamonds), and by varying pH (ref. 75; E. huxleyi—organic carbon isotopic compositions only: small yellow filled circles). Further data are superimposed in Supplementary Fig. 3. Data points are all the average of two replicates, with the ends of bars representing the value of each replicate. Solid lines and shaded envelopes represent respectively the mean and standard deviation of isotopic compositions predicted by the model for the 1,000 repeated calibrations. Values for compartmental pH are taken from published measurements33. Different species are represented by a representative set of parameters, which for ease of illustration are held constant with varying [CO2]. Representative values are taken from across a range of sources in the literature and from our own unpublished data (cell radius, division rate, PIC:POC are as follows: E. huxleyi – 2.3 μm, 1 day−1, 0.5; G. oceanica - 3.5 μm, 1 day−1, 1.1; C. pelagicus – 8 μm, 0.9 day−1, 1.2; C. leptoporus – 6 μm, 0.65 day−1, 2.2; and P. placolithoides – 7 μm, 0.8 day−1, 0.3). Strong co-variance between cell size or growth rate with CO2 would not be well represented in this projection. Note: for the pH manipulated model output uptake is enhanced at low CO2, manifest as a more rapid isotopic enrichment at low CO2 (steep yellow model curve). All species curves are a similar shape controlled by a number of trade-offs, and with the point of inflexion and maximum δ13Ccalcite determined by varying combinations of cell size, PIC:POC and growth rate. There are some regions of parameter space that are theoretically possible, but may not be observed in reality due to poor growing conditions. Regions of model space not populated by data should be treated as hypothetical.

Figure 3

Effect of cellular utilization and PIC:POC on isotopes.

Summary of model output—isotopic compositions as a function of PIC:POC and of the utilization parameter, τ. τ approximates utilization, and is a dimensionless linear function of growth rate, cell size [DIC], as defined in equation (1). (a) Calcite carbon isotopic compositions depend strongly on PIC:POC—if PIC:POC is high, calcite becomes isotopically light as utilization increases, if PIC:POC is low, calcite becomes isotopically heavy as utilization increases. (b) The carbon isotopic composition of organic matter is largely independent of PIC:POC, and is a strong function of the utilization parameter.

The observed phenomena are due to competing Rayleigh-type fractionation processes within the cell, which affect the isotopic composition of the bicarbonate pool in the coccolith vesicle, from which we assume calcite is precipitated. In preferentially incorporating the heavy isotope of carbon, the process of calcification drives its own substrate pool towards isotopically light values, with a greater isotopic drift manifest when the ratio of carbon usage to supply (utilization) in the coccolith vesicle is high. In the chloroplast meanwhile, the preferential removal of isotopically light carbon from the CO2 pool during carbon fixation catalysed by RuBisCO, causes the isotopic composition of CO2 to drift towards heavier values. The leakage of carbon from the chloroplast therefore becomes isotopically heavier, but also volumetrically decreased as the fraction of carbon supplied to the CO2 pool in the chloroplast is fixed and the concentration of CO2 decreases. The maximum influence of carbon fixation on the isotopic composition of the cytosolic pool, and thus on the coccolith vesicle and calcite, therefore occurs at intermediate values of CO2 utilization in the chloroplast (Fig. 4). Taken to either limit; zero; and complete utilization, carbon fixation in the chloroplast has no isotopic effects on the rest of the cell.

Figure 4

Summary of causes of vital effects.

Summary of causes of carbon isotope vital effects in coccolith calcite. (a) Coccolithophore cell with isotopically depleted POC produced in the chloroplast and isotopically enriched PIC produced in the coccolith vesicle. (b) Schematic representation of the effect of three dominant processes affecting the isotopic composition of calcite (δ13Ccalcite). (c) Net result of carbon isotopic vital effects in coccolith calcite. At low τ (that is, high CO2, low cell size or low growth rate), the flux of carbon into the cell far exceeds intracellular processes, and all species plateau at the same value (I). The value of this plateau is heavier than that of the carbon entering the cell due to intracellular interconversion of carbon species and preferential loss of light carbon dioxide. Low PIC:POC: for a cell of low PIC:POC, as τ increases (I–II), the flux of isotopically heavy carbon leaving the chloroplast (depleted in light carbon by carbon fixation catalysed by RuBisCO) influences the isotopic composition of the coccolith vesicle. As τ increases further (II–III), the flux of carbon from the chloroplast, although isotopically increasingly heavy, tends towards zero, and Rayleigh-type fractionation within the coccolith vesicle itself (which removes isotopically heavy carbon) takes over, driving the pool to light values. High PIC:POC: For a cell of high PIC:POC, as τ increases (I–IV) the Rayleigh-type fractionation due to calcite precipitation drives the isotopic composition of the coccolith vesicle to light values. The effect of heavy carbon leaking from the chloroplast is obscured. Across all cells of all PIC:POC values, the trends in vital effects are moderated by the biological up-regulation of HCO3− transport proteins, which increases overall permeability to carbon and increases the and isotopic composition of carbon entering the cell with increasing τ, and which dampens the Rayleigh fractionation-type effects when τ is high.

As carbon fixation and calcification have opposing effects on the isotopic composition of the bicarbonate pool in the coccolith vesicle, the net effect on the isotopic composition of calcite is a function of the ratio of these fluxes. When the calcification to net carbon fixation ratio, recorded physically in biogenic material as the molar PIC:POC ratio, is high, the effect of calcification dominates causing vital effects to be negative. When the PIC:POC ratio is low, the influence of isotopically heavy carbon leaking from the chloroplast is greatest, resulting in positive vital effects. Vital effects are greatest when carbon utilization is high (at low [CO2], or high μ or r), because the effective intracompartmental isotopic drift is greatest.

Cellular compartmentation and the non-infinite permeability of cellular membranes introduce the importance of intracellular heterogeneity, and the disproportionately large isotopic effect that any fractionating removal process has on its own substrate pool. Calcification, having a small fractionation factor, has very little effect on the carbon isotopic composition of the cell as a whole, and thus organic matter, but due to vicinity, does have a large effect on the pool in the coccolith vesicle. By contrast, the effect of carbon fixation on the isotopic composition of calcite is significant, because although the calcification substrate pool and carbon removal by fixation are separated by two membranes, the magnitude of the fractionation factor is comparatively very large.

At the limit of chemical and isotopic equilibrium, the model behaviour becomes a function of just two variables; the utilization index (τ), a compound parameter defined as the rate of cellular carbon usage to CO2 supply, and PIC:POC (see ‘Methods' section). At PIC:POC=0, as in non-calcifying eukaryotic phytoplankton such as diatoms, τ is the sole variable describing the isotopic composition of organic matter. τ, defined as:

where r=cell radius, ρ=cellular carbon density, μ=division rate, PC=membrane permeability to CO2, and Ce=extracellular CO2 concentration. τ is analogous to the utilization parameter described in the classic literature192232. PIC:POC and τ are parameters that can be mathematically constrained without assuming either chemical or isotopic equilibrium, and although they do not describe the system in its entirety when not at equilibrium, they are nevertheless useful for describing the system in a meaningful two-dimensional space (Fig. 3).

Discussion

Carbon isotopic compositions of coccolith calcite, and of organic matter, can be fully explained via competing Rayleigh-type fractionation processes within the compartmentalised cell (Fig. 4). The recently observed vanishing coccolith vital effect at high [CO2]79 naturally emerges as a result of the model, as Rayleigh-type fractionation is at a minimum when the rate of supply of carbon far exceeds usage.

There is an evolving consensus that forms a larger proportion of the carbon entering the cell at low [CO2]264748, a conclusion supported by a transcriptional upregulation of a putative cotransporter protein at low [DIC]42, and directly measured CO2 and uptake rates in E. huxleyi83738. Experimental and molecular evidence suggests that in coccolithophores, the transmembrane supply of is facilitated by transporters belonging to the solute carrier 4 (SLC4) family of cotransporter proteins40495051, which facilitate Cl−/ antiport, or Na+/ symport. Inside the cell, we assume that the transmembrane concentration gradients of Cl− and Na+ ions are negligible. In this scenario, transmembrane fluxes are proportional to the concentration on the proximal side of the membrane; the concentrations of cotransported ion species are incorporated in the effective permeability to . However, at the cell membrane, we allow the extra-cellular to intra-cellular concentration ratio of the cotransported ion species to vary, thus allowing up-gradient movement of . Due to the concentrations of both Cl− and Na+ likely being higher outside the cell than in the cytosol, we assume that symporting proteins dominate transport.

Other mechanisms for transport have been proposed which have implications for the isotopic composition of coccolith calcite. It has been suggested that negative carbon isotopic vital effects in coccolith calcite are due to an increased pumping of into the chloroplast at low [CO2], at the expense of moving to the coccolith vesicle with no diffusive efflux of from any compartment, leaving the cytosolic pool, and thus the coccolith vesicle, isotopically light in carbon2. We are aware of no molecular evidence to date of putative primary active transport proteins associated with bicarbonate transport in this group of organisms however. Furthermore, prohibiting efflux from any cellular compartment has significant implications for the isotopic composition of intracellular carbon pools, and imposes an artificially high importance on the ratio of inward fluxes of to CO2. Although effective membrane permeabilities to CO2 or have not been measured directly in coccolithophores, the values constrained by our model (permeability to CO2=∼9.3 × 10−4 ms−1, and to being 3–4 orders of magnitude lower) are consistent with those calculated for the closely related group, the diatoms46 (see ‘Methods' section).

In some unicellular algae, CA is emitted from the cell to accelerate sluggish kinetics of dehydration in the boundary layer. Despite reports that E. huxleyi lacks external CA (CAe)37, more recent work argues for the presence52 and possibly upregulation of CAe at low [CO2]4142. CAe was not however detected in either Gephyrocapsa oceanica nor C. pelagicus, even when cells were carbon limited49, though more research on species other than E. huxleyi is required. Here we assume that CAe is present in high enough concentrations that the external surface of the cell is in chemical and isotopic equilibrium with the bulk seawater. If future work concludes that CAe is absent in one or more species, it may be necessary to allow chemical and isotopic disequilibrium in the cellular boundary layer (for implications of this scenario see Supplementary Fig. 4 and Supplementary Note 1).

The model described here does not separately represent the pyrenoid: a compartment of low pH nested within the relatively alkaline chloroplast where CO2 and RuBisCO are thought to be highly concentrated. The value of the effective modelled isotopic discrimination factor for RuBisCO, describing the difference between the isotopic compositon of CO2 in the chloroplast and that of organic matter (fe), is therefore effectively a pyrenoid/RuBisCO black-box fractionation for the given set-up. The extent to which fe is representative of the enzymatic RuBisCO fractionation, , is a critical question regarding the validity of using in vitro measured values of in models that do not have a detailed representation of the pyrenoid. The introduction of an additional CO2—permeable membrane between the ambient medium and the site of carbon fixation, may allow a further isotopic discrimination to take place by the preferential uptake of CO2 over , which would allow the value of to be even smaller than the value of fe predicted by our model—thus providing an explanation when there is a discrepancy between low values of and high values of . However, conversely, if compartmental utilization of the pyrenoid were very high, may be far higher than the value of fe predicted by the model. Although isotopic fractionation by RuBisCO has been extensively studied in plants, only two studies to date have estimated the in vitro isotopic fractionation of carbon associated with carbon fixation catalysed by the type of RuBisCO (form ID) found in the ‘red-lineage algae' such as the coccolithophores and the diatoms3953, and these values are surprisingly low (11 and 19‰ respectively). Despite this, recent modelling work has continued to use the markedly larger fractionation factors associated with plants and cyanobacteria (Table 1). The fitted value for fe output by our model is ∼14–15‰.

The calcite saturation state in the coccolith vesicle (Ωcalcite) is controlled by [] and [Ca2+], and by the solubility product of calcite in seawater (Ksp), which is a function of salinity, temperature and pressure4454. At most CO2 concentrations, the model suggests that the coccolith vesicle is oversaturated with respect to calcite (Supplementary Fig. 5). However, at CO2 concentrations below around 30 μM, the model predicts that Ωcalcite falls below unity, which would usually render calcite thermodynamically unstable. A number of possibilities may individually, or in combination, lead to a precipitation-promoting Ωcalcite inside the CV at this lower CO2 range. First, measuring pH in an intracellular compartment is difficult, and is fraught with uncertainties, so the pH in the coccolith vesicle may be significantly higher than the previously measured value of 7.1 (ref. 33), increasing the fraction of DIC present in the form . Second, transport to the coccolith vesicle may be enhanced at low CO2− this could be an upregulation of transport proteins42 specific to the coccolith vesicle membrane, or, as suggested elsewhere, an energy-driven net Ca2+/ co-transport mechanism48. Third, the solubility product (Ksp) of calcite is two orders of magnitude higher (that is, Ωcalcite is lower and calcite precipitation is less likely to occur) in water with a salinity typical of seawater compared with distilled water54. Although the salinity of the CV is likely to be close to that of seawater due to the osmotic pressure that is associated with large salinity differences55, maintaining a lower salinity in the coccolith vesicle is a possible mechanism of elevating Ωcalcite. Fourth, through forming temporary bonds to Ca2+ ions in solution, acidic polysaccharides, which are known to be present in the CV, and throughout the coccolith calcite lattice, strongly influence calcification, and may cause the localized Ωcalcite to be significantly higher than that of the bulk CV. Values of [Ca2+] in the CV have not been measured directly with much success3356, but hypothetical estimates vary across an order of magnitude (0.5–4 mM)475758. The localized calcite precipitation-inducing effect on carbonate chemistry may even allow coccoliths to be built in a mostly calcite under-saturated vesicle. Little work has been undertaken to understand the localized influences of polysaccharides on carbonate chemistry however.

Under the modelled range of [CO2], CO2 of the chloroplast is at a level comparable with the Michaelis–Menten half-saturation constant of RuBisCO typical for a coccolithophore (∼30–40 μM CO2; Supplementary Fig. 5)3439, as might be expected for maximum efficiency. At low [CO2], CO2 in the chloroplast does drop significantly below this value, however, the pyrenoid/thylakoid complex is a small region of very low pH where RuBisCO and CO2 are in far higher concentration than in the rest of the chloroplast3459. Locally therefore, the low [CO2] inferred here for the chloroplast is unlikely to be reflective of the environment immediately in the vicinity of the RuBisCO enzyme.

Oxygen kinetic isotopic fractionation factors (KIFs) for the hydration and dehydration of CO2 and respectively are not known, and theoretical predictions vary widely6061. As a result, this model cannot currently be fully resolved with respect to oxygen isotopes. An additional complexity over carbon isotopes arises with oxygen because oxygen in DIC exchanges with oxygen in water, which even at high DIC concentrations is many orders of magnitude more abundant446263. The rate of oxygen exchange with water and the residence time of DIC species in a compartment is therefore critical to the preservation of vital effects in calcite. As coccolith calcite exhibits vital effects in oxygen, the cell cannot be in complete isotopic equilibrium with respect to oxygen. Across coccolithophores, a positive correlation is seen between the carbon and oxygen isotopic compositions of calcite (Supplementary Fig. 6a)9. However, although is isotopically heavier than CO2 in carbon, it is lighter than CO2 in oxygen, so the observed trends in carbon and oxygen isotopes cannot be simultaneously explained solely by a shift in cellular carbon source from dominantly to CO2 or vice versa. We hypothesize that the increased difference in oxygen isotopic vital effects between large and small cells observed at at low [CO2]23459 is due to a difference in the oxygen isotopic composition of DIC entering the cell following the increased influx of bicarbonate at high utilization (Supplementary Fig. 6b). We further hypothesize that this signal may be superimposed on a Rayleigh-type fractionation, which drives the combined / pool in the coccolith vesicle to heavy values at high utilization, as the isotopically lighter (in oxygen) is taken up into calcite, and by CA-catalysed hydration and dehydration which pushes the system further towards equilibrium at high DIC (Supplementary Fig. 6c).

The model of vital effects in coccolith calcite presented here not only accounts for previously unexplained variations in isotopic compositions, but also provides the theoretical framework for a new suite of proxies. The ultimate ambitious aim is to accurately quantify pCO2atm in the geological past, but if the whole system can be constrained, this would additionally yield information aspects of contemporaneous coccolithophore physiology. Coccolithophores include the only family of organisms known to be responsible for producing alkenones—the most widely used organic molecule extracted from ancient sediments for estimating past pCO2atm via isotopic analysis. In these analyses, traditionally isotopic compositions have been related to [CO2], and thus pCO2atm, through univariate and empirical relationships2728293031. Here we have shown that the system is more complicated than this, and the isotopic fractionation of carbon into calcite (and to a lesser degree organic material) is strongly bivariate in PIC:POC and τ space even when chemical and isotopic equilibrium can be assumed (Fig. 3). With paired measurements of C and O from analysis of ancient sediments, in theory our model can be used to iteratively search parameter space to minimize the misfit between observed and predicted isotopic compositions, and thus simultaneously predict the most likely values of PIC:POC and τ of these ancient organisms. We recently extracted acidic polysaccharides from within the calcite lattice of large ancient coccoliths64, which opens up the possibility of isotopically characterizing non-alkenone producing species, and may, in the future, supersede alkenones as the target molecule for organic carbon isotopic analyses. If, in addition, an independent estimate of PIC:POC were obtained65, the difference between the carbon isotopic compositions of organic matter and calcite alone could be used to constrain the system, circumventing the need for an independent estimate of the isotopic composition of DIC (Fig. 5). To extract an estimate of [CO2] from the constrained value of τ, growth rate and cell size must be known. The difference in isotopic composition of calcite produced by different sizes of coccolithophore cells growing in the same seawater may also be used to constrain [CO2], but independent estimates of PIC:POC, growth rate and cell size for both coccolithophore cell size fractions would be necessary (Fig. 5). For some species of coccolithophore, cell size can be estimated from coccolith size1165, and PIC:POC can be estimated from the cross-sectional aspect ratio of a coccolith65. Strontium/calcium ratios in coccolith calcite have also been proposed to reflect growth rate6667. At present, however, there are large statistical uncertainties associated with all of these estimates, which are multiplicative when propagated. Our novel mechanistic insight provides a biologically grounded theoretical basis, but improvements in estimates of these other integral parameters are necessary before accurate quantitative coccolith-based paleo-CO2 barometry will be possible.

Figure 5

Concept for interpretation of sedimentary data.

(a) Once calibrated, the model provides values of isotopic compositions of organic matter and of calcite (δ13Corg and δ13Ccalcite respectively) as outputs, given a number of input parameters (δ13CDIC, r, μ, [CO2] and PIC:POC). δ13CDIC can be estimated from foraminifera. r, μ and [CO2] form a single compound variable for each size fraction (semi-transparent ovals), τ (equation (1)), which cannot be deconvolved without additional evidence. The model can be inverted and used to iteratively search parameter space for the values of τ and PIC:POC that minimize the misfit between model-predicted and measured isotopic compositions. r can be estimated from coccolith size. Boxes without a division denote a shared variable between measured isotopic compositions. O is only weakly dependent on PIC:POC so to a first approximation, τ may be estimated from O alone (as shown by shaded box). The assumption included in [CO2] concerns the ambient concentration of , but the model output is highly insensitive to this assumption. (b–d) Model output describing theoretical carbon isotopic compositions of inorganic (b) and organic (c) material relative to δ13CDIC, and δ13Ccalcite—δ13Corg (d) over a range of PIC:POC and τ. The double headed arrows point to the difference between a cell with a PIC:POC typical of E. huxleyi (left) and a typical C. pelagicus (right); the number of δ13Ccalcite contours crossed decreases with decreasing τ.

Methods

Culture experiments

Culture experiments were undertaken in the Department of Earth Sciences at the University of Oxford. Duplicate monoclonal batch cultures of four strains of coccolithophore belonging to the family Noëlaerhabdaceae were grown in sterile filtered (0.2 μm) artificial seawater prepared according to ESAW68 adapted for a range of DIC concentrations ((DIC)=1.380 mM, 2.147 mM, 3.067 mM and 6.135 mM) at constant pH total scale (8.2) by varying sodium bicarbonate addition and titration with HCl and with nitrate (442 μM), phosphate (5.00 μM), vitamins, trace metals and ethylenediaminetetraacetic acid according to K/2 (ref. 69). Cultures were maintained at 15 °C with an incident photon flux of 250 μE and a 12/12 light/dark cycle. Cells were acclimated for >20 generations in dilute batch culture for each experimental condition before inoculation. Cells were inoculated in 2.4l polycarbonate flasks, with no headspace and sealed off to the air with teflon lined caps. Removal of medium during the experiment was unavoidable due to the need to count and measure cells, and resulted in a maximum headspace of 20 cm3 at harvest. To minimize the drift in culture conditions throughout the course of the experiment, cells were harvested at ∼1–2% (and never >4%) of maximum cell density, which was determined for each experimental condition and strain combination via prelimenary experimentation. Strains were AC478 (RCC1211 G. oceanica from Portuguese coast in Atlantic Ocean), AC472 (RCC1216 E. huxleyi, from Tasman Sea in Pacific Ocean), AC448 (RCC1256 E. huxleyi, Icelandic coast in Atlantic Ocean) and AC279 (RCC1314 G. oceanica, French coast in Atlantic Ocean) from the Roscoff culture collection (RCC). Particulate material was harvested by dry filtration onto pre-weighed membranes with 0.2 μm pore-size, and rinsed of salt with a minimal amount of neutralized deionised water. Coccolithophore size and concentration were obtained using a Beckman Z2 Coulter Counter. Coccosphere and cell size were measured three times each respectively pre- and post-decalcification both morning and evening on the harvest day and the preceding day. Cells were decalcified by reducing the pH of the suspension with HCl addition to 5.0 with for around 20 min. The Coulter counter was calibrated to use ESAW+K/2 medium as an electrolyte, and for use with the acidified electrolyte, to accommodate for the difference in ionic strength. Cell division is synchronized under the light/dark cycle and cell size was assumed to increase linearly throughout the day70. By measuring cell and coccosphere size morning and evening, the bias introduced due to the time of day of measurement can be removed by interpolation to the same time of day. This also removes the daily variation of cell size70. Culture health was monitored by cell counts and microscope inspection on alternate days. Molar PIC and POC were measured with a Rock Eval analyser in the Earth Sciences department at Oxford University.

Isotope measurements

Carbon and oxygen isotopic compositions of calcite were measured at the University of Oxford using a VG Isogas Prism II mass spectrometer with an on-line VG Isocarb common acid bath preparation system. Samples were first rinsed with neutralized deionised water to remove any salt. Samples were then dosed with acetone and dried at 60 °C for at least 30 min. In the instrument they were reacted with purified phosphoric acid at 90 °C. Calibration to PDB standard was via the international standard NBS-19 using the Oxford in-house (NOCZ) Carrara marble standard. Reproducibility of replicated standards was better than 0.1‰ (1σ) for δ13C and δ18O expressed relative to the V-PDB standard. Carbon isotopic composition of organic material was measured on an automated carbon and nitrogen elemental analyzer (Carlo Erba EA1108) at the Research Laboratory for Archaeology and the History of Art at the University of Oxford. Samples were decalcified with HCl, and rinsed with MilliQ water at least three times before being weighed into tin capsules. The internal alanine standard reproducibility is ∼0.16‰ expressed relative to the V-PDB standard.

Assumptions of the model

The assumptions on which this model is based, are:

Cells are in steady-state growth. The concentration of every carbon pool (that is, each carbon species in each intracellular compartment) is constant, which implies that all fluxes associated with that pool sum to zero.

The cell and compartments are assumed to be isomorphic across all strains and conditions.

CO2 and ions move across membranes via passive diffusion and facilitated diffusion respectively. Both chemical species are assumed to move across membranes with a flux proportional to surface area of the membrane and the concentration of the carbon species on the proximal side of the membrane. There is no concentration-independent active uptake of bicarbonate or carbon dioxide. Membranes are impermeable to carbonate ions. There is no other source of carbon to the cell.

Loss of carbon from the cell is via carbon fixation, calcification, passive diffusion of CO2, and co-transport of .

The membrane permeability to CO2 is constant.

The effective membrane permeability to , which incorporates the concentrations of coported ions, is a linear function of background carbon utilization (U 0; defined in equation (2)), representing the ability of cells to upregulate transport proteins. The special cases whereby permeability is non-adaptive, or is zero are explicitly allowed.

The constants describing the permeability of membranes to CO2 and to are identical across all species, and for all membranes except the cell membrane where the strong gradient of Na+ from the extra-cellular to intra-cellular environment is assumed to dominate transport via Na+ coport, thus introducing an asymmetric permeability. As this is an electroneutral process, the membrane potential does does affect this flux, which is assumed to be a product of mass-action.

pH in each compartment is constant, and is prescribed a priori, according to the values measured by Anning et al. 33.

Organic matter is not remobilised from the organic pool to the cytosol via mitochondrial respiration.

Movement of carbon across membranes does not impart a kinetic isotopic fractionation.

All cells are assumed to exhibit active external CA (see Supplementary Note 1).

The chloroplast is assumed to consist of a single compartment with no pyrenoid. The large isotopic fractionation of carbon that occurs during fixation catalysed by the enzyme RuBisCO , is here factored into the model as an effective pyrenoid/RuBisCO black box fractionation ( fe).

Membrane permeabilities

A putative bicarbonate transport protein has been observed to be upregulated of at low [CO2]42. Here these results are interpreted as consistent with an approximately linear increase in transcript abundance with the degree of carbon utilization (Supplementary Fig. 2). This physiological result is assumed to reflect an increase in synthesis of transport proteins, and thus the density of transport proteins in the membrane, manifest as an increase in the effective permeability of membranes to . The flux of through membranes facilitated by SCL4-like exchange proteins is proportional the the product of on the proximal side of the membrane and Na+ on the proximal side of the membrane for coport, and that of Cl− on the distal side of the membrane for antiport. For simplicity, we assume transport to be driven by Na+ coport, and the gradients of the coported ion to be negligible across all membranes, except the cellular membrane where it is substantial. This assumption is factored into the model via an additional universal constant, which describes the extra-cellular to intra-cellular concentration gradient of the coported ion. This assumption is more coherent with the current biological literature than equivalent assumptions of previous models, including concentration-independent implicitly ATP-driven active uptake of , assumed by Keller et al.24 to scale with growth rate, or by Bolton and Stoll2 to be independent of all other parameters.

In this model, the permeability of membranes to CO2 (PCcell) is assumed to be constant. The permeability of membranes to (PBcell) is assumed to increase linearly with utilization of background (that is, before upregulation of anion exchange proteins) carbon supply (U0), which is defined as:

where Cin0 is the net carbon supply when membrane permeabilities to CO2 and to are at their background values, and FFIX+FCAL is the rate of fixation of carbon into organic and inorganic matter. The membrane permeability to bicarbonate is therefore:

where T0 is the base-level membrane permeability to and TU is the gradient of the increase in membrane permeability to with increasing utilization; FFIX and FCAL are the rates of carbon fixation and calcification respectively, Ce and e are the concentrations of CO2 and in the ambient medium respectively, and SAcell is the surface area of the cell. The large Na+ gradient at the cell membrane is included as an additional constant, G. This line of reasoning is a set of assumptions factored into the model via four universal constants describing the membrane permeability to CO2 and to (PCcell, T0, TU and G). These four constants will be constrained by the data, explicitly leaving open the option of non-upregulated transport (through TU=0), complete impermeability of membranes to (through T0=TU=0), and through up- and down-gradient movement of dependent on the inferred Na+ ion gradient.

Compartment shapes and sizes

Each intracellular compartment is assumed to be an oblate spheroid. The equatorial axis of the spheroid, a, is assumed to have a constant ratio with the cell radius. This is referred to as the scaling factor (s fc × r=ac). The polar axis, c, is assumed to have a constant ratio with the equatorial axis. This ratio is referred to as the aspect ratio factor (a fc × ac=cc). Isometry is assumed across species, and in this way, two constants are used to describe the relative size and shape of each compartment, and therefore their volumes and surface areas. The surface area of an oblate spheroid is given by:

The volume of a spheroid is given by:

Intracellular carbonic anhydrase

The CAs are a family of zinc-containing metalloenzymes responsible for catalysing the hydration and dehydration of CO2 and respectively. Their behaviour is well described by Michaelis–Menten kinetics45. Given the general form of an enzymatically catalysed reaction:

where E, S, ES and P denote enzyme, substrate, enzyme–substrate complex and product concentrations, respectively. kf, kr and kcat are rate constants respectively for the binding and unbinding of substrate to the enzyme, and for the maximum catalytic throughput of the reaction.

According to the Michaelis–Menten equations, the velocity of a reaction such as equation (7) is given by:

where Km is the Michaelis constant, or half-saturation constant , and describes the substrate concentration when v=. When [S]<

and is the specific activity of the enzyme in units of M−1 s−1. This is analogous to a rate constant (kp), or activity, when multiplied by the concentration of the enzyme. The rate constant for CA-catalysed hydration of CO2 is then:

As enzymes catalyse reactions but do not alter the position of equilibrium, the rate constant for the reverse reaction, CA-catalysed dehydration of , is given by the hydration rate constant divided by the equilibrium constant:

A value of 2.7 × 107 M−1 s−1 was determined by Uchikawa et al.45 for the of the hydration reaction catalysed by bovine erythrocyte CA.

Model derivation

Fluxes of CO2 through membranes are the product of membrane surface area (SA), CO2 concentration on the source side of the membrane, and membrane permeability to the specific carbon species. From Fig. 1 expressions for transmembrane CO2 fluxes are:

and likewise for fluxes:

where all fluxes are defined in Fig. 1. G is the constant that describes the effective asymmetry of the cell membrane to due to the transmembrane Na+ gradient.

In each compartment, carbon dioxide and bicarbonate ions are interconverted by the reversible hydration and hydroxylation reactions, and by hydration in the presence of CA:

where arrow annotations denote rate constants (see section for CA rate constants). Combining the rate constants of equation (14) at constant pH and [CA] in a generic compartment, a, gives:

where kCBa denotes the reaction rate constant for the conversion of CO2 to , and kBCa denotes the rate constant for the reverse reaction. Ha, OHa and CAa refer to the concentrations of H+ ions, OH− ions and CA in compartment, a. As compartment-specific pH and [CA] are defined apriori, kCBa and kBCa can be treated as pH and [CA] dependent compound rate constants. The rates of the reactions described by equation (14) in each compartment are therefore:

Carbon fluxes throughout the cell can be fully described with equations (12), (13) and (17), and with two additional output fluxes; the rate of calcification (FCAL) and the rate of photosynthetic carbon fixation (FFIX; Fig. 1). Assuming steady state, the rate of change of the amount of each carbon species in each compartment is zero, giving:

At steady state these differential equations become linear functions of the concentration of the different inorganic carbon species. Substituting in equations (12), (13) and (17):

This set of interdependent equations can be written as a linear system of the form:

where A is the coefficient matrix of the linear system, N is the nonhomogeneous term and Φ is the unknown vector containing the carbon species concentrations in each compartment. Equation (20) can be solved for Φ, where A and N, defined by the dynamic carbon species concentration equations, are as follows:

[13C] is very low compared to [12C], carbon fluxes are assumed to represent 12C. 13C dynamics are therefore determined by prefixing each C flux with corresponding R, where:

By definition,

so the assumption of balanced growth allows both offset Rstandard and scale factor Rstandard to be eliminated while preserving the same system of equations. Effectively, we can thus prefix each flux directly with the associated δ13C value. For fractionation fluxes, this leads to a prefix with the sum of the δ13C of the source pool and the process-specific fractionation factor (). In the following, δΘa refers to the δ13C of carbon species, Θ, in compartment, a.

The following fractionation factors (in ‰)44, are assumed constant:

where the subscripts are consistent with those from equation (14). Analogous to the compound rate constants of equations (15) and (16), the interconversion of carbon dioxide and bicarbonate has a pH and [CA] dependent, compartment-specific, compound isotopic fractionation factor.

When the linear system described for carbon fluxes is solved, all carbon species concentrations in all compartments are known. Dynamic equations for the isotopic composition of each compartment, assuming balanced growth, and expressed in terms of carbon dioxide, and bicarbonate fluxes and fractionation factors are thus:

At steady state the differential equations become linear functions of δ13C of the different carbon species in the different compartments. Substituting in equations (12), (13) and (17):

This set of interdependent equations can also be written as a linear system of the form:

where A is the coefficient matrix of the linear system, N is the nonhomogeneous term and Φ is the unknown vector containing the carbon isotopic compositions of each carbon species in each compartment. Equation (36) can be solved for Φ, where A and N, defined by the dynamic carbon species concentration equations, are as follows:

Equilibrium limit special case

At the limit where CA activity is infinitely high in all compartments, complete chemical and isotopic equilibrium with respect to carbon can be assumed. In this scenario, the total [DIC] (D) in each compartment and its isotopic composition, along with (prescribed) pH, completely specifies the concentrations and isotopic compositions of the different inorganic carbon species in each compartment.

The dynamic equations for total DIC (D) in the three different compartments are:

In steady state all three equations must equal 0, which implies:

Inserting these into equation (37), the expression for cytoplasm total DIC, gives:

The four fluxes across the cell membrane equal:

where the permeability for relates to the permeability for CO2 through background utilization:

As a result, all cross-membrane fluxes are proportional to PCSA and equation (42) can be divided through by this factor. With chemical equilibrium and prescribed pHc in the cytoplasm, Bc is proportional to Cc for any compartment c∈{e, i, x, v}:

where is the temperature and salinity-dependent chemical equilibrium constant, and Hc the H+ concentration (Hc=). Thus equation (42) can be rewritten as:

Isolating Ci and divide through by Ce gives:

Thus, concentrations and fluxes of carbon species in all compartments are known.

Following the same process, two additional fluxes across the chloroplast membrane equal:

where fx is the ratio of the product of cellular membrane permeability and surface area to that of the chloroplast. Equation (40) can therefore be rewritten:

and inserting equation (50) becomes:

Similarly, two additional fluxes across the coccolith vesicle membrane equal:

where fv is the ratio of the product of cellular membrane permeability and surface area to that of the coccolith vesicle. Equation (41) can therefore be rewritten:

and inserting equation (50) likewise becomes:

The values fx and fv are constants when all membranes are assumed to have the same permeability, and isometry is assumed.

Membrane permeability to bicarbonate is reliant on the background utilization of carbon, U0. The value of U0 is given by:

At equilibrium, has a constant isotopic offset from CO2, denoted here as E. As both and CO2 in each compartment are known if one is known, the whole system can be described by CO2 alone. The dynamic equations for isotopic fluxes are:

As with chemical fluxes, substituting in expressions for concentration driven diffusion, membrane permeability to in terms of background utilization, concentration in terms of CO2 and H+, and rearranging gives:

The chemical carbon flux equations (50), (54) and (58) can be solved simultaneously. Further, the carbon isotope flux equations (63)–(65), , can also be solved simultaneously, using the chemical carbon flux equations. In these expressions, including U0, the rates of fixation and calcification always occur divided by PCSACe, and this is also the only appearance of these parameters. Thus,

are compound varibles that describe the entire system when constant and intracellular pH is assumed and known a priori. In this model, the fixation rate is related to cell size via the equation,

where ρ is the cellular carbon density and V and μ are the volume and division rate of a cell respectively. As PIC:POC is equal to , the compound parameters become:

and:

Therefore, at the limit, where intracellular CA is infinitely high such that all compartments are in chemical and isotopic equilibrium, the orthogonal variables that describe the entire system are:

τ is a parameter that is similar to the ratios that evolved through early work describing carbon fluxes in phytoplankton1617182123242532. The more complicated disequilibrium model collapses to a state completely consistent with this earlier work when equilibrium is assumed and PIC:POC=0. When the cell is close to chemical and isotopic equilibrium with respect to carbon therefore, changes in cell size, growth rate and extracellular [CO2] cannot be decoupled.

Fitted parameters

The results from culture were used to constrain PCcell (the permeability of membranes to CO2), T0, TU (the constants required to define the linear dependence of membrane permeability to on utilization (equation (4)), G (the asymmetry of the effective cell membrane permeability to , which may reflect the extra-cellular to intra-cellular concentration ratio of the coported ion), and fe (the effective enzymatic isotopic fractionation associated with carbon fixation by RuBisCO). Universal constants were fitted to the data by minimization of a misfit function, which quantified the sum of squared deviations of the modelled results from the measured data. Minimization of the misfit function was achieved using an R implementation71 of the constrained DIRECT_L algorithm72 within a defined parameter space. The model was prescribed to be undefined when concentrations are negative, or when the saturation state of calcite in the coccolith vesicle is below 0.1. Final parameter values were checked with an alternative implementation of the model in Python, using a combination of a broad search with the Differential Evolution algorithm73, followed by refinement of its result with the Nelder-Mead Simplex algorithm74. Only complete data sets were used to computationally constrain the universal parameters of the model. The required parameters are: cell size, PIC:POC, growth rate, medium pH and [DIC], δ13CDIC, δ13Ccalcite and δ13Corg. The model was fitted to these data multiple times with pseudo-replicate datasets generated by resampling each value from an assumed gaussian distribution defined by the expected value and uncertainty associated with each measurement. All input parameters were resampled in this way. The fitted value of PCcell is ∼9.3 × 10−4 ms−1, which is at the high end of the range suggested by Hopkinson et al.46 for diatoms. T0 is ∼7.8 × 10−8 ms−1; also consistent with values given by Hopkinson et al.46 Based on these results, at low utilization, the membrane is 3–4 orders of magnitude less permeable to as to CO2. The constant TU describing the increase in membrane permeability to with increased utilization is ∼5.1 × 10−6 ms−1 per unit utilization. The parameter G, which describes the asymmetry of the cell membrane permeability to , is constrained to be 2.7. This value partially reflects the Na+ gradient, but is somewhat too low, so probably also represents Cl−-coupled anion exchange, which would drive bicarbonate out of the cell. When utilization is low, either at high DIC or in small cells, ∼10% of carbon enters the cell in the form of . When cells experience high utilization, at low DIC or in larger cells, the contribution of to total DIC supply can be as high as 50%. CA in all compartments is between 0.1 and 1 mM in all intracellular compartments, based on a specific activity of 2.7 × 107 M−1 s−1 (ref. 45). Values below 0.1 mM lead to extremely unstable model output and increasing this value above 0.1 has little effect on model output suggesting that compartments are approximately in equilibrium with respect to carbon. As it takes an order of magnitude longer to reach equilibrium with respect to oxygen, and we know that the cell cannot be in equilibrium with respect to oxygen due to the existence of oxygen isotopic vital effects in coccolith calcite, the CA concentration must exist within this small range. fe is ∼−14.3‰, which is far smaller than that used in the literature to date221242532 but is closer to the in vitro measured value of given for E. huxleyi (−11‰ (ref. 39)) and for the diatom Skeletonema costatum (−19‰ (ref. 53)). Parameters associated with the model input and output are summarized in Table 2.

Table 2

Parameters from the model, their definitions, derivation and units.

Symbol	Definition	Derivation	Quantity	Units
Measured parameters
Ce, Be	External CO2 and concentration	Measured [DIC] and pH	Independent variable	mol m−3
δCe, δBe	δ13C of external CO2 &	Calculated from δ13C of DIC and pH	Independent variable	‰PDB
FFIX	Average carbon fixation rate	Division rate, mol C org per cell	Measured variable	mol s−1
FCAL	Average calcification rate	Division rate, mol calcite per cell	Measured variable	mol s−1
PIC:POC	Particulate inorganic to particulate organic carbon ratio of biomass	FCAL/FFIX	Measured variable	Molar ratio
TCR	Total carbon assimilation rate	FCAL+FFIX	Measured variable	mol s−1
δCorg	δ13C of organic material	Measured directly	Measured variable	‰PDB
δCcal	δ13C of calcite	Measured directly	Measured variable	‰PDB
rcell	Cell radius	Measured directly	Measured variable	m
SAcell, Vcell	Cell surface area and volume	Functions of rcell	Measured variable	m2&m3
Rho	Cellular organic carbon density	This study, consistent with Reinfelder77	20 × 10−15	mol m−3
Parameters from literature
	Carbon isotopic fractionation during calcification	Zeebe and Wolf-Gladrow44	+1 (from )	‰PDB
, , , , ,	Kinetic carbon isotopic fractionation factors	Zeebe and Wolf-Gladrow44	−13, −22, −11, −20, −1, −10	‰PDB
kp1, km1, kp4, km4	Kinetic rate constants	Calculated from Zeebe and Wolf-Gladrow44	T and S dependent.	various—dependent on reaction
,	Ion product of water	Calculated from (ref. 78)	T and S dependent.	mol2 m−6
,	1st acidity constant of carbonic acid	Calculated from Lueker et al.79	T and S dependent.	mol m−3
	Specific activity of CA for hydration reaction	Inferred from bovine erythrocyte CA45	2.7 × 107	M−1 s−1
pHi, pHx, pHv	pH in compartments i, x & v.	ref. 33	7.0, 7.9 & 7.1.	free scale
fx, fv	Scale factor for compartments x & v	TEM images8081	1.1 & 0.6	dimensionless
Sx, Sv	Shape factor for compartments x & y	TEM images8081	6 & 4	dimensionless
Fitted parameters
fe	Effective carbon isotopic fractionation during C-fixation by Rubisco	Fitted parameter	−14.3	‰PDB
PCei	Cell membrane permeability to CO2	Fitted parameter	9.3 × 10−4	ms−1
PBei (T0, TU)	Cell membrane permeability to HCO3−	Fitted parameter	7.8 × 10−8, 5.1 × 10−6	ms−1
G	Effective concentration ratio of the coported ion at the cell membrane	Fitted parameter	2.7	ratio
[CA]i, [CA]x, [CA]v	[CA] in compartments i, x & v.	Inferred	0.1–1	mol m−3

See supplementary information for values for measured parameters generated in this study.

Data availability

The culture data that support the findings of this study are contained in the Supplementary Material, and in the literature cited throughout the manuscript.

Additional information

How to cite this article: McClelland, H. L. O. et al. The origin of carbon isotope vital effects in coccolith calcite. Nat. Commun. 8, 14511 doi: 10.1038/ncomms14511 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.