青花椒種質資源重金屬鎘積累的差異分析

摘要：【目的】分析青花椒種質資源鎘積累的差異，探討鎘積累特性及富集機制，為青花椒品種選育和鎘防控提供理論依據。【方法】以重慶市江津區花椒國家林木種質資源庫中生產性狀較好的九葉青花椒、早熟九葉青花椒、關嶺藤椒、四川藤椒、云南藤椒和云南大麻椒6個種質為試驗材料，測定葉片、枝干、枝皮、果皮、種子和根系6個部位的鎘含量，分析鎘富集與轉運特征；選取鎘富集能力差異較大的種質，對其根系代謝物進行非靶向代謝組學分析。【結果】不同青花椒種質同一部位的鎘含量存在差異，其中云南藤椒各部位鎘含量均較高于其他種質。同一青花椒種質各部位鎘含量也存在差異，除關嶺藤椒外，其他種質各部位鎘含量均以枝皮為最高，其次是葉片、根系，且均高于種子、枝干、果皮。各青花椒種質根系對土壤鎘均有富集能力，富集系數為3.470～12.780，其中云南藤椒的富集能力顯著強于其他種質（Plt;0.05，下同），早熟九葉青花椒顯著弱于其他種質。根系到果皮、根系到葉片、根系到枝干、根系到枝皮、根系到種子的轉運系數分別為0.264～0.530、1.494～3.294、0.324～0.642、1.408～3.628、0.330～0.943。對鎘富集能力差異最大的早熟九葉青花椒與云南藤椒根系代謝物進行非靶向代謝組學分析，共鑒定出520個差異代謝物，其中196個代謝物上調表達，324個代謝物下調表達，2個種質間存在明顯差異；與早熟九葉青花椒根系代謝物相比，云南藤椒根系代謝物中有機酸、生物堿呈上調表達趨勢；KEGG代謝通路富集分析結果表明，差異代謝物主要富集在異黃酮生物合成、苯丙氨酸代謝、精氨酸生物合成、牛磺酸和亞牛磺酸代謝、亞油酸代謝等通路。【結論】青花椒種質資源對土壤鎘的積累存在差異，早熟九葉青花椒富集能力較弱，云南藤椒富集能力較強，其富集能力差異可能與根系代謝物有關。

關鍵詞：青花椒；鎘；富集系數；轉運系數；根系代謝物

中圖分類號：S573.9文獻標志碼：A文章編號：2095-1191（2025）02-0343-12

Differential analysis on heavy metal cadmium accumulation of Zanthoxylum armatum germplasm resources

LUO You-jin，CHEN Xia，WU Chun-qing，HU Jia-yu，ZHANG Yi-gang，CHENG Yue-qing，HUANG Ming，XIE Yong-hong*

（Fruit Research Institute，Chongqing Academy of Agricultural Sciences，Chongqing 401329，China）

Abstract：【Objective】To analyze the differences in cadmium（Cd）accumulation among Zanthoxylum armatum germ-plasm resources，explore the characteristics of cadmium accumulation and enrichment mechanisms，which could provide theoretical reference for the breeding of Z.armatum varieties and cadmium control.【Method】Six germplasms with good production traits from the National Forest Tree Germplasm Resource Repository of Zanthoxylum armatum in Jiangjin Dis-trict，Chongqing，including Z.armatum‘Jiuyeqinghuajiao’，Z.armatum‘Zaoshu Jiuyeqinghuajiao’，Guanling Z.arma-tum，Sichuan Z.armatum，Yunnan Z.armatum，Yunnan Zanthoxylum sp.，were used as experimental materials.The Cd content in 6 parts，including leaves，stems，branch bark，pericarp，seeds and roots，was measured to analyze the charac-teristics of Cd enrichment and transport.And germplasms with great differencesin Cd enrichment ability were selected for untargeted metabolomics analysis on root metabolites.【Result】Differencesin Cd content were observed in the same part among different Z.armatum germplasms，with Yunnan Z.armatum showing higher cadmium content in all parts com-pared to other germplasms.Within the same germplasm，variations in Cd content were also found among different parts.Except for Guanling Z.armatum，the branch bark had the highest Cd content in other germplasms，followed by leaves and roots，which were higher than those in seeds，stems and pericarp.All Z.armatum germplasms had soil Cd enrichment ability in their roots，with enrichment coefficients ranging from 3.470 to 12.780.Yunnan Z.armatum showed significantly stronger enrichment ability than other germplasms（Plt;0.05，the same below），while Z.armatum‘Zaoshu Jiuyeqinghua-jiao’was significantly weaker than other germplasms.The transfer coefficients from root to pericarp，root to leaf，root to stem，root to branch bark，and root to seed were 0.264-0.530，1.494-3.294，0.324-0.642，1.408-3.628 and 0.330-0.943 respectively.Untargeted metabolomics analysis was conducted on the root metabolites of Z.armatum‘Zaoshu Jiuyeqing-huajiao’and Yunnan Z.armatum，which had the greatest difference in Cd enrichment ability.A total of 520 differential metabolites were identified，including 196 up-regulated and 324 down-regulated metabolites，with significant differences between the 2 germplasms.Compared with the root metabolites of Z.armatum‘Zaoshu Jiuyeqinghuajiao’，the root me-tabolites of Yunnan Z.armatum showed an up-regulated trend in organic acids and alkaloids.The results of KEGG meta-bolic enrichment analysis indicated that the differential metabolites were mainly enriched in pathways such as isoflavonoid biosynthesis，phenylalanine metabolism，arginine biosynthesis，taurine and hypotaurine metabolism，and linoleic acid metabolism.【Conclusion】There are differences in the accumulation of soil Cd by Z.armatum germplasm resources.Specifically，Z.armatum‘Zaoshu Jiuyeqinghuajiao’demonstrated relatively low accumulation ability，whereas Yunnan Z.armatum shows high enrichment ability.This variation in Cd enrichment might be associated with root metabolites.

