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

Transcriptome analysis of liver and ileum reveals potential regulation of long non-coding RNA in pigs with divergent feed efficiency

OBJECTIVE

Long non-coding RNA (LncRNA) plays a significant role in regulating feed efficiency. This study aims to explore the key lncRNAs, associated genes, and pathways in pigs with extreme feed efficiencies.

METHODS

We screened pigs with extremely high and low residual feed intake through a 12-week animal growth trial and then conducted transcriptome analysis on their liver and ileum tissues. We analyzed the differential expressed lncRNAs, microRNAs (miRNAs), and messenger RNAs through target gene prediction and functional analysis. And we identified key lncRNAs and their potential regulatory genes associated with feed efficiency through the construction of competitive endogenous RNA network.

RESULTS

Differentially expressed lncRNAs were pinpointed in the liver, revealing 23 crucial target genes primarily associated with guanosine triphosphate metabolism and glycolipid biosynthesis. In the ileum, a screening identified 92 pivotal target genes, mainly linked to lipid and small molecule metabolism. Moreover, LOC106504303 and LOC102160805 emerged as potentially significant lncRNAs respectively, playing roles in mitochondrial oxidative phosphorylation in the liver, and lipid and cholesterol metabolism in the ileum.

CONCLUSION

The lncRNAs regulate energy metabolism and biosynthesis in the liver, and the digestive absorption capacity in the small intestine, affecting the feed efficiency of pigs.

Integration of epigenomic and genomic data to predict residual feed intake and the feed conversion ratio in dairy sheep via machine learning algorithms

BACKGROUND

Feed efficiency (FE) is an essential trait in livestock species because of the constant demand to increase the productivity and sustainability of livestock production systems. A better understanding of the biological mechanisms associated with FEs might help improve the estimation and selection of superior animals. In this work, differentially methylated regions (DMRs) were identified via genome-wide bisulfite sequencing (GWBS) by comparing the DNA methylation profiles of milk somatic cells from dairy ewes that were divergent in terms of residual feed intake. The DMRs were identified by comparing divergent groups for residual feed intake (RFI), the feed conversion ratio (FCR), and the consensus between both metrics (Cons). Additionally, the predictive performance of these DMRs and genetic variants mapped within these regions was evaluated via three machine learning (ML) models (xgboost, random forest (RF), and multilayer feedforward artificial neural network (deeplearning)). The average performance of each model was based on the root mean squared error (RMSE) and squared Spearman correlation (rho2). Finally, the best model for each scenario was selected on the basis of the highest ratio between rho2 and RMSE.

RESULTS

In total, 12,257, 9,328, and 6,723 genes were annotated for DMRs detected in the RFI, FCR, and Cons groups, respectively. These genes are associated with important pathways for regulating FE in dairy sheep, such as protein digestion and absorption, hormone synthesis and secretion, control of energy availability, cellular signaling, and feed behavior pathways. With respect to the ML predictions, the smallest mean RMSE (0.17) was obtained using RF, which was used to predict RFI. The highest mean rho2 (0.20) was obtained when the RFI was predicted via the mean methylation within the DMRs identified, the consensus groups were compared, and the genetic variants mapped within these DMRs were included. The best overall models were obtained for the prediction of RFI using the DMRs obtained in the comparison of RFI groups (RMSE = 0.10, rho2 = 0.86) using xgboost and the DMRs plus the genetic variants identified via the Cons groups (RMSE = 0.07, rho2 = 0.62) using RF.

CONCLUSIONS

The results provide new insights into the biological mechanisms associated with FE and the control of these processes through epigenetic mechanisms. Additionally, the potential use of epigenetic information as a biomarker for the prediction of FE can be suggested based on the obtained results.

Multi-omics analysis reveals critical cis-regulatory roles of transposable elements in livestock genomes

Transposable elements (TEs) are important sources of genetic and regulatory variation, yet their functional roles in domesticated animals remain insufficiently explored. To address this gap, we comprehensively annotated TE types, ages, and distributions in the genomes of pig (Sus scrofa), cattle (Bos taurus), and chicken (Gallus gallus). Our analysis revealed species-specific patterns in TE abundance, amplification, and activity in modern genomes. By integrating transcriptomic and epigenomic data, we explored the impact of specific TE types on cis-regulatory elements (CREs) and constructed a TE expression atlas across five tissues in all three species. Our findings underscored the critical roles of tissue-specific TE expression and chromatin accessibility in regulating tissue-specific biological processes. Most notably, we developed a computational framework to uncover TE-mediated gene regulatory networks (TE-GRNs). Our findings provide valuable insights into the regulatory functions of TEs in livestock and offer a robust approach for studying TE-GRNs in diverse biological contexts.

Review: Feed efficiency and metabolic flexibility in livestock

Improving the conversion of feed into product has been a key focus of genetic improvement in all livestock species. Livestock feed efficiency is the amount of product produced per unit of feed intake. Feed efficiency also depends on processes that are not directly related to economically important phenotypes, which can be considered 'waste' from a production point of view but are vital maintenance-related functions that are closely associated with environmental flexibility and adaptation. Resource allocation theory suggests that an animal's resource budget is narrowed when production efficiency is improved through an increase in productive output, along with a decrease in feed intake (capacity) and body reserves (improved leanness). The resulting trade-offs between productivity and vital functions may render the animal less capable of responding to unexpected challenges, potentially leading to negative side effects that are not directly related to economically important phenotypes. However, selection for feed efficiency may not narrow the metabolic space and result in trade-offs if the increase in feed efficiency is the result of increased metabolic flexibility in fuel substrate choice (carbohydrates, lipids, and/or proteins) and other energy-saving strategies. This review evaluates the relationship between metabolic flexibility and feed efficiency during anabolism (growth), fasting, immune activation, general stress, and heat stress, with a focus on pig production. We start with a brief overview of energy processes and substrate metabolism of carbohydrates, lipids, and protein. During muscle metabolism, the type of fuel used depends on fibre type characteristics of the muscle. Selection for improved meat production has resulted in pigs with a greater abundance of fast-twitch fibres with lower energy expenditure and higher metabolic efficiency. Metabolic flexibility for adaptation to disease, and response to regular stress implies that a more reactive immune response and reduced fear response results in higher feed efficiency. The examples presented in this review show that selection for improved feed efficiency does not necessarily narrow the metabolic space and result in trade-offs between productivity and vital functions because of energy-sparing mechanisms.

