MCP|LQ|DIAlignR provides precise retention time alignment across distant runs in DIA and targeted proteomics
文献名: DIAlignR provides precise retention time alignment across distant runs in DIA and targeted proteomics(DIAlignR 为DIA和靶向蛋白质组学提供了准确的保留时间对齐的算法)
期刊名:MCP
发表时间:(2019年1月)
单位:
- 多伦多大学
- 斯坦福大学
一、 概述:(用精炼的语言描述文章的整体思路及结果)
我们提出了一种新的保留时间算法,直接比对原始swathm – ms的二级色谱图的数据,采用了混合动态规划方法。该算法不强制按洗脱时间顺序排列色谱峰洗脱,峰对齐过程中允许峰位置的交换。粗略的全局拟合使其具有抗噪声的稳定性。在一个手动验证数据集,该策略优于当前最先进的方法。局部对其使得其精度更高。
二、 研究背景:(简要介绍研究进展动态、研究目的和意义)
尽管现在已经有大量的建立保留时间一致的方法。然而,目前在代谢组学和蛋白质组学中,保留时间对齐的算法是在SWATH-MS方法出现之前设计的。
这些算法通常使用动态规划算法在原始的MS1色谱图上或者特征表中寻找一个全局配对的对齐函数。然而大部分方法依赖于MS1数据,并且以此生成的对齐函数受到所含肽段的影响。 这些方法要么使用双向匹配,或者利用特征通过LOESS或者核密度函数来计算全局函数,从而对齐二级质谱,这在高噪声、缺少特征或者特征检测算法出错的情况下提供了一个次优解。全局单调函数也是如此,由于是单调函数,不考虑任何两个肽段之间保留时间的翻转。
本文提出了DIAlignR保留时间对齐算法, 该算法可以克服其他方法的缺点。不需要特征,就能直接将来自靶向蛋白质组学的原始多重MS2色谱峰图对齐。本方法使用动态规划来获得两者之间的最佳映射,含有局部信息的色谱图,如多个邻近峰进行对齐。该方法还能够使用全局运行对齐指导,使其抗噪声能力强。DIAlignR可以灵活处理用户的偏好,用户在全局对齐和局部对齐的两个极端之间进行选择。
三、实验设计:

四、研究成果:(重点图表展示)
靶向蛋白组学MS2色谱对齐算法。
a:多肽离子碎片的色谱图,两次测得的数据,run A在上面,run B在下面。
b:两次运行的色谱图相似度用矢量的点积计算
c:外点积色谱图相似度评分矩阵(S)用作对齐的近似路径。
d:基于特征的对齐,通过相似矩阵计算出最佳得分路径,偏离一定范围的点使用负分数来惩罚。
e:仿射罚分策略通过相似矩阵计算最高分的路径,该动态规划方法使用三维矩阵递归打分。对齐的路线使用黑色线表示
f:将色谱峰强度映射到对齐的时间路径上

使用全局优先的方法对不同相似性测量的方法,技术参数,相似性罚分对S. Pyogenes 数据集对齐准确性进行比较a:不同相似性测量方法的比较。b间隙罚分的选取对对齐结果在一定的RT差异容忍度内的峰百分比的影响。c, DGSVSVADSGR/2肽段在run11和run12的罚分相似性矩阵。d, 提取的肽段XIC图的端点。

