Characteristics of and Factors Affecting Late-stage Inflorescence on Pear Trees
-
摘要:目的 为加深对梨晚期花利用价值认识和为开展梨晚期花形成机理及促生技术研究提供借鉴，方法 采取定点跟踪调查方法，对梨晚期花发生时间、发生部位、花果特征特性及影响发生因素，展开较为系统观察调查分析，结果 梨晚期花发生于春季开花后至9月之间，其中果实迅速膨大期、果实缓慢增大期、果实迅速增大期、贮藏营养积累期分别占70.72%、16.89%、7.93%、4.46%；发生部位包括当年抽生果台副梢、中长梢和多年生枝不定芽梢顶端，占比分别为94.64 %、4.29%、1.07%。与梨正常花相比，晚期花序轴较长，似呈基部多带3～4片叶伞房花序；花着生间距大、花柄长、花托小；开花顺序为中花先开，边花后开。梨晚期花具有一定结果能力，‘黄花’梨发生于果实迅速膨大期、果实缓慢增大期、果实迅速增大期、贮藏营养积累期的坐果率分别为43.23%、16.61、29.44%、29.48%，尤其是早期发生于果台副梢上的果实可正常成熟为可食用的商品果。不同品种间比较，以花芽易形成且开花期长的品种发生晚期花比率高、数量多；不同年际间比较，以在开花幼果期气候或病虫等灾害较重的年份其发生率高、数量多；不同树势比较，以树势强且树体贮藏营养丰富的植株发生率高、数量多。结论 梨晚期花发生时间与发生部位分别以果实迅速膨大期和当年抽生的果台副梢占比大；花序轴较长、花朵着生间距大、花柄长、花托小，开花顺序与正常花相反；具有一定的结果能力，早期发生于果台副梢上的所结的果能正常成熟，对开展灾害性产量损失补救及果实熟期调节具有较大的潜在研发利用前景；大量促生的措施应以丰富树体营养为基础，应用环境胁迫刺激为诱导。Abstract:Objective Regularity and characteristics of as well as factors affecting the late-stage inflorescence of pear trees relating to fruit bearing of the plants were studied.Method Occurrence time and location of late-stage flowers on pear trees as well as physiological and environmental factors affecting the inflorescence were systematically monitored in the field.Result The late-stage flowers started to appear on a pear tree 10 days or more after regular inflorescence in spring and lasted till September. They showed up 70.27% of the time during the period when pears formed from regular flowers were rapidly expanding on a plant, 16.89% when the fruits were slowly enlarging, 7.93% when the fruits were rapidly enlarging, and 4.46% when the fruits were storing and accumulating nutrients. They were found mostly on tip of annual auxiliary shoots (94.64% of total), some on tip of mid-long shoots at 4.29%, and a few on tip of perennial adventitious buds at 1.07%. Compared with the regular flowers, the late comers were longer on the inflorescence axis with 3 or 4 small leaves on the base, bigger on the spacing between individual flowers, lengthier on the stalks, and smaller on the receptacle. They also started blooming from the center rather than the peripheral of a inflorescence. Of all, 43.23% of the fruit-setting from the late-stage flowers occurred in the rapid expanding period of regular fruits, 16.61% in the slow fruit enlarging period, 29.44% in the rapid fruit enlarging period, and 29.48% in time of nutrient storing and accumulating on Huanghua pear plants. The late-stage flowers found on the auxiliary shoots would generally develop and form fruits that ripened on the tree. Such result was observed on the pear plants with ample nutrient supply, of the varieties that showed a high rate of bud formation and long florescence period, and/or during the years encountering stresses such as serious diseases, pest infestation and/or climatic changes that affected fruit development.Conclusion The late-stage flowers on a pear tree were largely found in the rapid expanding time of the regular fruits and on the tip of annual auxiliary shoots. They had longer inflorescence axis, wider spacing between flowers, lengthier stalks, and smaller receptacles than the regular flowers on a same plant. Their blooming sequence was opposite to the regular inflorescence. The late-stage flowers on the auxiliary shoots that bloomed early generally developed fully into fruits that matured into edible product for marketing. At time of yield decline due to diseases, the late fruit bearing compensated somewhat the loss. The delay could also inadvertently alter the fruit harvesting. Therefore, if late-stage inflorescence was desired, applying fertilization coupled with external stress to induce changes on the plant physiology could be a plausible approach.
-
Keywords:
- Pear /
- late-stage inflorescence /
- occurrence regularity /
- characteristics /
- influence factors
-
0. 引言
【研究意义】 蝴蝶兰（Phalaenopsis equestris）花型独特、花色淡雅、花期持久，具有较高的观赏价值和市场经济价值，是我国销售量最大的盆花花卉，是目前市场上非常有竞争力和发展潜力的一种价值较高的花卉。但是病毒性病害严重影响和危害了蝴蝶兰的发展前景，已经成为蝴蝶兰产业持续、健康和快速发展的一个巨大绊脚石。福建漳州蝴蝶兰花卉产业在全国具有举足轻重的重要地位。福建漳州作为全国最大的兰花培育集散中心，年产出几千万盆兰花，其中蝴蝶兰约500万～800万盆，投放市场的成品兰占全国兰花市场份额的10%，同时也已成为向美国、日本、韩国和欧洲等国出口蝴蝶兰的重要基地[1]。植物病毒病被称“植物癌症”，一旦发病，尚无有效方法挽救，防治非常困难。近几年，由于全球气候变化、生态环境恶化、农业种植结构、耕种模式改变等原因，蓟马和蚜虫等虫媒传播的花卉病毒病日益严重，对我国蝴蝶兰等重要花卉的安全生产造成了严重的威胁。因此，亟需加强病毒-蝴蝶兰的互作研究，从而为制定防治策略提供理论依据。【前人研究进展】目前，研究表明感染蝴蝶兰植株的病毒种类较多，主要包括建兰花叶病毒（cymbidium mosaic virus, CymMV）[2]、齿兰环斑病毒（odontoglossum ringspot virus, ORSV）[3]、黄瓜嵌纹病毒（cucumber mosaic virus, CMV）[4]、番茄斑萎病毒（tomato spotted wilt virus, TSWV）[5]、建兰环斑病毒（cymbidium ringspot virus, CymRV）、辣椒褪绿病毒（capsicum chlorosis virus, CaCV）[6]、凤仙花坏死斑病毒（impetiens necrotic spot virus, INSV）[7]、兰花斑点病毒（ovrchid fleck rhabdovirus, OFV）[8]、蝴蝶兰黄化斑点病毒 （phalaenopsis chlorosisspot virus, PhCSV）[9]、康乃馨斑驳病毒（carnation mottle virus, CaMV）[10]和石斛兰叶脉坏疽病毒（dendrobium veinnecrosis closterovirus, DVNV） [11]等。其中CymMV和ORSV是蝴蝶兰生产中的主要病毒，且两者经常复合侵染蝴蝶兰植株。CymMV和ORSV分别属于甲型线形病毒科（Alphaflexiviridae） 和帚状病毒科（Virgaviridae）的正单链RNA病毒[12]。病毒感染蝴蝶兰会影响植株的生长发育，具体表现为生长速度减缓，叶片出现规则或者不规则的褪绿黄化、炭疽坏死斑块，花器出现畸形、色泽不均的斑块、提早凋萎等不同症状[13−14]，对蝴蝶兰影响巨大。RNA干扰（RNA interference, RNAi）是植物、昆虫和真菌等生物保守的抵御病毒等病原物侵染的免疫防御机制[15−16]。研究表明，病毒侵染模式植物拟南芥后，植物的Dicer-like (DCL)蛋白识别并切割长的病毒来源的双链基因组为21、22、24 nt等不同长度的小干扰RNA（small interference RNA, siRNA），这些病毒来源的siRNA（vsiRNAs）与Argonaute结合形成RNA诱导的沉默复合体，随后复合物中vsiRNAs通过碱基匹配进行基因定位，Argonaute蛋白通过核酸切割活性沉默基因[17]，从而沉默外源病毒等病原物基因的表达，由此可知，RNAi免疫途径被激活或者病毒与寄主植物发生互作的一个典型特征是病毒侵染寄主植物细胞后能够诱导产生大量病毒来源的vsiRNAs。由于病原系统的差异，病毒vsiRNAs在丰度、长度分布、5′端碱基偏好性、正负义来源、基因组高低频切割位点等方面呈现广泛的多样性[15]。不同种的病毒侵染不同寄主细胞后产生的 vsiRNAs 丰度存在较大差异，丰度一般在0.1%～20%[15]。vsiRNAs的长度分布主要集中在21 nt和22 nt[15]。vsiRNAs 5′端碱基偏好性根据vsiRNAs结合的AGO蛋白而有所不同[15]。vsiRNAs近乎等比例来源于正负义基因组或者主要来源于正义或负义基因组[15]。vsiRNAs在病毒基因组上呈现连续但不均匀的分布，即形成高低频切割位点[15]。 【本研究切入点】运用深度测序技术，已经在众多不同的植物细胞中鉴定到大量病毒来源的vsiRNAs，vsiRNAs的鉴定已成为研究病毒和寄主互作的有效途径或方法[18−19]。然而关于蝴蝶兰与病毒病原之间的互作尚有待深入探讨。【拟解决的关键问题】通过电镜和PCR等技术鉴定CymMV和ORSV感染蝴蝶兰，并通过深度测序技术鉴定分析病毒来源的vsiRNAs特征来分析CymMV、ORSV和蝴蝶兰寄主之间的相互作用，以期为防治蝴蝶兰病毒病提供依据。
1. 材料与方法
1.1 试验材料
病毒侵染和健康的蝴蝶兰样品均由福建漳州台湾农民创业园绿乐园艺有限公司提供。多糖多酚植物总RNA提取试剂盒（RNAprep pure plant plus kit）、一步法RT-PCR试剂盒（fastking one step RT-PCR kit）和DNA Marker分子量标准（MD102）均购自天根生化科技（北京）有限公司。
1.2 RT-PCR
根据CymMV和ORSV外壳蛋白基因序列（Accession登录号分别为: AB693982和AB693988）设计特异扩增引物CymMV-F：ATGGGAGAGCCCACTCCAACT，CymMV-R：TATCCTCCTGGAAACCAGCCT；ORSV-F：ATGTTTTACACTATTACAGA，ORSV-R：TGAGACTTGATTGTACATACCA，两个基因片段大小均约400 bp。用多糖多酚植物总RNA提取试剂盒（RNAprep pure plant plus kit）提取病毒侵染和健康蝴蝶兰叶片总RNA后，根据一步法RT-PCR试剂盒（fastking one step RT-PCR kit）说明进行RT-PCR检测。一步RT-PCR体系（50 μL）为：2×fastking one step RT-PCR mastermix 25 μL，25×RT-PCR Enzyme Mix 2 μL，上下游引物（10 μmol·L−1）各1.25 μL，RNA模版0.01～1.00 μg，ddH2O补至50 μL。RT-PCR程序：42 ℃，30 min；95 ℃，3 min；94 ℃，30 s，55～65 ℃，30 s，72 ℃，30 s，35个循环；72 ℃，5 min。RT-PCR扩增产物用琼脂糖凝胶电泳检测。
1.3 电镜观察
电镜负染观察主要步骤：获取病毒侵染和健康蝴蝶兰叶片的汁液，将汁液滴加在覆有Formvar膜的铜网上，2 min后用滤纸吸去铜网膜上多余的液体，用1%磷钨酸（pH 7.2）进行负染，用透射电子显微镜H-
7650 Hitachi（日立）观察病毒粒体的形态和大小等特征。电镜超薄切片观察主要步骤：病毒侵染和健康蝴蝶兰叶片切成0.5 cm的小块后，依次用2.5%戊二醛进行固定，用酒精梯度脱水，用Spurr试剂进行渗透，用环氧树脂进行包埋；随后在聚合器中70 ℃进行聚合；聚合后的样品在超薄切片机Leica UC7（德国莱卡）上进行切片制备；超薄切片转移至铜网进行染色；最后将铜网置于透射电镜H-
7650 Hitachi （日立）观察病毒粒体的形态和大小等特征。1.4 siRNAs深度测序和特征分析
参照Lan的方法[18]进行siRNAs的深度测序和鉴定分析。主要步骤为：用多糖多酚植物总RNA提取试剂盒（RNA preppure plant plus kit）提取病毒侵染和健康蝴蝶兰叶片总RNA，10 μg总RNA经8%聚丙烯酰胺凝胶电泳后，回收纯化18～32 nt的siRNAs，用T4 RNA酶分别连接3′接头和5′接头后，SuperScript Ⅱ酶反转录后进行PCR反应，通过聚丙烯酰胺凝胶电泳回收纯化PCR产物，Illumina HiSeq 2000测序仪（Illumina，美国）进行深度测序。用软件Bowtie将siRNAs与CymMV和ORSV基因组进行匹配，再进行病毒来源vsiRNAs丰度、长度、碱基偏好性和正负义链分布等特征分析。
2. 结果与分析
2.1 CymMV和ORSV在蝴蝶兰中的侵染
采集蝴蝶兰感病样品（图1A），以蝴蝶兰叶片总RNA为模板，CymMV-F/R和ORSV-F/R为引物，RT-PCR扩增两种病毒的外壳蛋白基因部分片段，经琼脂糖凝胶电泳检测。结果显示，RT-PCR扩增获得了两条特异的目的片段，片段大小约400 bp（图1B），与引物设计预期结果一致，初步表明CymMV和ORSV两种病毒侵染蝴蝶兰。
![]()
图 1 CymMV和ORSV在蝴蝶兰叶片的共侵染鉴定A：病毒侵染的蝴蝶兰样品；B：CymMV和ORSV外壳蛋白基因的RT-PCR扩增产物的琼脂糖凝胶电泳， M为DNA marker MD102；1和2分别为 CymMV和ORSV外壳蛋白基因产物； C：CymMV和ORSV病毒电镜观察的负染图；D：CymMV和ORSV病毒电镜观察的超薄切片图。Figure 1. Identification of CymMV/ORSV co-infections in P. equestrisA: P. equestris plant infected by viruses; B: agarose gel electrophoresis of RT-PCR products of CP genes of CymMV and ORSV; M: DNA marker MD102; 1 and 2: RT-PCR products of CP genes of CymMV and ORSV, respectively; C: electron microscopic negative staining images of CymMV and ORSV; D: electron microscopic images on thin sections of CymMV and ORSV.为了进一步验证RT-PCR结果的准确性和可靠性，采集蝴蝶兰感病叶片进行负染和超薄切片制备，然后用电镜观察结果。电镜结果表明，负染和超薄切片中都能够观察到长约300 nm线状的病毒粒子和长约500 nm杆状的病毒粒子（图1C、D），300 nm线状的病毒粒子和500 nm杆状的病毒粒子分别代表CymMV和ORSV两种病毒，在健康蝴蝶兰叶片样品中并未观察到该种大小和形态特征的病毒，因此电镜结果进一步证实CymMV和ORSV两种病毒侵染了蝴蝶兰。
2.2 CymMV和ORSV来源的vsiRNAs特征分析
2.2.1 siRNAs深度测序总体概况
病毒侵染证实后，提取感病和健康植株的总RNA进行siRNAs的深度测序，鉴定和分析病毒来源的vsiRNAs特征。siRNAs深度测序结果显示，健康和病毒侵染的蝴蝶兰植株分别获得7 563 892和6 133 689个读数的原始种类的小RNA（total siRNAs）；将这些total siRNAs匹配到CymMV和ORSV两种病毒基因组后，结果表明，健康和病毒侵染的蝴蝶兰植株中能匹配到病毒基因组的total vsiRNAs读数分别为906和826 682；在病毒侵染的蝴蝶兰样品中，匹配到CymMV和ORSV病毒基因组的total vsiRNAs读数分别为369 524（约44.7%）和457 104（55.3%）。
2.2.2 vsiRNAs的长度特征
获得siRNAs读数后进行特征分析，长度特征分析表明，在健康蝴蝶兰样品中，24、21、22 、23 nt长度的total siRNAs丰度依次较高（图2A），然而在病毒侵染的蝴蝶兰样品中，21、24、22、23 nt长度的total siRNAs丰度依次较高（图2A），total siRNAs最大丰度的差异表明病毒侵染改变了total siRNAs的长度分布形式或特征。另外，在病毒侵染的蝴蝶兰样品中，CymMV和ORSV来源的vsiRNAs长度丰度均主要为21 nt和22 nt（图2B），其中21 nt vsiRNAs丰度均为最高，22 nt vsiRNAs次之，该结果与Pai等的研究结果一致[19]。21 nt vsiRNAs高丰度不仅表明21 nt vsiRNAs在RNAi免疫过程中起主要的抵御病毒作用，还暗示在蝴蝶兰细胞中DCL4和DCL2蛋白可能是vsiRNAs分子形成过程中的主要参与者。vsiRNAs是双链小RNA分子，vsiRNAs 5′或3′端突出的碱基数量称为悬挂（overhang）。对于21 nt的vsiRNAs分子，悬挂的碱基数量=21−碱基匹配数量。21 nt vsiRNAs的overhang分析表明，CymMV和ORSV来源的21 nt的vsiRNAs中，具有2个碱基overhang的vsiRNAs丰度都最高（图2C），暗示在蝴蝶兰细胞中，CymMV和ORSV来源的21 nt的vsiRNAs中具有2个碱基overhang的vsiRNAs可能是RNAi免疫反应最有效的激发者。
![]()
图 2 蝴蝶兰细胞中CymMV和ORSV来源的vsiRNAs特征A：健康蝴蝶兰和病毒侵染蝴蝶兰中siRNAs长度分布；B： 病毒侵染蝴蝶兰中vsiRNAs长度分布； C：病毒侵染的蝴蝶兰中5'端距离为1~21 nt vsiRNAs读数；D： Total vsiRNAs（左）和21 nt vsiRNAs（右） 5'端第一个碱基偏好性。Figure 2. Characterization of CymMV- and ORSV-derived vsiRNAs in P. equestrisA: size distribution of total small RNAs in virus-infected and viruses-free P. equestris plants; B: size distribution of vsiRNAs matching CymMV and ORSV genomes in virus-infected P. equestris plants; C: reads of 21-nt vsiRNAs with 1-21 nt distance between 5′ ends from CymMV and ORSV genome in P. equestris plants; D: relative frequency of 5′-terminal nucleotide of total vsiRNAs (left) and 21 nt vsiRNAs (right).2.2.3 vsiRNAs的碱基偏好性
vsiRNAs 5′端的第一个碱基种类会影响vsiRNAs结合到哪种Argonaute蛋白。Total vsiRNAs 5′端碱基偏好性分析表明，CymMV和ORSV来源的vsiRNAs 5′端碱基偏好性都是依次偏好于U、A、C和G（图2D，左Panel），该结果表明CymMV和ORSV来源的vsiRNAs因第一个碱基的差异其可能结合各自不同的Argonaute蛋白，进而介导RNAi免疫效应。长度丰度分析表明，21 nt vsiRNAs丰度最高，其可能是蝴蝶兰细胞RNAi免疫途径的主要激发者，因此，进一步分析了21 nt vsiRNAs的5′端第一个碱基偏好性分析，分析结果表明，21 nt vsiRNAs的5′端第一个碱基偏好性也是依次偏好于U、A、C和G（图2D，右Panel），该结果进一步暗示vsiRNAs 5′端第一个碱基的类型与其结合的Argonaute蛋白之间具有密切的联系。
2.2.4 vsiRNAs的正负义链来源
total vsiRNAs正负义链来源（极性）分布表明，CymMV来源的vsiRNAs中，分别有
326662 个读数（88.4%）和42 862个读数（11.6%）的vsiRNAs来自于正义链和负义链病毒基因组（图3A）；与此类似，ORSV来源的vsiRNAs中，分别有425 548个读数（93.1%）和31 520个读数（6.9%）的vsiRNAs来自于正义链和负义链病毒基因组（图3A），该结果表明蝴蝶兰细胞中CymMV和ORSV vsiRNAs主要来源都是病毒基因组的正义链。![]()
图 3 CymMV和ORSV vsiRNAs的生成前体分析A：CymMV和ORSV vsiRNAs的正负义链分布；B：vsiRNAs在CymMV基因组上高低频切割位点；C：vsiRNAs在ORSV基因组上高低频切割位点； D：CymMV正链基因组的二级结构预测；E：ORSV正链基因组的二级结构预测。方框I为图3B和图3D中高频切割位点和基因组茎环结构的对应关系。方框II为图3C和图3E中高频切割位点和基因组茎环结构的对应关系。Figure 3. Biogenesis precursors of CymMV and ORSV vsiRNAsA: polarity distribution of vsiRNAs matching CymMV and ORSV genomes from co-infected P. equestris plants; B and C: vsiRNAs hotspots and cold spots distribution along CymMV and ORSV genomes, respectively; D and E: secondary structures of CymMV and ORSV RNAs predicted with RNAfold server. Box I shows the corresponding relationship between the high frequency cutting sites in Figures 3B and the stem ring structure in Figure 3D. Box II shows the corresponding relationship between the high frequency cutting sites in Figures 3C and the stem ring structure in Figure 3E.2.2.5 高低频切割位点及其与基因组结构的关系分析
将total vsiRNAs定位到CymMV和ORSV病毒基因组具体位点后，高低频切割位点分析表明，vsiRNAs沿着CymMV和ORSV病毒基因组呈现连续但不均匀的分布特征（图3B、C）；大部分（95%）total vsiRNAs产生位点位于CymMV和ORSV病毒基因组内部，而不是基因组的5′和3′两端（图3B、3C）；在病毒基因组不同位点形成了高频切割点（hotspot）和低频切割点（coldspot）（图3B、C）。
用RNAfold server在线平台预测CymMV和ORSV病毒基因二级结构，结构预测结果表明在CymMV和ORSV病毒基因内部形成了大量的茎环二级结构（图3D、E），且茎环二级结构的形成与高频切割位点是相对应的，表明CymMV和ORSV病毒基因内部茎环二级结构是RNAi免疫途径关键因子DCL蛋白的切割底物，从而形成不同大小的vsiRNAs，即茎环结构是vsiRNAs的主要生成前体。
3. 讨论
蝴蝶兰是具有重要观赏和经济价值的花卉之一。但是病毒性病害已经成为制约蝴蝶兰产业持续健康发展的较大因素。福建漳州蝴蝶兰花卉产业在全国具有举足轻重的重要地位，但是CymMV和ORSV两种病毒对蝴蝶兰危害严重[12−14]。因此，亟需加强蝴蝶兰病毒性病害尤其是CymMV和ORSV两种病毒的检测调查。本研究采集了漳州漳浦农民产业园区的健康和疑似病毒侵染的蝴蝶兰植株，依次采用RT-PCR和传统的电镜技术明确了漳州蝴蝶兰植株感染了CymMV和ORSV两种主要病毒。
植物病毒素有“植物癌症”之称，目前并无有效的防治办法。为了探索有效的防治措施，务必加强病毒和寄主之间相互作用的研究。RNAi是植物、昆虫和真菌等生物保守的抵御病毒等病原物侵染的免疫防御机制[15−16]。RNAi免疫途径被激活的一个典型特征是病毒侵染寄主植物细胞后能够诱导产生大量病毒来源的vsiRNAs。因此，病毒来源的vsiRNAs的鉴定和特征分析已经成为加深理解病毒和寄主相互作用的重要途径和方法。运用高通量深度技术已经在众多的病毒-植物病原系统中鉴定了大量的vsiRNAs[18]。在本研究中，利用siRNAs深度测序技术在CymMV和ORSV复合侵染的蝴蝶兰叶片中获得了大量（约13.48%）病毒来源的vsiRNAs，该结果表明CymMV和ORSV复合侵染蝴蝶兰叶片细胞后激活了寄主植物RNAi免疫途径，病毒与寄主确实存在相互作用。
病毒与寄主互作的具体机制可以通过分析vsiRNAs特征加以理解。研究表明24 nt siRNA主要在miRNA途径中调控寄主本身的生长发育等生命活动中发挥作用，而21 nt的siRNAs主要在siRNA介导的抗病毒免疫途径中发挥作用[15]。本研究中，高丰度值21 nt siRNAs的出现表明CymMV和ORSV的侵染激活了蝴蝶兰细胞中siRNAs介导的抗病毒免疫途径。在拟南芥等模式植物细胞中，21、22、24 nt的vsiRNAs分别由DCL4、DCL2和DCL1蛋白切割底物产生[17]，因此推测在蝴蝶兰细胞中该长度的vsiRNAs可能也有其同源蛋白切割底物产生。研究表明许多病原系统产生的21 nt的vsiRNAs中，其5′或3′端具有2个碱基的突出（overhang）siRNAs丰度最高[20−22]，与此相似，本研究的蝴蝶兰-病毒系统中，21 nt的vsiRNAs中，5′或3′端具有2个碱基的突出（overhang）vsiRNAs丰度也最高，该结果暗示CymMV和ORSV来源的21 nt的vsiRNAs中具有2个碱基overhang的vsiRNAs可能是RNAi免疫反应最有效的激发者。
研究表明siRNAs 5′端的碱基在介导siRNAs与Argonaute蛋白的结合密切相关[23]。比如，在拟南芥模式植物中，5′端为碱基A的siRNAs 更易与Argonaute2和Argonaute4结合，5′端为碱基U的siRNAs 更易与Argonaute1结合，5′端为碱基C的siRNAs 更易与Argonaute5结合，5′端为碱基G的siRNAs 更易与其他Argonaute结合。本研究中vsiRNAs 5′端碱基偏好性分析表明，CymMV和ORSV来源的vsiRNAs中，5’端碱基偏好性都是依次偏好于U、A、C和G；另外21 nt vsiRNAs 5′端碱基偏好性同样是依次偏好于U、A、C和G，该碱基偏好性结果与其他病原系统的vsiRNAs碱基偏好性一致[24−25]；该结果表明CymMV和ORSV来源的vsiRNAs可能结合各自不同的Argonaute蛋白，进而介导RNAi免疫效应，但是由于目前蝴蝶兰基因组无法获得，也就无法获得蝴蝶兰相关的Argonaute蛋白信息，因此后续获得蝴蝶兰基因组中编码Argonaute的基因序列是一项非常重要的工作。
在不同病毒激活寄主RNAi免疫途径中，vsiRNAs产生的前体各不相同[15]。vsiRNAs的形成前体主要通过分析其在病毒基因组正负义链分布、核苷酸产生的位点以及频率分布来推测。在本研究发现，CymMV和ORSV复合侵染的蝴蝶兰细胞中，CymMV和ORSV来源的vsiRNAs主要来自于病毒基因组的正义链，该结果与其他一些正单链RNA病毒在寄主细胞中形成的vsRNAs的正负义偏好性相反，这些病原系统中vsiRNAs大约等比例来源于病毒基因组正负义链，因此他们推测病毒基因组复制过程中形成的双链为vsiRNAs的形成前体[26−29]。进一步将 vsiRNAs定位到CymMV和ORSV病毒基因组具体位点后，高低频切割位点分析表明，vsiRNAs沿着CymMV和ORSV病毒基因组呈现持续但不均匀的分布特征；大部分（95%）total vsiRNAs产生位点位于CymMV和ORSV病毒基因组内部，而不是基因组的5′和3′两端；在病毒基因组不同位点形成了高频切割点和低频切割点。随后用RNAfold server在线平台预测CymMV和ORSV病毒基因二级结构，结构预测结果表明在CymMV和ORSV病毒基因内部形成了大量的茎环二级结构，且茎环二级结构的形成与高频切割位点相对应，因此推测CymMV和ORSV病毒基因组正义链通过碱基匹配自身折叠产生的 pre-miRNA（Precursor miRNA）的茎环结构区域为DCL蛋白的切割底物，即vsiRNAs产生的前体，该预测结果与Flynt和Molnár 等 [30−31] 的研究结果一致。
4. 结论
本研究首先通过电镜和PCR技术鉴定了CymMV和ORSV两种病毒在蝴蝶兰植株中的复合侵染；然后通过深度测序技术鉴定分析了来源于CymMV和ORSV两种病毒基因组的vsiRNAs的丰度、长度、悬挂、碱基偏好性、正负义链分布、Hotspot和Coldspot等特征，从而加深了蝴蝶兰抗病毒RNAi免疫途径防御CymMV和ORSV两种病毒的机制的理解，对研究蝴蝶兰病毒性病害的防治措施具有重要的理论意义。
-
![]()
图 2 ‘翠玉’梨晚期花与正常花特征特性比较
Figure 2. Characteristics of late-stage and regular inflorescences on Cuiyu pear trees
表 1 梨晚期花不同生长发育期发生数量与比例
Table 1 Number and occurrence rate of late-stage inflorescences at pear fruit development stages
年份
Year发生总量
Total number/
（朵·株−1）果实迅速膨大期
Fruit rapidly
expanding period果实缓慢增大期
Fruit slowly
enlargement period果实迅速增大期
Fruit rapidly
enlargement period贮藏营养积累期
Storage nutrition
accumulation period发生量
Number/
（朵·株−1）比例
Ratio/%发生量
Number/
（朵·株−1）比例
Ratio/%发生量
Number/
（朵·株−1）比例
Ratio/%发生量
Number/
（朵·株−1）比例
Ratio/%2010 217.9 153.2 70.31 41.0 18.81 16.4 7.53 7.3 3.35 2011 52.7 38.3 72.68 8.5 16.13 3.3 6.26 2.6 4.93 2012 41.9 26.6 63.48 7.8 18.62 4.5 10.74 3.0 7.16 2013 65.6 46.0 70.12 10.7 16.31 6.1 9.30 2.8 4.27 2014 44.9 31.4 69.93 4.8 10.69 5.1 11.36 3.6 8.02 2015 91.0 71.7 78.79 13.3 14.61 3.8 4.18 2.2 2.42 2016 38.3 23.4 60.10 7.2 18.80 4.6 12.01 3.1 8.09 合计Total 552.3 390.6 70.72 93.3 16.89 43.8 7.93 24.6 4.46 
下载: 导出CSV
表 2 ‘黄花’梨不同部位发生晚期花花序数量与比例
Table 2 Number and occurrence rate of late-stage inflorescences on parts of Huanghua pear tree
年份
Year花序总量
Total number of inflorescences/个当年生果台副梢
Annual auxiliary shoot当年生中长梢
Annual mid-long shoot多年生枝上不定芽梢
Adventitious buds of perennial shoot花序数
Number of inflorescences/个比例
Ratio/%花序数
Number of inflorescencesv/个比例
Ratio/%花序数
Number of inflorescences/个比例
Ratio/%2014 96 92 95.83 3 3.13 1 1.04 2015 116 107 92.24 7 6.03 2 1.73 2016 68 66 97.06 2 2.94 0 0.00 合计 Total 280 265 94.64 12 4.29 3 1.07
下载: 导出CSV
表 3 ‘黄花’不同生育期发生的晚期花坐果率差异比较
Table 3 Fruit-setting rate of late-stage inflorescences on Huanghua pear trees at fruit development stages
时期 Period 坐果率 Fruit setting rate/% 2014年 Year 2014 2015年 Year 2015 2016年 Year 2016 平均值 Mean 果实迅速膨大期 Fruit rapidly expanding period 42.38 47.61 39.70 43.23 Aa 果实缓慢增大期 Fruit slowly enlargement period 18.22 15.48 16.13 16.61 Cc 果实迅速增大期 Fruit rapidly enlargement period 34.14 27.87 26.32 29.44 Bb 贮藏营养积累期 Storage nutrition accumulation period 26.56 33.03 28.86 29.48 Bb 注：不同小写、大写字母分别表示0.05和0.01水平差异显著,采用费舍尔的LSD多重比较，表4-6同。
Note: Different uppercase or lowercase letters represent extreme or significant difference at 0.01 or 0.05 levels. ANOVA Fisher’s LSD test. The same as table 4-6.
下载: 导出CSV
表 4 ‘黄花’梨正常开花果实与晚期花果实品质鉴评测定差异分析
Table 4 Quality differences between Huanghua pears developed from late-stage and regular flowers
果实类型
Fruit type单果重
Weight per
fruit /g果形指数
Fruit shape index萼片残存率
Ratio of remnant sepals /%“公梨”比例
ratio of the male pear /%肉质
Flesh texture可溶性固形
物含量
Soluble solids content /%固酸比
Ratio of total soluble solids and titratable acidity正常果
Normal fruit246.2 Aa 8.96 Bb 42.44 Bb 38.26 Bb 质地较紧、果汁中
Denser with medium juiciness10.67 Bb 60.81:1 Aa 晚期花果
Late-stage flowers’ fruit178.8 Bb 9.98 Aa 87.61 Aa 79.93 Aa 质地松脆、果汁多
Crisp to soft with juicy11.94 Aa 51.54:1 Bb
下载: 导出CSV
表 5 不同梨品种晚期花发生能力差异分析
Table 5 Abilities of late-stage flower formation of different pear varieties
品种
Variety年均发
生率
Annual average rate /%年均株发生率
Annual average rate per plant /%株年均花序发生量
Annual average number of inflorescences per plant株年均花朵发生量
Annual average number of flowers per plant黄花 ‘Huanghua’ 100.00 Aa 92.43 Aa 9.51 Aa 32.35 Aa 翠玉 ‘Cuiyu’ 100.00 Aa 71.61 Bb 3.78 Bb 14.36 Bb 蜜雪梨 ‘Mixue’ 100.00 Aa 42.40 Cc 1.92 Cc 7.53 Cc 翠冠 ‘Cuiguan’ 100.00 Aa 26.82 Ee 1.20 Ee 3.36 Dd 清香 ‘Qingxiang’ 25.00 Bb 0.67 Dd 0.01 Dd 0.02 Ee
下载: 导出CSV
表 6 ‘黄花’梨年际之间晚期花发生能力差异分析
Table 6 Abilities of late-stage flower formation of Huanghua pear trees in different years
年份
Year株发生率
Rate per plant /%株均花序发生量
Average number of inflorescences per plant株均花朵发生量
Average number of flowers per plant2010 100.00 71.21 Aa 273.45 Aa 2011 91.86 8.16 Dd 28.56 Dd 2012 89.53 6.90 Ef 22.43 Ef 2013 95.35 7.58 DEe 27.83 Dd 2014 88.37 8.67 Dd 32.45 Dd 2015 100.00 28.44 Bb 98.69 Bb 2016 93.51 7.69 Dd 26.92 Dd 2017 100.00 18.22 Cc 80.63 Cc
下载: 导出CSV
表 7 ‘黄花’树势强弱不同树体晚期花发生能力差异分析
Table 7 Abilities of late-stage flower formation of Huanghua pear trees with different degrees of vigor
树体类型
Tree type株发生率
Ratio per plant/%
株均花序发生量
Average number
of inflorescences
per plant株均花朵发生量
Average number of flowers per plant树势强 Vigorous tree 95.00 12.25 44.65 树势弱 Weak tree 17.50 0.78 2.68
下载: 导出CSV
-
[1] 戈登(Gordon, R. F.), 乔丹(Jordan, F. T. W.), 狄伯雄等译. 禽病, 第二版[M]. 北京: 农业出版社, 1991. Poultry diseases[M]., 1991. [2] 王晴芳, 周绂, 肖春桥, 等. 砂梨资源返青返花调查及其预防 [J]. 中国南方果树, 2000, 29(2):44. DOI: 10.3969/j.issn.1007-1431.2000.02.043 WANG Q F, ZHOU F, XIAO C Q, et al. Investigation and prevention of reflorescence and refoliation in Pyrus pyrifolia Nakai [J]. South China Fruits, 2000, 29(2): 44.(in Chinese) DOI: 10.3969/j.issn.1007-1431.2000.02.043
[3] 黄新忠, 张长和, 曾少敏, 等. 梨叶部病害与梨叶早衰脱落相关性分析 [J]. 植物保护, 2015, 41(4):160−164, 183. DOI: 10.3969/j.issn.0529-1542.2015.04.031 HUANG X Z, ZHANG C H, ZENG S M, et al. Correlation analysis between pear leaf diseases and early-senescence and defoliation [J]. Plant Protection, 2015, 41(4): 160−164, 183.(in Chinese) DOI: 10.3969/j.issn.0529-1542.2015.04.031
[4] 黄新忠, 雷龑, 陈小明, 等. 2010年福建省梨树霜害调查 [J]. 中国果树, 2010(4):57−59, 80. HUANG X Z, LEI Y, CHEN X M, et al. Investigation on frost-bite in pear trees in Fujian province in 2010 [J]. China Fruits, 2010(4): 57−59, 80.(in Chinese)
[5] 郗荣庭. 果树栽培学总论(第三版)[M]. 北京: 中国农业出版社, 1997: 64-67. [6] 周俊辉. 果树的多次开花结实现象 [J]. 江西农业学报, 1999, 11(4):64−68. ZHOU J H. Unseasonable flowering and bearing phenomenon of fruit trees [J]. Acta Agriculturae Jiangxi, 1999, 11(4): 64−68.(in Chinese)
[7] 陈志参. 黑珍珠莲雾在福建沿海的反季节生产技术 [J]. 福建农业科技, 2014(8):51−52. DOI: 10.3969/j.issn.0253-2301.2014.08.023 CHEN Z C. Off-season cultivation techniques of wax apple cv. “Heizhenzhu” in Coast area of Fujian Province [J]. Fujian Agricultural Science and Technology, 2014(8): 51−52.(in Chinese) DOI: 10.3969/j.issn.0253-2301.2014.08.023
[8] 庄泽. 毛叶枣 [J]. 中国水土保持, 1991(3):39−40. ZHUANG Z. Zizyphus mauritiana [J]. Soil and Water Conservation in China, 1991(3): 39−40.(in Chinese)
[9] 梁广成, 张志群, 曾妮. 新加坡甜杨桃的特征特性及栽培技术 [J]. 广东农业科学, 2008, 35(4):90−91. DOI: 10.3969/j.issn.1004-874X.2008.04.033 LIANG G C, ZHANG Z Q, ZENG N. Characteristics and Cultivation Techniques of Sweet Star Fruit in Singapore [J]. Guangdong Agricultural Sciences, 2008, 35(4): 90−91.(in Chinese) DOI: 10.3969/j.issn.1004-874X.2008.04.033
[10] 肖邦森, 谢江辉, 雷新涛, 等. 番荔枝优质高效栽培技术[M]. 2001, 北京: 中国农业出版社: 19－23. [11] 林有兴. 果树产期调节技术及其在热带果树上的应用 [J]. 热带农业科技, 2004, 27(1):35−39, 54. DOI: 10.3969/j.issn.1672-450X.2004.01.012 LIN Y X. Adjusting technique of fruiter confinement and the application on the tropical fruiter [J]. Tropical Agricultural Science & Technology, 2004, 27(1): 35−39, 54.(in Chinese) DOI: 10.3969/j.issn.1672-450X.2004.01.012
[12] 邵则夏, 陆斌. 玫瑰香葡萄的一年多熟栽培 [J]. 云南农业科技, 1995(1):16−17. SHAO Z X, LU B. Multi-harvest a year of grape (cv. Muscat Hamburg) [J]. Yunnan Agricultural Science and Technology, 1995(1): 16−17.(in Chinese)
[13] 施金全, 王道平, 李以训, 等. 巨峰葡萄促成避雨一年两熟栽培技术 [J]. 农业科技通讯, 2013(1):135−137. DOI: 10.3969/j.issn.1000-6400.2013.01.056 SHI J Q, WANG D P, LI Y X, et al. Technology of rain-shelter, forcing and annual two-maturity cultivation of Kyoho grape [J]. Bulletin of Agricultural Science and Technology, 2013(1): 135−137.(in Chinese) DOI: 10.3969/j.issn.1000-6400.2013.01.056
[14] 陈英敏, 徐国辉, 王贺新, 等. 不同蓝莓品种花芽形成及二次开花特征 [J]. 北方果树, 2016(4):6−10. [15] 肖新超, 金群英, 彭华正, 等. 竹子分子生物学研究进展 [J]. 竹子研究汇刊, 2007, 26(2):1−5. XIAO X C, JIN Q Y, PENG H Z, et al. Advances in the research of bamboo molecular biology [J]. Journal of Bamboo Research, 2007, 26(2): 1−5.(in Chinese)
[16] 王小红. 环境因子对竹子开花影响研究 [J]. 四川动物, 2009, 28(4):618−621. WANG X H. Study on the influence of environmental factors on bamboo bloom [J]. Sichuan Journal of Zoology, 2009, 28(4): 618−621.(in Chinese)
[17] 黄新忠, 陈义挺, 雷龑, 等. 福建梨早期大量落叶诱因与防控策略 [J]. 中国农学通报, 2010, 26(2):91−95. HUANG X Z, CHEN Y T, LEI Y, et al. Causes and control strategies of a large number of early falling leaves of pear in Fujian [J]. Chinese Agricultural Science Bulletin, 2010, 26(2): 91−95.(in Chinese)
[18] 林庆光, 陶挺燕, 符碧海. 莲雾产期调控技术 [J]. 热带农业工程, 2013, 37(3):25−27. LIN Q G, TAO T Y, FU B H. A sketch of production regulation in wax-apples [J]. Tropical Agricultural Engineering, 2013, 37(3): 25−27.(in Chinese)
[19] 黄胜, 陈雄进. 促进莲雾提早开花调节产期措施研究 [J]. 热带农业工程, 2016, 40(2):38−41. HUANG S, CHEN X J. Measures to advance florescence and regulate yield of Syzygium samarangense [J]. Tropical Agricultural Engineering, 2016, 40(2): 38−41.(in Chinese)
[20] 武萍萍, 周碧燕. 杨桃新梢花芽分化及其碳水化合物含量的变化 [J]. 园艺学报, 2007, 34(5):1151−1156. DOI: 10.3321/j.issn:0513-353x.2007.05.014 WU P P, ZHOU B Y. Flower bud differentiation and changes of carbohydrate in new shoots in Averrhoa carambola [J]. Acta Horticulturae Sinica, 2007, 34(5): 1151−1156.(in Chinese) DOI: 10.3321/j.issn:0513-353x.2007.05.014
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 981
- HTML全文浏览量: 419
- PDF下载量: 31
- 被引次数: 1
