澳洲坚果MiSTPP3和MiSTPP7基因克隆与表达分析

蔡元保 杨祥燕 庞新华 李季东 林玉虹

摘  要：以澳洲坚果品种‘O.C.’为试材，利用RT-PCR技术克隆获得2个丝氨酸/苏氨酸磷酸酶基因，分别命名为MiSTPP3和MiSTPP7，其编码序列长度分别为2095 bp和1700 bp。序列分析和系统进化分析表明，MiSTPP3含有MPP_PP5_C结构域，属于PP5亚家族;MiSTPP7含有MPP_PP7结构域，属于PP7亚家族。蛋白理化性质和亚细胞定位分析表明，MiSTPP3是一个稳定的亲水性蛋白，可能定位于质膜;MiSTPP7是一个不稳定的亲水性蛋白，可能定位于线粒体基质。生物信息学分析表明，MiSTPP3和MiSTPP7都是非分泌跨膜蛋白，其二级和三级结构的主要元件都是α螺旋、无规则卷曲和延伸链。RT-qPCR分析表明，MiSTPP3和MiSTPP7在澳洲坚果根、茎、叶、花和小果中都有不同的表达量;MiSTPP3受低温胁迫下调表达，而MiSTPP7上调表达，且二者都受干旱胁迫上调表达。因此，推测MiSTPP3和MiSTPP7可能参与澳洲坚果的生长发育及逆境胁迫反应。

关键词：澳洲坚果;PP5;PP7;基因克隆;表达分析

中图分类号：S664.9      文献标识码：A

Abstract： Two serine / threonine phosphatase genes were named MiSTPP3 and MiSTPP7， with coding sequences 2095 bp and 1700 bp in length， respectively， which were cloned by RT-PCR from macadamia cultivar ‘O.C.’ （Macadamia integrifolia）. Sequence analysis and evolution analysis indicated that MiSTPP3 with the MPP_PP5_C domain belonged to PP5 subfamily， and MiSTPP7 with the MPP_PP7 domain belonged to PP7 subfamily. Physicochemical properties of protein and subcellular localization analysis showed that MiSTPP3 as a stable hydrophilic protein may be located in the plasma membrane， and MiSTPP7 as an unstable hydrophilic protein may be located in the mitochondrial matrix. Bioinformatics analysis indicated that MiSTPP3 and MiSTPP7 were non-secreted transmembrane proteins， α-helix， random coil and extended strand were the backbone of the secondary and 3D protein structures. RT-qPCR analysis indicated that MiSTPP3 and MiSTPP7 expressed differently in roots， stems， leaves， flowers and small fruits of M. integrifolia; MiSTPP3 was down regulated by low temperature stress， while MiSTPP7 was up-regulated， and both of them were up-regulated by drought stress. The results suggested that MiSTPP3 and MiSTPP7 may be involved in the growth and development of M. integrifolia and the stress response.

Keywords： Macadamia ternifolia; PP5; PP7; gene cloning; expression analysis

DOI： 10.3969/j.issn.1000-2561.2021.11.007

蛋白質可逆磷酸化是通过蛋白激酶和蛋白磷酸酶催化蛋白质磷酸化和去磷酸化[1-2]，这一过程几乎涉及细胞所有的生命活动，如生长发育、代谢调控、信号感知与转导、分化与衰老等[3-6]。蛋白磷酸酶催化的去磷酸化负调控蛋白激酶的生化作用，从而抑制蛋白激酶的细胞学功能，并在细胞应答信号通路中起到重要的作用，特别对于调控植物生长发育和响应逆境胁迫具有重要意义。

蛋白磷酸酶是一个庞大的蛋白家族，主要分为丝氨酸/苏氨酸蛋白磷酸酶（PSP）和酪氨酸蛋白磷酸酶（PTP）两大类，根据蛋白质系统演化分析PSP大家族又分为PPP家族和PPM家族，其中PPP家族主要包括PP1、PP2A、PP2B（PP3）、PP4、PP5、PP6、PP7、SLP、Kelch、PPP Unique共10个亚家族[7-8]。这10个亚家族中，只有PP1、PP2A和PP2B的功能研究比较多，而其他亚家族的功能研究鲜有报道[9-12]。PP5和PP7的蛋白组分与结构、理化特性在进化史上的亲缘关系较近[13-14]，与其他亚家族存在明显的差异。PP5蛋白通过TPR结构域，与分子伴侣Hsp90或Hsp70结合从而激活PP5蛋白的活性[15]。PP5蛋白在调控植物耐热性、抗病性等逆境胁迫中发挥重要的作用[16-17]。PP7蛋白是植物中第一个报道与钙调素（CaM）相互作用的丝氨酸/苏氨酸蛋白磷酸酶，在植物耐热性和氧化还原反应调控、光敏色素信号转导及胁迫信号转导中发挥重要作用[13， 18-20]。但是，目前有关PP5和PP7亚家族的功能研究并不多，有待深入研究。

