鄭玲 代琪芳 陳雪冬

摘要：bHLH轉錄因子在植物的生長發(fā)育、生物合成及脅迫反應等方面起著重要作用。利用生物信息學方法對森林草莓全基因組序列進行bHLH家族的搜索和鑒定，共鑒定出112個森林草莓的bHLH基因，這些基因分別被定位在7條染色體上，但分布不均勻，在5號染色體上分布得最多，3號染色體上分布得最少。進化樹結果顯示，森林草莓的bHLH家族可以進一步被分為12個亞家族，其中第Ⅰ、Ⅲ亞族的成員較多。森林草莓bHLH家族基因的表達模式分析結果顯示，有73個FvbHLH基因與花的發(fā)育相關，有28個FvbHLH基因與果實的發(fā)育相關，其中FvbHLH97/36/55/39/42/52/26/44/68/54對森林草莓花的整個發(fā)育過程都尤為重要，而FvbHLH57/59/66/78/100/15/17/46則在果實發(fā)育中起到主要的調控作用。

關鍵詞：森林草莓;bHLH;進化樹;花發(fā)育;果實發(fā)育

中圖分類號： S668.401? 文獻標志碼： A

文章編號：1002-1302（2022）02-0047-06

收稿日期：2021-03-28

基金項目：河南省高等學校重點科研項目（編號：20A180018）;洛陽師范學院博士科研啟動項目（編號：180141051235）。

作者簡介：鄭 玲（1981—），女，河南南陽人，碩士，副教授，主要從事植物學研究。E-mail：linda2266@qq.com。

轉錄因子又稱為反式作用因子，是一類能夠結合在順式作用元件上調控基因表達的DNA結合蛋白[1]。目前，人們已經從植物中鑒定出很多轉錄因子，其中bHLH（堿性螺旋-環(huán)-螺旋）是第二大類轉錄因子[2]。bHLH轉錄因子因含有高度保守的bHLH結構域而被命名，bHLH結構域約含有60個保守氨基酸序列，由basic區(qū)域、HLH（helix-loop-helix）區(qū)域組成[3]。位于N端的basic區(qū)域是 DNA結合的堿性區(qū)域，約含有15個氨基酸，其中有6個是堿性氨基酸殘基;HLH區(qū)域位于C末端，主要由疏水氨基酸組成，并參與形成同源或異源二聚體[4]。

bHLH轉錄因子在真核生物中最為廣泛存在，其功能也具有多樣性。在酵母和其他單細胞真核生物中，bHLH參與染色體分離和代謝調控等過程[5]。bHLH轉錄因子可參與植物的物質合成和信號轉導過程[6]，并在植物生長發(fā)育[7]及抗旱[8]、耐鹽[9]等脅迫應答反應中起著關鍵作用。植物中最早被發(fā)現(xiàn)的bHLH轉錄因子（即玉米中由R基因編碼的Lc蛋白）的主要功能為參與花青素的生物合成[10]。隨后，擬南芥的AtbHLH1/2/42、水稻的OsbHLH013/016等基因都被證明與花青素的合成調控有關[11-13]，擬南芥AtbHLH40則與赤霉素的合成轉導有關[14]。在響應植物脅迫方面，擬南芥AtbHLH112是鹽、干旱及滲透脅迫的正調控因子[15];AtbHLH116基因是耐寒性的調節(jié)子，能夠提高植物的抗凍能力[16]。在植物生長發(fā)育方面，擬南芥的AtbHLH9、AtbHLH65基因均是調控植物生長發(fā)育的關鍵信號樞紐[17];AtbHLH95基因主要在胚乳中表達，控制胚胎發(fā)育，從而導致植物生長遲緩[18];擬南芥的AtbHLH21和水稻的OsbHLH164對花粉粒的形成及絨氈層細胞的發(fā)育有著關鍵的調控作用[19-20]。

森林草莓（Fragaria vesca）也稱野草莓，是薔薇科草莓屬的多年生草本植物，屬于二倍體植物[21]，具有相當高的食用價值和藥用價值。目前，國內外在分子生物學方面對擬南芥、水稻等植物bHLH家族的研究已經較為成熟，然而還沒有關于森林草莓bHLH轉錄因子鑒定和分析的報道。本研究利用生物信息學方法對森林草莓bHLH（FvbHLH）家族的成員進行鑒定，并對其理化性質、系統(tǒng)進化、基因在染色體上的分布、基因結構及表達模式等方面進行全面分析，以期為將來深入研究該家族基因的結構和功能提供重要參考。