Expression genome-wide association study identifies key regulatory variants enriched with metabolic and immune functions in four porcine tissues

BACKGROUND

Integration of high throughput DNA genotyping and RNA-sequencing data enables the discovery of genomic regions that regulate gene expression, known as expression quantitative trait loci (eQTL). In pigs, efforts to date have been mainly focused on purebred lines for traits with commercial relevance as such growth and meat quality. However, little is known on genetic variants and mechanisms associated with the robustness of an animal, thus its overall health status. Here, the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred female finishers were studied, with the aim of identifying novel eQTL regulatory regions and transcription factors (TFs) associated with regulation of porcine metabolism and health-related traits.

RESULTS

An expression genome-wide association study with 535,896 genotypes and the expression of 12,680 genes in liver, 13,310 genes in lung, 12,650 genes in spleen, and 12,595 genes in muscle resulted in 4,293, 10,630, 4,533, and 6,871 eQTL regions for each of these tissues, respectively. Although only a small fraction of the eQTLs were annotated as cis-eQTLs, these presented a higher number of polymorphisms per region and significantly stronger associations with their target gene compared to trans-eQTLs. Between 20 and 115 eQTL hotspots were identified across the four tissues. Interestingly, these were all enriched for immune-related biological processes. In spleen, two TFs were identified: ERF and ZNF45, with key roles in regulation of gene expression.

CONCLUSIONS

This study provides a comprehensive analysis with more than 26,000 eQTL regions identified that are now publicly available. The genomic regions and their variants were mostly associated with tissue-specific regulatory roles. However, some shared regions provide new insights into the complex regulation of genes and their interactions that are involved with important traits related to metabolism and immunity.

The impact of dietary Laminaria digitata and alginate lyase supplementation on the weaned piglet liver: A comprehensive proteomics and metabolomics approach

The brown seaweed Laminaria digitata, a novel feedstuff for weaned piglets, has potentially beneficial prebiotic properties. However, its recalcitrant cell wall challenges digestion in monogastrics. Alginate lyase is a promising supplement to mitigate this issue. This study's aim was to investigate the impact of incorporating 10% dietary Laminaria digitata, supplemented with alginate lyase, on the hepatic proteome and metabolome of weaned piglets. These diets introduced minor variations to the metabolome and caused significant shifts in the proteome. Dietary seaweed provided a rich source of n-3 PUFAs that could signal hepatic fatty acid oxidation (FABP, ACADSB and ALDH1B1). This may have affected the oxidative stability of the tissue, requiring an elevated abundance of GST for regulation. The presence of reactive oxygen species likely inflicted protein damage, triggering increased proteolytic activity (LAPTM4B and PSMD4). Alginate lyase supplementation augmented the number of differentially abundant proteins, which included GBE1 and LDHC, contributing to maintain circulating glucose levels by mobilizing glycogen stores and branched-chain amino acids. The enzymatic supplementation with alginate lyase amplified the effects of the seaweed-only diet. An additional filter was employed to test the effect of missing values on the proteomics analysis, which is discussed from a technical perspective. SIGNIFICANCE: Brown seaweeds such as Laminaria digitata have prebiotic and immune-modulatory components, such as laminarin, that can improve weaned piglet health. However, they have recalcitrant cell wall polysaccharides, such as alginate, that can elicit antinutritional effects on the monogastric digestive system. The aim of this study was to evaluate the effect of a high level of dietary L. digitata and alginate lyase supplementation on the hepatic metabolism of weaned piglets, using high throughput Omics approaches.

1H-NMR based-metabolomics reveals alterations in the metabolite profiles of chickens infected with ascarids and concurrent histomonosis infection

BACKGROUND

Gut infections of chickens caused by Ascaridia galli and Heterakis gallinarum are associated with impaired host performance, particularly in high-performing genotypes. Heterakis gallinarum is also a vector of Histomonas meleagridis that is often co-involved with ascarid infections. Here, we provide a first insight into the alteration of the chicken plasma and liver metabolome as a result of gastrointestinal nematode infections with concomitant histomonosis. 1H nuclear magnetic resonance (1H-NMR) based-metabolomics coupled with a bioinformatics analysis was applied to explore the variation in the metabolite profiles of the liver (N = 105) and plasma samples from chickens (N = 108) experimentally infected with A. galli and H. gallinarum (+H. meleagridis). This was compared with uninfected chickens at different weeks post-infection (wpi 2, 4, 6, 10, 14, 18) representing different developmental stages of the worms.

RESULTS

A total of 31 and 54 metabolites were quantified in plasma and aqueous liver extracts, respectively. Statistical analysis showed no significant differences (P > 0.05) in any of the 54 identified liver metabolites between infected and uninfected hens. In contrast, 20 plasma metabolites including, amino acids, sugars, and organic acids showed significantly elevated concentrations in the infected hens (P < 0.05). Alterations of plasma metabolites occurred particularly in wpi 2, 6 and 10, covering the pre-patent period of worm infections. Plasma metabolites with the highest variation at these time points included glutamate, succinate, trimethylamine-N-oxide, myo-inositol, and acetate. Differential pathway analysis suggested that infection induced changes in (1) phenylalanine, tyrosine, and tryptophan metabolism, (2) alanine, aspartate and glutamate metabolism; and 3) arginine and proline metabolism (Pathway impact > 0.1 with FDR adjusted P-value < 0.05).

CONCLUSION

In conclusion, 1H-NMR based-metabolomics revealed significant alterations in the plasma metabolome of high performing chickens infected with gut pathogens-A. galli and H. gallinarum. The alterations suggested upregulation of key metabolic pathways mainly during the patency of infections. This approach extends our understanding of host interactions with gastrointestinal nematodes at the metabolic level.

