林顯祖，肖小虎，陽江華，秦云霞，龍翔宇，方永軍

摘? 要：鈣調(diào)蛋白結合轉(zhuǎn)錄激活因子（Calmodulin-binding transcription activator， CAMTA）是一類存在于真核生物中結構保守的重要鈣調(diào)素結合蛋白，廣泛參與植物的脅迫應答、生長發(fā)育及激素信號轉(zhuǎn)導等生物學過程。本文以巴西橡膠樹（Hevea brasiliensis）為研究對象，從全基因組水平鑒定到8個橡膠樹HbCAMTA基因，并對這些基因的理化性質(zhì)、基因結構、表達模式、啟動子區(qū)調(diào)控元件進行分析。結果顯示，該基因家族所有成員均定位于細胞核內(nèi)，編碼蛋白均為親水性蛋白，分子量在104.34～121.99 kD之間，等電點在5.32～7.83之間;家族基因根據(jù)蛋白序列同源性可以分成4個亞組，亞組成員在基因結構、保守基序數(shù)目和相對位置等方面具有較高相似性。調(diào)控HbCAMTA的miRNA預測分析顯示，每個家族成員受到至少3個miRNA調(diào)控。表達分析結果顯示，HbCAMTA在不同組織的表達有明顯的差異，在葉片中隨著葉片的發(fā)育表達不斷上調(diào)，在成熟期或變色期達到最高的表達水平;除HbCAMTA2外，乙烯刺激未能顯著誘導橡膠樹其他CAMTA基因的表達，而低溫可以顯著誘導HbCAMTA的上調(diào)表達。HbCAMTA1是膠乳、葉片、根、樹皮、雄花、雌花以及葉片整個發(fā)育過程（古銅期、變色期、淡綠期、成熟期）中最高表達的基因家族成員，表明其具有非組織特異性的基礎生物學功能。HbCAMTA3在種子中表達水平最高且受低溫脅迫誘導顯著上調(diào)表達，表明其在響應低溫脅迫過程中具有重要的作用。HbCAMTA7在大多數(shù)組織、乙烯處理、低溫脅迫處理以及葉片不同發(fā)育時期均是最低表達。本研究初步揭示了橡膠樹CAMTA家族成員的理化特征和表達特性，為進一步闡明該基因功能奠定了重要的基礎。

關鍵詞：巴西橡膠樹;轉(zhuǎn)錄因子;CAMTA;基因表達

中圖分類號：S794.1? ? ? 文獻標識碼：A

Genome-wide Identification and Expression Analysis of the CAMTA Family in Rubber Tree （Hevea brasiliensis）

LIN Xianzu1，2， XIAO Xiaohu2， YANG Jianghua2， QIN Yunxia2， LONG Xiangyu2， FANG Yongjun2*

1. College of Tropical Crops， Hainan University， Haikou， Hainan 570228， China; 2. Key Laboratory of Biology and Genetic Resources of Rubber Tree， Ministry of Agriculture and Rural Affairs / Institute of Rubber Research， Chinese Academy of Tropical Agricultural Sciences， Haikou， Hainan 571101， China

Abstract： Calmodulin-binding transcription activators （CAMTAs） are important proteins with a conserved calmodulin- binding domain， which are widely involved in biological processes such as stress response， growth and development， and hormone signalling in plant. In this study， we identified eight HbCAMTAs screening the genome-wide of the rubber tree and analyzed physicochemical properties， gene structures， regulatory elements， gene expression patterns， and so on. The results showed that all HbCAMTAs were located in the nucleus， and predicted as hydrophilic proteins with molecular weight from 104.34 to 121.99 kD and isoelectric points between 5.32 and 7.83. The sequence homology analysis showed that HbCAMTAs could be divided into four subclades， and each one had similar gene structure， conserved domains. Furthermore， each HbCAMTA was potentially regulated by at least three miRNAs. Expression analysis indicated that the expression of HbCAMTAs differed significantly in different tissues， and the expression was continuously up-regulated in leaves with the development of leaves and reached the highest expression level at mature or palegreen. Except for HbCAMTA2， ethylene stimulation failed to significantly induce the expression of other CAMTA genes in rubber tree while low-temperature could significantly induce the up-regulated expression of HbCAMTAs. HbCAMTA1 was constitutively expressed in all samples， and also mostly expressed in latex， leaves， roots， bark， male and female flowers， and four different leaf developmental stages （Bronzite， Changecolor， Palegreen， and Mature）， which suggesting its fundamental function in rubber tree. HbCAMTA3 was expressed at the highest level in the seeds and was significantly up-regulated by low temperature stress， indicating its important role in response to low temperature stress. HbCAMTA7 was undetectably express in all investigated samples including different tissues， leaf developments， and treatments. Finally， this study would provide important information for further revealing the function of HbCAMTAs in rubber tree.

Keywords： Hevea brasiliensis; transcription factor; CAMTA; gene expression

