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

Metagenomic insights into the relationship between intestinal flora and residual feed intake of meat ducks

In this study, we sought to determine the effects of intestinal flora on the feed efficiency of meat ducks by evaluating the correlation between intestinal flora and residual feed intake. The F2 generation of Cherry Valley ducks × Runzhou Crested White ducks was used as the study subjects, and feed consumption being recorded from d 21 to 42. RFI was calculated based on growth performance, and 20 low RFI and 20 high RFI ducks were randomly selected to characterize the effect of RFI on growth performance. To analyze the intestinal flora affecting RFI, 16s rDNA sequencing was performed on the contents of 5 intestinal segments from the HR and LR groups, and macrogenomic sequencing was performed on the cecal contents. Feed intake, average daily feed intake, feed conversion ratio, and residual feed intake were lower in low RFI. Analysis of the intestinal flora revealed the cecum to be more highly enriched in the carbohydrate metabolism pathway and less enriched with potentially pathogenic taxa than the other assessed intestinal regions. Further analysis of the cecal microbiota identified nine significantly differentially enriched intestinal flora. In this study, we accordingly identified a basis for the mechanisms underlying the effects of the intestinal flora on meat duck feed efficiency.

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.

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.

Gene coexpression network analysis reveals perirenal adipose tissue as an important target of prenatal malnutrition in sheep

We have previously demonstrated that pre- and early postnatal malnutrition in sheep induced depot- and sex-specific changes in adipose morphological features, metabolic outcomes, and transcriptome in adulthood, with perirenal (PER) as the major target followed by subcutaneous (SUB) adipose tissue. We aimed to identify coexpressed and hub genes in SUB and PER to identify the underlying molecular mechanisms contributing to the early nutritional programming of adipose-related phenotypic outcomes. Transcriptomes of SUB and PER of male and female adult sheep with different pre- and early postnatal nutrition histories were used to construct networks of coexpressed genes likely to be functionally associated with pre- and early postnatal nutrition histories and phenotypic traits using weighted gene coexpression network analysis. The modules from PER showed enrichment of cell cycle regulation, gene expression, transmembrane transport, and metabolic processes associated with both sexes' prenatal nutrition. In SUB (only males), a module of enriched adenosine diphosphate metabolism and development correlated with prenatal nutrition. Sex-specific module enrichments were found in PER, such as chromatin modification in the male network but histone modification and mitochondria- and oxidative phosphorylation-related functions in the female network. These sex-specific modules correlated with prenatal nutrition and adipocyte size distribution patterns. Our results point to PER as a primary target of prenatal malnutrition compared to SUB, which played only a minor role. The prenatal programming of gene expression and cell cycle, potentially through epigenetic modifications, might be underlying mechanisms responsible for observed changes in PER expandability and adipocyte-size distribution patterns in adulthood in both sexes.

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.

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.

Thigh muscle metabolic response is linked to feed efficiency and meat characteristics in slow-growing chicken

The Korat chicken (KR) is a slow-growing Thai chicken breed with relatively poor feed efficiency (FE) but very tasty meat with high protein and low fat contents, and a unique texture. To enhance the competitiveness of KR, its FE should be improved. However, selecting for FE has an unknown effect on meat characteristics. Thus, understanding the genetic basis underlying FE traits and meat characteristics is needed. In this study, 75 male KR birds were raised up to 10 wk of age. For each bird, the feed conversion ratio (FCR), residual feed intake (RFI), and physicochemical properties, flavor precursors, and biological compounds in the thigh meat were evaluated. At 10 wk of age, thigh muscle samples from 6 birds (3 with high FCR and 3 with low FCR values) were selected, and their proteomes were investigated using a label-free proteomic method. Weighted gene coexpression network analysis (WGCNA) was used to screen the key protein modules and pathways. The WGCNA results revealed that FE and meat characteristics significantly correlated with the same protein module. However, the correlation was unfavorable; improving FE may result in a decrease in meat quality through the alteration in biological processes including glycolysis/gluconeogenesis, metabolic pathway, carbon metabolism, biosynthesis of amino acids, pyruvate metabolism, and protein processing in the endoplasmic reticulum. The hub proteins of the significant module (TNNT1, TNNT3, TNNI2, TNNC2, MYLPF, MYH10, GADPH, PGK1, LDHA, and GPI) were also identified to be associated with energy metabolism, and muscle growth and development. Given that the same proteins and pathways are present in FE and meat characteristics but in opposite directions, selection practices for KR should simultaneously consider both trait groups to maintain the high meat quality of KR while improving FE.

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) > 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.

