陳微　潘美紅　惠林沖　李威亞　郇國磊　何林玉　繆美華　陳振泰　楊海峰

摘要：bHLH轉(zhuǎn)錄因子家族是真核生物中重要的轉(zhuǎn)錄因子家族，在植物生長發(fā)育、次生代謝和逆境應(yīng)答中發(fā)揮重要作用。本研究利用轉(zhuǎn)錄組數(shù)據(jù)對洋蔥基因家族成員進(jìn)行篩選和鑒定，并對家族成員進(jìn)行理化性質(zhì)、保守基序、亞細(xì)胞定位、系統(tǒng)發(fā)育和蛋白互作生物信息學(xué)分析，同時對bHLH基因家族在洋蔥幼苗中的表達(dá)情況進(jìn)行分析。結(jié)果表明，從洋蔥轉(zhuǎn)錄組共鑒定到36個bHLH基因，預(yù)測編碼蛋白質(zhì)含有的氨基酸數(shù)量為127～610個，且均為親水性蛋白，系統(tǒng)發(fā)育分析表明，洋蔥bHLH基因家族分為12個亞族，同一亞族的大多數(shù)基因具有相似的保守基序，大多數(shù)bHLH 基因在洋蔥幼苗中表達(dá)，不同基因的表達(dá)量水平差異較大。本研究結(jié)果為進(jìn)一步研究洋蔥bHLH 轉(zhuǎn)錄因子家族的生物學(xué)功能奠定了基礎(chǔ)。

關(guān)鍵詞：洋蔥；bHLH；轉(zhuǎn)錄因子；生物信息學(xué)分析

中圖分類號：S633.201文獻(xiàn)標(biāo)志碼：A

文章編號：1002-1302（2023）05-0011-08

轉(zhuǎn)錄因子是真核生物中重要的調(diào)節(jié)因子，可以與順式元件相互作用，調(diào)節(jié)特定基因的表達(dá)，來響應(yīng)環(huán)境脅迫。bHLH（ basic helix-loop-helix）轉(zhuǎn)錄因子家族是植物中最大的轉(zhuǎn)錄因子家族之一［1］，具有高度保守的 bHLH 結(jié)構(gòu)域，由大約50～60個保守氨基酸組成，包含2個功能區(qū)：一個是N端的堿性氨基酸區(qū)，長度約為15個氨基酸，另一個是 HLH 區(qū)［2］。bHLH 蛋白調(diào)節(jié)許多生理過程，在黃酮類生物合成［3-4］、信號轉(zhuǎn)導(dǎo)［5］和植物生長發(fā)育［6-8］等方面發(fā)揮著重要作用。

隨著不同物種全基因組和轉(zhuǎn)錄組測序的完成，多個物種中鑒定到bHLH 基因家族。目前，已從擬南芥［9］、毛竹［10］、胡蘿卜［11］、人參［12］、梨［13］、蘋果［14］、葡萄［15］等植物物種鑒定和分析了大量bHLH基因家族。隨著 bHLH 蛋白功能的確定，bHLH 的分類不斷得到改進(jìn)。根據(jù)序列相似性和進(jìn)化關(guān)系以及結(jié)合 DNA的能力，bHLH通常分為從A到F的6個主要組［16-17］。

洋蔥是全球范圍內(nèi)種植的蔬菜之一，具有豐富的營養(yǎng)成分，很多國家視其為一種非常重要的蔬菜。迄今為止，bHLH轉(zhuǎn)錄因子的研究主要在擬南芥、葡萄、水稻等植物中，在洋蔥中未見bHLH轉(zhuǎn)錄因子功能及理化性質(zhì)的報道。本研究利用洋蔥轉(zhuǎn)錄組測序的數(shù)據(jù)庫，在轉(zhuǎn)錄水平上對bHLH基因家族進(jìn)行鑒定，并對其基本理化性質(zhì)、保守基序、蛋白互作、系統(tǒng)進(jìn)化以及基因表達(dá)情況等進(jìn)行生物信息學(xué)分析，此外，還篩選到一些可能參與洋蔥黃酮生物合成途徑的bHLH基因，結(jié)果可為洋蔥AcbHLH基因功能的研究提供理論參考。

