王江英 朱朋波 湯雪燕 孫明偉 葛金濤 趙統(tǒng)利 李紀(jì)元 惠林沖 邵小斌

摘要：對(duì)0 mg/L、800 mg/L赤霉素處理的矮生山茶恨天高葉片、莖段和莖尖進(jìn)行了Illumina HiSeq高通量轉(zhuǎn)錄組測(cè)序，共獲得476 635個(gè)Unigenes，其中181 115（38.00%）個(gè)Unigenes獲得了注釋;差異基因分析發(fā)現(xiàn)了43 624個(gè)差異基因，其中包括43 295個(gè)特異的和329個(gè)共同的差異基因;GO和KEGG富集分析結(jié)果顯示，恨天高葉片增大、莖段伸長(zhǎng)及莖尖生長(zhǎng)與糖代謝、苯丙烷生物合成代謝相關(guān);莖尖生長(zhǎng)還與植物激素信號(hào)傳導(dǎo)和二萜生物合成有關(guān)。通過(guò)赤霉素處理后恨天高葉片、莖段、莖尖轉(zhuǎn)錄組分析發(fā)現(xiàn)，莖尖組織中參與赤霉素生物合成的KAO、GA20ox、GA3ox和GA2ox基因表達(dá)量發(fā)生變化，赤霉素信號(hào)傳導(dǎo)過(guò)程中DELLA蛋白調(diào)控以及赤霉素和生長(zhǎng)素、細(xì)胞分裂素等其他激素間的互作促進(jìn)恨天高植株生長(zhǎng)。

關(guān)鍵詞：矮生恨天高;赤霉素;轉(zhuǎn)錄組;差異表達(dá)基因

中圖分類(lèi)號(hào)：S685.14文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-4440（2020）01-0047-10

Abstract：Illumina HiSeq high-throughput sequencing was performed on the leaves， stem segments and shoot tips of Camellia reticulata Hentiangao treated with ?mg/L0 and 800 mg/L gibberellin， respectively. A total of 476 635 unigenes were obtained， and 181 115 （38.00%） of them were annotated. In addition， 43 624 differentially expressed genes （DEGs）， including 43 295 unique and 329 common， were identified. GO and KEGG enrichment analysis results showed that leaf enlargement， stem elongation and shoot tip growth were related to glucose metabdism and phenylpropanoid biosynthesis. Furthermore， many DEGs in shoot tip growth were involved in plant hormone signal transduction and diterpenoid biosynthesis. According to the transcriptome analysis， three factors were found to promote the growth of Hentiangao. The first was the expression changes of KAO， GA20ox， GA3ox and GA2ox genes involved in gibberellin biosynthesis in shoot tips. The second was the regulation of DELLA protein in gibberellin signal transduction. The third was the interaction of gibberellin with auxin， cytokinin and other hormones.

Key words：Camellia reticulata Hentiangao;gibberellin;transcriptome;differentially expressed genes

赤霉素（Gibberellins，GAs），二萜類(lèi)植物激素，主要促進(jìn)植株的營(yíng)養(yǎng)生長(zhǎng)，包括節(jié)間的延長(zhǎng)、根的伸長(zhǎng)和葉的擴(kuò)展[1-3]。在赤霉素合成途徑中，萜烯合成酶TPSs、細(xì)胞色素P450s和GA氧化酶（GAoxs）將前體牻牛兒牻牛兒焦磷酸GGDP催化成活性赤霉素，其中貝殼杉烯酸氧化酶（KAO）屬于細(xì)胞色素P450，催化赤霉素合成前期階段。GA氧化酶包括GA2ox、GA20ox和GA3ox，是赤霉素合成和代謝后期階段的主要調(diào)控點(diǎn)，將非活性赤霉素催化成活性赤霉素[4-7]。前期研究發(fā)現(xiàn)，多種植物的KAO突變能誘導(dǎo)植株矮化，如水稻d35突變體[8]。梨矮化砧木中KAO表達(dá)量低于喬化和半矮化的品種[9]。相對(duì)于早期合成酶，后期生物合成酶易受植物生長(zhǎng)發(fā)育和環(huán)境因素調(diào)節(jié)，這些酶在生物活性赤霉素水平的調(diào)控中起著關(guān)鍵作用。GA20ox和GA3ox過(guò)量表達(dá)后提高活性赤霉素的含量，加快植物生長(zhǎng)，節(jié)間伸長(zhǎng)，GA20ox和GA3ox功能缺失會(huì)引起植株矮化[10-11]。GA2ox的作用則是降低活性赤霉素的水平[11-12]。

