陆稻早期幼苗陆生适应性鉴定及综合评价

张德宇，CHALI Biru Alemu，PHYU Phyu Thin，路承凯，陆炳悦，陆定科，王效宁，万金鹏，徐 鹏

(1 中国科学院 西双版纳热带植物园 热带植物资源可持续利用重点实验室,云南 昆明 650223；2 中国科学院大学,北京 101408；3 广南县农业技术推广服务中心,云南 广南 663300；4 海南省农业科学院粮食作物研究所/海南省农作物遗传育种重点实验室,海南 海口 570216)

世界一半以上的人口以稻米为主食[1]。陆稻(Upland rice)和水稻(Lowland rice)是亚洲栽培稻OryzasativaL.的2种生态型[2]。陆稻，也称旱稻或者旱谷，通常是指在旱地、坡地及干旱生态环境下正常生长的栽培稻类型[3]，具有节约水资源、轻简宜机、劳动强度低和少甲烷排放等优点。陆稻种植不需灌溉水层，水分供给主要依赖于自然降雨或少量灌溉[4]。陆稻在适应陆生旱地环境过程中，形成了对陆生环境土壤养分、水分、季节和日长变化等的适应性[5-6]。与水稻相比，典型的陆稻具备低水分条件下种子萌发、萌发后出土、苗期与杂草竞争、趋避不利环境(如快速卷叶、快速灌浆成熟)、耐旱，以及适应于水分供给的光周期反应等能力，这些对陆生环境的适应表现被称为陆生适应性(Aerobic adaptation)[6]。

陆稻在选择驯化过程中，形成大量地方种质[7]。最近发现，陆稻可能先于水稻更早从野生稻中驯化而来[8]。基于水田生态系统育成的水稻品种普遍不适应通气有氧的陆生旱地环境[9-10]。传统的陆稻生产主要分布在我国云南、贵州、海南等山区及东南亚国家，具有深播出土、抗(耐)旱等特性[11-12]。虽然传统陆稻的产量低，但其蕴含的优异基因是稻作遗传育种的宝库。在全球气候变化背景下，水资源短缺已成为传统水稻生产最为主要的限制因素[13]。将传统陆稻与现代水稻改良品种相结合进行遗传改良，是应对水资源短缺和干旱等挑战的有效途径[14-15]。因此，有待进一步加强特异陆稻种质的挖掘和创新利用，但目前缺乏对陆稻陆生适应性关键节点性状简便、高效的系统鉴定方法，极大地限制了陆稻在遗传育种中的利用价值。

目前关于陆稻早期幼苗陆生适应性的综合性评价，特别是与陆稻生产实践紧密结合的低水分胁迫条件下的萌发和深土出土能力，仍缺乏系统深入地研究。现阶段对于水、陆稻种子低水分胁迫条件下萌发的研究，主要以聚乙二醇 (Polyethylene glycol 6000,PEG6000)作为诱导条件[16]。李自超等[17]对比了陆稻和水稻在不同质量浓度PEG6000处理条件下幼苗的生长情况，陆稻生长势比水稻更好。在200 g/L PEG6000溶液中培养时，陆稻仍有一定的生长量，而水稻均不能发芽。胡德勇等[18]调查了不同水分处理对陆稻种子萌发的影响，发现陆稻品种‘IAPAR’种子萌发的临界含水率可低至田间持水量的30%。徐建欣等[19]发现，在150 g/L PEG6000处理下，陆稻的相对发芽率、相对幼根长度等性状指标优于水稻，可视为陆稻对低水分环境萌发的适应特征。基于此方法，科学家已从陆稻中鉴定出许多抗旱基因[11,20-21]。

深播被认为是旱直播栽培中抗旱保墒的重要途径[22]。胚芽鞘具有保护胚芽中幼叶和生长锥的作用[23]，在其下方有一个短茎，称为中胚轴。胚芽鞘和中胚轴在幼苗出土过程中起到关键作用。不同作物的种子在萌发出土过程中延伸器官可能不同。高粱深播中胚轴伸长使其幼苗顶出土壤表面，大麦和小麦的中胚轴几乎没有伸长，而它们的胚芽鞘显著伸长[24-25]。有研究表明，胚芽鞘和中胚轴伸长被认为是水稻幼苗出土的主要动力源[26-28]。当前研究主要用出苗率来表示出土能力，出苗率越高其出土能力越强[29-30]，而大多研究均以水稻为主[31-32]，缺乏对特有陆稻种质资源出土能力相关指标的综合评价[33]。通过开展陆稻早期幼苗陆生适应性的综合评价，筛选具有较强低水分胁迫萌发及深土出苗能力的种质资源，可为后续陆生适应性遗传机理解析、节水栽培和耕作研究及陆稻育种奠定基础。

