质谱数据在质量评估中的分析方法、装置和存储介质与流程

1.本发明涉及数据处理和质谱技术领域，特别涉及一种质谱数据在质量评估中的分析方法、装置和存储介质。


背景技术：

2.基于质谱的定量蛋白质组学研究，通过表征蛋白质组丰度的动态变化规律，对生命科学及人类许多疾病的发生、发展的诊断有重大意义。相比于标记发现定量蛋白质组学技术，非标记的发现定量蛋白质组学技术中的样本前处理方法简单、通量高，更适合大队列发现蛋白质组学分析。而在非标记发现定量蛋白质组学中，基于数据非依赖模式（data independent acquisition, dia) 的定量蛋白质组学技术，相比于数据依赖模式（data dependent acquisition, dda）的定量蛋白质组学技术，在保证通量和样本前处理简单便捷的基础上，还具有较好的定量重现性和准确性。随着数据分析算法的提升，数据非依赖模式的定量蛋白质组学的数据分析复杂度较高的问题也逐渐改善，数据非依赖模式的定量蛋白质组学近几年开始逐渐应用在基于质谱的定量蛋白质组学中。基于质谱的大队列定量蛋白质组学的数据采集中，质谱仪器的稳定性和高灵敏度是保证质谱产生高质量数据的关键。虽然基于质谱的非标记定量蛋白质组学技术日新月异，但是在大队列的实验分析中，质谱运行异常也会带来原始文件重复性低的问题，因此建立自动化的仪器质量控制流程来保证仪器稳定性也显得尤为重要。
3.因数据非依赖模式产生的原始数据分析困难，目前多数实验室使用通过数据依赖模式的方法产生的原始文件的评估流程对非数据依赖模式产生的原始文件进行质量控制，例如，通过原始文件汇总肽段和蛋白鉴定数目。
4.目前，前人文献报道对基于数据依赖模式方法产生的原始文件的评价方案包含部分参数，但是存在输入参数过多导致用户使用不易上手现象。此外这些方法没有利用实际产出的结果形成一个定量结果，而这一定量结果是仪器性能状态变化的直接体现。
5.在相关技术中，“实时质量控制分析”（qc-art）工具能够动态标记仪器状态或者样本的潜在问题，是对dda原始文件建立的质控标准。在qc-art技术中，研究人员提取dda原始文件中的肽段鉴定数目、母离子质荷比/碎片离子质荷比的偏斜率、标记蛋白质组学中报告离子、鉴定到的肽段的tic（total ion chromatogram，总离子流图）数值中位数、二级谱图采集的四分位数等。然后，用户先自行定义一个基线参数作为标准，通过pca的方式去除冗余的参数加快后续数据处理进程。最后使用rpca（robust principal component analysis，鲁棒的主成分分析）算法结合马氏距离对文件中每个参数给出一个具体的打分。该方式具有与标准人工分析方法类似的准确度，同时还具有实时分析的优势。
6.而对于dia文件的质量控制，现阶段尚无公认的标准，因此很大程度上是依赖人工进行质量控制。


技术实现要素：

