The Arctic picoeukaryote Micromonas pusilla benefits synergistically from warming and ocean acidification

In the Arctic Ocean, climate change effects such as warming and ocean acidification (OA) are manifesting faster than in other regions. Yet, we are lacking a mechanistic understanding of the interactive effects of these drivers on Arctic primary producers. In the current study, one of the most abundant species of the Arctic Ocean, the prasinophyte Micromonas pusilla, was exposed to a range of different pCO2 levels at two temperatures representing realistic current and future scenarios for nutrient-replete conditions. We observed that warming and OA synergistically increased growth rates at intermediate to high pCO2 levels. Furthermore, elevated temperatures shifted the pCO2 optimum of biomass production to higher levels. Based on changes in cellular composition and photophysiology, we hypothesise that the observed synergies can be explained by beneficial effects of warming on carbon fixation in combination with facilitated carbon acquisition under OA. Our findings help to understand the higher abundances of picoeukaryotes such as M. pusilla under OA, as has been observed in many mesocosm studies.


Introduction
With the progress in using molecular tools to describe marine biodiversity in recent decades, the scientific community has become increasingly aware of the underestimated importance of picoeukaryotes, for both primary and export production of the world's oceans (Richardson and Jackson, 2007;Worden and Not, 2008).Larger phytoplankton such as diatoms are efficient vectors for carbon export due to aggregate formation and ingestion by large zooplankton, leading to the production of fast-settling faecal pellets (Sherr et al., 2003).In contrast, picoeukaryotes are mainly grazed by smaller heterotrophic protists such as ciliates, which have a low carbon retention, excrete relatively more dissolved material and thus fuel recycled production (Sherr and Sherr, 2002).Hence, changes in the relative abundance of pico-and nanoeukaryotes can have large implications for food webs and biogeochemistry (Worden et al., 2015).
Picoeukaryotes tend to dominate low nutrient environments, which is often attributed to their high surface : volume ratios and mixotrophic capacities (Raven, 1998;McKie-Krisberg and Sanders, 2014).The low nutrient concentrations in the Arctic surface ocean, for example, cause picoeukaryotes to be particularly successful in this region.In fact, the globally occurring prasinophyte Micromonas pusilla is considered the most abundant species in the Arctic ocean (Šlapeta et al., 2006;Lovejoy et al., 2007;Marquardt et al., 2016).In this environment, strong stratification causes low nutrient concentrations throughout the summer and autumn months (Tremblay et al., 2015), and the occurrence of the polar night requires organisms to either form resting stages or to have heterotrophic capacities (Tremblay et al., 2009;Lovejoy, 2014;Berge et al., 2015;Vader et al., 2015).
Climate change effects manifest faster in the Arctic than anywhere else on the planet (Stocker, 2014).In this region, for example, temperatures are rising more than twice as fast as throughout the rest of the globe (Miller et al., 2010).The concurrent rapid reduction in ice cover allows for more light penetration and longer growing seasons, while Published by Copernicus Publications on behalf of the European Geosciences Union.
increased stratification due to ice melt and warming constrain nutrient supply to surface waters, both of which will change the dynamics of primary production (Arrigo et al., 2008;Wassmann and Reigstad, 2011).Ocean acidification (OA) is also especially pronounced in the Arctic Ocean, because low temperatures and alkalinity make the system sensitive to anthropogenic CO 2 loading (AMAP, 2013;Qi et al., 2017).Picoeukaryotes such as M. pusilla may benefit from these changes and are considered potential winners of climate change.In the Canadian Arctic, for example, picoeukaryote abundances are increasing as surface waters get warmer, fresher and more oligotrophic (Li et al., 2009).Regarding OA effects, the majority of studies on natural phytoplankton assemblages have shown that picoeukaryotes, particularly M. pusilla, increase in relative abundance with increasing pCO 2 levels (Engel et al., 2008;Meakin and Wyman, 2011;Newbold et al., 2012;Brussaard et al., 2013;Hussherr et al., 2017;Schulz et al., 2017).Despite the evident sensitivity of M. pusilla to changes in pCO 2 levels, a detailed assessment of the OA effects, their interaction with warming as well as the underlying mechanisms in this important species is still missing.
Like all photosynthetic organisms, cells of M. pusilla need to maintain a balance between energy sources (i.e.light harvesting by the photosynthetic apparatus) and sinks (most importantly carbon fixation in the Calvin cycle) to prevent harmful levels of excitation pressure on the photosynthetic electron transport chain (Behrenfeld et al., 2008).Light harvesting and electron transport in the photosystems are largely independent of changes in temperature and pCO 2 (Mock and Hoch, 2005;Hoppe et al., 2015), but the impact of these drivers on energy sinks can potentially affect the energy balance of the cell: the beneficial effects of elevated pCO 2 observed in phytoplankton are thought to be caused by increased diffusive CO 2 supply, reduced CO 2 leakage or lowered costs to operate their CO 2 concentrating mechanisms (Rost et al., 2008;Bach et al., 2013).Elevated temperatures, however, can change enzyme kinetics, including those involved in the Calvin cycle, thus leading to a larger sink of excitation energy (Maxwell et al., 1994;Toseland et al., 2013).Hence, both ocean warming and acidification potentially increase the efficiency of photosynthesis and biomass production, at least up to the organisms' respective optimum levels.Above these levels, temperatures and proton concentrations start to disrupt enzymatic processes, increase the investment into pH homeostasis, and impair the delicate regulation of cellular processes (Levitt, 1980;Taylor et al., 2001;Flynn et al., 2012).Thus, the complex balance between beneficial and detrimental effects will determine whether the combination of warming and OA will synergistically promote or deteriorate phytoplankton growth and biomass build-up.
In the current study, we aim to investigate the responses of an Arctic M. pusilla strain to warming and OA.To this end, M. pusilla was grown at four pCO 2 levels ranging from preindustrial to future scenarios (180-1400 µatm) under 2 and 6 • C, which represent the magnitude of the projected future temperature increase in this region (Collins et al., 2013), but also the current spring and summer temperatures in the environment where the strain was isolated (Hegseth et al., 2018).

