Proteomic analysis of serum extracellular vesicles reveals Fibulin-3 as a new marker predicting liver-related events in MASLD.

BACKGROUND: There is a need for novel noninvasive markers for metabolic dysfunction-associated steatotic liver disease (MASLD) to stratify patients at high risk for liver-related events including liver cancer and decompensation. In the present study, we used proteomic analysis of proteins in extracellular vesicles (EVs) to identify new biomarkers that change with fibrosis progression and can predict the development of liver-related events. METHODS: We analyzed serum EVs from 50 patients with MASLD assessed for liver fibrosis by biopsy and identified proteins that altered with advanced fibrosis. A further evaluation was conducted on another cohort of 463 patients with MASLD with biopsy. RESULTS: Eight candidate proteins were identified by proteomic analysis of serum EVs. Among them, serum levels of Fibulin-3, Fibulin-1, and Ficolin 1 correlated with their EV levels. In addition, serum Fibulin-3 and serum Fibulin-1 levels changed significantly with advanced fibrosis. Using another cohort with biopsy, we found that the serum Fibulin-3 concentration was significantly greater in those with advanced fibrosis but that the serum Fibulin-1 concentration was not significantly different. Multivariate Cox proportional hazards analysis revealed that a higher Fibrosis-4 (FIB-4) index and higher serum Fibulin-3 concentration were independent risk factors for liver-related events. When the cutoff value for the serum Fibulin-3 concentration was 6.0 µg/mL according to the Youden index of AUROCs, patients with high serum Fibulin-3 significantly more frequently developed liver-related events than did other patients. Validation using another cohort of 226 patients with clinically diagnosed MASLD confirmed that high serum Fibulin-3 levels are associated with a greater frequency of liver-related events. CONCLUSIONS: Serum Fibulin-3 was identified as a biomarker for predicting liver-related events in patients with MASLD.

Integrative analysis of cancer dependency data and comprehensive phosphoproteomics data revealed the EPHA2-PARD3 axis as a cancer vulnerability in KRAS-mutant colorectal cancer.

Colorectal cancer (CRC), a common malignant tumour of the gastrointestinal tract, is a life-threatening cancer worldwide. Mutations in KRAS and BRAF, the major driver mutation subtypes in CRC, activate the RAS pathway, contribute to tumorigenesis in CRC and are being investigated as potential therapeutic targets. Despite recent advances in clinical trials targeting KRASG12C or RAS downstream signalling molecules for KRAS-mutant CRC, there is a lack of effective therapeutic interventions. Therefore, understanding the unique molecular characteristics of KRAS-mutant CRC is essential for identifying molecular targets and developing novel therapeutic interventions. We obtained in-depth proteomics and phosphoproteomics quantitative data for over 7900 proteins and 38 700 phosphorylation sites in cells from 35 CRC cell lines and performed informatic analyses, including proteomics-based coexpression analysis and correlation analysis between phosphoproteomics data and cancer dependency scores of the corresponding phosphoproteins. Our results revealed novel dysregulated protein-protein associations enriched specifically in KRAS-mutant cells. Our phosphoproteomics analysis revealed activation of EPHA2 kinase and downstream tight junction signalling in KRAS-mutant cells. Furthermore, the results implicate the phosphorylation site Y378 in the tight junction protein PARD3 as a cancer vulnerability in KRAS-mutant cells. Together, our large-scale phosphoproteomics and proteomics data across 35 steady-state CRC cell lines represent a valuable resource for understanding the molecular characteristics of oncogenic mutations. Our approach to predicting cancer dependency from phosphoproteomics data identified the EPHA2-PARD3 axis as a cancer vulnerability in KRAS-mutant CRC.

Integration of pharmacoproteomic and computational approaches reveals the cellular signal transduction pathways affected by apatinib in gastric cancer cell lines.