1 材料與方法

1.1 植物材料

本試驗所用的材料來自江蘇省連云港市農(nóng)業(yè)科學(xué)院東辛基地，在苗期，選取長勢相近的3株紫皮洋蔥和3株黃皮洋蔥，選取幼嫩鱗莖的表皮，液氮速凍后，-80 ℃冰箱保存?zhèn)溆茫糜谵D(zhuǎn)錄組測序（轉(zhuǎn)錄組數(shù)據(jù)已發(fā)表）。

1.2 數(shù)據(jù)來源及bHLH基因的鑒定

將測序得到的原始序列進(jìn)行過濾，去除接頭污染、低質(zhì)量、N比例大于5%的Reads，得到質(zhì)量較高的Clean Reads，利用Trinity軟件對高質(zhì)量序列進(jìn)行混合組裝，共獲到77 474條Unigene［18］。分別從擬南芥數(shù)據(jù)庫 （ https：//www.arabidopsis.org）和水稻數(shù)據(jù)庫（http：//rapdb.dna.affrc.go.jp/download/irgsp1.html）下載擬南芥bHLH家族蛋白序列和水稻bHLH家族蛋白序列。然后使用HMMER 3.0軟件對獲取的已知bHLH蛋白家族序列構(gòu)建隱馬爾可夫模型，以此模型搜尋洋蔥轉(zhuǎn)錄組所有編碼蛋白序列，找出洋蔥轉(zhuǎn)錄組蛋白序列中潛在的bHLH家族序列。使用BLASTP把洋蔥轉(zhuǎn)錄組所有蛋白的序列比對獲取的bHLH家族參考序列，e-value設(shè)為1e-10，比對上的序列用作候選的bHLH家族序列。

1.3 bHLH 生物信息學(xué)分析

利用在線軟件ExPASy-ProtParam（http：//web.expasy.org/protparam/）對洋蔥bHLH蛋白序列的氨基酸數(shù)目、分子量、等電點、不穩(wěn)定系數(shù)等理化性質(zhì)進(jìn)行預(yù)測。利用MEME軟件（http：//meme.nbcr.net/meme）分析洋蔥bHLH家族的保守基序，保守基序的數(shù)目選擇15，其他參數(shù)均采用默認(rèn)值。利用ProtComp（http：//www.softberry.com/berry.phtml？topic=protcomppl&group=programs &subgroup=proloc）網(wǎng)站對洋蔥 bHLH家族成員的亞細(xì)胞定位進(jìn)行預(yù)測。

1.4 洋蔥bHLH 基因的系統(tǒng)發(fā)育

構(gòu)建方法：用MAFFT進(jìn)行多序列比對，然后用MEGA6.0軟件對鑒定得到的洋蔥、擬南芥的bHLH蛋白家族序列構(gòu)建系統(tǒng)進(jìn)化樹。使用鄰接法（neighbor-joining）和P-distance 模型構(gòu)建進(jìn)化樹，校驗參數(shù)Bootstrap設(shè)置為重復(fù)1 000次。

1.5 洋蔥bHLH 蛋白功能聯(lián)系網(wǎng)絡(luò)

利用 String 蛋白相互作用的數(shù)據(jù)庫（http：//string-db.org/）分析洋蔥bHLH 蛋白的互作網(wǎng)絡(luò)，以模式植物擬南芥作為參考物種。

1.6 洋蔥bHLH 基因的表達(dá)分析

利用洋蔥轉(zhuǎn)錄組測序數(shù)據(jù)分析bHLH基因的表達(dá)特征。以FPKM 值表示轉(zhuǎn)錄本豐度，通過對數(shù)據(jù)均一化處理，利用Heml 工具做熱圖。

2 結(jié)果與分析

2.1 洋蔥bHLH 基因家族的鑒定及理化性質(zhì)分析

