-
土壤盐渍化现象几乎遍布所有的气候区域以及地理环境,是影响农业和林业生产的主要问题。土壤盐渍化使植物受到渗透胁迫、离子毒害、膜透性改变及生理代谢紊乱等危害[1],严重制约了植物的生长,是在植物生产与栽培上一个非常严峻的问题。因此,通过对植物抗盐的分子调控机制进行探析,挖掘耐盐相关基因,并应用于新种质创制和新品种选育,对于缓解土壤盐渍化问题具有重要的现实意义。
长期生长在滨海区域潮间带的红树林,由于其特殊的、有别于陆地和淡水植物的生长环境,演化出了一套特殊的盐逆境适应机制。植物的耐盐机制涉及渗透物质积累、蛋白积累、抗氧化酶的诱导、植物激素的诱导和光合作用路径的变化5个方面[2]。对于红树植物,不同的红树物种具体的耐盐机制有所不同。无瓣海桑(Sonneratia apetala Buch.-Ham.)是海桑科海桑属红树乔木,具有生长快、耐淹浸、耐贫瘠、耐寒强等特点,是滨海滩涂的先锋造林树种,在维持海岸带生态平衡方面起着重要作用[3]。先前有关其耐盐特性的研究主要集中在形态、生理生态层面,如廖岩等[4]研究了不同盐度处理3个月下无瓣海桑根、茎、叶中的蛋白质、可溶性糖、丙二醛含量、抗氧化酶活性的变化情况,并探讨了它们的相互关系,为耐盐树种的选育提供了参考。然而,目前鲜见关于无瓣海桑耐盐分子机制方面的研究报道。
转录组测序是转录组学研究的核心技术,被广泛应用于新基因的挖掘等[5]。目前,最新的测序技术已经发展到了第三代。三代测序技术具有单分子、超长读长、无需PCR扩增等优点,已在基因组学、转录组学、表观遗传学等领域上得到了应用[6]。在转录组学上,三代全长转录组测序(Iso-Seq)的应用,可弥补二代测序技术读段较短、无法覆盖整个转录本的缺点,从而做到从全长转录本水平上研究转录组[7]。但三代测序技术存在碱基错误率高、成本高等缺点,目前应用并不广泛[8]。二代测序有通量高、成本低的特点,且发展至今已取得了较好的工作基础,因此,二代和三代测序结合使用成为很多研究者的策略[9]。如梅瑜等[10]利用基于SMRT技术原理的PacBio平台三代测序技术,对甘葛藤(Pueraria thomsonii Benth)进行全长转录组测序及分析,从转录数据中鉴定了黄酮类生物合成基因、转录因子、R基因和SSR标记,为甘葛藤的选育和利用提供了基础。然而,对于红树植物,当前研究仅限于对白骨壤(Avicennia marina (Forsk.) Vierh.)[11]、秋茄(Kandelia candel (Linn.) Druce / Kandelia obovata)[12]、角果木(Ceriops tagal (perr.) C. B. Rob)[13]、竹节树(Carallia brachiate (Lour.) Merr.)[14]和老鼠簕(Acanthus ilicifolius L.)[15]等个别物种开展了转录组学相关研究。特别是,至今对无瓣海桑仅开展了以三代全长转录组数据集为基础的盐胁迫应答转录因子挖掘工作[16]。因此,本研究利用二代与三代结合的高通量测序技术,对盐胁迫处理前后的无瓣海桑根系进行转录组测序分析,从中鉴定并初步筛选出无瓣海桑根系响应盐胁迫的候选基因,可为深入研究无瓣海桑的耐盐基因组学、分子生物学和进一步揭示红树植物的耐盐机制提供理论基础。
-
通过SOAPnuke软件对原始数据进行过滤后,如表1所示,有效读序百分率均在98%以上,Q20(碱基量 ≥ 20%)均大于96%且Q30(碱基量 ≥ 30%)均 > 91%,比对率均在87%以上。6个样品转录本reads覆盖度的峰值在90 ~ 100;空白对照RCK组和盐处理RT组的三次生物学重复间的相关系数极高(> 0.98,图1)。以上数据均说明6个样品的测序质量较高,测序数据可靠,可用于后续的差异表达分析。
表 1 高通量转录组测序数据统计表
Table 1. Statistical of high-throughput transcriptome sequencing data
样本编号
Sample ID总原始读数
Total Raw Reads /M总有效读数
Total Clean Reads /M比对率
Mapping ratio/%Q20/% Q30/% 有效读序百分率
Clean Reads Ratio/%RCK_a 43.82 43.07 87.49 96.89 91.97 98.28 RCK_b 43.82 43.26 87.89 96.92 92.03 98.71 RCK_c 43.82 43.23 87.94 96.89 91.97 98.66 RT_a 43.82 43.35 88.08 97.53 93.21 98.93 RT_b 43.82 43.23 87.62 97.78 93.88 98.65 RT_c 43.82 43.25 87.62 97.77 93.85 98.71 注:RCK,空白对照组,RT,盐胁迫处理10 d的处理组,a、b、 c表示3个生物学重复;Q20、Q30分别表示质量值≥20、30的碱基数占总碱基数的比例。 Notes: RCK represents blank control group, RT represents samples treated with salt for 10 days, and a, b, c represent the three biological repetitions; Q20, Q30, the proportion of bases with value ≥ 20, 30. -
无瓣海桑根系在500 mmol·L−1的NaCl胁迫10 d后,共检测到20 289个基因表达量发生了改变(图2)。2组样本间共有14 401个基因差异表达,其中,7 153个差异基因上调表达,log2FC最大值为11.32(isoform_13 051);7 248个差异基因下调表达,log2FC最低值为−12.69(isoform_11 580)。结果表明,下调表达的差异基因数量多于上调的数量。
-
将差异基因比对到Gene Ontology (GO)数据库,共有11 068个差异基因注释到47条GO条目,注释率为76.9%,涉及3大主要功能分类:生物过程、细胞组分和分子功能,占比分别为25.02%,43.76%,31.21%,由此可知,大多数差异基因的功能与细胞组分具有较大的相关性。根据GO功能注释分类结果(图3),这些差异基因在参与生物过程中,主要集中在细胞过程和代谢过程;在细胞组分中,主要集中在细胞、膜、膜部分;在分子功能中,主要集中在连接、催化活性。表明在高盐胁迫下,上述的亚类中涉及到的基因可能在无瓣海桑根系中响应强烈。
-
在KEGG分析结果中,共有6 189个差异基因富集到134条通路,以Q-value<0.05作为显著富集标准,共有14条显著富集的通路(图4),其中,富集差异基因数量最多的是植物激素信号转导(ko04075),涉及512个基因,最显著富集的植物-病原体相互作用(ko04626)涉及367个基因。富集分析结果表明,无瓣海桑根系基因主要通过14条显著通路参与盐胁迫的响应,并且其中涉及的基因可能对无瓣海桑根系耐盐性的调控发挥作用,关注这些显著通路,有助于对无瓣海桑根系耐盐候选基因的筛选。
-
通过进一步对差异基因进行功能注释分析,本研究共筛选出89个无瓣海桑根系抗盐候选功能基因,涉及12条KEGG通路(表2,图5)。
表 2 无瓣海桑根系部分抗盐候选功能基因
Table 2. List of partial salt-resistance candidate genes in root of S. apetala
基因ID log2FC KEGG功能注释 参与通路 活性氧清除 isoform_61079 6.46 超氧化物歧化酶,铜锌家族
superoxide dismutase, Cu-Zn family过氧物酶体
Peroxisomeisoform_173361 1.47 超氧化物歧化酶,铁锰家族
superoxide dismutase, Fe-Mn familyisoform_287367 6.10 过氧化物酶
peroxidase苯丙素生物合成
Phenylpropanoid biosynthesisisoform_57012 −4.69 isoform_25838 3.16 过氧化氢酶
catalaseMAPK信号通路-植物
MAPK signaling pathway-plantisoform_58072 −1.83 谷胱甘肽过氧化物酶
Glutathione peroxidase谷胱甘肽代谢
Glutathione metabolismisoform_67198 5.39 黄酮醇合成酶
flavonol synthase类黄酮生物合成
Flavonoid biosynthesis渗透调节 isoform_17235 7.54 海藻糖-6-磷酸磷酸酶
trehalose 6-phosphate phosphatase淀粉和蔗糖代谢
Starch and sucrose metabolismisoform_134633 −4.24 海藻糖-6-磷酸合成酶/磷酸酶
trehalose 6-phosphate synthase/phosphataseisoform_9570 3.75 水苏糖合成酶
stachyose synthetase谷胱甘肽代谢
Galactose metabolismisoform_49176 7.98 棉子糖合成酶
raffinose synthaseisoform_16906 5.82 δ-1-吡咯啉-5-羧酸合成酶
delta-1-pyrroline-5-carboxylate synthetase氨基酸的生物合成
Biosynthesis of amino acidsisoform_31419 −1.20 isoform_272673 5.53 甜菜碱醛脱氢酶
betaine-aldehyde dehydrogenase甘氨酸、丝氨酸和苏氨酸代谢
Glycine, serine and threonine metabolismisoform_228012 2.21 脯氨酸脱氢酶
proline dehydrogenase精氨酸和脯氨酸代谢
Arginine and proline metabolismisoform_8893 −2.72 植物激素 isoform_6363 8.35 生长素应答因子
auxin response factor植物激素信号转导
Plant hormone signal transductionisoform_119200 6.23 生长素应答GH3基因家族
auxin responsive GH3 gene familyisoform_126871 3.48 生长素响应蛋白
auxin-responsive protein IAAisoform_278367 4.36 SAUR家族蛋白
SAUR family proteinisoform_50365 −7.94 isoform_48396 6.73 脱落酸受体PYR/PYL家族
abscisic acid receptor PYR/PYL familyisoform_49285 −6.23 isoform_221018 −2.92 蛋白磷酸酶2C
protein phosphatase 2Cisoform_267229 4.96 脱落酸应答元件结合因子
ABA responsive element binding factorisoform_43536 −2.46 蛋白激酶 isoform_216952 4.44 促分裂原活化蛋白激酶激酶激酶1
mitogen-activated protein kinase kinase kinase1MAPK信号通路-植物
MAPK signaling pathway-plantisoform_182097 4.28 促分裂原活化蛋白激酶4
mitogen-activated protein kinase4isoform_140969 1.38 促分裂原活化蛋白激酶6
mitogen-activated protein kinase 6isoform_161961 −7.52 钙依赖性蛋白激酶
calcium-dependent protein kinase植物-病原体相互作用
Plant-pathogen interaction -
正常情况下,植物体内存在着活性氧(reactive oxygen species, ROS)产生和清除的动态平衡,但外界的不良环境会打破这个平衡,导致ROS的积累,产生毒害[22]。植物体内存在2种ROS清除机制——酶促清除和非酶促清除,前者主要涉及超氧化物歧化酶(SOD)、过氧化氢酶(CAT)等,后者主要涉及抗坏血酸、类黄酮物质等[23]。本研究发现6类共24个活性氧清除基因(图6),其中,上调表达基因数量最多的是过氧化物酶基因,总数为7,表达量差异变化最大的为铜锌家族超氧化物歧化酶基因,差异倍数达到6.46。这些上调表达的基因有利于无瓣海桑在盐胁迫前期,合成多种抗氧化酶类物质,提高体内抗氧化和ROS清除的能力,增强耐受性。
-
植物在对盐胁迫的适应过程中,细胞会积累一定量的可溶性有机物质,作为渗透调节物质参与渗透调节,以适应外界的低水势[2]。在本研究中,发现10种共22个编码生成有机渗透调节物质的基因(图6),涉及水苏糖、棉子糖、蔗糖、海藻糖、脯氨酸、甜菜碱、热激蛋白等渗透调节物质。在这些基因中,部分基因存在上调表达,表明在盐胁迫下,无瓣海桑可通过合成一些渗透调节物质,以调节盐渍环境下植物细胞内的渗透势。
-
植物体内的多种激素在其生长发育、代谢等生理过程中发挥着重要的调控作用[24]。本研究发现7种共19个涉及2种植物激素(IAA与ABA)的基因(图6),其中,IAA途径相关的有生长素应答因子(auxin response factor, ARF)、生长素应答GH3家族基因(auxin responsive GH3 gene family, GH3)、SAUR家族蛋白(SAUR family protein, SAUR)和生长素响应蛋白(auxin-responsive protein IAA, IAA)共11个基因;涉及ABA途径的有脱落酸应答元件结合因子(ABA responsive element binding factor, ABF)、脱落酸受体PYR/PYL家族(abscisic acid receptor PYR/PYL family, PYL)和蛋白磷酸酶2C(protein phosphatase 2C, PP2C)共8个基因。IAA相关基因除GH3 基因均上调表达外, 其他同时存在上下调表达,而ABA相关基因大部分下调表达。上述结果表明,植物激素能够在无瓣海桑的信号转导途径中发挥作用,不同激素共同调控盐胁迫响应过程。
-
蛋白激酶可调控植物的生长,在逆境中起到调控感受胁迫信号、启动各种非生物逆境响应的作用[25]。 本研究发现6种共10个蛋白激酶相关基因(图6),其中,2种促分裂原活化蛋白激酶基因(mitogen-activated protein kinases, MAPK)MPK4、MPK6和2种丝氨酸/苏氨酸蛋白激酶基因PBS1、CTR1均上调表达,1种促分裂原活化蛋白激酶基因(mitogen activated protein kinase, MEKK)MEKK1和1种钙依赖性蛋白激酶基因(calcium-dependent protein kinase, CDPK)中均有2个上调表达,1个下调表达。可见这些蛋白激酶在无瓣海桑根系响应高盐胁迫时,发挥着尤为重要的调控作用。
-
转录因子是一类功能蛋白质,具有调控基因表达的作用。当植物受到逆境胁迫时,转录因子通过与相应的顺式作用元件相互作用,调控并减轻逆境胁迫带来的伤害,以适应各种逆境[26]。本研究发现8种共14个差异表达的转录因子(图6),涉及5个家族,分别为WRKY家族、MYB家族、GRAS家族、bHLH家族和EREBP家族,其中,WRKY家族发现的成员数量最多,暗示WRKY家族在无瓣海桑响应盐胁迫中起着重要的调控作用。
无瓣海桑根响应盐胁迫的转录组分析
Transcriptome Analysis of Sonneratia apetala Root in Response to Salt Stress
-
摘要:
目的 初步探究无瓣海桑的耐盐分子机制,筛选出无瓣海桑抗盐候选基因,为后续功能验证实验及林木抗盐性遗传育种奠定分子基础。 方法 以1年生无瓣海桑幼苗为材料,用500 mmol·L−1 NaCl分别处理0 d(对照组)和10 d(处理组),取不同条件下的根部组织进行转录组测序,并结合三代全长转录组数据进行后续生物信息学分析。 结果 (1)与对照组相比,NaCl处理10 d后,无瓣海桑幼苗根系中共有14401个差异表达基因,其中,7153个上调,7248个下调。(2)GO分析发现,共有11068个差异基因在47个GO 条目得到注释。(3)在KEGG富集分析中,共有6189个差异基因富集到134条通路,其中,共有14条通路显著富集(P值<0.01,Q值<0.05)。(4)通过进一步对差异基因进行功能注释分析,共筛选出抗盐候选基因89个,其中,抗氧化基因24个,渗透调节物质基因22个,植物激素基因19个,蛋白激酶基因10个,转录因子基因14个。 结论 活性氧清除、渗透调节、植物激素、蛋白激酶及转录因子相关基因参与调控无瓣海桑盐逆境适应过程。 Abstract:Objective To lay a molecular basis for functional verification experiment and genetic breeding of tree salinity tolerance, the genetic mechanisms was explored, and salt-related genes were identified for Sonneratia apetala Buch.-Ham. Method In this study, the root tissues were collected from a 1-year-old S. apetala treated with 500 mmol•L−1 NaCl for 0 (control group) and 10 d (treatment group). Then transcriptome sequencing and bioinformatics analysis was performed based on the three-generation full-length transcriptome dataset of S. apetala. Result (1) Compared with the control group, 14401 genes were differentially expressed after salt treatment, of which 7153 were up-regulated and 7248 were down-regulated. (2) GO analysis found that a total of 11068 differential genes were annotated in 47 GO items. (3) For KEGG enrichment analysis, a total of 6189 differential expression genes were enriched to 134 pathways, of which 14 were significantly enriched (P-value <0.01, Q-value <0.05). (4) Further functional annotation analysis of the differentially expressed genes revealed a total of 89 genes was potential salt-related candidate genes. Among these, 24, 22, 19, 10 and 14 genes encoded enzymes or functional proteins referred to antioxidation, osmotic adjustment substances, plant hormones, protein kinase and transcription factors, respectively. Conclusion Genes relating to active oxygen scavenging, osmotic regulation, plant hormones, protein kinases and transcription factors participate in the regulation of salt stress adaptation in S. apetala. -
Key words:
- Sonneratia apetala
- / salt stress
- / transcriptome
- / differentially expressed genes (DEGs)
-
表 1 高通量转录组测序数据统计表
Table 1. Statistical of high-throughput transcriptome sequencing data
样本编号
Sample ID总原始读数
Total Raw Reads /M总有效读数
Total Clean Reads /M比对率
Mapping ratio/%Q20/% Q30/% 有效读序百分率
Clean Reads Ratio/%RCK_a 43.82 43.07 87.49 96.89 91.97 98.28 RCK_b 43.82 43.26 87.89 96.92 92.03 98.71 RCK_c 43.82 43.23 87.94 96.89 91.97 98.66 RT_a 43.82 43.35 88.08 97.53 93.21 98.93 RT_b 43.82 43.23 87.62 97.78 93.88 98.65 RT_c 43.82 43.25 87.62 97.77 93.85 98.71 注:RCK,空白对照组,RT,盐胁迫处理10 d的处理组,a、b、 c表示3个生物学重复;Q20、Q30分别表示质量值≥20、30的碱基数占总碱基数的比例。 Notes: RCK represents blank control group, RT represents samples treated with salt for 10 days, and a, b, c represent the three biological repetitions; Q20, Q30, the proportion of bases with value ≥ 20, 30. 表 2 无瓣海桑根系部分抗盐候选功能基因
Table 2. List of partial salt-resistance candidate genes in root of S. apetala
基因ID log2FC KEGG功能注释 参与通路 活性氧清除 isoform_61079 6.46 超氧化物歧化酶,铜锌家族
superoxide dismutase, Cu-Zn family过氧物酶体
Peroxisomeisoform_173361 1.47 超氧化物歧化酶,铁锰家族
superoxide dismutase, Fe-Mn familyisoform_287367 6.10 过氧化物酶
peroxidase苯丙素生物合成
Phenylpropanoid biosynthesisisoform_57012 −4.69 isoform_25838 3.16 过氧化氢酶
catalaseMAPK信号通路-植物
MAPK signaling pathway-plantisoform_58072 −1.83 谷胱甘肽过氧化物酶
Glutathione peroxidase谷胱甘肽代谢
Glutathione metabolismisoform_67198 5.39 黄酮醇合成酶
flavonol synthase类黄酮生物合成
Flavonoid biosynthesis渗透调节 isoform_17235 7.54 海藻糖-6-磷酸磷酸酶
trehalose 6-phosphate phosphatase淀粉和蔗糖代谢
Starch and sucrose metabolismisoform_134633 −4.24 海藻糖-6-磷酸合成酶/磷酸酶
trehalose 6-phosphate synthase/phosphataseisoform_9570 3.75 水苏糖合成酶
stachyose synthetase谷胱甘肽代谢
Galactose metabolismisoform_49176 7.98 棉子糖合成酶
raffinose synthaseisoform_16906 5.82 δ-1-吡咯啉-5-羧酸合成酶
delta-1-pyrroline-5-carboxylate synthetase氨基酸的生物合成
Biosynthesis of amino acidsisoform_31419 −1.20 isoform_272673 5.53 甜菜碱醛脱氢酶
betaine-aldehyde dehydrogenase甘氨酸、丝氨酸和苏氨酸代谢
Glycine, serine and threonine metabolismisoform_228012 2.21 脯氨酸脱氢酶
proline dehydrogenase精氨酸和脯氨酸代谢
Arginine and proline metabolismisoform_8893 −2.72 植物激素 isoform_6363 8.35 生长素应答因子
auxin response factor植物激素信号转导
Plant hormone signal transductionisoform_119200 6.23 生长素应答GH3基因家族
auxin responsive GH3 gene familyisoform_126871 3.48 生长素响应蛋白
auxin-responsive protein IAAisoform_278367 4.36 SAUR家族蛋白
SAUR family proteinisoform_50365 −7.94 isoform_48396 6.73 脱落酸受体PYR/PYL家族
abscisic acid receptor PYR/PYL familyisoform_49285 −6.23 isoform_221018 −2.92 蛋白磷酸酶2C
protein phosphatase 2Cisoform_267229 4.96 脱落酸应答元件结合因子
ABA responsive element binding factorisoform_43536 −2.46 蛋白激酶 isoform_216952 4.44 促分裂原活化蛋白激酶激酶激酶1
mitogen-activated protein kinase kinase kinase1MAPK信号通路-植物
MAPK signaling pathway-plantisoform_182097 4.28 促分裂原活化蛋白激酶4
mitogen-activated protein kinase4isoform_140969 1.38 促分裂原活化蛋白激酶6
mitogen-activated protein kinase 6isoform_161961 −7.52 钙依赖性蛋白激酶
calcium-dependent protein kinase植物-病原体相互作用
Plant-pathogen interaction -
[1] 胡 涛, 张鸽香, 郑福超, 等. 植物盐胁迫响应的研究进展[J]. 分子植物育种, 2018, 16(9):3006-3015. [2] 牛 恋. 植物耐盐生理机制及抗盐性[J]. 河北农机, 2018(2):56. [3] 欧阳少明. 无瓣海桑的造林技术和管理措施[J]. 农业科技与信息, 2016(15):139-139. doi: 10.3969/j.issn.1003-6997.2016.15.116 [4] 廖 岩, 赵 肖, 陈桂珠. 盐胁迫对无瓣海桑幼苗根茎叶膜保护系统的影响[J]. 海洋环境科学, 2009, 28(2):154-158. doi: 10.3969/j.issn.1007-6336.2009.02.010 [5] 高慧娟, 吕昕培, 王润娟, 等. 转录组测序在林草植物抗逆性研究中的应用[J]. 草业学报, 2019, 28(12):184-196. doi: 10.11686/cyxb2019352 [6] 赵陆滟, 曹绍玉, 龙云树, 等. 全长转录组测序在植物中的应用研究进展[J]. 植物遗传资源学报, 2019, 20(6):1390-1398. [7] 李玉梅, 李书娴, 李向上, 等. 第三代测序技术在转录组学研究中的应用[J]. 生命科学仪器, 2018, 16(C1):114-121,113. [8] 周卫星, 石海鹤. 高通量测序中序列拼接算法的研究进展[J]. 计算机科学, 2019, 46(5):36-43. doi: 10.11896/j.issn.1002-137X.2019.05.005 [9] 马东娜, 张兴坦, 魏柳锋, 等. 基因组二代测序数据与三代测序数据的混合校正和组装[J]. 基因组学与应用生物学, 2018, 37(4):1547-1555. [10] 梅 瑜, 李向荣, 蔡时可, 等. 药食同源植物甘葛藤的全长转录组分析[J]. 华北农学报, 2021, 36(5):10-17. doi: 10.7668/hbnxb.20192213 [11] 郝赛琦. 红树植物白骨壤根系结构与转录组特征研究[D]. 厦门: 厦门大学, 2017: 11-57. [12] SU W, YE C, ZHANG Y, et al. Identification of putative key genes for coastal environments and cold adaptation in mangrove Kandelia obovata through transcriptome analysis[J]. Science of the Total Environment, 2019, 681: 191-201. doi: 10.1016/j.scitotenv.2019.05.127 [13] YANG Y, YANG S, LI J, et al. De novo assembly of the transcriptomes of two yellow mangroves, Ceriops tagal and C. zippeliana, and one of their terrestrial relatives, Pellacalyx yunnanensis[J]. Marine Genomics, 2015, 23: 33-36. doi: 10.1016/j.margen.2015.04.003 [14] GOU W, WU H, ZHANG Z, et al. Comparative analysis of transcriptomes in Rhizophoraceae provides insights into the origin and adaptive evolution of mangrove plants in intertidal environments[J]. Frontiers in Plant Science, 2017, 8: 795. doi: 10.3389/fpls.2017.00795 [15] YANG Y, YANG S, LI J, et al. Transcriptome analysis of the Holly mangrove Acanthus ilicifolius and its terrestrial relative, Acanthus leucostachyus, provides insights into adaptation to intertidal zones[J]. BMC Genomics, 2015, 16(1): 605. doi: 10.1186/1471-2164-16-1 [16] 刘婷婷, 莫玉剑, 欧成川, 等. 盐胁迫下无瓣海桑差异表达转录因子的转录组信息分析[J]. 分子植物育种, 2022, 20(10): 3210-3222. [17] CHEN Y, CHEN Y, SHI C, et al. SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data[J]. GigaScience, 2018, 7(1): gix120. [18] LANGMEAD B, SALZBERG S L. Fast gapped-read alignment with Bowtie 2[J]. Nature Methods, 2012, 9(4): 357-359. doi: 10.1038/nmeth.1923 [19] VUONG H, TRUONG T, TRAN T, et al. A revisit of RSEM generative model and its EM algorithm for quantifying transcript abundances[J]. bioRxiv, 2018: 503672. [20] LOVE M I, HUBER W, ANDERS S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biology, 2014, 15(12): 550. [21] KANEHISA M, ARAKI M, GOTO S, et al. KEGG for linking genomes to life and the environment[J]. Nucleic Acids Research, 2007, 36(suppl_1): D480-D484. [22] 吴 寒. 活性氧在植物体内的作用及其清除体制[J]. 广东蚕业, 2018, 52(3):18. [23] 郭明欣, 刘佳佳, 侯琳琳, 等. 植物体内活性氧的产生及清除机制研究进展[J]. 科技视界, 2021(8):104-106. [24] 黎 家, 李传友. 新中国成立70年来植物激素研究进展[J]. 中国科学:生命科学, 2019, 49(10):1227-1281. [25] 朱婷婷, 王彦霞, 裴丽丽, 等. 植物蛋白激酶与作物非生物胁迫抗性的研究[J]. 植物遗传资源学报, 2017, 18(4):763-770. [26] WANG Y, XU W, CHEN Z, et al. Gene structure, expression pattern and interaction of Nuclear Factor-Y family in castor bean (Ricinus communis)[J]. Planta, 2018, 247(3): 559-572. doi: 10.1007/s00425-017-2809-2 [27] ZHU J K. Abiotic stress signaling and responses in plants[J]. Cell, 2016, 167(2): 313-324. [28] 付晴晴, 李商锐, 刘小瑞, 等. NaCl胁迫对不同耐盐性葡萄株系根系活性氧代谢的影响[J]. 中外葡萄与葡萄酒, 2017(3):12-15. [29] 鲁 琳, 杨尚谕, 刘维东, 等. 基于转录组测序花烟草响应盐胁迫活性氧清除相关基因的挖掘[J]. 生物技术通报, 2020, 36(12):42-53. [30] KOSAR F, AKRAM N A, SADIQ M, et al. Trehalose: a key organic osmolyte effectively involved in plant abiotic stress tolerance[J]. Journal of Plant Growth Regulation, 2019, 38(2): 606-618. doi: 10.1007/s00344-018-9876-x [31] 李 辉, 李德芳, 邓 勇, 等. 红麻海藻糖生物合成关键酶基因HcTPPJ的克隆及响应逆境的表达分析[J]. 作物学报, 2020, 46(12):1914-1922. [32] 张莹莹, 吕 燕, 宣雯燕, 等. 龙须菜(Gracilariopsis lemaneiformis)海藻糖-6-磷酸合成酶(TPS)对逆境胁迫的响应[J]. 海洋与湖沼, 2021, 52(3):777-785. doi: 10.11693/hyhz20200700225 [33] ZHANG H, DONG J, ZHAO X, et al. Research progress in membrane lipid metabolism and molecular mechanism in peanut cold tolerance[J]. Frontiers in Plant Science, 2019, 10: 838. doi: 10.3389/fpls.2019.00838 [34] FUNCK D, BAUMGARTEN L, STIFT M, et al. Differential contribution of P5CS isoforms to stress tolerance in Arabidopsis[J]. Frontiers in Plant Science, 2020, 11: 1483. [35] GONG Z, XIONG L, SHI H, et al. Plant abiotic stress response and nutrient use efficiency[J]. Science China Life Sciences, 2020, 63(5): 635-674. doi: 10.1007/s11427-020-1683-x [36] 刘海洋, 陈玉珍. ABA受体结构及功能与ABA信号通路研究进展[J]. 中国农学通报, 2019, 35(21):75-81. doi: 10.11924/j.issn.1000-6850.casb19010079 [37] 张 弛, 蔚静玲, 储谟立, 等. 拟南芥蛋白磷酸酶PP2C31的盐胁迫响应功能研究[J]. 中国科技论文, 2018, 13(18):2070-2075. doi: 10.3969/j.issn.2095-2783.2018.18.004 [38] ZHANG M, SU J, ZHANG Y, et al. Conveying endogenous and exogenous signals: MAPK cascades in plant growth and defense[J]. Current Opinion in Plant Biology, 2018(45): 1-10. [39] ABULFARAJ A A. Stepwise signal transduction cascades under salt stress in leaves of wild barley (Hordeum spontaneum)[J]. Biotechnology & Biotechnological Equipment, 2020, 34(1): 860-872. [40] KUMAR K, RAINA S K, SULTAN S M. Arabidopsis MAPK signaling pathways and their cross talks in abiotic stress response[J]. Journal of Plant Biochemistry and Biotechnology, 2020, 29(4): 700-714. doi: 10.1007/s13562-020-00596-3 [41] 汪芳珍, 杨成行, 何子华, 等. 盐处理下旱生植物沙芥蛋白激酶相关基因的差异表达分析[J]. 草业学报, 2021, 30(10):116-124. doi: 10.11686/cyxb2020368 [42] 潘凌云, 马家冀, 李建民, 等. 植物盐胁迫应答转录因子的研究进展[J]. 生物工程学报, 2022, 38(1):50-65. [43] DU C, MA B, WU Z, et al. Reaumuria trigyna transcription factor RtWRKY23 enhances salt stress tolerance and delays flowering in plants[J]. Journal of Plant Physiology, 2019, 239: 38-51. [44] LIN J, DANG F, CHEN Y, et al. CaWRKY27 negatively regulates salt and osmotic stress responses in pepper[J]. Plant Physiology and Biochemistry, 2019, 145(C.): 43-51. doi: 10.1016/j.plaphy.2019.08.013