Introduction

Reproduction is an energetically expensive process. Successful reproduction requires investment of time and energy sometimes in detriment of other essential functions that maintain organisms’ fitness1. For organisms with multiple reproduction events, sacrificing reproduction over other processes, through prioritization of energy toward vital processes and functions, can increase their chances of survival in a stressing environment. However, disregarding reproduction can put at risk the long-term survival of the species2. In the two-spotted goby, Gobiusculus flavescens, used in this study, reproduction success includes physiological3 and behavioural4 processes. Females can allocate energy for partner selection, courtship, egg production, and offspring development5, while males energy can be allocated for partner selection, courtship, nest construction and defence6,7. The ability of a couple to successfully mate and produce clutches, as well as their potential investment in taking care of the clutches (parental care), is determinant aspects for their offspring’s fitness4.

Environmental factors can constraint reproduction of marine organisms8,9. Projections for the end of the century suggest an average increase in temperature and pCO2 in surface ocean up to + 2.5 °C and + 500µatm, respectively10. The increase of pCO2 in the water leads to a process known as ocean acidification, a perturbation of the carbonate chemistry including a decrease in pH. However, the situation is more complex in the coastal zone where natural and anthropogenic events, such as upwelling11 and eutrophication12, lead to high variability in temperature and pCO213.

Changes in temperature, one of the most important cues for the onset of the reproductive season, can be responsible for delaying or even interrupting reproductive processes14. Fluctuations in the reproductive timings, can in turn, cause a mismatch effect between offspring, food availability and predators15. For example, elevated temperatures can disrupt the metabolic activity and foraging behaviours of predators16,17 and preys18. pH has also been shown to impact reproduction in fish. Some studies showed that reproductive outputs can be up to 50% higher at lower pH, deviating from present environmental pH variability3,19,20. Other studies show reduced reproductive output in fish exposed to lower pH, as well as eggs’ abnormalities20,21. Impairing reproduction can have severe effects on the offspring quality, however, increasing the investment in reproduction might require using energy that otherwise would be allocated to maintain the homeostasis of the adult organisms22,23, and therefore make them more vulnerable to these stressors.

Studies looking at the combined impact of temperature and pCO2 have showed complex responses depending on the level of exposure24, life stages and habitat25,26,27,28,29. To date, there are only a few studies on the combined effect of temperature and pH on fish reproductive behaviour and performance. For example, in the tropical fish Amphiprion melanopus, clutch survival decreased at higher temperature, with no survival when combined with low pH treatments30.

The present study aims at evaluating the effect of temperature and pH on the reproductive performance of a temperate marine fish species, the two-spotted goby, G. flavescens. This semi-pelagic species inhabits shallow rocky waters of the European coasts. In Northern Europe, it plays a key ecological role as food of other important commercial fish species, such as Gadus morhua31. This species has sexual dimorphism and parental care during the embryonic phase, as the males defend the nest and attend to the eggs6. The two-spotted goby is also an excellent experimental model. It is easy to maintain and breed in captivity, with minimum stress for the fish3.

More specifically, this study evaluates the effect of temperature and pH on (1) the number of clutches laid and successfully hatched; (2) the size and number of eggs per clutch; (3) the larval length at hatch; (4) the parental care behaviour, namely fanning, cleaning, defending and chasing behaviours in the nest; and (5) hepatosomatic and gonadosomatic indexes of the breeding couples, as a proxy for parental investment.

Results

Seawater chemistry

A summary of seawater chemistry parameters is presented in Table 1. Temperature and pH were significantly different among treatments (F1:106 = 169.98, p < 0.001; F1:106 = 1412.22, p < 0.001, respectively), as expected for the different temperature and pH levels tested. No significant effects were observed for salinity (Supplementary Table S1 online on supplementary materials). As expected, time had a significant effect for temperature and salinity (F2:106 = 65.87, p < 0.001; F2:106 = 9.36, p < 0.001, respectively), reflecting seasonal variations in the fjord.

Table 1 Seawater chemistry parameters (mean ± s.d.) over the experimental period (June–August) in the 4 treatments: environment (EVM), temperature (T°), pH and temperature*pH (T°xpH). pH on the total scale (pHT), temperature (T; °C), salinity (PSU) and total alkalinity (TA; µmol kg−1) were used to calculate CO2 partial pressure (pCO2; µatm).

