In-Depth Duodenal Transcriptome Survey in Chickens with Divergent Feed Efficiency Using RNA-Seq

Comparative Analysis of Runs of Homozygosity Islands in Indigenous and Commercial Chickens Revealed Candidate Loci for Disease Resistance and Production Traits

Runs of homozygosity (ROH) are contiguous stretches of identical genomic regions inherited from both parents. Assessment of ROH in livestock species contributes significantly to our understanding of genetic health, population genetic structure, selective pressure and conservation efforts. In this study, whole genome re-sequencing data from 140 birds of 10 Iranian indigenous chicken ecotypes, 3 commercial chicken breeds and 1 red junglefowl (RJF) population were used to investigate their population genetic structure, ROH characteristics (length and frequency) and genomic inbreeding coefficients (FROH). Additionally, we examined ROH islands for selection footprints in the indigenous chicken populations. Our results revealed distinct genetic backgrounds, among which the White Leghorn breed exhibited the greatest genetic distance from other populations, while the gamecock populations formed a separate cluster. We observed significant differences in ROH characteristics, in which the commercial breeds showed a higher number of ROH compared to indigenous chickens and red junglefowls. Short ROH ranging from 0.1 to 1 Mb were dominant among the populations. The Arian line had the highest mean length of ROH, while the White Leghorn breed showed the highest number of ROH. Among indigenous chickens, the Lari-Afghani ecotype exhibited the highest FROH, but the Sari inherited the richest genetic diversity. Interestingly, GGA16 carried no ROH in the red junglefowls, whereas GGA22 had the highest FROH across all populations, except in the Isfahan ecotype. We also identified ROH islands associated with genetic adaptations in indigenous ecotypes. These islands harboured immune-related genes contributing to disease resistance (TLR2, TICAM1, IL22RA1, NOS2, CCL20 and IFNLR1), heat tolerance and oxidative stress response (NFKB1, HSF4, OSGIN1 and BDNF), and muscle development, lipid metabolism and reproduction (MEOX2, CEBPB, CDS2 and GnRH-I). Overall, this study highlights the genetic potential of indigenous chickens to survive and adapt to their respective environments.

Novel insight into the feed conversion ratio in laying hens and construction of its prediction model

Feed efficiency (FE) is an important economic factor in poultry production, and feed conversion ratio (FCR) is one of the most widely used measures of FE. Factors associated with FCR include genetics, the environment, and other factors. However, the mechanisms responsible for FCR in chickens are still less well appreciated. In this study, we examined the pattern changes of FCR, then delved into understanding the mechanisms behind these variations from both genetic and environmental perspectives. Most interestingly, the FCR at the front section of henhouse exhibited the lowest value. Further investigation revealed that laying rate in the high FCR (HFCR) group was lower than that in the low FCR (LFCR) group (P < 0.05). Cortisol, total antioxidant capacity (TAOC), and IgG levels in the LFCR group were significantly lower than those in the HFCR group (P < 0.05), while BUN level was significantly higher than that in the HFCR group (P < 0.05). We identified a total of 67 and 10 differentially expressed genes (DEGs) associated with FCR in ovarian and small intestine tissues, respectively. Functional enrichment analysis of DEGs revealed that they might affect FCR by modulating genes associated with salivary secretion, ferroptosis, and mineral absorption. Moreover, values for relative humidity (RH), air velocity (AV), PM2.5, ammonia (NH3), and carbon dioxide (CO2) in the LFCR group were significantly lower than those in the HFCR group (P < 0.05). Conversely, value for light intensity (LI) in the LFCR group was significantly higher than that in the HFCR group (P < 0.05). Correlation analysis revealed a positive correlation between FCR and RH, AV, PM2.5, NH3, and CO2, and a negative correlation with LI. Finally, the FCR prediction model was successfully constructed based on multiple environmental variables using the random forest algorithm, providing a valuable tool for predicting FCR in chickens.

Multi-tissue transcriptome profiling linked the association between tissue-specific circRNAs and the heterosis for feed intake and efficiency in chicken

