RNA-seq of muscle from pigs divergent in feed efficiency and product quality identifies differences in immune response, growth, and macronutrient and connective tissue metabolism

Genes involved in the cholecystokinin receptor signaling map were differentially expressed in the jejunum of steers with variation in residual feed intake

The jejunum is a critical site for nutrient digestion and absorption, and variation in its ability to take up nutrients within the jejunum is likely to affect feed efficiency. The purpose of this study was to determine differences in gene expression in the jejunum of beef steers divergent for residual feed intake (RFI) in one cohort of steers (Year 1), and to validate those genes in animals from a second study (Year 2). Steers from Year 1 (n = 16) were selected for high and low RFI. Jejunum mucosal tissue was obtained for RNA-seq. Thirty-two genes were differentially expressed (PFDR≤0.15), and five were over-represented in pathways including inflammatory mediator, cholecystokinin receptor (CCKR) signaling, and p38 MAPK pathways. Several differentially expressed genes (ALOX12, ALPI, FABP6, FABP7, FLT1, GSTA2, MEF2B, PDK4, SPP1, and TTF2) have been previously associated with RFI in other studies. Real-time qPCR was used to validate nine differentially expressed genes in the Year 1 steers used for RNA-seq, and in the Year 2 validation cohort. Six genes were validated as differentially expressed (P < 0.1) using RT-qPCR in the Year 1 population. In the Year 2 population, five genes displayed the same direction of expression as the Year 1 population and 3 were differentially expressed (P < 0.1). The CCKR pathway is involved in digestion, appetite control, and regulation of body weight making it a compelling candidate for feed efficiency in cattle, and the validation of these genes in a second population of cattle is suggestive of a role in feed efficiency.

Population dynamics of potentially harmful haplotypes: a pedigree analysis

BACKGROUND

The identification of low-frequency haplotypes, never observed in homozygous state in a population, is considered informative on the presence of potentially harmful alleles (candidate alleles), putatively involved in inbreeding depression. Although identification of candidate alleles is challenging, studies analyzing the dynamics of potentially harmful alleles are lacking. A pedigree of the highly endangered Gochu Asturcelta pig breed, including 471 individuals belonging to 51 different families with at least 5 offspring each, was genotyped using the Axiom PigHDv1 Array (658,692 SNPs). Analyses were carried out on four different cohorts defined according to pedigree depth and at the whole population (WP) level.

RESULTS

The 4,470 Linkage Blocks (LB) identified in the Base Population (10 individuals), gathered a total of 16,981 alleles in the WP. Up to 5,466 (32%) haplotypes were statistically considered candidate alleles, 3,995 of them (73%) having one copy only. The number of alleles and candidate alleles varied across cohorts according to sample size. Up to 4,610 of the alleles identified in the WP (27% of the total) were present in one cohort only. Parentage analysis identified a total of 67,742 parent-offspring incompatibilities. The number of mismatches varied according to family size. Parent-offspring inconsistencies were identified in 98.2% of the candidate alleles and 100% of the LB in which they were located. Segregation analyses informed that most potential candidate alleles appeared de novo in the pedigree. Only 17 candidate alleles were identified in the boar, sow, and paternal and maternal grandparents and were considered segregants.

CONCLUSIONS

Our results suggest that neither mutation nor recombination are the major forces causing the apparition of candidate alleles. Their occurrence is more likely caused by Allele-Drop-In events due to SNP calling errors. New alleles appear when wrongly called SNPs are used to construct haplotypes. The presence of candidate alleles in either parents or grandparents of the carrier individuals does not ensure that they are true alleles. Minimum Allele Frequency thresholds may remove informative alleles. Only fully segregant candidate alleles should be considered potentially harmful alleles. A set of 16 candidate genes, potentially involved in inbreeding depression, is described.

Candidate genes for longitudinal traits under sequential sampling in beef cattle

Both the measurement age of a longitudinal trait and the common pre-sampling procedures used in beef cattle herds may affect the identification of a functional candidate gene (FCG) that is potentially associated with a trait. To identify the FCG that takes part in the genetic control of body weight at five different ages in a beef cattle population with and without sequential sampling, the animals were weighed at different measurement events, around 330, 385, 440, 495 and 550 days old. Genetic parameters were estimated for body weight at each age using a single trait (STM) and a random regression model (RRM). In addition, two different databases were used to estimate the genetic parameters: the first (DB100) was formed by all animals that were weighed in the five measurement events, and the second (DB70) has records of the same population, considering that 70% of the heaviest animals were selected after each measurement event. For DB100, genome-wide association studies (GWAS) were performed with 21,667 SNP markers to identify genomic windows that explained at least 1% of the genetic variance. Additionally, prioritization analyses were performed and FCGs were selected. We associated seven different FCGs with body weight at different ages. Among them, the gene DUSP10 was suggested as FCG in all five ages evaluated. Genetic parameters estimated for body weight using DB100 were similar when STM and RRM were applied. However, when DB70 was used as phenotypic data, there were differences between the two models. When the STM was applied, there were differences between the genetic parameters estimated for body weight when DB100 or DB70 were used as sources of phenotypes, but not for the estimates obtained with RRM. The importance of each gene for animal growth can change at different ages, and different genes may be more relevant to body weight at each different growth stage for beef cattle. Besides, sequential sampling can affect the GWAS results of a longitudinal trait. The age of the animal when a longitudinal trait is measured and pre-sampling can also contribute to inconsistencies in GWAS results for body weight in beef cattle, depending on the time when that data were collected, and consequently on the identification of FCG between studies, even when models that consider a covariance structure are used.

Molecular mechanisms contributing to the development of beef sensory texture and flavour traits and related biomarkers: Insights from early post-mortem muscle using label-free proteomics

Beef sensory quality comprises a suite of traits, each of which manifests its ultimate phenotype through interaction of muscle physiology with environment, both in vivo and post-mortem. Understanding variability in meat quality remains a persistent challenge, but omics studies to uncover biological connections between natural variability in proteome and phenotype could provide validation for exploratory studies and offer new insights. Multivariate analysis of proteome and meat quality data from Longissimus thoracis et lumborum muscle samples taken early post-mortem from 34 Limousin-sired bulls was conducted. Using for the first-time label-free shotgun proteomics combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS), 85 proteins were found to be related with tenderness, chewiness, stringiness and flavour sensory traits. The putative biomarkers were classified in five interconnected biological pathways; i) muscle contraction, ii) energy metabolism, iii) heat shock proteins, iv) oxidative stress, v) regulation of cellular processes and binding. Among the proteins, PHKA1 and STBD1 correlated with all four traits, as did the GO biological process 'generation of precursor metabolites and energy'. Optimal regression models explained a high level (58-71%) of phenotypic variability with proteomic data for each quality trait. The results of this study propose several regression equations and biomarkers to explain the variability of multiple beef eating quality traits. Thanks to annotation and network analyses, they further suggest protein interactions and mechanisms underpinning the physiological processes regulating these key quality traits. SIGNIFICANCE: The proteomic profiles of animals with divergent quality profiles have been compared in numerous studies; however, a wide range of phenotypic variation is required to better understand the mechanisms underpinning the complex biological pathways correlated with beef quality and protein interactions. We used multivariate regression analyses and bioinformatics to analyse shotgun proteomics data to decipher the molecular signatures involved in beef texture and flavour variations with a focus on multiple quality traits. We developed multiple regression equations to explain beef texture and flavour. Additionally, potential candidate biomarkers correlated with multiple beef quality traits are suggested, which could have utility as indicators of beef overall sensory quality. This study explained the biological process responsible for determining key quality traits such as tenderness, chewiness, stringiness, and flavour in beef, which will provide support for future beef proteomics studies.

