利用离子组学和代谢组学研究硒肥对烤烟元素及化学成分的影响

范晓苏; 沈方科; 韦建玉; 石保峰; 贾海江; 路丹; 顾明华

doi:10.12101/j.issn.1004-390X(n).202209050

利用离子组学和代谢组学研究硒肥对烤烟元素及化学成分的影响

1.
广西大学农学院，植物科学国家级实验教学示范中心，广西南宁 530004
2.
广西中烟工业有限责任公司，广西南宁 530001
3.
广西壮族自治区烟草公司贺州市公司，广西贺州 542899

基金项目: 广西重点研发项目(桂科AB16380160)；广西中烟工业有限公司科技项目(20170206)。

详细信息

作者简介:
范晓苏(1988—)，女，广西来宾人，硕士，实验师，主要从事质谱和光谱仪器分析研究。E-mail：fanxs@gxu.edu.cn

通信作者:
沈方科(1973—)，男，广西博白人，硕士，高级农艺师，主要从事植物营养生理与施肥技术研究。E-mail：shenfk@gxu.edu.cn

中图分类号: S184
计量
- 文章访问数: 1528
- HTML全文浏览量: 349
- PDF下载量: 13
出版历程
- 收稿日期: 2022-09-26
- 修回日期: 2023-05-25
- 录用日期: 2023-05-28
- 网络出版日期: 2023-05-30
- 发布日期: 2023-05-29

摘要:

目的探究不同外源硒肥在不同施肥方式下对富硒烤烟元素及化学活性成分含量的影响。

方法以云烟87为供试材料，设置不施用外源硒(CK)、根际施用纳米硒(S₁)、叶面喷施纳米硒33.75 g/hm² (S₂)和叶面喷施纳米硒67.50 g/hm² (S₃)处理，利用电感耦合等离子体质谱仪、电感耦合发射光谱仪和超高效液相色谱仪—四级杆飞行时间质谱联用仪研究烤烟的组学差异。

结果4组烤烟的元素分布存在明显差异，主要做出贡献的元素是S、Mg、Mn、Al、P和Na等。4组烤烟的代谢组学也有明显差异，鉴定出29种差异代谢物，其中植物甾醇及其衍生物4种，酚类化合物9种，氨基酸类化合物7种，醛类化合物3种及其他物质6种。叶面喷硒肥对烤烟元素及化学活性成分的影响大于根际施硒肥。

结论研究结果为富硒烤烟元素及化学活性成分吸收与积累的机制研究提供了参考。

关键词:

硒肥 /
烤烟 /
离子组 /
代谢组

Effects of Selenium Fertilizer on Elements and Chemical Components of Flue-cured Tobacco Using Metallomics and Metabolomics Method

1.
College of Agriculture, Guangxi University, National Demonstration Center for Experimental Plant Science Education, Nanning 530004, China
2.
China Tobacoo Guangxi Industrial Company Limited, Nanning 530001, China
3.
Guangxi Zhuang Automomous Region Tobacoo Company, Hezhou Branch, Hezhou 542899, China

Abstract:

Purpose To study the effects of different fertilization methods and different selenium enrichment levels on the contents of elements and active components in flue-cured tobacco.

MethodsYunyan87 was used as test material, and the following treatments were applied: no application of exogenous selenium application (CK), application of selenium nanoparticles in the rhizosphere (S₁), foliar spray application of selenium nanoparticles with 33.75 g/hm² (S₂) and foliar spray application of selenium nanoparticles with 67.50 g/hm² (S₃). Compositional and multiple omics differences between flue-cured tobacco under different treatments were studied using inductively coupled plasma mass spectrometer, inductively coupled optical emission spectrometer, and ultra-performance liquid chromatography-quadruple time-of-flight mass spectrometer.

ResultsThere were obvious differences of element distribution among the four groups of flue-cured tobacco, with the main contributing elements being S, Mg, Mn, Al, P, and Na. The metabolomics of the four groups also displayed obvious differences, with a total of 29 differential metabolites identified, including four phytosterols and their derivatives, nine phenolic compounds, seven amino acid compounds, three aldehyde compounds, and six other compounds. The effect of foliar application of selenium fertilizer on the elements and chemically active components of flue-cured tobacco was found to be more significant than that of rhizosphere application.

ConclusionThese results may serve as a reference for investigating the absorption and bioaccumulation mechanism of selenium-enriched flue-cured tobacco’s elements and chemically active components.

Keywords:

HTML全文

干旱和半干旱地区的流动沙地是一种地表裸露疏松的生态系统，植被覆盖度低，风沙活动强烈，降水稀少，蒸发强烈^[1-2]。土壤水分被认为是控制沙漠化的约束因素^[3-4]，是气候、植被、地形及土壤等自然条件的综合反映，其时空分异性在大尺度上受气象和水文状况影响，局部地貌特征对土壤水分分异规律的影响也不可小觑^[5]。国内外学者就土壤水分时空动态及其运动机理^[6]，降水、植被类型和生物结皮对土壤水分的影响等内容进行了大量研究^[7-9]。土壤水分在季节性积雪、地形地貌、植被盖度等因素的影响下显现出明显的时空分异性，表现为在较小尺度上水分随坡度变幅明显^[10-11]。前人关于乌兰布和沙漠和毛乌素沙地的研究^[12-13]表明：沙丘上土壤含水量呈现出由丘间地向丘顶逐渐减小的趋势。在同样的地貌部位，决定土壤含水量的主导因素则变为植被类型及其盖度，灌木类植物的保水能力更强^[14]。植被长势较好的固定沙丘，由于生物结皮的存在，导致水分难以下渗而集中在表层，深层土壤含水量偏低，随沙地植被生长时间的增加，沙地水分逐渐降低，但当植被的生长时间超过40年时，沙地浅层保水能力显著提升^[15-16]。

