-
兰科(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]。目前,对兰科植物中CaM及CML基因家族的研究尚未见报道,本文利用生物信息学的方法预测了小兰屿蝴蝶兰和铁皮石斛中的CaM及CML基因家族成员,同时分析了这些成员的蛋白保守结构域,系统进化关系以及基因表达模式,为兰科植物CaM及CML基因家族功能验证提供数据基础。
-
从小兰屿蝴蝶兰和铁皮石斛基因组中均筛选获得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个手性分子结构域;小兰屿蝴蝶兰CaM及CML家族基因中39个基因没有内含子,19个基因具有内含子;铁皮石斛CaM及CML基因家族中有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/aaEF 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/aaEF 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位点,功能是形成氢键与其他手性分子环结合成对。
在兰科植物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氨基酸序列系统进化树(图 2、3)。依据拟南芥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%也较多。说明构建的系统进化树可信值较高。
-
为确定CaM及CML基因在兰科植物不同组织中的表达情况,本文利用小兰屿蝴蝶兰花、叶、茎、根4个组织的转录组测序数据进行基因表达分析,结果(图 4)表明:OPeCML13、OPeCML47、OPeCML51、OPeCML23、OPeCML30、OPeCML26、OPeCML08、OPeCML07和OPeCML15在叶中的表达相对于花、茎、根为上调表达,表达量最高的是OPeCML08(RPKM为1 178);而OPeCML05和OPeCML54在叶中的表达相对于花、茎、根却为下调表达;OPeCaM1、OPeCML16和OPeCML38在花中的表达相对于叶、茎、根为上调表达;OPeCaM4、OPeCML07和OPeCML15在花和叶中的表达相对于茎和根为上调表达。
图 4 小兰屿蝴蝶兰CaM及CML基因在4个组织部位表达水平热图
Figure 4. Heatmap showing the expression levels of OPeCaM and OPeCML genes in four tissues
为确定CaM和CML基因在铁皮石斛种子与菌根菌共生萌发过程中的表达水平,本文利用铁皮石斛体外非共生萌发种子(ODc)和共生萌发种子(ODcTcg)转录组数据进行ODcCaM和ODcCML基因的表达分析, 结果(图 5)表明:ODcCaM1、ODcCML05、ODcCML14和ODcCML49在共生萌发文库中的表达相对于非共生文库为上调表达,表达量最高的是ODcCaM1(RPKM为8 007);而ODcCaM3、ODcCML13、ODcCML27和ODcCML52在共生萌发文库中的表达相对于非共生文库为下调表达。
两种兰科植物CaM及CML基因家族全基因组分析
Genome-wide Analysis of CaM/CML Gene Family in Two Orchidaceae Species
-
摘要:
目的 探讨CaM和CML基因家族在兰科植物不同组织和与菌根真菌建立共生关系过程中的潜在功能。 方法 依据已发表的拟南芥CaM和CML基因家族蛋白序列,利用生物信息学方法对小兰屿蝴蝶兰和铁皮石斛全基因组进行搜索,通过多种生物软件或在线工具对CaM和CML基因家族蛋白序列进行筛选及确定、蛋白结构分析、系统进化分析及结构域预测;利用转录组数据绘制heatmap图,对小兰屿蝴蝶兰不同组织(花、叶、茎、根)及石斛种子与美胞胶膜菌共生条件下CaM和CML基因家族的表达情况进行分析。 结果 在小兰屿蝴蝶兰和铁皮石斛全基因组中均预测到4个CaM蛋白和54个CML蛋白;小兰屿蝴蝶兰CaM及CML家族基因中39个基因没有内含子,19个基因具有内含子;铁皮石斛CaM及CML基因家族中有41个基因没有内含子,17个基因具有内含子;系统进化分析将小兰屿蝴蝶兰和铁皮石斛CaM和CML蛋白家族各分为10个亚家族;小兰屿蝴蝶兰中9个基因在叶中的表达相对于花、茎、根为上调表达,2个基因在叶中的表达相对于花、茎、根却为下调表达;3个基因在花中的表达相对于叶、茎、根为上调表达;3个基因在花和叶中的表达相对于茎和根为上调表达;铁皮石斛中4个基因在与美胞胶膜菌共生萌发文库中的表达相对于非共生文库为上调表达;4个基因在共生萌发文库中的表达相对于非共生文库为下调表达。 结论 CaM及CML家族基因在兰科植物不同组织具有重要作用,且参与铁皮石斛种子与美孢胶膜菌共生萌发的生物学调控过程。 Abstract:Objective To research the CaM/CML gene family's characteristics in different tissues and symbiosis with mycorrhizal fungi of Orchidaceae. Method Genome-wides of Phalaenopsis equestris and Dendrobium officinale were separately predicted and determined with bioinformatics methods based on the publishing protein sequences of CaM/CML gene family of Arabidopsis. The protein structures, phylogenetic trees, exon-intron structures of CaM/CML gene family in P. equestris and D. officinale were analyzed using different biological softwares. Gene expression analysis of different tissues (flower, leave, stem and root) in P. equestris and symbiotic with Tulasnella calospora and asymbiotic germination seeds from D. officinale were conducted by using RNA-seq data. Result Four CaMs and fifty-four CMLs were separately predicted from P. equestris and D. officinale genomes. In P. equestris, forty-one genes did not have intron, seventeen genes had intron in CaM/CML gene family; and thirty-nine genes did not have intron, nineteen genes had intron in D. officinale. Phylogenetic analysis results showed that the CaMs and CMLs could be separately divided into ten subfamilies. The results of CaMs and CMLs genes expression analysis in different tissues of P. equestris indicated that nine genes were up-regulated in the leaf relative to the flower, stem and root, two genes were down-regulated; three genes were up-regulated in the flower relative to the leaf, stem and root; three genes were up-regulated in the flower and leaf relative to the stem and root. The results of gene expression analysis in symbiotic with T. calospora and asymbiotic germinated seeds from D. officinale by using RNA-seq data indicated that four genes were up-regulated, and four genes down-regulated. Conclusion These results will lay the basis for the further studying of potential biological function and molecular mechanism of CaMs and CMLs genes in orchid. -
表 1 小兰屿蝴蝶兰OPeCaM及OPeCML蛋白序列特征
Table 1. The characteristics of OPeCaM/OPeCML protein sequence in P.equestris
名称
Name基因ID
Gene ID亚家族
subfamily蛋白序列长度
Length of protein/aaEF 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/aaEF 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. -
[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] 王雪霁.杨树CaM及CML基因家族生物信息学分析[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