Heterosis has been widely utilized in chickens. The nonadditive inheritance of genes contributes to this biological phenomenon. However, the role of circRNAs played in the heterosis is poorly determined. In this study, we observed divergent heterosis for residual feed intake (RFI) between 2 crossbreds derived from a reciprocal cross between White Leghorns and Beijing You chickens. Then, circRNA landscape for 120 samples covering the hypothalamus, liver, duodenum mucosa and ovary were profiled to elucidate the regulatory mechanisms of heterosis. We detected that a small proportion of circRNAs (7.83-20.35%) were additively and non-additively expressed, in which non-additivity was a major inheritance of circRNAs in the crossbreds. Tissue-specific expression of circRNAs was prevalent across 4 tissues. Weighted gene co-expression network analysis revealed circRNA-mRNA co-expression modules associated with feed intake and RFI in the hypothalamus and liver, and the co-expressed genes were enriched in oxidative phosphorylation pathway. We further identified 8 nonadditive circRNAs highly correlated with 16 nonadditive genes regulating negative heterosis for RFI in the 2 tissues. Circ-ITSN2 was validated in the liver tissue for its significantly positive correlation with PGPEP1L. Moreover, the bioinformatic analysis indicated that candidate circRNAs might be functioned by binding the microRNAs and interacting with the RNA binding proteins. The integration of multi-tissue transcriptome firstly linked the association between tissue-specific circRNAs and the heterosis for feed intake and efficiency in chicken, which provide novel insights into the molecular mechanism underlying heterosis for feed efficiency. The validated circRNAs can act as potential biomarkers for predicting RFI and its heterosis.

The Omics Revolution in Understanding Chicken Reproduction: A Comprehensive Review

Omics approaches have significantly contributed to our understanding of several aspects of chicken reproduction. This review paper gives an overview of the use of omics technologies such as genomics, transcriptomics, proteomics, and metabolomics to elucidate the mechanisms of chicken reproduction. Genomics has transformed the study of chicken reproduction by allowing the examination of the full genetic makeup of chickens, resulting in the discovery of genes associated with reproductive features and disorders. Transcriptomics has provided insights into the gene expression patterns and regulatory mechanisms involved in reproductive processes, allowing for a better knowledge of developmental stages and hormone regulation. Furthermore, proteomics has made it easier to identify and quantify the proteins involved in reproductive physiology to better understand the molecular mechanisms driving fertility, embryonic development, and egg quality. Metabolomics has emerged as a useful technique for understanding the metabolic pathways and biomarkers linked to reproductive performance, providing vital insights for enhancing breeding tactics and reproductive health. The integration of omics data has resulted in the identification of critical molecular pathways and biomarkers linked with chicken reproductive features, providing the opportunity for targeted genetic selection and improved reproductive management approaches. Furthermore, omics technologies have helped to create biomarkers for fertility and embryonic viability, providing the poultry sector with tools for effective breeding and reproductive health management. Finally, omics technologies have greatly improved our understanding of chicken reproduction by revealing the molecular complexities that underpin reproductive processes.

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

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

Analysis of liver miRNA in Hu sheep with different residual feed intake

Feed efficiency (FE), an important economic trait in sheep production, is indirectly assessed by residual feed intake (RFI). However, RFI in sheep is varied, and the molecular processes that regulate RFI are unclear. It is thus vital to investigate the molecular mechanism of RFI to developing a feed-efficient sheep. The miRNA-sequencing (RNA-Seq) was utilized to investigate miRNAs in liver tissue of 6 out of 137 sheep with extreme RFI phenotypic values. In these animals, as a typical metric of FE, RFI was used to distinguish differentially expressed miRNAs (DE_miRNAs) between animals with high (n = 3) and low (n = 3) phenotypic values. A total of 247 miRNAs were discovered in sheep, with four differentially expressed miRNAs (DE_miRNAs) detected. Among these DE_miRNAs, three were found to be upregulated and one was downregulated in animals with low residual feed intake (Low_RFI) compared to those with high residual feed intake (High_RFI). The target genes of DE_miRNAs were primarily associated with metabolic processes and biosynthetic process regulation. Furthermore, they were also considerably enriched in the FE related to glycolysis, protein synthesis and degradation, and amino acid biosynthesis pathways. Six genes were identified by co-expression analysis of DE_miRNAs target with DE_mRNAs. These results provide a theoretical basis for us to understand the sheep liver miRNAs in RFI molecular regulation.

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.

Comparative analysis of the characteristics of digestive organs in broiler chickens with different feed efficiencies

