# Impacts of Long-Term Exposure to Ocean Acidification and Warming on Three-Spined Stickleback (Gasterosteus aculeatus) Growth and Reproduction

Jimmy Devergne, Véronique Loizeau, Christophe Lebigre, Anne Bado-Nilles, Sophie Collet, Olivier Mouchel, Ugo Iaria, Marie-Madeleine Le Gall, Lauriane Madec, Cyril Turiès, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Jimmy Devergne, Véronique Loizeau, Christophe Lebigre, Anne Bado-Nilles, Sophie Collet, et al. Impacts of Long-Term Exposure to Ocean Acidification and Warming on Three-Spined Stickleback (Gasterosteus aculeatus) Growth and Reproduction. Fishes, 2023, 8 (10), pp.523. 10.3390/fishes8100523. hal-04315890

HAL Id: hal-04315890

https://hal.univ-brest.fr/hal-04315890

Submitted on 1 Dec 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

# Impacts of Long-Term Exposure to Ocean Acidification and Warming on Three-Spined Stickleback (*Gasterosteus aculeatus*) Growth and Reproduction

Jimmy Devergne<fnref n="1" />,*, Véronique Loizeau<fnref n="1" />, Christophe Lebigre<fnref n="2" />, Anne Bado-Nilles<fnref n="3" />, Sophie Collet<fnref n="1" />, Olivier Mouchel<fnref n="1" />, Ugo Iaria<fnref n="3" />, Marie-Madeleine Le Gall<fnref n="1" />, Lauriane Madec<fnref n="1" />, Cyril Turiès<fnref n="3" /> and Arianna Servili<fnref n="1" />

<figure><img src="image_4.png" /></figure>

check for
updates

**Citation:** Devergne, J.; Loizeau, V.; Lebigre, C.; Bado-Nilles, A.; Collet, S.; Mouchel, O.; Iaria, U.; Le Gall, M.-M.; Madec, L.; Turiès, C.; et al. Impacts of Long-Term Exposure to Ocean Acidification and Warming on Three-Spined Stickleback (*Gasterosteus aculeatus*) Growth and Reproduction. *Fishes* **2023**, *8*, 523.
https://doi.org/10.3390/fishes8100523

Academic Editor: Li Li

Received: 30 August 2023

Revised: 6 October 2023

Accepted: 15 October 2023

Published: 21 October 2023

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** The warming and acidification of surface waters as predicted by the IPCC leads aquatic species to face major multifaceted changes in their environment. Although teleosts have efficient regulatory systems to cope with these changes, such changes clearly have the potential to impact their physiological functions. Hence, it is crucial to estimate the ability of teleost fishes to cope with multi-stresses to predict how they will deal with future environments. In this context, we investigated the joint effect of warming and acidification on three-spined stickleback (*Gasterosteus aculeatus*) from the juvenile stage to adulthood, focusing on parameters linked to growth, sexual maturation, and reproduction. Juvenile sticklebacks were split in 2 climate scenarios: a “Current” scenario corresponding to the current seasonal physico-chemical parameters of the water of the “Rade de Brest” in France, and a “RCP8.5” scenario with a warming of 3 °C and an acidification of 0.4 pH units. After 7 months, fish in the RCP8.5 scenario reached the same size and mass as those in the Current scenario, but they needed greater amounts of food to reach satiety. Furthermore, the mortality rate over the experiment was higher in the RCP8.5 scenario. Muscle lipid content, an indicator of energy reserves, was lower in females in the RCP8.5 scenario, suggesting an increased need for energy to maintain homeostasis and other physiological functions or a divergence in energy allocation strategy. Moreover, females exhibited lower sexual maturation and egg quality under the RCP8.5 scenario, which could have contributed to the lower fertilisation rate observed. Males were more resilient to the RCP8.5 scenario, exhibiting only a trend for lower kidney somatic index scores. Altogether, these results suggest a delay and/or an inhibition of gametogenesis and maturation in fish in warmed and acidified waters. The analysis of blood sex steroid concentrations, brain gene expression profiles, and physiological indexes did not allow us to discriminate between a delay and an inhibition of maturation in the RCP8.5 scenario. Overall, these findings clearly indicate that there is a long-term global impact of combined acidification and warming on the mortality and reproductive performance of three-spined stickleback.

**Keywords:** *Gasterosteus aculeatus*; climate change; life cycle; sexual maturation; long-term acclimation; lipids; multi-stress

**Key Contribution:** (1) In a warming and acidification experiment (+3 °C; -0.4 pH), conducted over 7 months, with an *ad libitum* feeding, the growth of stickleback (*G. aculeatus*) was not affected but the mortality rate was higher. (2) The lower fertilisation rate was associated with lower muscle lipid content in females and poorer egg quality.

<footnote n="1">¹ University of Brest, CNRS, IRD, Ifremer, UMR 6539, LEMAR, 29280 Plouzané, France</footnote>

<footnote n="2">² UMR DECOD (Ecosystem Dynamics and Sustainability), IFREMER, INRAE, Institut Agro, 29280 Plouzané, France</footnote>

<footnote n="3">³ Unité d'Ecotoxicologie In Vitro et In Vivo (ECOT), Institut National de l'Environnement Industriel et des Risques (INERIS), 60550 Verneuil-en-Halatte, France</footnote>

<footnote>* Correspondence: jimmy.devergne@ifremer.fr</footnote>
# 1. Introduction

Since the beginning of the industrial revolution in the 19th century that generated the expansion of anthropic activities, the levels of greenhouse gas emissions have never stopped increasing, with for instance, the rise in the partial pressure of carbon dioxide (pCO₂) in the atmosphere from 280 µatm during the preindustrial period to 380 µatm nowadays [1]. The increased atmospheric CO₂ has led to alterations in multiple physicochemical parameters of seawater (temperature, pH, oxygen and salinity [1,2]). The most pessimistic scenario of the Intergovernmental Panel on Climate Change (IPCC), named representative concentration pathway 8.5 (RCP8.5), predicts a warming of the sea surface temperature of 3.04 ± 0.62 °C and an acidification of 0.38 ± 0.005 pH units by 2100 [1,2].

Such changes in environmental conditions can affect the physiology and behavior of aquatic organisms as most of them are poikilothermic (i.e., unable to regulate their internal temperature), including a majority of fish species [3–5]. As a result, all of their physiological functions can be directly influenced by environmental temperature fluctuations [6]. The effects of warming on fish physiology have been widely studied, showing, for example, that even a rise of few degrees in the environment can lead to an increase in growth and metabolic rates in several fish species [7–12]. Warming is known to also impact fish sexual maturation [4,13,14], steroidogenesis [14–17] and development [3,18] by modifying molecular structures and rates of biochemical reactions with an end impact on enzymatic activities [19,20], gene expression profiles [14,16,17,21] and plasma hormone concentrations [4,15,22], as reviewed by many studies [23–25]. Ocean acidification is another major process directly related to increased pCO₂. This environmental hypercapnia leads to acidosis of the body/fluids of aquatic fauna by accumulation of carbonic acid (H₂CO₃) [26]. In teleosts, who possess an efficient acid-base regulation, body acidosis is compensated by increased gill bicarbonate concentration (HCO₃⁻), which could potentially drive an alteration in internal ionic concentrations, with possible consequences on neurotransmission [27,28]. The end result could be impacts on key physiological functions (basal metabolism, growth, reproduction) [29–32], behaviour (swimming activity, foraging, etc.) [31,33], and ultimately increase the vulnerability of the species [30,31].

As ocean warming and acidification will occur simultaneously under future climate change circumstances, it is essential to consider the combination of these factors and determine how organisms will cope with future environments. Results from previous studies evaluating the influence of several stressors combined are increasing. Most of these studies reported the effects of ocean warming and acidification on the early stages of physiological functions of the commercially important species. For instance, several studies about the futures projections of warming (4–5 °C) and acidification (0.3–0.4 pH unit) showed negative impacts on growth [34,35], survival [34], malformations [34,35], and digestive enzymes activity [35] in *Solea senegalensis*. In similar experimental designs, decreases in larval survival of Atlantic cod [36], gilthead seabream, and meagre [31] were also reported, associated in the two latter cases with increased malformation rates and behavioural impairment [31]. Regarding the combined effect of warming and acidification on the metabolism of early stages of marine fish, the reported results exhibit different physiological responses depending on the species [35,37,38], indicating inter-specific variations in sensitivity. Altogether, this reveals the need of further studies on different species (i.e., including non-commercially important species) and longer exposure duration also covering the most advanced life stages to enable more realistic predictions of the impact of climate change on fish physiology. Indeed, (i) over focus on a few model species makes it difficult to extrapolate results to other context and (ii) limiting the evaluation to short-term effects at early stages can be misleading since the main consequences of warming and acidification may occur in later life stages.

The three-spined stickleback (*Gasterosteus aculeatus*, Linnaeus, 1758) is an amphihaline teleost fish which is widely distributed over the northern part of the globe. It has been heavily used as a model species, providing a strong background knowledge of its physiology and ecology [39–41]. Furthermore, the full sequencing of its genome provides
tools to investigate the effects of environmental variations in its physiology at the molec-
ular level [42,43]. Its small size and short life cycle make this species easy to rear in lab
facilities, enabling the possibility to follow all life stages in a relatively short time span
(sexual maturity at one year). These characteristics make the stickleback a good model for
long-term studies of environmental variations in fish physiology.

Several studies have already shown direct effects of increased temperatures in stick-
leback, such as lower growth and mass [44], lower fecundity [10], smaller egg diameter
and faster egg development [45–47]. Metabolic and redox processes were also shown to
be upregulated at higher temperatures, as was the expression of certain genes involved
in reproduction and growth, such as smtlb, pomca and tshba [10,48]. Previous studies have
also reported the impacts of acidified water, with higher fecundity under high acidifica-
tion [32] and no effect of pH on egg survival and diameter [45]. However, to our knowledge,
the resilience of this species to multiple stressors of climatic and chemical origin is still
unknown.

In this context, we aimed at testing whether ocean acidification and warming together
could negatively impact the growth and reproductive performance of three-spined stick-
lebacks. For this purpose, juvenile sticklebacks were subjected to multi-stress (pH and
temperature), as predicted by the IPCC RCP8.5 to occur by 2100, for 7 months, until adult-
hood. According to the new set of scenarios, the shared socio-economic pathways (SSP),
which combine pathways of atmospheric greenhouse gas concentration (as in the previous
RCP scenarios) to socio-economic projections [49], the experimental scenario tested in the
present study, named RCP8.5, would correspond to both the SSP3-7.0 and the SSP5-8.5
scenarios [50].

During this multi-stress experiment, and in addition to estimating mortality rates, we evaluated the impact of acidification and warming on growth by estimating size, mass and physiological indexes as well as the feed conversion ratio in each experimental scenario. The potential disruption of the individuals' sexual maturation was tested by combining hormonal, histological, and molecular approaches (at the final sampling point, after 7 months of the experiment). Finally, we investigated the impact of ocean warming and acidification on the individuals' reproductive success by comparing eggs and sperm quality as well as the fertilisation rates. We are aware that without a complete cross-factorial design, it is not possible to distinguish whether the observed effects on fish physiology are due to just one factor or to their interaction. However, in the present study, the goal is to evaluate the global impact of ocean acidification and warming together to obtain a more realistic understanding of what will occur in nature.

**2. Materials and Methods**

2.1. *Animals and Experimental Design*

