雨生紅球藻鈣依賴蛋白激酶（CDPK）家族鑒定與表達分析

李亞男 張豪杰 梁夢靜 羅濤 李旺寧 張春輝 季春麗 李潤植薛金愛 崔紅利，2

（1.山西農業大學農學院 分子農業與生物能源研究所，太谷 030801；2.中國科學院煙臺海岸帶研究所 海岸帶生物學與生物資源利用重點實驗室，煙臺 264003）

鈣離子（Ca2+）是植物體內普遍存在的第二信使，通過瞬時改變植物生長發育過程中的Ca2+濃度，以及對外界環境刺激的反應，在信號轉導途徑中發揮重要作用［1-3］。在植物受到脅迫時，一些鈣感受器或鈣結合蛋白可以感知鈣離子濃度的變化并識別鈣信號，然后將其轉導到一系列下游反應，以應對環境脅迫損傷，提高自身抗性［4-5］。迄今已在植物細胞內發現多種鈣結合蛋白，如鈣調素（calmodulin,CaM）、類鈣調神經素B 亞基蛋白（calcineurin B-like proteins,CBLs）和鈣依賴蛋白激酶（calcium-dependent protein kinases,CDPKs 或CPKs）等［6］。

CDPK 是獨特的新型Ca2+傳感器。與其他已知的Ca2+傳感器相比，它具有一個催化功能的Ser/Thr激酶結構域和一個類鈣調素結合域（能結合鈣的EF-hand 基序）。因此，CDPK 同時具有Ca2+結合活性和激酶活性，能夠直接結合Ca2+，將鈣信號轉化為下游的磷酸化信號［7-8］。典型的CDPK 蛋白結構主要由N 端可變區、Ser/Thr 激酶域、自抑制域（連接域）以及類鈣調素結合域組成［6,9-10］。N 端可變區具有肉豆蔻?；妥貦磅；稽c，對CDPK 的亞細胞定位及調控植物生長發育十分重要［11］。自抑制結構域與激酶域相鄰，作為假底物與其結合，以使CDPK 失活［12］。目前已在萊茵衣藻（Chlamydomonas reinhardtii）［13］、擬南芥（Arabidopsis thaliana）［6］、煙 草（Nicotiana tabacum）［14］、水 稻（Oryza sativa）［15］、番茄（Solanum lycopersicum）［16］、玉米（Zea mays）［17］等中發現CDPK。

作為Ca2+的直接“傳感器”，CDPK 在調節植物生長和各種脅迫的應答反應中發揮重要作用。擬南芥Atcpk1 突變體表現出對鹽和干旱脅迫敏感，而其過表達株可顯著增強對鹽/干旱的抗性［18］。AtCDPK1 和AtCDPK2 基因在高鹽脅迫下快速上調表達，但在低溫下不受影響［19］。AtCDPK4/5/6/11 可以調節MAMP 或PAMP 觸發的免疫反應，并磷酸化轉錄因子WRKY 以調節下游基因表達［20］。過表達OsCDPK7 可通過上調一些脅迫響應基因的表達，進而增強植物對鹽堿或干旱的耐受性，但不參與冷脅迫響應。在萊茵衣藻中，大多數CDPK 基因在缺氮和鹽脅迫下顯著誘導，其中有6 個CrCDPK 基因在萊茵衣藻缺氮脅迫下的油脂積累中發揮積極作用［13］。大量研究表明，來自不同植物的CDPK 數量不同且具有功能多樣性。目前，關于藻類CDPK 基因家族鑒定及其生物學功能研究較少，有待進一步解析。

