摘要" 為了解蓖麻bZIP基因家族成員的組成和家族進化關(guān)系及特征，本研究利用生物信息學(xué)鑒定家族成員，進行基因結(jié)構(gòu)、蛋白質(zhì)保守基序分析及系統(tǒng)發(fā)育分析，同時分析其蛋白的理化性質(zhì)、磷酸化位點和互作關(guān)系，以及該基因在組織中的表達模式和表達量。結(jié)果表明，共鑒定出46個蓖麻bZIP基因家族成員，位于10個染色體上；基因結(jié)構(gòu)與蛋白質(zhì)保守基序較為復(fù)雜，可分為11個亞族，包含25個基序；與擬南芥進化相似，較為保守；bZIP蛋白為親水性蛋白，位于細胞核、細胞膜和葉綠體，主要通過絲氨酸的磷酸化來發(fā)揮生物學(xué)作用，蛋白RcbZIP18、RcbZIP27、RcbZIP2、RcbZIP17和RcbZIP46之間的聯(lián)系較多，通過蛋白間的相互作用參與蓖麻各個生長時期；該基因家族均參與葉、莖、根和胚等生長時期的發(fā)育調(diào)控，其中RcbZIP10、RcbZIP25和RcbZIP44在葉片和種子中表現(xiàn)出明顯的組織差異性。上述結(jié)果為蓖麻bZIP基因家族成員的功能研究提供參考。

關(guān)鍵詞" 蓖麻；bZIP家族；生物信息學(xué)分析；系統(tǒng)進化；表達分析

中圖分類號" S565.6" " " "文獻標(biāo)識碼" A" " " "文章編號" 1007-7731（2025）02-0027-09

DOI號" 10.16377/j.cnki.issn1007-7731.2025.02.006

Genome-wide identification， bioinformatics and expression analysis of bZIP gene family of Ricinus communis

WANG Dongyan1" " ZHANG Chunling2" " FENG Yan1" " WANG Lixin1" " ZAN Rui1" " DI Jianjun1，3，4，5，6

（1College of Life Science and Food Engineering， Inner Mongolia Minzu University， Tongliao 028000， China;

2Tongliao Agricultural and Animal Husbandry Academy， Tongliao 028015， China;

3Inner Mongolia Key Laboratory of Castor Breeding and Comprehensive Utilization， Tongliao 028000， China;

4Key Laboratory of Castor Breeding of the State Et hnic Affairs Commission， Tongliao 028000， China;

5Inner Mongolia Industrial Engineering Research Center of Universities for Castor， Tongliao 028000， China;

6Inner Mongolia Engineering Research Center of Industrial Technology Innovation of Castor， Tongliao 028000， China）

Abstract" To understand the composition， evolutionary relationship and characteristics of Ricinus communis bZIP gene family members， bioinformatics was used to identify family members， analyze gene structure， protein conserved motifs， phylogenetic analysis， and analysis of protein physicochemical properties， phosphorylation sites， and interactions， as well as gene expression patterns and expression levels in tissues. The results showed that 46 members of the Ricinus communis bZIP gene family were identified， which were located on 10 chromosomes.The gene structure and protein conserved motifs were complex， which can be divided into 11 subfamilies， including 25 motifs. Similar to Arabidopsis， the evolution was more conservative; bZIP protein was a hydrophilic protein located in the nucleus， cell membrane and chloroplast. It played a biological role mainly through the phosphorylation of serine.The proteins RcbZIP18， RcbZIP27， RcbZIP2， RcbZIP17 and RcbZIP46 were associated with each other， participating in each growth cycle of Ricinus communis through the interaction between proteins. They are involved in the developmental regulation of leaf， stem， root and embryo during growth，RcbZIP10， RcbZIP25 and RcbZIP44 showed significant organizational differences in leaves and seeds These results provide references for the functional study of Ricinus communis bZIP gene family members .

Keywords" Ricinus communis L.; bZIP family; bioinformatics analysis; phylogeny; expression analysis

