基于HS-SPME-GC-MS和OPLS-DA模型探究不同嫩度驴肉的关键挥发性物质成分差异

旨在检测驴肉挥发性物质成分,探究与广灵驴嫩度相关的关键差异风味物质并解析关键挥发性物质与嫩度基因之间的功能与联系。本研究选用30头生长环境和饲养条件相同、年龄相近的雌性广灵驴为研究对象,进行剪切力和肌内脂肪的测定并依据含量差异将其分为高嫩度组(HT,n=4)和低嫩度组(LT,n=4)。通过HS-SPME-GC-MS技术检测广灵驴背最长肌挥发性物质成分,利用香气活性值(odor activity value, OAV)筛选驴肉关键风味物质,并结合多元统计分析方法获得变量重要性投影(variable importance in the projection, VIP)筛选与嫩度相关的关键差异风味物质,后基于Pearson系数与转录组学进行联合分析,寻找出驴肉嫩度关键风味物质与差异基因之间的联系。结果,在广灵驴背最长肌中共鉴定出41种挥发性物质,包括醇类、醛类、烃类、酮类、酯类以及2种其他物质。通过OAV筛选出13种关键风味物质,发现影响驴肉的主要挥发性物质和贡献者是醛类。通过VIP和OAV值筛选出1-辛烯-3-醇、1-辛醇以及月桂醛既是嫩度的差异物质又是对驴肉风味有贡献的关键风味物质。利用Pearson系数筛选与关键风味物质相关的差异基因并对其进行KEGG功能分析,发现1-辛烯-3-醇的相关基因主要富集在果糖和甘露糖代谢、糖酵解/糖异生、胰高血糖素信号通路等;1-辛醇相关基因主要富集在胰岛素信号通路、丁酸代谢、Ⅱ型糖尿病、多种脂肪细胞因子信号通路等;月桂醛相关基因则多富集在次生代谢物的生物合成、胰岛素信号通路、糖酵解/糖异生、戊糖磷酸途径等。这些通路可能参与了驴肉关键风味物质1-辛烯-3-醇、1-辛醇以及月桂醛的形成。本研究利用SPME-GC-MS技术检测了驴肉的风味挥发性物质并分析了驴肉不同嫩度风味差异,筛选出1-辛烯-3-醇、1-辛醇和月桂醛3种嫩度差异物质和关键风味物质。KEGG富集结果分析发现,酵解/糖异生、胰高血糖素信号通路、PPAR信号通路、果糖和甘露糖代谢、MAPK信号通路等可能参与了1-辛烯-3-醇、1-辛醇和月桂醛3的形成,这些都为广灵驴肉质嫩度和风味的分子改良与育种提供新的方向。

基于转录组学和代谢组学研究调控驴背最长肌嫩度的分子机制

旨在基于转录组与代谢组联合分析,探究驴肉嫩度的分子调控机制。本研究以30头生长环境和饲养条件相同、33~36月龄的雌性广灵驴为研究对象,进行剪切力和肌内脂肪含量的测定。依据剪切力和肌内脂肪含量选择出8头驴并将其分为高嫩度组(HT,n=4)与低嫩度组(LT,n=4),通过转录组和代谢组分析筛选差异表达基因和差异代谢物,之后联合KEGG富集分析,构建相关网络互作图。转录组结果表明,在HT组和LT组中共发现有1253个差异表达基因,其中832个基因上调,421个基因下调。KEGG分析表明,差异表达基因主要富集在碳水化合物代谢、脂质代谢、内分泌系统、信号转导以及细胞过程等多种途径;代谢组结果表明,在HT组和LT组中鉴定出225个差异代谢物,其中上调的有154个,下调的有71个。KEGG通路分析表明,差异代谢物主要富集到碳水化合物代谢、脂质代谢、氨基酸代谢、核苷酸代谢以及信号转导等多个途径;联合分析表明,差异表达基因与差异代谢物显著富集在甘油磷脂代谢,磷酸戊糖途径,丙氨酸、天冬氨酸和谷氨酸代谢,精氨酸和脯氨酸代谢以及胰高血糖素信号通路中。显著差异表达基因和代谢物在上述代谢通路中起到了关键调控作用,可作为驴肉嫩度的潜在候选基因及代谢物,为今后探究驴肉嫩度的分子调控机制和驴肉分子育种提供理论基础。

