Early-life establishment of the swine gut microbiome and impact on host phenotypes

Early bacterial colonization and succession within the gastrointestinal tract has been suggested to be crucial in the establishment of specific microbiota composition and the shaping of host phenotype. Here, the composition and dynamics of faecal microbiomes were studied for 31 healthy piglets across five age strata (days 14, 36, 48, 60 and 70 after birth) together with their mothers. Faecal microbiome composition was assessed by 16S rRNA gene 454-pyrosequencing. Bacteroidetes and Firmicutes were the predominant phyla present at each age. For all piglets, luminal secretory IgA concentration was measured at day 70, and body weight was recorded until day 70. The microbiota of suckling piglets was mainly represented by Bacteroides, Oscillibacter, Escherichia/Shigella, Lactobacillus and unclassified Ruminococcaceae genera. This pattern contrasted with that of Acetivibrio, Dialister, Oribacterium, Succinivibrio and Prevotella genera, which appeared increased after weaning. Lactobacillus fermentum might be vertically transferred via breast milk or faeces. The microbiota composition coevolved with their hosts towards two different clusters after weaning, primarily distinguished by unclassified Ruminococcaceae and Prevotella abundances. Prevotella was positively correlated with luminal secretory IgA concentrations, and body weight. Our study opens up new possibilities for health and feed efficiency manipulation via genetic selection and nutrition in the agricultural domain.

Genetic parameters for residual feed intake in growing pigs, with emphasis on genetic relationships with carcass and meat quality traits

Data were collected over the first 4 generations of a divergent selection experiment for residual feed intake of Large White pigs having ad libitum access to feed. This data set was used to obtain estimates of heritability for residual feed intake and genetic correlations (r(a)) between this trait and growth, carcass, and meat quality traits. Individual feed intake of group-housed animals was measured by single-space electronic feeders. Upward and downward selection lines were maintained contemporarily, with 6 boars and 35 to 40 sows per line and generation. Numbers of records were 793 for residual feed intake (RFI1) of boar candidates for selection issued from first-parity (P1) litters and tested over a fixed BW range (35 to 95 kg) and 657 for residual feed intake (RFI2) and growth, carcass, and meat quality traits of castrated males and females issued from second-parity (P2) litters and tested from 28 to 107 kg of BW. Variance and covariance components were estimated using REML methodology applied to a series of multitrait animal models, which always included the criterion for selection as 1 of the traits. Estimates of heritability for RFI1 and RFI2 were 0.14 +/- 0.03 and 0.24 +/- 0.03, respectively, whereas the estimate of r(a) between the 2 traits was 0.91 +/- 0.08. Estimates of r(a) indicated that selection for low residual feed intake has the potential to improve feed conversion ratio and reduce daily feed intake, with minimal correlated effect for ADG of P2 animals. Estimates of r(a) between RFI2 and body composition traits of P2 animals were positive for traits related to the amount of fat depots (r(a) = 0.44 +/- 0.16 for carcass backfat thickness) and negative for carcass lean meat content (r(a) = -0.55 +/- 0.14). There was a tendency for a negative genetic correlation between RFI2 and carcass dressing percent (r(a) = -0.36 +/- 0.21). Moreover, selection for low residual feed intake is expected, through lower ultimate pH and lighter color, to decrease pork quality (r(a) = 0.77 +/- 0.14 between RFI2 and a meat quality index intended to predict the ratio of the weight of ham after curing and cooking to the weight of defatted and boneless fresh ham).

Between-breed variability of stillbirth and its relationship with sow and piglet characteristics

Litter characteristics at birth were recorded in 4 genetic types of sows with differing maternal abilities. Eighty-two litters from F(1) Duroc x Large White sows, 651 litters from Large White sows, 63 litters from Meishan sows, and 173 litters from Laconie sows were considered. Statistical models included random effects of sow, litter, or both; fixed effects of sow genetic type, parity, birth assistance, and piglet sex, as well as gestation length, farrowing duration, piglet birth weight, and litter size as linear covariates. The quadratic components of the last 2 factors were also considered. For statistical analyses, GLM were first considered, assuming a binomial distribution of stillbirth. Hierarchical models were also fitted to the data to take into account correlations among piglets from the same litter. Model selection was performed based on deviance and deviance information criterion. Finally, standard and robust generalized estimating equations (GEE) procedures were applied to quantify the importance of each effect on a piglet's probability of stillbirth. The 5 most important factors involved were, in decreasing order (contribution of each effect to variance reduction): difference between piglet birth weight and the litter mean (2.36%), individual birth weight (2.25%), piglet sex (1.01%), farrowing duration (0.99%), and sow genetic type (0.94%). Probability of stillbirth was greater for lighter piglets, for male piglets, and for piglets from small or very large litters. Probability of stillbirth increased with sow parity number and with farrowing duration. Piglets born from Meishan sows had a lower risk of stillbirth (P < 0.0001) and were little affected by the sources of variation mentioned above compared with the 3 other sow genetic types. Standard and robust GEE approaches gave similar results despite some disequilibrium in the data set structure highlighted with the robust GEE approach.

