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

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.

Pork carcass composition, meat and belly qualities as influenced by feed efficiency selection in replacement boars from Large White sire and dam lines

This study evaluated carcass attributes, meat and belly qualities in finisher boars (n = 79) selected for feed efficiency (low, intermediate and high) based on estimated breeding value for feed conversion ratio within a Large White dam and sire genetic lines. The sire line had lower trimmed fat proportions and higher lean than the dam line (P < 0.01). Genetic lines expressed slight colour changes and drip losses (P < 0.05), with no differences in pH, marbling and cooking traits (P > 0.05). High-efficient animals presented the highest lean yield (P < 0.01), the lowest trimmed fat proportion (P < 0.01) and no effect on meat and belly quality attributes (P > 0.05) compared with other efficient groups. Interaction between efficiency group and genetic line was only detected for belly weight and thickness (P < 0.01). High-efficient animals offer a greater leanness level, with minimal impact on meat and belly quality traits.

Transcriptomic analysis of intestinal organoids, derived from pigs divergent in feed efficiency, and their response to Escherichia coli

BACKGROUND

There is increasing interest in using intestinal organoids to study complex traits like feed efficiency (FE) and host-microbe interactions. The aim of this study was to investigate differences in the molecular phenotype of organoids derived from pigs divergent for FE as well as their responses to challenge with adherent and invasive Escherichia coli (E. coli).

RESULTS

Colon and ileum tissue from low and high FE pigs was used to generate 3D organoids and two dimensional (2D) monolayers of organoid cells for E. coli challenge. Genome-wide gene expression was used to investigate molecular differences between pigs that were phenotypically divergent for FE and to study the difference in gene expression after challenge with E. coli. We showed, (1) minor differences in gene expression of colon organoids from pigs with low and high FE phenotypes, (2) that an E. coli challenge results in a strong innate immune gene response in both colon and ileum organoids, (3) that the immune response seems to be less pronounced in the colon organoids of high FE pigs and (4) a slightly stronger immune response was observed in ileum than in colon organoids.

CONCLUSIONS

These findings demonstrate the potential for using organoids to gain insights into complex biological mechanisms such as FE.

Small networks of expressed genes in the whole blood and relationships to profiles in circulating metabolites provide insights in inter-individual variability of feed efficiency in growing pigs

BACKGROUND

Feed efficiency is a research priority to support a sustainable meat production. It is recognized as a complex trait that integrates multiple biological pathways orchestrated in and by various tissues. This study aims to determine networks between biological entities to explain inter-individual variation of feed efficiency in growing pigs.

RESULTS

The feed conversion ratio (FCR), a measure of feed efficiency, and its two component traits, average daily gain and average daily feed intake, were obtained from 47 growing pigs from a divergent selection for residual feed intake and fed high-starch or high-fat high-fiber diets during 58 days. Datasets of transcriptomics (60 k porcine microarray) in the whole blood and metabolomics (1H-NMR analysis and target gas chromatography) in plasma were available for all pigs at the end of the trial. A weighted gene co-expression network was built from the transcriptomics dataset, resulting in 33 modules of co-expressed molecular probes. The eigengenes of eight of these modules were significantly ([Formula: see text]) or tended to be ([Formula: see text]) correlated to FCR. Great homogeneity in the enriched biological pathways was observed in these modules, suggesting co-expressed and co-regulated constitutive genes. They were mainly enriched in genes participating to immune and defense-related processes, and to a lesser extent, to translation, cell development or learning. They were also generally associated with growth rate and percentage of lean mass. In the whole network, only one module composed of genes participating to the response to substances, was significantly associated with daily feed intake and body adiposity. The plasma profiles in circulating metabolites and in fatty acids were summarized by weighted linear combinations using a dimensionality reduction method. Close association was thus found between a module composed of co-expressed genes participating to T cell receptor signaling and cell development process in the whole blood and related to FCR, and the circulating concentrations of polyunsaturated fatty acids in plasma.

CONCLUSION

These systemic approaches have highlighted networks of entities driving key biological processes involved in the phenotypic difference in feed efficiency between animals. Connecting transcriptomics and metabolic levels together had some additional benefits.

The mRNA-lncRNA landscape of multiple tissues uncovers key regulators and molecular pathways that underlie heterosis for feed intake and efficiency in laying chickens

BACKGROUND

Heterosis is routinely exploited to improve animal performance. However, heterosis and its underlying molecular mechanism for feed intake and efficiency have been rarely explored in chickens. Feed efficiency continues to be an important breeding goal trait since feed accounts for 60 to 70% of the total production costs in poultry. Here, we profiled the mRNA-lncRNA landscape of 96 samples of the hypothalamus, liver and duodenum mucosa from White Leghorn (WL), Beijing-You chicken (YY), and their reciprocal crosses (WY and YW) to elucidate the regulatory mechanisms of heterosis.

RESULTS

We observed negative heterosis for both feed intake and residual feed intake (RFI) in YW during the laying period from 43 to 46 weeks of age. Analysis of the global expression pattern showed that non-additivity was a major component of the inheritance of gene expression in the three tissues for YW but not for WY. The YW-specific non-additively expressed genes (YWG) and lncRNA (YWL) dominated the total number of non-additively expressed genes and lncRNA in the hypothalamus and duodenum mucosa. Enrichment analysis of YWG showed that mitochondria components and oxidation phosphorylation (OXPHOS) pathways were shared among the three tissues. The OXPHOS pathway was enriched by target genes for YWL with non-additive inheritance of expression in the liver and duodenum mucosa. Weighted gene co-expression network analysis revealed divergent co-expression modules associated with feed intake and RFI in the three tissues from WL, YW, and YY. Among the negatively related modules, the OXPHOS pathway was enriched by hub genes in the three tissues, which supports the critical role of oxidative phosphorylation. Furthermore, protein quantification of ATP5I was highly consistent with ATP5I expression in the liver, which suggests that, in crossbred YW, non-additive gene expression is down-regulated and decreases ATP production through oxidative phosphorylation, resulting in negative heterosis for feed intake and efficiency.

CONCLUSIONS

Our results demonstrate that non-additively expressed genes and lncRNA involved in oxidative phosphorylation in the hypothalamus, liver, and duodenum mucosa are key regulators of the negative heterosis for feed intake and RFI in layer chickens. These findings should facilitate the rational choice of suitable parents for producing crossbred chickens.

Characterization of feed efficiency-related key signatures molecular in different cattle breeds

Feed efficiency is a major constraint in the beef industry and has a significant negative correlation with residual feed intake (RFI). RFI is widely used as a measure of feed efficiency in beef cattle and is independent of economic traits such as body weight and average daily gain. However, key traits with commonality or specificity among beef cattle breeds at the same level of RFI have not been reported. Accordingly, the present study hypothesized that signatures associated with feed efficiency would have commonality or specificity in the liver of cattle breeds at the same RFI level. By comparing and integrating liver transcriptome data, we investigated the critical signatures closely associated with RFI in beef cattle using weighted co-expression network analysis, consensus module analysis, functional enrichment analysis and protein network interaction analysis. The results showed that the consensus modules in Angus and Charolais cattle were negatively correlated, and four (turquoise, red, tan, yellow) were significantly positively correlated in Angus liver, while (turquoise, red) were significantly negatively correlated in Charolais liver. These consensus modules were found to be primarily involved in biological processes such as substance metabolism, energy metabolism and gene transcription, which may be one of the possible explanations for the difference in feed efficiency between the two beef breeds. This research also identified five key candidate genes, PLA2G12B, LCAT, MTTP, LCAT, ABCA1 and FADS1, which are closely associated with hepatic lipid metabolism. The present study has identified some modules, genes and pathways that may be the major contributors to the variation in feed efficiency among different cattle breeds, providing a new perspective on the molecular mechanisms of feed efficiency in beef cattle and a research basis for investigating molecular markers associated with feed efficiency in beef cattle.

Combined transcriptome and metabolome analysis reveal key regulatory genes and pathways of feed conversion efficiency of oriental river prawn Macrobrachium nipponense

BACKGROUND

