水氮互作對陜南烤煙煙堿合成及代謝的影響

Effects of Water and Nitrogen Interaction on Nicotine Synthesis and Metabolism of Flue-Cured Tobacco in Southern Shaanxi

HAO Sujuan'，ZHANG Aojie2， WANG Pingping^，ZHANG Mingcheng2，AI Suilong2，ZHANG Lixin2 （'Lianhu DistrictMeteorologicalServiceCenter，Xi'an7o;CollegeofLifeScience，NorhwestAamp;FUniversity， Xianyang712ioo，Shaanxi;YantaDistrict InstituteofTobaccoScience，Xi'an710000）

煙草是我國的主要葉用經(jīng)濟作物之一，是卷煙生產(chǎn)的主要原料，還可用于提取尼古丁、木質(zhì)素等，具有較高的商業(yè)價值[1]。此外，煙草在國家稅收中有著舉足輕重的地位，同時在精準扶貧方面發(fā)揮著重要作用?？緹熒a(chǎn)對提高農(nóng)民收人、推動地方經(jīng)濟發(fā)展有重要意義。

水分和肥料是影響煙草生長的最為重要的兩大栽培因素。水分是肥料發(fā)揮效用、被作物利用的基礎，肥料是水土系統(tǒng)中影響作物生產(chǎn)效率的關鍵[3]。氮是作物生產(chǎn)的主要驅(qū)動因素之一，與土壤水分條件有明顯互作[4-5]。施氮和灌水在提高烤煙水分利用率上存在顯著的互作效應，在土壤水分不足時，施加氮肥能夠促進烤煙對水分的利用，而在施氮量較多時，灌水又能促進烤煙對肥料的吸收利用。水氮之間的互作效應顯著促進了烤煙干物質(zhì)的積累，影響烤煙株高、莖圍、葉面積等農(nóng)藝性狀，從而提高煙葉產(chǎn)量。 ¹⁵N 穩(wěn)定同位素試驗表明，水氮互作可以顯著提高烤煙的氮肥利用率[。合理的水氮組合能協(xié)調(diào)烤煙各種化學成分含量，提高煙葉品質(zhì)。此外，施氮量增加可以提升烤煙的氧化防御系統(tǒng)，減輕水分脅迫下烤煙的氧化損傷[7]。

代謝物是細胞調(diào)控終產(chǎn)物，參與多種生化反應及生理調(diào)控[，最接近生物的表型或者生理狀態(tài)，其種類、含量變化可以看作是生物體對表型或基因變化響應的結(jié)果，更能反映生物的基因、環(huán)境與表型三者之間的關系。代謝組學應運而生，主要通過高通量檢測對生物體在限定條件下的代謝物進行定性、定量分析，探究其代謝產(chǎn)物差異及代謝途徑的變化規(guī)律來揭示內(nèi)在機制[9-10]。代謝組學作為系統(tǒng)生物學新的分支，目前已經(jīng)成為繼基因組學、轉(zhuǎn)錄組學和蛋白質(zhì)組學之后新的研究熱點。Zhao等[運用代謝組學方法分析不同產(chǎn)地來源的烤煙代謝物差異，發(fā)現(xiàn)不同產(chǎn)地的烤煙代謝物呈現(xiàn)明顯的分離，不同生態(tài)區(qū)域可影響烤煙氮代謝等多種代謝途徑，進一步篩選出多種與不同產(chǎn)地煙葉香味特征有關的差異代謝產(chǎn)物，為煙葉香型的研究提供了理論依據(jù)。隨著相關技術的發(fā)展，在玉米、水稻等作物上通過代謝組學方法進行機理探索的研究逐步深入。代謝組學技術在煙草研究中的應用也逐步展開[12-15]，成為全面、系統(tǒng)研究煙草對環(huán)境或外界刺激響應機理的重要手段[12]。

本研究針對陜南煙區(qū)烤煙生產(chǎn)中的水氮供應匹配性，以云煙99為材料，探究水氮互作對烤煙煙堿含量、煙堿合成過程中關鍵酶活性的影響，同時采用代謝組學方法探究不同水氮處理影響煙堿合成的內(nèi)在機制，旨在為烤煙產(chǎn)質(zhì)量調(diào)控提供參考。

1材料與方法

1.1試驗地概況試驗于2021年4-8月在西北農(nóng)林科技大學北校區(qū)曹新莊試驗站（ 34^°17^′23^′′N 108^°4^′12^′′E ，海拔 520m ）日光大棚內(nèi)完成。種植土壤為陜西省商洛市洛南縣植煙黃棕壤土，基本理化性質(zhì)： pH 值7.22，有機質(zhì)含量 10.25g/kg ，速效氮含量 30.25mg/kg ，速效磷含量 34.25mg/kg ，速效鉀含量 127.48mg/kg 。

1.2試驗設計烤煙品種為云煙99，由陜西省煙草科學研究所提供。本試驗涉及水、氮2個因素，其中土壤水分條件3個水平：W1（土壤相對含水量75%～85% ）W2（土壤相對含水量 55%～65% ）、W3（土壤相對含水量 35%～45% ），施氮量3個水平：N1（ 4g/ 株）N2（ 5g/ 株）N3（ 6g/ 株），雙因素完全組合設計共計9個處理（表1。煙苗成苗后移入桶中，每桶定植1株，移栽前保持桶中土壤相對含水量60% ，移栽后36d進人旺長期，開始進行水分處理。采用稱重法補水，控制土壤相對含水量，控水處理進行 30d 。按照試驗處理，稱取相應克數(shù)的氮肥給每株煙苗進行基施。

1.3樣品采集于移栽后36d（旺長初期）51d（旺長中期）66d（旺長末期）85d（成熟期）分別對根部及葉片進行取樣。葉片于 105% 殺青 30min 65°C 烘干至恒重，研磨并過100目篩，用于煙堿含量測定。取根尖3cm迅速放入液氮，后保存于 -80% 冰箱，用于測定煙堿合成途徑中關鍵酶的活性、代謝產(chǎn)物含量。取樣方式為各處理隨機3株樣品混合，每個處理3次重復。

1.4測定項目與方法煙堿含量測定參考鄭璞帆[的方法并稍作改動，采用HCI浸提法。稱取粉碎過篩烘干的煙葉樣品 0.03g ，加人 0.5mol/L 鹽酸溶液 25mL ，活性炭1/3角匙，震蕩 10min 后過濾，定容至 100mL ，用紫外分光光度計測定 236nm 、259nm 和 282nm 波長的吸光度。

煙堿合成途徑中關鍵酶活性測定煙草根樣添加液氮后進行粉碎，取 0.1g 加人 0.05mol/L 磷酸緩沖液（pH值7.4） 1.5mL ，充分混勻后 4^°C 下8000r/min 離心 20min ，上清液用于精氨酸脫羧酶（ADC）鳥氨酸脫羧酶（ODC）腐胺 N_- 甲基轉(zhuǎn)移酶（PMT）的測定，其測定步驟按照上海通蔚生物科技有限公司的酶聯(lián)免疫試劑盒說明書進行。

代謝產(chǎn)物含量測定樣品前處理：煙草根樣添加液氮后進行粉碎，取根樣 0.2g ，加入 75% 甲醇溶液 400μL ，充分混勻后超聲提取 30min，4^°C 下 14000r/min 離心 20min ，上清液用 0.22μm 有機相過濾器過濾后加入 2mL 進樣瓶中，用于高分辨離子淌度液質(zhì)聯(lián)用儀HRLCMS（ LC-30A+ TripleTOF5 5600+ ）上機進行正、負離子模式分析。

數(shù)據(jù)庫匹配將提取的峰信息在MassBank、Respect、GNPS這3個庫進行全庫檢索與比對。數(shù)據(jù)庫匹配參數(shù)見表2。

1.5數(shù)據(jù)統(tǒng)計與分析采用Excel2013、IBMSPSSStatistics26和RStudio1.3.1093.0對數(shù)據(jù)進行整理、統(tǒng)計分析和做圖。

2 結(jié)果與分析

2.1水氮互作對烤煙煙堿含量的影響由表3可以看出，不同水氮處理的烤煙中部葉煙堿含量在不同時間段均逐漸上升。移栽后36d，各處理的中部葉煙堿含量差異不大，此后水分控制開始，各處理的中部葉煙堿含量差異逐漸增大，在移栽后66d各處理中部葉煙堿含量差異最為顯著。整體來看，同一時期相同氮肥處理下，中部葉煙堿含量隨著土壤相對含水量的增加而降低；W1N1和W1N2處理在整個過程中保持相對平穩(wěn)的增長，在移栽后85d煙堿含量均較低；W3水分條件下，各時期中部葉煙堿含量均較高。同一時期相同水分處理下，中部葉煙堿含量隨著氮肥的增加而增加。

2.2水氮互作對煙堿合成途徑關鍵酶活性的影響

由圖1可知，移栽后36d即控水之前各處理的精氨酸脫羧酶（ADC）活性無顯著差異。移栽后51d，相同氮肥處理下ADC活性呈現(xiàn)隨土壤相對含水量降低而升高的趨勢，且W3處理下ADC活性達到W1條件下的 1.64～1.81 倍；W1水分條件下，W1N2處理的ADC活性略低于其余處理，但無顯著差異；W2、W3水分條件下，ADC活性隨施氮量增加而增加，W2N2、W2N3處理的ADC活性較W2N1處理分別提高 10.7%.21.1% ;W3N2、W3N3處理的ADC活性較W3N1處理分別提高 0.5%.0.8% ；因此，隨土壤相對含水量降低施氮量對ADC活性的提升效果減弱。移栽后66d，W2水分條件下ADC活性最高，且隨施氮量增長ADC活性提高；W3水分條件下，ADC活性較移栽后51d下降，且隨施氮量增加呈現(xiàn)先增加后降低的趨勢。移栽后85d，所有處理的ADC活性均處于較高水平，施氮量對ADC活性影響效果減弱。

由圖2可知，W1、W2水分條件下，各處理的鳥氨酸脫羧酶（ODC）活性隨移栽后時間增加而增加，W3水分條件下，各處理ODC活性呈現(xiàn)隨移栽后時間增加先增加后降低再增加的趨勢。移栽后36d，各處理的ODC活性無顯著差異；移栽后51d，W2、W3條件下ODC活性隨施氮量增加而升高，其中W2條件下增幅較高，W2N2、W2N3較W2N1分別提高了 17.6%.28.6% ;移栽后66d，W3水分條件下，各處理較移栽后51d分別降低了 2.1%.6.9%.2.5% 移栽后85d，W1水分條件下ODC活性仍隨施氮量增加而提高，W2、W3水分條件下則隨施氮量增加呈現(xiàn)先增加后降低的趨勢。

由圖3可知，所有處理的腐胺N-甲基轉(zhuǎn)移酶（PMT）活性均隨著移栽時間的增加而增加，移栽后36d各處理之間的PMT活性無顯著差異；進行水分控制后，各處理PMT活性出現(xiàn)差異，整體表現(xiàn)為隨土壤相對含水量的降低PMT活性升高。移栽后36＼～66d，均以W3N3處理的PMT活性最高；移栽后85d，W1、W3水分條件下不同施氮量處理的PMT活性相近；W2水分條件下，W2N3處理的PMT活性最高，比W2N1提高了 17.6% 。

2.3水氮互作下烤煙代謝組學分析

2.3.1總體代謝物層次聚類分析對煙堿含量差距較大的W3N3和W1N1兩組處理進行代謝組學分析，由PCA主成分分析結(jié)果可知進樣過程中儀器運行較為穩(wěn)定，可以用于后續(xù)的分析。對鑒定出的全部代謝物的相對定量值歸一轉(zhuǎn)化后進行層次聚類得到圖4，列為所有樣本，行為全部代謝物，紅色代表高表達代謝物，藍色代表低表達代謝物。由圖4可以看出，有部分代謝物高、低表達交互在一起，但高、低表達區(qū)域較明確。同一組別樣品的高、低表達區(qū)域較為相近，不同組別樣品間高、低表達區(qū)域差異明確。

產(chǎn)物中包括10種有機酸，7種胺類，5種醇類，3種酚類，2種酮類，2種酯類，1種氨基酸類，1種堿類，其他類8種。

2.3.3代謝通路富集分析如圖6所示，差異代謝產(chǎn)物主要富集的代謝通路為：氨酰tRNA生物合成，ABC轉(zhuǎn)運器，氨基酸生物合成，嘌呤代謝，代謝途徑，咖啡因代謝，纈氨酸、亮氨酸和異亮氨酸生物合成，丙氨酸、天冬氨酸、谷氨酸代謝，谷胱甘肽代謝，單酰胺類生物合成代射通路。氣泡圖中，氣泡顏色越深代表該通路差異顯著性越高，直徑越大代表富集到的代謝物越多。因此，在富集到的通路中，氨酰tRNA生物合成途徑富集的代謝物較少，但差異顯著性較高；代謝途徑富集的代謝物最多。

2.3.4煙堿代謝相關產(chǎn)物將已識別的代謝產(chǎn)物中與煙堿合成相關的代謝產(chǎn)物進行聚類分析，結(jié)果如圖7，有部分代謝物高、低表達交互在一起，但有相對明確的高、低表達區(qū)域。

a為柱狀圖，b為氣泡圖；Valine，leucineandisoleucinebiosynthesis為氨酸亮氨酸和異亮氨酸生物合成，Purinemetabolism為嘌呤代謝，Monobactam biosynthesis 為單酰胺類生物合成代謝通路，Metabolic pathways為代謝途徑，Glutathione metabolism為谷胱甘肽代謝，

Cafeine metabolism為咖啡因代謝，Biosynthesis ofaminoacids為氨基酸生物合成，Aminoacy1-tRNAbiosynthesis為氨酰tRNA生物合成，Alanine，aspartate and glutamate metabolism為丙氨酸、天冬氨酸和谷氨酸代謝，ABC transporter為ABC轉(zhuǎn)運器

對煙堿合成相關代謝產(chǎn)物含量進行相關性分析，結(jié)果如圖8。有9對代謝產(chǎn)物呈顯著正相關，有9對代謝產(chǎn)物呈顯著負相關。天冬氨酸、精氨酸含量與煙堿含量呈顯著正相關。初步判斷水氮處理對烤煙煙堿含量造成的影響可能是由于天冬氨酸、精氨酸含量的變化。

3討論

常規(guī)化學成分是判斷煙葉品質(zhì)好壞的關鍵因素[17]，煙堿是煙葉最重要的化學成分之一，煙堿是煙草特有的生物堿，影響烤煙的勁頭和刺激性，其含量高低是衡量烤煙品質(zhì)的重要指標[18-22]。土壤水分條件會影響煙堿含量，本研究中中部葉煙堿含量表現(xiàn)為土壤相對含水量 35%～45% 最高， 75%～85%最低，這與汪耀富等[3研究結(jié)果一致。氮作為煙堿的重要組成成分，其施用量與煙堿含量密切相關。水分和氮肥在作物生長中存在相互促進、相互制約的關系。馬振奇[23]進行水氮耦合研究發(fā)現(xiàn)，不同水氮處理下中部葉的煙堿含量均呈現(xiàn)隨土壤相對含水量降低而升高，隨施氮量增加而升高的規(guī)律。本研究中中部葉煙堿含量整體變化規(guī)律與之一致。

張瑞霞[24]測定了不同滲透脅迫程度下24h內(nèi)烤煙葉片精氨酸脫羧酶和鳥氨酸脫羧酶的酶活性變化，發(fā)現(xiàn)精氨酸脫羧酶活性隨脅迫程度的加劇而顯著升高，這與本研究一致；不同的是，在該研究中隨脅迫程度加劇、時間延長，烤煙鳥氨酸脫羧酶活性無顯著變化，而在本研究中移栽后36d、51d烤煙根部鳥氨酸脫羧酶活性受土壤相對含水量影響較小，在移栽后66d土壤相對含水量 35%～45% 條件下其活性降低，這種差異是否是由取樣部位不同、處理時長不同或者是由植株生長水平不同帶來的尚需進一步研究。本研究中，不同土壤水分條件下施氮量的增加不同程度上提高了精氨酸脫羧酶、鳥氨酸脫羧酶活性，但其作用機制尚不明確。黃化剛等[25]研究表明，同等用氮量下有機氮相比于無機氮可以提高鳥氨酸脫羧酶的表達量，導致烤煙上部葉煙堿含量增加。不同氮素形態(tài)也會影響精氨酸脫羧酶、鳥氨酸脫羧酶活性，從而導致煙堿積累量的差異[26]?？緹煾蹈種-甲基轉(zhuǎn)移酶活性與葉片煙堿含量呈顯著正相關[27-28]。

近年來，代謝組學不斷發(fā)展，為煙草的提質(zhì)增香及代謝機理研究奠定了基礎。本研究中，精氨酸含量與煙堿含量呈正相關。對比分析發(fā)現(xiàn)，這種差異的形成可能是由于取樣部位不同造成，本研究所用材料為合成煙堿的根系部位，而向海英等[29]研究對象為烤煙葉片部位，根部與葉片部位的這種差異與煙堿合成、轉(zhuǎn)運之間的關系仍需進一步研究；此外，研究共識別出4086個含有二級碎片的一級質(zhì)譜峰，鑒別出245個差異代謝產(chǎn)物。本研究中，對不同水氮處理后煙堿含量差異較高的兩個處理進行代謝物提取及組學分析，共鑒別出190種代謝物，這種差異可能是由于取樣部位、提取效率、儀器分辨率等不同導致。識別出的代謝產(chǎn)物數(shù)量及精準度是代謝組學精確分析的基礎，完善植物組織中代謝產(chǎn)物提取方法，提高代謝產(chǎn)物提取效率及比對精準度是代謝組學發(fā)展的重要內(nèi)容。

4結(jié)論

施氮量對煙堿合成的作用效果因土壤相對含水量不同而有所差異，中部葉煙堿含量在土壤相對含水量 35%～45% 條件下最高，且以施氮量 6g/ 株處理煙堿含量最高，煙堿含量以W1N1處理最低。不同水分條件下適宜的施氮量可以提高煙堿合成關鍵酶活性，代謝物差異分析結(jié)果表明，差異代謝物在氨酰tRNA生物合成途徑、ABC轉(zhuǎn)運器通路、氨基酸生物合成途徑及嘌呤代謝中顯著富集；煙堿合成途徑代謝產(chǎn)物存在顯著差異，其中天冬氨酸、精氨酸含量與煙堿含量呈顯著正相關。

參考文獻

[1]陳雪嬌，李云飛，李剛，宗凱，鄭海松，孫娟娟，姚劍，余曉峰，徐業(yè)平.煙草主要病害及其相關檢測方法研究進展.河南農(nóng)業(yè)科學，2013，42（11）：1-6

[2]周繼洲.基于GIS的陜西省環(huán)秦嶺生態(tài)煙區(qū)煙草生態(tài)適宜性區(qū)劃.楊凌：西北農(nóng)林科技大學，2014

[3]汪耀富，楊天旭，孫德梅，葉紅潮.灌水及氮素形態(tài)對烤煙營養(yǎng)元素含量及產(chǎn)量品質(zhì)的影響.河南農(nóng)業(yè)大學學報，2006，40（5）：477-481

[4]ZhaoSM，ZhengH，LuX，LiuY，SuBL，XuGW.Metabonomicsand phospholipid metabolic profilingof abnormal glucose metabolismbasedon high performance liquidchromatography-electrospraymassspectrometry.Chinese Journal of Chromatography，2011，29（4）：307-313

[5]LiQH，ZhaoCX，LiY，ChangYW，WuZM，PangT，LuX，WuY，XuGW.Liquid chromatography/mass spectrometry-based metabolicprofiling to elucidate chemical differences of tobacco leaves betweenZimbabwe and China.Journal of Separation Science，2010，34（2）：119-126

[6]袁仕豪.水氮互作對烤煙氮素吸收利用及煙葉產(chǎn)量和品質(zhì)的影響.鄭州：河南農(nóng)業(yè)大學，2008

[7]崔保偉，陸引罡，張振中，任琨.不同生育期水分脅迫對烤煙生理特性及化學品質(zhì)的影響.中國煙草科學，2009，30（3）：19-23

[8]高婭北，李成剛，陳二龍，張明剛，孔德輝，孫占偉，史龍飛，宋朝鵬.烤煙煙葉類黃酮合成關鍵酶活性及其調(diào)控基因的表達分析.中國煙草學報，2019，25（1）：86-92

[9]謝良文，路曉崇，裴曉東，蔣博文，李生棟，魏碩，李帆.烤煙烘烤過程中氨基酸含量變化對主要含氮化合物的影響.湖南農(nóng)業(yè)大學學報：自然科學版，2016，42（2）：152-156

[10]PiornosJA，BalagiannisDP，Methven L，KoussissiE，BrouwerE，ParkerJK.Elucidating the odor-activearoma compoundsinalcohol-freebeerand theircontributiontotheworty flavor.JournalofAgricultural andFoodChemistry，2020，68（37）：10088-10096

[11]ZhaoJY，HuCX，ZengJ，ZhaoYN，ZhangJJ，ChangYW，LiLL，ZhaoCX，LuX，XuGW.Studyofpolarmetabolitesintobacco from different geographical originsby usingcapillaryelectrophoresis-mass spectrometry.Metabolomics，2014，10（5）：805-815

[12]王小莉，付博，趙銘欽，賀凡，王鵬澤，劉鵬飛.代謝組學技術在煙草研究中的應用進展.中國煙草科學，2016，37（1）：89-96

[13]梁珂馨，奚家勤，陳若星，牟文君，劉萍萍，鄭慶霞，王得強，賈中林，李曉闖，張艷玲.基于非靶向代謝組學的片煙醇化進程判斷及代謝物差異分析.煙草科技，2023，56（4）：17-25

[14]虞昕磊，何結(jié)望，林國平，李金海，王大愛，袁躍斌，劉圣高，李志豪，陶德欣.夏冬兩季發(fā)酵雪茄煙葉的代謝組差異分析.生物技術通報，2024，40（6）：260-270

[15]陸岸芷，張樂毅，王平平，陳新穎，杜姍姍，李先鋒，王輝，張強，艾綏龍，張立新.半干旱區(qū)2種雪茄煙葉發(fā)酵后差異代謝物分析.干旱地區(qū)農(nóng)業(yè)研究，2024，42（4）：22-28

[16]鄭璞帆.陜南地區(qū)不同成熟度煙葉生理特性和烤后品質(zhì)研究.楊凌：西北農(nóng)林科技大學，2017

[17]康崢，唐衛(wèi)迪，楊繼利，曾皓，趙曉軍，李俊花，孔云鑫，張廷金，鄧綱，張海平.不同生物菌劑對根結(jié)線蟲病害地區(qū)煙草生長及品質(zhì)的影響.中國種業(yè)，2024（8）：79-86

[18]鄧金華.三種有機物料配施對土壤水分特性及烤煙生理與品質(zhì)的影響.楊凌：西北農(nóng)林科技大學，2022

[19]雷亞芳.不同酶制劑與香料處理對煙葉品質(zhì)及微生物多樣性的影響.楊凌：西北農(nóng)林科技大學，2021

[20]向羽，周成，黃楚天，何輝，柳洪入.結(jié)構(gòu)優(yōu)化對不同葉位煙葉品質(zhì)及揮發(fā)性物質(zhì)含量的影響.安徽農(nóng)業(yè)科學，2024，52（10）：15-20

[21]吳鑫穎.煙葉醇化過程的動態(tài)解析及功能微生物研究.無錫：江南大學，2022

[22]陳愛國，劉光亮，周道金，謝鳳標，周培祿，蔣代兵，周東新，林天然，沈燦柱，曾文龍.龍巖鮮煙葉特征香氣成分前體物代謝組學的通路分析.安徽農(nóng)業(yè)大學學報，2018，45（3）：545-550

[23]馬振奇.水氮耦合對烤煙生長發(fā)育以及土壤養(yǎng)分的影響.鄭州：河南農(nóng)業(yè)大學，2016

[24]張瑞霞.烤煙激素與多胺對水分迫和氮肥形態(tài)的響應及其調(diào)節(jié)研究.鄭州：河南農(nóng)業(yè)大學，2004

[25]黃化剛，申燕，吳飛躍，楊振智，劉世碧，喻奇?zhèn)?，賈宏昉.不同氮素來源對烤煙上部葉碳氮代謝及其基因表達的影響.南方農(nóng)業(yè)學報，2018，49（3）：462-468

[26]王學龍.礦質(zhì)營養(yǎng)對烤煙煙堿的調(diào)控及其生理機制研究.武漢：華中農(nóng)業(yè)大學，2012

[27]蒲媛媛.煙堿合成途徑中關鍵酶活性對雜種優(yōu)勢形成的影響.貴陽：貴州大學，2020

[28]MizusakiS，TanabeY，NoguchiM，TamakiE.Changesintheactivitiesof ornithine decarboxylase，putrescine N-methyltransferaseandN-methylputrescine oxidase in tobaccorootsin relationtonicotinebiosynthesis.PlantandCellPhysiology，1973，14（1）：103-110

[29]向海英，劉欣，逢濤，曾婉俐，許力，高茜，鄧樂樂，張建鐸，李雪梅，李勇，孔維松.基于非靶向代謝組學的紅花大金元不同生育期煙葉代謝圖譜差異分析.中國煙草學報，2022，28（4）：77-84

（收稿日期：2025-02-19）