Molecular Mechanism of m6A Methylation Modification Genes METTL3 and FTO in Regulating Heat Stress in Sheep

Heat stress is an important environmental factor affecting livestock production worldwide. Primary hepatocytes and preadipocytes derived from Hu sheep were used to establish a heat stress model. Quantitative reverse transcriptase-PCR (qRT-PCR) analysis showed that heat induction significantly increased the expression levels of heat stress protein (HSP) genes and the N6-methyladenosine (m6A) methylation modification genes: methyltransferase-like protein 3 (METTL3), methyltransferase-like protein 14 (METTL14), and fat mass and obesity associated protein (FTO). Heat stress simultaneously promoted cell apoptosis. Transcriptome sequencing identified 3980 upregulated genes and 2420 downregulated genes related to heat stress. A pathway enrichment analysis of these genes revealed significant enrichment in fatty acid biosynthesis, degradation, and the PI3K-Akt and peroxisome proliferator-activated receptor (PPAR) signaling pathways. Overexpression of METTL3 in primary hepatocytes led to significant downregulation of HSP60, HSP70, and HSP110, and significantly increased mRNA m6A methylation; FTO interference generated the opposite results. Primary adipocytes showed similar results. Transcriptome analysis of cells under METTL3 (or FTO) inference and overexpression revealed differentially expressed genes enriched in the mitogen-activated protein kinase (MAPK) signaling pathways, as well as the PI3K-Akt and Ras signaling pathways. We speculate that METTL3 may increase the level of m6A methylation to inhibit fat deposition and/or inhibit the expression of HSP genes to enhance the body's resistance to heat stress, while the FTO gene generated the opposite molecular mechanism. This study provides a scientific basis and theoretical support for sheep feeding and management practices during heat stress.

Exon-intron split analysis reveals posttranscriptional regulatory signals induced by high and low n-6/n-3 polyunsaturated fatty acid ratio diets in piglets

Polyunsaturated fatty acids (PUFA), such as omega-6 (n-6) and omega-3 (n-3), play a vital role in nutrient metabolism, inflammatory response, and gene regulation. microRNAs (miRNA), which can potentially degrade targeted messenger RNAs (mRNA) and/or inhibit their translation, might play a relevant role in PUFA-related changes in gene expression. Although differential expression analyses can provide a comprehensive picture of gene expression variation, they are unable to disentangle when in the mRNA life cycle the regulation of expression is taking place, including any putative functional miRNA-driven repression. To capture this, we used an exon-intron split analysis (EISA) approach to account for posttranscriptional changes in response to extreme values of n-6/n-3 PUFA ratio. Longissimus dorsi muscle samples of male and female piglets from sows fed with n-6/n-3 PUFA ratio of 13:1 (SOY) or 4:1 (LIN), were analyzed in a bidirectional contrast (LIN vs. SOY, SOY vs. LIN). Our results allowed the identification of genes showing strong posttranscriptional downregulation signals putatively targeted by significantly upregulated miRNA. Moreover, we identified genes primarily involved in the regulation of lipid-related metabolism and immune response, which may be associated with the pro- and anti-inflammatory functions of the n-6 and n-3 PUFA, respectively. EISA allowed us to uncover regulatory networks complementing canonical differential expression analyses, thus providing a more comprehensive view of muscle metabolic changes in response to PUFA concentration.

Hygiene of housing conditions and proinflammatory signals alter gene expressions in porcine adipose tissues and blood cells

Adipose tissue is an organ with metabolic, endocrine and immune functions. In this tissue, the expressions of genes associated with several metabolic pathways, including lipid metabolism, have been shown to be affected by genetic selection for feed efficiency, an important trait to consider in livestock. We hypothesized that the stimulation of immune system caused by poor hygiene conditions of housing impacts the molecular and cellular features of adipose tissue and that the impact may differ between pigs that diverge in feed efficiency. At the age of 12 weeks, Large White pigs from two genetic lines divergent for residual feed intake (RFI) were housed in two contrasting hygiene conditions (good vs poor). After six weeks of exposure, pigs were slaughtered (n = 36). Samples of blood, subcutaneous (SCAT) and perirenal (PRAT) adipose tissues were collected for cell response and gene expression investigations. The decrease in the relative weight of PRAT was associated with a decline in mRNA levels of FASN, ME, LCN2 and TLR4 (P < 0.05) in pigs housed in poor conditions compared with pigs housed in good conditions for both RFI lines. In SCAT, the expressions of only two key genes (PPARG and TLR4) were significantly affected by the hygiene of housing conditions. Besides, the mRNA levels of both LCN2 and GPX3 were influenced by the RFI line (P < 0.05). Because we suspected an effect of poor hygiene at the cellular levels, we investigated the differentiation of stromal vascular cells isolated from SCAT in vitro in the absence or presence of a pro-inflammatory cytokine, Tumor Necrosis Factor-α (TNF-α). The ability of these cells to differentiate in the absence or presence of TNF-α did not differ among the four groups of animals (P > 0.05). We also investigated the expressions of genes involved in the immune response and lipid metabolism in whole blood cells cultured in the absence and presence of LPS. The hygiene conditions had no effect but, the relative expression of the GPX3 gene was higher (P < 0.001) in high RFI than in low RFI pigs while the expressions of IL-10 (P = 0.027), TGFβ1 (P = 0.023) and ADIPOR2 (P = 0.05) genes were lower in high RFI than in low RFI pigs. Overall, the current study indicates that the hygiene of housing had similar effects on both RFI lines on the expression of genes in adipose tissues and on the features of SCAT adipose cells and whole blood cells in response to TNF-α and LPS. It further demonstrates that the number of genes with expression impacted by housing conditions was higher in PRAT than in SCAT. It suggests a depot-specific response of adipose tissue to the current challenge.

Modeling microRNA-driven post-transcriptional regulation using exon-intron split analysis in pigs