Apatinib is known to be a highly selective vascular endothelial growth factor receptor 2 (VEGFR2) inhibitor with anti-angiogenic and anti-tumor properties. In a phase III study, the objective response rate to apatinib was low. It remains unclear why the effectivity of apatinib varies among patients and what type of patients are candidates for the treatment. In this study, we investigated the anti-tumor efficacy of apatinib against 13 gastric cancer cell lines and found that it differed depending on the cell line. Using integrated wet and dry approaches, we showed that apatinib was a multi-kinase inhibitor of c-Kit, RAF1, VEGFR1, VEGFR2, and VEGFR3, predominantly inhibiting c-Kit. Notably, KATO-III, which was the most apatinib-sensitive among the gastric cancer cell lines investigated, was the only cell line expressing c-Kit, RAF1, VEGFR1, and VEGFR3 but not VEGFR2. Furthermore, we identified SNW1 as a molecule affected by apatinib that plays an important role in cell survival. Finally, we identified the molecular network related to SNW1 that was affected by treatment with apatinib. These results suggest that the mechanism of action of apatinib in KATO-III cells is independent of VEGFR2 and that the differential efficacy of apatinib was due to differences in expression patterns of receptor tyrosine kinases. Furthermore, our results suggest that the differential efficacy of apatinib in gastric cell lines may be attributed to SNW1 phosphorylation levels at a steady state. These findings contribute to a deeper understanding of the mechanism of action of apatinib in gastric cancer cells.

A large-scale targeted proteomics of plasma extracellular vesicles shows utility for prognosis prediction subtyping in colorectal cancer.

PURPOSE: The pathogenesis of cancers depends on the molecular background of each individual patient. Therefore, verifying as many biomarkers as possible and clarifying their relationships with each disease status would be very valuable. We performed a large-scale targeted proteomics analysis of plasma extracellular vesicles (EVs) that may affect tumor progression and/or therapeutic resistance. EXPERIMENTAL DESIGN: Plasma EVs from 59 were collected patients with colorectal cancer (CRC) and 59 healthy controls (HC) in cohort 1, and 150 patients with CRC in cohort 2 for the large-scale targeted proteomics analysis of 457 proteins as candidate CRC markers. The Mann-Whitney-Wilcoxon test and random forest model were applied in cohort 1 to select promising markers. Consensus clustering was applied to classify patients with CRC in cohort 2. The Kaplan-Meier method and Cox regression analysis were performed to identify potential molecular factors contributing to the overall survival (OS) of patients. RESULTS: In the analysis of cohort 1, 99 proteins were associated with CRC. The analysis of cohort 2 revealed two clusters showing significant differences in OS (p = 0.017). Twelve proteins, including alpha-1-acid glycoprotein 1 (ORM1), were suggested to be associated with the identified CRC subtypes, and ORM1 was shown to significantly contribute to OS, suggesting that ORM1 might be one of the factors closely related to the OS. CONCLUSIONS: The study identified two novel subtypes of CRC, which exhibit differences in OS, as well as important biomarker proteins that are closely related to the identified subtypes. Liquid biopsy assessment with targeted proteomics analysis was proposed to be crucial for predicting the CRC prognosis.

EphA2 on urinary extracellular vesicles as a novel biomarker for bladder cancer diagnosis and its effect on the invasiveness of bladder cancer.

BACKGROUND: Urinary extracellular vesicles (uEVs) secreted from bladder cancer contain cancer-specific proteins that are potential diagnostic biomarkers. We identified and evaluated a uEV-based protein biomarker for bladder cancer diagnosis and analysed its functions. METHODS: Biomarker candidates, selected by shotgun proteomics, were validated using targeted proteomics of uEVs obtained from 49 patients with and 48 individuals without bladder cancer, including patients with non-malignant haematuria. We developed an enzyme-linked immunosorbent assay (ELISA) for quantifying the uEV protein biomarker without ultracentrifugation and evaluated urine samples from 36 patients with and 36 patients without bladder cancer. RESULTS: Thirteen membrane proteins were significantly upregulated in the uEVs from patients with bladder cancer in shotgun proteomics. Among them, eight proteins were validated by target proteomics, and Ephrin type-A receptor 2 (EphA2) was the only protein significantly upregulated in the uEVs of patients with bladder cancer, compared with that of patients with non-malignant haematuria. The EV-EphA2-CD9 ELISA demonstrated good diagnostic performance (sensitivity: 61.1%, specificity: 97.2%). We showed that EphA2 promotes proliferation, invasion and migration and EV-EphA2 promotes the invasion and migration of bladder cancer cells. CONCLUSIONS: We established EV-EphA2-CD9 ELISA for uEV-EphA2 detection for the non-invasive early clinical diagnosis of bladder cancer.

Systematic identification of ALK substrates by integrated phosphoproteome and interactome analysis.

