RNA-seq analysis of bovine adipose tissue in heifers fed diets differing in energy and protein content

A review of the role of transcription factors in regulating adipogenesis and lipogenesis in beef cattle

In the past few decades, genomic selection and other refined strategies have been used to increase the growth rate and lean meat production of beef cattle. Nevertheless, the fast growth rates of cattle breeds are often accompanied by a reduction in intramuscular fat (IMF) deposition, impairing meat quality. Transcription factors play vital roles in regulating adipogenesis and lipogenesis in beef cattle. Meanwhile, understanding the role of transcription factors in regulating adipogenesis and lipogenesis in beef cattle has gained significant attention to increase IMF deposition and meat quality. Therefore, the aim of this paper was to provide a comprehensive summary and valuable insight into the complex role of transcription factors in adipogenesis and lipogenesis in beef cattle. This review summarizes the contemporary studies in transcription factors in adipogenesis and lipogenesis, genome-wide analysis of transcription factors, epigenetic regulation of transcription factors, nutritional regulation of transcription factors, metabolic signalling pathways, functional genomics methods, transcriptomic profiling of adipose tissues, transcription factors and meat quality and comparative genomics with other livestock species. In conclusion, transcription factors play a crucial role in promoting adipocyte development and fatty acid biosynthesis in beef cattle. They control adipose tissue formation and metabolism, thereby improving meat quality and maintaining metabolic balance. Understanding the processes by which these transcription factors regulate adipose tissue deposition and lipid metabolism will simplify the development of marbling or IMF composition in beef cattle.

Gene expression of free fatty acids-sensing G protein-coupled receptors in beef cattle

Many physiological functions are regulated by free fatty acids (FFA). Recently, the discovery of FFA-specific G protein-coupled receptors (FFARs) has added to the complexity of their actions at the cellular level. The study of FFAR in cattle is still in its earliest stages focusing mainly on dairy cows. In this study, we set out to map the expression of genes encoding FFARs in 6 tissues of beef cattle. We also investigated the potential effect of dietary forage nature on FFAR gene expression. To this end, 16 purebred Charolais bulls were fed a grass silage ration or a maize silage ration (n = 8/group) with a forage/concentrate ratio close to 60:40 for 196 d. The animals were then slaughtered at 485 ± 42 d and liver, spleen, ileum, rectum, perirenal adipose tissue (PRAT), and Longissimus Thoracis muscle were collected. FFAR gene expression was determined by real-time quantitative PCR. Our results showed that of the five FFARs investigated, FFAR1, FFAR2, FFAR3, and GPR84 are expressed (Ct < 30) in all six tissues, whereas FFAR4 was only expressed (Ct < 30) in PRAT, ileum, and rectum. In addition, our results showed that the nature of the forage, i.e., grass silage or maize silage, had no effect on the relative abundance of FFAR in any of the tissues studied (P value > 0.05). Taken together, these results open new perspectives for studying the physiological role of these receptors in beef cattle, particularly in nutrient partitioning during growth.

Fat deposition, fatty acid profiles, antioxidant capacity and differentially expressed genes in subcutaneous fat of Tibetan sheep fed wheat-based diets with and without xylanase supplementation

Xylanase, an exogenous enzyme that plays an essential role in energy metabolism by hydrolysing xylan into xylose, has been shown to positively influence nutrient digestion and utilisation in ruminants. The objective of this study was to evaluate the effects of xylanase supplementation on the back-fat thickness, fatty acid profiles, antioxidant capacity, and differentially expressed genes (DEGs) in the subcutaneous fat of Tibetan sheep. Sixty three-month-old rams with an average weight of 19.35 ± 2.18 kg were randomly assigned to control (no enzyme added, WH group) and xylanase (0.2% of diet on a dry matter basis, WE group) treatments. The experiment was conducted over 97 d, including 7 d of adaption to the diets. The results showed that xylanase supplementation in the diet increased adipocyte volume of subcutaneous fat (p < 0.05), shown by hematoxylin and eosin (H&E) staining. Gas chromatography showed greater concentrations of C14:0 and C16:0 in the subcutaneous fat of controls compared with the enzyme-treated group (p < 0.05), while opposite trend was seen for the absolute contents of C18:1n9t, C20:1, C18:2n6c, C18:3, and C18:3n3 (p < 0.05). Compared with controls, supplementation with xylanase increased the activity of T-AOC significantly (p < 0.05). Transcriptomic analysis showed the presence of 1630 DEGs between the two groups, of which 1023 were up-regulated and 607 were down-regulated, with enrichment in 4833 Gene Ontology terms, and significant enrichment in 31 terms (p < 0.05). The common DEGs were enriched in 295 pathways and significantly enriched in 26 pathways. Additionally, the expression of lipid-related genes, including fatty acid synthase, superoxide dismutase, fatty acid binding protein 5, carnitine palmytoyltransferase 1 A, and peroxisome proliferator-activated receptor A were verified via quantitative reverse-transcription polymerase chain reaction. In conclusion, dietary xylanase supplementation was found to reduce subcutaneous fat deposition in Tibetan sheep, likely through modulating the expression of lipid-related genes.

Identification of the hub genes related to adipose tissue metabolism of bovine

Due to the demand for high-quality animal protein, there has been consistent interest in how to obtain more high-quality beef. As well-known, the adipose content of beef has a close connection with the taste and quality of beef, and cattle with different energy or protein diet have corresponding effects on the lipid metabolism of beef. Thus, we performed weighted gene co-expression network analysis (WGCNA) with subcutaneous adipose genes from Norwegian red heifers fed different diets to identify hub genes regulating bovine lipid metabolism. For this purpose, the RNA sequencing data of subcutaneous adipose tissue of 12-month-old Norwegian red heifers (n = 48) with different energy or protein levels were selected from the GEO database, and 7,630 genes with the largest variation were selected for WGCNA analysis. Then, three modules were selected as hub genes candidate modules according to the correlation between modules and phenotypes, including pink, magenta and grey60 modules. GO and KEGG enrichment analysis showed that genes were related to metabolism, and participated in Rap, MAPK, AMPK, VEGF signaling pathways, and so forth. Combined gene interaction network analysis using Cytoscape software, eight hub genes of lipid metabolism were identified, including TIA1, LOC516108, SNAPC4, CPSF2, ZNF574, CLASRP, MED15 and U2AF2. Further, the expression levels of hub genes in the cattle tissue were also measured to verify the results, and we found hub genes in higher expression in muscle and adipose tissue in adult cattle. In summary, we predicted the key genes of lipid metabolism in the subcutaneous adipose tissue that were affected by the intake of various energy diets to find the hub genes that coordinate lipid metabolism, which provide a theoretical basis for regulating beef quality.

Analyzing the interactions of mRNAs, miRNAs and lncRNAs to predict ceRNA networks in bovine cystic follicular granulosa cells

Cross-talk between competitive endogenous RNAs (ceRNAs) may play a critical role in revealing potential mechanism of bovine follicular cysts. Ovarian cyst has always been an intractable scientific problem and has led to considerable economic losses to bovine breeding industry. However, its pathogenesis and molecular mechanisms are still not well understood. Here, this study aimed to investigate the role of non-coding RNAs (ncRNAs) and the ceRNA networks in bovine follicular cyst. Whole transcriptome sequencing of bovine follicular granulosa cells (GCs) was conducted to obtain the expression profiles of mRNAs, lncRNAs and miRNAs. The results for the identified expressions of 8,003 mRNAs, 579 lncRNAs and 205 miRNAs were often altered between cystic and normal follicular GCs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed on these differentially expressed mRNAs. Furthermore, the ceRNA network combining mRNAs, miRNAs, and lncRNAs using several bioinformatics methods based on co-expression analysis between the differentially expressed RNAs was conducted. Finally, the lncRNA NONBTAT027373.1-miR-664b-HSD17B7 pathway was verified by dual-luciferase reporting assay and RNA binding protein immunoprecipitation (RIP) assay. LncRNA NONBTAT027373.1 sponged miR-664b in GCs and prevented miR-664b from binding to the HSD17B7 3′-UTR. These results indicated that genes and lncRNAs related to steroid hormone synthesis and energy metabolism could play important roles in the formation of bovine cystic follicles through the ceRNA mechanism and represent candidate targets for further research. This can be used as a practical guideline for promoting healthy and highly efficient development in the bovine industry.

Comparative transcriptome analysis of goat (Capra hircus) adipose tissue reveals physiological regulation of body reserve recovery after the peak of lactation

Adipose tissue is the energy storage organ providing energy to other tissues, including mammary gland, that supports the achievement of successive lactation cycles. Our objective was to investigate the ability of goats to restore body fat reserves by comparing lipogenic enzyme activities and by transcriptomic RNA-Seq data at two different physiological stages, mid- and post-lactation. Key lipogenic enzyme activities were higher in goat omental adipose tissue during mid-lactation (74 days in milk) than during the post-lactation period (300 days postpartum). RNA-Sequencing analysis revealed 19,271 expressed genes in the omental adipose tissue. The comparison between adipose transcriptome analysis from mid- and post-lactation goats highlighted 252 differentially expressed genes (p_adj < 0.05) between these two physiological stages. The differential expression of 11 genes was confirmed by RT-qPCR. Functional genomic analysis revealed that 31% were involved in metabolic processes among which 38% in lipid metabolism. Most of the genes involved in lipid synthesis and those in lipid transport and storage were upregulated in adipose tissue of mid- compared to post-lactation goats. In addition, adipose tissue plasticity was emphasized by genes involved in cellular signaling and tissue integrity. Network analyses also highlighted three key regulators of lipid metabolism (LEP, APOE and HNF4A) and a key target gene (VCAM1). The greatest lipogenic enzyme activities with the upregulation of genes involved in lipid metabolism highlighted a higher recovery of lipid reserves after the lactation peak than 4 months post-lactation. This study contributes to a better understanding of the molecular mechanisms controlling the body lipid reserves management during the successive lactations.

Gene co-expression and alternative splicing analysis of key metabolic tissues to unravel the regulatory signatures of fatty acid composition in cattle

Increasing the healthy/unhealthy fatty acid (FA) ratio in meat is one of the urgent tasks required to address consumer concerns. However, the regulatory mechanisms ultimately resulting in FA profiles vary among animals and remain largely unknown. In this study, using ~1.2 Tb high-quality RNA-Seq-based transcriptomic data of 188 samples from four key metabolic tissues (rumen, liver, muscle, and backfat) together with the contents of 49 FAs in backfat, the molecular regulatory mechanisms of these tissues contributing to FA formation in cattle were explored. Using this large dataset, the alternative splicing (AS) events, one of the transcriptional regulatory mechanisms in four tissues were identified. The highly conserved and absent AS events were detected in rumen tissue, which may contribute to its functional differences compared with the other three tissues. In addition, the healthy/unhealthy FA ratio related AS events, differential expressed (DE) genes, co-expressed genes, and their functions in four tissues were analysed. Eight key genes were identified from the integrated analysis of DE, co-expressed, and AS genes between animals with high and low healthy/unhealthy FA ratios. This study provides an applicable pipeline for AS events based on comprehensive RNA-Seq analysis and improves our understanding of the regulatory mechanism of FAs in beef cattle.

Diet induced the change of mtDNA copy number and metabolism in Angus cattle

Background

Grass-fed and grain-fed Angus cattle differ in the diet regimes. However, the intricate mechanisms of different beef quality and other phenotypes induced by diet differences are still unclear. Diet affects mitochondrial function and dynamic behavior in response to changes in energy demand and supply. In this study, we examined the mtDNA copy number, mitochondria-related genes expression, and metabolic biomarkers in grass-fed and grain-fed Angus cattle.

Results

We found that the grass-fed group had a higher mtDNA copy number than the grain-fed group. Among different tissues, the mtDNA copy number was the highest in the liver than muscle, rumen, and spleen. Based on the transcriptome of the four tissues, a lower expression of mtDNA-encoded genes in the grass-fed group compared to the grain-fed group was discovered. For the mitochondria-related nuclear genes, however, most of them were significantly down-regulated in the muscle of the grass-fed group and up-regulated in the other three tissues. In which, COX6A2, POLG2, PPIF, DCN, and NDUFA12, involving in ATP synthesis, mitochondrial replication, transcription, and maintenance, might contribute to the alterations of mtDNA copy number and gene expression. Meanwhile, 40 and 23 metabolic biomarkers were identified in the blood and muscle of the grain-fed group compared to a grass-fed group, respectively. Integrated analysis of the altered metabolites and gene expression revealed the high expression level of MDH1 in the grain-fed group might contribute to the mitochondrial NADH oxidation and spermidine metabolism for adapting the deletion mtDNA copy number.

Conclusions

Overall, the study may provide further deep insight into the adaptive and regulatory modulations of the mitochondrial function in response to different feeding systems in Angus cattle.

Transcriptomic changes in bovine skeletal muscle cells after resveratrol treatment

Skeletal muscles constitute a high proportion of the cellular mass that is essential for the growth traits in cattle. Resveratrol (RSV) is a natural polyphenol compound involved in pleiotropic biological activities of muscle. Therefore, the aim of our study was to investigate the transcriptome-level effects of RSV on bovine primary myoblast to reveal differentially expressed genes (DEGs). We treated three replicates of primary myoblasts with 20 μM mother solution containing RSV, whereas three other replicates without RSV were used as control group. Then, we conducted genome-wide transcriptome analysis for the two groups. The results of expression analysis identified 3856 DEGs of which 1805 genes were up-regulated and 2051 genes were down-regulated (adjusted P < 0.05). In addition, qRT-PCR analysis of 19 selected DEGs were consistent with the expression levels observed in the transcriptome data. Gene Ontology (GO) and pathway enrichment analysis showed 72 and 66 significant GO terms and KEGG pathways, respectively (adjusted P < 0.05). The most significant GO term was actin cytoskeleton organization (GO:0030036). The top significant KEGG pathway was focal adhesion (bta04510). Predicted protein-protein interactions (PPIs) showed that CDKN1A encoding cyclindependent kinase inhibitor 1A connects several larger protein complexes. In conclusion, our results found a list of DEGs, significant GO terms and pathways, and provided an improved and expanded understanding of the impact of RSV on cattle muscle cells at the transcriptomic level. The study elucidates the potential of using the genes enriched in pathways mediating resveratrol effects as targets in genomic selection for muscle development and growth in beef cattle.

From Beef to Bees: High-Throughput Kinome Analysis to Understand Host Responses of Livestock Species to Infectious Diseases and Industry-Associated Stress

Within human health research, the remarkable utility of kinase inhibitors as therapeutics has motivated efforts to understand biology at the level of global cellular kinase activity (the kinome). In contrast, the diminished potential for using kinase inhibitors in food animals has dampened efforts to translate this research approach to livestock species. This, in our opinion, was a lost opportunity for livestock researchers given the unique potential of kinome analysis to offer insight into complex biology. To remedy this situation, our lab developed user-friendly, cost-effective approaches for kinome analysis that can be readily incorporated into most research programs but with a specific priority to enable the technology to livestock researchers. These contributions include the development of custom software programs for the creation of species-specific kinome arrays as well as comprehensive deconvolution and analysis of kinome array data. Presented in this review are examples of the application of kinome analysis to highlight the utility of the technology to further our understanding of two key complex biological events of priority to the livestock industry: host immune responses to infectious diseases and animal stress responses. These advances and examples of application aim to provide both mechanisms and motivation for researchers, particularly livestock researchers, to incorporate kinome analysis into their research programs.

