基于丹參基因組的蛋白磷酸酶2C家族的系統分析△

徐志超，浦香東，宋經元

(中國醫學科學院 北京協和醫學院 藥用植物研究所/國家中醫藥管理局 中藥資源保護重點研究室，北京 100193)

蛋白激酶和蛋白磷酸酶所介導的可逆磷酸化反應在植物信號傳遞和逆境脅迫響應進程中發揮重要調控作用。蛋白磷酸酶按照底物特異性分為絲氨酸/蘇氨酸蛋白磷酸酶類和酪氨酸蛋白磷酸酶類，PP2蛋白磷酸酶作為絲氨酸/蘇氨酸蛋白磷酸酶類的代表，根據亞基結構等分為PP2A、PP2B和PP2C 3類[1]。擬南芥Arabidopsisthaliana的PP2CA亞家族已證實負調控ABA信號通路?；瘜W遺傳學和蛋白互作研究揭示PYR/PYL/RCAR復合體為ABA的受體，當ABA信號缺失時，PP2C通過去磷酸化使SnRK2激酶失去活性；而響應ABA信號時，受體PYR/PYL/RCAR與ABA結合后抑制PP2C的活性，SnRK2激酶活性釋放并磷酸化下游因子實現信號的傳遞及細胞對ABA脅迫的響應[2-3]。除了ABA信號途徑，PP2C還參與多種植物信號轉導途徑的調控，如鉀通道蛋白、MAPK創傷信號途徑等。擬南芥和水稻基因組分別注釋到80和78個PP2C家族成員，劃分為13和11個亞家族，負調控的PP2CA亞家族基因響應不同的脅迫刺激，如ABA脅迫響應，而PP2CD亞家族可能正調控ABA介導的信號通路[4]?；诨蚪M系統鑒定植物PP2C基因家族，有助于揭示其在植物逆境脅迫下的響應機制。

次生代謝產物是藥用植物主要的藥效成分，同時次生代謝產物的產生和積累也是植物對環境脅迫響應的一種適應。環境脅迫影響藥用植物的品質，因此藥用植物逆境脅迫的分子機制解析將為中藥材種植環境及栽培技術的規范提供依據，推動中藥材可持續發展。研究證實逆境脅迫雖抑制植物生長，但促進藥用植物次生代謝產物的合成和積累，適當的逆境脅迫提升藥用植物的品質[5-6]。如適度的干旱脅迫誘導甘草酸合成途徑關鍵酶基因的上調，并促進甘草酸和甘草苷的積累[7]；輕度干旱可顯著提高丹參活性成分丹參酮類和丹酚酸類化合物的合成與積累，同時增強植株對磷、鉀、鈣、等礦質元素的吸收與積累[8]。

丹參為常用大宗藥材，遺傳多樣性豐富，是藥用植物次生代謝研究理想的模式植物[9-12]。不同產地丹參品質參差不齊，研究逆境脅迫影響丹參品質的分子機制，將有助于丹參優良品種的種植農藝措施的改良。本研究基于丹參基因組注釋PP2C基因家族，揭示PP2C在丹參不同組織部位的差異表達及對茉莉酸甲酯MeJA信號的響應，為丹參逆境脅迫分子機制的研究奠定基礎。

1 材料和方法

1.1 丹參基因組、轉錄組數據

丹參基因組組裝、注釋及轉錄組測序工作由本課題組完成并已釋放[10]。丹參基因組測序原始數據(Roche 454、PacBio RS、Illumina平臺)以及不同組織器官及處理的轉錄組數據(Illumina HiSeq2500測序平臺)已上傳至NCBI的SRA數據庫，登錄號：SRX753381、SRR1640458、SRP028388、SRP051564[10-11，13]。

1.2 基于蛋白序列比對的丹參PP2C基因家族成員鑒定

