不同切丝宽度烟丝在干燥过程中挥发性香味物质的变化
Change of the Volatile Aroma Components of Cut Tobacco in Different Cutting Widths in Drying Process
-
STRA8基因(stimulated by retinoic acid gene 8)是生殖细胞有丝分裂转变为减数分裂前的特异性表达基因[1-3],它在精子发生过程中的减数分裂前期起重要作用[4-5]。STRA8基因的表达受视黄酸(retinoic acid, RA)诱导因素的影响[6]。RA对于脊椎动物的生殖是必需的[1-2, 6],它对睾丸的正常结构和功能起维护作用[4]。RA含量缺乏时,会导致精子发生停止和细精管退化[1, 3, 6],因此它对精子发生的过程和细精管的正常功能非常重要;若RA活性较低,则其无法发挥作用,使STRA8基因无法表达,导致雄性生殖细胞一直停滞在有丝分裂阶段而无法进入减数分裂前期,更无法形成生殖细胞,从而影响精子发生过程[7]。STRA8在减数分裂过程中行使转录因子功能[4],是减数分裂起始的调节蛋白,参与雌性胚胎期和雄性出生后生殖细胞的减数分裂起始过程[8-11],其表达是雄性和雌性生殖腺及正常精子发生过程的重要条件,直接关系到RA是否有效[3, 7, 10],对生殖细胞中的RA的转导起重要作用[2, 7, 12]。
版纳微型猪近交系(Banna Mini-pig Inbred Line,BMI)是基于云南省优良地方猪种——滇南小耳猪培育成功的世界上第一个猪近交系,因其具有很多其他猪种不具有的重要特征,且基因高度纯合、遗传背景清楚,生理学、解剖学和疾病发生机理等方面与人类极为相似,可以作为良好的人类异种移植的器官供体,对其深入研究具有重要的科研价值和应用前景[13-17]。本研究克隆了BMI STRA8基因,通过测序等方法确定了STRA8基因的核苷酸序列,然后分析了该基因的核酸序列特征;利用生物信息学和比较基因组学对STRA8蛋白质进行了功能分析;通过分析多种组织的mRNA表达谱确定其组织表达水平,进而明确其发挥功能的重要组织,为今后深入研究STRA8基因的功能提供了参考。
1. 材料与方法
1.1 试验动物
版纳微型猪近交系第29世代的成年公猪,屠宰取样,液氮速冻保存。
1.2 方法
1.2.1 RNA提取和cDNA合成
利用RNA提取试剂抽提各组织总RNA,用分光光度计测定纯度、浓度,选择合格的RNA反转录为第1链cDNA。
1.2.2 设计并合成PCR引物
根据在GenBank上下载的水牛、绵羊、犬、人、仓鼠和大鼠等物种的STRA8 mRNA序列,以及猪STRA8 mRNA序列(登录号:JQ965783)和猪EST序列,使用Oligo 7软件设计特异性引物F1/R1来扩增STRA8基因;设计F2/R2特异性引物进行荧光定量,并以GAPDH基因为内参[18],扩增的产物长度为183 bp。引物序列及扩增信息如表1。
表 1 STRA8引物序列及扩增片段信息Table 1. STRA8 primer sequence and amplified fragment information名称
name序列 (5′→3′)
primer sequences产物大小/bp
product sizeSTRA8-F1 TGTGGGGAAGGCTGACTGAA 1 316 STRA8-R1 ACAGCCCACTCCAAAACATC STRA8-F2 AAGTACATCTAGCCCCAGTTCCG 143 STRA8-R2 GAAGGTCCGGTTCCTGTTTGA 1.2.3 PCR体系及程序
PCR反应体系为:ddH2O 5.75 μL、2×GC Buffer Ι 12.5 μL、dNTP 4 μL、睾丸cDNA模板1.5 μL、上下游引物各0.5 μL及Ex Taq酶0.25 μL,总体系25 μL。运行程序的退火温度为60 ℃,扩增产物经电泳后测序。
1.2.4 生物信息学分析
利用DNAStar软件进行序列组装、CDS预测、蛋白序列推导,利用ProtParam tool程序(http://web.expasy.org/protparam)预测蛋白质的分子量(Mw)和等电点(pI)。利用软件及网站SOPMA程序(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=/NPSA/npsa_sopma.html)预测推导出STRA8蛋白质的二级结构,SMART程序(http://smart.embl-heidelberg.de)预测出功能域,ProtScale程序(http://web.expasy.org/cgi-bin/protscale/protscale.pl)预测出疏水性,TMHMM 2.0程序(http://www.cbs.dtu.dk/services/TMHMM)预测跨膜结构,SignalP 4.1程序(http://www.cbs.dtu.dk/services/SignalP/)预测信号肽,PSORT Ⅱ程序(http://psort.hgc.jp/form2.html)进行亚细胞定位,NetPhos 3.1 Server程序(http://www.cbs.dtu.dk/services/NetPhos/)预测磷酸化位点、Protfun 2.2软件预测蛋白质功能。通过NCBI BLASTP程序(https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastp&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome)搜索其他物种的STRA8蛋白质序列,利用DNAStar软件中的MegAlign程序和Clustal X软件与BMI STRA8蛋白序列(KU705622)进行同源比对;基于多物种氨基酸比对结果,利用MEGA 5.2构建系统进化树。
1.2.5 STRA8基因组、mRNA分析
利用Ensembl (http://asia.ensembl.org/index.html?redirect=no)及NCBI网站(https://www.ncbi.nlm.nih.gov/)查找STRA8基因的外显子、内含子的数目、长度以及在基因组中的位置。
1.2.6 多组织荧光定量表达分析
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
2. 结果与分析
2.1 STRA8基因PCR结果
以F1/R1为引物,睾丸组织cDNA为模板,扩增的STRA8产物片段长度为1 316 bp (图1)。
2.2 STRA8基因序列及氨基酸序列
测序后拼接出BMI STRA8基因的CDS区为1 101 bp,编码366个氨基酸,已提交至NCBI GenBank,mRNA登录号和氨基酸登录号分别为KU705622和AOC89040 (图2)。
2.3 STRA8蛋白质基本特性
BMI STRA8蛋白质的一级结构信息和理化性质结果见表2。
表 2 STRA8蛋白质理化特性Table 2. The physics and chemistry properties of STRA8一级结构特性
primary structure characteristics预测结果
prediction results编码的氨基酸数number of amino acids encoded 366 等电点number of amino acids encoded (pI) 4.52 分子量/ku the molecular weight (Mw) 41.06 分子式molecular formula C1793H2821N479O591S16 正电荷残基(Arg+Lys) positive charge residue 37 负电荷残基(Asp+Glu) negative charge residue 71 平均疏水性(GRAVY) average hydrophobicity −0.537 不稳定系数(II) coefficient of instability 61.78 脂肪系数(AI) coefficient of fat 79.18 2.4 STRA8蛋白质二级结构
无规则卷曲106个氨基酸,占全部氨基酸序列的28.96%;α-螺旋191个氨基酸,占全部氨基酸序列的52.19%;延伸链结构42个氨基酸,占全部氨基酸序列的11.48%;β转角27个氨基酸,占全部氨基酸序列的7.38%。
2.5 STRA8蛋白质功能域
STRA8蛋白质有1个HLH蛋白功能域(螺旋—环—螺旋蛋白质功能域),位于34~84 AA处。在氨基端有1个HLH结构域和1个卷曲螺旋区域,在羧基端有2个低复杂度区域,HLH结构域包括2个保守性较高的α螺旋区和1个袢环(loop),袢环长度不定。
2.6 STRA8蛋白质疏水性
STRA8蛋白质的N末端、C末端疏水值均小于0,表明两端均亲水,而在第213 AA处具有最大疏水值1.944,第162 AA处具有最小疏水−3.589。
2.7 STRA8蛋白质跨膜螺旋、信号肽及亚细胞定位
STRA8蛋白质无跨膜螺旋和信号肽。其主要定位在细胞核,概率为87%;其次定位在细胞质,概率为8.7%。
2.8 STRA8蛋白质磷酸化位点分析
共发现34个磷酸化位点,包括23个丝氨酸位点、8个苏氨酸位点,3个酪氨酸位点(表3)。
表 3 预测的BMI STRA8蛋白质磷酸化位点Table 3. The prediction results of BMI STRA8 protein phosphorylated sites磷酸化氨基酸phosphorylated
amino acids位置position 分值score 位置position 分值score 位置position 分值score 丝氨酸serine 4 0.525 11 0.605 35 0.997 48 0.533 91 0.649 101 0.991 111 0.503 118 0.962 121 0.700 129 0.971 172 0.547 204 0.604 205 0.638 231 0.803 241 0.567 259 0.994 263 0.995 267 0.957 276 0.903 291 0.615 300 0.988 303 0.760 306 0.973 苏氨酸threonine 42 0.641 61 0.775 81 0.683 193 0.548 223 0.897 277 0.702 301 0.681 360 0.601 酪氨酸tyrosine 109 0.933 113 0.632 186 0.936 2.9 STRA8功能分析
Protfun 2.2预测结果(表4)显示:STRA8蛋白在复制和转录过程中发挥作用的可能性最大,概率达0.362。其次,在调节功能、转录、转录调节、中间代谢中枢等过程中发挥作用的可能性也较大。可见,此蛋白可能在复制和转录等方面具有重要作用。
表 4 STRA8功能预测Table 4. Prediction of STRA8 function功能分类functional category 概率probability 氨基酸合成amino acid biosynthesis 0.072 辅因子生物合成biosynthesis of cofactors 0.059 细胞被摸cell envelope 0.036 细胞过程cellular processes 0.031 中间代谢中枢central intermediary metabolism 0.144 能量代谢energy metabolism 0.032 脂肪酸代谢fatty acid metabolism 0.016 嘌呤和嘧啶purine and pyrimidine 0.171 调节功能regulatory functions $\geqslant$ 0.268复制和转录replication and transcription 0.362 翻译translation 0.115 转运和结合transport and binding 0.017 基因本体分类gene ontology category 概率probability 信号转导signal transducer 0.087 接收器receptor 0.003 激素hormone 0.001 结构蛋白structural protein 0.001 转运transporter 0.025 离子通道ion channel 0.011 电压门控离子通道voltage-gated ion channel 0.004 阳离子通道cation channel 0.010 转录transcription $\geqslant$ 0.168转录调节transcription regulation 0.163 应激反应stress response 0.006 免疫反应immune response 0.010 生长因子growth factor 0.005 金属离子传输metal ion transport 0.009 2.10 STRA8蛋白质系统进化分析
BMI STRA8氨基酸序列与家犬(XP_005629800)、绵羊(XP_014950882)、水牛(XP_006079922)、人(XP_011514439)、仓鼠(XP_003508946)、大鼠(XP_006224926) 6个物种STRA8氨基酸序列经同源性比对后,构建的分子系统进化树(图3、4)。
2.11 STRA8基因组、mRNA分析
在Ensembl数据库中搜索STRA8基因组DNA,结果显示该基因位于猪18号染色体NC_010460.3的14861386~14880354位置,其全长18 968 bp,包含8个外显子和7个内含子(图5)。
2.12 STRA8 组织mRNA表达分析
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
图 6 BMI STRA8基因多组织表达谱注:1.心;2.肝;3.脾;4.肺;5.肾;6.胃;7.脑;8.肌肉;9.十二指肠;10.结肠;11.精囊腺;12.前列腺;13.尿道球腺;14.睾丸;15.附睾;不同字母表示差异极显著(P<0.01)。Figure 6. The multi-tissue expression profile of BMI STRA8Note: 1. heart; 2. liver; 3. spleen; 4. lung; 5. kidney; 6. stomach; 7. brain; 8. muscle; 9. duodenum; 10. colon; 11. seminal vesicle; 12. prostate; 13. urethral ball glands; 14. testis; 15. epididymis; different letters mean extremely significant difference (P<0.01).3. 讨论
近年来对于STRA8基因的结构和功能的研究越来越深入[11],目前已有人、小鼠、山羊和鸡的研究报道[19-23],研究表明:STRA8基因在雄性和雌性动物的生殖生理中起重要作用,是生殖细胞有丝分裂转变为减数分裂前的特异表达基因。本研究采用的试验动物为BMI猪,采用电子克隆方式拼接出猪的STRA8基因序列,并设计引物,然后以BMI为试验材料进行了验证。
本研究克隆了猪STRA8基因的全长编码区序列1 316 bp,包括编码区1 101 bp;编码366个氨基酸,包含117个疏水性氨基酸、89个极性氨基酸、37个强碱性氨基酸、71个强酸性氨基酸。化学修饰是蛋白质翻译后加工的重要内容,糖基化和磷酸化是蛋白质化学修饰的两种主要方式,蛋白质磷酸化可以大大提高酶蛋白的活性,本研究发现BMI STRA8蛋白磷酸化位点共有34个,为以后的功能研究提供了参考。
蛋白质的结构是其发挥功能的基础。通过生物信息学分析获悉,BMI STRA8蛋白含有1个蛋白功能域HLH,HLH结构域是1个家族类转录因子,这类转录因子对脊椎动物和无脊椎动物发育过程中细胞的增殖、分化、生长以及细胞的生存起着重要的调控作用,该类转录因子能够形成同二聚体或者异二聚体,然后与下游目的基因调控区的DNA序列相结合,从而发挥转录调控作用[24-25]。在RA存在的情况下,HLH结构域可以开启精原细胞的减数分裂[19, 22]。在STRA8基因的二级结构中α螺旋占大多数,从而保证了STRA8蛋白的稳定性。预测结果表明:STRA8蛋白质是疏水性蛋白,通过分析氨基酸序列的亲/疏水性,从而预测蛋白质的跨膜区,这对预测蛋白是否为膜蛋白有重要意义。STRA8蛋白无信号肽序列,说明它属于非分泌性蛋白,主要在细胞质中发挥作用。STRA8蛋白无跨膜螺旋结构,说明它属于非跨膜蛋白。STRA8的功能分析显示该蛋白在复制和转录过程中发挥功能的可能性最大,在调节功能、转录、转录调节、中间代谢中枢等中发挥功能的可能性都较大,显著高于其他功能。可见,此蛋白可能在复制和转录等方面具有重要作用。
BMI与其他6个物种的STRA8蛋白序列相似性均在70%以上,说明STRA8基因在物种进化过程中保守性较高。系统进化树显示水牛与绵羊聚在一起,为一小类,这两者皆为牛科动物;然后和猪聚为一类,三者皆为偶蹄目;再与家犬聚为一类,四者都为哺乳纲动物;之后与人聚为一类,此5个物种皆为脊索动物门;中国仓鼠与大鼠聚为一类,两者皆为啮齿目;两个大的分支聚为一类,两类都为脊索动物门,符合系统进化树的分类标准。
STRA8基因在雌性和雄性成年动物的性腺中都有很高的表达[2]。在成年小鼠中,STRA8基因RNA的表达似乎被限制在睾丸中,且这个表达仅限于生殖细胞,而不是体细胞[4]。米美玲等[26]研究表明:STRA8基因仅在成年小鼠的睾丸中表达;也有研究表明:STRA8基因表达仅限于减数分裂前成人睾丸的生殖细胞[1]。MIYAMOTO等[4]研究表明:人STRA8基因也仅在睾丸中表达。本研究对BMI实时荧光定量多组织表达谱分析表明:STRA8基因在性腺睾丸组织的表达量显著高于其他14种组织,在睾丸中的表达量是其他组织表达量的8~2 302倍,推断该基因主要在猪的睾丸中发挥重要功能。在生产实践中,公猪的繁殖力高低会直接影响其后代的数量和质量,精液品质是衡量公猪繁殖力高低的重要指标。而在近交系猪生产中,公猪的繁殖力直接影响近交系的继代繁育,因为近交系相互交配的公母猪要求有血缘上的高度近亲关系,与配母猪的公猪选择余地较小,一旦公猪繁殖力有问题,会直接影响某个亚系甚至家系的继代繁育。STRA8基因在BMI睾丸组织中高表达,表明该基因与精子生成有关。研究结果将为深入研究STRA8基因在猪中的功能提供重要依据。
-
图 1 不同烟丝样品挥发性香味物质前两个主成分得分在二维平面上的投影图
注:◆ 干燥前烟丝;■ 0.85 mm宽度烟丝干燥后;▲ 0.95 mm宽度烟丝干燥后;● 1.05 mm宽度烟丝干燥后
Figure 1. Projection of the first two principal components score in the two-dimensional plane of volatile aroma components of different cut tobacco samples
Note: ◆ cut tobacco before drying; ■ 0.85 mm-width cut tobacco; ▲ 0.95 mm-width cut tobacco; ● 1.05 mm-width cut tobacco
图 2 不同烟丝样品挥发性香味物质主成分综合得分的箱线图
注:A. 干燥前烟丝;B. 0.85 mm切丝宽度烟丝干燥后;C. 0.95 mm切丝宽度烟丝干燥后;D. 1.05 mm切丝宽度烟丝干燥后。
Figure 2. Boxplot of principal components comprehensive score of volatile aroma components of different cut tobacco samples
Note: A. cut tobacco before drying; B. cut tobacco after drying with cutting width at 0.85 mm; C. cut tobacco after drying with cutting width at 0.95 mm; D. cut tobacco after drying with cutting width at 1.05 mm.
表 1 干燥前及不同切丝宽度烟丝干燥后的挥发性香味物质测定结果
Table 1 Determination results of volatile aroma components of cut tobacco before drying and cut tobacco with different cutting width after drying
致香前体物及其降解产物
aroma precursor and its degradation products干燥前烟丝
cut tobacco before drying0.85 mm切丝宽度烟丝干燥后
0.85 mm-width cut tobacco after dryingGQ-1 GQ-2 GQ-3 GQ-4 GQ-5 GQ-6 GH0.85-1 GH0.85-2 GH0.85-3 GH0.85-4 GH0.85-5 GH0.85-6 美拉德反应产物
maillard reaction products糠醛
furfural1.76 1.79 1.78 1.81 1.80 1.84 1.61 1.58 1.61 1.63 1.64 1.63 糠醇
furfuryl alcohol0.99 1.00 1.01 1.03 1.04 1.05 0.95 0.94 0.94 0.96 0.97 0.96 5-甲基糠醛
5-methylfurfural0.13 0.12 0.13 0.13 0.14 0.12 0.11 0.12 0.11 0.11 0.12 0.11 苯丙氨酸
phenylalanine苯甲醛
benzaldehyde0.27 0.25 0.23 0.25 0.25 0.26 0.21 0.20 0.20 0.20 0.19 0.20 苯甲醇
benzyl alcohol6.63 6.53 6.48 6.62 6.45 6.70 5.46 5.58 5.34 5.52 5.71 5.26 苯乙醛
phenylacetaldehyde0.76 0.75 0.77 0.77 0.76 0.80 0.60 0.65 0.65 0.65 0.68 0.66 苯乙醇
benzyl ethanol4.42 3.88 4.09 4.30 4.31 4.58 3.52 3.44 3.59 3.57 3.35 3.77 西柏烷类
cembranoids茄酮
solanone16.56 16.40 16.72 16.79 16.32 17.12 15.56 15.65 15.82 15.91 15.80 15.93 螺岩兰草酮
solavetivone3.67 3.63 3.70 3.77 3.79 3.89 4.75 4.65 4.72 4.68 4.60 4.46 降茄二酮
norsolandione0.89 0.91 0.91 0.92 0.92 0.94 1.00 0.99 1.03 1.01 0.98 1.03 胡萝卜素
carotenoidsβ-大马酮
β-damascenone6.83 6.77 6.58 6.84 6.66 7.02 6.53 6.59 6.68 6.75 6.89 6.71 β-二氢大马酮
β-2H-damascenone2.45 2.41 2.33 2.43 2.39 2.44 2.48 2.48 2.48 2.56 2.58 2.67 β-紫罗兰酮
β-ionone1.38 1.52 1.38 1.43 1.34 1.46 1.27 1.25 1.26 1.26 1.30 1.18 3-氧代-α-紫罗兰醇
3- oxo-α-ionol0.84 0.84 0.89 0.86 0.87 0.83 0.83 0.84 0.77 0.81 0.78 0.79 香叶基丙酮
geranyl acetone3.88 3.94 3.86 3.87 3.58 3.91 3.30 3.30 3.40 3.40 3.36 3.47 二氢猕猴桃内酯
dihydroactinidiolide2.34 2.54 2.61 2.55 2.60 2.56 2.72 2.79 2.65 2.75 2.75 2.71 巨豆三烯酮A
megastigmatrienone A2.95 3.08 2.93 3.00 2.90 3.00 2.99 2.91 3.06 3.04 2.91 3.17 巨豆三烯酮B
megastigmatrienone B12.75 12.71 11.81 12.68 12.55 12.97 13.23 12.99 13.67 13.38 13.29 13.07 巨豆三烯酮C
megastigmatrienone C2.51 2.52 2.67 2.64 2.69 2.70 2.80 2.81 2.87 2.90 2.80 3.07 巨豆三烯酮D
megastigmatrienone D11.87 11.93 11.75 12.22 12.31 12.64 12.61 12.84 12.85 13.08 12.91 13.52 金合欢基丙酮
farnesylacetone12.93 13.18 12.78 13.10 12.89 13.10 12.05 12.26 12.74 12.50 12.42 12.39 叶绿素
chlorophyll新植二烯
neophytadiene544.94 539.37 528.05 545.45 552.11 535.80 502.29 493.36 506.64 507.47 498.41 511.53 表 2 不同切丝宽度烟丝挥发性香味物质的多重比较结果
Table 2 Multiple comparison results of volatile aroma components of cut tobacco of different cutting width
挥发性香味物质
volatile aroma components干燥前烟丝
cut tobacco before drying不同宽度烟丝干燥后
cut tobacco of different cutting width after drying0.85 mm 0.95 mm 1.05 mm 糠醛
furfural1.80 Aa 1.62 Bb 1.64 Bb 1.08 Cc 糠醇
furfuryl alcohol1.02 Aa 0.95 Bb 0.92 Bc 0.57 Cd 5-甲基糠醛
5- methylfurfural0.13 Aa 0.11 Bb 0.10 Bc 0.08 Cd 苯甲醛
benzaldehyde0.25 Aa 0.20 Bb 0.19 Cc 0.17 Dd 苯甲醇
benzyl alcohol6.57 Aa 5.48 Bb 4.65 Cc 3.30 Dd 苯乙醛
phenylacetaldehyde0.77 Aa 0.65 Bb 0.44 Cc 0.29 Dd 苯乙醇
benzyl ethanol4.26 Aa 3.54 Bb 3.19 Cc 2.53 Dd 茄酮
solanone16.65 Aa 15.78 Bb 15.16 Cc 13.67 Dd 螺岩兰草酮
solavetivone3.74 Cc 4.64 Aa 4.00 Bb 2.59 Dd 降茄二酮
norsolandione0.91 Bb 1.01 Aa 0.91 Bb 0.68 Cc β-大马酮
β-damascenone6.78 Aa 6.69 Aa 6.07 Bb 4.83 Cc β-二氢大马酮
β-2H-damascenone2.41 Bb 2.54 Aa 2.22 Cc 1.74 Dd β-紫罗兰酮
β-ionone1.42 Aa 1.25 Bb 1.24 BCb 1.18 Cc 3-氧代-α-紫罗兰醇
3- oxo-α-ionol0.85 Aa 0.80 Ab 0.83 Aab 0.55 Bc 香叶基丙酮
geranyl acetone3.84 Aa 3.37 Bb 2.90 Cc 3.00 Cc 二氢猕猴桃内酯
dihydroactinidiolide2.53 Bb 2.73 Aa 2.48 Bb 1.92 Cc 巨豆三烯酮A
megastigmatrienone A2.98 Aab 3.01 Aa 2.87 Ab 2.56 Bc 巨豆三烯酮B
megastigmatrienone B12.58 Bb 13.27 Aa 11.52 Cc 10.57 Dd 巨豆三烯酮C
megastigmatrienone C2.62 Bb 2.88 Aa 2.47 Cc 1.95 Dd 巨豆三烯酮D
megastigmatrienone D12.12 Bb 12.97 Aa 11.68 Bb 10.04 Cc 金合欢基丙酮
farnesylacetone13.00 Aa 12.40 Ab 11.03 Bc 9.26 Cd 新植二烯
neophytadiene540.95 Aa 503.28 Bb 467.38 Cc 402.99 Dd 注:表中同一行数字后面小写字母、大写字母不同的分别表示差异达到0.05和0.01显著水平。
Note: Different small letters and capital letters in same row mean significant difference at the level of 0.05 and 0.01.表 3 挥发性香味物质主成分分析前的KMO及Bartlett球形度检验
Table 3 KMO and Bartlett sphericity test before principal component analysis of volatile aroma components
化学指标 chemical index KMO Bartlett球形度检验 Bartlett sphericity test 近似卡方 approximate chi square 自由度 degree of freedom 显著性水平 significant level 挥发性香味物质 volatile aroma component 0.511 1 536.906 231 0.000 表 4 挥发性香味物质主成分分析的特征值及方差贡献率
Table 4 Eigenvalue and variance contribution rate of principal component analysis of volatile aroma components
主成分
principal component特征值
eigenvalue方差贡献率/%
variance contribution rate累积方差贡献率/%
cumulative variance contribution ratePC1 18.200 82.729 93.244 PC2 2.313 10.515 表 5 前两个主成分的载荷矩阵和特征向量
Table 5 Load matrix and feature vector of the first two principal components
挥发性香味物质
volatile aroma components主成分载荷矩阵
load matrix of principal component主成分特征向量
feature vector of principal componentPC1 PC2 PC1 PC2 糠醛 furfural 0.943 −0.023 0.221 −0.015 糠醇 furfuryl alcohol 0.966 −0.082 0.226 −0.054 5−甲基糠醛 5− methylfurfural 0.953 0.189 0.223 0.124 苯甲醛 benzaldehyde 0.815 0.550 0.191 0.362 苯甲醇 benzyl alcohol 0.959 0.257 0.225 0.169 苯乙醛 phenylacetaldehyde 0.951 0.232 0.223 0.153 苯乙醇 benzyl ethanol 0.922 0.321 0.216 0.211 茄酮 solanone 0.972 0.172 0.228 0.113 螺岩兰草酮 solavetivone 0.812 −0.562 0.190 −0.369 降茄二酮 norsolandione 0.903 −0.411 0.212 −0.270 β−大马酮 β−damascenone 0.989 −0.058 0.232 −0.038 β−二氢大马酮 β−2H−damascenone 0.949 −0.237 0.222 −0.156 β−紫罗兰酮 β−ionone 0.731 0.581 0.171 0.382 3−氧代−α−紫罗兰醇 3− oxo−α−ionol 0.730 0.583 0.171 0.383 香叶基丙酮 eranyl acetone 0.902 −0.365 0.211 −0.240 二氢猕猴桃内酯 dihydroactinidiolide 0.904 −0.185 0.212 −0.122 巨豆三烯酮A megastigmatrienone A 0.891 −0.186 0.209 −0.122 巨豆三烯酮B megastigmatrienone B 0.917 −0.350 0.215 −0.230 巨豆三烯酮C megastigmatrienone C 0.904 −0.348 0.212 −0.229 巨豆三烯酮D megastigmatrienone D 0.897 −0.134 0.210 −0.088 金合欢基丙酮 farnesylacetone 0.971 0.067 0.228 0.044 新植二烯 neophytadiene 0.966 0.160 0.226 0.105 表 6 不同试验烟支样品感官质量的三点法检验结果
Table 6 Test result of sensory evaluation of different test cigarette samples by three point method
样品对子
sample pair参与评吸人数①
number of smoking正确识别人数
number of correct identification正确识别率/%
correct recognition rate整体品质特征差异
characteristic difference of whole quality0.85 mm与0.95 mm 9 4 44.44 不显著 0.95 mm与1.05 mm 9 7 77.78 显著 0.85 mm与1.05 mm 9 8 88.89 显著 注:① α=0.05时,三点检验临界值为6。
Note: ① Critical value of three point test is six when α is equal to 0.05.表 7 不同试验烟支样品的感官质量变化趋势
Table 7 Sensory quality changes of different test cigarette samples
样品编号
sample number香气特性 aroma characteristics 烟气特性 smoke characteristics 口感特性 taste characteristics 香气质
aroma quality香气量
aroma quantity杂气
offensive odor劲头
vigor浓度
concentration细腻程度
delicate degree刺激性
irritation干燥感
dry sense干净程度
clean degreeCK (0.85 mm) 0 0 0 0 0 0 0 0 0 0.95 mm −1 0 −1 0 +1 −1 −1 0 0 1.05 mm −2 +1 −1 +1 +2 −2 −2 −1 −1 注:表中的数值分别表示与对照样品相比指标质量的变化方向及变化程度。
Note: Different values in the table respectively express change direction and degree of quality indexes compared with the control sample. -
[1] 罗登山, 姚光明, 刘朝贤. 中式卷烟加工工艺技术探讨[J]. 烟草科技, 2005(5): 4. DOI: 10.1007/s12028-009-9223-y. [2] 卢彦华, 于录, 阮晓明, 等. 分组加工工艺过程研究[J]. 烟草科技, 2008(1): 12. DOI: 10.3969/j.issn.1002-0861.2008.01.004. [3] 魏玉玲, 阴耕耘, 李绍臣, 等. 几个重要制丝工序对烤烟烟丝填充值和碎丝率的影响[J]. 云南大学学报(自然科学版), 2010, 32(S1): 183. [4] 于存峰, 于录, 杨雷, 等. 模块加工中叶丝宽度组合对卷烟质量的影响[J]. 烟草科技, 2011(6): 9. DOI: 10.3969/j.issn.1002-0861.2011.06.002. [5] 邱玉春, 林志平, 黄朝章, 等. 切丝宽度对卷烟主流烟气7种Hoffmann成分释放量的影响[J]. 郑州轻工业学院学报(自然科学版), 2013, 28(6): 9. DOI: 10.3969/j.issn.2095-476X.2013.06.003. [6] 薛芳, 李东亮, 陈昆燕, 等. 卷烟加工重点工序工艺参数与卷烟主流烟气中苯酚释放量的关系研究[J]. 江西农业大学学报, 2010, 32(6): 1307. DOI: 10.3969/j.issn.1000-2286.2010.06.040. [7] 谭兰兰, 施丰成, 薛芳, 等. 卷烟制丝工艺参数对主流烟气中苯并[a]芘释放量的影响[J]. 安徽农业科学, 2010, 38(1): 165. [8] 张涛, 苏明亮, 赵伟, 等. 制丝工艺参数对9 种主流烟气成分释放量的影响[J]. 烟草科技, 2014(8): 32. DOI: 10.3969/j.issn.1002-0861.2014.08.008. [9] 田忠, 陈闯, 许宗保, 等. 制丝关键工序对细支卷烟燃烧温度及主流烟气成分的影响[J]. 中国烟草学报, 2015, 21(6): 19. DOI: 10.16472/j.issn.2015.107. [10] 王廷, 李斌, 朱文魁. 叶丝滚筒干燥过程中化学成分的动态变化特征研究[J]. 烟草科技, 2010(11): 9. DOI: 10.3969/j.issn.1002-0861.2010.11.002. [11] 王宗英, 王建民, 王海妮, 等. 基于均匀设计分析滚筒烘丝工艺参数对叶丝质量的影响[J]. 烟草科技, 2011(8): 5. DOI: 10.3969/j.issn.1002-0861.2011.08.001. [12] 张强, 董高峰, 李红武, 等. 滚筒烘丝机工艺参数对烤烟感官质量的影响[J]. 烟草科技, 2011(11): 10. DOI: 10.3969/j.issn.1002-0861.2011.11.002. [13] 刘泽, 何邦华, 杨蕾, 等. " 干头干尾”烟丝挥发性香味物质的主成分聚类分析[J]. 烟草科技, 2015(2): 47. DOI: 10.16135/j.issn1002-0861.20150210. [14] 李鹏飞, 周冀衡, 张建平, 等. 烤烟成熟期土壤水分状况对烟叶挥发性香气物质及主要化学成分的影响[J]. 中国烟草学报, 2009, 15(3): 44. DOI: 10.3969/j.issn.1004-5708.2009.03.009. [15] 薛超群, 王建伟, 奚家勤, 等. Agri土壤调理剂用量对烟叶香气质量的影响[J]. 烟草科技, 2012(7): 86. DOI: 10.3969/j.issn.1002-0861.2012.07.019. [16] 阴耕耘, 徐世涛, 阮昕, 等. 复烤过程中B4F烤烟烟叶类胡萝卜素和西柏烷类降解产物的变化[J]. 烟草科技, 2012(10): 61. DOI: 10.3969/j.issn.1002-0861.2012.10.014. [17] GB/T 17321—1998. 感官分析方法 二、三点检验[S]. [18] YC/T 138—1998. 烟草及烟草制品 感官评价方法[S]. [19] 时立文. SPSS 19.0统计分析从入门到精通[M]. 北京: 清华大学出版社, 2013: 290. [20] 汤银才. R语言与统计分析[M]. 北京: 高等教育出版社, 2008: 332. [21] 于建军. 卷烟工艺学[M]. 北京: 中国农业出版社, 2009: 179. [22] 张梦楚, 王龙, 朱波, 等. 红花大金元烟叶中性致香成分含量与感官质量预测分析[J]. 湖北农业科学, 2014, 53(22): 5438. DOI: 10.14088/j.cnki.issn0439-8114.2014.22.027. -
期刊类型引用(1)
1. 易志飞,付春,江纳,杨瑶君,肖开兴. 枯草芽孢杆菌表面活性素和伊枯草菌素的生物信息学分析. 基因组学与应用生物学. 2021(Z1): 2177-2200 . 百度学术
其他类型引用(0)