近年来，宁夏河东沙地特别是白芨滩地区防沙治沙工作取得了重大成果，但是缺乏对固沙植被建设以来土壤—植被相互作用的深入研究，尤其是随着固沙植被的生长，土壤水分在植被—土壤系统演替中的机理还不明确；对于该地区土壤水分的研究也仅是短期观测，缺乏长期定位监测与系统分析研究。本研究在该地区过去研究的基础上，进一步在野外定位观测土壤水分动态变化及其格局，揭示宁夏白芨滩防沙治沙区浅层土壤水分的空间变异特征，同时也可为该区域沙漠化治理生态效益评估提供基础数据和科学依据。

1. 材料与方法

1.1 研究区概况

如图1所示：宁夏河东沙地白芨滩自然保护区的大泉沙漠化治理区位于灵武市以南 40 km，N37°54′~38°22′，E106°23′~106°48′，海拔 1 150~1 650 m。属中温带干旱区，大陆性气候特征显著，全年太阳总辐射量为 6 056.20 MJ/m²，全年日照时间为 3 004.80 h，平均气温为 8.8 ℃，年平均降雨量为 198.9 mm，年内分配不均，多集中在 7、8 月，蒸发量大，年平均蒸发量为 1 928.4 mm。白芨滩位于毛乌素沙地南缘，鄂尔多斯台地西南缘，分布着干旱沙地、荒漠和流动沙丘等属荒漠类型的生态系统。白芨滩沙丘类型以流动沙丘为主，呈新月形沙丘链状分布，占 83%；半固定沙丘多呈椭圆形，占10%；流动浮沙土集中分布在流动沙丘的丘间平地，约占 5%；淡灰钙土集中分布在北部流动沙丘间的丘间平地，约占 2%。植被以沙生植被为主，如沙蒿(Artemisia desertorum )、沙米(Agriophyllum squarrosum)、苦豆子(Sophora alopecuroidesL.)、柠条锦鸡儿(Caragana korshinskii Kom.)和沙柳(Salix psammophila)等，植被盖度为 2%~20%^[17-18]，该沙区主要固沙植物种为柠条锦鸡儿(Caragana korshinskii Kom.)和花棒(Hedysarum scoparium Fisch.)等灌木植物。

图 1 研究区位置图和地形图

Figure 1. Location map and topographic map of the study area

下载: 全尺寸图片幻灯片

1.2 研究方法

2014年4月，在宁夏灵武白芨滩大泉林场防沙治沙区，选择已采用草方格和植物措施进行固定的沙丘(固沙时间为2008年)，用测线圈定100 m×100 m的范围作为样地，N38°0′18″, E106°24′15″。样地中包括迎风坡、丘间地和背风坡等地形，在其中布设 11 条样带，均匀设置10 m×10 m的样方 120个，取约0~30 cm的浅层风沙土。从2014年4月28日开始，分别在4月底、6月底、8月底和9月底测定121个样点的土壤含水量。其中表层10 cm深处的土壤含水量采用烘干法测定，30 cm深处的土壤含水量采用TDR300便携式土壤水分速测仪测定。前者是重量含水量，后者是体积含水量。

1.3 数据处理

采用SPSS 20.0对试验数据进行描述性统计，利用GS+10.0地统计学分析软件分析土壤含水量空间异质性，建立半方差函数模型，并进行克里格插值。块金值(C₀)与基台值(C₀+C)的比值即块金系数可表示系统变量空间相关程度，块金系数< 25 %表明空间相关性很强，块金系数在25%~50%之间说明变量有明显的空间自相关，50%~75%说明具有中等空间相关性，块金系数 > 75%表明空间相关性很弱^[19]。

2. 结果与分析

2.1 土壤含水量均值及变异系数

从2014年对白芨滩治沙区土壤含水量的观测发现：受降雨、蒸发、植物生长及微地形的影响，白芨滩防沙区表层以及30 cm深处的土壤含水量具有较明显的季节变化趋势(表1、2)，而且表层与30 cm深处的土壤含水量均表现出明显的空间分布特征(图2、3)。从表1可见：4月底、6月底和8月底土壤含水量较低，而9月底土壤含水量较高。30 cm的土壤含水量表现出4月底、6月底小于8月底、9月底。该研究区土壤含水量主要受降水和蒸发影响，因此季节变化较为明显，8月和9月降水较4月和6月多，该地区土壤水分主要来源于降水。变异系数CV的大小反映了土壤含水量的空间变异性大小，一般认为CV<0.01为弱变异性，0.01<CV<1为中等变异性，CV>1为强变异性^[20]。从表1中可以看出，10 cm和30 cm的土壤含水量变异系数均在0.01~1之间，属中等变异。

表 1 土壤含水量(白芨滩)描述性统计结果

Table 1. Descriptive statistical results on soil moisture content (Baijitan)