Reproductive success

Couples breed throughout the entire breeding season. A total of 134 clutches were produced by 29 breeding pairs: 36 laid and 20 hatched in the environment treatment; 30 laid and 3 hatched at high temperature; 38 laid and 24 hatched at low pH; 30 laid and 1 hatched in combined high temperature and low pH. Temperature and pH had no significant impact on the number of clutches laid per breeding pair (Supplementary Table S2 online on supplementary material). The number of clutches that successfully hatched was not influenced by pH (Supplementary Table S2 online on supplementary material). High temperature clutches had 3 to 10% of hatching success (Fig. 1), and this low number of clutches hatched at high temperature prevented a statistical comparison of the hatching success.

Figure 1
figure 1

(a) Effect of temperature (+ 0 °C, + 3 °C) and pH (8.0, 7.6) on the percentage (%), of hatched clutches. (b) Effect of pH (8.0, 7.6) on larva length (cm) at hatching (mean ± S.E.)

High temperature significantly decreased the average number of eggs per couple (30% less eggs in comparison to the environmental treatment), the egg area (15 to 20% smaller in comparison to the environmental treatment) and the overall reproductive output (number of eggs*mean egg area) (F1:104 = 14.24, p < 0.001; F1:104 = 73.92, p < 0.001; F1:104 = 41.53, p < 0.001, respectively Fig. 2), while low pH and the interaction of temperature and pH had no significant effect (Supplementary Table S2 online on supplementary material). The average number of eggs laid and the overall reproductive output also differed among parental pairs, highlighting intraspecific variability among the couples (Supplementary Table S2 online on supplementary material). The embryonic phase of successfully hatched clutches lasted for 5–6 days, in the environment and pH treatments (Supplementary Table S2 online on supplementary material). pH did not influence the size of the larvae at hatch (Fig. 1), but the length and duration of the embryonic phase differed among parental pairs (Supplementary Table S2 online on supplementary material). The low number of clutches that successfully hatched in the temperature and temperature × pH treatments prevented a statistical comparison of duration of embryonic phase and larval size at hatch under these conditions.

Figure 2
figure 2

Effect of temperature (+ 0 °C, + 3 °C) and pH (8.0, 7.6) on the number of eggs per clutch (a), egg area (b), and the reproductive output (c) mean ± S.E.

Paternal investment, assessed through the measure of time spent on fanning, cleaning, defending and chasing, was not significantly affected by temperature, pH or their interaction (Fig. 3, Supplementary Table S3 online on supplementary material). For all parameters of the behaviour, significant differences among parental pairs were registered (Supplementary Table S3 online on supplementary material), highlighting intraspecific variability among the couples. Males spent more 55% to 65% of the time inside than outside the nest (F1:169 = 90,484, p < 0.001), as would be expected during the embryonic phase.

Figure 3
figure 3

Time (in second; mean ± S.E.) spent by males fanning the eggs; cleaning the eggs; defending the nest; chasing females; and spent inside and outside the nest. All the measures were made from day 2 to day 4 after spawning in the 4 treatments.

No treatment difference was found in males and females across treatments for Fulton’s K (Supplementary Table S4 online on supplementary material). The gonadosomatic index was significantly affected by temperature in males (F1:24 = 32.07, p < 0.001), but no effect was found in females (Supplementary Table S4 online on supplementary material). No significant results were found in females’ hepatosomatic indexes (Supplementary Table S4 online on supplementary material), while pH had a significant effect in the males’ hepatosomatic index (F1:24 = 8.052, p = 0.009).

Discussion

Our results show that reproduction in two-spotted goby is strongly affected by elevated temperature. Couples in all treatments were able to breed, however hatching success was negatively affected by elevated temperature, with only 3–10% of the clutches successful hatched, when compared with the environmental group. High temperature also affected the average number of eggs per couple (less 30% than in environmental couples), and size of the eggs (eggs were 15–20% smaller compared to the environmental treatment).

