金保鋒，扈 強，葉為民，阮 志，李曉娜，陳富彩，勾 薇，崔志燕，4，張 玲，高 梅，張立新*

1.廣東中煙工業(yè)有限責任公司技術(shù)中心，廣州市天河區(qū)林和西橫路186號 510385

2.西北農(nóng)林科技大學生命科學學院，陜西省楊凌示范區(qū)西農(nóng)路22號 712100

3.西北農(nóng)林科技大學資源環(huán)境學院，陜西省楊凌示范區(qū)邰城路3號 712100

4.商洛市煙草公司鎮(zhèn)安縣分公司，陜西省商洛市鎮(zhèn)安縣宏化街17號 715000

秦巴不同生態(tài)區(qū)烤煙基因表達譜分析

金保鋒1，扈 強1，葉為民1，阮 志2，李曉娜3，陳富彩2，勾 薇2，崔志燕2，4，張 玲2，高 梅2，張立新*2

1.廣東中煙工業(yè)有限責任公司技術(shù)中心，廣州市天河區(qū)林和西橫路186號 510385

2.西北農(nóng)林科技大學生命科學學院，陜西省楊凌示范區(qū)西農(nóng)路22號 712100

3.西北農(nóng)林科技大學資源環(huán)境學院，陜西省楊凌示范區(qū)邰城路3號 712100

4.商洛市煙草公司鎮(zhèn)安縣分公司，陜西省商洛市鎮(zhèn)安縣宏化街17號 715000

為了揭示不同生態(tài)區(qū)烤煙風格特色形成的機理，采用Agilent煙草寡聚芯片技術(shù)，對商洛市洛南縣和寶雞市隴縣（對照）兩個生態(tài)區(qū)烤煙打頂后15 d同一品種（秦煙96）中部葉的基因表達譜進行了初步分析。結(jié)果表明：與隴縣相比，在44 000個基因中洛南煙區(qū)煙葉有2 058個基因差異表達，其中上調(diào)基因962個，占46.74%；下調(diào)基因1 096個，占53.26%。按功能富集分類，這些差異基因主要富集在代謝、響應(yīng)外界刺激及生理調(diào)控等生物學過程，大部分位于胞外區(qū)與細胞內(nèi)膜系統(tǒng)。其中洛南煙區(qū)碳水化合物合成代謝較為活躍，隴縣煙區(qū)萜類化合物代謝、苯丙烷代謝、芳香族化合物合成代謝較為活躍?；虼x通路分析表明，細胞信號轉(zhuǎn)導(dǎo)、轉(zhuǎn)錄因子及參與植物非生物脅迫相關(guān)基因在兩個煙草種植區(qū)均為高表達。

烤煙；生態(tài)區(qū)；秦煙96；基因；基因芯片；表達譜

基因芯片技術(shù)可用于高通量基因表達平行分析、大規(guī)?；虬l(fā)現(xiàn)及序列分析、基因多態(tài)性分析和基因組研究等［1］。目前，基因芯片技術(shù)已在水稻、玉米、小麥等［2-4］作物上廣泛應(yīng)用，而以煙草為材料的研究報道較少。楊天旭等［5］發(fā)現(xiàn)滲透脅迫48 h后，煙草根系上調(diào)表達基因主要是激素響應(yīng)蛋白基因、鈣信號響應(yīng)蛋白基因及蛋白激酶基因，而下調(diào)表達基因則主要是赤霉素響應(yīng)蛋白類基因和泛素連接酶類基因。利用基因芯片技術(shù)研究不同生態(tài)區(qū)煙草葉片基因表達情況僅有崔紅等［6］對烤煙移栽后60 d葉面腺毛表達譜進行過初步分析，而有關(guān)不同生態(tài)區(qū)烤煙打頂后煙葉基因表達譜變化的研究尚未見報道。因此，采用煙草寡聚芯片技術(shù)對兩個生態(tài)區(qū)烤煙打頂后煙葉基因表達變化進行試驗，旨在從分子水平上揭示生態(tài)因子對相同基因型烤煙物質(zhì)代謝的影響。

1 材料與方法

1.1 試驗材料