澳洲坚果又称夏威夷果，原产于亚热带雨林，对生长环境条件要求比较高[21-22]。低温、干旱等逆境条件会影响澳洲坚果的正常生长发育，甚至造成澳洲坚果产量下降或死亡[23]。蛋白磷酸酶

PP5和PP7在调控植物生长发育和响应逆境胁迫中发挥关键的作用[14， 20]，但目前國内外关于澳洲坚果PP5和PP7基因的功能研究还未见报道。因此，本研究从澳洲坚果光壳种中克隆获得PP5基因MiSTPP3和PP7基因MiSTPP7，利用生物信息学分析其编码蛋白的结构特征，并利用RT-qPCR分析其在不同组织和不同胁迫中的表达特性，以期为阐明澳洲坚果PP5和PP7基因在澳洲坚果生长发育及响应冷害、干旱逆境胁迫的分子机制提供科学依据。

1  材料与方法

1.1  材料

植物试材为澳洲坚果光壳种（Macadamia integrifolia）品种‘O.C.’，取自广西农业科学院广西亚热带作物研究所澳洲坚果种质资源圃。采集树龄为5～6年的‘O.C.’植株幼根、幼茎、嫩叶、刚开放的小花及谢花后30～45 d的小果。分别对处于苗期（接穗为4个月苗龄）的嫁接苗进行4 ℃低温处理（0.5、1、3、6、12、24 h）和20%浓度的PEG6000干旱处理（1、3、6、12、24、36 h），并采集处理后的嫩叶。其中，无胁迫处理的植株作为对照（CK），每个试验重复3次。

1.2  方法

1.2.1  MiSTPP3和MiSTPP7基因全长cDNA克隆根据蔡元保等[24]的CTAB改良法提取澳洲坚果总RNA，并参照通用型逆转录试剂盒进行澳洲坚果cDNA第1链合成。根据本课题组从澳洲坚果转录组获得的2个蛋白磷酸酶（PP5和PP7）基因，设计扩增这2个基因开放阅读框ORF的PCR引物（表1），其中PP3-F和PP3-R扩增MiSTPP3，以及PP7-F和PP7-R扩增MiSTPP7。取1.5 μL叶片的cDNA模板进行PCR扩增：94 ℃预变性2 min;94 ℃变性35 s，53.5 ℃（MiSTPP3基因）或52.6 ℃（MiSTPP7基因）退火35 s，72 ℃延伸100 s，33个循环;72 ℃再延伸7 min，10 ℃保温。

1.2.2  蛋白磷酸酶的生物信息学分析  利用NCBI网站的BLAST程序检索蛋白磷酸酶的同源序列，用其保守结构域库（CDD）分析蛋白磷酸酶的结构域，用ORF Finder分析蛋白磷酸酶的开放阅读框ORF，用PROSITE软件预测基元和功能结构域，用ProtParam程序分析蛋白磷酸酶的基本理化性质，用NetPhos 3.1在线程序分析蛋白质的磷酸化位点，用PSORT软件预测亚细胞定位，用TMpred软件和SignalP 4.1软件分别在线预测蛋白磷酸酶的跨膜结构和信号肽，用SOPMA软件预测蛋白质的二级结构组分及其含量，采用同源建模使用SWISS-MODEL在线工具构建蛋白质核心区域的三级结构。在NCBI网站GenBank数据库中搜索其他植物蛋白磷酸酶的同源蛋白，采用DNAMAN 6.0软件进行蛋白质多序列比对，并构建蛋白磷酸酶与其他同源蛋白的系统进化树。