测定日期 detection date	深度/cm depth	(平均值±标准差)/% average ± standard deviation	最小值/% minimum	最大值/% maximum	极差/% range	变异系数/% coefficient of variation	分布类型 distribution type
2014−04−30	10	3.22±1.42	0.99	4.38	3.39	0.44	正态分布 normal distribution
2014−04−30	30	15.4±1.41	10.00	17.40	7.40	0.09	正态分布 normal distribution
2014−06−30	10	2.15±1.55	0.38	4.82	4.44	0.72	正态分布 normal distribution
2014−06−30	30	10.9±0.42	8.20	13.50	5.30	0.04	正态分布 normal distribution
2014−08−30	10	3.41±0.87	1.06	4.02	2.96	0.26	正态分布 normal distribution
2014−08−30	30	16.4±0.707	13.00	18.40	5.40	0.04	正态分布 normal distribution
2014−09−30	10	7.47±0.56	2.09	10.16	8.07	0.07	正态分布 normal distribution
2014−09−30	30	20.35±2.05	10.20	23.80	13.60	0.10	正态分布 normal distribution

下载: 导出CSV

| 显示表格

表 2 土壤含水量变异函数模型

Table 2. Variation function models of soil moisture content

测定日期 detection date	深度/cm depth	理论模型 theoretical model	块金值/C₀ nugget	基台值/C₀+C sill	变程/m variable range	块金系数/%(块金值/基台值) coefficient of nugget (nugget/ sill)	决定系数/(r²) coefficient of determination	残差/RSS residual
2014−04−30	10	球状模型 spherical model	0.016 00	0.516 000	17.00	3.10	0.61	1.850E-03
2014−04−30	30	指数模型 index model	0.238 00	2.664 000	19.20	8.94	0.56	0.08
2014−06−30	10	球状模型 spherical model	0.001 00	0.584 000	16.40	0.17	0.29	6.690E-03
2014−06−30	30	指数模型 index model	0.064 00	0.626 000	24.60	10.22	0.53	0.01
2014−08−30	10	指数模型 index model	0.019 00	0.343 000	16.50	5.54	0.41	1.210E-03
2014−08−30	30	指数模型 index model	0.109 00	1.090 000	14.10	10.00	0.08	0.04
2014−09−30	10	球状模型 spherical model	0.046 00	1.190 000	14.90	0.96	0.09	0.03
2014−09−30	30	指数模型 index model	0.203 00	2.800 000	20.40	0.93	0.58	0.18

下载: 导出CSV

| 显示表格

图 2 土壤含水量(10 cm)分布半方差图

Figure 2. Semi-variance of soil moisture content (10 cm) distribution

下载: 全尺寸图片幻灯片

图 3 土壤含水量(30 cm)分布半方差图

Figure 3. Semi-variance of soil moisture content (30 cm) distribution

下载: 全尺寸图片幻灯片

本研究所选取的沙丘类型为固定沙丘，表层土壤含水量表现出明显的空间分异特征。李森等^[21]对腾格里沙漠东南缘固定沙丘表层土壤水分时空分异的研究表明：固定沙丘表层含水量时空分异之所以表现明显有两方面原因：降水的分布格局和植物蒸腾及土壤表面蒸散耗水。研究结果中10 cm和30 cm处4月底、6月底土壤含水量均值小于8月底和9月底，表明4月底和6月底是土壤水分损耗期，8月底和9月底是土壤水分积累期。朱海等^[22]在古尔班通古特沙漠南缘对固定沙丘土壤水分的研究认为3—5月是土壤水分补给期，6—10月是土壤水分耗损期。这一研究结果与白芨滩呈相反趋势，白芨滩年内降雨集中在7、8月，4、6月耗水源主要是植物生长和土壤表层蒸散，而古尔班通古特沙漠土壤含水量补给方式为春季冰雪融水和降雨，5月以后，温度回升，蒸发加剧，再加上大量短命植物生长耗水，土壤水分损失严重。对于乌兰布和流动沙地的研究得出，其水分含量主要随降雨量变化^[2]。本研究发现：除9月底以外其他时间段10 cm处的变异系数均大于30 cm处，土壤含水量的变异系数有随土层深度增加而减小的趋势，由于降水和蒸发的影响浅层土壤水分的变异性较大^[5]。

2.2 土壤含水量空间分布特征

利用GS+10.0地统计学软件做出的结果(表2和图2、3)表明：4、6、9月的土壤水分(10 cm)变异函数分别为球状模型，8月(10 cm)为指数模型；30 cm的土壤水分变异函数均为指数模型。由表2可知：不同时间的表层土壤含水量及30 cm深处的块金系数均小于25%，土壤含水量均表现出强烈的空间自相关。

从表2和图2可见：在研究区内10 cm的土壤含水量空间自相关范围在各月较接近，说明不同季节土壤水分的生态过程变化不明显，原因是该地区土壤水分主要受到降雨、蒸发的影响，比植物耗水更大。变程是指变异函数达到基台值所对应的距离，它反映的是属性因子空间自相关范围的大小^[18]。由表2可知：30 cm的土壤水分空间自相关范围在14.1~24.6 m之间，不同季节30 cm深处的土壤含水量变程范围明显不同，4、6和9月变程均大于10 cm深处的变程，表明30 cm深处在4、6和9月的土壤含水等空间变异程度低于10 cm深处的变程。其中8月最小，6月最大，变程逐渐减小，自相关范围的减小表明土壤含水量分布的均一性减弱，即8月30 cm深处的土壤含水量受降雨影响变化最明显。