DOI： 10.3969/j.issn.1000-2561.2021.10.015

鈣離子（Ca2+）是一種廣泛存在于植物中的重要的第二信使，在植物的生長發(fā)育和脅迫響應的信號轉(zhuǎn)導中發(fā)揮著重要的作用。鈣調(diào)蛋白（Calmodulin，CaM）作為植物中重要的鈣離子（Ca2+）結合蛋白之一，廣泛存在于植物當中。CaM在細胞中主要負責感知Ca2+濃度變化，進而達到傳遞Ca2+信號的目的，從而在諸多生理生化調(diào)節(jié)過程中扮演重要的角色[1]。Ca2+-CaM信號系統(tǒng)在植物調(diào)控抗病抗逆過程中起著重要作用，參與活性氧的產(chǎn)生[2]、過敏性壞死反應[3]、蛋白質(zhì)磷酸化過程及PR基因表達[4]的調(diào)節(jié)。

目前，植物中已經(jīng)鑒定出90多個轉(zhuǎn)錄因子為Ca2+/CaM結合蛋白，其中包括CAMTA、MYB、WRKY、NAC、bZIP等基因家族[5-9]。鈣調(diào)蛋白結合的轉(zhuǎn)錄激活子（CaM-binding transcription activator，CAMTA）是一類能夠與CaM結合的轉(zhuǎn)錄因子[4， 10-11]，廣泛存在于植物中，在生長發(fā)育、生物脅迫和非生物脅迫響應，特別是在冷脅迫方面起重要作用。模式植物擬南芥中研究表明AtCAMTA3具有包括CBF2在內(nèi)的10個相互作用因子，正調(diào)控參與冷應激基因CBF2的表達[12]，AtCAMTA1和AtCAMTA3雙突變體對冷凍的耐受性明顯降低[13]。玉米中研究表明ZmCAMTA4a、ZmCAMTA7a和ZmCAMTA7b通過冷脅迫處理在芽中顯著上調(diào)表達，只有ZmCAMTA4a在根中受冷脅迫誘導表達[14]。

巴西橡膠樹是目前商業(yè)天然橡膠的重要來源，在生產(chǎn)過程中受到多種脅迫威脅，其中低溫寒害是我國作為非傳統(tǒng)植膠區(qū)面臨的最大危害。低溫會引起橡膠樹長時間排膠，形成“長流”現(xiàn)象?！伴L流”易誘發(fā)橡膠樹發(fā)生死皮，不利于天然橡膠的高產(chǎn)、穩(wěn)產(chǎn)。在模式植物擬南芥中的研究表明，CAMTA基因家族成員參與了低溫脅迫的應答過程并調(diào)控相關冷應激基因的表達。而在橡膠樹中并未見相關研究報道，因而對于橡膠樹CAMTA基因研究有助于揭示CAMTA基因在橡膠樹低溫脅迫響應中的作用，所以了解橡膠樹抗寒機制具有重要的意義。本文以橡膠樹基因組和轉(zhuǎn)錄組數(shù)據(jù)為基礎，全面鑒定了橡膠樹8個CAMTA基因家族成員，利用生物信息學方法分析CAMTA家族成員的基因結構、保守結構域、啟動子元件、亞細胞定位、蛋白質(zhì)的理化性質(zhì)、家族成員的系統(tǒng)進化和表達水平。本研究將為進一步深入解析橡膠樹CAMTA基因功能奠定重要的基礎。

1? 材料與方法

1.1? HbCAMTA基因家族成員全基因組鑒定

參考PlantTFdb提供的CAMTA轉(zhuǎn)錄因子家族定義規(guī)則（http：//planttfdb.gao-lab.org/help_ famsc-hema.php），獲得用于該家族成員鑒定的DNA結合域CG1（PF03859），然后從PFAM數(shù)據(jù)庫（http：//pfam.xfam.org/）下載該結合域的隱馬爾可夫模型（hidden Markov model， HMM）文件，利用HMMER3.3（http：//hmmer.org/）程序搜索橡膠樹基因組數(shù)據(jù)庫中所有蛋白序列，從全基因組水平預測橡膠樹所有CAMTA基因家族成員。同時，從PlantTFdb數(shù)據(jù)庫下載擬南芥CAMTA蛋白序列（AtCAMTA1：AT5G09410.1、AT5G09410.2、AT5G09410.3;AtCAMTA2：AT5G64220.1、AT5G64220.2;AtCAMTA3：AT2G22300.1、AT2G22300.2;AtCAMTA4：AT1G67310.1;AtCAMTA5：AT4G16150.1;AtCAMTA6：AT3G16940.1）作為參考序列，采用blastp程序到HeveaDB數(shù)據(jù)庫（http：//hevea.catas. cn/home/index）進行比對檢索同源基因。將2種方法得到的基因集進行去冗余后，進一步通過NCBI的Batch Web CD-search Tool （https：//www. ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi）進行保守結構域分析，最終將含有該家族典型的結構域的蛋白序列鑒定為橡膠樹CAMTA家族成員（HbCAMTA）。

