作物杂志, 2025, 41(2): 196-206 doi: 10.16035/j.issn.1001-7283.2025.02.027

生理生化·植物营养·栽培耕作

叶面施硒对燕麦籽粒营养品质的影响

张丽娜,1, 杨文平2, 苏苗1, 张之玄1, 李军辉1, 陈杰1, 高志强1, 杨珍平,1

1山西农业大学农学院/黄土高原特色作物优质高效生产省部共建协同创新中心,030801,山西晋中

2华北理工大学生命科学学院,063210,河北唐山

Effect of Foliar Application of Selenium onNutritional Quality of Oat Grain

Zhang Lina,1, Yang Wenping2, Su Miao1, Zhang Zhixuan1, Li Junhui1, Aamir Ali1, Chen Jie1, Gao Zhiqiang1, Yang Zhenping,1

1College of Agriculture, Shanxi Agricultural University / Collaborative Innovation Center for High-Quality and Efficient Production of Characteristic Crops on the Loess Plateau, Jointly Built by the Province and the Ministry, Jinzhong 030801, Shanxi, China

2School of Life Sciences, North China University of Technology, Tangshan 063210, Hebei, China

通讯作者: 杨珍平,主要从事旱作栽培与作物生态及天然植物功能成分挖掘与利用研究,E-mail:yangzp.2@163.com

收稿日期: 2023-12-11   修回日期: 2024-01-30  

基金资助: 山西省回国留学人员科研资助项目(2023-093)
黄土高原特色作物优质高效生产省部共建协同创新中心(SBG JXTZX-21)
山西农谷建设科研专项项目(SXNGJSKYZX201701)
山西省现代农业杂粮产业技术体系(2024CYJSTX03-03)

Received: 2023-12-11   Revised: 2024-01-30  

作者简介 About authors

张丽娜,主要从事旱作栽培与作物生态及天然植物功能成分挖掘与利用研究,E-mail:15369317536@163.com

摘要

为探究叶面施硒(Se)对燕麦籽粒硒积累和营养品质的作用,以坝莜6号、坝莜12号、白燕2号3个燕麦品种为研究对象设置田间试验,在扬花期叶面喷施不同浓度(Se1:12 g/hm2;Se2:36 g/hm2)的硒营养液(SeO32-),同时以不喷硒为对照,测定灌浆期籽粒硒含量以及成熟期营养品质(氨基酸和脂肪酸)。结果表明,随灌浆时间延长,燕麦籽粒硒含量呈先升高后降低的趋势。其中,成熟期施硒各处理籽粒硒含量均在0.10~0.30 mg/kg的国家标准范围内,显著高于对照,且在Se2处理下籽粒硒含量最高,为0.16 mg/kg。与对照相比,施硒明显提高了3个燕麦品种籽粒中总氨基酸的含量,低氨基酸品种坝莜6号增加了2.13%~10.53%,高氨基酸品种白燕2号增加了7.44%~9.37%。与对照相比,施硒提高了3个燕麦品种籽粒中不同脂肪酸含量,且品种类型、硒肥浓度及二者的互作效应的影响均达到极显著水平(P<0.01)。结果说明叶面喷施硒肥能显著促进燕麦籽粒硒积累,强化籽粒氨基酸和不饱和脂肪酸的含量,其中以Se2处理效果更好。

关键词: 燕麦; 硒肥; 营养品质; 氨基酸; 脂肪酸

Abstract

In order to explore the effects of foliar selenium (Se) application on selenium accumulation and nutritional quality in oat grains, a field experiment was conducted using three oat varieties (Bayou 6, Bayou 12, and Baiyan 2). At the flowering stage, selenium nutrient solution (SeO32-) with different concentrations (Se1: 12 g/ha, Se2: 36 g/ha) were sprayed on the leaves, and selenium-free was used as control (CK). The selenium content in grains at grain-filling stage and the nutritional quality (amino acids and fatty acids) at maturity stage were determined. The findings demonstrated that the selenium concentration in oat grains first increased and then decreased with the extension of the filling period. Specifically, at maturity stage, selenium contents of all selenium-treated treatments was within the national standard range of 0.10-0.30 mg/kg, which was significantly higher than that of the CK. The highest selenium content (0.16 mg/kg) was observed under the Se2 treatment. Compared with the CK, selenium application significantly increased the total amino acid content in oat grains of all three varieties: the content in the low amino acid variety ʻBayou 6ʼ increased by 2.13%-10.53%, while the content in the high amino acid variety ʻBaiyan 2ʼ increased by 7.44%-9.37%. Compared with the CK, selenium application increased the content of different fatty acids in grains of the three oat varieties. The influences of variety type, selenium fertilizer concentration, and their interaction effect all reached a very significant level (P < 0.01). Foliar selenium spraying can significantly promote the accumulation of selenium and enhance the contents of amino acids and unsaturated fatty acids in oat grains, the effect of the Se2 treatment is better.

Keywords: Oat; Selenium fertilizer; Nutritional quality; Amino acid; Fatty acid