The contribution of microRNAs (miRNAs) to mRNA post-transcriptional regulation has often been explored by the post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes in mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Göttingen minipigs. Exon-intron split analysis of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of miRNAs upregulated in fed pigs. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, exon-intron split analysis provides a valuable view of the regulation of gene expression, complementary to canonical differential expression analyses.

Application of Genetic, Genomic and Biological Pathways in Improvement of Swine Feed Efficiency

Despite the significant improvement of feed efficiency (FE) in pigs over the past decades, feed costs remain a major challenge for producers profitability. Improving FE is a top priority for the global swine industry. A deeper understanding of the biology underlying FE is crucial for making progress in genetic improvement of FE traits. This review comprehensively discusses the topics related to the FE in pigs including: measurements, genetics, genomics, biological pathways and the advanced technologies and methods involved in FE improvement. We first provide an update of heritability for different FE indicators and then characterize the correlations of FE traits with other economically important traits. Moreover, we present the quantitative trait loci (QTL) and possible candidate genes associated with FE in pigs and outline the most important biological pathways related to the FE traits in pigs. Finally, we present possible ways to improve FE in swine including the implementation of genomic selection, new technologies for measuring the FE traits, and the potential use of genome editing and omics technologies.

Poor hygiene of housing conditions influences energy metabolism in a muscle type-dependent manner in growing pigs differing in feed efficiency

The ability of pigs to cope with inflammatory challenges may by modified by selection for residual feed intake (RFI), a measure of feed efficiency. In the current study, we evaluated skeletal muscle metabolic responses to degraded hygiene conditions in pigs divergently selected for RFI. At 82 d of age, low RFI and high RFI pigs were housed in either poor or good hygiene conditions. After a 6-week challenge, the poor hygiene conditions induced a decrease in growth performance (P < 0.001) and in plasma IGF-I concentrations (P < 0.003) in both lines. In the slow-twitch oxidative semispinalis muscle, poor hygiene conditions induced a shift towards a more oxidative metabolism and an activation of the AMPK pathway in pigs of both RFI lines. In the fast-twitch glycolytic longississimus muscle, poor hygiene conditions were associated to a less glycolytic metabolism in the HRFI line only. Poor hygiene conditions also increased the protein level of lipidation of microtubule-associated protein 1 light-chain 3β (LC3-II) in both RFI lines, suggesting an activation of the autophagy pathway. Altogether, the data revealed muscle-type specific metabolic adaptations to poor hygiene conditions, which may be related to different strategies to fuel the activated immune system.

Comparative proteomics revealed duodenal metabolic function associated with feed efficiency in slow-growing chicken

The Korat chicken (KR), developed in Thailand, is a slow-growing breed developed as an alternative breed for Thai chicken producers. The growing interest in slow-growing chicken meat, due to its unique taste, distinct texture, health benefits, and higher broiler welfare have led to higher market demand for KR. However, its low feed efficiency (FE) has a significant negative impact on farm profitability. Understanding the molecular mechanism regulating FE allows for designing a suitable selection program and contributing to breeding more efficient chicken for poultry production. Thus, the objective of our study was to investigate the proteome differences and possible pathways associated with FE in male KR using a label-free quantitative proteomic approach. Seventy-five KR males were individually evaluated for FE, and duodenum samples from 6 animals (3 high-FE and 3 low-FE chickens) were collected at 10 wk of age for differential abundant proteins (DAPs), protein networks, functional enrichment, and pathway analyses. In this study, we found 40 DAPs significantly associated with FE pathways, including glycolysis/gluconeogenesis, peroxisome, oxidative phosphorylation, tight junction, and cysteine and methionine metabolism. Thus, variations in observed DAPs or genes related to DAPs could be interesting biomarker candidates for selection for higher feed utilization efficiency in chicken.

Large-scale association study on daily weight gain in pigs reveals overlap of genetic factors for growth in humans

Background

Imputation from genotyping array to whole-genome sequence variants using resequencing of representative reference populations enhances our ability to map genetic factors affecting complex phenotypes in livestock species. The accumulation of knowledge about gene function in human and laboratory animals can provide substantial advantage for genomic research in livestock species.

Results

In this study, 201,388 pigs from three commercial Danish breeds genotyped with low to medium (8.5k to 70k) SNP arrays were imputed to whole genome sequence variants using a two-step approach. Both imputation steps achieved high accuracies, and in total this yielded 26,447,434 markers on 18 autosomes. The average estimated imputation accuracy of markers with minor allele frequency ≥ 0.05 was 0.94. To overcome the memory consumption of running genome-wide association study (GWAS) for each breed, we performed within-breed subpopulation GWAS then within-breed meta-analysis for average daily weight gain (ADG), followed by a multi-breed meta-analysis of GWAS summary statistics. We identified 15 quantitative trait loci (QTL). Our post-GWAS analysis strategy to prioritize of candidate genes including information like gene ontology, mammalian phenotype database, differential expression gene analysis of high and low feed efficiency pig and human GWAS catalog for height, obesity, and body mass index, we proposed MRAP2, LEPROT, PMAIP1, ENSSSCG00000036234, BMP2, ELFN1, LIG4 and FAM155A as the candidate genes with biological support for ADG in pigs.

Conclusion

Our post-GWAS analysis strategy helped to identify candidate genes not just by distance to the lead SNP but also by multiple sources of biological evidence. Besides, the identified QTL overlap with genes which are known for their association with human growth-related traits. The GWAS with this large data set showed the power to map the genetic factors associated with ADG in pigs and have added to our understanding of the genetics of growth across mammalian species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08373-3.

Jejunal Transcriptomic Profiling for Differences in Feed Conversion Ratio in Slow-Growing Chickens

Simple Summary

The slow-growing Korat chicken (KR) is economically attractive, as KR meat has a high selling price and has thus been used in Thailand to support smallholder farmers. However, low feed efficiency in KR stockbreeding makes the product less competitive and improving KR feed efficiency is central to increasing KR profitability. Using RNA sequencing, we compared the jejunal transcriptomic profiles of low- and high-feed conversion ratio (FCR) KR chickens, to identify FCR-related transcriptional variation and biological pathways. Gene Ontology and Kyoto Encyclopedia of Gene and Genome analysis revealed that the main pathways involved in KR FCR variation are related to immune response, glutathione metabolism, vitamin transport and metabolism, lipid metabolism, and neuronal and cardiac maturation, development, and growth. This is the first study to investigate, in the jejunum, the molecular genetic mechanisms affecting the FCR of slow-growing chickens. These findings will be useful in line-breeding programs to improve feed efficiency and profitability in slow-growing chicken stockbreeding.

Abstract