Similar results were reported in a study where high temperature decreased the hatching success of Clupea harengus32. Clutches produced under higher temperature had fewer and smaller eggs, and these couples had overall a smaller reproductive output. On the contrary low pH and the interaction of high temperature low pH did not have any effect on the number and size of the eggs. This is consistent with the study on Amphiprion melanopus, that found temperature had a greater impact on reproductive performance than CO230. Another study on two-spotted gobies, sampled at the same field site, did not find any significant effect of acute exposure to decreased pH (pH 7.6) on spawning or the clutch size21. For this population, acclimation times, acute (few days in Forsgren et al.21) or chronic (1 month, this study) did not modulate the response to decreased pH. Different results were found for the same species by Faria et al.3. Reproduction was stimulated under low pH with an increase in number of clutches and number of eggs. However, this came at a cost, with smaller eggs and larvae at hatch3. This is likely a consequence of different experimental conditions, as Faria et al.3 exposed fish to more extreme pCO2 levels than the levels in the current study (~ 2300 vs. ~ 1300 µatm, respectively). However, it could also be attributed to a population effect, as Faria et al.3 tested fish from the Portuguese coast (the southern limit distribution of two-spotted goby), while in the current study fish were sampled from a population that inhabits the Gullmarn fjord, in the Northern Sea. Different populations of a given species can present different sensitivities to pH through local adaptation; populations inhabiting more variable environments being able to tolerate lower pH33. The Gullmarn fjord has variable pH with variations up to 0.9 pH units within a single month and pH as low as 7.6 being documented34. Our tested scenario (pH 7.6) being within the natural range of variability, it is then not surprising that only mild effects were observed.

