全基因組水平蒺藜苜蓿CPP 基因家族的鑒定及表達模式分析

田驕陽，王秋霞，鄭淑文，劉文獻

（蘭州大學草地農業生態系統國家重點實驗室，蘭州大學農業農村部草牧業創新重點實驗室，蘭州大學草地農業教育部工程研究中心，蘭州大學草地農業科技學院，甘肅 蘭州 730020）

轉錄因子（transcription factors）是一類能夠與DNA 特異序列結合并調節翻譯過程的蛋白質，由轉錄調控區、核定位信號區、寡聚位點和DNA 結合區組成，通過4 個區域相互協同，激活或抑制下游基因轉錄［1］。轉錄因子能夠參與環境脅迫響應、代謝調節、物質生成等生命過程［2］，在動植物生長發育及應激反應中具有重要作用，是深入解析生命活動分子調控機制的關鍵環節。CPP（cysteine-rich polycomb-like protein）基因家族是一類小轉錄因子，廣泛存在于除酵母與原核生物外的各種生物中［3］。CPP 轉錄因子中包含一個或兩個富含半胱氨酸的CXC 功能域（PF03638），保守序列的結構為CXCX4CX3YCXCX6CX3CXCX2C，并被一個含有短的保守序列RNPXAFXPK 的基序分隔［3］。CXC 結構域可以通過結合DNA 對靶基因起調控作用［4］。突變的CXC 保守結構域失去與DNA 結合活性，抑制細胞周期進程，導致細胞核形態異常［5］。CXC 結構域在不同植物中高度保守，可通過DNA 結合來誘導基因的抑制或表達［6］。隨著植物基因組測序的相繼完成，在全基因組水平對CPP基因家族的鑒定和分析成為可能。截至目前，已在擬南芥（Arabidopsis thaliana）［7］、水稻（Oryza sativa）［7］、大豆（Glycine max）［8］、玉米（Zea mays）［9］、黃瓜（Cucumis sativus）［10］、茶樹（Camellia sinensis）［11］等物種中分別鑒定到8、11、20、13、5、11 個CPP基因，從而為深入解析該基因家族的生物學功能奠定了基礎。

前人研究發現，CPP基因家族在植物生長發育及響應脅迫過程中具有重要作用。例如，TSO1是在擬南芥中發現的第一個CPP 轉錄因子，在調節細胞分裂以及開花中起到重要作用［12-13］。AtTSO1編碼一個花特異性的細胞分裂成分，該基因突變會導致花分生組織細胞分裂過程受到影響，表現為細胞壁形成受阻以及DNA 倍性增加，導致突變體花的細胞核大小和形狀不規則［12］。另外，TSO1 和MYB3R1 形成一個調節模塊，協調芽和根組織細胞增殖與分化，造成莖和根分生組織收縮，從而導致地上部分和根部發育嚴重異常［11］。大豆中的GmCPP1能夠與血紅蛋白基因Gmlbc3的啟動子相互作用，參與共生根瘤中大豆血紅蛋白基因的調控［14］。GhCPP基因在陸地棉（Gossypium hirsutum）的雌蕊中高度表達，可能參與了花器官和雌蕊發育的調控過程。此外，CPP基因家族也在響應多種應激脅迫過程中發揮功能。通過qRT-PCR 對13 個ZmCPP基因的表達分析表明，大多數ZmCPP基因的表達在熱應激24 h 和冷應激12 h 后達到峰值［9］。大豆中GmCPP 轉錄因子除GmCPP03和GmCPP07外，其余18 個基因均在調節熱應激反應中起重要作用［8］。另外，黃瓜CsCPP基因能夠以依賴脫落酸的方式提高植物對逆境的耐性［10］。

