文献名: DIAlignR provides precise retention time alignment across distant runs in DIA and targeted proteomicsDIAlignR 为DIA和靶向蛋白质组学提供了准确的保留时间对齐的算法)

期刊名:MCP

发表时间:(2019年1月)

单位:

  1. 多伦多大学
  2. 斯坦福大学

 

一、 概述:(用精炼的语言描述文章的整体思路及结果)

我们提出了一种新的保留时间算法,直接比对原始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的更多相关文章

  1. InfluxDB学习之InfluxDB数据保留策略(Retention Policies)

    InfluxDB每秒可以处理成千上万条数据,要将这些数据全部保存下来会占用大量的存储空间,有时我们可能并不需要将所有历史数据进行存储,因此,InfluxDB推出了数据保留策略(Retention Po ...

  2. Alignment trap 解决方法  【转 结合上一篇

    前几天交叉编译crtmpserver到arm9下.编译通过,但是运行的时候,总是提示Alignment trap,但是并不影响程序的运行.这依然很令人不爽,因为不知道是什么原因引起的,这就像一颗定时炸 ...

  3. ARMLinux下Alignment trap的一些测试 【转自 李迟的专栏 CSDN http://blog.csdn.net/subfate/article/details/7847356

    项目中有时会遇到字节对齐的问题,英文为“Alignment trap”,如果直译,意思为“对齐陷阱”,不过这个说法不太好理解,还是直接用英文来表达. ARM平台下一般是4字节对齐,可以参考文后的给出的 ...

  4. Multiple sequence alignment Benchmark Data set

    Multiple sequence alignment Benchmark Data set 1. 汇总: 序列比对标准数据集: http://www.drive5.com/bench/ This i ...

  5. Retention、Documented、Inherited三种注解

    Retention注解 Retention(保留)注解说明,这种类型的注解会被保留到那个阶段. 有三个值:1.RetentionPolicy.SOURCE —— 这种类型的Annotations只在源 ...

  6. POJ 1836 Alignment

    Alignment Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 11450 Accepted: 3647 Descriptio ...

  7. cf.295.C.DNA Alignment(数学推导)

    DNA Alignment time limit per test 2 seconds memory limit per test 256 megabytes input standard input ...

  8. annotation(@Retention@Target)详解

    一.注解:深入理解JAVA注解 要深入学习注解,我们就必须能定义自己的注解,并使用注解,在定义自己的注解之前,我们就必须要了解Java为我们提供的元注解和相关定义注解的语法. 1.元注解(meta-a ...

  9. Rman备份的保留策略(retention policy)

    什么是备份的保留策略(retention policy) 保留策略就是指某份备份是否要保留以及保留多久.可以使用configure retention policy命令来创建一个一致的.自动的备份保留 ...

随机推荐

  1. c语言代码风格

    简要:主要介绍了K&R风格和Allman(BSD)风格之间的缩进大小和大括号位置等区别 关于其它的代码风格,详见:Indent style - Wikipedia. 1.K&R sty ...

  2. Java企业微信开发_01_接收消息服务器配置

    一.准备阶段 需要准备事项: 1.一个能在公网上访问的项目: 见:Java微信公众平台开发_01_本地服务器映射外网 2.一个企业微信账号: 去注册:(https://work.weixin.qq.c ...

  3. KVM-快照管理

    1.kvm克隆 kvm 虚拟机有两部分组成:img镜像文件和xml配置文件(/etc/libvirt/qemu 克隆命令:virt-clone -o rhel6- 71 -n xuegod63-kvm ...

  4. [BZOJ2962][清华集训]序列操作

    bzoj luogu 题意 有一个长度为\(n\) 的序列,有三个操作: \(I \ \ a\ b\ c\ :\)表示将\([a,b]\)这一段区间的元素集体增加\(c\): \(R \ \ a\ b ...

  5. C#添加修改控件css样式

    一.添加属性 MyStyleSheet.Attributes.Add("href","/css/flostyle.css") 二.改变css样式 if (use ...

  6. java流类

    总结:new FileInputStream package com.ds; import java.io.*; import com.da.fgbv; public class rter { pub ...

  7. Linux系统主流架构一

    Linux系统主流架构一 随着IT运维的不断发展,尤其的Linux的飞速发展,越来越多的企业开始使用Linux操作系统,例如:Centos.Red Hat.Ubuntu.Fadora等等,成千上亿个网 ...

  8. DataTable列查询加排序

    DataTable列查询加排序 DataRow[] drArray = dt.Select("ANLYCOM_ID='" + chSPrdtStblAnly.AnlyComId + ...

  9. 基于OpenCV之视频读取,处理和显示框架的搭建(一)

    主要包括以下内容: 1.使用的主要函数的说明. 2.两个实例:视频读取和显示.搭建视频读取和处理框架,调用canny函数提取边缘并显示. 3.一些注意事项和代码说明. 一.使用的主要函数 1.延时函数 ...

  10. 关于SCANF接受的一些总结

    当C语言中SCANF函数接受到非指定类型的参数时,如:定义为整形输入为字符串,函数会返回一个失败的BOOL类型,以此可以判断是否输入的为指定类型. 当发生一次输入错误的时候,如果想再次进行输入则需要使 ...