轉(zhuǎn)錄因子（Transcription factor，TF）是一類在特定的時間與空間下，目的基因與特定的序列進行專一性結(jié)合的蛋白質(zhì)分子[1]。其作為細胞過程（如分化和發(fā)育）的關(guān)鍵調(diào)節(jié)因子，通過信號通路對外部干擾做出反應(yīng)，在植物生長發(fā)育、代謝反應(yīng)和抗逆反應(yīng)等多種過程中發(fā)揮重要作用[2-3]。堿性亮氨酸拉鏈（Basic leucine zipper，bZIP）是真核生物中較為熟知的轉(zhuǎn)錄因子家族之一[4]。bZIP基因家族廣泛存在于植物中，在多種生物過程中起著重要作用，如孫耀國等[5]研究表明，西洋梨的bZIP基因家族參與該植物的生長發(fā)育、種子成熟、光信號響應(yīng)和環(huán)境脅迫等過程；曹麗茹等[6]研究發(fā)現(xiàn)，玉米bZIP轉(zhuǎn)錄因子ZmbZIP26在多個部位表達，且在氮脅迫和高鹽條件下表達量上調(diào)20倍左右，在高溫和干旱脅迫下表達量也呈上調(diào)趨勢；潘婷婷[7]研究表明，在水稻中共鑒定出89個bZIP基因，其中調(diào)控水稻開花的新基因OsbZIP65在多種部位表達，在根、莖和穗中表達量較高，在葉片等部位表達量較低。游冉[8]研究發(fā)現(xiàn)，過表達OsbZIP34基因能夠延遲植物開花，并參與調(diào)控水稻抽穗。目前，在野菊、水稻、玉米、大豆和睡蓮等多種植物中，已有bZIP轉(zhuǎn)錄因子的基因組鑒定及生物學(xué)分析[9]，但關(guān)于bZIP轉(zhuǎn)錄因子在蓖麻（Ricinus communis）中的分子進化研究相對較少。

蓖麻是大戟科油料植物[10]，其種子含油量較高，蓖麻油中富含的蓖麻油酸在生活與工業(yè)中有著廣泛的用途[11-12]，為再生能源提供優(yōu)質(zhì)原料[13]。蓖麻油及其衍生物可代替部分石油產(chǎn)品，被稱為能夠再生的“綠色石油”[14]。然而，受到環(huán)境和種植效益的影響，蓖麻產(chǎn)量逐年下降，使產(chǎn)品供不應(yīng)求。因此，利用基因工程和生物技術(shù)改良蓖麻的產(chǎn)量與性狀對生產(chǎn)具有重要意義[15]。

本研究主要通過生物信息學(xué)技術(shù)對蓖麻bZIP轉(zhuǎn)錄因子家族的蛋白和基因的保守結(jié)構(gòu)、系統(tǒng)進化及其基因表達模式等進行分析，以了解bZIP轉(zhuǎn)錄因子家族在蓖麻中的進化過程[16]，為蓖麻bZIP轉(zhuǎn)錄因子家族在油脂代謝調(diào)控機制中的功能研究提供參考，有利于提高蓖麻育種效率和培育高產(chǎn)優(yōu)質(zhì)抗逆品種，對蓖麻穩(wěn)產(chǎn)增產(chǎn)有重要作用。

1 材料與方法

1.1 數(shù)據(jù)來源

使用的蓖麻全基因組序列、編碼序列、蛋白質(zhì)序列和注釋文件均下載自網(wǎng)站Oil plant database（https：//oilplants.iflora.cn），蓖麻基因組文件下載自蓖麻基因組數(shù)據(jù)庫 （https：//castorbean.jcvi.org./index.php），擬南芥bZIP基因家族的基因與蛋白質(zhì)序列下載自TAIR數(shù)據(jù)庫（Araport11，https：//www.arabidopsis.org/）[17]，蓖麻的各組織相關(guān)表達量數(shù)據(jù)下載自GEO數(shù)據(jù)庫（https：//www.ncbi.nlm.nih.gov）。

1.2 家族成員鑒定

