METHOD OF IDENTIFICATION OF ENTITIES FROM MASS SPECTRA

Information

  • Patent Application
  • 20210241851
  • Publication Number
    20210241851
  • Date Filed
    July 19, 2019
    5 years ago
  • Date Published
    August 05, 2021
    3 years ago
Abstract
A method for determination of identity of at least one entity from a mass spectrum of at least one entity and from additional data from chemical, physical, biochemical or biological analysis of at least one entity for each entity having the steps of a) collecting analytical data from mass spectrum of the entity, b) obtaining a plurality and prevalences of candidate identities, c) calculation of its score, and d) determining the identity of an entity as the candidate identity with the score closest to the score that would correspond to the true identity of the entity. The entity may be any chemical or biological entity, in particular a peptide, a protein, a lipid, a nucleic acid, a metabolite or a small molecule.
Description
FIELD OF THE INVENTION

The invention relates to a method of determination of identity of an entity from mass spectra. The method is useful in proteomics, metabolomics and its applications in proteomics, metabolomics, genomics and transcriptomics.


BACKGROUND ART

Discovery proteomics contains wealth of rare information, often resistant to reliable interpretation. In shotgun proteomics, discovery-oriented subfield of bottom-up proteomics, proteins are enzymatically cleaved into peptides and the digested samples gradually introduced into an analyser, most commonly into a mass spectrometer using liquid chromatography. In mass spectrometry analysis, typically, in each cycle, masses of intact molecules are analyzed, those of further interest isolated, fragmented and second mass analyses performed on fragments, giving MS/MS spectra. The goal of identification is the peptide producing observed MS/MS spectrum and mapping of peptides to proteins concludes the protein identification task.


Because of immense proteome complexity, number of potential interpretations of each fragment spectrum is enormous. Size apart, even unbiased consideration of all interpretations is of limited use as high spectral homology draws many interpretations equal—a phenomenon often witnessed in peptide sequencing de novo. In practice, identification is intentionally biased towards much smaller database of reference peptides, commonly using search engines such as Sequest or X!Tandem. Reference search leaves many spectra uninterpreted, around 75% as recently evaluated on projects in PRIDE repository. Although many unmatched spectra were confidently found to contain peptides with post-translational modifications, majority still remains without interpretation. Several explanations remain open, e.g., splicing alterations, re-arranged or novel genes, however unmatched spectra often exhibit properties of peptides. As human proteomes differ from reference in around ten thousand peptide amino acid sites, some spectra might be naturally attributed to variants and their identification is the aim of the present invention.


The same problems with interpretation of outcomes of analyses of complex samples are encountered in many other fields of chemical, biochemical and biological research dealing with complex mixtures, irrespective of the analytical methods used and irrespective of the structural character of the samples.


SUMMARY OF THE INVENTION

The present invention relates to a method for determination of identity of at least one entity from a mass spectrum of said at least one entity and optionally from additional data from chemical, biochemical or biological analysis of said at least one entity, for each entity comprising the steps of:

    • a) collecting analytical data from mass spectrum of the entity, and optionally collecting additional analytical data from a chemical, biochemical or biological analysis of the entity,
    • b) obtaining a plurality of candidate identities of the entity and obtaining the prevalences of said candidate identities of the entity, whereas for each candidate identity it applies that all candidate identities with a higher prevalence are included in the plurality of candidate identities;
    • c) for each candidate identity of an entity, calculation of its score, said calculation involving at least prevalence of entity, or at least prevalence of entity and agreement with mass spectrum,
    • d) determining the identity of an entity as the candidate identity with the score closest to the score which would correspond to the true identity of the entity.


The candidate identities selected in step b) comprise candidate identities which are a possible or admissible interpretation of the mass spectrum and optionally of the additional data.


The score calculated in step c) and used for finally determining the identity in step d) may have a form of a numerical value (then in step d), usually the highest value of the score determines the identity which is finally determined to be the correct one for the analyzed entity), or another form, such as an interval of numbers, a non-numerical entity, entities with established order, a number with probabilistic interpretation. The skilled person will appreciate that when a form of the score is selected, also the score which would correspond to the true identity of the entity (the ideal score) is selected or determined by the form of the score or by its calculation. E.g., for a number with probabilistic interpretation, 100% probability (or value 1) corresponds to the true identity of the entity. The “true identity” is meant as the real identity of the entity, which is however unknown at the beginning of the process.


Preferably, in the step c), the calculation involves calculating maximal probability of candidate identity. The maximal probability may be the score, or it may be a variable in the calculation of the score.


Preferably, in the step c), the calculation involves calculating probability of candidate identity. The probability may be the score, or it may be a variable in the calculation of the score.


Preferably, in the step c), the calculation involves calculating probability of candidate identity using Bayes' Theorem.


In one preferred embodiment, in the step b) the value of prevalence is calculated based on at least one of population frequency of said entity, probability of modification of said entity in the environment, probability of modification of said entity during the analysis step.


Preferably, in the steps b) and c), the value of prevalence is expressed as prior probability or as prior-like probability.


In a preferred embodiment, in the step d) the determination of identity comprises evaluating whether multiple forms of isotopically labeled peptides were present.


Preferably, the entity is selected from a molecule having the molecular weight of up to 2000 mol/g, a peptide, a protein, a lipid, a nucleic acid, a metabolite.


Preferably, the entity is a peptide, and the method used to obtain the mass spectrum is tandem mass spectrometry (also referred to as MS/MS).


In preferred embodiments which are further referred to as comprising “enumeration”, in step b), the obtaining of the candidate entities and/or of prevalence of the candidate identities comprises enumeration which comprises the steps of:

    • b.a) selecting initial candidate identities with initial prevalences;
    • b.b) transferring said initial candidate identities into a base of candidate identities;
    • b.c) producing new candidate identities by application of events to said base of candidate identities, and incorporating said new candidate identities into said base of candidate identity identities and continuing said producing unless a limiting condition is met;
    • b.d) transforming the base of candidate identities obtained in step b.c) into candidate identities with associated prevalence.


In the enumeration embodiments, preferably said candidate identities are peptides; said prevalence is expressed as a prior-like probability; said initial entities are N-terminally-cleaved linear subsequences of reference proteins; said applicable events comprise modification, substitution and cleavage; said limiting condition is minimal prior-like probability of given form of peptide.


In the enumeration embodiments, preferably said candidate identities are proteins; said prevalence is expressed as a prior-like probability; said initial entities are reference exon-based protein models; said applicable events comprise exon exclusion and exon inclusion; said limiting condition is minimal prior-like probability of exon-based model; said transformation of entities into hypotheses is concatenation of exons into protein-coding sequence and translation in silico.


The method of the present invention has a number of potential utilizations, which may involve additional steps upstream or downstream, or may involve the utilization of the determined identity of one or more entities by the method of the present invention in known methods.


The method of the present invention, wherein the entities are proteins, wherein the step of obtaining the candidate identities of an entity in step b) includes database search in database of peptide variants, may be used for identification of mutant and polymorphic proteins from mass spectra of proteome, with alterations already observed globally on nucleotide level.


The method of the present invention, wherein the entities are peptides, further comprising the steps of: e) matching of entities determined as polymorphic peptides to database of origins, may be used for determination of identity on the basis of variability of known prevalence, in particular for authentication of cell lines or identification of an individual from mass spectra of proteome.


The method of the present invention, wherein the entities are non-host peptides, wherein in the step b) the prevalence is expressed as prior or prior-like probability and wherein the prevalence of non-host peptides is scaled down according to prevalence of non-host organism, may be used for identification of non-host organism of known prevalence from mass spectra of proteome of host organism.


The method of the present invention, wherein the entities are non-host peptides, wherein in the step b) in obtaining the candidate identities, peptides uniquely mapping to non-host organism are added to enumerated peptides of host organism and prevalence of non-host peptides is lower than of any host peptide, may be used for identification of non-host organism of unknown prevalence from mass spectra of proteome of host organism.


The method of the present invention, wherein the entities are donor peptides, wherein in the step b) the prevalence of donor peptides is scaled according to their prevalence among recipient peptides, may be used for identification of proteins originating from grafted tissue in allograft or xenograft.


The method of the present invention, wherein the entities are peptides, the method further comprising the step of: e) selecting somatic mutant peptides attributable to tumour, may be used for identification of presence of a tumour and estimation of its stage from mass spectra of circulating proteins.


The method of the present invention, wherein the entities are peptides, the method further comprising the step of: e) selection and quantification of polymorphic peptides attributable to donor, may be used for monitoring organ transplantation and early detection of transplant rejection from mass spectra of blood plasma or serum of recipient.


The method of the present invention, wherein the entities are peptides, said method further comprising the step of: e) calculating significance of match between two individuals based on polymorphic peptides, may be used for determination of presence of genetic relationship between two individuals from measured mass spectra of proteome.


Furthermore, the invention encompasses a data processing system comprising means for carrying out the steps of the method of any one of the preceding claims.


The invention also encompasses a computer program comprising instructions which, when the program is executed by a computer, cause the computer to carry out the steps of the method of any one of the preceding claims.


Yet furthermore, the invention encompasses a computer-readable medium comprising instructions which, when executed by a computer, cause the computer to carry out the steps of the method of any one of the preceding claims.





BRIEF DESCRIPTION OF THE DRAWINGS

In the drawings which schematically represent the method or its steps and sub-steps, lines with arrows refer to direct or indirect connection between individual units. Dotted lines with arrows correspond, in general, to alternative embodiments. Alternative embodiments are further indicated with addition of alphabetical letters grouping particular alternative embodiment. Reference numbers of subunits within units are formed as a concatenation of the reference number of the main unit, period and reference number of the subunit. Units depicted on drawings are assumed to be either standalone, or part of some larger units. Dotted outlines of blocks correspond to steps.



FIG. 1 is a schematic representation of incorporation of prevalence models into identification methods.



FIG. 2 is a schematic representation of incorporation of prevalence model for reevaluation.



FIG. 3 is a schematic representation of incorporation of prevalence model within identification system.



FIG. 4 is a schematic representation of incorporation of prevalence model for influencing selection of candidate identities.



FIG. 5 is a schematic representation of enumeration.



FIGS. 6a, 6b and 6c are a schematic representation of enumeration of peptides in shotgun proteomics.



FIG. 7 illustrates use of variants for identification of origin.



FIG. 8 is a schematic representation for evaluation of correspondence between entities.



FIG. 9 illustrates an MS/MS spectrum of particular precursor measured using tandem mass spectrometry.



FIG. 10 illustrates behaviour of particular agreement model in shotgun proteomics.



FIG. 11 illustrates behaviour of particular agreement model of true interpretations in shotgun proteomics.



FIG. 12 illustrates behaviour of particular agreement model of random interpretations in shotgun proteomics.



FIG. 13 is an example of distribution of precursor mass difference for true matches.



FIG. 14 is an example of experimental distributions of retention time at given theoretical retention time.



FIG. 15 shows a selection of true matches based on extreme behaviour of retention time.



FIG. 16 shows a distribution of differences of theoretical and experimental isotopic distributions.



FIG. 17 shows an example of combination of precursor mass difference and retention time into one value.



FIG. 18 shows the power of filtering when precursor mass difference, isotopic distribution difference, retention time and protein evidence are combined into single criterium.



FIG. 19 is a schematic representation of particular example of incorporation of prevalence model in shotgun proteomics.



FIG. 20 shows likely incompleteness of exome sequencing data for areas of low sequencing coverage.



FIG. 21 illustrates the family structure for calculation of correspondence.



FIG. 22 illustrates the behaviour of coverage of reference proteins in pairwise comparison.



FIG. 23 illustrates calculation of at least as good match at random between family members.



FIG. 24 illustrates results of identification of tumour-specific circulating proteins.



FIG. 25 illustrates identification of human mutant biomarkers in murine xenograft models.



FIG. 26 illustrates identification of microbial peptides and demonstrates practical use for diagnostics of microbial pathogens in human and animal materials.



FIG. 27 is a schematic representation of enumeration of splice variants in proteomics.



FIGS. 28a, 28b and 28c illustrate correspondence of tumour size versus proportion of somatic variants among identified peptides.





DETAILED DESCRIPTION OF THE INVENTION

“Entity” herein refers to a chemical or biological entity, such as a molecule, substance or organelle. In particular, entity may be selected from a substance, a compound, a lipid, a metabolite, a peptide, a protein and a nucleic acid.


“Prevalence” herein refers to frequency of occurrence of an entity. The frequency of occurrence of an entity refers to its frequency of occurrence in the nature, or in a specific part of the nature which was the source of the measured sample, such as organism, part of organism, specific environment, etc.


