Expression of candidate genes for residual feed intake in Angus cattle

Scans for Signatures of Selection in Genomes of Wagyu and Buryat Cattle Breeds Reveal Candidate Genes and Genetic Variants for Adaptive Phenotypes and Production Traits

Past and ongoing selection shapes the genomes of livestock breeds. Identifying such signatures of selection allows for uncovering the genetic bases of affected phenotypes, including economically important traits and environmental adaptations, for the further improvement of breed genetics to respond to climate and economic challenges. Turano-Mongolian cattle are a group of taurine breeds known for their adaptation to extreme environmental conditions and outstanding production performance. Buryat Turano-Mongolian cattle are among the few breeds adapted to cold climates and poor forage. Wagyu, on the other hand, is famous for high productivity and unique top-quality marbled meat. We used hapFLK, the de-correlated composite of multiple signals (DCMS), PBS, and FST methods to search for signatures of selection in their genomes. The scans revealed signals in genes related to cold adaptation (e.g., STAT3, DOCK5, GSTM3, and CXCL8) and food digestibility (SI) in the Buryat breed, and growth and development traits (e.g., RBFOX2 and SHOX2) and marbling (e.g., DGAT1, IQGAP2, RSRC1, and DIP2B) in Wagyu. Several putatively selected genes associated with reproduction, immunity, and resistance to pathogens were found in both breed genomes. The results of our work could be used for creating new productive adapted breeds or improving the extant breeds.

Identification of SNPs and expression patterns of ALB, AHSG and GC genes and their association with growth traits in Hu sheep

Growth traits are economically important traits in sheep breeding. This study was conducted to evaluate the polymorphisms of ALB, AHSG and GC genes and their association with growth traits in Hu sheep. We measured and recorded the body weight (BW), body height (BH), body length (BL) and feed conversion ratio (FCR) of 1418 male Hu sheep raised in the same environment from 80 to 180 days of age. The total of four SNPs in the ALB, AHSG and GC genes were identified by direct sequencing technology. The results of association analysis showed that two loci (g.8699 A>T and g.13458 T>C) of ALB gene significantly affect average daily gain (ADG; P < 0.05). The genotypes of SNP g.2454 T>C in AHSG gene were significantly associated with ADG and FCR (P < 0.05). There were significant associations between GC g.19484 A>C and BW, BH and BL (P < 0.05). The results of qRT-PCR showed that ALB, AHSG, and GC genes were extremely significantly higher in H_BW sheep compared with those in the L_BW sheep (P < 0.01). These results revealed that ALB-1 g.8699 A>T, ALB-2 g.13458 T>C, AHSG g.2454 T>C and GC g.19484 A>C loci are potential molecular markers for Hu sheep breeding.

Effects of Feeding and Drinking Behavior on Performance and Carcass Traits in Beef Cattle

Feed and water efficiency are important traits to improve beef cattle production’s economic and environmental sustainability. This study evaluated residual feed intake (RFI) and residual water intake (RWI) and their relationship with performance, ingestive behavior, and carcass traits in Caracu beef cattle. The data were analyzed using a generalized linear model with least squares means. The ingestive behavior, performance, and carcass traits were influenced by sex (p < 0.05). Males showed higher dry matter intake (DMI), average daily gain (ADG), mid-test metabolic weight (BW0.75), rib eye area, and rump fat thickness than females, besides spending more time drinking and eating. Low RFI animals exhibited higher DMI than high RFI animals. Low RWI animals ingested 3.89 L/d of water further than high RWI animals. The interaction between sex and RWI influenced the DMI, BW0.75, and backfat thickness. The ingestive behavior of low and high RFI animals was similar, although high RWI animals visited a smaller number of drinkers than low RWI animals. Water intake positively affects productive efficiency, and the combined use of RWI and RFI may help improve the selection of more efficient animals contributing to reducing the costs of beef cattle production and improving environmental sustainability.

Proteomic and parallel reaction monitoring approaches to evaluate biomarkers of mutton tenderness

