荔枝霜疫病由荔枝霜疫霉Peronophythora litchii引起，是造成荔枝果实褐变腐烂、产量降低和产量不稳定的主要病害，在我国荔枝栽培区普遍发生^[1-2]；该病为害荔枝新梢、嫩叶及花穗，在气候条件适宜发病的年份，可引起荔枝产业经济损失高达80%以上，严重限制了荔枝产业的发展^[3]。

关于荔枝霜疫霉对温湿度适应性的研究多为温度或者湿度单独对荔枝霜疫霉孢子囊萌发和侵染的影响，Chen等^[4]研究台湾地区的荔枝霜疫霉在V8汁培养基中，菌丝生长的温度范围为8~32 ℃，最适生长温度为24~28 ℃；Liu等^[5]研究华南地区的荔枝霜疫霉在PDA培养基中，菌丝生长温度为12~32 ℃，最适生长温度为28~30 ℃；蔡学清等^[6]研究认为在24~28 ℃条件下，果实最易被荔枝霜疫霉侵染且发病最重；黄河等^[7]认为荔枝霜疫霉孢子囊遇水立即脱落。荔枝霜疫病的发生和流行与日照时数、降雨天数有关，日照时数和降雨天数直接影响果园的温度和寄主露时，连续降雨数天即可造成该病害流行，且病害来势凶猛，一旦发生便迅速扩展蔓延^[8]。因此，监测种植区环境条件，结合物候期对病害发生进行有效的预测预报，在其发病前期和初期采取有效的防治措施，对遏制病害的流行为害至关重要。Montone 等^[9]通过监测草莓叶片露时得到相应的预测模型应用于草莓病害的预警系统中；Roubal等^[10]依据降雨量、温度及相对湿度等环境因素预测Fusicladium oleagineum对橄榄的侵染及潜伏期；Uysal等^[11]探索了温度、露时及叶龄对菠菜炭疽病菌侵染及扩展的影响。而有关荔枝霜疫病发生的流行条件、预测预报方法等报道较少^[12]。

农作物病虫害预测预报是植物保护工作的基础，是实现防灾减害的重要保证。本研究在人工控温、控湿条件下，研究温度、露时(即寄主湿润持续时间)对荔枝霜疫霉孢子囊萌发和侵染的影响，以期为荔枝霜疫病监测预警提供依据，提高预测预报准确性。

1.   材料和方法

1.1   材料

叶片和果实接种试验选取的荔枝品种均为‘怀枝’。

荔枝霜疫霉菌株(编号SHS3)由华南农业大学植物病理学系真菌研究室分离、鉴定和保存。菌株培养用PDA培养基。

1.2   方法

1.2.1   温度和露时对孢子囊萌发的影响

叶片和果实上孢子囊萌发的观察参考蒋继志等^[13]的方法并稍作改动：选取长势一致的嫩叶和成熟果实置于保鲜盒中，将玻璃纸剪成均等大小(2 cm×2 cm)，于无菌水中浸湿后附着于叶片或果实表面；然后吸取10 μL的孢子囊悬浮液(1×10⁴ mL⁻¹)滴于玻璃纸上，每个叶片或果实滴5滴，将叶片或果实分别放置在15、20、25和30 ℃条件下培养，待达到相应的露时2、4、8、12和24 h后，置玻璃纸于载玻片上，在显微镜下观察孢子囊的萌发情况。每个处理设置3次重复，每次重复至少计数100个孢子囊。

1.2.2   温度和露时对叶片和果实感病的影响

荔枝叶片的接种方法参考Neufeld 等^[14]和Arauz 等^[15]的方法并稍作改动。将孢子囊悬浮液喷雾接种于嫩叶上，至有液滴滴落。荔枝果实的接种方法：取新鲜的荔枝成熟果实用清水清洗后晾干，将果实在孢子囊悬浮液中浸泡1 min后取出。

将接种过的叶片分别于15、20、25和30 ℃条件下培养，接种过的果实分别于20、25和30 ℃条件下培养，待达到相应的露时2、4、8、12和24 h后，用强风吹干叶片和果实，重新置于各自的培养温度下，24 h后将所有的处理移至白天24 ℃、晚上20 ℃、12 h光暗交替、75%~90% 相对湿度下继续培养7 d。每天观察叶片或果实的发病情况并记录发病等级，计算接种3 d的病情指数。每个处理设置3次重复，每次重复至少接种50片叶片或30颗果实。

叶片发病分级标准：0级，无病斑；1级，病斑面积占叶片面积的1%~25%；2级，病斑面积占叶片面积的26%~50%；3级，病斑面积占叶片面积的51%~75%；4级，病斑面积占叶片面积的75%以上。

果实发病分级标准参考Qu等^[16]的方法并稍作修改。0级：果皮未出现病斑；1级：病斑面积占整个果皮面积的5%以下；3级：病斑面积占整个果皮面积的6%~15%；5级：病斑面积占整个果皮面积的16%~25%；7级：病斑面积占整个果皮面积的26%~50%；9级：病斑面积占整个果皮面积的50%以上。

病害发生严重度计算公式如下：

1.3   数据模型拟合

模型拟合采用SAS(9.4版本)软件进行分析。在Neufeld等^[14]分析方法的基础上略作修改，首先用初步方差分析法分析孢子囊萌发试验的3组重复、孢子囊侵染叶片试验的3组重复、侵染果实的2组重复是否有明显差异，若无明显差异，后续则取各组重复试验的平均值进行分析。其次对孢子囊在叶片、果实上萌发及侵染试验分两步进行分析，第一步使用SAS中的PROC GLM过程分析比较温度(θ)、露时(t)以及它们的交互效应(θ×t)对孢子囊萌发和寄主发病严重度的影响；第二步则以方程f(t, θ) ={1−exp[−(B×t)²]}/cosh[(θ−F)×G/2] 拟合模型。该模型是韦布尔函数的修正式，其中斜体小写字母为自变量，大写斜体字母为参数。B代表变量y响应自变量露时(t)的固有增长率，F和G分别代表变量y响应最佳温度和温度(θ)的固有变化率。使用SAS中的PROC NLIN过程，将孢子囊平均萌发率和寄主平均发病严重度分别拟合至非线性模型。使用SAS PROC REG过程分析孢子囊萌发率和寄主发病严重度观察值和预测值的相关性，评估模型拟合的可靠性。根据参数估计的显著性、渐近标准误差的大小评估模型的拟合优度，使用MATLAB MESH制作孢子囊在叶片、果实上萌发和侵染的三维曲面图。

1.4   病害风险预测模型

通过拟合模型f (t, θ)产生t和θ的函数，预测发病严重程度。根据最终模型参数估算温度和露时的范围，使用 MATLAB CONTOUR 制作等高线图，预测荔枝嫩叶和成熟果实霜疫病发生的严重度。

2.   结果与分析

2.1   孢子囊萌发及寄主发病程度响应温度和露时的方差分析

方差分析结果(表1)表明孢子囊萌发、孢子囊侵染嫩叶和孢子囊侵染成熟果实响应温度和露时的重复试验之间的差异不显著(P>0.05)，因此可使用试验数据平均值进行后续分析。温度、露时极显著影响孢子囊的萌发及叶片、果实接种3 d后的病害严重度(P<0.000 1)；温度和露时的交互作用也可以显著影响孢子囊的萌发及寄主发病的严重度(P<0.01)。

表  1  温度和露时对荔枝霜疫霉孢子囊萌发及荔枝叶片、果实发病程度影响的方差分析

Table  1.  Analysis of variance for the effects of temperature and wetness duration on sporangia germination of Peronophythora litchii and disease severities of litchi leaf and fruit

来源 Source	叶片 Leaf								果实 Fruit
	孢子囊萌发率 Sporangia germination rate				病害严重度 Disease severity				孢子囊萌发率 Sporangia germination rate				病害严重度 Disease severity
	df	SS	P		df	SS	P		df	SS	P		df	SS	P
重复 Replication	2	0	0.961 6		2	0.08	0.398 2		2	0	0.991 5		1	0	0.991 0
温度(θ) Temperature	3	0.89	0.000 1		3	0.61	0.000 1		3	0.38	0.000 1		2	0.06	0.000 1
露时(t) Wetness duration	4	1.41	0.000 1		4	0.46	0.000 1		4	6.62	0.000 1		4	0.71	0.000 1
θ×t	12	0.23	0.000 1		12	0.35	0.008 0		12	0.20	0.000 1		8	0.09	0.007 1

下载: 导出CSV 

| 显示表格

2.2   不同温度和露时对孢子囊萌发的影响情况

