Suspension culture dynamics of embryogenic callus from Torreya grandis ‘Merrillii’ and its response to gibberellin

以浙江农林大学省部共建亚热带森林培育国家重点实验室培育的香榧胚性愈伤组织为试验材料。使用SH液体培养基，其中添加脱落酸(ABA) 1.0 mg·L⁻¹、水解酪蛋白(CH) 500.0 mg·L⁻¹、蔗糖30.0 g·L⁻¹、活性炭2.0 g·L⁻¹和GA₃ 0.5 mg·L⁻¹，设置pH 5.7，摇床转速为110 r·min⁻¹，25 ℃暗培养。

在100 mL的锥形瓶中添置30 mL的SH液体培养基(配方同1.1)，香榧胚性愈伤组织接种量为30 g·L⁻¹，设置pH 5.7，转速为110 r·min⁻¹，每3 d取样测定细胞鲜质量及细胞活力(氯化三苯四氮还原法，TTC)，即波长为485 nm处的吸光度。

在100 mL锥形瓶中添置30 mL的SH液体培养基(配方及培养条件同1.1)，设置转接愈伤组织接种量为10、20、30、40、50 g·L⁻¹；pH为5.5、5.6、5.7、5.8、5.9；摇床转速为70、90、110、130、150 r·min⁻¹。25 ℃摇床中振荡暗培养。每个处理设置3个重复，21 d时统计鲜质量，计算细胞增长率，细胞增长率=(收获量−接种量)/接种量×100%。

超氧化物歧化酶(SOD)活性采用氮蓝四唑光化还原法测定；过氧化物酶(POD)活性采用愈创木酚法测定；过氧化氢酶(CAT)活性采用紫外吸收法测定^[15]。

将香榧胚性愈伤组织接种于SH液体培养基中(配方及培养条件同1.1)，设置对照组(ck)和实验组(GA₃+)，对照组培养基中不添加GA₃，实验组培养基中施加GA₃。18 d后收集材料，使用多糖多酚植物总RNA提取试剂盒(北京天根生化科技有限公司)，按说明书提取RNA，分别设置3个生物学重复，使用分光光度计检测浓度与纯度。

测序工作由北京诺禾致源科技有限公司完成，采用Illumina HiSeqTM 2000测序平台进行测序。由于香榧没有参考基因组，将聚类分析结果中最长的转录本作为参考序列单基因簇进行后续分析。

采用RSEM软件将每个样品过滤后的数据(clean reads)在参考序列上做比对。采用DESeq分析差异表达基因(DEG)，以|log₂F_c|＞1 (F_c为表达量的差异倍数)且P＜0.005为筛选条件^[16]。将差异表达基因注释到基因本体论数据库(GO)和京都基因与基因组百科全书(KEGG)，并进行GO和KEGG富集分析。

使用Primer 3.0 input在线软件设计引物(表1)，以TgActin为内参基因^[16]。用Takara生物科技有限公司的反转录试剂将提取的RNA反转录成cDNA进行RT-qPCR验证，每个样品3个生物学重复。反应体系为10.0 μL，包括5.0 μL的TB GREEN，0.2 μL的上下游引物，0.4 μL的模板cDNA，4.2 μL的ddH₂O。PCR反应程序为：95 ℃ 10 min，95 ℃ 10 s，60 ℃ 1 min，95 ℃ 30 s，60 ℃ 15 s共39个循环。

使用Excel处理数据，采用SPSS 26进行单因素方差分析(one-way ANOVA)和多重比较(邓肯法)，显著性水平为0.05，使用GraphPad Prism 9软件作图。

香榧胚性愈伤组织悬浮培养时，细胞鲜质量变化呈“S”型曲线。在0~6 d鲜质量增加缓慢；7~18 d鲜质量快速增加；18 d时达到最高值，为44.81 g；18~21 d进入稳定期，鲜质量基本不再增加；21 d后细胞进入衰退期，鲜质量下降。同时，在培养周期内香榧胚性细胞活力在0~3 d时快速增长，并在3 d时达到最高值(0.32)，之后持续下降，并且随着培养时间的延长，变化趋势趋向平缓(图1)，表明9~12 d为香榧胚性愈伤组织悬浮培养的最佳继代时期。

Figure 1.  Changes in fresh weight and viability of cells during incubation　　　　