Culture conditions
Monoclonal cultures of the picoeukaryote Micromonas pusilla (Butcher) Manton & Parke (isolated in 2014 by Klara Wolf in Kongsfjorden, Svalbard, 79 • N; local temperature range −1.5 to 8 • C) taxonomic identification confirmed by rDNA sequencing of SSU, LSU and ITS sequences) were grown in 1 L glass bottles in semi-continuous dilute-batch cultures (maximum 129 000 cells mL −1 ; diluted every 3-4 days) under constant irradiances of 150 ± 26 µmol photons m −2 s −1 .Media consisted of 0.2 µm sterile-filtered Arctic seawater with a salinity of 32.7 enriched with macronutrients, trace metals and vitamins according to F /2 R medium (Guillard and Ryther, 1962).Light intensities were provided by daylight lamps (Philips Master TL-D 18 W; emission peaks at wavelengths of 440, 560 and 635 nm), adjusted by neutral density screens and monitored using a LI-1400 data logger (Li-Cor) equipped with a 4π sensor (Walz).Cells were growing at four different CO 2 partial pressures (pCO 2 ; 180, 380, 1000 and 1400 µatm) and two temperatures (2.2 ± 0.3 and 6.3 ± 0.2 • C).Cultures were acclimated to these conditions for at least seven generations prior to sampling.
Different pCO 2 conditions were achieved by aeration of the incubation bottles with air of the respective pCO 2 levels delivered through sterile 0.2 µm air filters (Midisart 2000, Sartorius Stedim) for 24 h prior to inoculation.Gas mixtures were generated using a gas flow controller (CGM 2000 MCZ Umwelttechnik), in which CO 2 -free air (< 1 ppmv CO 2 ; Dominick Hunter) was mixed with pure CO 2 (Air Liquide Deutschland).The pCO 2 levels in the gas mixtures were regularly monitored with a non-dispersive infrared analyzer system (LI6252, LI-COR Biosciences), calibrated with CO 2 -free air and purchased gas mixtures of 150 ± 10 and 1000 ± 20 ppmv CO 2 (Air Liquide Deutschland).