1 材料與方法

1.1 森林草莓bHLH家族的鑒定

本研究對象是森林草莓（Fragaria vesca）。2020年3—6月在洛陽師范學院生命科學學院生物信息學實驗室（214）電腦上進行軟件操作和信息處理。森林草莓bHLH家族的蛋白質序列數據來自Phytozome（https：//phytozome.jgi.doe.gov/pz/portal.html？）與PlantTFDB（http：//planttfdb.cbi.pku.edu.cn/）數據庫，由于2個數據庫中搜索到的結果會存在差異，因此將2類搜索結果進行序列比對后，手動去除冗余序列，隨后利用InterPro（http：//www.ebi.ac.uk/interpro/）在線軟件鑒定搜索的候選基因是否具有bHLH結構域，進而獲得bHLH家族的所有成員。

1.2 森林草莓bHLH家族的理化性質分析

通過在線軟件ExPASy（https：//web.expasy.org/protparam/）得到森林草莓bHLH蛋白序列的基本理化性質信息，如理論等電點、氨基酸數量、分子量及總水平疏水性等。

1.3 森林草莓bHLH家族系統(tǒng)進化樹的構建

從PlantTFDB數據庫下載森林草莓和擬南芥的bHLH序列，隨后通過MEGA 6.0的Clustal W工具進行多序列比對，將獲得的比對結果用MEGA 6.0軟件構建系統(tǒng)進化樹[22]。同時采用相同方法單獨構建森林草莓的bHLH家族進化樹。

1.4 森林草莓bHLH基因在染色體上的定位

根據Phytozome數據庫中搜索到的FvbHLH基因在染色體上的位置信息，利用MapInspect軟件繪制基因在染色體上的分布圖[23]。

1.5 森林草莓bHLH基因在花、果實發(fā)育過程中的作用分析

從文獻[24-25]的附件中分別下載森林草莓花發(fā)育、果實發(fā)育過程中所有轉錄因子的表達量，進一步篩選出bHLH基因在各組織中的轉錄表達量。用MeV4軟件分別對檢測到的FvbHLH基因表達量進行聚類分析[26]。

2 結果與分析

2.1 森林草莓bHLH家族基因的鑒定及理化特性分析

通過PlantTFDB、Phytozome兩大數據庫共鑒定出112個候選bHLH家族基因，利用InterPro進行驗證，結果顯示，這112個基因編碼的蛋白質都具有完整的bHLH結構域。同時利用ExPASy在線分析工具得到這112個FvbHLH家族基因編碼蛋白的理論等電點、氨基酸數、原子組成及脂溶指數等理化信息。根據基因在染色體上的位置，將這112個轉錄因子依次命名為FvbHLH1～FvbHLH112。

結果顯示，各蛋白質分子中的氨基酸數量在90（FvbHLH64）～1 496（FvbHLH27）個之間，平均為434個;蛋白質分子量在10 296.6 u（FvbHLH64）～166 739.0 u（FvbHLH27）之間，平均為47 959.6 u;理論等電點在4.35（FvbHLH9）～10.94（FvbHLH106）之間，并且有70個蛋白質的等電點在4～7之間，說明森林草莓bHLH家族蛋白的氨基酸序列的等電點大多在酸性范圍內。有73個蛋白質分子的不穩(wěn)定指數均大于50，說明bHLH蛋白質的狀態(tài)不穩(wěn)定。對蛋白質的親疏水性進行分析發(fā)現(xiàn)，除FvbHLH76外，其余蛋白質的總水平疏水性均為負值，表明森林草莓bHLH家族蛋白大多屬于親水蛋白。

2.2 森林草莓bHLH家族基因進化樹的構建

選取111個擬南芥bHLH家族成員作為參考來構建森林草莓bHLH家族的系統(tǒng)進化樹（圖1）。根據擬南芥bHLH轉錄因子的分類系統(tǒng)，可將森林草莓的112個bHLH成員劃分為12個亞族。其中第Ⅵ亞族成員最少，只有1個FvbHLH、2個AtbHLH;其次是第Ⅱ亞族，有1個FvbHLH、3個AtbHLH;第Ⅲ亞族的成員數最多，包含19個FvbHLH、17個AtbHLH;第Ⅰ亞族的基因數也較多，有15個FvbHLH、17個AtbHLH。其余亞家族的成員數量不一。

2.3 森林草莓bHLH家族基因在染色體上的定位

