- 1Centro de Ciências do Mar, Faro, Portugal
- 2Energy and Environment Institute, University of Hull, Hull, United Kingdom
- 3Department of Biological and Marine Science, University of Hull, Hull, United Kingdom
- 4E.A. Milne Centre for Astrophysics and G.W. Gray Centre for Advanced Material, University of Hull, Hull, United Kingdom
The effects of ocean acidification on fish are only partially understood. Studies on olfaction are mostly limited to behavioral alterations of coral reef fish; studies on temperate species and/or with economic importance are scarce. The current study evaluated the effects of short- and medium-term exposure to ocean acidification on the olfactory system of gilthead seabream (Sparus aurata), and attempted to explain observed differences in sensitivity by changes in the protonation state of amino acid odorants. Short-term exposure to elevated PCO2 decreased olfactory sensitivity to some odorants, such as L-serine, L-leucine, L-arginine, L-glutamate, and conspecific intestinal fluid, but not to others, such as L-glutamine and conspecific bile fluid. Seabream were unable to compensate for high PCO2 levels in the medium term; after 4 weeks exposure to high PCO2, the olfactory sensitivity remained lower in elevated PCO2 water. The decrease in olfactory sensitivity in high PCO2 water could be partly attributed to changes in the protonation state of the odorants and/or their receptor(s); we illustrate how protonation due to reduced pH causes changes in the charge distribution of odorant molecules, an essential component for ligand-receptor interaction. However, there are other mechanisms involved. At a histological level, the olfactory epithelium contained higher densities of mucus cells in fish kept in high CO2 water, and a shift in pH of the mucus they produced to more neutral. These differences suggest a physiological response of the olfactory epithelium to lower pH and/or high CO2 levels, but an inability to fully counteract the effects of acidification on olfactory sensitivity. Therefore, the current study provides evidence for a direct, medium term, global effect of ocean acidification on olfactory sensitivity in fish, and possibly other marine organisms, and suggests a partial explanatory mechanism.
Introduction
Since the Industrial Revolution, atmospheric carbon dioxide (CO2) levels have been steadily increasing and the oceans of the world act as a sink for much of this. As atmospheric CO2 has increased from pre-industrial (about 280 μatm CO2) to present day values (about 390 μatm CO2) (Doney et al., 2009), this has led to a fall in oceanic pH of 0.1 unit, with a further fall of 0.77 units predicted by the year 2300 (Doney et al., 2009). Global mean oceanic CO2 levels are expected to reach 1,000 μatm by the year 2100, and 1900 μatm CO2 by 2300 (Caldeira and Wickett, 2003). Studies prior to 2009 suggested that fish, as effective acid-base regulators, were not particularly sensitive to changes in PCO2. However, more recent studies have documented notable impacts on neurosensory and behavioral endpoints in fish at PCO2 levels predicted to occur before the end of the 21st century (reviewed in: Heuer and Grosell, 2014; Clements and Hunt, 2015; Ashur et al., 2017). Consistent disturbances have been noted across a range of sensory systems, including olfaction (Dixson et al., 2010; Cripps et al., 2011; Devine et al., 2012a,b), hearing (Simpson et al., 2011), and vision (Chung et al., 2014), and have also been implicated in processes linked to general cognitive function including changes in lateralization (Domenici et al., 2012) and learning (Ferrari et al., 2012). The implications of these disturbances are expected to be substantial and include changes to dispersal, recruitment, connectivity, social interactions, predator-prey dynamics, population replenishment, biodiversity, habitat preference, and settlement timing; all of which could drastically affect population and ecosystem dynamics (Heuer and Grosell, 2014).
Evidence suggests that predicted end-of-the-century PCO2 levels will have adverse effects on the olfactory-mediated behavior of most fish species, including reef species (Dixson et al., 2010; Munday et al., 2016a,b), temperate species (Jutfelt et al., 2013), and sharks (Green and Jutfelt, 2014; Dixson et al., 2015). For example, clownfish larvae exposed to 1,000 μatm CO2 did not respond appropriately to the odor of predatory fishes (Dixson et al., 2010); moreover, juvenile damselfish exposed to increased PCO2 and then released into the wild suffered an 8-fold increase in predation (Munday et al., 2010). In both studies, control juveniles avoided predator odors, whereas CO2-treated juveniles were attracted to them (Dixson et al., 2010; Munday et al., 2010), suggesting that mal-adaptive behavioral responses may result from altered sensory processing by the central nervous system. An important step in understanding the basis for neurosensory disruption in fish was reported by Nilsson and colleagues in 2012, who demonstrated that the GABA and GABAA receptors were associated with disruptions in olfaction and lateralization during CO2 exposure (Nilsson et al., 2012). The addition of gabazine, a GABAA receptor antagonist, was found to re-establish adaptive olfactory-mediated behavior and lateralization in CO2-exposed fish (Nilsson et al., 2012). Since then, other studies have demonstrated the ability of gabazine to restore much of the high PCO2 related altered behaviors (Chivers et al., 2014; Hamilton et al., 2014; Lai et al., 2015). To date, disruption in brain neurotransmitter function has been proposed as the main mechanism involved in the behavioral effects of high CO2 in fish (critically reviewed by Tresguerres and Hamilton, 2017). However, the direct effect of high PCO2/low pH on the olfactory system may have been underestimated. Recently, acute exposure to acidified seawater was shown to directly reduce olfactory sensitivity of the sea bass (Dicentrarchus labrax) to some odorants (Porteus et al., 2018); however, the mechanism remains unknown. In addition, neither medium- or long-term effects on sensitivity were not investigated, nor whether the olfactory system was able to respond at the structural level and counteract the effects of acidification. Assessing the occurrence of morphological and/or physiological alterations in response to increased PCO2 is crucial to predict the effects of ocean acidification in fish. Therefore, the main objective of the current work was to evaluate the effects of both acute and longer-term exposure to low pH/high PCO2 on the olfactory sensitivity of the gilthead bream (Sparus aurata), an economically important species. The observed effects on olfactory sensitivity were considered in the context of predicted effects of pH change on the protonation of odorant molecules. Finally, we evaluated whether exposure to high PCO2 provoked structural changes in the olfactory epithelium that might be indicative of adaptation of the olfactory epithelium to lower water pH.
Materials And Methods
Fish Maintenance
Animal maintenance and experimentation were carried out in certified experimental facilities and followed Portuguese national legislation (DL 113/2013) under a “group-1” license by the Veterinary General Directorate, Ministry of Agriculture, Rural Development and Fisheries of Portugal. Gilthead seabream (S. aurata: 180 ± 7 g; 22.8 ± 0.6 cm), hereafter “seabream,” were obtained from a commercial supplier (Maresa–Mariscos de Esteros, SA, Huelva, Spain) and maintained in 1,000 L tanks with continuously running natural seawater, under natural photoperiod and temperature, and fed with commercial pellets (Sparos, Olhão, Portugal). For experimental purposes, seabream were randomly distributed into six 100 L tanks in an open circuit system kept at natural temperature and photoperiod. Three tanks were kept in control conditions (“control fish”), while the other three were kept at low pH by bubbling CO2 gas (see below) into the water and maintained for 4–8 weeks before electrophysiology experiments (“high CO2 fish”). No mortality was observed during this time.
Seawater Chemistry
Seawater was pumped from the ocean into two (2,000 L) header tanks, where it was aerated with ambient air (control) or CO2 to achieve the desired pH (elevated-CO2 treatment). The PCO2 of the CO2 treatment header tank was maintained at the target value of 1,000 μatm using a pH probe connected to an internal controller (EXAxt PH450G, Yokogawa Iberia, Portugal).
Control and high CO2 seawater was supplied to 100 L tanks at 2 L min−1. Seawater pH (Orion star A221, Thermo Scientific, Portugal), temperature (Orion star A221, Thermo Scientific, Portugal), and salinity (WTW, cond3310, Spain) were recorded daily in each tank. Water samples were analyzed for total alkalinity twice a week by Gran titration (DL15 titrator, Mettler Toledo, Portugal) using a certified acid titrant (0.1 M HCl, Fluka Analytical, Sigma-Aldrich). Carbonate chemistry parameters (Table 1) were calculated in CO2SYS (Pierrot et al., 2006) using the constants K1 and K2 from Mehrbach et al. (1973) refit by Dickson and Millero (1987) and Dickson (1990) for KHSO4.
Selection of Odorants
Odorant, in the context of the current study, can be defined as any chemical detected by the olfactory system which, in vertebrates, is the chemosensory system associated with cranial nerve I (discussed in Mollo et al., 2014, 2017). The effect of acute (minutes) and medium-term (4 weeks) exposure to high PCO2 on olfactory sensitivity was tested for amino acids, bile fluid, and intestinal fluid. Olfactory sensitivity to amino acids in fish is believed to be mainly involved in food detection (Hara, 1994; Velez et al., 2007). Bile fluid and intestinal fluid are known to be potent olfactory stimuli in fish, possibly involved in chemically-mediated interactions both within and between species (Velez et al., 2009; Huertas et al., 2010; Buchinger et al., 2014; Fatsini et al., 2017; Hubbard et al., 2017). Therefore, we tested the effect of elevated CO2/low pH on the olfactory responses of seabream to conspecific bile, intestinal fluid, and a range of amino acids in fish maintained under control conditions (pH nominally 8.2; PCO2 nominally 350 μatm) and high CO2 conditions (pH nominally 7.7, PCO2 nominally 1,000 μatm) for 4–8 weeks.
Collection of Body Fluids
Body fluids (bile and intestinal fluid) were collected from five seabream (144 ± 16 g; 18.7 ± 1.0 cm; all males) kept in normal seawater. Fish were anesthetized with 2-phenoxy ethanol (Sigma-Aldrich, Portugal; 2.0 ml L−1) and quickly sacrificed by decapitation. Bile fluid was taken directly from the gall bladder and the intestinal fluid was extracted from the posterior 10 cm of intestine. Samples were pooled, diluted in distilled water (1:2), mixed thoroughly, centrifuged, aliquoted, and frozen (−20°C) until use.
Olfactory Nerve Recording
Seabream were anesthetized in aerated natural seawater containing 300 mg L−1 MS222 (ethyl-3-aminobenzoate methanesulfonate salt, Sigma-Aldrich, Portugal) until the response to a tail pinch had stopped; an intramuscular injection of the neuromuscular blocker gallamine triethiodide (Sigma-Aldrich, Portugal; 10 mg/kg in 0.9% NaCl) was then given. Fish were then placed in a padded V-support and the gills flushed with aerated seawater containing 150 mg l−1 MS222. The pH of the water flushing the gills was around 8.2 for control fish and 7.7 for fish kept in high CO2. Therefore, for control fish, only the olfactory epithelium experienced high CO2 conditions during electrophysiology experiments; water chemistry parameters during experiments are summarized in Table 2.
Table 2. Water chemistry parameters for control and elevated CO2 fish in the electrophysiology set up.
The olfactory rosette was exposed by cutting the skin and connective tissue overlying the nasal cavity. The nostril was constantly irrigated with filtered sea water (without anesthetic) under gravity (flow rate: 6 ml min−1) via a glass tube. Test solutions were delivered to the tube irrigating the nasal cavity via a computer-operated three-way solenoid valve for 4 s. The olfactory nerves were exposed by removal of the skin, connective tissue, and overlying bone (Hubbard et al., 2000, 2003a). Charcoal-filtered sea water was used to make up the odorants solutions and to irrigate the olfactory rosette during experiments. This water was either bubbled with air (control) or CO2 (low pH) until the desired pHNBS was reached. Amino acid solutions were prepared from frozen aliquots of 10−2 M; all stimuli were diluted in charcoal-filtered sea water (either control or high CO2 as appropriate) immediately prior to use. The order in which odorants were given was randomized, but each odorant was always given from the lowest to the highest concentration.
The electrodes were placed in the olfactory nerve in a position that gave maximal response to 10−3 M L-serine. Olfactory nerve activity was recorded using tungsten micro-electrodes (0.1 MΩ, World Precision Instruments, UK) as previously described (Hubbard et al., 2000). Fish were connected to earth via a copper wire inserted in the flank. The raw signal was amplified (x20,000; AC pre-amplifier, Neurolog NL104; Digitimer Ltd., Welwyn Garden City, UK), filtered (high pass: 200 Hz, low pass: 3,000 Hz; Neurolog NL125, Digitimer Ltd.) and integrated (time constant 1 s; Neurolog NL703, Digitimer Ltd.). Raw and integrated signals were digitized (Digidata 1440A, Molecular Devices, San Jose, CA, USA) and recorded on a PC running AxoScope™ software (version 10.6, Molecular Devices).
All integrated response amplitudes were normalized to the amplitude of the integrated response to 10−3 M L-serine (the “standard”). Responses to the standard were recorded regularly at the beginning and end of each group of samples (every three to five samples) throughout the recording session. Each stimulus was applied for 4 s, with at least 1 min between odorants to allow complete recovery of the receptors.
Calculation of Amino Acid Protonation State
Odorant molecules, such as amino acids, can be subjected to ionization processes depending on the pH of their surroundings. Amino acids possess two or more ionizable functional groups that are de-protonated at different pH levels. The pKa value of each ionizable group of a molecule represents the pH at which 50% of the molecules in solution are protonated at this group and 50% remain unchanged. The Henderson-Hasselbalch equation relates pH and the group-specific pKa values for each amino acid, and therefore calculates the abundance of different protonation states across the entire pH range of interest (for details, see Po and Senozan, 2001, and references therein). At any given pH, a mixture of different protonation states can be present. For the current study, we assumed that the amino acid forms with deprotonated head group to be the form binding to the olfactory receptors. We, therefore, calculated the concentration of this form (= effective concentrations) of different amino acids used in the experiments for the control (pH 8.2) and treatment (pH 7.7) conditions based on the concentration used during the bioassay and the group-specific pKa values for each amino acid obtained from the CRC Handbook of Chemistry and Physics (Lide, 2004).
Comparison of pH-Induced Changes in Amino Acid Protonation State to Olfactory Potency
For direct comparison of the change in protonation state and the olfactory response for different concentrations and pH conditions, the respective effective stimulus concentrations were inserted into the linear regression equation/Hill plot obtained from the integrated nerve response at a given pH (see “Data Analysis” for details). The resulting points were plotted in the same figure to visualize the extent of change caused by stimulus protonation vs. the extent of change observed for the nerve response.
Conformers and Charge Distribution of Different Protonation States by Quantum Chemical Calculations
We used quantum chemical calculations to obtain the energetically most favorable model conformers for the two protonation states of interest and assess differences in their molecular electrostatic potential (MEP), which describes the charge distribution around the molecule. Starting from the structures published by Rak et al. (2001), the protonation state models of L-arginine were optimized using the PBE0 exchange correlation functional (Adamo and Barone, 1999) with a pc-2 basis set (Jensen, 2001, 2002a,b) and water as implicit solvent using COSMO (Klamt and Schüürmann, 1993) implemented in the ORCA suite of programs (Neese, 2012) (Version 3.0.0). We used the RIJ-COSX approximation (Neese et al., 2009) with a def2-TZVPP/J auxiliary basis set (Weigend and Ahlrichs, 2005) and included D3 dispersion corrections following Grimme et al. (2010, 2011). The VeryTightSCF and TightOpt criteria implemented in ORCA were used to stop the SCF gradient and the optimization at a total energy change of <10−8 Eh, respectively. Full optimization was confirmed by ensuring that there were no negative numerical frequencies for each conformer. The calculation of the molecular electrostatic potential (MEP) was performed with the GAMESS program (vJan122009R1) using the PBE exchange correlation functional (Perdew et al., 1996) in conjunction with a STO-3G basis set (Hehre et al., 1969; Collins et al., 1976). Three-dimensional electron density iso-surface was visualized with 100 grid points, a medium grid size, and a contour value of 0.03 e·a0−3 using the MacMolPlt program (Bode and Gordon, 1998) (v7.5141). The density iso-surface was colored according to the MEP with a maximum value of 0.25 Eh·e−1 and the RGB color scheme with red representing positive, green neutral, and blue negative charge.
Histology
Seabream were killed rapidly using an anesthetic overdose (500 mg/l MS222). Olfactory rosettes were carefully dissected from 12 fish (six controls, and six in high CO2 water for at least 4 weeks) under a stereo-microscope. Olfactory epithelia were fixed in Bouin’s solution (Sigma, Portugal) for 20 h and washed with 70% ethanol. Tissues were then embedded in paraffin and sectioned longitudinally (5 μm sections). Serial sections of the tissues were collected until the complete olfactory organ was visible (including the central raphe). The dewaxed and rehydrated sections were stained with Mayer’s hematoxylin/eosin and neutral and acidic mucus was identified using a combined periodic acid-Schiff (PAS)-Alcian blue stain. Slides were examined under a microscope (Leica DM 2000, Famalicão, Portugal) and photographic images were obtained using a digital camera (Leica DFC 480, IM50-software). For analysis, we focused on the three central lamellae (lamellae 7, 8, and 9) of the olfactory rosette. The ratio between nonsensory vs. sensory epithelium within an olfactory lamella was assessed by dividing the length of the apical nonsensory epithelium by the length of the olfactory lamella (from the top to the central raphe). Images were analyzed using the software Fiji Image J. Statistically significant differences between the ratio of nonsensory vs. sensory epithelium in control animals and experimental animals was assessed using Student’s t test (Prism 6). The goblet (mucus producing) cells in the apical surface of the lamellae were quantified using Fiji Image J.
Data and Statistical Analysis
All statistical analyses of electrophysiological results were carried out on normalized data. To determine whether there were effects of tank pH (i.e., control vs. high CO2), odorant pH and concentration, an ANOVA analysis was preformed (SPSS statistics). Factors that showed an effect were further analyzed by linear regression or fit to a three-parameter Hill plot.
Olfactory nerve responses to L-serine, L-leucine, L-arginine, bile fluid, and intestinal fluid were analyzed by linear regression; comparison between nerve responses at different tank, and odorant pH was assessed by linear regression of log-transformed data (Hubbard et al., 2003b) and comparing both the slopes and elevations (Prism 6). The slope of each curve indicates how a given increase in odorant concentration evokes a higher amplitude response; this depends on the binding affinity between ligand (odorant in this case) and receptor. Detection thresholds, which also depend on ligand-receptor affinity, were determined for each independent experiment, from the intercept with the x-axis of linear regression fit to individual concentration-response curves for each odorant and each treatment (Prism 6). Differences between detection thresholds were tested using Student’s t test for paired data (log). Responses to L-glutamine and L-glutamic acid were described by a three-parameter Hill curve as previously described (Hubbard et al., 2000). EC50 and Emax were calculated for each fish. These data were then compared using Student’s t-test (Prism 6).
To evaluate whether normalized data could be used to compare olfactory responses between control and high CO2 fish, we first compared raw integrated responses to L-serine at pH 8.2 in control and high CO2 fish. There were no differences in the slope (p = 0.71) or elevation (p = 0.88) of the concentration-response curves to L-serine in control and experimental animals (data not shown); therefore, we compared normalized olfactory nerve responses in control and experimental fish. A significance level of p < 0.05 was used throughout.
Results
For ANOVA analysis, the factors analyzed were the effects of tank pH (i.e., the pH at which the fish were kept: pH 7.7 or 8.2), odorant pH (i.e., the pH of the odorant solutions delivered to the olfactory epithelium: pH 7.7 or 8.2) and odorant concentration (10−3 to 10−6 M for the amino acids). Results are summarized in Table 3. While tank and odorant pH only showed significant effects for some of the odorants, the concentration always had a significant effect. Factors that showed an effect were further analyzed by linear regression or fit to a Hill three-parameter concentration-response curve.
Table 3. ANOVA analysis (Univariate analysis of variance, SPSS) for factors: tank pH, odorant pH, and odorant concentration on the detection of different odorants.
Effects of Odorant pH on Olfactory Sensitivity in Control Fish
In control seabream (kept at pH 8.2 prior to electrophysiological recording), the elevation of the concentration-response curves to L-serine (F45 = 9.363, p < 0.01), L-leucine (F53 = 25.225, p < 0.0001), and L-arginine (F52 = 0.1736, p = 0.68) tested at pH 8.2 was significantly higher than at pH 7.7. However, there were no differences between the slope of the concentration-response curves at pH 8.2 and 7.7 (Figures 1A–C). In addition, exposure to low pH water significantly increased the threshold of detection (Figure 1F) for L-serine (T5 = 3.170, p < 0.05), L-leucine (T6 = 2.848, p < 0.05), and L-arginine (T6 = 2.862, p < 0.05).
Figure 1. Normalized olfactory nerve responses of control fish (kept at pH 8.2) to (A) L-serine, (B) L-leucine, (C) L-arginine, (D) L-glutamine, and (E) L-glutamic acid; under control (red: odorant pH 8.2) and elevated PCO2 conditions (blue: odorant pH 7.7). (F) Effects of acute exposure to elevated PCO2 in the olfactory detection threshold of L-serine, L-leucine, and L-arginine. Values are shown as mean ± S.E.M. *p < 0.05, **p < 0.01, ***p < 0.001; n = 6.
In contrast to L-serine, L-arginine, and L-leucine, the concentration-response curves of L-glutamine and L-glutamate did not fit linear regression well; the best-fitting model was the three-parameter Hill curve (Figures 1D,E). Olfactory nerve responses of control fish to L-glutamine at pH 8.2 and 7.7 were statistically equivalent (Figure 1D). There was no difference (T5 = 2.205, p = 0.08) between the EC50 for L-glutamate of control fish tested at pH 8.2 (1.142 × 10−4 M) and pH 7.7 (3.89 × 10−4 M), but the maximum response (Emax) was significantly (T5 = 2.886, p < 0.05) higher at odorant pH 8.2 (0.88 ± 0.05) than at pH 7.7 (0.78 ± 0.05; Figure 1E).
For bile fluid, there were no differences in slope (F36 = 1.419, p = 0.24) or elevation (F37 = 0.022, p = 0.83) at pH 8.2 and 7.7 (Figure 2A). For intestinal fluid, the slope of the curves at pH 8.2 and 7.7 was not different (F44 = 0.072, p = 0.79); however, the elevation of the curves was statistically higher (F45 = 6.374, p < 0.05) at pH 8.2 than 7.7 (Figure 2B).
Figure 2. Normalized olfactory nerve responses of control fish (kept at pH 8.2) to conspecific bile fluid (A) and intestinal fluid (B) with odorant pH 8.2 (red) and pH 7.7 (blue). Values are represented as mean ± S.E.M; *p < 0.05, n = 6.
Effect of Odorant pH on Olfactory Sensitivity in High CO2 Fish
As for controls, in high CO2 fish (kept at pH 7.7 for at least 4 weeks prior to electrophysiological recording) the elevation of the concentration-response curves to L-serine (F45 = 42.156, p < 0.0001), L-leucine (F45 = 46.971, p < 0.0001), and L-arginine (F43 = 9.888, p < 0.01) tested at pH 8.2 was significantly higher than at pH 7.7 (Figures 3A–C). However, there were no differences between the slopes of the concentration-response curves at pH 8.2 and 7.7 (Figures 3A–C). In addition, exposure to low pH water significantly increased the threshold of detection (Figure 3F) to L-serine (T5 = 3.464, p < 0.05) and L-leucine (T5 = 3.084, p < 0.05) but not to L-arginine (T10 = 1.578, p = 0.07).
Figure 3. Normalized olfactory nerve responses of high CO2 fish (kept at pH 7.7) to (A) L-serine, (B) L-leucine, (C) L-arginine, (D) L-glutamine, and (E) L-glutamic acid; under odorant pH 8.2 (red) and odorant pH 7.7 (blue) conditions. (F) Effects of exposure to elevated PCO2 in the olfactory detection threshold of L-serine, L-leucine, and L-arginine. Values are shown as mean ± S.E.M. *p < 0.05, **p < 0.01, ***p < 0.001; n = 6.
As in the case of control fish, in the high CO2 animals the best-fitting model for the concentration-response curve to L-glutamine and L-glutamate was the three-parameter Hill curve (Figures 3D,E). Sensitivity of high CO2 fish to L-glutamine and L-glutamate at pH 8.2 and 7.7 were statistically equivalent; there were no differences between the EC50 or the Emax for these odorants (Figures 3D,E).
The olfactory nerve responses of high CO2 fish to bile fluid and intestinal fluid at pH 8.2 and 7.7 were statistically equivalent; there were no differences in the slope or elevation of the concentration-response curves (Figure 4).
Figure 4. Normalized olfactory nerve responses of high CO2 fish (kept at pH 7.7) to conspecific (A) bile fluid and (B) intestinal fluid with odorant pH 8.2 (red) and pH 7.7 (blue). Values are shown as mean ± S.E.M; n = 6.
Short-Term (Acute) vs. Medium-Term Effects
When tested at odorant pH 8.2 or 7.7, responses to L-serine, L-leucine, L-arginine, L-glutamine, and bile fluid were similar in control and high CO2 fish; that is, medium-term exposure to high CO2 did not influence the effects of acute changes in odorant pH. As there were significant effects of medium-term exposure to high PCO2 on the olfactory response to L-glutamic acid and intestinal fluid (Table 3), we compared responses of control and high CO2 fish at the same odorant pH. Comparing sensitivity to L-glutamic acid between control and high CO2 fish at odorant pH 8.2 showed that there were no significant differences between EC50 (T9 = 0.543, p = 0.60) or Emax (T9 = 1.814, p = 0.10) (Figure 5A). At odorant pH 7.7, however, the EC50 of control fish was significantly (T10 = 2.357, p < 0.05) higher than that of high CO2 fish. There was no difference between the Emax (T10 = 1.317, p = 0.22) of control and high CO2 fish (Figure 5B).
Figure 5. Normalized olfactory nerve responses of control and high CO2 fish to L-glutamate (A,B) and intestinal fluid (C,D) with the nostril conditioned with water pH 8.2 (A,C) and pH 7.7 (B,D). Values are shown as mean ± S.E.M. *p < 0.05; ***p < 0.001; n = 6.
For intestinal fluid at odorant pH 8.2, differences between the slopes of control and high CO2 fish were highly significant (F44 = 8.15, p < 0.01); it was therefore not possible to compare the intercepts (Figure 5C). For odorant pH 7.7, there was no differences in the slope (F44 = 1.70, p = 0.20) of control and high CO2 fish; however, the elevation was significantly different (F45 = 4.51, p < 0.05); the controls were more sensitive than the high CO2 fish (Figure 5D).
Comparison of pH-Induced Effects on Protonation State and Olfactory Sensitivity
The effect of high CO2 water on the olfactory sensitivity was reversed by exposing the olfactory epithelium to normal seawater for 10 min, with the exception of the decrease in amplitude of response to intestinal fluid seen in high CO2 fish. The transient nature of the observed changes in olfactory sensitivity at low pH, and that the extent of the changes depended on the odorant, raised the hypothesis that the mechanism involved is related to pH-driven changes in protonation of odorants and/or receptor(s), due to the reduction in water pH. To test this hypothesis, and understand the real impact of pH on olfactory sensitivity, the proportion of each protonation state over the pH range was calculated based on the pKa using the Henderson-Hasselbalch equation. Assuming that the protonation state of the ligand (odorant) and/or the binding site of its receptor affect(s) binding affinity (reviewed in Strasser et al., 2017), we calculated the stimulus concentrations at different states of protonation (Table 4).
Table 4. Calculated effective stimulus concentrations during bioassays at pH 8.2 and 7.7 as well as the respective change in stimulus concentration in percentage between pH-conditions.
These values were introduced into the linear regression or three-parameter Hill equation obtained experimentally from the response to odorants at a pH 8.2. The resulting points were plotted to simulate the extent of change caused by stimulus protonation state and compared with the extent of change observed for the olfactory nerve response at pH 7.7.
For all amino acids tested, a reduction in pH caused a reduction in the concentration of non-protonated ligand. The concentrations of active ligand at pH 7.7 were 32–35% lower than at pH 8.2 (Table 4). For L-serine, L-leucine, and L-arginine, the reduction in the olfactory response recorded was similar to that predicted by the theoretical protonation model (Figures 6A–C). However, for L-glutamine, a decrease in olfactory sensitivity was predicted (Figure 6D) but there was little observed difference between pH 8.2 and 7.7. Finally, for L-glutamic acid the decrease in the olfactory nerve response was lower than predicted by the protonation model (Figure 6E).
Figure 6. Normalized olfactory nerve responses of control fish (red) with the nostril conditioned with water pH 7.7 and expected response (blue), based on the effective odorant concentration at pH 7.7 and the linear regression equation or three-parameter Hill equation fit (A) L-serine, (B) L-leucine, (C) L-arginine, (D) L-glutamine, and (E) L-glutamic acid. Values are shown as mean ± S.E.M.
It is well-known that pH affects the protonation state of amino acids. This protonation changes the overall molecular charge and the charge distribution, which could affect interaction of the odorant with its receptor. To assess this, quantum chemical calculations were performed to obtain the energetically most favorable conformers and then the charge distribution for the two relevant protonation states of L-arginine were obtained (Figure 7) as an example. While the positive charge (red) of the side chain remains, the distribution of the negative charge (blue) around the head group shifts with protonation of the amino group.
Figure 7. Active (ArginineD, left) and protonated (ArginineZ, right) conformers of L-arginine with their molecular electrostatic potential mapped onto an iso-electron density surface. Negative charge is colored in blue, neutral in green, and positive charge in red. The orange circles highlight differences in the charge distribution between the two conformer models.
Histology
The olfactory epithelia of fish kept for 4–8 weeks in water of pH 7.7 showed pronounced morphological changes, including an increase in the area of the non-sensory epithelium, and an increase in the number of mucus cells. The oval olfactory rosette (OR) consists of 32–38 lamellae (OL) that radiate from a median raphe (MR). There were no marked differences in the gross anatomy of the olfactory organ between fish kept in control and high CO2 water. Each lamella has two layers of epithelium that enclose the central core (CC). The central core is separated from the epithelium by the basement membrane (BM). The non-sensory area consists of a stratified epithelium covering the outer margins of the lamella (NSE). The sensory area consists of pseudostratified epithelium covering lateral sides of the lamella (SE) (Figures 8C,D).
Figure 8. Representative histological sections of the olfactory lamellae of the gilthead seabream stained with periodic acid-Schiff/Alcian blue to show the effects of exposure to high PCO2/low pH on tissue morphology. (A,B) control fish; (C,D) high CO2 fish. CC, Central core; MC, Mucus cell; NSE, non-sensory epithelium; SE, sensory epithelium.
In high CO2 fish (pH 7.7), the ratio between non-sensory and sensory epithelium (0.34 ± 0.02) is significantly higher (T8 = 5.583, p < 0.001) than in controls (0.21 ± 0.01). Furthermore, the number of mucus cells per lamellae was significantly higher (T8 = 3.786, p < 0.01) in high CO2 fish (21.6 ± 1.7) than controls (11.8 ± 2.0) (Figure 8). Interestingly, mucus cells are present only in the nonsensory epithelium; the increased number of these cells in high CO2 seabream may be related to the augmented area of non-sensory epithelium. In addition, to the higher number of mucus cells, the color of these cells, when examined for neutral and acidic mucus by combined Periodic Acid-Schiff (PAS)-Alcian blue staining, was deep pink in high CO2 fish and colorless or light pink in controls (Figure 8); these differences suggest a shift in pH of the mucus from acidic to more neutral.
Discussion
The current study shows that, in the seabream, increased oceanic PCO2 and/or reduction of pH predicted for the end of the century cause(s) a significant reduction of olfactory sensitivity to some, but not all, odorants. Significantly, this reduction persisted in seabream kept in high PCO2 water for several weeks; fish were unable to compensate. The effect of exposure to high CO2/low pH was complex and varied with the odorant type. The change in olfactory sensitivity can be partly – but not entirely – explained by a shift in the protonation state of the odorant molecule. Moreover, the current study shows that after several weeks exposure to high PCO2 there were morphological changes in the olfactory epithelium indicative that exposure to high CO2/low pH evokes a structural re-organization.
Overall, acute exposure of the olfactory epithelium to high CO2/low pH water reduced the olfactory response for five out of seven odorants tested, and increased the detection threshold in three out of those five, similar to previously reported in bass (Porteus et al., 2018). Importantly, however, the current study also showed that this effect persisted, even after 4 weeks exposure to acidified water. There was an increase in the threshold of detection and a decrease in the amplitude of responses to L-serine, L-leucine and L-arginine, but not to bile fluid. This means that the fish would have to be closer to an odor source in order to detect it. This would have obvious repercussions in food-search and predator-avoidance behaviors, but may also help explain disrupted olfactory-guided homing and/or migration (Munday et al., 2009; Devine et al., 2012a). There were differences in the amplitude of response to intestinal fluid, but the detection threshold did not change. As intestinal fluid is likely a mixture of different odorants, this may reflect a reduction in all or just some; if the latter, then this may imply a qualitative effect on olfactory sensitivity, as well as a quantitative effect. When tested at odorant pH 8.2 or 7.7, responses to L-serine, L-leucine, L-arginine, L-glutamine, and bile fluid were similar in control and high CO2 fish. This suggests that medium term exposure to high CO2 does not cause changes in the olfactory sensitivity to these odorants per se; rather, the seabream is unable to compensate for the effect of low pH/high PCO2 on acute sensitivity. However, at odorant pH 7.7, olfactory sensitivity to glutamate was higher in high CO2 fish than in controls, suggesting that, in this specific case, the sensitivity to L-glutamate can somehow acclimate to high CO2/low pH. This may involve changes in the expression and/or type of receptor(s) in the olfactory epithelium; downregulation, but not upregulation of olfactory receptors was seen in the bass (Porteus et al., 2018), whereas no such changes were seen in the salmon (Williams et al., 2019). In contrast, the sensitivity to intestinal fluid when tested in seawater at pH 8.2, was higher in control fish than in high CO2 fish. This suggests that medium-term exposure to CO2 impairs the olfactory detection of intestinal fluid and, unlike the other stimuli tested, this was not reversible by bathing the olfactory epithelium with normal seawater. Again, this is likely to reflect a change in receptor expression. Thus, the effect of increased PCO2 on olfactory sensitivity varies with the odorant type.
The sensitivity to some odorants decreased as soon as the olfactory epithelium was exposed to high CO2 water; this effect occurred in both control and high CO2 fish and was reversed by bathing the epithelium with seawater at pH 8.2. We propose that the effect can be directly linked to the effects of pH on the odorant molecules and/or their receptors in fish as well as crabs (Roggatz et al., 2016). We show the changes in charge distribution can alter the conformation of protonated L-arginine, for example. Changes in protonation state are reversible and much more rapid (<micro-seconds) than any physiological and/or behavioral change which may take several days (for example, see Munday et al., 2009, 2010; Hamilton et al., 2014). Olfactory sensitivity in vertebrates, including fish, is mediated by G-protein coupled receptors (GPCRs) (Buck and Axel, 1991; Ngai et al., 1993). Changes in pH have been shown to alter GPCR-ligand affinity in vitro (Ruan et al., 2009), and charged amino acid residues in the receptor play multiple roles in receptor-ligand interaction (Strasser et al., 2017).
Another effect of medium term exposure to high CO2 water was the reduction of olfactory sensitivity to odorants in intestinal fluid. This may be due to qualitative and/or quantitative changes in the olfactory epithelium, as observed in the present study, which we hypothesize may cause a change in the expression of some odorant receptors (as in the bass; Porteus et al., 2018). Finally, in the case of L-glutamate, the decrease in olfactory sensitivity in high CO2 water was lower in animals chronically exposed than in those acutely exposed. This may show a degree of adaptation to high PCO2 and may be a result of morphological re-organization and/or changes in olfactory receptor expression. Further studies are required to evaluate whether exposure to high CO2 water induces changes in gene expression in the olfactory epithelium. Nonetheless, it should be noted that it is difficult to separate the direct effect of increased PCO2 and the consequent reduction in pH (Ou et al., 2015). A recent study in coho salmon (Williams et al., 2019), while showing that increased PCO2 affects both olfactory-mediated behavior and gene expression in the olfactory system, suggested that olfactory sensitivity is unaffected by increased PCO2/low pH at the level of the olfactory epithelium. However, this study used a different experimental approach (i.e., olfactory responses at different PCO2/pH levels were not compared in the same fish), and the high concentration (10−2 M) of odorants used, coupled with the lower sensitivity of the electro-olfactogram (EOG) when used in seawater (Silver et al., 1976; Hubbard et al., 2011) makes this negative finding difficult to interpret.
Does Exposure to High PCO2 Modify the Olfactory Epithelium Organization?
Medium term (4 weeks) exposure to high CO2 induced an increase in the ratio between nonsensory vs. sensory epithelium and an augmentation in the number of mucus cell per lamellae. Furthermore, there was a shift in pH of the mucus from acidic to more neutral. This suggests that high CO2 seawater evokes modification of the epithelium and causes increased proliferation of neutral mucus-producing goblet cells to reinforce the outer extrinsic barrier. This may be an adaptive response, as the mucus will protect the underlying epithelial cells, including the olfactory receptor neurones, but the charge, viscosity and the thickness of the layer may also reduce the odorant/receptor interaction and thus lower sensitivity (Rygg et al., 2013; Corfield, 2015). However, the olfactory sensitivity of seabream exposed for 4 weeks to high PCO2 was similar to that of controls, when tested at the same pH. The enlarged area of non-sensory epithelium might also represent an adaptation to high CO2; an augmented area of non-sensory epithelium positioned in the upper and outer margin of the olfactory lamellae results in a localization of the sensory area deeper in the olfactory lamellae, and as a consequence it is less exposed. On the other hand, it could represent impairment of olfactory epithelium formation, growth and/or renewal, and ultimately it might contribute to decreased olfactory sensitivity given that the sensory area is reduced. Initial studies clearly indicate changes in gene expression induced by ocean acidification in the olfactory system (Porteus et al., 2018; Williams et al., 2019). Further studies are required to evaluate the mucin forms expressed, the factors driving the increased number of goblet cells, necrosis and/or apoptosis of the epithelium and changes in the number, type and distribution of olfactory receptor neurones and/or receptors and ion pumps and channels within the olfactory epithelium.
Does Seawater pH Change the Protonation of Odorant Molecules and Thereby Modify Olfactory Sensitivity?
The transient nature of the observed changes in olfactory sensitivity raised the hypothesis it is due to changes in the protonation state of odorant and/or receptor, due to changes in water pH, as previously suggested for the crab pheromone (Roggatz et al., 2016). Fish odorants are mostly low molecular weight organic molecules with high water solubility (Laberge and Hara, 2004), such as amines (Rolen et al., 2003), amino acids, nucleotides, nucleosides and organic acids (reviewed in Hara, 1994). Amino acids contain both functional groups attached to the same carbon (α-carbon) and the various alpha amino acids differ in the side chain attached to the α-carbon. Amino acid detection by the olfactory epithelium probably involves different receptors that are relatively ligand-specific and for activation require specific molecular features (Friedrich and Korsching, 1997, 1998; Fuss and Korsching, 2001). Our hypothesis is that change of pH in future oceans could lead to changes in protonation states of amino acids (and other odorants) which, in turn, may lead to significant alterations in their structure and affinity for their receptor(s). Herein, we determined that the difference in pH from 8.1 to 7.7 would reduce the concentration of non-protonated amino acids by up to 35%. This coincided with a reduction in olfactory response, except for L-glutamine. This may be because, for L-glutamine, the binding site of the receptor is not dependent on the protonation state of the ligand. Furthermore, for L-glutamic acid, reductions in olfactory responses were greater than predicted by the change in protonated/non protonated ratio. This could be explained by the occurrence of changes of charge distribution within the binding site of the receptor (Tierney and Atema, 1988; Strasser et al., 2017); thus, conformational changes in both odorant and receptor may contribute to the decrease in olfactory sensitivity. However, other mechanisms – such as effects on the transduction pathways (Rong et al., 2018) – may also play a role.
In conclusion, the current study shows that predicted future oceanic pH will influence fish olfactory sensitivity by acting directly on the olfactory epithelium. Exposure to high CO2 water decreased olfactory nerve responses to most of the odorants tested. Furthermore, seabream were unable to compensate for high CO2 levels and, even after 4 weeks exposure to high, the olfactory sensitivity remained reduced. Our hypothesis is that changes of charge, structure and consequently function of odorant molecules and/or their receptors due to pH change reduces receptor-ligand affinity, as has previously been shown in crabs (Roggatz et al., 2016). Thus, we propose that this is a purely physicochemical issue and explains why fish are unable to compensate. However, changes in the protonation state of odorants is not the only mechanism involved in olfactory disruption due to high CO2 water; long exposure of seabream led to morphological changes in the olfactory epithelium and diminution of olfactory sensitivity to intestinal fluid. Importantly, the protonation effect and the modification of epithelial structure presumably act in combination with the reported alterations in brain neurotransmitter function (Nilsson et al., 2012) and suggest the increase in CO2 in seawater acts at multiple levels in sensory systems and will give rise to complex and diverse outcomes (Kelley et al., 2018). Further studies are required to evaluate whether changes in olfactory receptor expression or modifications in the olfactory epithelium are contributing to the long term effects of CO2 on the olfactory sensitivity of seabream and, probably, other fish.
Data Availability
The datasets generated for this study are available on request to the corresponding author.
Ethics Statement
Animal maintenance and experimentation was carried out in certified experimental facilities and followed Portuguese national legislation (DL 113/2013) under a “group-1” license by the Veterinary General Directorate, Ministry of Agriculture, Rural Development and Fisheries of Portugal.
Author Contributions
ZV, JH, and PH contributed to conception and design of the study. ZV and PH designed and carried out the electrophysiology experiments. CR and DB designed and carried out the chemical modeling. ZV did the histology. ZV, CR, and PH wrote the first draft of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This work was supported by grants from the Fundação para a Ciência e a Tecnologia (FCT), Portugal (PTDC/BIA-BMA/30262/2017, SFRH/BPD/100409/2014 to ZV and UID/Multi/04326/2019).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We acknowledge the Viper High Performance Computing facility of the University of Hull and its support team and funding for CR via the Energy and Environment Institute, University of Hull. Furthermore, thanks are due to Victoria Clark (University of Hull) for initial molecular optimizations of the L-arginine models used here.
References
Adamo, C., and Barone, V. (1999). Toward reliable density functional methods without adjustable parameters: the PBE0 model. J. Chem. Phys. 110, 6158–6170. doi: 10.1063/1.478522
Ashur, M. M., Johnston, N. K., and Dixson, D. L. (2017). Impacts of ocean acidification on sensory function in marine organisms. Integr. Comp. Biol. 57, 63–80. doi: 10.1093/icb/icx010
Bode, B. M., and Gordon, M. S. (1998). MacMolPlt: a graphical user interface for GAMESS. J. Mol. Graph. Model. 16, 133–138. doi: 10.1016/S1093-3263(99)00002-9
Buchinger, T. J., Li, W., and Johnson, N. S. (2014). Bile salts as semiochemicals in fish. Chem. Senses 39, 647–654. doi: 10.1093/chemse/bju039
Buck, L., and Axel, R. (1991). A novel multigene family may encode odorant receptors: a molecular basis for odor recognition. Cell 65, 175–187. doi: 10.1016/0092-8674(91)90418-X
Caldeira, K., and Wickett, M. E. (2003). Anthropogenic carbon and ocean pH. Nature 425:365. doi: 10.1038/425365a
Chivers, D. P., McCormick, M. I., Nilsson, G. E., Munday, P. L., Watson, S.-A., Meekan, M. G., et al. (2014). Impaired learning of predators and lower prey survival under elevated CO2: a consequence of neurotransmitter interference. Glob. Chang. Biol. 20, 515–522. doi: 10.1111/gcb.12291
Chung, W.-S., Marshall, N. J., Watson, S.-A., Munday, P. L., and Nilsson, G. E. (2014). Ocean acidification slows retinal function in a damselfish through interference with GABAA receptors. J. Exp. Biol. 217, 323–326. doi: 10.1242/jeb.092478
Clements, J. C., and Hunt, H. L. (2015). Marine animal behaviour in a high CO2 ocean. Mar. Ecol. Prog. Ser. 536, 259–279. doi: 10.3354/meps11426
Collins, J. B., Von, R., Schleyer, P., Binkley, J. S., and Pople, J. A. (1976). Self-consistent molecular orbital methods. XVII. Geometries and binding energies of second-row molecules. A comparison of three basis sets. J. Chem. Phys. 64, 5142–5151. doi: 10.1063/1.432189
Corfield, A. P. (2015). Mucins: a biologically relevant glycan barrier in mucosal protection. Biochim. Biophys. Acta 1850, 236–252. doi: 10.1016/j.bbagen.2014.05.003
Cripps, I. L., Munday, P. L., and McCormick, M. I. (2011). Ocean acidification affects prey detection by a predatory reef fish. PLoS One 6:e22736. doi: 10.1371/journal.pone.0022736
Devine, B. M., Munday, P. L., and Jones, G. P. (2012a). Homing ability of adult cardinalfish is affected by elevated carbon dioxide. Oecologia 168, 269–276. doi: 10.1007/s00442-011-2081-2
Devine, B. M., Munday, P. L., and Jones, G. P. (2012b). Rising CO2 concentrations affect settlement behaviour of larval damselfishes. Coral Reefs 31, 229–238. doi: 10.1007/s00338-011-0837-0
Dickson, A. G. (1990). Standard potential of the reaction: AgCl(s) +½H, (g) = Ag(s) + HCl(aq), and the standard acidity constant of the ion HSO4 in synthetic sea water from 273.15 to 318.15 K. J. Chem. Thermodyn. 22, 113–127. doi: 10.1016/0021-9614(90)90074-z
Dickson, A. G., and Millero, F. J. (1987). A comparison of the equilibrium constants for the dissociation of carbonic acid in seawater media. Deep Sea Res. Part A Oceanograph. Res. Pap. 34, 1733–1743. doi: 10.1016/0198-0149(87)90021-5
Dixson, D. L., Jennings, A. R., Atema, J., and Munday, P. L. (2015). Odor tracking in sharks is reduced under future ocean acidification conditions. Glob. Chang. Biol. 21, 1454–1462. doi: 10.1111/gcb.12678
Dixson, D. L., Munday, P. L., and Jones, G. P. (2010). Ocean acidification disrupts the innate ability of fish to detect predator olfactory cues. Ecol. Lett. 13, 68–75. doi: 10.1111/j.1461-0248.2009.01400.x
Domenici, P., Bridie, A., McCormick, M. I., and Munday, P. L. (2012). Elevated carbon dioxide affects behavioural lateralization in a coral reef fish. Biol. Lett. 8, 78–81. doi: 10.1098/rsbl.2011.0591
Doney, S. C., Fabry, V. J., Feely, R. A., and Kleypas, J. A. (2009). Ocean acidification: the other CO2 problem. Annu. Rev. Mar. Sci. 1, 169–192. doi: 10.1146/annurev.marine.010908.163834
Fatsini, E., Carazo, I., Chauvigné, F., Manchado, M., Cerdà, J., Hubbard, P. C., et al. (2017). Olfactory sensitivity of the marine flatfish Solea senegalensis to conspecific body fluids. J. Exp. Biol. 220, 2057–2065. doi: 10.1242/jeb.150318
Ferrari, M. C. O., Manassa, R. P., Dixson, D. L., Munday, P. L., McCormick, M. I., Meekan, M. G., et al. (2012). Effects of ocean acidification on learning in coral reef fishes. PLoS One 7:e31478. doi: 10.1371/journal.pone.0031478
Friedrich, R. W., and Korsching, S. I. (1997). Combinatorial and chemotopic odorant coding in the zebrafish olfactory bulb visualized by optical imaging. Neuron 18, 737–752. doi: 10.1016/S0896-6273(00)80314-1
Friedrich, R. W., and Korsching, S. I. (1998). Chemotopic, combinatorial, and noncombinatorial odorant representations in the olfactory bulb revealed using a voltage-sensitive axon tracer. J. Neurosci. 18, 9977–9988.
Fuss, S. H., and Korsching, S. I. (2001). Odorant feature detection: activity mapping of structure response relationships in the zebrafish olfactory bulb. J. Neurosci. 21, 8396–8407. doi: 10.1523/JNEUROSCI.21-21-08396.2001
Green, L., and Jutfelt, F. (2014). Elevated carbon dioxide alters the plasma composition and behaviour of a shark. Biol. Lett. 10:20140538. doi: 10.1098/rsbl.2014.0538
Grimme, S., Antony, J., Ehrlich, S., and Krieg, H. (2010). A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 132, 154104. doi: 10.1063/1.3382344
Grimme, S., Ehrlich, S., and Goerigk, L. (2011). Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 32, 1456–1465. doi: 10.1002/jcc.21759
Hamilton, T. J., Holcombe, A., and Tresguerres, M. (2014). CO2-induced ocean acidification increases anxiety in rockfish via alteration of GABAA receptor functioning. Proc. R. Soc. B 281:20132509. doi: 10.1098/rspb.2013.2509
Hara, T. J. (1994). The diversity of chemical stimulation in fish olfaction and gustation. Rev. Fish Biol. Fish. 4, 1–35. doi: 10.1007/BF00043259
Hehre, W. J., Stewart, R. F., and Pople, J. A. (1969). Self-consistent molecular-orbital methods. I. Use of Gaussian expansions of slater-type atomic orbitals. J. Chem. Phys. 51, 2657–2664. doi: 10.1063/1.1672392
Heuer, R. M., and Grosell, M. (2014). Physiological impacts of elevated carbon dioxide and ocean acidification on fish. Am. J. Phys. Regul. Integr. Comp. Phys. 307, R1061–R1084. doi: 10.1152/ajpregu.00064.2014
Hubbard, P. C., Baduy, F., Saraiva, J. L., Guerreiro, P. M., and Canário, A. V. M. (2017). High olfactory sensitivity to conspecific intestinal fluid in the chameleon cichlid Australoheros facetus: could faeces signal dominance? J. Fish Biol. 90, 2148–2156. doi: 10.1111/jfb.13297
Hubbard, P. C., Barata, E. N., and Canário, A. V. M. (2000). Olfactory sensitivity to changes in environmental [Ca2+] in the marine teleost Sparus aurata. J. Exp. Biol. 203, 3821–3829.
Hubbard, P. C., Barata, E. N., and Canário, A. V. M. (2003a). Olfactory sensitivity of the gilthead seabream (Sparus auratus L) to conspecific body fluids. J. Chem. Ecol. 29, 2481–2498. doi: 10.1023/A:1026357917887
Hubbard, P. C., Barata, E. N., and Canário, A. V. M. (2003b). Olfactory sensitivity to catecholamines and their metabolites in the goldfish. Chem. Senses 28, 207–218. doi: 10.1093/chemse/28.3.207
Hubbard, P. C., Barata, E. N., Ozório, R. O. A., Valente, L. M. P., and Canário, A. V. M. (2011). Olfactory sensitivity to amino acids in the blackspot seabream (Pagellus bogaraveo): a comparison between olfactory receptor recording techniques in seawater. J. Comp. Physiol. A. 197, 839–849. doi: 10.1007/s00359-011-0646-5
Huertas, M., Hagey, L., Hofmann, A. F., Cerdà, J., Canário, A. V. M., and Hubbard, P. C. (2010). Olfactory sensitivity to bile fluids and bile acids in eel (Anguilla anguilla), goldfish (Carassius auratus) and Mozambique tilapia (Oreochromis mossambicus) suggests a ‘broad range’ of sensitivity not confined to those produced by con-specifics alone. J. Exp. Biol. 213, 308–317. doi: 10.1242/jeb.033142
Jensen, F. (2001). Polarization consistent basis sets: principles. J. Chem. Phys. 115, 9113–9125. doi: 10.1063/1.1413524
Jensen, F. (2002a). Polarization consistent basis sets. II. Estimating the Kohn-Sham basis set limit. J. Chem. Phys. 116, 7372–7379. doi: 10.1063/1.1465405
Jensen, F. (2002b). Polarization consistent basis sets. III. The importance of diffuse functions. J. Chem. Phys. 117, 9234–9240. doi: 10.1063/1.1515484
Jutfelt, F., De Souza, K. B., Vuylsteke, A., and Sturve, J. (2013). Behavioural disturbances in a temperate fish exposed to sustained high-CO2 levels. PLoS One 8:e65825. doi: 10.1371/journal.pone.0065825
Kelley, J. L., Chapuis, L., Davies, W. I. L., and Collin, S. P. (2018). Sensory system responses to human-induced environmental change. Front. Ecol. Evol. 6:95. doi: 10.3389/fevo.2018.00095
Klamt, A., and Schüürmann, G. (1993). COSMO: a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient. J. Chem. Soc. Perkin Trans. 2, 799–805. doi: 10.1039/p29930000799
Laberge, F., and Hara, T. J. (2004). Electrophysiological demonstration of independent olfactory receptor types and associated neuronal responses in the trout olfactory bulb. Comp. Biochem. Physiol. 137A, 397–408. doi: 10.1055/s-2004-861558
Lai, F., Jutfelt, F., and Nilsson, G. E. (2015). Altered neurotransmitter function in CO2-exposed stickleback (Gasterosteus aculeatus): a temperate model species for ocean acidification research. Conserv. Physiol. 3:cov018. doi: 10.1093/conphys/cov018
Lide, D. R. (2004). CRC Handbook of chemistry and physics. 85th edn. (Boca Raton, Florida: CRC Press).
Mehrbach, C., Culberson, C. H., Hawley, J. E., and Pytkowicz, R. M. (1973). Measurement of apparent dissociation constants of carbonic acid in seawater at atmospheric pressure. Limnol. Oceanogr. 18, 897–907. doi: 10.4319/lo.1973.18.6.0897
Mollo, E., Fontana, A., Roussis, V., Polese, G., Amodeo, P., and Ghiselin, M. T. (2014). Sensing marine biomolecules: smell, taste, and the evolutionary transition from aquatic to terrestrial life. Front. Chem. 2:92. doi: 10.3389/fchem.2014.00092
Mollo, E., Garson, M. J., Polese, G., Amodeo, P., and Ghiselin, M. T. (2017). Taste and smell in aquatic and terrestrial environments. Nat. Prod. Rep. 34, 496–513. doi: 10.1039/c7np00008a
Munday, P. L., Dixson, D. L., Donelson, J. M., Jones, G. P., Pratchett, M. S., Devitsina, G. V., et al. (2009). Ocean acidification impairs olfactory discrimination and homing ability in a marine fish. Proc. Natl. Acad. Sci. USA 106, 1848–1852. doi: 10.1073/pnas.0809996106
Munday, P. L., Dixson, D. L., Mccormick, M. I., Meekan, M., Ferrari, M. C. O., and Chivers, D. P. (2010). Replenishment of fish populations is threatened by ocean acidification. Proc. Natl. Acad. Sci. USA 107, 12930–12934. doi: 10.1073/pnas.1004519107
Munday, P. L., Watson, S.-A., Parsons, D. M., King, A., Barr, N. G., Mcleod, I. M., et al. (2016a). Effects of elevated CO2 on early life history development of the yellowtail kingfish, Seriola lalandi, a large pelagic fish. ICES J. Mar. Sci. 73, 641–649. doi: 10.1093/icesjms/fsv210
Munday, P. L., Welch, M. J., Allan, B. J. M., Watson, S.-A., McMahon, S. J., and McCormick, M. I. (2016b). Effects of elevated CO2 on predator avoidance behaviour by reef fishes is not altered by experimental test water. Peerj 4:e2501. doi: 10.7717/peerj.2501
Neese, F. (2012). The ORCA program system. Wiley Interdisci. Rev. Comput. Mol. Sci. 2, 73–78. doi: 10.1002/wcms.81
Neese, F., Wennmohs, F., Hansen, A., and Becker, U. (2009). Efficient, approximate and parallel Hartree-Fock and hybrid DFT calculations. A ‘chain-of-spheres’ algorithm for the Hartree-Fock exchange. Chem. Phys. 356, 98–109. doi: 10.1016/j.chemphys.2008.10.036
Ngai, J., Dowling, M. M., Buck, L., Axel, R., and Chess, A. (1993). The family of genes encoding odorant receptors in the channel catfish. Cell 72, 657–666. doi: 10.1016/0092-8674(93)90395-7
Nilsson, G. E., Dixson, D. L., Domenici, P., McCormick, M. I., Sorensen, C., Watson, S.-A., et al. (2012). Near-future carbon dioxide levels alter fish behaviour by interfering with neurotransmitter function. Nat. Clim. Chang. 2, 201–204. doi: 10.1038/nclimate1352
Ou, M., Hamilton, T. J., Eom, J., Lyall, E. M., Gallup, J., Jiang, A., et al. (2015). Responses of pink salmon to CO2-induced aquatic acidification. Nat. Clim. Chang. 5, 950–955. doi: 10.1038/nclimate2694
Perdew, J. P., Burke, K., and Ernzerhof, M. (1996). Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868. doi: 10.1103/PhysRevLett.77.3865
Pierrot, D., Lewis, E., and Wallace, D. (2006). “MS Excel programme developed for CO2 system calculations” in ORNL/CDIAC-105a, Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, US Department of Energy, Oak Ridge, Tennessee. (Oak Ridge, Tennesee, USA).
Po, H. N., and Senozan, N. M. (2001). The Henderson-Hasselbalch equation: its history and limitations. J. Chem. Educ. 78, 1499–1503. doi: 10.1021/ed078p1499
Porteus, C. S., Hubbard, P. C., Uren Webster, T. M., Van Aerle, R., Canário, A. V. M., Santos, E. M., et al. (2018). Near-future CO2 levels impair the olfactory system of a marine fish. Nat. Clim. Chang. 8, 737–743. doi: 10.1038/s41558-018-0224-8
Rak, J., Skurski, P., Simons, J., and Gutowski, M. (2001). Low-energy tautomers and conformers of neutral and protonated arginine. J. Am. Chem. Soc. 123, 11695–11707. doi: 10.1021/ja011357l
Roggatz, C. C., Lorch, M., Hardege, J. D., and Beniot, D. M. (2016). Ocean acidification affects marine chemical communication by changing structure and function of peptide signalling molecules. Glob. Chang. Biol. 22, 3914–3926. doi: 10.1111/gcb.13354
Rolen, S. H., Sorensen, P. W., Mattson, D., and Caprio, J. (2003). Polyamines as olfactory stimuli in the goldfish Carassius auratus. J. Exp. Biol. 206, 1683–1696. doi: 10.1242/jeb.00338
Rong, J., Su, W., Guan, X., Shi, W., Zha, S., He, M., et al. (2018). Ocean acidification impairs foraging behavior by interfering with olfactory neural signal transduction in black sea bream, Acanthopagrus schlegelii. Front. Physiol. 9:1592. doi: 10.3389/fphys.2018.01592
Ruan, K.-H., Cervantes, V., and Wu, J. (2009). Ligand-specific conformation determines agonist activation and antagonist blockade in purified human thromboxane A2 receptor. Biochemistry 48, 3157–3165. doi: 10.1021/bi801443g
Rygg, A. D., Van Duin, A. C. T., and Craven, B. A. (2013). Molecular dynamics simulations of water/mucus partition coefficients for feeding stimulants in fish and the implications for olfaction. PLoS One 8:e72271. doi: 10.1371/journal.pone.0072271
Silver, W. L., Caprio, J., Blackwell, J. F., and Tucker, D. (1976). The underwater electro-olfactogram: a tool for the study of the sense of smell of marine fishes. Exp. Dermatol. 32, 1216–1217.
Simpson, S. D., Munday, P. L., Wittenrich, M. L., Manassa, R., Dixson, D. L., Gagliano, M., et al. (2011). Ocean acidification erodes crucial auditory behaviour in a marine fish. Biol. Lett. 7, 917–920. doi: 10.1098/rsbl.2011.0293
Strasser, A., Wittmann, H.-J., and Seifert, R. (2017). Binding kinetics and pathways of ligands to GPCRs. Trends Pharmacol. Sci. 38, 717–732. doi: 10.1016/j.tips.2017.05.005
Tierney, A. J., and Atema, J. (1988). Amino acid chemoreception: effects of pH on receptors and stimuli. J. Chem. Ecol. 14, 135–141. doi: 10.1007/BF01022537
Tresguerres, M., and Hamilton, T. J. (2017). Acid base physiology, neurobiology and behaviour in relation to CO2-induced ocean acidification. J. Exp. Biol. 220, 2136–2148. doi: 10.1242/jeb.144113
Velez, Z., Hubbard, P. C., Hardege, J. D., Barata, E. N., and Canário, A. V. M. (2007). The contribution of amino acids to the odour of a prey species in the Senegalese sole (Solea senegalensis). Aquaculture 265, 336–342. doi: 10.1016/j.aquaculture.2007.02.029
Velez, Z., Hubbard, P. C., Welham, K., Hardege, J. D., Barata, E. N., and Canário, A. V. M. (2009). Identification, release and olfactory detection of bile acids in the intestinal fluid of the Senegalese sole (Solea senegalensis). J. Comp. Physiol. A. 195, 691–698. doi: 10.1007/s00359-009-0444-5
Weigend, F., and Ahlrichs, R. (2005). Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: design and assessment of accuracy. Phys. Chem. Chem. Phys. 7, 3297–3305. doi: 10.1039/b508541a
Williams, C. R., Dittman, A. H., McElhany, P., Busch, D. S., Maher, M. T., Bammler, T. K., et al. (2019). Elevated CO2 impairs olfactory-mediated neural and behavioral responses and gene expression in ocean-phase coho salmon (Oncorhynchus kisutch). Glob. Chang. Biol. 25, 963–977. doi: 10.1111/gcb.14532
Keywords: olfaction, ocean acidification, fish, amino acid, receptor, olfactory epithelium, carbon dioxide, protonation
Citation: Velez Z, Roggatz CC, Benoit DM, Hardege JD and Hubbard PC (2019) Short- and Medium-Term Exposure to Ocean Acidification Reduces Olfactory Sensitivity in Gilthead Seabream. Front. Physiol. 10:731. doi: 10.3389/fphys.2019.00731
Edited by:
Eduardo Luís Cupertino Ballester, Federal University of Paraná, BrazilReviewed by:
Gianluca Polese, University of Naples Federico II, ItalyLorenzo Gallus, University of Genoa, Italy
Copyright © 2019 Velez, Roggatz, Benoit, Hardege and Hubbard. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peter C. Hubbard, cGh1YmJhcmRAdWFsZy5wdA==