Muscle characteristics and meat quality traits are affected by divergent selection on residual feed intake in pigs

Residual feed intake (RFI) is defined as the difference between the observed feed intake and that expected based on requirements for maintenance and production. A divergent selection was conducted during 4 generations in Large White male pigs to produce low and high RFI lines. The present study aims at determining the influence of this selection on biochemical and histological traits of skeletal muscle, and relating these changes to correlated effects on growth, carcass composition, and meat quality traits. At 8 d preslaughter, biopsies from the LM were taken in the fed state on 14 females from each RFI line fed ad libitum. Animals were slaughtered at 107.8 ± 8.0 kg of BW without any previous fasting. Samples of LM, semimembranosus (SM), biceps femoris (BFM), and rhomboideus muscles were taken at both 30 min and 24 h postmortem. Myofiber typing was only assessed in LM. Low RFI pigs ("efficient") had leaner carcasses with greater muscle content (P < 0.001), less backfat thickness (P < 0.001), and less intramuscular fat content in all 4 muscles (P < 0.01 to P = 0.04). Their greater muscle content was associated with hypertrophy of all fast-twitch fibers. Glycogen content in all glycolytic muscles (i.e., LM, SM and BFM), was greater in low than high RFI pigs. The greater accumulation of glycogen in LM of low RFI pigs was specifically located in the fast-twitch glycolytic IIBW fibers, which correspond to fibers containing IIb, IIb + IIx, or IIx myosin heavy chains. The difference in muscle glycogen content between RFI line pigs was more significant in the living animals (P = 0.0003) than at 30 min postmortem (P = 0.08). This was associated with a decreased ultimate pH (P = 0.001), and greater lightness of color (P = 0.002) and drip loss (P = 0.04) in LM of low than high RFI line pigs, suggesting that selection for reduced RFI may impair some meat quality traits, such as water-holding capacity. Pigs from the low RFI line exhibited a greater (P = 0.02) percentage of IIBW fibers in LM and tended (P < 0.10) to have less lipid β-oxidative capacity in LM, SM, and BFM. In contrast, no difference (P > 0.10) between lines was found for citrate synthase and lactate dehydrogenase activities, mitochondrial activity, and expression of genes coding for uncoupling proteins 2 and 3. Differences between RFI pigs in plasma leptin, cortisol, and thyroid hormone concentrations are presented and discussed. In conclusion, selection for low RFI influenced muscle properties in a way favoring muscle mass, but likely impairing meat quality.

Review: divergent selection for residual feed intake in the growing pig

This review summarizes the results from the INRA (Institut National de la Recherche Agronomique) divergent selection experiment on residual feed intake (RFI) in growing Large White pigs during nine generations of selection. It discusses the remaining challenges and perspectives for the improvement of feed efficiency in growing pigs. The impacts on growing pigs raised under standard conditions and in alternative situations such as heat stress, inflammatory challenges or lactation have been studied. After nine generations of selection, the divergent selection for RFI led to highly significant (P<0.001) line differences for RFI (-165 g/day in the low RFI (LRFI) line compared with high RFI line) and daily feed intake (-270 g/day). Low responses were observed on growth rate (-12.8 g/day, P<0.05) and body composition (+0.9 mm backfat thickness, P=0.57; -2.64% lean meat content, P<0.001) with a marked response on feed conversion ratio (-0.32 kg feed/kg gain, P<0.001). Reduced ultimate pH and increased lightness of the meat (P<0.001) were observed in LRFI pigs with minor impact on the sensory quality of the meat. These changes in meat quality were associated with changes of the muscular energy metabolism. Reduced maintenance energy requirements (-10% after five generations of selection) and activity (-21% of time standing after six generations of selection) of LRFI pigs greatly contributed to the gain in energy efficiency. However, the impact of selection for RFI on the protein metabolism of the pig remains unclear. Digestibility of energy and nutrients was not affected by selection, neither for pigs fed conventional diets nor for pigs fed high-fibre diets. A significant improvement of digestive efficiency could likely be achieved by selecting pigs on fibre diets. No convincing genetic or blood biomarker has been identified for explaining the differences in RFI, suggesting that pigs have various ways to achieve an efficient use of feed. No deleterious impact of the selection on the sow reproduction performance was observed. The resource allocation theory states that low RFI may reduce the ability to cope with stressors, via the reduction of a buffer compartment dedicated to responses to stress. None of the experiments focussed on the response of pigs to stress or challenges could confirm this theory. Understanding the relationships between RFI and responses to stress and energy demanding processes, as such immunity and lactation, remains a major challenge for a better understanding of the underlying biological mechanisms of the trait and to reconcile the experimental results with the resource allocation theory.

