
2022年4月7日,扬州大学教育部农业与农产品安全国际联合研究实验室在期刊Foods(IF:5.561)发表题为“Preliminary Results about Lamb Meat Tenderness Based on the Study of Novel Isoforms and Alternative Splicing Regulation Pathways Using Iso-seq, RNA-seq and CTCF ChIP-seq Data”的研究论文。该研究以羊肉质不同的两个杂交品种(Dorper × Hu × Hu,DHH 和 Dorper × Dorper × Hu,DDH)为研究对象。联合Iso-seq、RNA-seq和CTCF ChIP-seq等技术,探讨了由CTCF调控的肉嫩度相关的新异构体,并揭示了选择性剪接产生异构体的可能调控机制。爱基百客为该研究提供ChIP-seq的技术支持。

肉嫩度由几个内在因素(如蛋白水解活性、糖原量、纤维类型和结缔组织)和外在因素(如动物品种和营养条件)决定。在这些因素中,钙蛋白酶系统(calpain system)在肉的嫩度中起着重要的作用。此前有研究表明,CAPN3(Calpain3)蛋白和mRNA的表达水平与肉嫩度显著相关。此外,越多越多证据表明由于选择性剪接,从一个基因转录的不同异构体可能发挥不同的生物学作用。PacBio长度亚型测序(Iso-seq),三代测序技术,它能直接获得全长的异构体。联合传统的RNA-seq,Iso-seq可以促进与复杂性状相关新异构体的鉴定。尽管差异表达基因的表达研究已经在绵羊的肉质性状中研究过,但是目前对绵羊肌肉中新型异构体的研究还很少。
6只6月大公羊,3只Dorper ×Hu × Hu(DHH)和3只Dorper × Dorper × Hu (DDH)。

在本研究中,屠宰后共测量了3个肉质性状。3个肉类性状的均值和标准差见表1。这些性状中,DDH品种的剪力显著高于DHH,其他肉类性状在DDH和DHH的性状没差异。

为了更新绵羊基因组注释(Generic Feature Format, GFF),这对准确量化异构体丰度至关重要,作者汇集了6个RNA样本,实施Iso-seq测序。对基因位点、isoforms、可变剪接事件进行统计分析(图2)。

6份RNA-seq数据比对到更新的基因组注释文件,用以鉴定DDH和DHH之间的差异表达异构体。共鉴定到624个差异表达异构体。这624个差异表达异构体从492个基因转录而来,暗示一些基因可以产生超过1个差异表达异构体。所有检测到的差异表达异构体显著富集于280个GO项;显著富集于17个KEGG通路(图4)。最重要的途径是胰高血糖素信号通路(ko049222)。

图3 差异表达异构体(DEIs)火山图


为了验证差异表达异构体的结果,利用RT-PCR技术检测5个差异表达异构体(DEIs)。这5个异构体趋势与转录组数据一致。

CTCF ChIP-seq选择两个肉质不同的肌肉样本。经过质量控制和比对,有一半以上的比对reads位于基因间区域(图6),其次是内含子、启动子、外显子、3’UTR和5’UTR。在两个肉嫩度不同的样品检测到4388个差异peaks。差异峰数量最多分布在基因间的基因组区域,其次是内含子区域、启动子区域、外显子区域、5’UTR和3’UTR(图7)。作者还对这些差异峰的motif进行预测,共得到了66个motif。

图6 CTCF染色质免疫沉淀测序(ChIP-seq)在绵羊基因组中的reads分布

为了调查CTCF在羊肌肉组织中调控可变剪接的作用,对转录DEIs中的基因与位于差异peaks中的基因进行重叠分析。共找到了86个重叠基因(图8)

图8 差异peaks中的基因和转录差异表达异构体的基因的韦恩图在该研究中,DDH和DHH之间共鉴定到624个DEIs。比如,转录于CAPN3 异构体7.524.18可能与肉嫩度有关。转录DEIs与CTCF ChIP-seq得到的差异峰重叠基因有86个。这些重叠中,ANKRD23产生不同的异构体,可能受甲基化介导的CTCF调控。作为初步研究,研究结果发现了与肉嫩度相关的新型异构体,并通过生物信息学分析揭示了选择性剪接产生亚型的可能调控机制。