薏苡bZIP基因家族全基因組鑒定與非生物脅迫表達分析

周 宇，陳傳敏，王 健，周明強，楊小雨，班秀文，劉凡值，楊成龍*

薏苡bZIP基因家族全基因組鑒定與非生物脅迫表達分析

周 宇1,2，陳傳敏3，王 健2，周明強2，楊小雨2，班秀文2，劉凡值2，楊成龍2*

1. 貴州大學農學院，貴州貴陽 550025；2. 貴州省亞熱帶作物研究所，貴州興義 562400；3. 貴陽市第五中學，貴州貴陽 550003

薏苡是傳統藥食兼用經濟作物，具有極高的營養價值和重要的藥用價值，越來越受人們重視。轉錄因子是一類能夠與調控基因的順式作用元件或者功能基因特異性結合的重要調控蛋白質分子，在植物的許多生物學活動過程中起著重要的調控作用。堿性亮氨酸拉鏈（bZIP）轉錄因子家族在調控植物的生長發育及響應生物與非生物脅迫中起著重要作用。本研究利用北京薏苡（2n=20）基因組信息鑒定出83個bZIP基因家族成員，命名為～，2個bZIP基因（和）不能最終定位到任何連鎖群通過生物信息學方法，分析了基因結構、理化性質、染色體分布，研究了其與其他物種的系統進化關系，并利用轉錄組技術研究了家族成員在栽培種興仁小白殼苗期響應非生物脅迫的基因表達規律。結果表明：薏苡全基因組中bZIP成員分別編碼氨基酸97～1008個，最小的是5（97個氨基酸），而最大的是28（1008個氨基酸）；蛋白質分子量在11.33～110.03 kDa之間，pI在4.15（44）到12.33（50）之間；將其劃分為A～I、K、S等11個亞組，不均勻地分布于基因組的10條染色體上；除了廣泛分布的bZIP結構域的motif 1或motif 2外，每一類都具有相似基因結構與基序；然而，不同成員間存在外顯子數目差異，最多為15個外顯子（），最少為0個；干旱、高溫、干旱和高溫脅迫條件下，83個ClbZIP家族基因在干旱和高溫脅迫后表達量的發生明顯的變化，這說明這些基因可能在干旱和高溫中的調控作用存在著不同。上述結果為更深入研究薏苡bZIP家族基因的生物學功能奠定基礎及其在調控薏苡生長發育與響應生物與非生物脅迫機制提供參考。

薏苡；bZIP轉錄因子；生物信息學；干旱脅迫；高溫脅迫；基因表達

薏苡（L.）又名薏仁米、薏苡仁和六谷子等，為禾本科玉蜀黍族薏苡屬，一年生或多年生草本植物，起源于亞洲[1]，我國西南部也是薏苡屬植物的起源、演化和遷移的初生中心之一[2]。貴州是薏苡的主要栽培區，種植面積和產量均居全國第一，生產量占全國2/3，全省薏苡總產量達到28.8萬t，全省種植和粗加工產值達43億元，已建成全國最大的薏苡加工集聚區和產品集散地[3]。薏苡的營養和藥用價值非常高，不僅被譽為“世界禾本科植物之王”，也被東亞和東南亞許多國家廣泛用作藥材和保健品[4]。薏苡仁中富含蛋白質、脂肪、糖等營養物質，具有較高的食用價值。薏苡整株均可入藥，在醫學方面具有健脾利濕、除脾止瀉、清熱解毒等功效[5]。因此，薏苡作為食物與中草藥越來越受人們重視。

轉錄因子是一類能夠與調控基因的順式作用元件或者功能基因特異性結合的重要調控蛋白質分子，在植物的許多生物活動過程中起著重要的調控作用[6]。MYB、WRKY、AP2/EREBP、bZIP等轉錄因子是目前人們在植物中了解較早且研究較多的幾大類轉錄因子家族。堿性亮氨酸拉鏈（basic leucine zipper, bZIP）轉錄因子家族是目前在植物中發現的64個基因家族中數量最大、最多樣化且高度保守的基因家族之一[7]。bZI轉錄因子不僅參與調控植物的生長發育，而且在植物應對高溫、干旱、病原防御等生物與非生物脅迫中也起著重要調控作用。bZIP轉錄因子家族成員已在多種植物中被鑒定，廣泛存在于擬南芥（75個）[8]、玉米（125個）[9]、水稻（89個）[10]、谷子（85個）[11]等高等植物基因組。JAKOBY等[8]研究發現擬南芥基因通過與ABRE元件相結合，參與ABA信號傳導通路，調控擬南芥植株的生長發育。KANG等[12]發現擬南芥的ABF3和ABF4蛋白通過調控ABA信號，可以增強擬南芥的耐旱性。NIJHAWAN等[10]通過過量表達水稻基因，可以提高水稻抗干旱與耐熱能力。此外，水稻與基因共同過量表達，可以降低轉基因植株的ABA和氯化鈉的敏感性，從而提高水稻的耐受性[13-16]。WEI等[17]研究表明，通過誘導玉米基因表達，可以增強玉米耐鹽性。LIM等[18]研究辣椒被沉默后，通過ABA誘導基因表達、ABA含量改變和關閉氣孔，使其表現出干旱敏感和耐旱性。盧平等[11]研究谷子、和時，發現谷子在干旱處理下的根中表達量較高，參與根干旱脅迫應答。綜上所述，植物bZIP轉錄因子家族在調控植物生長發育和植物響應非生物脅迫并提高其耐受性起著重要的調控作用。

