Differentially expressed mRNAs, proteins and miRNAs associated to energy metabolism in skeletal muscle of beef cattle identified for low and high residual feed intake

Plasma proteomic analysis reveals key pathways associated with divergent residual body weight gain phenotype in beef steers

We utilized plasma proteomics profiling to explore metabolic pathways and key proteins associated with divergent residual body weight gain (RADG) phenotype in crossbred (Angus × Hereford) beef steers. A group of 108 crossbred growing beef steers (average BW = 282.87 ± 30 kg; age = 253 ± 28 days) were fed a high-forage total mixed ration for 49 days in five dry lot pens (20-22 beef steers per pen), each equipped with two GrowSafe8000 intake nodes to determine their RADG phenotype. After RADG identification, blood samples were collected from the beef steers with the highest RADG (most efficient; n = 15; 0.76 kg/d) and lowest RADG (least efficient; n = 15; -0.65 kg/d). Plasma proteomics analysis was conducted on all plasma samples using a nano LC-MS/MS platform. Proteins with FC ≥ 1.2 and false-discovery rate-adjusted p-values (FDR) ≤ 0.05 were considered significantly differentially abundant. The analysis identified 435 proteins, with 59 differentially abundant proteins (DAPs) between positive and negative-RADG beef steers. Plasma abundance of 38 proteins, such as macrophage stimulating 1 and peptidase D was upregulated (FC ≥ 1.2, FDR ≤ 0.05) in positive-RADG beef steers, while 21 proteins, including fibronectin and ALB protein were greater (FC < 1.2, FDR ≤ 0.05) in negative-RADG beef steers. The results of the Gene Ontology (GO) analysis of all the DAPs showed enrichment of pathways such as metabolic processes, biological regulation, and catalytic activity in positive-RADG beef steers. Results of the EuKaryotic Orthologous Groups (KOG) analysis revealed increased abundance of DAPs involved in energy production and conversion, amino acid transport and metabolism, and lipid transport and metabolism in positive-RADG beef steers. The results of this study revealed key metabolic pathways and proteins associated with divergent RADG phenotype in beef cattle which give more insight into the biological basis of feed efficiency in crossbred beef cattle.

Compilations and updates on residual feed intake in sheep

The increasing global demand for food and the strong effect of climate change have forced animal science to advance regarding new methods of selection in search of more efficient animals in production systems. Feed consumption represents more than 70% of the costs of sheep farms, and more efficient animals can increase the farmers' profitability. One of the main measures of feed efficiency is estimated residual feed intake (RFI), created in 1963 by Robert Koch for estimation in cattle and later adapted for sheep. Animals with negative RFI values (RFI-) are more efficient than animals with positive values (RFI+), with influence on the variables of performance, carcass quality and production of enteric gases. The RFI is the most common and accepted metric of the feed efficiency trait for genetic selection, since it is independent of growth traits, unlike the feed conversion ratio. The purpose of this review article was to present updated literature information on the relationship of RFI estimates with performance measures, molecular markers, greenhouse gas production and feed efficiency, the technical aspects and physiological basis of metabolic in sheep.

Polymorphisms of PLIN1 and MOGAT1 genes and their association with feed efficiency in Hu sheep

Feed cost accounts for a high proportion of sheep production, and improving sheep's utilization of feed will reduce production costs and improve economic benefits. The purpose of this study was to investigate the expression characteristics of PLIN1 and MOGAT1 genes and the relationship between their polymorphisms and feed efficiency traits in Hu sheep, and to find molecular Genetic marker that can be used in breeding. The expression levels of PLIN1 and MOGAT1 genes in various tissues were determined using quantitative real-time PCR (qRT-PCR). The results showed that PLIN1 and MOGAT1 genes were widely expressed in heart, liver, spleen, lungs, kidneys, rumen, duodenum, muscle, lymph, and tail fat. The PLIN1 gene had the highest expression level in in the tail fat compared to the other nine tissues. The expression levels of MOGAT1 gene in liver, tail fat, lung and heart was significantly higher than in kidney, muscle and lymph. The expression level of MOGAT1 was lowest in muscle compared to the other tissues (heart, liver, spleen, lung, rumen and tail fat). We recorded the body weight (BW80 and BW180) and feed intake (FI) information of 985 male Hu sheep at 80 and 180 days of age, and calculated the daily average feed intake (ADFI), average daily gain (ADG), and feed conversion rate (FCR) from 80 to 180 days of age. Two intronic mutations, g.18517910 A > G and g.224856118 G > C, were identified in PLIN1 and MOGAT1 genes by PCR amplification and Sanger sequencing. MassARRAY ® SNP detection technology was used to genotype the DNA of 985 Hu sheep and analyze its association with feed efficiency traits. The results showed that the SNP g.18517910 A > G was significantly associated with BW80, BW180, FI, ADFI and FCR (P < 0.05), while SNP g.2248561118 G > C was significantly associated with FCR (P < 0.05). Meanwhile, significant differences were also observed in different combinations of genotypes (P < 0.05). Therefore, these two polymorphic loci can serve as candidate molecular markers for improving feed utilization efficiency in Hu sheep.