Intensive fattening usually results in the changes of meat quality. Tenderness is a central attribute for mutton sensory qualities and consumers' choice. Here, we reported that intensive fattening mutton was more tender than that of traditionally raised sheep. By proteomic approach, we found 49 differentially expressed proteins in longissimus dorsi muscle. After bioinformatics analysis, 5 cytoskeletal proteins, 3 protein binding proteins and 7 metabolic enzymes were identified as potential biomarkers for mutton tenderness. Finally, we verified the expression of these abundant proteins by parallel reaction monitoring (PRM). Collectively, our results reveal that the mutton of sheep raised by intensive fattening is more tender than that of traditionally raised sheep. Myosin-2, myosin-13, vimentin, carbonic anhydrase, carbonic anhydrase-2, Glutathione S-transferase and Microtubule-associated protein 4 isoform X1 can be candidate biomarkers for mutton tenderness. Our data also indicate a central role of cytoskeletal proteins and metabolic enzymes in determining mutton tenderness.

Physiological responses and adaptations to high methane production in Japanese Black cattle

In this study, using enteric methane emissions, we investigated the metabolic characteristics of Japanese Black cattle. Their methane emissions were measured at early (age 13 months), middle (20 months), and late fattening phases (28 months). Cattle with the highest and lowest methane emissions were selected based on the residual methane emission values, and their liver transcriptome, blood metabolites, hormones, and rumen fermentation characteristics were analyzed. Blood β-hydroxybutyric acid and insulin levels were high, whereas blood amino acid levels were low in cattle with high methane emissions. Further, propionate and butyrate levels differed depending on the enteric methane emissions. Hepatic genes, such as SERPINI2, SLC7A5, ATP6, and RRAD, which were related to amino acid transport and glucose metabolism, were upregulated or downregulated during the late fattening phase. The above mentioned metabolites and liver transcriptomes could be used to evaluate enteric methanogenesis in Japanese Black cattle.

Proteomic analysis of adipose tissue revealing differentially abundant proteins in highly efficient mid-lactating dairy cows

The improvement of nutrient utilization efficiency in dairy cows represents an important task in view of the current rising demand for animal products and sustainable resource usage. In this perspective, the identification of appropriate markers to identify the most efficient animals for dairy production becomes a crucial factor. Residual feed intake (RFI), which represents the difference between predicted and actual intake, is used to define the efficiency of cows. In this study, subcutaneous adipose tissue (AT) was collected from five high efficient (HEF) and five low efficient (LEF) mid-lactation Holstein dairy cows, that represented subgroups of the 20% lowest RFI values (HEF) and highest 20% RFI values (LEF), out of a cohort of 155 cows that were examined for feed efficiency at the individual dairy barn at Volcani Institute, Israel. Adipose samples were examined for proteomic analysis by nano-LC/MS–MS and gene expression by RT-PCR. A total of 101 differential proteins (P ≤ 0.05 and fold change ± 1.5) and two protein networks related to feed efficiency were found between HEF and LEF cows. Among the enriched top canonical pathways, FAT10 signaling, EIF2 signaling, Sirtuin signaling, Acute phase response signaling, Protein ubiquitination and mTOR signaling pathways were related to feed efficiency in AT. Furthermore, abundance of transferrin (TF; FC = 78.35, P = 0.02) enriched pathways, including mTOR signaling, LXR/RXR and FXR/RXR activation was found in AT of HEF cows. Relative mRNA expression of RBM39, which is involved in energy metabolism, was decreased in AT of HEF versus LEF. The relationship found between the AT proteins and/or metabolic pathways and the feed efficiency demonstrates that AT may reflect metabolic adaptations to high efficiency, and suggests that these proteins together with their metabolic mechanisms are suitable candidates as biomarkers to identify efficient cows for dairy production.

Genome-wide Association Study for Carcass Primal Cut Yields Using Single-step Bayesian Approach in Hanwoo Cattle

