植物激素作为植物体内广泛存在的一类痕量信号分子，对于调节植物的各种生长发育过程和环境应答具有非常重要的作用^[1]。赤霉素(Gibberellin，GA)是调控植物体生长发育必不可缺的一类植物激素，最早发现于水稻的“恶苗病”，患病水稻体内赤霉素过表达并且会出现疯长的现象^[2]。目前已经鉴定的赤霉素有136种，按照它们被发现的顺序依次命名为GA₁、GA₂、GA₃、GA₄等^[3-4]，但其中只有GA₁、GA₃、GA₄和GA₇具有生物活性，其余均为GAs分子的代谢物和中间产物^[5]。赤霉素在植物生长发育的各个阶段都发挥着重要作用，它参与调控植物生长发育的各种生理过程：如打破种子休眠并促进其萌发、植物生长、开花和果实发育等。其中最显著的作用是促进节间的伸长，从而促进植物的生长^[6-7]。赤霉素作用可以通过不同的途径实现，主要包括赤霉素的合成代谢和极性运输，以及赤霉素信号途径，包括受体感知和信号转导。

赤霉素的主要合成部位为植株的幼嫩组织，如根尖、茎尖以及生长中的果实和种子等^[8]。在高等植物中，根据赤霉素合成途径的亚细胞区域和该途径中所参与的酶的种类，赤霉素的生物合成分为3个阶段：1)在质体中，4个类异戊二烯形成一个C₂₀的线性分子，牦牛儿牦牛儿焦磷酸(Geranylgeranyl diphosphate，GGDP)环化成赤霉素的前体——内根-贝壳杉烯酸(Ent-kaurenoic acid)；2)在质体折叠区和内质网中，内根-贝壳杉烯酸逐渐转变为第一种形式的GA——GA₁₂；3)在细胞溶质中，GA₁₂先转变为C₂₀-GA，然后转变为C₁₉-GA，包括有活性的赤霉素。赤霉素的生物合成过程需要多种酶的参与，主要包括柯巴焦磷酸合酶(Ent-copalyl diphosphate synthase，CPS)、内根-贝壳杉烯合成酶(Ent-kaurene synthase，KS)、内根-贝壳杉烯氧化酶(Ent-kaurene oxidase，KO)、GA20-氧化酶(GA20-oxidase，GA20ox)、GA3-氧化酶(GA3-oxidase，GA3ox)和GA2-氧化酶(GA2-oxidase，GA2ox)等^[8-9]。尽管赤霉素代谢途径早期代谢通量较大，但赤霉素后期的代谢酶(GA20ox、GA3ox、GA2ox)对活性GA的精确调控具有关键作用^[10]。活性赤霉素前体的合成过程都需要GA20ox的参与；GA20ox催化合成具有生物活性的赤霉素的倒数第二步，可以将GA₁₂和GA₅₃ (C₂₀类赤霉素)转化为GA₉和GA₂₀ (C₁₉类赤霉素)^[11]。GA3ox是赤霉素合成途径中最后一步的关键酶，进而将GA₉和GA₂₀催化合成为具有生物活性的GA₁和GA₄^[12]。GA2ox是植物体内赤霉素降解过程的关键酶，它可以将具有生物活性的GA₁和GA₄分别转化为无活性的GA₈和GA₃₄，进而维持植物体内具有生物活性的GAs与其中间产物之间的动态平衡^[13]。GA20ox和GA3ox的缺少，会导致合成有生物活性赤霉素的量减少；而过表达的GA2ox则会降解植物体内生物活性的赤霉素^[14]。