Improving feed efficiency is an important breeding target for the poultry industry; to achieve this, it is necessary to understand the molecular basis of feed efficiency. We compared the jejunal transcriptomes of low- and high-feed conversion ratio (FCR) slow-growing Korat chickens (KRs). Using an original sample of 75 isolated 10-week-old KR males, we took jejunal samples from six individuals in two groups: those with extremely low FCR (n = 3; FCR = 1.93 ± 0.05) and those with extremely high FCR (n = 3; FCR = 3.29 ± 0.06). Jejunal transcriptome profiling via RNA sequencing revealed 56 genes that were differentially expressed (p < 0.01, FC > 2): 31 were upregulated, and 25 were downregulated, in the low-FCR group relative to the high-FCR group. Functional annotation revealed that these differentially expressed genes were enriched in biological processes related to immune response, glutathione metabolism, vitamin transport and metabolism, lipid metabolism, and neuronal and cardiac maturation, development, and growth, suggesting that these are important mechanisms governing jejunal feed conversion. These findings provide an important molecular basis for future breeding strategies to improve slow-growing chicken feed efficiency.

GWAS and Post-GWAS High-Resolution Mapping Analyses Identify Strong Novel Candidate Genes Influencing the Fatty Acid Composition of the Longissimus dorsi Muscle in Pigs

Fatty acid (FA) composition is one of the most important parameters for the assessment of meat quality in pigs. The FA composition in pork can also affect human health. Our aim was to identify quantitative trait loci (QTLs) and positional candidate genes affecting the FA profile of the longissimus dorsi muscle in a large F2 intercross between Landrace and Korean native pigs comprising 1105 F2 progeny by genome-wide association studies (GWAS) and post-GWAS high-resolution mapping analyses. We performed GWAS using the PorcineSNP60K BeadChip and a linear mixed model. Four genome-wide significant QTL regions in SSC8, SSC12, SSC14, and SSC16 were detected (p < 2.53 × 10⁻⁷). Several co-localizations of QTLs in SSC12 for oleic acid, linoleic acid, arachidonic acid, monounsaturated FAs, polyunsaturated FAs, and the polyunsaturated/saturated FA ratio were observed. To refine the QTL region in SSC12, a linkage and linkage disequilibrium analysis was applied and could narrow down the critical region to a 0.749 Mb region. Of the genes in this region, GAS7, MYH2, and MYH3 were identified as strong novel candidate genes based on further conditional association analyses. These findings provide a novel insight into the genetic basis of FA composition in pork and could contribute to the improvement of pork quality.

Biological mechanisms of growth performance and meat quality in porcine muscle tissue

Growth performance and meat quality are important traits for pig production. The aim of the present study was to investigate the molecular mechanisms underlying growth performance and meat quality, and to identify novel target molecules for predicting the growth performance and meat quality. The differentially expressed genes (DEGs) in Diannan small ears pigs (DSP) and Landrace pigs (LP) were assessed by RNA-sequencing analyzing technology. A total of 339 DEGs were obtained between DSP and LP. 146 DEGs were upregulated in LP compared with DSP and 193 DEGs were upregulated in DSP compared with LP. The DEGs were significantly enriched in 26 GO and 3 KEGG pathways. The protein-protein interaction (PPI) network with 201 nodes and 382 edges was constructed and 5 modules were extracted from the entire network. The identified upregulated expression of genes involved in glycolysis and myogenesis as well as extracellular matrix may be associated with fast body and muscle deposition rates in LP. Increased expression of genes involved in PPAR signaling pathway and fatty acid metabolism as well as oxidative phosphate processes could be related to the intramuscular fat deposition and meat quality in DSP. The present study may provide an improved understanding of the growth performance and meat quality.

Effects of mineral methionine hydroxy analog chelate in sow diets on epigenetic modification and growth of progeny

The study was conducted to determine the effects of mineral methionine hydroxy analog chelate (MMHAC) partially replacing inorganic trace minerals in sow diets on epigenetic and transcriptional changes in the muscle and jejunum of progeny. The MMHAC is zinc (Zn), manganese (Mn), and copper (Cu) chelated with methionine hydroxy analog (Zn-, Mn-, and Cu-methionine hydroxy analog chelate [MHAC]). On day 35 of gestation, 60 pregnant sows were allotted to two dietary treatments in a randomized completed block design using parity as a block: 1) ITM: inorganic trace minerals with zinc sulfate (ZnSO₄), manganese oxide (MnO), and copper sulfate (CuSO₄) and 2) CTM: 50% of ITM was replaced with MMHAC (MINTREX trace minerals, Novus International Inc., St Charles, MO). Gestation and lactation diets were formulated to meet or exceed NRC requirements. On days 1 and 18 of lactation, milk samples from 16 sows per treatment were collected to measure immunoglobulins (immunoglobulin G, immunoglobulin A, and immunoglobulin M) and micromineral concentrations. Two pigs per litter were selected to collect blood to measure the concentration of immunoglobulins in the serum, and then euthanized to collect jejunal mucosa, jejunum tissues, and longissimus muscle to measure global deoxyribonucleic acid methylation, histone acetylation, cytokines, and jejunal histomorphology at birth and day 18 of lactation. Data were analyzed using Proc MIXED of SAS. Supplementation of MMHAC tended to decrease (P = 0.059) body weight (BW) loss of sows during lactation and tended to increase (P = 0.098) piglet BW on day 18 of lactation. Supplementation of MMHAC increased (P < 0.05) global histone acetylation and tended to decrease myogenic regulatory factor 4 messenger ribonucleic acid (mRNA; P = 0.068) and delta 4-desaturase sphingolipid1 (DEGS1) mRNA (P = 0.086) in longissimus muscle of piglets at birth. Supplementation of MMHAC decreased (P < 0.05) nuclear factor kappa B mRNA in the jejunum and DEGS1 mRNA in longissimus muscle and tended to decrease mucin-2 (MUC2) mRNA (P = 0.057) and transforming growth factor-beta 1 (TGF-β1) mRNA (P = 0.057) in the jejunum of piglets on day 18 of lactation. There were, however, no changes in the amounts of tumor necrosis factor-alpha, interleukin-8, TGF-β, MUC2, and myogenic factor 6 in the tissues by MMHAC. In conclusion, maternal supplementation of MMHAC could contribute to histone acetylation and programming in the fetus, which potentially regulates intestinal health and skeletal muscle development of piglets at birth and weaning, possibly leading to enhanced growth of their piglets.

Brain Transcriptome Analysis Reveals Potential Transcription Factors and Biological Pathways Associated with Feed Efficiency in Commercial DLY Pigs