随着接种量的增加，胚性愈伤组织增长率呈先上升后下降的趋势。当愈伤组织接种量为10 g·L⁻¹时，细胞增长率最低(21.78%)；当接种量为30 g·L⁻¹时，细胞增长率显著升高(P＜0.05)，达到最高值，为49.94%。SOD、POD和CAT活性均随接种量的增加而呈先升高后降低的趋势。当接种量为30 g·L⁻¹时，SOD、POD和CAT活性均最高，显著高于其他接种量时的酶活性(P＜0.05)(图2A)。

Figure 2.  Effect of different culture conditions on the liquid suspension culture of embryogenic callus of T. grandis ‘Merrillii’

随着pH的增加，香榧悬浮细胞增长率呈先上升后下降的趋势。当培养基pH为5.7时，细胞鲜质量增殖效果最佳，细胞增长率最大，为22.95%，显著高于其他处理(P＜0.05)。SOD、POD和CAT活性均随接种量的增加而呈先升高后降低的趋势。在pH为5.7时，悬浮细胞SOD、POD和CAT活性均最大，显著高于其他处理(P＜0.05)(图2B)。

随着转速的增加，香榧悬浮细胞增长率呈逐渐上升的趋势，当转速为70 r·min⁻¹时，细胞增长率较低，为20.33%；当转速达110 r·min⁻¹时，细胞增长率最大，为50.84%，显著高于其他处理(P＜0.05)。SOD、POD和CAT活性均随摇床转速的增加而呈先升高后降低的趋势。当转速为110 r·min⁻¹时，SOD、POD和CAT活性最大，显著高于其他处理(P＜0.05)(图2C)。

所有检测样品纯净碱基数均达6.66 Gb以上，质量为Q20 (每项碱基质量大于20的碱基数占总碱基数的比例)的碱基比例均大于97.47%，达Q30 (每项碱基质量大于30的碱基数占总碱基数的比例)以上的碱基比例大于92.70%，每个样本碱基GC含量为43.53%~44.45%，较为一致，且测序碱基平均错误率均在0.1%以下，表明转录组数据较为准确。

本研究共筛选到428个差异表达基因，其中236个基因上调表达，192个基因下调表达(图3A)。对照组与实验组共有表达的基因数为46 831个，对照组特有表达的基因数为10 580个，实验组特有表达的基因数为8 308个(图3B)。

Figure 3.  Differential gene analysis chart

GO功能注释与分类统计发现：差异表达基因被归类到2个功能分类组中，其中：分子功能中有17个基因涉及水解邻糖基化合物的水解酶活性；生物过程中细胞壁组织或生物合成过程注释到10个差异基因(图3C)。

将差异表达基因进行KEGG通路富集分析，对富集较多的前20个通路进行散点图绘制。差异基因富集最多的通路为苯丙烷生物合成；其次为淀粉和蔗糖代谢；再次为角质、木栓质和蜡的生物合成。而氰基氨基酸代谢、戊糖和葡萄糖醛酸相互转化、丙酮酸代谢、氨基糖和核苷酸糖代谢等通路富集较少(图3D)。

本研究有8个与苯丙烷生物合成相关的差异基因被富集到且呈上调表达，4个与过氧化物酶(POX)相关；2个与莽草酸-O-羟基肉桂酰转移酶(HCT)相关；2个与咖啡酸-O甲基转移酶(COMT)相关(图4)。本研究表明：GA₃处理会使得香榧胚性愈伤组织细胞增加对苯丙氨酸的合成与利用，调控HCT、COMT、POX活性的基因上调表达，促进了对香豆酰辅酶A参与合成咖啡酰莽草酸、愈创木酰木质素等产物的过程，同时生成酚类物质参与其他途径物质合成。该过程产生的酚类物质还可能参与清除部分愈伤组织和培养基产生的有害活性氧，提高香榧胚性愈伤组织的抗性^[17−19]。

Figure 4.  Pathways and differentially expressed genes related to phenylpropanoid biosynthesis and suberin synthesis

本研究有3个与木栓质生物合成相关的差异基因被富集到且呈上调表达，其中1个与脂肪酸ω-羟化酶(CYPB6B1)相关，2个与脂肪酰辅酶A还原酶(FAR)相关(图4)。与对照组相比，GA₃处理会使FAR与CYPB6B1上调表达，参与ω-羟基脂肪酸下游α,ω-二羧酸和α,ω-二醇合成途径，促进下游单酰基甘油酯的合成，共同参与形成木栓质单体或低聚物，同时产生酚类物质^[20]。该过程可结合苯丙烷合成途径产生的酚类物质，促进木栓质单体或低聚物的合成，应对环境变化并进行调整^[21−22]。