The sensitivity of phosphorylation site identification by mass spectrometry has improved markedly. However, the lack of kinase-substrate relationship (KSR) data hinders the improvement of the range and accuracy of kinase activity prediction. In this study, we aimed to develop a method for acquiring systematic KSR data on anaplastic lymphoma kinase (ALK) using mass spectrometry and to apply this method to the prediction of kinase activity. Thirty-seven ALK substrate candidates, including 34 phosphorylation sites not annotated in the PhosphoSitePlus database, were identified by integrated analysis of the phosphoproteome and crosslinking interactome of HEK 293 cells with doxycycline-induced ALK overexpression. Furthermore, KSRs of ALK were validated by an in vitro kinase assay. Finally, using phosphoproteomic data from ALK mutant cell lines and patient-derived cells treated with ALK inhibitors, we found that the prediction of ALK activity was improved when the KSRs identified in this study were used instead of the public KSR dataset. Our approach is applicable to other kinases, and future identification of KSRs will facilitate more accurate estimations of kinase activity and elucidation of phosphorylation signals.

Identification of CD14 and lipopolysaccharide-binding protein as novel biomarkers for sarcoidosis using proteomics of serum extracellular vesicles.

Sarcoidosis is a complex, polygenic, inflammatory granulomatous multi-organ disease of unknown cause. The granulomatous inflammation in sarcoidosis is driven by the interplay between T cells and macrophages. Extracellular vesicles (EVs) play important roles in intercellular communication. We subjected serum EVs, isolated by size exclusion chromatography, from seven patients with sarcoidosis and five control subjects to non-targeted proteomics analysis. Non-targeted, label-free proteomics analysis detected 2292 proteins in serum EVs; 42 proteins were up-regulated in patients with sarcoidosis relative to control subjects; and 324 proteins were down-regulated. The protein signature of EVs from patients with sarcoidosis reflected disease characteristics such as antigen presentation and immunological disease. Candidate biomarkers were further verified by targeted proteomics analysis (selected reaction monitoring) in 46 patients and 10 control subjects. Notably, CD14 and lipopolysaccharide-binding protein (LBP) were validated by targeted proteomics analysis. Up-regulation of these proteins was further confirmed by immunoblotting, and their expression was strongly increased in macrophages of lung granulomatous lesions. Consistent with these findings, CD14 levels were increased in lipopolysaccharide-stimulated macrophages during multinucleation, concomitant with increased levels of CD14 and LBP in EVs. The area under the curve values of CD14 and LBP were 0.81 and 0.84, respectively, and further increased to 0.98 in combination with angiotensin-converting enzyme and soluble interleukin-2 receptor. These findings suggest that CD14 and LBP in serum EVs, which are associated with granulomatous pathogenesis, can improve the diagnostic accuracy in patients with sarcoidosis.

GSK3 inhibition circumvents and overcomes acquired lorlatinib resistance in ALK-rearranged non-small-cell lung cancer.

Anaplastic lymphoma kinase (ALK) fusion is found in ~3%-5% of patients with non-small-cell lung cancers (NSCLCs). Although the third-generation ALK tyrosine kinase inhibitor (TKI) lorlatinib shows high clinical efficacy in ALK-positive NSCLC, most of the patients eventually relapse with acquired resistance. Recently, drug-tolerant persister (DTP) cells have been considered an important seed of acquired resistance cells. In this study, we established lorlatinib intermediate resistant cells from a patient-derived cell model. Glycogen synthase kinase 3 (GSK3) inhibitions significantly suppressed lorlatinib intermediate resistant cell growth. GSK3 inhibition also sensitized acquired resistance cells derived from alectinib-treated patients with or without secondary mutations to lorlatinib. Therefore, GSK3 plays a crucial role in developing acquired resistance against lorlatinib in ALK-positive NSCLC mediated by lorlatinib intermediate resistant cells and could be a potential molecular target to prevent acquired lorlatinib resistance and overcome ALK-TKI resistance.

Versatile and multiplexed mass spectrometry-based absolute quantification with cell-free-synthesized internal standard peptides.

Preparation of stable isotope-labeled internal standard peptides is crucial for mass spectrometry (MS)-based targeted proteomics. Herein, we developed versatile and multiplexed absolute protein quantification method using MS. A previously developed method based on the cell-free peptide synthesis system, termed MS-based quantification by isotope-labeled cell-free products (MS-QBiC), was improved for multiple peptide synthesis in one-pot reaction. We pluralized the quantification tags used for the quantification of synthesized peptides and thus, made it possible to use cell-free synthesized isotope-labeled peptides as mixtures for the absolute quantification. The improved multiplexed MS-QBiC method was proved to be applied to clarify ribosomal proteins stoichiometry in the ribosomal subunit, one of the largest cellular complexes. The study demonstrates that the developed method enables the preparation of several dozens and even several hundreds of internal standard peptides within a few days for quantification of multiple proteins with only a single-run of MS analysis. SIGNIFICANCE: The developed method can be applied for the preparation of internal standard peptides without limiting the number of peptides to be synthesized, which may result in more practical screening of quantitatively reliable peptides, one of the fundamental steps in the reliable absolute quantification using MS. Furthermore, the method is highly versatile for proteome analysis of any organisms or species without any cDNA or SIL peptide libraries. The quantification can be finished in a few days including design and preparation of appropriate SIL peptides using small-scale batch cell-free reactions, which has a potential to be a part of the standard methodology in a field of quantitative proteomics.

Proteomics of serum extracellular vesicles identifies a novel COPD biomarker, fibulin-3 from elastic fibres.

There is an unmet need for novel biomarkers in the diagnosis of multifactorial COPD. We applied next-generation proteomics to serum extracellular vesicles (EVs) to discover novel COPD biomarkers. EVs from 10 patients with COPD and six healthy controls were analysed by tandem mass tag-based non-targeted proteomics, and those from elastase-treated mouse models of emphysema were also analysed by non-targeted proteomics. For validation, EVs from 23 patients with COPD and 20 healthy controls were validated by targeted proteomics. Using non-targeted proteomics, we identified 406 proteins, 34 of which were significantly upregulated in patients with COPD. Of note, the EV protein signature from patients with COPD reflected inflammation and remodelling. We also identified 63 upregulated candidates from 1956 proteins by analysing EVs isolated from mouse models. Combining human and mouse biomarker candidates, we validated 45 proteins by targeted proteomics, selected reaction monitoring. Notably, levels of fibulin-3, tripeptidyl-peptidase 2, fibulin-1, and soluble scavenger receptor cysteine-rich domain-containing protein were significantly higher in patients with COPD. Moreover, six proteins; fibulin-3, tripeptidyl-peptidase 2, UTP-glucose-1-phosphate uridylyl transferase, CD81, CD177, and oncoprotein-induced transcript 3, were correlated with emphysema. Upregulation of fibulin-3 was confirmed by immunoblotting of EVs and immunohistochemistry in lungs. Strikingly, fibulin-3 knockout mice spontaneously developed emphysema with age, as evidenced by alveolar enlargement and elastin destruction. We discovered potential pathogenic biomarkers for COPD using next-generation proteomics of EVs. This is a novel strategy for biomarker discovery and precision medicine.

Proteomic analysis of urinary and tissue-exudative extracellular vesicles to discover novel bladder cancer biomarkers.

Proteomic analysis of urinary extracellular vesicles (EVs) is a powerful approach to discover potential bladder cancer (BCa) biomarkers, however urine contains numerous EVs derived from the kidney and normal urothelial epithelium, which can obfuscate information related to BCa cell-derived EVs. In this study, we combined proteomic analysis of urinary EVs and tissue-exudative EVs (Te-EVs), which were isolated from culture medium of freshly resected viable BCa tissues. Urinary EVs were isolated from urine samples of 11 individuals (7 BCa patients and 4 healthy individuals), and Te-EVs were isolated from 7 BCa tissues. We performed tandem mass tag (TMT)-labeling liquid chromatography (LC-MS/MS) analysis for both urinary EVs and Te-EVs and identified 1960 proteins in urinary EVs and 1538 proteins in Te-EVs. Most of the proteins identified in Te-EVs were also present in urinary EVs (82.4%), with 55 of these proteins showing upregulated levels in the urine of BCa patients (fold change > 2.0; P < .1). Among them, we selected 22 membrane proteins as BCa biomarker candidates for validation using selected reaction monitoring/multiple reaction monitoring (SRM/MRM) analysis on urine samples from 70 individuals (40 BCa patients and 30 healthy individuals). Six urinary EV proteins (heat-shock protein 90, syndecan-1, myristoylated alanine-rich C-kinase substrate (MARCKS), MARCKS-related protein, tight junction protein ZO-2, and complement decay-accelerating factor) were quantified using SRM/MRM analysis and validated as significantly upregulated in BCa patients (P < .05). In conclusion, the novel strategy that combined proteomic analysis of urinary EVs and Te-EVs enabled selective detection of urinary BCa biomarkers.

Hepatitis C virus enhances Rubicon expression, leading to autophagy inhibition and intracellular innate immune activation.