植物体内GAs平衡的调控过程，信号通路同样发挥着重要的作用。赤霉素信号感知是植物体将细胞外赤霉素信号传入细胞内的过程。赤霉素GA信号可以被赤霉素受体GID1(Gibberellin insensitive dwarf1)感知到，当GID1感知到GA信号后，能够激活信号通路，从而调控下游基因的表达，进而影响植物的生长发育和形态建成。赤霉素信号转导的基因路径是GA-GID1-DELLA信号通路，组分包括赤霉素受体GID1、DELLA蛋白家族成员(DELLAs)以及与DELLA蛋白降解有关的带有F-Box基序的同源蛋白^[15]。GA的信号转导是通过降解DELLA蛋白来实现的。当GA处于低水平时, GID1不与GA结合, 使DELLA蛋白与赤霉素应答基因结合并抑制其活性，进而阻遏植物的生长；当赤霉素处于高水平时，GID1可感知GA信号并与之结合形成GA-GID1，然后结合到DELLA蛋白中形成GID1-GA-DELLA复合三聚体^[16]。当形成GA-GID1-DELLA复合体时，F-Box蛋白中的泛素连接酶SCF能结合到DELLA蛋白的GRAS区域，通过泛素蛋白体通道迅速降解DELLA蛋白，导致其阻遏作用被解除，植物体表现出正常的赤霉素反应^[15]。

近年来，在拟南芥Arabidopsis thaliana^[17]、水稻Oryza sativa^[18]、南瓜Cucurbita moschata^[19]、玉米Zea mays^[20]、大麦Hordeum vulgare^[21]、苹果Malus pumila^[22]、小桐子Jatropha curcas^[10]等多种植物中，赤霉素合成、降解与信号转导的基因家族已经被大量鉴定和克隆。而对于多年生木本植物毛竹Phyllostachys edulis，这方面的研究还甚少。毛竹作为世界上最重要的非木材林产品之一，以其茎秆的快速拔节生长而备受关注^[23]。目前毛竹高生长的机制尚未完全了解，许多研究者认为植物内源激素，特别是生长素和赤霉素，在调节毛竹节间细胞伸长方面发挥着重要的作用^[24]。在毛竹茎秆生长的过程中，赤霉素含量的动态变化与竹笋的快速生长有着密切的联系，毛竹笋生长初期，GA₃含量较低；而到了毛竹快速生长阶段，GA₃的含量急剧上升到了一个峰值^[24]。因此，了解毛竹中赤霉素合成代谢途径与信号转导途径关键基因可以从分子水平上更深刻地理解赤霉素与毛竹快速伸长的作用机制。同时，毛竹的基因组序列已经被发布^[25]，为进行赤霉素基因家族相关基因全面分析提供了一个绝佳的机会。本研究对赤霉素的合成代谢以及信号转导相关基因进行鉴定，比较了毛竹与水稻和拟南芥在赤霉素合成代谢与信号转导相关基因的进化关系；通过实时荧光定量PCR检测了这些基因在外施赤霉素的毛竹实生苗和毛竹笋不同部位的表达特征，为深入探讨赤霉素调控竹子生长发育的机制奠定基础。

1 材料与方法 1.1 毛竹赤霉素生物合成与代谢及信号转导相关基因的数据获取及进化分析 1.1.1 同源基因的数据获取