下載蓖麻的蛋白質(zhì)序列和70個擬南芥bZIP基因家族成員序列，然后進行本地BLASTP比對（E-value≤1e-10）[18]。再從Pfam數(shù)據(jù)庫（http：//pfam.xfam.org/）下載bZIP基因家族典型結(jié)構(gòu)域的蛋白質(zhì)序列，進行本地HMMER比對[19]。采用這兩種方法進行家族成員的鑒定，將比對結(jié)果進行整合，刪除多余的重復(fù)序列。分別使用NCBI CDD工具（https：//www.ncbi.nlm.nih.gov/cdd/）與EBI InterPro工具（https：//www.ebi.ac.uk/interpro/）對候選序列進行結(jié)構(gòu)域檢測[20-21]，使用Clustal Omega軟件進行多重比對[22-23]，從而確定蓖麻bZIP基因家族的成員。

1.3 基因結(jié)構(gòu)與蛋白質(zhì)保守基序分析

根據(jù)下載的蓖麻基因組注釋文件編寫腳本，獲得bZIP基因家族成員的基因、內(nèi)含子與外顯子位點等信息。使用MEME在線工具（https：//meme-suite.org/tools/meme）對蓖麻bZIP蛋白進行保守基序分析[24]，基序數(shù)目設(shè)置為10個。使用Evoview在線工具對基因結(jié)構(gòu)、蛋白質(zhì)保守基序進行可視化分析[25]。

1.4 系統(tǒng)發(fā)育進化分析

在植物轉(zhuǎn)錄因子數(shù)據(jù)庫（http：//planttfdb.gao-lab.org/）下載模式植物擬南芥bZIP轉(zhuǎn)錄因子的家族成員基因序列。運用MEGA7軟件對擬南芥和蓖麻bZIP基因家族成員采用Clustal W方法進行多序列比對[26]。選用相鄰連接法（Neighbor-joining，NJ），設(shè)定Bootstrap值1 000進行校正，構(gòu)建進化樹。導(dǎo)出Newick Export文件。用EvolView（https：//evolgenius.info//evolview-v2/）進行進化樹美化[27]。

1.5 蛋白的理化性質(zhì)分析

利用在線工具ExPASy的ProtParam工具（https：//web.expasy.org/protparam/）對基因家族成員蛋白的氨基酸數(shù)目、分子量、理論等電點和親水性等理化性質(zhì)進行分析[28]。在Plant-mPLoc[25]中進行蛋白質(zhì)的亞細胞定位分析。

1.6 磷酸位點預(yù)測

利用在線網(wǎng)站NetPhos3.1（NetPhos 3.1-DTU Health Tech -Bioinformatic Services /netphos-3.1/）對潛在的蛋白磷酸化位點數(shù)量進行預(yù)測[29]，主要為蘇氨酸（Threonine，Thr）、酪氨酸（Tyrosine，Tyr）和絲氨酸（Serine，Ser）。

1.7 蛋白互作網(wǎng)絡(luò)分析

參考蓖麻蛋白質(zhì)數(shù)據(jù)庫，利用在線網(wǎng)站String （https：//cn.string-db.org/）對蓖麻bZIP家族成員蛋白進行蛋白質(zhì)網(wǎng)絡(luò)互作建模，從而預(yù)測該基因家族蛋白質(zhì)間的相互作用[30]。

1.8 表達模式分析

在NCBI GEO數(shù)據(jù)庫（https：//www.ncbi.nlm.nih.gov/geo/）獲取蓖麻在Ⅱ/Ⅲ期胚乳、Ⅴ/Ⅵ期胚乳、萌芽種子、葉片和雄花5個組織中基因的表達量[28]。得到的reads count 值進行FPKM標(biāo)準(zhǔn)化處理。對蓖麻的bZIP基因家族進行層次聚類，使用R4.1.0軟件（https：//www.r-project.org/）的Pheatmap包繪制基因表達量熱圖。

1.9 熒光定量分析

選擇具有代表性的差異基因，測定其在蓖麻種子和葉片中的表達量。提取蓖麻種子和葉片RNA，使用反轉(zhuǎn)錄試劑盒PrimeScript TMRT Master Mix（TaKaRa）反轉(zhuǎn)錄合成cDNA。用TB Green Premix Ex TaqTM Ⅱ試劑（TaKaRa）進行熒光定量PCR試驗，待測樣本設(shè)3個重復(fù)。參照2-ΔΔCt方法計算bZIP在蓖麻不同組織基因的相對表達量。

2 結(jié)果與分析

2.1 家族成員鑒定

由表1可知，通過比對鑒定共得到46個蓖麻bZIP基因家族成員，按照基因在染色體上的定位，命名為RcbZIP1～RcbZIP46，基因長度在476～19 184 bp。對蓖麻基因組注釋文件進行分析，發(fā)現(xiàn)46個蓖麻bZIP基因分布在10條染色體上。其中，2號染色體上分布的bZIP基因最多，有12個；3和6號染色體分布較少，各2個基因；1、4、5、8和9號染色體上均分布4個基因；7和10號染色體上各分布5個基因。

2.2 蓖麻bZIP基因結(jié)構(gòu)與蛋白質(zhì)保守基序

由表2可知，46個蓖麻bZIP基因家族成員分為11個亞族，第Ⅶ亞組包含了9個基因，是包含基因最多的亞族，Ⅲ組、Ⅷ組、Ⅸ組和Ⅹ組bZIP基因最少，各組僅有2個基因。除了Ⅰ組和Ⅱ組的部分基因無內(nèi)含子外，其余bZIP基因家族成員均含有內(nèi)含子和外顯子，內(nèi)含子數(shù)0～25個，外顯子數(shù)1～26個；其中，Ⅱ組和Ⅲ組組內(nèi)基因內(nèi)含子數(shù)和外顯子數(shù)完全一致。

基序是介于生物大分子二級和三級結(jié)構(gòu)之間的一種結(jié)構(gòu)層次，也被稱為模體（Motif）[31]。11個亞族共25個Motif，單個成員含有的Motif范圍在2～10個，其中Motif1為所有基因共有，經(jīng)比較發(fā)現(xiàn)，Motif1是典型的類鋅指結(jié)構(gòu)CX2CX6CX3C，是bZIP基因家族的標(biāo)志；Ⅰ～Ⅴ組內(nèi)家族成員同時具有Motif1和Motif3，Ⅵ～Ⅺ組內(nèi)家族成員同時具有Motif1和Motif2；其中，Ⅴ組組內(nèi)家族成員Motif種類及排布完全一致，Ⅱ組除了bZIP11外，其他家族成員Motif種類及排布完全一致，各組內(nèi)家族成員Motif種類差別較小。表明蓖麻bZIP家族的蛋白質(zhì)保守基序結(jié)構(gòu)較復(fù)雜，但與對應(yīng)的基因結(jié)構(gòu)基本一致（圖1）。

2.3 系統(tǒng)發(fā)育進化樹

由圖2可知，蓖麻轉(zhuǎn)錄因子保守結(jié)構(gòu)域和系統(tǒng)進化樹將進化樹分為5個亞族，蓖麻轉(zhuǎn)錄因子家族成員在各個亞族中呈不均勻分布。Ⅴ組中含有蓖麻bZIP基因最多為14個，Ⅰ、Ⅱ、Ⅲ和IV組中分別含有9、8、6和9個蓖麻bZIP基因。蓖麻和擬南芥的bZIP基因在多個分支中均有分布，表明二者可能具有相似的進化方式。

2.4 蛋白質(zhì)理化性質(zhì)

蛋白質(zhì)理化性質(zhì)方面，46個蓖麻bZIP蛋白的氨基酸數(shù)目和分子量有較大差異，其中有9個蛋白的氨基酸數(shù)小于200，分子量較?。籖cbZIP28只有136個氨基酸，分子量最小，為15 781.28 D；RcbZIP34則有1 223個氨基酸，分子量最高，為135 988.36 D。RcbZIP17等電點最低，為5.13；RcbZIP15等電點最高，為10.06。另外，46個蓖麻bZIP蛋白的疏水系數(shù)均為負值，即蓖麻bZIP蛋白均為親水性蛋白。在亞細胞定位中，除了RcbZIP16位于葉綠體，其余45個蓖麻bZIP蛋白全部分布于細胞核內(nèi)，而RcbZIP25在細胞核內(nèi)和細胞膜中均含有（表3）。

2.5 蛋白磷酸化位點預(yù)測

