Next Article in Journal
Arabidopsis V-ATPase d2 Subunit Plays a Role in Plant Responses to Oxidative Stress
Previous Article in Journal
Identification of QTLs for Resistance to Fusarium Head Blight Using a Doubled Haploid Population Derived from Southeastern United States Soft Red Winter Wheat Varieties AGS 2060 and AGS 2035
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 6
10.3390/genes11060700
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differentially Expressed miRNA-Gene Targets Related to Intramuscular Fat in Musculus Longissimus Dorsi of Charolais × Holstein F2-Crossbred Bulls

by
Bilal Ahmad Mir
1,*,
Henry Reyer
2,
Katrin Komolka
1,
Siriluck Ponsuksili
2,
Christa Kühn
2 and
Steffen Maak
1
1
Institute of Muscle Biology and Growth, Leibniz Institute for Farm Animal Biology (FBN), D-18196 Dummerstorf, Germany
2
Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), D-18196 Dummerstorf, Germany
*
Author to whom correspondence should be addressed.
Genes 2020, 11(6), 700; https://0-doi-org.brum.beds.ac.uk/10.3390/genes11060700
Submission received: 28 May 2020 / Revised: 22 June 2020 / Accepted: 23 June 2020 / Published: 25 June 2020
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Intramuscular fat (IMF) is a meat quality indicator associated with taste and juiciness. IMF deposition, influenced by genetic and non-genetic factors, occurs through a transcriptionally coordinated process of adipogenesis. MicroRNAs (miRNAs) are transcriptional regulators of vital biological processes, including lipid metabolism and adipogenesis. However, in bovines, limited data on miRNA profiling and association with divergent intramuscular fat content, regulated exclusively by genetic parameters, have been reported. Here, a microarray experiment was performed to identify and characterize the miRNA expression pattern in the Musculus longissimus dorsi of F2-cross (Charolais × German Holstein) bulls with high and low IMF. A total of 38 differentially expressed miRNAs (DE miRNAs), including 33 upregulated and 5 downregulated (corrected p-value ≤ 0.05, FC ≥ ±1.2), were reported. Among DE miRNAs, the upregulated miRNAs miR-105a/b, miR-695, miR-1193, miR-1284, miR-1287-5p, miR-3128, miR-3178, miR-3910, miR-4443, miR-4445 and miR-4745, and the downregulated miRNAs miR-877-5p, miR-4487 and miR-4706 were identified as novel fat deposition regulators. DE miRNAs were further analyzed, along with previously identified differentially expressed genes (DEGs) from the same samples and predicted target genes, using multiple bioinformatic approaches, including target prediction tools and co-expression networks, as well as Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. We identified DE miRNAs and their gene targets associated with bovine intramuscular adipogenesis, and we provide a basis for further functional investigations.

1. Introduction

Intramuscular fat (IMF) content is an important meat quality trait that attributes to taste, tenderness and juiciness [1,2]. Improved beef palatability and quality grade are linked to increased IMF abundance [1]. Deposition of IMF, also known as marbling, is influenced both by genetic (gender, breed and genotype) and non-genetic factors (age, weight, nutrition, castration and stressors) [3,4,5]. Hence, genetic factors could be solely responsible for marbling when non-genetic determinants are identical. Furthermore, marbling is associated with the size and number of adipocytes as well as the fine balance between lipogenesis and the lipolysis rate in muscle [6,7,8]. Adipogenesis (the process of preadipocyte differentiation into mature adipocytes with lipid accumulation) involves hypertrophy and hyperplasia, i.e., the subsequent increase in the size and number of adipocytes [9]. Transcriptional control of adipogenesis is a highly complex and orchestrated process that is mediated by a cascade of expression events of molecules, including coding and non-coding RNA [9,10]. Thus, understanding the molecular regulators of marbling would help in delineating avenues to increase intramuscular fat deposition.
MicroRNAs (miRNAs) are small, endogenous and non-coding RNAs that bind to target mRNA, thereby modulating gene expression by translation inhibition or mRNA breakdown [11,12]. Multiple genes can be regulated by a single miRNA, while a single gene can be modulated by more than one miRNA [13]. MicroRNAs are involved in the fat deposition by regulating the expression of key genes in adipogenesis including PPARG, CEBPα, FABP7, INSIG1, RXR, RYR, FASN, SCD and IGF1R, among others. [14,15,16,17]. Intramuscular adipocyte proliferation and differentiation are influenced by miR-125-5p and miR-23a [15,18]. To date, most of the miRNA functional exploration findings are related to human and murine cell models, and reports of miRNA profiling in bovine fat accretion in general, and IMF in particular, are limited. For instance, profiling of fat tissue in cattle revealed 54 fat-specific miRNAs [19], while 18 miRNAs were concomitant with subcutaneous backfat thickness, with miR-378 showing the highest association [20]. Regarding marbling correlated miRNA profiling, only a few RNA-sequencing studies [21,22,23,24] and a microarray experiment [25] identified divergent miRNAs related to intramuscular fat content. Therefore, profiling of miRNAs in animals with different IMF content is imperative to gain further insight into their potential role in beef fat accretion. Moreover, miRNA regulation of differential IMF deposition in bulls while keeping non-genetic parameters constant has not been reported before.
In a previous study conducted by our group on Musculus longissimus dorsi (MLD) of Charolais × Holstein F2-crossbred bulls with low and high IMF content, mRNAs were profiled and differentially expressed genes (DEGs) determined [26]. Here, the same MLD samples were used to identify, through a microarray analysis, differentially expressed miRNAs (DE miRNA) and whether these target and regulate the expression of DEGs. Furthermore, DE miRNA–DEG co-expression network and functional classification was performed, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment elucidation of predicted gene targets were carried out. The bioinformatic analyses revealed miRNA profiles that may be specifically involved in IMF accretion. The variation in IMF of the bulls was exclusively due to genetic differences as non-genetic parameters were kept standardized and identical.

2. Materials and Methods

2.1. Animals

A total of 20 crossbred bulls from a F2-population (Charolais × German Holstein) were assigned to the groups low IMF (1.9 ± 0.6%; n = 10) or high IMF (7.0 ± 0.6%; n = 10). Three F1 sires had an identical number of sons in both groups. The bulls were kept and fed under the standardized condition as described elsewhere [27,28] and slaughtered at an age of 18 months. Samples of Musculus longissimus dorsi (MLD) were shock-frozen in liquid nitrogen and stored until further analyses. The experiment was approved by the Animal Care Committee of the State Mecklenburg-Western Pomerania, Germany (State Office for Agriculture, Food Safety and Fishery; LALLF M-V/TSD/7221.3-2.1-010/03), and husbandry and slaughter conformed to the German Law of Animal Protection.

2.2. Samples and RNA Preparation

