草烏幼苗形成期葉綠體發(fā)育的轉(zhuǎn)錄組分析

摘要：為確定草烏（Aconitum kusnezoffii Reichb.）苗期生長發(fā)育過程中的轉(zhuǎn)錄組特征，本研究利用Illumina測序?qū)Σ轂鯊姆N子萌發(fā)到幼苗建成涉及到的四個時期（T1，T2，T3，T4）進行轉(zhuǎn)錄分析，挖掘草烏幼苗形成期葉綠體發(fā)育關(guān)鍵基因。結(jié)果顯示，4個時期的6個差異分類檢測到的差異基因數(shù)量依次為：16 379，28 155，16 100，34 368，38 250和44 568個。GO富集分析顯示T3vsT4差異表達基因主要富集到對幾丁質(zhì)的反應(yīng)，光合系統(tǒng)Ⅰ等。KEGG富集分析顯示差異表達基因顯著富集的通路主要是光合生物中的碳固定，植物激素信號轉(zhuǎn)導(dǎo)等。草烏幼苗形成期葉綠體發(fā)育過程中，19個差異表達基因注釋光信號調(diào)控葉綠體發(fā)育通路上，如GLK1/2，PIFs，HY5基因等。44個差異表達基因注釋植物激素信號轉(zhuǎn)導(dǎo)通路上，如AUX/IAA，ARF，DELLA基因等。本研究挖掘到一部分草烏幼苗形成期葉綠體發(fā)育相關(guān)基因，為進一步研究草烏生長發(fā)育相關(guān)分子機制提供理論基礎(chǔ)。

關(guān)鍵詞：草烏；轉(zhuǎn)錄組；幼苗形成；葉綠體發(fā)育

中圖分類號：Q943.2""" 文獻標識碼：A"""" 文章編號：1007-0435（2024）06-1779-10

Transcriptome Analysis of Chloroplast Development during the Formation

Stage of Aconitum kusnezoffii Seedlings

JIA Chun-he1， LYU Li-juan2， LI Xiao-jie1， ZHANG Xiao-ming1， MU Ying-tong1， FAN Dong-chang1，

LYU Yan-min3， SHI Jia-yi1， WANG Jun-jie1*

（1.Key Laboratory of Grassland Resources， College of Grassland and Resources Environment， Inner Mongolia Agricultural

University， Ministry of Education， Hohhot， Inner Mongolia 010019; 2.Hohhot Forestry and Grassland Bureau， Hohhot， Inner

Mongolia 010019; 3.Horqin Right Wing Middle Banner forestry and grassland Bureau， Horqin Right Wing Middle， Inner Mongolia 029400）

Abstract：In order to obtain the transcriptome characteristics of Aconitum kusnezoffii Reichb. during seedling growth and development in this study，Illumina sequencing was used to transcriptionally analyse key genes for chloroplast during the four periods （T1，T2，T3，T4） involved in the development of Aconitum kusnezoffii from seed germination to seedling establishment. The results showed that the number of differentially expressed genes detected in the six differential classifications of the four periods was 16 379，28 155，16 100，34 368，38 250 and 44 568，respectively. GO enrichment analysis showed that the differentially expressed genes in T3vsT4 were enriched in responses to chitin，photosynthetic system I，etc. KEGG enrichment analysis showed that the pathways significantly enriched in differentially expressed genes were mainly carbon fixation in photosynthetic organisms，plant hormone signal transduction，etc. For the development of chloroplasts during the formative stage of Aconitum kusnezoffii seedlings，nineteen differentially expressed genes，such as GLK1/2，PIFs， HY5 genes，are annotated to the light signaling pathway regulating chloroplast development. Forty-four differentially expressed genes annotated were in plant hormone signaling pathways，such as AUX/IAA，ARF，DELLA genes，etc. This study identified some genes related to chloroplast development during the formation stage of Aconitum kusnezoffii seedlings，providing a theoretical basis for further research on the molecular mechanisms related to the growth and development of Aconitum kusnezoffii.

Key words：Aconitum kusnezoffii Reichb;Transcriptome;Seedling formation;Chloroplast development

草烏（Aconitum kusnezoffii Reichb.）為毛茛科的多年生草本，別名北烏頭，廣泛分布于我國東北、華北地區(qū)。傳統(tǒng)蒙醫(yī)藥中，草烏塊根、葉、花及芽等部位分開入藥［1］。附子是北烏頭的干燥塊根，可入藥，具有祛風除濕、溫經(jīng)止痛之功效［2］。草烏的干燥葉、幼芽、花等用于治療風寒濕痹、關(guān)節(jié)疼痛、心腹冷痛及麻醉止痛等［3］。