赤霉素信號(hào)傳導(dǎo)途徑中，DELLA蛋白作為植物生長(zhǎng)抑制劑可被活性赤霉素抵消[13]，赤霉素信號(hào)啟動(dòng)DELLA蛋白的降解，調(diào)控赤霉素反應(yīng)基因的表達(dá)，從而促進(jìn)植物的生長(zhǎng)，所以DELLA蛋白對(duì)植物的生長(zhǎng)發(fā)育起負(fù)反饋調(diào)節(jié)[14-16]。目前DELLA基因在擬南芥[17]、葡萄[18]、草莓[19]和甘藍(lán)型油菜[20]等植物中陸續(xù)被克隆及分析，為植物矮化性狀研究提供依據(jù)。

恨天高（Camellia reticulata Hentiangao）是云南山茶最珍貴的品種之一，其株形矮小，生長(zhǎng)緩慢，花瓣平展，完全重瓣，形姿優(yōu)美，葉濃花艷，具有極高的觀賞價(jià)值[21]。隨著城市的快速發(fā)展，矮生和盆栽木本花卉植物在園林綠化、居室美化等方面有著重要作用。通過(guò)轉(zhuǎn)錄組分析挖掘高效調(diào)控恨天高株高的關(guān)鍵基因，為培育微型山茶新品種奠定基礎(chǔ)。

本研究以云南山茶矮生品種恨天高為材料，利用轉(zhuǎn)錄組測(cè)序分析恨天高植株矮生與赤霉素合成代謝和信號(hào)傳導(dǎo)間的關(guān)系，以期獲得調(diào)控植物株高的重要基因，并揭示其矮生機(jī)制。

1材料與方法

1.1試驗(yàn)材料

以立地條件相同且生長(zhǎng)年齡相近的恨天高山茶（Camellia reticulate Hentiangao）為材料，利用不同質(zhì)量濃度梯度[0 mg/L（CK）、200 mg/L（T1）、400 mg/L（T2）、800 mg/L（T3）和1 000 mg/L（T4）]的赤霉素（GA3）進(jìn)行莖尖涂抹。75 d后采集各質(zhì)量濃度處理后的葉片、莖段及莖尖組織進(jìn)行液氮速凍，-80 ℃保存?zhèn)溆谩?/p>

1.2莖尖和葉片總RNA提取

參照北京艾德萊生物科技有限公司RN53-EASYspin Plus植物試劑盒說(shuō)明書(shū)提取樣品總RNA。利用1%的瓊脂糖凝膠電泳和Nanodrop 2000分光光度計(jì)檢測(cè)RNA完整度、濃度、A260/280、A260/230及核酸峰度是否正常，之后再經(jīng)Agilent bioanalyzer 2100進(jìn)行RNA質(zhì)量檢測(cè)。

1.3赤霉素誘導(dǎo)恨天高植株生長(zhǎng)的質(zhì)量濃度篩選

分別利用0 mg/L、200 mg/L、400 mg/L、800 mg/L和1 000 mg/L的GA3涂抹5年生恨天高植株莖尖，每2 d涂1次，每個(gè)質(zhì)量濃度處理3株，每株分別處理5個(gè)莖尖，30 d和75 d后分別測(cè)量莖尖生長(zhǎng)值。

1.4文庫(kù)構(gòu)建和轉(zhuǎn)錄組測(cè)序

恨天高葉片、莖段和莖尖總RNA質(zhì)量分別檢測(cè)合格后，利用NEBNext Poly（A） mRNA Magnetic Isolation Module進(jìn)行mRNA捕獲，利用NEBNext Ultra RNA Library Prep Kit for Illumina和Beckman Agencourt AMPure XP beads分別進(jìn)行文庫(kù)構(gòu)建和文庫(kù)純化，利用TruSeq PE Cluster Kit V4完成cBOT自動(dòng)成簇，最后通過(guò)TrSeq SBS Kit v4-HS進(jìn)行HiSeq測(cè)序。以上文庫(kù)構(gòu)建和轉(zhuǎn)錄組測(cè)序均由蘇州金唯智生物科技有限公司完成。

1.5生物信息學(xué)分析

對(duì)測(cè)序結(jié)果原始圖像數(shù)據(jù)利用軟件Bcl2fastq（v2.17.1.14）進(jìn)行圖像堿基識(shí)別，初步質(zhì)量分析后得到原始測(cè)序數(shù)據(jù)。測(cè)序數(shù)據(jù)質(zhì)量評(píng)估采用FastQC（v0.10.1）進(jìn)行分析，使用Cutadapt（version1.9.1）對(duì)測(cè)序原始數(shù)據(jù)去除接頭以及低質(zhì)量序列等，得到后續(xù)信息分析用的高質(zhì)量有效序列。采用Trinity對(duì)樣品數(shù)據(jù)從頭組裝，組裝結(jié)果通過(guò)序列聚類(lèi)進(jìn)行進(jìn)一步序列拼接和去冗余處理，得到長(zhǎng)的非冗余Unigenes，利用TransDecoder軟件進(jìn)行ORF分析[22]。