从毛竹基因组数据库(BambooGDB，http://www.bamboogdb.org/)下载毛竹的基因注释和基因组序列。为了鉴定毛竹赤霉素途径中涉及的基因家族，从NCBI (https://www.ncbi.nlm.nih.gov/)、拟南芥信息资源(TAIR，https://www.arabidopsis.org/)、水稻基因数据库(RiceData，http://www.ricedata.cn/gene/)和植物基因组数据库(http: //www.phytozome. net/)中下载相应GAs合成代谢和信号转导相关的氨基酸序列，作为BLAST软件搜索毛竹数据库的查询序列，搜索毛竹基因组数据库获得第一循环的候选基因；同时具有e≤10^–20的和score≥100的序列被用作第二循环搜索的新查询以避免丢失额外的候选基因。利用Pfam (http://pfam.sanger.ac.uk/)和SMART (http://smart.embl-heidelberg.de)两个在线网站中验证这些候选基因，E值设定为 < e^-15，其他参数为默认值；删除不含特征保守结构域的基因，通过基因注释最终确定目的基因。

1.1.2 同源基因的进化分析和基序分布

将水稻和拟南芥中已经鉴定的上述基因的蛋白序列与上述毛竹中最终确定的同源基因的蛋白序列利用ClustalX 1.81进行多重比对，参数设定为默认值，将生成的“aln.”文件导入MEGA 7.0^[26]，并通过该软件中的邻接法(Neighbor-Joining，NJ)构建系统进化树。采用以下参数：模式采用“Poisson model”，缺口设置“Complete deletion”，重复次数为Bootstrap = 1 000次。利用MEME在线软件(http://meme-suite.org/)进行蛋白保守基序的分析，E值设定为 < e^-20，保守结构域最大基序数量设置为10个，其余设置为默认参数。

1.2 毛竹种子与GA₃处理 1.2.1 GA₃处理种子的萌发率

试验所使用的毛竹种子是半同胞，2017年采自广西壮族自治区灵川县大境乡同一株开花毛竹，挑选大小一致的饱满种子作为试验材料。用70%酒精消毒种子表面60 s，灭菌水冲洗3遍；再用2% NaClO溶液浸泡10 min，灭菌水冲洗3遍，用不同浓度GA₃ (50 mg/L、100 mg/L、200 mg/L)溶液浸种，灭菌水浸种作为对照。每组选择100颗种子，黑暗浸泡24 h后用无菌水冲洗2–3次。置于灭菌的双层滤纸上萌发，萌发环境25 ℃并全程黑暗处理。每天检查并记录种子萌发数量，2周后统计种子发芽率，并进行差异显著性分析。

1.2.2 幼苗生长实验

为了确定GA₃浓度对毛竹生长的影响，将毛竹种子在0、50、100、200 mg/L的赤霉素溶液中浸泡24 h并黑暗萌发2周后，每组选择萌发状态良好且一致的幼芽20颗，转移到光照培养室，培养条件为温度25 ℃，光周期为光照16 h、黑暗8 h。4周后测量幼苗的根、茎、叶长度并进行差异显著性分析。根长为伸直的主根长度，茎长为根部以上到最近的第一片伸展叶片之间的距离，叶长为第二片伸展叶片的长度。

1.2.3 毛竹中赤霉素作用相关的基因对GA₃的响应

在25 ℃培养室中16 h光照/8 h黑暗处理的条件下，在土壤中培育毛竹幼苗4周。用不同浓度的GA₃ (0、50、100、200 mg/L)喷施叶片，以2 d的间隔喷施2周，在最后一次喷施后立即剪取毛竹幼苗的茎和叶组织用于下一步的RT-PCR分析；对于时间梯度的试验，4周龄幼苗用100 mg/L赤霉素喷施叶片，在不同的时间点(1、3、6、12、24 h)剪取毛竹幼苗的茎和叶组织，以4周龄喷施清水的毛竹幼苗作为对照^[27]。试验的所有材料均采用3次生物学重复，剪取的幼苗叶片立即置于液氮速冻，-80 ℃保存。

采用Trizol (TaKaRa公司)法提取毛竹幼苗的茎和叶组织的总RNA，用Prime Script^TM RT Master Mix (TaKaRa公司)去除基因组DNA和反转录合成cDNA。引物的设计选择在每个基因非保守区域，扩增片段均为150 bp左右，PCR检测引物特异性(表 1)。反应体系(10 μL)：5 μL SYBR^® Premix Ex Taq^TM Ⅱ (TaKaRa公司)，0.8 μL cDNA，0.4 μL primer-F，0.4 μL primer-R，3.4 μL无菌水。反应条件为：95 ℃，7 min；预扩增95 ℃ 10 s，58 ℃ 10 s，72 ℃ 15 s，共40个循环。选取NTB为内参基因^[28]，以喷施清水的材料作为对照组设定为1，处理样品中每个基因的表达水平均以此为参照，以3次生物学重复结果用2^-△△Ct的方法计算相对表达量^[29]，并进行差异显著性分析。

1.3 毛竹不同节间中赤霉素相关基因的表达分析

试验所采集的毛竹笋均来自于浙江农林大学亚热带森林培育国家重点实验室的翠竹园中，选取当年正处茎秆快速伸长的地径接近且株高150 cm左右的毛竹笋3株作为材料重复，以地面附近没有根的第一节定义为竹笋的第1节。每株笋按形态学位置自下而上分别取材第7节、第15节、第27节作为毛竹基部、中部和顶部的代表节间。RNA的提取和RT-PCR同1.2.3。采用均一化法比较每个基因在毛竹笋不同部位的表达量差异。将每个基因的表达量进行Z-score处理^[30]，使得每一个基因在所有样本的表达量都变成了均值为0、标准差为1的一组值。通过在线网站http://www.omicshare.com/tools/绘制表达量热图。

2 结果与分析 2.1 毛竹基因组中赤霉素作用关键基因家族的鉴定

为了揭示赤霉素影响毛竹茎秆生长的作用，本研究搜索了毛竹全基因组，鉴定毛竹中的赤霉素关键基因。利用已知的拟南芥中26个GAs相关基因和水稻中12个GAs相关基因作为模板，经过Blastp检索，基于E值设定 < e^-20和去除冗余序列，最终在毛竹基因组中共确定了23个候选基因。其中，17个赤霉素氧化酶类基因，包括8个GA20ox基因、1个GA3ox基因和8个GA2ox基因，按照它们与水稻同源蛋白的相似性高低分别命名为PhGA20ox1–8、PhGA3ox1和PhGA2ox1–8；2个赤霉素受体基因GID1，命名为PhGID1A和PhGID1B；2个正调控F-Box类基因，命名为PhGID2A和PhGID2B，以及2个负调控DELLA类基因，命名为PhSLR1和PhGAI (表 2)。

2.1.1 赤霉素的生物合成基因

GA20ox和GA3ox家族基因调控赤霉素生物合成，对植物细胞中的赤霉素积累起着关键作用。它们都是由小基因家族编码，在不同的物种中有着各自的保守基序。图 1为GA20ox蛋白家族在拟南芥、水稻和毛竹的进化树，显示各个基因在这3个物种中都比较保守，而毛竹与同为禾本科的水稻亲缘关系更为接近。其中PhGA20ox4、PhGA20ox6、PhGA20ox7、PhGA20ox8这4个蛋白序列的保守性比较低，并且MEME的基序搜索结果显示这4个基因的序列长度更短并且有保守基序的缺少现象。Motif 1、Motif 2和Motif 9在各个物种中高度保守，其中包含了与GA底物结合有关的LPWKET基序、H和D残基等(图 1，表 3)。

GA3ox催化合成具有生物活性的GA₁和GA₄，是活性赤霉素合成途径的关键一步。进化关系显示GA3ox在不同的物种中高度保守，毛竹的蛋白序列与拟南芥和水稻的同源蛋白的相似度分别达到了76.77%和81.92%，并且它们都具有共同的保守基序(图 2，表 4)。

2.1.2 赤霉素的降解基因

调控GA降解的关键基因家族GA2ox在拟南芥、水稻和毛竹中的保守性较差，但共同含有几个重要的保守基序，如Motif 1和Motif 6存在于所有的基因中(图 3，表 5)。拟南芥GA2ox家族的基因普遍存在着Motif 10缺少现象，而水稻的基因家族缺少了Motif 9；同样在毛竹的PhGA2ox5、PhGA2ox7、PhGA2ox8蛋白存在着许多片段的丢失，表明该组基因在不同的物种中可能存在着部分功能的异化。

GA20ox、GA3ox、GA2ox这3类赤霉素氧化酶都属于2OG-Fe (Ⅱ) oxygenase亚家族，蛋白序列都含有2OG-Fe保守结构域。研究赤霉素氧化酶基因家族在拟南芥、水稻和毛竹这3个物种中的进化属性发现，功能不同的GA20ox、GA3ox、GA2ox聚类分布在不同的亚家族中(图 4)。在3个物种中，GA20ox和GA2ox的基因数目均多于GA3ox，表明GA20ox和GA2ox经历了动态进化路线，导致更多功能的冗余。

2.1.3 赤霉素的信号转导基因

赤霉素信号通路的发现来源于对水稻赤霉素受体OsGID1的研究；相对于水稻，拟南芥中含有3种GID1的同源类似蛋白，分别是AtGID1A、AtGID1B和AtGID1C。通过蛋白序列同源比对，在毛竹基因组文库种找到2个GID1成员，命名为PhGID1A和PhGID1B，与水稻同源蛋白的相似性分别为89.34%和80.87%。系统进化树与保守基序的分析均显示GID1在不同的物种中高度保守(图 5，表 6)。

在毛竹基因组文库中找到了2个DELLA的成员，命名为PhSLR1和PhGAI。构建的系统演化树(图 6)显示毛竹的DELLA类基因与水稻更加近缘；并且在毛竹的2个DELLA蛋白都存在着保守序列的缺少现象，其中PhSLR1蛋白中缺少了Motif 1、5、8、10；而PhGAI蛋白中缺少了Motif 4、6、7、10，可能它们上功能上出现了异化和互补(表 7)。

F-Box蛋白在毛竹基因组文库中有了2个F-Box的成员，命名为PhGID2A和PhGID2B，其中PhGID2A与水稻同源性更高而PhGID2B与拟南芥更加近缘。F-Box家族蛋白序列长度短且保守基序数目少，但是在不同的物种中具有着高度的保守性(图 7，表 8)。

2.2 毛竹幼苗对外源GA₃的响应 2.2.1 生长参数统计

将记录的种子萌发数据进行统计，分别计算不同浓度处理后的种子萌发率，结果显示GA₃处理显著提高了毛竹种子的萌发率(图 8A)。100 mg/L处理组的萌发率最高，达到了67%，50 mg/L处理组和200 mg/L处理组的萌发率分别是59%和54%，均高于对照组的43%，说明在一定范围内不同浓度的赤霉素浸泡毛竹种子均能提高其种子萌发率且种子萌发的最佳赤霉素浓度在100 mg/L左右。200 mg/L处理组的萌发率虽高于对照组但已呈现下降趋势，说明过高浓度的赤霉素反而会造成生长逆境，抑制种子的萌发。

同时，用不同浓度的GA₃浸泡毛竹种子，萌发后的毛竹幼苗的根茎叶长度都有所增加(图 8B–D。茎长和叶长增加效果显著，3个处理组的主茎长度分别提高了44.0%、55.1%、50.1%；主叶长度分别提高了41.4%、62.9%、55.9%。赤霉素浓度在100 mg/L时茎与叶达到最佳的生长效果，200 mg/L赤霉素处理虽然促进效果仍显著，但是相比100 mg/L已经有了一定程度的抑制。而根系的增长效果不明显，3个处理组仅分别增加了1.8%、4.8%、5.4%。表明外源赤霉素可以有效刺激毛竹的营养生长，明显促进了茎和叶的生长，但对根的生长基本无促进作用。另外，赤霉素处理的毛竹幼苗，“疯长”现象仅体现在垂直方向上的高度，茎秆细长没有粗度上的增加，并且大量出现倒伏现象。

2.2.2 外源GA3处理下赤霉素作用相关基因的表达变化

通过qRT-PCR检测不同浓度(0、50、100、200 mg/L)赤霉素GA₃处理下PhGA20ox1、PhGA3ox1、PhGA2ox1、PhGID1A、PhGID2A和PhSLR1的相对表达量(图 9)。

应用外源性赤霉素GA₃处理后，PhGA20ox1的表达量显著下调，且外源赤霉素浓度越高其表达水平越低。相一致的是，参与赤霉素生物合成的基因PhGA3ox1在各个浓度处理组中也都显示下调。而PhGA2ox1在GA₃处理后的幼苗中表达量升高，并且在高浓度处理组的表达量更高。至于信号转导的相关基因，PhSLR1的相对表达量受到了GA₃严格的负反馈，各个浓度的处理组与对照组相比全部下调表达；而GA₃处理却使PhGID2A表现为极显著的上调表达，200 mg/L处理组的PhGID2A表达量达到了对照组的2.7倍。赤霉素受体PhGID1A对外源GA₃的浓度变化不敏感，在各个处理组的表达量都略高于对照组。

对4周龄毛竹幼苗喷施100 mg/L GA₃，通过qRT-PCR检测其在处理后不同时期(0、1、3、6、12、24 h) PhGA20ox1、PhGA3ox1、PhGA2ox1、PhGID1A、PhGID2A和PhSLR1的相对表达量(图 10)。结果显示这些基因在GA₃处理后的24 h内表达量都发现了显著的变化。赤霉素生物合成相关基因PhGA20ox1和PhGA3ox1整体呈表达量下调趋势。其中PhGA20ox1在处理后3 h表达量达到最低，仅为对照的43.8%，之后开始升高但仍低于未处理组，处理后24 h的表达量为未处理时的80.1%；而PhGA3ox1对GA₃的响应有一个更为复杂的过程，其表达量也是在3 h时达到最低，只有对照组的32.9%，而6 h时恢复到了未处理时的表达水平，随后在12 h和24 h又出现了明显的降低和升高，最终在24 h时的表达量开始相比未处理时提高了38.9%。赤霉素代谢基因PhGA20ox1在GA₃处理后表达量开始显著上调，1 h时的表达量比对照提高了67.0%，随后开始持续下降，12 h后的表达量下降到了未处理时的水平，24 h后表达量为对照的111.5%。赤霉素受体基因PhGID1A在外源赤霉素处理过程中变化相对温和，在最初的1 h和3 h，表达量比对照下降了18.2%和30.8%，6 h时的表达量升高到了未处理时的111.6%，12 h开始略微下降，24 h时又逐渐升高到了未处理时的125.7%。PhGID2A在处理后的各个时间段均表达上调，1 h时表达量上调10.9%，3 h时已急剧上升为对照的237.8%，达到最高，随后逐渐下降但所有处理时间段表达量均高于未处理组，24 h表达量降至对照的121.6%。DELLA类基因PhSLR1表达量随时间不断呈现下调-上调-下调的趋势，然而在每个时间段的表达量均未高于对照组，最低的表达量出现在6 h，为对照组的45.4%。

2.3 毛竹笋不同节间的表达分析

赤霉素关键基因在1.5 m高毛竹笋3个节间均检测到表达，并且各个基因的表达趋势在3株重复中具有较好的一致性，推测这些基因在毛竹快速生长的过程都发挥着一定的作用(图 11)。赤霉素生物合成相关的基因家族PhGA20ox和PhGA3ox在毛竹笋的不同部位都高表达，尤其是在顶部和中部；而调控赤霉素失活的PhGA2ox基因家族也在毛竹笋的形态学上端高表达。赤霉素受体PhGID1表达主要集中在顶部和中部，基部基本上没有检测到表达。PhGID2在3个部位也都有表达，其中顶部的表达量最高，基部表达量最低，但不同部位的表达差异不是很明显；而DELLA相关基因在根部的表达量显著高于顶部及中部。可以看出，PhGID1与PhGID2的表达具有一致性，而与DELLA的表达趋势相反。

3 讨论

毛竹是我国分布面积最广且经济价值最高的竹种，毛竹笋的生长速度很快，高生长时期每天净高度增长量高达50 cm以上，这种快速生长的机理具有重要的学术研究价值。前人的大量研究发现，毛竹的高生长过程伴随着赤霉素等植物激素含量的急剧增加^[24,31]，因此推测赤霉素在毛竹的快速生长中发挥着重要的作用。在竹笋快速生长过程中，对赤霉素动态变化及其关键基因表达量变化的研究有助于我们了解毛竹茎秆快速伸长的机制。

通过生物信息学分析，从毛竹基因组中共鉴定出了17个赤霉素合成代谢基因和6个信号转导的相关基因。相比于拟南芥、水稻、玉米等植物，毛竹基因组更为庞大，但是参与赤霉素生物合成和信号转导的基因成员数量却较少，尤其是GA₃氧化酶家族，仅检索到了1位成员。这可能是由于本次研究删去了大量序列相似性较低的同源基因以及冗余基因，并且当前的毛竹基因组存在着个别序列的注释信息不完全和CDS缺失等组装问题。另一种假设是由于物种分化后，草本植物中竹子在蛋白质直系同源簇数据库(COGs)的构建和扩增中种类数量最多^[25]，因此导致毛竹中的基因家族比其他物种更加分散和极化。拟南芥、水稻和毛竹的赤霉素相关基因的系统发育分析显示赤霉素的合成代谢与信号传导途径基因在这些物种中是高度保守的。而毛竹的大多数基因与水稻的同源基因相似性更高，可能毛竹与水稻的赤霉素关键基因在进化关系上有着共同的祖先^[32]，而序列的高保守性预示着功能的保守性。

赤霉素参与调控植物生理发育的各个过程，包括种子萌发、根的生长、茎的伸长和叶片伸展等。王巍等^[33]用不同浓度的GA₃处理核桃种子，发现300 mg/L的GA₃处理组发芽率最高达到了87.2%。高春智等^[34]在樟子松上的研究也表明，喷施GA₃能显著提高樟子松种子的发芽率。毛竹自然结实率低且种子采集困难，室温贮藏还会使种子丧失发芽能力，这些都在一定程度上限制了毛竹实生苗的培育工作。本研究用不同浓度的GA₃处理毛竹种子的实验也同样表明，用不同浓度的GA₃浸种后促进了毛竹种子的萌发并且100 mg/L GA₃的促进效果最显著，并且有着最佳的浓度范围。吴建明等^[35]研究发现赤霉素处理可以诱导甘蔗节间伸长，增加株高；杨益善等^[36]在水稻上的研究也发现，赤霉素可以促进水稻茎的伸长。我们同时比较了不同浓度GA₃浸种处理对毛竹幼苗根、茎、叶生长的影响，结果与GA₃对发芽的效应一致，赤霉素处理后的毛竹茎和叶的伸长促进效果明显，并且最适浓度在100 mg/L左右，而对根的影响效果不明显。总之，赤霉素处理对毛竹种子的萌发与幼苗的生长均有明显的促进作用。

参与毛竹生物合成的相关基因对于外源GA₃的响应有着协调一致性，外施GA₃使毛竹幼苗内的赤霉素生物合成基因PhGA20ox1和PhGA3ox1表达量明显下调，并且浓度越高对基因表达的抑制作用越明显；相反，催化植物中活性GAs灭活的基因PhGA2ox1在高浓度的处理组表达量更高。同时对GA₃处理后的5个时期的表达量进行测定，发现PhGA20ox1表达量先下调后上调，整体表现下调趋势，而PhGA2ox1的表达量整体呈上调趋势。这与前人关于赤霉素处理植物的结论相同，Carrera等^[37]研究表明，马铃薯Solanum tuberosum L.的gal矮化突变体具有较高的GA-20氧化酶转录水平，外施GA₃后，3种GA-20氧化酶(StGA20ox1、StGA20ox2和StGA20ox3)在马铃薯gal突变体内表达明显减少；董凤等^[22]对苹果的研究也发现外施赤霉素会使GA20ox表达量减少而GA2ox表达增高。与PhGA20ox1的表达趋势不同，PhGA3ox1的表达量在处理后的某段时期还出现了升高的现象，可能是因为GA3ox作为赤霉素生物合成的最后一步，在赤霉素含量骤然升高或降低时，GA3ox起着更重要的反调节作用，从而维持植物体内赤霉素的平衡^[38]。并且，在外界GA₃的刺激下，大量赤霉素受体基因与赤霉素相结合，导致赤霉素受体基因的表达量升高，结果与李俊等^[39]东方山羊豆在GA₃处理下的基因表达水平相一致。总之，当外施GA₃时，一方面毛竹幼苗通过抑制体内赤霉素合成反应并加强体内活性赤霉素灭活基因GA2ox的表达，来实现植物体内GA水平的稳态平衡；另一方面，赤霉素受体基因GID1及正调控因子基因F-Box的表达水平显著提高，解除了DELLA蛋白的积累对毛竹生长的阻遏作用，植株表现出正常的赤霉素反应。

本研究选取了快速生长时期毛竹笋的基部、中部和顶部节间进行qRT-PCR分析，通过组织特异性分析发现，在毛竹笋快速生长阶段的3个节间，赤霉素合成代谢与信号转导途径关键酶基因的表达量均发生了一定变化。生物合成相关基因家族PhGA20ox和PhGA3ox在毛竹笋的形态学上部都过量表达，为顶端节间的伸长提供了充足的内源激素；而调控赤霉素代谢的基因家族PhGA2ox在中部和顶部也高表达，毛竹笋幼嫩部位的生命活动最旺盛，大量合成赤霉素，PhGA2ox的高表达能够将过量合成的赤霉素进行降解^[13]。丁兴萃^[31]曾分析了毛竹笋各部位的内源激素含量及动态变化，发现在竹笋生长初期，赤霉素在笋尖合成，通过极性运输到中部，一般赤霉素在各个部位的含量为顶部 > 中部 > 基部。而本研究从合成代谢关键基因进行分析，发现基因表达模式和赤霉素体内含量分布一致。信号转导途径的基因PhGID1和PhGID2在笋的各个部位也都有表达，并且在笋顶部和中部表达要明显高于基部。而DELLA类蛋白相关基因表达趋势相反，在毛竹笋形态学上部的中部和顶部表达量非常低甚至没有检测到表达。这三类基因已在其他物种中广泛报道，DELLA的大量表达抑制植物的生长伸长，是赤霉素的负调控基因；GID1是赤霉素的受体，能够感知赤霉素信号并介导DELLA类蛋白的泛素化降解；GID2是降解DELLA蛋白的泛素酶体复合物的主要成分，DELLA的降解使赤霉素得以大量表达并诱导植物快速生长^[15,40]。毛竹笋中赤霉素信号转导途径相关基因的表达模式与其生物学功能也是相吻合的：毛竹笋进入快速生长阶段，笋尖大量合成赤霉素，PhGID1大量表达并接收赤霉素信号，然后与DELLA蛋白形成GID1-GA-DELLA复合体，阻止赤霉素信号的传递；此时大量表达的PhGID2将DELLA蛋白泛素化标记并介导其降解，DELLA的阻遏作用解除，赤霉素大量表达，笋开始快速生长。

4 结论

本研究对赤霉素生物合成与信号转导途径的相关基因家族在毛竹中的同源基因进行了系统的鉴定，找到了6个多基因家族编码的蛋白，包括参与赤霉素合成代谢的GA20ox、GA3ox和GA2ox以及信号转导途径的GID1、GID2和DELLA蛋白。并初步探索了外源赤霉素处理对毛竹种子萌发的影响，证明了赤霉素在毛竹中也有着促进种子萌发和促进茎秆伸长等作用，且有最佳的浓度范围；同时，外施赤霉素会引起毛竹体内赤霉素相关基因表达模式变化，植物体通过反馈调节的方式控制体内活性GA的平衡。另外，本研究还发现了大量的赤霉素类基因参与到毛竹的快速生长这一复杂的生理过程中，并且基因的表达模式与毛竹体内的赤霉素含量分布一致，也与其生物学功能一致，为深入探讨赤霉素调控竹子生长发育的机制奠定了基础。但是，本研究尚未进行基因克隆等方面的分析，在后续的试验中应当重点关注。