Transcriptome analysis of divergent residual feed intake phenotypes in the M. longissimus thoracis et lumborum of Wannan Yellow rabbits

Introduction: Feed efficiency is an important economic trait in rabbit meat production. The identification of molecular mechanisms and candidate genes for feed efficiency may improve the economic and environmental benefits of the rabbit meat industry. As an alternative to the conventional feed conversion ratio, residual feed intake (RFI) can be used as an accurate indicator of feed efficiency. Methods: RNA sequencing was used to identify the differentially expressed genes (DEGs) in the M. longissimus thoracis et lumborum of eight Wannan Yellow rabbits with excessively high or low RFIs (HRFI or LRFI, respectively). Thereafter, Gene Ontology (GO) analysis, enrichment using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, gene set enrichment analysis (GSEA), and protein-protein interaction (PPI) network analysis was conducted. Results: In total, 445 DEGs were identified in the M. longissimus thoracis et lumborum of rabbits with high and low RFIs. The significantly enriched GO terms identified in these two groups were primarily involved in energy and mitochondrial metabolism and oxidation-reduction processes. KEGG analysis identified 11 significantly enriched pathways, including oxidative phosphorylation, PI3K-Akt signaling, and extracellular matrix-receptor interaction pathways. According to GSEA, the expressions of genes and pathways related to mitochondrial function were upregulated in HRFI rabbits, whereas genes with upregulated expressions in LRFI rabbits were related to immune response and energy metabolism. Additionally, PPI network analysis revealed five potential candidate genetic markers. Conclusion: Comparative analysis of the M. longissimus thoracis et lumborum transcriptomes in HRFI and LRFI rabbits revealed FOS, MYC, PRKACB, ITGA2, and FN1 as potential candidate genes that affect feed efficiency in rabbits. In addition, key signaling pathways involved in oxidative phosphorylation and PI3K-Akt and ECM-receptor interaction signaling impact rabbit feed efficiency. These findings will aid in breeding programs to improve feed efficiency and optimize RFI selection of rabbits for meat production.

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.

Growth performance, carcass traits, meat quality, and blood variables of small-sized meat ducks with different feed efficiency phenotypes

The study investigated the effects of feed efficiency (residual feed intake, RFI and residual intake and gain, RIG) on the production performance of small-sized meat ducks. Ninety ducks with intermediate and extreme (high and low) RFI values were selected from 1,083 male ducks of similar body weight, and the 3 groups were then redivided according to RIG. For both efficiency measures, the feed conversion ratio (FCR) and average daily feed intake (ADFI) of efficient ducks were significantly lower than those of inefficient ducks (P < 0.05), while the residual body weight gain (RG) was significantly higher in efficient ducks (P < 0.05). Inefficient-RFI animals showed greater skin fat yield (P < 0.05), but no other differences in carcass traits were observed (P > 0.05). RIG had positive effects on the pH1 value of the breast muscle (P < 0.05), but feed efficiency did not affect the other meat quality traits (P > 0.05). With regard to blood biochemical parameters, efficient ducks had significantly lower triglycerides (TG) (P < 0.05). Correlation analysis demonstrated that RFI was positively correlated with average daily feed intake and feed conversion ratio (P < 0.05), while RIG exhibited a strong negative correlation with both (P < 0.05). The average daily body weight gain was positively correlated with RIG (P < 0.05). RIG had a positive effect on the pH1 value of the breast muscle (P < 0.05). Furthermore, triglyceride and high-density lipoprotein cholesterol levels correlated with both efficiency classifications (P < 0.05). Overall, the efficiency measures did not affect the carcass and meat quality of small-sized meat ducks but could identify ducks with lower feed consumption and fast growth.

Profiling of Chromatin Accessibility in Pigs across Multiple Tissues and Developmental Stages

The study of chromatin accessibility across tissues and developmental stages is essential for elucidating the transcriptional regulation of various phenotypes and biological processes. However, the chromatin accessibility profiles of multiple tissues in newborn pigs and across porcine liver development remain poorly investigated. Here, we used ATAC-seq and rRNA-depleted RNA-seq to profile open chromatin maps and transcriptional features of heart, kidney, liver, lung, skeletal muscle, and spleen in newborn pigs and porcine liver tissue in the suckling and adult stages, respectively. Specifically, by analyzing a union set of protein-coding genes (PCGs) and two types of transcripts (lncRNAs and TUCPs), we obtained a comprehensive annotation of consensus ATAC-seq peaks for each tissue and developmental stage. As expected, the PCGs with tissue-specific accessible promoters had active transcription and were relevant to tissue-specific functions. In addition, other non-coding tissue-specific peaks were involved in both physical activity and the morphogenesis of neonatal tissues. We also characterized stage-specific peaks and observed a close association between dynamic chromatin accessibility and hepatic function transition during liver postnatal development. Overall, this study expands our current understanding of epigenetic regulation in mammalian tissues and organ development, which can benefit both economic trait improvement and improve the biomedical usage of pigs.

The components of tumor microenvironment as biomarker for immunotherapy in metastatic renal cell carcinoma

Substantial improvement in prognosis among metastatic renal cell carcinoma (mRCC) patients has been achieved, owing to the rapid development and utilization of immunotherapy. In particular, immune checkpoint inhibitors (ICIs) have been considered the backbone of systemic therapy for patients with mRCC alongside multi-targeted tyrosine kinase inhibitors (TKIs) in the latest clinical practice guidelines. However, controversies and challenges in optimal individualized treatment regarding immunotherapy remains still About 2/3 of the patients presented non-response or acquired resistance to ICIs. Besides, immune-related toxicities, namely immune-related adverse events, are still elusive and life-threatening. Thus, reliable biomarkers to predict immunotherapeutic outcomes for mRCC patients are needed urgently. Tumor microenvironment (TME), consisting of immune cells, vasculature, signaling molecules, and extracellular matrix and regulates tumor immune surveillance and immunological evasion through complex interplay, plays a critical role in tumor immune escape and consequently manipulates the efficacy of immunotherapy. Various studied have identified the different TME components are significantly associated with the outcome of mRCC patients receiving immunotherapy, making them potential valuable biomarkers in therapeutic guidance. The present review aims to summarize the latest evidence on the associations between the components of TME including immune cells, cytokines and extracellular matrix, and the therapeutic responses among mRCC patients with ICI-based treatment. We further discuss the feasibility and limitation of these components as biomarkers.

