Carbonate Disequilibrium in the External ... - ACS Publications

Michel Lavoie*† , Jérôme F. L. Duval‡§ , John A. Raven∥⊥ , Frédéric Maps† , Béchir Béjaoui# , David J. Kieber∇ , and Warwick F. Vin...
1 downloads 0 Views 2MB Size
Subscriber access provided by - UC Berkeley Library

Ecotoxicology and Human Environmental Health

Carbonate disequilibrium in the boundary layer of freshwater chrysophytes: Implications for contaminant uptake Michel Lavoie, Jerome F.L. Duval, J. A. Raven, Frédéric Maps, Béchir Béjaoui, David John Kieber, and Warwick Vincent Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b00843 • Publication Date (Web): 17 Jul 2018 Downloaded from http://pubs.acs.org on July 26, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 37

Environmental Science & Technology

1

Carbonate disequilibrium in the boundary layer of freshwater chrysophytes:

2

Implications for contaminant uptake

3

Michel Lavoie1*, Jérôme F. L. Duval2,3, John A. Raven4,5, Frédéric Maps1, Béchir Béjaoui6,

4 5

David J. Kieber7, Warwick F. Vincent1

6

1

7

Département de biologie, Université Laval, Québec, Québec, Canada G1V 0A6

8

2

9

UMR7360, Vandoeuvre-lès-Nancy F-54501, France.

Québec-Océan, Takuvik Unité Mixte Internationale (Université Laval-CNRS) and

CNRS, LIEC (Laboratoire Interdisciplinaire des Environnements Continentaux),

10

3

Université de Lorraine, LIEC, UMR7360, Vandoeuvre-lès-Nancy, F-54501, France

11

4

Division of Plant Science, University of Dundee at the James Hutton Institute, Invergowrie,

12

Dundee DD2 5DA, United Kingdom (permanent address of JAR)

13

5

14

Sydney, Ultimo, NSW 2007, Australia

15

6

16

rue du 2 mars 1934, 2025, Salammbô, Tunisie.

17

7

18

and Forestry, Syracuse, New York, USA 13210

19 20

*

21

Word count: 5633 words (from the Abstract to the end of the conclusion including

22

acknowledgments and supporting information abstract) + 2 figures x 600 words + 1 figure x

23

300 words = 7133 words-equivalent

24

Running head: Phytoplankton boundary layer and metal speciation

25

Keywords: pH, carbonate, boundary layer, metals, bioavailability, phytoplankton

Functional Plant Biology and Climate Change Cluster (C3), University of Technology

Laboratoire du Milieu Marin, Institut National des Sciences et Technologies de la Mer, 28,

Department of Chemistry, State University of New York, College of Environmental Science

corresponding author: [email protected]

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 37

26

ACS Paragon Plus Environment

2

Page 3 of 37

Environmental Science & Technology

27 28

Abstract

29

phycosphere and the resulting modulations of contaminant speciation and uptake is poorly

30

characterized. Here we modeled the effect of algal C and N uptake on carbonate cycling and

31

speciation of selected contaminants in the phycosphere (external boundary layer) of

32

chrysophytes, a key phytoplankton group in oligotrophic systems. We calculated enrichments

33

in H+ concentration relative to that in the bulk solution (pH 7.0) of approximately 40% for

34

NH4+-grown cells or depletions of approximately 500% for NO3--grown cells at the algal

35

membrane surface of a 5-µm or 30-µm radius cell or colony, respectively. Such changes are

36

mainly due to direct H+ uptake or release at the plasmalemma if NO3- or NH4+ is the N source,

37

respectively. Due to these pH changes in the boundary layer, competition between H+ and

38

metals for uptake is enhanced, which contributes to a decrease in potential metal uptake. Our

39

model suggests that the uptake of protonated weakly acidic organic acids (HA) is greater in

40

NH4+-grown cells compared to NO3--grown cells. The account of chemical reactions in the

41

algal boundary layer could improve ecological risk assessments for a wide range of

42

contaminants.

The interplay between biological and chemical reactions in the freshwater phytoplankton

43

ACS Paragon Plus Environment

3

Environmental Science & Technology

44

Page 4 of 37

Abstract Art efflux

HA

uptake

NO3- OHMez+

efflux

HA

uptake

Reactive boundary layer Mez+

NH4+ H+

45 46 47

1 Introduction

48

Chrysophyte species in the Chrysophyceae and the closely-related class Synurophyceae