由表4可知，各基因家族成員的磷酸化位點均是絲氨酸（Ser）占比最大，且其磷酸位點數(shù)量明顯高于蘇氨酸（Thr）和酪氨酸（Tyr）。其中bZIP32具有94個磷酸化位點，而RcbZIP28磷酸化位點最少，僅11個。RcbZIP32有最多的絲氨酸（Ser）位點，為75個，RcbZIP16有最多的蘇氨酸（Thr）位點，為27個。說明RcbZIP家族可能主要通過絲氨酸的磷酸化來發(fā)揮生物學(xué)作用。

2.6 蛋白互作網(wǎng)絡(luò)分析

由圖3可知，在bZIP蛋白基因家族中，各家族成員相互作用，其中以RcbZIP18、RcbZIP27、RcbZIP2、RcbZIP17和RcbZIP46之間的聯(lián)系較多，其他蛋白之間也相互聯(lián)系。

2.7 蓖麻bZIP基因家族表達模式

在GEO數(shù)據(jù)庫中下載蓖麻各個組織的bZIP基因表達量數(shù)據(jù)，其結(jié)果見圖4。蓖麻bZIP基因在不同組織中的表達量差異較大，具有組織特異性。其中多數(shù)的蓖麻bZIP基因僅在特定的組織部位中高表達，例如RcbZIP13、RcbZIP31、RcbZIP16、RcbZIP3和RcbZIP28在蓖麻各個生長時期的表達量較高，而RcbZIP20、RcbZIP24、RcbZIP15、RcbZIP11、RcbZIP5、RcbZIP38和RcbZIP39則在蓖麻各個生長時期的表達量均相對較低；其中，RcbZIP3和RcbZIP28，以及RcbZIP38和RcbZIP39在蓖麻各個生長時期表達量基本一致。

2.8 bZIP基因的熒光定量表達

選擇具有代表性的差異基因，在蓖麻種子和葉片中測定表達量，結(jié)果如圖5所示。4個基因的熒光定量表達量與GEO數(shù)據(jù)庫表達量大致相同。RcbZIP13在葉片和種子中均為高表達，而RcbZIP10、RcbZIP25和RcbZIP44在葉片和種子中具有組織差異性表達。RcbZIP44基因在種子中的表達量顯著高于在葉片中的表達量。

3" 結(jié)論與討論

植物中bZIP轉(zhuǎn)錄因子家族成員較多，是分布較廣且高度保守的轉(zhuǎn)錄因子，參與植物的生長調(diào)控過程。安爽等[32]研究發(fā)現(xiàn)，不同基因bZIP在藍莓不同組織或是不同pH條件下表達模式有差異，具有抗逆調(diào)控作用；靳正偉等[33]通過對蓖麻bZIP基因家族成員在不同組織和種子發(fā)育階段的表達規(guī)律的研究，鑒別出RcbZIP8、RcbZIP31和RcbZIP40基因的表達可能與種子油脂累積具有較強的相關(guān)性。本研究在蓖麻全基因組水平上共鑒定了46個bZIP轉(zhuǎn)錄因子家族成員。種內(nèi)系統(tǒng)發(fā)育樹顯示，根據(jù)基因結(jié)構(gòu)和蛋白質(zhì)保守基序，46個RcbZIP基因被分為11個亞族，這11個亞族的蛋白質(zhì)保守基序結(jié)構(gòu)較復(fù)雜，但與對應(yīng)的基因結(jié)構(gòu)基本一致，其中外顯子與內(nèi)含子結(jié)構(gòu)與分布，可能攜帶著基因家族的進化痕跡。蓖麻和擬南芥種間系統(tǒng)發(fā)育樹表明，RcbZIP基因家族成員與擬南芥家族成員在進化上相似，均較為保守。通過基因的組織表達分析可以直觀了解其功能，本研究發(fā)現(xiàn)，基因在組織中的表達量具有明顯的組織特異性。多數(shù)蓖麻bZIP基因家族成員在蓖麻發(fā)育過程均存在較高的表達，推測該基因家族在蓖麻發(fā)育過程中發(fā)揮了重要的作用。蓖麻是重要的油料作物，胚乳是其油脂積累的重要部位，部分bZIP基因在胚組織中的表達量低，推測這些基因與油脂積累無關(guān)，或是油脂積累的負調(diào)控因子。部分基因在胚組織中的表達量較高，推測這些基因可能與蓖麻油脂積累與代謝通路密切相關(guān)，可能是甘油三酯合成代謝通路中的關(guān)鍵轉(zhuǎn)錄調(diào)控因子或關(guān)鍵酶，在油脂形成與累積過程中發(fā)揮關(guān)鍵作用。蛋白磷酸化分析結(jié)果顯示，絲氨酸是該基因家族的主要磷酸化位點，預(yù)測絲氨酸的磷酸化修飾對蓖麻bZIP基因家族的生物功能起著至關(guān)重要的作用。同時，該基因家族成員之間相互作用，進而參與蓖麻的各個生長時期。

綜上，本研究在全基因組水平上共篩選與鑒定出了46個蓖麻bZIP基因家族成員，采用生物信息學(xué)手段對其進行分析，發(fā)現(xiàn)46個蓖麻bZIP基因家族成員，可分為11個亞族，進化保守，位于細胞核、細胞膜和葉綠體中，主要通過絲氨酸的磷酸化來發(fā)揮生物學(xué)作用；參與葉、莖、根和胚等生長時期的發(fā)育調(diào)控，其中bZIP44基因在種子中的表達量明顯高于在葉片中的表達量，推測在種子油脂形成中發(fā)揮作用。研究結(jié)果為bZIP基因家族功能研究提供理論基礎(chǔ)，并為改良種質(zhì)及優(yōu)良品種選育提供參考。

參考文獻

[1] 宿明星，孫穎顥，施鶴，等. 植物非生物脅迫相關(guān)轉(zhuǎn)錄因子研究方法[J]. 生物技術(shù)通報，2015，31（1）：51-60.

[2] 岳曉琦，楊成坤，張月，等. 荔枝RAV家族轉(zhuǎn)錄因子的全基因組鑒定及表達分析[J]. 植物生理學(xué)報，2022，58（6）：1092-1104.

[3] 代夢媛，高梅，李文昌. 蓖麻GRF轉(zhuǎn)錄因子家族生物信息學(xué)鑒定及表達分析[J]. 分子植物育種，2021，19（22）：7383-7390.

[4] 代春艷，張笑晗，王曉麗，等. 芥菜bZIP基因家族的鑒定和表達分析[J]. 分子植物育種，2024，22（7）：2104-2118.

[5] 孫耀國，蔡天潤，姬行舟，等. 西洋梨全基因組bZIP基因家族生物信息學(xué)分析[J]. 林業(yè)與生態(tài)科學(xué)，2021，36（1）：24-34.

[6] 曹麗茹，馬晨晨，龐蕓蕓，等. 玉米bZIP轉(zhuǎn)錄因子ZmbZIP26的分子特征、亞細胞定位及響應(yīng)非生物脅迫的表達分析[J]. 華北農(nóng)學(xué)報，2023，38（4）：1-10.

[7] 潘婷婷. OsbZIP65的克隆與調(diào)控水稻開花的機制解析[D]. 大慶：黑龍江八一農(nóng)墾大學(xué)，2023.

[8] 游冉. 水稻轉(zhuǎn)錄因子bZIP34調(diào)控胚乳淀粉合成和抽穗期的機制研究[D]. 揚州：揚州大學(xué)，2023.

[9] LI H Y，LI L X，SHANGGUAN G D，et al. Genome-wide identification and expression analysis of bZIP gene family in Carthamus tinctorius L.[J]. Scientific reports，2020，10：15521.

[10] 叢嬌嬌，王曉宇，李平，等. 蓖麻耐鹽基因HKT的克隆及表達載體構(gòu)建[J]. 分子植物育種，2018，16（14）：4648-4657.

[11] 李威，李國瑞，黃鳳蘭，等. 蓖麻RcDELLA基因的克隆及序列分析[J]. 華北農(nóng)學(xué)報，2018，33（1）：14-19.

[12] 楊麗鳳，常如慧，李國瑞，等. 噴施草銨膦與草甘膦對蓖麻葉片生理生化的影響研究[J]. 內(nèi)蒙古民族大學(xué)學(xué)報（自然科學(xué)版），2018，33（5）：411-416.

[13] 鄭鷺. 蓖麻種質(zhì)資源鑒定評價及遺傳多樣性ISSR與SRAP分析[D]. 福州：福建農(nóng)林大學(xué)，2007.

