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.

High-throughput phenotyping to characterise range use behaviour in broiler chickens

A key characteristic of free-range chicken farming is to enable chickens to spend time outdoors. However, each chicken may use the available areas for roaming in variable ways. To check if, and how, broilers use their outdoor range at an individual level, we need to reliably characterise range use behaviour. Traditional methods relying on visual scans require significant time investment and only provide discontinuous information. Passive RFID (Radio Frequency Identification) systems enable tracking individually tagged chickens' when they go through pop-holes; hence, they only provide partial information on the movements of individual chickens. Here, we describe a new method to measure chickens' range use and test its reliability on three ranges each containing a different breed. We used an active RFID system to localise chickens in their barn, or in one of nine zones of their range, every 30 seconds and assessed range-use behaviour in 600 chickens belonging to three breeds of slow- or medium-growing broilers used for outdoor production (all <40 g daily weight gain). From those real-time locations, we determined five measures to describe daily range use: time spent in the barn, number of outdoor accesses, number of zones visited in a day, gregariousness (an index that increases when birds spend time in zones where other birds are), and numbers of zone changes. Principal Component Analyses (PCAs) were performed on those measures, in each production system, to create two synthetic indicators of chickens' range use behaviour. The first two PCA axes represented over 90% of the variance of the five measures and were both consistent over time and correlated with independent visual scans. Contributions of the five measures to the PCAs were similar among breeds, except for the correlation between the number of outdoor accesses and the four other measures. PC1 correlated with time spent inside the barn and zone changes frequency, whilst PC2 was explained by exploration of the range. Taken together, PC1 and PC2 indicators showed that range use increased with age, outdoor temperature (in spring), and did not differ between males and females. Importantly, daily scores for both indicators were repeatable among individuals - particularly in PC1 - showing inter-individual variability on range-use. The characterisation of broiler behaviour around their range with these reliable and repeatable indicators provides novel tools to help understand individual variations of range-use in free-range farming.

Links between fecal microbiota and the response to vaccination against influenza A virus in pigs

This study describes the associations between fecal microbiota and vaccine response variability in pigs, using 98 piglets vaccinated against the influenza A virus at 28 days of age (D28) with a booster at D49. Immune response to the vaccine is measured at D49, D56, D63, and D146 by serum levels of IAV-specific IgG and assays of hemagglutination inhibition (HAI). Analysis of the pre-vaccination microbiota characterized by 16S rRNA gene sequencing of fecal DNA reveals a higher vaccine response in piglets with a richer microbiota, and shows that 23 operational taxonomic units (OTUs) are differentially abundant between high and low IAV-specific IgG producers at D63. A stronger immune response is linked with OTUs assigned to the genus Prevotella and family Muribaculaceae, and a weaker response is linked with OTUs assigned to the genera Helicobacter and Escherichia-Shigella. A set of 81 OTUs accurately predicts IAV-specific IgG and HAI titer levels at all time points, highlighting early and late associations between pre-vaccination fecal microbiota composition and immune response to the vaccine.

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.

Thermoregulation at birth differs between piglets from two genetic lines divergent for residual feed intake

Thermoregulation is essential to piglets' neonatal survival. This study used infrared thermography (IRT) to assess thermoregulation abilities of piglets from two lines divergent for residual feed intake (RFI). At birth, morphology (weight, length, width and circumference), vigour (respiration, mobility and vocalisation), and rectal temperature were recorded from piglets of the 11th generation of the low RFI (LRFI, more efficient; n = 34) and the high RFI (HRFI, less efficient; n = 28) lines. Infrared thermography images were taken at 8, 15, 30 and 60 min post partum. Temperatures of the ear base and tip, and of the back (i.e. shoulders to rumps) were extracted (Thermacam Researcher Pro 2.0) and analysed with linear mixed models (SAS 9.4). Piglets had different average hourly weight gain (HRFI = 7.1 ± 1.3 g/h, LRFI = 3.6 ± 1.3 g/h; P < 0,001) but did not differ in morphology or vigour. All temperatures increased overtime. At birth, piglets' rectal temperature was correlated with the initial temperature of the ear base and the maximum back temperature (0.37 and 0.33, respectively; P < 0.05). High residual feed intake piglets had lower ear tip temperatures than LRFI piglets at 15 (24.7 ± 0.37 °C vs. 26.3 ± 0.36 °C, respectively; F1, 63.5 = 9.11, P < 0.005) and 30 min post partum (26.2 ± 0.47 °C vs. 27.6 ± 0.44 °C, respectively; F1, 66.9 = 4.52, P < 0.05). Moreover, thermal pattern of the ear tip differed between the two genetic lines. In conclusion, IRT allowed non-invasive assessment of piglets' thermoregulation abilities and indicated an influence of genetic selection for RFI on neonatal thermoregulation abilities.

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.