Gene Co-expression Analysis Indicates Potential Pathways and Regulators of Beef Tenderness in Nellore Cattle

Expression network analysis of bovine skin infested with Rhipicephalus australis identifies pro-inflammatory genes contributing to tick susceptibility

The skin is the primary feeding site of ticks that infest livestock animals such as cattle. The highly specialised functions of skin at the molecular level may be a factor contributing to variation in susceptibility to tick infestation; but these remain to be well defined. The aim of this study was to investigate the bovine skin transcriptomic profiles of tick-naïve and tick-infested cattle and to uncover the gene expression networks that influence contrasting phenotypes of host resistance to ticks. RNA-Seq data was obtained from skin of Brangus cattle with high (n = 5) and low (n = 6) host resistance at 0 and 12 weeks following artificial tick challenge with Rhipicephalus australis larvae. No differentially expressed genes were detected pre-infestation between high and low resistance groups, but at 12-weeks there were 229 differentially expressed genes (DEGs; FDR < 0.05), of which 212 were the target of at least 1866 transcription factors (TFs) expressed in skin. Regulatory impact factor (RIF) analysis identified 158 significant TFs (P < 0.05) of which GRHL3, and DTX1 were also DEGs in the experiment. Gene term enrichment showed the significant TFs and DEGs were enriched in processes related to immune response and biological pathways related to host response to infectious diseases. Interferon Type 1-stimulated genes, including MX2, ISG15, MX1, OAS2 were upregulated in low host resistance steers after repeated tick challenge, suggesting dysregulated wound healing and chronic inflammatory skin processes contributing to host susceptibility to ticks. The present study provides an assessment of the bovine skin transcriptome before and after repeated tick challenge and shows that the up-regulation of pro-inflammatory genes is a prominent feature in the skin of tick-susceptible animals. In addition, the identification of transcription factors with high regulatory impact provides insights into the potentially meaningful gene-gene interactions involved in the variation of phenotypes of bovine host resistance to ticks.

A network-based approach to understanding gene-biological processes affecting economically important traits of Nelore cattle

This study aimed to build gene-biological process networks with differentially expressed genes associated with economically important traits of Nelore cattle from 17 previous studies. The genes were clustered into three groups by evaluated traits: group 1, production traits; group 2, carcass traits; and group 3, meat quality traits. For each group, a gene-biological process network analysis was performed with the differentially expressed genes in common. For production traits, 37 genes were found in common, of which 13 genes were enriched for six Gene Ontology (GO) terms; these terms were not functionally grouped. However, the enriched GO terms were related to homeostasis, the development of muscles and the immune system. For carcass traits, four genes were found in common. Thus, it was not possible to functionally group these genes into a network. For meat quality traits, the analysis revealed 222 genes in common. CSRP3 was the only gene differentially expressed in all three groups. Non-redundant biological terms for clusters of genes were functionally grouped networks, reflecting the cross-talk between all biological processes and genes involved. Many biological processes and pathways related to muscles, the immune system and lipid metabolism were enriched, such as striated muscle cell development and triglyceride metabolic processes. This study provides insights into the genetic mechanisms of production, carcass and meat quality traits of Nelore cattle. This information is fundamental for a better understanding of the complex traits and could help in planning strategies for the production and selection systems of Nelore cattle.

Molecular mechanisms contributing to the development of beef sensory texture and flavour traits and related biomarkers: Insights from early post-mortem muscle using label-free proteomics

Beef sensory quality comprises a suite of traits, each of which manifests its ultimate phenotype through interaction of muscle physiology with environment, both in vivo and post-mortem. Understanding variability in meat quality remains a persistent challenge, but omics studies to uncover biological connections between natural variability in proteome and phenotype could provide validation for exploratory studies and offer new insights. Multivariate analysis of proteome and meat quality data from Longissimus thoracis et lumborum muscle samples taken early post-mortem from 34 Limousin-sired bulls was conducted. Using for the first-time label-free shotgun proteomics combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS), 85 proteins were found to be related with tenderness, chewiness, stringiness and flavour sensory traits. The putative biomarkers were classified in five interconnected biological pathways; i) muscle contraction, ii) energy metabolism, iii) heat shock proteins, iv) oxidative stress, v) regulation of cellular processes and binding. Among the proteins, PHKA1 and STBD1 correlated with all four traits, as did the GO biological process 'generation of precursor metabolites and energy'. Optimal regression models explained a high level (58-71%) of phenotypic variability with proteomic data for each quality trait. The results of this study propose several regression equations and biomarkers to explain the variability of multiple beef eating quality traits. Thanks to annotation and network analyses, they further suggest protein interactions and mechanisms underpinning the physiological processes regulating these key quality traits. SIGNIFICANCE: The proteomic profiles of animals with divergent quality profiles have been compared in numerous studies; however, a wide range of phenotypic variation is required to better understand the mechanisms underpinning the complex biological pathways correlated with beef quality and protein interactions. We used multivariate regression analyses and bioinformatics to analyse shotgun proteomics data to decipher the molecular signatures involved in beef texture and flavour variations with a focus on multiple quality traits. We developed multiple regression equations to explain beef texture and flavour. Additionally, potential candidate biomarkers correlated with multiple beef quality traits are suggested, which could have utility as indicators of beef overall sensory quality. This study explained the biological process responsible for determining key quality traits such as tenderness, chewiness, stringiness, and flavour in beef, which will provide support for future beef proteomics studies.

