Next Article in Journal
Bioelectrical Phase Angle in Patients with Breast Cancer: A Systematic Review
Previous Article in Journal
T-Cells Expressing a Highly Potent PRAME-Specific T-Cell Receptor in Combination with a Chimeric PD1-41BB Co-Stimulatory Receptor Show a Favorable Preclinical Safety Profile and Strong Anti-Tumor Reactivity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Proteogenomics to Urine Analysis towards the Identification of Novel Biomarkers of Prostate Cancer: An Exploratory Study

1
Department of Medical Sciences, Institute of Biomedicine-iBiMED, University of Aveiro, 3810-193 Aveiro, Portugal
2
Cancer Biology and Epigenetics Group, Research Center of Portuguese Oncology Institute of Porto (GEBC CI-IPOP) & Porto Comprehensive Cancer Center (P.CCC), 4200-072 Porto, Portugal
3
UnIC@RISE, Department of Surgery and Physiology, Faculty of Medicine of the University of Porto, 4200-319 Porto, Portugal
4
LAQV/REQUIMTE, Department of Chemistry, University of Aveiro, 3810-193 Aveiro, Portugal
5
Department of Pathology, Portuguese Oncology Institute of Porto (IPO Porto) & Porto Comprehensive Cancer Center (P.CCC), 4200-072 Porto, Portugal
6
Department of Pathology and Molecular Immunology, Institute of Biomedical Sciences Abel Salazar, University of Porto (ICBAS-UP), 4050-513 Porto, Portugal
7
Department of Clinical Pathology, Portuguese Oncology Institute of Porto (IPO Porto) & Porto Comprehensive Cancer Center (P.CCC), 4200-072 Porto, Portugal
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2022, 14(8), 2001; https://doi.org/10.3390/cancers14082001
Submission received: 22 February 2022 / Revised: 12 April 2022 / Accepted: 13 April 2022 / Published: 15 April 2022
(This article belongs to the Topic Prostate Cancer: Symptoms, Diagnosis & Treatment)

Abstract

:

Simple Summary

Prostate cancer (PCa) is one of the most common cancers. Due to the limited and invasive approaches for PCa diagnosis, it is crucial to identify more accurate and non-invasive biomarkers for its detection. The aim of our study was to non-invasively uncover new protein targets for detecting PCa using a proteomics and proteogenomics approach. This work identified several dysregulated mutant protein isoforms in urine from PCa patients, some of them predicted to have a protective or an adverse role in these patients. These results are promising given urine’s non-invasive nature and offers an auspicious opportunity for research and development of PCa biomarkers.

Abstract

To identify new protein targets for PCa detection, first, a shotgun discovery experiment was performed to characterize the urinary proteome of PCa patients. This revealed 18 differentially abundant urinary proteins in PCa patients. Second, selected targets were clinically tested by immunoblot, and the soluble E-cadherin fragment was detected for the first time in the urine of PCa patients. Third, the proteogenome landscape of these PCa patients was characterized, revealing 1665 mutant protein isoforms. Statistical analysis revealed 6 differentially abundant mutant protein isoforms in PCa patients. Analysis of the likely effects of mutations on protein function and PPIs involving the dysregulated mutant protein isoforms suggests a protective role of mutations HSPG2*Q1062H and VASN*R161Q and an adverse role of AMBP*A286G and CD55*S162L in PCa patients. This work originally characterized the urinary proteome, focusing on the proteogenome profile of PCa patients, which is usually overlooked in the analysis of PCa and body fluids. Combined analysis of mass spectrometry data using two different software packages was performed for the first time in the context of PCa, which increased the robustness of the data analysis. The application of proteogenomics to urine proteomic analysis can be very enriching in mutation-related diseases such as cancer.

Graphical Abstract

1. Introduction

Prostate cancer (PCa) is one of the most prevalent cancers among men and the fifth leading cause of cancer-related death [1]. When detected at early stages, PCa can be treated. However, PCa diagnosis is challenging, largely due to the low specificity of PSA tests, particularly in the diagnostic window of 4–10 ng/mL [2], which underscores the need to identify new and more accurate biomarkers.
An ideal biomarker for PCa should be non-invasively assessed, inexpensive, highly sensitive, and specific [3]. For anatomical reasons, urine is enriched in prostatic secretions and better reflects the molecular changes associated with the prostate than blood, which contains markers and confounding factors from the whole body. Urine can be serially collected, requiring minimal processing steps, and presents a simpler matrix with more stability than blood [4].
The phenotype role of proteins combined with the variety of techniques available for proteome analysis makes the search for protein markers in cancer a very attractive strategy [5]. Some promising single-protein biomarkers have been reported, such as AMBP [6] and zinc-alpha-2-glycoprotein (AZGP1) [7,8]. AMBP discriminated PCa and benign prostatic hyperplasia (BPH) patients with a highest accuracy than that estimated for PSA [9], using 2D-DIGE MALDI-TOF/TOF and immunoturbidimetry as discovery and validation approaches, respectively. AZGP1 significantly improved the prediction of PCa in a cohort of candidates for a prostatic biopsy, using isobaric stable isotope labeling and 2D-LC-MS/MS as the discovery method and Western Blot as the validation approach. Multi-marker panels have been shown to improve performance because they better reflect the cancer complexity and heterogeneity, addressing the limitations of single biomarkers. Although promising, no urine protein panel is available for clinical practice due partly to failure in clinical validation, reflecting the need to discover new biomarkers and/or new combinations of biomarkers [7,8]. Interestingly, and to the best of our knowledge, only one assay (Promark®) that quantifies a protein panel in prostate tissue by Mass Spectrometry (MS) is commercially available [10] and, to date, only four mRNA-based urine tests—PCA3 [11], SelectMDX [12], ExoDx Prostate(IntelliScore) [13], and MyProstateScore [14]—have been commercialized.
Cancer is driven by accumulated mutations and other genomic alterations [15]. Mutations on proteins can affect their structure, function, and stability, which may increase their susceptibility to being degraded [16]. As in other types of cancer, in PCa, a weak correlation between RNA and proteins expression is observed. Therefore, the effect of mutations should also be directly investigated at the protein level [17]. To address this inference problem, integration of genome and proteome data (proteogenome) analyses has been performed to identify mutant protein isoforms. Integrated proteogenome analysis can provide new insights into PCa pathophysiology and unveil powerful clinically applicable biomarkers. A shotgun proteomics approach combined with a mutation database has been used to detect mutated peptides related to various types of cancer, such as breast [18], colon [19], and rectal cancer [20]. Still, in PCa, it is mostly unexplored. In 2018, Kwon et al. first applied a proteogenome approach to identify six mutated peptides in the conditioned media from human PCa cell lines related to androgen-independent PCa, which are specific markers for PCa and for metastasis sites [21]. More recently, the same team identified seventy mutant peptides in PCa cell lines, of which seven were differentially expressed in PCa compared to normal tissues [22].
To identify a panel of putative protein markers to be evaluated in a non-invasively collected body fluid for PCa screening, the urine proteome and proteogenome of PCa patients were characterized by an MS-based approach. The integration of results was used to select candidate targets for small-scale clinical testing. MS is widely used to discover urinary protein biomarkers for cancer, including PCa [23]. Usually, biomarker discovery relies on a shotgun proteomics approach, followed by a validation phase using antibody-based techniques or targeted MS. Considering the complex mixture of proteins in urine, separation methodologies are important to increase sensitivity. Thus, a combination of gel-based and gel-free methods, such as GeLC-MS/MS, appears to be a robust and reproducible method for proteome analysis [24], warranting its application in the present work.
This work aims to improve the diagnosis of PCa by investigating the effect of new mutations in proteins that can be detected in urine, a non-invasively collected fluid. Additipnally, it overcomes the limitations of prior studies by using a combination of two software packages for MS data analysis, a proteogenome approach, and a detailed revision and integration of other exploratory proteome analyses to select protein targets.

2. Materials and Methods

2.1. Urine Proteome Profile of PCa Patients and Cancer-Free Subjects

2.1.1. Patients and Sample Collection

Urine samples were collected, without a prior prostate massage, from patients diagnosed with PCa at the Portuguese Oncology Institute of Porto (IPO Porto, Porto, Portugal), before surgery or therapy. Patients with other types of cancer, obesity, or autoimmune diseases were excluded, and cancer-free subjects had no clinically apparent prostatic disease. All available clinical data of the subjects enrolled in this study (discovery (d) and testing cohorts) is depicted in Tables S1 and S2. The discovery cohort comprised five PCa patients and five cancer-free subjects (controls). The testing cohort comprised thirty patients and thirty cancer-free subjects, not considering benign prostate diseases, such as BPH, due to the unavailability of samples.

2.1.2. Urine Sample Preparation

Urine samples were kept at 4 °C and centrifugated at 4000× g for 20 min at 4 °C. The supernatant (4.5 mL per sample) was collected and stored at −80 °C until laboratory analysis. Each urine sample was concentrated using a filter device (10 kDa cut-off, Vivaspin 500 Sartorius Biotech) by sequential centrifugations at 10,000× g for 10 min at 10 °C. Afterward, the retentate was resuspended in 0.5 M Tris pH 6.8 and 4% SDS and protein concentration were assessed by DCTM kit (Bio-Rad, Hercules, CA, USA).

2.1.3. SDS-PAGE

The volume equivalent to 50 µg of protein was precipitated overnight with cold acetone (−20 °C) and centrifugated at 14,000× g for 30 min at 4 °C. Then, the precipitated protein was mixed 1:1 with sample Laemmli loading buffer (0.5 M Tris-HCl pH 6.8, 15% glycerol, 4% SDS, 20% 2-mercaptoethanol, bromophenol blue), heated to 100 °C for 5 min, and separated on 12% Tris-Glycine gels. Following electrophoretic separation, gels were fixed in methanol:acetic acid:water (4:1:5; for 30 min) and stained with Colloidal Coomassie Blue G250 (overnight). Gels were distained with 20% methanol until optimal contrast was achieved.

2.1.4. Liquid Chromatography Tandem-Mass-Spectrometry (LC-MS/MS)

Tryptic digestion was performed according to Shevchenko et al. [25], with a few modifications. All protein bands were manually excised from the gels and sliced into ten sections. The gel pieces were washed with ammonium bicarbonate (NH4HCO3) (25 mM) and ACN (acetonitrile). Proteins were reduced with dithiothreitol (10 mM, 30 min, 60 °C) and alkylated in the dark with iodoacetamide (55 mM, 30 min, 25 °C). The gel pieces were washed with 100 mM NH4HCO3 and then with ACN. Gel pieces were vacuum-dried (SpeedVac, Thermo Savant) and proteins digested with trypsin (Thermo Scientific™, Waltham, MA, USA. Pierce™ Trypsin Protease, MS Grade) in 50 mM NH4HCO3 to a final protease: protein ratio of 1:25 (w/w). After 30 min on ice, 50 μL of 50 mM NH4HCO3 was added, and the samples were incubated for 16 h at 37 °C. The extraction of tryptic peptides was performed by the serial addition of 10% formic acid (FA), 10% FA:ACN (1:1) twice, and 90% ACN. Tryptic peptides were lyophilized and resuspended in 1% FA upon HPLC injection. The samples were analyzed with an Orbitrap Q Exactive (Thermo Fisher Scientific, Bremen, Germany) through the EASY-spray nano ESI source (Thermo Fisher Scientific, Bremen) that was coupled to an Ultimate 3000 (Dionex, Sunnyvale, CA, USA) HPLC system. The trap (5 mm × 300 µm inner diameter) and the EASY-spray analytical (150 mm × 75 µm) columns used were C18 Pepmap100 (Dionex, LC Packings, Sunnyvale, CA, USA), having a particle size of 3 µm. One analytical replicate was performed for each sample and blank runs were acquired between samples. For quality control of the performance of the nano-LC system, the acquisition of cytochrome C digest (1 pmol/μL) (cytochrome c digest lyophilized P/N 161089-thermo scientific) was routinely performed. Peptides were trapped at 30 μL/min in 96% of solvent A (0.1% FA). Elution was achieved with the solvent B (0.1% FA/80% acetonitrile v/v) at 300 nL/min. The 92 min gradient used was as follows: 0–3 min, 4 solvent B; 3–70 min, 4–25% solvent B; 70–90 min, 25–40% solvent B; 90–92 min, 40–90% solvent B; 92–100 min, 90% solvent B; 100–101 min, 90–4% solvent B; 101–120 min, 4% solvent B. The mass spectrometer was operated at 2.2 kV in the data-dependent acquisition mode. An MS2 method was used with an FT survey scan from 400 to 1600 m/z (resolution 70,000; auto-gain control target 1 × 106). The 10 most intense peaks were subjected to high collision dissociation fragmentation (resolution 17,500; auto-gain control target 5 × 104, normalized collision energy 28%, max. injection time 100 ms, dynamic exclusion 35 s).

2.1.5. Protein Identification and Quantification