為了從洋蔥轉(zhuǎn)錄組中獲得bHLH家族基因，使用擬南芥和水稻的bHLH蛋白序列進(jìn)行隱馬爾可夫模型和BLASTP搜索，對獲得的候選序列利用軟件pfamscan和Pfam A數(shù)據(jù)庫對目標(biāo)序列進(jìn)行結(jié)構(gòu)域注釋，確定將只含有保守結(jié)構(gòu)域bHLH 基序（PF00010）的序列作為最終的bHLH序列，最終獲得36個bHLH基因，根據(jù)轉(zhuǎn)錄組編號順序?qū)⑵涿麨锳cbHLH1～AcbHLH36（表1）。

對洋蔥36 個bHLH基因編碼蛋白質(zhì)的基本理化性質(zhì)進(jìn)行分析（表1），結(jié)果顯示，bHLH 基因編碼的氨基酸數(shù)目為127～610個，平均長度為297個氨基酸，預(yù)測其相對分子量大小為14 718.73～66 957.14 u；等電點在4.84～9.44范圍內(nèi)。從整體看，有67%的洋蔥bHLH家族蛋白質(zhì)等電點小于7，等電點處于酸性范圍內(nèi)，蛋白質(zhì)分子含有豐富的酸性氨基酸。蛋白質(zhì)疏水性分析結(jié)果表明，洋蔥中36個bHLH家族蛋白總平均親水性均小于0，說明該家族蛋白質(zhì)均屬于親水性蛋白質(zhì)，但是不同蛋白質(zhì)之間親水性存在差異。不穩(wěn)定指數(shù)在38.84～80.53之間，其中94.44%是不穩(wěn)定蛋白（Ⅱ>40），只有2個屬于穩(wěn)定蛋白（Ⅱ<40），脂肪系數(shù)介于55.70～94.32之間。亞細(xì)胞定位預(yù)測結(jié)果顯示，AcbHLH1等27個蛋白定位于細(xì)胞核，其余的9個定位于細(xì)胞外。

2.2 洋蔥bHLH蛋白保守基序分析

洋蔥bHLH蛋白中的保守基序用在線軟件MEME分析（圖1），結(jié)果表明，36個bHLH基因家族成員中保守基序的數(shù)目及種類有差異，除AcbHLH5外， 所有的bHLH都含有基序motif 1和motif 2， 而且motif1和motif2 總是相鄰，二者共同構(gòu)成了bHLH 結(jié)構(gòu)域。9個成員含有motif 3，6個成員基序數(shù)目和構(gòu)成一樣，都是motif 1、motif 2 和motif 4，AcbHLH8、AcbHLH12、AcbHLH33和AcbHLH34含有7個基序，是基序最多的bHLH蛋白。

2.3 洋蔥bHLH 蛋白家族進(jìn)化關(guān)系分析

對擬南芥的bHLH蛋白序列和篩選的36條洋蔥bHLH蛋白序列構(gòu)建系統(tǒng)進(jìn)化樹（圖2）。根據(jù)擬南芥中bHLH 轉(zhuǎn)錄因子的分類系統(tǒng)和洋蔥bHLH轉(zhuǎn)錄因子的進(jìn)化樹分類情況，36個bHLH可以劃分為12個亞族，其中數(shù)量最多的是H亞族，共13個，J亞族含有的bHLH數(shù)量次之，有9個，A、C、E和G亞族僅含有1個bHLH。

2.4 洋蔥bHLH蛋白互作關(guān)系

利用在線軟件STRING，預(yù)測洋蔥bHLH家族成員在擬南芥中同源蛋白的互作關(guān)系（圖3、圖4）。

結(jié)果顯示，洋蔥bHLH家族成員和多個蛋白相互作用。其中光敏色素互作蛋白PIL5和PIF3共同作用來調(diào)節(jié)光形態(tài)發(fā)生。

2.5 洋蔥bHLH基因家族表達(dá)分析

利用洋蔥（紫皮和黃皮）的轉(zhuǎn)錄組數(shù)據(jù)庫，找到候選的 36個bHLH基因?qū)?yīng)轉(zhuǎn)錄本的RPKM 值，通過對RPKM 值取對數(shù)值轉(zhuǎn)換，利用Heml 熱圖軟件作熱圖（圖5）。結(jié)果顯示，36個bHLH基因在紫皮和黃皮洋蔥中展現(xiàn)出不同的表達(dá)水平，AcbHLH7、AcbHLH32、AcbHLH33、AcbHLH22、AcbHLH6、AcbHLH35、AcbHLH5、AcbHLH19和AcbHLH26在紫皮洋蔥中的表達(dá)量較高。