Carbonate chemistry
Samples for total alkalinity (A T ) were filtered through 0.7 µm glass fibre filters (GF/F, Whatman) and stored in borosilicate bottles at 3 • C. A T was estimated from duplicate potentiometric titration (Brewer et al., 1986) using a TitroLine alpha plus (Schott Instruments).A T values were corrected for systematic errors based on measurements of certified reference materials (CRMs provided by Andrew Dickson, Scripps, USA; batch no.111; reproducibility ± 5 µmol kg −1 ).Total dissolved inorganic carbon (C T ) samples were filtered through 0.2 µm cellulose-acetate filters (Sartorius Ste-dim) and stored in gas-tight borosilicate bottles at 3 • C. C T was measured colorimetrically in triplicates with a QuAAtro autoanalyzer (Seal; Stoll et al., 2001).The analyser was calibrated with NaHCO 3 solutions (with a salinity of 35, achieved by addition of NaCl) to achieve concentrations ranging from 1800 to 2300 µmol C T kg −1 .CRMs were used for corrections of errors in instrument performance such as baseline drifts (reproducibility ± 8 µmol kg −1 ).Seawater pH total was measured potentiometrically with a two-point calibrated glass reference electrode (IOline, Schott Instruments).An internal TRIS-based reference standard (Dickson et al., 2007) was used to correct for variability on electrode performance (reproducibility ± 0.015 pH units).Following recommendations by Hoppe et al. (2012), seawater carbonate chemistry including pCO 2 was calculated from A T and pH using CO 2SYS (Pierrot et al., 2006).The dissociation constants of carbonic acid of Mehrbach et al. (1973), as refitted by Dickson and Millero (1987), were used for calculations.Dissociation constants for KHSO 4 were taken from Dickson (1990).

Growth, elemental composition and production rates
Samples for cell counts were fixed with glutaraldehyde (0.5 % final concentration).After gentle mixing, samples were stored at room temperature in the dark for 15 min, and subsequently frozen in liquid nitrogen and stored at −80 • C. Prior to analysis, samples were thawed on ice and mixed thoroughly.After addition of 10 µL SybrGreen working solution (dissolved in DMSO) and 10 µL YG beads working solution (1 µm Fluoresbrite calibration beads grade YG, Polyscience), samples were counted on an Accuri C6 flow cytometer (BD Biosciences) equipped with a blue solid-state laser (488 nm excitation wavelength) run on medium fluidics settings (35 µL min −1 ; 16 µm core size) with a limit of 50 000 events or 250 µL.Analysis was performed based on red (FL3 channel, > 670 nm) and green (FL1 channel, 533 ± 30 nm) fluorescence, as well as sideward and forward light scattering.Specific growth rates constants (µ) were determined from exponential fits of cell counts over 4 consecutive days.Particulate organic carbon (POC) and nitrogen (PON) were measured after filtration onto precombusted (15 h, 500 • C) GF/F filters (Whatman), stored at −20 • C and dried for at least 12 h at 60 • C prior to sample preparation.Analysis was performed using a CHNS-O elemental analyser (Euro EA 3000, HEKAtech).Contents of POC and PON were corrected for blank measurements and normalised to filtered volume and cell densities to yield cellular quotas.Production rates of POC were calculated by multiplying the cellular quota with the division rate constant k of the respective incubation.Samples for determination of chlorophyll a (Chl a) were filtered onto GF/F filters (Whatman), immediately placed into liquid nitrogen and stored at −80 • C until analysis.Chl a was subsequently extracted in 8 mL 90 % acetone at 4 • C over night.Chl a concentrations were determined on a fluorometer (TD-700, Turner Designs), using an acidification step (1 M HCl) to determine phaeopigments (Knap et al., 1996).