1.2? HbCAMTA基因家族成員基本理化性質(zhì)預測

利用在線網(wǎng)站Plant-mPLoc（http：//www.csbio. sjtu.edu.cn/bioinf/plant-multi/）預測HbCAMTA的亞細胞定位;利用在線工具ExPASy ProtParam（https：//web.expasy.org/protparam/）計算蛋白質(zhì)的等電點和分子量。

1.3? HbCAMTA基因結構、保守結構域和Motif的分析

根據(jù)橡膠樹基因組的注釋信息，利用TBtools（https：//github.com/CJ-Chen/TBtools）繪制HbCAMTA基因結構圖;利用NCBI Batch CD- search（https：//www.ncbi.nlm.nih.gov/Structure/ bwrpsb/bwrpsb.cgi）搜索序列保守結構域，然后通過IBS （http：//ibs.biocuckoo.org/）進行圖形化展示;利用MEME（http：//meme-suite.org/tools/ meme）在線服務預測CAMTA家族基因的保守motif。

1.4? 多序列聯(lián)配和系統(tǒng)進化分析

從PlantTFdb（http：//planttfdb.gao-lab.org/ index.php）數(shù)據(jù)庫和NCBI（https：//www.ncbi.nlm. nih.gov/）數(shù)據(jù)庫下載擬南芥、蓖麻、木薯、胡楊、毛果楊、麻風樹和亞麻的CAMTA蛋白序列與HbCAMTA組成分析蛋白序列集，先用Cluster軟件進行多序列比對，再利用MEGA X軟件采用鄰接法（neighbor-joining method，NJ）進行系統(tǒng)進化樹構建，bootstrap取1 000。

1.5? 調(diào)控HbCAMTA基因的miRNA預測

通過psRNATarget在線軟件（http：//plantgrn. noble.org/psRNATarget/）對調(diào)控橡膠樹CAMTA基因家族成員的miRNA進行預測，設定參數(shù)Expect值為3.5，UPE值為15。預測可能參與調(diào)控HbCAMTA的miRNA。

1.6? HbCAMTA家族成員表達分析

利用本實驗室已有的Solexa轉(zhuǎn)錄組數(shù)據(jù)庫，對不同組織、乙烯處理不同時間、低溫脅迫處理不同時間和葉片不同發(fā)育時期的HbCAMTA基因家族成員的表達模式進行分析，并以熱圖形式對結果進行展示。

1.7? HbCAMTA家族成員啟動子順式作用元件分析

截取HbCAMTA轉(zhuǎn)錄起始位點上游2 kb的序列作為啟動子序列進行分析，利用在線軟件PlantCARE（http：//bioinformatics.psb.ugent.be/webt?o-ols/plantcare/html/）預測啟動子序列的順式作用元件，通過TBtools軟件對預測的結果進行可視化展示。

2? 結果與分析

2.1? HbCAMTA家族成員鑒定與理化性質(zhì)分析

用隱馬爾可夫模型預測得到8個橡膠樹候選CAMTA基因家族成員，利用蛋白序列同源比對獲得候選HbCAMTA基因家族成員14個，去冗余后共16個，進一步通過保守結構域分析，最終鑒定獲得8個橡膠樹CAMTA成員基因，分別命名為HbCAMTA1-8。橡膠樹CAMTA編碼的蛋白序列長度基本相似，在923～1086個氨基酸之間，HbCAMTA5編碼的蛋白序列最長，達到了1086個氨基酸，而HbCAMTA1編碼的蛋白序列最短，只有923個氨基酸。利用ExPASy ProtParam在線工具對橡膠樹CAMTA基因編碼蛋白的理化性質(zhì)進行分析，結果表明HbCAMTA蛋白的理化特性具有多樣性，其分子量在104.34～121.99 kD之間，等電點在5.32～7.83之間，其中7個HbCAMTA蛋白等電點小于7.0，只有HbCAMTA2的等電點大于7，說明多數(shù)HbCAMTA基因編碼弱酸性蛋白（表1）。亞細胞定位的預測結果顯示所有CAMTA蛋白均不含信號肽，定位在細胞核。親水性平均系數(shù)是衡量蛋白親疏水性的重要指標，預測結果顯示8個橡膠樹HbCAMTA蛋白的親水性平均系數(shù)在-0.574～-0.384之間，表明橡膠樹所有CAMTAs蛋白都是親水性蛋白。

2.2? HbCAMTA基因結構和保守結構域分析

根據(jù)Tang等[15]發(fā)表的基因組結構注釋信息，利用TBtools軟件分析HbCAMTA家族成員的基因結構并繪制該基因家族的基因結構圖（圖1），我們發(fā)現(xiàn)鑒定到HbCAMTA成員均含有11～12個內(nèi)含子不等。HbCAMTA1和HbCAMTA7，?HbCAMTA3和HbCAMTA8，HbCAMTA4和HbCAMTA6組成了3個同源基因?qū)Γ總€同源基因?qū)邢嗤耐怙@子數(shù)量和內(nèi)含子數(shù)量，且每個外顯子長度相似，而內(nèi)含子長度差異明顯，主要體現(xiàn)在第1個和第3個內(nèi)含子，同源基因?qū)χg在基因結構上具有較高的相似性，暗示著CAMTA基因結構具有保守性。保守結構域分析表明，HbCAMTA1、HbCAMTA2、HbCAMTA4、HbCAMTA6、HbCAMTA8五個基因均含有保守的CG-1結構域、IQ結構域、Ank_2結構域和TIG結構域，HbCAMTA3含有保守的CG-1結構域、Ank_2結構域和TIG結構域，不含有IQ結構域，HbCAMTA5含有保守的CG-1結構域、IQ結構域和TIG結構域，不含有Ank_2結構域，HbCAMTA7含有保守的CG-1結構域、IQ結構域和ANKYR結構域，不含有TIG結構域（圖2）。

2.3? HbCAMTA家族成員的保守基序分析

利用在線工具MEME Suite 5.1.1對HbCAMTA蛋白的保守基序進行分析的結果顯示，HbCAMTA蛋白包含了10個保守的基序，分別命名為Motif 1-10。在10個基序中，Motif 6是未知的，Pfam數(shù)據(jù)庫未能找到對應結構域，該基序的功能還有待進一步的實驗證明。Motif 1、Motif 4和Motif 10與CG-1結構域相關，而Motif 3與ANK結構域相關。Motif 2和Motif 7與IQ結構域相關，而Motif 5和Motif 8與TIG結構域相關。圖3顯示了所有已確定的基序的序列標志圖。

2.4? HbCAMTA家族成員進化分析

從PlantTFdb和NCBI數(shù)據(jù)庫分別下載木薯（6個）、蓖麻（5個）、胡楊（7個）、毛果楊（7個）、麻風樹（4個）、擬南芥（6個）和亞麻（9個）的CAMTA氨基酸序列，結合本研究鑒定到的8個橡膠樹CAMTA氨基酸序列，組成共計52個植物CAMTA蛋白的序列集，利用MEGA X構建系統(tǒng)進化樹，結果如圖4顯示，CAMTA基因家族可以分為Ⅰ、Ⅱ、Ⅲ和Ⅳ 4個亞組。HbCAMTA家族成員在4個亞組均有分布，且在同一分枝上與木薯MeCAMTA基因聚在一起，這與物種系統(tǒng)進化上橡膠與木薯親緣關系最近相一致。Ⅰ亞組主要包括HbCAMTA4、MeCAMTA4、AtCAMTA4等;Ⅱ亞組主要包括HbCAMTA1、HbCAMTA7、MeCAMTA1、AtCAMTA5、AtCAMTA6等;Ⅲ亞組主要包括HbCAMTA2、MeCAMTA2等;Ⅳ亞組主要包括HbCAMTA3、HbCAMTA8、HbCAMTA5、MeCAMTA3、MeCAMTA5、AtCAMTA1、AtCAMTA2、AtCAMTA3等。另外，系統(tǒng)進化顯示存在HbCAMTA1和HbCAMTA7，HbCAMTA3和HbCAMTA8，HbCAMTA4和HbCAMTA6三個同源基因?qū)Γ@與基因結構和序列同源性分析結果一致。

2.5? 調(diào)控HbCAMTA基因的miRNA分析

利用在線軟件psRNATarget，以HbCAMTA 基因的cDNA序列作為靶基因，對植物的small RNA數(shù)據(jù)庫進行檢索，結果見表2，HbCAMTA基因家族成員均有3個及以上的miRNA的結合位點。HbCAMTA6只有3個miRNA的結合位點，而HbCAMTA2和HbCAMTA8各有10個miRNA的結合位點。HbCAMTA3受到9個miRNA調(diào)控，?擬南芥：AtCAMTA1-6;橡膠樹：HbCAMTA1-8;蓖麻：RCOM_0213350、RCOM_0999870、RCOM_1002900、RCOM_1079700、RCOM_1296350;毛果楊：PtCAMTA1-7;胡楊：CCG000518.1、CCG020006.1、CCG020717.1、CCG023491.1、CCG025222.1、CCG025313.1、CCG031532.1;木薯：MeCAMTA1-6;亞麻：LuCAMTA1-9;麻風樹：JCGZ_19069、JCGZ_18971、JCGZ_24032、JCGZ_20177。

其中最多的是屬于miR159家族的成員。結果暗示了miRNA可能在參與調(diào)控HbCAMTA的表達方面具有重要作用。

2.6? HbCAMTA家族成員轉(zhuǎn)錄表達分析