Association of plasma miRNAs with early life performance and aging in dairy cattle

Early life performance traits in dairy cattle can have important influences on lifetime productivity. Poor health and fertility are of great economical and animal welfare concern. Circulating miRNAs have been linked to several livestock traits, including resistance to infection, fertility, and muscle development. This study aimed to identify circulating miRNAs associated with early life performance traits and aging in dairy cattle. Plasma samples from female calves (n = 12) identified retrospectively as differing in health, growth, and fertility outcomes prior to first calving were analyzed using PCR arrays detecting 378 miRNAs. Levels of 6 miRNAs differed significantly in calves with poor growth/fertility relative to controls (t-test: P<0.05). Additionally, general(ized) (non)linear mixed models identified 1 miRNA associated with average daily gain until weaning, 22 with live bodyweight at one year of age, 47 with age at first service, and 19 with number of infections before first calving. Out of 85 distinct miRNAs that were associated with at least one animal trait, 9 miRNAs were validated by RT-qPCR in a larger cohort (n = 91 animals), which included longitudinal plasma samples (calf, heifer, first lactation cow). Significant associations (P<0.05) involving individual miRNAs or ratios between miRNAs and early-life performance traits were identified, but did not retain significance after multiple testing adjustment. However, levels of 8 plasma miRNAs (miR-126-3p, miR-127, miR-142-5p, miR-154b, miR-27b, miR-30c-5p, miR-34a, miR-363) changed significantly with age, most prominently during the calf-to-heifer transition. Comparative RT-qPCR analyses of these miRNAs across 19 calf tissues showed that most were ubiquitously expressed. Online database mining identified several pathways involved in metabolism and cell signaling as putative biological targets of these miRNAs. These results suggest that miR-126-3p, miR-127, miR-142-5p, miR-154b, miR-27b, miR-30c-5p, miR-34a, miR-363 are involved in regulating growth and development from birth to first lactation (~2 years old) and could provide useful biomarkers of aging in cattle.

Feed efficiency in dairy sheep: An insight from the milk transcriptome

Introduction

As higher feed efficiency in dairy ruminants means a higher capability to transform feed nutrients into milk and milk components, differences in feed efficiency are expected to be partly linked to changes in the physiology of the mammary glands. Therefore, this study aimed to determine the biological functions and key regulatory genes associated with feed efficiency in dairy sheep using the milk somatic cell transcriptome.

Material and methods

RNA-Seq data from high (H-FE, n = 8) and low (L-FE, n = 8) feed efficiency ewes were compared through differential expression analysis (DEA) and sparse Partial Least Square-Discriminant analysis (sPLS-DA).

Results

In the DEA, 79 genes were identified as differentially expressed between both conditions, while the sPLS-DA identified 261 predictive genes [variable importance in projection (VIP) &gt; 2] that discriminated H-FE and L-FE sheep.

Discussion

The DEA between sheep with divergent feed efficiency allowed the identification of genes associated with the immune system and stress in L-FE animals. In addition, the sPLS-DA approach revealed the importance of genes involved in cell division (e.g., KIF4A and PRC1) and cellular lipid metabolic process (e.g., LPL, SCD, GPAM, and ACOX3) for the H-FE sheep in the lactating mammary gland transcriptome. A set of discriminant genes, commonly identified by the two statistical approaches, was also detected, including some involved in cell proliferation (e.g., SESN2, KIF20A, or TOP2A) or encoding heat-shock proteins (HSPB1). These results provide novel insights into the biological basis of feed efficiency in dairy sheep, highlighting the informative potential of the mammary gland transcriptome as a target tissue and revealing the usefulness of combining univariate and multivariate analysis approaches to elucidate the molecular mechanisms controlling complex traits.

