Divergent selection for feed efficiency in pigs altered the duodenum transcriptomic response to feed intake and its DNA methylation profiles.

Feed efficiency is a trait of interest in pigs as it contributes to lowering the ecological and economical costs of pig production. A divergent genetic selection experiment from a Large White pig population was performed for 10 generations, leading to pig lines with relatively low- (LRFI) and high- (HRFI) residual feed intake (RFI). Feeding behavior and metabolic differences have been previously reported between the two lines. We hypothesized that part of these differences could be related to differential sensing and absorption of nutrients in the proximal intestine. We investigated the duodenum transcriptome and DNA methylation profiles comparing overnight fasting with ad libitum feeding in LRFI and HRFI pigs (n = 24). We identified 1,106 differentially expressed genes between the two lines, notably affecting pathways of the transmembrane transport activity and related to mitosis or chromosome separation. The LRFI line showed a greater transcriptomic response to feed intake than the HRFI line. Feed intake affected genes from both anabolic and catabolic pathways in the pig duodenum, such as rRNA production and autophagy. Several nutrient transporter and tight junction genes were differentially expressed between lines and/or by short-term feed intake. We also identified 409 differentially methylated regions in the duodenum mucosa between the two lines, while this epigenetic mark was less affected by feeding. Our findings highlighted that the genetic selection for feed efficiency in pigs changed the transcriptome profiles of the duodenum, and notably its response to feed intake, suggesting key roles for this proximal gut segment in mechanisms underlying feed efficiency.NEW & NOTEWORTHY The duodenum is a key organ for the hunger/satiety loop and nutrient sensing. We investigated how the duodenum transcriptome and DNA methylation profiles are affected by feed intakes in pigs. We observed thousands of changes in gene expression levels between overnight-fasted and fed pigs in high-feed efficiency pig lines, but almost none in the related low-feed efficiency pig line.

The potential of microbiota information to better predict efficiency traits in growing pigs fed a conventional and a high-fiber diet.

BACKGROUND: Improving pigs' ability to digest diets with an increased dietary fiber content is a lever to improve feed efficiency and limit feed costs in pig production. The aim of this study was to determine whether information on the gut microbiota and host genetics can contribute to predict digestive efficiency (DE, i.e. digestibility coefficients of energy, organic matter, and nitrogen), feed efficiency (FE, i.e. feed conversion ratio and residual feed intake), average daily gain, and daily feed intake phenotypes. Data were available for 1082 pigs fed a conventional or high-fiber diet. Fecal samples were collected at 16 weeks, and DE was estimated using near­infrared spectrometry. A cross-validation approach was used to predict traits within the same diet, for the opposite diet, and for a combination of both diets, by implementing three models, i.e. with only genomic (Gen), only microbiota (Micro), and both genomic and microbiota information (Micro+Gen). The predictive ability with and without sharing common sires and breeding environment was also evaluated. Prediction accuracy of the phenotypes was calculated as the correlation between model prediction and phenotype adjusted for fixed effects. RESULTS: Prediction accuracies of the three models were low to moderate (< 0.47) for growth and FE traits and not significantly different between models. In contrast, for DE traits, prediction accuracies of model Gen were low (< 0.30) and those of models Micro and Micro+Gen were moderate to high (> 0.52). Prediction accuracies were not affected by the stratification of diets in the reference and validation sets and were in the same order of magnitude within the same diet, for the opposite diet, and for the combination of both diets. Prediction accuracies of the three models were significantly higher when pigs in the reference and validation populations shared common sires and breeding environment than when they did not (P < 0.001). CONCLUSIONS: The microbiota is a relevant source of information to predict DE regardless of the diet, but not to predict growth and FE traits for which prediction accuracies were similar to those obtained with genomic information only. Further analyses on larger datasets and more diverse diets should be carried out to complement and consolidate these results.

A genome-wide epistatic network underlies the molecular architecture of continuous color variation of body extremities.