1.2.3  MiSTPP3和MiSTPP7基因的表达分析  采用实时荧光定量RT-qPCR方法分析澳洲坚果MiSTPP3和MiSTPP7基因的表达情况，以澳洲坚果苹果酸脱氢酶基因MDH（GenBank登录号：MN627210）为内参基因[25]，引物为QMDH-F和QMDH-R;MiSTPP3基因引物为QPP3-F和QPP3-R;MiSTPP7基因引物为QPP7-F和QPP7-R（表1）。利用RT-qPCR分析澳洲坚果MiSTPP3和MiSTPP7基因在不同组织（根、茎、叶、花和小果）和不同胁迫下（低温和干旱胁迫）的表达情况。RT-qPCR反应程序为：94 ℃预变性3 min;94 ℃变性20 s，56.3～58.0 ℃（表1）退火20 s，72 ℃延伸30 s，40个循环。每个试验重复3次，根据2‒ΔΔCT方法处理数据以确定MiSTPP3和MiSTPP7基因的相对表达量。

1.3  数据处理

利用Excel软件对基因的相对表达量作图，应用SPSS 18.0软件进行数据的差异显著性分析。

2  结果与分析

2.1  MiSTPP3和MiSTPP7基因的全长序列克隆

PCR扩增、核苷酸测序及蛋白同源检索结果显示（图1），获得的澳洲坚果PP3基因序列全长2095 bp，包含1455 bp的完整编码区CDS序列，共编码含有484个氨基酸的多肽;PP7基因序列全长1700 bp，包含1338 bp的完整编码区CDS序列，共编码含有445个氨基酸的多肽，其中MiSTPP3蛋白和丝氨酸/苏氨酸磷酸酶PP5，以及MiSTPP7蛋白和丝氨酸/苏氨酸磷酸酶PP7都具有极高的氨基酸序列相似性。因此，将这2个基因分别命名为MiSTPP3和MiSTPP7，GenBank注册号分别为MT374550和MT374554。

2.2  MiSTPP3和MiSTPP7蛋白氨基酸序列比对

采用DNAMAN 6.0软件对澳洲坚果和其他植物的蛋白磷酸酶进行蛋白质多序列比对，并用PROSITE软件分析其功能结构域。分析结果显示（图2），澳洲坚果MiSTPP3与葡萄VvPP5、玫瑰木属植物RaPP5和辣椒CaPP5的氨基酸序列同源性分别为90.68%、89.63%和89.23%，都含有PP5亚家族的典型结构域MPP_PP5_C和N端的TPR结构域;MiSTPP7与橡胶树HbPP7、麻风树JcPP7和水芙蓉NnPP7的氨基酸序列同源性分别为77.48%、77.48%和77.38%，都含有PP7亚家族的典型结构域MPP_PP7。

2.3  MiSTPP3和MiSTPP7与其他蛋白磷酸酶的系统发育分析

选择澳洲坚果MiSTPP3和MiSTPP7蛋白与不同物种的代表性蛋白磷酸酶构建系统发育树。结果显示（图3），25个丝氨酸/苏氨酸蛋白磷酸酶被聚成PP1、PP2A、PP2B（PP3）、PP5和PP7共5个亚家族，其中澳洲坚果MiSTPP3蛋白属于PP5亚家族，与同一亚家族的玫瑰木属植物RaPP5和葡萄VvPP5的亲缘关系最近;MiSTPP7蛋白属于PP7亚家族，与同一亚家族的水芙蓉NnPP7和枣ZjPP7的亲缘关系最近;而且，在这5个亚家族中，PP5和PP7亚家族的亲缘关系最近，PP5和PP7亚家族又与PP1亚家族的亲缘关系最近，推测具有相同的进化起源和相似的生物学功能。

2.4  MiSTPP3和MiSTPP7蛋白的一级结构和理化性质分析

MiSTPP3蛋白的一级结构和理化性质分析结果显示，其分子质量为54.58 kDa，分子式为C2449H3786N646O719S25，理论等电点（pI）为5.81，带负电荷残基数（Asp+Glu）为69个，带正电荷残基数（Arg+Lys）为58个，不稳定系数为35.42，为稳定蛋白。脂肪族指数为79.79，亲水性系数GRAVY预测为负值（‒0.287），为亲水性蛋白。因此，MiSTPP3蛋白是一个稳定的亲水性蛋白。

MiSTPP7蛋白的一级结构和理化性质分析结果显示，其分子质量为49.30 kDa，分子式为C2211H3433N585O665S14，理论等电点（pI）为5.50，带负电荷残基数（Asp+Glu）为57个，带正电荷残基数（Arg+Lys）为47个，不稳定系数为41.81，为不稳定蛋白。脂肪族指数为81.46，亲水性系数GRAVY预测为负值（‒0.319），为亲水性蛋白。因此，MiSTPP7蛋白是一个不稳定的亲水性蛋白。

2.5  MiSTPP3和MiSTPP7蛋白的基本功能分析