Prevalence can be expressed in relative terms, e.g., entity A being more prevalent than entity B, or in absolute terms, such as percentages or amounts of the entity per unit of the sample or of the part of the nature. Prevalence also includes prior probabilities of entities. Prevalence also includes relative probabilistic terms, referred here as prior-like probabilities, wherein the relative differences between entities are the same as for prior probabilities of entities.


“Identity of an entity” herein refers to the determination of structural information about the entity, such as its chemical structure, sequence of amino acids or nucleotides. The structural information may refer to assigning a known structure to the entity, or determining its structure or part of its structure even though previously unknown.


“Candidate identity” herein refers to possible or admissible explanation (or interpretation) of the observed mass spectra and optionally additional chemical or biological data.


“Enumeration” herein refers to method of construction of candidate identities and their prevalence which is based on initial candidate entities and events for their combination. Such events include modifications of the initial entities which may have occurred.


“Score” is a value calculated for each candidate identity. The score may have a form of a numerical value, vector or array of numerical values, interval of numbers, non-numerical entity, entities with established order. Score also includes number with probabilistic interpretation, for example, probability of correctness, p-value, E-value, q-value, maximal probability and their intervals. The skilled person will appreciate that when determining the form of the score, its value which would correspond to the true identity of the entity is also determined. E.g., for a score corresponding to probability, the value corresponding to the true identity of the entity is from 1 to 100%.


“Mass spectrum” refers to mass spectrum (MS) obtained by introducing the entity into a mass spectrometer and performing the mass spectrum measurement, or to MS/MS spectrum. Analytical data from mass spectrum are typically the data about the fragment peaks shown in the spectrum (m/z values, intensities). Additional criteria from the mass spectrum may also be used, such as precursor mass difference, isotopic distribution difference, protein evidence.


“Chemical, physical, biochemical or biological analysis” include any analytical methods allowing to obtain data useful for the determination of identity of the entity. Such methods include spectroscopic analytical methods such as NMR spectroscopy, X-ray diffraction spectrometry, IR spectroscopy; immunochemical methods; optical observation methods; methods relying on interaction with further agents, such as antibodies, labels.


“Explanation” and “interpretation” are herein used to designate the assignment of the identity of at least one entity to an analytical method outcome, i.e. to a mass spectrum and optionally additional data.


The present invention describes a method of determination of identity of an entity based on their mass spectrum data, and optionally additional data from other analytical methods, said method utilizing prevalence data and prevalence or probabilistic calculations. Use of prevalence provides additional layer of discrimination and thus helps in resolution of otherwise indistinguishable situations. For instance, it is often the case that there are many explanations which agree equally well with an observed mass spectrum and additional data. The use of prevalence models might enable between these explanations if one explanation is much more prevalent than the rest. In effect, the utilization of prevalence reduces the complexity of the identification task.


If the way the prevalence is modelled, which is herein referred to as a “prevalence model”, exhibits desirable properties (e.g., the prevalence is expressed in relative probabilistic terms), the candidate interpretations might be often assigned a probability of correctness, or a maximal probability of correctness. Probability of correctness of an explanation has in turn the advantage of being usable in real-life scenarios as it enables long-term modelling of decision-making processes. Similarly, maximal probability of correctness provides strong grounds on which to rule out candidate explanations with direct real-life applicability. This might be shown in contrast to statistical significance of agreement (e.g., p-value or E-value), which does not possess such quality and even highly significant agreement might be often assigned to incorrect interpretations. This behaviour is practically impossible for a properly derived probability of correctness and thus gives practically much more desirable guarantees. The utilization of prevalence models especially helps in reliable identification of rare events. In case the prevalence of candidate explanations varies significantly (e.g., spanning multiple orders of magnitude as is the case in bottom-up proteomics), incorporation of prevalence might be necessary to obtain reliable results.



FIG. 1 represents several basic configurations of incorporation of prevalence model into identification system. In some embodiments 101, the prevalence model 101.2 is integrated into the identification system 101.1. Such incorporation is preferable for derivation of probability of correctness of candidate identity. More specific embodiments are illustrated on FIG. 3.


In other embodiments 102, the identification system 102.1 is separate from a system 102.2 comprising prevalence model, in this configuration the system 102.2 comprising prevalence model process results from said identification system 102.1. Such embodiments are usable for instance to derive maximal probability of candidate identity or probability of candidate identity. More specific embodiments of this kind are further illustrated on FIG. 2.


Yet other embodiments 103 comprise an identification system 103.2 and a system 103.1 comprising prevalence model in which the identification system 103.2 works with the selection of candidate identities influenced by the prevalence model 103.1. Such embodiments can be used to preselect candidate identities in a way which improves the behavior of the identification system. More specific embodiments of this kind are illustrated on FIG. 4.



FIG. 2 represents incorporation of a prevalence model for reevaluation of candidate identities. In this embodiment, evaluated candidate identities 201 pass through a system comprising the prevalence model 202. There are various possible alternatives. In an alternative embodiment 203.A, candidate identities are evaluated utilizing information from prevalence model. In such reevaluation, for instance, a new information can be added; for example the number of candidate identities which are at least as prevalent as the hypothesis and have the same agreement with the observed data as the candidate identity. In another embodiment 203.B, the candidate identities are assigned the maximal probability of their correctness. Specific embodiment of this kind in shotgun proteomics is illustrated on FIG. 27. In some embodiments 203.C, the candidate identities are assigned the probability of their correctness.



FIG. 3 represents embodiments used for determination of identity in identification wherein the prevalence model is integrated within the identification system. This configuration is in general suitable for scoring and derivation of probability of correctness of the candidate identity. In some embodiments, the identification system 302A comprises a true agreement model 302A.1, a random agreement model 302A.2 and a prevalence model 302A.3. Such configuration is particularly suitable for derivation of probability using Bayes' Theorem. In some embodiments, the identification system 302B comprises an agreement model 302B.1 and a prevalence model 302B.2 to obtain score or probability of candidate identity 303.



FIG. 4 represents incorporation of prevalence model to influence selection of tested candidate identities. In some embodiments 403.A, the selection of candidate identities is influenced based on their prevalence. An example in shotgun proteomics is selection of peptides more prevalent than peptides with some modification (for example, methylation), or amino acid substitution, or peptides more prevalent than peptides resulting from splicing alteration. An example in top-down proteomics is selection of proteins similarly prevalent as non-modified proteins.


In alternative embodiments 403.B, the selected candidate identities are at least as prevalent as the candidate identities accepted initially for testing (hypotheses 401). An example in bottom-up proteomics is when candidate identities 401 for testing correspond to variant peptides, and candidate identities which are at least as prevalent as the variant peptides 402 are selected in step 403.B (based on particular assumptions over prevalence of individual candidate identities).


The first step of the present invention comprises collecting analytical data. The methods for collecting analytical data, in particular mass spectrometry data, are well known to those skilled in the art. For example, in shotgun proteomics, the sample preparation protocols are well established and in general process samples into mixture of proteolytic peptides; see for instance an article comparing three protocols FASP, SP3 and iST (Sielaff et al. (2017): Journal of Proteome Research, 16(11): 4060-4072). Identification of the substances starts with their physical separation using liquid chromatography, coupled with mass spectrometer. Substances eluting at a particular time (retention time) enter mass spectrometer, undergo ionization and their masses are measured, giving a precursor spectrum. In data dependent acquisition, after each such measured precursor spectrum, masses of several most abundant precursors are selected, the ions are isolated, fragmented and MS/MS spectra (also fragment, or product spectra) acquired. These fragment spectra are of interest because calculation of the score may comprise evaluation of the agreement of theoretical spectra of candidate identities and the observed spectra of the entities.


In the second step of the method of the invention, candidate identities of the analyzed entity are obtained. This step can be performed obtained in multiple ways.


In the usual scenario, candidate identities are obtained through a database search of entities for the given samples. For example, the search may be for peptides or nucleic acids or lipids or compounds or metabolites for given analyzed organism. Often, candidate identities are obtained through reference database search containing reference entities (e.g., peptides) for the analyzed organism. Examples of such databases are UniProt and ENSEMBL. If the analyzed entities are proteins or peptides, then proteins from these databases are in silico digested with a protease used in the experiment. As reference proteolytic peptides are of highest prevalence, they are self-contained in the sense that all more prevalent peptides (than the peptide of the lowest prevalence) are considered as well. However, if some modifications of the reference entities are considered, care must be taken such that all modifications at least as prevalent as the modification of lowest prevalence are considered as well.


In some embodiments, the candidate identities may be obtained using enumeration of candidate identities.



FIG. 5 schematically shows the general process of enumeration, in which initial entities and events applicable to entities (e.g., chemical modifications occurring in the nature) are used for construction of prevalence model. In the beginning, initial entities 501 with associated prevalence are transferred to the base 502 of entities. The base 502 of entities is a part of a cycle, unlike initial entities. Entities from the base 502 are subjected to events 503 (in silico) which create additional entities, which are incorporated into the base 502. This continues until a pre-defined criterion 504 is met. Then the process stops, entities in the base 502 are optionally transformed in step 505 (if necessary) into the final form which then constitutes the prevalence model 506. This process has an important advantage, when coupled with prevalence: for each candidate identity e enumerated, all candidate identities which are at least as prevalent as said candidate identity e are enumerated as well.


Example of embodiments related to shotgun proteomics, i.e., wherein the analyzed entities are proteins, is used to describe some preferred features of the enumeration. The enumeration shown in FIG. 6a is used for assignment of prevalence to reference peptides, variant and modified peptides of varying cleavage specificity. The enumeration is performed over each reference protein independently and the behaviour for particular reference protein is described as follows. As initial candidate identities for a reference protein, all N-terminally cleaved sequences of said protein are used. The prevalence of these candidate identities depends on the probability of cleavage after residue just before the cleavage point (herein, a0, in FIGS. 6a, 6b and 6c). For example, in case of tryptic digestion, the initial prevalence will be usually large in case of lysine and arginine. In this example, if it is at N-terminus of the protein, the initial prevalence equals 1 (no cleavage needed). These initial candidate identities are transferred to the base of candidate identities. The events applicable to candidate identities are as follows: extension, modification and cleavage. Extension refers to the event of incorporation of next residue in reference amino acid chain and the probability of extension is derived as a complementary event of cleaving. Cleavage is modelled as cleaving after a specific amino acid and each candidate identity needs an exactly one cleavage to become a fully formed candidate identity (such cleavage does not need to happen at the C-terminus of the protein). Modifications (m1, . . . , mj) with respective prior-like probabilities (p1, . . . , pj) are applicable to each amino acid. Further, statistical independence of events is assumed, which enables assignment of prevalence in the form of prior-like probability to every peptide by multiplication of prior-like probabilities of events. The process continues until the minimal prior-like probability is met, which constitutes the stopping criterion. Herein, the entities themselves were candidate identities, therefore there is no need for any transformation step, and thus the base of entities is then taken for the prevalence model.


In some embodiments, prior-like probabilities are involved in the prevalence models and/or in the calculation of the score. Prior-like probabilities are also referred to in literature as relative probabilities. For prior-like probabilities, the relative proportions between individual prior-like probabilities are the same as for prior probabilities. Thus the following holds: suppose selection of large number n of outcomes (e.g., MS/MS spectra) and their interpretation by candidate identities q and r with prior-like probabilities Prq and Prr, respectively. The proportion of correct interpretations by candidate identity q as compared to r is then






Pr
q
·Pr
r
−1.  (1)


Thus prior-like probabilities preserve relative differences between prior probabilities.


As an example of preferred features and prior-like probabilities, the establishment of prior-like probabilities in proteomics is described. Prior-like probabilities can be derived from experimental data under these assumptions: the measured data represent the whole population; and the subset of data which is assumed to be correctly interpreted does not change the distribution.


Function Γ is used from measurements to peptides





Γcustom-charactercustom-character  (2)


which gives for each measurement s (MS/MS spectrum) a true interpretation q (peptide) of the spectrum (assuming that there is exactly one true interpretation).


For including the modifications of the peptides, the following is used: Suppose peptide q represented as a tuple






q=
custom-character
a
1
⊕m
1
, . . . ,a
n
⊕m
n
custom-character  (3)


where each ai is a coded amino acid residue and each mi is a modification applicable to the residue ai. The set of applicable modifications to ai is denoted as Φ(ai) and for technical brevity, existence of empty modification is considered.


