大豆根系响应早期缺铁胁迫的转录组分析

[目的]本文旨在探明2个大豆品种根系对早期缺铁胁迫的响应差异.[方法]以铁高效品种'吉育99'和铁低效品种'吉育93'大豆为材料进行水培试验,将对照(25μmol?L-1 Fe)和早期(1和6 h)缺铁胁迫(1μmol?L-1 Fe)处理的大豆根系进行转录组测序.[结果]与对照相比,铁高效大豆品种差异基因数随着缺铁胁迫...     

萱草根茎低温胁迫转录组分析

低温胁迫是萱草育种和生产中最严重的非生物胁迫之一.该研究以耐冷型和冷敏感型萱草低温处理前后的根茎为材料,利用转录组测序(RNA-seq)技术,比较了耐冷型和冷敏感型萱草在低温胁迫下的基因表达差异.结果表明,4个样品获得的高质量clean reads均大于35880734条;低温处理后,耐冷型萱草差异表达基因数量多于冷敏...     

铜胁迫下小麦幼根转录组学及蛋白质组学研究

为探讨外源铜胁迫对小麦幼根的影响,采用水培方法,测定0~60 mg/L不同质量浓度铜处理下小麦幼苗生长发育、转录组表达和蛋白质组表达的差异,探索其抵御重金属毒害的分子机制。结果表明:高质量浓度的铜影响小麦生长,随胁迫时间的延长,小麦幼苗萎蔫发黄,幼根生长受抑。小麦幼根转录组测序检测,共得到2 283个差异表达基因(DEGs),其中826个DEGs表达上调,1 457个DEGs表达下调;将DEGs进行KEGG pathway注释,被注释到31个代谢途径中。双向电泳试验检测,小麦幼根中有2 049个蛋白质点,30 mg/L铜处理96 h后,与对照(0 mg/L处理)相比,约130个蛋白质点表达出现显著差异,其中56个丰度增加,74个丰度降低;质谱鉴定部分差异表达蛋白,抗性蛋白如谷胱甘肽转移酶、27K蛋白等在胁迫下表达上升,而生理代谢相关蛋白表达降低。     

基于转录组学的植物响应盐胁迫调控机制研究进展

盐胁迫是限制植物生长和产量的重要环境因子之一.经过长期的进化,植物已形成了一套响应盐胁迫的调控机制.转录组学可以从植物mRNA整体转录水平揭示植物响应盐胁迫的调控机制,对研究植物抗盐、耐盐具有重要意义.本文针对转录组学在植物响应盐胁迫调控机制中的研究,简述了植物体内的信号传导、渗透调节、内源激素合成、光合作用、活性氧清...     

基于转录组分析筛选牛心朴子响应低温胁迫的转录因子家族

[目的]从转录水平分析牛心朴子在低温胁迫下的差异表达基因,筛选响应低温胁迫的转录因子家族,鉴定出牛心朴子低温胁迫响应的关键调控基因,为全面解析逆境胁迫响应分子调控网络及有效挖掘关键调控基因提供理论参考.[方法]通过高通量测序技术对低温胁迫(CT组)和常温处理(对照,CK组)的牛心朴子cDNA文库进行转录组测序分析,对差异表达基因进行功能注释和富集,鉴定出响应低温胁迫的转录因子家族,并选取4个差异表达的转录因子基因进行实时荧光定量PCR检测,以验证转录组测序结果的可信度.[结果]从CK组和CT组共获得30.50 Gb的原始数据,Cycle Q20平均值在96%以上,经数据过滤及去冗余后,拼接组装获得100006条Unigenes,但其功能注释率较低,有47082条至少在一个数据库中被功能注释,占Unigenes总数的47.07%;而在NR、NT、KO、SwissProt、PFAM、GO和KOG数据库中均被注释的Unigene有7070条,仅占Unigenes总数的7.06%.GO功能富集分析筛选到5545个差异表达基因,其中,上调表达基因2039个,下调表达基因3506个,分别富集到生物过程、细胞组分和分子功能三大类别中.从牛心朴子Unigene中共鉴定到83个转录因子家族的1826个转录因子,其中,以MYB转录因子家族成员数目最多,为136个(占7.45%).从差异表达基因中筛选到与低温胁迫有关的66个转录因子家族的550个转录因子,其中MYB、C3H、bHLH、AP2-EREBP、C2H2、NAC、bZIP、CCAAT和WRKY等转录因子家族均有大量转录因子能被低温胁迫诱导表达.基于实时荧光定量PCR的牛心朴子低温胁迫下转录因子基因表达水平检测结果与转录组测序分析结果基本一致.[结论]MYB、C3H、bHLH、AP2-EREBP、C2H2、NAC、bZIP、CCAAT和WRKY等转录因子家族成员在牛心朴子响应低温胁迫时发挥主导作用,同时各家族转录因子间存在共表达性或协同作用,通过复杂的转录调控网络发挥重要调节作用,进而提高牛心朴子对低温胁迫的耐受性.     

