Salatta Bruna Maria, Muniz Maria Malane Magalhães, Fonseca Larissa Fernanda Simielli, Mota Lucio Flavio Macedo, de Souza Teixeira Caio, Frezarim Gabriela Bonfá, Serna-García Marta, Dos Santos Silva Danielly Beraldo, Pereira Angélica Simone Cravo, Baldi Fernando, de Albuquerque Lucia Galvão
Animal Science Department, School of Agricultural and Veterinary Sciences, São Paulo State University (Unesp), Via de Acesso Paulo Donato Castellane S/N, Departamento de Zootecnia, Jaboticabal, 14884-900, SP, Brazil.
Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada.
Sci Rep. 2025 Jul 18;15(1):26109. doi: 10.1038/s41598-025-11179-4.
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.
本研究旨在鉴定根据脂肪酸谱聚类的内洛尔牛肌肉组织中差异表达的长链非编码RNA(lncRNA)。采集48头年轻公牛的胸最长肌样本,用于定量脂肪酸(FA)(肉豆蔻酸、棕榈酸、硬脂酸、油酸、亚油酸、共轭亚油酸(CLA)、α-亚麻酸以及饱和脂肪酸(SFA)、单不饱和脂肪酸(MUFA)、多不饱和脂肪酸(PUFA)、ω3、ω6、PUFA/SFA比值和ω6/ω3组),并生成用于转录组分析的RNA测序数据。采用K均值分析根据FA模式将48只动物分为三个聚类。聚类1(C1)的PUFA、ω3、ω6、亚油酸和α-亚麻酸含量显著(p≤0.05)更高。聚类3(C3)中SFA、肉豆蔻酸、棕榈酸和硬脂酸的比例显著(p≤0.05)更高,而聚类2(C2)呈现中间特征。对三种不同的比较(C1与C2、C1与C3以及C2与C3)进行差异分析,分别发现22、28和22个差异lncRNA(折叠变化> | 2 |,p值<0.01且错误发现率(FDR)<0.05)。对于这三种比较,新的差异转录本lncRNA_15786.3、lncRNA_13894.1和lincRNA_17393.3分别与CCN1、BNIP3和CNOT2基因相互作用,似乎有助于形成富含PUFA的脂肪酸谱。这些基因负责调节生脂基因、脂质代谢、免疫反应和脂质合成。同时,基因间差异lncRNA(lincRNA_18394.1、lincRNA_2526.3和lincRNA_17681.1)与DDX1、EIF4E和APOL3基因相关,似乎有助于形成富含SFA的脂肪酸谱。基因DDX1通过与RNA剪接相关的基因本体论(GO)术语(GO:0008380)富集,而其他基因(如EIF4E和APOL3)则通过与脂质转运(GO:0006869)、定位(GO:0010876)以及细胞对脂质的反应(GO:0071396)相关的GO术语富集。这些发现为牛肉脂肪酸组成基因调控的生物学机制提供了新的见解,并可能为进一步研究lncRNA与mRNA之间的相互作用及其对肉质的潜在影响提供有价值的基础。