Physiological responses of Skeletonema costatum to the interactions of seawater acidification and combination of photoperiod and temperature

Ocean acidification (OA), which is a major environmental change caused by increasing atmospheric CO2, has considerable influences on marine phytoplankton. But few studies have investigated interactions of OA and seasonal changes in temperature and photoperiod on marine diatoms. In the present study, a marine diatom Skeletonema costatum was cultured under two different CO2 levels (LC, 400 μatm; HC, 1000 μatm) and three different combinations of temperature and photoperiod length (8:16 L:D with 5 °C, 12:12 L:D with 15 °C, 16:8 L:D with 25 °C), simulating different seasons in typical 15 temperate oceans, to investigate the combined effects of these factors. The results showed that specific growth rate of S. costatum increased with increasing temperature and daylength. However, OA showed contrasting effects on growth and photosynthesis under different combinations of temperature and daylength: while positive effects of OA were observed under spring and autumn conditions, it significantly decreased growth (11 %) and photosynthesis (21 %) in winter. In addition, low temperature and short daylength decreased the proteins of PSII (D1, CP47 and RbcL) at ambient pCO2 level, while OA 20 alleviated the negative effect. These data indicated that future ocean acidification may show differential effects on diatoms in different cluster of other factors.

Abstract. Ocean acidification (OA), which is a major environmental change caused by increasing atmospheric CO 2 , has considerable influences on marine phytoplankton. But few studies have investigated interactions of OA and seasonal changes in temperature and photoperiod on marine diatoms. In the present study, a marine diatom Skeletonema costatum was cultured under two different CO 2 levels (LC, 400 µatm; HC, 1000 µatm) and three different combinations of temperature and photoperiod length (8 : 16 L : D with 5 • C, 12 : 12 L : D with 15 • C, 16 : 8 L : D with 25 • C), simulating different seasons in typical temperate oceans, to investigate the combined effects of these factors. The results showed that specific growth rate of S. costatum increased with increasing temperature and day length. However, OA showed contrasting effects on growth and photosynthesis under different combinations of temperature and day length: while positive effects of OA were observed under spring and autumn conditions, it significantly decreased growth (11 %) and photosynthesis (21 %) in winter. In addition, OA alleviated the negative effect of low temperature and short day length on the abundance of RbcL and key photosystem II (PSII) proteins (D1 and D2). These data indicated that future ocean acidification may show differential effects on diatoms in different clusters of other factors.