分析空间局部插值估计进行克里格插值制图可以精确了解变量的空间分布格局。由图4可知：4月底表层土壤含水量整体较低，分布相对较均匀。6月底研究区土壤含水量较其他月都低，整体分布相对较均匀，该时期无降水补充，蒸发强烈，对表层土壤含水量影响非常大，土壤含水量较低。8月底土壤含水量较6月增加，空间分布表现为从北向南降低，其中西北角土壤含水量最高，此处地形较低，为丘间地，容易集水。9月底研究区内土壤表层含水量较其他时期都高，分布较均匀。该区域降水集中于8、9月，该时期植被生长较茂盛，可以减少土壤表层蒸发。不同时期表层土壤含水量的空间异质性明显不同，这表明不同季节的空气温度、降雨、地形以及植被生长情况都会影响表层土壤含水量的空间分布特征。

图 4 不同季节表层(10 cm)土壤含水量空间分布图

Figure 4. Spatial distribution of soil moisture content in surface layer (10 cm) in different seasons

下载: 全尺寸图片幻灯片

由图5可见：不同季节30 cm的土壤含水量空间分布特征，显示该深度土壤含水量的空间变异受地形、季节影响明显，特别是4月底和6月底的土壤含水量空间分布图相似。30 cm深处的土壤含水量随季节发生变化，其空间分布主要受地形影响，丘间地等低洼地土壤含水量较高，与表层土壤含水量空间分布特征完全不同。

图 5 不同季节土壤含水量(30 cm)空间分布图

Figure 5. Spatial distribution of soil moisture content (30 cm) in different seasons

下载: 全尺寸图片幻灯片

3. 讨论

土壤含水量受众多因素影响，包括降水、蒸发、土壤质地、地形等。通常在干旱区浅层土壤含水量受降水和蒸发影响更明显。研究表明：0~60 cm是土壤水分受降水影响的活跃层^[2]。我们详细研究了白芨滩防沙治沙区浅层土壤不同深度水分在不同季节的空间变化特征，对该地区生态恢复有重要的理论指导意义。不同季节表层(10 cm)土壤水分空间分布特征明显不同。本研究区属于沙区，降水是土壤水分的主要来源，因此土壤含水量的变化与降水的时空分布密切相关，降水的季节性分布导致土壤水分的季节性变化。研究发现：降水时沙地含水量高于平均值，降水后含水量迅速下降，但降水只对0~40 cm的土壤含水量影响显著^{[2, 5, 23-24]}。降水时，水分既要克服蒸发消耗又要克服表土层的吸附截留，才能下渗到植物根系层，有效地补充该层土壤水分^[25]。除了降水，气象因素中还包括土壤温度和风况，都对沙地水分含量产生影响^[16]。其次，地形也是影响土壤含水量空间分布异质性的一个重要因素，研究得出该区域30 cm深处的土壤含水量空间分布特征受地形影响明显。另外，植被也对土壤水分空间分异影响比较明显。对毛乌素沙地与科尔沁沙地土壤水分的研究表明：地貌部位、植被类型及覆盖度是影响沙地土壤含水量的重要因素^{[14, 26-28]}。与盐池县同期结果对比^[20]，表层土壤含水量的空间变化趋势完全相反，白芨滩的空间分布特征表现出的异质性较小，分布相对均匀，30 cm处的土壤含水量两地都主要受地形影响。造成两地表层土壤含水量的空间分布差异的主要原因是降水量不同，两个研究区植被类型和覆盖度不同，选取样地地形不同。

本研究中表层土壤含水量为质量分数，30 cm深处土壤含水量为体积分数，因此，本研究中没有分析垂直分布规律，仅是探讨浅层土壤含水量空间分布特征，在后续工作中急需深入开展关于土壤含水量深层垂直分布规律的研究。此外，本研究属小尺度空间分布研究，土壤含水量的空间异质性在较小尺度上，主要受到植被状况、土壤特性和微地形等因素的影响^[3]。有关土壤水分空间变异的研究较多，而且各区域不同尺度研究结果也存在一定差异。最理想的采样是非常大的采样幅度、非常小的采样间距以及采样体积^[29]。因此，对本区域还需进一步开展不同尺度的空间变异特征研究。

4. 结论

(1)研究区10 cm土壤含水量和30 cm土壤含水量的变异程度中等，不同季节的土壤含水量表现出强烈的空间自相关，而且10~30 cm含水量的分布格局均受降水影响。10 cm和30 cm的土壤含水量8、9月相对较多，9月含水量是所有研究时期中最高的。

(2)研究区不同时期表层(10 cm)土壤含水量空间分布明显不同，30 cm深处的土壤含水量随季节发生变化，而且其空间分布特征受地形、季节影响明显。在丘间低地，不同深度的土壤含水量均表现出增加趋势。

图 1 烤烟元素主成分分析得分图

注：CK. 不施用外源硒，S₁. 根际施用纳米硒，S₂. 叶面喷施纳米硒33.75 g/hm²，S₃. 叶面喷施纳米硒67.50 g/hm²；下同。

Figure 1. Principal component analysis (PCA) score plot of flue-cured tobacco