Feed efficiency (FE) is one of the most important economic traits in the porcine industry. In this study, high-throughput RNA sequencing (RNA-seq) was first utilized for brain tissue transcriptome analysis in pigs to indicate the potential genes and biological pathways related to FE in pigs. A total of 8 pigs with either extremely high-FE group (HE-group) or low-FE group (LE-group) were selected from 225 Duroc × (Landrace × Yorkshire) (DLY) pigs for transcriptomic analysis. RNA-seq analysis was performed to determine differentially expressed genes (DEGs) between the HE- and LE-group, and 430 DEGs were identified in brain tissues of pigs (|log₂(FoldChange)| > 1; adjusted p-values <0.05). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the DEGs were mainly enriched in synaptic signaling or transmission, and hormone secretion pathways, in which insulin secretion, and oxytocin signaling pathways were closely associated with FE by regulating feeding behavior and energy metabolism (adjusted p-values <0.05). Further, the transcription factors (TFs) analysis and gene co-expression network analysis indicated three hub differentially expressed TFs (NR2F2, TFAP2D, and HNF1B) that affected FE by mainly regulating feeding behavior, insulin sensitivity, or energy metabolism. Our findings suggest several potential TFs and biological pathways for further investigations of FE in pigs.

Mitochondrial protein gene expression and the oxidative phosphorylation pathway associated with feed efficiency and energy balance in dairy cattle

Feed efficiency and energy balance are important traits underpinning profitability and environmental sustainability in animal production. They are complex traits, and our understanding of their underlying biology is currently limited. One measure of feed efficiency is residual feed intake (RFI), which is the difference between actual and predicted intake. Variation in RFI among individuals is attributable to the metabolic efficiency of energy utilization. High RFI (H_RFI) animals require more energy per unit of weight gain or milk produced compared with low RFI (L_RFI) animals. Energy balance (EB) is a closely related trait calculated very similarly to RFI. Cellular energy metabolism in mitochondria involves mitochondrial protein (MiP) encoded by both nuclear (NuMiP) and mitochondrial (MtMiP) genomes. We hypothesized that MiP genes are differentially expressed (DE) between H_RFI and L_RFI animal groups and similarly between negative and positive EB groups. Our study aimed to characterize MiP gene expression in white blood cells of H_RFI and L_RFI cows using RNA sequencing to identify genes and biological pathways associated with feed efficiency in dairy cattle. We used the top and bottom 14 cows ranked for RFI and EB out of 109 animals as H_RFI and L_RFI, and positive and negative EB groups, respectively. The gene expression counts across all nuclear and mitochondrial genes for animals in each group were used for differential gene expression analyses, weighted gene correlation network analysis, functional enrichment, and identification of hub genes. Out of 244 DE genes between RFI groups, 38 were MiP genes. The DE genes were enriched for the oxidative phosphorylation (OXPHOS) and ribosome pathways. The DE MiP genes were underexpressed in L_RFI (and negative EB) compared with the H_RFI (and positive EB) groups, suggestive of reduced mitochondrial activity in the L_RFI group. None of the MtMiP genes were among the DE MiP genes between the groups, which suggests a non-rate limiting role of MtMiP genes in feed efficiency and warrants further investigation. The role of MiP, particularly the NuMiP and OXPHOS pathways in RFI, was also supported by our gene correlation network analysis and the hub gene identification. We validated the findings in an independent data set. Overall, our study suggested that differences in feed efficiency in dairy cows may be linked to differences in cellular energy demand. This study broadens our knowledge of the biology of feed efficiency in dairy cattle.

Differences in the Loin Tenderness of Iberian Pigs Explained through Dissimilarities in Their Transcriptome Expression Profile

Simple Summary

The Iberian pig is the most representative autochthonous breed of the Mediterranean region with unique genetic and phenotypic characteristics. The breed has been successfully preserved by its high-quality meat and high-priced products. Tenderness is one of the most relevant meat quality traits, and meat tenderization is influenced by genetic and environmental effects such as pre-slaughter handling and post-mortem conditions. Tenderness could be included in Iberian pig breeding programs, mainly focused on the improvement of premium-cuts percentage, in order to avoid the meat quality decline. A better biological understanding of this trait is needed. In the current study, we analyze the transcriptome of pigs divergent for Warner–Bratzler shear force through RNA-seq technique for the identification, characterization and quantification of candidate genes involved in biological pathways, networks and functions affecting meat tenderness.

Abstract

Tenderness is one of the most important meat quality traits and it can be measured through shear force with the Warner–Bratzler test. In the current study, we use the RNA-seq technique to analyze the transcriptome of Longissimus dorsi (LD) muscle in two groups of Iberian pigs (Tough and Tender) divergent for shear force breeding values. We identified 200 annotated differentially expressed genes (DEGs) and 245 newly predicted isoforms. The RNAseq expression results of 10 genes were validated with quantitative PCR (qPCR). Functional analyses showed an enrichment of DE genes in biological processes related to proteolysis (CTSC, RHOD, MYH8, ACTC1, GADD45B, CASQ2, CHRNA9 and ANKRD1), skeletal muscle tissue development (ANKRD1, DMD, FOS and MSTN), lipid metabolism (FABP3 and PPARGC1A) and collagen metabolism (COL14A1). The upstream analysis revealed a total of 11 transcription regulatory factors that could regulate the expression of some DEGs. Among them, IGF1, VGLL3 and PPARG can be highlighted since they regulate the expression of genes involved in biological pathways that could affect tenderness. The experiment revealed a set of candidate genes and regulatory factors suggestive to search polymorphisms that could be incorporated in a breeding program for improving meat tenderness.

Genetic variations (eQTLs) in muscle transcriptome and mitochondrial genes, and trans-eQTL molecular pathways in feed efficiency from Danish breeding pigs

Feed efficiency (FE) is a key trait in pig production, as improvement in FE has positive economic and environmental impact. FE is a complex phenotype and testing animals for FE is costly. Therefore, there has been a desire to find functionally relevant single nucleotide polymorphisms (SNPs) as biomarkers, to improve our biological understanding of FE as well as accuracy of genomic prediction for FE. We have performed a cis- and trans- eQTL (expression quantitative trait loci) analysis, in a population of Danbred Durocs (N = 11) and Danbred Landrace (N = 27) using both a linear and ANOVA model based on muscle tissue RNA-seq. We analyzed a total of 1425x19179 or 2.7x10⁷ Gene-SNP combinations in eQTL detection models for FE. The 1425 genes were from RNA-Seq based differential gene expression analyses using 25880 genes related to FE and additionally combined with mitochondrial genes. The 19179 SNPs were from applying stringent quality control and linkage disequilibrium filtering on genotype data using a GGP Porcine HD 70k SNP array. We applied 1000 fold bootstrapping and enrichment analysis to further validate and analyze our detected eQTLs. We identified 13 eQTLs with FDR < 0.1, affecting several genes found in previous studies of commercial pig breeds. Examples include MYO19, CPT1B, ACSL1, IER5L, CPT1A, SUCLA2, CSRNP1, PARK7 and MFF. The bootstrapping results showed statistically significant enrichment (p-value<2.2x10^-16) of eQTLs with p-value < 0.01 in both cis and trans-eQTLs. Enrichment analysis of top trans-eQTLs revealed high enrichment for gene categories and gene ontologies associated with genomic context and expression regulation. This included transcription factors (p-value = 1.0x10^-13), DNA-binding (GO:0003677, p-value = 8.9x10^-14), DNA-binding transcription factor activity (GO:0003700,) nucleus gene (GO:0005634, p-value<2.2x10^-16), negative regulation of expression (GO:0010629, p-value<2.2x10^-16). These results would be useful for future genome assisted breeding of pigs to improve FE, and in the improved understanding of the functional mechanism of trans eQTLs.