Oriental river prawn Macrobrachium nipponense is an economically important aquaculture species in China, Japan, and Vietnam. In commercial prawn farming, feed cost constitutes about 50 to 65% of the actual variable cost. Improving feed conversion efficiency in prawn culture will not only increase economic benefit, but also save food and protect the environment. The common indicators used for feed conversion efficiency include feed conversion ratio (FCR), feed efficiency ratio (FER), and residual feed intake (RFI). Among these, RFI is much more suitable than FCR and FER during the genetic improvement of feed conversion efficiency for aquaculture species.

RESULTS

In this study, the transcriptome and metabolome of hepatopancreas and muscle of M. nipponense from high RFI low RFI groups, which identified after culture for 75 days, were characterized using combined transcriptomic and metabolomic analysis. A total of 4540 differentially expressed genes (DEGs) in hepatopancreas, and 3894 DEGs in muscle were identified, respectively. The DEGs in hepatopancreas were mainly enriched in KEGG pathways including the metabolism of xenobiotics by cytochrome P450 (down-regulated), fat digestion and absorption (down-regulated) and aminoacyl-tRNA biosynthesis (up-regulated), etc. The DEGs in muscle were mainly enriched in KEGG pathways including the protein digestion and absorption (down-regulated), glycolysis/gluconeogenesis (down-regulated), and glutathione metabolism (up-regulated), etc. At the transcriptome level, the RFI of M. nipponense was mainly controlled in biological pathways such as the high immune expression and the reduction of nutrients absorption capacity. A total of 445 and 247 differently expressed metabolites (DEMs) were identified in the hepatopancreas and muscle, respectively. At the metabolome level, the RFI of M. nipponense was affected considerably by amino acid and lipid metabolism.

CONCLUSIONS

M. nipponense from higher and lower RFI groups have various physiological and metabolic capability processes. The down-regulated genes, such as carboxypeptidase A1, 6-phosphofructokinase, long-chain-acyl-CoA dehydrogenase, et. al., in digestion and absorption of nutrients, and the up-regulated metabolites, such as aspirin, lysine, et. al., in response to immunity could be potential candidate factors contributed to RFI variation for M. nipponense. Overall, these results would provide new insights into the molecular mechanism of feed conversion efficiency and assist in selective breeding to improve feed conversion efficiency in M. nipponense.

Comparative Serum Proteome Analysis Indicates a Negative Correlation between a Higher Immune Level and Feed Efficiency in Pigs

Identifying and verifying appropriate biomarkers is instrumental in improving the prediction of early-stage pig production performance while reducing the cost of breeding and production. The main factor that affects the production cost and environmental protection cost of the pig industry is the feed efficiency of pigs. This study aimed to detect the differentially expressed proteins in the early blood index determination serum between high-feed efficiency and low-feed efficiency pigs and to provide a basis for further identification of biomarkers using the isobaric tandem mass tag and parallel reaction monitoring approach. In total, 350 (age, 90 ± 2 d; body weight, 41.20 ± 4.60 kg) purebred Yorkshire pigs were included in the study, and their serum samples were obtained during the early blood index determination. The pigs were then arranged based on their feed efficiency; 24 pigs with extreme phenotypes were grouped as high-feed efficiency and low-feed efficiency, with 12 pigs in each group. A total of 1364 proteins were found in the serum, and 137 of them showed differential expression between the groups with high- and low-feed efficiency, with 44 of them being upregulated and 93 being downregulated. PRM (parallel reaction monitoring) was used to verify 10 randomly chosen differentially expressed proteins. The proteins that were differentially expressed were shown to be involved in nine pathways, including the immune system, digestive system, human diseases, metabolism, cellular processing, and genetic information processing, according to the KEGG and GO analyses. Moreover, all of the proteins enriched in the immune system were downregulated in the high-feed efficiency pigs, suggesting that a higher immune level may not be conducive to improving feed efficiency in pigs. This study provides insights into the important feed efficiency proteins and pathways in pigs, promoting the further development of protein biomarkers for predicting and improving porcine feed efficiency.

Feed efficiency in dairy sheep: An insight from the milk transcriptome

Introduction

As higher feed efficiency in dairy ruminants means a higher capability to transform feed nutrients into milk and milk components, differences in feed efficiency are expected to be partly linked to changes in the physiology of the mammary glands. Therefore, this study aimed to determine the biological functions and key regulatory genes associated with feed efficiency in dairy sheep using the milk somatic cell transcriptome.

Material and methods

RNA-Seq data from high (H-FE, n = 8) and low (L-FE, n = 8) feed efficiency ewes were compared through differential expression analysis (DEA) and sparse Partial Least Square-Discriminant analysis (sPLS-DA).

Results

In the DEA, 79 genes were identified as differentially expressed between both conditions, while the sPLS-DA identified 261 predictive genes [variable importance in projection (VIP) &gt; 2] that discriminated H-FE and L-FE sheep.

Discussion

The DEA between sheep with divergent feed efficiency allowed the identification of genes associated with the immune system and stress in L-FE animals. In addition, the sPLS-DA approach revealed the importance of genes involved in cell division (e.g., KIF4A and PRC1) and cellular lipid metabolic process (e.g., LPL, SCD, GPAM, and ACOX3) for the H-FE sheep in the lactating mammary gland transcriptome. A set of discriminant genes, commonly identified by the two statistical approaches, was also detected, including some involved in cell proliferation (e.g., SESN2, KIF20A, or TOP2A) or encoding heat-shock proteins (HSPB1). These results provide novel insights into the biological basis of feed efficiency in dairy sheep, highlighting the informative potential of the mammary gland transcriptome as a target tissue and revealing the usefulness of combining univariate and multivariate analysis approaches to elucidate the molecular mechanisms controlling complex traits.

Transcriptome Profiling of the Liver in Nellore Cattle Phenotypically Divergent for RFI in Two Genetic Groups

The identification and selection of genetically superior animals for residual feed intake (RFI) could enhance productivity and minimize environmental impacts. The aim of this study was to use RNA-seq data to identify the differentially expressed genes (DEGs), known non-coding RNAs (ncRNAs), specific biomarkers and enriched biological processes associated with RFI of the liver in Nellore cattle in two genetic groups. In genetic group 1 (G1), 24 extreme RFI animals (12 low RFI (LRFI) versus 12 high RFI (HRFI)) were selected from a population of 60 Nellore bulls. The RNA-seq of the samples from their liver tissues was performed using an Illumina HiSeq 2000. In genetic group 2 (G2), 20 samples of liver tissue of Nellore bulls divergent for RFI (LRFI, n = 10 versus HRFI, n = 10) were selected from 83 animals. The raw data of the G2 were chosen from the ENA repository. A total of 1811 DEGs were found for the G1 and 2054 for the G2 (p-value ≤ 0.05). We detected 88 common genes in both genetic groups, of which 33 were involved in the immune response and in blocking oxidative stress. In addition, seven (B2M, ADSS, SNX2, TUBA4A, ARHGAP18, MECR, and ABCF3) possible gene biomarkers were identified through a receiver operating characteristic analysis (ROC) considering an AUC > 0.70. The B2M gene was overexpressed in the LRFI group. This gene regulates the lipid metabolism protein turnover and inhibits cell death. We also found non-coding RNAs in both groups. MIR25 was up-regulated and SNORD16 was down-regulated in the LRFI for G1. For G2, up-regulated RNase_MRP and SCARNA10 were found. We highlight MIR25 as being able to act by blocking cytotoxicity and oxidative stress and RMRP as a blocker of mitochondrial damage. The biological pathways associated with RFI of the liver in Nellore cattle in the two genetic groups were for energy metabolism, protein turnover, redox homeostasis and the immune response. The common transcripts, biomarkers and metabolic pathways found in the two genetic groups make this unprecedented work even more relevant, since the results are valid for different herds raised in different ways. The results reinforce the biological importance of these known processes but also reveal new insights into the complexity of the liver tissue transcriptome of Nellore cattle.

Multi-transcriptomics reveals RLMF axis-mediated signaling molecules associated with bovine feed efficiency