Key words：Zanthoxylum armatum；cadmium；enrichment coefficient；transfer coefficient；root metabolites

Foundation items：National Key Research and Development Program of China（2021YFD1000405）；Chongqing Academy of Agricultural Sciences Municipal Financial Youth Special Fund Project（cqaas2023sjczqn002）；Chongqing Agricultural Development Fund Project（NKY-2021AC001）

0引言

【研究意義】花椒作為香辛料在食品加工及家庭烹飪中被廣泛使用，按其顏色可分為青花椒和紅花椒，具有賦香、掩蓋異味、防腐、著色、保健、增加食欲等作用，且具有一定藥用價值，有溫中散寒、止痛、殺蟲等功效（趙秀玲，2012；國家藥典委員會，2015；席少陽等，2021）。目前，我國是世界花椒生產、消費和出口第一大國；重慶是青花椒的發源地，其青花椒種植規模達7萬ha，年產值超50億元。根據2014年《全國土壤污染狀況調查公報》，我國耕地土壤質量不容樂觀，耕地土壤污染點位超標率為19.4%，無機污染物中鎘的點位超標率最高，且主要分布在西南地區（劉霈珈，2017）。花椒作為典型的鎘富集植物，在重慶土壤重金屬背景值較高的情況下，青花椒鎘含量通常為0.05～0.45 mg/kg，超標較為普遍，已成為制約其產業可持續發展的主要因素。因此，探究青花椒種質資源的鎘富集特性及其機理對青花椒產業可持續發展具有重要意義。【前人研究進展】目前，已有較多關于花椒種質資源、有效成分及綜合利用的研究。杜凌等（2018）對不同品種花椒的果皮微量元素進行分析，結果發現銅、錳、鉀3個元素呈正態分布，可分為4級，將花椒樣品區分開來；敖厚豫等（2020）研究發現，通過頂壇花椒麻味成分的指紋圖譜可將其與其他花椒品種進行區分；馬堯（2021）對72份花椒種質資源的品種類型進行鑒定，結果發現不同花椒品種類型的品質性狀存在差異；王林等（2022）對南方4個主產區的紅花椒主要揮發性物質進行分析，結果發現其氣味差異的主要來源是非共有物質；尹慶豐等（2023）對36份花椒種質資源進行簡化基因組測序，結果發現其存在較大程度的遺傳分化，能為優良品種選育提供材料選擇。狄科（2011）對花椒多酚的最佳浸提條件進行研究，并證實花椒多酚具有較強的抗氧化性；徐雙美等（2020）對四川不同產地竹葉花椒的揮發性成分進行測定，結果發現22個種植區的竹葉花椒樣品成分相似但含量存在差異；慕鈺文等（2021）通過研究大紅袍花椒籽不同溶劑提取物的活性成分差異，發現不同溶劑對花椒籽的總多酚、總黃酮含量及DPPH自由基清除率和鐵還原能力影響較大。于霞（2011）研究發現，以花椒籽廢渣制備的改性活性炭對硝基苯酚、鉛有較強的吸附作用；王紀磊等（2020）研究發現，施用花椒籽（餅）能有效防治番茄根結線蟲。此外，也有部分學者對花椒鎘含量分布情況進行研究。余曉琴（2010）對36批花椒樣品進行重金屬含量測定，結果表明鎘含量為0.07～0.52 mg/kg；李美鳳等（2016）對四川漢源、茂縣的青花椒和紅花椒進行重金屬含量測定，結果發現鎘含量為0.08～0.20 mg/kg；李聯隊等（2019）對17個花椒林地采集的花椒樣品進行測定，結果發現部分樣地的花椒果皮鎘含量超標，并認為鎘含量超標可能與花椒吸收和富集特性有關；熊汝琴等（2019）對云南昭通4個縣所產的青花椒樣品進行測定，結果發現鐵、錳含量豐富且均不含鎘。【本研究切入點】目前，花椒的相關研究主要集中在種質資源評價與品種選育、有效成分及綜合利用等方面，而有關青花椒種質資源鎘積累特性及富集機制的研究鮮見報道。【擬解決的關鍵問題】以重慶市江津區花椒國家林木種質資源庫中生產性狀較好的青花椒種質資源為試驗材料，測定各部位的鎘含量，分析鎘富集與轉運特征；選取鎘富集能力差異較大的種質，對其根系代謝物進行非靶向代謝組學分析。分析青花椒種質資源鎘積累的差異，探討鎘積累特性及富集機制，為青花椒品種選育和鎘防控提供理論依據。

1材料與方法

1.1試驗材料

供試種質為九葉青花椒（重慶江津）、早熟九葉青花椒（重慶江津）、關嶺藤椒（貴州關嶺）、四川藤椒（四川洪雅）、云南藤椒（云南施甸）和云南大麻椒（云南寶山），均種植保存于重慶市江津區花椒國家林木種質資源庫（29°6′54.38″N～29°7′14.95″N，106°13′5.84″E～106°13′26.06″E），按當地栽培管理方式進行日常管理，根據病蟲害發生情況及時噴藥防治。該地為北高南低的淺丘地貌，海拔295.5～335.7 m，屬北半球中亞熱帶濕潤季風氣候區，常年平均氣溫18.4℃，年平均日照時數1141.0 h，年平均降水量1001.2 mm。土壤為沙溪廟組紫色砂壤土，排水良好，土壤理化性質：pH 5.48、有機質含量11.1 g/kg、全氮含量0.074%、全磷含量1.23 g/kg、全鉀含量26.2 g/kg、陽離子交換量10.3 cmol/kg，土壤重金屬：總砷含量3.91 mg/kg、總鎘含量0.17 mg/kg、總鉻含量46.4 mg/kg、總汞含量0.126 mg/kg、總鉛含量28.0 mg/kg、有效鎘含量0.029 mg/kg。

