(en)Methods are provided for determining the expression level of target proteins in a subject. A plurality of respective peptide transitions are generated from a plurality of proteins obtained from a biological sample from the subject, wherein the plurality of proteins comprises both target and normalizing proteins. A mass spectroscopy (MS) signal intensity is measured from the plurality of respective peptide transitions and a plurality of corresponding stable isotope-labeled internal standard (SIS) peptide transitions. For each of the plurality of proteins, a response ratio is calculated between the MS signal intensity of the respective peptide transition and the corresponding SIS peptide transition. The response ratio for each target protein is normalized by a sample-dependent normalization factor calculated from the response ratio for each normalizing protein, wherein the normalized response ratios provide a determination of the expression level of the target proteins.

1.ApplicationNumber: US-201514612959-A
1.PublishNumber: US-9594085-B2
2.Date Publish: 20170314
3.Inventor: LI XIAO-JUN
HUNSUCKER STEPHEN W.
HAYWARD CLIVE
KEARNEY PAUL EDWARD
LEE LIK WEE

4.Inventor Harmonized: LI XIAO-JUN(US)
HUNSUCKER STEPHEN W(US)
HAYWARD CLIVE(US)
KEARNEY PAUL EDWARD(CA)
LEE LIK WEE(US)

5.Country: US
6.Claims:
(en)Methods are provided for determining the expression level of target proteins in a subject. A plurality of respective peptide transitions are generated from a plurality of proteins obtained from a biological sample from the subject, wherein the plurality of proteins comprises both target and normalizing proteins. A mass spectroscopy (MS) signal intensity is measured from the plurality of respective peptide transitions and a plurality of corresponding stable isotope-labeled internal standard (SIS) peptide transitions. For each of the plurality of proteins, a response ratio is calculated between the MS signal intensity of the respective peptide transition and the corresponding SIS peptide transition. The response ratio for each target protein is normalized by a sample-dependent normalization factor calculated from the response ratio for each normalizing protein, wherein the normalized response ratios provide a determination of the expression level of the target proteins.

7.Description:
(en)RELATED APPLICATIONS
This application claims the benefit of, and priority to, U.S. Provisional Application No. 61/935,061 filed Feb. 3, 2014, the content of which is incorporated herein by reference in its entirety.


INCORPORATION-BY-REFERENCE OF SEQUENCE LISTING
The contents of the text file named “IDIA_012_001US_ST25.txt”, which was created on Feb. 20, 2015 and is 211 KB in size, are hereby incorporated by reference in their entireties.
BACKGROUND
Multiple reaction monitoring (MRM), also known as selected reaction monitoring (SRM) mass spectrometry (MS), allows for fast and reproducible measurement of tens to hundreds of proteins in complex biological samples such as bio-fluids, tissues, and cultured cells. There is tremendous interest in applying the technology to develop blood-based clinical tests for the diagnosis, prognosis or treatment stratification of various diseases. Due to the high complexity of the human blood proteome, proteomic analysis of blood samples typically consists of multiple experimental steps and is prone to variation ( FIG. 1A ). In addition, changes in laboratory conditions (e.g., operators, instruments, reagents) are expected during routine laboratory operations in clinical testing. Therefore, controlling analytical variability to satisfy rigorous quality control requirements for blood-based clinical testing using MRM-MS platforms has been challenging.
The principle of stable isotope labeling (SIL) is currently used in MS-based quantitative proteomics to control experimental variability. Protein abundance is measured by comparing MS signal intensities of endogenous peptides with those of their corresponding stable isotope-labeled internal standard (SIS) peptides. Three SIL approaches are potentially suitable for clinical testing ( FIG. 1B ). The first approach utilizes SIS peptides for protein quantification (SISQuan) and is the simplest one for implementation. SIS peptides are synthesized, optimized for MS analysis and spiked into samples before or after protein digestion to control variation in post-digestion procedures. However, variation occurring before or during digestion is not controlled. The second approach spikes full-length SIS proteins into samples before any analytical procedure takes place. While this approach offers the best control of analytical variability, it is applicable only to soluble proteins. Quality control of the production, the storage, etc., of SIS proteins as standards is challenging for routine laboratory operations. The third approach spikes either artificial or truncated SIS proteins into samples before protein digestion. It controls most variation in protein digestion and variation in subsequent procedures. However, it cannot control variation occurring before digestion and faces similar implementation challenges as the second approach. None of the above SIL approaches can control pre-analytical variability associated with sample collection and handling.
Thus, there is a need for a simple and robust method that provides sufficient control of pre-analytical and analytical variability for routine clinical testing on MS-based proteomics platforms. The present invention addresses that need.
SUMMARY
In one embodiment, the present invention provides a method for determining the expression level of target proteins in a subject. A plurality of respective peptide transitions are generated from a plurality of proteins obtained from a biological sample from the subject, wherein the plurality of proteins comprises both target and normalizing proteins. A mass spectroscopy (MS) signal intensity is measured from the plurality of respective peptide transitions and a plurality of corresponding stable isotope-labeled internal standard (SIS) peptide transitions. For each of the plurality of proteins, a response ratio is calculated between the MS signal intensity of the respective peptide transition and the corresponding SIS peptide transition. The response ratio for each target protein is normalized by a sample-dependent normalization factor calculated from the response ratio for each normalizing protein, wherein the normalized response ratios provide a determination of the expression level of the target proteins.
According to another embodiment, the determination of the expression level of the target proteins provides a diagnosis of lung disease for the subject.
According to yet another embodiment, the normalizing proteins are selected based on their ability to reduce intensity drift (D) of each of the plurality of respective peptide transitions, wherein intensity drift evaluates the deviation in abundance of each peptide transition from the overall median abundance of each peptide transition.
According to one embodiment, the normalizing proteins are selected based on their ability to reduce the median technical coefficient of variation (CV) of the plurality of proteins.
According to another embodiment, the plurality of proteins comprise at least two normalizing proteins selected from the group consisting of PEDF (Pigment epithelium-derived factor), MASP1 (Mannan-binding lectin serine protease 1), GELS (Gelsolin), LUM (Lumican), C163A (Scavenger receptor cysteine-rich type 1 protein M130), and PTPRJ (Receptor-type tyrosine-protein phosphatase eta).
According to yet another embodiment, the plurality of proteins comprise six normalizing proteins including PEDF_HUMAN, MASP1_HUMAN, GELS_HUMAN, LUM_HUMAN, C163A_HUMAN, and PTPRJ_HUMAN.
According to one embodiment, the plurality of respective peptide transitions comprise LQSLFDSPDFSK (SEQ ID NO: 25) (692.34, 593.30), TGVITSPDFPNPYPK (SEQ ID NO: 26) (816.92, 258.10), TASDFITK (SEQ ID NO: 27)(441.73, 710.40), SLEDLQLTHNK (SEQ ID NO: 28) (433.23, 499.30), INPASLDK (SEQ ID NO: 29)(429.24, 630.30), and VITEPIPVSDLR (SEQ ID NO: 30)(669.89, 896.50).
According to another embodiment, the target proteins comprise at least five of KIT_HUMAN, FRIL_HUMAN, COIA1_HUMAN, PRDX1_HUMAN, TENX_HUMAN, ENPL_HUMAN, GRP78_HUMAN, BGH3_HUMAN, ALDOA_HUMAN, GGH_HUMAN, CD14_HUMAN, LG3BP_HUMAN, TSP1_HUMAN, IBP3_HUMAN, TETN_HUMAN, and ISLR_HUMAN.
According to yet another embodiment, the target proteins comprise ALDOA_HUMAN, FRIL_HUMAN, COIA1_HUMAN, LG3BP_HUMAN, and TSP1_HUMAN.
According to yet another embodiment, the target proteins comprise ALDOA_HUMAN, FRIL_HUMAN, KIT_HUMAN, GGH_HUMAN, and TSP1_HUMAN.
According to one embodiment, the target proteins comprise KIT_HUMAN, FRIL_HUMAN, COIA1_HUMAN, PRDX_HUMAN, TENX_HUMAN, ENPL_HUMAN, TENX_HUMAN, ENPL_HUMAN, GRP78_HUMAN, BGH3_HUMAN, ALDOA_HUMAN, GGH_HUMAN, CD14_HUMAN, LG3BP_HUMAN, IBP3_HUMAN, TETN_HUMAN, and ISLR_HUMAN.
According to another embodiment, the biological sample is selected from the group consisting of tissue, blood, plasma, serum, whole blood, urine, saliva, genital secretion, cerebrospinal fluid, sweat, excreta, or bronchoalveolar lavage.
According to yet another embodiment, the plurality of proteins are obtained by immunoaffinity depletion.
According to one embodiment, the measuring step is performed by selected reaction monitoring mass spectrometry (SRM-MS).
According to another embodiment, the plurality of respective peptide transitions are generated by enzymatically digesting the plurality of proteins.
The methods of the current invention may include generation of the respective peptide transitions by fragmentation of the plurality of proteins by trypsin digestion. The methods of the current invention can include various manners to assess the expression of the plurality of proteins, including mass spectrometry (MS), liquid chromatography-selected reaction monitoring/mass spectrometry (LC-SRM-MS), reverse transcriptase-polymerase chain reaction (RT-PCR), microarray, serial analysis of gene expression (SAGE), gene expression analysis by massively parallel signature sequencing (MPSS), immunoassays, immunohistochemistry (IHC), transcriptomics, or proteomics.
According to yet another embodiment, the response ratio (R) is defined as:

R p,s =A p,s /Â p,s

where A p,s is the peak area of the signal intensity of peptide transition p in the biological sample s, and  p,s is the peak area of the signal intensity the corresponding SIS peptide transition.

According to one embodiment, the sample-dependent normalization factor (S) is defined as:
S s I = median ⁡ ( R 1 , s R ⋓ 1 , R 2 , s R ⋓ 2 , … ⁢ , R N , s R ⋓ N )
where S s I is the sample-dependent normalization factor calculated from the response ratios (R) of the N peptide transitions from normalizing proteins in the biological sample s, where R n,s is response ratio of peptide transitions from normalizing proteins n in the sample and {hacek over (R)} n is a scaling constant for the peptide normalizer that ensures values of {R n,s /{hacek over (R)} n } among all peptide transitions from normalizing proteins to be same on average.

According to another embodiment, the normalized response ratio ({tilde over (R)}) is defined as:

{tilde over (R)} p,s =R p,s /S s I ,

where p is peptide transition, s is the biological sample, R p,s is the response ratio of peptide p in biological sample s, and S s I is the sample-dependent normalization factor.

According to yet another embodiment, the intensity drift (D) is defined by:

D p,s =( I p,s −{hacek over (I)} p )/ {hacek over (I)} p

where p is peptide transition, s is the biological sample, I p,s is the abundance of peptide transition p in the sample s, and {hacek over (I)} p is the corresponding median value in all technical replica.

According to one embodiment, the plurality of respective peptide transitions and the plurality of corresponding stable isotope-labeled internal standard (SIS) peptide transitions are mixed together within a sample before the measuring step.
According to another embodiment, the determination of the expression level of the target proteins is independent of the volume of the sample.
In one embodiment, the subject has a pulmonary nodule, wherein the pulmonary nodule has a diameter of 30 mm or less. Preferably, the pulmonary nodule has a diameter of about 8 and 30 mm.
In one embodiment, the subject has a lung condition. The lung condition may be cancer or a non-cancerous lung condition. In another embodiment, the lung cancer is non-small cell lung cancer. The non-cancerous lung conditions include chronic obstructive pulmonary disease, hamartoma, fibroma, neurofibroma, granuloma, sarcoidosis, bacterial infection or fungal infection.
The subject can be a mammal. Preferably, the subject is a human.
One embodiment of the current invention includes assessing the expression of a plurality of proteins by liquid chromatography-selected reaction monitoring/mass spectrometry (LC-SRM-MS). In another aspect of the invention, at least one transition for each peptide is determined by liquid chromatography-selected reaction monitoring/mass spectrometry (LC-SRM-MS).
Although methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present invention, suitable methods and materials are described below. All publications, patent applications, patents, and other references mentioned herein are incorporated by reference in their entirety. The references cited herein are not admitted to be prior art to the claimed invention. In the case of conflict, the present specification, including definitions, will control. In addition, the materials, methods, and examples are illustrative only and are not intended to be limiting. Other features and advantages of the invention will be apparent from the following detailed description and claim.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a chart providing an overview of experimental variations and control methods. Panel (A) describes exemplar sources of variations. Systematic variations affect all proteins or peptides similarly. Protein- or peptide-specific variations affect only particular proteins or peptides. Random variations are not listed because they are not controllable. Panel (B) describes control of variations by different quantification methods in the analysis of plasma samples on a depletion-MRM-MS platform. Panel (C) provides an overview of three assessment studies and a combined HPS dataset.
FIG. 2 is a set of chromatograms showing the transitions of both endogenous and SIS peptides of individual proteins. Proteins are plotted in the same order as in Table 1. Transitions of endogenous peptides are plotted in the first, third, and fifth rows. Transitions of the corresponding SIS peptides are plotted in the second, fourth and sixth rows.
FIG. 3 is a set of graphs showing intensity drift of 16 target proteins and the six normalizing proteins as measured on 10 HPS samples in Study I described in the Examples. Results are shown in different panels for (A) raw MS data, (B) EPN, (C) SISQuan, and (D) InteQuan. The insert in panels (A)-(D) shows the mean (solid line) and 95% confidence internal (shaded band) of the protein drifts.
FIG. 4 is a graph showing receiver operating characteristic curves of the panel of all 16 target proteins, evaluated with Monte Carlo cross validation (MCCV) on clinical samples in Study I described in the Examples. Proteins are quantified using InteQuan (solid line) or SISQuan (dotted line).
FIG. 5 is a pair of charts showing principal variance component analysis (PVCA) of protein abundance in 15 clinical samples having three replicate measurements as described in the Examples. Protein abundance is shown in using (A) InteQuan and (B) SISQuan.
FIG. 6 is a set of graphs showing examples of SISQuan abundance versus loading volume in Study III. The top row shows two outlier measurements while the bottom row shows the corresponding linear behavior of the same proteins. Panel (A) shows PRDX1 in sample S3. The outlier measurement (Batch 2 with the loading volume of 15 μl) was removed from analysis. Panel (B) shows FRIL in sample S1. The outlier measurement (Batch 2 with the loading volume of 30 μl) was removed from analysis. Panel (C) shows PRDX1 in sample S2. Panel (D) shows FRIL in sample S6.
FIG. 7 is a set of graphs showing the calculation of generalized coefficient of variation (CV). Panels (A-D) show the results of all six clinical samples in Study III described in the Examples. Panels (E-H) show the results of the 29 HPS samples across all three studies described in the Examples. Panels (A, E) show the average InteQuan abundance versus experimental InteQuan abundance of individual proteins in individual samples. Panels (B, F) show the fitted SISQuan abundance versus experimental SISQuan abundance of individual proteins in individual samples. Panels (C, G) show the standard CV versus the generalized CV of InteQuan abundance. Panels (D, H) show the standard CV versus the generalized CV of SISQuan abundance.
FIG. 8 is a pair of graphs showing an assessment of alternative normalization panels for InteQuan. Panel (A) shows the median coefficient of variation (CV) of the target proteins after normalization by all possible combination of six normalizing proteins. Panel (B) shows the CV of the target proteins as evaluated using the “optimal” panel versus using the full panel.