光是影響植物生長和發(fā)育的關(guān)鍵環(huán)境因子之一，光作為能量來源為植物生長發(fā)育提供營養(yǎng)，還作為信號參與植物生命周期的調(diào)控，在植物幼苗形成過程中經(jīng)歷了光形態(tài)建成及葉綠體發(fā)育。光參與調(diào)控植物細胞的生長、分化，使植物細胞結(jié)構(gòu)、功能發(fā)生變化，匯聚成植物組織和器官，這個過程稱為植物的光形態(tài)建成［4］。植物的光形態(tài)構(gòu)建過程可人為按照最有效的光譜范圍分成兩種：紅光反應(yīng)和藍光反應(yīng)。植物葉子的增大主要靠紅光反應(yīng)，而氣孔打開、葉綠體分化與運動是藍光反應(yīng)［5］。在感知紅光時，光敏色素通過觸發(fā)磷酸化和蛋白酶體降解來抑制（Phytochrome-interacting factors，PIFs）的功能，PIFs是一類可以與光敏色素相互作用的螺旋-環(huán)-螺旋（helix-loop-helix）家族轉(zhuǎn)錄因子，進而減輕PIFs對葉綠素生物合成和光合系統(tǒng)相關(guān)基因的抑制［6］。此外，光敏色素還抑制光信號傳導(dǎo)的另一種抑制因子，即E3泛素連接酶復(fù)合物COP1-SPA1，從而驅(qū)動PhANGs的轉(zhuǎn)錄［7-8］。葉綠體正常發(fā)育的重要條件是光，一系列能夠感受特定波長的蛋白可感知光，蛋白通過構(gòu)象變化會與下游信號分子互相作用。光敏色素（Phytochrome A/B，phyA/B）、隱花色素（Cryptochrome，cry）被認為是調(diào)節(jié)光形態(tài)建成的兩類光受體［9］。當幼苗從土壤覆蓋中出現(xiàn)時，光引發(fā)了從暗生長的光形態(tài)建成到光生長的光形態(tài)建成的形態(tài)和生理巨變［10-11］，包括下胚軸伸長抑制、頂端彎鉤展開、子葉張開和膨大以及黃化體到葉綠體的轉(zhuǎn)變［12-14］。在光下PhyB能夠轉(zhuǎn)到核內(nèi)并與PIF3結(jié)合，以此自身磷酸化開啟，調(diào)控光形態(tài)建成相關(guān)基因的轉(zhuǎn)錄［9，15］，PIF3負調(diào)節(jié)編碼葉綠素合成調(diào)節(jié)酶基因HEMA1、GUN5以及光合系統(tǒng)PSI中LHCA1、PsaE1基因的表達［16-17］。光形態(tài)建成受PIFs家族基因抑制，尤其是葉綠體的發(fā)育［18］。葉綠體的生物發(fā)生通過大量信號轉(zhuǎn)導(dǎo)在細胞核和葉綠體之間進行協(xié)調(diào)，以此來促進功能葉綠體的組裝和生物合成［19］。

激素在調(diào)控葉綠體發(fā)育過程中同樣起著重要的作用。PIFs作為光信號調(diào)控葉綠體發(fā)育的重要調(diào)控因子，BZR1與PIF4的N端結(jié)構(gòu)域結(jié)合形成異源二聚體，從而共同抑制GLK1、GLK2和幾個葉綠素生物合成基因的表達。HY5與去磷酸化狀態(tài)BZR1形成異源二聚體，減弱了BZR1的活性，進而促進了葉綠體發(fā)育相關(guān)基因的表達，使葉綠體的發(fā)育不再受到抑制。DELLA蛋白是GA信號轉(zhuǎn)導(dǎo)的抑制子，當感受到GA信號時，DELLA蛋白經(jīng)26S蛋白酶體降解失活。當DELLA蛋白不響應(yīng)GA信號時，DELLA蛋白結(jié)合PIF并通過泛素-蛋白酶體途徑誘導(dǎo)其降解，以此降低PIF活性，進而削弱了PIF對葉綠體發(fā)育相關(guān)基因的轉(zhuǎn)錄抑制［20］。

目前草烏的研究多集中在栽培技術(shù)和生物堿含量等方面，關(guān)于草烏葉綠體發(fā)育方面的研究未見報道。葉綠體不光是植物細胞的重要器官，還是維持植物生命活動的關(guān)鍵因子，植物進行光合作用的主要場所是葉綠體，而光合作用的進行依賴葉綠體的正常發(fā)育，因此，葉綠體的正常發(fā)育對草烏生長發(fā)育具有重要意義。本文通過對草烏種子萌發(fā)及幼苗建成涉及到的四個時期進行轉(zhuǎn)錄測序，挖掘草烏幼苗形成期葉綠體發(fā)育相關(guān)基因，為進一步研究草烏生長發(fā)育調(diào)控過程以及烏頭屬植物進化提供依據(jù)。