3 討論與結(jié)論

作為植物中最大的轉(zhuǎn)錄因子家族之一，bHLH基因家族在植物的生長發(fā)育、新陳代謝以及應(yīng)對脅迫方面起著重要的作用?；谌蚪M研究的深入，越來越多物種的bHLH基因家族被陸續(xù)鑒定出來，其中擬南芥中鑒定出162個bHLH基因［19］，水稻中167個，番茄中152個［20］，葡萄中發(fā)現(xiàn)94個bHLH基因［15］，黃花蒿中鑒定出99個bHLH基因［21］，桃中有95個bHLH基因［22］，陳紅霖等從綠豆轉(zhuǎn)錄組中鑒定到122個bHLH轉(zhuǎn)錄因子［23］，Zhu等在龍血樹中發(fā)現(xiàn)了62個bHLH基因［24］。因為缺乏基因組，目前尚未見洋蔥bHLH 基因家族鑒定及其功能分析的報道，本研究利用洋蔥幼苗轉(zhuǎn)錄組數(shù)據(jù)庫，對洋蔥中bHLH基因家族成員進(jìn)行鑒定，并對AcbHLH基因家族中各蛋白的基本理化性質(zhì)、保守基序、亞細(xì)胞定位、系統(tǒng)進(jìn)化、基因表達(dá)等方面進(jìn)行了分析，為bHLH基因家族成員在結(jié)構(gòu)和功能以及基因表達(dá)調(diào)控等方面的研究提供參考。

通過分析鑒定，去除冗余，在洋蔥中鑒定出36個bHLH基因家族成員，根據(jù)轉(zhuǎn)錄組編號順序，將洋蔥36個bHLH基因分別命名為AcbHLH1～AcbHLH36。通過轉(zhuǎn)錄組在洋蔥中鑒定到的bHLH蛋白數(shù)量遠(yuǎn)遠(yuǎn)低于擬南芥、水稻和番茄中bHLH蛋白數(shù)量，這種情況可能與洋蔥基因組數(shù)據(jù)未獲得和轉(zhuǎn)錄組取樣時期有關(guān)，又或者測序組裝的條件限制有尚未發(fā)現(xiàn)的洋蔥bHLH基因家族成員。

洋蔥中36個bHLH基因的保守基序利用線軟件 MEME進(jìn)行分析，共找到15個高度保守基序，motif 1和motif 2 幾乎存在于所有的AcbHLH 蛋白中，它們總是彼此相鄰，共同構(gòu)成bHLH 結(jié)構(gòu)域，具有相似motif組成的洋蔥bHLH基因家族成員在聚類分析中距離較近，推測這些bHLH 保守基序在同一亞家族具有保守性及功能獨特性，而其他的保守基序也在每個亞家族發(fā)揮著不同的功能。

擬南芥和洋蔥bHLH蛋白的系統(tǒng)進(jìn)化結(jié)果表明，36個bHLH劃分為12個亞族，而擬南芥、水稻中的bHLH蛋白被劃分為21、22個亞族［24-25］，由于目前對洋蔥bHLH蛋白的功能了解較少，無法更深入的分析其原因。

目前，對洋蔥bHLH基因家族的研究遠(yuǎn)不如擬南芥、水稻等植物深入，bHLH 轉(zhuǎn)錄因子可以作為轉(zhuǎn)錄激活因子或抑制因子，在代謝和發(fā)育過程中起重要作用［26］。在擬南芥中已對大多數(shù)鑒定的 bHLH 蛋白進(jìn)行了功能鑒定。如轉(zhuǎn)錄因子PIF3 作為光感受器信號的直接光敏色素反應(yīng)伙伴，參與控制光調(diào)節(jié)基因的表達(dá)［9］。花青素是類黃酮中最大的一類［27］，是最重要的次級代謝產(chǎn)物之一。相關(guān)研究表明，MYB 和 WD40 轉(zhuǎn)錄因子可以與 bHLH 相互作用，參與花青素和黃酮醇的生物合成。bHLH2 和 bHLH46 在龍血樹中可能參與茉莉酸（JA）調(diào)節(jié)的黃酮類生物合成［28］。本研究中，通過網(wǎng)絡(luò)共表達(dá)發(fā)現(xiàn)bHLH 轉(zhuǎn)錄因子PIF3與HY5相互作用，HY5在花青素的合成中具有一定的作用，因此，推測PIF3可能在洋蔥花青素的合成中起到一定的作用，但具體的作用仍需進(jìn)一步驗證。

由于缺少基因組，目前在洋蔥分子生物學(xué)方面的研究相對滯后，未見洋蔥中bHLH 家族基因的報道，對該家族基因的功能和分類還缺乏系統(tǒng)深入的研究。本研究通過對轉(zhuǎn)錄組測序得到的洋蔥bHLH家族基因進(jìn)行鑒定，并對其進(jìn)行生物信息學(xué)、保守基序、系統(tǒng)進(jìn)化分析以及基因表達(dá)分析，為進(jìn)一步篩選洋蔥bHLH家族基因，揭示該家族基因參與植物生長發(fā)育的調(diào)控機(jī)制奠定理論基礎(chǔ)。

參考文獻(xiàn)：

［1］Hong Y Q，Ahmad N，Tian Y Y，et al. Genome-wide identification，expression analysis，and subcellular localization of Carthamus tinctorius bHLH transcription factors［J］. International Journal of Molecular Sciences，2019，20（12）：3044.

［2］Wang J Y，Hu Z Z，Zhao T M，et al. Genome-wide analysis of bHLH transcription factor and involvement in the infection by yellow leaf curl virus in tomato （Solanum lycopersicum）［J］. BMC Genomics，2015，16（1）：1-14.

［3］Xu W J，Dubos C，Lepiniec L. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes［J］. Trends in Plant Science，2015，20（3）：176-185.

［4］Zhao M，Li J，Zhu L，et al. Identification and characterization of MYB-bHLH-WD40 regulatory complex members controlling anthocyanidin biosynthesis in blueberry fruits development［J］. Genes，2019，10（7）：496.

［5］Zhu E G，You C J，Wang S S，et al. The DYT1-interacting proteins bHLH010，bHLH089 and bHLH091 are redundantly required for Arabidopsis anther development and transcriptome［J］. The Plant Journal，2015，83（6）：976-990.

［6］Lu R，Zhang J，Liu D，et al. Characterization of bHLH/HLH genes that are involved in brassinosteroid （BR） signaling in fiber development of cotton （Gossypium hirsutum）［J］. BMC Plant Biology，2018，18（1）：304.

［7］Xing B C，Yang D F，Yu H Z，et al. Overexpression of SmbHLH10 enhances tanshinones biosynthesis in Salvia miltiorrhiza hairy roots［J］. Plant Science，2018，276：229-238.

［8］Zhu Z G，Liang H L，Chen G P，et al. The bHLH transcription factor SlPRE2 regulates tomato fruit development and modulates plant response to gibberellin［J］. Plant Cell Reports，2019，38（9）：1053-1064.

［9］Toledo-Ortiz G，Huq E，Quail P H. The Arabidopsis basic/Helix-loop-Helix transcription factor family［J］. The Plant Cell，2003，15（8）：1749-1770.

［10］Cheng X R，Xiong R，Liu H L，et al. Basic helix-loop-helix gene family：genome wide identification，phylogeny，and expression in Moso bamboo［J］. Plant Physiology and Biochemistry，2018，132：104-119.

［11］Chen Y Y，Li M Y，Wu X J，et al. Genome-wide analysis of basic helix-loop-helix family transcription factors and their role in responses to abiotic stress in carrot［J］. Molecular Breeding，2015，35（5）：125.

［12］Chu Y，Xiao S M，Su H，et al. Genome-wide characterization and analysis of bHLH transcription factors in Panax ginseng［J］. Acta Pharmaceutica Sinica B，2018，8（4）：666-677.