Identification and characterization of hypothalamic circular RNAs associated with bovine residual feed intake

Residual feed intake (RFI) is crucial economic indicator used for calculating the feed efficiency of growing beef cattle. circRNA plays an important biological role in gene transcriptional regulation, but little is known about its potential functional regulation underlying RFI phenotypic variation. As the core center of regulation of animal feeding, the hypothalamus is closely associated with RFI. Therefore, the present study aimed to identify the key genes and functional pathways contributing to variance in cattle RFI phenotypes using RNA sequencing from hypothalamic tissue samples, in order to gain insight into the potential regulatory role of circRNAs in bovine RFI phenotypic variation. Differentially expressed genes were detected by RNA sequencing for beef cattle in the high and low RFI groups, followed by GO, KEGG enrichment, and circRNA-miRNA co-expression network analysis. A total of 257 circRNAs were differentially expressed between the two groups, with 128 significantly upregulated and 129 significantly downregulated genes in H group compared to L group. Among them, 9 unique circRNAs were present in group L and 4 unique circRNAs were present in group H. GO and KEGG enrichment analysis of the source genes of the differentially expressed circRNAs revealed that they were mainly involved in metabolic processes, such as cellular metabolic processes, cellular macromolecular metabolic processes, and regulatory pathways related to nutrient metabolism, including protein and amino acid metabolism, as well as vitamin metabolism and pancreatic secretion associated with the animal feeding behavior. The circRNAs detected in this study were mostly novel, and have not been investigated directly to be associated with the RFI phenotype. Interestingly, most miRNAs of differentially expressed circRNAs predicted based on the circRNA-miRNA co-expression network analysis by using top 50 differentially expressed circRNAs and 13 unique circRNAs, have been reported to be related to animal RFIs, implying that circRNAs in bovine hypothalamic tissue may regulate phenotypic variation in RFI through miRNAs. The study results illustrate the complex biological functions of the hypothalamus in regulating feed efficiency and showing the potential role of circRNAs in the feeding behavior regulation of livestock, which would contributing to expanding the understanding of circRNA.

Meta-analysis across Nellore cattle populations identifies common metabolic mechanisms that regulate feed efficiency-related traits

Background

Feed efficiency (FE) related traits play a key role in the economy and sustainability of beef cattle production systems. The accurate knowledge of the physiologic background for FE-related traits can help the development of more efficient selection strategies for them. Hence, multi-trait weighted GWAS (MTwGWAS) and meta-analyze were used to find genomic regions associated with average daily gain (ADG), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), and residual feed intake (RFI). The FE-related traits and genomic information belong to two breeding programs that perform the FE test at different ages: post-weaning (1,024 animals IZ population) and post-yearling (918 animals for the QLT population).

Results

The meta-analyze MTwGWAS identified 14 genomic regions (-log10(p -value) > 5) regions mapped on BTA 1, 2, 3, 4, 7, 8, 11, 14, 15, 18, 21, and 29. These regions explained a large proportion of the total genetic variance for FE-related traits across-population ranging from 20% (FCR) to 36% (DMI) in the IZ population and from 22% (RFI) to 28% (ADG) in the QLT population. Relevant candidate genes within these regions (LIPE, LPL, IGF1R, IGF1, IGFBP5, IGF2, INS, INSR, LEPR, LEPROT, POMC, NPY, AGRP, TGFB1, GHSR, JAK1, LYN, MOS, PLAG1, CHCD7, LCAT, and PLA2G15) highlighted that the physiological mechanisms related to neuropeptides and the metabolic signals controlling the body's energy balance are responsible for leading to greater feed efficiency. Integrated meta-analysis results and functional pathway enrichment analysis highlighted the major effect of biological functions linked to energy, lipid metabolism, and hormone signaling that mediates the effects of peptide signals in the hypothalamus and whole-body energy homeostasis affecting the genetic control of FE-related traits in Nellore cattle.

Conclusions