First insights into the dynamic protein changes in goat Semitendinosus muscle during the post-mortem period using high-throughput proteomics

Proteomics plays a key and insightful role in meat research in the post-genomic era. This study aimed to unveil using a shotgun proteomics approach the temporal dynamic changes in early post-mortem proteome of goat Semitendinosus muscle. Therefore, the evolution and comparison of the muscle proteome over three post-mortem times (1, 8, and 24 h) was assessed. The temporal proteomics profiling quantified 748 proteins, from which 174 were differentially abundant (DAPs): n = 55 between 1 h versus 8 h, n = 52 between 8 h versus 24 h, and n = 154 between 1 h versus 24 h. The DAPs belong to myriad interconnected pathways. Binding, transport and calcium homeostasis, as well as muscle contraction and structure, exhibited an equivalent contribution during post-mortem, demonstrating their central role. Catalytic, metabolism and ATP metabolic process, and proteolysis were active pathways from the first hours of animal bleeding. Conversely, oxidative stress, response to hypoxia and cell redox homeostasis along chaperones and heat shock proteins accounted for the large proportion of the biochemical processes, more importantly after 8 h post-mortem. Overall, the conversion of muscle into meat is largely orchestrated by energy production as well as mitochondrial metabolism and homeostasis through calcium and permeability transition regulation. The study further evidenced the role of ribosomal proteins in goat post-mortem muscle, signifying that several proteins experiencing changes during storage, also undergo splicing modifications, which is for instance a mechanism known for mitochondrial proteins. Overall, temporal proteomics profiling of early post-mortem muscle proteome offers an unparalleled view of the sophisticated post-mortem biochemical and proteolytic events associated with goat meat quality determination.

Differential Allele-Specific Expression Revealed Functional Variants and Candidate Genes Related to Meat Quality Traits in B. indicus Muscle

Traditional transcriptomics approaches have been used to identify candidate genes affecting economically important livestock traits. Regulatory variants affecting these traits, however, remain under covered. Genomic regions showing allele-specific expression (ASE) are under the effect of cis-regulatory variants, being useful for improving the accuracy of genomic selection models. Taking advantage of the better of these two methods, we investigated single nucleotide polymorphisms (SNPs) in regions showing differential ASE (DASE SNPs) between contrasting groups for beef quality traits. For these analyses, we used RNA sequencing data, imputed genotypes and genomic estimated breeding values of muscle-related traits from 190 Nelore (Bos indicus) steers. We selected 40 contrasting unrelated samples for the analysis (N = 20 animals per contrasting group) and used a beta-binomial model to identify ASE SNPs in only one group (i.e., DASE SNPs). We found 1479 DASE SNPs (FDR ≤ 0.05) associated with 55 beef-quality traits. Most DASE genes were involved with tenderness and muscle homeostasis, presenting a co-expression module enriched for the protein ubiquitination process. The results overlapped with epigenetics and phenotype-associated data, suggesting that DASE SNPs are potentially linked to cis-regulatory variants affecting simultaneously the transcription and phenotype through chromatin state modulation.

RNA-Seq exploration of the influence of stress on meat quality in Spanish goats

Studies exploring the transcriptome of stress and its effects on meat quality are very limited, particularly in goats. Fifty-four male Spanish goats (8-mo old; BW = 29.7 ± 2.03 kg) were randomly subjected to one of three treatments (TRT; n = 18 goats/treatment): (1) transported for 180 min, (2) transported for 30 min, or (3) held in pens (control) to analyze the transcriptome of stress and meat quality in goats using RNA-seq technology. Blood samples were collected before and after treatment, and meat samples were collected after humane slaughter for stress hormone, meat quality (Longissimus dorsi), and transcriptomic analysis. Plasma epinephrine concentrations were higher (P < 0.01) in 180 min and 30 min groups compared to the control group; however, norepinephrine concentrations were not affected by the treatment. Muscle glycogen concentrations (15 min postmortem) were lower (P < 0.01) in both 30 min and 180 min groups compared to the control group. Calpastatin levels were higher (P < 0.01) in 180 min and 30 min groups than the control group. Warner-Bratzler shear force values of loin chops were the highest in the 180 min group (4 ± 0.15, kg), lowest in the control group (3.51 ± 0.10, kg), and intermediate in the 30 min group (3.78 ± 0.09, kg; P < 0.01) both at day 1 and day 6 aging time. Additionally, desmin levels of day 6 samples were lowest in the control group, highest in 180 min group, and intermediate in 30 min group (P < 0.05). RNA-seq results showed that a total of 10,633 genes were differentially expressed (5194 up regulated; 5439 down regulated) among all comparisons (blood and day 1 and day 6 muscle samples). Among these differentially expressed genes (DEGs), KLF9, AMPK, FOXO3, PTX3, GADD45, PTPN1, CASP7, MAPK4, HSPA12A, and JAK-STAT were probably associated with the effects of stress on skeletal muscle proteins and involved in biological process such as cellular response to corticosteroid stimulus, endoplasmic reticulum stress, insulin resistance, DNA repair, apoptosis, MAPK cascade and regulation of proteolysis. The KEGG analysis revealed that AMPK and JAK-SAT signaling pathways and autophagy were among the top 20 enriched pathways in our treatment comparisons. The results provide an understanding of the genes and pathways involved in stress responses and related changes in postmortem muscle metabolism and meat quality characteristics in goats.