小麦条锈菌产孢转录组分析

为研究不同毒力小麦条锈菌侵染过程中基因表达差异,通过RNA-seq测定其转录组并分析差异表达基因。结果表明,每个小麦条锈菌转录组文库产生8 000万～9 000万条读长,其中,至少71.15%的读长可以比对到参考基因组。将读长转换计算出FPKM值并比较多个基因表达水平,FPKM值小于15的基因数分别占对应文库的79.17%、79.24%、80.01%和79.41%,FPKM值大于60的基因数所占比例均低于7.0%,而FPKM值大于10 000的基因在CY31、CY32、CY33和V26菌系样品文库中分别有7、8、7、7个。从|log₂Fold change|≥1且FDR<0.005这2个水平评估小麦条锈菌CY31、CY32、CY33和V26菌系间差异表达基因,结果显示,差异水平变化主要集中在6倍以内,差异倍数最大的一个基因为13倍。将差异表达基因与KEGG数据库比对发现,共有53个基因被注释到了KEGG中的36个通路中,其中,显著富集通路(P<0.05)有5个,主要参与到氧化磷酸化、淀粉和糖代谢等通路中。可见,小麦条锈菌CY31、CY32、CY33和V...     

波纹龙虾对亚硝酸盐胁迫响应的转录组学分析

为了探究波纹龙虾Panulirus homarus对亚硝酸盐的胁迫响应机理,采用转录组学测序分析方法,进行了正常养殖条件下波纹龙虾(体质量为39.00 g±5.31 g)肝胰腺和质量浓度为80 mg/L亚硝酸盐胁迫7 d的肝胰腺转录组比较研究。结果表明：转录组共筛选得到264个差异表达基因,其中,86个基因上调表达,178个基因下调表达;GO功能注释显示,被注释的差异表达基因主要与结合、催化和代谢等功能有关;KEGG通路富集分析显示,差异表达基因在黏着斑、白细胞跨内皮层迁移、胰岛素抵抗、紧密连接、淀粉与蔗糖代谢等通路中显著富集;随机挑选9个差异表达基因进行qRT-PCR验证,发现荧光定量结果与测序结果基本一致。研究表明,在亚硝酸盐胁迫下,通过GO和KEGG分析发现差异表达基因主要富集在能量代谢、生长和免疫等通路中,挖掘出的数据可为后续研究提供科学参考。     

不同桑树品种响应干旱胁迫的比较转录组学分析

【目的】发掘响应干旱胁迫的关键抗旱基因,从转录水平上揭示桑树的抗旱分子机制,为后续开展桑树分子抗旱性育种工作提供科学依据。【方法】以干旱敏感性品种德果1号和耐旱性品种湖桑32号为研究对象,通过盆栽种植方式开展干旱胁迫及复水处理试验,采集24个桑叶样本提取总RNA后构建cDNA文库,在Illumina HiSeqTM 4000测序平台上进行高通量测序,并结合生物信息学对相关基因进行注释分析。【结果】经转录组测序,各样本的Clean reads范围为45723096~67280168条,有效碱基数（Clean bases）集中在6.86~10.09 Gb,GC含量在44.84%~46.48%（平均为45.61%）,Q30在90.36%~93.07%。差异表达基因（DEGs）筛选结果显示,6个差异分组（A2 vs A1,B2 vs B₁,A3 vsA1,B3 vs B₁,A4 vs A1,B4 vs B₁）分别筛选出3510、3399、5677、5507、5124和2734个差异表达基因;经干旱胁迫处理后,桑叶功能组基因中呈下调表达的差异表达基因明显多于呈上调表达的差异表达基因,说明桑树在生长过程中存在不同的功能基因以控制桑叶生长发育。不同抗旱性桑叶转录组测序数据中以涉及生物学过程的差异表达基因最多,且主要集中在小分子代谢过程（Small molecule metabolic process）、跨膜运输（Transmembrane transport）及碳水化合物代谢过程（Carbohydrate metabolic process）等方面;与分子功能相关的差异表达基因次之,主要涉及转移酶活性（Transferase activity）和水解酶（Hydrolase activity）等。干旱胁迫下不同抗旱性桑叶差异表达基因主要富集到59条KEGG信号通路上,可划分为代谢、遗传信息处理、环境信息处理、细胞过程和生物系统五大信号通路;其中,抗旱性桑树品种通过提高能量代谢、碳水化合物代谢、增强光合作用及脂质代谢来更好地适应干旱胁迫。综合GO功能注释分析和KEGG信号通路富集分析,得到以下可能与干旱相关的基因: LOC21410404、LOC21404884、LOC21409623、LOC21401352、LOC21398977、LOC21388561、LOC21401447、LOC21398764、LOC21385254、LOC21384661、LOC-21410404及LOC21408971。【结论】不同桑树品种的耐旱性存在显著差异,其中抗旱性桑树品种是通过提高能量代谢、碳水化合物代谢、脂质代谢及光合作用共同应对干旱胁迫。     

