盐胁迫下7个刺榆无性系的叶片转录组分析

摘要：[目的]从转录组水平上探究刺榆叶片对盐胁迫的响应机制，筛选耐盐基因。[方法]本研究以7个不同耐盐能力刺榆无性系的扦插苗为试验材料，对0、100、150、200 mmol·L-1 NaCl溶液处理15 d的叶片进行了转录组测序与分析。[结果]（1）除DW1上调数量略多于下调数量外，其他无性系核心差异基因下调数量均多于上调数量，且各无性系在200 mmol·L-1盐溶液处理时差异表达基因数量最多，对盐胁迫的响应最为强烈。Go富集分析显示细胞壁等功能类别基因对盐胁迫的响应较强。KEGG富集分析显示植物激素信号转导等代谢通路在7个无性系响应盐胁迫时发挥了重要作用。（2）强耐盐型的无性系DJ1、DJ3与弱耐盐型的无性系DW1、DY1在200 mmol·L-1盐溶液处理时的4个对比组合共有的差异基因数量为1 406个，共同上调基因393个、共同下调基因996个，其中脱落酸调节因子蛋白磷酸酶PP2C基因AHG1、ABI2，乙烯受体基因ERS1、ETR2，以及环核苷酸门控离子通道基因C～GC10和C～GC/5可能与DJ1、DJ3较强的耐盐能力有关。[结论]盐胁迫引起刺榆转录水平上功能类别基因与代谢通路的响应，其中6个显著差异表达基因将作为候选耐盐基因。本研究为揭示刺榆响应盐胁迫的转录组特征与筛选耐盐基因提供了理论参考，为进一步开展刺榆种质资源的研究与应用奠定了基础。

关键词：刺榆；盐胁迫；差异表达基因；实时荧光定量PCR；耐盐基因

中图分类号：S722 文献标识码：A 文章编号：1001-1498（2024）05-0136-12

盐渍土在我国分布广泛，其类型多、面积大，且在不断扩张，严重威胁生态安全与农业可持续发展。诸多研究表明，植物具有改良土质、改善生态环境、实现有效治理与开发利用盐渍化土地的重要途径，开展耐盐植物种质资源的挖掘、筛选与应用研究对耐盐植物栽培利用和盐碱地改良具有重要意义。植物在受到盐胁迫时首先会在基因和转录组水平上进行调控，进而控制代谢以应对逆境。转录组学研究技术可以从植物整体转录水平揭示植物体中复杂的生物学通路和调控网络分子，对研究植物响应盐胁迫的复杂机制具有重要意义，目前在金丝楸（Catalpa bungei C.A.Mey）、多枝柽柳（Tamarix ramosissima Ledeb.）、平邑甜茶（Malus hupehensis Rehd.）等树种响应盐胁迫机理研究中得到较广泛的应用。

刺榆（Hemiptelea davidii （Hance） Planch.）为落叶小乔木，榆科（Ulmaceae）刺榆属（Hemiptelea）单种，耐瘠薄、耐低温、适应性强，在我国西北、西南、东北及华东等地分布较广泛，其中在山东主要分布在济南南部山区、淄博鲁山、胶东山区丘陵等地，近年来的调查发现在刘公岛等近海岛屿中也有分布。目前，刺榆的研究主要集中在种苗繁育与生态作用方面，抗逆研究报道较少。其中，耐盐性研究报道仅见杨颖对刺榆种子萌发及其幼苗进行的盐胁迫研究，该研究显示在NaCl＜30 mmol·L-1的盐胁迫下，刺榆种子能够正常发芽，幼苗内MDA含量和SOD、POD酶活性未发生显著变化，表明刺榆种子萌发与幼苗生长对低浓度盐环境具有一定的耐受性，为刺榆在盐碱地的种植和推广提供了理论依据。刺榆响应胁迫的分子机制研究报道亦较少，无性系响应盐胁迫的转录组学特征与耐盐基因筛选目前尚未见报道，刘华波利用干旱模拟试验和转录组数据探讨了刺榆核质基因响应干旱胁迫的分子机制，研究表明核质基因中光合作用相关通路、植物激素信号转导通路和次生代谢物合成通路受干旱胁迫影响最大，筛选并初步证实TrxG-like基因参与刺榆干旱胁迫响应、提高植株的抗旱性。本研究以7个刺榆无性系为试验材料，对不同浓度NaCl溶液处理下各无性系叶片的转录组学特征开展研究，以期为刺榆对盐胁迫的响应机制研究及耐盐基因的筛选与利用提供参考。

1 材料与方法

1.1 试验材料与处理

试验于山东省林草种质资源中心温室中进行。7个供试刺榆无性系（DJ1、DJ2、DJ3、DJ4、DW1、DY1、DH1）为2021年春季采集野生母株1年生枝条进行扦插繁育得到的无性系，其中DJ1～DJ4采自济南市柳埠街道锦阳川、DW1采自潍坊仰天山、DY1采自烟台昆嵛山、DH1采自威海刘公岛，前期研究表明DJ1与DJ3为强耐盐型，DJ2为较强耐盐型，DH1与DJ4为中等耐盐型，DW1与DY1为弱耐盐型。

试验苗以泥炭土、蛭石、珍珠岩按10：4：1的体积比混匀配制的营养土为基质，栽植于直径18 cm、高度15 cm的花盆中。2022年7月，将根系发育良好、长势基本一致的扦插苗置于温室内，室内温湿度与自然条件基本一致：试验期间实测夜间温度19～22℃、相对湿度73%～90%，日温23～32℃、相对湿度50%～74%。

2022年8月，参照杨颖的工作及在前期预试验的基础上，设置清水（0 mmol·L-1）为对照，50、100、150、200 mmol·L-14种浓度NaCL溶液。为避免产生盐激效应，在试验开展前先用100 mmol·L-1 NaCl溶液预处理150、200 mmol·L-1两个梯度的试验苗3d，预处理后与其他浓度一起开始盐胁迫处理。处理时，花盆置于托盘上，每盆分2次浇施1L相应浓度盐溶液，浇施的盐溶液从花盆自然渗出到托盘中，保留渗出的盐溶液使其浸泡花盆1夜（12 h）以使栽培基质充分吸收盐溶液，次日取出后用塑料薄膜包裹花盆减少盆口水分蒸发，每5d处理1次，各处理重复3次，共处理105株。试验处理至15 d时，各无性系叶片出现不同程度的盐害症状，即进行盐害情况调查，并取相同叶位的叶片，去离子水冲洗后用滤纸吸干水分，采集后立即放到液氮中速冻，然后置于-80℃冰箱中保存，用于转录组测序。

1.2 指标测定方法

1.2.1 盐害情况调查

参照白榆耐盐性评价研究中对盐害程度分级的方法，根据刺榆叶片损伤症状，将刺榆盐害分为0、Ⅰ、Ⅱ、Ⅲ4个等级（图1），盐害等级评价标准为：0级：叶片无盐害症状；Ⅰ级：叶缘轻度萎蔫，萎蔫面积不足1/5；Ⅱ级：叶缘卷曲干枯，叶片萎蔫面积1/5～1/2：Ⅲ级：叶缘卷曲干枯，叶片萎蔫面积1/2以上至干枯脱落。

1.2.2 转录组测序和分析

使用Trizol试剂（Invitrogen，CA，USA）提取样本总RNA，利用Bioanalyzer 2100和RNA 1000 Nano LabChipKit（Agilent，CA，USA）对提取的样本总RNA质量进行检测，质检合格后，利用附着有poly-T寡核苷酸的磁珠纯化mRNA，然后使用Fragmentation Buffer对其进行片段化，利用mRNASeq sample preparation kit（Illumina，SanDiego，USA）逆转录裂解的RNA片段创建最终的cDNA文库，文库质检合格后采用IlluminaNovaseqTM6000进行测序，测序读长为双端2x150 bp。

