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

张　蕾，阮羽萱，潘　娇，等．甘蓝型油菜响应干旱胁迫的转录组分析［Ｊ］．江苏农业科学，２０２３，５１（５）：５１－５５．
ｄｏｉ：１０．１５８８９／ｊ．ｉｓｓｎ．１００２－１３０２．２０２３．０５．００７
甘蓝型油菜响应干旱胁迫的转录组分析
张　蕾，阮羽萱，潘　娇，陈　敏，刘博宇，阮　颖，黄　勇
（湖南农业大学生物科学技术学院／湖南省作物表观遗传与发育重点实验室，湖南长沙４１０１２８）
摘要：甘蓝型油菜（ＢｒａｓｓｉｃａｎａｐｕｓＬ．）是我国主要油料作物之一，为国家粮油安全提供了重要保障。

干旱严重限制了油菜种植规模扩大、产量与品质提升。

为挖掘甘蓝型油菜抗旱的关键基因，阐明油菜耐旱机制，本研究以甘蓝型油菜湘油１
５为对象，通过实验室模拟干旱处理后的叶片为材料进行转录组测序（ＲＮＡ－Ｓｅｑ）、生物信息分析并ｑＲＴ－ＰＣＲ验证。

发现转录组样本平均Ｃｌｅａｎｂａｓｅａｓ约为９．６Ｇ，且各样品碱基百分比（Ｑ３０）大于９３．２１％，比对到参考基因组上的Ｒｅａｄｓ在８９％以上。

干旱处理４ｈ后，湘油１５出现了４４０个差异表达基因（ＤＥＧｓ），其中１４８个基因上调，２９２个基因下调。

随机选取１０个基因（上调与下调各５个）进行ｑＲＴ－ＰＣＲ验证，其中９个基因的ｑＲＴ－ＰＣＲ结果与ＲＮＡ－Ｓｅｑ结果相符，仅１个基因表达趋势相反，相似度可达９０％，证实转录组测序结果的可靠性。

结果表明甘蓝型油菜湘油１５主要通过光合作用、碳代谢和渗透调节等途径调节抗旱能力，为进一步阐明油菜耐旱机制奠定了基础。

关键词：甘蓝型油菜；干旱胁迫；转录组分析；差异表达基因
中图分类号：Ｓ６３４．３０１ 文献标志码：Ａ 文章编号：１００２－１３０２（２０２３）０５－００５１－０５
收稿日期：２０２２－０６－２１
基金项目：国家自然科学基金面上项目（编号：３１９７１８３４）；湖南省教育厅优秀青年项目（编号：
１９Ｂ２６６）。

作者简介：张　蕾（１９９８—），女，湖南衡阳人，硕士研究生，从事植物分子生物学研究。

Ｅ－ｍａｉｌ：ｚｌ８３９８＠１６３．ｃｏｍ。

通信作者：黄　勇，博士，教授，从事分子生物学研究。

Ｅ－ｍａｉｌ：ｙｏｎｇｈｕａｎｇ＠ｈｕｎａｕ．ｅｄｕ．ｃｎ。

干旱是全球发生频率最高、影响范围最广、对
国民经济造成巨大损失的自然灾害［１－２］。

油菜是我国种植面积最广、产量最高的油料作物之一［
３］
，其中甘蓝型油菜是广泛种植于长江中下游区域的油菜品种，但其整个生长期对水的依赖性较强，耐旱性较差。

干旱极大地影响了油菜的发育、产量和品
质，严重威胁人类的粮油安全［４］。

当遭受干旱胁迫时，植物在一定程度上能够通过自身调节适应环境的变化。

气孔是植物散失水分的主要途径，当植物处于干旱胁迫下，将调控叶片气孔闭合以降低蒸腾作用，从而减少植物体内水
分流失，进而影响光合作用与碳代谢［５］。

