Altered microRNA expression in bovine skeletal muscle with age

Characterization of Undiscovered miRNA Involved in Tumor Necrosis Factor Alpha-Induced Atrophy in Mouse Skeletal Muscle Cell Line

Muscular atrophy is a complex catabolic condition that develops due to several inflammatory-related disorders, resulting in muscle loss. Tumor necrosis factor alpha (TNF-α) is believed to be one of the leading factors that drive inflammatory response and its progression. Until now, the link between inflammation and muscle wasting has been thoroughly investigated, and the non-coding RNA machinery is a potential connection between the candidates. This study aimed to identify specific miRNAs for muscular atrophy induced by TNF-α in the C2C12 murine myotube model. The difference in expression of fourteen known miRNAs and two newly identified miRNAs was recorded by next-generation sequencing between normal muscle cells and treated myotubes. After validation, we confirmed the difference in the expression of one novel murine miRNA (nov-mmu-miRNA-1) under different TNF-α-inducing conditions. Functional bioinformatic analyses of nov-mmu-miRNA-1 revealed the potential association with inflammation and muscle atrophy. Our results suggest that nov-mmu-miRNA-1 may trigger inflammation and muscle wasting by the downregulation of LIN28A/B, an anti-inflammatory factor in the let-7 family. Therefore, TNF-α is involved in muscle atrophy through the modulation of the miRNA cellular machinery. Here, we describe for the first time and propose a mechanism for the newly discovered miRNA, nov-mmu-miRNA-1, which may regulate inflammation and promote muscle atrophy.

The role of non-coding RNAs in muscle aging: regulatory mechanisms and therapeutic potential

Muscle aging is a complex physiological process that leads to the progressive decline in muscle mass and function, contributing to debilitating conditions in the elderly such as sarcopenia. In recent years, non-coding RNAs (ncRNAs) have been increasingly recognized as major regulators of muscle aging and related cellular processes. Here, we comprehensively review the emerging role of ncRNAs, including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), in the regulation of muscle aging. We also discuss how targeting these ncRNAs can be explored for the development of novel interventions to combat age-related muscle decline. The insights provided in this review offer a promising avenue for future research and therapeutic strategies aimed at improving muscle health during aging.

Bta-miR-181d and Bta-miR-196a mediated proliferation, differentiation, and apoptosis in Bovine Myogenic Cells

Skeletal muscle is an important component of livestock and poultry organisms. The proliferation and differentiation of myoblasts are highly coordinated processes, which rely on the regulation of miRNA. MiRNAs are widely present in organisms and play roles in various biological processes, including cell proliferation, differentiation, and apoptosis. MiR-181d and miR-196a, identified as tumor suppressors, have been found to be involved in cell proliferation, apoptosis, directed differentiation, and cancer cell invasion. However, their role in beef cattle skeletal muscle metabolism remains unclear. In this study, we discovered that overexpression of bta-miR-181d and bta-miR-196a in Qinchuan cattle myoblasts inhibited proliferation and apoptosis while promoting myogenic differentiation through EDU staining, flow cytometry analysis, immunofluorescence staining, and Western blotting. RNA-seq analysis of differential gene expression revealed that after overexpression of bta-miR-181d and bta-miR-196a, the differentially expressed genes were mainly enriched in the PI3K-Akt and MAPK signaling pathways. Furthermore, the phosphorylation levels of key proteins p-AKT in the PI3K signaling pathway and p-MAPK in the MAPK signaling pathway were significantly decreased after overexpression of bta-miR-181d and bta-miR-196a. Overall, this study provides preliminary evidence that bta-miR-181d and bta-miR-196a may regulate proliferation, apoptosis, and differentiation processes in Qinchuan cattle myoblasts by affecting the phosphorylation status of key proteins in PI3K-Akt and MAPK-ERK signaling pathways.