阅读人:刘权
MCP|LQ|DIAlignR provides precise retention time alignment across distant runs in DIA and targeted proteomics的更多相关文章
- InfluxDB学习之InfluxDB数据保留策略(Retention Policies)
InfluxDB每秒可以处理成千上万条数据,要将这些数据全部保存下来会占用大量的存储空间,有时我们可能并不需要将所有历史数据进行存储,因此,InfluxDB推出了数据保留策略(Retention Po ...
- Alignment trap 解决方法 【转 结合上一篇
前几天交叉编译crtmpserver到arm9下.编译通过,但是运行的时候,总是提示Alignment trap,但是并不影响程序的运行.这依然很令人不爽,因为不知道是什么原因引起的,这就像一颗定时炸 ...
- ARMLinux下Alignment trap的一些测试 【转自 李迟的专栏 CSDN http://blog.csdn.net/subfate/article/details/7847356
项目中有时会遇到字节对齐的问题,英文为“Alignment trap”,如果直译,意思为“对齐陷阱”,不过这个说法不太好理解,还是直接用英文来表达. ARM平台下一般是4字节对齐,可以参考文后的给出的 ...
- Multiple sequence alignment Benchmark Data set
Multiple sequence alignment Benchmark Data set 1. 汇总: 序列比对标准数据集: http://www.drive5.com/bench/ This i ...
- Retention、Documented、Inherited三种注解
Retention注解 Retention(保留)注解说明,这种类型的注解会被保留到那个阶段. 有三个值:1.RetentionPolicy.SOURCE —— 这种类型的Annotations只在源 ...
- POJ 1836 Alignment
Alignment Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 11450 Accepted: 3647 Descriptio ...
- cf.295.C.DNA Alignment(数学推导)
DNA Alignment time limit per test 2 seconds memory limit per test 256 megabytes input standard input ...
- annotation(@Retention@Target)详解
一.注解:深入理解JAVA注解 要深入学习注解,我们就必须能定义自己的注解,并使用注解,在定义自己的注解之前,我们就必须要了解Java为我们提供的元注解和相关定义注解的语法. 1.元注解(meta-a ...
- Rman备份的保留策略(retention policy)
什么是备份的保留策略(retention policy) 保留策略就是指某份备份是否要保留以及保留多久.可以使用configure retention policy命令来创建一个一致的.自动的备份保留 ...
随机推荐
- php获取客户端IP地址的几种方法(转)
[php] view plain copy php获取客户端IP地址的几种方法 方法一 <?php $iipp=$_SERVER["REMOTE_ADDR"]; echo $ ...
- codeforces 553A A. Kyoya and Colored Balls(组合数学+dp)
题目链接: A. Kyoya and Colored Balls time limit per test 2 seconds memory limit per test 256 megabytes i ...
- uoj problem 10
uoj problem 10 题目大意: 给定任务若干,每个任务在\(t_i\)收到,需要\(s_i\)秒去完成,优先级为\(p_i\) 你采用如下策略: 每一秒开始时,先收到所有在该秒出现的任务,然 ...
- [CERC 2008] Suffix reconstruction
[题目链接] https://www.lydsy.com/JudgeOnline/problem.php?id=4319 [算法] 首先 , 我们可以求出这个字符串的rank数组 按照SA逐位枚举 , ...
- UILabel的富文本显示选项
UILabel的富文本格式设置 1.实例化方法和使用方法 实例化方法: 使用字符串初始化 - (id)initWithString:(NSString *)str; 例: NSMutableAttri ...
- github 在线 创建文件夹
参考: http://webapps.stackexchange.com/questions/36411/create-a-folder-in-github-via-the-web-interface
- BZOJ4695:最假女选手
浅谈区间最值操作和历史最值问题:https://www.cnblogs.com/AKMer/p/10225100.html 题目传送门:https://lydsy.com/JudgeOnline/pr ...
- 洛谷 P4546 & bzoj 5020 在美妙的数学王国中畅游 —— LCT+泰勒展开
题目:https://www.luogu.org/problemnew/show/P4546 先写了个55分的部分分,直接用LCT维护即可,在洛谷上拿了60分: 注意各处 pushup,而且 spla ...
- Servlet的生命周期以及简单工作原理的讲解
Servlet生命周期分为三个阶段: 1,初始化阶段 调用init()方法 2,响应客户请求阶段 调用service()方法 3,终止阶段 调用destr ...
- shell入门-系统和用户的配置文件
变量分为系统变量和用户自定义变量,他们的配置文件当然分为系统环境变量配置文件和用户环境变量控制文件 1.规定环境变量的配置文件 [root@wangshaojun 111]# vim /etc/pro ...