Evaluation of Shade Tolerance of Peanut with Different Genotypes and Screening of Identification Indexes
HU TingHui, CHENG LiangQiang, WANG Jun,, Lü JianWei, RAO QingLinGuizhou Oil Research Institute, Guizhou Academy of Agricultural Sciences, Guiyang 550006通讯作者:
责任编辑: 杨鑫浩
收稿日期:2019-08-12接受日期:2019-10-8网络出版日期:2020-03-16
基金资助: |
Received:2019-08-12Accepted:2019-10-8Online:2020-03-16
作者简介 About authors
胡廷会,E-mail:hutinghui2010@163.com。
成良强,E-mail:872586048@qq.com。
摘要
关键词:
Abstract
Keywords:
PDF (552KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
胡廷会, 成良强, 王军, 吕建伟, 饶庆琳. 不同基因型花生耐荫性评价及其鉴定指标的筛选[J]. 中国农业科学, 2020, 53(6): 1140-1153 doi:10.3864/j.issn.0578-1752.2020.06.006
HU TingHui, CHENG LiangQiang, WANG Jun, Lü JianWei, RAO QingLin.
0 引言
【研究意义】花生(Arachis hypogaea L.)是我国主要出口农产品之一,在国民经济中具有重要地位,因其具有矮秆、固氮、投入少且产出高等特点,已逐渐成为贵州等西南山区与幼龄茶园、果园、药园、辣椒和玉米等作物进行间套作的理想先锋作物[1,2,3,4,5]。但是,由于间作套种条件下高秆作物遮荫加上耐荫品种缺乏,不可避免会对花生植株正常生长发育产生影响[6,7]。因此,弄清不同基因型花生的耐荫性,筛选出合理鉴定花生耐荫性的评价指标,不仅对花生耐荫资源筛选、耐荫品种(系)选育以及花生耐荫性提高具有重要意义,而且可有效提高土地资源利用率、优化农业产业结构和增加农民经济收益。【前人研究进展】近年来,关于花生耐荫性的研究已有相关报道。前人研究指出不同生育时期遮荫胁迫对花生生长、产量和品质的影响不尽相同。苗期遮荫对花生产量和品质影响较小[8];花针期遮光花生单株秕果数增加,饱果数减少,植株干物质生产显著降低,产量下降;在结荚期和饱果期遮荫花生含油量降低,蛋白质含量和可溶性总糖含量稍有增加[9];全生育期遮荫花生产量和品质均显著降低,叶片净光合速率降低,生物产量减少[10,11]。不同遮荫程度对花生的影响不同,中度遮荫胁迫下花生会通过增大叶绿体的受光面积、增加基粒数和基粒片层数来提高对光能的捕获和利用能力,但重度弱光胁迫会导致叶绿体基粒发育不完全,基粒片层破损,捕获光能的能力降低[12]。目前,关于花生耐荫性研究大部分集中于遮荫胁迫对花生植株生产、产量、品质以及光合生理等方面。针对植物的耐荫性鉴定科学家们提出了多种评价方法,孙艳等[13]以各单项指标的耐荫系数作为衡量耐荫性的指标,利用主成分分析、回归分析和聚类分析法对菊花耐荫性进行综合评价;梁颖等[14]根据阴天测定的叶绿素含量、净光合速率、叶绿素a/b和比叶重对甘蓝型油菜耐荫性进行鉴定;袁刘正等[15]以综合耐荫系数为指标,对不同玉米自交系耐荫性进行了评价;武晓玲等[16]采用多元统计法对大豆苗期耐荫性进行评价;李春红[17]、赵银月[18]等利用逐步回归分析法筛选出主茎节数、分枝数、节间长度、抗倒性、单株荚数、百粒重和单株粒重可作为大豆耐荫性的鉴定指标。【本研究切入点】目前关于作物耐荫性评价研究较多,而关于花生耐荫性评价及其鉴定指标的研究鲜有报道。不同作物对遮荫胁迫的响应不同,即使同一作物对荫蔽胁迫的响应也因品种(系)、生态环境不同而存在差异[19]。【拟解决的关键问题】本研究拟采用多元统计分析法,以花生形态、叶片光合特性、产量以及产量构成等指标的耐荫系数计算耐荫性综合评价值,对不同基因型花生耐荫性进行评价,建立花生耐荫性数学评价模型,明确花生耐荫性鉴定指标,为花生耐荫资源的鉴定以及耐荫品种的选育提供理论支撑。1 材料与方法
1.1 试验地概况及试验材料
试验于2017—2018年在贵州省农业科学院油料研究所贵阳试验基地(26°34' N,106°42'E)进行。其海拔高度为1 140 m,土壤类型为黄壤,地势平坦,无灌溉设施,试验地土壤肥力为中等,全氮含量3.17 g·kg-1,碱解氮108.26 mg·kg-1,全磷1.1 g·kg-1,速效磷47.50 mg·kg-1,全钾6.00 g·kg-1,速效钾351 mg·kg-1,有机质36.6 g·kg-1,交换性钙14.9 cmol·kg-1,pH 7.12。试验供试花生品种(系)为连续多年开展遮荫网试验,并对其主茎高、侧枝长、分枝数、结果枝数、单株产量等进行比较分析,初步筛选出耐荫性不同的30个花生品种(系)(表1)。玉米品种为黔糯868。
Table 1
表1
表1参试材料信息
Table 1
材料 Material | 品种类型 Cultivar type | 选育单位 Breeding institute | 生育期 Growth period (d) |
---|---|---|---|
201240401A | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 132 |
201060401C | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 134 |
201070406 | 普通型Virginia type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 130 |
2010605 | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 130 |
201240413 | 普通型Virginia type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 133 |
201070407 | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 132 |
徐花9号Xuhua 9 | 中间型Irregular type | 江苏省徐州市农业科学院Institute of Agricultural Sciences of Xuzhou | 129 |
201150119 | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 130 |
花育21号Huayu 21 | 普通型Virginia type | 山东省花生研究所Shandong Peanut Research Institute | 133 |
201241001 | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 132 |
天府29号Tianfu 29 | 珍珠豆型Spanish type | 四川省南充市农业科学院Nanchong Institute of Agricultural Sciences | 130 |
6-2 | 普通型Virginia type | 山东省花生研究所Shandong Peanut Research Institute | 128 |
201060105 | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 133 |
201150104 | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 130 |
冀花9号Jihua 9 | 普通型Virginia type | 河北省农林科学院粮油作物研究所 Institute of Cereal and Oil Crops, Hebei Academy of Agricultural and Forestry Sciences | 129 |
花育33号Huayu 33 | 普通型Virginia type | 山东省花生研究所Shangdong Peanut Research Institute | 128 |
201060106B1 | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 132 |
201150118A | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 130 |
徐花13号Xuhua 13 | 中间型Irregular type | 江苏省徐州市农业科学院Institute of Agricultural Sciences of Xuzhou | 128 |
天府22号Tianfu 22 | 中间型Irregular type | 四川省南充市农业科学院Nanchong Institute of Agricultural Sciences | 130 |
201150201 | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 135 |
远杂6号Yuanza 6 | 珍珠豆型Spanish type | 河南省农业科学院Henan Academy of Agricultural Sciences | 131 |
闽花6号Minhua 6 | 珍珠豆型Spanish type | 福建农林大学作物科学学院 College of Crop Science, Fujian Agriculture and Forestry University | 129 |
天府18号Tianfu 18 | 普通型Virginia type | 四川省南充市农业科学院Nanchong Institute of Agricultural Sciences | 130 |
201150114A | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 131 |
中花6号Zhonghua 6 | 珍珠豆型Spanish type | 中国农业科学院油料作物研究所 Oil Crops Research Institute, Chinese Academy of Agricultural Sciences | 123 |
201240404 | 珍珠豆型Spanish type | 贵州省农业科学院油料研究所 Oil Research Institute, Guizhou Academy of Agricultural Sciences | 132 |
花育17号Huayu 17 | 普通型Virginia type | 山东省花生研究所Shandong Peanut Research Institute | 130 |
花育19号Huayu 19 | 普通型Virginia type | 山东省花生研究所Shandong Peanut Research Institute | 130 |
冀花10号Jihua 10 | 普通型Virginia type | 河北省农林科学院粮油作物研究所 Institute of Cereal and Oil Crops,Hebei Academy of Agricultural and Forestry Sciences | 129 |
新窗口打开|下载CSV
1.2 试验设计
试验为裂区设计,将光照设为主区,30个花生品种(系)设为副区。其中光照设遮荫(ZY)和正常光照(CK)2个水平。遮荫按照玉米与花生各种植1行(1﹕1)的比例进行间作(其遮光率在玉米开花期后达到72%—75%),正常光照为净作花生。玉米间作花生和净作花生均为平作,小区行长3 m,4行区,行距为0.6 m,玉米株距为0.15 m,每穴播2粒,留1苗,花生穴距为0.2 m,每穴播种3粒,留2苗。品种(系)间不留过道,随机排列,每个处理重复3次。2017年和2018年分别于4月24日和4月27日播种,9月26日和9月28日收获,收获时按田间实收进行计产。每公顷施用750 kg三元复合肥作基肥,随播种前整地施入土壤中,其他田间管理和病虫草害防治方法与普通大田花生生产一致。1.3 测定项目和方法
1.3.1 形态指标及产量性状 在花生成熟期每个处理选取具有代表性的10株,测定主茎高(X1)、侧枝长(X2)、总分枝数(X3)、有效分枝数(X4),单株结果数(X5)、单株饱果数(X6)和单株产量(X11);收获晾晒干后分别测定百果重(X7)、百果仁重(X8)、百仁重(X9)、出仁率(X10)和小区产量(X12)。1.3.2 光合参数测定 在天气晴朗的上午9:30—11:30,采用美国产的Li-6400便携式光合仪,在花生结荚期(2017年8月15—16日,2018年8月19—20日)选择具有代表性的植株,测定冠层上部倒数第3片叶的净光合速率(X13)、气孔导度(X14)、胞间CO2浓度(X15)和蒸腾速率(X16)。
1.4 数据处理与分析
采用Excel 2007软件进行数据输入、整理,采用SPSS 20.0软件进行方差分析、主成分分析、聚类分析和逐步回归等分析。计算各单项指标的耐荫系数(shade tolerance coefficient):RVS=遮荫处理下性状值/净作处理性状值、各综合指标的隶属函数值u(Xj)=(Xj-Xmin)/(Xmax-Xmin)、权重$w_{j}=p_{j}/\sum^{n}_{j=1}p_{j}$
式中,j=1,2,3,…,n,Xmin表示第j个综合指标的最小值,Xmax表示第j个综合指标的最大值,Xj表示第j个综合指标。wj表示第j个综合指标在所有综合指标中所占的权重;pj为各花生品种(系)第j个综合指标的贡献率。D值为各基因型花生在遮荫条件下耐荫性综合评价值。
2 结果
2.1 不同基因型花生各单项指标的耐荫系数及其相关性分析
采用2年的平均值,对30个花生品种(系)16个单项指标的耐荫系数进行计算分析(表2)可知,在各指标中,蒸腾速率和单株饱果数的变异系数大,分别为40.84%和36.56%,说明遮荫处理对蒸腾速率和单株饱果数的影响大;其次为单株结果数、百果重、百果仁重、单株产量、小区产量、净光合速率、气孔导度和胞间CO2浓度,其变异系数为22.38%—33.71%,说明遮荫胁迫对花生单株结果数、百果重、百果仁重、单株产量、小区产量、净光合速率、气孔导度和胞间CO2浓度的影响较大。不同花生品种(系)各指标的耐荫系数表现各异,遮荫条件下其主茎高、气孔导度、胞间CO2浓度和蒸腾速率总体上较正常光照有所上升(RVS>1),有效分枝数、单株结果数、单株饱果数、百果重、百果仁重、百仁重、出仁率、单株产量和小区产量与对照相比有所下降(RVS<1),而侧枝长、总分枝数和净光合速率则表现出种质差异,有的上升(RVS>1),有的下降(RVS<1)。可见,不同花生材料各单项指标的变异幅度不同,因此,难以采用单一指标的耐荫系数来判断花生种质的耐荫性。Table 2
表2
表2花生各单项指标的耐荫系数
Table 2
品种(系) Cultivar (Line) | X1 | X2 | X3 | X4 | X5 | X6 | X7 | X8 | X9 | X10 | X11 | X12 | X13 | X14 | X15 | X16 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
6-2 | 1.04 | 0.89 | 0.60 | 0.51 | 0.48 | 0.47 | 0.77 | 0.64 | 0.94 | 0.83 | 0.41 | 0.40 | 1.06 | 1.85 | 1.73 | 0.92 |
花育17号Huayu 17 | 1.18 | 1.10 | 0.66 | 0.52 | 0.48 | 0.46 | 0.62 | 0.52 | 0.68 | 0.82 | 0.26 | 0.32 | 0.74 | 1.29 | 1.78 | 0.60 |
花育19号Huayu 19 | 0.98 | 0.93 | 0.72 | 0.63 | 0.50 | 0.48 | 0.75 | 0.59 | 0.81 | 0.79 | 0.37 | 0.28 | 0.82 | 1.30 | 1.66 | 0.75 |
花育21号Huayu 21 | 1.24 | 1.12 | 0.76 | 0.63 | 0.52 | 0.41 | 0.69 | 0.56 | 0.88 | 0.81 | 0.42 | 0.48 | 0.72 | 1.45 | 1.62 | 1.15 |
花育33号Huayu 33 | 1.16 | 1.01 | 0.76 | 0.63 | 0.47 | 0.51 | 0.58 | 0.55 | 0.78 | 0.93 | 0.43 | 0.36 | 0.72 | 1.26 | 1.98 | 0.90 |
冀花10号Jihua 10 | 1.03 | 0.89 | 0.71 | 0.56 | 0.43 | 0.39 | 0.65 | 0.58 | 0.82 | 0.89 | 0.32 | 0.28 | 0.80 | 1.51 | 2.14 | 0.83 |
冀花9号Jihua 9 | 1.08 | 0.94 | 0.72 | 0.57 | 0.51 | 0.55 | 0.59 | 0.51 | 0.72 | 0.82 | 0.28 | 0.32 | 0.68 | 1.25 | 1.61 | 0.52 |
闽花6号Minhua 6 | 1.15 | 0.99 | 0.61 | 0.51 | 0.31 | 0.30 | 0.52 | 0.49 | 0.68 | 0.94 | 0.28 | 0.32 | 0.79 | 1.48 | 3.56 | 0.97 |
天府18号Tianfu 18 | 0.97 | 0.77 | 0.77 | 0.53 | 0.39 | 0.32 | 0.60 | 0.51 | 0.69 | 0.85 | 0.32 | 0.30 | 0.64 | 1.29 | 3.52 | 1.74 |
天府22号Tianfu 22 | 1.06 | 0.70 | 0.71 | 0.57 | 0.34 | 0.39 | 0.78 | 0.67 | 0.94 | 0.85 | 0.40 | 0.36 | 0.95 | 1.65 | 2.73 | 1.22 |
天府29号Tianfu 29 | 0.97 | 0.79 | 0.72 | 0.49 | 0.39 | 0.37 | 0.81 | 0.69 | 0.96 | 0.86 | 0.44 | 0.36 | 0.97 | 2.02 | 2.80 | 1.26 |
徐花13号Xuhua 13 | 1 | 0.83 | 0.71 | 0.53 | 0.35 | 0.39 | 0.64 | 0.56 | 0.76 | 0.85 | 0.37 | 0.33 | 0.77 | 1.44 | 3.34 | 1.44 |
徐花9号Xuhua 9 | 1.06 | 0.80 | 0.73 | 0.58 | 0.40 | 0.38 | 0.70 | 0.61 | 0.87 | 0.87 | 0.40 | 0.46 | 0.83 | 1.51 | 3.16 | 1.13 |
远杂6号Yuanza 6 | 0.92 | 0.87 | 0.67 | 0.59 | 0.43 | 0.35 | 0.65 | 0.53 | 0.63 | 0.78 | 0.29 | 0.34 | 0.76 | 1.07 | 2.90 | 0.95 |
中花6号Zhonghua 6 | 0.93 | 0.78 | 0.59 | 0.45 | 0.30 | 0.28 | 0.57 | 0.48 | 0.60 | 0.81 | 0.28 | 0.31 | 0.73 | 1.33 | 3.40 | 1.27 |
201060105 | 1.05 | 0.83 | 0.74 | 0.63 | 0.58 | 0.58 | 0.62 | 0.54 | 0.74 | 0.87 | 0.36 | 0.35 | 0.69 | 1.23 | 1.57 | 1.02 |
201060106B1 | 0.89 | 0.75 | 0.79 | 0.66 | 0.52 | 0.52 | 0.75 | 0.68 | 0.90 | 0.91 | 0.37 | 0.46 | 0.86 | 1.40 | 2.03 | 0.91 |
201060401C | 0.97 | 0.83 | 0.72 | 0.60 | 0.39 | 0.25 | 0.61 | 0.51 | 0.67 | 0.79 | 0.31 | 0.37 | 0.69 | 1.10 | 2.22 | 1.68 |
2010605 | 1.15 | 0.92 | 0.78 | 0.68 | 0.44 | 0.37 | 0.52 | 0.50 | 0.62 | 0.89 | 0.36 | 0.31 | 0.58 | 0.98 | 2.48 | 1.95 |
201070406 | 1.05 | 0.90 | 0.76 | 0.55 | 0.37 | 0.35 | 0.50 | 0.42 | 0.54 | 0.78 | 0.21 | 0.29 | 0.48 | 0.91 | 2.68 | 1.01 |
201070407 | 0.93 | 0.77 | 0.64 | 0.54 | 0.29 | 0.24 | 0.53 | 0.46 | 0.57 | 0.82 | 0.27 | 0.43 | 0.63 | 1.01 | 3.52 | 1.60 |
201150104 | 1.22 | 1.03 | 1.05 | 0.82 | 0.61 | 0.57 | 0.65 | 0.54 | 0.76 | 0.82 | 0.29 | 0.33 | 0.50 | 0.98 | 1.85 | 0.73 |
201150114A | 0.9 | 0.73 | 0.73 | 0.55 | 0.33 | 0.30 | 0.59 | 0.51 | 0.66 | 0.83 | 0.29 | 0.27 | 0.72 | 1.20 | 2.74 | 1.17 |
201150118A | 0.99 | 0.79 | 0.69 | 0.54 | 0.35 | 0.37 | 0.54 | 0.46 | 0.60 | 0.81 | 0.48 | 0.59 | 0.65 | 1.07 | 2.58 | 2.10 |
201150119 | 1.14 | 0.93 | 0.69 | 0.69 | 0.56 | 0.51 | 0.57 | 0.49 | 0.64 | 0.83 | 0.48 | 0.36 | 0.64 | 0.96 | 1.69 | 1.30 |
201150201 | 1.18 | 0.95 | 0.64 | 0.57 | 0.44 | 0.55 | 0.71 | 0.64 | 0.86 | 0.92 | 0.36 | 0.45 | 1.03 | 1.61 | 2.28 | 0.86 |
201240404 | 1.11 | 0.88 | 0.84 | 0.63 | 0.44 | 0.58 | 0.56 | 0.53 | 0.72 | 0.91 | 0.34 | 0.33 | 0.64 | 1.14 | 2.21 | 0.68 |
201241001 | 1.04 | 0.84 | 0.80 | 0.59 | 0.38 | 0.45 | 0.77 | 0.66 | 0.95 | 0.85 | 0.39 | 0.45 | 0.83 | 1.62 | 2.27 | 0.86 |
201240401A | 1.07 | 0.96 | 0.85 | 0.76 | 0.56 | 0.64 | 0.70 | 0.61 | 0.81 | 0.87 | 0.40 | 0.59 | 0.70 | 1.17 | 2.14 | 1.18 |
201240413 | 1.42 | 1.03 | 0.74 | 0.61 | 0.44 | 0.49 | 0.71 | 0.56 | 0.81 | 0.79 | 0.47 | 0.40 | 0.76 | 1.50 | 2.15 | 1.00 |
平均值Average | 1.06 | 0.88 | 0.73 | 0.59 | 0.43 | 0.43 | 0.64 | 0.55 | 0.75 | 0.85 | 0.35 | 0.37 | 0.75 | 1.32 | 2.41 | 1.12 |
标准差Stdev | 0.18 | 0.15 | 0.12 | 0.12 | 0.12 | 0.15 | 0.10 | 0.09 | 0.14 | 0.06 | 0.11 | 0.12 | 0.19 | 0.38 | 0.77 | 0.50 |
变异系数CV(%) | 17.00 | 17.72 | 16.03 | 20.39 | 30.47 | 36.56 | 22.38 | 23.10 | 19.92 | 7.05 | 29.76 | 33.71 | 29.06 | 28.99 | 29.38 | 40.84 |
新窗口打开|下载CSV
对16个单项指标进行相关性分析,获得相关系数矩阵(表3),除出仁率外,其他各个单项指标间均存在一定的相关性,其中,单株饱果数与9个指标存在显著或极显著相关性,单株结果数与8个指标存在显著或极显著相关性。由于各单项指标耐荫系数间或多或少存在一定相关性,不可避免使得各个单项指标所提供的信息会发生重叠,表明花生耐荫性是一个复杂的综合性状。另外,各单项指标在花生耐荫中所起的作用也可能不同,直接利用各单项指标不能直观、准确地评价花生耐荫性。因此,为弥补单项指标耐荫性评价所存在的不足,需在此基础上采用多元统计方法进行综合评价分析。
Table 3
表3
表3各单项指标耐荫系数的相关系数矩阵
Table 3
指标Index | X1 | X2 | X3 | X4 | X5 | X6 | X7 | X8 | X9 | X10 | X11 | X12 | X13 | X14 | X15 | X16 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
X1 | 1 | |||||||||||||||
X2 | 0.799** | 1 | ||||||||||||||
X3 | 0.245 | 0.178 | 1 | |||||||||||||
X4 | 0.357 | 0.356 | 0.789** | 1 | ||||||||||||
X5 | 0.398* | 0.525** | 0.520** | 0.753** | 1 | |||||||||||
X6 | 0.455* | 0.429* | 0.487** | 0.621** | 0.827** | 1 | ||||||||||
X7 | 0.008 | -0.089 | 0.079 | 0.053 | 0.236 | 0.303 | 1 | |||||||||
X8 | 0.000 | -0.129 | 0.118 | 0.095 | 0.235 | 0.380* | 0.935** | 1 | ||||||||
X9 | 0.168 | 0.024 | 0.151 | 0.088 | 0.275 | 0.403* | 0.914** | 0.951** | 1 | |||||||
X10 | 0.089 | -0.009 | 0.062 | 0.083 | 0.035 | 0.262 | -0.002 | 0.330 | 0.289 | 1 | ||||||
X11 | 0.294 | 0.035 | 0.030 | 0.210 | 0.253 | 0.340 | 0.449* | 0.455* | 0.511** | 0.111 | 1 | |||||
X12 | 0.085 | 0.010 | 0.067 | 0.224 | 0.129 | 0.252 | 0.261 | 0.290 | 0.288 | 0.062 | 0.547** | 1 | ||||
X13 | -0.106 | -0.178 | -0.478** | -0.397* | -0.107 | 0.062 | 0.762** | 0.791** | 0.745** | 0.263 | 0.355 | 0.175 | 1 | |||
X14 | 0.037 | -0.102 | -0.313 | -0.457* | -0.154 | 0.016 | 0.753** | 0.770** | 0.815** | 0.263 | 0.356 | 0.099 | 0.878** | 1 | ||
X15 | -0.406* | -0.522** | -0.312 | -0.522** | -0.827** | -0.719** | -0.269 | -0.230 | -0.316 | 0.098 | -0.294 | -0.099 | -0.030 | 0.041 | 1 | |
X16 | -0.282 | -0.440* | -0.147 | -0.124 | -0.429* | -0.554** | -0.309 | -0.319 | -0.371* | -0.121 | 0.230 | 0.259 | -0.243 | -0.239 | 0.476** | 1 |
新窗口打开|下载CSV
2.2 主成分分析
以30个花生品种(系)的16个单项指标耐荫系数为基础进行主成分分析,根据特征值>1的原则,共提取得到5个主成分,并将16个单项指标转换为5个综合指标(CI1—CI5),各主成分的特征向量及贡献率如表4所示。由主成分分析特征值可知,前5个主成分的累计贡献率高达87.127%,具有较强的信息代表性,代表了原始指标所携带的绝大部分信息。其中,第1主成分可反映原始数据信息量的33.860%,代表了5.418个原始指标,包括单株产量、单株饱果数、百果重、百果仁重、百仁重等指标,主要反映花生产量及产量构成因素等特征;第2主成分可反映原始数据信息量的26.666%,代表了4.267个原始指标,包括侧枝长、总分枝数、有效分枝数,单株结果数、净光合速率和气孔导度等指标,主要反映花生的形态和光合特征;第3主成分可反映原始数据信息量的11.176%,代表了1.788个原始指标,包括小区产量和蒸腾速率;第4主成分可反映原始数据信息量的8.471%,代表了1.355个原始指标,起主要贡献的是株高,反映花生形态特性;第5主成分可反映原始数据信息量的6.954%,代表了1.113个原始指标,起主要贡献的是出仁率,主要反映花生荚果饱满程度。Table 4
表4
表4各性状主成分的特征向量及贡献率
Table 4
主成分Principle factor | CI1 | CI2 | CI3 | CI4 | CI5 | |
---|---|---|---|---|---|---|
特征向量Eigenvector | X1 | 0.413 | -0.477 | -0.201 | 0.587* | 0.203 |
X2 | 0.331 | -0.598* | -0.400 | 0.505 | 0.054 | |
X3 | 0.310 | -0.614* | 0.283 | -0.473 | 0.113 | |
X4 | 0.408 | -0.746* | 0.353 | -0.238 | 0.063 | |
X5 | 0.645 | -0.646* | 0.000 | -0.105 | -0.140 | |
X6 | 0.751* | -0.489 | -0.015 | -0.081 | 0.100 | |
X7 | 0.792* | 0.462 | 0.048 | -0.164 | -0.281 | |
X8 | 0.831* | 0.474 | 0.077 | -0.232 | 0.031 | |
X9 | 0.879* | 0.406 | 0.013 | -0.092 | 0.016 | |
X10 | 0.274 | 0.132 | -0.015 | -0.110 | 0.911* | |
X11 | 0.569* | 0.136 | 0.550 | 0.426 | -0.030 | |
X12 | 0.347 | 0.062 | 0.680* | 0.323 | 0.005 | |
X13 | 0.554 | 0.768* | -0.164 | 0.088 | -0.021 | |
X14 | 0.556 | 0.751* | -0.208 | 0.120 | 0.056 | |
X15 | -.638* | 0.512 | 0.132 | -0.059 | 0.327 | |
X16 | -0.514 | 0.181 | 0.722* | 0.216 | 0.017 | |
特征值Eigenvalue | 5.418 | 4.267 | 1.788 | 1.355 | 1.113 | |
贡献率Contribution rate | 33.860 | 26.666 | 11.176 | 8.471 | 6.954 | |
累计贡献率Cumulative contribution rate | 33.860 | 60.526 | 71.702 | 80.173 | 87.127 |
新窗口打开|下载CSV
2.3 耐荫性综合评价
2.3.1 隶属函数分析 对各花生品种(系)综合指标进行隶属函数分析,获得隶属函数值(表5),对于综合指标CI1,主要反映花生产量及其产量构成因素特征,其隶属函数值最大的是6-2,为1.00,最小的是201070407,为0.00,可见在该综合指标下,6-2的耐荫性好,而201070407的耐荫性差;对于综合指标CI2,主要反映花生的形态、光合特性,其中以天府29的隶属函数值最大,201150104的隶属函数值最小,可见在天府29在 CI2这一综合指标表现强耐荫性,而耐荫性最差的是201150104。Table 5
表5
表5各花生品种(系)耐荫性综合评价
Table 5
品种(系) Cultivar (Line) | CI1 | CI2 | CI3 | CI4 | CI5 | U(X1) | U(X2) | U(X3) | U(X4) | U(X5) | D值 D value | 综合评价 Comprehensive valuation |
---|---|---|---|---|---|---|---|---|---|---|---|---|
6-2 | 1.266 | 1.358 | -0.700 | 0.767 | -1.211 | 1.000 | 0.851 | 0.298 | 0.640 | 0.140 | 0.761 | 耐荫型 Shade tolerant |
花育17号Huayu 17 | 0.021 | -0.676 | -2.157 | 0.988 | -0.564 | 0.615 | 0.412 | 0.000 | 0.693 | 0.284 | 0.455 | 中度耐荫 Medium shade tolerant |
花育19号Huayu 19 | 0.651 | -0.123 | -0.777 | -0.750 | -1.838 | 0.810 | 0.531 | 0.282 | 0.279 | 0.000 | 0.541 | 中度耐荫 Medium shade tolerant |
花育21号Huayu 21 | 0.991 | -0.646 | 0.222 | 1.771 | -0.776 | 0.915 | 0.419 | 0.486 | 0.880 | 0.237 | 0.651 | 耐荫型Shade tolerant |
花育33号Huayu 33 | 0.571 | -0.687 | -0.189 | 0.676 | 1.619 | 0.785 | 0.410 | 0.402 | 0.619 | 0.771 | 0.604 | 耐荫型Shade tolerant |
冀花10号Jihua 10 | 0.193 | 0.405 | -1.125 | -0.537 | 0.536 | 0.668 | 0.645 | 0.211 | 0.330 | 0.530 | 0.559 | 中度耐荫 Medium shade tolerant |
冀花9号Jihua 9 | 0.094 | -0.865 | -1.423 | -0.260 | -0.777 | 0.638 | 0.371 | 0.150 | 0.396 | 0.237 | 0.438 | 中度耐荫 Medium shade tolerant |
闽花6号Minhua 6 | -0.963 | 0.588 | -1.535 | 1.162 | 2.644 | 0.311 | 0.685 | 0.127 | 0.735 | 1.000 | 0.498 | 中度耐荫 Medium shade tolerant |
天府18号Tianfu 18 | -1.158 | 0.423 | 0.580 | -0.753 | 0.568 | 0.251 | 0.649 | 0.560 | 0.278 | 0.537 | 0.438 | 中度耐荫 Medium shade tolerant |
天府22号Tianfu 22 | 0.683 | 1.596 | 0.416 | -0.663 | -0.096 | 0.820 | 0.902 | 0.526 | 0.300 | 0.389 | 0.722 | 耐荫型Shade tolerant |
天府29号Tianfu 29 | 0.950 | 2.051 | 0.186 | -0.322 | -0.251 | 0.902 | 1.000 | 0.479 | 0.381 | 0.354 | 0.783 | 耐荫型Shade tolerant |
徐花13号Xuhua 13 | -0.489 | 0.822 | 0.253 | -0.069 | 0.442 | 0.457 | 0.735 | 0.493 | 0.441 | 0.509 | 0.549 | 中度耐荫 Medium shade tolerant |
徐花9号Xuhua 9 | 0.354 | 0.933 | 0.764 | -0.034 | 0.616 | 0.718 | 0.759 | 0.597 | 0.450 | 0.547 | 0.675 | 耐荫型Shade tolerant |
远杂6号Yuanza 6 | -0.756 | 0.074 | -0.397 | -0.648 | -1.339 | 0.375 | 0.574 | 0.360 | 0.303 | 0.111 | 0.406 | 中度耐荫 Medium shade tolerant |
中花6号Zhonghua 6 | -1.647 | 0.981 | -0.639 | 0.116 | -0.448 | 0.099 | 0.769 | 0.310 | 0.485 | 0.310 | 0.386 | 中度耐荫 Medium shade tolerant |
201060105 | 0.420 | -0.803 | -0.031 | -0.653 | -0.017 | 0.738 | 0.385 | 0.435 | 0.302 | 0.406 | 0.522 | 中度耐荫 Medium shade tolerant |
201060106B1 | 1.226 | 0.441 | 0.940 | -1.920 | 0.521 | 0.988 | 0.653 | 0.633 | 0.000 | 0.526 | 0.707 | 耐荫型Shade tolerant |
201060401C | -0.996 | 0.003 | 0.677 | -0.105 | -1.237 | 0.301 | 0.559 | 0.579 | 0.433 | 0.134 | 0.415 | 中度耐荫 Medium shade tolerant |
2010605 | -0.849 | -1.023 | 0.939 | 0.304 | 1.367 | 0.346 | 0.337 | 0.633 | 0.530 | 0.715 | 0.428 | 中度耐荫 Medium shade tolerant |
201070406 | -1.735 | -1.017 | -0.922 | -0.365 | -0.645 | 0.072 | 0.339 | 0.253 | 0.371 | 0.266 | 0.221 | 敏感型 Shade susceptible |
201070407 | -1.968 | 0.405 | 0.717 | 0.045 | -0.065 | 0.000 | 0.645 | 0.588 | 0.468 | 0.396 | 0.350 | 中度耐荫 Medium shade tolerant |
201150104 | 0.784 | -2.589 | 0.071 | -1.614 | -0.101 | 0.851 | 0.000 | 0.456 | 0.073 | 0.388 | 0.427 | 中度耐荫 Medium shade tolerant |
201150114A | -1.171 | 0.441 | -0.256 | -1.276 | -0.282 | 0.246 | 0.653 | 0.389 | 0.153 | 0.347 | 0.388 | 中度耐荫 Medium shade tolerant |
201150118A | -1.068 | 0.106 | 2.734 | 1.762 | -0.591 | 0.278 | 0.581 | 1.000 | 0.878 | 0.278 | 0.522 | 中度耐荫 Medium shade tolerant |
201150119 | 0.061 | -1.305 | 0.755 | 1.060 | -0.596 | 0.627 | 0.277 | 0.596 | 0.710 | 0.277 | 0.496 | 中度耐荫 Medium shade tolerant |
201150201 | 1.200 | 0.749 | -0.692 | 0.867 | 1.246 | 0.980 | 0.720 | 0.300 | 0.664 | 0.688 | 0.759 | 耐荫型Shade tolerant |
201240404 | 0.132 | -0.944 | -0.398 | -0.894 | 1.628 | 0.649 | 0.354 | 0.360 | 0.245 | 0.773 | 0.492 | 中度耐荫 Medium shade tolerant |
201241001 | 1.036 | 0.793 | 0.333 | -0.630 | -0.262 | 0.929 | 0.729 | 0.509 | 0.307 | 0.352 | 0.707 | 耐荫型Shade tolerant |
201240401A | 1.244 | -1.063 | 1.814 | -0.301 | 0.506 | 0.993 | 0.329 | 0.812 | 0.386 | 0.523 | 0.670 | 耐荫型Shade tolerant |
201240413 | 0.923 | -0.429 | -0.161 | 2.275 | -0.597 | 0.894 | 0.465 | 0.408 | 1.000 | 0.277 | 0.662 | 耐荫型Shade tolerant |
权重Index weight | 0.389 | 0.306 | 0.128 | 0.097 | 0.080 |
新窗口打开|下载CSV
2.3.2 权重确定 根据各个综合指标贡献率大小,对其权重进行计算(表5)。经计算获得5个综合指标的权重分别为0.389、0.306、0.128、0.097和0.080。2.3.3 综合评价及分类 对各花生品种(系)的耐荫综合评价值D值进行计算,并根据D 值大小对各花生材料耐荫性进行排序(表5)。其中,D值最大的是天府29号,为0.783,说明其耐荫性最强,D值最小的是201070406,为0.221,说明其耐荫性最差。根据D值采用欧式距离法对其进行聚类分析(图1),可将30个不同基因型花生材料划分为3 类,其中,6-2、201150201、天府29号、天府22号、201060106B1、201241001、201240401A、201240413、花育21号、花育33号和徐花9号为第一类,占供试材料的36.67%,其耐荫综合评价值大,属于耐荫型品种(系);花育17号、花育19号、冀花10号、冀花9号、闽花6号、天府18号、徐花13号、远杂6号、中花6号、201060105、201060401C、2010605、201150119、201150104、201150114A、201150118A、201070407、201240404属于第二类,占供试材料的60.00%,其耐荫综合评价值较大,属于中度耐荫型品种(系);201070406属于第三类,占供试材料的3.33%,其耐荫综合评价最小,属于敏感型品种(系)。综上所述,参试的30个花生品种(系)中,有11个属于耐荫品种(系),18个属于中度耐荫品种(系),而敏感型品种(系)较少,只有1个。
图1
新窗口打开|下载原图ZIP|生成PPT图130个花生品种(系)聚类树状图
Fig. 1The dendrogram of clusters for 30 peanut cultivars (line)
2.4 回归分析及耐荫鉴定指标的选择
为更好的对花生耐荫性进行预测,弄清各个单项指标与耐荫性之间的关系,筛选出鉴定花生耐荫性的重要指标,获得耐荫性评价数学模型,以耐荫性综合评价值(D值)作因变量,各单项指标的耐荫系数为自变量进行逐步回归分析,建立最优回归方程:D=-0.741+0.576X9+0.507X11+0.298X13+0.272X12+0.406X10(R2=0.990,F=483.497,P=0.000),并对其估计精度进行分析可知(表6),各个花生品种(系)的估计精度均在93.18%以上,表明方程中5个指标所代表的百仁重、出仁率、单株产量、小区产量和净光合速率对遮荫条件响应较为明显,可作为花生耐荫性评价指标。因此,可在相同条件下测定其他花生品种(系)的上述5个指标,并求得其耐荫系数,利用该方程计算耐荫性综合评价D值,便可预测其耐荫性。Table 6
表6
表6回归方程的估计精度分析
Table 6
品种(系) Cultivar (Line) | D值 D value | 回归值 Regression | 拟合误差 Fitting error | 估计精度 Evaluation accuracy (%) |
---|---|---|---|---|
6-2 | 0.7606 | 0.7688 | -0.0081 | 98.93 |
花育17号Huayu 17 | 0.4553 | 0.4243 | 0.0310 | 93.18 |
花育19号Huayu 19 | 0.5407 | 0.5541 | -0.0133 | 97.53 |
花育21号Huayu 21 | 0.6506 | 0.6570 | -0.0064 | 99.02 |
花育33号Huayu 33 | 0.6039 | 0.6181 | -0.0142 | 97.65 |
冀花10号Jihua 10 | 0.5586 | 0.5693 | -0.0107 | 98.08 |
冀花9号Jihua 9 | 0.4381 | 0.4374 | 0.0007 | 99.85 |
闽花6号Minhua 6 | 0.4979 | 0.4912 | 0.0067 | 98.66 |
天府18号Tianfu 18 | 0.4377 | 0.4346 | 0.0031 | 99.30 |
天府22号Tianfu 22 | 0.7222 | 0.7290 | -0.0067 | 99.07 |
天府29号Tianfu 29 | 0.7835 | 0.7663 | 0.0172 | 97.81 |
徐花13号Xuhua 13 | 0.5494 | 0.5482 | 0.0012 | 99.79 |
徐花9号Xuhua 9 | 0.6754 | 0.6868 | -0.0114 | 98.32 |
远杂6号Yuanza 6 | 0.4059 | 0.4081 | -0.0022 | 99.46 |
中花6号Zhonghua 6 | 0.3858 | 0.3790 | 0.0068 | 98.23 |
201060105 | 0.5223 | 0.5177 | 0.0046 | 99.11 |
201060106B1 | 0.7070 | 0.7194 | -0.0124 | 98.24 |
201060401C | 0.4149 | 0.4279 | -0.0130 | 96.87 |
2010605 | 0.4276 | 0.4122 | 0.0153 | 96.41 |
201070406 | 0.2214 | 0.2179 | 0.0035 | 98.43 |
201070407 | 0.3500 | 0.3595 | -0.0095 | 97.29 |
201150104 | 0.4272 | 0.4189 | 0.0082 | 98.07 |
201150114A | 0.3881 | 0.4101 | -0.0220 | 94.33 |
201150118A | 0.5217 | 0.5323 | -0.0106 | 97.98 |
201150119 | 0.4961 | 0.4999 | -0.0038 | 99.23 |
201150201 | 0.7589 | 0.7374 | 0.0216 | 97.16 |
201240404 | 0.4925 | 0.4974 | -0.0049 | 99.00 |
201241001 | 0.7074 | 0.7133 | -0.0058 | 99.17 |
201240401A | 0.6700 | 0.6539 | 0.0161 | 97.59 |
201240413 | 0.6615 | 0.6225 | 0.0391 | 94.10 |
新窗口打开|下载CSV
2.5 各品种(系)耐荫类别的特征分析
结合聚类和逐步回归分析结果,将不同基因型花生在对照和遮荫条件下各主要性状的平均值列于表7中,分析遮荫处理条件下不同耐荫性花生的表现特征:第一类(耐荫型):净光合速率较高,百仁重、单株产量、小区产量和出仁率高;
第二类(中度耐荫型):百仁重、单株产量、小区产量和出仁率均较高,净光合速率高;
第三类(敏感型):百仁重、单株产量、净光合速率、小区产量和出仁率均最低。
Table 7
表7
表7不同耐荫类型各性状的表现特征
Table 7
类别 Cluster | 正常光照CK | 遮荫Shade | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
X9 | X11 | X13 | X12 | X10 | X9 | X11 | X13 | X12 | X10 | ||
第一类 First cluster | 48.59 | 15.05 | 18.15 | 308.18 | 0.62 | 43.43 | 5.89 | 10.45 | 126.89 | 0.54 | |
第二类 Second cluster | 55.40 | 14.92 | 18.22 | 296.86 | 0.63 | 35.17 | 4.87 | 10.91 | 100.05 | 0.53 | |
第三类 Third cluster | 63.38 | 17.29 | 17.12 | 305.41 | 0.63 | 26.80 | 3.80 | 9.34 | 88.55 | 0.49 |
新窗口打开|下载CSV
3 讨论
3.1 花生耐荫性评价分析
花生耐荫性是一个十分复杂的综合性状,由其自身遗传性状、生物和非生物等[20]诸多因素共同决定。花生在适当荫蔽胁迫下具有自我调节适应弱光的能力。在间套作复合种植体系中,除了减少高秆作物对花生的遮荫外,选择耐荫性强的花生品种,是解决间作套种环境下花生植株瘦弱、产量低、品质差的有效措施。因此,开展花生耐荫性研究,找到科学合理评价花生耐荫性的方法,对筛选、挖掘适宜间作套种花生资源以及耐荫花生品种的选育具有重要作用。本研究采用玉米与花生(1﹕1)间作对花生进行自然遮荫处理,选择能直观反映花生生长性状的16个容易测定指标,通过对其耐荫系数及其相关性分析发现,不同花生品种(系)各单项指标的耐荫系数变异幅度不同,除出仁率外各个单项指标间存在显著或极显著的相关性,可见各指标间不仅会发生不同程度的重叠,而且各个指标对遮荫环境的响应程度不同。因此,花生的耐荫性不能采用单一指标进行评价,应该对所测指标进行判断分析后再进行综合评价才更加科学有效[21]。很多****为弄清作物的耐逆性,利用综合评价法对作物的耐荫性、耐寒性、和耐盐性[16,22-23]等进行评价分析。目前,关于作物的耐荫性研究目前较多的是大豆,科学家们采用综合评价法分别从苗期特性[16]、根系[24]、产量、产量构成因素[17,18]、茎秆木质素含量[25,26,27]以及茎秆抗倒性[28,29]等方面对大豆耐荫性进行评价研究。本研究采用多元分析方法,通过主成分分析,提取得到5个主成分,其累计贡献率高达87.127%,具有较强的信息代表性,可将16个单项指标转换为5个彼此独立的综合指标;利用隶属函数分析法获得花生的耐荫综合评价D值;对D值进行聚类分析,将30个花生品种(系)分为3种类型,分别为耐荫型、中度耐荫型和敏感型,其中有11个属于耐荫型,占供试材料的36.67%,有18个属于中度耐荫型,占供试材料的60.00%,而敏感型品种(系)较少,只有1个。可见绝大部分花生属于耐荫和中度耐荫型品种(系),少部分属于敏感型品种(系),说明在遮荫环境下,花生会对一定程度的荫蔽胁迫产生适应性,是一种既喜光又耐荫的作物[30]。
3.2 花生耐荫性鉴定指标选择
光是作物获得能量的来源和形态建成的重要环境因子,在遮荫条件下,作物的形态、光合生理、产量[31]以及品质[32]等特性会发生变化,但并不是所有指标都能有效、准确鉴定花生耐荫性,花生的耐荫性因遮荫环境、作物特性和品种差异等而存在不同。因此,建立可靠的花生耐荫评价体系十分必要。本研究利用逐步回归分析方法建立数学模型D= -0.741+0.576X9+0.507X11+0.298X13+0.272X12+0.406X10,分析可知各个花生品种(系)的估计精度均在93.18%以上,可见,该数学模型所筛选出的5个指标(百仁重、出仁率、单株产量、小区产量和净光合速率)对遮荫环境的响应较为明显,可作为鉴定花生耐荫性的指标。遮荫会对花生正常生长产生影响,有研究已表明,荫蔽导致花生净光合速率、产量和成熟饱满度降低[9,33],而在本研究中,净光合速率、单株产量、百仁重和出仁率均被筛选出来,可见这些指标对弱光胁迫较为敏感,可作为鉴定花生耐荫性的指标。4 结论
将30个花生品种(系)划分为耐荫型花生品种(系)11个、中等耐荫型花生品种(系)18个和敏感型花生品种(系)1个。百仁重、出仁率、单株产量、小区产量和净光合速率可作为鉴定花生品种(系)耐荫性的指标。耐荫型花生在遮荫条件下表现为净光合速率较高,百仁重、单株产量、小区产量和出仁率高。参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
[本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
[本文引用: 2]
[本文引用: 2]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
[本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
URL [本文引用: 3]
URL [本文引用: 3]
,
URL [本文引用: 2]
URL [本文引用: 2]
,
[本文引用: 2]
[本文引用: 2]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]
,
URL [本文引用: 1]
URL [本文引用: 1]