49

constitute a dominant or sub-dominant portion of phytoplankton biomass in lakes, with

50

moderate to low productivity.1 With diverse carbon acquisition mechanisms, ranging from

51

photoautotrophy in all Synurophyceae to phagomixotrophy in many Chrysophyceae, these

52

algal classes play an important role in determining the carbon flow through the food web,

53

including in Canadian Shield water bodies that represent ~50% of the world’s freshwater

54

lakes.2 Several freshwater chrysophyte habitats are impacted by metals released from mining

55

industries3,4 and organic contaminants discharged from urban areas.5,6 Anthropogenic

56

eutrophication of freshwater systems as well as global change related-CO2 increases and

57

concomitant acidification7 can also affect the physiology and diversity of chrysophytes.8 Such

58

pH changes may further impact the uptake and toxicity of several inorganic and hydrophobic

59

contaminants in phytoplankton.9-11 For freshwater phytoplankton, chemical reactions in the

60

phycosphere (external boundary layer) surrounding cells rather than those in the well-mixed

61

bulk phase dictate the overall interactions between algae and their external aquatic medium.

62

However, due to practical and technical limitations, in-depth microscale analysis of these

63

chemical interactions has been proven difficult and has received little attention.12,13

ACS Paragon Plus Environment

4

Page 5 of 37

Environmental Science & Technology

64

A major uncertainty is how photosynthesis affects carbonate chemistry in the immediate

65

vicinity of freshwater algal cells, and what the consequences are in terms of speciation and

66

uptake of nutrients and contaminants. This issue remain difficult to evaluate for cyanobacteria

67

and algae because of their carbon concentrating mechanisms (CCM), which increases algal

68

photosynthesis at low inorganic carbon concentrations.14 In contrast, chrysophytes lack such a

69

CCM and hence soley rely on diffusive CO2 transport for photosynthesis.15-18 They are

70

therefore an attractive model for the analysis of carbonate system and associated contaminant

71

disequilibria in the critical reactive layer adjacent to the algal cell surface.

72

Whereas modeling work and experiments have shown that marine phytoplankton can

73

substantially deplete carbon dioxide at the cell surface, which is particularly true for large

74

phytoplankton species,19-21 little information is available for phytoplankton species growing in

75

freshwater environments, many of which are poorly buffered.22,23 A better knowledge of the

76

ecophysiology of CCM-lacking algae is needed to assess the impacts of global climate change

77

and contaminant pollution on organic carbon production by freshwater algae. Such a coupling

78

between algal cell physiology and ecotoxicology is of increasing interest especially given the

79

ongoing development of metal toxicity models that account for physiological regulation of

80

metal transport systems,24-26 and given the increasing recognition of the importance of algal

81

physiological state on pesticide toxicity.27,28

82

The aim of the present study is to model and evaluate the effect of photosynthesis on

83

carbonate chemistry within the boundary layer of chrysophyte cells or colonies of different

84

sizes. The proposed model simulates reactive carbon transfers with an explicit account of

85

species diffusion and intertwined acid/base disequilibria as coupled to C and N nutrient

86

uptake, and the associated effect on the uptake of inorganic and organic contaminants.

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 37

87 88

2. Theory

89

2.1. Modeling the influence of algal metabolism on the carbonate system in the boundary

90

layer

91

For strictly phototrophic cells grown under continuous light conditions as experienced at high

92

latitudes in summer, adopting a cell size scaling derived for algal growth from laboratory

93

experiments, the net carbon fixation rate (d[CO2] /dt) can be written as the product of the

94

carbon cell quota (Qc, pgC/cell) and the specific growth rate (µ):

95

d[CO2]/dt = Qc µ

1

96 97

The cell carbon quota in freshwater phytoplankton grown under continuous light scales

98

linearly with cell volume with a proportionality coefficient of 0.279 pg C µm-3 according to

99

the following equation:29

100

Qc = 0.279 Vc

2

101 102

where Vc is the cell volume (µm3). This relationship is valid for cells with biovolumes

103

between ~2 and 100,000 µm3, i.e., a cell radius in the range 1.3 to 30 µm assuming spherical

104

geometry for chrysophytes such as Dinobryon and Mallomonas.30 However, Dinobryon

105

generally occur as loose colonies of cells unlike the chrysophyte Synura, which forms more

106

compact colonies.

107

The nutritional mode of several freshwater phytoplankton species lacking CCM varies widely

108