RT-qPCR验证结果(图5)表明：5个基因与转录组测序结果呈现相同的表达趋势，2个基因上调表达，3个基因下调表达，说明转录组测序结果可靠。

Figure 5.  RT-qPCR validation of candidate genes

植物细胞悬浮培养是一个持续变化的过程。在悬浮细胞培养过程中，接种量直接影响细胞的生长。细胞需要足够的空间和营养物质支持其增殖，选择适当的接种量对于细胞的生长至关重要。已有研究表明：落叶松胚性愈伤组织的最佳接种量为4 g·L^−1[6]。油茶Camellia oleifera愈伤细胞悬浮培养的最佳初始接种量为33 g· L^−1[23]。本研究中，胚性愈伤组织接种量为30 g·L⁻¹，有利于香榧悬浮细胞的生长，此时细胞增长率最大。在细胞生长过程中，细胞内部很多起关键调控作用的酶对pH有一定的要求。陈继光等^[24]研究表明：青钱柳在pH为5时细胞活力较强，细胞质量增量达最大，当pH达6时，细胞活力和质量增量都开始下降。本研究中，pH为5.7时，香榧胚性愈伤组织悬浮细胞抗逆性最强，细胞增长率最大。悬浮培养中摇床转速对细胞生长也有较大的影响。较低的转速会使愈伤组织堆积导致褐化^[4]；转速较高时会使细胞间产生碰撞导致破裂死亡，只有合适的转速才能使得细胞快速增殖生长^[25]。本研究表明：转速为110 r·min⁻¹时，香榧悬浮细胞生长最快。悬浮培养过程中对多种条件进行优化，能显著提升细胞的生长能力及抗性。

植物细胞在离体悬浮培养过程中，还须通过外施激素来促进细胞的生长。赤霉素在植物体的生命过程中起着重要的作用。当植物体处于胁迫条件时，外源赤霉素可以有效缓解胁迫的危害^[26]。有研究表明：加入外源GA₃可以有效缓解牡丹Paeonia × suffruticosa愈伤组织的褐化，并且提高愈伤组织诱导率^[14]。本研究表明：在差异基因的GO富集分析中，GA₃处理后上调表达基因主要参与催化活性、氧化还原过程、氧化还原酶活性等过程，说明活性氧的清除是香榧胚性愈伤组织增殖发育时的一个重要过程。在代谢通路中，差异表达基因主要富集在苯丙烷生物合成、淀粉和蔗糖代谢、木栓质合成等途径中。KEGG通路富集结果表明：苯丙烷生物合成途径中注释的差异表达基因最多，说明GA₃处理后，香榧胚性愈伤组织的次级代谢途径产生了一定的变化。

苯丙烷代谢途径是植物体内多种防御性次生代谢产物(如黄酮类、木质素及水杨酸等)的主要来源^[27]。有关木栓质相关合成途径的研究集中于植株根系以及部分受创伤部位^[28]，对于其在植物愈伤组织内的合成与作用研究较少。本研究中，木栓质合成途径也有7个差异基因被显著富集，且呈上调表达。已有研究表明：细胞色素P450单加氧酶(CYP)家族与FAR主要参与木栓质的合成^[29]。木栓质大分子的酚类组分已证明主要有阿魏酸、香豆酸等物质，通过苯丙烷生物合成途径产生，可继续生成木栓质单体中相应的底物物质，在一定程度上促进了木栓质单体的合成，影响细胞生长，以应对环境的胁迫^[30]。这与本研究结果相似，推测外源GA₃处理下香榧愈伤组织苯丙烷生物合成途径所产生的酚类物质不仅参与了活性氧的清除，还参与了木栓质单体的合成，对愈伤组织的生长及胁迫适应性产生一定的影响^[31]。

香榧胚性愈伤组织悬浮细胞生长曲线呈“S”型，接种后9~12 d为最佳继代时期，设置接种量为30 g·L⁻¹，pH为 5.7，摇床转速为110 r·min⁻¹，可获得生长良好的香榧胚性愈伤组织细胞，同时也可以提高悬浮细胞的抗氧化能力。通过转录组分析发现：外源GA₃处理下香榧胚性愈伤组织细胞苯丙烷生物合成、木栓质生物合成等相关代谢通路为差异基因的主要参与部分，且其中部分调空愈伤组织形成、生长和发育的关键基因表达量产生变化，可能在香榧愈伤组织的生长发育、适应环境变化以及应对胁迫相关过程中发挥重要作用。