Integration of Non-Coding RNA and mRNA Profiles Reveals the Mechanisms of Rumen Development Induced by Different Types of Diet in Calves

Selecting suitable feed types and understanding the gastrointestinal digestive mechanism are helpful for the growth and health of calves in intensive dairy farming. However, the effects on rumen development of changing the molecular genetic basis and the regulatory mechanism by using different feed types are still unclear. Nine 7-day-old Holstein bull calves were randomly divided into GF (concentrate), GFF (alfalfa: oat grass = 3:2) and TMR (concentrate: alfalfa grass: oat grass: water = 0.30:0.12:0.08:0.50) diet experiment groups. Rumen tissue and serum samples were collected for physiological and transcriptomic analysis after 80 days. The results showed that serum α-amylase content and ceruloplasmin activity were significantly higher in the TMR group, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis ncRNAs and mRNAs were significantly enriched in the pathways of rumen epithelial development and stimulated rumen cell growth, including the Hippo signaling pathway, Wnt signaling pathway, thyroid hormone signaling pathway, ECM-receptor interaction and the absorption of protein and fat. The circRNAs/lncRNA-miRNAs-mRNA networks constructed, including novel_circ_0002471, novel_circ_0012104, TCONS_00946152, TCONS_00960915, bta-miR-11975, bta-miR-2890, PADI3 and CLEC6A, participated in metabolic pathways of lipid, immune system, oxidative stress and muscle development. In conclusion, the TMR diet could improve rumen digestive enzyme activities, stimulate rumen nutrient absorption and stimulate the DEGs related to energy homeostasis and microenvironment balance, and is thus better than the GF and GFF diets for promoting rumen growth and development.

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.

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.

Effective Quality Breeding Directions-Comparison and Conservative Analysis of Hepatic Super-Enhancers between Chinese and Western Pig Breeds

The transcriptional initiation of genes is closely bound to the functions of cis-regulatory elements, including promoters, typical enhancers (TEs), and recently-identified super-enhancers (SEs). In this study, we identified these cis-regulatory elements in the livers of two Chinese (Meishan and Enshi Black) and two Western (Duroc and Large White) pig breeds using ChIP-seq data, then explored their similarities and differences. In addition, we analyzed the conservation of SEs among different tissues and species (pig, human, and mouse). We observed that SEs were more significantly enriched by transcriptional initiation regions, TF binding sites, and SNPs than other cis-elements. Western breeds included fewer SEs in number, while more growth-related QTLs were associated with these SEs. Additionally, the SEs were highly tissue-specific, and were conserved in the liver among humans, pigs, and mice. We concluded that intense selection could concentrate functional SEs; thus, SEs could be applied as effective detection regions in genomic selection breeding.

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.

Effects of residual feed intake divergence on growth performance, carcass traits, meat quality, and blood biochemical parameters in small-sized meat ducks

Feed efficiency (FE) is a major economic trait of meat duck. This study aimed to evaluate the effects of residual feed intake (RFI) divergence on growth performance, carcass traits, meat quality, and blood biochemical parameters in small-sized meat ducks. A total of 500 healthy 21-day-old male ducks were housed in individual cages until slaughter at 63 d of age. The growth performance was determined for all the ducks. The carcass yield, meat quality, and blood biochemical parameters were determined for the selected 30 high-RFI (HRFI) and 30 low-RFI (LRFI) ducks. In terms of growth performance, the RFI, feed conversion ratio (FCR), and average daily feed intake (ADFI) were found to be significantly lower in the LRFI group (P < 0.01), whereas no differences were observed in the BW and body weight gain (P > 0.05). For slaughter performance, no differences were observed in the carcass traits between the LRFI and HRFI groups (P > 0.05). For meat quality, the shear force of breast muscle was significantly lower in the LRFI group (P < 0.05), while the other meat quality traits of breast and thigh muscles demonstrated no differences (P > 0.05). For blood biochemical parameters, the serum concentrations of triglycerides (TG) and glucose (GLU) were significantly lower in the LRFI group (P < 0.05), while the other parameters showed no differences (P > 0.05). The correlation analysis demonstrated a high positive correlation between RFI, FCR, and ADFI (P < 0.01). The RFI demonstrated a negative effect on the breast muscle and lean meat yields, but a positive effect on the shear force of breast muscle (P < 0.05). Further, the RFI demonstrated a positive effect on the TG and GLU levels (P < 0.05). These results indicate that the selection for low RFI could improve the FE of small-sized meat ducks without affecting the production performance. This study provides valuable insight into the biological processes underlying the variations in FE in small-sized meat ducks.

Residual feed intake in peripartal dairy cows is associated with differences in milk fat yield, ruminal bacteria, biopolymer hydrolyzing enzymes, and circulating biomarkers of immunometabolism

Residual feed intake (RFI) measures feed efficiency independent of milk production level, and is typically calculated using data past peak lactation. In the current study, we retrospectively classified multiparous Holstein cows (n = 320) from 5 of our published studies into most feed-efficient (M-eff) or least feed-efficient (L-eff) groups using performance data collected during the peripartal period. Objectives were to assess differences in profiles of plasma biomarkers of immunometabolism, relative abundance of key ruminal bacteria, and activities of digestive enzymes in ruminal digesta between M-eff and L-eff cows. Individual data from cows with ad libitum access to a total mixed ration from d -28 to d +28 relative to calving were used. A linear regression model including dry matter intake (DMI), energy-corrected milk (ECM), changes in body weight (BW), and metabolic BW was used to classify cows based on RFI divergence into L-eff (n = 158) and M-eff (n = 162). Plasma collected from the coccygeal vessel at various times around parturition (L-eff = 60 cows; M-eff = 47 cows) was used for analyses of 30 biomarkers of immunometabolism. Ruminal digesta collected via esophageal tube (L-eff = 19 cows; M-eff = 29 cows) was used for DNA extraction and assessment of relative abundance (%) of 17 major bacteria using real-time PCR, as well as activity of cellulase, amylase, xylanase, and protease. The UNIVARIATE procedure of SAS 9.4 (SAS Institute Inc.) was used for analyses of RFI coefficients. The MIXED procedure of SAS was used for repeated measures analysis of performance, milk yield and composition, plasma immunometabolic biomarkers, ruminal bacteria, and enzyme activities. The M-eff cows consumed less DMI during the peripartal period compared with L-eff cows. In the larger cohort of cows, despite greater overall BW for M-eff cows especially in the prepartum (788 vs. 764 kg), no difference in body condition score was detected due to RFI or the interaction of RFI × time. Milk fat content (4.14 vs. 3.75 ± 0.06%) and milk fat yield (1.75 vs. 1.62 ± 0.04 kg) were greater in M-eff cows. Although cumulative ECM yield did not differ due to RFI (1,138 vs. 1,091 ± 21 kg), an RFI × time interaction due to greater ECM yield was found in M-eff cows. Among plasma biomarkers studied, concentrations of nonesterified fatty acids, β-hydroxybutyrate, bilirubin, ceruloplasmin, haptoglobin, myeloperoxidase, and reactive oxygen metabolites were overall greater, and glucose, paraoxonase, and IL-6 were lower in M-eff compared with L-eff cows. Among bacteria studied, abundance of Ruminobacter amylophilus and Prevotella ruminicola were more than 2-fold greater in M-eff cows. Despite lower ruminal activity of amylase in M-eff cows in the prepartum, regardless of RFI, we observed a marked linear increase after calving in amylase, cellulase, and xylanase activities. Protease activity did not differ due to RFI, time, or RFI × time. Despite greater concentrations of biomarkers reflective of negative energy balance and inflammation, higher feed efficiency measured as RFI in peripartal dairy cows might be associated with shifts in ruminal bacteria and amylase enzyme activity. Further studies could help address such factors, including the roles of the liver and the mammary gland.

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.

Differential Gene Expression Associated with Soybean Oil Level in the Diet of Pigs

The aim of this study was to identify the differentially expressed genes (DEG) from the skeletal muscle and liver samples of animal models for metabolic diseases in humans. To perform the study, the fatty acid (FA) profile and RNA sequencing (RNA-Seq) data of 35 samples of liver tissue (SOY1.5, n = 17 and SOY3.0, n = 18) and 36 samples of skeletal muscle (SOY1.5, n = 18 and SOY3.0, n = 18) of Large White pigs were analyzed. The FA profile of the tissues was modified by the diet, mainly those related to monounsaturated (MUFA) and polyunsaturated (PUFA) FA. The skeletal muscle transcriptome analysis revealed 45 DEG (FDR 10%), and the functional enrichment analysis identified network maps related to inflammation, immune processes, and pathways associated with oxidative stress, type 2 diabetes, and metabolic dysfunction. For the liver tissue, the transcriptome profile analysis revealed 281 DEG, which participate in network maps related to neurodegenerative diseases. With this nutrigenomics study, we verified that different levels of soybean oil in the pig diet, an animal model for metabolic diseases in humans, affected the transcriptome profile of skeletal muscle and liver tissue. These findings may help to better understand the biological mechanisms that can be modulated by the diet.

Identification and characterization of circular RNAs in association with the feed efficiency in Hu lambs

Background

Circular RNA (circRNA), as a new members of noncoding RNA family, have vital functions in many biological processes by as microRNA sponges or competing endogenous RNAs (ceRNAs). However, little has been reported about the genetic mechanism of circRNAs regulation of feed efficiency in sheep.

Results

This study aimed to explore the expression of circRNAs in the liver of Hu sheep with High-RFI (High residual feed intake) and Low-RFI (Low residual feed intake) using transcriptome sequencing. A total of 20,729 circRNAs were identified in two groups, in which 219 circRNAs were found as significantly differentially expressed. Several circRNAs were validated by using RT-PCR, sanger sequencing and RT-qPCR methods. These results demonstrated that the RNA-seq result and expression level of circRNAs identified are reliable. Subsequently, GO and KEGG enrichment analysis of the parental genes of the differentially expressed (DE) circRNAs were mainly involved in immunity response and metabolic process. Finally, the ceRNA regulatory networks analysis showed that the target binding sites for miRNA such as novel_41, novel_115, novel_171 and oar-miR-485-3p in the identified DE cirRNAs. Importantly, two metabolic (SHISA3 and PLEKHH2) and four (RTP4, CD274, OAS1, and RFC3) immune-related target mRNAs were identified from 4 miRNAs. Association analysis showed that the polymorphism (RTP4 c.399 A > G) in the target gene RTP4 were significantly associated with RFI (P < 0.05).

Conclusions

Analysis of sequencing data showed some candidate ceRNAs that may play key roles in the feed efficiency in sheep by regulating animal immune and metabolic. These results provide the basis data for further study of the biological functions of circRNAs in regulating sheep feed efficiency.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08517-5.

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

Background

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

Results

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

Conclusion

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

Supplementary Information

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

Differences in Liver microRNA profiling in pigs with low and high feed efficiency