試驗于2013年進行，分別設(shè)置在寶雞市隴縣（LX，對照）曹家灣鎮(zhèn)流渠村（海拔1 055 m，N 34°53'57.3"，E 106°42'46.6"）和商洛市洛南縣（LN）城關(guān)鎮(zhèn)樊村（海拔 920 m，N 34°01'58''，E 110°18'05"）。供試材料為烤煙品種秦煙96。2013年洛南煙區(qū)大田期（5月—9月）平均氣溫12.51～25.56℃，月降水量為29.80～45.00 mm，大田期光照時數(shù)累計值為1 183.14 h；隴縣煙區(qū)大田期平均溫度為13.00～26.37℃，月降水量為65.20～259.40 mm，大田期光照時數(shù)累計值為1 220.55 h。兩生態(tài)煙區(qū)烤煙均于5月7日移栽，7月23日打頂。按當?shù)責熑~生產(chǎn)技術(shù)規(guī)范進行煙田管理。移栽后15 d選擇標記長勢長相均勻一致的煙株10株，分別在烤煙打頂后15 d（8月6日），用打孔器打下中部葉片第10葉位的葉片樣品，置于液氮速凍，帶回實驗室，置于-80℃冰箱中，保存?zhèn)溆谩?/p>

1.2 試驗方法

1.2.1 總RNA提取及純化

采用Trizol（美國Gibco BRL公司）提取煙草葉片總RNA，并進一步用RNA去除試劑盒（德國MACHEREY-NAGEL公司）對總RNA進行純化。

1.2.2 cRNA的合成及熒光標記

以總RNA起始，含有T7噬菌體RNA聚合酶啟動子序列的寡聚dT為引物，使用CbcScript酶合成cDNA第一條鏈。用核糖核酸酶H（RNaseH）將雜合鏈中的RNA切成隨機短片段，DNA聚合酶以RNA短片段作為引物進行延伸，合成與cDNA第一條鏈互補配對的第二條鏈，并純化雙鏈cDNA。以cDNA為模板，利用T7酶混合物合成cRNA，然后用RNA去除試劑盒純化。取5 μg cRNA，用CbcScriptⅡ酶利用隨機引物進行反轉(zhuǎn)錄，反轉(zhuǎn)錄產(chǎn)物用PCR NucleoSpinExtract II Kit（德 國 MACHEREYNAGEL公司）純化。取上述反轉(zhuǎn)錄產(chǎn)物，以隨機引物進行Klenow酶標記，之后用PCR NucleoSpin Extract II Kit純化經(jīng)Klenow酶標記的產(chǎn)物，純化后干燥。用熒光染料Cy5-dCTP和Cy3-dCTP（美國GE Healthcare公司，型號：No.PA 55021/PA53021）分別標記試驗樣品（一種為洛南樣品，另一種為隴縣樣品）。洛南與隴縣兩地樣品設(shè)為一組，組內(nèi)試驗采用熒光正反標的方法進行標記，每個樣品進行2組4次獨立雜交試驗，每組3次生物學重復(fù)。

1.2.3 芯片的檢測與數(shù)據(jù)分析

將標記的DNA溶于2×GEx Hyb Buffer+25%甲酰胺的雜交液中，45℃條件下雜交過夜。雜交結(jié)束后，在42℃左右含0.2%SDS+2×SSC的液體中漂洗5 min，之后置于25℃下，0.2×SSC中漂洗5 min。玻片干燥后用微陣列掃描儀（G2565CA型，美國Agilent公司）對芯片進行掃描，得到雜交圖。

采用Feature Extraction圖像分析軟件將芯片圖像信號轉(zhuǎn)為數(shù)字信號。運用批Percentile shift方法，將原始數(shù)據(jù)導(dǎo)入至GeneSpring GX軟件中，對數(shù)字信號值進行歸一化處理。以絕對差異倍數(shù)（Absolute Fold change）大于或等于2、同時用Flag數(shù)據(jù)標簽標記成可檢測的標準形式進行差異基因篩選。絕對差異倍數(shù)越大，說明兩個樣品之間的差異越大。

1.2.4 生物信息學分析

查閱NCBI數(shù)據(jù)庫對差異表達基因進行分析。采用AgriGO和WEGO等［7］網(wǎng)站對差異表達基因進行GO注釋和富集分析，并運用BLASTN對基因代謝途徑進行分析。