7.发明人对相关技术进行分析后发现，qc-art工具存在以下局限性：1）输入参数过多，限制了自身应用的推广；2）输入参数存在局限性，仅适用于orbitrap仪器产生的数据；3）qc-art的参数都是在质谱参数的基础上做一个简单的统计分析，并且其输出结果不包含任何的定量结果，这一结果是评价仪器状态的直接体现，因此使用qc-art需要引入人工的二次评价，需要实验室技术人员具备强大的专业知识，而这有会限制这些指标对系统行能诊断的潜力。
8.对于dia文件的质控，现阶段没有公认的标准，因此很大程度上是人工进行质控。
9.本发明实施例所要解决的一个技术问题是：如何更准确地评估质谱文件的状态。
10.根据本发明一些实施例的第一个方面，提供一种质谱数据在质量评估中的分析方法，包括：确定质谱文件的一个或多个液相色谱系统（liquid chromatography，简称：lc）特征和一个或多个质谱系统（mass spectrometer，简称：ms）特征；根据lc特征、ms特征和产生质谱文件的仪器的信息，生成每个特征对应的输入特征；利用每个特征对应的预测模块对相应的输入特征进行处理，获得每个特征对应的预测值；以及，分别根据lc特征对应的预测值和ms特征对应的预测值，确定质谱文件的lc状态和ms状态。
11.在一些实施例中：lc特征包括总离子流图；或者，ms特征包括肽段鉴定数目、蛋白鉴定数目、一级质谱信号值或二级质谱信号值中的至少一种。
12.在一些实施例中：lc特征还包括目标肽段离子色谱峰保留时间；或者，ms特征还包括目标肽段离子色谱峰采集点数、目标肽段离子色谱峰一级质谱峰面积或目标肽段离子色谱峰二级质谱峰面积中的至少一种。
13.在一些实施例中：lc特征还包括目标肽段离子色谱峰峰形、或者目标肽段离子色谱峰峰宽；或者，ms特征还包括鉴定到的离子电荷分布、离子源喷雾稳定性、肽段母离子质谱鉴定准确度、肽段碎片离子质谱鉴定准确度或一级质谱或二级质谱信号差异度中的至少一种。
14.在一些实施例中，根据lc特征、ms特征和产生质谱文件的仪器的信息，生成每个特征对应的输入特征包括：获取产生质谱文件的仪器对应的编码特征；对于lc特征和ms特征中的每个特征，利用特征和仪器对应的编码特征，生成特征对应的输入特征。
15.在一些实施例中，对于lc特征和ms特征中的每个特征：在特征的维度数为1的情况下，特征对应的预测模块包括一个第一全连接层；或者，在特征的维度数大于1的情况下，特征对应的预测模块包括多层感知机、或者包括lstm网络和第二全连接层。
16.在一些实施例中，在特征为总离子流图的情况下，特征对应的预测模块包括lstm网络和全连接层。
17.在一些实施例中，分别根据lc特征对应的预测值和ms特征对应的预测值，确定质谱文件的lc状态和ms状态包括：利用第三全连接层，对将每个lc特征对应的预测值进行拼接的结果进行处理，获得第三全连接层输出的、质谱文件的lc状态；利用第四全连接层，对将每个ms特征对应的预测值进行拼接的结果进行处理，获得第四全连接层输出的、质谱文件的ms状态。
18.在一些实施例中，分析方法还包括：利用训练文件，对每个特征对应的预测模块进行训练，其中，训练文件具有lc状态的标记值和ms状态的标记值、以及训练文件中每个特征
对应的标记值。
19.在一些实施例中，对每个特征对应的预测模块进行训练包括：利用lc状态的预测误差、ms状态的预测误差、每个特征对应的预测误差以及预设一组或多组特征的一致性误差，计算损失值，其中，特征一致性误差与同一组中的特征对应的预测值之间的差异成正相关关系；根据损失值，对每个特征对应的预测模块进行训练。
20.在一些实施例中：一级质谱信号值和目标肽段离子色谱峰一级质谱峰面积属于同一组特征；或者二级质谱信号值和目标肽段离子色谱峰二级质谱峰面积属于同一组特征；或者肽段鉴定数目和蛋白鉴定数目属于同一组特征。
21.在一些实施例中，在质谱文件的lc特征和ms特征中，第一部分的特征是根据质谱文件的整体确定的，第二部分的特征是根据质谱文件的目标肽段确定的。
22.在一些实施例中，将质谱文件的数据集中，缺失率低于预设阈值的肽段确定为目标肽段。
23.在一些实施例中，质谱文件是数据非依赖模式的文件。
24.根据本发明一些实施例的第二个方面，提供一种质谱数据在质量评估中的分析装置，包括：确定模块，被配置为确定质谱文件的一个或多个液相色谱系统lc特征和一个或多个质谱系统ms特征；输入特征生成模块，被配置为根据lc特征、ms特征和产生质谱文件的仪器的信息，生成每个特征对应的输入特征；预测模块，被配置为利用每个特征对应的预测模块对相应的输入特征进行处理，获得每个特征对应的预测值；以及，状态确定模块，被配置为分别根据lc特征对应的预测值和ms特征对应的预测值，确定质谱文件的lc状态和ms状态。
25.在一些实施例中，分析装置还包括：训练模块，被配置为利用训练文件，对每个特征对应的预测模块进行训练，其中，训练文件具有lc状态的标记值和ms状态的标记值、以及训练文件中每个特征对应的标记值。
26.根据本发明一些实施例的第三个方面，提供一种质谱数据在质量评估中的分析装置，包括：存储器；以及耦接至存储器的处理器，处理器被配置为基于存储在存储器中的指令，执行前述任意一种质谱数据在质量评估中的分析方法。
27.根据本发明一些实施例的第四个方面，提供一种计算机可读存储介质，其上存储有计算机程序，其中，该程序被处理器执行时实现前述任意一种质谱数据在质量评估中的分析方法。
28.上述发明中的一些实施例具有如下优点或有益效果。本发明的实施例在评估质谱文件的过程中，能够从lc和ms两个维度评估质谱文件的质谱质量，也能够根据需要从更细致的lc特征和ms特征维度获得定量的评价结果，从而能够更合理地进行质量控制。该方案能够自动化地执行，从而更高效、准确。并且，上述实施例的通用性强，适用于多种质谱仪器产生的文件的评估。
29.通过以下参照附图对本发明的示例性实施例的详细描述，本发明的其它特征及其优点将会变得清楚。
附图说明
30.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现
有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
31.图1示出了根据本发明一些实施例的质谱数据在质量评估中的分析方法的流程示意图。
32.图2a示例性地示出了第一机器学习模型的结构示意图。
33.图2b示例性地示出了第二机器学习模型的结构示意图。
34.图3示出了根据本发明一些实施例的训练方法的流程示意图。
35.图4示出了根据本发明一些实施例的目标肽段筛选方法的流程示意图。
36.图5示出了根据本发明一些实施例的质谱数据在质量评估中的分析装置的结构示意图。
37.图6示出了根据本发明另一些实施例的质谱数据在质量评估中的分析装置的结构示意图。
38.图7示出了根据本发明又一些实施例的质谱数据在质量评估中的分析装置的结构示意图。
具体实施方式
39.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。以下对至少一个示例性实施例的描述实际上仅仅是说明性的，决不作为对本发明及其应用或使用的任何限制。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
40.除非另外具体说明，否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本发明的范围。
41.同时，应当明白，为了便于描述，附图中所示出的各个部分的尺寸并不是按照实际的比例关系绘制的。
42.对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论，但在适当情况下，所述技术、方法和设备应当被视为授权说明书的一部分。
43.在这里示出和讨论的所有示例中，任何具体值应被解释为仅仅是示例性的，而不是作为限制。因此，示例性实施例的其它示例可以具有不同的值。
44.应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步讨论。
45.图1示出了根据本发明一些实施例的质谱数据在质量评估中的分析方法的流程示意图。如图1所示，该实施例的质谱数据在质量评估中的分析方法包括步骤s102~s108。
46.在步骤s102中，确定质谱文件的一个或多个lc特征和一个或多个ms特征。
47.在一些实施例中，质谱文件是通过dia方式采集的文件，即，是dia模式的文件。
48.在一些实施例中，lc特征包括总离子流图；或者，ms特征包括肽段鉴定数目、蛋白鉴定数目、一级质谱信号值或二级质谱信号值中的至少一种。经过发明人的实验，总离子流图、肽段鉴定数目、蛋白鉴定数目、一级质谱信号值或二级质谱信号值能够更好地评价质谱
文件的状态。
49.在一些实施例中，lc特征包括目标肽段离子色谱峰保留时间；ms特征包括目标肽段离子色谱峰采集点数、目标肽段离子色谱峰一级质谱峰面积或目标肽段离子色谱峰二级质谱峰面积中的至少一种。经过发明人的实验，目标肽段离子色谱峰保留时间、目标肽段离子色谱峰采集点数、目标肽段离子色谱峰一级质谱峰面积和目标肽段离子色谱峰二级质谱峰面积也能够较准确地评价质谱文件的状态。
50.在一些实施例中，lc特征包括目标肽段离子色谱峰峰形、或者目标肽段离子色谱峰峰宽；ms特征包括鉴定到的离子电荷分布、离子源喷雾稳定性、肽段母离子质谱鉴定准确度、肽段碎片离子质谱鉴定准确度或一级质谱或二级质谱信号差异度中的至少一种。经过发明人的实验，目标肽段离子色谱峰峰形、目标肽段离子色谱峰峰宽、离子电荷分布、离子源喷雾稳定性、肽段母离子质谱鉴定准确度、肽段碎片离子质谱鉴定准确度和一级质谱或二级质谱信号差异度也能够较准确地评价质谱文件的状态。
51.下面示例性地描述各个特征的含义。f1：目标肽段离子色谱峰峰形表征（peak shape of targeted peptide precursor）。
52.提取目标离子的分析软件导出的结果，并根据该结果计算每个目标离子的峰形信息，例如不对称因子和拖尾因子。
53.特征f1在相关技术中并不应用于dia领域。
54.f2：目标肽段离子色谱峰峰宽（peak width of targeted peptide precursor）。
55.使用确定f1的流程获得母离子的轮廓图，并对该轮廓图进行峰检测，计算母离子半峰宽信息。
56.f3：目标肽段离子色谱峰保留时间（retention time of targeted peptide precursor）。
57.基于提取目标离子的分析软件与dia数据分析工具软件的分析结果提取每个文件的目标肽段离子色谱图的保留时间。
58.f4：一级离子流图（ms1 chromatogram）。
59.通过质谱文件（mzxml格式）提取一级离子质子数与电荷数的比值（m/z）与强度（intensity），使用样条插值生成n维的离子数值，以生成一级离子流图。
60.在相关技术中，特征f4从未应用于现有的质谱文件质量评估中。
61.f5：鉴定到的离子电荷分布（charge state distribution of identified peptide precursors）。
62.提取dia数据分析工具软件的分析结果，统计每个文件鉴定到的母离子数目、不同电荷数的母离子的数目、母离子的电荷平均值、不同电荷数的母离子的数目占总离子数目的比例。
63.f6：离子源喷雾稳定性（stability of ionization）。
64.提取mzxml文件中的信息，统计所有连续间隔几张一级质谱信号超过n倍的次数，计算出现次数。
65.特征f6在相关技术中并不应用于dia领域。
66.f7：肽段母离子质谱鉴定准确度(mass accuracy of peptide precursor)。
67.通过dia数据分析工具软件的分析结果统计每个文件的质荷比准确度的中位数值。
68.f8：一级质谱信号（ms1 signal）值。
69.对mzxml文件的数据进行插值，统计每个原始文件ms1信号的面积数值。
70.f9：目标肽段离子色谱峰采集点数（data points of each identified peptide precursor）。
71.从搜库软件中或根据采集质谱采集频率和色谱出峰时间计算目标离子在每个原始文件采集的点数。
72.f10：肽段碎片离子质谱鉴定准确度（mass accuracy of peptide fragment ion）。
73.通过dia数据分析工具软件结果统计每个文件的所有肽段二级碎片离子的质荷比准确度的中位数值。
74.f11：二级质谱信号（ms2 signal）值。
75.对mzxml的数据进行插值，统计每个原始文件ms2信号的面积数值。
76.f12：一级质谱和二级质谱信号差异度（ratio of ms1 signal to ms2 signal）。
77.每个文件鉴定的tic ms1/tic ms2比值，即f8/f11的数值。
78.特征f12在相关技术中并不应用于dia领域。
79.f13：肽段鉴定数目（number of identified peptide）。
80.通过dia数据分析工具软件 统计每个原始文件鉴定到的肽段数目。
81.f14：蛋白鉴定数目（number of identified protein）。
82.通过dia数据分析工具软件 统计每个原始文件鉴定到的蛋白数目。
83.f15：目标肽段离子色谱峰一级质谱峰面积（ms1 signal of targeted peptide precursor）。
84.统计每个原始文件鉴定到的目标肽段离子的一级质谱峰面积及在所有原始文件的分布。
85.f16：目标肽段离子色谱峰二级质谱峰面积（ms2 signal of targeted peptide precursor）。
86.统计每个原始文件鉴定到的目标肽段离子的二级质谱峰面积及在所有原始文件的分布。
87.在步骤s104中，根据lc特征、ms特征和产生质谱文件的仪器的信息，生成每个特征对应的输入特征。即，生成的输入特征既有原始特征的信息、又有仪器的信息。
88.在一些实施例中，获取产生质谱文件的仪器对应的编码特征；对于lc特征和ms特征中的每个特征，利用特征和仪器对应的编码特征，生成特征对应的输入特征。
89.仪器对应的编码特征例如通过独热（one-hot）编码方式表示。例如，对于仪器a、b和c，分别使用[1, 0, 0]、[0, 1, 0]、[0, 0, 1]表示。
[0090]
在步骤s106中，利用每个特征对应的预测模块对相应的输入特征进行处理，获得每个特征对应的预测值。预测模块为根据输入特征进行计算、并输出预测值的模块，其可以包括机器学习模型、子网络等等，例如包括全连接层（fc）、多层感知机（mlp）或神经网络模型等。
[0091]
例如，将每个输入特征输入到相应的预测模块中，获得预测模块输出的结果，作为预测值。
[0092]
在一些实施例中，在特征的维度数为1的情况下，特征对应的预测模块包括一个第一全连接层。
[0093]
在一些实施例中，在特征的维度数大于1的情况下，特征对应的预测模块包括mlp、或者包括长短期记忆网络（long short-term memory，简称：lstm）网络和第二全连接层。例如，在特征为总离子流图的情况下，特征对应的预测模块包括lstm网络和全连接层。总离子流图包括时序信息，因此通过lstm网络处理能够得到更准确的预测结果。
[0094]
在步骤s108中，分别根据lc特征对应的预测值和ms特征对应的预测值，确定质谱文件的lc状态和ms状态。lc状态和ms状态能够用于评估质谱文件的质量，其状态值例如包括用于表示“合格的”“不合格的”的数值。根据需要，也可以采用其他的状态值，这里不再赘述。此外，步骤s106中获得的每个特征对应的预测值也能够作为质谱文件评估的定量结果。
[0095]
在一些实施例中，利用第三全连接层，对将每个lc特征对应的预测值进行拼接（concat）的结果进行处理，获得第三全连接层输出的、质谱文件的lc状态；利用第四全连接层，对将每个ms特征对应的预测值进行拼接的结果进行处理，获得第四全连接层输出的、质谱文件的ms状态。
[0096]
在一些实施例中，利用两个机器学习模型实现步骤s106和s108的处理，两个机器学习包括第一机器学习模型和第二机器学习模型。第一机器学习模型包括lc特征对应的预测模块和第三全连接层，第二机器学习模型包括ms特征对应的预测模块和第四全连接层。
[0097]
图2a示例性地示出了第一机器学习模型的结构示意图。
[0098]
在图2a中，特征f1、f2和f3各对应一个mlp，特征f4对应一个lstm和fc。将特征f1~f4对应的输入特征输入到相应的预测模块后，分别获得对应的预测值o1~o4。然后，将o1~o4输入到fc中，获得输出的lc状态。
[0099]
图2b示例性地示出了第二机器学习模型的结构示意图。
[0100]
在图2b中，特征f5、f15和f16各对应一个mlp，特征f6~f14各对应一个fc。将特征f5~f16对应的输入特征输入到相应的预测模块后，分别获得对应的预测值o5~o16。然后，将o5~o16输入到fc中，获得ms状态。
[0101]
本发明的实施例在评估的过程中，能够从lc和ms两个维度评估质谱文件的质谱质量，也能够根据需要从更细致的lc特征和ms特征维度获得定量的评价结果，从而能够更合理地进行质量控制。该方案能够自动化地执行，从而更高效、准确。并且，上述实施例的通用性强，适用于多种质谱仪器产生的文件的评估。
[0102]
本发明的实施例能够应用于对dia文件的评估，从而能够高效、准确地评估dia文件。然而，本发明的实施例也能够对dda文件进行质量评价。
[0103]
在一些实施例中，可以预先对各个特征对应的预测模块进行训练：利用训练文件，对每个特征对应的预测模块进行训练，其中，训练文件具有lc状态的标记值和ms状态的标记值、以及训练文件中每个特征对应的标记值。训练文件是指用于训练的质谱文件。下面参考图3描述本发明训练方法的实施例。
[0104]
图3示出了根据本发明一些实施例的训练方法的流程示意图。如图3所示，该实施例的训练方法包括步骤s302~s310。
[0105]
在步骤s302中，获取训练文件，训练文件为质谱文件，具有lc状态的标记值和ms状态的标记值、以及训练文件的一个或多个lc特征和一个或多个ms特征中每个特征对应的标记值。这些标记值可以是由有经验的从业者根据经验评估的。
[0106]
在步骤s304中，根据训练文件的lc特征、ms特征和产生质谱文件的仪器的信息，生成每个特征对应的输入特征。
[0107]
在步骤s306中，将lc特征对应的输入特征输入到第一机器学习模型，将ms特征对应的输入特征输入到第二机器学习模型。第一机器学习模型和第二机器学习模型的架构可以参见前述实施例，这里不再赘述。根据需要，可以对第一机器学习模型和第二机器学习模型进行一些变形处理。
[0108]
在第一机器学习模型和第二机器学习模型中，利用每个特征对应的预测模块对相应的输入特征进行处理，获得每个特征对应的预测值。并且，第一机器学习模型还根据lc特征对应的预测值确定质谱文件的lc状态，第二机器学习模型还根据ms特征对应的预测值确定质谱文件的ms状态。
[0109]
在步骤s308中，根据损失函数计算损失值。
[0110]
在一些实施例中，利用lc状态的预测误差、ms状态的预测误差、每个特征对应的预测误差以及预设一组或多组特征的一致性误差，计算损失值，其中，特征一致性误差与同一组中的特征对应的预测值之间的差异成正相关关系。
[0111]
一些特征对应的预测值在理论上应当是相同的。例如，设特征t1与t2属于同一组。理论上，如果t1的预测值为1，则t2的预测值也应当为1；如果t1的预测值为0，则t2的预测值也应当为0。
[0112]
在一些实施例中，一级质谱信号值和目标肽段离子色谱峰一级质谱峰面积属于同一组特征；二级质谱信号值和目标肽段离子色谱峰二级质谱峰面积属于同一组特征；肽段鉴定数目和蛋白鉴定数目属于同一组特征。
[0113]
一个示例性的损失函数如以下公式所示：其中，表示损失值，例如交叉熵损失（cross entropy loss）；表示每个特征对应的预测误差的总和，表示特征标识，表示第个特征对应的系数，表示第个特征对应的预测误差（例如为标记值和预测值的差距）；表示一致性误差的总和，其中，和表示属于同一组的特征的特征标识，表示第个特征对应的系数，表示属于同一组的特征的预测值的差异带来的损失，例如，第个特征对应的预测值和第个特征对应的预测值一致时，为0，否则为1；表示lc状态的预测误差、ms状态的预测误差的总和，lc状态和ms状态的预测误差例如分别根据相应状态的预测值和标记值直接的差距确定。
[0114]
从而，可以从多个维度衡量模型的损失，使得训练后的模型具有更高的准确率。
[0115]
在步骤s310中，根据损失值，对每个特征对应的预测模块进行训练。例如，通过反向传播算法对第一机器学习模型和第二机器学习模型的参数进行联合调整，则模型中的预
测模块的参数也得到了调整。
[0116]
通过上述实施例，能够通过预先进行的训练过程获得预测模块的参数，从而在利用预测模块确定质谱文件的状态时，能够具有更高的准确性。
[0117]
在一些实施例中，在质谱文件的lc特征和ms特征中，第一部分的特征是根据质谱文件的整体确定的，第二部分的特征是根据质谱文件的目标肽段确定的。目标肽段是指通过筛选而确定的特定肽段，这些肽段相比于其他肽段更具有代表性，可以体现质谱文件的状态。
[0118]
例如，在特征f1~f16中，f4~f14根据原始的质谱文件的整体确定，f1、f2、f15和f16根据该质谱文件的目标肽段确定。
[0119]
图4示出了根据本发明一些实施例的目标肽段筛选方法的流程示意图。如图4所示，该实施例的目标肽段筛选方法包括步骤s402~s404。
[0120]
在步骤s402中，获取质谱文件的数据集。
[0121]
在一些实施例中，数据集中的质谱文件来自于目前常用于发现蛋白质组学的所有类型的质谱仪。从而，筛选的目标肽段能够适用于更多种类的质谱仪器。
[0122]
在步骤s404中，将质谱文件的数据集中，缺失率低于预设阈值的肽段确定为目标肽段。
[0123]
在一些实施例中，预设阈值为80%。
[0124]
通过上述实施例，能够将出现频率高、有代表性的肽段确定为目标肽段，从而基于目标肽段提取的特征能够更准确地反映质谱文件的特点，使得预测的准确性得到了提升。
[0125]
图5示出了根据本发明一些实施例的质谱数据在质量评估中的分析装置的结构示意图。如图5所示，该实施例的分析装置50包括：确定模块510，被配置为确定质谱文件的一个或多个液相色谱系统lc特征和一个或多个质谱系统ms特征；输入特征生成模块520，被配置为根据lc特征、ms特征和产生质谱文件的仪器的信息，生成每个特征对应的输入特征；预测模块530，被配置为利用每个特征对应的预测模块对相应的输入特征进行处理，获得每个特征对应的预测值；以及，状态确定模块540，被配置为分别根据lc特征对应的预测值和ms特征对应的预测值，确定质谱文件的lc状态和ms状态。
[0126]
在一些实施例中：lc特征包括总离子流图；或者，ms特征包括肽段鉴定数目、蛋白鉴定数目、一级质谱信号值或二级质谱信号值中的至少一种。
[0127]
在一些实施例中：lc特征还包括目标肽段离子色谱峰保留时间；或者，ms特征还包括目标肽段离子色谱峰采集点数、目标肽段离子色谱峰一级质谱峰面积或目标肽段离子色谱峰二级质谱峰面积中的至少一种。
[0128]
在一些实施例中：lc特征还包括目标肽段离子色谱峰峰形、或者目标肽段离子色谱峰峰宽；或者，ms特征还包括鉴定到的离子电荷分布、离子源喷雾稳定性、肽段母离子质谱鉴定准确度、肽段碎片离子质谱鉴定准确度或一级质谱或二级质谱信号差异度中的至少一种。
[0129]
在一些实施例中，输入特征生成模块520进一步被配置为获取产生质谱文件的仪器对应的编码特征；对于lc特征和ms特征中的每个特征，利用特征和仪器对应的编码特征，生成特征对应的输入特征。
[0130]
在一些实施例中，对于lc特征和ms特征中的每个特征：在特征的维度数为1的情况
下，特征对应的预测模块包括一个第一全连接层；或者，在特征的维度数大于1的情况下，特征对应的预测模块包括多层感知机、或者包括lstm网络和第二全连接层。
[0131]
在一些实施例中，在特征为总离子流图的情况下，特征对应的预测模块包括lstm网络和全连接层。
[0132]
在一些实施例中，状态确定模块540进一步被配置为利用第三全连接层，对将每个lc特征对应的预测值进行拼接的结果进行处理，获得第三全连接层输出的、质谱文件的lc状态；利用第四全连接层，对将每个ms特征对应的预测值进行拼接的结果进行处理，获得第四全连接层输出的、质谱文件的ms状态。
[0133]
在一些实施例中，分析装置50还包括：训练模块550，被配置为利用训练文件，对每个特征对应的预测模块进行训练，其中，训练文件具有lc状态的标记值和ms状态的标记值、以及训练文件中每个特征对应的标记值。
[0134]
在一些实施例中，训练模块550进一步被配置为利用lc状态的预测误差、ms状态的预测误差、每个特征对应的预测误差以及预设一组或多组特征的一致性误差，计算损失值，其中，特征一致性误差与同一组中的特征对应的预测值之间的差异成正相关关系；根据损失值，对每个特征对应的预测模块进行训练。
[0135]
在一些实施例中：一级质谱信号值和目标肽段离子色谱峰一级质谱峰面积属于同一组特征；或者二级质谱信号值和目标肽段离子色谱峰二级质谱峰面积属于同一组特征；或者肽段鉴定数目和蛋白鉴定数目属于同一组特征。
[0136]
在一些实施例中，在质谱文件的lc特征和ms特征中，第一部分的特征是根据质谱文件的整体确定的，第二部分的特征是根据质谱文件的目标肽段确定的。
[0137]
在一些实施例中，分析装置50还包括：目标肽段确定模块660，被配置为将质谱文件的数据集中，缺失率低于预设阈值的肽段确定为目标肽段。
[0138]
在一些实施例中，质谱文件是dia模式的文件。
[0139]
在一些实施例中，质谱数据在质量评估中的分析装置可以部署在后端，即服务器侧。用户通过前端设备输入原始的质谱文件，前端设备将质谱文件发送给后端的分析装置进行处理。然后，分析装置将对质谱文件的评估结果返回给前端设备。
[0140]
图6示出了根据本发明另一些实施例的质谱数据在质量评估中的分析装置的结构示意图。如图6所示，该实施例的分析装置60包括：存储器610以及耦接至该存储器610的处理器620，处理器620被配置为基于存储在存储器610中的指令，执行前述任意一个实施例中的质谱数据在质量评估中的分析方法。
[0141]
其中，存储器610例如可以包括系统存储器、固定非易失性存储介质等。系统存储器例如存储有操作系统、应用程序、引导装载程序（boot loader）以及其他程序等。
[0142]
图7示出了根据本发明又一些实施例的质谱数据在质量评估中的分析装置的结构示意图。如图7所示，该实施例的分析装置70包括：存储器710以及处理器720，还可以包括输入输出接口730、网络接口740、存储接口750等。这些接口730，740，750以及存储器710和处理器720之间例如可以通过总线760连接。其中，输入输出接口730为显示器、鼠标、键盘、触摸屏等输入输出设备提供连接接口。网络接口740为各种联网设备提供连接接口。存储接口750为sd卡、u盘等外置存储设备提供连接接口。
[0143]
本发明的实施例还提供一种计算机可读存储介质，其上存储有计算机程序，其特
征在于，该程序被处理器执行时实现前述任意一种质谱数据在质量评估中的分析方法。
[0144]
本领域内的技术人员应当明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用非瞬时性存储介质（包括但不限于磁盘存储器、cd-rom、光学存储器等）上实施的计算机程序产品的形式。
[0145]
本发明是参照根据本发明实施例的方法、设备（系统）、和计算机程序产品的流程图和／或方框图来描述的。应理解为可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的装置。
[0146]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能。
[0147]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的步骤。
[0148]
以上仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。