［13］Jin C，Huang X S，Li K Q，et al. Overexpression of a bHLH1 transcription factor of Pyrus ussuriensis confers enhanced cold tolerance and increases expression of stress-responsive genes［J］. Frontiers in Plant Science，2016，7：441.

［14］Mao K，Dong Q L，Li C，et al. Genome wide identification and characterization of apple bHLH transcription factors and expression analysis in response to drought and salt stress［J］. Frontiers in Plant Science，2017，8（12）：480.

［15］Wang P F，Su L，Gao H H，et al. Genome-wide characterization of bHLH genes in grape and analysis of their potential relevance to abiotic stress tolerance and secondary metabolite biosynthesis［J］. Frontiers in Plant Science，2018，9：64.

［16］Atchley W R，F(xiàn)itch W M. A natural classification of the basic helix-loop-helix class of transcription factors［J］. Proceedings of the National Academy of Sciences of the United States of America，1997，94（10）：5172-5176.

［17］Simionato E，Ledent V，Richards G，et al. Origin and diversification of the basic helix-loop-helix gene family in metazoans：insights from comparative genomics［J］. BMC Evolutionary Biology，2007，7：33.

［18］陳 微，潘美紅，惠林沖，等. 洋蔥轉(zhuǎn)錄組測序及黃酮類化合物合成相關(guān)基因的分析［J/OL］.分子植物育種：1-6［2022-10-17］.http：//kns.cnki.net/kcms/detail/46.1068.S.20210816.1540.010.html.

［19］Carretero-Paulet L，Galstyan A，Roig-Villanova I，et al. Genome-wide classification and evolutionary analysis of the bHLH family of transcription factors in Arabidopsis，poplar，rice，moss，and algae［J］. Plant Physiology，2010，153（3）：1398-1412.

［20］Sun H，F(xiàn)an H J，Ling H Q. Genome-wide identification and characterization of the bHLH gene family in tomato［J］. BMC Genomics，2015，16（1）：9.

［21］甘 雨，吳 端，張 棟，等. 黃花蒿bHLH轉(zhuǎn)錄因子基因家族鑒定及光調(diào)控分析［J］. 中國現(xiàn)代中藥，2021，23（3）：441-452.

［22］Zhang C H，F(xiàn)eng R C，Ma R J，et al. Genome-wide analysis of basic helix-loop-helix superfamily members in peach［J］. PLoS One，2018，13（4）：e0195974.

［23］陳紅霖，胡亮亮，王麗俠，等. 綠豆bHLH轉(zhuǎn)錄因子家族的鑒定與生物信息學(xué)分析［J］. 植物遺傳資源學(xué)報，2017，18（6）：1159-1167.

［24］Zhu J H，Xia D N，Xu J，et al. Identification of the bHLH gene family in Dracaena cambodiana reveals candidate genes involved in flavonoid biosynthesis［J］. Industrial Crops and Products，2020，150：112407.

［25］Heim M A. The basic helix-loop-helix transcription factor family in plants：a genome-wide study of protein structure and functional diversity［J］. Molecular Biology and Evolution，2003，20（5）：735-747.

［26］Li X X，Duan X P，Jiang H X，et al. Genome-wide analysis of basic/helix-loop-helix transcription factor family in rice and Arabidopsis［J］. Plant Physiology，2006，141（4）：1167-1184.

［27］Feller A，Machemer K，Braun E L，et al. Evolutionary and comparative analysis of MYB and bHLH plant transcription factors［J］. The Plant Journal，2011，66（1）：94-116.

［28］Welch C R，Wu Q L，Simon J E. Recent advances in anthocyanin analysis and characterization［J］. Current Analytical Chemistry，2008，4（2）：75-101.

收稿日期：2022-04-21

基金項目：江蘇省農(nóng)業(yè)科技自主創(chuàng)新資金［編號：CX（20）3164］；連云港市財政專項（編號：QNJJ2005）。

作者簡介：陳 微（1990—），女，山東菏澤人，碩士，助理研究員，主要從事蔬菜育種和栽培研究。E-mail：1135050599@qq.com。

通信作者：楊海峰，副研究員，主要從事蔬菜育種和栽培研究。E-mail：13775492959@163.com。