Unigenes表達(dá)計(jì)算使用RSEM軟件，該軟件使用FPKM（Fragment per kilo bases per million reads）方法計(jì)算基因表達(dá)量。利用EdgeR軟件進(jìn)行基因差異表達(dá)分析，以變化2倍以上、Pvalue≤0.05且錯(cuò)誤發(fā)現(xiàn)率（FDR）<0.01為差異表達(dá)基因篩選標(biāo)準(zhǔn)。最后基于差異表達(dá)基因進(jìn)行GO功能注釋、KEGG代謝通路分析和COG功能注釋。

1.6實(shí)時(shí)熒光定量PCR驗(yàn)證

利用實(shí)時(shí)熒光定量PCR驗(yàn)證轉(zhuǎn)錄組測(cè)序的準(zhǔn)確性，挑選5個(gè)差異表達(dá)基因設(shè)計(jì)引物，以山茶18S RNA為內(nèi)標(biāo)基因（所用引物見(jiàn)表1），生物學(xué)試驗(yàn)重復(fù)3次。運(yùn)用2-△△Ct計(jì)算方法求得差異基因的表達(dá)量。

2結(jié)果與分析

2.1外施赤霉素促進(jìn)恨天高莖尖伸長(zhǎng)

選擇立地條件相同、生長(zhǎng)狀態(tài)良好且一致的三年生恨天高（Camellia reticulate Hentiangao）幼苗為材料（圖1A），利用不同質(zhì)量濃度梯度的GA3進(jìn)行莖尖涂抹（圖1B）。研究結(jié)果發(fā)現(xiàn)，GA3處理30 d后，對(duì)照組和試驗(yàn)組莖尖長(zhǎng)度差異顯著（P<0.05）（圖1C），其中CK的莖尖長(zhǎng)度最短，T1、T2、T4處理分別是CK的1.23倍、1.58倍和1.92倍，T3為CK的2.53倍（圖2A）;GA3處理75 d后，對(duì)照組和試驗(yàn)組莖尖長(zhǎng)度差異顯著（P<0.05）（圖1D），其中CK莖尖長(zhǎng)度依舊最短，其次是T1處理和T2處理，為CK的1.59倍和1.93倍，而T4處理和T3處理分別為CK的2.27倍和3.61倍（圖2B）。比較發(fā)現(xiàn)，對(duì)照組和試驗(yàn)組在2次測(cè)量的45 d間隔期，CK的莖尖少許伸長(zhǎng)，試驗(yàn)組除了T1處理長(zhǎng)勢(shì)變化不明顯外，T2、T3、T4處理長(zhǎng)勢(shì)旺，尤其T3處理（圖2C）。由此可以推測(cè)，外源GA3能誘導(dǎo)恨天高山茶的植株生長(zhǎng)，并且隨著外施質(zhì)量濃度的增高（0 mg/L、200 mg/L、400 mg/L），莖尖生長(zhǎng)值也隨之增加，當(dāng)質(zhì)量濃度達(dá)800 mg/L時(shí)，莖尖生長(zhǎng)最快，然而質(zhì)量濃度達(dá)1 000 mg/L時(shí)，莖尖增加幅度卻有所下降，從而表明800 mg/L是誘導(dǎo)恨天高山茶植株生長(zhǎng)的較佳質(zhì)量濃度。

2.2轉(zhuǎn)錄組序列組裝及Unigene注釋分析

利用最佳赤霉素誘導(dǎo)質(zhì)量濃度800 mg/L處理恨天高莖尖組織，0 mg/L的GA3處理作對(duì)照，75 d后采集葉片、莖段、莖尖進(jìn)行轉(zhuǎn)錄組測(cè)序，其中葉片和莖段分別取2組樣品進(jìn)行測(cè)試，莖尖只取1組樣品進(jìn)行測(cè)試。