Autophagy, a degradation system, works to maintain cellular homeostasis. However, as the impact of Hepatitis C virus (HCV) infection on hepatocyte autophagy and its effect on HCV replication remain unclear, we examined them. HCV infection suppressed late-stage autophagy and increased Rubicon. siRNA-mediated knockdown of Rubicon promoted autophagy in HCV-infected cells. In Huh-7 cells harbouring the HCV replicon, Rubicon knockdown downregulated the expression of type 1 interferon (IFN)-related genes and upregulated HCV replication. Rubicon overexpression or administration of bafilomycin A1 or chloroquine, an inhibitor of late-stage autophagy, suppressed autophagy and activated the type 1 IFN pathway. On the other hand, Atg7 knockout suppressed early-stage autophagy and did not activate the type 1 IFN pathway. In livers of humanized liver chimeric mice, HCV infection increased Rubicon and enhanced type 1 IFN signalling. Elimination of HCV in the mice reduced the increase in Rubicon due to HCV infection. The expression levels of Rubicon and IFN-stimulated genes in chronic hepatitis C patients were higher than those in non-B, non-C hepatitis patients. HCV infection increased Rubicon and suppressed hepatocyte autophagy, leading to activation of the intracellular immune response. Rubicon induction is involved in HCV replication via activation of the intracellular immune response.

The role of actinin-4 (ACTN4) in exosomes as a potential novel therapeutic target in castration-resistant prostate cancer.

Prostate cancer is the second leading cause of cancer death in men in the United States. Several novel therapeutic agents have been developed for castration-resistant prostate cancer (CRPC), but the prognosis for patients with CRPC remains poor. The identification of novel therapeutic targets for CRPC is an urgent issue. Exosomes are small vesicles secreted by a variety of cells, and exosomes derived from cancer cells have been reported to circulate in the patient's bodily fluids, promoting metastasis and invasion. We aimed to identify novel therapeutic targets for CRPC by proteomic analysis of serum exosomes. Exosomes were isolated by ultracentrifugation of sera from 36 men with metastatic prostate cancer: untreated (n = 8), well-controlled with primary androgen deprivation therapy (ADT) (n = 8), and CRPC (n = 20). We identified 823 proteins in the serum exosomes. Six proteins were increased in CRPC patients compared with untreated patients. In contrast, only ACTN4 was increased in the CRPC patients compared to the ADT patients. We focused on ACTN4 as a candidate for targeted therapeutics. ACTN4 was highly expressed in the prostate cancer cell line DU145 as well as exosomes from this line. RNA interference-mediated downregulation of ACTN4 significantly attenuated cell proliferation and invasion in DU145 cells. ACTN4 could be a potential therapeutic target for CRPC.

Adenovirus vector-based incorporation of a photo-cross-linkable amino acid into proteins in human primary cells and cancerous cell lines.

The site-specific incorporation of cross-linkable designer amino acids into proteins is useful for covalently bonding protein complexes upon exposure to light. This technology can be used to study networks of protein-protein interactions in living cells; however, to date it has only been applicable for use with a narrow range of cell types, due to the limited availability of plasmid-based transfection protocols. In the present study, we achieved adenovirus-based expression of a variant of an archaeal pyrrolysyl-tRNA synthetase and UAG-recognising tRNA pair, which was used to incorporate unnatural amino acids into proteins at sites defined by in-frame UAG codons within genes. As such, the site-specific photo-cross-linking method is now applicable to a wide variety of mammalian cells. In addition, we repositioned the reactive substituent of a useful photo-cross-linker, Nε-(para-trifluoromethyl-diazirinyl-benzyloxycarbonyl)-l-lysine (pTmdZLys), to the meta position, which improved its availability at low concentration. Finally, we successfully applied this system to analyse the formation of a protein complex in response to a growth signal in human cancerous cells and human umbilical vein endothelial cells. This adenovirus-based system, together with the newly designed cross-linkable amino acid, will facilitate studies on molecular interactions in various cell lines of medical interest.

Targeted Phosphoproteome Analysis Using Selected/Multiple Reaction Monitoring (SRM/MRM).

Mass spectrometry-based phosphoproteomics has been rapidly spread based on the advancement of mass spectrometry and development of efficient enrichment techniques for phosphorylated proteins or peptides. Non-targeted approach has been employed in most of the studies for phosphoproteome analysis. However, targeted approach using selected/multiple reaction monitoring (SRM/MRM) is an indispensible technique used for the quantitation of known targets especially when we have many samples to quantitate phosphorylation events on proteins in biological or clinical samples. We herein describe the application of a large-scale phosphoproteome analysis and SRM-based quantitation for the systematic discovery and validation of biomarkers.

