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 F(2)-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.