1. Introduction
In 21st century, Atlantic salmon (
Salmo salar) has become one of the most intensively farmed finfish in the world and the industry continues to expand as consumer demand and global population grow [
1]. Unlike many intensively farmed fin-fish, Atlantic salmon is a carnivore and has a high dietary protein demand [
2]. Sourcing protein for inclusion in salmon aquafeeds has become a major challenge and bottleneck for the aquaculture sector [
3]. Historically fish meal dominated the aquafeed market, but increasing and unsustainable pressure on wild-capture fisheries and rising prices have resulted in the commercial uptake of soy protein sources [
3,
4]. Soybean meal (SBM) gained initial popularity because it is highly available, cheap and has a favourable amino acid profile for farmed fish [
5]. However, for Atlantic salmon, SBM can have adverse consequences to gut health, and triggers soybean meal induced enteritis (SBMIE) attributed to the presence of saponins, non-starch polysaccharides (NSP) and other anti-nutritional factors (ANFs) [
6,
7,
8,
9,
10]. These ANFs can adversely alter the distal gut morphology of salmon [
11,
12], reduce nutrient absorption and efficient digestion [
13], adversely alter the gut microbiome community [
14] and consequently compromises the immune resilience and welfare of fish [
15,
16]. To address this problem, the salmon industry now primarily uses soy protein concentrate (SPC) an alcohol extracted soy product that reduces the inflammatory effect of SBM but has still been observed to have negative impact on gut health and gut microbiome of Atlantic salmon [
17,
18,
19]. Applications of such treatments increase the cost of the aquafeed sector, thus there is a need to ensure tangible benefits for salmon production and welfare [
20].
The gut microbiome of fish has been linked to key traditional measures of aquaculture productivity, including, nutrient digestibility and availability, growth performance [
21], the metabolism, immune development, and disease resistance [
22,
23,
24], the metabolism [
25,
26,
27]. The composition of the gut microbiome is in turn influenced by the feed, host´s environment [
28], and genetics [
29,
30]. In juvenile salmonids, the gut microbiome is particularly malleable to alterations in dietary protein sources [
31,
32] and this can influence the establishment and development of the gut microbiome into the adult developmental stage [
33]. When salmonids are fed diets containing SBM or SPC, the composition of the gut microbiome is distinct from fish fed marine diets. There is a greater presence of bacteria associated with inflammation, and an imbalance in the community which has been linked to poor intestinal health [
17,
18]. Increasingly, aquaculture research is highlighting the potential of dietary treatments that can have functional benefits to the host by altering the gut microbiome. Such treatments include the application of probiotics, prebiotics, synbiotics, fermentation and enzyme-treatment of ingredients [
34,
35,
36,
37,
38]. Prebiotics are non-digestible fibers that are not directly used by the host but act as a nutritional source to support the proliferation of desirable gut microbiota [
39], they are widely used in human foods, terrestrial farming and have high acceptance for application to aquaculture diets [
40]. In aquaculture, fructose oligosaccharides (FOS), and mannan oligosaccharides (MOS) have been used in feeds to provide gut microbiome-enhancing substrates and even growth performance [
41,
42]. Enzyme treatment has also become a popular processing method to improve the value of proteins for animal and aquaculture feeds [
35]. Enzyme treatment of plant proteins such as SBM has been used to breakdown long-chain carbohydrates such as NSPs to improve their nutritional value for dietary inclusion [
37,
43]. This has a potential added benefit, that the breakdown product of long-chain carbohydrates includes shorter chain oligosaccharides, which in fish diets are not utilised by the host but may have a prebiotic value to the gut microbiome of the host [
44].
There is a growing understanding of the role of gut microbiome not just in traditional measures of productivity for aquaculture but also in the observation of the bi-directional link between the gut microbiome and the central nervous system of animals including fish, which has been dubbed, the gut-brain-axis [
45,
46]. In fish there is evidence that the gut microbiome influences swimming behaviour [
47], feeding behaviour [
21] and social behaviour [
48]. Behaviour is an early indicator of changes in environment [
49]. Thus, monitoring behaviour is a valuable tool for managing welfare [
50], a growing concern for the aquaculture sector [
51,
52].
In fish, five behavioural types relevant for welfare in aquaculture have been identified along five major axes: boldness-shyness, exploration-avoidance, sociability, activity, and aggression [
53,
54]. When such behavioural differences are consistent over time and context, they can be referred to as personality [
55,
56]. Boldness is a measure of predisposition to take risks [
57] and individuals can be categorised on a continuum from bold to shy. Exploration is a measure of predisposition to engage with a novel environment or object [
49] and individuals can be categorised on a continuum from active explorers to cautious explorers [
58]. Coping style, closely related to personality, refers to a coherent set of behavioural and physiological stress responses, which is consistent over time and which is characteristic to a certain group of individuals [
55]. This approach allows to classify the individuals of a given population into reactive and proactive categories. Proactive fish are bold individuals taking more risks and exploring their environment faster (less cautiously) when exposed to novelty [
59,
60] and they are more aggressive, dominant [
60,
61], and display low variability and flexibility in their behavioural responses with environmental changes [
62]. By contrast, reactive fish are shyer individuals tending to be risk-averse and are generally neophobic, show higher behavioural flexibility and are more responsive to their environment [
63], and are more sociable than proactive individuals [
64]. If there is a correlation between different behaviour traits at the individual level, this is known as the behavioural syndrome which can be associated with welfare and freedom to express natural behaviour [
51,
54,
65]. Behavioural traits can be altered by dietary ingredients making them important indicators to monitor [
15,
47,
66]. In juvenile Senegalese sole (
Solea senegalensis), replacing fish oil with vegetable oils resulted in significantly more reactive individuals [
67]. It will be important to understand how dietary treatments affect the coping style and stress response of juvenile Atlantic salmon since it could have impacts on the productivity, fitness and welfare of farmed fish [
53].
The objective of this study was to optimise the application of SBM in farmed Atlantic salmon diets and to elucidate the impact of novel enhancements of SBM (enzyme pre-treatment and addition of prebiotic respectively) on growth performance, gut microbiome, and behaviour. This study hypothesises that these novel enhancements will benefit salmonid aquaculture by supporting good growth performance, increasing the presence of desirable bacteria in the gastro-intestinal tract and improve fish welfare compared with unenhanced SBM. This study addresses important knowledge gaps for both the aquaculture industry and the field of applied microbiology in a commercially relevant model organism.
2. Materials and Methods
2.1. Experimental animals and study design
Atlantic salmon (Salmo salar) hatched by Stofnfiskur Ltd. (Vogar, Iceland) were incubated at 5.5°C and eyed eggs were transferred to Laxar ehf. (Kópavogur, Iceland). Fish were raised to first feeding using standard commercial techniques and start-feed diet BioMar Inicio-plus (United Kingdom) of 0.5mm pellet size and 12°C water temperature. Fry were transferred to Verid aquaculture station of Hólar University (Saudárkrókur, Iceland) where they were acclimated for 6 weeks prior to the start of the experiment. All fish within the experiment were individually weighed and measured under anaesthetic (2-phenyoxyethanol of 300ppm) following a 24-hour fasting period. Fish were split into 12 identical 20L-White cylindrical PVC tanks, in triplicate for each feed treatment. Each tank contained 40 fish with initial weight of 2.2±0.4g. Fish were maintained at 11±1.9°C under continuous light of 250±50 lux and 90±10% oxygen saturation. Fish were fed with experimental feed treatments for 70 days. The experiment was performed following Icelandic guidelines and within the permits and licenses of Verid aquaculture station.
2.2. Experimental feeds and feeding
A fish meal-based control (FM) diet, an untreated soybean meal (US) diet, an enzyme treated soybean meal (ETS), and an untreated soybean meal with added fructose oligosaccharide (USP) were used. The chemical composition of protein sources used in this study is presented in
Table 1. For each of the soy diets, US, ETS and USP were included at 25%, replacing fish meal, the diet formulation and nutritional composition of each diet are presented in
Table 2. The enzyme treated soybean for the ETS diet was produced by treatment with a 50:50 blend of Hostazyme X (main activity: Endo – 1,4 – β-xylanase, side enzyme activity: End – 1,4 -β-glucanase (cellulase), endo – 1,3(4) – β-glucanase, α-amylase, protease) and Hostazyme C (main activity; End-1,4 – β-glucanase, side enzyme activity: Endo – 1,4 – β-xylanase, endo – 1,3(4) – β-glucanase, α-amylase). The enzyme preparation was blended in 50°C water. The solution was added to 2L of water and 1kg of soybean meal and mixed for 15 minutes. The mix was then incubated for 3 hours at 50°C and manually mixed every 30 minutes. Following incubation, excess water was removed through pressing. The material was then further dried for 80°C. The diets were produced by cold pelletisation at Matís ohf. (Iceland). All dry ingredients were milled to homogenise particle size (IPHARMACHINE, Germany). Dry ingredients were homogenised in a standard food mixer (KitchenAid, USA) and the mix was milled again to improve the homogeneity of the feed. The dry mix was returned to the food mixer and fish oil was added and 200ml of water to produce the optimal consistency for processing in a pasta machine set to 0.5mm strings (ADE, Germany). Strings were dried in a commercial food dryer (Kreuzmayr, Austria) to <10% moisture content. During the 70-day feeding trial, fish were fed continuously by electric belt-feed with identical feed volumes at 15% excess fed based on the feed requirements for this developmental stage.
2.3. Growth performance
The fish were not fed for 24hours prior to measuring. After 70-days of feeding, all individual fish from each replicate tank within each treatment were anesthetised (2-phenoxyethanol 300ppm) and their wet weight (g) and total length (cm) measured. The specific growth rate (SGR) (%) over the study period was calculated: SGR = ((Ln(Final Weight)-Ln(Initial Weight)) x 100)/t, where t is the number of days of the experimental period. Mortality was monitored daily during the trial period.
2.4. Gut sampling
Following the final measurements of growth performance, all fish were left for one week to recover where they were maintained on their respective experimental diets. Fish were then fasted for 12-hours. Three fish per tank n = 9 per treatment were randomly selected for gut microbiome analysis. Fish were humanely euthanised with a lethal dose of anaesthetic (2-phenoxyethanol 600ppm). The outside of the fish was washed in 90% ethanol followed by sterile distilled water. The gastro-intestinal (GI) tract and contents (from the start of the mid-gut, just below the pyloric caeca to the end of the distal gut) was sampled under sterile conditions and directly frozen at -80°C.
2.5. DNA extraction, PCR amplification and sequencing
GI tract samples were individually manually homogenised in a sterile petri dish with a sterile scalpel to physically break up the material. Sample material was transferred to a sterile 2ml Eppendorf tube with 300µl of sterile 1mm diameter silica beads (BioSpec Products, United States) and 800µl of CD1 solution from the QIAamp PowerFecal Pro DNA kit (QIAGEN, Germany). Samples were vortexed for 5 seconds and shaken at maximum speed (30Hz) in a laboratory mixer mill (Retsch MM400) for 1 minute. The supernatant (~800µl) was transferred to the PowerBead Pro tube from the QIAGEN QIAamp PowerFecal Pro DNA kit. The protocol for this DNA extraction kit was then followed and final DNA was eluted with 80µl of C6 solution. A DNA negative (no material added) was also run to ensure no contamination occurred during DNA extraction. DNA concentration was measured for each sample with 2µl with Invitrogen Qubit dsDNA BR Assay kit (Invitrogen, Carlsbad, CA, USA). DNA was diluted to 4ng µl-1 in 50µl aliquots. Samples were subjected to PCR of the V3-V4 regions of the 16S rRNA gene with a universal bacterial primer pair S-D-Bact-0341-b-S-17 (5′−CCTACGGGNGGCWGCAG-3′)/S-D-Bact-0785-a-A-21(5′-GACT- ACHVGGGTATCTAATCC-3′) (Klindworth et al., 2013). PCR mater mix included diluted DNA, nuclease-free water, Q5 High Fidelity DNA polymerase (New England Biolabs, Ipswich, USA), Q5 GC Enhancer, 0.5µM of each primer containing Illumina overhand adapters, and 1 x Q5 Reaction buffer, 200µM dNTPs (New England Biolabs, Ipswich, USA). Both positive and negative samples were also run in the PCR to monitor for absence of contamination and successful amplification. The thermocycling protocol had an initial denaturation step (90°C for 30s), followed by 35 cycles of: denaturation (90°C for 10s), annealing (52°C for 30s), and extension step (72°C for 30s), with a final extension (72°C for 2 min). Libraries were multiplexed with Nextera XT v2 barcodes (Illumina, USA), normalised using Sequel-Prep normalisation places (ThermoFisher Scientific, USA) and sequenced on MiSeq desktop sequencer (Illumina, USA) using v3 chemistry and 2 x 300 cycles.
2.6. Behaviour
Following gut microbiome sampling, eight (8) fish from the remaining tank population were randomly selected to be uniquely tagged with Visible Implant Elastomer (VIE) colour tags (Leblanc and Noakes, 2012) under anaesthetic (2-phenoxyethanol 300ppm). Each fish was injected parallel to the dorsal fin, below the surface skin layer with two 1cm strips of colour with unique combinations so each of the 8 selected fish per tank (n=24 per treatment) could be identified and followed through all behavioural observations. All remaining untagged fish per tank were retained in tank to maintain a consistent stocking density but were not included in behavioural observations. Each tagged fish was observed in two different behaviour contexts (a swimming test and an open field test with shelter) and each of these tests was performed two times. There was always a one-week recovery period between each behavioural observation. Behaviour observations were performed in a dedicated room to minimise external disturbance. Fish were individually tested in all behavioural observations. Behaviour was recorded using a monochrome camera (Basler Ace acA1920-150um, Germany) with a frame rate of 30Hz and resolution of 1280 x 1024 pixels placed 112cm (swimming test arena) and 110cm (open-field test) above the water level of the respective behavioural arena. Video recordings were analysed with EthoVision XT software (Noldus, The Netherlands), which was used to track fish in each arena. The VIE tags were recorded for each fish after each behaviour observation so their individual data could be tracked across the different time and contexts of the behavioural observations.
2.6.1. Swimming test for assessing response to stress and exploration.
A swimming test was used to collect behaviour observation data for both stress response and exploration-avoidance axis. The apparatus for assessing these traits consisted of 4 circular arenas (diameter = 25cm, water depth 7.8cm (4L), height= 15cm). The arena was illuminated from below to provide uniform light intensity (260 lux) (
Figure 1). Tagged fish from each tank were collected from the trial tank and transferred in identical white transfer buckets with closed tops and then placed into individual arenas in random order. All four arenas were filmed simultaneously and continuously for 20 minutes. For the purposes of this analysis, this was separated into four virtual periods. The 0-5 minutes was the acclimation period (AC), 5-10 minutes was the normal swimming period (NS). At the 10-minute mark a stress was elicited (bottom lighting was switched off for 3 seconds and then switched back on again), minutes 10-15 was thus the post light stress period (PLS), and 15-20 minutes was the recovery period (RC). Behaviour was not recorded while the light was off.
Variables collected through the swimming tests were, mean distance from centre-point (cm): the distance of the centre-point of the fish body from arena centre (DisCent). Total distance moved (cm): the distance travelled by a fish measured from the centre-point of body between consecutive X-Y coordinates (TotDis). Mean velocity (BL s-1): the distance moved by the centre-point of body per unit time between consecutive X-Y coordinates normalised to body lengths per second (Vel). Absolute angular velocity (° s-1): expressed in degrees per second was calculated by EthoVision XT software. Vangn = RTAn/tn – (tn-1), where RTAn represents the relative turn angle for sample n, and tn-(tn-1) is the time difference between the current and previous sample (AngVel).). Mobility state, the cumulative duration for which the fish body is changing, highly mobile (s) when cumulative duration is 60% (HiMob). Moderate mobile (s) when cumulative duration was between 20-60% (MedMob). Immobile (s) when cumulative duration was below 20% (Immob).
2.6.2. Open-field test for assessing boldness
An open-field test (OFT) with a shelter was used to assess boldness [
68,
69,
70]. The apparatus for assessing this trait consisted of 4 rectangular arenas (40 x 30 x 25 cm), water depth 6cm (7L) with an identical shelter placed in the bottom-right of each arena (14 x 6.5 x 6.5 cm). The arenas were illuminated from below to provide uniform light intensity (260 lux) (
Figure 2A). Tagged fish from each tank were collected from the trial tank and transferred in identical white transfer buckets with closed tops and then placed in random order into individual arena shelters through a top compartment (4cm diameter) which was then closed, and the main shelter door (a sliding opaque trapdoor) was kept shut, closing off access to the arena. All four arenas were filmed simultaneously and continuously for 25 minutes. The first 5 minutes the fish remained shut inside the shelter acting as an acclimation period, after which the door to each arena shelter was lifted simultaneously and filming continued for a further 20 minutes uninterrupted. Each individual arena was virtually divided into four zones using the EthoVision XT software. These zones were called Shelter, Entry, Border and Centre (
Figure 2B). The centre zone is considered high-risk and staying close to the border edges of a space is considered more cautious and an indicator of a shyer individual [
57,
71]. Variables used to characterise OFT behaviour for the trait of boldness were as follows: Latency (s) to emerge-time taken to exit the shelter (Lat). Time spent in each zone (s) (Shelter, Entry, Border, Centre) respectively (Shelt, Ent, Bord, Cent). Mean distance from shelter (cm) (DisShelt). Number of returns to shelter (Ret). Cumulative duration spent highly mobile (s) (HiMob). Distance moved: the distance travelled by a fish measured from the centre-point of body between consecutive X-Y coordinates. Mean velocity: the distance moved by the centre-point of body per unit time between consecutive X-Y coordinates normalised to body lengths per second. Absolute angular velocity: expressed in degrees per second.
2.7. Statistical methods
Statistical analyses were performed in Rstudio version 3.6.1 (2019-07-05). Tests for growth performance, gut microbiome and behavioural characteristics were two-tailed with a significance level set to α = 0.05. We fitted multivariate linear mixed models to estimate the repeatabilities and the correlations between the personality estimates using Bayesian statistical modelling.
2.7.1. Growth performance
SGR (%) was used to assess growth performance across the different dietary treatments in this study. A Linear Mixed Model (LMM) with the package lme4 [
72] was selected defining feed treatment (FM, US, ETS, USP) as a fixed factor and Tank as a random nested factor of feed treatment. Since the random nested factor of tank did not cause significant variation in SGR (%), a simplified statistical test was adopted. A Linear Model (LM) with the package nmle [
73] was used, where feed treatment was a fixed factor and the model residuals had a normal distribution. A Tukey post-hoc test was applied to assess pairwise differences between feed treatments.
2.7.2. Gut microbiome
To assess the gut microbiome of fish fed different dietary treatments in this study, demultiplexed FASTQ files from Illumina were processed to produce amplicon sequence variants (ASVs) with DADA2 package version 1.16.0 [
74]. The filterAndTrim variables used were as follows: truncLen=c(280,250), trimLeft=21, maxN=0, maxEE=c(2,2), truncQ=2. Taxonomy was assigned to ASVs from version 138 of the SILVA database and the function assignTaxonomy [
75]. The R packages, phyloseq [
76], microbiome [
77] and vegan [
78] were used to analyse the microbiome community and ggplot2 was sued to visualise key data [
79]. The average read number output from the DADA2 pipeline was 20643±12702, and two samples were removed as they did not contain any reads after processing. Five PCR negative samples were also sequenced to control for any contamination that may have occurred during sample amplification. These controls were used to remove suspected contamination from samples using the decontam package with prevalence method and threshold of 0.5 [
80]. Read depth was normalised across all samples using the function rarefy_even_depth. Raw 16S rRNA gene amplicon reads can be found in the Sequence Read Archive following submission.
Alpha and beta diversity indices were used to quantitatively analyse the gut microbiome. Alpha diversity measures selected were, observed richness of AVSs, Shannon diversity, Chao1 diversity and Pielou´s evenness. A GLMM was used to assess if there was a significant difference in the alpha diversity measures between the different feed treatments in this study. In this model, dietary treatment was a fixed factor and tank was a random nested factor of feed treatment (which was tested by a Likelihood Ratio Test [
81]). A Tukey test was applied for post-hoc testing. Beta diversity was assessed by transforming microbiome community data using a Bray-Curtis dissimilarity matrix. An Analysis of Similarity (ANOSIM) was performed to determine if there was any significant difference between and within the gut microbiome community of fish fed the different diets. Relative abundance of taxa at the phylum level and genera level as a proportion was calculated. Taxa present in the community at <1% relative abundance were grouped into a category labelled “Other”.
2.7.3. Behavioural characteristics
To assess stress response, all swimming test variables were analysed using Linear mixed effect models (LMM). For each analysis, the explanatory variables included in the full model were, feed treatment (FM, US, ETS, USP), period (AC, NS, PLS, RC), and TL, and interactions between feed treatment and period were considered. Random effects considered in the model were trial replicate number (1, 2), and Tank (1,2,3) nested in feed treatment, and individual ID.
To assess the behaviour trait of exploration, TotDis recorded for AC and NS only, was used as a proxy [
82]. Both traits i.e., boldness and exploration were then compared between feed treatments in two respective LMM models where boldness score or TotDis were the response variables. The full models were reduced by backward selection based on the Akaike Information Criterion (AIC) [
83]. Diagnostics based on residuals of the model were performed to assess the adequacy of the reduced model and compliance with the underlying assumptions. Dependent variables were transformed whenever necessary to ensure that the residuals followed the assumed error distribution. Finally, the effects of the independent variables were estimated from the reduced models and their significance was tested by likelihood ratio tests (LRT) between models respecting marginality of the effects that are supposed to follow a χ2 distribution under the null hypothesis (type II tests; [
81]). This analysis was followed by a post-hoc multiple comparison test [
84] to assess pairwise differences between models.
2.7.4. Repeatability of behaviour traits and correlations
A Bayesian multivariate linear model using Stan [
85] was run using the ‘brms’ package [
86]. We used one model for each feed treatment (four models in total). Each model simultaneously regressed each dependent variable (i.e., boldness score and TotDis) against a set of fixed factors (Tank, Trial, SGR and TL). The random-effects structure included individual fish identity (ID) as a grouping variable, allowing us to calculate the repeatability of Boldness and TotDis as the ratio of the among individual variance and the sum of the among individual and residual-level variances (i.e., personality) [
87]. Fish that did not exit the shelter for at least one of two the trial were removed from the analysis. Moreover, the model estimated covariances between Boldness and TotDis at both the ID and residual levels. The among individual covariance quantified the degree to which Boldness and TotDis was correlated among individuals across multiple trials (i.e., behavioural syndrome). The model was run for 64000 iterations (32000 for warmup and 32000 for sampling), four chains, adapt delta was set to 0.9, the max tree depth to 40, and All other parameters were set to their defaults. Convergence was assessed using the standard diagnostics provided by Stan, including the potential scale reduction factor (̂R), effective sample size, and visual inspection of trace plots and histograms for each model parameter. Unless otherwise noted, we used posterior modes for point estimates and higher posterior density with 95 % coverage for uncertainty intervals (UI95 %), respectively calculated using the map_estimate and hdi functions from the “coda” package [
88].
4. Discussion
The objective of the present multi-disciplinary study was to optimise the application of soybean meal (SBM) in formulated aquafeeds for salmonids. This study reports the impact of two novel enhancements of SBM, an enzyme pre-treated SBM (ETS), and a SBM with addition of the prebiotic (fructose oligosaccharides, FOS) (USP) compared with an untreated SBM (US) and a fish meal control (FM) on the growth performance, gut microbiome, and individual behaviour traits relevant to the culture of juvenile Atlantic salmon (Salmo salar).
The growth performance of Atlantic salmon in this study across all treatments was within a normal range for the freshwater stage at the experiment water temperature [
89,
90]. The findings were consistent with the existing literature, showing that FM diets support higher growth rates than untreated SBM for Atlantic salmon, both at inclusion levels of 20% SBM [
12], and at higher levels of 31% SBM [
91]. However, salmon parr fed diets containing 16.7% SMB [
92] or 40% SBM [
93] performed comparably to those fed a FM diets. In seawater phase, inclusion levels of 20% SBM and higher show reduced growth rate in Atlantic salmon compared to those fed FM diets [
11], although inconsistencies in the literature exist [
94]. The two novel enhancements of SBM applied in this study supported comparable growth performance to FM, suggesting both enzyme pre-treatment and addition of FOS have potential benefits for improving SBM application for growth performance of salmon during the freshwater stage. Similar improvements have also been seen with SBM treated to reduce the content of oligosaccharides and ANFs compared with an untreated SBM at a 40% crude protein supply [
95]. In Japanese seabass (Lateolabrax japonicus) NSP enzymes applied to feed also had beneficial effects on growth performance [
96] and in the white-spotted snapper (Lutjanus stellatus) NSP targeted enzyme treatment of Gracilaria lemaneiformis growth performance was also enhanced [
97]. On the contrary, [
98] found no growth performance benefits of NSP-targeted enzyme treatment of SBM for Atlantic salmon. Furthermore, phytase pre-treatment on SPC replacing 60% of FM in SW Atlantic salmon did not improve growth performance [
98]. The addition of FOS to diets improved growth performance of Rainbow trout (Oncorhynchus mykiss) [
99], however, when FOS was added to FM diets for Atlantic salmon, no effect on growth performance was observed [
42].
The differences in the gut microbiome community assemblage between fish fed different diets have been observed in many studies with salmonids, between FM and SBM [
30,
32] and between different treatments of conventional protein sources [
18]. Other studies have also found FM based diets to support greater diversity and community evenness in fish gut microbiome compared with other protein sources [
100,
101]. The findings of this study that USP fed fish had comparable diversity to US, and that ETS fed fish had lower diversity than all other treatments, were similar to what was found in white sea bream (Diplodus sargus) fed FOS [
102]. However, research has shown increased gut microbiome diversity with the additions of prebiotics for juvenile hybrid Tilapia (Oreochromis niloticus♀ × Oreochromis aureus ♂) [
103]. The lack of difference in the present study may be driven by the need for longer term continuous administration to detect positive gut health impacts of feed additives in fish [
104]. It has been noted that diets with SBM which contains natural oligosaccharide sources may mask any beneficial effects of additional prebiotics [
102,
105]. The impact of low values for diversity and evenness found for ETS fed fish will be important to elucidate, since high gut microbiome community evenness and diversity values have been associated with improved productivity in aquaculture [
106]. In the present study, it seems that this lower diversity and evenness is driven by a dominance of the lactic acid bacteria (LAB) Enterococcus, a genus that has been associated with growth and immune benefits to fish [
107]. Because the presence of this genus has not excluded other taxa establishing in the gut, reduction in evenness and diversity may not adversely affect fish gut health. Although community dominance at such an early development stage should be monitored to ensure that an undesirable dysbiosis does not establish itself [
21]. This study suggests that the application of FOS has quite different impacts than the enzyme-treatment at the genus taxonomic level, both may be having a prebiotic effect, but stimulate the growth of different communities. The increased presence of the LABs in SBM diets compared to FM is consistent with the existing literature for salmonids [
30], as is the addition of prebiotic ingredients increasing levels of LAB in the common carp (Cyprinus carpio) [
108]. LABs have been associated with improved digestive function, improved gut health, and disease resistance [
38,
109,
110]. Since the levels are lower for these two taxa in the ETS fed fish, it suggests that the enzyme treatment may have broken down part of the SBM favoured by these two taxa, a phenomenon that warrants further investigation. The greater presence of Clostridium_senu_stricto_7, Clostridium_senu_stricto_18 and Peptostreptococcus in FM fed fish have been found in faster growing individual Rainbow trout [
111], which is consistent with the growth performance findings for this trial. These bacteria have been linked to fermentation of different amino acids [
112], suggesting their growth may be facilitated by amino acids present in FM that are not present or less present in SBM. There is a prevalence of Hathewaya and Anaerosalibacter in FM fed fish compared with SBM fed fish in the study, these taxa have been found in Chinook salmon (Oncorhynchus tshawytscha) gut microbiomes associated with presence in feed [
113] suggesting an association with FM. There is a need for future research to address the functional role of LAB for Atlantic salmon and to establish taxa associated with growth performance and gut health benefits to improve the application of feed additives and dietary enhancements.
This study is the first to investigate behavioural traits relevant to welfare related to diets containing different protein sources, and enhancements of dietary protein for Atlantic salmon. The light stress did not appear to produce a response in the fish tested, suggesting the need for a different stress source, as the fish may have been habituated to light stress as they were kept under continuous light [
114]. It is also possible that broodstock selection and process of domestication of Atlantic salmon genetic lines could already have selected for more proactive styles [
62,
115] which may also explain why there was low responsiveness to the elicited stress. While the scores of the traits of exploration and boldness behaviour traits did not differ significantly between the diets in this study, there was an interesting trend that FM fed fish displayed a more proactive coping strategy i.e., higher boldness and exploration values than fish fed any of the SBM diets. For ETS fed fish, there was indication of more proactive-type behaviour compared to US fed fish. This is consistent with results from juvenile Rainbow trout, where fish fed a plant-based diet, exhibit an increase in apathetic behaviour and an increased stress response (both traits of reactive coping styles [
53] compared to a marine diet [
15]. This goes along with other results in juvenile Senegalese sole (Solea senegalensis) individuals fed diets enriched with fish oil from Cod liver, compared with vegetable oils from Linseed, Soybean and Olive, that were more proactive [
116].
For boldness, the repeatability was highest in FM fed fish showing the existence of personality trait which could indicate that they may be more suited with the aquaculture environment. However, regarding the exploration, high repeatability was shown in all treatments, which could indicate lesser impact of the diet treatment on this trait i.e., highly canalised trait [
117]. Existing studies have indicated that more proactive individuals do show a higher degree of repeatability for behaviour traits which matches well with the findings of the present study. Interestingly, there was a positive link between growth and exploration in the FM treatment only which is in accordance with previous studies showing higher fitness in proactive fish e.g., higher reproductive success in gilthead seabream (Sparus aurata) [
67], higher growth rates for salmon [
118], faster feeding recovery after stress [
119] and lesser sensitivity to environmental stress [
120] for Rainbow trout, all of which may make them better suited to the intensive farming environment of salmonid aquaculture [
115] as these individuals may be more productive and also experience better welfare. Repeatability of both traits was also high in the ETS which could indicate a benefit of the enzyme treatment whereas the US (and the USP to a lesser extent) showed no evidence for any pronounced personality trait. This could indicate an impact of these diets on welfare. In summary, high level of repeatability of both traits may show consistent freedom to express natural behaviour which may indicate a good welfare condition [
65,
115].
In the present study a clear behavioural syndrome was detected only for ETS but surprisingly not for FM which may indicate the fish fed ETS experiencing different conditions or pressures than the fish fed other diets [
58]. Domestication of salmon may also reduce the presence of behavioural syndromes. In a study of urban song sparrows no correlation between traits of boldness and aggression could be found, whereas there was a correlation, and thus behaviour syndrome in their wild counterparts [
121].This also may suggest that the behaviour of US fed fish in this study was more impacted by the culture conditions and have adjusted by becoming more behaviourally plastic or showing adapted behaviours to cope [
58]. There is a need for further research attention on the impact of dietary proteins and treatments on behaviour and welfare for farmed Atlantic salmon to support the optimisation of existing feed materials and the industry uptake of alternative proteins. This study also highlights the need for greater understanding on how domestication may alter salmon behaviour and therefore welfare in intensive farms which will help inform selection programs to optimise coping styles in farmed fish [
122,
123].
Although, this study did not specifically address the link between gut-microbiome and behaviour at the individual level, our results suggest potential impacts of diets on the gut microbiome that in turn influence behaviour. Indeed, fish fed SBM diets displayed an altered gut-microbiome, lower boldness repeatability and tended to have lower mean boldness values. Studies have shown that manipulating fish microbiome can modulate the behaviour via the gut-brain axis [
47] and a recent study on farmed Arctic charr (Salvelinus alpinus) [
124] showed evidence of bacterial strains influencing the response to stress and growth where fish fed with E. thailandicus 04-394 and L.brevis ISCAR-07433 displayed an increase in motility and slower growth which can be interpreted as lower stress coping ability. Using personality has been suggested to be a valuable tool to reduce chronic stress in captive fish where bolder individuals are usually more resistant to chronic stress [
125]. Therefore, the gut alteration by SBM diets could directly impact boldness in Atlantic salmon, which in turn could affect reaction to stress and ultimately the growth. In this study, FM fed fish were the treatment group that grew the fastest, had a distinctly different gut microbiome community visible at the taxonomic level and tended toward more proactive behaviour and potentially better welfare compared to any of the SBM treatment. Conversely, the US fish differed the most from the FM fish in all aspects of this multi-discipline approach, showing poorer growth, a lower diversity and evenness of the gut microbiome community and behaviours that indicated differing impact compared with any other diet. While enhancements of SBM ameliorated some of these traits, there were still difference with FM fish, highlighting the importance of such multi-disciplinary studies to elucidate the impact both to traditional productivity measures for aquaculture but also to measures that explore indirect impacts of nutrition on immune development and welfare. Future studies should aim to establish functional links between these different fields to optimise Atlantic salmon aquaculture. This will be especially important as new feed enhancements and alternative proteins enter the market to ensure they are safe and effective for the aquaculture sector and to support optimal broodstock selection programs.
Author Contributions
“Funding Acquisition, Methodology and Conceptualization, A.L, D.B, W.K, and M.Ø.; Project administration A.L.; investigation and formal analysis, A.L, D.B, C.S, S.K, S.S, B.B and G.L.; writing—original draft preparation, A.L, D.B, M.Ø, and C.S.; writing—review and editing, A.L, D.B, MØ; W.P, R.B, and G.L.; supervision, D.B, M.Ø, W.P and R.B.; All authors have read and agreed to the published version of the manuscript.
Figure 1.
Swimming test arenas used for observation of behavioural traits of exploration and stress response.
Figure 1.
Swimming test arenas used for observation of behavioural traits of exploration and stress response.
Figure 2.
(A). Open Field test arenas used for observation of the behavioural trait of boldness. (B). An arena with four virtual zones outlined in the EthoVision XT software, Shelter, Entry, Border, and Centre.
Figure 2.
(A). Open Field test arenas used for observation of the behavioural trait of boldness. (B). An arena with four virtual zones outlined in the EthoVision XT software, Shelter, Entry, Border, and Centre.
Figure 3.
Average Individual Specific Growth Rate (%) for each dietary treatment. Different letters above the bars indicate significant differences (Tukey post-hoc test, p<0.05). FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 3.
Average Individual Specific Growth Rate (%) for each dietary treatment. Different letters above the bars indicate significant differences (Tukey post-hoc test, p<0.05). FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 4.
Alpha diversity measures of fish fed different dietary treatments: (A) observed richness of ASVs, (B) Shannon diversity, (C) Chao1 diversity, (D) Pielou evenness. Different letters above the bars indicate significant differences (Tukey post-hoc test, p<0.05). FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 4.
Alpha diversity measures of fish fed different dietary treatments: (A) observed richness of ASVs, (B) Shannon diversity, (C) Chao1 diversity, (D) Pielou evenness. Different letters above the bars indicate significant differences (Tukey post-hoc test, p<0.05). FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 5.
NMDS plot of individual fish gut microbiome composition for each feed treatment. FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 5.
NMDS plot of individual fish gut microbiome composition for each feed treatment. FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 6.
Relative abundance of genera of taxa present in gut microbiome community of fish fed each dietary treatment. FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 6.
Relative abundance of genera of taxa present in gut microbiome community of fish fed each dietary treatment. FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 7.
Distance travelled (i.e., proxy for exploration) of fish fed different dietary treatments. FM is fish-meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 7.
Distance travelled (i.e., proxy for exploration) of fish fed different dietary treatments. FM is fish-meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 8.
Boldness score of fish fed different feed treatments where lower values indicate shyer individuals and high values indicate bolder individuals. FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 8.
Boldness score of fish fed different feed treatments where lower values indicate shyer individuals and high values indicate bolder individuals. FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 9.
Left side: Posterior modes and their corresponding uncertainty intervals for the variances and pairwise covariance of boldness and exploration for each diet treatment, on both the ID level and residual level. The black dashed line indicates zero. Right side: Posterior modes and their corresponding 95 % uncertainty intervals for the fixed effects SGR, TL, Tank and Trial.TL, SGR and Trial. Dashed line indicates zero. Because the non-categorial fixed effects were z-scored, the coefficient values are presented in units of standard deviation. FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Figure 9.
Left side: Posterior modes and their corresponding uncertainty intervals for the variances and pairwise covariance of boldness and exploration for each diet treatment, on both the ID level and residual level. The black dashed line indicates zero. Right side: Posterior modes and their corresponding 95 % uncertainty intervals for the fixed effects SGR, TL, Tank and Trial.TL, SGR and Trial. Dashed line indicates zero. Because the non-categorial fixed effects were z-scored, the coefficient values are presented in units of standard deviation. FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide.
Table 1.
Chemical composition of the protein sources used in feed treatments for this study.
Table 1.
Chemical composition of the protein sources used in feed treatments for this study.
Table 2.
Feed formulation and chemical composition for feed treatments for this study.
Table 2.
Feed formulation and chemical composition for feed treatments for this study.
Table 3.
Repeatability estimates (R) in the four dietary treatments. FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide. Boldness scores and exploration were obtained from two consecutive open field tests repeated 7 days apart. CI: confidence interval.
Table 3.
Repeatability estimates (R) in the four dietary treatments. FM is fish meal-based control diet, US is an untreated soybean meal diet, ETS is an enzyme treated soybean meal, and USP is an untreated soybean meal with fructose oligosaccharide. Boldness scores and exploration were obtained from two consecutive open field tests repeated 7 days apart. CI: confidence interval.
Diet treatment |
R |
CI |
FM |
|
|
Boldness |
0.84 |
[0.49-0.94] |
Exploration |
0.82 |
[0.52-0.94] |
ETS |
|
|
Boldness |
0.49 |
[0.16-0.77] |
Exploration |
0.84 |
[0.64-0.94] |
US |
|
|
Boldness |
0.00 |
[0.00-0.62] |
Exploration |
0.72 |
[0.30-0.89] |
USP |
|
|
Boldness |
0.30 |
[0.00-0.61] |
Exploration |
0.85 |
[0.60-0.94] |