The importance of meat and carcass quality is growing in beef cattle production to meet both producer and consumer demands. Primal cut yields, which reflect the body compositions of carcass, could determine the carcass grade and, consequently, command premium prices. Despite its importance, there have been few genome-wide association studies on these traits. This study aimed to identify genomic regions and putative candidate genes related to 10 primal cut traits, including tenderloin, sirloin, striploin, chuck, brisket, top round, bottom round, shank, flank, and rib in Hanwoo cattle using a single-step Bayesian regression (ssBR) approach. After genomic data quality control, 43,987 SNPs from 3,745 genotyped animals were available, of which 3,467 had phenotypic records for the analyzed traits. A total of 16 significant genomic regions (1-Mb window) were identified, of which five large-effect quantitative trait loci (QTLs) located on chromosomes 6 at 38–39 Mb, 11 at 21–22 Mb, 14 at 6–7 Mb and 26–27 Mb, and 19 at 26–27 Mb were associated with more than one trait, while the remaining 11 QTLs were trait-specific. These significant regions were harbored by 154 genes, among which TOX, FAM184B, SPP1, IBSP, PKD2, SDCBP, PIGY, LCORL, NCAPG, and ABCG2 were noteworthy. Enrichment analysis revealed biological processes and functional terms involved in growth and lipid metabolism, such as growth (GO:0040007), muscle structure development (GO:0061061), skeletal system development (GO:0001501), animal organ development (GO:0048513), lipid metabolic process (GO:0006629), response to lipid (GO:0033993), metabolic pathways (bta01100), focal adhesion (bta04510), ECM–receptor interaction (bta04512), fat digestion and absorption (bta04975), and Rap1 signaling pathway (bta04015) being the most significant for the carcass primal cut traits. Thus, identification of quantitative trait loci regions and plausible candidate genes will aid in a better understanding of the genetic and biological mechanisms regulating carcass primal cut yields.

Exploring Biological Impacts of Prenatal Nutrition and Selection for Residual Feed Intake on Beef Cattle Using Omics Technologies: A Review

Approximately 70% of the cost of beef production is impacted by dietary intake. Maximizing production efficiency of beef cattle requires not only genetic selection to maximize feed efficiency (i.e., residual feed intake (RFI)), but also adequate nutrition throughout all stages of growth and development to maximize efficiency of growth and reproductive capacity, even during gestation. RFI as a measure of feed efficiency in cattle has been recently accepted and used in the beef industry, but the effect of selection for RFI upon the dynamics of gestation has not been extensively studied, especially in the context of fluctuating energy supply to the dam and fetus. Nutrient restriction during gestation has been shown to negatively affect postnatal growth and development as well as fertility of beef cattle offspring. This, when combined with the genetic potential for RFI, may significantly affect energy partitioning in the offspring and subsequently important performance traits. In this review, we discuss: 1) the importance of RFI as a measure of feed efficiency and how it can affect other economic traits in beef cattle; 2) the influence of prenatal nutrition on physiological phenotypes in calves; 3) the benefits of investigating the interaction of genetic selection for RFI and prenatal nutrition; 4) how metabolomics, transcriptomics, and epigenomics have been employed to investigate the underlying biology associated with prenatal nutrition, RFI, or their interactions in beef cattle; and 5) how the integration of omics information is adding a level of deeper understanding of the genetic architecture of phenotypic traits in cattle.

Expression of candidate genes for residual feed intake in tropically adapted Bos taurus and Bos indicus bulls under thermoneutral and heat stress environmental conditions

The objectives of this study were to measure the relative expression of the ATP1A1, NR3C1, POMC, NPY, and LEP genes in Caracu (Bos taurus) and Nelore (Bos indicus) bulls submitted to feed efficiency tests at high environmental temperatures, and to evaluate differences in adaptability to tropical conditions between breeds. Thirty-five Caracu and 30 Nelore bulls were submitted to a feed efficiency test using automated feeding stations. At the end of the test, the animals were subjected to thermoneutral (TN) and heat stress (HS) conditions. Blood samples were collected after the exposure to the TN and HS conditions and the relative expression of genes was measured by qPCR. The bulls exhibited lower expression of ATP1A1 in the HS condition than in the TN condition (1.98 ± 0.27 and 2.86 ± 0.26, P = 0.02), while the relative expression of NR3C1, POMC, and LEP did not differ (P > 0.05) between climatic conditions. The breed and feed intake influenced NPY and LEP expression levels (P < 0.05). Different climate conditions associated with residual feed intake can modify the gene expression patterns of ATP1A1 and NPY. The association observed among all genes studied shows that they are involved in appetite control. Bos taurus and Bos indicus bulls exhibited similar adaptability to tropical climate conditions.