Identification of QTL with effects on intramuscular fat content and fatty acid composition in a Duroc x Large White cross

BACKGROUND

Improving pork quality can be done by increasing intramuscular fat (IMF) content. This trait is influenced by quantitative trait loci (QTL) sought out in different pig populations. Considering the high IMF content observed in the Duroc pig, it was appealing to determine whether favourable alleles at a major gene or QTL could be found. The detection was performed in an experimental F2 Duroc x Large White population first by segregation analysis, then by QTL mapping using additional molecular information.

RESULTS

Segregation analysis provided evidence for a major gene, with a recessive Duroc allele increasing IMF by 1.8% in Duroc homozygous pigs. However, results depended on whether data were normalised or not. After Box-Cox transformation, likelihood ratio was indeed 12 times lower and no longer significant. The QTL detection results were partly consistent with the segregation analysis. Three QTL significant at the chromosome wide level were evidenced. Two QTL, located on chromosomes 13 and 15, showed a high IMF Duroc recessive allele with an overall effect slightly lower than that expected from segregation analysis (+0.4 g/100 g muscle). The third QTL was located on chromosome 1, with a dominant Large White allele inducing high IMF content (+0.5 g/100 g muscle). Additional QTL were detected for muscular fatty acid composition.

CONCLUSION

The study presented results from two complementary approaches, a segregation analysis and a QTL detection, to seek out genes involved in the higher IMF content observed in the Duroc population. Discrepancies between both methods might be partially explained by the existence of at least two QTL with similar characteristics located on two different chromosomes for which different boars were heterozygous. The favourable and dominant allele detected in the Large White population was unexpected. Obviously, in both populations, the favourable alleles inducing high IMF content were not fixed and improving IMF by fixing favourable alleles using markers can then be applied both in Duroc and LW populations. With QTL affecting fatty acid composition, combining an increase of IMF content enhancing monounsaturated fatty acid percentage would be of great interest.

Effects of quantitative trait loci on chromosomes 1, 2, 4, and 7 on growth, carcass, and meat quality traits in backcross Meishan x Large White pigs

The aim of this work was to estimate whether genetic dissection of QTL on chromosomes 1, 2, 4, and 7, detected in an F2 Meishan x Large White population, can be achieved with a recombinant back-cross progeny test approach. For this purpose, a first generation of backcross (BC1) was produced by using frozen semen of F1 Large White x Meishan boars with Large White females. Four BC1 boars were selected because of their heterozygosity for at least 1 of the 4 regions. The BC1 boars were crossed with Large White sows, and the resulting BC2 offspring were measured for several growth and body composition traits. Contrary to the F2 animals, BC2 animals were also measured for meat quality traits in adductor, gluteus superficialis (GS), longissimus dorsi, and biceps femoris (BF) muscles. Each BC1 boar was tested for a total of 39 traits and for the 4 regions with statistical interval mapping analyses. The QTL effects obtained in BC1 families showed some differences compared with those described in F1 families. However, we confirmed QTL effects for growth in the SW1301-SW2512 markers interval on chromosome 1 and also for body composition in the SW1828-SW2512 markers interval on chromosome 1, in the SW2443-SWR783 markers interval on chromosome 2, and in the SW1369-SW632 markers interval on chromosome 7. In addition, we detected new QTL for growth traits on chromosome 2 and for meat quality traits on chromosomes 1 and 2. Growth of animals from weaning to the end of the test was influenced by the IGF2 gene region on chromosome 2. Concerning meat quality, ultimate pH of adductor, longissimus dorsi, and BF were affected by the interval delimited by UMNP3000 and SW2512 markers on chromosome 1, and a* of GS, L* of BF, and water-holding capacity of GS were affected by QTL located between marker loci SW2443 and SWR783 on chromosome 2. Recombinant progeny testing appeared to be a suitable strategy for the genetic dissection of the QTL investigated.

Correlated responses in sow appetite, residual feed intake, body composition, and reproduction after divergent selection for residual feed intake in the growing pig