Deciphering the molecular architecture of coat coloration for a better understanding of the biological mechanisms underlying pigmentation still remains a challenge. We took advantage of a rabbit French experimental population in which both a pattern and a gradient of coloration from white to brown segregated within the himalayan phenotype. The whole experimental design was genotyped using the high density Affymetrix® AxiomOrcun™ SNP Array and phenotyped into 6 different groups ordered from the lighter to the darker. Genome-wide association analyses pinpointed an oligogenic determinism, under recessive and additive inheritance, involving genes already known in melanogenesis (ASIP, KIT, MC1R, TYR), and likely processed pseudogenes linked to ribosomal function, RPS20 and RPS14. We also identified (i) gene-gene interactions through ASIP:MC1R affecting light cream/beige phenotypes while KIT:RPS responsible of dark chocolate/brown colors and (ii) a genome-wide epistatic network involving several others coloration genes such as POT1 or HPS5. Finally, we determined the recessive inheritance of the English spotting phenotype likely involving a copy number variation affecting at least the end of the coding sequence of the KIT gene. Our analyses of coloration as a continuous trait allowed us to go beyond much of the established knowledge through the detection of additional genes and gene-gene interactions that may contribute to the molecular architecture of the coloration phenotype.

Accounting for overlapping annotations in genomic prediction models of complex traits.

BACKGROUND: It is now widespread in livestock and plant breeding to use genotyping data to predict phenotypes with genomic prediction models. In parallel, genomic annotations related to a variety of traits are increasing in number and granularity, providing valuable insight into potentially important positions in the genome. The BayesRC model integrates this prior biological information by factorizing the genome according to disjoint annotation categories, in some cases enabling improved prediction of heritable traits. However, BayesRC is not adapted to cases where markers may have multiple annotations. RESULTS: We propose two novel Bayesian approaches to account for multi-annotated markers through a cumulative (BayesRC+) or preferential (BayesRC[Formula: see text]) model of the contribution of multiple annotation categories. We illustrate their performance on simulated data with various genetic architectures and types of annotations. We also explore their use on data from a backcross population of growing pigs in conjunction with annotations constructed using the PigQTLdb. In both simulated and real data, we observed a modest improvement in prediction quality with our models when used with informative annotations. In addition, our results show that BayesRC+ successfully prioritizes multi-annotated markers according to their posterior variance, while BayesRC[Formula: see text] provides a useful interpretation of informative annotations for multi-annotated markers. Finally, we explore several strategies for constructing annotations from a public database, highlighting the importance of careful consideration of this step. CONCLUSION: When used with annotations that are relevant to the trait under study, BayesRC[Formula: see text] and BayesRC+ allow for improved prediction and prioritization of multi-annotated markers, and can provide useful biological insight into the genetic architecture of traits.

Microbiability and microbiome-wide association analyses of feed efficiency and performance traits in pigs.

BACKGROUND: The objective of the present study was to investigate how variation in the faecal microbial composition is associated with variation in average daily gain (ADG), backfat thickness (BFT), daily feed intake (DFI), feed conversion ratio (FCR), and residual feed intake (RFI), using data from two experimental pig lines that were divergent for feed efficiency. Estimates of microbiability were obtained by a Bayesian approach using animal mixed models. Microbiome-wide association analyses (MWAS) were conducted by single-operational taxonomic units (OTU) regression and by back-solving solutions of best linear unbiased prediction using a microbiome covariance matrix. In addition, accuracy of microbiome predictions of phenotypes using the microbiome covariance matrix was evaluated. RESULTS: Estimates of heritability ranged from 0.31 ± 0.13 for FCR to 0.51 ± 0.10 for BFT. Estimates of microbiability were lower than those of heritability for all traits and were 0.11 ± 0.09 for RFI, 0.20 ± 0.11 for FCR, 0.04 ± 0.03 for DFI, 0.03 ± 0.03 for ADG, and 0.02 ± 0.03 for BFT. Bivariate analyses showed a high microbial correlation of 0.70 ± 0.34 between RFI and FCR. The two approaches used for MWAS showed similar results. Overall, eight OTU with significant or suggestive effects on the five traits were identified. They belonged to the genera and families that are mainly involved in producing short-chain fatty acids and digestive enzymes. Prediction accuracy of phenotypes using a full model including the genetic and microbiota components ranged from 0.60 ± 0.19 to 0.78 ± 0.05. Similar accuracies of predictions of the microbial component were observed using models that did or did not include an additive animal effect, suggesting no interaction with the genetic effect. CONCLUSIONS: Our results showed substantial associations of the faecal microbiome with feed efficiency related traits but negligible effects with growth traits. Microbiome data incorporated as a covariance matrix can be used to predict phenotypes of animals that do not (yet) have phenotypic information. Connecting breeding environment between training sets and predicted populations could be necessary to obtain reliable microbiome predictions.