雨生紅球藻（Haematococcus pluvialis）是一種具有雙鞭毛的單細胞淡水綠藻，是目前已知天然蝦青素含量最高的生物，可在逆境條件下大量積累蝦青素［21］。蝦青素作為一種高價值的酮類胡蘿卜素，具有很強的抗氧化能力，已廣泛應用于食品、制藥和化妝品等領域［22-23］。溫度、光照和營養鹽等非生物脅迫可促進雨生紅球藻合成積累大量蝦青素［24］。干旱、鹽分和低溫等非生物脅迫都會引起植物細胞內Ca2+濃度的變化，進而誘導下游信號的傳遞途徑［25］。CDPK 可迅速感受瞬時Ca2+信號的變化，識別并磷酸化特異性底物后，引發各種生理反應，從而調控植物的生長發育及對多種脅迫的響應。而雨生紅球藻在非生物脅迫下，CDPK 的轉錄水平如何？具有怎樣的功能？是否參與蝦青素等類胡蘿卜素代謝的調控？回答這些問題有助于深入解析CDPK 基因參與雨生紅球藻脅迫應答和蝦青素的合成機制，以及建立優化的雨生紅球藻培養體系，高效生產蝦青素。

基于雨生紅球藻轉錄組數據，本研究擬對雨生紅球藻HpCDPK 基因家族進行鑒定和系統的生物信息學分析，進一步檢測UV-B、高白光、高藍光和α-酮戊二酸處理條件下HpCDPK 各成員的表達譜，以期為深入解析雨生紅球藻HpCDPK 基因的生物學功能，尤其是介導雨生紅球藻生長發育及蝦青素生物合成積累的分子機制提供新的科學依據。

1 材料與方法

1.1 材料

本試驗所用藻株為雨生紅球藻［Haematococcus pluvialis（Flotow,1844）］藻株,現存于山西農業大學分子農業與生物能源研究所。

1.2 方法

1.2.1 雨生紅球藻CDPK 基因家族成員鑒定及基本參數分析 雨生紅球藻氨基酸序列來自本地轉錄組數據庫。從Phytozome v13（https://phytozome-next.jgi.doe.gov/）下載擬南芥和萊茵衣藻CDPK 蛋白序列，與雨生紅球藻轉錄組中對應的蛋白序列進行同源比對，篩選雨生紅球藻中的CDPK 蛋白，其中同源比對E 值的閾值設置為1e-5。利用Interpro（http://www.ebi.ac.uk/interpro/）網站中 的search 命令檢測雨生紅球藻蛋白序列中同時包含Pkinase 結構域和EF-hand 基序的序列。對Interpro 序列比對的結果進行篩選。綜合以上結果，使用SMART 數據庫（http://smart.embl-heidelberg.de）對所含結構域進行驗證。利用在線軟件SWISS-MODEL（https://swissmodel.expasy.org/interactive）預測蛋白的三級結 構。利 用ExPASy Proteomics Server（http：//cn.expasy.org/tools）分析雨生紅球藻CDPK 蛋白序列的長度、等電點、分子質量等參數。肉豆蔻?；稽c由ExPASy 中 的Myristoylator 程 序預測（http://web.expasy.org/myristoylator/）。棕櫚?；稽c通過CSSPALM 4.0（http://csspalm.biocuckoo.org/）預 測。利用在線工具CELLO（http://cello.life.nctu.edu.tw/）和Softberry（http://linux1.softberry.com/berry.phtml）預測雨生紅球藻HpCDPK 的亞細胞定位。

1.2.2 雨生紅球藻CDPK 基因系統進化分析 利用MEGA7 軟件中的 Clustal W 對鑒定所得雨生紅球藻CDPK 蛋白序列與擬南芥、萊茵衣藻的CDPK 蛋白序列進行比對?；谛蛄斜葘Y果，采用鄰近法構建系統進化樹，其中bootstrap 值設定為1 000。

1.2.3 雨生紅球藻CDPK 蛋白的保守基序和功能結構域分析 利用在線網站MEME Suite v 5.5.0（https://meme-suite.org）預測雨生紅球藻CDPK 基因家族成員蛋白序列中的保守基序（motif），其中motif 的最大數量設為15，最優的motif 寬度為10-100 個氨基酸，其余參數為默認數值。利用SMART 數據庫（http://smart.emblheidelberg.de/）分析各CDPK 蛋白中的功能結構域。并利用TBtools（https://github.com/CJ-Chen/TBtools）將雨生紅球藻 CDPK 基因的系統進化樹和保守基序數據進行合并。