世界各地的陆稻，包括低纬度低海拔(含赤道附近)的陆稻地方品种及改良品种以粳稻为主，水稻与陆稻的遗传分化主要发生在粳稻亚种内[15]。云南省是中国稻种遗传和生态多样性中心之一，其南部雨热同季、多山的环境孕育了丰富的陆稻资源，为筛选具有优良陆生适应性的种质提供了大量的资源[34-35]。本研究利用150 g/L PEG6000模拟低水分胁迫和8 cm深土播种进行典型水、陆稻早期幼苗陆生适应性的差异指标筛选，同时对246份地方陆稻种质资源进行综合评价，以期明确评价标准，筛选出具有优良陆生适应性的代表性品种，为进一步创新利用云南陆稻种质资源提供理论依据。

1 材料与方法

1.1 试验材料

试验材料选用典型水稻和改良陆稻品种。24份典型的水稻品种主要来源于长江中下游及东北地区；24份典型陆稻品种主要来源于国内外改良品种，246份地方陆稻品种均来源于云南省，由中国科学院西双版纳热带植物园作物保护与育种基地扩繁保育提供试验。

1.2 试验方法

1.2.1 陆稻低水分胁迫萌发能力鉴定 试验于2022年9—12月在中国科学院西双版纳热带植物园温室内进行。参试材料选取籽粒饱满、大小一致的种子60粒，使用质量分数为1%的次氯酸钠消毒液消毒15 min，清洗干净后将种子置于垫有滤纸的9 cm培养皿中。设置对照组(蒸馏水)和低水分胁迫组(150 g/L PEG6000溶液) 2个处理，每处理3个重复。处理组加入5 mL的150 g/L PEG6000溶液，对照组加入等体积的蒸馏水，培养皿加盖置于28 ℃恒温、光周期为14 h光照/10 h黑暗的温室中培养。从种子置床之日起开始观察，以胚根或胚芽达种子长的1/2为发芽标准，每天定时统计发芽种子数[36]。处理7 d后，每个处理取10株幼苗，测量芽长、根长、根数等指标。

1.2.2 陆稻深土出苗能力鉴定 每个参试材料选取45粒饱满的种子，使用孔径为8 mm的筛子除去土壤大颗粒，填充于11 cm 深、底部钻有9个孔的黑色塑料培养盒中，每个塑料盒先平铺3 cm土壤，播种后将盒子填满(播种深度8 cm)，以模拟田间土壤播种条件。然后将培养盒放入3 cm水深的容器中，反渗吸水确保土壤充分浸润，放置24 h后取出置于苗床上培养。每盒播种3个材料，设置3个重复。7 d后调查幼苗生长情况，取样测量胚芽鞘长、中胚轴长、根数和根长等指标。

1.3 数据处理与分析

试验数据采用Excel 2016录入并绘制图表，利用SPSS软件进行描述性统计、方差和主成分分析。采用模糊数学隶属函数法进行综合鉴定，并进行聚类分析[37]。具体计算公式如下：

隶属函数值计算公式：

权重(Wj)计算公式：

综合低水分胁迫萌发/出土能力(D)计算公式：

式中，Xj为某项指标的测定值，Xmax为全部供试材料某项指标的最大值，Xmin为全部供试材料某项指标的最小值；Pj表示第j个综合指标贡献率；uj表示第j个综合指标的隶属函数，Wj表示第j个权重。

2 结果与分析

2.1 典型水稻和改良陆稻品种早期幼苗的陆生适应性评价