此外，植
物细胞中的一些离子和渗透调节物质将累积以维持细胞内离子稳态和渗透调节平衡，从而使自身保存水分的能力保持在较高的水平
［６］。

当植物处于
逆境胁迫下，细胞代谢发生紊乱，大量积累的活性氧破坏了生物膜的结构，影响其功能的正常发挥，
严重时危害细胞生长，导致植株死亡［７］。

为减轻活
性氧积累带来的影响，作物进化出复杂的酶和抗氧
化剂系统，如低分子抗氧化剂和活性氧清除酶［
８］。

转录水平上的调控是生物基因表达的主要调控方式之一，是对生物进行研究的关键组成部
分［９］。

高通量测序技术的出现，使研究者能够更有
效地研究不同处理条件或不同发育过程中基因的
差异表达谱，筛选新的功能基因［１０－１２］。

在本研究
中，用２０％聚乙二醇（ＰＥＧ－６０００）模拟干旱胁迫处理甘蓝型油菜湘油１５，借助高通量测序技术进行分析，筛选干旱处理诱导表达基因信息，初步挖掘耐旱基因和代谢通路，构建抗旱性相关的调控网络。

为甘蓝型油菜耐旱机制提供理论依据，为培育耐旱品种提供参考。

１　材料与方法
１．１　植物材料和干旱处理
植物材料为甘蓝型油菜湘油１５，先将其种子进行表面消毒，后播种于含有霍格兰营养液的蛭石中，温室条件下正常生长４５ｄ（２５℃，１６ｈ／８ｈ光照／黑暗，光照强度为１００００ｌｘ，湿度为４０％～６０％），长出３～４片真叶后小心洗去蛭石，置于霍格兰营养液中并充气，光照培养箱中预培养３ｄ，根据营养液体积，加入２０％聚乙二醇（ＰＥＧ－６０００）模拟干旱胁迫，处理时间为４ｈ。

选取生长状态基本一致的５株植株进行混合取样，处理前的样品编号为
Copyright ©博看网. All Rights Reserved.
Ｃ１、Ｃ２、Ｃ３，处理后的样品编号为Ｔ１、Ｔ２、Ｔ３。

样本于液氮中速冻，保存于－８０℃中，用于后续试验。

１．２　试验时间与地点
２０２１年７月２１日于湖南省作物表观遗传与发育重点实验室进行。

１．３　转录组测序
样品由百迈克生物科技（长沙）有限公司进行ＲＮＡ提取、质量分析、ｃＤＮＡ文库构建和转录组测序分析［１３］。

１．４　差异表达基因的筛选和分析
首先对样品中ＭａｐｐｅｄＲｅａｄｓ数目与转录本长度归一化，以ＦＰＫＭ方法作为检测转录本或基因表达水平的指标［１４］。

时间和空间特异性是基因表达的２个特性，当处于不同的条件下，表达水平具有显著差异的基因或转录本称为差异表达基因（ＤＥＧｓ）。

根据不同样品间的相对表达水平，可将ＤＥＧｓ划分为下调基因和上调基因。

利用Ｒ语言ＤＥＳｅｑ２程序包对过滤后的数据进行筛选来获得ＤＥＧｓ，筛选标准为：差异倍数（ｆｏｌｄｃｈａｎｇｅ，简称ＦＣ）大于２，即｜ｌｏｇ
２
ＦＣ｜≥１；错误发现率（ｆａｌｓｅｄｉｓｃｏｖｅｒｙｒａｔｅ，简称ＦＤＲ）小于０．０１［１５］，ＦＤＲ是通过对差异显著性Ｐ值（ｐ－ｖａｌｕｅ）进行校正得到的，即ｐ－ｖａｌｕｅ＜０．０１。

使用Ｂｌａｓｔ２ＧＯ软件对差异基因进行基因功能注释。