目前栽培薏苡與野生薏苡基因組測序已完成，為薏苡抗性相關基因的挖掘與利用提供了重要的基因組資源。本研究利用栽培薏苡基因組信息和前期研究獲得‘興仁小白殼’轉錄組數據，通過生物信息學方法對薏苡bZIP轉錄因子家族進行全基因組鑒定，并分析了bZIP基因家族在薏苡苗期葉片響應干旱、高溫和干旱高溫脅迫的表達規律，為研究薏苡bZIP轉錄因子基因家族功能奠定基礎，旨在挖掘薏苡抗旱候選基因，為解析薏苡抗旱和耐高溫分子機理提供理論基礎。

1 材料與方法

1.1 ClbZIP基因家族的鑒定

在薏苡中鑒定bZIP蛋白，本研究從NCBI（https://www.ncbi.nlm.nih.gov/bioproject/PRJNA544872）網站下載薏苡基因組序列，以及‘興仁小白殼’干旱、高溫和干旱高溫脅迫轉錄組數據中編碼序列（coding sequences, CDS）。另外，利用Pfam（http://pfam.xfam.org/）蛋白質數據庫下載bZIP轉錄因子的隱馬爾科夫模型文件（PF00170），利用HMMER 3.0軟件中的Hmmsearch程序對薏苡蛋白質序列進行搜索，將得到的結果去除重復，并通過Pfam、SMART和CDD三大數據庫進一步鑒定和篩選，最后得到薏苡bZIP家族蛋白質序列。薏苡bZIP轉錄因子基因的染色體定位和基因組位置信息來源于Phytozome數據庫。然后，對錯誤預測的基因進行人工整理，其中一些基因通過PCR擴增和測序進一步驗證。冗余序列通過手動剔除，經過全面整理，最終在薏苡基因組中鑒定出83個bZIP基因模型（ClbZIP）。使用ExPasy（http://web.expasy.org/protparam/）在線軟件對ClbZIP蛋白的序列長度、分子量、等電點等進行理化性質預測。

1.2 ClbZIP基因家族的結構分析和染色體定位

根據薏苡基因組數據庫中的物理位置信息，利用Mapchart 2.3軟件（https://www.wur.nl/en）將所有ClbZIP基因定位到薏苡染色體上。采用多重共線性掃描工具包（MCScanX）以默認參數分析基因復制事件。為了展示從薏苡和其他選定物種獲得的同源bZIP基因的共線關系，使用Dual-Systeny Plotter軟件構建了共線分析圖。使用TBtools在線軟件繪制基因結構圖；利用BLASTP對轉錄組鑒定的序列與薏苡基因組中的序列比較分析。

1.3 系統進化分析

利用MEGA 7.0軟件把ClbZIP轉錄因子的氨基酸序列與擬南芥、水稻的bZIP蛋白序列進行多序列比對，并用鄰接法（NJ）構建系統進化樹，參數Bootstrap值設為1000，其他參數使用系統默認值。ClbZIP基因結構通過Gene StructureDisplay Server 2.0（http://gsds.cbi.pku.edu.cn/）網站預測，蛋白結構域通過MEME在線工具查詢[19-21]。

1.4 ClbZIP家族成員啟動子順式作用元件分析

從薏苡全基因組數據庫中提取每個薏苡bZIP基因啟動子區域（上游1500 bp），并利用Plant CARE（http://bioinformatics.psb.ugent.be/webtools/ plantcare/html/）網站分析薏苡bZIP基因啟動子的順式作用元件種類、數目及功能。

1.5 基因表達分析