蒺藜苜蓿（Medicago truncatula）是一種豆科模式植物，具有基因組?。?70 Mb）、自花授粉、植株再生時間短以及突變體多樣的特點［15］。另外，蒺藜苜蓿與大多數豆科植物有著較高的遺傳相似性，通過研究蒺藜苜?？梢詾槠渌箍浦参铮貏e是紫花苜蓿（Medicago sativa）提供理論依據［16］。隨著蒺藜苜?；蚪M測序的完成，多種蒺藜苜?；蚣易澹鏛BD 家族［17］、AQPs 家族［18］、FAD 家族［19］被鑒定和系統研究，從而為深入解析蒺藜苜蓿生長發育調控及逆境脅迫響應過程提供了重要基因資源。截至目前，蒺藜苜蓿CPP基因家族的相關研究尚未見報道。本研究以蒺藜苜蓿基因組數據為基礎，利用生物信息學方法對全基因組水平的蒺藜苜蓿CPP基因家族進行鑒定，進而對不同家族成員的理化性質、結構特征、順式作用元件、進化關系及壓力、染色體位置等進行系統分析。此外，對CPP基因家族在蒺藜苜蓿不同組織及非生物脅迫下的表達模式進行研究，以期為后期深入研究蒺藜苜蓿CPP基因家族功能以及通過基因工程技術創制高抗逆苜蓿新種質提供理論依據。

1 材料與方法

1.1 MtCPP 基因家族成員鑒定

本試驗于2020 年7 月-2021 年3 月進行。根據Pfam（http：//pfam. sanger. ac. uk/search）［20］數據庫獲得的CPP 的保守結構域（PF03638），從Phytozome 網站（https：//phytozome.jgi.doe.gov/pz/portal.html）下載蒺藜苜蓿CPP基因家族的蛋白質編碼區(coding sequence, CDS)和蛋白序列。通過HMMER（http：//www.ebi.ac.uk/）［21］驗 證 保 守 結 構 域。在ExPASy 網 站（http：//www. expasy. org/）分 析 等 電 點（isoelectric point），分 子 量（molecular weight）和親水指數（grand average of hydropathicity）等參數［22］。通過WoLF-PSORT（https：//www.genscript.com/wolf-psort.html）網站對MtCPP 蛋白亞細胞定位進行預測分析［23］。

1.2 多序列比對、系統發育及基因進化壓力分析

擬南芥、水稻及大豆的CPP 蛋白序列分別從http：//www. arabidopsis. org/，http：//rice. plantbiology. msu.edu/，http：//phyozome. jgi. doe. gov/soybean 下載。利用MEGA 7.0 軟件（http：//www. megasoftware. net/）及Bioedit 軟件構建系統發育樹，方法為最大似然法（maximum likelihood estimate），設置最大循環數為1000［24］。利用MCScanX 和TBTools 軟件進行基因區塊識別并計算同義替換率(the ratio of the number of synonymous substitutions per synonymous site,Ks)、非同義替換率（the ratio of the number of nonsynonymous substitutions per nonsynonymous site,Ka）［25-26］。

1.3 染色體定位、結構分析、順式作用元件及保守結構域鑒定

利用Mapchart 軟件對MtCPP基因染色體定位進行作圖。通過GSDS 2.0 網站（http：//gsds.cbi.pku.edu.cn）分析蒺藜苜蓿的外顯子-內含子結構。MtCPP基因的啟動子區（翻譯起始密碼子上游2000 bp）從蒺藜苜?；蚪M數據中獲得，利用PlantCARE（http：//bioinformatics.psb.ugent.be/webtools/plantcare/html/）對啟動子順式作用元件進行分析。利用MEME（http：//meme.nbcr.net/meme/cgibin/meme.cgi）分析蒺藜苜蓿保守結構域，采用TBtools 工具進行作圖［27］。

1.4 MtCPP 基因家族表達模式及響應非生物脅迫分析

為分析MtCPP基因在不同器官及非生物脅迫下的表達模式，從Noble 網站（https：//mtgea. noble.org/v3/blastt_search_form.php）中獲取不同MtCPP基因探針，并獲得了對應探針在葉、葉柄、葉芽、幼苗、莖、花、莢果、種衣、種子、根等組織中以及干旱脅迫和鹽脅迫下的表達量。干旱處理是對幼苗和根組織干旱脅迫2、3、4、7、10、14 d 以及解除脅迫后1、2、4 d。鹽處理是利用200 mmol·L-1NaCl 對根 組織處理1、2、5、10 和24 h。最后使用TBtools 軟件構建表達熱圖。