The regulatory axis plays a vital role in interpreting the information exchange and interactions among mammal organs. In this study on feed efficiency, it was hypothesized that a rumen-liver-muscle-fat (RLMF) regulatory axis exists and scrutinized the flow of energy along the RLMF axis employing consensus network analysis from a spatial transcriptomic standpoint. Based on enrichment analysis and protein-protein interaction analysis of the consensus network and tissue-specific genes, it was discovered that carbohydrate metabolism, energy metabolism, immune and inflammatory responses were likely to be the biological processes that contribute most to feed efficiency variation on the RLMF regulatory axis. In addition, clusters of genes related to the electron respiratory chain, including ND (2,3,4,4L,5,6), NDUF (A13, A7, S6, B3, B6), COX (1,3), CYTB, UQCR11, ATP (6,8), clusters of genes related to fatty acid metabolism including APO (A1, A2, A4, B, C3), ALB, FG (A, G), as well as clusters of the ribosomal-related gene including RPL (8,18A,18,15,13, P1), the RPS (23,27A,3A,4X), and the PSM (A1-A7, B6, C1, C3, D2-D4, D8 D9, E1) could be the primary effector genes responsible for feed efficiency variation. The findings demonstrate that high feed efficiency cattle, through the synergistic action of the regulatory axis RLMF, may improve the efficiency of biological processes (carbohydrate metabolism, protein ubiquitination, and energy metabolism). Meanwhile, high feed efficiency cattle might enhance the ability to respond to immunity and inflammation, allowing nutrients to be efficiently distributed across these organs associated with digestion and absorption, energy-producing, and energy-storing organs. Elucidating the distribution of nutrients on the RLMF regulatory axis could facilitate an understanding of feed efficiency variation and achieve the study on its molecular regulation.

Identification of transcriptional regulatory variants in pig duodenum, liver, and muscle tissues

BACKGROUND

In humans and livestock species, genome-wide association studies (GWAS) have been applied to study the association between variants distributed across the genome and a phenotype of interest. To discover genetic polymorphisms affecting the duodenum, liver, and muscle transcriptomes of 300 pigs from 3 different breeds (Duroc, Landrace, and Large White), we performed expression GWAS between 25,315,878 polymorphisms and the expression of 13,891 genes in duodenum, 12,748 genes in liver, and 11,617 genes in muscle.

RESULTS

More than 9.68 × 1011 association tests were performed, yielding 14,096,080 significantly associated variants, which were grouped in 26,414 expression quantitative trait locus (eQTL) regions. Over 56% of the variants were within 1 Mb of their associated gene. In addition to the 100-kb region upstream of the transcription start site, we identified the importance of the 100-kb region downstream of the 3'UTR for gene regulation, as most of the cis-regulatory variants were located within these 2 regions. We also observed 39,874 hotspot regulatory polymorphisms associated with the expression of 10 or more genes that could modify the protein structure or the expression of a regulator gene. In addition, 2 motifs (5'-GATCCNGYGTTGCYG-3' and a poly(A) sequence) were enriched across the 3 tissues within the neighboring sequences of the most significant single-nucleotide polymorphisms in each cis-eQTL region.

CONCLUSIONS

The 14 million significant associations obtained in this study are publicly available and have enabled the identification of expression-associated cis-, trans-, and hotspot regulatory variants within and across tissues, thus shedding light on the molecular mechanisms of regulatory variations that shape end-trait phenotypes.

Identification of Differentially Expressed miRNAs in Porcine Adipose Tissues and Evaluation of Their Effects on Feed Efficiency

Feed efficiency (FE) is a very important trait affecting the economic benefits of pig breeding enterprises. Adipose tissue can modulate a variety of processes such as feed intake, energy metabolism and systemic physiological processes. However, the mechanism by which microRNAs (miRNAs) in adipose tissues regulate FE remains largely unknown. Therefore, this study aimed to screen potential miRNAs related to FE through miRNA sequencing. The miRNA profiles in porcine adipose tissues were obtained and 14 miRNAs were identified differentially expressed in adipose tissues of pigs with extreme differences in FE, of which 9 were down-regulated and 5 were up-regulated. GO and KEGG analyses indicated that these miRNAs were significantly related to lipid metabolism and these miRNAs modulated FE by regulating lipid metabolism. Subsequently, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) of five randomly selected DEMs was used to verify the reliability of miRNA-seq data. Furthermore, 39 differentially expressed target genes of these DEMs were obtained, and DEMs-target mRNA interaction networks were constructed. In addition, the most significantly down-regulated miRNAs, ssc-miR-122-5p and ssc-miR-192, might be the key miRNAs for FE. Our results reveal the mechanism by which adipose miRNAs regulate feed efficiency in pigs. This study provides a theoretical basis for the further study of swine feed efficiency improvement.

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.

Transcription Analysis of Liver and Muscle Tissues from Landrace Finishing Pigs with Different Feed Conversion Ratios

The efficiency of feed utilization determines the cost and economic benefits of pig production. In the present study, two pairs of full-sibling and two pairs of half-sibling female Landrace finishing pigs were selected, with each pair including individuals with different feed conversion rates, with liver and longissimus muscle tissue samples collected from each group for transcriptome analysis. A total of 561 differentially expressed genes (DEGs), among which 224 were up-regulated and 337 were down-regulated, were detected in the liver transcriptomes in the high-feed efficiency group compared to the low-feed efficiency group. The DEGs related to phosphorus and phosphate metabolism, arginine biosynthesis, chemical carcinogenesis, cytokine-cytokine receptor interaction, the biosynthesis of amino acids, and drug metabolism-cytochrome P450 in liver tissue were also associated with feed efficiency. In total, 215 DEGs were screened in the longissimus muscle tissue and were mainly related to disease and immune regulation, including complement and coagulation cascades, systemic lupus erythematosus, and prion diseases. The combination of gene expression and functional annotation results led to the identification of candidate feed efficiency-related biomarkers, such as ARG1, ARG2, GOT1, GPT2, ACAA2, ACADM, and ANGPTL4, members of cytochrome P450 family, and complement component family genes. Although the novel feed efficiency-related candidate genes need to be further evaluated by a larger sample size and functional studies, the present study identifies novel candidate biomarkers for the identification of functional SNPs underlying porcine feed efficiency.

Gut transcriptome reveals differential gene expression and enriched pathways linked to immune activation in response to weaning in pigs

Weaning represents one of the most critical periods in pig production associated with increase in disease risk, reduction in performance and economic loss. Physiological changes faced by piglets during the weaning period have been well characterised, however little is currently known about the underlying molecular pathways involved in these processes. As pig meat remains one of the most consumed sources of protein worldwide, understanding how these changes are mediated is critical to improve pig production and consequently sustainable food production globally. In this study, we evaluated the effect of weaning on transcriptomic changes in the colon of healthy piglets over time using an RNA-sequencing approach. The findings revealed a complex and coordinated response to weaning with the majority of genes found to be rapidly differentially expressed within 1 day post weaning. Multiple genes and pathways affected by weaning in the colon were associated with immune regulation, cell signalling and bacterial defence. NOD-like receptors, Toll-like receptor and JAK-STAT signalling pathways were amongst the pathways significantly enriched. Immune activation was evidenced by the enrichment of pathways involved in interferon response, cytokines interactions, oxidoreductase activities and response to microbial invasion. Biosynthesis of amino acids, in particular arginine, was also amongst the most enriched KEGG pathways in weaned pigs, reinforcing the critical role of arginine in gut homeostasis under stress conditions. Overall, transcriptomic and physiological results suggest that pigs going through the weaning transition undergo a transient period of inflammatory state with a temporary breakdown of barrier functions in the gut. These findings could provide valuable tools to monitor host response post weaning, and may be of particular relevance for the investigation and development of intervention strategies aimed to reduce antibiotic use and improve pig health and performance.

Genome-Wide Association Analysis and Genetic Parameters for Feed Efficiency and Related Traits in Yorkshire and Duroc Pigs

Simple Summary

Genetic improvements in feed efficiency (FE) and related traits could considerably reduce pig production costs and energy consumption. Thus, we performed a genetic parameter estimation and genome-wide association study of four FE and FE-related traits, namely, average daily feed intake, average daily gain, the feed conversion ratio, and residual feed intake, of two pig breeds, Yorkshire and Duroc. The results demonstrate the genetic relationships of FE and FE-related traits with two growth traits, age and backfat thickness at 100 kg. We also identified many single-nucleotide polymorphisms (SNPs) and novel candidate genes related to these traits. In addition, we found many pathways significantly associated with FE and FE-related traits, and they are generally involved in digestive and metabolic processes. The results of this study are expected to provide a valuable reference for the genomic selection of FE and FE-related traits in pigs.