Quality of breeding value predictions from longitudinal analyses, with application to residual feed intake in pigs.

BACKGROUND: An important goal in animal breeding is to improve longitudinal traits. The objective of this study was to explore for longitudinal residual feed intake (RFI) data, which estimated breeding value (EBV), or combination of EBV, to use in a breeding program. Linear combinations of EBV (summarized breeding values, SBV) or phenotypes (summarized phenotypes) derived from the eigenvectors of the genetic covariance matrix over time were considered, and the linear regression method (LR method) was used to facilitate the evaluation of their prediction accuracy. RESULTS: Weekly feed intake, average daily gain, metabolic body weight, and backfat thickness measured on 2435 growing French Large White pigs over a 10-week period were analysed using a random regression model. In this population, the 544 dams of the phenotyped animals were genotyped. These dams did not have own phenotypes. The quality of the predictions of SBV and breeding values from summarized phenotypes of these females was evaluated. On average, predictions of SBV at the time of selection were unbiased, slightly over-dispersed and less accurate than those obtained with additional phenotypic information. The use of genomic information did not improve the quality of predictions. The use of summarized instead of longitudinal phenotypes resulted in predictions of breeding values of similar quality. CONCLUSIONS: For practical selection on longitudinal data, the results obtained with this specific design suggest that the use of summarized phenotypes could facilitate routine genetic evaluation of longitudinal traits.

Analysis of the causal structure of traits involved in sow lactation feed efficiency.

BACKGROUND: Feed efficiency during lactation involves a set of phenotypic traits that form a complex system, with some traits exerting causal effects on the others. Information regarding such interrelationships can be used to predict the effect of external interventions on the system, and ultimately to optimize management practices and multi-trait selection strategies. Structural equation models can be used to infer the magnitude of the different causes of such interrelationships. The causal network necessary to fit structural equation models can be inferred using the inductive causation (IC) algorithm. By implementing these statistical tools, we inferred the causal association between the main energy sources and sinks involved in sow lactation feed efficiency for the first time, i.e., daily lactation feed intake (dLFI) in kg/day, daily sow weight balance (dSWB) in kg/day, daily litter weight gain (dLWG) in kg/day, daily back fat thickness balance (dBFTB) in mm/day, and sow metabolic body weight (SMBW) in kg0.75. Then, we tested several selection strategies based on selection indices, with or without dLFI records, to improve sow efficiency during lactation. RESULTS: The IC algorithm using 95% highest posterior density (HPD95%) intervals resulted in a fully directed acyclic graph, in which dLFI and dLWG affected dSWB, the posterior mean of the corresponding structural coefficients (PMλ) being 0.12 and - 0.03, respectively. In turn, dSWB influenced dBFTB and SMBW, with PMλ equal to 0.70 and - 1.22, respectively. Multiple indirect effects contributed to the variances and covariances among the analyzed traits, with the most relevant indirect effects being those involved in the association between dSWB and dBFTB and between dSWB and SMBW. Selection strategies with or without phenotypic information on dLFI, or that hold this trait constant, led to the same pattern and similar responses in dLFI, dSWB, and dLWG. CONCLUSIONS: Selection based on an index including only dBFTB and dLWG records can reduce dLFI, keep dSWB constant or increase it, and increase dLWG. However, a favorable response for all three traits is probably not achievable. Holding the amount of feed provided to the sows constant did not offer an advantage in terms of response over the other strategies.

Gut microbiota and host genetics contribute to the phenotypic variation of digestive and feed efficiency traits in growing pigs fed a conventional and a high fiber diet.

