基于软件的图像分析系统为研究匍匐茎和根状茎物种的形态提供了一种替代方法。该方案可用于测量匍匐茎和根状茎的长度与直径,适用于生物量较大的样品,并可应用于多种不同物种。
方法文章
基于软件的图像分析系统为研究匍匐茎和根状茎物种的形态提供了一种替代方法。该方案可用于测量匍匐茎和根状茎的长度与直径,适用于生物量较大的样品,并可应用于多种不同物种。
通常使用简单的直尺和卡尺测量匍匐茎或根状茎的长度和直径。该方法耗时且费力,因此通常仅用于有限数量的匍匐茎或根状茎。正因如此,这些性状在植物形态学表征中的应用受到限制。采用数字图像分析软件技术可克服因人为操作错误导致的测量误差,而此类误差往往随着样本数量和尺寸的增加而上升。该方案适用于各类作物,但尤其适用于植株较小且数量较多的牧草或禾本科植物。草坪样品包括地上生物量以及深度达到最大根状茎发育范围的表层土壤,具体深度取决于目标物种。在研究中,样品需从土壤中清洗出来,并在进行数字图像分析软件分析前人工清理匍匐茎/根状茎。随后将样品置于实验室烘箱中烘干以测定干重;因此,每个样品所得数据包括总长度、总干重和平均直径。在分析前,可通过排除扫描图像中可见的多余部分(如清洗过程中未完全去除的残留根系或叶片)对图像进行校正。事实上,这些碎片的直径通常远小于匍匐茎或根状茎,因此可通过设定最小直径阈值(低于该值的对象不纳入分析)轻松将其排除在分析之外。随后可根据样品面积计算单位面积内的匍匐茎或根状茎密度。该方法的优势在于能够快速高效地测量大量匍匐茎或根状茎样本的长度和平均直径。
植物形态学的研究广泛涉及植物科学的各个学科,包括生态学、农学、生物学和生理学。植物根系因其在胁迫耐受性、土壤稳定性、植物生长和生产力方面的重要性而被广泛研究。匍匐茎和根状茎也因其在植物繁殖策略、恢复能力以及碳水化合物储存中的作用而受到广泛关注。匍匐茎和根状茎是横向生长的变态茎,分别生长在地表之上(匍匐茎)或地表之下(根状茎)。匍匐茎和根状茎上还具有规律排列的节和节间,以及具有分生能力的节,能够产生新的根和芽1。已有大量研究针对不同植物的根系、匍匐茎和根状茎开展了多方面的探讨2,3,4,5,6,7,8。草坪草的根系、匍匐茎和根状茎因其对草坪质量9、冬季休眠后的春季返青能力10以及耐磨性和恢复能力的重要影响而被深入研究11。此外,这些器官也在其他作物中被研究,例如水稻12、大豆4和玉米13等草坪草类作物,以及在牧草系统中,其中侧生茎在控制土壤侵蚀方面发挥关键作用5。
根长密度(单位土壤体积内的根长)和平均直径通常使用扫描软件进行测量3,4,5,9,14,15,16,17,18。相反,匍匐茎或根状茎的长度和直径通常使用直尺和游标卡尺进行测量3,19,20,且需要耗费大量时间和人力21,22,23,24。因此,这类测量通常仅限于少量匍匐茎或根状茎11,20,25,且往往仅用于孤立植株的形态学表征。对成熟冠层中匍匐茎和根状茎性状的研究需要采集大量生物量,因此通常仅测定匍匐茎和根状茎的干重密度(单位面积上的干重)7,11,26,27。事实上,通过将样品在烘箱中烘干,比测量长度和直径更容易获得匍匐茎干重。然而,匍匐茎长度是一个重要的物种和品种特征,与干重的相关性较差。最近一项关于匍匐型多年生黑麦草(Lolium perenne)的研究表明,具有高匍匐茎长度密度的样品不一定具有高的匍匐茎重量密度6。
图像分析系统相比传统的手工方法31,32,33,能够更快28,29、更准确地分析根系,且更不易受人为误差影响30,21。此外,这些系统提供了高度的灵活性和易于使用的工具,包括光源、光学设置和分辨率,通常针对每幅扫描图像进行校准34。Pornaro et al.24 证明,WinRHIZO 系统是一种专为清洗后根系测量设计的图像分析系统,可通过克服人为错误导致的测量偏差,更全面地分析匍匐茎和根状茎性状,从而提供一种替代现有方法的新途径。对于匍匐茎和根状茎生长的形态学描述及定量信息获取,图像分析系统能够快速分析大量样本,即使生物量较大时也能实现,从而提高统计精度。因此,根系分析软件为研究不同植物物种的匍匐茎和根状茎的生长与形态提供了一种可靠、快速的替代方法24。
本研究在意大利东北部开展,旨在探讨四种狗牙根(Cynodon spp.)品种的匍匐茎和根状茎发育情况。该研究旨在增进对种子繁殖型("LaPaloma" 和 "Yukon")与营养繁殖型("Patriot" 和 "Tifway")狗牙根匍匐茎和根状茎发育过程的理解。试验于2013年5月建立,自2013年秋季至2015年夏季,每年分三次取样(分别为3月(返青前)、7月(生长旺盛期)和10月(冬季休眠前))。在本方法的描述与说明中,我们采用了第二个生长季夏季(2014年7月)采集的样品,因为此时样品生物量较大,适合进行快速分析。采用WinRHIZO软件——一种专为清洗后根系测量设计的数字图像分析工具——测定匍匐茎长度密度和平均直径。
1. 生物量样品的采集
2. 生物质样品的清洗
3. 样品的扫描与图像分析
4. 干重的测定
5. 数据校正及长度与重量密度的计算
于2013年秋季开展了一项田间试验,旨在比较四种狗牙根品种的匍匐茎和根状茎发育情况,包括两个种子繁殖类型(“LaPaloma”和“Yukon”)以及两个不育的营养型杂交种(“Patriot”和“Tifway”)。试验采用随机完全区组设计,设三个重复,共12个小区(每个小区2 m × 2 m)。
从每种草坪型品种以及野生狗牙根的试验小区内,以及小区附近生长的野生狗牙根植株中,随机采集14个匍匐茎和14个根状茎,共计70个匍匐茎和70个根状茎。所有匍匐茎和根状茎在进一步测量前均按照实验方案(步骤2)进行清洗。使用游标卡尺测量节间直径,使用直尺测量节间长度,并统计每个匍匐茎或根状茎的节数。同时记录使用直尺和游标卡尺清洗及测量匍匐茎和根状茎样本所需的时间。匍匐茎和根状茎的直径计算为所测所有节间直径的平均值。总匍匐茎长度和总根状茎长度计算为所有节间长度之和。此外,使用数字图像分析系统测量每个匍匐茎和根状茎的总扫描长度和扫描直径,具体方法见步骤3和步骤5。同时记录使用数字图像分析系统测量匍匐茎和根状茎性状所需的时间。随后用剪刀将每个匍匐茎和根状茎剪开,以分离节间与节,节间部分用于估算扫描得到的节间直径,方法见步骤3和步骤5。针对匍匐茎和根状茎(n = 70个匍匐茎,n = 70个根状茎),计算以下指标间的皮尔逊相关系数:实测长度与扫描长度、实测直径与扫描直径、节数与实测直径和扫描直径差值的绝对值、实测直径与扫描节间直径。使用直尺测量的长度用于校准数字图像分析系统估算的长度。
回归分析表明,匍匐茎扫描长度与实测长度之间具有高度相关性(图1a),斜率为1.03,截距为-4.22;根状茎扫描长度与实测长度之间也具有高度相关性(图1b),斜率为1.03,截距为4.22。人工清理14条匍匐茎和14条根状茎的平均耗时分别为21分钟24秒和11分钟12秒。使用直尺和游标卡尺测量长度和直径的平均耗时,匍匐茎为14分钟6秒,根状茎为13分钟35秒。使用WinRHIZO软件对样品进行扫描及分析的平均耗时,匍匐茎为11分钟,根状茎为12分钟4秒。
匍匐茎和根状茎的实测直径与扫描直径也呈显著相关。实测直径与扫描直径之间的关系接近1:1,表明数据拟合良好(图2a 和 2b)。然而,回归截距表明,数字图像分析系统对实测直径存在高估现象,尤其在较小直径值时更为明显,而对较大根状茎直径值则存在低估。这种高估可能是由于软件在扫描过程中包含了匍匐茎的节,从而影响了用于计算直径的总投影表面积(总投影表面积与总长度的比值),而使用游标卡尺测量时则排除了节的影响。节点数量与两种方法(实测与扫描)所得直径值差异之间的相关性仅在匍匐茎中达到显著水平(图3a);此外,节点数量的变化仅能解释该差异变异的一小部分(R2 = 14%)。扫描所得节间直径与实测直径之间存在显著相关性(匍匐茎和根状茎的斜率分别为1.01和0.98,截距接近零)(图4a 和 4b),表明只要去除节点,即可通过数字图像分析系统准确估算节间直径。因此,由多个匍匐茎或根状茎组成的样品的总长度和平均直径均可通过数字图像分析系统进行快速且准确的量化。
作为一项持续实验的一部分,从2013年秋季至2015年夏季,每个小区每季采集一个草坪样品(20 × 20 × 15 cm 深度),并按照本方案所述方法进行处理。2014年7月采集的样品中,单位面积匍匐茎和根状茎的长度(长度密度)及其单位面积重量(重量密度)如图5所示。营养繁殖品种(“Patriot”和“Tifway”)与种子繁殖品种(“La Paloma”和“Yukon”)在匍匐茎长度密度上存在差异。“Patriot”的根状茎长度密度最高,其次为“Tifway”和种子繁殖品种。所有品种的匍匐茎重量密度均不同,其中“Patriot”最高,随后依次为“Tifway”、“La Paloma”和“Yukon”。营养繁殖品种的根状茎重量密度也高于种子繁殖品种。品种“Patriot”在整个研究期间单位面积匍匐茎和根状茎长度(长度密度)及其单位面积重量(重量密度)的变化情况如图6所示。匍匐茎长度密度从2014年3月至2014年7月有所增加,此后从2014年7月至2015年7月保持稳定。2013年10月和2014年3月采集的样品中仅发现少量根状茎。根状茎长度密度在2014年7月增加并达到最高值,但在2014年10月再次下降。匍匐茎重量密度从2014年3月至7月略有上升;然而,从2014年7月至10月上升更为迅速,随后在2015年3月有所下降。根状茎重量密度的变化趋势与根状茎长度密度相似,在2014年7月达到最高值。
该软件在分析中包含扫描图像中的所有对象。图中展示了 WinRHIZO 软件生成的数字图像分析布局示例(图7),其中不同颜色的线条覆盖了不同直径的匍匐茎对象,以计算每个直径等级的总长度。我们可以观察到,该分析包含了根或叶的碎片。如步骤3.9所述,可以限制被分析的直径等级的宽度和数量。直方图显示了长度在所选直径等级中的分布情况(图7)。该直方图可用于评估应排除的最小直径等级。通过目视观察屏幕图像上部的该图可知,长度在平均直径等级周围呈正态分布,但前两个等级的数值明显高于符合正态分布的预期值。即使样品已仔细清洗,若包含这些较小的直径等级,数据分析仍可能影响结果,导致长度密度被高估,而平均直径被低估。我们的结果显示,在软件分析得出的总根状茎长度中,较小直径等级(直径 < 0.2 mm)所占比例为13–32%(表1)。此外,平均直径被低估了2–17%(表1)。