Abstract

Feed efficiency (FE) traits are key factors that can influence the economic benefits of pig production. However, little is known about the genetic architecture of FE and FE-related traits. This study aimed to identify SNPs and candidate genes associated with FE and FE-related traits, namely, average daily feed intake (ADFI), average daily gain (ADG), the feed conversion ratio (FCR), and residual feed intake (RFI). The phenotypes of 5823 boars with genotyped data (50 K BeadChip) from 1365 boars from a nucleus farm were used to perform a genome-wide association study (GWAS) of two breeds, Duroc and Yorkshire. Moreover, we performed a genetic parameter estimation for four FE and FE-related traits. The heritabilities of the FE and FE-related traits ranged from 0.13 to 0.36, and there were significant genetic correlations (−0.69 to 0.52) of the FE and FE-related traits with two growth traits (age at 100 kg and backfat thickness at 100 kg). A total of 61 significant SNPs located on eight different chromosomes associated with the four FE and FE-related traits were identified. We further identified four regions associated with FE and FE-related traits that have not been previously reported, and they may be potential novel QTLs for FE. Considering their biological functions, we finally identified 35 candidate genes relevant for FE and FE-related traits, such as the widely reported MC4R and INSR genes. A gene enrichment analysis showed that FE and FE-related traits were highly enriched in the biosynthesis, digestion, and metabolism of biomolecules. This study deepens our understanding of the genetic mechanisms of FE in pigs and provides valuable information for using marker-assisted selection in pigs to improve FE.

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.

Protein metabolism in the liver and white muscle is associated with feed efficiency in Chinook salmon (Oncorhynchus tshawytscha) reared in seawater: Evidence from proteomic analysis

Understanding the molecular mechanisms that underlie differences in feed efficiency (FE) is an important step toward optimising growth and achieving sustainable salmonid aquaculture. In this study, the liver and white muscle proteomes of feed efficient (EFF) and inefficient (INEFF) Chinook salmon (Oncorhynchus tshawytscha) reared in seawater were investigated by liquid chromatography-tandem mass spectrometry (LC-MS/MS). In total, 2746 liver and 702 white muscle proteins were quantified and compared between 21 EFF and 22 INEFF fish. GSEA showed that gene sets related to protein synthesis were enriched in the liver and white muscle of the EFF group, while conversely, pathways related to protein degradation (amino acid catabolism and proteolysis, respectively) were the most affected processes in the liver and white muscle of INEFF fish. Estimates of individual daily feed intake and share of the meal within tank were significantly higher in the INEFF than the EFF fish showing INEFF fish were likely more dominant during feeding and overfed. Overeating by the INEFF fish was associated with an increase in protein catabolism. This study found that fish with different FE values had expression differences in the gene sets related to protein turnover, and this result supports the hypothesis that protein metabolism plays a role in FE.

Genetic and genomic analysis of oxygen consumption in mice

We estimated genetic parameters for oxygen consumption (OC), OC per metabolic body weight (OCMBW) and body weight at three through 8 weeks of age in divergently selected mice populations, with an animal model considering maternal genetic, common litter environmental and cytoplasmic inheritance effects. Cytoplasmic inheritance was considered based on maternal lineage information. With respect to OC, estimated direct heritability was moderate (0.32) and the estimated proportion of the variance of cytoplasmic inheritance effects to the phenotypic variance was very low (0.01), implying that causal genes for OC could be located on autosomes. To assess this hypothesis, we attempted to identify possible candidate causal genes through selective signature detection with the results of pooled whole-genome resequencing using pooled DNA samples from high and low OC mice. We made a list of possible candidate causal genes for OC, including those relating to electron transport chain and ATP-binding proteins (Ndufa12, Sdhc, Atp10b, etc.), Prr16 encoding Largen protein, Cry1 encoding a key component of the circadian core oscillator and so on. The results, although careful interpretation must be required, could contribute to elucidate the genetic mechanism of OC, an indicator for maintenance energy requirement, and therefore feed efficiency.

Transcriptome analysis reveals the potential roles of long non-coding RNAs in feed efficiency of chicken

Feed efficiency is an important economic trait and reduces the production costs per unit of animal product. Up to now, few studies have conducted transcriptome profiling of liver tissue in feed efficiency-divergent chickens (Ross vs native breeds). Also, molecular mechanisms contributing to differences in feed efficiency are not fully understood, especially in terms of long non-coding RNAs (lncRNAs). Hence, transcriptome profiles of liver tissue in commercial and native chicken breeds were analyzed. RNA-Seq data along with bioinformatics approaches were applied and a series of lncRNAs and target genes were identified. Furthermore, protein-protein interaction network construction, co-expression analysis, co-localization analysis of QTLs and functional enrichment analysis were used to functionally annotate the identified lncRNAs. In total, 2,290 lncRNAs were found (including 1,110 annotated, 593 known and 587 novel), of which 53 (including 39 known and 14 novel), were identified as differentially expressed genes between two breeds. The expression profile of lncRNAs was validated by RT-qPCR. The identified novel lncRNAs showed a number of characteristics similar to those of known lncRNAs. Target prediction analysis showed that these lncRNAs have the potential to act in cis or trans mode. Functional enrichment analysis of the predicted target genes revealed that they might affect the differences in feed efficiency of chicken by modulating genes associated with lipid metabolism, carbohydrate metabolism, growth, energy homeostasis and glucose metabolism. Some gene members of significant modules in the constructed co-expression networks were reported as important genes related to feed efficiency. Co-localization analysis of QTLs related to feed efficiency and the identified lncRNAs suggested several candidates to be involved in residual feed intake. The findings of this study provided valuable resources to further clarify the genetic basis of regulation of feed efficiency in chicken from the perspective of lncRNAs.

Emerging Roles of Non-Coding RNAs in the Feed Efficiency of Livestock Species

A global population of already more than seven billion people has led to an increased demand for food and water, and especially the demand for meat. Moreover, the cost of feed used in animal production has also increased dramatically, which requires animal breeders to find alternatives to reduce feed consumption. Understanding the biology underlying feed efficiency (FE) allows for a better selection of feed-efficient animals. Non-coding RNAs (ncRNAs), especially micro RNAs (miRNAs) and long non-coding RNAs (lncRNAs), play important roles in the regulation of bio-logical processes and disease development. The functions of ncRNAs in the biology of FE have emerged as they participate in the regulation of many genes and pathways related to the major FE indicators, such as residual feed intake and feed conversion ratio. This review provides the state of the art studies related to the ncRNAs associated with FE in livestock species. The contribution of ncRNAs to FE in the liver, muscle, and adipose tissues were summarized. The research gap of the function of ncRNAs in key processes for improved FE, such as the nutrition, heat stress, and gut–brain axis, was examined. Finally, the potential uses of ncRNAs for the improvement of FE were discussed.

Characterization and Duodenal Transcriptome Analysis of Chinese Beef Cattle With Divergent Feed Efficiency Using RNA-Seq

Residual feed intake (RFI) is an important measure of feed efficiency for agricultural animals. Factors associated with cattle RFI include physiology, dietary factors, and the environment. However, a precise genetic mechanism underlying cattle RFI variations in duodenal tissue is currently unavailable. The present study aimed to identify the key genes and functional pathways contributing to variance in cattle RFI phenotypes using RNA sequencing (RNA-seq). Six bulls with extremely high or low RFIs were selected for detecting differentially expressed genes (DEGs) by RNA-seq, followed by conducting GO, KEGG enrichment, protein-protein interaction (PPI), and co-expression network (WGCNA, n = 10) analysis. A total of 380 differentially expressed genes was obtained from high and low RFI groups, including genes related to energy metabolism (ALDOA, HADHB, INPPL1), mitochondrial function (NDUFS1, RFN4, CUL1), and feed intake behavior (CCK). Two key sub-networks and 26 key genes were detected using GO analysis of DEGs and PPI analysis, such as TPM1 and TPM2, which are involved in mitochondrial pathways and protein synthesis. Through WGCNA, a gene network was built, and genes were sorted into 27 modules, among which the blue (r = 0.72, p = 0.03) and salmon modules (r = -0.87, p = 0.002) were most closely related with RFI. DEGs and genes from the main sub-networks and closely related modules were largely involved in metabolism; oxidative phosphorylation; glucagon, ribosome, and N-glycan biosynthesis, and the MAPK and PI3K-Akt signaling pathways. Through WGCNA, five key genes, including FN1 and TPM2, associated with the biological regulation of oxidative processes and skeletal muscle development were identified. Taken together, our data suggest that the duodenum has specific biological functions in regulating feed intake. Our findings provide broad-scale perspectives for identifying potential pathways and key genes involved in the regulation of feed efficiency in beef cattle.

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