下載擬南芥PP2C基因家族A-L亞家族的蛋白，采用BLASTP比對丹參基因組蛋白序列，Evalue設置為1E-5。從基因組數據庫中提取所有PP2C序列，根據NCBI 網站BLASTP比對的結果，利用APOLLO軟件校正PP2C成員的CDS、DNA、蛋白序列信息。在蛋白質數據庫PROSITE(http：//prosite.expasy.org/)中分析所有蛋白的保守結構域序列。

1.3 PP2C家族成員系統發育、基因結構和保守基序分析

將丹參和擬南芥PP2C家族成員的全長蛋白序列采用MUSCLE做多序列比對，構建NJ(Neighbour-Joining)樹，選擇Jones-Taylor-Thornton(JTT)模型，重復次數1000次?；贑DS序列及對應基因組序列，利用在線服務器Gene Structure Display Server(GSDS 2.0)(http：//gsds.cbi.pku.edu.cn/index.php)對基因結構進行分析。利用MEME(Suite version 4.9.1)(http：//meme-suite.org/)分析PP2C基因家族成員的保守基序，具體參數：E值小于2×10-30，基序重復次數不限，氨基酸個數10～200個。提取PP2C家族成員啟動子區1500 bp通過NEW PLACE在線軟件(https：//sogo.dna.affrc.go.jp/cgi-bin/)分析順式作用元件。

1.4 PP2C家族基因的差異表達分析

采用已報道[10-11，13]的丹參根組織(周皮、韌皮部、木質部)、根、莖、葉、花以及MeJA處理的轉錄組數據，通過HISAT2比對到丹參基因組，通過cufflinks計算丹參編碼基因的FPKM(Fragments Per Kilobase of exon model per Million mapped reads)值，篩選出PP2C基因的表達量。FPKM值通過Z-score(Standard score)標準化處理后，采用R語言的gplots下的heatmap來繪制熱圖。采用皮爾森相關系數(Pearson’s Correlation Tests，r)分析基因在不同組織部位的共表達情況。

2 結果和分析

2.1 丹參PP2C基因家族的注釋及系統進化分析

在丹參基因組中共預測到83個PP2C基因成員，基因數量與高等植物中已報道的PP2C基因數量近乎一致，如擬南芥(80)[4]、水稻Oryzasativa(78)[4]、二穗短柄草Brachypodiumdistachyon(86)[14]、谷子Setariaitalica(80)[15]。丹參PP2C的長度范圍在145 aa(SMil_00028881)到1048 aa(SMil_00011829)，覆蓋基因組的長度在698 bp(SMil_00028881)到8891 kb(SMil_00011829)。

用丹參和擬南芥PP2C的氨基酸序列構建系統發育樹，根據擬南芥的亞家族分類，丹參PP2C基因家族被分為13個亞家族A～L(A、B、C、D、E、F1、F2、G、H、I、J、K、L)，各亞家族成員個數分別為11、7、6、8、9、6、4、9、5、2、1、7、2個，其中PP2CA亞家族成員數量最多，如圖1所示。除此之外，6個成員聚為其他獨立的分枝，暫未命名。丹參PP2C基因家族成員兩兩比對計算同源關系，發現PP2CF1亞家族的SMil_00030307和SMil_00028881的同源關系最近為98.57%，而SMil_00011829(PP2CL)與SMil_00021239(PP2CC)的同源關系最低為16.46%。

2.2 丹參PP2C基因家族的基因結構分析

PP2C基因家族的基因結構顯示亞家族內和不同亞家族間基因結構差異較大，內含子數量變異較大，其中PP2CL亞家族的SMil_00011829內含子數量最多，為15個。丹參PP2CA亞家族的11個成員，命名為SmPP2CA1-SmPP2CA11，內含子數量為2～7個，多數為3個內含子，基因結構相似，如圖2所示。

PP2CA亞家族各基因氨基酸大小在286～541 aa，等電點在4.78～6.84，分子量大小31～58 kD。對PP2CA亞家族的11個基序motif進行預測，Motif 1在11個亞家族成員中保守存在，大部分PP2CA成員保守存在motif 2、4、5、6、7。PP2A亞家族的系統進化分析顯示3個同源分枝，其中同源分枝2和3的兩對亞家族成員的基因結構和保守結構域基本相同。

2.3 PP2C基因的差異表達分析

丹參酮是丹參重要的脂溶性活性成分，丹參根及根周皮是丹參酮的合成和積累部位[11]。本研究對丹參不同組織部位，如根、莖、葉、花、根周皮、根韌皮部和根木質部的轉錄組數據比對到丹參基因組，計算FPKM值，將FPKM值用scale函數做Z-score標準化處理并聚類，如圖3所示。

圖1 丹參PP2C基因家族系統進化分析

差異表達結果顯示28.91%(24)的基因低表達(FPKM<10)，其中7個基因顯示沉默表達(FPKM<1)。丹參根3個組織部位(周皮、韌皮部、木質部)的PP2C基因表達聚在一起。組織特異性表達分析發現：PP2CE亞家族的SMil_00028462和SMil_00013051兩個基因在花中特異表達，可能與花發育的調控相關；PP2CE亞家族的SMil_00004715在丹參根周皮特異高表達，在莖、葉、花等器官表達極低，推測間接調控丹參酮的生物合成。3個PP2CA亞家族基因在丹參組織器官中表達低(FPKM<10)，PP2CA亞家族的同源分枝3的兩個基因在丹參不同組織部位中表達最高，且在根、莖、葉、花中表達差異不顯著，呈共表達(r=0.81)。茉莉酸甲酯MeJA作為重要的信號轉導因子能夠誘導植物次生代謝產物的產生，在丹參中MeJA能夠正調控丹參酮的生物合成。對于丹參葉片MeJA處理12 h后的轉錄組數據，分析丹參PP2C家族成員編碼基因的表達，發現13個PP2C基因在MeJA處理下顯著降低，其中包含6個PP2CA和3個PP2CK亞家族的基因。PP2CA亞家族同源分枝2的兩個基因在丹參所有轉錄組(不同組織器官及MeJA處理)的表達中呈顯著共表達(r=0.95)。

分析丹參PP2CA基因的啟動子中的順式作用元件發現：PP2CA的啟動子中ABA響應元件(ABRE，ACGTGG)均有分布(見表1)，大部分啟動子區呈密集分布，與響應ABA信號密切相關；在PP2CA的啟動子中富含WRKY、MYB和MYC 3類轉錄因子的作用元件(見表1)，推測在響應逆境脅迫時受轉錄因子的調控；干旱和低溫響應的元件部分分布。

表1 丹參PP2CA亞家族成員啟動子區順式作用元件分析

注：ABRE、DRE、LTRE代表響應ABA信號、干旱、低溫的作用元件；WRKY、MYB、MYC、GCC代表轉錄因子WRKY、MYB、MYC、AP2的結合元件。

3 討論

本研究基于丹參全基因組測序，系統鑒定蛋白磷酸激酶PP2C家族成員，丹參基因組注釋的PP2C基因數量(83)與已報道的植物中相似；丹參與擬南芥PP2C序列的系統進化顯示丹參PP2C亞組的分類與擬南芥相似，亞家族成員數量也大多相似，保守結構域相似度高，證明PP2C基因家族及其基因結構在植物中進化較為保守。低等維管植物(石松綱)卷柏Selaginellatamariscina[16]中的PP2C基因家族鑒定到46個成員(未發表)，遠少于已報道的高等植物，PP2C在進化過程中應對復雜逆境環境時顯著的基因擴張反映PP2C基因家族在植物逆境響應進程中重要的調控作用。

PP2CA亞家族成員作為ABA的受體之一，負調控ABA信號轉導途徑，響應多種逆境脅迫，如干旱、低溫等[17-19]。丹參PP2CA基因啟動子含有多種逆境脅迫應答的順式作用元件如ABRE、DRE、LTRE等，大多數的PP2CA基因響應MeJA信號轉導因子。PP2CA的同源分枝2、3的基因結構相似，motif高度保守，且顯著共表達，推測為基因組復制導致的相似或相同功能基因的分化。有報道指出植物中ABA信號與花發育相關，如ABA信號途徑中一個正調控因子ABI4(乙烯響應的轉錄因子)正調控開花關鍵基因FLC的轉錄，進而負調控開花時間[20]。丹參PP2CE亞家族的兩個基因在花中特異表達，一個基因在根中特異高表達，預測PP2CE亞家族與植物組織器官發育相關，如根和花的發育，但需要進一步的實驗驗證。

丹參是最常用的大宗藥材之一，生態適應性較強，種植面積廣泛分布在華東、華北、西北、西南、中南等地區，野生資源大幅減少，主要來源于人工栽培品，然而不同產地的丹參品質差異明顯。研究指出輕度干旱可顯著提高丹參活性成分丹參酮類和丹酚酸類化合物的合成與積累，同時增強植株對礦質元素的吸收，因此丹參對不同環境的適應導致活性成分合成和積累的差異，揭示丹參對逆境脅迫的分子機制有助于指導丹參優良品種的栽培與繁育。丹參PP2C基因家族的系統分析將為丹參響應逆境脅迫分子機制的解析奠定基礎，同時丹參作為藥用模式植物，將對其他藥用植物的栽培育種及道地性的理解具有指導意義。

[1] Shi Y.Serine/Threonine Phosphatases：Mechanism through Structure[J].Cell，2009，139(3)：468-484.

[2] Zhu J K.Abiotic Stress Signaling and Responses in Plants[J].Cell，2016，167(2)：313.

[3] Danquah A，Zelicourt A D，Colcombet J，et al.The role of ABA and MAPK signaling pathways in plant abiotic stress responses[J].Biotechnol Adv，2014，32(1)：40.

[4] Xue T，Wang D，Zhang S，et al.Genome-wide and expression analysis of protein phosphatase 2C in rice and Arabidopsis[J].BMC Genomics，2008，9(1)：550.

[5] 司燦，張君毅，徐護朝.藥用植物在干旱脅迫下生長代謝變化規律及應答機制的研究進展[J].中國中藥雜志，2014，39(13)：2432-2437.

[6] 楊枝中，馬逾英，盧曉琳，等.干旱和溫度脅迫在藥用植物研究中的應用現狀與展望[J].時珍國醫國藥，2010，21(5)：1223-1225.

[7] 張春榮，桑雪雨，渠萌，等.基于轉錄組測序揭示適度干旱脅迫對甘草根基因表達的調控[J].中國中藥雜志，2015，40(24)：4817-4823.

[8] 劉大會，郭蘭萍，黃璐琦，等.土壤水分含量對丹參幼苗生長及有效成分的影響[J].中國中藥雜志，2011，36(3)：321-325.

[9] 宋經元，羅紅梅，李春芳，等.丹參藥用模式植物研究探討[J].藥學學報，2013(7)：1099-1106.

[10] Xu H，Song J，Luo H，et al.Analysis of the genome sequence of the medicinal plantSalviamiltiorrhiza[J].Mol Plant，2016，9(6)：949.

[11] Xu Z，Peters R J，Weirather J，et al.Full-length transcriptome sequences and splice variants obtained by a combination of sequencing platforms applied to different root tissues ofSalviamiltiorrhizaand tanshinone biosynthesis[J].Plant J，2015，82(6)：951.

[12] Xu Z，Song J.The 2-oxoglutarate-dependent dioxygenase superfamily participates in tanshinone production inSalviamiltiorrhiza[J].J Exp Bot，2017，68(9)：2299-2308.

[13] Luo H，Zhu Y，Song J，et al.Transcriptional data mining ofSalviamiltiorrhizain response to methyl jasmonate to examine the mechanism of bioactive compound biosynthesis and regulation[J].Physiol Plantarum，2014，152(2)：241.

[14] Cao J，Jiang M，Li P，et al.Genome-wide identification and evolutionary analyses of the PP2C gene family with their expression profiling in response to multiple stresses inBrachypodiumdistachyon[J].BMC Genomics，2016，17(1)：175.

[15] 閔東紅，薛飛洋，馬亞男，等.谷子PP2C基因家族的特性[J].作物學報，2013，39(12)：2135-2144.

[16] Banks J A.Selaginella and 400 million years of separation.[J].Annu Rev Plant Biol，2009，60(60)：223.

[17] Zelicourt A D，Colcombet J，Hirt H.The Role of MAPK Modules and ABA during Abiotic Stress Signaling[J].Trends Plant Sci，2016，21(8)：677.

[18] Zhu J K.Salt and Drought Stress Signal Transduction in Plants[J].Annu Rev Plant Biol，2002，53(53)：247.

[19] Schweighofer A，Hirt H，Meskiene I.Plant PP2C phosphatases：emerging functions in stress signaling.[J].Trends Plant Sci，2004，9(5)：236-243.

[20] Shu K，Chen Q，Wu Y，et al.ABSCISIC ACID-INSENSITIVE 4 negatively regulates flowering through directly promoting Arabidopsis FLOWERING LOCUS C transcription[J].J Exp Bot，2016，67(1)：195.