We observed a shorter embryonic duration (5–6 days, under all treatments) as compared to the average 10 days of embryonic development reported in other studies3,21. This may be a consequence of a particularly warm temperature in Sweden in 2018. The two-spotted goby experiences an average temperature of 16 °C during their reproductive season21,35, and in 2018 the maximum recorded temperature during the breeding period reached, 20 °C (Kristineberg Weather Station—https://www.weather.loven.gu.se/kristineberg/en/data.shtml). These record temperatures might have contributed to the observed lower reproductive success under the high temperature treatments.

The two-spotted goby, as most goby species, provide paternal care during the embryonic phase, and, as expected, the time the males spent inside the nest, taking care of the eggs, in all treatments was significantly higher than time spent outside6,36. Our results suggest that neither high temperature, low pH or the interaction of both, influences the time males spend on the different parental care behaviours. All these behaviours were recorded from day 2 to day 4 of the embryonic phase. Although males spent time and energy providing care for the eggs, in the end this investment did not pay off in some of the treatments, as the eggs most likely died and were removed by the male. The increase of temperature can increase fungus, infections and general organic matter decomposition can affect the eggs36, and would explain why clutches under elevated temperature did not successfully hatch6.

Body condition of the two-spotted goby was similar in all treatments, as suggested by the lack of significant differences in Fulton’s K index. A previous study on the effects of low pH on reproductive performance of the two-spotted goby also showed no difference in the K index3. Although fish in these treatments can have higher metabolic rates as an exposure to higher temperature, the ad libitum feeding ensured they can compensate for the extra energetically costs and maintain their physical condition under elevated temperature and low pH29. The ad libitum food may also have prevented males from eating clutch’s eggs37, allowing the focus to be solely on the tested stressors.

GSI was not affected by the treatments for females but was negatively affected by high temperature for males. The gonadosomatic index is used to assess the development of the gonads. Higher values within a reproductive season are usually attributed to the reproduction peak. While low values are attributed to the end of the breeding season38. The two-spotted goby have a short life-cycle and reproductive season. Studies on the Gullmarn fjord population mating behaviour showed a transition from strong male–male competition at the beginning of the season, to strong female–female competition later in the season39. Our experiment ended in August, close to the end of the breeding season4. This can explain the different values of GSI between sexes. Moreover, the significantly lower investment in male gonads at elevated temperature explain the poor quality of the clutches, or why clutches failed to hatch.

pH had a significant effect on males HSI while temperature had no effect. Acanthochromis polyacanthus subjected to an increase of 3 °C had larger livers40. Larger livers are likely related to higher needs for energy storage under metabolic stress at higher temperature41. However, in the two-spotted goby, neither high temperature nor the interaction of high temperature and low pH had an effect on the HSI. Low pH had a significant effect, that might be related to the stress caused by pH. There is no evidence suggesting that liver size directly compromises reproduction, nonetheless indirectly it might influence the size of the gonads (due to space limitation in the abdominal cavity) and overall the reproductive and metabolic processes.

Metabolic demands increase with temperature32. This could explain the lower investment in reproductive output at high temperature observed in this study. This is consistent with previous studies showing that temperature influences reproduction in fish14,42. Exposure to high temperature within the thermal range fishes are used to endure will make them spend more energy to maintain cellular functions and impair other behaviours not directly essential for their survival30. The effect of temperature changes on reproduction is well documented in fish from warmer regions43. Studies on species exposed to wider thermal ranges are more scarce. Reproduction performance is an important part of fish development and it is then crucial to understand the impacts of multiple stressors such as pH and temperature in cross-factorial studies28,44. Out study shows that exposure to extreme present conditions of temperature and pH already have negative effects on multiple aspects of reproduction. Projections show that such conditions will be more frequent and be more extreme in the near-future. For example, temperature up to 21 °C have been documented in June 2020 in the Gullmarn fjord (Kristineberg Weather Station—https://www.weather.loven.gu.se/kristineberg/en/data.shtml). Without drastic CO2 mitigation strategies, local population of two-spotted gobies will then be exposed to more challenging conditions eventually compromising their survival.

Reproduction involves a series of hormonal, physiological and behavioural processes. Future work should consider not only reproductive performance, but also maternal and paternal investment, offspring development, and acclimatization factors in transgenerational experiments. Future works should look into local adaptation in different populations, using the same experimental designs and conditions to test different populations, so results can be compared, as well as other modulating factors, such as food availability. Additionally, transgenerational experiments, where more than one generation is bread in the lab under different environmental conditions, would be important to evaluate adaptation and acclimatization potential to future conditions. Overall, these studies are necessary to better comprehend how reproduction will be affected by global changes, and how will different species, and populations cope with those changes.

Methods

Fish collection

The two-spotted gobies were collected in Sweden in the Gullmarsfjord (58° 6 N, 11° 10E), from 25th May to 12th June 2018, by snorkelers using hand-dip nets. In the Gullmarsfjord, the two-spotted goby inhabits seaweed beds on the fjord rocky shores4. During breeding season, they seek protected and moderately exposed algal areas in shallow waters21.

Aquarium system

In the laboratory, single sex groups of 8 fish were housed in eight 55L aquaria maintained at 17 °C and pHT 8.00, for a 3-day quarantine period. Flow-through aquaria were used with seawater pumped directly from the fjord at 5 m of depth. After laboratory acclimation, fish were exposed to a total of four treatments: 1) Environment (EVM)- + 0.0 °C × pHT 8.0; 2) pH- + 0.0 °C × pHT 7.6; 3) Temperature (T°)- + 3.0 °C × pHT 8.0; 4) Temperature × pH (T°xpH)- + 3.0 °C × pHT 7.6. pH 7.6 was selected as the average pH in the fjord by 2100 and the extreme of the present natural variability34. Eight replicates per treatment were used. Temperature and pH were slowly changed in elevated temperature and low pH aquaria (+ 1 °C/per week; − 0.1 pH/per week), until temperature was + 3 °C as compared to ambient conditions and pH was 7.6. pH was maintained through bubbling of pure CO2 and using pH-computers (Aqua Medic, Bissendorf Germany) in each aquarium. In the other aquaria, temperature and pH were maintained at the natural variability of the fjord, throughout the duration of the experiment. High temperature treatments were adjusted throughout the experiment to always maintain a + 3 °C difference from the in situ temperature. Fish were kept under a light cycle of 16/8, and were fed Artemia twice a day, ad libitum.

After 1 month of acclimation to the new conditions, all 64 fish were measured and weighed. A male and female were chosen at random from the same treatment, placed in one 30L aquaria, with temperature and pH conditions matching the ones from the acclimation aquaria, and allowed to breed.

Carbonate chemistry