Background

Improving feed efficiency (FE) is an important goal due to its economic and environmental significance for farm animal production. The FE phenotype is complex and based on the measurements of the individual feed consumption and average daily gain during a test period, which is costly and time-consuming. The identification of reliable predictors of FE is a strategy to reduce phenotyping efforts.

Results

Gene expression data of the whole blood from three independent experiments were combined and analyzed by machine learning algorithms to propose molecular biomarkers of FE traits in growing pigs. These datasets included Large White pigs from two lines divergently selected for residual feed intake (RFI), a measure of net FE, and in which individual feed conversion ratio (FCR) and blood microarray data were available. Merging the three datasets allowed considering FCR values (Mean = 2.85; Min = 1.92; Max = 5.00) for a total of n = 148 pigs, with a large range of body weight (15 to 115 kg) and different test period duration (2 to 9 weeks). Random forest (RF) and gradient tree boosting (GTB) were applied on the whole blood transcripts (26,687 annotated molecular probes) to identify the most important variables for binary classification on RFI groups and a quantitative prediction of FCR, respectively. The dataset was split into learning (n = 74) and validation sets (n = 74). With iterative steps for variable selection, about three hundred’s (328 to 391) molecular probes participating in various biological pathways, were identified as important predictors of RFI or FCR. With the GTB algorithm, simpler models were proposed combining 34 expressed unique genes to classify pigs into RFI groups (100% of success), and 25 expressed unique genes to predict FCR values (R² = 0.80, RMSE = 8%). The accuracy performance of RF models was slightly lower in classification and markedly lower in regression.

Conclusion

From small subsets of genes expressed in the whole blood, it is possible to predict the binary class and the individual value of feed efficiency. These predictive models offer good perspectives to identify animals with higher feed efficiency in precision farming applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07843-4.

RNA-Seq Data for Reliable SNP Detection and Genotype Calling: Interest for Coding Variant Characterization and Cis-Regulation Analysis by Allele-Specific Expression in Livestock Species

In addition to their common usages to study gene expression, RNA-seq data accumulated over the last 10 years are a yet-unexploited resource of SNPs in numerous individuals from different populations. SNP detection by RNA-seq is particularly interesting for livestock species since whole genome sequencing is expensive and exome sequencing tools are unavailable. These SNPs detected in expressed regions can be used to characterize variants affecting protein functions, and to study cis-regulated genes by analyzing allele-specific expression (ASE) in the tissue of interest. However, gene expression can be highly variable, and filters for SNP detection using the popular GATK toolkit are not yet standardized, making SNP detection and genotype calling by RNA-seq a challenging endeavor. We compared SNP calling results using GATK suggested filters, on two chicken populations for which both RNA-seq and DNA-seq data were available for the same samples of the same tissue. We showed, in expressed regions, a RNA-seq precision of 91% (SNPs detected by RNA-seq and shared by DNA-seq) and we characterized the remaining 9% of SNPs. We then studied the genotype (GT) obtained by RNA-seq and the impact of two factors (GT call-rate and read number per GT) on the concordance of GT with DNA-seq; we proposed thresholds for them leading to a 95% concordance. Applying these thresholds to 767 multi-tissue RNA-seq of 382 birds of 11 chicken populations, we found 9.5 M SNPs in total, of which ∼550,000 SNPs per tissue and population with a reliable GT (call rate ≥ 50%) and among them, ∼340,000 with a MAF ≥ 10%. We showed that such RNA-seq data from one tissue can be used to (i) detect SNPs with a strong predicted impact on proteins, despite their scarcity in each population (16,307 SIFT deleterious missenses and 590 stop-gained), (ii) study, on a large scale, cis-regulations of gene expression, with ∼81% of protein-coding and 68% of long non-coding genes (TPM ≥ 1) that can be analyzed for ASE, and with ∼29% of them that were cis-regulated, and (iii) analyze population genetic using such SNPs located in expressed regions. This work shows that RNA-seq data can be used with good confidence to detect SNPs and associated GT within various populations and used them for different analyses as GTEx studies.

Identification of the relationship between the gut microbiome and feed efficiency in a commercial pig cohort

Feed efficiency (FE) is an economically important trait in pig production. Gut microbiota plays an important role in energy harvest, nutrient metabolism, and fermentation of dietary indigestible components. Whether and which gut microbes affect FE in pigs are largely unknown. Here, a total of 208 healthy Duroc pigs were used as experimental materials. Feces and serum samples were collected at the age of 140 d. We first performed 16S rRNA gene and metagenomic sequencing analysis to investigate the relationship between the gut microbiome and porcine residual feed intake (RFI). 16S rRNA gene sequencing analysis detected 21 operational taxonomic units showing the tendency to correlation with the RFI (P < 0.01). Metagenomic sequencing further identified that the members of Clostridiales, e.g., Ruminococcus flavefaoiens, Lachnospiraceae bacterium 28-4, and Lachnospiraceae phytofermentans, were enriched in pigs with low RFI (high-FE), while 11 bacterial species including 5 Prevotella spp., especially, the Prevotella copri, had higher abundance in pigs with high RFI. Functional capacity analysis suggested that the gut microbiome of low RFI pigs had a high abundance of the pathways related to amino acid metabolism and biosynthesis, but a low abundance of the pathways associated with monosaccharide metabolism and lipopolysaccharide biosynthesis. Serum metabolome and fecal short-chain fatty acids were determined by UPLC-QTOF/MS and gas chromatography, respectively. Propionic acid in feces and the serum metabolites related to amino acid metabolism were negatively correlated with the RFI. The results from this study may provide potential gut microbial biomarkers that could be used for improving FE in pig production industry.

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.

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.

Transcriptome Profiling across Five Tissues of Giant Panda

Gene differential expression studies can serve to explore and understand the laws and characteristics of animal life activities, and the difference in gene expression between different animal tissues has been well demonstrated and studied. However, for the world-famous rare and protected species giant panda (Ailuropoda melanoleuca), only the transcriptome of the blood and spleen has been reported separately. Here, in order to explore the transcriptome differences between the different tissues of the giant panda, transcriptome profiles of the heart, liver, spleen, lung, and kidney from five captive giant pandas were constructed with Illumina HiSeq 2500 platform. The comparative analysis of the intertissue gene expression patterns was carried out based on the generated RNA sequencing datasets. Analyses of Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and protein-protein interaction (PPI) network were performed according to the identified differentially expressed genes (DEGs). We generated 194.52 GB clean base data from twenty-five sequencing libraries and identified 18,701 genes, including 3492 novel genes. With corrected p value <0.05 and |log₂FoldChange| >2, we finally obtained 921, 553, 574, 457, and 638 tissue-specific DEGs in the heart, liver, spleen, lung, and kidney, respectively. In addition, we identified TTN, CAV3, LDB3, TRDN, and ACTN2 in the heart; FGA, AHSG, and SERPINC1 in the liver; CD19, CD79B, and IL21R in the spleen; NKX2-4 and SFTPB in the lung; GC and HRG in the kidney as hub genes in the PPI network. The results of the analyses showed a similar gene expression pattern between the spleen and lung. This study provided for the first time the heart, liver, lung, and kidney's transcriptome resources of the giant panda, and it provided a valuable resource for further genetic research or other potential research.

Exploring the Regulatory Potential of Long Non-Coding RNA in Feed Efficiency of Indicine Cattle