[14] 張偉，龔文照，趙廣，等. 蓖麻油及其下游產(chǎn)業(yè)市場分析與發(fā)展探討[J]. 山西化工，2016，36（3）：24-27.

[15] OGUNNIYI D S. Castor oil：a vital industrial raw material[J]. Bioresource technology，2006，97（9）：1086-1091.

[16] 楊晗. 基于隱馬爾科夫模型的噬菌體尾部模塊識別[D]. 哈爾濱：哈爾濱工業(yè)大學(xué)，2022.

[17] 陳渝，鄧潔，陳君愉，等. 番茄FAR1/FHY3轉(zhuǎn)錄因子家族的全基因組鑒定及表達分析[J]. 植物生理學(xué)報，2021，57（10）：1983-1995.

[18] MOUNT D W. Using the basic local alignment search tool （BLAST）[J]. Cold spring harbor protocols，2007，2007（7）：pdb.top17.

[19] FINN R D，CLEMENTS J，EDDY S R. HMMER web server：interactive sequence similarity searching[J]. Nucleic acids research，2011，39（2）：29-37.

[20] LU S N，WANG J Y，CHITSAZ F，et al. CDD/SPARCLE：the conserved domain database in 2020[J]. Nucleic acids research，2020，48（1）：265-268.

[21] ZDOBNOV E M，APWEILER R. InterProScan-an integration platform for the signature-recognition methods in InterPro[J]. Bioinformatics，2001，17（9）：847-848.

[22] SIEVERS F，HIGGINS D G. Clustal omega，accurate alignment of very large numbers of sequences[M]//Methods in molecular biology. Totowa，NJ：Humana Press，2013.

[23] WATERHOUSE A M，PROCTER J B，MARTIN D M A，et al. Jalview Version 2-a multiple sequence alignment editor and analysis workbench[J]. Bioinformatics，2009，25（9）：1189-1191.

[24] BAILEY T L，ELKAN C. Fitting a mixture model by expectation maximization to discover motifs in biopolymers[J]. Proceedings. international conference on intelligent systems for molecular biology，1994，2：28-36.

[25] SUBRAMANIAN B，GAO S H，LERCHER M J，et al. Evolview v3：a webserver for visualization，annotation，and management of phylogenetic trees[J]. Nucleic acids research，2019，47（1）：270-275.

[26] CHEN J W. Dynamic research of growth，N，P，K and active ingredients in Angelica dahurica var. Formosana from Sichuan province[J]. China journal of Chinese materia medica，2010，35（21）：2812-2817.

[27] GUO C L，GUO R R，XU X Z，et al. Evolution and expression analysis of the grape （Vitis vinifera L.） WRKY gene family[J]. Journal of experimental botany，2014，65（6）：1513-1528.

[28] GASTEIGER E，HOOGLAND C，GATTIKER A，et al. Protein identification and analysis tools on the ExPASy server[M]//The proteomics protocols handbook. Totowa，NJ：Humana Press，2005.

[29] 練從龍，賴鐘雄，盧秉國，等. 荔枝古樹胚性愈傷組織Fe-SOD成員基因克隆及生物信息學(xué)分析[J]. 熱帶作物學(xué)報，2014，35（1）：74-81.

[30] SZKLARCZYK D，F(xiàn)RANCESCHINI A，KUHN M，et al. The STRING database in 2011：functional interaction networks of proteins，globally integrated and scored[J]. Nucleic acids research，2011，39（1）：561-568.

[31] 段強，何智彪，李國瑞，等. 蓖麻LBD基因家族全基因組鑒定、進化和表達分析[J]. 北方農(nóng)業(yè)學(xué)報，2021，49（6）：1-14.

[32] 安爽，劉金蓮，占曉慧，等. 藍莓bZIP轉(zhuǎn)錄因子基因預(yù)測及其在根中表達分析[J]. 浙江農(nóng)業(yè)科學(xué)，2021，62（10）：1977-1983.

[33] 靳正偉. 蓖麻bZIP家族鑒定及與油脂基因互作分析[D]. 合肥：中國科學(xué)技術(shù)大學(xué)，2014.

（責(zé)任編輯：胡立萍）