１．５　通过ｑＲＴ－ＰＣＲ验证差异表达
随机选取１０个ＤＥＧｓ（上、下调基因各５个）进行实时荧光定量分析（ｑＲＴ－ＰＣＲ）验证转录组测序中差异基因表达的准确性。

通过Ｔｒｉｚｏｌ法提取转录组测序备份样品的总ＲＮＡ，每样品取１ｎｇＲＮＡ通过擎科生物ｃＤＮＡ合成试剂盒反转录，最终ｃＤＮＡ体积为２０μＬ。

ｑＰＣＲ引物设计见表１。

ｑＲＴ－ＰＣＲ使用ＡｐｐｌｉｅｄＢｉｏｓｙｓｔｅｍｓ７３００实时ＰＣＲ系统，设置３次技术重复，使用２－ΔΔＣＴ法计算相对定量数据。

表１　ｑＰＣＲ的引物及内参引物
基因正向引物（５′→３′） 反向引物（５′→３′） 产物长度（ｂｐ）
ＡｃｔｉｎＧＧＴＴＧＧＧＡＴＧＧＡＣＣＡＧＡＡＧＧＴＣＡＧＧＡＧＣＡＡＴＡＣＧＧＡＧＣ１６５ＢｎａＣ０６ｇ３８９１０ＤＧＡＴＧＡＴＴＡＣＴＧＡＴＧＧＧＴＴＡＧＧＧＡＴＧＡＧＡＧＧＧＣＡＴＧＧＴＧＧＡ１１４ＢｎａＡ０８ｇ１３８６０ＤＴＴＡＣＡＧＣＧＧＧＡＴＧＧＧＧＴＧＡＡＧＧＧＴＴＴＣＡＴＧＧＧＡＣＡＡＧ９５ＢｎａＡｎｎｇ１３０４０ＤＴＧＡＡＴＧＴＴＧＴＣＡＧＴＧＡＧＧＡＧＧＴＴＴＧＧＡＧＧＡＡＡＡＡＡＧＣＧＧ１５７ＬＯＣ１０６４２９２１１ＡＣＧＣＴＣＴＴＴＡＴＡＴＣＡＴＣＡＣＣＡＡＧＧＧＡＣＧＡＡＡＣＡＧＣＴＧＴＴＣＣＣＴＣ８７ＢｎａＡ０７ｇ１１４５０ＤＣＴＣＡＡＧＡＡＡＣＧＡＣＣＡＣＡＣＣＧＴＧＣＴＴＣＧＣＣＡＣＡＡＡＣＧＣ８２ＢｎａＡｎｎｇ３０６８０ＤＣＧＧＡＧＡＧＡＴＧＧＡＡＡＡＡＧＣＧＧＧＡＡＧＡＧＧＡＴＧＣＧＡＧＡＣＴＧＧ１６７ＢｎａＣ０２ｇ４６７００ＤＣＣＴＴＴＧＣＧＧＴＴＴＴＣＧＧＧＣＴＴＧＧＴＣＡＣＴＣＴＧＣＴＣＴＴＧ１１８ＢｎａＣ０５ｇ００８４０ＤＡＧＧＣＧＴＣＡＡＡＴＧＴＧＧＡＧＡＧＧＴＡＧＴＡＧＣＡＧＡＧＧＡＴＧＴＴＧＴＧＧＴＡ１４８ＬＯＣ１０６４２７８８４ＣＴＣＴＣＴＣＴＡＡＡＣＡＣＴＧＴＣＧＣＴＡＴＣＣＴＴＧＴＣＧＴＧＣＴＣＴＴＣＣＴＣＣ１８４ＢｎａＣ０８ｇ１２６６０ＤＧＣＴＣＴＡＧＧＴＡＴＣＣＴＣＧＧＴＧＣＴＣＴＡＴＧＡＧＴＧＴＣＴＧＴＣＧＴＧＴＣＧ１５７ 注：基因名称来源于ＮＣＢＩ和ｇｒａｍｅｎｅ数据库比对的结果。