RNA sequencing analysis of the longissimus dorsi to identify candidate genes underlying the intramuscular fat content in Anqing Six-end-white pigs

Intramuscular fat (IMF) is a significant marker for pork quality. The Anqing Six-end-white pig has the characteristics of high meat quality and IMF content. Owing to the influence of European commercial pigs and a late start in resource conservation, the IMF content within local populations varies between individuals. This study analyzed the longissimus dorsi transcriptome of purebred Anqing Six-end-white pigs with varying IMF content to recognize differentially expressed genes. We identified 1528 differentially expressed genes between the pigs with high (H) and low (L) IMF content. Based on these data, 1775 Gene Ontology terms were significantly enriched, including lipid metabolism, modification and storage, and regulation of lipid biosynthesis. Pathway analysis revealed 79 significantly enriched pathways, including the Peroxisome proliferator-activated receptor and mitogen-activated protein kinase signaling pathways. Moreover, gene set enrichment analysis indicated that the L group had increased the expression of genes related to ribosome function. Additionally, the protein-protein interaction network analyses revealed that VEGFA, KDR, LEP, IRS1, IGF1R, FLT1 and FLT4 were promising candidate genes associated with the IMF content. Our study identified the candidate genes and pathways involved in IMF deposition and lipid metabolism and provides data for developing local pig germplasm resources.

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.

Genetic and microbiome analysis of feed efficiency in laying hens

Improving feed efficiency is an important target for poultry breeding. Feed efficiency is affected by host genetics and the gut microbiota, but many of the mechanisms remain elusive in laying hens, especially in the late laying period. In this study, we measured feed intake, body weight, and egg mass of 714 hens from a pedigreed line from 69 to 72 wk of age and calculated the residual feed intake (RFI) and feed conversion ratio (FCR). In addition, fecal samples were also collected for 16S ribosomal RNA gene sequencing (V4 region). Genetic analysis was then conducted in DMU packages by using AI-REML with animal model. Moderate heritability estimates for FCR (h² = 0.31) and RFI (h² = 0.52) were observed, suggesting that proper selection programs can directly improve feed efficiency. Genetically, RFI was less correlated with body weight and egg mass than that of FCR. The phenotypic variance explained by gut microbial variance is defined as the microbiability (m²). The microbiability estimates for FCR (m² = 0.03) and RFI (m² = 0.16) suggested the gut microbiota was also involved in the regulation of feed efficiency. In addition, our results showed that the effect of host genetics on fecal microbiota was minor in three aspects: 1) microbial diversity indexes had low heritability estimates, and genera with heritability estimates more than 0.1 accounted for only 1.07% of the tested fecal microbiota; 2) the genetic relationship correlations between host genetics and different microbial distance were very weak, ranging from -0.0057 to -0.0003; 3) the microbial distance between different kinships showed no significant difference. Since the RFI has the highest microbiability, we further screened out three genera, including Anaerosporobacter, Candidatus Stoquefichus, and Fournierella, which were negatively correlated with RFI and played positive roles in improving the feed efficiency. These findings contribute to a great understanding of the genetic background and microbial influences on feed efficiency.

Identification of biomarkers associated with the feed efficiency by metabolomics profiling: results from the broiler lines divergent for high or low abdominal fat content

Background

