馬李廣 張賀龍 龐小可 王浩 李廣 汪承剛 袁凌云 侯金鋒 唐小燕 陳國(guó)戶

摘要： 為研究白菜（ Brassica rapa ）堿性亮氨酸拉鏈（Basic leucine zipper，bZIP）轉(zhuǎn)錄因子家族基因（ BrbZIP ）的相關(guān)功能，通過(guò)生物信息學(xué)分析技術(shù)，鑒定全基因組的 BrbZIP 基因家族成員，并對(duì)其染色體分布、進(jìn)化關(guān)系、表達(dá)模式及其應(yīng)答春化反應(yīng)等進(jìn)行了分析。白菜bZIP轉(zhuǎn)錄因子基因家族共有118個(gè)成員，在染色體上不均等分布。白菜組織轉(zhuǎn)錄組分析結(jié)果顯示，大部分 BrbZIP 基因在根、莖、葉、花及莢中均有較高的表達(dá)豐度，且具有組織表達(dá)特異性;春化反應(yīng)轉(zhuǎn)錄組、熒光定量PCR及相關(guān)基因互作網(wǎng)絡(luò)分析結(jié)果表明，白菜 bZIP 基因家族中應(yīng)答春化反應(yīng)相關(guān)基因上調(diào)與下調(diào)表達(dá)基因數(shù)量差異不大，相互之間存在復(fù)雜的互作網(wǎng)絡(luò)，其中4個(gè)成員（Bra039631、Bra020620、Bra004550與Bra020471）是應(yīng)答春化反應(yīng)的中心節(jié)點(diǎn)。

關(guān)鍵詞： 白菜; bZIP轉(zhuǎn)錄因子; 基因家族; 春化反應(yīng); 生物信息學(xué)分析

中圖分類號(hào)： S634.3?? 文獻(xiàn)標(biāo)識(shí)碼： A?? 文章編號(hào)： 1000-4440（2022）03-0765-10

Genome-wide identification of bZIP transcription factor gene family in ?Brassica rapa ?and its association with vernalization

MA Li-guang 1 ， ZHANG He-long 1 ， PANG Xiao-ke 1 ， WANG Hao 1 ， LI Guang 1，2 ， WANG Cheng-gang 1，2，3 ，YUAN Ling-yun 1，2 ， HOU Jin-feng 1，2 ， TANG Xiao-yan 1，2 ， CHEN Guo-hu 1，2，3

（1.Anhui Provincial Engineering Laboratory of Horticultural Crop Breeding， School of Horticulture， Anhui Agricultural University， Hefei 230036， China; 2.Wanjiang Vegetable Industrial Technology Institute， Maanshan 238200， China; 3.Anhui Southwest Comprehensive Experimental Station of Anhui Agricultural University， Huaining 246100， China）

Abstract： In order to explore the function of basic leucine zipper （bZIP） transcription factor in ?Brassica rapa ?and its association with vernalization， bioinformation methods were used to identify the ?B. rapa ?bZIP gene family members. The chromosome location， evolutionary relationship， expression pattern and vernalization response were analyzed through a variety of biological information softwares. The results revealed that 118 ?BrbZIP ?transcription factors were unevenly distributed on the chromosomes of ?B. rapa . Transcriptome analysis showed that most of the ?BrbZIP ?genes had high expression levels in roots， stems， leaves， flowers and pods， and had tissue expression specificity. The results of transcriptome， fluorescence quantitative PCR and interaction network analysis suggested that there was n o significant difference in the number of up-regulated an d down-regulated genes responding to vernalization in bZIP gene family of ?B. rapa ， and the interaction network was complex. Four ?BrbZIP ?genes were the central nodes in response to vernalization.

Key words： ?Brassica rapa ; bZIP transcription factor; gene family; vernalization; bioinformatics analysis

轉(zhuǎn)錄因子（Transcription factors，TFs），是一類通過(guò)與下游基因啟動(dòng)子順式元件特異性結(jié)合或與其他蛋白質(zhì)互作來(lái)調(diào)控基因表達(dá)的一類蛋白質(zhì)。在只存在于真核生物中的轉(zhuǎn)錄因子家族中，堿性亮氨酸拉鏈（Basic leucine zipper，bZIP）家族是轉(zhuǎn)錄因子中成員數(shù)量最多、最保守的家族之一，在植物生長(zhǎng)發(fā)育和應(yīng)答脅迫反應(yīng)等生物過(guò)程中起著重要的調(diào)控作用 [1] 。