Denoting by uppercase Q the set of all modified forms of peptide q, the proportion rq of a particular form of peptide q to all its forms Q is expressed as










r
q

=



P


(

q
=

Γ


(
s
)



)






q

Q




P


(

q
=

Γ


(
s
)



)




.





(
4
)







In practice, there is not enough data to model each peptide separately. It is assumed that their behaviour is independent of peptide sequence and depends only on the modification m itself. For the purpose, peptides can be considered with exactly one residue to which m is applicable, denoting such set as Hm. Denote






custom-character
m
+custom-characterm  (5)


the peptides with the modification. Then, for a particular modification m, the










r
m

=


P


(


Γ


(
s
)





m
+


)



P


(


Γ


(
s
)





m


)







(
6
)







is equal to (4) for each peptide q and can be calculated from larger set of data.


Furthermore, the approach can be extended to account for peptides with varying number of modifiable residues. Such extension behaves in the same way on peptides with exactly one residue and enables the utilization of the whole set of interpretations. Specifically, the proportion for modification m on residue a is derived as the total number of a residues modified with m to the total number of a residues with any applicable modification (also the empty one). Suppose






n(a⊖m)  (7)


refers to the number of a residues with modification m. Then the proportion rm can be derived as










r
m

=


P


(


Γ


(
s
)





m
+


)



P


(


Γ


(
s
)





m


)







(
6
)







As another example, probabilities of DNA/RNA substitution are derived. Derivation is analogous as for modifications but with the following difference in the modelling approach. Due to the low ratio of substitutions in the data, the substitution event is modelled in an aggregated manner (independent of the residue). Specifically, the proportion r of all altered residues to all residues is obtained










r
m

=



n


(

a

m

)






m





(
a
)






n


(

a

m

)




.





(
8
)







and interpret it as a probability of amino acid substitution.


As yet another example, for derivation of cleavage probabilities (after particular amino acid) peptides with missed cleavages and semi-specific cleavage (specific at N-terminus and not specific at C-terminus) were utilized. Denoting ncleavage(a) as the number of residues a followed by cleavage, and n(a) as the total number of residues a, the proportion ra of cleavages after residue a is










r
a

=




n
cleavage



(
a
)



n


(
a
)



.





(
10
)







As yet another example, calculation of prevalence of peptides in allografts or xenografts is shown, for identification of peptides originating from grafted tissue in tumour xenograft model. In this case, the prevalence of peptides from different organisms (donor and recipient) is different and needs to be taken into account. Often the recipient would be an animal model and the donor would be a human; peptides from the animal are expected to be more prevalent than those of the human, for instance to monitor tissue/organ grafting or rejection or identification of donor's peptides in the recipient. Alternatively, the recipient can be a patient undergoing tissue transplantation and the donor can be a tissue/organ donor. In the following, it is shown how to estimate the difference in the prevalence of donor and recipient peptides.


The relative difference in the prevalences of donor peptides to recipient peptides can be estimated through derivation of the origin of homologous peptides of the donor and the recipient. Suppose a homologous peptide attributable to both donor and recipient was identified. The interest is in knowing whether the peptide is from the donor or the recipient. For this purpose, protein evidence (of donor proteins and of recipient proteins) of a given peptide can be used which provides the evidence of the origin of the peptide. The proportion p is estimated as the proportion of the homologous peptides with the donor protein evidence as compared to those with the recipient protein evidence. In the construction of the protein evidence, the protein evidence is restricted to heterologous peptides only. In another approach, the proportion p is estimated as the ratio of detected heterologous peptides. Both approaches can be used when there is a limited homology between donor and recipient, which is often the case in xenografts. In allografts, the proportion can be set equal. From a practical perspective, the relative difference between prevalence of donor and recipient peptides is rather small; for instance the number of the donor peptides is in the order of tens of percent of those from recipient. This is important to note as it simplifies identification of donor peptides as there is no other organism (other than donor) expected to be of higher prevalence than that of recipient.


As another option, determination of prevalence of peptides of non-host organisms is described. The situations when identification of non-host organism is of interest include for example detection of microbial presence in an organism, for example for diagnosis of microbial infection.


In some situations, prevalence is known. In general, it is necessary to take the prevalence (or its estimate) of the non-host organism into account. The situation is partially similar to allografts or xenografts, however with the difference that the prevalence of peptides of the non-host organism is generally lower than that of a grafted tissue and non-host peptides are phylogenetically more distant. This has some consequences, notably that all non-host organisms of higher prevalence need to be considered as well (among other at least as prevalent peptides). If the estimate of prevalence (po) of the non-host organism o is known and also the prevalences of all organisms q at least as prevalent as the organism o, po≤pq, then the prevalence model can be easily configured as follows. The prevalence should be expressed in prior or prior-like probabilities and then the prevalence of non-host peptides of the organism o is multiplied with the value of prevalence po.


Often, the prevalence of the non-host organism is unknown, which complicates the identification task. It is nevertheless possible to increase the specificity of identification, if it can be ascertained that all at least as prevalent organisms are considered as well. In case of unknown prevalence, one possible resolution is to consider all known organisms. This is because it is rather unlikely that organisms which were not described yet are more prevalent than the one we are interested in identification.


In the strictest scenario, all peptides of the non-host organism of interest are assumed of the lowest prevalence among all peptides of all considered organisms. An example of realization of such scenario is to enumerate all organisms with a limiting prior-like probability r (e.g., and estimated prior-like probability, such as r=4·10−6) and further scale down prevalences of non-host organisms. In particular, the prevalence of the host remains identical, but the prevalence of all non-host organisms except the one of interest is multiplied by r, and the prevalence of the non-host organism of interest is multiplied by r2. In such case, the non-host peptides of interest are of strictly the lowest prevalence. As the prior-like probabilities are just very roughly established, it is preferable to calculate the maximal probability Pmax instead of the probability; furthermore in the calculation of Pmax, only the relative order of prior-like probabilities is preferred to be assumed. This scenario requires enumeration of peptides for all known organisms. In such case, the database for r≈4·10−6 will likely have thousands of terabytes.


A more preferred approach to the problem of identifying organisms of unknown prevalence was developed. The computational aspect can be significantly simplified under the following assumptions: i. if a peptide of the non-host organism was measured, it was a fully specific (e.g., tryptic) reference peptide, and ii. all non-host organisms are of equal prevalence which is lower than that of any enumerated peptide of the host. Under these circumstances, the peptides exclusively attributable to the non-host organism of interest can be preselected and added to obtained peptides in enumeration of host.


The third step of the method of the invention, the score is calculated for each candidate entity.


First, agreement metrics used for observed and theoretical spectra is described. Various metrics of agreement between observed and theoretical spectra exist, such as simple number of matching peaks consisting of singly charged ions (b, y ions for CID and HCD). The matching of peaks occurs for prespecified fragment mass tolerance depending on instrumentation used (e.g., 0.3 Th for linear ion trap). The skilled person is aware of the available options and is capable of choosing the suitable one. Derivation of maximal probability of candidate identity can be based on adequacy of the agreement model. The following corresponds to the assertion that a higher agreement of predicted and observed spectrum leads to an increase in the probability of interpretation being correct. Thus for every spectrum o ∈O, given two agreements





α≤β  (11)





it is assumed that






P(p=Γ(o)|Φ(p,o)=α)≤P(q=Γ(o)|Φ)(q,o)=β)   (12)





where





Γ:custom-charactercustom-character  (13)


is a function giving correct interpretation for a spectrum and





Φ:custom-character×custom-charactercustom-character  (14)


an agreement model (for some set X, on which there is order). Furthermore, it is assumed that the same agreement for different candidate identities gives an equal probability of being the correct interpretation from the viewpoint of agreement:






P(p=Γ(o)|Φ(p,o)=α)=P(q=Γ(o)|Φ(q,o)=α)  (15)


For the candidate interpretation q of o, with the prior-like probability Prq, there is a need to obtain a set R of all at least as good interpretations. R is of the following form






R={
custom-character
p,Pr
p
custom-character
|pcustom-character and Prp≥Prq}  (16)


wherein p is a candidate identity and PrP its prior-like (or prior) probability. Thus for the establishment of Pmax, only candidate identities which are a priori equal or more likely need to be considered.


If the correct order of prior-like probabilities is considered, then the maximal probability Pmax of q is inversely proportional to the number of the at least as good interpretations, thus






P≤|R|
−1
=P
max.  (17)


If numerical aspects of prior-like probabilities are assumed, then Pmax is the proportion of Prq among all the at least as good interpretations, thus:










P



Pr
q








p
,

Pr
p





R




Pr
p




=


P
max






R



-
1


.






(
18
)







The value of Pmax is independent of search space size.


Preferred methods for derivation of the probability of candidate identity, with prior-like (or prior) probabilities, are described herein. The probability that a particular candidate identity h is the true candidate identity h=δ(o) of the spectrum o, given its agreement d=Φ(h,o) is:






P(h=Γ(o)|Φ(h,o)=d).  (19)


Such probability can be derived, for example, using Bayes' Theorem, wherein (19) equals to












P


(


Φ


(

h
,
o

)


=


d
|
h

=

Γ


(
o
)




)


·

P


(

h
=

Γ


(
o
)



)




P


(


Φ


(

h
,
o

)


=
d

)



.




(
20
)







The equation (20) can be easily changed to incorporate prior-like probabilities, instead of prior probabilities. In case of prior probabilities, the following holds for each o ∈O, if exhaustively all candidate identities h ∈H are considered:













h






P


(

h
=

Γ


(
o
)



)



=
1.




(
21
)







In case of prior-like probabilities, the sum (21) might be different, but because prior-like probabilities preserve relative differences, they can be always normalized by resealing.


Prior-like probabilities are easier to establish, however it might be not clear how they should be rescaled. If candidate identities are selected, such that the true candidate identity is among them, then prior-like probabilities can be rescaled to sum to 1 and then are equivalent to prior probabilities. If the true identity of the analyzed entity is not within the given candidate identities






custom-character
0custom-character  (22)





and therefore













h






P


(

h
=

Γ


(
o
)



)



=

c
<
1.





(
23
)







The variable c in (23) then corresponds to the probability that the true identity of the analyzed entity is within the selected candidate identities H0. Then, prior-like probabilities of selected candidate identities H0 can be rescaled (their sum) to c and will be equal to prior probabilities.


For instance, in shotgun proteomics, not all candidate interpretations of the spectrum are considered and therefore the true interpretation might be not among them. However, around 25% of spectra are usually correctly interpreted in standard bottom-up proteomics experiment. Thus the value of c≥0.25 and c≤1. This also limits the possible range of prior probabilities in shotgun proteomics and in turn the range of probability of correctness.


In mass spectrometry, multiple additional (supporting) criteria, e.g., precursor mass difference, can be directly used in identification. Furthermore, these criteria are useful for identification of rare events, for instance variant peptides. For practical purposes, it is often desirable for these criteria to have statistical interpretation for simplification of making decisions. In particular, the probability was modelled that the true interpretation of a spectrum has a specific additional/supporting criterion at least as extreme as was observed. This in effect enables removal of interpretations.


An example of selection of true interpretations in proteomics field in experiments over a particular organism (for example, human) is mentioned here. There is no need for the interpretations to be correct, rather that the distribution of these criteria to be same as for correct interpretations. Therefore some degree of incorrect results (e.g., 10%) will very likely not invalidate the results. It is assumed that all tryptic reference peptides of expected organism (e.g., human) with sufficient spectral significance (herein, E-Value of 0.1 in X!Tandem in the reference peptide database search) selects suitable interpretations for modelling of supporting criteria.


In some embodiments, precursor mass difference is used as an additional criterion. Distribution of differences between observed and calculated mass of peptide for true interpretations can be readily calculated. Further, association of probabilistic interpretation to differences enables their direct use in identification.


It is assumed that the distribution D of n differences between observed and calculated precursor mass for correct spectral interpretations matches






D=
custom-character
d
1
, . . . ,d
n
custom-character.  (24)


The number n is often rather large (order of thousands, or tens of thousands) for a particular sample, or even for a single run on modern instruments (such as Orbitrap). Therefore it is not even necessary to model the distribution and thus it is possible to work directly with data, e.g., through percentiles. For probabilistic interpretation of the difference d, D is utilized to calculate pd as proportion of true matches having at least as extreme difference as is d. Thus for difference d, we are interested in










p
d

=





{


e

D

|

e

d


}





D



.





(
25
)







Thus if it is unlikely (e.g., pd at most 0.01) for true match to have at least as extreme precursor mass difference, it provides probabilistic grounds for removal of interpretation.


