Identification of CIPK Family in Rice and qRT-PCR Analysis on OsCIPK5 Induced by Magnaporthe oryzae
-
摘要:目的 植物通过启动一系列信号传导过程来应对外部环境，这些过程通常涉及多种蛋白激酶，包括钙调神经磷酸酶B样蛋白互作激酶（calcineurin B-like protein-interacting protein kinases, CIPKs）。为更加清晰全面地了解水稻CIPK基因家族，本研究根据最新的基因组测序数据对水稻基因组中的CIPKs进行了鉴定。方法 通过探讨拟南芥和水稻中CIPKs蛋白家族的结构特点，结合生物信息学和qRT-PCR技术系统分析水稻中CIPKs家族蛋白的结构。结合转录组数据，比较了粳稻云引受稻瘟病菌诱导后的表达情况。结果 根据最新的水稻基因组数据，鉴定出31个水稻OsCIPK基因。系统发育树分析结果表明，31个OsCIPK基因可分为5个亚家族，这些亚家族具有不同的外显子-内含子和UTR的结构特点。从广谱抗稻瘟病品种粳稻云引受稻瘟病菌诱导的基因表达谱的趋势聚类中筛选出了OsCIPK5基因并对其进行了表达分析，结果表明，云引中OsCIPK5基因受稻瘟病菌的诱导表达。结论 内含子缺失和片段重复在水稻OsCIPK基因家族的扩展中起到重要作用，同时OsCIPK5受到稻瘟病菌的诱导表达。Abstract:Objective The family of calcineurin B-like protein-interacting protein kinase genes (CIPKs) in rice was studied using the latest sequencing data to further understand the signal transduction involving a variety of kinases in plants in respond to environmental changes.Methods The structures of CIPKs in Arabidopsis thaliana and rice were analyzed. Combining the bioinformatics and qRT-PCR technology, expressions of the CIPK familie and that of the genes induced by M. grisea in Japonica rice cv. 'Yunyin' were compared.Results A total of 31 OsCIPK genes were identified in the rice genome databank. The phylogenetic tree analysis showed that these 31 OsCIPKs could be divided into 5 subfamilies, which had different structural characteristics of exon-intron and UTR. The expression of OsCIPK5, which was selected from a trend cluster of the gene expressiong profile of 'Yunyin' induced by M. grisea, could be induced by M. grisea.Conclusions It appeared that the intron deletion and fragment duplication played an important role in the expansion of OsCIPK family in rice, and that OsCIPK5 expression was induced in 'Yunyin' by M. grisea.
-
0. 引言
【研究意义】植物利用多种信号转导途径对外界环境做出响应，细胞内钙浓度的瞬时升高是机体对非生物胁迫和生物胁迫等多种信号响应的主要事件[1-5]。高等植物表达几个主要的钙传感器家族，包括钙调蛋白（calmodulin, CaM）、CaM相关蛋白、CBL蛋白以及钙依赖蛋白激酶（calcium-dependent protein kinases, CDPKs）[6-10]。CBL蛋白通过激活靶蛋白，如CBL互作蛋白激酶（CBL-interacting protein kinases, CIPKs）来触发Ca2+介导的信号网络。植物CIPK的N端含有一个Ser/Thr蛋白激酶结构域，C端含有一个NAF结构域。拟南芥中的CIPK蛋白最初被称为SOS2蛋白激酶(Sail overly sensitive 2, SOS2)，激活的CIPK可通过磷酸化下游信号元件来转导钙信号[11]。在拟南芥中，含有SOS3（CBL蛋白）[12-13]、SOS2（CIPK蛋白）[12-14]和SOS1（被SOS2激活的Na+/H+反向运输载体）[15-17]的SOS通路在Ca2+信号转导中的作用已被详细解析，其可以在盐胁迫条件下维持离子稳态[11, 13, 15, 18- 19]。【前人研究进展】前人基于基因组的分析表明，拟南芥中存在25个AtCIPK基因，水稻中存在30个OsCIPK基因[20-21]。Kolukisaoglu等[20]报道了30个水稻基因组中的OsCIPK基因，并将它们分别命名为OsCIPK1至OsCIPK30。研究表明，干旱、盐碱化、低温、聚乙二醇（PEG）或脱落酸（abscisic acid, ABA）等非生物胁迫可诱导20个OsCIPK基因的转录[21]。3个OsCIPK基因（OsCIPK3、OsCIPK12和OsCIPK15）的过表达分别显著提高了转基因植物对低温、干旱和盐胁迫的耐受性[21]。Chen等[22]分离鉴定了15个OsCIPK基因，基因表达分析表明，这些基因可不同程度地被生物胁迫，如：白叶枯病菌和非生物胁迫（重金属、高盐度、低温和ABA）诱导表达。其中，经白叶枯病菌感染后，5个OsCIPK基因（OsCIPK1、2、10、11和12）转录上调，另有4个OsCIPK基因（OsCIPK2、10、11和14）受所有胁迫处理诱导，表明这些OsCIPK基因参与了多个胁迫反应途径和/或植物的抗病反应。【本研究切入点】虽然CIPK基因家族在许多植物中得到了广泛的研究，但是目前对水稻CIPK蛋白的认识还很有限。由于CIPK基因在植物多种生理功能中发挥重要作用，因此，对水稻CIPK家族成员进行系统的研究具有重要意义。【拟解决的关键问题】在本研究中，我们利用Ensembl Plants(http://plants.ensembl.org/)数据库和UniProt(https://www.uniprot.org/uniprot/) 蛋白数据库的最新基因组数据鉴定出了31个水稻CIPK基因，进一步研究了该基因家族成员的外显子-内含子的结构、基序组成、基因复制、染色体分布、系统发育和同源性等，并分析了粳稻云引中OsCIPK5基因受稻瘟病菌诱导的表达情况，为水稻CIPK基因家族成员的功能鉴定提供了有价值的线索。
1. 材料与方法
1.1 水稻CIPK基因家族的鉴定
在Ensembl Plants(http://plants.ensembl.org/)数据库下载水稻和拟南芥基因组序列文件，用TBtools软件从基因组序列中提取CDS序列，再转换成蛋白质序列。
在 UniProt(https://www.uniprot.org/uniprot/) 蛋白数据库下载水稻或拟南芥CIPK基因家族蛋白序列，用HMMER 3.0软件建模并搜索水稻或拟南芥中CIPK基因家族的成员，去掉序列重复的成员，筛选后的家族成员在NCBI (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/) 预测保守结构域，根据CIPK基因家族的保守结构域，验证基因家族中序列异常的成员，并对其进行手动校正。等电点、分子量等数据由ExPasy (http://web.expasy.org/protparam/)网站在线工具计算得来。
1.2 水稻CIPK基因家族结构分析
用MEGAX软件的极大似然估计法(Maximum Likelihood)构建31个水稻CIPK基因家族的系统发育树。bootstrap replications参数为1 000，模型为LG+G。CIPK蛋白序列用MEME（http://meme.nbcr.net/meme/intro.html）在线工具分析水稻CIPK基因家族保守motif，motif数为10。用TBtools软件对水稻CIPK基因结构、保守结构域可视化，再将上述分析结果合并为一张图。
1.3 水稻CIPK基因家族进化分析
1.3.1 CIPK基因在染色体上的分布
通过NCBI(ftp://ftp.ncbi.nlm.nih.gov/blast/)下载的blast软件对水稻基因组进行自身比对，再用MCScanX(Multiple Collinearity Scan toolkit)软件对比对结果进行共线性分析，提取水稻CIPK基因的直系同源基因ID。用TBtools软件对CIPK基因在水稻12条染色体上的分布以及CIPK基因中直系同源基因关系进行可视化。
1.3.2 水稻和拟南芥CIPK基因的进化与分类
用MEGAX软件的极大似然估计法(Maximum Likelihood)构建31个水稻CIPK基因和26个拟南芥CIPK基因的系统发育树，设置参数同上。导出文件，用Figtree软件对水稻和拟南芥CIPK基因分类进行注释与可视化。
1.4 qRT-PCR分析
利用TriPure Isolation Reagent提取总RNA，选择Actin为内参基因，下载Actin和OsCIPK5的CDs序列，设计引物如表1所示。
表 1 OsCIPK5和Actin的引物序列Table 1. The primer sequence of OsCIPK5 and Actin基因 Genes q-F q-R OsCIPK5 CAGAGCGTCGCCATCAAGGTC GTGACAGAAGTCCACAGCCCCTATC Actin CCTCGTCTGCGATAATGGAACTG CCCTGGGCGCATCATCTC 采用喷雾接种法将稻瘟病菌孢子液分别喷洒到粳稻云引三叶一心期幼苗的叶片上，分别于0、12、24、36、48、72、96、120、144、168 h共10个时间点取云引叶片样品，用液氮速冻，保存于−80 ℃冰箱中备用。提取的叶片RNA用High-Capacity cDNA Archive Kit反转成单链cDNA，将cDNA稀释20倍用于qRT-PCR分析。通过ABI Prism® 7500 型荧光定量PCR仪进行数据处理，之后通过仪器自带软件生成柱形图。
2. 结果与分析
2.1 水稻CIPK基因的鉴定、染色体分布及同源性分析
通过BLAST和Hidden Markov Model搜索，对OsCIPK蛋白进行了鉴定。从水稻基因组中共鉴定出31个非冗余CIPK。保守域分析证实，所有的OsCIPK都能识别出S_TKc结构域，这是CIPK家族的基本特征。此外，对CIPK蛋白的多重序列比对表明，本研究中鉴定的所有OsCIPK在N端均含有激酶结构域，29个OsCIPK在C端均含有CIPK_C结构域。由于部分基因含有选择性mRNA剪接，因此选择每个基因最长的蛋白进行进一步分析。31个预测的OsCIPK蛋白的氨基酸残基范围为405～541个，相对分子质量范围为43.97～59.81 kDa（表2）。
表 2 水稻基因组CIPK家族的特征Table 2. Characteristics of CIPK family in rice genome基因名称
Name of gene转录本ID(RAP-DB)
Transcriptome染色体位置
Locus of chromosome氨基酸长度
Length of amino acid等电点(PI)
Isoelectric point分子量(WM)
Weight of molecular/kDaOsCIPK1 Os01t0292200-01 chr1:10622951-10627532 462 6.26 52.20 OsCIPK2 Os07t0678600-01 chr7:28726817-28729527 444 9.18 50.27 OsCIPK3 Os07t0687000-01 chr7:29191197-29194187 446 6.9 50.95 OsCIPK5 Os01t0206700-03 chr1:5809589-5811389 462 9.28 51.97 OsCIPK6 Os08t0441100-01 chr8:21467453-21468593 305 5.44 32.48 OsCIPK7 Os03t0634400-01 chr3:24226372-24227930 448 9.35 48.37 OsCIPK8 Os01t0536000-01 chr1:19469262-19477590 447 6.45 50.60 OsCIPK9 Os03t0126800-00 chr3:1515355-1518947 479 8.31 53.78 OsCIPK10 Os03t0339900-01 chr3:12630705-12635182 452 9.01 51..49 OsCIPK11 Os01t0824600-01 chr1:35228913-35231725 503 8.48 56.56 OsCIPK12 Os01t0759400-02 chr1:31943529-31948580 541 8.3 59.81 OsCIPK13 Os01t0206300-00 chr1:5799607-5801142 512 8.05 56.01 OsCIPK14 Os12t0113500-01 chr12:679079-682143 440 9.44 50.32 OsCIPK15 Os11t0113700-01 chr11:630585-633622 435 9.53 49.70 OsCIPK16 Os09t0418000-03 chr9:15009182-15010716 457 8.84 50.50 OsCIPK17 Os05t0136200-01 chr5:2113595-2117110 455 6.93 50.92 OsCIPK18 Os05t0332300-01 chr5:15556085-15560404 458 8.84 51.50 OsCIPK19 Os05t0514200-01 chr5:25489349-25491319 509 7.26 56.91 OsCIPK20 Os05t0208100-01 chr5:6701053-6702933 467 8.29 51.35 OsCIPK21 Os07t0637000-01 chr7:26464444-26468858 445 7.96 50.10 OsCIPK22 Os05t0334750-00 chr5:15643309-15644732 452 8.03 49.32 OsCIPK23 Os07t0150700-01 chr7:2643540-2647773 451 9.23 50.71 OsCIPK24 Os06t0606000-01 chr6:24038959-24044785 454 8.52 50.94 OsCIPK25 Os06t0543400-01 chr6:20490388-20492271 515 8.78 56.98 OsCIPK26 Os02t0161000-01 chr2:3283264-3286366 494 9.13 55.91 OsCIPK27 Os09t0418500-00 chr9:15027311-15034300 809 5.53 87.30 OsCIPK28 Os05t0476350-00 chr5:23430891-23432198 436 9.32 49.33 OsCIPK29 Os07t0678300-01 chr7:28710749-28712470 444 8.69 48.21 OsCIPK30 Os01t0759200-01 chr1:31938559-31940184 477 9.34 53.56 OsCIPK31 Os03t0319400-01 chr3:11526273-11530306 450 7.98 50.95 OsCIPK33 Os11t0134300-01 chr11:1605182-1609130 406 6.13 46.60 图1显示，除连锁群（linkage groups, LG）4和10无OsCIPK基因外，OsCIPK基因不均匀地分布于水稻10个LG。LG1的OsCIPK基因数量最多，共有7个。有些连锁群基因较多（如：LG1、LG5），而另一些连锁群基因较少；有些LG只有一个基因（如：LG2）。LG长度与OsCIPK基因数量无显著正相关。根据Holub[23]的描述，在200 kb内包含两个或多个基因的染色体区域被定义为串联重复事件。在LG1、5、7和9上，有10个OsCIPK基因（OsCIPK12/OsCIPK30、OsCIPK5/OsCIPK13、OsCIPK22/OsCIPK18、OsCIPK2/OsCIPK29和OsCIPK16/OsCIPK27）聚集成5个串联重复区域。LG1有两个簇，表明其为OsCIPK基因的一个分布热点。
![]()
图 1 水稻CIPK基因染色体分布及染色体间关系注：红色线表示直系同源的CIPK基因对。染色体数显示在每条染色体的底部Figure 1. Chromosomal distribution and inter-chromosomal relationship of CIPKs in riceNote: Red line indicates homology pair of CIPKs. Number of chromosomes is shown at bottom of each chromosome.为了更好地理解OsCIPK基因家族的进化程度，计算了OsCIPK基因对的Ka/Ks比值（表3）。所有直系同源的OsCIPK基因对Ka/Ks＜1，说明水稻OsCIPK基因家族在进化过程中可能经历了较强的纯化选择性压力。
表 3 水稻基因组中的直系同源CIPK基因对的Ka/Ks比值Table 3. Ka/Ks value of orthologous CIPK gene pairs in rice genome直系同源CIPK基因对
The homologous CIPK gene pairKa Ks Ka/Ks OsCIPK31/OsCIPK3 0.113 65 1.115 82 0.101 85 OsCIPK15/OsCIPK14 0.004 17 0.035 81 0.116 45 OsCIPK11/OsCIPK28 0.161 62 0.880 47 0.183 55 OsCIPK5/OsCIPK20 0.221 63 1.297 74 0.170 78 OsCIPK6/OsCIPK27 0.348 63 0.529 36 0.658 59 2.2 水稻CIPK基因的结构和保守基序
为了更好地了解水稻CIPK基因家族的结构特征，本研究根据系统发育关系分析了内含子/外显子分布、保守基序和顺式元件。系统发育分析表明，31个OsCIPK家族成员被分为5组（A、B、C、D和E；图2）。含有10～21个外显子的多外显子OsCIPK基因被聚集到A组，而只有1～4个外显子的缺外显子成员被聚集到其他4组（B、C、D和E；图2）。各类群保守的外显子数目支持了它们之间的密切进化关系和类群的分类。在拟南芥中，CIPK基因也被分为外显子丰富组和外显子贫乏组[20]。目前在玉米、杨树等单子叶植物和双子叶植物中发现的CIPK基因也存在类似的现象[24-25]，说明了CIPK基因家族结构的保守性。在植物进化的早期，内含子的产生率较高[26]。此外，根据Nuruzzaman等[27]的研究，分段复制后内含子丢失的速度快于内含子获得的速度。因此，A组可能代表CIPK家族的原始基因。此外，除了OsCIPK27，所有B组的缺外显子OsCIPK含有极短或没有5′-UTR，3′-UTR也短；除了密切相关的基因OsCIPK5/OsCIPK20以及OsCIPK28，C组OsCIPK的5′-UTR长且分段；除了OsCIPK13，D组的5′-UTR也极短，但3′-UTR长；E组的UTR虽短，但两端都有。总的来说，同一组中的基因结构都基本一致，说明CIPK组别的分类与基因结构的特性有关。
![]()
图 2 水稻CIPK基因的系统发育分类、基因结构及保守蛋白基序的特征注：① a为系统发育树，b为蛋白基序结构特征，c为基因结构；② OsCIPK系统发育树分为五组：A、B、C、D和E，用不同的颜色标出；③b图标尺为氨基酸长度，c图标尺为核苷酸长度；④绿色方框表示UTR，黄色方框表示CDS，黑线表示内含子，其他颜色方框代表保守结构域。Figure 2. Phylogenetic relationships, gene structure and architecture of conserved protein motifs in CIPKs from riceNote: ① a, phylogenetic tree; b, protein motif structure; and c, gene structure. ② OsCIPK phylogenetic tree was divided into 5 groups -A, B, C, D, and E- shown in different colors. ③ The ruler of b is for amino acid length; and the ruler of c is for nucleotide length. ④ Green box shows UTR; yellow box shows CDs; black line shows introns; and boxes of other colors show conservative domains.利用基序分析进一步支持了水稻OsCIPK的系统发育关系和分类（图2）。利用MEME软件对水稻OsCIPK的10个保守基序进行了特征分析。值得注意的是，所有水稻OsCIPK蛋白都含有标记为S_TKc domain的1、2、3、4、5、9 motif，而除OsCIPK22和OsCIPK29外的29个水稻OsCIPK蛋白都含有与CIPK_C domain密切相关的6、7、8、10 motif。虽然OsCIPK29缺失了motif 6，但它也含有其他基序赋予的S_TKc domain和CIPK_C domain。虽然OsCIPK22 C端为AMPKA_C_like domain，但OsCIPK22 C端基序排列和结构与OsCIPK29一样，说明AMPKA_C_like domain与CIPK_C domain的功能是一样的。CIPK蛋白具有N端催化激酶结构域和含有NAF基序的C端调控结构域的保守特性[25]。因此，本研究鉴定的所有水稻OsCIPK均具有CIPK家族的保守特征。保守的NAF基序介导了CBL-CIPK之间的物理互作，表明所鉴定的OsCIPK可能在CBL-CIPK网络中存在功能上的连接[28]。另外，在OsCIPK22和OsCIPK29这对旁系同源序列中，缺乏同样注释为NAF域的motif 6。
2.3 水稻和拟南芥CIPK基因家族的进化与分类
为了研究水稻OsCIPK蛋白与其他物种已知CIPK蛋白的进化关系，以水稻和拟南芥CIPK家族蛋白的全氨基酸为基础，构建了系统发育树（图3）。结果表明，根据序列相似性，57个CIPK蛋白可分为5个不同的类群（A～E），这与前人对拟南芥和胡杨CIPK的报道一致[24]。这31个OsCIPK被分配到5个具有较高自举值（bootstrap values，BV）的组，同时还有26个AtCIPK。油菜中CIPK的分类也出现了类似的现象，其中油菜CIPK与拟南芥、水稻、玉米等物种的CIPK共分为4组[29]。
![]()
图 3 水稻和拟南芥中CIPK蛋白的系统发育分析Figure 3. Phylogenetic trees constructed with 26 CIPKs from A. thaliana and 31 CIPKs from rice系统发育分析还发现了水稻与拟南芥亲缘关系密切的同源CIPK（A组中OsCIPK24与AtCIPK24、OsCIPK8与AtCIPK8、OsCIPK23与AtCIPK23、OsCIPK9与AtCIPK9；B组中OsCIPK16与AtCIPK16；C组中OsCIPK6与AtCIPK6、OsCIPK20与AtCIPK20；D组中OsCIPK7与AtCIPK7），表明在水稻和拟南芥分化之前就存在着一组古老的CIPK基因。AtCIPK24/AtSOS2与OsCIPK24高度相似，可以与AtCBL4/AtSOS3相互作用，作用于Na+/H+的反向运输载体AtSOS1/AtNHX7，增强拟南芥根系的耐盐胁迫能力[30]。AtCIPK8是OsCIPK8的同源基因，参与调控初级硝酸盐反应的低亲和阶段[31]。这些证据表明CIPK基因可能在水稻中发挥作用。
2.4 稻瘟病菌对水稻OsCIPK5基因表达的诱导
在植物中，钙在各种信号转导途径中起着普遍的信使作用，包括对生物[32-33]和非生物胁迫[3,34]的反应，以及对细胞和发育过程的调控[35-36]。在水稻中，已有一些基于表达模式和生化数据的报道推测了OsCIPK可能的功能。OsCIPK23对授粉和干旱胁迫均有一定的作用[37]；OsCIPK19响应光和营养物质[38]；OsCIPK31的表达受低温、盐、光、细胞分裂素和糖等多种信号的影响[39]；OsCIPK15在水稻缺O2耐受性中发挥了重要作用[40]；OsCIPK1、OsCIPK2、OsCIPK10、OsCIPK11和OsCIPK12在接种白叶枯病菌后表达显著上调[22]。这些研究表明OsCIPK可能参与多种逆境信号的转导，但到目前为止还未有报道OsCIPK基因受稻瘟病菌的影响。本课题组利用Self-Organizing Map (SOM)方法对抗稻瘟病粳稻品种云引在不同侵染时间下的基因表达差异进行了分析，发现OsCIPK5基因显著上调 [41]，因此从中筛选出OsCIPK5基因并鉴定其表达是否受到稻瘟病菌的影响。随着侵染时间延长，OsCIPK5基因的的表达水平呈升高趋势。进一步进行qRT-PCR验证，分析其在稻瘟病菌侵染后10个时间点[0、12、24、36、48、72、96、120、144、168 hai（hours after infection）]的表达模式（图4）。结果表明，OsCIPK5基因在稻瘟病菌接种前表达量最低，在接种稻瘟病菌后表达量开始上调，但24 hai后基因的表达量下降，到36 hai时，基因的表达丰度甚至略低于0 hai，之后基因的表达量又上调，48 hai的表达水平高于前四个时间点，但之后表达又下调，至72 hai到达低值，但高于0、12、24、36 hai的表达量，后基因的表达水平又上调，并在120 hai时表达丰度达到最大，之后OsCIPK5基因的表达急剧下降，在168 hai时几乎不表达。
![]()
图 4 OsCIPK5在稻瘟病菌侵染后不同时间段的表达情况Figure 4. Expressions of OsCIPK5 at times after M. grisea induction3. 讨论与结论
水稻是最重要的粮食作物之一，也是单子叶模式植物。在本研究中，我们对水稻OsCIPK家族进行了全基因组的分析。在水稻基因组数据库的计算机分析中，通过对保守结构域和多重序列的比对分析，鉴定出了31个OsCIPK基因。通过对水稻和拟南芥的系统发育分析，将其分为5个亚家族。多重序列比对显示5组OsCIPK蛋白在N端均含有激酶结构域，在C端均含有CIPK_C/NAF结构域，说明水稻CIPK蛋白的高度保守性。在这些保守结构域的作用下，CIPK蛋白表现出对相应CBL蛋白的结合偏好。根据以往的研究，为响应各种刺激，CIPK可特异性地靶向结合CBL，并转导捕捉到的钙信号，这属于Ca2+介导的CBL-CIPK网络[29, 42]。
内含子的得失是CIPK基因家族扩展的一种发散力。基因结构分析表明，31个基因中有10个OsCIPK具有10～21个内含子，被聚类为同一组，表明其具有较高的保守性。在拟南芥CIPK家族中，有8个基因含有9～12个保守的内含子，其中7个是保持不变的[20]。同样，杨树中有9个CIPK基因含有12～14个内含子，其中7个是非常保守的[24]。因此，基于内含子数量的不同，水稻、拟南芥和杨树中CIPK家族成员可分为两个不同的基因簇，这进一步说明了植物基因组具有共同的祖先，进而分化为单子叶和双子叶谱系。在本研究中，水稻OsCIPK缺外显子成员又被分为4组，它们之间的区别主要集中在5′-和3′-UTR。
连锁群和进化选择压分析表明，串联和分段复制事件在水稻OsCIPK家族的扩展中也发挥了重要作用（图1，表3）。全基因组复制事件是被子植物进化过程中常见的现象，通常会导致基因家族的扩大[43]。此外，LG1和LG5中OsCIPK基因的密集分布也是导致其基因家族规模多样性的潜在原因，这说明这两个连锁群为OsCIPK基因家族进化中最具活力的群体。3个串联重复基因对（OsCIPK12/OsCIPK30、OsCIPK5/OsCIPK13和OsCIPK22/OsCIPK18）和2个直系同源基因对（OsCIPK11/OsCIPK28和OsCIPK5/OsCIPK20）分布在这两个连锁群上，这可以为OsCIPK基因家族的进化历史提供更多线索。在本研究中，我们还发现亲缘关系近的水稻OsCIPK成员倾向于散在分布，说明植物CIPK基因的多样性低。
CBL-CIPK通路作为植物钙信号转导网络的重要组成部分，进一步分析拟南芥、水稻等物种的应激诱导的CBL和CIPK家族成员的功能，有助于我们更好地理解应激反应中钙信号的分子机制[21, 44-47]。本研究发现，OsCIPK5基因受到稻瘟病菌的诱导表达。目前，尚未有报道描述OsCIPK5基因与抗病性相关，本研究结果将为水稻抗稻瘟病的信号通路中CBL-CIPK信号通路的研究提供新的思路。将抗病基因导入改良品种以减少因稻瘟病造成的水稻产量的损失，是目前认为的最经济便捷的方法。但是，当培育品种的抗病性是基于单一的R基因时，它可迅速被新出现的相容性病原菌株崩溃[48]。因此，Ca2+介导的CBL-CIPK信号通路将大大提高作物的广谱抗性。
-
![]()
图 1 水稻CIPK基因染色体分布及染色体间关系
注：红色线表示直系同源的CIPK基因对。染色体数显示在每条染色体的底部
Figure 1. Chromosomal distribution and inter-chromosomal relationship of CIPKs in rice
Note: Red line indicates homology pair of CIPKs. Number of chromosomes is shown at bottom of each chromosome.
![]()
图 2 水稻CIPK基因的系统发育分类、基因结构及保守蛋白基序的特征
注：① a为系统发育树，b为蛋白基序结构特征，c为基因结构；② OsCIPK系统发育树分为五组：A、B、C、D和E，用不同的颜色标出；③b图标尺为氨基酸长度，c图标尺为核苷酸长度；④绿色方框表示UTR，黄色方框表示CDS，黑线表示内含子，其他颜色方框代表保守结构域。
Figure 2. Phylogenetic relationships, gene structure and architecture of conserved protein motifs in CIPKs from rice
Note: ① a, phylogenetic tree; b, protein motif structure; and c, gene structure. ② OsCIPK phylogenetic tree was divided into 5 groups -A, B, C, D, and E- shown in different colors. ③ The ruler of b is for amino acid length; and the ruler of c is for nucleotide length. ④ Green box shows UTR; yellow box shows CDs; black line shows introns; and boxes of other colors show conservative domains.
![]()
图 3 水稻和拟南芥中CIPK蛋白的系统发育分析
Figure 3. Phylogenetic trees constructed with 26 CIPKs from A. thaliana and 31 CIPKs from rice
![]()
图 4 OsCIPK5在稻瘟病菌侵染后不同时间段的表达情况
Figure 4. Expressions of OsCIPK5 at times after M. grisea induction
表 1 OsCIPK5和Actin的引物序列
Table 1 The primer sequence of OsCIPK5 and Actin
基因 Genes q-F q-R OsCIPK5 CAGAGCGTCGCCATCAAGGTC GTGACAGAAGTCCACAGCCCCTATC Actin CCTCGTCTGCGATAATGGAACTG CCCTGGGCGCATCATCTC 
下载: 导出CSV
表 2 水稻基因组CIPK家族的特征
Table 2 Characteristics of CIPK family in rice genome
基因名称
Name of gene转录本ID(RAP-DB)
Transcriptome染色体位置
Locus of chromosome氨基酸长度
Length of amino acid等电点(PI)
Isoelectric point分子量(WM)
Weight of molecular/kDaOsCIPK1 Os01t0292200-01 chr1:10622951-10627532 462 6.26 52.20 OsCIPK2 Os07t0678600-01 chr7:28726817-28729527 444 9.18 50.27 OsCIPK3 Os07t0687000-01 chr7:29191197-29194187 446 6.9 50.95 OsCIPK5 Os01t0206700-03 chr1:5809589-5811389 462 9.28 51.97 OsCIPK6 Os08t0441100-01 chr8:21467453-21468593 305 5.44 32.48 OsCIPK7 Os03t0634400-01 chr3:24226372-24227930 448 9.35 48.37 OsCIPK8 Os01t0536000-01 chr1:19469262-19477590 447 6.45 50.60 OsCIPK9 Os03t0126800-00 chr3:1515355-1518947 479 8.31 53.78 OsCIPK10 Os03t0339900-01 chr3:12630705-12635182 452 9.01 51..49 OsCIPK11 Os01t0824600-01 chr1:35228913-35231725 503 8.48 56.56 OsCIPK12 Os01t0759400-02 chr1:31943529-31948580 541 8.3 59.81 OsCIPK13 Os01t0206300-00 chr1:5799607-5801142 512 8.05 56.01 OsCIPK14 Os12t0113500-01 chr12:679079-682143 440 9.44 50.32 OsCIPK15 Os11t0113700-01 chr11:630585-633622 435 9.53 49.70 OsCIPK16 Os09t0418000-03 chr9:15009182-15010716 457 8.84 50.50 OsCIPK17 Os05t0136200-01 chr5:2113595-2117110 455 6.93 50.92 OsCIPK18 Os05t0332300-01 chr5:15556085-15560404 458 8.84 51.50 OsCIPK19 Os05t0514200-01 chr5:25489349-25491319 509 7.26 56.91 OsCIPK20 Os05t0208100-01 chr5:6701053-6702933 467 8.29 51.35 OsCIPK21 Os07t0637000-01 chr7:26464444-26468858 445 7.96 50.10 OsCIPK22 Os05t0334750-00 chr5:15643309-15644732 452 8.03 49.32 OsCIPK23 Os07t0150700-01 chr7:2643540-2647773 451 9.23 50.71 OsCIPK24 Os06t0606000-01 chr6:24038959-24044785 454 8.52 50.94 OsCIPK25 Os06t0543400-01 chr6:20490388-20492271 515 8.78 56.98 OsCIPK26 Os02t0161000-01 chr2:3283264-3286366 494 9.13 55.91 OsCIPK27 Os09t0418500-00 chr9:15027311-15034300 809 5.53 87.30 OsCIPK28 Os05t0476350-00 chr5:23430891-23432198 436 9.32 49.33 OsCIPK29 Os07t0678300-01 chr7:28710749-28712470 444 8.69 48.21 OsCIPK30 Os01t0759200-01 chr1:31938559-31940184 477 9.34 53.56 OsCIPK31 Os03t0319400-01 chr3:11526273-11530306 450 7.98 50.95 OsCIPK33 Os11t0134300-01 chr11:1605182-1609130 406 6.13 46.60
下载: 导出CSV
表 3 水稻基因组中的直系同源CIPK基因对的Ka/Ks比值
Table 3 Ka/Ks value of orthologous CIPK gene pairs in rice genome
直系同源CIPK基因对
The homologous CIPK gene pairKa Ks Ka/Ks OsCIPK31/OsCIPK3 0.113 65 1.115 82 0.101 85 OsCIPK15/OsCIPK14 0.004 17 0.035 81 0.116 45 OsCIPK11/OsCIPK28 0.161 62 0.880 47 0.183 55 OsCIPK5/OsCIPK20 0.221 63 1.297 74 0.170 78 OsCIPK6/OsCIPK27 0.348 63 0.529 36 0.658 59
下载: 导出CSV
-
[1] BAUM G, LONG J C, JENKINS G I, et al. Stimulation of the blue light phototropic receptor NPH1 causes a transient increase in cytosolic Ca2+ [J]. Proceedings of the National Academy of Sciences of the United States of America, 1999, 96(23): 13554−13559. DOI: 10.1073/pnas.96.23.13554
[2] EVANS N H, MCAINSH M R, HETHERINGTON A M. Calcium oscillations in higher plants [J]. Current Opinion in Plant Biology, 2001, 4(5): 415−420. DOI: 10.1016/S1369-5266(00)00194-1
[3] KNIGHT H, KNIGHT M R. Abiotic stress signalling pathways: specificity and cross-talk [J]. Trends in Plant Science, 2001, 6(6): 262−267. DOI: 10.1016/S1360-1385(01)01946-X
[4] MACROBBIE E A C. ABA activates multiple Ca2+ fluxes in stomatal guard cells, triggering vacuolar K+(Rb+) release [J]. Proceedings of the National Academy of Sciences of the United States of America, 2000, 97(22): 12361−12368. DOI: 10.1073/pnas.220417197
[5] SANDERS D, BROWNLEE C, HARPER J F. Communicating with calcium [J]. The Plant Cell, 1999, 11(4): 691. DOI: 10.1105/tpc.11.4.691
[6] CHENG S H, WILLMANN M R, CHEN H C, et al. Calcium signaling through protein kinases. The Arabidopsis calcium-dependent protein kinase gene family [J]. Plant Physiology, 2002, 129(2): 469−485. DOI: 10.1104/pp.005645
[7] HARMON A C, GRIBSKOV M, GUBRIUM E, et al. The CDPK superfamily of protein kinases [J]. New Phytologist, 2001, 151(1): 175−183. DOI: 10.1046/j.1469-8137.2001.00171.x
[8] LUAN S, RODRIGUEZCONCEPCION M, YALOVSKY S, et al. Calmodulins and calcineurin B-like proteins: calcium sensors for specific signal response coupling in plants [J]. The Plant Cell, 2002, 14: S389−400. DOI: 10.1105/tpc.001115
[9] SNEDDEN W A, FROMM H. Calmodulin as a versatile calcium signal transducer in plants [J]. New Phytologist, 2001, 151(1): 35−66. DOI: 10.1046/j.1469-8137.2001.00154.x
[10] ZIELINSKI R E. Calmodulin and calmodulin-binding proteins in plants [J]. Annual Review of Plant Physiology and Plant Molecular Biology, 1998, 49(1): 697−725. DOI: 10.1146/annurev.arplant.49.1.697
[11] LIU J, ISHITANI M, HALFTER U, et al. The Arabidopsis thaliana SOS2 gene encodes a protein kinase that is required for salt tolerance [J]. Proceedings of the National Academy of Sciences of the United States of America, 2000, 97(7): 3730−3734. DOI: 10.1073/pnas.97.7.3730
[12] ISHITANI M, LIU J P, HALFTER U, et al. SOS3 function in plant salt tolerance requires N-myristoylation and calcium binding [J]. The Plant Cell, 2000, 12(9): 1667−1677. DOI: 10.1105/tpc.12.9.1667
[13] SÁNCHEZ-BARRENA M J, MARTÍNEZ-RIPOLL M, ZHU J K, et al. The structure of the Arabidopsis thaliana SOS3: molecular mechanism of sensing calcium for salt stress response [J]. Journal of Molecular Biology, 2005, 345(5): 1253−1264. DOI: 10.1016/j.jmb.2004.11.025
[14] GONG D M, GUO Y, SCHUMAKER K S, et al. The SOS3 family of calcium sensors and SOS2 family of protein kinases in Arabidopsis [J]. Plant Physiology, 2004, 134(3): 919−926. DOI: 10.1104/pp.103.037440
[15] SHI H Z, ISHITANI M, KIM C, et al. The Arabidopsis thaliana salt tolerance gene SOS1 encodes a putative Na+/H+ antiporter [J]. Proceedings of the National Academy of Sciences of the United States of America, 2000, 97(12): 6896−6901. DOI: 10.1073/pnas.120170197
[16] QIU Q S, GUO Y, DIETRICH M A, et al. Regulation of SOS1, a plasma membrane Na+/H+ exchanger in Arabidopsis thaliana, by SOS2 and SOS3 [J]. Proceedings of the National Academy of Sciences of the United States of America, 2002, 99(12): 8436−8441. DOI: 10.1073/pnas.122224699
[17] QIU Q S, GUO Y, QUINTERO F J, et al. Regulation of vacuolar Na+/H+ exchange in Arabidopsis thaliana by the salt-overly-sensitive (SOS) pathway [J]. Journal of Biological Chemistry, 2004, 279(1): 207−215. DOI: 10.1074/jbc.M307982200
[18] XIONG L M, SCHUMAKER K S, ZHU J K. Cell signaling during cold, drought, and salt stress [J]. The Plant Cell, 2002, 14(S1): S165−183.
[19] ZHU J K. Salt and drought stress signal transduction in plants [J]. Annual Review of Plant Biology, 2002, 53(1): 247−273. DOI: 10.1146/annurev.arplant.53.091401.143329
[20] KOLUKISAOGLUÜ, WEINL S, BLAZEVIC D, et al. Calcium sensors and their interacting protein kinases: genomics of the Arabidopsis and rice CBL-CIPK signaling networks [J]. Plant Physiology, 2004, 134(1): 43−58. DOI: 10.1104/pp.103.033068
[21] XIANG Y, HUANG Y M, XIONG L Z. Characterization of stress-responsive CIPK genes in rice for stress tolerance improvement [J]. Plant Physiology, 2007, 144(3): 1416−1428. DOI: 10.1104/pp.107.101295
[22] CHEN X F, GU Z M, LIU F, et al. Molecular analysis of rice CIPKs involved in both biotic and abiotic stress responses [J]. Rice Science, 2011, 18(1): 1−9. DOI: 10.1016/S1672-6308(11)60001-2
[23] HOLUB E B. The arms race is ancient history in Arabidopsis, the wildflower [J]. Nature Reviews Genetics, 2001, 2(7): 516−527. DOI: 10.1038/35080508
[24] YU Y H, XIA X L, YIN W L, et al. Comparative genomic analysis of CIPK gene family in Arabidopsis and Populus [J]. Plant Growth Regulation, 2007, 52(2): 101−110. DOI: 10.1007/s10725-007-9165-3
[25] KANWAR P, SANYAL S K, TOKAS I, et al. Comprehensive structural, interaction and expression analysis of CBL and CIPK complement during abiotic stresses and development in rice [J]. Cell Calcium, 2014, 56(2): 81−95. DOI: 10.1016/j.ceca.2014.05.003
[26] ROY S W, PENNY D. Patterns of intron loss and gain in plants: intron loss–dominated evolution and genome-wide comparison of O. sativa and A. thaliana [J]. Molecular Biology and Evolution, 2007, 24(1): 171−181.
[27] NURUZZAMAN M, MANIMEKALAI R, SHARONI A M, et al. Genome-wide analysis of NAC transcription factor family in rice [J]. Gene, 2010, 465(1/2): 30−44.
[28] KLEIST T J, SPENCLEY A L, LUAN S. Comparative phylogenomics of the CBL-CIPK calcium-decoding network in the moss Physcomitrella, Arabidopsis, and other green lineages [J]. Frontiers in Plant Science, 2014(5): 187.
[29] ZHANG H F, YANG B, LIU W Z, et al. Identification and characterization of CBL and CIPK gene families in canola (Brassica napus L.) [J]. BMC Plant Biology, 2014, 14(1): 8. DOI: 10.1186/1471-2229-14-8
[30] KUDLA J, BATISTIČ O, HASHIMOTO K. Calcium signals: the lead currency of plant information processing [J]. The Plant Cell, 2010, 22(3): 541−563. DOI: 10.1105/tpc.109.072686
[31] HU H C, WANG Y Y, TSAY Y F. AtCIPK8, a CBL‐interacting protein kinase, regulates the low-affinity phase of the primary nitrate response [J]. Plant Journal, 2009, 57(2): 264−278. DOI: 10.1111/j.1365-313X.2008.03685.x
[32] CHUNG E, PARK J M, OH S K, et al. Molecular and biochemical characterization of the Capsicum annuum calcium-dependent protein kinase 3 (CaCDPK3) gene induced by abiotic and biotic stresses [J]. Planta, 2004, 220(2): 286−295. DOI: 10.1007/s00425-004-1372-9
[33] LUDWIG A A, SAITOH H, FELIX G, et al. Ethylene-mediated cross-talk between calcium-dependent protein kinase and MAPK signaling controls stress responses in plants [J]. Proceedings of the National Academy of Sciences of the United States of America, 2005, 102(30): 10736−10741. DOI: 10.1073/pnas.0502954102
[34] FROHNMEYER H, LOYALL L, BLATT M R, et al. Millisecond UV-B irradiation evokes prolonged elevation of cytosolic-free Ca2+ and stimulates gene expression in transgenic parsley cell cultures [J]. Plant Journal, 1999, 20(1): 109−117. DOI: 10.1046/j.1365-313X.1999.00584.x
[35] TREWAVAS A J, MALHÓ R. Ca2+ signalling in plant cells: the big network! [J]. Current Opinion in Plant Biology, 1998, 1(5): 428−443. DOI: 10.1016/S1369-5266(98)80268-9
[36] TRAN P O T, HINMAN L E, UNGER G M, et al. A wound-induced[Ca2+]i increase and its transcriptional activation of immediate early genes is important in the regulation of motility [J]. Experimental Cell Research, 1999, 246(2): 319−326. DOI: 10.1006/excr.1998.4239
[37] YANG W Q, KONG Z S, OMO-IKERODAH E, et al. Calcineurin B-like interacting protein kinase OsCIPK23 functions in pollination and drought stress responses in rice (Oryza sativa L.) [J]. Journal of Genetics and Genomics, 2008, 35(9): 531−543. DOI: 10.1016/S1673-8527(08)60073-9
[38] OHBA H, STEWARD N, KAWASAKI S, et al. Diverse response of rice and maize genes encoding homologs of WPK4, an SNF1-related protein kinase from wheat, to light, nutrients, low temperature and cytokinins [J]. Molecular and General Genetics, 2000, 263(2): 359−366. DOI: 10.1007/s004380051179
[39] KIM K N, LEE J S, HAN H E, et al. Isolation and characterization of a novel rice Ca2+-regulated protein kinase gene involved in responses to diverse signals including cold, light, cytokinins, sugars and salts [J]. Plant Molecular Biology, 2003, 52(6): 1191−1202. DOI: 10.1023/B:PLAN.0000004330.62660.a2
[40] LEE K W, CHEN P W, LU C A, et al. Coordinated responses to oxygen and sugar deficiency allow rice seedlings to tolerate flooding [J]. Science Signaling, 2009, 2(91): ra61.
[41] HUANG F, LIAN L, HE W, et al. Genome-wide profiling of changes in gene expression in response to infection of the japonica rice variety Yunyin by Magnaporthe oryzae [J]. Molecular Breeding, 2014, 34(4): 1965−1974. DOI: 10.1007/s11032-014-0155-7
[42] DENG X M, HU W, WEI S Y, et al. TaCIPK29, a CBL-interacting protein kinase gene from wheat, confers salt stress tolerance in transgenic tobacco [J]. PLoS One, 2013, 8(7): e69881. DOI: 10.1371/journal.pone.0069881
[43] CANNON S B, MITRA A, BAUMGARTEN A, et al. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana [J]. BMC Plant Biology, 2004, 4(1): 10. DOI: 10.1186/1471-2229-4-10
[44] ZHU J K, LIU J P, XIONG L M. Genetic analysis of salt tolerance in Arabidopsis: evidence for a critical role of potassium nutrition [J]. The Plant Cell, 1998, 10(7): 1181−1191. DOI: 10.1105/tpc.10.7.1181
[45] LI L, KIM B G, CHEONG Y H, et al. A Ca2+ signaling pathway regulates a K+ channel for low-K response in Arabidopsis [J]. Proceedings of the National Academy of Sciences of the United States of America, 2006, 103(33): 12625−12630. DOI: 10.1073/pnas.0605129103
[46] XU J, LI H D, CHEN L Q, et al. A protein kinase, interacting with two calcineurin B-like proteins, regulates K+ transporter AKT1 in Arabidopsis [J]. Cell, 2006, 125(7): 1347−1360. DOI: 10.1016/j.cell.2006.06.011
[47] MARTÍNEZ -ATIENZA J, JIANG X Y, GARCIADEBLAS B, et al. Conservation of the salt overly sensitive pathway in rice [J]. Plant Physiology, 2007, 143(2): 1001−1012. DOI: 10.1104/pp.106.092635
[48] HAYASHI K, YOSHIDA H, ASHIKAWA I. Development of PCR-based allele-specific and InDel marker sets for nine rice blast resistance genes [J]. Theoretical and Applied Genetics, 2006, 113(2): 251−260. DOI: 10.1007/s00122-006-0290-6
-
期刊类型引用(5)
1. 高姗姗，冯浩，黄鑫，张鑫月，张泽阳，刘贵河. 大麦CIPK基因家族鉴定及表达分析. 现代农业科技. 2024(04): 156-163 . 
百度学术
2. 杨宏羽，李瑞瑞，何雷，余斌，师桂英. 马铃薯StCIPK基因家族鉴定及表达特征分析. 植物生理学报. 2022(10): 1919-1934 . 百度学术
3. 张和桥，廖元德，韦鸿，卢晓惠. 关于优质杂交稻的育种创新探讨. 农村实用技术. 2021(04): 64-65 . 百度学术
4. 王晓宇，刘栩铭，韩孟良，李敏，吴颖，靳亚楠，张继星. 蓖麻RcCIPKs基因家族鉴定与冷胁迫下的表达模式分析. 农业生物技术学报. 2021(07): 1260-1273 . 百度学术
5. 林悦龙，王昱澎，朱永生，郑燕梅，蔡秋华，何炜，陈丽萍，谢鸿光，魏毅东，谢华安，张建福. 粳稻云引稻瘟病抗性近等基因系与免疫应答相关的LRR-RLKs差异基因分析. 科学通报. 2021(30): 3873-3885 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 3760
- HTML全文浏览量: 1034
- PDF下载量: 144
- 被引次数: 6