1.2.4 雨生紅球藻CDPK 基因表達特性分析 高藍光和高白光處理組：所用培養基為BG11，在光強25 μmol/（m2·s）（LED 白光）、溫度（22.5±1）℃、光/暗周期（12 h/12 h）的條件下靜置培養，每8 h 搖1 次。將對數期雨生紅球藻在黑暗環境下處理24 h，經離心（5 000 r/min，4℃，5 min）收集，用新鮮培養基重懸細胞，分別于正常培養（N,25 μmol/（m2·s））、高白光（W,500 μmol/（m2·s））和高藍光（B,500 μmol/（m2·s））條件下處理72 h。最后離心收藻，將樣品于液氮速凍后存于-80℃，每個處理3 個重復，并進行轉錄組測序［26］。

UV-B 處理組：采用 MCM 培養基培養藻細胞，培養條件為：光強為25 μmol/（m2·s）、溫度 25℃、光暗比為12 h/12 h，靜置培養，每天手搖2-3 次。使用 Philips TL20W/01RS 窄帶 UV-B 燈管（15 W）提供UV-B 光源，通過移動錐形瓶與紫外燈的距離調節強度，光強使用勒克斯計測量，在試驗前連續照射24 h 使系統穩定［27］。離心收集正常光照組（L,25 μmol/（m2·s））、低UV-B 輻射組（LU200,200 W）、高UV-B 輻射組（LU400,400 W）以及對應的黑暗條件下處理組（D）、低UV-B 輻射組（DU200,200 W）、高UV-B 輻射組（DU400,400 W）培養的雨生紅球藻藻株經液氮速凍后存于-80℃，每個樣本3 個重復，用于轉錄組測序分析。

α-酮戊二酸（α-ketoglutaric acid,OG）處理組：將培養至對數生長期的雨生紅球藻以5 000 r/min 離心5 min，收集藻細胞，然后重新懸浮在正常培養基和缺氮培養基中，缺氮（ND）培養基為不添加硝酸鈉的BG11 培養基。對照為高光處理（CK），光照強度為270 μmol/（m2·s）。OG 濃度設置為10.0 mg/L。離心收集CK、OG、ND 和ND-OG 組培養的雨生紅球藻藻株經液氮速凍后存于-80℃，每個樣本3 個重復，用于轉錄組測序分析。

1.2.5 雨生紅球藻HpCDPK 基因與類胡蘿卜素合成基因的共表達分析 為預測HpCDPK2、HpCDPK3和HpCDPK5 的功能，使用從轉錄組中鑒定的類胡蘿卜素合成相關基因，包括AACT、HMGS、MVD、IDI、DXS、DXR、CMT、CMK、MDS、HDS、HDR、GGPS、PSY、PDS、Z-ISO、ZDS、CRTISO、LYCB、LYCE、BKT、BCH、CYP97A、CYP97B、CYP97C、ZEP、VDE 和NCED/CCD。結合熱圖分析中使用的轉錄組數據，使用GraphPad Prism 9.0 中的ggcor 函數和默認的Pearson 相關系數計算HpCDPK 與類胡蘿卜素合成相關基因表達的相關性。

2 結果

2.1 雨生紅球藻CDPK基因及編碼蛋白的基本信息

以擬南芥和萊茵衣藻CDPK 蛋白序列作為查詢序列，通過BlastP 程序從雨生紅球藻轉錄組中檢索到候選CDPK 蛋白序列。通過SMART 和序列比對驗證后，總共有7個非冗余的HpCDPK基因被鑒定，并命名為HpCDPK1-HpCDPK7（表1）。這些HpCDPK基因編碼氨基酸長度主要集中在478-1 479 aa，所編碼蛋白質分子質量大小在53 280.35-153 784.74 Da之間，等電點變化區間為5.71-7.29，其中除HpCDPK6 蛋白等電點大于7 顯堿性，其余6 個蛋白等電點小于7 顯酸性。三級結構預測結果表明，HpCDPK 蛋白主要由α 螺旋組成（圖1）。在鑒定的7 個HpCDPK 中僅HpCDPK7 預測到在N 末端有肉豆蔻化基序。利用CELLO（http://cello.life.nctu.edu.tw/）和Softberry（http://linux1.softberry.com/berry.phtml）預測HpCDPK 蛋白亞細胞定位發現，HpCDPK1、3、4 均定位于細胞質，該定位更易識別Ca2+內流信號。