Net Phos 3.1软件预测表明，MiSTPP3蛋白有41个磷酸化位点，包括22个丝氨酸，8个苏氨酸，11个酪氨酸。PSORT软件预测表明，MiSTPP3蛋白定位于质膜可能性最高，为70%。SignalP4.1和TMpred软件预测显示，MiSTPP3蛋白不存在信号肽，有2个跨膜结构域（262～286位，311～331位氨基酸），故MiSTPP3为非分泌跨膜蛋白。

Net Phos 3.1软件预测表明，MiSTPP7蛋白有50个磷酸化位点，包括32个丝氨酸，11个苏氨酸，7个酪氨酸。PSORT软件预测表明，MiSTPP7蛋白定位于线粒体基质可能性最高，为46.8%。SignalP和TMpred软件预测显示，MiSTPP7蛋白不存在信号肽，有3个跨膜结构域（76～99位，155～177位，225～242位氨基酸），故MiSTPP7为非分泌跨膜蛋白。

2.6  MiSTPP3和MiSTPP7蛋白的二级和三级结构预测

SOPMA软件预测显示，澳洲坚果MiSTPP3蛋白二级结构的主要结构元件是α-螺旋、无规则卷曲、延伸链和β-转角，分别占46.28%、30.99%、15.70%和7.02%，α-螺旋主要集中在N端（图4A）;MiSTPP7蛋白二级结构的主要结构元件是无规则卷曲、α-螺旋、延伸链和β-转角，分别占46.74%、30.34%、16.40%和6.52%（图4B）。

以丝氨酸/苏氨酸磷酸酶PP5的晶体结构（PDB ID：5obl.1.A）为同源模板，利用在线工具SWISS-MODEL预测MiSTPP3蛋白的三级结构。结果表明，MiSTPP3蛋白的氨基酸序列与模板序列一致性为85.74%，GMQE值为0.98，同源建模区域的组成以α-螺旋和无规卷曲为主，以延伸链为辅（图5A）。目前，SWISS-MODEL软件的蛋白数据库中还没有丝氨酸/苏氨酸磷酸酶PP7的晶体结构，只能以亲缘关系最近的PP5晶体结构（PDB ID：5hpe.1.A）为同源模板，利用SWISS-MODEL在线工具构建MiSTPP7蛋白的三级结构。结果表明，MiSTPP7蛋白的氨基酸序列与模板序列一致性為36.94%，GMQE值为0.55，同源建模区域主要由无规卷曲和α-螺旋组成（图5B）。与MiSTPP3和MiSTPP7蛋白的二级结构预测结果相符合。

2.7  MiSTPP3和MiSTPP7基因的表达分析

MiSTPP3和MiSTPP7基因在不同组织的表达分析显示（图6），MiSTPP3和MiSTPP7在澳洲坚果根、茎、叶、花和小果中呈现显著性差异表达，其中MiSTPP3在小果中的表达量最高，MiSTPP7在花中的表达量最高。4 ℃低温胁迫的表达分析显示（图7），澳洲坚果苗期受低温胁迫后MiSTPP3在叶片中表达显著性下调;而MiSTPP7表达显著上调，且胁迫3 h后的表达量达到峰值。PEG模拟干旱胁迫的表达分析显示（图8），澳洲坚果苗期受干旱胁迫后叶片中MiSTPP3 和MiSTPP7的表达都显著上调，且在胁迫处理12 h时达到峰值。

3  讨论

目前，PP5和PP7亚家族基因在动物中的研究较多，在细胞生长发育调控、响应逆境胁迫以及细胞应答信号通路中扮演重要角色[26-27]，但在植物中这2个亚家族基因的鉴定、结构及功能研究较少报道。PP5和PP7是PPP家族的独特成员[7-8]，MiSTPP3和MiSTPP7分别作为PP5和PP7亚家族的新成员，分别含有典型结构域MPP_PP5_C和MPP_PP7。

蛋白去磷酸化是蛋白可逆磷酸化修饰的2个关键步骤之一[5， 28]，MiSTPP3和MiSTPP7蛋白在发挥去磷酸化作用的过程中主要以丝氨酸进行磷酸化修饰。植物PP5与PP7蛋白的亚细胞定位取决于mRNA前体的选择性剪接[29]。MiSTPP3和MiSTPP7作为非分泌跨膜蛋白，分别定位于质膜和线粒体基质，至于其亚细胞定位的mRNA前体剪接机制还有待于深入研究。

