富含半胱氨酸类受体蛋白激酶（cysteinerich receptor-like kinase, CRK）在植物类受体蛋白激酶（receptor-like kinase, RLK）家族中具有重要地位,在植物生长发育、激素信号传导、生物和非生物胁迫以及植物抗病性细胞过敏性死亡中发挥着至关重要的作用[1,2]。
国内外研究发现CRK家族基因可调控植物生长发育及抵抗逆境。在对拟南芥的CRK基因研究中：CRK45通过调控ABA的生物合成而调控着对非生物逆境胁迫,如干旱、高盐以及ABA的响应,对植物耐药性有正向调控作用,即与抗病性相关,参与了水杨酸介导的抗病信号途径[3];CRK5是参与调节基本生理过程的生化机制的重要组成部分,CRKs突变体的功能特性表现出气孔导度降低和衰老加速,与活性氧种类的积累有关,CRK5可以在紫外线照射下表现出细胞死亡和氧化损伤的增强[4],转基因拟南芥对紫外线强度的适应性可通过过表达CRK5来增强[4];CRK17受水杨酸诱导表达参与植物的抗病防御反应[5];用诱导型启动子过量表达CRK4、CRK5、CRK19、CRK20能引起过敏反应导致细胞程序化死亡[4,5,6];CRK13可通过依赖水杨酸抑制病原菌[7]。从番茄分离的CRK1基因在花粉中可以特异性表达[8];菜豆中CRK基因PvPK20-1可调控病原菌、共生菌与肿瘤因子[9];水稻中OsCRK1基因沉默后严重影响水稻对白叶枯病的抗性,可证明其与叶枯病抗性密切相关[10];在甜瓜中,大多数CRK基因在白粉病和ABA胁迫后下调表达,在盐和低温处理后上调表达,并对冷害胁迫表现较强烈[11]。
龙眼（Dimocarpus longan Lour.）是我国重要的亚热带果树,种植面积及产量均居世界首位[12]。但高温、干旱、盐碱、寒冷等非生物胁迫以及生物胁迫均会对龙眼的生长发育造成不良影响,且现今对龙眼生物和非生物胁迫耐受的分子调控机制了解较为匮乏。因此,通过对CRK基因家族的分析,在龙眼生长发育以及抗逆性研究中具有重大意义。本文通过对龙眼CRK家族成员进行分析鉴定以及表达调控研究（进化树分析、基因结构分析、启动子顺式作用元件分析、蛋白保守结构域分析以及在不同体胚发生过程中和生长发育阶段不同组织器官特异性表达的FPKM值分析）,了解了CRK家族基因及其在龙眼中的表达情况,可为CRK家族在龙眼生长发育及对其他植物调控的研究中提供理论基础。
1 材料与方法
1.1 材料
1.1.1 实验材料 实验材料为福建农林大学园艺植物生物工程研究所提供的‘红核子’品种龙眼体胚发生早期胚性愈伤组织（embryonic callus, EC）、非胚性愈伤组织（non-embryoniccallus, NEC）、不完全胚性紧实结构（incomplete embryonic compact structure, ICpEC）、球形胚（globular embryo, GE）[13]。
1.1.2 数据来源 实验室构建的龙眼基因组数据库（NCBI登录号：BioProject PRJNA305337）;龙眼转录组数据库（SRA050205）;拟南芥（Arabidopsis thaliana）CRK基因家族成员序列下载于https://phytozome.jgi.doe.gov/pz/portal. html。
1.2 方法
1.2.1 龙眼CRK基因家族成员的鉴定及其基本理化性质分析 对龙眼CRK基因家族成员进行鉴定,经龙眼基因组数据库同源对比,筛选获得98条具有完整Open Reading Frame（ORF）的龙眼CRK候选序列。根据基因注释,采用DNAMAN 6.0软件对CRK的gDNA、CDS以及启动子序列进行比对分析,结合NCBI Blast同源比对分析,初步确定龙眼CRK家族成员。其次,根据其基因座位置[14],将龙眼CRK家族成员命名为CRK1-CRK98（表1）。此外,从ExPASy（www.expasy.org）中获取CRK家族的基本理化性质,并对其分子量（Mw）、氨基酸个数（aa）、等电点（pI）、信号肽（signal peptide）及其糖基化位点（N-glycosyl sites）进行系统的分析。
表1 DlCRK家族基因编码的蛋白理化性质分析Tab. 1 Physicochemical properties analysis of longan DlCRK gene family members coding proteins |
| 基因名称Gene name | 基因ID Gene ID | 位置 Location | 氨基酸个数 Number of acmino aid | 分子量 Mw/ku | 等电点 pI | 信号肽 Signal peptide | 糖基化位点 N-glycosyl sites | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| DlCRK1 | Dlo_008262.1 | scaffold17:1418468:1422193 | 649 | 72.84 | 5.31 | YES | 11 | |||||||||||
| DlCRK2 | Dlo_011267.1 | scaffold21:1649322:1653952 | 915 | 101.30 | 7.45 | NO | 9 | |||||||||||
| DlCRK3 | Dlo_011268.1 | scaffold21:1656389:1659473 | 670 | 74.56 | 5.25 | YES | 8 | |||||||||||
| DlCRK4 | Dlo_011270.1 | scaffold21:1700403:1703482 | 664 | 74.14 | 5.76 | YES | 6 | |||||||||||
| DlCRK5 | Dlo_012981.1 | scaffold24:431411:434686 | 564 | 64.81 | 9.49 | NO | 5 | |||||||||||
| DlCRK6 | Dlo_015369.1 | scaffold29:482575:485078 | 641 | 70.49 | 8.22 | NO | 6 | |||||||||||
| DlCRK7 | Dlo_016006.1 | scaffold30:248448:253318 | 636 | 71.63 | 7.79 | YES | 8 | |||||||||||
| DlCRK8 | Dlo_016054.1 | scaffold30:1053225:1056696 | 724 | 81.94 | 5.51 | YES | 11 | |||||||||||
| DlCRK9 | Dlo_016055.1 | scaffold30:1072291:1079097 | 793 | 88.37 | 5.16 | NO | 10 | |||||||||||
| DlCRK10 | Dlo_016056.3 | scaffold30:1107155:1110588 | 745 | 84.11 | 5.47 | YES | 11 | |||||||||||
| DlCRK11 | Dlo_017641.1 | scaffold34:1498040:1501933 | 777 | 87.86 | 6.11 | YES | 9 | |||||||||||
| DlCRK12 | Dlo_018615.1 | scaffold37:839212:841936 | 651 | 72.51 | 8.28 | YES | 9 | |||||||||||
| DlCRK13 | Dlo_018616.2 | scaffold37:845764:848709 | 596 | 66.03 | 8.66 | NO | 11 | |||||||||||
| DlCRK14 | Dlo_019866.1 | scaffold40:207589:211282 | 864 | 97.16 | 8.17 | YES | 9 | |||||||||||
| DlCRK15 | Dlo_019870.1 | scaffold40:251674:262633 | 1063 | 118.54 | 8.56 | NO | 10 | |||||||||||
| DlCRK16 | Dlo_019871.1 | scaffold40:272347:280591 | 771 | 86.36 | 5.91 | YES | 10 | |||||||||||
| DlCRK17 | Dlo_019872.1 | scaffold40:338615:342561 | 755 | 85.00 | 6.36 | YES | 11 | |||||||||||
| 基因名称Gene name | 基因ID Gene ID | 位置 Location | 氨基酸个数 Number of amino aid | 分子量 Mw/ku | 等电点 pI | 信号肽 Signal peptide | 糖基化位点 N-glycosyl sites | |||||||||||
| DlCRK18 | Dlo_019878.1 | scaffold40:390413:392855 | 430 | 48.66 | 5.33 | NO | 2 | |||||||||||
| DlCRK19 | Dlo_019879.1 | scaffold40:397709:400689 | 772 | 86.85 | 5.82 | YES | 8 | |||||||||||
| DlCRK20 | Dlo_019880.1 | scaffold40:455598:458720 | 699 | 79.00 | 5.25 | NO | 11 | |||||||||||
| DlCRK21 | Dlo_019888.1 | scaffold40:581049:584203 | 769 | 86.61 | 5.09 | NO | 10 | |||||||||||
| DlCRK22 | Dlo_032127.1 | scaffold87:321002:324107 | 637 | 72.03 | 8.39 | NO | 11 | |||||||||||
| DlCRK23 | Dlo_032129.1 | scaffold87:340063:347128 | 693 | 76.83 | 6.99 | NO | 8 | |||||||||||
| DlCRK24 | Dlo_032130.1 | scaffold87:359813:362708 | 590 | 65.82 | 6.41 | YES | 7 | |||||||||||
| DlCRK25 | Dlo_032131.1 | scaffold87:386688:389578 | 666 | 74.73 | 5.94 | YES | 9 | |||||||||||
| DlCRK26 | Dlo_032137.1 | scaffold87:472245:475145 | 643 | 71.83 | 6.11 | YES | 9 | |||||||||||
| DlCRK27 | Dlo_032138.2 | scaffold87:484271:504722 | 893 | 100.00 | 5.45 | NO | 13 | |||||||||||
| DlCRK28 | Dlo_032139.1 | scaffold87:525085:527717 | 554 | 61.93 | 6.31 | NO | 7 | |||||||||||
| DlCRK29 | Dlo_032140.1 | scaffold87:558625:561567 | 667 | 74.51 | 8.01 | YES | 9 | |||||||||||
| DlCRK30 | Dlo_032141.1 | scaffold87:565433:568357 | 685 | 76.81 | 6.52 | NO | 7 | |||||||||||
| DlCRK31 | Dlo_032143.1 | scaffold87:583584:586963 | 755 | 84.15 | 7.46 | NO | 8 | |||||||||||
| DlCRK32 | Dlo_033594.1 | scaffold96:72103:75322 | 782 | 88.64 | 5.77 | YES | 12 | |||||||||||
| DlCRK33 | Dlo_033616.1 | scaffold96:432411:442605 | 772 | 87.12 | 6.56 | NO | 8 | |||||||||||
| DlCRK34 | Dlo_033618.1 | scaffold96:470428:473420 | 772 | 87.65 | 6.90 | YES | 10 | |||||||||||
| DlCRK35 | Dlo_033624.1 | scaffold96:557800:560981 | 416 | 47.66 | 5.41 | NO | 4 | |||||||||||
| DlCRK36 | Dlo_033631.1 | scaffold96:680833:683852 | 770 | 86.73 | 8.49 | YES | 12 | |||||||||||
| DlCRK37 | Dlo_033632.1 | scaffold96:704371:709405 | 963 | 108.96 | 5.70 | YES | 12 | |||||||||||
| DlCRK38 | Dlo_033633.1 | scaffold96:726239:729602 | 792 | 90.18 | 6.48 | NO | 13 | |||||||||||
| DlCRK39 | Dlo_033635.1 | scaffold96:745310:749077 | 727 | 81.75 | 6.72 | YES | 6 | |||||||||||
| DlCRK40 | Dlo_033641.1 | scaffold96:821207:825408 | 764 | 85.48 | 5.37 | YES | 10 | |||||||||||
| DlCRK41 | Dlo_001229.1 | scaffold104:636374:639626 | 793 | 89.96 | 5.01 | YES | 12 | |||||||||||
| DlCRK42 | Dlo_004774.1 | scaffold132:53035:57043 | 845 | 94.31 | 7.88 | YES | 11 | |||||||||||
| DlCRK43 | Dlo_004800.1 | scaffold132:504512:508496 | 686 | 76.82 | 7.78 | NO | 11 | |||||||||||
| DlCRK44 | Dlo_004802.1 | scaffold132:576147:580124 | 629 | 70.39 | 7.47 | YES | 10 | |||||||||||
| DlCRK45 | Dlo_004804.1 | scaffold132:603079:607481 | 844 | 94.28 | 7.19 | NO | 14 | |||||||||||
| DlCRK46 | Dlo_004805.1 | scaffold132:664877:668837 | 907 | 101.53 | 6.71 | YES | 14 | |||||||||||
| DlCRK47 | Dlo_006631.1 | scaffold149:562502:566459 | 664 | 75.09 | 7.88 | YES | 11 | |||||||||||
| DlCRK48 | Dlo_008543.1 | scaffold174:544129:545811 | 337 | 42.64 | 9.29 | NO | 2 | |||||||||||
| DlCRK49 | Dlo_008834.1 | scaffold179:1110932:1113998 | 776 | 87.91 | 5.15 | YES | 9 | |||||||||||
| DlCRK50 | Dlo_010001.1 | scaffold194:895848:898551 | 572 | 65.10 | 5.49 | NO | 4 | |||||||||||
| DlCRK51 | Dlo_011279.1 | scaffold210:240423:243279 | 696 | 78.93 | 6.41 | NO | 8 | |||||||||||
| DlCRK52 | Dlo_011280.1 | scaffold210:305315:308440 | 778 | 88.28 | 6.62 | YES | 11 | |||||||||||
| DlCRK53 | Dlo_011283.1 | scaffold210:341608:342490 | 229 | 26.09 | 8.51 | NO | 1 | |||||||||||
| DlCRK54 | Dlo_012661.1 | scaffold233:426012:431645 | 812 | 90.06 | 6.05 | YES | 14 | |||||||||||
| DlCRK55 | dlo_037774.1 | scaffold283:4344:44240 | 607 | 68.83 | 6.79 | NO | 6 | |||||||||||
| DlCRK56 | dlo_038798.1 | scaffold283:31197:44309 | 637 | 71.67 | 7.42 | YES | 6 | |||||||||||
| DlCRK57 | Dlo_018035.1 | scaffold352:67438:70571 | 666 | 74.19 | 7.80 | YES | 10 | |||||||||||
| DlCRK58 | Dlo_018037.1 | scaffold352:95501:102355 | 775 | 85.78 | 7.49 | YES | 8 | |||||||||||
| DlCRK59 | Dlo_018039.1 | scaffold352:112731:116236 | 656 | 72.96 | 8.61 | YES | 6 | |||||||||||
| DlCRK60 | Dlo_018042.1 | scaffold352:147934:154803 | 1126 | 125.45 | 7.80 | YES | 17 | |||||||||||
| DlCRK61 | Dlo_018043.1 | scaffold352:157768:160745 | 660 | 73.84 | 5.96 | YES | 9 | |||||||||||
| DlCRK62 | Dlo_018044.1 | scaffold352:182596:188678 | 616 | 68.75 | 5.53 | NO | 12 | |||||||||||
| 基因名称Gene name | 基因ID Gene ID | 位置 Location | 氨基酸个数 Number of amino aid | 分子量 Mw/ku | 等电点 pI | 信号肽 Signal peptide | 糖基化位点 N-glycosyl sites | |||||||||||
| DlCRK63 | Dlo_018045.1 | scaffold352:192448:197453 | 633 | 70.60 | 7.16 | YES | 6 | |||||||||||
| DlCRK64 | Dlo_018048.1 | scaffold352:222999:225762 | 593 | 66.58 | 8.19 | NO | 5 | |||||||||||
| DlCRK65 | Dlo_018049.1 | scaffold352:228776:232713 | 654 | 73.30 | 5.60 | YES | 8 | |||||||||||
| DlCRK66 | Dlo_018052.1 | scaffold352:263577:266686 | 637 | 71.92 | 7.15 | NO | 10 | |||||||||||
| DlCRK67 | Dlo_018053.1 | scaffold352:270544:287481 | 1103 | 121.89 | 7.16 | NO | 18 | |||||||||||
| DlCRK68 | Dlo_018696.1 | scaffold370:8789:13475 | 755 | 84.49 | 7.95 | YES | 10 | |||||||||||
| DlCRK69 | Dlo_018751.1 | scaffold372:120492:123631 | 624 | 68.54 | 7.76 | YES | 5 | |||||||||||
| DlCRK70 | Dlo_018752.1 | scaffold372:125612:129420 | 635 | 71.27 | 8.44 | YES | 3 | |||||||||||
| DlCRK71 | Dlo_019197.1 | scaffold382:209164:213134 | 651 | 71.48 | 6.44 | YES | 8 | |||||||||||
| DlCRK72 | Dlo_027776.1 | scaffold660:178964:182585 | 777 | 87.73 | 6.03 | YES | 9 | |||||||||||
| DlCRK73 | Dlo_028832.1 | scaffold707:132315:136458 | 401 | 45.12 | 9.40 | NO | 2 | |||||||||||
| DlCRK74 | Dlo_030058.1 | scaffold764:5410:9160 | 674 | 74.49 | 6.17 | YES | 9 | |||||||||||
| DlCRK75 | Dlo_030059.1 | scaffold764:24974:38644 | 1201 | 133.92 | 7.16 | YES | 14 | |||||||||||
| DlCRK76 | Dlo_030062.1 | scaffold764:77998:82065 | 675 | 73.00 | 5.47 | YES | 9 | |||||||||||
| DlCRK77 | Dlo_030064.1 | scaffold764:100652:115239 | 1152 | 129.01 | 6.04 | YES | 10 | |||||||||||
| DlCRK78 | Dlo_030104.1 | scaffold769:160423:163023 | 397 | 43.58 | 6.32 | NO | 2 | |||||||||||
| DlCRK79 | Dlo_030449.1 | scaffold784:27483:30714 | 748 | 84.64 | 6.01 | YES | 8 | |||||||||||
| DlCRK80 | Dlo_030450.1 | scaffold784:67704:79677 | 931 | 105.77 | 5.49 | YES | 11 | |||||||||||
| DlCRK81 | Dlo_031218.3 | scaffold811:54705:56903 | 523 | 58.90 | 7.15 | YES | 4 | |||||||||||
| DlCRK82 | Dlo_031222.1 | scaffold811:85102:88290 | 619 | 69.16 | 5.94 | YES | 7 | |||||||||||
| DlCRK83 | Dlo_031223.1 | scaffold811:94091:97522 | 627 | 69.64 | 6.40 | YES | 10 | |||||||||||
| DlCRK84 | Dlo_031224.1 | scaffold811:106448:109735 | 691 | 77.75 | 7.06 | YES | 10 | |||||||||||
| DlCRK85 | Dlo_031226.1 | scaffold811:136398:141354 | 652 | 72.30 | 7.09 | YES | 10 | |||||||||||
| DlCRK86 | Dlo_031227.2 | scaffold811:148465:151140 | 600 | 67.44 | 5.95 | NO | 9 | |||||||||||
| DlCRK87 | Dlo_031228.1 | scaffold811:154152:157003 | 629 | 70.46 | 5.24 | NO | 7 | |||||||||||
| DlCRK88 | Dlo_032861.1 | scaffold900:78995:83404 | 420 | 47.49 | 5.59 | NO | 3 | |||||||||||
| DlCRK89 | Dlo_033363.1 | scaffold938:54771:66654 | 900 | 101.04 | 5.68 | NO | 6 | |||||||||||
| DlCRK90 | Dlo_002474.1 | scaffold1116:5932:8771 | 616 | 68.87 | 6.48 | YES | 10 | |||||||||||
| DlCRK91 | Dlo_002475.1 | scaffold1116:11063:14292 | 420 | 47.09 | 5.96 | NO | 2 | |||||||||||
| DlCRK92 | Dlo_006918.1 | scaffold1521:2848:6164 | 665 | 74.49 | 6.58 | YES | 11 | |||||||||||
| DlCRK93 | Dlo_006923.1 | scaffold1521:37767:51591 | 1322 | 147.80 | 6.55 | YES | 20 | |||||||||||
| DlCRK94 | Dlo_009244.1 | scaffold1844:28894:32878 | 896 | 99.82 | 6.42 | YES | 9 | |||||||||||
| DlCRK95 | Dlo_013355.1 | scaffold2454:4643:7529 | 606 | 67.85 | 8.24 | NO | 7 | |||||||||||
| DlCRK96 | Dlo_013357.1 | scaffold2454:16566:20078 | 614 | 68.49 | 8.50 | YES | 14 | |||||||||||
| DlCRK97 | Dlo_013358.1 | scaffold2454:23468:26828 | 541 | 61.00 | 6.48 | NO | 7 | |||||||||||
| DlCRK98 | Dlo_013681.1 | scaffold2528:5202:8328 | 744 | 84.25 | 6.06 | NO | 9 | |||||||||||
1.2.2 CRK基因家族进化树构建 采用MEGA 6.06软件中的邻近法Neighbor-joining method对龙眼98个基因进行进化树分析,且将龙眼、拟南芥2个物种共141条CRK家族氨基酸序列构建系统进化树,手动删除同源差异较大的序列后进行Bootstrap分析,设置参数为1000次重复性检验。最后使用iTOL（itol.embl.de）对进化树进行修饰编辑。
1.2.3 龙眼CRK基因家族启动子顺式作用元件及蛋白结构域分析 本研究利用PlantCARE对龙眼CRK家族基因各成员ATG上游2000 bp序列进行分析以及利用TBtools结合龙眼CRK家族基因本身进化树进行启动子顺式作用元件分析。为了进一步了解不同龙眼CRK蛋白间motif组成的相似性和多样性,可利用MEME（http://meme）对龙眼CRK蛋白序列中的motif进行预测,其中参数为：保守基序最小宽度≥6,最大宽度为50,保守基序数目设置为25,其他参数为默认值,最后利用TBtools结合龙眼CRK家族自身进化树绘制分析图,对蛋白结构域进行预测分析。
1.2.4 龙眼CRK基因家族不同成员的不同组织器官特异性表达分析 为了解各成员在龙眼不同组织器官或在胚胎中可能发挥的功能特点,通过结合龙眼基因组数据库提取的CRK家族基因在不同体胚发生阶段NEC、EC、ICpEC、GE以及不同组织器官（种子、根、茎、叶、花、花蕾、果肉、幼果、果皮）中特异表达的FPKM值,分析龙眼CRK基因家族各成员在不同组织器官中特异性表达情况。
2 结果与分析
2.1 龙眼CRK基因家族成员的鉴定及其基本理化性质分析
根据CRK家族基因成员基因座的位置将98个成员命名为DlCRK1~DlCRK98。为进一步分析龙眼CRK基因家族的基因结构特征,了解其生物学功能,因此对龙眼98条基因的基因结构进行内含子、外显子数目及位置进行分析。结果显示,龙眼CRK基因长度介于2~39 kb,其中大部分CRK基因长度小于5 kb;外显子数目5~14个,其中大多数外显子个数为6~8个。
对龙眼DlCRK家族氨基酸序列的基本理化性质进行分析,结果如表1所示,龙眼CRK家族成员的蛋白质分子量为26.09 ku（DlCRK53）~147.8 ku（DlCRK93）,其中大部分基因蛋白分子量分布在60~90 ku之间,平均值为79.34 ku;等电点在5.01（DlCRK41）~9.49（DlCRK5）之间,平均值为6.74,其中pI值大于7的有80条,即大部分为酸性;DlCRK5、DlCRK21、DlCRK35、DlCRK43、DlCRK51、DlCRK66等39个成员不含有信号肽,剩余59个成员均含信号肽,可以引导蛋白质在细胞内运输;糖基化位点有1个（CRK53）~20个（CRK93）,大部分分布6~14个,糖基化是一种蛋白转录翻译后的修饰方式,在生命活动中起重要作用,糖基化位点的多少可能严重影响细胞的生命活动[15]。
2.2 CRK基因家族进化树构建
采用龙眼CRK氨基酸序列进行进化树构建时手动删除一个同源差异较大的基因Dlo- 011283.1（DlCRK53）,将该基因的长度、蛋白质分子量以及糖基化位点与其他基因相对比,其数值明显较低,且与其他基因的同源差异大,龙眼CRK家族基因与其他模式作物的CRK家族基因亲缘关系较远。利用龙眼和拟南芥2个物种共141条基因在手动删除3条同源差异较大的基因（Dlo-011283.1 DlCRK53、Dlo-032861.1 CRK88、AT4G23170）后构建进化树,可将龙眼和拟南芥2个物种的CRK蛋白家族分为7类（图1）：Group1中只含有龙眼CRK家族基因,这部分基因与拟南芥CRK基因同源差异较大,Group4中含有大部分ATCRK基因家族成员与DlCRK46、DlCRK63等聚类相对较近,其余各组中ATCRK基因相对较少且同物种间同源性高。
2.3 龙眼CRK基因家族启动子顺式作用元件分析
通过对龙眼CRK基因家族的启动子预测结果分析表明,龙眼CRK家族的启动子区域富含光响应、逆境胁迫响应以及植物激素响应的顺式作用元件。逆境胁迫中包括低温胁迫、干旱胁迫、创伤修复等非生物胁迫。龙眼CRK家族98个成员中可将97个进行聚类,绘制进化树与其启动子顺式作用元件结合,其中DlCRK53同源差异较大,对其单独作图（图2）。其中所有CRK基因启动子均包含至少1个光响应元件;92%的基因响应脱落酸、赤霉素、茉莉酸甲酯、生长素等激素应答;此外,86%的基因序列包含厌氧诱导所必需的元件,其中DlCRK44、DlCRK51、DlCRK97、DlCRK12这4个基因还参与缺氧特异性诱导。含有MBY转录因子是龙眼CRK基因的另一个值得注意的特点,其中72%的龙眼CRK基因包含MYB光响应以及干旱诱导的结合位点,DlCRK98、DlCRK69、DlCRK23、DlCRK31、DlCRK34中包含黄酮类化合物生物合成相关基因调控的MYB转录因子结合位点;还有部分基因（DlCRK41、DlCRK7、DlCRK67、DlCRK68、DlCRK15、DlCRK81、DlCRK31、DlCRK40）包含与细胞周期相关的顺式作用元件。总之,龙眼CRK基因家族含有多个顺式作用元件可推测其在功能上的多样性。
2.4 龙眼CRK基因家族蛋白结构域分析
为解龙眼CRK家族基因在龙眼中的蛋白保守结构域特点以及分布情况,使用MEME搜索出25个motif。大部分序列（图3）包含motif10（81）、motif6（74）、motif8（82）、motif17(73),说明这4个基序在龙眼CRK基因家族中比较保守。龙眼CRK基因具有相对保守的蛋白功能结构域（图4）,功能结构域分析显示68条基因具有STKc_IRAK结构域,只有DlCRK80和DlCRK86具有B_lecyin结构域。而不能与其他成员进行聚类的CRK53其长度较短,相应的蛋白结构域也仅有一个,其功能特性有缺失。通过对龙眼CRK蛋白家族系统发育树的分析结果可知,位于进化树同一个分支的同源基因有相似的功能结构域,进一步推测这些基因具有相似的生物学功能。
2.5 龙眼CRK基因家族不同成员在体胚发生早期不同阶段和不同组织器官特异性表达分析
为了解龙眼CRK基因家族成员们在龙眼体胚发生的不同时期,生长发育过程中不同组织器官中所起到的作用,利用龙眼全基因组数据库中获得DlCRK基因家族成员在体胚发生早期不同阶段,以及生长发育阶段不同组织器官的FPKM值,制作热图对其进行分析（图5）。除了DlCRK36、DlCRK97、DlCRK80、DlCRK88、DlCRK67这5个DlCRK成员在整个体胚发生过程以及生长发育阶段的不同器官中都不表达,其余93个DlCRK成员在体胚发生早期不同阶段或不同组织器官中都有所表达;在不同体胚发生阶段,98个成员中有22个成员在体胚发生阶段完全不表达,其余76个成员大部分在非胚性阶段（NEC）显著表达,少量成员在CE、ICPEC、GE阶段显著表达;另外,龙眼CRK家族成员在不同组织器官形态建成方面（幼果、种子、花、花蕾、叶子、果皮、茎）也发挥着重要作用,但在种子和果实阶段几乎未表达。
在不同体胚发生过程有11种表达模式：DlCRK31、DlCRK7、DlCRK46、DlCRK29、DlCRK56、DlCRK18、DlCRK22、DlCRK19、DlCRK34、DlCRK37、DlCRK35、DlCRK21、DlCRK11、DlCRK98、DlCRK79、DlCRK6、DlCRK53这17个基因在体胚发生早期的4个阶段（NEC、EC、ICpEC、GE）均未表达,不同组织器官中少量表达;在NEC阶段上调表达而其他阶段都下调表达（DlCRK25、71、76、64、1、50、2、43、85、83、55、44、66、20、89、17、40、32、8、59、77、12、61、49、74、60、45、48、5、68、93、3、96）;在EC阶段上调表达其他阶段都下调表达（DlCRK58、41、72）;在ICpEC阶段上调表达在其他阶段都下调表达（DlCRK30、47、65、9）;在NEC和EC阶段上调表达其他阶段都下调表达（DlCRK14且在EC阶段上调表达不明显）;在NEC和ICPEC阶段中表达其他阶段下调表达（DlCRK24、91、92、69、4、26、84、95、73、63）;在NEC和GE阶段上调表达其他阶段下调表达（DlCRK10、39、13）;在EC和ICPEC阶段上调表达其他阶段下调表达（DlCRK94、42、33、16、78、62、51、96、92）;在EC和GE阶段上调表达其他阶段下调表达（DlCRK28）;在ICpEC和GE阶段上调表达其他阶段下调表达（DlCRK38、54、57、27、90、75、87、73）;DlCRK81、94、15在EC、ICpEC、GE阶段均上调。总体来看龙眼CRK家族成员在体胚发生过程中主要呈下调趋势。
在生长发育过程中不同组织器官中FPKM表达情况总体而言仍是下调为主,其中种子中部分成员未表达整体呈下调趋势;果肉中只有DlCRK5显著上调以及CRK2少量上调;在幼果以及叶子中有部分家族成员显著上调表达其余均下调表达;DlCRK96、95、90在幼果、种子、花、花蕾、叶子、果皮、果实、根和茎中均未表达;DlCRK35、37、79只在叶子部分上调其余部位均未表达;DlCRK19、66、22只有在根部分上调其余部位均未表达。在不同组织器官中龙眼CRK家族成员表达情况较分散,但总体仍呈下调趋势。
3 讨论
3.1 龙眼CRK家族成员可能具有功能上的多样性
类受体蛋白激酶（receptor-like protein kinases, RLKs）定位在细胞膜上,在植物体内可以参与信号转导、抗逆反应以及病原反应等。其主要生物学功能有：参与抗病防御反应、调控植物生长与发育、参与抗非生物胁迫反应、参与胞内信号转导、调控花粉自交不亲和性[16]。而CRK基因家族（富含半胱氨酸的类受体蛋白激酶又被称为DUF26（domain of unknown function 26）类受体激酶）即类受体蛋白激酶中富含半胱氨酸的那一大类,其生物学功能即与RLK相类似。相关研究指出CRKs在植物中生物与非生物胁迫、生长发育过程以及细胞信号转导过程中起着重要的作用[17]。在其他植物中（拟南芥[3,4,5,6,7]、番茄[8]、甜瓜[11]、陆地棉[1]等）已有研究,龙眼与拟南芥通过进化树分析,同源差异较小,其生物学功能也很可能相似。在龙眼CRK基因家族中有98个成员,通过进化树分析将其中1个基因剔除（DlCRK53）后聚类可将97个成员分为7大类。其中每个成员本身又是独特的,例如在进化树中亲缘关系非常近的DlCRK8与DlCRK9,DlCRK9基因中含有MYB光响应以及与干旱诱导相关的顺式作用元件而DLCRK8基因中不含有,DlCRK8在NEC中显著表达而DlCRK9在ICpEC中显著表达;DlCRK20和DlCRK89亲缘关系非常相近,DlCRK20基因中含有控制细胞代谢周期的顺式作用元件而DlCRK89基因中不含有,DlCRK20中有20个motif而DlCRK89中有21个motif。在龙眼CRK家族中即使亲缘关系相近,其功能特性仍有差异。
龙眼CRK基因家族成员在龙眼体胚发生过程以及幼果、种子、花、花蕾、叶子、果皮、茎这些组织器官中起着重要作用。其中包含26个顺式作用元件：低温胁迫、干旱胁迫、创伤修复、缺氧特异性诱导、厌氧诱导、响应各种激素应答等。植物受渍害的本质是根系缺氧,从而影响其生长发育[18]。龙眼CRK基因家族中含有缺氧特异性诱导和厌氧诱导的顺式作用元件,推测龙眼CRK基因可以通过厌氧诱导来提高植物的耐渍能力;MYB转录因子在植物形态建成、生长发育以及初生和次级代谢产物合成等生命过程中起着重要调控作用[19],表示龙眼CRK基因家族成员在植物生长发育过程中可能起着重要作用,其中72%的龙眼CRK基因包含MYB光响应以及干旱诱导的结合位点,还有DlCRK98、DlCRK69、DlCRK23、DlCRK31、DlCRK34中包含黄酮类化合物生物合成相关基因调控的MYB转录因子结合位点,这些又暗示龙眼CRK基因可能参与龙眼生长发育过程的色素合成途径[20,21,22]。总之,98个家族成员在功能上具有多样性每个成员之间具有功能特异性,在调控植物生长发育以及在抵御逆境胁迫中发挥重要作用。
3.2 龙眼CRK家族成员参与胚胎发育和组织器官生长发育的调控
根据龙眼CRK基因家族不同成员在不组织器官特异性表达分析结果可以看出,龙眼在参与胚胎发育和组织器官生长发育的表达模式具有多样性。DlCRK家族成员在胚胎发育和组织器官生长发育中呈下调趋势,龙眼胚发生早期过程中有大量的CRK基因在NEC阶段高度表达显著,龙眼CRK基因家族成员在非胚性阶段和胚性阶段影响显著。根据观察细胞表明龙眼的胚性愈伤组织比非胚性愈伤组织细胞更加规则,并且有更多的内含物[23]。所以在龙眼中主要对其胚性愈伤组织进行研究,龙眼胚性细胞的建立与保持为龙眼研究提供便利[24]。在DlCRK41、DlCRK7、DlCRK67、DlCRK68、DlCRK15、DlCRK81、DlCRK31、DlCRK40中含有与细胞周期相关的顺式作用元件,而植物细胞周期可以调控体胚的发育[25],所以龙眼CRK家族成员可能通过调节细胞周期来参与胚胎发育的调控。
在研究中发现番茄的SlCRK1基因在花这一器官中表达具有组织特异性（具有花粉特异性的顺式作用元件）[8]。根据龙眼CRK基因家族成员在不同组织器官中特异性表达分析结果来看,部分基因在花中显著表达,所以很可能龙眼CRK家族成员在花这一器官中也具有组织特异性的表达。在DlCRK64、90、74、81中含有种子特异性的顺式作用元件,但DlCRK家族成员在龙眼种子中并未表达,表明虽含有与某种功能相关的顺式作用元件,但最终并不一定会在组织器官中表达出来。拟南芥CRK5可以通过控制植物光和作用来控制植物的生长发育[4],在龙眼CRK家族成员中有与光响应相关的顺式作用元件。对于龙眼来说与光合作用相关性最大的组织器官为叶子,与之相对应有一部分DlCRK基因在叶子中显著表达,所以龙眼CRK家族成员可能通过对植物色素合成、叶片上气孔数量[26]等因素的调控来影响植物光合作用以及抵抗光胁迫[18,19,20]。龙眼CRK基因家族成员在幼果中表达量呈上调趋势,说明龙眼CRK基因家族成员在花粉的特异性以及在幼果的形成过程中起着重要作用。龙眼CRK家族成员在胚胎与组织表达特性具有多样性与特异性,对整个龙眼生长发育过程有着重要影响。
3.3 大多数基因受植物激素和逆境调控
龙眼CRK家族成员中,大部分具有抵抗各种逆境以及受不同激素调控的顺式作用元件,其中包含低温胁迫、干旱胁迫、创伤修复等非生物胁迫;响应脱落酸、赤霉素、茉莉酸甲酯、生长素等激素应答。脱落酸在种子和芽休眠的起始和维持过程中起着重要作用[27],脱落酸可能通过调控龙眼CRK家族基因的表达来影响种子与芽的休眠过程;赤霉素能调节种子中淀粉酶的合成,来打破种子休眠,进而促进茎尖分生组织的细胞分裂,增强顶端优势,也可使果实增大[28],所以龙眼CRK家族基因可能与淀粉酶的合成相关,以此影响种子休眠以及植物生长发育。植物激素在整个植物生长发育过程中起着重要作用,而龙眼CRK家族成员大多数受多种植物激素调控,进一步说明了龙眼CRK基因家族对龙眼生长发育过程中起重要作用。
在龙眼中已发现多个基因家族中含有抵抗逆境相关的功能：龙眼漆酶家族基因可以通过影响细胞壁组成来响应非生物胁[29];DlUGD6参与龙眼中细胞壁多糖合成,可影响植物对逆境的抵抗[30];龙眼Hsf可能通过DNA甲基化来影响植物对逆境的抵抗[13];龙眼ERF家族基因与乙烯之间存在负调控[31]。那么龙眼中CRK基因家族成员抵抗逆境的机理是什么呢?蛋白激酶主要是通过激活不同磷酸化途径来介导外界环境信号的感知与传递,调控下游抗逆基因的转录表达来启动相应的生理生化等适应性的反应来降低或者消除危害[22]。有研究发现拟南芥中CRK45在植物应答非生物胁迫（干旱、高盐以及ABA胁迫）中发挥着重要作用[3],参与水杨酸介导的抗病信号途径的转变[3];AtCRK45和AtCRK36在植物细胞内共同作用,并且参与应答脱落酸和渗透胁迫信号的负调控反应[3, 10]。以上可说明龙眼CRK家族基因可能通过感知和传递信号来调控抗逆基因的表达,来启动适应逆境的能力。另外,植物激素可以对拟南芥CRK基因家族成员进行转录调控[32];而龙眼CRK家族成员中有许多响应激素应答的顺式作用元件,其必然与植物抗逆性有直接或间接相关。在植物的生长发育过程中植物激素起着重要作用,而植物抵抗逆境的能力直接决定了植物能否正常的生长发育,因此龙眼CRK家族成员对植物激素的响应与抵抗逆境的能力息息相关。