基于生物质谱数据鉴定单核苷酸变异的生物信息学方法,生物化学论文.docx

《基于生物质谱数据鉴定单核苷酸变异的生物信息学方法,生物化学论文.docx》由会员分享，可在线阅读，更多相关《基于生物质谱数据鉴定单核苷酸变异的生物信息学方法,生物化学论文.docx（25页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、基于生物质谱数据鉴定单核苷酸变异的生物信息学方法,生物化学论文单核苷酸变异(single nucleotide variations, SNVs)是由 DNA 序列上单个碱基变异产生的, 包括碱基的缺失、插入、转换及颠换等. SNVs 是基因组序列变异的主要形式1, 同时也是生物体生理和病理变异的遗传基础2. 从遗传学的角度看, SNVs 既能够存在于具有遗传性的生殖细胞中, 可以以存在于不具有遗传性的体细胞中. 华而不实, 只要位于基因编码区的SNVs 能够影响蛋白的编码. 位于编码区的 SNVs 能够分为 3 类: () 同义 SNVs, 不改变相应的氨基酸种类; () 无义 SNVs,

2、突变成为终止密码子, 提早结束编码; () 非同义 SNVs(nonsynonymous SNVs,nsSNVs), 改变氨基酸的种类. nsSNVs 能够改变蛋白的构造、功能、表示出以及亚细胞定位等3, 进而对多种遗传性的特征、疾病以及癌症等产生影响49, 如人类耳垢的类型6、腋窝的气味7、癌症与肿瘤的发生8、阿尔茨海默病9以及镰刀形红细胞贫血症10等. 因而, 对 SNVs 展开研究能够揭示出基因与表型多样性和基因与疾病间的关系, 并且有可能研发出治疗疾病的新方式方法. 当前, 全基因组关联研究(genome-wide association studies, GWAS)11固然在基因变异

3、与表型多样性的研究中产出了很多能够用来解释特异性疾病分子途径的结果, 但是仍然难以对绝大部分具有复杂特征的分子机制以及SNVs与复杂疾病表型间的关系进行解释12. 在这种情况下, 对突变蛋白的研究提供了另一种了解基因型与表型间关联的方式方法13. 由 SNVs 引起的单个氨基酸的变异称为单氨基酸变异(single amino acid variations, SAVs), 因而SAVs 是 SNVs 在蛋白水平上的表现. 对 SAVs 的研究, 有助于了解基因型与表型间的关系, 进而从本质上了解基因是如何在蛋白水平上影响生物体的生命经过的14. 当前, 基于串联质谱的鸟枪法蛋白质组学(shot

4、gun proteomics)技术由于其自动化、高通量、高灵敏度和高分辨率等特点, 已成为大规模蛋白质研究的主要方式方法. 序列数据库搜索算法由于具有较高的可靠性以及灵敏度而成为当今鸟枪法蛋白质组学中蛋白鉴定的主要生物信息学方式方法. 然而, 通常蛋白质数据库在构建时为了减小数据库的冗余程度, 往往有意压缩对 SAVs 信息的收录(如 Swiss-Prot 数据库15,16, IPI 数据库17等), 进而使得常用的数据库搜索策略不能有效地鉴定出样本中的氨基酸突变信息. 为此, 研究人员提出了一系列鉴定突变蛋白的方式方法,如构建包含有突变信息的蛋白质数据库、构建类似性图谱库等.在基于串联质谱进

5、行 SAVs 鉴定时, 能够采用与蛋白质翻译后修饰(post-translational modifications,PTMs)鉴定18一样的方式方法, 这是由于肽段的突变和修饰在质谱图中的表现都是质量迁移, 如甲硫氨酸(Met)氧化与丙氨酸(Ala)突变为丝氨酸(Ser)在质量上都是增加 16 Da19, 所以鉴定 PTMs 的算法和流程通常也能够鉴定 SAVs(如 Bonanza 算法20). 固然 PTMs和 SAVs 的质谱鉴定方式方法非常类似, 但由于其来源上的差异不同, 在实际的鉴定策略中有所不同. () PTMs的种类远比 SAVs 要多, 鉴定 PTMs 所需的搜索空间一般会比

6、鉴定 SAVs 所需的大, 在质量控制方面具有更大的挑战; () 蛋白水平的 SAVs 大部分是从基因组或转录组延续过来的, 充分利用 SNVs 的数据能大大降低搜索空间, 进而得到更可靠的结果. 因而在计算方式方法与策略方面, SAVs 和 PTMs 的鉴定具有一定的类似性, 也有其独有的特点. 本文从序列数据库搜索算法、序列标签搜索算法以及图谱库搜索算法 3 个大方面, 具体地介绍了当前基于生物质谱数据鉴定SAVs的各种生物信息学方式方法,并分析了各种突变鉴定方式方法的缺乏之处, 最后介绍了基于生物质谱的 SAVs 鉴定研究现在状况及其发展方向. 1 氨基酸突变鉴定的算法 当下基于生物质谱