An examination of skeletal muscle and hepatic tissue transcriptomes from beef cattle divergent for residual feed intake

The selection of cattle with enhanced feed efficiency is of importance with regard to reducing feed costs in the beef industry. Global transcriptome profiling was undertaken on liver and skeletal muscle biopsies from Simmental heifers and bulls divergent for residual feed intake (RFI), a widely acknowledged feed efficiency phenotype, in order to identify genes that may be associated with this trait. We identified 5 genes (adj. p < 0.1) to be differentially expressed in skeletal muscle between high and low RFI heifers with all transcripts involved in oxidative phosphorylation and mitochondrial homeostasis. A total of 11 genes (adj. p < 0. 1) were differentially expressed in liver tissue between high and low RFI bulls with differentially expressed genes related to amino and nucleotide metabolism as well as endoplasmic reticulum protein processing. No genes were identified as differentially expressed in either heifer liver or bull muscle analyses. Results from this study show that the molecular control of RFI in young cattle is modified according to gender, which may be attributable to differences in physiological maturity between heifers and bulls of the same age. Despite this we have highlighted a number of genes that may hold potential as molecular biomarkers for RFI cattle.

Ancestral Haplotype Mapping for GWAS and Detection of Signatures of Selection in Admixed Dairy Cattle of Kenya

Understanding the genetic structure of adaptation and productivity in challenging environments is necessary for designing breeding programs that suit such conditions. Crossbred dairy cattle in East Africa resulting from over 60 years of crossing exotic dairy breeds with indigenous cattle plus inter se matings form a highly variable admixed population. This population has been subject to natural selection in response to environmental stresses, such as harsh climate, low-quality feeds, poor management, and strong disease challenge. Here, we combine two complementary sets of analyses, genome-wide association (GWA) and signatures of selection (SoS), to identify genomic regions that contribute to variation in milk yield and/or contribute to adaptation in admixed dairy cattle of Kenya. Our GWA separates SNP effects due to ancestral origin of alleles from effects due to within-population linkage disequilibrium. The results indicate that many genomic regions contributed to the high milk production potential of modern dairy breeds with no region having an exceptional effect. For SoS, we used two haplotype-based tests to compare haplotype length variation within admixed and between admixed and East African Shorthorn Zebu cattle populations. The integrated haplotype score (iHS) analysis identified 16 candidate regions for positive selection in the admixed cattle while the between population Rsb test detected 24 divergently selected regions in the admixed cattle compared to East African Shorthorn Zebu. We compare the results from GWA and SoS in an attempt to validate the most significant SoS results. Only four candidate regions for SoS intersect with GWA regions using a low stringency test. The identified SoS candidate regions harbored genes in several enriched annotation clusters and overlapped with previously found QTLs and associations for different traits in cattle. If validated, the GWA and SoS results indicate potential for SNP-based genomic selection for genetic improvement of smallholder crossbred cattle.

Testis-specific changes in gene expression of post-pubertal beef bulls divergent for residual feed intake and exposure to different pre-natal diets

Context Selection for residual feed intake (RFI) and its impact on male reproductive development has had mixed reviews in the past. Our previous studies demonstrated earlier puberty, larger testes and greater percentage of progressively motile sperm in high-RFI bulls. However, the molecular mechanisms within testes of bulls with varying RFI remain unclear. Aims To determine the effect of RFI and pre-natal diet on the expression patterns of testicular genes and use this information to explain differences observed across RFI. Methods The study included 25 purebred-Angus bulls with a genetic background of either high or low RFI and fed either normal or low pre-natal nutrition from 30 to 150 days post conception. After slaughter (17 months), testicular tissue was recovered, and RNA was extracted and sequenced. Key results Of 19218 expressed genes, 17 were differentially expressed for RFI (including PLCD1, INPP4B), with no differences being observed for pre-natal diet or diet × RFI interaction (false discovery rate) &lt; 0.1%). KEGG pathway analysis indicated that differentially expressed genes were associated with inositol phosphate metabolism, and phosphatidylinositol signalling. On the basis of a candidate gene-expression study, IGF1R was upregulated in high-RFI bulls (P &lt; 0.1). Conclusions Increased expression of IGF1R and lowered PLCD1 and INPP4B expression could activate PI3K–Akt signalling responsible for cell growth, proliferation and steroid metabolism in high-RFI bulls. Implications Selecting bulls for feed efficiency might affect molecular networks associated with reproduction and fertility.