在去除测序接头、低质量等不合格序列后，使用FastQC验证序列质量，包括有效数据的Q20、Q30和GC含量，所有的下游分析均基于高质量的有效数据。使用Trinity 2.4.0进行转录组的从头组装，归一化后得到非重复基因序列（unigenes）。使用DIAMOND将组装完的unigenes与六大数据库（NCBI_NR、GO、Swiss-Prot、KEGG、Pfam、eggNOG）进行比对、分类和相应的功能注释，E-value阈值小于0.000 01。

通过Salmon计算样品基因TPM值以评估unigenes的表达水平。使用R包edgeR对样本之间进行显著差异分析，将llog2 fold changel＞1，且P值＜ 0.05的unigenes作为差异表达基因，并对其进行GO和KEGG富集分析。

1.2.3 实时荧光定量PCR分析

随机选择4个差异表达基因进行实时荧光定量PCR分析，进一步验证转录组数据的可靠性。利用Trizol法提取样品的总RNA，使用NovoScript Plus All-in-one 1stStrand cDNA Synthesis SuperMix（gDNAPurge）反转录试剂盒将RNA反转录为cDNA。以较为稳定的GAPDH作为内参基因，按照TBGreen⑩Premix Ex TaqTM（Tli RNaseH Plus）荧光定量试剂盒的操作方法，进行实时荧光定量PCR，测定样品中相关基因的表达量，按照2-△△CT法计算基因的相对表达量，基因扩增引物序列见表1。

2 结果与分析

2.1 盐胁迫对各株系叶片表型的影响

如图2所示，在50 mmol·L-1浓度处理下各无性系均未表现出盐害症状；在100 mmol·L-1盐溶液处理时，各无性系均表现为Ⅰ级盐害症状；在150 mmol·L-1盐溶液处理时，各无性系进一步表现出受害症状，且受盐害程度出现较大分化，其中DJ1、DJ2、DJ3、DJ44个无性系盐害程度为Ⅰ级，DW1、DH1盐害程度为Ⅱ级，DY1无性系受害情况最重、达到Ⅲ级；在200 mmol·L-1盐溶液处理时，各无性系健康分值均为处理间最低，不同无性系的受害程度也进一步分化，其中DJ1无性系、DJ3无性系盐害程度虽然加重但仍为I级，DJ2、DJ4、DH1 3个无性系盐害程度为Ⅱ级，DW1、DY1两个无性系受害情况最严重、达到Ⅲ级。

2.2 测序质量

测序得到的原始序列（Raw Reads）数量介于39 993 532条～44 393 618条，共1 169 469262条，过滤后得到的有效序列（Valid Reads）数量介于37 318 336条～43 091 940条，共1 123 282 434条，有效序列占比介于92.28%～97.07%，平均占比96.05%；原始碱基（RawBases）数量6.00 G～6.66 G，共175.42 G，有效碱基（Valid Bases）数量5.19 G～5.99 G，共156.03 G；Q20介于96.64%～97.14%，平均97.01%，Q30介于90.82%-91 .88%，平均91.62%；GC含量介于43.61%～44.61%，平均44.14%。有效数据、Q20、Q30及GC含量表明测序质量可靠，可进行深层次分析。

2.3 基因组装结果

将获得的有效数据使用Trinity 2.4.0进行组装，共组装转录本174 852条、碱基192 45681 1个，其长度介于201～16 721 bp，中值为715 bp，200～500 bp占比为3g%，GC含量为39.97%，N50为1799 bp；归一化后共得到基因87 712条、碱基78 334 647个，长度范围与转录本相同，中值为483 bp，200～500 bp占比达51%，GC含量为40.03%，N50为1590 bp，组装结果较好。

2.4 基因功能注释结果

使用DIAMOND对组装好的Unigenes进行功能注释，6种数据库的注释结果统计见表2，各数据库注释的基因数量由多到少依次为eggNOG、GO、Pfam、Swissprot、NR、KEGG，总注释87712个基因。