在早期幼苗陆生适应性的相关性状中，典型水稻、改良陆稻品种间观测的几乎所有指标均表现出极显著差异，陆稻的相对发芽率与相对发芽势大于0.700，而水稻则小于0.150。陆稻的相对根长、相对根数分别为水稻的5倍和3倍；此外，水稻的各相对指标的变异系数均大于200%，而陆稻除相对根长外均低于70%，表现出更稳定、更强的适应性。对这些材料的深土出苗观察分析发现，除根数外，其余性状指标在水、陆稻间均存在显著的差异，特别是根长、中胚轴长显示出极显著差异，根长、根数、胚芽鞘与芽长较水稻的变异程度更低，陆稻表现出更强的适应深播的能力(表1)。

表1 典型水、陆稻陆生适应性差异性分析Table 1 Difference analysis of aerobic adaptability between typical lowland rice and upland rice

在150 g/L PEG6000处理下，有3份典型水稻品种的芽长、根数和根长均为0，而5份改良陆稻品种的芽长和根数为典型水稻品种的10倍，根长为典型水稻品种的30～50倍，均显著强于典型的水稻品种(图1)。

图1 不同水、陆稻品种在150 g/L PEG6000条件下培养7 d后的生长情况Fig.1 Growth of different lowland rice and upland rice varieties after cultured for 7 days under 150 g/L PEG6000 conditions

在深土萌发下，陆稻品种的各项指标表现更好，同时发现部分水稻品种的根数、根长上也有不错的表现，可进一步用于发掘相关优良基因(图2)。以上结果表明，低水分胁迫萌发、出土能力等相关指标除根数外，均可用于陆稻早期幼苗陆生适应性的筛选。

图2 不同水、陆稻品种在深土8 cm深土播种条件下培养7 d后的生长情况Fig.2 Growth of different lowland rice and upland rice varieties after cultured for 7 days under 8 cm deep-seeding conditions

2.2 早期幼苗陆生适应性性状的相关性分析

在150 g/L PEG6000胁迫条件下，对294份供试材料的相对芽长、相对根数、相对发芽率等5个指标的相对值进行相关性分析，结果表明各性状均呈现显著相关，其中相对芽长与相对根数之间的相关系数达0.611，相对芽长与相对发芽率相关系数达0.465。这些结果说明陆稻苗期的陆生适应性可以在多个性状指标上得以体现，需要综合考虑多个指标，才能进行更准确地鉴定(表2)。

表2 150 g/L PEG6000 胁迫下生长性状的相关性分析1)Table 2 Correlation analysis of growth traits under 150 g/L PEG6000 conditions

在8 cm深土播种条件下，除中胚轴与根长、芽长之间无显著相关外，其余各性状指标之间均显著相关。芽长与胚芽鞘长相关系数达0.691，这与在幼苗出土过程中，胚芽鞘保护幼芽穿过土壤的生物学过程一致(表3)。此外，胚芽鞘长与中胚轴长间存在极显著相关性，与前人关于中胚轴和胚芽鞘的伸长是种子出苗的主要动力源的观点相吻合[26]。根数在陆稻中与各指标均存在显著相关性，而在典型水陆稻中不存在显著差异(表3)。这些结果表明数量性状的复杂性，说明陆稻幼苗的出土能力受到多个性状的共同调控。

表3 8 cm深土播种条件下生长性状的相关性分析1)Table 3 Correlation analysis of growth traits under 8 cm deep-seeding conditions

2.3 早期幼苗陆生适应性的主成分分析

主成分分析能够提取出数据中的主要特征和结构，保留数据中最重要的信息，可方便快捷地对供试样品的相似性进行科学评价。对早期幼苗陆生适应性的低水分胁迫萌发和出土能力指标分别进行主成分分析，以累计贡献率大于80%为原则选择主成分。在150 g/L PEG6000处理条件下，因子分析提取了3个主成分，其贡献率分别为50.038%、19.683%和 12.547%，累计贡献率达82.268%，表明可以解释原始数据82.268%的总变异，涵盖了所有指标的大部分信息。在8 cm深土条件下，由于典型水陆稻间根数性状不存在显著差异，因而只对根长、中胚轴、胚芽鞘、芽长进行主成分分析，该过程将原来的4个单项指标组合为2个新的互相独立的综合指标，2个主成分贡献率分别为55.082%和26.351%，累计贡献率为81.433%，可以代替大多数指标。