2 結(jié)果與分析

2.1 RNA提取與芯片雜交結(jié)果

經(jīng)甲醛變性膠電泳檢測，RNA樣品電泳條帶清晰，見圖1。28S∶18S rRNA條帶亮度大于或接近 2∶1，且 A260/A280≥1.80，RNA 總量≥8 μg（表 1），滿足了表達譜芯片實驗要求。表明RNA純度及完整性可滿足后續(xù)實驗要求。

圖1RNA電泳圖Fig.1 Electrophoretogram of RNA

表1 樣品質(zhì)檢結(jié)果Tab.1 Quality of samples

烤煙打頂后15 d，在44 000個芯片探針中，檢測到的雜交信號有34 080個，檢出率為77.45%。說明在兩個煙區(qū)樣品中共檢測到34 080個基因的表達量發(fā)生了變化。將兩組實驗的熒光數(shù)據(jù)校正、整合后制作散點圖，見圖2。從圖2可以看出，多數(shù)數(shù)據(jù)點集中在上下兩條斜線之間，說明絕大多數(shù)基因在兩個生態(tài)區(qū)中表達差異不大。

2.2 烤煙基因表達差異分析

圖3顯示，與隴縣相比，洛南煙區(qū)煙葉表達量差異在2倍以上的有2 058個，其中上調(diào)基因有962個，占46.74%；下調(diào)基因有1 096個，占53.26%。

2.3 烤煙差異表達基因富集分析

應(yīng)用AgriGO和WEGO數(shù)據(jù)庫進行基因本體注釋和功能富集分析。圖4顯示，烤煙打頂后15 d，與隴縣相比，在洛南煙區(qū)烤煙葉片2 058個差異表達基因中有763個（37.07%）基因有功能注釋。差異表達基因顯著富集在代謝、響應(yīng)刺激及生理調(diào)控等生物學過程中，定位在細胞膜、胞內(nèi)及細胞器中。同時這些差異表達基因翻譯的蛋白質(zhì)主要具有催化活性、分子轉(zhuǎn)運及轉(zhuǎn)錄調(diào)控活性。

圖2 洛南煙區(qū)烤煙煙葉差異表達基因散點圖Fig.2 Differentially expressed gene scatterplot of flue-cured tobacco leaves from Luonan

圖3 洛南煙區(qū)烤煙煙葉差異表達基因Fig.3 Differentially expressed genes of flue-curedtobacco leaves from Luonan

由表2可知，與隴縣相比，洛南煙區(qū)上調(diào)差異表達基因達到顯著水平的GO條目有10個［錯誤發(fā)現(xiàn)率（FDR）＜0.05］，其中生物學過程5個、分子功能4個、細胞組分1個，分別占50%、40%和10%。生物學過程主要為節(jié)律過程（GO:0048511，P）、晝夜節(jié)律（GO:0007623，P）、響應(yīng)二糖刺激（GO:0034285，P）、蔗糖刺激（GO:0009744，P）、果糖刺激（GO:0009750，P）。分子功能主要為原葉綠素酸酯還原酶活性（GO:0016630，F(xiàn)）、氧化還原酶活性（GO:0016628、GO:0016627，F(xiàn)）和水解酶活性（GO:0016798，F(xiàn)）。差異表達基因富集的細胞組分主要涉及胞外區(qū)（GO:0005576，C）。

圖4 洛南煙區(qū)烤煙煙葉差異表達基因富集分析Fig.4 Analysis of differentially expressed genes of flue-cured tobacco leaves from Luonan

表2 洛南煙區(qū)烤煙煙葉上調(diào)差異表達基因GO功能分類Tab.2 GO function classification of differentially expressed up-regulated genes of flue-cured tobacco from Luonan