Genes and pathways associated with common signals for feed efficiency-related traits provide better knowledge about regions with biological relevance in physiological mechanisms associated with differences in energy metabolism and hypothalamus signaling. These pleiotropic regions would support the selection for feed efficiency-related traits, incorporating and pondering causal variations assigning prior weights in genomic selection approaches.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08671-w.

Integrated analysis of expression profiles with meat quality traits in cattle

MicroRNAs (miRNAs) play a vital role in improving meat quality by binding to messenger RNAs (mRNAs). We performed an integrated analysis of miRNA and mRNA expression profiling between bulls and steers based on the differences in meat quality traits. Fat and fatty acids are the major phenotypic indices of meat quality traits to estimate between-group variance. In the present study, 90 differentially expressed mRNAs (DEGs) and 18 differentially expressed miRNAs (DEMs) were identified. Eighty-three potential DEG targets and 18 DEMs were used to structure a negative interaction network, and 75 matching target genes were shown in this network. Twenty-six target genes were designated as intersection genes, screened from 18 DEMs, and overlapped with the DEGs. Seventeen of these genes enriched to 19 terms involved in lipid metabolism. Subsequently, 13 DEGs and nine DEMs were validated using quantitative real-time PCR, and seven critical genes were selected to explore the influence of fat and fatty acids through hub genes and predict functional association. A dual-luciferase reporter and Western blot assays confirmed a predicted miRNA target (bta-miR-409a and PLIN5). These findings provide substantial evidence for molecular genetic controls and interaction among genes in cattle.

Differences in Liver microRNA profiling in pigs with low and high feed efficiency

Feed cost is the main factor affecting the economic benefits of pig industry. Improving the feed efficiency (FE) can reduce the feed cost and improve the economic benefits of pig breeding enterprises. Liver is a complex metabolic organ which affects the distribution of nutrients and regulates the efficiency of energy conversion from nutrients to muscle or fat, thereby affecting feed efficiency. MicroRNAs (miRNAs) are small non-coding RNAs that can regulate feed efficiency through the modulation of gene expression at the post-transcriptional level. In this study, we analyzed miRNA profiling of liver tissues in High-FE and Low-FE pigs for the purpose of identifying key miRNAs related to feed efficiency. A total 212~221 annotated porcine miRNAs and 136~281 novel miRNAs were identified in the pig liver. Among them, 188 annotated miRNAs were co-expressed in High-FE and Low-FE pigs. The 14 miRNAs were significantly differentially expressed (DE) in the livers of high-FE pigs and low-FE pigs, of which 5 were downregulated and 9 were upregulated. Kyoto Encyclopedia of Genes and Genomes analysis of liver DE miRNAs in high-FE pigs and low-FE pigs indicated that the target genes of DE miRNAs were significantly enriched in insulin signaling pathway, Gonadotropin-releasing hormone signaling pathway, and mammalian target of rapamycin signaling pathway. To verify the reliability of sequencing results, 5 DE miRNAs were randomly selected for quantitative reverse transcription-polymerase chain reaction (qRT-PCR). The qRT-PCR results of miRNAs were confirmed to be consistent with sequencing data. DE miRNA data indicated that liver-specific miRNAs synergistically acted with mRNAs to improve feed efficiency. The liver miRNAs expression analysis revealed the metabolic pathways by which the liver miRNAs regulate pig feed efficiency.

Emerging Roles of Non-Coding RNAs in the Feed Efficiency of Livestock Species

A global population of already more than seven billion people has led to an increased demand for food and water, and especially the demand for meat. Moreover, the cost of feed used in animal production has also increased dramatically, which requires animal breeders to find alternatives to reduce feed consumption. Understanding the biology underlying feed efficiency (FE) allows for a better selection of feed-efficient animals. Non-coding RNAs (ncRNAs), especially micro RNAs (miRNAs) and long non-coding RNAs (lncRNAs), play important roles in the regulation of bio-logical processes and disease development. The functions of ncRNAs in the biology of FE have emerged as they participate in the regulation of many genes and pathways related to the major FE indicators, such as residual feed intake and feed conversion ratio. This review provides the state of the art studies related to the ncRNAs associated with FE in livestock species. The contribution of ncRNAs to FE in the liver, muscle, and adipose tissues were summarized. The research gap of the function of ncRNAs in key processes for improved FE, such as the nutrition, heat stress, and gut–brain axis, was examined. Finally, the potential uses of ncRNAs for the improvement of FE were discussed.

