Biological material
The experiment was carried out at the Laboratorio de Fisiología en Recursos Acuáticos (LAFIRA) in DACBiol-UJAT. Spawning was induced to obtain A. tropicus larvae, using one female (3 kg) and six males (1.5 kg). Female was injected with GnRH (Ourofino Saúde animal, Ltda) (35 µg kgˉ¹) in the base of the pelvic fin. Fish were placed in a 2000 L circular pond with raffia thread to simulate the natural spawning site.
After hatching (72 h post fertilization), larvae were placed in a 70 l plastic tanks for maintenance. Once the nutrients have been absorbed from the yolk sac and after opening of the mouth (3 days post hatching), the larvae were fed 5 times per day (8:00, 11:00, 13:00, 15:00 and 18:00 h) with artemia nauplii following the feeding scheme proposed by De la Cruz-Alvarado et al. (2021) until reaching juvenile size. For this experiment, 180 juveniles (3.65 ± 0.12 g mean weight and 10.38 ± 0.10 cm total length) were used.
Experimental design
Experimental design consisted of a completely randomized single factor experiment. Four treatments were evaluated by triplicate with different concentrations of SP (0.5, 1.0, 1.5 and 2.0%) and a control group (0%). For each experimental unit, 12 organisms were placed in 15 circular tanks of 70 L connected to a recirculation system with a 1500 L reservoir for solids sedimentation and a biological filter driven by a 1 HP water pump (Jacuzzi, JWPA5D-230A, Delavan WI, USA). Water quality parameters were monitored daily until the end of the experiment. Temperature (29.5 ± 0.6°C) and dissolved oxygen (5.4 ± 0.3 mg/l) were measured with an oximeter (YSI 85, Ohio, USA), and pH (7.2 ± 0.1) with a potentiometer (HANNA HI 991001, Romania). The organisms were fed 4 times a day (8:00, 11:00, 14:00 and 18:00 h) at 6% of biomass weight during 61 days of feeding trial.
Experimental diets
A modified diet proposed by Frías-Quintana et al. (2016) was used as base diet. Different concentrations of SP were added according to treatments (Table 1). For the preparation of experimental diets, the methodology described by Álvarez-González et al. (2001) was followed. Macronutrients, micronutrients, and sodium propionate were weighed separately. Macronutrients were mixed for 15 min using an industrial mixer (Bathamex, 178716, Mexico), then the micronutrients and sodium propionate were added and mixed for another 15 min. Finally, liquid ingredients (soy lecithin) and water (400 ml kgˉ¹) were weighted, incorporated, and mixed for 15 min. The resulting mixture was processed in a grinder (Torrey, M-22RI, Monterrey, Mexico) to obtain 5 mm pellets, which were dried in an oven (Coriat, HC-35-D, CDMX, México) at 45°C. for 12 h and stored at -20°C until further analysis.
Table 1
Composition of experimental diets with diferent concentrations of SP
Sodium propionate (%)
|
Ingredients (g kg− 1)
|
0
|
0.5
|
1.0
|
1.5
|
2.0
|
Pork mealᵃ
Poultry mealᵃ
Soybean mealᵇ
Corn starchᶜ
Fish mealᵃ
Wheat mealᵇ
Sodium propionateᵈ
Soybean lecithinᵉ
Gelatinᶠ
Vitamins and minerals premixᶢ
Vitamin Cʰ
|
250
218.4
161.1
150
100
55.5
0
30
20
10
5
|
250
218.4
161.1
150
100
50.5
5
30
20
10
5
|
250
218.4
161.1
150
100
45.5
10
30
20
10
5
|
250
218.4
161.1
150
100
40.5
15
30
20
10
5
|
250
218.4
161.1
150
100
35.5
20
30
20
10
5
|
Proximal composition (g/100g dry matter)
|
Protein
Lipids
Ash
NFE
|
44.32
15.32
12.77
28.49
|
44.81
15.74
12.74
28.43
|
43.99
15.43
13.23
28.43
|
44.83
15.53
13.71
28.48
|
44.74
15.3
13.53
28.46
|
ᵃProteínas Marinas y Agropecuarias, S. A. de C. V., Guadalajara, Jalisco
ᵇGALMEX Comercializadora de Insumos Agrícolas, Villahermosa, Tabasco, México
ᶜMSA Industrializadora de Maíz, Guadalajara, Jalisco, México
ᵈMi granero Ingredientes Selectos. San Pedro Cholula, Puebla, México
ᵉPronat Ultra, Mérida, Yucatán, México
ᶠD’gari Productos alimenticios y dietéticos Relámpago, Tlalpan, CDMX, México
ᶢPedregal (para trucha Silver Cup), Toluca, Edo. de Méx. México
ʰROVIMIX® C-EC (Roche) active agent of 35%
Growth and feed indexes
Every 15 days, weight and total length of the organisms were measured until the end of the experiment (day 61). Weight was recorded with an electronic balance (Ohaus HH120, precision 120 ± 0.01 g, Shenzhen, China) and length was measured using a digital vernier. Survival was calculated by counting the juveniles at the end of the feeding trial. With the obtained data, the following indexes were calculated:
- $$Survival \left(S\%\right)=\frac{final fish number}{initial fish number }*100$$
- $$Absolute weight gain \left(AWG\right)=final weight \left(g\right)-initial weight\left(g\right)$$
- $$Specific growth rate \left(SGR\right)=\frac{(\text{ln} final weight-\text{ln} initial weight)}{days of feeding}*100$$
- $$Feed conversion rate \left(FCR\right)=\frac{feed intake \left(g\right)}{weight gain \left(g\right)}$$
- $$Condition factor \left(K\right)=\frac{final weight \left(g\right)}{total final lenght \left(cm\right)}*100$$
Biological samples
All procedures were performed according to Mexican Norm (NOM-062-ZOO-1999, 2001) of Animal Welfare. At the end of the feeding trial, 15 organisms per treatment were dissected to obtain the visceral package and register biometric data of the organs. For the digestive enzyme analysis, intestine and stomach from three fish per tank were dissected and stored at -80°C until further analysis. For gene expression analysis, intestine and liver from two fish per tank were dissected and stored in RNAlater at -80°C until processing.
Somatic indexes
With the data obtained from the visceral package, somatic indexes were calculated:
- $$Hepatosomatic index \left(HSI\right)=\frac{liver weight \left(g\right)*100}{body weight \left(g\right)}$$
- $$Viscerosomatic index \left(VSI\right)=\frac{viscera weight \left(g\right)*100}{body weight \left(g\right)}$$
- $$Relative intestine length \left(RIL\right)=\frac{intestine length \left(cm\right)*100}{body length \left(cm\right)}$$
Digestive enzyme activity
The crude enzymatic extract was obtained from the stomach and intestine, which were macerated and homogenized separately with distilled water in a 1:10 ratio (weight: volume) under cold conditions (4°C). The mixture was centrifuged at 12,000 rpm for 15 min at 4°C and the supernatant was aliquoted and frozen at -80°C until analysis. Soluble protein concentration was determined following Bradford (1976) protocol using bovine serum albumin (BSA) as the standard protein. Acid proteases (stomach) were determined according to Anson (1938) methodology, using hemoglobin (0.25%) in L–1 glycine-HCl buffer (100 mM, pH 2) as substrate. The activity of alkaline proteases in intestine was determined by the method described by Walter (1984) using casein (0.25%) as substrate in 100 mM Tris-HCL buffer + 10 mM Ca Cl₂ (pH 9). For both techniques, absorbance was measured at 280 nm using a molar extinction coefficient (MOC) of 0.005 ml/µM cm. The methodology described by Erlanger et al. (1961) was followed to determine trypsin activity using 2 mM α-Benzoyl-DL-Arginine- ƿ -nitroanilide (BAPNA) as substrate diluted in 50 mM Tris-HCL + 10 mM CaCl₂ buffer (pH 8.2), absorbance was measured at 410 nm. Chymotrypsin activity was determined according to DelMar et al. (1979) using 1.25 mM SAAPNA substrate in 100 mM Tris-HCL + 100 mM CaCl₂ buffer (pH 7.8), absorbance was measured at 410 nm. Amylase activity was determined according to Robyt & Whelan (1968) using 2% starch as substrate in 100 mM citrate-phosphate buffer + 50 mM NaCl, pH 7.5; the reaction product was measured at 600 nm. Lipase activity was determined according Versaw et al. (1989) technique using 100 mM β-naphthyl acetate as substrate in 50 mM Tris-HCL buffer (pH 7.5) with 100 mM sodium cholate, the reaction was measured at 540 nm. Alkaline and acid phosphatases were determined following methodology described by Bergmeyer (1974) using 4-nitrophenyl phosphate (2.4%), dissolved in 100 mM L–1 glycine-NaOH (pH 10.1) buffer as substrate for alkaline phosphatases and 100 mM citric acid (pH 5.5) for acid phosphatases; the absorbances was measured at 410 nm. The enzyme activity was determined using the following equations: units by ml (U/ml) = [Δabs × final reaction volume (ml)]/[ε × time (min) × extract volume (ml)] − 1; specific activity (U/mg protein) = U ml/ mg of soluble protein.
RNA extraction and qPCR
RNA was extracted from intestine and liver samples using Trizol Reagent according to the manufacturer's protocol (Invitrogen, Waltham, MA). RNA concentration and purity were measured with a spectrophotometer (NanoDrop2000, Thermo Fisher Scientific, Madrid, Spain) using an absorbance radius of 260 and 280 nm. The cDNA synthesis was performed with the Superscript II kit (Invitrogen), using 1 µg of RNA for a final volume of 20 µl. Reactions were performed in a thermocycler (Mastercycler nexus GSX1, Eppendorf AG, Hamburg, Germany) using the following protocol: 5 min at 65°C, 10 min at 25°C (alignment), 50 min at 42°C (cDNA extension), 15 min at 70°C (RT activation) and finally 20 min at 37°C.
Genes related to intestinal health occludin (ocln), nucleotide-binding oligomerization domain-2 (nod2) and mucin 2 (muc2) were analyzed (Table 2).
The RT-qPCR reaction was performed using a CFX96 Real-Time System (BioRad, Hercules, CA) using 5µl of SyberGreen, 3µl of primer mix and 2µl of cDNA for a final volume of 10µl. The protocol followed was: 95°C for 10 min, 40 cycles of 15 s at 95°C, 60°C for 30 s, 70°C for 5 s and finally a melting curve with increments of 0.5°C from 75°C to 95°C. The rpl8 gene was used as the reference gene. Relative gene expression was calculated as fold-change compared with control and using 2−ΔΔCt formula (Livak & Schmittgen, 2001).
Table 2
Designed primers for qPCR of genes from intestinal barrier of A. tropicus
Gen
|
Primer sequence (5’-3’)
|
Size (bp)
|
Amplification efficiency (%)
|
Reference
|
Ribosomal protein (rpl8)
|
F: TGTGCTGCCTGGAAGAGAAG
R: TTTCGGGGTTGTGGGAGATG
|
90
|
99.82
|
Aranda-Morales et al. (2021)
|
Occludin (ocln)
|
F: TGACGAATACCACAGACTGAAG
R: CGATCATAGTCGCTGACCATC
|
123
|
|
(Sepúlveda-Quiroz et al., 2021)
|
Mucin (muc2)
|
F: GGCCTCCTCAAGAGCACGGTG
R: TCTGCACGCTGGAGCACTCAATG
|
100
|
90.94
|
(Nieves-Rodríguez et al., 2018)
|
Nucleotide-binding oligomerization domain-2 (nod2)
|
F: GTAGTGAACAAGGAGGCGGAC
R: TGAGCTCATCCAGGCCATCG
|
295
|
|
(Sepúlveda-Quiroz et al., 2021)
|
Statistical Analysis
Normality (Kolmogorov-Smirnov) and homoscedasticity (Barttlet) tests were performed for all variables. One-way ANOVA was performed to analyze the variables (growth and digestive enzyme activity). A posteriori test (Tukey) was performed when finding significant differences. Gene expression variables were analyzed using the Kruskal-Wallis and Dunn´s nonparametric test. Statistical analyses were performed in the GraphPad Prism software (v.8.0.2) with a significance value of 0.05.