２　结果与分析
２．１　转录组测序质量分析
转录组测序分析对照组和干旱处理试验组各３个样本。

高通量测序后，为保证数据分析的质量，对原始数据进行过滤得到ｃｌｅａｎｒｅａｄｓ。

整体上分析，干旱处理组和对照组的３个重复平均ｒｅａｄｓ数分别约为３４０９４３６９、３０１９０７５７，样本平均碱基数约为９．６Ｇ，且各样品Ｑ３０碱基百分比在９３．２１％及以上，ＧＣ含量为４７．５７％～４８．８％。

使用ＨＩＳＡＴ２比对系统分析得到，各样品ｒｅａｄｓ与参考基因组的比对效率在８９．１８％～９１．６５％（表２）。

整体测序数据质量良好，可用于后续分析。

２．２　差异表达基因分析
利用ＤＥＳｅｑ２对甘蓝型油菜各样本的转录组数据进行分析，以Ｐ值＜０．０１且差异倍数｜ｌｏｇ
２
ＦＣ｜≥１作为显著差异表达基因的筛选条件。

从图１可知，对照组和试验组（ＣｖｓＴ）中鉴定到ＤＥＧｓ共有４４０个，其中１４８个ＤＥＧｓ表达上调，２９２个ＤＥＧｓ表达下调。

２．３　差异表达基因ＧＯ分析
以Ｐ≤０．０１为标准将差异表达基因进行ＧＯ注释。

根据差异基因对应的生物过程、分子功能和细胞组分的不同，分别进行ＧＯ富集分析。

数据结果统计发现，生物过程主要集中在韧皮部蔗糖的装载、
Copyright©博看网. All Rights Reserved.
表２　样本测序数据质量控制数据统计
样品干净序列读数干净序列碱基数
（Ｇ）
Ｇ和Ｃ含量
（％）
Ｑ３０碱基比例
（％）
比对率
（％）Ｃ１３５６７９１９８１０．６７Ｇ４８．８９５．０８３２６３０８７７（９１．４６）Ｃ２２３２４０４０５６．９３Ｇ４７．５７９３．２１２０９４２０００（９０．１１）Ｃ３３１６５２６６９９．４７Ｇ４８．５８９５．０５２９００９２８６（９１．６５）Ｔ１４０４８２７２５１２．１０Ｇ４８．３２９４．８９２９７３５１７８（９０．８８）Ｔ２３７６９１６４１１１．２７Ｇ４８．２８９４．６６３２８９３１２５（９１．２５）Ｔ３２４１０８７４１７．２１Ｇ４７．７９９３．７０２１５０１１４１（８９．１８）
总计
１９２８５５３７９
５７．６４Ｇ
≥４
７．５７≥９
３．２１≥８
９．１
８应对镉离子、磷酸戊糖分流、糖酵解、叶绿体移动等（图２－
Ａ）；分子功能主要集中在果糖二磷酸醛缩酶活性、磷酸核酮糖激酶活性、吡哆醛磷酸的结合等（图２－Ｂ）；细胞组分主要集中在胞质核糖体、质膜、质外体、叶绿体被膜等（图２－Ｃ）。

结果表明，甘蓝型油菜在干旱胁迫下的差异表达基因主要与植物生长发育如光合作用和糖代谢过程有关。

甘蓝型油菜可以通过体内细胞全方位的协调作用来减少外部恶劣环境对植物体所造成的伤害，同时调节内部代谢活动以保证植物的基本生长。

２．４　差异表达基因ＫＥＧＧ分析
为分析在干旱胁迫下甘蓝型油菜发育过程中体内代谢途径的变化，对差异基因进行ＫＥＧＧ分类统计（Ｐ＜０．０１）。

这些差异基因主要富集通路包括碳代谢、乙醛酸和二羧酸代谢、氨基酸生物合成和代谢、磷酸戊糖途径（图３）。