铜胁迫下的西瓜食酸菌转录组分析

【目的】分析西瓜食酸菌(Acidovorax citrulli)在铜胁迫下转录水平的响应特征,研究对其抗铜系统和抗铜机制,为进一步解析西瓜食酸菌的铜耐受分子机制积累数据。【方法】以西瓜食酸菌亚群I的抗铜菌株FC440铜胁迫和无胁迫培养下的菌体为材料,提取RNA并利用RNA-Seq技术获取细菌对铜胁迫的应答数据,分析数据质量、差异表达基因的类型、KEGG富集与Pathway注释、GO富集分析等转录组水平特征,并进行qPCR验证。【结果】共获得6个转录组文库,总数据量23.6 Gb;该细菌共有4 344个基因表达,其中有466个基因差异表达(上调266个,下调200个);差异表达的基因显著富集于ABC转运系统的跨膜转运、双组分系统通路的信号转导以及氨基酸代谢等约9大通路;基因表达趋势与转录组测序数据的趋势一致。【结论】铜胁迫下谷氨酸、谷胱甘肽等氨基酸代谢在西瓜食酸菌抗铜胁迫中发挥重要作用;双组分系统参与西瓜食酸菌对铜离子胁迫的响应且受到负调控;该菌中可能存在包含ABC转运体的跨膜转运、氧化还原反应、胞内束缚沉淀等多种系统互作的抗铜机制。     

基于RNA-Seq技术分析植物激素信号途径在水稻幼苗中对低温胁迫的应答规律

水稻对低温敏感,幼苗时期的低温胁迫会影响水稻的生长,甚至致死,从而导致减产。植物激素在植物低温胁迫响应过程中发挥着重要作用。在本研究中,我们以籼稻品种特青和粳稻品种02428为试验材料,在不同低温处理时间下(10℃处理0、3或24 h),运用RNA-Seq技术开展水稻苗期低温应答转录组分析。结果表明:粳稻02428的耐低温性显著强于籼稻特青,低温处理后,在两个品种中均检测到大量差异表达基因。在差异表达基因的KEGG通路分析中,植物激素信号通路是低温胁迫下差异表达基因最显著富集的通路,说明植物激素通路在水稻苗期耐低温中具有重要作用。我们通过对生长素、细胞分裂素等8种激素信号途径中差异表达基因的数量及其表达量进行分析,解析了各激素信号通路中的关键基因在低温胁迫下表达量的变化趋势,最终发现在水稻苗期耐低温中可能起正调控作用的激素是乙烯、脱落酸和茉莉酸,负调控耐低温的是细胞分裂素和赤霉素。本研究为深入了解植物激素参与水稻苗期耐低温调控的机制提供了方向与依据。     

干旱胁迫下割手密根系转录组差异表达分析

