• 中国中文核心期刊
  • 中国科学引文数据库(CSCD)核心库来源期刊
  • 中国科技论文统计源期刊(CJCR)
  • 第二届国家期刊奖提名奖

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

优先发表

优先发表栏目展示本刊经同行评议确定正式录用的文章,这些文章目前处在编校过程,尚未确定卷期及页码,但可以根据DOI进行引用。
显示方式:
基于连清样地数据的全国杉木人工林平均木树高-胸径模型
牛思圆, 刘鹏举, 雷相东, 任怡, 高影
, doi: 10.12403/j.1001-1498.20220322
摘要:
目的 基于我国森林资源连续清查(简称“连清”)样地数据,分省区研建全国杉木人工林平均木树高-胸径的最优基础模型,以期为全国各省区杉木人工林的树高预测提供基础模型。 方法 研究范围为杉木人工林分布的15个省份,数据来自第六次、第七次连清样地数据的树高调查表,总样本数为23 239个。选取18种基础生长方程作为候选模型,分别拟合各省区杉木平均木树高与胸径的关系,根据模型的决定系数(R2)、平均绝对误差(MAE)、平均相对误差(MRE)、均方根误差(RMSE)和平均预估误差(MPE),并结合模型残差分布图,确定各省区最优模型,同时采用5折法验证各省区最优模型的预测能力,最终决定各省区最优树高-胸径模型。 结果 15个省区的杉木最优树高-胸径模型并不相同,四川、云南、重庆、陕西、浙江、江西、湖南、广西的最优模型为模型18(Mitscherlich方程),江苏、安徽、河南和福建的最优模型为模型16(Hossfeld方程),广东、湖北、贵州的最优模型分别为模型10(双曲线方程)、模型11(Logistic方程)和模型13(Gompertz方程),R2分布在0.602~0.807之间,MAE分布在0.94~1.53 m之间,MRE分布在−2.93%~−4.72%之间,RMSE分布在1.23~2.00 m之间,MPE分布在0.50%~2.77%之间。模型拟合效果较好,满足精度要求,且参数具有生物学意义,可作为全国各省区杉木人工林平均木树高-胸径基础模型。 结论 本研究构建全国杉木人工林分布的15个省区的最优树高-胸径基础模型,能较好的模拟各省区的杉木平均木树高随胸径的变化规律,可以作为全国各省区基本的杉木人工林平均木树高-胸径模型,为各省区杉木人工林的树高预测提供依据。
银腺杨PagWOX11/12a基因对茎生长发育的影响
文爽爽, 王留强, 卢孟柱
, doi: 10.12403/j.1001-1498.20230264
摘要:
目的 分析PagWOX11/12a基因对杨树生长发育的影响,为进一步研究该基因在木本植物中的调控机制奠定基础。 方法 利用生物信息学方法及相关软件进行系统发育进化树构建、序列比对及生化特征分析。利用qRT-PCR分析其组织特异性表达模式。利用转基因植株35S::PagWOX11/12a-SRDX(DR),分析显性抑制PagWOX11/12a后杨树的表型。 结果 PagWOX11/12a基因可编码含255个氨基酸的蛋白质,该基因在84K的不同组织中均有表达。对DR转基因植株的表型分析表明,与非转基因植株相比,显性抑制该基因可导致茎中韧皮部细胞、髓心细胞及木质部纤维细胞长度变小,节间伸长受到抑制,株高明显降低。 结论 PagWOX11/12a基因通过影响节间伸长参与调控了杨树的高生长,该研究为进一步揭示PagWOX11/12a基因参与杨树生长发育的调控机制提供了参考。
无瓣海桑根响应盐胁迫的转录组分析
梁锐涛, 韩维栋, 杨少瑕, 陈蓓蓓
, doi: 10.12403/j.1001-1498.20230142
摘要:
摘要 : 目的 初步探究无瓣海桑的耐盐分子机制,筛选出无瓣海桑抗盐候选基因,为后续功能验证实验及林木抗盐性遗传育种奠定分子基础。 方法 以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个。 结论 活性氧清除、渗透调节、植物激素、蛋白激酶及转录因子相关基因参与调控无瓣海桑盐逆境适应过程。
九叶青花椒叶绿体基因组结构及系统进化分析
刘霞, 孙冲, 黄勤琴, 谢润泸, 刘浩文, 陈泽雄
, doi: 10.12403/j.1001-1498.20220277
摘要:
目的 揭示九叶青花椒叶绿体基因组的结构特征及其与其他物种的系统进化关系,为花椒属种质资源的鉴定、新品种的培育及品种间的遗传分析提供参考。 方法 使用改良CTAB法提取九叶青花椒叶片总DNA。利用BGISeq-500平台进行高通量测序,SPAdes v3.13.0软件组装叶绿体基因组,通过GeSeq注释九叶青花椒的全叶绿体基因组信息,对其叶绿体基因组序列的结构特征、重复序列、密码子偏好性及系统进化关系进行分析。 结果 九叶青花椒叶绿体基因组为典型的四分体组成结构,全长序列为158 558 bp,编码133个基因;测到19个串联重复序列,49个长片段重复,70个简单重复序列;九叶青花椒叶绿体基因组中的蛋白编码基因共有26 398个密码子(不包括终止密码子),在密码子第三位碱基上有较强的A/T碱基偏好性。 结论 本研究首次组装了九叶青花椒叶绿体基因组全序列,明确了花椒属为单系类群,九叶青花椒与野花椒(Z. simulans)关系密切,这将为花椒的遗传信息、花椒种质资源评价、分子育种、cpSSR分子标记的开发及遗传多样性等研究提供了重要参考。
马尾松次生维管组织的变化特征
杨超, 范付华, 徐刚
, doi: 10.12403/j.1001-1498.20220483
摘要:
目的 旨在揭示马尾松次生维管组织在发育过程中细胞结构、成分和功能的变化规律。 方法 通过组织化学分析和细胞离析等手段,研究了马尾松维管组织的细胞形态、结构、木质素和纤维素的变化。 结果 维管形成层分化形成新生韧皮部和木质部时会逐渐富集纤维素;韧皮部发育过程中筛胞发生形变、细胞壁木质化,韧皮部纤维素含量占比下降。木质部管胞由早材向晚材的发育过程中,纹孔直径降低,纹孔数量下降,次生壁木质化程度增加,胞腔面积下降,纤维素含量占比下降,弦向壁及径向壁增厚。 结论 次生维管组织内木质部与韧皮部的细胞结构和成分随发育进程会发生规律性变化。其中,木质部细胞的胞内及胞间输导能力下降,机械强度上升;次生韧皮部丧失输导能力,机械强度上升。
杉木连栽氨氧化古菌群落结构与硝态氮含量的关系
陈雯雯, 王淑真, 姜宇杰, 周垂帆
, doi: 10.12403/j.1001-1498.20220369
摘要:
目的 分析不同代数杉木林的土壤有效氮含量变化,探讨硝态氮与不同分类水平的氨氧化古菌(AOA)群落结构、多样性之间的联系,为人工林土壤氮素的有效性和杉木林土壤质量评估提供参考。 方法 在福建南平选取4个代数的杉木人工林,采用高通量测序技术对 PCR 扩增的 amoA 基因进行测定,运用 Mantel_r 相关分析、随机森林模型和偏最小二乘路径模型等分析方法研究不同代数的土壤有效氮含量与氨氧化古菌群落丰度、多样性之间的关系。 结果 随杉木连栽代数增加,硝态氮(NO3-N)含量显著下降,铵态氮(NH4+-N)和微生物量氮(MBN)含量变化不明显。土壤酶活性、氨氧化古菌群落丰度与多样性总体呈现降低趋势,土壤有效氮含量与氨氧化古菌群落以及酶活性之间联系密切,其中氨氧化古菌丰度与多样性指数是影响硝态氮含量的最主要因素。 结论 随着连栽代数增加,氨氧化古菌群落的丰富度与多样性一定程度上降低,除第四代土壤脲酶和氨单加氧酶酶活性略有上升,土壤氮循环相关酶活性基本呈现降低趋势,从而导致土壤硝态氮含量显著降低,限制了杉木连栽林土壤氮素的有效性。
秋茄低温胁迫转录组分析及脱落酸信号途径基因挖掘
郭晋敏, 杨升, 刘星, 王金旺, 王文卿, 陈秋夏
, doi: 10.12403/j.1001-1498.20220417
摘要:
目的 深入了解秋茄响应低温胁迫的分子机制以及培育秋茄的抗寒新品种。 方法 以耐寒红树植物‘龙港’秋茄1年生容器苗为实验材料,15 ℃处理12 h为对照组(CK),−5 ℃处理12 h为低温组(LT),采用 Illumina HiSeq测序平台进行转录组测序,挖掘脱落酸信号途径相关基因。 结果 转录组测序共鉴定到148个转录因子,分属于25个转录因子家族,其中ERF、NAC、WRKY、bHLH、MYB、bZIP、HB-other和MYB-related等家族所包含的基因数量较多,分别为17、14、12、12、10、9、6和6;差异组共筛选到1 330个差异表达基因(Differentially Expressed Genes, DEGs),其中698(52.48%)个上调表达,632(47.52%)个下调表达; KEGG通路富集分析发现,DEGs显著富集在植物激素信号转导、苯丙素生物合成、半乳糖代谢、光合作用-天线蛋白和α-亚麻酸代谢等通路中; 脱落酸信号途径相关基因 KoPYL1KoABF1KoABF2上调表达,KoPP2C1KoABF3下调表达,且这些基因表达情况与qRT-PCR验证结果一致。 结论 ERF、NAC、WRKY、bHLH、MYB、bZIP、HB-other和MYB-related等家族转录因子对秋茄响应低温胁迫起重要调控作用。植物激素信号转导、苯丙素生物合成、半乳糖代谢、光合作用-天线蛋白和α-亚麻酸代谢等是秋茄响应低温胁迫的重要KEGG通路。脱落酸信号途径中的KoPYL1KoPP2C1KoABF1KoABF2KoABF3基因可作为后期研究秋茄响应低温胁迫的重要候选基因。
白桦BpAMT基因家族鉴定及表达模式分析
杨海昕, 刘晓莹, 詹亚光, 范桂枝
, doi: 10.12403/j.1001-1498.20220242
摘要:
目的 鉴定白桦AMT基因家族成员,分析AMT基因的表达模式。 方法 本研究利用生物信息学方法鉴定了该家族成员并通过实时荧光定量技术分析其基因表达模式。 结果 从白桦基因组中鉴定了9个AMT家族成员,并将其分为AMT1AMT2两个亚家族,分别命名为BpAMT1.1~1.4BpAMT2.1~2.5;BpAMTs氨基酸残基数为384~522个,等电点为4.61~8.16,均定位于质膜与细胞器膜上;BpAMT基因不均匀的分布在5条染色体上,且成员间存在串联重复现象。BpAMT基因的表达具有组织部位特异性,呈现叶>根>茎趋势;同时发现,硝酸钾、氯化铵、茉莉酸甲酯、赤霉素、脱落酸、氯化镉、干旱、4 ℃和日变化等均可以影响BpAMT基因表达,且不同处理下家族成员的响应程度存在差异。 结论 从白桦中鉴定9个BpAMT基因,聚类为两个亚家族,在调节氮素吸收转运、响应激素信号及非生物胁迫中发挥不同的作用,该研究结果为深入解析BpAMT基因在白桦生长发育和抵抗逆境中的功能奠定基础。
PEG和NaCl胁迫下毛竹萌发种子的MicroRNAs表达谱分析
王晓静, 王涛, 杨凯, 李潞滨
, doi: 10.12403/j.1001-1498.20220211
摘要:
目的 鉴定和分析不同干旱和盐胁迫下毛竹种子在露白时期的microRNAs(miRNAs)及其表达模式。 方法 分别使用聚乙二醇(PEG6000)和氯化钠(NaCl)模拟干旱和盐胁迫,构建H2O、10% PEG、15% PEG、50 mmol·L−1 NaCl和100 mmol·L−1 NaCl处理下毛竹露白时期种子的small RNA文库,通过高通量测序和生物信息学分析揭示样本中的miRNA及其表达谱。 结果 通过small RNA测序共鉴定了246条已知miRNAs的成熟体序列,并预测得到262条novel miRNAs的成熟体序列;在毛竹种皮破裂阶段的种子中,丰度最高的已知miRNA为miR166,其次是miR159、miR6478、miR319等;根据miRNA靶基因预测结果,靶基因数量最多的已知miRNA属于MIR396家族,PH02Gene13935GAMYB)被预测受到MIR159、MIR319、MIR396家族共28个miRNAs的调控;在6个比较组中共鉴定了181个差异表达miRNA,与对照组相比,10% PEG、15% PEG、50 mmol·L−1 NaCl和100 mmol·L−1 NaCl胁迫下表达水平最高且显著差异表达的已知miRNA分别为phe-miR171e-5p、phe-miR3630-3p、phe-miR171e-5p和phe-miR159a.1;差异表达miRNA靶基因能够显著富集在不同的GO和KEGG途径中;对10个差异表达miRNAs进行qRT-PCR验证,荧光定量结果的整体趋势和测序数据一致。 结论 毛竹种皮破裂阶段种子中积累了大量的miR159、miR6478、miR319等已知miRNA,且在对照和胁迫处理组中均具有高表达水平,可能在毛竹种子萌发中具有保守调控作用;与对照组相比,phe-miR171e-5p,phe-miR3630-3p,phe-miR171e-5p和phe-miR159a.1在 10% PEG、15% PEG、50 mmol·L−1 NaCl、100 mmol·L−1 NaCl胁迫下分别显著差异表达,能够在毛竹种子露白阶段响应PEG或NaCl胁迫。
杜仲种子内生微生物群落组成及生态功能分析
张青青, 董醇波, 邵秋雨, 陆莹霞, 董旋, 梁宗琦, 韩燕峰
, doi: 10.12403/j.1001-1498.20220239
摘要:
目的 揭示杜仲种子内的微生物群落组成和多样性,丰富杜仲相关微生物资源。 方法 本研究采用传统培养法和高通量测序技术对杜仲种子内生微生物群落组成和多样性进行分析,同时基于FAPROTAX 和FUNGuild数据库对其进行功能注释。 结果 从杜仲种子中共分离获得真菌40株(3门11属),优势属为曲霉属(Aspergillus)(30%);细菌142株(1门6属),其中土壤芽孢杆菌属(Solibacillus)(47.18%)为主要优势属。高通量测序获得141个真菌OTUs,隶属于6门101属;442个细菌OTUs,隶属于24门313属。在属水平,种子优势真菌属为Apiotrichum(31.28%)和德巴利氏酵母属(Debaryomyces)(26.07%);优势细菌属主要为假单胞菌属(Pseudomonas)(16.66%)和乳杆菌属(Lactobacillus)(9.68%)。功能注释结果显示,真菌以降解有机物、促进宿主植物生长、促进宿主对土壤养分的吸收、改善土壤环境微生物群落结构以及致病菌群为主;细菌以化能异氧、有氧化能异养以及具有发酵和固氮作用的菌群为主。 结论 杜仲种子内携带多种有益功能菌和潜在致病菌,这些微生物对植物的生长发育具有重要影响。同时本研究还分离到多株有益功能菌,这将为后续人为调控种子微生物组以增强植物健康提供微生物资源。
梵净山典型植被类型的根系真菌群落结构及功能类群特征研究
刘敏, 李龙, 李海波, 王祖华
, doi: 10.12403/j.1001-1498.20220099
摘要:
目的 解析世界自然遗产地梵净山地区典型植被类型对真菌群落结构及功能的影响,为深入探讨可见的地上植物群落和未知的地下真菌群落结构间的关联机制提供数据支撑。 方法 以梵净山常绿阔叶林、常绿落叶阔叶混交林和亚高山矮林3种植被类型为研究对象,高通量测序技术探讨根系真菌群落结构在植被类型间的差异,FUNGuild对真菌进行功能预测,并分析真菌群落与土壤性质、细根性状和叶性状间的关系。 结果 子囊菌门和担子菌门为优势真菌门,undefined saprotroph、丛枝菌根真菌、外生菌根真菌、植物病原菌、内生真菌、寄生性真菌和杜鹃花类菌根真菌为主要功能类群。常绿阔叶林、常绿落叶阔叶混交林和亚高山矮林的显著差异功能类群分别是丛枝菌根真菌、外生菌根真菌和杜鹃花类菌根真菌。真菌群落组成在植被类型间存在显著差异,常绿阔叶林真菌多样性和Shanon指数显著高于常绿落叶阔叶混交林和亚高山矮林,而Bray-Curtis相异性指数在亚高山矮林显著低于常绿阔叶林和常绿落叶阔叶混交林。ABT分析表明,根碳含量是影响真菌丰富度的首要因子,叶干物质含量是影响真菌Shannon指数、Simpson指数和Pielou指数的首要因子。冗余分析揭示植被类型、土壤pH、根碳含量和叶碳含量显著影响真菌群落;方差分解分析表明植被类型对真菌群落的影响最大,其次是细根性状和土壤性质,叶性状影响最小。 结论 植被类型能显著影响真菌群落结构,真菌功能类群与植被类型间存在某种特定联系,植物性状对根系真菌群落的构建具有重要意义。
8种寄主幼树对美国白蛾幼虫生长发育及营养效应的影响
张书曼, 方静, 赵佳齐, 王楷, 任英杰, 孟昭军, 严善春
, doi: 10.12403/j.1001-1498.20220019
摘要:
目的 探究美国白蛾(Hyphantria cunea)对不同寄主幼树的适生性,科学保护林木、防控美国白蛾。 方法 以东北地区广泛分布及栽植的8种树(家榆、胡桃楸、旱柳、白桦、水曲柳、长白落叶松、山槐、紫椴)为实验对象,室内饲养并测定美国白蛾幼虫取食不同寄主幼树叶片后的生长发育和营养效应指标。 结果 取食家榆和胡桃楸的美国白蛾幼虫存活率较高,其次是水曲柳、旱柳和紫椴;取食长白落叶松的较低,其次是山槐。取食家榆和胡桃楸的美国白蛾幼虫体重、体长、头壳宽数值较高,其次是水曲柳、紫椴;取食山槐的较低,其次是白桦。取食胡桃楸、家榆、水曲柳的美国白蛾幼虫化蛹率、羽化率数值较高,发育历期较短,其次是紫椴;取食长白落叶松和山槐不能化蛹,取食白桦、山槐、长白落叶松的不能羽化且发育历期较长。美国白蛾幼虫对家榆、胡桃楸、水曲柳取食量较大,其次是紫椴、旱柳;取食家榆和胡桃楸的幼虫食物利用率和食物转换率均处于较高水平,其次是紫椴、水曲柳;幼虫对山槐和白桦的取食量较低,其食物利用率和食物转换率均处于较低水平;幼虫食物消耗率与其食物利用率和食物转换率结果相反。 结论 家榆和胡桃楸较适合美国白蛾幼虫生长发育,其次是水曲柳、紫椴;而长白落叶松、山槐不适于其生长;表明美国白蛾对多种寄主有不同的适生性和营养利用策略。
基于UHPLC–QTOF–MS技术的五倍子单宁成分分析
倪冰冰, 刘宏, 于立洋, 张国昀, 何彩云, 张建国
, doi: 10.12403/j.1001-1498.20210511
摘要:
目的 基于超高效液相色谱–四极杆飞行时间–质谱(UHPLC–QTOF–MS)技术分析五倍子中单宁成分,建立一种有效的五倍子单宁成分测定方法,并解析单宁的裂解规律。 方法 使用高温和超声波辅助方法提取五倍子中的酚类化合物,基于二级质谱推定化合物的裂解途径。 结果 共鉴定到20种化合物,包含14种没食子单宁、3种酚酸、2种没食子酸酯和1种黄烷醇。以高温121℃为提取条件,均未鉴定到分子质量高于1 500 Da的没食子单宁,而在超声波辅助提取的五倍子水提取液中,成功地检测到1–O–没食子酰基葡萄糖至14–O–没食子酰基葡萄糖,其中,超声功率1 500 W,提取温度65℃条件下,鉴定到更多的没食子单宁异构体。MS/MS图谱显示,五倍子单宁主要通过没食子酸、没食子酰基、糖苷、水、羧基的损失进行裂解。 结论 本研究证明,超声波辅助提取条件下,应用UHPLC–QTOF–MS技术能够全面高效地分析五倍子中的单宁成分,本研究为富含单宁的植物的化学成分研究提供参考依据。
西藏色季拉山西藏红杉径向生长对气温和降水波动的响应
于德水, 卢杰, 张萌, 张新生
, doi: 10.12403/j.1001-1498.20220144
摘要:
目的 基于色季拉山地区特有树种西藏红杉年轮样芯,研究西藏红杉径向生长对气温和降水量的波动响应。 方法 利用树木年代学的方法,结合TASP-Win及ARSTAN程序建立标准年表,采用Pearson及bootstrap的计算方式将年轮宽度指数与1961—2020年气温及降水分别进行相关分析。 结果 径向生长在气温波动前后存在较大的差异。在气温波动前西藏红杉与前一年8月,当年1月及6月平均气温呈正相关(P<0.05),与前一年9—10月及12月平均降水呈正相关,与前一年8月降水呈负相关,与1月至12月的相对湿度呈显著正相关(P<0.01)。在气温发生波动后,与当年3—12月平均气温呈负相关,与前一年8月及当年6月降水为显著正相关。 结论 藏东南地区的气候条件极为特殊,气候的波动导致水热条件发生变化是限制西藏红杉径向生长的主要原因,在非生长季,气温短暂的变化对树木径向生长同样有着显著的影响。
六盘山华北落叶松林的结构随林龄变化及其水文影响
徐小勤, 于澎涛, 王彦辉, 徐丽宏, 于艺鹏, 刘泽彬
, doi: 10.12403/j.1001-1498.20220218
摘要:
目的 探讨林分结构随林龄的变化及其所带来的水文影响,为六盘山半干旱区森林植被建设和经营管理提供理论依据。 方法 在六盘山半干旱区叠叠沟小流域,以主要人工林类型—华北落叶松人工林为研究对象,基于1993—2002年的日气象数据,应用林分结构耦合模型和BROOK90生态水文模型,通过情景模拟,得到植被结构(平均树高、胸径、郁闭度、LAI)随林龄的变化,及其对水分平衡分量(截留量、蒸腾量、土壤蒸发量、产流量)的影响。 结果 华北落叶松林随林龄的增加可分为快速生长期(林龄5~15 a)、缓慢生长期(林龄15~40 a)和稳定期(林龄大于40 a)。在快速生长期,林分平均树高、平均胸径、郁闭度和冠层LAI均随林龄的增加而快速增大,年均增长率分别为0.48 m·a−1、0.54 cm·a−1、0.03、0.19;在缓慢生长期,林分平均树高、平均胸径和郁闭度呈缓慢增大趋势,而冠层LAI则呈先增大后减小的趋势;在稳定期,林分结构基本稳定。在快速生长期,年均截留量和蒸腾量随林龄增加而快速增大,年均增长率分别为1.91、24.13 mm·a−1;年均土壤蒸发量和产流量随林龄增加快速减小,其年均减小速率为10.58、14.88 mm·a−1。在缓慢生长期,年均蒸腾量、土壤蒸发量和产流量随林龄的变化减慢,并在林龄30 a时趋于稳定,其年均变化率分别为0.62、−0.75、−0.76 mm·a−1 结论 六盘山半干旱区华北落叶松人工林在林龄小于30 a,特别是林龄小于15 a时,林分结构不断变化,显著影响人工林的生态水文功能。在林龄30 a后,林分结构稳定,水分平衡分量也趋于稳定,在林水综合管理中可不必再考虑林龄的影响。
金丝楸幼苗响应盐碱胁迫的生理和转录组分析
郜新强, 王小艳, 焦伟, 王静, 郑丽月, 王丹宁, 王兴云, 侯立江, 尚增振, 刘彦珍
摘要:
目的 研究不同盐碱胁迫对金丝楸幼苗生长、光合和生理指标的影响,并结合转录组测序分析,探究楸树耐盐碱的生理机制和分子机制。 方法 采用盆栽法对金丝楸幼苗进行不同盐碱胁迫处理,分析其生物量、光合及生理指标对不同盐碱响应的差异,采用Illumina高通量测序技术进行转录组测序,通过生物信息学分析盐碱胁迫对转录水平的影响。 结果 不同盐碱胁迫下,金丝楸幼苗叶片受伤害程度为Na2CO3>混合盐碱>NaCl;新增株高和地径、地上部和根的干质量和鲜质量、生物量、根冠比均受到明显抑制,并随盐碱浓度增加而抑制加强,但生长胁迫指数均随浓度的增加而降低;丙二醛(MDA)含量和相对电导率都随胁迫浓度增加而不同程度上升,超氧化物歧化酶(SOD)活性、可溶性糖含量、脯氨酸(Pro)含量、叶绿素总量和光合速率都呈现先上升后下降的趋势。转录组测序共产生约60.4 Gb原始数据,组装得到55 793个Unigenes,其中29 534(52.93%)个Unigenes获得了注释;通过差异表达基因(DEGs)分析,3个比较组(CK vs NaCl,CK vs Na2CO3 和 CK vs 混合盐碱)分别筛选出1 779、2 835和4 059 个DEGs;DEGs GO富集分析表明,膜的整体成分、膜的内在成分、催化活性、类异戊二烯代谢和合成过程、氧化还原酶活性等条目被显著富集;DEGs KEGG分析表明,苯丙素生物合成、淀粉和蔗糖代谢、植物激素信号转导、萜类主干生物合成和精氨酸代谢等通路被显著富集;此外,在DEGs中鉴定的bHLH、ERF、MYB-related、NAC、C2H2、WRKY、MYB和bZIP转录因子家族成员最多。 结论 金丝楸主要通过积累可溶性糖和Pro,提高SOD酶活和光合作用来抵御盐碱胁迫,但都呈现“低促高抑”的现象,说明其具有一定阈值。金丝楸通过调节膜成分、催化活性、类异戊二烯代谢和生物合成过程、苯丙素生物合成、淀粉和蔗糖代谢、植物激素信号转导等生物过程和代谢途径,并结合有关转录因子共同响应盐碱胁迫。本研究为深入研究楸树耐盐碱生理机制和分子机制提供科学的理论依据。