PDF (1552KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

张丽娜, 杨文平, 苏苗, 张之玄, 李军辉, 陈杰, 高志强, 杨珍平. 叶面施硒对燕麦籽粒营养品质的影响. 作物杂志, 2025, 41(2): 196-206 doi:10.16035/j.issn.1001-7283.2025.02.027

Zhang Lina, Yang Wenping, Su Miao, Zhang Zhixuan, Li Junhui, Aamir Ali, Chen Jie, Gao Zhiqiang, Yang Zhenping. Effect of Foliar Application of Selenium onNutritional Quality of Oat Grain. Crops, 2025, 41(2): 196-206 doi:10.16035/j.issn.1001-7283.2025.02.027

硒(Se)是人体必需微量元素,具有提高免疫力、增强机体新陈代谢等积极作用[1-3]。但土壤硒元素地理分布不均匀,我国大约72%的地区土壤存在不同程度的缺硒或低硒情况,由此大部分粮食作物硒含量也较低,不能满足人体日常硒摄入量[4-5]。燕麦(Avena sativa)被誉为“全营养谷类作物”[6-8],不仅富含膳食纤维,而且具有良好的氨基酸构成、丰富的不饱和脂肪酸和大量抗氧化活性物质,但其对硒的积累能力在禾谷类作物中相对较弱[9],因此,强化燕麦籽粒硒积累,明确其对氨基酸构成以及不饱和脂肪酸组成的影响,对助推营养健康型农业发展具有重要意义[10]

外源施加硒肥是提高粮食作物硒含量的一种主要方式,常用的2种途径是叶面喷施和土壤施用[11-12],其中叶面喷施是最直接有效的途径[13],其富集效果与作物品种、施硒时期、施硒量密切相关。研究[14]表明,对14个燕麦品种进行叶面施硒处理后,不同品种的硒富集能力存在差异。其中,白燕2号硒富集效果最好(0.30 mg/kg2),晋燕8号和晋燕18号富集效果最差(0.20 mg/kg2)。但在其他研究[15]结果中,硒富集效果最好的品种则是白燕7号。关于喷施时期及施用量的研究[14-18]发现,在燕麦抽穗期喷施Na2SeO3籽粒硒富集效果最佳,但最佳施用量还存在争议,从21.92到70.00 g/hm2不等。在叶面施硒对作物籽粒营养品质影响的研究中发现,硒处理提高了燕麦蛋白质(3.29%~18.62%)、β-葡聚糖(50.37%~186.54%)[14]和粗蛋白的含量(18.90%)[16],降低了燕麦籽粒中铜、铁、锌、镁含量,提高了锰含量,但对钙和硫没有影响[17]。此外,叶面喷硒同样有助于提高小麦营养品质,增加其抗氧化能力[19],促进水稻的食用和糊化品质[20],增加谷子粗蛋白(13.90%)、赖氨酸(17.90%)和叶酸(7.50%)的含量[21]。目前,关于燕麦硒强化的研究[14-18]主要集中在籽粒硒含量、产量、蛋白质及其组分和矿质元素等方面的影响,而对灌浆期籽粒硒积累规律,燕麦氨基酸组成和脂肪酸成分强化方面的研究较少。

本研究以优质专用型燕麦品种坝莜12号、坝莜6号和白燕2号为研究对象,在山西省大同市左云县进行种植,于燕麦扬花期喷施不同浓度叶面硒营养液(SeO32-形式),通过测定灌浆期燕麦籽粒硒积累量和成熟期籽粒氨基酸和不饱和脂肪酸组成及含量,分析不同浓度硒处理对燕麦籽粒营养品质是否具有强化作用,为高寒生态区富硒燕麦优质生产提供理论和实践依据。

1 材料与方法

1.1 试验地概况

试验于2021年在山西省大同市左云县安烟墩村(112°44′22″ E,40°8′33″ N)进行。该地区位于山西省最北部,属于晋西北高寒生态区,平均海拔1200.00 m以上,是典型的温带半干旱大陆性季风气候,年均气温6.11 ℃,年均日照时数2763.80 h,年均降水量399.00 mm,年均无霜期125 d。土壤类型为淡栗褐土,播前0~20 cm土层土壤有机质、碱解氮、速效磷、速效钾和全氮含量分别为11.04 g/kg、94.40 mg/kg、13.80 mg/kg、63.00 mg/kg和0.92 g/kg,硒含量为0.08 mg/kg,pH 8.35。燕麦生长季为每年4-10月,生长季平均气温14.61 ℃,降水量374.90 mm(图1)。

图1

图1   2021年4-10月降水量和平均气温

Fig.1   Rainfall and average temperature from April to October in 2021


1.2 试验材料

供试燕麦品种为坝莜12号(BY12)、坝莜6号(BY6)和白燕2号(BY2),由山西农业大学高寒区作物研究所与国家燕麦荞麦产业技术体系右玉县农业试验站提供。坝莜12号具有低纤维、易脱粒的特点,主要是燕麦米专用加工品种,为试验当地主栽品种。坝莜6号具有生育期短、低脂肪、面粉保存时间长、出粉率高、面粉色泽白的特点,主要是救灾备荒及大众化磨粉专用加工品种。白燕2号具有高蛋白和高脂肪等特点,主要是大粒型燕麦片专用加工品种。

硒肥为叶面硒营养液(含SeO32- 8.00 mg/mL),由山西省太原市伊丰富硒农业产品研发有限公司提供。

1.3 试验处理

试验采用两因素裂区设计。主区为3个燕麦品种,副区为3个施硒浓度:0(CK)、12.00(Se1)和36.00 g/hm2(Se2)。各处理3个小区(3.00 m×6.00 m)重复,共27个小区。各小区间隔0.50 m,防止喷雾漂移造成交叉污染。

播前整地,机械化条播,行距0.25 m,播种量120.00 kg/hm2。在燕麦扬花期喷施硒肥,按照试验设计,各小区分别取硒原液0、1.44、4.32 mL,对水量450.00 L/hm2,对照组喷相同体积蒸馏水。其余田间管理措施按常规操作进行。喷施硒肥7 d后(7月23日)开始每7 d取样1次,直至成熟期(9月3日)结束,共采样6次。每次取样时各小区剪10穗装入自封袋,标明处理和取样时间,带回实验室。

1.4 指标测定与方法

1.4.1 籽粒硒含量

将每次采集的燕麦穗用大量去离子水冲洗干净后,于烘箱105 ℃杀青30 min,然后调至80 ℃烘干至恒重。脱粒后用电动研磨机(HCP100,中国浙江)磨成面粉,过0.15 mm筛后,参照GB 5009.93-2017采用原子荧光光谱法测定籽粒硒含量。

1.4.2 成熟期籽粒营养品质

采用上述1.4.1中研磨过筛的燕麦籽粒粉末进行下述指标测定,参照GB/T 5009采用酸水解法测定16种氨基酸含量。

按照以下公式计算[22]氨基酸相关指标值:

氨基酸指数(ratio of amino acid,RAA)=某必需氨基酸含量/标准模式(FAO/WHO)中相对应的氨基酸含量;

氨基酸比值系数(ratio coefficient of amino acid,RC)=某必需氨基酸RAA/各种氨基酸RAA均值;

氨基酸比值系数分(score ratio coefficient amino acid,SRC)=100-CV×100;

必需氨基酸指数(essential amino acid index,EAAI)=na1b1×a2b2×a3b3××anbn

式中,n为比较的必需氨基酸(EAA)的数

目,n≤8;a1a2a3,…,an为试验中EAA的含量[g/(100 g)干物质];b1b2b3,…,bn为标准蛋白质[FAO/WHO评分模式]中与an相应的必需氨基酸的含量[g/(100 g)干物质]。

采用脂肪酸甲酯化法进行提取脂肪酸,再采用气相色谱仪测定脂肪酸含量[23]

1.5 数据处理

采用Excel 2010软件对数据进行处理,用SPSS 26软件进行方差分析、多重比较、相关性及通径分析,用Origin 2022软件进行作图,用R studio进行相关性绘图。

2 结果与分析

2.1 叶面施硒对燕麦灌浆期籽粒硒积累动态的影响

图2表明,叶面喷硒显著提高了各燕麦品种籽粒硒含量,且随硒肥浓度增加,提高的幅度越明显。成熟期各品种的籽粒硒含量均在0.10~0.30 mg/kg标准范围(GB/T 22499-2008),其中白燕2号籽粒总硒积累量最多(0.17 mg/kg),坝莜6号中最少(0.15 mg/kg)。

图2

图2   叶面喷硒下燕麦灌浆期籽粒硒含量积累的动态

Ov代表燕麦品种,Se代表硒肥处理,Ov×Se代表燕麦品种和硒肥处理互作。不同小写字母表示在同一品种不同硒肥浓度间差异显著(P < 0.05);大写字母不同表示在同一硒肥浓度不同品种间差异显著(P < 0.05)。“*”、“**”、“***”和ns分别表示P < 0.05、P < 0.01、P < 0.001和P > 0.05。下同。

Fig.2   The dynamics of selenium accumulation in grains of oat at filling stage under foliar selenium spraying

Ov means oat variety, Se means foliar selenium treatment, Ov×Se signify interaction between oat varieties and foliar selenium. Different lowercase letters indicate significant difference at P < 0.05 level among different treatments under the same variety; Different capital letters indicate that the difference between different varieties under the same treatment at P < 0.05 level.“*”,“**”,“***”and“ns”indicate P < 0.05, P < 0.01, P < 0.001 and P > 0.05, respectively. The same below.


花后籽粒中硒积累呈现先增后减的趋势,但同一浓度下不同品种间变化幅度不同,各时期出现的时间点也不同。相对于坝莜12号来说(峰值出现期在花后28 d),坝莜6号和白燕2号需要较长的积累时期(峰值出现期在花后35 d),但高硒处理会缩短籽粒硒积累的时间,使籽粒提前进入硒饱和阶段,具体表现为坝莜12号积累峰值出现在花后21 d,坝莜6号和白燕2号峰值则出现在花后28 d,均缩短了7 d。

同一品种对于不同浓度的叶面硒响应也不同。相对于CK处理来说,低硒品种坝莜6号对外源硒处理更为敏感,使籽粒在形成初期快速累积硒元素,其含量可提高67.37%~125.40%,大大缩短籽粒硒积累时间;而高硒品种白燕2号对外源硒处理敏感程度低,快速积累时期(花后7~14 d)籽粒硒积累幅度小,为9.33%~15.25%,但在积累起始期(花后7 d),其含量在3个品种中最高(图2)。

不同品种对于硒响应的阈值也不同。当地主推品种坝莜12号阈值最低,呈现出明显的低促进高抑制的趋势,CK处理下,快速积累时期籽粒硒积累幅度最大(66.67%),Se2处理最小,为9.80%;坝莜6号虽然品种原始籽粒硒含量低,但它的阈值较坝莜12号高,与CK处理相比,处理Se1和Se2籽粒硒积累幅度分别为125.40%和67.37%;白燕2号在高浓度处理(Se2)时,籽粒硒积累幅度为12.57%。

2.2 叶面施硒对燕麦成熟期籽粒氨基酸组成及含量的影响

根据所测燕麦籽粒的16种氨基酸含量(表1)可知,3个燕麦品种中白燕2号的总氨基酸含量及大多数氨基酸含量最高,其次是坝莜12号和坝莜6号。施硒明显提高了各品种的总氨基酸含量。与CK相较,低氨基酸品种坝莜6号的总氨基酸含量增加2.13%~10.53%;高氨基酸品种白燕2号总氨基酸含量增加7.44%~9.37%,而坝莜12号较对照增加1.45%~2.08%。

表1   叶面喷硒下燕麦成熟期籽粒总氨基酸和必需氨基酸含量

Table 1  Essential and total amino acid contents in grains of oat at maturity stage under foliar selenium spraying mg/g

氨基酸
Amino acid
品种
Variety
硒肥浓度Se concentrationFF-value
CKSe1Se2均值Mean品种Ov硒肥Se品种×硒肥Ov×Se
总氨基酸TAA坝莜6号109.16cC111.48bC120.66aB113.77*********
坝莜12号115.20bB116.87aB117.60aC116.56
白燕2号122.35cA131.45bA133.81aA129.20
异亮氨酸Ile坝莜6号13.46bA13.49bA15.35aA14.10****ns
坝莜12号12.88aA13.50aA13.74aB13.37
白燕2号14.06aA14.72aA15.15aA14.64
甲硫氨酸Met坝莜6号11.36cA11.79bA12.73aA11.96ns*ns
坝莜12号12.11aA12.55aA13.18aA12.61
白燕2号11.36aA11.98aA12.66aA12.00
亮氨酸Leu坝莜6号9.87bC9.95bB11.63aB10.48*********
坝莜12号10.50aB10.57aB10.79aC10.62
白燕2号11.39bA12.78aA13.50aA12.56
缬氨酸Val坝莜6号9.78aA9.92aA10.44aA10.05*nsns
坝莜12号9.98aA10.00aA10.17aA10.05
白燕2号10.12aA10.57aA11.43aA10.71
苯丙氨酸Phe坝莜6号4.73cB5.20bC5.58aB5.17******ns
坝莜12号5.70bA6.09aB6.25aA6.01
白燕2号5.91bA6.40aA6.56aA6.29
赖氨酸Lys坝莜6号3.94bB4.15abC4.39aB4.16******ns
坝莜12号4.29bA4.55aB4.58aB4.47
白燕2号4.60bA4.92aA5.01aA4.84
苏氨酸Thr坝莜6号3.43aB3.45aB3.61aB3.50******ns
坝莜12号3.73aA3.95aA3.99aAB3.89
白燕2号3.66bA4.24aA4.34aA4.08

不同小写字母表示在同一品种不同硒肥浓度间差异显著(P < 0.05);不同大写字母表示在同一硒肥浓度不同品种间差异显著(P < 0.05)。“*”、“**”、“***”和“ns”分别表示P < 0.05、P < 0.01、P < 0.001和P > 0.05。下同。

Different lowercase letters indicate significant difference at P < 0.05 level among different treatments under the same variety; Different capital letters indicate significant difference between different varieties at P < 0.05 level under the same treatment.“*”,“**”,“***”and“ns”indicate P < 0.05, P < 0.01, P < 0.001 and P > 0.05, respectively. The same below.

新窗口打开| 下载CSV


在必需氨基酸(色氨酸除外)中,异亮氨酸含量最高,其次是甲硫氨酸、亮氨酸、缬氨酸,这4种氨基酸含量平均达到10.00 mg/g以上,但各氨基酸丰度在不同品种中表现不同。除甲硫氨酸(坝莜12号最高)外,白燕2号中3种氨基酸含量最高。施硒对不同品种籽粒中缬氨酸含量没有明显的促进作用,但对苯丙氨酸和赖氨酸的含量有明显的促进作用。不同品种中必需氨基酸含量对硒响应不同,与低浓度处理相比,坝莜6号中甲硫氨酸含量显著提高,异亮氨酸和亮氨酸含量对高浓度硒响应更为显著,分别提高13.83%和16.86%;坝莜12号中则是赖氨酸对硒浓度响应明显;白燕2号中赖氨酸、苏氨酸以及亮氨酸含量显著提高,分别提高7.13%~8.99%、15.95%~18.72%和12.21%~ 18.51%。

在非必需氨基酸中(表2),谷氨酸含量最高,可达15.00 mg/g以上,白燕2号含量最高。硒处理显著提高了坝莜6号中谷氨酸、天冬氨酸含量和白燕2号中丝氨酸含量;与CK相比,高浓度硒处理下白燕2号中酪氨酸、甘氨酸、谷氨酸、丙氨酸、天冬氨酸、胱氨酸以及组氨酸7种非必需氨基酸含量显著提高,分别增加10.96%、15.13%、9.40%、11.00%、5.49%、12.20%和9.25%;坝莜6号酪氨酸、甘氨酸、丙氨酸、组氨酸和精氨酸含量显著提高,而在坝莜12号中仅有酪氨酸、精氨酸和天冬氨酸3种氨基酸显著提高,分别较CK增加14.13%、32.35%和6.46%。其他氨基酸含量随施硒量的变化因品种而异。总体来看,施硒提高了坝莜6号12种氨基酸含量、坝莜12号5种氨基酸含量和白燕2号12种氨基酸含量。

表2   叶面硒肥下燕麦成熟期籽粒非必需氨基酸含量

Table 2  Non-essential amino acid contents in grains of oat at maturity stage under foliar selenium spraying mg/g

氨基酸
Amino acid
品种
Varietiy
硒肥浓度Se concentrationFF-value
CKSe1Se2均值Mean品种Ov硒肥Se品种×硒肥Ov×Se
谷氨酸Glu坝莜6号16.20cB17.71bC18.79aC17.57*********
坝莜12号20.03aA20.30aB20.51aB20.28
白燕2号20.20bA22.03aA22.10aA21.44
天冬氨酸Asp坝莜6号7.68cC8.04bC8.42aC8.05******ns
坝莜12号8.36bB8.66abB8.90aB8.64
白燕2号9.16bA9.26bA9.67aA9.36
精氨酸Arg坝莜6号6.05bB6.11bB6.87aB6.34*****ns
坝莜12号3.06bC3.23bC4.05aC3.45
白燕2号7.69aA7.90aA8.09aA7.89
酪氨酸Tyr坝莜6号4.18bC4.36bB4.98aA4.51******ns
坝莜12号4.46bB4.61bB5.09aA4.72
白燕2号5.02bA5.24abA5.57aA5.28
甘氨酸Gly坝莜6号4.18bB4.25bC4.53aC4.32*********
坝莜12号4.64aA4.68aB4.76aB4.69
白燕2号4.56bA5.12aA5.25aA4.98
丝氨酸Ser坝莜6号4.21aB4.31aB4.35aC4.29*********
坝莜12号4.67aA4.79aA4.82aB4.76
白燕2号4.17cB4.96bA5.27aA4.80
丙氨酸Ala坝莜6号3.56bC3.73abB3.91aB3.73******ns
坝莜12号3.85aB3.98aB4.08aB3.97
白燕2号4.19bA4.65aA4.65aA4.50
胱氨酸Cys坝莜6号2.68aA2.88aA2.99aA2.85ns**ns
坝莜12号2.53aA2.59aB2.85aA2.66
白燕2号2.75bA2.80bAB3.08aA2.88
组氨酸His坝莜6号2.44bB2.80abB2.87aB2.70******ns
坝莜12号2.61aAB2.61aC2.64aC2.62
白燕2号2.82bA2.97aA3.08aA2.96

新窗口打开| 下载CSV


燕麦中必需氨基酸分别与FAO/WHO氨基酸模式谱和理想蛋白质标准比较(表3),异亮氨酸、甲硫氨酸+胱氨酸含量与模式蛋白组成较一致,而亮氨酸、苯丙氨酸+酪氨酸、酪氨酸和缬氨酸含量不足。E/T与E/N值明显高于FAO/WHO氨基酸模式值,均高于理想蛋白标准,说明燕麦具有较高食用价值。对于同一品种来说,高浓度处理Se2可强化燕麦籽粒氨基酸含量,提高其食用价值。对于不同品种来说,喷硒对白燕2号营养强化效果显著。

表3   FAO/WHO氨基酸模型和理想蛋白质标准下不同处理燕麦籽粒必需氨基酸含量比较

Table 3  Comparison of contents of essential amino acids of oats in different treatments with the FAO/WHO amino acid model and the ideal protein standard

品种
Variety
处理
Treatment
氨基酸类型Amino acid type (mg/g)理想蛋白标准Ideal protein standard
IleLeuLysMet+CysPhe+TyrThrValE/TE/N
FAO/WHO标准模式4.007.005.503.506.004.005.004060
坝莜6号CK3.37bA1.41bC0.72bB4.01cA1.49cC0.86aA1.96aA52aA111aA
Se13.37bA1.42bB0.75abC4.19bA1.59bC0.86aB1.98aA52aA107aB
Se23.84aA1.66aB0.80aB4.49aA1.76aB0.90aB2.09aA53aA110aA
坝莜12号CK3.22aA1.50aB0.78bA4.18aA1.69cB0.93aA2.00aA51aAB109aA
Se13.38aA1.51aB0.83aB4.33aA1.78bB0.99aA2.00aA52aA109aA
Se23.43aB1.54aC0.83aB4.58aA1.89aAB1.00aAB2.03aA53aA110aA
白燕2号CK3.52aA1.63bA0.84bA4.03bA1.82bA0.91bA2.02aA50aB101aB
Se13.68aA1.83aA0.90aA4.22abA1.94abA1.06aA2.11aA50aB101aC
Se23.79aA1.93aA0.91aA4.50aA2.02aA1.09aA2.29aA51aA103aB

新窗口打开| 下载CSV


为进一步说明不同燕麦品种氨基酸的营养价值,通过RC、SRC和EAAI指标进行了营养价值评价。各材料7种必需氨基酸的RC、SRC和EAAI值见表4

表4   不同处理下燕麦必需氨基酸RC、SRC和EAAI的比较

Table 4  Comparison of RC, SRC, and EAAI of essential amino acids in oat varieties under different treatments

品种
Variety
处理
Treatment
RCSRCEAAI
LysThrLeuValIleMet+CysPhe+Tyr
坝莜6号CK0.36bB0.44aA0.71bB0.99bA1.71aA2.03cB0.75bC40.99bC0.83cB
Se10.38aB0.44aB0.72bB1.01bB1.67bA2.12bB0.81aC41.26abC0.85bB
Se20.40aB0.46aB0.84aB1.06aB1.73aA2.28aB0.89aB42.33aB0.93aB
坝莜12号CK0.40aA0.47bA0.76aB1.01aA1.58aC2.12bA0.86bB44.21aB0.87bB
Se10.42aA0.50aA0.77aB1.01aB1.60aA2.19bA0.90bB44.47aB0.91aA
Se20.42aB0.51aA0.78aC1.03aB1.57aC2.32aA0.96aB43.28aB0.93aB
白燕2号CK0.42bA0.46bA0.83bA1.03bA1.67aB2.04cB0.92bA47.14bA0.91cA
Se10.45aA0.54aA0.93aA1.07bA1.64aA2.14bB0.98aA49.96aA0.98bA
Se20.46aA0.55aA0.98aA1.16aA1.61bB2.28aB1.03aA49.53aA1.02aA

新窗口打开| 下载CSV


3种燕麦材料的第一限制性氨基酸均为赖氨酸,其次是苏氨酸和亮氨酸,但叶面施硒可不同程度提高必需氨基酸的含量,且强化效果均为高浓度硒处理下最好。不同燕麦材料的SRC值排序为白燕2号>坝莜12号>坝莜6号。

对这些材料必需氨基酸指数进行计算,根据EAAI值实用评价标准,白燕2号为优质氨基酸源;其他均属于良好氨基酸源。且从整体来看,硒处理均能不同程度提高必需氨基酸含量、氨基酸比值、氨基酸比值系数、氨基酸比值系数分数和必需氨基酸指数,说明硒处理可以适当提高燕麦籽粒氨基酸的营养价值(表4)。

2.3 叶面施硒对不饱和脂肪酸组成及含量的影响

燕麦籽粒以不饱和脂肪酸为主要成分。由表5可知,3个燕麦品种主要以油酸和亚油酸2个不饱和脂肪酸成分为主,分别约占脂肪含量的40.00%和35.00%;其次是棕榈酸,约为15.00%;另外还包含约1.00%的亚麻酸;饱和脂肪酸(硬脂酸)含量极低,约2.00%。施硒对各类脂肪酸含量均有明显的促进作用。白燕2号各脂肪酸含量最高,高浓度处理Se2对各脂肪酸成分强化效果最好,与CK相比,Se2处理下白燕2号籽粒中油酸、亚油酸、棕榈酸和亚麻酸含量分别提高3.29%、5.11%、7.35%和46.36%。坝莜6号各脂肪酸含量最低。与对照相比,Se2处理下坝莜6号强化效果最好,籽粒中油酸、亚油酸、棕榈酸和亚麻酸含量分别提高0.36%、0.93%、1.66%和53.66%。不同燕麦品种籽粒脂肪酸成分对硒响应不同。不同浓度叶面硒对坝莜12号籽粒中棕榈酸强化效果最好(11.51%~ 19.40%),对坝莜6号籽粒中的亚麻酸强化效果最好(26.83%~53.66%)。品种、硒肥及其互作效应均达到极显著水平。

表5   叶面喷硒对燕麦成熟期籽粒脂肪酸含量的影响

Table 5  Effects of foliar selenium spraying the fatty acid contents in oat grains at maturity stage %

指标
Index
品种
Variety
处理TreatmentFF-Value
CKSe1Se2平均值Mean品种Ov硒肥Se品种×硒肥Ov×Se
脂肪Fat坝莜6号4.32cC4.55bC4.83aC4.57*********
坝莜12号6.81cB6.97bB7.19aB6.99
白燕2号7.62cA8.09bA8.17aA7.96
油酸Oleic acid坝莜6号39.23bC39.26bC39.37aC39.29*********
坝莜12号39.66cB40.00bA40.36aB40.01
白燕2号39.84bA39.86bB41.15aA40.28
亚油酸Linoleic acid坝莜6号34.25bC34.43aC34.57aC34.42*********
坝莜12号36.10cA36.23bB36.96aB36.43
白燕2号35.64cB37.31bA37.46aA36.80
棕榈酸Palmitic acid坝莜6号15.62bA15.69abA15.88aA15.73*********
坝莜12号12.94cC14.43bB15.45aB14.27
白燕2号14.82bB14.95bB15.91aA15.23
硬脂酸Stearic acid坝莜6号1.27cC1.68bC2.19aB1.71*********
坝莜12号1.87cA2.00bB2.16aB2.01
白燕2号1.80bB2.51aA2.55aA2.29
亚油酸Linoleic acid坝莜6号0.41cC0.52bC0.63aC0.52********
坝莜12号0.63bB0.75abB0.87aB0.75
白燕2号1.10cA1.26bA1.61aA1.32

新窗口打开| 下载CSV


2.4 叶面施硒对燕麦籽粒氨基酸和脂肪酸强化作用分析

叶面施硒可以显著强化成熟期燕麦籽粒硒、必需氨基酸、棕榈酸和硬脂酸的含量(图3a),且通过对其进行通径分析发现,硒浓度与成熟期硒含量、必需氨基酸、棕榈酸和硬脂酸存在显著的线性关系(P=0.000、0.004、0.031和0.004),其影响系数分别为0.891、0.689、0.584和0.695(图3b),进一步说明叶面施硒可显著提高上述物质的含量。

图3

图3   叶面施硒对燕麦籽粒氨基酸和脂肪酸的强化效果分析

(a)和(b)分别为籽粒氨基酸和脂肪酸的相关性分析和通径分析;(c)和(d)为籽粒中必需、非必需氨基酸和脂肪酸的相关性分析。SeIG:成熟期燕麦籽粒硒含量;TAA:总氨基酸;EAA:必需氨基酸;NEE:非必需氨基酸。

Fig.3   Analysis of strengthening effects of foliar selenium spraying on amino acids and fatty acids of oat grains

(a) and (b): the correlation and path analysis of amino acids and fatty acids in grains by foliar selenium, respectively; (c) and (d): correlation analysis of essential and non-essential amino acids and fatty acids in grains. SeIG: selenium content in oat grains at maturity stage; TAA: total amino acids; EAA: essential amino acid; NEE: non-essential amino acid.


此外,必需氨基酸、非必需氨基酸和脂肪酸成分(除棕榈酸外)间存在显著的相关性(图3a)。进一步分析发现,在必需氨基酸中,油酸和亚油酸与苯丙氨酸、赖氨酸和苏氨酸极显著相关,与亮氨酸和缬氨酸显著相关。除上述氨基酸外,硬脂酸还与异亮氨酸和蛋氨酸显著相关,亚麻酸还与异亮氨酸显著相关(图3c)。在非必需氨基酸中,油酸和亚油酸与谷氨酸、天冬氨酸、酪氨酸、甘氨酸、丙氨酸和丝氨酸显著相关,除上述氨基酸外,硬脂酸和亚麻酸还与组氨酸相关(图3d)。

进一步对其进行通径分析(表6)发现,必需氨基酸(P=0.017,拟合优度指数=0.935)和非必需氨基酸(P=0.042,拟合优度指数=0.969)模型拟合良好。必需氨基酸中(表7),缬氨酸对油酸和亚麻酸存在显著的线性关系(P=0.001和0.005),蛋氨酸对硬脂酸和亚麻酸存在显著线性关系(P= 0.041和0.017)。异亮氨酸与亚油酸线性相关(P= 0.019),其影响系数(标准化系数)为-0.404。苏氨酸与亚油酸也线性相关(P=0.003),影响系数为0.577。非必需氨基酸中,天冬氨酸与油酸、硬脂酸和亚麻酸线性相关(P=0.000、0.009和0.000),其影响系数(标准化系数)分别为1.515、-0.763和0.904。丝氨酸与油酸也存在显著线性相关(P=0.050)。

表6   通径分析模型拟合指标

Table 6  Fitting index of path analysis model

模型
Model
X2dfP
P-value
卡方
自由度比
Chi-square/df
拟合优
度指数
GFI
近似误差
均方根
RMSEA
均方根残差
RMSR
比较
拟合指数
CFI
规范
拟合系数
NFI
非规范
拟合系数
NNFI
必需氨基酸EAA15.44860.017**2.5750.9350.4440.0025480.9480.9350.526
非必需氨基酸NAA9.20360.042*1.5340.9690.2580.0024640.9860.9690.812

新窗口打开| 下载CSV


表7   通径分析模型回归系数

Table 7  Regression coefficient of the path analysis model

模型
Model
X→Y非标准化系数
USC
标准化系数
SC
标准误
SE
临界比值
CR
P
P-value
必需氨基酸EAA缬氨酸→油酸1.4931.2290.4583.2580.001***
异亮氨酸→亚油酸-0.592-0.4040.252-2.3470.019**
苏氨酸→亚油酸2.1830.5770.7322.9820.003***
蛋氨酸→硬脂酸0.2430.3790.1192.0430.041**
蛋氨酸→亚麻酸-0.218-0.3540.092-2.3810.017**
缬氨酸→亚麻酸0.4590.5890.1642.8010.005***
非必需氨基酸NAA天冬氨酸→油酸1.4651.5150.3833.8300.000***
丝氨酸→油酸1.2010.7450.6121.9620.050**
天冬氨酸→硬脂酸-0.491-0.7630.187-2.6310.009***
天冬氨酸→亚麻酸0.5610.9040.1593.5210.000***

新窗口打开| 下载CSV


此外,节点协方差数据(表8)表明,必需氨基酸中,苯丙氨酸与赖氨酸、苏氨酸,亮氨酸与赖氨酸、缬氨酸,赖氨酸与苏氨酸均有显著的关联性。非必需氨基酸中,丝氨酸与甘氨酸,丙氨酸与谷氨酸、甘氨酸、天冬氨酸和酪氨酸,谷氨酸与甘氨酸,甘氨酸与天冬氨酸和天冬氨酸,天冬氨酸与酪氨酸也均具有显著的关联。

表8   通径分析节点协方差关系

Table 8  Path analysis node covariance relationship

模型
Model
X↔Y非标准化估计系数
USEC
标准化估计系数
SEC
标准误
SE
z检验值
z-test value
P
P-value
必需氨基酸EAA苯丙氨酸↔赖氨酸0.1710.9490.0832.0650.039**
苯丙氨酸↔苏氨酸0.1590.9260.0782.0380.042**
亮氨酸↔赖氨酸0.3420.8950.1712.0010.045**
亮氨酸↔缬氨酸0.5180.9330.2532.0460.041**
赖氨酸↔苏氨酸0.0930.9200.0462.0310.042**
非必需氨基酸AA丝氨酸↔甘氨酸0.1060.8980.0532.0040.045**
丙氨酸↔谷氨酸0.5860.9030.2912.0110.044**
丙氨酸↔甘氨酸0.1110.9350.0542.0490.040**
丙氨酸↔天冬氨酸0.2000.9430.0972.0580.040**
丙氨酸↔酪氨酸0.1400.9060.0702.0150.044**
谷氨酸↔甘氨酸0.5670.9410.2762.0560.040**
谷氨酸↔天冬氨酸0.9960.9250.4892.0380.042**
甘氨酸↔天冬氨酸0.1780.8990.0892.0050.045**
天冬氨酸↔酪氨酸0.2390.9280.1172.0410.041**

新窗口打开| 下载CSV


3 讨论

3.1 叶面施硒对燕麦灌浆期籽粒硒积累的强化作用

叶片喷施外源硒是提高作物可食用部位硒含量的有效途径[13]。研究[24-25]发现,植物叶片以被动方式吸收亚硒酸盐,吸收量随喷施溶液中硒浓度的提高而持续增加,但硒强化需在国家富硒安全标准内,否则会对人体产生毒害作用。本研究中,Se2处理下强化了成熟期籽粒硒含量,且均在安全标准内。燕麦灌浆期籽粒硒积累整体呈现先增后减的趋势,可分为积累起始期、快速积累期、缓慢增加期、峰值出现期和缓慢回落期。但同一浓度下不同品种变化幅度不同,各时期出现的时间点也不同。这可能是因为品种硒积累能力不同有关[15]。且在花后7、21和28 d,燕麦品种和叶面硒肥互作达到显著差异水平。因此,在晋西北高寒生态区种植燕麦在灌浆初期和中后期进行叶面喷施硒肥可达到最佳效果。但品种间籽粒硒积累存在差异,成熟期籽粒硒含量下降的原因有待进一步了解。

3.2 叶面施硒对燕麦成熟期籽粒氨基酸的强化作用

燕麦是谷物中氨基酸最平衡的谷物之一,包含人体所必需的8种氨基酸,各种氨基酸含量均接近或高于FAO/WHO所规定的标准,明显优于大米、面粉等谷物蛋白质中必需氨基酸的构成比例[26]。Hu等[27]表明,适量施硒能促进水稻中氨基酸的合成。Liang等[28]发现,硒增加了谷子中苏氨酸合成酶和天冬氨酸合成酶的活性,并促进了相关氨基酸和总氨基酸的含量。本研究发现,Se2处理下3个燕麦品种总氨基酸及各种氨基酸含量均有提高,改善了燕麦的氨基酸营养。

导致氨基酸含量增加的重要因素可能是硒元素引起了植物体内硫和氮代谢的变化。利用无机硒可以提高蛋氨酸和胱氨酸的含量,以合成硒蛋氨酸和硒胱氨酸,原因可能是由于硒和硫同化互作,氮化合物合成过程中氮代谢途径可能发生了变化。可能是外源硒作为激活因子,激活了原本沉默的基因,从而直接产生新的酶或者增加原有通路中关键酶的含量和活性;或者通过干扰其他物质的吸收来间接影响原有代谢;或者通过干扰代谢途径中原本的酶、激活因子或辅助因子,从而使代谢途径流入其他方向,减少目标物质的合成;或者是通过调节促进或抑制目标物质生成的“协同物”或“拮抗物”的合成来调控目标物质的量[29-30]

本研究表明,燕麦籽粒中必需氨基酸对硒均有不同程度的响应。其中,品种和硒肥互作亮氨酸差异显著,硒肥单效应下异亮氨酸、赖氨酸、苏氨酸和苯丙氨酸差异显著,甲硫氨酸仅对硒肥单效应差异显著。且在强化分析中发现,苯丙氨酸与赖氨酸、苏氨酸,亮氨酸与赖氨酸、缬氨酸,赖氨酸与苏氨酸均有显著的关联性。所以,叶面硒强化必需氨基酸的原因可能是上述氨基酸间的“串联效应”。

本研究中,非必需氨基酸对硒肥单效应响应均达到极显著水平。其中,品种和硒肥互作谷氨酸、甘氨酸和丝氨酸差异显著,半胱氨酸仅对硒肥单效应差异显著。合成通路上,它们合成起始物质均来自糖代谢产生的关键中间体(α-酮戊二酸、草酰乙酸、3-磷酸甘油酸、丙酮酸、磷酸烯醇式丙酮酸和4-磷酸赤藓糖等),合成所需的氨基基团大部分来源于谷氨酸的转氨基反应。叶面喷硒提高燕麦各种氨基酸含量的原因之一可能是因为叶面硒提高了谷氨酸的含量。不同氨基酸对硒元素响应的具体机制有待进一步研究。

3.3 叶面施硒对燕麦成熟期籽粒脂肪酸的强化作用

燕麦虽为非油料作物,但裸燕麦脂肪平均含量可达6.30%,且不饱和脂肪酸含量占80.00%以上[31]。本研究结果表明,施硒对3个燕麦品种的脂肪酸含量均有明显的促进作用,品种、硒肥及其互作效应均达到极显著水平。Liu等[32]在谷子研究中发现,参与脂肪酸合成的还原蛋白Os04g0403500(K3ZVQ6)在富硒小米中上调,而还原剂的供应可以限制脂质的积累。Wang等[33]在稻谷中也发现施硒后能促进还原当量的增加,一些参与脂质降解的蛋白质也被下调,包括磷脂酶、D(K3ZQM7)和Os02g0654100(K3YUG6),进而促进脂肪酸的生物合成。

本研究中,脂肪酸(除棕榈酸外)与缬氨酸、异亮氨酸、苏氨酸、天冬氨酸、甲硫氨酸、天冬氨酸和丝氨酸存在不同程度的线性关系。推测脂肪酸的强化除了因叶面施硒引起的直接效应外,还有因氨基酸引起的间接效应。但叶面硒促进燕麦籽粒脂肪酸的机制尚不清楚,在碳链去饱和或延长的过程中,有哪些基因或酶受硒调控有待进一步研究验证。

4 结论

叶面施硒提高了燕麦籽粒中硒含量,优化了必需氨基酸构成比例,促进了不饱和脂肪酸含量,从而强化了燕麦营养品质。

参考文献

Dumont E, Vanhaecke F, Cornelis R.

Selenium speciation from food source to metabolites: a critical review

Analytical and Bioanalytical Chemistry, 2006, 385(7):1304-1323.

[本文引用: 1]

Jones G D, Droz B, Greve P, et al.

Selenium deficiency risk predicted to increase under future climate change

Proceedings of the National Academy of Sciences, 2017, 114(11):2848-2853.

[本文引用: 1]

Williame P N, Lombi E, Sun G, et al.

Selenium characterization in the global rice supply chain

Environmental Science & Technology, 2009, 43(15):6024-6030.

[本文引用: 1]

尹雪斌, 赵其国, 印遇龙, .

功能农业关键科学问题研究进展与“十四五”发展建议

科学通报, 2022, 67(6):497-510.

[本文引用: 1]

Schiavon M, Nardi S, Dalla Vecchia F, et al.

Selenium biofortification in the 21st century: status and challenges for healthy human nutrition

Plant and Soil, 2020,453:245-270.

[本文引用: 1]

Kaur S, Bhardwaj R D, Kapoor R, et al.

Biochemical characterization of oat (Avena sativa L.) genotypes with high nutritional potential

LWT-Food Science and Technology, 2019,110:32-39.

[本文引用: 1]

Mushtaq A, Gul-Zaffar Z. A D, et al.

A review on Oat (Avena sativa L.) as a dual-purpose crop

Scientific Research and Essays, 2014, 9(4):52-59.

[本文引用: 1]

Singh M, Kumar S. Broadening the Genetic Base of Grain Cereals. New Delhi: Springer (India) Private Limited, 2016.

[本文引用: 1]

Lyons G, Stangoulis J, Graham R.

High-selenium wheat: biofortification for better health

Nutrition Research Reviews, 2003, 16(1):45-60.

[本文引用: 1]

汤超华, 赵青余, 张凯, .

富硒农产品研究开发助力我国营养型农业发展

中国农业科学, 2019, 52(18):3122-3133.

[本文引用: 1]

Winkel L, Vriens B, Jones G, et al.

Selenium cycling across soil-plant-atmosphere interfaces: A critical review

Nutrients, 2015, 7(6):4199-4239.

[本文引用: 1]

Winkel L H E, Johnson C A, Lenz M, et al.

Environmental selenium research: from microscopic processes to global understanding

Environmental Science & Technology, 2012, 46 (2):571-579.

[本文引用: 1]

Ros G H, Rotterdam A M D, Bussink D W, et al.

Selenium fertilization strategies for bio-fortification of food: an agro- ecosystem approach

Plant and Soil, 2016,404:99-112.

[本文引用: 2]

李越, 宁丹, 吕玉峰, .

抽穗期叶面喷硒对燕麦品种硒积累影响及品种综合性状评价

作物杂志, 2023(4):1-11.

[本文引用: 4]

王敏. 燕麦种质资源主要微量元素含量的多样性及硒富集效应研究. 呼和浩特: 内蒙古农业大学, 2020.

[本文引用: 4]

张新军, 杨才, 曾昭海, .

叶面喷施硒肥对裸燕麦产量和品质的影响

麦类作物学报, 2015, 35(3):408-412.

[本文引用: 3]

孙发宇, 李长成, 王安, .

叶面喷施硒酸钠对不同小麦品种(系)籽粒硒及其他矿质元素含量的影响

麦类作物学报, 2017, 37(4):559-564.

[本文引用: 3]

Hao S N, Liu P F, Qin J, et al.

Effects of applying different doses of selenite to soil and foliar at different growth stage on selenium content and yield of different oat varieties

Plants, 2022, 11(14):1810.

[本文引用: 2]

Delaqua D, Carnier R, Berton R S, et al.

Increase of selenium concentration in wheat grains through foliar application of sodium selenate

Journal of Food Composition and Analysis, 2021,99:103886.

[本文引用: 1]

Yuan Z Q, Long W X, Liang T, et al.

Effect of foliar spraying of organic and inorganic selenium fertilizers during different growth stages on selenium accumulation and speciation in rice

Plant and Soil, 2023,486:87-101.

[本文引用: 1]

穆婷婷, 杜慧玲, 张福耀, .

外源硒对谷子生理特性、硒含量及其产量和品质的影响

中国农业科学, 2017, 50(1):51-63.

[本文引用: 1]

托列霍加·加吾提, 吾买尔夏提·塔汉, 隋晓青.

20份糜子材料的氨基酸含量分析及营养价值评价

种子, 2020, 39(7):31-36.

[本文引用: 1]

王超群. 基于脂肪酸组成对燕麦—小麦复配粉中燕麦粉的定量分析. 无锡: 江南大学, 2016.

[本文引用: 1]

张联合, 赵巍, 郁飞燕, .

水稻离体叶片吸收亚硒酸盐的生理特性

土壤学报, 2012, 49(1):189-193.

[本文引用: 1]

曹昌林, 吕慧卿, 郝志萍, .

外源硒对‘晋苦荞5号’光合特性、籽粒硒积累和产量及品质的影响

中国土壤与肥料, 2021(3):207-213.

[本文引用: 1]

Yue J X, Gu Z X, Zhu Z B, et al.

Impact of defatting treatment and oat varieties on structural, functional properties, and aromatic profile of oat protein

Food Hydrocolloids, 2021,112:106368.

[本文引用: 1]

Hu Z, Cheng Y, Suzuki N, et al.

Speciation of selenium in brown rice fertilized with selenite and effects of selenium fertilization on rice proteins

International Journal of Molecular Sciences, 2018, 19(11):3494.

[本文引用: 1]

Liang K H, Liang S, Zhu H.

Comparative proteomics analysis of the effect of selenium treatment on the quality of foxtail millet

LWT, 2020,131:109691.

[本文引用: 1]

Bielecka M, Watanabe M, Morcuende R, et al.

Transcriptome and metabolome analysis of plant sulfate starvation and resupply provides novel information on transcriptional regulation of metabolism associated with sulfur, nitrogen, and phosphorus nutritional responses in arabidopsis

Frontiers in Plant Science, 2015,5:805.

[本文引用: 1]

Harris J, Schneberg K A, Pilon-Smits E A.

Sulfur-selenium- molybdenum interactions distinguish selenium hyperaccumulator stanleya pinnata from non-hyperaccumulator Brassica Juncea (Brassicaceae)

Plant, 2014, 239(2):479-491.

[本文引用: 1]

相玉婷, 王晓龙, 胡新中, .

燕麦品种间脂肪酶活性差异及低脂肪酶优质品种的预测

中国农业科学, 2022, 55(21):4104-4117.

[本文引用: 1]

Liu M X, Zhang Z W, Ren G X, et al.

Evaluation of selenium and carotenoid concentrations of 200 foxtail millet accessions from china and their correlations with agronomic performance

Journal of Integrative Agriculture, 2016, 15(7):1449-1457.

[本文引用: 1]

Wang Y D, Wang X, Ngai S M, et al.

Comparative proteomics analysis of selenium responses in selenium-enriched rice grains

Journal of Proteome Research, 2013, 12(2):808-820.

[本文引用: 1]

/