Variable Chl a fluorescence
Photophysiological characteristics, based on photosystem II (PSII) variable Chl a fluorescence, were measured using a fast repetition rate fluorometer (FRRf; FastOcean PTX, Chelsea Technologies) in combination with a Fas-tAct Laboratory system (Chelsea Technologies).The excitation wavelength of the fluorometer's light-emitting diodes (LEDs) was 450 nm, and the applied light intensity was 21 587 µmol photons m −2 s −1 .The FRRf was used in single turnover mode, with a saturation phase comprising 100 flashlets on a 2 µs pitch and a relaxation phase comprising 40 flashlets on a 50 µs pitch.Measurements from all replicates (n = 3) were conducted in a temperature-controlled chamber (± 0.2 • C) at the respective treatment temperature.
After subtraction of a blank value, the minimum (F 0 and F 0 for light-and dark-acclimated measurements, respectively) and maximum Chl a fluorescence (F m , and F m for light-and dark-acclimated measurements, respectively) were estimated from iterative algorithms for the induction (Kolber et al., 1998) and relaxation phase (Oxborough, 2012) after 15 min of dark acclimation, which was sufficient to achieve a dark-acclimated state (data not shown).All fluorescence parameters were calculated by standard equations (Genty et al., 1989;Maxwell and Johnson, 2000).Maximum quantum yields of PSII (apparent PSII photochemical quantum efficiency; F v /F m ) were calculated as (1) Fluorescence-based photosynthesis-irradiance curves (PI) were conducted at six irradiances (I ) between 33 and 672 µmol photons m −2 s −1 , with an acclimation time of 10 min per light step.Electron transfer rate (ETR, measured in mol e − (mol RCII) −1 s −1 ) through PSII for each light step was calculated as follows: Following the suggestion by Silsbe and Kromkamp (2012), the light-use efficiency (α, mol e − m 2 (mol RCII) −1 (mol photons) −1 ) and the maximum electron transfer rates per RCII (ETR max , mol e − (mol RCII) −1 s −1 ) were estimated by fitting the data to the model by Webb et al. (1974): ]. (3) The light saturation index (E k , µmol photons m −2 s −1 ) was then calculated as ETR max / α.
where F q is the differences between measured and maximal fluorescence (Suggett et al., 2010).F 0 was measured after each light step (with a duration of 90 s).

Statistics
All data are given as the mean of three biological replicates with ± 1 standard deviation.To test for significant differences between the treatments, two-way analyses of variance (ANOVA) with additional normality (Kolmogorov-Smirnov) and post hoc (Holm-Šidák) tests were performed.The significance level was set to 0.05.Statistical analyses were performed with the program SigmaPlot (Systat Software Inc, Version 12.5).

Carbonate chemistry
Regular dilution of cultures with a pre-aerated seawater medium kept carbonate chemistry stable over the course of the experiment.More specifically, in each bottle the drift in A T and C T compared to initial values was ≤ 3 and ≤ 4 %, respectively (data not shown).Final carbonate chemistry in the 2 • C treatments yielded pCO 2 levels of 197 ± 3, 323 ± 12, 959 ± 22 and 1380 ± 53 µatm (Table 1).In the 6  Results from statistical analysis can be found in Table S2.
Molar C : N ratios of the biomass (Tables 2, S1) increased with temperature (two-way ANOVA, F = 14, p = 0.002), yet this overall difference was mainly driven by results at  S2.
The ratios of C : Chl a (Fig. 2d, Tables 2, S1) were elevated under high compared to low temperature conditions (two-way ANOVA, F = 14, p = 0.002), an effect that was most pronounced at pCO 2 levels of 180 µatm (post hoc test, t = 5.5, p < 0.001).While there was no effect of pCO 2 on C : Chl a at low temperature, C : Chl a decreased with increasing pCO 2 at high temperature (two-way ANOVA, interaction term, F = 6, p = 0.007; 180 vs. 1400 µatm at 6 • C post hoc test, t = 3.9, p = 0.008).

Chl a fluorescence-based photophysiology
The effects of the applied treatments on photophysiology were studied by means of FRRf, which investigates photochemistry at photosystem II (PSII).No effects of the applied treatments were observed in most parameters investigated (Tables 3, S1).This was true for the dark-acclimated quantum yield efficiency of PSII (F v /F m ), which was similar in all treatments with values of 0.45 ± 0.06, as well as for absorption cross sections of PSII light harvesting (σ PSII ).
Furthermore, the fitted parameters of FRRf-based PI curves (α, ETR max and E K ) were independent of the experimental treatments (Tables 3, S1).In contrast, the rate constant of the reopening of PSII reaction centres (τ ES ; Tables 3,  S1) was slightly yet significantly smaller under high temperatures (two-way ANOVA, F = 6, p = 0.029), even though this overall response also depended on the applied pCO 2 levels (two-way ANOVA, interaction term, F = 4, p = 0.033).