使用軟件Cutadapt對(duì)測(cè)序原始數(shù)據(jù)去除接頭以及低質(zhì)量序列，獲得后續(xù)信息分析用的過(guò)濾數(shù)據(jù)。恨天高3種組織及對(duì)照共8個(gè)文庫(kù)產(chǎn)生了411 464 870條序列，約60.96 Gb轉(zhuǎn)錄組數(shù)據(jù)，Unigene平均長(zhǎng)度為147.38～148.61 bp，G+C堿基含量為44.58%～46.84%（表2）。由表3可以看出，轉(zhuǎn)錄組數(shù)據(jù)經(jīng)組裝產(chǎn)生了43 073 816個(gè)重疊群，476 635個(gè)Unigenes，平均長(zhǎng)度為497.69 bp。在所有Unigenes序列中，序列長(zhǎng)度小于500 bp的占71.60%。Unigene注釋分析顯示有181 115條Unigenes獲得注釋?zhuān)糢nigenes總數(shù)的38.00%，其中獲得注釋較多的數(shù)據(jù)庫(kù)分別是Nr注釋172 121條（36.11%），COG注釋74 120條（15.55%），Swissport注釋106 818條（22.41%）和KEGG注釋25 837條（5.42%）（圖3）。

2.3差異表達(dá)基因篩選及分析

將Clean data比對(duì)到Unigene上，利用bowtie2軟件進(jìn)行短reads比對(duì)，結(jié)果顯示，唯一比對(duì)到參考基因的reads數(shù)量在2.5×107～3.6×107，用于基因表達(dá)分析。在葉片、莖段和莖尖3種組織中，發(fā)現(xiàn)了特異的和共同的差異表達(dá)基因，共發(fā)現(xiàn)43 624個(gè)差異表達(dá)基因，包括43 295個(gè)特異差異表達(dá)基因和329個(gè)共同差異表達(dá)基因（圖4A）。

CK1（對(duì)照葉片樣品）與L1（800 mg/L外源赤霉素處理后葉片樣品1）比較中發(fā)現(xiàn)27 424個(gè)差異表達(dá)基因，其中上調(diào)表達(dá)基因16 495個(gè)，下調(diào)表達(dá)基因10 929個(gè);CK1與L2（800 mg/L外源赤霉素處理后葉片樣品2）比較中發(fā)現(xiàn)23 450個(gè)差異表達(dá)基因，其中上調(diào)表達(dá)基因12 956個(gè)，下調(diào)表達(dá)基因10 494個(gè);CK2（對(duì)照莖段樣品）與SS1（800 mg/L外源赤霉素處理后莖段樣品1）比較中發(fā)現(xiàn)20 773個(gè)差異表達(dá)基因，其中上調(diào)表達(dá)基因11 168個(gè)，下調(diào)表達(dá)基因9 605個(gè);CK2與SS2（800 mg/L外源赤霉素處理后莖段樣品2）比較中發(fā)現(xiàn)21 040個(gè)差異表達(dá)基因，其中上調(diào)表達(dá)基因11 706個(gè)，下調(diào)表達(dá)基因9 334個(gè);CK3（對(duì)照莖尖樣品）與ST1（800 mg/L外源赤霉素處理后莖尖樣品）比較中發(fā)現(xiàn)4 256個(gè)差異表達(dá)基因，其中上調(diào)表達(dá)基因1 734個(gè)，下調(diào)表達(dá)基因2 522個(gè)（圖4B）。在5組比較中，CK3與ST1比較獲得的上調(diào)和下調(diào)基因數(shù)量均最低，可能是由于莖尖組織相對(duì)于葉片和莖段而言，與植株伸長(zhǎng)生長(zhǎng)更為密切，獲得的差異表達(dá)基因控制植株生長(zhǎng)的作用更大。

根據(jù)FPKM值將8個(gè)文庫(kù)中的基因表達(dá)水平分為極高、高、中、低和極低5類(lèi)（圖4C），其中差異表達(dá)基因中占比例最大的是極低表達(dá)基因，F(xiàn)PKM為0～1，占比例最小的是極高表達(dá)基因，F(xiàn)PKM>60。另外，在葉片、莖段和莖尖等組織中檢測(cè)到共同差異表達(dá)基因329個(gè)，通過(guò)層次聚類(lèi)分析將其劃分成A、B 2個(gè)簇（圖4D）。

2.4差異基因GO富集分析

通過(guò)GO富集分析對(duì)差異表達(dá)基因進(jìn)行生物學(xué)功能研究，在錯(cuò)誤發(fā)現(xiàn)率FDR<0.01時(shí)，43個(gè)GO條目富集在CK1/L1，42個(gè)GO條目富集在CK1/L2，39個(gè)GO條目富集在CK2/SS1，39個(gè)GO條目富集在CK2/SS2，37個(gè)GO條目富集在CK3/ST1?；谒蠫O富集分析，確定37個(gè)共同GO類(lèi)別和8個(gè)差異GO類(lèi)別，在37個(gè)共同GO條目中，催化活性、結(jié)合、代謝過(guò)程、細(xì)胞過(guò)程、細(xì)胞組分和細(xì)胞器是富集基因最多的GO類(lèi)別。在8個(gè)差異GO類(lèi)別中，CK3/ST1無(wú)分子功能調(diào)節(jié)劑富集;鳥(niǎo)苷酸交換因子活性顯著富集在CK1/L1和CK1/L2;CK2/SS1和CK2/SS2無(wú)細(xì)胞外基質(zhì)部分和細(xì)胞外區(qū)域部分富集;CK3/ST1無(wú)類(lèi)核和細(xì)胞外基質(zhì)富集;轉(zhuǎn)運(yùn)狀態(tài)顯著富集在CK1/L1;CK3/ST1無(wú)生物附著富集。由此可以看出，在外源赤霉素處理下葉片、莖段和莖尖組織中的GO類(lèi)別存在差異，誘導(dǎo)植株生長(zhǎng)。

2.5差異基因KEGG富集分析

KEGG富集分析結(jié)果顯示，分別有30、30、18、14和27條代謝途徑富集在CK1/L1、CK1/L2、CK2/SS1、CK2/SS2和CK3/ST1，顯著富集的則為5、4、5、4和5條代謝途徑（圖5）。苯丙烷生物合成、戊糖和葡萄糖醛酸相互轉(zhuǎn)化過(guò)程是共同顯著富集的代謝途徑（圖5），倍半萜和三萜生物合成富集在CK1/L1、CK1/L2和CK2/SS2（圖5A，圖5B，圖5D），細(xì)胞色素P450對(duì)外源化合物代謝的影響、藥物代謝-細(xì)胞色素P450只富集在CK1/L1（圖5A），糖酵解/糖異生只富集在CK1/L2（圖5B），光合-天線(xiàn)蛋白、油菜素內(nèi)酯生物合成只富集在CK2/SS1（圖5C），玉米素生物合成只富集在CK2/SS1和CK2/SS2（圖5C、圖5D），CK3/ST1顯著富集植物激素信號(hào)轉(zhuǎn)導(dǎo)、二萜生物合成和類(lèi)胡蘿卜素生物合成（圖5E）。由此可以看出，CK3/ST1中所富集的植物激素信號(hào)轉(zhuǎn)導(dǎo)和二萜生物合成兩條途徑與恨天高植株生長(zhǎng)關(guān)系較為密切，其他途徑則間接調(diào)控植株生長(zhǎng)。

KEGG富集研究結(jié)果顯示，CK3/ST1所富集的植物激素信號(hào)轉(zhuǎn)導(dǎo)途徑中色氨酸代謝CH3、SAUR下調(diào)表達(dá)，玉米素生物合成AHP上調(diào)表達(dá)、A-ARR下調(diào)表達(dá)，二萜生物合成途徑DELLA下調(diào)表達(dá)，油菜素內(nèi)酯合成TCH4上調(diào)表達(dá)、BSK和CYCD3下調(diào)表達(dá)（圖6A）。表明，植株生長(zhǎng)受多種激素互作的影響，本研究中可以認(rèn)為二萜生物合成途徑對(duì)植株生長(zhǎng)首先起作用。KEGG富集結(jié)果顯示，二萜生物赤霉素生物合成途徑中KAO上調(diào)表達(dá)、GA20ox和GA3ox下調(diào)表達(dá)、GA2ox既有上調(diào)也有下調(diào)表達(dá)（CK1/L1上調(diào)，CK3/ST1下調(diào)）（圖6B）。說(shuō)明，二萜生物合成途徑中DELLA、KAO、GA20ox、GA3ox和GA2ox調(diào)控恨天高植株生長(zhǎng)。

2.6轉(zhuǎn)錄組的實(shí)時(shí)熒光定量PCR驗(yàn)證

挑選5個(gè)差異表達(dá)基因，利用實(shí)時(shí)熒光定量PCR進(jìn)行對(duì)照組和試驗(yàn)組（800 mg/L GA3處理75 d）莖尖基因的表達(dá)量測(cè)定。結(jié)果顯示，實(shí)時(shí)熒光定量PCR（qPCR）與轉(zhuǎn)錄組測(cè)序（RNA-Seq）的結(jié)果在基因表達(dá)量變化幅度上具有一定差異，但基因的表達(dá)趨勢(shì)是一致的（圖7）。利用SPSS Pearson對(duì)qPCR和RNA-Seq進(jìn)行雙變量相關(guān)性分析，結(jié)果顯示5個(gè)差異表達(dá)基因的qPCR和RNA-Seq在P<0.01水平（雙側(cè)）上顯著相關(guān)，相關(guān)系數(shù)為0.978（表4）。由此可以說(shuō)明，轉(zhuǎn)錄組測(cè)序獲得的基因差異表達(dá)結(jié)果可靠。

3討論

恨天高作為天然矮生資源，是山茶矮化育種的優(yōu)良親本材料。本研究采用Illumina HiSeq高通量技術(shù)對(duì)0 mg/L和800 mg/L GA3處理75 d的葉片、莖段和莖尖組織進(jìn)行測(cè)序，組裝得到的181 115條Unigenes獲得了注釋。通過(guò)實(shí)時(shí)熒光定量PCR檢驗(yàn)，qPCR和RNA-Seq結(jié)果顯著相關(guān)，因此，本研究轉(zhuǎn)錄組測(cè)序結(jié)果是可靠的。

GO富集分析發(fā)現(xiàn)，大量差異表達(dá)基因在催化活性、結(jié)合、代謝過(guò)程、細(xì)胞過(guò)程、細(xì)胞組分和細(xì)胞器等類(lèi)別中顯著富集，在差異條目中，葉片、莖段和莖尖所富集的類(lèi)別存在差異，可能是外源赤霉素對(duì)不同組織的誘導(dǎo)作用存在差異，也可能是組織本身存在差異所造成。KEGG通路分析發(fā)現(xiàn)，葉片、莖段和莖尖顯著富集的代謝途徑不完全相同，其中苯丙烷生物合成、戊糖和葡萄糖醛酸相互轉(zhuǎn)化過(guò)程是共同顯著富集的代謝途徑。苯丙烷生物合成是植物最重要的次生代謝途徑之一[23]，而戊糖和葡萄糖醛酸相互轉(zhuǎn)化也是糖代謝的重要環(huán)節(jié)，參與植物多種發(fā)育過(guò)程[24]，表明在外源赤霉素誘導(dǎo)過(guò)程中，苯丙烷生物合成、戊糖和葡萄糖醛酸相互轉(zhuǎn)化途徑中差異基因?qū)尢旄咧仓晟L(zhǎng)具有重要作用。然而，相比于其他4組，CK3/ST1顯著富集的植物激素信號(hào)轉(zhuǎn)導(dǎo)、二萜生物合成和類(lèi)胡蘿卜素生物合成途徑中差異表達(dá)基因直接參與調(diào)控恨天高植株生長(zhǎng)。

本研究中，莖尖組織中外源赤霉素誘導(dǎo)的差異表達(dá)基因包括赤霉素信號(hào)傳導(dǎo)途徑中的DELLA，赤霉素生物合成途徑中的KAO、GA20ox、GA3ox和GA2ox等基因，DELLA、GA20ox、GA3ox和GA2ox下調(diào)表達(dá)，KAO上調(diào)表達(dá)。研究結(jié)果表明，DELLA蛋白在赤霉素信號(hào)傳導(dǎo)途徑中起中心開(kāi)關(guān)的作用，對(duì)赤霉素誘導(dǎo)的反應(yīng)非常靈敏，赤霉素處理5～10 min后蛋白水平顯著下降[25-26]，本研究中外源赤霉素的誘導(dǎo)信號(hào)傳遞給DELLA蛋白，DELLA蛋白呈下調(diào)表達(dá)，從而誘導(dǎo)下游赤霉素代謝基因如KAO的上調(diào)表達(dá)。

高等植物赤霉素水平通過(guò)正反饋和負(fù)反饋來(lái)實(shí)現(xiàn)調(diào)節(jié)體內(nèi)激素的平衡。當(dāng)體內(nèi)活性赤霉素水平低時(shí)，擬南芥GA20ox1、GA20ox2、GA20ox3和GA3ox1的表達(dá)量均顯著提高，外施GA3后，其表達(dá)量又顯著降低[27-28]。本研究中，外施赤霉素后恨天高莖尖中GA20ox、GA3ox也呈下調(diào)表達(dá)。當(dāng)活性赤霉素水平增加時(shí)，擬南芥GA2ox表達(dá)量上調(diào)[29]，在植物中過(guò)量表達(dá)GA2ox能誘導(dǎo)植株矮化，本研究中通過(guò)差異基因表達(dá)水平測(cè)定發(fā)現(xiàn)，莖尖中GA2ox下調(diào)表達(dá)，葉片中GA2ox呈上調(diào)表達(dá)，可能由于莖尖在外源赤霉素的誘導(dǎo)下，快速生長(zhǎng)，積累較多的活性赤霉素，當(dāng)活性赤霉素由莖尖傳輸?shù)饺~片時(shí)，則誘導(dǎo)葉片中的GA2ox上調(diào)表達(dá)，以抑制植株進(jìn)一步生長(zhǎng)，起到反饋調(diào)節(jié)作用。

另外，我們挑選了DELLA、GA20ox、GA3ox、GA2ox和KAO 5個(gè)差異表達(dá)基因進(jìn)行了實(shí)時(shí)熒光定量PCR驗(yàn)證，在莖尖組織中，實(shí)時(shí)熒光定量PCR（qPCR）與轉(zhuǎn)錄組測(cè)序（RNA-Seq）結(jié)果一致，DELLA、GA20ox、GA3ox、GA2ox呈下調(diào)表達(dá)，KAO上調(diào)表達(dá)。由此可以推測(cè)，800 mg/L外源赤霉素莖尖涂抹處理后，赤霉素質(zhì)量濃度升高，DELLA蛋白接受信號(hào)呈下調(diào)表達(dá)，誘導(dǎo)下游KAO上調(diào)表達(dá)，GA20ox、GA3ox、GA2ox呈下調(diào)表達(dá)，多個(gè)基因相互作用，促進(jìn)植株生長(zhǎng)，當(dāng)體內(nèi)赤霉素質(zhì)量濃度積累到一定程度時(shí)，GA2ox則反饋調(diào)節(jié)抑制植株繼續(xù)生長(zhǎng)。外源赤霉素誘導(dǎo)植株生長(zhǎng)的機(jī)理，以本文為基礎(chǔ)，后續(xù)還需要進(jìn)一步深入研究。

本研究利用高通量轉(zhuǎn)錄組測(cè)序技術(shù)建立矮生山茶恨天高赤霉素誘導(dǎo)前后葉片、莖段和莖尖轉(zhuǎn)錄組數(shù)據(jù)庫(kù)，分析比較3種組織在外源赤霉素處理后生物反應(yīng)途徑的差異，發(fā)現(xiàn)相對(duì)于葉片和莖段，莖尖對(duì)植株生長(zhǎng)影響較大，進(jìn)一步篩選出莖尖組織中的差異表達(dá)基因，為研究恨天高矮生機(jī)制以及培育矮生山茶品種提供基礎(chǔ)。

參考文獻(xiàn)：

[1]RICHAEDS D E， KING K E， AIT-ALI T， et al. How gibberellin regulates plant growth and development： a molecular genetic analysis of gibberellin signaling [J]. Annual Review of Plant Biology， 2001， 52： 67-88.

[2]張旭紅，孫美玉，李靖銳，等. 東方百合‘索邦GA20ox的克隆及表達(dá)分析[J]. 園藝學(xué)報(bào)， 2019， 46（1）： 74-86.

[3]YAMAGUCHI S. Gibberellin metabolism and its regulation [J]. Annual Review of Plant Biology， 2008， 59： 225-251.

[4]DAVIDSON S E， ELLIOTT R C， HELLIWELL C A， et al. The pea gene NA encodes ent-kaurenoic acid oxidase [J]. Plant Physiology， 2003， 131（1）： 335-344.

[5]LANGE T， HEDDEN P， GRAEBE J E. Expression cloning of a gibberellin 20-oxidase， a multifunctional enzyme involved in gibberellin biosynthesis [J]. Proceedings of the National Academy of Sciences， 1994， 91（18）： 8522-8566.

[6]APPLEFORD N E J， LENTON J R. Gibberellins and leaf expansion in near-isogenic wheat lines containing Rht1 and Rhe3 dwarfing alleles [J]. Planta， 1991， 183（2）： 229-236.

[7]COLES J P， PHILLIPS A L， CROKER S J， et al. Modification of gibberellin production and plant development in Arabidopsis by sense and antisense expression of gibberellin 20-oxidase genes [J]. The Plant Journal， 1999， 17（5）： 547-563.

[8]SAKAMOTO T， MIURA K， ITOH H， et al. An overview of gibberellin metabolism enzyme genes and their related mutants in rice [J]. Plant Physiology， 2004， 134（4）： 1642-1653.

[9]歐春青，姜淑苓，王斐，等. 梨貝殼杉烯酸氧化酶基因PcKAO1的克隆與表達(dá)分析[J]. 園藝學(xué)報(bào)， 2013， 40（5）： 849-858.

[10]WANG J Y， WU B， LI J Y， et al. Overexpression and silent expression of CrGA20ox1 from Camellia reticulate ‘Hentiangao and its effect on morphological alterations in transgenic tobacco plants [J]. Plant Breeding， 2018， 9（3）： 1-9.

[11]HUANG Y， WANG X， GE S， et al. Divergence and adaptive evolution of the gibberellin oxidase genes in plants [J]. BMC Evolutionary Biology， 2015， 15： 207-221.

[12]THOMAS S G， PHILLIPS A L， HEDDEN P. Molecular cloning and functional expression of gibberellin 2-oxidases， multifunctional enzymes involved in gibberellin deactivation [J]. Proceedings of the National Academy of Sciences， 1999， 96（8）： 4698-4703.

[13]李巧峽，張麗，王玉，等. 赤霉素調(diào)控植物開(kāi)花及器官發(fā)育的研究進(jìn)展[J]. 中國(guó)細(xì)胞生物學(xué)學(xué)報(bào)， 2019， 41（4）： 746-758.

[14]WILLIGE B C， GHOSH S， NILL C， et al. The DELLA domain of GA INSENSITIVE mediates the interaction with the GA INSENSITIVE DWARF1A gibberellin receptor of Arabidopsis [J]. The Plant Cell， 2007， 19（4）： 1209-1220.

[15]DILL A， THOMAS S G， HU J， et al. The Arabidopsis F-Box protein SLEEPY1 targets gibberellin signaling repressors for gibberellin-induced degradation [J]. The Plant Cell， 2004， 16（6）： 1392-1405.

[16]SUN T P. Gibberellin-GID1-DELLA： a pivotal regulatory module for plant growth and development [J]. Plant Physiology， 2010， 154（2）： 567-570.

[17]CHEN J H， CHENG T L， WANG P K， et al. Genome-wide bioinformatics analysis of DELLA-family proteins from plants [J]. Plant Omics， 2013， 6（3）： 201-207.

[18]張文穎，王晨，朱旭東. 葡萄全基因組DELLA 蛋白基因家族鑒定及其應(yīng)答外源赤霉素調(diào)控葡萄果實(shí)發(fā)育的特征[J]. 中國(guó)農(nóng)業(yè)科學(xué)， 2018， 51（16）： 3130-3146.

[19]LI W J， ZHANG J X， SUN H Y. FveRGA1， encoding a DELLA protein， negatively regulates runner production in Fragaria vesca [J]. Planta， 2018， 247（4）： 941-951.

[20]賈永鵬，李開(kāi)祥，昝領(lǐng)兄，等. 甘藍(lán)型油菜全基因組DELLA蛋白基因家族的鑒定和表達(dá)分析[J]. 中國(guó)油料作物學(xué)報(bào)， 2019， 41（3）： 360-368.

[21]管開(kāi)云.中國(guó)茶花圖鑒[M]. 杭州： 浙江科學(xué)技術(shù)出版社， 2014：53-54.

[22]張恩亮，馬玲玲，楊如同，等. IBA誘導(dǎo)楸樹(shù)嫩枝扦插不定根發(fā)育的轉(zhuǎn)錄組分析[J]. 林業(yè)科學(xué)， 2018， 54（5）： 48-61.

[23]WANG Y， WEATHERS P J. Sugar proportionately affect artemisinin production [J]. Plant Cell Report， 2007， 26（7）： 1073-1081.

[24]HU Y L， GAI Y， YIN L， et al. Crystal structures of a Populus tomentosa 4- coumarate： CoA ligase shed light on its enzymatic mechanisms[J]. Plant Cell， 2010，22（9）：3093-3104.

[25]ITOH H， UEGUCHI-TANAKA M， SATO Y， et al. The gibberellin signaling pathway is regulated by the appearance and disappearance of SLENDER RICE1 in nuclei [J]. Plant Cell， 2002， 14（1）： 57-70.

[26]TYLER L， THOMAS S G， HU J， et al. DELLA proteins and gibberellin regulated seed germination and floral development in Arabidopsis [J]. Plant Physiology， 2004， 135（2）： 1008-1019.

[27]MITCHUM M G， YAMAGUCHI S， HANADA A， et al. Distinct and overlapping roles of two gibberellin 3-oxidases in Arabidopsis development [J]. The Plant Journal， 2006， 45（5）： 804-818.

[28]RIEU I， RUIZ-RIVERO O， FERNANDEZ-GARCIA N， et al. The gibberellin biosynthetic genes At GA20ox1 and At GA20ox2 act， partially redundantly， to promote growth and development throughout the Arabidopsis life cycle [J]. The Plant Journal， 2008， 53（3）： 488-504.

[29]RIEU I， THOMAS S， POWERS S J， et al. Genetic analysis reveals that C19-GA 2-oxidation is a major gibberellin inactivation pathway in Arabidopsis [J]. Plant Cell， 2008， 20（9）： 2420-2436.

（責(zé)任編輯：陳海霞）