Residual feed intake (RFI) has been explored as an alternative selection criterion to feed conversion ratio to capture the fraction of feed intake not explained by expected production and maintenance requirements. Selection experiments have found that low RFI in the growing pig is genetically correlated with reduced fatness and feed intake. Selection for feed conversion ratio also reduces sow appetite and fatness, which, together with increased prolificacy, has been seen as a hindrance for sow lifetime performance. The aims of our study were to derive equations for sow RFI during lactation (SRFI) and to evaluate the effect of selection for RFI during growth on sow traits during lactation. Data were obtained on 2 divergent lines selected for 7 generations for low and high RFI during growth in purebred Large Whites. The RFI was measured on candidates for selection (1,065 pigs), and sow performance data were available for 480 sows having from 1 to 3 parities (1,071 parities). Traits measured were sow daily feed intake (SDFI); sow BW and body composition before farrowing and at weaning (28.4 ± 1.7d); number of piglets born total, born alive, and surviving at weaning; and litter weight, average piglet BW, and within-litter SD of piglet BW at birth, 21 d of age (when creep feeding was available), and weaning. Sow RFI was defined as the difference between observed SDFI and SDFI predicted for sow maintenance and production. Daily production requirements were quantified by litter size and daily litter BW gain as well as daily changes in sow body reserves. The SRFI represented 24% of the phenotypic variability of SDFI. Heritability estimates for RFI and SRFI were both 0.14. The genetic correlation between RFI and SRFI was 0.29 ± 0.23. Genetic correlations of RFI with sow traits were low to moderate, consistent with responses to selection; selection for low RFI during growth reduced SDFI and increased number of piglets and litter growth, but also increased mobilization of body reserves. No effect on rebreeding performance was found. Metabolic changes previously observed during growth in response to selection might explain part of the better efficiency of the low-RFI sows, decreasing basal metabolism and favoring rapid allocation of resources to lactation. We propose to consider SRFI as an alternative to SDFI to select for efficient sows with reduced input demands during lactation.

Effect of chronic and acute heat challenges on fecal microbiota composition, production, and thermoregulation traits in growing pigs1,2

The present study aimed at investigating the impact of heat challenges on gut microbiota composition in growing pigs and its relationship with pigs’ performance and thermoregulation responses. From a total of 10 F1 sire families, 558 and 564 backcross Large White × Créole pigs were raised and phenotyped from 11 to 23 wk of age in temperate (TEMP) and in tropical (TROP) climates, respectively. In TEMP, all pigs were subjected to an acute heat challenge (3 wk at 29 °C) from 23 to 26 wk of age. Feces samples were collected at 23 wk of age both in TEMP and TROP climate (TEMP23 and TROP23 samples, respectively) and at 26 wk of age in TEMP climate (TEMP26 samples) for 16S rRNA analyses of fecal microbiota composition. The fecal microbiota composition significantly differed between the 3 environments. Using a generalized linear model on microbiota composition, 182 operational taxonomic units (OTU) and 2 pathways were differentially abundant between TEMP23 and TEMP26, and 1,296 OTU and 20 pathways between TEMP23 and TROP23. Using fecal samples collected at 23 wk of age, pigs raised under the 2 climates were discriminated with 36 OTU using a sparse partial least square discriminant analysis that had a mean classification error-rate of 1.7%. In contrast, pigs in TEMP before the acute heat challenge could be discriminated from the pigs in TEMP after the heat challenge with 32 OTU and 9.3% error rate. The microbiota can be used as biomarker of heat stress exposition. Microbiota composition revealed that pigs were separated into 2 enterotypes. The enterotypes were represented in both climates. Whatever the climate, animals belonging to the Turicibacter–Sarcina–Clostridium sensu stricto dominated enterotype were 3.3 kg heavier (P < 0.05) at 11 wk of age than those belonging to the Lactobacillus-dominated enterotype. This latter enterotype was related to a 0.3 °C lower skin temperature (P < 0.05) at 23 wk of age. Following the acute heat challenge in TEMP, this enterotype had a less-stable rectal temperature (0.34 vs. 0.25 °C variation between weeks 23 and 24, P < 0.05) without affecting growth performance (P > 0.05). Instability of the enterotypes was observed in 34% of the pigs, switching from an enterotype to another between 23 and 26 wk of age after heat stress. Despite a lower microbial diversity, the Turicibacter–Sarcina–Clostridium sensu stricto dominated enterotype was better adapted to heat stress conditions with lower thermoregulation variations.

Extensive expression differences along porcine small intestine evidenced by transcriptome sequencing

The aim of this study was to analyse gene expression along the small intestine (duodenum, jejunum, ileum) and in the ileal Peyer's patches in four young pigs with no clinical signs of disease by transcriptome sequencing. Multidimensional scaling evidenced that samples clustered by tissue type rather than by individual, thus prefiguring a relevant scenario to draw tissue-specific gene expression profiles. Accordingly, 1,349 genes were found differentially expressed between duodenum and jejunum, and up to 3,455 genes between duodenum and ileum. Additionally, a considerable number of differentially expressed genes were found by comparing duodenum (7,027 genes), jejunum (6,122 genes), and ileum (6,991 genes) with ileal Peyer's patches tissue. Functional analyses revealed that most of the significant differentially expressed genes along small intestinal tissues were involved in the regulation of general biological processes such as cell development, signalling, growth and proliferation, death and survival or cell function and maintenance. These results suggest that the intrinsic large turnover of intestinal tissues would have local specificities at duodenum, ileum and jejunum. In addition, in concordance with their biological function, enteric innate immune pathways were overrepresented in ileal Peyer's patches. The reported data provide an expression map of the cell pathway variation in the different small intestinal tissues. Furthermore, expression levels measured in healthy individuals could help to understand changes in gene expression that occur in dysbiosis or pathological states.