Improving feed efficiency is one of the main goals of chicken breeding and production. The function of the digestive system, where feed is digested and nutrients are absorbed, is closely related to feed efficiency. However, the association between feed efficiency and the development of different digestive organs in chickens remains unclear. Here, we investigated the individual feed efficiency of 207 broilers during the fast-growing period with an electronic feeder and examined the characteristics of 8 organs of their digestive system (the liver, bile, proventriculus, gizzard, duodenum, jejunum, ileum, and cecum) at market age. Both the feed conversion ratio (FCR) and residual feed intake (RFI) were significantly negatively correlated with the gizzard weight (GW) and significantly positively correlated with the relative weight of the liver (RLW). Additionally, we found an obvious negative relationship between the FCR and cecal length (CL). A two-tailed t test further confirmed these correlation analysis results. Specifically, compared to birds with the lowest feed efficiencies, the GW of broilers with the highest feed efficiencies (the lowest FCR or RFI) was 22.74% and 17.97% higher, respectively. The RLW of chickens with the highest feed efficiencies was 10.82 to 13.73% less than that of chickens with the lowest feed efficiencies. In addition, we found that increased CL (5.42–12.09%) was significantly associated with better feed efficiency. Thus, our study showed that the feed efficiency of broilers was related to the development of the gizzard, liver, and cecum. These findings provide new insight into the genetic and physiological regulation of feed efficiency in broilers.

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.

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.

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.

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.

Transcriptome Analysis Identifies Candidate Genes and Signaling Pathways Associated With Feed Efficiency in Xiayan Chicken

Feed efficiency is an important economic factor in poultry production, and the rate of feed efficiency is generally evaluated using residual feed intake (RFI). The molecular regulatory mechanisms of RFI remain unknown. Therefore, the objective of this study was to identify candidate genes and signaling pathways related to RFI using RNA-sequencing for low RFI (LRFI) and high RFI (HRFI) in the Xiayan chicken, a native chicken of the Guangxi province. Chickens were divided into four groups based on FE and sex: LRFI and HRFI for males and females, respectively. We identified a total of 1,015 and 742 differentially expressed genes associated with RFI in males and females, respectively. The 32 and 7 Gene Ontology (GO) enrichment terms, respectively, identified in males and females chiefly involved carbohydrate, amino acid, and energy metabolism. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified 11 and 5 significantly enriched signaling pathways, including those for nutrient metabolism, insulin signaling, and MAPK signaling, respectively. Protein-protein interaction (PPI) network analysis showed that the pathways involving CAT, ACSL1, ECI2, ABCD2, ACOX1, PCK1, HSPA2, and HSP90AA1 may have an effect on feed efficiency, and these genes are mainly involved in the biological processes of fat metabolism and heat stress. Gene set enrichment analysis indicated that the increased expression of genes in LRFI chickens was related to intestinal microvilli structure and function, and to the fat metabolism process in males. In females, the highly expressed set of genes in the LRFI group was primarily associated with nervous system and cell development. Our findings provide further insight into RFI regulation mechanisms in chickens.

RNA seq analyses of chicken reveals biological pathways involved in acclimation into different geographical locations

Transcriptome expression reflects genetic response in diverse conditions. In this study, RNA sequencing was utilized to profile multiple tissues such as liver, breast, caecum, and gizzard of Korean commercial chicken raised in Korea and Kyrgyzstan. We analyzed ten samples per tissue from each location to identify candidate genes which are involved in the adaptation of Korean commercial chicken to Kyrgyzstan. At false discovery rate (FDR) < 0.05 and fold change (FC) > 2, we found 315, 196, 167 and 198 genes in liver, breast, cecum, and gizzard respectively as differentially expressed between the two locations. GO enrichment analysis showed that these genes were highly enriched for cellular and metabolic processes, catalytic activity, and biological regulations. Similarly, KEGG pathways analysis indicated metabolic, PPAR signaling, FoxO, glycolysis/gluconeogenesis, biosynthesis, MAPK signaling, CAMs, citrate cycles pathways were differentially enriched. Enriched genes like TSKU, VTG1, SGK, CDK2 etc. in these pathways might be involved in acclimation of organisms into diverse climatic conditions. The qRT-PCR result also corroborated the RNA-Seq findings with R² of 0.76, 0.80, 0.81, and 0.93 for liver, breast, caecum, and gizzard respectively. Our findings can improve the understanding of environmental acclimation process in chicken.

Estimation of the genetic parameters of traits relevant to feed efficiency: result from broiler lines divergent for high or low abdominal fat content