Note: CK. no exogenous selenium application, S₁. application of selenium nanoparticles in the rhizosphere, S₂. foliar spray application of selenium nanoparticles with 33.75 g/hm², S₃. foliar spray application of selenium nanoparticles with 67.50 g/hm²; the same as below.

下载: 全尺寸图片幻灯片

图 2 烤烟元素主成分分析载荷图

Figure 2. PCA loading plot of flue-cured tobacco

下载: 全尺寸图片幻灯片

图 3 对照组(a)、根际施肥组(b)以及叶面喷施组(c和d)样本的负离子模式总离子流图

Figure 3. Negative ionize mode total ion current chromatography of control (a), rhizosphere fertilizer treatment (b) and foliar spray treatments (c and d)

下载: 全尺寸图片幻灯片

图 4 对照组(a)、根际施肥组(b)以及叶面喷施组(c和d)样本的正离子模式总离子流图

Figure 4. Positive ionize mode total ion current chromatography of control (a), rhizosphere fertilizer treatment (b) and foliar spray treatments (c and d)

下载: 全尺寸图片幻灯片

图 5 不同硒肥处理下各处理间PCA图

注：a) 正离子模式；b) 负离子模式。

Figure 5. PCA plot of different treatments under different selenium fertilization

Note: a) positive ionize mode; b) negative ionize mode.

下载: 全尺寸图片幻灯片

表 1 微波消解升温程序

Table 1 Temperature program of microwave digestion

步骤 step	设定温度/℃ setting temperature	恒温时间/min constant temperature time	功率/W power
1	120	6	1200
2	150	3	1200
3	180	5	1200
4	200	30	1200

下载: 导出CSV

表 2 待测元素同位素信息表

Table 2 Information of element to be measured and internal standard isotopes

序号 No.	元素 element	测量同位素 isotope internal standards	序号 No.	元素 element	测量同位素 isotope internal standards
1	Li	⁷Li	32	Pd	¹⁰⁶Pd
2	Be	⁹Be	33	Ag	¹⁰⁷Ag
3	Ti	⁴⁷Ti	34	Sn	¹¹⁸Sn
4	V	⁵¹V	35	Te	¹³⁰Te
5	Cr	⁵²Cr	36	Cs	¹³³Cs
6	Mn	⁵⁵Mn	37	Ba	¹³⁸Ba
7	Fe	⁵⁷Fe	38	La	¹³⁹La
8	Co	⁵⁹Co	39	Ce	¹⁴⁰Ce
9	Ni	⁶⁰Ni	40	Pr	¹⁴¹Pr
10	Cu	⁶³Cu	41	Nd	¹⁴²Nd
11	Zn	⁶⁶Zn	42	Sm	¹⁵²Sm
12	As	⁷⁵As	43	Eu	¹⁵³Eu
13	Ru	¹⁰²Ru	44	Gd	¹⁵⁸Gd
14	Sr	⁸⁸Sr	45	Dy	¹⁶⁴Dy
15	Mo	⁹⁸Mo	46	Ho	¹⁶⁵Ho
16	Cd	¹¹¹Cd	47	Er	¹⁶⁶Er
17	Sb	¹²¹Sb	48	Tm	¹⁶⁹Tm
18	Tl	²⁰⁵Tl	49	Yb	¹⁷⁴Yb
19	Pb	²⁰⁸Pb	50	Lu	¹⁷⁵Lu
20	B	¹¹B	51	Hf	¹⁸⁰Hf
21	Na	²³Na	52	Ta	¹⁸¹Ta
22	Mg	²⁴Mg	53	W	¹⁸⁴W
23	Al	²⁷Al	54	Re	¹⁸⁷Re
24	P	³¹P	55	Os	¹⁹²Os
25	S	³²S	56	Ir	¹⁹³Ir
26	Ca	⁴³Ca	57	Pt	¹⁹⁵Pt
27	Ga	⁶⁹Ga	58	Au	¹⁹⁷Au
28	Rb	⁸⁵Rb	59	Hg	²⁰²Hg
29	Zr	⁹⁰Zr	60	Pa	²³¹Pa
30	Nb	⁹³Nb	61	Th	²³²Th
31	Rh	¹⁰³Rh	62	U	²³⁸U

下载: 导出CSV

表 3 标准曲线参数

Table 3 Standard curve parameters

元素 element	质量浓度范围/(mg·L⁻¹) range of mass concentration	方程式 regression equation	相关系数 correlation index
S	1~50	y=0.01925x−0.58511	0.9901
Mg	1~50	y=5.554 08e⁻⁷x−0.05930	0.9999
Al	1~50	y=1.637 02e⁻⁶x−0.04165	0.9998
P	1~50	y=0.22530x−3.94446	0.9986

下载: 导出CSV

表 4 不同处理的硒含量

Table 4 Selenium content in different treatments

处理 treatments	硒含量/(mg·kg⁻¹) selenium content
CK	0.031±0.001 Aa
S₁	0.240±0.031 Bb
S₂	0.231±0.048 Bb
S₃	0.441±0.084 Cc
注：CK. 不施用外源硒，S₁. 根际施用纳米硒，S₂. 叶面喷施纳米硒33.75 g/hm²，S₃. 叶面喷施纳米硒67.50 g/hm²；下同。不同大写字母表示差异极显著 (P<0.01)，不同小写字母表示差异显著 (P<0.05)。 Note: CK. no exogenous selenium application, S₁. application of selenium nanoparticles in the rhizosphere, S₂. foliar spray application of selenium nanoparticles with 33.75 g/hm², S₃. foliar spray application of selenium nanoparticles with 67.50 g/hm²; the same as below. Different uppercase letters indicate extremely significant differences (P<0.01), different lowercase letters indicate significant differences (P<0.05).

