Research progress in secretory cavity morphogenesis in Citrus
-
摘要:
分泌囊是柑橘属Citrus植物普遍存在的结构，是药用成分合成、积累的主要场所。目前研究显示柑橘属植物分泌囊以裂溶生方式发生，其分泌细胞的发育过程是细胞程序化死亡的过程。虽然如此，但是分泌囊发生发育过程的研究已经有百年的历程。本文重点介绍了柑橘属植物分泌囊的起源、发生发育方式和细胞程序化死亡过程，以期为分泌囊形态建成的未来研究提供参考。
Abstract:The secretory cavity is a ubiquitous structure in Citrus plants, which is the main place for the synthesis and accumulation of medicinal components. Current studies have shown that the secretory cavity of citrus plants occurs in schizolysogenesis, and the development of secretory cells is a process of programmed cell death. However, the study of the development of secretory cavity has experienced a hundred years. This review focuses on the origin, development and programmed cell death of secretory cavity in Citrus, in order to provide a reference for future studies on secretory cavity morphogenesis.
-
Keywords:
- Citrus /
- Secretory cavity /
- Schizolysigenous /
- PCD /
- Development
-
随着经济全球化的加深，我国进口货物和进口国家数量不断增长，同时，许多有害杂草伴随着国际流通进入我国[1]。研究表明，进口额每增长1个百分点，植物疫情在海关截获量约增长2个百分点[2]。近年来有多起南方三棘果Emex australis在海关截获的报告。例如，2012—2015年广东黄埔口岸在主要进境粮食中截获的有害生物中，南方三棘果截获高达90批次排第7位。此外在广东的肇庆和汕头口岸以及江苏口岸也都有截获报道[3-4]。
南方三棘果是蓼科亦模属的冬型一年生杂草，喜生长在湿润和半湿润热带、亚热带和温带气候区，被我国列入进口植物检疫对象。目前在我国内地尚未发现，仅在台湾有分布记录。南方三棘果单株植株可以产生1000粒以上的果实，种子在土壤中休眠多年，仍有繁殖能力，且具有繁殖速度快、传播迅速、危害性大和防治难度高的特点，一旦大面积入侵我国，难以根治，将会给农业生产造成巨大破坏[5-7]。
随着信息化技术的发展，目前用于预测物种潜在分布区的模型有CLIMEX、GARP、DIVA-GIS、BIOCLIM和MaxEnt等。相比于其他模型，MaxEnt生态位模型在分布数据较少的情况下，预测精确度和准确度高[8-9]。其通过将研究区域的环境因子转换为栅格数据，形成多个栅格单元，把物种已知分布区域的栅格单位作为样点，然后以分布点区域的环境变量数据对所有栅格单元进行数理统计，得到当前和未来一段时间内研究区域的每一个栅格单元分布的概率以及环境变量数值与分布概率之间的关系[10-12]。
由于各物种的研究和关注程度不同，物种在世界范围内详实的分布记录存在差异，因此物种的分布数据获得量存在差异[12]。南方三棘果属于数据获得量较少的物种，因此选择MaxEnt作为预测南方三棘果在我国潜在适生区的模型。全球气温呈现升高趋势，会影响植物的潜在分布区域变化[13-14]。气候系统模式能够很好地反映气候变化的趋势，对未来一段时间的气候变化进行准确预测[15-17]。因此本文采用国家气候中心4种气候情景下的未来气候数据对南方三棘果在我国的适生区进行分析预测。
本研究将MaxEnt模型与ArcGIS软件相结合，得到南方三棘果潜在分布区的可视化地图，以期能够为有关部门在潜在分布区的相关地区开展入侵风险评估、检验检疫和监测预警等提供理论参考，对防控南方三棘果传入、定殖和扩散，保障农业生产和生物多样性健康以及生态安全和经济贸易安全等具有重要意义[18]。
1. 材料与方法
1.1 技术路线
技术路线包括查找和处理数据、模型优化与构建以及制作可视化地图3部分，按照“物种分布点查找、处理→查找影响分布的环境因子以及SPSS分析→MaxEnt模型参数优化→构建MaxEnt模型→ArcGIS制作可视化图”的研究方案进行，数据分析贯穿于整个研究过程。
1.2 相关软件和程序
MaxEnt软件3.4.1版、ArcGIS软件10.4.1版、SPSS Statistics软件20.0版、R语言3.6.3版、ENMeval程序包以及Microsoft Excel 2016。
1.3 环境数据及处理
从世界气候数据网站( https://www.worldclim.org/)下载2.1版、分辨率为10 min(340 km2)的1970—2000年和国家气候中心气候系统模式(BCC-CSM2-MR)4个气候情景模式(SSP1-2.6、SSP 2-4.5、SSP 3-7.0和SSP5-8.5)2021—2040年世界范围内19个生物气候变量数据(Bio1~Bio19)以及高程海拔数据(Elev)。采用ArcGIS软件从高程海拔数据中提取出坡度(Slo)和坡向(Asp)[19-20]。一共得到温度、降水和海拔特征的环境变量数据22个(表1)。
表 1 环境变量因子Table 1. Environmental variable factors变量 Variable 描述 Description 变量 Variable 描述 Description Bio1 年平均温度/℃ Annual mean temperature Bio12 年均降水量/mm Annual mean precipitation Bio2 昼夜温差月均值/℃ Monthly mean diurnal temperature range Bio13 最湿月降水量/mm Precipitation of the wettest month Bio3 等温性/% (Bio2/Bio7) ×100 Isothermality Bio14 最干月降水量/mm Precipitation of the driest month Bio4 温度季节变化标准差/℃ Standard deviation of seasonal variation of temperature Bio15 降水量变异系数 Standard deviation of the precipitation Bio5 最暖月最高温/℃ Maximum temperature of the warmest month Bio16 最湿季度降水量/mm Precipitation of the wettest quarter Bio6 最冷月最低温/℃ Minimum temperature of the coldest month Bio17 最干季度降水量/mm Precipitation of the driest quarter Bio7 年均温度变化范围/℃ Temperature annual range Bio18 最暖季度降水量/mm Precipitation of the warmest quarter Bio8 最湿季度平均温度/℃ Mean temperature of the wettest quarter Bio19 最冷季度降水量/mm Precipitation of the coldest quarter Bio9 最干季度平均温度/℃ Mean temperature of the driest quarter Elev 高程海拔 Elevation Bio10 最暖季度平均温度/℃ Mean temperature of the warmest quarter Slo 坡度/° Slope Bio11 最冷季度平均温度/℃ Mean temperature of the coldest quarter Asp 坡向 Aspect 参考文献[21-22]方法，用ArcGIS将所有环境变量坐标系及基准面统一定义，并转出为ASCⅡ格式，以待MaxEnt模型预测分布；使用R语言将TIFF格式的环境变量因子统一转换为ASC格式，以待使用R语言ENMeval 程序包对MaxEnt模型进行优化。为了避免变量间的高度共线性导致的过拟合和减少运算量，对环境数据进行相关性分析。将物种分布点与环境变量输入MaxEnt模型运算10次，获得各个环境变量的贡献率，将环境数据和物种分布点数据在ArcGIS处理为“*．txt”格式文件，将文件利用Excel转化为“*．csv”格式文件，最后将分析数据导入SPSS进行双变量相关性分析，作双尾检验。对相关系数(R)>0.8的环境变量，选择贡献率较高的作为主要变量进行下一步分析，舍去贡献率较低项Bio5、Bio6、Bio7、Bio10、Bio12、Bio16和Bio17。舍去贡献率为0的环境变量Bio8。由于南方三棘果为冬型一年生草本植物，在最湿月和最暖季度是种子形态[23-24]，因此舍去Bio13和Bio18。最终得到12个环境变量因子(表1)，分别为Bio1、Bio2、Bio3、Bio4、Bio9、Bio11、Bio14、Bio15、Bio19、Elev、Slo和Asp。
1.4 南方三棘果分布数据及处理
在全球生物多样性信息交换平台( https://www. gbif.org/)等平台数据库检索物种分布相关记录，通过发表的文献、专著、文件、报告等查询相关分布信息及记录[25-27]。一共得到1473个分布数据点，去除重复、无效、缺失、错误分布数据后，得到1205个分布点数据。利用下载好的气候数据用ENMeval程序包对南方三棘果分布点进行过拟合分析，最终得到南方三棘果自然分布点647个。
1.5 地图数据
从标准地图服务系统( http://bzdt.ch.mnr.gov.cn/)下载GS(2022)4312号标准地图。
1.6 模型优化
使用R语言ENMeval程序包将南方三棘果分布点分为2部分，75%随机样本数据用于训练集，25%用于测试集，对MaxEnt软件中调控倍频(Regularization multiplier，RM)和特征组合(Feature combination，FC)参数进行优化，利用数据包设定RM从0.1到4，每次增加0.1，共40个调控倍频，同时设定29个特征组合，即H、L 、LH、LP、LPH、LPT、LPTH、LQ、LQH、LQP、LQPH、LQPT、LQPTH、LQT、LQTH、LT、P、PH、PT、PTH、Q、QH、QP、QPH、QPT、QT、QTH、T、TH，其中L为线性(Linear)，Q为二次型(Quadratic)，H为片段化(Hinge)，P为乘积型(Product)，T为阈值型(Threshold)，ENMeval 语言包将上述 1160 种参数组合进行检验，最终采用最小信息准则(The minimum information criterion AICc value)模型评估受试者工作特征曲线(Receiver operating curve，ROC)下的面积(Area under curve，AUC)差异检验模型的拟合度与复杂度，当AICc值最低(AICc=0)时模型参数组合最佳[28]，用于MaxEnt模型建模。
1.7 构建 MaxEnt 预测模型
参考前人所述方法[29-33]，将处理后的当前环境变量和物种分布点数据输入MaxEnt模型，模型参数设置随机选择25%的样本作为测试集，75%的样本数据为训练集，同时勾选制作响应曲线功能以获得分布概率与气候因子的关系，勾选刀切法检验环境变量的重要程度，RM与FC取优化后的参数，其他参数为默认设置。以AUC评价模型的准确性。AUC越大，表明模型的可靠性越高，0.9以上表示模型结果优秀可靠，为降低模型的随机性和模拟结果的不确定性，模型重复运行10次，以运行的均值作为模型运算结果。
1.8 绘制物种分布图
将MaxEnt模型生成的ASC格式物种分布数据导入ArcGIS软件，设置坐标系为地理坐标系GCS_WGS_1984，通过掩膜提取获得南方三棘果在我国的分布图，对预测结果利用自然间断点分级法进行重分类，得到南方三棘果的非适生区、低适生区、中适生区和高适生区。
2. 结果与分析
2.1 MaxEnt生态位模型优化结果
基于12个环境变量因子和南方三棘果的647个分布点数据，MaxEnt软件初始默认RM=1，FC为LQPH，AICc=142.939251，经ENMeval语言包优化后RM=1，FC为LQPT，AICc=0，表明其为最优参数组合。在优化参数下进行模型预测，采用ROC分析法进行模型准确度检验。经10次重复训练得到的AUC > 0.9( 图1)，模拟结果优秀可靠。
2.2 影响分布环境变量因子贡献率及适宜范围
使用12个环境因子建立模型，模型参数为ENMeval语言包优化结果，RM=1，FC为LQPT，采用刀切法评估各环境因子对南方三棘果地理分布的重要性。贡献率超过15%的包括最冷季度平均气温(Bio11，27.7%)、等温性(Bio3，24.9%)、温度季节变化标准差(Bio4，18.2%)，这3个因子的贡献率合计达70.8%。贡献率低于15%的包括年平均温度(Bio1，14.0%)、最干月降水量(Bio14，6.5%)、降水量变异系数(Bio15，6.5%)、最冷季度降水量(Bio19，2.7%)、高程海拔(Elev，1.4%)、坡度(Slo，0.6%)、最干季度平均温度(Bio9，0.3%)、昼夜温差月均值(Bio2，0.2%)、坡向(Asp，0.2%)。
刀切法检验表明，仅使用单一环境变量时，正则化训练增益和测试增益都表明最重要的环境因子为最冷季度平均温度(Bio11)(图2)。
![]()
图 2 环境变量重要性的刀切法检验结果Asp: 坡向；Slo: 坡度；Elev: 高程海拔；Bio9: 最干季度平均温度；Bio4: 温度季节变化标准差；Bio3: 等温性 ( Bio2/ Bio7) ×100；Bio2: 昼夜温差月均值；Bio19: 最冷季度降水量；Bio15: 降水量变异系数；Bio14: 最干月降水量；Bio11: 最冷季度平均温度；Bio1: 年平均温度Figure 2. Results of Jackknife for the importance of the environmental variablesAsp: Aspect; Slo: Slope; Elev: Elevation; Bio9: Mean temperature of the driest quarter; Bio4: Standard deviation of seasonal variation of temperature; Bio3: Isothermality; Bio2: Mean diurnal range; Bio19: Precipitation of the coldest quarter; Bio15: Precipitation seasonality; Bio14: Precipitation of the driest month; Bio11: Mean temperature of the coldest quarter; Bio1: Annual mean temperature一般认为存在概率≥0.5对应的环境因子范围是物种适宜分布的阈值。根据响应曲线确定影响南方三棘果分布环境因子的阈值：最冷季度平均温度9.35~12.76 ℃、等温性45.2%~53.72%以及温度季节变化标准差337.06~616.40 ℃(图3)。
2.3 南方三棘果在当前气候下的潜在适生分布区
将南方三棘果分布点和环境变量数据导入MaxEnt模型进行模拟，得到南方三棘果在当前(2022年)气候下的潜在分布区(图4)。总适生面积约69.13万km2，高适生区、中适生区、低适生区分别为11.64万、13.37万、44.02万km2。南方三棘果的适生范围主要集中在我国云南、广西、广东、福建以及西藏局部，其中，高适生区和中适生区主要位于云南、海南以及西藏局部，低适生区主要位于广西、广东、台湾、福建和四川。
![]()
图 4 南方三棘果在中国的适生区(2022年)Figure 4. Current suitable regions of Emex australis in China (2022)2.4 南方三棘果在未来气候模型下的潜在适生分布区
如表4所示，SSP1-2.6气候情景下，2040年南方三棘果总适生区面积约63.54万km2，较2022年分布面积减少约5.48万km2，其中，低适生区面积减少5.93万km2，中适生区面积增加 6.57万km2，高适生区面积减少6.12万km2；SSP2-4.5气候情景下，2040年南方三棘果总适生区面积约65.10万km2，较2022年分布面积减少约3.92万km2，其中，低适生区面积减少2.31万km2，中适生区面积增加2.31万km2，高适生区面积减少3.92万km2；SSP3-7.0气候情景下， 2040年南方三棘果总适生区面积约63.34万km2，较2022年分布面积减少约5.68万km2，其中，低适生区面积减少9.74万km2，中适生区面积增加6.51万km2，高适生区面积减少2.45万km2；SSP5-8.5情景下，2040年南方三棘果总适生区面积约62.90万km2，较2022年分布面积减少约6.12万km2，其中，低适生区面积减少5.42万km2，中适生区面积增加3.39万km2，高适生区面积减少4.09万km2。
表 4 不同时期气候情景下南方三棘果在中国潜在适生区的面积Table 4. Areas of the potential suitable region for Emex australis in China under different climate scenarios in different periods×104 km2 年份 Year 气候情景 Climate scenario 低适生区 Low suitable region 中适生区 Middle suitable region 高适生区 High suitable region 总适生区 Total suitable region 2022 44.02 13.36 11.63 69.02 2040 SSP1-2.6 38.09 19.93 5.51 63.54 2040 SSP2-4.5 41.71 15.67 7.71 65.10 2040 SSP3-7.0 34.27 19.88 9.18 63.34 2040 SSP5-8.5 38.59 16.76 7.54 62.90 如图5所示，对不同气候情景下南方三棘果适生区的空间格局进行对比分析，结果表明，2040年，未来4种情景模式下的适生区都呈现轻微缩小趋势，广东和广西的部分低适生区变为非适生区，云南的高适生区变为中适生区，海南的高、中、低适生区都呈减少趋势。从单个情景模式下看也有一些地区变得更加适应南方三棘果分布，SSP1-2.6气候情景模式下西藏出现局部高适生区，四川的局部低适生区转变为中适生区；SSP2-4.5气候情景模式下福建的一些非适生区变为低适生区；SSP3-7.0气候情景模式下四川部分低适生区转变为中适生区，局部地区变为高适生区，西藏出现局部高适生区；SSP5-8.5气候情景模式下，四川和西藏的一些非适生区变为低适生区，其中，西藏出现局部高适生区。这些变化区域为南方三棘果分布的敏感区。
![]()
图 5 4种气候情景模式下南方三棘果在中国的适生区Figure 5. The suitable growth reqion of Emex australis in China under four climate scenarios3. 结论与讨论
本模型的AUC测试值 > 0.9，模型结果优秀可靠。最冷季度平均气温(Bio11)、等温性(Bio3)和温度季节变化标准差(Bio4)是影响南方三棘果地理分布的主要环境因子。正则化训练增益和测试增益以及贡献率都表明最冷季度平均气温(Bio11)对南方三棘果分布影响最大，贡献率为27.7%。环境因子响应曲线表明，当最冷季度平均气温为9.35~12.76 ℃时，南方三棘果发生概率≥0.5。南方三棘果的适宜分布区总面积约为69.02万km 2，适生区域主要集中于云南、广东、广西和福建地区，其中高、中适生区主要集中于云南。
通过对2022—2040年国家气候中心气候系统模式(BCC-CSM2-MR)的4种未来气候情景下的南方三棘果潜在分布区研究发现，总体上适生区呈现轻微减少趋势。四川、西藏、广东和广西属于南方三棘果分布的敏感区，易随气候变化而导致适生区等级以及适生面积发生变化。广东、广西、海南和云南的一些适生区在未来一段时间会出现转变为非适生区或降低适生等级的情况。福建部分地区在SSP2-4.5情景模式下会更适宜南方三棘果生存，西藏和四川部分地区在SSP1-2.6、SSP3-7.0和SSP5-8.5情景下会更适宜南方三棘果生存。
其中，西藏的适生区尤为特殊，主要集中在藏南地区的错那县和墨脱县以及朗县最南部。这一部分地区由于地形原因导致气候复杂，高原气候和亚热带、热带季风气候同时存在，加上地处边疆，适生区西边与不丹接壤，南边与印度接壤，导致藏南外来物种入侵严重，值得相关部门进一步关注。中国科学院植物研究所的植物物种信息系统( http://www.iplant.cn/)显示三叶鬼针草Bidens pilosa、胜红蓟Ageratum conyzoides等原产于热带地区的外来杂草在墨脱县都有分布记录[34]。国家气候中心气候系统模式数据显示未来这一地区的温度季节变化标准差(Bio4)为390~550 ℃，等温性(Bio3)为39%~45%，最冷季度平均温度(Bio11)为0.9~19 ℃，这种气候适宜南方三棘果分布，因此，对这里需要加大监测力度，阻止南方三棘果进入藏南地区，并防止由藏南地区入侵我国西南和华南地区。
南方三棘果是具有严重危害的检疫杂草，对农业生产会造成巨大破坏[35-36]。云南是我国葡萄十大产地之一，南方三棘果在我国的高、中适生区主要集中于云南，南方三棘果一旦定殖于葡萄园，可能会有部分种子在机械化收割时与葡萄一起收割，影响葡萄的品相和食品安全性[37]。藏南是西藏的主要农区，藏南地区生产的小麦是西藏的主要粮食作物之一，研究表明南方三棘果茎的水提取物和植株残留物会抑制小麦种子萌发和早期幼苗的生长，降低小麦产量[38-39]。油菜是广西、广东和福建的主要油料作物，在农田里南方三棘果与油料作物存在种间竞争关系[40]。
外来入侵杂草传入与人员和货物的流动密切相关，作为国与国之间管控人员和货物流通的重要节点，海关和动植物检验检疫局实施的植物检疫是防止外来杂草传入的最重要方式。在适生区域内以及邻近适生区域的海关和动植物检验检疫局应当加大对南方三棘果的检疫力度，防止南方三棘果的传入[41]。尤其值得注意的是，广东、广西和福建的海关每年都有大量的货物进口，加上这些地区存在大面积的南方三棘果适生区，海关和动植物检验检疫局的植物疫情数据也表明在广东和福建的海关有多起南方三棘果的截获记录，因此，这些地区是南方三棘果传入的高风险地区。
有害杂草入侵的定殖和扩散阶段在适生区范围内进行，其中，外部气候环境的变化是影响其扩散速度的重要因素，试验表明最冷季度平均温度是影响南方三棘果分布的最重要环境因子，在适生区气候适宜其生存的重要年份要加大南方三棘果的监测，有效防止其定殖、扩散，以免影响农作物的质量和产量[42]。生物入侵杂草的分布与气候环境、地形因素、土壤类型以及植被类型等多种因素相关，在使用MaxEnt生态位模型预测入侵杂草潜在分布区时，目前尚未有统一的影响因子选择。朱煜建等[43]在预测疣果匙荠Bunias orientalis适生区时选择了气候环境、地形因素以及土壤类型的环境变量因子。郭燕青等[44]在预测假臭草Praxelis clematidea时仅选择了气候因素有关因子。本研究选用与环境和地形有关的变量因子，缺少土壤类型、物种类型等与南方三棘果分布有关的其他变量，这可能对适生区的预测带来一定的偏差。伴随着我国与世界各国日益频繁的贸易往来，生物入侵难以避免。近些年假高粱Sorghum halepense和豚草Ambrosia artemisiifolia等入侵生物对我国农业生产和生态安全造成了巨大破坏。采用MaxEnt生态位模型将信息学技术与生物的生态学特性结合，对各种入侵杂草的潜在分布区进行预测，对我国的口岸检疫和入侵杂草的防控具有指导意义。
-
![]()
图 1 柑橘属分泌囊发育过程模式图
a：原始细胞早期；b：原始细胞中期；c：原始细胞晚期；d：腔隙形成期；e：空腔扩大期；f：成熟期
Figure 1. Secretory cavity development model in Citrus
a: The early stage of initial cell; b: The middle stage of initial cell; c: The late stage of initial cell; d: The lumen formation stage; e: The lumen expansion stage; f: The mature stage
-
[1] VAN TIEQHEM P. Deuxieme memoire sur les canaux secreteurs[J]. Annales des Sciences Naturelles Botanique, 1885: 1-96.
[2] SCHULZE H, BLATTANAT B. Rutaceae[M]. Heidelberg: Diss, 1902: 50.
[3] LIANG S J, WU H, LUN X, et al. Secretory cavity development and its relationship with the accumulation of essential oil in fruits of Citrus medica L. var. sarcodactylis (Noot. ) Swingle[J]. Journal of Integrative Plant Biology, 2006, 48(5): 573-583. doi: 10.1111/j.1744-7909.2006.00230.x
[4] WU Z, LI H, YANG Y, et al. Variation in the components and antioxidant activity of Citrus medica L. var. sarcodactylis essential oils at different stages of maturity[J]. Industrial Crops and Products, 2013, 46: 311-316. doi: 10.1016/j.indcrop.2013.02.015
[5] MARTINET J. Organes de secretion des vegetaux[J]. Annales Des Sciences Naturelles Botanique, 1872, 14(5): 91-232.
[6] HABERLANDT G. Physioloqical plant anotomy[M]. London: Macmillan Press, 1914: 777.
[7] FOHN M. Zurentstehunq und weiterbildlung der exkreträume von Citrus medica L. und Eucalyptus globulus Lab[J]. Ö sterreichische Botanische Zeitschrift, 1935, 84(3): 198-209.
[8] KISSER J G. Die Ausscheidung von atherischen glen und harzen[J]. Handbuch der Pflanzenphysiologie, 1958, 10: 91-131.
[9] AMELUNXEN F, ARBEITER H. Untersuchungen an den Spritzdrusen von Dictamnus albus L.[J]. Zeitschrift Pflanzenphysiologie, 1967, 58: 49-69.
[10] HEINRICH G. Licht- und elektronenmikroskopische untersuchungen zur genese der exkrete in den lysigenen exkretrumen von Citrus medica[J]. Flora Oder Allgemeine Botanische Zeitung. Abt. A, Physiologie und Biochemie, 1966, 156(5): 451-460. doi: 10.1016/S0367-1836(17)30213-6
[11] HEINRICH G. Elektronenmikroskopische beobachtungen zur entstehungsweise der exkretbehlter von Ruta graveolens, Citrus limon und Poncirus trifoliata[J]. Ö sterreichische Botanische Zeitschrift, 1969, 117(4): 397-403.
[12] HEINRICH G. Elektronenmikroskopische beobachtungen an den drüsenzellen von Poncirus trifoliata; zugleich ein beitrag zur wirkung ätherischeröle auf pflanzenzellen und eine method zur unterscheidung flüchtiger von nichtflüchtigen lipophilen komponenten[J]. Protoplasma, 1970, 69(1): 15-36. doi: 10.1007/BF01276649
[13] ESAU K. Plant anatomy[M]. 2nd ed. New York: Wiley, 1965.
[14] ESAU K. Anatomy of seed plants[M]. 2nd ed. New York: Wiley, 1977.
[15] BOSABALIDIS A, TSEKOS I. Ultrastructural studies on the secretory cavities of Citrus deliciosa ten. I. Early stages of the gland cells differentiation[J]. Protoplasma, 1982, 112: 55-62. doi: 10.1007/BF01280215
[16] BOSABALIDIS A, TSEKOS I. Ultrastructural studies on the secretory cavities of Citrus deliciosa ten. Ⅱ. Early stages of the gland cells differentiation[J]. Protoplasma, 1982, 112: 63-70. doi: 10.1007/BF01280216
[17] TSCHIRCH A, STOCK E. Die harze und die harzbehälter[M]. Beilin: Gebrüder Borntraeger, 1933.
[18] THOMSON W W, PLATT-ALOIA K A, ENDRESS A G. Ultrastructure of oil gland development in the leaf of Citrus sinensis L[J]. Botanical Gazette, 1976, 137: 330-340. doi: 10.1086/336880
[19] FAHN A. Secretory tissues in plants[M]. London: Academic Press, 1979: 176-209.
[20] TURNER G W, BERRY A M, GIFFORD E M. Schizogenous secretory cavities of Citrus limon (L. ) Burm. F. and a reevaluation of the lysigenous gland concept[J]. International Journal of Plant Sciences, 1998, 159: 75-88. doi: 10.1086/297523
[21] TURNER G W. A brief history of the lysigenous gland hypothesis[J]. The Botanical Review, 1999, 65(1): 76-88. doi: 10.1007/BF02856558
[22] KNIGHT T G, KLIEBER A, SEDGLEY M. The relationship between oil gland and fruit development in Washington navel orange (Citrus sinensis L. Osbeck)[J]. Annals of Botany, 2001, 88(6): 1039-1047. doi: 10.1006/anbo.2001.1546
[23] KNIGHT T G, KLIEBER A, SEDGLEY M. Structural basis of the rind disorder oleocellosis in Washington navel orange (Citrus sinensis L. Osbeck)[J]. Annals of Botany, 2002, 90(6): 765-773. doi: 10.1093/aob/mcf258
[24] BROCHERIOU J, BELIN-DEPOUX M. Contribution a I etude ontogenique des poches des feuilles de quelques Myrtaceae[J]. Phytumurphology, 1974, 24: 321-338.
[25] SIECK W. Die schizolysigenen secretbehalter[J]. Jahrb Wiss Botanik, 1895, 27: 197-242.
[26] BUVAT R. Ontogeny, cell differentiation and structure of vascular plants[M]. Berlin: Springer-Verlag, 1989.
[27] SPRECHER E. BeitrRge zur frage der biogenese sekundirer pflanzenstoffe der weinraute (Ruta graveolens L.)[J]. Planta, 1956, 47(4): 323-358. doi: 10.1007/BF01911861
[28] BENNICI A, TANI C. Anatomical and ultrastructural study of the secretory cavity development of Citrus sinensis and Citrus limon: Evaluation of schizolysigenous ontogeny[J]. Flora-Morphology, Distribution, Functional Ecology of Plants, 2004, 199(6): 464-475. doi: 10.1078/0367-2530-00174
[29] CHEN Y, WU H. Programmed cell death involved in the schizolysigenous formation of the secretory cavity in Citrus sinensis L. (Osbeck)[J]. Chinese Science Bulletin, 2010, 55(20): 2160-2168. doi: 10.1007/s11434-010-3275-4
[30] LIU P, LIANG S, YAO N, et al. Programmed cell death of secretory cavity cells in fruits of Citrus grandis cv. Tomentosa is associated with activation of caspase 3-like protease[J]. Trees, 2012, 26(6): 1821-1835. doi: 10.1007/s00468-012-0752-1
[31] BAI M, LIANG M, HUAI B, et al. Ca2+-dependent nuclease is involved in DNA degradation during the formation of the secretory cavity by programmed cell death in fruit of Citrus grandis ‘Tomentosa’[J]. Journal of Experimental Botany, 2020, 71(16): 4812-4827. doi: 10.1093/jxb/eraa199
[32] HUAI B, BAI M, TONG P P, et al. CgPBA1 may be involved in nuclear degradation during secretory cavity formation by programmed cell death in Citrus grandis ‘Tomentosa’ fruits[J]. Plant Physiology and Biochemistry, 2021, 160: 306-314. doi: 10.1016/j.plaphy.2021.01.018
[33] LIANG M, BAI M, WU H. Zn2+-dependent nuclease is involved in nuclear degradation during the programmed cell death of secretory cavity formation in Citrus grandis ‘Tomentosa’ fruits[J]. Cells, 2021, 10(11): 3222. doi: 10.3390/cells10113222
[34] TONG P, HUAI B, CHEN Y, et al. CisPG21 and CisCEL16 are involved in the regulation of the degradation of cell walls during secretory cavity cell programmed cell death in the fruits of Citrus sinensis (L. ) Osbeck[J]. Plant Science, 2020, 297: 110540. doi: 10.1016/j.plantsci.2020.110540
[35] 仝盼盼. 化橘红和甜橙果实分泌囊发育过程中细胞壁降解的细胞和分子特性研究[D]. 广州: 华南农业大学, 2020. [36] CUTTER E G. Plant anatomy: Experiment and interpretation part I: Cells and tissues[M]. London: Addison Wesley, 1978: 106-143.
[37] FAHN A, RACHMILEVITZ T. Ultrastructure and nectar secretion in Lonicera japonica[J]. Botanical Journal of the Linnean Society, 1970, 63: 51-56.
[38] CHATIN J. Etudes histologiques et histogeniques sur les glandes foliairesinterieures et quelques productions analogues[J]. Annales Des Sciences Naturelles. Botanique, 1875, 2(6): 199-221.
[39] PETERSON R L, SCOTT M G, ELLS B E. Structure of astem-derived callus of Rutagravelens; meristems, leaves, and secretory structures[J]. Canadian Journal of Botany-Revue Cannadiennede Botanique, 1978, 56: 271-283.
[40] CARR D J, CARR S G M. Oil glands and ducts in Eucalyptus L'Herit. II: Development and structure of oil glands in the embryo[J]. Australian Journal of Botany, 1970, 18(2): 191-212. doi: 10.1071/BT9700191
[41] FREY-WYSSLING, A. Elimination processes in higher plants[J]. Saussurea, 1972, 3: 79-90.
[42] 胡正海, 余刚. 枳分泌囊的结构与发育的研究[J]. 植物学报, 1993, 35(6): 447-452. [43] LIANG S J, WANG H Y, YANG M, et al. Sequential actions of pectinases and cellulases during secretory cavity formation in Citrus fruits[J]. Trees, 2009, 23(1): 19-27. doi: 10.1007/s00468-008-0250-7
[44] BOZHKOV P V, LAM E. Green death: Revealing programmed cell death in plants[J]. Cell Death & Differentiation, 2011, 18(8): 1239-1240. doi: 10.1038/cdd.2011.86
[45] VAN HAUTEGEM T, WATERS AJ, GOODRICH J, et al. Only in dying, life: Programmed cell death during plant development[J]. Trends in Plant Science, 2015, 20(2): 102-113. doi: 10.1016/j.tplants.2014.10.003
[46] SPERRY J S. Evolution of water transport and xylem structure[J]. International Journal of Plant Sciences, 2003, 164(S3): S115-S127. doi: 10.1086/368398
[47] ESCAMEZ S, TUOMINEN H. Programmes of cell death and autolysis in tracheary elements: When a suicidal cell arranges its own corpse removal[J]. Journal of Experimental Botany, 2014, 65(5): 1313-1321. doi: 10.1093/jxb/eru057
[48] YAMAUCHI T, SHIMAMURA S, NAKAZONO M, et al. Aerenchyma formation in crop species: A review[J]. Field Crops Research, 2013, 152: 8-16. doi: 10.1016/j.fcr.2012.12.008
[49] GUNAWARDENA A H L A N. Programmed cell death and tissue remodelling in plants[J]. Journal of Experimental Botany, 2008, 59(3): 445-451. doi: 10.1093/jxb/erm189
[50] GUNAWARDENA A H L A N, GREENWOOD J S, DENGLER N G. Programmed cell death remodels lace plant leaf shape during development[J]. Plant Cell, 2004, 16(1): 60-73. doi: 10.1105/tpc.016188
[51] LORD C E N, DAUPHINEE A N, WATTS R L, et al. Unveiling interactions among mitochondria, caspase-like proteases, and the actin cytoskeleton during plant programmed cell death (PCD)[J]. PLoS One, 2013, 8(3): e57110. doi: 10.1371/journal.pone.0057110
[52] PENNELL R I, LAMB C. Programmed cell death in plants[J]. The Plant Cell, 1997, 9(7): 1157-1168.
[53] LOMBARDI L, CASANI S, CECCARELLI N, et al. Programmed cell death of thenucellus during Sechium edule Sw. seed development is associated with activation of caspase-like proteases[J]. Journal of Experimental Botany, 2007, 58(11): 2949-2958. doi: 10.1093/jxb/erm137
[54] CARLILE G W, SMITH D H, WIEDMANN M. Caspase-3 has a nonapoptotic functionin erythroid maturation[J]. Blood, 2004, 103: 4310-4316. doi: 10.1182/blood-2003-09-3362
[55] HATSUGAI N, IWASAKI S, TAMURA K, et al. A novel membrane fusion-mediated plant immunity against bacterial pathogens[J]. Genes and Development, 2009, 23(21): 2496-2506. doi: 10.1101/gad.1825209
[56] ZHANG L R, XU Q, XING D, et al. Real-time detection of caspase-3-like proteases activation in vivo using fluorescence resonance energy transfer during plant programmed cell death induced by ultraviolet-C overexposure[J]. Plant Physiology, 2009, 150(4): 1773-1783. doi: 10.1104/pp.108.125625
[57] CASANI S, FONTANINI D, CAPOCCHI A, et al. Investigation on cell death in themegagametophyte of Araucaria bidwillii Hook. post-germinated seeds[J]. Plant Physiology and Biochemistry, 2009, 47(7): 599-607. doi: 10.1016/j.plaphy.2009.02.012
[58] HAN J J, LIN W, ODA Y, et al. The proteasome is responsible for caspase-3-like activity during xylem development[J]. The Plant Journal, 2012, 72(1): 129-141. doi: 10.1111/j.1365-313X.2012.05070.x
[59] SUGIYAMA M, ITO J, AOYAGI S, et al. Endonucleases[J]. Plant Molecular Biology, 2000, 44(3): 387-397. doi: 10.1023/A:1026504911786
[60] ZHENG P, BAI M, CHEN Y, et al. Programmed cell death of secretory cavity cells of Citrus fruits is associated with Ca2+ accumulation in the nucleus[J]. Trees, 2014, 28(4): 1137-1144. doi: 10.1007/s00468-014-1024-z
-
期刊类型引用(2)
1. 徐慧娟，邝智祥，左珂菁，吴承洋，郑泓昊，谢青梅，陈伟国. tva受体基因起始密码子突变对清远麻鸡感染K亚群禽白血病病毒的影响. 华南农业大学学报. 2025(03): 351-357 . 
本站查看
2. 郑文静，田亚东，李转见，李东华. 禽白血病抗病育种研究进展. 中国畜禽种业. 2024(01): 42-46 . 百度学术
其他类型引用(0)
下载:
计量
- 文章访问数: 121
- HTML全文浏览量: 9
- PDF下载量: 355
- 被引次数: 2