Evolutionary and Functional Analysis of Korean Native Pig Using Single Nucleotide Polymorphisms

Time and cost-effective production of next-generation sequencing data has enabled the performance of population-scale comparative and evolutionary studies for various species, which are essential for obtaining the comprehensive insight into molecular mechanisms underlying species- or breed-specific traits. In this study, the evolutionary and functional analysis of Korean native pig (KNP) was performed using single nucleotide polymorphism (SNP) data by comparative and population genomic approaches with six different mammalian species and five pig breeds. We examined the evolutionary history of KNP SNPs, and the specific genes of KNP based on the uniqueness of non-synonymous SNPs among the used species and pig breeds. We discovered the evolutionary trajectory of KNP SNPs within the used mammalian species as well as pig breeds. We also found olfaction-associated functions that have been characterized and diversified during evolution, and quantitative trait loci associated with the unique traits of KNP. Our study provides new insight into the evolution of KNP and serves as a good example for a better understanding of domestic animals in terms of evolution and domestication using the combined approaches of comparative and population genomics.

Genome Regulation and Gene Interaction Networks Inferred From Muscle Transcriptome Underlying Feed Efficiency in Pigs

Improvement of feed efficiency (FE) is key for Sustainability and cost reduction in pig production. Our aim was to characterize the muscle transcriptomic profiles in Danbred Duroc (Duroc; n = 13) and Danbred Landrace (Landrace; n = 28), in relation to FE for identifying potential biomarkers. RNA-seq data on the 41 pigs was analyzed employing differential gene expression methods, gene-gene interaction and network analysis, including pathway and functional analysis. We also compared the results with genome regulation in human exercise data, hypothesizing that increased FE mimics processes triggered in exercised muscle. In the differential expression analysis, 13 genes were differentially expressed, including: MRPS11, MTRF1, TRIM63, MGAT4A, KLH30. Based on a novel gene selection method, the divergent count, we performed pathway enrichment analysis. We found five significantly enriched pathways related to feed conversion ratio (FCR). These pathways were mainly related to mitochondria, and summarized in the mitochondrial translation elongation (MTR) pathway. In the gene interaction analysis, the most interesting genes included the mitochondrial genes: PPIF, MRPL35, NDUFS4 and the fat metabolism and obesity genes: AACS, SMPDL3B, CTNNBL1, NDUFS4, and LIMD2. In the network analysis, we identified two modules significantly correlated with FCR. Pathway enrichment of module genes identified MTR, electron transport chain and DNA repair as enriched pathways. The network analysis revealed the mitochondrial gene group NDUF as key network hub genes, showing their potential as biomarkers. Results show that genes related to human exercise were enriched in identified FCR related genes. We conclude that mitochondrial activity is a key driver for FCR in muscle tissue, and mitochondrial genes could be potential biomarkers for FCR in pigs.

Identification of key genes and pathways associated with feed efficiency of native chickens based on transcriptome data via bioinformatics analysis

BACKGROUND

Improving feed efficiency is one of the important breeding targets for poultry industry. The aim of current study was to investigate the breast muscle transcriptome data of native chickens divergent for feed efficiency. Residual feed intake (RFI) value was calculated for 1008 closely related chickens. The 5 most efficient (LRFI) and 5 least efficient (HRFI) birds were selected for further analysis. Transcriptomic data were generated from breast muscle collected post-slaughter.

RESULTS

The differently expressed genes (DEGs) analysis showed that 24 and 325 known genes were significantly up- and down-regulated in LRFI birds. An enrichment analysis of DEGs showed that the genes and pathways related to inflammatory response and immune response were up-regulated in HRFI chickens. Moreover, Gene Set Enrichment Analysis (GSEA) was also employed, which indicated that LRFI chickens increased expression of genes related to mitochondrial function. Furthermore, protein network interaction and function analyses revealed ND2, ND4, CYTB, RAC2, VCAM1, CTSS and TLR4 were key genes for feed efficiency. And the 'phagosome', 'cell adhesion molecules (CAMs)', 'citrate cycle (TCA cycle)' and 'oxidative phosphorylation' were key pathways contributing to the difference in feed efficiency.

CONCLUSIONS

In summary, a series of key genes and pathways were identified via bioinformatics analysis. These key genes may influence feed efficiency through deep involvement in ROS production and inflammatory response. Our results suggested that LRFI chickens may synthesize ATP more efficiently and control reactive oxygen species (ROS) production more strictly by enhancing the mitochondrial function in skeletal muscle compared with HRFI chickens. These findings provide some clues for understanding the molecular mechanism of feed efficiency in birds and will be a useful reference data for native chicken breeding.

Potential Biomarkers for Feed Efficiency-Related Traits in Nelore Cattle Identified by Co-expression Network and Integrative Genomics Analyses

Feed efficiency helps to reduce environmental impacts from livestock production, improving beef cattle profitability. We identified potential biomarkers (hub genes) for feed efficiency, by applying co-expression analysis in Longissimus thoracis RNA-Seq data from 180 Nelore steers. Six co-expression modules were associated with six feed efficiency-related traits (p-value ≤ 0.05). Within these modules, 391 hub genes were enriched for pathways as protein synthesis, muscle growth, and immune response. Trait-associated transcription factors (TFs) ELF1, ELK3, ETS1, FLI1, and TCF4, were identified with binding sites in at least one hub gene. Gene expression of CCDC80, FBLN5, SERPINF1, and OGN was associated with multiple feed efficiency-related traits (FDR ≤ 0.05) and were previously related to glucose homeostasis, oxidative stress, fat mass, and osteoblastogenesis, respectively. Potential regulatory elements were identified, integrating the hub genes with previous studies from our research group, such as the putative cis-regulatory elements (eQTLs) inferred as affecting the PCDH18 and SPARCL1 hub genes related to immune system and adipogenesis, respectively. Therefore, our analyses contribute to a better understanding of the biological mechanisms underlying feed efficiency in bovine and the hub genes disclosed can be used as biomarkers for feed efficiency-related traits in Nelore cattle.

Review: Enhancing intramuscular fat development via targeting fibro-adipogenic progenitor cells in meat animals