Micromonas pusilla benefits from warming
We observed a strong stimulation of growth rates and biomass build-up with increasing temperature (Fig. 1, Table 2).Even though the isolate stems from 1.8 • C water temperature, the beneficial effects of warming from 2 to 6 • C are not surprising: M. pusilla is known to dominate Arctic phytoplankton assemblages in the summer and autumn situations (Lovejoy et al., 2007;Marquardt et al., 2016) when surface temperatures of 6 • C or more can be reached (Hegseth et al., 2018).Our results are also in line with mesocosm experiments that indicate stimulatory effects of warming on picoplankton abundances (Daufresne et al., 2009;Sommer et al., 2015) as well as with the temperature optimum of 6-8 • C observed for another Arctic strain of M. pusilla (Lovejoy et al., 2007).
Table 3. FRR fluorometrical PSII photochemistry measurements -PSII quantum yield efficiency F v /F m (dimensionless), functional absorption cross section (σ PSII ) (nm −2 PSII −1 ), rate of PSII re-opening (τ ES , ms), maximum nonphotochemical quenching at 672 µmol photons m −2 s −1 (NPQ max , dimensionless), maximum light-use efficiency (initial slope α, mol e − m 2 (mol RCII) −1 (mol photons) −1 ), maximal absolute electron transfer rates through PSII (ETR max , mol e − (mol RCII) −1 s −1 ), and the light saturation index (E K , µmol photons m −2 s −1 ) under the different temperature and pCO 2 treatments (n = 3; mean ± 1 SD).Results from statistical analyses can be found in Table S2 Below the temperature optimum of a cell, warming causes an acceleration of the entire metabolism, as enzymatic reactions run faster under these conditions (Eppley, 1972;Brown et al., 2004).In this study, warming caused higher growth rates, POC quotas and biomass production (Fig. 2, Tables 2,  S1), indicating that particularly the fixation and storage of carbon was facilitated by increasing temperature.Electron transport processes, however, were largely independent of temperature (Tables 3, S1).Thus, temperature affected the balance between electron transport ("light reaction") and carbon fixation in the Calvin cycle ("dark reactions").Especially under relatively low temperatures, as investigated here, warming can decrease the excitation pressure on the electron transport chain of the photosystems by increasing the temperature-limited turnover rates of enzyme reactions such as RuBisCO (Mock and Hoch, 2005).Thus, cells grown under low temperatures need to invest relatively more energy into biosynthesis than into photochemistry compared to cells grown under high temperatures (Toseland et al., 2013).While it has been shown that Antarctic diatoms can compensate for slow RuBisCO kinetics by increasing the expression of this enzyme (Young et al., 2014), it is unknown whether such acclimation responses also occur in prasinophytes.Regarding the C : Chl a ratio, this can be taken as an indicator of how many resources the cell retains as carbon biomass (e.g.structural and storage compounds) relative to how much is invested into its light harvesting capacities (Halsey and Jones, 2015).In this study, the strong temperature-dependent increase in C : Chl a (Fig. 2, Table S1) under potentially limiting pCO 2 levels of 180 µatm suggests that under warming, the balance between light harvesting and carbon fixation was indeed more beneficial for biomass build-up.Furthermore, elevated temperature significantly decreased τ ES (Tables 3,  S1), which can serve as a proxy of the rate at which downstream processes can remove electrons from PSII (Kolber et al., 1998).Thus, our results indicate that the drainage of electrons into carbon fixation was faster under warmer conditions, explaining the higher growth and biomass production under these conditions.
4.2 Warming shifts CO 2 optima towards higher pCO 2 levels Under 6 • C and pCO 2 levels expected to be reached by the end of this century, OA had a significantly positive effect on growth and biomass build-up (Fig. 1).This finding is in line with previous studies, which have shown that picoeukaryotes can benefit strongly from OA in both laboratory and mesocosm studies (Meakin and Wyman, 2011;Newbold et al., 2012;Schaum et al., 2012;Brussaard et al., 2013;Maat et al., 2014;Schulz et al., 2017).Such positive response to OA could indicate that picoeukaryotes such as M. pusilla are mainly dependent on diffusive CO 2 supply and thus directly benefit from higher CO 2 concentrations (Brussaard et al., 2013;Schulz et al., 2013Schulz et al., , 2017)).While this could be related to the large surface-to-volume ratio of small cells, opposite trends within the group of diatoms (i.e. higher sensitivity of larger compared to smaller diatoms; Wu et al., 2014;Sett et al., 2018) suggest more complex underlying mechanisms at play.Despite this overall effect, growth rates of M. pusilla tended to follow a non-linear response curve over the tested range of glacial to elevated future pCO 2 levels (i.e. 180 to 1400 µatm), i.e. growth increased with increasing pCO 2 from low to intermediate, but decreased again under the highest pCO 2 levels (Fig. 1).Such an optimum behaviour can be expected for most environmental drivers (Harley et al., 2017) and has previously been observed in response to OA (Sett et al., 2014;Wolf et al., 2018).The response patterns in these studies were attributed to a combination of beneficial effects of rising pCO 2 under potentially carbon-limiting conditions for photosynthesis, and negative effects of declining pH on cellular homeostasis and enzyme performance, which man- ifest mainly at high pCO 2 (Bach et al., 2013).This nonlinearity in the observed pCO 2 effects emphasises the importance of experiments with more than two pCO 2 levels in order to properly describe OA-response patterns of organisms.
On a more general level, apparent discrepancies between OA studies can be attributed to actual differences in the environmental settings and their interactive effects with pCO 2 (Riebesell and Gattuso, 2015).When comparing the two most commonly applied pCO 2 levels, i.e. the present-day and the anticipated end-of-century situation, the effects of OA on most of the investigated physiological parameters are reversed under 6 compared to 2 • C (Fig. 3).This illustrates how difficult it is to infer responses to OA from experiments applying only one set of environmental conditions.It is also noteworthy that the combination of OA and warming led to more densely packed cells (no change in cell size based on flow cytometric measurements; data not shown) with similar stoichiometry compared to the control treatment (Table 2).This indicates that cells managed to cope well with the experienced future conditions.Furthermore, warming altered the OA-dependent change in most of the investigated parameters in a direction that indicates higher fitness compared to low temperatures (e.g. higher growth rates and higher elemental quota; Fig. 3; Hoppe et al., 2018b).Thus, the increase in growth under future compared to ambient conditions was larger than what would be expected by the respective responses to warming and OA in isolation, indicating synergistic beneficial effects of both drivers.