bZIP蛋白最顯著的特征是含有一個(gè)保守的60～ 80個(gè)氨基酸組成的bZIP結(jié)構(gòu)域（ N -×7-R/K-×9- L -×6- L ），包括1個(gè)堿性氨基酸區(qū)域和1個(gè)相鄰的亮氨酸拉鏈區(qū);堿性氨基酸區(qū)域位于C端，由16～20個(gè)氨基酸殘基組成 N -×7-R/K基序，負(fù)責(zé)與DNA啟動(dòng)子區(qū)域內(nèi)特定序列結(jié)合;亮氨酸（Leu）拉鏈區(qū)位于 N 端，由亮氨酸或異亮氨酸（Ile）、纈氨酸（Val）組成的七肽重復(fù)序列（×9- L -×6- L -×6- L ），主要作用為bZIP蛋白特定識(shí)別和形成同源二聚體 [2] 。植物bZIP蛋白優(yōu)先與含有ACGT核心元件的DNA序列結(jié)合，特別是G-box（CACGTG）、C-box（GACGTC）以及A-box（TACGTA）等元件 [3] 。

目前bZIP轉(zhuǎn)錄因子基因家族已在多個(gè)植物基因組中被鑒定出來(lái)，如擬南芥（ Arabidopsis thaliana ） [4] 、油菜（ B.napus ） [5] 、麻瘋樹（ Jatropha curcas ） [6] 、藜麥（ Chenopodium quinoa ） [7] 、楊樹（ Populus tremula ） [8] 、獼猴桃（ Actinidia chinensis ） [9] 等。大量研究結(jié)果表明，bZIP轉(zhuǎn)錄因子在植物光信號(hào)、花的誘導(dǎo)與發(fā)育、種子成熟與萌發(fā)、非生物和生物脅迫等多個(gè)關(guān)鍵生物學(xué)過(guò)程中起著重要的作用。例如，在擬南芥中，bZIP11基因可影響根系發(fā)育 [10] ，HY5基因參與調(diào)控根和下胚軸的發(fā)育 [11] ;而bZIP1基因與花藥發(fā)育過(guò)程中的水分運(yùn)動(dòng)密切相關(guān)，其過(guò)量表達(dá)不利于花粉發(fā)育 [12] 。小麥（ Triticum aestivum ）bZIP轉(zhuǎn)錄因子 TaGBF 基因調(diào)控 CO 及 FT 等基因的表達(dá)，影響植物開花 [13] 。大豆（ Glycine max ）植株中FDL19蛋白與FT2a、FT5a相互作用，可影響大豆的開花時(shí)間 [14-15] 。風(fēng)疹卷尾草（ Capsella rubella ）TGA4基因的過(guò)表達(dá)，可延遲擬南芥開花時(shí)間 [16] 。編碼bZIP轉(zhuǎn)錄因子的HY5基因，在植物幼苗光形態(tài)建成期間，對(duì)激活光反應(yīng)基因起著重要的作用 [17] 。水稻（ Oryza sativa ）應(yīng)激激活蛋白激酶SAPK100磷酸化OsbZIP77，導(dǎo)致開花時(shí)間提前 [18] ;OsFD4基因通過(guò)作用于莖尖分生組織中成花整合因子復(fù)合物的形成，促進(jìn)水稻開花 [19] 。因此，bZIP轉(zhuǎn)錄因子在植物開花過(guò)程中具有重要的作用。

目前，bZIP轉(zhuǎn)錄因子在植物開花過(guò)程中的功能研究，已有一些進(jìn)展。然而迄今為止，未有白菜bZIP轉(zhuǎn)錄因子參與春化反應(yīng)的研究報(bào)道。本研究擬對(duì)白菜（ B. rapa ） bZIP 基因家族進(jìn)行系統(tǒng)研究，包括基因家族成員鑒定、蛋白質(zhì)氨基酸序列分析和系統(tǒng)發(fā)育關(guān)系、基因的染色體分布、基因組串聯(lián)復(fù)制和片段復(fù)制、基因共線性分析等，此外，基于轉(zhuǎn)錄組數(shù)據(jù)及熒光定量PCR分析，探討白菜 bZIP 基因在不同組織中的表達(dá)模式及應(yīng)答春化反應(yīng)相關(guān)基因的表達(dá)，為進(jìn)一步解析白菜 bZIP 基因在春化過(guò)程中的功能奠定基礎(chǔ)。

1 材料與方法

1.1 數(shù)據(jù)來(lái)源

從白菜基因組數(shù)據(jù)庫(kù)（BRAD;http：//brassicadb.cn/）下載白菜全基因組數(shù)據(jù)（V1.5）及其注釋文件，從NCBI Sequence Read Archive（SRA）數(shù)據(jù)庫(kù)下載白菜不同組織（根、莖、葉、花、莢等） [20] 及春化反應(yīng)（JWW材料春化0 d、25 d、35 d;XBJ材料春化0 d、10 d、25 d） [21] RNA-Seq原始數(shù)據(jù)。RNA-Seq原始數(shù)據(jù)經(jīng)過(guò)質(zhì)量控制過(guò)濾后，比對(duì)至白菜基因組（V1.5）序列，再利用FeatureCounts進(jìn)行有參轉(zhuǎn)錄組定量?；虮磉_(dá)豐度采用FPKM值進(jìn)行計(jì)算，差異表達(dá)基因采用DESeq軟件包進(jìn)行計(jì)算 [22-23] 。

1.2 白菜 bZIP 基因家族成員的鑒定

從白菜基因組數(shù)據(jù)庫(kù)（BRAD）下載 bZIP 基因家族候選成員ID，利用TBtools軟件 [24] 提取白菜 bZIP 基因家族候選基因序列、蛋白質(zhì)氨基酸序列等信息。使用Pfam（http：//pfam.xfam.org/）、SwissProt（https：//www.expasy.org/resources/uniprotkb-swiss-prot）、SMART（http：//smart.embl-heidelberg.de/）和InterProScan（https：//www.ebi.ac.uk/interpro/about/interproscan/）等數(shù)據(jù)庫(kù)鑒定所有候選蛋白質(zhì)氨基酸序列。對(duì)確定BrbZIP成員的蛋白質(zhì)氨基酸序列，利用MEGA 7.0軟件采用鄰接法以及泊松校正與成對(duì)刪除法，構(gòu)建系統(tǒng)進(jìn)化樹 [25] 。

1.3 白菜bZIP家族理化性質(zhì)及 bZIP 基因組信息分析

采用ExPasy工具（https：//web.expasy.org/compute_pi/），分析白菜bZIP家族的相對(duì)分子質(zhì)量、等電點(diǎn)和長(zhǎng)度等理化性質(zhì)。提取白菜 bZIP 基因家族注釋信息，并定位到染色體上，利用TBtools工具進(jìn)行可視化，獲得 BrbZIP 基因家族成員染色體定位圖譜。將白菜 bZIP 基因家族編碼的蛋白質(zhì)氨基酸序列進(jìn)行all-BLASTP-all分析（ top <5， E <1e -5 ）。采用TBtools軟件中MCScanX工具提取共線性基因?qū)?，再利用Circos工具進(jìn)行可視化。

1.4 白菜 bZIP 基因家族組織表達(dá)分析及其應(yīng)答春化反應(yīng)分析

利用白菜不同組織轉(zhuǎn)錄組數(shù)據(jù) [20] ，分析 BrbZIP 基因家族的組織表達(dá)特征。利用白菜春化反應(yīng)材料轉(zhuǎn)錄組數(shù)據(jù) [21] ，分析白菜 bZIP 基因家族應(yīng)答春化反應(yīng)。利用STRING數(shù)據(jù)庫(kù)（https：//string-db.org/）分析白菜應(yīng)答春化相關(guān) bZIP 基因互作網(wǎng)絡(luò)關(guān)系。對(duì) BrbZIP 基因家族中應(yīng)答春化反應(yīng)關(guān)鍵基因，利用熒光定量PCR技術(shù)分析其在JWW不同春化時(shí)間材料中的相對(duì)表達(dá)量 [21] 。

2 結(jié)果與分析

2.1 白菜 bZIP 基因家族鑒定及其編碼的蛋白質(zhì)理化特性分析

