The experiment was reviewed and approved by the Committee of Ethics and Animal Welfare of the Universitat Politècnica de València (UPV), following the Spanish Royal Decree 53/2013 on the protection of animals used for scientific purposes [31].
A total number of 240 juveniles of gilthead seabream (averagare weigth 7,5 g and 60 days) were obtained from the fish farm BERSOLAZ (Bersolaz Spain, S.L.U, Culmarex Group) located in Port de Sagunt (Valencia, Spain) and transported to the facilities at the Universitat Politècnica de València, where the growth trial was conducted after 15 days of adaptation to experimental conditions. Features of the system and water parameters set were described in previous growth trials carried out in these facilities [32,33]. Lighting conditions were determined by the natural photoperiod. Temperature, pH, oxygen, ammonia, nitrite and nitrate concentrations were monitored along the experiment. The fish were daily fed by hand to apparent satiation two times per day (9:00h and 17:00h). The pellets were slowly distributed, allowing fish to eat, in a weekly regime of six day of feeding and one day of fasting.
Diets were prepared by cooking extrusion process using a semi-industrial twin-screw extruder (CLEXTRAL BC-45, St. Etienne, France). A fish meal based control diet (FM), in which most of the protein was provided by fish meal (59%), and a plant protein based diet (PP), in which all the fish meal was replaced by plant sources and synthetic amino acid were added to meet the minimum amino acid requirement for gilthead seabream juveniles [34]. Ingredients and proximate composition are shown in Table 1. Prior to diet formulation, dry matter, crude protein, crude lipid, ashes and crude fibre (CF) of different sources and ingredients used were analysed according to AOAC procedures [35]. All analyses were performed in triplicate. Amino acids of raw diets were also analysed by reverse phase – high performance liquid chromatography [36]. Macronutrients and essential amino acid content were determined in the experimental diets, and they are shown in Table 1.
Table 1
Ingredients, proximal composition and essential amino acids profile of experimental diets
|
FM
|
PP
|
Ingredients (g kg-1)
|
|
|
Fish meal
|
589
|
|
Wheat meal
|
260
|
|
Wheat gluten
|
|
295
|
Broad bean meal
|
|
41
|
Soybean meal
|
|
182
|
Pea meal
|
|
41
|
Sunflower meal
|
|
158
|
Krill meal
|
|
|
Squid meal
|
|
|
Fish oil
|
38.1
|
90
|
Soybean oil
|
92.9
|
90
|
Soy Lecithin
|
10
|
10
|
Vitamin-mineral mix*
|
10
|
10
|
Calcium phosphate
|
|
38
|
Arginine
|
|
5
|
Lysine
|
|
10
|
Methionine
|
|
7
|
Taurine
|
|
20
|
Threonine
|
|
3
|
Proximate composition (% dry weight)
|
|
|
Dry matter
|
88.1
|
93.9
|
Ashes
|
10.1
|
7.4
|
Crude lipid
|
18.5
|
19.8
|
Crude fiber
|
0.8
|
4.3
|
Crude protein
|
44.2
|
45.0
|
Essential aminoacids (g 100 g-1)
|
|
|
Arginine
|
3.39
|
3.30
|
Histidine
|
1.00
|
0.82
|
Isoleucine
|
1.47
|
1.17
|
Leucine
|
3.24
|
2.98
|
Lysine
|
3.68
|
2.26
|
Methionine
|
1.16
|
1.06
|
Phenylalanine
|
1.80
|
1.87
|
Threonine
|
1.98
|
1.44
|
Valine
|
2.01
|
1.47
|
Bacterial strains and growth conditions
Cultures of Pseudomonas anguilliseptica CECT 901, Vibrio algynoliticus CECT 521 and Photobacterium damselae subsp. piscicida CECT 7198 were obtained from the Colección Española de Cultivos Tipo (CECT, Valencia, Spain). The culture medium used for P. anguilliseptica was Tryptic Soy Broth (containing per litre of water, 15.0 g tryptone, 5.0 g soy peptone and 5.0 g NaCl, pH 7.3), and V. algynoliticus and P. damselae subsp. piscicida were grown on Marine Broth (containing per litre of water, 5.0 g Bacto peptone, 1.0 g yeast extract, 0.10 g Fe(III) citrate, 19.45 g NaCl, 0.16 g Na2CO3, 3.24 g NaSO4, 1.80 g CaCl2, 8.80 g MgCl2, 0.55 g KCl, 0.08 g KBr, 34.00 mg SrCl2, 22.00 mg H3BO3, 4.00 mg Na-silicate, 2.40 mg NaF, 1.60 mg (NH4)NO3 and 8.00 mg Na2HPO4, pH 7.6). Bacteria were grown under agitation at 26º C for 2 days (P. anguilliseptica and P. damselae subsp. piscicida) and at 30º C for 1 day (V. algynoliticus). Then, 1.5 g/L of bacteriological agar were added to these media to prepare solid medium in petri dishes. For the bacterial challenge, optical density (600 nm) of the bacterial cultures was determined and bacterial cell number was estimated using the standard curves established for each strain. Then, bacterial cultures were centrifuged at 4.000 g for 20 min, washed once with PBS, and re-suspended in CO2-independent cell culture medium (Gibco, ThermoFisher) to a final concentration of 3·107 ufc/mL in the case of V. algynoliticus, and 1·107 ufc/mL of P. damselae subsp. piscicida and P. anguilliseptica.
Experimental design
The aim of this work was to evaluate the impact of dietary fish protein substitution by plant protein on the intestinal health status and its immune response capacity. For this purpose, fish were fed with fish meal (FM) or plant protein (PP) based diets, and animals were sacrificed and processed at two critical growth phases of sea bream [16,17]: Phase I (90 days; from 12g to~ 68g) and Phase II (305 days; up to 250g). Each diet was assayed in tanks per triplicate. Fig 1. illustrates the experimental design. At each sampling time intestine fragments were collected, part of them were used to determine basal expression of selected genetic markers and other fragments were used in the culture explant procedure (ex vivo assay) to test the immune competence of the intestinal mucosa when challenged with different bacterial pathogens.
Phase I: Up to 68 g
Fish were fed with FM or PP diets, in tanks per triplicated (40 fish per tank), up to 90 days, being scattering 2 fish per group to submit bacterial challenge (68±37.8g). For basal gene expression two fragments from foregut (FG) and hindgut (HG) from each fish were placed in an eppendorf tube containing 500 µl of RNA Later® (Qiagen, Valencia, CA) for subsequent total RNA (tRNA) extraction. Additionally, four fragments from of FG and HG fragments were used for the ex vivo assay and exposed to pathogens challenge (see below). Gene expression was determined in all samples to evaluate the inflammatory and immune status of the intestinal mucosa due to changes in the diet and the bacterial challenge.
Phase II: Up to 250 g
Fish, ~30 fish per tank, were fed with the same diets, FM and PP, up to 305 days when the mean weight was 252±70.1g (Fig 1). In this second phase, in addition to the impact of long-term feeding with 100% of PP diets, also a short-term exposure (15 days) of total fish meal substitution was evaluated. For this purpose, fish bred with FM diet (n=~15 fish per tank) were changed to a PP diet two weeks before the termination of the experiment (from day 290 to day 305) (PP* group). Three fish from the FM group and two from the PP and PP* were sacrificed to obtain FG and HG explants for ex vivo assays. As in the previous assay, for basal gene expression two fragments from each fish were placed in RNA Later® for subsequent total RNA (tRNA) extraction, and four pieces for ex vivo study.
Ex vivo assays and bacterial challenge
Before tissue preparation, fish were sacrificed by immersion in benzocaine (60 ppm) during 15 min. Then, they were dissected and the intestine was obtained and separated in two sections (foregut and hindgut). Each section was cut with a scalpel in small pieces (4 mm x mm), which were immediately placed in culture filter plates (15 mm diameter wells with 500 µm bottom-mesh, Netwell culture systems, Costar, Cambridge, MA) with the epithelial surface facing up. Filters were placed into wells containing 1 mL of the different bacterial solutions (one of them was preserved without bacteria as control; Ex vivo Unchallenged group) in CO2-independent cell culture medium (Gibco, ThermoFisher). 100 µL of the corresponding bacterial solutions were finally added to epithelial surface to ensure that samples were completely submerged. At the end of the incubation time, samples were carefully collected from the culture filter plates and stored in 100 mM Tris-HCl at 4ºC or RNA later at -80º C for lactate dehydrogenase (LDH) activity evaluation or RNA isolation, respectively. Changes in pH of the explant culture medium due to different bacterial treatments were monitored. Explants of foregut (FG) and hindgut (HG) from two fish per group were incubated during 4 and 6 h at 22º C in independent CO2 atmosphere, depending on the experiment. The bacterial species used in the pathogen challenge were: Photobacterium damselae subsp, Pseudomonas anguilliseptica and Vibrio algynoliticus. Pseudomonas anguilliseptica was discarded in phase II, because bacterial concentration could not be determined due to aggregate formation. After explant assay, the samples were placed into RNA Later (Qiagen, Valencia, CA) for subsequent tRNA extraction. All conditions (fish/section/stimuli) were assayed in duplicate and gene expression was determined in all samples to evaluate the intestinal inflammatory and immune status based on experimental diet and bacterial challenge.
LDH activity assay
In order to determine the tissue integrity, the LDH activity was determined [37] in the tissue (U/mg protein) and explant culture medium (U/L) at different times of the incubation (0, 4, 6 and 24 h). LDH activity was analysed measuring the nicotinamide adenine dinucleotide (NADH) absorbance at 340 nm using the commercial kit (BioSystems S. A., Barcelona, Spain). Tissue was weighed, homogenised in Tris-Hcl 100 mM while maintaining the tubes on ice, centrifuged at 12.000 rpm and 4º C for 15 min and supernatant was collected for LDH assessment. Total protein in tissue extracts was determined using Bradford [38].
Gene expression assay of intestinal inflammatory and immune markers
Based on the gene expression analysis used in previous studies in this species (Sparus aurata) to evaluate the intestinal inflammatory and immune status [32] tRNA was extracted from intestinal tissue samples using the phenol/chloroform method with Trizol Reagent (Invitrogen, Spain) and treated with DNAse I (Roche) to remove DNases. Total RNA concentration, quality and integrity were assessed using a NanoDrop 2000C Spectrophotometer (Fisher Scientific SL, Spain). The integrity of 28S/18S was also determined by gel electrophoresis. 1 µg of total RNA was used for cDNA synthesis reaction using the qScript cDNA synthesis kit (Quanta BioScience), according to the manufacturer’s instructions. An Applied Biosystems 2720 Thermal Cycler was used with the following cycling conditions: 22 ºC for 5 min, 42 ºC for 30 min, and 85 ºC for 5 min. cDNA samples were stored at -20º C until gene expression was analysed.
Four housekeeping candidate genes (Table 2) were tested to be used as reference genes, and for assessing RNA integrity along the assay. The Cq of the four genes was determined in six pooled samples from Experiment 1 (two dietary groups: FM and PP; three times: 0, 4 and 6 h). Relative gene expression of six genes was determined in the foregut and hindgut samples. The genetic markers monitored in this assay were three pro-inflammatory markers, IL1-β, IL-6 and COX-2, the main immunoglobulin, IgM, the main intestinal mucin, I-Muc, and the occludin gene, Ocl, with primers listed in Table 2.
Table 2
Primer sequences of candidate genes (reference and target genes) in the RT-qPCR assay
Gene
|
Abbreviation
|
GeneBank ID
|
Primer Forward
|
Primer Reverse
|
Lenght
|
Reference
|
REFERENCE GENES
|
|
Elongation Factor 1α
|
EF-1α
|
AF184170
|
CTGTCAAGGAAATCCGTCGT
|
TGACCTGAGCGTTGAAGTTG
|
87
|
[35, 36]
|
Glyceraldehide 3-phosphate dehydrogenase
|
GAPDH
|
DQ641630
|
CCAACGTGTCAGTGGTTGAC
|
AGCCTTGACGACCTTCTTGA
|
80
|
[37]
|
Ribosomal Protein S18
|
Rps18
|
AM490061
|
AGGGTGTTGGCAGACGTTAC
|
CGCTCAACCTCCTCATCAGT
|
97
|
[37]
|
β-Actin
|
β-Act
|
X89920
|
TCTGTCTGGATCGGAGGCTC
|
AAGCATTTGCGGTGGACG
|
113
|
[38]
|
TARGET GENES
|
|
Interleukin 1β
|
IL-1β
|
AJ277166
|
GCGACCTACCTGCCACCTACACC
|
TCGTCCACCGCCTCCAGATGC
|
131
|
[37]
|
Interleukin 6
|
IL-6
|
AM749958
|
AGGCAGGAGTTTGAAGCTGA
|
ATGCTGAAGTTGGTGGAAGG
|
101
|
[35]
|
Cyclooxygenase 2
|
COX-2
|
AM296029
|
GAGTACTGGAAGCCGAGCAC
|
GATATCACTGCCGCCTGAGT
|
192
|
[35, 36]
|
Intestinal Mucin
|
I-Muc
|
JQ277712
|
GTGTGACCTCTTCCGTTA
|
GCAATGACAGCAATGACA
|
102
|
[38]
|
Immunoglobulin M
|
IgM
|
JQ811851
|
TCAGCGTCCTTCAGTGTTTATGATGCC
|
CAGCGTCGTCGTCAACAAGCCAAGC
|
131
|
[39]
|
Occludin
|
Ocl
|
JK692876
|
GTGCGCTCAGTACCAGCAG
|
TGAGGCTCCACCACACAGTA
|
81
|
[35, 36]
|
All qPCR assays and expression analyses were performed using the Applied Biosystems 7500 Real-Time PCR with SYBR® Green PCR Master Mix (ThermoFisher Scientific, Waltham, Massachusetts, USA). After an initial Taq activation of polymerase at 95 ºC for 10 min, 42 cycles of PCR were performed with the following cycling conditions: 95 ºC for 10 s and 60 ºC for 30 s in all genes. In order to evaluate assay specificity, a melting curve analysis was directly performed after PCR cycles by slowly increasing the temperature (1º C/min) from 60 to 95 ºC, with continuous registration of changes in fluorescent emission intensity. The total volume for every PCR reaction was 20 μl, performed from diluted (1:20) cDNA template (5 μl), forward and reverse primers (10 μM, 1 μL), SYBR® Green PCR Master Mix (10 μl), ROX (2 μL, 10 nM) and nuclease-free water up to 20 μl. The analysis of the results was carried out using the 2-∆∆Ct method. The target gene expression quantification was expressed relative to the expression of the selected reference gene. A cDNA pool from all the samples was included in each run and acted as a calibrator, and a non-template control for each primer pair, in which cDNA was replaced by water, was run on all plates. Reference and target genes in all samples were run in duplicate PCR reactions.
Statistics
Statistical data analysis were performed with Statgraphics© Centurion XVI software (Statistical Graphics Corp., Rockville, MO, USA).
LDH enzymatic activity in tissues and the supernatant was statistically analysed by one-way analysis of variance (ANOVA) using Newman-Keuls test to determine possible differences across the assay (0, 4, 6 and 24 hours) in FG and HG.
The expression stability of reference genes was assessed using the BestKeeper program, basing on the arithmetic means of the Cq values [39]. Lower deviation in the expression is related to better stability.
The evaluation of intestinal inflammatory and immune status was performed through the gene expression of the target genes both in vivo and ex vivo conditions. The relative gene expression was statistically analysed by ANOVA. Gene expression of cultured pieces was normalised with the expression of ex vivo unchallenged samples at 4 and 6 hours. Multifactorial analysis was used to determine the significance (p<0.05) of different factors considered (dietary treatment: FM/PP, intestinal section: FG/HG and bacterial stimuli: P. damselae subsp. piscicida/P. anguilliseptica/V. algynoliticus) at different times and to determine differences in normalised gene expression between dietary groups, sections and bacterial stimuli, using Newman-Keuls test. Data was expressed with the mean and the standard error of the normalised expression values, and differences were considered statistically significant when p<0.05.
Additionally, with the aim to evaluate if the bacterial challenge is individually inducing the target genes or a combination of a set of genes, a correlation analysis was carried out and the Pearson product-moment coefficient was obtained for each pair of genes.
Finally, in order to confirm the assay reproducibility, the gene expression of the biological replicate samples was randomly assigned to different variables (x and y). Data consistency was evaluated for each gene by simple regression analysis using the model y = ax. 95% confidence intervals for a (a±1.96σ) were obtained for each gene to validate the hypothesis a=1 (y=x).