Potential mechanism underlying the interaction between warming and OA
The observed synergistic effects could be explained by their specific impacts on carbon acquisition and fixation.As outlined in the introduction, light and dark reaction of photosynthesis need to be balanced to achieve high biomass production while avoiding photodamage (Behrenfeld et al., 2008).According to our data, this balance is shifted towards higher biomass production rates under warming and OA.At higher temperatures, seawater CO 2 concentrations were lower than under colder conditions (Table 1; Zeebe and Wolf-Gladrow, 2001).At the same time, warming from 2 to 6 • C caused up to 60 % higher growth and 110 % higher biomass build-up rates (Fig. 1, Table 2).Furthermore, the decrease in τ ES indicates a faster transfer of photochemical energy into downstream processes such as RuBisCO activity (Table 3).Increased carbon demand in concert with lower carbon supply at higher temperatures thus increases the risk of CO 2 shortage in the cell, which in turn causes OA to have larger effects than at colder temperatures.Moreover, warming changes the kinetics of carbon fixation, with RuBisCO increasing its maximum turnover rates but decreasing its affinity for CO 2 (Young et al., 2014).At higher temperatures, cells thus have the potential for higher carboxylation rates provided sufficient CO 2 is available (Kranz et al., 2015).Under elevated pCO 2 levels, diffusive CO 2 supply increases and/or costs for active carbon acquisition decrease.Consequently, the positive effect of increasing temperature on the carbon fixation rate can develop its full potential under OA.
In conclusion, elevated catabolic activity under warmer conditions can explain the observed upward shift in the CO 2 optimum of growth with increasing temperature (Fig. 1), as the corresponding higher carbon demand causes CO 2 fixation to saturate under higher pCO 2 levels.In combination with a faster and more efficient machinery for pH homeostasis at elevated temperatures (Morgan-Kiss et al., 2006), this could explain why declining growth rates were only observed at relatively higher pCO 2 levels compared to those under low temperature conditions (Fig. 2).