Improving feed efficiency (FE) is one of the main objectives in broiler breeding. It is difficult to directly measure FE traits, and breeders hence have been trying to identify biomarkers for the indirect selection and improvement of FE traits. Metabolome is the "bridge" between genome and phenome. The metabolites may potentially account for more of the phenotypic variation and can suitably serve as biomarkers for selecting FE traits. This study aimed to identify plasma metabolite markers for selecting high-FE broilers. A total of 441 birds from Northeast Agricultural University broiler lines divergently selected for abdominal fat content were used to analyze plasma metabolome and estimate the genetic parameters of differentially expressed metabolites.

Results

The results identified 124 differentially expressed plasma metabolites (P < 0.05) between the lean line (high-FE birds) and the fat line (low-FE birds). Among these differentially expressed plasma metabolites, 44 were found to have higher positive or negative genetic correlations with FE traits (|r_g| ≥ 0.30). Of these 44 metabolites, 14 were found to display moderate to high heritability estimates (h² ≥ 0.20). However, among the 14 metabolites, 4 metabolites whose physiological functions have not been reported were excluded. Ultimately, 10 metabolites were suggested to serve as the potential biomarkers for breeding the high-FE broilers. Based on the physiological functions of these metabolites, reducing inflammatory and improving immunity were proposed to improve FE and increase production efficiency.

Conclusions

According to the pipeline for the selection of the metabolite markers established in this study, it was suggested that 10 metabolites including 7-ketocholesterol, dimethyl sulfone, epsilon-(gamma-glutamyl)-lysine, gamma-glutamyltyrosine, 2-oxoadipic acid, L-homoarginine, testosterone, adenosine 5'-monophosphate, adrenic acid, and calcitriol could be used as the potential biomarkers for breeding the "food-saving broilers".

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.

Transcriptomic analysis of the liver in aged laying hens with different eggshell strength

Eggshell is composed of a very ordered and mineralized structure and is important for egg quality. Eggshell strength is particularly important because of its direct association with economic outcomes and egg safety. Various factors related to laying hens and their environment affects eggshell strength. However, the molecular mechanisms of liver functions related to decreased eggshell strength of aged laying hens are largely unknown. Therefore, this study aimed to identify potential factors affecting eggshell strength in aged laying hens at the hepatic transcriptomic level. A total of five hundred 92-wk-old Hy-line Brown laying hens were screened to select those exhibiting the greatest variation in eggshell strength. Based on the final eggshell strength, 12 hens producing eggs with strong eggshell strength (SES) and weak eggshell strength (WES) were finally selected (n = 6) for liver tissue sampling. The RNA-sequencing was performed to identify differentially expressed genes (DEGs) between the 2 groups. We identified a total of 2,084 DEGs, of which 1,358 genes were upregulated and 726 genes were downregulated in the WES group compared with SES group. According to the Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis, the DEGs indicated the mammalian target of rapamycin signaling pathway, the Janus kinase-signal transducer and activator of transcription pathway, the mitogen‑activated protein kinase signaling pathway, and the insulin resistance pathways. Genes related to fatty liver disease were upregulated in WES group compared with SES group. In addition, expression of several genes associated with oxidative stress and bone resorption activity was altered in aged laying hens with different eggshell strength. Overall, these findings contribute to the identification of genes involved in different intensity of eggshell strength, enabling more understanding of the hepatic molecular mechanism underlying in decreased eggshell strength of aged laying hens.

EEF1A1 transcription cofactor gene polymorphism is associated with muscle gene expression and residual feed intake in Nelore cattle

Cis-acting effects of noncoding variants on gene expression and regulatory molecules constitute a significant factor for phenotypic variation in complex traits. To provide new insights into the impacts of single-nucleotide polymorphisms (SNPs) on transcription factors (TFs) and transcription cofactors (TcoF) coding genes, we carried out a multi-omic analysis to identify cis-regulatory effects of SNPs on these genes' expression in muscle and describe their association with feed efficiency-related traits in Nelore cattle. As a result, we identified one SNP, the rs137256008C > T, predicted to impact the EEF1A1 gene expression (β = 3.02; P-value = 3.51E-03) and the residual feed intake trait (β = - 3.47; P-value = 0.02). This SNP was predicted to modify transcription factor sites and overlaps with several QTL for feed efficiency traits. In addition, co-expression network analyses showed that animals containing the T allele of the rs137256008 SNP may be triggering changes in the gene network. Therefore, our analyses reinforce and contribute to a better understanding of the biological mechanisms underlying gene expression control of feed efficiency traits in bovines. The cis-regulatory SNP can be used as biomarker for feed efficiency in Nelore cattle.