Deep Small RNA Sequencing Reveals Important miRNAs Related to Muscle Development and Intramuscular Fat Deposition in Longissimus dorsi Muscle From Different Goat Breeds

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have been shown to play important post-transcriptional regulatory roles in the growth and development of skeletal muscle tissues. However, limited research into the effect of miRNAs on muscle development in goats has been reported. In this study, Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with significant phenotype difference in meat production performance were selected and the difference in Longissimus dorsi muscle tissue expression profile of miRNAs between the two goat breeds was then compared using small RNA sequencing. A total of 1,623 miRNAs were identified in Longissimus dorsi muscle tissues of the two goat breeds, including 410 known caprine miRNAs, 928 known species-conserved miRNAs and 285 novel miRNAs. Of these, 1,142 were co-expressed in both breeds, while 230 and 251 miRNAs were only expressed in LC and ZB goats, respectively. Compared with ZB goats, 24 up-regulated miRNAs and 135 miRNAs down-regulated were screened in LC goats. A miRNA-mRNA interaction network showed that the differentially expressed miRNAs would target important functional genes associated with muscle development and intramuscular fat deposition. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the target genes of differentially expressed miRNAs were significantly enriched in Ras, Rap 1, FoxO, and Hippo signaling pathways. This study suggested that these differentially expressed miRNAs may be responsible for the phenotype differences in meat production performance between the two goat breeds, thereby providing an improved understanding of the roles of miRNAs in muscle tissue of goats.

miR-152 Regulates Bovine Myoblast Proliferation by Targeting KLF6

Though miRNAs have been reported to regulate bovine myoblast proliferation, but many miRNAs still need to be further explored. Specifically, miR-152 is a highly expressed miRNA in cattle skeletal muscle tissues, but its function in skeletal muscle development is unknown. Herein, we aimed to investigate the role of miR-152 in regulating bovine myoblast proliferation. Functionally, RT-qPCR, Western blotting, EdU assay, and flow cytometry detection results showed that miR-152 inhibited bovine myoblast proliferation. Mechanistically, we demonstrated transcription factor KLF6 was a target gene of miR-152 by means of bioinformatics software prediction and dual-luciferase report analysis, which had been demonstrated to be favorable for myoblast proliferation. Collectively, our research suggested that miR-152 inhibits bovine myoblast proliferation via targeting KLF6.

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.

Bidirectional regulation of genistein on the proliferation and differentiation of C2C12 myoblasts

Genistein is a widely studied phytoestrogen. The effects of genistein on myoblasts were reported long ago, but the conclusions are controversial. In this study, we evaluated the effects of different concentrations of genistein on C2C12 myoblasts. Genistein treatment promoted myoblast proliferation in a dose-dependent manner in the concentration range of 0-2 µM/L, reaching its maximum effect at 2 µM/L. Proliferation then declined, and a concentration higher than 20 µM/L showed significant inhibition. In addition, genistein treatment promoted myoblast differentiation at a dose of 10 µM/L. However, at treatment concentrations higher than 10 µM/L, the effect on myoblast differentiation was rapidly inhibited as the concentration increased. Genistein treatment also down-regulated the expression of miR-222, resulting in increased expression of its target genes, MyoG, MyoD, and ERα and thereby promoting myoblast differentiation. Our results suggest that genistein has a dose-dependent and bidirectional regulation effect on myoblast proliferation and differentiation. We also found that genistein is a miRNA inducer, and it specifically affects the expression of miR-222 to regulate myoblast differentiation.

Expression of corticosteroid hormone receptors, prereceptors, and molecular chaperones in hypothalamic-pituitary-adrenal axis and adipose tissue after the administration of growth promoters in veal calves

The action of glucocorticoids on target tissues is regulated by the glucocorticoid and mineralocorticoid receptors (codified by the NR3C1 and NR3C2 gene, respectively). Moreover, the prereceptor system, represented by the hydroxysteroid 11-beta dehydrogenases (HSD11Bs), catalyzes the interconversion from active glucocorticoids into inactive compounds. This study aimed to determine whether the expression of the prereceptor system, the corticosteroid receptors, and the molecules regulating their intracellular trafficking (FKBP prolyl isomerase 4 and FKBP prolyl isomerase 5) could be regulated in the hypothalamic-pituitary-adrenal axis and in different type of adipose tissue of calves by the administration of dexamethasone in combination with estradiol or prednisolone. Research about the glucocorticoid effects on bovine target tissues may allow development of new diagnostic methods that use potential molecular biomarkers of glucocorticoid treatment. The administration of dexamethasone in combination with estradiol increased the gene expression of HSD11B1 (P < 0.01), HSD11B2 (P < 0.05), NR3C1 (P < 0.01), and NR3C2 (P < 0.01) in the adrenal glands; NR3C2 in the intramuscular adipose tissue (P < 0.01), and HSD11B1 in the subcutaneous adipose tissue (P < 0.01). Prednisolone administration increased the gene expression of HSD11B1 (P < 0.01), NR3C1 (P < 0.05), and NR3C2 (P < 0.05) in the adrenal glands and HSD11B1 (P < 0.01) in the subcutaneous adipose tissue. Interestingly, most of the examined tissues/organs showed a significant variation of FKBP5 gene expression after the administration of dexamethasone in combination with estradiol. So, these changes suggest that the FKBP5 gene expression could be a possible biomarker of the illegal dexamethasone administration in calves.

Expression profiles of the MXD3 gene and association of sequence variants with growth traits in Xianan and Qinchuan cattle

Max dimerization protein 3 (MXD3) belongs to the MYC superfamily of basic helix‐loop‐helix leucine zipper transcription factors, and MXD3‐MAX heterodimers can bind to promoters of target genes to modulate their expression. The aim of this study was to determine the MXD3 mRNA expression levels in various cattle tissues comprising heart, liver, spleen, lung, kidney, Longissimus dorsi muscle and subcutaneous fat in Chinese Qinchuan and Xianan cattle breeds. The RT‐qPCR data showed that the MXD3 gene was variably expressed between all tissues and at levels that were significantly different between two breeds (p < .05). We used the polymerase chain reaction‐restriction fragment length polymorphism (PCR‐RFLP) method to investigate the possible association between single‐nucleotide polymorphisms (SNP) within the MXD3 gene and five different growth traits in cattle. We found two intronic SNPs (g.2694 C>T and g.3801 T>C) and one SNP in 3′untranslated region (3′UTR) (g.6263 G>A) of MXD3 gene. Association analysis revealed strong associations between pairwise and triple SNP combinations and the growth traits. Based on these results, we suggest that MXD3 polymorphisms could be useful as molecular markers in the Chinese beef cattle breeding program.

Altered miRNA and mRNA Expression in Sika Deer Skeletal Muscle with Age

Studies of the gene and miRNA expression profiles associated with the postnatal late growth, development, and aging of skeletal muscle are lacking in sika deer. To understand the molecular mechanisms of the growth and development of sika deer skeletal muscle, we used de novo RNA sequencing (RNA-seq) and microRNA sequencing (miRNA-seq) analyses to determine the differentially expressed (DE) unigenes and miRNAs from skeletal muscle tissues at 1, 3, 5, and 10 years in sika deer. A total of 51,716 unigenes, 171 known miRNAs, and 60 novel miRNAs were identified based on four mRNA and small RNA libraries. A total of 2,044 unigenes and 11 miRNAs were differentially expressed between adolescence and juvenile sika deer, 1,946 unigenes and 4 miRNAs were differentially expressed between adult and adolescent sika deer, and 2,209 unigenes and 1 miRNAs were differentially expressed between aged and adult sika deer. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that DE unigenes and miRNA were mainly related to energy and substance metabolism, processes that are closely associate with the growth, development, and aging of skeletal muscle. We also constructed mRNA–mRNA and miRNA–mRNA interaction networks related to the growth, development, and aging of skeletal muscle. The results show that mRNA (Myh1, Myh2, Myh7, ACTN3, etc.) and miRNAs (miR-133a, miR-133c, miR-192, miR-151-3p, etc.) may play important roles in muscle growth and development, and mRNA (WWP1, DEK, UCP3, FUS, etc.) and miRNAs (miR-17-5p, miR-378b, miR-199a-5p, miR-7, etc.) may have key roles in muscle aging. In this study, we determined the dynamic miRNA and unigenes transcriptome in muscle tissue for the first time in sika deer. The age-dependent miRNAs and unigenes identified will offer insights into the molecular mechanism underlying muscle development, growth, and maintenance and will also provide valuable information for sika deer genetic breeding.