由表3可知，與隴縣相比，洛南煙區(qū)下調(diào)差異表達基因達到顯著水平的GO條目有33個（FDR＜0.05），其中生物學過程19個、分子功能13個、細胞組分1個，分別占57.58%、39.39%和3.03%。生物學過程主要為次級代謝（GO:0019748，P）、倍半萜類代謝（GO:0016106、GO:0051762、GO:0051761、GO:0006714，P）、萜烯代謝（GO:0042214，P）、萜類化合物代謝（GO:0016114、GO:0042214，P）、營養(yǎng)生長向生殖生長轉(zhuǎn)變時間調(diào)控（GO:0048510，P）、分生組織相變時間調(diào)控（GO:0048506，P）、發(fā)育調(diào)控（GO:0040034，P）、類異戊二烯合成（GO:0008299、GO:0006720，P）、脂類合成（GO:0008610，P）、苯丙烷代謝（GO:0009698，P）、氨基酸衍生物合成（GO:0006575，P）、芳香族化合物合成（GO:0019438，P）過程。分子功能主要為α-蛇麻烯合成酶活性（GO:0080017，F(xiàn)）、反式-E-β-石竹烯合成酶活性（GO:0080016，F(xiàn)）、碳-氧裂解酶活性（GO:0016838、GO:0016835，F(xiàn)）、腺苷酸環(huán)化酶活性（GO:0009975，F(xiàn)）、催化活性（GO:0003824，F(xiàn)）、主動跨膜轉(zhuǎn)運蛋白活性（GO:0015291、GO:0022804，F(xiàn)）、裂解酶（GO:0016829，F(xiàn)）、陽離子轉(zhuǎn)運蛋白活性（GO:0015294，F(xiàn)）、轉(zhuǎn)運蛋白活性（GO:0015293，F(xiàn)）及轉(zhuǎn)移酶活性（GO:0016765、GO:0016740，F(xiàn)）。細胞組分主要位于細胞內(nèi)膜系統(tǒng)（GO:0012505，C）。

表3 洛南煙區(qū)烤煙煙葉下調(diào)差異表達基因GO功能分類Tab.3 GO function classification of differentially expressed down-regulated genes of flue-cured tobacco leaves from Luonan

2.4 烤煙差異表達基因代謝途徑分析

由表4可知，與隴縣相比，在洛南煙區(qū)煙草葉片中差異表達20倍以上的上調(diào)基因共有152個，有13個基因功能已知，其余均為未知功能蛋白。

已知的上調(diào)表達功能基因主要分為4類：①光合色素相關(guān)基因。原葉綠素酸酯還原酶在開花植物依賴于光的顯綠過程中起重要作用［8］。捕光色素葉綠素a/b結(jié)合蛋白為光合系統(tǒng)蛋白，是與色素形成的色素蛋白復(fù)合體（LHC），LHC能將捕獲的光能迅速傳至反應(yīng)中心并引起光化學反應(yīng)，在光合作用和對各種環(huán)境的適應(yīng)過程中發(fā)揮作用［9］。②細胞結(jié)構(gòu)形成相關(guān)基因。果膠甲基酯酶抑制蛋白，果膠甲基酯酶（PME）是一種果膠酶，在細胞壁降解過程中起著重要作用［10］。堿性α-半乳糖苷酶在植物種子、幼苗、葉片及果實中均存在，參與種子萌發(fā)、葉片發(fā)育、果實成熟及低溫耐受等生理過程［11-12］。類萌發(fā)素蛋白是存在于植物體中的重要防御酶，既能調(diào)控植物特定的生長發(fā)育過程，又能在植物響應(yīng)生物及非生物脅迫中發(fā)揮作用［13］。過氧化物酶前體在細胞器過氧化物酶體中被催化生成過氧化物酶，過氧化物酶與植物的抗逆性有關(guān)，是植物保護酶系的重要保護酶之一［14］。富含甘氨酸蛋白互作蛋白和富含脯氨酸蛋白（PRP）是一類細胞壁結(jié)構(gòu)蛋白，在細胞壁的形成和細胞的防衛(wèi)過程中有重要作用，其表達具有組織特異性并受發(fā)育階段特異性及多種環(huán)境因素的調(diào)控［15］。③細胞信號轉(zhuǎn)導(dǎo)相關(guān)基因。絲氨酸/蘇氨酸蛋白激酶SAPK8類蛋白是一大類特異性催化蛋白質(zhì)絲氨酸和蘇氨酸殘基磷酸化的激酶家族，參與細胞中許多信號轉(zhuǎn)導(dǎo)途徑［16］。含NB-ARC結(jié)構(gòu)域蛋白在植物超敏反應(yīng)中發(fā)揮作用［17］。④轉(zhuǎn)錄因子類：富含甘氨酸結(jié)構(gòu)域RNA結(jié)合蛋白通過與RNA相互作用來調(diào)節(jié)細胞的功能，從而適應(yīng)多種逆境條件［18］。Drm3類蛋白、生長素抑制蛋白均參與生長素誘導(dǎo)基因的下調(diào)表達調(diào)控［19］。