2.5 差异表达基因分析

2.5.1 各无性系核心差异表达基因筛选与分析

2.5.1.1 各无性系核心差异表达基因筛选

为明确获得刺榆各株系响应盐胁迫的差异表达基因，对各无性系中与对照组差异明显的100、150、200 mmol·L-13个处理浓度（依次标记为2、3、4）与对照0 mmol·L-1（标记为1）进行比较，以llog2 fold changel＞1且P值＜0.05为标准筛选差异表达基因，将3个对比组共有的差异基因作为该无性系核心差异表达基因。

如图3（A-G）所示，在DJ1株系共筛选出703个核心差异基因，其中在3个处理浓度对比组中均上调的基因184个、均下调的514个；DJ2株系共筛选出430个核心差异基因，其中在3个处理浓度对比组中均上调的基因154个、均下调的263个；DJ3株系共筛选出1 021个核心差异基因，其中在3个处理浓度对比组中均上调的基因178个、均下调的830个；DJ4株系共筛选出960个核心差异基因，其中在3个处理浓度对比组中均上调的基因166个、均下调的780个；DW1株系共筛选出702个核心差异基因，其中在3个处理浓度对比组中均上调的基因351个、均下调的319个；DY1株系共筛选出1 757个核心差异基因，其中在3个处理浓度对比组中均上调的基因629个、均下调的1 097个；DH1株系共筛选出323个核心差异基因，其中在3个处理浓度对比组中均上调的基因108个、均下调的194个。除DW1株系上调基因数量略多于下调基因数量外，其他无性系下调基因数量均多于上调基因数量。此外，各无性系的3个对比组中差异表达的基因数量均为200 mmol·L-1VS对照0 mmol·L叫的对比组最多，表明各无性系在200 mmol·L-1高浓度盐溶液处理下受到的胁迫压力最重，对盐胁迫的响应也最强烈。

2.5.1.2 各无性系差异表达基因富集分析

如图4（A-G）所示，GO功能富集分析显示DJ1株系核心差异表达基因功能主要被显著富集在细胞壁、细胞外区、细胞壁组织、钙调素结合、对过氧化氢的反应、质膜等类别。DJ2株系核心差异表达基因功能主要被显著富集在叶绿体类囊体膜、叶绿体包膜、类囊体、叶绿体类囊体、细胞壁、叶绿体等类别。DJ3株系核心差异表达基因功能主要被显著富集在细胞壁、细胞外区、细胞壁组织、植物型细胞壁、细胞分裂、胞间连丝等类别。DJ4株系核心差异表达基因功能主要被显著富集在细胞壁、DNA结合、细胞外区、核、细胞分裂、DNA结合转录因子活性等类别。DW1株系核心差异表达基因功能主要被显著富集在细胞外泌体、膜、细胞壁、核糖体的结构成分、细胞外区、蛋白质结合等类别。DY1株系核心差异表达基因功能主要被显著富集在DNA结合转录因子活性、质膜、细胞外区、转录与DNA模板、蛋白质丝氨酸／苏氨酸激酶活性、转录调控与DNA模板等类别。DH1株系核心差异表达基因功能主要被显著富集在对缺水的反应、对过氧化氢的反应、细胞外区、序列特异性DNA结合、细胞壁等类别。其中，细胞壁、细胞外区、细胞壁组织、植物型细胞壁、细胞壁生物发生、膜锚定构件、质膜锚定构件、α-L-阿拉伯呋喃糖苷酶活性、DNA结合转录因子活性在7个无性系中均被显著富集，表明上述功能类别基因，尤其是与细胞壁和细胞膜相关的基因对盐胁迫的响应较强。