from strictly phototrophic to phagomixotrophic so that cells may acquire carbon not only

109

from photosynthesis but also by the phagotrophy variant of heterotrophy. The current study

ACS Paragon Plus Environment

6

Page 7 of 37

Environmental Science & Technology

110

considers spherical cells without CCM growing photoautotrophically without or with (see

111

below) phagocytosis, which is important in several protist species. Photosynthesis has been

112

shown to provide more than 75% of cell carbon in the mixotroph Dinobryon31 whereas

113

several other chrysophytes are strict (e.g., Mallomonas and Synura) or possibly facultative

114

(Hydrurus) photoautotrophs.15

115

Based on eq. 2, the mean cellular organic carbon concentration in a range of freshwater algae

116

growing in nutrient-sufficient conditions is 2.3 × 104 mol C per m3 of cell biovolume. This is

117

a good approximation of the cellular carbon concentration of 2.7 × 104 mol C m-3 in

118

Dinobryon divergens 32 and of 2.6 × 104 mol C m-3 in Dinobryon cylindrum .31

119

Adopting a theoretical power law relating cell volume to growth rate (µ = constant Vc-0.25) (33

120

and references therein), which is at the lower end of the range of values of the exponent for

121

nano- and micro-phytoplankton,34 and a growth rate of 0.4 d-1 as measured for Dinobryon

122

(equivalent spherical radius = 4 µm) in an axenic culture,31 the specific growth rate (µ in d-1)

123

as a function of cell radius (rc in m) can be evaluated from:

124

µ = (2.90 x 10-5)rc-0.75

3

125 126

It should be noted that this theoretical cell size-scaling relationship minimizes growth rates

127

(C and N nutrient uptake rates) and hence boundary layer perturbations in the carbonate

128

system and contaminant speciation. This choice adds robustness to our modeling approach.

129 130

Combining eq.1 for the net CO2 uptake flux (mol CO2 m-2 d-1) by an alga with eq. 3, and

131

assuming a spherical cell with a mean intracellular C concentration per cell surface area of 2.3

132

× 104 mol C per m3 (see calculation above), with rc in m, it comes:

ACS Paragon Plus Environment

7

Environmental Science & Technology

133

d[CO2]/dt = (3.16 x 10-6) rc0.25

Page 8 of 37

4

134 135

The carbon uptake rate (d[CO2]/dt) is assumed to be independent of bulk-solution pH from 5

136

to 8, even though photosynthesis in Mallomonas sp. has been shown to vary by ~30%

137

between pH 5 and 7 and by 3- to 4-fold between pH 6 and 8.18 The effect of pH on carbon

138

uptake is specifically considered in the discussion section.

139

Calculations were also performed considering mixotrophic growth based on data from

140

Dinobryon cultures with added bacteria. This chrysophyte can grow at 0.7 d-1 and acquire

141

75% or more of its internal cell C through photosynthesis, with the remaining carbon

142

provided by ingestion of bacteria.31,35 Further assuming an identical C:N ratio in bacteria and

143

algae of 116:16, it follows that ≤ 25% of algal cell N can be provided by ingestion of bacteria.

144

Since the cellular fate of bacterial N and its associated coupling with H+/OH- exchange at the

145

plasmalemma is unknown to our knowledge, no model corrections were made regarding the

146

impact of bacterial N uptake on H+/OH- exchange at the plasmalemma. We further considered

147

the amount of bacterial carbon respired as CO2 by assuming a 1.05 mean molar ratio of CO2

148

production to organic C taken up through heterotrophy in the mixotrophic chrysophyte

149

Ochromonas malhamensis.32 Therefore, in the case of mixotrophic growth, the net CO2

150

uptake rate for a given cell radius decreased by 25% due to a decrease in photosynthesis, and

151

by 25% due to bacterial organic carbon respiration, but increased by 75% due to an increase

152

in growth rate relative to the phototrophic case. As a result, the net CO2 uptake rate (mol m-2

153

d-1) for mixotrophic growth was only 14% lower than that for a strict phototroph. This uptake

154

can thus be expressed as a function of rc according to:

155

d[CO2]/dt = (2.77 x 10-6)rc0.25

5

ACS Paragon Plus Environment

8

Page 9 of 37

Environmental Science & Technology

156 157

2.2 Modeling the carbonate distribution in the algal boundary layer

158

We modeled the effect of pH (and carbonate speciation) in the bulk solution on the pH and

159