表4 洛南煙區(qū)烤煙煙葉差異表達20倍以上的基因功能注釋Tab.4 Function of genes differentially expressed more than 20 times of flue-cured tobacco from Luonan

下調(diào)差異表達20倍以上的基因有86個，其中有8個基因功能已知，其余均為未知功能蛋白。半胱氨酸蛋白酶參與脅迫響應(yīng)、組織分化、衰老等過程［20］。MYB類轉(zhuǎn)錄因子廣泛參與植物的生長發(fā)育和代謝調(diào)控過程［21］。RNA酶H家族蛋白，該酶具有核酸酶活性，催化DNA-RNA雜合體的RNA部分的核內(nèi)降解，對細胞的增殖和分化起調(diào)控作用［22］。馬兜鈴烯合成酶是一種萜烯合成酶，能催化不同類異戊二烯中間產(chǎn)物使其環(huán)化［23］。肉桂酰輔酶A還原酶（CCR）是木質(zhì)素特異途徑的第一個關(guān)鍵酶，對木質(zhì)素單體的生物合成起重要作用，而木質(zhì)素對植物的結(jié)構(gòu)與防御功能都具有重要作用［24-25］。Patatin 類蛋白1（PLP1）為一類酯酰基水解酶，參與植物生物脅迫與非生物脅迫過程［26-27］。生長素類蛋白是一類參與植物生長素產(chǎn)生和降解的重要蛋白［28］。肌醇-3-磷酸合成酶是植酸合成代謝中催化葡糖糖轉(zhuǎn)化為肌醇-3-磷酸的功能蛋白，植酸可以防止植物葉片褐變，提高植物的抗氧化性［29-30］。

3 結(jié)論與討論

在本試驗中通過對兩生態(tài)煙區(qū)烤煙品種秦煙96在打頂后15 d煙葉基因表達譜分析顯示，洛南生態(tài)煙區(qū)差異表達基因為2 058個，其中上調(diào)差異表達基因為962個，下調(diào)差異表達基因為1 096個。這些差異基因主要在植物生長發(fā)育、細胞信號轉(zhuǎn)導(dǎo)、物質(zhì)代謝與運輸、抵御逆境等過程中發(fā)揮作用。葉綠素合成、碳水化合物代謝與萜類化合物代謝相關(guān)基因的差異表達可能是造成兩個生態(tài)區(qū)煙葉風格特色差異的內(nèi)因。

適宜的溫度能滿足烤煙生長發(fā)育，有利于物質(zhì)的積累和煙葉充分成熟，同時適宜的光照可以保證碳水化合物的合成積累及煙葉的色澤［31］。葉綠素在高溫和強光照下易發(fā)生降解，煙葉成熟期溫度過高，會破壞葉綠素，影響光合作用，從而使植物新陳代謝失調(diào)，明顯地影響煙株的生長、成熟和煙葉品質(zhì)［32-33］。洛南煙區(qū)在烤煙打頂后煙葉葉綠素合成相關(guān)基因，如光合色素原葉綠素酸酯還原酶、捕光色素葉綠素a/b結(jié)合蛋白前體相關(guān)基因高度表達。洛南煙區(qū)的溫度有利于煙葉中葉綠素生物合成，從而促進了光合作用。整個生育期的平均氣溫與煙葉總糖和還原糖含量密切相關(guān)，平均氣溫適宜時總糖和還原糖含量較高，過高則不利于總糖和還原糖合成［34］。推測洛南煙區(qū)烤煙成熟前期溫度較為適宜，促進了碳代謝化合物相關(guān)基因高表達，提高了碳水化合物累積量，進而為后期次生代謝提供了充足的前體物質(zhì)和能量，為萜類化合物合成創(chuàng)造了良好的前提條件。不同氣候條件如降雨量和降雨分布會影響土壤水分狀況與空氣濕度，導(dǎo)致煙草葉面腺毛分泌物發(fā)生變化，進而對煙葉香氣品質(zhì)產(chǎn)生影響［35］。隴縣煙區(qū)煙葉萜類化合物、芳香族化合物基因表達量較高，如馬兜鈴烯合成酶基因。這可能是由于洛南煙區(qū)8月上旬平均氣溫與降雨量要高于隴縣煙區(qū)所致［36］?？梢姡鷳B(tài)因子對兩個生態(tài)煙區(qū)葉綠素合成、碳代謝和萜類代謝相關(guān)基因表達產(chǎn)生較大影響，導(dǎo)致兩地煙葉化學成分的差異。張松濤等［37］研究提出生態(tài)條件差異對細胞內(nèi)核苷酸代謝、激素代謝、DNA復(fù)制和修復(fù)以及信號轉(zhuǎn)導(dǎo)相關(guān)基因的表達有一定影響。在本研究中也有類似發(fā)現(xiàn)，兩個生態(tài)區(qū)差異表達基因均涉及細胞結(jié)構(gòu)（細胞壁）的形成、細胞信號轉(zhuǎn)導(dǎo)及參與非生物脅迫等過程。表明不同生態(tài)區(qū)煙葉風格特色的差異是由多基因參與調(diào)控的結(jié)果。

