《Food Chemistry: Molecular Sciences》：Integrated analysis of phenotypes and whole-transcriptome data for the preliminary exploration of molecular regulatory networks influencing beef quality
编辑推荐：
固原牛是中国西北地区特有的本土遗传资源，具有高粗饲料利用率和产肉潜力。然而，其肉质性状的分子机制仍不清楚。在这项探索性研究中，研究人员表征了固原牛和和牛的肉质表型和转录组谱，这两个品种在肌内脂肪（IMF）含量上存在显著差异。表型比较（n = 6）揭示了屠宰性能
固原牛是中国西北地区特有的本土遗传资源，具有高粗饲料利用率和产肉潜力。然而，其肉质性状的分子机制仍不清楚。在这项探索性研究中，研究人员表征了固原牛和和牛的肉质表型和转录组谱，这两个品种在肌内脂肪（IMF）含量上存在显著差异。表型比较（n = 6）揭示了屠宰性能、IMF含量、滴水损失、水分含量、剪切力、脂肪酸组成、氨基酸谱和肌纤维类型方面的显著种间差异。通过整合全转录组测序（n = 3）和表型分析，并基于差异表达分析结合表型相关性筛选，研究人员初步鉴定出37个与脂肪沉积（16个）、剪切力（11个）和肌肉发育（10个）相关的差异基因。基于种子区域互补性、表型相关性和信号通路富集分析，研究人员预测了潜在靶基因，并构建了16个共表达网络。需要注意的是，本研究受限于小样本量，且多组学数据本质上是关联性的，无法进行因果推断。因此，上述基因和网络具有预测性，尚未经过独立的功能验证。这些发现为固原牛关键肉质性状的潜在分子机制提供了初步的探索性见解，并为未来的功能基因验证、本土黄牛遗传改良及相关工业应用提供了参考。
论文解读
研究背景与意义
固原牛是中国西北特有的黄牛遗传资源，具备卓越的攀山能力、抗逆性与繁殖性能，但其肉质，特别是决定风味、多汁性和嫩度的关键性状（如肌内脂肪（Intramuscular Fat, IMF）含量和剪切力（Shear Force, SF）），相较于经过长期定向选育的专业肉牛品种（如和牛）存在显著缺陷。肉质是一个复杂的量化性状，受高度网络化的转录后调控。现有的竞争性内源RNA（competing endogenous RNA, ceRNA）网络研究往往过于简化，且大多局限于单一肌肉组织采样，忽略了脂肪与肌肉组织间的信号串扰及跨组织通讯对调控网络的影响。因此，从多组织视角系统剖析ceRNA网络在牛肉品质调控中的作用具有重要的理论价值。该研究旨在通过极端表型比较设计，初步描绘与肉质性状紧密相关的ceRNA调控网络景观，为固原牛及其他本土黄牛品种的肉质形成机制研究与遗传改良提供理论依据。
主要关键技术方法
研究人员选取了6头低IMF含量的固原牛（Guyuan cattle, GY）和6头高IMF含量的和牛（Wagyu cattle, WA）作为实验模型。研究采集了背最长肌（Longissimus Dorsi Muscle, LDM）和皮下脂肪组织（Subcutaneous Adipose Tissue, SAT）。通过对肉质性状（包括IMF、剪切力、氨基酸和脂肪酸组成等）进行生化测定，并结合苏木精-伊红（Hematoxylin and Eosin, HE）染色、三磷酸腺苷酶（ATPase）染色及油红O（Oil Red O）染色进行组织形态学观察。随后，研究人员对上述两种组织的全转录组（mRNA、长链非编码RNA（long non-coding RNA, lncRNA）、环状RNA（circular RNA, circRNA）和微小RNA（microRNA, miRNA））进行了测序。数据分析整合了差异表达分析（Differential Expression Analysis）与加权基因共表达网络分析（Weighted Gene Co-expression Network Analysis, WGCNA），并通过实时荧光定量PCR（quantitative Real-Time PCR, qRT-PCR）对测序结果进行了验证，最终利用生物信息学工具预测了潜在的ceRNA调控网络。
研究结果
3.1. 屠宰性能与肉质性状分析
表型评估显示，和牛组的宰前活重、胴体重和屠宰率均显著高于固原牛组。在肉质方面，和牛表现出更高的IMF含量、更低的含水量和剪切力，以及更高的脂肪沉积标志物基因（如PPARγ、FABP4、LPL等）表达水平。脂肪酸分析表明，和牛的单不饱和脂肪酸（Monounsaturated Fatty Acid, MUFA）/饱和脂肪酸（Saturated Fatty Acid, SFA）比率更高，而固原牛的多不饱和脂肪酸（Polyunsaturated Fatty Acid, PUFA）/SFA比率更高。氨基酸分析则显示，两组在总必需氨基酸（Essential Amino Acid, EAA）、非必需氨基酸（Non-Essential Amino Acid, NEAA）、鲜味氨基酸（Delicious Amino Acids, DAA）和甜味氨基酸（Sweet Amino Acids, SAA）含量上存在显著差异。
3.2. 肌纤维特征分析
组织学分析表明，固原牛的肌纤维直径和横截面积大于和牛。ATPase染色显示，固原牛I型肌纤维的密度、数量比例和面积比小于和牛，而II型肌纤维的直径、面积和数量比例更大。基因表达分析进一步证实，固原牛肌肉组织中MyHC IIa、MyHC IIb和MyHC IIx的表达显著上调，而MyHC I的表达显著下调，表明固原牛肌肉代谢更倾向于糖酵解型，这与其较高的剪切力（较差的嫩度）直接相关。
3.3. 肉质性状相关mRNA的筛选
通过差异表达分析与WGCNA的交叉筛选，研究人员分别鉴定出与皮下脂肪组织（SAT）、IMF、剪切力（SF）和肌肉发育相关的候选mRNA。其中，STC1、ADAMTS15和TNNI2被预测参与脂肪细胞分化和脂质积累；CRYAB与肉质嫩度和大理石花纹特征相关；MYLPF、TNNI2和MYH2则在肌肉发育中起关键作用。STRING数据库互作分析进一步锁定了VCAM1、ITGAM等10个连接度最高的核心候选基因。
3.4. miRNA-mRNA关联的筛选
研究人员筛选出与IMF相关的bta-miR-486，与SAT相关的bta-miR-206、bta-miR-196a和bta-miR-6120-3p，以及与SF相关的bta-miR-486、bta-miR-22-5p和bta-miR-504。这些miRNA的表达趋势经qRT-PCR验证与测序结果一致，并进一步构建了miRNA-mRNA调控网络示意图。
3.5. lncRNA-miRNA-mRNA调控网络的构建
在lncRNA层面，研究人员筛选出23个与脂肪生成相关的候选lncRNA和12个与剪切力相关的候选lncRNA。选取模块成员资格（Module Membership, MM）> 0.9且基因显著性（Gene Significance, GS）> 0.9的核心lncRNA（如TCONS_00130221、TCONS_0053529等）构建了lncRNA介导的竞争性内源RNA调控网络。
3.6. circRNA-miRNA-mRNA调控网络的构建
对于circRNA，研究人员预测chr12_50842014_50873294+和chr25_28389214_28452131+可能是与IMF相关的候选分子；chr12_50842014_50873294+、chr25_28389214_28452131+、chr7_83426962_83432176+和chr8_38717329_38741975+可能与SAT相关；chr8_38717329_38741975_+可能与SF相关。Sanger测序验证了circRNA的环状结构，qRT-PCR验证了其在不同组织中的表达差异，并据此绘制了circRNA-miRNA-mRNA网络图。
讨论与结论
讨论部分指出，和牛优异的生长性能和高IMF含量伴随着特定的脂质标记基因高表达及偏向氧化代谢的肌纤维类型，这与其卓越的嫩度和风味物质基础相符。通过整合差异表达基因（Differentially Expressed Genes, DEGs）与WGCNA，本研究成功提高了候选基因筛选的准确性，并初步预测了涉及MAPK信号通路、Wnt信号通路等核心生物学过程的ceRNA调控框架。尽管小样本量和极端表型选择可能限制了统计效力和普适性，且相关性分析无法替代因果关系论证，但这些发现为后续实验提供了可检验的假设。
结论部分重申，通过整合全转录组测序与表型数据，本研究筛选出了与SAT沉积、IMF沉积、剪切力和肌肉发育相关的一系列潜在mRNA和ncRNA，并预测了多个调控网络。这些发现是探索性的，完全基于计算预测，旨在为未来的功能验证奠定基础。
