Skeletal muscle lncRNA profile associated with fatty acids in Nellore beef cattle

Abstract

This study aimed to identify differentially expressed (DE) long non-coding RNAs (lncRNAs) in muscle tissue of Nellore cattle clustered by their fatty acid profile. Longissimus thoracis muscle samples from 48 young bulls were used to quantify fatty acid (FA) (myristic, palmitic, stearic, oleic, linoleic, conjugated linoleic (CLA), α-linolenic and the groups of saturated fatty acids (SFA), monounsaturated (MUFA), polyunsaturated (PUFA), ω3, ω6, PUFA/SFA ratio and ω6/ω3) and to generate RNA-Sequencing data for transcriptomic analyses. The K-means analysis was used to classify the 48 animals into three clusters based on their FA patterns. The C1 had significantly (p ≤ 0.05) higher PUFA, ω3, ω6, linoleic and α-linolenic content. The proportion of SFA, myristic, palmitic and stearic were significantly (p ≤ 0.05) higher in C3, while C2 presented an intermediate profile. DE analyses were performed on three different comparisons, C1 vs. C2, C1 vs. C3 and C2 vs. C3, and 22, 28 and 22 DE lncRNAs (fold change > | 2 |, p-value < 0.01 and false discovery rate (FDR) < 0.05) were found, respectively. For three comparisons, the novel DE transcripts, lncRNA_15786.3, lncRNA_13894.1 and lincRNA_17393.3 interacted with CCN1, BNIP3, and CNOT2 genes, respectively, and appeared to contribute to a PUFA-enriched fatty acid profile. These genes are responsible for regulating the lipogenic genes, lipid metabolism, immune response and lipid synthesis. Meanwhile, the intergenic DE lncRNAs (lincRNA_18394.1, lincRNA_2526.3 and lincRNA_17681.1) were associated with the genes DDX1, EIF4E and APOL3, and appeared to contribute to a SFA-enriched fatty acid profile. The gene DDX1 was enriched by GO terms related to RNA splicing (GO:0008380), while the other genes (e.g., EIF4E and APOL3) were enriched to GO terms related to lipid transport (GO:0006869), localization (GO:0010876) and to cellular response to lipid (GO:0071396). These findings offer new insights into the biological mechanisms underlying the gene regulation of FA composition in beef and may provide a valuable foundation for further investigations regarding the interactions between lncRNAs and mRNAs, as well as their potential impact on meat quality.