【目的】探明云南割手密82-114受干旱胁迫的内在分子机制,挖掘与抗旱密切相关的功能基因,提高割手密抗旱基因的研究与利用。【方法】以干旱胁迫24、48和72 h的云南割手密82-114的根系进行Illumina Hi Seq TM 4000高通量转录组测序分析,将组装得到的Unigene分别与Swiss-Prot、Nr、KOG、Pfam和KEGG数据库进行比对,利用RPKM来衡量各样本间的基因表达丰度,以FDR≤0.05和|log2 fold change|≥1来评估样本间的差异表达基因,通过Gene Ontology(GO)数据库、KEGG pathway数据库对不同干旱胁迫样本差异表达基因的功能和参与的调控路径进行分析。【结果】未胁迫处理的CK与胁迫24、48和72 h后的样本转录组分别检测到134 724、130 368、133 564和131 321个表达基因,与CK相比,各胁迫样本显著上调(或下调)的差异表达基因分别有3 061(1 302)、2 304(2 841)和3 236(2 525)个;以FDR值=0为筛选条件,3个样本的极显著差异表达基因主要参与转录激活、水分运输、DNA结合、ATP结合及细胞膜、跨膜运输和防御反应等有关代谢活动。GO富集分析表明:3个胁迫处理的样本在生物学途径中富集最多的类别并不完全相同,其中,胁迫24 h与48 h的样本富集最多的类别是与DNA依赖型转录和转录调节子相关的基因,胁迫72 h样本富集最多的类别是与翻译相关的基因;而在细胞组件和分子功能中,3个胁迫样本富集最多的基因类别均是与内在的膜、核和ATP结合相关的基因。KEGG富集分析表明:胁迫24 h的样本有2 248个差异表达基因参与了107条代谢途径,胁迫48 h的样本有2 114个差异表达基因参与了130条代谢途径,胁迫72 h的样本有2 392个差异表达基因参与了144条代谢途径;经KEGG显著性富集分析,筛选到鞘糖脂生物合成途径、MAP激酶信号途径和ABC转运蛋白等与非生物胁迫或逆境相关。【结论】获得3个干旱胁迫样本与CK样本间差异表达基因的变化及其功能信息,发现参与云南割手密82-114干旱胁迫响应相关的细胞外信号调节激酶、磷脂酶D、β-氨基己糖苷酶和ATP结合盒B亚族(MDR/TAP)-1等,获得在整个干旱胁迫时期稳定上调表达基因70个,稳定下调表达基因11个,通过同源基因序列的比对分析,阐明这些基因在各自的代谢通路中被强烈诱导或抑制表达,显示与干旱胁迫存在密切关系。     

甘蓝型油菜响应干旱胁迫的转录组分析

甘蓝型油菜(Brassica napus L.)是我国主要油料作物之一,为国家粮油安全提供了重要保障。干旱严重限制了油菜种植规模扩大、产量与品质提升。为挖掘甘蓝型油菜抗旱的关键基因,阐明油菜耐旱机制,本研究以甘蓝型油菜湘油15为对象,通过实验室模拟干旱处理后的叶片为材料进行转录组测序(RNA-Seq)、生物信息分析并qRT-PCR验证。发现转录组样本平均Clean baseas约为9.6 G,且各样品碱基百分比(Q30)大于93.21%,比对到参考基因组上的Reads在89%以上。干旱处理4 h后,湘油15出现了440个差异表达基因(DEGs),其中148个基因上调,292个基因下调。随机选取10个基因(上调与下调各5个)进行qRT-PCR验证,其中9个基因的qRT-PCR结果与RNA-Seq结果相符,仅1个基因表达趋势相反,相似度可达90%,证实转录组测序结果的可靠性。结果表明甘蓝型油菜湘油15主要通过光合作用、碳代谢和渗透调节等途径调节抗旱能力,为进一步阐明油菜耐旱机制奠定了基础。     

猴樟叶片转录组测序及盐碱胁迫响应基因表达模式研究

为了解猴樟响应盐碱胁迫机制,以2个耐盐碱猴樟品系和1个盐碱敏感型猴樟品系3年生种苗叶片为材料,利用Illumina Hiseq^TM平台对其进行转录组测序,对6个精氨酸和脯氨酸通路相关基因进行qRT−PCR分析。结果表明：研究共获得72369条Unigenes,其中39156条Unigenes在NR、NT、KEGG、SwissProt、PFAM、GO、KOG数据库中获得注释;HDHZ vs. HJHZ 组有5415个DEGs,其中4003个得到GO注释;YYHZ vs. HJHZ 组有4665个DEGs,其中331个得到GO注释。此外,KEGG通路分析发现,HDHZ vs. HJHZ组上调基因主要富集在氨基酸代谢和生物碱合成相关通路,下调基因主要富集在倍半萜和三萜的生物合成,以及光合作用相关通路;YYHZ vs. HJHZ组上调基因主要富集在倍半萜和三萜的生物合成及类黄酮生物合成通路;下调基因主要富集在色素代谢和酚类物质代谢相关通路。结合6个差异表达基因的qRT−PCR结果,生物碱类、精氨酸和脯氨酸代谢可能在猴樟响应碱胁迫过程中起主要作用,所得基因的表达模式与RNA−Seq数据结果一致,证实RNA−Seq数据准确。     

