中国是全球最大的猪肉生产和消费国,近年来生产高品质猪肉已成为养猪业发展目标。脂肪组织主要包括皮下脂肪、内脏脂肪和肌内脂肪,其中皮下脂肪和内脏脂肪在不同部位起着保温、支撑、储能、保护等功能,而肌内脂肪分布在肌肉组织内部,是形成猪肉大理石纹的关键,对猪肉的嫩度、风味和多汁性等感官品质具有重要影响
[1]。脂肪沉积过程受多种信号调控,例如:脂肪酸结合蛋白(FABPs)家族中的心脏型(H-FABP)和脂肪细胞型(A-FABP)可以促进心脏和脂肪中甘油三酯的沉积,进而促进肌内脂肪的积累
[2];过氧化物激活物受体γ(PPARγ)和CCAAT-增强子结合蛋白(C/EBP)是控制脂肪细胞形成的关键因子
[3];另外,环磷酸腺苷-蛋白激酶A(cAMP-PKA)信号通路通过下游的cAMP反应元件结合蛋白(cAMP response element-binding protein,CREB)和其共激活剂转录辅活化因子(CRTCs)发挥作用调控能量平衡
[4]。肌肉组织中的非编码RNA也参与调节肌内脂肪沉积过程,影响脂肪形成
[5]。因此,脂肪沉积是营养与遗传因素二者相互作用的结果,对动物生产效益和肉质具有重要影响。
CRTC2(CREB-regulated transcription coactivator 2)是CREB的转录共激活因子之一,CREB家族有
CRTC1、
CRTC2和
CRTC3成员。
CRTC1在大脑组织中高度表达,
CRTC2主要在肌肉和肝脏中表达,
CRTC3主要在肺组织和淋巴细胞中表达。除了作为CREB的共激活因子外,
CRTC2还参与能量代谢的调控,
CRTC2在空腹和饥饿信号(如胰高血糖素、肾上腺素下降等)传导下,结合相应受体,激活cAMP-PKA信号通路,导致
CRTC2被激活并进入细胞核与CREB结合,从而促进胰岛素基因的转录,调节机体的血糖平衡
[6]。此外,
CRTC2还与其他转录因子相互作用,调控脂肪酸合成、糖原合成和分解等代谢通路的基因转录
[7]。前期研究发现,与脂质代谢相关的miR-149-5p过表达显著降低了
CRTC2的mRNA表达,表明miR-149-5p可能阻止
CRTC2的表达,从而抑制Ca
2+信号通路和脂质沉积
[8]。
非同义单核苷酸多态性(Non-synonymous single nucleotide polymorphisms,nsSNPs)是发生在编码区(Coding sequence,CDS)的单个核苷酸突变,导致氨基酸序列改变,这种变异一般是由单个碱基的转换、颠换引起
[9]。因其可能对蛋白质的空间结构和功能产生影响,通常被认为是动植物表型和基因型变异的重要原因之一。目前,关于猪
CRTC2基因的nsSNPs在调控猪脂肪沉积和代谢方面的作用机制尚不明确,因而本研究通过生物信息学分析,挖掘猪
CRTC2基因CDS区的nsSNPs,筛选可能影响猪肉品质的功能性突变位点,预测其可能的作用机理,并对其相互作用基因和蛋白进行KEGG(Kyoto encyclopedia of genes and genomes)富集分析,随后通过qRT-PCR检测猪
CRTC2基因在不同组织中的相对表达水平。本研究有助于理解猪
CRTC2基因调控脂质代谢的作用机制,为挖掘脂肪沉积和代谢机制提供参考依据。
1 材料与方法
1.1 猪CRTC2基因SNPs的检索及m6A甲基化位点预测
从Ensembl数据库(
http://asia.ensembl.org/index.html)中检索位于猪
CRTC2基因外显子区、内含子区、3'UTR(Untranslated region)和5'UTR上的SNPs,并筛选出在外显子区域上的nsSNPs。使用SRAMP在线网站(
https://www.cuilab.cn/sramp)预测
CRTC2的CDS区的N6-甲基腺苷(m6A)修饰位点,为确保低错误率,以高置信度(high confidence)预测的m6A位点作为其是否具有m6A修饰的标准
[10]。
1.2 猪CRTC2基因nsSNPs功能预测
本研究利用SNAP(
http://www.rostlab.org/services/SNAP)、SIFT(
http://blocks.fhcrc.org/sift/SIFT.html)和Polyphen-2(
http://genetics.bwh.harvard.edu/pph2)方法预测位于
CRTC2基因CDS区的nsSNPs的功能。SNAP主要预测SNPs的功能效应,其预测结果可分为中性(Neutral)和有影响(Effect)。SIFT基于氨基酸的物理特性和序列同源,主要用于检测氨基酸替换是否会导致蛋白质功能的变化,其得分越低,表明引起蛋白质功能改变的可能性越大
[11]。Polyphen-2主要进行氨基酸替换改变蛋白质结构和功能可能性的预测,分数越接近1.0,改变蛋白质结构和功能的可能性越大,阈值为0.50。
1.3 猪CRTC2基因编码蛋白质的进化保守性位点预测和进化树分析
使用MEGA11对
CRTC2的核苷酸和其编码蛋白质的氨基酸组成进行分析,并对
CRTC2编码蛋白质的氨基酸序列进行同源比对。同时,利用在线工具ConSurf(
https://consurf.tau.ac.il/)分析
CRTC2编码蛋白质的氨基酸序列中各位点的进化保守水平。
1.4 猪CRTC2编码的蛋白质稳定性预测
利用在线网站I-Mutant 2.0和Mupro(
http://mupro.proteomics.ics.uci.edu/)预测
CRTC2单点突变后编码蛋白质的稳定性,评估结果中自由能变化值
<0表示蛋白质的稳定性降低,
>0则相反。
1.5 猪CRTC2编码的蛋白质二级结构预测
1.6 CRTC2及其互作基因/蛋白KEGG通路富集分析
通过GeneMANIA(
http://genemania.org/)和STRING(
https://cn.string-db.org/)软件分别对
CRTC2基因以及其互作基因
,CRTC2编码蛋白(F1SFW7)及其互作蛋白进行预测和分析,进而预测这些通路可能发挥的作用。
1.7 组织表达谱分析
使用TRizol法提取各组织总RNA,使用逆转录试剂盒、SYBR Green Premix Pro Taq HS qPCR试剂盒(湖南艾科瑞生物工程有限公司)对RNA样品进行反转录和qPCR检测,具体操作参照试剂盒说明书进行。采用两步法进行扩增,qPCR反应体系为10 μL,生物学重复设置为3。通过2-ΔΔCt法计算基因的相对表达量,并以GAPDH作为内参对照。
2 结果与分析
2.1 猪CRTC2序列组成分析及m6A甲基化位点预测
利用MEGA11对
CRTC2基因的氨基酸和碱基组成进行分析。结果发现,Ser(丝氨酸)、Pro(脯氨酸)、Leu(亮氨酸)在氨基酸组成中占比较大;碱基C在碱基组成类别中占比最高,而T(U)最低(
表1和
图1)。此外,
CRTC2编码的蛋白质分子量为73 143.85 g/mol,等电点为7.038 5。通过SRAMP网站对
CRTC2的CDS区的N6-甲基腺苷(m6A)修饰位点进行预测,发现其CDS区仅存在3个中度可信的位点(Moderate confidence),此结果表明
CRTC2 的CDS区不存在m6A修饰(
图2)。
2.2 猪CRTC2基因上的SNPs及nsSNPs分布
从Ensembl数据库中检索到猪
CRTC2基因的转录本
CRTC2-202对应的基因座上存在190个SNPs,其中17个位于3' UTR区、1个位于5' UTR区、142个位于内含子区、30个位于外显子区。其中CDS区中有10个错义突变,即nsSNPs(
图3)。
2.3 猪CRTC2基因功能性nsSNPs预测
本研究利用SNAP、SIFT和Polyphen-2方法预测
CRTC2基因上的nsSNPs的功能。这3个网站分别预测到8个、3个和1个非中性的nsSNPs;共预测到9个潜在影响
CRTC2编码蛋白功能的nsSNPs,分别为rs709245010(P148L)、rs337849183(S199I)、rs1107669551(P567S)、rs699432551(V591M)、rs709060088(S592G)、rs695272651(E596D)、rs705078240(L600F)、rs694414401(R610H)和rs318352528(G621R)(
表2)。
2.4 猪CRTC2基因编码蛋白质的进化保守性位点预测和进化树分析
2.4.1 CRTC2基因编码蛋白质的进化保守性位点预测
利用在线工具ConSurf预测转录本
CRTC2-202对应基因座上的nsSNPs位点的保守性,结果显示,rs709245010(P148L)是高度保守并暴露的功能性位点,因此推测该保守性强的位点突变可能对其蛋白质功能产生一定影响(
表3和
图4)。
2.4.2 CRTC2基因编码蛋白质的进化树分析
在NCBI中的Protein Blast数据库检索发现该序列在人(
Homo sapiens)、小鼠(
Mus musculus)、牛(
Bos taurus)、猪(
Sus scrofa)、野兔(
Oryctolagus cuniculus)、绵羊(
Ovis aries)之间高度同源,利用MEGA11对
CRTC2编码蛋白质的氨基酸序列进行比对,结果显示,在本研究中rs709245010(P148L)多态位点在除人外的物种中均表现出其保守性(
图5(a))。进化树分析显示,猪和牛、绵羊最近,但是与人、野兔和小鼠相对较远(
图5(b))。
2.5 猪CRTC2编码蛋白质稳定性预测
使用Mupro工具评估nsSNPs对蛋白质稳定性的影响,发现除rs337849183(S199I)突变位点外,rs709245010(P148L)等9个突变位点均会使
CRTC2编码的蛋白质稳定性降低(
表4)。
2.6 rs709245010(P148L)突变位点蛋白质二级结构分析
根据上述预测分析,rs709245010(P148L)位点可能是重要的nsSNPs突变位点。利用在线软件Sopma对这个突变位点的蛋白质二级结构进行预测。结果显示,
CRTC2编码蛋白质的野生型二级结构中
α螺旋占19.19%、延伸链占6.93%、
β转角占7.07%、无规则卷曲占66.81%。其中rs709245010(P148L)突变位点导致无规则卷曲比率降至65.66%,
β转角比率提高至7.22%,延伸链百分比提高至7.36%,并且导致
α螺旋百分比提高至19.77%。因此rs709245010(P148L)位点的突变造成了
CRTC2编码的蛋白质二级结构的改变(
表5)。
2.7 CRTC2互作基因和蛋白质预测及其KEGG通路富集分析
通过GeneMANIA和STRING软件分别对
CRTC2基因及其互作基因,以及
CRTC2编码蛋白及其互作蛋白进行预测和分析,结果发现,
CRTC3、
CRTC1与
CRTC2基因存在较强的互作关系(
图6(a));CREBBP与
CRTC2编码的蛋白质存在较强的互作关系(
图6(b))。
为进一步确定
CRTC2基因的功能,本研究通过对与
CRTC2互作基因及蛋白质进行KEGG通路富集分析。结果发现,
CRTC2及其互作基因主要参与胰高血糖素信号通路和脂质和动脉粥样硬化等脂质代谢和肌肉发育相关通路(
图7(a));
CRTC2及其互作蛋白质主要参与胰高血糖素信号通路和cGMP-PKG信号通路等脂质代谢和肌肉发育相关通路(
图7(b))。因此,推测
CRTC2可能对脂质代谢有重要作用。
2.8 CRTC2基因在猪肌肉和脂肪相关组织中的表达谱分析
为进一步探索
CRTC2在脂肪沉积和肌肉发育中的作用,本研究通过qRT-PCR分析其在长白猪的不同组织中的表达情况,结果显示,
CRTC2在长白猪的肝脏中表达量显著高于其他组织(
P<0.05),其次是腹脂和背脂,其表达量均显著高于背肌和腿肌(
P<0.05)(
图8)。
3 讨论
3.1 CRTC2生物信息学分析
m6A是一种转录后甲基化修饰,广泛存在于RNA转录本的腺苷碱基上。这种修饰被认为参与RNA转录本的降解、亚细胞定位和剪接等的调节。为进一步探究
CRTC2是否存在m6A修饰,本研究通过预测发现其CDS区仅存在3个中度可信的位点,因此认为其CDS区不存在m6A修饰。本研究运用软件对猪
CRTC2基因外显子区的功能性nsSNPs进行预测,预测结果显示,rs709245010(P148L)突变导致
CRTC2编码的蛋白质结构发生改变,并且蛋白质的稳定性降低。由于保守序列在生物体生存中起着重要作用,错义突变可能降低蛋白质的稳定性,从而导致疾病发生
[12]。本研究对
CRTC2编码的蛋白质的稳定性和二级结构进行预测,蛋白质的二级结构变化会影响其稳定性,进而调控蛋白的翻译过程。此外,rs709245010(P148L)突变将脯氨酸替换为亮氨酸,亮氨酸作为支链氨基酸,在脂肪细胞分化和脂肪积累中扮演关键角色。据报道,亮氨酸对脂质代谢具有调节作用,能够促进脂肪分解和抑制脂肪合成
[13]。因此,本研究推测rs709245010(P148L)位点可能是影响猪脂肪沉积和代谢的重要功能性突变位点。本研究为后续深入研究rs709245010(P148L)位点突变对
CRTC2基因功能的影响提供参考。
3.2 CRTC2可能参与脂质代谢调节过程
脂质代谢是一个复杂的过程,涉及脂肪合成和分解代谢、脂滴变化、脂肪沉积等多个关键过程。
CRTC2是CREB家族的成员,在CREB介导的转录激活中扮演着关键角色。早期研究表明,
CRTC2在调节葡萄糖代谢和脂质代谢方面具有重要功能,它通过影响cAMP信号传导途径调节葡萄糖代谢、甾醇调节元件结合蛋白1(SREBP1)和甾醇调节元件结合蛋白2(SREBP2)来影响脂质代谢
[14],并且可以通过去磷酸化的方式影响碳水化合物代谢
[15]。此外,
CRTC2对脂肪沉积和脂质代谢产生影响。研究发现,
CRTC2参与调控脂肪细胞发育和脂肪合成的关键信号通路,并与其他调控脂质代谢的关键蛋白质如AMPK和盐诱导激酶SIK2等相互作用,共同调节脂肪沉积过程。
CRTC2在细胞中受到多种信号通路的调控,包括cAMP-PKA和Ca
2+-calmodulin依赖性信号传导。当细胞内的cAMP或Ca
2+水平增加,激活PKA或Ca
2+-calmodulin依赖性激酶,进而导致
CRTC2从细胞质转移到细胞核,在细胞核内,
CRTC2与CREB结合并激活下游基因转录
[16]。
为进一步探讨其作用,对与
CRTC2基因的互作基因进行预测和分析,结果显示,
CRTC2与
CRTC3、
CRTC1存在较强的相互作用,相关研究表明
CRTC3通过Ca
2+-cAMP信号通路在调节脂质代谢和脂肪生成起着重要作用,同时,
CRTC3还通过调节与脂质代谢相关的基因,促进磷脂的形成增加,从而扩大脂肪细胞以增加脂质储存
[17];
CRTC1选择性介导肝脏中cAMP信号传导的有益抗脂肪变性作用,且
CRTC1广泛抑制促进肝脏脂质积累基因的表达
[18]。此外,结果发现
CRTC2编码的蛋白质与CREBBP蛋白相互作用,CREBBP作为共激活因子,研究发现CREB部分通过刺激转录激活因子3(ATF3)的表达来损害脂肪细胞功能,ATF3是一种与脂联素和葡萄糖转运蛋白-4(GLUT4)结合并抑制转录的阻遏蛋白。除了对ATF3的影响外,CREB还可以调节影响脂肪细胞脂质代谢的其他细胞基因的表达
[19-20]。例如,CREB通过激活
C/EBPβ引发
C/EBPα和
PPARγ基因的表达,从而启动脂肪细胞表型相关基因的转录
[21]。通过KEGG富集分析发现,与
CRTC2相互作用的基因和蛋白参与胰高血糖素信号通路和Wnt信号通路等多条脂质代谢相关通路,其中胰高血糖素信号通路通过AMPK信号通路激活减弱脂质积累
[22],Wnt信号通路是脂肪生成中最重要的调控通路,Wnt信号转导通过抑制成脂转录因子C/EBPα和PPARγ抑制脂肪细胞的分化和脂肪生成
[23]。肝脏是脂肪沉积的关键器官
[24],且
CRTC2在调节脂质代谢和控制肝脏胆固醇合成中起着至关重要的作用
[25-26]。
qRT-PCR结果表明
CRTC2在长白猪的肝脏、背脂和腹脂等脂肪代谢相关组织中表达量较高,进一步表明
CRTC2可能在脂肪沉积和代谢过程中发挥重要作用。据报道,
CRTC2在未受刺激的肝细胞中高度磷酸化,并且主要定位于细胞质中。在用胰高血糖素或诱导细胞cAMP升高的刺激物处理后,
CRTC2被去磷酸化并易位到细胞核
[27-28]。SIK2是AMPK相关激酶家族的成员,在肝脏中,SIK2作为LKB1(Liver kinase B1)的下游糖异生抑制因子,其抑制剂HG-9-91-01能够促进转录共激活因子
CRTC2的去磷酸化,从而增强肝细胞中糖异生基因的表达和葡萄糖产生
[29]。因此推测
CRTC2可能在调节脂质代谢方面发挥作用,其为解析猪
CRTC2基因在脂质代谢中的功能提供了参考依据。
4 结论
本研究通过对猪CRTC2基因进行生物学特征和功能预测分析,发现了猪CRTC2基因上的rs709245010错义突变位点(P148L),其可能影响CRTC2编码的蛋白质结构及其稳定性;CRTC2基因参与多条与脂肪代谢相关的通路,且在肝脏中的表达量显著高于其他组织。本研究为进一步挖掘猪CRTC2基因功能提供参考依据,对解析猪脂质代谢调控机制具有参考意义。
科技创新2030-重大项目(2023ZD04046)
国家“十四五”重点研发计划项目(2021YFD1301202)
河南省农业良种联合攻关项目(2022020101)