Introduction
Ocean acidification (OA) is one of major environmental changes caused by increasing atmospheric CO 2 , which has directly increased from 280 ppm in the preindustrial era to higher than 400 ppm at present (Friedlingstein et al., 2019). It is predicted that surface seawater pH would drop 0.3-0.5 and 0.5-0.7 units by the years 2100 and 2300, respectively (Caldeira and Wickett, 2003). It has been suggested that calcifying organisms, such as coral reefs and coccolithophores, are vulnerable to OA due to the decreased calcification at elevated CO 2 (Albright et al., 2016). The responses of noncalcifying organisms such as diatoms to OA vary widely among taxonomic groups and may be detrimental, negligible or even beneficial (Gao and Campbell, 2014). Consequently, the abundance of marine phytoplankton and community structure might be altered by OA (Gattuso et al., 2015).
Diatoms are ubiquitous photosynthetic phytoplankton which account for about 20 % of global primary productivity and thus play a crucial role in the global cycling of carbon and silicon (Falkowski et al., 2004). To overcome the limited aqueous CO 2 concentration in seawater, they have developed CO 2 -concentrating mechanisms (CCMs) (Spalding, 2007). Decreased photosynthetic affinity for dissolved inorganic carbon (DIC) and activity of CCM-related enzymes are generally found under increased CO 2 conditions (Raven and Beardall, 2014). For phytoplankton assemblages, elevated CO 2 could lead to increases in chlorophyll a concentrations and the abundance of diatoms (Johnson et al., 2013). Species and strain specificity are observed in studies on physiological responses of diatoms to OA, which might be caused by the balance between positive effects of elevated CO 2 and negative effects of decreased pH (Langer et al., 2009;Li et al., 2016). In addition, acclimation and adaptation processes, i.e., the timescale of diatoms exposed to OA, could also influence the physiological effects of OA (Wu et al., 2014;Li et al., 2017). Moreover, other environmental factors, such as temperature (Seebah et al., 2014), light (Gao et al., 2012), nutrients (W. Hartmann et al., 2013) and clusters of multiple factors (Xu et al., 2014a, b), are shown to have interaction with OA on diatoms.
Diatoms are widespread across oceans; thus they would experience different photoperiods. Photoperiod controls the total light dose received by phytoplankton and thus could remarkably influence the physiological performance such as growth and lipid content of microalgae (Wahidin et al., 2013). For Antarctic sea ice microalgae Chlamydomonas sp., continuous illumination stimulates higher growth and nutrient absorption rates than successive darkness conditions (Xu et al., 2014c). Growth rate of Chlamydomonas reinhardtii is gradually enhanced following an increasing photoperiod (Hsieh et al., 2017). In contrast, Alexandrium minutum grows faster under short day length relative to longer and even continuous day length (Wang et al., 2019). Moreover, different photoperiods could influence intracellular carbon demand of microalgae, which has a stronger regulation effect on CCMs compared with effects of changes in CO 2 supply (Rost et al., 2006).
Under the combined influence of photoperiod and OA, physiological performance of phytoplankton might be different from that under a single factor. For example, continuous light moderates the negative effect of OA on coccolithophore growth, although species isolated from different regions show diverse responses (Bretherton et al., 2019). The changes of photoperiod are often accompanied by an increase or decrease in temperature, and impacts of OA on diatoms can also be changed by temperature. For example, under OA condition, decreased metabolic activity is observed in Phaeodactylum tricornutum when temperature is elevated (Bautista-Chamizo et al., 2018), while elevated temperature enhances the growth rate of Nitzschia lecointei (Torstensson et al., 2013).
But limited studies have investigated interactions between OA and the combination of temperature and photoperiod (i.e., conditions in different seasons) on diatoms. Skeletonema costatum is a widespread eurythermal and euryhaline diatom species, which frequently causes red tide. We hypothesized the effect of OA on S. costatum may be modulated by photoperiod and temperature. In the present study, we investigate the physiological performance of marine diatom Skeletonema costatum under two different CO 2 levels and three combinations of temperature and photoperiod, which simulated different seasons in typical temperate oceans (win-ter, 5 • C with 8 : 16 L : D; spring or autumn, 15 • C with 12 : 12 L : D; summer, 25 • C with 16 : 8 L : D).

Culture conditions
The diatom Skeletonema costatum in this study was isolated from Gaogong Island, Lianyungang, Jiangsu Province (34 • 70 74.95 N, 119 • 49 26.47 E). Before being used in experiments, the cells were cultured in autoclaved natural seawater enriched with f/2 medium (Guillard and Ryther, 1962). Semi-continuous cultures were maintained in 500 mL Erlenmeyer flasks with a filter unit (Millex-GP, Merck, USA) in order to aerate sterile air. Triplicate independent cultures were set for each treatment at the light intensity of 150 µmol photons m −2 s −1 .

Experimental setup
In order to evaluate effects of pCO 2 levels and different combinations of temperature and photoperiod on S. costatum, cells were cultured under winter (5 • C with light: dark cycle of 8 : 16 h), spring-autumn (15 • C with 12 : 12 h) and summer (25 • C with 16 : 8 h) conditions independently with two pCO 2 levels (400 ppm, LC; 1000 ppm, HC), simulating temperature and day length conditions of different seasons in typical temperate oceans. Temperatures and light intensity (150 µmol photons m −2 s −1 ) were controlled by illumination incubators (GXZ-500B, Ningbo, China). Cells were inoculated in cultures with fresh medium which was aerated with ambient air (400 ppm) or CO 2 -enriched air (1000 ppm). The high pCO 2 level was manipulated by a CO 2 plant incubator (HP 1000 G-D, Ruihua Instruments, Wuhan, China). Cultures were kept at the exponential phase by diluting every 3 d, and cell concentrations were controlled below 2 × 10 5 cell mL −1 in order to minimize the effect of cell metabolism on carbonate chemistry in the medium. The changes in culture pH were less than 0.05 during the 3 d (8.10 ± 0.01 for LC and 7.85 ± 0.01 for HC in winter; 8.14±0.01 for LC and 7.85±0.01 for HC in spring-autumn; 8.19±0.02 for LC and 7.89±0.02 for HC in summer). After acclimating to different treatments for at least 40 generations, the following parameters were measured.

Growth measurement
To estimate the growth of S. costatum, triplicate samples (1 mL each) were collected from each treatment at 48 and 72 h after dilution and fixed with 10 µL Lugol's solution; then a plankton counting chamber (DSJ-01, Xundeng Instruments, Xiamen, China) was used to count cells directly under an optical microscope (DM500, Leica, Germany). The specific growth rate was calculated as µ = (ln N t − ln N 0 )/(t − t 0 ), where N t represents the cell concentration (cells mL −1 ) at time t; N 0 represents the cell concentration at time t 0 , t − t 0 = 1 d. The growth rates were averaged from three dilution processes within each growth condition.

Chlorophyll a and BSi measurements
Samples were filtered onto GF/F filters (25 mm, Whatman, UK), and chlorophyll a was extracted with 4 mL of methanol at 4 • C for 24 h in darkness. An ultraviolet spectrophotometer (Ultrospect 3300 pro, Amersham Bioscience, Sweden) was used to detect the absorption values of supernatant under 632, 665 and 750 nm after centrifuging (Biofuge primo R, Thermo, Germany). The chlorophyll a concentration (pg per cell) of S. costatum was calculated by the equation of Ritchie (2006).
Samples (200 mL) for biogenic silica (BSi) measurement (pmol per cell) were filtered onto polycarbonate filters (0.8 µm, Merck Millipore, Germany) by a polysulfone filter funnel (25 mm, Pall Corporation, UK), and filters were then dried at 80 • C for 24 h. BSi on the filter was digested by 4 mL of 0.2 M NaOH in a boiling bath for 40 min and was neutralized with 1 mL of 1 M HCl when cooled. The supernatant (1 mL) was diluted with 4 mL of Milli-Q water, and then 2 mL of molybdate solution and 3 mL of reducing agent were added into tubes. The absorption was measured at 810 nm by an ultraviolet spectrophotometer (Ultrospect 3300 pro, Amersham Bioscience, Sweden) after the color developed for 2-3 h (Brzezinski and Nelson, 1995).

Photosynthesis and respiration measurements
The net photosynthetic rate under culture conditions and photosynthetic oxygen evolution rate vs. light intensity (P -I ) curve of S. costatum were measured through a Clark-type oxygen electrode (Oxygraph+, Hansatech, UK), in which temperature was controlled by a thermostatic water bath (DHX-2005, China).
For measurement of net photosynthesis under culture conditions, light intensity was set as 150 µmol photons m −2 s −1 (light intensity during culture), provided by a halogen lamp (QVF135, Philips, Netherlands). A sample of 50 mL was filtered (< 0.02 MPa) onto a cellulose acetate membrane (Xinya Instruments, Shanghai, China). Cells were filtered and resuspended in 5 mL pre-aerated fresh medium under cultured conditions, which was then used to determine oxygen evolution rate and cell concentration. Oxygen consumption was measured under darkness which was realized by covering the reaction chamber with an opaque box.

P -I curve measurement
For the P -I curve, oxygen consumption rates in darkness and net photosynthetic oxygen evolution at seven different light intensities (0, 10, 20, 50, 100, 200, 500, 1000 µmol photons m −2 s −1 ) were identified. Light intensity was achieved by adjusting the distance between the halo-gen lamp and oxygen electrode chamber. Cells were also filtered and resuspended in pre-aerated fresh medium in a similar way to photosynthesis measurement under culture conditions. Photosynthetic rate and light intensity data were fitted according to Henley (1993): P = P m ×tanh(α ×PAR/P m )+ R d , where PAR is irradiance, P is photosynthetic rate, P m is light-saturated photosynthetic rate, α is initial slope of P -I curve and R d is dark respiration rate. I k (saturating irradiance for photosynthesis) and I c (light compensation point) were also calculated by I k = P m /α, I c = R d /α.

Chlorophyll fluorescence measurement
Cells were concentrated by gentle vacuum filtration (< 0.02 MPa) for measurement of rapid light curves (RLCs) under eight different PAR levels (0, 10, 20, 50, 100, 200, 500, 1000 µmol photons m −2 s −1 ) lasting for 10 s each using a handheld fluorometer (AquaPen-C AP-P 100, Chech). Relative electron transport rates (rETRs) of S. costatum were measured and were estimated as in the study of Wu et al. (2010): rETR = PAR × Y (II) × 0.5, where PAR represents the photon flux density of actinic light, Y (II) represents the effective quantum yield of PSII and 0.5 is based on the assumption that PSII receives half of all absorbed quanta. RLCs were fitted as P = PAR/(a × PAR 2 + b × PAR + c), where P represents rETR; a, b and c are model parameters. The relative photoinhibition ratio of rETR was calculated as Inh(%) = (rETR max − rETR x )/rETR max × 100 %, where rETR x is the value of rETR at 1000 µmol photons m −2 s −1 .

Protein measurements
The abundance of ribulose-1,5-bisphosphate carboxylase/oxygenase (RubisCO), large subunit binding protein (RbcL) and PSII proteins (PsbA (D1), PsbD (D2) and PsbB (CP47)) was measured under different pCO 2 levels and combinations of temperature and day length. D1 and D2 proteins are located in the reaction center, while CP47 is a junction of the antenna, and RbcL is a component of RubisCO which is a key enzyme in the CO 2 fixation process. For the relative value of protein measurements, cells were filtered and resuspended in 2 mL of extracting medium (50 mM Tris-HCl, pH 7.6, 5.0 mM MgCl 2 , 10 mM NaCl, 0.4 M sucrose and 0.1 % BSA) according to Ma et al. (2020). After cell disruption and centrifugation, supernatant liquid was used to measure chlorophyll concentration (µg mL −1 ) according to Arnon (1949): C = D 652 × 1000/34.5 × T , where C represents total chlorophyll concentration, D 652 represents the absorption value at 652 nm, T represents the dilution ratio, and 1000 and 34.5 are constants. The same concentration of chlorophyll (2.4 micrograms) per lane was used for 12 % sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE, Mini PROTEAN, Bio-rad, America) at 150 V for 1 h, and the proteins were transferred into polyvinylidene difluoride (PVDF) membranes which were then immersed in blocking solution with antibodies (D1, D2, CP47, RbcL and Actin; Agrisera) for 1 h, and successively goat anti-rabbit secondary antibodies were used. Blots were developed by using enhanced chemiluminescence luminescence (ECL) reagent and were quantified with a chemiluminescence detection system (Tanon 5500, Shanghai, China). Actin was used as the internal control in order to correct the experimental error in the process of quantitative sample loading of protein, to ensure the accuracy of the experimental results. And the data provide us with a general trend, not accurate concentrations, among different treatments.

Data analyses
Data were analyzed with IBM SPSS Statistics 24 and are presented as mean ±SD (standard deviation). One-way ANOVA was used to compare differences among combinations of temperature and photoperiod treatments. The independentsample t test was applied to compare differences between two pCO 2 levels. A general linear model was used to assess the interactive effects of CO 2 level and combination of temperature and photoperiod on growth rate, rETR, photosynthesis and respiration, contents of chlorophyll a and BSi, and proteins. When P values were under 0.05, the Tukey test was used for post hoc analysis.

Specific growth rate
The growth rate of S. costatum ranged from 0.47 ± 0.01 to 3.22 ± 0.08 d −1 under different treatments and increased significantly with increasing temperature and day length (P < 0.05) regardless of the pCO 2 level (Fig. 1). In the summer season, elevated pCO 2 showed no significant effects; however, it remarkably influenced the growth rate in other seasons. Elevated pCO 2 enhanced the growth rate by 11 % in spring and autumn (P < 0.001), while the reverse pattern was found in winter (P < 0.001). The general linear model indicated that the season and CO 2 level had a notable interaction on specific growth rate (P < 0.001, Table 1).

Chlorophyll a and BSi contents
Under ambient CO 2 conditions, chlorophyll a content was enhanced by increased temperature and day length (Table 2), and the content was 22 % higher in summer compared with winter (P = 0.008). When CO 2 was elevated, chlorophyll a content in winter was 42 % and 32 % lower than that in spring and summer, respectively (P = 0.001, P = 0.004). Elevated pCO 2 decreased chlorophyll a content in winter (P = 0.022) while enhancing it in spring (P = 0.002) and had no significant impact in summer. A significant interaction between season and CO 2 can be found (P < 0.001) ( Table 1). Figure 1. Specific growth rate of S. costatum acclimated to ambient (LC, black bars) and elevated pCO 2 (HC, white bars) under different combination of temperature and photoperiod conditions. The data are mean ±SD values of triplicate cultures (n = 3). Asterisks represent significant differences (P < 0.05) between two CO 2 levels under the same season conditions (t test).
A different trend was detected for BSi content (Table 2). Under ambient pCO 2 , BSi content decreased with higher temperature along with longer day length, and the value in winter was significantly higher than that in spring and summer (P = 0.005, 0.002, respectively). Higher pCO 2 decreased BSi significantly in winter (P = 0.016) and spring (P = 0.007), while it had no significant influence on the content in summer (P = 0.3). There is a significant interaction between season and CO 2 on BSi content (P < 0.05) (Table 1).

Photosynthesis and respiration
Net photosynthetic oxygen evolution and dark respiration rates showed similar patterns under the same CO 2 conditions (Fig. 2). The lowest photosynthesis and respiration rates were observed under winter conditions, and maximal rates were observed in summer at each pCO 2 level. Both photosynthesis and respiration rates increased with increasing temperature and day length (P < 0.05). Elevated pCO 2 inhibited net photosynthetic rate under winter conditions (P = 0.0013) while photosynthesis was enhanced by elevated pCO 2 in spring and autumn (P = 0.006). In addition, higher pCO 2 stimulated dark respiration rate in spring and autumn (P < 0.001). Both photosynthesis and respiration were not impacted by higher pCO 2 in summer. Interaction between season and CO 2 on net photosynthetic rate was detected (P = 0.035). Positive relationships of dark respiration or net photosynthetic rates and growth rate were observed ( Fig. 4a and b). Table 1. Significance test results of growth, BSi, photosynthetic and respiration rates, parameters of P -I (P max is the maximum net photosynthetic rate, α is the photosynthetic efficiency, R d is the dark respiration rate), and RLC curves (rETR max , which is the maximum relative electron rate) for a combination of temperature and photoperiod (season), CO 2 , and their interactions.

P -I curve
Net photosynthetic oxygen evolution rates increased with increasing light intensity initially and then reached the plateaus in all seasons; the curves under winter conditions reached the plateaus much earlier than the other two seasons (Fig. 3). Higher temperature and prolonged day length had a main effect on P max , R d , I k and I c . However, elevated pCO 2 instead of season had the main effect on α (Table 1). P max was enhanced when temperature increased with prolonged day length (P < 0.05) except when the summer season was compared with spring and autumn conditions at elevated pCO 2 ( Table 3). Effects of elevated pCO 2 were only observed in spring and autumn (P = 0.03). I k increased remarkably in summer under both pCO 2 treatments, which was similar to R d at elevated pCO 2 (P < 0.05). There was a significant interaction between CO 2 and the combination of temperature and photoperiod on R d (Table 1). At higher pCO 2 levels, rETR max values were significantly different among different seasons (P < 0.05), and the highest value was found in spring and autumn (52 % and 14 % higher than winter and summer, Table 3). Elevated pCO 2 decreased rETR max in winter and summer (P < 0.001 in winter, P = 0.01 in summer). rETR max had a positive correlation with growth rate when temperature and day length increased from winter to spring-autumn conditions. However, as temperature and day length continuously increased from spring-autumn to summer, the correlation became negative (Fig. 4c). Interactions between the two factors were detected on rETR max (P < 0.001, Table 1). Photoinhibitions were found in RLCs of cells under all treatments. As the season proceeded from winter to summer, photoinhibitions were alleviated significantly at both pCO 2 levels (P < 0.05). In spring and autumn, the inhibition at higher pCO 2 was significantly decreased compared with ambient pCO 2 conditions (P = 0.039, Table 3).

PSII protein concentrations
The relative values of RbcL and key PSII proteins (PsbA (D1), PsbD (D2) and PsbB (CP47)) were quantified in different seasons under ambient and elevated CO 2 conditions (Fig. 5). At the ambient pCO 2 level, the highest contents of all four proteins were detected in spring and autumn. Elevated pCO 2 significantly enhanced RbcL, D1 and D2 protein contents in winter (P < 0.05). In addition, higher pCO 2 led to an increase in CP47 in summer (P = 0.006).

Discussion
Phytoplankton like diatoms have already evolved several strategies to cope with different temperatures and day lengths in temperate oceans, where variations in season are evident. However, the ongoing elevated pCO 2 combined with Table 3. Photosynthetic parameters fitted from P (O 2 ) -I and rapid light curves for S. costatum acclimated to ambient and elevated pCO 2 in different seasons. P max (pmol O 2 per cell per hour) is the maximum net photosynthetic rate, α is the photosynthetic efficiency, R d (pmol O 2 per cell per hour) is the dark respiration rate, I k (µmol photons m −2 s −1 ) is the photosynthetic saturated light intensity, I c (µmol photons m −2 s −1 ) is light compensation point, rETR max is the maximum relative electron rate and Inh is the relative photoinhibition ratio of rETR. Different superscript letters represent significant differences (P < 0.05) between two CO 2 levels under the same season (t test).

Physiological responses of S. costatum to different combinations of temperature and photoperiod
In the present study, the growth rate of S. costatum increased with increasing temperature and day length regardless of the pCO 2 level (Fig. 1). Previous studies showed that most phytoplankton, such as Chlamydomonas reinhardtii, Trichodesmium or Alexandrium catenella, grew faster under prolonged photoperiods (Cai and Gao, 2015), although Alexandrium minutum grew faster under shorter photoperiods (Wang et al., 2019). For S. costatum, a remarkably higher contribution of HCO − 3 to the overall carbon uptake was ob- served under light-dark cycles compared with continuous light, and a shorter photoperiod led to lower photosynthetic affinity for inorganic carbon (Rost et al., 2006). Basically, the chl a quota in microalgae increases with decreasing day length; however, our results exhibited an inverse pattern (Table 2). The inconsistency might be caused by the different temperatures set in studies, which is another main environmental factor affecting the growth of diatoms. Increasing temperature may lead to various changes in growth rate depending on whether the temperature is optimal for the species. For S. costatum, its growth rate has been shown to increase with temperature up to 30 • C (Ebrahimi and Salarzadeh, 2016). Zhang et al. (2020) also found that the growth rate of S. costatum could increase with temperature from 5 to 30 • C and then drop sharply. The underlying mechanism is that elevated temperature promotes S. costatum metabolic rates when nutrients are abundant. This could be shown by the relationship between respiration and growth. Higher mitochondrial respiration can result in higher growth rate theoretically (Geider and Osborne, 1989), since this process provides ATP and carbon skeletons (Raven et al., 2017). It seems that temperature shows the dominant effect compared with day length.
RubisCO is an important enzyme for carbon fixation; psychrophilic diatoms utilize increasing abundance of RbcL protein to maintain high cellular enzymatic rates and growth rate Figure 5. (a) Immunoblot using antibodies against RbcL and key PSII proteins (PsbA (D1), PsbD (D2) and PsbB (CP47)) isolated from S. costatum acclimated to ambient and elevated pCO 2 in different season levels (W for winter, A for autumn, S for summer). Each line was loaded with similar amounts of proteins. (b) Quantitative analysis of proteins. The relative abundance of each band was estimated by densitometric scanning of the exposed films. Asterisks represent significant differences (P < 0.05) between two CO 2 levels under same season (t test). at low temperature (Young et al., 2015). However, for diatoms in temperate areas, the amount of RbcL protein decreased in low temperatures along with short day lengths compared with higher temperatures and longer day lengths in our results. Low abundance of RbcL was in line with slower growth rate in winter conditions.