Feed consumption represents a major cost in poultry production and improving feed efficiency is one of the important goals in breeding strategies. The present study aimed to analyze the relationship between feed efficiency and relevant traits and find the proper selection method for improving feed efficiency by using the Northeast Agricultural University High and Low Fat broiler lines that were divergently selected for abdominal fat content. A total of 899 birds were used to measure the feed intake (FI), abdominal fat weight (AFW), and body weight traits. The abdominal fat percentage (AFP), feed conversion ratio (FCR), and the residual feed intake (RFI) were calculated for each individual broiler. The differences in the AFW, AFP, and in traits relevant to feed efficiency, such as FCR and RFI, between the fat line and the lean line were analyzed, and the genetic parameters were estimated for AFW, AFP, and feed efficiency relevant traits. The results showed that AFW, AFP, body weight gain (BWG), FI, FCR, and RFI were significantly higher in the fat line compared with the lean line. The heritability of FI, BWG, FCR, RFI, AFW, and AFP were 0.45, 0.28, 0.36, 0.38, 0.33, and 0.30, respectively. Both FCR and RFI showed high positive genetic correlations with FI, AFW, and AFP and relatively low, negative genetic correlations with BWG. The RFI showed much higher positive genetic correlation with the abdominal fat traits than FCR. In addition, the FCR showed negative genetic correlation with body weight of 4 wk (BW4) and 7 wk (BW7), whereas RFI showed positive genetic correlation with BW4 and BW7. The results showed that both RFI and FCR could be used for improving feed efficiency. When selecting against RFI, the AFP could be significantly reduced, and by selecting against FCR, the body weight could be improved simultaneously.

Quercetin Improving Lipid Metabolism by Regulating Lipid Metabolism Pathway of Ileum Mucosa in Broilers

This study is aimed at evaluating the regulatory mechanism of quercetin on lipid metabolism in the ileum of broilers to better understand these pathways decreasing abdominal fat. 480 chickens were randomly divided into 4 groups (control, 0.02% quercetin, 0.04% quercetin, and 0.06% quercetin). Breast muscle, thigh muscle, and abdominal fat pad were removed and weighed at 42 d of age. Serum was obtained by centrifuging blood samples from the jugular vein (10 ml) to determine high-density lipoprotein (HDL), total cholesterol (TC), low-density lipoprotein (LDL), triglyceride (TG), leptin, and adiponectin using ELISA. About 5 g of the ileum was harvested and immediately frozen in liquid nitrogen for RNA-seq. Then, the confirmation of RNA-seq results by the Real-Time Quantitative PCR (RT-qPCR) method was evaluated using Pearson's correlation. Compared with control, abdominal fat percentage was significantly decreased with increasing quercetin supplementation, and the best result was obtained at 0.06% dietary quercetin supplementation (P < 0.01). Breast muscle percentage was significantly decreased at 0.02% quercetin (P < 0.01), and thigh muscle percentage tended to increase (P = 0.078). Meanwhile, 0.04% and 0.06% quercetin significantly decreased TG (P < 0.01), TC (P < 0.01), and LDL content (P < 0.05) in serum. Serum leptin and adiponectin contents were significantly increased by 0.04% and 0.06% dietary quercetin supplementation, compared with the control (P < 0.01). Analyses of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database were used to identify differently expressed genes and lipid metabolism pathways. Quercetin decreased abdominal fat percentage through regulating fat digestion and absorption, glycerophospholipid metabolism, AMPK signaling pathway, fatty acid degradation, and cholesterol metabolism.

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 re-sequencing and transcriptome analysis reveal candidate genes associated with the pendulous comb phenotype in domestic chickens

To determine the causative variations associated with two chicken comb phenotypes, pendulous comb (PC) or upright comb (UC), two pooled genomic DNA samples from PC and UC chickens were re-sequenced by Next-Generation Sequencer, and genome-wide Single nucleotide polymorphisms (SNPs) were detected. Using three selective sweep approaches, F_ST , θπ, and Tajima's D, with top 5% window values serving as the threshold, a total of 84 positively selective genes (PSGs) were identified. There were no SNPs in exons of the PSGs with significant differences in allele frequencies between the two comb phenotype groups. Then, 515 differentially expressed genes (DEGs) between the PC and UC were identified by RNA-seq. Three genes including CD36 (CD36 molecule), ADAMTSL3 (ADAMTS-like 3), and AOX1 (aldehyde oxidases 1) are overlapped between PSGs and DEGs. After genotyping seven candidate SNPs in the regulatory regions of the three overlapping genes in 120 chickens from two other breeds, two variants (rs14607046 and rs731818051) in the regulatory regions of AOX1 and ADAMTSL3 were found to have significant differences in allele frequency between the PC and UC, suggesting that the two variants may be causative mutations for PC. Overall, our study shed light on the genetic basis underlying the PC phenotype in chickens.

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

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

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

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

Expression of genes related to lipid transport in meat-type ducks divergent for low or high residual feed intake

Objective

This study examined the effects of divergence in residual feed intake (RFI) on expression profiles of key genes related to lipid transport in the liver and duodenal epithelium and their associations with feed efficiency traits in meat-type ducks.