1.2樣品采集

于2022年6月22日，分別選取長勢良好、無病蟲害、具有代表性的2株各種質植株，在其中部南北兩側各剪取1枝枝條，分解成葉片、枝干、枝皮和果實。將葉片、枝干和枝皮洗凈后烘干、稱重。將果實風干、去雜后分解成果皮和種子。所有樣品分別裝入密封袋保存，用于鎘含量測定。同時，在所選取的2株植株南北側滴水線處采集根系，其中一部分經清洗、烘干后裝入密封袋保存，用于鎘含量測定，另一部分用液氮速凍后于-80℃保存，用于根系代謝物分析。

1.3土壤理化性質及青花椒植株鎘含量測定

土壤pH、有機質、全氮、全磷、全鉀、陽離子交換量、總砷含量、總鎘含量、總鉻含量、總汞含量、總鉛含量、有效鎘含量參照鮑士旦（2000）的方法進行測定。青花椒葉片、枝干、枝皮、果皮、種子和根系樣品的鎘含量測定參照GB 5009.268—2016《食品安全國家標準食品中多元素的測定》，樣品經微波消解后，用電感耦合等離子體質譜法（ICP-MS）測定，所用儀器為NexION 1000電感耦合等離子體質譜儀［珀金埃爾默企業管理（上海）有限公司］。計算各種質的根系鎘富集系數和各部位鎘轉運系數，公式如下：

鎘富集系數=某部位鎘含量/土壤總鎘含量

鎘轉運系數=某部位鎘含量/根系鎘含量

1.4根系代謝物分析

選取鎘富集能力差異較大的種質，對根系代謝物進行非靶向代謝組學分析，委托上海派森諾生物科技股份有限公司進行，數據預處理后使用派森諾基因云（https：//www.genescloud.cn/）進行分析。以變量投影重要性（VIP）gt;1且rlt;0.05為篩選條件，進行差異代謝物鑒定，并對差異代謝物進行KEGG代謝通路富集分析。

1.5統計分析

采用SPSS 19.0對試驗數據進行整理，對青花椒種質資源農藝性狀及鎘含量、鎘富集系數、鎘轉運系數進行單因素方差分析（One-way ANOVA）。利用Excel 2019作圖。

2結果與分析

2.1青花椒種質資源農藝性狀分析結果

6個青花椒種質的果皮、種子、枝干、枝皮和葉片重量分別為17.98～62.29 g/枝、15.10～52.62 g/枝、20.95～48.96 g/枝、15.17～35.45 g/枝、3.87～20.36 g/枝（表1）。果皮、種子、枝干和枝皮重量均以早熟九葉青花椒為最高，且與其他種質差異顯著（Plt;0.05，下同），而均以四川藤椒為最低。葉片重量則以四川藤椒為最高，其后依次為早熟九葉青花椒、云南藤椒、云南大麻椒、九葉青花椒和關嶺藤椒，其中四川藤椒顯著高于其他種質，早熟九葉青花椒與云南藤椒也顯著高于云南大麻椒、九葉青花椒和關嶺藤椒。上述農藝性狀差異可能與各種質的抽梢能力、結實特性等有關。

2.2青花椒種質資源鎘含量測定結果

如圖1所示，6個青花椒種質的果皮、種子、葉片、枝干、枝皮、根系鎘含量分別為0.279～0.579 mg/kg、0.483～0.934 mg/kg、1.413～3.927 mg/kg、0.350～0.829 mg/kg、2.064～4.148 mg/kg、0.590～2.173 mg/kg。不同青花椒種質同一部位的鎘含量存在差異。其中，果皮鎘含量以早熟九葉青花椒為最低，但與四川藤椒無顯著差異（Pgt;0.05，下同）；種子鎘含量以關嶺藤椒為最低，但與早熟九葉青花椒和四川藤椒均無顯著差異；葉片鎘含量以九葉青花椒為最低，且與其他種質差異顯著；枝干鎘含量以九葉青花椒為最低，但與早熟九葉青花椒、關嶺藤椒、四川藤椒均無顯著差異；枝皮鎘含量以關嶺藤椒為最低，但與九葉青花椒、早熟九葉青花椒、四川藤椒、云南大麻椒均無顯著差異；根系鎘含量以早熟九葉青花椒為最低，且與其他種質差異顯著。同時，云南藤椒各部位鎘含量均高于其他種質的同一部位。

此外，同一種質青花椒不同部位的鎘含量存在差異。除關嶺藤椒外，其他種質各部位鎘含量均以枝皮為最高，其次是葉片、根系，且均高于種子、枝干、果皮。關嶺藤椒各部位的鎘含量中，以葉片為最高，其次是枝皮，但二者間無顯著差異。除九葉青花椒外，其他種質各部位鎘含量均以果皮為最低。

2.3青花椒種質資源的鎘富集與轉運特征