利用MapInspect軟件繪制出基因在染色體上的分布情況，由圖2可以看出，在森林草莓的全部bHLH基因中，F(xiàn)vbHLH104～FvbHLH112尚未定位到染色體上，其余103個基因在7條染色體上均有分布，但分布不均勻。其中5號染色體上分布的基因最多，有22個（FvbHLH55～FvbHLH76）;其次是2、4號染色體，都各有17個基因（FvbHLH12～FvbHLH28和FvbHLH38～FvbHLH54）;3號染色體上分布的基因較少，為9個（FvbHLH29～FvbHLH37），其余染色體上的基因分布數量為11～15個不等。多條染色體上出現(xiàn)了由2個或多個基因聚集在一起形成的基因簇。

2.4 森林草莓bHLH基因在花發(fā)育過程中的作用分析

用MeV4軟件對檢測到的FvbHLH家族成員在花發(fā)育過程中相關基因的表達量進行聚類分析。由圖3可以看出，在森林草莓的花瓣（flower）、花被（perianth）、花托（receptacle）、花藥（anther）及小孢子（micrsopore）、花粉（pollen）和雌蕊（carpel）發(fā)育的過程中，共有73個FvbHLH家族基因表達。

結果顯示，有10個基因（FvbHLH97/36/55/39/42/52/26/44/68/54）在所有組織中的表達量均較高。FvbHLH90/19/9/107/80/67/22/1/11/18/15/51/24/46這14個基因除了在花粉中不表達之外，在其余組織中均大量表達。同樣的，F(xiàn)vbHLH57/93/37/43/78/71也僅在小孢子、花粉中的表達量較低。此外，有些基因的表達模式顯示出組織特異性，例如FvbHLH53在花藥、花粉和雌蕊中高表達;FvbHLH13在小孢子、花藥中高表達;FvbHLH86在雌蕊、花藥中高表達;FvbHLH89/34在花藥、花粉中高表達;FvbHLH16只在花粉中大量表達;FvbHLH12/4只在花藥發(fā)育過程中的表達水平較高;FvbHLH10/

112/88/63/29只在雌蕊中高表達。

2.5 森林草莓bHLH家族基因在果實發(fā)育過程中的作用分析

通過MeV4軟件對檢測到的FvbHLH家族成員在果實發(fā)育過程中相關基因的表達量進行聚類分析。由圖4可以看出，共有28個FvbHLH家族基因在森林草莓胚胎（embryo）、胚珠（ovule）、種皮（ghost）、種子（seed）及子房壁（wall）、果皮（cortex）、木髓（pith）的發(fā)育過程中表達。

研究發(fā)現(xiàn)，在胚珠、種子和子房壁發(fā)育的過程中，有14個基因即FvbHLH90/96/58/40/57/91/16/66/78/100/81/15/17/46均顯示出高表達，其中FvbHLH96/57/91/16/66/78/100/81/15/17/46在種皮整個發(fā)育階段中的表達量也很高，還有部分基因的表達水平則在種皮的不同發(fā)育階段存在較大差異，例如FvbHLH10/91/33的表達量持續(xù)增加，而FvbHLH58/66/15的表達量卻逐漸減少。在果皮、木髓中，F(xiàn)vbHLH77/57/59/112/99/66/78/100/79/64/30/15/17/46共14個基因的表達水平都很高，有些基因則只在這2種組織的特定發(fā)育階段表達，例如FvbHLH90、FvbHLH84分別只在果皮、木髓發(fā)育的 1～2階段顯示出較高的表達量，F(xiàn)vbHLH72只在果皮發(fā)育的4～5階段高表達。在胚胎發(fā)育過程中，bHLH基因的表達量較少，只有FvbHLH90/10/59/101/33/88/81/72共8個基因的表達量較高。

3 討論

目前， 還未見關于森林草莓bHLH轉錄因子的

報道，本研究利用生物信息學方法鑒定得到112個森林草莓bHLH轉錄因子，對其理化性質進行全面分析發(fā)現(xiàn)，bHLH轉錄因子大多是不穩(wěn)定的親水蛋白，且等電點大多處于酸性范圍。在擬南芥、水稻中的bHLH蛋白等電點也存在這種現(xiàn)象，表明bHLH轉錄因子的保守結構域中雖然存在basic堿性區(qū)域，但酸性氨基酸也同樣富集于蛋白質分子中[27]。

