Ethical statement
All fish handling procedures were conducted in accordance with relevant guidelines on the care and use of animals for scientific purposes developed by the Singapore Institutional Animal Care and Use Committee (approval number: TLL(F)-23-001).
Fish and system
In two separate trials, two cohorts consisting of 150 barramundi (48.9 ± 4.2 g, 95 days post-hatching [dph]) and 150 tilapia (48.7 ± 11.4 g, 219 dph) were fed one of two experimental diets in triplicate tanks for 23 and 29 days, respectively. The barramundi cohort originated from the Barramundi Group hatchery, while the tilapia cohort came from a 6th generation of the Temasek Life Sciences Laboratory (TLL) breeding program that was introduced to Singapore from South Africa (University of Stellenbosch) in 2009. Fish were grown in nine round conical-bottom 300-L tanks, with 25 fish per tank, and subjected to a natural 12:12-hour photoperiod. The tanks were supplied with brackish water (5 ppt salinity) by a recirculated aquaculture system (RAS). The RAS included UV treatment and sand and biological filtration steps and supplied water to each tank at a flow rate of 12 L per minute (resulting in 250% tank renewal per hour). Parameters such as dissolved oxygen, salinity, temperature, pH, nitrate, nitrite, and total ammonia in the system water were measured daily. After 14 days of acclimatization in the system, the fish were individually weighed (± 0.1 g) following light sedation with MS-222 (100–200 mg L− 1), and the same weighing procedure was repeated on the last day of the trial. On the penultimate day of each trial, two fish were sampled from each tank, euthanized in MS-222 (500 mg L− 1) 1 hour after feeding, and their whole digestive tracts were collected and stored at − 80°C for proteomic analysis.
Diet production
The BPM was a commercial product (China) made of a dried and autolyzed mixture of bacteria grown on brewery by-products. A control diet (CTL) was formulated with Wittaya software (Wittaya Aqua International Inc., 2022) to be as compatible as possible with the nutritional requirements of barramundi 10 and saline tilapia in terms of digestible proteins, energy, and lipids. A CTL premix was prepared at the R&D feed mill of the Marine Aquaculture Centre (Singapore). Ingredients were hammer milled (Mikro AP-1SH, Hosokawa Micron, USA) and sieved to a particle size of 500 µm using a vibratory sifter. Subsequently, the ingredients were thoroughly dry-mixed by utilizing a powder mixer (KSE-PM100, Kong Shiang Engineering, Singapore), and yttrium oxide (Y2O3) was added as an indigestible marker. BPM diets were prepared by incorporating 30% of BPM into the CTL premix (Table 1). In a 60-L mixer (HL-600, Hobart, Pinkenba, Australia), premixes were mixed with fish oil and screw-pressed through a 3.175-mm diameter die using a meat mincer (HL-400, Hobart, Australia). The diets were subsequently dried at 65°C in an oven (Memmert VO200) for 24 hours, bagged, and stored at 4°C until use in the trials. Tested ingredients and diets were subjected to proximate composition analysis (Table 2).
Table 1. Formulation (g kg-1) of experimental diets for both tilapia and barramundi1
|
CTL
|
BPM
|
Fish meal2
|
259.5
|
181.6
|
BPM3
|
|
300.0
|
Soybean meal
|
381.3
|
266.9
|
SPC4
|
187.2
|
131.0
|
Soy lecithin
|
30.0
|
21.0
|
Rice flour
|
89.8
|
62.9
|
Sardine oil5
|
43.8
|
30.7
|
Vitamin premix6
|
5.0
|
3.5
|
Mineral premix6
|
2.5
|
1.7
|
Yttrium oxide7
|
1.0
|
0.7
|
1 CTL: Control diet; BPM: Bacterial protein meal
2 FF-Classic, Skagen, Denmark, 70% crude protein
3 Bacterial protein meal made of a mix of autolyzed bacteria grown on brewery by-products (China)
4 Soy protein concentrate, X-Soy 200 (C.J. Selecta, Brazil), 60% crude protein
5 Sardine oil contained no antioxidants
6 ANA Fish Vitamin Premix-199 and ANA Fish Mineral Premix-199 (Zagro, Singapore)
7 Regent Chemicals Pte. Ltd., Singapore
Feeding and waste collection
Operators manually distributed the experimental diets (Table 1) at 11:00 AM five days a week at a feeding rate of 1.8% of body weight for barramundi and 2.5% of body weight for tilapia. Feed amounts were progressively adjusted based on the known thermal-unit growth coefficient (TGC) for each species and refined further based on weight gained during the acclimatization period. The uneaten pellets were removed and quantified after 30 min of feeding. During the last 6 days of the acclimatization period, fish were fed the control diet to minimize the influence of the commercial diet. As per the recommendation of Blyth et al. (2015), fecal collection started on day 4. Feces within the conical-shaped-bottom tanks accumulated in a customized collector at the base of each tank, which was collected daily over 13 days. The fecal samples were centrifuged at 3,900 × g for 10 minutes at 4°C, with the supernatant being discarded, and the pellets were frozen at − 20°C before being analyzed for their chemical content.
Chemical analysis
Chemical analyses of diets and feces were conducted by an accredited analytical service laboratory (Eurofins, Singapore) following the Association of Official Analytical Chemists methods 12. Dry Matter (DM) content was calculated by gravimetric analysis (AOAC 925.09) following oven drying at 100°C for 5 hours and desiccator cooling. Crude lipid (CL) content was estimated by acid hydrolysis (AOAC 922.06). Ash content was determined gravimetrically following the loss of mass after combustion of the sample at 550°C (AOAC 923.03). Total nitrogen (N) was estimated after pyrolysis and combustion (AOAC 968.06), and crude protein (CP) levels were based on N × 6.25. Crude fiber (fiber) was estimated after the digestion of the sample with sulfuric acid and sodium hydroxide (AOAC 962.09). Nitrogen-free extract (NFE) was calculated using the DM-CP-CL-fiber-ash formula, and gross energy (GE) was determined using the formula 23.6 × CP + 39 × CL + 17.6 × (fiber + NFE), as detailed by the National Research Council (1993). Fatty acids (FAs) in the samples were methylated to FA methyl-esters and quantitively measured by capillary gas chromatography (AOAC 996.06). The complete list of measured FAs used to compute omega-3 (n-3), omega-6 (n-6), omega-7 (n-7), omega-9 (n-9), saturated fatty acid (SFA), mono-unsaturated fatty acid (MFA), poly-unsaturated fatty acid (PUFA), and long-chain poly-unsaturated fatty acid (LC-PUFA) levels is provided in Table A1. Following digestion in hydrofluoric acid at 80°C, yttrium was quantified by a separate analytical service laboratory (Latech, Singapore) using inductively coupled plasma mass spectrometry. The amino acid composition of samples was determined by acid hydrolysis prior to separation using high-performance liquid chromatography, as explained by Schuster (1988).
Table 2. Composition (g kg-1 DM unless otherwise indicated) of ingredient and experimental diets
|
Ingredient
|
|
Tilapia trial
|
|
Barramundi trial
|
|
BPM
|
|
CTL
|
BPM
|
|
CTL
|
BPM
|
DM1
|
936.2
|
|
939.6
|
927
|
|
948.0
|
939.7
|
CL
|
32
|
|
111.7
|
93.9
|
|
115.0
|
87.3
|
Fiber
|
2.6
|
|
18.4
|
11.8
|
|
23.0
|
18.9
|
CP
|
461.4
|
|
520.4
|
504.9
|
|
526.4
|
508.7
|
Ash
|
132.5
|
|
89.8
|
102.7
|
|
90.2
|
103.2
|
NFE2
|
371.5
|
|
259.6
|
286.8
|
|
245.5
|
281.9
|
GE3
|
18.7
|
|
21.5
|
20.8
|
|
21.6
|
20.7
|
Yttrium oxide
|
|
|
7.1
|
4.5
|
|
11.0
|
5.9
|
|
|
|
|
|
|
|
|
Essential amino acids
|
|
|
|
|
|
|
|
Lysine
|
13.6
|
|
43
|
36.9
|
|
32.0
|
26.7
|
Arginine
|
24.8
|
|
36.9
|
33.4
|
|
36.4
|
30.9
|
Histidine
|
6.6
|
|
13.8
|
10.9
|
|
13.4
|
10.7
|
Isoleucine
|
18.1
|
|
23.1
|
20.9
|
|
24.8
|
22.5
|
Leucine
|
30.0
|
|
41.7
|
37.9
|
|
43.0
|
38.2
|
Valine
|
26.1
|
|
25.6
|
25.0
|
|
27.4
|
26.8
|
Methionine
|
10.4
|
|
11.2
|
11.0
|
|
12.0
|
11.4
|
Phenylalanine
|
17.2
|
|
25.4
|
22.9
|
|
26.2
|
23.1
|
Threonine
|
19.3
|
|
23.4
|
22.0
|
|
23.3
|
21.5
|
|
|
|
|
|
|
|
|
Non-essential amino acids
|
|
|
|
|
|
|
|
Alanine
|
32.9
|
|
29.4
|
30.7
|
|
29.2
|
30.3
|
Aspartic acid
|
37.0
|
|
59.2
|
52.2
|
|
58.3
|
50.3
|
Cystine
|
2.5
|
|
7.3
|
4.7
|
|
6.6
|
5.4
|
Glutamic acid
|
45.5
|
|
86.8
|
73.6
|
|
90.9
|
74.6
|
Glycine
|
19.7
|
|
31.1
|
30.0
|
|
26.1
|
24.6
|
Proline
|
16.3
|
|
33.6
|
32.6
|
|
25.4
|
22.0
|
Serine
|
13.9
|
|
26.3
|
22.3
|
|
26.2
|
21.3
|
Tyrosine
|
14.4
|
|
18.4
|
17.4
|
|
19.5
|
15.0
|
|
|
|
|
|
|
|
|
Fatty acids4
|
|
|
|
|
|
|
|
SFA
|
2.7
|
|
30.9
|
21.7
|
|
30.5
|
19.5
|
MFA
|
8.2
|
|
21.7
|
18.7
|
|
21.2
|
16.4
|
PUFA
|
0.0
|
|
38.8
|
25.8
|
|
38.2
|
24.5
|
LC-PUFA
|
0.0
|
|
20.3
|
12.9
|
|
19.6
|
12.8
|
Omega-3
|
0.0
|
|
21.9
|
14.1
|
|
21.3
|
13.8
|
Omega-6
|
0.0
|
|
16.8
|
11.7
|
|
16.9
|
10.7
|
Omega-7
|
2.0
|
|
5.4
|
4.5
|
|
5.4
|
4.1
|
Omega-9
|
0.0
|
|
13.8
|
9.9
|
|
13.2
|
8.8
|
TFA
|
10.9
|
|
104.6
|
76.2
|
|
102.9
|
68.9
|
1 Dry matter in g kg− 1.
2 NFE: nitrogen-free extracts
3 Gross energy in MJ kg− 1 DM
4 SFA: saturated fatty acids, MFA: mono-unsaturated fatty acids, PUFA: poly-unsaturated fatty acids, LC-PUFA: long-chain poly-unsaturated fatty acids, TFA: total fatty acids. Refer to Table A1 for the composition of fatty acid groups.
Enzyme assays
Three fish from each tank were sampled, and their entire digestive tracts were mixed with 10 ml of deionized water. Intestinal enzymes were extracted by homogenizing the mixture in a 15 ml Dounce homogenizer and centrifuging it at 21,000 × g for 30 minutes at 4°C. The supernatant was stored at − 80°C for further enzyme activity analyses. The protein concentration (in mg ml− 1) of each sample was determined via the Bradford assay using standards of bovine serum albumin (Quick Start Bradford Protein Assay Kit, Bio-Rad 5000201). Protease activities were measured using standard colorimetric substrates and the methods described by Natalia et al. (2004), adapted for 96-well format using a microplate reader. Assays were performed at room temperature by mixing the crude extracts with substrates in 50 mM Tris, pH 7.4, in a total volume of 350 µl. The substrates used with trypsin, chymotrypsin, and leucine aminopeptidase were, respectively, 0.1 mM Nα-benzoyl-DL-arginine 4-nitroanilide hydrochloride (BAPNA) + 20 mM CaCl2, 0.1 mM N-succinyl-L-alanyl-L-alanyl-L-prolyl-L-phenylalanine 4-nitroanilide (SAPNA) + 20 mM CaCl2, and 1.0 mM L-leucine-4-nitroanilide (LeuPNA) + 20 mM MgCl2. Proteolytic reactions were followed over time by recording the absorbance at 410 nm. Lipase activity was measured using a commercial kit (Sigma-Aldrich, MAK046-1KT), where the production of glycerol from triglyceride by lipase was monitored using absorbance at 570 nm. One enzyme unit (U) was defined as one nmol of product produced per minute. For trypsin, chymotrypsin, and leucine aminopeptidase, the conversion of absorbance to nmol of p-nitroaniline was calculated using the molar extinction coefficient of 8,800 M− 1cm− 1 14. For lipase, the molar amount of glycerol was derived from a standard curve obtained by measuring the absorbance at 570 nm of known concentrations of glycerol. Specific activity was expressed as U mg− 1 total protein present in the crude extract. The enzyme kinetic parameters, namely the maximum velocity (Vmax, in U mg− 1) and the Michaelis–Menten constant (Km), were determined for trypsin, chymotrypsin, and leucine aminopeptidase. This was accomplished using concentration gradients of BAPNA, SAPNA, and LeuPNA as substrates.
Data collection
The mean values for the initial and final fish body weights (FBW, in g) were employed to compute both the body weight gain (BWG, in g) and the TGC for each tank. TGC served as a standardized measure of growth 15 and was assumed to be unaffected by variations in body weight, time intervals, and water temperature 16,17. In this analysis, the TGC was calculated with a base value of 20 to enhance its stability, which is especially important in high-temperature conditions 18:
\(\:{TGC}_{a}=1000\text{*}\frac{{FBW}_{a}^{(1/3)}-\:{IBW}_{a}^{(1/3)}}{{\sum\:}_{\text{i}=0}^{51}({\text{T}}_{\text{i}}-20)}\) (Eq. 1)
where FBWa and IBWa are the final and initial average FBW (g) in tank a, and Ti (°C) is the water temperature on day i. The feed conversion ratio (FCR) in tank a was calculated as follows:
\(\:{FCR}_{a}=\left(\sum\:_{i=1}^{51}{TI}_{ai}/{qty\:of\:fish}_{ai})/{(FBW}_{a}-{IBW}_{a}\right)\) (Eq. 2)
where TIai is the tank feed intake (g day− 1) in tank a on day i. The apparent digestibility coefficients (ADCs) (%) of DM, CP, CL, GE, nutrients, or nutrient groups in the experimental diets were calculated as follows:
\(\:{ADC}_{parameter}=\left[1-\left(\frac{{Y}_{Diet}\times\:\:{Parameter}_{feces}}{{Y}_{feces}\times\:\:{Parameter}_{diet}}\right)\right]\times\:100\) (Eq. 3)
where \(\:{Y}_{diet}\) and \(\:{Y}_{feces}\) are the yttrium content of the diet and feces, respectively, and \(\:{Parameter}_{diet}\) and \(\:{Parameter}_{feces}\) are the content of the nutritional parameter of interest (DM, CP, CL, GE, nutrient, or nutrient group) in the diet and feces, respectively. The ADCs of a given nutrient from BPM ingredient (ADCingr) were calculated as follows:
\(\:{ADC}_{ingr}=\left({ADC}_{test}{\times\:Nutr}_{test}-\left({ADC}_{ref}{\times\:Nutr}_{ref}\times\:Rb\right)\right)/\left(\left(1-Rb\right){\times\:Nutr}_{ingr})\right)\) (Eq. 4)
where ADCtest is the ADC of the diet containing the ingredient, ADCref is the ADC of the reference diet. Nutrtest, Nutrref, Nutringr are the level of the nutrient of interest in the test diet, reference diet and the ingredient respectively 19. The ratio 1 − Rb at which the test ingredient was included in the test diet was determined based on the effective dilution of yttrium between the basal and test diets, according 20.
Statistical analysis
To analyze fish performance and ADCs, one-way analysis of variance (ANOVA) was employed, and pairwise comparisons between species or diets were conducted using Tukey's honestly significant difference (HSD) multiple range test. An arcsine-square-root transformation was applied to the survival rate and ADCs. All kinetic parameters were obtained by subjecting the data to the non-linear regression methods described by the Michaelis–Menten equation: V = Vmax([S]/ (Km + [S])) (Eq. 5), where V is the velocity of the reaction, [S] is the substrate concentration, Vmax is the maximum velocity, and Km is the Michaelis–Menten constant. This methodology was preferred over the Lineweaver-Burk method to improve the accuracy of kinetic parameter estimates 21. For All statistical analyses were performed using R software 22, and effects were considered statistically significant at a significance level of P < 0.05.