2021 年 1 月
雪特兰
第44卷第1期湖南师范大学自然科学学报Journai of Naturai Science of Hunan Normai University Vvl44 No.1Jan. ,2021DOI ： 104612/j4ssn4096W281402141405
藜麦叶片响应低氮胁迫基因的可变剪接事件分析
杨利艳1，杨小兰1，朱满喜1，张玉荣1，杨雅舒1，王创云"*
*收稿日期：2020-08-09
基金项目：山西省应用基础研究计划项目(201901D111285) +中国农业科学院科技创新工程协同任务资助项目(CAAS-
XTCX20190025-3)
* 通信作者,E-mail ：wrocy@
(1.山西师范大学生命科学学院，中国临汾041000 + 2.山西农业大学农学院，中国太谷030800)
摘 要 缺氮环境对藜麦生长不利，为鉴定和分析藜麦响应低氮的可变剪接基因，采用Lymma Hi S eq TM
分别对 低氮处理(氮浓度丄42 mmol-L g )的藜麦叶片和对照组(氮浓度：7.1 mmoi-L-1)藜麦叶片的转录组进行测序。结 果表明，对照组(CK )中26 094个基因发生可变剪接事件68 764个，低氮处理组(LN )中25 562个基因发生可变剪 接事件66 611个。两组处理的可变剪接事件均以第一个外显子可变剪切！TSS )和最后一个外显子可变剪切(TTS ) 两种类型为主。对比两组发生可变剪接的基因发现LN 组中有1 041个新基因发生可变剪接，存在147个显著差异 表达基因。GO 功能富集分析表明，发生可变剪接的差异基因主要富集于代谢过程、刺激反应、转移酶活性及水解 酶活性等功能上。KEGG 通路分析表明，可变剪接差异基因显著富集到DNA 复制、异源末端链接、抗坏血酸和醛酸 代谢、错配修复系统等通路上，表明这些差异可变剪接基因参与了藜麦低氮胁迫的响应。
俚语是什么意思关键词藜麦；低氮；可变剪接事件
中图分类号S519 文献标识码 A 文章编号2096-5281( 2021) 01-0041-06
Alternative Splicing Events Analysis of Low Nitrogen
Strbss Respon Genes 0 Quinoa Leaves
E4NG Li-yan 1, E4NG Xiao-lan 1, ZHU Man-xi 1, ZHANG Yu-rong l , E4NG Ya-shu 1, WANG CCuang-yun 1 *
(1. Colleac or Life Science , Shanxi Normai University , Linfen 041000, China ；
2. Colleac of Agriculture , Shanxi Agriculturai University , Taiyu 030800, China )
Abstract The lack of nitrogen is unfaverabie for quinog growth. To identify and anHyze alternative splicer genes in respon to low nitrogen, Ulimina Hiq™ was applied to quence transcoptome of quinod leaves treated with low nitrogen ( N concentration : 1.42 mmco • L 1 ) and tho treated witU low nitrogen ( N concentration : 7.1 mmol ■ L 1) , respectively. Our results show that 68 764 alternative splicing events occurred in 26 094 genes in the conieoageoup ( C_) and 66 611 aaieenaioeeCpaocongeeenioccu e ed on 25 562 geneCon iheaownoieogen ieeaimeni group ( LN) . The two groups of splicing event were mainly of tUe alternative 5’ first exon (TSS) and alternative 3’ last exon (TTS ) . Comparing tUe two groups of genes with vayablv splicing , 1 041 new genes witU vayablv spli cing were found in LN group , and 147 genes were significantly dlferential i y exposd. GO functional enochment a nalysis results unveil tUat tUe dmarentially expresd genes witU vayablv splicing were mainly concentrated in the metabolic process , stiniulHion reaction , transfera activity and hydrola activity. KEGG patUway analys results indicate that tUe split-dOferentialiy expresd genes were signigcantly enoched in DNA replication , Hlocenic termi - 会计属于什么行业类别
naaaonks , asFo eb o and aadehydemeiaboaosm , and mosmaih eepaoesysiems , suggesiongihaiihespaoi-do t eeenioa a y expresd genes were involved in tUe respon to low nitrocen stress in quinoc.
Key wordt quinoc ； low nitrocen + alternative splicing events
42湖南师范大学自然科学学报第44卷
藜麦！Chenopodium quinoa Willd）为一年生双子叶植物，苋科藜亚科藜属，是一种起源于南美洲安第斯区的古老物种⑴（藜麦种子富含无谷蛋白以及人体无法合成的8种必需氨基酸，且铁、钙、磷含量丰富，具有良好的营养健康效益［2,3］（藜麦不仅有极高的营养价值，还具有很强的抗逆性，因此受到世界的关注并被广泛引种，是未来最具潜力的农作物之一⑷。
氮元素（N）是细胞结构不可或缺的一部分，担负着植物的生长发育及机体调控的重要责任［5］（氮素约占植物干重的0.3%~5%，缺氮或低氮是导致农作物生长不良、减产的重要原因.句（氮胁迫条件下，植物会通过根的伸长来增强氮的吸收［7］，在甘蓝型油菜.8/、小麦［9］、水稻［10］中皆有报道。然而目前对藜麦低氮响应的研究鲜有报道。
可变剪接（altewative splicing,AS）是一种转录后调节机制。该机制可使蛋白的最终产物表现出不同
的功能和结构特性，或者在相同的细胞中可能由于表达水平的不同而导致表型发生变化，因此对增强机体内转录组的可塑性、基因的表达以及蛋白质结构的多样性起重要作用［11,12］（目前，利用RNA-q技术对可变剪接的分析有不少报道，可变剪接基因（alteoativela spliced genes,ASG）在一些生物中占较高比例，如人（95%）、线虫（71%）［13,14］，在禽类.15/、甲壳类动物［16〕中也发现了新的剪接变体。作为模式植物，拟南芥AS 事件的发生比例也较高，为61%"］，约有42%的内含子参与拟南芥的生长发育及胁迫响应［18,19］O KHOKHAR等.20］发现在不同生态类型的拟南芥中，许多反式剪接数量性状位点（sqtl）在昼夜节律钟、开花和应激反应基因中富集，表明差异可变剪接在调节这些重要过程中发挥了作用。KANNO等.21］在拟南芥中通过对一个po4Ca突变体的定量磷酸化蛋白组学分析，检测到几个富含丝氨酸/精氨酸的蛋白的磷酸化变化，这些蛋白调节组成型和选择性剪接以及其他剪接相关的因素。因此，研究可变剪接有利于对基因的深入了解。本研究对不同氮处理的藜麦叶片进行转录组测序，统计并分析可变剪接基因，旨在挖掘响应低氮胁迫的ASG，从而探究AS对植物低氮响应基因的调控作用。
1材料与方法
1.1材料
藜麦ZK7种子来自于山西省农业科学院作物科学研究所。将种子浸泡在10%NaClC溶液中消毒15 mC，
用蒸馏水冲洗3~5次，置于放有湿润滤纸的培养皿中，在25l下培养发芽。待芽长约0.5cm时，将其移入干净的砂子中置于光照培养箱育苗，光暗周期16h/8h，温度18-22l。待藜麦幼苗生长至2叶，采用霍格兰全营养液培养，每3天浇一次，长至6叶期进行以下处理。
1.2实验设置huge
选取长势相同的幼苗进行两组处理，即对照（CK）：全营养液培养（N浓度为7.1mmol-L-1）；低氮培养（LN）：氮元素浓度为全营养液的20%（N浓度1.42mmomL-）o每个处理设置3个重复，处理5天后取相同叶位材料，液氮速冻后保存。提取总RNA，转录组数据委托生物公司测定。
1.3总RNA的提取与质量检测
采用试剂盒RNAprep Pure Plant Kit（TIANGEN）的方法，分别提取藜麦叶片的总RNA。用Nanodop ND-1000检测RNA的浓度Labchip GX精确检测RNA的纯度和完整性（RIN）（
1.4 RNA-q文库构建和测序
RNA-q文库构建和测序由北京百迈客生物科技公司完成，将对照组！CK）和低氮胁迫组（LN）两个样品检测合格后构建文库，之后进行Alumina Hi Se q T M测序。测序去除含有接头和低质量的Reads，获得高质量Clean数据。本研究通过分析已获得的Clean数据鉴定在低氮胁迫条件下发生AS事
件的差异表达基因，从中找寻CK和LN中ASG的共性和差异性，进而理解AS对低氮胁迫的调控作用。
1.5可变剪接事件的鉴定
采用StUngTie.22］对Hisat2的比对结果进行拼接，通过ASpofile［23］软件分析可变剪接类型，细分为12类，并获取CK和LN中相应的可变剪接类型及表达量.24］（
1.6低氮胁迫相关基因GO分类及KEGG功能富集分析
为筛选响应低氮胁迫的ASG，对LN中特有ASG进行差异分析［25］O并以Fold Change（表达量差异倍数）%2,FDR（错误发现率）<0.01且q<0.05（q为校正后的P值）为筛选标准，由于不同处理组基因存在差
tear什么意思
第1期杨利艳等：藜麦叶片响应低氮胁迫基因的可变剪接事件分析43异，根据CK 和LN 组中基因的表达量值绘制各样本中的基因表达量热图，采用Heatmap 软件作图。GO 分类 根据 GO 数据库(Genc Ontology Consortium "获得分子功能(moleculaT function )、生物学过程(biolooical process )和细胞组成！ cellular component )的 GO 编号，功能富集采用 cluster Profiler 分析！ q<0.05 )。KEGG 通 路分析以Pathway 为单位，应用超几何检验，统计并分析其与整个基因组背景相比显著富集的Pathway o 2结果与分析
2.1 RNA 的提取与检测
将提取的 RNA （ RIN 值％7.0）采用 Calipcr Lab Chip GX 检测,28S/18S 均在 2.0〜2.1 ,A 16Q /A^Q 均在 2.1〜
2.2（图 1 ）（
图1 CK 与LN 总RNA 的GX 结果
Fig. 1 The GX results of total RNA in CK and LN
2.2可变剪接事件的鉴定
对CK 和LN 两个样品的转录组数据分析后发现CK 中26 094个基因发生了 AS ，共发生68 764个,LN
中25 562个基因发生AS ，共发生66 611个。
在各种AS 类型中，CK 组以TSS 和TTS 两种类型为主，分别占所有可变剪接事件的40.95%和40.87%； 此外AE,IR 及SKIP 三种类型分别占所有剪接事件的7.06%,3.97%及3.98%；剩下7种类型共占3.17%。低 氮处理组中,LN 组的可变剪接事件仍以TSS 和TTS 两种类型为主，但较CK 组有所增加，分别占41.45%和 41.11%；AE,IR 及SKI 三种类型分别占6.18% ,4.58%及3.68% ，其中IR 剪接类型较CK 组增加15.4% ；剩 下7种类型共占3.05%(表1 )(
表1 CK 和LN 基因组中鉴定的可变剪接事件类型及数量
Tab. 1 Types and number of alternative splicing events identified in CK and LN
组别
AE IR MIR MSKIP SKIP TSS TTS XAE XIR XSKIP XMSKIP XMIR CK 4 8562 7322823142 740
28 15628 1063825145805646LN
4 1183 0502742942 42027 60927 3873295085106646对CK 与LN 中发生的ASG 进行整理对比后发现，有24 520个
AS 发生重合。然而，有1 573个新基因仅在CK 中发生AS o 同时，
LN 组中有1 041个新基因发生AS ，在CK 中却没有发生(图2)(
2.3 LN 特有可变剪接基因差异表达分析
对LN 中1 041个ASG 进行分析，差异表达显著的有147个基
因，对其聚类分析见图3(这147个基因中AS 类型以TSS 和TTS 为
主，其中，上调基因有72个，下调基因有74个，还有一个基因(NW_
018743732.1)在TSS 事件中下调，在TTS 事件中上调。表明藜麦在
低氮条件下的基因表达与这两种可变剪接事件密切相关。同时，这
147个可变剪接基因的功能主要有无机离子运输和传递，次生代谢
物的合成(NW_018742204.1 及 NW_018742461.1 )、DNA 解旋酶及
修复蛋白(NW_018742412.1)、翻译后修饰(NW_018742673.1 )、信 号转导(NW_018743014.1)、细胞骨架的形成(NW_018742814.1)及 Fig. 2 Number of overlapped and specific alterna tive splicing genes in CK and LN
44湖南师范大学自然科学学报第44卷细胞周期调控(NW_018742966.1)
2.4 LN 中差异可变剪接基因GO 注释及功能分析
对LN 基因组中147个可变剪接基因进行GO 注释
和功能分析，发现39个GO 注释的基因(图4)。其中，
显著富集(q<0. 05)的GO 项在生物学过程中有81个,
在细胞构成中有26个，分子功能中有60个。在生物学
过程中，发生可变剪接的差异基因主要集中于代谢过
程(GO ：0008152)和刺激反应(GO ： 0050896)；在细胞
构成中，发生可变剪接的差异基因主要集中于核(GO ：
0005634)、细胞质(GO ： 0005737 )、细胞壁(GO ：
0005618)以及质外体(GO ： 0048046)上；在分子功能
中，发生可变剪接的差异基因主要富集于转移酶活性
(GO : 0016740)和水解酶活性(GO ： 0016787 )相关
GO 项。
2.5差异可变剪接基因KEGG 通路分析
对LN 中147个差异可变剪接基因进行KEGG 通
路分析(表2)，发现显著富集到DNA 复制、异源末端链
接、抗坏血酸和醛酸代谢、错配修复系统等通路上。
DNA 复制、异源末端链接参与新细胞的合成；抗坏血酸 和醛酸代谢参与植物逆境响应、细胞生长；错配修复系
CK LN 图3 LN 差异可变剪接基因表达聚类热图
3 Cluster thermography of dimeontiaVy expresd genes 统参与DNA 双链上错配碱基对的矫正，泛素介导的蛋白水解调控细胞周期，囊泡运输的相互作用参与分子 的转运，谷胱甘肽代谢参与抗氧化和解毒作用(
A Biological Process
B Cellular Component
C Molecular Function ■ multicellular organismal process (2)■respon to oxidative stress (2)■ cellular process (8)
lowell广州新东方英语■ respon to stimulus (8)■ biological regulation (6)■ reproductive process (6)■ biosynthetic process (6)■ nucleus (4)■ cytoplasm (3)
■ cell wall (3)■ apoplast (3)■ integral component of membrane (2)■ chloroplast envelope (2)■ transfera activity (7)■ hydrola activity (6)■ oxidoreducta activity
(4)
DNAbinding (4)
■ metal ion binding (4)
■ heme binding (3)
■ methylation (2)■ cell killing (1)■ localization (1)■ single-organism process (1)
■ protein phosphorylation (1)
旗开得胜英文版■ regulation of tanscription,
DNA-templated (1)
■ signaling (1)■oxidation-reduction process (6)■ developmental process (4)
■ cell cycle (4)
■ reproduction (4)■ cellular component organization or biogenesis (4)■ extracellular transport (3)■ immune system process (2)■ macromolecular complex (2)1■ chloroplast (1)■ cytoplasmic part (1)■ spindle (1)■ intracellular part (1)■ plasma membrane (1)■ mitochondrion (1)■ membra
ne (1)■ methyltransfera
activity (1) electron carrier activity
(1)■ liga activity (1)■ endonuclea activity (1)■exonuclea activity (1)■ cation binding (1)■ iron ion binding (1)■ nucleotide binding (2)■ ATPbinding (2)■ protein kina activity (2)■ peroxida activity (2)
plating■ isomera activity (2)■flavin adenine dinucleotide binding (2〕■ lya activity (1)■ mangane ion binding (1) ■ transporter activity (1)
注：数字代表在该功能下注释的基因数;A ：生物学过程；B ：细胞构成；C ：分子功能。
图4 LN 中特有可变剪接差异基因的GO 注释
Fig. 4 GO annotations for specific dimeontiat expresd alternative splicing genes in LN
3讨论
选择性剪接在广大真核生物产生潜在的基因组信息方面发挥着重要作用，它是一个重要的转录后调节
机制，可以增加蛋白质多样性，影响mRNA 的稳定性：26/ o 然而，关于AS 如何参与藜麦对低氮胁迫的响应, 人们知之甚少。在本实验中，笔者在LN 组中共识别到66 611个AS 事件，包括TSS ,TTS ,AE,IR ,SKIA 等12
第1期杨利艳等：藜麦叶片响应低氮胁迫基因的可变剪接事件分析45
种。其中TSS及TTS事件最为丰富，其次是AE,IR和SKIP。与CK组相比，鉴别出低氮处理组(LN)基因组中有1041个新基因发生AS，这些基因的生物学过程富集在代谢过程和刺激反应等过程中。LEVIUTAN 等［27］发现拟南芥在低温胁迫下，其表达水平保持相对不变，但其剪接方式在异构体的相对丰度上有明显变化。大多数低温调控的选择性剪接将提前终止密码子(PTC)引入到转录本中，通过mRNA衰变(NMD)过程产生潜在的降解目标，这可能导致新蛋白质功能改变或产生负面影响。FANG等［28］发现水稻中11个OsAAT选择性剪接基因，其中1个基因剪接变异的表达受氮(N)的调控,5个基因在各种形态氮浓度较高时表达上调,9个基因在硝酸盐浓度较高时表达上调。剪接变异体的mRNA水平在水稻叶和穗发育过程中也受到调控，结果表明水稻OsAAT家族的可变剪接变异体在不同氮素条件下的表达率均有变化。本实验发现，藜麦响应低氮的差异可变剪接基因富集于与转移酶活性和水解酶活性相关的通路。GIURETTA等［29］研究表明由GGT1和GGT2编码的Apoplastic-谷氨酰胺转移酶活性对植物营养和生殖发育具有重要意义，因此笔者推测转移酶对藜麦缓解低氮胁迫起一定作用。对LN中差异可变剪接基因进行KEGG分析，发现显著富集到DNA复制、异源末端链接、抗坏血酸
和醛酸代谢、错配修复系统等通路上，这些信号通路与藜麦氮吸收密切相关，表明在藜麦生长过程中可能通过差异可变剪接基因调节氮吸收来响应低氮胁迫。
表2LN与CK差异可变剪接基因KEGG通路分析
Tab.2Analysis of KEGG p athway in dimeynhaHy expresd altemaHve splicine genes between LN and CK KEGG通路基因编号3值核苷酸切除修复NW_018742314.1,NW_018744348.10.06 DNA复制NW_018742556.1,NW_018742616.1,NW_018744348.10.0039
异源末端链接NW_018742616.1,NW_018742733.10.00046
维生素B6代谢NW_018742684.10.072同源重组NW_018742733.1,NW_018744348.10.049
嘌吟代谢NW_018742738.1,NW_018744784.10.13 ABC转运蛋白NW_018743212.10.093苯丙素的生物合成NW_018743975.10.60
DNA错配修复系统NW_018744066.1,NW_018744348.10.032
抗坏血酸和醛酸代谢NW_018744240.1,NW_018745650.10.025谷胱甘肽代谢NW_018744240.10.34
脂肪酸生物合成NW_018744578.10.17丙酮酸代谢NW_018744578.10.32碳代谢NW_018744578.10.67脂肪酸代谢NW_018744578.10.25
二萜生物合成NW_018744718.10.087囊泡运输的相互作用NW_018744755.10.17
泛素介导的蛋白水解NW_018744955.10.13
植物激素信号转导NW018745291.10.57
4结论
本试验利用转录组样品的高通量测序数据鉴定了藜麦响应低氮的特异可变剪接的差异表达基因，藜麦叶片主要通过TSS和TTS两种剪接方式形成差异基因，其中NW_018742556.1,NW_018742616.1,NW_ 018744348.1,NW_018742616.1,NW_018742733.1等基因通过DNA复制和异源末端链接等通路来响应低氮。
参考文献：
［1］林春，刘正杰，董玉梅，等.藜麦的驯化栽培与遗传育种［J］.遗传,2019,41(11)：1009-1022.
［2］HINOJOSA L,SANAD M,JARVIS D,et al.Impact of heat and drought stress on peroxisome proliferation in quinoa.J/.Plant J,2019,99(6):1144-1158.
皮肤美白小窍门.3］申瑞玲，张文杰，董吉林，等.藜麦的营养成分、健康促进作用及其在食品工业中的应用［J］.中国粮油学报，2016,31