The MaxQuant (version 1.6.5.0, Thermo software) and Proteome Discoverer (version 2.2, Thermo Fisher Scientific) software packages were used for peptide identification and label-free quantification. In MaxQuant, the Andromeda, and Proteome Discoverer, the MS Amanda, and Sequest HT search engines were used to search the MS/MS spectra against the Uniprot (TrEMBL and Swiss-Prot) protein sequence database under Homo Sapiens (version December 2018). Both database search parameters were as follows: methionine oxidation, protein N-term acetylation and phosphorylation, as variable modifications, and cysteine carbamidomethylation as a fixed modification. The mass tolerance of precursor mass was 20 ppm for MaxQuant and 10 ppm for Proteome Discoverer, and fragment ion mass tolerance was 0.15 Da (MaxQuant) and 0.02 Da (Proteome Discoverer). Minimal peptide length was set to 7 amino acids and, at most, 2 missed cleavages were allowed for both software. The false discovery rate (FDR) for identification was set to 1% at peptide and protein levels. Only the top-ranking protein of each group (master proteins), identified with at least two peptides, were considered. Exclusion of contaminants relied on those identified by the MaxQuant software and the cRAP protein sequences—THE GPM (https://www.thegpm.org/crap/) (accessed on 2 April 2019).
The MS proteome data have been deposited on the ProteomeXchange Consortium via the PRIDE [26] partner repository with the data set identifier PXD017902.

2.1.6. Exploratory Analysis of Urine Proteome Data

The protein abundances in Proteome Discoverer (normalized to the respective median) and normalized LFQ intensities in MaxQuant were log 2-transformed. In an exploratory analysis of proteome data, the proteins identified in all individuals were used as variables to perform Principal Component Analysis (PCA) and Heatmap analyses. These analyses were performed on MetaboAnalyst 5.0 [27]. To identify dysregulated proteins in PCa patients, the fold-change in protein abundance between PCa patients and cancer-free subjects was then calculated from the average log2 difference of protein intensities. Student’s t-test assessed the statistical significance of this difference.

2.1.7. Comparison with a Previous Bioinformatic Analysis of Putative Urinary Markers of PCa and Selection of Candidate Protein Targets for the Testing Phase

Dysregulated proteins were compared with the results of a bioinformatic analysis focused on comparing and mining the proteome profile of tumor prostate tissue and urine from PCa patients reported by several MS studies [28]. The bioinformatic analysis reported 2641 and 616 dysregulated proteins in tumor prostate tissue and urine from PCa patients, respectively. To place urine proteome as a reflection of events taking place in prostate tissue and to identify specific urinary protein targets for PCa, the dysregulated proteins identified in tumor prostate tissue and urine from PCa patients were compared, resulting in 339 overlapping proteins. In this sense, the dysregulated proteins identified by MS in the present work, common to the 2641 dysregulated proteins expressed in tumor prostate tissue or to the 339 urinary proteins with prostate expression, correspond to the selection criteria of candidate proteins to be tested. Then, the selected proteins were compared with the normal human urinary proteome [29].

2.1.8. Measurement of Candidate Protein Targets in Urine Using Immunoblot

The selected protein targets from the discovery phase were tested by slot blot or Western blot immunoassays. In slot blot analysis, performed according to Caseiro et al. [30], the urine protein concentrated fraction was diluted in TBS to a final protein concentration of 0.01 μg/μL and slot-blotted onto a nitrocellulose membrane (Amersham Protran NC 0.45; Amersham Pharmacia Biotech, Buckinghamshire, UK). Antibodies specificity, selectivity, and sensitivity were assessed previously through Western blot by the bands appearing at the expected molecular weights without evidence of non-specific binding of the antibodies. The blocking and incubation conditions were optimized as follows: EFEMP1 (GTX111657: 1:1000, 1 h; GE Amersham-NA934: HRP-linked donkey anti-rabbit 1:10,000); AMBP (sc-81948: 1:1000, 1 h; GE Amersham-NA931: HRP-linked sheep anti-mouse 1:5000); LMAN2 (sc-130026, 1 h; GE Amersham-NA931: HRP-linked sheep anti-mouse 1:5000). Regarding Western blot, 20 µg of protein from each sample was separated on a 12% SDS-PAGE gel and transferred onto nitrocellulose membranes. In both immunoblot experiments, Ponceau S staining was used to normalize the antibody signal to total protein levels. In any case, the membranes were washed with TBS-T (TBS 25 Mm Tris−HCl, pH 7.4, 150 Mm NaCl, 0.1% Tween 20) and imaged in a ChemiDocTM Touch imaging system (Bio-Rad) using the Enhanced Chemiluminescence kit (ECL Select Western Blotting Detection Reagent, RPN2235, Amersham). Optical density was assessed with Image Lab Software (Bio-Rad) and normalized to a loading control sample. Western blot conditions were: CDH1 (GTX629691: 1:1000, 1 h; GE Amersham-NA931: HRP-linked sheep anti-mouse 1:5000); TTR (GTX100577: 1:500, 1 h; GE Amersham-NA934: HRP-linked donkey anti-rabbit 1:10,000).

2.1.9. Measurement of Urinary PSA Levels

Urinary PSA levels were determined using the same method (Elecsys total PSA, 08791732500) used to determine serum PSA levels. This electrochemiluminescence assay is used in the clinical routine of IPO Porto. It quantifies total PSA (free + complexed PSA) using a Cobas e 801 module, a member of Roche Cobas 8000 Modular Analyzer (Roche, Woerden, The Netherlands).

2.2. Urine Proteogenome Profile of PCa Patients and Cancer-Free Subjects

2.2.1. Identification of Cancer-Associated Mutations

Considering the high impact of mutations on cancer progression, the proteogenome profile of urine from PCa patients was explored. For this, mass spectra resulting from the MS analysis were searched against a database built into the Pinnacle software (https://rimuhc.ca/-/protein-quantification-software-pinnacle?redirect=%2Fproteomics-software, accessed on 5 January 2022). This type of analysis aimed to investigate the existence of cancer-associated mutations that were translated in proteins present in the urine from PCa patients. To select high-confidence urinary proteins with a very likely origin in the prostate, only mutations on proteins present in all samples and with known prostate expression were considered. The prostate proteome was searched in the HPA database and in the above-mentioned bioinformatic analysis [28]. The prostate proteome in the HPA consisted of proteins with evidence at the protein level and its last access was on 8 November 2021.

2.2.2. Exploratory Analysis of Urine Proteogenome Data

The abundances of proteins with known prostate expression in Pinnacle were log 2-transformed. In an exploratory analysis of proteogenome data, the levels of mutant protein isoforms identified in all individuals were used as variables to perform Principal Component Analysis (PCA) and Heatmap analyses. These analyses were performed on MetaboAnalyst 5.0 [27]. To identify dysregulated proteins with mutations in PCa patients, the fold-change in protein abundance between PCa patients and cancer-free subjects was then calculated from the average log2 difference of protein intensities. Student’s t-test assessed the statistical significance of this difference.

2.2.3. Integration with the Cancer Genome Atlas (TCGA), DisGeNET and Literature Data

To investigate whether mutations identified in proteins with known prostate expression were already described in PCa, TCGA, DisGeNET (v7.0), and literature data were searched.
TCGA is a cancer genomics consortium that generates data (https://www.cancer.gov/tcga, accessed on 12 January 2022) encompassing the profiling of over 20,000 primary tumors and matched non-tumoral samples related to various human cancers, including PCa. The characterization of PCa samples disclosed 20,237 mutated genes and 33,334 mutations. DisGeNET is one of the largest repositories of Gene-Disease (GDA) and Variant-Disease (VDA) Associations [31]. The latest version of DisGeNET contains 1,134,942 GDAs and 369,554 VDAs. In the present work, variants associated with PCa were extracted from the Prostate Carcinoma C0600139 (January 2022).

2.2.4. Comparison of the Levels of Native and Mutant Forms of Proteins in the Urine from PCa Patients

To investigate the influence of mutations on the abundance of proteins with known expression in the prostate, the levels of their native and mutant forms were compared.

2.2.5. Prediction of the Likely Impact of Single-Residue Substitutions in Proteins

The PolyPhen-2 (Polymorphism Phenotyping v2) web tool was used to predict the likely impact of each amino acid substitution on the structure and function of the proteins with known prostate expression [32]. Each mutation is assigned a score, which is the probability of the substitution being damaging, in addition to a sensitivity and specificity value of the prediction confidence. According to the PolyPhen-2 tool, single-residue substitutions in the protein sequence can be classified as benign (score: 0–0.4), possibly damaging (score: 0.4–0.9), or probably damaging (score: 0.9–1) [33].

2.2.6. Protein–Protein Interaction Analysis

Due to the pivotal role of Protein–Protein interactions (PPIs) in cancer and the possible effect of mutations on its dynamics, the interactions between proteins in which point mutations has been identified were explored. For this, the STRING database v 11.5 was sourced on 12 January 2022, and only protein interactions with a confidence score of ≥0.4 were considered [34]. However, we must be cautious when extrapolating the significance of these PPIs to biological fluids such as urine, as most PPIs are identified or predicted from studies in cells and tissues.

2.2.7. Prediction of the Likely Impact of Single-Residue Substitutions in Protein–Protein Affinity

Considering the impact of mutations on PPIs, the SAAMBE-SEQ Web Server was used to predict the effect of point mutations detected in this work on protein binding affinity [35].

2.3. Statistical Data Analysis

Statistical analyses were carried out in R software for Windows version 3.6.2 and GraphPad Prism version 6.0 (GraphPad Software, Inc.; San Diego, CA, USA). The Shapiro normality test and visual inspection of the histograms were used to assess the data distribution. To evaluate the effect size of the dysregulated proteins when comparing the tested groups, Cohen’s d was determined. Differences were considered statistically significant if p-value was ≤ 0.05. The clinical parameters and protein levels are expressed as mean ± standard deviation (SD).

3. Results

3.1. Urine Proteome Profile of PCa Patients and Cancer-Free Subjects

To identify potential protein targets for PCa prediction, shotgun proteomics was performed in urine collected from PCa patients and cancer-free subjects. To boost MS data analysis, a combination of two different software packages, MaxQuant and Proteome Discover, sourcing three databases (Andromeda, Amanda, and Sequest HT) in total, was used.
Considering only the top-ranking protein of each group identified with at least two peptides and filtering out identifications from reversed sequences and contaminants, 605 and 592 urinary proteins were identified by MaxQuant and Proteome Discoverer, respectively. In total, 732 proteins were identified, excluding those common to both software.

3.1.1. Exploratory Analysis of Urine Proteome Data

Aiming to select and identify proteins of interest for PCa monitoring, only proteins present in all samples analyzed by MaxQuant (82 proteins) and by Proteome Discoverer (84 proteins) were considered for further analysis. These high-confidence proteins were separately used for Principal Component Analysis (PCA) (Figure 1A and Figure 2A) and Heatmap analyses (Figure 1B and Figure 2B). In both software, no separation of groups was observed in the PCa analysis. However, the proteins identified by the MaxQuant software alone seem to provide a discrimination between PCa patients and non-cancer subjects based on two protein clusters, depicted in the heatmap: AZGP1(zinc-alpha-2-glycoprotein)-SPP1 (Osteopontin); CD14 (Monocyte differentiation antigen CD14)-MASP2 (Mannan-binding lectin serine protease 2) (Figure 1B). In the first cluster, proteins are mostly upregulated in PCa patients compared to non-cancer subjects, while in the second cluster proteins are predominantly downregulated in PCa patients.
Then, differential protein analysis revealed 18 dysregulated proteins in PCa, with 4 proteins (p-value ≤ 0.05) identified only by Proteome Discoverer, 9 proteins only by MaxQuant analysis, and 5 proteins (Cadherin-1 (CDH1), EGF-containing fibulin-like extracellular matrix protein 1 (EFEMP1), Prostate-specific antigen (PSA) (KLK3), Secreted and transmembrane protein 1 (SECTM1), and Transthyretin (TTR)) discovered by both software. Altogether, 11 proteins were significantly downregulated (fold change less than 1), and 7 proteins were significantly upregulated (fold change greater than 1) in PCa patients (Table 1 and Table 2). Reassuringly, the most widely used biomarker for PCa diagnosis, PSA, was one of the dysregulated proteins in common in the analysis by both software packages. When the tested groups were compared, proteins showing significant differences (p-value ≤ 0.05) and revealed a “large” effect-size (|Cohen’s d|) > 0.8 (Table 1 and Table 2). Besides a large effect-size, dysregulated proteins identified by both software presented a consistent direction of dysregulation. It is noteworthy that in the heatmap of MaxQuant data, seven proteins (TTR, KLK3, SECTM1, CDH13, AMY2A, EFEMP1, ITIH4, HSPG2, PTGDS, CDH1, and LMAN2) responsible for the separation of groups were also found dysregulated in PCa patients. It was observed that the decreased levels of SECTM1, CDH13, AMY2A, EFEMP1, ITIH4, HSPG2, PTGDS, CDH1, and LMAN2 and increased levels of TTR and KLK3 characterized the urine proteome of PCa patients.

3.1.2. Comparison with a Previous Bioinformatic Analysis of Putative Urinary Markers of PCa and Selection of Candidate Protein Targets for the Testing Phase

To select the most promising proteins for further analysis, dysregulated proteins revealed by MS analysis were compared with proteins resulting from a bioinformatic analysis integrating urine and tumor tissue proteomes of PCa from several MS studies [28]. From this comparison, some common proteins emerged, such as AMBP, CDH1, EFEMP1, KLK3, SECTM1, LMAN2, and TTR.
From the previous study of our group, the dysregulated proteins AMBP, KLK3, LMAN2, and TTR were found dysregulated in urine and tumor tissue from PCa patients, while SECTM1 was only found in urine from PCa patients, and CDH1 and EFEMP1 were only in PCa tissue.
Taken together, and keeping in mind that candidate targets should be urinary proteins with prostate expression, AMBP, CDH1, EFEMP1, KLK3, LMAN2, and TTR were selected for testing in an independent cohort. The presence of these proteins in the urine was already expected, because they are characteristic of the normal human urine proteome [29].

3.1.3. Measurement of Candidate Protein Targets in Urine

Five protein targets, AMBP, CDH1, EFEMP1, LMAN2, and TTR were selected for immunoblot-based testing in a larger and independent cohort (testing group). However, none of the MS findings could be reproduced (Table S3, Figure S1). Measurement of urinary PSA levels in the testing cohort did not agree with the MS findings (p = 0.29, Mann–Whiney test). The results are shown in Figure 3.

3.2. Urine Proteogenome Profile of PCa Patients and Cancer-Free Subjects

3.2.1. Identification of Cancer-Associated Mutations

To characterize the proteogenome landscape of urine from PCa patients, MS/MS spectra were searched against a repository of information from a wide variety of databases encompassing somatic mutations. This search resulted in identifying 6418 mutated peptides corresponding to 1665 mutant protein isoforms. Of these, 609 mutated peptides, which correspond to 417 mutant protein isoforms, were associated with cancer. Only mutant protein isoforms that occurred in all urine samples (322 proteins) were selected for further analysis. Immunoglobulins and highly abundant urinary proteins (serum albumin, uromodulin, serotransferrin) were excluded due to their high abundance in biological samples and the lack of specificity for cancer, resulting in 170 proteins. These 170 proteins corresponded to 122 proteins after filtering out duplicates. As our focus was high confidence proteins with mutations whose origin was very likely the prostate, these data were integrated with the prostate proteome searched in the HPA database and in a bioinformatic analysis [28], resulting in 86 proteins with known expression in the prostate (Table S4). Among these proteins are some of known relevance for PCa, namely Acid ceramidase (ASAH1), Extracellular superoxide dismutase [Cu-Zn] (SOD3), Glutathione S-transferase P (GSTP1), Osteopontin (SPP1), Prostatic acid phosphatase (PAP), and Zinc-alpha-2-glycoprotein (ZAG).

3.2.2. Exploratory Analysis of Urine Proteogenome Data

The levels of the mutant protein isoforms were used for PCA (Principal Component Analysis) (Figure 4A) and Heatmap analyses (Figure 4B). No group separation was observed in the PCA of the proteogenome profile of PCa patients. However, the heatmap indicates a discrimination between PCa patients and non-cancer subjects based on two protein clusters: ITIH4*G893S (Inter-alpha-trypsin inhibitor heavy chain H4)-LMAN2*D222N (Vesicular integral-membrane protein VIP36); KLK3*C209Y (PSA)-MVB12B*T198M (Multivesicular body subunit 12B) (Figure 4B). In the first cluster, mutant forms of proteins are mostly downregulated in PCa patients compared to non-cancer subjects, while in the second cluster mutant forms of proteins are upregulated predominantly in PCa patients.

3.2.3. Integration with the Cancer Genome Atlas (TCGA), DisGeNET and Literature Data

According to TCGA, DisGeNET, and the literature, only three of the mutations identified in the 86 proteins with known prostate expression have already been described. These mutations (rs17632542, rs1695, rs7041) were mapped on KLK3 (PSA) [36], GSTP1 (Glutathione S-transferase P) [37,38], and GC (Vitamin D-binding protein) [39], respectively. To the best of our knowledge, there is no association of the remaining mutant protein isoforms with PCa. Especially notable are the proteins SPP1, VASN, ASAH1, RBP4, and ASS1, which, until now, have had no mutation related to PCa described in the literature.

3.2.4. Comparison of the Levels of Native and Mutant Forms of Proteins in the Urine from PCa Patients

The analysis of proteogenome data revealed 6 differentially abundant mutant protein isoforms in PCa patients compared with cancer-free individuals, namely Protein AMBP (AMBP*A286G), Sodium/hydrogen exchanger 9B1 (SLC9B1*N70S), Basement membrane-specific heparan sulfate proteoglycan core protein (HSPG2*Q1062H), Zinc finger protein 624 (ZNF624*S207F), Vasorin (VASN*R161Q), and Complement decay-accelerating factor (CD55*S162L) (Table S4, Figure S2). Mutant AMBP isoform was upregulated in PCa patients, while the remaining 5 differentially abundant mutant protein isoforms were downregulated.
Comparing the proteome profile analysis of MaxQuant and Proteome Discoverer with the proteogenome profile of PCa patients resulted in 30 and 31 common proteins, respectively. Of these common proteins, AMBP, CDH1, EFEMP1, HSPG2, ITIH4, KLK3, LMAN2, PTGDS, VASN, and CD55 proteins stood out. The native form of AMBP, CDH1, EFEMP1, HSPG2, ITIH4, KLK3, LMAN2, and PTGDS proteins was found dysregulated in urine from PCa patients, but only the mutant protein isoforms (AMBP*A286G; HSPG2*Q1062H) were found dysregulated (Figure S2). In the remaining common proteins, the presence of mutations did not affect their abundance in urine. The native form of VASN and CD55 proteins was not found dysregulated in the urine from PCa patients, but their mutant protein isoforms (VASN*R161Q; CD55*S162L) were.
The mutations identified in these proteins and in those with recognized relevance to PCa are summarized in Table 3.

3.2.5. Prediction of the Likely Impact of Single-Residue Substitutions in Proteins

With the purpose of determining the potential impact of point mutations on protein function, PolyPhen-2 tool was used. It is worthy of mention that AMBP*A286G and CD55*S162L mutant protein isoforms were predicted to be probably damaging, while SLC9B1*N70S, ZNF624*S207F, VASN*R161Q, and HSPG2*Q1062H were predicted to be benign. Most point mutations were predicted to be possibly or probably damaging. The results are presented in Table 4 and Table S5.

3.2.6. Protein–Protein Interaction Analysis

In addition to impacting the function of proteins, mutations can also affect interactions between proteins and, consequently, important biological processes and signaling pathways. To predict interactions between the proteins in which point mutations were identified, the STRING search tool was used. As shown in Figure 5, the network consisted of 86 connected proteins (nodes) through 214 edges with different confidence levels. The protein–protein interaction enrichment p-value was <1.0 × 10−16. Reactome enrichment analysis showed 12 pathways enriched in this network (Table S6). Regulation of Insulin-like Growth Factor (IGF) transport and uptake by Insulin-like Growth Factor Binding Proteins (IGFBPs) was the third most important pathway in this network, while Extracellular matrix (ECM) organization was the tenth. This network shows predicted interactions between most of the proteins.

3.2.7. Prediction of the Likely Impact of Single-Residue Substitutions in Protein–Protein Affinity

To predict the impact of point mutations on PPIs, the SAAMBE-SEQ tool was used. The likely effect of AMBP*A286G, HSPG2*Q1062H, VASN*R161Q, and CD55*S162L point mutations on protein–protein interactions was scrutinized. Point mutations detected on SLC9B1 and ZNF624 were not examined as these proteins do not interact with any proteins in the network. Additionally, the impact of point mutations on proteins involved in the IGF pathway was also explored. This analysis revealed that the likely effect of these point mutations is destabilizing for PPIs (Table S7).

4. Discussion

The limitations and the invasive nature of serum PCa screening have driven the discovery of new candidate urinary biomarkers, especially protein markers. However, so far, none has translated into clinically useful tools, reflecting the need to discover novel biomarkers and/or new combinations of biomarkers. Thus, this study aimed to take advantage of a non-invasively collected biofluid, urine, and a high throughput approach, proteomics, to identify new protein targets for predicting the risk of developing PCa. This work was divided into three stages: characterization of the urine proteome profile and selection of protein targets; testing of shortlisted protein targets in a larger, independent cohort; and characterization of the urine proteogenome profile. The urine proteome profile of PCa and cancer-free subjects was analyzed by two software packages and 18 dysregulated proteins, of which 5 (TTR, EFEMP1, CDH1, SECTM1, KLK3) common to both software, were found. The integration of the urine proteome profile of PCa patients with proteome data from other studies reviewed by us [28] supported the selection of potential discriminatory protein targets. As a result, AMBP, CDH1, EFEMP1, LMAN2, and TTR stood out as potential targets and were tested in an independent cohort of patients. In this testing phase, incubation with anti-E-cadherin did not result in a band around 120 kDa (full-length protein), but rather a band about 80 kDa. We realized that this 80 kDa fragment corresponded to soluble E-cadherin (sE-cadherin) and has been previously identified in tissue and serum from PCa patients [93,94] and in urine from patients with other cancers [95,96], using antibody-based techniques. Concerning PCa, as far as we know, here we present the first report of the detection of sE-cadherin fragment in the urine. Kuefer et al. [93] suggested that the 80 kDa fragment is originated from the extracellular domain of full-length E-cadherin. Increased levels of sE-cadherin have been reported in serum and tumor prostate tissue from PCa patients and are correlated with disease stage [94,97,98]. Differential abundances of these MS-detected proteins were tested in an independent cohort using immunoblot, but different variations were observed. Additionally, urinary PSA levels were also assessed in this independent cohort, but did not distinguish PCa patients from controls, which agrees with other studies [99].
The proteogenome landscape of urine from PCa patients was then characterized and 1665 mutant protein isoforms were disclosed, of which 417 were cancer-related mutations. After considering only mutations present in all urine samples and proteins with known prostate expression, 86 mutant protein isoforms emerged. Among these proteins are some of known relevance for PCa, namely Acid ceramidase (ASAH1), Extracellular superoxide dismutase [Cu-Zn] (SOD3), Glutathione S-transferase P (GSTP1), Osteopontin (SPP1), Prostatic Acid Phosphatase (PAP), and Zinc-Alpha-2-Glycoprotein (ZAG). PAP is gaining renewed interest due to its superior predictive role of cause-specific survival and GS compared to serum PSA in men with high risk PCa [100,101]. Remarkably, it was recently suggested that a form of PAP (PLPAcP) associates with early PCa [102]. Identifying a new mutation in this protein in a non-invasive biological fluid, adding to the prediction of PAP mutation to be probably damaging, strengthens the renewed interest in its study in PCa. Mutations found on the 86 proteins were searched for in databases and the literature and, to the best of our knowledge, only rs17632542 [36,103,104,105], rs1695 [37,38,106,107], and rs7041 [39] mutations mapped on PSA, GSTP1, and GC proteins have been described in the PCa context. In that vein, these results validate the proteogenome analysis performed in the present study.
The analysis of the urine proteogenome profile of PCa patients revealed 6 differentially abundant mutant protein isoforms, namely AMBP*A286G, SLC9B1*N70S, HSPG2*Q1062H, ZNF624*S207F, VASN*R161Q, and CD55*S162L. From the comparison of the proteome and proteogenome profile of PCa patients, AMBP, CDH1, EFEMP1, KLK3, and LMAN2 proteins stood out. Their native form was found dysregulated in urine from PCa patients, but the same was not observed with their mutant form, with the exception of AMBP*A286G and HSPG2*Q1062H. These results may explain the differences between MS and immunoblot data, because the antibodies either do not recognize the mutated peptides or do not specifically recognize them.
PPIs play a pivotal role in most biological processes. Dysregulation of these protein interactions may result in pathological conditions, such as cancer, being involved in tumor progression, invasion, and metastasis [108,109]. In this sense, PPIs have been claimed as promising therapeutic targets for numerous types of cancer, including for PCa. For this type of cancer, 28 small molecules and 14 peptides have been proposed to disrupt PPIs with relevance to PCa progression [110]. To explore PPIs between proteins with known prostate expression and the pathways in which these interactions were involved, the STRING tool was used. In this analysis, the IGF transport and uptake by IGFBPs proved to be the third most important pathway in the network. The IGF axis is a network of ligands (GF1, IGF2, insulin) and IGFBP receptors (IGF1R, IGF2R, INSR), the latter being responsible for mediating the activity of IGFs [111]. IGFs are oncogenic regulators, promoting prostate tumor growth, survival, and proliferation, and the role of IGF axis has been well documented in PCa. For instance, IGFBP-2 enhanced proliferation of androgen-independent prostate cancer cells [112] and IGF-I levels were found raised in serum and prostate tissue from PCa patients, being a predictor of risk for this type of cancer [113,114]. In accordance with this, IGF1R and INSR act as oncogenes in PCa, enhancing tumor growth, proliferation, invasion, and angiogenesis [115]. Considering the relevance of the IGF pathway in PCa, the impact of mutations on the interaction of proteins involved in this pathway was predicted. According to SAAMBE-SEQ, the mutations were predicted to destabilize all PPIs involved in the IGF pathway, which naturally could affect this pathway and consequently the progression of PCa.
To investigate the likely impact of each amino acid substitution on protein function and PPIs involving the dysregulated mutant protein isoforms (AMBP*A286G, SLC9B1*N70S, HSPG2*Q1062H, ZNF624*S207F, VASN*R161Q, and CD55*S162L), the PolyPhen-2 and SAAMBE-SEQ prediction tools were used. The role of the SLC9B1 and ZNF624 proteins on cancer is completely unknown, so the downregulation of their mutant protein isoforms and the prediction of their benign impact do not allow conclusions to be drawn. HSPG2, in its intact form, is a well-described pro-angiogenic molecule, being correlated with GS and increased cell proliferation and viability [55,56,116]. The intact form of this protein was found increased in tumor prostate tissue, but in sera from PCa patients raised levels of HSPG2-derived fragments resulting from matrix metalloproteinase 7 (MMP7) degradation were observed. These fragments were mostly originated from domain IV and were not present in sera from non-cancer subjects, suggesting that HSPG2 cleavage occurs during metastasis and before the protein enters the bloodstream. Using an in silico analysis, Grindel et al. predicted that domains III and V of HSPG2 are the most prone to cleavage by MMP-7 and generate new peptides for other extracellular proteases to digest [55]. Curiously, in this work, the mutated peptide identified in the mutant HSPG2 isoform is located on domain III. The cleavage of HSPG2 and other components of basement membrane occurs during PCa cell invasion and is orchestrated by proteases such as MMPs, cathepsin L, and BMP1/Tolloid-like proteases. Both Cathepsin L and BMP1/Tolloid-like proteases cleave HSPG2 in domain V, originating the Endorepellin [117] and LG3 [118] peptides, respectively. Unlike the intact form, cleaved Endorepellin and LG3 peptides behave as powerful anti-angiogenic factors, being claimed as potential therapeutic targets for cancer [118]. In fact, the administration of endorepellin to mice with squamous cell carcinomas and lung carcinomas resulted in mitigation of tumor growth, angiogenesis and metabolism and promotion of tumor hypoxia [119]. Accordingly, LG3-diminished levels were noticed in breast cancer cells and in plasma from breast cancer patients [120]. Only the LG3 peptide has been detected in urine [121,122]. In PCa, both the existence and the role of these peptides are unknown, and the only recognized HSPG2 protease is MMP7. A complex network between HSPG2 and other basement membrane components, such as collagens, laminin, and nidogen is responsible for ECM integrity. When this integrity is disturbed, the metastatic process is compromised [123]. In the present work, mutations were identified in HSPG2, collagens, nidogen, and in other proteins involved in ECM organization. When the impact of these mutations on PPIs was predicted, they all proved to be destabilizing, which eventually affects ECM dynamics and tumor progression. All these results, together with the fact that the HSPG2*Q1062H point mutation was predicted to be benign and the mutant peptide was downregulated in PCa patients, suggest that this mutant peptide may have beneficial effects in patients with PCa and opens doors for its study in PCa treatment. Concerning the AMBP protein, it is cleaved into three chains, namely Alpha-1-microglobulin, Bikunin, and Trypstatin. The function of the AMBP protein in cancer remains undisclosed. However, it has been claimed that the AMBP-derived product bikunin is underexpressed in oral squamous cell carcinoma and plays an antitumor role [40]. In line with this, there is evidence that bikunin significantly prevented tumor invasion and metastasis in Lewis lung carcinoma and ovarian carcinoma cells [124,125]. Curiously, in this work, the mutant peptide identified in the AMBP isoform is located on the bikunin fragment. The mutation identified in AMBP was predicted to be probably damaging, destabilized all PPIs in which AMBP was involved, and resulted in an upregulation of mutant AMBP isoform in PCa patients. This may suggest a detrimental role of this mutation on PCa patients. Regarding CD55, it blocks complement response by accelerating the decay of C3 and C5 convertases [126] and is involved in PCa cell survival and metastasis [92]. This interplay between CD55 and C3 is visible by their interaction in the STRING network. The mutation detected on the CD55 protein was predicted to be probably damaging and destabilizing for CD55-C3 interaction. With these findings, it seems reasonable to suspect the detrimental role of this mutation on PCa patients. Regarding VASN, it is a known inhibitor of TGF-β signaling [127]. The TGF-β pathway has a dual role in cancer, because it prevents cell proliferation in early stages and in advanced stages stimulates proliferation, epithelial-to-mesenchymal transition (EMT) and evasion of immune surveillance, and attenuates apoptosis [128]. The mechanism involved in this inhibitory action of VASN on TGF-beta was revealed in breast cancer cell lines. It was demonstrated that a soluble form of VASN resulting from the proteolytic shedding of its extracellular domain by Metalloprotease domain 17 (ADAM17) is responsible for controlling the TGFβ pathway [129]. In PCa, the role of VASN is largely unexplored, including the interplay between the VASN and TGFβ pathways. However, overexpression of VASN in prostate tumor tissue and in serum from PCa patients and the subsequent promotion of cell proliferation and PCa progression have already been reported, in agreement with other types of cancer [90]. Interestingly, in this work, the mutated peptide identified in the VASN protein is located on the extracellular domain of the protein, the domain cleaved by ADAM17. The mutation identified in VASN resulted in a downregulation of this mutant protein isoform in PCa patients and was predicted to be benign, which may suggest a protective role of this mutation on PCa patients.
These findings indicate that, in mutational diseases such as cancer and in biofluids with high proteolytic activity, such as urine, the application of proteogenomics to urine analysis and the study of peptides can be very enriching because point mutations can go unnoticed at the protein level but are detected at the peptide level. This may sharpen or renew interest in underexplored targets, as observed in this work. We hope to address some of these questions in future work. Furthermore, it would be interesting to test these mutant peptides by an MS-targeted approach such as MRM, but this is beyond the scope of this work. This work’s novelty lies in the proteogenome characterization of urine from PCa patients and the combined analysis of MS data using two different software packages, increasing certainty in the identification of urinary proteins modulated by PCa.

5. Conclusions

The majority of mutations identified in this work have never been associated with PCa, and some are predicted to be damaging, which offers an auspicious opportunity for research and development of PCa biomarkers, especially in the HSPG2 context. Additionally, the discovery of cancer-associated mutations in PCa-related proteins in urine is promising given this biofluid’s non-invasive and dynamic nature.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14082001/s1: Table S1, Clinical data of subjects included in the discovery cohort group; Table S2, Clinical data of subjects included in the testing cohort group; Table S3, Summary of statistical analysis results of shortlisted proteins evaluated in the testing group; Table S4, List of mutant protein isoforms identified in the 86 proteins with known prostate expression; Table S5, Prediction of likely impact of point mutations on protein function using PolyPhen-2 tool; Table S6, Reactome pathway enrichment analysis of the network; Table S7, Prediction of likely impact of point mutations on protein–protein interactions using SAAMBE-SEQ tool; Figure S1, Original Western blots figures; Figure S2, Levels of AMBP*A286G, SLC9B1*N70S, HSPG2*Q1062H, ZNF624*S207F, VASN*R161Q, and CD55*S162L mutant protein isoforms and respective levels of native form (when applicable) in the urine from PCa patients.

Author Contributions

Conceptualization, T.L., R.H., R.V. and M.F.; methodology, T.L., R.H., R.V. and M.F.; validation, T.L., R.H., R.V. and M.F.; formal analysis, T.L., A.S.B., F.T., R.F., C.J., A.L.-M., R.H., R.V. and M.F.; investigation, T.L.; resources, D.B.-S., C.J., L.A., A.L.-M., R.H., R.V. and M.F.; data curation, A.S.B., F.T. and R.F.; writing—original draft preparation, T.L.; writing—review and Editing, T.L., A.S.B., F.T., R.F., A.L.-M., D.B.-S., L.A., C.J., R.H., R.V. and M.F.; visualization, T.L., R.H., R.V. and M.F.; supervision, R.H., R.V. and M.F.; project administration, T.L.; funding acquisition, T.L., R.H., R.V. and M.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Portuguese Foundation for Science and Technology (FCT), European Union, QREN, FEDER, and COMPETE for the Unidade de Investigação Cardiovascular (UIDB/IC/00051/2020 and UIDP/00051/2020), Institute of Biomedicine (iBiMED) (UIDB/04501/2020, POCI-01-0145-FEDER-007628), and FCT QOPNA ((FCT UID/QUI/00062/2019) and LAQV/REQUIMTE (UIDB/50006/2020) research units and the DOCnet (NORTE-01-0145-FEDER-000003), by Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF). T.L. is supported by an individual scholarship (SFRH/BD/136904/2018), F.T. by a post-doctoral research grant by UnIC (UIDP/00051/2020), and R.V. by individual fellowship grant (IF/00286/2015).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the IPO-Porto Ethics Committee (Comissão de Ética para a Saúde, Reference 282R/2017).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data generated during Mass Spectrometry analysis is available in the ProteomeXchange Consortium via the PRIDE partner repository with the data set identifier PXD017902.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  2. Saini, S. PSA and beyond: Alternative prostate cancer biomarkers. Cell. Oncol. 2016, 39, 97–106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Jakobsen, N.A.; Hamdy, F.C.; Bryant, R.J. Novel biomarkers for the detection of prostate cancer. J. Clin. Urol. 2016, 9, 3–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Eskra, J.N.; Daniel, R.; Pavlovich, C.P.; Catalona, W.J.; Luo, J. Approaches to urinary detection of prostate cancer. Prostate Cancer Prostatic Dis. 2019, 22, 362–381. [Google Scholar] [CrossRef]
  5. Van der Burgt, Y.E.M. Protein biomarker discovery is still relevant and has entered a new phase. eBioMedicine 2019, 43, 15. [Google Scholar] [CrossRef] [Green Version]
  6. Davalieva, K.; Kiprijanovska, S.; Komina, S.; Petrusevska, G.; Zografska, C.C.; Polenakovic, M. Proteomics analysis of urine reveals acute phase response proteins as candidate diagnostic biomarkers for prostate cancer. Proteome Sci. 2015, 13, 2. [Google Scholar] [CrossRef] [Green Version]
  7. Garbis, S.D.; Tyritzis, S.I.; Roumeliotis, T.; Zerefos, P.; Giannopoulou, E.G.; Vlahou, A.; Kossida, S.; Diaz, J.; Vourekas, S.; Tamvakopoulos, C.; et al. Search for potential markers for prostate cancer diagnosis, prognosis and treatment in clinical tissue specimens using amine-specific isobaric tagging (iTRAQ) with two-dimensional liquid chromatography and tandem mass spectrometry. J. Proteome Res. 2008, 7, 3146–3158. [Google Scholar] [CrossRef]
  8. Katafigiotis, I.; Tyritzis, S.I.; Stravodimos, K.G.; Alamanis, C.; Pavlakis, K.; Vlahou, A.; Makridakis, M.; Katafigioti, A.; Garbis, S.D.; Constantinides, C.A. Zinc α2-glycoprotein as a potential novel urine biomarker for the early diagnosis of prostate cancer. BJU Int. 2012, 110, E688–E693. [Google Scholar] [CrossRef]
  9. Louie, K.S.; Seigneurin, A.; Cathcart, P.; Sasieni, P. Do prostate cancer risk models improve the predictive accuracy of PSA screening? A meta-analysis. Ann. Oncol. 2015, 26, 848–864. [Google Scholar] [CrossRef]
  10. Blume-Jensen, P.; Berman, D.M.; Rimm, D.L.; Shipitsin, M.; Putzi, M.; Nifong, T.P.; Small, C.; Choudhury, S.; Capela, T.; Coupal, L.; et al. Biology of human tumors development and clinical validation of an in situ biopsy-based multimarker assay for risk stratification in prostate cancer. Clin. Cancer Res. 2015, 17, 372–381. [Google Scholar] [CrossRef] [Green Version]
  11. Deras, I.L.; Aubin, S.M.J.; Blase, A.; Day, J.R.; Koo, S.; Partin, A.W.; Ellis, W.J.; Marks, L.S.; Fradet, Y.; Rittenhouse, H.; et al. PCA3: A molecular urine assay for predicting prostate biopsy outcome. J. Urol. 2008, 179, 1587–1592. [Google Scholar] [CrossRef] [PubMed]
  12. Van Neste, L.; Hendriks, R.J.; Dijkstra, S.; Trooskens, G.; Cornel, E.B.; Jannink, S.A.; de Jong, H.; Hessels, D.; Smit, F.P.; Melchers, W.J.G.; et al. Detection of high-grade prostate cancer using a urinary molecular biomarker—Based risk score. Eur. Urol. 2016, 70, 740–748. [Google Scholar] [CrossRef] [PubMed]
  13. McKiernan, J.; Donovan, M.J.; O’Neill, V.; Bentink, S.; Noerholm, M.; Belzer, S.; Skog, J.; Kattan, M.W.; Partin, A.; Andriole, G.; et al. A novel urine exosome gene expression assay to predict high-grade prostate cancer at initial biopsy. JAMA Oncol. 2016, 2, 882–889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Tosoian, J.J.; Dunn, R.L.; Niknafs, Y.S.; Saha, A.; Vince, R.A.; St. Sauver, J.L.; Jacobson, D.J.; McGree, M.E.; Siddiqui, J.; Groskopf, J.; et al. Association of urinary myprostatescore, age, and prostate volume in a longitudinal cohort of healthy men: Long-term findings from the olmsted county study. Eur. Urol. Open Sci. 2021, 29, 30–35. [Google Scholar] [CrossRef]
  15. Lin, T.; Zhang, T.; Kitata, R.B.; Liu, T.; Smith, R.D.; Qian, W.; Shi, T. Mass spectrometry-based targeted proteomics for analysis of protein mutations. Mass Spectrom. Rev. 2021, e21741. [Google Scholar] [CrossRef]
  16. Reva, B.; Antipin, Y.; Sander, C. Predicting the functional impact of protein mutations: Application to cancer genomics. Nucleic Acids Res. 2011, 39, e118. [Google Scholar] [CrossRef] [Green Version]
  17. Sinha, A.; Huang, V.; Livingstone, J.; Wang, J.; Fox, N.S.; Kurganovs, N.; Ignatchenko, V.; Fritsch, K.; Donmez, N.; Heisler, L.E.; et al. The proteogenomic landscape of curable prostate cancer. Cancer Cell 2019, 35, 414–427.e6. [Google Scholar] [CrossRef] [Green Version]
  18. Lazar, I.M.; Karcini, A.; Ahuja, S.; Estrada-Palma, C. Proteogenomic analysis of protein sequence alterations in breast cancer cells. Sci. Rep. 2019, 9, 10381. [Google Scholar] [CrossRef] [Green Version]
  19. Mathivanan, S.; Ji, H.; Tauro, B.J.; Chen, Y.S.; Simpson, R.J. Identifying mutated proteins secreted by colon cancer cell lines using mass spectrometry. J. Proteom. 2012, 76, 141–149. [Google Scholar] [CrossRef]
  20. Zhang, B.; Wang, J.; Wang, X.; Zhu, J.; Liu, Q.; Shi, Z.; Chambers, M.C.; Zimmerman, L.J.; Shaddox, K.F.; Kim, S.; et al. Proteogenomic characterization of human colon and rectal cancer. Nature 2014, 513, 382–387. [Google Scholar] [CrossRef] [Green Version]
  21. Kwon, O.K.; Jeon, J.M.I.; Sung, E.; Na, A.Y.; Kim, S.J.; Lee, S. Comparative secretome profiling and mutant protein identification in metastatic prostate cancer cells by quantitative mass spectrometry-based proteomics. Cancer Genom. Proteom. 2018, 15, 279–290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Kwon, O.K.; Ha, Y.-S.; Lee, J.N.; Kim, S.; Lee, H.; Chun, S.Y.; Kwon, T.G.; Lee, S. Comparative proteome profiling and mutant protein identification in metastatic prostate cancer cells by quantitative mass spectrometry-based proteogenomics. Cancer Genom. Proteom. 2019, 16, 273–286. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Yang, J.; Roy, R.; Jedinak, A.; Moses, M.A. Mining the human urinary proteome biomarker discovery for human cancer and its metastases. Cancer J. 2015, 21, 327–336. [Google Scholar] [CrossRef] [PubMed]
  24. Dzieciatkowska, M.; Hill, R.C.; Hansen, K. GeLC-MS/MS analysis of complex protein mixtures. Methods Mol. Biol. 2014, 1156, 53–66. [Google Scholar] [CrossRef]
  25. Shevchenko, A.; Tomas, H.; Havliš, J.; Olsen, J.V.; Mann, M. In-gel digestion for mass spectrometric characterization of proteins and proteomes. Nat. Protoc. 2007, 1, 2856–2860. [Google Scholar] [CrossRef]
  26. Perez-Riverol, Y.; Csordas, A.; Bai, J.; Bernal-Llinares, M.; Hewapathirana, S.; Kundu, D.J.; Inuganti, A.; Griss, J.; Mayer, G.; Eisenacher, M.; et al. The PRIDE database and related tools and resources in 2019: Improving support for quantification data. Nucleic Acids Res. 2019, 47, D442–D450. [Google Scholar] [CrossRef]
  27. Pang, Z.; Chong, J.; Zhou, G.; de Lima Morais, D.A.; Chang, L.; Barrette, M.; Gauthier, C.; Jacques, P.É.; Li, S.; Xia, J. MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 2021, 49, W388–W396. [Google Scholar] [CrossRef]
  28. Lima, T.; Henrique, R.; Vitorino, R.; Fardilha, M. Bioinformatic analysis of dysregulated proteins in prostate cancer patients reveals putative urinary biomarkers and key biological pathways. Med. Oncol. 2021, 38, 10381. [Google Scholar] [CrossRef]
  29. Zhao, M.; Li, M.; Yang, Y.; Guo, Z.; Sun, Y.; Shao, C.; Li, M.; Sun, W.; Gao, Y. A comprehensive analysis and annotation of human normal urinary proteome. Sci. Rep. 2017, 7, 3024. [Google Scholar] [CrossRef]
  30. Caseiro, A.; Barros, A.; Ferreira, R.; Padrão, A.; Aroso, M.; Quintaneiro, C.; Pereira, A.; Marinheiro, R.; Vitorino, R.; Amado, F. Pursuing type 1 diabetes mellitus and related complications through urinary proteomics. Transl. Res. 2014, 163, 188–199. [Google Scholar] [CrossRef]
  31. Piñero, J.; Ramírez-Anguita, J.M.; Saüch-Pitarch, J.; Ronzano, F.; Centeno, E.; Sanz, F.; Furlong, L.I. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020, 48, D845–D855. [Google Scholar] [CrossRef] [Green Version]
  32. Adzhubei, I.A.; Schmidt, S.; Peshkin, L.; Ramensky, V.E.; Gerasimova, A.; Bork, P.; Kondrashov, A.S.; Sunyaev, S.R. A method and server for predicting damaging missense mutations. Nat. Methods 2010, 7, 248–249. [Google Scholar] [CrossRef] [Green Version]
  33. Adzhubei, I.; Jordan, D.M.; Sunyaev, S.R. Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Protoc. Hum. Genet. 2013, 76, 7–20. [Google Scholar] [CrossRef] [Green Version]
  34. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef]
  35. Li, G.; Pahari, S.; Murthy, A.K.; Liang, S.; Fragoza, R.; Yu, H.; Alexov, E. SAAMBE-SEQ: A sequence-based method for predicting mutation effect on protein-protein binding affinity. Bioinformatics 2021, 37, 992–999. [Google Scholar] [CrossRef]
  36. Otto, J.J.; Correll, V.L.; Engstroem, H.A.; Hitefield, N.L.; Main, B.P.; Albracht, B.; Johnson-Pais, T.; Yang, L.F.; Liss, M.; Boutros, P.C.; et al. Targeted mass spectrometry of a clinically relevant psa variant from post-DRE urines for quantitation and genotype determination. Proteom. Clin. Appl. 2020, 14, e2000012. [Google Scholar] [CrossRef]
  37. Santric, V.; Djokic, M.; Suvakov, S.; Pljesa-Ercegovac, M.; Nikitovic, M.; Radic, T.; Acimovic, M.; Stankovic, V.; Bumbasirevic, U.; Milojevic, B.; et al. GSTP1 rs1138272 polymorphism affects prostate cancer risk. Medicina 2020, 56, 128. [Google Scholar] [CrossRef] [Green Version]
  38. Zhang, Y.; Yuan, Y.; Chen, Y.; Wang, Z.; Li, F.; Zhao, Q. Association between GSTP1 Ile105Val polymorphism and urinary system cancer risk: Evidence from 51 studies. Onco. Targets. Ther. 2016, 9, 3565–3569. [Google Scholar] [CrossRef] [Green Version]
  39. Gilbert, R.; Bonilla, C.; Metcalfe, C.; Lewis, S.; Evans, D.M.; Fraser, W.D.; Kemp, J.P.; Donovan, J.L.; Hamdy, F.C.; Neal, D.E.; et al. Associations of vitamin D pathway genes with circulating 25-hydroxyvitamin-D, 1, 25-dihydroxyvitamin-D, and prostate cancer: A nested case–control study. Cancer Causes Control. 2015, 26, 205–218. [Google Scholar] [CrossRef] [Green Version]
  40. Sekikawa, S.; Onda, T.; Miura, N.; Nomura, T.; Takano, N.; Shibahara, T.; Honda, K. Underexpression of α-1-microglobulin/bikunin precursor predicts a poor prognosis in oral squamous cell carcinoma. Int. J. Oncol. 2018, 53, 2605–2614. [Google Scholar] [CrossRef] [Green Version]
  41. Huang, H.; Han, Y.; Gao, J.; Feng, J.; Zhu, L.; Qu, L.; Shen, L.; Shou, C. High level of serum AMBP is associated with poor response to paclitaxel-capecitabine chemotherapy in advanced gastric cancer patients. Med. Oncol. 2013, 30, 748. [Google Scholar] [CrossRef]
  42. Fujita, K.; Kume, H.; Matsuzaki, K.; Kawashima, A.; Ujike, T. Proteomic analysis of urinary extracellular vesicles from high Gleason score prostate cancer. Sci. Rep. 2017, 7, 42961. [Google Scholar] [CrossRef] [Green Version]
  43. Hamm, A.; Veeck, J.; Bektas, N.; Wild, P.J.; Hartmann, A.; Heindrichs, U.; Kristiansen, G.; Werbowetski-Ogilvie, T.; Del Maestro, R.; Knuechel, R.; et al. Frequent expression loss of Inter-alpha-trypsin inhibitor heavy chain (ITIH) genes in multiple human solid tumors: A systematic expression analysis. BMC Cancer 2008, 8, 25. [Google Scholar] [CrossRef] [Green Version]
  44. Davalieva, K.; Kiprijanovska, S.; Kostovska, I.M.; Stavridis, S.; Stankov, O.; Komina, S.; Petrusevska, G.; Polenakovic, M. Comparative proteomics analysis of urine reveals down-regulation of acute phase response signaling and LXR/RXR activation pathways in prostate cancer. Proteomes 2017, 6, 1. [Google Scholar] [CrossRef] [Green Version]
  45. Zhao, H.; Zhao, X.; Lei, T.; Zhang, M. Screening, identification of prostate cancer urinary biomarkers and verification of important spots. Invest. New Drugs 2019, 37, 935–947. [Google Scholar] [CrossRef]
  46. Rodríguez-Blanco, G.; Zeneyedpour, L.; Duijvesz, D.; Marije Hoogland, A.; Verhoef, E.I.; Kweldam, C.F.; Burgers, P.C.; Smitt, P.S.; Bangma, C.H.; Jenster, G.; et al. Tissue proteomics outlines AGR2 AND LOX5 as markers for biochemical recurrence of prostate cancer. Oncotarget 2018, 9, 36444–36456. [Google Scholar] [CrossRef]
  47. Quanico, J.; Franck, J.; Gimeno, J.P.; Sabbagh, R.; Salzet, M.; Day, R.; Fournier, I. Parafilm-assisted microdissection: A sampling method for mass spectrometry-based identification of differentially expressed prostate cancer protein biomarkers. Chem. Commun. 2015, 51, 4564–4567. [Google Scholar] [CrossRef]
  48. Khan, A.P.; Poisson, L.M.; Bhat, V.B.; Fermin, D.; Zhao, R.; Kalyana-Sundaram, S.; Michailidis, G.; Nesvizhskii, A.I.; Omenn, G.S.; Chinnaiyan, A.M.; et al. Quantitative proteomic profiling of prostate cancer reveals a role for miR-128 in prostate cancer. Mol. Cell. Proteom. 2010, 9, 298–312. [Google Scholar] [CrossRef] [Green Version]
  49. Debelec-Butuner, B.; Alapinar, C.; Ertunc, N.; Gonen-Korkmaz, C.; Yörükoğlu, K.; Korkmaz, K.S. TNFα-mediated loss of β-catenin/E-cadherin association and subsequent increase in cell migration is partially restored by NKX3.1 expression in prostate cells. PLoS ONE 2014, 9, e109868. [Google Scholar] [CrossRef] [Green Version]
  50. Tsui, K.-H.; Lin, Y.-H.; Chung, L.-C.; Chuang, S.-T.; Feng, T.-H.; Chiang, K.-C.; Chang, P.-L.; Yehg, C.-J.; Juang, H.H. Prostate-derived ets factor represses tumorigenesis and modulatesepithelial-to-mesenchymal transition in bladder carcinoma cells. Cancer Lett. 2016, 375, 142–151. [Google Scholar] [CrossRef]
  51. Li, F.; Pascal, L.E.; Stolz, D.B.; Wang, K.; Zhou, Y.; Chen, W.; Xu, Y.; Chen, Y.; Dhir, R.; Parwani, A.V.; et al. E-cadherin is down-regulated in benign prostate hyperplasia and required for tight junction formation and permeability barrier in prostatic epithelial cell monolayer. Prostate 2020, 79, 1226–1237. [Google Scholar] [CrossRef] [PubMed]
  52. Shen, H.; Zhang, L.; Zhou, J.; Chen, Z.; Yang, G.; Liao, Y.; Zhu, M. Epidermal growth factor-containing fibulin-like extracellular matrix protein 1 (EFEMP1) acts as a potential diagnostic biomarker for prostate cancer. Med. Sci. Monit. 2017, 23, 216–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Kim, Y.J.; Yoon, H.Y.; Kim, S.K.; Kim, Y.W.; Kim, E.J.; Kim, I.Y.; Kim, W.J. EFEMP1 as a novel DNA methylation marker for prostate cancer: Array-based DNA methylation and expression profiling. Clin. Cancer Res. 2011, 17, 4523–4530. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Almeida, M.; Costa, V.L.; Costa, N.R.; Ramalho-Carvalho, J.; Baptista, T.; Ribeiro, F.R.; Paulo, P.; Teixeira, M.R.; Oliveira, J.; Lothe, R.A.; et al. Epigenetic regulation of EFEMP1 in prostate cancer: Biological relevance and clinical potential. J. Cell. Mol. Med. 2014, 18, 2287–2297. [Google Scholar] [CrossRef] [PubMed]
  55. Grindel, B.; Li, Q.; Arnold, R.; Petros, J.; Zayzafoon, M.; Muldoon, M.; Stave, J.; Chung, L.W.K.; Farach-Carson, M.C. Perlecan/HSPG2 and matrilysin/MMP-7 as indices of tissue invasion: Tissue localization and circulating perlecan fragments in a cohort of 288 radical prostatectomy patients. Oncotarget 2016, 7, 10433–10447. [Google Scholar] [CrossRef] [Green Version]
  56. Savorè, C.; Zhang, C.; Muir, C.; Liu, R.; Wyrwa, J.; Shu, J.; Zhau, H.E.; Chung, L.W.K.; Carson, D.D.; Farach-Carson, M.C. Perlecan knockdown in metastatic prostate cancer cells reduces heparin-binding growth factor responses in vitro and tumor growth in vivo. Clin. Exp. Metastasis 2005, 22, 377–390. [Google Scholar] [CrossRef]
  57. Whitelock, J.; Melrose, J.; Iozzo, R.V. Diverse cell signaling events modulated by perlecan. Biochemistry 2008, 47, 11174–11183. [Google Scholar] [CrossRef] [Green Version]
  58. Piñeiro, M.; Andrés, M.; Iturralde, M.; Carmona, S.; Hirvonen, J.; Pyörälä, S.; Heegaard, P.M.H.; Tjørnehøj, K.; Lampreave, F.; Piñeiro, A.; et al. ITIH4 (inter-alpha-trypsin inhibitor heavy chain 4) is a new acute-phase protein isolated from cattle during experimental infection. Infect. Immun. 2004, 72, 3777–3782. [Google Scholar] [CrossRef] [Green Version]
  59. Huang, M.; Zhang, W.; Zhao, B.; Li, L. Relationship between inter-α-trypsin inhibitor heavy chain 4 and ovarian cancer. Chin. J. Cancer Res. 2019, 31, 955–964. [Google Scholar] [CrossRef]
  60. Duffy, M.J. Biomarkers for prostate cancer: Prostate-specific antigen and beyond. Clin. Chem. Lab. Med. 2020, 58, 326–339. [Google Scholar] [CrossRef] [Green Version]
  61. Söllner, T. Adaptor protein CD2AP and L-type lectin LMAN2 regulate exosome cargo protein trafficking through the golgi. J. Biol. Chem. 2016, 291, 25462–25475. [Google Scholar] [CrossRef] [Green Version]
  62. The L-Type Lectin LMAN2 Is Over-Expressed in Brain Metastatic Breast Cancer. Available online: https://www.sciencegate.app/app/document/full-text#10.31219/osf.io/s3k2n. (accessed on 1 January 2022).
  63. Potapenko, I.O.; Haakensen, V.D.; Lüders, T.; Helland, Å.; Bukholm, I.; Sørlie, T.; Kristensen, V.N.; Lingjærde, O.C.; Børresen-Dale, A.L. Glycan gene expression signatures in normal and malignant breast tissue; possible role in diagnosis and progression. Mol. Oncol. 2010, 4, 98–118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. L’Espérance, S.; Popa, I.; Bachvarova, M.; Plante, M.; Patten, N.; Wu, L.; Têtu, B.; Bachvarov, D. Gene expression profiling of paired ovarian tumors obtained prior to and following adjuvant chemotherapy: Molecular signatures of chemoresistant tumors. Int. J. Oncol. 2006, 29, 5–24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Thompson, V.C.; Day, T.K.; Bianco-Miotto, T.; Selth, L.A.; Han, G.; Thomas, M.; Buchanan, G.; Scher, H.I.; Nelson, C.C.; Greenberg, N.M.; et al. A gene signature identified using a mouse model of androgen receptor-dependent prostate cancer predicts biochemical relapse in human disease. Int. J. Cancer 2012, 131, 662–672. [Google Scholar] [CrossRef] [PubMed]
  66. Kim, J.; Yang, P.; Suraokar, M.; Sabichi, A.L.; Llansa, N.D.; Mendoza, G.; Subbarayan, V.; Logothetis, C.J.; Newman, R.A.; Lippman, S.M.; et al. Suppression of prostate tumor cell growth by stromal cell prostaglandin D synthase-derived products. Cancer Res. 2005, 65, 6189–6198. [Google Scholar] [CrossRef] [Green Version]
  67. Ogretmen, B.; Hannun, Y.A. Biologically active sphingolipids in cancer pathogenesis and treatment. Nat. Rev. Cancer 2004, 4, 604–616. [Google Scholar] [CrossRef]
  68. Seelan, R.S.; Qian, C.; Yokomizo, A.; Bostwick, D.G.; Smith, D.I.; Liu, W. Human acid ceramidase is overexpressed but not mutated in prostate cancer. Genes Chromosom. Cancer 2000, 29, 137–146. [Google Scholar] [CrossRef]
  69. Camacho, L.; Meca-Cortés, Ó.; Abad, J.L.; García, S.; Rubio, N.; Díaz, A.; Celiá-Terrassa, T.; Cingolani, F.; Bermudo, R.; Fernández, P.L.; et al. Acid ceramidase as a therapeutic target in metastatic prostate cancer. J. Lipid Res. 2013, 54, 1207–1220. [Google Scholar] [CrossRef] [Green Version]
  70. Saad, A.F.; Meacham, W.D.; Bai, A.; Anelli, V.; Elojeimy, S.; Mahdy, A.E.M.; Turner, L.S.; Cheng, J.; Bielawska, A.; Bielawski, J.; et al. The functional effects of acid ceramidase overexpression in prostate cancer progression and resistance to chemotherapy. Cancer Biol. Ther. 2007, 6, 1455–1460. [Google Scholar] [CrossRef] [Green Version]
  71. McCord, J.M.; Edeas, M.A. SOD, oxidative stress and human pathologies: A brief history and a future vision. Biomed. Pharmacother. 2005, 59, 139–142. [Google Scholar] [CrossRef]
  72. Zheng, S.; Lin, X.; Gan, X.; Wang, X. The impact of SOD3 on prostatic diseases: Elevated SOD3 is a novel biomarker for the diagnosis of chronic nonbacterial prostatiti. Res. Sq. 2021. [Google Scholar] [CrossRef]
  73. Kim, J.; Mizokami, A.; Shin, M.; Izumi, K.; Konaka, H.; Kadono, Y.; Kitagawa, Y.; Keller, E.T.; Zhang, J.; Namiki, M. SOD3 acts as a tumor suppressor in PC-3 prostate cancer cells via hydrogen peroxide accumulation. Anticancer Res. 2014, 34, 2821–2832. [Google Scholar]
  74. Allocati, N.; Masulli, M.; Di Ilio, C.; Federici, L. Glutathione transferases: Substrates, inihibitors and pro-drugs in cancer and neurodegenerative diseases. Oncogenesis 2018, 7, 8. [Google Scholar] [CrossRef]
  75. Lin, X.; Tascilar, M.; Lee, W.H.; Vles, W.J.; Lee, B.H.; Veeraswamy, R.; Asgari, K.; Freije, D.; Van Rees, B.; Gage, W.R.; et al. GSTP1 CpG island hypermethylation is responsible for the absence of GSTP1 expression in human prostate cancer cells. Am. J. Pathol. 2001, 159, 1815–1826. [Google Scholar] [CrossRef] [Green Version]
  76. Nakayama, M.; Bennett, C.J.; Hicks, J.L.; Epstein, J.I.; Platz, E.A.; Nelson, W.G.; De Marzo, A.M. Hypermethylation of the human glutathione S-transferase-π gene (GSTP1) CpG island is present in a subset of proliferative inflammatory atrophy lesions but not in normal or hyperplastic epithelium of the prostate: A detailed study using laser-capture micro. Am. J. Pathol. 2003, 163, 923–933. [Google Scholar] [CrossRef]
  77. Brooks, J.D.; Weinstein, M.; Lin, X.; Sun, Y.; Pin, S.S.; Bova, G.S.; Epstein, J.I.; Isaacs, W.B.; Nelson, W.G. CG island methylation changes near the GSTP1 gene in prostatic intraepithelial neoplasia. Cancer Epidemiol. Biomark. Prev. 1998, 7, 531–536. [Google Scholar]
  78. Badowska-Kozakiewicz, A.M.; Budzik, M.P. The multidirectional role of osteopontin in cancer. Nowotwory 2018, 68, 176–183. [Google Scholar] [CrossRef]
  79. Khodavirdi, A.C.; Song, Z.; Yang, S.; Zhong, C.; Wang, S.; Wu, H.; Pritchard, C.; Nelson, P.S.; Roy-Burman, P. Increased expression of osteopontin contributes to the progression of prostate cancer. Cancer Res. 2006, 66, 883–888. [Google Scholar] [CrossRef] [Green Version]
  80. Yu, A.; Guo, K.; Qin, Q.; Xing, C.; Zu, X. Clinicopathological and prognostic significance of osteopontin expression in patients with prostate cancer: A systematic review and meta-analysis. Biosci. Rep. 2021, 41, BSR20203531. [Google Scholar] [CrossRef]
  81. Wiśniewski, T.; Zyromska, A.; Makarewicz, R.; Zekanowska, E. Osteopontin and angiogenic factors as new biomarkers of prostate cancer. Urol. J. 2019, 16, 134–140. [Google Scholar] [CrossRef]
  82. Prager, A.J.; Peng, C.R.; Lita, E.; Mcnally, D.; Kaushal, A.; Sproull, M.; Compton, K.; Dahut, W.L.; Figg, W.D.; Citrin, D.; et al. Urinary aHGF, IGFBP3 and OPN as diagnostic and prognostic biomarkers for prostate cancer. Biomark. Med. 2013, 7, 831–841. [Google Scholar] [CrossRef]
  83. Graddis, T.J.; McMahan, C.J.; Tamman, J.; Page, K.J.; Trager, J.B. Prostatic acid phosphatase expression in human tissues. Int. J. Clin. Exp. Pathol. 2011, 4, 295–306. [Google Scholar] [CrossRef] [Green Version]
  84. Gutman, B.Y.A.B.; Gutman, E.B. An “acid” phosphatase occurring in the serum of patients with metastasizing carcinoma of the prostate gland. J. Clin. Investig. 1938, 17, 473–478. [Google Scholar] [CrossRef]
  85. Bing, C.; Bao, Y.; Jenkins, J.; Sanders, P.; Manieri, M.; Cinti, S.; Tisdale, M.J.; Trayhurn, P. Zinc-α2-glycoprotein, a lipid mobilizing factor, is expressed in adipocytes and is up-regulated in mice with cancer cachexia. Proc. Natl. Acad. Sci. USA 2004, 101, 2500–2505. [Google Scholar] [CrossRef] [Green Version]
  86. Bondar, O.P.; Barnidge, D.R.; Klee, E.W.; Davis, B.J.; Klee, G.G. LC-MS/MS quantification of Zn-α2 glycoprotein: A potential serum biomarker for prostate cancer. Clin. Chem. 2007, 53, 673–678. [Google Scholar] [CrossRef] [Green Version]
  87. Hale, L.P.; Price, D.T.; Sanchez, L.M.; Demark-Wahnefried, W.; Madden, J.F. Zinc α-2-glycoprotein is expressed by malignant prostatic epithelium and may serve as a potential serum marker for prostate cancer. Clin. Cancer Res. 2001, 7, 846–853. [Google Scholar]
  88. Chintapalli, V.R.; Kato, A.; Henderson, L.; Hirata, T.; Woods, D.J.; Overend, G.; Davies, S.A.; Romero, M.F.; Dow, J.A.T. Transport proteins NHA1 and NHA2 are essential for survival, but have distinct transport modalities. Proc. Natl. Acad. Sci. USA 2015, 112, 11720–11725. [Google Scholar] [CrossRef] [Green Version]
  89. Shi, W.; Dong, F.; Jiang, Y.; Lu, L.; Wang, C.; Tan, J.; Yang, W.; Guo, H.; Ming, J.; Huang, T. Construction of prognostic microRNA signature for human invasive breast cancer by integrated analysis. Onco. Targets. Ther. 2019, 12, 1979–2010. [Google Scholar] [CrossRef] [Green Version]
  90. Cui, F.L.; Mahmud, A.N.; Xu, Z.P.; Wang, Z.Y.; Hu, J.P. VASN promotes proliferation of prostate cancer through the YAP/TAZ axis. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 6589–6596. [Google Scholar] [CrossRef]
  91. Dho, S.H.; Lim, J.C.; Kim, L.K. Beyond the role of CD55 as a complement component. Immune Netw. 2018, 18, e11. [Google Scholar] [CrossRef] [Green Version]
  92. Loberg, R.D.; Day, L.S.L.; Dunn, R.; Kalikin, L.M.; Pienta, K.J. Inhibition of decay-accelerating factor (CD55) attenuates prostate cancer growth and survival in vivo. Neoplasia 2006, 8, 69–78. [Google Scholar] [CrossRef] [Green Version]
  93. Kuefer, R.; Hofer, M.D.; Gschwend, J.E.; Pienta, K.J.; Sanda, M.G.; Chinnaiyan, A.M.; Rubin, M.A.; Day, M.L. The role of an 80 kDa fragment of E-cadherin in the metastatic progression of prostate cancer. Clin. Cancer Res. 2003, 9, 6447–6452. [Google Scholar]
  94. Iacopino, F.; Pinto, F.; Bertaccini, A.; Calarco, A.; Proietti, G.; Totaro, A.; Martorana, G.; Bassi, P.; Sica, G. Soluble E-cadherin and IL-6 serum levels in patients affected by prostate cancer before and after prostatectomy. Oncol. Rep. 2012, 28, 370–374. [Google Scholar] [CrossRef] [Green Version]
  95. Banks, R.E.; Porter, W.H.; Whelan, P.; Smith, P.H.; Selby, P.J. Soluble forms of the adhesion molecule E-cadherin in urine. J. Clin. Pathol. 1995, 48, 179–180. [Google Scholar] [CrossRef] [Green Version]
  96. Katayama, M.; Hirai, S.; Yasumoto, M.; Nishikawa, K.; Nagata, S.; Otsuka, M.; Kato, I.; Kamihagi, K. Soluble fragments of E-cadherin cell adhesion molecule increase in urinary excretion of cancer patients, potentially indicating its shedding from epithelial tumor cells. Int. J. Oncol. 1994, 5, 1049–1057. [Google Scholar] [CrossRef]
  97. Tsaur, I.; Thurn, K.; Juengel, E.; Gust, K.M.; Borgmann, H.; Mager, R.; Bartsch, G.; Oppermann, E.; Ackermann, H.; Nelson, K.; et al. sE-cadherin serves as a diagnostic and predictive parameter in prostate cancer patients. J. Exp. Clin. Cancer Res. 2015, 34, 43. [Google Scholar] [CrossRef] [Green Version]
  98. Kuefer, R.; Hofer, M.D.; Zorn, C.S.M.; Engel, O.; Volkmer, B.G.; Juarez-Brito, M.A.; Eggel, M.; Gschwend, J.E.; Rubin, M.A.; Day, M.L. Assessment of a fragment of e-cadherin as a serum biomarker with predictive value for prostate cancer. Br. J. Cancer 2005, 92, 2018–2023. [Google Scholar] [CrossRef] [Green Version]
  99. Duijvesz, D.; Versluis, C.Y.L.; Van Der Fels, C.A.M.; Vredenbregt-Van Den Berg, M.S.; Leivo, J.; Peltola, M.T.; Bangma, C.H.; Pettersson, K.S.I.; Jenster, G. Immuno-based detection of extracellular vesicles in urine as diagnostic marker for prostate cancer. Int. J. Cancer 2015, 137, 2869–2878. [Google Scholar] [CrossRef]
  100. Xu, H.; Wang, F.; Li, H.; Ji, J.; Cao, Z.; Lyu, J.; Shi, X.; Zhu, Y.; Zhang, C.; Guo, F. Prostatic acid phosphatase (PAP) predicts prostate cancer progress in a population-based study: The renewal of PAP? Dis. Markers. 2019, 2019, 7090545. [Google Scholar] [CrossRef]
  101. Fang, L.C.; Dattoli, M.; Taira, A.; True, L.; Sorace, R.; Wallner, K. Prostatic Acid Phosphatase Adversely Affects Cause-Specific Survival in Patients with Intermediate to High-Risk Prostate Cancer Treated with Brachytherapy. Urology 2008, 71, 146–150. [Google Scholar] [CrossRef]
  102. Alpert, E.; Akhavan, A.; Gruzman, A.; Hansen, W.J.; Lehrer-Graiwer, J.; Hall, S.C.; Johansen, E.; McAllister, S.; Gulati, M.; Lin, M.-F.; et al. Multifunctionality of prostatic acid phosphatase in prostate cancer pathogenesis. Biosci. Rep. 2021, 41, BSR20211646. [Google Scholar] [CrossRef]
  103. Kote-Jarai, Z.; Amin Al Olama, A.; Leongamornlert, D.; Tymrakiewicz, M.; Saunders, E.; Guy, M.; Giles, G.G.; Severi, G.; Southey, M.; Hopper, J.L.; et al. Identification of a novel prostate cancer susceptibility variant in the KLK3 gene transcript. Hum. Genet. 2011, 129, 687–694. [Google Scholar] [CrossRef] [Green Version]
  104. Parikh, H.; Wang, Z.; Pettigrew, K.A.; Jia, J.; Daugherty, S.; Yeager, M.; Jacobs, K.B.; Hutchinson, A.; Burdett, L.; Cullen, M.; et al. Fine mapping the KLK3 locus on chromosome 19q13.33 associated with prostate cancer susceptibility and PSA levels. Hum. Genet. 2011, 129, 675–685. [Google Scholar] [CrossRef] [Green Version]
  105. Gudmundsson, J.; Besenbacher, S.; Sulem, P.; Gudbjartsson, D.F.; Olafsson, I.; Arinbjarnarson, S.; Agnarsson, B.A.; Benediktsdottir, K.R.; Isaksson, H.J.; Kostic, J.P.; et al. Genetic correction of PSA values using sequence variants associated with PSA levels. Sci. Transl. Med. 2010, 2, 62ra92. [Google Scholar] [CrossRef] [Green Version]
  106. Cotignola, J.; Leonardi, D.B.; Shahabi, A.; Acuña, A.D.; Stern, M.C.; Navone, N.; Scorticati, C.; De Siervi, A.; Mazza, O.; Vazquez, E. Glutathione-S-transferase (GST) polymorphisms are associated with relapse after radical prostatectomy. Prostate Cancer Prostatic Dis. 2013, 16, 28–34. [Google Scholar] [CrossRef]
  107. Oskina, N.A.; Ermolenko, N.A.; Boyarskih, U.A.; Lazarev, A.F.; Petrova, V.D.; Ganov, D.I.; Tonacheva, O.G.; Lifschitz, G.I.; Filipenko, M.L. Associations between SNPs within antioxidant genes and the risk of prostate cancer in the Siberian region of Russia. Pathol. Oncol. Res. 2014, 20, 635–640. [Google Scholar] [CrossRef]
  108. Gulfidan, G.; Turanli, B.; Beklen, H.; Sinha, R.; Arga, K.Y. Pan-cancer mapping of differential protein-protein interactions. Sci. Rep. 2020, 10, 3272. [Google Scholar] [CrossRef]
  109. Engin, B.H.; Guney, E.; Keskin, O.; Oliva, B.; Gursoy, A. Integrating structure to protein-protein interaction networks that drive metastasis to brain and lung in breast cancer. PLoS ONE 2013, 8, e81035. [Google Scholar] [CrossRef] [Green Version]
  110. Matos, B.; Howl, J.; Jerónimo, C.; Fardilha, M. The disruption of protein-protein interactions as a therapeutic strategy for prostate cancer. Pharmacol. Res. 2020, 161, 105145. [Google Scholar] [CrossRef]
  111. Foulstone, E.; Prince, S.; Zaccheo, O.; Burns, J.L.; Harper, J.; Jacobs, C.; Church, D.; Hassan, A.B. Insulin-like growth factor ligands, receptors, and binding proteins in cancer. J. Pathol. 2005, 205, 145–153. [Google Scholar] [CrossRef]
  112. Chatterjee, S.; Sung Park, E.; Soloff, M.S. Proliferation of DU145 prostate cancer cells is inhibited by suppressing insulin-like growth factor binding protein-2. Int. J. Urol. 2004, 11, 876–884. [Google Scholar] [CrossRef]
  113. Hellawell, G.O.; Turner, G.D.H.; Davies, D.R.; Poulsom, R.; Brewster, S.F.; Macaulay, V.M. Expression of the type 1 insulin-like growth factor receptor is up-regulated in primary prostate cancer and commonly persists in metastatic disease. Cancer Res. 2002, 62, 2942–2950. [Google Scholar]
  114. Wolk, A.; Mantzoros, C.S.; Andersson, S.O.; Bergström, R.; Signorello, L.B.; Lagiou, P.; Adami, H.O.; Trichopoulos, D. Insulin-like growth factor 1 and prostate cancer risk: A population- based, case-control study. J. Natl. Cancer Inst. 1998, 90, 911–915. [Google Scholar] [CrossRef] [Green Version]
  115. Heidegger, I.; Kern, J.; Ofer, P.; Klocker, H.; Massoner, P. Oncogenic functions of IGF1R and INSR in prostate cancer include enhanced tumor growth, cell migration and angiogenesis. Oncotarget 2014, 5, 2723–2735. [Google Scholar] [CrossRef] [Green Version]
  116. Datta, M.W.; Hernandez, A.M.; Schlicht, M.J.; Kahler, A.J.; DeGueme, A.M.; Dhir, R.; Shah, R.B.; Farach-Carson, C.; Barrett, A.; Datta, S. Perlecan, a candidate gene for the CAPB locus, regulates prostate cancer cell growth via the Sonic Hedgehog pathway. Mol. Cancer 2006, 5, 9. [Google Scholar] [CrossRef] [Green Version]
  117. Mongiat, M.; Sweeney, S.M.; San Antonio, J.D.; Fu, J.; Iozzo, R.V. Endorepellin, a novel inhibitor of angiogenesis derived from the C terminus of perlecan. J. Biol. Chem. 2003, 278, 4238–4249. [Google Scholar] [CrossRef] [Green Version]
  118. Gonzalez, E.M.; Reed, C.C.; Bix, G.; Fu, J.; Zhang, Y.; Gopalakrishnan, B.; Greenspan, D.S.; Iozzo, R.V. BMP-1/Tolloid-like metalloproteases process endorepellin, the angiostatic C-terminal fragment of perlecan. J. Biol. Chem. 2005, 280, 7080–7087. [Google Scholar] [CrossRef] [Green Version]
  119. Bix, G.; Castello, R.; Burrows, M.; Zoeller, J.J.; Weech, M.; Iozzo, R.A.; Cardi, C.; Thakur, M.L.; Barker, C.A.; Camphausen, K.; et al. Endorepellin in vivo: Targeting the tumor vasculature and retarding cancer growth and metabolism. J. Natl. Cancer Inst. 2006, 98, 1634–1646. [Google Scholar] [CrossRef]
  120. Chang, J.W.; Kang, U.B.; Kim, D.H.; Yi, J.K.; Lee, J.W.; Noh, D.Y.; Lee, C.; Yu, M.H. Identification of circulating endorepellin LG3 fragment: Potential use as a serological biomarker for breast cancer. Proteom.-Clin. Appl. 2008, 2, 23–32. [Google Scholar] [CrossRef]
  121. Parker, T.J.; Sampson, D.L.; Broszczak, D.; Chng, Y.L.; Carter, S.L.; Leavesley, D.I.; Parker, A.W.; Upton, Z. A fragment of the LG3 peptide of endorepellin is present in the urine of physically active mining workers: A potential marker of physical activity. PLoS ONE 2012, 7, e33714. [Google Scholar] [CrossRef] [Green Version]
  122. Oda, O.; Shinzato, T.; Ohbayashi, K.; Takai, I.; Kunimatsu, M.; Maeda, K.; Yamanaka, N. Purification and characterization of perlecan fragment in urine of end-stage renal failure patients. Clin. Chim. Acta 1996, 255, 119–132. [Google Scholar] [CrossRef]
  123. Venning, F.A.; Wullkopf, L.; Erler, J.T. Targeting ECM disrupts cancer progression. Front. Oncol. 2015, 5, 224. [Google Scholar] [CrossRef] [Green Version]
  124. Kobayashi, H.; Shinohara, H.; Fujie, M.; Gotoh, J.; Itoh, M.; Takeuchi, K.; Terao, T. Inhibition of metastasis of lewis lung carcinoma by urinary trypsin inhibitor in experimental and spontaneous metastasis models. Int. J. Cancer 1995, 63, 455–462. [Google Scholar] [CrossRef]
  125. Suzuki, M.; Kobayashi, H.; Tanaka, Y.; Hirashima, Y.; Kanayama, N.; Takei, Y.; Saga, Y.; Suzuki, M.; Itoh, H.; Terao, T. Suppression of invasion and peritoneal carcinomatosis of ovarian cancer cell line by overexpression of bikunin. Int. J. Cancer 2003, 104, 289–302. [Google Scholar] [CrossRef] [Green Version]
  126. Geller, A.; Yan, J. The role of membrane bound complement regulatory proteins in tumor development and cancer immunotherapy. Front. Immunol. 2019, 10, 1074. [Google Scholar] [CrossRef] [Green Version]
  127. Ikeda, Y.; Imai, Y.; Kumagai, H.; Nosaka, T.; Morikawa, Y.; Hisaoka, T.; Manabe, I.; Maemura, K.; Nakaoka, T.; Imamura, T.; et al. Vasorin, a transforming growth factor β-binding protein expressed in vascular smooth muscle cells, modulates the arterial response to injury in vivo. Proc. Natl. Acad. Sci. USA 2004, 101, 10732–10737. [Google Scholar] [CrossRef] [Green Version]
  128. Zhang, M.; Zhang, Y.Y.; Chen, Y.; Wang, J.; Wang, Q.; Lu, H. TGF-β Signaling and Resistance to Cancer Therapy. Front. Cell Dev. Biol. 2021, 9, 786728. [Google Scholar] [CrossRef]
  129. Malapeira, J.; Esselens, C.; Bech-Serra, J.J.; Canals, F.; Arribas, J. ADAM17 (TACE) regulates TGFΒ signaling through the cleavage of vasorin. Oncogene 2011, 30, 1912–1922. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Exploratory analysis of proteome data from MaxQuant. (A) Principal Component Analysis of the urine proteome of the two groups. (B) The heatmap of proteins identified in all individuals. Samples are represented in columns and proteins in rows. Proteins whose gene name is not available are indicated by their UniProt accession number. The dashed line on the heatmap indicates the two clusters of proteins.
Figure 1. Exploratory analysis of proteome data from MaxQuant. (A) Principal Component Analysis of the urine proteome of the two groups. (B) The heatmap of proteins identified in all individuals. Samples are represented in columns and proteins in rows. Proteins whose gene name is not available are indicated by their UniProt accession number. The dashed line on the heatmap indicates the two clusters of proteins.
Cancers 14 02001 g001
Figure 2. Exploratory analysis of proteome data from Proteome discoverer. (A) Principal Component Analysis of the urine proteome of the two groups. (B) The heatmap of proteins identified in all individuals. Samples are represented in columns and proteins in rows.
Figure 2. Exploratory analysis of proteome data from Proteome discoverer. (A) Principal Component Analysis of the urine proteome of the two groups. (B) The heatmap of proteins identified in all individuals. Samples are represented in columns and proteins in rows.
Cancers 14 02001 g002
Figure 3. Urinary protein levels of the candidate targets for PCa in the discovery group (using MS) and in the testing group (using immunoblot and immunoassay). MS: mass spectrometry.
Figure 3. Urinary protein levels of the candidate targets for PCa in the discovery group (using MS) and in the testing group (using immunoblot and immunoassay). MS: mass spectrometry.
Cancers 14 02001 g003
Figure 4. Exploratory analysis of proteogenome data from Pinnacle. (A) Principal Component Analysis of the urine proteogenome of the two groups. (B) The heatmap of mutant proteins identified in all individuals. Samples are represented in columns and proteins in rows. Proteins are identified by their gene name, and the mutation identified. The dashed line on the heatmap indicates the two clusters of proteins.
Figure 4. Exploratory analysis of proteogenome data from Pinnacle. (A) Principal Component Analysis of the urine proteogenome of the two groups. (B) The heatmap of mutant proteins identified in all individuals. Samples are represented in columns and proteins in rows. Proteins are identified by their gene name, and the mutation identified. The dashed line on the heatmap indicates the two clusters of proteins.
Cancers 14 02001 g004
Figure 5. PPI network of 86 mutated proteins with known expression in the prostate.
Figure 5. PPI network of 86 mutated proteins with known expression in the prostate.
Cancers 14 02001 g005
Table 1. Dysregulated proteins between PCa patients and cancer-free subjects (Proteome Discoverer).
Table 1. Dysregulated proteins between PCa patients and cancer-free subjects (Proteome Discoverer).
Uniprot IDProtein NameGene Namep-ValueCohen’s d
[Lower; Upper 95% CI]
P07288Prostate-specific antigenKLK30.004.21 (3.50; 4.91)
Q8WVN6Secreted and transmembrane protein 1SECTM10.01−2.16 (−2.39; −1.93)
P12830Cadherin-1CDH10.03−1.73 (−2.05; −1.41)
P0DOX5Immunoglobulin gamma-1 heavy chainN/A0.031.73 (1.39; 2.07)
Q12805EGF-containing fibulin-like extracellular matrix protein 1EFEMP10.03−1.68 (−2.25; −1.12)
P02766TransthyretinTTR0.031.66 (0.86; 2.46)
P01861Immunoglobulin heavy constant gamma 4IGHG40.041.52 (0.90; 2.15)
P01034Cystatin-CCST30.051.50 (0.91; 2.08)
Q01459Di-N-acetylchitobiaseCTBS0.05−1.44 (−1.86; −1.02)
The protein identification and label-free quantification performed by the Proteome Discoverer software revealed nine dysregulated proteins (p-value ≤ 0.05) between the tested groups. These proteins are shown in this table along with their p-value and effect size. The Cohen’s d for individual proteins is presented together with the lower and upper 95% confidence interval (CI). Abbreviation: Confidence interval (CI).
Table 2. Dysregulated proteins between PCa patients and cancer-free subjects (MaxQuant).
Table 2. Dysregulated proteins between PCa patients and cancer-free subjects (MaxQuant).
Uniprot IDProtein NameGene Namep-ValueCohen’s d
[Lower; Upper 95% CI]
Q8WVN6Secreted and transmembrane protein 1SECTM10.01−2.10 (−2.48; −1.73)
P07288Prostate-specific antigenKLK30.012.01 (1.08; 2.95)
P41222Prostaglandin-H2 D-isomerasePTGDS0.01−1.97 (−2.44; −1.49)
Q14624Inter-alpha-trypsin inhibitor heavy chain H4ITIH40.01−1.96 (−2.32; −1.60)
Q12805EGF-containing fibulin-like extracellular matrix protein 1EFEMP10.01−1.84 (−2.33; −1.35)
P55290Cadherin-13CDH130.02−1.75 (−2.11; −1.40)
P98160Basement membrane-specific heparan sulfate proteoglycan core proteinHSPG20.03−1.63 (−2.07; −1.19)
P04746Pancreatic alpha -amylaseAMY2A0.03−1.57 (−1.95; −1.19)
P01876Immunoglobulin heavy constant alpha 1IGHA10.041.55 (1.32; 1.78)
P02760Protein AMBPAMBP0.04−1.51 (−1.88; −1.13)
P12830Cadherin-1CDH10.05−1.48 (−1.90; −1.07)
Q12907Vesicular integral-membrane protein VIP36LMAN20.05−1.46 (−2.10; −0.83)
Q9NPP6Immunoglobulin heavy chain variantN/A0.041.58 (1.22; 1.93)
P02766TransthyretinTTR0.051.42 (0.97; 1.87)
The protein identification and label-free quantification performed by the MaxQuant software revealed fourteen dysregulated proteins (p-value ≤ 0.05) between the tested groups. These proteins are shown in this table along with their p-value and effect size. The Cohen’s d fof individual proteins is presented together with the lower and upper 95% confidence interval (CI). Abbreviation: Confidence interval (CI).
Table 3. List of mutations mapped on some proteins and respective mutant peptides identified in urine from PCa patients.
Table 3. List of mutations mapped on some proteins and respective mutant peptides identified in urine from PCa patients.
Uniprot IDProtein NameGene NameMutation DescriptionMutation TypeProtein Role in PCa or Other Types of Cancer
P02760Protein AMBPAMBPG238S; E192G; V69M; A286G;
P197S; R185Q;
G338S; G341A;
I198T; V313I;
G186R; R185Q
missenseAMBP is an inflammation-regulating protein, associated with human cancers [40,41], including PCa [42,43]. Increased urinary levels [6,42,44,45] but diminished levels in tumor prostate tissue have been reported in PCa patients [46,47,48].
P12830Cadherin-1CDH1H233R; A408EmissenseCDH1 is a protein implicated in cell adhesion, migration, and epithelial-mesenchymal transition [49,50] and its downregulation is correlated with a poor prognosis in PCa patients [51].
Q12805EGF-containing fibulin-like extracellular matrix protein 1EFEMP1V463MmissenseEFMP1 plays a role in cell adhesion and migration, acting as a tumor suppressor in PCa. Diminished EFEMP1 mRNA and protein levels [52] and EFEMP1 promoter hypermethylation were observed in PCa patients [53,54].
P98160Basement membrane-specific heparan sulfate proteoglycan core proteinHSPG2V4332I; A1503V;
S970F; M638V;
Q1062H
missenseHSPG2, found predominantly in the ECM and bone marrow, modulates tumor angiogenesis, proliferation, and differentiation. It is overexpressed in PCa tissues compared to non-malignant tissues, correlating with high GS and PCa cell proliferation and viability [55,56,57].
Q14624Inter-alpha-trypsin inhibitor heavy chain H4ITIH4R866C; G893SmissenseITIH4 is an acute-phase response protein whose function remains unclear [58]. Research points to a tumor suppressor activity of ITIH4 in human cancers and dysregulation in PCa [43,59].
P07288Prostate-specific antigen (PSA)KLK3C209Y; V55M; G156V; AVCG (47–50);
S117P; G87R;
L124F; A154T;
I179T
Missense; inframe_insertionPSA is widely used as serum biomarker for PCa. It was approved by the US Food and Drug Administration (FDA) in 1994 [60].
Q12907Vesicular integral-membrane protein VIP36LMAN2G250S; D229NmissenseLMAN2 protein is involved in endoplasmic reticulum to Golgi trafficking of some glycoproteins [61]. Dysregulation of the LMAN2 gene has been indicated in some cancers [62,63,64], while the role in PCa remains obscure. However, raised LMAN2 urinary levels were detected in PCa patients [44].
P41222Prostaglandin-H2 D-isomerasePTGDSL130MmissensePTGDS is involved in prostaglandins metabolism and lipid transport. The PTGDS gene is downregulated in malignant prostate tissues compared to non-malignant tissues and integrates a signature that predicts relapse after prostatectomy. In vitro, its overexpression increased death and suppressed the growth of PCa cells [65,66].
Q13510Acid ceramidaseASAH1V246AmissenseASAH1 hydrolyzes ceramide to sphingosine and fatty acid [67] and its protein levels are elevated in tumor prostate tissue [68]. Its increased levels have been suggested as a therapeutic target in PCa as they have been correlated with metastasis establishment and resistance to chemotherapy [69,70].
P08294Extracellular superoxide dismutase [Cu-Zn]SOD3A58TmissenseSOD3 is a known tumor suppressor gene in PCa. It is an antioxidant enzyme that catalyzes the dismutation of the superoxide radical anion [71]. SOD3-reduced levels were reported in PCa patients, and its overexpression in PCa cells prevented cell proliferation, migration, and invasion, suggesting a role as a therapeutic target and predictive marker [72,73].
P09211Glutathione S-transferase PGSTP1I105VmissenseGSTP1 is a known tumor suppressor gene in PCa and is responsible for cellular detoxification through glutathione conjugation [74]. PCa is characterized by loss of GSTP1 function, mostly due to hypermethylation of its regulatory CpG island [75], and it is purported to occur early in prostatic carcinogenesis [76,77].
P10451OsteopontinSPP1A22GmissenseSPP1 is a bone matrix protein involved in bone remodeling, modulation of inflammation, cell adhesion, and migration and angiogenesis [78]. In PCa, SPP1 is associated with metastasis and proliferation [79], lower overall survival and biochemical relapse-free survival, and high GS [80]. Higher SPP1 levels were reported in PCa patients [80,81,82].
P15309Prostatic acid phosphatasePAPG68DmissensePAP is one of the main secreted proteins by the prostate cells and was the first serum screening marker for PCa. PAP was latter replaced by PSA [83,84].
P25311Zinc-alpha-2-glycoproteinZAGP187L; A46TmissenseZAG promotes adipocyte lipolysis, resulting in cancer cachexia [85]. Elevated levels of this protein have been proposed as a serum marker for PCa [86,87], and a significant predictive ability was found for urinary ZAG [8].
Q4ZJI4Sodium/hydrogen exchanger 9B1SLC9B1N70SmissenseSLC9B1 is a Na+/H+ transporter responsible for preserving cellular homeostasis [88], but this transporter has not yet been correlated with any type of cancer.
Q9P2J8Zinc finger protein 624ZNF624S207FmissenseZNF624 has not been well studied yet, but in breast cancer was one of the target genes of a microRNA found to be significantly and independently correlated with patient prognosis [89].
Q6EMK4VasorinVASNR161QmissenseVASN, an inhibitor of TGF-beta signaling, is upregulated in PCa tissues and stimulates PCa proliferation [90].
P08174Complement decay-accelerating factorCD55S162LmissenseCD55 inhibits the complement system [91]. In PCa, CD55 mediates tumor cells survival and growth [92].
This table shows the UniProt IDs, protein and gene names, mutation site/description and type, and the role of proteins in PCa.
Table 4. Results of Polyphen-2 score and prediction for the mapped mutations.
Table 4. Results of Polyphen-2 score and prediction for the mapped mutations.
Gene NameMutationPredictionScoreSensitivitySpecificity
AMBPG238SProbably damaging1.0000.001.00
AMBPE192GProbably damaging0.750.9810.96
AMBPV69MPossibly damaging0.7580.850.92
AMBPA286GProbably damaging1.0000.001.00
AMBPP197SBenign0.0510.940.83
AMBPG338SProbably damaging0.9940.690.97
AMBPG341AProbably damaging0.9580.780.95
AMBPV313IBenign0.0250.950.81
AMBPG186RProbably damaging1.0000.001.00
AMBPR185QProbably damaging0.9920.700.97
CDH1H233RPossibly damaging0.8310.840.93
CDH1A408EPossibly damaging0.7980.840.93
EFEMP1V463MProbably damaging0.9990.140.99
HSPG2V4332IBenign0.0010.990.15
HSPG2A1503VProbably damaging1.000.001.00
HSPG2S970FPossibly damaging0.4980.880.90
HSPG2M638VBenign0.001.000.00
HSPG2Q1062HBenign0.001.000.00
ITIH4R866CProbably damaging10.001.00
ITIH4G893SBenign0.001.000.00
KLK3C209YProbably damaging1.0000.001.00
KLK3G156VProbably damaging1.0000.001.00
KLK3V55MProbably damaging0.9720.770.96
KLK3S117PPossibly damaging0.6210.870.91
KLK3G87RBenign0.1280.930.86
KLK3L124FProbably damaging1.0000.001.00
KLK3A154TPossibly damaging0.6570.860.91
KLK3I 179TPossibly damaging0.8000.840.93
LMAN2G250SProbably damaging1.000.001.00
LMAN2D229NProbably damaging0.9830.740.96
PTGDSL130MProbably damaging1.000.001.00
ASAH1V246ABenign0.001.000.00
SOD3A58TBenign0.1880.920.87
GSTP1I105VBenign0.001.000.00
SPP1A22GPossibly damaging0.6110.870.91
ACP3G68DProbably damaging1.000.001.00
AZGP1P187LProbably damaging0.940.690.97
AZGP1A46TBenign0.0020.990.30
SLC9B1N70SBenign0.0360.940.82
ZNF624S207FBenign0.2140.920.88
VASNR161QBenign0.0190.950.80
CD55S162LProbably damaging0.9900.720.97
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lima, T.; Barros, A.S.; Trindade, F.; Ferreira, R.; Leite-Moreira, A.; Barros-Silva, D.; Jerónimo, C.; Araújo, L.; Henrique, R.; Vitorino, R.; et al. Application of Proteogenomics to Urine Analysis towards the Identification of Novel Biomarkers of Prostate Cancer: An Exploratory Study. Cancers 2022, 14, 2001. https://doi.org/10.3390/cancers14082001

AMA Style

Lima T, Barros AS, Trindade F, Ferreira R, Leite-Moreira A, Barros-Silva D, Jerónimo C, Araújo L, Henrique R, Vitorino R, et al. Application of Proteogenomics to Urine Analysis towards the Identification of Novel Biomarkers of Prostate Cancer: An Exploratory Study. Cancers. 2022; 14(8):2001. https://doi.org/10.3390/cancers14082001

Chicago/Turabian Style

Lima, Tânia, António S. Barros, Fábio Trindade, Rita Ferreira, Adelino Leite-Moreira, Daniela Barros-Silva, Carmen Jerónimo, Luís Araújo, Rui Henrique, Rui Vitorino, and et al. 2022. "Application of Proteogenomics to Urine Analysis towards the Identification of Novel Biomarkers of Prostate Cancer: An Exploratory Study" Cancers 14, no. 8: 2001. https://doi.org/10.3390/cancers14082001

APA Style

Lima, T., Barros, A. S., Trindade, F., Ferreira, R., Leite-Moreira, A., Barros-Silva, D., Jerónimo, C., Araújo, L., Henrique, R., Vitorino, R., & Fardilha, M. (2022). Application of Proteogenomics to Urine Analysis towards the Identification of Novel Biomarkers of Prostate Cancer: An Exploratory Study. Cancers, 14(8), 2001. https://doi.org/10.3390/cancers14082001

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop