Use of serial pig body weights for genetic evaluation of daily gain

Gut microbiome composition differences among breeds impact feed efficiency in swine

BACKGROUND

Feed efficiency is a crucial parameter in swine production, given both its economic and environmental impact. The gut microbiota plays an essential role in nutrient digestibility and is, therefore, likely to affect feed efficiency. This study aimed to characterize feed efficiency, fatness traits, and gut microbiome composition in three major breeds of domesticated swine and investigate a possible link between feed efficiency and gut microbiota composition.

RESULTS

Average daily feed intake (ADFI), average daily gain (ADG), feed conversion ratio (FCR), residual feed intake (RFI), backfat, loin depth, and intramuscular fat of 615 pigs belonging to the Duroc (DR), Landrace (LR), and Large White (LW) breeds were measured. Gut microbiota composition was characterized by 16S rRNA gene sequencing. Orthogonal contrasts between paternal line (DR) and maternal lines (LR+LW) and between the two maternal lines (LR versus LW) were performed. Average daily feed intake and ADG were statistically different with DR having lower ADFI and ADG compared to LR and LW. Landrace and LW had a similar ADG and RFI, with higher ADFI and FCR for LW. Alpha diversity was higher in the fecal microbial communities of LR pigs than in those of DR and LW pigs for all time points considered. Duroc communities had significantly higher proportional representation of the Catenibacterium and Clostridium genera compared to LR and LW, while LR pigs had significantly higher proportions of Bacteroides than LW for all time points considered. Amplicon sequence variants from multiple genera (including Anaerovibrio, Bacteroides, Blautia, Clostridium, Dorea, Eubacterium, Faecalibacterium, Lactobacillus, Oscillibacter, and Ruminococcus) were found to be significantly associated with feed efficiency, regardless of the time point considered.

CONCLUSIONS

In this study, we characterized differences in the composition of the fecal microbiota of three commercially relevant breeds of swine, both over time and between breeds. Correlations between different microbiome compositions and feed efficiency were established. This suggests that the microbial community may contribute to shaping host productive parameters. Moreover, our study provides important insights into how the intestinal microbial community might influence host energy harvesting capacity. A deeper understanding of this process may allow us to modulate the gut microbiome in order to raise more efficient animals. Video Abstract.

Longitudinal analysis of weight showed little relationship with age at puberty in gilts

Continued selection for increased gilt body weight could negatively impact selection for age at puberty (AP) in gilts. The purpose of this study was to compare the genetic potential for growth to that for reducing age of puberty in swine. Females utilized (n = 1,079) were produced over a 6-yr period from a population developed to determine the impact of energy restrictions and genetic influences on sow development and longevity. From 120 to 235 d of age, BW was collected every 14 d and attainment of puberty tested. Age at puberty was defined as the first observed standing estrus in the presence of a mature boar. All females were genotyped with the Porcine SNP60K BeadChip and genotypes were used to construct a genomic relationship matrix. Univariate (AP), repeatability (BW), and random regression (BW; RR) models were fitted. Univariate analysis included the fixed effects of contemporary group (CG) and age at first boar exposure, and random direct additive and common litter effects. Repeatability analysis included the fixed effects of CG and random effects of direct additive, common litter, and permanent environmental (PE) effects. Random regression analysis included fixed effects of CG, and random direct additive, common litter, and PE effects. Proportion of phenotypic variation due to direct additive and common litter variance for AP were 0.33 and 0.06, respectively. Proportion of phenotypic variation due to direct additive, common litter, and PE variance estimated from the repeatability model for BW were 0.26, 0.11, and 0.40, respectively. Proportion of phenotypic variation due to direct additive, common litter, and PE variance estimated from the RR for BW ranged from (mean) 0.19 to 0.30 (0.27), 0.08 to 0.31 (0.19), and 0.42 to 0.62 (0.50), respectively. Direct additive correlations between test days for BW from RR ranged from 0.30 to 0.99. Rank correlations between estimated breeding values (EBV) for AP and BW from the RR were near zero across all age points ranging from -0.03 to 0.09. Rank correlations were higher (0.63) when BW was considered at the age at which puberty was reached. Correlations between AP and BW EBV from the repeatability model were low (-0.11). Growth appears to be less related to AP than previously reported, suggesting the need to either directly measure AP or investigate alternative indicator traits. Selecting gilts with most desirable BW EBV alone would not result in improvement in AP, at least in the current population.

The relationship between different measures of feed efficiency and feeding behavior traits in Duroc pigs

Utilization of feed in livestock species consists of a wide range of biological processes, and therefore, its efficiency can be expressed in various ways, including direct measurement, such as daily feed intake, as well as indicator measures, such as feeding behavior. Measuring feed efficiency is important to the swine industry, and its accuracy can be enhanced by using automated feeding systems, which record feed intake and associated feeding behavior of individual animals. Each automated feeder space is often shared among several pigs and therefore raises concerns about social interactions among pen mates with regard to feeding behavior. The study herein used a data set of 14,901 Duroc boars with individual records on feed intake, feeding behavior, and other off-test traits. These traits were modeled with and without the random spatial effect of Pen_Room, a concatenation of room and pen, or random social interaction among pen mates. The nonheritable spatial effect of common Pen-Room was observed for traits directly measuring feed intake and accounted for up to 13% of the total phenotypic variance in the average daily feeding rate. The social interaction effect explained larger proportions of phenotypic variation in all the traits studied, with the highest being 59% for ADFI in the group of feeding behaviors, 73% for residual feed intake (RFI; RFI4 and RFI6) in the feed efficiency traits, and 69% for intramuscular fat percentage in the off-test traits. After accounting for the social interaction effect, residual BW gain and RFI and BW gain (RIG) were found to have the heritability of 0.38 and 0.18, respectively, and had strong genetic correlations with growth and off-test traits. Feeding behavior traits were found to be moderately heritable, ranging from 0.14 (ADFI) to 0.52 (average daily occupation time), and some of them were strongly correlated with feed efficiency measures; for example, there was a genetic correlation of 0.88 between ADFI and RFI6. Our work suggested that accounting for the social common pen effect was important for estimating genetic parameters of traits recorded by the automated feeding system. Residual BW gain and RIG appeared to be two robust measures of feed efficiency. Feeding behavior measures are worth further investigation as indicators of feed efficiency.

How to improve breeding value prediction for feed conversion ratio in the case of incomplete longitudinal body weights

With the development of automatic self-feeders, repeated measurements of feed intake are becoming easier in an increasing number of species. However, the corresponding BW are not always recorded, and these missing values complicate the longitudinal analysis of the feed conversion ratio (FCR). Our aim was to evaluate the impact of missing BW data on estimations of the genetic parameters of FCR and ways to improve the estimations. On the basis of the missing BW profile in French Large White pigs (male pigs weighed weekly, females and castrated males weighed monthly), we compared 2 different ways of predicting missing BW, 1 using a Gompertz model and 1 using a linear interpolation. For the first part of the study, we used 17,398 weekly records of BW and feed intake recorded over 16 consecutive weeks in 1,222 growing male pigs. We performed a simulation study on this data set to mimic missing BW values according to the pattern of weekly proportions of incomplete BW data in females and castrated males. The FCR was then computed for each week using observed data (obser_FCR), data with missing BW (miss_FCR), data with BW predicted using a Gompertz model (Gomp_FCR), and data with BW predicted by linear interpolation (interp_FCR). Heritability (h) was estimated, and the EBV was predicted for each repeated FCR using a random regression model. In the second part of the study, the full data set (males with their complete BW records, castrated males and females with missing BW) was analyzed using the same methods (miss_FCR, Gomp_FCR, and interp_FCR). Results of the simulation study showed that h were overestimated in the case of missing BW and that predicting BW using a linear interpolation provided a more accurate estimation of h and of EBV than a Gompertz model. Over 100 simulations, the correlation between obser_EBV and interp_EBV, Gomp_EBV, and miss_EBV was 0.93 ± 0.02, 0.91 ± 0.01, and 0.79 ± 0.04, respectively. The heritabilities obtained with the full data set were quite similar for miss_FCR, Gomp_FCR, and interp_FCR. In conclusion, when the proportion of missing BW is high, genetic parameters of FCR are not well estimated. In French Large White pigs, in the growing period extending from d 65 to 168, prediction of missing BW using a Gompertz growth model slightly improved the estimations, but the linear interpolation improved the estimation to a greater extent. This result is due to the linear rather than sigmoidal increase in BW over the study period.

Genome-wide association study on legendre random regression coefficients for the growth and feed intake trajectory on Duroc Boars

Background