下载: 导出CSV

表 5 组间差异元素比较

Table 5 Comparison of different elements between groups

组间比较 comparison between groups	差异元素 different element	变量投影重要度 variable importance for the projection
CK vs S₁	S	6.88073
	Mg	3.05575
	P	1.13191
CK vs S₂	Mg	5.75246
	S	4.52645
	Mn	1.91015
	Na	1.25467
	Al	1.23085
CK vs S₃	S	5.64746
	Mg	4.77327
	Al	1.58081
	Mn	1.53345
S₁ vs S₂	Mg	6.28388
	S	3.50569
	Mn	2.03379
	Na	1.53989
	Al	1.33968
S₁ vs S₃	Mg	5.47132
	S	4.73899
	Al	1.77892
	Mn	1.67301
S₂ vs S₃	Mg	6.13279
	S	3.80198
	Mn	2.09729
	Na	1.74293

下载: 导出CSV

表 6 差异元素在富硒烤烟样本中的含量

Table 6 Content of different element in tobacco samples

元素 element	Mn/(mg·kg⁻¹)	S/(g·kg⁻¹)	Mg/(g·kg⁻¹)	Al/(g·kg⁻¹)	P/(g·kg⁻¹)	Na/(mg·kg⁻¹)
CK	40.20	1.80	1.95	0.079	1.81	64.98
S₁	33.35	2.29	1.78	0.107	1.59	61.64
S₂	123.00	3.09	2.76	0.177	1.42	115.70
S₃	82.49	2.51	2.36	0.116	1.33	70.92

下载: 导出CSV

表 7 不同处理间差异代谢物的相对丰度

Table 7 Relative abundance of differential metabolites among treatments

代谢物类别 class of metabolites	化合物名称 name of compounds	处理 treatments
代谢物类别 class of metabolites	化合物名称 name of compounds	CK	S₁	S₂	S₃
植物甾醇及其衍生物 phytosterols and their derivatives	麦角固醇 ergosterol	1.000±0.027	0.866±0.057 ↓	3.339±0.070 ↑	3.081±0.041 ↑
	环阿屯醇 cycloartenol	1.000±0.017	0.953±0.030 ↓	4.417±0.045 ↑	4.310±0.028 ↑
	胡萝卜甾醇 daucosterol	1.000±0.176	0.948±0.180 ↓	4.651±0.168 ↑	4.291±0.048 ↑
	丁香树脂醇 (+)-syringaresinol	1.000±0.022	1.110±0.010 ↑	3.269±0.068 ↑	3.319±0.055 ↑
酚类化合物 phenolic compounds	黄氧素 xanthoxin	1.000±0.017	1.089±0.031 ↑	0.898±0.052 ↓	0.978±0.034
	木香烃内酯 (+)-costunolide	1.000±0.011	1.053±0.011 ↑	1.012±0.026	1.054±0.018 ↑
	金丝桃苷 hyperoside	1.000±0.036	0.970±0.021	1.007±0.032	0.976±0.027 ↓
	玉米黄质 zeaxanthin	1.000±0.013	1.089±0.027 ↑	0.983±0.038	1.054±0.018 ↑
	叶黄素 lutein	1.000±0.016	1.057±0.040 ↑	0.853±0.046 ↓	0.921±0.041 ↓
	环巴胺 cyclopamine	1.000±0.006	1.068±0.029 ↑	0.977±0.031	1.042±0.031 ↑
	眉藻黄素 caloxanthin	1.000±0.008	1.037±0.015 ↑	0.978±0.011 ↓	1.020±0.012 ↑
	苦木素 quassin	1.000±0.024	1.026±0.018	0.960±0.013 ↓	0.957±0.019 ↓
	7-羟基-8-甲氧基香豆素 7-hydroxy-8-methoxycoumarin	1.000±0.010	1.132±0.038 ↑	1.703±0.031 ↑	1.132±0.052 ↑
氨基酸类化合物 amino acid compounds	缬氨酸 valine	1.000±0.008	1.047±0.022 ↑	0.959±0.022 ↓	1.003±0.028
	水合酪氨酸 hydroxytyrosol	1.000±0.006	1.068±0.029 ↑	0.977±0.031	1.042±0.031 ↑
	亮氨酸 L-(+)-leucine	1.000±0.009	1.145±0.070 ↑	1.043±0.030 ↑	1.222±0.068 ↑
	苯基丙氨酸 D-(+)-phenylalanine	1.000±0.008	1.071±0.023 ↑	1.004±0.029	1.039±0.061
	色氨酸 D-(+)-tryptophan	1.000±0.008	1.077±0.017 ↑	1.023±0.031 ↑	1.023±0.050 ↑
	乙酰基精氨酸 acetylarginine	1.000±0.081	0.977±0.048 ↓	0.892±0.050 ↓	1.072±0.051
	甲基酪氨酸 o-methyltyrosine	1.000±0.014	0.906±0.036 ↓	0.902±0.055 ↓	0.895±0.057 ↓