Feed cost is the main factor affecting the economic benefits of pig industry. Improving the feed efficiency (FE) can reduce the feed cost and improve the economic benefits of pig breeding enterprises. Liver is a complex metabolic organ which affects the distribution of nutrients and regulates the efficiency of energy conversion from nutrients to muscle or fat, thereby affecting feed efficiency. MicroRNAs (miRNAs) are small non-coding RNAs that can regulate feed efficiency through the modulation of gene expression at the post-transcriptional level. In this study, we analyzed miRNA profiling of liver tissues in High-FE and Low-FE pigs for the purpose of identifying key miRNAs related to feed efficiency. A total 212~221 annotated porcine miRNAs and 136~281 novel miRNAs were identified in the pig liver. Among them, 188 annotated miRNAs were co-expressed in High-FE and Low-FE pigs. The 14 miRNAs were significantly differentially expressed (DE) in the livers of high-FE pigs and low-FE pigs, of which 5 were downregulated and 9 were upregulated. Kyoto Encyclopedia of Genes and Genomes analysis of liver DE miRNAs in high-FE pigs and low-FE pigs indicated that the target genes of DE miRNAs were significantly enriched in insulin signaling pathway, Gonadotropin-releasing hormone signaling pathway, and mammalian target of rapamycin signaling pathway. To verify the reliability of sequencing results, 5 DE miRNAs were randomly selected for quantitative reverse transcription-polymerase chain reaction (qRT-PCR). The qRT-PCR results of miRNAs were confirmed to be consistent with sequencing data. DE miRNA data indicated that liver-specific miRNAs synergistically acted with mRNAs to improve feed efficiency. The liver miRNAs expression analysis revealed the metabolic pathways by which the liver miRNAs regulate pig feed efficiency.

Bioinformatics and Network Pharmacology-Based Approaches to Explore the Potential Mechanism of the Antidepressant Effect of Cyperi Rhizoma through Soothing the Liver

Major depressive disorder (MDD) has become the second most common disease worldwide, making it a threat to human health. Cyperi Rhizoma (CR) is a traditional herbal medicine with antidepressant properties. Traditional Chinese medicine theory states that CR relieves MDD by dispersing stagnated liver qi to soothe the liver, but the material basis and underlying mechanism have not been elucidated. In this study, we identified the active compounds and potential anti-MDD targets of CR by network pharmacology-based approaches. Through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, we hypothesized that the anti-MDD effect of CR may be mediated by an altered response of the liver to lipopolysaccharide (LPS) and glucose metabolism. Through bioinformatics analysis, comparing normal and MDD liver tissue in rats with spontaneous diabetes, we identified differentially expressed genes (DEGs) and selected PAI-1 (SERPINE1) as a target of CR in combating MDD. Molecular docking and molecular dynamics analysis also verified the binding of the active compound quercetin to PAI-1. It can be concluded that quercetin is the active compound of CR that acts against MDD by targeting PAI-1 to enhance the liver response to LPS and glucose metabolism. This study not only reveals the material basis and underlying mechanism of CR against MDD through soothing the liver but also provides evidence for PAI-1 as a potential target and quercetin as a potential agent for MDD treatment.

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.

Factors Influencing Proteolysis and Protein Utilization in the Intestine of Pigs: A Review

Pigs are among the most important farm animals for meat production worldwide. In order to meet the amino acid requirements of the animals, pigs rely on the regular intake of proteins and amino acids with their feed. Unfortunately, pigs excrete about two thirds of the used protein, and production of pork is currently associated with a high emission of nitrogen compounds resulting in negative impacts on the environment. Thus, improving protein efficiency in pigs is a central aim to decrease the usage of protein carriers in feed and to lower nitrogen emissions. This is necessary as the supply of plant protein sources is limited by the yield and the cultivable acreage for protein plants. Strategies to increase protein efficiency that go beyond the known feeding options have to be investigated considering the characteristics of the individual animals. This requires a deep understanding of the intestinal processes including enzymatic activities, capacities of amino acid transporters and the microbiome. This review provides an overview of these physiological factors and the respective analyses methods.

The Impact of Maternal and Piglet Low Protein Diet and Their Interaction on the Porcine Liver Transcriptome around the Time of Weaning

Maternal diet during early gestation affects offspring phenotype, but it is unclear whether maternal diet during late gestation influences piglet metabolism. We evaluated the impact of two dietary protein levels in sow late gestation diet and piglet nursery diet on piglet metabolism. Diets met or exceeded the crude protein and amino acid requirements. Sows received either 12% (Lower, L) or 17% (Higher, H) crude protein (CP) during the last five weeks of gestation, and piglets received 16.5% (L) or 21% (H) CP from weaning at age 3.5 weeks. This resulted in a 2 × 2 factorial design with four sow/piglet diet treatment groups: HH and LL (match), HL and LH (mismatch). Piglet hepatic tissues were sampled and differentially expressed genes (DEGs) were determined by RNA sequencing. At age 4.5 weeks, 25 genes were downregulated and 22 genes were upregulated in the mismatch compared to match groups. Several genes involved in catabolic pathways were upregulated in the mismatch compared to match groups, as were genes involved in lipid metabolism and inflammation. The results show a distinct interaction effect between maternal and nursery diets, implying that sow late gestation diet could be used to optimize piglet metabolism.

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.

Transcriptome Analysis Reveals Key Genes and Pathways Associated with Mummify Piglets

China is the country with the largest pork consumption in the world. However, the incidence of high mummify piglets (3-5%) is one of the important factors that cause the slow improvement of pig reproductive capacity, and the genetic mechanism is still unclear. This study aimed to identify candidate genes related to high mummify piglets. RNA-seq technology was used to comparative transcriptome profiling of blood from high piglets mummified and healthy sow at different stages of pregnancy (35d, 56d, 77d and 98d). A total of 137 to 420 DEGs were detected in each stage. Seven differentially expressed genes were significantly differentially expressed at various stages. IL-9R, TLR8, ABLIM3, FSH-α, ASCC1, PRKCZ, and GCK may play an important role in course of mummify piglets. The differential genes we identified between the groups were mainly enriched in immune and inflammation regulation, and others were mainly enriched in reproduction. Considering the function of candidate genes, IL-9R and TLR8 were suggested as the most promising candidate genes involved in mummify piglet traits. We speculate that during pregnancy, it may be the combined effects of the above-mentioned inflammation, immune response, and reproduction-related signal pathways that affect the occurrence of mummifying piglets, and further affect pig reproduction.

A hierarchical dynamic model used for investigating feed efficiency and its relationship with hepatic gene expression in APOE*3-Leiden.CETP mice

BACKGROUND

Feed efficiency (FE) is an important trait for livestock and humans. While the livestock industry focuses on increasing FE, in the current obesogenic society it is more of interest to decrease FE. Hence, understanding mechanisms involved in the regulation of FE and particularly how it can be decreased would help tremendously in counteracting the obesity pandemic. However, it is difficult to accurately measure or calculate FE in humans. In this study, we aimed to address this challenge by developing a hierarchical dynamic model based on humanized mouse data.