鎘富集系數是植物某部位鎘含量與土壤總鎘含量的比值，比值越大說明植物對鎘的富集能力越強。從圖2可知，6個青花椒種質根系的鎘富集系數為3.470～12.780，表明各青花椒種質根系對土壤鎘均有富集能力，富集能力排序為云南藤椒gt;關嶺藤椒gt;云南大麻椒gt;四川藤椒gt;九葉青花椒gt;早熟九葉青花椒，其中云南藤椒顯著強于其他種質，早熟九葉青花椒顯著弱于其他種質。

鎘轉運系數是植物某部位鎘含量與根系鎘含量的比值，比值越大說明植物體內對鎘的轉運能力越強。從圖3-A可知，6個青花椒種質根系到果皮的鎘轉運系數為0.264～0.530，轉運能力排序為九葉青花椒gt;早熟九葉青花椒gt;云南大麻椒gt;四川藤椒gt;云南藤椒gt;關嶺藤椒，且九葉青花椒和早熟九葉青花椒顯著強于其他種質。圖3-B顯示，6個青花椒種質根系到葉片的鎘轉運系數為1.494～3.294，轉運能力排序為早熟九葉青花椒gt;云南大麻椒gt;四川藤椒gt;云南藤椒gt;九葉青花椒gt;關嶺藤椒，其中早熟九葉青花椒顯著強于其他種質資源，關嶺藤椒弱于其他種質。從圖3-C可知，6個青花椒種質根系到枝干的鎘轉運系數為0.324～0.642，轉運能力排序為早熟九葉青花椒gt;云南大麻椒gt;四川藤椒gt;九葉青花椒gt;云南藤椒gt;關嶺藤椒，其中早熟九葉青花椒和云南大麻椒顯著強于其他種質，關嶺藤椒弱于其他種質。圖3-D顯示，6個青花椒種質根系到枝皮的鎘轉運系數為1.408～3.628，轉運能力排序為早熟九葉青花椒gt;九葉青花椒gt;四川藤椒gt;云南大麻椒gt;云南藤椒gt;關嶺藤椒，其中早熟九葉青花椒顯著強于其他種質，關嶺藤椒顯著弱于其他種質。從圖3-E可知，6個青花椒種質根系到種子的鎘轉運系數為0.330～0.943，轉運能力排序為早熟九葉青花椒gt;九葉青花椒gt;四川藤椒gt;云南大麻椒gt;云南藤椒gt;關嶺藤椒，其中早熟九葉青花椒顯著強于其他種質，關嶺藤椒顯著弱于其他種質。

2.4青花椒種質資源根系代謝物分析結果

通過比較各青花椒種質根系的鎘富集系數，發現早熟九葉青花椒與云南藤椒根系對鎘的富集能力差異最大，為探究其富集機制，基于非靶向代謝組學對早熟九葉青花椒與云南藤椒根系代謝物進行分析。正離子模式下的主成分分析（PCA）結果如圖4-A所示，可明顯區分早熟九葉青花椒與云南藤椒樣品，第一主成分（PC1）貢獻率為23.1%，第二主成分（PC2）貢獻率為15.1%；負離子模式下的PCA分析結果如圖4-B所示，早熟九葉青花椒與云南藤椒樣品間差異明顯，PC1貢獻率為24.1%，PC2貢獻率為14.4%。

為進一步驗證早熟九葉青花椒與云南藤椒根系代謝物樣品間的分離情況，采用正交偏最小二乘判別分析（OPLS-DA）對數據進行處理，結果如圖5所示。正離子模式下OPLS-DA模型R2 Y和Q2分別為0.992、0.696，負離子模式下OPLS-DA模型R2 Y和Q2分別為0.999、0.767，均大于0.5，表明模型穩定可靠。與PCA分析結果相比，OPLS-DA分析結果顯示同一種質的樣品聚集更明顯。

以VIPgt;1且Plt;0.05為篩選條件，對早熟九葉青花椒與云南藤椒間的差異代謝物進行鑒定，并將正、負離子模式下鑒定出的差異代謝物進行聚類分析。如圖6-A所示，正離子模式下早熟九葉青花椒與云南藤椒間共鑒定出292個差異代謝物，其中124個代謝物上調表達，168個代謝物下調表達；如圖6-B所示，負離子模式下共鑒定出228個代謝物，其中72個代謝物上調表達，156個代謝物下調表達；合計有520個差異代謝物，其中196個代謝物上調表達，324個代謝物下調表達。對前20種差異代謝物進行分析，發現18種代謝物上調表達，2種代謝物下調表達，其中有2種氨基酸、4種生物堿、6種有機酸、4種糖類、2種脂質、2種酰胺。與早熟九葉青花椒根系代謝物相比，云南藤椒根系代謝物中有機酸、生物堿呈上調表達趨勢。

對早熟九葉青花椒與云南藤椒間的根系差異代謝物進行KEGG代謝通路富集分析，結果表明差異代謝物共注釋到213條通路，其中正離子模式下48條，負離子模式下165條。在正離子模式下，差異代謝物主要富集在異黃酮生物合成（Isoflavonoid biosynthesis）、苯丙氨酸代謝（Phenylalanine metabo-lism）等通路（圖7-A）；在負離子模式下，差異代謝物主要富集在精氨酸生物合成（Arginine biosynthe-sis）、牛磺酸和亞牛磺酸代謝（Taurine and hypotau-rine metabolism）、亞油酸代謝（Linoleic acid metabo-lism）等通路（圖7-B）。

3討論

3.1青花椒種質資源的鎘富集特征

不同植物對鎘的吸收存在較大差異，同一植物的不同品種也存在差異，水稻（Yan et al.，2010；林小兵等，2018）、甘薯（溫國燦，2020）、小麥（潘建清等，2021）、高粱（邵建均等，2022）、玉米（劉奇等，2024）等作物中均存在上述差異。本研究結果表明，九葉青花椒、早熟九葉青花椒、關嶺藤椒、四川藤椒、云南藤椒和云南大麻椒等6個青花椒種質對鎘的吸收存在較大差異，其中云南藤椒的鎘含量高于其他種質。

