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.

Sex and fetal genome influence gene expression in pig endometrium at the end of gestation.

BACKGROUND: A fine balance of feto-maternal resource allocation is required to support pregnancy, which depends on interactions between maternal and fetal genetic potential, maternal nutrition and environment, endometrial and placental functions. In particular, some imprinted genes have a role in regulating maternal-fetal nutrient exchange, but few have been documented in the endometrium. The aim of this study is to describe the expression of 42 genes, with parental expression, in the endometrium comparing two extreme breeds: Large White (LW); Meishan (MS) with contrasting neonatal mortality and maturity at two days of gestation (D90-D110). We investigated their potential contribution to fetal maturation exploring genes-fetal phenotypes relationships. Last, we hypothesized that the fetal genome and sex influence their endometrial expression. For this purpose, pure and reciprocally crossbred fetuses were produced using LW and MS breeds. Thus, in the same uterus, endometrial samples were associated with its purebred or crossbred fetuses. RESULTS: Among the 22 differentially expressed genes (DEGs), 14 DEGs were differentially regulated between the two days of gestation. More gestational changes were described in LW (11 DEGs) than in MS (2 DEGs). Nine DEGs were differentially regulated between the two extreme breeds, highlighting differences in the regulation of endometrial angiogenesis, nutrient transport and energy metabolism. We identified DEGs that showed high correlations with indicators of fetal maturation, such as ponderal index at D90 and fetal blood fructose level and placental weight at D110. We pointed out for the first time the influence of fetal sex and genome on endometrial expression at D90, highlighting AMPD3, CITED1 and H19 genes. We demonstrated that fetal sex affects the expression of five imprinted genes in LW endometrium. Fetal genome influenced the expression of four genes in LW endometrium but not in MS endometrium. Interestingly, both fetal sex and fetal genome interact to influence endometrial gene expression. CONCLUSIONS: These data provide evidence for some sexual dimorphism in the pregnant endometrium and for the contribution of the fetal genome to feto-maternal interactions at the end of gestation. They suggest that the paternal genome may contribute significantly to piglet survival, especially in crossbreeding production systems.

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.

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.

Influence of genetics and the pre-vaccination blood transcriptome on the variability of antibody levels after vaccination against Mycoplasma hyopneumoniae in pigs.

BACKGROUND: The impact of individual genetic and genomic variations on immune responses is an emerging lever investigated in vaccination strategies. In our study, we used genetic and pre-vaccination blood transcriptomic data to study vaccine effectiveness in pigs. RESULTS: A cohort of 182 Large White pigs was vaccinated against Mycoplasma hyopneumoniae (M. hyo) at weaning (28 days of age), with a booster 21 days later. Vaccine response was assessed by measuring seric M. hyo antibodies (Ab) at 0 (vaccination day), 21 (booster day), 28, 35, and 118 days post-vaccination (dpv). Inter-individual variability of M. hyo Ab levels was observed at all time points and the corresponding heritabilities ranged from 0.46 to 0.57. Ab persistence was higher in females than in males. Genome-wide association studies with a 658 K SNP panel revealed two genomic regions associated with variations of M. hyo Ab levels at 21 dpv at positions where immunity-related genes have been mapped, DAB2IP on chromosome 1, and ASAP1, CYRIB and GSDMC on chromosome 4. We studied covariations of Ab responses with the pre-vaccination blood transcriptome obtained by RNA-Seq for a subset of 82 pigs. Weighted gene correlation network and differential expression analyses between pigs that differed in Ab responses highlighted biological functions that were enriched in heme biosynthesis and platelet activation for low response at 21 dpv, innate antiviral immunity and dendritic cells for high response at 28 and 35 dpv, and cell adhesion and extracellular matrix for high response at 118 dpv. Sparse partial least squares discriminant analysis identified 101 genes that efficiently predicted divergent responders at all time points. We found weak negative correlations of M. hyo Ab levels with body weight traits, which revealed a trade-off that needs to be further explored. CONCLUSIONS: We confirmed the influence of the host genetics on vaccine effectiveness to M. hyo and provided evidence that the pre-vaccination blood transcriptome co-varies with the Ab response. Our results highlight that both genetic markers and blood biomarkers could be used as potential predictors of vaccine response levels and more studies are required to assess whether they can be exploited in breeding programs.

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.

Adipokines expression profiles in both plasma and peri renal adipose tissue in Large White and Meishan sows: A possible involvement in the fattening and the onset of puberty.