METHODS

We analyzed existing experimental data derived from 105 APOE*3-Leiden.CETP (E3L.CETP) mice fed a high-fat high-cholesterol (HFHC) diet for 1 (N = 20), 2 (N = 19), 3 (N = 20), and 6 (N = 46) month. We developed an ordinary differential equation (ODE) based model to estimate the FE based on the longitudinal data of body weight and food intake. Since the liver plays an important role in maintaining metabolic homeostasis, we evaluated associations between FE and hepatic gene expression levels. Depending on the feeding duration, we observed different relationships between FE and hepatic gene expression levels.

RESULTS

After 1-month feeding of HFHC diet, we observed that FE was associated with vitamin A metabolism, arachidonic acid metabolism, and the PPAR signaling pathway. After 3- and 6-month feeding of HFHC diet, we observed that FE was associated most strongly with expression levels of Spink1 and H19, genes involved in cell proliferation and glucose metabolism, respectively.

CONCLUSIONS

In conclusion, our analysis suggests that various biological processes such as vitamin A metabolism, hepatic response to inflammation, and cell proliferation associate with FE at different stages of diet-induced obesity.

Identification of functional candidate variants and genes for feed efficiency in Holstein and Jersey cattle breeds using RNA-sequencing

The identification of functional genetic variants and associated candidate genes linked to feed efficiency may help improve selection for feed efficiency in dairy cattle, providing economic and environmental benefits for the dairy industry. This study used RNA-sequencing data obtained from liver tissue from 9 Holstein cows [n = 5 low residual feed intake (RFI), n = 4 high RFI] and 10 Jersey cows (n = 5 low RFI, n = 5 high RFI), which were selected from a single population of 200 animals. Using RNA-sequencing, 3 analyses were performed to identify: (1) variants within low or high RFI Holstein cattle; (2) variants within low or high RFI Jersey cattle; and (3) variants within low or high RFI groups, which are common across both Holstein and Jersey cattle breeds. From each analysis, all variants were filtered for moderate, modifier, or high functional effect, and co-localized quantitative trait loci (QTL) classes, enriched biological processes, and co-localized genes related to these variants, were identified. The overlapping of the resulting genes co-localized with functional SNP from each analysis in both breeds for low or high RFI groups were compared. For the first two analyses, the total number of candidate genes associated with moderate, modifier, or high functional effect variants fixed within low or high RFI groups were 2,810 and 3,390 for Holstein and Jersey breeds, respectively. The major QTL classes co-localized with these variants included milk and reproduction QTL for the Holstein breed, and milk, production, and reproduction QTL for the Jersey breed. For the third analysis, the common variants across both Holstein and Jersey breeds, uniquely fixed within low or high RFI groups were identified, revealing a total of 86,209 and 111,126 functional variants in low and high RFI groups, respectively. Across all 3 analyses for low and high RFI cattle, 12 and 31 co-localized genes were overlapping, respectively. Among the overlapping genes across breeds, 9 were commonly detected in both the low and high RFI groups (INSRR, CSK, DYNC1H1, GAB1, KAT2B, RXRA, SHC1, TRRAP, PIK3CB), which are known to play a key role in the regulation of biological processes that have high metabolic demand and are related to cell growth and regeneration, metabolism, and immune function. The genes identified and their associated functional variants may serve as candidate genetic markers and can be implemented into breeding programs to help improve the selection for feed efficiency in dairy cattle.

Association of residual feed intake with growth performance, carcass traits, meat quality, and blood variables in native chickens

Improving feed efficiency is a primary goal in poultry breeding strategies. Residual feed intake (RFI) in chickens typically calculated during the growing period is a measure of feed efficiency that is independent of the level of production. The objective of this study was to evaluate phenotypic correlations of growth performance, carcass traits, meat quality, and blood variables with RFI in growing native chickens. A total of 1,008 chickens were selected for the experiment to derive RFI. After the RFI measurement period of 42 d, 25 chickens with low RFI values, 25 chickens with medium RFI values, and 25 chickens with high RFI values were selected. The RFI was significantly positively correlated with feed conversion ratio and average daily feed intake, while it was not significantly correlated with initial body weight (BW), final BW, average daily body weight gain, and metabolic BW0.75. The abdominal fat weight and yield of high RFI group were significantly greater than those of medium and low RFI groups, and the abdominal fat yield was significantly positively correlated with RFI. Moreover, the plasma insulin-like growth factor 1 (IGF-1) content of low RFI group was significantly greater than those of high and medium RFI groups. The plasma concentrations of adrenocorticotropic hormone, triiodothyronine (T3), and cortisol of high RFI birds were significantly greater than that of low RFI birds. RFI was significantly positively correlated with plasma concentrations of T3 and cortisol, while it was significantly negatively correlated with plasma concentration of IGF-1. In addition, the serum levels of glucose and triglyceride of high RFI birds were significantly lower than that of low RFI birds. The serum low-density lipoprotein cholesterol (LDL-C) content of high RFI group was significantly greater than that of medium and low RFI groups, and it was significantly positively correlated with RFI. Our data suggested that selection of chickens with low RFI values may be beneficial to reduce fat deposition in native chickens without affecting the meat quality. Circulating IGF-1, T3, cortisol, and LDL-C concentrations can be used as indirect selection indicators of feed efficiency in native chickens. The effect of IGF-1, T3, cortisol, and LDL-C on feed efficiency of native chickens should be carefully examined and validated in future breeding programs.

Integrative Analysis of Metabolomic and Transcriptomic Profiles Uncovers Biological Pathways of Feed Efficiency in Pigs