主成分与所有相关性状的因子负荷量反映了彼此间的相关性。在150 g/L PEG6000条件下，构成载荷分析表明，在第1主成分中所有性状指标均起着重要作用；相对根数的特征向量(0.814)最高，表明了根系对于水分吸收的重要性；在第2主成分中，相对发芽势特征向量(0.851)数值最大，说明种子的自身活力同样对低水分萌发起着重要的作用；在第3主成分中，相对根长对应特征向量(0.637)数值最大，表明发达的根系对于适应低水分的重要性。对深土出苗评价发现，在第1主成分当中，除中胚轴外，其余指标的特征向量均较高，它们起着主要作用；在第2主成分中，中胚轴特征向量较高，对深土出苗发挥着主要作用，表明在深土条件下这些指标均需要考虑作为幼苗出土能力的重要评价指标(表4)。

表4 在150 g/L PEG6000和8 cm深土播种条件下主成分性状特征值、贡献率和累积贡献率1)Table 4 The characteristic value,contribution rate and cumulative contribution rate of principal component traits under 150 g/L PEG6000 and 8 cm deep-seeding conditions

2.4 早期幼苗陆生适应性的聚类分析

通过主成分分析得到因子得分值，作为鉴定陆生适应性的综合指标。利用公式(4)计算各综合指标的隶属函数值，再用公式(5)通过综合指标的贡献率计算权重。在150 g/L PEG6000胁迫下，3个主成分的权重依次是0.608、0.239和0.153。在8 cm深土条件下，2个主成分的权重依次是0.608和0.392。最后采用公式(5)计算得到综合评价值D。对D均值采用欧式距离平方法进行系统聚类分析(表5)。

表5 294份参试材料的陆生适应性等级及分类Table 5 Aerobic adaptation level and classification of 294 tested germplasms

结果表明，在欧式距离等于5时，可将294份供试材料的陆生适应性划分为极强型(Ⅰ类)、强型(Ⅱ类)、中间型(Ⅲ类)、弱型(Ⅳ类)、极弱型(Ⅴ类)5类。第I类群包含‘N90-253大长谷’‘黑花谷’‘夺哥谷2’‘大香糯’‘大麻谷’‘红壳谷’和‘黄壳谷’等7份种质，为来自于云南的地方陆稻品种，说明地方种质具有优异的陆生适应性特性，对地方资源的挖掘评价具有重大意义。第II类群包含96份种质，除地方品种外，还涵盖有4份改良的陆稻品种。第Ⅲ类群包含157份种质，除了3份水稻外，绝大多数为地方及改良陆稻品种。第Ⅳ类群包含29份种质，大部分为典型水稻品种，但也有5份改良的陆稻品种，这表明在陆稻改良品种时可能注重于单个性状，而忽视了对其陆生适应性的方面。第Ⅴ类群包含5份种质，全为典型水稻品种，说明典型水稻品种的陆生适应性确有很大的提升潜力(表5)。这些结果表明，陆生适应性的评价是一个全面综合的评价过程，一些改良的陆稻由于育种过程可能仅注重单方面抗性的提高，忽略了出土能力等性状的改良，导致其在陆生适应性综合评价中依旧表现不佳。参试云南地方陆稻品种的陆生适应性强弱参差不齐，说明该省陆稻地方种质资源类型多样、遗传丰富，筛选出的陆生适应性强的地方陆稻品种可深入挖掘并加以利用。

3 讨论与结论

水、陆稻在陆生适应性相关性状方面存在明显的遗传差异，陆稻适应性关键遗传差异是水、陆稻遗传差异的基础[38-39]；陆稻的陆生适应性需在生长发育过程中与土壤水分供应特征保持一致，并具备对干旱逆境的抵抗、忍耐与恢复能力，以及较强的出土能力[40]。目前，缺乏有效的陆生适应性系统鉴定评价方法，极大地限制了优良陆稻种质资源在育种中的利用。在幼苗的抗旱性鉴定研究中，主要利用PEG6000溶液模拟低水分胁迫环境，将萌发相关的生长指标作为抗旱鉴定的指标，已成为幼苗抗旱鉴定研究的重要手段[41]。杨瑰丽等[42]研究发现，干旱胁迫严重影响了水稻的萌发，明显抑制根的生长，与芽相比，根中的干物质积累易受干旱胁迫的抑制。李其勇等[43]研究模拟干旱对水稻育种材料‘川香29B’近等基因导入系种子萌发的影响，发现从轻度到重度的水分胁迫对种子发芽的抑制作用呈现增强的趋势。丁国华等[44]发现，干旱胁迫条件下杂草稻的发芽率、茎长、根长、根数、茎干质量、根干质量和根冠比等指标均显著降低，与水稻对同样胁迫条件的表现没有显著差异。