Common and diet-specific metabolic pathways underlying residual feed intake in fattening Charolais yearling bulls

Residual feed intake (RFI) is one of the preferred traits for feed efficiency animal breeding. However, RFI measurement is expensive and time-consuming and animal ranking may depend on the nature of the diets. We aimed to explore RFI plasma biomarkers and to unravel the underlying metabolic pathways in yearling bulls fed either a corn-silage diet rich in starch (corn diet) or a grass-silage diet rich in fiber (grass diet). Forty-eight extreme RFI animals (Low-RFI, n = 24, versus High-RFI, n = 24, balanced per diet) were selected from a population of 364 Charolais bulls and their plasma was subjected to a targeted LC-MS metabolomic approach together with classical metabolite and hormonal plasma analyses. Greater lean body mass and nitrogen use efficiency, and lower protein turnover were identified as common mechanisms underlying RFI irrespective of the diet. On the other hand, greater adiposity and plasma concentrations of branched-chain amino acids (BCAA) together with lower insulin sensitivity in High-RFI animals were only observed with corn diet. Conversely, greater plasma concentrations of BCAA and total triglycerides, but similar insulin concentrations were noted in efficient RFI cattle with grass diet. Our data suggest that there are diet-specific mechanisms explaining RFI differences in fattening Charolais yearling bulls.

Dark-cutting beef: A brief review and an integromics meta-analysis at the proteome level to decipher the underlying pathways

Comprehensive characterization of the post-mortem muscle proteome defines a fundamental goal in meat proteomics. During the last decade, proteomics tools have been applied in the field of foodomics to help decipher factors underpinning meat quality variations and to enlighten us, through data-driven methods, on the underlying mechanisms leading to meat quality defects such as dark-cutting meat known also as dark, firm and dry (DFD) meat. In cattle, several proteomics studies have focused on the extent to which changes in the post-mortem muscle proteome relate to dark-cutting beef development. The present data-mining study firstly reviews proteomics studies which investigated dark-cutting beef, and secondly, gathers the protein biomarkers that differ between dark-cutting versus beef with normal-pH in a unique repertoire. A list of 130 proteins from eight eligible studies was curated and mined through bioinformatics for Gene Ontology annotations, molecular pathways enrichments, secretome analysis and biological pathways comparisons to normal beef color from a previous meta-analysis. The major biological pathways underpinning dark-cutting beef at the proteome level have been described and deeply discussed in this integromics study.

MicroRNA-mRNA Regulatory Networking Fine-Tunes Polyunsaturated Fatty Acid Synthesis and Metabolism in the Inner Mongolia Cashmere Goat

Fatty acid composition is an important aspect of meat quality in ruminants. Improving the beneficial fatty acid level in cashmere goat meat is important to its economic value. To investigate microRNAs (miRNAs) and mRNAs that regulate or coregulate polyunsaturated fatty acid (PUFA) synthesis and metabolism in the Inner Mongolia cashmere goat, we used longissimus dorsi muscle (WLM) and biceps femoris muscle (WBM) for transcript-level sequencing. RT-qPCR was used to evaluate the expression of mRNAs and miRNAs associated with PUFA synthesis and metabolism. The total PUFA content in the WBM was significantly higher than that in the WLM (P < 0.05). Our study is the first to systematically report miRNAs in cashmere goat meat. At the mRNA level, 20,375 genes were identified. ACSL1, CD36 and TECRL were at the center of a gene regulatory network and contributed significantly to the accumulation and metabolic regulation of fatty acids. At the miRNA level, 426 known miRNAs and 30 novel miRNAs were identified. KEGG analysis revealed that the miRNA target genes were involved mainly in the PPAR signaling pathway. The mRNA-miRNA coregulation analysis showed that ACSL1 was negatively targeted by nine miRNAs: chi-miR-10a-5p, chi-miR-10b-5p, chi-miR-130b-5p, chi-miR-15a-5p_R-1, chi-miR-15b-5p, chi-miR-16a-5p, chi-miR-16b-5p, chi-miR-181c-5p_R+1, and chi-miR-26b-5p. Finally, we speculated that the simultaneous silencing of ACSL1 by one or more of these nine miRNAs through PPAR signaling led to low ACSL1 expression in the WLM and, ultimately to high PUFA content in the WBM. Our study helps elucidate the metabolic regulation of fatty acids in Inner Mongolia cashmere goats.