1 材料與方法

1.1 材料來源

草烏（Aconitum kusnezoffii Reichb.）種子于2020年在內(nèi)蒙古自治區(qū)赤峰市巴林右旗賽罕烏拉自然保護區(qū)采集。選取顆粒飽滿、大小一致的種子用清水浸泡48 h，種子與砂土按1∶3拌勻，置于容器中封好，砂土濕度（15%～20%），置于-5℃/4℃變溫環(huán)境下每 7 d交替一次，相對濕度 50%低溫培養(yǎng)箱中層積。第0 d（CK）、14 d（種子吸脹）、42 d（種子露白）取樣，待種子露白后，放置于15℃光照培養(yǎng)箱中催芽處理，14 d后取樣（種子子葉露出，并且變?yōu)榫G色），各時期稱取0.5 g新鮮樣品，3次重復(fù)，分別用T1，T2，T3和T4表示。

1.2 RNA提取和轉(zhuǎn)錄組測序

采用Trizol法提取樣品總RNA［21］，得到草烏總RNA后用Agilent 2100檢測RNA的質(zhì)量，質(zhì)量合格后一部分-70℃保存用于qRT-PCR，一部分送往上海歐易生物科技有限公司進行上機測序。原始數(shù)據(jù)（raw data）用Fastqc質(zhì)量評估后用trim-galore對原始數(shù)據(jù)進行剪切得到clean data。使用Trinity［22］軟件采用paired-end方法對clean data進行無參轉(zhuǎn)錄組從頭組裝（de novo）得到轉(zhuǎn)錄本，選擇最長的轉(zhuǎn)錄本作為Unigene，用CD-HIT［23］軟件去冗余得到最終的Unigene，作為后續(xù)分析的參考序列。

1.3 草烏的unigenes的功能注釋和葉綠體發(fā)育相關(guān)基因的鑒定

為獲得草烏的unigene功能信息，將測序得到的Unigene比對Nr、Nt、KOG/COG、Swiss-Prot、KO和GO數(shù)據(jù)庫［24］，使用HMMER［25］軟件對Pfam數(shù)據(jù)庫進行Unigene的功能分析。利用FPKM法計算基因的表達水平，利用DESeq2軟件分析T3vsT4時期差異表達，篩選條件為qlt;0.05 且差異倍數(shù)（foldChange）gt; 2為差異表達基因。

1.4 實時熒光定量qRT-PCR驗證

在本研究所有差異表達基因中隨機選取9個與草烏葉綠體發(fā)育相關(guān)基因進行驗證，使用Primer3plus網(wǎng)站設(shè)計PCR引物（表1），以草烏中actin為內(nèi)參基因，采用Trizol法提取樣品總RNA，使用PrimeScript TM 1 st Strand cDNA Synthesis Kit（TaKaRa）反轉(zhuǎn)錄試劑盒進行cDNA合成，以cDNA為模板進行qRT-PCR分析。qRT-PCR反應(yīng)體系見表2，試驗技術(shù)設(shè)置3次重復(fù)，基因表達量采用2-ΔΔCT法計算［26］。

2 結(jié)果與分析

2.1 草烏轉(zhuǎn)錄組數(shù)據(jù)組裝與功能注釋

通過用Trinity軟件對所有reads進行組裝，我們在草烏中發(fā)現(xiàn)了79 251個unigene，平均長度為881.74 bp，組裝的unigene的N50長度為1 221 bp，平均GC含量為46.72%。為了獲得本研究中草烏葉綠體的基因功能信息，我們采用diamon軟件將Unigene與6大數(shù)據(jù)庫進行比對，選擇具有最高序列相似性的蛋白，獲得功能注釋信息（表3）。

2.2 葉綠體發(fā)育過程差異表達基因及其功能富集分析

