• 中国中文核心期刊
  • 中国科学引文数据库(CSCD)核心库来源期刊
  • 中国科技论文统计源期刊(CJCR)
  • 第二届国家期刊奖提名奖

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

两种兰科植物CaMCML基因家族全基因组分析

杨前宇 王涛 梁立雄 李潞滨 刘蕾

引用本文:
Citation:

两种兰科植物CaMCML基因家族全基因组分析

    作者简介: 杨前宇(1985-), 男, 博士研究生, 研究方向为园林植物与观赏园艺.E-mail:yangdayu0055@163.com.
    通讯作者: 刘蕾, liulei0828@126.com
  • 基金项目:

    北京市植物资源研究开发重点实验室开放课题 PRRD-2016-ZD1

    国家高技术研究发展计划(863计划)子课题 2013AA102607-3

  • 中图分类号: S718.46

Genome-wide Analysis of CaM/CML Gene Family in Two Orchidaceae Species

    Corresponding author: LIU Lei, liulei0828@126.com
  • CLC number: S718.46

  • 摘要: 目的 探讨CaMCML基因家族在兰科植物不同组织和与菌根真菌建立共生关系过程中的潜在功能。 方法 依据已发表的拟南芥CaMCML基因家族蛋白序列,利用生物信息学方法对小兰屿蝴蝶兰和铁皮石斛全基因组进行搜索,通过多种生物软件或在线工具对CaMCML基因家族蛋白序列进行筛选及确定、蛋白结构分析、系统进化分析及结构域预测;利用转录组数据绘制heatmap图,对小兰屿蝴蝶兰不同组织(花、叶、茎、根)及石斛种子与美胞胶膜菌共生条件下CaMCML基因家族的表达情况进行分析。 结果 在小兰屿蝴蝶兰和铁皮石斛全基因组中均预测到4个CaM蛋白和54个CML蛋白;小兰屿蝴蝶兰CaMCML家族基因中39个基因没有内含子,19个基因具有内含子;铁皮石斛CaMCML基因家族中有41个基因没有内含子,17个基因具有内含子;系统进化分析将小兰屿蝴蝶兰和铁皮石斛CaM和CML蛋白家族各分为10个亚家族;小兰屿蝴蝶兰中9个基因在叶中的表达相对于花、茎、根为上调表达,2个基因在叶中的表达相对于花、茎、根却为下调表达;3个基因在花中的表达相对于叶、茎、根为上调表达;3个基因在花和叶中的表达相对于茎和根为上调表达;铁皮石斛中4个基因在与美胞胶膜菌共生萌发文库中的表达相对于非共生文库为上调表达;4个基因在共生萌发文库中的表达相对于非共生文库为下调表达。 结论 CaMCML家族基因在兰科植物不同组织具有重要作用,且参与铁皮石斛种子与美孢胶膜菌共生萌发的生物学调控过程。
  • 图 1  兰科植物与拟南芥CaM蛋白序列比对

    Figure 1.  The protein sequence alignment of OPeCaMs, ODcCaMs and AtCaM1

    图 2  小兰屿蝴蝶兰CaM及CML家族系统进化树

    Figure 2.  The phylogenetic tree of OPeCaM and OPeCML

    图 3  铁皮石斛CaM及CML家族系统进化树

    Figure 3.  The phylogenetic tree of ODcCaM and ODcCML

    图 4  小兰屿蝴蝶兰CaMCML基因在4个组织部位表达水平热图

    Figure 4.  Heatmap showing the expression levels of OPeCaM and OPeCML genes in four tissues

    图 5  CaMCML基因在共生和非共生萌发的铁皮石斛种子中的表达水平热图

    Figure 5.  Heatmap showing the expression levels of ODcCaM and ODcCML genes in symbiotic and asymbiotic germinated seeds

    表 1  小兰屿蝴蝶兰OPeCaM及OPeCML蛋白序列特征

    Table 1.  The characteristics of OPeCaM/OPeCML protein sequence in P.equestris

    名称
    Name
    基因ID
    Gene ID
    亚家族
    subfamily
    蛋白序列长度
    Length of protein/aa
    EF hands 分子量
    Mol Wt/kDa
    等电点
    pI
    内含子
    Intron
    甲硫氨酸的含量
    percentage methionine/%
    OPeCaM1 PEQU_18438 1 149 4 16.96 3.89 2 6.02
    OPeCaM2 PEQU_10717 1 149 4 16.96 3.89 2 6.01
    OPeCaM3 PEQU_12537 1 148 4 16.96 3.88 2 6.00
    OPeCaM4 PEQU_26335 1 149 4 17.96 3.95 1 6.05
    OPeCML01 PEQU_04925 1 148 4 12.96 3.63 4 4.82
    OPeCML02 PEQU_07491 10 155 4 16.96 3.95 1 5.01
    OPeCML03 PEQU_06469 10 166 4 16.96 4.02 0 5.90
    OPeCML04 PEQU_16519 1 135 3 20.96 3.80 0 4.11
    OPeCML05 PEQU_34003 10 199 4 30.98 4.21 6 4.67
    OPeCML06 PEQU_29891 1 120 3 14.96 3.52 2 2.83
    OPeCML07 PEQU_25947 7 190 4 24.95 4.16 1 6.36
    OPeCML08 PEQU_22407 6 155 4 19.97 3.98 0 3.42
    OPeCML09 PEQU_40659 5 154 4 18.95 3.95 0 3.71
    OPeCML10 PEQU_35611 1 72 2 8.96 3.73 1 2.44
    OPeCML11 PEQU_15751 10 141 3 17.96 3.83 4 2.12
    OPeCML12 PEQU_06469 10 166 4 16.97 4.19 0 6.23
    OPeCML13 PEQU_41809 6 219 4 24.95 4.38 0 5.22
    OPeCML14 PEQU_30031 8 137 4 20.96 3.88 0 5.01
    OPeCML15 PEQU_40424 2 193 4 23.97 4.07 0 4.93
    OPeCML16 PEQU_18134 6 199 4 25.97 4.18 0 5.22
    OPeCML17 PEQU_24120 8 134 4 18.96 3.80 0 5.33
    OPeCML18 PEQU_32857 8 140 4 17.95 3.81 0 4.31
    OPeCML19 PEQU_32861 8 141 4 16.97 3.82 0 5.46
    OPeCML20 PEQU_11970 8 140 4 17.96 3.87 0 6.82
    OPeCML21 PEQU_32860 8 141 4 15.96 3.87 0 5.11
    OPeCML22 PEQU_11969 8 237 5 25.95 4.28 1 6.30
    OPeCML23 PEQU_25577 8 143 3 18.95 4.38 0 4.41
    OPeCML24 PEQU_28290 8 143 4 15.96 3.77 0 5.02
    OPeCML25 PEQU_05240 3 143 3 16.97 3.88 0 4.92
    OPeCML26 PEQU_19202 5 160 3 18.97 3.96 0 3.83
    OPeCML27 PEQU_17666 5 197 3 22.95 4.19 0 3.77
    OPeCML28 PEQU_35797 6 110 3 11.96 3.81 0 3.61
    OPeCML29 PEQU_11304 8 152 3 19.96 3.85 0 4.16
    OPeCML30 PEQU_20039 6 114 3 15.96 3.77 0 3.03
    OPeCML31 PEQU_28291 8 110 3 11.96 3.78 0 4.23
    OPeCML32 PEQU_09901 8 111 3 16.96 3.79 0 3.19
    OPeCML33 PEQU_29894 2 101 3 6.95 3.78 0 2.68
    OPeCML34 PEQU_32859 8 115 3 14.96 3.78 0 3.46
    OPeCML35 PEQU_09640 9 183 4 38.96 3.87 0 3.10
    OPeCML36 PEQU_15088 4 166 2 23.95 3.85 0 2.89
    OPeCML37 PEQU_28312 4 176 2 20.97 3.83 0 2.54
    OPeCML38 PEQU_00748 7 168 3 7.96 3.85 0 3.88
    OPeCML39 PEQU_39569 10 108 2 14.95 3.81 0 3.16
    OPeCML40 PEQU_40143 10 127 2 18.96 3.77 0 1.93
    OPeCML41 PEQU_09077 4 103 2 13.96 3.88 0 2.42
    OPeCML42 PEQU_09518 4 158 2 21.96 3.89 0 2.02
    OPeCML43 PEQU_24484 10 225 4 34.95 4.28 3 5.33
    OPeCML44 PEQU_00516 7 106 2 7.97 3.77 0 2.11
    OPeCML45 PEQU_11010 1 304 2 30.96 4.75 4 1.66
    OPeCML46 PEQU_04326 1 253 2 23.96 4.56 0 3.35
    OPeCML47 PEQU_41809 6 219 4 24.96 4.23 5 4.26
    OPeCML48 PEQU_17666 5 197 3 22.97 4.06 0 3.28
    OPeCML49 PEQU_01306 10 153 3 18.95 3.99 1 3.66
    OPeCML50 PEQU_12715 10 147 3 17.97 4.20 2 2.82
    OPeCML51 PEQU_17236 10 211 4 17.96 4.12 1 4.21
    OPeCML52 PEQU_00999 1 113 4 13.96 4.32 0 3.16
    OPeCML53 PEQU_12714 10 153 2 19.95 4.41 1 3.19
    OPeCML54 PEQU_01183 1 121 3 13.96 3.52 0 2.47
    下载: 导出CSV

    表 2  铁皮石斛ODcCaM及ODcCML蛋白序列特征

    Table 2.  The characteristics of ODcCaM/ODcCMLprotein sequence in D. officinale

    名称
    Name
    基因ID
    Gene ID
    亚家族
    subfamily
    蛋白序列长度
    Length of protein/aa
    EF hands 分子量
    Mol Wt/kDa
    等电点
    pI
    内含子
    Intron
    甲硫氨酸的含量
    percentage methionine /%
    ODcCaM1 DG_10048944 1 149 4 16.96 3.89 1 6.01
    ODcCaM2 DG_10042177 1 149 4 16.96 3.88 1 6.03
    ODcCaM3 DG_10042221 1 149 4 16.96 3.88 3 6.00
    ODcCaM4 DG_10069718 1 148 4 11.96 3.67 3 6.02
    ODcCML01 DG_10099596 8 147 4 16.96 3.63 3 5.11
    ODcCML02 DG_10026132 6 154 4 22.96 4.20 0 3.96
    ODcCML03 DG_10032884 6 152 4 18.96 4.12 0 4.88
    ODcCML04 DG_10102530 6 154 4 18.97 4.32 0 4.46
    ODcCML05 DG_10001096 6 165 4 22.96 4.41 0 5.52
    ODcCML06 DG_10086740 6 142 4 12.95 3.52 0 4.26
    ODcCML07 DG_10112859 6 154 4 21.96 4.16 0 4.14
    ODcCML08 DG_10069652 6 140 4 20.15 3.78 0 3.82
    ODcCML09 DG_10048052 6 143 4 16.96 3.75 0 6.26
    ODcCML10 DG_10013983 6 145 4 12.96 3.73 0 3.93
    ODcCML11 DG_10087227 6 145 4 11.96 3.73 0 3.74
    ODcCML12 DG_10120493 10 113 3 13.95 3.79 0 2.43
    ODcCML13 DG_10127394 2 112 3 13.96 3.78 0 2.16
    ODcCML14 DG_10045310 2 177 3 17.97 4.88 1 6.29
    ODcCML15 DG_10099614 9 124 3 7.96 3.77 2 5.26
    ODcCML16 DG_10006683 10 86 2 9.98 3.88 1 1.19
    ODcCML17 DG_10099594 9 130 3 14.96 3.89 2 3.95
    ODcCML18 DG_10012949 8 107 2 11.97 3.80 1 2.27
    ODcCML19 DG_10018482 6 110 3 13.96 3.82 0 3.36
    ODcCML20 DG_10068885 6 206 4 25.96 4.67 0 4.42
    ODcCML21 DG_10080246 6 206 4 25.96 4.67 0 5.16
    ODcCML22 DG_10107219 7 195 4 27.95 4.28 0 4.82
    ODcCML23 DG_10115967 6 219 4 24.96 4.38 0 5.56
    ODcCML24 DG_10045416 7 165 4 25.96 3.77 0 6.01
    ODcCML25 DG_10131684 7 151 2 21.95 3.88 0 2.44
    ODcCML26 DG_10013708 8 87 2 6.98 3.72 1 2.22
    ODcCML27 DG_10086551 8 147 3 21.96 3.89 0 4.01
    ODcCML28 DG_10057311 4 141 3 10.97 3.85 0 2.82
    ODcCML29 DG_10047181 7 104 2 10.96 3.78 0 1.71
    ODcCML30 DG_10112504 3 231 3 39.95 4.67 0 3.29
    ODcCML31 DG_10114938 8 303 3 37.97 5.88 6 4.15
    ODcCML32 DG_10061210 5 107 3 12.96 3.79 3 3.23
    ODcCML33 DG_10067799 7 100 3 7.96 3.78 0 4.53
    ODcCML34 DG_10089833 7 195 4 28.97 4.58 0 5.14
    ODcCML35 DG_10069653 6 116 3 18.95 3.87 0 2.96
    ODcCML36 DG_10112846 6 133 3 21.96 3.85 0 3.54
    ODcCML37 DG_10112841 6 128 3 20.96 3.83 0 3.26
    ODcCML38 DG_10122957 6 133 3 21.95 3.85 0 2.80
    ODcCML39 DG_10012468 7 190 4 24.97 4.58 0 2.51
    ODcCML40 DG_10113507 6 185 4 25.97 4.17 0 3.82
    ODcCML41 DG_10112861 6 337 4 48.95 5.88 2 3.99
    ODcCML42 DG_10140024 6 148 2 16.96 3.89 1 1.58
    ODcCML43 DG_10008482 6 184 3 30.96 4.28 0 3.44
    ODcCML44 DG_10088170 6 175 3 31.97 4.02 0 2.15
    ODcCML45 DG_10078503 7 189 4 17.96 4.27 0 5.52
    ODcCML46 DG_10098264 7 189 4 17.95 4.27 0 4.32
    ODcCML47 DG_10078224 7 195 3 24.96 4.99 0 3.66
    ODcCML48 DG_10002690 7 107 2 11.95 3.72 0 2.32
    ]ODcCML49 DG_10049524 7 176 2 23.93 4.25 0 1.20
    ODcCML50 DG_10045416 7 165 2 25.97 4.17 0 2.21
    ODcCML51 DG_10024582 7 181 2 24.96 4.28 0 3.03
    ODcCML52 DG_10116046 7 327 2 29.97 5.88 4 2.52
    ODcCML53 DG_10064998 7 691 2 76.94 8.85 10 2.23
    ODcCML54 DG_10098265 5 191 4 25.95 4.58 0 4.19
    下载: 导出CSV

    表 3  小兰屿蝴蝶兰CaM蛋白序列比对(OPeCaMs)

    Table 3.  Alignment of protein sequence of OPeCaMs

    OPeCaM1 OPeCaM2 OPeCaM3 OPeCaM4
    OPeCaM1 6 9 12
    OPeCaM2 95.97 10 11
    OPeCaM3 93.96 93.29 12
    OPeCaM4 91.95 92.62 91.95
      注:数字分别表示相似度(100%)和氨基酸差异数。
      Note: The numbersaresimilarity(100%) and quantity difference amino acids.
    下载: 导出CSV

    表 4  铁皮石斛CaM蛋白序列比对(ODcCaMs)

    Table 4.  Alignment of protein sequence of ODcCaMs

    ODcCaM1 ODcCaM2 ODcCaM3 ODcCaM4
    ODcCaM1 3 3 12
    ODcCaM2 97.99 2 11
    ODcCaM3 97.99 98.66 11
    ODcCaM4 92.00 92.67 92.67
      注:数字分别表示相似度(100%)和氨基酸差异数。
      Note: The numbersaresimilarity(100%) and quantity difference amino acids.
    下载: 导出CSV
  • [1]

    Cribb P J, Keu S P, Dixon K W, et al. Ochid Conservation[J]. Sabah:Natural History Publications, 2003, 4(1):21-30.
    [2]

    Dressler R L. Phylogeny and classification of the orchid family[J]. Cambridge:University Press, 1993, 7(2):55-62.
    [3] 郎楷永.兰科植物区系中一些有意义属的地理分布格局的研究[J].植物分类学报, 1994, 32(4):328-339.

    [4] 张殷波, 杜昊东, 金效华, 等.中国野生兰科植物物种多样性与地理分布[J].科学通报, 2015, 60(2):179-188.

    [5]

    Rasmussen H N. Terrestrial orchids from seed to mycotrophic plant[J]. Cambridge University Press, Cambridge, 1995, 8(1):68-75.
    [6] 陈心启, 吉占和.中国兰花全书.第二版[M].北京:中国林业出版社, 2003.

    [7]

    Harley J L. The significance of mycorrhiza[J]. Mycological Research, 1989, 92(2):129-139. doi: 10.1016/S0953-7562(89)80001-2
    [8]

    Cai J, Liu X, Vanneste K, et al. The genome sequence of the orchid Phalaenopsis equestris[J]. Nature Genetics, 2015, 47(1):65-72. doi: 10.1038/ng.3149
    [9]

    Zhang G Q, Xu Q, Bian C, et al. The Dendrobium catenatum Lindl. genome sequence provides insights into polysaccharide synthase, floral development and adaptive evolution[J]. Scientific reports, 2016, 6:19029-19038. doi: 10.1038/srep19029
    [10]

    Yan L, Wang X, Liu H, et al. The Genome of Dendrobium officinale Illuminates the Biology of the Important Traditional Chinese Orchid Herb[J].Molecular Plant.2015, 8(6):922-934. doi: 10.1016/j.molp.2014.12.011
    [11]

    Defalco T A, Bender K W, Snedden W A. Breaking the code:Ca2+ sensors in plant signalling[J]. Biochemical Journal, 2010, 425(1):27-40. doi: 10.1042/BJ20091147
    [12] 梁秋芬, 刘宽灿, 徐碧玉, 等.钙调蛋白在植物发育中的功能[J].生命科学研究, 2005, 9(S2):1-5.

    [13]

    Popescu S C, Popescu G V, Bachan S, et al. Differential binding of calmodulin-related proteins to their targets revealed through high-density Arabidopsis protein microarrays[J]. Procee-dings of the National Academy of Sciences of the United States of America, 2007, 104(11):4730-4735. doi: 10.1073/pnas.0611615104
    [14]

    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
    [15] 毛国红, 宋林霞, 孙大业.植物钙调素结合蛋白研究进展[J].植物生理与分子生物学学报, 2004, 30(5):481-488.

    [16] 张成, 王喆之.豇豆钙调蛋白cDNA的克隆及序列分析[J].陕西师范大学学报, 2006, 34(2):88-91.

    [17]

    Larkin M A, Blackshields G, Brown N P, et al. Clustal W and clustal X version 2.0[J]. Bioinformatics, 2007, 23(21):2947-2948. doi: 10.1093/bioinformatics/btm404
    [18]

    Bologna G, Yvon C, Duvaud S, et al. N-Terminal myristoylation predictions by ensembles of neural networks[J]. Proteomics, 2004, 4(6):1626-1632. doi: 10.1002/(ISSN)1615-9861
    [19]

    Saitou N, Nei M. The neighbor-joining method:a new method for reconstructing phylogenetic trees[J]. Molecular Biology and Evolution, 1987, 4(4):406-425.
    [20]

    Hu B, Jin J, Guo A, et al. GSDS 2.0:an upgraded gene feature visualization server[J]. Bioinformatics, 2015, 31(8):1296-1297. doi: 10.1093/bioinformatics/btu817
    [21] 宋争.石斛与美胞胶膜菌共生关系研究[D].北京: 中国林业科学研究院.2017.

    [22]

    Deng W, Wang Y, Liu Z, et al. HemI:a toolkit for illustrating heatmaps[J]. PLOS ONE, 2014, 9(11):e111988. doi: 10.1371/journal.pone.0111988
    [23]

    Zielinski R E. Calmodulin and calmodulin-binding proteins in plants[J]. Annual Review of Plant Biology, 1998, 49(1):697-725. doi: 10.1146/annurev.arplant.49.1.697
    [24] 张春梅, 谢晓蓉, 刘金荣, 等.Ca2+和钙调素拮抗剂W7对PEG胁迫下不同耐旱型紫花苜蓿抗氧化系统的影响[J].草地学报, 2012, 20(6):1072-1076.

    [25]

    Mccormack E, Braam J. Calmodulins and related potential Calcium sensors of Arabidopsis[J]. New Phytologist, 2003, 159(3):585-598. doi: 10.1046/j.1469-8137.2003.00845.x
    [26]

    Boonburapong B, Buaboocha T. Genome-wide identification and analyses of the rice calmodulin and related potential Calcium sensor proteins[J]. BMC Plant Biology, 2007, 7:4-5. doi: 10.1186/1471-2229-7-4
    [27] 王雪霁.杨树CaMCML基因家族生物信息学分析[D].北京: 中国林业科学研究院, 2016.

    [28]

    Zhao Y, Liu W, Xu Y, et al. Genome-wide identification and functional analyses of calmodulin genes in Solanaceous species[J]. BMC Plant Biology, 2013, 13(1):70.
    [29]

    Mccormack E, Tsai Y, Braam J. Handling Calcium signaling:Arabidopsis CaMs and CMLs[J]. Trends in Plant Science, 2005, 10(8):383-389. doi: 10.1016/j.tplants.2005.07.001
    [30]

    Zielinski R E. Calmodulin and calmodulin-binding proteins in plants[J]. Annual Review of Plant Biology, 1998, 49(1):697-725. doi: 10.1146/annurev.arplant.49.1.697
    [31] 庞俊峰, 于卓, 张占路, 等.木榄CaM基因的克隆及序列分析[J].生物技术通报, 2010(12):127-131.

    [32]

    Luan S, Kudla J, Rodriguezconcepcion M, et al. Calmodulins and calcineurin B-like proteins Calcium sensors for specific signal response coupling in plants[J]. The Plant Cell, 2002, 14(Suppl):S389-S400.
    [33] 孔海燕, 贾桂霞, 温跃戈.钙在植物花发育过程中的作用[J].植物学通报, 2003, 20(2):168-177. doi: 10.3969/j.issn.1674-3466.2003.02.006

    [34]

    Lemaire F, Sintes G, Morel P. Mineral needs of Begonia×elatior during the growing period and the flowering time[J]. Hydroponics and Transplant Production, 1995(396):219-226.
    [35] 褚孟螈, 黄金城.梅树花芽形态分化及其物质代谢的研究[J].北京林业大学学报, 1995, 17(1):68-74.

    [36]

    Galon Y, Finkler A, Fromm H. Calcium-regulated transcription in plants[J]. Molecular Plant, 2010, 3(4):653-669. doi: 10.1093/mp/ssq019
    [37]

    Yang T, Poovaiah B W. Hydrogen peroxide homeostasis:Activation of plant catalase by calcium/calmodulin[J]. Proceedings of the National Academy of Sciences of the United States of America, 2002, 275(49):38467-38473.
    [38]

    Heo WD, Lee SH, et al. Involvement of specific calmodulin isoforms in salicylic acid-independent activation of plant disease resistance responses[J]. Proceedings of the National Academy of Sciences of the United States of America, 1999, 96(2):766-771. doi: 10.1073/pnas.96.2.766
    [39]

    Cao J, Li X, Lv Y, et al. Comparative analysis of the phytocyanin gene family in 10 plant species a focus on Zea mays[J]. Frontiers in Plant Science, 2015, 6:515-527.
    [40]

    Yamakawa H, Mitsuhara I, Ito N, et al. Transcriptionally and post-transcriptionally regulated response of 13 calmodulin genes to tobacco Mosaic virus-induced cell death and wounding in tobacco plan[J]t. European Journal of Biochemistry, 2001, 268(14):3916-3929. doi: 10.1046/j.1432-1327.2001.02301.x
    [41] 徐亮, 刘振宇, 吕全, 等.马尾松钙调素基因的克隆及响应松材线虫侵染的表达特征[J].林业科学研究, 2014, 27(3):323-328.

    [42]

    Leba LJ, Cheval C, et al. CML9 an Arabidopsis calmodulin-like protein, contributes to plant innate immunity through a flagellin-dependent signalling pathway[J]. Plant journal, 2014, 10(2):19-25.
    [43]

    Harrison M J. Development of the arbuscular mycorrhizal symbiosis[J]. Current Opinion in Plant Biology, 1998, 1(4):360-365. doi: 10.1016/1369-5266(88)80060-8
  • [1] 宋争李潞滨梁立雄王涛 . 铁皮石斛phytocyanin基因家族全基因组分析. 林业科学研究, 2018, 31(2): 98-106. doi: 10.13275/j.cnki.lykxyj.2018.02.014
    [2] 王雪霁梁立雄李潞滨王涛 . 小兰屿蝴蝶兰R2R3-MYB转录因子分析. 林业科学研究, 2018, 31(3): 104-113. doi: 10.13275/j.cnki.lykxyj.2018.03.014
    [3] 郑宝强乔红娟李柏君KAOTachung张燕王雁 . 铁皮石斛(Dendrobium catenatum)叶艺形成机制的初探. 林业科学研究, 2020, 33(4): 177-183. doi: 10.13275/j.cnki.lykxyj.2020.04.022
    [4] 吕中睿刘宏张国昀于立洋罗红梅何彩云 . 沙棘UGT基因家族的全基因组鉴定与表达分析. 林业科学研究, 2021, 34(6): 9-19. doi: 10.13275/j.cnki.lykxyj.2021.06.002
    [5] 刘聪张洋夏德安陈雪冰魏志刚 . 毛果杨PLD基因家族全基因组水平鉴定及其盐胁迫下的表达分析. 林业科学研究, 2021, 34(3): 23-36. doi: 10.13275/j.cnki.lykxyj.2021.03.003
    [6] 陈鸿鹏朱凤云吴志华谢耀坚 . 赤桉GAPDH家族基因的克隆及其序列分析. 林业科学研究, 2014, 27(1): 120-127.
    [7] 张进李建波刘伯斌陈军卢孟柱 . 杨树CDPK基因家族的表达分析及功能预测. 林业科学研究, 2014, 27(5): 604-611.
    [8] 杨海昕刘晓莹詹亚光范桂枝 . 白桦BpAMT基因家族鉴定及表达模式分析. 林业科学研究, 2023, 36(2): 133-143. doi: 10.12403/j.1001-1498.20220242
    [9] 李煜张进王丽娟卢孟柱 . 杨树PtROP家族基因的表达分析与功能预测. 林业科学研究, 2017, 30(1): 1-9. doi: 10.13275/j.cnki.lykxyj.2017.01.001
    [10] 孙化雨娄永峰李利超赵韩生高志民 . 毛竹TIPs基因家族成员组织表达模式研究. 林业科学研究, 2016, 29(4): 521-528.
    [11] 胡梦璇宋学勤刘颖丽赵树堂 . 杨树MYC基因家族成员表达模式研究. 林业科学研究, 2023, 36(3): 32-40. doi: 10.12403/j.1001-1498.20220503
    [12] 徐向东任逸秋张利李煜王丽娟卢孟柱 . 杨树PIF基因家族成员表达模式研究. 林业科学研究, 2018, 31(2): 19-25. doi: 10.13275/j.cnki.lykxyj.2018.02.003
    [13] 王瑶孔祥波张苏芳刘福张真严善春 . 云南松毛虫线粒体基因组全序列测定和分析. 林业科学研究, 2019, 32(5): 11-20. doi: 10.13275/j.cnki.lykxyj.2019.05.002
    [14] 饶国栋睢金凯张建国 . 毛白杨α-微管蛋白基因家族的克隆与序列分析. 林业科学研究, 2015, 28(1): 44-49.
    [15] 宁坤杨洋马述山李慧玉 . 11条白桦BpSPL家族基因的生物信息学和表达分析. 林业科学研究, 2016, 29(5): 646-653.
    [16] 薛丽刘晓霞罗莹吕中睿张建国饶国栋 . 毛果杨MAP65基因家族的扩张与表达分析. 林业科学研究, 2021, 34(2): 92-101. doi: 10.13275/j.cnki.lykxyj.2021.02.010
    [17] 张恺恺杨立莹丰美静张林凤陈段芬邱德有杨艳芳 . 曼地亚红豆杉NAC基因家族鉴定及表达分析. 林业科学研究, 2022, 35(2): 97-103. doi: 10.13275/j.cnki.lykxyj.2022.02.011
    [18] 袁婷婷朱成磊杨克彬宋新章高志民 . 毛竹硝态氮转运蛋白家族PeNPFs基因鉴定及其表达特性分析. 林业科学研究, 2021, 34(3): 1-12. doi: 10.13275/j.cnki.lykxyj.2021.03.001
    [19] 胡景涛李彦杰段艳艳阮宇顾欣肖国生 . 桑树HD-Zip I亚家族基因的鉴定及表达分析. 林业科学研究, 2022, 35(4): 130-142. doi: 10.13275/j.cnki.lykxyj.2022.004.014
    [20] 詹妮谢耀坚陈鸿鹏刘果 . 巨桉EXP基因家族的生物信息学分析. 林业科学研究, 2018, 31(6): 39-46. doi: 10.13275/j.cnki.lykxyj.2018.06.006
  • 加载中
图(5) / 表(4)
计量
  • 文章访问数:  5676
  • HTML全文浏览量:  2685
  • PDF下载量:  294
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-08-22
  • 录用日期:  2018-04-20
  • 刊出日期:  2018-12-01

两种兰科植物CaMCML基因家族全基因组分析

    通讯作者: 刘蕾, liulei0828@126.com
    作者简介: 杨前宇(1985-), 男, 博士研究生, 研究方向为园林植物与观赏园艺.E-mail:yangdayu0055@163.com
  • 1. 中国林业科学研究院林业研究所, 林木遗传育种国家重点实验室, 国家林业局林木培育重点实验室, 北京 100091
  • 2. 沈阳农业大学林学院, 辽宁 沈阳 110866
基金项目:  北京市植物资源研究开发重点实验室开放课题 PRRD-2016-ZD1国家高技术研究发展计划(863计划)子课题 2013AA102607-3

摘要:  目的 探讨CaMCML基因家族在兰科植物不同组织和与菌根真菌建立共生关系过程中的潜在功能。 方法 依据已发表的拟南芥CaMCML基因家族蛋白序列,利用生物信息学方法对小兰屿蝴蝶兰和铁皮石斛全基因组进行搜索,通过多种生物软件或在线工具对CaMCML基因家族蛋白序列进行筛选及确定、蛋白结构分析、系统进化分析及结构域预测;利用转录组数据绘制heatmap图,对小兰屿蝴蝶兰不同组织(花、叶、茎、根)及石斛种子与美胞胶膜菌共生条件下CaMCML基因家族的表达情况进行分析。 结果 在小兰屿蝴蝶兰和铁皮石斛全基因组中均预测到4个CaM蛋白和54个CML蛋白;小兰屿蝴蝶兰CaMCML家族基因中39个基因没有内含子,19个基因具有内含子;铁皮石斛CaMCML基因家族中有41个基因没有内含子,17个基因具有内含子;系统进化分析将小兰屿蝴蝶兰和铁皮石斛CaM和CML蛋白家族各分为10个亚家族;小兰屿蝴蝶兰中9个基因在叶中的表达相对于花、茎、根为上调表达,2个基因在叶中的表达相对于花、茎、根却为下调表达;3个基因在花中的表达相对于叶、茎、根为上调表达;3个基因在花和叶中的表达相对于茎和根为上调表达;铁皮石斛中4个基因在与美胞胶膜菌共生萌发文库中的表达相对于非共生文库为上调表达;4个基因在共生萌发文库中的表达相对于非共生文库为下调表达。 结论 CaMCML家族基因在兰科植物不同组织具有重要作用,且参与铁皮石斛种子与美孢胶膜菌共生萌发的生物学调控过程。