Long non-coding RNA (lncRNA) can regulate several aspects of gene expression, being associated with complex phenotypes in humans and livestock species. In taurine beef cattle, recent evidence points to the involvement of lncRNA in feed efficiency (FE), a proxy for increased productivity and sustainability. Here, we hypothesized specific regulatory roles of lncRNA in FE of indicine cattle. Using RNA-Seq data from the liver, muscle, hypothalamus, pituitary gland and adrenal gland from Nellore bulls with divergent FE, we submitted new transcripts to a series of filters to confidently predict lncRNA. Then, we identified lncRNA that were differentially expressed (DE) and/or key regulators of FE. Finally, we explored lncRNA genomic location and interactions with miRNA and mRNA to infer potential function. We were able to identify 126 relevant lncRNA for FE in Bos indicus, some with high homology to previously identified lncRNA in Bos taurus and some possible specific regulators of FE in indicine cattle. Moreover, lncRNA identified here were linked to previously described mechanisms related to FE in hypothalamus-pituitary-adrenal axis and are expected to help elucidate this complex phenotype. This study contributes to expanding the catalogue of lncRNA, particularly in indicine cattle, and identifies candidates for further studies in animal selection and management.

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.

A Transcriptome Analysis Reveals that Hepatic Glycolysis and Lipid Synthesis Are Negatively Associated with Feed Efficiency in DLY Pigs

Feed efficiency (FE) is an important trait in the porcine industry. Therefore, understanding the molecular mechanisms of FE is vital for the improvement of this trait. In this study, 6 extreme high-FE and 6 low-FE pigs were selected from 225 Duroc × (Landrace × Yorkshire) (DLY) pigs for transcriptomic analysis. RNA-seq analysis was performed to determine differentially expressed genes (DEGs) in the liver tissues of the 12 individuals, and 507 DEGs were identified between high-FE pigs (HE- group) and low-FE pigs (LE- group). A gene ontology (GO) enrichment and pathway enrichment analysis were performed and revealed that glycolytic metabolism and lipid synthesis-related pathways were significantly enriched within DEGs; all of these DEGs were downregulated in the HE- group. Moreover, Weighted gene co-expression analysis (WGCNA) revealed that oxidative phosphorylation, thermogenesis, and energy metabolism-related pathways were negatively related to HE- group, which might result in lower energy consumption in higher efficiency pigs. These results implied that the higher FE in the HE- group may be attributed to a lower glycolytic, energy consumption and lipid synthesizing potential in the liver. Furthermore, our findings suggested that the inhibition of lipid synthesis and glucose metabolic activity in the liver may be strategies for improving the FE of DLY pigs.

Landscape of multi-tissue global gene expression reveals the regulatory signatures of feed efficiency in beef cattle

MOTIVATION

Feed efficiency is an important trait for sustainable beef production that is regulated by the complex biological process, but the mode of action behinds it has not been clearly defined. Here, we aimed to elucidate the regulatory mechanisms of this trait through studying the landscape of the genome-wide gene expression of rumen, liver, muscle and backfat tissues, the key ones involved in the energy metabolism.

RESULTS

The transcriptome of 189 samples across four tissues from 48 beef steers with varied feed efficiency were generated using Illumina HiSeq4000. The analysis of global gene expression profiles of four tissues, functional analysis of tissue-shared and -unique genes, co-expressed network construction of tissue-shared genes, weighted correlations analysis between gene modules and feed efficiency-related traits in each tissue were performed. Among four tissues, the transcriptome of muscle tissue was distinctive from others, while those of rumen and backfat tissues were similar. The associations between co-expressed genes and feed efficiency related traits at single or all tissues level exhibited that the gene expression in the rumen, liver, muscle and backfat were the most correlated with feed conversion ratio, dry matter intake, average daily gain and residual feed intake, respectively. The 19 overlapped genes identified from the strongest module-trait relationships in four tissues are potential generic gene markers for feed efficiency.

AVAILABILITY AND IMPLEMENTATION

The distribution of gene expression data can be accessed at https://www.cattleomics.com/transcriptome.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome Wide Assessment of Genetic Variation and Population Distinctiveness of the Pig Family in South Africa

Genetic diversity is of great importance and a prerequisite for genetic improvement and conservation programs in pigs and other livestock populations. The present study provides a genome wide analysis of the genetic variability and population structure of pig populations from different production systems in South Africa relative to global populations. A total of 234 pigs sampled in South Africa and consisting of village (n = 91), commercial (n = 60), indigenous (n = 40), Asian (n = 5) and wild (n = 38) populations were genotyped using Porcine SNP60K BeadChip. In addition, 389 genotypes representing village and commercial pigs from America, Europe, and Asia were accessed from a previous study and used to compare population clustering and relationships of South African pigs with global populations. Moderate heterozygosity levels, ranging from 0.204 for Warthogs to 0.371 for village pigs sampled from Capricorn municipality in Eastern Cape province of South Africa were observed. Principal Component Analysis of the South African pigs resulted in four distinct clusters of (i) Duroc; (ii) Vietnamese; (iii) Bush pig and Warthog and (iv) a cluster with the rest of the commercial (SA Large White and Landrace), village, Wild Boar and indigenous breeds of Koelbroek and Windsnyer. The clustering demonstrated alignment with genetic similarities, geographic location and production systems. The PCA with the global populations also resulted in four clusters that where populated with (i) all the village populations, wild boars, SA indigenous and the large white and landraces; (ii) Durocs (iii) Chinese and Vietnamese pigs and (iv) Warthog and Bush pig. K = 10 (The number of population units) was the most probable ADMIXTURE based clustering, which grouped animals according to their populations with the exception of the village pigs that showed presence of admixture. AMOVA reported 19.92%-98.62% of the genetic variation to be within populations. Sub structuring was observed between South African commercial populations as well as between Indigenous and commercial breeds. Population pairwise F ST analysis showed genetic differentiation (P ≤ 0.05) between the village, commercial and wild populations. A per marker per population pairwise F ST analysis revealed SNPs associated with QTLs for traits such as meat quality, cytoskeletal and muscle development, glucose metabolism processes and growth factors between both domestic populations as well as between wild and domestic breeds. Overall, the study provided a baseline understanding of porcine diversity and an important foundation for porcine genomics of South African populations.

Identification and Annotation of Potential Function of Regulatory Antisense Long Non-Coding RNAs Related to Feed Efficiency in Bos taurus Bulls

Long non-coding RNAs (lncRNAs) can influence transcriptional and translational processes in mammalian cells and are associated with various developmental, physiological and phenotypic conditions. However, they remain poorly understood and annotated in livestock species. We combined phenotypic, metabolomics and liver transcriptomic data of bulls divergent for residual feed intake (RFI) and fat accretion. Based on a project-specific transcriptome annotation for the bovine reference genome ARS-UCD.1.2 and multiple-tissue total RNA sequencing data, we predicted 3590 loci to be lncRNAs. To identify lncRNAs with potential regulatory influence on phenotype and gene expression, we applied the regulatory impact factor algorithm on a functionally prioritized set of loci (n = 4666). Applying the algorithm of partial correlation and information theory, significant and independent pairwise correlations were calculated and co-expression networks were established, including plasma metabolites correlated with lncRNAs. The network hub lncRNAs were assessed for potential cis-actions and subjected to biological pathway enrichment analyses. Our results reveal a prevalence of antisense lncRNAs positively correlated with adjacent protein-coding genes and suggest their participation in mitochondrial function, acute phase response signalling, TCA-cycle, fatty acid β-oxidation and presumably gluconeogenesis. These antisense lncRNAs indicate a stabilizing function for their cis-correlated genes and a putative regulatory role in gene expression.

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.

Genes associated with body weight gain and feed intake identified by meta-analysis of the mesenteric fat from crossbred beef steers