生態(tài)因子影響著細胞內(nèi)眾多基因的表達，造成這些基因表達差異的原因尚不完全明晰。另外相關(guān)數(shù)據(jù)庫中基因的信息及功能注釋不太全面，也影響了數(shù)據(jù)分析的深度。此外，烤煙葉片基因差異表達還需通過熒光定量PCR進一步驗證。在差異表達20倍以上的基因中，還有一部分并未進行功能注釋，它們的差異表達對于兩個生態(tài)煙區(qū)烤煙風格特色形成的影響尚不清楚，需要進行更深入的研究。

［1］ 滕曉坤，肖華勝.基因芯片與高通量DNA測序技術(shù)前景分析［J］.中國科學（C輯:生命科學），2008，38（10）：891-899.

［2］ 肖軍，楊濤，楊鎮(zhèn)，等.應(yīng)用基因芯片分析植物內(nèi)生菌醇提取物處理條件下水稻的基因表達圖譜［J］.生物技術(shù)通報，2013，29（1）：89-95.

［3］ 袁廣勝，杜娟，高健，等.玉米穗粒腐病差異表達基因的生物信息學分析［J］.玉米科學，2013，21（5）：46-51.

［4］ 董劍，劉迎團，劉宏偉，等.基因芯片分析殺雄劑SQ-1誘導(dǎo)小麥生理型雄性不育的基因差異表達［J］.中國生物化學與分子生物學報，2013，29（10）：948-955.

［5］ 楊天旭，汪耀富，劉國順，等.滲透脅迫下煙草葉片基因的差異表達研究［J］.作物學報，2007，33（6）：914-920.

［6］ 崔紅，冀浩，楊惠絹，等.不同生態(tài)區(qū)煙草的葉面腺毛基因表達［J］.生態(tài)學報，2011，31（11）：3143-3149.

［7］ Ye J，F(xiàn)ang L，Zheng H，et al.WEGO：a web tool for plotting GO annotations［J］.Nucleic Acids Research，2006，34（2）：293-297.

［8］ Lebedev N，Timko M P.Protochlorophyllide oxidoreductase B-catalyzed protochlorophyllide photoreduction in vitro：insight into the mechanism of chlorophyll formation in light-adapted plants［J］.Proceedings of the National Academy of Sciences of the United States of America，1999，96（17）：9954-9959.

［9］ 未麗，徐秉良，雷江麗，等.胡楊葉綠素a/b結(jié)合蛋白基因的克隆及序列特性分析［J］.中國農(nóng)業(yè)科技導(dǎo)報，2008，10（4）：63-69.

［10］董章勇，王振中.香蕉枯萎病菌果膠甲基酯酶基因的克隆與生物信息學分析［J］.生物技術(shù)通報，2012，28（6）：111-115.

［11］ Marraccini P，Rogers W J，Caillet V，et al.Biochemical and molecular characterization of alpha-D-galactosidase from coffee beans［J］. Plant Physiology and Biochemistry，2005，43（43）：909-920.