BACKGROUND: Breeding pigs that can efficiently digest alternative diets with increased fiber content is a viable strategy to mitigate the feed cost in pig production. This study aimed at determining the contribution of the gut microbiota and host genetics to the phenotypic variability of digestive efficiency (DE) traits, such as digestibility coefficients of energy, organic matter and nitrogen, feed efficiency (FE) traits (feed conversion ratio and residual feed intake) and growth traits (average daily gain and daily feed intake). Data were available for 791 pigs fed a conventional diet and 735 of their full-sibs fed a high-fiber diet. Fecal samples were collected at 16 weeks of age to sequence the V3-V4 regions of the 16S ribosomal RNA gene and predict DE with near-infrared spectrometry. The proportions of phenotypic variance explained by the microbiota (microbiability) were estimated under three OTU filtering scenarios. Then, microbiability and heritability were estimated independently (models Micro and Gen) and jointly (model Micro+Gen) using a Bayesian approach for all traits. Breeding values were estimated in models Gen and Micro+Gen. RESULTS: Differences in microbiability estimates were significant between the two extreme filtering scenarios (14,366 and 803 OTU) within diets, but only for all DE. With the intermediate filtering scenario (2399 OTU) and for DE, microbiability was higher (> 0.44) than heritability (< 0.32) under both diets. For two of the DE traits, microbiability was significantly higher under the high-fiber diet (0.67 ± 0.06 and 0.68 ± 0.06) than under the conventional diet (0.44 ± 0.06). For growth and FE, heritability was higher (from 0.26 ± 0.06 to 0.44 ± 0.07) than microbiability (from 0.17 ± 0.05 to 0.35 ± 0.06). Microbiability and heritability estimates obtained with the Micro+Gen model did not significantly differ from those with the Micro and Gen models for all traits. Finally, based on their estimated breeding values, pigs ranked differently between the Gen and Micro+Gen models, only for the DE traits under both diets. CONCLUSIONS: The microbiota explained a significant proportion of the phenotypic variance of the DE traits, which was even larger than that explained by the host genetics. Thus, the use of microbiota information could improve the selection of DE traits, and to a lesser extent, of growth and FE traits. In addition, our results show that, at least for DE traits, filtering OTU is an important step and influences the microbiability.

Physiological response to the weaning in two pig lines divergently selected for residual feed intake.

Breeding efficient pigs is a way to reduce dietary costs and environmental waste. However, optimization of feed efficiency must not be linked to a decrease of the ability of animals to cope with stress, such as the weaning. This study characterizes the response after weaning of pigs from two lines divergently selected for residual feed intake (RFI) during growth. Animals of the low (L) RFI line are more efficient than animals from the high (H) RFI line. Thirty-six piglets from each line weaned at 28 days of age were individually housed and fed a conventional dietary sequence. Their performance, behaviour, health and oxidative status, immune and nutritional parameters were followed during three weeks. Daily feed intake and growth rate of pigs from the LRFI line were 35% and 40% lower compared with HRFI (p < 0.001). Pigs from the LRFI-line had lower total tract apparent digestibility (-6% for OM) and suffered more from undernutrition with a 167 and 55% higher plasmatic concentration of NEFA and urea compared with HRFI (p < 0.01). In the first week after the weaning, they had more diarrhoea and had a higher inflammatory status with concentration of haptoglobin 52% higher (p < 0.001). These piglets then seemed to adapt to the weaning conditions and to recover during the second and third weeks. Both lines had similar zootechnical performance and physiological characteristics at the end of the post-weaning period. To conclude, the physiological responses to the weaning differed between lines. Pigs from the LRFI line, selected for greater feed efficiency, were more sensitive to the weaning stress. They were also more resilient as they finally adapted to the new condition and recovered to show similar performance results as pigs of the HRFI line.

Analysis of merged whole blood transcriptomic datasets to identify circulating molecular biomarkers of feed efficiency in growing pigs.

BACKGROUND: Improving feed efficiency (FE) is an important goal due to its economic and environmental significance for farm animal production. The FE phenotype is complex and based on the measurements of the individual feed consumption and average daily gain during a test period, which is costly and time-consuming. The identification of reliable predictors of FE is a strategy to reduce phenotyping efforts. RESULTS: Gene expression data of the whole blood from three independent experiments were combined and analyzed by machine learning algorithms to propose molecular biomarkers of FE traits in growing pigs. These datasets included Large White pigs from two lines divergently selected for residual feed intake (RFI), a measure of net FE, and in which individual feed conversion ratio (FCR) and blood microarray data were available. Merging the three datasets allowed considering FCR values (Mean = 2.85; Min = 1.92; Max = 5.00) for a total of n = 148 pigs, with a large range of body weight (15 to 115 kg) and different test period duration (2 to 9 weeks). Random forest (RF) and gradient tree boosting (GTB) were applied on the whole blood transcripts (26,687 annotated molecular probes) to identify the most important variables for binary classification on RFI groups and a quantitative prediction of FCR, respectively. The dataset was split into learning (n = 74) and validation sets (n = 74). With iterative steps for variable selection, about three hundred's (328 to 391) molecular probes participating in various biological pathways, were identified as important predictors of RFI or FCR. With the GTB algorithm, simpler models were proposed combining 34 expressed unique genes to classify pigs into RFI groups (100% of success), and 25 expressed unique genes to predict FCR values (R2 = 0.80, RMSE = 8%). The accuracy performance of RF models was slightly lower in classification and markedly lower in regression. CONCLUSION: From small subsets of genes expressed in the whole blood, it is possible to predict the binary class and the individual value of feed efficiency. These predictive models offer good perspectives to identify animals with higher feed efficiency in precision farming applications.