Feed intake and growth are economically important traits in swine production. Previous genome wide association studies (GWAS) have utilized average daily gain or daily feed intake to identify regions that impact growth and feed intake across time. The use of longitudinal models in GWAS studies, such as random regression, allows for SNPs having a heterogeneous effect across the trajectory to be characterized. The objective of this study is therefore to conduct a single step GWAS (ssGWAS) on the animal polynomial coefficients for feed intake and growth.

Results

Corrected daily feed intake (DFI_Adj) and average daily weight measurements (DBW_Avg) on 8981 (n = 525,240 observations) and 5643 (n = 283,607 observations) animals were utilized in a random regression model using Legendre polynomials (order = 2) and a relationship matrix that included genotyped and un-genotyped animals. A ssGWAS was conducted on the animal polynomials coefficients (intercept, linear and quadratic) for animals with genotypes (DFI_Adj: n = 855; DBW_Avg: n = 590). Regions were characterized based on the variance of 10-SNP sliding windows GEBV (WGEBV). A bootstrap analysis (n =1000) was conducted to declare significance. Heritability estimates for the traits trajectory ranged from 0.34-0.52 to 0.07-0.23 for DBW_Avg and DFI_Adj, respectively. Genetic correlations across age classes were large and positive for both DBW_Avg and DFI_Adj, albeit age classes at the beginning had a small to moderate genetic correlation with age classes towards the end of the trajectory for both traits. The WGEBV variance explained by significant regions (P < 0.001) for each polynomial coefficient ranged from 0.2-0.9 to 0.3-1.01 % for DBW_Avg and DFI_Adj, respectively. The WGEBV variance explained by significant regions for the trajectory was 1.54 and 1.95 % for DBW_Avg and DFI_Adj. Both traits identified candidate genes with functions related to metabolite and energy homeostasis, glucose and insulin signaling and behavior.

Conclusions

We have identified regions of the genome that have an impact on the intercept, linear and quadratic terms for DBW_Avg and DFI_Adj. These results provide preliminary evidence that individual growth and feed intake trajectories are impacted by different regions of the genome at different times.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0218-8) contains supplementary material, which is available to authorized users.

Feed intake, average daily gain, feed efficiency, and real-time ultrasound traits in Duroc pigs: I. Genetic parameter estimation and accuracy of genomic prediction

The efficiency of producing salable products in the pork industry is largely determined by costs associated with feed and by the amount and quality of lean meat produced. The objectives of this paper were 1) to explore heritability and genetic correlations for growth, feed efficiency, and real-time ultrasound traits using both pedigree and marker information and 2) to assess accuracy of genomic prediction for those traits using Bayes A prediction models in a Duroc terminal sire population. Body weight at birth (BW at birth) and weaning (BW at weaning) and real-time ultrasound traits, including back fat thickness (BF), muscle depth (MD), and intramuscular fat content (IMF), were collected on the basis of farm protocol. Individual feed intake and serial BW records of 1,563 boars obtained from feed intake recording equipment (FIRE; Osborne Industries Inc., Osborne, KS) were edited to obtain growth, feed intake, and feed efficiency traits, including ADG, ADFI, feed conversion ratio (FCR), and residual feed intake (RFI). Correspondingly, 1,047 boars were genotyped using the Illumina PorcineSNP60 BeadChip. The remaining 516 boars, as an independent sample, were genotyped with a low-density GGP-Porcine BeadChip and imputed to 60K. Magnitudes of heritability from pedigree analysis were moderate for growth, feed intake, and ultrasound traits (ranging from 0.44 ± 0.11 for ADG to 0.58 ± 0.09 for BF); heritability estimates were 0.32 ± 0.09 for FCR but only 0.10 ± 0.05 for RFI. Comparatively, heritability estimates using marker information by Bayes A models were about half of those from pedigree analysis, suggesting "missing heritability." Moderate positive genetic correlations between growth and feed intake (0.32 ± 0.05) and back fat (0.22 ± 0.04), as well as negative genetic correlations between growth and feed efficiency traits (-0.21 ± 0.08, -0.05 ± 0.07), indicate selection solely on growth traits may lead to an undesirable increase in feed intake, back fat, and reduced feed efficiency. Genetic correlations among growth, feed intake, and FCR assessed by a multiple-trait Bayes A model resulted in increased genetic correlation between ADG and ADFI, a negative correlation between ADFI and FCR, and a positive correlation between ADG and FCR. Accuracies of genomic prediction for the traits investigated, ranging from 9.4% for RFI to 36.5% for BF, were reported that might provide new insight into pig breeding and future selection programs using genomic information.