A Specific microRNA Targets an Elongase of Very Long Chain Fatty Acids to Regulate Fatty Acid Composition and Mitochondrial Morphology of Skeletal Muscle Cells

Recently, miR-22 has been suggested to be an important microRNA (miRNA) affecting meat quality. Studies have shown that muscle fatty acid composition and mitochondrial function are closely related to meat quality. The regulatory mechanism of miR-22 on skeletal muscle fatty acid composition and mitochondrial function is not well characterized. Therefore, we aimed to explore the effects of miR-22 on fatty acid composition and mitochondrial function in C2C12 cells. Here, it demonstrate that elevated expression of miR-22 significantly repressed fatty acid elongation and mitochondrial morphology in C2C12 myoblasts, while the knockdown of miR-22 showed opposite results. Furthermore, miR-22 targets the elongase of very long chain fatty acids 6 (ELOVL6) and represses its expression in muscle cells. Knockdown of ELOVL6 mimicked the effect of miR-22 on fatty acid composition and mitochondrial function, while overexpression of ELOVL6 restored the effects of miR-22. These findings indicate that miR-22 downregulates the elongation of fatty acids and mitochondrial morphology by inhibiting ELOVL6 expression in muscle cells, which may provide some useful information for controlling muscle lipid accumulation and mitochondrial function in livestock in the future.

Network analysis reveals different hub genes and molecular pathways for pig in vitro fertilized early embryos and parthenogenotes

Maternal-to-zygotic transition (MZT) occurs when maternal transcripts decay and zygotic genome is activated gradually at early stage of embryo development. Previously, single cell RNA-seq (scRNA-seq) has helped us to uncover the MZT-associated mRNA dynamics of in vitro produced pig early embryos. Here, to further investigate functional modules and hub genes associated with MZT process, the weighted gene-coexpression network analysis (WGCNA) was performed on our previously generated 45 scRNA-seq datasets. For the in vitro fertilized embryo (IVF) group, 5 significant modules were identified (midnightblue/black/red and blue/brown modules, positively correlated with 1-cell (IVF1) and 8-cell (IVF8), respectively), containing genes mainly enriched in signaling pathways such as Wnt, regulation of RNA transcription, fatty acid metabolic process, poly(A) RNA binding and lysosome. For the parthenogenetically activated embryo (PA) group, 9 significant modules were identified (black/purple/red, brown/turquoise/yellow, and magenta/blue/green modules, positively correlated with MII oocytes, 1-cell (PA1), and 8-cell (PA8), respectively), mainly enriched in extracellular exosome, poly(A) RNA binding, mitochondrion, transcription factor activity. Moreover, some of identified hub genes within 3 IVF and 9 PA significant modules, including ADCY2, DHX34, KDM4A, GDF10, ABCC10, PAFAH2, HEXIM2, COQ9, DCAF11, SGK1, ESRRB etc., have been reported to play vital roles in different biological processes. Our findings provide information and resources for subsequent in-depth study on the regulation and function of MZT in pig embryos.

Genome-Wide Associative Study of Phenotypic Parameters of the 3D Body Model of Aberdeen Angus Cattle with Multiple Depth Cameras

Simple Summary

This article aims to develop a new approach to the lifetime evaluation of cattle by 3-D visualization of economic-biological and genetic features. The following indicators were selected as phenotypic features: chest width and chest girth retrieved by 3-D model and meat output on the bones. Correlation analysis showed a reliable positive relationship between chest width and meat output on the bones, which can potentially be used for lifetime evaluation of meat productivity of animals. Genome-wide associations analysis revealed the following potential loci of quantitative traits on cattle chromosomes for chest width, chest girth, and meat output on bones.

Abstract

In beef cattle breeding, genome-wide association studies (GWAS) using single nucleotide polymorphisms (SNPs) arrays can reveal many loci of various production traits, such as growth, productivity, and meat quality. With the development of genome sequencing technologies, new opportunities are opening up for more accurate identification of areas associated with these traits. This article aims to develop a novel approach to the lifetime evaluation of cattle by 3-D visualization of economic-biological and genetic features. The purpose of this study was to identify significant variants underlying differences in the qualitative characteristics of meat, using imputed data on the sequence of the entire genome. Samples of biomaterial of young Aberdeen-Angus breed cattle (n = 96) were the material for carrying out genome-wide SNP genotyping. Genotyping was performed using a high-density DNA chip Bovine GPU HD BeadChip (Illumina Inc., San Diego, CA, USA), containing ~150 thousand SNPs. The following indicators were selected as phenotypic features: chest width and chest girth retrieved by 3-D model and meat output on the bones. Correlation analysis showed a reliable positive relationship between chest width and meat output on the bones, which can potentially be used for lifetime evaluation of meat productivity of animals.

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.