醛类化合物 aldehydes	芸香糖 rutinose (beta-anomer)	1.000±0.012	1.080±0.056 ↑	0.997±0.039	1.110±0.033 ↑
	对叔丁基苯丙醛 bourgeonal	1.000±0.008	1.037±0.015 ↑	0.978±0.011 ↓	1.020±0.012 ↑
	香草素 vanillin	1.000±0.020	1.079±0.034 ↑	1.026±0.026	1.049±0.039
其他 other	伞形酸 umbellic acid	1.000±0.008	1.078±0.043 ↑	0.979±0.040 ↓	1.053±0.032 ↑
	己糖酸 hexonic acid	1.000±0.007	1.105±0.032 ↑	1.131±0.036 ↑	1.103±0.049 ↑
	3-乙氧基苯甲酰胺 3-ethoxybenzamide	1.000±0.115	2.092±0.057 ↑	4.106±0.037 ↑	3.853±0.050 ↑
	乙氧酰胺苯甲酯 ethopabate	1.000±0.031	1.117±0.063 ↑	1.002±0.030	1.174±0.048 ↑
	泛酸 pantothenic acid	1.000±0.024	1.049±0.025 ↑	1.064±0.019 ↑	1.039±0.020
	乙酰胺 acetamide	1.000±0.044	1.264±0.029 ↑	1.508±0.043 ↑	1.407±0.046 ↑
注： ↑ 表示显著上调， ↓ 表示显著下调。 Note: ↑ means significant increase, ↓ means significant decrease.

下载: 导出CSV

参考文献(26)

[1]	LIU A M, YUAN K L, LI Q, et al. Metabolomics and proteomics revealed the synthesis difference of aroma precursors in tobacco leaves at various growth stages[J]. Plant Physiology and Biochemistry, 2022, 192: 308. DOI: 10.1016/j.plaphy.2022.10.016.
[2]	王忠, 宋春燕. 关于对烟草种植技术的研究与探讨[J]. 农业与技术, 2015, 35(12): 129. DOI: 10.11974/nyyjs.20150633116.
[3]	韦建玉, 金亚波, 吴峰, 等. 烤烟品种K326、云烟85及云烟87的适应性研究[J]. 安徽农业科学, 2008, 36(6): 2362. DOI: 10.13989/j.cnki.0517-6611.2008.06.096.
[4]	曾昭松, 吴才源, 龙立汪, 等. 烤烟品种云烟87生产技术研究进展[J]. 安徽农学通报, 2018, 24(9): 38. DOI: 10.16377/j.cnki.issn1007-7731.2018.09.017.
[5]	黎妍妍, 李锡宏, 王林, 等. 富硒烟叶对烟气有害成分释放量作用分析[J]. 中国烟草学报, 2013, 19(2): 7. DOI: 10.3969/j.issn.1004-5708.2013.02.002.
[6]	韩丹. 硒在烤烟中的累积、形态转化及缓解砷毒害的机理研究[D]. 武汉: 华中农业大学, 2015.
[7]	韩丹, 熊双莲, 许自成, 等. 硒对烤烟生长、化学指标及矿质营养元素含量的影响[J]. 核农学报, 2017, 31(10): 2072. DOI: 10.11869/j.issn.100-8551.2017.10.2072.
[8]	杨应粉, 詹军, 罗华元, 等. 增施硒肥对连作地块烤烟生长发育及内在品质的影响[J]. 甘肃农业大学学报, 2016, 51(3): 83. DOI: 10.13432/j.cnki.jgsau.2016.03.014.
[9]	刘晓迪, 景延秋, 张豹林, 等. 钼、硒营养对烤烟化学成分和香气物质的影响[J]. 安徽农业大学学报, 2016, 43(4): 610. DOI: 10.13610/j.cnki.1672-352x.20160712.031.
[10]	GARCIA-SEVILLANO M A, GARCIA-BARRERA T, NAVARRO-ROLDAN F, et al. A combination of metallomics and metabolomics studies to evaluate the effects of metal interactions in mammals. Application to Mus musculus mice under arsenic/cadmium exposure[J]. Journal of Proteomics, 2014, 104: 66. DOI: 10.1016/j.jprot.2014.02.011.
[11]	SZPUNAR J. Metallomics: a new frontier in analytical chemistry[J]. Analytical and Bioanalytical Chemistry, 2004, 378(1): 54. DOI: 10.1007/s00216-003-2333-z.
[12]	FIEHN O. Metabolomics-the link between genotypes and phenotypes[J]. Plant Molecular Biology, 2002, 48(1/2): 155. DOI: 10.1023/A:1013713905833.
[13]	NICHOLSON J K, LINDON J C, HOLMES E. ‘Metabonomics’: understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data[J]. Xenobiotica, 1999, 29(11): 1181. DOI: 10.1080/004982599238047.
[14]	MISHRA P, KUMAR R, DWIVEDI A, et al. Analysis of constituents present in smokeless tobacco (Nicotiana tabacum) using spectroscopic techniques[J]. Methods and Applications in Fluorescence, 2022, 10(3): 1. DOI: 10.1088/2050-6120/ac5e11.
[15]	HU Z R, YAN W J, YANG C K, et al. Integrative analysis of transcriptome and metabolome provides insights into the underlying mechanism of cold stress response and recovery in two tobacco cultivars[J]. Environmental and Experimental Botany, 2022, 200: 104920. DOI: 10.1016/j.envexpbot.2022.104920.
[16]	梁珂馨, 奚家勤, 陈若星, 等. 基于非靶向代谢组学的片烟醇化进程判断及代谢物差异分析[J]. 烟草科技, 2023, 56(4): 17. DOI: 10.16135/j.issn1002-0861.2022.0582.
[17]	WANG Z J, WU J, SUN Z P, et al. ICP-MS based metallomics and GC-MS based metabolomics reveals the physiological and metabolic responses of Dendrobium huoshanense plants exposed to Fe₃O₄ nanoparticles[J]. Frontiers in Nutrition, 2022, 9: 1013756. DOI: 10.3389/fnut.2022.1013756.
[18]	谢瑞, 易忠经, 申修贤, 等. 不同磷肥施用量对植烟土壤磷有效性及烟叶产量和品质的影响[J]. 山地农业生物学报, 2017, 36(6): 38. DOI: 10.15958/j.cnki.sdnyswxb.2017.06.007.
[19]	蔡海洋, 李玉梅, 熊德中. 不同硫肥用量对烤烟光合作用和产量的影响[J]. 亚热带农业研究, 2009, 5(3): 155. DOI: 10.13321/j.cnki.subtrop.agric.res.2009.03.010.
[20]	马友华, 丁瑞兴, 张继榛, 等. 硒硫相互作用对烤烟(Nicotina tabacum L.)吸收硒的影响[J]. 南京农业大学学报, 2001, 24(1): 55. DOI: 10.7685/j.issn.1000-2030.2004.04.012.
[21]	BLOEM E, HANEKLAUS S H, SCHNUG E. Influence of nitrogen and sulfur fertilization on the alliin content of onions and garlic[J]. Journal of Plant Nutrition, 2004, 27(10): 1827. DOI: 10.1081/LPLA-200026433.
[22]	JIANG C Q, ZU C L, SHEN J, et al. Effects of selenium on the growth and photosynthetic characteristics of flue-cured tobacco (Nicotiana tabacum L.)[J]. Acta Societatis Botanicorum Poloniae, 2015, 84(1): 71. DOI: 10.5586/asbp.2015.006.
[23]	FAN J, WAN H, WANG R Y, et al. Growth and physiological characteristics of flue-cured tobacco (Nicotiana tabacum L.) under selenate-Se stress[J]. Communications in Soil Science and Plant Analysis, 2018, 49(18): 2238. DOI: 10.1080/00103624.2018.1499759.
[24]	LIU Y, SANZ-SÁEZ A, MENG X Y, et al. Promoting effect of nano-Se on tobacco growth and reactive oxygen species metabolism[J]. Journal of Plant Nutrition, 2020, 43(18): 2810. DOI: 10.1080/01904167.2020.1793187.
[25]	LOCHMAN J, MIKES V. Ergosterol treatment leads to the expression of a specific set of defence-related genes in tobacco[J]. Plant Molecular Biology, 2006, 62(1/2): 43. DOI: 10.1007/s11103-006-9002-5.
[26]	OUYANG M A, WEIN Y S, ZHANG Z K, et al. Inhibitory activity against tobacco mosaic virus (TMV) replication of pinoresinol and syringaresinol lignans and their glycosides from the root of Rhus javanica var. roxburghiana[J]. Journal of Agricultural and Food Chemistry, 2007, 55(16): 6460. DOI: 10.1021/jf0709808.