Effect of ocean acidification under different seasons
Our results showed that the impacts of elevated pCO 2 on S. costatum depended on seasonal changes in temperature and photoperiod. High pCO 2 had enhanced the growth of S. costatum in spring-autumn and reduced it in winter, while no significant effects were detected in summer (Fig. 1). CO 2 concentration mechanisms (CCMs) are energy-dependent, and high pCO 2 downregulates CCMs of most phytoplankton including S. costatum, so the saved energy could be used for growth (Raven et al., 2017). Higher initial slope of the P -I curve at elevated pCO 2 might be partly responsible for the higher growth rate compared with that at ambient pCO 2 (Table 3). But in winter, growth decreased under OA conditions, although respiration and P max in P -I curves had no significant changes. This is because the combination of biochemical and biophysical CCMs may cause the lack of a positive response to elevated pCO 2 under near-optimal growth conditions (Passow and Laws, 2015). In addition, when other environmental factors are stressful, the sensitivity of diatoms to CO 2 and temperature is prominent (Taucher et al., 2015). The shorter day length and low temperature simulating winter conditions in the present study can be seen as stressors, under which S. costatum was more sensitive to elevated pCO 2 . When temperature and photoperiod are optimal, positive or neutral effects of higher pCO 2 were observed. However, different patterns were reported for other species, such as E. huxleyi and the macroalgae Ulva linza (Bretherton et al., 2019;Yue et al., 2019). For these two species, reduced growth rates at elevated pCO 2 were found when the day length was longest. High temperature might accelerate nutrient uptake and metabolic rates, which may alleviate the negative effects of longer day length under a higher pCO 2 environment (Bretherton et al., 2019;Yue et al., 2019). The maximum photosynthetic rate increased significantly under higher pCO 2 in spring and autumn conditions. This was consistent with the higher photosynthetic efficiency and growth rate (Table 3, Fig. 1).
Silicification directly relates to cell division and growth and is independent of photoperiod (Brzezinski, 1992). BSi contents generally increased with decreasing growth rate when any limiting factors such as temperature, light or ammonium exist (Martin et al., 2000). In the present study, elevated pCO 2 mitigated the negative effects of temperature and photoperiod limitation on BSi content (Table 2). Higher BSi contents in winter under ambient CO 2 conditions can intensify the ballasting effects and thus impact the sinking rate of organic matter produced by diatoms. Lavaud et al. (2016) indicated that PSII activity and phosphorylation of thylakoid protein may play a crucial role in controlling the change of the photosynthetic activity. The contents of D1 and D2 proteins decreased in winter because of the negative effect of lower temperature (Mock and Valentin, 2004), while elevated pCO 2 level increased their contents (Fig. 5). The protein contents included both photochemically active PSII and PSII that is inactivated but retains D1 and D2 subunits (G. . Proteins will be degraded and synthesized rapidly after damage. However, the degraded rate of photosynthetic proteins from photoinactivated PSII complexes could be different according to culture conditions. For example, the removal rate of D1 protein increased with growth light in the diatom Thalassiosira pseudonana (Campbell et al., 2013). Therefore, the photosynthetic rate might decouple with the contents of proteins. The decline in growth under winter conditions might result from the increased metabolic costs of photoprotection and elevated D1 turnover under the combination of short day length limitation and low temperature (Hoppe et al., 2015). Although RbcL decreased with elevated CO 2 in some studies that might be caused by the decrease in RuBP concentration for diatoms (Endo et al., 2015), McCarthy et al. (2012) observed an increase in RubisCO concentration with higher CO 2 , which is in line with the present study under winter conditions. Light and temperature could affect RbcL amount in phytoplankton. RbcL increased slightly with higher CO 2 at low-light-intensity conditions, while it decreased slightly at higher light intensity in the research of Levitan et al. (2010). In the present study, the combination of low temperature and short day length might lead to a complex trend of RbcL: elevated CO 2 might stimulate the sensitivity of S. costatum to low temperature, which leads to a steep rise of RbcL in winter.
For diatoms, an increasing number of studies have paid attention to the effect of interactions of ocean acidification and other environmental factors such as light intensity, UV, temperature, nutrient limitation, salinity or photoperiod (Gao et al., 2012;Yue et al., 2019), and the results showed positive, negative or neutral effects (Xu et al., 2014c;Li et al., 2018). However, few studies combined elevated pCO 2 with seasonal changes in seawater physical and chemical characters on marine diatoms. In this study, temperature and photoperiod were chosen as seasonal factors to investigate the combined effects of pCO 2 and these two factors on physiology of S. costatum. Our results suggested temperature and photoperiod could mediate effects of elevated pCO 2 on the typical diatom S. costatum. Positive effects of OA on growth and photosynthesis were observed in spring and autumn, while negative effects were found in winter conditions. To better understand how global climate changes would affect marine diatoms in the future, it is necessary to explore the interactive effects of ocean acidification with seasonal changes in seawater characters.
Data availability. All data included in this study are available upon request via contact with the corresponding author.
Author contributions. JX and FL conceived and designed the experiments. HL and TX carried out the experiments. HL, FL and JX analyzed data. HL wrote the draft of the paper. FL, JM and JX revised the manuscript and approved this version for submission.