基于转录组测序数据识别黑猩猩RNA编辑位点

使用转录组测序(RNA-Seq)数据识别黑猩猩RNA编辑位点,探索了RNA编辑的识别机制以及潜在的功能影响.基于黑猩猩RNA-Seq数据与基因组序列的比对信息发现RNA-DNA错配位点,并构建编辑位点候选集.从中滤除基因组或转录组测序质量低的位点,其他的过滤条件包括3′端测不准、覆盖度、SNP位点以及估算的编辑水平.构建二项分布统计模型和Bonferroni多重检验滤除候选集中的随机错误,得到RNA编辑位点.选取落在已知基因上的编辑位点进行功能分析,并用Two Sample Logo软件分析编辑位点上下游序列的特征.识别出黑猩猩12种碱基替换型RNA编辑位点8 334个,其中有41个编辑位点改变原有的氨基酸,另有3个编辑位点落在microRNA(miRNA)潜在靶基因的种子结合区.统计学分析表明,分别有640和872个RNA编辑位点存在组织和性别差异.上下游碱基频率分析表明,多种类型的编辑位点紧邻碱基具有显著偏好.结果显示,RNA编辑在黑猩猩体内大量存在,且潜在具有重要的生物学功能,为进一步深入研究灵长类RNA编辑的机制奠定了基础.

基于RNA-Seq数据识别果蝇剪接位点和可变剪接事件

完整基因结构的预测是当前生命科学研究的一个重要基础课题,其中一个关键环节是剪接位点和各种可变剪接事件的精确识别.基于转录组测序(RNA-seq)数据,识别剪接位点和可变剪接事件是近几年随着新一代测序技术发展起来的新技术策略和方法.本工作基于黑腹果蝇睾丸RNA-seq数据,使用TopHat软件成功识别出39718个果蝇剪接位点,其中有10584个新剪接位点.同时,基于剪接位点的不同组合,针对各类型可变剪接特征开发出计算识别算法,成功识别了8477个可变剪接事件(其中新识别的可变剪接事件3922个),包括可变供体位点、可变受体位点、内含子保留和外显子缺失4种类型.RT-PCR实验验证了2个果蝇基因上新识别的可变剪接事件,发现了全新的剪接异构体.进一步表明,RNA-seq数据可有效应用于识别剪接位点和可变剪接事件,为深入揭示剪接机制及可变剪接生物学功能提供新思路和新手段.