同一作物不同部位對鎘的富集也存在差異，唐茜等（2008）研究發現，在鎘脅迫下，茶樹各部位的鎘含量排序為根gt;莖gt;新梢（1芽2葉）gt;枝條gt;老葉。本研究結果顯示，6個青花椒種質各部位的鎘含量分布規律不一致；除關嶺藤椒外，其他種質均以枝皮的鎘含量為最高，其次是葉片、根系，且均高于種子、枝干、果皮；同時，除九葉青花椒外，其他種質均以果皮的鎘含量為最低。鎘含量分布規律不一致的原因可能與不同青花椒種質的鎘運輸途徑和積累差異有關。植物吸收鎘后，會根據各自的應激策略將其隔離在細胞的不同部位，如通過與細胞壁的結合及液泡區室化等途徑阻滯鎘的遷移運輸（于輝和向言詞，2018）。彭秋等（2019）研究發現，辣椒根、莖、葉和果各亞細胞組分的鎘含量排序均為細胞壁gt;細胞器gt;細胞可溶性組分，鎘被限制在細胞壁中。唐敏等（2021）研究發現，欒樹、臭椿、銀杏3種喬木葉片中的鎘多存在于細胞壁或可溶性組分中，根系中的鎘多存在于細胞壁中，鎘在葉綠體、線粒體和細胞核等細胞器中分布較少。

鎘富集系數能反映鎘在植物體內的富集情況，富集系數越大，說明植物吸收鎘的能力越強。本研究結果顯示，6個青花椒種質的根系鎘富集系數均大于1，表明各種質根系對土壤鎘有富集能力；同時，各種質間的富集能力存在差異，云南藤椒顯著強于其他種質，早熟九葉青花椒顯著弱于其他種質資源。邵建均等（2022）對不同高粱品種的鎘富集差異進行研究，結果發現6個高粱品種根部鎘富集系數為2.80～5.05，均具有鎘富集特性，本研究結果與之相似。但鎘的運輸和積累是長期持續的過程，且不同發育期對鎘的積累規律可能發生變化，還需進一步對青花椒不同發育期進行系統研究。

3.2青花椒種質資源的鎘富集機制

植物對土壤中鎘的吸收機制與植物種類、品種密切相關。Grant等（2008）研究表明，不同植物品種的鎘積累差異與根系生理活動有關。根系分泌物是植物在生長過程中根系向根際環境釋放或分泌的各類有機化合物，能通過改變根際土壤的酸堿狀況或氧化還原狀況、與重金屬發生螯合或絡合沉淀等化學反應、影響土壤微生物數量和活性等，直接或間接地影響重金屬在土壤中的結合形態及生物有效性，從而導致作物對重金屬積累的差異（陳美靜等，2015）。Cie?liński等（1998）研究發現，不同植物品種的根系分泌物存在差異。謝婕等（2021）對有高、低鎘積累特征的2個蕹菜品種根系分泌物組分進行檢測分析，發現2個蕹菜品種根系分泌物種類和相對含量存在差異，并推測與不同品種間的鎘積累能力差異有關。本研究對早熟九葉青花椒和云南藤椒根系代謝物進行非靶向代謝組學分析，共鑒定出520個差異代謝物，其中196個代謝物上調表達，324個代謝物下調表達，2個種質間存在明顯差異；與早熟九葉青花椒根系代謝物相比，云南藤椒根系代謝物中有機酸、生物堿呈上調表達趨勢。KEGG代謝通路富集分析結果表明，差異代謝物主要富集在異黃酮生物合成、苯丙氨酸代謝、精氨酸生物合成、牛磺酸和亞牛磺酸代謝、亞油酸代謝等通路。相關研究表明，植物根系分泌的有機酸會增加其附近鎘的溶解度和流動性，從而促進植物對鎘的吸收、富集。羅慶等（2016）研究發現，不同生態型東南景天可通過調節根系分泌來耐受或超富集鎘，其可能通過根系分泌更多的癸酸、月桂酸、苯甲酸等來減少鎘的活化從而提升耐受性。Fu等（2018）研究發現，鎘脅迫會促進水稻根系主動分泌有機酸，且在超鎘積累水稻中有機酸的分泌量比正常水稻高出1～2倍。同時，王沛琦等（2021）認為，蓖麻體內積累的草酸、蘋果酸及檸檬酸可能通過根系分泌作用活化環境中的鎘，進而促進對鎘的富集。

3.3青花椒副產物的安全利用

隨著青花椒種植規模和產量的快速增長，青花椒產業鏈不斷延長，其副產物也通過不同方式得到開發利用。如花椒籽可作為肉雞飼料添加劑，提升肉雞生產性能、改良肉質（付興周等，2016）；花椒葉可制成花椒葉茶、口服液、豆瓣醬、醬菜、沐浴液、芳香保健肥皂、空氣清新劑、香煙等（孟資寬等，2021）；花椒枝梢可作為食用菌栽培基質（李曉林等，2014）。本研究結果表明，6個青花椒種質的葉片、枝皮鎘含量相對較高，分別為1.413～3.927 mg/kg、2.064～4.148 mg/kg，種子的鎘含量相對較低，但也達0.483～0.934 mg/kg。因此，在進行青花椒副產物開發利用時，應考慮其相關產品的安全性。

4結論

青花椒種質資源對土壤鎘的積累存在差異，早熟九葉青花椒富集能力較弱，云南藤椒富集能力較強，其富集能力差異可能與根系代謝物有關。

參考文獻（References）：

敖厚豫，李欣，余天華，楊慶雄.2020.基于麻味成分的頂壇花椒HPLC指紋圖譜研究［J］.中國調味品，45（5）：19-26.［Ao H Y，Li X，Yu T H，Yang Q X.2020.Research on HPLC fingerprint spectra of Zanthoxylum planispinum var.dingtanensis based on numb-taste components［J］.China Condiment，45（5）：19-26.］doi：10.3969/j.issn.1000-9973.2020.05.004.

鮑士旦.2000.土壤農化分析［M］.第3版.北京：中國農業出版社.［Bao S D.2000.Soil agrochemical analysis［M］.The 3rd Edition.Beijing：China Agriculture Press.］

陳美靜，劉倩雯，譚佳緣，李雪梅.2015.重金屬脅迫對植物有機酸代謝影響研究進展［J］.廣東農業科學，42（24）：86-91.［Chen M J，Liu Q W，Tan J Y，Li X M.2015.Research advance on effects of heavy metal stress on plant organic acid metabolism［J］.Guangdong Agricultural Sciences，42（24）：86-91.］doi：10.16768/j.issn.1004-874x.2015.24.059.

狄科.2011.花椒中多酚類物質的提取純化及活性研究［D］.南京：南京農業大學.［Di K.2011.Study on extraction，purification and activity of polyphenols from Zathoxlum bungeanum［D］.Nanjing：Nanjing Agricultural University.］doi：10.7666/d.Y2361585.

杜凌，陳應福，侯娜.2018.基于微量元素分析的花椒種質資源評價［J］.種子，37（9）：69-72.［Du L，Chen Y F，Hou N.2018.Germplasm resource evaluation of Zanthoxylum bun-geanum based on microelement analysis［J］.Seed，37（9）：69-72.］doi：10.16590/j.cnki.1001-4705.2018.09.069.

付興周，路志芳，申海燕.2016.花椒籽粉對肉雞生產性能、免疫器官指數及肉品質的影響［J］.河南農業科學，45（9）：120-124.［Fu X Z，Lu Z F，Shen H Y.2016.Effect of Chi-nese prickly ash seed powder on growth performance，immune organ index and meat quality of broilers［J］.Jour-nal of Henan Agricultural Sciences，45（9）：120-124.］doi：10.15933/j.cnki.1004-3268.2016.09.024.

國家藥典委員會.2015.中華人民共和國藥典（一部）［M］.北京：中國醫藥科技出版社：159-160.［National Pharmaco-poeia Commission.2015.Pharmacopoeia of the People’s Republic of China（Volume I）［M］.Beijing：China Medi-cal Science and Technology Press：159-160.］

李聯隊，彌云，謝毓芬，慈順，謝斌，李孟樓.2019.花椒器官中6種重金屬的含量與分布研究［J］.西南林業大學學報（自然科學），39（4）：24-33.［Li L D，Mi Y，Xie Y F，Ci S，Xie B，Li M L.2019.Content and distribution of 6 heavy metals in Zanthoxylum bungeanum［J］.Journal of South-west Forestry University（Natural Sciences），39（4）：24–33.］doi：10.11929/j.swfu.201812010.

李美鳳，陳艷，蔣麗施，孟曉.2016.漢源、茂縣花椒中重金屬的測定［J］.輕工科技，32（5）：14-15.［Li M F，Chen Y，Jiang L S，Meng X.2016.Determination of heavy metals in Zanthoxylum bungeanum Maxim.from Hanyuan County and Maoxian County［J］.Light Industry Science and Tech-nology，32（5）：14-15.］

李曉林，孫海艷，梁國魯，何橋，郭啟高，夏清清，趙申清玉.2014.花椒枝屑栽培金針菇的研究［J］.中國蔬菜，（12）：30-33.［Li X L，Sun H Y，Liang G L，He Q，Guo Q G，Xia Q Q，Zhao S Q Y.2014.Studies on planting Flammulina velutipes at Chinese prickly ash branch［J］.China Vegeta-bles，（12）：30-33.］doi：10.3969/j.issn.1000-6346.2014.12.011.

林小兵，周利軍，王惠明，劉暉，武琳，俞瑩，胡敏，何波，周青輝，黃欠如.2018.不同水稻品種對重金屬的積累特性［J］.環境科學，39（11）：5198-5206.［Lin X B，Zhou L J，Wang H M，Liu H，Wu L，Yu Y，Hu M，He B，Zhou Q H，Huang Q R.2018.Accumulation of heavy metals in diffe-rent rice varieties［J］.Environmental Science，39（11）：5198-5206.］doi：10.13227/j.hjkx.201804027.

劉霈珈.2017.縣域重金屬污染農用地土壤安全利用評價與調控［D］.北京：中國地質大學.［Liu P J.2017.Safe utili-zation evaluation and regulation of agricultural land soil contaminated by heavy mental elements in a Chinese county［D］.China University of Geosciences.］

劉奇，王晟，何濤，包立，張乃明.2024.不同玉米品種對Pb、Cd積累與轉運差異研究［J］.生態與農村環境學報，40（1）：130-137.［Liu Q，Wang S，He T，Bao L，Zhang N M.2024.Differences in lead and cadmium accumulation and translocation in different varieties of Zea mays［J］.Journal of Ecology and rural environment，40（1）：130-137.］doi：10.19741/j.issn.1673-4831.2022.1008.