2 結果與分析

2.1 MtCPP 基因家族的鑒定和分析

通過Pfam 結構域搜索最終在蒺藜苜?；蚪M共鑒定出9 個CPP基因家族成員，分別重新命名為MtCPP1~MtCPP9。利用在線網站對MtCPP基因編碼蛋白的長度、分子量、理論等電點、親水性和亞細胞定位進行分析。MtCPP基因的編碼序列長度介于924（MtCPP5）~2604 bp（MtCPP1），蛋白質由307（MtCPP5）~867 AA（MtCPP1）個 氨基酸 組 成，分 子 量 為34610.89~93893.24 Da，等 電 點介于4.60~8.93，其 中MtCPP1 和MtCPP3 等電點小于6.5，為酸性氨基酸；MtCPP5、MtCPP6 和MtCPP9 等電點大于8，為堿性氨基酸。此外，MtCPP家族的蛋白質親水性范圍為-0.802（MtCPP4）~-0.169（MtCPP5），9 個蛋白質均為親水性蛋白質。根據WoLF PSORT 預測，除MtCPP5為葉綠體基因外，其他8 個基因都為核基因（表1）。

表1 MtCPP 家族基本信息Table 1 The basic information of MtCPP gene family in M.truncatula

2.2 MtCPP 基因系統進化分析

為了研究MtCPP基因家族成員的系統進化特征，利用MEGA 7.0 軟件構建了蒺藜苜蓿與擬南芥、水稻和大豆CPP 家族之間的系統發育樹（圖1）。基于系統發育樹，可將CPP 蛋白分為A、B、C 三組。A 組中包含7 個GmCPP 蛋白、3 個MtCPP 蛋白及1 個OsCPP 蛋白。在A 組中，除了OsCPP4 來源于單子葉植物外，其余蛋白都來源于雙 子葉植物。B 組包含3 個MtCPP 蛋 白、6 個GmCPP 蛋 白、5 個OsCPP 蛋白及4 個AtCPP 蛋白。C 組含有3 個MtCPP 蛋白、7 個GmCPP 蛋白、5 個OsCPP 蛋白及4 個AtCPP 蛋白。由以上 結果發現，MtCPP基因均勻地分布于3 組中；3 組中均含有單子葉與雙子葉基因，表明CPP基因是在單子葉和雙子葉分化前形成的［7］。在B和C 組中，蒺藜苜蓿的CPP 基因與大豆的遺傳距離較近，如B 組中MtCPP6 與GmCPP5、GmCPP7，C 組中MtCPP1 與GmCPP16、GmCPP17、GmCPP18，表明雙子葉和單子葉植物在分化后沿著不同方向進化［7］。研究表明，A 組中的GmCPP1 可以調節豆血紅蛋白的合成［14］，因此，同組的蒺藜苜蓿MtCPP基因可能具有類似的功能。C 組中的AtCPP4 和AtCPP5 分別對應SOL1、TSO1基因的蛋白序列，在花器官中表達，并調控花的發育［3］，因此C 組中的MtCPP基因功能可能與花的發育有關。

圖1 MtCPP 基因家族進化分析Fig. 1 The phylogenetic analysis of MtCPP gene family

2.3 MtCPP 基因進化壓力分析

為了探究MtCPPs 的進化模式、選擇壓力，對MCScanX 篩選出的兩對MtCPP同源基因（MtCPP6和MtCPP9、MtCPP8和MtCPP9）計算非同義/同義替換比（Ka/Ks），結果如表2 所示，兩對MtCPP基因的Ka都明顯小于Ks，Ka/Ks都小于1，表明在MtCPP基因進化過程中，主要是凈化選擇發揮作用［28］。

表2 MtCPP 基因進化壓力分析Table 2 Evolutionary stress analysis of MtCPP genes

2.4 保守結構域及基因結構分析

根據蛋白序列分析發現（圖2），大部分MtCPP 中有兩個CXC 保守結構域以及在兩個保守結構域之間的R 基序；兩個CXC 結構域中都富含9 個保守的半胱氨酸（cysteine，Cys），在R 基序中含有相對保守的RNPXAFXPK 序列。有研究表明CXC 結構域可以結合特定DNA 序列［13］。另外，并不是所有MtCPP 都具有兩個CXC 保守結構域，其中MtCPP2 缺少CXC1 結構域，并且在CXC1 區域具有較多的冗余氨基酸序列，這與前人研究結果相同［7］。MtCPP5 缺少CXC2 部分序列；MtCPP3 的CXC2 保守結構域較短，其原因為N 端丟失了含有3 個保守半胱氨酸的11 個氨基酸序列；缺失和缺少CXC 結構域可能會影響基因功能的發揮。蒺藜苜蓿MtCPP 家族的CXC1 結構域保守性高于CXC2 結構域，這一結果與水稻［29］、擬南芥［29］、茶樹［11］等的研究結果一致。

圖2 MtCPP 氨基酸序列比對Fig. 2 Multiple alignments of the amino acid sequences of MtCPPs

為了進一步驗證CPP 轉錄因子在蒺藜苜蓿中的保守結構域，使用MEME 軟件對蒺藜苜蓿的7 個結構域進行分析（圖3）。結果表明，9 個MtCPP 蛋白的結構域長度介于11~49 AA，分別含有2~7 個保守結構域。除MtCPP2 以外，其他家族成員均含有保守結構域1（位于CXC1 的N 端）；除MtCPP5 以外，其他家族成員均含有保守結構域5（位于CXC1 的C 端）。保守結構域2 和7 在多個MtCPP 蛋白中存在。雖然MtCPP基因家族在進化過程中發生分化，但一些基序在特定的亞家族中具有保守性，例如，MtCPP6、MtCPP8、MtCPP9 三個成員之間以及MtCPP1、MtCPP4、MtCPP7 三個成員之間均具有相同的保守結構域。

圖3 基于綜合系統發育樹的MtCPPs 特征分析Fig. 3 Characteristic features of the MtCPPs based on the comprehensive phylogenetic tree

為進一步了解MtCPP基因的結構差異，進一步比較了MtCPP的全長基因序列及結構差異（圖4）。結果表明，MtCPP基因家族片段長度差異較大，MtCPP5的片段長度最短，MtCPP2的片段長度最長。MtCPP基因內含子數目在6~16 個之間，外顯子數目在7~17 個之間。同一亞組的內含子數目相近，如MtCPP6與MtCPP9都具有6 個內含子；MtCPP1和MtCPP7基因結構特征類似，但MtCPP7較MtCPP1缺少1 個內含子。

圖4 MtCPP 系統發育關系及基因結構分析Fig. 4 Phylogenetic relationships and gene structure of MtCPP genes

2.5 MtCPP 基因順式作用元件分析

轉錄因子通過與啟動子特定序列結合對基因轉錄的起始與表達進行調控，進而在調控植物生長發育及響應逆境脅迫過程中發揮重要作用。為了預測MtCPP基因家族成員在響應逆境脅迫方面的潛在功能，通過PlantCare 網站（http：//www.bioinformatics.psb.ugent.be/webtools/plantcare/html/）分析了不同MtCPP家族成員2000 bp 長度啟動子中的順式作用元件［30］。結果如圖5 所示，以激素調節和逆境脅迫為篩選條件，共在9 個基因啟動子中篩選了15 種順式作用元件，分別為無氧誘導（ARE、GC-motif）、水楊酸反應（TCA-element）、防御和壓力反應（TC-rich repeats）、茉莉酸甲酯反應（CGTCA-motif、TGACG-motif）、干旱誘導（MBS）、類黃酮生物合成基因調控（MBSI）、脫落酸反應（ABRE）、赤霉素反應（GARE-motif、P-box、TATC-box）、創傷反應（WUN-motif）以及生長素反應（AuxRR-core、TGA-element）。結果表明，平均每個基因中含有7 個順式作用元件。其中，MtCPP5和MtCPP6含有最多與激素調節和逆境脅迫相關的順式作用元件（13 個），MtCPP3中順式作用元件最少（2 個）；MtCPP1、MtCPP2、MtCPP5、MtCPP6中有較多的無氧誘導元件，MtCPP5和MtCPP7中有較多的創傷誘導元件?？刂祁慄S酮生物合成基因啟動子調控元件僅在MtCPP3中存在。具有干旱誘導元件的基因有MtCPP2、MtCPP5、MtCPP6、MtCPP7；MtCPP4、MtCPP5、MtCPP9具有脫落酸誘導元件，而MtCPP7和MtCPP8啟動子中具有生長素誘導元件。