PP5与PP7蛋白在进化过程中都具有保守的氨基酸序列，二者亲缘关系相对于其他亚家族最为相近[14]，具有相同的进化起源和相似的生物学功能。已有研究表明，这2个亚家族基因在调控组织的生长发育中发挥重要作用，如调控拟南芥气孔的生长发育[29]。MiSTPP3和MiSTPP7可能参与澳洲坚果不同组织（如根、茎、叶、花和小果等）的生长发育，尤其在调控花和小果生长发育中发挥重要作用。

分子伴侣Hsp90是PP5蛋白催化的去磷酸化的关键酶，两者的结合可调控胁迫诱导的细胞信号转导[14]。如热激胁迫条件下，拟南芥AtPP5与底物MDH相互作用，与AtHsp90结合形成复合体从而增强植株的耐热性[16];番茄R蛋白I-2通过与PP5和Hsp90的直接相互作用，增强了植株对尖孢镰刀菌的抗性[30]。MiSTPP3受低温诱导后显著下调表达，受干旱诱导后显著上调表达，这2种不同的表达模式都证实该基因参与澳洲坚果逆境胁迫反应。但是，MiSTPP3是否与Hsp90结合从而响应这些逆境胁迫反应还需深入研究。

植物PP7蛋白在光和胁迫信号传递途径的交叉点上发挥了重要的调控作用[20]。拟南芥AtPP7除了调控氧化还原反应外[13]，还与钙调素（CaM）结合从而增强植株的耐热性[18]，且作为正调控因子调控其光敏色素信号转导[19，31]。与已知的PP7亚家族基因一样，MiSTPP7受低温和干旱诱导后都显著上调表达。因此，推测MiSTPP7可能作为正调控因子，参与澳洲坚果逆境胁迫的信号转导。

本研究结果表明，MiSTPP3和MiSTPP7可能参与了澳洲坚果植株生长发育和逆境胁迫响应。在后续的工作中，利用转基因技术对澳洲坚果MiSTPP3和MiSTPP7基因功能的深层次研究，有助于阐明这2个基因在响应逆境胁迫的生物学功能和作用机制，尤其是对于揭示植物PP5和PP7亚家族基因的信号转导网络及其亚家族间的相互关系具有重要作用。

参考文献

[1] Macek B， Mann M， Olsen J V. Global and site-specific quantitative phosphoproteomics： principles and applica-tions[J]. Annual Review of Pharmacology and Toxicology， 2009， 49（1）： 199-221.

[2] Goel R， Harsha H C， Pandey A， et al. Human protein refer-ence database and human proteinpedia as resources for phosphoproteome analysis[J]. Molecular BioSystems， 2012， 8（2）： 453-463.

[3] Ptacek J， Devgan G， Michaud G， et al. Global analysis of protein phosphorylation in yeast[J]. Nature， 2005， 438（7068）： 679-684.

[4] Aivaliotis M， Macek B， Gnad F， et al. Ser/thr/Tyr protein phosphorylation in the archaeon Halobacterium salina-rum— A representative of the third domain of life[J]. PLoS One， 2009， 4（3）： e4777.

[5] Barrero-Gil J， Salinas J. Post-translational regulation of cold acclimation response[J]. Plant Science， 2013， 205/206： 48- 54.

[6] 肖  冬， 崔燕娇， 王宁宁. 叶片衰老过程中的蛋白激酶和蛋白磷酸酶[J]. 植物生理學报， 2014， 50（9）： 1267-1273.

[7] Cohen P. The structure and regulation of protein phospha-tases[J]. Annual Review of Biochemistry， 1989， 58（1）： 453- 508.

[8] Owens D M， Keyse S M. Differential regulation of MAP kinase signalling by dual-specificity protein phosphatases[J]. Oncogene， 2007， 26（22）： 3203-3213.

[9] Durian G， Rahikainen M， Alegre S， et al. Protein phospha-tase 2A in the regulatory network underlying biotic stress resistance in plants[J]. Frontiers in Plant Science， 2016， 7： 812.

[10] Khan Z H， Agarwal S， Rai A， et al. Co-expression network analysis of protein phosphatase 2A （PP2A） genes with stress-responsive genes in Arabidopsis thaliana reveals 13 key regulators[J]. Scientific Reports， 2020， 10（1）： 1-16.

[11] Tan S T， Abas M， Verstraeten I， et al. Salicylic acid targets protein phosphatase 2A to attenuate growth in plants[J]. Current Biology， 2020， 30（3）： 381-395.e8.

