Comprehensive Evaluation of Nitrogen Efficiency and Screening of Varieties with High Grain Yield and High Nitrogen Efficiency of Inbred Middle-Ripe Japonica Rice in the Middle and Lower Reaches of Yangtze River
LIU QiuYuan1,2, ZHOU Lei1, TIAN JinYu1, CHENG Shuang1, TAO Yu1, XING ZhiPeng1, LIU GuoDong1, WEI HaiYan,1, ZHANG HongCheng,1通讯作者:
责任编辑: 杨鑫浩
收稿日期:2020-08-30接受日期:2020-11-23网络出版日期:2021-04-01
基金资助: |
Received:2020-08-30Accepted:2020-11-23Online:2021-04-01
作者简介 About authors
刘秋员,E-mail:
摘要
关键词:
Abstract
Keywords:
PDF (554KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
刘秋员, 周磊, 田晋钰, 程爽, 陶钰, 邢志鹏, 刘国栋, 魏海燕, 张洪程. 长江中下游地区常规中熟粳稻氮效率综合评价及高产氮高效品种筛选[J]. 中国农业科学, 2021, 54(7): 1397-1409 doi:10.3864/j.issn.0578-1752.2021.07.007
LIU QiuYuan, ZHOU Lei, TIAN JinYu, CHENG Shuang, TAO Yu, XING ZhiPeng, LIU GuoDong, WEI HaiYan, ZHANG HongCheng.
开放科学(资源服务)标识码(OSID):
0 引言
【研究意义】产量是水稻生产追求的永恒主题,增加产量是保证粮食安全的重要途径。氮素是决定水稻产量重要的养分因子,同时也是水稻生产肥料成本投入的重要组成部分[1]。由于氮肥的过量投入,所引起的环境污染、生产效益降低等问题,正逐步成为限制水稻持续健康发展的重要障碍[2,3,4,5]。筛选和培育高产氮高效型水稻品种,能够保障产量与氮素吸收利用效率的协同提高,对促进水稻产业健康发展有重要意义。【前人研究进展】科学的筛选评价方法是筛选出高产氮高效类型水稻品种的关键。水稻产量较直观,易于比较筛选,而氮效率涉及的评价指标众多,由于筛选对象与目的的不同,所选用的氮效率评价指标也各有差异,由此也形成了较多的氮效率筛选与评价方法。有****认为在高氮和低氮条件下均高于平均产量水平的水稻品种是氮高效型品种[6,7]。也有****认为可以把氮素农学利用效率[8]、氮肥生理利用效率[9]、氮素稻谷生产效率[10]等作为水稻氮效率评价指标。值得注意的是,不同氮效率评价指标反映的是水稻氮素吸收利用的不同方面,而且同一基因型在不同氮素吸收利用效率评价指标下的排序差异较大[11,12]。因此,利用单个氮效率评价指标对不同基因型的氮效率特性进行评价,并不能全面反映其氮素吸收利用整体能力,可能还会造成一些优良基因型的丢失。氮效率是氮素吸收、同化、运转、再利用等多个生理过程综合作用的结果,较理想的氮高效基因型应该既有较高的氮素利用效率,同时又有较大的吸收效率[13,14]。因此,有必要基于多项氮效率评价指标对其进行综合评价,这有利于筛选结果同时兼顾氮素吸收与利用两个方面。近年来,有研究者也开始探索基于多指标对水稻氮效率进行综合评价。陈琛等[15]基于成熟期吸氮量和产量把114个株系分为了低吸氮低产、低吸氮中产、中吸氮低产、中吸氮中产、中吸氮高产、高吸氮高产6种类型。陈晨等[16]基于隶属函数法和聚类分析,以整株生物量、茎叶生物量、根系生物量、茎叶氮累积量作为评价指标,将55个供试水稻品种分为了氮高效型、氮中效型、氮低效型3个类型。【本研究切入点】自20世纪30年代作物氮效率基因型差异被提出后,尽管前人对水稻氮效率开展了大量研究,但在水稻氮高效品种的筛选方法上仍以单一指标筛选方法为主,缺乏基于多指标对水稻氮效率进行综合评价的研究。当前主流的综合评价方法中,熵权模糊隶属函数法将信息论中熵的概念引入,实现了对影响事物的多个指标的客观赋权,是一种有效的多因素决策方法,在多个领域已被广泛应用[17,18,19,20],但目前尚未见有应用熵权模糊隶属函数法对水稻氮效率进行综合评价的研究报道。【拟解决的关键问题】本研究于2017、2018年分别收集90和105个长江中下游地区常规中熟粳稻品种(系)为材料,基于熵权模糊隶属函数方法,对不同品种(系)的氮素吸收利用效率进行综合评价,不仅可以丰富水稻氮效率评价方法,而且还可以筛选出高产氮高效的粳稻品种(系),为该地区高产氮高效品种的选育与生产应用提供参考。1 材料与方法
1.1 供试材料与种植方法
供试材料为长江中下游粳稻生产地区育成的主要中熟粳稻品种(系),2017年供试品种(系)共90份,2018年在2017年基础上新增15个品种(系),共105份,具体见表1。试验地点位于泰州市姜堰区沈高镇河横村农业生态园。试验田土壤类型为潴育型水稻土,质地黏性。0—20 cm土层含有机质含量31.72 g·kg-1、全氮1.96 g·kg-1、速效钾165.26 mg·kg-1和速效磷32.54 mg·kg-1。Table 1
表1
表1供试品种(系)清单
Table 1
年度 Year | 品种(系)Variety (line) |
---|---|
2017 | JD6602、JD6614、JD6619、常软07-1 Changruan07-1、常软07-11 Changruan07-11、常软07-2 Changruan07-2、常软07-3 Changruan07-3、常软07-4 Changruan07-4、常软07-5 Changruan07-5、常软07-6 Changruan07-6、丰粳1606 Fengjing1606、沪香粳165 Huxiangjing165、沪早软粳 Huzaoruanjing、沪早香软1号 Huzaoxiangruan1、沪早香软2号 Huzaoxiangruan2、华丰1502 Huafeng1502、华粳295 Huajing295、华粳5号 Huajing5、华粳8号 Huajing8、淮330 Huai330、淮稻5号 Huaidao5、连粳11号 Lianjing11、连粳12号 Lianjing12、连粳13 Lianjing13、连粳13264 Lianjing13264、连粳15 Lianjing15、连粳15113 Lianjing15113、连粳7号 Lianjing7、南繁1604 Nanfan1604、南繁1605 Nanfan1605、南繁1609 Nanfan1609、南繁1610 Nanfan1610、南粳2728 Nanjing2728、南粳505 Nanjing505、南粳5711 Nanjing5711、南粳5833 Nanjing5833、南粳9108 Nanjing9108、南粳3818 Nanjing3818、南粳5718 Nanjing5718、宁5720 Ning5720、宁9003 Ning9003、宁9022 Ning9022、宁9039 Ning9039、宁粳040 Ningjing040、宁粳4号Ningjing4、宁粳7号Ningjing7、圣稻1647 Shendao1647、圣稻18-15 Shendao18-15、圣稻18-4 Shendao18-4、圣稻19 Shendao19、圣稻20 Shendao20、圣稻22 Shendao22、圣稻2620 Shendao2620、泗15-234 Si15-234、泗15-301 Si15-301、泗稻14-211 Sidao14-211、泗稻15号 Sidao15、松早香1号 Songzaoxiang1、苏1795 Su1795、苏粳815 Sujing815、苏香粳3号 Suxiangjing3、苏秀867 Suxiu867、泰粳1152 Taijing1152、泰粳2340 Taijing2340、皖垦粳3号 Wankenjing3、武4610 Wu4610、武6267 Wu6267、武粳004 Wujing004、武育粳3号 Wuyujing3、武运5020 Wuyun5020、武运5051 Wuyun5051、武运粳21 Wuyunjing21、武运粳27号 Wuyunjing27、武运粳32号 Wuyunjing32、武运粳80 Wuyunjing80、新稻22 Xindao22、新科稻31 Xinkedao31、徐36618 Xu36618、徐41368 Xu41368、徐稻9号 Xudao9、徐农33202 Xunong33202、盐粳16号Yanjing16、扬粳1612 Yangjing1612、扬粳239 Yangjing239、扬粳3012 Yangjing3012、扬粳3491 Yangjing3491、扬粳5515 Yangjing5515、扬育粳2号Yangyujing2、镇9471 Zhendao9471、镇稻99 Zhendao99 |
2018 (新增 Newly added) | 早香粳1号 Zaoxiangjing1、福粳1601 Fujing1601、福粳1608 Fujing1608、沪早香181 Huzaoxiang181、宁9036 Ning9036、申粳1221 Shenjing1221、圣稻18 Shendao18、圣稻23 Shendao23、圣稻24 Shendao24、圣香66 Shenxiang66、苏1785 Su1785、武育粳36号Wuyujing36、武运4326 Wuyun4326、徐40398 Xu40398、徐稻10号 Xudao10 |
新窗口打开|下载CSV
试验设施氮肥(N18)和不施氮肥(N0)2个处理,其中施氮处理纯氮用量为270 kg·hm-2,按照基肥﹕分蘖肥﹕穗肥=3﹕4﹕3施用,分蘖肥于4叶1心施用,穗肥于倒4叶期施用,氮﹕磷﹕钾比例为2﹕1﹕2,磷肥一次性基施,钾肥分别于耕翻前、拔节期等量施入。不施氮肥处理的磷、钾肥用量及施用方法同施氮处理。2年均于6月13日播种,6月23日至24日移栽,行距25 cm,穴距6 cm,每穴2苗,小区面积15 m2(3 m×5 m),重复2次。水分管理及病虫草害防治等相关栽培措施均按照当地高产栽培技术要求进行。
1.2 测定内容与方法
1.2.1 产量及其构成因素 在成熟期,各供试品种(系)每个小区普查60穴茎糵数,计算有效穗数。每个小区按照每穴平均有效穗数选取长势相近植株5穴,调查每穗粒数、结实率和千粒重。在每个小区中间割8 m2(2 m×4 m)进行实际产量测定。1.2.2 干物质及氮素积累测定 在成熟期,各供试品种(系)按照每穴平均茎糵数选取长势相近植株3穴,分茎鞘、叶、穗置于105℃杀青30 min,80℃烘至恒重后称重,计算各供试品种(系)单位面积茎鞘、叶、穗干物质积累量及总干物质积累量(茎鞘、叶、穗干物质积累量之和)。将各器官干物质样品粉碎,采用FOSS凯氏定氮仪测定茎、叶、穗各器官含氮率。按照各器官干物质量×各器官含氮率计算成熟期茎、叶、穗各器官氮素吸收量,总吸氮量为茎、叶、穗3个器官氮素吸收量之和。
1.3 数据处理与统计方法
1.3.1 氮素吸收利用效率评价指标计算方法氮肥农学利用率(nitrogen agronomic efficiency,NAE,kg·kg-1)=(施氮区籽粒产量-无氮区籽粒产量)/氮肥施用量;
氮素生理利用率(nitrogen physiological efficiency,NPE,kg·kg-1)=(施氮区籽粒产量-无氮区籽粒产量)/(施氮区总吸氮量-无氮区总吸氮量);
氮肥回收效率(nitrogen recovery efficiency,NRE,%)=(施氮区总吸氮量-无氮区总吸氮量)/氮肥施用量×100;
氮素干物质生产效率(nitrogen dry matter production efficiency,DME,kg·kg-1)=施氮区植株干物质积累量/植株总吸氮量;
氮素籽粒生产效率(nitrogen grain yield production efficiency,GYE,kg·kg-1)=施氮区水稻籽粒产量/植株总吸氮量。
1.3.2 氮效率综合评价方法 参照张先起等[18]、余健等[20]提供的方法,具体步骤为:
(1)构建从优隶属度矩阵 按照从优隶属度原则,各单项评价指标相应的模糊量值隶属于标准方案对应评价指标相应模糊量值的程度称为从优隶属度。本研究各评价指标按照越大越优型指标计算。rij=$\frac{x_{ij}}{max x_{ij}}$,式中rij为从优隶属度,maxxij为各评价指标的最大值。
(2)熵权法确定指标权重 首先原始数据归一化处理。按照越大越优型指标计算处理,计算公式为$b_{ij}=\frac{x_{ij}-min x_{ij}}{max x_{ij}-min x_{ij}}$,式中,max xij和min xij分别为各评价指标的最大值和最小值。其次计算熵值。m个评价事物n个评价指标的熵计算公式为Hi=$- \frac{\sum{j=1}{m}f_{ij}In f_{ij}}{In m}$,式中i=1,2,···,n,j=1,2,···,m,fij=$\frac{1+b_{ij}}{\sum{j-1}{m}(1+b_{ij})}$。最后计算评价指标权重W,计算公式为$W_{i}=\frac{1-H_{i}}{n-\sum{i-1}{n}H_{i}}$,且满足$\sum{i=1}{n}W_{i}$=1。
(3)计算模糊评价的综合指数 把各评价指标的
权重值wi与相应评价指标的从优隶属度值相乘并累加,可得模糊评价的综合指数值Z,计算公式为Zj=$\sum{i-1}{n}w_{i}r_{ij}$。该值越大,说明该品种综合评价越优。
采用 Microsoft Excel 2013整理数据和绘制图表,SPSS 22.0进行相关统计分析。
2 结果
2.1 供试品种(系)产量、干物质、氮素吸收与利用效率
施氮肥处理(N18)下,供试品种(系)2年产量变化范围为6.80—10.28 t·hm-2,平均产量为8.50 t·hm-2;总干物质积累量和总吸氮量变化范围为14.56— 21.25 t·hm-2和175.72—270.38 kg·hm-2,平均值分别为18.27 t·hm-2和225.53 kg·hm-2(表 2)。不管是施氮肥处理还是不施氮肥处理,产量构成因素中均是以穗粒数的变异系数最大,各器官的干物质积累量则是以叶的变异系数最大,2年规律一致。此外,与不施氮肥处理相比,施氮肥处理增大了穗粒数、产量、叶干重、穗干重、总干重等指标在品种间的变异系数,而降低了总吸氮量的变异系数,说明供试品种(系)的产量以及干物质生产对氮肥的响应程度存在差异。氮素吸收利用效率指标除了氮素籽粒生产效率和氮素干物质生产效率变异系数较小(2年均在5%以下)外,氮肥回收效率、氮肥农学利用率、氮素生理利用率的变异系数均在10%以上,其中氮肥农学利用率的变异系数2年均在30%以上,由此也说明供试品种(系)在氮素吸收与利用效率上存在较大差异。Table 2
表2
表2供试品种(系)产量、干物质、氮素吸收与利用性状
Table 2
年份 Year | 指标 Trait | N18 | N0 | ||||
---|---|---|---|---|---|---|---|
变幅 Range | 平均值 Mean | 变异系数 CV (%) | 变幅 Range | 平均值 Mean | 变异系数 CV (%) | ||
2017 | 有效穗数 Number of effective panicles (×104 hm-2) | 302.58-467.52 | 365.20 | 8.16 | 230.10-426.90 | 303.97 | 9.31 |
千粒重 1000-grain weight (g) | 20.60-29.88 | 26.38 | 5.87 | 20.93-30.25 | 26.88 | 5.75 | |
穗粒数 Spikelet per panicle | 75.38-136.08 | 103.21 | 10.48 | 63.88-114.3 | 94.47 | 10.05 | |
结实率 Percentage of filled grains (%) | 82.43-95.77 | 89.43 | 3.39 | 79.34-98.43 | 92.37 | 4.27 | |
产量 Grain yield (t·hm-2) | 6.80-9.89 | 8.42 | 8.11 | 5.25-7.76 | 6.69 | 6.91 | |
茎干重 SDM (t·hm-2) | 4.62-7.44 | 6.01 | 9.29 | 3.36-5.30 | 4.36 | 7.82 | |
叶干重 LDM (t·hm-2) | 2.24-4.37 | 3.28 | 12.34 | 1.66-2.99 | 2.39 | 10.72 | |
穗干重 PDM (t·hm-2) | 6.94-10.99 | 9.00 | 8.93 | 5.41-8.36 | 7.13 | 7.38 | |
总干重 TDM (t·hm-2) | 14.56-21.25 | 18.29 | 7.89 | 10.75-15.89 | 13.89 | 6.65 | |
茎吸氮量 SNC (kg·hm-2) | 43.92-84.95 | 61.50 | 12.71 | 20.28-37.56 | 25.29 | 11.40 | |
叶吸氮量 LNC (kg·hm-2) | 36.62-70.98 | 52.00 | 12.57 | 21.60-41.75 | 31.84 | 13.11 | |
穗吸氮量 PNC (kg·hm-2) | 89.19-149.51 | 121.92 | 8.60 | 61.64-98.58 | 83.74 | 8.73 | |
总吸氮量 TNC (kg·hm-2) | 194.77-270.38 | 235.42 | 6.50 | 117.17-165.88 | 140.87 | 7.53 | |
氮肥回收效率 NRE (%) | 24.55-47.83 | 35.02 | 13.73 | — | — | — | |
氮肥农学利用率 NAE (kg·kg-1) | 1.70-13.58 | 6.39 | 35.17 | — | — | — | |
氮素生理利用率 NPE (kg·kg-1) | 6.38-32.16 | 17.92 | 25.81 | — | — | — | |
氮素籽粒生产效率 GYE (kg·kg-1) | 31.52-39.38 | 35.75 | 3.97 | — | — | — | |
氮素干物质生产效率 DME (kg·kg-1) | 69.24-83.74 | 77.66 | 3.54 | — | — | — | |
2018 | 有效穗数 Number of effective panicles (×104 hm-2) | 268.51-456.36 | 350.82 | 8.58 | 239.04-403.73 | 288.24 | 8.80 |
千粒重 1000-grain weight (g) | 20.78-30.17 | 26.62 | 6.57 | 22.05-30.88 | 27.86 | 6.19 | |
穗粒数 Spikelet per panicle | 74.19-134.88 | 103.31 | 11.08 | 67.22-119.13 | 91.33 | 9.85 | |
结实率 Percentage of filled grains (%) | 86.84-98.56 | 94.60 | 2.62 | 91.69-98.57 | 96.42 | 1.59 | |
产量 Grain yield (t·hm-2) | 6.82-10.28 | 8.58 | 8.49 | 5.04-8.01 | 6.60 | 8.03 | |
茎干重 SDM (t·hm-2) | 4.78-7.09 | 5.99 | 8.34 | 3.39-5.46 | 4.32 | 8.95 | |
叶干重 LDM (t·hm-2) | 2.15-3.83 | 3.05 | 11.43 | 1.45-2.43 | 1.98 | 11.08 | |
穗干重 PDM (t·hm-2) | 7.20-11.63 | 9.21 | 9.72 | 5.08-8.70 | 7.06 | 8.62 | |
总干重 TDM (t·hm-2) | 14.77-21.23 | 18.25 | 7.94 | 10.84-15.93 | 13.36 | 7.40 | |
茎吸氮量 SNC (kg·hm-2) | 37.94-73.85 | 55.97 | 11.60 | 15.61-30.53 | 21.67 | 13.71 | |
叶吸氮量 LNC (kg·hm-2) | 25.96-58.01 | 42.33 | 13.61 | 14.79-28.33 | 21.09 | 15.84 | |
穗吸氮量 PNC (kg·hm-2) | 89.37-144.78 | 117.32 | 9.73 | 61.02-95.06 | 78.76 | 8.27 | |
总吸氮量 TNC (kg·hm-2) | 175.72-246.13 | 215.63 | 6.80 | 100.01-146.65 | 121.51 | 7.32 | |
氮肥回收效率 NRE (%) | 21.06-45.79 | 34.86 | 15.17 | — | — | — | |
氮肥农学利用率 NAE (kg·kg-1) | 2.27-13.97 | 7.33 | 32.29 | — | — | — | |
氮素生理利用率 NPE (kg·kg-1) | 10.76-30.96 | 20.65 | 22.31 | — | — | — | |
氮素籽粒生产效率 GYE (kg·kg-1) | 36.28-44.26 | 39.77 | 4.27 | — | — | — | |
氮素干物质生产效率 DME (kg·kg-1) | 77.15-91.93 | 84.59 | 3.60 | — | — | — |
新窗口打开|下载CSV
2.2 氮素吸收利用效率综合评价
不同氮素吸收利用效率评价指标反映的是氮素吸收与利用效率的不同方面,其中氮肥回收效率、氮肥农学利用率反映的是不同水稻品种对施入氮肥的吸收和增产效率,氮素生理利用率、氮素籽粒生产效率、氮素干物质生产效率反映的是植株吸收的氮素形成生物产量和经济产量能力。本研究以该5项指标,应用熵权模糊隶属函数法对供试品种(系)的氮素吸收与利用效率进行综合评价。2017年权重最大的指标为氮肥回收效率,最小的是氮素籽粒生产效率,2018年权重最大的指标为氮素生理利用率,最小的是氮素干物质生产效率(表3)。各供试品种(系)氮效率综合得分可知,2年得分范围分别为0.485—0.952(2017)和0.535—0.977(2018),平均得分0.707(2017)和0.752(2018)(表4)。Table 3
表3
表3氮素吸收利用效率评价指标权重
Table 3
年份 Year | 氮肥回收效率 NRE | 氮肥农学利用率 NAE | 氮素生理利用率 NPE | 氮素籽粒生产效率 GYE | 氮素干物质生产效率 DME |
---|---|---|---|---|---|
2017 | 0.247 | 0.219 | 0.186 | 0.169 | 0.178 |
2018 | 0.186 | 0.192 | 0.230 | 0.211 | 0.182 |
新窗口打开|下载CSV
Table 4
表4
表4供试品种(系)的氮效率综合值和产量-氮效率综合指数
Table 4
指标 Trait | 年份 Year | 变幅 Range | 变异系数 CV (%) | 平均值 Mean |
---|---|---|---|---|
氮效率综合值CEV | 2017 | 0.485-0.952 | 12.28 | 0.707 |
2018 | 0.535-0.977 | 11.63 | 0.752 | |
产量-氮效率综合指数GYNEI | 2017 | 0.569-1.583 | 19.45 | 1.008 |
2018 | 0.628-1.541 | 18.80 | 1.007 |
新窗口打开|下载CSV
2.3 高产氮高效品种(系)筛选
本研究在参照FAGERIA等[21]提供的产量效率指数计算公式的基础上进行了改进,改为产量-氮效率综合指数(grain yield-nitrogen efficiency index, GYNEI),GYNEI=(产量/群体平均产量)×(氮效率综合指数/群体氮效率综合指数平均值)。供试品种(系)GYNEI 2017年变化范围为0.569—1.583,平均值为1.008,2018年变化范围为0.628—1.541,平均值为1.007(表 4)。采用系统聚类方法(组内联结-平方欧式距离),基于GYNEI把供试品种(系)分为3类(表5)。根据分类结果,供试品种(系)类型Ⅰ 2017年有23个,2018年有27个,其中2年相同的品种(系)有19个;类型Ⅱ 2017年有42个,2018年有35个,其中2年相同的品种(系)有23个;类型Ⅲ 2017年有25个,2018年有43个,其中2年相同的品种(系)有24个。对不同类型的产量以及氮效率指标进行统计分析(表6),GYNEI大小依次为类型Ⅰ>类型Ⅱ>类型Ⅲ,其中产量、氮肥回收效率、氮肥农学利用率、氮素生理利用率、氮素籽粒生产效率、氮素干物质生产效率在3个类型间的大小顺序也表现为类型Ⅰ>类型Ⅱ>类型Ⅲ,且类型Ⅰ与类型Ⅱ、类型Ⅲ之间的差异均达到了显著水平,说明类型Ⅰ为高产氮高效类型,类型Ⅲ为低产氮低效类型。
Table 5
表5
表5供试品种的产量-氮效率类型
Table 5
品种(系) Variety (line) | 产量-氮效率类型 GYNET | 品种(系) Variety (line) | 产量-氮效率类型 GYNET | 品种(系) Variety (line) | 产量-氮效率类型 GYNET | |||||
---|---|---|---|---|---|---|---|---|---|---|
2017 | 2018 | 2017 | 2018 | 2017 | 2018 | |||||
常软07-3 Changruan07-3 | Ⅰ | Ⅰ | 连粳13264 Lianjing13264 | Ⅱ | Ⅲ | 沪早软粳Huzaoruanjing | Ⅲ | Ⅱ | ||
华丰1502 Huafeng1502 | Ⅰ | Ⅱ | 连粳15113 Lianjing15113 | Ⅱ | Ⅱ | 沪早香软1号Huzaoxiangruan1 | Ⅲ | Ⅲ | ||
华粳8号Huajing8 | Ⅰ | Ⅰ | 南繁1610 Nanfan1610 | Ⅱ | Ⅲ | 华粳295 Huajing295 | Ⅲ | Ⅲ | ||
南繁1605 Nanfan1605 | Ⅰ | Ⅰ | 南粳505 Nanjing505 | Ⅱ | Ⅲ | 华粳5号Huajing5 | Ⅲ | Ⅲ | ||
南粳5711 Nanjing5711 | Ⅰ | Ⅰ | 南粳5833 Nanjing5833 | Ⅱ | Ⅱ | 连粳15 Lianjing15 | Ⅲ | Ⅲ | ||
南粳9108 Nanjing9108 | Ⅰ | Ⅰ | 南粳3818 Nanjing3818 | Ⅱ | Ⅱ | 连粳7号 Lianjing7 | Ⅲ | Ⅲ | ||
南粳5718 Nanjing5718 | Ⅰ | Ⅰ | 宁5720 Ning5720 | Ⅱ | Ⅲ | 南繁1604 Nanfan1604 | Ⅲ | Ⅲ | ||
宁粳040 Ningjing040 | Ⅰ | Ⅰ | 宁9003 Ning9003 | Ⅱ | Ⅱ | 南繁1609 Nanfan1609 | Ⅲ | Ⅲ | ||
宁粳4号Ningjing4 | Ⅰ | Ⅰ | 宁9022 Ning9022 | Ⅱ | Ⅱ | 南粳2728 Nanjing2728 | Ⅲ | Ⅲ | ||
宁粳7号Ningjing7 | Ⅰ | Ⅰ | 宁9039 Ning9039 | Ⅱ | Ⅰ | 圣稻18-4 Shendao18-4 | Ⅲ | Ⅲ | ||
泗稻15号Sidao15 | Ⅰ | Ⅰ | 圣稻1647 Shendao1647 | Ⅱ | Ⅱ | 圣稻2620 Shendao2620 | Ⅲ | Ⅲ | ||
苏1795 Su1795 | Ⅰ | Ⅰ | 圣稻18-15 Shendao18-15 | Ⅱ | Ⅱ | 泗稻14-211 Sidao14-211 | Ⅲ | Ⅲ | ||
泰粳1152 Taijing1152 | Ⅰ | Ⅰ | 圣稻19 Shendao19 | Ⅱ | Ⅱ | 苏粳815 Sujing815 | Ⅲ | Ⅲ | ||
泰粳2340 Taijing2340 | Ⅰ | Ⅰ | 圣稻20 Shendao20 | Ⅱ | Ⅱ | 苏香粳3号 Suxiangjing3 | Ⅲ | Ⅲ | ||
武粳004 Wujing004 | Ⅰ | Ⅱ | 圣稻22 Shendao22 | Ⅱ | Ⅱ | 苏秀867 Suxiu867 | Ⅲ | Ⅲ | ||
武运5020 Wuyun5020 | Ⅰ | Ⅰ | 泗15-234 Si15-234 | Ⅱ | Ⅲ | 武育粳3号 Wuyujing3 | Ⅲ | Ⅲ | ||
武运5051 Wuyun5051 | Ⅰ | Ⅱ | 泗15-301 Si15-301 | Ⅱ | Ⅱ | 武运粳27号 Wuyunjing27 | Ⅲ | Ⅲ | ||
徐36618 Xu36618 | Ⅰ | Ⅱ | 松早香1号 Songzaoxiang1 | Ⅱ | Ⅲ | 新稻22 Xindao22 | Ⅲ | Ⅲ | ||
徐41368 Xu41368 | Ⅰ | Ⅰ | 皖垦粳3号 Wankenjing3 | Ⅱ | Ⅱ | 新科稻31 Xinkedao31 | Ⅲ | Ⅲ | ||
扬粳1612 Yangjing1612 | Ⅰ | Ⅰ | 武4610 Wu4610 | Ⅱ | Ⅱ | 徐稻9号Xudao9 | Ⅲ | Ⅲ | ||
扬粳239 Yangjing239 | Ⅰ | Ⅰ | 武6267 Wu6267 | Ⅱ | Ⅱ | 福粳1601 Fujing1601 | — | Ⅱ | ||
扬粳3491 Yangjing3491 | Ⅰ | Ⅰ | 武运粳21 Wuyunjing21 | Ⅱ | Ⅲ | 福粳1608 Fujing1608 | — | Ⅱ | ||
扬粳5515 Yangjing5515 | Ⅰ | Ⅰ | 武运粳32号 Wuyunjing32 | Ⅱ | Ⅱ | 沪早香181 Huzaoxiang181 | — | Ⅱ | ||
JD6602 | Ⅱ | Ⅱ | 武运粳80 Wuyunjing80 | Ⅱ | Ⅰ | 宁9036 Ning9036 | — | Ⅲ | ||
常软07-1 Changruan07-1 | Ⅱ | Ⅰ | 徐农33202 Xunong33202 | Ⅱ | Ⅲ | 申粳1221 Shenjing1221 | — | Ⅲ | ||
常软07-2 Changruan07-2 | Ⅱ | Ⅱ | 盐粳16号 Yanjing16 | Ⅱ | Ⅲ | 圣稻18 Shendao18 | — | Ⅲ | ||
常软07-4 Changruan07-4 | Ⅱ | Ⅲ | 扬粳3012 Yangjing3012 | Ⅱ | Ⅱ | 圣稻23 Shendao23 | — | Ⅰ | ||
常软07-6 Changruan07-6 | Ⅱ | Ⅰ | 扬育粳2号 Yangyujing2 | Ⅱ | Ⅱ | 圣稻24 Shendao24 | — | Ⅲ | ||
丰粳1606 Fengjing1606 | Ⅱ | Ⅱ | 镇9471 Zhendao9471 | Ⅱ | Ⅱ | 圣香66 Shenxiang66 | — | Ⅱ | ||
沪早香软2号 Huzaoxiangruan2 | Ⅱ | Ⅱ | 镇稻99 Zhendao99 | Ⅱ | Ⅲ | 苏1785 Su1785 | — | Ⅰ | ||
淮330 Huai330 | Ⅱ | Ⅱ | JD6614 | Ⅲ | Ⅲ | 武育粳36号 Wuyujin36 | — | Ⅲ | ||
淮稻5号 Huaidao5 | Ⅱ | Ⅱ | JD6619 | Ⅲ | Ⅲ | 武运4326 Wuyun4326 | — | Ⅰ | ||
连粳11号 Lianjing11 | Ⅱ | Ⅱ | 常软07-11 Changruan07-11 | Ⅲ | Ⅲ | 徐40398 Xu40398 | — | Ⅱ | ||
连粳12号 Lianjing12 | Ⅱ | Ⅲ | 常软07-5 Changruan07-5 | Ⅲ | Ⅲ | 徐稻10号Xudao10 | — | Ⅰ | ||
连粳13 Lianjing13 | Ⅱ | Ⅲ | 沪香粳165 Huxiangjing165 | Ⅲ | Ⅲ | 早香粳1号 Zaoxiangjing1 | — | Ⅲ |
新窗口打开|下载CSV
Table 6
表6
表6不同产量-氮效率类型的产量和氮素吸收利用效率
Table 6
年份 Year | 产量-氮 效率类型GYNET | 指标 Trait | 产量-氮效率 综合指数 GYNEI | 产量 Grain yield (t·hm-2) | 氮肥回收效率 NRE (%) | 氮肥农学 利用率 NAE (kg·kg-1) | 氮素生理 利用率 NPE (kg·kg-1) | 氮素籽粒 生产效率GYE (kg·kg-1) | 氮素干物质 生产效率 DME (kg·kg-1) |
---|---|---|---|---|---|---|---|---|---|
2017 | Ⅰ | 变幅 Range | 1.142-1.583 | 8.51-9.89 | 30.65-47.83 | 7.41-13.58 | 18.53-32.16 | 35.40-39.38 | 76.91-83.74 |
变异系数 CV (%) | 9.36 | 3.88 | 10.46 | 18.56 | 13.46 | 2.91 | 2.59 | ||
平均值 Mean | 1.271a | 9.27a | 39.18a | 9.15a | 23.33a | 37.19a | 79.97a | ||
Ⅱ | 变幅Range | 0.882-1.107 | 7.41-9.16 | 26.90-42.48 | 4.37-8.52 | 12.60-24.42 | 33.75-37.56 | 72.97-81.41 | |
变异系数 CV (%) | 7.09 | 4.22 | 9.88 | 17.95 | 15.01 | 2.70 | 2.55 | ||
平均值 Mean | 0.993b | 8.38b | 35.51b | 6.30b | 17.74b | 35.47b | 77.32b | ||
Ⅲ | 变幅 Range | 0.569-0.867 | 6.80-8.32 | 24.55-36.85 | 1.70-5.27 | 6.38-17.72 | 31.52-36.93 | 69.24-81.52 | |
变异系数 CV (%) | 8.84 | 5.25 | 10.21 | 21.01 | 20.59 | 3.97 | 4.08 | ||
平均值 Mean | 0.791c | 7.70c | 30.36c | 4.01c | 13.25c | 34.88b | 76.10b | ||
2018 | Ⅰ | 变幅 Range | 1.156-1.541 | 8.90-10.28 | 30.10-45.79 | 7.71-13.97 | 19.83-30.96 | 39.43-44.26 | 83.75-91.58 |
变异系数 CV (%) | 7.16 | 3.59 | 9.78 | 16.09 | 11.49 | 2.47 | 1.81 | ||
平均值 Mean | 1.256a | 9.50a | 39.45a | 9.95a | 25.20a | 41.30a | 87.02a | ||
Ⅱ | 变幅 Range | 0.940-1.130 | 7.47-9.51 | 26.92-41.97 | 5.76-10.12 | 15.83-26.32 | 36.28-42.81 | 79.93-88.64 | |
变异系数 CV (%) | 5.54 | 4.19 | 9.96 | 13.69 | 11.54 | 3.57 | 2.89 | ||
平均值 Mean | 1.042b | 8.65b | 36.51b | 8.02b | 22.02b | 39.69b | 84.35b | ||
Ⅲ | 变幅 Range | 0.628-0.924 | 6.82-8.71 | 21.06-35.58 | 2.27-7.31 | 10.76-22.47 | 36.53-43.13 | 77.15-91.93 | |
变异系数 CV (%) | 9.12 | 5.16 | 12.71 | 24.33 | 19.98 | 4.09 | 3.98 | ||
平均值 Mean | 0.823c | 7.94c | 30.62c | 5.13c | 16.68c | 38.87b | 83.27b |
新窗口打开|下载CSV
2.4 不同产量-氮效率类型的产量构成因素
穗粒数在3个类型间的大小顺序均表现为类型Ⅰ>类型Ⅱ>类型Ⅲ,且类型间的差异均达到显著水平,2年规律一致。有效穗数、千粒重以及结实率在3个类型间的差异均没有达到显著水平,2年规律一致(表7)。由此可见,较高的穗粒数能够促进产量与氮效率的协同提高。Table 7
表7
表7不同产量-氮效率类型的产量构成因素
Table 7
年份 Year | 产量-氮效率类型 GYNET | 指标 Trait | 有效穗数 Number of effective panicles (×104 hm-2) | 千粒重 1000-grain weight (g) | 穗粒数 Spikelet per panicle | 结实率 Percentage of filled grains (%) |
---|---|---|---|---|---|---|
2017 | Ⅰ | 变幅 Range | 313.50-410.25 | 24.15-29.88 | 94.55-136.08 | 83.55-93.44 |
变异系数 CV (%) | 7.21 | 5.63 | 8.85 | 3.21 | ||
平均值 Mean | 363.20a | 26.81a | 111.18a | 88.77a | ||
Ⅱ | 变幅 Range | 317.34-429.58 | 24.00-29.13 | 83.32-128.31 | 82.43-95.61 | |
变异系数 CV (%) | 6.61 | 5.01 | 9.48 | 3.25 | ||
平均值 Mean | 369.49a | 26.46a | 102.36b | 88.99a | ||
Ⅲ | 变幅 Range | 302.58-467.52 | 20.60-29.38 | 75.38-112.30 | 84.52-95.77 | |
变异系数 CV (%) | 11.03 | 7.04 | 9.54 | 3.43 | ||
平均值 Mean | 359.82a | 25.86a | 97.31c | 90.80a | ||
2018 | Ⅰ | 变幅 Range | 286.46-386.83 | 24.15-30.17 | 95.95-134.88 | 88.93-97.03 |
变异系数 CV (%) | 7.84 | 5.91 | 8.99 | 2.43 | ||
平均值 Mean | 347.67a | 27.09a | 112.31a | 94.31a | ||
Ⅱ | 变幅 Range | 308.83-425.53 | 23.98-30.06 | 90.77-119.47 | 90.22-97.92 | |
变异系数 CV (%) | 7.46 | 6.83 | 8.42 | 2.48 | ||
平均值 Mean | 350.85a | 26.69a | 103.94b | 94.42a | ||
Ⅲ | 变幅 Range | 268.51-456.36 | 20.78-29.65 | 74.19-126.02 | 86.84-98.56 | |
变异系数 CV (%) | 9.88 | 6.61 | 10.70 | 2.85 | ||
平均值 Mean | 352.78a | 26.26a | 97.14c | 94.93a |
新窗口打开|下载CSV
2.5 不同产量-氮效率类型的干物质及氮素积累
茎干重、叶干重、穗干重以及总干重均表现为类型Ⅰ>类型Ⅱ>类型Ⅲ,且类型间的差异均达到显著水平,2年规律一致。此外,类型Ⅲ总干重的最大值分别为18.31 t·hm-2(2017)和18.50 t·hm-2(2018),而类型Ⅰ的总干重的最小值分别为18.56 t·hm-2(2017)和18.63 t·hm-2(2018),说明类型Ⅰ的品种(系)的总干重积累量均要高于类型Ⅲ的品种(系)。茎吸氮量、叶吸氮量、穗吸氮量以及总吸氮量也都表现出类型Ⅰ>类型Ⅱ>类型Ⅲ的规律,但只有穗吸氮量和总吸氮量在3个类型间差异显著(表8)。Table 8
表8
表8不同产量-氮效率类型的干物质和氮素积累
Table 8
年份 Year | 产量-氮 效率类型 GYNET | 指标 Trait | 茎干重 SDM (t·hm-2) | 叶干重 LDM (t·hm-2) | 穗干重 PDM (t·hm-2) | 总干重 TDM (t·hm-2) | 茎吸氮量 SNC (kg·hm-2) | 叶吸氮量 LNC (kg·hm-2) | 穗吸氮量 PNC (kg·hm-2) | 总吸氮量 TNC (kg·hm-2) |
---|---|---|---|---|---|---|---|---|---|---|
2017 | Ⅰ | 变幅 Range | 5.63-7.44 | 2.96-4.37 | 8.68-10.99 | 18.56-21.25 | 52.45-84.95 | 43.89-67.41 | 113.42-149.51 | 227.38-270.38 |
变异系数 CV (%) | 7.33 | 9.14 | 6.01 | 3.98 | 12.20 | 11.69 | 7.38 | 4.72 | ||
平均值 Mean | 6.44a | 3.54a | 9.95a | 19.93a | 64.35a | 53.48a | 131.60a | 249.42a | ||
Ⅱ | 变幅 Range | 4.79-6.97 | 2.49-4.14 | 7.35-9.76 | 15.28-20.11 | 49.90-76.20 | 36.62-67.06 | 99.58-135.63 | 207.82-253.72 | |
变异系数 CV (%) | 7.41 | 10.16 | 5.54 | 4.72 | 11.22 | 12.47 | 5.97 | 4.44 | ||
平均值 Mean | 6.06b | 3.32b | 8.90b | 18.27b | 63.44a | 52.04a | 120.88b | 236.37b | ||
Ⅲ | 变幅 Range | 4.62-6.28 | 2.24-3.89 | 6.93-9.02 | 14.56-18.31 | 43.92-66.82 | 38.06-70.98 | 89.19-129.54 | 194.77-242.98 | |
变异系数 CV (%) | 7.93 | 13.00 | 6.08 | 5.94 | 10.31 | 13.48 | 8.14 | 5.44 | ||
平均值 Mean | 5.54c | 2.97c | 8.30c | 16.80c | 55.60b | 50.58a | 114.75c | 220.93c | ||
2018 | Ⅰ | 变幅 Range | 5.61-7.11 | 3.02-3.83 | 9.29-11.63 | 18.63-21.23 | 46.63-68.80 | 37.03-55.14 | 108.10-144.78 | 208.18-246.13 |
变异系数 CV (%) | 5.78 | 6.15 | 5.16 | 3.18 | 9.65 | 8.92 | 7.37 | 3.79 | ||
平均值 Mean | 6.37a | 3.33a | 10.32a | 20.02a | 57.48 a | 44.64 a | 128.00 a | 230.12 a | ||
Ⅱ | 变幅 Range | 5.16-6.94 | 2.42-3.83 | 7.78-10.08 | 16.54-20.12 | 44.90-73.85 | 35.27-58.01 | 99.49-133.55 | 193.36-239.16 | |
变异系数 CV (%) | 7.68 | 11.09 | 5.16 | 4.49 | 11.62 | 15.57 | 7.68 | 4.93 | ||
平均值 Mean | 6.03 b | 3.10 b | 9.25 b | 18.39 b | 57.20 a | 43.44 a | 117.48 b | 218.11 b | ||
Ⅲ | 变幅 Range | 4.78-6.64 | 2.15-3.30 | 7.20-9.73 | 14.77-18.50 | 37.95-66.34 | 25.96-49.08 | 89.37-125.45 | 175.72-220.35 | |
变异系数 CV (%) | 7.67 | 9.89 | 6.44 | 5.30 | 12.13 | 12.56 | 8.14 | 5.46 | ||
平均值 Mean | 5.70 c | 2.83 c | 8.48 c | 17.01 c | 54.03 a | 39.97 b | 110.49 c | 204.50 c |
新窗口打开|下载CSV
3 讨论
3.1 关于水稻氮效率的综合评价方法
作物氮效率是氮素吸收、同化、运转、再利用等多个生理过程综合作用的结果,其中氮素吸收效率与氮素物质生产效率之间还存在一定的相互抑制关系[22]。因此,单一的氮效率评价指标难以全面反映不同水稻品种的整体氮效率水平,需要基于多项氮效率评价指标才能对不同水稻品种的氮效率进行准确的评价。前人多采用模糊隶属函数法、主成分分析法等综合评价方法对作物品种进行筛选和评价[19, 23-24]。应用主成分分析法进行多指标综合评价时,其优点是可以消除评价指标间的相关影响,有助于保证客观性,但其缺点是所提取主成分的方差贡献率一般很难达到100%,不能完整的代表评价对象的所有信息[25]。模糊隶属函数法虽然不存在上述问题,但如何确定其评价指标的权重是保证其评价结果合理性的关键[23]。本研究在模糊隶属函数法的基础上,引入熵权理论,从评价指标本身提供的有效信息量的多寡程度来计算指标权重[18],保证了评价指标赋权的客观性。本研究得到的5项氮效率评价指标的权重大小在年际间并不一致,这可能一方面与2年供试品种(系)的样本量不同有关,另一方面可能与不同氮效率指标在年际间的稳定程度不同有关[11]。由此说明熵权法能够充分依据评价指标的变化规律来合理确定指标权重,从而保证其赋权的客观性。指标权重以及样本量在年际间存在差异,可能是导致一些品种(系)的产量-氮效率类型在年度间存在差异的主要原因。高产条件下的氮高效水稻品种筛选对缓解粮食需求、环境恶化等面临的压力具有现实意义。因此,本研究在得到了各供试品种(系)的氮效率综合值后,然后根据氮效率综合值和产量计算了产量-氮效率综合指数,并采用系统聚类方法基于产量-氮效率综合指数对供试品种(系)进行了类型划分。分类结果表明,高产氮高效类型(类型Ⅰ)的产量、氮肥回收效率、氮肥农学利用率、氮素生理利用率、氮素籽粒生产效率、氮素干物质生产效率的平均值均要显著高于低产氮低效类型(类型Ⅲ),说明本研究所应用的筛选方法有效可行。我们注意到不同产量-氮效率类型的产量、氮素吸收利用效率、干物质积累及吸氮量等指标的变异范围之间均存在不同程度的交叉,这可能一方面与不同氮效率评价指标反映的是水稻氮素吸收与利用效率的不同方面,且同一品种在不同氮素吸收与利用效率指标上的排序存在较大差异有关[11];另一方面可能与本研究是基于产量和多项氮素吸收与利用效率指标的综合值进行水稻品种类型的划分,且产量与氮素吸收利用效率之间也并不是简单的线性相关关系有关[26,27]。由此可见,采用单一氮效率评价指标对水稻品种进行筛选,可能会导致一些具有其他优良性状的基因型丢失。虽然熵权模糊隶属函数法在水稻氮效率评价领域的应用尚属尝试,还存在入选的评价指标是否合理等问题需要进一步研究探讨,但我们认为这是一种有益的尝试,可以进一步丰富高产氮高效水稻品种筛选评价方法。3.2 关于水稻产量与氮效率协同的途径
实现水稻高产氮高效,对提高水稻产量和氮肥利用效率,促进水稻产业持续健康发展具有重要意义。前人研究表明,氮高效的水稻品种,其产量、生物量以及氮素吸收量均要显著高于氮低效品种[9, 28-29]。本研究也有一致的发现,随着产量-氮效率综合指数的提高,产量、总干重以及总吸氮量也均呈增加趋势。增加干物质积累量是增加氮素吸收量的有效途径[30],而且在增加生物产量的基础上,进一步优化产量结构,是夺取水稻超高产的最佳路径[31]。本研究结果也发现,高产氮高效类型(类型Ⅰ)总干重的最小值均要比低产氮低效类型(类型Ⅲ)的最大值大,说明增加干物质积累量是促进产量和氮素吸收协同提高的关键所在。水稻地上部干物质积累量可以理解为群体茎蘖数与单茎重的乘积。本研究结果表明有效穗数在不同产量-氮效率类型之间差异不显著,与魏海燕等[9]研究结果一致,说明通过提高群体茎糵数来提高干物质生产量的途径可能行不通。因此,笔者推测增加水稻干物质积累量最有效的方式可能是通过提高单茎重,特别是通过提高单茎茎鞘重,这是因为较高的单茎茎鞘重,有利于形成较大的叶面积,同时促进大穗的形成,提高产量[32]。在较高的生物产量基础上,进一步增加库容,不仅能够增加氮素吸收量,同时还能促进氮素向穗部转移,提高氮素的物质和籽粒生产效率[33,34,35]。水稻库容量大小由单位面积有效穗数、每穗粒数和千粒重三个因素决定。从本研究结果看,产量构成因素中只有穗粒数在不同产量-氮效率类型之间的差异达到显著水平,且以高产氮高效类型(类型Ⅰ)最高,低产氮低效类型(类型Ⅲ)最低,与前人研究结果基本一致[36],说明选择穗粒数较多的品种,有利于实现产量与氮效率的协同提升。值得注意的是,增加生物产量能够增加吸氮量,这虽然有利于提高氮肥农学利用率、氮素生理利用率、氮肥吸收效率,但不利于氮素籽粒生产效率及氮素干物质生产效率的提升[37]。本研究也有类似的结果,氮肥农学利用率、氮素生理利用率、氮肥回收效率在3个不同产量-氮效率类型间均存在显著差异,而氮素籽粒生产效率和氮素干物质生产效率在类型Ⅱ与类型Ⅲ之间差异不显著,这可能与吸收的氮素更多地滞留在秸秆内,导致氮素利用效率降低有关[38]。氮素滞留在秸秆内可能是导致本研究中茎、叶吸氮量在筛选出来的3个产量-氮效率类型之间基本没有显著差异的主要原因之一。由此可见,促进茎、叶中的氮素更多地向穗转移,提高穗的吸氮量,可能是进一步提高水稻品种物质生产效率,增加产量和氮素利用效率的有效方式[27,28]。综上所述,结合本研究结果,通过以较大的单茎茎鞘重形成较高的地上部生物量,同时保证较多的穗粒数和较高的茎、叶氮素转运效率,能够实现产量与氮效率的协同提升。4 结论
运用熵权模糊隶属函数法,以氮肥回收效率、氮肥农学利用率、氮素生理利用率、氮素籽粒生产效率、氮素干物质生产效率为氮效率评价指标,建立了水稻氮效率综合评价方法,并基于产量-氮效率综合指数,分别从2017年90个和2018年105个供试品种(系)中筛选出了23和27个高产氮高效品种(系),其中南粳5718、南粳9108、宁粳7号、泗稻15号、扬粳239等19个品种(系)2年均被划分为高产氮高效类型。与低产氮低效类型品种(系)相比,高产氮高效类型品种(系)主要表现出生物量大、穗粒数多、穗吸氮量以及总吸氮量高等特征。上述研究结果可以为长江中下游地区常规中熟粳稻品种选择提供参考,同时也可为产量与氮效率协同机理的研究提供材料。参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 3]
[本文引用: 3]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 3]
[本文引用: 3]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 3]
[本文引用: 3]
,
[本文引用: 2]
[本文引用: 2]
,
[本文引用: 2]
[本文引用: 2]
,
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 2]
[本文引用: 2]
,
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 2]
[本文引用: 2]
,
[本文引用: 2]
[本文引用: 2]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]
,
[本文引用: 1]
[本文引用: 1]