Mass spectrometry is in modern settings coupled to liquid chromatography which enables utilization of the predicted and observed retention time, similarly as the precursor mass difference. In practice, it is also beneficial to have a statistical interpretation of the difference between those two. In the simplest cases, the retention time difference can be modelled exactly as a precursor mass difference explained above. The prediction of the retention time can be done, for example, via BioLCCC (Liquid Chromatography of Biomacromolecules at Limiting Conditions; http://theorchromo.ru/).


It is preferred, however, to model the behaviour more locally, especially because the observed retention time often shows a non-linear behaviour as related to predicted time. The modelling can be performed for each predicted time ti separately and for each ti construct the distribution Di. Each Di consists of an experimental times (experimental counterparts of tj, tj being neighbors of ti). Each Di contains 2·w neighbors, where 2·w is the window size (the preferred size is 500):






D
i
=
custom-character
e
i−w
, . . . ,e
i
, . . . ,e
i+w
custom-character.  (26)


Then for some theoretical time t and experimental time e the position of e within distribution Di is obtained wherein Di is selected such that its corresponding ti is closest to t. A percentile q is obtained










q
=




{


f


D
i


|

f

e


}






D
i





,




(
27
)







representing where e is within distribution Di. It is of interest to remove an expected proportion of true results based on their retention time. It is assumed that a symmetric removal of those which are on both sides of the distribution is needed and then q is converted to






p=2·(½−|½−q|),  (28)


which gives the required proportion. Thus, for instance, selecting results with p≤0.1 is then expected to give 10% of results with the largest differences (on both sides).


In tandem mass spectrometry, precursor spectra are often also measured and thus difference between the theoretical isotopic distribution and the observed one can be readily calculated as well. The difference can be also associated to statistical interpretation, analogously as for precursor mass difference. The software Isotopic Pattern Calculator (http://isotopatcalc.sourceforge.net/) can be used for the prediction of theoretical isotopic distributions.


Although multiple ways to calculate difference between distributions exist, we will utilize very simple one. For calculation of difference, peaks of theoretical and experimental distributions are initially mass-matched to some precursor mass tolerance (for instance 5 ppm on Orbitrap). Intensities of both distributions (experimental and theoretical) are normed to the sum to one and from this alignment the sum of squares is calculated from difference of intensities. The proportion of true results having at least as extreme difference as expected is then calculated. The calculation can be carried out in the same way as for precursor mass difference.


In the specific example of bottom-up proteomics, proteins are enzymatically digested into peptides and therefore in the resulting mixture it is expected that all peptides (of a particular protein) are present. This is called “protein evidence”. It is therefore unlikely that just one peptide of a protein is identified, and this behavior may be modelled. Although multiple options for modelling of protein evidence exist, the modelling is restrained just to the presence or absence of different protein evidence (e.g., by assigning zero and one, respectively). Thus the probability of true match having p as the extreme protein evidence is:









p
=





{


r

R

|


e
r


e


}





R



.





(
29
)







In practice, the p≤0.1 for no protein evidence and p=1 for protein evidence. The task can be performed even before the step of protein inference, stating whether there exists a particular reference protein isoform for which there exists another peptide.


Additional/supporting criteria (e.g. precursor mass difference, retention time, isotopic distribution difference, protein evidence) can be combined to obtain a single criterium with a desirable statistical interpretation. This criterium is built in a way that it is expected to remove a desired proportion of true matches.


For the particular criteria ci of interest (e.g. precursor mass difference and retention time) of some peptide-spectrum match






custom-character
c
1
, . . . ,c
n
custom-character  (30)





a single value is calculated






c=Πc
i.  (31)


For each c the proportion f(c) of the retained true matches among the results is calculated:










f


(
c
)


=




{


r

R

|


c
r


c


}





R







(
32
)







For a new instance of peptide-spectrum match with criteria






custom-character
d
1
, . . . d
n
custom-character.  (33)


its single value d is calculated as in (31), and then its value of f(d). This can be done by looking up the closest value c to d, and obtaining f(c) (this, for example, equals the behavior of k-nearest neighbor with one neighbor). Of course, in this step various options exist, however, as often a lot of data is available, k-nearest neighbor with one neighbor is preferred due to simplicity and obvious interpretation. Then we can filter out the results with the expected loss of a desired proportion of true interpretations; based, however, on all additional/supporting criteria used.


The fourth step of the method of the invention relates to determining the identity of the analyzed entity.


Maximal probability Pmax of interpretations provides a rationale for removal of unlikely matches with a predictable long-term behaviour. For instance, selection of a large number n of candidate interpretations with Pmax=p is expected to result with at most n·p correct interpretations. Therefore, the knowledge that a particular interpretation is less likely than a pre-determined value provides a rationale for its removal.


The most strict settings of Pmax corresponds to interpretations of Pmax=1, and are a preferred embodiment.


The probability P of candidate interpretations provides a rationale for selection of matches with a predictable long-term behaviour. For instance, selection of a large number n of candidate interpretations with a probability higher than p, is expected to result in at least n·p correct interpretations.


The probabilistic interpretation for the additional/supporting criteria is built in a way to express how likely it is that the true interpretation has the supporting criteria as extreme as observed. If it is therefore unlikely (e.g., up to 10%) that the true interpretation would have as extreme criteria, then by removal of these interpretations it is expected that the same proportion (e.g., up to 10%) of correct matches can be removed.


In the following paragraphs, several uses of the method of the invention are described.


In one preferred embodiment, the method of the present invention can be utilized for matching to databases of origins. The following section describes matching of identified peptide or nucleic acid variants of known prevalence to database of origins, with each origin containing set of variants; FIG. 7 schematically describes the process.


For an analyzed sample s, we are interested in its true origin Γ(s) and agreement Φ(s,Ci) of the sample s and a candidate origin Ci can be used for its establishment. Further, the sample s is considered as a set of variants {v1, . . . , vk} identified in the sample s, denoted as






v(s)={v1, . . . ,vk}.  (34)


The agreement Φ(s,Ci) can be, for instance, a number of matching variants. However, it is more preferable to define the agreement as











Φ


(



i

,
s

)


=




v




i



v


(
s
)












ϕ


(
v
)



-
1




,




(
35
)







wherein φ represents the prevalence of a variant, because it provides the probabilistic interpretation over the obtained results. Rescaling of (35) such that the sum over all Ci equals one (the true origin is within the considered origins), gives the probability of the origin Ci being the true origin:






P(Γ(s)=custom-character|Φ(custom-character,s)=x).  (36)


Herein, a rationale for the use of the formula (35) for derivation of (36) is provided. Foremost, consider the agreement Φ(s,Ci) to be the actual matching variants:





Φ(custom-characteri,s)=custom-characteri∩v(s).  (37)


Said agreement (37) is used to enable derivation of the probability of determination of origin. Suppose, v(s)={va,vb} variants identified in sample s and suppose they are identified correctly (the true origin has these variants). Consider two origins






custom-character
a
=
custom-character
∪{v
a},custom-characterb=custom-character∪{vb},vacustom-character,vbcustom-character,va≠vb.  (38)


Let us first set up the following notation: denote C+ all origins in a population having variants C. Probability of the true origin being within the sets of origins custom-charactera+ and Cb+ might be considered equal, given the same number of matching variants






P(Γ(s)∈custom-charactera+|Φ(custom-charactera,s)={va})=P(Γ(s)∈custom-characterb+|Φ(custom-characterb,s)={vb})  (39)


Thus if the true origin is within C+a or C+b with the same probability, the size of both sets can be investigated. Suppose that the individual origins (within C+a and C+b) are equally likely a priori. Then the relative differences in the probability of the origins Ca, Cb are inverse to number of entities having the corresponding variants:











P


(


Γ


(
s
)


=


a


)



P


(


Γ


(
s
)


=


b


)



=


(





a
+







b
+




)


-
1






(
40
)







Furthermore, the relative difference between the number of the entities might be derived using population frequencies as:











(





a
+







b
+




)


-
1


=






v



a









ϕ


(
v
)







v



b









ϕ


(
v
)




.





(
41
)







The relative differences between the individual origins are preserved, and if the true origin is within the considered origins, they can be normalized to the value of one and thus (36) is readily established. Although, the assumption (39) was based on same number of matching variants, it might be assumed to hold, in general, for any sets of matching variants.


Another use of the method of the invention is in diagnosis of cancer by identifying somatic mutant peptides attributable to tumour in a sample taken from the body of a patient, e.g., blood or other fluids. The identification of the somatic mutant peptides attributable uniquely to tumour can be used for non-invasive diagnosis and monitoring of progression and recurrence of the disease.


For determination of the status of a variant (somatic or germline), various criteria can be used. Herein, global nucleotide alterations for the purpose are used.


Germline variants are considered as follows: a variant is present in dbSNP (v. 147), or ExAC (version of ExAC compilation without TCGA) and is preferably of a population frequency higher than 1·10−4 (in any of dbSNP or ExAC). Somatic variants are defined as those present in COSMIC, ICGC or TCGA, but not present in dbSNP and also not present in ExAC.


The presence of the somatic mutant protein variants (e.g., in blood of individual) can be in itself a sign of presence of a tumour. This is especially true for tumours with a high mutation rate, e.g., a melanoma.


For a more precise identification of the somatic mutant proteins exclusively attributable to tumour, the sample (e.g. of blood or blood plasma or blood serum or tears or urine or saliva or stool or breath condensate or lavages or effusions or liquors, etc.) of the patients may be analyzed before and after treatment (e.g., surgery, radiation, chemotherapy, biological therapy, immunotherapy, etc.). The drop in somatic mutant proteins after treatment establishes their exclusive correspondence to tumour and ultimately the tumour response. This can be done for establishing standards for such measurements or for the monitoring of a patient.


Another possible use of the method of the invention is in the monitoring the response of a recipient after transplantation by selection and quantification of peptides of the donor in a sample taken from the body of the recipient. Identification of increasing quantities of donor peptides in the sample (e.g. of blood or blood plasma or blood serum or tears or urine or saliva or stool or breath condensate or lavages or effusions or liquors, etc.) of the recipient is a sign of rejection or a risk of rejection of the transplanted organ.


Analysis of polymorphisms in the proteins is performed both on donor and recipient separately. Once these polymorphisms are associated with donors and recipients, identification of polymorphic peptides is uniquely linked to donor and recipient as well. Non-exclusive polymorphisms are not considered.


The quantification can be done using any label-free quantification method, for example by integration of area under curve in LC/MS spectra. For precise quantification, targeted quantitative methods such as SRM/MRM can be used. Once polymorphisms are established and transitions of polymorphic peptides available, their monitoring in individuals can be readily performed.


Furthermore, calculation of pair-wise variant-based agreement between individuals is performed. The method is schematically shown in FIG. 8; the correspondence 806 between entities 802 and 804 based on their variants 803 and 805 is determined. Two samples sa,sb from the tested individuals are analyzed and their agreement Φ(sa,sb) is determined using the method of the invention, optionally with a probabilistic interpretation. In the following description of the method, φ denotes function from variant to its population frequency (such function might be for instance derived from population frequencies in dbSNP database).


The agreement may be based plainly on number of matching variants identified using particular methods ma,mb, for example, as follows:





Φ(sa,sb)=|ma(sa)∩mb(sb)|  (42)


Alternatively, the agreement may be in probabilistic terms. For the purpose, Γ is a function from sample to its true origin, wherein the origin e is a subset of all variants (neglecting the possibility that two distinct origins have the same variants). The probability of two samples having the same origin, given the observed agreement, is:






Pa(sa)=Γb(sb|Φ(sa,sb)=x).  (43)


Yet alternatively, the probability of at least as extreme match as x at random may be used:






P(Φ(sa,sb)≥x).  (44)


Method m of identification of variants applied to a sample may identify exactly the variants in origin, and the origin is the same if and only if the identified variants are equal in both samples, however such situation is less likely in practice.


In some embodiments, the method m applied to the sample identifies a proportion r of variants in the sample. The proportion might be unknown in advance (or it might depend on the concentration of the sample, etc.), but the fact that samples are drawn from a known population may be utilized for its derivation. In such case, the expected number of variants in the sample with known population frequencies is





Σϕ(v)=n.  (45)


The probability of both the presence of variant v and its identification using method m is designated as Pm+(v).


In some embodiments, identification of a variant might be independent of the variant itself, and therefore the probability of the identification is equal for each variant. In other embodiments, the probabilities may be different. Nevertheless, if n variants were identified using the method m in a sample, then the probabilities of identification can be expressed as the expected number of identified variants being the actual number of identified variants:





ΣPm+(v)=n.  (46)


In shotgun proteomics, it is beneficial to model the probability of identification of a variant. Such identification might be modelled as a function of variant abundance because more abundant proteins are more likely to be measured. This is also important for establishment of agreement due to chance, because identification of variants of high population frequency in highly abundant proteins is quite likely even for random individuals. The probability of identification of a variant can be modelled on per-protein basis as its abundance. However, it is preferred to model it as a coverage C(p) of protein p (by identified reference peptides), which effectively normalizes the abundance within zero-one range. Proteins which are highly covered by peptides give a high probability of identification of a variant (if the variant is present) as opposed to proteins with a low coverage. For further simplification, the coverage might be calculated over genes and restricted to peptides uniquely alignable to genes (around 90%). The coverage for a gene might be then defined as the average of coverages of proteins (corresponding to the gene). This is followed by a further normalization of the probabilities of identification, such that (46) holds, as follows:












P

m
+




(
v
)


=




C


(
v
)


·

ϕ


(
v
)






v



(


C


(
v
)


·

ϕ


(
v
)



)



·
n


,




(
47
)







Depending on the agreement model, at least as good match (44) may be calculated using different approaches. In general, the probability may be numerically calculated using viable methods, for example, Monte Carlo simulation. The following paragraphs focus on number of matching variants (42).


In general, the probability of a match of particular variant v at random using methods ma,mb is






P
m

a
(vPmb(v),  (48)


if these events are statistically independent. The situation is, however, more complicated in shotgun proteomics, because some peptides are more suitable for identification (e.g. due to ionisation characteristics) than others. In other words, identification of variant in first sample changes, usually increases, the probability of identification of a variant in the second sample.


This effect can be modelled on pairwise comparison of coverage of matching proteins in both samples. Suppose Ca(p) is a coverage of protein p in sample a and Cb(p) in sample b. The expected shared coverage if distributed uniformly is Ca(p)·Cb(p). However, as there are preferences for individual peptides, the actual shared coverage is usually higher. The relationship could be modelled in variety of ways. Given the large set of available data, it can be modelled also using k-nearest neighbor regression. Here, the regression model is represented as a function k (5 neighbors, Euclidean distance). Therefore instead of multiplication of probabilities as in (48), it is calculated as






k(Pma+(v),Pmb+(v)).  (49)


If the agreement model is a number of matching variants (42), the calculation of (44) may be approximated, for example, using Binomial distribution, with the probability of success in one trial being equal to mean value of (49) over all variants.


The determination of the identity of the entity solves the problems of interpretation of mass spectra commonly encountered in shotgun proteomics and many other fields.


The method of the present invention may be also used for determination of identity, in particular for authentication of cell lines or identification of an individual from mass spectra of proteome.


The method may be also used for identification of non-host organism from mass spectra of proteome of host organism, in particular for diagnosis of microbial infection or colonization.


The method may be also used for identification of presence of a tumour from mass spectra of body fluid proteins or estimation of tumour characteristics through presence or absence of somatic mutations.


The method may be also used for monitoring of organ transplantation and early detection of transplant rejection from mass spectra of biological materials of recipient.


EXAMPLES OF CARRYING OUT THE INVENTION
Example 1—Determination of Identity of an Entity
Collecting Analytical Data

This example illustrates fragment mass spectrum of analytical data collected for unknown peptide in shotgun proteomics. This particular example of MS/MS spectrum is shown on FIG. 9 and the process of determination of entity is further illustrated on it.


Obtaining Candidate Identity by Enumeration

The candidate identities for spectrum on FIG. 9 are obtained through enumeration, whose description follows. The probabilities for cleavage after particular amino acid (even modified) are specified in FIG. 6b. Probabilities of few modifications were set as in FIG. 6c. The rest of modifications (but not substitutions), were set prior-like probability of 0.001. The prior-like probability of amino acid substitution was set to be dependent on the least number of nucleotide substitutions of given codon. If the substitution can happen (for any combination of codons) in one nucleotide alteration, it is 0.0002=q, otherwise its powers; thus if n is the minimum number of nucleotide substitutions, then the prior-like probability is nq. The prior-like probabilities of coded amino acids and terminals were set such that sum over their prior-like probability of amino acid and all its modifications equals one. Small partial list of modifications of amino-acids, along with their prior-like probabilities is illustrated in the following table:














modification
amino acids
prior-like








[C]
0.0000



[M]
0.6828



[K]
0.8186



[N-term]
0.9080



[S]
0.9156



[T]
0.9314



[Y]
0.9338



[N]
0.9496



[R]
0.9566



[H]
0.9626



[C-term]
0.9670



[D]
0.9686



[E]
0.9688



[W]
0.9800



[P]
0.9826



[F]
0.9848



[Q]
0.9858



[L]
0.9900



[G]
0.9914



[I]
0.9942



[V]
0.9954



[A]
0.9966


Methyl
[N-term, E, D, C-term, L, I,
0.0010



R, Q, N, K, H, C]



Methyl
[T, S]
0.0100


Oxidation
[W, H, C, R, Y, F, P, N, K,
0.0010



D, G]



Oxidation
[M]
0.3000


Phospho
[R, C, D, Y, H, T, S]
0.0010


. . .
. . .
. . .









In the next step, all candidate identities of prevalence higher than 4·10−6 were obtained, and only candidates within 5 ppm (parts-per-million) of calculated precursor mass were considered and are ordered from highest prevalence, an excerpt shown in the following table:














sequence
prevalence
ppm

















R.-FYPHDVLSL-.P
1.1 . 10-3
2.859482


(SEQ ID NO. 1)







S.-MVELWAWR-.E
8.8 . 10-4
−3.124660


(SEQ ID NO. 2)







R.-SPS(Methyl)GFGDPGKK-.D
8.5 . 10-4
−0.839305


(SEQ ID NO. 3)







R.-MTQALALQAGS-.L
8.0 . 10-4
2.272081


(SEQ ID NO. 4)







D.-NTVVMEEIR-.R
7.8 . 10-4
2.272081


(SEQ ID NO. 5)







T.-AISNASDVWK-.K
7.3 . 10-4
−0.830127


(SEQ ID NO. 6)







R.-(Formyl)ANAFLEELR-.P
5.6 . 10-4
−0.834716


(SEQ ID NO. 7)









The precursor mass difference of 5 ppm was selected due to accuracy of the employed mass spectrometer (Orbitrap Elite). Depending on experimental conditions, the precursor tolerance can be much wider (e.g., 500 Da) as is the case of open search or total (all candidate identities independent of precursor mass are considered). In these cases, the mass difference is further localized (or decomposed into multiple modifications and their localization) as is usual in open search, but prevalences of candidate identities with localized masses are further updated by corresponding prevalences of modifications.


The limiting prevalence of 4·10−6 should be low enough for most experiments at present. The number of spectra per experiment is in order of hundreds of thousands; in such case we would expect at most few peptides identified if lower minimal prevalence was taken into account.


Calculating Score
Agreement

This section describes agreement of theoretical spectrum of peptide and experimental (measured) spectrum. Number of matching peaks (of experimental and theoretical spectrum) is used as particular agreement model (FIG. 10). In this example, only singly charged ions (b, y) are used for prediction of theoretical spectrum. The agreement in FIG. 10 is shown for two peptides (from those enumerated in previous step), placed at top and bottom separately. Prefix (b) ions are shown closer to the MS/MS spectrum and suffix (y) ions are shown further (both on top and bottom). Ions matching experimental spectrum (fragment tolerance of 0.3 Da), are thicker. The agreement corresponds to total number of matching peaks. The agreement of individual peptides is illustrated in the following table, in which the first few peptides are ordered from highest spectral match.















sequence
agreement
prevalence
ppm


















K.-M(A→T) TQALEELR-.S (SEQ ID NO. 8)
12
5.7 . 10-5
2.272081





R.-(Carboxymethyl)TSGGAGGLGSLR-.A (SEQ ID NO. 9)
11
4.8 . 10-4
−4.529832





K.-(Methyl)M(Oxidation)AQALEELR-.S (SEQ ID NO. 10)
11
1.5 . 10-4
2.267492





R.-SL(R→M) MVDIAEGR-.K (SEQ ID NO. 11)
11
6.8 . 10-5
2.262903





K.-LV(K→M)MDSIAEGR-.I (SEQ ID NO. 12)
11
5.6 . 10-5
2.262903





R.-TMQALEIE(L→R) R-.Q (SEQ ID NO. 13)
10
6.9 . 10-5
2.272081





.-MAS(Acetyl) S(Methyl) LLAGER-.L (SEQ ID NO. 14)
10
4.5 . 10-6
2.267492









Maximal Probability

The following table illustrates determination of identity using maximal probability of candidate identities (Pmax column), calculated from agreement and prior-like probabilities.















sequence
agreement
prevalence
Pmax


















K.-M(A→T) TQALEELR-.S (SEQ ID NO. 8)
12
5.7 . 10-5
1





R.-(Carboxymethyl)TSGGAGGLGSLR-.A (SEQ ID NO. 9)
11
4.8 . 10-4
1





K.-(Methyl) M(Oxidation) AQALEELR-.S (SEQ ID NO. 10)
11
1.5 . 10-4
0.235





R.-SL(R→M) MVDIAEGR-.K (SEQ ID NO. 11)
11
6.8 . 10-5
0.098





K.-LV(K→M) MDSIAEGR-.I (SEQ ID NO. 12)
11
5.6 . 10-5
0.069





K.-DC(GlyGly)VLNATLK-.Q (SEQ ID NO. 15)
10
2.9 . 10-4
0.377









Probability

The following description explains the association of probability of correctness (P) using Bayes' Theorem. For this purpose, model of true agreement and random agreement is defined.


Agreement of true interpretations is modelled as follows. Agreement is evaluated on interpreted spectra from spectral database of X!Hunter, which are assumed to be true interpretations. The behaviour (FIG. 11) is shown only for doubly charged fragment mass spectra. In this example, the model is taken as an average behaviour over number of residues. This is meaningful, as the behaviour over number of residues is quite independent on the length of peptide.


Agreement of random interpretations for this spectrum is visualized on FIG. 12. In this example, random agreement is modelled in such way that the probability of random match decreases 20× with increasing number of matching peaks.


In derivation of prior probability from prior-like probabilities, the c in (23) equals 1.0 (this is assumption that true interpretation is within the candidates), which enables establishment of probability. The associated probability of candidate identities is shown in the following table:


















sequence
agreement
prevalence
Ptrue
Prandom
Pprior
P







K.-M(A→T) TQALEELR-.S
12
5.7 . 10-5
0.066353
8.871817e−09
0.002357
0.520089


(SEQ ID NO. 8)











R.-(Carboxymethyl)
11
4.8 . 10-4
0.082515
1.774363e−07
0.019913
0.273163


TSGGAGGLGSLR-.A (SEQ ID








NO. 9)











K.-(Methyl)
11
1.5 . 10-4
0.082515
1.774363e−07
0.006103
0.083718


M(Oxidation)AQALEELR-.S








(SEQ ID NO. 10)











R.-SL(R→M)MVDIAEGR-.K
11
6.8 . 10-5
0.082515
1.774363e−07
0.002823
0.038721


(SEQ ID NO. 11)











K.-LV(K→M) MDSIAEGR-.I
11
5.6 . 10-5
0.082515
1.774363e−07
0.002325
0.031888


(SEQ ID NO. 12)











K.-DC(GlyGly)VLNATLK-.Q
10
2.9 . 10-4
0.103277
3.548727e−06
0.012038
0.010334


(SEQ ID NO. 15)











K.-DCVLNATLK-.Q (SEQ ID
10
2.7 . 10-4
0.103277
3.548727e−06
0.011107
0.009535


NO. 16)









Additional/Supporting Criteria

The run corresponding to the analytical data was used for selection of true interpretation as explained earlier (tryptic reference human peptides, 5 ppm precursor mass difference, statistical significance of 0.1).


Precursor Mass Difference

The FIG. 13 shows the distribution of precursor mass differences for true interpretations.


Retention Time

The FIG. 14 shows the distribution of experimental time for particular predicted theoretical time and its neighbors (the theoretical time is significantly shifted). With the assumption of symmetric difference, the FIG. 15 shows selection of interpretations near tails of the distribution (<5%) and those in the center (>95%).


Isotopic Distribution Difference

Similarly as for precursor mass difference, FIG. 16 shows the distribution of differences between theoretical and experimental isotopic distribution.


Protein Evidence

In case of protein evidence, 8.129% of assumed true interpretations were without existence of other peptide from the same protein.


Combination of Additional/Supporting Criteria

The FIG. 17 shows combination of supporting evidence, split for unlikely (≤5%) and likely (≥95%) results. For instance, in case of likely results, it can be seen that as retention time is getting closer to the center of the distribution (p closer to one), the precursor mass difference can be higher to still obtain probability greater than 95%. Therefore the figure captures the numerical relationship between these supporting criteria and the resulting probability.


The ROC curve on FIG. 18 shows capabilities of removal of incorrect interpretations with the use of supporting evidence. The filtering is evaluated on interpretations of variant peptides (statistical spectral significance E-Value of 0.1 in X!Tandem). In the ROC curve, true interpretations are assumed to be those which have sequencing support (variant also found in sequencing). It is clear that supporting evidence can help in removal of incorrect interpretations. For instance, here around 50% of sequencing unsupported results is removed, while retaining around 90% of sequencing supported results.


The following table contains associated evidence from supporting criteria and their combined values.
























pre-
reten-







agree-

cursor
tion
protein
isotopic



sequence
prevalence
ment
Pmax
mass p
p
p
p
comb. p























K.-M(A→T)TQALEELR-.S (SEQ
5.7 . 10-5
12
1
0.219550
0.712
1.000000
0.389195
0.479031


ID NO. 8)













K.-










(Methyl)M(Oxidation)AQALEELR-
1.5 . 10-4
11
0.235
0.221310
0.552
1.000000
0.389195
0.427736


.S (SEQ ID NO. 10)













R.-SL(R→M) MVDIAEGR-.K
6.8 . 10-5
11
0.098
0.222536
0.116
1.000000
0.389195
0.196714


(SEQ ID NO. 11)













K.-LV(K→M)MDSIAEGR-.I (SEQ
5.6 . 10-5
11
0.069
0.222536
0.436
0.081207
0.389195
0.099081


ID NO. 12)













R.-










(Carboxymethyl)TSGGAGGLGSLR-
4.8 . 10-4
11
1
0.001813
0.604
1.000000
0.970927
0.053857


.A (SEQ ID NO. 9)













K.-M(Oxidation) AQ(Methyl)
1.3 . 10-4
10
0.086
0.221310
0.552
1.000000
0.389195
0.427736


ALEELR-.S (SEQ ID NO. 17)













R.-TS(G→D) DGAGGLGSLR-.A
8.5 . 10-5
10
0.049
0.001813
0.636
1.000000
0.970927
0.055138


(SEQ ID NO. 18)









The combined p can be used for removal of matches which are unlikely to be correct. In this case, selecting expected removal of 10% of correct results, the first interpretation (highest scoring from viewpoint of spectral match) is not removed.


Determining Identity
Maximal Probability

Determination of identity in this example is based on selecting interpretation which is of maximal agreement and has Pmax=1. Such interpretation can be at most one, and it is the best candidate for given agreement and prevalence modeling; it is the first one in the previous table.


Probability

Determination of identity in this example is based on selection of the interpretation which has higher probability than 0.5; such interpretation can be at most one and it is the most likely interpretation. In this example, it is again the first one and the identity determined is the same as in previous example using Pmax and highest agreement.


Example 2—Identification of Observed Mutant and Polymorphic Proteins

The system overview (FIG. 19) represents example 102 of incorporation of prevalence model (FIG. 1) in shotgun proteomics for identification of variant peptides. In this example, candidate identities are first scored in database search and further reevaluated with the use of prevalence to obtain maximal probability of their correctness.


In general, the identification system 1901 corresponds to 102.1 and the rejection system 1902 corresponds to system comprising prevalence model 102.2. The search database for X!Tandem is represented in form of variant peptide fasta file constructed by translation of variant mRNA, and excerpt from it looks as follows:











>ID-00000000 none



(SEQ ID NO. 19)



NEIPIR







>ID-00000001 none



(SEQ ID NO. 20)



AAVAAITQALVGR







>ID-00000002 none



(SEQ ID NO. 21)



SPPLPGDLGGPSK







>ID-00000003 none



(SEQ ID NO. 22)



LSAAQTNGGGSAGMEGIMNPYTALPTPQQLLAIEQSVYSSDPFR







>ID-00000004 none



(SEQ ID NO. 23)



NTEILTGSWSDQTYPEGTHAIYK






The deep database 1902.1 corresponds to prevalence model, was obtained through enumeration (FIGS. 6a, 6b and 6c) and was stored as a peptide database, along with prior-like probabilities. It is preferred to store the database and index it by precursor mass, because the interpretations will be loaded for given precursor mass range; an excerpt of such record is shown here:
















precursor

prior-




mass
position
like
protein
sequence







1322.659981
540
0.000004
ENSP00000364543
QPT(Methyl)HQGGPAQDR(Methyl:2H(3)13C(1))






(SEQ ID NO. 24)





1322.659984
539
0.000058
ENSP00000292377
GGGGSARY(iTRAQ4plex114)NQGR (SEQ ID






NO. 25)





1322.659984
539
0.000059
ENSP00000292377
(iTRAQ4plex114)GGGGSARYNQGR (SEQ ID






NO. 26)





1322.659984
512
0.000484
ENSP00000400142
(iTRAQ4plex114)GNPSGGGNVQHR (SEQ ID






NO. 27)


1322.659991
179
0.000479
ENSP00000211998
QQELTHQEHR(Methyl:2H(3)13C(1)) (SEQ ID






NO. 28)









Further, it is preferred to construct the database first for a wide range of masses (for example, 700-2500 Da) and further index peptides into smaller ranges (for example, 0.01 Da), to save computational time.


The rejection system 1902 is an example of incorporation of prevalence model for reevaluation of candidate identities (FIG. 2), corresponding to 203.B. Rejection designates the reevaluation of candidate identities, in which maximal probability of correctness of candidate identity is evaluated and used for rejection of candidates.


Obtaining Candidate Identities, Database Search in Peptide Database of Variants

The process is illustrated step by step on identification of variant peptides on samples measured on colorectal cancer cell line HCT116. The steps can be split into three phases: i) spectral match using database search, ii) assignment of additional information, iii) obtaining additional candidate identities. In the first step, the variant peptide database is searched using database search method, herein using X!Tandem. Matching of spectra and variant peptides gives initial results, an example is illustrated in the following table, ordered by most significant matches first (E-Value).