野生罗非鱼响应低温胁迫的脑组织转录组测序分析

【目的】比较低温胁迫下野生罗非鱼和养殖吉富罗非鱼的基因表达模式,揭示野生罗非鱼低温应答的特有分子调控机制,为后续开展野生罗非鱼耐冷亲本大范围筛选打下基础。【方法】挑选规格相近的野生罗非鱼和养殖吉富罗非鱼进行梯度降温试验,水温先以1℃/12 h的速度从26℃降至14℃,并保持288 h（12 d）,于26℃、20℃及14℃保持0、120和288 h等5个时间点分别解剖罗非鱼采集脑组织样品,构建cDNA文库后以Illumina HiSeq×Ten测序平台进行双端测序及比较转录组分析,并采用实时荧光定量PCR对具有重要功能的差异表达基因（DEGs）进行表达验证。【结果】野生罗非鱼在11℃时开始出现死亡个体,但在8℃时仍有50.0%的存活个体,说明野生罗非鱼较养殖吉富罗非鱼具有更强的低温耐受能力。在14℃的长时间低温胁迫过程中,养殖吉富罗非鱼脑组织中表达显著上调的差异表达基因数量约是野生罗非鱼的10倍,即养殖吉富罗非鱼的应激反应远比野生罗非鱼强烈。KEGG信号通路富集分析结果显示,养殖吉富罗非鱼和野生罗非鱼在14℃保持120和288 h的上调差异表达基因均富集到核糖体发生、内质网蛋白加工及剪接体信号等信号通路上,野生罗非鱼的上调差异表达基因还额外富集到核苷酸切除修复、N-糖基化生物合成和DNA复制信号等通路上。与养殖吉富罗非鱼相比,在野生罗非鱼中发现577个特有的差异表达基因,主要富集在NOD受体信号通路、凋亡和内吞等通路上,且表现为NOD受体信号通路被启动,而细胞凋亡受抑制。在野生罗非鱼中,参与NOD受体信号通路的关键功能基因（Nemo、NFκB、p38、JNK和IL-1β）在长时间低温胁迫中均维持在一个较高的表达水平,而内吞途径关键基因的表达上升倍数也明显高于养殖吉富罗非鱼。【结论】野生罗非鱼通过避免过度的应激反应和维持低水平代谢的细胞稳定以减轻低温胁迫对机体的损伤,并持续启动NOD受体信号通路及内吞途径以维持其免疫能力,而具有较养殖罗非鱼更强的低温耐受力。     

大豆PP2C家族基因鉴定与响应盐胁迫的转录组分析

植物蛋白磷酸酶2C(PP2C)能够通过影响植物体内多种生物学过程在环境胁迫响应中发挥重要作用.为研究栽培大豆PP2C家族成员在盐胁迫中的功能,通过生物信息学手段结合转录组学分析对栽培大豆PP2C家族进行了系统探究.在栽培大豆中共鉴定出126个PP2C家族成员,并将其在进化树上分为11个亚家族,同一亚家族的成员具有类似的...     

花生幼苗响应低温胁迫的转录组分析

低温冷害是限制我国东北地区花生产量和品质提升的主要环境因素,鉴定花生耐寒的关键功能基因,并揭示其调控机制是创制耐寒花生种质,提高花生耐冷性的重要途径之一。以耐寒型花生品种农花5号为试验材料,对6℃低温胁迫下的花生叶片进行RNA-Seq转录组测序分析。结果表明:与对照相比,共有3910个基因持续差异表达;GO功能注释将差异表达基因富集到43个功能组,其中生物学过程中的代谢过程富集到的基因数目最多;KEGG代谢通路分析将持续差异表达基因富集到116个代谢通路中,其中植物激素信号转导、淀粉和蔗糖代谢、植物-病原菌互作、植物昼夜节律和α-亚麻酸代谢途径富集到的基因数目最多,且大部分为上调表达,在低温胁迫中具有重要作用。从转录组数据库中随机挑选6个差异表达基因进行qRT-PCR验证,其表达趋势与高通量测序结果相一致,证明了转录组数据的可靠性。研究结果将为揭示花生响应低温胁迫的分子机制提供理论依据。     