New residual feed intake criterion for longitudinal data.

BACKGROUND: Residual feed intake (RFI) is one measure of feed efficiency, which is usually obtained by multiple regression of feed intake (FI) on measures of production, body weight gain and tissue composition. If phenotypic regression is used, the resulting RFI is generally not genetically independent of production traits, whereas if RFI is computed using genetic regression coefficients, RFI and production traits are independent at the genetic level. The corresponding regression coefficients can be easily derived from the result of a multiple trait model that includes FI and production traits. However, this approach is difficult to apply in the case of multiple repeated measurements of FI and production traits. To overcome this difficulty, we used a structured antedependence approach to account for the longitudinality of the data with a phenotypic regression model or with different genetic and environmental regression coefficients [multi- structured antedependence model (SAD) regression model]. RESULTS: After demonstrating the properties of RFI obtained by the multi-SAD regression model, we applied the two models to FI and production traits that were recorded for 2435 French Large White pigs over a 10-week period. Heritability estimates were moderate with both models. With the multi-SAD regression model, heritability estimates were quite stable over time, ranging from 0.14 ± 0.04 to 0.16 ± 0.05, while heritability estimates showed a U-shaped profile with the phenotypic regression model (ranging from 0.19 ± 0.06 to 0.28 ± 0.06). Estimates of genetic correlations between RFI at different time points followed the same pattern for the two models but higher estimates were obtained with the phenotypic regression model. Estimates of breeding values that can be used for selection were obtained by eigen-decomposition of the genetic covariance matrix. Correlations between these estimated breeding values obtained with the two models ranged from 0.66 to 0.83. CONCLUSIONS: The multi-SAD model is preferred for the genetic analysis of longitudinal RFI because, compared to the phenotypic regression model, it provides RFI that are genetically independent of production traits at all time points. Furthermore, it can be applied even when production records are missing at certain time points.

Identification of genomic regions affecting production traits in pigs divergently selected for feed efficiency.

BACKGROUND: Feed efficiency is a major driver of the sustainability of pig production systems. Understanding the biological mechanisms that underlie these agronomic traits is an important issue for environment questions and farms' economy. This study aimed at identifying genomic regions that affect residual feed intake (RFI) and other production traits in two pig lines divergently selected for RFI during nine generations (LRFI, low RFI; HRFI, high RFI). RESULTS: We built a whole dataset of 570,447 single nucleotide polymorphisms (SNPs) in 2426 pigs with records for 24 production traits after both imputation and prediction of genotypes using pedigree information. Genome-wide association studies (GWAS) were performed including both lines (global-GWAS) or each line independently (LRFI-GWAS and HRFI-GWAS). Forty-five chromosomal regions were detected in the global-GWAS, whereas 28 and 42 regions were detected in the HRFI-GWAS and LRFI-GWAS, respectively. Among these 45 regions, only 13 were shared between at least two analyses, and only one was common between the three GWAS but it affects different traits. Among the five quantitative trait loci (QTL) detected for RFI, two were close to QTL for meat quality traits and two pinpointed novel genomic regions that harbor candidate genes involved in cell proliferation and differentiation processes of gastrointestinal tissues or in lipid metabolism-related signaling pathways. In most cases, different QTL regions were detected between the three designs, which suggests a strong impact of the dataset structure on the detection power and could be due to the changes in allelic frequencies during the establishment of lines. CONCLUSIONS: In addition to efficiently detecting known and new QTL regions for feed efficiency, the combination of GWAS carried out per line or simultaneously using all individuals highlighted chromosomal regions that affect production traits and presented significant changes in allelic frequencies across generations. Further analyses are needed to estimate whether these regions correspond to traces of selection or result from genetic drift.

Genetics of digestive efficiency in growing pigs fed a conventional or a high-fibre diet.

The use of diets with increased dietary fibre content (HF) from alternative feedstuffs is a solution to limit the impact of increased feed costs on pig production. This study aimed at determining the impact of an alternative HF diet on pig digestibility and at estimating genetic parameters of this trait. Digestibility coefficients (DC) of energy, organic matter and nitrogen were predicted from faecal samples analysed with near infrared spectrometry for 1,242 samples, and it represented 654 Large White pigs fed a conventional (CO) diet and 588 fed a HF diet. Growth and feed efficiency traits, carcass composition and meat quality traits were recorded. Pigs fed the HF diet had significantly lower DC than pigs fed the CO diet (-4.5 to 6.0 points). The DC were moderately to highly heritable (about 0.26 ± 0.12 and 0.54 ± 0.15 in the CO and the HF diet, respectively). Genetic correlations were favourable with feed conversion ratio, daily feed intake and residual feed intake, but unfavourable with average daily gain (ADG) and carcass yield (CY). To conclude, DC could be an interesting trait to include in future breeding objectives if pigs were fed diet with HF diets, but adverse genetic trends with ADG and CY would have to be taken into account.

Genetic relationships between feed efficiency and gut microbiome in pig lines selected for residual feed intake.

This study aimed to evaluate the genetic relationship between faecal microbial composition and five feed efficiency (FE) and production traits, residual feed intake (RFI), feed conversion ratio (FCR), daily feed intake (DFI), average daily gain (ADG) and backfat thickness (BFT). A total of 588 samples from two experimental pig lines developed by divergent selection for RFI were sequenced for the 16 rRNA hypervariable V3-V4 region. The 75 genera with less than 20% zero values (97% of the counts) and two α-diversity indexes were analysed. Line comparison of the microbiota traits and estimations of heritability (h2 ) and genetic correlations (rg ) were analysed. A non-metric multidimensional scaling showed line differences between genera. The α-diversity indexes were higher in the LRFI line than in the HRFI line (p < .01), with h2 estimates of 0.19 ± 0.08 (Shannon) and 0.12 ± 0.06 (Simpson). Forty-eight genera had a significant h2 (>0.125). The rg of the α-diversities indexes with production traits were negative. Some rg of genera belonging to the Lachnospiraceae, Ruminococcaceae, Prevotellaceae, Lactobacillaceae, Streptococcaceae, Rikenellaceae and Desulfovibrionaceae families significantly differed from zero (p < .05) with FE traits, RFI (3), DFI (7) and BFT (11). These results suggest that a sizable part of the variability of the gut microbial community is under genetic control and has genetic relationships with FE, including diversity indicators. It offers promising perspectives for selection for feed efficiency using gut microbiome composition in pigs.

The impact of training on data from genetically-related lines on the accuracy of genomic predictions for feed efficiency traits in pigs.

BACKGROUND: Most genomic predictions use a unique population that is split into a training and a validation set. However, genomic prediction using genetically heterogeneous training sets could provide more flexibility when constructing the training sets in small populations. The aim of our study was to investigate the potential of genomic prediction of feed efficiency related traits using training sets that combine animals from two different, but genetically-related lines. We compared realized prediction accuracy and prediction bias for different training set compositions for five production traits. RESULTS: Genomic breeding values (GEBV) were predicted using the single-step genomic best linear unbiased prediction method in six scenarios applied iteratively to two genetically-related lines (i.e. 12 scenarios). The objective for all scenarios was to predict GEBV of pigs in the last three generations (~ 400 pigs, G7 to G9) of a given line. For each line, a control scenario was set up with a training set that included only animals from that line (target line). For all traits, adding more animals from the other line to the training set did not increase prediction accuracy compared to the control scenario. A small decrease in prediction accuracies was found for average daily gain, backfat thickness, and daily feed intake as the number of animals from the target line decreased in the training set. Including more animals from the other line did not decrease prediction accuracy for feed conversion ratio and residual feed intake, which were both highly affected by selection within lines. However, prediction biases were systematic for these cases and might be reduced with bivariate analyses. CONCLUSIONS: Our results show that genomic prediction using a training set that includes animals from genetically-related lines can be as accurate as genomic prediction using a training set from the target population. With combined reference sets, accuracy increased for traits that were highly affected by selection. Our results provide insights into the design of reference populations, especially to initiate genomic selection in small-sized lines, for which the number of historical samples is small and that are developed simultaneously. This applies especially to poultry and pig breeding and to other crossbreeding schemes.