在模式植物擬南芥中，研究者將133個bHLH基因劃分為12個亞家族[28];水稻中的167個bHLH轉錄因子被分為22個亞家族[29];159個馬鈴薯[30]、230個大白菜[31]bHLH家族分別被劃分為21、24個亞族。本研究參照擬南芥的分類方式，將森林草莓的112個bHLH轉錄因子分成12個亞族，其中第Ⅱ、Ⅵ亞族的成員最少，均只有1個FvbHLH基因;第Ⅲ亞族的成員最多，為19個。在擬南芥bHLH亞家族中較少的也是第Ⅱ、Ⅵ亞族，第Ⅲ亞族的bHLH成員也較多。通過染色體定位分析發(fā)現(xiàn)，bHLH基因在各號染色體上的分布都不均勻，并且形成了基因簇，說明它們可能屬于串聯(lián)復制。

已有大量研究發(fā)現(xiàn)，bHLH轉錄因子與植物的脅迫反應有關，包括抗鹽、干旱、缺鐵等抗脅迫過程。進化樹分析發(fā)現(xiàn)，森林草莓FvbHLH30與擬南芥AtbHLH112[15]聚在同一分支，推測其在抗鹽、耐旱和抗?jié)B透脅迫方面也有重要作用;FvbHLH59與AtbHLH116[16]位于同一分支，因此FvbHLH59可能與植物的抗低溫調控有關。bHLH家族在植物信號轉導、物質合成方面也有關鍵作用。森林草莓FvbHLH93/57/21與擬南芥AtbHLH1/2/42[11-12]的進化關系較近，推測它們與花青素的合成有關;FvbHLH6與AtbHLH40[14]聚在同一分支，推測FvbHLH6也參與了赤霉素的合成轉導。

表達模式分析結果顯示，在花的發(fā)育過程中顯示出高表達的部分基因主要集中于第Ⅲ、Ⅳ、Ⅸ及Ⅺ、Ⅻ亞族，其中FvbHLH97/36/55/39/42/52/26/44/68/54的表達量在森林草莓花的整個發(fā)育過程中都很高，表明這些基因對于花的發(fā)育尤為重要。FvbHLH57/59/66/78/100/15/17/46幾乎在果實發(fā)育的全過程表達，說明這8個基因與森林草莓果實的發(fā)育調控有關。草莓是一種重要的小型漿果，其果實主要由花托發(fā)育而來[24]。由此推測，在花托發(fā)育中起較大調控作用的FvbHLH19/21/57/9等基因可能對改善森林草莓果實的品質起著重要的作用。

本研究對森林草莓bHLH家族成員的理化性質、進化關系、染色體分布及組織表達模式等方面進行了全面分析，初步分析了bHLH基因在森林草莓花和果實發(fā)育過程中的作用，可為提高其經濟價值提供理論基礎和研究方向，并為進一步了解森林草莓bHLH轉錄因子結構與功能之間的關系和在植物生長發(fā)育過程中的調控機制提供重要的依據。

參考文獻：

[1]Wu S，Gallagher K L. Transcription factors on the move[J]. Current Opinion in Plant Biology，2012，15（6）：645-651.

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

[3]Feng X M，Zhao Q，Zhao L L，et al. The cold-induced basichelix-loop-helix transcription factor gene MdCIbHLH1 encodes an ICE-like protein in apple[J]. BMC Plant Biology，2012，12（1）：1-14.

[4]Massari M E，Murre C. Helix-loop-helix proteins：regulators of transcription in eucaryotic organisms[J]. Molecular and Cellular Biology，2000，20（2）：429-440.

[5]Robinson K A，Lopes J M.SURVEY AND SUMMARY：Saccharomyces cerevisiae basic helix-loop-helix proteins regulate diverse biological processes[J]. Nucleic Acids Research，2000，28（7）：1499-1505.

[6]Fairchild C D，Schumaker M A，Quail P H.HFR1 encodes an atypical bHLH protein that acts in phytochrome A signal transduction[J]. Gene Dev，2000，14（18）：2377-2391.

[7]Penfield S，Josse E M，Kannangara R，et al. Cold and light control seed germination through the bHLH transcription factor SPATULA[J]. Curr Biol，2005，15（22）：1998-2006.

[8]Abe H，Yamaguchi-Shinozaki K，Urao T，et al. Role of Arabidopsis MYC and MYB homologs in drought-and abscisic acid-regulated gene expression[J]. Plant Cell，1997，9（10）：1859-1868.

[9]Jiang Y Q，Deyholos M K. Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes[J]. BMC Plant Biology，2006，6：25.