carbonate speciation prevailing in the algal boundary layer. To that end and for simplicity, we

160

considered that the bulk-solution dissolved CO2 concentration was in equilibrium with the

161

atmosphere, as an approximation of carbonate system for several oligotrophic freshwater

162

lakes. We varied the bulk-solution pH from 5 to 8, encompassing the range of pH values

163

typically measured in freshwater systems from acidic Canadian Shield lakes to lakes with

164

limestone in the catchment basin. We assumed that the bulk-solution concentrations of all

165

inorganic carbon species and the bulk-solution pH did not change over time as expected under

166

non-bloom, low algal cell density conditions. The steady-state CO2, HCO3-, CO32-, OH- and

167

H+ concentration gradients in the external boundary layer were evaluated as a function of

168

cell/colony size. All calculations were performed at 1 mol m-3 ionic strength (I), and model

169

outcomes were evaluated to verify that the model changed by less than 30% percent over a

170

range of I applicable in freshwater systems (from 0.5 to 5 mol m-3). The rate of change in

171

concentration of each chemical species in the boundary layer (dc/dt) was further given by the

172

following diffusion-reaction equation:20,36

173

 

 + = 

6

174 175

 is the diffusion Laplacian operator in spherical geometry and R is the sum of all where 

176

reaction rates involving the chemical species of interest. The diffusion rate of each inorganic

177

carbon species is defined by:

178

 c = D

  

 

 (     )

7

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 37

179 180

where D and [c] are the diffusion coefficient and concentration of a given chemical species,

181

respectively, and r is the radial distance from the center of the spherical cell. A common value

182

for D of 1.18 x 10-9 m2 s-1 was adopted for all carbonate species to maintain electroneutrality

183

(SI section 1). The extent of disequilibrium under steady-state diffusion conditions in the

184

carbonate system due to algal nutrient C uptake was estimated taking into account the rate

185

constants of the relevant chemical reactions involving CO2, HCO3-, CO32- and H+/OH- in the

186

external boundary layer. Direct H+ release or uptake at the plasmalemma due to NH4+ or NO3-

187

uptake, respectively, was further considered (0.15 mol H+ per mol C taken up, assuming a 1:1

188

N:H+ co-transport and a Redfield C:N ratio).37 At steady-state, the time rate of change in the

189

concentration of inorganic carbon species and H+/OH- are zero, and eq. 6 can be solved for

190

the concentration of CO2, HCO3- and CO32-, H+ and OH- as a function of the distance from the

191

cell surface using two boundary conditions: 1) the equilibrium concentration of chemical

192

species in the bulk solution beyond the external boundary layer, which is determined from the

193

atmospheric CO2 pressure and bulk-solution pH, and 2) the uptake flux of CO2, which is the

194

only inorganic carbon species taken up by the cell. Note that CO2 uptake is not coupled to net

195

H+ transport at the cell membrane, since CO2 assimilation produces neutral carbohydrates. A

196

proton is taken up or released, however, at the cell membrane due to NO3- and NH4+

197

assimilation, respectively (see section 2.1). All chemical equations and modeling procedures

198

used to determine R in eq. 6, are detailed in the SI, section 2.

199 200

2.3 Modeling metal speciation in the boundary layer as a function of pH

ACS Paragon Plus Environment

10

Page 11 of 37

Environmental Science & Technology

201

To test the effect of external boundary-layer pH on metal uptake, we focussed on cadmium

202

Cd(II), which has been extensively studied in ecotoxicology laboratory experiments using

203

well-defined inorganic media. Cd(II) was chosen for its high toxicity10 as well as its relatively

204

simple speciation, with a large fraction of Cd(II) present as the free, hydrated ion or as

205

inorganic complexes at pHs greater than 7.5; typically less than 25% of total Cd is complexed

206

to dissolved organic matter (DOM) in low DOM freshwater environments that are often

207

associated with oligotrophic conditions.4,38 Cadmium speciation was evaluated in the external

208

algal boundary layer considering a total Cd concentration of 1 × 10-6 mol m-3 and

209

complexation with OH- and CO32- (SI, section 3). It was further assumed that DOM was not

210

present in the external boundary layer for simplification in order to mimic experimental media

211

routinely used in short-term metal uptake experiments. Cadmium uptake rates by freshwater

212

phytoplankton are generally much slower (ca. by one order of magnitude) than metal

213

diffusion from the bulk solution to the algal membrane surface. Therefore, Cd uptake is not

214

expected to be limited by diffusion, but rather by the internalisation step particularly at

215

nanomolar Cd concentrations, provided that Cd uptake kinetics in green algae is comparable

216

to that in chrysophytes.10,24 Furthermore, since metal internalisation kinetics are much slower

217

than that for other chemical species considered in this work (e.g., CO2), algal Cd uptake was

218

ignored by setting the corresponding flux at the algal surface to zero as a boundary condition;

219

this is detailed in the SI section 3. As a second boundary condition, it was assumed that Cd

220

speciation was at equilibrium in the bulk solution.

221 222

2.4 Modeling the speciation of ionisable weak acid/base organic contaminants in the

223

boundary layer

ACS Paragon Plus Environment

11

Environmental Science & Technology

Page 12 of 37

224

Protonation/deprotonation of weak acid/base organic contaminants is important from an

225

ecotoxicological viewpoint, since ionic species do not freely pass through the cell membrane

226

and therefore are generally much less bioavailable (or toxic) for living cells than neutral

227

species that are able to cross algal membranes via simple diffusion.23 Here we considered two

228

examples of ionisable organic contaminants whose uptake may not be limited by diffusion in

229

the external boundary layer but rather may be controlled by changes in their speciation.

230

Contaminants were selected in such a way that their modeled membrane permeability (Pm)

231

was much smaller than the modeled boundary layer permeability (Pbl). The Pm term was

232

calculated using the Xiang and Anderson empirical equations39 and Pbl was determined by

233

dividing the diffusion coefficient of the organic contaminant by the external boundary layer

234

thickness (see SI section 4.1).

235

As a first example, we calculated the boundary layer speciation of the polar herbicide maleic

236

hydrazide (MA) (Kow = 0.72, pKa = 5.64), which contaminates freshwater systems and may

237

affect phytoplankton.40 Secondly, we estimated the speciation of a mildly hydrophobic

238

antifungal pharmaceutical clindamycin (Clind) (pKa = 7.79, Kow = 101.83), which is released in

239

freshwater and is of potential concern for phytoplankton.41 For each case, the calculated

240

exchange rate constant of the neutral and charged species with H+ were considered. It should

241

be noted that the proportion of the neutral species relative to the charged contaminant

242

concentration dominates at pH < pKa or pH > pKa, for MA or Clind, respectively. Two

243

boundary conditions were considered: 1) organic contaminant uptake was set to 0 since the

244

uptake is limited by the membrane transport process; and 2) organic contamination speciation

245

in the bulk solution is at equilibrium. The modeling procedure for these contaminants is

246

further detailed in SI section 4.

ACS Paragon Plus Environment

12

Page 13 of 37

Environmental Science & Technology

247 248

2.5 Model assumptions

249

Since the thickness of the cell wall in unicellular algae is less than 10% of the thickness of the

250

effective external boundary layer in small to giant-celled algae, diffusion of inorganic carbon

251

species and H+ were assumed to be controlled by the external boundary layer,23 and thus the

252

potential effect of the cell wall on diffusive fluxes was not considered in the model.

253

Therefore, throughout this study, enrichment at the cell or colony surface refers to the nearest

254

point of the external boundary layer at the surface of the algal plasmalemma. In all

255

calculations, it was assumed that dissolved CO2 and N concentrations sustained maximal algal

256

growth. We also assumed that C and N uptake rates were not significantly affected by the

257

presence of contaminants as expected at low contaminant concentrations or under short-term

258

exposure. Algal cells or colonies were assumed spherical and interstitial spaces within loose

259

colony of cells were not considered for simplicity. The influence of turbulence was neglected

260

at the micrometer scale (see SI section 5)42 as was the effect of cell swimming and sinking

261

rates on diffusion in the boundary layer; potential influences of turbulence, swimming and

262

gravitational sinking are discussed in SI section 6. The role of heterotrophic bacteria in CO2

263

production within the algal boundary layer was also not considered in the model based on

264

calculations performed in SI section 7. The role of other processes that may regulate cellular

265

acid/base reactions (e.g., production of cell intermediary metabolites, nucleic acids and

266

proteins as well as organic acid excretion) are neglected, since these processes are not

267

constrained for chrysophytes. The role of cell-surface charge on the enrichment of ions was

268

also not considered because of the insignificant or small differences in ionic bio-partitioning

269