圖1 雨生紅球藻CDPK 蛋白三級結構預測Fig.1 Prediction of tertiary structures of CDPK proteins in H.pluvialis

表1 雨生紅球藻CDPK 基因家族成員信息Table 1 Members of CDPK gene family in H.pluvialis

2.2 雨生紅球藻CDPK基因家族系統進化分析

為了剖析CDPK 基因家族成員的進化關系，我們使用雙子葉模式植物擬南芥、單細胞模式生物萊茵衣藻和雨生紅球藻CDPK 基因家族成員構建了系統發育樹，這56 個成員可分為4 組（圖2）。第1組成員最多，包含了31 個AtCDPK，皆為擬南芥成員。第2 組包括4 個CrCDPK 和2 個HpCDPK。第3組是由AtCDPK18、AtCDPK16和AtCDPK28聚集而成。第4 組則是由11 個CrCDPK 和5 個HpCDPK 組成。此外，雨生紅球藻CDPK 基因家族成員與高等植物擬南芥的AtCDPK 基因家族成員進化關系相對較遠，與同為低等植物的萊茵衣藻CrCDPK 基因家族成員親緣關系較近。

圖2 雨生紅球藻、萊茵衣藻和擬南芥CDPK 基因系統發育樹Fig.2 Phylogenetic tree of CDPK genes in H.pluvialis,C.reinhardtii and A.thaliana

2.3 雨生紅球藻CDPK蛋白特征分析

保守基序（motif）是同源蛋白中的潛在功能區。為了研究HpCDPK 的蛋白結構和功能多樣性，我們使用全長的HpCDPK 蛋白序列構建了一個單獨的無根系統發育樹，并將其與相應基因序列的保守基序和結構域組合進行了比較分析（表1 和圖3）。利用在線MEME 工具分析獲得HpCDPK 蛋白中最保守的9 個基序（圖4），以此幫助我們研究它們的功能。無根系統發育樹表明，HpCDPK 大致分為3 個主要亞組。motif7、motif2、motif3、motif1、motif4 為蛋白激酶結構域，該區域序列在不同物種間的同源性較高；motif5、motif8 和motif6 為EF-hand 基序，其可分為N 端和C 端兩部分，且C 端的鈣親和力高于N 端（圖3，5 和6）?；蚍治鼋Y果表明，每個亞群內部保守的基序結構支持它們之間密切的進化關系，但不同亞群之間可能存在功能上的差異。此外，使用Interpro 和SMART 蛋白質數據庫和序列比對預測了HpCDPK 的蛋白質結構（圖3-6）。結果表明，本研究鑒定的所有CDPK 都具有4 個特征結構，即一個可變的N-末端結構域，一個催化Ser/Thr 蛋白激酶結構域，一個自抑制區和一個類鈣調素結合域。此外，對HpCDPK 的保守結構域和重要位點分析發現，7 個HpCDPK 蛋白中，僅有HpCDPK1 具有C2 保守結構域。在激酶保守結構區域，所有HpCDPK 具有Ser/Thr 活性位點；除HpCDPK5 外，其余蛋白都含有激酶ATP 結合位點。此外，在EF 手型保守區域，7 個HpCDPK 都含有Ca2+結合位點，但各蛋白所含Ca2+結合位點的數量不同（圖7）。

圖3 雨生紅球藻HpCDPK 家族成員的保守基序和功能結構域Fig.3 Conserved motif and functional domain of HpCDPK family members of H.pluvialis