[10]Ludwig S R，Habera L F，Dellaporta S L，et al. Lc，a member of the maize R gene family responsible for tissue-specific anthocyanin production，encodes a protein similar to transcriptional activators and contains the myc-homology region[J]. Proceedings of the National Academy of Sciences of the United States of America，1989，86（18）：7092-7096.

[11]Nesi N，Debeaujon I，Jond C，et al. The TT8 gene encodes a basic helix-loop-helix domain protein required for expression of DFR and BAN genes in Arabidopsis siliques[J]. Plant Cell，2000，12（10）：1863-1878.

[12]Ramsay N A，Walker A R，Mooney M，et al. Two basic-helix-loop-helix genes （MYC-146 and GL3） from Arabidopsis can activate anthocyanin biosynthesis in a white-flowered Matthiola incana mutant[J]. Plant Molecular Biology，2003，52（3）：679-688.

[13]Sakamoto W，Ohmori T，Kageyama K，et al. The purple leaf （Pl） locus of rice：the Plw allele has a complex organization and includes two genes encoding basic helix-loop-helix proteins involved in anthocyanin biosynthesis[J]. Plant and Cell Physiology，2001，42（9）：982-991.

[14]Arnaud N，Girin T，Sorefan K，et al. Gibberellins control fruit patterning in Arabidopsis thaliana[J]. Genes & Development，2010，24（19）：2127-2132.

[15]Liu W，Tai H，Li S，et al. bHLH122 is important for drought and osmotic stress resistance in? Arabidopsis and in the repression of ABA catabolism[J]. New Phytologist，2014，201（4）：1192-1204.

[16]Chinnusamy V，Ohta M，Kanrar S，et al. ICE1：a regulator of cold-induced transcriptome and freezing tolerance in Arabidopsis[J]. Genes Development，2003，17（8）：1043-1054.

[17]Paik I，Kathare P K，Kim J I，et al. Expanding roles of PIFs in signal integration from multiple processes[J]. Molecular Plant，2017，10（8）：1035-1046.

[18]Kondou Y，Nakazawa M，Kawashima M，et al. RETARDED GROWTH OF EMBRYO1，a new basic helix-loop-helix protein，expresses in endosperm to control embryo growth[J]. Plant Physiology，2008，147（4）：1924-1935.

[19]Sorensen A M，Krber S，Unte U S，et al. The Arabidopsis ABORTED MICROSPORES （AMS） gene encodes a MYC class transcription factor[J]. Plant Journal for Cell & Molecular Biology，2003，33（2）：413-423.

[20]Jung K H，Han M J，Lee Y S，et al. Rice undeveloped Tapetum1 is a major regulator of early tapetum development[J]. The Plant Cell，2005，17（10）：2705-2722.

[21]周鶴瑩，張 瑋，張 卿，等. 森林草莓‘Hawaii4’高效遺傳轉化系統(tǒng)的建立[J]. 北京農學院學報，2015，30（1）：10-14.

[22]Tamura K，Stecher G，Peterson D，et al. MEGA6：molecular evolutionary genetics analysis version 6.0[J]. Molecular Biology and Evolution，2013，30（12）：2725-2729.

[23]Zhang S P，Liu M M，Miao H，et al. Chromosomal mapping and QTL analysis of resistance to downy mildew in Cucumis sativus[J]. Plant Disease，2013，97 （2）：245-251.

[24]Hollender C A，Kang C，Darwish O，et al. Floral transcriptomes in woodland strawberry uncover developing receptacle and anther gene networks[J]. Plant Physiol，2014，165（3）：1062-1075.

[25]Kang C，Darwish O，Geretz A，et al. Genome-scale transcriptomic insights into early-stage fruit development in woodland strawberry Fragaria vesca[J]. Plant Cell，2013，25（6）：1960–1978.

[26]Howe E，Holton K，Nair S，et al. MeV：multi experiment viewer[M]. Boston：Springer，2010：267–277.

[27]應炎標，朱友銀，郭衛(wèi)東，等. 櫻桃bHLH轉錄因子家族基因鑒定及表達分析[J]. 分子植物育種，2018，16（14）：4559-4568.

[28]Heim M A，Jakoby M，Werber M，et al. 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.

[29]Li X，Duan X，Jiang H，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.

[30]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.

[31]Song X，Huang Z，Duan W，et al. Genome-wide analysis of the bHLH transcription factor family in Chinese cabbage（Brassica rapa ssp. pekinensis）[J]. Molecular Genetics and Genomics，2014，289（1）：77-91.