圖5 MtCPP 基因啟動子順式作用元件分析Fig. 5 The analysis of cis-elements in the promoters of MtCPP genes

2.6 MtCPP 基因染色體定位

通過分析蒺藜苜蓿GFF3 文件獲取9 個MtCPP基因的染色位置進行染色體定位。結果如圖6 所示，9 個MtCPP定位在6 條染色體上。其中，4 個基因（MtCPP1、MtCPP2、MtCPP3、MtCPP4）定位在1 號染色體上；2、3、5、6、8 號染色體上都包含一個基因。相較于其他染色體，染色體1 上存在密度較大的基因簇，其中，MtCPP3與MtCPP4為串聯重復。基于基因重復分析，在MtCPP6與MtCPP9、MtCPP8和MtCPP9之間發生了兩個片段重復事件，預測這些基因在進化過程中可能來源于片段重復與串聯重復。

圖6 MtCPP 基因家族共線性分析Fig. 6 The collinearity analysis of MtCPP gene family

2.7 MtCPP 基因家族表達模式分析

為了進一步解析MtCPP基因的表達模式，基于蒺藜苜蓿探針提取了9 個MtCPP基因不同組織、相同組織不同時期以及逆境脅迫下的56 個表達芯片數據［31］并進行熱圖制作，基因表達強度用紅-黃-藍聚類圖指示。從紅色到黃色最后到藍色，基因的表達量逐漸降低，結果如圖7 所示。在營養器官及發育時期特異性方面，根據聚類分析可將9 個基因劃分為兩類（圖7A），MtCPP1~MtCPP4為Ⅰ類，MtCPP5~MtCPP9為Ⅱ類。Ⅰ類基因在根中相對表達量較高，表明Ⅰ類基因主要為根表達基因。Ⅰ類中只有MtCPP1基因在種衣中表達量較高，Ⅱ類中除MtCPP5基因外，其余4 個基因在種衣中表達量均較高。在莢果中，Ⅰ類基因表達量較高，Ⅱ類基因中除了MtCPP5，其余基因表達量均較低。MtCPP1、MtCPP2、MtCPP5和MtCPP8基因在種子授粉后表達量在種子中下降，MtCPP2在種子授粉后20 d，達到最低值；MtCPP1、MtCPP5和MtCPP8在種子授粉后36 d 達到最低值，表明這4 個基因可能主要在植物的營養生長過程中發揮作用。MtCPP3表達量在授粉后上升，可能在植物的生殖生長中發揮作用。MtCPP5基因主要在葉芽和種子中表達；MtCPP8基因隨著莖的發育表達量增加；MtCPP9基因在花、葉以及莢果中表達量較高。

圖7 MtCPP 基因不同組織及響應脅迫表達模式分析Fig. 7 Analysis of MtCPP gene expression profiles in different tissues and stress

水分脅迫處理基因表達數據分析結果顯示（圖7B），對幼苗和根脅迫處理后，MtCPP6和MtCPP7基因被誘導表達，并在解除脅迫后表達量降低，表明MtCPP6和MtCPP7基因在響應干旱脅迫方面可能發揮著重要作用。對根進行干旱脅迫后，MtCPP2和MtCPP5基因表達量升高，可能為根特異性干旱脅迫響應基因。MtCPP4、MtCPP5、MtCPP8、MtCPP9基因在受脅迫后表達量降低，解除脅迫后表達量升高，根中MtCPP1基因受干旱脅迫影響，表明這些基因受干旱脅迫的影響。MtCPP1和MtCPP2基因在幼苗中對干旱脅迫不敏感，MtCPP8基因對根的干旱脅迫不敏感，表明MtCPP基因在響應干旱過程中具有組織特異性。MtCPP2、MtCPP8在鹽脅迫誘導的條件下表達量升高，表明MtCPP2、MtCPP8是植物響應鹽脅迫的重要基因。MtCPP1、MtCPP5、MtCPP6、MtCPP7基因在鹽脅迫條件下表達量降低，表明鹽脅迫影響著這4 個基因的表達（圖7C）。

3 討論

CPP基因普遍存在于植物和動物基因組中，截至目前，CPP基因家族的一些成員已經在模式植物中進行了廣泛的研究，例如大豆和擬南芥［12，14］。并且發現CPP家族的基因在調控生長發育及響應逆境脅迫過程中具有重要的生物學功能［11］。然而，蒺藜苜蓿MtCPP基因家族尚未進行鑒定，其生物學功能仍不清晰。本研究利用生物信息學方法，在蒺藜苜蓿中共鑒定出9 個CPP家族基因，不均勻分布在蒺藜苜蓿6 條染色體上。其中MtCPP3與MtCPP4為串聯重復，MtCPP6與MtCPP9、MtCPP8與MtCPP9為片段重復。因此，蒺藜苜蓿CPP基因家族進化過程中，串聯重復以及片段重復可能起到了重要的促進作用。

3.1 MtCPP 保守結構域及進化關系

CPP基因家族廣泛存在著兩個CXC 保守結構域及保守結構域之間的保守R 基序［4］。其中，CXC 保守結構域廣泛存在于動物和植物中，但動物中的CPP同源基因缺失兩個CXC 保守基序之間的保守R 基序［3］。在兩個CXC結構域中具有9 個保守的半胱氨酸，C1 結構域的保守性大于C2 結構域，該特征廣泛存在于蒺藜苜蓿、水稻［7］、擬南芥［7］、黃瓜［10］等CPP基因家族中。多重序列比對結果顯示，并非所有的CPP 家族蛋白都有3 個完整的跨膜結構域。蒺 藜 苜 蓿MtCPP5 只 含 有C1 保 守 結 構 域，與 水 稻OsCPP4 相 同［10］；MtCPP7 只 含 有C2 結 構 域，與 小 麥（Triticum aestivum）中TaCPP3 和TaCPP4［32］相同；其余7 個MtCPP 均含有C1 和C2 結構域。較多的冗余氨基酸出現在MtCPP2 兩種蛋白C1 結構域的第3 和4 個半胱氨酸殘基之間，而MtCPP3 在C1 的N 端則缺失16 個氨基酸，包括3 個半胱氨酸殘基，與AtCPP6 研究結果類似［7］。目前，已有研究表明兩個CXC 結構域和兩個結構域之間的保守結構域間存在著高度相關，并共同進化，兩個結構域相互協同發揮特定功能，一個結構域中氨基酸序列的變化將導致另一個結構域中的反選擇或補償變化［33］，缺失某部分CXC 可能會導致兩個CXC 結構域之間的相互作用及其功能喪失。本研究中MtCPP6 和MtCPP9 來自片段重復，有著相同的保守結構域，并且它們的基因結構也非常相似，推測它們可能具有相同的功能。

為了進一步探究蒺藜苜蓿的進化關系，本研究通過與其他3 個已鑒定物種的CPP 蛋白的系統進化進行了分析，進而可將MtCPP分為A、B、C 三個類群。在A 組中，除1 個水稻基因外，其他都為大豆和蒺藜苜蓿CPP基因，表明同為豆科植物的CPP家族之間的親緣關系更為相近。其他兩組中均含有4 個物種基因，以上結果說明單子葉植物與雙子葉植物擁有共同的進化起源以及不同方向的進化選擇［7］。這種進化關系在水稻的SAUR基因家族［34］、擬南芥和水稻的Dof基因家族［29］等中都有類似報道。另外，進化關系相近的蛋白之間表現出更為相近的生物學功能。本研究中，A 組中的大豆CPP1 蛋白具有調節豆血紅蛋白合成的功能［14］，C 組中擬南芥的SOL1及TSO1在調節花器官的發育方面具有重要功能［12-13］，從而為與這些功能已知CPP基因進化關系更近的蒺藜苜蓿CPP基因的進一步生物學功能解析提供了參考和指導。