Inclusive inheritance for residual feed intake in pigs and rabbits.

Non-genetic information (epigenetic, microbiota, behaviour) that results in different phenotypes in animals can be transmitted from one generation to the next and thus is potentially involved in the inheritance of traits. However, in livestock species, animals are selected based on genetic inheritance only. The objective of the present study was to determine whether non-genetic inherited effects play a role in the inheritance of residual feed intake (RFI) in two species: pigs and rabbits. If so, the path coefficients of the information transmitted from sire and dam to offspring would differ from the expected transmission factor of 0.5 that occurs if inherited information is of genetic origin only. Two pigs (pig1, pig2) and two rabbits (rabbit1, rabbit2) datasets were used in this study (1,603, 3,901, 5,213 and 4,584 records, respectively). The test of the path coefficients to 0.5 was performed for each dataset using likelihood ratio tests (null model: transmissibility model with both path coefficients equal to 0.5, full model: unconstrained transmissibility model). The path coefficients differed significantly from 0.5 for one of the pig datasets (pig2). Although not significant, we observed, as a general trend, that sire path coefficients of transmission were lower than dam path coefficients in three of the datasets (0.46 vs 0.53 for pig1, 0.39 vs 0.44 for pig2 and 0.38 vs 0.50 for rabbit1). These results suggest that phenomena other than genetic sources of inheritance explain the phenotypic resemblance between relatives for RFI, with a higher transmission from the dam's side than from the sire's side.

Prospects for sustainability of pig production in relation to climate change and novel feed resources.

Pig production systems provide multiple benefits to humans. However, the global increase in meat consumption has profound consequences for our earth. This perspective describes two alternative scenarios for improving the sustainability of future pig production systems. The first scenario is a high input-high output system based on sustainable intensification, maximizing animal protein production efficiency on a limited land surface at the same time as minimizing environmental impacts. The second scenario is a reduced input-reduced output system based on selecting animals that are more robust to climate change and are better adapted to transform low quality feed (local feeds, feedstuff co-products, food waste) into meat. However, in contrast to the first scenario, the latter scenario results in reduced predicted yields, reduced production efficiency and possibly increased costs to the consumer. National evaluation of the availability of local feed and feedstuff co-product alternatives, determination of limits to feed sourced from international markets, available land for crop and livestock production, desired production levels, and a willingness to politically enforce policies through subsidies and/or penalties are some of the considerations to combine these two scenarios. Given future novel sustainable alternatives to livestock animal protein, it may become reasonable to move towards an added general premium price on 'protein from livestock animals' to the benefit of promoting higher incomes to farmers at the same time as covering the extra costs of, politically enforced, welfare of livestock animals in sustainable production systems. © 2020 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Investigation of muscle transcriptomes using gradient boosting machine learning identifies molecular predictors of feed efficiency in growing pigs.

BACKGROUND: Improving feed efficiency (FE) is a major challenge in pig production. This complex trait is characterized by a high variability. Therefore, the identification of predictors of FE may be a relevant strategy to reduce phenotyping efforts in breeding and selection programs. The aim of this study was to investigate the suitability of expressed muscle genes in prediction of FE traits in growing pigs. The approach considered different transcriptomics experiments to cover a large range of FE values and identify reliable predictors. RESULTS: Microarrays data were obtained from longissimus muscles of two lines divergently selected for residual feed intake (RFI). Pigs (n = 71) from three experiments belonged to generations 6 to 8 of selection, were fed either a diet with a standard composition or a diet rich in fiber and lipids, received feed ad libitum or at restricted level, and weighed between 80 and 115 kg at slaughter. For each pig, breeding value for RFI was estimated (RFI-BV), and feed conversion ratio (FCR) and energy-based feed conversion ratio (FCRe) were calculated during the test periods. Gradient boosting algorithms were used on the merged muscle transcriptomes to identify very important predictors of FE traits. About 20,405 annotated molecular probes were commonly expressed in longissimus muscle across experiments. Six to 267 expressed muscle genes covering a variety of biological processes were found as important predictors for RFI-BV (R2 = 0.63-0.65), FCR (R2 = 0.61-0.70) and FCRe (R2 = 0.49-0.52). The error of prediction was less than 8% for FCR. Altogether, 56 predictors were common to RFI-BV and FCR. Expression levels of 24 target genes were further measured by qPCR. Linear regression confirmed the good accuracy of combining mRNA levels of these genes to fit FE traits (RFI-BV: R2 = 0.73, FRC: R2 = 0.76; FCRe: R2 = 0.75). Stepwise regression procedure highlighted 10 genes (FKBP5, MUM1, AKAP12, FYN, TMED3, PHKB, TGF, SOCS6, ILR4, and FRAS1) in a linear combination predicting FCR and FCRe. In addition, FKBP5 and expression levels of five other genes (IGF2, SERINC3, CSRNP3, EZR and RPL16) significantly contributed to RFI-BV. CONCLUSION: It was possible to identify few genes expressed in muscle that might be reliable predictors of feed efficiency.