Small RNAs were isolated from MLD samples of the 20 bulls using a miReasy Mini kit (Qiagen, Hilden, Germany) and subsequently enriched with an RNeasy MinElute Cleanup kit (Qiagen) according to the manufacturer’s instructions. Assessment of the amount and integrity of the isolated small RNAs was done on an Agilent 2100 Bioanalyser (Agilent, Santa Clara, CA, USA) with the Agilent small RNA kit.

2.3. MicroRNA Microarray Assay

We used the GeneChip miRNA 3.0 array (Affymetrix, Santa Clara, CA, USA) for miRNA expression profiling essentially as described by the manufacturer. This array included all microRNAs listed in release 17 of miRBase. The microarray is multispecies, covering almost 20,000 miRNAs, including 676 mature miRNAs assigned to Bos taurus. Labelling of small RNAs (200 ng/sample) was done with a FlashTag Biotin RNA labelling kit for Affymetrix GeneChip miRNA arrays (Genisphere, Hatfield, PA, USA). Hybridzation of the labelled small RNA with the array lasted 16 h. Subsequent washing and staining were performed in a Fluidics Station 450, and the arrays were then scanned on a G3000 GeneArray Scanner (Affymetrix). Robust multi-array average (RMA) background correction, log2 transformations and quantile normalization methods implemented in JMP Genomics 6 (SAS Institute, Cary, NC, USA) were applied.

2.4. Database and Identification of Differentially Expressed miRNAs

The data from the microarray were deposited to Gene Expression Omnibus (GEO) database and are accessible through a GEO Series accession number GSE147622. Here, a similar approach was taken as for the mRNA transcriptome data from the same samples [26] to identify DEGs, but with slight modifications. To identify DE miRNAs, a t-test was performed and p values were corrected for multiple testing using instructions detailed on the biostathandbook website [29]). MicroRNAs with FC ≥ 1.2 cutoff and a corrected p-value of ≤ 0.05 were considered to be statistically different, where FC ≥ 1.2 represented upregulated and FC ≥ −1.2 downregulated miRNAs. Only miRNAs with present calls in at least 7 out of 10 animals per group were selected for further analyses.

2.5. Prediction of Potential Target Genes

The miRWalk [30] is an open-source platform that provides predicted and validated miRNA binding sites of human, mouse, rat, dog and cow genes. This database is updated twice a year and provides updated information on miRNA–gene interaction. It uses random-forest-based approach software TarPmiR, searching the complete transcript sequence, including the 5′-UTR, CDS and 3′-UTR [30]. For the determination of mRNA targets of DE miRNAs, 32 previously identified DEGs from the same samples [26] were used. According to established miRNA regulation mechanisms, the upregulated miRNAs were considered to target downregulated DEGs and downregulated miRNAs to target upregulated DEGs. In addition to miRWalk, miRDB [31] and microRNA.org [32] were utilized to strengthen the prediction of targets. Hence, common target genes of functional DE miRNAs from DEGs were determined through these algorithms. However, DE miRNAs can potentially target other genes involved in adipogenesis as well.

2.6. Construction of miRNA-Target Co-Expression Interaction Network

The MicroRNA Target Filter tool and core analyses of QIAGEN’s Ingenuity Pathway Analysis (IPA®) [33] were used to find co-expression interactions and functional enrichment associations of DE miRNAs and DEGs. The IPA uses TargetScan, miRecords and TarBase as the miRNA target gene database. Given the fact that mammalian miRNAs, as well as their targets sites, are highly evolutionarily conserved, this information is exploited by the IPA® of predicted targets of microRNAs [34,35]. Core analyses and functional enrichment were performed to identify enriched networks and functions associated with fat deposition.

2.7. GO Function and KEGG Pathway Enrichment Analysis

Gene Ontology (GO) is a knowledge base used to annotate genes and classify characteristic biological aspects for transcriptome and high-throughput data. Kyoto Encyclopedia of Genes and Genomes (KEGG) is a database that is employed to conduct searches regarding biological pathways, genomes, chemical substances, diseases and drugs. The Database for Annotation, Visualization and Integrated Discovery (DAVID) is an online analysis tool that is utilized to provide a comprehensive functional understanding and biological meaning of large lists of genes [36,37]. List of predicted target genes of DE miRNAs obtained from IPA was filtered by removing those not involved in any pathway or involved in cancer. The remaining target genes in the list were imported to DAVID 6.8. The final target gene list had approximately 60% of DEGs present in it. Key biological processes (BP), cellular components (CC), molecular functions (MF) and KEGG pathways among target genes were analyzed by DAVID.

3. Results

3.1. Differentially Expressed miRNAs and Their Predicted mRNA Targets

We combined the results of the miRNA array with data on DEGs [26] obtained from the same MLD samples. The analysis was performed based on the criteria mentioned in Section 2.1. A total of 38 DE miRNAs were identified (Corr. p-value ≤ 0.05). Among these 38 DE miRNAs, 33 were upregulated (FC ≥ 1.2) and the remaining 5 were downregulated (FC ≥ −1.2) in bulls with high IMF (7%) compared to low IMF (1.9%) content (Table 1). The array was multispecies, hence the miRNAs chosen were from cows, humans and mice because they are highly conserved across mammalian species [34,38]. When predicting mRNA targets, 21 DE miRNAs, including 19 upregulated and 2 downregulated, were identified to interact with predicted DEG targets (Table 2). The predicted mRNA targets were chosen from 32 DEGs identified from a microarray analysis of the same MLD samples [26]. Out of 21 downregulated DEGs, a maximum of 16 were predicted targets of upregulated miRNAs miR-1224 and -1287, while only three targets were predicted for miR-421 (Table 2). However, 8 and 7 upregulated genes were predicted to be targeted by downregulated miRNAs let-7e and miR-877, respectively (Table 2). Other miRNAs were predicted to target 6 to 15 DEGs. The DEGs have been related to IMF deposition previously [26,39], hence DE miRNAs could potentially regulate fat development and accretion by modulating these and other potential gene targets as well as associated pathways.

3.2. DE miRNA and DEGs Co-Expression Network and Functional Classification

A miRNA–gene co-expression network was constructed to reveal the relationship of DE miRNAs with all DEGs that are basically involved in fat development [26] (Figure 1). The results revealed that, among upregulated miRNAs, miR-27b-3p, miR-31-5p, miR-100-5p, miR-149-5p, miR-421-3p, miR-504, miR-874-3p, miR-1287-5p, miR-1284, miR-4443 and miR-3128 may potentially inhibit the expression of several downregulated DEGs, including IGF1R. Contrarily, the downregulated miRNAs let-7e-5p, miR-4706, miR-877-5p and miR-4487 may modulate upregulated DEGs including THRSP and SCD (Figure 1). Apart from miR-331-3p in upregulated and miR-877-5p in downregulated, the remaining candidates were revealed as nodal miRNAs interacting with two to nine DEGs. The DEGs in the co-expression network are related to lipid metabolism and fat accretion [26]. Similarly, functional enrichment analysis revealed that DE miRNAs regulate DEGs associated with the biological functions “quantity of adipose tissue”, “differentiation of adipocytes”, “concentration of lipid” and “quantity of carbohydrate” (Figure 2).

3.3. Gene Ontology (GO) and KEGG Pathway Enrichment Analyses of Target Genes of the DE miRNAs

The predicted gene target list used for Gene Ontology (GO) and KEGG enrichment was obtained from IPA and is available in the supplementary data (Table S1). GO enrichment analysis indicated that DE miRNA gene targets belonged to 1002 biological processes (including cell proliferation regulation, fatty acid metabolic process, regulation of lipid metabolic and cholesterol process), 107 cellular components (including protein-lipid complex, microsome, receptor complex and cell fraction) and 185 molecular functions (including lipid binding, steroid hormone receptor activity, cytokine activity and transcription factor binding) (Table 3). KEGG pathway enrichment analysis revealed that the predicted gene targets of DE miRNAs belonged to 35 pathways, including fatty acid metabolism, Peroxisome proliferator-activated receptor (PPAR) signaling, adipocytokine signaling pathway, glycerolipid and retinol metabolism and type II diabetes mellitus (Table 4). The full list of GO and KEGG enrichment analyses is available in the supplementary data (Table S2).

4. Discussion

In our study, we evaluated the differences in miRNA expression in MLD with low and high IMF content, and we identified putative potential target genes from the mRNA microarray from our previous study [26]. The DE miRNAs were related to DEGs and to further gene targets predicted by IPA that may modulate pathways involved in fat development and deposition.
As bulls with varied IMF content were chosen despite the identical husbandry and feeding conditions, in addition to having the same sires, the causal factors for the differences should exclusively be genetic. Hence, differentially expressed genes and miRNAs need to be identified to understand their different regulatory capacity for fat deposition in MLD. Among the DE miRNAs identified in the microarray, some miRNAs are known to have a significant regulatory role in adipogenesis, lipid metabolism and fat tissue development. However, upregulated miRNAs, miR-105a/b, miR-695, miR-1193, miR-1284, miR-1287-5p, miR-3128, miR-3178, miR-3910, miR-4443, miR-4445 and miR-4745, and downregulated miRNAs, miR-877-5p, miR-4487 and miR-4706 are novel to the list of IMF regulators. Further functional investigation and validation of all these miRNAs, especially those interacting with multiple DEGs or involved in the fat-associated biological functions, may reveal whether they play an important role in marbling, lipid and fat biology. This study confirms the expression pattern of various miRNAs that have already been implicated in adipogenesis and fat tissue. For example, increased levels of let-7a in IMF and subcutaneous fat tissue [25,40], miR-101 in serum and adipose tissue of obese individuals [41], miR-1282 in nonalcoholic fat deposition in liver [42], miR-149 in rumen and fat tissue compared to liver, spleen, lungs, kidney and small intestine tissues [43], miR-504 in adipose-derived mesenchymal stem cells of obese individuals [44], miR-132 in the omental fat of obese individuals [45] and miR-421 in fatty hepatic tissues of high fat diet (HFD)-diet mice [46] have been reported. Moreover, miR-10b [47], miR-30c [48,49] miR-132 [50], miR-149 [51], miR-331 [52], miR-335 [53] and miR-425 [54] were all upregulated in adipogenesis and accelerated the process of adipogenic differentiation. Interestingly, miR-10b through targeting PPARA [55], miR-101 by modulating phosphatidylserine synthase 1 (PTDSS1) and fatty acid elongase 5 (ELOVL5) genes [56], miR-100 curtailing IGF1R and SREBF1 [57,58] and miR-135a regulating IRS2 play important roles in fat accretion, lipid and fatty acid metabolism and insulin signaling. While miR-135a [59] and miR-31 [60,61] are downregulated in 3T3L-1 cellular adipogenesis, they have also been associated with increased glycerol/cholesterol levels [62] and fat deposition through angiotensin [63], respectively. The expression of miR-874 has been confirmed and found increased in cattle with moderate intramuscular fat [24]. To date, concerning the downregulated DE miRNAs, no functional analyses or expression confirmation have been reported. Therefore, further investigation on these miRNAs would reveal their role in and regulation of intramuscular fat deposition.
Target genes of DE miRNAs were enriched for GO terms associated with biological processes, molecular functions and cellular components involved in adipose tissue regulation and lipid metabolism. These findings suggest that these DE miRNAs may play a significant role in the regulation and deposition of fat tissue in bovines. Many of these DE miRNAs, such as miR-27a/b, miR-10b-5p, miR-31-5p, miR-135a and miR-1224-5p, play a range of roles in lipid and fatty acid metabolism, and in fat development by regulating multiple genes [16,17,64,65]. The microRNA-27 family targets RXRA, FASN, SREBP1/2, PPARA and PPARG to regulate fat and lipid metabolism [14,66,67]. Fat accumulation is promoted by miR-31 inhibition through regulation of CEBPA and PPARG. MicroRNA-1224 curtails AMPKA1 to promote hepatic lipogenesis [65] and is also a negative regulator of TNFA [68], which may corroborate its potential role in fat development and biology [69]. In this study, predicted targets of miRNAs were also found to be involved in enriched KEGG pathways, including fatty acid metabolism, PPAR signaling, adipocytokine signaling pathway, glycerolipid and retinol metabolism and type II diabetes mellitus. PPAR signaling and downstream targets are essential for differentiation of both in vivo and in vitro differentiation of adipose tissue [70,71]. Moreover, adipocytokine signaling pathway targets have been related to insulin signaling, fat tissue and inflammation [72,73]. Together, GO terms and KEGG pathway enrichment suggest that DE miRNAs target mRNA genes, which are involved in bovine fat development, but this would need to be experimentally verified. Nevertheless, bioinformatic analyses of an informative data set imply the role which DE miRNAs may play in fat accretion and lipid metabolism. Further investigation into DE miRNA expression validation and functional exploration in the fat-deposition-related biological processes is fundamental to elucidate how they regulate IMF development.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/11/6/700/s1. Table S1: DE miRNAs gene targets predicted by IPA; Table S2: GO term and KEGG pathway analyses (complete list).

Author Contributions

Conceptualization, S.M., C.K. and S.P.; writing—original draft preparation, B.A.M. and S.M.; bioinformatics analyses, B.A.M. and H.R.; sample preparation and analyses, K.K. and S.P.; editing, S.M., H.R., S.P. and C.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study was in part supported by a grant from the Leibniz Association (grant no. SAW-2013-FBN-3). The publication of this article was funded by the Open Access Fund of the Leibniz Association and the Open Access Fund of the Leibniz Institute for Farm Animal Biology (FBN).

Acknowledgments