施引文献(8)

期刊类型引用(4)

1.	王军川，逄红，王瑞霞，马国栋，靳磊，张萍. 基于GEE和Canoco的沙丘移动速率及其内外控制因子贡献率分析. 中国水土保持科学(中英文). 2024(03): 93-99 . 百度学术
2.	陈嘉嘉，黄磊，贾鸿飞，杨利贞，杨贵森，姚军，张艳艳. 腾格里沙漠人工植被区不同植被类型土壤水分垂直变异特征及其动态变化模型. 干旱区资源与环境. 2022(04): 126-133 . 百度学术
3.	刘秋丽. 霍泉灌区小尺度农田土壤含水率时变性研究. 水资源开发与管理. 2021(10): 6-11 . 百度学术
4.	陈晓莹，陈林，宋乃平，李敏岚，吴婷，杨新国，王磊，王兴. 荒漠草原两种类型土壤的水分动态对比. 应用生态学报. 2020(05): 1560-1570 . 百度学术

其他类型引用(4)

资源附件(0)

图(5) / 表(7)

计量

文章访问数: 1528
PDF下载量: 13
被引次数: 8

1. 材料与方法
1.1 研究区概况
1.2 研究方法
1.3 数据处理
2. 结果与分析
2.1 土壤含水量均值及变异系数
2.2 土壤含水量空间分布特征
3. 讨论
4. 结论

利用离子组学和代谢组学研究硒肥对烤烟元素及化学成分的影响

作者简介: 范晓苏(1988—)，女，广西来宾人，硕士，实验师，主要从事质谱和光谱仪器分析研究。E-mail：fanxs@gxu.edu.cn

通信作者: 沈方科(1973—)，男，广西博白人，硕士，高级农艺师，主要从事植物营养生理与施肥技术研究。E-mail：shenfk@gxu.edu.cn

计量

出版历程