Methods

A total of 1,000 male ducks with similar body weight (1,042.1±87.2 g) were used in this study, and their individual RFI was calculated from 21 to 42 d of age. Finally, the 10 highest RFI (HRFI) and 10 lowest RFI (LRFI) ducks were chosen for examining the expression of key genes related to lipid transport in the liver and duodenal epithelium using quantitative polymerase chain reaction.

Results

In the liver, expression levels of albumin (ALB), CD36 molecule (CD36), fatty acid hydroxylase domain containing 2 (FAXDC2), and choline kinase alpha (CHKA) were significantly higher in LRFI ducks than in HRFI ducks (p<0.01); negative correlations (p<0.05) between expression levels of ALB, CD36, FAXDC2, and CHKA and RFI were detected in the liver. Additionally, ALB expression was strongly positively correlated (p<0.05) with CD36, FAXDC2, CHKA, and apolipoprotein H (APOH) expression in the liver. In duodenal epithelium, we found that mRNA levels of ALB, CD36, FAXDC2, and APOH were significantly higher in LRFI ducks than in HRFI ducks (p<0.01); RFI was strongly negatively correlated (p<0.05) with ALB, FAXDC2, and APOH expression, while ALB expression was strongly positively correlated with APOH expression (p<0.01) in duodenal epithelium. Furthermore, expression levels of both ALB and FAXDC2 genes were significantly associated with feed conversion ratio and RFI in both liver and duodenal epithelium (p<0.05).

Conclusion

Our findings therefore suggest that ALB and FAXDC2 genes might be used as potential gene markers designed to improve feed efficiency in future meat-type duck breeding programs.

A Nutrigenomics Approach Using RNA Sequencing Technology to Study Nutrient-Gene Interactions in Agricultural Animals

Thorough understanding of animal gene expression driven by dietary nutrients can be regarded as a bottom line of advanced animal nutrition research. Nutrigenomics (including transcriptomics) studies the effects of dietary nutrients on cellular gene expression and, ultimately, phenotypic changes in living organisms. Transcriptomics can be applied to investigate animal tissue transcriptomes at a defined nutritional state, which can provide a holistic view of intracellular RNA expression. As a novel transcriptomics approach, RNA sequencing (RNA-Seq) technology can monitor all gene expressions simultaneously in response to dietary intervention. The principle and history of RNA-Seq are briefly reviewed, and its 3 principal steps are described in this article. Application of RNA-Seq in different areas of animal nutrition research is summarized. Lastly, the application of RNA-Seq in swine science and nutrition is also reviewed. In short, RNA-Seq holds significant potential to be employed for better understanding the nutrient-gene interactions in agricultural animals.

Feed efficiency measures and their relationships with production and meat quality traits in slower growing broilers

Feed consumption accounts for the major cost of broiler production. Improving the efficiency of feed utilization is a primary goal in breeding strategies, although few studies have focused on slower growing broilers. Here, we recorded the feed intake (FI) during the fast-growing period (d 56 to 76) and measured the live weight, body measurements, carcass characteristics, and intramuscular fat (IMF) content of Chinese yellow broilers. Then, the residual feed intake (RFI) and feed conversion ratio (FCR) were calculated for each individual. Pair-wise phenotypic correlations were subsequently calculated between feed efficiency traits and others. Finally, we separately selected the more efficient individuals based on RFI and FCR values to evaluate the impacts on the traits of FI, growth, carcass characteristics, and meat quality. The results showed higher correlations between FCR and production traits than with RFI, while RFI showed a moderate and positive phenotypic correlation with abdominal fat. FCR was weakly correlated with FI and slightly positively correlated with IMF content. The correlation coefficient between RFI and FI was 0.62, and that between RFI and IMF content was close to zero. Without increasing FI, decreasing FCR could effectively enhance the growth rate and market weight with no adverse effect on meat quality. In contrast, by improving RFI, FI and abdominal fat mass were significantly reduced and thus increased the yield with no unfavorable effects on meat quality. In consideration of consumer preference and overall economical benefits, RFI is a more suitable index to improve feed efficiency in slower growing broilers.

Transcriptomic Profiles of Brain Provide Insights into Molecular Mechanism of Feed Conversion Efficiency in Crucian Carp (Carassius auratus)