A population of freshwater juvenile three-spined stickleback (*Gasterosteus aculeatus*), born in July 2020, were obtained on the 7th of October 2020 from the experimental freshwater mesocosms of the *Institut national de l'environnement industriel et des risques* (E60-769-02, INERIS, Verneuil-en-Halatte, France). The stabling of the fish was performed within Ifremer's experimental facilities in Brest (Brittany, France), at the aquaculture facility of the laboratory of Ecophysiology and Life History Traits of Marine Organisms (agreement number B29-212-05). On arrival at the Brest experimental facility, these fish were placed for a week in water at a salinity level of 10‰, then gradually acclimated to seawater (increasing the salinity by 10‰ per week) until the conditions reached the 34‰ salinity level of Brest Bay. The fish were fed once a day with frozen Artemia and bloodworms (proportion 1:4, Antinea, France). The food was distributed progressively into the tank until satiety was reached. The satiety of the fish population in a tank was determined by observation as a cessation of eating behaviour. The theoretical food consumption was calculated as the difference between mass of food available and mass of food not distributed, once satiety was reached. This amount was used to adjust the ration per tank for the next day. One month after salinity acclimation, the juvenile population, then around 4 months old, was
randomly divided into 2 tanks corresponding to 2 experimental groups of 557 fish per tank. Then, on 12 December 2020, after acclimating to these new conditions for around 3 weeks, during which no mortality was observed, the 5-month-old fish were subjected to two experimental scenarios. The rearing conditions of the 1st group (labelled "Current scenario") correspond to the current temperature and pH of the Bay of Brest. The second group (labelled "RCP8.5 scenario") was subjected to warming (+3 °C) and acidification (-0.4 pH units). These scenarios simulate the IPCC RCP8.5 scenario predicted for 2100 [26]. The experimental tanks were supplied with flow-through seawater directly pumped at a location 500 m from the coastline from a depth of 20 m in the Bay of Brest. The water was filtered by passing through a sand filter (Waterco Lacron, Australia), heated via a tungsten plate heat exchanger (Vicarb, Sweden), degassed through a column packed with plastic rings, refiltered with a 2 µm membrane, and sterilized by a UV lamp (PZ50, 75W, Ocene, France). The RCP8.5 scenario was supplied with sea water via a header tank where gas CO₂ was diffused in. The water temperature and pH were monitored daily (Figure 1) with a WTW 3110 pH meter (Xylem Analytics Germany, Weilheim, Germany; with electrode: WTW Sentix 41, NBS scale). The devices were daily calibrated with NBS certified WTW technical buffers pH 4.01 and pH 7.00 (Xylem Analytics Germany, Weilheim, Germany; the values are summarized in Table 1). Juveniles and adults experienced a temperature range from 11.3 to 17.9 °C in the Current scenario and from 13.2 to 20.6 °C in the RCP8.5 scenario, following the natural temperature fluctuation of the Bay of Brest. The photoperiod regime was also set to simulate the natural photoperiod of the Bay of Brest (Halogen lamp at 42W, 55–60 lux).

<figure><img src="image_6.png" /><figcaption>**Figure 1.** Temperature (°C) in the Current (in grey) and RCP8.5 (in black) scenarios throughout the experiment.</figcaption></figure>

<table><caption>**Table 1.** Mean physico-chemical conditions recorded twice a month during the total rearing period (mean ± S.D.). TA: total alkalinity, pCO₂: partial pressure of CO₂, mean ± se calculated on the monthly (7 months) measure of alkalinity.</caption><thead><tr><th rowspan="2">Rearing Scenarios</th><th rowspan="2">Salinity (‰)</th><th rowspan="2">O₂ (%)</th><th colspan="2">Temperature (°C)</th><th rowspan="2">pH Total</th><th rowspan="2">TA (µmol·kg<sub>sw</sub><sup>-1</sup>)</th><th rowspan="2">pCO₂ (µatm)</th></tr><tr><th>Min</th><th>Max</th></tr></thead><tbody><tr><td>Current</td><td>32.70 ± 0.98</td><td>96.54 ± 4.49</td><td>11.90</td><td>17.90</td><td>7.95 ± 0.06</td><td>2184.07 ± 86.59</td><td>527.78 ± 135.91</td></tr><tr><td>RCP8.5</td><td>32.79 ± 0.99</td><td>96.67 ± 3.82</td><td>15.40</td><td>20.30</td><td>7.53 ± 0.05</td><td>2193.56 ± 138.72</td><td>1509.03 ± 265.52</td></tr></tbody></table>
The total alkalinity of the seawater was measured twice a month according to the protocol described in Cohen-Rengifo et al. [51] adapted from the protocol of Anderson and Robinson [52] and Strickland and Parsons [53]. Salinity and Oxygen (O₂) were measured simultaneously with a WTW LF 340 salinometer (Xylem Analytics Germany, Weilheim, Germany; with electrode: WTW Tetracon® 325/C, ‰ scale) and a WTW Oxi 340 oximeter (Xylem Analytics Germany, Weilheim, Germany; with electrode: WTW CellOx 325, ‰ scale). Then, the pCO₂ was calculated through the CO₂sys macro [54] with the constants from Mehrbach et al. [55] and the refit by Dickson and Millero [56]. The total scale pH was afterwards calculated using pCO₂ values from the same macro (Table 1).

## 2.2. Non-Invasive Biometric Approaches

During the 7-month rearing period, a sample of 144 fish were collected (and returned) from each experimental scenario every two months to assess their mass and size and subsequently calculate the thermal growth coefficient of each experimental group. Se- dated (0.05 g·L⁻¹ of Tricaine methane sulfonate 222 (MS222), Pharmaq, Fordingbridge, Hampshire, UK) fish were weighed (Mettler-Toledo, MS4002TSDR, Switzerland) and pho- tographed (NIKON, D7200, Japan) in 12 batches of 12 fish. Because MS222 residues are eliminated rapidly from the tissues of fish medicated by bath administration, even when used repeatedly [57], and due to the low doses used here, we exclude any long-term impact of MS222 on fish physiology. To measure the individuals' size, pictures of each batch were taken in front of a scale and measurements were carried out using the ImageJ software (Version 1.53e, Java 1.8.0_172; [58]). The feed conversion ratio (FCR; [59]), mortal- ity rate, and thermal growth coefficients (TGC; [60]) were measured using the following formulas: FCR = (Total mass of food consumed in the tank)/(Biomass gain in the tank); Mortality rate = (Total number of dead fish)/(Total number of fish); TGC = (Mass<sub>final</sub><sup>1/3</sup> – Mass<sub>initial</sub><sup>1/3</sup>)/(Σ(Temperature)), where Σ(Temperature) corresponds to the accumulation of temperature in degrees Celsius experienced by the stickleback.

## 2.3. Final End-Points

### 2.3.1. Individual Physiological Indexes

In July, after 7 months of acclimation to the climate conditions, a sampling was carried out on 33 fish in total. The remaining adult sticklebacks were preserved for longer exposure for use in further studies. The fish were fasted for 24 h to reduce faecal contamination of measures and risk of regurgitation. Before sampling, the fish were anaesthetised with 0.125 g·L⁻¹ MS222 in the tank, and then fished and euthanized with a lethal concentration of 0.250 g·L⁻¹ MS222. Individual biometric data, including the standard body length, body mass, Fulton index corrected from the mass of the gonads (K<sub>corr</sub>), gonado-somatic index (GSI), hepato-somatic index (HSI) and kidney-somatic index (KSI), were estimated (Sartorius, MC BA 1000, Germany) using the following formulas [61]: K<sub>corr</sub> = ((Mass of the fish – Mass of the gonads)/(Size of the fish)<sup>3</sup>) × 100; GSI = ((Mass of the gonads)/(Mass of the whole body)) × 100; HSI = ((Mass of the liver)/(Mass of the whole body)) × 100; and KSI = ((Mass of the kidney)/(Mass of the whole body)) × 100.

### 2.3.2. Lipids Extraction and Analysis

The total lipid content of fish muscle was assessed individually. During the extraction, the atmosphere in the tube was replaced by azote gas as far as possible to preserve the lipids. White muscle tissue was flash frozen in liquid nitrogen after sampling and kept at -80 °C until extraction. The total lipid content of each muscle sample was extracted following a modified Folch method [62]. The micro-extraction was initiated by homogenising 100 to 500 mg of tissues, with the addition of an equivalent amount of Na₂So₄, using a polytron in Dichloromethane:Methanol (98:2) solution [62]. Afterword, the following extraction process was repeated thrice: after shaking for 3 min with a vortex, the homogenate was run through an ultrasonic tank for 5 min and centrifuged at 4 °C for 10 min at 1.500 G to collect the supernatant in a separate tube. The samples were concentrated in a SpeedVac
at 45 °C for 1 h. The concentrates were then transferred to clean, tared, round-bottomed reaction tubes and all the remaining solvent was removed by impinging with a stream of nitrogen. The total lipid extracts (TL) were measured to the nearest 0.01 mg using an analytical balance (Sartorius, MSE225P-100-DU, Germany).

### 2.3.3. RNA Extraction and Purification

The fish brains were quickly removed and the median region containing the hypothalamus was dissected and stored in an RNA stabilization reagent (RNAlater, Qiagen, Hilden, Germany) at 4 °C for 24 h and then placed at -20 °C until the extraction. The total RNA extraction and purification of the median brain was performed with the Nucleospin® RNA<sub>XS</sub> kit (Macherey Nagel, Germany) according the protocol of the supplier (Version January 2020/Rev.10). The concentration and purity of the RNA were assessed with Nanodrop<sup>TM</sup> 2000 (Thermo Scientific Inc., Waltham, MA, USA) and the integrity of the RNA was assessed using Agilent Tapstation 4150 (Agilent Technologies Inc., Santa Clara, CA, USA). All the samples showed an integrity number (RIN) higher than 7.4 and were stored at -80 °C until used in reverse-transcription quantitative PCR (RT-qPCR) analysis.

### 2.3.4. RT-qPCR Analysis

The reverse transcription (RT) of cDNA was performed on all the RNA samples with positive and negative (without enzyme) reaction by following the protocol recommended by the supplier, with 500 ng of RNA, using the iScript™ cDNA Synthesis kit (Bio-Rad Laboratories Inc., Hercules, CA, USA) and run in ThermalCycler T100TM (Bio-Rad Laboratories Inc.). Then, cDNA was stored at -20 °C until qPCR.

The qPCR was performed on the following genes of interest: Gonadotropin Releasing Hormone 2 and 3 (gnrh2, gnrh3), aromatase (cyp19a1b), Kisspeptin 2 (kiss2), Gonadotropin-Inhibitory Hormone (gnih), and 2 housekeeping genes (Ribosomal Protein L8 (rpl8), and Ribosomal Protein L13A (rpl13a)). The primers were specially designed (Table 2) for this study with the primer 3Plus software [63].

<table><caption>**Table 2.** Designed primers used in the qPCR analysis in the median brain.</caption><thead><tr><td>Gene Symbol</td><td>Forward Primer</td><td>Reverse Primer</td><td>GenBank Accession Number</td></tr></thead><tbody><tr><td>RPL8</td><td>5'GTGAAGGACATCATCCACGA3'</td><td>5'CTCCTCCACACAGCAGATGA3'</td><td>ENSGACT00000002668.1</td></tr><tr><td>RPL13A</td><td>5'CGCCCTACGACAAGAGGAAG3'</td><td>5'CTCTGCCACCTTGGTCAACT3'</td><td>ENSGACT00000012382.1</td></tr><tr><td>GNRH2</td><td>5'TGTGTTGGAGCTCAGCTGTC3'</td><td>5'CTCTGGCTAAGGCATCCAA3'</td><td>ENSGACT00000011943.1</td></tr><tr><td>GNRH3</td><td>5'GTTGGTGGTCCAGGTCUCT3'</td><td>5'CTCTCTTGGGTCTGGGCACT3'</td><td>ENSGACT00000012668.1</td></tr><tr><td>CYP19A1B</td><td>5'ATACCCGGTCCGTGCTTCT3'</td><td>5'CCCGAATCTGGCTGTGTAGT3'</td><td>ENSGACT00000007929.1</td></tr><tr><td>KISS2</td><td>5'GGGTCAGTTCTTTCTACGCTCA3'</td><td>5'AATGTAGCGTTTCCCAAAGC3'</td><td>KT202354.1</td></tr><tr><td>GNIH</td><td>5'AGAGCTTCCGCATACTCTCG3'</td><td>5'ATGTTGTGGGTCGAACTGGT3'</td><td>KT202315.1</td></tr></tbody></table>

The qPCR was performed with the CFX96 Touch Real-Time PCR Detection system (Bio-Rad Laboratories Inc.). The efficiency of the primer pairs was calculated by using serial dilution (1/5, 1/10, 1/20 and 1/40) of a cDNA pool. Each well was filled with 5 µL of cDNA (1/20 dilution) to conduct the qPCR in duplicate, followed by 10 µL of reaction mix made up of 333 nM of each primer pair, 7.5 µL of SsoAdvanced Universal SYBR® Green Supermix (Bio-Rad Laboratories Inc.) and completed with of RNase/DNase-free water. The qPCR program initial activation step started at 98 °C for 2 min, followed by 40 cycles of 5 s at 95 °C and 20 s at 60 °C. The protocol ended with the generation of a high-resolution melting curve after the amplification phase in order to confirm the amplification of a single product in each reaction. The Gene Expression Module of the CFX Manager program (Bio-Rad Laboratories Inc.) was used to automatically determine the related quantification cycle (Cq) value for each sample. *rpl8* and *rpl13a* were used as reference genes in the ΔΔCt method to normalize the relative quantity of transcripts. No significant differences in Cq values were observed for housekeeping genes between experimental scenarios.
### 2.3.5. Histological Assessment of the Gonadal Maturation Stage