MZ
sequence
[XTandem]scoreT
E-Value


















945.122864
GSGDPSSSSSSGNPIVYLDVDANGKPLGR (SEQ ID
107.099998
1.000000e−15



NO. 29)







912.987976
IGLIGGTGLDDPEILEGR (SEQ ID NO. 30)
77.800003
1.000000e−15





1156.580444
VEGTEPTTAFNIFVGNLNFNK (SEQ ID NO. 31)
72.400002
2.700000e−15





978.494324
GEELGGGQDPVQLISGFPR (SEQ ID NO. 32)
98.500000
3.900000e−15





593.311707
HPLHVTYAGAAVDELGK (SEQ ID NO. 33)
77.900002
5.500000e−15





975.565063
VPIASQGLGPGSTVLLVVDK (SEQ ID NO. 34)
80.599998
8.800000e−15





1045.49340
(Gln-¿pyro-Glu)QSGGSSQAGAVTVSDVQELMR (SEQ
70.900002
7.900000e−14



ID NO. 35)







893.930054
DPGGLTAGSTDEPPMLTK (SEQ ID NO. 36)
75.400002
1.500000e−13





893.929443
DPGGLTAGSTDEPPMLTK (SEQ ID NO. 37)
83.699997
6.000000e−13





1174.115967
NFYGGNGIVGAQVPLGAGIALAC(Carbamidomethyl)K
57.900002
1.100000e−12