Impact of the severity of negative energy balance on gene expression in the subcutaneous adipose tissue of periparturient primiparous Holstein dairy cows: Identification of potential novel metabolic signals for the reproductive system

The severity of negative energy balance (NEB) in high-producing dairy cows has a high incidence among health diseases. The cow’s energy status during early lactation critically affects metabolic and reproductive parameters. The first objective of this study was to investigate by RNA-seq analysis and RT-qPCR the gene expression profile in white adipose tissue and by gene ontology and upstream regulation tools the relationships with energy metabolism and reproduction in two groups of primiparous dairy cows with extreme NEB statuses (NEB < -9 Mcal/day vs. NEB > -9 Mcal/day) around parturition. The second objective was to determine the potential involvement of a new adipokine identified as a candidate for the regulation of ovarian function in our RNA-seq analysis by using bovine primary granulosa culture, thymidine incorporation to determine cell proliferation and ELISA assays to measure progesterone secretion. The RNA-seq analysis revealed that 514 genes were over-expressed and 695 were under-expressed in the adipose tissue of cows with severe NEB (SNEB) and cows with moderate NEB (MNEB) during the -4 and 16 wkpp period. In addition, 491 genes were over-expressed and 705 genes were under-expressed in the adipose tissue of SNEB cows compared to MNEB cows. Among these differently expressed genes (DEGs), 298 were related to metabolic functions and 264 to reproductive traits. A set of 19 DEGs were validated by RT-qPCR, including CCL21 (C-C motif chemokine ligand 21). Moreover, CCL21, a gene known to be secreted by adipose tissue, was chosen for further analysis in plasma and ovaries. The use of next-generation sequencing technologies allowed us to characterise the transcriptome of white adipose tissue from primiparous cows with different levels of NEB during lactation. This study highlighted the alteration of the expression of genes related to lipid metabolism, including CCL21, which is released in the bloodstream and associated with the in vitro regulation of ovarian functions.

Long-term effects of prior diets, dietary transition and pregnancy on adipose gene expression in dairy heifers

Adipose tissue is highly involved in whole-body metabolism and is the main site for lipid synthesis, storage and mobilization in ruminants. Therefore, knowledge about adipose tissue responses to different diets is important, especially in growing heifers as the feeding regimes of replacement heifers affect their future success as dairy cows. However, at gene expression level such knowledge is limited. As part of a larger feed trial, adipose tissue biopsies from 24 Norwegian Red heifers were collected at 12 months of age (12MO) and at month seven of gestation (PREG) and analyzed by next-generation mRNA sequencing. Between these two sampling points, all heifers had gone through a successful conception and a feed change from four dietary treatments of high or low energy (HE/LE) and protein (HP/LP) content (treatments LPHE, HPHE, LPLE and HPLE) to a low-energy, low-protein pregnancy feed given to all animals. Gene expression differences between different feed treatments at 12MO are described in an earlier publication from our group. The main objectives of this study were to investigate the long-term effects of diets differing in protein and energy density level on gene expression in adipose tissue of growing replacement dairy heifers. To achieve this, we examined the post-treatment effects between the treatment groups at month seven of gestation; 6 months after the termination of experimental feeding, and the long-term gene expression changes occurring in the adipose tissue between 12MO and PREG. Post-treatment group comparisons showed evidence of long-term effects of dietary treatment on adipose gene expression. Differences between protein treatments were smaller than between energy treatments. Adipose gene expression changes from 12MO to PREG were much larger for the HE than the LE treatments and seemed to mostly be explained by the characteristics of the diet change. 97 genes displayed a unidirectional expression change for all groups from 12MO to PREG, and are considered to be treatment-independent, possibly caused by pregnancy or increased age. This study provides candidate genes and key regulators for further studies on pregnancy preservation (TGFB1, CFD) and metabolic regulation and efficiency (PI3K, RICTOR, MAP4K4,) in dairy cattle.