圖4 MEME 預測的雨生紅球藻HpCDPK 保守基序Fig.4 Conserved motifs of HpCDPKs in H.pluvialis predicted by MEME

圖5 雨生紅球藻HpCDPK 蛋白保守PDK 結構域的多序列比對Fig.5 Multiple sequence alignments of conserved PDK domains of HpCDPK proteins in H.pluvialis

圖6 雨生紅球藻HpCDPK 蛋白保守EF 基序的多序列比對Fig.6 Multiple sequence alignments of conserved EF motifs of HpCDPK proteins in H.pluvialis

圖7 雨生紅球藻HpCDPK 的保守結構域和重要位點分析Fig.7 Conserved domain and important site analysis of HpCDPKs in H.pluvialis

2.4 雨生紅球藻CDPK基因在非生物脅迫下的表達模式

為研究HpCDPK 基因是否在雨生紅球藻脅迫響應過程中發揮功能，我們檢測3 種不同脅迫下雨生紅球藻HpCDPK 基因的表達模式。在正常培養條件（N）下（圖8-A），雨生紅球藻HpCDPK 基因均表達，且HpCDPK1 表達量最高。在高白光和高藍光處理下，HpCDPK1 和HpCDPK6 表達量均顯著降低，而且HpCDPK1 表達量減少幅度遠低于HpCDPK6。

圖8 雨生紅球藻HpCDPK 基因在非生物脅迫下的表達模式Fig.8 Expression patterns of HpCDPK genes in H.pluvialis in response to abiotic stresses

在正常光照條件下（圖8-B），UV-B（200 W 和400 W）處理的雨生紅球藻HpCDPK2 和HpCDPK5基因的表達量均提高；HpCDPK3 在UV-B（200 W）輻射時上調，當UV-B 輻射增至400 W 時HpCDPK3的表達量下調。在黑暗條件進行UV-B（200 W和400 W）處理下，HpCDPK3 的表達 量下調，HpCDPK5 的表達量均上調；HpCDPK2 在UV-B（200 W）和UV-B（400 W）輻射時表達量均上調。

在高光條件（CK）下（圖8-C），HpCDPK3 的基因表達量最高，HpCDPK5 的表達量最低。與CK相比，在OG、ND 和ND-OG 條件下，HpCDPK6 的表達量均下調；HpCDPK1、HpCDPK2、HpCDPK4、HpCDPK5 和HpCDPK7 的表達量均差異性上調，且在ND 和ND-OG 條件下的基因表達量高于OG 條件下的基因表達量，但其基因表達量在ND 和ND-OG條件下差異不顯著。此外，在ND 和ND-OG 條件下，HpCDPK3、HpCDPK5 和HpCDPK6 的表達量差異同樣不明顯。

2.5 HpCDPK基因與類胡蘿卜素及蝦青素合成基因表達量的相關性分析

有研究表明，在UV-B 脅迫條件下可顯著誘導雨生紅球藻蝦青素積累。為了探索HpCDPK 是否參與蝦青素合成積累的調控，我們根據課題組已有的UV-B 誘導條件下雨生紅球藻轉錄組數據，對不同強度UV-B 誘導處理下的雨生紅球藻中HpCDPK2、HpCDPK3 和HpCDPK5 與蝦青素合成及積累的關鍵基因轉錄水平進行相關性分析。結果表明，HpCDPK2、HpCDPK3 和HpCDPK5 與蝦青素合成通路基因表達高度相關，暗示HpCDPK 可能介導雨生紅球藻蝦青素生物合成（圖9）。進一步相關性分析（圖10）顯示，HpCDPK2 與蝦青素合成關鍵基因PSY、BCH 和BKT 和類胡蘿卜素合成關鍵基因Z-ISO（圖10-A、C、D 和B）轉錄水平呈正相關（P＜0.05）。HpCDPK3 與類胡蘿卜素合成基因關鍵基因ZDS、CYP97B 和HMGS（圖10-E、F 和H）轉錄水平呈正相關（P＜0.05）；而與其降解基因NCED/CCD（圖10-G）轉錄水平呈負相關（P＜0.05）。HpCDPK5 與類胡蘿卜素合成關鍵基因DXS、DXR 和LYCB 以及其降解基因NCED/CCD（圖10-I、J、L 和K）轉錄水平呈正相關（P＜0.05）。由此推測HpCDPK 可能在UV-B 條件下參與蝦青素生物合成的調控，亦可能在類胡蘿卜素的積累中發揮類似功能，具體分子機制還需進一步試驗驗證。