Polymorphisms in ovine ME1 and CA1 genes and their association with feed efficiency in Hu sheep

Improving feed efficiency would increase profitability for producers. The objective of this study was to detect the expression levels of ME1 and CA1 and the polymorphisms of ME1 and CA1 associated with the feed conversion ratio (FCR) and residual feed intake (RFI) of Hu sheep by using qRT-PCR, pooled DNA sequencing and KASPar assay. The qRT-PCR results indicated that the expression levels of ME1 and CA1 were significantly higher in the liver tissues of low-RFI sheep than in those of the high-RFI sheep (p <.01). Association analysis demonstrated that the polymorphism ME1 g.453 C > T was significantly associated with FCR and RFI (p <.05). The polymorphism CA1 g.199 C > T had a significantly associated with FCR (p <.05) and no association with RFI (p >.05). Significant differences were observed between the combined genotypes and FCR and RFI at different measurement periods (p <.05). Thus, we propose the use of these two polymorphisms as new candidate molecular markers for improving feed efficiency in sheep populations.

Bovine hepatic miRNAome profiling and differential miRNA expression analyses between beef steers with divergent feed efficiency phenotypes

MicroRNAs (miRNAs) are small RNA molecules involved in regulation of multiple biological processes through modulating expression of their target genes. Here we employed RNAseq to profile liver tissue miRNAome of 60 steers from Angus, Charolais, and Kinsella Composite (KC) populations. Of these animals, 36 animals (n = 12 for each breed) were utilized to identify differentially expressed (DE) miRNAs between animals with high (n = 6) or low (n = 6) phenotypic values of residual feed intake (RFI), a common measurement of feed efficiency. At a threshold of fold-change > 1.5 and P-value < 0.05, we detected 12 (7 up- and 5 downregulated in low-RFI animals), 18 (12 up- and 6 downregulated), and 13 (8 up- and 5 downregulated) DE miRNAs for Angus, Charolais, and KC steers, respectively. Most of the DE miRNAs were breed specific, with bta-miR-449a and bta-miR-AB-2 being differentially expressed in all three breeds. The predicted target genes of the identified DE miRNA are mainly involved in cell cycle, cell death and survival, cell signaling, cellular growth and proliferation, protein trafficking, cell morphology, cell-to-cell signaling and interaction, cellular development, molecular transport, post-translational modification, as well as nutrient metabolism (lipids, carbohydrates, protein and amino acid). Our results provide insights into the bovine hepatic miRNAome and their potential roles in molecular regulation of RFI in beef cattle.

The Role of MicroRNAs in Muscle Tissue Development in Beef Cattle

In this review, we highlight information on microRNA (miRNA) identification and functional characterization in the beef for muscle and carcass composition traits, with an emphasis on Qinchuan beef cattle, and discuss the current challenges and future directions for the use of miRNA as a biomarker in cattle for breeding programs to improve meat quality and carcass traits. MicroRNAs are endogenous and non-coding RNA that have the function of making post-transcriptional modifications during the process of preadipocyte differentiation in mammals. Many studies claim that diverse miRNAs have an impact on adipogenesis. Furthermore, their target genes are associated with every phase of adipocyte differentiation. It has been confirmed that, during adipogenesis, several miRNAs are differentially expressed, including miR-204, miR-224, and miR-33. The development of mammalian skeletal muscle is sequentially controlled by somite commitment into progenitor cells, followed by their fusion and migration, the proliferation of myoblasts, and final modification into fast- and slow-twitch muscle fibers. It has been reported that miRNA in the bovine MEG3-DIO3 locus has a regulatory function for myoblast differentiation. Likewise, miR-224 has been associated with controlling the differentiation of bovine adipocytes by targeting lipoprotein lipase. Through the posttranscriptional downregulation of KLF6, miR-148a-3p disrupts the proliferation of bovine myoblasts and stimulates apoptosis while the miR-23a~27a~24-2 cluster represses adipogenesis. Additional to influences on muscle and fat, bta-mir-182, bta-mir-183, and bta-mir-338 represent regulators of proteolysis in muscle, which influences meat tenderness.