羅慶，孫麗娜，王輝.2016.基于代謝組學方法的兩種生態型東南景天根系分泌物差異研究［J］.環境化學，35（11）：2353-2361.［Luo Q，Sun LN，Wang H.2016.Study on the differences of root exudates between two ecotypes of Sedum alfredii based on metabolomics［J］.Environmental Chemistry，35（11）：2353-2361.］

馬堯.2021.不同種質資源花椒主要品質性狀差異及其影響因素研究［D］.楊凌：西北農林科技大學.［Ma Y.2021.Differences and influence factors for pericarp quality of different Zanthoxylum germplasm resources［D］.Yangling：Northwest Aamp;F University.］doi：10.27409/d.cnki.gxbnu.2021.001885.

孟資寬，任元元，鄒育，張鑫，張星燦，王擁軍.2021.花椒副產物的研究現狀及展望［J］.食品與發酵科技，57（3）：140-144.［Meng Z K，Ren Y Y，Zou Y，Zhang X，Zhang X C，Wang Y J.2021.Summary of the research status and future prospects of Zanthoxylum bungeanum by-products［J］.Food and Fermentation Sciencesamp;Technology，57（3）：140-144.］doi：10.3969/j.issn.1674-506X.2021.03-022.

慕鈺文，黃玉龍，張輝元.2021.花椒籽不同溶劑提取物的活性成分及抗氧化性比較［J］.包裝與食品機械，39（3）：21-26.［Mu Y W，Huang Y L，Zhang H Y.2021.Bioactive compounds extracted with different solvents from Zan-thoxylum bungeanum Maxim.seeds and comparison of their antioxidant activities［J］.Packaging and Food Machi-nery，39（3）：21-26.］doi：10.3969/j.issn.1005-1295.2021.03.004.

潘建清，陸敏，楊肖娥.2021.不同小麥品種灌漿期生長和鎘積累的差異研究［J］.農業環境科學學報，40（4）：756-765.［Pan J Q，Lu M，Yang X E.2021.Variations in growth and cadmium accumulation among 30 wheat（Triticum aesti-vum L.）cultivars during the grain filling stage［J］.Journal of Agro-Environment Science，40（4）：756-765.］doi：10.11654/jaes.2020-0967.

彭秋，李桃，徐衛紅，焦璐琛，鄧繼寶.2019.不同品種辣椒鎘亞細胞分布和化學形態特征差異［J］.環境科學，40（7）：3347-3354.［Peng Q，Li T，Xu W H，Jiao L C，Deng J B.2019.Differences in the cadmium-enrichment capacity and subcellular distribution and chemical form of cadmium in different varieties of pepper［J］.Environmental Science，40（7）：3347-3354.］doi：10.13227/j.hjkx.20190 1112.

邵建均，樓振邦，徐立軍，謝煒，盧華兵，喻曼.2022.不同高粱品種對土壤重金屬Cd的積累差異［J］.浙江農業科學，63（11）：2500-2502.［Shao J J，Lou Z B，Xu L J，Xie W，Lu H B，Yu M.2022.Study on the accumulation diversity of Cd in different sorghum varieties［J］.Journal of Zhejiang Agricultural Sciences，63（11）：2500-2502.］doi：10.16178/j.issn.0528-9017.20220751.

唐敏，張欣，劉燕，王美仙.2021.鎘在3種喬木中的積累及其亞細胞分布和化學形態研究［J］.環境科學學報，41（6）：2440-2447.［Tang M，Zhang X，Liu Y，Wang M X.2021.Accumulation，subcellular distribution and chemical forms of cadmium in three tree species［J］.Acta Scientiae Cir-cumstantiae，41（6）：2440-2447.］doi：10.13671/j.hjkxxb.2020.0413.

唐茜，葉善蓉，陳能武，張蘇敏，譚和平.2008.茶樹對鉻、鎘的吸收積累特性研究［J］.茶葉科學，28（5）：339-347.［Tang Q，Ye S R，Chen N W，Zhang S M，Tan H P.2008.Study on the absorption and accumulation of Cr3+and Cd2+in tea plants［J］.Journal of Tea Science，28（5）：339-347.］doi：10.3969/j.issn.1000-369X.2008.05.005.

王紀磊，李霜，史慶華.2020.花椒籽（餅）熏蒸對溫室根結線蟲防效及番茄生長發育的影響［J］.山東農業科學，52（1）：120-125.［Wang J L，Li S，Shi Q H.2020.Effect of Zanthoxylum bungeanum seeds（cake）fumigation on con-trol of root-knot nematode as well as growth and develop-ment of tomato in greenhouse［J］.Shandong Agricultural Sciences，52（1）：120-125.］doi：10.14083/j.issn.1001-4942.2020.01.022.

王林，胡金祥，喬明鋒，何蓮，魯道云，吳華昌，鄧靜，易宇文.2022.南方4產區紅花椒揮發性物質鑒定及差異研究［J］.中國調味品，47（6）：165-170.［Wang L，Hu J X，Qiao M F，He L，Lu D Y，Wu H C，Deng J，Yi Y W.2022.Research on identification and differences of volatile com-pounds in Zanthoxylum bungeanum Maxim.from four pro-ducing areas in South China［J］.China Condiment，47（6）：165-170.］doi：10.3969/j.issn.1000-9973.2022.06.031.

王沛琦，胡尊紅，胡學禮，劉旭云，楊謹，李文昌.2021.鎘脅迫對蓖麻有機酸含量及鎘吸收的影響［J］.山西農業科學，49（7）：822-827.［Wang P Q，Hu Z H，Hu X L，Liu X Y，Yang J，Li W C.2021.Effect of cadmium stress on organic acid content and cadmium absorption of Ricinus communis L.［J］.Journal of Shanxi Agricultural Sciences，49（7）：822-827.］doi：10.3969/j.issn.1002-2481.2021.07.06.