Towards the quantitative characterisation of piglets' robustness to weaning: a modelling approach

Weaning is a critical transition phase in swine production in which piglets must cope with different stressors that may affect their health. During this period, the prophylactic use of antibiotics is still frequent to limit piglet morbidity, which raises both economic and public health concerns such as the appearance of antimicrobial-resistant microbes. With the interest of developing tools for assisting health and management decisions around weaning, it is key to provide robustness indexes that inform on the animals' capacity to endure the challenges associated with weaning. This work aimed at developing a modelling approach for facilitating the quantification of piglet resilience to weaning. A total of 325 Large White pigs weaned at 28 days of age were monitored and further housed and fed conventionally during the post-weaning period without antibiotic administration. Body weight and diarrhoea scores were recorded before and after weaning, and blood was sampled at weaning and 1 week later for collecting haematological data. A dynamic model was constructed based on the Gompertz-Makeham law to describe live weight trajectories during the first 75 days after weaning, following the rationale that the animal response is partitioned in two time windows (a perturbation and a recovery window). Model calibration was performed for each animal. Our results show that the transition time between the two time windows, as well as the weight trajectories are characteristic for each individual. The model captured the weight dynamics of animals at different degrees of perturbation, with an average coefficient of determination of 0.99, and a concordance correlation coefficient of 0.99. The utility of the model is that it provides biologically meaningful parameters that inform on the amplitude and length of perturbation, and the rate of animal recovery. Our rationale is that the dynamics of weight inform on the capability of the animal to cope with the weaning disturbance. Indeed, there were significant correlations between model parameters and individual diarrhoea scores and haematological traits. Overall, the parameters of our model can be useful for constructing weaning robustness indexes by using exclusively the growth curves. We foresee that this modelling approach will provide a step forward in the quantitative characterisation of robustness.

Muscle transcriptomic investigation of late fetal development identifies candidate genes for piglet maturity

Background

In pigs, the perinatal period is the most critical time for survival. Piglet maturation, which occurs at the end of gestation, leads to a state of full development after birth. Therefore, maturity is an important determinant of early survival. Skeletal muscle plays a key role in adaptation to extra-uterine life, e.g. glycogen storage and thermoregulation. In this study, we performed microarray analysis to identify the genes and biological processes involved in piglet muscle maturity. Progeny from two breeds with extreme muscle maturity phenotypes were analyzed at two time points during gestation (gestational days 90 and 110). The Large White (LW) breed is a selected breed with an increased rate of mortality at birth, whereas the Meishan (MS) breed produces piglets with extremely low mortality at birth. The impact of the parental genome was analyzed with reciprocal crossed fetuses.

Results

Microarray analysis identified 12,326 differentially expressed probes for gestational age and genotype. Such a high number reflects an important transcriptomic change that occurs between 90 and 110 days of gestation. 2,000 probes, corresponding to 1,120 unique annotated genes, involved more particularly in the maturation process were further studied. Functional enrichment and graph inference studies underlined genes involved in muscular development around 90 days of gestation, and genes involved in metabolic functions, such as gluconeogenesis, around 110 days of gestation. Moreover, a difference in the expression of key genes, e.g. PCK2, LDHA or PGK1, was detected between MS and LW just before birth. Reciprocal crossing analysis resulted in the identification of 472 genes with an expression preferentially regulated by one parental genome. Most of these genes (366) were regulated by the paternal genome. Among these paternally regulated genes, some known imprinted genes, such as MAGEL2 or IGF2, were identified and could have a key role in the maturation process.

Conclusion

These results reveal the biological mechanisms that regulate muscle maturity in piglets. Maturity is also under the conflicting regulation of the parental genomes. Crucial genes, which could explain the biological differences in maturity observed between LW and MS breeds, were identified. These genes could be excellent candidates for a key role in the maturity.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-797) contains supplementary material, which is available to authorized users.

Meiotic recombination analyses of individual chromosomes in male domestic pigs (Sus scrofa domestica)