miR-148a-3p regulates proliferation and apoptosis of bovine muscle cells by targeting KLF6

Skeletal muscle development is regulated by a series of regulatory factors, and also including noncoding RNA, especially microRNAs (miRNAs). Recently, miR-148a has been found to be involved in murine C2C12 differentiation by targeting ROCK1. However, the function of miR-148a-3p for the proliferation and apoptosis of bovine muscle cells has not been determined. In this study, we found that miR-148a-3p was highly expressed in fetal bovine skeletal muscle and exhibited a decreasing trend in muscle cells during its growth phase. Functional studies indicated that gain of miR-148a-3p inhibited the proliferation of bovine muscle cells and promoted apoptosis. Conversely, interference with miR-148a-3p inhibitor promoted muscle cell proliferation and inhibited its apoptosis. Mechanistically, KLF6 was confirmed as a new potential target gene of miR-148a-3p by TargetScan software prediction and the dual-luciferase assay verification. Additionally, after a gain or loss of KLF6, the function of KLF6 for muscle cell proliferation and apoptosis was opposite to that of miR-148a-3p. Collectively, these findings proposed a novel avenue whereby miR-148a-3p impeded bovine myoblast cell proliferation and promoted apoptosis through the posttranscriptional downregulation of KLF6.

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

Background

Feed efficiency is one of the most important parameters that affect beef production costs. The energy metabolism of skeletal muscle greatly contributes to variations in feed efficiency. However, information regarding differences in proteins involved in the energy metabolism of the skeletal muscle in beef cattle divergently identified for feed efficiency is scarce. In this study, we aimed to investigate energy metabolism of skeletal muscle of Nellore beef cattle, identified for low and high residual feed intake using a proteomics approach. We further assessed the expression of candidate microRNAs as a one of the possible mechanisms controlling the biosynthesis of the proteins involved in energy metabolism that were differentially abundant between high and low residual feed intake animals.

Results

A greater abundance of 14–3-3 protein epsilon (P = 0.01) was observed in skeletal muscle of residual feed intake (RFI) high animals (RFI-High). Conversely, a greater abundance of Heat Shock Protein Beta 1 (P < 0.01) was observed in the skeletal muscle of RFI-Low cattle. A greater mRNA expression of YWHAE, which encodes the 14–3-3 protein epsilon, was also observed in the skeletal muscle of RFI-High animals (P = 0.01). A lower mRNA expression of HSPB1, which encodes the Heat Shock Protein Beta 1, was observed in the skeletal muscle of RFI-High animals (P = 0.01). The miR-665 was identified as a potential regulator of the 14–3-3 protein epsilon, and its expression was greater in RFI-Low animals (P < .001). A greater expression of miR-34a (P = 0.01) and miR-2899 (P < .001) was observed in the skeletal muscle of RFI-High animals, as both miRNAs were identified as potential regulators of HSPB1 expression.

Conclusion