驴肌内脂肪沉积关键调控因子研究

旨在基于转录组(RNA-Seq)和代谢组(UPLC-MS/MS)关联分析,筛选驴肌内脂肪沉积的关键调控因子。本研究选用饲养条件相同、平均体重为236.10 kg的雌性广灵驴30头,对其背最长肌进行肌内脂肪(IMF)含量的测定,选择年龄一致的驴,并根据其IMF含量的高低分为两组:低肌内脂肪组(L组,每组3头)与高肌内脂肪组(H组,每组3头),进行转录组学和代谢组学测序分析,对差异表达基因(DEGs)和差异代谢物(DAMs)进行KEGG关联分析和相关性分析,并构建相关性网络图。结果表明,在L组和H组之间共鉴定出72种差异代谢物,其中27种代谢物上调,45种代谢物下调;关联分析表明,在差异表达基因和差异代谢物之间共有35条共富集通路,其中甘油脂代谢、甘油磷脂代谢、花生四烯酸代谢、亚油酸代谢和不饱和脂肪酸的生物合成5条通路与IMF沉积相关,有8个差异表达基因和15种差异代谢物在这5条通路中富集,包括DGKA、DGAT2、PLA2G3、GPCPD1和SCD等差异表达基因以及甘油-3-磷酸(glycerol-3-phosphate,G3P)、溶血卵磷脂酸16:0(lysophosphatidic acid,LPA)、溶血卵磷脂酰胆碱16:0(lysophosphatidylcholine,LPC)、二十二碳六烯酸(docosahexaenoic acid,DHA)和花生四烯酸(arachidonic acid,AA)等差异代谢物。另外,通过网络图分析发现,甘油脂代谢通路和甘油磷脂代谢通路中的DGAT2、PLA2G3、GPCPD1和DGKA基因与它们相应的代谢产物具有高度相关性。本研究所鉴定差异表达基因可作为IMF沉积相关的潜在候选基因,为今后进一步探究驴IMF沉积的分子调控机制和驴的肉品质分子育种提供理论依据。

广灵驴DGAT 2基因克隆、生物信息学分析及组织表达研究

本研究旨在对广灵驴二酰基甘油酰基转移酶2(DGAT2)基因进行克隆,生物信息学分析和检测其在不同组织中的表达情况,为探究DGAT2基因在广灵驴脂肪沉积和提高乳脂率等方面的作用机制提供理论参考。根据GenBank上公布的马(登录号:XM_023645689.1)、双峰驼(登录号:XM_010973154.1)、绵羊(登录号:XM_027979550.1)等物种的DGAT 2基因mRNA序列,利用Primer Premier 3.0在线工具设计同源引物,应用RT-PCR法扩增DGAT 2基因序列,用生物信息学方法分析DGAT 2基因编码序列,用实时荧光定量PCR技术检测DGAT 2基因在广灵驴心脏、肝脏、脾脏、肺脏、肾脏、背最长肌、肌间脂肪、皮下脂肪组织中的表达。结果显示,广灵驴DGAT 2基因CDS序列1086 bp,编码361个氨基酸,提交到GenBank,获得登录号:MT993643,其编码序列与马、牛、双峰驼、猪、绵羊、人、小鼠的同源性分别为99.0%、92.0%、93.5%、92.0%、92.7%、85.3%、84.1%。系统进化树分析表明,驴与马的亲缘关系最近,和小鼠的关系最远。DGAT2蛋白分子质量40.96 ku,脂肪系数92.85,等电点9.16,是一种具有跨膜区的稳定碱性疏水蛋白。DGAT2蛋白有28个磷酸化修饰位点,2个糖基化修饰位点,没有信号肽,主要定位在内质网,α-螺旋(39.89%)和无规则卷曲(36.01%)是主要的二级结构。DGAT 2基因在检测的8个组织中都有表达,其中皮下脂肪中的表达量显著高于其余组织(P<0.05),其次是心脏、肝脏和肾脏,背最长肌中的表达量最低。本试验结果为探究DGAT 2基因在广灵驴脂肪沉积和提高乳品质性状的作用奠定基础。