reported for Chlorophyceae with increasing pH (see SI section 8 and discussions for further

ACS Paragon Plus Environment

13

Environmental Science & Technology

Page 14 of 37

270

details).43,44 Reactions involving reactive oxygen species were not considered due to

271

insufficient knowledge regarding these reactions in the external boundary layer. Finally, since

272

the steady-state concentration gradient and chemical equilibria for CO2 and other chemical

273

species are achieved within milliseconds to minutes, much shorter than algal doubling times,

274

any effect of cell-number changes on model results due to algal growth were not included (see

275

SI section 9).

276 277

2.6 Model implementation

278

All simulations were run and all figures plotted using R software (v 1.0.153).45 All scripts

279

including functions are freely available on github (DOI: 10.5281/zenodo.1252222). The R-

280

package “Reactran” (v 1.4.3.1) was used to solve the differential diffusion-reaction equation

281

(eq. 6) in spherical coordinates using the cell radius and bulk-solution pH as model input

282

variables, and a set of model parameters (C and N uptake rates, equilibrium and rate

283

constants, ionic strength and bulk contaminant speciation).46 Steady-state concentrations of

284

carbonate species or contaminants were computed in the external boundary layer using

285

unidimensional model grids over a distance from the algal cell membrane of 60 µm (for cell

286

radius of 5 µm) or 100 µm (for cell or colony radius of 30-um) with a model grid resolution

287

of 10,000 elements.

288 289

3. RESULTS AND DISCUSSION

290

3.1 Carbonate system disequilibrium and associated pH gradients in the boundary layer

291

as a function of nitrogen source

ACS Paragon Plus Environment

14

Page 15 of 37

Environmental Science & Technology

292

Photosynthesis in model chrysophyte cells or colonies of 5 or 30-µm radius led to significant

293

gradients in CO2, H+/OH-, HCO3-, CO32- in the external boundary layer (Fig. 1 and 2). Slow

294

CO2 hydration or dehydration kinetics relative to diffusion generated a CO2 drawdown in the

295

boundary layer; N uptake altered H+/OH- concentrations at the membrane surface (Fig. 1 and

296

2 and SI section 9). At pH 7, for a 5 µm-radius chrysophyte cell (representative of small

297

Chromulina or Mallomonas), algal metabolism resulted in an increase or decrease of cell-

298

surface [H+] by around 30-40% for cells grown on NH4+ or NO3-, respectively (Fig 1 and 2).

299

This change was almost entirely accounted for by direct H+/OH- exchange at the

300

plasmalemma due to N uptake and not a decrease in carbonate system buffering due to CO2

301

uptake (SI section 10).

302

The proton gradient was even more pronounced for larger cells. The H+ concentration at the

303

surface of a 30-µm radius cell/colony was depleted by approximately 5-fold if NO3- was the N

304

source, whereas a 3.5-fold cell-surface enrichment in H+ was obtained for the same

305

cell/colony size when NH4+ was the N source (Fig. 1 and 2). However, cell surface depletion

306

in H+ due to C uptake (assuming no inorganic N uptake) only reached ~25% for a 30 µm-

307

radius cell/colony despite a decrease in cell surface CO2 concentrations of more than 3.5-fold,

308

indicating that HCO3- buffering remained high (SI section 10). Indeed, over the range of cell

309

radii tested between 5 and 30 µm, and at pH 7, HCO3- concentrations changed by less than

310

1%, and therefore, as for the 5-µm radius cells, the proton gradient was mostly due to nitrogen

311

uptake (and direct H+ exchange at the cell membrane) rather than to a change in carbonate

312

buffering. Uptake of N influences more the cell surface pH than CO2 uptake since N uptake

313

and assimilation, unlike CO2 uptake, is coupled to net H+ transport at the cell membrane,

314

which directly affects cell surface pH, and this effect may be significant depending on

ACS Paragon Plus Environment

15

Environmental Science & Technology

Page 16 of 37

315

boundary layer pH buffering, algal cell size and N uptake rate. By contrast, CO2 uptake

316

increases the cell surface pH indirectly through H+ consumption via thereaction (H+ + HCO3-

317

-> CO2 + H2O, see SI section 2), a reaction that is slow relative to diffusion leading to kinetic

318

limitation of CO2 replenishment at the cell surface and minimizing the associated pH rise. As

319

stated previously, whereas bicarbonate concentrations did not change, carbonate

320

concentrations did increase, but concentrations of this species were quite low and it therefore

321

only played a minor role (ca.