Temperature and pH (on the total scale—pHT) were measured every 2 days. pH probe was calibrated every week using TRIS (Tris/HCL) and AMP (2-aminopyridine/HCl) buffers (salinity 25) following best practices45. Total alkalinity (TA) was measured weekly in each aquarium from filtered water samples by titration. Salinity data was obtained from the continuous monitoring of the flow-through system at the Kristineberg station. Other parameters of the carbonate systems parameters were calculated from temperature, salinity, pHT and TA using CO2SYS46.

Endpoints

Breeding couples were provided with a pvc tube (10 cm, Ø3.2 cm) for shelter and nesting. Each pipe was line with a removable acetate sheet, where the spawning females could attach their eggs. Presence of an egg clutch was checked daily. In the presence of eggs, the acetate sheet was carefully removed, placed in a Petri dish filled with seawater from the same aquarium and then photographed for later counting. Approximately 10 eggs were gently removed from the clutch, placed in a smaller Petri dish and photographed under a stereomicroscope for measurement of their surface area. Couples were not affected by the brief removal of the acetate sheet, as the fish resumed normal activity after the acetate sheet was put back in the nest. Males with clutches were recorded with a Panasonic HDC-SD90 camera, for 10 min in the morning 1 h after being fed, as it has been shown that food availability has no effect on parental care for the two spotted goby6. Studies showed records of parental care averaging 8 to 20 min are enough to observe fish parental care36,47, which corroborates our personal observations. All recordings were done from day 2 to 4 after spawning. Videos were subsequently analysed using the Observer XT software (Noldus technology version 5, Netherlands) to quantify the time spent: (1) fanning, (2) cleaning, (3) defending, and (4) chasing the female, as well as (5) total time spent outside and inside the nest. Fanning behaviour consisted of swimming movements with the caudal, dorsal and pelvic fins, to aerate the eggs (adapted from Baklow48). Cleaning behaviours included the removal of dead, unfertilized or infected eggs, and debris adapted from Blumer49). Defending was defined as the guarding position the male assumed on the nest entry, to avoid intruders. Chasing the female behaviours consisted in charge and chasing the female, outside the nest, to prevent her entrance in the nest.

Reproductive output, as a proxy for maternal investment, was calculated as the overall area (mm2) of the clutch by multiplying the total number of eggs per clutch by the average of individual egg area in the same clutch. Close to hatching, the acetate sheet was transferred to a maternity cage (16.5 × 12.5 × 11.5 cm) inside a 55L aquaria matching the temperature and pH conditions of the parental treatment. Ten to 15 larvae were collected at hatching, euthanized using an overdose of MS222, and immediately photographed under a stereomicroscope to measure standard length.

At the end of the experiment, adults were euthanized with an overdose of MS222, weighed and measured, to determine the Fulton’s K condition, K = 100*(W/SL3), where W is wet mass in grams and SL is standard length in cm. Fish were immediately dissected and gonads and liver were weighed to determine the gonadosomatic (GSI) and hepatosomatic (HSI) index using the formula: GSI/HSI = (gonad or liver mass (g)/fish mass (g)) * 100.

Statistical analysis

Data normality was assessed using a Shapiro–Wilk test and homogeneity of variances using a Barlett test. Effect of the different treatments (temperature, pH and time) on the physico-chemical parameters (temperature, pH and salinity) were assessed using a 3 factor linear model Anova (temperature, pH and month). Impact of temperature and pH on the number of clutches laid and hatched, and body condition of the adult fishes (Fulton’s K, gonadosomatic and hepatosomatic indexes) were tested using a 2 factor linear model Anova (temperature and pH). For the number of eggs per clutch, area of the eggs, reproductive output and male reproductive investment, a 3 factor linear model Anova (temperature, pH, and parental pair nested within the treatment) was used. Aquariums with one couple (parental pair) each were replicates within their treatment. Three couples did not reproduce during the experiment, and were excluded from the statistical analysis.

All the statistical analysis was performed with R statistics, RStudio (Version 1.1.463), and results were considered statistically significant at p < 0.01.

Ethics statement

The research was performed in strict accordance with the current Swedish national legislation on Animal Welfare, the experimental procedures were approved by the Swedish Board of Agriculture and the ethical permission was given by the Gothenburg Animal experiment Ethical board under the licence Dnr 103-2014.