The authors are grateful to H. Tychsen and E. Schwitulla for expert technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hausman, G.J.; Dodson, M.V.; Ajuwon, K.; Azain, M.; Barnes, K.M.; Guan, L.L.; Jiang, Z.; Poulos, S.; Sainz, R.D.; Smith, S.; et al. BOARD-INVITED REVIEW: The biology and regulation of preadipocytes and adipocytes in meat animals1,2. J. Anim. Sci. 2009, 87, 1218–1246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Fan, Y.; Fu, G.; Fu, C.; Zan, L.; Tian, W. A missense mutant of the PPAR-γgene associated with carcass and meat quality traits in Chinese cattle breeds. Genet. Mol. Res. 2012, 11, 3781–3788. [Google Scholar] [CrossRef] [PubMed]
  3. Pethick, D.W.; Harper, G.; Oddy, V.H. Growth, development and nutritional manipulation of marbling in cattle: A review. Aust. J. Exp. Agric. 2004, 44, 705–715. [Google Scholar] [CrossRef]
  4. Bong, J.J.; Jeong, J.Y.; Rajasekar, P.; Cho, Y.M.; Kwon, E.G.; Kim, H.C.; Paek, B.H.; Baik, M. Differential expression of genes associated with lipid metabolism in longissimus dorsi of Korean bulls and steers. Meat Sci. 2012, 91, 284–293. [Google Scholar] [CrossRef]
  5. Albrecht, E.; Schering, L.; Liu, Y.; Komolka, K.; Kühn, C.; Wimmers, K.; Gotoh, T.; Maak, S. TRIENNIAL GROWTH AND DEVELOPMENT SYMPOSIUM: Factors influencing bovine intramuscular adipose tissue development and cellularity. J. Anim. Sci. 2017, 95, 2244–2254. [Google Scholar] [CrossRef]
  6. Cianzio, D.S.; Topel, D.G.; Whitehurst, G.B.; Beitz, D.C.; Self, H.L. Adipose Tissue Growth and Cellularity: Changes in Bovine Adipocyte Size and Number. J. Anim. Sci. 1985, 60, 970–976. [Google Scholar] [CrossRef]
  7. Hocquette, J.; Cassar-Malek, I.; Jurie, C.; Bauchart, D.; Picard, B.; Renand, G. Relationships between muscle growth potential, intramuscular fat content and different indicators of muscle fibre types in young Charolais bulls. Anim. Sci. J. 2012, 83, 750–758. [Google Scholar] [CrossRef]
  8. Martínez-Álvaro, M.; Agha, S.; Blasco, A.; Hernández, P. Muscle lipid metabolism in two rabbit lines divergently selected for intramuscular fat. J. Anim. Sci. 2017, 95, 2576–2584. [Google Scholar] [CrossRef]
  9. Lowe, C.E.; O’Rahilly, S.; Rochford, J.J. Adipogenesis at a glance. J. Cell. Sci. 2011, 124, 2681–2686. [Google Scholar] [CrossRef] [Green Version]
  10. De Sá, P.M.; Richard, A.J.; Hang, H.; Stephens, J.M. Transcriptional Regulation of Adipogenesis. Compr. Physiol. 2017, 272, 635–674. [Google Scholar] [CrossRef]
  11. O’Brien, J.; Hayder, H.; Zayed, Y.; Peng, C. Overview of MicroRNA Biogenesis, Mechanisms of Actions, and Circulation. Front. Endocrinol. 2018, 9, 402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Bartel, B. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  13. Bartel, B. MicroRNAs: Target Recognition and Regulatory Functions. Cell 2009, 136, 215–233. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Yang, Z.; Cappello, T.; Wang, L. Emerging role of microRNAs in lipid metabolism. Acta Pharm. Sin. B 2015, 5, 145–150. [Google Scholar] [CrossRef] [Green Version]
  15. Du, J.; Xu, Y.; Zhang, P.; Zhao, X.; Gan, M.; Li, Q.; Ma, J.; Tang, G.; Jiang, Y.; Wang, J.; et al. MicroRNA-125a-5p Affects Adipocytes Proliferation, Differentiation and Fatty Acid Composition of Porcine Intramuscular Fat. Int. J. Mol. Sci. 2018, 19, 501. [Google Scholar] [CrossRef] [Green Version]
  16. Engin, A.B. MicroRNA and Adipogenesis. Adv. Exp. Med. Biol. 2017, 960, 489–509. [Google Scholar] [CrossRef]
  17. Son, Y.H.; Ka, S.; Kim, A.Y.; Kim, J.B. Regulation of Adipocyte Differentiation via MicroRNAs. Endocrinol. Metab. 2014, 29, 122–135. [Google Scholar] [CrossRef] [Green Version]
  18. Guan, L.; Hu, X.; Liu, L.; Xing, Y.; Zhou, Z.; Liang, X.; Yang, Q.; Jin, S.; Bao, J.; Gao, H.; et al. bta-miR-23a involves in adipogenesis of progenitor cells derived from fetal bovine skeletal muscle. Sci. Rep. 2017, 7, 43716. [Google Scholar] [CrossRef]
  19. Gu, Z.; Eleswarapu, S.; Jiang, H. Identification and characterization of microRNAs from the bovine adipose tissue and mammary gland. FEBS Lett. 2007, 581, 981–988. [Google Scholar] [CrossRef] [Green Version]
  20. Jin, W.; Dodson, M.V.; Moore, S.; Basarab, J.A.; Guan, L.L. Characterization of microRNA expression in bovine adipose tissues: A potential regulatory mechanism of subcutaneous adipose tissue development. BMC Mol. Boil. 2010, 11, 29. [Google Scholar] [CrossRef] [Green Version]
  21. Li, N.; Zhang, Y.; Li, H.-P.; Han, L.; Yan, X.-M.; Li, H.-B.; Du, W.; Zhang, J.-S.; Yu, Q.-L. Differential expression of mRNA-miRNAs related to intramuscular fat content in the longissimus dorsi in Xinjiang brown cattle. PLoS ONE 2018, 13, e0206757. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Huang, J.; Wang, S.; Feng, X.; Liu, X.; Zhao, J.; Zheng, Q.; Wei, X.; Ma, Y. miRNA transcriptome comparison between muscle and adipose tissues indicates potential miRNAs associated with intramuscular fat in Chinese swamp buffalo. Genome 2019, 62, 729–738. [Google Scholar] [CrossRef] [PubMed]
  23. Oliveira, G.; Regitano, L.C.D.A.; Cesar, A.S.M.; Reecy, J.M.; Degaki, K.Y.; Poleti, M.; Felício, A.M.; Koltes, J.E.; Coutinho, L.L. Integrative analysis of microRNAs and mRNAs revealed regulation of composition and metabolism in Nelore cattle. BMC Genom. 2018, 19, 126. [Google Scholar] [CrossRef]
  24. Guo, Y.; Zhang, X.; Huang, W.; Miao, X. Identification and characterization of differentially expressed miRNAs in subcutaneous adipose between Wagyu and Holstein cattle. Sci. Rep. 2017, 7, 44026. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wang, H.; Zheng, Y.; Wang, G.; Li, H. Identification of microRNA and bioinformatics target gene analysis in beef cattle intramuscular fat and subcutaneous fat. Mol. BioSyst. 2013, 9, 2154. [Google Scholar] [CrossRef]
  26. Komolka, K.; Ponsuksili, S.; Albrecht, E.; Kühn, C.; Wimmers, K.; Maak, S. Gene expression profile of Musculus longissimus dorsi in bulls of a Charolais × Holstein F2-cross with divergent intramuscular fat content. Genomics data 2015, 7, 131–133. [Google Scholar] [CrossRef] [Green Version]
  27. Kühn, C.; Bellmann, O.; Voigt, J.; Wegner, J.; Guiard, V.; Ender, K. An experimental approach for studying the genetic and physiological background of nutrient transformation in cattle with respect to nutrient secretion and accretion type. Arch. Anim. Breed. 2002, 45, 317–330. [Google Scholar] [CrossRef]
  28. Pfuhl, R.; Bellmann, O.; Teuscher, F.; Ender, K.; Wegner, J.; Kühn, C. Beef versus dairy cattle: A comparison of feed conversion, carcass composition, and meat quality. Arch. Anim. Breed. 2007, 50, 59–70. [Google Scholar] [CrossRef]
  29. McDonald, J.H. Handbook of Biological Statistics, 3rd ed.; Sparky House Publishing: Baltimore, MD, USA, 2014; pp. 254–260. [Google Scholar]
  30. Sticht, C.; De La Torre, C.; Parveen, A.; Gretz, N. miRWalk: An online resource for prediction of microRNA binding sites. PLoS ONE 2018, 13, e0206239. [Google Scholar] [CrossRef]
  31. Chen, Y.; Wang, X. miRDB: An online database for prediction of functional microRNA targets. Nucleic Acids Res. 2019, 48, D127–D131. [Google Scholar] [CrossRef] [Green Version]
  32. Betel, D.; Wilson, M.; Gabow, A.; Marks, S.D.; Sander, C. The microRNA.org resource: Targets and expression. Nucleic Acids Res. 2008, 36, D149–D153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Krämer, A.; Green, J.; Pollard, J., Jr.; Tugendreich, S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics 2014, 30, 523–530. [Google Scholar] [CrossRef] [PubMed]
  34. Friedman, R.; Farh, K.K.-H.; Burge, C.B.; Bartel, B. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2008, 19, 92–105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Naeem, A.; Zhong, K.; Moisá, S.; Drackley, J.; Moyes, K.; Loor, J. Bioinformatics analysis of microRNA and putative target genes in bovine mammary tissue infected with Streptococcus uberis. J. Dairy Sci. 2012, 95, 6397–6408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  37. Huang, D.W.; Sherman, B.T.; Lempicki, R. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2008, 37, 1–13. [Google Scholar] [CrossRef] [Green Version]
  38. Weber, M.J. New human and mouse microRNA genes found by homology search. FEBS J. 2004, 272, 59–73. [Google Scholar] [CrossRef]
  39. Hudson, N.; Reverter, A.; Greenwood, P.; Guo, B.; Cafe, L.M.; Dalrymple, B.P. Longitudinal muscle gene expression patterns associated with differential intramuscular fat in cattle. Animal 2014, 9, 650–659. [Google Scholar] [CrossRef] [Green Version]
  40. Nardelli, C.; Iaffaldano, L.; Pilone, V.; Labruna, G.; Ferrigno, M.; Carlomagno, N.; Dodaro, C.A.; Forestieri, P.; Buono, P.; Salvatore, F.; et al. Changes in the MicroRNA Profile Observed in the Subcutaneous Adipose Tissue of Obese Patients after Laparoscopic Adjustable Gastric Banding. J. Obes. 2017, 2017, 1–6. [Google Scholar] [CrossRef]
  41. Higuchi, C.; Nakatsuka, A.; Eguchi, J.; Teshigawara, S.; Kanzaki, M.; Katayama, A.; Yamaguchi, S.; Takahashi, N.; Murakami, K.; Ogawa, D.; et al. Identification of Circulating miR-101, miR-375 and miR-802 as Biomarkers for Type 2 Diabetes. Metabolism 2015, 64, 489–497. [Google Scholar] [CrossRef] [Green Version]
  42. Soronen, J.; Yki-Jarvinen, H.; Zhou, Y.; Sädevirta, S.; Sarin, A.; Leivonen, M.; Sevastianova, K.; Perttilä, J.; Laurila, P.; Sigruener, A.; et al. Novel hepatic microRNAs upregulated in human nonalcoholic fatty liver disease. Physiol. Rep. 2016, 4, e12661. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Khan, R.; Raza, S.H.A.; Junjvlieke, Z.; Wang, X.; Wang, H.; Cheng, G.; Mei, C.; Elnour, I.E.; Zan, L. Bta-miR-149-5p inhibits proliferation and differentiation of bovine adipocytes through targeting CRTCs at both transcriptional and posttranscriptional levels. J. Cell. Physiol. 2020, 235, 5796–5810. [Google Scholar] [CrossRef] [PubMed]
  44. Meng, Y.; Eirin, A.; Zhu, X.-Y.; Tang, H.; Chanana, P.; Lerman, L.O.; Van Wijnen, A.J.; Lerman, L.O. Obesity-induced mitochondrial dysfunction in porcine adipose tissue-derived mesenchymal stem cells. J. Cell. Physiol. 2018, 233, 5926–5936. [Google Scholar] [CrossRef] [PubMed]
  45. Klöting, N.; Berthold, S.; Kovács, P.; Schön, M.P.; Fasshauer, M.; Ruschke, K.; Stumvoll, M.; Blüher, M. MicroRNA Expression in Human Omental and Subcutaneous Adipose Tissue. PLoS ONE 2009, 4, e4699. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Cheng, Y.; Mai, J.; Hou, T.; Ping, J. MicroRNA-421 induces hepatic mitochondrial dysfunction in non-alcoholic fatty liver disease mice by inhibiting sirtuin 3. Biochem. Biophys. Res. Commun. 2016, 474, 57–63. [Google Scholar] [CrossRef] [PubMed]
  47. Tsiloulis, T.; Pike, J.; Powell, D.; Rossello, F.J.; Canny, B.J.; Meex, R.C.R.; Watt, M.J. Impact of endurance exercise training on adipocyte microRNA expression in overweight men. FASEB J. 2016, 31, 161–171. [Google Scholar] [CrossRef] [Green Version]
  48. Irani, S.; Hussain, M.M. Role of microRNA-30c in lipid metabolism, adipogenesis, cardiac remodeling and cancer. Curr. Opin. Lipidol. 2015, 26, 139–146. [Google Scholar] [CrossRef]
  49. Karbiener, M.; Neuhold, C.; Opriessnig, P.; Prokesch, A.; Bogner-Strauss, J.G.; Scheideler, M. MicroRNA-30c promotes human adipocyte differentiation and co-repressesPAI-1andALK2. RNA Boil. 2011, 8, 850–860. [Google Scholar] [CrossRef] [Green Version]
  50. Shi, C.; Huang, F.; Gu, X.; Zhang, M.; Wen, J.; Wang, X.; You, L.; Cui, X.; Ji, C.; Guo, X. Adipogenic miRNA and meta-signature miRNAs involved in human adipocyte differentiation and obesity. Oncotarget 2016, 7, 40830–40845. [Google Scholar] [CrossRef] [Green Version]
  51. Yi, X.; Liu, J.; Wu, P.; Gong, Y.; Xu, X.; Li, W. The key microRNA on lipid droplet formation during adipogenesis from human mesenchymal stem cells. J. Cell. Physiol. 2019, 235, 328–338. [Google Scholar] [CrossRef]
  52. Chen, T.; Cui, J.; Ma, L.; Zeng, Y.; Chen, W. The Effect of MicroRNA-331-3p on Preadipocytes Proliferation and Differentiation and Fatty Acid Accumulation in Laiwu Pigs. BioMed Res. Int. 2019, 2019, 9287804-13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Zhu, L.; Chen, L.; Shi, C.-M.; Xu, G.-F.; Xu, L.-L.; Zhu, L.-L.; Guo, X.-R.; Ni, Y.; Cui, Y.; Ji, C. MiR-335, an Adipogenesis-Related MicroRNA, is Involved in Adipose Tissue Inflammation. Cell Biophys. 2013, 68, 283–290. [Google Scholar] [CrossRef] [PubMed]
  54. Qi, R.; Wang, J.; Wang, Q.; Qiu, X.; Yang, F.; Liu, Z.; Huang, J. MicroRNA-425 controls lipogenesis and lipolysis in adipocytes. Biochim. Biophys. Acta (BBA)-Mol. Cell Boil. Lipids 2019, 1864, 744–755. [Google Scholar] [CrossRef] [PubMed]
  55. Zheng, L.; Lv, G.-C.; Sheng, J.; Yang, Y. Effect of miRNA-10b in regulating cellular steatosis level by targeting PPAR-α expression, a novel mechanism for the pathogenesis of NAFLD. J. Gastroenterol. Hepatol. 2010, 25, 156–163. [Google Scholar] [CrossRef]
  56. Li, H.; Ma, Z.; Jia, L.; Li, Y.; Xu, C.; Wang, T.; Han, R.; Jiang, R.; Li, Z.; Sun, G.; et al. Systematic analysis of the regulatory functions of microRNAs in chicken hepatic lipid metabolism. Sci. Rep. 2016, 6, 31766. [Google Scholar] [CrossRef] [Green Version]
  57. Pek, S.L.T.; Sum, C.F.; Lin, M.X.; Cheng, A.K.S.; Wong, M.T.K.; Lim, S.C.; Tavintharan, S. Circulating and visceral adipose miR-100 is down-regulated in patients with obesity and Type 2 diabetes. Mol. Cell. Endocrinol. 2016, 427, 112–123. [Google Scholar] [CrossRef]
  58. Hu, Y.; Zhang, R.; Zhang, Y.; Li, J.; Grossmann, R.; Zhao, R. In ovo leptin administration affects hepatic lipid metabolism and microRNA expression in newly hatched broiler chickens. J. Anim. Sci. Biotechnol. 2012, 3, 16. [Google Scholar] [CrossRef] [Green Version]
  59. Chen, C.; Peng, Y.; Peng, J.; Jiang, S. miR-135a-5p inhibits 3T3-L1 adipogenesis through activation of canonical Wnt/ -catenin signaling. J. Mol. Endocrinol. 2014, 52, 311–320. [Google Scholar] [CrossRef] [Green Version]
  60. Sun, F.; Wang, J.; Pan, Q.; Yu, Y.; Zhang, Y.; Wan, Y.; Wang, J.; Li, X.; Hong, A. Characterization of function and regulation of miR-24-1 and miR-31. Biochem. Biophys. Res. Commun. 2009, 380, 660–665. [Google Scholar] [CrossRef]
  61. Tang, Y.-F.; Zhang, Y.; Li, X.-Y.; Li, C.; Tian, W.; Liu, L. Expression of miR-31, miR-125b-5p, and miR-326 in the Adipogenic Differentiation Process of Adipose-Derived Stem Cells. OMICS A J. Integr. Boil. 2009, 13, 331–336. [Google Scholar] [CrossRef]
  62. Monfared, Y.K.; Honardoost, M.; Sarookhani, M.R.; Farzam, S.A. Circulating miR-135 May Serve as a Novel Co-biomarker of HbA1c in Type 2 Diabetes. Appl. Biochem. Biotechnol. 2019, 191, 623–630. [Google Scholar] [CrossRef] [PubMed]
  63. Machal, J.; Novak, J.; Hezova, R.; Zlamal, F.; Vasku, A.; Slaby, O.; Bienertová-Vašků, J. Polymorphism in miR-31 and miR-584 binding site in the angiotensinogen gene differentially influences body fat distribution in both sexes. Genes Nutr. 2015, 10, 37. [Google Scholar] [CrossRef] [PubMed]
  64. Xie, H.; Lim, B.; Lodish, H.F. MicroRNAs Induced During Adipogenesis that Accelerate Fat Cell Development Are Downregulated in Obesity. Diabetes 2009, 58, 1050–1057. [Google Scholar] [CrossRef] [Green Version]
  65. Chen, T.; Yan, N.; Cheng, X.; Ji, X.; Yin, W. MiR-1224-5p enhances hepatic lipogenesis by targeting Adenosine monophosphate activated protein kinase? 1 in male mice. Endocrinology 2018, 159, 2008–2021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Ji, J.; Zhang, J.; Huang, G.; Qian, J.; Wang, X.; Mei, S. Over-expressed microRNA-27a and 27b influence fat accumulation and cell proliferation during rat hepatic stellate cell activation. FEBS Lett. 2009, 583, 759–766. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Jennewein, C.; von Knethen, A.; Schmid, T.; Brüne, B. MicroRNA-27b contributes to lipopolysaccharide-mediated peroxisome proliferator-activated receptor gamma (PPARgamma) mRNA destabilization. J. Biol. Chem. 2010, 285, 11846–11853. [Google Scholar] [CrossRef] [Green Version]
  68. Niu, Y.; Mo, D.; Qin, L.; Wang, C.; Li, A.; Zhao, X.; Wang, X.; Xiao, S.; Wang, Q.; Xie, Y.; et al. Lipopolysaccharide-induced miR-1224 negatively regulates tumour necrosis factor-α gene expression by modulating Sp1. Immunology 2011, 133, 8–20. [Google Scholar] [CrossRef]
  69. Cawthorn, W.P.; Sethi, J.K. TNF-alpha and adipocyte biology. FEBS Lett. 2007, 582, 117–131. [Google Scholar] [CrossRef] [Green Version]
  70. Rosen, E.D.; Sarraf, P.; E Troy, A.; Bradwin, G.; Moore, K.; Milstone, D.S.; Spiegelman, B.M.; Mortensen, R.M. PPARγ Is Required for the Differentiation of Adipose Tissue In Vivo and In Vitro. Mol. Cell 1999, 4, 611–617. [Google Scholar] [CrossRef]
  71. Ma, X.; Wang, N.; Zhao, W.; Xu, L. Deciphering the Roles of PPARγ in Adipocytes via Dynamic Change of Transcription Complex. Front. Endocrinol. 2018, 9, 473. [Google Scholar] [CrossRef] [Green Version]
  72. Ogunyemi, D.; Xu, J.; Mahesan, A.M.; Rad, S.; Kim, E.; Yano, J.; Alexander, C.; Rotter, J.I.; Chen, Y.-D.I. Differentially expressed genes in adipocytokine signaling pathway of adipose tissue in pregnancy. J. Diabetes Mellit. 2013, 3, 86–95. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Tilg, H.; Moschen, A.R. Adipocytokines: Mediators linking adipose tissue, inflammation and immunity. Nat. Rev. Immunol. 2006, 6, 772–783. [Google Scholar] [CrossRef] [PubMed]