DNA methylation may affect beef tenderness through signal transduction in Bos indicus

BACKGROUND

Beef tenderness is a complex trait of economic importance for the beef industry. Understanding the epigenetic mechanisms underlying this trait may help improve the accuracy of breeding programs. However, little is known about epigenetic effects on Bos taurus muscle and their implications in tenderness, and no studies have been conducted in Bos indicus.

RESULTS

Comparing methylation profile of Bos indicus skeletal muscle with contrasting beef tenderness at 14 days after slaughter, we identified differentially methylated cytosines and regions associated with this trait. Interestingly, muscle that became tender beef had higher levels of hypermethylation compared to the tough group. Enrichment analysis of predicted target genes suggested that differences in methylation between tender and tough beef may affect signal transduction pathways, among which G protein signaling was a key pathway. In addition, different methylation levels were found associated with expression levels of GNAS, PDE4B, EPCAM and EBF3 genes. The differentially methylated elements correlated with EBF3 and GNAS genes overlapped CpG islands and regulatory elements. GNAS, a complex imprinted gene, has a key role on G protein signaling pathways. Moreover, both G protein signaling pathway and the EBF3 gene regulate muscle homeostasis, relaxation, and muscle cell-specificity.

CONCLUSIONS

We present differentially methylated loci that may be of interest to decipher the epigenetic mechanisms affecting tenderness. Supported by the previous knowledge about regulatory elements and gene function, the methylation data suggests EBF3 and GNAS as potential candidate genes and G protein signaling as potential candidate pathway associated with beef tenderness via methylation.

Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle

Background

A comprehensive analysis of gene expression profiling across tissues can provide necessary information for an in-depth understanding of their biological functions. We performed a large-scale gene expression analysis and generated a high-resolution atlas of the transcriptome in beef cattle.

Results

Our transcriptome atlas was generated from 135 bovine tissues in adult beef cattle, covering 51 tissue types of major organ systems (e.g., muscular system, digestive system, immune system, reproductive system). Approximately 94.76% of sequencing reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We detected a total of 60,488 transcripts, and 32% of them were not reported before. We identified 2654 housekeeping genes (HKGs) and 477 tissue-specific genes (TSGs) across tissues. Using weighted gene co-expression network analysis, we obtained 24 modules with 237 hub genes (HUBGs). Functional enrichment analysis showed that HKGs mainly maintain the basic biological activities of cells, while TSGs were involved in tissue differentiation and specific physiological processes. HKGs in bovine tissues were more conserved in terms of expression pattern as compared to TSGs and HUBGs among multiple species. Finally, we obtained a subset of tissue-specific differentially expressed genes (DEGs) between beef and dairy cattle and several functional pathways, which may be involved in production and health traits.

Conclusions

We generated a large-scale gene expression atlas across the major tissues in beef cattle, providing valuable information for enhancing genome assembly and annotation. HKGs, TSGs, and HUBGs further contribute to better understanding the biology and evolution of multiple tissues in cattle. DEGs between beef and dairy cattle also fill in the knowledge gaps about differential transcriptome regulation of bovine tissues underlying economically important traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01269-4.

Characterization of novel lncRNA muscle expression profiles associated with meat quality in beef cattle

The aim of this study was to identify novel lncRNA differentially expressed (DE) between divergent animals for beef tenderness and marbling traits in Nellore cattle. Longissimus thoracis muscle samples from the 20 most extreme bulls (of 80 bulls set) for tenderness, tender (n = 10) and tough (n = 10) groups, and marbling trait, high (n = 10) and low (n = 10) groups were used to perform transcriptomic analysis using RNA‐Sequencing. For tenderness, 29 lncRNA were DE (p‐value ≤ 0.01) in tough beef animals in relation to tender beef animals. We observed that genic lncRNAs, for example, lncRNA_595.1, were overlapping exonic part of the PICK gene, while lncRNA_3097.2 and lncRNA_3129.5 overlapped intronic part of the genes GADL1 and PSMD6. The lncRNA associated with PICK1, GADL1, and PMD6 genes were enriched in the pathways associated with the ionotropic glutamate receptor, gamma‐aminobutyric acid synthesis, and the ubiquitin–proteasome pathway. For marbling, 50 lncRNA were DE (p‐value ≤ 0.01) in high marbling group compared with low marbling animals. The genic lncRNAs, such as lncRNA_3191.1, were overlapped exonic part of the ITGAL gene, and the lncRNA_512.1, lncRNA_3721.1, and lncRNA_41.4 overlapped intronic parts of the KRAS and MASP1 genes. The KRAS and ITGAL genes were enriched in pathways associated with integrin signaling, which is involved in intracellular signals in response to the extracellular matrix, including cell form, mobility, and mediates progression through the cell cycle. In addition, the lincRNAs identified to marbling trait were associated with several genes related to calcium binding, muscle hypertrophy, skeletal muscle, lipase, and oxidative stress response pathways that seem to play a role important in the physiological processes related to meat quality. These findings bring new insights to better understand the biology mechanisms involved in the gene regulation of these traits, which will be valuable for a further investigation of the interactions between lncRNA and mRNAs, and of how these interactions may affect meat quality traits.

Genome-wide association study for meat tenderness in beef cattle identifies patterns of the genetic contribution in different post-mortem stages

The beef tenderization process during the post-mortem period is one of the most important sensorial attributes and it is well-established. The aim of this study was to identify the genetic contribution pattern to meat tenderness at 7-(LMD7), 14-(LMD14), and 21-(LMD21) days post-mortem. The heritabilities for LMD7 (0.194), LMD14 (0.142) and LMD21 (0.048) are well established in the population evaluated here. However, its genetic contribution in terms of genomic candidate regions is still poorly understood. Tenderness was measured in the Longissiums thoracis using Warner-Bratzler shear force in the three post-mortem periods. A total of 4323 crossbred beef cattle were phenotyped and genotyped using the Illumina BovineSNP50K. The percentage of the total genetic variance was estimated using the weighted single-step genomic best linear unbiased prediction method. The main candidate windows for LMD7 were associated with proteolysis of myofibrillar structures and the weakening endomysium and perimysium. Candidate windows for LMD14 and LMD21 were mapped in bovine QTLs for body composition, height and growth. Results presented herein highlight, the largest contribution of proteolysis related processes before 14-days post-mortem and body composition characteristics in later stages for meat tenderness.

Genes and proteins associated with ribeye area and meat tenderness in a commercial Nellore cattle population

Despite several studies on genetic markers and differently expressed genes related to ribeye area (REA) and tenderness traits in beef cattle, there is divergence in the results regarding the genes associated with these traits. Thirteen genes that had been associated or have biological functions that may influence such phenotypes were included in this study. A total of five genes for REA (IGF-1, IGF-2, MSTN, NEDD4, and UBE4A) and eight genes for meat tenderness (CAPN1, CAPN2, CAST, HSPB1, DNAJA1, FABP4, SCD, and PRKAG3) were selected from previously studies in beef cattle. Genes and its respective proteins expression were validated in a commercial population of Nellore cattle using quantitative real-time PCR (RT-qPCR) and advanced mass spectrometry (LC / MS-MS) techniques, respectively. MSTN gene was upregulated in animals with low REA. CAPN1, CAPN2, CAST, HSPB1, and DNAJA1 genes were upregulated in animals with tougher meat. The proteins translated by these genes were not differentially expressed. Our results could confirm the potential of some studied genes as biomarkers for carcass and meat quality in Nellore cattle.

Muscle allele-specific expression QTLs may affect meat quality traits in Bos indicus

Single nucleotide polymorphisms (SNPs) located in transcript sequences showing allele-specific expression (ASE SNPs) were previously identified in the Longissimus thoracis muscle of a Nelore (Bos indicus) population consisting of 190 steers. Given that the allele-specific expression pattern may result from cis-regulatory SNPs, called allele-specific expression quantitative trait loci (aseQTLs), in this study, we searched for aseQTLs in a window of 1 Mb upstream and downstream from each ASE SNP. After this initial analysis, aiming to investigate variants with a potential regulatory role, we further screened our aseQTL data for sequence similarity with transcription factor binding sites and microRNA (miRNA) binding sites. These aseQTLs were overlapped with methylation data from reduced representation bisulfite sequencing (RRBS) obtained from 12 animals of the same population. We identified 1134 aseQTLs associated with 126 different ASE SNPs. For 215 aseQTLs, one allele potentially affected the affinity of a muscle-expressed transcription factor to its binding site. 162 aseQTLs were predicted to affect 149 miRNA binding sites, from which 114 miRNAs were expressed in muscle. Also, 16 aseQTLs were methylated in our population. Integration of aseQTL with GWAS data revealed enrichment for traits such as meat tenderness, ribeye area, and intramuscular fat . To our knowledge, this is the first report of aseQTLs identification in bovine muscle. Our findings indicate that various cis-regulatory and epigenetic mechanisms can affect multiple variants to modulate the allelic expression. Some of the potential regulatory variants described here were associated with the expression pattern of genes related to interesting phenotypes for livestock. Thus, these variants might be useful for the comprehension of the genetic control of these phenotypes.

USP2-Related Cellular Signaling and Consequent Pathophysiological Outcomes

Ubiquitin specific protease (USP) 2 is a multifunctional deubiquitinating enzyme. USP2 modulates cell cycle progression, and therefore carcinogenesis, via the deubiquitination of cyclins and Aurora-A. Other tumorigenic molecules, including epidermal growth factor and fatty acid synthase, are also targets for USP2. USP2 additionally prevents p53 signaling. On the other hand, USP2 functions as a key component of the CLOCK/BMAL1 complex and participates in rhythmic gene expression in the suprachiasmatic nucleus and liver. USP2 variants influence energy metabolism by controlling hepatic gluconeogenesis, hepatic cholesterol uptake, adipose tissue inflammation, and subsequent systemic insulin sensitivity. USP2 also has the potential to promote surface expression of ion channels in renal and intestinal epithelial cells. In addition to modifying the production of cytokines in immune cells, USP2 also modulates the signaling molecules that are involved in cytokine signaling in the target cells. Usp2 knockout mice exhibit changes in locomotion and male fertility, which suggest roles for USP2 in the central nervous system and male genital tract, respectively. In this review, we summarize the cellular events with USP2 contributions and list the signaling molecules that are upstream or downstream of USP2. Additionally, we describe phenotypic differences found in the in vitro and in vivo experimental models.

Differential expression of miRNAs in skeletal muscles of Indian sheep with diverse carcass and muscle traits

The study presents the miRNA profiles of two Indian sheep populations with divergent carcass and muscle traits. The RNA sequencing of longissimus thoracis muscles from the two populations revealed a total of 400 known miRNAs. Myomirs or miRNAs specific to skeletal muscles identified in our data included oar-miR-1, oar-miR-133b, oar-miR-206 and oar-miR-486. Comparison of the two populations led to identification of 100 differentially expressed miRNAs (p < 0.05). A total of 45 miRNAs exhibited a log₂ fold change of ≥ ( ±) 3.0. Gene Ontology analysis revealed cell proliferation, epithelial to mesenchymal transition, apoptosis, immune response and cell differentiation as the most significant functions of the differentially expressed miRNAs. The differential expression of some miRNAs was validated by qRT-PCR analysis. Enriched pathways included metabolism of proteins and lipids, PI3K-Akt, EGFR and cellular response to stress. The microRNA-gene interaction network revealed miR-21, miR-155, miR-143, miR-221 and miR-23a as the nodal miRNAs, with multiple targets. MicroRNA-21 formed the focal point of the network with 42 interactions. The hub miRNAs identified in our study form putative regulatory candidates for future research on meat quality traits in Indian sheep. Our results provide insight into the biological pathways and regulatory molecules implicated in muscling traits of sheep.

The Profiling of DNA Methylation and Its Regulation on Divergent Tenderness in Angus Beef Cattle

Beef is an essential food source in the world. Beef quality, especially tenderness, has a significant impact on consumer satisfaction and industry profit. Many types of research to date have focused on the exploration of physiological and developmental mechanisms of beef tenderness. Still, the role and impact of DNA methylation status on beef tenderness have yet to be elucidated. In this study, we exhaustively analyzed the DNA methylation status in divergent tenderness observed in Angus beef. We characterized the methylation profiles related to beef tenderness and explored methylation distributions on the whole genome. As a result, differentially methylated regions (DMRs) associated with tenderness and toughness of beef were identified. Importantly, we annotated these DMRs on the bovine genome and explored bio-pathways of underlying genes and methylation biomarkers in beef quality. Specifically, we observed that the ATP binding cassette subfamily and myosin-related genes were highly methylated gene sets, and generation of neurons, regulation of GTPase activity, ion transport and anion transport, etc., were the significant pathways related with beef tenderness. Moreover, we explored the relationship between DNA methylation and gene expression in DMRs. Some methylated genes were identified as candidate biomarkers for beef tenderness. These results provide not only novel epigenetic information associated with beef quality but offer more significant insights into meat science, which will further help us explore the mechanism of muscle biology.

The application of omics in ruminant production: a review in the tropical and sub-tropical animal production context

The demand for animal products (e.g. dairy and beef) in tropical regions is expected to increase in parallel with the public demand for sustainable practices, due to factors such as population growth and climate change. The necessity to increase animal production output must be achieved with better management and production technologies. For this to happen, novel research methodologies, animal selection and postgenomic tools play a pivotal role. Indeed, improving breeder selection programs, the quality of meat and dairy products as well as animal health will contribute to higher sustainability and productivity. This would surely benefit regions where resource quality and quantity are increasingly unstable, and research is still very incipient, which is the case of many regions in the tropics. The purpose of this review is to demonstrate how omics-based approaches play a major role in animal science, particularly concerning ruminant production systems and research associated to the tropics and developing countries. SIGNIFICANCE: Environmental conditions in the tropics make livestock production harder, compared to temperate regions. Due to global warming, the sustainability of livestock production will become increasingly problematic. The use of novel omics technologies could generate useful information to understand adaptation mechanisms of resilient breeds and/or species. The application of omics to tropical animal production is still residual in the currently available literature. With this review, we aim to summarize the most notable results in the field whilst encouraging further research to deal with the future challenges that animal production in the tropics will need to face.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

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.

Signatures of positive selection underlying beef production traits in Korean cattle breeds

The difference in the breeding programs and population history may have diversely shaped the genomes of Korean native cattle breeds. In the absence of phenotypic data, comparisons of breeds that have been subjected to different selective pressures can aid to identify genomic regions and genes controlling qualitative and complex traits. In this study to decipher genetic variation and identify evidence of divergent selection, 3 Korean cattle breeds were genotyped using the recently developed high-density GeneSeek Genomic Profiler F250 (GGP-F250) array. The three Korean cattle breeds clustered according to their coat color phenotypes and breeding programs. The Heugu breed reliably showed smaller effective population size at all generations considered. Across the autosomal chromosomes, 113 and 83 annotated genes were identified from Hanwoo-Chikso and Hanwoo-Heugu comparisons, respectively of which 16 genes were shared between the two pairwise comparisons. The most important signals of selection were detected on bovine chromosomes 14 (24.39–25.13 Mb) and 18 (13.34–15.07 Mb), containing genes related to body size, and coat color (XKR4, LYN, PLAG1, SDR16C5, TMEM68, CDH15, MC1R, and GALNS). Some of the candidate genes are also associated with meat quality traits (ACSF3, EIF2B1, BANP, APCDD1, and GALM) and harbor quantitative trait locus (QTL) for beef production traits. Further functional analysis revealed that the candidate genes (DBI, ACSF3, HINT2, GBA2, AGPAT5, SCAP, ELP6, APOB, and RBL1) were involved in gene ontology (GO) terms relevant to meat quality including fatty acid oxidation, biosynthesis, and lipid storage. Candidate genes previously known to affect beef production and quality traits could be used in the beef cattle selection strategies.