2.2.1 差異表達基因的統(tǒng)計分析 以qlt;0.05 且差異倍數(shù)（fold Change）gt; 2作為篩選條件，將4個不同時期草烏兩兩比對，即T1vsT2、T1vsT3、T1vsT4、T2vsT3、T2vsT4、T3vsT4篩選差異表達基因（DEGs），結(jié)果表明，在T1vsT2中獲得16 100個差異表達基因，上調(diào)基因數(shù)13 119個，下調(diào)基因數(shù)2 981個。在T1vsT3中獲得 38 250個差異表達基因，上調(diào)基因數(shù)20 985個，下調(diào)基因數(shù)17 265個。在T1vsT4中獲得44 568個差異表達基因，上調(diào)基因數(shù)23 594個，下調(diào)基因數(shù)20 974個。在T2vsT3中獲得28 155個差異表達基因，上調(diào)基因數(shù)11 728個，下調(diào)基因數(shù)16 427個。在T2vsT4中獲得34 268個差異表達基因，上調(diào)基因數(shù)13 491個，下調(diào)基因數(shù)20 877個。在T3vsT4中獲得16 379個差異表達基因，上調(diào)基因數(shù)6 433個，下調(diào)基因數(shù)9 946個（圖1）。

2.2.2 差異表達基因的GO富集分析 為了研究草烏幼苗形成期葉綠體發(fā)育過程中差異表達基因的主要生物學功能，對T3vsT4差異表達基因進行GO功能分析。圖2展示了T3vsT4對比組合在生物過程、細胞組分、分子功能中顯著性前10的GO條目，從圖中可以看出差異基因富集在生物過程和細胞組分較多。T3vsT4差異基因在生物過程中多富集在對幾丁質(zhì)的反應(yīng)、防御反應(yīng)過程等；在細胞組分中差異基因主要富集在光合系統(tǒng)Ⅰ、細胞壁等；在分子功能多富集在脂肪酶活性、鄰苯二酚氧化酶活性等。對所有差異基因進行GO富集分析，結(jié)果顯示草烏在光形態(tài)建成和葉綠體發(fā)育過程中在生物過程富集到參與光合作用（GO：0015979）、防御反應(yīng)（GO：0006952）。在分子功能富集到DNA結(jié)合轉(zhuǎn)錄因子活性（GO：0003700）。在細胞組分富集到葉綠體類囊體膜（GO：0009535）、葉綠體（GO：0009507）、類囊體（GO：0009579）、葉綠體基質(zhì)（GO：0009570）、葉綠體類囊體（GO：0009534）、植物型細胞壁（GO：0009505）、質(zhì)外體（GO：004804）。以上結(jié)果表明這些差異表達基因主要參與光合作用和葉綠體發(fā)育相關(guān)代謝途徑，參與調(diào)控葉綠體的發(fā)育，影響光合作用（表4）。

2.2.3 差異表達基因的KEGG富集分析 通過對T3，T4實驗組間比較的差異表達基因進行KEGG代謝通路顯著性分析，選擇前20條（篩選對應(yīng)差異基因數(shù)目大于 2的Pathway條目，按照每個條目對應(yīng)的-lgPvalue由大到小排序）通路做氣泡圖進行分析（圖3），T3vsT4時期差異表達基因顯著富集的通路主要是核糖體、植物-病原體相互作用、植物激素信號轉(zhuǎn)導(dǎo)、光合生物中的碳固定、苯丙類生物合成等。核糖體、植物-病原體相互作用、植物激素信號轉(zhuǎn)導(dǎo)、光合生物中的碳固定等代謝通路可能在草烏幼苗形成期葉綠體發(fā)育過程中起到一定作用。

2.3 光信號調(diào)控草烏葉綠體發(fā)育過程中相關(guān)基因表達模式分析

在光對葉綠體發(fā)育的影響途徑中，本研究鑒定出了1個編碼GOLDEN2-LIKE（GLK）蛋白的DEGs，是促進葉綠體發(fā)育和葉綠素合成的核心轉(zhuǎn)錄因子，常以GLK1/GLK2基因?qū)Φ男问酱嬖?，該基因在植物中的功能是保守的，葉綠體發(fā)育受它正調(diào)控，GLKs調(diào)節(jié)一些-Cp和PhANGs（葉綠體和光合作用相關(guān)的核基因）的表達［20］。在草烏葉綠體發(fā)育過程中，（GLK）在幼苗形成階段上調(diào)表達，說明其調(diào)控葉綠體生物發(fā)生和葉綠體分裂過程（圖4）。

在草烏葉綠體發(fā)育過程中的光調(diào)控模型中共鑒定出了18個顯著差異表達的基因；在T4時期，2個DEGs（Phy）上調(diào)表達，1個DEGs（Phy）下調(diào)表達，2個DEGs（Cry）下調(diào)表達，2個DEGs（SPA1）下調(diào)表達，2個DEGs（COP1）上調(diào)表達，1個DEGs（COP1）下調(diào)表達，4個DEGs（PIFs）下調(diào)表達，2個DEGs（HY5）上調(diào)表達；在T3時期，1個DEGs（Phy）上調(diào)表達；在T2時期，1個DEGs（SPA1）上調(diào)表達（圖5）。