干旱胁迫棉花转录组DNA损伤修复相关基因的分析

【目的】研究干旱胁迫下,棉花DNA损伤修复基因的表达情况,从整体水平探讨棉花DNA损伤修复相关基因表达与抗旱的相关性。【方法】用2.5%PEG6000处理棉花幼苗,利用RNA-Seq技术对干旱胁迫下棉花幼苗的转录组进行测序。【结果】从棉花干旱胁迫响应的转录组中,筛选获得差异表达的DNA损伤修复相关基因共51个,其中差异表达的上调基因23个,差异表达的下调基因28个,干旱胁迫能够影响棉花DNA损伤修复相关基因的表达。选取4个差异基因进行生物信息学分析及qRT-PCR验证,HMGB1、rec A1、UDGs和GMP synthase基因的表达变幅有一定的差异,但基因的表达趋势一致,棉花转录组测序结果通过qRTPCR验证是可靠的。【结论】DNA损伤修复相关基因可能与干旱胁迫有一定的相关性,干旱胁迫能够影响棉花DNA损伤修复相关基因的表达。     

盐胁迫下栽培大豆与野生大豆根系离子流及转录组分析

为揭示栽培大豆和野生大豆盐胁迫下根系离子流变化及内在分子响应机制,挖掘大豆耐盐基因,以普通栽培大豆(Glycine max, Gm)和滩涂野生大豆(Glycine soja, Gs)为试验材料,通过设置100 mmol·L^-1NaCl盐胁迫处理,用非损伤微测技术(NMT)检测大豆根系离子流对盐胁迫的响应,并用高通量测序技术对栽培大豆及野生大豆盐胁迫24 h后根系的转录组进行分析,以寻找离子转运体相关的耐盐基因。结果表明：与对照相比,100 mmol·L^-1NaCl胁迫10 d后栽培大豆鲜重下降41.5%,而野生大豆鲜重下降29.2%。根系离子流测定结果显示,盐胁迫24 h后野生大豆比栽培大豆Na⁺外排量高46.3%, K⁺损失量少33.0%,且H⁺吸收量高93.9%。转录组分析结果显示,野生大豆阳离子/H⁺反向转运体CHX20和钠钾根缺陷NaKR2显著下调,K⁺/H⁺逆向转运体KEA2、KEA3和K^+<...     

海马齿根系响应盐胁迫的转录组分析

【目的】对盐胁迫下海马齿根系进行转录组测序分析,挖掘海马齿根系耐盐相关基因,为揭示海马齿耐盐的分子机制提供参考。【方法】利用Illumina测序技术对0 mmol/L NaCl （对照组）和400 mmol/L NaCl胁迫处理（盐胁迫处理组）下的海马齿根系进行转录组测序分析,从中筛选出差异表达基因,选取13个基因进行实时荧光定量PCR （qRTPCR）检测,以验证转录组数据的可靠性。【结果】在海马齿根系转录组中共鉴定出305145个转录本,平均长度为622 bp,其中,对照组有146177个长度>300 bp的转录本,盐胁迫处理组有72173个长度>300 bp的转录本;共有65535条Unigenes在Nr、GO、Swiss-Prot、COG和KEGG五大数据库注释成功,占Unigenes总数的52.36%。对照组和盐胁迫处理组共有65535个差异Unigenes,其中,有182个热休克蛋白基因。对照组和盐胁迫处理组间共有24042个差异表达基因,从中选取13个基因进行qRT-PCR检测,结果显示,9个基因表达上调,其余4个基因表达下调,与转录组测序结果一致。24042个差异表达基因中,共有10106个显著差异基因富集到129条代谢通路,其中富集程度排名前10的代谢途径为核糖体、次级代谢生物合成、RNA转运、内吞作用、剪接体、甘油磷脂代谢、内质网加工、吞噬、醚脂类代谢和植物-病原体相互作用,参与盐胁迫相关的硫代谢、脯氨酸积累、活性氧（ROS）代谢、与盐胁迫相关的钙信号通路和过氧化氢代谢等途径的差异基因上调。【结论】在盐胁迫下海马齿差异表达基因如小分子量热激蛋白基因、抗氧化酶相关基因及与离子交换相关基因发挥了重要调控作用。