Implications for the current and future Arctic pelagic ecosystem
Picoeukaryotes such as M. pusilla are considered to be potential winners of climate change: they are not only thriving in warmer, more stratified environments, which are predicted to further expand in the future, but also seem to benefit from OA (Li et al., 2009;Schulz et al., 2017).Our results for M. pusilla confirm beneficial effects of warming and OA on growth and biomass production under nutrientreplete conditions (Fig. 1, Table 2).Hence, in warmer spring conditions of the future, this species may experience growth stimulation under OA, potentially increasing its importance early in the growing season.Currently, M. pusilla already dominates Arctic phytoplankton assemblages in the nutrientlimited summer and autumn situations, which were not investigated in this study.Regarding the importance of nutrient availability, laboratory experiments found beneficial OA effects on M. pusilla primary production to persist also under P limitation (Maat et al., 2014), while in a mesocosm community, OA-dependent increases in M. pusilla abundances disappeared when the system ran into P and N co-limitation (Engel et al., 2008).Thus, it remains to be seen how the combined effects of warming and OA manifest under low nutrient conditions as well as how the responses may depend on sources and types of nutrients (e.g.mixing-delivered nitrate vs. regenerated ammonium).
A species' success in the environment does not only depend on individual performance, but also on how it compares to that of competing species.When we compare our results with the responses of the Arctic diatom Thalassiosira hyalina, isolated from the same location and exposed to the same experimental conditions (Wolf et al., 2018), the diatom had higher growth rates than the picoeukaryote under most treatment conditions, as can be expected for nutrient-replete conditions (Sarthou et al., 2005).The relative increase in growth rates from ambient (2 or 3 • C and 380 µatm pCO 2 ) to future conditions (6 • C and 1000 µatm pCO 2 ) was, however, much higher for M. pusilla than for T. hyalina.The fact that our were conducted under nutrientreplete conditions, which typically favour diatoms over picoeukaryotes, may indicate an even stronger increase in fitness (Collins et al., 2014) and could mean that M. pusilla gains another competitive advantage over phytoplankton like diatoms in the future, in addition to those resulting from changes in stratification (Li et al., 2009).Thus, our findings suggest higher picoplankton contribution to future Arctic phytoplankton assemblages under non-limiting conditions, e.g.early in the growing season when picoeukaryotes can already contribute quite substantially to the phytoplankton standing stocks (Marquardt et al., 2016;Paulsen et al., 2015).How such competition between diatoms and picoeukaryotes would manifest under nutrient-depleted conditions that strongly favour M. pusilla is currently unknown.
Even though picoeukaryotes seem to contribute more to the downward export of organic matter than previously assumed (Waite et al., 2000;Richardson and Jackson, 2007), in comparison to diatoms for example, they are less efficient vectors for carbon export to depth and have a lower energy transfer along trophic levels (Sherr et al., 2003).Consequently, Arctic food webs dominated by picoeukaryotes would look very different from those fuelled by diatom production (Sherr et al., 2003;Paulsen et al., 2015).Due to its motility and capability to grow mixotrophically, M. pusilla is characterized by an exceptionally high cellular C : N ratio compared to other Arctic phytoplankton (Table 2; Halsey et al., 2014;McKie-Krisberg and Sanders, 2014).An increased importance of this species would thus not only affect the food web due to its small size and concurrent grazer pref-erences, but also in terms of food quality (van de Waal and Boersma, 2012).The expected higher growth rates and thus abundances of this species may thus strengthen the Arctic microbial food web.Together with a concurrent weakening of the classical diatom-fuelled food web, this could have severe implications for the flow of energy and nutrients through future marine Arctic ecosystems (Post, 2016).