2.4 激素調(diào)控草烏葉綠體發(fā)育過程中相關(guān)基因表達模式分析

BR對光形態(tài)發(fā)生和葉綠體發(fā)育的調(diào)控主要由BR信號轉(zhuǎn)導(dǎo)通路的核心轉(zhuǎn)錄因子（BR-activated transcription factor 1，BZR1）及其同源物BES1介導(dǎo)的。BZR1是BR信號轉(zhuǎn)導(dǎo)過程中必不可少的轉(zhuǎn)錄因子，許多靶基因的表達都受它直接調(diào)控［27］，在BR調(diào)控葉綠體發(fā)育途徑中共鑒定出了5個DGEs，2個DEGs（TCH4）在T4時期上調(diào)表達，1個DEGs（TCH4）在T4時期下調(diào)表達，2個DEGs（CYCD3）在T3時期上調(diào)表達（圖6）。

在細胞分裂素調(diào)控葉綠體發(fā)育過程中，1個細胞分裂素受體的His激酶（Arabidopsis His Kinase2）（AHK2）、AHK3和CRE1/AHK4在T4時期上調(diào)表達。與細胞分裂素結(jié)合后，這些受體通過組氨酸磷酸轉(zhuǎn)移蛋白（AHPs）將信號傳遞至B型應(yīng)答調(diào)節(jié)劑B-ARR，ARR（Arabidopsis Response Regulator）激活下游靶基因的轉(zhuǎn)錄［28］，在T4時期，1個DEGs（B-ARR）下調(diào)表達，1個DEGs（GLK1/2）上調(diào)表達，以調(diào)節(jié)葉綠體的發(fā)育（圖6）。

暗形態(tài)建成到光形態(tài)建成期間，DELLA蛋白介導(dǎo)赤霉素信號與光信號之間的互作。本研究中，在T3時期，4個DEGs（DELLA）上調(diào)表達；在T4時期，4個DEGs（PIF4）上調(diào)表達，促進了葉綠體的發(fā)育和下胚軸伸長。本研究中，2個POR酶的DEGs在T4時期上調(diào)表達（圖6）。

在細胞水平上，ARF （Auxin Response Factor） 和AUX/IAA （Auxin/indole-3-Acetic Acid Inducible）蛋白因子調(diào)控生長素，如圖6（D）中的a過程所示，本研究中鑒定到1個DEGs（Aux/IAA）在T3時期上調(diào)表達，7個DEGs（ARF）在T4時期上調(diào)表達，2個DEGs（ARF）在T4時期下調(diào)表達，1個DEGs （GLK1/2）在T4時期上調(diào)表達。如圖6（D）中的b過程所示，本研究中，在T4時期鑒定到9個DEGs（Aux）上調(diào)表達，6個DEGs（ARF）下調(diào)表達，1個DEGs（GLK1/2）上調(diào)表達，因此，GNC和GNL的表達受生長素和赤霉素在轉(zhuǎn)錄水平上的交叉調(diào)控，進而影響葉綠體發(fā)育和葉綠素合成［29］ （圖6）。

2.5 關(guān)鍵基因相對表達量分析

草烏幼苗形成期葉綠體發(fā)育過程中關(guān)鍵基因表達情況如圖7所示，結(jié)果表明，除了PHY、ARF外，DELLA、POR、HY5、AUX/IAA、COPA、CYCD3和CRY的轉(zhuǎn)錄組分析與qRT-PCR分析結(jié)果基本一致。PHY在轉(zhuǎn)錄組分析T4時期表達量下調(diào)，qRT-PCR分析結(jié)果中表達量也下調(diào)，但表達量與轉(zhuǎn)錄組分析數(shù)量有差異。ARF在轉(zhuǎn)錄組分析T3時期表達量上調(diào)，qRT-PCR分析結(jié)果中表達量也上調(diào)，但表達量與轉(zhuǎn)錄組分析數(shù)量間差異較大。

3 討論