Discovery of colorectal cancer biomarker candidates by membrane proteomic analysis and subsequent verification using selected reaction monitoring (SRM) and tissue microarray (TMA) analysis.

Recent advances in quantitative proteomic technology have enabled the large-scale validation of biomarkers. We here performed a quantitative proteomic analysis of membrane fractions from colorectal cancer tissue to discover biomarker candidates, and then extensively validated the candidate proteins identified. A total of 5566 proteins were identified in six tissue samples, each of which was obtained from polyps and cancer with and without metastasis. GO cellular component analysis predicted that 3087 of these proteins were membrane proteins, whereas TMHMM algorithm predicted that 1567 proteins had a transmembrane domain. Differences were observed in the expression of 159 membrane proteins and 55 extracellular proteins between polyps and cancer without metastasis, while the expression of 32 membrane proteins and 17 extracellular proteins differed between cancer with and without metastasis. A total of 105 of these biomarker candidates were quantitated using selected (or multiple) reaction monitoring (SRM/MRM) with stable synthetic isotope-labeled peptides as an internal control. The results obtained revealed differences in the expression of 69 of these proteins, and this was subsequently verified in an independent set of patient samples (polyps (n = 10), cancer without metastasis (n = 10), cancer with metastasis (n = 10)). Significant differences were observed in the expression of 44 of these proteins, including ITGA5, GPRC5A, PDGFRB, and TFRC, which have already been shown to be overexpressed in colorectal cancer, as well as proteins with unknown function, such as C8orf55. The expression of C8orf55 was also shown to be high not only in colorectal cancer, but also in several cancer tissues using a multicancer tissue microarray, which included 1150 cores from 14 cancer tissues. This is the largest verification study of biomarker candidate membrane proteins to date; our methods for biomarker discovery and subsequent validation using SRM/MRM will contribute to the identification of useful biomarker candidates for various cancers. Data are available via ProteomeXchange with identifier PXD000851.

Establishment of TSH ß real-time monitoring system in mammalian photoperiodism.

Organisms have seasonal physiological changes in response to day length. Long-day stimulation induces thyroid-stimulating hormone beta subunit (TSHß) in the pars tuberalis (PT), which mediates photoperiodic reactions like day-length measurement and physiological adaptation. However, the mechanism of TSHß induction for day-length measurement is largely unknown. To screen candidate upstream molecules of TSHß, which convey light information to the PT, we generated Luciferase knock-in mice, which quantitatively report the dynamics of TSHß expression. We cultured brain slices containing the PT region from adult and neonatal mice and measured the bioluminescence activities from each slice over several days. A decrease in the bioluminescence activities was observed after melatonin treatment in adult and neonatal slices. These observations indicate that the experimental system possesses responsiveness of the TSHß expression to melatonin. Thus, we concluded that our experimental system monitors TSHß expression dynamics in response to external stimuli.

A strategy for large-scale phosphoproteomics and SRM-based validation of human breast cancer tissue samples.

Protein phosphorylation is a key mechanism of cellular signaling pathways and aberrant phosphorylation has been implicated in a number of human diseases. Thus, approaches in phosphoproteomics can contribute to the identification of key biomarkers to assess disease pathogenesis and drug targets. Moreover, careful validation of large-scale phosphoproteome analysis, which is lacking in the current protein-based biomarker discovery, significantly increases the value of identified biomarkers. Here, we performed large-scale differential phosphoproteome analysis using IMAC coupled with the isobaric tag for relative quantification (iTRAQ) technique and subsequent validation by selected/multiple reaction monitoring (SRM/MRM) of human breast cancer tissues in high- and low-risk recurrence groups. We identified 8309 phosphorylation sites on 3401 proteins, of which 3766 phosphopeptides (1927 phosphoproteins) were able to be quantified and 133 phosphopeptides (117 phosphoproteins) were differentially expressed between the two groups. Among them, 19 phosphopeptides were selected for further verification and 15 were successfully quantified by SRM using stable isotope peptides as a reference. The ratio of phosphopeptides between high- and low-risk groups quantified by SRM was well correlated with iTRAQ-based quantification with a few exceptions. These results suggest that large-scale phosphoproteome quantification coupled with SRM-based validation is a powerful tool for biomarker discovery using clinical samples.