图1:百慕大草匍匐茎24(a)和根状茎(b)使用直尺测量的长度值与通过数字图像分析系统估算的长度值之间的回归分析。 虚线表示1:1的比例关系。图A经Pornaro 等24修改。 请点击此处查看该图的放大版本。

图2:利用游标卡尺测量的狗牙根匍匐茎直径值与数字图像分析系统估算值的回归分析24 (a)和根茎(b)。 虚线代表1:1的比例。图A已根据Pornaro修改 等24. 请点击此处以查看此图的放大版本。

图3:狗牙根匍匐茎24(a)和根状茎(b)的节点数量与通过数字图像分析系统估算值和游标卡尺实测值之间的直径差绝对值的回归分析。 虚线表示1:1比例关系。图A已根据Pornaro et al.24 进行修改。请点击此处查看此图的放大版本。

图4:百慕大草匍匐茎24(a)和根状茎(b)节间直径的游标卡尺实测值与数字图像分析系统估算值之间的回归分析。 虚线表示1:1比例关系。图A经Pornaro et al.24修改。 请点击此处查看该图的放大版本。

图 5:四种草坪狗牙根品种(Patriot、Tifway、La Paloma、Yukon)匍匐茎和根状茎长度与重量密度的田间试验示例结果。 匍匐茎长度密度(a)、根状茎长度密度(b)、匍匐茎重量密度(c)和根状茎重量密度(d)。垂直条形表示六次重复的标准误差。请点击此处查看此图的放大版本。