In the livestock industry, subcutaneous and visceral fat pads are considered as wastes, while intramuscular fat or marbling fat is essential for improving flavor and palatability of meat. Thus, strategies for optimizing fat deposition are needed. Intramuscular adipocytes provide sites for lipid deposition and marbling formation. In the present article, we addressed the origin and markers of intramuscular adipocyte progenitors - fibro-adipogenic progenitors (FAPs), as well as the latest progresses in mechanisms regulating the proliferation and differentiation of intramuscular FAPs. Finally, by targeting intramuscular FAPs, possible nutritional manipulations to improve marbling fat deposition are discussed. Despite recent progresses, the properties and regulation of intramuscular FAPs in livestock remain poorly understood and deserve further investigation.

Integrating genome-wide co-association and gene expression to identify putative regulators and predictors of feed efficiency in pigs

BACKGROUND

Feed efficiency (FE) has a major impact on the economic sustainability of pig production. We used a systems-based approach that integrates single nucleotide polymorphism (SNP) co-association and gene-expression data to identify candidate genes, biological pathways, and potential predictors of FE in a Duroc pig population.

RESULTS

We applied an association weight matrix (AWM) approach to analyse the results from genome-wide association studies (GWAS) for nine FE associated and production traits using 31K SNPs by defining residual feed intake (RFI) as the target phenotype. The resulting co-association network was formed by 829 SNPs. Additive effects of this SNP panel explained 61% of the phenotypic variance of RFI, and the resulting phenotype prediction accuracy estimated by cross-validation was 0.65 (vs. 0.20 using pedigree-based best linear unbiased prediction and 0.12 using the 31K SNPs). Sixty-eight transcription factor (TF) genes were identified in the co-association network; based on the lossless approach, the putative main regulators were COPS5, GTF2H5, RUNX1, HDAC4, ESR1, USP16, SMARCA2 and GTF2F2. Furthermore, gene expression data of the gluteus medius muscle was explored through differential expression and multivariate analyses. A list of candidate genes showing functional and/or structural associations with FE was elaborated based on results from both AWM and gene expression analyses, and included the aforementioned TF genes and other ones that have key roles in metabolism, e.g. ESRRG, RXRG, PPARGC1A, TCF7L2, LHX4, MAML2, NFATC3, NFKBIZ, TCEA1, CDCA7L, LZTFL1 or CBFB. The most enriched biological pathways in this list were associated with behaviour, immunity, nervous system, and neurotransmitters, including melatonin, glutamate receptor, and gustation pathways. Finally, an expression GWAS allowed identifying 269 SNPs associated with the candidate genes' expression (eSNPs). Addition of these eSNPs to the AWM panel of 829 SNPs did not improve the accuracy of genomic predictions.

CONCLUSIONS

Candidate genes that have a direct or indirect effect on FE-related traits belong to various biological processes that are mainly related to immunity, behaviour, energy metabolism, and the nervous system. The pituitary gland, hypothalamus and thyroid axis, and estrogen signalling play fundamental roles in the regulation of FE in pigs. The 829 selected SNPs explained 61% of the phenotypic variance of RFI, which constitutes a promising perspective for applying genetic selection on FE relying on molecular-based prediction.

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

Background

Improving feed efficiency (FE) is a major challenge in pig production. This complex trait is characterized by a high variability. Therefore, the identification of predictors of FE may be a relevant strategy to reduce phenotyping efforts in breeding and selection programs. The aim of this study was to investigate the suitability of expressed muscle genes in prediction of FE traits in growing pigs. The approach considered different transcriptomics experiments to cover a large range of FE values and identify reliable predictors.

Results

Microarrays data were obtained from longissimus muscles of two lines divergently selected for residual feed intake (RFI). Pigs (n = 71) from three experiments belonged to generations 6 to 8 of selection, were fed either a diet with a standard composition or a diet rich in fiber and lipids, received feed ad libitum or at restricted level, and weighed between 80 and 115 kg at slaughter. For each pig, breeding value for RFI was estimated (RFI-BV), and feed conversion ratio (FCR) and energy-based feed conversion ratio (FCRe) were calculated during the test periods. Gradient boosting algorithms were used on the merged muscle transcriptomes to identify very important predictors of FE traits. About 20,405 annotated molecular probes were commonly expressed in longissimus muscle across experiments. Six to 267 expressed muscle genes covering a variety of biological processes were found as important predictors for RFI-BV (R² = 0.63–0.65), FCR (R² = 0.61–0.70) and FCRe (R² = 0.49–0.52). The error of prediction was less than 8% for FCR. Altogether, 56 predictors were common to RFI-BV and FCR. Expression levels of 24 target genes were further measured by qPCR. Linear regression confirmed the good accuracy of combining mRNA levels of these genes to fit FE traits (RFI-BV: R² = 0.73, FRC: R² = 0.76; FCRe: R² = 0.75). Stepwise regression procedure highlighted 10 genes (FKBP5, MUM1, AKAP12, FYN, TMED3, PHKB, TGF, SOCS6, ILR4, and FRAS1) in a linear combination predicting FCR and FCRe. In addition, FKBP5 and expression levels of five other genes (IGF2, SERINC3, CSRNP3, EZR and RPL16) significantly contributed to RFI-BV.

Conclusion

It was possible to identify few genes expressed in muscle that might be reliable predictors of feed efficiency.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6010-9) contains supplementary material, which is available to authorized users.

Transcriptomic Analysis Provides Novel Insights into Heat Stress Responses in Sheep

Simple Summary

The general increase in global temperatures has meant that heat stress has become an increasingly significant problem for sheep. This has both direct and indirect impact on their physiological functions, productivity, and health of sheep. Sheep generally live in high-temperature environments; however, the genes and pathways that play regulatory roles in the heat stress responses of sheep remain unclear. In this study, we applied RNA-Seq technology to analyze liver tissues of sheep from heat-stressed and control groups, and screened genes and pathways related to sheep heat stress. This work provides a theoretical foundation for the breeding and production of heat-resistant sheep.

Abstract

With the intensified and large-scale development of sheep husbandry and global warming, sheep heat stress has become an increasingly important issue. However, little is known about the molecular mechanisms related to sheep responses to heat stress. In this study, transcriptomic analysis of liver tissues of sheep in the presence and absence of heat stress was conducted, with the goal of identifying genes and pathways related to regulation when under such stress. After a comparison with the sheep reference genome, 440,226,436 clean reads were obtained from eight libraries. A p-value ≤ 0.05 and fold change ≥ 2 were taken as thresholds for categorizing differentially expressed genes, of which 1137 were identified. The accuracy and reliability of the RNA-Seq results were confirmed by qRT-PCR. The identified differentially expressed genes were significantly associated with 419 GO terms and 51 KEGG pathways, which suggested their participation in biological processes such as response to stress, immunoreaction, and fat metabolism. This study’s results provide a comprehensive overview of sheep heat stress-induced transcriptional expression patterns, laying a foundation for further analysis of the molecular mechanisms of sheep heat stress.

A comparative analysis of the transcriptome profiles of liver and muscle tissue in pigs divergent for feed efficiency

BACKGROUND