荔枝霜疫霉孢子囊在温度为15~30 ℃时，露时为2 h时即开始萌发，且萌发率随着露时的延长而逐渐增大；在温度为25 ℃时，所有露时处理条件下的孢子囊萌发率达到最大值。孢子囊在叶片上萌发率大于0.2的温度为25 ℃，露时≥2 h或温度20和30 ℃，露时≥4 h(图1A)；在果实上孢子囊萌发率大于0.6的温度为20和30 ℃，露时≥12 h或温度25 ℃，露时≥8 h(图1C)。

图  1  不同温度和露时处理下荔枝霜疫霉孢子囊在荔枝叶片和果实上的萌发率

Figure  1.  Sporangia germination rates of Peronophythora litchii on litchi leaf and fruit under the conditions of different temperature and wetness duration

下载: 全尺寸图片 幻灯片

2.3   不同温度和露时对叶片和果实发病程度的影响

露时相同的条件下，当温度为25 ℃时，孢子囊为害荔枝嫩叶和成熟果实最为严重；所有的温度处理下，叶片和果实的发病程度随着露时延长而加重。相同的处理条件下，孢子囊为害果实比叶片更为严重和迅速(图2)。荔枝叶片接种孢子囊后，温度为25 ℃，露时≥2 h时，病害严重度可达0.2；或者温度为20或30 ℃时，露时≥4 h；而温度为15 ℃时，则需露时≥12 h(图2A)。果实接种孢子囊后，温度为25 ℃、露时≥2 h或者温度为20和30 ℃、露时≥4 h，果实的病情指数即可超过0.6(图2C)。

图  2  荔枝叶片和果实接种荔枝霜疫霉孢子囊3 d后的病情指数

Figure  2.  Disease severity of litchi leaf and fruit at three days after inoculation with Peronophythora litchii sporangia

下载: 全尺寸图片 幻灯片

2.4   模型的参数估计和检验

选用韦布尔模型的修正式，在SAS软件中拟合得到了准确度较好的孢子囊萌发和发病严重度模型f(t, θ) ={1−exp[−(B×t)²]}/cosh[(θ−F)×G/2]。模型很好地描述了温度和露时对荔枝霜疫霉孢子囊萌发和侵染的影响。孢子囊萌发率和病害严重度随露时变化呈S形曲线，随温度变化呈单峰形曲线(图1B、1D和图2B、2D)。模型拟合的渐进参数估计值均有效且差异显著(P<0.05)，近似标准误差均较小(表2)。

表  2  模型参数估计¹⁾

Table  2.  Parameter estimates for the model

参数 Parameter	荔枝叶片 Litchi leaf										荔枝果实 Litchi fruit
	孢子囊萌发率 Sporangia germination rate					病害严重度 Disease severity					孢子囊萌发率 Sporangia germination rate					病害严重度 Disease severity
	e	ASE	ACIL	ACIU		e	ASE	ACIL	ACIU		e	ASE	ACIL	ACIU		e	ASE	ACIL	ACIU
B	0.13	0.01	0.10	0.16		0.19	0.05	0.09	0.30		0.12	0.01	0.10	0.15		0.47	0.03	0.39	0.54
F	24.91	0.61	23.63	26.19		27.77	0.29	27.16	28.38		26.46	11.56	2.08	50.84		24.56	0.96	22.47	26.64
G	0.39	0.05	0.30	0.49		0.93	0.12	0.67	1.18		0.08	0.10	−0.14	0.30		0.17	0.03	0.10	0.25
1) B代表变量y响应自变量露时(t)的固有增长率，F和G分别代表变量y响应最佳温度和温度(θ)的固有变化率，e为参数估计值，ASE为参数近似标准误差，ACIL、ACIU为参数在95%置信区间的近似下限和上限 　1) B represents the inherent growth rate of variable y responding to the independent variable(t), F and G represent the inherent changing rates of variable y responding to the best temperature and temperature(θ)respectively, e represents the estimate value of parameter, ASE is asymptotic standard error, ACIL and ACIU are the asymptotic lower and upper limits of estimates in 95% confidence interval respectively

下载: 导出CSV 

| 显示表格

从预测值和观测值相关性分析来看，孢子囊对叶片和果实萌发率的预测值和观测值差异不显著(P> 0.05)，即可认为预测值与观测值无差异，模型拟合效果较好，其决定系数(r²)分别为0.931 2和0.948 4；孢子囊侵染叶片和果实病害严重度的预测值和观测值差异不显著(P>0.05)，即可认为预测值与观测值无差异，模型拟合效果较好，决定系数分别为0.814 7和0.991 6(表3)。

表  3  模型检验与拟合公式

