2.
3.
4.
Differential Expression of LncRNAs in Nicotiana tabacum var. Samsun NN Infected by TMV at Different Temperatures
JIA HaiYan1, SONG LiYun1, XU Xiang1, XIE Yi1, ZHANG ChaoQun2, LIU TianBo3, ZHAO CunXiao4, SHEN LiLi1, WANG Jie1, LI Ying1, WANG FengLong1, YANG JinGuang12.
3.
4.
责任编辑: 岳梅
收稿日期:2019-09-16接受日期:2019-11-6网络出版日期:2020-04-01
基金资助: |
Editor:
Received:2019-09-16Accepted:2019-11-6Online:2020-04-01
作者简介 About authors
贾海燕,E-mail:
宋丽云,E-mail:jiayoulily2009@126.com。贾海燕和宋丽云为同等贡献作者。
![](https://www.chinaagrisci.com/richhtml/0578-1752/richHtml_jats1_1/images/down.png)
摘要
关键词:
Abstract
Keywords:
PDF (2478KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
贾海燕, 宋丽云, 徐翔, 解屹, 张超群, 刘天波, 赵存孝, 申莉莉, 王杰, 李莹, 王凤龙, 杨金广. 不同温度下TMV侵染枯斑三生烟的LncRNA差异表达[J]. 中国农业科学, 2020, 53(7): 1381-1396 doi:10.3864/j.issn.0578-1752.2020.07.008
JIA HaiYan, SONG LiYun, XU Xiang, XIE Yi, ZHANG ChaoQun, LIU TianBo, ZHAO CunXiao, SHEN LiLi, WANG Jie, LI Ying, WANG FengLong, YANG JinGuang.
开放科学(资源服务)标识码(OSID):
![](https://www.chinaagrisci.com/article/2020/0578-1752/0578-1752-53-7-1381/img_1.jpg)
0 引言
【研究意义】烟草花叶病毒(Tobacco mosaic virus,TMV)对主要经济作物生产造成极大危害且难以防治,来源于野生种黏烟草(Nicotiana glutinosa)的N基因对TMV具有抗性[1]。研究长链非编码RNA(long non-coding RNA,lncRNA)在N-TMV互作引起的枯斑三生烟(N. tabacum var. Samsun NN)抗性反应中的作用,揭示lncRNA在植物系统获得性抗性中的调控功能,可为新型抗病毒技术开发打下基础。【前人研究进展】非编码RNA即不编码蛋白质的RNA,在RNA水平上行使调控功能。LncRNA是长度超过200 nt的非编码RNA,按其与蛋白质编码基因在基因组上的位置,分为基因间长链非编码RNA(long intergenetic ncRNA,lincRNA)、内含子长链非编码RNA(intronic long ncRNA,incRNA)、反义长非编码RNA(antisense lncRNA)等[2]。LncRNA具有同源性较低、组织特异性高等特点[3]。N基因介导的TMV抗性具备温度敏感性,环境温度低于28℃时含有N基因的烟草在TMV侵染后发生过敏反应(hypersensitive reaction,HR)、建立系统获得性抗性(systemic acquired resistance,SAR),环境温度高于28℃时,N基因对TMV的抗性丧失,TMV在植株体内系统侵染[4,5]。SAR的建立依赖多种激素相互协作或拮抗[6,7,8],进而引发一系列信号转导级联反应[9,10],其建立标志是病程相关(pathogenicity-related,PR)基因的大量表达[11]。肖万福对感染TMV后烟草差异表达的基因与蛋白质展开研究,发现它们主要富集在光合作用、能量代谢、植物免疫以及与病原物互作等途径[12]。与此同时非编码RNA的调控作用已得到诸多揭示,最近的研究表明许多lncRNA与植物同病原物的互作相关。番茄黄化曲叶病毒(Tomato yellow leaf curl virus,TYLCV)侵染易感番茄后,2 056个lncRNA参与到植株内的激素信号传导、内质网中的蛋白质加工和植物-病原物互作等多种生物过程中;利用病毒诱导的基因沉默技术降低lncRNA S-slylnc0957的表达后,易感番茄增强了对TYLCV的抗性[13]。此外,寄主lncRNA与病毒产生的小干扰RNA(viral small- interfering RNA,vsRNA)也可能发生相互作用。番茄长链非编码RNA SlLNR1可以通过与TYLCV产生的vsRNA互作来控制病害发展。VsRNA可以诱导易感番茄的SlLNR1沉默,SlLNR1表达下调时番茄叶片呈现发育迟缓和卷曲,而SlLNR1的过表达使TYLCV积累减少[14]。SEO等研究发现,敲除拟南芥长链非编码RNA ELENA1后PR1表达降低,对丁香假单胞菌(Pseudomonas syringe pv. tomato)入侵更敏感;过表达ELENA1的植物PR1表达也升高,对丁香假单胞菌产生抗性[15]。进一步研究揭示ELENA1可以解离FIB2/MED19a复合物,从PR1启动子上的激活因子(MED19a)中释放阻遏物FIB2,最终导致PR1的表达增加[16]。LncRNA在寄主植物与病原物互作中发挥着不可或缺的作用,但对于lncRNA在烟草响应病原物胁迫中的作用少有揭示。【本研究切入点】N-TMV互作引发寄主抗性反应,此过程中差异表达基因和蛋白的功能已有研究,而关于lncRNA在其中的作用鲜见报道。【拟解决的关键问题】利用长链非编码核糖核酸测序(long non-coding RNA sequencing,lncRNA-seq)技术,分析不同温度(25℃和31℃)TMV侵染后枯斑三生烟系统叶lncRNA的表达情况。确定差异表达lncRNA,注释其功能及富集的相关途径,通过对比预测在抗性反应中起作用的lncRNA及作用途径,探究lncRNA在植物系统获得性抗性中的调控功能。1 材料与方法
试验于2018年在中国农业科学院烟草研究所烟草行业病虫害监测与综合治理重点开放实验室完成。1.1 样品准备
供试烟草品种枯斑三生烟,由中国烟草种质资源平台提供;烟草花叶病毒为普通株系(TMV-C),由中国农业科学院烟草研究所植物保护研究中心提供。三生烟生长在昼夜周期为16 h/8 h、25℃或31℃、光合有效辐射为100 μmol·m-2·s-1和相对湿度约为80%的人工气候室中。1.2 样品处理
以磷酸盐缓冲盐水(phosphate buffered saline,PBS)处理为对照,10 mg·mL-1 TMV接种六叶期枯斑三生烟上部第3叶片(图1)。设置4个处理,每个处理包含3个生物学重复,其中处理Ⅰ:TMV+25℃;处理Ⅱ:TMV+31℃;处理Ⅲ:PBS+25℃;处理Ⅳ:PBS+31℃;TMV和PBS处理48 h后,分别将12个样本的上部第2叶片送于北京诺禾致源科技股份有限公司进行lncRNA-seq。图1
![](https://www.chinaagrisci.com/article/2020/0578-1752/0578-1752-53-7-1381/thumbnail/img_2.png)
图1不同温度下TMV接种前后的枯斑三生烟
A:25℃条件下TMV接种前的枯斑三生烟Samsun NN before TMV inoculation at 25℃;B:25℃条件下TMV接种48 h后的枯斑三生烟Samsun NN after 48 hours of TMV inoculation at 25℃;C:31℃条件下TMV接种前的枯斑三生烟Samsun NN before inoculation of TMV at 31℃;D:31℃条件下TMV接种48 h后的枯斑三生烟Samsun NN after 48 hours of TMV inoculation at 31℃。标注为小写字母的图片展示相应枯斑三生烟的接种叶The pictures marked in lowercase letters show the inoculated leaves of the corresponding tobacco
Fig. 1Samsun NN before and after TMV inoculation at different temperatures
1.3 文库制备与高通量测序
利用TRIzol regent(Invitrogen,美国)对4个处理的12个样品进行总RNA提取,分析样品RNA完整性及是否存在DNA污染,每个样本使用3 μg左右RNA,除去rRNA(ribosome RNA)后逆转录合成cDNA第一条链。链特异性文库与普通建库方法的区别在于将第二链cDNA合成反应缓冲液中的dTTP替换为dUTP,对cDNA末端进行修饰后,利用USER(Uracil-Specific Excision Reagent)酶(New England Biolabs,英国)降解含U的cDNA第二链,从而保留第一条cDNA信息进行PCR扩增并获得文库。LncRNA的来源具备链特异性[3],通过与参考基因组比对,链特异性建库可以明确lncRNA正负链来源,获得的基因信息更准确。使用Illumina NovaSeq6000 sequencer(Illumina,美国)按读取长度150 bp进行测序,获得原始序列数据(raw reads)后,经测序错误率分布检查、测序数据过滤,得到有效读段(clean reads)用于随后的比对分析。将测序数据上传至NCBI SRA数据库,BioProject号:PRJNA588726。
1.4 序列比对与拼接
由于枯斑三生烟未进行基因组测序,而近缘品种TN90的基因组信息与枯斑三生烟测序数据比对率高、符合要求,故将其作为参考基因组。参考基因组和基因模型注释文件下载自NCBI(1.5 LncRNA的识别与筛选
利用Cuffmerge(V2.0.2)对拼接得到的转录本进行合并,去除其中链方向不确定的转录本后得到完整转录组信息。从转录组中筛选出外显子数量≥1且长度>200 bp的转录本,用于已知lncRNA(annotated lncRNA,即与数据库注释exon区域有重叠的转录本)的识别和新lncRNA(novel lncRNA)的筛选。利用Cuffcompare将这些转录本与数据库比对,得到已知lncRNA。对剩余转录本,通过Cuffquant筛选出FPKM(expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced,每百万fragments中来自某一基因每千碱基长度的fragments数目)≥2的转录本后[19],利用CNCI(coding-non-coding index,V2.0.0)、CPC2(coding potential calculator)、Pfam-scan(V1.3.0)预测其中没有编码潜能的转录本[20,21,22],取3种工具预测结果的交集作为新lncRNA,而至少被一种软件预测能编码蛋白的转录本作为TUCP(transcripts of uncertain coding potential,编码潜力不确定的转录本)[23]。
1.6 差异表达lncRNA的筛选和聚类分析
通过计算某一lncRNA在不同处理后的FPKM值是否存在差异,判定该lncRNA是否为差异表达lncRNA(differentially-expressed lncRNA,DElncRNA)。通过 StringTie-eB(V1.3.1)定量lncRNA在12个样本中的FPKM值[18],对同一处理的3个重复样本取均值得到lncRNA在该处理下的FPKM信息。筛选枯斑三生烟在25℃环境中差异表达的lncRNA时,将处理Ⅰ设为试验组、将处理Ⅲ设为对照组(TMV_25_vs_PBS_25),利用edgeR(V3.0.8)[24]筛选DElncRNA:当某一lncRNA对应的 |log2 (TMV_25_FPKM / PBS_25_FPKM)|>2、且 校正后的P值(P-adjust)<0.05时,认为该lncRNA在25℃、TMV处理后差异表达。31℃条件下DElncRNA的筛选依据同样方法,以处理Ⅱ为试验组、以处理Ⅳ为对照组(TMV_31_vs_PBS_31)。
为判断不同处理DElncRNA的表达模式,根据DElncRNA在4个处理中的FPKM值,进行表达量层次聚类(hierarchical clustering)分析。
1.7 DElncRNA的实时荧光定量PCR(qRT-PCR)验证
使用Primer Premier(V6.0.0)对各处理DElncRNA进行荧光定量引物设计(表1、表2)。使用ChamQ Universal SYBR qPCR Master Mix(Vazyme,中国)实时荧光定量试剂盒,利用Applied Biosystems 7500 Fast Real-Time PCR System(ABI,美国)进行qRT-PCR反应。总反应体系为20 μL:SYBR qPCR Master Mix 10 μL,PCR正反向引物各0.4 μL,作为模板的逆转录所得cDNA 2 μL,RNase-free water 7.2 μL。扩增条件为95℃预变性30 s;95℃反应10 s和60℃反应30 s,循环40次;95℃反应15 s和60℃反应1 min,95℃反应15 s分析融解曲线。依据融解曲线分析引物的特异性。根据红花烟草肌动蛋白(actin)序列(序列号:GQ339768),设计荧光定量引物(正向引物:5′-TGAGACATTCAACGTTCCGGC-3′,反向引物:5′-GTTAGGTCACGGCCAGCAAG-3′)。以actin为内参,根据2-ΔΔCt法[25]计算基因相对定量的表达,借助SPSS(V24.0.0)进行数据分析,利用Excel 2016绘制图表,得到DElncRNA表达情况,对测序结果进行验证。Table 1
表1
表125℃下差异表达lncRNA引物序列
Table 1
转录本ID Transcript ID | 基因ID Gene ID | 正向引物序列 Forward primer sequence (5′ to 3′) | 反向引物序列 Reverse primer sequence (5′ to 3′) |
---|---|---|---|
LNC_000098 | XLOC_000472 | CCCTCCACGCAGTTCTTCTGA | TGCCGGGAGTTTCGGGTTTA |
LNC_000799 | XLOC_003050 | GCCTTCTGCCGAATTGTTGG | GCTCAAAGGAGGCCAAGTCC |
LNC_001093 | XLOC_004138 | ATGCGTCAGTTTGCATGGGA | GAAATCTCAAGGGTAGCTGTACCA |
LNC_006130 | XLOC_023441 | GCCTTCGACTTGTCGTTTGCT | ATCCGCAAATCAGCCCTTCC |
LNC_012986 | XLOC_049416 | CCTCTGCCTCAGGCTGTCTCG | TTCTGCCGCTTCCGTTATTGCTG |
LNC_019285 | XLOC_073034 | GGATTGAATCTAAGATGAATGCTTGGT | ACAGTCACAAGGGTCTAGAAGG |
LNC_028742 | XLOC_108482 | ACGCTGGCGAACGATAATAGAACC | ACGGTCGAGATTCTCCTGGTCAG |
LNC_037710 | XLOC_141968 | TCTCTTTCTACCGGGAATTTAAAAAGT | CCATTACGTATTGATTGGTGAGCT |
XR_001643474.1 | 107765638 | CCACTTGCATGGGCAGTGAT | GCCACCACACCAACTTCCTC |
XR_001643621.1 | 107766213 | CGGAGAGGTGGAATGGCAAGTG | TCCGAGATGTGCCTCCAAGACC |
XR_001644424.1 | 107769516 | AACATTCACGGCCAGCAGAACC | CTGAAGGCGCGACTGAGATTAGC |
XR_001645141.1 | 107772701 | TCATTCTCAGGTCCTCCGTTGT | TCTCTCTGTGCTTTCCGCCT |
XR_001650509.1 | 107796801 | CGCCCATGTACCCGATTCTG | CTTGTTCCAGATGCGCCAGT |
XR_001654053.1 | 107812073 | CACTGCCTACACAACCTACACGTC | AACAACATGCTGCGAAGAGACTCC |
XR_001657220.1 | 107826167 | AATGGCTTTGTCTGCCTCGG | AGGCGCGTTGATGAAAGAGG |
LNC_001159 | XLOC_004417 | GACACAGGTACTTGGTTCGCTC | AAGAACCGGCACAATACGCC |
LNC_003116 | XLOC_011996 | GGTTTGCAGCGATAACTCGG | TGGCTTCAACTATTCCCGCA |
LNC_004594 | XLOC_017602 | AGTGAATGTAACGGAGGCAGCAAC | ACACAGAAGAGGATCGAGGTCGTC |
LNC_015076 | XLOC_057288 | CTCTTGCAGGAAGGTTGGCT | CCTCATACGGCTCCTCCCTT |
LNC_015077 | XLOC_057288 | CTCTTGCAGGAAGGTTGGCT | CCTCATACGGCTCCTCCCTT |
LNC_028743 | XLOC_108482 | ACGCTGGCGAACGATAATAGAACC | ACGGTCGAGATTCTCCTGGTCAG |
LNC_030105 | XLOC_113591 | TCCCAACGTGTTACATAAGGCA | GCACGGTGTTCTTCGACTCA |
LNC_036321 | XLOC_136711 | CTCTGTCATTCTTCGGTGCCGATG | TGAGCTGCCATGCCATGATACAAG |
LNC_038081 | XLOC_143307 | GCTTAGTGTGTGACTCGTTGGT | TGCGAAGTGATGGGTTGGTT |
XR_001643840.1 | 107767123 | GTTAGTTCGAGGCTGCGGTCAAG | CATCGGCGGCGGACAATTACC |
XR_001646643.1 | 107779549 | GGTTGGATCAGCAGCAGCAGTAG | AGCAGCAGTCTCATTAGCAGCAAC |
XR_001648910.1 | 107789767 | ATCCGAACGACCACTCCCAG | TTGCAAGTCATCACCGTCCG |
XR_001649452.1 | 107792124 | TGTTGGAGAAGGCTGGGACT | CGGCATCTGCTGCTTCTAGT |
XR_001651797.1 | 107802252 | GGCATCAATCACACATGCCGT | AGCAGCCACTGACTTGGAAAC |
XR_001654265.1 | 107813134 | GCACTCGAACCAGAAGGCAA | CAGCAGAAGCTCGACCACAA |
XR_001654587.1 | 107814577 | GTTCATCGCCTTCTTCTGCCTCTC | AAGCACCGAGAGCAACCAACATAG |
XR_001655076.1 | 107816953 | AATGGCTTTGTCTGCCTCGG | AGGCGCGTTGATGAAAGAGG |
XR_001657217.1 | 107826167 | GTGGATGGCGAATGCGATCT | CTGAGGTCCAGCTGCCAATG |
XR_001657841.1 | 107829108 | GTGGTGGTGCATGGCCGTTC | TAGCAGGCTGAGGTCTCGTTCG |
XR_001658459.1 | 107831913 | GTGGTGGTGCATGGCCGTTC | TAGCAGGCTGAGGTCTCGTTCG |
新窗口打开|下载CSV
Table 2
表2
表231℃下差异表达lncRNA引物序列
Table 2
转录本ID Transcript ID | 基因ID Gene ID | 正向引物序列 Forward primer sequence (5′ to 3′) | 反向引物序列 Reverse primer sequence (5′ to 3′) |
---|---|---|---|
LNC_008396 | XLOC_031876 | GCGATCTGCCGAAGCTGTGG | TGCGTTACTCAAGCCGACATTCTC |
LNC_008562 | XLOC_032544 | TCCTGTTGCTTCATTGCTGCT | ACAGATGAACACAGCGCAGG |
LNC_008815 | XLOC_033553 | GACTGTGAAACTGCGAATGGC | GCATCCCTTCCAGAAGTCGG |
LNC_010207 | XLOC_038635 | AGACGAACAACTGCGAAAGCA | CTGGTCGGCATCGTTTATGGT |
LNC_013541 | XLOC_051547 | AGCACCATCCTCACTTCCACCAG | CGCAGGCTGCCTCACAACTTG |
LNC_016042 | XLOC_060869 | GATTAAGACAGCAGGACGGTGGTC | GGCTAGTTGATTCGGCAGGTGAG |
LNC_018096 | XLOC_068634 | AGCCAAGCGTTCATAGCGAC | AACCCAGCTCACGTTCCCTA |
LNC_024626 | XLOC_093025 | ACCACCTGTGGCTCCAGTTACC | TAGGCGGACCTCGCAGAATCTTAG |
LNC_033085 | XLOC_124537 | AGCGAGCAGTCTGGACTCCTATG | ACGAGTAACCAGCGCAATTGGAG |
LNC_033563 | XLOC_126355 | GTGGTGGTGCATGGCCGTTC | TAGCAGGCTGAGGTCTCGTTCG |
LNC_036321 | XLOC_136711 | CTCTGTCATTCTTCGGTGCCGATG | TGAGCTGCCATGCCATGATACAAG |
LNC_037998 | XLOC_143050 | GATTAAGACAGCAGGACGGTGGTC | GGCTAGTTGATTCGGCAGGTGAG |
LNC_038712 | XLOC_145744 | CAACCGCTCAGCCATCTCTC | GGGGTATCCACCCCTATGGC |
XR_001648128.1 | 107786374 | GGAGGGAAGCTGGGTCGTAT | GGCAATTCTCCATCGGCTCC |
XR_001655221.1 | 107817578 | GAGAGTCGGGCTGCAACTTC | GGGAAGCATGGCACCAACAA |
XR_001655293.1 | 107817783 | GTCTGAGGTTGCGGTGAAGGC | ATTGGCGGCAGAGAATTGACAGAG |
XR_001656793.1 | 107824145 | TCTCATGTGGGATGGCACGA | CATTCCACTGCGCACCTCTC |
XR_001657266.1 | 107826393 | TTGAGGCCACTGTTCAAGACTTGG | CGCGAGCAAGGTCATCAGAGC |
LNC_018541 | XLOC_070302 | ACTTCACAAGCAGCAGCTAGTTCC | GTAATGCGCCAGGTGCCGTAG |
LNC_034892 | XLOC_131403 | GTGGGTTCCAGACCACAACC | AGATGGAAGGGGCAGGTGTT |
XR_001644489.1 | 107769825 | CCGTGGACGTGACAACATTGGAG | CGGTGACTGTCGCTGAGATACTTG |
XR_001646112.1 | 107776988 | TCGCTGCTGATTGTTGCTGTCTC | GACTCGTCGGCAACCTCAACTG |
XR_001647877.1 | 107785495 | ATCAATGTGGGAGCTTGCCT | GCTTGTGCGTGTCTCGACTA |
XR_001647999.1 | 107785872 | GCGGATGAGGTATGGTTCACAGC | AGCTGCTTCCATAGTCTGTTGCTG |
XR_001648126.1 | 107786374 | ATGAGAGTGGAGTGGGGCAG | TGCGCCCAATAGGTTATGCG |
XR_001649356.1 | 107791865 | ACGTTGCTTTACACTTTACCCTG | TGGAGGAAGAAACCAAGCGT |
XR_001654070.1 | 107812109 | TCCGTTCAAGTCCACCTGAAGTTC | GCAGCTCACCGTGGAAGTCTC |
XR_001656841.1 | 107824342 | CAAGCCATCATGCCTCACAGT | CCGGAGTCTACCACCCATTCT |
XR_001656939.1 | 107824792 | CTGCTGGTCCAGAAGCCGTTAAG | AATGTCGTCACGAGTTCGGTCATG |
XR_001658441.1 | 107831745 | CAGACGAGGGGTGGAGATGT | GCAAATGCCAACAAACAGGGA |
新窗口打开|下载CSV
1.8 DELncRNA靶基因预测及注释
LncRNA的调控方式分为作用于相邻靶基因的顺式作用模式,以及在其表达水平上识别其他基因的反式作用模式[26]。DElncRNA顺式作用靶基因通过共定位(co-location)分析预测得到,筛选条件为mRNA位于lncRNA上下游100 kb以内。通过共表达(co-expression)分析预测反式作用靶基因[27],采用Pearson相关系数法研究DElncRNA与mRNA的相关性,相关系数绝对值>0.95(P<0.05)时作为反式作用靶基因。研究参考基因组对靶基因的注释,分析其生物学功能。1.9 GO与KEGG富集分析
对DElncRNA的靶基因分别进行GO(Gene Ontology,基因本体论)和KEGG(Kyoto Encyclopedia of Genes and Genomes,京都基因和基因组百科全书)途经分析,GO富集分析通过GOseq(V3.3.2)实现[28],使用KOBAS(V2.0.0)进行KEGG途径分析[29,30],当P<0.05时认为鉴定的GO条目或者KEGG途径被显著富集。2 结果
2.1 不同温度下TMV侵染的枯斑三生烟系统叶lncRNA测序概述
4个处理12个样本获得clean reads均为8 000万条左右(表3),碱基测序错误率(error rate)最高为0.02%。与参考基因组(total mapped)比对发现,95%以上的clean reads被比对到参考基因组,约46%在参考序列上有多个比对位置(multiple mapped),有唯一比对位置的clean reads(uniquely mapped)占49%左右。Table 3
表3
表3测序数据概览
Table 3
样本名称 Sample name | 原始读段 Raw reads | 有效读段 Clean reads | 测序错误率 Error rate (%) | 可定位序列总数及占比 Total mapped | 多定位序列计数及占比 Multiple mapped | 单定位序列计数及占比 Uniquely mapped |
---|---|---|---|---|---|---|
PBS_25_1 | 81285872 | 80011202 | 0.02 | 76771937 (95.95%) | 36818251 (46.02%) | 39953686 (49.94%) |
PBS_25_2 | 79744462 | 78438974 | 0.02 | 75316854 (96.02%) | 36751852 (46.85%) | 38565002 (49.17%) |
PBS_25_3 | 82086686 | 81323058 | 0.02 | 78074493 (96.01%) | 38211722 (46.99%) | 39862771 (49.02%) |
PBS_31_1 | 112775324 | 110071878 | 0.01 | 106555280 (96.81%) | 51925797 (47.17%) | 54629483 (49.63%) |
PBS_31_2 | 110922928 | 108381110 | 0.01 | 105100267 (96.97%) | 51096062 (47.14%) | 54004205 (49.83%) |
PBS_31_3 | 111456946 | 109144098 | 0.01 | 105971578 (97.09%) | 51922061 (47.57%) | 54049517 (49.52%) |
TMV_25_1 | 82779880 | 81084374 | 0.02 | 77865206 (96.03%) | 38330300 (47.27%) | 39534906 (48.76%) |
TMV_25_2 | 81397982 | 79800752 | 0.02 | 76525369 (95.90%) | 36994028 (46.36%) | 39531341 (49.54%) |
TMV_25_3 | 84779126 | 83206166 | 0.02 | 79834756 (95.95%) | 40768694 (49.00%) | 39066062 (46.95%) |
TMV_31_1 | 86454182 | 85145676 | 0.02 | 81734255 (95.99%) | 39993324 (46.97%) | 41740931 (49.02%) |
TMV_31_2 | 82671238 | 81754880 | 0.01 | 78713311 (96.28%) | 35893618 (43.90%) | 42819693 (52.38%) |
TMV_31_3 | 96258450 | 93993096 | 0.02 | 90712104 (96.51%) | 44629779 (47.48%) | 46082325 (49.03%) |
新窗口打开|下载CSV
为研究样本基因表达情况,进行已知类型基因定量分析,约36%的clean reads为protein_coding,lncRNA约占1.4%。经过筛选共获得4 737条已知lncRNA,40 169条新lncRNA;其中lincRNA占47.7%,antisense lncRNA较少,为14.8%,剩余37.5%的lncRNA为intronic lncRNA。
密度分布显示已知lncRNA、新lncRNA及mRNA之间存在一定差异。如图2-A所示,已知lncRNA、新lncRNA相对mRNA长度偏短,已知lncRNA长度集中在600 bp,而新lncRNA在300 bp左右。多数新lncRNA有2个外显子,而已知lncRNA多有3个外显子(图2-B)。两种lncRNA的开放阅读框(open reading frame,ORF)长度都在250 bp以下,新lncRNA相对更短(图2-C)。
3种转录本表达水平差异明显(图2-D),mRNA表达量显著高于TUCP和lncRNA,lncRNA较TUCP更低。对同一处理的重复样品取平均值获得该处理FPKM,图2-E显示4种处理表达水平稍有不同,31℃两个处理转录本表达量高于25℃的两个处理,同一温度条件下TMV和PBS处理间基因表达水平差异不大。
图2
![](https://www.chinaagrisci.com/article/2020/0578-1752/0578-1752-53-7-1381/thumbnail/img_3.png)
图2LncRNA-seq结果分析
A:lncRNA和mRNA序列长度的密度分布Density distribution of sequence lengths of lncRNA and mRNA;B:lncRNA和mRNA外显子个数的密度分布Density distribution of exon number of lncRNA and mRNA;C:lncRNA和mRNA ORF长度的密度分布Density distribution of ORF length of lncRNA and mRNA;D:不同类型转录本表达水平比较Comparison of expression levels of different types of transcripts;E:不同处理转录本表达水平比较Comparison of expression levels of different treatments of transcripts;F:25℃和31℃比较组中差异表达lncRNA数目统计The numbers of differentially expressed lncRNAs in the comparison groups of 25℃ and 31℃
Fig. 2LncRNA-seq results analysis
2.2 DElncRNA的筛选
为确定系统获得性抗性中相关的lncRNA,以25℃和31℃两个温度条件下PBS处理的枯斑三生烟系统叶为对照,筛选TMV处理后差异表达lncRNA,共测得64个DElncRNA(图2-F)。25℃下TMV处理后15个lncRNA上调表达,20个lncRNA下调表达(TMV_25_vs_PBS_25);18个lncRNA在31℃及TMV侵染的情况下表达上调,12个lncRNA下调(TMV_31_vs_PBS_31)。其中LNC_036321在25℃条件下接种TMV后表达下调,在31℃表达上调。结合差异表达lncRNA整体分布情况(图3-A、3-B),DElncRNA的表达在不同处理下存在差异。聚类分析(图3-C)将表达模式相同或相近的lncRNA聚成一组,它们可能具有相同功能或参与同一生物学过程。
图3
![](https://www.chinaagrisci.com/article/2020/0578-1752/0578-1752-53-7-1381/thumbnail/img_4.png)
图3差异表达lncRNA分析
A:25℃差异表达lncRNA火山图,横坐标代表转录本在不同样本中表达倍数变化,纵坐标代表转录本表达量变化差异的统计学显著性,q-value即P-adjust;图中红色点表示有显著性差异表达的上调转录本,蓝色点表示有显著性差异表达的下调转录本Volcano plot of differentially expressed lncRNAs at 25℃. The abscissa represents expression fold change in different samples, the ordinate represents the statistical significance of difference in transcript expression, and the q-value is P-adjust. The red dot in the figure indicates the up-regulated transcript with significant differential expression, and the blue dot indicates the down-regulated transcript with significant differential expression;B:31℃差异表达lncRNA火山图Volcano plot of differentially expressed lncRNAs at 31℃。C:差异表达lncRNA表达量层次聚类,以 lg (FPKM+1) 值进行聚类,红色代表基因高表达,蓝色表示低表达基因,颜色从红到蓝表示lg (FPKM+1)从大到小DElncRNA expression level hierarchical clustering, clustering was performed with lg (FPKM+1) values, with red representing high expression of genes and blue indicating low expressed genes, colors from red to blue indicate lg (FPKM+1) from large to small
Fig. 3DElncRNA analysis
2.3 DElncRNA的表达验证
为验证lncRNA-seq结果的可靠性,利用qRT-PCR对不同温度(25℃和31℃)TMV接种后三生烟的DElncRNA进行相对表达量测定。如图4-A所示,25℃的15个uDElncRNA(up-regulated DElncRNA)中有11个上调表达,其中LNC_000799、XR_001643621.1显著高表达(P<0.05)。对25℃的20个dDElncRNA(down-regulated DElncRNA)定量发现(图4-B)LNC_004594、LNC_038081、XR_001651797.1、XR_001657217.1表达显著下调,共17个与测序结果一致。18个31℃的uDElncRNA中14个表达上调(图5-A),LNC_016042、LNC_024626、LNC_033085、XR_001648128.1显著上调。定量验证的31℃ 12个dDElncRNA中10个lncRNA与lncRNA-seq结果一致(图5-B),XR_001647999.1、XR_001656939.1、XR_001658441.1表达显著降低。80%的DElncRNA的定量结果与lncRNA-seq结果一致,对于不一致的定量结果,除了表达量过低lncRNA(FPKM<1)的测序错误率较高外,与引物设置不合理或者相反链上的基因转录干扰结果也有一定关系[31,32]。图4
![](https://www.chinaagrisci.com/article/2020/0578-1752/0578-1752-53-7-1381/thumbnail/img_5.png)
图425℃条件下差异表达lncRNA 的qRT-PCR验证
A:25℃和TMV处理下uDElncRNA的相对表达量Relative expression of up-regulated DElncRNAs after TMV treatment at 25℃;B:25℃和TMV处理下dDElncRNA的相对表达量Relative expression of down-regulated DElncRNAs after TMV treatment at 25℃ 误差线表示3组生物学重复的正负标准差,采用t-检验分析差异显著性 Error bars indicate the ±SD of three biological replicates. Student’s t-test. *:P<0.1; **:P<0.01; ***:P<0.001。
Fig. 4Validation of differentially expressed lncRNAs at 25℃ by qRT-PCR
图5
![](https://www.chinaagrisci.com/article/2020/0578-1752/0578-1752-53-7-1381/thumbnail/img_6.png)
图531℃条件下差异表达lncRNA 的qRT-PCR验证
A:31℃和TMV处理下uDElncRNA的相对表达量Relative expression of up-regulated DElncRNAs after TMV treatment at 31℃;B:31℃和TMV处理下dDElncRNA的相对表达量Relative expression of down-regulated DElncRNAs after TMV treatment at 31℃
Fig. 5Validation of differentially expressed lncRNAs at 31℃ by qRT-PCR
2.4 DElncRNA靶基因预测
通过共定位和共表达分析预测DElncRNA靶基因,lncRNA可能有多个靶基因,同一基因也可能被不同lncRNA靶向。预测得到的DElncRNA靶基因注释功能丰富,除了与光合作用、植物抗性相关外,还与多种植物激素通路密切联系,也作用于TMV侵染、寄主代谢过程等。例如LNC_037710的靶基因LOC107811086,预测编码抗病性蛋白RPP13(resistance pathogen protein 13);XR_001658459.1靶向的LOC107761822,预测编码丝裂原活化蛋白激酶2(mitogen-activated protein kinase 2)[33];XR_001648910.1的靶基因LOC107789767预测编码过氧化物酶9(peroxidase 9)。LNC_006130靶向的LOC107793779预测编码TMV抗性蛋白N(TMV resistance protein N),另一个靶基因LOC107821288预测编码烟草花叶病毒增殖蛋白3(tobamovirus multiplication protein 3)。值得关注的是,部分可能与激素通路相关的lncRNA在25℃和TMV侵染条件下呈现下调趋势,而在31℃和TMV侵染条件下呈现上调趋势。例如XR_001648910.1(TMV_25_vs_PBS_25)和LNC_008815(TMV_31_vs_PBS_31)同时靶向的LOC107800521,预测编码乙烯响应转录因子ERF003(ethylene-responsive transcription factor ERF003)。两个温度lncRNA的差异表达可能引起靶蛋白相反的积累,进而在25℃的植物抗性建立中发挥作用,而在31℃丧失对TMV的抗性,后续研究中需进一步验证靶基因在不同处理中的表达变化。
2.5 DElncRNA靶基因富集分析
如表4所示,为进一步研究差异表达lncRNA的功能,对其靶基因进行了GO富集分析(P<0.05)。结果显示,膜结合的细胞器(GO:0043227)、质膜部分(GO:0044459)等细胞组分被富集,表明膜组织在枯斑三生烟抗病毒应激反应中具有重要作用。分子功能方面,一些靶基因具备转移酶活性、离子通道抑制剂活性等分子功能,使相应物质得以转运,参与后续应激过程。靶基因也参与到化合物代谢、抗原加工代谢等过程。Table 4
表4
表4GO条目富集的 DElncRNA靶基因数目
Table 4
条目类型 Term type | GO条目 GO term | 比较组合 Comparison group | 基因数 Number of genes |
---|---|---|---|
生物过程 Biological process (BP) | 单组织进程Single-organism process | TMV_25_vs_PBS_25 | 133 |
单组织细胞进程Single-organism cellular process | TMV_25_vs_PBS_25 | 113 | |
氮化合物代谢过程Nitrogen compound metabolic process | TMV_25_vs_PBS_25 | 108 | |
有机环状化合物代谢过程Organic cyclic compound metabolic process | TMV_25_vs_PBS_25 | 105 | |
细胞氮化合物代谢过程Cellular nitrogen compound metabolic process | TMV_25_vs_PBS_25 | 104 | |
单组织进程Single-organism process | TMV_31_vs_PBS_31 | 21 | |
囊泡介导的运输Vesicle-mediated transport | TMV_31_vs_PBS_31 | 20 | |
细胞对刺激的反应Cellular response to stimulus | TMV_31_vs_PBS_31 | 9 | |
自RNA聚合酶II启动子转录Transcription from RNA polymerase II promoter | TMV_31_vs_PBS_31 | 8 | |
端粒维持Telomere maintenance | TMV_31_vs_PBS_31 | 7 | |
细胞组分 Cellular component (CC) | 细胞Cell | TMV_25_vs_PBS_25 | 188 |
细胞组分Cell part | TMV_25_vs_PBS_25 | 188 | |
细胞内组分Intracellular part | TMV_25_vs_PBS_25 | 166 | |
细胞内膜结合细胞器Intracellular membrane-bounded organelle | TMV_25_vs_PBS_25 | 114 | |
膜结合细胞器Membrane-bounded organelle | TMV_25_vs_PBS_25 | 114 | |
质膜组分Plasma membrane part | TMV_31_vs_PBS_31 | 8 | |
微管组织中心Microtubule organizing center | TMV_31_vs_PBS_31 | 4 | |
反式高尔基体转运囊泡膜Trans-Golgi network transport vesicle membrane | TMV_31_vs_PBS_31 | 2 | |
网格蛋白囊膜Clathrin vesicle coat | TMV_31_vs_PBS_31 | 2 | |
网格蛋白外壳Clathrin coat | TMV_31_vs_PBS_31 | 3 | |
分子功能 Molecular function (MF) | 阳离子结合Cation binding | TMV_25_vs_PBS_25 | 69 |
锌离子结合Zinc ion binding | TMV_25_vs_PBS_25 | 41 | |
糖基转移酶活性Transferase activity, transferring glycosyl groups | TMV_25_vs_PBS_25 | 31 | |
己糖基转移酶活性Transferase activity, transferring hexosyl groups | TMV_25_vs_PBS_25 | 23 | |
DNA解旋酶活性DNA helicase activity | TMV_25_vs_PBS_25 | 13 | |
UDP-糖基转移酶活性UDP-glycosyltransferase activity | TMV_31_vs_PBS_31 | 8 | |
抗氧化活性Antioxidant activity | TMV_31_vs_PBS_31 | 8 | |
过氧化物酶活性Peroxidase activity | TMV_31_vs_PBS_31 | 7 | |
作为受体作用于过氧化物的氧化还原酶活性 Oxidoreductase activity, acting on peroxide as acceptor | TMV_31_vs_PBS_31 | 7 | |
翻译监管活动Translation regulator activity | TMV_31_vs_PBS_31 | 5 |
新窗口打开|下载CSV
25℃条件下DElncRNA靶基因富集在色素沉积(GO:0043473)、参与凋亡过程的半胱氨酸型内肽酶活性激活(GO:0006919)等过程,可能与HR形成有关;同时这些靶基因参与蛋白质解聚(GO:0051261)、细胞骨架组织的负调节(GO:0051494)等过程,使细胞膜透性增加、信号分子得以释放。25℃下DElncRNA靶向的基因参与对生长素刺激的反应(GO:0009733),而细胞分裂素代谢过程(GO:00096909)被31℃的靶基因富集。LOC107764060预测被31℃条件下差异表达的lncRNA靶向,其富集在肌醇磷酸代谢过程(GO: 0043647)参与信号转导[34]。两个温度uDElncRNA靶基因都参与到铁载体的合成(GO:0019290)、转运(GO:0015343)、代谢(GO:0009237)等过程,表明铁载体可能在植物应激过程中发挥作用[35]。
KEGG分析(表5)发现,除22条25℃ dDElncRNA靶基因极显著富集(P-adjust<0.05)到植物激素信号传导通路(sly04075)外,31℃ uDElncRNA的15条靶基因也富集(P<0.05)到该通路;分析发现除赤霉素外,水杨酸、乙烯、茉莉酸等激素信号通路均得到富集。进一步研究发现植物病原PAMP(pathogen-associated molecular patterns,病原体相关分子模式)触发免疫通路上的CDPK(sly:101261579)节点,同时被这两个分组lncRNA靶向的基因富集,引起HR反应和细胞壁加固。基于25℃ dDElncRNA和31℃ uDElncRNA存在较多相同靶基因,对比分析发现它们共同富集在ABC运输蛋白,缬氨酸、亮氨酸和异亮氨酸降解,苯丙烷类生物合成,以及mRNA监测途径、嘧啶代谢通路。此外,DElncRNA靶基因也参与到二萜生物合成(sly00904)、磷脂酰肌醇信号系统(sly04070)、水泡运输中的SNARE相互作用(sly04130)、酮体的合成和降解(sly00072)等过程。
Table 5
表5
表5KEGG 通路富集的DElncRNA靶基因数目比较
Table 5
通路 Pathway | TMV_25_vs_PBS_25 | TMV_31_vs_PBS_31 |
---|---|---|
植物激素信号转导Plant hormone signal transduction | 22 | 15 |
苯丙烷类生物合成Phenylpropanoid biosynthesis | 12 | 9 |
mRNA监测途径mRNA surveillance pathway | 11 | 6 |
嘧啶代谢Pyrimidine metabolism | 10 | 6 |
ABC运输蛋白ABC transporters | 3 | 3 |
非同源末端连接Non-homologous end-joining | 3 | 1 |
RNA转运RNA transport | 6 | _ |
碱基切除修复Base excision repair | 5 | _ |
错配修复Mismatch repair | 4 | _ |
氨酰基-tRNA生物合成Aminoacyl-tRNA biosynthesis | 4 | _ |
缬氨酸、亮氨酸和异亮氨酸的生物合成Valine, leucine and isoleucine biosynthesis | 3 | _ |
缬氨酸、亮氨酸和异亮氨酸降解Valine, leucine and isoleucine degradation | 3 | _ |
同源重组Homologous recombination | 3 | _ |
一个叶酸碳库One carbon pool by folate | 3 | _ |
二萜生物合成Diterpenoid biosynthesis | 1 | _ |
苯丙氨酸代谢Phenylalanine metabolism | _ | 7 |
水泡运输中的SNARE相互作用SNARE interactions in vesicular transport | _ | 1 |
磷脂酰肌醇信号系统Phosphatidylinositol signaling system | _ | 1 |
肌醇磷酸代谢Inositol phosphate metabolism | _ | 1 |
酮体的合成和降解Synthesis and degradation of ketone bodies | _ | 1 |
丁酸代谢Butanoate metabolism | _ | 1 |
鞘脂代谢Sphingolipid metabolism | _ | 1 |
新窗口打开|下载CSV
3 讨论
LncRNA已被证明参与植物对病毒、细菌、真菌、根结线虫[36,37,38,39]等多种病原物的抗性反应,作用于植物激素信号转导、系统抗性诱导、过敏反应等多个过程。系统获得性抗性使植物对多种病原物的入侵具备广谱抗性,PR是系统获得性抗性的标志基因[11],一定程度上PR的表达变化揭示系统获得性抗性的建立过程[40]。珊西烟(N. tabacum cv. Xanthi NN)为枯斑三生烟的近缘品种,接种TMV后,珊西烟系统叶SA含量在48 hpi(hours post infection)开始增加、PR1a蛋白在72 hpi开始积累[41]。基于N基因的温度敏感性选定25℃和31℃两个温度条件,在TMV侵染枯斑三生烟48 h后,对其系统叶进行长链非编码RNA测序。LncRNA-seq发现4 737条已知lncRNA、40 169条新lncRNA表达,其中lincRNA和antisense lncRNA占多数。64个lncRNA在不同温度(25℃和31℃)TMV处理后差异表达,qRT-PCR显示这些DElncRNA上 调/下调情况与高通量测序结果基本一致。本研究发现的DElncRNA数目相对于测得的新lncRNA总量而言较少,可能与此时抗性反应较弱有关。DElncRNA在4种处理下的表达模式不同,25℃ TMV处理后15个lncRNA上调表达,20个lncRNA下调表达;18个lncRNA在31℃ TMV处理后高表达,同时12个lncRNA表达降低。对比较得到的DElncRNA顺式作用和反式作用的靶基因展开研究,预测它们编码抗性蛋白、激素信号通路相关蛋白等。部分基因被25℃差异下调lncRNA和31℃差异上调lncRNA同时靶向。考虑到两个温度下枯斑三生烟抗性的差异,认为这些靶基因通过相应lncRNA的调控在烟草抗性反应中发挥作用。靶基因富集的GO条目符合植物感知到病原入侵而启动的一系列应激过程,包括启动离子转运、细胞膜透性改变、激素合成转运等[42]。KEGG途径分析显示两组靶基因都富集到植物激素信号传导,ABC运输蛋白,缬氨酸、亮氨酸和异亮氨酸降解,苯丙烷类生物合成等通路。
DElncRNA富集到水杨酸、茉莉酸、乙烯等多个激素信号通路,水杨酸与HR的形成密切相关[6,7],是重要的内源信号分子,茉莉酸、乙烯也作用于SAR的建立[43,44]。DElncRNA参与调控多种激素通路,使烟草得以建立多重抗性机制,抵御多种病原物的再次入侵[11]。ABC运输蛋白与分子运输过膜有关[45],亮氨酸、异亮氨酸、缬氨酸响应脱落酸介导的渗透胁迫[46],苯丙烷类化合物也被证明与烟草抗病性相关[47],这些通路在两个温度条件下同时富集表明它们可能在枯斑三生烟的抗性反应中发挥作用。
LNC036321在25℃及31℃条件下同时差异表达,它共定位的靶基因LOC107805368预测编码温度诱导的载脂蛋白(temperature-induced lipocalin-1,TIL1),同样富集在激素信号传导通路。进一步的研究发现接种TMV 6 d内LNC_036321与TIL1 mRNA的相对含量呈现此消彼长的趋势,表明二者可能存在负调控关系。LncRNA调控方式多样,依照功能可将其分为4种:吸附RNA或蛋白与之结合的诱饵分子(decoys)、通过改变自身转录调节相关基因表达的信号分子(signals)、结合蛋白形成复合物指导其移动到特定位置的导向分子(guides)以及作为RNA及蛋白等分子结合平台的骨架分子(scaffolds)[27]。LncRNA16397在致病疫霉抗/感番茄品种间差异表达,SlGRX22是它的共定位靶基因,lncRNA16397作为SlGRX22反义转录物调节其表达,使活性氧积累减少、细胞膜损伤减轻,增强番茄对致病疫霉的抗性[38]。正义转录本(sense transcript)的表达水平随正义和反义对之间重叠的增加而显著降低[31]。lncRNA16397作为“信号分子”,对SlGRX22表达进行顺势调控。而LNC_036321与LOC107805368序列相似性不高,考虑到LNC_036321为lincRNA,其更有可能以“诱饵分子”或“导向分子”的方式对LOC107805368进行调控[48]。
载脂蛋白一般在植物对温度等非生物胁迫的应激中发挥作用[49],但SADE等研究发现番茄对病毒抗性/易感载脂蛋白(Solanum lycopersicum virus resistant/ susceptible lipocalin,SlVRSLip)与TYLCV抗性相关[50],需要通过进一步的实验明确TIL1是否参与枯斑三生烟对TMV的抗性反应。
本研究仅对接种TMV 48 h的枯斑三生烟系统叶进行了lncRNA测序和分析,此时DElncRNA数目较少,这些lncRNA通过调控植物激素信号传导等过程参与枯斑三生烟抗性反应。为找到更多与系统获得性抗性相关的lncRNA,需要对TMV侵染更长时间的系统叶进行lncRNA测序分析,配合接种叶的lncRNA表达情况,对枯斑三生烟抗性反应中的lncRNA进行更全面深入的分析。
4 结论
不同温度(25℃和31℃)条件下对枯斑三生烟接种TMV后,4 737条已知lncRNA、40 169条新lncRNA在系统叶表达,其中64个lncRNA在不同处理间差异表达。利用qRT-PCR技术验证了DElncRNA的特异性与可靠性,并通过靶基因预测和功能分析,发现这些DElncRNA广泛参与到激素信号传导、物质转运等过程,在枯斑三生烟的抗性反应中发挥作用。研究结果可为N-TMV互作过程中关键 lncRNA的筛选和功能研究提供理论支持。参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子
[本文引用: 1]
[本文引用: 1]
[本文引用: 2]
[本文引用: 1]
[本文引用: 1]
[本文引用: 2]
[本文引用: 2]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 3]
[D].
[本文引用: 1]
[D].
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
URL [本文引用: 1]
[本文引用: 2]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 2]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 2]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 2]
[本文引用: 1]
[D].
[本文引用: 1]
[D].
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]