DETAILED DESCRIPTION
The disclosed invention derives from the surprising discovery that endogenous protein normalization (EPN) and SISQuan may act synergistically in the control of different types of variations present in clinical proteomics. By selecting endogenous normalizing proteins for their ability to reduce both technical drift and technical variations of other proteins, rather than for their “housekeeping” properties, EPN may be combined with SISQuan in a unique fashion to provide numerous benefits and advantages over prior methods.
In summary, the quantification methods of the present invention demonstrate superiority to prior quantification methods such as SISQuan. The methods of the present invention are robust, simple to implement, capable of reducing pre-analytical and analytical variability, and able to improve the measurement of biological differences. All of these features make the methods of the present invention an ideal technique for MS-based quantitative proteomics, e.g., for applications in biomarker research and in routine clinical testing.
A. Definitions
The term “pulmonary nodules” (PNs) as used herein refers to lung lesions that can be visualized by radiographic techniques. A pulmonary nodule is any nodules less than or equal to three centimeters in diameter. In one example a pulmonary nodule has a diameter of about 0.8 cm to 3 cm.
The term “masses” or “pulmonary masses” as used herein refers to lung nodules that are greater than three centimeters maximal diameter.
The term “blood biopsy” as used herein refers to a diagnostic study of the blood to determine whether a patient presenting with a nodule has a condition that may be classified as either benign or malignant.
The term “acceptance criteria” as used herein refers to the set of criteria to which an assay, test, diagnostic or product should conform to be considered acceptable for its intended use. As used herein, acceptance criteria are a list of tests, references to analytical procedures, and appropriate measures, which are defined for an assay or product that will be used in a diagnostic. For example, the acceptance criteria for the classifier refer to a set of predetermined ranges of coefficients.
The term “partial AUC factor or pAUC factor” as used herein is greater than expected by random prediction. At specificity=0.80 the pAUC factor is the trapezoidal area under the ROC curve from 0.0 to 0.2 False Positive Rate/(0.2*0.2/2).
The term “robust” as used herein refers to a test or procedure that is not seriously disturbed by violations of the assumptions on which it is based. For the present invention, a robust test is a test wherein the proteins or transitions of the mass spectrometry chromatograms have been manually reviewed and are “generally” free of interfering signals.
The term “coefficients” as used herein refers to the weight assigned to each protein used to in the logistic regression model to score a sample.
In certain embodiments of the invention, it is contemplated that in terms of the logistic regression model of MC CV, the model coefficient and the coefficient of variation (CV) of each protein's model coefficient may increase or decrease, dependent upon the method (or model) of measurement of the protein classifier. For each of the listed proteins in the panels, there is about, at least, at least about, or at most about a 2-, 3-, 4-, 5-, 6-, 7-, 8-, 9-, or 10-, -fold or any range derivable therein for each of the coefficient and CV. Alternatively, it is contemplated that quantitative embodiments of the invention may be discussed in terms of as about, at least, at least about, or at most about 10, 20, 30, 40, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99% or more, or any range derivable therein.
The term “classifying” as used herein with regard to a lung condition refers to the act of compiling and analyzing expression data for using statistical techniques to provide a classification to aid in diagnosis of a lung condition, particularly lung cancer.
The term “classifier” as used herein refers to an algorithm that discriminates between disease states with a predetermined level of statistical significance. A two-class classifier is an algorithm that uses data points from measurements from a sample and classifies the data into one of two groups. In certain embodiments, the data used in the classifier is the relative expression of proteins in a biological sample. Protein expression levels in a subject can be compared to levels in patients previously diagnosed as disease free or with a specified condition. Table 5 lists representative rule-in classifiers (e.g., panels 1, 4, and 5).
The “classifier” maximizes the probability of distinguishing a randomly selected cancer sample from a randomly selected benign sample, i.e., the AUC of ROC curve.
In addition to the classifier's constituent proteins with differential expression, it may also include proteins with minimal or no biologic variation to enable assessment of variability, or the lack thereof, within or between clinical specimens; these proteins may be termed endogenous proteins and serve as internal controls for the other classifier proteins.
The term “normalization” or “normalizer” as used herein refers to the expression of a differential value in terms of a standard value to adjust for effects which arise from technical variation due to sample handling, sample preparation and mass spectrometry measurement rather than biological variation of protein concentration in a sample. For example, when measuring the expression of a differentially expressed protein, the absolute value for the expression of the protein can be expressed in terms of an absolute value for the expression of a standard protein that is substantially constant in expression. This prevents the technical variation of sample preparation and mass spectrometry measurement from impeding the measurement of protein concentration levels in the sample. A skilled artisan could readily recognize that any normalization methods and/or normalizers suitable for the present invention can be utilized.
The term “condition” as used herein refers generally to a disease, event, or change in health status.
The term “treatment protocol” as used herein includes further diagnostic testing typically performed to determine whether a pulmonary nodule is benign or malignant. Treatment protocols include diagnostic tests typically used to diagnose pulmonary nodules or masses such as for example, CT scan, positron emission tomography (PET) scan, bronchoscopy or tissue biopsy. Treatment protocol as used herein is also meant to include therapeutic treatments typically used to treat malignant pulmonary nodules and/or lung cancer such as for example, chemotherapy, radiation or surgery.
The terms “diagnosis” and “diagnostics” also encompass the terms “prognosis” and “prognostics”, respectively, as well as the applications of such procedures over two or more time points to monitor the diagnosis and/or prognosis over time, and statistical modeling based thereupon. Furthermore the term diagnosis includes: a. prediction (determining if a patient will likely develop a hyperproliferative disease); b. prognosis (predicting whether a patient will likely have a better or worse outcome at a pre-selected time in the future); c. therapy selection; d. therapeutic drug monitoring; and e. relapse monitoring.
In some embodiments, for example, classification of a biological sample as being derived from a subject with a lung condition may refer to the results and related reports generated by a laboratory, while diagnosis may refer to the act of a medical professional in using the classification to identify or verify the lung condition.
The term “providing” as used herein with regard to a biological sample refers to directly or indirectly obtaining the biological sample from a subject. For example, “providing” may refer to the act of directly obtaining the biological sample from a subject (e.g., by a blood draw, tissue biopsy, lavage and the like). Likewise, “providing” may refer to the act of indirectly obtaining the biological sample. For example, providing may refer to the act of a laboratory receiving the sample from the party that directly obtained the sample, or to the act of obtaining the sample from an archive.
As used herein, “lung cancer” preferably refers to cancers of the lung, but may include any disease or other disorder of the respiratory system of a human or other mammal. Respiratory neoplastic disorders include, for example small cell carcinoma or small cell lung cancer (SCLC), non-small cell carcinoma or non-small cell lung cancer (NSCLC), squamous cell carcinoma, adenocarcinoma, broncho-alveolar carcinoma, mixed pulmonary carcinoma, malignant pleural mesothelioma, undifferentiated large cell carcinoma, giant cell carcinoma, synchronous tumors, large cell neuroendocrine carcinoma, adenosquamous carcinoma, undifferentiated carcinoma; and small cell carcinoma, including oat cell cancer, mixed small cell/large cell carcinoma, and combined small cell carcinoma; as well as adenoid cystic carcinoma, hamartomas, mucoepidermoid tumors, typical carcinoid lung tumors, atypical carcinoid lung tumors, peripheral carcinoid lung tumors, central carcinoid lung tumors, pleural mesotheliomas, and undifferentiated pulmonary carcinoma and cancers that originate outside the lungs such as secondary cancers that have metastasized to the lungs from other parts of the body. Lung cancers may be of any stage or grade. Preferably the term may be used to refer collectively to any dysplasia, hyperplasia, neoplasia, or metastasis in which the protein biomarkers expressed above normal levels as may be determined, for example, by comparison to adjacent healthy tissue.
Examples of non-cancerous lung condition include chronic obstructive pulmonary disease (COPD), benign tumors or masses of cells (e.g., hamartoma, fibroma, neurofibroma), granuloma, sarcoidosis, and infections caused by bacterial (e.g., tuberculosis) or fungal (e.g., histoplasmosis) pathogens. In certain embodiments, a lung condition may be associated with the appearance of radiographic PNs.
As used herein, “lung tissue” and “lung cancer” refer to tissue or cancer, respectively, of the lungs themselves, as well as the tissue adjacent to and/or within the strata underlying the lungs and supporting structures such as the pleura, intercostal muscles, ribs, and other elements of the respiratory system. The respiratory system itself is taken in this context as representing nasal cavity, sinuses, pharynx, larynx, trachea, bronchi, lungs, lung lobes, aveoli, aveolar ducts, aveolar sacs, aveolar capillaries, bronchioles, respiratory bronchioles, visceral pleura, parietal pleura, pleural cavity, diaphragm, epiglottis, adenoids, tonsils, mouth and tongue, and the like. The tissue or cancer may be from a mammal and is preferably from a human, although monkeys, apes, cats, dogs, cows, horses and rabbits are within the scope of the present invention. The term “lung condition” as used herein refers to a disease, event, or change in health status relating to the lung, including for example lung cancer and various non-cancerous conditions.
“Accuracy” refers to the degree of conformity of a measured or calculated quantity (a test reported value) to its actual (or true) value. Clinical accuracy relates to the proportion of true outcomes (true positives (TP) or true negatives (TN)) versus misclassified outcomes (false positives (FP) or false negatives (FN)), and may be stated as a sensitivity, specificity, positive predictive values (PPV) or negative predictive values (NPV), or as a likelihood, odds ratio, among other measures. The term “biological sample” as used herein refers to any sample of biological origin potentially containing one or more biomarker proteins. Examples of biological samples include tissue, organs, or bodily fluids such as whole blood, plasma, serum, tissue, lavage or any other specimen used for detection of disease.
The term “subject” as used herein refers to a mammal, preferably a human.
The term “biomarker protein” as used herein refers to a polypeptide in a biological sample from a subject with a lung condition versus a biological sample from a control subject. A biomarker protein includes not only the polypeptide itself, but also minor variations thereof, including for example one or more amino acid substitutions or modifications such as glycosylation or phosphorylation.
The term “biomarker protein panel” as used herein refers to a plurality of biomarker proteins. In certain embodiments, the expression levels of the proteins in the panels can be correlated with the existence of a lung condition in a subject. In certain embodiments, biomarker protein panels comprise 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 60, 70, 80, 90 or 100 proteins. In certain embodiments, the biomarker proteins panels comprise 2-5 proteins, 5-10 proteins, 10-20 proteins or more.
“Treating” or “treatment” as used herein with regard to a condition may refer to preventing the condition, slowing the onset or rate of development of the condition, reducing the risk of developing the condition, preventing or delaying the development of symptoms associated with the condition, reducing or ending symptoms associated with the condition, generating a complete or partial regression of the condition, or some combination thereof.
Biomarker levels may change due to treatment of the disease. The changes in biomarker levels may be measured by the present invention. Changes in biomarker levels may be used to monitor the progression of disease or therapy.
“Altered”, “changed” or “significantly different” refer to a detectable change or difference from a reasonably comparable state, profile, measurement, or the like. One skilled in the art should be able to determine a reasonable measurable change. Such changes may be all or none. They may be incremental and need not be linear. They may be by orders of magnitude. A change may be an increase or decrease by 1%, 5%, 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 99%, 100%, or more, or any value in between 0% and 100%. Alternatively the change may be 1-fold, 1.5-fold 2-fold, 3-fold, 4-fold, 5-fold or more, or any values in between 1-fold and five-fold. The change may be statistically significant with a p value of 0.1, 0.05, 0.001, or 0.0001.
Using the methods of the current invention, a clinical assessment of a patient is first performed. If there exists is a higher likelihood for cancer, the clinician may rule in the disease which will require the pursuit of diagnostic testing options yielding data which increase and/or substantiate the likelihood of the diagnosis. “Rule in” of a disease requires a test with a high specificity.
“FN” is false negative, which for a disease state test means classifying a disease subject incorrectly as non-disease or normal.
“FP” is false positive, which for a disease state test means classifying a normal subject incorrectly as having disease.
The term “rule in” refers to a diagnostic test with high specificity that optionally coupled with a clinical assessment indicates a higher likelihood for cancer. If the clinical assessment is a lower likelihood for cancer, the clinician may adopt a stance to rule out the disease, which will require diagnostic tests which yield data that decrease the likelihood of the diagnosis. “Rule out” requires a test with a high sensitivity. Accordingly, the term “ruling in” as used herein is meant that the subject is selected to receive a treatment protocol.
The term “rule out” refers to a diagnostic test with high sensitivity that optionally coupled with a clinical assessment indicates a lower likelihood for cancer. Accordingly, the term “ruling out” as used herein is meant that the subject is selected not to receive a treatment protocol.
The term “sensitivity of a test” refers to the probability that a patient with the disease will have a positive test result. This is derived from the number of patients with the disease who have a positive test result (true positive) divided by the total number of patients with the disease, including those with true positive results and those patients with the disease who have a negative result, i.e., false negative.
The term “specificity of a test” refers to the probability that a patient without the disease will have a negative test result. This is derived from the number of patients without the disease who have a negative test result (true negative) divided by all patients without the disease, including those with a true negative result and those patients without the disease who have a positive test result, e.g., false positive. While the sensitivity, specificity, true or false positive rate, and true or false negative rate of a test provide an indication of a test's performance, e.g., relative to other tests, to make a clinical decision for an individual patient based on the test's result, the clinician requires performance parameters of the test with respect to a given population.
The term “positive predictive value” (PPV) refers to the probability that a positive result correctly identifies a patient who has the disease, which is the number of true positives divided by the sum of true positives and false positives.
The term “negative predictive value” or “NPV” is calculated by TN/(TN+FN) or the true negative fraction of all negative test results. It also is inherently impacted by the prevalence of the disease and pre-test probability of the population intended to be tested. The term NPV refers to the probability that a negative test correctly identifies a patient without the disease, which is the number of true negatives divided by the sum of true negatives and false negatives. A positive result from a test with a sufficient PPV can be used to rule in the disease for a patient, while a negative result from a test with a sufficient NPV can be used to rule out the disease, if the disease prevalence for the given population, of which the patient can be considered a part, is known.
The term “disease prevalence” refers to the number of all new and old cases of a disease or occurrences of an event during a particular period. Prevalence is expressed as a ratio in which the number of events is the numerator and the population at risk is the denominator.
The term disease incidence refers to a measure of the risk of developing some new condition within a specified period of time; the number of new cases during some time period, it is better expressed as a proportion or a rate with a denominator.
Lung cancer risk according to the “National Lung Screening Trial” is classified by age and smoking history. High risk—age ≧55 and ≧30 pack-years smoking history; Moderate risk—age ≧50 and ≧20 pack-years smoking history; Low risk—<age 50 or <20 pack-years smoking history.
The clinician must decide on using a diagnostic test based on its intrinsic performance parameters, including sensitivity and specificity, and on its extrinsic performance parameters, such as positive predictive value and negative predictive value, which depend upon the disease's prevalence in a given population.
Additional parameters which may influence clinical assessment of disease likelihood include the prior frequency and closeness of a patient to a known agent, e.g., exposure risk, that directly or indirectly is associated with disease causation, e.g., second hand smoke, radiation, etc., and also the radiographic appearance or characterization of the pulmonary nodule exclusive of size. A nodule's description may include solid, semi-solid or ground glass which characterizes it based on the spectrum of relative gray scale density employed by the CT scan technology.
“Mass spectrometry” refers to a method comprising employing an ionization source to generate gas phase ions from an analyte presented on a sample presenting surface of a probe and detecting the gas phase ions with a mass spectrometer.
B. Development of InteQuan
Hundreds of human plasma samples were analyzed in a prior label-free discovery study of lung cancer biomarkers, using immunoaffinity-based protein depletion coupled with MRM-MS (depletion-MRM-MS) (Li, X J et al, Sci. Tranl. Med. 2013, 5:207ra142). Six endogenous normalizing proteins were identified in this study out of 371 protein candidates. Since the normalizing proteins were processed and analyzed together with target proteins of interest, it was expected that they would serve as monitors for systematic variation in both pre-analytical and analytical procedures ( FIG. 1B ). It has been verified that experimental variability was reduced after normalization by a panel of the six proteins.
Although this method of endogenous protein normalization (EPN) and similar approaches have been reported previously in quantitative proteomics, the approach taken by the present invention is unique. The six normalizing proteins were selected by their ability to compensate both the drift of depletion columns and the technical variation of other proteins, rather than their “housekeeping” properties as utilized in other approaches. However, the EPN quantification method used in our previous study is a label-free approach and cannot control analytical variability as narrowly as the three SIL approaches.
The present invention comprises more accurate MRM assays for target proteins of interest, using the SIL approach of SIS peptides. As described herein, neither SIL nor EPN is ideal for protein quantification in clinical proteomics. To deal with this challenge, the methods of the present invention describe a new quantification method, named integrated quantification (InteQuan), which uniquely combines the advantageous features of two methods: EPN and SISQuan. As further described in the Examples below, six endogenous normalizing proteins were used to compensate systematic variation in pre-analytical procedures and in depletion and digestion, while SIS peptides were used to control variation in desalting and MS analysis ( FIG. 1B ).
To mimic an actual clinical testing scenario, the details of the methods of the present invention were defined based on data from a different study (Vachani et al., Validation of a Multi-Protein Plasma Classifier to Identify Benign Lung Nodules, J. Thoracic Onc . (doi: 10.1097/JTO.0000000000000447)) before assessing its performance in three independent studies ( FIG. 1C ), Studies I-III, which are further described in the Examples. No one has previously described such a method to quantify individual proteins before. A recent work used “sparse reference labeling” to anchor protein abundance that can be valuable for biomarker discovery (Chang, C T et al, Targeted protein quantification using spare reference labeling, Nat. Methods 2014, 11:301-304). Nevertheless, individual proteins were essentially quantified in a label-free approach in the study, leaving peptide-specific variation in MS analysis uncontrolled and thus reducing its validity for routine clinical testing. The methods of the present invention, utilizing aspects of the newly developed InteQuan, demonstrate increased precision, robustness, and resolution of protein measurement in the three independent assessment studies, Studies I-III ( FIG. 1 C).
C. Advantages of InteQuan
At least three aspects of InteQuan and the methods of the present invention enhance their relevance towards the development of blood-based laboratory-developed tests using MRM-MS platforms. First, the target proteins used in the development of InteQuan were all potential lung cancer biomarkers. Second, endogenous proteins in low ng/ml to low μg/ml plasma concentrations were quantified in both clinical plasma samples and in the standard HPS samples. Third, the longitudinal assessment on the robustness of the depletion-MRM-MS platform was performed in settings similar to actual laboratory operations for clinical testing. In contrast, medium- to high-abundant endogenous proteins, spike-in peptides or spike-in proteins were quantified on single plasma samples, in settings of academic research rather than clinical testing, and using different MRM-MS platforms in other studies. Nevertheless, the precision obtained from the development of InteQuan was comparable to the precision reported in those studies. Furthermore, the precision of the whole depletion-MRM-MS platform was assessed in the development of InteQuan, not just the precision of MRM-MS platforms as in some studies.
In one aspect, the selection of proper endogenous normalizing proteins provides the advantageous performance of InteQuan relative to other quantification methods. According to error propagation theory, normalization by endogenous proteins has two opposite effects. On one hand, it reduces the overall variability in protein measurement by cancelling out systematic variation that similarly affects target and normalizing proteins. On the other hand, it increases the overall variability by transferring protein-specific and random variation of normalizing proteins to target proteins.
Three strategies were applied to ensure a favorable outcome from the normalization process utilized in various aspects of the methods of the present invention. First, a large dataset was generated to capture both technical variability on the platform and biological variability of the intended patient population. Owing to considerations of cost, a label-free approach was used to quantify proteins in the study. Second, six normalizing proteins were selected for their performance in reducing column drift and technical CV of other proteins. In other words, the proteins were specifically selected to fulfill the role of normalizers. Third, the six normalizing proteins were used as a panel that was more stable compared to individual proteins. In addition, although plasma concentration was not used as a selection criterion, the wide concentration range (three orders of magnitude) of the six normalizing proteins was likely beneficial. Similar strategies should be used for selecting endogenous normalizing proteins on other MS platforms and/or for other proteomics projects. In our case, the six normalizing proteins were selected from a pool of 371 protein candidates based on a set of label-free depletion-MRM-MS data.
It is possible that other proteins outside the pool may serve as normalizers for on other MS platforms. As shown later in FIG. 8 , different combinations of normalizers may be used. Accordingly, the methods of endogenous normalizing protein selection described herein may be used to select normalizing proteins for different MS platform techniques and/or focused on quantifying biomarkers for diseases other than lung disease.
InteQuan measures the abundance of the target proteins relative to the abundance of the endogenous normalizing proteins, which explains its high tolerance against variation in the total protein concentration. When testing actual clinical samples, pre-analytical variability (due to differences in patient posture, diurnal cycle, sample collection, and/or sample handling, etc.) and analytical variability (due to differences in sample loading volume, instrument performance, and/or operator, etc.) are hard to avoid and all contribute to the overall variability of the assay. Thus, a high tolerance against such variation is a desirable feature that will increase the reproducibility of clinical tests.
In summary, InteQuan has been developed as a quantification method for, e.g, MS-based quantitative proteomics, and has demonstrated its superiority to SISQuan in three independent studies and on a combined HPS dataset described in the Examples. The InteQuan method is robust, simple to implement, capable of reducing pre-analytical and analytical variability, and able to improve the measurement of biological differences. All these features make the method an ideal technique for MS-based quantitative proteomics, e.g., for applications in biomarker research and in routine clinical testing.
D. Applications of InteQuan
According to one aspect of the methods of the present invention, the abundance or expression level of protein biomarkers for lung disease may be measured by MS. MS analyzes the mass spectrum produced by an ion after its production by the vaporization of its parent protein and its separation from other ions based on its mass-to-charge ratio. The most common modes of acquiring MS data are 1) full scan acquisition resulting in the typical total ion current plot (TIC), 2) selected ion monitoring (SIM), and 3) selected reaction monitoring (SRM).
In certain embodiments of the methods provided herein, biomarker protein expression levels are measured by LC-SRM-MS. LC-SRM-MS is a highly selective method of tandem mass spectrometry which has the potential to effectively filter out all molecules and contaminants except the desired analyte(s). This is particularly beneficial if the analysis sample is a complex mixture which may comprise several isobaric species within a defined analytical window. LC-SRM-MS methods may utilize a triple quadrupole mass spectrometer which, as is known in the art, includes three quadrupole rod sets. A first stage of mass selection is performed in the first quadrupole rod set, and the selectively transmitted ions are fragmented in the second quadrupole rod set. The resultant transition (product) ions are conveyed to the third quadrupole rod set, which performs a second stage of mass selection. The product ions transmitted through the third quadrupole rod set are measured by a detector, which generates a signal representative of the numbers of selectively transmitted product ions. The RF and DC potentials applied to the first and third quadrupoles are tuned to select (respectively) precursor and product ions that have m/z values lying within narrow specified ranges. By specifying the appropriate transitions (m/z values of precursor and product ions), a peptide corresponding to a targeted protein may be measured with high degrees of sensitivity and selectivity. Signal-to-noise ratio is superior to conventional tandem mass spectrometry (MS/MS) experiments, which select one mass window in the first quadrupole and then measure all generated transitions in the ion detector. LC-SRM-MS.
The expression level of a biomarker protein can be measured using any suitable method known in the art, including but not limited to mass spectrometry (MS), reverse transcriptase-polymerase chain reaction (RT-PCR), microarray, serial analysis of gene expression (SAGE), gene expression analysis by massively parallel signature sequencing (MPSS), immunoassays (e.g., ELISA), immunohistochemistry (IHC), transcriptomics, and proteomics.
According to one aspect of the methods, the measuring step is performed by selected reaction monitoring mass spectrometry, using a compound that specifically binds the protein being detected or a peptide transition. In one embodiment, the compound that specifically binds to the protein being measured is an antibody or an aptamer.
According to one aspect of the methods, the biological sample includes such as for example tissue, blood, plasma, serum, whole blood, urine, saliva, genital secretion, cerebrospinal fluid, sweat and excreta.
In some embodiments of the methods, the determination of the abundance or expression level of the target proteins provides a diagnosis of lung disease. Accordingly, the methods of the present invention may serve as diagnostic methods.
In certain embodiments, the diagnostic methods disclosed herein can be used in combination with other clinical assessment methods, including for example various radiographic and/or invasive methods. Similarly, in certain embodiments, the diagnostic methods disclosed herein can be used to identify candidates for other clinical assessment methods, or to assess the likelihood that a subject will benefit from other clinical assessment methods.
To evaluate the diagnostic performance of a particular set of peptide transitions, a ROC curve is generated for each significant transition.
An “ROC curve” as used herein refers to a plot of the true positive rate (sensitivity) against the false positive rate (specificity) for a binary classifier system as its discrimination threshold is varied. A ROC curve can be represented equivalently by plotting the fraction of true positives out of the positives (TPR=true positive rate) versus the fraction of false positives out of the negatives (FPR=false positive rate). Each point on the ROC curve represents a sensitivity/specificity pair corresponding to a particular decision threshold.
AUC represents the area under the ROC curve. The AUC is an overall indication of the diagnostic accuracy of 1) a biomarker or a panel of biomarkers and 2) a ROC curve. AUC is determined by the “trapezoidal rule.” For a given curve, the data points are connected by straight line segments, perpendiculars are erected from the abscissa to each data point, and the sum of the areas of the triangles and trapezoids so constructed is computed. In certain embodiments of the methods provided herein, a biomarker protein has an AUC in the range of about 0.75 to 1.0. In certain of these embodiments, the AUC is in the range of about 0.8 to 0.85, 0.85 to 0.9, 0.9 to 0.95, or 0.95 to 1.0.
According to one aspect of the methods of the present invention, the subject has or is suspected of having a pulmonary nodule or a pulmonary mass. The pulmonary nodule has a diameter of less than or equal to 3.0 cm. The pulmonary mass has a diameter of greater than 3.0 cm. In some embodiments, the pulmonary nodule has a diameter of about 0.8 cm to 3.0 cm. The subject may have stage IA lung cancer (i.e., the tumor is smaller than 3 cm).
The high abundance of certain proteins in a biological sample such as plasma or serum can hinder the ability to assay a protein of interest, particularly where the protein of interest is expressed at relatively low concentrations. Several methods are available to circumvent this issue, including enrichment, separation, and depletion. Enrichment uses an affinity agent to extract proteins from the sample by class, e.g., removal of glycosylated proteins by glycocapture. Separation uses methods such as gel electrophoresis or isoelectric focusing to divide the sample into multiple fractions that largely do not overlap in protein content. Depletion typically uses affinity columns to remove the most abundant proteins in blood, such as albumin, by utilizing advanced technologies such as IgY14/Supermix (Sigma St. Louis, Mo.) that enable the removal of the majority of the most abundant proteins.
In certain embodiments of the methods provided herein, a biological sample may be subjected to enrichment, separation, and/or depletion prior to assaying biomarker or putative biomarker protein expression levels. In certain of these embodiments, blood proteins may be initially processed by a glycocapture method, which enriches for glycosylated proteins, allowing quantification assays to detect proteins in the high pg/ml to low ng/ml concentration range. Exemplary methods of glycocapture are well known in the art (see, e.g., U.S. Pat. No. 7,183,188; U.S. Patent Appl. Publ. No. 2007/0099251; U.S. Patent Appl. Publ. No. 2007/0202539; U.S. Patent Appl. Publ. No. 2007/0269895; and U.S. Patent Appl. Publ. No. 2010/0279382). In other embodiments, blood proteins may be initially processed by a protein depletion method, which allows for detection of commonly obscured biomarkers in samples by removing abundant proteins. In one such embodiment, the protein depletion method is a Supermix (Sigma) depletion method.
The following examples are provided to better illustrate the claimed invention and are not to be interpreted as limiting the scope of the invention. To the extent that specific materials are mentioned, it is merely for purposes of illustration and is not intended to limit the invention. One skilled in the art may develop equivalent means or reactants without the exercise of inventive capacity and without departing from the scope of the invention
EXAMPLES
Example 1
Protein Quantification in Human Plasma Samples
A total of 21 lung cancer biomarker candidates were identified in a prior discovery study (Li, X J et al, Sci. Tranl. Med. 2013, 5:207ra142). Two of the 21 candidates (GSLG1 and EF1A1) were eliminated from this study due to weak signals on a new MRM-MS platform. Another candidate (FIBA) was eliminated due to its affinity to the depletion column. With reference to Table 1, the remaining 18 candidates were targeted for quantification in human plasma samples.