Our results show that Nellore cattle divergently identified for feed efficiency by RFI present changes in the abundance of proteins involved in energy expenditure in skeletal muscle. Moreover, our data point towards that miR-665, miR34a and miR-2899 are likely involved in controlling both 14-3-3 epsilon and HSPB1 proteins identified as differentially abundant in the skeletal muscle of RFI-High and RFI-Low Nellore cattle.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5890-z) contains supplementary material, which is available to authorized users.

Survey of allele specific expression in bovine muscle

Allelic imbalance is a common phenomenon in mammals that plays an important role in gene regulation. An Allele Specific Expression (ASE) approach can be used to detect variants with a cis-regulatory effect on gene expression. In cattle, this type of study has only been done once in Holstein. In our study we performed a genome-wide analysis of ASE in 19 Limousine muscle samples. We identified 5,658 ASE SNPs (Single Nucleotide Polymorphisms showing allele specific expression) in 13% of genes with detectable expression in the Longissimus thoraci muscle. Interestingly we found allelic imbalance in AOX1, PALLD and CAST genes. We also found 2,107 ASE SNPs located within genomic regions associated with meat or carcass traits. In order to identify causative cis-regulatory variants explaining ASE we searched for SNPs altering binding sites of transcription factors or microRNAs. We identified one SNP in the 3’UTR region of PRNP that could be a causal regulatory variant modifying binding sites of several miRNAs. We showed that ASE is frequent within our muscle samples. Our data could be used to elucidate the molecular mechanisms underlying gene expression imbalance.

Association of Circulating Transfer RNA fragments with antibody response to Mycoplasma bovis in beef cattle

Background

High throughput sequencing allows identification of small non-coding RNAs. Transfer RNA Fragments are a class of small non-coding RNAs, and have been identified as being involved in inhibition of gene expression. Given their role, it is possible they may be involved in mediating the infection-induced defense response in the host. Therefore, the objective of this study was to identify 5′ transfer RNA fragments (tRF5s) associated with a serum antibody response to M. bovis in beef cattle.

Results

The tRF5s encoding alanine, glutamic acid, glycine, lysine, proline, selenocysteine, threonine, and valine were associated (P < 0.05) with antibody response against M. bovis. tRF5s encoding alanine, glutamine, glutamic acid, glycine, histidine, lysine, proline, selenocysteine, threonine, and valine were associated (P < 0.05) with season, which could be attributed to calf growth. There were interactions (P < 0.05) between antibody response to M. bovis and season for tRF5 encoding selenocysteine (anticodon UGA), proline (anticodon CGG), and glutamine (anticodon TTG). Selenocysteine is a rarely used amino acid that is incorporated into proteins by the opal stop codon (UGA), and its function is not well understood.

Conclusions

Differential expression of tRF5s was identified between ELISA-positive and negative animals. Production of tRF5s may be associated with a host defense mechanism triggered by bacterial infection, or it may provide some advantage to a pathogen during infection of a host. Further studies are needed to establish if tRF5s could be used as a diagnostic marker of chronic exposure.

microRNA-200a-3p enhances mitochondrial elongation by targeting mitochondrial fission factor

Mitochondria play pivotal roles in the ATP production, apoptosis and generation of reactive oxygen species. Although dynamic regulation of mitochondria morphology is a critical step to maintain cellular homeostasis, the regulatory mechanisms are not yet fully elucidated. In this study, we identified miR-200a-3p as a novel regulator of mitochondrial dynamics by targeting mitochondrial fission factor (MFF). We demonstrated that the ectopic expression of miR-200a-3p enhanced mitochondrial elongation, mitochondrial ATP synthesis, mitochondrial membrane potential and oxygen consumption rate. These results indicate that miR-200a-3p positively regulates mitochondrial elongation by downregulating MFF expression.

Physiological and Molecular Mechanisms of Methionine Restriction

Methionine restriction (MR) has been studied extensively over the last 25 years for its role in altering metabolic hallmarks of disease. Animals subjected to MR, display changes in metabolic flexibility demonstrated by increases in energy expenditure, glucose tolerance, and lifespan. These changes have been well characterized in a number of model systems and significant progress has been made in understanding how hepatic fibroblast growth factor 21 links MR to several components of its metabolic phenotype. Despite these advances, a complete understanding of mechanisms engaged by dietary MR remains elusive. In this review, we offer a brief history of MR and its known mechanisms associated with stress, metabolism, and lifespan extension. We consider the role of epigenetics in the response of animals to MR and propose a novel epigenetic pathway involving the regulation of microRNAs during MR.

PlanHab (Planetary Habitat Simulation): the combined and separate effects of 21 days bed rest and hypoxic confinement on human skeletal muscle miRNA expression

The study concerns effects of 21 days of sustained bedrest and hypoxia, alone and in combination, on skeletal muscle microRNA (miRNA) expression. It is expected that astronauts undertaking long‐duration missions will be exposed not only to microgravity but also to a hypoxic environment. The molecular machinery underlying microgravity‐induced alterations in skeletal muscle structure and function is still largely unknown. One possible regulatory mechanism is altered expression of miRNAs, a group of noncoding RNAs which down‐regulate many different target genes through increased degradation or translation of their messenger RNA. Thirteen healthy men underwent three 21‐day interventions, interspersed by 4‐month washout periods: horizontal bedrest in normoxia, bedrest in hypoxia, ambulation in hypoxia. The level of hypoxia corresponded to 4000 m altitude. miRNAs from v. lateralis muscle biopsies were analyzed using a microarray covering ≈4000 human miRNAs. Sixteen mature miRNAs were up‐regulated and three down‐regulated after bedrest. The magnitudes of these changes were small and a large portion of the miRNAs affected by bedrest was also differentially expressed after washout periods. In fact, the number of differentially expressed probe sets over time was substantially larger than what could be detected after bedrest. Still, the majority of the miRNAs (let‐7, miR‐15, miR‐25, miR‐199, miR‐133) that were differentially expressed following bedrest, belong to miRNA families previously reported in the context of muscle physiology, in particular to respond to changes in mechanical loading. Since only minor changes in miRNA expression could be detected after bedrest, our data indicate miRNA to play only a minor role in the substantial change in muscle phenotype seen with unloading.

Association of MicroRNAs with Antibody Response to Mycoplasma bovis in Beef Cattle

The objective of this study was to identify microRNAs associated with a serum antibody response to Mycoplasma bovis in beef cattle. Serum from sixteen beef calves was collected at three points: in summer after calves were born, in fall at weaning, and in the following spring. All sera collected in the summer were ELISA-negative for anti-M. bovis. By the fall, eight animals were seropositive for IgG (positive group), while eight remained negative (negative group). By spring, all animals in both groups were seropositive. MicroRNAs were extracted from sera and sequenced on the Illumina HiSeq next-generation sequencer. A total of 1,374,697 sequences mapped to microRNAs in the bovine genome. Of these, 82% of the sequences corresponded to 27 microRNAs, each represented by a minimum of 10,000 sequences. There was a statistically significant interaction between ELISA response and season for bta-miR-24-3p (P = 0.0268). All sera collected at the initial summer had a similar number of copies of this microRNA (P = 0.773). In the fall, the positive group had an increased number of copies when compared to the negative group (P = 0.021), and this grew more significant by the following spring (P = 0.0001). There were 21 microRNAs associated (P< 0.05) with season. These microRNAs could be evaluated further as candidates to potentially improve productivity in cattle. The microRNAs bta-let-7b, bta-miR- 24-3p, bta-miR- 92a, and bta-miR-423-5p, were significatly associated with ELISA status (P< 0.05). These microRNAs have been recognized as playing a role in the host defense against bacteria in humans, mice, and dairy cattle. Further studies are needed to establish if these microRNAs could be used as diagnostic marker or indicator of exposure, or whether intervention strategies could be developed as an alternative to antibiotics for controlling disease due to M. bovis.