(SEQ ID NO. 38)







889.464172
HPLHVTYAGAAVDELGK (SEQ ID NO. 39)
77.500000
1.400000e−12





908.914551
TVSGTC(Carbamidomethyl)GPGQPASSSGGPGR (SEQ
79.199997
1.600000e−12



ID NO. 40)







1079.080811
NSNPVIAELSQAINSGTLLSK (SEQ ID NO. 41)
77.599998
1.700000e−12





971.964539
VWLDPNETNELANANSR (SEQ ID NO. 42)
73.500000
1.900000e−12









There are many highly significant matches of variant peptides, which however, does not mean they are correct interpretations.


In the second step, variant peptides are aligned to reference protein-coding sequences (ENSEMBL, human genome), their distance to reference genome is calculated and additional information attached. Only reference peptides which can be result of a single nucleotide variation are considered in this example (this is also because prevalence of such peptides is much higher and simplifies identification task). Furthermore, here, only peptides which can be aligned to one genomic location are considered (such decision has some benefits, for example, it is easier to establish peptide-derived nucleotide variation, which has further benefits of deriving population frequency, or calculating correspondence to nucleotide sequencing of matching sample). Excerpt of results of this procedure is illustrated by the following table:















E-Value
sequence
substitution
candidate SNVs







1.000000e−15
GSGDPSSSSSSGNPIVYLDVDANGKP
(L, I)
[(chr10, 81107442, 81107442, C, A)]



LGR (SEQ ID NO. 43)







2.700000e−15
VEGTEPTTAFNIFVGNLNFNK (SEQ
(L, I)
[(chr2, 232324977, 232324977, C, A)]



ID NO. 44)







3.900000e−15
GEELGGGQDPVQLISGFPR (SEQ ID
(L, I)
[(chr11, 62444347, 62444347, C, A)]



NO. 45)







8.800000e−15
VPIASQGLGPGSTVLLVVDK (SEQ ID
(L, I)
[(chr1, 949537, 949537, C, A)]



NO. 46)







1.500000e−13
DPGGLTAGSTDEPPMLTK (SEQ ID
(I, L)
[(chr16, 30900254, 30900254, A, C),



NO. 47)

(chr16, 30900254, 30900254, A, T)]





1.700000e−11
VEHLEEGPMLEQLSK (SEQ ID NO. 48)
(I, L)
[(chr12, 93881377, 93881377, A, C),





(chr12, 93881377, 93881377, A, T)]





1.800000e−11
IVVECVMNNVTCTR (SEQ ID NO. 49)
(L, I)
[(chr8, 82196200, 82196200, T, A)]





9.100000e−11
GVEAANVTGPGGVPVQGSK (SEQ ID
(A, V)
[(chr1, 43162316, 43162316, C, T)]



NO. 50)







1.000000e−10
DMLLEIEEQLAESR (SEQ ID NO. 51)
(L, I)
[(chr2, 238662071, 238662071, C, A)]





1.100000e−10
IVVECVMNNVTCTR (SEQ ID NO. 52)
(L, I)
[(chr8, 82196200, 82196200, T, A)]





1.200000e−10
VNNSSLIGLGYTQTIKPGIK (SEQ ID
(L, I)
[(chr5, 133309490, 133309490, C, A)]



NO. 53)







1.300000e−10
SPASDTYIVFGEVK (SEQ ID NO. 54)
(A, V)
[(chr12, 57107324, 57107324, C, T)]





1.400000e−10
IVVECVMNNVTCTR (SEQ ID NO. 55)
(L, I)
[(chr8, 82196200, 82196200, T, A)]





1.500000e−10
IVVECVMNNVTCTR (SEQ ID NO. 56)
(L, I)
[(chr8, 82196200, 82196200, T, A)]





1.700000e−10
IEAGEGETAVLNQLQEK (SEQ ID
(Q, E)
[(chr12, 93213151, 93213151, C, G)]



NO. 57)







2.000000e−10
DMDLTGCLESGGSEEPGGIGVGEK
(N, D)
[(chr11, 57076574, 57076574, A, G)]



(SEQ ID NO. 58)







2.200000e−10
LSYDQQQQQQQQQQQQQQAIQSR
(L, I)
[(chr17, 17697138, 17697138, C, A)]



(SEQ ID NO. 59)









In the third step, for each spectrum with variant peptide interpretation, all candidate identities (minimal prior-like probability of 4·10−6) are enumerated (as described previously).


Calculating Score—Maximal Probability

Maximal probability of interpretations is established. Results of this procedure are illustrated in the following table:



















sequencing



sequence
E-Value
substitution
support
Pmax







GEELGGGQDPVQLISGFPR (SEQ ID NO. 60)
3.900000e−15
(L, I)

0.000201





VPIASQGLGPGSTVLLVVDK (SEQ ID NO. 61)
8.800000e−15
(L, I)

0.000201





DPGGLTAGSTDEPPMLTK (SEQ ID NO. 62)
1.500000e−13
(I, L)

0.000199





VWLDPNETNELANANSR (SEQ ID NO. 63)
1.900000e−12
(I, L)

0.000199





ALDVGSGSGLLTAC(Carbamidomethyl)FAR
4.400000e−12
(I, L)

0.000199


(SEQ ID NO. 64)









NSGPIGTEM(Oxidation)NTGFSSEVK (SEQ ID
6.900000e−12
(L, I)

0.000201


NO. 65)









FLQEHGSDSFIAEHK (SEQ ID NO. 66)
9.000000e−09
(L, I)

0.000201





TAMNVNEIFMALAK (SEQ ID NO. 67)
9.800000e−09
(I, L)

0.000199





SC(Carbamidomethyl)QTALVEILDVIAR (SEQ ID
1.100000e−08
(V, A)
X
1.000000


NO. 68)









WSGPLSLQEVDEQPK (SEQ ID NO. 69)
1.200000e−08
(Q, K)
X
1.000000





NLTNPNTVIIIIGNK (SEQ ID NO. 70)
1.200000e−08
(L, I)

0.000201





YGYTHISAGELLR (SEQ ID NO. 71)
1.300000e−08
(L, I)

0.000201





TVNELQNLTATEVVVPR (SEQ ID NO. 72)
6.600000e−08
(A, T)

0.000204





TILDQHGQYPIWMNQR (SEQ ID NO. 73)
1.000000e−07
(L, I)

0.000201





AVFVDLEPTVIDDIR (SEQ ID NO. 74)
1.300000e−07
(E, D)

0.000197





AAELLANSLATAGDGLIELR (SEQ ID NO. 75)
1.700000e−07
(I, L)

0.000199









Further, column “support” contains cell-line matched sequencing support; showing whether such variant was observed in sequencing. It can be seen that many of statistically significant results are not supported by sequencing. However, they have also low Pmax and can be rejected. The only two sequencing supported results in the example have Pmax=1.


Determining Identity—Maximal Probability

The criterium for determination of identity used here was maximal spectral match and Pmax=1.


Results—Maximal Probability
Identified Variants

The method was used for identification of variants in human family members (FIG. 21). The following table contains numbers of identified variant peptides and their sequencing support (evaluated against exome sequencing), separately for each family member.


















family
germline
germline
exome



member
variants [a]
variants [b]
support





















father
38
  30 (78.9%)
100.0%



mother
55
  42 (76.4%)
 91.3%



daughter 1
57
  43 (75.4%)
 91.5%



daughter 2
61
  53 (86.9%)
 93.0%



daughter 3
64
  53 (82.8%)
 94.6%



son 1 (twin)
35
  24 (68.6%)
100.0%



son 2 (twin)
42
  36 (85.7%)
 94.7%



average
50.3
40.1 (79.8%)
 95.0%










Note that the exome sequencing of particular sample was not used in construction of the global database. The evaluation of sequencing support against exome sequencing is most meaningful for germline variants as those are always present in substantial proportion.


Comparison with Proteome from Translated Exome


The previous table also shows comparison of number of identified variants, if knowledge of exome sequencing was used to create proteome with all variants. In such case, germline variants were based on exome sequencing in a following way: variant was found in at least one parent and one kid. The results suggest that even if sequencing of sample is available, its benefits are limited as around 80% of germline variants are identified with the use of global nucleotide database (at around 95% of sequencing correspondence).


Incompleteness of Exome Sequencing Support

In some cases when sequencing support is evaluated, the results might contain some of correct, but sequencing unsupported interpretation. This is because some variants are likely not sequencing supported because of low sequencing coverage of surrounding area, as illustrated on FIG. 20. For that reason, in previous comparison areas with coverage lower than 10 reads were removed from comparison.


Example 3—Authentication of Cell Lines

This example shows utilization of claimed method for identification of cell line. The analysis is performed on publicly available data of NCI60 panel (Gholami et al. (2013) Cell Reports, 4(3): 609-620). Variants were identified as in previous example (the system architecture on FIG. 19). For establishment of genetic origin, only variants of high population frequency, (above 1% as specified in dbSNP) were considered; variants of this kind are likely germline variants, are easy to identify (not many more likely interpretations a priori, and statistical significance of E-Value <0.1 is often enough) and are suitable for identification of origin.


Matching to Database of Origins

The identification of origin was performed against NCI60 exome database (Shankavaram et al. (2009) BMC Genomics, 10(1):277) and it was assumed that the true origin is within considered origins (thus within NCI60 exome database).


The process of matching of database of origins is further illustrated on a specific sample (P0001751) from NCI60 proteomes. An excerpt of results of variant peptide identification restricted to polymorphic peptides is shown in the following table:
















population


sequence
candidate SNVs
frequency

















G[R→P]EVDVNLPK (SEQ ID NO. 76)
[(chr11, 62298073, 62298073,
0.01338



G, C)]






EIENIQTY[Q→R] (SEQ ID NO. 77)
[(chr4, 186379734,
0.2133



186379734, A, G)]






GAAGSGD[A→S]AAAAEWIR (SEQ ID NO. 78)
[(chr19, 460641, 460641, G,
0.02097



T)]






TILTLTGVSTL[G→R] (SEQ ID NO. 79)
[(chr16, 2263835, 2263835, G,
0.01058



A), (chr16, 2263835, 2263835,




G, C)]






IIIE[Y→D]LLEATR (SEQ ID NO. 80)
[(chr12, 120575548,
0.9994



120575548, T, G)]






GS[S→L]GDRPEASMTPDAK (SEQ ID NO. 81)
[(chr16, 89927150, 89927150,
0.03235



C, T)]






GLVEPVN[V→M](Ox)VDNGDGTHTVTYTPSQEGPYM(Ox)VSVK
[(chr3, 58118554, 58118554,
0.6082


(SEQ ID NO. 82)
G, A)]






VGEAGLLSV[D→N]C(Carb)SEAGPGALGLEAVSDSGTK (SEQ ID
[(chr3, 58109161, 58109161,
0.5461


NO. 83)
G, A)]






QVVSAVT[T→A]LVEAAER (SEQ ID NO. 84)
[(chr8, 144942902,
0.6853



144942902, A, G)]









Polymorphic peptides were used to calculate match with exome sequencing data and used for calculation of probability of correct determination of origin, with an excerpt of results illustrated in the following table:
















product of
number of



cell line
matching
matching
probability







PR:PC 3
3.447000e−13
24
9.999987e−01


CNS:SNB 19
5.679920e−07
17
6.068739e−07


LE:SR
1.853638e−06
19
1.859584e−07


ME:MALME 3M
2.686513e−06
18
1.283074e−07


RE:ACHN
4.735100e−06
17
7.279668e−08


CNS:U251
5.081215e−06
18
6.783801e−08


OV:OVCAR 8
6.737296e−06
14
5.116289e−08









The data indicates that the most likely cell line within considered origins is PR:PC 3. In the metadata of P0001751, it can be seen that the cell line is claimed as RE:SN12C, therefore it can be concluded that it is likely, that the cell line is mislabeled.


The procedure was further performed over all proteomes in the dataset; the results are visualized in the following table and point out two likely mislabeled cell lines:



















Best






Claimed
matching
Claimed:
Best:
Best:


sample
cell line
cell line
SNPs
SNPs
Error Prob.







P001891 
ME:UACC
ME:UACC
41
41
2.14 · 10 − 12



62
62





P003198 
LC:NCI H23
LC:NCI H23
26
26
2.52 · 10 − 12


P003381 
LC:NCI
LC:NCI
42
42
7.86 · 10 − 11



H460
H460





P003207 
RE:ACHN
RE:ACHN
32
32
4.41 · 10 − 8 


P0001568
CO:HCC
CO:HCC
35
35
2.63 · 10 − 7 



2998
2998





P0001751
RE:SN12C
PR:PC 3
15
24
1.28 · 10 − 6 


P003820 
CO:KM12
CO:SW 620
32
41
2.29 · 10 − 6 


P001897 
ME:UACC
ME:UACC
31
31
2.64 · 10 − 6 



257
257





P001561 
CNS:SNB 75
CNS:SNB 75
33
33
7.56 · 10 − 6 









Identification of mislabeled cell lines is of high importance, because it prevents propagation of incorrectly drawn conclusions and contributes to reproducibility of science outputs.


Example 4—Identification of a Person

This example shows utilization of the method for identification of a person. The analysis is performed on in-house data of family of particular structure (FIG. 21).


Matching to Database of Origins

The example is analogous to matching of cell lines. The database of origins corresponds to sequencing database of family members. The same methods are used for assignment.


Results


The results of identification are shown in the following table:


















Claimed
Best matching
Claimed:
Best:
Best:


sample
person
person
SNPs
SNPs
Error Prob.







FAMD3
daughter 3
daughter 3
81
81
1.89 · 10 − 8


FAMM
mother
mother
67
67
1.56 · 10 − 7


FAMD2
daughter 2
daughter 2
62
62
6.27 · 10 − 4


FAMF
father
father
41
41
2.01 · 10 − 2


FAMD1
daughter 1
daughter 1
64
64
7.38 · 10 − 2


FAMS1
son 1 (twin)
son 1 (twin)
35
35
0.25


FAMS2
son 2 (twin)
son 1 (twin)
53
55
0.40









The only misidentification happened in case of monozygous twin and was indicated by high probability of error (0.4).


Example 5—Presence of Genetic Relationship Between Two Individuals

This example illustrates the utilization of the identified variants for determination of genetic relationship. For this purpose, variants were identified in blood lymphocytes of family members (FIG. 21) in proteomics data (the architecture of variant identification was done as in FIG. 19, with determination of identity as maximal agreement among candidate identities and Pmax=1).


Calculating Significance of Match

The calculation of significance of match is further illustrated here. Similarly as in identification of genetic origin, only variants more prevalent than 1% (as population frequency in dbSNP) were used for calculation of match (these variants are more likely germline variants).


For calculation of pairwise match, it is beneficial to establish tabular structure, which helps in organization of data. Small subset of such tabular structure is visualized in the following table:






























k(Pma + (v),




Chr
Ref
Alt
Start
p+
gene
Pma + (v)
Pmb + (v)
Pmb + (v))
a:has
b:has







chr11
G
T
62292881
1.0000
ENSG00000124942
0.389718
0.349291
0.294993




chr12
T
G
53070120
0.8011
ENSG00000167768
0.347364
0.446011
0.222365




chr1
T
G
40533346
1.0000
ENSG00000131236
0.225182
0.267302
0.213775
X
X


chr1
T
G
40533314
1.0000
ENSG00000131236
0.225182
0.267302
0.213775




chr1
T
G
40533265
1.0000
ENSG00000131236
0.225182
0.267302
0.213775




chr1
T
G
40533319
1.0000
ENSG00000131236
0.225182
0.267302
0.213775




chr1
T
G
40533286
1.0000
ENSG00000131236
0.225182
0.267302
0.213775









In the table, each line corresponds to particular variant. The “p+” refers to population frequency of a variant from database, and Pma+(v) to multiplication of coverage of individual gene in sample a and population frequency. Analogously for b. The k(Pma+(v), Pmb+(v)) refers to probability of its identification in both samples.


The calculation of significance of match was approximated by Binomial distribution with mean value of






k(Pma+(v),Pmb+(v)).  (50)


The results of this method (the approximation by Binomial distribution) are further illustrated here:


















k
n
p *
p-value





















0
35433
0.000474
1.000000e+00



1
35433
0.000474
9.999999e−01



. . .
. . .
. . .
. . .



23
35433
0.000474
8.624331e−02



24
35433
  0.000474 *
5.665743e−02



25
35433
0.000474
3.596854e−02



26
35433
0.000474
2.208026e−02



. . .
. . .
. . .
. . .










Results

The method was further applied in pairwise manner to all family members (FIG. 21) and the heatmap (FIG. 23) shows results of at least as good match between them. The interest was thus in calculation of probability of at least as good match at random (43), here visualized as a color on the heatmap.


Example 6—Presence of a Tumour in Human

This example illustrates embodiment for identification of tumour-specific circulating proteins in blood serum. In this example, publicly available data accessible on PRIDE (identifiers: PXD004624, PXD004625, PXD004626) were used for identification of mutant proteins. For identification of variants, the same method which corresponds to FIG. 19 was used.


Determination of Variant Status
Selection of Tumour-Attributable Mutations

In this example, mutations attributable to tumour were assumed to be all somatic mutations identified.