如图5（A-G）所示，KEGG富集分析显示DJ1株系核心差异表达基因主要被显著富集在淀粉和蔗糖代谢、植物激素信号转导、苯丙烷生物合成、戊糖和葡萄糖醛酸相互转化、植物MAPK信号通路、氧化磷酸化等通路。DJ2株系核心差异表达基因主要被显著富集在光合作用、植物MAPK信号通路、植物激素信号转导、苯丙烷生物合成、戊糖和葡萄糖醛酸相互转化、植物与病原体的相互作用等通路。DJ3株系核心差异表达基因主要被显著富集在戊糖和葡萄糖醛酸相互转化、淀粉和蔗糖代谢、植物激素信号转导、氨基糖与核苷酸糖代谢、半乳糖代谢、苯丙烷生物合成等通路。DJ4株系核心差异表达基因主要被显著富集在植物激素信号转导、DNA复制、植物MAPK信号通路、半乳糖代谢、淀粉和蔗糖代谢、嘧啶代谢等通路。DW1株系核心差异表达基因主要被显著富集在核糖体、氧化磷酸化、吞噬体、内质网中的蛋白质加工、植物激素信号转导、淀粉和蔗糖代谢等通路。DY1株系核心差异表达基因主要被显著富集在植物激素信号转导、苯丙烷生物合成、植物MAPK信号通路、淀粉和蔗糖代谢、戊糖和葡萄糖醛酸相互转化、半乳糖代谢等通路。DH1株系核心差异表达基因主要被显著富集在植物MAPK信号通路、淀粉和蔗糖代谢、半乳糖代谢、植物激素信号转导、氨基糖与核苷酸糖代谢、核糖体等通路。其中，植物激素信号转导、淀粉和蔗糖代谢、苯丙烷生物合成、植物与病原体互作、植物MAPK信号通路在各无性系中富集基因数量较多且在多个无性系中呈显著富集，表明盐胁迫对上述代谢通路产生了较大影响。

2.5.2 强耐盐型与弱耐盐型的无性系之间差异表达基因筛选与分析

由2.5.1可知，各无性系在最高浓度200 mmol·L-1盐溶液处理下对盐胁迫响应最强烈、差异表达基因最多，为进一步了解耐盐性好的刺榆无性系响应盐胁迫的机制及筛选耐盐基因，对强耐盐型无性系DJ1、DJ3与弱耐盐型无性系DW1、DY1在最高浓度200 mmol·L-1盐溶液处理时的叶片转录组数据进行比较分析。

如图3（H）所示，DJ14 vs DW14有6 248个差异显著表达基因，3 209个上调、3 039个下调；DJ14 vs DY14有2 954个差异表达基因，1 057个上调、1 897个下调；DJ34 vs DW14有6 930个差异显著表达基因，3 704个上调、3 226个下调；DJ34 vs DY14有3 913个差异表达基因，1 687个上调、2 226个下调；4个对比组中共有的差异基因数量为1 406个，其中共同上调基因393个、共同下调基因996个。

对4个对比组共有的1 406个差异基因进行富集分析，如图4（H） GO富集分析所示，差异基因主要富集在细胞外泌体、DNA结合转录因子活性、膜、转录，DNA模板、转录调控-DNA模板、质膜、细胞壁等功能类别。图5 （H） KEGG富集分析显示，差异基因主要富集在植物激素信号转导、淀粉和蔗糖代谢、内质网中的蛋白质加工、半乳糖代谢、植物MAPK信号通路、苯丙烷生物合成等代谢通路。

植物激素信号转导通路在2.5.2各无性系差异基因分析和本节强耐盐型与弱耐盐型的无性系间差异表达基因分析中均为主要富集通路，且富集基因数量较多，反映其与刺榆响应盐胁迫有较大关系，本研究对其在本节富集的22个差异基因进行分析。如图6代谢通路图所示，生长素信号转导途径中富集6个基因，其中上调节点中2个基因注释为生长素响应蛋白IAA、1个基因注释为生长素响应蛋白SAUR，下调节点中1个基因注释为生长素响应因子ARF、2个基因注释到GH3生长素响应基因家族；细胞分裂素信号转导途径富集2个下调基因，分别注释到细胞分裂素受体CRE1和转录因子ARR-B；赤霉素信号转导途径富集的1个上调基因注释为赤霉素受体GID1，富集的2个下调基因均注释为DELLA蛋白；脱落酸信号转导途径富集到5个下调基因，均注释为蛋白磷酸酶PP2C;乙烯信号转导途径富集的2个下调基因均注释为乙烯受体ETR；油菜素内酯富集的2个上调基因分别注释到油菜素内酯受体BRI1和转录因子CYCD3，富集的1个下调基因注释为木葡聚糖：木葡糖基转移酶TCH4；茉莉酸信号转导途径富集到1个下调基因，注释为转录因子MYC2。上述途径中相关基因的差异表达表明多个激素参与了DJ1与DJ3对盐胁迫的响应，且这些差异基因的表达可能与两个无性系有较好的耐盐性有关。

此外，在对4个样本中基因表达量的分析中发现一个在耐盐性好的无性系中大幅表达而在耐盐性差的无性系中表达量较低的基因C～GC15，注释显示其为植物环核苷酸门控通道（cyclic nucleotidegated channel，plant），KEGG分析显示该基因在植物与病原体互作途径中通过调控Ca2+进而影响下游防御响应相关基因表达，其大幅表达可能与刺榆DJ1和DJ3株系有较好的耐盐性有关。

2.6 差异表达基因qRT-PCR验证

如图7所示，TRINITY_DN61310_c1_g2、TRINITY_DN38595_c1_g1、TRINITY DN63193c0_g2、TRINITY_DN65673_c0_g24个基因的qRT-PCR结果与转录组测序数据趋势一致、相关性较强，表明本研究的测序结果具有可靠性。

3 讨论

通过转录组测序得到87 712条基因，对各无性系及处理浓度之间进行比较显示，除DW1株系差异基因中下调基因数量略低于上调基因数量外，其他无性系下调差异基因数量均多于上调基因数量，表明在刺榆响应盐胁迫的过程中下调基因发挥了重要作用。盐胁迫会引起植物细胞壁与细胞膜发生变化，本研究中各无性系核心差异基因GO富集分析亦显示与细胞壁和细胞膜相关的功能类别在7个无性系均被富集到。KEGG富集分析显示植物激素信号转导、淀粉和蔗糖代谢、苯丙烷生物合成、植物与病原体互作、植物MAPK信号通路在各无性系中富集基因数量较多且在多个无性系中呈显著富集，与其他学者对金丝楸（C．bungeiC.A.Mey）、无瓣海桑（Sonneratia apetalaBuchanan-Hamilton）等植物响应盐胁迫的研究结果相似，表明上述通路在7个刺榆无性系响应盐胁迫时发挥了重要作用。

对耐盐性好的DJ1、DJ3株系与耐盐性差的DW1、DY1株系在高盐胁迫下的叶片转录组进行比较，共有的差异基因数量为1 406个，其中共同上调基因393个、共同下调基因996个。对4个对比组共有的1 406个差异基因进行富集分析，GO富集分析显示，差异基因主要富集在细胞外泌体、DNA结合转录因子活性、膜、转录与DNA模板、转录调控与DNA模板、质膜、细胞壁等功能类别；KEGG富集分析显示，差异基因主要富集在植物激素信号转导、淀粉和蔗糖代谢、内质网中的蛋白质加工、半乳糖代谢、植物MAPK信号通路、苯丙烷生物合成等代谢通路，GO富集分析与KEGG富集分析表明200 mmol·L-1高浓度盐溶液胁迫下DJ1与DJ3株系有更强的调节信号转导、蛋白质、碳水化合物、次生代谢物与维护细胞壁与细胞膜结构和功能的能力，上述过程对DJ1和DJ3株系响应盐胁迫、提高抵御盐害的能力具有重要作用。