Common Gene Modules Identified for Chicken Adiposity by Network Construction and Comparison

Excessive fat deposition can cause chicken health problem, and affect production efficiency by causing great economic losses to the industry. However, the molecular underpinnings of the complex adiposity trait remain elusive. In the current study, we constructed and compared the gene co-expression networks on four transcriptome profiling datasets, from two chicken lines under divergent selection for abdominal fat contents, in an attempt to dissect network compositions underlying adipose tissue growth and development. After functional enrichment analysis, nine network modules important to adipogenesis were discovered to be involved in lipid metabolism, PPAR and insulin signaling pathways, and contained hub genes related to adipogenesis, cell cycle, inflammation, and protein synthesis. Moreover, after additional functional annotation and network module comparisons, common sub-modules of similar functionality for chicken fat deposition were identified for different chicken lines, apart from modules specific to each chicken line. We further validated the lysosome pathway, and found TFEB and its downstream target genes showed similar expression patterns along with chicken preadipocyte differentiation. Our findings could provide novel insights into the genetic basis of complex adiposity traits, as well as human obesity and related metabolic diseases.

Proteome alterations associated with the oleic acid and cis-9, trans-11 conjugated linoleic acid content in bovine skeletal muscle

Oleic acid (OA) and cis-9, trans-11 conjugated linoleic acid (c9t11-CLA) are fatty acids found in beef with beneficial effects in human health. This study investigated differentially abundant proteins (DAPs) in skeletal muscle of bovines with extreme values of OA, and c9t11-CLA. For each one of the fatty acids, twenty muscle samples were divided into two groups (N = 10_High; N = 10_Low) and analyzed by high definition mass spectrometry. We identified 103 and 133 DAPs between the groups for each fatty acid. We found 64 and 45 up-regulated and 39 and 68 down-regulated proteins for OA and c9t11-CLA, respectively. Comparative analysis between proteomic and transcriptomic data revealed eight and ten genes with a consistent between mRNA expression levels and protein abundance for OA and c9t11-CLA, respectively. Unconventional myosin-Id (MYO1D), mineralocorticoid receptor (NR3C2), geranylgeranyl transferase type-2 subunit-alpha (RABGGTA), and uveal autoantigen with coiled-coil domains and ankyrin repeats (UACA) were found as putative candidate proteins for OA content. Fatty acid synthase (FASN), tubulin alpha-4A chain (TUBA4A), vinculin (VCL), NADH dehydrogenase 1 alpha subcomplex 5 (NDUFA5), and prefoldin subunit 6 (PFDN6) for c9t11-CLA. Our findings contribute to a deeper understanding of the molecular mechanisms behind the regulation of the OA and c9t11-CLA content in cattle skeletal muscle. SIGNIFICANCE: Questions about the association between meat intake and disease incidence in humans has driven animal scientist to pursue a better understanding of the biological processes associated with differences in the intramuscular fat composition. The beneficial effects of oleic acid and conjugated linoleic acid in human health have been demonstrated by improving the immune system and preventing atherosclerosis, different types of cancers, hypertension, and diabetes. Previous genome-wide association and gene expression studies identified genomic regions and differentially expressed genes associated with the fatty acid profile in skeletal muscle. In this work, differences were evaluated at the protein level. The use of a label-free quantitative proteomic approach, compared with muscle transcriptome results obtained by RNA-sequencing, allowed us to earn new insights into the variability in fatty acid deposition in skeletal muscle of farm animals. This study opens new avenues to explore the effect of the fatty acids in the skeletal muscle of livestock animals, which is associated with nutritional values of the meat, and perhaps to understand the mechanisms correlated with metabolic diseases in other species.

Cross talk between mineral metabolism and meat quality: a systems biology overview

Meat quality has an inherent complexity because of the multiple interrelated causative factors and layers of feedback regulation. Understanding the key factors and their interactions has been challenging, despite the availability of remarkable high-throughput tools and techniques that have provided insights on muscle metabolism and the genetic basis of meat quality. Likewise, we have deepened our knowledge about mineral metabolism and its role in cell functioning. Regardless of these facts, complex traits like mineral content and meat quality have been studied under reductionist approaches. However, as these phenotypes arise from complex interactions among different biological layers (genome, transcriptome, proteome, epigenome, etc.), along with environmental effects, a holistic view and systemic-level understanding of the genetic basis of complex phenotypes are in demand. Based on the state of the art, we addressed some of the questions regarding the interdependence of meat quality traits and mineral content. Furthermore, we sought to highlight potential regulatory mechanisms arising from the genes, miRNAs, and mineral interactions, as well as the pathways modulated by this interplay affecting muscle, mineral metabolism, and meat quality. By answering these questions, we did not intend to give an exhaustive review but to identify the key biological points, the challenges, and benefits of integrative genomic approaches.