English Abstract

  • 兰科(Orchidaceae)植物为多年生草本种子植物[1-2],约6个亚科800属35 000种[3-4]。兰科植物不仅有很高的观赏价值,还具有很高的药用和经济价值[5-7]。小兰屿蝴蝶兰(Phalaenopsis equestris (Schauer) Rchb.)[8]和铁皮石斛(Dendrobium officinale Kimura et Migo)[9-10]全基因组测序工作的完成,为科学工作者深入研究兰科植物在形态进化、生理适应和共生关系等方面的分子机制提供了海量的信息资源和前所未有的机遇。

    Ca2+是植物细胞在发育和应答外界信号中功能最广泛的第二信使。植物经长期进化,在细胞内形成一种独特的蛋白,这种蛋白拥有高度保守性结构(EF-hand motif),可以结合Ca2+,行使生物学功能,即:钙调蛋白(CaMs)和类钙调蛋白(CMLs)[11]。在植物细胞中,CaM和CML两类蛋白与Ca2+结合后形成的复合物参与多种生物学反应[12-13],是目前被研究最多的Ca2+受体蛋白[14],且在结构上非常保守,具有调控植物生长发育基因和生物胁迫条件下应答基因的表达作用[15-16]。目前,对兰科植物中CaMCML基因家族的研究尚未见报道,本文利用生物信息学的方法预测了小兰屿蝴蝶兰和铁皮石斛中的CaMCML基因家族成员,同时分析了这些成员的蛋白保守结构域,系统进化关系以及基因表达模式,为兰科植物CaMCML基因家族功能验证提供数据基础。

    • 通过OrchidBase数据库[8](http://orchidbase.itps.ncku.edu.tw/)下载小兰屿蝴蝶兰和铁皮石斛的全基因蛋白序列文件;通过TAIR基因组序列数据库[17](http://www.arabidopsis.org/)下载拟南芥(Arabidopsis thaliana)CaM和CML蛋白序列。将拟南芥的CaM和CML蛋白序列作为查询序列,利用BLASTP程序对小兰屿蝴蝶兰和铁皮石斛蛋白序列文件进行本地化比对检索,手动去除错配序列后,依据与拟南芥序列的同源性对其进行命名[18];利用DNAMAN软件[19]计算CaM及CML蛋白序列的分子量及等电点等;通过在线分析服务GSDS[20]对小兰屿蝴蝶兰和铁皮石斛CaMCML基因家族内含子及外显子结构进行分析。

    • 将小兰屿蝴蝶兰和铁皮石斛的候选蛋白序列与拟南芥CaM蛋白序列通过软件ClustalX 2.0.10 [17]进行序列比对。利用在线数据库SMART[18]和Pfam[19]对候选蛋白序列的保守结构域进行分析和鉴定。通过在线软件Myristoylator[20](http://web.expasy.org/myristoylator/)确定是否存在潜在豆蔻酰化位点。

    • 利用软件MEGA6.05[19]分别构建小兰屿蝴蝶兰和铁皮石斛CaM及CML蛋白家族系统进化树,建树方法为最大似然法(Maximum Likelihood,ML)[19],Bootstrap运行次数设置为1 000,其他参数均为默认。

    • 从兰科植物基因组数据库(http://www.nature.com/ng/journal/v47/n1/full/ng.3149.html)下载小兰屿蝴蝶兰不同组织(花、叶、茎、根)[8]以及利用本实验室获得的石斛种子与美胞胶膜菌互作转录组数据[21]。将RPKM(Reads Per Kilo bases per Million reads)值通过HemI(Heatmap Illustrator)1.0.3.3软件进行作图分析[22],分别统计CaMCML家族基因表达情况。

    • 从小兰屿蝴蝶兰和铁皮石斛基因组中均筛选获得4个CaM(OPeCaM1~4,ODcCaM1~4)和54个CML(OPeCML1~54,ODcCML1~54),其中,CaM蛋白序列全长为148~149 aa,分子量为11.96~17.96 kDa,等电点为3.67~3.95,都具有4个手性分子结构域(EF-hand);CML蛋白序列全长为72~337 aa(除 ODcCML53 为691aa),分子量为6.95~76.94 kDa,等电点为3.52~8.85,分别具有1~5个手性分子结构域;小兰屿蝴蝶兰CaMCML家族基因中39个基因没有内含子,19个基因具有内含子;铁皮石斛CaMCML基因家族中有41个基因没有内含子,17个基因具有内含子(表 1表 2)。小兰屿蝴蝶兰OPeCaM1~4的相似度为91.95%~95.97%,相差的氨基酸个数为6~12个;铁皮石斛ODcCaM1~4的相似度为92.00%~98.66%,相差的氨基酸个数为2~12个(表 3表 4)。

      表 1  小兰屿蝴蝶兰OPeCaM及OPeCML蛋白序列特征

      Table 1.  The characteristics of OPeCaM/OPeCML protein sequence in P.equestris

      名称
      Name
      基因ID
      Gene ID
      亚家族
      subfamily
      蛋白序列长度
      Length of protein/aa
      EF hands 分子量
      Mol Wt/kDa
      等电点
      pI
      内含子
      Intron
      甲硫氨酸的含量
      percentage methionine/%
      OPeCaM1 PEQU_18438 1 149 4 16.96 3.89 2 6.02
      OPeCaM2 PEQU_10717 1 149 4 16.96 3.89 2 6.01
      OPeCaM3 PEQU_12537 1 148 4 16.96 3.88 2 6.00
      OPeCaM4 PEQU_26335 1 149 4 17.96 3.95 1 6.05
      OPeCML01 PEQU_04925 1 148 4 12.96 3.63 4 4.82
      OPeCML02 PEQU_07491 10 155 4 16.96 3.95 1 5.01
      OPeCML03 PEQU_06469 10 166 4 16.96 4.02 0 5.90
      OPeCML04 PEQU_16519 1 135 3 20.96 3.80 0 4.11
      OPeCML05 PEQU_34003 10 199 4 30.98 4.21 6 4.67
      OPeCML06 PEQU_29891 1 120 3 14.96 3.52 2 2.83
      OPeCML07 PEQU_25947 7 190 4 24.95 4.16 1 6.36
      OPeCML08 PEQU_22407 6 155 4 19.97 3.98 0 3.42
      OPeCML09 PEQU_40659 5 154 4 18.95 3.95 0 3.71
      OPeCML10 PEQU_35611 1 72 2 8.96 3.73 1 2.44
      OPeCML11 PEQU_15751 10 141 3 17.96 3.83 4 2.12
      OPeCML12 PEQU_06469 10 166 4 16.97 4.19 0 6.23
      OPeCML13 PEQU_41809 6 219 4 24.95 4.38 0 5.22
      OPeCML14 PEQU_30031 8 137 4 20.96 3.88 0 5.01
      OPeCML15 PEQU_40424 2 193 4 23.97 4.07 0 4.93
      OPeCML16 PEQU_18134 6 199 4 25.97 4.18 0 5.22
      OPeCML17 PEQU_24120 8 134 4 18.96 3.80 0 5.33
      OPeCML18 PEQU_32857 8 140 4 17.95 3.81 0 4.31
      OPeCML19 PEQU_32861 8 141 4 16.97 3.82 0 5.46
      OPeCML20 PEQU_11970 8 140 4 17.96 3.87 0 6.82
      OPeCML21 PEQU_32860 8 141 4 15.96 3.87 0 5.11
      OPeCML22 PEQU_11969 8 237 5 25.95 4.28 1 6.30
      OPeCML23 PEQU_25577 8 143 3 18.95 4.38 0 4.41
      OPeCML24 PEQU_28290 8 143 4 15.96 3.77 0 5.02
      OPeCML25 PEQU_05240 3 143 3 16.97 3.88 0 4.92
      OPeCML26 PEQU_19202 5 160 3 18.97 3.96 0 3.83
      OPeCML27 PEQU_17666 5 197 3 22.95 4.19 0 3.77
      OPeCML28 PEQU_35797 6 110 3 11.96 3.81 0 3.61
      OPeCML29 PEQU_11304 8 152 3 19.96 3.85 0 4.16
      OPeCML30 PEQU_20039 6 114 3 15.96 3.77 0 3.03
      OPeCML31 PEQU_28291 8 110 3 11.96 3.78 0 4.23
      OPeCML32 PEQU_09901 8 111 3 16.96 3.79 0 3.19
      OPeCML33 PEQU_29894 2 101 3 6.95 3.78 0 2.68
      OPeCML34 PEQU_32859 8 115 3 14.96 3.78 0 3.46
      OPeCML35 PEQU_09640 9 183 4 38.96 3.87 0 3.10
      OPeCML36 PEQU_15088 4 166 2 23.95 3.85 0 2.89
      OPeCML37 PEQU_28312 4 176 2 20.97 3.83 0 2.54
      OPeCML38 PEQU_00748 7 168 3 7.96 3.85 0 3.88
      OPeCML39 PEQU_39569 10 108 2 14.95 3.81 0 3.16
      OPeCML40 PEQU_40143 10 127 2 18.96 3.77 0 1.93
      OPeCML41 PEQU_09077 4 103 2 13.96 3.88 0 2.42
      OPeCML42 PEQU_09518 4 158 2 21.96 3.89 0 2.02
      OPeCML43 PEQU_24484 10 225 4 34.95 4.28 3 5.33
      OPeCML44 PEQU_00516 7 106 2 7.97 3.77 0 2.11
      OPeCML45 PEQU_11010 1 304 2 30.96 4.75 4 1.66
      OPeCML46 PEQU_04326 1 253 2 23.96 4.56 0 3.35
      OPeCML47 PEQU_41809 6 219 4 24.96 4.23 5 4.26
      OPeCML48 PEQU_17666 5 197 3 22.97 4.06 0 3.28
      OPeCML49 PEQU_01306 10 153 3 18.95 3.99 1 3.66
      OPeCML50 PEQU_12715 10 147 3 17.97 4.20 2 2.82
      OPeCML51 PEQU_17236 10 211 4 17.96 4.12 1 4.21
      OPeCML52 PEQU_00999 1 113 4 13.96 4.32 0 3.16
      OPeCML53 PEQU_12714 10 153 2 19.95 4.41 1 3.19
      OPeCML54 PEQU_01183 1 121 3 13.96 3.52 0 2.47

      表 2  铁皮石斛ODcCaM及ODcCML蛋白序列特征

      Table 2.  The characteristics of ODcCaM/ODcCMLprotein sequence in D. officinale

      名称
      Name
      基因ID
      Gene ID
      亚家族
      subfamily
      蛋白序列长度
      Length of protein/aa
      EF hands 分子量
      Mol Wt/kDa
      等电点
      pI
      内含子
      Intron
      甲硫氨酸的含量
      percentage methionine /%
      ODcCaM1 DG_10048944 1 149 4 16.96 3.89 1 6.01
      ODcCaM2 DG_10042177 1 149 4 16.96 3.88 1 6.03
      ODcCaM3 DG_10042221 1 149 4 16.96 3.88 3 6.00
      ODcCaM4 DG_10069718 1 148 4 11.96 3.67 3 6.02
      ODcCML01 DG_10099596 8 147 4 16.96 3.63 3 5.11
      ODcCML02 DG_10026132 6 154 4 22.96 4.20 0 3.96
      ODcCML03 DG_10032884 6 152 4 18.96 4.12 0 4.88
      ODcCML04 DG_10102530 6 154 4 18.97 4.32 0 4.46
      ODcCML05 DG_10001096 6 165 4 22.96 4.41 0 5.52
      ODcCML06 DG_10086740 6 142 4 12.95 3.52 0 4.26
      ODcCML07 DG_10112859 6 154 4 21.96 4.16 0 4.14
      ODcCML08 DG_10069652 6 140 4 20.15 3.78 0 3.82
      ODcCML09 DG_10048052 6 143 4 16.96 3.75 0 6.26
      ODcCML10 DG_10013983 6 145 4 12.96 3.73 0 3.93
      ODcCML11 DG_10087227 6 145 4 11.96 3.73 0 3.74
      ODcCML12 DG_10120493 10 113 3 13.95 3.79 0 2.43
      ODcCML13 DG_10127394 2 112 3 13.96 3.78 0 2.16
      ODcCML14 DG_10045310 2 177 3 17.97 4.88 1 6.29
      ODcCML15 DG_10099614 9 124 3 7.96 3.77 2 5.26
      ODcCML16 DG_10006683 10 86 2 9.98 3.88 1 1.19
      ODcCML17 DG_10099594 9 130 3 14.96 3.89 2 3.95
      ODcCML18 DG_10012949 8 107 2 11.97 3.80 1 2.27
      ODcCML19 DG_10018482 6 110 3 13.96 3.82 0 3.36
      ODcCML20 DG_10068885 6 206 4 25.96 4.67 0 4.42
      ODcCML21 DG_10080246 6 206 4 25.96 4.67 0 5.16
      ODcCML22 DG_10107219 7 195 4 27.95 4.28 0 4.82
      ODcCML23 DG_10115967 6 219 4 24.96 4.38 0 5.56
      ODcCML24 DG_10045416 7 165 4 25.96 3.77 0 6.01
      ODcCML25 DG_10131684 7 151 2 21.95 3.88 0 2.44
      ODcCML26 DG_10013708 8 87 2 6.98 3.72 1 2.22
      ODcCML27 DG_10086551 8 147 3 21.96 3.89 0 4.01
      ODcCML28 DG_10057311 4 141 3 10.97 3.85 0 2.82
      ODcCML29 DG_10047181 7 104 2 10.96 3.78 0 1.71
      ODcCML30 DG_10112504 3 231 3 39.95 4.67 0 3.29
      ODcCML31 DG_10114938 8 303 3 37.97 5.88 6 4.15
      ODcCML32 DG_10061210 5 107 3 12.96 3.79 3 3.23
      ODcCML33 DG_10067799 7 100 3 7.96 3.78 0 4.53
      ODcCML34 DG_10089833 7 195 4 28.97 4.58 0 5.14
      ODcCML35 DG_10069653 6 116 3 18.95 3.87 0 2.96
      ODcCML36 DG_10112846 6 133 3 21.96 3.85 0 3.54
      ODcCML37 DG_10112841 6 128 3 20.96 3.83 0 3.26
      ODcCML38 DG_10122957 6 133 3 21.95 3.85 0 2.80
      ODcCML39 DG_10012468 7 190 4 24.97 4.58 0 2.51
      ODcCML40 DG_10113507 6 185 4 25.97 4.17 0 3.82
      ODcCML41 DG_10112861 6 337 4 48.95 5.88 2 3.99
      ODcCML42 DG_10140024 6 148 2 16.96 3.89 1 1.58
      ODcCML43 DG_10008482 6 184 3 30.96 4.28 0 3.44
      ODcCML44 DG_10088170 6 175 3 31.97 4.02 0 2.15
      ODcCML45 DG_10078503 7 189 4 17.96 4.27 0 5.52
      ODcCML46 DG_10098264 7 189 4 17.95 4.27 0 4.32
      ODcCML47 DG_10078224 7 195 3 24.96 4.99 0 3.66
      ODcCML48 DG_10002690 7 107 2 11.95 3.72 0 2.32
      ]ODcCML49 DG_10049524 7 176 2 23.93 4.25 0 1.20
      ODcCML50 DG_10045416 7 165 2 25.97 4.17 0 2.21
      ODcCML51 DG_10024582 7 181 2 24.96 4.28 0 3.03
      ODcCML52 DG_10116046 7 327 2 29.97 5.88 4 2.52
      ODcCML53 DG_10064998 7 691 2 76.94 8.85 10 2.23
      ODcCML54 DG_10098265 5 191 4 25.95 4.58 0 4.19

      表 3  小兰屿蝴蝶兰CaM蛋白序列比对(OPeCaMs)

      Table 3.  Alignment of protein sequence of OPeCaMs

      OPeCaM1 OPeCaM2 OPeCaM3 OPeCaM4
      OPeCaM1 6 9 12
      OPeCaM2 95.97 10 11
      OPeCaM3 93.96 93.29 12
      OPeCaM4 91.95 92.62 91.95
        注:数字分别表示相似度(100%)和氨基酸差异数。
        Note: The numbersaresimilarity(100%) and quantity difference amino acids.

      表 4  铁皮石斛CaM蛋白序列比对(ODcCaMs)

      Table 4.  Alignment of protein sequence of ODcCaMs

      ODcCaM1 ODcCaM2 ODcCaM3 ODcCaM4
      ODcCaM1 3 3 12
      ODcCaM2 97.99 2 11
      ODcCaM3 97.99 98.66 11
      ODcCaM4 92.00 92.67 92.67
        注:数字分别表示相似度(100%)和氨基酸差异数。
        Note: The numbersaresimilarity(100%) and quantity difference amino acids.
    • 利用ClustalX 2.0.10软件将OPeCaM和ODcCaM蛋白序列与拟南芥AtCaM1进行多重序列比对(图 1)。兰科植物CaM蛋白具有能与Ca2+结合的4个手性分子结构域(EF-hand),每个EF-hand都是由30个氨基酸组成的螺旋-环-螺旋结构(Helix-Loop-Helix)。第1至9位氨基酸组成E螺旋(E-helix),第22至30位氨基酸组成的F螺旋(F-helix),中间的12个氨基酸为结合Ca2+的环(Loop),均为9个氨基酸长度的E螺旋与F螺旋结合于钙离子环的两侧。E螺旋由1个谷氨酸(E)起始,疏水氨基酸依“h**h**h”的规则在序列分布,其中,“h”代表疏水氨基酸,“*”为任意一种氨基酸。这种特殊的分布模式同样适用于EF-hand1和EF-hand3手性分子的F螺旋,而EF-hand2和EF-hand4手性分子的F螺旋的氨基酸分布模式与该模式有细微的差别。中间Loop环上的1、3位点是非常保守的天冬氨酸(D);第5位点通常为天冬氨酸(D)或者天冬酰胺(N),但占据第7、9位点的氨基酸有多种可能;谷氨酸(E)在0位点,功能是连接2个配位点;甘氨酸(G)在6位点,功能是保持环的结构;异亮氨酸(I)在8位点,功能是形成氢键与其他手性分子环结合成对。

      图  1  兰科植物与拟南芥CaM蛋白序列比对

      Figure 1.  The protein sequence alignment of OPeCaMs, ODcCaMs and AtCaM1

      在兰科植物CaM蛋白序列27位上存在1个半胱氨酸(C),常见于植物钙调蛋白中,在动物钙调蛋白中却没有被发现[23];在116位上存在1个位于第3个和第4个手性分子之间的赖氨酸(K),这一位点被认为是三甲基作用的潜在位点,存在于除酵母CaM以外的所有真核生物的CaM中[24]

      另外,OPeCaM和ODcCaM的甲硫氨酸的含量均为6%左右。OPeCML5、OPeCML15、OPeCML33、OPeCML35、ODcCML6、ODcCML7、ODcCML8、ODcCML18、ODcCML28、ODcCML42都包含潜在的豆蔻酰化序列(表 1表 2)。完整的CaM结构是CaM结合Ca2+和靶蛋白结合进而行使功能的前提。

    • 利用最大似然法(ML)构建小兰屿蝴蝶兰和铁皮石斛的CaM及CML氨基酸序列系统进化树(图 23)。依据拟南芥CaM及CML亚家族分类方法(9个亚家族)[25],将小兰屿蝴蝶兰和铁皮石斛CaM及CML蛋白各自分成10个亚家族,分别命名为S1~S10,其中,OPeCaM1~4与OPeCML01、OPeCML04、OPeCML06、OPeCML10、OPeCML45、OPeCML46、OPeCML52、OPeCML54聚类在同一个亚家族里,而ODcCaM1~4聚类的亚家族里没有ODcCMLn;OPeCML25、OPeCML35、ODcCML28、ODcCML30均单独聚类到一个亚家族。将Bootstrap≥75%分支加粗,在小兰屿蝴蝶兰系统进化树中S2、S3、S5的Bootstrap值均在75%~100%,并且各分支中Bootstrap值≥75%较多;铁皮石斛系统进化树中S5的Bootstrap值为97%,其他各亚群中Bootstrap值≥75%也较多。说明构建的系统进化树可信值较高。

      图  2  小兰屿蝴蝶兰CaM及CML家族系统进化树

      Figure 2.  The phylogenetic tree of OPeCaM and OPeCML

      图  3  铁皮石斛CaM及CML家族系统进化树

      Figure 3.  The phylogenetic tree of ODcCaM and ODcCML

    • 为确定CaMCML基因在兰科植物不同组织中的表达情况,本文利用小兰屿蝴蝶兰花、叶、茎、根4个组织的转录组测序数据进行基因表达分析,结果(图 4)表明:OPeCML13OPeCML47OPeCML51OPeCML23OPeCML30OPeCML26OPeCML08OPeCML07OPeCML15在叶中的表达相对于花、茎、根为上调表达,表达量最高的是OPeCML08(RPKM为1 178);而OPeCML05OPeCML54在叶中的表达相对于花、茎、根却为下调表达;OPeCaM1OPeCML16OPeCML38在花中的表达相对于叶、茎、根为上调表达;OPeCaM4OPeCML07OPeCML15在花和叶中的表达相对于茎和根为上调表达。

      图  4  小兰屿蝴蝶兰CaMCML基因在4个组织部位表达水平热图

      Figure 4.  Heatmap showing the expression levels of OPeCaM and OPeCML genes in four tissues

      为确定CaMCML基因在铁皮石斛种子与菌根菌共生萌发过程中的表达水平,本文利用铁皮石斛体外非共生萌发种子(ODc)和共生萌发种子(ODcTcg)转录组数据进行ODcCaMODcCML基因的表达分析, 结果(图 5)表明:ODcCaM1ODcCML05ODcCML14ODcCML49在共生萌发文库中的表达相对于非共生文库为上调表达,表达量最高的是ODcCaM1(RPKM为8 007);而ODcCaM3ODcCML13ODcCML27ODcCML52在共生萌发文库中的表达相对于非共生文库为下调表达。

      图  5  CaMCML基因在共生和非共生萌发的铁皮石斛种子中的表达水平热图

      Figure 5.  Heatmap showing the expression levels of ODcCaM and ODcCML genes in symbiotic and asymbiotic germinated seeds

    • 钙是植物生长发育所必需的营养元素,也是一种转导多种生理过程的重要信号分子。钙调蛋白(CaM)以及类钙调蛋白(CML)是细胞内两种重要的钙结合靶蛋白,开展植物中钙调蛋白和类钙调蛋白的研究,对了解上述基因在植物生长发育、细胞代谢、环境胁迫以及生物胁迫等方面的作用具有重要意义。目前,CaM基因和CML基因已在多种植物中开展研究,如拟南芥(Arabidopsis thaliana (L.) Heynh.)基因组中发现6个AtCaM基因和50个AtCML基因[25];水稻(Oryza sativa L.)基因组中发现5个OsCaM基因和32个OsCML基因[26];毛果杨(Populus trichocarpa Torr. & Gray)基因组中发现9个PtCaM基因和50个PtCML基因[27];胡杨(Populus euphratica)基因组中发现8个PeCaM基因和51个PeCML基因[27]。Zhao等[28]于2013年完成了对茄属(Solanum L.)植物CaM基因的鉴别和功能分析,分别于番茄(Lycopersicon esculentum Miller),烟草(Nicotiana tabacum L.)和马铃薯(Solanum tuberosum L.)的基因组中鉴定到6、7和4个CaM基因。本研究利用生物信息学方法预测分析了小兰屿蝴蝶兰和铁皮石斛的CaMCML基因家族,发现2种兰科植物基因组中均含有4个CaM基因和54个CML基因,较上述植物中CaMCML基因数量少。对比分析NCBI所发表的基因组数据发现,拟南芥、水稻、毛果杨和胡杨分别有33 583、30 534、45 555和34 279个蛋白编码基因,而小兰屿蝴蝶兰和铁皮石斛基因组分别有29 431个和35 567个蛋白编码基因,推测该基因家族的数量可能不与基因组蛋白编码基因的数量成正比。有研究表明,在拟南芥、水稻和杨树中,片段复制和串联复制可能是导致CaMCML基因家族扩张的主要原因[25-27]

      2003年Mccormack等[25]对拟南芥的CaMCML基因家族进行了分析,共找到了6个CaM基因和50个CML编码基因,其中,6个CaM基因编码的蛋白均含有4个手性分子结构域,推测6个CaM基因并不足以维持植株的功能,分布于50个CML中的手性分子结构域与CaM中手性分子作用类似;同时与CaM基因家族相比,CML基因家族在手性分子的数量和部分功能上有改变。本研究发现,在2种兰科植物基因组中均含有4个CaM基因和54个CML基因,CML基因数量远大于CaM基因数量,且CML基因具有1~5个数量不等的手性分子结构域。推测CML基因家族是CaM基因在进化过程中演变产生的,这一演变结果体现出明显的基因进化作用,同时CML中手性分子在序列上的差异可能反映了不同的CML基因家族成员间对Ca2+信号响应的微小差别[29]

      甲硫氨酸通过侧链异常的灵活性和极化性影响钙调蛋白的结构和功能,当CaM结合钙离子时,CaM形成开放构象,使结构更加稳定[30]。本研究发现,小兰屿蝴蝶兰和铁皮石斛中的CaM与拟南芥、水稻、杨树等植物的CaM甲硫氨酸含量相同,均为6%左右[25-27];豆蔻酰化序列的存在影响到对信号转导的调节,在拟南芥中就发现有4个AtCML蛋白序列存在潜在的豆蔻酰化位点[25],而水稻中只有1个[26],本研究在小兰屿蝴蝶兰和铁皮石斛中分别找到了4个和6个。

      Mccormack等[25]通过对拟南芥CaM及CML蛋白家族进行系统进化分析,将CaM及CML蛋白分成了9个亚家族;Boonburapong等[26]在2007年对水稻的CaM及CML蛋白进行了分析,将其分为6个亚家族。本研究依据拟南芥CaM及CML亚家族分类方法,将小兰屿蝴蝶兰和铁皮石斛分别分成了10个亚家族, 并通过Bootstrap值≥75%分支较多,推测出构建的系统进化树可信度较高。研究发现,AtCML49/50通过影响脱落酸介导(ABA-mediated)途径调控拟南芥在盐胁迫下的耐受性[25],OPeCML11和ODcCML05/30与AtCML49/50聚类到同一分支,推测它们在功能上存在相似性。拟南芥AtCML37通过茉莉酸信号途径结合钙离子从而调节其抗虫害胁迫的耐受性[25],同时OPeCML40/45和ODcCML52/53与AtCML37聚类到同一分支,由此推测, 他们可能调控兰花在抗虫害胁迫方面的耐受性。庞俊峰等[31]于2010年通过生物信息学分析了木榄(Bruguiera gymnorrhiza (L.) Poir.)的CaMCML基因家族,并与拟南芥、花生、大麦、榉木及甘薯等植物进行了比较,发现聚类到同一个亚家族的蛋白,不仅在结构上基本相同,所含有的外显子和内含子数也基本一致。通过本文CaMCML基因家族系统进化树和结构图分析,也可以得到相同的推测。

      大量研究表明,CaMCML家族基因参与植物不同组织在生长发育过程中的响应机制[32-34]。芒果(Mangifera indica L.)在花瓣分化、雄蕊分化、雌蕊分化时期时花中Ca2+的含量相比其它组织中的含量增加,CaM/CML含量也随之升高[33];梅花(Armeniaca mume Sieb.)中花芽的CaM/CML含量相比于其他组织也较高[34]。杨树在叶发育过程中,高表达量的PtCaMPtCML家族基因均参与钙调蛋白的积累过程[27]。在拟南芥矮小突变体和正常植株的研究中发现,AtCaMAtCML的含量不同,说明CaM/CML参与调控植物茎的生长发育[35-36]。Yang等[37]在种子生根过程的研究中发现,CaM/CML表达量随着根的发育逐渐降低。本研究对CaMCML基因在小兰屿蝴蝶兰花、叶、茎、根中表达水平分析发现,9个基因在叶中的表达相对于花、茎、根为上调表达,其中,OPeCML08表达量最高,RPKM为1 178;而OPeCML05OPeCML54在叶中的表达相对于花、茎、根却为下调表达。推测它们可能在小兰屿蝴蝶兰的叶中发挥着某些重要作用。OPeCaM1OPeCML16OPeCML38在花中的表达相对于叶、茎、根为上调表达;OPeCaM4OPeCML07OPeCML15在花和叶中的表达相对于茎和根为上调表达。推测这些基因可能在兰科植物不同组织过程中起到关键性作用。

      前期研究结果发现,在植物受到生物及非生物胁迫时,CaMCML家族基因具有重要作用[24]。Boonburapong等[26]针对水稻叶、根、花、未成熟的种子和愈伤组织PCR扩增反转录得到的cDNA,结果发现生物胁迫下的水稻叶子中OsCaM1-1基因表达量显著增多。这与“细胞中Ca2+的变化是高等植物对病虫害侵染的早期反应之一[38]”的结论吻合。大豆SCaM4和SCaM5的表达受到病原真菌的强烈诱导,这一诱导事件也可以通过外施Ca2+来实现,然而用Ca2+螯合剂处理后抑制了大豆钙调蛋白基因在病原菌诱导下的表达[39]。Yamakawa等[40]通过沉默NtCaM2NtCaM6基因,发现烟草(Nicotiana tabacum L.)对脆裂病毒、瓜果腐霉菌和真菌病原体这三种常见病害的抵抗力显著降低,推测NtCaM基因家族在抗生物胁迫方面发挥着显著的作用。烟草NtCaM4NtCaM5基因在正常环境下表达水平很低,但在病原体侵入的条件下表达量明显提高,表现出高诱导性[40]。马尾松(Pinus massoniana Lamb.)接种松材线虫后,马尾松苗为响应松材线虫侵染,根、茎、叶器官中的PmCaM基因均不同程度地上调表达[41]。Leba等[42]对拟南芥CML基因在多种胁迫条件下研究发现,AtCML9受到致病菌、鞭毛蛋白和水杨酸的诱导后上调表达。Harrison等[43]研究植物在受到丛枝菌根真菌侵染后,胞内Ca2+浓度迅速升高,推测Ca2+-CaM/CML信号通路被激活,CaMCML基因活跃表达。兰科植物在自然条件下必须从特定的真菌中获得营养才能成功萌发,并在随后的生活史中与真菌建立稳定的共生关系[5-6]。这些特定的真菌被定义为兰科菌根真菌,均属于丛枝菌根真菌[7]。因此,对兰科植物与兰科菌根真菌的共生关系建立过程中CaMCML基因家族的的表达情况进行分析,对进一步阐明共生机制具有重要意义。本研究发现,在非共生和共生的铁皮石斛种子中ODcCaM1ODcCML05ODcCML14ODcCML49在共生萌发文库中的表达相对于非共生文库为上调表达;而 ODcCaM3ODcCML13ODcCML27和ODcCML52在共生萌发文库中的表达相对于非共生文库为下调表达。推测这些基因可能参与铁皮石斛种子与美孢胶膜菌共生萌发中的生物学调控。

    • 本研究依据已发表的拟南芥CaMCML基因家族蛋白序列,利用生物信息学方法对小兰屿蝴蝶兰和铁皮石斛CaMCML基因家族进行分析。在小兰屿蝴蝶兰和铁皮石斛全基因组中均预测到4个CaM蛋白和54个CML蛋白;小兰屿蝴蝶兰CaMCML家族基因中39个基因没有内含子,19个基因有内含子;铁皮石斛CaMCML基因家族中有41个基因没有内含子,17个基因有内含子;系统进化分析将其各分为10个亚家族;小兰屿蝴蝶兰中9个基因在叶中的表达相对于花、茎、根为上调表达,2个基因在叶中的表达相对于花、茎、根却为下调表达,3个基因在花中的表达相对于叶、茎、根为上调表达,3个基因在花和叶中的表达相对于茎和根为上调表达。说明CaMCML家族基因在兰科植物不同组织中具有重要作用;铁皮石斛中4个基因在共生萌发文库中的表达相对于非共生文库为上调表达,4个基因在共生萌发文库中的表达相对于非共生文库为下调表达。说明CaMCML基因家族可能参与铁皮石斛种子与美孢胶膜菌共生萌发中的生物学调控。本研究为进一步研究兰科植物中CaMCML基因家族的潜在生物学功能提供数据基础。

参考文献 (43)

目录

    /

    返回文章
    返回