為了解橡膠樹CAMTA基因的表達特性，我們利用本實驗室已有的Solexa轉(zhuǎn)錄組測序數(shù)據(jù)，對不同組織、葉片不同發(fā)育時期、低溫脅迫處理以及乙烯處理膠乳中的家族成員進行表達分析。結果顯示（圖5），HbCAMTA1是橡膠樹CAMTA家族主要功能基因，在所檢測的7個組織（膠乳、葉片、根、雄花、雌花、樹皮、種子）、葉片的4個發(fā)育時期（古銅期、變色期、淡綠期、成熟期）、低溫脅迫處理3個時間點以及乙烯處理4個時間點均有表達，且整體表達水平較高。特別是在葉片發(fā)育過程中和低溫脅迫處理下明顯上調(diào)表達，推測該基因為組成型表達，且在葉片發(fā)育過程中和響應低溫脅迫起重要作用。HbCAMTA3在種子中表達水平最高，其次是樹皮中，其受低溫脅迫誘導明顯地持續(xù)上調(diào)表達且最顯著。在乙烯處理膠乳中HbCAMTA整體表達水平相對較低，HbCAMTA1呈現(xiàn)下調(diào)表達的趨勢，HbCAMTA2和HbCAMTA8呈現(xiàn)明顯的上調(diào)表達。低溫脅迫下，HbCAMTA1、HbCAMTA2、HbCAMTA3、HbCAMTA4、HbCAMTA5、HbCAMTA6、HbCAMTA8均受低溫誘導上調(diào)表達。HbCAMTA7基因在所有檢測樣品中表達水平均為最低甚至不表達，從該結果可以推測HbCAMTA1和HbCAMTA7在進化過程中出現(xiàn)了功能分化，HbCAMTA1在橡膠樹生長發(fā)育等方面具有重要功能，而HbCAMTA7出現(xiàn)了去功能化。

綜上結果表明，橡膠樹HbCAMTA基因家族成員在組織表達、乙烯處理的響應和葉片不同發(fā)育時期的表達方面存在明顯差異，反映了橡膠樹HbCAMTA基因功能的多樣性，而在低溫脅迫的表達方面具有一定的一致性。

2.7? HbCAMTA家族成員啟動子順式作用元件的分析

分別截取8個HbCAMTA基因轉(zhuǎn)錄起始位點上游2000 bp序列作為啟動子區(qū)序列，利用PlantCARE在線工具對啟動子區(qū)上可能存在的順式作用元件進行預測，得到位于啟動子區(qū)域的多種順式作用元件（圖6），包括脫落酸響應元件（14個）、逆境脅迫誘導（2個）、赤霉素誘導（8個）、茉莉酸甲酯響應元件（18個）、光響應元件（58個）、水楊酸響應元件（2個）、干旱響應元件（3個）、厭氧誘導元件（14個）等，表明HbCAMTA參與多種植物激素和環(huán)境脅迫應答。HbCAMTA中的6個基因（HbCAMTA1、HbCAMTA2、HbCAMTA3、HbCAMTA4、HbCAMTA7和HbCAMTA8）的啟動子區(qū)上擁有茉莉酸甲酯響應元件。除了HbCAMTA7之外，HbCAMTA中的7個基因在啟動子區(qū)上都有厭氧誘導元件。所有的HbCAMTA都擁有光響應元件。HbCAMTA2、HbCAMTA4、HbCAMTA5、HbCAMTA7和HbCAMTA8五個基因擁有脫落酸響應元件;HbCAMTA2和HbCAMTA5擁有干旱響應元件;HbCAMTA6和HbCAMTA7擁有水楊酸響應元件。

3? 討論

CaM是一種重要的Ca2+結合蛋白，通過生物化學、細胞生物學和分子生物學手段證明具有明確的生物學作用，是一種與多種靶蛋白結合的調(diào)節(jié)因子[11， 16-17]。CAMTA、WRKY、MYB和bZIP等轉(zhuǎn)錄因子已經(jīng)被證明通過與CaMs/CMLs相互作用來響應不同的生物和非生物脅迫[18-20]。

Latex：膠乳;Leaf：葉片;Root：根;Bark：樹皮;Female flower：雌花;Male flower：雄花;Seed：種子;ET-0h：乙烯處理0 h;ET-3h：乙烯處理3 h;ET-12h：乙烯處理12 h;ET-24h：乙烯處理24 h;Bronzite：古銅期;Changecolor：變色期;Palegreen：淡綠期;Mature：成熟期;4C-0h：低溫脅迫（4 ℃）處理0 h;4C-6 h：低溫脅迫（4 ℃）處理6 h;4C-12 h：低溫脅迫（4 ℃）處理12 h。