代谢途径分析表明，干旱胁迫对甘蓝型油菜的碳代谢、氨基酸代谢有影响。

２．５　差异表达基因的ｑＲＴ－ＰＣＲ验证
随机选取上调与下调各５个差异表达基因验证转录组测序结果的可靠性，通过ｑＲＴ－ＰＣＲ验证干旱处理前后的表达水平，将所选基因的相对表达与ＲＮＡ－ｓｅｑ分析的相对表达进行比较，结果表明，
ｑＲＴ－ＰＣＲ结果的相对趋势与ＲＮＡ－ｓｅｑ数据基本一致，证明本研究ＲＮＡ－ｓｅｑ结果的可靠性（图４）。

３　讨论与结论
光合作用直接影响作物生长发育及产量，干旱环境中光合作用强度是衡量作物抗旱能力的重要
依据［
１６－１７］。

笔者发现，差异表达基因ＧＯ和ＫＥＧＧ注释到的结果与光合作用、碳代谢、渗透调节有关，研究结果与前人在不同物种的发现较为一致。

研
究发现，棉花［
１８－１９］和苏丹草［２０］
干旱转录组的ＫＥＧＧ代谢通路显著富集于光合作用中。

对梭
梭［２１］和甘蓝型油菜［２２］
的干旱转录组分析表明，其
ＫＥＧＧ富集于碳代谢等过程。

碳水化合物可以为逆境下的植物提供基本能量以维持正常生命活动，同时还可以作为体内信号分子调节一系列代谢反应
以适应不利环境［２３］。

植物受到水分胁迫时，渗透调
节系统会改变植物的物质代谢过程，一些大分子物质如蛋白质、淀粉等会分解成小分子有机物如可溶
性糖和氨基酸等［
２４］
，它们能作为渗透调节因子保持植物体内渗透压的平衡，保护逆境环境中的植物细胞。

对干旱胁迫棉花的转录组分析显示，差异表达
基因显著富集在渗透调节系统上［１３］。

在干旱胁迫期间，为通过调节细胞过程在干旱
条件下生存，植物中的基因表达发生了显著变化［
２５］。

通过差异基因表达方法阐明植物的耐旱性，为确定可能的耐旱机制提供了有价值的信息。

对低温胁迫下云雾贡茶进行转录组分析，筛选并鉴定了差异
基因ＣｓＰＰＯ具有提高植物耐寒性的功能［２６］。

对盐
碱胁迫下的紫花苜蓿转录组分析筛选到差异基因ＭｓＮＡＣ４７，对其功能研究发现ＭｓＮＡＣ４７通过调控钙
离子的转运使植物抵御盐碱胁迫［２７］。

在本研究中，
ＰＥＧ－６０００模拟干旱处理前后湘油１５共有４４０个ＤＥＧｓ，其中１４８个表达上调，２９２个表达下调，并随
Copyright ©博看网. All Rights Reserved.
机对其中１０个ＤＥＧｓ进行ｑＲＴ－ＰＣＲ进行验证。

虽然结果与ＲＮＡ－ｓｅｑ检测的基因表达量上调或下调的倍数存在差异，但这两者的趋势基本一致，表明转录组测序富集到ＤＥＧｓ数据是准确的，推测此ＤＥＧｓ对于甘蓝型油菜抗旱有调控作用，可以进行基因克隆与功能验证等深入研究。

项目有效地筛选了甘蓝型油菜湘油１５干旱胁迫响应基因，这些基因主要调控光合作用、碳代谢
Copyright©博看网. All Rights Reserved.
和渗透调节等过程，研究结果为阐明甘蓝型油菜的耐旱机制提供了参考。

参考文献：
［１］韩兰英．气候变暖背景下中国农业干旱灾害致灾因子、风险性特征及其影响机制研究［Ｄ］．兰州：兰州大学，２０１６．
［２］李茂松，李　森，李育慧．中国近５０年旱灾灾情分析［Ｊ］．中国农业气象，２００３，２４（１）：７－１０．
［３］张治礼，郑学勤，吕应堂．内源细胞分裂素调控油菜叶片衰老进程的研究［Ｊ］．作物学报，２００５，３１（１）：１－６，１３７．
［４］黄　倩，赵永国，黄祥伟，等．甘蓝型油菜蕾薹期抗旱相关性状的ＱＴＬ分析［Ｊ］．干旱地区农业研究，２０１７，３５（６）：８８－９４．
［５］ＷｉｌｋｉｎｓｏｎＳ，ＤａｖｉｅｓＷＪ．Ｄｒｏｕｇｈｔ，ｏｚｏｎｅ，ＡＢＡａｎｄｅｔｈｙｌｅｎｅ：ｎｅｗｉｎｓｉｇｈｔｓｆｒｏｍｃｅｌｌｔｏｐｌａｎｔｔｏｃｏｍｍｕｎｉｔｙ［Ｊ］．Ｐｌａｎｔ，Ｃｅｌｌ＆Ｅｎｖｉｒｏｎｍｅｎｔ，２０１０，３３（４）：５１０－５２５．
［６］ＡｎｊｕｍＳ，ＸｉｅＸＹ，ＷａｎｇＬＣ，ｅｔａｌ．Ｍｏｒｐｈｏｌｏｇｉｃａｌ，ｐｈｙｓｉｏｌｏｇｉｃａｌａｎｄｂｉｏｃｈｅｍｉｃａｌｒｅｓｐｏｎｓｅｓｏｆｐｌａｎｔｓｔｏｄｒｏｕｇｈｔｓｔｒｅｓｓ［Ｊ］．ＡｆｒｉｃａｎＪｏｕｒｎａｌｏｆＡｇｒｉｃｕｌｔｕｒａｌＲｅｓｅａｒｃｈ，２０１１，６（９）：２０２６－２０３２．
［７］张　妍．外源ＥＢＲ和ＮＯ增强华北驼绒藜耐盐及耐旱性的生理机制研究［Ｄ］．呼和浩特：内蒙古大学，２０２１．
［８］谭筱玉．利用转录组解析油菜苗期响应干旱和盐胁迫的分子机制［Ｄ］．武汉：华中农业大学，２０２０．
［９］胡玉锋，金峰学，程云清，等．干旱胁迫下玉米苗期的转录组分析［Ｊ］．东北农业科学，２０２１，４７（６）：４８－５２．
［１０］ＭａｌｉｎｉｃｈＥＡ，ＷａｎｇＫ，ＭｕｋｈｅｒｊｅｅＰＫ，ｅｔａｌ．ＤｉｆｆｅｒｅｎｔｉａｌｅｘｐｒｅｓｓｉｏｎａｎａｌｙｓｉｓｏｆＴｒｉｃｈｏｄｅｒｍａｖｉｒｅｎｓＲＮＡｒｅｖｅａｌｓａｄｙｎａｍｉｃｔｒａｎｓｃｒｉｐｔｏｍｅｄｕｒｉｎｇｃｏｌｏｎｉｚａｔｉｏｎｏｆＺｅａｍａｙｓｒｏｏｔｓ［Ｊ］．ＢＭＣＧｅｎｏｍｉｃｓ，２０１９，２０（１）：１－１９．
［１１］ＫｕｍａｒＭ，ＫｉｍＩ，ＫｉｍＹＫ，ｅｔａｌ．ＳｔｒｉｇｏｌａｃｔｏｎｅｓｉｇｎａｌｉｎｇｇｅｎｅｓｓｈｏｗｉｎｇｄｉｆｆｅｒｅｎｔｉａｌｅｘｐｒｅｓｓｉｏｎｐａｔｔｅｒｎｓｉｎＡｒａｂｉｄｏｐｓｉｓｍａｘｍｕｔａｎｔｓ［Ｊ］．Ｐｌａｎｔｓ，２０１９，８（９）：３５２．
［１２］ＺｈｕＪＹ，ｅｔａｌ．Ｃｏｍｐｒｅｈｅｎｓｉｖｅｃｏ－ｅｘｐｒｅｓｓｉｏｎａｎａｌｙｓｉｓｐｒｏｖｉｄｅｓｎｏｖｅｌｉｎｓｉｇｈｔｓｉｎｔｏｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｏｆｆｌａｖｏｎｏｉｄｓｉｎｆｒｅｓｈｌｅａｖｅｓｏｆｔｈｅｔｅａｐｌａｎｔ（Ｃａｍｅｌｌｉａｓｉｎｅｎｓｉｓ）［Ｊ］．ＰｌａｎｔＳｃｉｅｎｃｅ，２０２０，２９０：１１０３０６．
［１３］包秋娟．干旱胁迫下棉花转录组分析［Ｄ］．乌鲁木齐：新疆大
学，２０１８．
［１４］ＭｏｒｔａｚａｖｉＡ，ＷｉｌｌｉａｍｓＢＡ，ＭｃＣｕｅＫ，ｅｔａｌ．ＭａｐｐｉｎｇａｎｄｑｕａｎｔｉｆｙｉｎｇｍａｍｍａｌｉａｎｔｒａｎｓｃｒｉｐｔｏｍｅｓｂｙＲＮＡ－Ｓｅｑ［Ｊ］．ＮａｔｕｒｅＭｅｔｈｏｄｓ，２００８，５（７）：６２１－６２８．
［１５］陈智勇．甘蓝型油菜叶片生长对干旱胁迫的响应分析［Ｄ］．福州：福建农林大学，２０２０．
［１６］ＧａｌｍéｓＪ，ＭｅｄｒａｎｏＨ，ＦｌｅｘａｓＪ．ＰｈｏｔｏｓｙｎｔｈｅｔｉｃｌｉｍｉｔａｔｉｏｎｓｉｎｒｅｓｐｏｎｓｅｔｏｗａｔｅｒｓｔｒｅｓｓａｎｄｒｅｃｏｖｅｒｙｉｎＭｅｄｉｔｅｒｒａｎｅａｎｐｌａｎｔｓｗｉｔｈｄｉｆｆｅｒｅｎｔｇｒｏｗｔｈｆｏｒｍｓ［Ｊ］．ＮｅｗＰｈｙｔｏｌｏｇｉｓｔ，２００７，１７５（１）：８１－９３．［１７］吴金芝，黄　明，王志敏，等．干旱对冬小麦旗叶光合参数、产量和水分利用效率的影响［Ｊ］．江苏农业学报，２０２１，３７（５）：１１０８－１１１８．［１８］ＳｉｎｇｈＲ，ＰａｎｄｅｙＮ，ＮａｓｋａｒＪ，ｅｔａｌ．ＰｈｙｓｉｏｌｏｇｉｃａｌｐｅｒｆｏｒｍａｎｃｅａｎｄｄｉｆｆｅｒｅｎｔｉａｌｅｘｐｒｅｓｓｉｏｎｐｒｏｆｉｌｉｎｇｏｆｇｅｎｅｓａｓｓｏｃｉａｔｅｄｗｉｔｈｄｒｏｕｇｈｔｔｏｌｅｒａｎｃｅｉｎｃｏｎｔｒａｓｔｉｎｇｖａｒｉｅｔｉｅｓｏｆｔｗｏＧｏｓｓｙｐｉｕｍｓｐｅｃｉｅｓ［Ｊ］．Ｐｒｏｔｏｐｌａｓｍａ，２０１５，２５２（２）：４２３－４３８．
［１９］ＣｈｅｎＹ，ＬｉｕＺＨ，ＦｅｎｇＬ，ｅｔａｌ．Ｇｅｎｏｍｅ－ｗｉｄｅｆｕｎｃｔｉｏｎａｌａｎａｌｙｓｉｓｏｆｃｏｔｔｏｎ（Ｇｏｓｓｙｐｉｕｍｈｉｒｓｕｔｕｍ）ｉｎｒｅｓｐｏｎｓｅｔｏｄｒｏｕｇｈｔ［Ｊ］．ＰＬｏＳＯｎｅ，２０１３，８（１１）：ｅ８０８７９．
［２０］ＺｈｕＹＱ，ＷａｎｇＸ，ＨｕａｎｇＬＫ，ｅｔａｌ．Ｔｒａｎｓｃｒｉｐｔｏｍｉｃｉｄｅｎｔｉｆｉｃａｔｉｏｎｏｆｄｒｏｕｇｈｔ－ｒｅｌａｔｅｄｇｅｎｅｓａｎｄＳＳＲｍａｒｋｅｒｓｉｎＳｕｄａｎｇｒａｓｓｂａｓｅｄｏｎＲＮＡ－ｓｅｑ［Ｊ］．ＦｒｏｎｔｉｅｒｓｉｎＰｌａｎｔＳｃｉｅｎｃｅ，２０１７，８：６８７．
［２１］ＬｏｎｇＹ，ＺｈａｎｇＪＷ，ＴｉａｎＸＪ，ｅｔａｌ．ＤｅｎｏｖｏａｓｓｅｍｂｌｙｏｆｔｈｅｄｅｓｅｒｔｔｒｅｅＨａｌｏｘｙｌｏｎａｍｍｏｄｅｎｄｒｏｎ（Ｃ．Ａ．Ｍｅｙ．）ｂａｓｅｄｏｎＲＮＡ－Ｓｅｑｄａｔａｐｒｏｖｉｄｅｓｉｎｓｉｇｈｔｉｎｔｏｄｒｏｕｇｈｔｒｅｓｐｏｎｓｅ，ｇｅｎｅｄｉｓｃｏｖｅｒｙａｎｄｍａｒｋｅｒｉｄｅｎｔｉｆｉｃａｔｉｏｎ［Ｊ］．ＢＭＣＧｅｎｏｍｉｃｓ，２０１４，１５（１）：１１１１．［２２］王力敏．甘蓝型油菜应答干旱胁迫的组学研究［Ｄ］．武汉：华中农业大学，２０１７．
［２３］李　亮．玉米雄穗响应干旱胁迫的主要代谢途径解析及耐旱ＳＮＡＣ基因筛选鉴别［Ｄ］．乌鲁木齐：新疆农业大学，２０１５．［２４］张　华．ＡＭ真菌对小叶锦鸡儿和蒙古扁桃生长、抗旱性及总黄酮含量的影响［Ｄ］．呼和浩特：内蒙古大学，２０１３．
［２５］ＳｈａｎｋｅｒＡＫ，ＭａｈｅｓｗａｒｉＭ，ＹａｄａｖＳＫ，ｅｔａｌ．Ｄｒｏｕｇｈｔｓｔｒｅｓｓｒｅｓｐｏｎｓｅｓｉｎｃｒｏｐｓ［Ｊ］．Ｆｕｎｃｔｉｏｎａｌ＆ＩｎｔｅｇｒａｔｉｖｅＧｅｎｏｍｉｃｓ，２０１４，１４（１）：１１－２２．
［２６］王　莹．云雾贡茶低温胁迫下的转录组分析及ＣｓＰＰＯ基因的功能鉴定［Ｄ］．贵阳：贵州大学，２０２１．
［２７］安逸民．紫花苜蓿响应盐碱胁迫的转录组分析及ＭｓＮＡＣ４７基因的功能研究［Ｄ］．哈尔滨：哈尔滨师范大学，２０２１．
Copyright©博看网. All Rights Reserved.。