After sampling, and weighing, gonads were fixed in paraformaldehyde (PFA) 4% for 24 h at 4 °C. Then, the tissues were washed in phosphate-buffer saline (PBS) solution and 0.85% sodium chloride (NaCl), each for 15 min. Next, 5 washes of increasing ethanol solution (70% to 100%) were performed at 4 °C to dehydrate the samples. The samples were lightly shaken overnight in 100% ethanol, then in 2 xylene baths. Finally, they were impregnated into warm filtrate paraffin. The gonads were then embedded into cassettes, the anterior part closest to the cut. A microtome (RM 2035, Leica, UK) was used to make 5 µm thick slices starting at 200 to 250 µm from the apical lobe of the gonad. The slices were dried on glass slides at 37 °C overnight and warmed at 53 °C for 20 min before the trichrome Prenant-Gabe staining. The slides were dewaxed into xylene and gradually rehydrated in decreasing ethanol solutions (100% to 70%) and distilled water. Staining of the nucleus, in purplish blue, was obtained with 1% Weigert's Haematoxylin (1 min), cleared by tap and distilled water. The cytoplasm was stained in red with 1% aqueous 1% Eosin Y (8 min) and the conjunctive tissue in green after an aqueous Phosphomolybdique acid (30 s) bath by 0.05% aqueous Light Green (30 s). The tissues were dehydrated by 2 baths of 100% ethanol and xylene. The slides were mounted with Eukitt® medium (ORSAtec, Royaume-Unis). Analysis and photography (EOS 700D, Canon, Japan) were performed under a light microscope (AXIO observer Z1, ZEISS, Germany). In total, 5 different regions of 8 gonads per sex and scenario were analysed to estimate the maturity of gametes identified. To evaluate maturation, the stages described by Sokolowska and Kulczykowska [64] and Furin et al. [65] were pooled into 5 and 4 groups for females and males, respectively. This resulted in the following classification of maturation stages in ovaries: "Early" (stages 1 and 2), "Midterm" (stages 3 to 4), "Late" (stages 5 to 6), "Mature" (stages 7 to 8), and "Regression" stages (stage 9). In the testis, the classification groups analysed were "Early" (stages 6 and 7), "Midterm" (stages 8 to 9), "Late" (stages 1 to 2), "Mature" (stages 3 to 4), and "Spawning" stages (stage 5). For statistical purposes, to evaluate the maturation of females, an advanced maturation oocytes stage corresponding to the "Late" and "Mature" groups was created. The "Spawning" group was used in the statistical analysis to test differences in maturation in males.

### 2.3.6. Steroid Blood Analysis

