用高通量测序技术研究松辽黑猪与长白猪背最长肌 mRNA和lncRNA的差异表达

【背景】肌肉生长和脂肪沉积是猪生产中非常重要的经济性状,同时也是复杂的数量性状,受到众多因素的影响。长白猪原产于丹麦, 是目前世界上使用最为广泛的瘦肉型猪种之一, 具有生长速度快, 瘦肉率高等特点,但肉质欠佳。松辽黑猪是以民猪、长白猪和杜洛克猪为素材经过三元杂交育成的黑色母系地方培育品种,具有繁殖率高、肉质优良、适应性强、耐粗饲等特性。二者在肉质方面存在较大差异。【目的】采用高通量测序技术对松辽黑猪和长白猪的背最长肌转录组进行分析,筛选影响猪肌肉生长、肉质和脂肪沉积的关键mRNA和长链非编码RNA（long non-coding RNA, lncRNA）,为猪肉品质研究提供新的参考信息。【方法】采集6头松辽黑猪和6头长白猪的背最长肌作为样本,提取其RNA,将每个品种内的个体取等量RNA混池,采用Illumina HiSeq 2500高通量测序技术对混池后的样本RNA进行测序,然后对所获得的Reads进行比对、注释和差异表达等分析,筛选出差异表达mRNA和lncRNA并进行相关生物学功能富集分析。而后在筛选出的差异表达mRNA和lncRNA中各挑选10个基因,包括5个上调基因和5个下调基因,采用荧光定量PCR（Real-time PCR, RT-PCR）方法鉴定所选mRNA和lncRNA的表达水平,并验证测序结果的可靠性。【结果】将得到的原始数据进行质量控制后,每个样本分别得到103 M和150 M reads,其中98%以上reads质量较高,满足进行下一步分析的要求。将reads比对到参考基因组上,共得到26 774个转录本,其中包括8 428个新的转录本。采用edgeR软件包进行差异表达分析,参考标准是错误发现率（false discovery rate, FDR）小于0.05和差异倍数大于2（fold change, FC）。在2个猪种中共得到4 239个显著差异表达的mRNA,其中在松辽黑猪肌肉中2 023个上调,2 216个下调。其中HLCS、BTD、DGKα、LPIN1、FGF1、FGF10、FGFR1和ZNF7等基因参与脂类代谢和肌肉发育相关的调控。差异表达的mRNA的生物学功能富集分析发现,其主要富集在细胞发育、生物代谢调控、肌肉发育、脂类代谢等相关的信号通路和生物代谢途径上。同时,在2个猪种中还检测到的178个差异表达显著的lncRNAs,其中在松辽黑猪肌肉组织中有84个上调,94个下调。联合差异表达lncRNAs的靶基因预测和差异表达基因结果分析,显示有4个差异表达的lncRNA与差异表达mRNA存在潜在的调控关系,分别是TCONS-00041713、TCONS-00041712、ENSSSCT00000034982、ENSSSCT00000034269。同时10个差异表达mRNA和10个差异表达lncRNA的RT-PCR结果显示：其表达水平在两组中同样差异显著,且与转录组测序中表达的趋势相似,表明高通量测序结果具有可靠性。 【结论】 获得了影响猪肌肉生长、肉质和脂肪沉积的关键mRNA和lncRNA,同时发现了几对差异表达lncRNA和mRNA之间可能存在着潜在的调控关系。这些潜在调控关系的发现,可能对于研究调控肌肉组织代谢活动的分子机制的具有一定意义,并为其提供新的参考信息。