The linkage between the underlying mechanisms behind the fish physiological adaptation to long-term stressors and reliable indicators able to measure unambiguously the stress condition of the fish, is a longstanding problematic in fish welfare research. Within this context, high-throughput technologies like proteomics are powerful techniques for the discovery of new potential biomarkers and lever the aquaculture industry with new approaches to its sustainability (28). To the best of our knowledge, few studies have used similar approaches to assess the potential of the fish protein-based adaptations as reliable signatures of stress, in contrast with the common and untrustworthy stress indicators i.e., cortisol, glucose and lactate. In this study we have shown that farmed gilthead seabream exposed to different chronic stressors exhibited a panoply of physiological responses, therefore, supporting the accepted idea that the commonly used indicators are far from the robust to assess the level of fish stress. Fish respond to acute and chronic challenges by the activation of a coordinated series of neuroendocrine pathways, in an attempt to regain homeostasis, but when this stability is surpassed, detrimental effects at the whole-animal level are visible (14). These responses or effects are commonly used as indicators of compromised welfare and are rather late indicators and usually a point of no return. By this reason, cortisol has been the most commonly used physiological indicator of the primary response to stress (21). However, there is a shred of evidence indicating that this corticosteroid is not a reliable biomarker of long-term stress (27,35–37), and our results are in agreement with this allegation. Gilthead seabream exposed to high stocking densities (45 kg/m3) during the 54 days of the OC trial showed to reconfigure the cortisol secretion by expressing a tendency decreased response of this metabolite as compared to unstressed fish. Such consequence is suggestive of habituation to the stressor whereas a potential negative feedback has occurred when unstressed fish become stressed (i.e., consequence of the sampling on fish that did not have to adjust physiologically to previous environmental challenges). The same outcome was observed in juvenile gilthead seabream confined for 14 days at 26 kg/m3 (38) and in meagre cultured at different stocking densities for 40 days (39). In the NET trial, contrarily, plasma cortisol levels were significantly higher in handled fish, suggesting that gilthead seabream was not able to adapt to the handling stressor and it might indicate a mechanism of sensitization of the HPI axis. The fact that it was not constant (like overcrowding) but repeated instead, summed to its unpredictability and severity, could have prevented the possibility of accommodation to stress. Within this same context, under hypoxia conditions, was observed a chronic down-regulation of the cortisol response and here, either habituation or exhaustion might have occurred by overstimulation of the HPI axis. Exhaustion is however speculative since differences in behaviour were not monitored despite the reduction of activity seen alongside with the reduction of oxygen. Other explanation could be the fact that the decrease in food intake by prolonged hypoxia, reduces energy needs and thus oxygen demand, therefore lowering cortisol levels. Overall, the aforementioned observations suggest that the cortisol response and the capacity of adaptation were modulated by the nature of the stressor.
A prolonged and constant stressor might, therefore, induce a different physiological reaction compared to repeated acute stress. Studies report contradictory results regarding whether fish adapt similarly to prolonged chronic stressors and repeated acute stressors, but factors like species, age, sex and individual coping mechanisms seem to be ubiquitous and impact their adaptive processes (24,40,41). However, the process of habituation, characterized by a progressively weaker stress response after a prolonged stimulus, was already suggested and demonstrated in other studies as a result of the acclimation of the hypothalamus-pituitary-interrenal (HPI) axis (37,42), a mechanism that is not completely understood so far. High individual variability was also found in every trial, for both stressed and undisturbed fish, most likely due to individual differences in the stress response related to intrinsic factors of the animal, namely, their coping styles, the cognitive perception of the stressor and the value given to such challenge (38,43,44). Additionally, values registered for control fish, in every trial, are higher than the reference values reported in the literature for this species (45). These discrepancies can have several causes and that is why cortisol should be used with caution when evaluating the magnitude of the stress response. The difficulty of measuring the resting levels of this hormone might also be one of these causes, since the lack of proper planning when sampling cortisol, or the manipulation needed to net and anaesthetize the fish, can result in high “control” cortisol levels that do not correspond to the “genuine” basal levels i.e., the non-manipulated fish levels. Also, it is well established that following the perception of an acute stressor, the levels of circulating stress markers increase within the first minutes or hours of stress response, returning to basal levels while time elapses, usually within 24 h (17,46,47).
Secondary physiological responses are characterized by an increase in glucose and lactate levels in blood plasma to satisfy the increased energy expenditure. Glucose release is primarily mediated by catecholamine-induced glycogenolysis and maintained by cortisol-induced hepatic gluconeogenesis (13). Thus, changes in glucose usually follow similar trends than cortisol after the stressor (14), which is corroborated by the levels of plasma glucose registered for all the three trials (Fig. 1). In this study, glucose levels, besides following the same trend as cortisol levels, are, in general, below the basal values for this species (45). This could be related to the fish’s inability to maintain the same levels of glucose in the blood due to the high demand for glucose mobilization to other tissues. The decrease of plasma glucose levels in OC is consonant with the decrease in the cortisol levels, supporting the hypothesis of habituation, or by exhaustion of the endocrine system to the chronic stressor and depletion of energy reserves (27). The significant increases in the plasma glucose levels of stressed fish from NET and HYP trials are consistent with previous studies, which show that glucose rises during air exposure or low oxygen levels, due to stimulation of muscle glycogenolysis and hepatic gluconeogenesis, where glucose is synthesized to maintain the energetic substrates’ demand (48). Similarly to cortisol, glucose and lactate circulating levels also return to basal levels within hours post-stressor, which also make of these metabolites weak reliable markers in case of prolonged stressors (49,50). Additionally, studies also demonstrate that glucose variations in the blood are not only hormonal-induced due to stressful practices. Factors like variations in the water temperature and pH, anaesthesia, diet composition or fasting can also affect plasma glucose levels (51,52).
When tissue oxygen demands are exceeded, in the case of exhaustive exercise, air exposure or hypoxia, for instance, the response to a stressor can become anaerobic and trigger instant glycolysis in order to produce ATP, from stored glycogen, to meet cellular requirements, which consequently leads to lactate accumulation in the muscle (19,53). In this study, changes in circulating lactate are not in agreement with cortisol and glucose variations. Statistically significant differences in the lactate levels were only observed in the OC trial. In this case, if the cortisol response is indeed lower due to HPI-axis acclimation, as suggested before, the lactate recycling rate in the hepatic glycogenolysis is reduced, explaining the significant plasma lactate increase in stressed fish. Additionally, previous studies show that during hypoxia and intense swimming activity, fish produces lactate in the muscle at a higher rate than it can be processed by other tissues (53).
Post-mortem muscle pH and rigor mortis have been used as tissue indicators of ante-mortem stress in numerous fish species (54–56). Stress at slaughter is the main contributor leading to an acceleration of the adenosine triphosphate (ATP) degradation rate before and after cessation of the blood circulation and oxygen supply, consequently influencing the post-mortem muscle pH and the onset and release of rigor mortis (57). The depletion of ATP reserves stimulates the anaerobic glycolysis in the muscle in order to maintain the energy expenditure. This process results in the accumulation of lactic acid, generating H+ ions and consequently lowering muscle pH (57). The magnitude and rate of this pH fall depend on the fish’s energy reserves prior to death, while these energy reserves are influenced by the intensity and duration of the pre-mortem stress. To our knowledge, no studies were performed in this species regarding the effects of long-term chronic stressors on the evolution of post-mortem biochemical processes. Results from this study (Fig. 2) followed the same pH trends as previous studies on gilthead seabream (58,59), however, comparing with the existent studies on pre-slaughter stress (55,60,61), muscle pH values immediately after death are below the ones found in this study, suggesting that slaughter stress was low in our fish. Poli et al 2005 state that in cases of exposure to a chronic stressor for a long time before death, the lactic acid produced can be gradually cleared from the muscle, but simultaneously the energy sources will likewise become gradually exhausted. Hence, when the fish is killed, muscle pH remains higher due to an early end of post-mortem anaerobic glycolysis caused by energy source scarcity. This might explain the significant differences found in the HYP trial, where the highest pH values were observed in the highly stressed fish (HYP15), suggesting that these fish had lower energy reserves. Nevertheless, pH values registered after the 24 HAD, in every treatment, are in agreement with the reported by previous studies in this species at the same sampling times (58,62).
The onset of rigor mortis occurs with ATP depletion. When ATP reaches low levels, actin and myosin in the muscle bind together forming the actomyosin complex and causing stiffness of the fish body (63). Rigor mortis is inextricably correlated with muscle ATP and pH decline. A strong relationship between low muscle pH immediately after death, resultant of increased muscular activity at slaughter, and a rapid onset of the rigor state was demonstrated in a range of fish species (57,60). In this study, the evolution of rigor mortis (Fig. 2) was similar between treatments and significant differences were only found in the NET and HYP trials at 8, and at 8 and 24 HAD, respectively. A delayed onset was observed, starting between 2 and 6 HAD in every trial and reaching the maximum rigor index between 24 and 48 HAD. This delay is consonant with the high muscle pH registered at the time of death, suggesting low energetic reserves. A delayed rigor onset in rested and anaesthetized fish was also shown in previous studies (60). Measuring the glycogen and ATP content in the fish muscle and liver would be a complementary assessment to infer about the energetic reserves a corroborate our hypothesis.
Plasma proteins were evaluated in this study since blood plasma is a very informative biological fluid as it acts as a mirror of the physiological condition of the organism. Also, it is well known that different kinds of stressors result in a variety of immune changes, and proteins with key roles in immunological systems are the main components of blood plasma (42,64). Stress and stress-related hormones are recognized as modulators of the fish immune system (65), however, reactions depend on the intensity and duration of the stressor. The innate immune system is a fundamental defence mechanism in fish for protection against potentially harmful situations (66). The acute phase response is part of this system and it is mainly regulated by cytokines and glucocorticoids (67). This response is characterized by the release of acute-phase proteins (APP), by the hepatocytes, into circulation (68). APP can be classified as “positive” or “negative” depending on their plasma concentration increases or decreases during an inflammatory event (69). The response profile of our fish demonstrated the same tendency of protein changes.
This proteomics-based study of protein changes in the plasma after exposure of farmed gilthead seabream to multiple chronic stressors indicated that the fish’s immune system was affected mainly by net handling and hypoxia stressors. Nevertheless, net handling was shown to be the most impacting. The levels of 20 different plasma proteins (distributed by 56 significantly differential spots), all related with immunological processes, were shown to be modulated by repetitive net handling, comparing to 2 proteins modulated by hypoxia. As mentioned, the same proteins were often detected from different spots on the 2D gels. Such a phenomenon can be due to existent isoforms or caused by adaptive changes of the proteome in an attempt to maintain cellular homeostasis under stress. This adaptation, in addition to protein abundance, may involve changes at the protein degradation level, localization, function and activity – all of which can be modulated by post-translational modifications (PTMs) (70). PTMs can regulate fundamental biochemical processes and be more energetically efficient than altering protein abundance, constituting potential interesting signatures of stress. Studies on PTMs in fish are still scarce.
The protein changes detected (listed in additional file 2), along with the network and GO enrichment analyses (Fig. 5) performed, confirmed the involvement of several components of the innate immune system in the physiological adaptation to these stressors. Proteins considered to be “positive” APP were likewise shown to be increased in abundance in the plasma of fish stressed by net handling (fibrinogen alpha-chain, complement component C3, haptoglobin, complement factor B, warm-temperature acclimation 65 kDa protein, alpha-1-antitrypsin), while proteins considered as “negative” were decreased (transferrin, inter-alpha-trypsin inhibitor, apolipoprotein A-I) (71). A diverse number of proteins involved in the APR was also previously found to be modulated in chronically stressed gilthead seabream (72).
Apolipoprotein A-I (Apo-AI) was only modulated by net handling stress and 17 proteoforms were identified in the plasma proteome map, being mostly decreased in abundance. Apo-AI is the main protein constituent of the high-density lipoprotein (HDL), playing a role in lipid metabolism and participating in the reverse transport of cholesterol from tissues to the liver (73,74). Apo-AI was also found to be decreased in abundance in crowded Atlantic salmon (75). In cod (Gadus morhua) it acted as a negative regulator of the complement system (76). Other two apolipoproteins were also found to be down-regulated in the plasma of fish from NET2 and NET4 groups (Apolipoprotein Eb and apolipoprotein B-100).
The complement system is an essential part of the innate immune system which can be activated through three pathways: the classical, alternative and lectin pathways (77). Fish display a plethora of complement components, mainly complement component C3 (C3), which may present around five proteoforms in a single species (78). C3 is one of the most abundant proteins in the plasma and plays a central role in the innate immune system, supporting the activation of all three pathways (77). In this study, C3, identified in 5 proteoforms, and complement factor B (Bf), identified in 4, were found to be increased in abundance by net handling. Contrarily, C3 was down-regulated in fish exposed to low oxygen levels. Bf also plays a role in complement activation by acting as the catalytic subunit of C3 convertase, an enzyme responsible for the proteolytic cleavage of C3, in the classical and alternative pathways (77).
Several metal-binding proteins, existent in the plasma of vertebrates, can chelate iron, zinc and copper, which are essential elements for the virulence of bacteria (79). Alpha-2-macroglobulin (A2M) is a multifunctional protein (80) found to be down-regulated in the plasma of fish submitted to handling stress. It is mostly known to act as a broad range serine proteinase inhibitor and to bind metal ions (79). Contrarily, haptoglobin, which is also responsible for the sequestration of iron by binding to hemoglobin, was found to be increased in the plasma of handled fish. Similarly, warm-temperature acclimation-related 65 kDa protein (Wap65), which is involved in the scavenging of free heme (81), was increased in abundance by net handling and hypoxia stressors. Wap65 in fish is the homologue of mammalian hemopexin (82) and in most teleosts presents two proteoforms (83). In this study, two spots were also matched to this protein suggesting the presence of these two proteoforms. Transferrin (Tf) decreased in abundance in the plasma of fish stressed by net handling. Tf is a plasma protein also capable of binding iron and an important constituent of the iron homeostasis (33).
In fish, antiproteases are important participants of the non-specific humoral immune defence mechanism (71). A2M is an important factor in this mechanism. Alpha-1-antitrypsin is a serine protease inhibitor, up-regulated in net-handled fish, which is responsible to negatively regulate blood clotting molecules to prevent thrombosis (84). Inter-alpha-trypsin inhibitor H3 is also a serine protease inhibitor, which was found to be down-regulated in the plasma of fish from NET groups. Same protein changes were verified for fetuin-B, a cysteine proteinase inhibitor recently described in teleosts (85). Fibrinogen alpha-chain, a beta-globulin involved in blood clotting, an integral part of innate immunity (84), was found to be up-regulated in the plasma of fish belonging to NET groups.