Cecal Microbial Hydrogen Cycling Potential Is Linked to Feed Efficiency Phenotypes in Chickens

In chickens, early life exposure to environmental microbes has long-lasting impacts on gastrointestinal (GI) microbiome development and host health and growth, via mechanisms that remain uncharacterized. In this study, we demonstrated that administrating a fecal microbiome transplant (FMT) from adults to day-of-hatch chicks results in significantly higher body mass of birds and decreased residual feed intake (RFI), implying enhanced feed efficiency, at 6 weeks of age. To assess the potential mechanisms through which FMT affects adult bird phenotype, we combined 16 S rRNA gene amplification, metagenomic, and comparative genomic approaches to survey the composition and predicted activities of the resident microbiome of various GI tract segments. Early life FMT exposure had a long-lasting significant effect on the microbial community composition and function of the ceca but not on other GI segments. Within the ceca of 6-week-old FMT birds, hydrogenotrophic microbial lineages and genes were most differentially enriched. The results suggest that thermodynamic regulation in the cecum, in this case via hydrogenotrophic methanogenic and sulfur-cycling lineages, potentially serving as hydrogen sinks, may enhance fermentative efficiency and dietary energy harvest capacity. Our study provides a specific mechanism of action through which early-life microbiome transplants modulate market-relevant phenotypes in poultry and, thereby, may represent a significant advance toward microbiome-focused sustainable agriculture.

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 of Breast Muscle Reveals Pathways Related to Protein Deposition in High Feed Efficiency of Native Turkeys

Feed efficiency is important due to the high cost of food, which accounts for about 70% of the total cost of a turkey breeding system. Native poultry are an important genetic resource in poultry breeding programs. This study aimed to conduct a global transcriptome analysis of native male turkeys which have been phenotyped for high and low feed efficiency. Feed efficiency traits were recorded during the experimental period. After slaughter, the three most efficient and three least efficient male turkeys were selected for RNA-Seq analysis. A total of 365 genes with different expressions in muscle tissue were identified between turkeys with a high feed efficiency compared to turkeys with a low feed efficiency. In the pathway analysis of up-regulated genes, major pathways included the "metabolism of glycine, serine, and threonine"; the "adipocytokine signaling pathway" and the "biosynthesis of amino acids". In the pathway analysis of down-regulated genes, the major pathways included "dorso-ventral axis formation" and "actin cytoskeleton regulation". In addition, gene set enrichment analyses were performed, which showed that high feed efficiency birds exhibit an increased expression of genes related to the biosynthesis of amino acids and low feed efficiency birds an increased expression of genes related to the immune response. Furthermore, functional analysis and protein network interaction analysis revealed that genes including GATM, PSAT1, PSPH, PHGDH, VCAM1, CD44, KRAS, SRC, CAV3, NEDD9, and PTPRQ were key genes for feed efficiency. These key genes may be good potential candidates for biomarkers of feed efficiency in genetic selection in turkeys.

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.

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

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

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.

Revealing Pathways Associated with Feed Efficiency and Meat Quality Traits in Slow-Growing Chickens

Simple Summary

Korat is a new chicken breed with high-protein meat, low fat, and low purine content. However, the effects of improving the breed’s feed efficiency, which would enhance production, on meat quality are unknown. Hence, understanding the genetic architecture underlying feed efficiency and meat quality traits in chicken offers new opportunities toward genetic improvement. Through a weighted gene co-expression network analysis on Korat chickens, the presented results provide new information on the molecular pathways that play important roles in FE and meat quality that could help achieve the optimum feed efficiency while maintaining meat quality in Korat chickens.

Abstract

