作者简介:马 轩(1978—),男,深圳大学博士后研究人员. 研究方向:植物miRNA介导的翻译抑制机理. E-mail: xuanma@genetics.ac.cn
中文责编:晨 兮; 英文责编:艾 琳
1)深圳大学生命与海洋科学学院,广东深圳518060; 2)中国科学院遗传与发育生物学研究所,北京100101; 3)中国农业科学院深圳农业基因组研究所,广东深圳518120; 4)美国加利福尼亚大学河滨分校植物科学系, 加利福尼亚河滨CA 92521,美国
Ma Xuan1,2, Li Shengben3, Mo Beixin1, Cao Xiaofeng2, and Chen Xuemei1,41)College of Life Sciences and Oceanography, Shenzhen University, Shenzhen 518060, Guangdong Province, P.R.China2)Institute of Genetics and Developmental Biology, Chinese Academy of Sciences,Beijing 100101, P.R.China3)Shenzhen Institute of Agricultural Genomics, Chinese Academy of Agricultural Sciences,Shenzhen 518120, Guangdong Province, P.R.China4)Department of Botany and Plant Sciences, Institute of Integrative Genome Biology, University of California, Riverside, California 92521, USA
phytobiochemistry; microribonucleic acid(miRNA); ago1-27; RNA-seq; degradome; target genes; Arabidopsis
DOI: 10.3724/SP.J.1249.2017.01027
miRNA在植物生长发育和环境适应等方面起着极为重要的作用.本研究对拟南芥ago1-27突变体进行了RNA-seq,鉴定到几千个差异表达基因,并且miRNA靶基因倾向于在ago1-27中表达上调.ago1-27 RNA-seq结合降解组分析鉴定到新的miRNA靶基因,如miR396靶向半胱氨酸蛋白酶基因和miR167靶向编码AT-hook DNA结合蛋白的基因等.
miRNAs play important roles in plant development and adaptation to environment. We conduct the RNA-seq of Arabidopsis ago1-27 mutant and has identified thousands of differentially expressed genes. The study shows that miRNA targets tend to be upregulated in the ago1-27 mutant. New miRNA targets are identified by ago1-27 RNA-seq combined with degradome analysis. For examples, miR396 targets cysteine proteinase genes and miR167 targets a gene encoding AT-hook DNA-binding protein.
小核糖核酸(microribonucleic acid, microRNA或miRNA)是一类长度约为20~24个核苷酸的内源非编码RNA.miRNA与Argonaute(AGO)蛋白相结合,对靶RNA进行RNA剪切或翻译抑制[1- 4].植物miRNA在植物形态建成、生长发育和环境适应等方面起着极为重要的作用[3- 4].自2002年首次证实植物中存在miRNA[5]以来,关于植物miRNA的研究进展异常迅速.植物miRNA和靶RNA之间互补紧密,因此,miRNA的靶基因预测较为容易.模式匹配程序scan_for_matches最早被应用于植物miRNA靶基因预测[6].之后,以序列匹配程序FASTA为基础的工具,如Targetfinder和psRNATarget能够快速预测植物靶基因[7-8].降解组(degradome)测序使miRNA靶基因的大规模鉴定成为可能,目前已广泛用于各种植物miRNA靶基因的鉴定中[9-10].
目前已经确认的miRNA靶基因,主要是互补紧密的类型,并集中在高度保守的miRNA上.本研究对模式植物拟南芥的ago1-27突变体进行了RNA-seq,得到了新的候选miRNA靶基因.结合降解组测序数据分析,鉴定出新的miRNA靶基因.
拟南芥野生型(wild type, WT)和ago1-27突变体的花序用于RNA-seq.使用TRIzol(Invitrogen)提取花序的总RNA.每个样品制备3个生物学重复.使用TruSeq试剂盒(Illumina)构建RNA-seq文库.在IlluminaHiSeq 2500平台上进行50-nt单端非链特异性RNA-seq.原始测序数据已经提交至NCBI SRA数据库(数据号GSE77211).
使用Tophat(v2.0.11,默认参数)将RNA-seq数据与拟南芥基因组(TAIR10)比对,唯一定位的读段(reads)用于后续的分析.应用Cuffdiff进行基因表达分析,采取的标准为:① FDR(false discovery rate)<0.05; ② 基因差异表达>2倍.拟南芥miRNA数据源于miRBase(v21),已验证的miRNA靶基因数据源于TarBase(v6.0).miRNA靶基因预测应用TarHunter-L(https://github.com/XMaBio/TarHunter-L).降解组数据应用CleaveLand4进行分析,采用的标准是:① 0类或1类; ② 罚分<5.0; ③ 剪切位点reads丰度≥10×TP10M,TP10M(transcripts per ten million)为剪切位点reads数占每千万总reads数的比例.
对拟南芥野生型和ago1-27突变体植物进行RNA-seq,每个样品进行3个生物学重复.首先比较3种分析RNA-seq数据的流程方法:Tophat2-Cuffdiff、Tophat2-HTSeq-edgeR和Tophat2-HTSeq-DESeq.采用的标准是:① FDR(或调整的P值)<0.05; ② 基因差异表达>2.0.Cuffdiff、DESeq和edgeR 3种流程分别鉴定到1 862、1 823和2 387 个差异表达基因,大多数差异表达基因都能被3种方法检测到,如图1(a).
其中,edgeR具有较高的灵敏度,DESeq具有最高的特异性.对这3种方法进行两两比较,基因差异表达的倍数(FC为ago1和WT基因差异表达倍数)变化相关性很高(Pearson相关系数> 0.99),如图1(b)和(c),说明该分析方法可靠.因为Cuffdiff方法平衡了灵敏性和特异性,以下结果将基于Cuffdiff的分析方法.同时对3个生物重复进行了比较,3个生物学样品具有高度可重复性(Pearson相关系数> 0.996,如图2(a)和(b),这充分说明了分析样品的可靠性.
观察到在ago1-27突变体中AGO1基因的表达值升高,如图3,RPKM(reads per kilobases per million reads)为基因表达值单位.对于这一现象,笔者的解释是:首先, AGO1是miR168的靶基因; 在ago1-27突变体中,miRNA与内质网的关联程度下降[11],而内质网膜是miRNA作用的重要场所[11-12],因此在ago1-27中miRNA靶基因的表达失调.应用Cuffdiff在ago1-27突变体中鉴定到689个上调基因和1 175个下调基因,上调和下调的基因中包含24和5个已经验证的miRNA靶基因.这个结果显示,已经确认的miRNA靶基因倾向于富集在ago1-27上调的基因中(Fisher's exact test,P值= 6.97×10-7).
应用TarHunter-L(见第1章“材料与方法”)预测拟南芥miRNA在编码序列上的靶位点.如图4(a),罚分≤3.0的靶基因在ago1-27突变体中的表达比在野生型的高,而非靶基因在ago1-27中的表达与野生型相当.在ago1-27中上调的预测的靶基因(罚分≤3)中,17个没有报道过(表1).这17个基因包括5个miR396靶基因,还有一些保守miRNA(如miR172和miR167等)的靶基因.
图4 ago1-27 RNA-seq结合降解组鉴定到新的miRNA靶基因
Fig.4 ago1-27 RNA-seq combined with degradome analysis identified novel miRNA targets
为进一步确认通过ago1-27 RNA-seq分析得到的靶基因,对已发表的降解组数据(包括花序和叶片等4个测序文库)进行了分析,鉴定到22个miR396靶基因,包括一些新的靶基因,如FLU、LSU2、RD21B和AT4G32250等.miRNA在靶位点的剪切效率以TP10M表示.在这22个miR396靶基因中,5个没有在ago1-27 RNA-seq中检测到,其余17个基因中的13个在ago1-27突变体中上调变化> 1.5倍,因此,ago1-27 RNA-seq结合降解组分析能有效鉴定miRNA靶基因.
通过ago1-27 RNA-seq分析,还鉴定到AHL11(AT3G61310)是miR167潜在的靶基因,因为它在ago1-27中的表达上调1.6倍.对降解组数据的分析显示,AHL11确实是miR167的靶基因,如图4(b).另外,鉴定到一个果胶甲酯酶(pectin methylesterase,PME)基因AT5G64640是miR156潜在的靶基因,因为它在ago1-27突变体中上调1.6倍,并且受4个降解组数据(成叶、幼叶、野生型花序和xrn4花序)的支持,如图4(c). AT5G64640被划分为A组PME,其在长角果中高度或独特地表达; 该结果说明,miR156除调节SQUAMOSA PROMOTER(SPL)外,还具有调控其他基因的新功能.
本研究对ago1-27突变体进行了RNA-seq,发现许多已确认的miRNA靶基因在ago1-27中表达上调.以往的研究对少数几个miRNA靶基因进行了分析,在RNA水平它们在ago1-27中表达上调; 在蛋白水平它们在ago1-27中上调更为明显[13].因此,ago1-27点突变可能主要影响了miRNA介导的翻译抑制过程.文献[11]的研究表明,在ago1-27中miRNA与内质网膜的关联度下降.由于内质网膜是蛋白翻译的重要场所[14-16],miRNA与膜的解离可能导致miRNA介导的翻译抑制过程的失调.
目前已确认的miRNA靶基因集中在保守miRNA和紧密互补的类型上[17],新的证据表明松弛配对类型的靶基因确实存在[18-20].本研究结合miRNA靶基因预测、降解组数据和ago1突变体RNA-seq分析,鉴定出新的miRNA靶基因,这有助于深入认识植物miRNA和靶基因的调控网络.
深圳大学学报理工版
JOURNAL OF SHENZHEN UNIVERSITY SCIENCE AND ENGINEERING
(1984年创刊 双月刊)
主 管 深圳大学
主 办 深圳大学
编辑出版 深圳大学学报理工版编辑部
主 编 阮双琛
国内发行 深圳市邮电局
国外发行 中国国际图书贸易集团有限公司(北京399信箱)
地 址 北京东黄城根北街16号
邮 编 100717
电 话 0755-26732266
0755-26538306
Email journal@szu.edu.cn
标准刊号 ISSN 1000-2618
CN 44-1401/N