CAMTA轉(zhuǎn)錄因子在鈣/鈣調(diào)素轉(zhuǎn)導信號通路和CAMTA介導的基因轉(zhuǎn)錄調(diào)控中發(fā)揮關鍵作用，這是植物對外源激素和非生物脅迫做出反應的關鍵過程[21-24]。橡膠樹發(fā)育具有明顯的特異性（乳管），生產(chǎn)過程受到多種生物及非生物脅迫以及激素刺激，推測CAMTA可能發(fā)揮重要生物學功能。2016年Tang等[15]完成橡膠樹基因組測序，使得我們能夠?qū)bCAMTA進行全面地分析。在本研究中，我們對橡膠樹的CAMTA基因進行了全基因組水平的鑒定，并從理化性質(zhì)、基因結構、系統(tǒng)進化、調(diào)控因子以及轉(zhuǎn)錄表達模式等方面進行分析。

橡膠樹中鑒定到8個CAMTA基因，略多于擬南芥（6個）[25]、木薯（6個）[26]和毛果楊（7個）[27]的，而比同屬金虎尾目下的亞麻（9個）[28]少一個。我們采用基于隱馬爾可夫模型的HMMER軟件和blastp程序從全基因組水平進行鑒定橡膠樹的CAMTA成員，并取二者輸出的結果交集進一步進行保守結構域分析，最終將含有CG-1結構域的確認橡膠樹的CAMTA成員。HMMER軟件是基于隱馬爾可夫模型構建的HMM文件，利用概率推理方法和隱馬爾可夫模型（HMM）進行同源搜索，常用于同源蛋白的識別。同源過度擴展[29]和收斂進化[30]是序列比對程序錯誤地識別蛋白質(zhì)同源關系的重要因素，Mistry等[31]使用Bateman等[32]收集的配置文件-HMMs測試HMMER3分配蛋白質(zhì)序列到同源家族的準確性，結果表明假陽性率在1%～2%之間。同時，我們用擬南芥的CAMTA家族成員蛋白序列進行blastp，然后取HMMER和blastp輸出結果的交集，可有效避免因為程序算法引起的假陽性結果和假陰性結果，并通過保守結構域分析再次驗證了鑒定結果的可靠性。

基因結構分析發(fā)現(xiàn)橡膠樹CAMTA同源基因?qū)哂休^高的相似性，暗示著橡膠樹同擬南芥、毛果楊、玉米和木薯等物種一樣表現(xiàn)為基因結構具有較高的保守性[14， 25-27]。內(nèi)含子相對位置在同一分支中也具有保守性，而不同分支中則存在差異。系統(tǒng)進化分析顯示CAMTA可以分為4個亞組，其他物種一般分為4個或3個亞組，例如亞麻（3個）[28]、玉米（4個）[14]、毛果楊（3個）[27]、小麥（3個）[33]，擬南芥CAMTA在其他的物種系統(tǒng)進化分組中的分布與在本次系統(tǒng)進化分組中的分布是一致的，因此橡膠樹CAMTA分為4個亞組具有合理性。橡膠樹和木薯的CAMTA高度同源具有較近的親緣關系，這與二者同屬大戟科植物，進化距離較近相一致。HbCAMTA3和MeCAMTA3處在同一分枝上并與HbCAMTA5、HbCAMTA8、AtCAMTA1、AtCAMTA2和AtCAMTA3處在同一個亞組。該分支中的擬南芥CAMTA成員生物學功能已經(jīng)進行了充分研究，AtCAMTA3是一個響應低溫脅迫的轉(zhuǎn)錄因子，在低溫脅迫下通過正調(diào)控在植物低溫脅迫發(fā)揮重要作用的CBF2基因[21]和調(diào)控水楊酸（salicylic acid， SA）合成[13]發(fā)揮作用，而橡膠樹中HbCAMTA3受低溫脅迫誘導顯著持續(xù)上調(diào)表達，推測HbCAMTA3可能通過類似AtCAMTA3的生化反應在橡膠樹的低溫脅迫中發(fā)揮作用。同時AtCAMTA3還是一個鹽脅迫[34]響應因子，暗示了HbCAMTA3可能也具有功能多樣性。Galon等[10]通過基因芯片進行研究發(fā)現(xiàn)， AtCAMTA1突變體中有17個與生長素相關的基因表達上調(diào)，AtCAMTA1突變體和RNAi介導的CAMTA1抑制型轉(zhuǎn)基因株系的下胚軸伸長都表現(xiàn)出對生長素超敏的表型，認為擬南芥AtCAMTA1參與生長素信號傳導途徑并調(diào)節(jié)植物的生長發(fā)育。Yang等[35]在番茄中發(fā)現(xiàn)7個CAMTA基因在果實發(fā)育和成熟過程中差異表達，可能參與番茄果實的發(fā)育和成熟過程的調(diào)控，推測橡膠樹中CAMTA同源基因參與了組織和器官的生長發(fā)育過程。擬南芥中的研究表明AtCAMTA2是 AtALMT1（鋁激活的蘋果酸轉(zhuǎn)運蛋白1）的激活轉(zhuǎn)錄因子[36]，AtCAMTA3在SA介導的植物防衛(wèi)反應中發(fā)揮負向調(diào)節(jié)作用[23， 37-38]，AtCAMTA1、AtCAMTA2和AtCAMTA3 喪失阻遏作用能誘導植物防御基因和系統(tǒng)抗性的啟動[39]，木薯中MeCAMTA3通過調(diào)控木薯-黃單胞菌相互作用過程中的多種免疫反應調(diào)節(jié)木薯細菌性枯萎病的抗性[26]。因此，我們推測橡膠樹HbCAMTA3、HbCAMTA5和HbCAMTA8基因具有功能多樣性，可能參與激素信號轉(zhuǎn)導、生長發(fā)育調(diào)節(jié)、低溫脅迫響應、鹽脅迫和防衛(wèi)反應等過程，而其中HbCAMTA3基因是我們進一步研究橡膠樹響應低溫脅迫分子機制的重要候選基因之一。本研究首次對橡膠樹中CAMTA基因進行系統(tǒng)分析，將為后續(xù)進一步研究HbCAMTA基因在橡膠樹中的功能奠定基礎。

參考文獻

[1] Chin D， Means A R. Calmodulin： a prototypical calcium sensor[J]. Trends in Cell Biology， 2000， 10（8）： 322-328.

[2] Grant M， Brown I， Adams S， et al. The RPM1 plant disease resistance gene facilitates a rapid and sustained increase in cytosolic calcium that is necessary for the oxidative burst and hypersensitive cell death[J]. Plant Journal， 2000， 23（4）： 441- 450.

[3] Xu H， Heath M C. Role of calcium in signal transduction during the hypersensitive response caused by basidiospore- derived infection of the cowpea rust fungus[J]. The Plant Cell， 1998， 10（4）： 585-597.

[4] Kim M C， Chung W S， Yun D J， et al. Calcium and calmodulin-mediated regulation of gene expression in plants[J]. Molecular Plant， 2009， 2（1）： 13-21.

[5] Kim Y， Gilmour S J， Chao L， et al. Arabidopsis CAMTA transcription factors regulate pipecolic acid biosynthesis and priming of immunity genes[J]. Molecular Plant， 2020， 13（1）： 157-168.

[6] Pei J， Flieder D B， Patchefsky A， et al. Detecting MYB and MYBL1 fusion genes in tracheobronchial adenoid cystic carcinoma by targeted RNA-sequencing[J]. Modern Pathology， 2019， 32（10）： 1416-1420.

[7] Zhu D， Hou L， Xiao P， et al. VvWRKY30， a grape WRKY transcription factor， plays a positive regulatory role under salinity stress[J]. Plant Science， 2019， 280： 132-142.

[8] Yang Y， Yoo C G， Rottmann W， et al. PdWND3A， a wood-associated NAC domain-containing protein， affects lignin biosynthesis and composition in Populus[J]. BMC Plant Biology， 2019， 19（1）： 486.

[9] Lorenzo O. bZIP edgetic mutations： at the frontier of plant metabolism， development and stress trade-off[J]. Journal of Experimental Botany， 2019， 70（20）： 5517-5520.

[10] Galon Y， Finkler A， Fromm H. Calcium-regulated transcription in plants[J]. Molecular Plant， 2010， 3（4）： 653-669.

[11] Reddy A S N， Ali G S， Celesnik H， et al. Coping with stresses： Roles of calcium-and calcium/calmodulin-regulated gene expression[J]. The Plant Cell， 2011， 23（6）： 2010.

[12] Thomashow M F. Molecular basis of plant cold acclimation： insights gained from studying the CBF cold response pathway[J]. Plant Physiology， 2010， 154（2）： 571-577.

[13] Kim Y， Park S， Gilmour S J， et al. Roles of CAMTA transcription factors and salicylic acid in configuring the low- temperature transcriptome and freezing tolerance of Arabidopsis[J]. The Plant Journal， 2013， 75（3）： 364-376.

[14] Yue R， Lu C， Sun T， et al. Identification and expression profiling analysis of calmodulin-binding transcription activator genes in maize （Zea mays L.） under abiotic and biotic stresses[J]. Frontiers in Plant Science， Frontiers， 2015， 6（576）： 576.

[15] Tang CR， Yang M， Fang Y， et al. The rubber tree genome reveals new insights into rubber production and species adaptation[J]. Nature Plants， 2016， 2（6）： 16073.

[16] Kudla J， Batisti? O， Hashimoto K. Calcium signals： The lead currency of plant information processing[J]. The Plant Cell， 2010， 22（3）： 541-563.

[17] DeFalco T A， Chiasson D， Munro K， et al. Characterization of GmCaMK1， a member of a soybean calmodulin-binding receptor-like kinase family[J]. FEBS Letters， 2010， 584（23）： 4717-4724.

[18] Yang T， Poovaiah B W. A calmodulin-binding/CGCG box DNA-binding protein family involved in multiple signaling pathways in plants[J]. Journal of Biological Chemistry， 2002， 277（47）： 45049-45058.