Co-Expression Networks Reveal Potential Regulatory Roles of miRNAs in Fatty Acid Composition of Nelore Cattle

Fatty acid (FA) content affects the sensorial and nutritional value of meat and plays a significant role in biological processes such as adipogenesis and immune response. It is well known that, in beef, the main FAs associated with these biological processes are oleic acid (C18:1 cis9, OA) and conjugated linoleic acid (CLA-c9t11), which may have beneficial effects on metabolic diseases such as type 2 diabetes and obesity. Here, we performed differential expression and co-expression analyses, weighted gene co-expression network analysis (WGCNA) and partial correlation with information theory (PCIT), to uncover the complex interactions between miRNAs and mRNAs expressed in skeletal muscle associated with FA content. miRNA and mRNA expression data were obtained from skeletal muscle of Nelore cattle that had extreme genomic breeding values for OA and CLA. Insulin and MAPK signaling pathways were identified by WGCNA as central pathways associated with both of these fatty acids. Co-expression network analysis identified bta-miR-33a/b, bta-miR-100, bta-miR-204, bta-miR-365-5p, bta-miR-660, bta-miR-411a, bta-miR-136, bta-miR-30-5p, bta-miR-146b, bta-let-7a-5p, bta-let-7f, bta-let-7, bta-miR 339, bta-miR-10b, bta-miR 486, and the genes ACTA1 and ALDOA as potential regulators of fatty acid synthesis. This study provides evidence and insights into the molecular mechanisms and potential target genes involved in fatty acid content differences in Nelore beef cattle, revealing new candidate pathways of phenotype modulation that could positively benefit beef production and human consumption.

Detection of Co-expressed Pathway Modules Associated With Mineral Concentration and Meat Quality in Nelore Cattle

Meat quality is a complex trait that is influenced by genetic and environmental factors, which includes mineral concentration. However, the association between mineral concentration and meat quality, and the specific molecular pathways underlying this association, are not well explored. We therefore analyzed gene expression as measured with RNA-seq in Longissimus thoracis muscle of 194 Nelore steers for association with three meat quality traits (intramuscular fat, meat pH, and tenderness) and the concentration of 13 minerals (Ca, Cr, Co, Cu, Fe, K, Mg, Mn, Na, P, S, Se, and Zn). We identified seven sets of co-expressed genes (modules) associated with at least two traits, which indicates that common pathways influence these traits. From pathway analysis of module hub genes, we further found an over-representation for energy and protein metabolism (AMPK and mTOR signaling pathways) in addition to muscle growth, and protein turnover pathways. Among the identified hub genes FASN, ELOV5, and PDE3B are involved with lipid metabolism and were affected by previously identified eQTLs associated to fat deposition. The reported hub genes and over-represented pathways provide evidence of interplay among gene expression, mineral concentration, and meat quality traits. Future studies investigating the effect of different levels of mineral supplementation in the gene expression and meat quality traits could help us to elucidate the regulatory mechanism by which the genes/pathways are affected.

MiRNAs differentially expressed in skeletal muscle of animals with divergent estimated breeding values for beef tenderness

Background

MicroRNAs (miRNAs) are small noncoding RNAs of approximately 22 nucleotides, highly conserved among species, which modulate gene expression by cleaving messenger RNA target or inhibiting translation. MiRNAs are involved in the regulation of many processes including cell proliferation, differentiation, neurogenesis, angiogenesis, and apoptosis. Beef tenderness is an organoleptic characteristic of great influence in the acceptance of meat by consumers. Previous studies have shown that collagen level, marbling, apoptosis and proteolysis are among the many factors that affect beef tenderness. Considering that miRNAs can modulate gene expression, this study was designed to identify differentially expressed miRNAs that could be modulating biological processes involved with beef tenderness.

Results

Deep sequence analysis of miRNA libraries from longissimus thoracis muscle allowed the identification of 42 novel and 308 known miRNAs. Among the known miRNAs, seven were specifically expressed in skeletal muscle. Differential expression analysis between animals with high (H) and low (L) estimated breeding values for shear force (EBVSF) revealed bta-mir-182 and bta-mir-183 are up-regulated (q value < 0.05) in animals with L EBVSF, and bta-mir-338 is up-regulated in animals with H EBVSF. The number of bovine predicted targets for bta-mir-182, bta-mir-183 and bta-mir-338 were 811, 281 and 222, respectively, which correspond to 1204 unique target genes. Among these, four of them, MEF2C, MAP3K2, MTDH and TNRC6B were common targets of the three differentially expressed miRNAs. The functional analysis identified important pathways related to tenderness such as apoptosis and the calpain–calpastatin system.

Conclusion

The results obtained indicate the importance of miRNAs in the regulatory mechanisms that influence muscle proteolysis and meat tenderness and contribute to our better understanding of the role of miRNAs in biological processes associated with beef tenderness.

Electronic supplementary material

The online version of this article (10.1186/s12867-018-0118-3) contains supplementary material, which is available to authorized users.