Mesenteric fat is a visceral fat depot that increases with cattle maturity and can be influenced by diet. There may be a relationship between the accumulation of mesenteric fat and feed efficiency in beef cattle. The purpose of this study was to identify genes that may be differentially expressed in steers with high and low BW gain and feed intake. RNA-Seq was used to evaluate the transcript abundance of genes in the mesenteric fat from a total of 78 steers collected over 5 different cohorts. A meta-analysis was used to identify genes involved with gain, feed intake or the interaction of both phenotypes. The interaction analysis identified 11 genes as differentially expressed. For the main effect of gain, a total of 87 differentially expressed genes (DEG) were identified (P_ADJ<0.05), and 24 were identified in the analysis for feed intake. Genes identified for gain were involved in functions and pathways including lipid metabolism, stress response/protein folding, cell proliferation/growth, axon guidance and inflammation. The genes for feed intake did not cluster into pathways, but some of the DEG for intake had functions related to inflammation, immunity, and/or signal transduction (JCHAIN, RIPK1, LY86, SPP1, LYZ, CD5, CD53, SRPX, and NF2). At P_ADJ<0.1, only 4 genes (OLFML3, LOC100300716, MRPL15, and PUS10) were identified as differentially expressed in two or more cohorts, highlighting the importance of evaluating the transcriptome of more than one group of animals and incorporating a meta-analysis. This meta-analysis has produced many mesenteric fat DEG that may be contributing to gain and feed intake in cattle.

Identification of genetic loci associated with growth traits at weaning in yak through a genome-wide association study

A multilocus GWAS was performed to explore the genetic architecture of four growth traits in yak. In total, 354 female yaks for which measurements of body weight (BW), withers height (WH), body length (BL) and chest girth (CG) at weaning were available underwent genotyping with the Illumina BovineHD BeadChip (770K). After quality control, we retained 98 688 SNPs and 354 animals for GWAS analysis. We identified seven, 18, seven and nine SNPs (corresponding to seven, 17, seven and eight candidate genes) associated with BW, WH, BL and CG at weaning respectively. Interestingly, most of these candidate genes were reported to be involved in growth-related processes such as muscle formation, lipid deposition, feed efficiency, carcass composition and development of the central and peripheral nervous system. Our results offer novel insight into the molecular architecture underpinning yak growth traits. Further functional analyses are thus warranted to explore the molecular mechanisms whereby these genes affect these traits of interest.

Effect of divergence in residual feed intake on expression of lipid metabolism-related genes in the liver of meat-type ducks1

Lipid metabolism is considered one of the important factors affecting residual feed intake (RFI). However, the relationship between RFI and expression of lipid metabolism-related genes is unknown in meat-type ducks. To address this issue, a total of 1,000 male meat-type ducks with similar body weight were randomly selected to measure body weight gain and feed intake from 21 to 42 d of age to estimate RFI. The 8 greatest- (high RFI [HRFI]) and lowest- (low RFI [LRFI]) ranking birds were then selected for the present study. Relative expressions of key genes, namely sirtuin 1 (Sirt1), forkhead box O1 (Foxo1), peroxisome proliferator-activated receptor gamma (PPARγ), sterol regulatory element-binding transcription factor 1c (SREBP-1c), fas cell surface death receptor (FAS), acetyl-CoA carboxylase alpha (ACC), carnitine palmitoyltransferase 1A (CPT1A), and acyl-CoA oxidase 1 (ACOX1), were then determined in the HRFI and LRFI ducks by quantitative PCR. The results showed that RFI, feed conversion ratio (FCR), and average daily feed intake (ADFI) were significantly lower (P < 0.05) in LRFI ducks than in HRFI ducks. In addition, expression of Sirt1, Foxo1, CPT1A, and ACOX1 were significantly higher in LRFI ducks than in HRFI ducks (P < 0.05), whereas PPARγ and FAS expression levels were significantly lower in LRFI ducks than in HRFI ducks (P < 0.01). Correlation analysis showed that Sirt1, CPT1A, and ACOX1 expressions were significantly negatively correlated with FCR (r = -0.81 to -0.93; P < 0.01), whereas PPARγ and FAS expressions were significantly positively correlated with FCR (r = 0.74 to 0.87; P < 0.01). PPARγ expression was significantly positively correlated with RFI (r = 0.83; P < 0.01), whereas CPT1A and ACOX1 expressions were significantly negatively correlated with RFI (r = -0.84 to -0.89; P < 0.01). Sirt1 mRNA expression was positively correlated with Foxo1, CPT1A, and ACOX1 mRNA expression (r = 0.78 to 0.92; P < 0.01). Association of Foxo1 with CPT1A and ACOX1 was positive (r = 0.88 to 0.96; P < 0.01). These results suggest that genes related to fatty acid oxidation are upregulated in the liver of ducks with high feed efficiency, while genes associated with lipid synthesis are downregulated. Furthermore, the inclusion of lipid metabolism-related genes in future breeding programs might be beneficial for selecting ducks with greater feed efficiency phenotype.

Combining information from genome-wide association and multi-tissue gene expression studies to elucidate factors underlying genetic variation for residual feed intake in Australian Angus cattle

BACKGROUND

Genome-wide association studies (GWAS) are extensively used to identify single nucleotide polymorphisms (SNP) underlying the genetic variation of complex traits. However, much uncertainly often still exists about the causal variants and genes at quantitative trait loci (QTL). The aim of this study was to identify QTL associated with residual feed intake (RFI) and genes in these regions whose expression is also associated with this trait. Angus cattle (2190 steers) with RFI records were genotyped and imputed to high density arrays (770 K) and used for a GWAS approach to identify QTL associated with RFI. RNA sequences from 126 Angus divergently selected for RFI were analyzed to identify the genes whose expression was significantly associated this trait with special attention to those genes residing in the QTL regions.

RESULTS

The heritability for RFI estimated for this Angus population was 0.3. In a GWAS, we identified 78 SNPs associated with RFI on six QTL (on BTA1, BTA6, BTA14, BTA17, BTA20 and BTA26). The most significant SNP was found on chromosome BTA20 (rs42662073) and explained 4% of the genetic variance. The minor allele frequencies of significant SNPs ranged from 0.05 to 0.49. All regions, except on BTA17, showed a significant dominance effect. In 1 Mb windows surrounding the six significant QTL, we found 149 genes from which OAS2, STC2, SHOX, XKR4, and SGMS1 were the closest to the most significant QTL on BTA17, BTA20, BTA1, BTA14, and BTA26, respectively. In a 2 Mb windows around the six significant QTL, we identified 15 genes whose expression was significantly associated with RFI: BTA20) NEURL1B and CPEB4; BTA17) RITA1, CCDC42B, OAS2, RPL6, and ERP29; BTA26) A1CF, SGMS1, PAPSS2, and PTEN; BTA1) MFSD1 and RARRES1; BTA14) ATP6V1H and MRPL15.

CONCLUSIONS

Our results showed six QTL regions associated with RFI in a beef Angus population where five of these QTL contained genes that have expression associated with this trait. Therefore, here we show that integrating information from gene expression and GWAS studies can help to better understand the genetic mechanisms that determine variation in complex traits.

Transcriptome Analysis Identifies Candidate Genes and Pathways Associated With Feed Efficiency in Hu Sheep

In the genetic improvement of livestock and poultry, residual feed intake (RFI) is an important economic trait. However, in sheep, the genetic regulatory mechanisms of RFI are unclear. In the present study, we measured the feed efficiency (FE)-related phenotypes of 137 male Hu lambs, and selected six lambs with very high (n = 3) and very low (n = 3) RFI values and analyzed their liver transcriptomes. A total of 101 differentially expressed genes were identified, of which 40 were upregulated and 61 were downregulated in the low-RFI group compared with that in the high-RFI group. The downregulated genes were mainly concentrated in immune function pathways, while the upregulated genes were mainly involved in energy metabolism pathways. Two differentially expressed genes, ADRA2A (encoding adrenoceptor alpha 2A) and RYR2 (ryanodine receptor 2), were selected as candidate genes for FE and subjected to single nucleotide polymorphism scanning and association analysis. Two synonymous mutations, ADRA2A g.1429 C > A and RYR2 g.1117 A > C, were detected, which were both significantly associated with the feed conversion rate. These findings provide a deeper understanding of the molecular mechanisms regulating FE, and reveal key genes and genetic variants that could be used to genetically improve FE in sheep.

Tissue-specific responses of antioxidant pathways to poor hygiene conditions in growing pigs divergently selected for feed efficiency

Background