Here, molecular pathways and genes involved in the feed efficiency (FE) and thigh-meat quality of slow-growing Korat chickens were investigated. Individual feed intake values and body weights were collected weekly to the calculate feed conversion ratios (FCR) and residual feed intake. The biochemical composition and meat quality parameters were also measured. On the basis of extreme FCR values at 10 weeks of age, 9 and 12 birds from the high and the low FCR groups, respectively, were selected, and their transcriptomes were investigated using the 8 × 60 K Agilent chicken microarray. A weighted gene co-expression network analysis was performed to determine the correlations between co-expressed gene modules and FE, thigh-meat quality, or both. Groups of birds with different FE values also had different nucleotide, lipid, and protein contents in their thigh muscles. In total, 38 modules of co-expressed genes were identified, and 12 were correlated with FE and some meat quality traits. A functional analysis highlighted several enriched functions, such as biological processes, metabolic processes, nucleotide metabolism, and immune responses. Several molecular factors were involved in the interactions between FE and meat quality, including the assembly competence domain, baculoviral IAP repeat containing 5, cytochrome c oxidase assembly factor 3, and myosin light chain 9 genes.

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.

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

Simple Summary

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

Abstract

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

Transcriptomics analysis of differentially expressed genes in subcutaneous and perirenal adipose tissue of sheep as affected by their pre- and early postnatal malnutrition histories

BACKGROUND

Early life malnutrition is known to target adipose tissue with varying impact depending on timing of the insult. This study aimed to identify differentially expressed genes in subcutaneous (SUB) and perirenal (PER) adipose tissue of 2.5-years old sheep to elucidate the biology underlying differential impacts of late gestation versus early postnatal malnutrition on functional development of adipose tissues. Adipose tissues were obtained from 37 adult sheep born as twins to dams fed either NORM (fulfilling energy and protein requirements), LOW (50% of NORM) or HIGH (110% of protein and 150% of energy requirements) diets in the last 6-weeks of gestation. From day 3 to 6 months of age, lambs were fed high-carbohydrate-high-fat (HCHF) or moderate low-fat (CONV) diets, and thereafter the same moderate low-fat diet.

RESULTS

The gene expression profile of SUB in the adult sheep was not affected by the pre- or early postnatal nutrition history. In PER, 993 and 186 differentially expressed genes (DEGs) were identified in LOW versus HIGH and NORM, respectively, but no DEG was found between HIGH and NORM. DEGs identified in the mismatched pre- and postnatal nutrition groups LOW-HCHF (101) and HIGH-HCHF (192) were largely downregulated compared to NORM-CONV. Out of 831 DEGs, 595 and 236 were up- and downregulated in HCHF versus CONV, respectively. The functional enrichment analyses revealed that transmembrane (ion) transport activities, motor activities related to cytoskeletal and spermatozoa function (microtubules and the cytoskeletal motor protein, dynein), and responsiveness to the (micro) environmental extracellular conditions, including endocrine and nervous stimuli were enriched in the DEGs of LOW versus HIGH and NORM. We confirmed that mismatched pre- and postnatal feeding was associated with long-term programming of adipose tissue remodeling and immunity-related pathways. In agreement with phenotypic measurements, early postnatal HCHF feeding targeted pathways involved in kidney cell differentiation, and mismatched LOW-HCHF sheep had specific impairments in cholesterol metabolism pathways.

CONCLUSIONS

Both pre- and postnatal malnutrition differentially programmed (patho-) physiological pathways with implications for adipose functional development associated with metabolic dysfunctions, and PER was a major target.

Identification of Predictor Genes for Feed Efficiency in Beef Cattle by Applying Machine Learning Methods to Multi-Tissue Transcriptome Data