本研究中多个差异基因被注释到生长素、赤霉素、细胞分裂素、油菜素内酯、脱落酸、乙烯等植物激素信号途径，研究表明这些激素不仅可以调节植物生长发育，在植物响应盐胁迫的过程中也发挥着重要的作用。其中脱落酸（ABA）是植物响应盐胁迫非常重要的应激激素，在ABA信号通路中包含多种关键的调节因子，蛋白磷酸酶2C（protein phosphatases 2C，PP2C）作为核心调节因子之一，以负调控的方式发挥作用，其机理为PP2C与ABA正调控因子SnRK2蛋白激酶结合后抑制了其酶活性，阻止了其对下游转录因子的激活，盐胁迫下，PP2C活性会被抑制从而激活ABA信号传导，调控植株耐盐性，李梦霜对盐碱胁迫下茶树（Camellia sinensis（L．）O．Ktze．）ABA信号转导系统PP2C基因家族的研究结果显示CsPP2CA6呈现显著的下调表达。本研究中强耐盐型无性系相较于弱耐盐型无性系的差异基因中AHG1、ABI2等5个注释为PP2C的基因均显著下调，与前述研究结果一致。AHG1、AB12是PP2C家族中的A组成员，作为核心负调控因子在ABA信号传导通路中发挥着重要作用，这2个基因可能与DJ1株系和DJ3株系有较好的耐盐性相关。乙烯是一种气态植物激素，在植物响应盐胁迫过程中同样起着非常重要的调控作用。张驶等阐述了在乙烯信号转导途径中ETR1、ERS1、ETR2等乙烯受体可与下游的蛋白激酶CTR1结合阻断乙烯信号传导，即以负反馈形式控制乙烯信号的输出，在对强耐盐型无性系与弱耐盐型无性系的差异基因分析中发现2个乙烯受体基因ERS1和ETR2表达量均显著下调，这2个基因表达量的下调可能会增强乙烯信号传导，进而影响了DJ1株系与DJ3株系对盐胁迫的的耐受能力。此外，注释为CNGC15的基因在强耐盐型无性系中大幅上调表达而在弱耐盐型无性系中表达量较小，其为环核苷酸门控离子通道（CNGCs）家族成员，CNGCs是植物细胞中非选择性阳离子通道和分子开关，参与植物的生长发育，抵御生物和非生物胁迫。Saand等以番茄（Solanum lycopersicum L．）沉默株为材料对C～GC15基因的进行了研究，结果显示其与抗旱性有关，与耐盐性相关性较小，但Saand也阐述了该基因在沉默株中有表达，认为需要进一步研究C～GC15基因在植物耐盐胁迫中的作用，在本研究中C～GC15基因在耐盐性好的DJ1株系、DJ3株系大幅上调表达，而在耐盐性较差的DY1株系、DW1株系中表达量较低，可能与刺榆DJ1株系、DJ3株系较好的耐盐性机制有关。CNGC10是CNGCs家族另一重要成员，Guo等、Jin等的研究均显示C～GC10对拟南芥（Arabidopsis thaliana （L.） Heynh.）的耐盐性具有负调控作用，在本研究中，C～GC10在强耐盐型无性系中表达量较弱耐盐型无性系显著下调，根据上述研究推测，差异基因C～GC10亦可能与DJ1株系、DJ3株系较好的耐盐能力相关。

4 结论

盐胁迫引起刺榆细胞壁、细胞外区、DNA结合转录因子活性等功能类别基因与植物激素信号转导、淀粉和蔗糖代谢、苯丙烷生物合成等代谢通路的响应，强耐盐型无性系较弱耐盐型无性系对上述过程具有更强的调节能力，其中6个显著差异表达基因将作为候选耐盐基因。本研究为揭示刺榆响应盐胁迫的转录组特征与筛选耐盐基因提供了理论参考，为进一步开展刺榆种质资源的研究与应用奠定了基础。

（责任编辑：张研）

基金项目：山东省重点研发计划（重大科技创新工程）项目“珍贵用材树种种质资源挖掘与精准鉴定”（2021LZGC02304）