广灵驴CAST基因的克隆与表达

对广灵驴钙蛋白酶抑制(CAST)基因进行克隆及序列分析,并用实时荧光定量PCR技术对其在心、肝、脾、肺、肾和背最长肌等6种组织中的表达水平进行分析,探究CAST基因对驴肉品质的影响.结果表明,广灵驴CAST基因的CDS编码区全长2427 bp,共编码808个氨基酸,将测序序列提交至NCBI,得到的GenBank登录号为MN231002.广灵驴CAST蛋白的氨基酸序列与马、牛、人、猪、山羊、绵羊、小鼠、鸡的同源性分别为98.6%、77.1%、70.8%、77.5%、70.1%、60.9%、65.1%和41.5%;CAST基因在广灵驴6种组织中均存在表达,其中,表达量最高的组织为背最长肌.本试验成功获得了广灵驴的CAST基因,并发现其在背最长肌中高度表达,为进一步研究CAST基因在肌肉嫩度方面的表达调控功能提供参考.

广灵驴PPARγ基因克隆与组织表达分析

以广灵驴作为试验材料,使用RT-PCR方法对广灵驴PPARγ基因CDS编码区进行克隆然后对序列进行分析,并利用实时荧光定量PCR技术检测PPARγ基因在广灵驴的心、肝、脾、肺、肾、背最长肌和皮下脂肪7种组织中的表达特征.结果表明,广灵驴PPARγ基因含有1个1425 bp的开放阅读框,可编码474个氨基酸残基,其核苷酸序列与马的同源性最高;PPARγ蛋白是一种酸性不稳定的亲水性蛋白,二级结构由α-螺旋、延伸链、β-转角和无规则卷曲4种形式组成,蛋白序列中具有3个结构域和48个磷酸化位点,无信号肽剪切位点以及跨膜区域,主要定位在细胞质和细胞核中;PPARγ基因在广灵驴的7种不同组织中均有表达但存在一定差异,在皮下脂肪中的表达水平最为丰富.本试验成功克隆了广灵驴PPARγ基因,并发现其在皮下脂肪中有较高的表达水平,说明PPARγ基因可能与广灵驴脂肪沉积有关.

广灵驴HSL基因克隆、序列分析与差异表达

试验旨在对广灵驴的激素敏感脂酶(hormone sensitive lipase,HSL)基因进行克隆和序列分析,并对HSL基因在广灵驴不同组织中的差异表达水平进行分析。使用RT-PCR法扩增并克隆广灵驴HSL基因CDS区部分序列,将序列拼接后得到HSL基因完整的CDS区全长序列,并对序列进行一系列生物信息学分析,通过实时荧光定量PCR检测HSL基因mRNA在广灵驴的心脏、肝脏、脾脏、肺脏、肾脏、背最长肌和皮下脂肪7个组织中的表达情况。结果显示,广灵驴HSL基因完整的CDS区全长为2286 bp,共编码761个氨基酸,序列已提交到NCBI,登录号:MN231003。广灵驴HSL基因的核苷酸序列与马、羊驼、骆驼、猪、牛、山羊、小鼠、绵羊相应序列的同源性分别为99.6%、88.9%、88.6%、88.1%、86.9%、85.6%、80.8%、79.1%。系统进化树预测表明,广灵驴HSL基因与马的亲缘关系最近,与小鼠的亲缘关系最远。生物信息学分析发现,HSL蛋白的理论等电点为6.51,不稳定指数为56.83,亲水性的总平均值为-0.048,说明HSL是酸性不稳定的水溶性蛋白。蛋白保守域中存在N-末端结构域、α/β水解酶折叠结构域以及调节结构域。蛋白序列中共存在88个磷酸化修饰位点、25个糖基化修饰位点。蛋白中存在较强的疏水性区域,没有信号肽及跨膜区域。二级结构显示此蛋白是由α-螺旋、延伸链、β-转角和无规则卷曲构成的,分别占45.33%、11.70%、5.65%、37.32%。实时荧光定量PCR检测结果显示,HSL基因mRNA在广灵驴的7种组织中均有表达但存在差异,在皮下脂肪中表达量最高,在心脏中表达量最低,说明广灵驴HSL基因可能在体内脂肪沉积中发挥着非常重要的作用。该试验为进一步研究HSL蛋白功能及其在广灵驴脂肪沉积中代谢调控机制提供了理论基础。