Feed efficiency (FE) is an economically important trait. Thus, reliable predictors would help to reduce the production cost and provide sustainability to the pig industry. We carried out metabolome-transcriptome integration analysis on 40 purebred Duroc and Landrace uncastrated male pigs to identify potential gene-metabolite interactions and explore the molecular mechanisms underlying FE. To this end, we applied untargeted metabolomics and RNA-seq approaches to the same animals. After data quality control, we used a linear model approach to integrate the data and find significant differently correlated gene-metabolite pairs separately for the breeds (Duroc and Landrace) and FE groups (low and high FE) followed by a pathway over-representation analysis. We identified 21 and 12 significant gene-metabolite pairs for each group. The valine-leucine-isoleucine biosynthesis/degradation and arginine-proline metabolism pathways were associated with unique metabolites. The unique genes obtained from significant metabolite-gene pairs were associated with sphingolipid catabolism, multicellular organismal process, cGMP, and purine metabolic processes. While some of the genes and metabolites identified were known for their association with FE, others are novel and provide new avenues for further research. Further validation of genes, metabolites, and gene-metabolite interactions in larger cohorts will elucidate the regulatory mechanisms and pathways underlying FE.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome-Wide Epistatic Interaction Networks Affecting Feed Efficiency in Duroc and Landrace Pigs

Interactions among genomic loci have often been overlooked in genome-wide association studies, revealing the combinatorial effects of variants on phenotype or disease manifestation. Unexplained genetic variance, interactions among causal genes of small effects, and biological pathways could be identified using a network biology approach. The main objective of this study was to determine the genome-wide epistatic variants affecting feed efficiency traits [feed conversion ratio (FCR) and residual feed intake (RFI)] based on weighted interaction SNP hub (WISH-R) method. Herein, we detected highly interconnected epistatic SNP modules, pathways, and potential biomarkers for the FCR and RFI in Duroc and Landrace purebreds considering the whole population, and separately for low and high feed efficient groups. Highly interacting SNP modules in Duroc (1,247 SNPs) and Landrace (1,215 SNPs) across the population and for low feed efficient (Duroc-80 SNPs, Landrace-146 SNPs) and high feed efficient group (Duroc-198 SNPs, Landrace-232 SNPs) for FCR and RFI were identified. Gene and pathway analyses identified ABL1, MAP3K4, MAP3K5, SEMA6A, KITLG, and KAT2B from chromosomes 1, 2, 5, and 13 underlying ErbB, Ras, Rap1, thyroid hormone, axon guidance pathways in Duroc. GABBR2, GNA12, and PRKCG genes from chromosomes 1, 3, and 6 pointed towards thyroid hormone, cGMP-PKG and cAMP pathways in Landrace. From Duroc low feed efficient group, the TPK1 gene was found involved with thiamine metabolism, whereas PARD6G, DLG2, CRB1 were involved with the hippo signaling pathway in high feed efficient group. PLOD1 and SETD7 genes were involved with lysine degradation in low feed efficient group in Landrace, while high feed efficient group pointed to genes underpinning valine, leucine, isoleucine degradation, and fatty acid elongation. Some SNPs and genes identified are known for their association with feed efficiency, others are novel and potentially provide new avenues for further research. Further validation of epistatic SNPs and genes identified here in a larger cohort would help to establish a framework for modelling epistatic variance in future methods of genomic prediction, increasing the accuracy of estimated genetic merit for FE and helping the pig breeding industry.

Prediction of nutrient digestibility in grower-finisher pigs based on faecal microbiota composition

Microbiota play an important role in total tract nutrient digestion, especially when fibrous diets are fed to pigs. This study aimed to use metagenomics to predict faecal nutrient digestibility in grower‐finisher pigs. The study design consisted of 160 three‐way crossbreed grower‐finisher pigs (80 female and 80 male) which were either fed a diet based on corn/soybean meal or a more fibrous diet based on wheat/barley/by‐products. On the day before slaughter, faecal samples were collected and used to determine faecal digestibility of dry matter, ash, organic matter, crude protein, crude fat, crude fibre and non‐starch polysaccharides. The faecal samples were also sequenced for the 16S hypervariable region of bacteria (V3/V4) to profile the faecal microbiome. With these data, we calculated the between‐animal variation in faecal nutrient digestibility associated with variation in the faecal microbiome, that is the “microbiability”. The microbiability values were significantly greater than zero for dry matter, organic matter, crude protein, crude fibre and non‐starch polysaccharides, ranging from 0.58 to 0.93, as well as for crude fat with a value of 0.37, but not significantly different from zero for ash. Using leave‐one‐out cross‐validation, we estimated the accuracy of predicting digestibility values of individual pigs based on their faecal microbiota composition. The accuracies of prediction for crude fat and ash digestibility were virtually 0, and for the other nutrients, the accuracies ranged from 0.42 to 0.63. In conclusion, the faecal microbiota composition gave high microbiability values for faecal digestibility of dry matter, organic matter, crude protein, crude fibre and non‐starch polysaccharides. The accuracies of prediction are relatively low if the interest is in precisely predicting faecal nutrient digestibility of individual pigs, but are promising from the perspective of ranking animals in a genetic selection context.

Identification and Conservation Analysis of Cis-Regulatory Elements in Pig Liver

The liver plays a key role in metabolism and affects pig production. However, the functional annotation of noncoding regions of the pig liver remains poorly understood. We revealed the landscape of cis-regulatory elements and their functional characterization in pig liver. We identified 102,373 cis-regulatory elements in the pig liver, including enhancers, promoters, super-enhancers, and broad H3K4me3 domains, and highlighted 26 core transcription regulatory factors in the pig liver as well. We found similarity of cis-regulatory elements among those of pigs, humans, and cattle. Despite the low proportion of functionally conserved enhancers (~30%) between pig and human liver tissue, ~78% of the pig liver enhancer orthologues sequence could play an enhancer role in other human tissues. Additionally, we observed that the ratio of consistent super-enhancer-associated genes was significantly higher than the ratio of functionally conserved super-enhancers. Approximately 54% of the core regulation factors driven by super-enhancers were consistent across the liver from these three species. Our pig liver annotation and functional characterization studies provide a system and resource for noncoding annotation for future gene regulatory studies in pigs. Furthermore, our study also showed the high level functional conservation of cis-regulatory elements in mammals; it also improved our understanding of regulation function of mammal cis-regulatory elements.