Liver transcriptome profiling of beef steers with divergent growth rate, feed intake, or metabolic body weight phenotypes1

Average daily gain (ADG) and daily dry matter intake (DMI) are key determinants of beef industry profitability. These traits together with metabolic body weight (MWT) are combined as component traits to calculate residual feed intake (RFI), a common measure of feed efficiency in beef cattle. Recently, there have been significant efforts towards molecular genetic characterization of RFI through transcriptomic studies in different breeds and tissues. However, molecular mechanisms of RFI component traits still remain predominately unexplored. Therefore, in the current study, we investigated the hepatic transcriptomic profiles and their associations with ADG, DMI, and MWT in Angus, Charolais, and Kinsella Composite (KC) populations through global RNAseq analyses. In each population and for each trait, 12 steers with extreme phenotypes (n = 6 low and n = 6 high) were analyzed for differential gene expression. These animals were from 20 beef steers of each Angus, Charolais, and KC breed population that were initially selected for a transcriptome study of RFI. At a false discovery rate <0.05 and fold change >1.5, we identified 123, 102, and 78 differentially expressed (DE) genes between high- and low-ADG animals of Angus, Charolais, and KC populations, respectively. For DMI, 108, 180, and 156 DE genes were identified between high- and low-DMI from Angus, Charolais, and KC populations, respectively, while for MWT, 80, 82, and 84 genes were differentially expressed between high- and low-MWT animals in Angus, Charolais, and KC populations, respectively. The identified DE genes were largely breed specific (81.7% for ADG, 82.7% for DMI, and 83% for MWT), but were largely involved in the same biological functions across the breeds. Among the most enriched biological functions included metabolism of major nutrients (lipids, carbohydrates, amino acids, vitamins, and minerals), small molecule biochemistry, cellular movement, cell morphology, and cell-to-cell signaling and interaction. Notably, we identified multiple DE genes that are involved in cholesterol biosynthesis, and immune response pathways for the 3 studied traits. Thus, our findings present potential molecular genetic mechanisms and candidate genes that influence feed intake, growth, and MWT of beef cattle.

Population structure and genomic inbreeding in nine Swiss dairy cattle populations

BACKGROUND

Domestication, breed formation and intensive selection have resulted in divergent cattle breeds that likely exhibit their own genomic signatures. In this study, we used genotypes from 27,612 autosomal single nucleotide polymorphisms to characterize population structure based on 9214 sires representing nine Swiss dairy cattle populations: Brown Swiss (BS), Braunvieh (BV), Original Braunvieh (OB), Holstein (HO), Red Holstein (RH), Swiss Fleckvieh (SF), Simmental (SI), Eringer (ER) and Evolèner (EV). Genomic inbreeding (F ROH) and signatures of selection were determined by calculating runs of homozygosity (ROH). The results build the basis for a better understanding of the genetic development of Swiss dairy cattle populations and highlight differences between the original populations (i.e. OB, SI, ER and EV) and those that have become more popular in Switzerland as currently reflected by their larger populations (i.e. BS, BV, HO, RH and SF).

RESULTS

The levels of genetic diversity were highest and lowest in the SF and BS breeds, respectively. Based on F ST values, we conclude that, among all pairwise comparisons, BS and HO (0.156) differ more than the other pairs of populations. The original Swiss cattle populations OB, SI, ER, and EV are clearly genetically separated from the Swiss cattle populations that are now more common and represented by larger numbers of cows. Mean levels of F ROH ranged from 0.027 (ER) to 0.091 (BS). Three of the original Swiss cattle populations, ER (F ROH: 0.027), OB (F ROH: 0.029), and SI (F ROH: 0.039), showed low levels of genomic inbreeding, whereas it was much higher in EV (F ROH: 0.074). Private signatures of selection for the original Swiss cattle populations are reported for BTA4, 5, 11 and 26.

CONCLUSIONS

The low levels of genomic inbreeding observed in the original Swiss cattle populations ER, OB and SI compared to the other breeds are explained by a lesser use of artificial insemination and greater use of natural service. Natural service results in more sires having progeny at each generation and thus this breeding practice is likely the major reason for the remarkable levels of genetic diversity retained within these populations. The fact that the EV population is regionally restricted and its small census size of herd-book cows explain its high level of genomic inbreeding.

Comparative transcriptome analysis of the Pacific White Shrimp (Litopenaeus vannamei) muscle reveals the molecular basis of residual feed intake

Feed efficiency is an economically important trait in genetic improvement programs of L. vannamei. Residual feed intake (RFI), an ideal measure of feed efficiency, is the difference between observed feed intake and expected feed requirement predicted from maintenance and production. Exploring the molecular basis of RFI is essential to facilitate the genetic breeding of feed efficiency in L. vannamei. However, few studies have been reported in this aspect. In this study, we sequenced muscle transcriptomes of a high-efficiency group, a low-efficiency group and a control group originating from two families, and compared the gene expression patterns between each extreme group and the control group. A total of 383 differentially expressed genes were identified, most of which were involved in cell proliferation, growth and signaling, glucose homeostasis, energy and nutrients metabolism. Functional enrichment analysis of these genes revealed 13 significantly enriched biological pathways, including signaling pathways such as PI3K-Akt signaling pathway, AMPK signaling pathway and mTOR signaling pathway, as well as some important pathways such as ubiquitin mediated proteolysis, cell cycle, pentose phosphate pathway and glycolysis/gluconeogenesis. These genes and pathways provide initial insight into the molecular mechanisms driving the feed efficiency in L. vannamei.

Searching for new loci and candidate genes for economically important traits through gene-based association analysis of Simmental cattle

Single-marker genome-wide association study (GWAS) is a convenient strategy of genetic analysis that has been successful in detecting the association of a number of single-nucleotide polymorphisms (SNPs) with quantitative traits. However, analysis of individual SNPs can only account for a small proportion of genetic variation and offers only limited knowledge of complex traits. This inadequacy may be overcome by employing a gene-based GWAS analytic approach, which can be considered complementary to the single-SNP association analysis. Here we performed an initial single-SNP GWAS for bone weight (BW) and meat pH value with a total of 770,000 SNPs in 1141 Simmental cattle. Additionally, 21836 cattle genes collected from the Ensembl Genes 83 database were analyzed to find supplementary evidence to support the importance of gene-based association study. Results of the single SNP-based association study showed that there were 11 SNPs significantly associated with bone weight (BW) and two SNPs associated with meat pH value. Interestingly, all of these SNPs were located in genes detected by the gene-based association study.

Gene expression differences in Longissimus muscle of Nelore steers genetically divergent for residual feed intake

Residual feed intake (RFI), a measure of feed efficiency (FE), is defined as the difference between the observed and the predictable feed intake considering size and growth of the animal. It is extremely important to beef production systems due to its impact on the allocation of land areas to alternative agricultural production, animal methane emissions, food demand and cost of production. Global differential gene expression analysis between high and low RFI groups (HRFI and LRFI: less and more efficient, respectively) revealed 73 differentially expressed (DE) annotated genes in Longissimus thoracis (LT) muscle of Nelore steers. These genes are involved in the overrepresented pathways Metabolism of Xenobiotics by Cytochrome P450 and Butanoate and Tryptophan Metabolism. Among the DE transcripts were several proteins related to mitochondrial function and the metabolism of lipids. Our findings indicate that observed gene expression differences are primarily related to metabolic processes underlying oxidative stress. Genes involved in the metabolism of xenobiotics and antioxidant mechanisms were primarily down-regulated, while genes responsible for lipid oxidation and ketogenesis were up-regulated in HRFI group. By using LT muscle, this study reinforces our previous findings using liver tissue and reveals new genes and likely tissue-specific regulators playing key-roles in these processes.

Post-weaning blood transcriptomic differences between Yorkshire pigs divergently selected for residual feed intake

Background

