- Open Access
Identification of a seven glycopeptide signature for malignant pleural mesothelioma in human serum by selected reaction monitoring
Clinical Proteomics volume 10, Article number: 16 (2013)
Serum biomarkers can improve diagnosis and treatment of malignant pleural mesothelioma (MPM). However, the evaluation of potential new serum biomarker candidates is hampered by a lack of assay technologies for their clinical evaluation. Here we followed a hypothesis-driven targeted proteomics strategy for the identification and clinical evaluation of MPM candidate biomarkers in serum of patient cohorts.
Based on the hypothesis that cell surface exposed glycoproteins are prone to be released from tumor-cells to the circulatory system, we screened the surfaceome of model cell lines for potential MPM candidate biomarkers. Selected Reaction Monitoring (SRM) assay technology allowed for the direct evaluation of the newly identified candidates in serum. Our evaluation of 51 candidate biomarkers in the context of a training and an independent validation set revealed a reproducible glycopeptide signature of MPM in serum which complemented the MPM biomarker mesothelin.
Our study shows that SRM assay technology enables the direct clinical evaluation of protein-derived candidate biomarker panels for which clinically reliable ELISA’s currently do not exist.
Malignant pleural mesothelioma (MPM) is a fatal cancer of the pleura induced by asbestos exposure. Treatments developed over the last decade have improved patient survival [1–5]. However, their efficacy is limited by the frequent detection of MPM only at advanced stages [6, 7]. Easily and longitudinally accessible blood biomarkers are expected to support diagnosis and therapy selection at early disease stages, when benefit from treatment is the highest . To date, the best available MPM biomarker in serum is mesothelin . While the protein is frequently elevated at advanced stages of the disease, its value for early detection remains limited . The search and evaluation of additional MPM biomarkers in serum remains thus a priority. Generally, this is approached applying enzyme linked immunosorbent assays (ELISA), which commonly allow for the reliable evaluation of only one biomarker candidate at the time, like in the case of the recently proposed fibulin-3 protein . An alternative would be the investigation of panels of simultaneously measured biomarkers. Such a multiplexed strategy would be a more efficient approach in terms of samples consumption and diagnostic accuracy . To achieve this goal, in our study we developed and applied a hypothesis-driven targeted proteomics strategy which enabled the parallel quantitative evaluation of potential MPM candidate biomarkers in serum through SRM assay technology.
SRM assay technology relies on the ability of a triple-quadrupole mass spectrometer (QQQ) to selectively isolate predefined peptides of interest in a complex protein mixture after enzymatic digestion (usually using trypsin) . SRM-assays encompass the analytical coordinates necessary for the unambiguous detection and quantification of the target candidates . They consist of selected peptide-transitions to monitor (the pairs of signals representing the precursor peptide-ion and a corresponding fragment-ion), best collision energies to apply for peptide-fragmentations in the mass spectrometer and retention times of the target peptides in a chromatographic separation column. In a single SRM analysis, dozens of peptides are simultaneously quantified in complex samples with high sensitivity and reproducibility as surrogates of their proteins . This multiplexing potential enables parallel testing and clinical evaluation of proposed candidate biomarkers in clinically relevant specimens [16, 17].
In our study, to identify new candidate biomarkers for MPM we performed quantitative discovery-driven screening of the surfaceome of MPM model cell lines. This has previously been proposed as a valuable source of MPM candidate biomarkers . Subsequently, we used SRM assay technology to clinically evaluate the surfaceome-derived MPM candidate biomarkers in serum samples of suitably collected cohorts of MPM subjects and controls (Figure 1A and B).
Quantitative analysis of mesothelin in serum
To verify the SRM based approach for MPM biomarkers in serum, we performed quantitative investigations of mesothelin in samples enriched for N-linked glycopeptides (N-glycopeptides). This enrichment step was chosen in our study to efficiently and reproducibly reduce the analytical complexity of the serum proteome and to focus our investigation on a sub-proteome particularly relevant for biomarker research, like represented by the N-linked glycoproteins (N-glycoproteins) [19–23]. To generate and optimize SRM-assays for the quantitative analysis of mesothelin, we followed the strategy proposed by Picotti et al.  using spectral libraries from chemically synthesized peptide-sequences. We established SRM-assays for the N-glycopeptides of membrane-bound mesothelin  (Additional file 1: Table S1) which is shed in serum [26, 27] and for which a commercial ELISA-kit approved by the US Food and Drug Administration (FDA) (Mesomark®) [28, 29] is available. We selected quantitative SRM-assays based on the N-glycopeptide KWDVTSLETLK (UniProt/Swiss-Prot: Q13421; D replaces the formerly glycosylated N after deamidation with the enzyme PNGaseF in the protocol for N-glycopeptide enrichment in serum ) which was routinely observed in our experiments (Additional file 2: Figure S1). We then quantified mesothelin using SRM assay technology in enriched serum samples from 75 subjects (23 MPM, 26 healthy donors, HD, and 26 non-small cell lung cancer, NSCLC, Additional file 3: Table S2). In parallel, in the same subjects we performed mesothelin quantification based on the FDA approved ELISA (Additional file 3: Table S2). The results showed firstly, a significant correlation between the quantification of mesothelin assessed by SRM at the glycopeptide level and by ELISA at the protein level (Figure 2A). Furthermore, differences in mean concentrations of mesothelin measured by ELISA among the groups of MPM, HD and NSCLC were reflected by SRM assay technology (Figure 2B). These results provided evidence for the ability of the SRM based approach in serum to accurately identify and quantify MPM biomarkers of clinical relevance.
Surfaceome derived serum candidate biomarkers for malignant pleural mesothelioma
To identify candidate biomarkers for MPM, we performed a quantitative discovery-driven surfaceome screening in cell lines. To do so, we applied the mass spectrometry (MS) based Cell Surface Capture (CSC) technology  in two epithelioid and two biphasic MPM cell lines in parallel with two NSCLC (lung adenocarcinoma) and two non-cancerous pleural cell lines. A total of 668 N-glycopeptides were confidently (PeptideProphet probability ≥ 0.9) detected from more than 350 N-glycoproteins, which could potentially be shed into the blood stream. 514 N-glycopeptides were from MPM and 557 from non-MPM cell lines with 403 N-glycopeptides in common between the two groups. We prioritized candidate biomarkers potentially specific for MPM by focusing on N-glycopeptides reproducibly detected in higher abundance or in strong association with MPM cell lines. This screen led to the selection of 125 N-glycopeptides candidate biomarkers for MPM (Additional file 4: Table S3).
Subsequently, we applied SRM assay technology for the multiplexed assessment of these candidate biomarkers in serum. To generate SRM-assays for the candidate biomarkers, we used again spectral libraries from chemically synthesized peptide-sequences and acquired spectra without confounding background matrices. The strategy enabled us to avoid potential MS platform related transition-selection biases  and sufficient peptide amounts for establishing optimized SRM-assays . We obtained fragment spectra suitable for the library from 112 out of the initial 125 MPM candidate biomarkers and for these peptides we established initial SRM-assays in the background matrix of serum samples enriched for N-glycopeptides. We then performed multiplexed verification of the in serum unknown detectability of these candidate biomarkers by analyzing enriched serum samples of five MPM subjects (Additional file 5: Skyline file). In this initial screen we detected a total of 51 N-glycopeptides belonging to 36 N-glycoproteins which in the literature are reported at concentrations in serum/plasma between 180 μg/ml (hemopexin, UniProt/Swiss-Prot: P02790, peptide SWPAVGNCSSALR) and 0.92 ng/ml (poliovirus receptor-related protein 1, UniProt/Swiss-Prot: Q15223, peptide NPNGTVTVISR) (Additional file 6: Figure S2 and Additional file 7: Table S4) . The majority of these proteins are reported at concentrations below 100 ng/ml and thus in the range of MPM biomarkers previously proposed in the literature [11, 33, 34] (Figure 3). Taken together, these observations supported the suitability and sensitivity of our targeted SRM-based approach for MPM biomarker investigation in serum.
Seven glycopeptide signature for malignant pleural mesothelioma
To assess the diagnostic potential of the 51 MPM candidate biomarkers, we relatively quantified them via optimized SRM-assays in an initial set of enriched serum samples from clinical cohorts of MPM (N = 25), HD (N = 25) and NSCLC subjects (N = 25) (clinical characteristics are reported in Additional file 8: Table S5). We used these samples as the training set for a predictive analysis. SRM assay technology allowed us to quantify all 51 MPM candidate biomarkers in parallel in one single injection per sample (Additional file 9: Table S6). Statistical significance analysis  identified four N-glycopeptides significantly higher in abundance in MPM as compared to HD and six N-glycopeptides of higher abundance in MPM as compared to NSCLC (Additional file 7: Table S4). We selected these ten peptides as candidate predictors in two separate logistic regression models for binary responses between MPM and HD or MPM and NSCLC. We than performed a stepwise selection of predictors and evaluation of their predictive ability using receiver operating characteristic (ROC) curves to identify a multiplexed panel of MPM candidate biomarkers with best discriminatory performance for MPM. The panel was composed by the six glycopeptides derived from the N-glycoproteins intercellular adhesion molecule 1 (UniProt/Swiss-Prot: P05362, peptide sequence ADLTVVLLR), basement membrane-specific heparan sulfate proteoglycan core protein (P98160, ALVDFTR), anthrax toxin receptor 1 (Q9H6X2, DFDETQLAR), serum paraoxonase/arylesterase 1 (P27169, HADWTLTPLK), hypoxia up-regulated protein 1 (Q9Y4L1, VIDETWAWK) and thrombospondin-1 (P07996, VVDSTTGPGEHLR). In the training set, the panel discriminated MPM from HD with an area under the curve (AUC) of 0.94 (95% confidence interval, CI, [0.87, 0.99]) and the highest accuracy (highest true positives and true negatives) of 90% (95% CI, [82, 100]) at a cut-off of 0.61 (Figure 4A). Subsequently, we confirmed the performance of the six glycopeptide panel in an independent validation set of 87 sera newly enriched for N-glycopeptides from MPM (N = 30), HD (N = 29) and NSCLC (N = 28) subjects (Additional file 8: Table S5). Here the panel discriminated MPM from HD with AUC of 0.94 (95% CI, [0.86, 0.99]) and accuracy of 86% (95% CI, [78, 95]) at the 0.61 cut-off (Figure 4B). However, the panel failed to discriminate MPM from NSCLC: AUC in the training set was 0.77 (95% CI, [0.62, 0.91]) with accuracy of 74% (cut-off 0.35; 95% CI, [62, 88]) (Figure 4C) and in the validation set was 0.56 (95% CI, [0.40, 0.71]) with accuracy of 59% (cut-off 0.35; 95% CI, [50, 88]) (Figure 4D).
We compared the discriminatory performance of the six glycopeptides panel to that of the FDA approved ELISA assay for mesothelin (Mesomark®)  (Additional file 10: Figure S3). In the 75 subjects above (23 MPM, 26 HD and 26 NSCLC, Additional file 3: Table S2) for which ELISA measurements were available, mesothelin ELISA discriminated MPM from HD with AUC of 0.92 (95% CI, [0.83, 0.99]) and accuracy of 82% (95% CI, [71, 94]) at the 2 nM cut-off proposed in the literature . In the same subjects, which were part of the validation set above, the six glycopeptides panel based on SRM assay technology had AUC of 0.94 (95% CI, [0.86, 0.99]) and accuracy of 86% (95% CI, [71, 96]) at the 0.61 cut-off, indicating a discriminatory power similar to that of the ELISA assay. For the discrimination between MPM and NSCLC, mesothelin ELISA was superior to SRM and had AUC of 0.84 (95% CI, [0.71, 0.94]) whereas the six glycopeptides panel had AUC of 0.54 (95% CI, [0.37, 0.71]).
From this last observation, we argued that the specificity of the SRM panel for MPM could be increased by including mesothelin. Thus, we added the SRM monitoring of mesothelin peptide KWDVTSLETLK to the signature (Table 1 and Additional file 11: Table S7) and tested it in the group of 75 subjects for which mesothelin ELISA measurements were available for comparison. We performed parameterization of the signature using a sub-training set of 12 MPM, 14 HD and 14 NSCLC and assessed its performance in a sub-validation set of 11 MPM, 12 HD and 12 NSCLC (Figures 5A and B). Here, the SRM based signature discriminated MPM from HD with AUC of 0.95 (95% CI, [0.83, 1.0]) and highest accuracy of 91% (95% CI, [83, 100]) and MPM from NSCLC with AUC of 0.84 (95%CI, [0.66, 0.97]) and highest accuracy of 74% (95% CI, [57, 91]). The performance of the signature was higher than for mesothelin ELISA, which in the same subjects had AUC of 0.93 (95% CI, [0.84, 1.0]) and accuracy of 78% (95% CI, [61, 96]; cut-off 2 nM) in discriminating MPM from HD and AUC of 0.80 (95% CI, [0.59, 0.95]) and accuracy of 74% (95% CI, [57, 100]; cut-off 2 nM) in discriminating MPM from NSCLC.
The development of protein biomarkers in serum requires the availability of reliable analytical tools for the unbiased prioritization and large scale clinical testing of novel candidates [38, 39]. ELISA’s assays for biomarker investigation are typically obtainable only for a subset of candidates and establishing new ELISA’s solely for the purpose of testing candidates is time consuming and too expensive [40, 41]. In our study, we presented the application of targeted proteomics SRM assay technology in serum for the investigation and clinical evaluation of candidate biomarker panels for MPM. The approach presented itself as an accurate alternative to immunoassays and allowed us to follow a hypothesis driven biomarker investigation independently of the still piecemeal availability of antibodies.
The underlying hypothesis of our investigation was that the surfaceome of MPM can reveal novel blood accessible biomarkers. This was suggested by the fact that the cell surface proteins are exposed to the tumor environment and thus prone to be shed or released to the stroma and finally collected into the blood [42, 43]. Indeed, many proposed blood tumor-markers like mesothelin, carcinoembryonic antigen (CEA) or cancer antigen 125 (CA-125) are glycoproteins of cell surface origin. Following this hypothesis, our SRM investigation in serum detected 51 out of the 112 surfaceome derived candidate biomarkers. This was a considerable fraction of the candidates, considering that they were selected without prior knowledge and from extremely simplified tumor models as represented by cell lines. Interestingly, the majority of the candidates were at concentrations in the range of proposed MPM biomarkers like fibulin-3, megakaryocyte potentiating factor (MPF) or osteopontin. This observation, together with the SRM detection of mesothelin, confirmed that our targeted proteomics strategy could reliably access that fraction of the serum proteome which seems to be of particular relevance for MPM biomarker investigations.
In our study, if available, at least two surfaceome detected N-glycopeptides per protein were initially investigated by SRM in serum. For our final MPM signature in serum we selected the best peptide for a particular protein, e.g. these peptides were most consistently detected. Potential variations in the detectability and response factor of peptides from the same proteins are related to a number of reasons. One reason is certainly related to the fact that each peptide has its peculiar physicochemical characteristics which will influence its mass spectrometric detection independently of a common protein of origin [44, 45]. At the same time it is also likely that the signature reflects the complexity of a natural tissue environment [46, 47]. Indeed, several biochemical and proteolytic processes are expected to take place in the tumor microenvironment which can modify the original structure of the cell surface proteins [48, 49]. It is thus likely that not intact proteins but rather only fragments of them will reach and pass the vessel barriers [50, 51]. This could at least in part explain the apparently asynchronous behavior in serum of peptides from the same protein.
Despite the confident discrimination between MPM and healthy, our study cannot conclusively answer the question if the candidate biomarkers of the signature are MPM specific or rather more generally cancer associated. Indeed, without mesothelin, the six biomarkers of the signature failed to discriminate MPM from NSCLC and their association with other tumors is reported [52–56]. Nevertheless, the SRM signature inclusive of mesothelin presented accuracies higher than the ELISA test for the single marker mesothelin. This indicated that the integration of the seven MPM biomarkers in the multiplexed SRM signature could complement the limited sensitivity of mesothelin, taking at the same time advantage of its specificity for MPM. Here, we have to point out that, because of the exploratory nature of our investigation, the majority of the patients included in our study were at advanced disease stages and that controls did not include confounding conditions like chronic inflammations or other non-malignant pathologies of the lung. As a consequence, the accuracy of the MPM signature could be lower if applied to more heterogeneous populations, like indirectly suggested by the higher AUCs of mesothelin ELISA of our study in respect to literature reports .
Finally, it has to be highlighted that the MPM signature includes hypoxia up-regulated protein 1 (also known as ORP150 or GRP170; UniProt/Swiss-Prot: Q9Y4L1, gene name HYOU1) which is a heat shock protein with chaperone function in the endoplasmic reticulum . This could arise some concern about the specificity of our approach. It is therefore worthwhile to mention here that, in accord with other groups [58, 59], in our surfaceome experiments we reproducibly observed the protein and it is known that heat shock proteins can be expressed on cell surfaces or be secreted to blood [60–63].
In conclusion, the SRM assay technology based approach chosen for our clinical MPM investigation allowed us to directly evaluate a larger set of candidate serum biomarkers resulting in a seven glycopeptide signature with diagnostic potential for MPM. Our results indicate that the SRM assay technology lends itself for the fast clinical evaluation of candidate biomarkers in serum. In this respect, larger SRM-assays repositories are currently being generated [64, 65], which will ultimately enable the quantitative evaluation of biomarker candidates of interest in the disease setting of choice.
The MPM cell lines ZL55, SDM4, SDM5 and SDM34 were from surgical tumor samples and the pleural cell line SDM104 was from a surgical biopsy of a patient with chronic pleuritis. Cell lines were established as previously described [66, 67] and were from patients with pathologically confirmed diagnosis and treated at the University Hospital Zürich. HCC4012 was from human mesothelial cells immortalized with hTERT (kind gift of Dr. A. Gazdar, The University of Texas, Southwestern Medical Center). ADCA cell lines Calu-3 and SK-LU-1 were from American Type Culture Collection (ATCC; Manassas, VA). Detailed growing conditions can be found in Additional file 12: Supplementary Methods.
CSC-based surfaceome analysis and MPM candidate biomarkers selection
CSC followed by MS analysis was performed as described previously . For label free relative-quantification, raw data of duplicate measurements were acquired in profile mode on a Fourier-Transform LTQ MS (FT-LTQ; Thermo Electron, San Jose, CA), converted to mzXML  and analyzed with the software Superhirn . For sequence identification MS/MS spectra of centroided raw files were converted to mzXML and searched against the IPI Human database v3.26 using the search algorithm SEQUEST v27 . Criteria for MPM candidate biomarker peptides were: 1. fully tryptic. 2. deamidation of asparagine in the consensus sequence NxS/T (x denotes any amino acid excluded proline) after treatment with PNGaseF. 3. PeptideProphet probability ≥ 0.9. 4. sequence proteotypic and unique for proteins reviewed in Uniprot  and with subcellular localization associated to membranes or secreted. 5. reproducibly higher abundant in MPM in at least two MPM vs non-MPM cell lines comparisons, or originating from the same protein of an higher abundant peptide, or deriving from a protein not observed in non-MPM cell lines but detected in MPM at least in two cell lines or with two peptides. Further details about quantitative CSC analysis are reported in Additional file 12: Supplementary Methods.
Generation of SRM-assays
To establish glycopeptide-specific SRM assays, spectra of MPM candidate biomarker glycopeptides were generated by using synthetic heavy isotope-labeled (heavy, with R 13C6/15 N4 and/or K 13C6/15 N2) peptides (SpikeTides_L™, JPT Peptide Technologies, Berlin, Germany) with aspartic acid (D) replacing the putative glycosylated asparagines (N) according to the mass modification introduced by treatment with the enzyme PNGaseF in the protocol for enrichment of N-glycopeptides from serum. Spectra were acquired on Quadrupole Time-of-Flight (QTOF) LC/MS series 6520 or 6550 instruments (Agilent Technologies, Santa Clara, CA) equipped with an HPLC-Chip Cube interface (Agilent Technologies) and operated in data dependent mode. MS/MS spectra were used to generate initial SRM-assays for MPM candidate biomarkers. They consisted of at least six transitions per peptide selected based on signal intensities of heavy peptides (SpikeTides_L™, JPT Peptide Technologies) spiked in the matrix of enriched serum. SRM-assays of candidate biomarkers detected in serum were further individually optimized and consisted of four transitions per peptide with at least three fixed transitions used for quantification. Details about spectra acquisition, MS settings, SRM-assays generation and optimization can be found in Additional file 12: Supplementary Methods. All assays developed can also be downloaded in form of a Skyline library file (Additional file 5: Skyline file).
Whole blood samples were obtained after written informed consent from therapy naïve patients with pathologically proven diagnosis of MPM or NSCLC and treated at the University Hospital Zürich. Staging was based on TNM-International Union Against Cancer (UICC, sixth edition) selecting the highest stage in case of ambiguous report. Whole blood samples from HD were from blood donors at the Blood Transfusion Service Zürich, SRC, Schlieren, Switzerland and judged healthy based on standardized medical questionnaire . The study was approved by the Ethics Committee of the University Hospital Zürich. Serum processing is reported in Additional file 12: Supplementary Methods.
Serum enrichment for N-glycopeptides and MS analysis
For SRM analysis, 100 μl of serum were enriched for N-glycopeptides using a modified protocol of the method for solid phase extraction of N-glycopeptides (SPEG) . 1.5 μl of peptide mixture were analyzed on a QQQ LC/MS 6460 series (Agilent Technologies) equipped with an HPLC-Chip Cube interface (Agilent technologies) and using a nano-flow gradient of 5 to 35% acetonitrile (ACN) /water, 0.1% formic acid (FA) over 30 min. The software Skyline  was used for SRM-traces visualization after Savitzky-Golay smoothing, SRM-methods building and calculation of peak transition-intensities. Details about serum processing and MS settings can be found in Additional file 12: Supplementary Methods.
Verification of MPM candidate biomarker peptides in serum
To verify the detectability of MPM candidate biomarker peptides in serum, samples from five MPM subjects were enriched for N-glycopeptides and analyzed on a QQQ LC/MS instrument using not-optimized SRM-assays. Sample processing and MS settings were as described above. Transitions were monitored in scheduled SRM-mode allowing for a maximum of 339 total transitions and 176 concurrent transitions per method. Cycle-times ranged from 2 to 4.1 s allowing for a minimal dwell time of 18.5 ms per transition. Delta retention time window was 4 or 5 min. Confident detection of MPM candidate biomarker peptides in serum was manually confirmed based on transition co-elution with simultaneously monitored heavy isotope-labeled synthetic peptides with matching sequences (SpikeTides_L™, JPT Peptide Technologies) spiked in the samples before MS analysis.
SRM analysis of candidate biomarker N-glycopeptides from clinical cohorts
Serum samples of training and validation sets were enriched for N-glycopeptides and analyzed using optimized SRM assays on a QQQ LC/MS instrument as described above. The two sets were processed and analyzed at separate time points. Samples of the same set were processed simultaneously in randomized order and analyzed in technical duplicates on the QQQ. Eleven samples (normalizing-samples) from the training set were re-processed and re-analyzed in parallel with the validation set and results were used for normalization of SRM signals between the two groups. These samples were subsequently excluded from the validation set. For relative quantification, a mix of heavy isotope-labeled synthetic peptides with sequences matching the MPM candidate biomarker peptides was used as internal standard (SpikeTides_L™, JPT Peptide Technologies, for mesothelin heavy isotope-labeled synthetic peptides were from Thermo Scientific) and spiked at fixed concentration in each sample before MS analysis at a volume ratio of 1:5 of heavy-peptide-mix to serum sample. To assess technical variations among runs, iRT peptides (Biognosys, Schlieren, Switzerland)  were spiked in each sample before MS analysis. MS analysis of serum samples from the training set was performed using a scheduled SRM method including a total of 468 light and heavy transitions. Cycle time was of 3.7 s allowing for the acquisition of at least eight data points per peptide elution profile. RT window was set to 5 min. Dwell time per transition ranged from a minimum of 16 ms to a maximum of 459 ms. The number of concurrent transitions ranged from 8 to maximal 190. Samples of the validation set were analyzed in scheduled SRM-mode monitoring for a total of 288 light and heavy transitions. Cycle time was set to 3 s for the acquisition of at least eight data points per peptide elution profile using a delta RT window of 5 min. Minimal dwell time per transition was 26.6 ms and maximal was 459 ms. Minimal and maximal number of concurrent transitions were 8 and 123 respectively. Both method included transitions from the iRT peptides and peptides of the serum proteins haptoglobin (UniProt entry P00738) and kininogen-1 (UniProt entry P01042) used as internal reference control for sample handling and MS performance. Confident detection of MPM candidate biomarker peptides was confirmed manually based on transition co-elution with heavy isotope-labeled internal standards.
Statistical significance analysis and prediction analysis
Statistical analysis of peptide differential abundance utilized SRMstats package in R [37, 75]. Ten peptides of higher abundance in training set in either comparison for MPM vs. HD or MPM vs. NSCLC were further used in two logistic regression models for MPM vs HD and MPM vs NSCLC. In order to account for relative experimental yield and reproducibility of sample preparation between training and validation sets, we developed a two-step normalization procedure based on the eleven normalizing-samples that were present in both sets. The first normalization step accounted for variations in the mass spectrometer performance, separately for the training and the validation sets, by equating median intensities of reference transitions between the runs. The second normalization step shifted the intensities of the endogenous transitions in the validation set to the scale of the training set. Specifically, for each endogenous transition we calculated the median difference of log-intensities among the eleven normalizing-samples in the validation and the training sets. The difference was then subtracted from the endogenous intensities in all the validation samples. All inputs for the logistic regressions are estimates of peptide abundance in each biological sample on a relative scale, which are summarized across multiple transitions and technical replicate runs. This summarization was performed in SRMstats fitting logistic regression in R. 'pROC’ package in R was used to draw ROC curves and to calculate AUCs and CI with bootstrap methods . Correlations and Mann-Whithney test were calculated and visualized using IBM SPSS Statistics Standard v17.0 (SPSS, Inc, Chicago, IL) or GraphPad Prism 5 (GraphPad Software, Inc, San Diego, CA).
Mesothelin ELISA in serum was performed in duplicates using the Mesomark-kit™ (Fujirebio Diagnostic, Malvern, PA) according to the manufacturer’s protocol. Averaged values were used for analysis. Samples with coefficient of variation > 15% were excluded.
Sugarbaker DJ, Flores RM, Jaklitsch MT, Richards WG, Strauss GM, Corson JM, DeCamp MM, Swanson SJ, Bueno R, Lukanich JM: Resection margins, extrapleural nodal status, and cell type determine postoperative long-term survival in trimodality therapy of malignant pleural mesothelioma: results in 183 patients. J Thorac Cardiovasc Surg. 1999, 117: 54-63. discussion 63–55,
Weder W, Kestenholz P, Taverna C, Bodis S, Lardinois D, Jerman M, Stahel RA: Neoadjuvant chemotherapy followed by extrapleural pneumonectomy in malignant pleural mesothelioma. J Clin Oncol. 2004, 22: 3451-3457.
Flores RM, Krug LM, Rosenzweig KE, Venkatraman E, Vincent A, Heelan R, Akhurst T, Rusch VW: Induction chemotherapy, extrapleural pneumonectomy, and postoperative high-dose radiotherapy for locally advanced malignant pleural mesothelioma: a phase II trial. J Thorac Oncol. 2006, 1: 289-295.
Weder W, Stahel RA, Bernhard J, Bodis S, Vogt P, Ballabeni P, Lardinois D, Betticher D, Schmid R, Stupp R: Multicenter trial of neo-adjuvant chemotherapy followed by extrapleural pneumonectomy in malignant pleural mesothelioma. Ann Oncol. 2007, 18: 1196-1202.
Krug LM, Pass HI, Rusch VW, Kindler HL, Sugarbaker DJ, Rosenzweig KE, Flores R, Friedberg JS, Pisters K, Monberg M: Multicenter phase II trial of neoadjuvant pemetrexed plus cisplatin followed by extrapleural pneumonectomy and radiation for malignant pleural mesothelioma. J Clin Oncol. 2009, 27: 3007-3013.
Wilcox BE, Subramaniam RM, Peller PJ, Aughenbaugh GL, Nichols Iii FC, Aubry MC, Jett JR: Utility of integrated computed tomography-positron emission tomography for selection of operable malignant pleural mesothelioma. Clin Lung Cancer. 2009, 10: 244-248.
Sorensen JB, Ravn J, Loft A, Brenoe J, Berthelsen AK: Preoperative staging of mesothelioma by 18F-fluoro-2-deoxy-D-glucose positron emission tomography/computed tomography fused imaging and mediastinoscopy compared to pathological findings after extrapleural pneumonectomy. Eur J Cardiothorac Surg. 2008, 34: 1090-1096.
de Perrot M, Feld R, Cho BC, Bezjak A, Anraku M, Burkes R, Roberts H, Tsao MS, Leighl N, Keshavjee S, Johnston MR: Trimodality therapy with induction chemotherapy followed by extrapleural pneumonectomy and adjuvant high-dose hemithoracic radiation for malignant pleural mesothelioma. J Clin Oncol. 2009, 27: 1413-1418.
Robinson BW, Creaney J, Lake R, Nowak A, Musk AW, de Klerk N, Winzell P, Hellstrom KE, Hellstrom I: Mesothelin-family proteins and diagnosis of mesothelioma. Lancet. 2003, 362: 1612-1616.
Hollevoet K, Reitsma JB, Creaney J, Grigoriu BD, Robinson BW, Scherpereel A, Cristaudo A, Pass HI, Nackaerts K, Rodriguez Portal JA: Serum mesothelin for diagnosing malignant pleural mesothelioma: an individual patient data meta-analysis. J Clin Oncol. 2012, 30: 1541-1549.
Pass HI, Levin SM, Harbut MR, Melamed J, Chiriboga L, Donington J, Huflejt M, Carbone M, Chia D, Goodglick L: Fibulin-3 as a blood and effusion biomarker for pleural mesothelioma. N Engl J Med. 2012, 367: 1417-1427.
Fu Q, Schoenhoff FS, Savage WJ, Zhang P, Van Eyk JE: Multiplex assays for biomarker research and clinical application: translational science coming of age. Proteomics Clin Appl. 2010, 4: 271-284.
Picotti P, Bodenmiller B, Aebersold R: Proteomics meets the scientific method. Nat Methods. 2013, 10: 24-27.
Lange V, Picotti P, Domon B, Aebersold R: Selected reaction monitoring for quantitative proteomics: a tutorial. Mol Syst Biol. 2008, 4: 222
Addona TA, Abbatiello SE, Schilling B, Skates SJ, Mani DR, Bunk DM, Spiegelman CH, Zimmerman LJ, Ham AJ, Keshishian H: Multi-site assessment of the precision and reproducibility of multiple reaction monitoring-based measurements of proteins in plasma. Nat Biotechnol. 2009, 27: 633-641.
Elschenbroich S, Ignatchenko V, Clarke B, Kalloger S, Boutros P, Gramolini A, Shaw P, Jurisica I, Kislinger T: In-depth proteomics of ovarian cancer ascites: combining shotgun proteomics and selected reaction monitoring mass spectrometry. J Proteome Res. 2011, 10: 2286-2299.
Whiteaker JR, Lin C, Kennedy J, Hou L, Trute M, Sokal I, Yan P, Schoenherr RM, Zhao L, Voytovich UJ: A targeted proteomics-based pipeline for verification of biomarkers in plasma. Nat Biotechnol. 2011, 29: 625-634.
Ziegler A, Cerciello F, Bigosch C, Bausch-Fluck D, Felley-Bosco E, Ossola R, Soltermann A, Stahel RA, Wollscheid B: Proteomic surfaceome analysis of mesothelioma. Lung Cancer. 2012, 75: 189-196.
Anderson NL, Anderson NG: The human plasma proteome: history, character, and diagnostic prospects. Mol Cell Proteomics. 2002, 1: 845-867.
Zhang H, Li XJ, Martin DB, Aebersold R: Identification and quantification of N-linked glycoproteins using hydrazide chemistry, stable isotope labeling and mass spectrometry. Nat Biotechnol. 2003, 21: 660-666.
Tian Y, Zhou Y, Elliott S, Aebersold R, Zhang H: Solid-phase extraction of N-linked glycopeptides. Nat Protoc. 2007, 2: 334-339.
Zhou Y, Aebersold R, Zhang H: Isolation of N-linked glycopeptides from plasma. Anal Chem. 2007, 79: 5826-5837.
Stahl-Zeng J, Lange V, Ossola R, Eckhardt K, Krek W, Aebersold R, Domon B: High sensitivity detection of plasma proteins by multiple reaction monitoring of N-glycosites. Mol Cell Proteomics. 2007, 6: 1809-1817.
Picotti P, Rinner O, Stallmach R, Dautel F, Farrah T, Domon B, Wenschuh H, Aebersold R: High-throughput generation of selected reaction-monitoring assays for proteins and proteomes. Nat Methods. 2010, 7: 43-46.
Chang K, Pastan I: Molecular cloning of mesothelin, a differentiation antigen present on mesothelium, mesotheliomas, and ovarian cancers. Proc Natl Acad Sci U S A. 1996, 93: 136-140.
Hellstrom I, Raycraft J, Kanan S, Sardesai NY, Verch T, Yang Y, Hellstrom KE: Mesothelin variant 1 is released from tumor cells as a diagnostic marker. Cancer Epidemiol Biomarkers Prev. 2006, 15: 1014-1020.
Ho M, Onda M, Wang QC, Hassan R, Pastan I, Lively MO: Mesothelin is shed from tumor cells. Cancer Epidemiol Biomarkers Prev. 2006, 15: 1751
Beyer HL, Geschwindt RD, Glover CL, Tran L, Hellstrom I, Hellstrom KE, Miller MC, Verch T, Allard WJ, Pass HI, Sardesai NY: MESOMARK: a potential test for malignant pleural mesothelioma. Clin Chem. 2007, 53: 666-672.
Scholler N, Fu N, Yang Y, Ye Z, Goodman GE, Hellstrom KE, Hellstrom I: Soluble member(s) of the mesothelin/megakaryocyte potentiating factor family are detectable in sera from patients with ovarian carcinoma. Proc Natl Acad Sci U S A. 1999, 96: 11531-11536.
Wollscheid B, Bausch-Fluck D, Henderson C, O’Brien R, Bibel M, Schiess R, Aebersold R, Watts JD: Mass-spectrometric identification and relative quantification of N-linked cell surface glycoproteins. Nat Biotechnol. 2009, 27: 378-386.
de Graaf E, Altelaar A, van Breukelen B, Mohammed S, Heck A: Improving SRM assay development: a global comparison between triple quadrupole, ion trap, and higher energy CID peptide fragmentation spectra. J Proteome Res. 2011, 10: 4334-4341.
Farrah T, Deutsch EW, Omenn GS, Campbell DS, Sun Z, Bletz JA, Mallick P, Katz JE, Malmstrom J, Ossola R: A high-confidence human plasma proteome reference set with estimated concentrations in PeptideAtlas. Mol Cell Proteomics. 2012, 10: M110 006353
Pass HI, Lott D, Lonardo F, Harbut M, Liu Z, Tang N, Carbone M, Webb C, Wali A: Asbestos exposure, pleural mesothelioma, and serum osteopontin levels. N Engl J Med. 2005, 353: 1564-1573.
Grigoriu B, Chahine B, Zerimech F, Gregoire M, Balduyck M, Copin MC, Devos P, Lassalle P, Scherpereel A: Serum mesothelin has a higher diagnostic utility than hyaluronic acid in malignant mesothelioma. Clin Biochem. 2009, 42: 1046-1050.
Hollevoet K, Nackaerts K, Thimpont J, Germonpre P, Bosquee L, De Vuyst P, Legrand C, Kellen E, Kishi Y, Delanghe JR, van Meerbeeck JP: Diagnostic performance of soluble mesothelin and megakaryocyte potentiating factor in mesothelioma. Am J Respir Crit Care Med. 2010, 181: 620-625.
Schouwink H, Korse CM, Bonfrer JM, Hart AA, Baas P: Prognostic value of the serum tumour markers Cyfra 21–1 and tissue polypeptide antigen in malignant mesothelioma. Lung Cancer. 1999, 25: 25-32.
Chang CY, Picotti P, Hüttenhain R, Heinzelmann-Schwarz V, Jovanovic M, Aebersold R, Vitek O: Protein significance analysis in selected reaction monitoring (SRM) measurements. Mol Cell Proteomics. 2012, 11: M111 014662
Pepe MS, Etzioni R, Feng Z, Potter JD, Thompson ML, Thornquist M, Winget M, Yasui Y: Phases of biomarker development for early detection of cancer. J Natl Cancer Inst. 2001, 93: 1054-1061.
Zhang Z, Chan DW: The road from discovery to clinical diagnostics: lessons learned from the first FDA-cleared in vitro diagnostic multivariate index assay of proteomic biomarkers. Cancer Epidemiol Biomarkers Prev. 2010, 19: 2995-2999.
Elschenbroich S, Kislinger T: Targeted proteomics by selected reaction monitoring mass spectrometry: applications to systems biology and biomarker discovery. Mol Biosyst. 2011, 7: 292-303.
Rifai N, Gillette MA, Carr SA: Protein biomarker discovery and validation: the long and uncertain path to clinical utility. Nat Biotechnol. 2006, 24: 971-983.
Roth J, Zuber C, Park S, Jang I, Lee Y, Kysela K, Fourn VL, Santimaria R, Guhl B, Cho J: Protein N-glycosylation, protein folding, and protein quality control. Mol Cells. 2010, 30: 497-506.
Zhang H, Liu AY, Loriaux P, Wollscheid B, Zhou Y, Watts JD, Aebersold R: Mass spectrometric detection of tissue proteins in plasma. Mol Cell Proteomics. 2007, 6: 64-71.
Fusaro VA, Mani DR, Mesirov JP, Carr SA: Prediction of high-responding peptides for targeted protein assays by mass spectrometry. Nat Biotechnol. 2009, 27: 190-198.
Duncan MW, Aebersold R, Caprioli RM: The pros and cons of peptide-centric proteomics. Nat Biotechnol. 2010, 28: 659-664.
Fang Q, Kani K, Faca VM, Zhang W, Zhang Q, Jain A, Hanash S, Agus DB, McIntosh MW, Mallick P: Impact of protein stability, cellular localization, and abundance on proteomic detection of tumor-derived proteins in plasma. PLoS One. 2011, 6: e23090
Wang J, Willumsen N, Zheng Q, Xue Y, Karsdal MA, Bay-Jensen AC: Bringing cancer serological diagnosis to a new level: focusing on HER2, protein ectodomain shedding and neoepitope technology. Future Oncol. 2013, 9: 35-44.
Craig SE, Brady-Kalnay SM: Cancer cells cut homophilic cell adhesion molecules and run. Cancer Res. 2011, 71: 303-309.
Mason SD, Joyce JA: Proteolytic networks in cancer. Trends Cell Biol. 2011, 21: 228-237.
Wilken JA, Perez-Torres M, Nieves-Alicea R, Cora EM, Christensen TA, Baron AT, Maihle NJ: Shedding of soluble epidermal growth factor receptor (sEGFR) is mediated by a metalloprotease/fibronectin/integrin axis and inhibited by cetuximab. Biochemistry. 2013, 52: 4531-4540.
Ebos JM, Bocci G, Man S, Thorpe PE, Hicklin DJ, Zhou D, Jia X, Kerbel RS: A naturally occurring soluble form of vascular endothelial growth factor receptor 2 detected in mouse and human plasma. Mol Cancer Res. 2004, 2: 315-326.
Shah N, Cabanillas F, McIntyre B, Feng L, McLaughlin P, Rodriguez MA, Romaguera J, Younes A, Hagemeister FB, Kwak L, Fayad L: Prognostic value of serum CD44, intercellular adhesion molecule-1 and vascular cell adhesion molecule-1 levels in patients with indolent non-Hodgkin lymphomas. Leuk Lymphoma. 2012, 53: 50-56.
Iozzo RV, Sanderson RD: Proteoglycans in cancer biology, tumour microenvironment and angiogenesis. J Cell Mol Med. 2011, 15: 1013-1031.
Nanda A, Carson-Walter EB, Seaman S, Barber TD, Stampfl J, Singh S, Vogelstein B, Kinzler KW, Croix BS: TEM8 Interacts with the cleaved C5 domain of collagen α3(VI). Cancer Res. 2004, 64: 817-820.
Miyagi T, Hori O, Koshida K, Egawa M, Kato H, Kitagawa Y, Ozawa K, Ogawa S, Namiki M: Antitumor effect of reduction of 150-kDa oxygen-regulated protein expression on human prostate cancer cells. Int J Urol. 2002, 9: 577-585.
Isenberg JS, Martin-Manso G, Maxhimer JB, Roberts DD: Regulation of nitric oxide signalling by thrombospondin 1: implications for anti-angiogenic therapies. Nat Rev Cancer. 2009, 9: 182-194.
Easton DP, Kaneko Y, Subjeck JR: The Hsp110 and Grp170 stress proteins: newly recognized relatives of the Hsp70s. Cell Stress Chaperones. 2000, 5: 276-290.
Naaby-Hansen S, Herr JC: Heat shock proteins on the human sperm surface. J Reprod Immunol. 2010, 84: 32-40.
Calvert ME, Digilio LC, Herr JC, Coonrod SA: Oolemmal proteomics – identification of highly abundant heat shock proteins and molecular chaperones in the mature mouse egg and their localization on the plasma membrane. Reprod Biol Endocrinol. 2003, 1: 27-37.
Beachy SH, Kisailus AJ, Repasky EA, Subjeck JR, Wang XY, Kazim AL: Engineering secretable forms of chaperones for immune modulation and vaccine development. Methods. 2007, 43: 184-193.
Altmeyer A, Maki R, Feldweg A, Heike M, Protopopov V, Masur S, Srivastava P: Tumor-specific cell surface expression of the-KDEL containing, endoplasmic reticular heat shock protein gp96. Int J Cancer. 1996, 69: 340-349.
Gastpar R, Gross C, Rossbacher L, Ellwart J, Riegger J, Multhoff G: The cell surface-localized heat shock protein 70 epitope TKD induces migration and cytolytic activity selectively in human NK cells. J Immunol. 2004, 172: 972-980.
Vega VL, Rodrı’guez-Silva M, Frey T, Gehrmann M, Diaz JC, Steinem C, Multhoff G, Arispe N, Maio AD: Hsp70 Translocates into the plasma membrane after stress and is released into the extracellular environment in a membrane-associated form that activates macrophages. J Immunol. 2008, 180: 4299-4307.
Hüttenhain R, Soste M, Selevsek N, Röst H, Sethi A, Carapito C, Farrah T, Deutsch EW, Kusebauch U, Moritz RL: Reproducible quantification of cancer-associated proteins in body fluids using targeted proteomics. Sci Transl Med. 2012, 4: 142ra194
The SRMAtlas. http://www.srmatlas.org,
Thurneysen C, Opitz I, Kurtz S, Weder W, Stahel RA, Felley-Bosco E: Functional inactivation of NF2/merlin in human mesothelioma. Lung Cancer. 2009, 64: 140-147.
Knobel PA, Kotov IN, Felley-Bosco E, Stahel RA, Marti TM: Inhibition of REV3 expression induces persistent DNA damage and growth arrest in cancer cells. Neoplasia. 2011, 13: 961-970.
Pedrioli PG, Eng JK, Hubley R, Vogelzang M, Deutsch EW, Raught B, Pratt B, Nilsson E, Angeletti RH, Apweiler R: A common open representation of mass spectrometry data and its application to proteomics research. Nat Biotechnol. 2004, 22: 1459-1466.
Mueller LN, Rinner O, Schmidt A, Letarte S, Bodenmiller B, Brusniak MY, Vitek O, Aebersold R, Müller M: SuperHirn - a novel tool for high resolution LC-MS-based peptide/protein profiling. Proteomics. 2007, 7: 3470-3480.
Eng JK, McCormack AL, Yates JR: An approach to correlate tandem mass spectral data of peptides with amino acid sequences in a protein database. J Am Soc Mass Spectrom. 1994, 5: 976-989.
The UniProt Knowledgebase. (UniProtKB) http://www.uniprot.org,
Blood Transfusion Service Zürich, SRC, Switzerland. http://www.zhbsd.ch,
MacLean B, Tomazela DM, Shulman N, Chambers M, Finney GL, Frewen B, Kern R, Tabb DL, Liebler DC, MacCoss MJ: Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics. 2010, 26: 966-968.
Escher C, Reiter L, Maclean B, Ossola R, Herzog F, Chilton J, Maccoss MJ, Rinner O: Using iRT, a normalized retention time for more targeted measurement of peptides. Proteomics. 2012, 12: 1111-1121.
Laboratory for Statistical Proteomics and Bioinformatics. http://www.stat.purdue.edu/~ovitek/Home.html,
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, Müller M: pROC: an open-source package for R and S + to analyze and compare ROC curves. BMC Bioinforma. 2011, 12: 77
We are grateful to the Wollscheid and Aebersold groups, to the Laboratory of Molecular Oncology and to Dr. Ralph Schiess and the team of ProteoMediX AG (Schlieren, Switzerland) for critical discussion and support. We are grateful to Colette Bigosch and Ubiratan Moura for their contribution to sample processing. FC was recipient of a fellowship under the Swiss national MD-PhD program of the Swiss National Science Foundation (SNSF) and enrolled in the MD-PhD and Cancer Biology program of the University of Zürich. This work was supported by the Cancer League of Zürich, the Huggenberger-Bischoff foundation for Cancer Research, the NCCR Neural Plasticity and Repair of the SNSF (to BW), BIP SystemsX.ch (to BW). The project was further supported in part by the SNSF (Grant# 3100A0-107679), the European Research Council (Grant# ERC-2008-AdG 233226) and by grant Prime-X (EC) to RA.
The authors declare that they have no competing interests.
FC, RA, EFB, RS and BW conceived the project, designed and interpreted experiments. FC and AN performed experiments. DBF provided support and interpreted CSC experiments. MC and OV performed statistical analysis in clinical serum cohorts. FC, EFB, RS and AZ coordinated clinical sample acquisition, FC and AZ processed and registered samples. FC, MC and BW wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3: Table S2: SRM and ELISA quantitative investigation of mesothelin in serum samples of 75 subjects. (XLS 156 KB)
Additional file 4: Table S3: A. Surfaceome peptides fulfilling requirements for MPM candidate biomarker. B. MPM candidate biomarkers selected for SRM analysis in serum. (XLS 128 KB)
Additional file 5: (Skyline, https://brendanx-uw1.gs.washington.edu/labkey/project/home/software/Skyline/begin.view). Skyline library file and serum verification of SRM-assays developed for surfaceome derived MPM candidate biomarkers. (ZIP 7 MB)
Additional file 6: Figure S2: MPM candidate biomarker peptides detected by SRM in the screening of five MPM sera after enrichment for N-glycopeptides. (PDF 1010 KB)
Additional file 7: Table S4: MPM candidate biomarker peptides identified by SRM in sera enriched for N-glycopeptides. (XLS 78 KB)
Additional file 8: Table S5: Clinical characteristics of the 55 MPM, 53 NSCLC and 54 healthy donors in the training and validation set. (XLS 93 KB)
About this article
Cite this article
Cerciello, F., Choi, M., Nicastri, A. et al. Identification of a seven glycopeptide signature for malignant pleural mesothelioma in human serum by selected reaction monitoring. Clin Proteom 10, 16 (2013). https://doi.org/10.1186/1559-0275-10-16
- Malignant pleural mesothelioma
- Selected reaction monitoring
- Targeted proteomics
- Serum biomarkers