For the first time in the domestic pig, meiotic recombination along the 18 porcine autosomes was directly studied by immunolocalization of MLH1 protein. In total, 7,848 synaptonemal complexes from 436 spermatocytes were analyzed, and 13,969 recombination sites were mapped. Individual chromosomes for 113 of the 436 cells (representing 2,034 synaptonemal complexes) were identified by immunostaining and fluorescence in situ hybridization (FISH). The average total length of autosomal synaptonemal complexes per cell was 190.3 µm, with 32.0 recombination sites (crossovers), on average, per cell. The number of crossovers and the lengths of the autosomal synaptonemal complexes showed significant intra- (i.e. between cells) and inter-individual variations. The distributions of recombination sites within each chromosomal category were similar: crossovers in metacentric and submetacentric chromosomes were concentrated in the telomeric regions of the p- and q-arms, whereas two hotspots were located near the centromere and in the telomeric region of acrocentrics. Lack of MLH1 foci was mainly observed in the smaller chromosomes, particularly chromosome 18 (SSC18) and the sex chromosomes. All autosomes displayed positive interference, with a large variability between the chromosomes.

Time course of the response to ACTH in pig: biological and transcriptomic study

Background

HPA axis plays a major role in physiological homeostasis. It is also involved in stress and adaptive response to the environment. In farm animals in general and specifically in pigs, breeding strategies have highly favored production traits such as lean growth rate, feed efficiency and prolificacy at the cost of robustness. On the hypothesis that the HPA axis could contribute to the trade-off between robustness and production traits, we have designed this experiment to explore individual variation in the biological response to the main stress hormone, cortisol, in pigs. We used ACTH injections to trigger production of cortisol in 120 juvenile Large White (LW) pigs from 28 litters and the kinetics of the response was measured with biological variables and whole blood gene expression at 4 time points. A multilevel statistical analysis was used to take into account the longitudinal aspect of the data.

Results

Cortisol level reached its peak 1 h after ACTH injection. White blood cell composition was modified with a decrease of lymphocytes and monocytes and an increase of granulocytes (FDR<0.05). Basal level of cortisol was correlated with birth and weaning weights. Microarray analysis identified 65 unique genes of which expression responded to the injection of ACTH (adjusted P<0.05). These genes were classified into 4 clusters with distinctive kinetics in response to ACTH injection. The first cluster identified genes strongly correlated to cortisol and previously reported as being regulated by glucocorticoids. In particular, DDIT4, DUSP1, FKBP5, IL7R, NFKBIA, PER1, RGS2 and RHOB were shown to be connected to each other by the glucocorticoid receptor NR3C1. Most of the differentially expressed genes that encode transcription factors have not been described yet as being important in transcription networks involved in stress response. Their co-expression may mean co-regulation and they could thus provide new patterns of biomarkers of the individual sensitivity to cortisol.

Conclusions

We identified 65 genes as biological markers of HPA axis activation at the gene expression level. These genes might be candidates for a better understanding of the molecular mechanisms of the stress response.

Electronic supplementary material

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

Exclusion of the swine leukocyte antigens as candidate region and reduction of the position interval for the Sus scrofa chromosome 7 QTL affecting growth and fatness1

Pig chromosome 7 (SSC 7) has been shown to be rich in QTL affecting performance and quality traits. Most studies mapped the QTL close to the swine leukocyte antigens (SLA), which has a large effect on adaptability and natural selection. Previous comparative mapping studies suggested that the 15-cM region limited by markers LRA1 (mapped at 55 cM) and S0102 (mapped at 70 cM) contains hundreds of genes. To decrease the number of candidate genes, we improved the mapping resolution with a genetic chromosome dissection through a backcross recombinant progeny test program between Meishan (MS) and European (EU; i.e., Large White or Landrace) breeds. Three first-generation backcross--(EU x MS) x EU--and two second-generation backcross--([EU x MS] x EU) x EU--sires carrying a recombination in the QTL mapping interval were progeny-tested (i.e., measured for a total of 44 growth, fatness, carcass and meat quality traits). Progeny family size varied from 29 to 119 pigs. Animals were genotyped for markers covering the region of interest. Progeny-test results allowed the QTL interval to be decreased from 15 to 20 cM down to 10 cM, and even less than 6 cM if we assumed that the EU pigs used in this study share only one QTL allele. Except for a putative QTL affecting some carcass composition traits, the SLA is excluded as a candidate region, suggesting that it might be possible to apply a marker-assisted selection strategy for this QTL, while controlling SLA allele diversity. The strong QTL effects remaining in animals with only 12.5% (issued from first-generation backcross boars) and 6.25% (issued from second-generation back-cross boars) Meishan genetic background shows that epistatic interactions are likely to be limited. Finally, the QTL does not have strong effects on meat quality traits.

Integrated Analysis of Proteomic and Transcriptomic Data Highlights Late Fetal Muscle Maturation Process