Improving feed efficiency (FE) of pigs by genetic selection is of economic and environmental significance. An increasingly accepted measure of feed efficiency is residual feed intake (RFI). Currently, the molecular mechanisms underlying RFI are largely unknown. Additionally, to incorporate RFI into animal breeding programs, feed intake must be recorded on individual pigs, which is costly and time-consuming. Thus, convenient and predictive biomarkers for RFI that can be measured at an early age are greatly desired. In this study, we aimed to explore whether differences exist in the global gene expression profiles of peripheral blood of 35 to 42 day-old pigs with extremely low (more efficient) and high RFI (less efficient) values from two lines that were divergently selected for RFI during the grow-finish phase, to use such information to explore the potential molecular basis of RFI differences, and to initiate development of predictive biomarkers for RFI.

Results

We identified 1972 differentially expressed genes (DEGs) (q ≤ 0.15) between the low (n = 15) and high (n = 16) RFI groups of animals by using RNA sequencing technology. We validated 24 of 37 selected DEGs by reverse transcription-quantitative PCR (RT-qPCR) in a joint analysis of 24 (12 per line) of the 31 samples already used for RNA-seq plus 24 (12 per line) novel samples from the same contemporary group of pigs. Using an analysis of the 24 novel samples alone, only nine of the 37 selected DEGs were validated. Genes involved in small molecule biosynthetic process, antigen processing and presentation of peptide antigen via major histocompatibility complex (MHC) class I, and steroid biosynthetic process were overrepresented among DEGs that had higher expression in the low versus high RFI animals. Genes known to function in the proteasome complex or mitochondrion were also significantly enriched among genes with higher expression in the low versus high RFI animals. Alternatively, genes involved in signal transduction, bone mineralization and regulation of phosphorylation were overrepresented among DEGs with lower expression in the low versus high RFI animals. The DEGs significantly overlapped with genes associated with disease, including hyperphagia, eating disorders and mitochondrial diseases (q < 1E-05). A weighted gene co-expression network analysis (WGCNA) identified four co-expression modules that were differentially expressed between the low and high RFI groups. Genes involved in lipid metabolism, regulation of bone mineralization, cellular immunity and response to stimulus were overrepresented within the two modules that were most significantly differentially expressed between the low and high RFI groups. We also found five of the DEGs and one of the co-expression modules were significantly associated with the RFI phenotype of individual animals (q < 0.05).

Conclusions

The post-weaning blood transcriptome was clearly different between the low and high RFI groups. The identified DEGs suggested potential differences in mitochondrial and proteasomal activities, small molecule biosynthetic process, and signal transduction between the two RFI groups and provided potential new insights into the molecular basis of RFI in pigs, although the observed relationship between the post-weaning blood gene expression and RFI phenotype measured during the grow-finish phase was not strong. DEGs and representative genes in co-expression modules that were associated with RFI phenotype provide a preliminary list for developing predictive biomarkers for RFI in pigs.

Electronic supplementary material

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

Liver transcriptomic networks reveal main biological processes associated with feed efficiency in beef cattle

BACKGROUND

The selection of beef cattle for feed efficiency (FE) traits is very important not only for productive and economic efficiency but also for reduced environmental impact of livestock. Considering that FE is multifactorial and expensive to measure, the aim of this study was to identify biological functions and regulatory genes associated with this phenotype.

RESULTS

Eight genes were differentially expressed between high and low feed efficient animals (HFE and LFE, respectively). Co-expression analyses identified 34 gene modules of which 4 were strongly associated with FE traits. They were mainly enriched for inflammatory response or inflammation-related terms. We also identified 463 differentially co-expressed genes which were functionally enriched for immune response and lipid metabolism. A total of 8 key regulators of gene expression profiles affecting FE were found. The LFE animals had higher feed intake and increased subcutaneous and visceral fat deposition. In addition, LFE animals showed higher levels of serum cholesterol and liver injury biomarker GGT. Histopathology of the liver showed higher percentage of periportal inflammation with mononuclear infiltrate.

CONCLUSION

Liver transcriptomic network analysis coupled with other results demonstrated that LFE animals present altered lipid metabolism and increased hepatic periportal lesions associated with an inflammatory response composed mainly by mononuclear cells. We are now focusing to identify the causes of increased liver lesions in LFE animals.