Animals and housing
All pigs were handled and treated in accordance with the guidelines of the Canadian Council on Animal Care (57). The protocol was approved by the Animal Use Ethics Committee of the Faculté de médecine vétérinaire of the Université de Montréal (certificate number 17-Rech-1858).
A total of 352 individually ear-tagged pigs (Landrace x Large White hybrid sows sired with Duroc x synthetic hybrid boar), with undocked tail and of 8 weeks of age, were randomly distributed into 32 growing-finishing pens (2.06 × 3.35 m) located in two separate rooms (16 pens/room) at a commercial farm in Ange-Gardien, Quebec, Canada (Additional file 1). Each pen housed 11 pigs (6 gilts and 5 barrows). To facilitate video identification, in each pen, a 30 cm number ranging from 0 to 10 was painted on each pig back. Each number was associated to each pig ear-tag identifier.
In accordance with the feeding practices in place at the commercial farm, pigs were fed a corn- and soy-based granulated diet. The feeding protocol consisted of a pre-fattening feed (supplemented with 0.5 kg/ton of salinomycin) for 2.5 weeks, a fattening feed (supplemented with 0.5 kg/ton of salinomycin) for 3.5 weeks, a first growth feed (supplemented with 0.21 kg/ton of salinomycin) for 3.5 weeks, a second growth feed (supplemented 0.21 kg/ton of salinomycin) for 10 days, and a finishing feed (supplemented with 0.15 kg/ton of narasin) (F. Ménard, Inc., Ange-Gardien, QC, Canada) until slaughter. Animals had free access to food and water at all times. Ambient conditions, ventilation and temperature, were established according to standard housing recommendations.
Behavior Assessment
Tail-biting behaviors were recorded on video cameras (1080p resolution, day/night) with a recording capacity of 7 days. A total of 16 cameras were installed above the 32 growing-finishing pens under study in a position that allowed the recording of two pens simultaneously. Recordings were done 24 h/24 h and transferred to a computer every 48 h for analysis of pig behavior. A 2-minute interval sampling on a 24 h video in four pens from which animals presented TB behavior was done as a preliminary test to determine the periods when TB occurred as recommended (58).
Subsequently, five pens where TB behavior was present were selected for developing an ethogram (Table 1) based on a pre-observation period of 48 h video recordings. The main behaviors presented in the ethogram included tail accessibility, nibbling of the tail by the biter, duration of the nibbling behavior, frequency of tail nibbling, time when the nibbling occurred and the receiver’s response.
In order to further define the groups of animals to be studied, the “target behavior sampling” method (58), centered on observations of the behavior of interest, here TB, was used on recordings of the 10 pens where TB behavior was present, based on the ethogram. Sampling observations were distributed among two periods: from 10:00 am to noon and from 8:00 pm to 10:00 pm, which corresponded to the peaks of TB behavior. Sampling survey lasted for eight consecutive days until the study groups included 12 biter pigs (distributed in 8 pens), 12 bitten pigs (distributed in 7 pen) and 12 negative control pigs (distributed in 3 pens). The control pigs (non-bitten and non-biter) were selected in pens where TB behavior was absent according to the criteria defined in the ethogram. For the group of bitten pigs, all animals were examined based on a growing scale of tail damage scores from 0 to 3, where 0: no lesion; 1: presence of scratches; 2: presence of a moderate bleeding lesion; and 3: significant lesion with sometimes loss of a tail Sect. (58). The bitten and biter pigs were selected by video according to the ethogram. For confirmation of the pig tail status, an examination was made on-farm, every two days in the pens where TB behavior was present among the 32 pens in the study. All selected pigs were kept in their respective pens to continue their monitoring in their original and controlled environment. The behavior of all groups was studied simultaneously throughout a TB behavior episode that lasted four weeks, during which blood and fresh feces samples were taken for cortisol and IM analysis, respectively.
A positive control group consisting of 12 non-biter/non-bitten antibiotic-treated pigs was also formed. Blood samples for serum cortisol analysis were not taken. Chlortetracycline at 1210 ppm was added to their diet seven days prior to each sampling date to induce conformational changes in their IM. This group served as a positive control in our later 16S rRNA gene amplification, sequencing and analysis.
The average bite number of biter pigs was plotted on a weekly basis, w1, w2, w3 and w4.
Blood And Feces Sampling
Blood and feces were sampled in the biter, bitten and negative control groups, during two different periods, t0 and t1. No sampling was done on the AT group. At t0, pigs were 12 weeks of age, and TB behavior episode was beginning. Four weeks later, at t1, the TB behavior episode was considered finished.
Blood (approx. 5 ml) was collected from the jugular vein of snared pigs by a single experienced animal technician in the pen. Blood samples were kept at room temperature for 2 h to allow blood to clot prior to centrifugation (15 min at 1,000 x g). The serum was transferred to 1.5-mL Eppendorf tubes and stored at -80ºC pending cortisol concentration analysis. Serum cortisol concentration was determined using the Cortisol ELISA Kit (Pig) (Abnova, Taiwan) according to the manufacturer’s recommendation. The minimum detectable concentration of cortisol was 0.2 ng/ml.
Fresh fecal material was collected directly from the rectum of animals. A one-gram fraction (from five 200 mg subsamples from the same animal) was frozen immediately in liquid nitrogen and subsequently stored at -80ºC until DNA extraction.
DNA extraction, 16S ribosomal RNA gene amplification, sequencing and analysis
Total DNA was extracted from 500 µg of each feces sample according to Thibodeau et al., (2016) (59). Briefly, samples were put in tubes containing 0.1 mm glass beads. Bacteria were lysed with 500 ml of lysis buffer (500 mM Tris-HCl, 200 mM EDTA, 1% SDS; Fisher Scientific, Ottawa, ON, Canada) and a FastPrep-24 5G™ High Speed Homogenizer (mpbio, Santa Ana, CA, USA) for 2 cycles of 40 seconds at 6 m/s. Samples were kept on ice between cycles. DNA was purified using phenol: chloroform: isoamyl alcohol 25:24:1 (Sigma-Aldrich, St. Louis, MO, USA). The phenol traces were removed using chloroform: isoamyl alcohol 24:1 (Sigma-Aldrich). The DNA was precipitated in 90% ethanol for 24 h at -20 °C and resuspended in 1 mM Tris-HCl:0.1 M EDTA, pH 8.0. Negative controls without feces and positive controls with a known bacterial community (ZymoBIOMICS ™ Microbial Community Standard; Zymo Research, Irvine, CA, USA) were processed in parallel with the fecal samples. The purified DNAs were quantified using a QFX Fluorometer (DeNovix, Wilmington, DE, USA) with Qubit BR reagents (Fisher Scientific). DNA extracts were stored at -80ºC.
The hypervariable V4 region of the 16S rRNA gene was amplified for each sample by PCR using the primer pair 515FP1-CS1F ACACTGACGACATGGTTCTACAGTGCCAGCMGCCGCGGTAA and 806RP1-CS2R TACGGTAGCAGAGACTTGGTCTGGACTACHVGGGTWTCTAAT (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) [57]. The Platinum SuperFi PCR Master Mix (Invitrogen, Thermo Fisher Scientific) was used with 12.5 ng of DNA in a total reaction volume of 25 µl. The amplification was carried out with an initial denaturation step at 95ºC for 5 min, followed by 23 cycles at 98ºC for 30 sec, 55ºC for 30 sec and 72ºC for 180 sec, and a final elongation step at 72ºC for 10 min (60). Amplification was confirmed by gel electrophoresis and amplicons sent to the McGill University and Génome Québec Innovation Centre (Montreal, QC, Canada) for barcoding and subsequent Illumina Miseq sequencing (250 paired-ends).
All sequences were analyzed using Mothur software version 1.35.5 (61) according to Larivière-Gauthier et al., (2017) (62). Taxonomic assignment of the sequences was made using the Ribosomal Database Project (RDP; https://rdp.cme.msu.edu) (63). The sequences with 97% similarity (equivalent to species level) were grouped into operational taxonomic units (OTUs). Alpha-diversity (number of OTUs per sample, Shannon-even and inverse Simpson indices) of fecal samples from pigs of different groups, at t0 and t1were calculated in Mothur, using a subsample of 32,374 sequences, the lowest number of samples returned in all samples. For beta-diversity analysis, the distance between all samples was measured by the Yue & Clayton and the Jaccard indices using the same subsampling. Structure of the bacterial communities were visualized by a non-metric multidimensional scaling (NMDS) graph, and each combination of two pig groups, biter/negative, bitten/negative and biter/bitten, was compared at t0 and t1 by the analysis of molecular variance (AMOVA) (64). In addition, a Linear discriminant analysis (LDA) effect size (LEfSe) (65) was used to discover bacterial taxa significantly associated with each group at each sampling time.
Real-time Quantitative Pcr Of Specific Bacterial Populations
To validate and quantify results obtained from sequencing, a quantitative PCR (qPCR) targeting lactobacilli was performed on all samples (66). Standard curves were made from amplicons derived from the control strain Lactobacillus acidophilus ATCC 314. Each well contained 4 µl of Evagreen (MBI Montreal Biotech, Kirkland, QC, Canada), 0.6 µl of forward primer, 0.6 µl of reverse primer, 12.8 µl of water and 20 ng of DNA. The amplification was done in a LightCycler 96 real-time PCR (Roche Diagnostics, Mannheim, Germany) using the following program: 50ºC for 120 sec, 95ºC for ten minutes, 45 amplification cycles of 95ºC for 15 seconds and 60ºC for 60 seconds, and a final high-resolution melt analysis. The results, the number of gene copies, were expressed in log per ng of DNA.
Statistical analysis
A linear model for repeated measures with time as within-subject factor (SAS v9.3, Cary, NC, USA) was used to compare variations in TB across time. For the comparison of the different alpha-diversity indices, a linear model for repeated measures with time as within-subject factor and group as between-subject factor was used. For quantitative PCR results and serum cortisol levels, a linear model was used with group as factor. The AMOVA test was used to compare the beta-diversity between the different groups (comparison of two groups at a time) using Mothur (67). Statistical tests were made on comparisons between the following pig groups: biter/negative, bitten/negative and biter/bitten, at t0 and t1. The alpha level for these comparisons was adjusted downwards with the sequential method of Benjamini and Hochberg (1995) (68).