In pigs, the perinatal period is the most critical time for survival. Piglet maturation, which occurs at the end of gestation, is an important determinant of early survival. Skeletal muscle plays a key role in adaptation to extra-uterine life, e.g. motor function and thermoregulation. Progeny from two breeds with extreme neonatal mortality rates were analyzed at 90 and 110 days of gestation (dg). The Large White breed is a highly selected breed for lean growth and exhibits a high rate of neonatal mortality, whereas the Meishan breed is fatter and more robust and has a low neonatal mortality. Our aim was to identify molecular signatures underlying late fetal longissimus muscle development. First, integrated analysis was used to explore relationships between co-expression network models built from a proteomic data set (bi-dimensional electrophoresis) and biological phenotypes. Second, correlations with a transcriptomic data set (microarrays) were investigated to combine different layers of expression with a focus on transcriptional regulation. Muscle glycogen content and myosin heavy chain polymorphism were good descriptors of muscle maturity and were used for further data integration analysis. Using 89 identified unique proteins, network inference, correlation with biological phenotypes and functional enrichment revealed that mitochondrial oxidative metabolism was a key determinant of neonatal muscle maturity. Some proteins, including ATP5A1 and CKMT2, were important nodes in the network related to muscle metabolism. Transcriptomic data suggest that overexpression of mitochondrial PCK2 was involved in the greater glycogen content of Meishan fetuses at 110 dg. GPD1, an enzyme involved in the mitochondrial oxidation of cytosolic NADH, was overexpressed in Meishan. Thirty-one proteins exhibited a positive correlation between mRNA and protein levels in both extreme fetal genotypes, suggesting transcriptional regulation. Gene ontology enrichment and Ingenuity analyses identified PPARGC1A and ESR1 as possible transcriptional factors positively involved in late fetal muscle maturation.

Genotype by environment interactions for performance and thermoregulation responses in growing pigs1,2

Heat stress affects pig health, welfare, and production, and thus the economic viability of the pig sector in many countries. Breeding for heat tolerance is a complex issue, increasingly important due to climate change and the development of pig production in tropical areas. Characterizing genetic determinism of heat tolerance would help building selection schemes dedicated to high performance in tropical areas. The main objective of our study was to estimate the genetic parameters for production and thermoregulation traits in two highly related growing pig populations reared in temperate (TEMP) or tropical humid (TROP) environment. Pigs came from a backcross population between Large White (LW, heat sensitive) and Creole (CR, heat tolerant) pigs. Phenotypic data were obtained on a total of 1,297 pigs using the same procedures in both environments, for body weight (BW, at weeks 11 and 23), daily feed intake (ADFI), backfat thickness (BFT, at weeks 19 and 23), cutaneous temperature (CT, at weeks 19 and 23), and rectal temperature (RT, at weeks 19, 21, and 23). Feed conversion ratio (FCR) and residual feed intake (RFI) were computed for the whole test period (11 to 23 wk). Criteria comparing the fits to the data revealed genotype × environment (G × E) interactions for most traits but not for FCR. The variance components were obtained using two different methods, a restricted maximum likelihood method and a Bayesian Markov chain Monte Carlo method, considering that traits are either similar or different in each environment. Regardless of the method, heritability estimates for production traits were moderate to high, except for FCR (lower than 0.18). Heritability estimates for RT were low to moderate, ranging from 0.04 to 0.34. The genetic correlations of each trait between environments generally differed from 1, except for FCR and ADG. For most thermoregulation traits, they also did not differ significantly from zero, suggesting that the main genetic bases of heat tolerance may vary in different environment. Within environments, the unfavorable genetic correlations between production traits and RT suggest an antagonism between the ability to maintain inner temperature and the ability to increase ADFI and ADG. However, greater RT were also associated to leaner pigs and better feed efficiency. Nevertheless, due to large inaccuracies of these estimations, larger cohorts would be needed to decide about the best breeding schemes to choose for tropical pig production.

1HNMR-Based metabolomic profiling method to develop plasma biomarkers for sensitivity to chronic heat stress in growing pigs