葉綠體的發(fā)育在植物萌發(fā)到成苗這一階段至關(guān)重要，影響著植物的整個生命周期。黑暗中，光敏色素相互作用因子（PIF）抑制與光合作用相關(guān)的核編碼（PhANG）轉(zhuǎn)錄，包括編碼葉綠素生物合成酶的基因和（GLK1）轉(zhuǎn)錄因子等。同時，E3泛素連接酶復(fù)合物也會降解大量的轉(zhuǎn)錄因子和活化的光敏色素，特別是（HY5）轉(zhuǎn)錄因子［30］。遺傳學研究發(fā)現(xiàn)，植物的子葉和真葉這兩種葉片里葉綠體的發(fā)育是不同的，二者有自己的發(fā)育通路，但任何一個部位受到影響，都會影響植物的生長［31］。了解葉綠體發(fā)育的調(diào)控機制意義重大。前人研究發(fā)現(xiàn)，在黑暗條件下，GNC和GNL（GATA家族轉(zhuǎn)錄因子）的過表達促進了原質(zhì)體向黃化質(zhì)體的分化；在光照下，促進葉綠體發(fā)育［32］。本研究中光信號調(diào)控葉綠體發(fā)育過程中發(fā)現(xiàn)1個DEGs（GLK1/2）在T4時期顯著下調(diào)，2個DEGs（CRY）在T4時期顯著下調(diào)，4個DEGs（PIFs）和2個DEGs（HY5）在T4時期顯著上調(diào)，激活下游轉(zhuǎn)錄因子促進葉綠體發(fā)育。

葉綠體的生物發(fā)生由核和葉綠體編碼的基因產(chǎn)物共同決定。除了光以外，一些激素如生長素、細胞分裂素、赤霉素等都參與調(diào)節(jié)葉綠體發(fā)育。植物激素可以通過控制細胞分裂和擴張之間的平衡，或通過影響光合作用相關(guān)基因的表達、葉綠體的生物合成和分裂來影響細胞中葉綠體的大小，從而影響其發(fā)育。在光形態(tài)建成過程中BR信號被植物接收后，BRZ1受體會結(jié)合GATA2的啟動子，因此GATA2的表達被抑制；光形態(tài)建成相關(guān)因子COP1的水解被光抑制，GATA2被激活跟光響應(yīng)基因的啟動子元件結(jié)合，光和BR響應(yīng)基因的表達被調(diào)控，因此植物幼苗生長被影響，在光形態(tài)建成過程中，GATA2起正調(diào)節(jié)因子的作用［33］。HY5是連接光信號與激素信號的重要節(jié)點，通過整合光信號以此來影響激素含量來促進光形態(tài)建成。例如，HY5可以與BZR1相互作用拮抗BZR1來調(diào)控子葉開放相關(guān)基因的轉(zhuǎn)錄活性［34］，本研究在BR調(diào)控葉綠體發(fā)育途徑中發(fā)現(xiàn)2個DEGs（TCH4）在T4時期上調(diào)表達，1個DEGs（TCH4）在T4時期下調(diào)表達，2個DEGs（CYCD3）在T3時期上調(diào)表達，在草烏葉綠體發(fā)育過程中起調(diào)控作用。

番茄（Solanum lycopersicum L.）中AtBZR1-1D與SlGLK2協(xié)同上調(diào)促進葉綠素合成與葉綠體發(fā)育［35］，與本研究的結(jié)果一致。CRF2作為細胞分裂素反應(yīng)因子家族的成員，屬于AP2/ERF類轉(zhuǎn)錄因子。CRF2過表達導(dǎo)致了PDV2蛋白水平的升高和葉綠體分裂的增加。作為參與葉綠體發(fā)育和分裂的重要調(diào)控因子CGA1和CRF2被發(fā)現(xiàn)在調(diào)控細胞分裂素信號上起著極其重要的作用［15］。本研究在CK調(diào)控葉綠體發(fā)育途徑中發(fā)現(xiàn)3個DGEs影響其發(fā)育，1個DEGs（AHK2/3/4）和1個DEGs（GLK1/2）在T4時期顯著上調(diào)，1個DEGs（B-ARR）在T4時期顯著下調(diào)。GA能抑制CGA1和GNC這兩個葉綠體發(fā)育促進因子的表達，進而抑制CGA1和GUC下游相關(guān)葉綠體發(fā)育基因的表達［20］。光可以誘導(dǎo)POR基因表達，POR酶可以促進葉綠素合成［30］，在GA調(diào)控葉綠體發(fā)育途徑中發(fā)現(xiàn)10個DGEs調(diào)控其發(fā)育，4個DEGs（PIF4）、2個DEGs（POR）在T4時期顯著上調(diào)，其余基因既有上調(diào)也有下調(diào)。原葉綠素酸酯氧化還原酶（POR）是編碼光合體中的一個關(guān)鍵酶，參與葉綠素a的合成過程［36］，本研究也證明POR對草烏葉綠體發(fā)育起到了一定的調(diào)控作用。Auxin調(diào)控葉綠體發(fā)育途徑中發(fā)現(xiàn)26個DGEs調(diào)控其發(fā)育，10個DEGs（AUX/IAA）和2個DEGs（GLK1/2）在T4時期顯著上調(diào)，9個DEGs（ARF）在T4時期顯著下調(diào)，5個DEGs（ARF）在T4時期顯著上調(diào)。綜上所述，光與其他刺激如激素有交叉作用，共同調(diào)節(jié)植物發(fā)育，我們應(yīng)把焦點聚集在將復(fù)雜的信號轉(zhuǎn)導(dǎo)網(wǎng)絡(luò)系統(tǒng)與功能基因組的信息整合起來，以此來正確認識植株水平上的光調(diào)控發(fā)育機理，有助于明晰環(huán)境因素如何調(diào)控植物的生長發(fā)育［37］。