[12] Casamayor A， Ariño J. Controlling Ser/Thr protein phosphatase PP1 activity and function through interaction with regulatory subunits[M]//Advances in Protein Chemistry and Structural Biology. Amsterdam： Elsevier， 2020： 231-288.

[13] Andreeva A V， Solov'eva O V， Kakuev D L， et al. Purifica-tion of plant protein phosphatase PP7 and evidence for its redox regulation[J]. Archives of Biochemistry and Biophysics， 2001， 396（1）： 65-70.

[14] Oberoi J， Dunn D M， Woodford M R， et al. Structural and functional basis of protein phosphatase 5 substrate specificity[J]. PNAS， 2016， 113（32）： 9009-9014.

[15] Connarn J N， Assimon V A， Reed R A， et al. The molecular chaperone Hsp70 activates protein phosphatase 5 （PP5） by binding the tetratricopeptide repeat （TPR） domain[J]. Journal of Biological Chemistry， 2014， 289（5）： 2908-2917.

[16] Park J H， Kim W Y， Chae H B， et al. Serine/threonine pro-tein phosphatase 5 （PP5） interacts with substrate under heat stress conditions and forms protein complex in Arabidop-sis[J]. Plant Signaling & Behavior， 2012， 7（5）： 535-538.

[17] Van Bentem S D L F， Vossen J H， Vries K J， et al. Heat shock protein 90 and its co-chaperone protein phosph¬atase 5 interact with distinct regions of the tomato I-2 disea¬se resistance protein[J]. Plant Journal， 2005， 43（2）： 284-298.

[18] Liu H T， Li G L， Chang H， et al. Calmodulin-binding protein phosphatase PP7 is involved in thermotolerance in Arabidopsis[J]. Plant Cell & Environment， 2007， 30（2）： 156-164.

[19] Genoud T， Treviño Santa Cruz M， Kulisic T， et al. The protein phosphatase 7 regulates phytochrome signaling in Arabidopsis[J]. PLoS One， 2008， 3（7）： e2699.

[20] Andreeva A V， Kutuzov M A. PPEF/PP7 protein Ser/Thr phosphatases[J]. Cellular and Molecular Life Sciences， 2009， 66（19）： 3103-3110.

[21] Buthelezi N M D， Magwaza L S， Tesfay S Z. Postharvest pre-storage processing improves antioxidants， nutritional and sensory quality of Macadamia nuts[J]. Scientia Horticulturae， 2019， 251： 197-208.

[22] Herbert S W， Walton D A， Wallace H M. The influence of pollen-parent and carbohydrate availability on Macadamia yield and nut size[J]. Scientia Horticulturae， 2019， 251： 241-246.

[23] 蔡元保， 楊祥燕， 陈显国， 等. 澳洲坚果SCoT反应体系的建立及应用[J]. 热带亚热带植物学报， 2013， 21（3）： 253- 258.

[24] 蔡元保， 杨祥燕， 孙光明， 等. 菠萝花发育相关基因AcMADS1的克隆与组织表达特性分析[J]. 植物学报， 2014， 49（6）： 692-703.

[25] 杨  倩， 杨子平， 周娅丽， 等. 澳洲坚果实时荧光定量PCR分析中内参基因的筛选[J]. 热带作物学报， 2020， 41（8）： 1505-1512.

[26] Kostich M， English J， Madison V， et al. Human members of the eukaryotic protein kinase family[J]. Genome Biology， 2002， 3（9）： research0043.1.

[27] 邹  奕， 闫彩燕， 邳  植. 甜菜蛋白磷酸酶超家族分析[J]. 中国糖料， 2018， 40（6）： 19-22.

[28] 阮班军， 代  鹏， 王  伟， 等. 蛋白质翻译后修饰研究进展[J]. 中国细胞生物学学报， 2014， 36（7）： 1027-1037.

[29] Van Bentem S D L F， Vossen J H， Vermeer J E M， et al. The subcellular localization of plant protein phosphatase 5 isoforms is determined by alternative splicing[J]. Plant Physiology， 2003， 133（2）： 702-712.

[30] Sun X D， Kang X J， Ni M. Hypersensitive to red and blue 1 and its modification by protein phosphatase 7 are im-plicated in the control of Arabidopsis stomatal aperture[J]. PLoS Genetics， 2012， 8（5）： e1002674.

[31] Møller S G， Kim Y S， Kunkel T， et al. PP7 is a positive regulator of blue light signaling in Arabidopsis[J]. Plant Cell， 2003， 15（5）： 1111-1119.

責任编辑：黄东杰