Poor hygiene of housing induces a systemic inflammatory response. Because inflammation and oxidative stress are processes that can sustain each other, the ways pigs are able to activate their antioxidant defenses are critical for production performance and health during periods when the immune system is solicited. Selection for production performance can also influence reactive oxygen species (ROS) production and expression levels of genes involved in cellular response to oxidative stress in different tissues. To establish the extent by which poor hygiene and selection for feed efficiency affected redox status, pigs divergently selected for residual feed intake (RFI) were housed in poor or good hygiene during 6 weeks. At the end, blood was collected in all pigs, and half of them were killed for tissue sampling. The remaining pigs were reared in good hygiene conditions during a recovery period of 7–8 weeks.

Results

At week 6, poor hygiene was associated with a lower total antioxidant capacity assessed by plasma ferric reducing ability in all pigs, and with greater plasma levels of hydrogen peroxides in the high RFI pigs (less efficient). Adipose tissue of high RFI pigs exhibited higher activities of catalase and glutathione reductase, and greater thiobarbituric acid reactive substances (TBARS) concentrations when compared with the low RFI pigs (more efficient). Poor hygiene conditions activated the antioxidant enzymes activities (glutathione reductase, superoxide dismutase and catalase) in adipose tissue of both lines, but led to higher ROS production by mature adipocytes isolated from the high RFI pigs only. In liver and muscle, there were only minor changes in antioxidant molecules due to genetics and hygiene conditions. After the resilience period, adipose tissue of pigs previously challenged by poor hygiene maintained higher antioxidant enzyme activities, and for the high RFI line, displayed higher TBARS concentrations.

Conclusions

Pigs selected for improved feed efficiency showed a lower susceptibility to oxidative stress induced by poor hygiene conditions. This could led to a lower inflammatory response and less impaired growth when these pigs are facing sanitary challenges during the production period.

Acute systemic inflammatory response to lipopolysaccharide stimulation in pigs divergently selected for residual feed intake

Background

It is unclear whether improving feed efficiency by selection for low residual feed intake (RFI) compromises pigs’ immunocompetence. Here, we aimed at investigating whether pig lines divergently selected for RFI had different inflammatory responses to lipopolysaccharide (LPS) exposure, regarding to clinical presentations and transcriptomic changes in peripheral blood cells.

Results

LPS injection induced acute systemic inflammation in both the low-RFI and high-RFI line (n = 8 per line). At 4 h post injection (hpi), the low-RFI line had a significantly lower (p = 0.0075) mean rectal temperature compared to the high-RFI line. However, no significant differences in complete blood count or levels of several plasma cytokines were detected between the two lines. Profiling blood transcriptomes at 0, 2, 6, and 24 hpi by RNA-sequencing revealed that LPS induced dramatic transcriptional changes, with 6296 genes differentially expressed at at least one time point post injection relative to baseline in at least one line (n = 4 per line) (|log₂(fold change)| ≥ log₂(1.2); q < 0.05). Furthermore, applying the same cutoffs, we detected 334 genes differentially expressed between the two lines at at least one time point, including 33 genes differentially expressed between the two lines at baseline. But no significant line-by-time interaction effects were detected. Genes involved in protein translation, defense response, immune response, and signaling were enriched in different co-expression clusters of genes responsive to LPS stimulation. The two lines were largely similar in their peripheral blood transcriptomic responses to LPS stimulation at the pathway level, although the low-RFI line had a slightly lower level of inflammatory response than the high-RFI line from 2 to 6 hpi and a slightly higher level of inflammatory response than the high-RFI line at 24 hpi.

Conclusions

The pig lines divergently selected for RFI had a largely similar response to LPS stimulation. However, the low-RFI line had a relatively lower-level, but longer-lasting, inflammatory response compared to the high-RFI line. Our results suggest selection for feed efficient pigs does not significantly compromise a pig’s acute systemic inflammatory response to LPS, although slight differences in intensity and duration may occur.

Host transcriptome and microbiome interaction modulates physiology of full-sibs broilers with divergent feed conversion ratio

Efficient livestock production relies on effective conversion of feed into body weight gain (BWG). High levels of feed conversion are especially important in production of broiler chickens, birds reared for meat, where economic margins are tight. Traits associated with improved broiler growth and feed efficiency have been subjected to intense genetic selection, but measures such as feed conversion ratio (FCR) remain variable, even between full siblings (sibs). Non-genetic factors such as the composition and function of microbial populations within different enteric compartments have been recognized to influence FCR, although the extent of interplay between hosts and their microbiomes is unclear. To examine host–microbiome interactions we investigated variation in the composition and functions of host intestinal-hepatic transcriptomes and the intestinal microbiota of full-sib broilers with divergent FCR. Progeny from 300 broiler families were assessed for divergent FCR set against shared genetic backgrounds and exposure to the same environmental factors. The seven most divergent full-sib pairs were chosen for analysis, exhibiting marked variation in transcription of genes as well as gut microbial diversity. Examination of enteric microbiota in low FCR sibs revealed variation in microbial community structure and function with no difference in feed intake compared to high FCR sibs. Gene transcription in low and high FCR sibs was significantly associated with the abundance of specific microbial taxa. Highly intertwined interactions between host transcriptomes and enteric microbiota are likely to modulate complex traits like FCR and may be amenable to selective modification with relevance to improving intestinal homeostasis and health.

A Transcriptome Analysis Identifies Biological Pathways and Candidate Genes for Feed Efficiency in DLY Pigs

Feed cost accounts for approximately 65–75% of overall commercial pork production costs. Therefore, improving the feed efficiency of pig production is important. In this study, 12 individuals with either extremely high (HE) or low (LE) feed efficiency were selected from 225 Duroc × (Landrace × Yorkshire) (DLY) pigs. After the pigs were slaughtered, we collected small intestine mucosal tissue. Next, RNA sequencing (RNA-seq) analysis was used to reveal the presence and quantity of genes expressed between these extremely HE- and LE-groups. We found 433 significantly differentially expressed genes (DEGs) between the HE- and LE-groups. Of these, 389 and 44 DEGs were upregulated and downregulated in the HE-group, respectively. An enrichment analysis showed that the DEGs were mainly enriched in functions related to apical plasma membrane composition, transporter activity, transport process and hormone regulation of digestion and absorption. Protein network interaction and gene function analyses revealed that SLC2A2 was an important candidate gene for FE in pigs, which may give us a deeper understanding of the mechanism of feed efficiency. Furthermore, some significant DEGs identified in the current study could be incorporated into artificial selection programs for increased feeding efficiency in pigs.

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.

The effect of breed and diet type on the global transcriptome of hepatic tissue in beef cattle divergent for feed efficiency

Background

Feed efficiency is an important economic and environmental trait in beef production, which can be measured in terms of residual feed intake (RFI). Cattle selected for low-RFI (feed efficient) have similar production levels but decreased feed intake, while also emitting less methane. RFI is difficult and expensive to measure and is not widely adopted in beef production systems. However, development of DNA-based biomarkers for RFI may facilitate its adoption in genomic-assisted breeding programmes. Cattle have been shown to re-rank in terms of RFI across diets and age, while also RFI varies by breed. Therefore, we used RNA-Seq technology to investigate the hepatic transcriptome of RFI-divergent Charolais (CH) and Holstein-Friesian (HF) steers across three dietary phases to identify genes and biological pathways associated with RFI regardless of diet or breed.

Results

Residual feed intake was measured during a high-concentrate phase, a zero-grazed grass phase and a final high-concentrate phase. In total, 322 and 33 differentially expressed genes (DEGs) were identified across all diets for CH and HF steers, respectively. Three genes, GADD45G, HP and MID1IP1, were differentially expressed in CH when both the high-concentrate zero-grazed grass diet were offered. Two canonical pathways were enriched across all diets for CH steers. These canonical pathways were related to immune function.

Conclusions

The absence of common differentially expressed genes across all dietary phases and breeds in this study supports previous reports of the re-ranking of animals in terms of RFI when offered differing diets over their lifetime. However, we have identified biological processes such as the immune response and lipid metabolism as potentially associated with RFI divergence emphasising the previously reported roles of these biological processes with respect to RFI.

Electronic supplementary material

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