After sampling, the blood samples were preserved at -80 °C in a storage solution (50% PBS, 30% heparin, 20% glycerol, 1% protein inhibitor). The quantification of blood 11-ketotestosterone and 17-β-estradiol in males was carried out by commercial Cayman ELISA kits (Cayman Chemical, Lexington KY, 11-ketotestosterone (Cat # 582751), 17-β-estradiol (Cat # 501890)), following the supplier's recommendations. These assays have been previously validated for fish blood samples [66,67]. The intra- and inter-assay variations reported were 13.4% and 7.7%, respectively, for a concentration of 12.5 pg·mL⁻¹ of 11-ketotestosterone and 33.8% and 27.8%, respectively, for a concentration of 2.4 pg·mL⁻¹ of 17-β-estradiol.

### 2.3.7. Gamete Quality, Reproduction and Fertilisation Rate

During the reproduction period (May-July), the fish exhibited secondary sexual traits such as a large belly showing the shape of the eggs and the red ventral pigmentation of the mature males [64]. To perform in vitro fertilisation (IVF), males were anesthetised and euthanised before extraction and dilaceration of their gonads in 200 µL of 300 mOsmol·kg⁻¹ Hank Balanced saline solution (HBSS). The semen was stored on ice until a gravid female fish was sampled, anaesthetised (0.125 g·L⁻¹ MS222) and smoothly stripped to extract the eggs. The eggs of each female were counted and fertilised with the sperm solution of 2 males. The sperm viability and spermatozoa concentration were carried out by flow cytometry (MACSQuant X, Miltenyi Biotec, Germany) by using the kit LIVE/DEAD™ Sperm Viability (L7011, ThermoFisher, USA). Briefly, 1 µL of SYBR™ 14 (50-fold diluted) and 1 µL of propidium iodide were incubated with 200 µL of spermatozoa suspension
(diluted 1/1000 with HBSS) for 10 min in the dark at room temperature. Then, the viable and dead cells expressed green (FL1) and red (FL3) fluorescence, respectively. The assessment of the fertilisation rate was carried out 24 h after fertilisation under a binocular magnifier (Discovery.V8 SteREO, ZEISS, Germany). The fertilised eggs were photographed (NIKON, D7200, Japan) with a binocular magnifier (Discovery.V8 SteREO, ZEISS, Germany) at 1 and 3 days post fertilisation (Dpf). Their diameter and the perivitelline index were measured individually using ImageJ software (Version 1.53e, Java 1.8.0_172, [58]) with the following formulas: Fertilisation rate = ((Fertilised eggs)/(Total of eggs)) × 100; and Perivitelline index = ((Perivitellin width)/(Diameter of the fertilised)) × 100.

### 3. Statistical Analysis

Data analyses were performed using the R software [68]. The significance threshold was set to *p*-value < 0.05 (See Supplementary Table S4). Prior to the ANOVA analysis, the assumptions of normality of the data and homogeneity of variance were tested by the Shapiro-Wilk test and Bartlett test. Levene's test was used to determine the homogeneity of variance when the assumption of normality was not fulfilled. If the assumption of homogeneity was met, the Kruskal-Wallis rank sum test followed by a pairwise Wilcoxon rank sum post hoc test was used, with Holm's *p*-value adjustment method [69]. If none of these assumptions for the ANOVA were met, the Welch's heteroscedastic F test, followed by the Games Howell post hoc test, was used.

Pearson's chi-squared test was used to determine the difference in mortality between the two scenarios.

The function glmmPQL from the r package "MASS" was used [70] to identify differences in gonadal maturation with the scenario as the fixed factor and the slide and the observer as random factors. Furthermore, the residuals normality and homogeneity of variance were checked graphically.

In addition, the analyses of the standard length, body mass, K<sub>corr</sub>, GSI, HIS and KSI were performed for the males and females separately.

### 4. Results

#### 4.1. Mortality and Biometric Study

Over the 7-month experimental period, mortality occurred homogeneously under both scenarios. A significantly higher mortality rate was observed in the RCP8.5 scenario (*p*-value = 0.01; Table 3). The final feed conversion ratio (FCR) was 24 mg<sub>food</sub>·mg<sub>body</sub><sup>-1</sup> in the RCP8.5 scenario and 20 mg<sub>food</sub>·mg<sub>body</sub><sup>-1</sup> in the Current scenario, with a mass thermal growth coefficient (mTGC) of 0.27 and 0.31 mg·degree days<sup>-1</sup>, respectively. Table S3 in Supplementary Materials shows the evolution of the stickleback FCR, TGC and mortality rates throughout the 7 months of the experiment.

<table><caption>**Table 3.** Stickleback mortality and biometric data after 7 months of acclimation to the RCP8.5 scenario. The asterisk indicates a significance between the scenarios.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Mortality (%)</th><th colspan="3">Standard Length (mm)</th><th colspan="3">Mass (g)</th></tr><tr><th>Mean ± S.D.</th><th>Min</th><th>Max</th><th>Mean ± S.D.</th><th>Min</th><th>Max</th></tr></thead><tbody><tr><td>Current</td><td>24.0</td><td>40.5 ± 5.8</td><td>22.1</td><td>55.0</td><td>1.4 ± 0.3</td><td>0.6</td><td>2.1</td></tr><tr><td>RCP8.5</td><td>31.8 *</td><td>41.3 ± 6.2</td><td>25.3</td><td>51.0</td><td>1.4 ± 0.4</td><td>0.8</td><td>2.5</td></tr></tbody></table>

At the beginning of the experiment (November), the mean standard length and the mass of the fish population were 26.3 ± 4.5 mm and 381.5 ± 191.7 mg, respectively (Figure 2). After 7 months of exposure to the RCP8.5 scenario, the mass of the fish increased significantly with time (*p*-value < 0.01) in a comparable manner in the two experimental scenarios (Supplementary Table S1), and no significant differences in the final body mass (*p*-value = 0.80) and standard length (*p*-value = 0.17) were found.
<figure><img src="image_7.png" /><figcaption>**Figure 2.** Stickleback mass values under Current (in light grey) and RCP8.5 (in dark grey) scenarios (*n* = 144 per scenario). The white box represents the mean mass of the population sample before the beginning of the experiment (*n* = 125). The 10th and 90th percentiles are represented by the whiskers, the 25th and 75th percentiles by the boxes, the median values by horizontal lines and outliers by points.</figcaption></figure>

At the sampling day, the individual body standard length and mass as well as the physiological indexes, including the corrected Fulton (Kcorr), gonado-somatic (GSI), hepato-somatic (HSI), and kidney-somatic (KSI) indexes, were estimated for each gender (Table 4). Male fish reared in the RCP8.5 scenario tended to have lower KSI (males: *p*-value = 0.06; females: *p*-value = 0.94) and GSI (males: *p*-value = 0.80; females: *p*-value = 0.37) than those reared in the Current scenario. There were no significant differences between scenario groups in the other biometric parameters measured, such as the standard length (males: *p*-value = 0.74; females: *p*-value = 0.09), body mass (males: *p*-value = 0.70; females: *p*-value = 0.35), Kcorr (males: *p*-value = 0.99; females: *p*-value = 0.25), and HSI (males: *p*-value = 0.32; females: *p*-value = 0.77).

<table><caption>**Table 4.** Stickleback standard body length and mass and physiological indexes at the sampling in July, after 7 months of acclimation to the RCP8.5 scenario. Kcorr: Fulton's index corrected of the mass of the gonads; GSI: gonado-somatic index; HSI: hepato-somatic index; KSI: kidney-somatic index (*n* = 5–9 per scenario and sex). Data are presented as mean ± S.D.</caption><thead><tr><th></th><th>Sex</th><th>Standard Length (mm)</th><th>Mass (g)</th><th>Kcorr</th><th>GSI</th><th>HSI</th><th>KSI</th></tr></thead><tbody><tr><td rowspan="2">Current</td><td>♀</td><td>45.0 ± 3.3</td><td>1.5 ± 0.2</td><td>1.4 ± 0.2</td><td>14.6 ± 7.4</td><td>4.9 ± 0.7</td><td>0.4 ± 0.3</td></tr><tr><td>♂</td><td>44.5 ± 2.4</td><td>1.3 ± 0.2</td><td>1.4 ± 0.1</td><td>0.6 ± 0.2</td><td>3.6 ± 0.7</td><td>2.8 ± 1.0</td></tr><tr><td rowspan="2">RCP8.5</td><td>♀</td><td>47.5 ± 2.1</td><td>1.6 ± 0.2</td><td>1.3 ± 0.1</td><td>11.1 ± 7.6</td><td>5.1 ± 1.8</td><td>0.4 ± 0.2</td></tr><tr><td>♂</td><td>44.9 ± 1.9</td><td>1.3 ± 0.2</td><td>1.4 ± 0.2</td><td>0.8 ± 0.5</td><td>4.4 ± 1.9</td><td>1.7 ± 0.7</td></tr></tbody></table>

## 4.2. Muscle Lipids Storage

The total lipid content of the white muscle of stickleback under each scenario (Figure 3) revealed a significantly lower percentage of lipids in the muscle of the females of the RCP8.5 scenario compared with females of the Current scenario and to males, irrespective of the scenario (*p*-value = 0.01; see Supplementary Table S2).
<figure><img src="image_8.png" /><figcaption>**Figure 3.** Stickleback muscle lipid storage after 7 months under Current (in light grey) and RCP8.5 (in dark grey) scenarios for female (white stripes) and male (no stripes) (n = 8). The 10th and 90th percentiles are represented by the whiskers, the 25th and 75th percentiles by the boxes, the median values by horizontal lines and outliers by points. Different letters indicate a significant difference between the boxes.</figcaption></figure>

### 4.3. Genes Expression Profiles in the Brain

The analysis of the gene expression in the median brain (Figure 4) revealed a downregulation of *cyp19a1b* in males under both Current and RCP8.5 scenarios (sex: *p*-value < 0.01). There were no significant differences between scenarios and sexes in all the other genes analysed (*gnrh2*: *p*-value = 0.24; *gnrh3*: *p*-value = 0.11; *kiss2*: *p*-value = 0.49; *gnih*: *p*-value = 0.42), although an overall trend suggesting an upregulation under RCP8.5 scenario is observed for all these genes related to reproductive processes. Furthermore, except for *cyp19a1b*, the analysis of the variance for data concerning the expressions of *gnrh2*, *gnrh3*, *kiss2* and *gnih* showed significantly higher interindividual variability in the RCP8.5 scenario (*gnrh2*: *p*-value = 0.02; *gnrh3*: *p*-value = 0.02; *kiss2*: *p*-value = 0.03; *gnih*: *p*-value = 0.01).

### 4.4. Gonad Maturation

The evaluation of the ovarian maturation by histology revealed a shift in the proportion of developmental stages between RCP8.5 and Current scenario (Figure 5A), with a significantly lower proportion of advanced maturation oocytes (Figure 5B) in RCP8.5 scenario (*p*-value = 0.01). In males, the analysis of testicular maturation (Figure 6B) did not reveal any difference in the proportion of tubules in the Spawning stage between the scenarios (*p*-value = 0.46). The Mature stage was predominant in the testis (Figure 6A) and in the ovary in both scenarios (Figure 5A). The histological structure of the gonadic tissues did not show any alterations.
<figure><img src="image_9.png" /><figcaption>**Figure 4.** Boxplots representing relative expression of *gnrh2* (A), *gnrh3* (B), *kiss2* (C), *gnih* (D) and *cyp19a1b* (E) in the median brain of sticklebacks after 7 months of exposure to Current (in light grey) and RCP8.5 (in dark grey) scenarios for female (white stripes) and male (no stripes) (females: *n* = 4 and 5 for Current and RCP8.5 scenarios, respectively; males: *n* = 5 and 7 for Current and RCP8.5, respectively). The 10th and 90th percentiles are represented by the whiskers, the 25th and 75th percentiles by the boxes, the median values by horizontal lines and outliers by points. Different letters indicate a significant difference between the boxes.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 5.** Proportions of oocytes in different maturation stages in stickleback ovaries (A) and proportion of advanced maturation oocytes (B) under Current (light grey) and RCP8.5 (dark grey) scenarios (*n* = 8). The 10th and 90th percentiles are represented by the whiskers, the 25th and 75th percentiles by the boxes, and the median values by horizontal lines. Different letters indicate a significant difference between the boxes.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 6.** Proportions of tubules in different maturation stages in the stickleback testis (A) and proportion of spawning tubules (B) under Current (light grey) and RCP8.5 (dark grey) scenarios (*n* = 8). The 10th and 90th percentiles are represented by the whiskers, the 25th and 75th percentiles by the boxes, the median values by horizontal lines and outliers by points.</figcaption></figure>

### 4.5. *Steroid Hormones*

The blood concentrations of 11-keto-testosterone (11-KT; Figure 7A) and 17-β-estradiol (E2; Figure 7B) at the July sampling point did not show any differences between scenarios in males (11-KT: *p*-value = 0.28; E2: *p*-value = 0.14). Furthermore, the variance of data revealed a trend in E2 blood concentration, showing higher interindividual variability in the RCP8.5 scenario (*p*-value = 0.06).

<figure><img src="image_12.png" /><figcaption>**Figure 7.** Blood concentrations of 11-keto-testosterone (11-KT; panel (A)) and 17-β-estradiol (E2; panel (B)) in stickleback males under Current (light grey; *n* = 17) and RCP8.5 (dark grey; *n* = 26) scenarios. The 10th and 90th percentiles are represented by the whiskers, the 25th and 75th percentiles by the boxes, the median values by horizontal lines and outliers by points.</figcaption></figure>

### 4.6. *Reproductive Success*

The number of eggs per fish obtained was 66 ± 23 for the Current scenario and 57 ± 18 for RCP8.5 (*p*-value = 0.28). The fertilisation rate was significantly lower in RCP8.5 scenario (63 ± 21.3%) than in the Current scenario (81.8 ± 4.5%; *p*-value = 0.02). In the scenario RCP8.5
at 1 and 3 Dpf, the size of the eggs was significantly smaller (scenario: *p*-value < 0.01, Dpf: *p*-value < 0.01; scenario/Dpf: *p*-value = 0.68; Figure 8A) and the perivitelline index significantly higher (scenario: *p*-value < 0.01, Dpf: *p*-value < 0.01; scenario/Dpf: *p*-value = 0.58; Figure 8B). The sperm quality was assessed, and no significant differences in spermatozoa concentration (*p*-value = 0.92) and mortality (*p*-value = 0.38) were observed between scenarios (Table 5).

<figure><img src="image_13.png" /><figcaption>**Figure 8.** Boxplots representing the diameter (A) and the perivitelline index (B) at 1 and 3 days post fertilisation (Dpf) under Current (light grey) and RCP8.5 (dark grey) scenarios (Current, Dpf 1 *n* = 19, Dpf 3 *n* = 14; RCP8.5, Dpf 1 *n* = 22, Dpf 3 *n* = 10). The 10th and 90th percentiles are represented by the whiskers, the 25th and 75th percentiles by the boxes, the median values by horizontal lines and outliers by points. Different letters indicate a significant difference between the boxes.</figcaption></figure>

<table><caption>**Table 5.** Summary of sperm quality markers of the three-spined stickleback after 7 months of acclimation to RCP8.5 (*n* = 8) or Current (*n* = 7) scenarios. Data are presented as mean ± S.D.</caption><thead><tr><th></th><th>Current</th><th>RCP8.5</th></tr></thead><tbody><tr><td>Spermatozoa concentration (10<sup>7</sup> cells/mL)</td><td>1.68 ± 1.25</td><td>1.62 ± 1.08</td></tr><tr><td>Spermatozoa mortality (%)</td><td>7.45 ± 2.70</td><td>5.81 ± 4.03</td></tr></tbody></table>

## 5. Discussion

Overall, the findings from this study indicate that a long acclimation to combined ocean warming and acidification, as predicted by the RCP8.5 IPCC scenario by 2100, could impact the survival and reproductive performance of three-spined stickleback.

Given the freshwater origin of the sticklebacks used in this study, one concern could be related to the additional stress the fish suffered in seawater, which could compromise their resilience to warming and acidification, although during the gradual acclimatisation to the salinity of Brest Harbour (34‰), no mortality was observed and the fish showed a behaviour and feeding rate consistent with and equivalent to individuals of the same stage in fresh water [71,72]. These observations are similar to the results obtained by Grotan et al. [73] on this same species. Their study, conducted on sticklebacks from populations locally adapted to saline, brackish and freshwater environments at three different salinities (0, 15 and 30‰), indicated that the standard metabolic rate did not increase in non-indigenous salinity environments. This means that marine sticklebacks paid no additional metabolic cost in freshwater, and freshwater sticklebacks paid no additional metabolic cost in marine waters. Altogether, it indicates that sticklebacks are able to rapidly acclimate to seawater, suggesting that it would not alter their physiological response to the exposure to a different climatic scenario.
In each scenario, the mortality was globally homogenous for the whole duration of the exposure, suggesting that it was likely natural mortality. No increase in mortality events was observed during the breeding season, meaning that the experimental set-up and fish density succeeded in avoiding the typical aggressive behaviour of dominant stickleback males. However, the cumulative mortality at the end of this experiment (from November to July) was higher in the RCP8.5 scenario, even if the lowest upper lethal temperature (21.6 °C) of stickleback (as previously reported) was not reached in this study [39]. Despite numerous studies on the effects of ocean acidification on the survival of early developmental stages of fish [74], no studies have examined the impact of acidification on stickleback survival at juvenile and adult stages. In our experimental conditions, we cannot attribute the increased vulnerability observed in sticklebacks directly to acidification conditions rather than to the interaction between acidification and warming. According to Pörtner [26], the physiological response to this multi-stress could require higher energy expenditure. Indeed, metabolic rates increase with temperature until the temperature level of intolerance for the species is reached [26]. This leads to an increase in food consumption and overall growth of individuals when food is not limiting [7,9,75], with a potential indirect impact on other key physiological functions.

In this study, biometric monitoring (body mass and standard length) during the 7 months of the experiment showed no difference between the fish subjected to the two experimental scenarios. However, the RCP8.5 scenario fish needed larger amounts of food to reach satiety. In this work, the RCP8.5 scenario represented a combination of warming and acidification, and the potential faster growth of fish reared in warmer water temperatures might have been counter-balanced by the costs associated with their simultaneous exposure to lower pH [76,77]. Moreover, and in contrast to many other studies, sticklebacks were not exposed to constant temperatures during this experiment [61], but rather to the thermal seasonal variations of the “Rade de Brest” in France, meaning that a positive effect at a given temperature might have been offset by a negative one (especially in warmers periods). This led to the same standard body size and mass in the two scenarios, but with a FCR of 24 $mg_{food}·mg_{body}^{-1}$ in the RCP8.5 scenario and of 20 $mg_{food}·mg_{body}^{-1}$ in the Current scenario. This finding suggests that at higher temperature and lower pH, the sticklebacks needed to increase their food intake to reach the same body size, indicating an extra energetic demand in the RCP8.5 scenario. For instance, Cominassi et al. [59] observed that in European sea bass, the decreased growth rate at higher temperature and $pCO_2$ levels was associated with a lower specific activity of digestive enzymes. Altogether, these findings suggest that sticklebacks facing warming and acidification would exhibit a lower FCR and/or higher metabolic cost. Previous studies hypothesised that a high $CO_2$ scenario leads to increased metabolic costs due to hypercapnia and the consequent acid-base regulation, which could compromise the energy available for growth and other physiological functions of the fish [26,78–80]. However, the effect of $CO_2$ on growth depends on the temperature [80–83]. To estimate how warming and acidification impact the energy absorption efficiency and metabolism, digestive enzyme activities should be determined, as well as oxygen consumption under swimming challenge.

The amount of reserves that fish can accumulate, particularly for reproduction, is likely to be impacted by such increases in FCR and metabolic costs. In fish, lipids, along with proteins, play a major role as energy sources for growth and reproduction [84]. Thus, their reproductive cycle includes consecutive phases of energy storage in the form of lipids, mobilisation, and transfer from the lipid deposition sites to the liver for metabolization and transport to the gonads [85]. In sticklebacks, the energetic requirements associated with reproduction may be met in two ways: stored energy or energy from food consumed as reproduction proceeds [86]. For instance, as the energy content of sticklebacks ovaries increases considerably between March and June, a parallel depletion in liver and carcass lipid and glycogen is observed [87]. In our study, the analysis of the total lipids in the muscle was carried out during the reproductive period. It showed that females under the RCP8.5 scenario had significantly lower lipid levels, suggesting a smaller amount of energy available for reproduction at the time of sampling.
This is in line with the observation of a lower proportion of mature oocytes in females under the RCP8.5 scenario. In July, during the spawning period, mature oocytes were the most common stage in both scenarios, but a smaller number of advanced maturation oocytes (late and mature stages) was observed in the RCP8.5 scenarios. It has already been reported that warmer environments can lead to a shift, a partial or even total inhibition in gametogenesis, depending on the species [24]. The slower oocyte maturation in the RCP8.5 scenario suggests that a delay and/or a partial inhibition of the maturation process might be occurring in these fish when warming occurred simultaneously with acidification. Keeping in mind our hypothesis concerning the increase in metabolic rate in the RCP8.5 scenario, we could assume that the reproductive capital of females under the RCP8.5 scenario is inadequate at the time of sampling to ensure the energy flow necessary for maturation [88,89]. What are the consequences of this potential delay? Fish exposed to the RCP8.5 scenario might reach their peak reproductive activity with a similar delay, which could have major consequences for their offspring if they hatch in suboptimal conditions. Future studies would be needed to understand the fate of the offspring of RCP8.5 fish, but such a mismatch (commonly reported in the literature [90–92]) and the decline in reproductive success, which itself comes after a significant decline in individuals’ survival, suggest that the RCP8.5 scenario could have major effects on the fitness of three-spined stickleback. Secondly, the quality of the eggs might be impaired, and this is what our data indicates. Indeed, egg quality also appeared to be affected in the RCP8.5 scenario, with eggs of smaller diameter and a higher perivitelline index at 1 and 3 Dpf. Similar observations, i.e., smaller diameters, have been reported by numerous studies on eggs of seawater sticklebacks with a constant warmer condition of 4 °C [46,47,93,94], and as reported by Glippa et al. [45], who showed a negative effect of temperature on stickleback egg diameter. Another study hypothesised lower egg quality in freshwater sticklebacks (i.e., small eggs) under a constant rise of 5 °C [61]. Other fish species also show smaller diameters after warming exposure (short winter), such as *Perca flavescens* [95] and *Gadus chalcogrammus*, under constant temperature (8 °C) increase [96]. Indeed, in many marine fish species, it is often observed that an increase in water temperature is associated with a decrease in egg size [97]. No effect of acidification has been found on the size of sticklebacks eggs [32,45], but in European eel (*Anguilla anguilla*), acidification induced a decline in both egg diameter and perivitelline space. The authors hypothesised that the entrance of acidified water into eggs could impact the development of the embryos and compromise their survival [98]. Then, if we assume the effects observed in European eel’s egg under acidification would occur also in stickleback, the reduction in egg diameter could be an additive or synergic effect of warming and acidification. Overall, our findings show smaller diameters of eggs with higher perivitelline index scores at 1 and 3 Dpf in the RCP8.5 scenario, suggesting lower availability of reserves for the embryo. This could potentially increase embryo vulnerability, which could result in suboptimal hatching conditions that substantially reduce the survival rate of larvae.

Males seemed to be less impacted by the RCP8.5 scenario, with muscle lipid levels and testicular maturation stable between scenarios. However, males' KSI scores tended to be lower in the RCP8.5 scenario after 7 months of exposure. Since no difference in body mass was observed between scenarios, our finding indicates a decrease in the mass of the kidney, which is in line with other studies showing similar results in male sticklebacks under warmer conditions [44,61,99]. Kidney hypertrophy is one of the secondary sexual characters of the males in the reproductive period of sticklebacks, which results from the kidney's physiological function in producing the spiggin, a protein glue that allows males to make the nest required for reproduction in nature. This hypertrophy is directly related to the androgen concentration in the blood and specifically 11-KT [100]. The potential hypotrophy of the male kidney under the RCP8.5 experimental scenario is not associated with a decrease in blood 11-KT, with blood 11-KT levels being similar in the two experimental groups at the sampling time. However, the same value for 11-KT concentration in the two scenarios may not correspond to the same point of the seasonal cycle of 11-KT levels in the blood. Therefore, we cannot exclude a delay or an inhibition of maturation in the RCP8.5 scenario, even in males. There were no significant differences in the concentration or mortality of
spermatozoa between scenarios. Nevertheless, we could not evaluate the direct or indirect effects of the climate scenario on the motility of spermatozoa, which is directly linked to the fertilisation potential [24,101,102]. An impact of the RCP8.5 scenario on sperm motility is likely, since reduction in fish sperm velocity has been already reported under higher temperature [13] and acidified water [102]. However, to our knowledge, the combined effect of acidification and warming on fish sperm motility is not clear yet. Overall, although the impacts of the RCP8.5 scenario on stickleback males were less clear than those on females, we cannot rule out a decline in their reproductive performance or at least a short delay in their maturation. To resolve this, monitoring of the kinetics of physiological indexes and hormone concentrations during the reproductive cycle would be needed. This would enable us to better decipher the impact of the RCP8.5 scenario in both sexes and, notably, to discriminate between a shift in and an inhibition of gonad maturation in this scenario. Similarly, kinetic assessments of lipids in the liver and the muscle could help to evaluate the energy status of fish in the context of climate change and contribute to distinguishing between potential delay or inhibition of maturation.

The neuroendocrine control of reproduction showed only a trend in the impact of the RCP8.5 scenario. The expression levels of genes involved in fish reproduction (*gnrh2*, *gnrh3*, *kiss2*, *gnih*) did not show any significant differences between the scenarios in males or females. However, there was a general trend suggesting a slight up-regulation in the RCP8.5 scenario in most genes. More specifically, the genes *gnrh2*, *gnrh3*, and *kiss2* in fish are all involved in the triggering of reproduction, and *gnrh3* is also involved in the maturation of fish [103–107]. Altogether, these data could suggest that sticklebacks under RCP8.5 would undergo central stimulation of the reproductive axis. This central stimulation could represent a positive feedback coming from the gonad and returning to the brain, aiming to increase the maturation of gonads in both sexes. This could be further evidence of a delay or a slight inhibition of fish maturation in this scenario in July at the time of the sampling. Interestingly, a previous study of the eurythermal sheepshead minnow (*Cyprinodon variegatus*) showed that the transcripts of genes encoding for *gnih* and *gnrh3* were higher in relative abundance in the hypothalamus of both sexes at higher temperature when a global inhibition of reproduction (reduced gonadal steroidogenesis and oogenesis) is observed [108]. The higher expression of *cyp19a1b* in female sticklebacks is not surprising, due to the role of *cyp19a1b* in converting androgens to oestrogens. Brain aromatase did not show any response to the experimental RCP8.5 scenario, which is in line with the lack of effect on blood sex steroid levels. This study also provided information on the reproductive success of sticklebacks, with lower fertilisation rates observed in the RCP8.5 scenario than in the Current scenario. This could be explained by a potential decline in gamete quality, for example, a possible decrease in spermatozoa motility and/or change in egg morphology (size and composition).

Interestingly, the RCP8.5 scenario exhibited increased variability in interindividual responses for most of the analysed end points related to reproduction, suggesting a potential amplification of the phenotypic variation in this scenario. Higher variance in the RCP8.5 scenario was observed for the expression levels of genes related to reproduction and blood sex steroids concentrations, notably 17-β-estradiol. Phenotypic plasticity is the capacity of a single genotype to express multiple phenotypes in response to environmental stimuli [109–111], allowing organisms to cope with environmental change [112]. Plasticity can occur within a generation and across generations (transgenerational plasticity). When rapid change in environmental conditions occurs, as in the context of climate change, within-generation plasticity can be important to mitigate the immediate impacts of the changing environment on fish individuals [94,113]. Most of the within-generation plasticity is driven by the developmental plasticity, defined as the phenotypic variation in later life stages explained by the environmental conditions experienced in early life [114,115]. In the present study, the exposure to environmental changes began when fish were already at the juvenile stage, excluding any chance of developmental plasticity. Nevertheless, the exposure lasted until adulthood and the gametogenesis process, which has been described as a
sensitivity window towards environmental changes that generates plasticity in fish [116]. The increased variation in hormonal responses and the neuroendocrine regulation of reproduction in the RCP8.5 scenario could lead, for instance, to a higher variation in spawning period starting and duration, with consequences for hatching time that could potentially impact offspring survival in the changed environment. Previous works have reported both increases [117] and decreases [118,119] in phenotypic variability in fish under warmer temperatures or acidification. According to the proposition of Rodriguez-Dominguez and collaborators [118], this discrepancy could be explained by a difference in phenotypic variation in response to changing environments in studies based on wild populations versus laboratory experiments. The authors noticed that most of environmental reduction in phenotype variation is observed in wild populations that face whole-life exposure to the changing environments, while increases in the phenotypic responses are observed mostly in laboratory experiments following short exposures. One explanation could be that in short-term laboratory experiments, the frequency of rare phenotypes that would unlikely survive in the wild is potentially increased [118]. The present study showed higher response variability in a laboratory experiment in which sticklebacks were facing warming and acidification at the same time after a relatively long exposure (7 months). The amplification of variations in phenotypes related to reproductive function is consistent with observations in experimental studies focusing only on warming impacts [117]. This could suggest that the stickleback population is still in the initial phase of selection for beneficial traits that provide advantages in the new environment. Evaluating the impact of changing environments on phenotypic plasticity is crucial to identifying fish acclimation ability in future environments.

## 6. Conclusions

The aim of this study was to estimate the global effect of the combination of warming (+3 °C) and acidification (–0.4 pH unit) on individuals’ mortality, growth, and reproductive performance in three-spined stickleback. The experimental design of this study enabled us to monitor the biometric parameters (size and mass) over 7 months and thus to estimate the growth of sticklebacks under the two scenarios from the juvenile stages to the first reproduction. Stickleback size and mass did not change with exposure to warming and acidification when feeding was not limited, but food consumption was greater under the RCP8.5 scenario. This suggests an increased need for energy to maintain homeostasis and other physiological functions. The measurement of lipid content in muscle, an indicator of energy reserves, revealed lower levels in females in the RCP8.5 scenario, supporting the hypothesis of a divergence in energy allocation strategy between scenarios. This experiment also enabled us to characterise the first reproduction of sticklebacks in the two different scenarios. After 7 months of exposure, the females showed lower maturation and egg quality under the RCP8.5 scenario, which could have contributed to the lower fertilisation rate observed. This suggests a delay and/or inhibition of gametogenesis and maturation in warmed and acidified waters. The analysis of blood sex steroid concentrations and physiological indexes at one single time point did not allow us to discriminate between a delay and an inhibition of maturation in the RCP8.5 scenario. Altogether, these results highlight the importance of studying the effects of combined stresses, since resilience and physiological responses of organisms to multi-stress can differ from those observed when organisms of the same species face just one stress at a time. This experiment is not based on a complete cross-factorial design and does not allow to distinguish the effects of single stressors and the potential interaction between them, but instead, shows the resulting global effect of the two stresses when applied together.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/fishes8100523/s1. Supplementary Table S1: Pairwise Games Howell post hoc test on stickleback growth in mass between Current and RCP8.5 condition through the 7 months of acclimation; Supplementary Table S2: Pairwise Wilcoxon post hoc test on stickleback fresh lipids content between sexes and experimental scenarios at the final sampling in July after 7 months
of acclimation to RCP8.5 scenario; Supplementary Table S3: Sticklebacks feed conversion ratio (FCR), thermal growth coefficients (TGC) and mortality rate through the 7 months of experiments; Supplementary Table S4: Statistical tests and scores of data analysis.

**Author Contributions:** Conceptualization, V.L. and A.S.; Data curation, J.D., A.B.-N., S.C., O.M., U.I., M.-M.L.G., L.M. and C.T.; Formal analysis, J.D., V.L., C.L. and A.S.; Funding acquisition, V.L. and A.S.; Methodology, J.D., V.L., C.L., A.B.-N., S.C., M.-M.L.G., C.T. and A.S.; Project administration, V.L. and A.S.; Validation, C.L.; Writing—original draft, J.D.; Writing—review and editing, V.L., C.L., A.B.-N., O.M. and A.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the French national program EC2CO-DYCOVI through the ECHANGE project to V.L. (2021–2023) and Ifremer (ICEfish project). J.D. PhD grant was funded by "Université de Bretagne Occidental" and by the ISblue project, Interdisciplinary graduate school for the blue planet (ANR-17-EURE-0015) and co-funded by a grant from the French government under the program "Investissements d'Avenir" embedded in France 2030.

**Institutional Review Board Statement:** All animal experiments in this study were carried out in accordance with French national laws and EU Directive 2010/63/EU and the European Commission recommendation 2007/526/EC regarding the welfare of animals used for scientific purposes. The experimental design was authorized by the Regional Ethics Committee for animal testing (CEFEA: Comité d'Éthique Finistérien en Expérimentation Animale, registering code C2EA-74) and by the French Ministère de l'Enseignement Supérieur de la Recherche et de l'Innovation (Authorization APAFIS #27510, permit number 2020100911422298_v3).

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available in the Supplementary Materials.

**Acknowledgments:** Special thanks to all the members of the PHYTNESS unit (Ifremer-LEMAR) for their assistance in the experimental work and fish maintenance.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Pörtner, H.-O.; Karl, D.M.; Boyd, P.W.; Cheung, W.W.L.; Lluch-Cota, S.E.; Nojiri, Y.; Schmidt, D.N.; Zavialov, P.O. Ocean Systems. In *Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change*; Field, C.B., Barros, V.R., Dokken, D.J., Mach, K.J., Mastrandrea, M.D., Bilir, T.E., Chatterjee, M., Ebi, K.L., Estrada, Y.O., Genova, R.C., et al., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2014; pp. 411–484, ISBN 978-1-107-41537-9.

2. Kwiatkowski, L.; Torres, O.; Bopp, L.; Aumont, O.; Chamberlain, M.; Christian, J.R.; Dunne, J.P.; Gehlen, M.; Ilyina, T.; John, J.G.; et al. Twenty-First Century Ocean Warming, Acidification, Deoxygenation, and Upper-Ocean Nutrient and Primary Production Decline from CMIP6 Model Projections. *Biogeosciences* **2020**, *17*, 3439–3470. [CrossRef]

3. Grorud-Colvert, K.; Sponaugle, S. Variability in Water Temperature Affects Trait-Mediated Survival of a Newly Settled Coral Reef Fish. *Oecologia* **2011**, *165*, 675–686. [CrossRef] [PubMed]

4. Hermelink, B.; Wuertz, S.; Trubiroha, A.; Rennert, B.; Kloas, W.; Schulz, C. Influence of Temperature on Puberty and Maturation of Pikeperch, Sander Lucioperca. *Gen. Comp. Endocrinol.* **2011**, *172*, 282–292. [CrossRef] [PubMed]

5. Tierney, K.B.; Patterson, D.A.; Kennedy, C.J. The Influence of Maternal Condition on Offspring Performance in Sockeye Salmon *Oncorhynchus nerka*. *J. Fish. Biol.* **2009**, *75*, 1244–1257. [CrossRef] [PubMed]

6. Fry, F.E.J. The Effect of Environmental Factors on the Physiology of Fish. In *Fish Physiology*; Elsevier: Amsterdam, The Netherlands, 1971; Volume 6, pp. 1–98, ISBN 978-0-12-350406-7.

7. Catalán, I.A.; Auch, D.; Kamermans, P.; Morales-Nin, B.; Angelopoulos, N.V.; Reglero, P.; Sandersfeld, T.; Peck, M.A. Critically Examining the Knowledge Base Required to Mechanistically Project Climate Impacts: A Case Study of Europe's Fish and Shellfish. *Fish. Fish.* **2019**, *20*, 501–517. [CrossRef]

8. Cominassi, L.; Moyano, M.; Claireaux, G.; Howald, S.; Mark, F.C.; Zambonino-Infante, J.-L.; Le Bayon, N.; Peck, M.A. Combined Effects of Ocean Acidification and Temperature on Larval and Juvenile Growth, Development and Swimming Performance of European Sea Bass (*Dicentrarchus labrax*). *PLoS ONE* **2019**, *14*, e0221283. [CrossRef]

9. DeCourten, B.; Romney, A.; Brander, S. The Heat Is On: Complexities of Aquatic Endocrine Disruption in a Changing Global Climate. In *Separation Science and Technology*; Elsevier: Amsterdam, The Netherlands, 2019; Volume 11, pp. 13–49, ISBN 978-0-12-815730-5.

10. Kim, S.-Y.; Metcalfe, N.B.; da Silva, A.; Velando, A. Thermal Conditions during Early Life Influence Seasonal Maternal Strategies in the Three-Spined Stickleback. *BMC Ecol.* **2017**, *17*, 34. [CrossRef]
11. Oliveira, H.; Maulvault, A.L.; Castanho, S.; Repolho, T.; Valente, L.M.P.; Pousão-Ferreira, P.; Rosa, R.; Marques, A.; Anacleto, P. Lack of Detrimental Effects of Ocean Acidification and Warming on Proximate Composition, Fitness and Energy Budget of Juvenile Senegalese Sole (*Solea senegalensis*). *Sci. Total Environ.* **2023**, *857*, 159491. [CrossRef]

12. Sampaio, E.; Santos, C.; Rosa, I.C.; Ferreira, V.; Pörtner, H.-O.; Duarte, C.M.; Levin, L.A.; Rosa, R. Impacts of Hypoxic Events Surpass Those of Future Ocean Warming and Acidification. *Nat. Ecol. Evol.* **2021**, *5*, 311–321. [CrossRef]

13. Fenkes, M.; Fitzpatrick, J.L.; Ozolina, K.; Shiels, H.A.; Nudds, R.L. Sperm in Hot Water: Direct and Indirect Thermal Challenges Interact to Impact on Brown Trout Sperm Quality. *J. Exp. Biol.* **2017**, *220*, 2513–2520. [CrossRef]

14. Lim, B.-S.; Kagawa, H.; Gen, K.; Okuzawa, K. Effects of Water Temperature on the Gonadal Development and Expression of Steroidogenic Enzymes in the Gonad of Juvenile Red Seabream, *Pagrus major*. *Fish. Physiol. Biochem.* **2003**, *28*, 161–162. [CrossRef]

15. Pankhurst, N.W.; Purser, G.J.; Van Der Kraak, G.; Thomas, P.M.; Forteath, G.N.R. Effect of Holding Temperature on Ovulation, Egg Fertility, Plasma Levels of Reproductive Hormones and In Vitro Ovarian Steroidogenesis in the Rainbow Trout *Oncorhynchus mykiss*. *Aquaculture* **1996**, *146*, 277–290. [CrossRef]

16. Elisio, M.; Chalde, T.; Miranda, L.A. Effects of Short Periods of Warm Water Fluctuations on Reproductive Endocrine Axis of the Pejerrey (*Odontesthes bonariensis*) Spawning. *Comp. Biochem. Physiol. Part. A Mol. Integr. Physiol.* **2012**, *163*, 47–55. [CrossRef] [PubMed]

17. Anderson, K.; Pankhurst, N.; King, H.; Elizur, A. Effect of Thermal Challenge on the Expression of Genes Involved in Ovarian Steroidogenesis in Tasmanian Atlantic Salmon (*Salmo salar*). *Aquaculture* **2017**, *479*, 474–478. [CrossRef]

18. Dionísio, G.; Campos, C.; Valente, L.M.P.; Conceição, L.E.C.; Cancela, M.L.; Gavaia, P.J. Effect of Egg Incubation Temperature on the Occurrence of Skeletal Deformities in Solea Senegalensis. *J. Appl. Ichthyol.* **2012**, *28*, 471–476. [CrossRef]

19. Schulte, P.M.; Healy, T.M.; Fangue, N.A. Thermal Performance Curves, Phenotypic Plasticity, and the Time Scales of Temperature Exposure. *Integr. Comp. Biol.* **2011**, *51*, 691–702. [CrossRef]

20. González, A.; Piferrer, F. Characterization of Aromatase Activity in the Sea Bass: Effects of Temperature and Different Catalytic Properties of Brain and Ovarian Homogenates and Microsomes. *J. Exp. Zool.* **2002**, *293*, 500–510. [CrossRef]

21. Anderson, K.; Luckenbach, J.A.; Yamamoto, Y.; Elizur, A. Impacts of Fsh, Igf1, and High Temperature on the Expression of Genes Involved in Steroidogenesis, Cell Communication, and Apoptosis in Isolated Coho Salmon Previtellogenic Ovarian Follicles. *Aquaculture* **2019**, *506*, 60–69. [CrossRef]

22. Anderson, K.; Swanson, P.; Pankhurst, N.; King, H.; Elizur, A. Effect of Thermal Challenge on Plasma Gonadotropin Levels and Ovarian Steroidogenesis in Female Maiden and Repeat Spawning Tasmanian Atlantic Salmon (*Salmo salar*). *Aquaculture* **2012**, *334*–337, 205–212. [CrossRef]

23. Tenugu, S.; Pranoty, A.; Mamta, S.-K.; Senthilkumaran, B. Development and Organisation of Gonadal Steroidogenesis in Bony Fishes—A Review. *Aquac. Fish.* **2021**, *6*, 223–246. [CrossRef]

24. Alix, M.; Kjesbu, O.S.; Anderson, K.C. From Gametogenesis to Spawning: How Climate-driven Warming Affects Teleost Reproductive Biology. *J. Fish. Biol.* **2020**, *97*, 607–632. [CrossRef] [PubMed]

25. Servili, A.; Canario, A.V.M.; Mouchel, O.; Muñoz-Cueto, J.A. Climate Change Impacts on Fish Reproduction Are Mediated at Multiple Levels of the Brain-Pituitary-Gonad Axis. *Gen. Comp. Endocrinol.* **2020**, *291*, 113439. [CrossRef] [PubMed]

26. Pörtner, H.O. Synergistic Effects of Temperature Extremes, Hypoxia, and Increases in CO₂ on Marine Animals: From Earth History to Global Change. *J. Geophys. Res.* **2005**, *110*, C09S10. [CrossRef]

27. Nilsson, G.E.; Dixson, D.L.; Domenici, P.; McCormick, M.I.; Sørensen, C.; Watson, S.-A.; Munday, P.L. Near-Future Carbon Dioxide Levels Alter Fish Behaviour by Interfering with Neurotransmitter Function. *Nat. Clim. Chang.* **2012**, *2*, 201–204. [CrossRef]

28. Tresguerres, M.; Hamilton, T.J. Acid–Base Physiology, Neurobiology and Behaviour in Relation to CO₂-Induced Ocean Acidification. *J. Exp. Biol.* **2017**, *220*, 2136–2148. [CrossRef]

29. Frommel, A.Y.; Maneja, R.; Lowe, D.; Malzahn, A.M.; Geffen, A.J.; Folkvord, A.; Piatkowski, U.; Reusch, T.B.H.; Clemmesen, C. Severe Tissue Damage in Atlantic Cod Larvae under Increasing Ocean Acidification. *Nat. Clim. Chang.* **2012**, *2*, 42–46. [CrossRef]

30. Enzor, L.A.; Zippay, M.L.; Place, S.P. High Latitude Fish in a High CO₂ World: Synergistic Effects of Elevated Temperature and Carbon Dioxide on the Metabolic Rates of Antarctic Notothenioids. *Comp. Biochem. Physiol. Part A Mol. Integr. Physiol.* **2013**, *164*, 154–161. [CrossRef]

31. Pimentel, M.S.; Faleiro, F.; Marques, T.; Bispo, R.; Dionísio, G.; Faria, A.M.; Machado, J.; Peck, M.A.; Pörtner, H.; Pousão-Ferreira, P.; et al. Foraging Behaviour, Swimming Performance and Malformations of Early Stages of Commercially Important Fishes under Ocean Acidification and Warming. *Clim. Chang.* **2016**, *137*, 495–509. [CrossRef]

32. Schade, F.M.; Clemmesen, C.; Mathias Wegner, K. Within- and Transgenerational Effects of Ocean Acidification on Life History of Marine Three-Spined Stickleback (*Gasterosteus aculeatus*). *Mar. Biol.* **2014**, *161*, 1667–1676. [CrossRef]

33. Munday, P.L.; Welch, M.J.; Allan, B.J.M.; Watson, S.-A.; McMahon, S.J.; McCormick, M.I. Effects of Elevated CO₂ on Predator Avoidance Behaviour by Reef Fishes Is Not Altered by Experimental Test Water. *PeerJ* **2016**, *4*, e2501. [CrossRef]

34. Pimentel, M.S.; Faleiro, F.; Dionísio, G.; Repolho, T.; Pousão-Ferreira, P.; Machado, J.; Rosa, R. Defective Skeletogenesis and Oversized Otoliths in Fish Early Stages in a Changing Ocean. *J. Exp. Biol.* **2014**, *217*, 2062–2070. [CrossRef] [PubMed]

35. Pimentel, M.S.; Faleiro, F.; Diniz, M.; Machado, J.; Pousão-Ferreira, P.; Peck, M.A.; Pörtner, H.O.; Rosa, R. Oxidative Stress and Digestive Enzyme Activity of Flatfish Larvae in a Changing Ocean. *PLoS ONE* **2015**, *10*, e0134082. [CrossRef] [PubMed]
36. Dahlke, F.; Lucassen, M.; Bickmeyer, U.; Wohlrab, S.; Puvanendran, V.; Mortensen, A.; Chierici, M.; Pörtner, H.-O.; Storch, D. Fish Embryo Vulnerability to Combined Acidification and Warming Coincides with a Low Capacity for Homeostatic Regulation. *J. Exp. Biol.* **2020**, 223, jeb212589. [CrossRef] [PubMed]

37. Villalobos, C.; Love, B.A.; Olson, M.B. Ocean Acidification and Ocean Warming Effects on Pacific Herring (*Clupea pallasi*) Early Life Stages. *Front. Mar. Sci.* **2020**, *7*, 597899. [CrossRef]

38. Pimentel, M.S.; Faleiro, F.; Machado, J.; Pousão-Ferreira, P.; Rosa, R. Seabream Larval Physiology under Ocean Warming and Acidification. *Fishes* **2020**, *5*, 1. [CrossRef]

39. Wootton, R.J. *A Functional Biology of Sticklebacks*; Springer US: Boston, MA, 1984; ISBN 978-1-4615-8515-2.

40. Blaker, E.; Sebire, M.; Ellis, T.; Katsiadaki, I. The Housing, Care, and Use of a Laboratory Three-Spined Stickleback Colony. In *Laboratory Fish in Biomedical Research*; Elsevier: Amsterdam, The Netherlands, 2022; pp. 349–371, ISBN 978-0-12-821099-4.

41. *Biology of the Three-Spined Stickleback*; Ostlund-Nilsson, S.; Mayer, I.; Huntingford, F.A. (Eds.) CRC Press: Boca Raton, FL, USA, 2006; ISBN 978-0-429-12283-5.

42. Jones, F.C.; Grabherr, M.G.; Chan, Y.F.; Russell, P.; Mauceli, E.; Johnson, J.; Swofford, R.; Pirun, M.; Zody, M.C.; White, S.; et al. The Genomic Basis of Adaptive Evolution in Threespine Sticklebacks. *Nature* **2012**, *484*, 55–61. [CrossRef]

43. Kingsley, D.M.; Zhu, B.; Osoegawa, K.; de Jong, P.; Schein, J.; Marra, M.; Peichel, C.; Amemiya, C.; Schlüter, D.; Balabhadra, S.; et al. New Genomic Tools for Molecular Studies of Evolutionary Change in Threespine Sticklebacks. *Behav. 2004*, *141*, 1331–1344. [CrossRef]

44. Guderley, H.; Lavoie, B.A.; Dubois, N. The Interaction among Age, Thermal Acclimation and Growth Rate in Determining Muscle Metabolic Capacities and Tissue Masses in the Threespine Stickleback, *Gasterosteus aculeatus*. *Fish. Physiol. Biochem.* **1994**, *13*, 419–431. [CrossRef]

45. Glippa, O.; Brutemark, A.; Johnson, J.; Spilling, K.; Candolin, U.; Engström-Öst, J. Early Development of the Threespine Stickleback in Relation to Water PH. *Front. Mar. Sci.* **2017**, *4*, 427. [CrossRef]

46. Shama, L.N.S. Bet Hedging in a Warming Ocean: Predictability of Maternal Environment Shapes Offspring Size Variation in Marine Sticklebacks. *Glob. Chang. Biol.* **2015**, *21*, 4387–4400. [CrossRef]

47. Shama, L.N.S. The Mean and Variance of Climate Change in the Oceans: Hidden Evolutionary Potential under Stochastic Environmental Variability in Marine Sticklebacks. *Sci. Rep.* **2017**, *7*, 8889. [CrossRef] [PubMed]

48. Kim, S.-Y.; Costa, M.M.; Esteve-Codina, A.; Velando, A. Transcriptional Mechanisms Underlying Life-History Responses to Climate Change in the Three-Spined Stickleback. *Evol. Appl.* **2017**, *10*, 718–730. [CrossRef] [PubMed]

49. Gidden, M.J.; Riahi, K.; Smith, S.J.; Fujimori, S.; Luderer, G.; Kriegler, E.; van Vuuren, D.P.; van den Berg, M.; Feng, L.; Klein, D.; et al. Global Emissions Pathways under Different Socioeconomic Scenarios for Use in CMIP6: A Dataset of Harmonized Emissions Trajectories through the End of the Century. *Geosci. Model Dev.* **2019**, *12*, 1443–1475. [CrossRef]

50. Meinshausen, M.; Nicholls, Z.R.J.; Lewis, J.; Gidden, M.J.; Vogel, E.; Freund, M.; Beyerle, U.; Gessner, C.; Nauels, A.; Bauer, N.; et al. The Shared Socio-Economic Pathway (SSP) Greenhouse Gas Concentrations and Their Extensions to 2500. *Geosci. Model Dev.* **2020**, *13*, 3571–3605. [CrossRef]

51. Cohen-Rengifo, M.; Danion, M.; Gonzalez, A.-A.; Bégout, M.-L.; Cormier, A.; Noël, C.; Cabon, J.; Vitré, T.; Mark, F.C.; Mazurais, D. The Extensive Transgenerational Transcriptomic Effects of Ocean Acidification on the Olfactory Epithelium of a Marine Fish Are Associated with a Better Viral Resistance. *BMC Genom.* **2022**, *23*, 448. [CrossRef]

52. Anderson, D.H.; Robinson, R.J. Rapid Electrometric Determination of Alkalinity of Sea Water Using Glass Electrode. *Ind. Eng. Chem. Anal. Ed.* **1946**, *18*, 767–769. [CrossRef]

53. Strickland, J.D.H.; Parsons, T.R. *A Practical Handbook of Seawater Analysis*, 2nd ed.; Fisheries Research Board of Canada: Ottawa, ON, Canada, 1972. [CrossRef]

54. Lewis, E.; Wallace, D.W.R. MS Excel Program Developed for CO₂ System Calculations; Carbon Dioxide Information Analysis Center (CDIAC), Environmental Sciences Division: Oak Ridge, TN, USA, 2012.

55. Mehrbach, C.; Culberson, C.H.; Hawley, J.E.; Pytkowicz, R.M. Measurement of the Apparent Dissociation Constants of Carbonic Acid in Seawater at Atmospheric Pressure. *Linnol. Oceanogr.* **1973**, *18*, 897–907. [CrossRef]

56. Dickson, A.G.; Millero, F.J. A Comparison of the Equilibrium Constants for the Dissociation of Carbonic Acid in Seawater Media. *Deep. Sea Res. Part A Oceanogr. Res. Pap.* **1987**, *34*, 1733–1743. [CrossRef]

57. Xie, C.; Li, Q.; Han, G.; Liu, H.; Yang, J.; Li, J. Stable Isotope Dilution Assay for the Accurate Determination of Tricaine in Fish Samples by HPLC-MS-MS. *Biomed. Chromatogr.* **2019**, *33*, e4512. [CrossRef]

58. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 Years of Image Analysis. *Nat. Methods* **2012**, *9*, 671–675. [CrossRef]

59. Cominassi, L.; Moyano, M.; Claireaux, G.; Howald, S.; Mark, F.C.; Zambonino-Infante, J.-L.; Peck, M.A. Food Availability Modulates the Combined Effects of Ocean Acidification and Warming on Fish Growth. *Sci. Rep.* **2020**, *10*, 2338. [CrossRef] [PubMed]

60. Iwama, G.K.; Tautz, A.F. A Simple Growth Model for Salmonids in Hatcheries. *Can. J. Fish. Aquat. Sci.* **1981**, *38*, 649–656. [CrossRef]

61. Hani, Y.M.I.; Turies, C.; Palluel, O.; Delahaut, L.; Bado-Nilles, A.; Geffard, A.; Dedourge-Geffard, O.; Porcher, J.-M. Effects of a Chronic Exposure to Different Water Temperatures and/or to an Environmental Cadmium Concentration on the Reproduction of the Threespine Stickleback (*Gasterosteus aculeatus*). *Ecotoxicol. Environ. Saf.* **2019**, *174*, 48–57. [CrossRef] [PubMed]
62. Folch, J.; Lees, M.; Sloane Stanley, G.H. A Simple Method for the Isolation and Purification of Total Lipides from Animal Tissues. *J. Biol. Chem.* **1957**, *226*, 497–509. [CrossRef] [PubMed]

63. Untergasser, A.; Nijveen, H.; Rao, X.; Bisseling, T.; Geurts, R.; Leunissen, J.A.M. Primer3Plus, an Enhanced Web Interface to Primer3. *Nucleic Acids Res.* **2007**, *35*, W71–W74. [CrossRef]

64. Sokolowska, E.; Kulczykowska, E. Annual Reproductive Cycle in Two Free Living Populations of Three-Spined Stickleback (*Gasterosteus aculeatus* L.): Patterns of Ovarian and Testicular Development. *Oceanologia* **2006**, *48*, 103–124.

65. Furin, C.G.; von Hippel, F.A.; Postlethwait, J.H.; Buck, C.L.; Cresko, W.A.; O'Hara, T.M. Developmental Timing of Sodium Perchlorate Exposure Alters Angiogenesis, Thyroid Follicle Proliferation and Sexual Maturation in Stickleback. *Gen. Comp. Endocrinol.* **2015**, *219*, 24–35. [CrossRef]

66. Gardell, A.M.; Dillon, D.M.; Smayda, L.C.; von Hippel, F.A.; Cresko, W.A.; Postlethwait, J.H.; Buck, C.L. Perchlorate Exposure Does Not Modulate Temporal Variation of Whole-Body Thyroid and Androgen Hormone Content in Threespine Stickleback. *Gen. Comp. Endocrinol.* **2015**, *219*, 45–52. [CrossRef]

67. Kitano, J.; Kakioka, R.; Ishikawa, A.; Toyoda, A.; Kusakabe, M. Differences in the Contributions of Sex Linkage and Androgen Regulation to Sex-biased Gene Expression in Juvenile and Adult Sticklebacks. *J. Evol. Biol.* **2020**, *33*, 1129–1138. [CrossRef]

68. R Core Team R: A Language and Environment for Statistical Computing 2018, Vienna, Austria. Available online: https://www.R-project.org/ (accessed on 30 September 2023).

69. Holm, S. A Simple Sequentially Rejective Multiple Test Procedure. *Scand. J. Stat.* **1979**, *6*, 65–70.

70. Venables, W.N.; Ripley, B.D. *Modern Applied Statistics with S*; Statistics and Computing; Springer: New York, NY, USA, 2002; ISBN 978-1-4419-3008-8.

71. Leloutre, C.; Péry, A.R.R.; Porcher, J.-M.; Beaudouin, R. A Bioenergetics Model of the Entire Life Cycle of the Three-Spined Stickleback, *Gasterosteus aculeatus*. *Ecol. Freshw. Fish.* **2018**, *27*, 116–127. [CrossRef]

72. Hani, Y.M.I.; Turies, C.; Palluel, O.; Delahaut, L.; Gaillet, V.; Bado-nilles, A.; Porcher, J.-M.; Geffard, A.; Dedourge-geffard, O. Effects of Chronic Exposure to Cadmium and Temperature, Alone or Combined, on the Threespine Stickleback (*Gasterosteus aculeatus*): Interest of Digestive Enzymes as Biomarkers. *Aquat. Toxicol.* **2018**, *199*, 252–262. [CrossRef] [PubMed]

73. Grøtan, K.; Østbye, K.; Taugbøl, A.; Vøllestad, L.A. No Short-Term Effect of Salinity on Oxygen Consumption in Threespine Stickleback (*Gasterosteus aculeatus*) from Fresh, Brackish, and Salt Water. *Can. J. Zool.* **2012**, *90*, 1386–1393. [CrossRef]

74. Esbaugh, A.J. Physiological Implications of Ocean Acidification for Marine Fish: Emerging Patterns and New Insights. *J. Comp. Physiol. B* **2018**, *188*, 1–13. [CrossRef] [PubMed]

75. Baag, S.; Mandal, S. Combined Effects of Ocean Warming and Acidification on Marine Fish and Shellfish: A Molecule to Ecosystem Perspective. *Sci. Total Environ.* **2022**, *802*, 149807. [CrossRef] [PubMed]

76. Gräns, A.; Jutfelt, F.; Sandblom, E.; Jönsson, E.; Wiklander, K.; Seth, H.; Olsson, C.; Dupont, S.; Ortega-Martinez, O.; Einarsdottir, I.; et al. Aerobic Scope Fails to Explain the Detrimental Effects on Growth Resulting from Warming and Elevated CO₂ in Atlantic Halibut. *J. Exp. Biol.* **2014**, *217*, 711–717. [CrossRef] [PubMed]

77. Pistevos, J.C.A.; Nagelkerken, I.; Rossi, T.; Olmos, M.; Connell, S.D. Ocean Acidification and Global Warming Impair Shark Hunting Behaviour and Growth. *Sci. Rep.* **2015**, *5*, 16293. [CrossRef] [PubMed]

78. Heuer, R.M.; Grosell, M. Physiological Impacts of Elevated Carbon Dioxide and Ocean Acidification on Fish. *Am. J. Physiol. -Regul. Integr. Comp. Physiol.* **2014**, *307*, R1061–R1084. [CrossRef]

79. Heuer, R.M.; Grosell, M. Elevated CO₂ Increases Energetic Cost and Ion Movement in the Marine Fish Intestine. *Sci. Rep.* **2016**, *6*, 34480. [CrossRef]

80. Munday, P.; Crawley, N.; Nilsson, G. Interacting Effects of Elevated Temperature and Ocean Acidification on the Aerobic Performance of Coral Reef Fishes. *Mar. Ecol. Prog. Ser.* **2009**, *388*, 235–242. [CrossRef]

81. Concannon, C.A.; Cross, E.L.; Jones, L.F.; Murray, C.S.; Matassa, C.M.; McBride, R.S.; Baumann, H. Temperature-Dependent Effects on Fecundity in a Serial Broadcast Spawning Fish after Whole-Life High CO₂ Exposure. *ICES J. Mar. Sci.* **2021**, *78*, 3724–3734. [CrossRef]

82. Gobler, C.J.; Merlo, L.R.; Morrell, B.K.; Griffith, A.W. Temperature, Acidification, and Food Supply Interact to Negatively Affect the Growth and Survival of the Forage Fish, Menidia Beryllina (Inland Silverside), and Cyprinodon Variegatus (Sheepshead Minnow). *Front. Mar. Sci.* **2018**, *5*, 86. [CrossRef]

83. Murray, C.; Baumann, H. You Better Repeat It: Complex CO₂ × Temperature Effects in Atlantic Silverside Offspring Revealed by Serial Experimentation. *Diversity* **2018**, *10*, 69. [CrossRef]

84. Tocher, D.R. Metabolism and Functions of Lipids and Fatty Acids in Teleost Fish. *Rev. Fish. Sci.* **2003**, *11*, 107–184. [CrossRef]

85. Craig, S.R.; MacKenzie, D.S.; Jones, G.; Gatlin, D.M. Seasonal Changes in the Reproductive Condition and Body Composition of Free-Ranging Red Drum, *Sciaenops ocellatus*. *Aquaculture* **2000**, *190*, 89–102. [CrossRef]

86. Wootton, R.J.; Fletcher, D.A. Effect of Spawning Number and Ration on Reproductive Performance of the Batch-Spawning Three-Spined Stickleback *Gasterosteus aculeatus*. *J. Fish. Biol.* **2009**, *75*, 618–629. [CrossRef]

87. Wootton, R.J. Energy Cost of Egg Production and Environmental Determinants of Fecundity in Teleost Fishes. *Symp. Zool. Soc. Lond.* **1979**, *44*, 133–155.

88. Baker, J.A.; Wund, M.A.; Heins, D.C.; King, R.W.; Reyes, M.L.; Foster, S.A. Life-History Plasticity in Female Threespine Stickleback. *Heredity* **2015**, *115*, 322–334. [CrossRef] [PubMed]
89. Madsen, T.; Shine, R. The Adjustment of Reproductive Threshold to Prey Abundance in a Capital Breeder. *J. Anim. Ecol.* **1999**, *68*, 571–580. [CrossRef]

90. Asch, R.G.; Stock, C.A.; Sarmiento, J.L. Climate Change Impacts on Mismatches between Phytoplankton Blooms and Fish Spawning Phenology. *Glob. Chang. Biol.* **2019**, *25*, 2544–2559. [CrossRef] [PubMed]

91. Cohen, J.M.; Lajeunesse, M.J.; Rohr, J.R. A Global Synthesis of Animal Phenological Responses to Climate Change. *Nat. Clim. Chang.* **2018**, *8*, 224–228. [CrossRef]

92. Durant, J.; Hjermann, D.; Falkenhaug, T.; Gifford, D.; Naustvoll, L.; Sullivan, B.; Beaugrand, G.; Stenseth, N. Extension of the Match-Mismatch Hypothesis to Predator-Controlled Systems. *Mar. Ecol. Prog. Ser.* **2013**, *474*, 43–52. [CrossRef]

93. Shama, L.N.S.; Wegner, K.M. Grandparental Effects in Marine Sticklebacks: Transgenerational Plasticity across Multiple Generations. *J. Evol. Biol.* **2014**, *27*, 2297–2307. [CrossRef] [PubMed]

94. Fuxjäger, L.; Wanzenböck, S.; Ringler, E.; Wegner, K.M.; Ahnelt, H.; Shama, L.N.S. Within-Generation and Transgenerational Plasticity of Mate Choice in Oceanic Stickleback under Climate Change. *Phil. Trans. R. Soc. B* **2019**, *374*, 20180183. [CrossRef]

95. Farmer, T.M.; Marschall, E.A.; Dabrowski, K.; Ludsin, S.A. Short Winters Threaten Temperate Fish Populations. *Nat. Commun.* **2015**, *6*, 7724. [CrossRef]

96. Tanaka, H.; Nakagawa, T.; Yokota, T.; Chimura, M.; Yamashita, Y.; Funamoto, T. Effects of Spawning Temperature on the Reproductive Characteristics of Walleye Pollock *Gadus chalcogrammus*. *Fish. Sci.* **2019**, *85*, 901–911. [CrossRef]

97. Chambers, R.C. Environmental Influences on Egg and Propagule Sizes in Marine Fishes. In *Early Life History and Recruitment in Fish Populations*; Chambers, R.C., Trippel, E.A., Eds.; Chapman & Hall Fish and Fisheries Series; Springer: Dordrecht, The Netherlands, 1997; pp. 63–102, ISBN 978-94-009-1439-1.

98. Sganga, D.E.; Dahlke, F.T.; Sørensen, S.R.; Butts, I.A.E.; Tomkiewicz, J.; Mazurais, D.; Servili, A.; Bertolini, F.; Politis, S.N. CO₂ Induced Seawater Acidification Impacts Survival and Development of European Eel Embryos. *PLoS ONE* **2022**, *17*, e0267228. [CrossRef]

99. Vézina, D.; Guderley, H. Anatomic and Enzymatic Responses of the Three-Spined Stickleback, *Gasterosteus aculeatus* to Thermal Acclimation and Acclimatization. *J. Exp. Zool.* **1991**, *258*, 277–287. [CrossRef]

100. Jakobsson, S.; Borg, B.; Haux, C.; Hyllner, S.J. An 11-Ketotestosterone Induced Kidney-Secreted Protein: The Nest Building Glue from Male Three-Spined Stickleback, *Gasterosteus aculeatus*. *Fish. Physiol. Biochem.* **1999**, *20*, 79–85. [CrossRef]

101. Alavi, S.; Cosson, J. Sperm Motility in Fishes. (II) Effects of Ions and Osmolality: A Review. *Cell Biol. Int.* **2006**, *30*, 1–14. [CrossRef]

102. Magnotti, C.; Figueroa, E.; Farias, J.G.; Merino, O.; Valdebenito, I.; Oliveira, R.P.S.; Cerqueira, V. Sperm Characteristics of Wild and Captive Lebranche Mullet Mugil Liza (Valenciennes, 1836), Subjected to Sperm Activation in Different PH and Salinity Conditions. *Anim. Reprod. Sci.* **2018**, *192*, 164–170. [CrossRef] [PubMed]

103. Roufidou, C.; Schmitz, M.; Mayer, I.; Sebire, M.; Katsiadaki, I.; Shao, Y.T.; Borg, B. Hormonal Changes over the Spawning Cycle in the Female Three-Spined Stickleback, *Gasterosteus aculeatus*. *Gen. Comp. Endocrinol.* **2017**, *257*, 97–105. [CrossRef]

104. Roufidou, C.; Sebire, M.; Katsiadaki, I.; Mustafa, A.; Schmitz, M.; Mayer, I.; Shao, Y.T.; Borg, B. Overripening of Eggs and Changes in Reproductive Hormones in the Threespine Stickleback. *Evol. Ecol. Res.* **2016**, *17*, 317–329.

105. Shao, Y.T.; Roufidou, C.; Chung, P.C.; Borg, B. Changes in Kisspeptin, GnRH, and Gonadotropin mRNA Levels in Male Threespine Stickleback (*Gasterosteus aculeatus*) during Photoperiod-Induced Sexual Maturation. *Evol. Ecol. Res.* **2019**, *20*, 317–329.

106. Shao, Y.T.; Tseng, Y.C.; Chang, C.-H.; Yan, H.Y.; Hwang, P.P.; Borg, B. GnRH mRNA Levels in Male Three-Spined Sticklebacks, *Gasterosteus aculeatus*, under Different Reproductive Conditions. *Comp. Biochem. Physiol. Part A Mol. Integr. Physiol.* **2015**, *180*, 6–17. [CrossRef]

107. Zohar, Y.; Muñoz-Cueto, J.A.; Elizur, A.; Kah, O. Neuroendocrinology of Reproduction in Teleost Fish. *Gen. Comp. Endocrinol.* **2010**, *165*, 438–455. [CrossRef]

108. Bock, S.L.; Chow, M.I.; Forsgren, K.L.; Lema, S.C. Widespread Alterations to Hypothalamic-Pituitary-Gonadal (HPG) Axis Signaling Underlie High Temperature Reproductive Inhibition in the Eurythermal Sheepshead Minnow (*Cyprinodon variegatus*). *Mol. Cell. Endocrinol.* **2021**, *537*, 111447. [CrossRef]

109. Pigliucci, M. Evolution of Phenotypic Plasticity: Where Are We Going Now? *Trends Ecol. Evol.* **2005**, *20*, 481–486. [CrossRef] [PubMed]

110. Scheiner, S.M. Genetics and Evolution of Phenotypic Plasticity. *Annu. Rev. Ecol. Syst.* **1993**, *24*, 35–68. [CrossRef]

111. Souza, M.L.; Duarte, A.A.; Lovato, M.B.; Fagundes, M.; Valladares, F.; Lemos-Filho, J.P. Climatic Factors Shaping Intraspecific Leaf Trait Variation of a Neotropical Tree along a Rainfall Gradient. *PLoS ONE* **2018**, *13*, e0208512. [CrossRef]

112. Bonamour, S.; Chevin, L.-M.; Charmantier, A.; Teplitsky, C. Phenotypic Plasticity in Response to Climate Change: The Importance of Cue Variation. *Phil. Trans. R. Soc. B* **2019**, *374*, 20180178. [CrossRef] [PubMed]

113. Crozier, L.G.; Hendry, A.P.; Lawson, P.W.; Quinn, T.P.; Mantua, N.J.; Battin, J.; Shaw, R.G.; Huey, R.B. Potential Responses to Climate Change in Organisms with Complex Life Histories: Evolution and Plasticity in Pacific Salmon. *Evol. Appl.* **2008**, *1*, 252–270. [CrossRef] [PubMed]

114. Metcalfe, N.B.; Monaghan, P. Compensation for a Bad Start: Grow Now, Pay Later? *Trends Ecol. Evol.* **2001**, *16*, 254–260. [CrossRef]

115. West-Eberhard, M.J. *Developmental Plasticity and Evolution*; OUP USA: New York, NY, USA, 2003; ISBN 978-0-19-512235-0.

116. Fellous, A.; Wegner, K.M.; John, U.; Mark, F.C.; Shama, L.N.S. Windows of Opportunity: Ocean Warming Shapes Temperature-Sensitive Epigenetic Reprogramming and Gene Expression across Gametogenesis and Embryogenesis in Marine Stickleback. *Glob. Chang. Biol.* **2022**, *28*, 54–71. [CrossRef]
117. O'Dea, R.E.; Lagisz, M.; Hendry, A.P.; Nakagawa, S. Developmental Temperature Affects Phenotypic Means and Variability: A Meta-analysis of Fish Data. *Fish. Fish.* **2019**, *20*, 1005–1022. [CrossRef]

118. Rodriguez-Dominguez, A.; Connell, S.D.; Coni, E.O.C.; Sasaki, M.; Booth, D.J.; Nagelkerken, I. Phenotypic Responses in Fish Behaviour Narrow as Climate Ramps Up. *Clim. Chang.* **2022**, *171*, 19. [CrossRef]

119. Gunderson, A.R.; Stillman, J.H. Plasticity in Thermal Tolerance Has Limited Potential to Buffer Ectotherms from Global Warming. *Proc. R. Soc. B Biol. Sci.* **2015**, *282*, 20150401. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.