圖9 HpCDPK 基因與蝦青素合成及積累基因的表達量相關性分析Fig.9 Correlation between HpCDPKs and gene expressions for astaxanthin synthesis and accumulation

圖10 HpCDPK2、HpCDPK3 和HpCDPK5 基因與類胡蘿卜素及蝦青素合成及積累基因的表達量相關性Fig.10 Correlation between HpCDPK2,HpCDPK3,HpCDPK5 genes and gene expressions for the biosynthesis and accumulation of carotenoid and astaxanthin

3 討論

CDPK 激酶在植物生長、發育和脅迫反應中發揮重要作用［10］。目前已在多種植物中闡明其生物學功能。雨生紅球藻作為高效積累天然蝦青素的優異種質，其CDPK 成員能否參與脅迫條件誘導蝦青素合成的調控尚不清楚。本研究共鑒定7 個HpCDPK家族成員。通過與擬南芥和萊茵衣藻CDPK 成員比對與進化分析，分為4 組。其中HpCDPK 與來自萊茵衣藻的CDPK 聚為一類，與高等植物（擬南芥）來源的AtCDPK 明顯分開，形成兩個姐妹分枝，暗示在進化過程中具有種屬特異性［13］。

N-肉豆蔻?；揎椯x予蛋白靈活可變的膜結合能力，對蛋白的膜靶向尤為重要［28-29］。CDPK 具有N-肉豆蔻?；?，暗示定位在質膜，在胞膜結合和鈣信號轉導中發揮作用［1,30-32］。在本研究中，僅有HpCDPK7 預測含有N-肉豆蔻?；?，但亞細胞定位預測顯示在線粒體。研究表明玉米ZmCK1 含有 N 末端肉豆蔻?；?，但ZmCK1:hGFP 融合蛋白定位于細胞質和細胞核［33］。CDPK 蛋白的亞細胞定位還可能受到其他因素影響，需要進一步的實驗驗證。

蛋白激酶結構域能將核苷三磷酸（ATP）的γ-磷酸轉移到蛋白質側鏈，導致構象變化進而影響其功能［34］。因此，含有ATP 結合位點的HpCDPK 成員可以利用ATP 作為磷酸基團的來源，以此來磷酸化靶底物［35］。盡管蛋白激酶結構域非常保守，但不同成員間的ATP 結合位點和Ser/Thr 蛋白激酶活性位點序列存在著大量的變異。例如，HpCDPK3 的激酶ATP 結合位點序列與其他成員此位點序列相比較短，HpCDPK5 不具有激酶ATP 結合位點序列。這可能是因為其保守性較低導致的，也可能是不同成員功能差異性所在。

為探究HpCDPK 是否參與雨生紅球藻脅迫響應，我們檢測了3 種不同脅迫下雨生紅球藻HpCDPK 基因的表達模式。結果顯示，HpCDPK1 和HpCDPK6基因表達受到光條件的影響，但HpCDPK1 對藍光更敏感。UV-B 輻射對植物的生長發育具有重要的作用，高強度的UV-B 作為逆境脅迫因子，影響植物正常的生理功能；而低強度的UV-B 作為信號調控因子，促進植物光形態的建成［36-39］。已有研究表明，UV-B輻射顯著影響長春花懸浮細胞培養中CDPK 蛋白激酶的活性［40］。HpCDPK5 無論在黑暗條件還是光照條件下UV-B（400 W）處理，其基因的表達量均上調，而HpCDPK3 表達則是下調。因此，HpCDPK3和HpCDPK5 在高強度UV-B（400 W）誘導下，其基因表達量具有相反模式。與之相似，OsCDPK2 的表達也與光誘導密切相關。暴露在光下的水稻綠葉中幾乎檢測不到OsCDPK2 蛋白，但當植株轉移到黑暗中時，該蛋白水平急劇上升［41］。雨生紅球藻在高光（HL）和缺氮（ND）雙重脅迫，以及添加α-酮戊二酸（OG）處理的HpCDPK 基因表達模式結果顯示，ND 和ND-OG 處理激活了更多差異表達的HpCDPK 基因，這與萊茵衣藻CDPK 在缺氮處理下的表達結果相一致［13］。其中HpCDPK7 基因轉錄水平在缺氮條件下上調最為明顯。此外，結果還表明，HpCDPK 基因表達更易受到ND 影響，而OG 處理對其影響較小。