Feed efficiency is an economically crucial trait for cultured animals, however, progress has been scarcely made in the genetic analyses of feed conversion efficiency (FCE) in fish because of the difficulties in measurement of trait phenotypes. In the present investigation, we present the first application of RNA sequencing (RNA-Seq) combined with differentially expressed genes (DEGs) analysis for identification of functional determinants related to FCE at the gene level in an aquaculture fish, crucian carp (Carassius auratus). Brain tissues of six crucian carp with extreme FCE performances were subjected to transcriptome analysis. A total of 544,612 unigenes with a mean size of 644.38 bp were obtained from Low- and High-FCE groups, and 246 DEGs that may be involved in FCE traits were identified in these two groups. qPCR confirmed that genes previously identified as up- or down-regulated by RNA-Seq were effectively up- or down-regulated under the studied conditions. Thirteen key genes, whose functions are associated with metabolism (Dgkk, Mgst3 and Guk1b), signal transduction (Vdnccsa1b, Tgfα, Nr4a1 and Tacr2) and growth (Endog, Crebrtc2, Myh7, Myh1,Myh14 and Igfbp7) were identified according to GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) annotations. Our novel findings provide useful pathway information and candidate genes for future studies of genetic mechanisms underlying FCE in crucian carp.

Identification of molecular pathways and candidate genes associated with cocks' comb size trait by genome-wide transcriptome analysis

The comb of the male is an important secondary sexual characteristic. Although quantitative trait loci (QTLs) related to comb size have been identified, molecular mechanisms underlying this trait remain mostly unknown. In this study, RNA sequencing (RNA-seq) was employed to compare whole transcriptomic differences between two groups of Partridge Shank chickens that are divergent in comb sizes. A total of 563 differentially expressed genes (DEGs) were identified, including 277 up-regulated and 286 down-regulated DEGs. According to the animal QTL database, eight DEGs including BMP2 and CHADL matching the reported QTLs were associated with the comb size. Functional annotation analysis revealed that DEGs were involved in cell communication and calcium signaling. Protein-protein interaction network analysis showed that STK32A, PIK3R1, EDN1, HSPA5, and HSPA8 have an impact on comb growth. Moreover, potential alternative splicing events and single nucleotide polymorphisms were also identified. Our data provide a source for identifying genes and pathways with functions critical to comb size and accelerate studies involving molecular mechanisms of this sexual ornament.

Exploring Genomic Variants Related to Residual Feed Intake in Local and Commercial Chickens by Whole Genomic Resequencing

Improving feed efficiency is a major goal in poultry production to reduce production costs and increase profitability. The genomic variants and possible molecular mechanisms responsible for residual feed intake (RFI) in chickens, however, remain poorly understood. In this study, using both local and commercial breeds, genome re-sequencing of low RFI and high RFI chickens was performed to elucidate the genomic variants underlying RFI. Results showed that 8,505,214 and 8,479,041 single nucleotide polymorphisms (SNPs) were detected in low and high RFI Beijing-You chickens, respectively; 8,352,008 and 8,372,769 SNPs were detected in low- and high-RFI Cobb chickens, respectively. Through a series of filtering processes, 3746 candidate SNPs assigned to 1137 genes in Beijing-You chickens and 575 candidate SNPs (448 genes) in Cobb chickens were found. The validation of the selected 191 SNPs showed that 46 SNPs were significantly associated with the RFI in an independent population of 779 Cobb chickens, suggesting that the method of screening associated SNPs with whole genome sequencing (WGS) strategy was reasonable. Functions annotation of RFI-related genes indicated that genes in Beijing-You were enriched in lipid and carbohydrate metabolism, as well as the phosphatase and tensin homolog (PTEN) signaling pathway. In Cobb, however, RFI-related genes were enriched in the feed behavior process and cAMP responsive element binding protein (CREB) signaling pathway. For both breeds, organismal development physiological processes were enriched. Correspondingly, NOS1, PHKG1, NEU3 and PIP5K1B were differentially expressed in Beijing-You, while CDC42, CSK, PIK3R3, CAMK4 and PLCB4 were differentially expressed in Cobb, suggesting that these might be key genes that contribute to RFI. The results of the present study identified numerous novel SNPs for RFI, which provide candidate biomarkers for use in the genetic selection for RFI. The study has improved knowledge of the genomic variants and potential biological pathways underlying RFI in chickens.

Thyroid transcriptome analysis reveals different adaptive responses to cold environmental conditions between two chicken breeds

