苹果(
Malus pumila Mill.)是蔷薇科苹果属落叶乔木,喜光、耐寒
[1]。中国现存苹果是从美国和日本引进的西洋品种
[2]。我国苹果产区分布极广,主要集中在环渤海、西北黄土高原、黄河故道和西南冷凉高地四大产区
[3]。‘长富2号’是1980年引入我国的着色系富士品种,经过20多年发展,现在已成为我国苹果主栽品种之一,果实综合品质优良。
Askenaly于1867年研究了酸樱桃的花芽分化时期,是果树方面最早的花芽分化研究,其后许多学者从形态、解剖和生理生化等多方面对果树的花芽分化进行了研究,而苹果花芽分化的研究也有100多年的历史。苹果花芽属于混合芽,即展开后可以成为叶芽或花芽
[4]。花芽数量、质量是衡量苹果树体结果潜力的关键、核心指标。果树生长包括营养生长和生殖生长,两者既相互促进,又相互抑制。花芽分化是从营养生长到生殖生长的一个重要的变化过程,这个过程是营养生长结束和生殖生长开始的标志。这种变化极大地影响了花期和花的数量
[5]。花芽分化的实质是植物经一定时期的营养生长后,由于受到光照、温度等外界因素的影响,产生某些与花芽分化相关的激素,在适宜的温度、光照和营养条件下,茎生长锥不再向叶芽方向转变形成叶原基和腋芽原基,而是发生花原基和花序原基,逐渐形成花芽
[6]。果树不同发育阶段的成花能力不同,其枝条生长状况、叶片质量状况都有所差异,这些差异的形成决定了其具有不同的成花能力
[7]。有研究发现,温度、光照、水分、矿物质含量等是影响花芽分化的关键因素
[5]。
甘肃陇东黄土高原区是我国红富士等晚熟苹果栽培的最佳适宜区,苹果产业已成为当地经济发展和乡村振兴的主导产业。该产区20世纪90年代至21世纪初栽植的苹果园大多数采用乔化密植栽培,虽然乔化密植栽培曾经起到过促进丰产的积极作用,但随着树龄的增长,果园郁闭程度的增加,出现枝量繁多、冠层光照条件恶化、病虫害防治难度增大、结果部位外移,年际结果不均衡,产量和果实品质下降等诸多问题
[8-10]。严重制约着苹果产业的进一步发展。为此,乔化成龄密闭苹果园的改造和树体结构的合理调整研究已成为当地苹果栽培中的热点课题。
果园间伐是乔化成龄密闭苹果园改造的主要措施,能够在较短时间内改善果园通风透光条件、优化果园群体结构、促进花芽分化、改善果实内外品质
[11-12],但是随着时间推移,果园密闭现象将再次显现,既不能持久拥有良好的通风条件,也不能持续具备合理的群体结构
[13-14]。缩小树冠、疏除大枝、提高主干、落去主头等改形措施,也可暂时缓解单个树体密闭的问题,进而提高果实内外品质,但也不能永久解决果园整体郁闭问题,不能持久具备优质高产的潜能
[15]。目前有研究发现,对老果园进行间伐改形后,花芽分化明显,成花率显著提高,第三年坐果率也大幅提高
[16],对密闭苹果园进行间伐改形后发现成花能力和花芽质量明显提高,冠层花芽和果实的分布空间也得到改善
[17]。
植物激素是调节高等植物发育和环境信号应答的重要分子,其在植物生长的各个阶段发挥作用,花芽分化亦不例外
[18]。研究发现生长素、脱落酸、多效唑、赤霉素等激素都能够促进花芽分化
[19-20]。有研究表明,高水平含量的脱落酸有利于花器官形成发育;反式玉米素核苷含量升高对菠萝花芽分化的整个过程都有促进作用,与花原基形成的数量正相关;低浓度的生长素和赤霉素有利于藤本月季“安吉拉”的花芽分化
[21]。本试验以陇东高原雨养区十五年生‘长富2号’红富士密闭果园为对象,通过间伐与树体改形相结合的改造试验,研究了不同处理对树体花芽分化的影响,通过不同分化时期的生长类相关激素和转录组分析,在激素和分子生物学层面解析间伐改形对乔化成龄密闭苹果园树体成花的影响,为乔化成龄密闭红富士苹果园合理改造和富士苹果品种的改良提供理论依据。
1 材料与方法
1.1 试验区概况和试验材料
试验于2022年在国家苹果产业技术体系平凉苹果综合试验站的静宁示范县城川镇试验园(35°24′16″ N,105° 47′15″ E)进行。试验园海拔1592 m,年降雨量460 mm,春季少雨较为干旱,秋季雨量丰沛且较集中,年无霜期160 d,年日照时数2242 h,平均气温7.21 ℃。土壤为黄绵土,粘性较弱,有机质含量为9.23 g·kg-1,pH为8.02。试验园面积为3620 m2,供试品种为十五年生‘长富2号’红富士苹果树,砧木为山定子(Malus baccata),行株距为4 m×3 m,密度为834 株·hm-2,改良纺锤形树形,树冠交接率大于140%,树高4.2~4.6 m,主干高度60~75 cm,主枝数6~9个,主枝角度70°~80°,授粉树品种为秦冠和金冠,占比为10%。果园为雨养果园,无灌溉条件,地面自然生草,栽培管理水平中等偏上,果园管理一致,以单个树体为单位进行施肥,即各处理每株树体施肥水平完全一致。
1.2 试验方法
参照牛军强等
[5]的方法,于2022年2月上中旬开始对试验园开展间伐改形改造试验,试验采用单因素随机设计的方法,只做2种处理,分别为试验T(间伐改形)和对照CK(不间伐不改形)。间伐改形采用隔株间伐,行株距由原来的4 m×3 m变为4 m×6 m,同时结合缩小树冠、疏除大枝、提高主干、落去主头等改形措施,使树高变为310~340 cm,干高100~110 cm,选留健康无病虫害、大小合理、分布方位均衡、间距适度的枝条4~6个作为永久主枝,开张角度为90°~100°,其他修剪采用富士苹果常用的手法进行修剪,总体修剪量为30%左右。对照(CK)为不间伐改形(既不间伐也不改形),其它修剪手法和修剪量与间伐改形处理完全相同。
1.3 取样时间与方法
每个处理区面积均为20 m×25 m=500 m2,均有3个重复,每个重复选取树势中庸、无明显病虫害、有代表性的树体5株作为供试树体。于2022年6月22日花芽生理分化前期(CK1和T1)和2022年8月22日花芽生理分化后期(CK2和T2),分别取样1次,取样部位为树体冠层中下部中短枝上的封顶芽体,要求芽体大小、饱满程度一致,以最快速度采取芽体样品,然后用锡箔纸将样品包裹严实并标记清楚,于液氮罐贮藏带回实验室,-80 ℃超低温冰箱保存,等待检测分析。
1.4 转录组学测定方法
用TRIzol (thermofisher, 15596018),根据厂商提供的操作方案对总样品的RNA进行分离和纯化。然后用NanoDrop ND-1000(NanoDrop,美国)对总RNA的量与纯度进行质控并通过Bioanalyzer 2100(Agilent,美国)对RNA的完整性进行检测。使用oligo(dT)磁珠(Dynabeads Oligo (dT), cat.25-61005, Thermo Fisher,美国)通过2轮纯化对其中的带有PolyA(多聚腺苷酸)的mRNA进行特异性捕获。将捕获到的mRNA在高温条件下利用镁离子打断试剂盒(NEBNext Magnesium RNA Fragmentation Module, cat.E6150S,美国)进行片段化,94 ℃保持5~7 min。将片段化的RNA在逆转录酶(Invitrogen SuperScript™ II Reverse Transcriptase, cat.1896649,美国)的作用下合成cDNA。然后使用E. coli DNA polymerase I(NEB, cat.m0209,美国)与RNase H(NEB, cat.m0297,美国)进行二链合成,将这些DNA与RNA的复合双链转化成DNA双链,同时在二链中掺入dUTP Solution(Thermo Fisher, cat.R0133,美国),将双链DNA的末端补齐为平末端。再在其两端各加上1个A碱基,使其能与末端带有T碱基的接头进行连接,并利用磁珠对其片段大小进行筛选和纯化。以UDG酶(NEB, cat.m0280,美国)消化二链,再通过PCR——预变性95 ℃保持3 min,98 ℃变性总计8个循环、每次15 s,退火到60 ℃保持15 s,72 ℃下延伸30 s,最后72 ℃保留5 min,使其形成片段大小为300 bp±50 bp的文库。最后,使用illumina Novaseq™ 6000按照标准操作对其进行双端测序,测序模式为PE150。
1.5 激素含量测定
通过ELISA试剂盒测定赤霉素、生长素、脱落酸含量。
1.6 RNA提取和基因表达分析
根据制造商的说明,使用TRIzol(Vazyme,中国)方法从叶中提取总RNA。使用反转录试盒(TransGen,中国)将总RNA逆转录成cDNA。使用TransStart® Tip Green qPCR SuperMix(TaKaRa,中国)和引物进行qRT-PCR。反应的总体积为20 μL,其中cDNA为200 ng,正向引物、反向引物和ROX Reference Dye均为0.4 μL,TB Green
Premix Ex Taq (Tli RnaseH Plus)(2×)qPCR super mix为10 μL,少于20 μL的部分用双蒸水补充。PCR循环条件设定如下:在94℃下进行30 s的预变性,接着94 ℃下5 s、60 ℃下30 s进行40个循环,最后使用溶解阶段检测荧光信号以验证目标产物是否是单一。qRT-PCR分析中使用的引物如
表1所示,
Actin为内参。
1.7 数据处理
数据采用Excel 2013软件进行统计处理,差异显著性分析使用SPSS 26 统计分析软件。
2 结果与分析
2.1 激素含量测定
对T1、T2、CK1、CK2进行赤霉素、吲哚乙酸、植物脱落酸含量测定,结果如
图1所示,T组的激素含量基本都高于CK组的激素含量。与CK1相比,T1的赤霉素含量、吲哚乙酸含量和脱落酸含量明显增加,与CK2相比,T2赤霉素含量和脱落酸含量显著增加,吲哚乙酸含量稍有下降。这说明,间伐改形能显著提高树体花芽赤霉素、吲哚乙酸和脱落酸等生长相关激素的含量,进而说明间伐改形对花芽分化有一定促进作用。
2.2 转录组学分析
2.2.1 相关性分析
为探索间伐改形对花芽分化的影响,对T1(间伐改形)和CK(不间伐改形)花芽进行转录组测序。样本之间的相关系数越大,样本聚类越好。Pearson相关系数可直观看出样本之间的可重复性,并有助于发现数据异常样本。
图2可见,本研究中样本相关性结果很好,足以支撑后续研究。
2.2.2 质量评估
T1与CK1相比,总共检测48个基因表达上调,54个基因表达下调,T2与CK2相比,总共检测353个基因表达上调,227个基因表达下调,如
图3所示。为了解差异表达基因的整体分布情况,以log(fc)为横坐标,-log(q value)为纵坐标,对差异表达分析中的所有基因绘制火山图,如
图4所示。
2.2.3 DEGs GO富集分析
GO富集分析结果以气泡图的形式呈现。根据Unigene在GO数据库中的比对结果,选取每组样品中差异表达量最高的20个基因进行泡绘,如
图5所示。
2.2.4 生长发育相关基因的挑选
根据GO富集结果,从3 diffExpression文件夹的各比较组中的文件-Gene differential expression中找到如下57个生长发育相关基因(
表2)。由
表2可以看出,gene:MD10G1048600、gene:MD01G1090300、gene:MD09G1176100、gene:MD16G1218900、gene:MD15G1393600等基因在花芽分化前期和后期在生长素相关GO富集通路检测到明显变化;gene:MD10G1048600、gene:MD01G1090300、gene:MD06G1003800、gene:MD15G1012500、gene:MD02G1124900、gene:MD01G1016100等基因在花芽分化前期和后期在脱落酸相关GO富集通路检测到明显变化;gene:MD01G1090300、gene:MD07G1280300、gene:MD15G1012500、gene:MD15G1109300等基因在花芽分化前期和后期在赤霉素相关GO富集通路检测到明显变化。
综合分析后挑选出表达量变化明显的基因5个,分别是:gene:MD10G1048600(生长素、脱落酸相关基因)、gene:MD01G1090300(生长素、脱落酸、赤霉素相关基因)、gene:MD02G1124900(脱落酸相关基因)、gene:MD01G1016100(脱落酸相关基因)、gene:MD07G1280300(赤霉素相关基因)。
2.3 差异表达基因qRT-PCR验证
为了验证转录组测序数据的可信度,将上述挑选出的5个基因:gene:MD10G1048600(NAC030)、gene:MD01G1090300(CCA1)、gene:MD02G1124900(PIP2-7)、gene:MD01G1016100(ABR1)、gene:MD07G1280300(WRKY27),进行qRT-PCR分析。qRT-PCR结果与转录组结果趋势基本一致,表明转录组数据的可信度。与CK1相比,T1处理的MD02G1124900表达量上调,MD10G1048600、MD01G1090300、MD01G1016100、MD07G1280300这4个基因的表达量都下调;与CK2相比,T2处理的MD01G1090300表达量上调,MD10G1048600、MD02G1124900、MD01G1016100、MD07G1280300这4个基因表达量都下调;与T1相比,T2处理的MD02G1124900表达量下调,MD10G1048600、MD01G1090300、MD01G1016100、MD07G1280300这4个基因表达量都上调。其中MD02G1124900上调比较明显。
3 讨论
很多苹果果园目前都存在果园郁闭的问题,果园郁闭会影响通风透光性、降低花芽质量、加大营养消耗、导致大小年现象严重与树势衰弱等
[22-23]。为了改变苹果果园郁闭的现象,通常采用间伐的方式缓解。前人研究表明,间伐可以提高植株叶片的光合速率、促进新梢生长和短枝的成枝、改善果树的枝类组成比例、促进果实着色,提高果面的光洁程度,提高单果重和可溶性固形物含量,可以解决果园郁闭的问题
[24-25]。
植物激素的研究于20世纪30年代因生长素的发现而拉开序幕,目前已基本明确植物激素的生物合成途径
[26]。植物激素是调节植物生长、发育和抗逆性的重要物质,其几乎参与了植物生长发育所有的生理过程
[27]。研究发现低水平生长素、赤霉素(GA3)、茉莉酸甲酯(MeJA)和高水平脱落酸(ABA)有利于花芽分化,且较高的mABA/mGA3和mABA/mMeJA比值有利于花芽分化;吲哚乙酸可以激活花原基发育,有利于花器官的形成;低水平含量的脱落酸有利于花原基形成,促进菠萝进入花芽分化
[28-29]。本试验中对苹果园进行间伐改形后,与对照相比花芽分化后期吲哚乙酸含量降低,花芽分化前期和后期的赤霉素和脱落酸含量及花芽分化前期吲哚乙酸含量均明显提高,说明间伐改形能显著提高生长相关激素的含量,这一结果与前人研究中生长素、赤霉素、脱落酸等激素有利于花芽分化的结果一致,进而说明间伐改形对花芽分化有一定促进作用。
转录组学是从RNA水平研究基因表达的情况,即一个活细胞所能转录出来的所有RNA的总和,是一种重要的生命科学研究手段,被广泛应用到多个物种的研究中
[30]。本试验中苹果园进行间伐改形后进行花芽转录组学分析,以期找到间伐改形促进苹果树体花芽分化的关键基因,通过转录组分析,在苹果花芽分化2个阶段的GO富集分析中找到激素介导的信号通路,包括生长素相关信号通路、对脱落酸的反应、对赤霉素的反应、吲哚乙酸生物合成过程等通路。在上述激素介导的信号通路中发现大量与吲哚乙酸、赤霉素和脱落酸等信号通路的相关基因,从中挑选出5个表达量变化明显的基因,分别是:gene:MD10G1048600(NAC030)、gene:MD01G1090300(CCA1)、gene:MD02G1124900(PIP2-7)、gene:MD01G1016100(ABR1)、gene:MD07G1280300(WRKY27)。对这5个基因进行qRT-PCR分析,其结果与转录组结果趋势基本一致,进一步证实转录组数据的可信度。
间伐改形是常用的提高苹果果实品质的物理方法,但是间伐改形通过哪些基因发挥作用仍未可知。本研究通过转录组分析获得的5个生长发育相关基因,为以后进一步研究如何促进花芽分化,提高苹果果实坐果率,提高果实品质提供理论基础。转基因已经成为现代育种的重要方向,我们下一步将通过对选出的5个基因进行转基因等领域研究,为培育高质量苹果品种做出贡献。
4 结论
综上所述,郁闭苹果园通过隔间伐改形可以有效提高赤霉素、吲哚乙酸和脱落酸等生长相关激素的含量,进而说明间伐改形对花芽分化有一定促进作用。这种作用可能是通过gene:MD10G1048600(NAC030)、gene:MD01G1090300(CCA1)、gene:MD02G1124900(PIP2-7)、gene:MD01G1016100(ABR1)、gene:MD07G1280300(WRKY27)等基因来调控的。
国家自然科学基金(32160683)
国家自然科学基金(31760556)
甘肃省技术创新引导计划项目(22CX8NA026)