Estimating direct genetic and maternal effects affecting rabbit growth and feed efficiency with a factorial design.

The aim of this experiment was to evaluate the significance of neonatal environment on feed efficiency. For that purpose, rabbits from a line selected for residual feed intake (RFI) during 10 generations (G10 kits) were cross-fostered with non-selected control does (i.e., G0 line), and reciprocally. In parallel, sibs were fostered by mothers from their original line. Nine hundred animals were raised in individual (N = 456) or collective (N = 320) cages. Traits analysed in this study were body weight at 32 days and at 63 days, average daily gain (ADG), feed intake between weaning and 63 days (FI), feed conversion ratio (FCR) and RFI. The maternal environment offered by does from the line selected for RFI deteriorated the FCR of the kits, independently of their line of origin, during fattening (+0.08 ± 0.02) compared to FCR of kits nursed by G0 does. The line, the type of housing and the batch were significant effects for all the measured traits: G10 kits were lighter than their G0 counterparts at 32 days (-82.9 ± 9 g, p < 0.0001) and at 63 days (-161 ± 16 g, p < 0.0001). They also had a lower ADG (-2.36 ± 0.36 g/day, p < 0.0001), RFI (-521 ± 24 g/day, p < 0.0001) and a lower FI (-855 ± 31 g, p < 0.0001), resulting in a more desirable feed efficiency (FCR: -0.35 ± 0.02). There was no significant difference in the contrast of G10 and G0 performances between collective and individual/digestive cages (p > 0.22): -2.35 g/day versus 2.94 g/day for ADG, -0.39 versus -0.40 for FCR, -577 g versus -565 g for RFI and -879 g versus -859 g for FI, respectively). Thus, no genotype-by-environment (housing) interaction is expected at the commercial level, that is, no re-ranking of the animals due to collective housing.

A transcriptome multi-tissue analysis identifies biological pathways and genes associated with variations in feed efficiency of growing pigs.

BACKGROUND: Animal's efficiency in converting feed into lean gain is a critical issue for the profitability of meat industries. This study aimed to describe shared and specific molecular responses in different tissues of pigs divergently selected over eight generations for residual feed intake (RFI). RESULTS: Pigs from the low RFI line had an improved gain-to-feed ratio during the test period and displayed higher leanness but similar adiposity when compared with pigs from the high RFI line at 132 days of age. Transcriptomics data were generated from longissimus muscle, liver and two adipose tissues using a porcine microarray and analyzed for the line effect (n = 24 pigs per line). The most apparent effect of the line was seen in muscle, whereas subcutaneous adipose tissue was the less affected tissue. Molecular data were analyzed by bioinformatics and subjected to multidimensional statistics to identify common biological processes across tissues and key genes participating to differences in the genetics of feed efficiency. Immune response, response to oxidative stress and protein metabolism were the main biological pathways shared by the four tissues that distinguished pigs from the low or high RFI lines. Many immune genes were under-expressed in the four tissues of the most efficient pigs. The main genes contributing to difference between pigs from the low vs high RFI lines were CD40, CTSC and NTN1. Different genes associated with energy use were modulated in a tissue-specific manner between the two lines. The gene expression program related to glycogen utilization was specifically up-regulated in muscle of pigs from the low RFI line (more efficient). Genes involved in fatty acid oxidation were down-regulated in muscle but were promoted in adipose tissues of the same pigs when compared with pigs from the high RFI line (less efficient). This underlined opposite line-associated strategies for energy use in skeletal muscle and adipose tissue. Genes related to cholesterol synthesis and efflux in liver and perirenal fat were also differentially regulated in pigs from the low vs high RFI lines. CONCLUSIONS: Non-productive functions such as immunity, defense against pathogens and oxidative stress contribute likely to inter-individual variations in feed efficiency.