4 結(jié)論

對草烏4個生長發(fā)育階段的轉(zhuǎn)錄組數(shù)據(jù)分析，共篩選出79 251個基因，T3vsT4時期差異表達基因在GO富集分析中主要富集到參與對幾丁質(zhì)的反應(yīng)、光合系統(tǒng)Ⅰ、脂肪酶活性等過程中。在KEGG富集分析中主要富集到植物-病原體相互作用、植物激素信號轉(zhuǎn)導(dǎo)、光合生物中的碳固定。同時，草烏的4個不同時期在光和激素調(diào)節(jié)葉綠體發(fā)育的不同階段其表達量不同，19個差異表達基因注釋光信號調(diào)控葉綠體發(fā)育通路上。44個差異表達基因注釋植物激素信號轉(zhuǎn)導(dǎo)通路上。本研究挖掘到一部分草烏幼苗形成期葉綠體發(fā)育相關(guān)基因，為進一步研究草烏生長發(fā)育相關(guān)分子機制提供理論基礎(chǔ)。

參考文獻

［1］ 張躍華. 草烏頭種子休眠及繁殖特性研究［D］. 呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學，2016：1

［2］ 彭劭.草烏中生物堿類化學成分的研究［D］. 長春：吉林大學，2014：1

［3］ 山丹. 不同種源草烏生物學特性及入藥部位生物堿含量的研究［D］. 呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學，2020：1

［4］ 王小菁. 植物生理學［M］. 第8版.北京：高等教育出版社，2019：236285

［5］ 鄭素玲. 談植物的光形態(tài)建成［J］. 唐山師專學報，1999（5）：75-76

［6］ MOON J，ZHU L，SHEN H，et al. PIF1 directly and indirectly regulates chlorophyll biosynthesis to optimize the greening process in Arabidopsis［J］. Proceedings of the National Academy of Sciences，2008，105（27）：9433-9438

［7］ LU X D，ZHOU C M，XU P B，et al. Red-light-dependent interaction of phyB with SPA1 promotes COP1-SPA1 dissociation and photomorphogenic development in Arabidopsis［J］. Molecular Plant，2015，8（3）：467-78

［8］ SHEERIN D J，MENON C，ZUR OVEN-KROCKHAUS S，et al. Light-activated phytochrome A and B interact with members of the SPA family to promote photomorphogenesis in Arabidopsis by reorganizing the COP1/SPA complex［J］. Plant Cell，2015，27（1）：189-201

［9］ BAUER D，VICZIN A，KIRCHER S，et al. Constitutive photomorphogenesis 1 and multiple photoreceptors control degradation of phytochrome interacting factor 3，a transcription factor required for light signaling in Arabidopsis［J］. Plant Cell，2004，16：1433-1445

［10］ZHONG S，ZHAO M，SHI T，et al. EIN3/EIL1 cooperate with PIF1 to prevent photo-oxidation and to promote greening of Arabidopsis seedlings［J］. Proceedings of the National Academy of Sciences，USA，2009，106（50）：21431-21436

［11］LEIVAR P，MONTE E，OKA Y，et al. Multiple phytochrome-interacting bHLH transcription factors repress premature seedling photomorphogenesis in darkness［J］. Current Biology， 2008，18（23）：1815-1823

［12］SHI H，SHEN X，LIU R，et al. The Red Light Receptor Phytochrome B Directly Enhances Substrate-E3 Ligase Interactions to Attenuate Ethylene Responses［J］. Developmental Cell， 2016，39（5）：597-610

［13］LIU X，LIU R，LI Y，et al. EIN3 and PIF3 Form an Interdependent Module That Represses Chloroplast Development in Buried Seedlings［J］. Plant Cell， 2017，29（12）：3051-3067

［14］JIAO Y，LAU OS，DENG XW. Light-regulated transcriptional networks in higher plants［J］. Nature Reviews Geneticst，2007（8）：217-230

［15］HUQ E，AL-SADY B，HUDSON M，et al. Phytochrome-interacting factor 1 is a critical bHLH regulator of chlorophyll biosynthesis［J］. Science，2004，305：1937-1941