經(jīng)篩選鑒定，在白菜基因組中共獲得118個(gè) bZIP 基因家族成員。理化性質(zhì)分析發(fā)現(xiàn)，白菜 bZIP 基因家族開放閱讀框（ORF）長(zhǎng)度為273～ 1 893 ?bp，編碼90～ 630個(gè)氨基酸（AA），預(yù)測(cè)相對(duì)分子質(zhì)量大小在1.062× 10 4 至6.874× 10 4 之間，等電點(diǎn)（pI）范圍為4.80～ 10.11。對(duì)白菜bZIP蛋白的等電點(diǎn)與相對(duì)分子質(zhì)量之間進(jìn)行比較分析，發(fā)現(xiàn)酸性與堿性BrbZIP蛋白數(shù)量相差不大，且相對(duì)分子質(zhì)量分布較規(guī)律。對(duì)白菜bZIP蛋白進(jìn)行亞細(xì)胞定位預(yù)測(cè)分析，結(jié)果顯示絕大部分白菜bZIP蛋白家族成員定位于細(xì)胞核內(nèi)，只有6個(gè)成員定位于細(xì)胞膜、線粒體或葉綠體中，表明白菜bZIP蛋白家族的功能主要發(fā)生在細(xì)胞核中。

2.2 白菜 bZIP 基因家族染色體定位分析

利用白菜基因組注釋信息，將白菜 bZIP 基因家族定位于染色體上。結(jié)果如圖1所示，118個(gè)白菜 bZIP 家族基因不均等分布在10條染色體上，各染色體上的基因數(shù)量與染色體大小無(wú)關(guān)。其中，染色體A09上的 BrbZIP 基因家族數(shù)量最多（23個(gè)），染色體A10上的 BrbZIP 基因家族數(shù)量最少（7個(gè)）。此外，白菜 bZIP 家族基因在A04（Bra016959與Bra016953）、A09（Bra007274與Bra007276）染色體上各有1對(duì)基因串聯(lián)重復(fù)，表明白菜 bZIP 家族基因 在進(jìn)化過(guò)程中發(fā)生了少量的串聯(lián)重復(fù)現(xiàn)象。

2.3 白菜與擬南芥 bZIP 基因家族系統(tǒng)發(fā)育分析

為了確定白菜 bZIP 基因家族成員的進(jìn)化關(guān)系，通過(guò)MEGA-X軟件中MUSCLE工具將118個(gè)白菜 bZIP 基因與73個(gè)擬南芥 bZIP 基因進(jìn)行序列分析比對(duì)，運(yùn)用鄰接法構(gòu)建了系統(tǒng)發(fā)育樹（圖2）。白菜與擬南芥 bZIP 基因家族分成了5個(gè)大簇（GroupA～ GroupE）、11個(gè)小簇。其中，Group E的 bZIP 基因數(shù)量最多，包括21個(gè) AtbZIP 基因和41個(gè) BrbZIP 基因;Group C的數(shù)量最少，僅有16個(gè)，包括7個(gè) AtbZIP 基因和9個(gè) BrbZIP 基因。系統(tǒng)發(fā)育樹分析結(jié)果表明，各簇中白菜和擬南芥 bZIP 基因具有較高的同源性。

2.4 白菜 bZIP 家族基因共線性分析

為了分析白菜 bZIP 家族基因的擴(kuò)張模式，對(duì)其基因復(fù)制事件進(jìn)行了分析。如圖3所示，118個(gè)白菜 bZIP 家族基因成員之間，有110對(duì)存在共線性基因?qū)?，其中存在著單個(gè)基因與多個(gè)基因?qū)?yīng)的情況，表明白菜 bZIP 家族在發(fā)育進(jìn)化過(guò)程中發(fā)生了大量的基因復(fù)制現(xiàn)象，是白菜 bZIP 基因多樣性的主要驅(qū)動(dòng)力。

2.5 白菜 bZIP 家族基因表達(dá)模式及響應(yīng)春化反應(yīng)分析

為了研究白菜 bZIP 基因在不同組織中的表達(dá)模式，利用白菜組織轉(zhuǎn)錄組數(shù)據(jù)分析了 BrbZIP 基因在根、莖、葉、花、莢等組織中的表達(dá)情況。結(jié)果如圖4A所示，大部分白菜 bZIP 基因在各組織中均有表達(dá)，表明這些 bZIP 基因在白菜整個(gè)生長(zhǎng)發(fā)育階段均有重要的作用。對(duì) FPKM 值大于500的 BrbZIP 基因進(jìn)行了維恩分析，發(fā)現(xiàn)各組織共有 BrbZIP 基因僅13個(gè)（圖4B）。此外，少數(shù) BrbZIP 基因具有較強(qiáng)的組織表達(dá)特異性，如Bra037809、Bra018634、Bra031871、Bra023540、Bra031622等主要在根系中表達(dá)，Bra010504、Bra004329等主要在莖中表達(dá)，Bra025144主要在花中表達(dá)，而Bra009288、Bra035957等在莢中表現(xiàn)出較高的轉(zhuǎn)錄豐度，表明 BrbZIP 基因的表達(dá)具有組織特異性（圖4）。

我們還利用白菜春化反應(yīng)轉(zhuǎn)錄組數(shù)據(jù)分析了應(yīng)答春化的相關(guān) BrbZIP 基因，結(jié)果如圖5所示，在白菜JWW材料與XBJ材料春化轉(zhuǎn)錄組數(shù)據(jù)中，共篩選出42個(gè)差異表達(dá) BrbZIP 基因（DEGs）（圖5A），這些DEGs中，JWW材料與XBJ材料共有18個(gè)（圖5B），其中8個(gè)DEGs上調(diào)表達(dá)，10個(gè)DEGs下調(diào)表達(dá)（圖5C）。對(duì)白菜JWW材料與XBJ材料中春化相關(guān) BrbZIP 基因分別進(jìn)行表達(dá)趨勢(shì)分析發(fā)現(xiàn)，共有5個(gè)顯著富集趨勢(shì)（ P < 0.05）。在JWW材料中，18個(gè) BrbZIP 基因隨春化時(shí)間延長(zhǎng)，轉(zhuǎn)錄豐度逐漸增大，其中12個(gè) BrbZIP 基因在春化25 d時(shí)表達(dá)量最大（Profile 6）、6個(gè) BrbZIP 基因在春化35 d時(shí)表達(dá)量最大（Profile 7）;10個(gè) BrbZIP 基因隨春化時(shí)間延長(zhǎng)轉(zhuǎn)錄豐度逐漸降低（Profile 1）。在XBJ材料中，16個(gè) BrbZIP 基因隨春化時(shí)間延長(zhǎng)，轉(zhuǎn)錄豐度逐漸增大，在春化10 d時(shí)表達(dá)量最大（Profile 6）;7個(gè) BrbZIP 基因隨春化時(shí)間延長(zhǎng)轉(zhuǎn)錄豐度逐漸降低（Profile 0）。在JWW材料與XBJ材料共有模塊6（Profile 6）中，共有9個(gè) BrbZIP 基因（圖6），表明這些 BrbZIP 基因在白菜春化反應(yīng)中具有重要的作用。

2.6 白菜 bZIP 家族應(yīng)答春化相關(guān)基因互作網(wǎng)絡(luò)分析

為了進(jìn)一步研究白菜春化應(yīng)答 BrbZIP 基因的相互作用，利用STRING數(shù)據(jù)庫(kù)對(duì)JWW材料與XBJ材料中共有18個(gè)應(yīng)答春化 BrbZIP 基因進(jìn)行互作網(wǎng)絡(luò)分析（圖7）。結(jié)果顯示18個(gè) BrbZIP 基因中有11個(gè)與其他基因有相互作用關(guān)系，且互作網(wǎng)絡(luò)表現(xiàn)出復(fù)雜的相關(guān)性，其中Bra039631（AtbZIP2）、Bra020620（AtbZIP63）、Bra004550（AtbZIP55）、Bra020471（AtbZIP9）基因在互作網(wǎng)絡(luò)中是應(yīng)答春化反應(yīng)的中心節(jié)點(diǎn)，表明其在應(yīng)答春化反應(yīng)調(diào)控網(wǎng)絡(luò)中具有重要的作用。

2.7 白菜 bZIP 家族應(yīng)答春化相關(guān)基因熒光定量PCR分析

利用熒光定量PCR技術(shù)，進(jìn)一步分析白菜 bZIP 家族應(yīng)答春化反應(yīng)關(guān)鍵基因的表達(dá)水平。結(jié)果（圖8）顯示，在春化35 d內(nèi)，Bra039631、Bra020620、Bra033464、Bra001742、Bra017316基因的表達(dá)量均低于春化前，而其他6個(gè)基因的表達(dá)水平均高于春化前。以上結(jié)果表明，上述 BrbZIP 基因?yàn)榘撞舜夯瘧?yīng)答基因，在白菜春化過(guò)程中具有重要的作用。

3 討 論

bZIP轉(zhuǎn)錄因子基因家族是植物中數(shù)量比較多的基因家族，不同植物家族成員數(shù)量不等，但在進(jìn)化上比較保守 [3] 。本研究經(jīng)過(guò)嚴(yán)格篩選鑒定，從白菜基因組中獲得118個(gè) bZIP 基因，數(shù)量是擬南芥 bZIP 成員的1.6倍，與白菜 PRX 家族基因數(shù)量比值類似 [23] ，表明 BrbZIP 家族基因成員在白菜基因組倍化過(guò)程中，同樣發(fā)生了基因丟失事件。通過(guò)系統(tǒng)發(fā)育樹分析將已鑒定的白菜與擬南芥 bZIP 基因成員分為5大簇11個(gè)小簇，與擬南芥根據(jù)結(jié)構(gòu)分組的結(jié)果類似 [4] ，表明白菜與擬南芥的 bZIP 基因在進(jìn)化上是保守的。

在物種進(jìn)化過(guò)程中，基因家族通常經(jīng)歷串聯(lián)復(fù)制或大規(guī)模的片段復(fù)制，以維持家族較大規(guī)模 [8， 26] 。白菜經(jīng)歷過(guò)2次全基因組復(fù)制，及1次全基因組三倍化事件 [27] 。然而，白菜118個(gè) bZIP 基因只有2對(duì)串聯(lián)重復(fù)，數(shù)量顯著少于油菜（247個(gè)成員，15個(gè)基因復(fù)制對(duì)） [5] ，但白菜基因組內(nèi)共線性基因?qū)Ω哌_(dá)110組，表明基因復(fù)制事件在蕓薹屬作物進(jìn)化中的復(fù)雜性，以及在 bZIP 基因家族擴(kuò)展中的重要性 [5， 28] 。此外，通過(guò)組織表達(dá)模式分析也發(fā)現(xiàn)，一些重復(fù)的 BrbZIP 的轉(zhuǎn)錄水平也相似，這可能與它們高度相似的結(jié)構(gòu)和順式調(diào)控元件有關(guān)。

對(duì)于多基因家族來(lái)說(shuō)，基因表達(dá)分析常常為功能預(yù)測(cè)提供有用的線索。本研究利用白菜組織轉(zhuǎn)錄組數(shù)據(jù) [20] 分析了 BrbZIP 在各組織中的表達(dá)模式。結(jié)果表明，大部分 BrbZIP 在白菜根、莖、葉、花、莢中均有較高的表達(dá)豐度，且在根組織中高表達(dá)豐度基因數(shù)量最多，這可能與根部應(yīng)答逆境相關(guān) [1] ;在各組織高表達(dá)豐度（ FPKM >500）的 BrbZIP 基因中，各組織僅共有13個(gè) BrbZIP 基因，說(shuō)明大部分 BrbZIP 基因具有組織表達(dá)特異性，且這些組成型表達(dá) BrbZIP 基因在白菜生長(zhǎng)發(fā)育過(guò)程中具有重要的作用 [8] 。

由于目前對(duì) bZIP 基因應(yīng)答春化反應(yīng)的調(diào)控作用知之甚少，我們探討了 bZIP 基因在白菜幼苗春化反應(yīng)中的表達(dá)模式。比較春化前后樣本的表達(dá)數(shù)據(jù) [21] ，結(jié)果顯示白菜 bZIP 家族基因中應(yīng)答春化反應(yīng)表達(dá)有上調(diào)也有下調(diào)，結(jié)果與大麥（ Hordeum vulgare ）中類似 [29] 。在白菜兩個(gè)材料中應(yīng)答春化反應(yīng) bZIP 基因有42個(gè)，共有差異表達(dá)基因18個(gè)，其中上調(diào)基因8個(gè)，下調(diào)基因10個(gè)，表明這些基因與春化反應(yīng)密切相關(guān)。為了了解它們的功能，對(duì)白菜18個(gè)春化應(yīng)答關(guān)鍵 bZIP 基因進(jìn)行了基因互作網(wǎng)絡(luò)分析及熒光定量PCR分析，發(fā)現(xiàn)這些基因表現(xiàn)出較復(fù)雜的互作關(guān)系，其中下調(diào)表達(dá)基因Bra039631與Bra020620，及上調(diào)表達(dá)基因Bra004550與Bra020471作為互作網(wǎng)絡(luò)的中心節(jié)點(diǎn)，在白菜應(yīng)答春化過(guò)程中發(fā)揮著重要作用。

參考文獻(xiàn)：

[1] 魏瑞敏，鄭井元，劉 峰，等. 辣椒 bZIP 家族基因的鑒定與表達(dá)分析[J]. 園藝學(xué)報(bào)， 2018， 45（8）： 1535-1550.

[2] LI D， FU F， ZHANG H， et al. Genome-wide systematic characterization of the ?bZIP ?transcriptional factor family in tomato （ Solanum lycopersicum ?L.） [J]. BMC Genomics， 2015， 16（1）： 771.

[3] WOLFGANG D L， SNOEK B， BEREND S， et al. The ?Arabidopsis bZIP ?transcription factor family-an update [J]. Current Opinion in Plant Biology， 2018， 45： 36-49.

[4] JAKOBY M， WEISSHAAR B， DRGE-LASER W， et al. ?bZIP ?transcription factors in ?Arabidopsis ?[J]. Trends in Plant Science， 2002， 7（3）： 106-111.

[5] YAN Z， XU D， JIA L， et al. Genome-wide identification and structural analysis of ?bZIP ?transcription factor genes in ?Brassica napus ?[J]. Genes， 2017， 8（11）： 288.

[6] WANG Z， ZHU J， YUAN W， et al. Genome-wide characterization of ?bZIP ?transcription factors and their expression patterns in response to drought and salinity stress in ?Jatropha curcas [J]. International Journal of Biological Macromolecules， 2021， 181： 1207-1223.

[7] LI F， LIU J， GUO X， et al. Genome-wide survey， characterization， and expression analysis of ?bZIP ?transcription factors in ?Chenopodium quinoa ?[J]. BMC Plant Biology， 2020， 20（1）： 405.

[8] ZHAO K， CHEN S， YAO W， et al. Genome-wide analysis and expression profile of the ?bZIP ?gene family in poplar [J]. BMC Plant Biology， 2021， 21（122）： 122.

[9] JIN M， GAN S， JIAO J， et al. Genome-wide analysis of the ?bZIP ?gene family and the role of AchnABF1 from postharvest kiwifruit （ Actinidia chinensis ?cv. Hongyang） in osmotic and freezing stress adaptations [J]. Plant Science， 2021， 308： 110927.

[10] CHRISTOPH W， LORENZO P， JEBASINGH S， et al. The Arabidopsis bZIP11 transcription factor links low-energy signalling to auxin-mediated control of primary root growth [J]. PLoS Genetics， 2017， 13（2）： e1006607.

[11] OYAMA T， SHIMURA Y， OKADA K. The Arabidopsis HY5 gene encodes a ?bZIP ?protein that regulates stimulus-induced development of root and hypocotyl [J]. Genes Development， 1997， 11（22）： 2983-2995.

[12] HOU X， HU W， SHEN L， et al. Global identification of ?DELLA ?target genes during ?Arabidopsis ?flower development [J]. Plant Physiology， 2008， 147（3）： 1126-1142.

[13] 徐 偉. 小麥 bZIP 基因 TaGBF 參與植物開花調(diào)控機(jī)制研究[D]. 濟(jì)南： 山東大學(xué)， 2015.

[14] NAN H， CAO D， ZHANG D， et al. GmFT2a and ?GmFT5a ?redundantly and differentially regulate flowering through interaction with and upregulation of the ?bZIP ?transcription factor GmFDL19 in soybean [J]. PLoS One， 2014， 9（5）： e97669.

[15] TAKESHIMA R， NAN H， HARIGAI K， et al. Functional divergence between soybean ?FLOWERING LOCUS T ?orthologues， FT2a and FT5a， in post-flowering stem growth [J]. Journal of Experimental Botany， 2019， 70（15）： 3941-3953.

[16] LI M， HUA W， YUAN Y， et al. Capsella rubella TGA4， a ?bZIP ?transcription factor， causes delayed flowering in ?Arabidopsis thaliana ?[J]. Archives of Biological Sciences， 2015， 68： 141.

[17] 周 波，王 宇，孫 梅，等. 津田蕪菁bZIP蛋白HY5 cDNA的克隆及表達(dá)特性[J]. 分子植物育種， 2008， 6（1）： 59-64.

[18] LIU X， LI Z， HOU Y， et al. Protein interactomic analysis of SAPKs and ABA-inducible ?bZIPs ?revealed key roles of SAPK10 in rice flowering [J]. International Journal of Molecular Sciences， 2019， 20（6）： 1427.

[19] CERISE M， GIAUME F， GALLI M， et al. OsFD4 promotes the rice floral transition via florigen activation complex formation in the shoot apical meristem [J]. New Phytologist， 2020， 229（1）： 429-443.

[20] TONG C， WANG X， YU J， et al. Comprehensive analysis of RNA-seq data reveals the complexity of the transcriptome in ?Brassica rapa ?[J]. BMC Genomics， 2013， 14： 689.

[21] DAI Y， ZHANG S， SUN X， et al. Comparative transcriptome analysis of gene expression and regulatory characteristics associated with different vernalization periods in ?Brassica rapa ?[J]. Genes， 2020， 11（4）： 392.

[22] CHEN G， WANG J， WANG H， et al. Genome-wide analysis of proline-rich extension-like receptor protein kinase （PERK） in ?Brassica rapa ?and its association with the pollen development [J]. BMC Genomics， 2020， 21： 401.

[23] 陳國(guó)戶，王 浩，李 廣，等. 白菜 PRX 基因家族的鑒定與生物信息學(xué)分析[J]. 浙江大學(xué)學(xué)報(bào)， 2020， 46（6）： 677-686.

[24] CHEN C， CHEN H， ZHANG Y， et al. TBtools： an integrative toolkit developed for interactive analyses of big biological data [J]. Molecular Plant， 2020， 13（8）： 1194-1202.

[25] 孫 宇，劉志鑫，葉 子，等. 杧果 RAV 基因家族的全基因組分析[J]. 江蘇農(nóng)業(yè)學(xué)報(bào)， 2021， 37（4）： 957-967.

[26] CANNON S， MITRA A， BAUMGARTEN A， et al. The roles of segmental and tandem gene duplication in the evolution of large gene families in ?Arabidopsis thaliana ?[J]. BMC Plant Biology， 2004， 4： 10.

[27] CHENG F， MANDKOV T， WU J， et al. Deciphering the diploid ancestral genome of the mesohexaploid ?Brassica rapa ?[J]. The Plant Cell， 2013， 25： 1541-1554.

[28] FLAGEL L， WENDEL J. Gene duplication and evolutionary novelty in plants [J]. New Phytologist， 2009， 183（3）： 557-564.

[29] POURABED E， GOLMOHAMADI F， MONFARED P， et al. Basic leucine zipper family in barley： genome-wide characterization of members and expression analysis [J]. Molecular Biotechnology， 2015， 57（1）： 12-26.

（責(zé)任編輯：張震林）

收稿日期：2021-09-10

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目（31801853）;國(guó)家大學(xué)生創(chuàng)新創(chuàng)業(yè)基金項(xiàng)目（202110364080）;安徽省大學(xué)生創(chuàng)新創(chuàng)業(yè)基金項(xiàng)目（S202010364238、S202010364239）;安徽省博士后基金項(xiàng)目（2019B320）

作者簡(jiǎn)介：馬李廣（2000-），男，廣東河源人，學(xué)士，研究方向?yàn)槭卟朔N質(zhì)資源與遺傳育種。（E-mail）1554727563@qq.com。張賀龍為共同第一作者。

通訊作者：陳國(guó)戶，（E-mail）cgh@ahau.edu.cn