溫國燦.2020.不同甘薯品種對土壤重金屬鎘轉移富集規律研究［J］.福建農業科技，7：20-24.［Wen G C.2020.Study on the transfer and enrichment regulation of different varie-ties of sweet potato to the heavy metal cadmium in soil［J］.Fujian Agricultural Science and Technology，7：20-24.］doi：10.13651/j.cnki.fjnykj.2020.07.004.

席少陽，郭延秀，馬曉輝，詹志來，晉玲.2021.花椒化學成分及藥理作用的研究進展［J］.華西藥學雜志，36（6）：717-722.［Xi S Y，Guo Y X，Ma X H，Zhan Z L，Jin L.2021.Research progress on chemical constituents and pharmaco-logical effects of Zanthoxylum［J］.West China Journal of Pharmaceutical Sciences，36（6）：717-722.］doi：10.13375/j.cnki.wcjps.2021.06.027.

謝婕，李建杏，龔玉蓮，曾碧健，余樂洹，黃慧婷，唐鍵.2021.蕹菜不同鎘積累品種根系分泌物化學組分分析［J］.中國蔬菜，（12）：73-78.［Xie J，Li J X，Gong Y L，Zeng B J，Yu L H，Huang H T，Tang J.2021.Analysis on root exudates of different cadmium-accumulating cultivars of water spi-nach［J］.China Vegetables，（12）：73-78.］doi：10.19928/j.cnki.1000-6346.2021.1050.

熊汝琴，路紹彩，王艷.2019.昭通地產青花椒微量元素含量對比分析［J］.浙江農業科學，60（5）：813-815.［Xiong R Q，Lu S C，Wang Y.2019.Contrastive analysis of trace ele-ments contents in Zanthoxylum schinifolium Sieb.et Zucc.in Zhaotong［J］.Journal of Zhejiang Agricultural Sciences，60（5）：813-815.］doi：10.16178/j.issn.0528-9017.201905 44.

徐雙美，楊森，李澤娜，楊丹，陳鴻平，劉友平.2020.GC-MS分析四川主產地竹葉花椒的揮發性成分［J］.中國調味品，45（12）：124-129.［Xu S M，Yang S，Li Z N，Yan D，Chen H P，Liu Y P.2020.Analysis of volatile components from Zanthoxylum armatum DC.in Sichuan Province by GC-MS［J］.China Condiment，45（12）：124-129.］doi：10.3969/j.issn.1000-9973.2020.12.026.

尹慶豐，王毅，呼延麗，郝佳波，蒙進芳，陸斌.2023.基于簡化基因組測序的野生竹葉花椒遺傳多樣性分析［J/OL］.分子植物育種.https：//link.cnki.net/urlid/46.1068.s.20230821.1644.015.［Yin Q F，Wang Y，Huyan L，Hao J B，Meng J F，Lu B.2023.Genetic diversity of wild Zanthoxylum armatum by ddRAD-seq［J/OL］.Molecular Plant Breeding.https：//link.cnki.net/urlid/46.1068.s.20230821.1644.015.］

于輝，向言詞.2018.鎘不同積累型作物品種對鎘吸收、轉運和積累特性的研究進展［J］.生態科學，37（2）：223-228.［Yu H，Xiang Y C.2018.Advance in study of cadmium uptake，translocation and accumulation in different crop varieties［J］.Ecological Science，37（2）：223-228.］doi：10.14108/j.cnki.1008-8873.2018.02.028.

于霞.2016.花椒籽廢渣活性炭的制備、改性及吸附性能研究［D］.蘭州：西北師范大學.［Yu X.2016.Synthesis and modification of activated carbon with the waste residue of Chinese prickly ash seeds and its adsorption properties［D］.Lanzhou：Northwest Normal University.］

余曉琴.2010.花椒品質評價方法及其應用研究［D］.重慶：西南大學.［Yu X Q.2010.Study on quality evaluation system of Pericarpium zanthoxyli and its application［D］.Chongqing：Southwest University.］

趙秀玲.2012.花椒的化學成分、藥理作用及其資源開發的研究進展［J］.中國調味品，37（3）：1-5.［Zhao X L.2012.Research progress in chemical constituents，biological activities and exploration utilization of Pericarpium zan-thoxyli［J］.China Condiment，37（3）：1-5.］doi：10.3969/j.issn.1000-9973.2012.03.001.

Cie?liński G，Van Rees K C J，Szmigielska A M，Krishnamurti G S R，Huang P M.1998.Low-molecular-weight organic acids in rhizosphere soils of durum wheat and their effect on cadmium bioaccumulation［J］.Plant and Soil，203：109-117.doi：10.1023/A：1004325817420.

Fu H J，Yu H Y，Li T X，Zhang X Z.2018.Influence of cad‐mium stress on root exudates of high cadmium accumula-ting rice line（Oryza sativa L.）［J］.Ecotoxicology and Environmental Safety，150：168-175.doi：10.1016/j.ecoenv.2017.12.014.

Grant C A，Clarke J M，Duguid S，Chaney R L.2008.Selection and breeding of plant cultivars to minimize cadmium accu‐mulation［J］.Science of the Total Environment，390（2-3）：301-310.doi：10.1016/j.scitotenv.2007.10.038.

Yan Y F，Choi D H，Kim D S，Lee B W.2010.Absorption，translocation，and remobilization of cadmium supplied at different growth stages of rice［J］.Journal of Crop Science and Biotechnology，13：113-119.doi：10.1007/s 12892-010-0045-4.

（責任編輯：陳燕，劉可丹）