為了研究基因在非生物脅迫下的表達變化規律，本研究利用‘興仁小白殼’作為試驗材料，約1周左右，待材料幼苗生長至3片葉時，選擇大小一致的材料分別進行干旱（15% PEG6000模擬干旱）、高溫（40℃）和干旱高溫（15%PEG6000+40℃）等脅迫處理不同時間（0、3、6、12、24 h）。脅迫后取葉片，每一個樣品分別取3個生物學重復，液氮中速凍后保存于–80℃超低溫冰箱中備用。每一個樣品分別提取RNA構建文庫并利用Illumina測序，經過注釋獲得差異表達基因。借助于Pheatmap對轉錄組數據進行均一化并鑒定基因差異表達，繪制基因表達熱圖，分析響應不同脅迫下基因表達的變化規律。‘興仁小白殼’干旱、高溫和高溫干旱的轉錄組數據已上傳至NCBI數據庫（https://www.ncbi. nlm.nih.gov/），編號為PRJNA812268。

2 結果與分析

2.1 ClbZIP基因家族的鑒定

從Pfam數據庫中獲得的bZIP轉錄因子的隱馬爾科夫模型文件（PF00170），結合bZIP蛋白保守序列，多重比對分析薏苡基因組數據庫，通過剔除不含保守結構域的蛋白質序列，手動去除重復和不完整的序列。最后，基于存在明顯完整的bZIP結構域，選擇83個基因模型并注釋為薏苡bZIP基因。最后，根據其在染色體上的位置依次命名為～（表1），2個ClbZIP基因（和）不能最終定位到任何連鎖群，分別命名為和。

在83個ClbZIP蛋白質中，ClbZIP5被鑒定為最小的蛋白質，含有97個氨基酸，而最大的是ClbZIP28（1008個氨基酸）。蛋白質分子量在11.33～110.03 kDa之間，pI在4.15（44）到12.33（50）之間（表1）。

表1 ClbZIP基因家族理化性質分析

續表1 ClbZIP基因家族理化性質分析

Tab. 1 Analysis of the physicochemical properties of ClbZIPs (continued)

由圖1可知，除了和外，薏苡bZIP基因家族成員不均勻地分布在薏苡10條染色體上，其中，1號染色體含有薏苡bZIP基因最多，有14個基因，占總數的16.87%；9號染色體含有薏苡bZIP基因最少，只有2個。

2.2 系統進化分析

對鑒定得到的83個ClbZIP基因進行基因結構分析，將薏ClbZIP轉錄因子的氨基酸序列進行多序列比對，并利用MEGA7.0軟件構建薏苡、擬南芥和水稻bZIP家族成員構建系統進化樹（圖2）。結果表明，83個薏苡bZIP家族成員物種之間存在差異，與擬南芥親緣關系較遠，與水稻親緣關系較近，并被劃分為A（17個）、B（1個）、C（9個）、D（14個）、E（2個）、F（9個）、G（8個）、H（0個）、I（4個）、K（14個）、S（5個）等11個亞組。除H亞組外，均有ClbZIP基因家族成員分布，A亞組中最多，17個，B、E亞組中分布最少（1個，2個），和水稻（和）聚到一起，其同源性較高，并聚在B亞組，和兩個成員親緣關系較近聚集在E亞組。

2.3 ClbZIP基因家族成員的基因結構和基序組成

對所有已鑒定的ClbZIP基因的外顯子-內含子組成特征進行分析，以進一步了解薏苡中bZIP家族的進化（圖2、圖3A）。結果表明，可劃分為3個大類群，類群Ⅰ有57個基因，類群Ⅱ有7個基因，類群Ⅲ有19個基因；類群Ⅰ的57個ClbZIP基因可進一步劃分為3個亞組，Ⅰa組有34個基因，Ⅰb組有9個基因，Ⅰc組有14個基因。由此可見，多組ClbZIP有較近的親緣關系，如等均為旁系同源基因。

圖1 ClbZIP基因家族的染色體分布

圖2 薏苡與水稻、擬南芥的bZIP基因家族系統進化樹

根據MEME motif分析結果構建了代表所有ClbZIP蛋白結構的示意圖。如圖3A所示，除了廣泛分布的bZIP結構域的motif 1或motif 2外，同一組中的ClbZIP成員通常具有相似的基序組成，除Ⅰa類中有一部分包括motif 1和motif 2。例如，motif 9是I組一部分所特有的，而motif 3、4、5、6是Ⅲ組特有的且大部分含有多個motif。聚集的ClbZIP對，即，顯示出高度相似的基序分布。亞組內ClbZIP蛋白之間相似的基序排列表明，該蛋白結構在亞家族中是相對保守。這些保守基序的大部分功能仍有待闡明。總體而言，同一類群中bZIP成員的保守基序組成和相似的基因結構，以及系統發育分析結果，可以有力地支持類群分類的可靠性。