水、陆稻对干旱胁迫响应存在着巨大的遗传多样性，特别是地方陆稻品种往往具有很强的干旱胁迫耐受性[45-46]。陆稻地方品种虽然存在较多不足，但其适应于干旱环境等潜在的遗传变异可以用于稻作新品种耐旱性的遗传改良。传统的陆稻主要分布于湿润旱地，生长季节与雨热同期，耐旱性可能不是陆稻驯化的主要推动力[47]。陆稻耐直播性使其适应低水分萌发、耐深埋、扎根快，有利于其适应陆生有氧的旱地环境并快速完成苗期形态建成，是陆稻陆生适应性的重要体现。发达的根系通过维持水分和矿物质，在赋予植物耐旱性方面发挥着重要作用，陆稻品种深、粗的根系使其能从深层土壤中获取水分来抵御水分亏缺[48]。Gao等[12]发现水稻根系发育的负调控因子RobustRootSystem1(RRS1)在抗旱性中起着重要作用。我们的研究发现相对根数、根长与低水分胁迫萌发能力显著相关，表明发达的根系对于适应旱生环境具有重要作用。

胚芽鞘和中胚轴的伸长被认为是幼苗出土的主要动力源[26-28]。中胚轴长是影响水稻种子出苗的主要因素，中胚轴伸长对深播种子的出苗起关键作用，长胚轴种质出土动力来自于胚芽鞘与中胚轴的共同作用；而短中胚轴出苗动力主要来源于胚芽鞘的伸长[28]。我们的研究发现，芽长与胚芽鞘长、胚芽鞘长与中胚轴长间存在显著相关性，说明胚芽鞘在保护嫩芽穿过土壤中起着重要作用[23,27]。已有研究通过数量性状位点(QTL)定位和全基因组关联分析的方法,发掘了50余个水稻中胚轴伸长相关位点，分布在水稻全基因组中[31]。Zhao等[49]发现控制中胚轴伸长的2个主效基因OsML1和OsML2，其优良单倍型组合可使水稻中胚轴长度增加。通过GWAS对208个水稻种质的中胚轴进行分析，发现16个与中胚轴伸长显着相关的基因座[31]。利用长中胚轴品种‘Changai’和短中胚轴品种‘IR 145’构建F2遗传分离群体，发现LOC_Os03g52450、LOC_Os03g56060和LOC_Os03g58290为水稻中胚轴伸长相关的候选基因，它们参与植物激素和细胞分裂调控[32]。通过研究播种深度与中胚轴和胚芽鞘长度相关的QTL，发现qMel-1和qMel-3增加品种‘Kasalath’中胚轴长，而qMel-6增加‘Nipponbare’的中胚轴长；另外，qCol-3和qCol-5增加‘Nipponbare’的胚芽鞘长[50]。本研究发掘7个优异陆稻种质资源将可用于胚芽鞘和中胚轴等关键性状调控机理的研究。

在我国的云南、海南、四川和贵州等山区，坡地种植陆稻的生长主要依赖于自然降雨，较强的陆生适应性是保障前期出苗率和后期齐苗、壮苗的关键环节。本研究紧密结合陆稻生产中种子低水分胁迫萌发和深播幼苗出土问题，通过对早期幼苗低水分胁迫萌发和出土能力鉴定，采用相关性状指标相对值的主成分、隶属函数等分析方法对陆生适应性进行综合评价。发现150 g/L PEG6000胁迫下的相对发芽率、相对发芽势和相对根长，以及8 cm深土下的胚芽鞘长和中胚轴长，可作为陆稻早期幼苗陆生适应性鉴定的主要指标。通过聚类分析，将294份种质资源分为陆生适应性极强型(7个)、强型(96个)、中间型(157个)、弱型(29个)和极弱型(5个)。本研究建立了陆稻早期幼苗陆生适应性的筛选评价体系，为后续陆稻种质资源挖掘、创新利用及优良陆生适应性新品种的培育提供了种质资源和理论依据。