TABLE 1



List of six normalizing proteins and eighteen target proteins of interest









Con-






centration a
Transition b




Protein name
(ng/ml)
(peptide (SEQ ID NO)_Q1_Q3)
FDR c
F n d
{hacek over (R)} n e










Pigment epithelium-derived
7200
LQSLFDSPDFSK (SEQ ID NO: 25)
1.40E−04
0.971
1.756

factor

_692.34_593.30




Mannan-binding lectin serine
240
TGVITSPDFPNPYPK (SEQ ID NO: 26)
5.75E−04
0.957
0.360

protease 1

_816.92_258.10




Gelsolin
16000
TASDFITK(SEQ ID NO: 27) _441.73_710.40
3.18E−04
0.852
0.502

Lumican
4000
SLEDLQLTHNK (SEQ ID NO: 28)_433.23_499.30
3.82E−04
0.838
10.846

Scavenger receptor cysteine-
94
INPASLDK (SEQ ID NO: 29)_429.24_630.30
1.19E−03
0.823
0.392

rich type 1 protein M130






Receptor-type tyrosine-protein
9.9
VITEPIPVSDLR (SEQ ID NO: 30)_669.89_896.50
1.44E−03
0.926
0.275

phosphatase eta












Apoptosis-inducing factor 1,
1.4
ELWFSDDPNVTK (SEQ ID NO: 31)_725.85_558.30
3.70E−02
Assay specificity not verified

mitochondrial












Mast/stem cell growth factor
8.2
YVSELHLTR (SEQ ID NO: 32) _373.21_428.30
2.40E−03
0.730


receptor






Ferritin light chain
12
LGGPEAGLGEYLFER(SEQ ID NO: 33)
4.30E−05
0.844




_804.40_1083.60










Prolow-density lipoprotein
20
TVLWPNGLSLDIPAGR(SEQ ID NO: 34)
1.40E−04
Assay specificity not verified

receptor-related protein 1

_855.00_1209.70










Collagen alpha-1(XVIII) chain
35
AVGLAGTFR(SEQ ID NO: 35) _446.26_721.40
6.70E−04
0.732


Peroxiredoxin-1
60
QITVNDLPVGR (SEQ ID NO: 36)_606.30_970.50
1.90E−05
1.714


Tenascin-X
70
YEVTVVSVR(SEQ ID NO: 37) _526.29_293.10
1.10E−03
0.699


Endoplasmin
88
SGYLLPDTK(SEQ ID NO: 38) _497.27_308.10
1.10E−03
0.649


78 kDa glucose-regulated
100
TWNDPSVQQDIK (SEQ ID NO: 39)
1.80E−03
1.140


protein

_715.85_288.10




Transforming growth factor-beta-
140
LTLLAPLNSVFK(SEQ ID NO: 40)
1.40E−04
0.779


induced protein ig-h3

_658.40_804.50




Fructose-bisphosphate aldolase
250
ALQASALK(SEQ ID NO: 41) _401.25_617.40
3.70E−05
0.777


A






Gamma-glutamyl hydrolase
250
YYIAASYVK(SEQ ID NO: 42) _539.28_638.40
1.70E−03
0.834


Monocyte differentiation
420
ATVNPSAPR(SEQ ID NO: 43) _456.80_527.30
4.30E−04
0.789


antigen CD14






Galectin-3-binding protein
440
VEIFYR (SEQ ID NO: 44)_413.73_598.30
2.80E−05
0.842


Thrombospondin-1
510
GFLLLASLR(SEQ ID NO: 45) _495.31_559.40
1.90E−05
0.625


Insulin-like growth factor-
5700
FLNVLSPR (SEQ ID NO: 46)_473.28_685.40
2.80E−05
0.790


binding protein 3






Tetranectin
58000
LDTLAQEVALLK(SEQ ID NO: 47) _657.39_871.50
3.70E−05
0.760


Immunoglobulin superfamily

ALPGTPVASSQPR(SEQ ID NO: 48)
4.40E−03
0.850


containing leucine-rich repeat

_640.85_841.50




protein



a Predicted plasma concentration [26].

b Thc transition that was used for quantification.

c False discovery rate for peptide MRM assay (peptide Q value) [6],

d Correction factor {F n } in Study II in which a new lot of SIS peptides were used.

e Scaling constant {{hacek over (R)} n } for InteQuan.

f Scaling constant {{hacek over (A)} n } for EPN.
















































































With further reference to Table 1, six endogenous normalizing proteins were selected from a pool of 371 protein candidates in a previous label-free discovery study (Li, X J et al, Sci. Tranl. Med. 2013, 5:207ra142). The predicted plasma concentration of the six proteins, estimated from the occurrence of protein detection in human plasma or serum samples by the proteomics community, ranged from 9.9 ng/ml (PTPRJ) to 16 μg/ml (GELS). All six proteins were used as normalizing proteins for InteQuan and for EPN in this study.
Human plasma samples were analyzed on a depletion-MRM-MS platform. SIS peptides of the target and the normalizing proteins were synthesized and spiked into peptide samples after digestion. The specificity of MRM assays to the corresponding proteins was verified for all proteins except for LRP1 and AIFM1. As shown in FIG. 2 , MRM signals of verified assays were well above the corresponding noise level; endogenous and SIS peptides coeluted and had comparable intensity ratios between different transitions. The highest false discovery rate (FDR) of the original assays was 0.44% (ISLR, see Table 1).
As shown later in FIG. 7B , proteins were also measured within the respective linear dynamic range of the assays. Two blank samples were processed and analyzed at the end of each experimental batch in Study II and III to monitor possible carryover from previous samples (see Table 2). MRM signals in those blank samples were just above noise level (data not shown), indicating that carryover was not a problem for the depletion-MRM-MS platform. After validating the MRM assays, LRP1 and AIFM1 were both eliminated from further analysis. The predicted plasma concentration of the 16 remaining target proteins spanned four orders of magnitude from 8.2 ng/ml (KIT) to 58 μg/ml (TETN).





TABLE 2



Sample layout in the three assessment studies a,b













a HPS, human plasma standard. NC, negative control (0.15M (NH 4 )HCO 3 ).

b Samples with no data were highlighted in box and were lost due to 1) unexpected

pump stoppage and restarting during depletion (the first HPS in Study I), 2) sample

exhaustion (the sample in Study III) or 3) sample contamination during solid-phase

extraction (all other samples).

c The loading volume for sample labeled as L0, L1, and L2 was 30, 20 or 15 μl,

respectively.





















The target proteins were quantified based on MRM-MS data using four different methods (raw MS data, EPN, SISQuan and InteQuan). In this study, the abundance of a protein was evaluated based on the MRM signal intensity of the strongest transition from the protein. Thus, no distinction was made between protein abundance, peptide abundance and transition abundance. For raw MS data, protein abundance was measured by the peak area of the strongest transition of the protein. For EPN, protein abundance using the raw MS data was divided by a sample-dependent normalization factor that was calculated from the peak areas of the six normalizing proteins. Six scaling constants, one for each of the six normalizing proteins, were used in the calculation of the normalization factor. For SISQuan, protein abundance was measured by the response ratio between the peak area of the strongest transition of the target protein and the peak area of the matching transition of the corresponding SIS peptide. For InteQuan, protein abundance using SISQuan was divided by a sample-dependent normalization factor that was calculated from the response ratios of the six normalizing proteins. As with EPN, six scaling constants were used in the calculation of the normalization factor. In the study, we mainly focused on comparing the new InteQuan method with the widely used SISQuan method.
With reference to Table 1, all scaling constants for InteQuan and for EPN were determined from a different study of 100 clinical samples and 20 aliquots of a human plasma standard (HPS) sample. See Vachani et al., Validation of a Multi-Protein Plasma Classifier to Identify Benign Lung Nodules, J. Thoracic Oncology (doi: 10.1097/JTO.0000000000000447). None of the scaling constants were modified in this study. Therefore, the assessment of the four quantification methods was based on independent datasets.
A. Clinical Samples
Archival K2-EDTA plasma samples were obtained from subjects that provided informed consent and with approval by either the Ethics Review Board at Institut Universitaire de Cardiologie et de Pneumologie de Quebec or the Institutional Review Boards at New York University and University of Pennsylvania. All samples were collected prior to surgery or from patients without surgery. Disease status of patients was histopathologically confirmed. All cancer patients were at Stage I or II. Clinical data associated with subjects were handled in accordance with the guidance established by the Health Insurance Portability and Accountability Act of 1996 to ensure subject privacy.
B. Selection of Endogenous Normalizing Proteins
In a prior discovery study, 72 cancer and 71 benign samples were analyzed in five experimental batches along with 15 aliquots of a pooled HPS sample that was purchased from Bioreclamation (Hicksville, N.Y.) (Li, X J et al, Sci. Tranl. Med. 2013, 5:207ra142). The HPS samples were embedded among clinical samples and analyzed repeatedly to monitor analytical variability in the experiment. The clinical samples were used to represent biological variability and possible pre-analytical variability.
Endogenous normalizing proteins were selected from proteins whose strongest transitions were detected in all samples. Each protein candidate was used to normalize the abundance of other proteins and evaluated based on the following criteria: (A) Its rank, as a normalizer, in reducing median technical CV of other proteins; (B) its rank in compensating median column drift, that is a technical variation associated with depletion; (C) its own median technical CV on HPS samples; and (D) its own median biological CV on clinical samples.
In the end, six endogenous normalizing proteins were selected. Owing to considerations of cost, the selection of endogenous normalizing proteins was performed in a label-free approach.
C. Immunoaffinity Chromatography
Experimental protocols for sample preparation were adapted and modified from a recent study (Li, X J. Sci. Transl. Med. 2013, 5:207ra142). Immunoaffinity columns containing a 2:1 ratio of IgY14 and SuperMix resins were purchased from Sigma Aldrich (St. Louis). Each column was conditioned with 0.15 M (NH 4 )HCO 3 at 0.5 ml/min for 45 min. Prior to immunoaffinity separation of each sample batch, column performance was assessed with replicate injections of aliquots of the HPS sample.
To isolate low abundance proteins, 45, 50, or 60 μl of plasma were diluted in 0.15M (NH 4 )HCO 3 to a final volume of 135, 150, or 180 μl, respectively, and filtered using a 0.45 μm AcroPrep 96-well filter plate (Pall Life Sciences). Immunoaffinity separation was conducted on a IgY14-SuperMix column connected to an high performance liquid chromatography (HPLC) system (Agilent 1260 Infinity Bioinert Quaternary liquid chromatography (LC)) using 3 buffers (loading/washing: 0.15 M (NH 4 )HCO 3 ; stripping/elution: 0.1 M glycine, pH 2.5; and neutralization: 0.01 M Tris-HCl and 0.15 M NaCl, pH 7.4) with a cycle comprised of load, wash, elute, neutralization, and re-equilibration lasting 36 min. The total plasma volume loaded onto the depletion column was 15, 20, or 30 μl, respectively. The unbound and bound fractions were monitored at 280 nm and were baseline resolved after separation. Unbound fractions (containing the low abundance proteins) were collected for downstream processing and analysis and lyophilized prior to enzymatic digestion. Every 24 samples were grouped as an experimental batch and were processed sequentially in a throughput of one batch per day.
D. Enzymatic Digestion
Lyophilized fractions containing low abundance proteins were digested with trypsin after being reconstituted under mild denaturing conditions in 200 μl of 1:1 0.1 M (NH 4 )HCO 3 trifluoroethanol (TFE) (v/v) and then allowed to incubate on an orbital shaker for 30 min at room temperature (RT). Samples were diluted in 800 μl of 0.1 M (NH 4 )HCO 3 and digested with 0.4 μg trypsin (Princeton Separations) per sample for 16+/−2 hours at 37° C. Following digestion samples were stored at −70° C. for 2 hours and then lyophilized. Samples within each study were digested in parallel.
E. Stable Isotope-labeled Standard Peptides
A total of 26 SIS peptides were purchased from New England Peptide (Gardner, Mass.), including one SIS peptide for each of the six normalizing proteins and the 18 target proteins in Table 1. SIS peptides of two additional proteins (S10A6 and PROF1) were included as potential biomarkers earlier on but were later eliminated. Each SIS peptide was purified to 95% or greater as determined by reversed phase HPLC; mass determination for each peptide was confirmed to be within 0.1% of the calculated mass by matrix-assisted laser desorption/ionization (MALDI)-time of flight (TOF) MS. The concentration of the stock solution for each peptide was determined by amino acid analysis. The SIS peptide mixture was produced per specified formulation in 10% acetonitrile, 0.1% formic acid final concentration with 100 fmol/μL BSA digest added for stability. Concentrations of individual SIS peptides were tailored so that their MRM-MS signal intensities were comparable to those of the corresponding endogenous peptides. The mixture was aliquoted into individual 300 μL single use microfuge tubes and stored at −80° C. Aliquots of the SIS peptide mixture were thawed on wet ice, mixed briefly and spiked into peptide samples after enzymatic digestion and lyophilization and during solubilization just prior to solid-phase extraction. Two different preparations (lots) of the SIS peptide mixture were prepared and used in this study. The stability of SIS peptides was monitored based on their MRM signal intensities. No evidence for the instability of SIS peptides was observed over a period of 20 months (data not shown).
F. Solid-Phase Extraction
Aliquots of the SIS peptide mixture were spiked into the lyophilized peptide samples, followed by reconstitution in 350 μl of 0.01 M (NH 4 )HCO 3 , incubation on an orbital shaker for 15 min at RT, reduction using 30 μl of 0.05 M TCEP, incubation for 1 hour at RT, and dilution in 375 μl of 90% water/10% acetonitrile/0.2% trifluoroacetic acid. The solid phase extraction plate (Empore C18, 3M Bioanalytical Technologies) was conditioned according to the manufacturer's protocol, and after sample loading were washed in 500 μl of 95% water/5% acetonitrile/0.1% trifluroacetic acid and eluted by 200 μl of 52% water/48% acetonitrile/0.1% trifluoroacetic acid into a collection plate. The eluate was split into 2 equal aliquots and was taken to dryness in a vacuum concentrator. One aliquot was used immediately for mass spectrometry, while the other was stored at −80° C. Samples were reconstituted in 12 μl of 90% water/10% acetonitrile/0.2% formic acid just prior to LC-MRM-MS analysis. Samples within each study were processed in parallel in this step.
G. Optimization of MRM Assays
MRM assays of endogenous peptides of the target and normalizing proteins were developed previously on a 5500 QTrap® reversed-phase LC-MRM-MS platform (AB Sciex). The specificity of the assays was verified with a FDR of 3.70% or lower. These assays, along with MRM assays of the corresponding SIS peptides, were transferred to and optimized on a 6490 Triple Quadrupole LC-MRM-MS platform (Agilent) based on the highly purified synthetic SIS peptides. The optimal assays were further tested on processed HPS samples to check for signal intensity and possible interference. Unless specified, the signal of the assays was well above noise and within the respective linear dynamic range. In addition to the low FDRs of the original assays, the specificity of the transitions to the corresponding proteins was further verified from the co-elution of endogenous and SIS peptides and from the consistency between the peptides on intensities of different transitions. Seventeen additional proteins were analyzed for exploratory purposes without optimizing their transitions or spiking in the corresponding SIS peptides. The 17 extra proteins were not analyzed in this study. A total of 302 transitions from 38 proteins were measured in this study.
H. MRM-MS Analysis
Peptide samples were separated using a capillary reversed-phase LC column (Agilent Poroshell 120 EC-C18; 2.1 mm×100 mm, particle size 2.7 μm) and an Agilent 1290 Infinity HPLC system. The mobile phases were (A) 0.1% formic acid in water and (B) 0.1% formic acid in acetonitrile. The samples were injected (8 μl) and separated using a linear gradient (98% A to 70% A) at 0.4 mL/minute for 21.7 min. Peptides were eluted directly into the electrospray source of the mass spectrometer (6490 Triple Quadrupole, Agilent) operating in scheduled MRM positive-ion mode (Q1 resolution: wide; Q3 resolution: unit; detection window variable: 124 to 240 seconds; cycle time: 1.0 seconds). Peak areas of transitions were integrated by MassHunter (Agilent) and manually curated to ensure quality. Samples within each experiment were analyzed sequentially.
I. Four Quantification Methods
In this study the abundance of a protein was evaluated based on the MRM signal intensity of the strongest transition of the protein and no distinctions between protein abundance, peptide abundance and transition abundance were made. Without losing generality, the four quantification methods were described in terms of peptide quantification as follows.
Raw MS data. In this label-free quantification approach, the abundance of peptide p in sample s was measured by its raw peak area (A p,s ) without normalization.
Endogenous protein normalization (EPN). In this label-free quantification approach, the abundance of peptide p in sample s was measured by its normalized peak area à p,s =A p,s /S s E , where S s E was a sample-dependent normalization factor and was calculated from the peak areas of a predetermined set of N=6 endogenous, normalizing peptides in the sample. More specifically,
S s E = median ⁡ ( A 1 , s A ⋓ 1 , A 2 , s A ⋓ 2 , … ⁢ , A N , s A ⋓ N ) ( EQN . ⁢ 1 )
where A p,s was the peak area of peptide normalizer n (with n=1, . . . , N) in the sample and {hacek over (A)} n was a scaling constant for the normalizer that ensured values of {A n,s /{hacek over (A)} n } among all normalizers to be the same on average. The scaling constants {{hacek over (A)} n } were determined as the median values (over all clinical samples) of {A n,s } in an independent study of 120 samples.