基于WGCNA技术研究驴肉嫩度基因及代谢物调控网络

旨在通过WGCNA技术解析转录组和代谢组数据,探究驴肉嫩度调控机制。试验动物采用24~36月龄的健康雌性广灵驴(平均体重236.10 kg),将驴肉剪切力和肌内脂肪含量作为表型数据,以每个样品3个重复进行表型数据测定。本试验基于前期研究的具有显著剪切力和肌内脂肪含量差异的14个广灵驴背最长肌样本的转录组和代谢组测序数据,运用WGCNA技术筛选与驴肉嫩度相关的基因及代谢物并进行转录组与代谢组联合分析,解析嫩度相关基因与代谢物。结果表明,利用WGCNA技术通过|r|≥0.5以及P≤0.05筛选标准得到3个与嫩度相关的关键基因模块Greenyellow、Darkgrey、Darkgreen以及2个关键代谢物模块Brown、Yellow。对关键基因模块进行GO富集分析发现,模块内基因主要在甘油磷脂的生物合成、脂质氧化、脂肪酸β-氧化、细胞大分子分解代谢过程、肌肉器官发育、钙离子结合等GO功能上富集。存在于关键模块内的基因及代谢物经KEGG功能富集分析发现,其大多集中在精氨酸和脯氨酸代谢、Wnt信号通路、蛋白质消化吸收、脂肪酸代谢、TCA循环、胰高血糖素信号通路、甘油磷脂代谢、嘌呤代谢、β-丙氨酸代谢等通路上。联合分析表明,丙氨酸、天冬氨酸和谷氨酸代谢、精氨酸和脯氨酸代谢、β-丙氨酸代谢以及PPAR信号通路可能调控驴肉嫩度。WGCNA及联合KEGG共富集分析筛选到的丙氨酸、天冬氨酸和谷氨酸代谢、精氨酸和脯氨酸代谢、β-丙氨酸代谢以及PPAR信号通路可能对调控驴肉嫩度有重要作用;而GAD1、PPAT、NIT2、AGMAT、CARNS1、ACOXL以及腺苷酸基琥珀酸、L-脯氨酸、L-谷氨酸、肌酸、高肌肽、肌肽、泛酸、(9S)-羟基十八碳二烯酸则可能是影响驴肉嫩度的候选基因及代谢物。本试验可为今后广灵驴肉质嫩度的分子调控与改良育种提供一定的理论基础。

广灵驴FTO基因克隆、序列分析及组织差异表达

【目的】对广灵驴脂肪和肥胖相关基因(fat mass and obesity-associated gene,FTO)进行克隆及序列分析,同时检测其在广灵驴各组织中的表达差异,以探究广灵驴FTO基因结构对其生理代谢功能的影响。【方法】运用RT-PCR技术扩增并克隆广灵驴FTO基因CDS序列,进行基因及蛋白质功能分析,同时运用实时荧光定量PCR方法检测FTO基因在广灵驴7种组织(心脏、肝脏、脾脏、肺脏、肾脏、背最长肌及皮下脂肪)中的差异表达。【结果】广灵驴FTO基因CDS区序列长1518 bp,编码505个氨基酸,序列提交至NCBI,登录号:MZ169553。广灵驴FTO基因与马、猪、牛、人、羊驼、绵羊和山羊的相似性为99.3%、90.3%、89.5%、90.8%、90.7%、89.2%和89.2%;系统进化树分析发现,广灵驴与马的亲缘关系最近。FTO蛋白分子质量为58.35 ku,理论等电点为5.07,脂肪系数为80.36,不稳定系数为48.82,平均疏水指数为-0.550,为不稳定的酸性亲水蛋白。FTO蛋白无信号肽和跨膜区,主要定位于细胞质中,具有34个磷酸化位点与5个糖基化位点。FTO蛋白二级结构主要以α-螺旋(43.96%)和无规则卷曲(37.82%)为主。实时荧光定量PCR分析发现,FTO基因在广灵驴7个组织中均有表达,其中在肺脏和皮下脂肪中的表达量极显著高于其他组织(P<0.01),在背最长肌中表达量最低。【结论】本试验结果可为下一步基因表达与调控脂肪沉积机制及改善驴肉品质的研究奠定基础。

广灵驴ADSL基因克隆、序列分析及组织表达研究

本研究旨在对广灵驴腺苷琥珀酸裂解酶(adenylosuccinatelyase,ADSL)基因进行克隆及生物信息学分析,并检测其在不同组织中的表达情况,为探究ADSL基因在广灵驴肌苷酸合成及风味形成中的作用机制提供理论参考。根据GenBank中公布的马(登录号:XM_001917207.5)、牛(登录号:NM_001102377.2)、猪(登录号:GU249574.1)等物种的ADSL基因mRNA序列,通过Primer Premier 3.0在线工具设计同源引物,RT-PCR法扩增并克隆ADSL基因序列,对编码序列进行结构与功能分析,最后用实时荧光定量PCR检测ADSL基因在广灵驴心脏、肝脏、脾脏、肺脏、肾脏、背最长肌组织中的表达水平。结果显示,广灵驴ADSL基因CDS长1473 bp,编码490个氨基酸,提交至NCBI,登录号:MW037837,其核苷酸序列与马、牛、双峰驼、猪、绵羊、人、小鼠的相似性分别为99.5%、90.8%、92.3%、90.4%、90.7%、90.5%和86.0%。进化树分析结果表明,广灵驴与马的种属关系最近,与小鼠的亲缘关系最远。ADSL蛋白分子质量为55.44 ku,等电点为6.52,平均疏水指数为-0.243,是一种不稳定的酸性亲水蛋白。ADSL蛋白有41个磷酸化修饰位点,6个糖基化修饰位点,没有信号肽和跨膜结构,有1个卷曲螺旋。ADSL蛋白主要定位在细胞质,α-螺旋(68.98%)是主要的二级结构。ADSL基因在广灵驴6个组织中都有表达,其中肺脏中的表达量最高,显著高于其他组织(P<0.05),其次是心脏和肝脏,脾脏、肾脏和背最长肌中的表达量最低。本研究结果为今后探究ADSL基因在广灵驴肌苷酸合成及风味形成的分子机制奠定基础。

广灵驴SCD基因克隆、生物信息学及组织差异表达分析

[目的]本试验旨在对广灵驴SCD基因CDS区进行生物信息学分析,并检测其在广灵驴组织中的表达,为探究SCD基因在广灵驴脂肪沉积中的作用提供理论参考。[方法]设计SCD基因同源引物,RT-PCR法克隆CDS序列,对SCD基因编码产物进行生物信息学预测,同时使用qRT-PCR法检测广灵驴7个组织中SCD基因的表达量。[结果]广灵驴SCD基因CDS长1080 bp,可编码359个氨基酸,提交到GenBank,获得登录号:MW132164。广灵驴编码序列与马、牛、双峰驼、猪、绵羊、人、小鼠、大鼠的同源性分别为99.4%、89.2%、90.7%、90.1%、89.5%、88.2%、81.2%、82.4%。SCD蛋白分子质量为41.55 ku,等电点为9.29,不稳定指数为45.54,平均疏水指数为-0.146,属于不稳定的碱性亲水蛋白。SCD蛋白既没有信号肽也没有卷曲螺旋,有4个跨膜区,主要定位在内质网。二级、三级结构显示,SCD蛋白主要由α-螺旋(40.67%)和无规则卷曲(40.11%)构成。SCD基因在广灵驴7个组织内都表达,在肝脏、肺脏和皮下脂肪中的表达量最高,心脏和背最长肌最低。[结论]本研究结果可为今后探究SCD基因在广灵驴脂肪沉积中的功能,以及提高广灵驴乳肉品质奠定基础。