[19] Finkler A， Ashery-Padan R， Fromm H. CAMTAs： calmodulin-binding transcription activators from plants to human[J]. FEBS Letters， 2007， 581（21）： 3893-3898.

[20] Yang T， Peng H， Whitaker B D， et al. Differential expression of calcium/calmodulin-regulated SlSRs in response to abiotic and biotic stresses in tomato fruit[J]. Physiologia Plantarum， 2013， 148（3）： 445-455.

[21] Doherty C J， Buskirk H A V， Myers S J， et al. Roles for arabidopsis CAMTA transcription factors in cold-regulated gene expression and freezing tolerance[J]. The Plant Cell， 2009， 21（3）： 972-984.

[22] Galon Y， Nave R， Boyce J M， et al. Calmodulin-binding transcription activator （CAMTA） 3 mediates biotic defense responses in Arabidopsis[J]. FEBS Letters， 2008， 582（6）： 943-948.

[23] Nie H， Zhao C， Wu G， et al. SR1， A calmodulin-binding transcription fFactor， modulates plant defense and ethylene-Induced senescence by directly regulating NDR1 and EIN3[J]. Plant Physiology， 2012， 158（4）： 1847-1859.

[24] Qiu Y， Xi J， Du L， et al. Coupling calcium/ calmodulin- mediated signaling and herbivore-induced plant response thro?ugh calmodulin-binding transcription factor AtSR1/ CAMT?A3[J]. Plant Molecular Biology， 2012， 79（1-2）： 89-99.

[25] Bouché N， Scharlat A， Snedden W， et al. A novel family of calmodulin-binding transcription activators in multicellular organisms[J]. Journal of Biological Chemistry， 2002， 277（24）： 21851-21861.

[26] Chang Y， Bai Y， Wei Y， et al. CAMTA3 negatively regulates disease resistance through modulating immune response and extensive transcriptional reprogramming in cassava[J]. Tree Physiology， 2020， 40（11）： 1520-1533.

[27] Wei M， Xu X， Li C. Identification and expression of CAMTA genes in Populus trichocarpa under biotic and abiotic stress[J]. Scientific Reports， 2017， 7（1）： 17910.

[28] Ali E， Raza M A， Cai M， et al. Calmodulin-binding transcription activator （CAMTA） genes family： Genome-wide survey and phylogenetic analysis in flax （Linum usitatissimum）[J]. PLoS One， 15（7）： e0236454.

[29] Gonzalez M W， Pearson W R. Homologous over-extension： a challenge for iterative similarity searches[J]. Nucleic Acids Research， 2010， 38（7）： 2177-2189.

[30] Pearson W R， Sierk M L. The limits of protein sequence comparison？[J]. Current Opinion in Structural Biology， 2005， 15（3）： 254-260.

[31] Mistry J， Finn R D， Eddy S R， et al. Challenges in homology search： HMMER3 and convergent evolution of coiled-coil regions[J]. Nucleic Acids Research， 2013， 41（12）： e121.

[32] Bateman A， Coin L， Durbin R， et al. The Pfam protein families database[J]. Nucleic Acids Research， 2004， 32（suppl_1）： D138-D141.

[33] Yang F， Dong F， Hu F， et al. Genome-wide identification and expression analysis of the calmodulin-binding transcription activator （CAMTA） gene family in wheat （Triticum aestivum L.）[J]. BMC Genetics， 2020， 21（1）： 105.

[34] Prasad K V S K， Abdel-Hameed A A E， Xing D， et al. Global gene expression analysis using RNA-seq uncovered a new role for SR1/CAMTA3 transcription factor in salt stress[J]. Scientific Reports， 2016， 6（1）： 27021.

[35] Yang T， Peng H， Whitaker B D， et al. Characterization of a calcium/calmodulin-regulated SR/CAMTA gene family during tomato fruit development and ripening[J]. BMC Plant Biology， 2012， 12（1）： 19.

[36] Tokizawa M， Kobayashi Y， Saito T， et al. STOP1， CAMTA2 and other transcription factors are involved in aluminum-inducible AtALMT1 expression[J]. Plant Physiology， 2015： 114.

[37] Du L， Ali G S， Simons K A， et al. Ca2+/calmodulin regulates salicylic-acid-mediated plant immunity[J]. Nature， 2009， 457（7233）： 1154-1158.

[38] Zhang L， Du L， Shen C， et al. Regulation of plant immunity through ubiquitin-mediated modulation of Ca2+-calmodulin- AtSR1/CAMTA3 signaling[J]. The Plant Journal， 2014， 78（2）： 269-281.

[39] Kim Y， Gilmour S J， Chao L， et al. Arabidopsis CAMTA transcription factors regulate pipecolic acid biosynthesis and priming of immunity Genes[J]. Molecular Plant， 2020， 13（1）： 157-168.

責任編輯：孫海燕