The negative impact of heat stress (HS) on the production performances in pig faming is of particular concern. Novel diagnostic methods are needed to predict the robustness of pigs to HS. Our study aimed to assess the reliability of blood metabolome to predict the sensitivity to chronic HS of 10 F1 (Large White × Creole) sire families (SF) reared in temperate (TEMP) and in tropical (TROP) regions (n = 56±5 offsprings/region/SF). Live body weight (BW) and rectal temperature (RT) were recorded at 23 weeks of age. Average daily feed intake (AFDI) and average daily gain were calculated from weeks 11 to 23 of age, together with feed conversion ratio. Plasma blood metabolome profiles were obtained by Nuclear Magnetic Resonance spectroscopy (^1HNMR) from blood samples collected at week 23 in TEMP. The sensitivity to hot climatic conditions of each SF was estimated by computing a composite index of sensitivity (I_sens) derived from a linear combination of t statistics applied to familial BW, ADFI and RT in TEMP and TROP climates. A model of prediction of sensitivity was established with sparse Partial Least Square Discriminant Analysis (sPLS-DA) between the two most robust SF (n = 102) and the two most sensitive ones (n = 121) using individual metabolomic profiles measured in TEMP. The sPLS-DA selected 29 buckets that enabled 78% of prediction accuracy by cross-validation. On the basis of this training, we predicted the proportion of sensitive pigs within the 6 remaining families (n = 337). This proportion was defined as the predicted membership of families to the sensitive category. The positive correlation between this proportion and I_sens (r = 0.97, P < 0.01) suggests that plasma metabolome can be used to predict the sensitivity of pigs to hot climate.

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 (h²) and genetic correlations (r_g) 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 h² estimates of 0.19 ± 0.08 (Shannon) and 0.12 ± 0.06 (Simpson). Forty‐eight genera had a significant h² (>0.125). The r_g of the α‐diversities indexes with production traits were negative. Some r_g 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.

Acute social stress-induced immunomodulation in pigs high and low responders to ACTH

Pig husbandry is known as an intensive breeding system, piglets being submitted to multiple stressful events such as early weaning, successive mixing, crowding and shipping. These stressors are thought to impair immune defences and might contribute, at least partly, to the prophylactic use of antibiotics. Robustness was recently defined as the ability of an individual to express a high-production potential in a wide variety of environmental conditions. Increasing robustness thus appears as a valuable option to improve resilience to stressors and could be obtained by selecting piglets upon their adrenocortical activity. In this study, we aimed at depicting the consequences of an acute social stress on the immune capacity of piglets genetically selected upon divergent hypothalamic-pituitary-adrenocortical (HPA) axis activity. For this purpose, we monitored neuroendocrine and immune parameters, in high- (HPA^hi) and low- (HPA^lo) responders to ACTH, just before and immediately after a one-hour mixing with unfamiliar conspecifics. As expected, stressed piglets displayed higher levels of circulating cortisol and norepinephrine. Blood cell count analysis combined to flow cytometry revealed a stress-induced leukocyte mobilization in the bloodstream with a specific recruitment of CD8α⁺ lymphocytes. Besides, one-hour mixing decreased LPS-induced IL-8 and TNFα secretions in whole-blood assays (WBA) and reduced mononuclear cell phagocytosis. Altogether, our data demonstrate that acute social stress alters immune competence of piglets from both groups, and bring new insights in favour of good farming practices. While for most parameters high- and low-responders to ACTH behaved similarly, HPA^hi piglets displayed higher number of CD4⁺ CD8α^- T cells, as well as increased cytokine production in WBA (LPS-induced TNFα and PIL-induced IL-8), which could confer them increased resistance to pathogens. Finally, a principal component analysis including all parameters highlighted that overall stress effects were less pronounced on piglets with a strong HPA axis. Thus, selection upon adrenocortical axis activity seems to reduce the magnitude of response to stress and appears as a good tool to increase piglet robustness.

Metabolic and histochemical characteristics of fat and muscle tissues in homozygous or heterozygous pigs for the body composition QTL located on chromosome 7

Quantitative trait loci (QTL) influencing many traits including backfat thickness and carcass composition have been detected on porcine chromosome 7 (SSC7) in an F2 cross between Large White (LW) and Meishan (MS) pigs. However, the genes and controlled pathways underlying the QTL effects on body phenotype remain unknown. This study aimed at investigating the tissue characteristics at metabolic and cellular levels in pigs that were either homozygous or heterozygous for a body composition SSC7 QTL. A backcross pig (BC3) was first progeny tested to confirm its heterozygoty for the SSC7 QTL; results on all offspring (n = 80) confirmed the QTL effects on body fatness. This boar was then mated with three sows known to be heterozygous for this QTL. In the subset of pigs per genotype, we found that heterozygous LW(QTL7)/MS(QTL7) pigs had smaller adipocytes in backfat, together with a lower basal rate of glucose incorporation into lipids and lower activities of selected lipogenic enzymes in backfat isolated cells, compared with homozygous LW(QTL7)/LW(QTL7) pigs. A higher number of adipocytes was also estimated in backfat of LW(QTL7)/MS(QTL7) animals compared with LW(QTL7)/LW(QTL7) pigs. The SSC7 QTL did not influence oxidative and glycolytic metabolisms of longissimus and trapezius muscles, as estimated by the activities of specific energy metabolism enzymes, or the myofiber type properties. Altogether, this study provides new evidence for an altered adipocyte cellularity in backfat of pigs carrying at least one MS allele for the SSC7 QTL. Some candidate genes known for their functions on adipocyte growth and differentiation are suggested.

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.