基于代谢组学分析明火熏烤烟叶的代谢差异及形成机制
运用代谢组学技术分析明火熏烤烟叶的代谢差异及其形成机制。
以烤烟品种K326为试验材料，进行明火熏烤(X)和常规烘烤(CK)对比试验，利用代谢组学技术分析2种调制方法的代谢差异，并探讨其形成机制。
KEGG富集结果表明：在X-38 vs CK-38比较组中，亚油酸代谢和植物激素信号转导是富集最显著的2条通路；在X-42 vs CK-42比较组中，糖酵解/糖异生和氨基酸的生物合成是富集最显著的2条通路；在CK-42 vs CK-38比较组和X-42 vs X-38比较组中，亚油酸代谢和核苷酸代谢是富集最显著的2条通路。在4个比较组中，上调表达的代谢物包括游离脂肪酸、氨基酸及其衍生物、酚酸类、糖类、有机酸、生物碱、核苷酸及其衍生物等。
明火熏烤烟叶的差异代谢物包含大量酸性代谢物。形成差异的原因是明火熏烤使烟叶额外受到烟雾和光热辐射胁迫，促使其形成更多的代谢物。
Metabolomic Analysis of Metabolic Differences and Formation Mechanisms in Open Fire-cured Tobacco Leaves
To analyze the metabolic differences and formation mechanisms in open fire-cured tobacco leaves by metabolomics techniques.
Tobacco cultivar K326 was used as the test material for a comparative study between open fire-curing (X) and conventional curing (CK). Metabolomic techniques were employed to analyze the metabolic differences between the two curing methods and to explore their formation mechanisms.
KEGG enrichment results indicated that: in the X-38 vs CK-38 comparison groups, linoleic acid metabolism and plant hormone signal transduction were the two most significantly enriched pathways; in the X-42 vs CK-42 comparison groups, glycolysis/gluconeogenesis and amino acid biosynthesis were the two most significantly enriched pathways; in the CK-42 vs CK-38 and X-42 vs X-38 comparison groups, linoleic acid metabolism and nucleotide metabolism were the two most significantly enriched pathways. The up-regulated metabolites in all four comparison groups included free fatty acids, amino acids and their derivatives, phenolic acids, sugars, organic acids, alkaloids, nucleotides and their derivatives.
The differential metabolites in open fire-cured tobacco leaves include a large number of acidic metabolites. The formation of these differences is due to additional environmental stresses such as smoke and thermal radiation experienced during open fire-curing, which lead to the production of more metabolites in open fire-cured tobacco leaves.
-
Keywords:
- metabolomics /
- open fire-curing /
- KEGG /
- differential metabolites
-
猪流行性腹泻病毒(Porcine epidemic diarrhea virus，PEDV)是一种猪肠道致病性冠状病毒，其引起的猪流行性腹泻是一种以严重呕吐、腹泻和脱水为特征的传染性疾病[1]。该病毒可通过消化道和呼吸道传播；能感染各年龄段的猪，对1周龄内的仔猪危害最为严重，可造成100%死亡[2]。临床上，PEDV多与猪传染性胃肠炎病毒、猪轮状病毒和猪圆环病毒等病毒混合感染[3]。因此，该病在许多国家造成了巨大的经济损失，并对全球生猪养殖业造成巨大威胁。
PEDV属于冠状病毒科(Coronaviridae)甲型冠状病毒属(Alphacoronavirus)成员，为有囊膜单股正链RNA病毒[4]。基因组全长约28 kb，编码4种结构蛋白[ 纤突糖蛋白(S)、基质蛋白(M)、囊膜蛋白(E)和核蛋白(N)]、1个辅助蛋白(ORF3)和2个非结构蛋白 (复制酶多聚蛋白pp1a和pp1b)。PEDV结构蛋白在介导病毒入侵、中和抗体产生、病毒基因组复制和转录以及病毒组装、出芽和拮抗宿主先天免疫反应中发挥着重要作用[5-7]；PEDV的非结构蛋白在病毒基因表达、RNA合成加工以及病毒粒子组装和释放等方面发挥着重要作用[8]。
通过将PEDV S蛋白基因进行遗传进化分析，可将PEDV分为传统PEDV毒株和新型PEDV毒株[9]，新型PEDV毒株又分为S non-INDEL型和S INDEL型2个亚群[10]。在生产中，S non-INDEL型导致的病死率通常更高[11]。S non-INDEL型毒株是近年来流行的优势毒株，其基因变异较快，并引起抗原表位发生变化，是导致当前商品化疫苗免疫效果不佳的重要原因之一。本研究病料组织采集于接种过PEDV疫苗的规模化猪场，利用该病料进行病毒分离鉴定、测定分离株的全基因序列并分析该毒株与疫苗株之间的遗传进化关系，以期为了解PEDV流行株的变异情况及疫苗开发提供参考。
1. 材料与方法
1.1 病料样本与细胞株
15份仔猪小肠组织样本采自接种过PEDV疫苗的四川某猪场；PEDV、TGEV (猪传染性胃肠炎病毒)和PoRV (猪轮状病毒)阳性血清以及猪肾细胞(PK-15)均由四川农业大学动检实验室保存。
1.2 主要试剂及试验动物
细胞培养液(DMEM)、胎牛血清(FBS)、胰蛋白酶和氨苄青霉素(Amp)等购自美国Gibco公司；二甲基亚砜(DMSO)、低熔点琼脂糖、Hanks液(10×)和DL2000 Marker购自上海生工生物工程有限公司；病毒RNA提取试剂盒购自北京天根生化科技有限公司；反转录试剂盒和2×Taq Master Mix购自大连宝生物工程有限公司。
6头1周龄哺乳仔猪(未免疫)由四川农业大学合作猪场提供。
1.3 引物设计
PEDV检测引物参考《猪流行性腹泻诊断技术》(NY/T 544—2015)[12]；TGEV和PoRV检测引物参考郭容利等[13]和卓秀萍[14]的研究(表1)；根据PEDV CV777株(登录号：AF353511)、CH/S株(登录号：JN547228)及DR13株(登录号：JQ023161)的全基因序列，运用Primer 5.0软件设计22对引物扩增分离株全基因组序列。以上引物均由上海生工生物工程有限公司合成。
表 1 鉴别引物序列Table 1. The identified primers sequences扩增片段基因
gene fragment引物名称
name of primers引物序列 (5′→3′)
primers sequence退火温度/℃
anneal temperature扩增片段长度/bp
products lengthPEDV-ORF3 PEDV1-F CCTAGACTTCAACCTTACGA 56 774 PEDV1-R CAGGAAAAAGAGTACGAAAA TGEV-M TGEV1-F GCGTCTGATTGTGAGTCATG 54 554 TGEV1-R AATAGTCCTGCTGGGTAATG PoRV-VP6 PoRV1-F CTATGTTGAAAATGCTCGAAC 52 229 PoRV1-R GAATACGAATCCAGTACGC G-R TACAACAAGCCCAATC 1.4 组织样本PCR鉴定
提取15份组织样本的总RNA并反转录成cDNA。根据《猪流行性腹泻诊断技术》(NY/T 544—2015)[12]及郭容利等[13]和卓秀萍[14]设计的TGEV与PoRV 检测方法进行RT-PCR扩增，PCR产物经1%琼脂糖凝胶进行电泳，对15份样品中PEDV的感染情况进行分析。
1.5 病毒分离、纯化及病毒滴度测定
1.5.1 病毒分离
取PEDV阳性小肠组织1 g制成原代病毒液，分装后于−80 ℃保存备用。将1 mL原代病毒液接种至状态良好的PK-15细胞，于37 ℃孵育1 h后吸弃病毒液，添加胰酶终质量浓度为4 μg/mL的维持液5 mL，待细胞病变达80%~90%时，将细胞置于液氮中冻融3次后离心取上清作为传代接种液。按上述方法进行病毒传代，并观察细胞病变情况。
1.5.2 病毒纯化
将1 mL病毒传代接种液用DMEM进行5个系列的10倍梯度稀释(10−1~10−5)，后接种单层PK-15细胞，另有一孔加1 mL维持液作空白对照。于培养箱孵育1 h后吸弃传代接种液，每孔添加2 mL含1%低熔点琼脂糖的DMEM (2×)覆盖液，铺平后继续培养。当出现细胞病变时，每孔添加2 mL含0.05%结晶紫的覆盖液，继续培养。眼观无色斑点为蚀斑，镜下观察蚀斑具体位置，用灭菌10 μL枪头挑取单个蚀斑进行扩大培养。
为进行分离株的RT-PCR检测，分别提取第3、7、10代及纯化病毒的细胞上清液总RNA，反转录成cDNA后进行PCR鉴定。PCR反应体系为25 μL，包括：Premix Taq 12 μL，RNase Free ddH2O 10 μL，上、下游引物各0.5 μL，模板2 μL。反应条件为：95 ℃预变性3 min；94 ℃变性40 s，56 ℃退火40 s，72 ℃延伸45 s，共35个循环；72 ℃延伸10 min。反应结束后，取7.5 μL PCR产物进行电泳。
1.5.3 病毒滴度测定
将纯化的病毒用DMEM进行10个系列的10倍梯度稀释(10−1~10−10)，后接种于长成单层细胞的96孔板中，每个稀释度设8孔，每孔接种100 μL；最后2列接种维持液作阴性对照。将细胞继续培养，3 d后观察并记录每孔细胞病变情况。按照Reed-Muench法计算病毒滴度(TCID50)。
1.6 仔猪感染试验及临床症状和病理学观察
将6头1周龄健康哺乳仔猪(未免疫)分为试验组和对照组，每组3头。饲养环境相同，分别隔离饲养。试验组将10 mL纯化PEDV病毒液拌料饲喂仔猪；对照组将10 mL细胞培养液以同样的方式饲喂仔猪。每隔4 h观察仔猪临床症状。将感染后死亡的仔猪剖检，观察各组织器官病理变化；取病变明显的肺脏、肝脏及小肠组织制作病理切片，利用实验室建立的方法测定各组织器官病毒载量。
1.7 分离株的全基因组测序及系统发育进化树构建
提取PEDV SNJ-P株的总RNA并反转录为cDNA。使用1.3节中设计的全基因组扩增引物扩增分离株全基因序列，并回收纯化。将纯化的DNA片段连接到pMD19-T (simple)载体，并转入DH5α感受态细胞进行克隆扩增。经PCR鉴定为阳性的菌液送至上海生工生物工程有限公司测序。使用MEGA 6.0软件将获得的22条序列进行拼接，获得分离株全基因序列。基于NCBI上38株国内外PEDV毒株的S蛋白基因序列信息，运用MegAlign软件构建系统发育进化树，并进行遗传进化分析。以传统疫苗株CV777为参考，分析SNJ-P株的各基因序列及S蛋白氨基酸序列遗传变异情况。
2. 结果与分析
2.1 病料PCR鉴定结果
由图1可知：PEDV、TGEV和PoRV的阴性和阳性对照均成立，阳性对照条带大小分别为774、554和229 bp，与预期相符。15份样本中仅有1份样本检测PEDV呈阳性，且该份样本TGEV和PoRV检测均呈阴性。PCR鉴定出的PEDV阳性样本为后续毒株的分离提供了材料。
![]()
图 1 病料RT-PCR检测结果注：M. Marker DL2000；1、4和7. PEDV、TGEV和PoRV阳性对照；2、5和8. PEDV、TGEV和PoRV阴性对照；3、6和9. PEDV阳性病料的PEDV、TGEV和PoRV检测结果。Figure 1. The results of RT-PCR detection of disease materialNote: M. Marker DL2000; 1, 4 and 7. PEDV, TGEV and PoRV positive control; 2, 5 and 8. PEDV, TGEV and PoRV negative control; 3, 6 and 9. detection results of PEDV, TGEV and PoRV in PEDV positive samples.2.2 病毒分离、纯化及病毒滴度测定结果
由图2可知：病毒感染后，盲传至第4代开始出现明显的细胞病变效应(cytopathic effect， CPE)；盲传至第7代，细胞产生稳定的CPE；盲传至第10代，细胞产生CPE的时间提前，在接毒后约12 h即观察到CPE产生，24~36 h CPE可达80%以上。细胞病变主要表现为细胞内颗粒增多、细胞圆缩脱落、大面积溶解和细胞间隙增大等。
![]()
图 2 PK-15细胞接种病毒液后各代病变情况注：a) 空白对照；b) 盲传第4代48 h；c) 盲传第7代48 h；d) 盲传第10代12 h。Figure 2. The pathological changes of different generations PK-15 cells after inoculation with virus fluidNote: a) blank control; b) 48 hours after blind passage the 4th generation; c) 48 hours after blind passage the 7th generation; d) 12 hours after blind passage the 10th generation.由图3可知：接毒后12 h细胞开始产生病变；3 d后覆盖含有结晶紫的覆盖液，2~3 d可以观察到蚀斑，空白对照组无蚀斑产生。
![]()
图 3 细胞接种不同稀释度病毒液后产生的蚀斑注：a) 细胞接种病毒液(稀释度为10−4)产生的蚀斑染色；b) 细胞接种病毒液(稀释度为 10−5)产生的蚀斑染色；c) 空白对照。Figure 3. Plaque formation after inoculation of different dilutions of virus fluidNote: a) plaque staining result of inoculate 10−4 virus solution; b) plaque staining result of inoculate 10−5 virus solution; c) blank control.由图4可知：第3、7、10代以及纯化病毒液PEDV检测均呈阳性，且条带亮度随传代数的增加而增大。这表明随传代数的增加，病毒增殖能力不断增强，不断适应培养环境。按照Reed-Muench法计算病毒滴度(TCID50)为1×107.5/100 μL。
![]()
图 4 各代病毒液RT-PCR检测结果注：M. DL2000 DNA Marker；1~3. 第3、7、10代病毒液；4. 纯化病毒液；5. 阴性对照；6. 阳性对照。Figure 4. The results of RT-PCR detection in each generation of virus fluidNote: M. DL2000 DNA Marker; 1-3. the 3rd, 7th and 10th generation virus liquid; 4. purified virus fluid; 5. negative control; 6. positive control.2.3 仔猪感染临床症状和病理变化
临床表现：试验组仔猪在感染后6 h开始出现精神沉郁；12 h后陆续出现腹泻症状；30~36 h后出现食欲减退、水样腹泻、脱水及黏膜苍白等症状；感染后40~48 h，3只仔猪均死亡。对照组仔猪表现正常。
病理学观察：病死仔猪消化道病变最严重，胃和小肠内有未消化的乳糜，并发酵产生大量气体；小肠壁薄而透明，内充满黄绿色水样内容物，气味腥臭；肠系膜充血，肿胀严重；其他组织器官病变不明显，脾脏轻微血肿，底部有淤血；肺脏轻微肿胀，右肺叶出现淤血点；心脏局部水肿；肾脏与肝脏无明显病变。
取病变明显的肺脏、肝脏和小肠组织制作病理切片进行观察发现：肺泡壁增厚，细胞成分增多，肺泡出现萎陷融合现象(图5a)；局部肝细胞肿胀，肝窦变窄或消失，肝索纹理不清，含铁血黄素沉积(图5b)；空肠绒毛高度普遍降低，并出现绒毛萎缩、断裂和脱落等现象(图5c)；十二指肠绒毛高度普遍降低，严重萎缩(图5d)。
![]()
图 5 病死仔猪肺脏、肝脏、空肠及十二指肠组织病理学变化(H.E，400×)注：a) 肺脏；b) 肝脏；c) 空肠；d) 十二指肠；箭头所指为病理变化位置。Figure 5. Pathological changes of lung, liver, jejunum and duodenum in dead pigletsNote: a) lung; b) liver; c) jejunum; d) duodenum; the arrow points to the location of pathological changes.利用本实验室建立的qPCR方法检测病死仔猪各组织器官的病毒载量，结果表明：肠道组织病毒载量最高，为3.82×106 copies/mL；其次是肺脏，为1.47×105 copies/mL；心脏、肝脏、脾脏和肾脏等组织器官病毒拷贝数均较低，分别为3.03×103、1.31×104、3.00×103和1.09×103 copies/mL。该毒株具有较强的肠道嗜性，符合PEDV感染特点。以上结果表明成功分离获得1株PEDV毒株，并将其命名为PEDV SNJ-P。
2.4 分离株的全基因组测序及遗传进化分析
PEDV SNJ-P株全基因共28003 bp，序列已上传NCBI，登录号为MK702008.1。将分离株与不同地区毒株及疫苗株的S蛋白基因进行遗传进化分析，结果(图6)表明：分离株属于S non-INDEL型，与同型的中国分离株PEDV-WS (登录号：KM609213.1)、CH/JLDH/2016 (登录号：MF346935.1)和CH/HNLH/2015 (登录号：KT199103.1)亲缘关系较近，与同型的疫苗株AJ1102-F12 (登录号：MK584552.1)和LW/L (登录号：MK392335.1)亲缘关系相对较远，与传统疫苗株CV777 (登录号：LT905450.1)和DR13 (登录号：JQ023162.1)亲缘关系较远。
![]()
图 6 基于SNJ-P株与参考株S蛋白基因构建的系统发育进化树Figure 6. Phylogenetic tree based on S protein gene construction of SNJ-P strain and reference strain2.5 各序列片段遗传与变异情况
SNJ-P株与传统疫苗CV777的全基因序列比对结果表明：SNJ-P株全基因序列碱基共发生980处突变、16处插入和11处缺失，两者的碱基变化数如表2所示。
表 2 SNJ-P株与CV777株基因组碱基变化数对比Table 2. Comparison of nucleotide changes between SNJ-P strain and CV777 strain毒株
strains基因组全长/bp
full length of genome5′非翻译区
5′UTR开放阅读框1
ORF1纤突糖蛋白
S开放阅读框3
ORF3囊膜蛋白
E基质蛋白
M核蛋白
N3′非翻译区
3′UTRSNJ-P 28003 261 20345 4161 675 231 681 1326 330 CV777 28033 296 20345 4152 675 231 681 1326 334 注/Note: UTR. untranslated region; ORF. open reading frame; S. spike glycoprotein; E. envelope protein; M. matrix protein; N. nucleoprotein. 2.6 S蛋白氨基酸位点变化情况
由图7可知：SNJ-P株S蛋白基因共有4161个碱基，与CV777株相比其碱基存在较多变化，主要集中在S1区 (1~789 aa)，共发生271处突变、6处缺失和15处插入。164~166 nt 插入TTG等3个碱基，177~185 nt 插入AGGTGTTAA等9个碱基，415~417 nt 插入GAT等3个碱基，分别导致55~56 aa、59~62 aa及139 aa位插入氨基酸IG (异亮氨酸和甘氨酸)、QGVN (谷氨酰胺、甘氨酸、缬氨酸和天冬酰胺) 和D (天冬氨酸)。480~486 nt缺失TGGAAA等6个碱基，导致160~162 aa位缺失氨基酸DGK (天冬氨酸、甘氨酸和赖氨酸)。79~87 nt由CAGTCTACT突变为TCAGCTAAC，导致27~29 aa位氨基酸由QST (谷氨酰胺、丝氨酸和苏氨酸)变为SAN (丝氨酸、丙氨酸和天冬酰胺) ；203~215 nt由GCACAGGCATTGA突变为CTGGCCAACATCC，导致68~72 aa位氨基酸由GTGIE (甘氨酸、苏氨酸、甘氨酸、异亮氨酸和谷氨酸) 变为AGQHP (丙氨酸、甘氨酸、谷氨酰胺、组氨酸和脯氨酸) ；601~608 nt位由CGCAGAAGT突变为AGTGGAGGT，导致201~203 aa位氨基酸由SGG (丝氨酸、甘氨酸和甘氨酸)变为RRS (精氨酸、精氨酸和丝氨酸)。其余突变位点均是由单碱基突变导致单个或2个间断氨基酸发生改变。S蛋白的插入与缺失均发生在S1区，同时S1区有63处氨基酸位点发生突变，S2区有19处氨基酸位点发生突变。
![]()
图 7 SNJ-P与CV777的S蛋白氨基酸序列比对Figure 7. Comparison of amino acid sequence of S protein between SNJ-P and CV7773. 讨论
猪流行性腹泻病毒是引起仔猪腹泻、脱水并导致死亡的重要病原之一，在中国猪群中传播时间长、波及范围广，对生猪养殖业造成了巨大的经济损失。PEDV毒株的分离仍然是病毒学研究及疫苗开发的重要内容。该病毒毒株在细胞增殖过程中存在细胞适应能力弱、易出现病毒滴度降解等问题[15]。因此，探索能使病毒稳定增殖的细胞系及培养条件，对于流行毒株的分离尤为重要。研究表明：PEDV可在多数猪源细胞系中生长，包括猪膀胱细胞[16]、猪肾细胞系[17]、猪睾丸细胞系[18]、猪肺泡巨噬细胞系[19]以及猪小肠上皮细胞系[20]等，但目前猪源细胞系分离PEDV流行株的培养条件仍需进一步探索。 1988年，HOFFMAN通过在Vero细胞培养基中加入胰酶成功分离获得1株PEDV毒株；随着学者对PEDV更加深入的研究，Vero细胞成为分离PEDV最常用的细胞种类，相关的试验方法也趋于完善，但Vero细胞系分离PEDV毒株存在增殖不稳定、滴度易降低等问题[20]。PK-15细胞具有生长速度快、环境耐受力强等特点，且可以用于体外分离PEDV毒株[21]，但相关记录较少，试验方法并不详尽。因此，本研究尝试建立以PK-15细胞为基础的毒株分离条件，利用PK-15细胞进行病毒分离，并在其维持培养基中添加终质量浓度为4 μg/mL的胰酶，盲传至第4代即产生明显的病变，传至第7代后产生稳定CPE。
PEDV能感染各年龄段的猪，尤以1周龄内的仔猪病死率最高。究其原因，一方面与初生仔猪肠绒毛更新修复能力较弱有关；另一方面与初生仔猪免疫系统不健全、抵抗病原入侵的能力较弱有关[22]。本研究利用SNJ-P株进行动物感染试验，结果显示：仔猪出现呕吐、腹泻和脱水等典型临床症状。病毒载量检测结果显示：小肠病毒载量最高，表明此分离毒株具有较高的肠道组织嗜性。以上研究结果与周书亭等[23]和卓秀萍等[24]的研究结果基本一致。
将分离株与不同地区毒株及疫苗株的S蛋白基因进行遗传进化分析，结果表明：分离株属于S non-INDEL型，与同型的中国分离株PEDV-WS、CH/JLDH/2016和CH/HNLH/2015亲缘关系较近，与同型的疫苗株AJ1102-F12和LW/L亲缘关系相对较远，与传统疫苗株CV777和DR13亲缘关系较远。通过对SNJ-P株与CV777的S蛋白氨基酸序列进行比对发现：与CV777相比，SNJ-P株S蛋白S1亚基存在多处氨基酸突变、插入和缺失等变化，导致SNJ-P株相应的抗原表位发生改变。以上结果从遗传进化层面说明传统毒株制备的疫苗免疫效果不佳，可能与当前流行毒株基因组变异较快有关，从而导致相应抗原表位发生改变。
目前，针对猪流行性腹泻尚无特效药，感染后只能通过隔离猪群或紧急免疫降低猪群的死亡率。疫苗免疫仍然是预防该病的主要措施，现有的疫苗主要以灭活苗和弱毒苗为主[25]。在中国，早在1999年就开发出针对PEDV经典毒株的疫苗[26]；近年来，针对S non-INDEI型毒株的商品化疫苗不断上市，如基于AJ1102株的PEDV和TGEV二价灭活疫苗等[27]，为PEDV的防控提供了一定的帮助。由于病毒变异较快，现有的疫苗保护力有限，一些规模化猪场通过制备自家组织灭活苗，免疫母猪产生抗体，仔猪通过吮吸初乳获得抗体产生免疫保护，在一定程度上提高了仔猪存活率。但自家组织灭活苗同样具有一定的局限性，如：仅针对本场常发毒株，不适合大规模推广应用；仅在疫情暴发后才能制备使用；疫苗制备过程不规范，易使病毒灭活不完全，导致猪场长期带毒，加大病毒传播风险等。
4. 结论
本研究利用PK-15细胞从疑似病料中成功分离获得PEDV SNJ-P株，该毒株与传统疫苗株和S non-INDEI型疫苗株亲缘关系均较远，说明当前PEDV毒株变异速度快，商品化疫苗无法提供有效保护。
-
![]()
图 1 密集烤房常规烘烤(左)和明火熏烤(右)工艺
Figure 1. Process of conventional curing in dense roasting room (left) and open fire-curing (right)
![]()
图 2 正(上)、负(下)离子模式下QC样品的总离子流
Figure 2. Total ion chromatography of QC samples under positive (above) and negative (below) ion mode
![]()
图 4 相同温度不同调制方式处理烟叶的差异代谢物KEGG分类图
Figure 4. KEGG taxonomic diagram of differential metabolites in tobacco leaves by different curing methods under the same temperature
![]()
图 5 相同调制方式不同温度处理烟叶的差异代谢物KEGG分类图
Figure 5. KEGG taxonomic diagram of differential metabolites in tobacco leaves under different temperatures by the same curing method
![]()
图 6 差异代谢物火山图
注：FC. 差异倍数。
Figure 6. Volcanic map of differential metabolites
Note: FC. fold change.
表 1 代谢组学液相色谱洗脱梯度
Table 1 Metabolomics liquid chromatography gradient elution
时间/min
time流速/(mL·min−1)
velocity of flow流动相A/%
mobile phase A流动相B/%
mobile phase B0.00 0.35 95 5 9.00 0.35 5 95 10.00 0.35 5 95 11.00 0.35 95 5 14.00 0.35 95 5 
下载: 导出CSV
表 2 质控样本(QC)的相关性
Table 2 Correlation of quality control (QC) samples
QC01 QC02 QC03 QC01 1.00 QC02 0.98 1.00 QC03 0.99 0.98 1.00 注：QC1~QC3分别表示3个质控样本。
Note: QC1-QC3 represent the three quality control samples respectively.
下载: 导出CSV
表 3 差异代谢物筛选结果
Table 3 Screening results of differential metabolites
比较组
comparison
group代谢物总数
total number of
metabolites差异显著代谢物数量
number of
differential metabolites上调代谢物数量
number of
up-regulated metabolites下调代谢物数量
number of
down-regulated metabolitesX-38 vs CK-38 1656 434 386 48 X-42 vs CK-42 1657 477 340 137 CK-42 vs CK-38 1638 309 275 34 X-42 vs X-38 1657 316 260 56
下载: 导出CSV
-
[1] 刘国顺. 烟草栽培学[M]. 北京: 中国农业出版社, 2016. [2] SIFOLA M I, CARRINO L, COZZOLINO E, et al. A survey of fertility program responses of kentucky dark fire-cured tobacco (Nicotiana tabacum L.) yield and quality for cigars manufacture in the Benevento Province (southern Italy)[J]. Beiträ ge zur Tabakforschung International, 2018, 28(1): 14. DOI: 10.2478/cttr-2018-0003.
[3] GIACOMO M D, PAOLINO M, SILVESTRO D, et al. Microbial community structure and dynamics of dark fire-cured tobacco fermentation[J]. Applied and Environmental Microbiology, 2007, 73(3): 825. DOI: 10.1128/AEM.02378-06.
[4] TANTON T W, ABINGTON J B. The effect of different cultural practices upon the major diseases of dark fire-cured tobacco in Malawi[J]. PANS, 1979, 25(2): 158. DOI: 10.1080/09670877909411691.
[5] ABINGTON J B, TANTON T W. Effects of variety, fertilizer and environment on fire-cured tobacco in Malawi[J]. Experimental Agriculture, 1979, 15(2): 153. DOI: 10.1017/S0014479700000545.
[6] TANTON T W, ABINGTON J B. Factors affecting yields of dark fire-cured tobacco in Malawi[J]. Experimental Agriculture, 1979, 15(4): 353. DOI: 10.1017/S00 14479700012990.
[7] 张兴伟, 邢丽敏, 齐义良, 等. 新型烟草制品未来发展探讨[J]. 中国烟草科学, 2015, 36(4): 110. DOI: 10.13496/j.issn.1007-5119.2015.04.020. [8] 蔡洁云, 王惠平, 李雪梅, 等. 无烟气烟草制品现状及发展趋势[J]. 云南民族大学学报(自然科学版), 2019, 28(2): 121. DOI: 10.3969/j.issn.1672-8513.2019.02.005. [9] LEFFINGWELL J C, ALFORD E D, LEFFINGWELL D. Aroma constituents of a supercritical CO2 extract of kentucky dark fire-cured tobacco[J]. Leffingwell Reports, 2013, 5(1): 1. DOI: 10.13140/RG.2.2.12294.37440.
[10] FIEHN O, KOPKA J, DÖRMANN P, et al. Metabolite profiling for plant functional genomics[J]. Nature Biotechnology, 2000, 18(11): 1157. DOI: 10.1038/81137.
[11] HARDY N W, HALL R D. Plant metabolomics: methods and protocols[M]. Totowa: Humana Press, 2012.
[12] 顾渝娟, 吴振先. 代谢组学在植物研究中的应用[J]. 广东农业科学, 2012, 39(4): 105. DOI: 10.16768/j.issn.10 04-874x.2012.04.003. [13] 彭翠云. 延迟采收对豫中烤烟上部叶生理指标和代谢组学的影响[D]. 郑州: 河南农业大学, 2021. [14] 田培, 李晓辉, 贺文俊, 等. 南方根结线虫侵染烟草抗感品种前后的代谢组学分析[J]. 中国烟草学报, 2019, 25(4): 81. DOI: 10.16472/j.chinatobacco.2018.261. [15] 李勇, 林茜, 逄涛, 等. 基于气相色谱—质谱的烟草根部代谢组学方法建立及应用[J]. 分析试验室, 2017, 36(8): 910. DOI: 10.13595/j.cnki.issn1000-0720.2017.0195. [16] 邵琪, 解群杰, 刘贺娟, 等. 接种枯萎病菌对甜瓜脂氧合酶活性及基因表达的影响[J]. 植物生理学报, 2016, 52(4): 479. DOI: 10.13592/j.cnki.ppj.2015.0473. [17] 马俊卿. 基于多组学方法研究AMF共生诱导的玉米根系分泌物变化及变化机理[D]. 南宁: 广西大学, 2022. [18] 张坤, 孙仕仙, 石傲傲, 等. 扑草净胁迫下香根草叶片的代谢组学研究[J]. 西南林业大学学报(自然科学), 2022, 42(3): 90. DOI: 10.11929/j.swfu.202108076. [19] 徐忠山. 利用多组学联合分析燕麦响应干旱胁迫及盐胁迫的分子机制[D]. 呼和浩特: 内蒙古农业大学, 2021. [20] 杜宪. 高温、CO2耦合作用对结果期温室黄瓜生理生化及代谢组学影响的研究[D]. 呼和浩特: 内蒙古农业大学, 2022. [21] LI Z Q, ZHANG Y, LI H, et al. Genome-wide characterization and expression analysis provide basis to the biological function of cotton FBA genes[J]. Frontiers in Plant Science, 2021, 12: 696698. DOI: 10.3389/FPLS.2021.696698.
[22] 张蕾, 阮羽萱, 潘娇, 等. 甘蓝型油菜响应干旱胁迫的转录组分析[J]. 江苏农业科学, 2023, 51(5): 51. DOI: 10.15889/j.issn.1002-1302.2023.05.007. [23] 王弛, 刘慧, 刘芳, 等. 基于代谢组学探究金线莲对铅胁迫的响应机制[J]. 南方农业学报, 2022, 53(12): 3478. DOI: 10.3969/j.issn.2095-1191.2022.12.019. [24] 甄熙. 野生罗布麻生态适应性及响应盐胁迫机制研究[D]. 呼和浩特: 内蒙古农业大学, 2022. [25] 赵景娜. 植物乳植杆菌P9对甲拌磷胁迫的响应机制研究[D]. 呼和浩特: 内蒙古农业大学, 2022. [26] 刘苏宁, 王力荣, 方伟超, 等. 基于广泛靶向代谢组学的桃芽抗寒代谢物的筛选与鉴定[J]. 果树学报, 2023, 40(1): 1. DOI: 10.13925/j.cnki.gsxb.20220371. [27] 罗素雅, 虞新磊, 陈琼琳, 等. 利用代谢组学揭示脱落酸对高温胁迫龙须菜(Gracilariopsis lemaneiformis)的作用机制[J]. 海洋与湖沼, 2022, 53(1): 224. DOI: 10.116 93/hyhz20210800173. [28] HAN R C, HE X F, PAN X H, et al. Enhancing xanthine dehydrogenase activity is an effective way to delay leaf senescence and increase rice yield[J]. Rice, 2020, 13(1): 16. DOI: 10.1186/s12284-020-00375-7.
[29] KAFER C, ZHOU L, SANTOSO D, et al. Regulation of pyrimidine metabolism in plants[J]. Frontiers in Bioscience, 2004, 9(2): 1611. DOI: 10.2741/1349.
-
期刊类型引用(1)
1. 尤永君，刘思桐，王幸，潘静，邱贞娜. 猪流行性腹泻病毒流行病学调查及毒株致病力研究. 天津科技. 2023(S1): 81-87 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 963
- PDF下载量: 38
- 被引次数: 1