7、的 SAVs 鉴定算法都是由常规鉴定算法改良而来的, 因而根据常规串联质谱鉴定算法中对数据库的依靠程度以及使用的数据库种类, 能够将基于生物质谱的 SAVs 鉴定算法分为 3 大类(表 1): () 完全依靠序列数据库的搜索算法, 即基于序列数据库搜索的氨基酸突变鉴定算法. 此算法利用前体离子质量从序列数据库中挑选出候选肽段, 然后将候选肽段的理论图谱与目的图谱进行比对, 进而鉴定出样品中的突变肽段; () 将从头测序算法(de novo)与序列比对结合的算法, 即基于序列标签的氨基酸突变鉴定算法. 此算法首先通过 de novo测序算法推导出目的图谱中的肽序列标签(peptidesequen

8、ce tags, PSTs), 然后利用 PSTs 过滤数据库挑选出候选肽段, 最后结合 PSTs 对理论谱图与目的图谱进行比拟打分, 进而鉴定出样品中的突变肽段; ()依靠于图谱库的搜索算法, 即基于图谱库的氨基酸突变鉴定算法. 此算法将实验图谱与图谱库中的一致性图谱进行比对, 进而鉴定出样品中的突变肽段.这 3 类方式方法和策略在施行经过中各有其优劣(表 1),互相之间暂无法替代, 因而在不同的目的下各有其适用性. 1.1 基于序列数据库搜索的氨基酸突变鉴定算法 基于序列数据库搜索的氨基酸突变鉴定算法,根据不同的数据库构建方式方法能够细分为 3 类: () 基于穷举法的氨基酸突变鉴定算法,

9、 即通过枚举数据库中氨基酸残基的所有可能突变种类进行突变肽段的鉴定; () 结合已经知道氨基酸突变信息对突变肽段进行鉴定, 即结合当下变异数据库(如 dbSNP 数据库21、COSMIC 数据库22等, 表 2 列举了常用的氨基酸与基因突变数据库)中的变异信息构建数据库进行突变肽段的鉴定; () 基于样本特异性的数据库鉴定突变肽段, 即结合样本数据中可能存在的突变肽段信息构建数据库进行突变肽段的鉴定. 下面将对这 3 种方式进行逐一具体地讲明. (1) 基于穷举法的氨基酸突变鉴定算法. 在序列数据库搜索中, 最早对突变肽段进行鉴定的自动化方式方法是穷举法, 此方式方法不仅原理简单而且理论上能够

10、鉴定出样品中所有可能的突变肽段. 这类算法的大体步骤是: 通过穷举法罗列出所有可能的突变肽段序列, 然后用常规鉴定方式方法进行比对打分挑选出最有可能的突变肽段序列. 此类算法的代表有SEQUEST-SNP 算法27和 Sipros v2.0 算法18等.Gatlin 等人27在2000年, 利用改良的 SEQUEST算法(SEQUEST-SNP)率先实现了利用自动化的数据库搜索对突变肽段进行鉴定. 此方式方法特点在于动态生成所有可能的核苷酸突变序列, 将其翻译成肽段并构建成一个数据库用于对突变肽段的鉴定. 此后, 通过穷举蛋白序列中所有可能的氨基酸突变进行肽段突变鉴定的方式方法在 Mascot

11、28和 X!Tandem29相继采用.2020 年, Hyatt 和 Pan18提出了不受数据库约束的穷举法突变肽段鉴定算法 Sipros v2.0, 此算法通过肽段产生模块和肽段打分模块实现对 CPU 和内存效率的优化以应对穷举法产生的大数据库. 理论上, 穷举法能够鉴定出样品中所有的突变肽段, 但肽段中的每一个氨基酸残基都有 18 种可能的突变, 因而利用此方式方法会大大增加搜索空间18,24, 延长搜索时间, 并且会增加假阳性风险进而降低结果的灵敏度. (2) 结合已经知道氨基酸突变信息对突变氨基酸进行鉴定. 为了避免穷举法引起搜索空间过大的问题,一些团队提出结合已经知道的编码SNVs信

12、息或是与疾病等有关的突变信息构建蛋白质数据库, 以减小突变肽段的搜索范围. 此类数据库的代表有 MSIPI17和MS-CanProVar24等. 2007 年, Schandorff 等人17将一些来自 dbSNP 数据库21的编码 SNP(single nucleotidepolymorphism) 以 及 与 IPI(the international proteinindex)数据库中数据有冲突的序列等整合到 IPI 数据库30中构建了质谱友好型的变异数据库 MSIPI. 其质谱友好型具体表现出在, 在保存原始 IPI 条目完好性的基础上, 将后加的肽段序列附加到原有序列中, 用不代表任

13、何氨基酸的字母 J 将原始条目与附加肽段区分开来, 并且将在原始条目的表头信息中参加附加肽段信息. 同年, Bunger等人31也利用dbSNP数据库中人类基因变异信息构建变异蛋白质数据库 K-SNPdb,并构建相应的常规数据库. 然后对分开搜库结果进行比对打分, 挑选出高可信的变异肽段. Li 等人24在2018 年基于人类癌症蛋白质变异数据库 CanProVar32构建了一个 MS-CanProVar 数据库, 此数据库中不仅包含了 dbSNP 数据库中的编码的 SNP 信息, 还包括了 COSMIC22和 OMIM23等数据库中与癌症相关的体细胞变异信息. 除了自定义构建突变数据库以外,

14、 氨基酸突变信息也被一些在线平台收录、整合, 如 Swiss-Var33,SysPIMP34和 RAId_DbS35等. Swiss-Var 网站搜集的是 Swiss-Prot 数据库36中突变肽段的信息, 主要为用户提供Swiss-Prot数据库中的突变肽段信息及其与疾病间的关系. SysPIMP 主要用于鉴定与人类疾病有关的突变肽段序列, 它的数据主要来源于 OMIM 数据库中等位基因突变信息、蛋白质突变数据库(proteinmutation database, PMD)37以及 Swiss-Prot 数据库中与 人 类 疾 病 和 多 态 性 有 关 的 序 列 信 息 . 而 在RAI

15、d_DbS 数据库中不仅整合了 SAVs 与疾病的信息,同时也收录了 PTMs 与疾病有关的信息. 2020 年, Mathivanan 等人25提出的 iMASp 策略即是利用现有的突变信息对突变肽段进行鉴定. 这种策略利用了分步搜索的方式方法, 即是第一次通过常规搜索鉴定出样本中的常规蛋白, 第二次利用突变数据库对第一次没有鉴定出的质谱图进行搜索鉴定样品中的突变肽段. 相比穷举法, 结合已经知道氨基酸突变信息对突变氨基酸进行鉴定的方式方法固然在一定程度上缩小了搜索空间, 但在数据库中添加的上万条突变肽段序列绝大部分不会在样品数据集中出现.因而, 这种方式方法并没有特别有效地躲避假阳性升高以

16、及鉴定结果灵敏性降低的缺点14. (3) 基于样本特异性的数据库鉴定突变肽段 .除了直接利用公共数据库中的突变数据外, 利用DNA/RNA 等信息提供的样本特异性突变构建的数据库能更好地贴合实际样本数据, 提高鉴定效率. 当前利用样本特异性鉴定突变肽段的方式方法有2种: 两次搜索数据库的方式方法以及利用转录组数据构建数据库的方式方法. 两次搜索数据库的方式方法与 iMASp 策略中所使用的分步搜索以及 Mascot 和 X!Tandem 中的容错搜索类似, 不同的地方在于两次搜索数据库中所使用的突变数据库依靠于样本特异性的 DAN/RAN 信息, 而 iMASp 策略中的突变数据库是整合所有已

17、经知道的蛋白突变信息, 不具有样本特异性; Mascot 和X!Tandem 则是对第一次搜索所得的蛋白序列进行穷举进而鉴定出突变或修饰肽段. Chernobrovkin 等人38提出的二次迭代法以及 Su 等人39构建样本特异性突变数据库的策略都是样本特异性的两次搜索方式方法的代表. 另一种方式方法是利用转录组数据构建样本特异性数据库用于突变肽段的鉴定. 相对于利用公共的突变数据库, 利用转录组数据构建蛋白质数据库能够由样品转录组数据直接推导样本中可能存在的蛋白及其突变序列并由其构建数据库40. 用此方式方法构建的数据库所包含的蛋白质信息愈加接近样品中真实信息, 因而这种无偏性的数据库能高效

18、地鉴定出样品中存在的突变序列16,41. 由于转录组数据特别庞大, 在现有的计算能力下要想利用转录组数据构建数据库就必需要对转录组数据进行压缩. 2007 年,Edwards16提出了一个压缩表示出序列标签(expressedsequence tags, ESTs)数据的策略, 实现了利用 EST 数据库进行常规化的肽段序列和变异位点的鉴定. 此压缩策略的特点在于选用某种方式方法来表示肽段, 确保绝大多数的重复肽段序列被消除, 并且不影响肽段序列的鉴定. 随着下一代测序(next generationsequencing, NGS) 技术的出现 , RNA 测序 (RNA-sequecing,

19、 RNA-Seq)的成本越来越低14, 并且克制了 EST 测序存在的克隆偏性和高花费等缺点42, 因而利用 RNA-Seq 数据构建样本特异性数据库逐步遭到人们的重视. Wang 等人41在 2020 年提出了一个利用 RNA-Seq 数据构建样本特异性数据库的策略, 此策略通过两步来实现: () 利用一个经历体验性的RPKM (reads per kilo bases per million reads)值排除不表示出或低表示出基因以减小数据库中的条目; ()将由 RNA-Seq 数据鉴定得来的高可靠性 SNVs 的相应肽段添加到数据库中, 以寻找变异肽段. 此后,Wang 和 Zhang

20、43为生成自定义RNA-Seq数据库编写了 R程序包 customProDB, 能够生成含有突变、插入、缺失等变异肽段的 RNA-Seq 数据库. 2020 年,Sheynkman 等人14实践了 Wang 和 Zhang43的方式方法,利用 Jurkat 细胞系的 RNA-Seq 数据构建一个自定义的变异蛋白质数据库, 并成功地应用在 Jurkat 细胞系的质谱数据突变鉴定中. 同年, Woo 等人44在尽量不影响鉴定结果灵敏性的基础上, 将秀丽隐杆线虫(Caenorhabditis elegans)的 RNA-Seq 数据压缩了近1000 倍, 并利用此数据库成功地鉴定到了新型蛋白. 由于

21、并不是所有的样本都同时拥有蛋白质数据和RNA-Seq 数据, 因而, Wang 和 Zhang43利用 64 个大肠癌的 RNA-Seq 数据构建了一致性蛋白质数据库,并成功地将此数据库应用在蛋白鉴定中. 样本特异性的数据库, 十分是利用 RNA-Seq 数据构建的样本数据库不仅能够有效地缩减搜索空间, 而且能够鉴定出样品中所有已经知道类型的蛋白种类以及新型的变异肽段序列. 随着计算方式方法的不断改良, 通过RNA-Seq 数据对样本进行突变肽段的鉴定方式方法有望成为常规的突变鉴定方式方法. (4) 基于序列数据库搜索的氨基酸突变鉴定算法的缺点. 在鉴定突变肽段的方式方法中, 固然通过构建含有

22、突变信息的序列数据库鉴定突变肽段的方式方法是当前被最广泛采用的方式方法, 但它的缺点也是不容忽视的. () 除了利用穷举法构建的突变数据库以外, 利用其他方式方法构建的突变数据库对突变信息包含得都不够全面, 如公共数据库通常会有意忽略对变异数据的收录, 而样本特异性数据库为了减小搜索空间通常也会去除低表示出的蛋白质; () 序列数据搜索中, 当图谱中的碎裂信息不够完好、信噪比拟低时, 搜索引擎就不能将候选肽段正确地区分开45,因此会增加假阳性的概率. 为了避免序列数据库的上述缺点, 提出了其他鉴定突变肽段的方式方法, 如序列标签算法、图谱库搜索算法等. 1.2 基于序列标签的氨基酸突变鉴定算法

23、 相比序列数据库搜索算法利用肽段母离子质量从数据库中挑选候选肽段, 序列标签算法利用 denovo 测序算法推导的 PSTs 能够更有效地过滤数据库,减少候选肽段的数目以缩小搜索空间, 使得更复杂和计算更密集的方式方法能够应用到对候选肽段的突变打分算法中45, 进而提高了突变鉴定结果的灵敏性并且减少了结果中的假阳性率. 下面从序列标签搜索算法与 de novo 测序算法之间的关系以及当下结合PSTs 进行氨基酸突变鉴定的主流工具两个方面对序列标签算法鉴定突变氨基酸进行介绍. (1) 序列标签搜索算法与 de novo 测序算法. 相比序列数据库搜索算法, de novo 算法在对质谱图进行氨基

24、酸序列推导时不依靠蛋白质数据库, 因而它在鉴定氨基酸突变方面有独特的优势4547. 当下使用 de novo 测序算法的代表性工具有 SHERENGA48,PEAKS4951以及 PepNovo52等. 这些工具所使用的算法都是通过生成前缀残基质量图谱(prefix residuemass spectra)重构整个图谱进行肽段序列推导的, 因而这些算法对质谱图的质量具有较高的要求45. 但通过诱导碰撞解离(collision-induced dissociation, CID)产生的串联图谱中不可避免地含有不完好的碎裂离子系列、噪音离子和精度较差的碎裂离子质量, 这使得 de novo 算法经

25、常产生一些不确定的序列区域, 导致 de novo 算法通常只能准确地推导出肽段序列中的部分序列46. 因而, 结合 de novo 算法鉴定的部分肽段序列进行数据库搜索的序列标签算法应运而生,这种算法不仅能够利用 de novo 推导出的 PSTs 作为挑选候选肽段时的过滤指标, 有效地减少搜索空间,而且能够通过改变PSTs与候选肽段匹配的打分算法,提高对突变肽段的鉴定效率. (2) 结合肽序列标签的氨基酸突变鉴定算法 .最早结合 PSTs 进行数据库搜索的方式方法是由 Mann 和Wilm53在 1994 年提出的, 此方式方法不仅能有效地对常规图谱进行鉴定, 而且能够鉴定出带有突变或修饰

26、图谱的肽段序列. 当下结合肽序列标签对氨基酸突变进行鉴定的算法或程序有 GutenTag 程序54,Opensea 工具55, SPIDER 程序56,57, InsPecT 搜索引擎45, DirecTag 算法58以及 MoDa 算法59等. 鉴定突变氨基酸常用的序列标签软件及其网址见表 3. GutenTag 是由 Yates 实验室开发出来的能够自动推导+2电荷母离子串联图谱PSTs用于数据库搜索的算法, 其特点是利用碎片离子峰强度经历体验模型并结合相邻氨基酸和碎片离子的相对质量对肽段碎裂的影响推导 PSTs, 之后用多个 PSTs 进行搜库, 同时放宽对 PST 两端质量匹配的限制,

27、 进而能够有效地进行突变肽段的鉴定. 但由于GutenTag算法没有考虑同源突变或修饰, 所以此算法只能对数据库中已存在的突变序列进行鉴定, 并且由于在打分方面存在漏洞55,所以鉴定出来的结果中存在较高的假阳性. 在 GutenTag 算法发表后的第 2 年, Searle 等人55初次将序列标签算法的思想应用于非限制翻译后修饰, 并提出了基于质量的序列比对算法工具Opensea.Opensea 的特点是利用基于质量的宽度优先的算法( breadth-first search algorithm)鉴定出突变位点或修饰位点. 但宽度优先的算法是一种贪心的匹配算法,并且在 Opensea 中没有考

28、虑在一个位点上同时存在de novo 的测序错误和同源突变的情况, 所以它不能保证最终结果的可靠性. SPIDER 方式方法与 Opensea 工具有类似的序列标签算法思想, 但与 Opensea 不同的是, 它能够在一个位置上同时考虑 de novo 的测序错误和同源突变的情况, 并且利用动态规划算法进行比对打分. SPIDER 算法已被整合进 PEAKS 软件中,专门用来对突变肽段和跨物种的同源性肽段进行鉴定. 在 GutenTag 算法推出后, Pevzner 实验室迅速推出了 InsPecT 序列标签算法搜索引擎45, 它是最早实现规模化鉴定翻译后修饰肽段的搜索工具, 如今仍然被广泛使

29、用. InsPect 搜索引擎推导 PSTs 的算法的特点在于利用改良的 de novo 算法推导出 PSTs 作为过滤器缩小候选肽段的范围, 并利用树状快速搜索方式方法(fast tree-based search)找出与 PSTs 匹配的候选肽段, 用基于动态规划算法(dynamic programming)的图谱比对方式方法鉴定修饰肽段, 并在打分算法中考虑肽段的碎裂形式. 在推导 PSTs 时, InsPecT 需要构建前缀残基质量图, 而 DirecTag 算法则是直接利用串联图谱的质核比值和峰强度信息对可能的标签进行打分. 由于 DirecTag 只能用来推导 PSTs, 因而其团

30、队后续开发了 TagRecon 算法47并将 DirecTag,TagRecon 和 IDPicker 工具60整合成鉴定突变和修饰肽段的流程, 其大致经过为: () 利用 DirecTag 生成PSTs; () TagRecon 利用 PSTs 对常规数据库进行候选肽段过滤, 并且定位数据集中的突变或修饰肽段;() 利用 IDPicker 工具对鉴定结果进行质量控制并且装配成蛋白. 此流程算法在 2020 年由 Abraham 等人19在鉴定胡杨树(Populus)单氨基酸多态性的实验中被成功地使用. 当前序列标签算法都依靠于 de novo 测序构建PSTs, 但是由 de novo 算法

31、测出的肽片段往往存在部分构建错误的序列56. MoDa 算法59在搜索候选肽段时, 由于采用序列标签链算法(tag chain algorithm)61,能有效地避免由 de novo 测序引起的错误匹配. 在MoDa 算法中, 将序列标签算法和动态规划算法结合,同时利用多条序列标签与候选肽段进行比对, 找出存在质量差的位点, 然后利用基于动态规划算法的图谱比对算法找出最佳的肽段序列. 此方式方法能够大规模地鉴定出存在多个修饰位点或突变位点的肽段. (3) 基于序列标签的氨基酸突变鉴定算法面临的问题. 基于肽段序列标签的氨基酸突变序列鉴定算法固然能够有效地利用 PSTs 过滤数据库, 弥补 d

32、enovo 测序算法的测序错误并且提高对突变或修饰肽段鉴定的效率和准确性, 但当前已有的 PSTs 算法仍然存在着很多缺乏, 如在 GutenTag 算法中没有考虑同源突变或修饰, 所以不能鉴定出数据库中不存在的突变序列, 而在 Opensea 软件中没有考虑到突变位点的出现可能是由 de novo 的测序错误引起的等. 但是图谱质量是限制序列标签算法的主要因素, 由于低能 CID 碎裂形式通常很难将质量一样或相近的碎裂离子区分开来, 如亮氨酸(Leu)和异亮氨酸(Ile)、赖氨酸(Lys)和谷氨酰胺(Gln)以及苯丙氨酸(Phe)和氧化的甲硫氨酸(Met)等46. 近年来, 随着电子转移解离

33、(electron transfer dissociation, ETD)和高能碰撞解离(high-energy collision induced dissociation, HCD)的出现, 越来越多的比 CID 质谱图质量高的、含有丰富的碎裂离子信息的高精度质谱图被产出, 这些高精度的质谱图能更好地适用于序列标签算法, 提高其准确性. 1.3 基于图谱库搜索的氨基酸突变鉴定算法 在肽段鉴定领域, 图谱库搜索是一种有望取代序列数据库搜索的鉴定策略62. 相比序列数据库搜索策略, 图谱库搜索策略有下面优点: () 直接利用图谱库中每一张真实图谱的各种不同的特征信息进行比对, 如碎片离子峰的峰

34、强度信息、碎裂形式等,使图谱比对算法具有更高层次的灵敏性; () 能够在一个更小、更精到准确的搜索空间内进行搜索, 能够比序列搜索速度快好几个数量级; () 能够轻松地鉴定出图谱库中已存在的变异肽段63. 对于依靠于图谱库搜索的蛋白突变鉴定来讲, 当前最大的限制来源于图谱库的覆盖范围, 尤其是对突变和修饰肽段图谱的包含63,64. 由于在类似的条件下, 肽段的图谱具有可再生性65并且类似序列的肽段通常能够产生类似的质谱图20,66, 因而一批利用图谱库中已收录的肽段图谱来扩大图谱库对肽段的覆盖范围, 以实现对氨基酸突变进行鉴定的算法或工具应运而生. 当前常用的图谱搜索软件及其网址见表 4. 在

35、蛋白质组学中, 图谱库搜索概念早在 1998 年就由 Yates 等人70率先提出, 但由于质谱仪通量不高、生物质谱数据缺乏以及质谱数据的自动化分析方式方法不完善等71原因使得图谱库搜索策略发展缓慢.直到近期 10 年, 随着质谱和计算机技术的快速发展,鉴定出的肽段图谱匹配对(peptide spectrum match,PSM)的数目与日俱增, 图谱库搜索策略才逐步被应用到大规模数据集和数据库中. 近期, 图谱库搜索策略更是被用于开掘样品中的突变肽段. 要用图谱搜索策略来鉴定样品中的突变肽段, 就必需要扩大图谱库对突变肽段的覆盖范围. 当前用于扩大图谱库覆盖范围的算法有 pMatch63、半

36、经历体验算法72,73以及Ji 等人66提出的类似性算法等. pMatch 在构建图谱时,利用肽段已经知道的实验图谱和理论图谱混合构建图谱,用来缓冲由修饰或突变氨基酸残基引起的肽段碎裂形式的变化64. 由 Hu 等人72在 2018 年提出的半经历体验方式方法通过利用图谱库中已收录的PSMs构建突变肽段的质谱图以扩大对突变肽段的覆盖范围. 这种算法把图谱库中图谱对应的肽段序列替换为相应的突变肽段序列, 并将突变肽段的碎裂离子的质核比值替换到图谱中. 2020 年, Ji 等人66提出的类似性算法通过利用类似序列肽段的图谱来推断目的肽段的图谱, 以到达扩大图谱库的覆盖范围的目的. 这种算法的特点

37、是, 通过加权 K 邻近类似算法66(weightedK-nearest neighbor method) 和支持向量机 (supportvector machine, SVM)74, 利用与目的肽段序列类似且长度相等的肽段的图谱来精到准确地预测目的肽段序列的优势碎裂离子(如 b, y 离子类型以及华而不实性丢失离子类型等)的峰强度, 并且利用 SpectraST64,67开创建立的模型构建诱饵数据库进行数据过滤. 同时, Ji 等人66指出, 将此算法应用于扩建美国国家标准与技术研究院(National Institute of Standards and Tech-nology, NIST

38、)图谱数据库, 能有效地将 NIST 图谱库的覆盖率提高 20%60%, 并且用此数据库能够鉴定到样品中更多的突变肽段. 除了通过扩大图谱库覆盖范围以提高图谱库搜索对样品突变肽段的鉴定率以外, 通过改善图谱-图谱匹配(spectrum-spectrum match, SSM)的打分算法也是一条有效提高突变肽段鉴定效率的途径. 当前点积法是 SSMs 打分的主流算法, 如 SpectraST 和X!Hunter68等主流的图谱搜索工具都是利用点积算法进行匹配打分的. 近年来, 一些基于点积法、用于搜索变异肽段的图谱库搜索工具或算法也逐步被开发出来, 如 pMatch 工具63, Bonanza

39、算法20等. pMatch 工具的特点在于, 利用电荷依靠型的质量位移进行离子峰匹配, 并且将常规的点积法与基于概率的模型相结合对图谱间的匹配进行打分. Bonanza算法特点在于, 在挑选候选图谱时, 不限制母离子质量, 能够将不同母离子质量的图谱聚在一起作为候选图谱; 在对离子峰进行匹配时, 不仅将质量相近的子离子峰考虑进去, 还将母离子间的质量差考虑进去; 最后, 利用改良的点积法进行打分. 近期, 考虑到点积法不能提供一个清楚明晰的统计学上的解释并且在打分中忽略了碎裂离子质核比值的差异等缺点,Dasari 等人62构建了一个利用概率评分标准对 SSMs的质量进行评估的搜索引擎 Pepi

40、tome, 并且在错误发现率(false discovery rate, FDR)为 2%的情况下, 成功地鉴定到比SpectraST多10%12%的肽段数目. 前面所提到的图谱鉴定方式方法都需要依靠图谱库, Fu 等人69在 2018 年提出了一个不需要搜索图谱库就能直接对突变肽段进行鉴定的统计学算法 DeltAMT, 此算法通过二维高斯混合模型利用高精到准确度的母离子质量差和保存时间信息对变异肽段进行鉴定. 总体而言, 基于图谱库的蛋白质突变鉴定算法能够有效地缩小搜索空间, 降低搜索时间, 提高搜索的灵敏度. 当前, 由于存在谱图搜索软件对质谱图的整体解析度无法到达传统数据库搜索策略的程度

41、以及谱库的覆盖范围小等原因, 谱图搜索更多的是作为传统数据库搜索策略的互补策略被使用. 但是随着算法的改良以及 PeptideAtlas75计划的进行62, 相信在不久的将来, 利用图谱库对串联质谱进行鉴定的方式方法会越来越广泛地被使用. 2 氨基酸突变鉴定的应用 当调节细胞增殖、分化、死亡的蛋白序列突变累积到一定程度就会引起癌变76. DNA 测序显示, 在复杂的癌症基因组中通常包含 40100 个可能的氨基酸突变位点25, 然而这些突变中只要小部分会对癌症与肿瘤的发生产生作用. 因而假如能够鉴定出与癌症或肿瘤发生有关的突变肽段, 进而对能够真正引发癌症的基因进行重注释就有时机从更深的层次上

42、了解癌症或肿瘤病发的机理, 找寻到治疗癌症或肿瘤的新方式方法. 所以, 提高图谱的解析率, 鉴定出更多的高质量的突变图谱是找寻突变肽段的关键. 受限于质谱数据的质量、计算能力以及当下已经知道SAVs 的覆盖范围等因素, SAVs 鉴定首先应用在小规模的样品数据集中. 2000 年, Gatlin 等人27通过动态构建人类血红蛋白变异数据库初次成功地对人类血红蛋白样品进行了突变鉴定. 随后, 2003 年 Tabb 等人54利用序列标签算法对 32950 张人类晶状体蛋白质样品(human lens sample)质谱数据中的突变肽段鉴定作出了尝试, 成功地鉴定出 742 条肽段, 华而不实 1

43、34条与突变有关. 随着科技的进步和算法的优化,SAVs 鉴定逐步被应用到大规模数据集的鉴定中. 2007 年, Bunger 等人31通过搜索结合 dbSNP 数据库的自建蛋白质数据库从 DU4475 乳腺肿瘤细胞样品质谱数据中鉴定出 629 个 nsSNVs. 同时他们指出,在大规模数据集中, 要想鉴定出高可信的 SAVs, 不仅要依靠鉴定算法还要对假阳性鉴定结果进行过滤,如通过诱饵数据库去除假阳性鉴定等. Tanner 等人77利用 InsPecT 对 1850 万张人类蛋白质样品 HEK293质谱图进行鉴定, 并结合 PTMfinder78算法对鉴定结果进行了假阳性过滤, 从中发现了与

44、 308 个 nsSNVs有关的肽段. 之后, SAVs 的鉴定被广泛地应用于组织、器官等复杂样品数据集中. 2020 年, Hyatt 和 Pan18将 Sipros v2.0 算法应用于鉴定酸性矿坑水(acid minedrainage)环境中的微生物群落蛋白质的突变氨基酸,在含有 57001 个蛋白的数据库中进行搜索, 鉴定出1683 张图谱对应的 755 个突变肽段. 同时, Hyatt 和Pan18指出, 氨基酸突变中有些可能来自于氨基酸的修饰作用, 如在鉴定出的频率最高的突变氨基酸中,谷氨酰胺(Gln)与谷氨酸(Glu)以及天冬酰胺(Asn)与天冬氨酸(Asp)之间都能经过脱氨基作

45、用进行转换. Su等人39利用自定义的人类血浆蛋白质突变数据库,从 290 个亚洲人血浆样品中鉴定出 2029 个 SAVs, 并挑选出 3 对与糖尿病和肥胖有关的 SAVs进行了绝对定量分析, 指出表型不仅和 SAVs 的浓度有关, 也和SAVs 变种的相对率有关系. Mathivanan 等人25通过构建人类蛋白质突变数据库(HPMD), 从直肠癌细胞系中鉴定出 2728 个蛋白, 华而不实有 57 个突变蛋白是初次在直肠癌中被鉴定出来的. 这些新鉴定出来的突变蛋白在发展新的直肠癌生物标志物和研究治疗直肠癌的靶蛋白方面将会发挥宏大的作用. 近年来, 利用 RNA-seq 数据对蛋白质组数据

46、进行鉴定逐步遭到人们的青睐. 2020 年, Wang 等人41基于RNA-Seq 数据构建了蛋白质数据库并对 2 个直肠癌细胞系 SW480 和 RKO 进行了鉴定, 分别鉴定出 18760和 22623 张质谱图. 这些图谱包含 23 条不存在于dbSNP54 中的变异肽段, 华而不实鉴定到的 TP53P309S突变能够增加SW480细胞的增殖能力, 并且能够加强对细胞抗癌药物的耐受性; HSP90AA1D393N突变对致癌蛋白的构象和稳定性有着宏大的影响. 固然利用质谱数据结合序列数据库搜索是当前主流的蛋白鉴定的策略, 但在传统的数据库搜索中,即便利用最好的质谱平台和最优的分析软件, 也

47、有相当一部分质谱图不能被解析出来79,80. 随着越来越多的 PSMs 被鉴定出来, 人们开场利用质谱图数据库来鉴定突变肽段, 并且成功地鉴定到了比序列数据库搜索更多的 SAVs. 在 FDR=0.0001%的条件下,Hu 等人72,73利用 SpectraST 搜索半经历体验图谱库并结合 PeptideProphet81对结果进行检验, 成功地从人类血浆样品中鉴定出了与 SAVs 有关的 2045 条肽段,而一样条件下, X!Tandem 则只从序列数据库中鉴定出来 623 条与 SAVs 有关的肽段. 3 结束语 随着 DNA 测序成本大幅降低, 越来越多个体的基因组序列被鉴定出来82.

48、但即使在知道人类全基因序列信息的情况下, 科学家们对基因型与分子表型间关系的了解也只是冰山一角83. 而对分子表型的了解有助于科学家们对人类疾病发生机理的理解,比方由 RNA、蛋白质以及翻译后修饰数据能够容易地推断出信号通路能否被激活. 固然当前出现了很多能够预测基因突变对蛋白分子构造及功能影响的软件和在线工具, 如 IntOGen84, SIFT85和 Poly-Phen-286等, 但这些预测工具只能辅助性地对突变氨基酸进行挑选和排序, 以便减少实验验证的候选者87. 而结合了变异蛋白信息的基因信息能够有效地帮助科学家对特定生物学经过的分子途径以及疾病发生的机制等进行理解, 进而增加预防、诊断、治疗疾病的手段88. 本文从数据库搜索、数据库搜索与 de novo 结合的序列标签搜索以及新兴的图谱比对搜索方式方法 3 个方面对大规模鉴定突变蛋白的方式方法作出了比拟全面的介绍. 当前, 无论哪一种搜索方式方法都遭到离子碎裂形式理解程度的深切进入、计算能力高低以及数据库覆盖范围大小等因素的限制, 而结合不同搜索方式方法能够实现不同方式方法间的互补, 能有效地提高鉴定结果的灵敏度. Dasari 等人62发现, 将序列数据库搜索和图谱库搜素结合起来