3.2 MtCPP 基因表達的時空特異性

赤霉素在響應干旱、水淹、鹽脅迫方面發揮著重要作用［35］。MtCPP5、MtCPP6和MtCPP7中含有多種與激素相關的順式作用元件，表明這3 個基因可能受多種激素調控響應逆境脅迫。CPP基因家族參與植物對多種非生物脅迫的反應。在本研究中，9 個MtCPP基因在干旱和鹽脅迫處理的不同時間點表現出不同的表達模式。在干旱脅迫下，MtCPP2、MtCPP5、MtCPP6和MtCPP7被誘導表達，其中MtCPP2和MtCPP5僅在根中誘導表達。前人研究中，干旱脅迫導致玉米中ZmCPP1.1、ZmCPP7.1、ZmCPP9、ZmCPP12在處理12 h 后上調表達［9］；大豆［8］、黃瓜［10］、茶樹［11］等研究中也表明CPP家族一些基因具有抗旱作用。啟動子是調控基因表達及響應逆境脅迫的重要順式調控元件。植物激素在響應脅迫方面發揮重要作用。干旱條件下，脫落酸調節植物氣孔的開閉并參與抗逆基因的表達調控。MtCPP5基因中含有脫落酸和干旱誘導順式作用元件，可能通過依賴脫落酸的方式發揮抗旱作用。對茶樹的CsCPP基因進行qRT-PCR 分析，6 個基因在干旱和脫落酸處理后上調表達，并且啟動子區分布有抗旱及脫落酸相關的順式作用元件［11］。黃瓜CsCPP基因的抗逆性依賴于脫落酸調節作用［10］。在外源脫落酸的作用下，橡膠樹（Hevea brasiliensis）葉片HbCPP1表達量升高［36］。這些結論表明，CPP基因家族中的部分基因可能通過脫落酸的調節誘導植物發揮抗旱功能。

CPP基因家族的一些基因與生殖相關。例如，擬南芥的TSO1影響著胚珠、花粉的生長；SOL1主要在花粉中表達［3］；黃瓜CsCPP4參與黃瓜的子房發育［10］；巴西橡膠樹HbCPP1基因隨著花器官的生長表達量增加［36］。本研究發現，MtCPP9在花器官及莢果中大量表達，并隨著生殖階段的變化而在不同生殖器官差異表達，MtCPP9可能參與蒺藜苜蓿花器官的發育以及生殖調控過程。B 類基因除MtCPP5以外，其他基因均可能參與種皮的形成，可能在生殖生長的后期發揮重要作用。另外，CPP基因在植物的營養生長方面也起到重要作用，例如，TSO1與MYB3R1結合共同調節擬南芥根和芽的分化［14］；北美鵝掌楸（Liriodendron tulipifera）中LtTCX2基因在葉中大量表達；大豆中GmCPP16在幼嫩葉片和根瘤中特異表達［8］。本研究發現MtCPP3、MtCPP8、MtCPP9分別在根、莖、葉中大量表達，暗示這些MtCPP基因可能對營養生長發育過程具有重要調控作用。

4 結論

本研究在全基因組水平首次鑒定分析了蒺藜苜蓿的9 個CPP基因，發現MtCPP基因中存在著CXC 保守結構域，可能為MtCPP基因發揮功能的關鍵因素。系統發育分析表明MtCPP基因與大豆CPP基因親緣關系較近。此外，9 個MtCPP基因表達模式具有一定的時空特異性，并具有響應干旱和鹽脅迫的特性。本研究結果可為進一步研究MtCPP在植物逆境響應中的調控機制和功能提供基礎，也可為在其他豆科植物，特別是在紫花苜蓿中鑒定和分析CPP家族基因功能提供一定參考。