Results

The results (FIG. 24) show presence of mutant peptides in melanoma cancer patients, with higher presence in advanced cachectic ones, lower presence in less advanced non-cachectic ones and almost none in controls. In this example, mutant peptides can be roughly associated with the presence of tumour and extent/stage of cancer.


Example 7—Identification of Grafted Peptides in Xenograft Model

In this example, human reference and variant proteins are identified in blood serum from murine xenografts. Configuration of the experiment is based on FIG. 19, with difference in enumeration of candidate identities explained further.


Obtaining Candidate Identities, Enumeration of Candidate Identities

In enumeration of candidates, peptides are enumerated for both organisms (here, mouse and human), limiting condition being prior-like probability of 4·10−6. Prior-like probabilities of peptides enumerated for human are multiplied (herein and in practice linearly scaled down) by relative difference of prevalence of human to mouse. The number is derived for particular experimental circumstances.


Prevalence of Peptides in Xenografts

The proportion in this example of p=0.25 was estimated from homologous peptides and heterologous protein evidence, explained earlier. Tabular structure corresponding to homologous peptides and their heterologous protein evidence is illustrated here:



















human
mouse


homologous sequence
human proteins
mouse proteins
evidence
evidence







YIETDPANR (SEQ ID NO. 85)
ENSP00000377882; . . . 
ENSMUSP00000028239; . . . 

X





VGDTYERPK (SEQ ID NO. 86)
ENSP00000346839; . . . 
ENSMUSP00000140816; . . . 

X





DIFTGLIGPMK (SEQ ID NO. 87)
ENSP00000420545; . . . 
ENSMUSP00000003714; . . . 

X





SGIPIVTSPYQIHFTK (SEQ ID
ENSP00000245907
ENSMUSP00000024988
X
X


NO. 88)









KLWAYLTINQLLAER (SEQ ID
ENSP00000351190; . . . 
ENSMUSP00000046530

X


NO. 89)









RFDEILEASDGIMVAR (SEQ ID
ENSP00000320171; . . . 
ENSMUSP00000128770; . . . 

X


NO. 90)









LLCNGDNDCGDQSDEANCRR
ENSP00000440113; . . . 
ENSMUSP00000066940; . . . 

X


(SEQ ID NO. 91)









IHGILSNTHR (SEQ ID NO. 92)
ENSP00000366460
ENSMUSP00000110350; . . . 

X





EVYGFNPEGK (SEQ ID NO. 93)
ENSP00000382200; . . . 
ENSMUSP00000127147; . . . 
X
X





VNLGVGAYR (SEQ ID NO. 94)
ENSP00000359539
ENSMUSP00000026196

X





VQFELHYQEVK (SEQ ID
ENSP00000351190; . . . 
ENSMUSPP0000046530; . . . 

X


NO. 95)









Determination of Identity

Peptides with Pmax=1 and maximal spectral agreement are retained.


Results

The identification method was used for identification of human protein biomarkers across wide range of cancer tissues transplanted to mice. The results (FIG. 25) show presence of human peptides and in general lack of such peptides in immunocompromised SCID mice, providing reliability of results.


Example 8—Identification of Non-Host Organisms

This example illustrates utilization of prevalence for diagnosis of mycoplasma in host organism. In this case, the prevalence of non-host organism is assumed to be unknown and thus refers to the more complicated situation as described earlier.


Obtaining Candidate Identities

For particular mass spectrum, peptides mapping exclusively to reference mycoplasmal peptides (among all organisms) and all human peptides (prior-like probability of 4·10−6) were obtained. As previously described, mycoplasmal peptides were defined to be of strictly lower prevalence than any enumerated human peptide.


Determining the Identity

Only mycoplasmal peptides of highest agreement and Pmax=1 were retained.


Results

This identification approach was applied to subset of projects in PRIDE repository and results of diagnostics of Mycoplasma are illustrated on FIG. 26 (SEQ ID NOs. 96-131).


Example 9—Use of Isotopic Labels for Identification of Mutations

The following example illustrates utility of presence of both light and heavy isotopic forms using Stable isotope labeling with amino acids in cell culture (SILAC) of sample for identification of variants. The variant identification was performed as in previous examples (FIG. 19).


Determining Identity

The additional criterium in this case is identification of both light and heavy forms of peptide of interest.


Results

The analysis of SILAC-pairing and its effect on sequencing support of variants was analyzed for two criteria (initial significance: E-Value ≤0.1) and after rejection (Pmax=1). The results in the following table show that peptides which were identified in pairs have much higher sequencing support (when evaluated against sequencing of matching sample).
















Sequencing support












Initial
After rejection




(E-Value <0.1)
(Pmax = 1)














SILAC

SILAC





Paired
Single
Paired
Single















nonSNP
HCT116
 2.3%
0.8%
 45.5%
 36.8%



CCRF-CEM
 3.2%
0.7%
 50.0%
 29.4%


SNP
HCT116
60.0%
3.5%
 96.0%
100.0%



CCRF-CEM
64.3%
2.9%
100.0%
 82.4%









This difference was most remarkable in variants of low population frequency (more likely somatic mutations). Therefore, isotopic labels can be utilized to increase specificity of identification of somatic mutations.


Example 10—Identification of Splicing Variants
Obtaining Candidate Identities, Enumeration of Candidate Identities

The scheme (FIG. 27) refers to enumeration, in which alternatively spliced proteins (and their prevalence) are constructed from reference exon-based protein models. This schema is in direct, unit-wise correspondence with generic enumeration (FIG. 5).


The enumeration for particular protein starts with reference exon-based protein model 2701, for which individual exons of corresponding gene are either present or not. Such a model can be represented by a binary vector representing presence of exon in the model. Different protein models 2702 are constructed by exon inclusion or exon exclusion events 2703 with associated effect on prevalence.


In some embodiments, prevalence might be expressed in prior-like probabilities and then exon inclusion, or exon exclusion are assigned probabilities of these events.


The enumeration process continues until limiting minimal prevalence condition 2704 is met. The protein models are transformed by concatenation of individual exons and translated into proteins 2705 with their associated prevalence, which then further constitutes the prevalence model 2706. Proteins constructed in this way might be used, for example, directly in top-down proteomics in identification, or proteins might be further digested for use in bottom-up proteomics.


Example 11—Identification of Tumour, Protein Variants, Correlation with Clinical Characteristics

This example illustrates embodiment for identification of tumor size and associated disease stage. Herein, publicly available data from Clinical Proteomic Tumor Analysis Consortium, specifically TCGA Colorectal Cancer were used for identification of mutant proteins using the same method which corresponds to FIG. 19.


Results

The parameters of determination of somatic and germline variants are as follows. Germline variants are considered as follows: a variant is present in dbSNP (v. 147), or ExAC (version of ExAC compilation without TCGA) and is preferably of a population frequency higher than 1·10−4 (in any of dbSNP or ExAC)


The results on FIGS. 28a, 28b and 28c show the behavior of identified variants. The proportion of identified somatic mutant peptides among all reference peptides is visualized on FIG. 28a and shows clear increase in proportion of mutations with increasing tumour stage. Therefore, given particular reference measurement system, the increase in somatic mutations shows strong correlation to the tumour stage. Similar, but more pronounced effect can be seen when derived using nucleotide sequencing (FIG. 28b). Finally, the proportion of germline variants derived using proteomics does not show association with the tumour stage, showing that it is the effect of somatic mutations which is increased due to the higher tumour heterogeneity in more advanced stages.

Claims
  • 1: A method for determination of identity of at least one entity from a mass spectrum of said at least one entity and optionally from additional data from chemical, physical, biochemical or biological analysis of said at least one entity, for each entity comprising the steps of: a) collecting analytical data from mass spectrum of the entity, and optionally collecting additional analytical data from a chemical, physical, biochemical or biological analysis of the entity,b) obtaining a plurality of candidate identities of the entity and obtaining the prevalences of said candidate identities of the entity, whereas for each candidate identity it applies that all candidate identities with a higher prevalence are included in the plurality of candidate identities;c) for each candidate identity of an entity, calculation of its score, said calculation involving at least prevalence of entity, or at least prevalence of entity and agreement with mass spectrum,d) determining the identity of an entity as the candidate identity with the score closest to the score which would correspond to the true identity of the entity.
  • 2: The method according to claim 1, wherein in the step c), the calculation involves calculating maximal probability of candidate identity or calculating probability of candidate identity, optionally using Bayes' Theorem.
  • 3: The method according to claim 1, wherein in the step b) the value of prevalence is calculated based on at least one of population frequency of said entity, probability of modification of said entity in the environment, probability of modification of said entity during the analysis step.
  • 4: The method according to claim 1, wherein in the steps b) and c), the value of prevalence is expressed as prior probability or as prior-like probability.
  • 5: The method according to claim 1, wherein the entity is selected from a peptide, a protein, a lipid, a nucleic acid, a metabolite and a molecule having the molecular weight of up to 2000 mol/g.
  • 6: The method according to claim 1, wherein in the step b), the obtaining of the candidate entities and/or of prevalences of the candidate identities comprises enumeration which comprises the steps of: b.a) selecting initial candidate identities with initial prevalences;b.b) transferring said initial candidate identities into a base of candidate identities;b.c) producing new candidate identities by application of events to said base of candidate identities, and incorporating said new candidate identities into said base of candidate identity identities and continuing said producing unless a limiting condition is met;b.d) transforming the base of candidate identities obtained in step b.c) into candidate identities with associated prevalence.
  • 7: The method according to claim 6, wherein said candidate identities are peptides; said prevalence is expressed as prior-like probability; said initial entities are N-terminally-cleaved linear subsequences of reference proteins; said applicable events comprise modification, substitution and cleavage; said limiting condition is minimal prior-like probability of given form of peptide; or wherein said candidate identities are proteins; said prevalence is expressed as prior-like probability; said initial entities are reference exon-based protein models; said applicable events comprise exon exclusion and exon inclusion; said limiting condition is minimal prior-like probability of exon-based model; said transformation of entities into hypotheses is concatenation of exons into protein-coding sequence and translation in silico.
  • 8: The method according to claim 1, wherein the entities are proteins, wherein the step of obtaining the candidate identities of an entity in step b) includes database search in database of peptide variants, and wherein the method is used for identification of mutant and polymorphic proteins from mass spectra of proteome, with alterations already observed globally on nucleotide level
  • 9: The method according to claim 1, wherein the entities are peptides, further comprising the steps of: e) matching of entities determined as polymorphic or germline peptides to database of origins,and wherein the method is used for authentication of cell lines or identification of a person from mass spectra of proteome.
  • 10: The method according to claim 1, wherein the entities are non-host peptides, wherein in the step b) the prevalence is expressed as prior or prior-like probability and prevalence of non-host peptides is scaled down according to prevalence of non-host organism, and wherein the method is used for identification of non-host organism of known prevalence from mass spectra of proteome of host organism, for instance to identify microbial infection or colonization of the host.
  • 11: The method according to claim 6, wherein the entities are non-host peptides, wherein in the step b) in obtaining the candidate identities, peptides uniquely mapping to non-host organism are added to enumerated peptides of host organism and prevalence of non-host peptides is lower than of any host peptide, and wherein the method is used for identification of non-host organism of unknown prevalence from mass spectra of proteome of host organism.
  • 12: The method according to claim 1, wherein the entities are donor peptides, wherein in the step b) the prevalence of donor peptides is scaled according to their prevalence among recipient peptides, and wherein the method is used for identification of proteins originating from grafted tissue in recipient.
  • 13: The method according to claim 1, wherein the entities are peptides, the method further comprising the step of: e) selecting somatic mutant peptides attributable to tumour,wherein the method is used for identification of presence of a tumour from mass spectra of circulating proteins or estimation of tumour biological characteristics through increase in number of somatic mutations.
  • 14: The method according to claim 1, wherein the entities are peptides, the method further comprising the step of: e) selection and quantification of polymorphic peptides attributable to donor,wherein the method is used for monitoring tissue or organ grafting and early detection of transplant rejection from mass spectra of biological material of recipient.
  • 15: The method according to claim 1, wherein the entities are peptides, said method further comprising the steps of: e) calculating significance of match between two individuals based on polymorphic peptides,wherein the method is used for determination of presence of genetic relationship between two or more individuals from measured mass spectra of proteome.
Priority Claims (1)
Number Date Country Kind
18184710.4 Jul 2018 EP regional
PCT Information
Filing Document Filing Date Country Kind
PCT/EP2019/069552 7/19/2019 WO 00