如圖3B所示，83個ClbZIP基因蛋白質保守基序均有所差異，含有0～7個不等。ClbZIP28中沒有蛋白質保守基序，62個ClbZIP只含有1個保守基序，占薏苡ClbZIP全基因組的75%。Ⅰa組中，只有ClbZIP12含有Motif 2和Motif 7，ClbZIP21ClbZIP40ClbZIP42只含有Motif 2，其余基因均含有Motif 1。Ⅰb組中ClbZIP75ClbZIP8只含有Motif2，其余基因均只含有Motif 1。Ⅰc組中ClbZIP79只含有Motif 2，其余基因均含有Motif 2和Motif 7。類群Ⅱ中所有基因均只含一個保守基序，其中只有ClbZIP5含有Motif7，其余基因只含Motif 2。ClbZIP69ClbZIP77ClbZIP9ClbZIP27ClbZIP39ClbZIP48ClbZIP54ClbZIP76ClbZIP81ClbZIP30含有保守基序最多，有7個，均存在于類群Ⅲ中，并且ClbZIP71ClbZIP45ClbZIP17ClbZIP58均只含Motif 2。親緣關系較近的基因蛋白質保守基序組成較相似，例如ClbZIP3/ ClbZIP4ClbZIP52/ClbZIP72均只含有Motif 2，且位置相近。

從圖3C可以看出，ClbZIP基因家族成員的每個基因長度、內含子與外顯子個數都有所不同，其內含子個數在0～15個之間，其中內含子數量最多的ClbZIP基因是有15個。均無內含子。

A：ClbZIP家族成員的系統進化樹，不同顏色的枝代表不同的亞族；B：ClbZIP蛋白的基序組成，不同的基序用不同顏色的方框表示；C：ClbZIP基因的結構，框表示外顯子，黑線代表內含子。

2.4 ClbZIP家族成員的染色體分布及共線性分析

圖4所示，ClbZIP家族成員基因在10個薏苡連鎖群（LG）上的分布不均勻。LG1、LG2、LG4含有最多的ClbZIP基因。一些連鎖群（如LG1、LG2、LG4）有更多的基因，而其他連鎖群則很少；有些LGs只有一個基因（例如LG9）。LG長度與ClbZIP基因數無顯著正相關。根據200 kb內包含2個或更多基因的染色體區域被定義為串聯復制事件，通過分析ClbZIP的基因復制事件，ClbZIP基因家族中不存在串聯復制現象，有29對大片段復制基因。

為了進一步推斷ClbZIP家族的系統發育機制，我們構建了2個代表物種相關的薏苡比較共線圖，包括1個雙子葉植物（擬南芥）和1個單子葉植物（水稻）（圖5）。在薏苡與擬南芥的bZIP基因中僅有一對基因存在共線性關系，而薏苡與水稻的bZIP基因中存在52對基因存在共線性關系。在薏苡和擬南芥之間發現的共線基因對，在薏苡和水稻之間沒有發現，這可能表明這些同源基因對是在雙子葉植物和單子葉植物分化之后形成的。

chr1～chr10表示薏苡的1～10號染色體；灰色部分表示薏苡全基因組的共線性基因；紅色線條表示薏苡bZIP基因家族中的共線性。

數字表示不同物種中的染色體；灰色部分表示2個物種間全基因組的共線性基因；紅色線條表示2個物種間bZIP基因家族中的共線性。

2.5 ClbZIP基因家族啟動子順式作用元件分析

將ClbZIP基因起始密碼子上游1500 bp序列用于順式作用元件分析，如圖6所示。ClbZIP基因啟動子區域富含多個與ABA響應元件（ABRE）、植物光周期響應相關的順式作用元件（AE-box、Box）。此外該基因家族還有與植物激素調節相關的順式作用元件及與逆境響應相關的順式作用元件。此結果表明ClbZIP基因在應對生物與非生物脅迫中起重要的調控作用。

圖6 ClbZIP基因家族順式作用元件

2.6 ClbZIP基因家族表達分析

根據本研究‘興仁小白殼’干旱、高溫和干旱高溫脅迫轉錄組中的FPKM數據，建立基因表達熱圖（圖7、圖8、圖9），以研究ClbZIP家族基因在干旱和高溫脅迫中所起的作用。從熱圖中發現，83個ClbZIP家族基因在干旱和高溫脅迫后的表達量發生明顯的變化，這說明這些基因可能在干旱和高溫中的調控作用存在著不同。

如圖7所示，在干旱脅迫3 h后，的表達量較高；在干旱6 h脅迫后，的表達量較高；在干旱脅迫12 h后，的表達量較高；在干旱脅迫24 h后，等基因具有較高的表達量。