［12］ PetersS， Mundree S G，Thomson JA，etal.Protection mechanisms in the resurrection plant Xerophyta viscosa（Baker）：both sucrose and raffinose family oligosaccharides（RFOs） accumulate in leaves in response to water deficit［J］. Journal of Experimental Botany，2007，58（8）：1947-1956.

［13］王沙沙，舒文濤，陳軍營，等.麥類作物萌發(fā)素及其相關(guān)蛋白研究進展［J］.河南農(nóng)業(yè)大學學報，2007，41（2）：242-246.

［14］蔣選利，李振岐，康振生.過氧化物酶與植物抗病性研究進展［J］.西北農(nóng)林科技大學學報：自然科學版，2001，29（6）：124-129.

［15］柴團耀，張玉秀.菜豆富含脯氨酸蛋白質(zhì)基因在生物和非生物脅迫下的表達［J］.植物學報，1999，41（1）：111-113.

［16］臧夢維.絲氨酸/蘇氨酸蛋白激酶［J］.基礎(chǔ)醫(yī)學與臨床，1996，16（2）：86-95.

［17］王文娟，張飛云.植物抗病分子機制研究進展［J］.生物技術(shù)通報，2007，23（1）：19-23.

［18］ Kim Y O，Kim J S，Kang H.Cold-inducible zinc finger-containing glycine-rich RNA-binding protein contributes to the enhancement of freezing tolerance in Arabidopsis thaliana［J］.The Plant Journal，2005，42（6）：890-900.

［19］Misra P，Pandey A，Tiwari M，et al.Modulation of transcriptome and metabolome of tobacco by Arabidopsistranscription factor，AtMYB12，leads to insect resistance［J］.Plant Physiology，2010，152（4）：2258-2268.

［20］閆龍鳳，楊青川，韓建國，等.植物半胱氨酸蛋白酶研究進展［J］.草業(yè)學報，2005，14（5）：11-19.

［21］ Martin C.Multifunctionality and diversity within the plant MYB-gene family［J］.Plant Molecular Biology，1999，41（5）：577-585.

［22］ SchultzS J，Champoux JJ.RNaseH activity：structure， specificity， and function in reverse transcription［J］.Virus Research，2008，134（1）：86-103.

［23］郝宏蕾，朱旭芬，曾云中.類異戊二烯的生物合成及調(diào)控［J］.浙江大學學報：農(nóng)業(yè)與生命科學版，2002，8（2）：108-114.

［24］李波，梁穎，柴友榮.植物肉桂酰輔酶A還原酶（CCR）基因的研究進展［J］.分子植物育種，2006，4（S1）：55-65.

［25］章霄云，郭安平，賀立卡，等.木質(zhì)素生物合成及其基因調(diào)控的研究進展［J］.分子植物育種，2006，4（3）：431-437.

［26］Camera S L，Balagué C，G?bel C，et al.The Arabidopsispatatin-like protein 2 （PLP2） plays an essential role in cell death execution and differentially affects biosynthesis of oxylipins and resistance to pathogens［J］. Molecular Plant-microbe Interactions，2009，22（4）：469-481.

［27］Matos A R，D'Arcy-Lameta A.A novel patatin-like gene stimulated by drought stress encodes a galactolipid acyl hydrolase［J］.Febs Letters，2001，491（3）：188-192.

［28］DreherK A，BrownJ，Saw R E，etal.The Arabidopsis Aux/IAA protein family has diversified in degradation and auxin responsiveness［J］.The Plant Cell，2006，18（3）：699-714.

［29］袁鳳杰，董德坤，李百權(quán)，等.兩個低植酸大豆突變體中肌醇-3-磷酸合成酶MIPS1基因表達特性的研究［J］.核農(nóng)學報，2013，27（03）：286-292.

［30］李偉，高玉珍，劉霞，等.植酸對煙葉破碎后幾個生理指標的影響［J］.鄭州輕工業(yè)學院學報：自然科學版，2007，22（4）：9-12.

［31］王玉璽，栗珂，韋成才，等.陜南優(yōu)質(zhì)烤煙氣候條件及區(qū)域劃分的研究［J］.陜西氣象，2001，42（5）：15-18.