Table  3.  Verification and fitting formula of model

指标 Index	预测值与观测值相关性 Relationship of observed value and predicted value		模型拟合公式 Model fitting formula
	r2	P
孢子囊在荔枝叶片上的萌发率 Sporangia germination rate on litchi leaf	0.931 2	0.601 5	f(t，θ)={1−exp[−(0.126 6×t)2]}/cosh[(θ− 24.908 4)×0.394 2/2]
孢子囊在荔枝果实上的萌发率 Sporangia germination rate on litchi fruit	0.948 4	0.140 0	f(t，θ)={1−exp[−(0.123 1×t)2]}/cosh[(θ− 26.457 4)×0.084 4/2]
荔枝叶片发病严重度 Disease severity of litchi leaf	0.814 7	0.052 8	f(t，θ)={1−exp[−(0.194 3×t)2]}/cosh[(θ− 27.769 6)×0.927 7/2]
荔枝果实发病严重度 Disease severity of litchi fruit	0.991 6	0.919 1	f(t，θ)={1−exp[−(0.469 3×t)2]}/cosh[(θ− 24.556)×0.170 9/2]

下载: 导出CSV 

| 显示表格

2.5   病害风险预测

根据荔枝嫩叶和成熟果实接种孢子囊3 d后的病害严重度的模型，制得相应的风险预测等高线图(图3)。依据荔枝霜疫霉侵染叶片风险预测图可以推算出病害严重度从0.1~0.9变化所对应的温度和露时。当温度为22~30 ℃、露时约为3~24 h时，叶片的病害严重度达到0.2以上(图3A)。依据侵染果实风险预测图可以推算出病害严重度从0.5~0.9变化所对应的温度和露时。由图3B可知，荔枝霜疫霉为害果实更易成灾，当温度为15~30 ℃、露时约2~3 h时，果实发病的严重度高达0.6。因此，在荔枝果实接近成熟时，必须严格监控温度和露时，及时发布霜疫病发生流行的预警信息。

图  3  荔枝霜疫霉侵染荔枝叶片和果实的风险预测图

Figure  3.  Risk prediction charts for Peronophythora litchi infecting litchi leaf and fruit

下载: 全尺寸图片 幻灯片

3.   讨论与结论

温度和湿度在病原孢子产生和病害的发生过程中起着至关重要的作用^[17]。据报道，病原真菌孢子的近距离传播在一定程度上依赖于局部的温度和湿度^[18]。近年来，多用定量分析模型研究温度和露时对病原菌孢子萌发侵染的影响。在卵菌病害的研究中，Arauz等^[15]通过定量分析模型研究了温度和露时对黄瓜霜霉病菌Pseudoperonospora cubensis孢子萌发和侵染的影响；Rakotonindraina等^[19]则结合分析模型研究了不同气象条件、土质等外界因素影响下马铃薯晚疫病菌Phytophthora infestans对马铃薯的为害。在果园病害的研究中，Li等^[20]定量研究了温度和相对湿度对梨黑星病菌Venturia nashicola分生孢子萌发和侵染的影响。本研究在人工控温、控湿条件下，定量分析了4组温度(15~30 ℃)、5组露时(2~24 h)对荔枝霜疫霉孢子囊萌发和侵染的影响，利用SAS软件构建了以温度和露时为自变量的函数。并在此基础上，制得荔枝霜疫霉侵染嫩叶和成熟果实相对应的风险预测图。在卵菌侵染叶片的试验中，Neufeld等^[14]研究认为病害严重度不大于0.1时，为轻度发病，病害严重度为0.1~0.2时属于中度发病，不小于0.2则为重度发病。参考此标准，本研究认为霜疫霉侵染荔枝嫩叶时，当温度为22~30 ℃、露时为3~24 h时，叶片的病害严重度达到0.2以上，认定为嫩叶发病严重，应及时防治。孢子囊为害成熟果实则更加迅速和严重，当温度为15~30 ℃、露时约2~3 h时，果实的病害严重度即高达0.6，足可为害成灾。

本研究获得的荔枝霜疫霉侵染荔枝嫩叶和成熟果实的风险预测图可为实际生产中的荔枝霜疫病预测预报提供参考，但仍存在不足，例如在构建模型后，使用SAS PROC REG过程分析孢子囊萌发率和寄主发病严重度观察值和预测值的相关性，认为预测值与观测值无差异，模型拟合效果较好。模型的准确性，还要结合生产实际中病害发生流行期的温度和露时进行进一步验证。