Selection for cold tolerance in chickens is important for improving production performance and animal welfare. The identification of chicken breeds with higher cold tolerance and production performance will help to target candidates for the selection. The thyroid gland plays important roles in thermal adaptation, and its function is influenced by breed differences and transcriptional plasticity, both of which remain largely unknown in the chicken thyroid transcriptome. In this study, we subjected Bashang Long-tail (BS) and Rhode Island Red (RIR) chickens to either cold or warm environments for 21 weeks and investigated egg production performance, body weight changes, serum thyroid hormone concentrations, and thyroid gland transcriptome profiles. RIR chickens had higher egg production than BS chickens under warm conditions, but BS chickens produced more eggs than RIRs under cold conditions. Furthermore, BS chickens showed stable body weight gain under cold conditions while RIRs did not. These results suggested that BS breed is a preferable candidate for cold-tolerance selection and that the cold adaptability of RIRs should be improved in the future. BS chickens had higher serum thyroid hormone concentrations than RIRs under both environments. RNA-Seq generated 344.3 million paired-end reads from 16 sequencing libraries, and about 90% of the processed reads were concordantly mapped to the chicken reference genome. Differential expression analysis identified 46-1,211 genes in the respective comparisons. With regard to breed differences in the thyroid transcriptome, BS chickens showed higher cell replication and development, and immune response-related activity, while RIR chickens showed higher carbohydrate and protein metabolism activity. The cold environment reduced breed differences in the thyroid transcriptome compared with the warm environment. Transcriptional plasticity analysis revealed different adaptive responses in BS and RIR chickens to cope with the cold, and showed higher responsiveness in BS compared with RIR chickens, suggesting greater adaptability of the thyroid in BS chickens. Moreover, 10,053 differential splicing events were revealed among the groups, with RNA splicing and processing, gene expression, transport, and metabolism being the main affected biological processes, identifying a valuable alternative splicing repertoire for the chicken thyroid. A short isoform of TPO (encoding thyroid peroxidase) containing multiple open reading frames was generated in both breeds by skipping exons 4 and 5 in the cold environment. These findings provide novel clues for future studies of the molecular mechanisms underlying cold adaptation and/or acclimation in chickens.

Quantitative trait loci mapping for feed conversion efficiency in crucian carp (Carassius auratus)

QTL is a chromosomal region including single gene or gene clusters that determine a quantitative trait. While feed efficiency is highly important in aquaculture fish, little genetic and genomic progresses have been made for this trait. In this study, we constructed a high-resolution genetic linkage map in a full-sib F1 family of crucian carp (Carassius auratus) consisting of 113 progenies with 8,460 SNP markers assigning onto 50 linkage groups (LGs). This genetic map spanned 4,047.824 cM (0.478 cM/marker) and covered 98.76% of the crucian carp genome. 35 chromosome-wide QTL affecting feed conversion efficiency (FCE, 8 QTL), relative growth rate (RGR, 9 QTL), average daily gain (ADG, 13 QTL) and average daily feed intake (ADFI, 5 QTL) were detected on 14 LGs, explaining 14.0–20.9% of the phenotypic variations. In LGs of LG16, LG25, LG36 and LG49, several QTL affecting different traits clustered together at the identical or close regions of the same linkage group. Seven candidate genes, whose biological functions may involve in the energy metabolism, digestion, biosynthesis and signal transduction, were identified from these QTL intervals by comparative genomics analysis. These results provide a basis for elucidating genetic mechanism of feed efficiency and potential marker-assisted selection in crucian carp.

Gene expression profiling reveals candidate genes related to residual feed intake in duodenum of laying ducks

Feed represents two-thirds of the total costs of poultry production, especially in developing countries. Improvement in feed efficiency would reduce the amount of feed required for production (growth or laying), the production cost, and the amount of nitrogenous waste. The most commonly used measures for feed efficiency are feed conversion ratio (FCR) and residual feed intake (RFI). As a more suitable indicator assessing feed efficiency, RFI is defined as the difference between observed and expected feed intake based on maintenance and growth or laying. However, the genetic and biological mechanisms regulating RFI are largely unknown. Identifying molecular mechanisms explaining divergence in RFI in laying ducks would lead to the development of early detection methods for the selection of more efficient breeding poultry. The objective of this study was to identify duodenum genes and pathways through transcriptional profiling in 2 extreme RFI phenotypes (HRFI and LRFI) of the duck population. Phenotypic aspects of feed efficiency showed that RFI was strongly positive with FCR and feed intake (FI). Transcriptomic analysis identified 35 differentially expressed genes between LRFI and HRFI ducks. These genes play an important role in metabolism, digestibility, secretion, and innate immunity including (), (), (), β (), and (). These results improve our knowledge of the biological basis underlying RFI, which would be useful for further investigations of key candidate genes for RFI and for the development of biomarkers.