［32］張武昌，王榮.光和溫度對葉綠素a和脫鎂葉綠酸a降解的影響［J］.海洋科學，2000，24（4）：50-52.

［33］李琦.烤煙優(yōu)劣質(zhì)年的氣象條件分析［J］.安徽農(nóng)業(yè)科學，1997，25（2）：32-35.

［34］韋成才，馬英明，艾綏龍，等.陜南烤煙質(zhì)量與氣候關(guān)系研究［J］.中國煙草科學，2004，25（3）：38-41.

［35］楊克偉.生態(tài)因子對煙葉香型風格劃分的影響［J］.北京農(nóng)業(yè)，2013，32（3）：19-20.

［36］郭兆夏，賀文麗，李星敏，等.基于GIS的陜西省烤煙氣候生態(tài)適宜性區(qū)劃［J］.中國煙草學報，2012，18（2）：21-24.

［37］張松濤，楊永霞，楊惠娟，等.不同生態(tài)區(qū)煙草葉片基因表達譜的分析［J］.農(nóng)業(yè)生物技術(shù)學報，2014，22（10）：1251-1260.

Analysis of Gene Expression Profiles of Flue-cured Tobaccos from Different Ecological Zones in Qinba Mountain Area

JIN Baofeng1,HU Qiang1,YE Weimin1,RUAN Zhi2,LI Xiaona3,CHEN Fucai2,GOU Wei2,CUI Zhiyan2,4,ZHANG Ling2,GAO Mei2,and ZHANG Lixin*2
1.Technology Center,China Tobacco Guangdong Industrial Co.,Ltd.,Guangzhou 510385,China
2.College of Life Sciences,Northwest A&F University,Yangling 712100,Shaanxi,China
3.College of Natural Resources&Environment,Northwest A&F University,Yangling 712100,Shaanxi,China
4.Zhen'an County Tobacco Company of Shaanxi Province,Zhen'an 715000,Shaanxi,China

To investigate the formation mechanism of style characteristics of flue-cured tobaccos from different ecological environment conditions,the middle leaves of flue-cured tobacco cv.Qinyan96 grown in two ecological zones,i.e.Luonan County in Shangluo City and Longxian County in Baoji City（the control）were sampled 15 days after topping,and their gene expression profiles were preliminarily analyzed with Agilent tobacco oligo microarray.The results showed that,comparing with those from Longxian,of the 44 000 genes of leaves from Luonan,there were 2 058 differentially expressed genes,including 962（46.74%）up-regulated genes and 1 096（53.26%）down-regulated genes.Categorizing byfunction,those genes mainly involved in biological processes,including metabolism,response to external environment stimulation,physiological regulation;and mostly distributed in extracellular regions and endomembrane systems.Synthesis metabolism of carbohydrate was more active in leaves from Luonan,while that of terpenoid,phenylpropanoid and aromatic compound was active in leaves from Longxian.Metabolic pathway analysis indicated that genes related to cell signal transduction,transcription factor and involved in plant abiotic stress were of higher expression in leaves from both Luonan and Longxian.

Flue-cured tobacco;Ecological zone;Qinyan96;Gene;Microarray;Expression profile

S572.01

A

1002-0861（2015）12-0001-08

10.16135/j.issn1002-0861.20151201

2015-04-17

2015-06-15

廣東中煙工業(yè)有限責任公司科技計劃項目“‘雙喜’品牌秦巴生態(tài)煙葉關(guān)鍵技術(shù)研究與集成推廣”［粵煙工05XM-QK（2014）015］。

金保鋒（1973—），碩士，高級工程師，主要從事特色煙葉開發(fā)與應(yīng)用研究。E-mail：jinbf@gdzygy.com；*

張立新，E-mail：zhanglixin@nwsuaf.edu.cn

金保鋒，扈強，葉為民，等.秦巴不同生態(tài)區(qū)烤煙基因表達譜分析［J］.煙草科技，2015，48（12）：1-8.JIN Baofeng,HU Qiang,YE Weimin,et al.Analysis of gene expression profiles of flue-cured tobaccos from different ecological zones in Qinba mountain area［J］.Tobacco Science&Technology,2015,48（12）：1-8.

責任編輯 董志堅