In pig, backfat deposition is strongly related to the growth and reproductive performance. However, the molecular regulatory mechanisms of adipose tissue are not clearly understood. Adipose tissue is now recognized as an important endocrine organ that secretes a variety of factors including adipokines. However, the regulation of expression pattern of these adipokines in both plasma and visceral white adipose tissue (WAT) in lean and fat pig is unclear. In the present study, we used two representative porcine breeds (Large White, LW; Meishan, MS) with contrasting backfat thickness and sexual maturity age. Using specific ELISA assays, we determined the plasma profile of eight adipokines, leptin, adiponectin, visfatin, apelin, chemerin, resistin, omentin and vaspin in LW and MS sows. By RT-qPCR and western-blot we also investigated the mRNA and protein levels of these adipokines and their cognate receptors (LEPR, ADIPOR1, ADIPOR2, CMKLR1, CCRL2, GPR1, APLNR, TLR4, ROR1, CAP1 and HSPA5) in the peri renal WAT, respectively. At both plasma and peri renal WAT level, we found that the amounts of leptin, chemerin, resistin and vaspin were higher whereas those of adiponectin and omentin were lower in MS than LW sows. Plasma and adipose tissue visfatin and apelin levels were not different between the two breeds. Moreover, we noted that the variations of peri renal WAT adipokines observed between MS and LW were similar at the protein and mRNA level except for chemerin and apelin suggesting post-transcriptional modifications for these two adipokines. Finally, among the eight adipokines studied, we showed that only the plasma concentrations of leptin and chemerin were positively and those of adiponectin, negatively associated with the thickness of fat and opposite correlation was found for the onset of puberty in both LW and MS animals. Taken together, these results support a potential involvement of adipokines in WAT regulation and its link with the onset of the puberty in sows.

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.

Effects of divergent selection upon adrenocortical activity on immune traits in pig.

BACKGROUND: The sustainability of farming and animal welfare requires the reconsideration of current selection schemes. In particular, implementation of new selection criteria related to animal health and welfare should help to produce more robust animals and to reduce anti-microbial use. The hypothalamo-pituitary-adrenocortical (HPA) axis plays a major role in metabolic regulation and adaptation processes and its activity is strongly influenced by genetic factors. A positive association between HPA axis activity and robustness was recently described. To explore whether selecting pigs upon HPA axis activity could increase their robustness, a divergent selection experiment was carried out in the Large White pig breed. This allowed the generation of low (HPAlo) and high (HPAhi) responders to adrenocorticotropic hormone administration. RESULTS: In this study, we compared 23 hematologic and immune parameters of 6-week-old, HPAlo and HPAhi piglets and analysed their response to a low dose of lipopolysaccharide (LPS) two weeks later. At six weeks of age, HPAhi piglets displayed greater red blood cell and leucocyte number including CD8α+ γδ cells, cytotoxic T lymphocytes, naive T helper (Th) cells and B lymphocytes as compared to HPAlo individuals. The ability of blood cells to secrete TNFα in response to LPS ex vivo was higher for HPAhi pigs. At week eight, the inflammatory response to the LPS in vivo challenge was poorly affected by the HPA axis activity. CONCLUSIONS: Divergent selection upon HPA axis activity modulated hematologic and immune parameters in 6-week-old pigs, which may confer an advantage to HPAhi pigs at weaning. However, HPAlo and HPAhi piglets did not exhibit major differences in the parameters analysed two weeks later, i. e. in 8-week-old pigs. In conclusion, chronic exposure to high cortisol levels in HPAhi pigs does not negatively impact immunity.

Time course study of the response to LPS targeting the pig immune gene networks.

BACKGROUND: Stress is a generic term used to describe non-specific responses of the body to all kinds of challenges. A very large variability in the response can be observed across individuals, depending on numerous conditioning factors like genetics, early influences and life history. As a result, there is a wide range of individual vulnerability and resilience to stress, also called robustness. The importance of robustness-related traits in breeding strategies is increasing progressively towards the production of animals with a high level of production under a wide range of climatic conditions and management systems, together with a lower environmental impact and a high level of animal welfare. The present study aims at describing blood transcriptomic, hormonal, and metabolic responses of pigs to a systemic challenge using lipopolysaccharide (LPS). The objective is to analyze the individual variation of the biological responses in relation to the activity of the HPA axis measured by the levels of plasma cortisol after LPS and ACTH in 120 juvenile Large White (LW) pigs. 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 reaches its peak 4 h after LPS injection. The characteristic changes of white blood cell count to LPS were observed, with a decrease of total count, maximal at t=+4 h, and the mirror changes in the respective proportions of lymphocytes and granulocytes. The lymphocytes / granulocytes ratio was maximal at t=+1 h. An integrative statistical approach was used and provided a set of candidate genes for kinetic studies and ongoing complementary studies focused on the LPS-stimulated inflammatory response. CONCLUSIONS: The present study demonstrates the specific biomarkers indicative of an inflammation in swine. Furthermore, these stress responses persist for prolonged periods of time and at significant expression levels, making them good candidate markers for evaluating the efficacy of anti-inflammatory drugs.

Comparing the intestinal transcriptome of Meishan and Large White piglets during late fetal development reveals genes involved in glucose and lipid metabolism and immunity as valuable clues of intestinal maturity.

BACKGROUND: Maturity of intestinal functions is critical for neonatal health and survival, but comprehensive description of mechanisms underlying intestinal maturation that occur during late gestation still remain poorly characterized. The aim of this study was to investigate biological processes specifically involved in intestinal maturation by comparing fetal jejunal transcriptomes of two representative porcine breeds (Large White, LW; Meishan, MS) with contrasting neonatal vitality and maturity, at two key time points during late gestation (gestational days 90 and 110). MS and LW sows inseminated with mixed semen (from breed LW and MS) gave birth to both purebred and crossbred fetuses. We hypothesized that part of the differences in neonatal maturity between the two breeds results from distinct developmental profiles of the fetal intestine during late gestation. Reciprocal crossed fetuses were used to analyze the effect of parental genome. Transcriptomic data and 23 phenotypic variables known to be associated with maturity trait were integrated using multivariate analysis with expectation of identifying relevant genes-phenotypic variable relationships involved in intestinal maturation. RESULTS: A moderate maternal genotype effect, but no paternal genotype effect, was observed on offspring intestinal maturation. Four hundred and four differentially expressed probes, corresponding to 274 differentially expressed genes (DEGs), more specifically involved in the maturation process were further studied. In day 110-MS fetuses, Ingenuity® functional enrichment analysis revealed that 46% of DEGs were involved in glucose and lipid metabolism, cell proliferation, vasculogenesis and hormone synthesis compared to day 90-MS fetuses. Expression of genes involved in immune pathways including phagocytosis, inflammation and defense processes was changed in day 110-LW compared to day 90-LW fetuses (corresponding to 13% of DEGs). The transcriptional regulator PPARGC1A was predicted to be an important regulator of differentially expressed genes in MS. Fetal blood fructose level, intestinal lactase activity and villous height were the best predicted phenotypic variables with probes mostly involved in lipid metabolism, carbohydrate metabolism and cellular movement biological pathways. CONCLUSIONS: Collectively, our findings indicate that the neonatal maturity of pig intestine may rely on functional development of glucose and lipid metabolisms, immune phagocyte differentiation and inflammatory pathways. This process may partially be governed by PPARGC1A.

Multiple-trait structured antedependence model to study the relationship between litter size and birth weight in pigs and rabbits.

BACKGROUND: Some genetic studies need to take into account correlations between traits that are repeatedly measured over time. Multiple-trait random regression models are commonly used to analyze repeated traits but suffer from several major drawbacks. In the present study, we developed a multiple-trait extension of the structured antedependence model (SAD) to overcome this issue and validated its usefulness by modeling the association between litter size (LS) and average birth weight (ABW) over parities in pigs and rabbits. METHODS: The single-trait SAD model assumes that a random effect at time [Formula: see text] can be explained by the previous values of the random effect (i.e. at previous times). The proposed multiple-trait extension of the SAD model consists in adding a cross-antedependence parameter to the single-trait SAD model. This model can be easily fitted using ASReml and the OWN Fortran program that we have developed. In comparison with the random regression model, we used our multiple-trait SAD model to analyze the LS and ABW of 4345 litters from 1817 Large White sows and 8706 litters from 2286 L-1777 does over a maximum of five successive parities. RESULTS: For both species, the multiple-trait SAD fitted the data better than the random regression model. The difference between AIC of the two models (AIC_random regression-AIC_SAD) were equal to 7 and 227 for pigs and rabbits, respectively. A similar pattern of heritability and correlation estimates was obtained for both species. Heritabilities were lower for LS (ranging from 0.09 to 0.29) than for ABW (ranging from 0.23 to 0.39). The general trend was a decrease of the genetic correlation for a given trait between more distant parities. Estimates of genetic correlations between LS and ABW were negative and ranged from -0.03 to -0.52 across parities. No correlation was observed between the permanent environmental effects, except between the permanent environmental effects of LS and ABW of the same parity, for which the estimate of the correlation was strongly negative (ranging from -0.57 to -0.67). CONCLUSIONS: We demonstrated that application of our multiple-trait SAD model is feasible for studying several traits with repeated measurements and showed that it provided a better fit to the data than the random regression model.

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 (F D R<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.

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.

Detection of DNA methylation signatures through the lens of genomic imprinting.

Genomic imprinting represents an original model of epigenetic regulation resulting in a parent-of-origin expression. Despite the critical role of imprinted genes in mammalian growth, metabolism and neuronal function, there is no molecular tool specifically targeting them for a systematic evaluation. We show here that enzymatic methyl-seq consistently outperforms the bisulfite-based standard in capturing 165 candidate regions for genomic imprinting in the pig. This highlights the potential for a turnkey, fully customizable and reliable capture tool of genomic regions regulated by cytosine methylation in any population of interest. For the field of genomic imprinting, it opens up the possibility of detecting multilocus imprinting variations across the genome, with implications for basic research, agrigenomics and clinical practice.

Driving gut microbiota enterotypes through host genetics.

BACKGROUND: Population stratification based on interindividual variability in gut microbiota composition has revealed the existence of several ecotypes named enterotypes in humans and various animal species. Enterotypes are often associated with environmental factors including diet, but knowledge of the role of host genetics remains scarce. Moreover, enterotypes harbor functionalities likely associated with varying abilities and susceptibilities of their host. Previously, we showed that under controlled conditions, 60-day-old pig populations consistently split into two enterotypes with either Prevotella and Mitsuokella (PM enterotype) or Ruminococcus and Treponema (RT enterotype) as keystone taxa. Here, our aim was to rely on pig as a model to study the influence of host genetics to assemble enterotypes, and to provide clues on enterotype functional differences and their links with growth traits. RESULTS: We established two pig lines contrasted for abundances of the genera pairs specifying each enterotype at 60 days of age and assessed them for fecal microbiota composition and growth throughout three consecutive generations. Response to selection across three generations revealed, per line, an increase in the prevalence of the selected enterotype and in the average relative abundances of directly and indirectly selected bacterial genera. The PM enterotype was found less diverse than the RT enterotype but more efficient for piglet growth during the post-weaning period. Shotgun metagenomics revealed differentially abundant bacterial species between the two enterotypes. By using the KEGG Orthology database, we show that functions related to starch degradation and polysaccharide metabolism are enriched in the PM enterotype, whereas functions related to general nucleoside transport and peptide/nickel transport are enriched in the RT enterotype. Our results also suggest that the PM and RT enterotypes might differ in the metabolism of valine, leucin, and isoleucine, favoring their biosynthesis and degradation, respectively. CONCLUSION: We experimentally demonstrated that enterotypes are functional ecosystems that can be selected as a whole by exerting pressure on the host genetics. We also highlight that holobionts should be considered as units of selection in breeding programs. These results pave the way for a holistic use of host genetics, microbiota diversity, and enterotype functionalities to understand holobiont shaping and adaptation. Video Abstract.

Secreted phosphoprotein 1 expression in endometrium and placental tissues of hyperprolific large white and meishan gilts.

Increased litter size and within-litter uniformity in birth weight would improve pig reproductive efficiency. This study compared the location and gene and protein expression of secreted phosphoprotein 1 in placental and uterine tissues supplying a normally sized and the smallest fetus carried by hyperprolific Large White and Meishan gilts on Days 41-42 of pregnancy. Immunohistochemistry and in situ hybridization showed that the protein and gene encoding secreted phosphoprotein 1 were located in the glandular and luminal epithelium of the endometrium and in the placenta. Secreted phosphoprotein 1 protein levels were higher in glandular epithelium, luminal epithelium, and placenta from Meishan gilts compared to corresponding tissues from hyperprolific Large White gilts. Reverse transcription quantitative PCR demonstrated secreted phosphoprotein 1 mRNA levels were higher in endometrium, but not placenta, from Meishan compared to hyperprolific Large White gilts. In hyperprolific Large White gilts, secreted phosphoprotein 1 protein levels were higher in glandular epithelium and placenta surrounding small fetuses than corresponding tissues supplying normal-sized fetuses. Similarly, in Meishan gilts, secreted phosphoprotein 1 protein levels were higher in luminal epithelium surrounding small compared to normal-sized fetuses. Within hyperprolific Large White, but not Meishan, gilts secreted phosphoprotein 1 mRNA was higher in endometrium surrounding the normal-sized fetus than the control fetus. The contradictory relationship between fetal size and secreted phosphoprotein 1 protein and mRNA in the hyperprolific Large White is intriguing and may reflect breed differences in posttranslational modification. The striking breed differences in secreted phospoprotein 1 expression suggest that SPP1 may be associated with placental efficiency.

Plasma 1H-NMR metabolic and amino acid profiles of newborn piglets from two lines divergently selected for residual feed intake.

Together with environmental factors, physiological maturity at birth is a major determinant for neonatal survival and postnatal development in mammalian species. Maturity at birth is the outcome of complex mechanisms of intra-uterine development and maturation during the end of gestation. In pig production, piglet preweaning mortality averages 20% of the litter and thus, maturity is a major welfare and economic concern. Here, we used both targeted and untargeted metabolomic approaches to provide a deeper understanding of the maturity in a model of lines of pigs divergently selected on residual feed intake (RFI), previously shown to have contrasted signs of maturity at birth. Analyses were conducted on plasma metabolome of piglets at birth and integrated with other phenotypic characteristics associated to maturity. We confirmed proline and myo-inositol, previously described for their association with delayed growth, as potential markers of maturity. Urea cycle and energy metabolism were found more regulated in piglets from high and low RFI lines, respectively, suggesting a better thermoregulation ability for the low RFI (with higher feed efficiency) piglets.

Rabbit targeted genomic sequences after heterologous hybridization using human exome.

OBJECTIVE: Causal mutations for major genes that underlie a broad range of morphological traits are often located within exons of genes that then affect protein functions. Non-model organism genetic studies are not easy to perform due to the lack of genome-wide molecular tools such as SNP genotyping array. Genotyping-By-Sequencing (GBS) methods offer an alternative. Consequently, we used this approach that is focused on the exome to target and identify major genes in rabbit populations. Data description We used a heterologous enrichment method before sequencing, allowing us to capture the rabbit exome using the marketed human panel since mammal protein coding genes are well conserved across the phylogenic tree of species. This targeted strategy was performed on 52 French rabbits from 5 different French strains (Californian, New-Zealand, Castor, Chinchilla and Laghmere). We generated 3.4 billion sequencing reads and approximately 29-140 million of reads per DNA sample. The expected exome coverage per sample ranged between 118 and 566X. The present dataset could be useful for the scientific community working on rabbit species in order to (i) improve the annotation of the rabbit reference genome Oryctolagus cuniculus (OryCun2.0), (ii) enrich the characterization of polymorphisms segregating in rabbits and (iii) evaluate the genetic biodiversity in different rabbit strains. Raw sequences were deposited in the European Nucleotide Archive (ENA) at the European Molecular Biology Laboratory- European Bioinformatics Institute (EMBL-EBI) data portal under bioproject accession number PRJEB37917.

Analysis of image-based sow activity patterns reveals several associations with piglet survival and early growth.

An activity pattern describes variations in activities over time. The objectives of this study are to automatically predict sow activity from computer vision over 11 days peripartum and estimate how sow behavior influences piglet's performance during early lactation. The analysis of video images used the convolutional neural network (CNN) YOLO for sow detection and posture classification of 21 Large White and 22 Meishan primiparous sows housed in individual farrowing pens. A longitudinal analysis and a clustering method were combined to identify groups of sows with a similar activity pattern. Traits under study are as follows: (i) the distribution of time spent daily in different postures and (ii) different activities while standing. Six postures were included along with three classes of standing activities, i.e., eating, drinking, and other, which can be in motion or not and root-pawing or not. They correspond to a postural budget and a standing-activity budget. Groups of sows with similar changes in their budget over the period (D-3 to D-1; D0 and D1-D7) were identified with the k-means clustering method. Next, behavioral traits (time spent daily in each posture, frequency of postural changes) were used as explanatory variables in the Cox proportional hazards model for survival and in the linear model for growth. Piglet survival was influenced by sow behavior on D-1 and during the period D1-D7. Piglets born from sows that were standing and doing an activity other than drinking and eating on D-1 had a 26% lower risk of dying than other piglets. Those born from sows that changed posture more frequently on D1-D7 had a 44% lower risk of dying. The number of postural changes, which illustrate sow restlessness, influenced piglet growth in the three periods. The average daily gain of piglets born from sows that were more restless on D1-D7 and that changed posture more frequently to hide their udder on D0 decreased by 22 and 45 g/d, respectively. Conversely, those born from sows that changed posture more frequently to hide their udder during the period of D1-D7 grew faster (+71 g/d) than the other piglets. Sow restlessness at different time periods influenced piglet performance.