Genes 11 00700 g001 550
Figure 1. MicroRNA-target gene co-expression network constructed using Ingenuity Pathway Analysis (IPA). Here, upregulation and downregulation of microRNA (miRNA)/gene are represented by red and green, respectively. The intensity of color represents the level of expression, e.g., intense red means higher expression, while light red is low expression of a miRNA or DEG. The gene inhibition and activation by miRNA is shown by blue and orange lines, respectively. Further, predicted targets of the miRNAs are connected with yellow lines. Transcription factors, receptors, mature miRNA etc. are represented by different shapes in a blue box on top right of the image.
Figure 1. MicroRNA-target gene co-expression network constructed using Ingenuity Pathway Analysis (IPA). Here, upregulation and downregulation of microRNA (miRNA)/gene are represented by red and green, respectively. The intensity of color represents the level of expression, e.g., intense red means higher expression, while light red is low expression of a miRNA or DEG. The gene inhibition and activation by miRNA is shown by blue and orange lines, respectively. Further, predicted targets of the miRNAs are connected with yellow lines. Transcription factors, receptors, mature miRNA etc. are represented by different shapes in a blue box on top right of the image.
Genes 11 00700 g001
Genes 11 00700 g002 550
Figure 2. The relationship between differentially expressed miRNAs and their target DEGs with associated biological function related to adipogenesis and fat development. Here, upregulation and downregulation of miRNA/gene are represented by red and green, respectively. The intensity of color represents the level of expression, e.g., intense red means higher expression, while light red is low expression of a miRNA or DEG. The gene inhibition and activation by miRNA is shown by blue and orange lines, respectively. Further, predicted targets of the miRNAs are connected with yellow lines. Transcription factors, receptors, mature miRNA etc. are represented by different shapes in a blue box on top right of the image.
Figure 2. The relationship between differentially expressed miRNAs and their target DEGs with associated biological function related to adipogenesis and fat development. Here, upregulation and downregulation of miRNA/gene are represented by red and green, respectively. The intensity of color represents the level of expression, e.g., intense red means higher expression, while light red is low expression of a miRNA or DEG. The gene inhibition and activation by miRNA is shown by blue and orange lines, respectively. Further, predicted targets of the miRNAs are connected with yellow lines. Transcription factors, receptors, mature miRNA etc. are represented by different shapes in a blue box on top right of the image.
Genes 11 00700 g002
Table
Table 1. Differentially expressed miRNAs (DE miRNAs) in bulls with high intramuscular fat (IMF) (7.0%) and low IMF (1.9%), corrected p-value ≤ 0.05 and the threshold of FC ≥ 1.2.
Table 1. Differentially expressed miRNAs (DE miRNAs) in bulls with high intramuscular fat (IMF) (7.0%) and low IMF (1.9%), corrected p-value ≤ 0.05 and the threshold of FC ≥ 1.2.
Upregulated miRNAsFCp-ValueCorr. p-Value
bta-let-7a-5p1.360.0010.017
bta-mir-10b-5p1.350.0070.02
bta-mir-100-5p1.320.0030.014
bta-mir-101-5p1.480.0060.02
bta-mir-105a/b1.390.0040.018
bta-mir-135a-5p1.310.0450.045
bta-mir-11931.330.0130.024
bta-mir-1224-5p1.260.0090.02
bta-mir-12821.380.0220.034
bta-mir-12841.310.0290.036
bta-mir-1287-5p1.200.03010.036
hsa-mir-27b1.270.0230.033
hsa-mir-30c-5p1.240.0270.036
hsa-mir-31-5p1.360.0440.045
hsa-mir-149-5p1.340.0080.02
hsa-mir-331-3p1.210.0350.04
hsa-mir-504-5p1.30.0210.032
hsa-mir-548a-3p1.270.0280.036
hsa-mir-874-3p1.210.0360.04
hsa-mir-31281.300.0050.018
hsa-mir-31781.210.0140.023
hsa-mir-39101.400.010.02
hsa-mir-44431.220.0080.019
hsa-mir-4445-5p1.470.0210.033
hsa-mir-4745-5p1.320.0430.045
mmu-mir-30c-5p1.210.0290.037
mmu-mir-31-5p1.440.0040.017
mmu-mir-132-3p1.220.0320.037
mmu-mir-149-5p1.450.0080.018
mmu-mir-335-5p1.350.0080.017
mmu-mir-421-3p1.250.0240.033
mmu-mir-425-5p1.260.0420.045
mmu-mir-6951.410.0020.012
Downregulated miRNAsFCp-ValueCorr. p-Value
bta-let-7e-5p−1.300.00140.015
bta-let-7f-5p−1.520.0020.011
hsa-mir-877-5p−1.300.01390.024
hsa-mir-4487−1.300.000040.001
hsa-mir-4706−1.350.0000060.0001
Table
Table 2. Predicted DE miRNA gene targets from differentially expressed genes (DEGs) with functions associated with lipid metabolism and fat deposition.
Table 2. Predicted DE miRNA gene targets from differentially expressed genes (DEGs) with functions associated with lipid metabolism and fat deposition.
Upregulated DE miRNACountPredictedDownregulated DEGTargets
miR-1224-5p16ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, BACH2, GNAS, RAD50, SMARCAD1, TAP2, CDC40, SESN1, ART3, PPDPFL, C2H2orf88, TBXAS1
miR-1287-5p16ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, BACH2, GNAS, RAD50, SMARCAD1, TAP2, CDC40, ART3, PPDPFL, C2H2orf88, TBXAS1, FRAT2
miR-31-5p15ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, BACH2, GNAS, RAD50, SMARCAD1, TAP2, CDC40, SESN1, ART3, PPDPFL, TBXAS1
miR-149-5p15ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, GNAS, RAD50, SMARCAD1, TAP2, CDC40, SESN1, ART3, PPDPFL, C2H2orf88, FRAT2
miR-331-3p15ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, BACH2, GNAS, RAD50, SMARCAD1, TAP2, SESN1, ART3, PPDPFL, C2H2orf88, CETN3
miR-874-3p15ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, BACH2, GNAS, SMARCAD1, CDC40, SESN1, ART3, PPDPFL, C2H2orf88, TBXAS1, CETN3
Let-7b14ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, BACH2, GNAS, RAD50, SMARCAD1, TAP2, CDC40, SESN1, ART3, TBXAS1
miR-425-5p14ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, BACH2, GNAS, RAD50, SMARCAD1, TAP2, CDC40, SESN1, PPDPFL, C2H2orf88
miR-128413ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, BACH2, RAD50, TAP2, CDC40, SESN1, ART3, CETN3, FRAT2
miR-128212ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, BACH2, GNAS, RAD50, SMARCAD1, CDC40, SESN1, PPDPFL
miR-504-5p12ATXN7L1, GPD2, IL4R, PROX1, SLC43A2, BACH2, GNAS, RAD50, SMARCAD1, TAP2, PPDPFL, TBXAS1
miR-105a9ART3, BACH2, BACH2, CDC40, GNAS, GPD2, GPD2, IL4R, TBXAS1
miR-105b9BACH2, PPDPFL, GNAS, GPD2, GPD2, IL4R, PROX1, RAD50, SLC43A2
miR-27b-3p8BACH2, CETN3, IL4R, SESN1, RAD50, SMARCAD1, TAP2, FRAT2
miR-132-3p8SLC43A2, BACH2, CDC40, FRAT2, GNAS, GNAS, SESN1, SESN1
miR-100-5p6CDC40, GPD2, SESN1, SLC43A2, SMARCAD1, CETN3
miR-10b-5p6ATXN7L1, GPD2, IL4R, RAD50, SLC43A2, TAP2
miR-11936BACH2, BACH2, GPD2, IL4R, PROX1, TBXAS1
miR-421-3p3GPD2, RAD50, SLC43A2
Downregulated DE miRNACountPredictedUpregulated DEGTargets
bta-let-7e-5p8ACLY, AGRN, CIDEC, THRSP, MRAS, SCD, INSIG1, MGST1
bta-miR-877-5p7ACLY, AGRN, CIDEC, THRSP, MRAS, SCD, INSIG1
Table
Table 3. Gene Ontology (GO) enrichment analyses of DE miRNA gene targets associated with lipid metabolism and fat accretion.
Table 3. Gene Ontology (GO) enrichment analyses of DE miRNA gene targets associated with lipid metabolism and fat accretion.
CategoryTermCount%p-ValueFold
Enrichment		FDR
GOTERM_BP_FATGO:0042127~regulation of cell proliferation15412.097417.95 × 10−222.2041268021.49 × 10−18
GOTERM_BP_FATGO:0006631~fatty acid metabolic process433.3778487.30 × 10−82.4462106491.37 × 10−4
GOTERM_BP_FATGO:0019216~regulation of lipid metabolic process352.7494117.02 × 10−113.5199833471.32 × 10−7
GOTERM_BP_FATGO:0008203~cholesterol metabolic process292.2780833.13 × 10−93.5505918985.86 × 10−6
GOTERM_CC_FATGO:0000267~cell fraction17914.061272.80 × 10−222.067152384.09 × 10−19
GOTERM_CC_FATGO:0005792~microsome473.6920661.22 × 10−82.4802612571.79 × 10−5
GOTERM_CC_FATGO:0043235~receptor complex342.6708565.61 × 10−113.6658006618.20 × 10−8
GOTERM_CC_FATGO:0032994~protein-lipid complex100.7855460.0013483.5733855191.952788
GOTERM_MF_FATGO:0005125~cytokine activity554.3205032.43 × 10−143.1298049534.04 × 10−11
GOTERM_MF_FATGO:0003707~steroid hormone receptor activity211.6496472.33 × 10−94.7556776563.87 × 10−6
GOTERM_MF_FATGO:0008289~lipid binding695.4202671.34 × 10−51.7014757830.022292
GOTERM_MF_FATGO:0008134~transcription factor binding846.5985867.99 × 10−81.8169840561.33 × 10−4
Table
Table 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DE miRNA gene targets associated with lipid and fatty acid metabolism, and fat development.
Table 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DE miRNA gene targets associated with lipid and fatty acid metabolism, and fat development.
CategoryTermCount%p-ValueFold
Enrichment		FDR
KEGG_PATHWAYhsa04060:Cytokine-cytokine receptor interaction866.7556951.72 × 10−112.0256342.13 × 10−8
KEGG_PATHWAYhsa00562:Inositol phosphate metabolism231.8067561.31 × 10−52.6284390.016234
KEGG_PATHWAYhsa03320:PPAR signaling pathway251.9638651.14 × 10−42.2359120.141802
KEGG_PATHWAYhsa00071:Fatty acid metabolism171.3354282.57 × 10−42.6227250.317883
KEGG_PATHWAYhsa04930:Type II diabetes mellitus171.3354280.0020122.2321062.465567
KEGG_PATHWAYhsa00561:Glycerolipid metabolism151.1783190.0092532.05703910.88421
KEGG_PATHWAYhsa00830:Retinol metabolism171.3354280.009321.94275910.95931
KEGG_PATHWAYhsa04920:Adipocytokine signaling pathway171.3354280.065051.56580656.55956

Share and Cite

MDPI and ACS Style

Mir, B.A.; Reyer, H.; Komolka, K.; Ponsuksili, S.; Kühn, C.; Maak, S. Differentially Expressed miRNA-Gene Targets Related to Intramuscular Fat in Musculus Longissimus Dorsi of Charolais × Holstein F2-Crossbred Bulls. Genes 2020, 11, 700. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11060700

AMA Style

Mir BA, Reyer H, Komolka K, Ponsuksili S, Kühn C, Maak S. Differentially Expressed miRNA-Gene Targets Related to Intramuscular Fat in Musculus Longissimus Dorsi of Charolais × Holstein F2-Crossbred Bulls. Genes. 2020; 11(6):700. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11060700

Chicago/Turabian Style

Mir, Bilal Ahmad, Henry Reyer, Katrin Komolka, Siriluck Ponsuksili, Christa Kühn, and Steffen Maak. 2020. "Differentially Expressed miRNA-Gene Targets Related to Intramuscular Fat in Musculus Longissimus Dorsi of Charolais × Holstein F2-Crossbred Bulls" Genes 11, no. 6: 700. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11060700

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop