郭由兵，程廷才，李開和，盧燕回，王衛(wèi)峰，夏慶友*

（1.家蠶基因組生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室，西南大學(xué)，重慶 400715；2.中國煙草總公司廣西壯族自治區(qū)公司，南寧 530022）

打頂期煙草組織基因表達(dá)分析及打頂對(duì)腋芽基因表達(dá)的影響

郭由兵1，程廷才1，李開和2，盧燕回2，王衛(wèi)峰2，夏慶友1*

（1.家蠶基因組生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室，西南大學(xué)，重慶 400715；2.中國煙草總公司廣西壯族自治區(qū)公司，南寧 530022）

打頂是煙葉生產(chǎn)重要的技術(shù)環(huán)節(jié)之一。為了研究煙草（Nicotiana tabacum）打頂期基因表達(dá)模式，選取廣西大寧地區(qū)曬黃煙現(xiàn)蕾期10個(gè)打頂或非打頂?shù)慕M織樣本，構(gòu)建鏈特異總RNA測(cè)序文庫，利用Illumina HiSeq 2000完成測(cè)序。通過質(zhì)量控制，獲得9千多萬條測(cè)序數(shù)據(jù)，分析獲得50 706個(gè)基因的表達(dá)信息。篩選組織特異和打頂差異表達(dá)基因，共富集687個(gè)組織特異表達(dá)基因和444個(gè)腋芽打頂差異表達(dá)基因，為煙草組織基因表達(dá)模式和打頂分子響應(yīng)機(jī)制研究提供數(shù)據(jù)支撐，也為調(diào)控?zé)煵轄I養(yǎng)和生殖生長提供分子靶標(biāo)。

煙草；打頂；轉(zhuǎn)錄組；生長調(diào)控

煙草（Nicotiana tabacum）是以收獲煙葉為主的經(jīng)濟(jì)作物，在其營養(yǎng)生長后期，莖頂形成生殖中心，大量營養(yǎng)物質(zhì)集中到頂端，影響煙葉產(chǎn)量和品質(zhì)[1]。在煙葉生產(chǎn)過程中，打頂抹杈是抑制生殖生長和促進(jìn)營養(yǎng)生長的有效措施，將引起煙草植株一系列生理生化響應(yīng)，如增加營養(yǎng)元素的吸收利用、代謝物質(zhì)的積累[1-2]。研究人員利用miRNA測(cè)序技術(shù)分析煙草根部對(duì)打頂響應(yīng)的分子機(jī)制，發(fā)現(xiàn)打頂差異基因主要涉及激素代謝、轉(zhuǎn)錄、應(yīng)激反應(yīng)等[3-4]。楊潔等[5]報(bào)道打頂可影響腋芽內(nèi)激素的含量和腋芽的生長，而楊惠娟等[6]報(bào)道了打頂差異蛋白能促進(jìn)光合作用、增強(qiáng)植株抗性。隨著測(cè)序技術(shù)的發(fā)展，3個(gè)煙草品種[Basma Xanthi (BX, Oriental)，K326 (Flue-cured)，TN90 (Burley)]的基因組測(cè)序完成，共鑒定到9萬多個(gè)基因[7]。但迄今為止，打頂期全基因組范圍的基因表達(dá)模式卻鮮有報(bào)道。本研究采用RNA-seq技術(shù)，對(duì)煙草根、莖、葉、腋芽、花蕾、頂芽等組織進(jìn)行轉(zhuǎn)錄組分析，比較打頂和非打頂樣品的基因表達(dá)差異，為揭示煙草打頂響應(yīng)分子機(jī)制和調(diào)控營養(yǎng)和生殖生長提供數(shù)據(jù)支撐。

1 材料與方法

1.1 試驗(yàn)材料

本試驗(yàn)選用廣西大寧地區(qū)曬黃煙作為試驗(yàn)材料。收集現(xiàn)蕾期打頂3 d的根、莖、葉、腋芽等4份樣品，同期未打頂根、莖、葉、腋芽、頂芽和花蕾等6份樣品，共10份組織樣本。所有樣本用液氮保存，用于提取總RNA。

1.2 Illumina測(cè)序文庫構(gòu)建

參照TRIZOL Reagent (Cat#15596-018, Life technologies, Carlsbad, CA, US)提取總RNA，采用RNase-Free DNase Set (Cat#79254, QIAGEN, GmBH, Germany)純化。使用Agilent 2100檢查RNA完整性和純度。采用TruSeq? Stranded Total RNA kit (Cat#RS-122-2201, Illumina, CA, US)構(gòu)建Paired-end測(cè)序文庫，使用Illumina HiSeq 2000完成測(cè)序。

1.3 轉(zhuǎn)錄組數(shù)據(jù)分析

使用軟件fqtrim (v0.9.4, http://ccb.jhu.edu/ software/fqtrim/index.shtml)去除ployA，Bowtie（version 2.2.5）[8]去除rRNA序列，Trimmamotic（Version 0.36）[9]去除接頭、低質(zhì)量堿基和長度較短的讀段(Reads)，fastqc (v0.11.2, http://www. bioinformatics.babraham.ac.uk/projects/fastqc/)檢測(cè)數(shù)據(jù)質(zhì)量。處理后的數(shù)據(jù)使用RSEM[10](v1.2.22)比對(duì)到煙草K326（Nicotiana tabacum K326）的參考基因集上，對(duì)樣本中的基因表達(dá)定量。煙草K326的參考基因集數(shù)據(jù)下載自煙草基因組數(shù)據(jù)庫(http://solgenomics.net)[11]。通過腳本程序abundance_estimates_to_matrix.pl將基因表達(dá)結(jié)果整合成表達(dá)矩陣，并采用TMM（M值的截尾均值，Trimmed Mean of M-values）方法均一化。采用軟件edgeR[12]鑒定樣本間差異表達(dá)的基因，腳本程序get_tissue_enriched_DE_one_vs_all.pl分析各樣本富集的高表達(dá)基因。

1.4 GO和KEGG注釋與富集分析

基因功能注釋使用Trinonate流程(v2.0.2, http:// trinotate.github.io/)完成，采用KAAS（KEGG自動(dòng)注釋系統(tǒng)，KEGG Automatic Annotation Server）[13]對(duì)基因進(jìn)行KEGG（京都基因和基因組百科全書，Kyoto Encyclopedia of Genes and Genomes）注釋。通過WEGO (在線基因本體注釋圖示，Web Gene Ontology Annotation Plot)[14]比較打頂差異表達(dá)基因的功能。GO（基因本體，Gene Ontology）和KEGG富集通過程序GO_seq[15]軟件完成。以上perl腳本來自Trinity[16]軟件，相關(guān)圖表由軟件R和Excel生成。

2 結(jié) 果

2.1 轉(zhuǎn)錄組數(shù)據(jù)統(tǒng)計(jì)

經(jīng)過質(zhì)量控制的測(cè)序數(shù)據(jù)比對(duì)到參考基因集上，共獲得50 706個(gè)基因的表達(dá)信息。由于參考序列僅為基因集，不能代表全基因組，因此各個(gè)樣本的比對(duì)率較低（表1）。樣本間比對(duì)率存在差異，如 打頂樣本中，葉和根的比對(duì)率為23.38%和31.94%，其余樣本的比對(duì)率都在35%以上。以FPKM（每百萬測(cè)序片段覆蓋每千堿基片段數(shù)，F(xiàn)ragments Per Kilobase per Million）值大于等于2為限，各樣本表達(dá)的基因占總基因數(shù)的2/3（66.84%）以上，表達(dá)的轉(zhuǎn)錄本占總轉(zhuǎn)錄本數(shù)的1/3（36.63%）以上。為避免高表達(dá)和零表達(dá)基因的干擾，對(duì)FPKM矩陣進(jìn)行l(wèi)og10(FPKM+1)轉(zhuǎn)化（圖1），結(jié)果顯示基因在各樣本中的表達(dá)相對(duì)均衡。上述結(jié)果表明，盡管各樣本測(cè)序數(shù)據(jù)和比對(duì)率存在差異，但基因表達(dá)相對(duì)均衡，其結(jié)果可用于下游分析。

表1 RNA-seq數(shù)據(jù)統(tǒng)計(jì)Table 1 Summary of RNA-seq data

圖1 各樣本中表達(dá)基因分布Fig. 1 Distribution of expressed genes in each sample

2.2 組織特異基因表達(dá)模式分析

通過差異表達(dá)基因分析，獲得煙草組織特異表達(dá)的基因，即在某一組織高表達(dá)且與其他所有組織都有顯著差異的基因。利用組織特異表達(dá)基因的FPKM值繪制熱圖（圖2），結(jié)果顯示根部含有數(shù)目較多的特異表達(dá)基因，暗示地上部分和地下部分在基因表達(dá)上的差異，在轉(zhuǎn)錄水平上揭示了煙草根部組織的特異性。葉片中存在一類高表達(dá)的基因，且在花蕾、腋芽、芽和莖中也有表達(dá)，功能注釋顯示這些基因可能與光合作用及色素合成代謝相關(guān)。

2.3 組織特異表達(dá)基因功能注釋

圖2 組織差異基因表達(dá)模式Fig. 2 Expression patterns of tissue-specific genes

通過軟件GO_seq對(duì)特異表達(dá)基因進(jìn)行GO富集。相較于其他組織，根部有較多的特異表達(dá)基因，主要涉及響應(yīng)刺激物（GO:0050896 response to stimulus）、物質(zhì)運(yùn)輸（GO:0006810 transport）、氧化還原（GO:0055114 oxidation- reduction process）、天然免疫（GO:0045087 innate immune response）和細(xì)胞程序性死亡（GO: 0012501 programmed cell death）等。其中煙草生物堿合成過程中的腐胺N甲基轉(zhuǎn)移酶全部在根部特異表達(dá)，和生物堿合成相關(guān)的腐胺、聚胺等物質(zhì)合成相關(guān)基因也在根部特異表達(dá)，表明煙草根系是生物堿合成的重要場(chǎng)所。氧化還原過程中的血紅素結(jié)合、氧化還原酶、谷胱甘肽轉(zhuǎn)移酶等相關(guān)酶類，在根部都特異表達(dá)，暗示根部存在特殊的氧化還原保護(hù)機(jī)制。這與研究人員報(bào)道的生物堿和抗氧化基因在根部上調(diào)表達(dá)是一致的[3,7]。葉片特異表達(dá)的基因主要和光合作用相關(guān)，組成包括葉綠體、類囊體膜以及光合體系等組分，涉及光合作用、光呼吸、丙酮酸代謝等生物學(xué)過程，包括果糖二磷酸醛縮酶、醛裂解酶、果糖1,6-二磷酸1-磷酸酶等功能基因。這些基因的特異表達(dá)是葉片進(jìn)行光合作用的分子基礎(chǔ)?；ɡ偬禺惐磉_(dá)的基因主要涉及苯酸鹽羥化降解、天冬氨酸代謝等生物學(xué)過程，相應(yīng)的分子功能為酵母賴氨酸脫氫酶（GO:0004753 saccharopine dehydrogenase activity）、丙氨酸乙醛酸轉(zhuǎn)氨酶（GO:0008453 alanine- glyoxylate transaminase activity）、alpha-L-阿拉伯呋喃糖苷酶（GO:0046556 alpha-L- arabinofuranosidase activity）和天冬氨酰合酶（GO: 0004066 asparagine synthase activity）等。這些基因的表達(dá)，暗示花蕾中存在復(fù)雜的物質(zhì)代謝過程。

腋芽位于煙草莖上葉的腋部，具有萌發(fā)成分枝的潛力。腋芽組織特異表達(dá)的基因主要是與環(huán)境響應(yīng)、轉(zhuǎn)錄和幾丁質(zhì)酶活性等功能相關(guān)，主要包括水份或壓力應(yīng)激相應(yīng)，DNA結(jié)合轉(zhuǎn)錄調(diào)控因子等功能，可能與基因表達(dá)調(diào)控相關(guān)。參與幾丁質(zhì)酶降解過程的幾丁質(zhì)酶也在腋芽中富集表達(dá)，幾丁質(zhì)是僅次于纖維素的細(xì)胞內(nèi)最豐富的有機(jī)化合物，其降解在細(xì)胞營養(yǎng)循環(huán)中發(fā)揮著重要作用。這些結(jié)果暗示打頂期的腋芽正加快物質(zhì)循化和功能分化。

利用KAAS對(duì)基因集進(jìn)行KEGG功能注釋，共注釋11 318個(gè)基因。組織特異表達(dá)基因的KEGG富集分析顯示，葉中富集到的ALDO、rbcS、GAPA、FBP等基因與光合作用相關(guān)（表2），GGAT參與碳固定和氨基酸代謝[17]，DOX1參與長鏈脂肪酸a氧化的起始步驟，與植物抗菌或其他病原反應(yīng)相關(guān)[18]。

表2 組織特異表達(dá)基因KEGG注釋富集Table 2 Enrichment of KEGG annotation of tissue-specific genes

2.4 打頂差異基因分析

利用edgeR軟件，以樣本間基因表達(dá)差異倍數(shù)大于等于4，且基因在樣本間差異顯著的P值小于等于0.001（即P-vlaue≤1e-3和|log2(FC)|≥2）為標(biāo)準(zhǔn)，計(jì)算樣本間的基因表達(dá)差異。比較葉、莖、腋芽、根四組打頂與未打頂樣本的基因表達(dá)差異，發(fā)現(xiàn)腋芽中差異基因最多，莖、根中存在少量差異基因，暗示打頂對(duì)腋芽的影響遠(yuǎn)大于其他組織。同時(shí)，打頂腋芽和頂芽僅有16個(gè)差異表達(dá)基因，暗示打頂后腋芽的頂端優(yōu)勢(shì)抑制解除，替代頂芽成為新的生長點(diǎn)。以未打頂?shù)囊秆繛閰⒄諛颖荆▓D3），打頂腋芽和未打頂腋芽相比，上調(diào)152個(gè)基因，其中105個(gè)在芽中同樣也上調(diào)，約占69.08%，下調(diào)292個(gè)基因，其中245個(gè)在芽中也同樣上調(diào)，約占83.90%，暗示這些基因在頂芽萌發(fā)和打頂后腋芽發(fā)育中起著關(guān)鍵作用。

圖3 打頂處理的腋芽和頂芽差異表達(dá)基因數(shù)量比較Fig. 3 Comparison of differentially expressed genes in apical buds and axillary buds upon topping

2.5 打頂差異基因的功能注釋

比較腋芽打頂差異表達(dá)基因的GO功能注釋（圖4），結(jié)果顯示打頂上調(diào)的基因主要分布在GO注釋的細(xì)胞組分類目下，涉及細(xì)胞壁、細(xì)胞質(zhì)部分、蛋白復(fù)合物、轉(zhuǎn)運(yùn)載體和形態(tài)學(xué)構(gòu)成等。打頂下調(diào)的基因主要分布在GO注釋的分子功能和生物學(xué)過程類目下，涉及的生物學(xué)過程有DNA結(jié)合，蛋白質(zhì)結(jié)合，RNA合成和轉(zhuǎn)錄，應(yīng)激和耐受性，激素響應(yīng)和抗毒素的合成等。RNA合成和轉(zhuǎn)錄受抑制，延緩細(xì)胞功能的進(jìn)一步分化；激素和抗毒素的合成下降，減緩頂端優(yōu)勢(shì)，促進(jìn)腋芽生長。腋芽中的細(xì)胞組分類蛋白上調(diào)則是腋芽向營養(yǎng)生長轉(zhuǎn)變的開始，暗示打頂后腋芽由先前的功能分化向營養(yǎng)生長轉(zhuǎn)變。

圖4 腋芽打頂差異表達(dá)基因GO功能注釋Fig. 4 GO functional output of differentially expressed genes after topping in axillary buds

表3 腋芽打頂差異表達(dá)基因KEGG注釋富集Table 3 Enrichment of KEGG annotation differentially expressed genes in axillary buds after topping

依據(jù)KAAS的注釋結(jié)果，對(duì)腋芽打頂差異表達(dá)基因進(jìn)行KEGG富集分析（表3）。與未打頂腋芽相比，打頂腋芽中的基因表達(dá)模式發(fā)生了顯著變化。上調(diào)的組蛋白H2A、H4為核小體組裝因子，可導(dǎo)致核內(nèi)染色質(zhì)結(jié)構(gòu)重塑[19-20]；RPS10為線粒體中核糖體的組成部分之一，上調(diào)表達(dá)導(dǎo)致線粒體大小亞基比例失調(diào)，影響植株的休眠與復(fù)蘇[21-22]；TIP和PIP為水通道蛋白，可維系細(xì)胞內(nèi)水分平衡，其上調(diào)表達(dá)可增進(jìn)細(xì)胞對(duì)水分的吸收[23-24]；LHCA3為光合作用觸角蛋白，以較慢的速率吸收紅光[25]，暗示打頂后腋芽中光合作用也受到影響。打頂腋芽下調(diào)的TPS和sacA參與海藻糖和蔗糖的代謝，有研究顯示TPS與冬小麥耐寒特性相關(guān)[26]；WRKY33和HSFF為轉(zhuǎn)錄因子，WRKY33參與植物與病原菌的相互作用[27]；PRR5抑制植物的避蔭反應(yīng)，影響植物周期節(jié)律[28]。

3 討 論

本研究采集打頂期10個(gè)組織樣本，基于RNA-Seq技術(shù)，獲得50 706個(gè)基因的表達(dá)信息。由于公共數(shù)據(jù)庫中煙草基因組的注釋信息并不完整，因此，將測(cè)序數(shù)據(jù)比對(duì)到基因集上，分析基因表達(dá)模式。分析中測(cè)序數(shù)據(jù)比對(duì)到基因集上的比對(duì)率較低，可能是由于參考基因集僅來自根和葉[7]，且根、葉與環(huán)境接觸較多，取樣過程中難免存在污染，同時(shí)基因表達(dá)具有的時(shí)期和條件的特異性，共同導(dǎo)致分析的組織中根和葉的比對(duì)率不高。即使如此，以FPKM值大于等于2為限，超過總基因和轉(zhuǎn)錄本的2/3（66.84%）和1/3（36.63%）是表達(dá)的。各樣本表達(dá)基因數(shù)目和FPKM分布相對(duì)均衡，數(shù)據(jù)可以用于比較分析。

分析組織特異表達(dá)基因發(fā)現(xiàn)，根和葉中含有較多的特異基因（351和208個(gè)），這與根和葉具有活躍代謝活動(dòng)的組織特性是相符的[7]。與其他組織相比，根部有較多的特異基因表達(dá)，在轉(zhuǎn)錄水平上揭示了地上組織與地下組織的區(qū)別。對(duì)特異表達(dá)基因功能注釋發(fā)現(xiàn)，根部特異表達(dá)的基因多參與環(huán)境響應(yīng)、氧化還原以及物質(zhì)運(yùn)輸，從轉(zhuǎn)錄水平反映了根與環(huán)境的緊密接觸和在植株攝取營養(yǎng)中的關(guān)鍵作用。另外，生物堿合成相關(guān)的基因在根部富集表達(dá)，暗示根部是生物堿合成的重要場(chǎng)所，這與先前的報(bào)道是一致的[7]。作為光合作用的主要場(chǎng)所，葉片中富集表達(dá)了一系列的酶、葉綠體組分、暗反應(yīng)等光合作用相關(guān)的基因。花蕾作為生殖器官，富集表達(dá)的基因中較多地涉及苯酸鹽羥化降解、天冬氨酸類氨基酸代謝等相關(guān)過程，暗示花蕾中存在較為復(fù)雜的生殖代謝過程[29]。與此對(duì)應(yīng)，腋芽中特異表達(dá)的基因除與環(huán)境相關(guān)，還包括轉(zhuǎn)錄因子，能調(diào)控相關(guān)基因的表達(dá)，暗示腋芽有進(jìn)一步分化或增殖的潛力[2,30]。

打頂抹杈是煙草生產(chǎn)的重要技術(shù)環(huán)節(jié)，可以打破生殖與營養(yǎng)生長的平衡，使較多的營養(yǎng)富集到葉片中，從而提高煙葉產(chǎn)量和品質(zhì)。比較打頂對(duì)各組織中基因表達(dá)的影響，發(fā)現(xiàn)腋芽中含有較多的差異基因，根、莖及葉等組織中基因表達(dá)差異不明顯，打頂后腋芽中的基因表達(dá)模式更趨向于頂芽。這些結(jié)果暗示打頂后如不及時(shí)抹杈，腋芽可能取代頂芽成為新的生長中心，根、葉等組織不再產(chǎn)生明顯變化，無法實(shí)現(xiàn)提高葉片產(chǎn)量的目標(biāo)，從轉(zhuǎn)錄水平證實(shí)了打頂后抹杈的必要[31-34]。

4 結(jié) 論

本研究采用轉(zhuǎn)錄組測(cè)序技術(shù)，對(duì)10個(gè)煙草組織樣本進(jìn)行轉(zhuǎn)錄組測(cè)序，獲得超過5萬個(gè)基因的表達(dá)信息，揭示各組織樣本的基因表達(dá)特征以及植株對(duì)打頂?shù)捻憫?yīng)機(jī)制。結(jié)果表明根部主要富集生物堿等物質(zhì)合成代謝相關(guān)基因，葉片富集光合作用相關(guān)基因。比較打頂差異表達(dá)基因，顯示腋芽具有440余差異表達(dá)基因，明顯高于其他組織樣本，暗示打頂主要影響煙草腋芽生長發(fā)育。

[1] 林桂華，周冀衡，范啟福，等. 打頂技術(shù)對(duì)烤煙產(chǎn)質(zhì)量和生物堿組成的影響[J]. 中國煙草科學(xué)，2002，23（4）：9-13.

[2] 劉進(jìn)平. 植物腋芽生長與頂端優(yōu)勢(shì)[J]. 植物生理學(xué)通訊，2007（3）：575-582.

[3] Qi Y, Guo H, Li K, et al. Comprehensive analysis of differential genes and miRNA profiles for discovery of topping-responsive genes in flue-cured tobacco roots[J]. FEBS Journal, 2012, 279(6): 1054-1070.

[4] Tang S, Wang Y, Li Z, et al. Identification of wounding and topping responsive small RNAs in tobacco (Nicotiana tabacum)[J]. BMC Plant Biol, 2012, 12: 28.

[5] 楊潔，胡日生，童建華，等. 打頂對(duì)煙草腋芽生長及植物激素含量的影響[J]. 煙草科技，2013（10）：72-75.

[6] 楊惠娟，王景，王紅麗，等. 打頂前后烤煙葉片microRNAs表達(dá)差異的研究[J]. 中國煙草學(xué)報(bào)，2014，20（5）：110-115.

[7] Sierro N, Battey J N, Ouadi S, et al. The tobacco genome sequence and its comparison with those of tomato and potato[J]. Nat Commun, 2014, 5: 3833.

[8] Langmead B, Trapnell C, Pop, M, et al. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome[J]. Genome Biology, 2009, 10(3): R25-R25.

[9] Bolger A M, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data[J]. Bioinformatics, 2014, 30(15): 2114-2120.

[10] Li B, Dewey C N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome[J]. BMC Bioinformatics, 2011, 12: 323.

[11] Fernandez-Pozo N, Menda N, Edwards J D, et al. The Sol Genomics Network (SGN)--from genotype to phenotype to breeding[J]. Nucleic Acids Res, 2015, 43(Database issue): D1036-1041.

[12] Robinson M D, McCarthy D J, Smyth G K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data[J]. Bioinformatics, 2010, 26(1): 139-140.

[13] Moriya Y, Itoh M, Okuda S, et al. KAAS: an automatic genome annotation and pathway reconstruction server[J]. Nucleic Acids Res, 2007, 35(Web Server issue): W182-185.

[14] Ye J, Fang L, Zheng H, et al. WEGO: a web tool for plotting GO annotations[J]. Nucleic Acids Res, 2006, 34(Web Server issue): W293-297.

[15] Young, M D, Wakefield M J, Smyth G K, et al. Gene ontology analysis for RNA-seq: accounting for selection bias[J]. Genome Biol, 2010, 11(2): R14.

[16] Haas B J, Papanicolaou A, Yassour M, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis[J]. Nat Protocols, 2013, 8(8): 1494-1512.

[17] Liepman A H, Olsen, L J Alanine aminotransferase homologs catalyze the glutamate: glyoxylate aminotransferase reaction in peroxisomes of Arabidopsis[J]. Plant Physiol, 2003, 131(1): 215-227.

[18] Shimada T L, Takano Y, Shimada T, et al. Leaf oil body functions as a subcellular factory for the production of a phytoalexin in Arabidopsis[J]. Plant Physiol, 2014, 164(1): 105-118.

[19] Tolstorukov M Y, Goldman, J A, Gilbert C, et al. Histone variant H2A. Bbd is associated with active transcription and mRNA processing in human cells[J]. Mol Cell, 2012, 47(4): 596-607.

[20] Miotto B, Struhl K. HBO1 histone acetylase activity is essential for DNA replication licensing and inhibited by Geminin[J]. Mol Cell, 2010, 37(1): 57-66.

[21] Kwasniak M, Majewski P, Skibior R, et al. Silencing of the nuclear RPS10 gene encoding mitochondrial ribosomal protein alters translation in Arabidopsis mitochondria[J]. Plant Cell, 2013, 25(5): 1855-1867.

[22] Majewski P, Woloszynska M, Janska H. Developmentally early and late onset of Rps10 silencing in Arabidopsis thaliana: genetic and environmental regulation[J]. J Exp Bot, 2009, 60(4): 1163-1178.

[23] Venkatesh J, Yu J W, Park S W. Genome-wide analysis and expression profiling of the Solanum tuberosum aquaporins[J]. Plant Physiol Biochem, 2013, 73: 392-404.

[24] Reuscher S, Akiyama M, Mori C, et al. Genome-Wide Identification and Expression Analysis of Aquaporins in Tomato[J]. PLoS ONE, 2013, 8(11): e79052.

[25] Wientjes E, van Stokkum I H, van Amerongen H, et al. The role of the individual Lhcas in photosystem I excitation energy trapping[J]. Biophys J, 2011, 101(3): 745-754.

[26] Xie D W, Wang X N, Fu L S, et al. Identification of the trehalose-6-phosphate synthase gene family in winter wheat and expression analysis under conditions of freezing stress[J]. J Genet, 2015, 94(1): 55-65.

[27] Eulgem T, Rushton P J, Robatzek S, et al. The WRKY superfamily of plant transcription factors[J]. Trends Plant Sci, 2000, 5(5): 199-206.

[28] Takase M, Mizoguchi T, Kozuka T, et al. The unique function of the Arabidopsis circadian clock gene PRR5 in the regulation of shade avoidance response[J]. Plant Signal Behav, 2013, 8(4): e23534.

[29] Bokvaj P, Hafidh S, Honys D. Transcriptome profiling of male gametophyte development in Nicotiana tabacum[J]. Genom Data, 2015, 3: 106-111.

[30] 楊潔. 植物激素與煙草腋芽生長的關(guān)系及其調(diào)控研究[D]. 長沙：湖南農(nóng)業(yè)大學(xué)，2013.

[31] 楊惠娟，許俐，王景，等. 打頂前后烤煙葉片蛋白質(zhì)表達(dá)差異的研究[J]. 中國煙草學(xué)報(bào)，2012，18（5）：73-78.

[32] 陳常偉，姚忠達(dá)，張新華，等. 打頂時(shí)期及仿生型信號(hào)分子施用時(shí)間對(duì)抑制煙草生物堿效應(yīng)的影響[J]. 中國煙草學(xué)報(bào)，2011，21（6）：32-36.

[33] 王茂. 打頂前后煙草miRNA表達(dá)譜的生物信息學(xué)分析及靶基因的電子克隆[D]. 鄭州：河南農(nóng)業(yè)大學(xué)，2010.

[34] 林中麟，石健林，周益. 煙草打頂研究進(jìn)展[J]. 江西農(nóng)業(yè)學(xué)報(bào)，2009（6）：32-36.

Analysis of Expressed Genes at the Topping Stage and Effects of Topping to Genes Expression in Axillary Buds in Nicotiana tabacum

GUO Youbing1, CHENG Tingcai1, LI Kaihe2, LU Yanhui2, WANG Weifeng2, XIA Qingyou1*

(1. State Key Laboratory of Silkworm Genome Biology, Southwest University, Chongqing 400715, China; 2. China Tobacco Guangxi Zhuang Autonomous Region Corporation, Nanning 530022, China)

Topping is one of the important technical procedures of Nicotiana tabacum production. In order to study the gene expression patterns at the topping stage in Nicotiana tabacum, we selected the Daning tobacco growing area in Guangxi and collected ten tissue samples at the topping stage. The samples were used to extract total RNAs and to construct stranded specific total RNA sequencing library and to sequence using the Illumina HiSeq 2000 sequencer. Through quality control, more than 91 million clean reads were obtained and a total of 50706 genes were detected to be expressed from ten tissues at the topping stage. After normalization, screening for tissue-specific and genes differentially expressed among tested tissues, a total of 687 genes were identified to be tissue-specifical and 444 genes showed differential expression before and after topping. Above data and results will not only provide some support for exploring the molecular responses of topping, but also provide molecular targets for the regulatiion of Nicotiana tabacum vegetative and reproductive growth.

Nicotiana tabacum ; topping; transcriptome; growth regulation

S572.03

1007-5119（2016）06-0014-07

10.13496/j.issn.1007-5119.2016.06.003

廣西壯族自治區(qū)煙草公司科技創(chuàng)新項(xiàng)目“廣西賀州大寧曬煙基因組測(cè)序、表達(dá)譜芯片和代謝組學(xué)研究”(GYK201226)

郭由兵（1989-），男，碩士研究生，研究方向?yàn)樯镄畔W(xué)。E-mail：gyoubing@yeah.net。*通信作者，E-mail：xiaqy@swu.edu.cn

2016-05-06

2016-09-30