在高溫脅迫3 h后，的表達量較高；在高溫脅迫6 h后，的表達量較高；在高溫脅迫12 h后，的表達量較高；在高溫脅迫24 h后，的表達量較高（圖8）。

在干旱和高溫共同脅迫3 h后，的表達量較高；在干旱和高溫共同脅迫6 h和12 h下有較高的表達量；干旱和高溫共同脅迫24 h后，的表達量較高（圖9）。

熱圖bar為log2；CK為未經干旱處理；Drought-3 h為干旱處理3 h；Drought-6 h為干旱處理6 h；Drought-12 h為干旱處理12 h；Drought-24 h為干旱處理24 h。

熱圖bar為log2；CK為未經高溫處理，Heat-3 h為高溫處理3 h；Heat-6 h為高溫處理6 h；Heat-12 h為高溫處理12 h；Heat-24 h為高溫處理24 h。

3 討論

bZIP基因家族普遍存在于植物中，其在植物生長、發育和抗逆等方面中起著重要作用，是最大的轉錄因子家族之一[22]。目前，bZIP基因家族已在多種植物中被鑒定出，但在薏苡的研究中尚未見報道。本研究通過對薏苡基因組數據中的蛋白質序列及興仁小白殼的轉錄組數據進行分析，去除冗余后鑒定83個bZIP轉錄因子家族成員。本研究對鑒定的83個ClbZIP家族成員在染色體位置分析，發現ClbZIP基因不均勻地分布在1～10號染色體上。目前在與薏苡同為禾本科的谷子與水稻中分別鑒定出89和85個bZIP基因，其數量上相近。薏苡bZIP基因具有一定的保守性，親緣關系相近的蛋白保守基序與基因結構也較為相似。薏苡bZIP轉錄因子家族的每個基因長度、內含子與外顯子個數都有所差異，基因長度在97～ 1008個氨基酸之間，內含子個數最多的有15個，結構存在著復雜性，從而說明結構與功能存在關聯性，復雜生物功能分化出不同類型基因結構。根據前人研究發現多種作物的bZIP家族基因中均普遍存在沒有內含子的bZIP基因，如水稻、木薯、西洋梨，它們分別有15.3%、20.8%、28.8%的bZIP基因沒有內含子[23-25]。基因復制事件在基因的進化過程中起著重要作用[26-27]。薏苡bZIP基因中存在29對大片段復制基因，基因片段復制是基因家族產生的重要因素，在大豆bZIP基因中有58%的基因存在片段復制現象[28]。薏苡與水稻的bZIP基因中存在52對基因存在共線性關系，而與擬南芥僅有1對基因存在共線性關系，說明薏苡與水稻親緣關系較近。

前人研究發現，藜麥bZIP基因在不同組織部位均有表達[29]。谷子在暗處理后和的表達量較高[11]。水稻在干旱、PEG、ABA處理后具有較高的表達量[30-31]。苧麻在高鹽脅迫下的表達量較高[32]。擬南芥正調控鹽誘導基因，參與鹽信號轉導途徑，并進行級聯反應調控擬南芥的響應鹽脅迫[33]。擬南芥和屬于bZIP家族的TGA亞族，在抗逆過程中起著重要作用[34-35]。通常認為，聚類越近越相似，功能和結構上相似度越高。因此，推測擬南芥中的聚為一類的和基因可能參與薏苡抗旱的過程，參與調控激素GA和ABA的代謝過程，是否參與類似功能有待進一步驗證。一般認為，基因表達量與其功能存在一定相關性，在轉錄組中發現薏苡bZIP家族多種基因在不同環境條件脅迫與不同脅迫時間條件下有不同程度的表達量，并且發現薏苡bZIP基因啟動子含有多種調控生長發育及響應生物與非生物脅迫的順式作用元件，但其是否在薏苡抗逆過程中起到重要的作用還需進一步驗證。

熱圖bar為log2；CK為未經處理；DH-3 h為干旱和高溫處理3 h；DH-6 h為干旱和高溫處理6 h；DH-12 h為干旱和高溫處理12 h；DH-24 h為干旱和高溫處理24 h。

薏苡作為一種糧藥兼用作物，具有較高的營養與藥用價值，但目前對其研究起步較晚并且尚未深入。bZIP轉錄因子是植物生長發育過程中重要的調控基因，并且在植物響應逆境環境脅迫中起重要的調控作用。尤其借助轉錄組與轉錄因子對薏苡抗性調控機制的研究，報道比較少。因此，對薏苡bZIP基因家族的研究是及有必要的，但在薏苡中研究抗逆性遠沒有擬南芥、水稻中的深入和清楚。本研究利用北京薏苡基因組和‘興仁小白殼’干旱和高溫脅迫的轉錄組對薏苡bZIP基因家族鑒定及干旱和高溫脅迫后bZIP基因表達規律的研究，旨在為更深入研究薏苡bZIP基因家族的生物學功能奠定基礎及其在調控薏苡生長發育與響應生物與非生物脅迫機制提供參考。

[1] ARORA R K. Job’s tears (L.) a minor food and fodder crop of Northeastern India[J]. Economic Botany, 1977, 31(3): 358-366.

[2] 陸 平, 左志明. 廣西水生薏苡種的發現與鑒定[J]. 廣西農業科學, 1996, 1: 18-20.

LU P, ZUO Z M. Discovery and identification ofRoxb. in Guangxi[J]. Guangxi Agricultural Science, 1996, 1: 18-20. (in Chinese)

[3] 李祥棟, 潘 虹, 陸秀娟, 魏心元, 陸 平, 石 明, 秦禮康. 薏苡種質的主要營養組分特征及綜合評價[J]. 中國農業科學, 2018, 51(5): 835-850.

LI X D, PAN H, LU X J, WEI X Y, LU P, SHI M, QIN L K. Characteristics and comprehensive assessment of principal nutritional components in adlay landraces[J]. Scientia Agricultura Sinica, 2018, 51(5): 835-850. (in Chinese)

[4] 汪 燦, 周棱波, 張國兵, 張立異, 徐 燕, 高 旭, 姜 訥, 邵明波. 薏苡種質資源萌發期抗旱性鑒定及抗旱指標篩選[J]. 植物遺傳資源學報, 2017, 18(5): 846-859.

WANG C, ZHOU L B, ZHANG G B, ZHANG L Y, XU Y, GAO X, JIANG N, SHAO M B. Identification and indices screening of drought resistance in job’s tears germplasm resources at germination stage[J]. Journal of Plant Genetic Resources, 2017, 18(5): 846-859. (in Chinese)

[5] 汪 燦, 周棱波, 張國兵, 張立異, 徐 燕, 高 旭, 姜 訥, 邵明波. 薏苡種質資源苗期抗旱性鑒定及抗旱指標篩選[J]. 中國農業科學, 2017, 50(15): 2872-2887.

WANG C, ZHOU L B, ZHANG G B, ZHANG L Y, XU Y, GAO X, JIANG N, SHAO M B. Drought resistance identification and drought resistance indices screening of job’s tears (L.) germplasm resources at seedling stage[J]. Scientia Agricultura Sinica, 2017, 50(15): 2872-2887. (in Chinese)

[6] 王金英, 丁 峰, 潘介春, 張樹偉, 楊亞涵, 黃 幸, 范志毅, 李 琳, 王 穎. 植物轉錄因子家族的研究進展[J]. 熱帶農業科學, 2019(6): 39-45.

WANG J Y, DING F, PAN J C, ZHANG S W, YANG Y H, HUANG X, FAN Z Y, LI L, WANG Y. Research progress oflineage transcription factors in plant[J]. Chinese Journal of Tropical Agriculture,2019(6): 39-45. (in Chinese)

[7] PAULINO P R, MAURICIO R, GUEDES C, RENSING S A, BIRGIT K, BERND M R. PlnTFDB: updated content and new features of the plant transcription factor database[J]. Nucleic Acids Research, 2010, 38(suppl 1): 822-827.

[8] JAKOBY M, WEISSHAAR B, DRGE-LASER W, VICENTE-CARBAJOSA J, TIEDEMANN J, KROJ T, PARCY F, GROUP B.transcription factors in[J]. Trends in Plant Science, 2002, 7(3): 106-111.

[9] 于 滔, 王成波, 曹士亮, 孫培元. 玉米轉錄因子的生物信息學分析[J]. 黑龍江農業科學, 2016(4): 1-5.

YU T, WANG C B, CAO S L, SUN P Y. Bioinformatics analysis of thetranscription factor in maize[J]. Heilongjiang Agricultural Sciences, 2016(4): 1-5. (in Chinese)

[10] NIJHAWAN A, JAIN M, TYAGI A K, KHURANA J P. Genomic survey and gene expression analysis of the basic leucine zipper transcription factor family in rice[J]. Plant Physiology.2008, 146(2): 333-350

[11] 盧 平, 武懿茂, 武強強, 李雪垠. 谷子轉錄因子家族的全基因組鑒定與生物信息學分析[J]. 山西農業科學, 2020, 48(9): 1361-1370, 1430.

LU P, WU Y M, WU Q Q, LI X Y. Genome-wide identification and bioinformatics analysis oftranscription factor family[J]. Journal of Shanxi Agricultural Sciences, 2020, 48(9): 1361-1370, 1430.

[12] KANG J Y, CHOI H I, KIM S Y, IM M Y.basic leucine zipper proteins that mediate stress responsive abscisic acid signaling[J].Plant Cell, 2002, 14: 343-357.

[13] ZONG W, TANG N, YANG J, PENG L, MA S, XU Y, LI G, XIONG L. Feedback regulation of ABA signaling and biosynthesis by atranscription factor targets drought- resistance-related genes[J]. Plant Physiology, 2016, 171(4): 2810-2825.

[14] PANDEY A S, SHARMA E, JAIN N, SINGH B, BURMAN N, KHURANA J P. A ricetranscription factor, Os16, regulates abiotic stress tolerance when over-expressed in[J].Journal of Plant Biochemistry & Biotechnology, 2018, 27(4): 393-400.

[15] 李紅婷. 擬南芥轉錄因子靶基因的篩選及鑒定[D]. 合肥: 安徽農業大學, 2020.

LI H T. Screening and identification of target genes ofin[D]. Hefei: Anhui Agricultural University, 2020. (in Chinese).

[16] 孫陽陽. 水稻轉錄因子Os18和OsBCAT家族基因的功能研究[D]. 武漢: 華中農業大學, 2020.

SUN Y Y. Function alanalysis of rice transcription factor Os18 and OsBCAT family genes[D]. Wuhan: Huazhong Agricultural University, 2020. (in Chinese).

[17] WEI K F, CHEN J , WANG Y M, CHEN Y H, CHEN S X, LIN Y N, PAN S, ZHONG X J, XIE D X. Genome-wide analysis of-Encoding genes in maize[J]. DNA Research, 2012, 19(6): 463-476.

[18] LIM C W, BAEK W, LEE S C. Roles of pepperprotein CaDILZ1 and its interacting partner RING-type E3 ligase CaDSR1 in modulation of drought tolerance[J].Plant Journal, 2018, 96(2): 452-467.

[19] JIN Z, WEI X, LIU A. Genomic surveys and expression analysis ofgene family in castor bean (L.)[J].Planta, 2014, 239(2).

[20] WANG M , YUAN F, HAO H , ZHANG Y, ZHAO H, GUO A, HU J, ZHOU X, XIE C G. BolOST1, an ortholog of open stomata 1 with alternative splicing products in, positively modulates drought responses in plants[J]. Biochemical & Biophysical Research Communications, 2013, 442(3/4): 214-220.

[21] YUAN F, WANG M, HAO H, ZHANG Y, ZHAO H, GUO A, HU J, ZHOU X, XIE C G. Negative regulation of abscisic acid signaling by theABI1 ortholog[J]. Biochemical & Biophysical Research Communications, 2013, 442(3/4): 202-208.

[22] DEPPMANN C D, TAPAROWSKY A E J. Cross-species annotation of basic leucine zipper factor interactions: insight into the evolution of closed interaction networks[J]. Molecular Biology and Evolution,2006, 23(8): 1480.

[23] HU W, YANG H, YAN Y, WEI Y, TIE W, DING Z, ZUO J, PENG M, LI K. Genome-wide characterization and analysis oftranscription factor gene family related to abiotic stress in cassava[J]. Scientific Reports, 2016, 6: 22783.

[24] 孫耀國, 蔡天潤, 姬行舟, 張 軍. 西洋梨全基因組基因家族生物信息學分析[J]. 林業與生態科學, 2021, 36(1): 24-34.

SUN Y G, CAI T R, JI X Z, ZHANG J. Genome-wide bioinformatics analysis ofgene family in[J]. Forestry and Ecological Sciences, 2021, 36(1): 24-34. (in Chinese)

[25] 金曉芬. 水稻類轉錄因子的功能分析[D]. 揚州: 揚州大學, 2008.

JIN X F. Functional analysis of, atranscription factor in rice (L)[D]. Yangzhou: Yangzhou University, 2008. (in Chinese).

[26] FLAGEL L E, WENDEL J F. Gene duplication and evolutionary novelty in plants[J]. New Phytologist, 2009, 183(3): 557-564.

[27] MOORE R C, MD P. The evolutionary dynamics of plant duplicate genes[J]. Current Opinion in Plant Biology, 2005, 8(2): 122-128.

[28] ZHANG M, LIU Y, SHI H, GUO M, CHAI M, HE Q, YAN M, CAO D, ZHAO L, CAI H. Evolutionary and expression analyses of soybean basic leucine zipper transcription factor family[J]. BMC Genomics, 2018, 19(1): 159

[29] LI F, LIU J, GUO X, YIN L, WEN R. Genome-wide survey, characterization, and expression analysis oftranscription factors in[J]. BMC Plant Biology, 2020, 20(1): 405.

[30] HAO C, WEI C, ZHOU J, HUANG H, CHEN L, CHEN H, XING W D. Basic leucine zipper transcription factor Os16 positively regulates drought resistance in rice[J]. Plant Science, 2012, 193/194: 8-17.

[31] HOSSAIN M A, CHO J I, HAN M, AHN C H, JEON J S, AN G, PARK P B. The ABRE-bindingtranscription factor OsABF2 is a positive regulator of abiotic stress and ABA signaling in rice[J]. Journal of Plant Physiology, 2010, 167(17): 1512-1520.

[32] HUANG C, ZHOU J, JIE Y, XING H, ZHONG Y, YU W, SHE W, MA Y, LIU Z, ZHANG Y. A ramietranscription factor Bn2 is involved in drought, salt, and heavy metal stress response[J]. DNA & Cell Biology, 2016, 35(12): 776-786.

[33] LIU J X, SRIVASTAVA R, CHE P, HOWELL S H. Salt stress responses inutilize a signal transduction pathway related to endoplasmic reticulum stress signaling[J]. The Plant Journal, 2007, 51(5): 897-909.

[34] MIAO Z H, LIU X J, LAM E. TGA3 is a distinct member of the TGA family oftranscription factors in[J]. Molecular Biology Reports, 1994, 25(1): 1-11.

[35] ZANDER M, LA CAMERA S, LAMOTTE O.class-II TGA transcription factors are essential activators of jasmonic acid/ethylene-induced defense responses[J]. The Plant Journal, 2010, 61(2): 200-210.

Genome-wide Identification and Expression Analysis of bZIP Gene Family under Abiotic Stress inL.

ZHOU Yu1,2, CHEN Chuanmin3, WANG Jian2, ZHOU Mingqiang2, YANG Xiaoyu2, BNA Xiuwen2, LIU Fanzhi2, YANG Chenglong2*

1. College of Agriculturre, Guizhou University, Guiyang, Guizhou 550025, China; 2. Guizhou Institute of Subtropical Crops, Xingyi, Guizhou 562400, China; 3. Guiyang Fifth Middle School, Guiyang, Guizhou 550003, China

Coix is a traditional economic crop that can be used both as food and medicine, and has high nutritional and medicinal value. As a kind of food and traditional Chinese medicine, coix is getting more and more attention. Transcription factors are important regulatory protein molecules that can specifically bind to cis-acting elements or functional genes of regulatory genes and play an important regulatory role in many biological processes of plants. The basic leucine zipper (bZIP) transcription factor family plays an important role in regulating plant growth and development and responding to biotic and abiotic stresses. In this study, 83 bZIP gene family members named1–83 were identified using the genomic information of Beijing Coix (2n=20). Twogenes (and) could not be located to any linkage group. Through bioinformatics methods, the gene structure, physical and chemical properties and chromosome distribution were analyzed, and its phylogenetic relationship with other species was studied. The gene expression rules of family members in response to abiotic stress at the seedling stage of Xingren small white shell were studied by the transcriptome technology. The results showed that the bZIP members in coix’s genome encoded 97–1008 amino acids, respectively. The smallest was(97 amino acids), and the largest was(1008 amino acids). The molecular weight of the protein was between 11.33–110.03 kDa, the pI was between 4.15 () and 12.33 (); They were divided into 11 subgroups such as A–I, K and s, which were unevenly distributed on 10 chromosomes of the genome. Except motif 1 or motif 2 of the widely distributed bZIP domain, each class had similar gene structures and motifs. However, there were differences in the number of exons among different members, with a maximum of 15 exons () and a minimum of 0. Under drought, high temperature, drought and high temperature stress, the expression of 83gene families changed significantly after drought and high temperature stress, indicating that the genes may have different regulatory roles in drought and high temperature. The above results would lay a foundation for further understanding the biological functions of the bZIP gene family inand provide reference for the regulation of growth and development and response to biological and abiotic stresses, and provide theoretical guidance for the study of the regulation of the bZIP gene family in drought-resistant and high-temperature environments.

L.; bZIP transcription factor; bioinformatics analysis; drought stress; high temperature stress; gene expression

S519

A

10.3969/j.issn.1000-2561.2022.10.006

2022-02-15；

2022-03-04

國家自然科學基金項目（No. 31701459）；貴州省農業科學院院專項（黔農科院院專項(2016)019號）；貴州省農業科學院青年基金項目（黔農科院青年科技基金[2022]22號）。

周 宇（1997—），男，碩士研究生，研究方向：作物遺傳育種。*通信作者（Corresponding author）：楊成龍（YANG Chenglong），E-mail：yangchenglong208@163.com。