RNA m6A位点预测
前面已经给大家介绍过m6A甲基化的一些基本概念以及m6A的检测方法
也给大家分享过一些m6A相关的发文套路
当然我们也会分享一些m6A数据分析相关的干货
☞m6a甲基化相关基因根据临床信息分组绘制boxplot并显示p值
在☞RNA m6A检测方法一文中我们介绍了MeRIP-qPCR这种方法。在第四步中我们提到IP下来的RNA需要用随机引物进行逆转录,而第四步中设计的引物用的是特异性引物,这个大家要注意区分和甄别。有人可能会比较好奇,这个特异性引物是怎么设计的。首先我们需要确定m6A究竟发生在mRNA上的什么位置,接下来才能确定究竟应该在哪里放特异性引物。
下图是一个m6A位点在mRNA上的分布图,可以看出m6A修饰倾向于富集在终止密码子附近。当然在5‘UTR和编码区也有。
以上只是一个大规模的统计结果,那么针对于不同的mRNA,m6A位点究竟存在于什么位置呢?今天就给大家介绍一个免费在线预测哺乳动物m6A修饰位点的网站SRAMP。
网址:http://www.cuilab.cn/sramp
1. 首先进入SRAMP网页
点击“Prediction”按钮
2. 选择预测的模型
左边的“Full transcript mode”在对编码和非编码RNA进行预测时,建议使用此模式。 注意,在该模式中,应使用完整转录物(具有内含子)而不是成熟mRNA / cDNA序列的基因组序列。而对于没有完整基因组序列的用户来说,可以选择右侧的“Mature mRNA mode”,这种预测模式是一种备选解决方案,相对来说预测的准确性不如左边的,该模型适用于成熟mRNA(cDNA)序列,不能预测内含子中的m6A位点。
我们用网站默认的test1(TROVE2)序列演示一下,后面的选项都是默认的,点击“Submit”就可以预测序列的m6A结合位点了。
3. 预测结果
有两个结合位点,在“Decision”列中可以看出预测的两个结合位点具有很高的可信度。SRAMP网站提供了几个阈值:99 %/95%/90%/85%,分别对应very high/high/moderate/low的自信度,当然我们一般都挑选“very high”啦。预测结果是按照“Position”位置排序的,所以还需要综合各个位点“Score”情况。
SRAMP网站还提供了RNA二级结构m6A位点结合预测功能,写材料,发文章,有个图不是更加美观?针对于这个功能,网站默认选项是“NO(default)”,我们选择“YES(much slower)”,然后点击“Submit”就行,计算过程需要花更多的时间,这个取决于提交序列的长度。
预测结果
点击“Draw”,RNA二级结构m6A位点预测图就出来了
SRAMP网站运算速度比较慢,如果需要预测的序列较多,有一台不错的电脑,可以自行下载SRAMP tool(压缩包大概2G),同时安装R语言和Perl语言就可以运算了。
SRAMP算是一个不错的网站,界面相当简洁,由北京大学创建的;下面两个m6A甲基化预测的网站都是由中山大学创建的。
(1)RMBase v2.0
网址:http://rna.sysu.edu.cn/rmbase/
(2)m6Avar
网址:http://m6avar.renlab.org
相关文章
- 数据挖掘项目一般多少钱_预测类数据挖掘项目
- Python中利用长短期记忆模型LSTM进行时间序列预测分析 - 预测电力负荷数据|附代码数据
- R语言Lasso回归模型变量选择和糖尿病发展预测模型|附代码数据
- 使用PyG进行图神经网络的节点分类、链路预测和异常检测
- 【视频】随机波动率SV模型原理和Python对标普SP500股票指数预测|数据分享|附代码数据
- Adv Sci | 首创引入蛋白动态信息的AI模型,高效预测药物-蛋白亲和力
- 【鹅厂网事】服务器硬盘故障预测实践
- [KDD 2022 | 论文简读] KPGT: 用于分子性质预测的知识指导的预训练图形变换模型
- [Bioinformatics | 论文简读] 用于从单细胞 RNA 测序数据中预测基因调控的图注意力网络
- 快手技术VP王仲远对于AI、大模型、深度学习的预测
- R语言预测人口死亡率:用李·卡特(Lee-Carter)模型、非线性模型进行平滑估计|附代码数据
- MATLAB用GARCH-EVT-Copula极值理论模型VaR预测分析股票投资组合|附代码数据
- R语言随机森林RandomForest、逻辑回归Logisitc预测心脏病数据和可视化分析|附代码数据
- Gartner预测:2021年全球5G网络基础设施收入将增长39%
- 凯文·凯利来华再谈人工智能,并预测未来 25 年的技术发展趋势