［16］GOSLINGS D，MESKAUSKIENE R，KIM C，et al. Concurrent interactions of heme and FLU with Glu tRNA reductase （HEMA1），the target of metabolic feedback inhibition of tetrapyrrole biosynthesis，in dark and light-grown Arabidopsis plants［J］. Plant Journal，2004，40：957-96

［17］SHIRR J，KIM K，KANG H，et al. Phytochromes promote seedling light responses by inhibiting four negatively-acting phytochrome-interacting factors［J］. Proceedings of the National Academy of Sciences，USA，2009，106（18）：7660-7665

［18］李保珠，趙孝亮，彭雷. 植物葉綠體發(fā)育及調(diào)控研究進展［J］. 植物學報，2014，49（3）：337-345

［19］唐鈺瑩，劉陽軒，潘婷，等. 葉綠體蛋白CV調(diào)節(jié)質(zhì)體逆向信號的研究［J］. 植物生理學報，2020，56（7）：1553-1560

［20］陳煒. 調(diào)控水稻葉綠體發(fā)育基因YGL3的圖位克隆及功能研究［D］. 南昌：江西農(nóng)業(yè)大學，2022：5-6

［21］YOUMING Y，CHUNIAN X，BAOMIN W，et al. Effects of plant growth regulators on secondary wall thickening of cotton fibres［J］. Plant Growth Regulation，2004，35（3）：233-237

［22］ROBERTS E H. Seed dormancy and oxidation processes［J］. Symposia of the Society for Experimental Biology，1969，23：161-192

［23］GRABHERR M G，HAAS B J，YASSOUR M，et al. Trinity：reconstructing a full-length transcriptome without a genome from RNA-Seq data［J］. Nature Biotechnology，2011，29（7）：644

［24］洪森榮，劉佳凝，袁昕，等. 蘇丹草和高丹草轉(zhuǎn)錄組測序及其差異基因表達分析［J］. 草地學報，2024，32（3）：714-725

［25］BENIAMIN B， XIE C， HUSON D H. Fast and sensitive protein alignment using DIAMOND［J］. Nature Methods，2015，12：59-60

［26］金忠民，李春月，劉本松，等. 菌株JB12影響鉛鎘脅迫下菊苣黃酮合成的轉(zhuǎn)錄組分析［J］. 草地學報，2023，31（6）：1648-1655

［27］王晉，王新宇，沈淵博，等. 番茄果實葉綠體發(fā)育調(diào)控及其應(yīng)用的研究進展［J］. 園藝學報，2022，49（12）：2669-2682

［28］衛(wèi)云豐. 擬南芥PGA37轉(zhuǎn)錄因子調(diào)控葉綠體發(fā)育的分子機制［D］. 晉中：山西農(nóng)業(yè)大學，2020：4

［29］SAKAMOTO W，MIYAGISHIMA S Y，JARVIS P. Chloroplast biogenesis：control of plastid development，protein import，division and inheritance ［J］. Arabidopsis Book，2008，6：e0110

［30］吳釩漳. miR319-TCP4調(diào)控葉綠體發(fā)育的分子機制研究［D］. 昆明：昆明理工大學，2022：4

［31］SAKAI A，TAKANO H，KUROIWA T. Organelle nuclei in higher plants：structure，composition，function，and evolution ［J］. International Review of Cytology，2004，238：59-118

［32］王娟，蘭海燕. GATA轉(zhuǎn)錄因子對植物發(fā)育和脅迫響應(yīng)調(diào)控的研究進展［J］. 植物生理學報，2016，52（12）：1785-1794

［33］CHIANG Y H，ZUBO Y O，TAPKEN W，et al. Functional characterization of the GATA transcription factors GNC and CGA1 reveals their key role in chloroplast development，growth，and division in Arabidopsis［J］. Plant Physiol，2012，160（1）：332-348

［34］詹為民. 利用轉(zhuǎn)錄組和代謝組解析玉米光信號與雜種優(yōu)勢的互作機制［D］. 鄭州：河南農(nóng)業(yè)大學，2023：5

［35］LIU L，JIA C，ZHANG M，et al. Ectopic expression of a BZR1-1D transcription factor in brassinosteroid signalling enhances carotenoid accumulation and fruit quality attributes in tomato［J］. Plant Biotechnology Journal，2014，12（1）：105-115

［36］LIANG M，GU D，LIE Z，et al. Regulation of chlorophyll biosynthesis by light-dependent acetylation of NADPH：protochlorophyll oxidoreductase A in Arabidopsis［J］. Plant Science，2023，330：111641

［37］王小菁. 我國光形態(tài)建成研究回顧［J］. 植物學通報，2003（4）：407-415

（責任編輯 彭露茜）