Conclusions
This study is the first to show synergistic effects of warming and OA on M. pusilla, one of the most abundant species of the worlds' oceans.Individually, both warming and OA cause more efficient biomass build-up under nutrient-replete conditions.Beneficial effects manifest, however, even more strongly in combination, when facilitated carbon acquisition (e.g.due to higher diffusive CO 2 supply) co-occurs with higher fixation rates (e.g.due to higher turnover-rates of Ru-BisCO).Our results provide an explanation for the observations of previous mesocosm studies, which indicated beneficial effects of OA and warming on M. pusilla and other picoeukaryotes.Characterising the responses of this key Arctic species to warming and OA will help to develop mechanistic phytoplankton functional types and more realistic model representation of phytoplankton assemblages as well as their responses to multiple drivers.Future studies are needed to elucidate further multifactorial environmental changes, addressing both abiotic (e.g.changes in light and nutrients) as well as biotic (e.g.heterotrophy, competition, grazers, viruses) interactions.
Data availability.The raw data of this study are available in the PANGAEA data repository with the identifier https://doi.org/10.1594/PANGAEA.892370(Hoppe et al., 2018a).
Author contributions.CJMH and BR designed the study.CJMH and CMF conducted the experiment.CJMH analysed the data and prepared the paper with contributions from BR and CMF.
Competing interests.The authors declare that they have no conflict of interest.identification.Furthermore, Laura Wischnewski, Anja Terbrüggen and Marcel Machnik are acknowledged for their help with sample analyses.
The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.
Edited by: Katja Fennel Reviewed by: Douglas Campbell and two anonymous referees

Figure 1 .
Figure 1.Specific growth rate constant µ (a) and POC production (b) of M. pusilla under low (filled symbols) and high (open symbols) temperatures as a function of pCO 2 (n = 3; mean ± 1 SD).Results from statistical analysis can be found in TableS2.

Figure 2 .
Figure 2. Cellular composition, i.e.POC (a), PON (b) and Chl a quota (c) as well as C : Chl a ratios (d), of M. pusilla under low (filled symbols) and high (open symbols) temperatures as a function of pCO 2 (n = 3; mean ± 1 SD).Results from statistical analysis can be found in TableS2.

Figure 3 .
Figure3.Schematic illustration of results for both temperatures over the entire range of pCO 2 levels as well as focusing on the responses between 380 and 1000 µatm (as the representation for commonly used OA treatments) and their modulation by temperature.

Table 1 .
(Pierrot et al., 2006)mistry at the end of the experiments (n = 3; mean ± 1 SD).CO 2 partial pressure (pCO 2 ) and dissolved CO 2 concentrations were calculated from total alkalinity (A T ) and pH total at 2 or 6 • C, a salinity of 32.7 using CO 2 SYS(Pierrot et al., 2006), and phosphate and silicate concentrations of 10 and 100 µmol kg −1 , respectively.NA indicates that values are not available for this specific treatment.

Table 2 .
Growth rate constants µ, division rate constants k, POC production rates and cellular quota of Chl a, and POC and PON as well as their ratios of M. pusilla at the end of the experiment under the different treatment conditions (n = 3; mean ± 1 SD).Results from statistical analysis can be found in TableS2.
• C treatments, pCO 2 levels were 198 ± 6, 394 ± 10, 1036 ± 31 and 1449 ± 18 µatm.Please note that the same pCO 2 level translates into differing dissolved CO 2 concentrations at different temperatures due to the temperature dependency of the carbonate system.Specifically, the treatment pCO 2 values translated into up to 13 % lower dissolved CO 2 concentrations in the 6 • C compared to the 2 • C treatment (Table1; see Fig.S1in the Supplement).Concurrently, the pCO 2 levels at 2 • C corresponded to pH total values of 8.30 ± 0.01, 8.11 ± 0.01, 7.68 ± 0.01 and 7.52 ± 0.02, respectively.In the 6 • C treatment, pH total values of the four pCO 2 treatments .