Iron Content Affects Lipogenic Gene Expression in the Muscle of Nelore Beef Cattle

Iron (Fe) is an essential mineral for metabolism and plays a central role in a range of biochemical processes. Therefore, this study aimed to identify differentially expressed (DE) genes and metabolic pathways in Longissimus dorsi (LD) muscle from cattle with divergent iron content, as well as to investigate the likely role of these DE genes in biological processes underlying beef quality parameters. Samples for RNA extraction for sequencing and iron, copper, manganese, and zinc determination were collected from LD muscles at slaughter. Eight Nelore steers, with extreme genomic estimated breeding values for iron content (Fe-GEBV), were selected from a reference population of 373 animals. From the 49 annotated DE genes (FDR<0.05) found between the two groups, 18 were up-regulated and 31 down-regulated for the animals in the low Fe-GEBV group. The functional enrichment analyses identified several biological processes, such as lipid transport and metabolism, and cell growth. Lipid metabolism was the main pathway observed in the analysis of metabolic and canonical signaling pathways for the genes identified as DE, including the genes FASN, FABP4, and THRSP, which are functional candidates for beef quality, suggesting reduced lipogenic activities with lower iron content. Our results indicate metabolic pathways that are partially influenced by iron, contributing to a better understanding of its participation in skeletal muscle physiology.

Deep sequencing shows microRNA involvement in bovine mammary gland adaptation to diets supplemented with linseed oil or safflower oil

BACKGROUND

Bovine milk fat composition is responsive to dietary manipulation providing an avenue to modify the content of fatty acids and especially some specific unsaturated fatty acid (USFA) isomers of benefit to human health. MicroRNAs (miRNAs) regulate gene expression but their specific roles in bovine mammary gland lipogenesis are unclear. The objective of this study was to determine the expression pattern of miRNAs following mammary gland adaptation to dietary supplementation with 5 % linseed or safflower oil using next generation RNA-sequencing.

METHODS

Twenty-four Canadian Holstein dairy cows (twelve per treatment) in mid lactation were fed a control diet (total mixed ration of corn:grass silages) for 28 days followed by a treatment period (control diet supplemented with 5 % linseed or safflower oil) of 28 days. Milk samples were collected weekly for fat and individual fatty acid determination. RNA from mammary gland biopsies harvested on day-14 (control period) and on days +7 and +28 (treatment period) from six randomly selected cows per treatment was subjected to small RNA sequencing.

RESULTS

Milk fat percentage decreased significantly (P < 0.001) during treatment with the two diets as compared to the control period. The individual saturated fatty acids C4:0, C6:0, C8:0, C14:0 and C16:0 decreased significantly (P < 0.05) while five USFAs (C14:1, C18:1n11t, C20:3n3, C20:5n3 and CLA:t10c12) increased remarkably (P < 0.05) in response to both treatments. Analysis of 361 million sequence reads generated 321 known bovine miRNAs and 176 novel miRNAs. The expression of fourteen and twenty-two miRNAs was affected (P < 0.05) by linseed and safflower oil treatments, respectively. Seven miRNAs including six up-regulated (bta-miR-199c, miR-199a-3p, miR-98, miR-378, miR-148b and miR-21-5p) and one down-regulated (bta-miR-200a) were found to be regulated (P < 0.05) by both treatments, and thus considered core differentially expressed (DE) miRNAs. The gene targets of core DE miRNAs have functions related to gene expression and general cellular metabolism (P < 0.05) and are enriched in four pathways of lipid metabolism (3-phosphoinositide biosynthesis, 3-phosphoinositide degradation, D-myo-inisitol-5-phosphate metabolism and the superpathway of inositol phosphate compounds).

CONCLUSION

Our results suggest that DE miRNAs in this study might be important regulators of bovine mammary lipogenesis and metabolism. The novel miRNAs identified in this study will further enrich the bovine miRNome repertoire and contribute to understanding mammary gland biology.