图6:来自田间试验的匍匐茎和根状茎长度密度与重量密度的示例结果,显示Patriot狗牙根品种的匍匐茎和根状茎发育情况。 匍匐茎长度密度(a)、根状茎长度密度(b)、匍匐茎重量密度(c)和根状茎重量密度(d)。垂直条形表示六次重复的标准误差。请点击此处查看此图的放大版本。

图 7:WinRHIZO 软件进行数字图像分析的示例布局。 前景中的扫描图像以及屏幕图像上方的柱状图显示了选定直径等级中的长度分布情况。彩色线条表示图像分析结果,每种颜色对应柱状图中报告的相应直径等级的颜色。请点击此处查看该图的放大版本。
| 品种 | 区组 | 根茎长度 (cm/dm2) | 平均直径 (mm) | ||||
| < 0.2 mm | 总计 | 比率a | < 0.2 mm | 总计 | 比率b | ||
| Patriot | 1 | 231 | 278 | 16.9 | 1.637846 | 1.5994 | 97.7 |
| Patriot | 2 | 304 | 349 | 12.8 | 1.620667 | 1.588371 | 98.0 |
| Patriot | 3 | 304 | 366 | 16.8 | 1.649918 | 1.621367 | 98.3 |
| Tifway | 1 | 184 | 231 | 20.6 | 2.149745 | 1.9951 | 92.8 |
| Tifway | 2 | 155 | 193 | 19.9 | 1.866253 | 1.76605 | 94.6 |
| Tifway | 3 | 119 | 150 | 20.9 | 1.877386 | 1.75865 | 93.7 |
| La Paloma | 1 | 17 | 23 | 24.4 | 2.139019 | 1.8904 | 88.4 |
| La Paloma | 2 | 26 | 38 | 31.6 | 2.101385 | 1.7455 | 83.1 |
| La Paloma | 3 | 34 | 47 | 27.5 | 2.033729 | 1.7354 | 85.3 |
| Yukon | 1 | 32 | 44 | 28.0 | 1.700155 | 1.4945 | 87.9 |
| Yukon | 2 | 17 | 25 | 33.2 | 1.68339 | 1.4284 | 84.9 |
| Yukon | 3 | 67 | 87 | 23.6 | 1.844721 | 1.6774 | 90.9 |
| a ≤0.2 mm 类别长度/总长度 | |||||||
| b 总直径/≤0.2 mm 类别直径 | |||||||
表1:包含与不包含较小直径等级的根茎长度密度及根茎平均直径。 包含与不包含直径小于0.2 mm等级时的长度密度及其比值(直径≤0.2 mm等级的长度/总长度);以及包含与不包含直径小于0.2 mm等级时的平均直径及其比值(包含直径<0.2 mm等级的平均直径 / 不包含直径<0.2 mm等级的平均直径)。
本文所述方案是为草坪草研究开发并评估的。然而,通过根据植物的形态特征、环境条件以及样品清洁精度进行必要的调整,该方案也可应用于其他多种匍匐茎或根状茎植物。
通过本方案估算的平均直径不能与使用游标卡尺测量的节间直径进行比较。数字图像分析在计算平均直径时包含了节和节间,该平均直径定义为总投影表面积与总长度的比值。Pornaro 等24指出,使用WinRHIZO系统测得的狗牙根匍匐茎平均直径比游标卡尺在节间处测得的实际直径值偏大。匍匐茎直径通常用于描述匍匐茎节间的直径,是植物形态描述中常用的一个参数18,25。因此,Pornaro 等24强调,通过WinRHIZO系统估算的平均直径与人工测量的节间直径反映的是两种不同的形态学特征。
完成本方案所需的时间仍是常规分析中的一个限制因素。最耗时的阶段是样品的清洗(步骤2.4)。根据我们的经验,清洗一份生物量较大的草皮样品即, 20 x 20 cm)的样本通常需要约三名人员工作2至4小时。如实验方案所述,无论使用数字分析系统,还是使用游标卡尺和直尺,清洁步骤都是必要的。当样本中匍匐茎或根状茎数量较少时,两种方法采集数据所需时间相近。然而,随着样本量增加,基于软件的方法不会显著增加所需时间,其唯一限制因素是扫描仪的扫描面积;相反,使用直尺和游标卡尺测量器官所需时间则随样本中匍匐茎或根状茎数量的增加而相应延长。
对成熟草坪草匍匐茎和根状茎性状的研究一直基于节间长度、直径以及干重的测量7,11,26,27。由于样品处理耗时较长,且随着样品数量增加,人工测量的准确性会下降,因此人工测量应仅限于少量匍匐茎或根状茎的测定11,20,25,可能仅适用于单株植物实验。与传统方法相比,图像分析系统的优势在于能够测量大量匍匐茎或根状茎样品的长度,并可同时计算长度密度和比重量(重量-长度比)。
本方案可用于测量具有较大生物量样品中的匍匐茎和根状茎长度,并计算其长度密度(对于此类样品,目前仅以匍匐茎或根状茎重量作为形态描述的参数)。在许多研究中,匍匐茎和/或根状茎的长度可能是一个重要参数,但现有技术无法对其进行准确估算。近期针对不同草坪草物种的研究6表明,匍匐茎的重量密度与长度密度并不总呈相关性,提示有必要测量多个参数以充分评估匍匐茎和根状茎系统。该方法特别适用于品种间或栽培管理措施的比较研究。
在实验方案中,有多个步骤对于准确估算匍匐茎和根状茎的长度及平均直径至关重要。由于植物在不同环境条件下的形态具有高度变异性,因此应仔细评估样本数量(样本量)以及采样地块的面积范围(采样尺度),使其尽可能代表整个种群,以降低数据的变异性。此外,在分析前需仔细清除匍匐茎上的根和叶片,此项工作较为精细,需特别注意以避免高估结果。最后,在处理图像之前,建议通过软件设置仔细选择直径分级的宽度以及最小直径,以排除所有非匍匐茎或根状茎的结构进入分析。由于直径大小随物种和环境条件(包括栽培措施)而变化,每个实验均需设定合适的最小直径阈值。
作者无任何利益冲突需要披露。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 实验室胶带 | 任意 | NA | 可使用标签标记样品 |
| 塑料袋 | 任意 | NA | 可使用任意塑料袋在样品清洗前保存样品 |
| 纸袋 | 任意 | NA | 可使用任意纸袋存放已清洗的样品,以避免霉菌形成 |
| 纸巾 | 任意 | NA | 样品用水清洗后、在用剪刀清理前,将其放置于纸巾上以吸收水分会更方便 |
| 剪刀 | 任意 | NA | 任何具有细尖头的剪刀 |
| 铝盒 | 任意 | NA | 任何足够容纳样品的铝盒 |
| 托盘 | 任意 | NA | 在清洗过程中使用塑料托盘盛放样品较为方便 |
| 孔径为 0.5–1.5 mm 的筛子 | 任意 | NA | 任意筛子 |
| 土壤钻取采样器 | 任意 | NA | 我们使用直径至少为 8 cm 的土壤钻取采样器进行土壤采集 |
| 方形采样框 | 任意 | NA | 采集大样品时,我们使用方形采样框(10 × 10 cm、或 15 × 15 cm、或 20 × 20 cm) |
| 铁铲 | 任意 | NA | 我们使用铁铲挖出由方形采样框划定的样品 |
| 精密电子天平 | 任意 | NA | 任意精密电子天平 |
| 实验室烘箱 | 任意 | NA | 任意实验室烘箱 |
| 冷冻柜 | 任意 | NA | 任意冷冻柜 |
| WinRHIZO 软件 | Regent Instruments Inc., Quebec | NA | 不包括 "basic" 版本 |
| WinRHIZO 扫描仪 | Regent Instruments Inc., Quebec | NA | WinRHIZO 系统包含一台与软件校准匹配的扫描仪 |
| WinRHIZO 扫描仪配件 | Regent Instruments Inc., Quebec | NA | WinRHIZO 系统包含配件,如塑料托盘和定位器,用于配合扫描仪使用 |
申请许可以重复使用本 JoVE 文章的文本或图表
申请许可