Machine learning (ML) methods have shown promising results in identifying genes when applied to large transcriptome datasets. However, no attempt has been made to compare the performance of combining different ML methods together in the prediction of high feed efficiency (HFE) and low feed efficiency (LFE) animals. In this study, using RNA sequencing data of five tissues (adrenal gland, hypothalamus, liver, skeletal muscle, and pituitary) from nine HFE and nine LFE Nellore bulls, we evaluated the prediction accuracies of five analytical methods in classifying FE animals. These included two conventional methods for differential gene expression (DGE) analysis (t-test and edgeR) as benchmarks, and three ML methods: Random Forests (RFs), Extreme Gradient Boosting (XGBoost), and combination of both RF and XGBoost (RX). Utility of a subset of candidate genes selected from each method for classification of FE animals was assessed by support vector machine (SVM). Among all methods, the smallest subsets of genes (117) identified by RX outperformed those chosen by t-test, edgeR, RF, or XGBoost in classification accuracy of animals. Gene co-expression network analysis confirmed the interactivity existing among these genes and their relevance within the network related to their prediction ranking based on ML. The results demonstrate a great potential for applying a combination of ML methods to large transcriptome datasets to identify biologically important genes for accurately classifying FE animals.

Transcriptomic analysis of the liver in aged laying hens with different intensity of brown eggshell color

OBJECTIVE

Eggshell color is an important indicator of egg quality for consumers, especially for brown eggs. Various factors related to laying hens and their environment affect brown eggshell coloration. However, there have been no studies investigating hepatic functions of laying hens with variable intensity of brown eggshell color. Therefore, this study was aimed to identify potential factors affecting brown eggshell coloration in aged laying hens at the hepatic transcriptomic level.

METHODS

Five hundred 92-wk-old Hy-line Brown laying hens were screened to select laying hens with different intensity of brown eggshell color based on eggshell color fans. Based on eggshell color scores, hens with dark brown eggshells (DBE; eggshell color fan score = 14.8) and hens with light brown eggshells (LBE; eggshell color fan score = 9.7) were finally selected for the liver sampling. We performed RNA-seq analysis using the liver samples through the paired-end sequencing libraries. Differentially expressed genes (DEGs) profiling was carried out to identify their biological meaning by bioinformatics.

RESULTS

A total of 290 DEGs were identified with 196 being up-regulated and 94 being down-regulated in DBE groups as compared to LBE groups. The Kyoto encyclopedia of genes and genomes (KEGG) analysis revealed that these DEGs belong to several biological pathways including herpes simplex infection (toll-like receptor 3 [TLR3], cyclin-dependent kinase 1, etc.) and influenza A (TLR3, radical S-adenosyl methionine domain containing 2, myxovirus [influenza virus] resistance 1, etc.). Genes related to stress response (ceremide kinase like) and nutrient metabolism (phosphoenolpyruvate carboxy-kinase 1, methylmalonic aciduria [cobalamin deficiency] cblB type, glycine receptor alpha 2, solute carrier family 7 member 11, etc.) were also identified to be differentially expressed.

CONCLUSION

The current results provide new insights regarding hepatic molecular functions related to different intensity of brown eggshell color in aged laying hens. These insights will contribute to future studies aiming to optimize brown eggshell coloration in aged laying hens.

Identification of hub genes involved in apparent metabolizable energy of chickens

Feed efficiency is one of the most economically significant traits in a breeding program. Apparent metabolizable energy is the most used method to evaluate energy utilization for feed efficiency. The purpose of this study was to identify candidate genes of chickens with divergent apparent metabolizable energy by bioinformatics analysis. The gene expression profile of duodenal of the highest and lowest apparent metabolizable energy-ranked birds were analyzed. Differentially expressed genes were picked out using GEO2R tool. Gene ontology and pathway analysis were performed using bioinformatics tools. Cytoscape software was used to visualize protein-protein network. There were 201 DEGs, including 99 up-regulated genes enriched in metabolic pathways, Cellular senescence and Focal adhesion, and 102 down-regulated genes enriched in metabolic pathways, Regulation of actin cytoskeleton, Neuroactive ligand-receptor interaction, Calcium signaling pathway and Focal adhesion. Two important modules were detected and pathway enrichment analysis showed that they were mainly associated with Focal adhesion, Regulation of actin cytoskeleton and RNA transport. Fifteen hub genes were selected and among them, ITGA8, CDC42 and GSK3B might be the core genes related to apparent metabolizable energy of chickens. These hub genes can be used as biomarkers for apparent metabolizable energy and feed efficiency in breeding program of chickens.