RNA-Seq analysis on chicken taste sensory organs: An ideal system to study organogenesis

RNA-Seq is a powerful tool in transcriptomic profiling of cells and tissues. We recently identified many more taste buds than previously appreciated in chickens using molecular markers to stain oral epithelial sheets of the palate, base of oral cavity, and posterior tongue. In this study, RNA-Seq was performed to understand the transcriptomic architecture of chicken gustatory tissues. Interestingly, taste sensation related genes and many more differentially expressed genes (DEGs) were found between the epithelium and mesenchyme in the base of oral cavity as compared to the palate and posterior tongue. Further RNA-Seq using specifically defined tissues of the base of oral cavity demonstrated that DEGs between gustatory (GE) and non-gustatory epithelium (NGE), and between GE and the underlying mesenchyme (GM) were enriched in multiple GO terms and KEGG pathways, including many biological processes. Well-known genes for taste sensation were highly expressed in the GE. Moreover, genes of signaling components important in organogenesis (Wnt, TGFβ/ BMP, FGF, Notch, SHH, Erbb) were differentially expressed between GE and GM. Combined with other features of chicken taste buds, e.g., uniquely patterned array and short turnover cycle, our data suggest that chicken gustatory tissue provides an ideal system for multidisciplinary studies, including organogenesis and regenerative medicine.

Rapid, broad-scale gene expression evolution in experimentally harvested fish populations

Gene expression changes potentially play an important role in adaptive evolution under human-induced selection pressures, but this has been challenging to demonstrate in natural populations. Fishing exhibits strong selection pressure against large body size, thus potentially inducing evolutionary changes in life history and other traits that may be slowly reversible once fishing ceases. However, there is a lack of convincing examples regarding the speed and magnitude of fisheries-induced evolution, and thus, the relevant underlying molecular-level effects remain elusive. We use wild-origin zebrafish (Danio rerio) as a model for harvest-induced evolution. We experimentally demonstrate broad-scale gene expression changes induced by just five generations of size-selective harvesting, and limited genetic convergence following the cessation of harvesting. We also demonstrate significant allele frequency changes in genes that were differentially expressed after five generations of size-selective harvesting. We further show that nine generations of captive breeding induced substantial gene expression changes in control stocks likely due to inadvertent selection in the captive environment. The large extent and rapid pace of the gene expression changes caused by both harvest-induced selection and captive breeding emphasizes the need for evolutionary enlightened management towards sustainable fisheries.

RNA-Seq transcriptomics and pathway analyses reveal potential regulatory genes and molecular mechanisms in high- and low-residual feed intake in Nordic dairy cattle

Background

The selective breeding of cattle with high-feed efficiencies (FE) is an important goal of beef and dairy cattle producers. Global gene expression patterns in relevant tissues can be used to study the functions of genes that are potentially involved in regulating FE. In the present study, high-throughput RNA sequencing data of liver biopsies from 19 dairy cows were used to identify differentially expressed genes (DEGs) between high- and low-FE groups of cows (based on Residual Feed Intake or RFI). Subsequently, a profile of the pathways connecting the DEGs to FE was generated, and a list of candidate genes and biomarkers was derived for their potential inclusion in breeding programmes to improve FE.

Results

The bovine RNA-Seq gene expression data from the liver was analysed to identify DEGs and, subsequently, identify the molecular mechanisms, pathways and possible candidate biomarkers of feed efficiency. On average, 57 million reads (short reads or short mRNA sequences < ~200 bases) were sequenced, 52 million reads were mapped, and 24,616 known transcripts were quantified according to the bovine reference genome. A comparison of the high- and low-RFI groups revealed 70 and 19 significantly DEGs in Holstein and Jersey cows, respectively. The interaction analysis (high vs. low RFI x control vs. high concentrate diet) showed no interaction effects in the Holstein cows, while two genes showed interaction effects in the Jersey cows. The analyses showed that DEGs act through certain pathways to affect or regulate FE, including steroid hormone biosynthesis, retinol metabolism, starch and sucrose metabolism, ether lipid metabolism, arachidonic acid metabolism and drug metabolism cytochrome P450.

Conclusion

We used RNA-Seq-based liver transcriptomic profiling of high- and low-RFI dairy cows in two breeds and identified significantly DEGs, their molecular mechanisms, their interactions with other genes and functional enrichments of different molecular pathways. The DEGs that were identified were the CYP’s and GIMAP genes for the Holstein and Jersey cows, respectively, which are related to the primary immunodeficiency pathway and play a major role in feed utilization and the metabolism of lipids, sugars and proteins.

Electronic supplementary material

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