The improvement of feed efficiency is a key economic goal within the pig production industry. The objective of this study was to examine transcriptomic differences in both the liver and muscle of pigs divergent for feed efficiency, thus improving our understanding of the molecular mechanisms influencing feed efficiency and enabling the identification of candidate biomarkers. Residual feed intake (RFI) was calculated for two populations of pigs from two different farms of origin/genotype. The 6 most efficient (LRFI) and 6 least efficient (HRFI) animals from each population were selected for further analysis of Longissimus Dorsi muscle (n = 22) and liver (n = 23). Transcriptomic data were generated from liver and muscle collected post-slaughter.

RESULTS

The transcriptomic data segregated based on the RFI value of the pig rather than genotype/farm of origin. A total of 6463 genes were identified as being differentially expressed (DE) in muscle, while 964 genes were identified as being DE in liver. Genes that were commonly DE between muscle and liver (n = 526) were used for the multi-tissue analysis. These 526 genes were associated with protein targeting to membrane, extracellular matrix organisation and immune function. In the muscle-only analysis, genes associated with RNA processing, protein synthesis and energy metabolism were down regulated in the LRFI animals while in the liver-only analysis, genes associated with cell signalling and lipid homeostasis were up regulated in the LRFI animals.

CONCLUSIONS

Differences in the transcriptome segregated on pig RFI value rather than the genotype/farm of origin. Multi-tissue analysis identified that genes associated with GO terms protein targeting to membrane, extracellular matrix organisation and a range of terms relating to immune function were over represented in the differentially expressed genes of both liver and muscle.

Machine learning applied to transcriptomic data to identify genes associated with feed efficiency in pigs

Background

To date, the molecular mechanisms that underlie residual feed intake (RFI) in pigs are unknown. Results from different genome-wide association studies and gene expression analyses are not always consistent. The aim of this research was to use machine learning to identify genes associated with feed efficiency (FE) using transcriptomic (RNA-Seq) data from pigs that are phenotypically extreme for RFI.

Methods

RFI was computed by considering within-sex regression on mean metabolic body weight, average daily gain, and average backfat gain. RNA-Seq analyses were performed on liver and duodenum tissue from 32 high and 33 low RFI pigs collected at 153 d of age. Machine-learning algorithms were used to predict RFI class based on gene expression levels in liver and duodenum after adjusting for batch effects. Genes were ranked according to their contribution to the classification using the permutation accuracy importance score in an unbiased random forest (RF) algorithm based on conditional inference. Support vector machine, RF, elastic net (ENET) and nearest shrunken centroid algorithms were tested using different subsets of the top rank genes. Nested resampling for hyperparameter tuning was implemented with tenfold cross-validation in the outer and inner loops.

Results

The best classification was obtained with ENET using the expression of 200 genes in liver [area under the receiver operating characteristic curve (AUROC): 0.85; accuracy: 0.78] and 100 genes in duodenum (AUROC: 0.76; accuracy: 0.69). Canonical pathways and candidate genes that were previously reported as associated with FE in several species were identified. The most remarkable pathways and genes identified were NRF2-mediated oxidative stress response and aldosterone signalling in epithelial cells, the DNAJC6, DNAJC1, MAPK8, PRKD3 genes in duodenum, and melatonin degradation II, PPARα/RXRα activation, and GPCR-mediated nutrient sensing in enteroendocrine cells and SMOX, IL4I1, PRKAR2B, CLOCK and CCK genes in liver.

Conclusions

ML algorithms and RNA-Seq expression data were found to provide good performance for classifying pigs into high or low RFI groups. Classification was better with gene expression data from liver than from duodenum. Genes associated with FE in liver and duodenum tissue that can be used as predictive biomarkers for this trait were identified.

Electronic supplementary material

The online version of this article (10.1186/s12711-019-0453-y) contains supplementary material, which is available to authorized users.

RNA-Seq of Liver From Pigs Divergent in Feed Efficiency Highlights Shifts in Macronutrient Metabolism, Hepatic Growth and Immune Response

Liver is a metabolically complex organ that influences nutrient partitioning and potentially modulates the efficiency of converting energy acquired from macronutrients ingestion into a muscle and/or adipose tissue (referred to as feed efficiency, FE). The objective of this study was to sequence the hepatic tissue transcriptome of closely related but differently feed efficient pigs (n = 16) and identify relevant biological processes that underpin the differences in liver phenotype between FE groups. Liver weight did not significantly differ between the FE groups, however, blood parameters showed that total protein, glucose, cholesterol and percentage of lymphocytes were significantly greater in high-FE pigs. Ontology analysis revealed carbohydrate, lipid and protein metabolism to be significantly enriched with differentially expressed genes. In particular, high-FE pigs exhibited gene expression patterns suggesting improved absorption of carbohydrates and cholesterol as well as enhanced reverse cholesterol transport. Furthermore, the inferred decrease in bile acid synthesis in high-FE pigs may contribute to the observed greater levels of serum glucose, which can be then delivered to cells and utilized for growth and maintenance. Gene ontology analysis also suggested that livers of more efficient pigs may be characterized by higher protein turnover and increased epithelial cell differentiation, whereby an enhanced quantity of invariant natural killer T-cells and viability of natural killer cells could induce a quicker and more effective hepatic response to inflammatory stimuli. Our findings suggest that this prompt hepatic response to inflammation in high-FE group may contribute to the more efficient utilization of nutrients for growth in these animals.

Transcriptome analysis of adipose tissue from pigs divergent in feed efficiency reveals alteration in gene networks related to adipose growth, lipid metabolism, extracellular matrix, and immune response

Adipose tissue is hypothesized to play a vital role in regulation of feed efficiency (FE; efficiency in converting energy and nutrients into tissue), of which improvement will simultaneously reduce environmental impact and feed cost per pig. The objective of the present study was to sequence the subcutaneous adipose tissue transcriptome in FE-divergent pigs (n = 16) and identify relevant biological processes underpinning observed differences in FE. We previously demonstrated that high-FE pigs were associated with lower fatness when compared to their counterparts. Here, ontology analysis of a total of 209 annotated genes that were differentially expressed at a p < 0.01 revealed establishment of a dense extracellular matrix and inhibition of capillary formation as one underlying mechanism to achieve suppressed adipogenesis. Moreover, mechanisms ensuring an efficient utilization of lipids in high-FE pigs might be orchestrated by upstream regulators including CEBPA and EGF. Consequently, high-FE adipose tissue could exhibit more efficient cholesterol disposal, whilst inhibition of inflammatory and immune response in high-FE pigs may be an indicator of an optimally functioning adipose tissue. Taken together, adipose tissue growth, extracellular matrix formation, lipid metabolism and inflammatory and immune response are key biological events underpinning the differences in FE. Further investigations focusing on elucidating these processes would assist the animal production industry in optimizing strategies related to nutrient utilization and product quality.