ABSTRACT
Fish in coastal ecosystems can be exposed to acute variations in CO2 of between 0.2 and 1 kPa CO2 (2000–10,000 µatm). Coping with this environmental challenge will depend on the ability to rapidly compensate for the internal acid–base disturbance caused by sudden exposure to high environmental CO2 (blood and tissue acidosis); however, studies about the speed of acid–base regulatory responses in marine fish are scarce. We observed that upon sudden exposure to ∼1 kPa CO2, European sea bass (Dicentrarchus labrax) completely regulate erythrocyte intracellular pH within ∼40 min, thus restoring haemoglobin–O2 affinity to pre-exposure levels. Moreover, blood pH returned to normal levels within ∼2 h, which is one of the fastest acid–base recoveries documented in any fish. This was achieved via a large upregulation of net acid excretion and accumulation of HCO3− in blood, which increased from ∼4 to ∼22 mmol l−1. While the abundance and intracellular localisation of gill Na+/K+-ATPase (NKA) and Na+/H+ exchanger 3 (NHE3) remained unchanged, the apical surface area of acid-excreting gill ionocytes doubled. This constitutes a novel mechanism for rapidly increasing acid excretion during sudden blood acidosis. Rapid acid–base regulation was completely prevented when the same high CO2 exposure occurred in seawater with experimentally reduced HCO3− and pH, probably because reduced environmental pH inhibited gill H+ excretion via NHE3. The rapid and robust acid–base regulatory responses identified will enable European sea bass to maintain physiological performance during large and sudden CO2 fluctuations that naturally occur in coastal environments.
INTRODUCTION
Increased CO2 in aquatic environments, or environmental hypercapnia, causes significant physiological challenges for water-breathing animals including fish. As environmental CO2 increases, there is a corresponding rise in CO2 within the fish's blood, which in turn induces a decrease in blood pH (Brauner et al., 2019; Cameron and Randall, 1972). This condition is referred to as respiratory acidosis and, depending on its magnitude, can disrupt multiple homeostatic processes including gas exchange (Crocker and Cech, 1998; Eddy et al., 1977; Perry and Kinkead, 1989) and cardiovascular function (Lee et al., 2003; Perry and McKendry, 2001; Perry et al., 1999). Globally, CO2 levels in the ocean are increasing as a result of anthropogenic climate change, and are predicted to reach ∼0.1 kPa (0.1% CO2, 1000 µatm) by 2100 under a ‘business as usual’ scenario (Orr et al., 2005; Doney et al., 2011; IPCC, 2014).
The increase in oceanic CO2 levels, known as ocean acidification (OA) (Doney et al., 2009), has renewed interest in acid–base regulatory mechanisms of aquatic organisms. However, coastal and estuarine environments already experience much larger variations in CO2 levels (Sunda and Cai, 2012; Wallace et al., 2014), which are likely to be exacerbated in the future (Cai et al., 2011; Melzner et al., 2013). High variation in PCO2 within coastal environments primarily results from diel changes in metabolic activity (i.e. a switch between net autotrophy during the day and net heterotrophy at night) (Baumann et al., 2015; Saderne et al., 2013), riverine inputs (Salisbury et al., 2008), coastal upwelling events (Feely et al., 2008; Saderne et al., 2013) and increased impacts of anthropogenic activities (Cai et al., 2011; Duarte et al., 2013). Fluctuations in PCO2 can occur rapidly over minutes to hours, and reach levels greater than 10 times higher than atmospheric CO2 concentrations (>0.4 kPa CO2) (Baumann et al., 2015; Hofmann et al., 2011; Melzner et al., 2013). In addition, estuarine areas in European coastal waters have been documented to reach PCO2 levels approaching 1 kPa (1% CO2, 10,000 µatm) (Borges et al., 2006), which may result in mobile species experiencing rapid changes in CO2 as they move through these habitats. This type of environmental hypercapnia implies a different physiological challenge to fish compared with OA. Firstly, as environmental CO2 levels increase above CO2 levels in venous blood of fish (typically ∼0.3 kPa), CO2 diffusion gradients are reversed, resulting in net uptake from the environment into the blood and inducing a much more pronounced respiratory acidosis (Tresguerres and Hamilton, 2017). Secondly, the sudden and extreme CO2 fluctuations must be met by equally fast, robust and reversible acid–base regulatory responses.
Fish have a great capacity to restore blood pH to compensate for CO2-induced respiratory acidosis, which is largely achieved by excreting excess H+ and absorbing HCO3− using their gills (Brauner et al., 2019; Claiborne et al., 2002; Esbaugh, 2017; Evans et al., 2005; Perry and Gilmour, 2006). At the cellular level, these processes take place in specialised ion-transporting cells, or ionocytes. However, the underlying ion-transporting proteins and regulatory mechanisms are intrinsically different between freshwater and marine fishes and may also vary between species (Brauner and Baker, 2009; Claiborne et al., 2002; Evans et al., 2005; Perry and Gilmour, 2006). The few studies that have investigated acid–base regulation after acute exposure to ∼1 kPa CO2 in marine fish have reported a large variation of responses, with full blood pH compensation occurring between ∼2 and 24 h depending on the species (Hayashi et al., 2004; Larsen et al., 1997; Perry, 1982; Toews et al., 1983). Given the exquisite sensitivity of most proteins to changes in pH, variation in the time course of acid–base regulatory responses between species has important implications for whole-organism performance. Haemoglobins (Hb) of fish species show strong Bohr and Root effects (Wells, 2009) which reduce Hb–O2 binding affinity and the O2 carrying capacity when erythrocyte intracellular pH (pHi) decreases. While fish have adaptations to regulate pHi of erythrocytes (Cossins and Richardson, 1985; Nikinmaa and Tufts, 1989; Thomas and Perry, 1992), erythrocyte pHi in many fish species (particularly marine fish) is closely linked to whole-blood pH (Brauner and Baker, 2009; Shartau et al., 2020). Therefore, adaptations which enhance the speed of whole-blood acid–base regulation will also provide faster restoration of O2 transport capacity and minimise disruption to energetically expensive activities such as foraging and digestion. However, little is known about why some species are able to compensate for respiratory acidosis faster than others.
The gill ionocytes of marine fish excrete H+ using apical Na+/H+ exchangers (NHEs) driven by basolateral Na+/K+-ATPase (NKA) (Brauner et al., 2019; Claiborne et al., 2002; Evans et al., 2005). Theoretically, H+ excretion during sudden exposure to hypercapnia could be upregulated by increased biosynthesis of NKA and NHE; however, transcriptional and translational responses typically take at least a few hours (e.g. Tresguerres et al., 2005, 2006). Furthermore, protein turnover is energetically expensive (Pan et al., 2015), so short-term regulation of H+ excretion by synthesis and degradation of ion-transporting proteins would not be particularly efficient. Alternatively, the rapid upregulation of H+ excretion may be mediated by post-translational regulatory modifications such as insertion of pre-existing proteins into the ionocyte membrane, or morphological adjustments of its apical area, as reported for a variety of fishes in response to other acid–base disturbances (reviewed in Brauner et al., 2019; Tresguerres et al., 2019).
In the present study, we investigated acid–base regulation of European sea bass, Dicentrarchus labrax, an active predator which seasonally inhabits shallow coastal, estuarine and saltmarsh environments (Doyle et al., 2017) where rapid and large fluctuations in CO2 levels occur (Baumann et al., 2015). Specifically, we characterised blood acid–base regulation, erythrocyte pHi and O2 affinity, effects of seawater chemistry on speed of acid–base regulation, and changes in ionocyte NKA and NHE3 abundance, intracellular localisation, and apical surface area during acute exposure to ∼1 kPa CO2.
MATERIALS AND METHODS
Capture and pre-experimental conditions
Juvenile European sea bass, Dicentrarchus labrax (Linnaeus 1758), were obtained by seine netting in estuaries and salt marshes from Dorset and the Isle of Wight on the south coast of the United Kingdom. Sea bass were transferred to the Aquatic Resources Centre of the University of Exeter where they were held in ∼500 l tanks in a recirculating aquaculture system (RAS, total volume ∼2500 l) at temperatures between 14 and 22°C. Sea bass were fed 3 times a week with commercial pellet (Horizon 80, Skretting) with a supplement of chopped frozen mussel (Mytilus edulis) once a week. For ∼6 months before experiments, sea bass were maintained at a temperature of 14°C and seawater CO2 levels of ∼0.05–0.06 kPa (pH ∼8.10). Prior to all experimental procedures, food was withheld from sea bass for a minimum of 72 h. Animal collections were conducted under appropriate permits (Marine Management Organisation permit no. 030/17 and Natural England permit no. OLD1002654) and all experimental procedures were carried out under UK Home Office licence P88687E07 and approved by the University of Exeter's Animal Welfare and Ethical Review Board.
Hypercapnia exposure
Individual sea bass were moved to isolation tanks (∼20 l) and left to acclimate overnight for a minimum of 12 h before exposure to hypercapnia. During the acclimation period, isolation tanks were fed by the RAS at a rate of ∼4 l min−1; with overflowing water recirculated back to the RAS. After overnight acclimation, hypercapnia exposure began by switching the inflow of water from low CO2 control conditions to high CO2 seawater delivered from a header tank (∼150 l) in which PCO2 levels had already been increased to ∼1 kPa using an Aqua Medic pH computer (AB Aqua Medic GmbH). The pH computer maintained header tank PCO2 levels using an electronic solenoid valve which fed CO2 to a diffuser in the header tank if the pH rose above 6.92 and stopped CO2 flow if the pH dropped below 6.88. Additionally, to reduce CO2 fluctuation in isolation tanks during exposures, the gas aerating each tank was switched from ambient air to a gas mix of 1% CO2, 21% O2 and 78% N2 (G400 Gas mixing system, Qubit Biology Inc.). During exposures, overflowing water from each isolation tank recirculated to the header tank creating an isolated experimental system of ∼250 l. The experimental system was maintained at 14°C using a heater/chiller unit (Grant TX150 R2, Grant Instruments) attached to a temperature exchange coil in the header tank. To characterise the time course of acid–base regulation, sea bass were exposed to hypercapnia (∼1 kPa CO2) for ∼10, ∼40 or ∼135 min before measurements were taken. The pH of isolation tanks was monitored with a separate pH probe and matched the header tank ∼5 min after initial exposure. Measurements of an additional group of sea bass were obtained at normocapnic CO2 levels (∼0.05 kPa CO2) to act as a pre-exposure control (hereafter this group is referred to as time=0). At the time of sampling, measurements of seawater pH (NBS scale), temperature and salinity, as well as samples of seawater to measure total CO2 (TCO2)/dissolved inorganic carbon (DIC) were taken from each isolation tank. DIC analysis was conducted using a custom-built system described in detail by Lewis et al. (2013). Measurements of pH, salinity, temperature and DIC were then input into CO2SYS (Pierrot et al., 2006) to calculate PCO2 and total alkalinity (TA) based on the equilibration constants refitted by Dickson and Millero (1987).
Blood sampling and analysis
Following hypercapnia exposures (Table 1), sea bass were individually anaesthetised in situ with a dose of 100 mg l−1 benzocaine. Anaesthetic was introduced into the isolation box without any visual or physical disturbance to sea bass, thus eliminating capture or handling stress prior to sampling. When sea bass were moderately anaesthetised (no response to physical stimuli and cessation of gill ventilation), they were immediately transferred to a gill irrigation tank (containing 40 mg l−1 benzocaine to maintain anaesthesia) in which gill ventilation was artificially maintained by a micro-pump. Once gill water flow was stable (gill operculum just open and exhalant water flow just visualised), blood was sampled via caudal vessel puncture using a heparinised 1 ml syringe. The gill irrigation tank used was filled with water from the header tank and maintained at an appropriate PCO2 level by aeration with the same gas mix feeding the isolation tanks. The water chemistry of the gill irrigation chamber was measured following the same procedures outlined for the isolation chambers, with one DIC sample taken at the end of blood sampling (Table S1).
Immediately after sampling, extracellular pH was measured on 30 µl of whole blood using an Accumet micro-pH electrode and Hanna HI8314 pH meter calibrated to 14°C using pHNBS 7.04 and 9.21 appropriate buffers. Measurements of blood pH were made in a temperature-controlled water bath. Three 75 µl microcapillary tubes were then filled with whole blood and anaerobically sealed with Critoseal capillary tube sealant (Fisher) and paraffin oil, and centrifuged for 2 min at 10,000 rpm. Haematocrit (Hct) was measured using a Hawksley micro-haematocrit reader. Plasma was then extracted from capillary tubes for analysis of TCO2 using a Mettler Toledo 965D carbon dioxide analyser. Plasma PCO2 and HCO3− were then calculated from TCO2, temperature and blood pH using the Henderson–Hasselbalch equation with values for solubility and pK1app based on Boutilier et al. (1984, 1985). The Hb content of 10 µl of whole blood was also assessed by the cyanmethaemoglobin method (after addition to 2.5 ml of Drabkin's reagent, Sigma). Half the remaining whole blood was centrifuged at 10,000 rpm for 2 min at 4°C. The resulting plasma was separated and 10 µl was diluted in ultrapure water, snap frozen in liquid N2, and stored at −80°C before later being used to measure plasma cation and anion concentrations using ion chromatography (Dionex ICS 1000 and 1100, Thermo Scientific, UK). The remaining plasma was snap frozen in liquid N2 and stored at −80°C before measurements of plasma lactate and glucose were made (YSI 2900D Biochemistry Analyzer, Xylem Analytics). After separating the plasma, the surface of the leftover erythrocyte pellet was blotted to remove the leukocyte layer. The erythrocyte pellet was then snap frozen in liquid nitrogen for 10 s and thawed in a 37°C water bath for 1 min prior to pHi measurements as described by Zeidler and Kim (1977), and validated by Baker et al. (2009b). All measurements or storage of blood occurred within 10 min of blood sampling. Finally, the P50 of whole blood was measured as a proxy for Hb–O2 affinity following the methods outlined in Montgomery et al. (2019) using a Blood Oxygen Binding System (BOBS, Loligo Systems), as detailed by Oellermann et al. (2014). The P50 point is the PO2 at which 50% of Hb is bound to O2 and a higher P50 indicates a right shift in the O2 dissociation curve of whole blood and a reduction in Hb–O2 affinity.
Flux measurements
The flux of acid–base relevant ions between sea bass and seawater was measured over a ∼135 min time period in normocapnic conditions (n=7) and immediately following exposure to hypercapnia (n=8, Table 2). At the start of the measurements the flow to the isolation tanks was stopped and water chemistry maintained at the desired PCO2 by gassing the tanks with either ambient air (control) or a 1% CO2 gas mix (hypercapnia). Seawater samples for measuring TA were taken at the beginning and end of the ∼135 min flux period, preserved by adding 40 µl of 4% (w/v) mercuric chloride per 10 ml of seawater, and stored at 4°C (Dickson et al., 2007) prior to analysis by double titration using an autotitrator (Metrohm 907 Titrando with 815 Robotic USB Sample Processor XL, Metrohm). TA measurements were made using a double titration method modified from Cooper et al. (2010) as detailed by Middlemiss et al. (2016). Briefly, 20 ml samples were titrated to pH 3.89 using 0.02 mol l−1 HCl whilst gassing with CO2-free N2; pH was then returned to starting values by titrating with 0.02 mol l−1 NaOH. Samples for measuring total ammonia were frozen at −20°C before ammonia concentration was measured using a modified version of the colorimetric method of Verdouw et al. (1978) at 660 nm using a microplate reader (NanoQuant infinite M200 pro, Tecan Life Sciences). A calibration curve was constructed using NH4Cl standards in seawater.
Exposure to low total alkalinity and pH
A group of sea bass was exposed to ∼1 kPa CO2 for ∼135 min in low total alkalinity and pH seawater in order to examine whether the mechanism for acid–base regulation was affected by these water chemistry parameters. Sufficient 1 mol l−1 HCl was added to ∼250 l seawater to reduce TA by over 90% from ∼2800 µmol l−1 to ∼200 µmol l−1, followed by overnight aeration to equilibrate CO2 with atmospheric levels. We then adjusted the PCO2 of the low TA seawater to the desired level of ∼1 kPa as described above, and a pH set point of 5.75. Sea bass were placed in the individual isolation boxes (fed by the RAS as detailed for normal TA hypercapnia exposures) and left to acclimate overnight before being exposed to the combined low TA and hypercapnia treatment. Flow to individual isolation boxes was stopped, and ∼75% of the water from the isolation box was drained and refilled with low TA, hypercapnic water. This process was repeated 3 times over a period of ∼5 min. The gas mix aerating each isolation box was switched from ambient air to a 1% CO2 gas mix to maintain the desired PCO2 levels. After ∼135 min exposure, each sea bass was anaesthetised and sampled for blood acid–base measurements as detailed previously. The water chemistry of isolation boxes and gill irrigation chambers was measured at the time of blood sampling (Table S2).
Gill sampling
Gill tissue was sampled from sea bass exposed to ambient CO2 conditions (n=5) and to hypercapnia for ∼135 min (n=5, taken immediately after the flux measurements) in normal TA seawater (Table S3). Mean water chemistry conditions during flux measurements (Table 2) and those experienced by sea bass prior to gill sampling (Table S3) differ because gill samples were only collected from 5 of the 8 sea bass used for flux measurements. Sea bass were euthanised in an anaesthetic bath (benzocaine, 1 g l−1), and gills were dissected and rinsed in phosphate-buffered saline (PBS). The first gill arch on the left side was flash frozen in liquid N2 and stored at −80°C for Western blotting, and the first gill arch on the right side was fixed in 4% paraformaldehyde in 0.1 mol l−1 phosphate buffer saline (PBS) (diluted from 16% electron microscopy grade paraformaldehyde; Electron Microscope Science catalogue no. 15710), overnight at 4°C for immunohistochemistry. Following ∼10 h fixation, gill samples were transferred to 50% ethanol for ∼10 h at 4°C, and then stored in 70% ethanol at 4°C.
Antibodies
NKA was immunodetected using α5, a mouse monoclonal antibody against the α-subunit of chicken NKA (a5, Developmental Studies Hybridoma Bank, Iowa City, IA, USA; Lebovitz et al., 1989). This antibody universally recognises NKA in teleost fishes including yellowfin tuna (Thunnus albacares; Kwan et al., 2019), Pacific chub mackerel (Scomber japonicus; Kwan et al., 2020) and California killifish (Fundulus parvipinnis; Nadler et al., 2021). Rabbit anti-NHE3 polyclonal antibodies were generously donated by Dr Junya Hiroi (St Marianna University School of Medicine, Kawaski, Japan); they target two epitope regions within rainbow trout (Oncorhynchus mykiss) NHE3b (GDEDFEFSEGDSASG and PSQRAQLRLPWTPSNLRRLAPL), and recognise NHE3 of multiple teleost species including European sea bass (D. labrax; Blondeau-Bidet et al., 2019). The secondary antibodies were goat anti-mouse HRP-linked and goat anti-rabbit HRP-linked (Bio-Rad, Hercules, CA, USA) for immunoblotting, and goat anti-mouse Alexa Fluor 546 and goat anti-rabbit Alexa Fluor 488 (Invitrogen, Grand Island, USA) for immunohistochemistry.
Western Blotting
Western blotting followed the procedures outlined in Kwan et al. (2019, 2020). While frozen on dry ice, the gill filament and lamellae were separated from the gill arch using a razor blade. The excised tissue was then immersed in liquid N2 and pulverised in a porcelain grinder, then submerged within an ice-cold, protease-inhibiting buffer (250 mmol l−1 sucrose, 1 mmol l−1 EDTA, 30 mmol l−1 Tris, 10 mmol l−1 benzamidine hydrochloride hydrate, 1 mmol l−1 phenylmethanesulfonyl fluoride, 1 mmol l−1 dithiothreitol, pH 7.5). Samples were further homogenised using a handheld VWR Pellet Mixer (VWR, Radnor, PA, USA) for 15 s intervals (3 times) while on ice. Next, samples were centrifuged (3000 g, 4°C; 10 min), and the resulting supernatant was considered the crude homogenate. An aliquot of the crude homogenate was further subjected to a higher speed centrifugation (21,130 g, 4°C; 30 min), and the pellet was saved as the membrane-enriched fraction. The Bradford assay was used to determine protein concentration (Bradford, 1976), which was used to normalise protein loading.
On the day of Western blotting, samples were mixed with an equal volume of 90% 2× Laemmli buffer and 10% β-mercaptoethanol. Samples were then heated at 70°C for 5 min, and the proteins (20 µg per lane) were loaded onto a polyacrylamide mini gel (4% stacking, 10% separating) – alternating between control and high CO2 treatments to avoid possible gel lane effects. The gel ran at 60 V for 15 min, then 100 V for 1.5 h, and proteins were transferred onto a polyvinylidene difluoride (PVDF) membrane using a wet transfer cell (Bio-Rad) at 70 V for 2 h at 4°C. PVDF membranes were then incubated in Tris-buffered saline with 1% Tween (TBS-T) with milk powder (0.1 g ml−1) at room temperature for 1 h, then incubated with primary antibody (NKA: 10.5 ng ml−1; NHE3: 1:1000) in blocking buffer at 4°C overnight. The following day, PVDF membranes were washed in TBS-T (3 times; 10 min each), incubated in blocking buffer with secondary antibodies (1:10,000) at room temperature for 1 h, and washed again in TBS-T (3 times; 10 min each). Bands were made visible through addition of ECL Prime Western Blotting Detection Reagent (GE Healthcare, Waukesha, WI, USA) and imaged with the Universal III Hood (Bio-Rad). Following imaging, the PVDF membrane was incubated in Ponceau stain (10 min, room temperature) to estimate protein loading. Relative NKA and NHE protein abundance (n=5 per treatment) was quantified using the Image Lab software (version 6.0.1; Bio-Rad) and normalised by the protein content in each lane. The anti-NKA antibodies detected a single band of ∼100 kDa (Fig. S2, left), which matches previous reports from multiple fish species ranging from ∼90 to ∼120 kDa (Chandrasekar et al., 2014; Dutta et al., 2019; Frommel et al., 2021; Kwan et al., 2020, 2021; McCormick et al., 2009). The anti-NHE3 antibodies detected a single band of ∼70 kDa (Fig. S2, right), which is smaller than what may be expected based on bioinformatic predictions for NHE3s from other animals. The discrepancy between predicted and apparent NKA and NHE3 sizes may be explained by intrinsic limitations in Western blotting for precisely determining protein size (Bass et al., 2017; Gwozdz and Dorey, 2017), together with the presence of post-translational modifications and isoforms (see Harter et al., 2021, for a case of a fish NHE with an alternative start codon that results in a smaller-than-predicted protein). For example, the predicted size of NHE3 is 97 kDa in climbing perch Anabas testudineus, 94 kDa in killifish Fundulus heteroclitus and 93 kDa in rabbit Oryctolagus cuniculus, yet their apparent Western blot sizes are ∼85 kDa (Chen et al., 2017), ∼77 kDa (Edwards et al., 2005) and ∼80 kDa (Biemesderfer et al., 1997), respectively. Although many studies have previously used the same anti-NHE3 antibodies as us for immunohistochemistry in other fish species (Blondeau-Bidet et al., 2019; Christensen et al., 2012; Hiroi and McCormick, 2012; Seo et al., 2013; Wright et al., 2016), to our knowledge only Wright et al. (2016) performed Western blots and did not report the predicted size of the immunoreactive band. Regardless, the presence of single Western blot bands close to previously reported protein sizes combined with the labelling of expected subcellular localisations of specific gill cells (see below) indicates that the anti-NKA and anti-NHE3 antibodies used in the current study specifically recognise their target proteins.
Immunohistochemistry
Immunohistochemistry was performed as described in Kwan et al., (2020). Fixed gill tissue stored in 70% ethanol was rehydrated in PBS+0.1% Tween (PBS-T) for 10 min, and gill filaments were dissected out to ease subsequent imaging. Autofluorescence was quenched by rinsing in ice-cold PBS with sodium borohydride (1.0–1.5 mg ml−1; 6 times; 10 min each), followed by incubation in blocking buffer (PBS-T, 0.02% normal goat serum, 0.0002% keyhole limpet haemocyanin) at room temperature for 1 h. Samples were incubated with blocking buffer containing primary antibodies (NKA: 40 ng ml−1; NHE3: 1:500, cf. Seo et al., 2013) at 4°C overnight. The following day, samples were washed in PBS-T (3 times at room temperature; 10 min each), and incubated with the fluorescent secondary antibodies (1:500) counterstained with DAPI (1 µg ml−1) at room temperature for 1 h. Samples were washed again in PBS-T as before, then placed on a concave slide for imaging using an inverted confocal microscope (Zeiss LSM 800 with Zeiss ZEN 2.6 blue edition software; Cambridge, UK). Samples incubated without primary antibodies had no signal (Fig. S1A).
Quantification of ionocyte apical surface area
The apical surface area of gill ionocytes were quantified through a combination of whole-mount imaging (40× objective lens with deionised water immersion), optical sectioning, and XZ- and YZ-plane analysis. Gill filaments immunostained for NKA and NHE3 were coded, then blind analysis was employed to eliminate experimental bias. A relatively flat surface on the gill filament was selected under 0.5× scanning confocal magnification to ensure imaging would be performed on ionocytes in an upright position, thus minimising errors in apical surface area quantification due to angle distortion. After locating an ionocyte by its distinctive NKA signal, the scanning confocal magnification was increased to 5.0× and the entire cell was Z-stack imaged (optimal interval automatically selected: 0.07–0.12 µm per slice). Subsequent viewing of the Z-stack from the X–Z and Y–Z planes allowed us to assess intracellular localisation, and to identify the image slice that captured the entire apical surface (typically, the second slice from the top of the cell). Next, the ionocyte apical surface area (identified by NHE3 immunofluorescence signal) was quantified using FIJI (Schindelin et al., 2012). The apical surface area for control and hypercapnia-exposed fish (n=5) was determined as the mean apical surface area of five ionocytes taken from different gill filaments.
Statistical analysis
All statistical analysis was performed using R v3.6.3 (http://www.R-project.org/). Changes in blood chemistry parameters over time in response to hypercapnia exposure were analysed using one-way ANOVA before assumptions of equal variance of data and normality of model residuals were checked. Post hoc tests were conducted on least-square means generated by package ‘emmeans’ (https://CRAN.R-project.org/package=emmeans), with Tukey adjusted P-values for multiple comparisons. Some data did not meet required assumptions for one-way ANOVA. Unequal variances were observed in measurements of plasma PCO2 between treatments; as such, we used Welch's ANOVA with Tukey's pairwise comparisons using Benjamini–Hochberg corrections for post hoc testing. Measurements of blood pH and P50 did not meet assumptions of normality and were analysed using the Kruskal–Wallis test with post hoc comparisons made with Dunn's test from package ‘FSA’ (https://github.com/droglenc/FSA), using Benjamini–Hochberg corrections for multiple comparisons. As a result of unusually high measurements of plasma [Cl−] and [Na+] in some samples, a ROUT test was conducted (Q=0.5%) in Graphpad Prism 9 to identify potential outliers (Motulsky and Brown, 2006). Samples in which plasma [Cl−] and [Na+] were both identified as outliers (n=6) by the ROUT test were excluded from the dataset prior to subsequent statistical analysis. The relationship between plasma [Cl−] and [HCO3−] as well as plasma [Na+] and [H+] across treatments was examined using Kendall's tau correlation. Comparisons of blood chemistry between fish exposed to hypercapnia for ∼135 min in normal alkalinity and low total alkalinity (and pH) seawater were made using two-sample t-tests. Measurements of plasma HCO3− were analysed using Welch's t-test as unequal variances were observed between fish exposed to hypercapnia in normal and low total alkalinity seawater. Flux measurements were analysed using Student's t-test after checking data met assumptions of normality and equal variance. Relative protein abundance and ionocyte apical area met both assumptions of normality and equal variance and were analysed using one-tailed t-test (control response<CO2-exposed response).
RESULTS
Blood chemistry
Exposure to environmental hypercapnia caused significant changes in blood pH of sea bass over time (Kruskal–Wallis test, χ2=25.0, d.f.=3, P<0.001). There was a pronounced acidosis of the blood from pH 7.84±0.02 (mean±s.e.m.) in control conditions (normocapnia, time=0) to 7.50±0.03 after exposure to hypercapnia for ∼10 min (Fig. 1A,D). Following this initial acidosis, sea bass completely restored blood pH to control levels after ∼135 min (Fig. 1A,D). Blood pH regulation was accompanied by a ∼5-fold increase in plasma HCO3−, from 4.5±0.3 to 21.9±0.7 mmol l−1, over the ∼135 min exposure (Fig. 1C,D; one-way ANOVA, F=203.3, d.f.=3, P<0.001).
Plasma PCO2 showed significant changes during exposure to hypercapnia (Fig. 1B,D; Welch's ANOVA, F=202.5, d.f.=3, P<0.001). The initial decrease in blood pH of sea bass was driven by a rapid and large (∼6-fold) increase in plasma PCO2, from 0.200±0.016 to 1.185±0.049 kPa CO2, within the first 10 min of exposure. There was a small but significant decline in plasma PCO2 between sea bass sampled ∼10 min after exposure and sea bass sampled ∼135 min after exposure (Fig. 1B). There were no significant differences in plasma glucose or lactate levels between any treatment groups, with values for all sea bass of 5.90±0.43 mmol l−1 and 0.45±0.05 mmol l−1, respectively.
To test the influence of total alkalinity and environmental pH on acid–base regulation, a group of sea bass were exposed to hypercapnia for ∼135 min in low alkalinity and pH seawater. These sea bass were unable to compensate for a respiratory acidosis when exposed to acute hypercapnia for ∼135 min. Blood pH was 0.37 units (95% confidence interval, CI 0.33–0.41) lower than that of sea bass exposed to hypercapnia in normal alkalinity seawater for the same length of time (Fig. 1A; two-sample t-test, d.f.=15, t=−13.551, P<0.001). Additionally, sea bass in low alkalinity seawater did not actively accumulate HCO3− when exposed to environmental hypercapnia for ∼135 min (Fig. 1C; Welch's t-test, d.f.=9.8996, t=−16.817, P<0.001). Indeed, the increase in mean plasma [HCO3−] followed the predicted non-bicarbonate buffering line (Fig. 1D).
Flux measurements
Sea bass switched from slight net base excretion under control normocapnic conditions to net acid excretion that was ∼2.5-fold larger in magnitude during 135 min of hypercapnia (Fig. 2C, Student's t-test, t=−2.25, d.f.=13, P=0.042). This was driven by a switch from a small apparent HCO3− excretion to a large apparent HCO3− uptake (Fig. 2A). There were no significant differences in Tamm excretion (Fig. 2B).
Oxygen transport capacity
The initial drop in blood pH during hypercapnia exposure was reflected in changes in erythrocyte pHi (Fig. 3A; one-way ANOVA, F=12.34, d.f.=3, P<0.001). However, erythrocyte pHi returned to control levels after ∼40 min of exposure to hypercapnia (Fig. 3A). As expected, the significant changes in erythrocyte pHi affected Hb–O2 binding affinity, leading to a ∼3-fold increase in P50 after 10 min, from 1.78±0.30 kPa O2 in control sea bass to 5.60±0.36 kPa O2 (Kruskall–Wallis test, χ2=17.4, d.f.=3, P<0.001). The rapid recovery of erythrocyte pHi after ∼40 min led to P50 returning to pre-exposure levels (Fig. 3B).
Sea bass exposed to hypercapnia also experienced a ∼25% increase in Hb levels (Fig. 3D), at ∼10 min and∼135 min compared with control sea bass (one-way ANOVA, F=4.60, d.f.=3, P=0.009). In addition, sea bass exposed to hypercapnia exhibited an ∼8–10% increase in Hct (Fig. 3C), although this increase was marginally non-significant (one-way ANOVA, F=2.40, d.f.=3, 0.086).
Plasma ion concentrations
Plasma [Cl−] significantly decreased by 13.1 mmol l−1 (95% CI 10.1–16.2 mmol l−1) from 134.9±5.1 mmol l−1 (mean±s.e.m.) in sea bass in normocapnia to 121.7±1.3 mmol l−1 in sea bass exposed to hypercapnia for ∼135 min (Kruskall–Wallis test, χ2=7.87, d.f.=3, P=0.049). Sea bass exposed to hypercapnia in low alkalinity seawater had significantly higher plasma [Cl−] when compared with fish in normal alkalinity seawater (Fig. 4A; two sample t-test, d.f.=11, t=5.7306, P<0.001). Decreases in plasma [Cl−] showed a correlation with increasing bicarbonate (Fig. 4A inset; Kendall's tau correlation, τ=−0.32, P=0.005). Plasma [Na+] showed no significant changes over the time course of hypercapnia exposure (one-way ANOVA, F=1.063, P=0.391), and there were no differences in [Na+] after ∼135 min of exposure to hypercapnia between sea bass in normal and low TA seawater (Fig. 4B; two-sample t-test, d.f.=11, t=−0.84857, P=0.414). As such, there was no correlation between plasma [Na+] and [H+] (Fig. 4B inset; Kendall's tau correlation, τ=0.083, P=0.47).
NKA and NHE3 protein abundance
Exposure to hypercapnia did not induce significant changes in the abundance of NKA or NHE3 in crude homogenates (indicative of total protein abundance) or the abundance of NKA and NHE3 in membrane-enriched fractions (indicative of protein that was present in the apical or basolateral plasma membranes) (one-tailed t-test, P>0.05; Table 3).
Ionocyte intracellular localisation and apical surface area
NKA-rich ionocytes were primarily localised on the gill filament trailing edges and the basal portion of the gill lamellae (Fig. S1B); all NKA-rich ionocytes also expressed NHE3 in their apical region (Fig. 5A). Despite analysis using high magnification imaging, optical sectioning and X–Z and Y–Z plane visualisation, we found no evidence of NHE3 intracellular localisation (Fig. 5B,Bi,C,Ci). Blind analysis revealed the ionocyte's apical surface area (based on the NHE3 signal) significantly increased, almost doubling from 3.34±0.17 to 6.36±0.49 µm2, after exposure to ∼135 min of hypercapnia (one-tailed t-test, t=5.873, d.f.=4.963, P=0.001; Fig. 5D).
DISCUSSION
Our results indicate that European sea bass are able to rapidly compensate for hypercapnia-induced blood acidosis when the environmental CO2 is at the extreme high end of the spectrum encountered in their natural habitat. Complete restoration of blood pH after exposure to ∼1 kPa CO2 was achieved within ∼2 h via a switch from net base excretion to net acid excretion and a subsequent accumulation of HCO3− in plasma. In addition, erythrocyte pHi and Hb–O2 affinity were restored to pre-exposure levels after just ∼40 min, and there was a 20% increase in the blood Hb concentration together with a trend for ∼10% Hct increase. These results suggest an adrenergic response that stimulates a β-NHE in the erythrocytes (Nikinmaa, 2012), and contracts the spleen, resulting in the release of erythrocytes into the circulation (Crocker and Cech, 1997; Lee et al., 2003; Perry and Kinkead, 1989; Vermette and Perry, 1988). The end result is a boost in blood O2 transport capacity that counteracts the reduced Hb–O2 affinity induced by the initial hypercapnia-induced acidosis.
Regulation of respiratory acidosis by sea bass exposed to hypercapnia in normal alkalinity sea water (TA ∼2800 µmol l−1) resulted in an elevation of plasma [HCO3−] by ∼18 mmol l−1, which was correlated with a decrease in plasma [Cl−] of ∼13 mmol l−1. In comparison, while we saw a slight rise in plasma [Na+] on initial exposure to hypercapnia, this was transient, and there was no overall correlation between plasma [Na+] and [H+] during the whole 135 min experiment. However, a lack of increase in plasma [Na+] during acid–base regulation does not preclude increased Na+ uptake (to facilitate H+ excretion) during acid–base regulation. Instead, a lack of increased plasma [Na+] may simply reflect upregulation of the hypo-ionoregulatory mechanism for NaCl excretion in marine teleosts, which depends on transcellular Cl− movement via basolateral Na+/K+/2Cl− cotransporters and apical Cl− channels, and on paracellular Na+ movement (Zadunaisky, 1996). In principle, this could compensate for the enhanced uptake of Na+ into the blood that is associated with H+ excretion by NHE, and would also help explain the observed reduction in plasma [Cl−] in fish exposed to hypercapnia.
In comparison, sea bass exposed to hypercapnia in low alkalinity sea water (TA ∼200 µmol l−1) showed no ability to accumulate HCO3−, to compensate for respiratory acidosis, and did not experience a decrease in plasma [Cl−]. At face value, these results may support potential direct uptake of HCO3− from seawater in exchange for blood Cl− through HCO3−/Cl− exchange across the gills (Esbaugh et al., 2012; Perry and Gilmour, 2006; Tovey and Brauner, 2018). However, the thermodynamics of this proposed mechanism are not clear, as [Cl−] is much higher in seawater than in internal fluids of fish, and the opposite is true for [HCO3−]. This implies that both counterions would have to be transported against their concentration gradients and, furthermore, these gradients would become increasingly unfavourable as blood acidosis is compensated for. Importantly, our hypercapnic low alkalinity seawater had a pH of ∼5.7, which was ∼1.2 pH units lower than hypercapnic normal alkalinity seawater (a 15-fold increase in [H+]). Based on nominal values of [Na+] and [H+] inside fish gill ionocytes and calculations in Parks et al. (2008), the low alkalinity seawater would not sustain H+ excretion via apical NHEs (Table S4). Interestingly, if the NHE3 is not activated, there would not be an associated increased Na+ load into the blood and, as Na+ excretion is coupled to Cl− excretion as explained above, there would be no need to upregulate NaCl excretion, potentially explaining the observed lack of decrease of plasma [Cl−]. Overall, we believe this evidence supports enhanced NHE3-mediated H+ excretion and associated retention of metabolically produced HCO3− in the blood as the primary mechanism underlying regulation of respiratory acidosis in sea bass. However, we cannot exclude a potential role for Cl−-dependent HCO3− absorption from sea water as an additional mechanism for acid–base compensation by sea bass. Evidence from some other teleost species after exposure to equivalent levels of environmental CO2 show changes in gene expression of ionocyte transporters and carbonic anhydrase, which may lead to more favourable conditions for apical HCO3− uptake (Esbaugh et al., 2012; Tseng et al., 2013). Unfortunately, experiments to produce unequivocal evidence to distinguish between apical H+ excretion and HCO3− uptake during acid–base regulation are not trivial because of the interlink and reciprocal compensation between ion-transporting processes for acid–base and ionoregulation in marine fish gills and intestine.
Immunostaining confirmed the presence of basolateral NKA and apical NHE3 in European sea bass gill ionocytes, indicating these cells play a role in H+ excretion. Moreover, all NKA-rich cells co-expressed NHE3. Noteworthy, a recent study reported a lack of NH3 immunostaining in gills of seawater-acclimated D. labrax (Blondeau-Bidet et al., 2019). This is puzzling because they used the same species, similar environmental salinity and the same antibodies as in our study. Additionally, that study was able to detect NHE3 mRNA in gills from the same specimens. We do not know the reasons for this discrepancy between positive gene expression but a lack of immunostaining for the protein itself, but we speculate this might be due to differences between D. labrax populations, feeding and exercise regimes, immunostaining protocols, or reagents, potentially coupled with a disconnect between mRNA and protein expression.
To investigate the mechanisms used by sea bass to enhance acid excretion, we examined whether changes in gill NKA and NHE3 occur after acute (∼135 min) exposure to hypercapnia. Gill NKA and NHE3 protein abundance did not change, ruling out increased protein synthesis as the mechanism responsible for the observed upregulation in acid excretion; this is not surprising considering the short time frame of our experiments. We also examined the potential translocation of pre-existing NKA and NHE3 to the ionocyte basolateral and apical membranes, respectively. Such mechanisms upregulate acid–base regulatory ion transport in elasmobranchs (Roa et al., 2014; Tresguerres et al., 2005, 2006, 2007b) and hagfish (Parks et al., 2007; Tresguerres et al., 2007a); however, NKA and NHE3 protein abundance in the gill membrane fraction of European sea bass gills was also unchanged, ruling out NKA and NHE3 translocation in our experiments. Finally, we hypothesised that sea bass could have remodelled the apical membrane of ionocytes to increase the sites for H+ excretion in contact with seawater. Indeed, this was the case as the surface area of the NHE3-abundant apical membrane of NKA-rich ionocytes roughly doubled in response to hypercapnia. Based on the accepted model of acid excretion in marine fish gill ionocytes (Brauner et al., 2019; Claiborne et al., 2002; Evans et al., 2005), the mechanism probably also includes enhanced carbonic anhydrase activity as well as decreased HCO3− excretion via apical anion exchangers and H+ absorption into the blood via basolateral NHE1. Given the velocity of the response (∼2 h), these changes in activity are likely to be mediated by changes in the availability of CO2 substrate and post-translational modifications rather than changes in protein abundance. Interestingly, marine Japanese flounder (Paralichthys olivaceus) exposed to 1 kPa CO2 for 24 h also showed evidence for an approximately doubling in gill ionocyte apical surface area (although the differences with control fish were not statistically significant, possibly as a result of the large and variable responses, as seen in fish exposed to 5 kPa CO2, which were part of the same analyses even though many of them died between 24 and 72 h after hypercapnic exposure) (Hayashi et al., 2013). Moreover, gill NKA activity remained unchanged upon 24 h exposure to 1 kPa CO2, suggesting that, like European sea bass, Japanese flounder are able to cope with this level of hypercapnia without having to upregulate the biosynthesis of ion-transporting proteins. Additional mechanisms for acid–base regulation, which were not measured in this study, may also contribute to compensation of acute respiratory acidosis in marine fish. For example, hyperventilation by red drum (Sciaenops ocellatus) in response to a 4 h exposure to 0.5 kPa CO2 has been estimated to reduce the need for metabolic compensation by ∼39% as a result of lower blood PCO2 (Ern and Esbaugh, 2016). Although blood sampling in our study involves forced gill ventilation, potentially eliminating the ventilatory responses observed by Ern and Esbaugh (2016), we believe it is unlikely that the rapid acid–base regulation we observed in sea bass is an artefact of this sampling method. Indeed, it has previously been demonstrated that this method provides equivalent blood acid–base measurements to those for cannulated fish in which natural ventilatory responses would be present (see fig. 6 in Brauner et al., 2019).
Previous studies on freshwater fishes have also documented morphological adjustments in gill ionocytes upon comparable hypercapnic exposures. However, the responses were the opposite to our study, as those freshwater fishes experienced a significant reduction in ionocyte apical surface area (Baker et al., 2009a; Goss et al., 1992b; Leino and McCormick, 1984). In some cases, the apical membrane retracted into a more pronounced apical pit (Goss et al., 1992a, 1994), which was suggested to create a microenvironment with higher [Na+] compared with the bulk freshwater and to facilitate Na+/H+ exchange (Kumai and Perry, 2012). In contrast, exposure to more pronounced hypercapnia (8 kPa CO2 over 4 days) induced an increase in gill ionocyte apical surface area in freshwater catfish (Ictalurus punctatus) (Cameron and Iwama, 1987). This response was similar to that in the sea bass in our study; however, the considerably longer exposures and higher PCO2 levels probably resulted in additional acid–base responses that were not investigated, such as increased synthesis of ion-transporting proteins, or a change in the H+ excreting mechanism (e.g. NHE versus V-type H+-ATPase). Regardless, the ability of sea bass to rapidly compensate for a blood respiratory acidosis by simply increasing gill ionocyte apical surface area is in large part possible as a result of the overabundance of Na+ in sea water, which establishes favourable conditions for NHE-mediated H+ excretion.
Freshwater species typically take from 24 h to >72 h to regulate blood pH after exposure to 1 kPa CO2 (Baker et al., 2009a; Claiborne and Heisler, 1984, 1986; Damsgaard et al., 2015; Larsen and Jensen, 1997; Perry, 1982; Perry et al., 1981; Smatresk and Cameron, 1982). While it is generally believed that marine teleosts can regulate their blood acid–base status at a faster rate than freshwater species (Brauner et al., 2019), relatively little research has been conducted to characterise the speed of acid–base regulation in marine fish. A bibliography search revealed four previous studies on five marine teleost species that characterised the time course of the acid–base regulatory response after exposure to 1 kPa CO2 (Fig. 6). Of these species, only the Japanese amberjack (Seriola quinqueradiata) was able to restore blood extracellular pH (pHe) faster than sea bass (∼60 min versus ∼135 min; Fig. 6B). The remaining four species regulated blood pHe between 3 and 24 h post-CO2 exposure (Fig. 6C–F).
The robust ability of sea bass to rapidly acid–base regulate in response hypercapnia probably plays a significant role in their natural environment. Specifically, sea bass feed in shallow coastal estuaries and salt marsh habitats during the summer (Doyle et al., 2017), and these habitats typically experience large fluctuations in CO2 levels over short time periods (Hofmann et al., 2011; Melzner et al., 2013; Wallace et al., 2014). For example, equivalent salt marshes on the US coast regularly experience CO2 fluctuations of ∼0.4 kPa across tide cycles during the summer (Baumann et al., 2015). We should note that natural CO2 levels equivalent to those used in our study (∼1 kPa) have previously only been reported in upper estuaries where salinity is reduced (Borges et al., 2006) and high CO2 levels in estuaries often negatively correlate with salinity (Frankignoulle and Borges, 2001). Future studies should examine the impacts of salinity on acid–base regulation in euryhaline marine fish which may encounter acute hypercapnia in combination with reduced salinity. Despite this, the fast acid–base regulatory response observed in our study indicates that sea bass would be able to rapidly correct the respiratory blood acidosis caused by environmentally relevant, acute CO2 variation in seawater in <1 h. Critically, regulation of erythrocyte pHi by sea bass occurred quicker than regulation of blood pHe. By rapidly restoring erythrocyte pHi, O2 transport capacity is maintained, which is crucial for active predatory teleosts. However, environmental CO2 variation cannot be the sole driver for enhanced acid–base regulatory capacities in all species. For example, Japanese amberjack show a similarly fast blood acid–base regulatory response (Hayashi et al., 2004) but primarily inhabit pelagic, offshore ecosystems in which large variation in environmental CO2 would be less likely to occur. Alternatively, it may be that active, predatory species have developed higher capacities for acid–base regulation to deal with large metabolic acidosis (as a result of anaerobic respiration used during intense exercise involved in prey capture). Understanding the mechanisms that determine species-specific differences in acid–base regulatory capacity will help us to understand the differential impacts of acute exposure to elevated CO2, both by itself and in combination with other stressors such as hypoxia. For example, we have recently reported that sea bass showed enhanced hypoxia tolerance when exposed to progressive and environmentally relevant hypercapnia and hypoxia over a 6 h period (Montgomery et al., 2019). In contrast, European plaice (Pleuronectes platessa) and European flounder (Platichthys flesus) exposed to the same conditions showed reduced hypoxia tolerance (Pcrit), which was associated with reduced Hb–O2 affinity and O2 uptake resulting from an uncompensated respiratory acidosis (Rogers, 2015; D.W.M. and R.W.W., unpublished observations). Thus, species with more robust acid–base regulatory mechanisms may be more resilient to interactive effects between hypercapnia and hypoxia.
Overall, our study highlights the capacity of European sea bass to rapidly (2 h) regulate blood and erythrocyte acid–base status and O2 transport capacity upon exposure to a pronounced and sudden increase in environmental CO2 levels. The ability of sea bass to rapidly upregulate H+ excretion appears to be mediated via the increased exposure of NHE3-containing apical surface area of gill ionocytes, rather than changes in NHE3 or NKA protein abundance or localisation. Additionally, sea bass erythrocyte pHi is regulated even more rapidly than blood pH (40 min), which enables them to quickly restore the affinity of Hb for O2, and therefore blood O2 transport capacity during exposure to elevated CO2. In conjunction, these acid–base regulatory responses will minimise the impact of pronounced and rapidly fluctuating CO2 in their natural environments, and so may prevent disruption of energetically costly activities such as foraging or digestion, and may make sea bass more resilient to the impacts of hypoxia and additional stressors during acute periods of hypercapnia. This is an avenue where we believe further research effort is necessary.
Acknowledgements
We thank Dr Junya Hiroi for providing the NHE3 antibodies, and Dr Andrew Esbaugh and Dr Till Harter for discussions related to data interpretation. We would similarly like to thank the aquarium staff, particularly Dr Gregory Paul, Paul Tyson, Rebecca Turner, Alex Bell and Richard Silcox, of the Aquatic Resource Centre at the University of Exeter for assistance with fish husbandry and maintenance of aquarium facilities.
Footnotes
Author contributions
Conceptualization: D.W.M., R.W.W.; Methodology: D.W.M., G.T.K., R.W.W.; Formal analysis: D.W.M., G.T.K., S.S., G.H.E., S.N.R.B., M.T., R.W.W.; Investigation: D.W.M., G.T.K., W.G.D., J.F., A.B., R.W.W.; Resources: M.T., R.W.W.; Data curation: D.W.M.; Writing - original draft: D.W.M., G.T.K., M.T., R.W.W.; Writing - review & editing: D.W.M., G.T.K., W.G.D., J.F., A.B., S.S., G.H.E., S.N.R.B., M.T., R.W.W.; Visualization: D.W.M., G.T.K.; Supervision: R.W.W.; Project administration: D.W.M.; Funding acquisition: D.W.M., G.T.K., M.T., R.W.W.
Funding
This work was supported by a NERC GW4+ Doctoral Training Partnership studentship from the Natural Environment Research Council (NE/L002434/1) with additional funding from CASE partner, The Centre for Environment, Fisheries and Aquaculture Science (Cefas) to D.W.M., and from the Biotechnology and Biological Sciences Research Council (BB/D005108/1 and BB/J00913X/1) and Natural Environment Research Council (NE/H017402/1) to R.W.W. G.T.K. was funded by the National Science Foundation (NSF) Postdoctoral Fellowship Program (#1907334). M.T. received funding from the National Science Foundation (IOS #1754994). Open Access funding provided by University of Exeter. Deposited in PMC for immediate release.
Data availability
Data are available from the University of Exeter's Open Research Exeter (ORE) depository: https://doi.org/10.24378/exe.3723
References
Competing interests
The authors declare no competing or financial interests.