Quantification using SIS peptides (SISQuan). In this labeled quantification approach, the abundance of peptide p in sample s was measured by the response ratio between the endogenous peptide to the corresponding SIS peptide, that is R p,s =A p,s / p,s where  p,s was the peak area of the SIS peptide.
Integrated quantification (InteQuan). In this labeled quantification approach, the abundance of peptide p in sample s was measured by its normalized response ratio {tilde over (R)} p,s =R p,s /S s I , where S s I was a sample-dependent normalization factor and was calculated from the response ratios of the N peptide normalizers in the sample. More specifically,
S s I = median ⁡ ( R 1 , s R ⋓ 1 , R 2 , s R ⋓ 2 , … ⁢ , R N , s R ⋓ N ) ( EQN . ⁢ 2 )
where R n,s was the response ratio of peptide normalizer n in the sample and {hacek over (R)} n was a scaling constant for the normalizer that ensured values of {R n,s /{hacek over (R)} n } among all normalizers to be same on average. Similar to {{hacek over (A)} n }, the scaling constants {{hacek over (R)} n } were determined as the median values (over all clinical samples) of {R n,s } in the same study of 120 samples.

J. Migration to new lot of SIS peptides
Six aliquots of the HPS sample (30 μl per aliquot) were processed and pooled together after digestion. The pooled sample was split into two identical aliquots. Two lots of SIS peptide mixtures (old and new) were each spiked into one of the two aliquots of HPS. The two aliquots of SIS peptide/HPS mixture were then each further split into three equal aliquots and lyophilized. The SIS peptide/HPS mixtures were reconstituted, desalted, lyophilized, and stored. The SIS peptide/HPS samples were then solubilized and analyzed by MRM-MS. A correction factor was calculated for each peptide as F p ={circumflex over (R)} p,old /{circumflex over (R)} p,new , where {circumflex over (R)} p,old ({circumflex over (R)} p,new ) was the median response ratio of peptide p as evaluated using the old (new) lot of SIS peptides. In Study II, the abundance ratio R′ p,s of peptide p in sample s as measured using the new lot was multiplied by the correction factor F p , that is R p,s =R′ p,s *F p . This correction was applied to both the target and the normalizing peptides. Afterwards, the evaluation of protein abundance using InteQuan and using SISQuan were both based on the corrected abundance ratios {R p,s }.
N. Differences Between Study I and a Previous Study
All clinical samples in Study I have been previously processed and analyzed by a contract research organization (CRO; Caprion, Montreal). Similar protocols were used in immunoaffinity depletion, protein digestion and desalting (Li, X J et al, Sci. Tranl. Med. 2013, 5:207ra142). Major differences between the two studies included: 1) Laboratory: Study I was carried out in-house but the discovery study was done by the CRO. 2) Depletion: The depletion column was ordered directly from vendor in Study I but packed by the CRO with a different lot of IgY14-Supermix resin beads in the discovery study. 3) Quantification: SIS peptides were used for quantification in Study I but not in the discovery study. 4) MS platform: Peptides were analyzed by an Agilent 6490 Triple Quadrupole LC/MS System in Study I but by an AB SCIEX QTrap® 5500 LC/MS system in the discovery study. 5) Monitored transitions: 302 transitions of 38 proteins were monitored in Study I. In comparison, 1550 transitions of 344 proteins were monitored in the discovery study.
O. Statistical Data Analysis
Data analysis was performed using the R statistical environment. Code for PVCA was adapted from: [http://www.niehs.nih.gov/research/resources/software/biostatist-ics/pvca/], setting the threshold to capture at least 90% of variance and a minimum of two principal components. The p value for comparing different quantification methods was based on the most-applicable, non-parametric paired sign test, assuming that measurements were independent and from a continuous population. The p value was evaluated using the function “SIGN.test” in the “BSDA” library. Functions “glm” and “predict” were used to train and test logistic regression models. Function “lm.fit” in the “stats” library was used to fit the linear relationship between the SISQuan abundances and the sample loading volumes.
Example 2
Demonstration of Complementary Control of Variation by SISQuan and EPN
In Study I, 60 clinical samples and 12 aliquots of the HPS sample were analyzed in three experimental batches using one depletion column ( FIG. 1C and Table 2). Clinical information of the patients is listed in Table 3. MRM-MS data was successfully collected on 55 clinical samples and 10 HPS samples while seven samples were lost during processing (Table 2). The normalization factors of the six normalizing proteins had a median coefficient of variation (CV) of 20.4% as evaluated from individual samples in the study.





TABLE 3



Clinical information of patients in Study I

















Nodule




Disease
Cancer



Size
Smoking
Pack-

Subject
Status
Stage
Histopathology of Cancer
Gender
Age
(mm)
Status
Year













S1
Benign


Male
69
33
Past
98

S2
Cancer
I
ADENOCARCINOMA
Male
71
6
Past
25




w/MIXED






S3
Cancer
I
SQUAMOUS
Male
53
45
Past
114

S4
Benign


Male
56
26
Past
10

S5
Benign


Female
63
35
Never
0

S6
Cancer
I
ADENOCARCINOMA
Female
64
39
Past
7.5

S7
Benign


Male
67
47.5
Current
55

S8
Cancer
II
SQUAMOUS
Male
72
43
Past
48

S9
Benign


Male
58
25
Current
40

S10
Cancer
I
ADENOCARCINOMA
Male
57
30
Past
40

S11
Cancer
I
ADENOCARCINOMA
Male
61
36
Past
20

S12
Benign


Male
55
22
Past
15

S13
Benign


Female
52
43
Never
0

S14
Benign


Female
46
40
Current
30

S15
Benign


Male
52
38
Past
9

S16
Cancer
II
MIXED
Male
48
40
Current
35

S17
Benign


Male
64
32
Past
12

S18
Cancer
II
SQUAMOUS
Male
62
24
Current
92

S19
Cancer
I
SQUAMOUS
Male
73
32
Past
130

S20
Benign


Male
74
22
Never
0

S21
Cancer
II
ADENOCARCINOMA
Female
67
40
Past
40

S22
Benign


Female
66
40
Never
0

S23
Cancer
I
ADENOCARCINOMA
Male
60
12
Past
17.5

S24
Benign


Male
62
40
Past
NA

S25
Benign


Female
76
40
Never
0

S26
Cancer
II
ADENOCARCINOMA
Female
75
44
Past
58

S27
Benign


Female
62
23
Past
NA

S28
Cancer
I
ADENOCARCINOMA
Female
58
22
Current
19

S29
Benign


Male
43
56
Current
30

S30
Cancer
II
SQUAMOUS
Male
51
58
Past
95

S31
Benign


Male
77
26
NA
45.6

S32
Cancer
I
ADENOCARCINOMA
Male
75
9
Past
91.5

S33
Cancer
II
ADENOCARCINOMA
Male
79
50
Past
12.5

S34
Benign


Male
78
48
Past
NA

S35
Cancer
II
NA
Female
56
24
Past
36

S36
Benign


Female
59
24
Past
47.5

S37
Benign


Female
64
26
Never
0

S38
Cancer
I
ADENOCARCINOMA
Female
64
33
Never
0

S39
Benign


Female
53
54
Past
30

S40
Cancer
II
ADENOCARCINOMA
Female
72
32
Current
50

S41
Benign


Female
61
35
Past
25

S42
Cancer
II
ADENOCARCINOMA
Female
60
25
Current
90




w/MIXED






S43
Benign


Female
38
22
Never
0

S44
Cancer
II
ADENOCARCINOMA
Female
75
45
Current
12.5

S45
Benign


Male
65
24
Past
30

S46
Cancer
I
ADENOCARCINOMA
Male
64
50
Past
NA

S47
Cancer
II
SQUAMOUS
Female
47
55
Never
0

S48
Benign


Female
52
24
Past
15

S49
Cancer
IIA
SQUAMOUS
Male
82
34
Past
30

S50
Benign


Male
84
32
Past
30

S51
Cancer
II
SQUAMOUS
Female
63
47
Past
15

S52
Benign


Male
44
45
Past
30

S53
Benign


Male
41
31
Past
20

S54
Cancer
II
ADENOCARCINOMA
Male
51
35
Past
30

S55
Cancer
II
SQUAMOUS
Female
67
52
Current
NA

S56
Benign


Female
65
42
Never
0

S57
Benign


Female
54
5
Past
20

S58
Cancer
II
ADENOCARCINOMA
Female
67
44
Past
30




w/MIXED






S59
Benign


Female
60
41
Current
40

S60
Cancer
II
ADENOCARCINOMA
Female
61
27
Past
40
























































































To understand how SISQuan and EPN controlled technical variability, intensity drift was defined as the relative deviation of protein intensity in individual samples from the corresponding median intensity in all samples, and was evaluated based on data of the 10 HPS samples ( FIG. 3 ), using the four quantification methods. Since the 10 HPS samples were identical, the deviation of protein drifts from zero represented the analytical variability in the experiment. The mean of protein drifts, plotted as a solid line in the inserts of FIG. 3 , measured the strength of variation that affected all proteins similarly (i.e., the strength of systematic variation). The 95% confidence interval (CI) of protein drifts, plotted as a shaded band in the inserts of FIG. 3 , measured the strength of variation that affected different proteins differently (i.e., the strength of protein-specific variation). In comparison with the protein drifts for the raw MS data ( FIG. 3A ), the protein drifts for EPN had a lower absolute mean but a comparable 95% CI ( FIG. 3B ) while the protein drifts for SISQuan had a lower 95% CI but a comparable absolute mean ( FIG. 3C ). Thus, EPN effectively controlled systematic variation and SISQuan effectively controlled protein-specific variation, illustrating the complementary nature of the two methods. The protein drifts for InteQuan had a lower absolute mean and a lower 95% CI ( FIG. 3D ), illustrating that InteQuan suppressed both systematic and protein-specific variation.
A. Intensity Drift
The intensity drift of peptide p in sample s was defined as

D p,s =( I p,s +{hacek over (I)} p )/ {hacek over (I)} p (EQN. 3)

where I p,s was the abundance of the peptide in the sample and {hacek over (I)} p was the corresponding median value in all technical replica. The intensity drift D p,s evaluated how far the abundance of the peptide in the sample deviated from the overall median abundance of the peptide. The median value of D p,s was zero by definition for all peptides.

Example 3
Improvement on Precision of Protein Measurement
To assess the precision of InteQuan and SISQuan, CVs of the target proteins were evaluated from the 10 HPS samples (Table 4). InteQuan had better precision than SISQuan on all proteins except for ISLR. The median CV of all proteins was 9.3% using InteQuan versus 13.3% using SISQuan. InteQuan was statistically more precise than SISQuan (P=5.2×10 −4 ) and lowered protein CV by a median value of 4.9%. Using InteQuan, the highest CV was 16.8% (FRIL, 12 ng/ml). CVs of the remaining 15 target proteins were all below 15%, including 10 proteins with a CV below 10% and two proteins with a CV at or below 5%.





TABLE 4



Coefficient of variation (CV) of protein abundance as evaluated

using InteQuan and using SISQuan











Study I
Study II
Study III


CV (%)
CV (%)
CV (%)









Protein
InteQuan
SISQuan
InteQuan
SISQuan
InteQuan
SISQuan



KIT
8.8
9.6
4.3
10.5
7.8
35.0

FRIL
16.8
25.2
5.3
11.8
7.3
33.6

COIA1
12.6
17.8
9.0
9.7
10.4
38.5

PRDX1
10.7
15.4
4.9
11.8
9.6
32.2

TENX
11.1
13.6
8.8
12.7
10.7
26.6

ENPL
13.1
18.9
11.2
8.6
11.3
34.9

GRP78
6.8
11.4
19.3
24.2
11.5
33.7

BGH3
5.0
12.3
5.7
9.9
12.4
42.3

ALDOA
6.6
13.6
9.1
17.5
15.1
35.4

GGH
6.9
7.1
9.0
13.7
13.9
38.9

CD14
4.1
8.0
4.6
12.0
4.6
35.3

LG3BP
8.8
13.0
5.9
10.0
5.6
31.1

TSP1
11.6
18.3
12.9
17.0
21.5
45.7

IBP3
5.7
11.6
6.3
13.5
13.5
41.5

TETN
9.9
17.8
9.9
12.4
25.5
52.3

ISLR
10.0
9.1
4.7
10.0
4.2
31.3

Median
9.3
13.3
7.6
11.9
11.0
35.1

CV (%)







Median
4.9

4.8

25.9


of CV







reduction a







(%)







Total
15
1
15
1
16
0

proteins







with







lower







CV












P value
5.2 × 10 −4
5.2 × 10 −5
3.1 × 10 −4

(paired




sign




test)




Comments
CV of 10 HPS
Median CV of 15
Median CV of 6 samples


aliquots
clinical samples








Study III
Combined HPS dataset







Generalized CV (%)
CV (%)
Generalized CV (%)











Protein
InteQuan
SISQuan
InteQuan
SISQuan
InteQuan
SISQuan




KIT
7.7
6.5
10.5
31.7
10.6
12.8


FRIL
7.3
12.0
15.9
25.4
16.7
18.8


COIA1
10.2
14.1
16.1
31.1
16.6
17.5


PRDX1
9.4
14.6
31.9
15.8
31.8
15.5


TENX
10.8
11.1
10.7
32.3
10.6
15.2


ENPL
11.3
13.6
11.7
32.9
11.8
12.6


GRP78
11.2
12.9
13.9
28.5
13.6
18.4


BGH3
12.4
13.8
9.8
33.7
9.6
14.9


ALDOA
14.9
19.9
11.0
35.4
11.2
14.9


GGH
13.7
16.5
11.0
31.7
10.5
16.6


CD14
4.7
6.8
7.9
30.5
7.8
11.9


LG3BP
5.6
7.5
8.6
30.5
8.5
13.9


TSP1
20.9
21.9
20.3
41.1
22.9
25.3


IBP3
14.0
15.5
19.8
26.7
20.5
21.6


TETN
26.9
29.7
33.1
47.6
33.3
37.7


ISLR
4.2
5.0
12.2
33.8
12.1
9.6


Median
11.0
13.7
11.9
31.7
11.9
15.3


CV (%)








Median
2.0

21.0

3.0



of CV








reduction a








(%)








Total
15
1
15
1
14
2


proteins








with








lower








CV














P value
5.2 × 10 −4
5.2 × 10 −4
4.2 × 10 −3


(paired





sign





test)










Comments
Median CV of
CV of 29 HPS aliquots



6 samples



a CV reduction was defined as CV using SISQuan minus CV using InteQuan.




















































































































Example 4
Improvement on Panel Performance in Disease Diagnosis
To illustrate the benefit of using InteQuan in clinical applications, a protein panel was constructed of all 16 target proteins and tested on the clinical samples in Study I using Monte Carlo cross validation (MCCV). Since the sample size was very small, the panel was not optimized for intended use, owing to concerns on both high false positive rate and high false negative rate. Using either InteQuan or SISQuan, the performance of the panel was summarized by the two receiver operating characteristic (ROC) curves in FIG. 4 . The corresponding AUC was 0.573 (95% CI 0.569-0.576) using InteQuan or 0.528 (95% CI 0.524-0.532) using SISQuan, respectively. The improvement by InteQuan was 0.045 (95% CI 0.042-0.048, P<0.0001). Thus, the panel had a significantly better performance using InteQuan than using SISQuan. More importantly, the ROC curve using InteQuan was consistently better than the ROC curve using SISQuan everywhere: See FIG. 4 . This comparative analysis demonstrated that InteQuan improved the performance of the 16-protein panel in disease diagnosis, illustrating its value for biomarker research, despite the fact that the panel was not optimized for clinical application. A protein panel comprising a subset of the 16 target proteins was recently optimized and validated using the InteQuan quantification method (Vachani et al., Validation of a Multi-Protein Plasma Classifier to Identify Benign Lung Nodules, J. of Thoracic. Onc. (doi: 10.1097/JTO.0000000000000447)).
A. Monte Carlo Cross Validation
Monte Carlo cross validation (MCCV) was performed as follows: First, all clinical samples in Study I were randomly assigned to a training group (including 24 benign and 24 cancer samples) or a test group (including 4 benign and 3 cancer samples). Second, two logistic regression models were developed to fit the disease status of the training samples, using either the InteQuan abundances or the SISQuan abundances of all 16 proteins in Table 4 as predictors. The first two steps were repeated if any one of the two models failed to converge. Third, the models were used to calculate scores of the test samples, evaluating their likelihood of being a cancer sample, based on protein InteQuan or SISQuan abundances, respectively. Fourth, the test samples were ranked by their scores from the InteQuan model or the SISQuan model, respectively. Fifth, the first four steps were repeated 10,000 times with different sample permutations. The ranking and the corresponding disease status of the test samples in all permutations were assembled under either InteQuan or SISQuan, respectively. Finally, comparison of ROC curves was carried out to compare the MCCV performance of the 16-protein panel using InteQuan with the corresponding performance using SISQuan. Due to small sample size, covariates in both the training samples and the test samples were unavoidable and difficult to adjust, which made it not meaningful to direct compare scores of the test samples between different permutations. Thus, the ranking instead of the score was combined for the ROC analysis, which effectively standardized the scores between different permutations. The ROC comparison analysis was performed by MedCalc (Ostend, Belgium), selecting “DeLong et al.” and “Binomial exact Confidence Interval for the AUC” as options.
Example 5
Better Control of Analytical Variability
To determine whether InteQuan can better control analytical variability during use of multiple depletion columns on clinical samples, 18 clinical samples in triplicate along with 12 aliquots of the HPS sample were analyzed in three experimental batches using three depletion columns in Study II ( FIG. 1C and Table 2). The three aliquots of the clinical samples were processed either using different depletion columns or using the same column but at different positions in the depletion sequence, monitoring analytical variability due to column or position difference. Out of the 66 samples, an HPS sample and three clinical samples were lost during processing (Table 2). A new lot of SIS peptide mixture was used in this study. The correction factors between the new and the old lots of SIS peptide mixture were determined from a migration experiment and are listed in Table 1. This dataset to compare different quantification methods.
The median CVs of the target proteins were evaluated from the 15 clinical samples having three replicate measurements (Table 4). InteQuan demonstrated better precision than SISQuan on all proteins except for ENPL. The median CV of all proteins was 7.6% using InteQuan versus 11.9% using SISQuan. InteQuan was statistically more precise than SISQuan (P=5.2×10 −4 ) and lowered protein CV by a median value of 4.8%. Using InteQuan, the highest CV was 19.3% (GRP78, 100 ng/ml). CVs of the remaining 15 target proteins were all below 15%, including 13 proteins with a CV below 10% and four proteins with a CV below 5%.
To assess whether InteQuan can better control analytical variability without compromising its ability to reveal biological difference among the clinical samples, principal variance component analysis (PVCA) was carried out to identify the major sources of variation in the experiment, including biological variation among individual patients (denoted as “patient”), analytical variation among depletion columns (denoted as “column”), and analytical variation among positions within a depletion sequence (denoted as “position”) ( FIG. 5 ). For InteQuan, “patient” alone contributed 97.6% to the total variability while other sources jointly contributed a negligible fraction of 2.4%. For SISQuan, “patient” alone contributed 87.2% to the total variability while other sources jointly contributed 12.8%. Thus, InteQuan enhanced the ability of measuring biological difference among the clinical samples, in agreement with the previous observation that InteQuan improved the performance of the 16-protein panel in Study I. In other words, InteQuan improved the resolution of protein measurement in clinical samples.
Example 6
High Tolerance Against Variation in Total Protein Concentration
To demonstrate InteQuan's ability to handle the variation in the total protein concentration, six samples were diluted into three concentrations (equivalent to the loading of 15, 20, or 30 μl of the original plasma samples: see Example 1, C.) and analyzed in duplicate using two depletion columns along with eight aliquots of the HPS sample in Study III ( FIG. 1C and Table 2). No data was collected on one of the 44 samples owing to sample exhaustion. Manual review of experimental data identified two erratic measurements ( FIG. 6A , B) that were eliminated from further analysis.
The median CVs of the target proteins were evaluated from the six samples using all valid measurements (Table 4). The median CV of all proteins was 11.0% using InteQuan and 35.1% using SISQuan. As a reference, the CV evaluated from the equivalent loading volumes (duplicates of 15, 20, and 30 μl) was 31.5%. While the median CV using SISQuan was higher than the CV of the loading volume, the median CV using InteQuan was much lower. Using InteQuan, all proteins had a median CV less than 20% except for TETN (25.5%) and TSP1 (21.5%), despite a two-fold difference in the total protein concentration.
Example 7
Usage of Generalized CV for Precision Evaluation in Study III
The high CVs of the target proteins using SISQuan in Study III reflected the large difference in the total protein concentration ( FIG. 6 ) rather than the precision of SISQuan. To compare the precision of InteQuan and SISQuan, a generalized method for CV calculation was developed. This method included two steps: First, the abundance of proteins in a sample was modeled either as linear functions of the loading volume (SISQuan) or as constants independent of the loading volume (InteQuan). Second, error propagation theory was applied to calculate the generalized CV as the standard deviation of differences between the modeled and the experimental abundances after logarithmic transformation. The modeled and the experimental abundances of all proteins in all samples collapsed nicely onto the respective diagonal line in FIGS. 7A and B, indicating that the method worked very well for both InteQuan and SISQuan. For SISQuan, it also demonstrated that proteins were measured within the respective linear dynamic range of the assays at all three concentrations. The generalized CVs and the standard CVs of InteQuan abundance were almost identical for all proteins in all samples ( FIG. 7C ). On the contrary, the generalized CVs of SISQuan abundance were uniformly lower than the corresponding standard CVs ( FIG. 7D ).
The median generalized CVs of the target proteins were evaluated from the six samples using all valid measurements (Table 4). InteQuan demonstrated better precision than SISQuan on all proteins except for KIT. The median generalized CV of all proteins was 11.0% using InteQuan versus 13.7% using SISQuan. InteQuan was statistically more precise than SISQuan (P=5.2×10 −4 ) and lowered protein generalized CV by a median value of 2.0%.
The generalized CV can be applied to analyze data from dilution experiments within the linear dynamic range and to provide an assessment on precision over the whole concentration range. Ideally, generalized CV should be evaluated on data covering three or more concentrations to avoid over-fitting.
A. Calculation of Generalized CV
The method consisted of two steps:
In the first step, protein InteQuan abundances in a sample were modeled as constants independent of the loading volume. Thus, the expected InteQuan abundances were assigned to the corresponding average values, that is









R
~


p
,
s

I

=


1

N

p
,
s



⁢


∑

i
=
1


N

p
,
s



⁢



⁢



R
~


p
,
s
,
i


.







(

EQN
.



⁢
4

)







Here {tilde over (R)} p,s,i was the InteQuan abundance of protein p in sample s at the loading volume v i =15, 20, or 30 and N p,s was the number of repeat measurements of the protein on the sample regardless of the loading volumes, that is N p,s =5 or 6 in Study III and N p,s =29 for the combined HPS dataset.
On the contrary, protein SISQuan abundances in a sample were modeled as linear functions of the loading volume. More specifically, the expected SISQuan abundances were fitted as linear functions of the loading volume such that

R p,s,i S =a p,s *v i +b p . (EQN. 5)

Here a p,s was proportional to the concentration of the protein in the sample and b y was common to all samples. Parameters {a p,s } and b p were evaluated from repeat measurements of the protein in all samples using maximum likelihood estimation. Ideally one should have at least three loading volumes to avoid over-fitting.
In the second step, error propagation theory was applied to evaluate the generalized CV. According to the theory, the CV of a quantity equals to the standard deviation of the same quantity after logarithmic transformation, that is
CV ⁡ ( x ) = σ ⁡ ( x ) x = σ ⁡ ( ln ⁡ ( x ) )
where σ(x) represents the standard deviation of x and ln(x) is the natural logarithmic function. Thus, the generalized CV of protein abundance was evaluated from differences between the expected and the experimental values after logarithmic transformation. More specifically, the generalized CV of InteQuan abundance was evaluated as:









CV

p
,
s

I

=





∑

i
=
1


N

p
,
s



⁢



⁢


[


ln
⁡

(


R
~


p
,
s
,
i


)


-

ln
⁡

(


R
~


p
,
s

I

)



]

2




N

p
,
s


-
1



.





(

EQN
.



⁢
6

)







And the generalized CV of SISQuan abundance was evaluated as:








CV

p
,
s

S

=





∑

i
=
1


N

p
,
s



⁢



⁢


[


ln
⁡

(

R

p
,
s
,
i


)


-

ln
⁡

(

R

p
,
s
,
i

S

)



]

2




N

p
,
s


-
1
-

1
/
K




.





(

EQN
.



⁢
7

)







Here K was the number of different samples used in the study and was needed to account for the fitting of b p . Thus, K=6 in Study III and K=1 for the combined HPS dataset.
Example 8
Robustness of the Depletion-MRM-MS Platform
The three assessment studies lasted over six months, were carried out by different operators, encountered major instrument repairs, required implementation of a protocol change in sample loading volume, and used different reagent lots (Table 5).





TABLE 5



Summary of main experimental differences among the three assessment

studies and major instrument services







Procedure
Description
Study I
Study III
Study II



Depletion
Date
18-21 Nov 2012
21-23 Feb 2013
15-17 Apr 2013


Operator
JM
JM
JM


Column lot
20093648
20093648 &
3315179




SLBD7383


HPLC machine
HPLC-LC00001
HPLC-LC00001
HPLC-LC-





00001: Batch 1





& 2





HPLC-LC-





00002: Batch 3


HPLC service a
18 Nov 2012
20 & 22 Feb 2013
HPLC-





LC00001: 26,





28 & 29 Mar;





15 & 17 Apr





2013





HPLC-





LC00002: 17





Apr 2013

Digestion
Date
26-27 Nov 2012
28 Feb-01 Mar
23-24 Apr 2013




2013


Operator
DAS
JM
DAS/JM

Solid-phase
Date
3-4 Dec 2012
03-04 Mar 2013
25-26 Apr 2013

extraction
Operator
DAS
JM
DAS/JM

MS analysis
Date
11 Dec 2012-9
6-10 Mar 2013
29 Apr-5 May



Jan 2013

2013


Operator
JT
DAS
DAS


SIS lot
Lot 1
Lot 1
Lot 2


RP column lot
Not Recorded
S/N USCGC02709
S/N





USCGC02348






RP HPLC service
HPLC Flush (17 Jan), analytical column change


in Year
(14 Feb 21 & 25 Mar), B side check valve replaced


2013
(15 Feb), A side pump head replaced (19 Feb), B side pump



head replaced (4 Mar), and PM service on LC unit



(11 Mar).


MS service
Hexa bore capillary change (18 Dec 2012, 6 Feb 2013,



4 Mar 2013) and Ion funnels cleaned & PM on MS



(11 Mar 2013).








Data curation
JT
DAS
DAS





a Changed depletion column, purge valve frit, 2 inline filters, and needle seat.































































To assess the robustness of the depletion-MRM-MS platform, the CVs and the generalized CVs of the target proteins were computed from the 29 HPS samples across all three studies ( FIG. 7E-H and Table 4). Using InteQuan, 13 of the 16 target proteins had a CV less than 20%, 10 had a CV less than 15%, and three had a CV less than 10%. Only three proteins had a CV greater than 20%, including TSP1 with a CV of 20.3%, PRDX1 with a CV of 31.9%, and TETN with a CV of 33.1%.
We investigated possible causes for the high CVs of PRDX1 and TETN. On PRDX1, we noticed that, despite a lower loading volume per sample in Study II that was only two thirds of the loading volume in the other two studies, its SISQuan abundance was almost the same in all three studies. As a result, its InteQuan abundance was about 77% higher in Study II than in the other two studies, which led to the large CV value, The CV of its EPN abundance was only 16.7% Thus, the large CV of its InteQuan abundance was likely due to issues on isotopic labeling rather than protein normalization. Possible causes for the inflated PRDX1 abundance in Study II include: (i) the correction factor for PRDX1 in Table 1 was incorrectly determined and/or (ii) the SIS peptide of QITVNDLPVGR (SEQ ID NO: 36) of PRDX1. was partially cyclized [31]in Study II. On TETN, we noticed that the generalized CV of its SISQuan abundance was even higher at 37.7%. It turns out that TETN partially binds to the IgY14-Supermix resin column [25]. Possible causes for the large CV of TETN include: (i) the binding affinity varied between different depletion columns and/or (ii) the binding affinity was sensitive to the loading volume. In both cases InieQuan, as a quantification method itself, was not the cause for the high CV values.
Based on generalized CV, InteQuan had better precision than SISQuan on all proteins except for PRDX1 and ISLR. The median generalized CV of all proteins was 11.9% using InteQuan versus 15.3% using SISQuan. The generalized CV likely overestimated the precision of SISQuan since linear functions were used to fit SISQuan abundances at only two different protein concentrations, instead of the desirable three or more concentrations to avoid over-fitting. Nevertheless, InteQuan was statistically more precise than SISQuan (P=4.2×10 −3 ) and lowered protein generalized CV by a median value of 3.0%. Based on standard CV, the superiority of InteQuan to SISQuan was even more significant (P=5.2×10 −4 ).
Example 9
Evaluation of Alternative Normalizing Panels
To assess whether a subset of the six normalizing proteins can form a normalization panel that is better than the full panel, all combinations of the six normalizing proteins were used to quantify the target proteins using InteQuan on the combined HPS dataset ( FIG. 8A ). The “best” panel, selected on the basis that the corresponding median CV of all proteins was the lowest, consisted of PEDF, MASP1, GELS and C163A. The “best” panel led to better precision than the full panel on 11 out of the 16 target proteins. The median CV of all proteins was 11.2% using the “best” panel versus 11.9% using the full panel. The “best” panel was marginally more precise than the full panel (P=0.21) and lowered protein CV by a median value of 0.2%.
To assess whether normalization panels that are tailored to individual target proteins are better than the full panel, an “optimal” panel was selected for each of the 16 target proteins on the basis that the corresponding CV was the lowest ( FIG. 8B ). To avoid over-fitting, the “optimal” panels were selected from the same dataset in which the scaling constants of InteQuan were determined and tested on the combined HPS dataset. The “optimal” panels led to better precision than the full panel on only 6 out of the 16 target proteins on the combined HPS dataset. The median CV of all proteins was 13.4% using the “optimal” panels versus 11.9% using the full panel. The “optimal” panels were marginally less precise than the full panel (P=0.45) and increased protein CV by a median value of 0.5%. The full panel had a better overall performance than “optimal” panels on the combined HPS dataset.
Example 10
EPN as an Alternative to InteQuan
To compare the precision of all four quantification methods (raw MS data, EPN, SISQuan, and InteQuan), CVs of protein abundance were evaluated from data of the 10 HPS samples in Study I (Table 6). Among the four methods, InteQuan was statistically more precise than SISQuan (P=5.2×10 −4 ), SISQuan was only marginally better than EPN (P=0.80), and EPN was significantly better than the raw data (P=3.1×10 −5 ). Thus, the four quantification methods were ranked by their precision in descending order as InteQuan, SISQuan, EPN, and the raw MS data.





TABLE 6



Coefficient of variation (CV) of protein abundance as evaluated on the 10 HPS

samples in Study I





CV reduction (%)







Protein
CV (%)
SISQuan −
EPN −











(HUMAN)
InteQuan
SISQuan
EPN
Raw
InteQuan
SISQuan
Raw − EPN












KIT
8.8
9.6
15.0
24.1
0.8
5.4
9.2

FRIL
16.8
25.2
21.9
28.0
8.4
−3.3
6.1

COIA1
12.6
17.8
20.3
27.2
5.2
2.5
7.0

PRDX1
10.7
15.4
12.1
19.7
4.7
−3.3
7.6

TENX
11.1
13.6
9.6
20.2
2.5
−4.0
10.6

ENPL
13.1
18.9
6.9
19.1
5.8
−12.0
12.2

GRP78
6.8
11.4
16.7
22.5
4.5
5.3
5.8

BGH3
5.0
12.3
16.7
23.3
7.3
4.4
6.5

ALDOA
6.6
13.6
17.7
28.6
7.0
4.1
11.0

GGH
6.9
7.1
6.8
17.2
0.2
−0.3
10.3

CD14
4.1
8.0
21.1
31.8
3.8
13.1
10.7

LG3BP
8.8
13.0
8.8
19.4
4.2
−4.2
10.6

TSP1
11.6
18.3
19.6
22.5
6.7
1.3
2.8

IBP3
5.7
11.6
5.7
12.4
5.9
−5.9
6.6

TETN
9.9
17.8
19.6
21.7
7.9
1.8
2.1

ISLR
10.0
9.1
13.5
20.5
−0.9
4.4
7.0

Median
9.3
13.3
15.9
22.1
4.9
1.6
7.3






P value (paired sign test)
5.2 × 10 −04
0.80
3.1 × 10 −05





















































The above results indicate that the label-free EPN could be an attractive alternative to InteQuan, especially when it was too costly to obtain SIS peptides for hundreds to thousands of proteins of interest in early-stage biomarker discovery studies. Using EPN, the median CV of all proteins was 15.9%. Three proteins had a CV just above 20%, including FRIL (12 ng/ml, 21.9%), CD14 (420 ng/ml, 21.1%), and COIA1 (35 ng/ml, 20.3%). CVs of the remaining 13 target proteins were all below 20%, including eight proteins with a CV at or below 15% and five proteins with a CV below 10%.
To further assess EPN, Pearson correlation coefficients of protein abundance as evaluated using different quantification methods were calculated on data of the 55 clinical samples in Study I (Table 7). The median Pearson correlation coefficient between InteQuan and EPN was 0.843. The lowest coefficient between them was 0.621 (CD14, P=4.3×10 −7 ). So the correlation between InteQuan and EPN was significant for all the target proteins.





TABLE 7



Pearson correlation coefficient of protein abundance as evaluated on the 55

clinical samples in Study I










InteQuan
InteQuan
InteQuan
EPN



Protein
vs.
vs.
vs.
vs.
EPN vs.
SISQuan vs.

(HUMAN)
EPN
SISQuan
Raw
SISQuan
Raw
Raw



KIT
0.789
0.669
0.502
0.630
0.625
0.883

FRIL
0.963
0.919
0.862
0.933
0.898
0.971

COIA1
0.801
0.735
0.630
0.712
0.784
0.888

PRDX1
0.965
0.979
0.958
0.977
0.986
0.990

TENX
0.818
0.811
0.663
0.641
0.742
0.805

ENPL
0.936
0.883
0.838
0.781
0.831
0.893

GRP78
0.850
0.802
0.630
0.741
0.774
0.841

BGH3
0.740
0.679
0.589
0.685
0.760
0.882

ALDOA
0.954
0.958
0.943
0.902
0.927
0.977

GGH
0.837
0.804
0.749
0.621
0.844
0.792

CD14
0.621
0.498
0.234
0.704
0.727
0.782

LG3BP
0.900
0.910
0.826
0.807
0.865
0.913

TSP1
0.972
0.951
0.954
0.928
0.945
0.992

IBP3
0.918
0.816
0.749
0.756
0.784
0.872

TETN
0.775
0.779
0.702
0.680
0.745
0.905

ISLR
0.737
0.637
0.518
0.634
0.678
0.890

Median
0.843
0.808
0.725
0.726
0.784
0.889








































All the 55 clinical samples in Study I had matching data from a previous label-free study. Major differences between the two studies were described in Example 1, N. Pearson correlation coefficients of protein abundance were computed on data from the 55 clinical samples, using InteQuan on data from Study I and using EPN on data from the discovery study (Table 8). The median Pearson correlation coefficient between the two studies was 0.821. All proteins had a correlation coefficient above 0.5 except for TETN (0.418, P=1.5×10 −3 ). Despite major differences between the two studies, the correlation between EPN and InteQuan was significant for all the target proteins. Correlations between all feasible quantification methods on the two datasets are also listed in Table 8. Based on this evidence, it was justified to use EPN as an economical alternative to InteQuan in early-stage biomarker discovery studies.





TABLE 8



Pearson correlation coefficient of protein abundance between Study I and a

discovery study a












InteQuan
InteQuan
EPN
EPN
SISQuan
SISQuan
Raw
Raw

Protein
vs.
vs.
vs.
vs.
vs.
vs.
vs.
vs.

(HUMAN)
EPN
Raw
EPN
Raw
EPN
Raw
EPN
Raw



KIT
0.711
0.536
0.560
0.404
0.270
0.434
0.196
0.290

FRIL
0.953
0.829
0.850
0.781
0.815
0.857
0.721
0.791

COIA1
0.770
0.610
0.715
0.605
0.679
0.735
0.626
0.646

PRDX1
0.978
0.971
0.946
0.943
0.950
0.975
0.937
0.956

TENX
0.831
0.690
0.800
0.711
0.607
0.683
0.551
0.615

ENPL
0.648
0.652
0.629
0.638
0.499
0.647
0.432
0.560

GRP78
0.649
0.618
0.620
0.656
0.444
0.624
0.367
0.597

BGH3
0.521
0.364
0.216
0.202
0.284
0.523
0.204
0.363

ALDOA
0.900
0.882
0.868
0.847
0.865
0.899
0.826
0.870

GGH
0.835
0.622
0.830
0.666
0.656
0.717
0.745
0.755

CD14
0.841
0.395
0.588
0.412
0.543
0.642
0.330
0.411

LG3BP
0.921
0.836
0.833
0.772
0.904
0.935
0.854
0.911

TSP1
0.909
0.802
0.876
0.750
0.918
0.902
0.918
0.884

IBP3
0.811
0.664
0.750
0.598
0.515
0.609
0.457
0.486

TETN
0.418
0.416
0.277
0.289
0.353
0.551
0.443
0.581

ISLR
0.783
0.651
0.700
0.637
0.554
0.712
0.500
0.624

Median
0.821
0.652
0.733
0.647
0.581
0.698
0.526
0.620



a Evaluated on the 55 common clinical samples between the two studies and labeled as method on data of Study I versus method on data of the discovery study.







































The amino acid sequences for all proteins of interest described herein are described in Table 9.





TABLE 9



Amino Acid Sequences for Proteins of Interest





Protein

Seq.

Name
Amino Acid Sequences
ID.



BGH3_
MALFVRLLALALALALGPAATLAGPAKSPYQLVLQHSRLRGRQHGPNVCAVQKVIGTNRKYFTNCKQWYQRKICGKSTVISYECC
1

HUMAN
PGYEKVPGEKGCPAALPLSNLYETLGVVGSTTTQLYTDRTEKLRPEMEGPGSFTIFAPSNEAWASLPAEVLDSLVSNVNIELLNA



LRYHMVGRRVLTDELKHGMTLTSMYQNSNIQIHHYPNGIVTVNCARLLKADHHATNGVVHLIDKVISTITNNIQQIIEIEDTFET



LRAAVAASGLNTMLEGNGQYTLLAPTNEAFEKIPSETLNRILGDPEALRDLLNNHILKSAMCAEAIVAGLSVETLEGTTLEVGCS



GDMLTINGKAIISNKDILATNGVIHYIDELLIPDSAKTLFELAAESDVSTAIDLFRQAGLGNHLSGSERLTLLAPLNSVFKDGTP



PIDAHTRNLLRNHIIKDQLASKYLYHGQTLETLGGKKLRVFVYRNSLCIENSCIAAHDKRGRYGTLFTMDRVLTPPMGTVMDVLK



GDNRFSMLVAAIQSAGLTETLNREGVYTVFAPTNEAFRALPPRERSRLLGDAKELANILKYHIGDEILVSGGIGALVRLKSLQGD



KLEVSLKNNVVSVNKEPVAEPDIMATNGVVHVITNVLQPPANRPQERGDELADSALEIFKQASAFSRASQRSVRLAPVYQKLLER



MKH




GGH_
MASPGCLLCVLGLLLCGAASLELSRPHGDTAKKPIIGILMQKCRNKVMKNYGRYYIAASYVKYLESAGARVVPVRLDLTEKDYEI
2

HUMAN
LFKSINGILFPGGSVDLRRSDYAKVAKIFYNLSIQSFDDGDYFPVWGTCLGFEELSLLISGECLLTATDTVDVAMPLNFTGGQLH



SRMFQNFPTELLLSLAVEPLTANFHKWSLSVKNFTMNEKLKKFFNVLTTNTDGKIEFISTMEGYKYPVYGVQWHPEKAPYEWKNL



DGISHAPNAVKTAFYLAEFFVNEARKNNHHFKSESEEEKALIYQFSPIYTGNISSFQQCYIFD




LG3BP_
MTPPRLFWVWLLVAGTQGVNDGDMRLADGGATNQGRVEIFYRGQWGTVCDNLWDLTDASVVCRALGFENATQALGRAAFGQGSGP
3

HUMAN
IMLDEVQCTGTEASLADCKSLGWLKSNCRHERDAGVVCTNETRSTHTLDLSRELSEALGQIFDSQRGCDLSISVNVQGEDALGFC



GHTVILTANLEAQALWKEPGSNVTMSVDAECVPMVRDLLRYFYSRRIDITLSSVKCFHKLASAYGARQLQGYCASLFAILLPQDP



SFQMPLDLYAYAVATGDALLEKLCLQFLAWNFEALTQAEAWPSVPTDLLQLLLPRSDLAVPSELALLKAVDTWSWGERASHEEVE



GLVEKIRFPMMLPEELFELQFNLSLYWSHEALFQKKTLQALEFHTVPFQLLARYKGLNLTEDTYKPRIYTSPTWSAFVTDSSWSA



RKSQLVYQSRRGPLVKYSSDYFQAPSDYRYYPYQSFQTPQHPSFLFQDKRVSWSLVYLPTIQSCWNYGFSCSSDELPVLGLTKSG



GSDRTIAYENKALMLCEGLFVADVTDFEGWKAAIPSALDTNSSKSTSSFPCPAGHFNGFRTVIRPFYLTNSSGVD




PRDX1_
MSSGNAKIGHPAPNFKATAVMPDGQFKDISLSDYKGKYVVFFFYPLDFTFVCPTEIIAFSDRAEEFKKLNCQVIGASVDSHFCHL
4

HUMAN
AWVNTPKKQGGLGPMNIPLVSDPKRTIAQDYGVLKADEGISFRGLFIIDDKGILRQITVNDLPVGRSVDETLRLVQAFQFTDKHG



EVCPAGWKPGSDTIKPDVQKSKEYFSKQK




TSP1_
MGLAWGLGVLFLMHVCGTNRIPESGGDNSVFDIFELTGAARKGSGRRLVKGPDPSSPAFRIEDANLIPPVPDDKFQDLVDAVRAE
5

HUMAN
KGFLLLASLRQMKKTRGTLLALERKDHSGQVFSVVSNGKAGTLDLSLTVQGKQHVVSVEEALLATGQWKSITLFVQEDRAQLYID



CEKMENAELDVPIQSVFTRDLASIARLRIAKGGVNDNFQGVLQNVRFVFGTTPEDILRNKGCSSSTSVLLTLDNNVVNGSSPAIR



TNYIGHKTKDLQAICGISCDELSSMVLELRGLRTIVTTLQDSIRKVTEENKELANELRRPPLCYHNGVQYRNNEEWTVDSCTECH



CQNSVTICKKVSCPIMPCSNATVPDGECCPRCWPSDSADDGWSPWSEWTSCSTSCGNGIQQRGRSCDSLNNRCEGSSVQTRTCHI



QECDKRFKQDGGWSHWSPWSSCSVTCGDGVITRIRLCNSPSPQMNGKPCEGEARETKACKKDACPINGGWGPWSPWDICSVTCGG



GVQKRSRLCNNPTPQFGGKDCVGDVTENQICNKQDCPIDGCLSNPCFAGVKCTSYPDGSWKCGACPPGYSGNGIQCTDVDECKEV



PDACFNHNGEHRCENTDPGYNCLPCPPRFTGSQPFGQGVEHATANKQVCKPRNPCTDGTHDCNKNAKCNYLGHYSDPMYRCECKP



GYAGNGIICGEDTDLDGWPNENLVCVANATYHCKKDNCPNLPNSGQEDYDKDGIGDACDDDDDNDKIPDDRDNCPFHYNPAQYDY



DRDDVGDRCDNCPYNHNPDQADTDNNGEGDACAADIDGDGILNERDNCQYVYNVDQRDTDMDGVGDQCDNCPLEHNPDQLDSDSD



RIGDTCDNNQDIDEDGHQNNLDNCPYVPNANQADHDKDGKGDACDHDDDNDGIPDDKDNCRLVPNPDQKDSDGDGRGDACKDDFD



HDSVPDIDDICPENVDISETDFRRFQMIPLDPKGTSQNDPNWVVRHQGKELVQTVNCDPGLAVGYDEFNAVDFSGTFFINTERDD



DYAGFVFGYQSSSRFYVVMWKQVTQSYWDTNPTRAQGYSGLSVKVVNSTTGPGEHLRNALWHTGNTPGQVRTLWHDPRHIGWKDF



TAYRWRLSHRPKTGFIRVVMYEGKKIMADSGPIYDKTYAGGRLGLFVFSQEMVFFSDLKYECRDP




CD44_
MDKFWWHAAWGLCLVPLSLAQIDLNITCRFAGVFHVEKNGRYSISRTEAADLCKAFNSTLPTMAQMEKALSIGFETCRYGFIEGH
6

HUMAN
VVIPRIHPNSICAANNTGVYILTSNTSQYDTYCFNASAPPEEDCTSVTDLPNAFDGPITITIVNRDGTRYVQKGEYRTNPEDIYP



SNPTDDDVSSGSSSERSSTSGGYIFYTFSTVHPIPDEDSPWITDSTDRIPATTLMSTSATATETATKRQETWDWFSWLFLPSESK



NHLHTTTQMAGTSSNTISAGWEPNEENEDERDRHLSFSGSGIDDDEDFISSTISTTPRAFDHTKQNQDWTQWNPSHSNPEVLLQT



TTRMTDVDRNGTTAYEGNWNPEAHPPLIHHEHHEEEETPHSTSTIQATPSSTTEETATQKEQWFGNRWHEGYRQTPKEDSHSTTG



TAAASAHTSHPMQGRTTPSPEDSSWTDFFNPISHPMGRGHQAGRRMDMDSSHSITLQPTANPNTGLVEDLDRTGPLSMTTQQSNS



QSFSTSHEGLEEDKDHPTTSTLTSSNRNDVTGGRRDPNHSEGSTTLLEGYTSHYPHTKESRTFIPVTSAKTGSFGVTAVTVGDSN



SNVNRSLSGDQDTFHPSGGSHTTHGSESDGHSHGSQEGGANTTSGPIRTPQIPEWLIILASLLALALILAVCIAVNSRRRCGQKK



KLVINSGNGAVEDRKPSGLNGEASKSQEMVHLVNKESSETPDQFMTADETRNLQNVDMKIGV




ENPL_
MRALWVLGLCCVLLTFGSVRADDEVDVDGTVEEDLGKSREGSRTDDEVVQREEEATQLDGLNASQIRELREKSEKFAFQAEVNRM
7

HUMAN
MKLIINSLYKNKEIFLRELISNASDALDKIRLISLTDENALSGNEELTVKIKCDKEKNLLHVTDTGVGMTREELVKNLGTIAKSG



TSEFLNKMTEAQEDGQSTSELIGQFGVGFYSAFLVADKVIVTSKHNNDTQHIWESDSNEFSVIADPRGNTLGRGTTITLVLKEEA



SDYLELDTIKNLVKKYSQFINFPIYVWSSKTETVEEPMEEEEAAKEEKEESDDEAAVEEEEEEKKPKTKKVEKTVWDWELMNDIK



PIWQRPSKEVEEDEYKAFYKSFSKESDDPMAYIHFTAEGEVTFKSILFVPTSAPRGLFDEYGSKKSDYIKLYVRRVFITDDFHDM



MPKYLNFVKGVVDSDDLPLNVSRETLQQHKLLKVIRKKLVRKTLDMIKKIADDKYNDTFWKEFGTNIKLGVIEDHSNRTRLAKLL



RFQSSHHPTDITSLDQYVERMKEKQDKIYFMAGSSRKEAESSPFVERLLKKGYEVIYLTEPVDEYCIQALPEFDGKRFQNVAKEG



VKFDESEKTKESREAVEKEFEPLLNWMKDKALKDKIEKAVVSQRLTESPCALVASQYGWSGNMERIMKAQAYQTGKDISTNYYAS



QKKTFEINPRHPLIRDMLRRIKEDEDDKTVLDLAVVLFETATLRSGYLLPDTKAYGDRIERMLRLSLNIDPDAKVEEEPEEEPEE



TAEDTTEDTEQDEDEEMDVGTDEEEETAKESTAEKDEL




TENX_
MMPAQYALTSSLVLLVLLSTARAGPFSSRSNVTLPAPRPPPQPGGHTVGAGVGSPSSQLYEHTVEGGEKQVVFTHRINLPPSTGC
8

HUMAN
GCPPGTEPPVLASEVQALRVRLEILEELVKGLKEQCTGGCCPASAQAGTGQTDVRTLCSLHGVFDLSRCTCSCEPGWGGPTCSDP



TDAEIPPSSPPSASGSCPDDCNDQGRCVRGRCVCFPGYTGPSCGWPSCPGDCQGRGRCVQGVCVCRAGFSGPDCSQRSCPRGCSQ



RGRCEGGRCVCDPGYTGDDCGMRSCPRGCSQRGRCENGRCVCNPGYTGEDCGVRSCPRGCSQRGRCKDGRCVCDPGYTGEDCGTR



SCPWDCGEGGRCVDGRCVCWPGYTGEDCSTRTCPRDCRGRGRCEDGECICDTGYSGDDCGVRSCPGDCNQRGRCEDGRCVCWPGY



TGTDCGSRACPRDCRGRGRCENGVCVCNAGYSGEDCGVRSCPGDCRGRGRCESGRCMCWPGYTGRDCGTRACPGDCRGRGRCVDG



RCVCNPGFTGEDCGSRRCPGDCRGHGLCEDGVCVCDAGYSGEDCSTRSCPGGCRGRGQCLDGRCVCEDGYSGEDCGVRQCPNDCS



QHGVCQDGVCICWEGYVSEDCSIRTCPSNCHGRGRCEEGRCLCDPGYTGPTCATRMCPADCRGRGRCVQGVCLCHVGYGGEDCGQ



EEPPASACPGGCGPRELCRAGQCVCVEGFRGPDCAIQTCPGDCRGRGECHDGSCVCKDGYAGEDCGEARVPSSASAYDQRGLAPG



QEYQVTVRALRGTSWGLPASKTITTMIDGPQDLRVVAVTPTTLELGWLRPQAEVDRFVVSYVSAGNQRVRLEVPPEADGTLLTDL



MPGVEYVVTVTAERGRAVSYPASVRANTEEREEESPPRPSLSQPPRRPWGNLTAELSRFRGTVQDLERHLRAHGYPLRANQTYTS



VARHIHEYLQRQVLGSSADGALLVSLDGLRGQFERVVLRWRPQPPAEGPGGELTVPGTTRTVSLPDLRPGTTYHVEVHGVRAGQT



SKSYAFITTTGPSTTQGAQAPLLQQRPQELGELRVLGRDETGRLRVVWTAQPDTFAYFQLRMRVPEGPGAHEEVLPGDVRQALVP



PPPPGTPYELSLHGVPPGGKPSDPIIYQGIMDKDEEKPGKSSGPPRLGELTVTDRTSDSLLLRWTVPEGEFDSFVIQYKDRDGQP



QVVPVEGPQRSAVITSLDPGRKYKFVLYGFVGKKRHGPLVAEAKILPQSDPSPGTPPHLGNLWVTDPTPDSLHLSWTVPEGQFDT



FMVQYRDRDGRPQVVPVEGPERSFVVSSLDPDHKYRFTLFGIANKKRYGPLTADGTTAPERKEEPPRPEFLEQPLLGELTVTGVT



PDSLRLSWTVAQGPFDSFMVQYKDAQGQPQAVPVAGDENEVTVPGLDPDRKYKMNLYGLRGRQRVGPESVVAKTAPQEDVDETPS



PTELGTEAPESPEEPLLGELTVTGSSPDSLSLFWTVPQGSFDSFTVQYKDRDGRPRAVRVGGKESEVTVGGLEPGHKYKMHLYGL



HEGQRVGPVSAVGVTAPQQEETPPATESPLEPRLGELTVTDVTPNSVGLSWTVPEGQFDSFIVQYKDKDGQPQVVPVAADQREVT



VYNLEPERKYKMNMYGLHDGQRMGPLSVVIVTAPATEASKPPLEPRLGELTVTDITPDSVGLSWTVPEGEFDSFVVQYKDRDGQP



QVVPVAADQREVTIPDLEPSRKYKFLLFGIQDGKRRSPVSVEAKTVARGDASPGAPPRLGELWVTDPTPDSLRLSWTVPEGQFDS



FVVQFKDKDGPQVVPVEGHERSVTVTPLDAGRKYRFLLYGLLGKKRHGPLTADGTTEARSAMDDTGTKRPPKPRLGEELQVTTVT



QNSVGLSWTVPEGQFDSFVVQYKDRDGQPQVVPVEGSLREVSVPGLDPAHRYKLLLYGLHHGKRVGPISAVAITAGREETETETT



APTPPAPEPHLGELTVEEATSHTLHLSWMVTEGEFDSFEIQYTDRDGQLQMVRIGGDRNDITLSGLESDHRYLVTLYGFSDGKHV



GPVHVEALTVPEEEKPSEPPTATPEPPIKPRLGELTVTDATPDSLSLSWTVPEGQFDHFLVQYRNGDGQPKAVRVPGHEEGVTIS



GLEPDHKYKMNLYGFHGGQRMGPVSVVGVTEPSMEAPEPAEEPLLGELTVTGSSPDSLSLSWTVPQGRFDSFTVQYKDRDGRPQV



VRVGGEESEVTVGGLEPGRKYKMHLYGLHEGRRVGPVSAVGVTAPEEESPDAPLAKLRLGQMTVRDITSDSLSLSWTVPEGQFDH



FLVQFKNGDGQPKAVRVPGHEDGVTISGLEPDHKYKMNLYGFHGGQRVGPVSAVGLTASTEPPTPEPPIKPRLEELTVTDATPDS



LSLSWTVPEGQFDHFLVQYKNGDGQPKATRVPGHEDRVTISGLEPDNKYKMNLYGFHGGQRVGPVSAIGVTEEETPSPTEPSMEA



PEPPEEPLLGELTVTGSSPDSLSLSWTVPQGRFDSFTVQYKDRDGRPQVVRVGGEESEVTVGGLEPGRKYKMHLYGLHEGRRVGP



VSTVGVTAPQEDVDETPSPTEPGTEAPGPPEEPLLGELTVTGSSPDSLSLSWTVPQGRFDSFTVQYKDRDGRPQAVRVGGQESKV



TVRGLEPGRKYKMHLYGLHEGRRLGPVSAVGVTEDEAETTQAVPTMTPEPPIKPRLGELTMTDATPDSLSLSWTVPEGQFDHFLV



QYRNGDGQPKAVRVPGHEDGVTISGLEPDHKYKMNLYGFHGGQRVGPISVIGVTEEETPSPTELSTEAPEPPEEPLLGELTVTGS



SPDSLSLSWTIPQGHFDSFTVQYKDRDGRPQVMRVRGEESEVTVGGLEPGRKYKMHLYGLHEGRRVGPVSTVGVTVPTTTPEPPN



KPRLGELTVTDATPDSLSLSWMVPEGQFDHFLVQYRNGDGQPKVVRVPGHEDGVTISGLEPDHKYKMNLYGFHGGQRVGPISVIG



VTEEETPAPTEPSTEAPEPPEEPLLGELTVTGSSPDSLSLSWTIPQGRFDSFTVQYKDRDGRPQVVRVRGEESEVTVGGLEPGCK



YKMHLYGLHEGQRVGPVSAVGVTAPKDEAETTQAVPTMTPEPPIKPRLGELTVTDATPDSLSLSWMVPEGQFDHFLVQYRNGDGQ



PKAVRVPGHEDGVTISGLEPDHKYKMNLYGFHGGQRVGPVSAIGVTEEETPSPTEPSTEAPEAPEEPLLGELTVTGSSPDSLSLS



WTVPQGRFDSFTVQYKDRDGQPQVVRVRGEESEVTVGGLEPGRKYKMHLYGLHEGQRVGPVSTVGITAPLPTPLPVEPRLGELAV



AAVTSDSVGLSWTVAQGPFDSFLVQYRDAQGQPQAVPVSGDLRAVAVSGLDPARKYKFLLFGLQNGKRHGPVPVEARTAPDTKPS



PRLGELTVTDATPDSVGLSWTVPEGEFDSFVVQYKDKDGRLQVVPVAANQREVTVQGLEPSRKYRFLLYGLSGRKRLGPISADST



TAPLEKELPPHLGELTVAEETSSSLRLSWTVAQGPFDSFVVQYRDTDGQPRAVPVAADQRTVTVEDLEPGKKYKFLLYGLLGGKR



LGPVSALGMTAPEEDTPAPELAPEAPEPPEEPRLGVLTVTDTTPDSMRLSWSVAQGPFDSFVVQYEDTNGQPQALLVDGDQSKIL



ISGLEPSTPYRFLLYGLHEGKRLGPLSAEGTTGLAPAGQTSEESRPRLSQLSVTDVTTSSLRLNWEAPPGAFDSFLLRFGVPSPS



TLEPHPRPLLQRELMVPGTRHSAVLRDLRSGTLYSLTLYGLRGPHKADSIQGTARTLSPVLESPRDLQFSEIRETSAKVNWMPPP



SRADSFKVSYQLADGGEPQSVQVDGQARTQKLQGLIPGARYEVTVVSVRGFEESEPLTGFLTTVPDGPTQLRALNLTEGFAVLHW



KPPQNPVDTYDVQVTAPGAPPLQAETPGSAVDYPLHDLVLHTNYTATVRGLRGPNLTSPASITFTTGLEAPRDLEAKEVTPRTAL



LTWTEPPVRPAGYLLSFHTPGGQNQEILLPGGITSHQLLGLFPSTSYNARLQAMWGQSLLPPVSTSFTTGGLRIPFPRDCGEEMQ



NGAGASRTSTIFLNGNRERPLNVFCDMETDGGGWLVFQRRMDGQTDFWRDWEDYAHGFGNISGEFWLGNEALHSLTQAGDYSMRV



DLRAGDEAVFAQYDSFHVDSAAEYYRLHLEGYHGTAGDSMSYHSGSVFSARDRDPNSLLISCAVSYRGAWWYRNCHYANLNGLYG



STVDHQGVSWYHWKGFEFSVPFTEMKLRPRNFRSPAGGG




AIFM1_
MFRCGGLAAGALKQKLVPLVRTVCVRSPRQRNRLPGNLFQRWHVPLELQMTRQMASSGASGGKIDNSVLVLIVGLSTVGAGAYAY
9

HUMAN
KTMKEDEKRYNERISGLGLTPEQKQKKAALSASEGEEVPQDKAPSHVPFLLIGGGTAAFAAARSIRARDPGARVLIVSEDPELPY



MRPPLSKELWFSDDPNVTKTLRFKQWNGKERSIYFQPPSFYVSAQDLPHIENGGVAVLTGKKVVQLDVRDNMVKLNDGSQITYEK



CLIATGGTPRSLSAIDRAGAEVKSRTTLFRKIGDFRSLEKISREVKSITIIGGGFLGSELACALGRKARALGTEVIQLFPEKGNM



GKILPEYLSNWTMEKVRREGVKVMPNAIVQSVGVSSGKLLIKLKDGRKVETDHIVAAVGLEPNVELAKTGGLEIDSDFGGFRVNA



ELQARSNIWVAGDAACFYDIKLGRRRVEHHDHAVVSGRLAGENMTGAAKPYWHQSMFWSDLGPDVGYEAIGLVDSSLPTVGVFAK



ATAQDNPKSATEQSGTGIRSESETESEASEITIPPSTPAVPQAPVQGEDYGKGVIFYLRDKVVVGIVLWNIFNRMPIARKIIKDG



EQHEDLNEVAKLFNIHED




IBP3_
MQRARPTLWAAALTLLVLLRGPPVARAGASSAGLGPVVRCEPCDARALAQCAPPPAVCAELVREPGCGCCLTCALSEGQPCGIYT
10

HUMAN
ERCGSGLRCQPSPDEARPLQALLDGRGLCVNASAVSRLRAYLLPAPPAPGNASESEEDRSAGSVESPSVSSTHRVSDPKFHPLHS



KIIIIKKGHAKDSQRYKVDYESQSTDTQNFSSESKRETEYGPCRREMEDTLNHLKFLNVLSPRGVHIPNCDKKGFYKKKQCRPSK



GRKRGFCWCVDKYGQPLPGYTTKGKEDVHCYSMQSK




GELS_
MAPHRPAPALLCALSLALCALSLPVRAATASRGASQAGAPQGRVPEARPNSMVVEHPEFLKAGKEPGLQIWRVEKFDLVPVPTNL
11

HUMAN
YGDFFTGDAYVILKTVQLRNGNLQYDLHYWLGNECSQDESGAAAIFTVQLDDYLNGRAVQHREVQGFESATFLGYFKSGLKYKKG



GVASGFKHVVPNEVVVQRLFQVKGRRVVRATEVPVSWESFNNGDCFILDLGNNIHQWCGSNSNRYERLKATQVSKGIRDNERSGR



ARVHVSEEGTEPEAMLQVLGPKPALPAGTEDTAKEDAANRKLAKLYKVSNGAGTMSVSLVADENPFAQGALKSEDCFILDHGKDG



KIFVWKGKQANTEERKAALKTASDFITKMDYPKQTQVSVLPEGGETPLFKQFFKNWRDPDQTDGLGLSYLSSHIANVERVPFDAA



TLHTSTAMAAQHGMDDDGTGQKQIWRIEGSNKVPVDPATYGQFYGGDSYIILYNYRHGGRQGQIIYNWQGAQSTQDEVAASAILT



AQLDEELGGTPVQSRVVQGKEPAHLMSLFGGKPMIIYKGGTSREGGQTAPASTRLFQVRANSAGATRAVEVLPKAGALNSNDAFV



LKTPSAAYLWVGTGASEAEKTGAQELLRVLRAQPVQVAEGSEPDGFWEALGGKAAYRTSPRLKDKKMDAHPPRLFACSNKIGRFV



IEEVPGELMQEDLATDDVMLLDTWDQVFVWVGKDSQEEEKTEALTSAKRYIETDPANRDRRTPITVVKQGFEPPSFVGWFLGWDD



DYWSVDPLDRAMAELAAGCGCGCCCGCTGCAGGCGCTGCTGGATGGCCGCGGCCTGTGCGTGAACGCGAGCGCGGTGAGCCGCCT



GCGCGCGTATCTGCTGCCGGCGCCGCCGGCGCCGGGCGAACCGCCGGCGCCGGGCAACGCGAGCGAAAGCGAAGAAGATCGCAGC



GCGGGCAGCGTGGAAAGCCCGAGCGTGAGCAGCACCCATCGCGTGAGCGATCCGAAATTTCATCCGCTGCATAGCAAAATTATTA



TTATTAAAAAAGGCCATGCGAAAGATAGCCAGCGCTATAAAGTGGATTATGAAAGCCAGAGCACCGATACCCAGAACTTTAGCAG



CGAAAGCAAACGCGAAACCGAATATGGCCCGTGCCGCCGCGAAATGGAAGATACCCTGAACCATCTGAAATTTCTGAACGTGCTG



AGCCCGCGCGGCGTGCATATTCCGAACTGCGATAAAAAAGGCTTTTATAAAAAAAAACAGTGCCGCCCGAGCAAAGGCCGCAAAC



GCGGCTTTTGCTGGTGCGTGGATAAATATGGCCAGCCGCTGCCGGGCTATACCACCAAAGGCAAAGAAGATGTGCATTGCTATAG



CATGCAGAGCAAA




MASP1_
MRWLLLYYALCFSLSKASAHTVELNNMFGQIQSPGYPDSYPSDSEVTWNITVPDGFRIKLYFMHFNLESSYLCEYDYVKVETEDQ
12

HUMAN
VLATFCGRETTDTEQTPGQEVVLSPGSFMSITFRSDFSNEERFTGFDAHYMAVDVDECKEREDEELSCDHYCHNYIGGYYCSCRF



GYILHTDNRTCRVECSDNLFTQRTGVITSPDFPNPYPKSSECLYTIELEEGFMVNLQFEDIFDIEDHPEVPCPYDYIKIKVGPKV



LGPFCGEKAPEPISTQSHSVLILFHSDNSGENRGWRLSYRAAGNECPELQPPVHGKIEPSQAKYFFKDQVLVSCDTGYKVLKDNV



EMDTFQIECLKDGTWSNKIPTCKIVDCRAPGELEHGLITFSTRNNLTTYKSEIKYSCQEPYYKMLNNNTGIYTCSAQGVWMNKVL



GRSLPTCLPVCGLPKFSRKLMARIFNGRPAQKGTTPWIAMLSHLNGQPFCGGSLLGSSWIVTAAHCLHQSLDPEDPTLRDSDLLS



PSDFKIILGKHWRLRSDENEQHLGVKHTTLHPQYDPNTFENDVALVELLESPVLNAFVMPICLPEGPQQEGAMVIVSGWGKQFLQ



RFPETLMEIEIPIVDHSTCQKAYAPLKKKVTRDMICAGEKEGGKDACAGDSGGPMVTLNRERGQWYLVGTVSWGDDCGKKDRYGV



YSYIHHNKDWIQRVTGVRN




COIA1_
MAPYPCGCHILLLLFCCLAAARANLLNLNWLWFNNEDTSHAATTIPEPQGPLPVQPTADTTTHVTPRNGSTEPATAPGSPEPPSE
13

HUMAN
LLEDGQDTPTSAESPDAPEENIAGVGAEILNVAKGIRSFVQLWNDTVPTESLARAETLVLETPVGPLALAGPSSTPQENGTTLWP



SRGIPSSPGAHTTEAGTLPAPTPSPPSLGRPWAPLTGPSVPPPSSGRASLSSLLGGAPPWGSLQDPDSQGLSPAAAAPSQQLQRP



DVRLRTPLLHPLVMGSLGKHAAPSAFSSGLPGALSQVAVTTLTRDSGAWVSHVANSVGPGLANNSALLGADPEAPAGRCLPLPPS



LPVCGHLGISRFWLPNHLHHESGEQVRAGARAWGGLLQTHCHPFLAWFFCLLLVPPCGSVPPPAPPPCCQFCEALQDACWSRLGG



GRLPVACASLPTQEDGYCVLIGPAAERISEEVGLLQLLGDPPPQQVTQTDDPDVGLAYVFGPDANSGQVARYHFPSLFFRDFSLL



FHIRPATEGPGVLFAITDSAQAMVLLGVKLSGVQDGHQDISLLYTEPGAGQTHTAASFRLPAFVGQWTHLALSVAGGFVALYVDC



EEFQRMPLARSSRGLELEPGAGLFVAQAGGADPDKFQGVIAELKVRRDPQVSPMHCLDEEGDDSDGASGDSGSGLGDARELLREE



TGAALKPRLPAPPPVTTPPLAGGSSTEDSRSEEVEEQTTVASLGAQTLPGSDSVSTWDGSVRTPGGRVKEGGLKGQKGEPGVPGP



PGRAGPPGSPCLPGPPGLPCPVSPLGPAGPALQTVPGPQGPPGPPGRDGTPGRDGEPGDPGEDGKPGDTGPQGFPGTPGDVGPKG



DKGDPGVGERGPPGPQGPPGPPGPSFRHDKLTFIDMEGSGFGGDLEALRGPRGFPGPPGPPGVPGLPGEPGRFGVNSSDVPGPAG



LPGVPGREGPPGFPGLPGPPGPPGREGPPGRTGQKGSLGEAGAPGHKGSKGAPGPAGARGESGLAGAPGPAGPPGPPGPPGPPGP



GLPAGFDDMEGSGGPFWSTARSADGPQGPPGLPGLKGDPGVPGLPGAKGEVGADGVPGFPGLPGREGIAGPQGPKGDRGSRGEKG



DPGKDGVGQPGLPGPPGPPGPVVYVSEQDGSVLSVPGPEGRPGFAGFPGPAGPKGNLGSKGERGSPGPKGEKGEPGSIFSPDGGA



LGPAQKGAKGEPGFRGPPGPYGRPGYKGEIGFPGRPGRPGMNGLKGEKGEPGDASLGFGMRGMPGPPGPPGPPGPPGTPVYDSNV



FAESSRPGPPGLPGNQGPPGPKGAKGEVGPPGPPGQFPFDFLQLEAEMKGEKGDRGDAGQKGERGEPGGGGFFGSSLPGPPGPPG



PPGPRGYPGIPGPKGESIRGQPGPPGPQGPPGIGYEGRQGPPGPPGPPGPPSFPGPHRQTISVPGPPGPPGPPGPPGTMGASSGV



RLWATRQAMLGQVHEVPEGWLIFVAEQEELYVRVQNGFRKVQLEARTPLPRGTDNEVAALQPPVVQLHDSNPYPRREHPHPTARP



WRADDILASPPRLPEPQPYPGAPHHSSYVHLRPARPTSPPAHSHRDFQPVLHLVALNSPLSGGMRGIRGADFQCFQQARAVGLAG



TFRAFLSSRLQDLYSIVRRADRAAVPIVNLKDELLFPSWEALFSGSEGPLKPGARIFSFDGKDVLRHPTWPQKSVWHGSDPNGRR



LTESYCETWRTEAPSATGQASSLLGGRLLGQSAASCHHAYIVLCIENSFMTASK




GRP78_
MKLSLVAAMLLLLSAARAEEEDKKEDVGTVVGIDLGTTYSCVGVFKNGRVEIIANDQGNRITPSYVAFTPEGERLIGDAAKNQLT
14

HUMAN
SNPENTVFDAKRLIGRTWNDPSVQQDIKFLPFKVVEKKTKPYIQVDIGGGQTKTFAPEEISAMVLTKMKETAEAYLGKKVTHAVV



TVPAYFNDAQRQATKDAGTIAGLNVMRIINEPTAAAIAYGLDKREGEKNILVFDLGGGTFDVSLLTIDNGVFEVVATNGDTHLGG



EDFDQRVMEHFIKLYKKKTGKDVRKDNRAVQKLRREVEKAKRALSSQHQARIEIESFYEGEDFSETLTRAKFEELNMDLFRSTMK



PVQKVLEDSDLKKSDIDEIVLVGGSTRIPKIQQLVKEFFNGKEPSRGINPDEAVAYGAAVQAGVLSGDQDTGDLVLLDVCPLTLG



IETVGGVMTKLIPRNTVVPTKKSQIFSTASDNQPTVTIKVYEGERPLTKDNHLLGTFDLTGIPPAPRGVPQIEVTFEIDVNGILR



VTAEDKGTGNKNKITITNDQNRLTPEEIERMVNDAEKFAEEDKKLKERIDTRNELESYAYSLKNQIGDKEKLGGKLSSEDKETME



KAVEEKIEWLESHQDADIEDFKAKKKELEEIVQPIISKLYGSAGPPPTGEEDTAEKDEL




KIT_
MRGARGAWDFLCVLLLLLRVQTGSSQPSVSPGEPSPPSIHPGKSDLIVRVGDEIRLLCTDPGFVKWTFEILDETNENKQNEWITE
15

HUMAN
KAEATNTGKYTCTNKHGLSNSIYVFVRDPAKLFLVDRSLYGKEDNDTLVRCPLTDPEVTNYSLKGCQGKPLPKDLRFIPDPKAGI



MIKSVKRAYHRLCLHCSVDQEGKSVLSEKFILKVRPAFKAVPVVSVSKASYLLREGEEFTVTCTIKDVSSSVYSTWKRENSQTKL



QEKYNSWHHGDFNYERQATLTISSARVNDSGVFMCYANNTFGSANVTTTLEVVDKGFINIFPMINTTVFVNDGENVDLIVEYEAF



PKPEHQQWIYMNRTFTDKWEDYPKSENESNIRYVSELHLTRLKGTEGGTYTFLVSNSDVNAAIAFNVYVNTKPEILTYDRLVNGM



LQCVAAGFPEPTIDWYFCPGTEQRCSASVLPVDVQTLNSSGPPFGKLVVQSSIDSSAFKHNGTVECKAYNDVGKTSAYFNFAFKG



NNKEQIHPHTLFTPLLIGFVIVAGMMCIIVMILTYKYLQKPMYEVQWKVVEEINGNNYVYIDPTQLPYDHKWEFPRNRLSFGKTL



GAGAFGKVVEATAYGLIKSDAAMTVAVKMLKPSAHLTEREALMSELKVLSYLGNHMNIVNLLGACTIGGPTLVITEYCCYGDLLN



FLRRKRDSFICSKQEDHAEAALYKNLLHSKESSCSDSTNEYMDMKPGVSYVVPTKADKRRSVRIGSYIERDVTPAIMEDDELALD



LEDLLSFSYQVAKGMAFLASKNCIHRDLAARNILLTHGRITKICDFGLARDIKNDSNYVVKGNARLPVKWMAPESIFNCVYTFES



DVWSYGIFLWELFSLGSSPYPGMPVDSKFYKMIKEGFRMLSPEHAPAEMYDIMKTCWDADPLKRPTFKQIVQLIEKQISESTNHI



YSNLANCSPNRQKPVVDHSVRINSVGSTASSSQPLLVHDDV




LRP1_
MLTPPLLLLLPLLSALVAAAIDAPKTCSPKQFACRDQITCISKGWRCDGERDCPDGSDEAPEICPQSKAQRCQPNEHNCLGTELC
16

HUMAN
VPMSRLCNGVQDCMDGSDEGPHCRELQGNCSRLGCQHHCVPTLDGPTCYCNSSFQLQADGKTCKDFDECSVYGTCSQLCTNTDGS



FICGCVEGYLLQPDNRSCKAKNEPVDRPPVLLIANSQNILATYLSGAQVSTITPTSTRQTTAMDFSYANETVCWVHVGDSAAQTQ



LKCARMPGLKGFVDEHTINISLSLHHVEQMAIDWLTGNFYFVDDIDDRIFVCNRNGDTCVTLLDLELYNPKGIALDPAMGKVFFT



DYGQIPKVERCDMDGQNRTKLVDSKIVFPHGITLDLVSRLVYWADAYLDYIEVVDYEGKGRQTIIQGILIEHLYGLTVFENYLYA



TNSDNANAQQKTSVIRVNRFNSTEYQVVTRVDKGGALHIYHQRRQPRVRSHACENDQYGKPGGCSDICLLANSHKARTCRCRSGF



SLGSDGKSCKKPEHELFLVYGKGRPGIIRGMDMGAKVPDEHMIPIENLMNPRALDFHAETGFIYFADTTSYLIGRQKIDGTERET



ILKDGIHNVEGVAVDWMGDNLYWTDDGPKKTISVARLEKAAQTRKTLIEGKMTHPRAIVVDPLNGWMYWTDWEEDPKDSRRGRLE



RAWMDGSHRDIFVTSKTVLWPNGLSLDIPAGRLYWVDAFYDRIETILLNGTDRKIVYEGPELNHAFGLCHHGNYLFWTEYRSGSV



YRLERGVGGAPPTVTLLRSERPPIFEIRMYDAQQQQVGTNKCRVNNGGCSSLCLATPGSRQCACAEDQVLDADGVTCLANPSYVP



PPQCQPGEFACANSRCIQERWKCDGDNDCLDNSDEAPALCHQHTCPSDRFKCENNRCIPNRWLCDGDNDCGNSEDESNATCSART



CPPNQFSCASGRCIPISWTCDLDDDCGDRSDESASCAYPTCFPLTQFTCNNGRCININWRCDNDNDCGDNSDEAGCSHSCSSTQF



KCNSGRCIPEHWTCDGDNDCGDYSDETHANCTNQATRPPGGCHTDEFQCRLDGLCIPLRWRCDGDTDCMDSSDEKSCEGVTHVCD



PSVKFGCKDSARCISKAWVCDGDNDCEDNSDEENCESLACRPPSHPCANNTSVCLPPDKLCDGNDDCGDGSDEGELCDQCSLNNG



GCSHNCSVAPGEGIVCSCPLGMELGPDNHTCQIQSYCAKHLKCSQKCDQNKFSVKCSCYEGWVLEPDGESCRSLDPFKPFIIFSN



RHEIRRIDLHKGDYSVLVPGLRNTIALDFHLSQSALYWTDVVEDKIYRGKLLDNGALTSFEVVIQYGLATPEGLAVDWIAGNIYW



VESNLDQIEVAKLDGTLRTTLLAGDIEHPRAIALDPRDGILFWTDWDASLPRIEAASMSGAGRRTVHRETGSGGWPNGLTVDYLE



KRILWIDARSDAIYSARYDGSGHMEVLRGHEFLSHPFAVTLYGGEVYWTDWRTNTLAKANKWTGHNVTVVQRTNTQPFDLQVYHP



SRQPMAPNPCEANGGQGPCSHLCLINYNRTVSCACPHLMKLHKDNTTCYEFKKFLLYARQMEIRGVDLDAPYYNYIISFTVPDID



NVTVLDYDAREQRVYWSDVRTQAIKRAFINGTGVETVVSADLPNAHGLAVDWVSRNLFWTSYDTNKKQINVARLDGSFKNAVVQG



LEQPHGLVVHPLRGKLYWTDGDNISMANMDGSNRTLLFSG



QKGPVGLAIDFPESKLYWISSGNHTINRCNLDGSGLEVIDAMRSQLGKATALAIMGDKLWWADQVSEKMGTCSKADGSGSVVLRN



STTLVMHMKVYDESIQLDHKGTNPCSVNNGDCSQLCLPTSETTRSCMCTAGYSLRSGQQACEGVGSFLLYSVHEGIRGIPLDPND



KSDALVPVSGTSLAVGIDFHAENDTIYWVDMGLSTISRAKRDQTWREDVVTNGIGRVEGIAVDWIAGNIYWTDQGFDVIEVARLN



GSFRYVVISQGLDKPRAITVHPEKGYLFWTEWGQYPRIERSRLDGTERVVLVNVSISWPNGISVDYQDGKLYWCDARTDKIERID



LETGENREVVLSSNNMDMFSVSVFEDFIYWSDRTHANGSIKRGSKDNATDSVPLRTGIGVQLKDIKVFNRDRQKGTNVCAVANGG



CQQLCLYRGRGQRACACAHGMLAEDGASCREYAGYLLYSERTILKSIHLSDERNLNAPVQPFEDPEHMKNVIALAFDYRAGTSPG



TPNRIFFSDIHFGNIQQINDDGSRRITIVENVGSVEGLAYHRGWDTLYWTSYTTSTITRHTVDQTRPGAFERETVITMSGDDHPR



AFVLDECQNLMFWTNWNEQHPSIMRAALSGANVLTLIEKDIRTPNGLAIDHRAEKLYFSDATLDKIERCEYDGSHRYVILKSEPV



HPFGLAVYGEHIFWTDWVRRAVQRANKHVGSNMKLLRVDIPQQPMGIIAVANDTNSCELSPCRINNGGCQDLCLLTHQGHVNCSC



RGGRILQDDLTCRAV



NSSCRAQDEFECANGECINFSLTCDGVPHCKDKSDEKPSYCNSRRCKKTFRQCSNGRCVSNMLWCNGADDCGDGSDEIPCNKTAC



GVGEFRCRDGTCIGNSSRCNQFVDCEDASDEMNCSATDCSSYFRLGVKGVLFQPCERTSLCYAPSWVCDGANDCGDYSDERDCPG



VKRPRCPLNYFACPSGRCIPMSWTCDKEDDCEHGEDETHCNKFCSEAQFECQNHRCISKQWLCDGSDDCGDGSDEAAHCEGKTCG



PSSFSCPGTHVCVPERWLCDGDKDCADGADESIAAGCLYNSTCDDREFMCQNRQCIPKHFVCDHDRDCADGSDESPECEYPTCGP



SEFRCANGRCLSSRQWECDGENDCHDQSDEAPKNPHCTSQEHKCNASSQFLCSSGRCVAEALLCNGQDDCGDSSDERGCHINECL



SRKLSGCSQDCEDLKIGFKCRCRPGFRLKDDGRTCADVDECSTTFPCSQRCINTHGSYKCLCVEGYAPRGGDPHSCKAVTDEEPF



LIFANRYYLRKLNLDGSNYTLLKQGLNNAVALDFDYREQMIYWTDVTTQGSMIRRMHLNGSNVQVLHRTGLSNPDGLAVDWVGGN



LYWCDKGRDTIEVSKLNGAYRTVLVSSGLREPRALVVDVQNGYLYWTDWGDHSLIGRIGMDGSSRSVIVDTKITWPNGLTLDYVT



ERIYWADAREDYIEFASLDGSNRHVVLSQDIPHIFALTLFEDYVYWTDWETKSINRAHKTTGTNKTLLISTLHRPMDLHVFHALR



QPDVPNHPCKVNNGGCSNLCLLSPGGGHKCACPTNFYLGSDGRTCVSNCTASQFVCKNDKCIPFWWKCDTEDDCGDHSDEPPDCP



EFKCRPGQFQCSTGICTNPAFICDGDNDCQDNSDEANCDIHVCLPSQFKCTNTNRCIPGIFRCNGQDNCGDGEDERDCPEVTCAP



NQFQCSITKRCIPRVWVCDRDNDCVDGSDEPANCTQMTCGVDEFRCKDSGRCIPARWKCDGEDDCGDGSDEPKEECDERTCEPYQ



FRCKNNRCVPGRWQCDYDNDCGDNSDEESCTPRPCSESEFSCANGRCIAGRWKCDGDHDCADGSDEKDCTPRCDMDQFQCKSGHC



IPLRWRCDADADCMDGSDEEACGTGVRTCPLDEFQCNNTLCKPLAWKCDGEDDCGDNSDENPEECARFVCPPNRPFRCKNDRVCL



WIGRQCDGTDNCGDGTDEEDCEPPTAHTTHCKDKKEFLCRNQRCLSSSLRCNMFDDCGDGSDEEDCSIDPKLTSCATNASICGDE



ARCVRTEKAAYCACRSGFHTVPGQPGCQDINECLRFGTCSQLCNNTKGGHLCSCARNFMKTHNTCKAEGSEYQVLYIADDNEIRS



LFPGHPHSAYEQAFQGDESVRIDAMDVHVKAGRVYWTNWHTGTISYRSLPPAAPPTTSNRHRRQIDRGVTHLNISGLKMPRGIAI



DWVAGNVYWTDSGRDVIEVAQMKGENRKTLISGMIDEPHAIVVDPLRGTMYWSDWGNHPKIETAAMDGTLRETLVQDNIQWPTGL



AVDYHNERLYWADAKLSVIGSIRLNGTDPIVAADSKRGLSHPFSIDVFEDYIYGVTYINNRVFKIHKFGHSPLVNLTGGLSHASD



VVLYHQHKQPEVTNPCDRKKCEWLCLLSPSGPVCTCPNGKRLDNGTCVPVPSPTPPPDAPRPGTCNLQCFNGGSCFLNARRQPKC



RCQPRYTGDKCELDQCWEHCRNGGTCAASPSGMPTCRCPT



GFTGPKCTQQVCAGYCANNSTCTVNQGNQPQCRCLPGFLGDRCQYRQCSGYCENFGTCQMAADGSRQCRCTAYFEGSRCEVNKCS



RCLEGACVVNKQSGDVTCNCTDGRVAPSCLTCVGHCSNGGSCTMNSKMMPECQCPPHMTGPRCEEHVFSQQQPGHIASILIPLLL



LLLLVLVAGVVFWYKRRVQGAKGFQHQRMTNGAMNVEIGNPTYKMYEGGEPDDVGGLLDADFALDPDKPTNFTNPVYATLYMGGH



GSRHSLASTDEKRELLGRGPEDEIGDPLA




PEDF_
MQALVLLLCIGALLGHSSCQNPASPPEEGSPDPDSTGALVEEEDPFFKVPVNKLAAAVSNFGYDLYRVRSSTSPTTNVLLSPLSV
17

HUMAN
ATALSALSLGAEQRTESIIHRALYYDLISSPDIHGTYKELLDTVTAPQKNLKSASRIVFEKKLRIKSSFVAPLEKSYGTRPRVLT



GNPRLDLQEINNWVQAQMKGKLARSTKEIPDEISILLLGVAHFKGQWVTKFDSRKTSLEDFYLDEERTVRVPMMSDPKAVLRYGL



DSDLSCKIAQLPLTGSMSIIFFLPLKVTQNLTLIEESLTSEFIHDIDRELKTVQAVLTVPKLKLSYEGEVTKSLQEMKLQSLFDS



PDFSKITGKPIKLTQVEHRAGFEWNEDGAGTTPSPGLQPAHLTFPLDYHLNQPFIFVLRDTDTGALLFIGKILDPRGP




LUM_
MSLSAFTLFLALIGGTSGQYYDYDFPLSIYGQSSPNCAPECNCPESYPSAMYCDELKLKSVPMVPPGIKYLYLRNNQIDHIDEKA
18

HUMAN
FENVTDLQWLILDHNLLENSKIKGRVFSKLKQLKKLHINHNNLTESVGPLPKSLEDLQLTHNKITKLGSFEGLVNLTFIHLQHNR



LKEDAVSAAFKGLKSLEYLDLSFNQIARLPSGLPVSLLTLYLDNNKISNIPDEYFKRFNALQYLRLSHNELADSGIPGNSFNVSS



LVELDLSYNKLKNIPTVNENLENYYLEVNQLEKFDIKSFCKILGPLSYSKIKHLRLDGNRISETSLPPDMYECLRVANEVTLN




C163A_
MSKLRMVLLEDSGSADFRRHFVNLSPFTITVVLLLSACFVTSSLGGTDKELRLVDGENKCSGRVEVKVQEEWGTVCNNGWSMEAV
19

HUMAN
SVICNQLGCPTAIKAPGWANSSAGSGRIWMDHVSCRGNESALWDCKHDGWGKHSNCTHQQDAGVTCSDGSNLEMRLTRGGNMCSG



RIEIKFQGRWGTVCDDNFNIDHASVICRQLECGSAVSFSGSSNFGEGSGPIWFDDLICNGNESALWNCKHQGWGKHNCDHAEDAG



VICSKGADLSLRLVDGVTECSGRLEVRFQGEWGTICDDGWDSYDAAVACKQLGCPTAVTAIGRVNASKGFGHIWLDSVSCQGHEP



AIWQCKHHEWGKHYCNHNEDAGVTCSDGSDLELRLRGGGSRCAGTVEVEIQRLLGKVCDRGWGLKEADVVCRQLGCGSALKTSYQ



VYSKIQATNTWLFLSSCNGNETSLWDCKNWQWGGLTCDHYEEAKITCSAHREPRLVGGDIPCSGRVEVKHGDTWGSICDSDFSLE



AASVLCRELQCGTVVSILGGAHFGEGNGQIWAEEFQCEGHESHLSLCPVAPRPEGTCSHSRDVGVVCSRYTEIRLVNGKTPCEGR



VELKTLGAWGSLCNSHWDIEDAHVLCQQLKCGVALSTPGGARFGKGNGQIWRHMFHCTGTEQHMGDCPVTALGASLCPSEQVASV



ICSGNQSQTLSSCNSSSLGPTRPTIPEESAVACIESGQLRLVNGGGRCAGRVEIYHEGSWGTICDDSWDLSDAHVVCRQLGCGEA



INATGSAHFGEGTGPIWLDEMKCNGKESRIWQCHSHGWGQQNCRHKEDAGVICSEFMSLRLTSEASREACAGRLEVFYNGAWGTV



GKSSMSETTVGVVCRQLGCADKGKINPASLDKAMSIPMWVDNVQCPKGPDTLWQCPSSPWEKRLASPSEETWITCDNKIRLQEGP



TSCSGRVEIWHGGSWGTVCDDSWDLDDAQVVCQQLGCGPALKAFKEAEFGQGTGPIWLNEVKCKGNESSLWDCPARRWGHSECGH



KEDAAVNCTDISVQKTPQKATTGRSSRQSSFIAVGILGVVLLAIFVALFFLTKKRRQRQRLAVSSRGENLVHQIQYREMNSCLNA



DDLDLMNSSENSHESADFSAAELISVSKFLPISGMEKEAILSHTEKENGNL




PTPRJ_
MKPAAREARLPPRSPGLRWALPLLLLLLRLGQILCAGGTPSPIPDPSVATVATGENGITQISSTAESFHKQNGTGTPQVETNTSE
20

HUMAN
DGESSGANDSLRTPEQGSNGTDGASQKTPSSTGPSPVFDIKAVSISPTNVILTWKSNDTAASEYKYVVKHKMENEKTITVVHQPW



CNITGLRPATSYVFSITPGIGNETWGDPRVIKVITEPIPVSDLRVALTGVRKAALSWSNGNGTASCRVLLESIGSHEELTQDSRL



QVNISGLKPGVQYNINPYLLQSNKTKGDPLGTEGGLDASNTERSRAGSPTAPVHDESLVGPVDPSSGQQSRDTEVLLVGLEPGTR



YNATVYSQAANGTEGQPQAIEFRTNAIQVFDVTAVNISATSLTLIWKVSDNESSSNYTYKIHVAGETDSSNLNVSEPRAVIPGLR



SSTFYNITVCPVLGDIEGTPGFLQVHTPPVPVSDFRVTVVSTTEIGLAWSSHDAESFQMHITQEGAGNSRVEITTNQSIIIGGLF



PGTKYCFEIVPKGPNGTEGASRTVCNRTVPSAVFDIHVVYVTTTEMWLDWKSPDGASEYVYHLVIESKHGSNHTSTYDKAITLQG



LIPGTLYNITISPEVDHVWGDPNSTAQYTRPSNVSNIDVSTNTTAATLSWQNFDDASPTYSYCLLIEKAGNSSNATQVVTDIGIT



DATVTELIPGSSYTVEIFAQVGDGIKSLEPGRKSFCTDPASMASFDCEVVPKEPALVLKWTCPPGANAGFELEVSSGAWNNATHL



ESCSSENGTEYRTEVTYLNFSTSYNISITTVSCGKMAAPTRNTCTTGITDPPPPDGSPNITSVSHNSVKVKFSGFEASHGPIKAY



AVILTTGEAGHPSADVLKYTYEDFKKGASDTYVTYLIRTEEKGRSQSLSEVLKYEIDVGNESTTLGYYNGKLEPLGSYRACVAGF



TNITFHPQNKGLIDGAESYVSFSRYSDAVSLPQDPGVICGAVFGCIFGALVIVTVGGFIFWRKKRKDAKNNEVSFSQIKPKKSKL



IRVENFEAYFKKQQADSNCGFAEEYEDLKLVGISQPKYAAELAENRGKNRYNNVLPYDISRVKLSVQTHSTDDYINANYMPGYHS



KKDFIATQGPLPNTLKDFWRMVWEKNVYAIIMLTKCVEQGRTKCEEYWPSKQAQDYGDITVAMTSEIVLPENTIRDFTVKNIQTS



ESHPLRQFHFTSWPDHGVPDTTDLLINFRYLVRDYMKQSPPESPILVHCSAGVGRTGTFIAIDRLIYQIENENTVDVYGIVYDLR



MHRPLMVQTEDQYVFLNQCVLDIVRSQKDSKVDLIYQNTTAMTIYENLAPVTTFGKTNGYIA




ALDOA_
MPYQYPALTPEQKKELSDIAHRIVAPGKGILAADESTGSIAKRLQSIGTENTEENRRFYRQLLLTADDRVNPCIGGVILFHETLY
21

HUMAN
QKADDGRPFPQVIKSKGGVVGIKVDKGVVPLAGTNGETTTQGLDGLSERCAQYKKDGADFAKWRCVLKIGEHTPSALAIMENANV



LARYASICQQNGIVPIVEPEILPDGDHDLKRCQYVTEKVLAAVYKALSDHHIYLEGTLLKPNMVTPGHACTQKFSHEEIAMATVT



ALRRTVPPAVTGITFLSGGQSEEEASINLNAINKCPLLKPWALTFSYGRALQASALKAWGGKKENLKAAQEEYVKRALANSLACQ



GKYTPSGQAGAAASESLFVSNHAY




FRIL_
MSSQIRQNYSTDVEAAVNSLVNLYLQASYTYLSLGFYFDRDDVALEGVSHFFRELAEEKREGYERLLKMQNQRGGRALFQDIKKP
22

HUMAN
AEDEWGKTPDAMKAAMALEKKLNQALLDLHALGSARTDPHLCDFLETHFLDEEVKLIKKMGDHLTNLHRLGGPEAGLGEYLFERL



TLKHD




TETN_
MELWGAYLLLCLFSLLTQVTTEPPTQKPKKIVNAKKDVVNTKMFEELKSRLDTLAQEVALLKEQQALQTVCLKGTKVHMKCFLAF
23

HUMAN
TQTKTFHEASEDCISRGGTLGTPQTGSENDALYEYLRQSVGNEAEIWLGLNDMAAEGTWVDMTGARIAYKNWETEITAQPDGGKT



ENCAVLSGAANGKWFDKRCRDQLPYICQFGIV




ISLR_
MQELHLLWWALLLGLAQACPEPCDCGEKYGFQIADCAYRDLESVPPGFPANVTTLSLSANRLPGLPEGAFREVPLLQSLWLAHNE
24

HUMAN
IRTVAAGALASLSHLKSLDLSHNLISDFAWSDLHNLSALQLLKMDSNELTFIPRDAFRSLRALRSLQLNHNRLHTLAEGTFTPLT



ALSHLQINENPFDCTCGIVWLKTWALTTAVSIPEQDNIACTSPHVLKGTPLSRLPPLPCSAPSVQLSYQPSQDGAELRPGFVLAL



HCDVDGQPAPQLHWHIQIPSGIVEITSPNVGTDGRALPGTPVASSQPRFQAFANGSLLIPDFGKLEEGTYSCLATNELGSAESSV



DVALATPGEGGEDTLGRRFHGKAVEGKGCYTVDNEVQPSGPEDNVVIIYLSRAGNPEAAVAEGVPGQLPPGLLLLGQSLLLFFFL



TSF