類胡蘿卜素合成的前體物質是異戊烯焦磷酸（IPP），在雨生紅球藻中，IPP 被認為只能由MEP途徑合成［42］。已有研究表明，UV-B 輻射可提高雨生紅球藻蝦青素含量［43］。因此，我們對不同強度UV-B 輻射處理下的雨生紅球藻中HpCDPK2、HpCDPK3 和HpCDPK5 與類胡蘿卜素合成相關基因轉錄水平進行相關性分析，結果顯示，HpCDPK5 和參與IPP 合成的酶基因DXS、DXR 和HDS 的轉錄水平呈正相關。IPP 在酶的作用下最終生成香葉基香葉基焦磷酸（GGPP），它是后續胡蘿卜素合成的直接前提物質。GGPP 在經PSY、PDS、Z-ISO、ZDS 等酶的催化反應后生成全反式番茄紅素，后者在兩種不同類型的環化酶（LYCB 和LYCE）協同作用下形成不同的胡蘿卜素［44-45］。其中，番茄紅素在LYCB的催化下生成β-胡蘿卜素，然后在BKT 和BCH 酶的參與反應下最終生成蝦青素。結果顯示，在GGPP合成全反式番茄紅素階段，HpCDPK2 與PSY 和Z-ISO基因的轉錄水平呈正相關，HpCDPK3 與ZDS 基因的轉錄水平呈正相關。在合成蝦青素階段，HpCDPK5與LYCB 基因的轉錄水平呈正相關，HpCDPK2 與蝦青素合成關鍵基因BKT 和BCH 轉錄水平正相關。此外，我們發現NCED/CCD 與HpCDPK3 基因的轉錄水平負相關，與HpCDPK5 的轉錄水平呈正相關。類胡蘿卜素裂解雙加氧酶（CCD）家族包括CDD 和NCED，是雙加氧酶酶促降解途徑中的重要酶，可特異性識別不同底物和不同的雙鍵位置進行斷裂，從而形成自然界中多種多樣的脫輔基類胡蘿卜素［44］。因此，HpCDPK3 和HpCDPK5 基因的可能發揮相反的功能。

已有研究表明，CDPK 可通過磷酸化調節轉錄因子、離子通道、轉運蛋白或一些其他信號通路相關的蛋白活性，從而調控基因的表達和維持ROS穩 態［46］。在本研究中，HpCDPK2、HpCDPK3 和HpCDPK5 基因的表達模式與類胡蘿卜素合成相關基因的表達具有相關性。因此，我們推測上述基因可能通過磷酸化調節某種底物蛋白或轉錄因子的活性，進而調節雨生紅球藻類胡蘿卜素合成相關基因的表達，從而增加其含量。

4 結論

本研究在雨生紅球藻轉錄組數據庫中鑒定到7個HpCDPK 家族成員，HpCDPK 基因表達量受到多種脅迫條件的影響，其中HpCDPK7 基因轉錄水平在缺氮條件下上調最為明顯。此外，HpCDPK2 與多個控制類胡蘿卜素和蝦青素合成的酶基因轉錄水平主要呈正相關，預示著HpCDPK2 可能在類胡蘿卜素及蝦青素合成積累中發揮積極調控作用。