Next Article in Journal
A Review of Fibraurea tinctoria and Its Component, Berberine, as an Antidiabetic and Antioxidant
Previous Article in Journal
Unveiling the Efficacy of Sesquiterpenes from Marine Sponge Dactylospongia elegans in Inhibiting Dihydrofolate Reductase Using Docking and Molecular Dynamic Studies
Previous Article in Special Issue
Lectin-Based Affinity Enrichment and Characterization of N-Glycoproteins from Human Tear Film by Mass Spectrometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Glycoproteomic Analysis of Urinary Extracellular Vesicles for Biomarkers of Hepatocellular Carcinoma

1
Center for Supramolecular Chemical Biology, School of Life Sciences, Jilin University, Changchun 130012, China
2
Prenatal Diagnosis Center, Reproductive Medicine Center, The First Hospital of Jilin University, Changchun 130021, China
3
Department of Hepatopancreatobiliary Medicine, The Second Hospital, Jilin University, Changchun 130041, China
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(3), 1293; https://doi.org/10.3390/molecules28031293
Submission received: 26 December 2022 / Revised: 21 January 2023 / Accepted: 24 January 2023 / Published: 29 January 2023
(This article belongs to the Special Issue Frontiers in Mass Spectrometry Based Glycomics and Glycoproteomics)

Abstract

:
Hepatocellular carcinoma (HCC) accounts for the most common form of primary liver cancer cases and constitutes a major health problem worldwide. The diagnosis of HCC is still challenging due to the low sensitivity and specificity of the serum α-fetoprotein (AFP) diagnostic method. Extracellular vesicles (EVs) are heterogeneous populations of phospholipid bilayer-enclosed vesicles that can be found in many biological fluids, and have great potential as circulating biomarkers for biomarker discovery and disease diagnosis. Protein glycosylation plays crucial roles in many biological processes and aberrant glycosylation is a hallmark of cancer. Herein, we performed a comprehensive glycoproteomic profiling of urinary EVs at the intact N-glycopeptide level to screen potential biomarkers for the diagnosis of HCC. With the control of the spectrum-level false discovery rate ≤1%, 756 intact N-glycopeptides with 154 N-glycosites, 158 peptide backbones, and 107 N-glycoproteins were identified. Out of 756 intact N-glycopeptides, 344 differentially expressed intact N-glycopeptides (DEGPs) were identified, corresponding to 308 upregulated and 36 downregulated N-glycopeptides, respectively. Compared to normal control (NC), the glycoproteins LG3BP, PIGR and KNG1 are upregulated in HCC-derived EVs, while ASPP2 is downregulated. The findings demonstrated that specific site-specific glycoforms in these glycoproteins from urinary EVs could be potential and efficient non-invasive candidate biomarkers for HCC diagnosis.

1. Introduction

Liver cancer is the sixth most commonly diagnosed cancer and the third leading cause of cancer death worldwide in 2020, with approximately 906,000 new cases and 830,000 deaths [1]. It is estimated that, by 2025, more than 1 million individuals will be affected by liver cancer annually [2]. Hepatocellular carcinoma (HCC) accounts for 75–85% of primary liver cancer cases and constitutes a major health problem worldwide [3]. Although patients with HCC can be treated with local ablation, surgical resection, and liver transplantation, due to the lack of specific clinical symptoms and efficient diagnostic method in the early stage, the majority of patients are already in advanced stages when they are diagnosed, thus having poor prognosis. For example, the 5-year survival rate for HCC dropped from 70% for early-stage patients to 16% for late-stage patients [4]. Therefore, early detection of HCC is of particular importance to reduce the mortality rates of HCC patients. Currently, serum α-fetoprotein (AFP) measurement and liver ultrasonography are widely accepted as the most effective and affordable tools to screen and diagnose HCC in the clinic [5]. However, not all patients with HCC have elevated serum AFP, and the sensitivity and specificity of AFP are limited to 40–62% [6]. Furthermore, tumors can only be detected by imaging technology if they are greater than 1 cm in diameter [7]. Therefore, there is an urgent need to develop non-invasive biomarkers that are more sensitive for HCC diagnosis at the early stage.
In recent years, liquid biopsy has gained much attention as an alternative biomarker to conventional biomarkers and tissue biopsy for cancer diagnosis and monitoring. Liquid biopsy represents a minimally invasive procedure that is usually performed on biofluid samples to detect circulating tumor by-products, which includes circulating tumor cells (CTCs), cell-free DNA (cfDNA), and extracellular vesicles (EVs) [8,9,10]. Compared with CTCs and cfDNA, EVs have attracted researchers’ interest because they are not only present in circulation at relatively early stages of disease, but also with high abundance and high stability [11,12]. EVs are heterogeneous populations of phospholipid bilayer-enclosed vesicles that are released from almost all types of cells and can be found in biological fluids such as blood, saliva, urine, breast milk, and cerebrospinal fluid (CSF) [13,14]. According to their biophysical properties and biogenesis pathways, EVs are generally grouped into two subclasses, exosomes and microvesicles. Exosomes are small EVs with a size range of ~30–150 nm that derived from the endosomal membrane pathway, while microvesicles are larger than exosomes (100–1000 nm) and are formed by budding directly from the cell membrane [11]. EVs play a central role in intercellular communication by selectively packaging and transporting bioactive cargoes such as proteins, RNA and metabolites to recipient cells, and participating in immune responses, cardiovascular diseases and cancer [15,16]. Recent studies have found that cancer cell-derived EVs released into the tumor microenvironment and circulation can promote tumor progression and metastasis by inducing matrix remodeling, angiogenesis, inflammation, and metastatic niche formation [17]. Furthermore, multiple studies on EV-based diagnosis of cancer have suggested that EV transcriptomic and proteomic biomarkers can increase the likelihood of early and sensitive detections of cancer. For example, plasma EV-derived miRNA-483-3 could be used as a potential biomarker for the diagnosis of early-stage small cell lung cancer, and miRNA-152-3p and miRNA-1277-5p could be used for the diagnosis of early-stage non-small cell lung cancer [18]. Glypican 1 (GCP1) and macrophage migration inhibitory factor (MIF) were identified as potential exosome-associated biomarkers for early diagnosis of pancreatic cancer [19,20]. As has been reported, tissue- and plasma-derived exosomal hsa-miR-483-5p could differentiate HCC and non-HCC cases and may represents a potential specific and sensitive biomarker for HCC diagnosis [21]. However, while these studies mainly focus on the utilization of proteins or RNAs in EVs as biomarkers for disease diagnosis, little attention is paid to researching the use of post-translational modifications (PTMs) of proteins in EVs as biomarker.
Protein glycosylation is one of the most important PTMs found in eukaryotic cells, which is involved in molecular recognition, immune response, and intercellular communication [22,23,24]. Notably, aberrant glycosylation of proteins is a hallmark of cancer, and glycoproteins account for the majority of FDA-approved biomarkers for cancer diagnosis and monitoring [25]. As has been reported, bifucosylated N-glycans with both core and antennary fucosylation were elevated in HCC patients as compared to patients with cirrhosis [26]. Moreover, AFP-L3, a corefucosylated variant of AFP which reacts with LCA fraction, was known to be more specific to HCC than AFP and was a useful biomarker for HCC diagnosis and assessment [27]. In recent years, studies on EV glycosylation have gradually increased. Because EVs originate from the plasma membrane or from multivesicular endosomes, their surface is highly glycosylated, and specific cancer-associated glycosignatures have already been identified in cancer EVs, thus highlighting EV glycosylation as a potential source of novel circulating biomarkers in cancer diagnosis [28,29]. Currently, studies have been reported on the glycosylation of EVs derived from different tumor cell lines and biofluids. For example, the proteoglycan glypican-1 positive (GPC1+) exosomes derived from the serum of pancreas cancer patients could distinguish early- and late-stage pancreas cancer as well as patients with a benign pancreas disease [30]. Furthermore, Sakaue et al. found that increased levels of sialic acids of CD133 glycoprotein from the pancreatic cancer patient’s ascites-derived exosomes, and indicated that highly glycosylated CD133 could be a potential biomarker for advanced pancreatic cancer [31]. Moreover, the proteoglycan biglycan (BGN) and the glycoprotein basigin (EMMPRIN) were found to be enriched in pancreatic tumor tissue-derived EVs [32]. In addition, high levels of EMMPRIN were observed from breast cancer patient serumderived MVs [33]. Notably, galectin-3-binding protein (LG3BP) was identified in several tumor-derived EVs, including ovarian cancer, uveal melanoma, and gliomas [34,35,36]. With respect to HCC, there are currently only two glycosylation studies on EVs that have been reported, one from serum-derived exosome, and the other from cell line-derived exosome. Lv et al. performed a reverse capture strategy to profile the N-glycome of exosomes isolated from HCC patients’ serum samples, and compared this to the N-glycans identified in exosomes from healthy samples. The majority of the N-glycans detected in exosomes from HCC patients were modified with sialic acids or fucoses [37]. Wang et al. found that compared to exosomes derived from naive HCC cells, α2,6-sialylation degradation abolished both the proliferation-promoting and migration-promoting effects of HCC-exo, which suggests that a loss of α2,6-sialylation decreases HCC progression through the loss of cancer cell-derived exosomes [38]. The above studies indicate that, compared with normal EVs, specific glycosylation modifications are existing on the surface of tumor cell-derived EVs, therefore they can be used as potential biomarkers for disease diagnosis and prognosis. Nevertheless, previously studies on glycosylation of EVs mainly focus on the identification of released glycans or deglycosylated peptides, losing the site-specific information of each individual glycan. Little is known about the information on the glycosylation of tumor-derived EVs at the intact glycopeptide (including both the glycosylation sites and glycan structures) level.
In this study, we investigated the glycoproteomic profiling of EVs from HCC patient urine samples to screen the potential candidate biomarker. Urinary EVs were first isolated by a previously reported EVTRAP method [39]. Then, intact-glycoproteomic analysis of EVs glycoproteins were detected using the C18-HILIC-MS/MS analysis with the N-glycan database engine GPSeeker search [40]. In total, 756 intact N-glycopeptides mapping to 154 N-glycosites on 107 intact N-glycoproteins were identified. Furthermore, variations in the N-glycans and site-specific glycans were elucidated. With the criteria of differentially expressed intact N-glycopeptides (DEGPs) as fold change cutoff ratio ≥2.0 or ≤0.5, 344 DEGPs were identified from EVs in HCC relative to the NC group, corresponding to 308 upregulated and 36 downregulated N-glycopeptides, respectively. Moreover, the relationship between aberrated glycoproteins in EVs and HCC was further discussed.

2. Results and Discussion

2.1. Workflow of EVs Proteomics and Glycoproteomics from Urine Samples

The workflow of EVs proteomics and glycoproteomics analysis from urine samples is shown in Figure 1. In brief, EVs were first isolated by EVTRAP beads using pooled urine samples from HCC patients and from NC, respectively. Then, urinary EVs were lysed and the extracted proteins were digested into peptides; thereafter, 5% of the digested peptides were directly used to perform proteomic analysis with nanoflow LC-MS/MS. The other 95% portion of peptides were labeled with isotopomeric dimethyl labels and mixed at a 1:1 ratio from HCC and NC, and the N-glycopeptides enrichment was performed by the HILIC method. Finally, the N-glycopeptides were analyzed by nanoflow LC-MS/MS and data were analyzed with the GPSeeker search engine.

2.2. General Description of Extracellular Vesicles Isolated by EVTRAP

Compared to the ultracentrifugation method, the EVTRAP method had a high recovery yield (over 95%) from urine samples [39], and had the capacity to uncover more low-abundant urinary biomarkers, especially EV proteins with PTMs such as phosphorylation and glycosylation. The quality of the isolated EVs was evaluated by proteins identified in the proteomic analysis with mass spectrometry. In three replicate tests, 1331 and 1275 proteins were identified in urinary EVs from the HCC and NC group, respectively (Supplementary Table S1). The quality of the isolated EVs was evaluated by comparing the proteins identified in EVs by direct mass spectrometry with that in the Vesiclepedia database [41], and up to 96% of the proteins were matched form the HCC group (Figure 2a). Moreover, 94 proteins were overlapped with the top 100 most common identified EV marker proteins, including EV markers such as CD9, CD63, Alix, TSG101, annexin A1 and syntenin-1 (Figure 2b). In addition, Gene Ontology (GO) analysis of cellular component revealed “extracellular exosome” as the top enriched GO item (Figure 2c). The quality evaluation results of EVs from the NC group was shown in Figure S1. Therefore, after quality assessment of the isolated EVs, we performed comprehensive glycoproteomic analysis of the EVs at the intact N-glycopeptide level and evaluated its potential role as a biomarker for the hepatocellular carcinoma.

2.3. Comprehensive Analysis of Intact N-glycopeptides in the Urinary EVs

With target and decoy database searches using the intact N-glycopeptide search engine GPSeeker, a total of 756 intact N-glycopeptides were identified, corresponding to 154 N-glycosites on 107 intact N-glycoproteins with 286 putative N-glycan linkages from 94 monosaccharide compositions (Figure 3a, Table S2). All of the identified N-glycosites were observed with N-X-S and N-X-T sequons, in which ~60% of the glycosites had the N-X-T sequons and 40% had N-X-S sequons (Figure 3b). In addition, the percentage of mannose, hybrid, and complex N-glycosylation among the 756 intact N-glycopeptide IDs is 30.3%, 14.9% and 54.8%, respectively. Furthermore, compared with the NC-derived EVs, the percentage of mannose N-glycosylation in HCC-derived EVs was increased while the hybrid and complex N-glycosylation showed the opposite trend (Figure 3c). Moreover, for the monosaccharide compositions of the N-glycan moieties among the identified intact N-glycopeptides, the composition containing two HexNAc was the most common, and the eight most common monosaccharide compositions and their occurrence numbers are shown in Figure S2. Macro-heterogeneity and micro-heterogeneity are two key features of protein glycosylation; therefore, we further performed the heterogeneity of glycosylation analysis with our data set. As for macro-heterogeneity, approximately 30% of the glycoproteins were identified to have more than one N-glycosite; whereas for micro-heterogeneity, ~46% of the glycosites were occupied by more than one N-linked glycan compositions (Figure 3d). For example, on probable serine carboxypeptidase (Q9H3G5, CPVL_HUMAN), the same N-glycan composition N2H5F0S0 was identified on three N-glycosites (N81, N307, and N346); while on the same N-glycosite N346, other N-glycan compositions (N2H3F0S0, N2H4F1S0 and N2H4F0S0) were also identified (Table 1).

2.4. Relatively Quantification Analysis of the Site-Specific Glycans in the Urinary EVs

To further analyze the site-specific glycans in urinary EVs, the abundance of the 756 intact N-glycopeptide IDs together with their isotopic pairs (4.0134 Da) in the corresponding MS spectra were searched with GPSeekerQuan. With the criteria of observation of all of the six most abundant isotopic peaks; and with the designated differentially expressed intact N-glycopeptides (DEGPs) as fold change cutoff ratio ≥2.0 or ≤0.5, a total of 344 DEGPs were identified corresponding to 308 upregulated and 36 downregulated N-glycopeptides, respectively (Table S3). The 308 upregulated intact N-glycopeptides come from 51 N-glycoproteins and the 36 downregulated intact N-glycopeptides come from 7 N-glycoproteins. For the 51 intact N-glycoproteins with up-regulation, 38 were quantified with more than one DEGPs. Among the upregulated N-glycopeptides, intact N-glycopeptide NVTLLSR_N2H9F0S0 from N-glycosite N828 of N-glycoprotein Glutamyl aminopeptidase (AMPE_HUMAN, Q07075) displayed the highest up-regulation, with a fold change of 40 in HCC relative to NC (Figure S3). For Galectin-3-binding protein (LG3BP_HUMAN, Q08380), on N-glycosite N398, intact N-glycopeptide GLNLTEDTYKPR_N2H3F1S0 was found to be 9.31-fold upregulated in HCC relative to NC (Figure 4), while on N-glycosite N551, intact N-glycopeptides AAIPSALDTNSSK with compositions of N2H3F1S0, N2H5F0S0 and N6H6F1S0 were identified with a range of upregulated (2.21–9.31). LG3BP is a large oligomeric glycoprotein present in human body fluids and originally identified as a tumor-secreted antigen from breast and lung cancer [42,43]. Upregulation of LG3BP was previously observed in the serum of patients with HCC, and high expression of LG3BP was as a poor prognostic biomarker for HCC [44,45]. Recently, proteomic analysis on serum exosomes revealed that the LG3BP has higher diagnostic capacity than AFP, and has the potential as a candidate biomarker for HCC detection [46,47]. In addition, upregulation of intact N-glycopeptideWNNTGCQALPSQDEGPSK_N2H5F0S0 at N-glycosite N499 of the polymeric immunoglobulin receptor (PIGR_HUMAN, P01833) was quantified with a fold change of 2.25 in HCC relative to NC (Figure S4). PIGR is universally expressed in epithelial cells and regulates transcytosis of dimeric IgA and pentameric IgM, which are the first-line antibodies in response to initial infection. Previously, studies have shown that PIGR was significantly overexpressed in tumor tissues and in serum samples in HCC patients, and the level of PIGR was statistically significantly associated with early recurrence in early-stage HCC and in hepatitis B surface antigen-positive HCC patients [48,49]. Previous studies have shown that elevated PIGR levels were found in the serum EVs from patients with HCC, and EV-PIGR could promote liver cancer cell aggressiveness through activating thePDK1/Akt/GSK3β/β-catenin signaling axis [46,50]. Therefore, EV-PIGR holds the potential as a diagnostic and prognostic marker in HCC. For Kininogen-1 (KNG1_HUMAN, P01042), intact N-glycopeptides with compositions of N6H3F0S0, N5H4F1S1 and N6H3F2S0 were identified with a range of upregulated (2.27–3.39) in HCC relative to NC. KNG1 is a precursor for kinins and plays a central role in the blood coagulation system and in the kinin–kallikrein system [51]. Previous studies have shown that the serum level of KNG1 was significantly elevated in patients with various cancers, including advanced colorectal adenoma (ACA), colorectal cancer (CRC), and HCC [52,53]. Hence, KNG1 is a potential serum biomarker for the early detection of ACA, CRC and HCC. Moreover, we found that the N-glycan of KNG1 possesses fucose modification, and this is consistent with the results from Wang et al. [54], who also revealed an increased amount of both core and outer-arm fucosylation of KNG1 in patients with HCC. In addition, they revealed that the increased reactivity with fucose-binding lectin observed in patients with HCC is the result of increased levels of fucosylation. Therefore, they proposed a biomarker algorithm that use the combination of fucosylated KNG1, AFP, and clinical characteristics to improve the detection of early-stage HCC.
For the seven intact N-glycoproteins with down-regulation, all of them were quantified with more than one DEGP. Among the downregulated N-glycopeptides, intact N-glycopeptide SVNVSLTQEELDSGLDELSVR_N4H5F0S0 from N-glycosite N463 of A-kinase anchor protein 2 (AKAP2_HUMAN, Q9Y2D5) displayed the lowest down-regulation, with the fold change of 0.18 in HCC relative to NC (Figure S5). Intact N-glycopeptideKNQSSEDILR_N4H5F1S1 at N-glycosite N478 of apoptosis stimulating of p53 protein 2 (ASPP2_HUMAN, Q13625) was downregulated with a fold change of 0.34 in HCC relative to NC (Figure 5). ASPP2, a member of the ankyrin-repeat, SH3-domain, and proline-rich region containing protein (ASPP) family, which also include ASPP1 and inhibitory ASPP (iASPP), is a key regulator of various cancer-relevant cellular phenotypes and behaviors, including cell proliferation, apoptosis, and polarity [55,56]. ASPP2 is a haplo-insufficient tumor suppressor, and is frequently downregulated in a variety of human cancers, including breast cancer, lung cancer and leukemia [57]. Furthermore, a previous study by Zhao et al. [58] also found that ASPP2 is downregulated in hepatocellular carcinoma owing to DNA methylation. Chen et al. [59] revealed that downregulation of ASPP2 may contribute to tumor progression and chemoresistance via promoting BECN1-dependent autophagy in HCC. Moreover, studies have shown that reduced ASPP2 expression results in EMT, and is associated with poor survival in hepatocellular carcinoma patients [60]. Therefore, ASPP2 plays important roles in the development of HCC, which may be a potential target for the treatment of HCC.
Interestingly, there were 5 glycoproteins with both upregulated and downregulated site-specific glycans. For example, on Alpha-2-macroglobulin-like protein 1 (A2ML1_HUMAN, A8K2U0), intact N-glycopeptide LGHI867NFTISTK_N2H5F0S0 was quantified with up-regulation while THHW857NITAVK_N2H8F0S0 was significantly downregulated. Furthermore, there were N-glycosites with both upregulated and downregulated glycans. For example, on glycoprotein Vasorin (VASN_HUMAN, Q6EMK4), on the same glycosylation site N500, the upregulated intact N-glycopeptide500NLSGPDK_N2H5F0S0 and 500NLSGPDK_N3H6F0S0 as well as downregulated 500NLSGPDK_N3H5F0S1 were observed (Figure 6). VASN was found to be highly expressed in the cells and tissues of HCC, and at a high level in the serum samples of HCC patients [61]. Being a transmembrane glycoprotein, VASN regulates cellular responses to vascular lesion, and exosomal VASN from HepG2 cells could promote the migration of human umbilical vein endothelial cells (HUVEC) [62]. These results suggest that VASN may play an important role in the HCC pathogenesis.

2.5. Gene Ontology Analysis of the Differentianl Expressed Intact N-glycopeptides

To further investigate the biological information of the differentially expressed N-glycoproteins identified in the urinary EVs, GO analysis was carried out using DAVID Bioinformatics Resources 6.8 for the differentially expressed intact N-glycopeptides. As is shown in Figure 7a, they are mainly localized on the cell, extracellular region, organelle and membrane regions in cellular components. However, upregulation was more circulated in the extracellular region and protein-containing complex region, down-regulation was also circulated in the cell junction region. As far as molecular function is concerned, it was shown that most of them participate in catalytic activity, molecular function regulator and binding, while up-regulation also played a role in transporter activity and molecular transducer activity (Figure 7b). In the aspect of biological process, both up-regulation and down-regulation are involved in the cellular process, metabolic process, cellular component organization or biogenesis, and biological regulation (Figure 7c). The effect of site-specific glycoforms on these important biological processes and molecular functions needs further investigation.
Despite the valuable results of this study, our research also has several limitations. Firstly, the HCC patient sample size was relatively small, the glycoproteins identified in urinary EVs as potential candidate biomarkers for HCC diagnosis should be further validated in larger samples. Secondly, the candidate glycoproteins were not validated by Western blot or other technologies, therefore, there is a need to evaluate the precision and robustness of the identified candidate glycoproteins as non-invasive biomarkers for HCC in any further study.

3. Materials and Methods

3.1. Patients and Sample Collection

This study was approved by the Ethics Committee of the Second Hospital of Jilin University and written informed consent was obtained from all the participants. Urine samples were obtained from 21 HCC patients (Table S1) as well as seven normal controls. The diagnosis of HCC was determined pathologically and immunohistochemically according to the WHO classification. The first morning urine samples were collected and kept for less than 1 h at room temperature followed by centrifugation at 2000× g for 10 min to remove cell debris and large apoptotic bodies, and then stored at − 80 °C until use.

3.2. Isolation of Extracellular Vesicles from Urine by EVTRAP

First, pooled urine samples were made by aspirating equal volumes of urine samples from HCC patients or from NC before EV isolation. Then, 10 mL of the pooled urine samples were used for EV isolation by EVTRAP method as we previously described [39]. In brief, EVTRAP beads was mixed with the urinary sample at a beads-to-urine ratio of 1:50 (v/v). The sample was incubated for 1 h with end-over-end rotation at room temperature, then a magnetic separator was used to remove the solution. Next, 150 μL of fresh 100 mM TEA solution was used twice to elute the EVs from beads, and the eluate was dried out in a vacuum centrifuge (Labconco, SpeedVac).

3.3. Enzyme Digestion and Protein Extraction of Urinary EVs

The dried urinary EVs were lysed to extract proteins using a phase-transfer- surfactant (PTS)-aided procedure as previously reported [63]. First, 100 μL of lysis buffer (50 mM Tris-HCl, 12 mM SDC, 12 mM SLS, 10 mM TCEP, 40 mM CAA, and 1% phosphatase inhibitor mixture in in 50 mM Tris·HCl) was used to solubilize the dried EVs. Second, the solution was boiled at 95 °C for 10 min and diluted five-fold with 50 mM TEAB. Then, the proteins were digested with Lys-C at an enzyme/proteome ratio of 1:100 (w/w) for 3 h at 37 °C. Afterwards, trypsin was added to a final 1:50 (w/w) enzyme/proteome ratio for overnight digestion at 37 °C. The digested peptides were acidified with trifluoroacetic acid to a final concentration of 0.5% trifluoroacetic acid, and 500 μL of ethyl acetate was added. After vortexing for 2 min it was centrifuged at 20,000× g for 2 min to obtain aqueous and organic phases. The aqueous phase was collected and desalted using the Sep-pak C18 column (Waters). Each sample was split into 5% and 95% aliquots for the proteomic and glycoproteomic experiments, respectively.

3.4. Stable Isotopic Dimethyl Labeling and Enrichment of N-glycopeptides from Urinary EVs

For glycoproteomic analysis, the 95% portion of each peptide sample was first labeled with isotopomeric dimethyl labels, then the N-glycopeptide enrichment was performed by the HILIC method. The stable isotope dimethyl labeling was performed as previously described with slight modifications [64]. Briefly, 200 μL of 100 mM TEAB (pH = 8.0) was first used to reconstitute the above desalted EV peptides, and 16 μL of CD2O (4%, v/v) and CH2O (4%, v/v) was added to the peptides derived from HCC and NC, respectively. After shaking for 1 min, 16 μL of 0.6 M NaBD3CN/NaBH3CN solution was added, and the solutions were incubated for 1 h in a shaker at room temperature. Next, 40 μL of 1M Tris (pH = 7.0) was used to terminate the reaction and 32 μL of formic acid was added to further quench the reaction. Afterward, the obtained light- and heavy-labeled peptides were mixed at a 1:1 ratio, desalted and dried in the SpeedVac.
The N-glycopeptides enrichment was performed with the centrifugation-assisted click maltose–HILIC approach as previously reported [65]. In brief, the dried dimethyl labeled peptide sample was reconstituted in loading buffer (80% ACN/1% TFA). Next, the peptide sample was mixed with the click maltose material in a new tube to a final volume of 200 μL and incubated for 1 h at room temperature. Then the mixture was pipetted into a 200 μL tip with a 1 mm Teflon disk blocking in the tip, and after centrifugation, the loading buffer was spun-out from HILIC tip, and the HILIC tip was washed twice with 200 μL 80% ACN/1% TFA. Finally, the enriched N-glycosylated peptides were eluted with 100 μL of 30% ACN/0.1% FA and dried with a vacuum freeze centrifuge.

3.5. LC−MS/MS Analysis

The dried peptides and glycopeptides were resuspended in 15 μL of 0.1% formic acid (FA) and 2% ACN solution. LC-MS/MS analysis was performed using a Dionex UltiMate 3000 RSLCnano system (Thermo Fisher Scientific, Bremen, Germany) coupled with an Orbitrap Fusion Lumos Tribrid spectrometer (Thermo Fisher Scientific, Bremen, Germany). The samples were separated on a 15 cm in-house packed column (150 μm i.d.) containing C18 AQ beads (1.9 μm, Dr. Maisch, GmbH, Germany). Mobile phase A was composed of 98% H2O and 2% ACN and mobile phase B was composed of 80% ACN and 20% H2O, both containing 0.1% FA. For the proteomic analysis, the flow rate was maintained at 600 nL/min with a linear 100 min gradient. The elution gradient was kept at 1% mobile phase B for 10 min, and then increased linearly from 1% to 34% mobile phase B over 60 min, followed by increasing linearly from 34% to 90% mobile phase B over 15 min. After remaining at 90% mobile phase B for 12 min, the content of mobile phase B was returned to 1% and maintained for 3 min. The mass spectrometer was operated in data-dependent mode. Full mass scan MS spectra (m/z 350−1800) were acquired by the Orbitrap mass analyzer with a resolution of 240,000. The AGC target for MS1 was set as 4 × 105 and a maximum injection time of 50 ms. The most intense ions above a threshold ion count of 4 × 103 and charge state ≥2 were selected and fragmented using HCD fragmentation with 30% normalized collision energy. The dynamic exclusion was set for 60 s with a 10 ppm mass window. Fragment ion spectra were acquired in the linear IT with an AGC of 2 × 104 and a maximum injection time of 35 ms for IT MS2 detection. For the glycoproteomic analysis, the elution gradient was kept at 4% mobile phase B for 3 min, and then increased linearly from 4% to 34% mobile phase B over 60 min, followed by increasing linearly from 34% to 90% mobile phase B over 23 min. After remaining at 90% mobile phase B for 10 min, the content of mobile phase B was returned to 1% and maintained for 4 min. The flow rate was maintained at 800 nL/min. The parent ion is selected in the Orbitrap cell (FTMS) at a resolution of 120,000. Up to the top 20 most abundant isotope patterns with a charge ≥ +2 from the survey scan (350–1500 m/z) were selected with an isolation window of 1.6 m/z and fragmented by HCD with normalized collision energies of 27. The AGC target value for MS1 and MS2 scan modes was set to 5 × 105 and 1 × 105, respectively. System control and data collection were carried out by Xcalibur software.

3.6. Protein Identification

RAW files were searched against the Uniprot human database (www.uniprot.org, accessed on 20 March 2022) using MaxQuant (version 2.2) and the default search parameter settings: 10 ppm and 20 ppm mass tolerance for precursor and fragment ions, respectively; fixed modifications: carbamidomethylation of cysteine (+57.0214 Da); variable modifications: oxidation of methionine (+15.9949 Da); fully Trypsin/P digestion with up to two missed cleavages. The false discovery rates (FDR) of proteins and peptides were set at 0.01.

3.7. Intact N-glycopeptideIdentification and Quantification

Database search and intact N-glycopeptides identification was performed using GPSeeker as previously described [40]. In brief, the human proteome DB (UniProt, 20,375 entries) and the human N-glycome DB (75,888 entries) were used to build the theoretical customized human intact N-glycopeptides databases [66], and the LC-MS/MS raw dataset was searched against the four databases independently. To search matching precursor and fragment ions, the isotope peak abundance cut-off, isotope peak mass-to-charge ratio (m/z) deviation and isotope peak abundance deviation were set to 40%, 20 ppm and 50%, respectively. The search of intact N-glycopeptide spectrum matches (GPSMs) included the following parameters: Y1 ion, Top5; the minimum percentage of matched fragment ions of the peptide backbone, ≥10%; the minimum matched fragment ion of the N-glycan moiety, ≥1; TopN hits, n = 2 with Top1 hit(s) having the lowest P score; G-bracket, ≥1; and GF score, ≥1. For each dataset, the GPSMs from the target and decoy searches were combined and ranked with the increasing order of the P score; a cutoff P score was chosen to obtain spectrum-level FDR ≤1%. Target GPSMs with FDR control were subjected to duplicate removal to obtain the final intact N-glycopeptide IDs. GPSeekerQuan software was used to quantify the differentially expressed intact N-glycopeptides (DEGPs) from HCC relative to NC samples. For each intact N-glycopeptide ID, GPSeekerQuan searched its paired precursor ion with a mass difference of 4.01344 Da and an isotopic peak m/z tolerance of 20 ppm. For each precursor ion, the summed abundance of Top 3 isotopic peaks was used for relative quantitation. For each intact N-glycopeptide ID, all the six isotopic peaks in the pair are required to be observed to obtain the relative ratio (HCC/NC).

3.8. Bioinformatics Analysis

The functional gene enrichment analysis and Vesiclepedia database search was performed with the FunRich software (version 3.1.3) [67]. Gene ontology data were analyzed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) bioinformatics tool (version 6.8, https://david.ncifcrf.gov, accessed on 20 March 2022) [68].

4. Conclusions

In this study, we performed a comprehensive glycoproteomic profiling of urinary EVs at the intact N-glycopeptide level to screen potential biomarkers for the diagnosis of HCC. Combined with DB search using GPSeeker, 756 intact N-glycopeptides with 154 N-glycosites, 158 peptide backbones, and 107 N-glycoproteins were identified. Out of 756 intact N-glycopeptides, 344 DEGPs were identified, corresponding to 308 upregulated and 36 downregulated N-glycopeptides, respectively. Compared to NC, several aberrated glycoproteins were identified in urinary EVs from patients with HCC. The glycoproteins LG3BP, PIGR and KNG1 are upregulated in HCC-derived EVs, while ASPP2 is downregulated. The findings demonstrated that specific site-specific glycoforms in these glycoproteins from urinary EVs could be potential and efficient non-invasive candidate biomarkers for HCC diagnosis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28031293/s1, Figure S1: Evaluation of the urinary EVs isolated by EVTRAP method derived from normal controls (NC); Figure S2: The eight most common N-glycan monosaccharide compositions vs. the occurrence number of intact N-glycopeptide IDs identified in this study; Figure S3: Upregulation of intact N-glycopeptide NVTLLSR_N2H9F0S0 in urinary EVs from patients with HCC relative to NC; the N-glycosite is N828 on Glutamyl aminopeptidase (AMPE_HUMAN, Q07075); Figure S4: Upregulation of intact N-glycopeptide WNNTGCQALPSQDEGPSK_N2H5F0S0 in urinary EVs from patients with HCC relative to NC; the N-glycosite is N499 on polymeric immunoglobulin receptor (PIGR_HUMAN, P01833); Figure S5: Downregulation of intact N-glycopeptide SVNVSLTQEELDSGLDELSVR_N4H5F0S0 in urinary EVs from patients with HCC relative to NC; the N-glycosite is N463 on A-kinase anchor protein 2 (AKAP2_HUMAN, Q9Y2D5). Table S1: Clinical characteristics of the HCC patient enrolled in this study; Table S2: Proteins identified in the urinary EVs by the proteomic study; Table S3: Detail information of the intact N-glycopeptides in the urinary EVs identified by the N-glycoproteomic study; Table S4: Differentially expressed intact N-glycopeptides quantitated with ≥2.0- or ≤0.5-fold change in urinary EVs form HCC (relative to normal control).

Author Contributions

Conceptualization, D.L. and S.J.; methodology, S.W.; validation, D.L. and S.J.; formal analysis, S.W.; investigation, L.H.; resources, L.H.; original draft preparation, D.L.; writing—review and editing, S.J.; visualization, S.W.; supervision, L.H.; project administration, L.H.; funding acquisition, L.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 21675061) and Jilin Province Health Project (No. 2020SCZT057 to SN.J.).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of The Second Hospital of Jilin University.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patient(s) to publish this paper.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

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. Llovet, J.M.; Kelley, R.K.; Villanueva, A.; Singal, A.G.; Pikarsky, E.; Roayaie, S.; Lencioni, R.; Koike, K.; Zucman-Rossi, J.; Finn, R.S. Hepatocellular carcinoma. Nat. Rev. Dis. Prim. 2021, 7, 6. [Google Scholar] [CrossRef] [PubMed]
  3. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Tsuchiya, N.; Sawada, Y.; Endo, I.; Saito, K.; Uemura, Y.; Nakatsura, T. Biomarkers for the early diagnosis of hepatocellular carcinoma. World J. Gastroenterol. 2015, 21, 10573–10583. [Google Scholar] [CrossRef] [PubMed]
  5. Marrero, J.A.; Kulik, L.M.; Sirlin, C.B.; Zhu, A.X.; Finn, R.S.; Abecassis, M.M.; Roberts, L.R.; Heimbach, J.K. Diagnosis, Staging, and Management of Hepatocellular Carcinoma: 2018 Practice Guidance by the American Association for the Study of Liver Diseases. Hepatology 2018, 68, 723–750. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Marrero, J.A.; Lok, A.S. Newer markers for hepatocellular carcinoma. Gastroenterology 2004, 127 (Suppl. 1), S113–S119. [Google Scholar] [CrossRef] [PubMed]
  7. Maluccio, M.; Covey, A. Recent progress in understanding, diagnosing, and treating hepatocellular carcinoma. CA Cancer J. Clin. 2012, 62, 394–399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Babayan, A.; Pantel, K. Advances in liquid biopsy approaches for early detection and monitoring of cancer. Genome Med. 2018, 10, 21. [Google Scholar] [CrossRef] [PubMed]
  9. Kosaka, N.; Kogure, A.; Yamamoto, T.; Urabe, F.; Usuba, W.; Prieto-Vila, M.; Ochiya, T. Exploiting the message from cancer: The diagnostic value of extracellular vesicles for clinical applications. Exp. Mol. Med. 2019, 51, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Krebs, M.G.; Metcalf, R.L.; Carter, L.; Brady, G.; Blackhall, F.H.; Dive, C. Molecular analysis of circulating tumour cells-biology and biomarkers. Nat. Rev. Clin. Oncol. 2014, 11, 129–144. [Google Scholar] [CrossRef]
  11. van Niel, G.; D’Angelo, G.; Raposo, G. Shedding light on the cell biology of extracellular vesicles. Nat. Rev. Mol. Cell Biol. 2018, 19, 213–228. [Google Scholar] [CrossRef] [PubMed]
  12. Pegtel, D.M.; Gould, S.J. Exosomes. Annu. Rev. Biochem. 2019, 88, 487–514. [Google Scholar] [CrossRef] [PubMed]
  13. Colombo, M.; Raposo, G.; Thery, C. Biogenesis, secretion, and intercellular interactions of exosomes and other extracellular vesicles. Annu. Rev. Cell Dev. Biol. 2014, 30, 255–289. [Google Scholar] [CrossRef] [PubMed]
  14. Yuana, Y.; Sturk, A.; Nieuwland, R. Extracellular vesicles in physiological and pathological conditions. Blood Rev. 2013, 27, 31–39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Raposo, G.; Stoorvogel, W. Extracellular vesicles: Exosomes, microvesicles, and friends. J. Cell Biol. 2013, 200, 373–383. [Google Scholar] [CrossRef] [Green Version]
  16. Valadi, H.; Ekstrom, K.; Bossios, A.; Sjostrand, M.; Lee, J.J.; Lotvall, J.O. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat. Cell Biol. 2007, 9, 654–659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Hu, T.; Wolfram, J.; Srivastava, S. Extracellular Vesicles in Cancer Detection: Hopes and Hypes. Trends Cancer 2021, 7, 122–133. [Google Scholar] [CrossRef] [PubMed]
  18. Jiang, Y.F.; Wei, S.N.; Geng, N.; Qin, W.W.; He, X.; Wang, X.H.; Qi, Y.P.; Song, S.; Wang, P. Evaluation of circulating small extracellular vesicle-derived miRNAs as diagnostic biomarkers for differentiating between different pathological types of early lung cancer. Sci. Rep. 2022, 12, 17201. [Google Scholar] [CrossRef]
  19. Herreros-Villanueva, M.; Bujanda, L. Glypican-1 in exosomes as biomarker for early detection of pancreatic cancer. Ann. Transl. Med. 2016, 4, 64. [Google Scholar]
  20. Costa-Silva, B.; Aiello, N.M.; Ocean, A.J.; Singh, S.; Zhang, H.; Thakur, B.K.; Becker, A.; Hoshino, A.; Mark, M.T.; Molina, H.; et al. Pancreatic cancer exosomes initiate pre-metastatic niche formation in the liver. Nat. Cell Biol. 2015, 17, 816–826. [Google Scholar] [CrossRef]
  21. Lin, J.; Lin, W.; Bai, Y.; Liao, Y.; Lin, Q.; Chen, L.; Wu, Y. Identification of exosomal hsa-miR-483-5p as a potential biomarker for hepatocellular carcinoma via microRNA expression profiling of tumor-derived exosomes. Exp. Cell Res. 2022, 417, 113232. [Google Scholar] [CrossRef] [PubMed]
  22. Tan, N.Y.; Bailey, U.M.; Jamaluddin, M.F.; Mahmud, S.H.; Raman, S.C.; Schulz, B.L. Sequence-based protein stabilization in the absence of glycosylation. Nat. Commun. 2014, 5, 3099. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Croci, D.O.; Cerliani, J.P.; Dalotto-Moreno, T.; Mendez-Huergo, S.P.; Mascanfroni, I.D.; Dergan-Dylon, S.; Toscano, M.A.; Caramelo, J.J.; Garcia-Vallejo, J.J.; Ouyang, J.; et al. Glycosylation-dependent lectin-receptor interactions preserve angiogenesis in anti-VEGF refractory tumors. Cell 2014, 156, 744–758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kolarich, D.; Lepenies, B.; Seeberger, P.H. Glycomics, glycoproteomics and the immune system. Curr. Opin. Chem. Biol. 2012, 16, 214–220. [Google Scholar] [CrossRef] [PubMed]
  25. Badr, H.A.; Alsadek, D.M.; Darwish, A.A.; Elsayed, A.I.; Bekmanov, B.O.; Khussainova, E.M.; Zhang, X.; Cho, W.C.; Djansugurova, L.B.; Li, C.Z. Lectin approaches for glycoproteomics in FDA-approved cancer biomarkers. Expert Rev. Proteom. 2014, 11, 227–236. [Google Scholar] [CrossRef] [PubMed]
  26. Zhu, J.; Lin, Z.; Wu, J.; Yin, H.; Dai, J.; Feng, Z.; Marrero, J.; Lubman, D.M. Analysis of serum haptoglobin fucosylation in hepatocellular carcinoma and liver cirrhosis of different etiologies. J. Proteome Res. 2014, 13, 2986–2997. [Google Scholar] [CrossRef] [PubMed]
  27. Choi, J.; Park, Y.; Kim, J.H.; Kim, H.S. Evaluation of revisited fucosylated alpha-fetoprotein (AFP-L3) with an autoanalyzer muTAS in a clinical laboratory. Clin. Chim. Acta 2012, 413, 170–174. [Google Scholar] [CrossRef] [PubMed]
  28. Costa, J. Glycoconjugates from extracellular vesicles: Structures, functions and emerging potential as cancer biomarkers. Biochim. Biophys. Acta Rev. Cancer 2017, 1868, 157–166. [Google Scholar] [CrossRef] [PubMed]
  29. Martins, A.M.; Ramos, C.C.; Freitas, D.; Reis, C.A. Glycosylation of Cancer Extracellular Vesicles: Capture Strategies, Functional Roles and Potential Clinical Applications. Cells 2021, 10, 109. [Google Scholar] [CrossRef]
  30. Melo, S.A.; Luecke, L.B.; Kahlert, C.; Fernandez, A.F.; Gammon, S.T.; Kaye, J.; LeBleu, V.S.; Mittendorf, E.A.; Weitz, J.; Rahbari, N.; et al. Glypican-1 identifies cancer exosomes and detects early pancreatic cancer. Nature 2015, 523, 177–182. [Google Scholar] [CrossRef] [Green Version]
  31. Sakaue, T.; Koga, H.; Iwamoto, H.; Nakamura, T.; Ikezono, Y.; Abe, M.; Wada, F.; Masuda, A.; Tanaka, T.; Fukahori, M.; et al. Glycosylation of ascites-derived exosomal CD133: A potential prognostic biomarker in patients with advanced pancreatic cancer. Med. Mol. Morphol. 2019, 52, 198–208. [Google Scholar] [CrossRef] [PubMed]
  32. Hoshino, A.; Kim, H.S.; Bojmar, L.; Gyan, K.E.; Cioffi, M.; Hernandez, J.; Zambirinis, C.P.; Rodrigues, G.; Molina, H.; Heissel, S.; et al. Extracellular Vesicle and Particle Biomarkers Define Multiple Human Cancers. Cell 2020, 182, 1044–1061.e18. [Google Scholar] [CrossRef]
  33. Menck, K.; Scharf, C.; Bleckmann, A.; Dyck, L.; Rost, U.; Wenzel, D.; Dhople, V.M.; Siam, L.; Pukrop, T.; Binder, C.; et al. Tumor-derived microvesicles mediate human breast cancer invasion through differentially glycosylated EMMPRIN. J. Mol. Cell Biol. 2015, 7, 143–153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Gomes, J.; Gomes-Alves, P.; Carvalho, S.B.; Peixoto, C.; Alves, P.M.; Altevogt, P.; Costa, J. Extracellular Vesicles from Ovarian Carcinoma Cells Display Specific Glycosignatures. Biomolecules 2015, 5, 1741–1761. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Surman, M.; Hoja-Lukowicz, D.; Szwed, S.; Kedracka-Krok, S.; Jankowska, U.; Kurtyka, M.; Drozdz, A.; Litynska, A.; Stepien, E.; Przybylo, M. An Insight into the Proteome of Uveal Melanoma-Derived Ectosomes Reveals the Presence of Potentially Useful Biomarkers. Int. J. Mol. Sci. 2019, 20, 3789. [Google Scholar] [CrossRef] [Green Version]
  36. Costa, J.; Gatermann, M.; Nimtz, M.; Kandzia, S.; Glatzel, M.; Conradt, H.S. N-Glycosylation of Extracellular Vesicles from HEK-293 and Glioma Cell Lines. Anal. Chem. 2018, 90, 7871–7879. [Google Scholar] [CrossRef] [PubMed]
  37. Lv, J.; Wang, Z.; Li, F.; Zhang, Y.; Lu, H. Reverse capture for selectively and sensitively revealing the N-glycome of serum exosomes. Chem. Commun. 2019, 55, 14339–14342. [Google Scholar] [CrossRef]
  38. Wang, L.; Chen, X.; Wang, L.; Wang, S.; Li, W.; Liu, Y.; Zhang, J. Knockdown of ST6Gal-I expression in human hepatocellular carcinoma cells inhibits their exosome-mediated proliferation- and migration-promoting effects. IUBMB Life 2021, 73, 1378–1391. [Google Scholar] [CrossRef]
  39. Wu, X.; Li, L.; Iliuk, A.; Tao, W.A. Highly Efficient Phosphoproteome Capture and Analysis from Urinary Extracellular Vesicles. J. Proteome Res. 2018, 17, 3308–3316. [Google Scholar] [CrossRef] [PubMed]
  40. Xiao, K.; Tian, Z. GPSeeker Enables Quantitative Structural N-Glycoproteomics for Site- and Structure-Specific Characterization of Differentially Expressed N-Glycosylation in Hepatocellular Carcinoma. J. Proteome Res. 2019, 18, 2885–2895. [Google Scholar] [CrossRef]
  41. Kalra, H.; Simpson, R.J.; Ji, H.; Aikawa, E.; Altevogt, P.; Askenase, P.; Bond, V.C.; Borras, F.E.; Breakefield, X.; Budnik, V.; et al. Vesiclepedia: A compendium for extracellular vesicles with continuous community annotation. PLoS Biol. 2012, 10, e1001450. [Google Scholar] [CrossRef] [PubMed]
  42. Iacobelli, S.; Arno, E.; D’Orazio, A.; Coletti, G. Detection of antigens recognized by a novel monoclonal antibody in tissue and serum from patients with breast cancer. Cancer Res. 1986, 46, 3005–3010. [Google Scholar] [PubMed]
  43. Linsley, P.S.; Horn, D.; Marquardt, H.; Brown, J.P.; Hellstrom, I.; Hellstrom, K.E.; Ochs, V.; Tolentino, E. Identification of a novel serum protein secreted by lung carcinoma cells. Biochemistry 1986, 25, 2978–2986. [Google Scholar] [CrossRef] [PubMed]
  44. Jiang, S.S.; Weng, D.S.; Wang, Q.J.; Pan, K.; Zhang, Y.J.; Li, Y.Q.; Li, J.J.; Zhao, J.J.; He, J.; Lv, L.; et al. Galectin-3 is associated with a poor prognosis in primary hepatocellular carcinoma. J. Transl. Med. 2014, 12, 273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Correale, M.; Giannuzzi, V.; Iacovazzi, P.A.; Valenza, M.A.; Lanzillotta, S.; Abbate, I.; Quaranta, M.; Caruso, M.L.; Elba, S.; Manghisi, O.G. Serum 90K/MAC-2BP glycoprotein levels in hepatocellular carcinoma and cirrhosis. Anticancer Res. 1999, 19, 3469–3472. [Google Scholar] [PubMed]
  46. Arbelaiz, A.; Azkargorta, M.; Krawczyk, M.; Santos-Laso, A.; Lapitz, A.; Perugorria, M.J.; Erice, O.; Gonzalez, E.; Jimenez-Aguero, R.; Lacasta, A.; et al. Serum extracellular vesicles contain protein biomarkers for primary sclerosing cholangitis and cholangiocarcinoma. Hepatology 2017, 66, 1125–1143. [Google Scholar] [CrossRef] [Green Version]
  47. Zhao, L.; Shi, J.; Chang, L.; Wang, Y.; Liu, S.; Li, Y.; Zhang, T.; Zuo, T.; Fu, B.; Wang, G.; et al. Serum-Derived Exosomal Proteins as Potential Candidate Biomarkers for Hepatocellular Carcinoma. ACS Omega 2021, 6, 827–835. [Google Scholar] [CrossRef]
  48. Zhang, Y.; Zhang, J.; Chen, X.; Yang, Z. Polymeric immunoglobulin receptor (PIGR) exerts oncogenic functions via activating ribosome pathway in hepatocellular carcinoma. Int. J. Med. Sci. 2021, 18, 364–371. [Google Scholar] [CrossRef] [PubMed]
  49. Ai, J.; Tang, Q.; Wu, Y.; Xu, Y.; Feng, T.; Zhou, R.; Chen, Y.; Gao, X.; Zhu, Q.; Yue, X.; et al. The role of polymeric immunoglobulin receptor in inflammation-induced tumor metastasis of human hepatocellular carcinoma. J. Natl. Cancer Inst. 2011, 103, 1696–1712. [Google Scholar] [CrossRef] [Green Version]
  50. Tey, S.K.; Wong, S.W.K.; Tung Chan, J.Y.; Mao, X.; Ng, T.H.; Yeung, C.L.S.; Leung, Z.; Fung, H.L.; Tang, A.H.N.; Wong, D.K.H.; et al. Patient pIgR-enriched extracellular vesicles drive cancer stemness, tumorigenesis and metastasis in hepatocellular carcinoma. J. Hepatol. 2022, 76, 883–895. [Google Scholar] [CrossRef]
  51. Sainz, I.M.; Pixley, R.A.; Colman, R.W. Fifty years of research on the plasma kallikrein-kinin system: From protein structure and function to cell biology and in-vivo pathophysiology. Thromb. Haemost. 2007, 98, 77–83. [Google Scholar] [CrossRef] [PubMed]
  52. Abdel Wahab, A.H.A.; El-Halawany, M.S.; Emam, A.A.; Elfiky, A.; Abd Elmageed, Z.Y. Identification of circulating protein biomarkers in patients with hepatocellular carcinoma concomitantly infected with chronic hepatitis C virus. Biomarkers 2017, 22, 621–628. [Google Scholar] [CrossRef]
  53. Wang, J.; Wang, X.; Lin, S.; Chen, C.; Wang, C.; Ma, Q.; Jiang, B. Identification of kininogen-1 as a serum biomarker for the early detection of advanced colorectal adenoma and colorectal cancer. PLoS ONE 2013, 8, e70519. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wang, M.; Sanda, M.; Comunale, M.A.; Herrera, H.; Swindell, C.; Kono, Y.; Singal, A.G.; Marrero, J.; Block, T.; Goldman, R.; et al. Changes in the Glycosylation of Kininogen and the Development of a Kininogen-Based Algorithm for the Early Detection of HCC. Cancer Epidemiol. Biomark. Prev. 2017, 26, 795–803. [Google Scholar] [CrossRef] [Green Version]
  55. Trigiante, G.; Lu, X. ASPP [corrected] and cancer. Nat. Rev. Cancer 2006, 6, 217–226. [Google Scholar] [CrossRef] [PubMed]
  56. Sottocornola, R.; Royer, C.; Vives, V.; Tordella, L.; Zhong, S.; Wang, Y.; Ratnayaka, I.; Shipman, M.; Cheung, A.; Gaston-Massuet, C.; et al. ASPP2 binds Par-3 and controls the polarity and proliferation of neural progenitors during CNS development. Dev. Cell 2010, 19, 126–137. [Google Scholar] [CrossRef] [Green Version]
  57. Sullivan, A.; Lu, X. ASPP: A new family of oncogenes and tumour suppressor genes. Br. J. Cancer 2007, 96, 196–200. [Google Scholar] [CrossRef] [Green Version]
  58. Zhao, J.; Wu, G.; Bu, F.; Lu, B.; Liang, A.; Cao, L.; Tong, X.; Lu, X.; Wu, M.; Guo, Y. Epigenetic silence of ankyrin-repeat-containing, SH3-domain-containing, and proline-rich-region- containing protein 1 (ASPP1) and ASPP2 genes promotes tumor growth in hepatitis B virus-positive hepatocellular carcinoma. Hepatology 2010, 51, 142–153. [Google Scholar] [CrossRef]
  59. Chen, R.; Wang, H.; Liang, B.; Liu, G.; Tang, M.; Jia, R.; Fan, X.; Jing, W.; Zhou, X.; Wang, H.; et al. Downregulation of ASPP2 improves hepatocellular carcinoma cells survival via promoting BECN1-dependent autophagy initiation. Cell Death Dis. 2016, 7, e2512. [Google Scholar] [CrossRef]
  60. Wang, Y.; Bu, F.; Royer, C.; Serres, S.; Larkin, J.R.; Soto, M.S.; Sibson, N.R.; Salter, V.; Fritzsche, F.; Turnquist, C.; et al. ASPP2 controls epithelial plasticity and inhibits metastasis through beta-catenin-dependent regulation of ZEB1. Nat. Cell Biol. 2014, 16, 1092–1104. [Google Scholar] [CrossRef]
  61. Li, S.; Li, H.; Yang, X.; Wang, W.; Huang, A.; Li, J.; Qin, X.; Li, F.; Lu, G.; Ding, H.; et al. Vasorin is a potential serum biomarker and drug target of hepatocarcinoma screened by subtractive-EMSA-SELEX to clinic patient serum. Oncotarget 2015, 6, 10045–10059. [Google Scholar] [CrossRef]
  62. Huang, A.; Dong, J.; Li, S.; Wang, C.; Ding, H.; Li, H.; Su, X.; Ge, X.; Sun, L.; Bai, C.; et al. Exosomal transfer of vasorin expressed in hepatocellular carcinoma cells promotes migration of human umbilical vein endothelial cells. Int. J. Biol. Sci. 2015, 11, 961–969. [Google Scholar] [CrossRef] [Green Version]
  63. Andaluz Aguilar, H.; Iliuk, A.B.; Chen, I.H.; Tao, W.A. Sequential phosphoproteomics and N-glycoproteomics of plasma-derived extracellular vesicles. Nat. Protoc. 2020, 15, 161–180. [Google Scholar] [CrossRef] [PubMed]
  64. Boersema, P.J.; Raijmakers, R.; Lemeer, S.; Mohammed, S.; Heck, A.J. Multiplex peptide stable isotope dimethyl labeling for quantitative proteomics. Nat. Protoc. 2009, 4, 484–494. [Google Scholar] [CrossRef] [PubMed]
  65. Yu, L.; Li, X.; Guo, Z.; Zhang, X.; Liang, X. Hydrophilic interaction chromatography based enrichment of glycopeptides by using click maltose: A matrix with high selectivity and glycosylation heterogeneity coverage. Chemistry 2009, 15, 12618–12626. [Google Scholar] [CrossRef]
  66. Xiao, K.; Wang, Y.; Shen, Y.; Han, Y.; Tian, Z. Large-scale identification and visualization of N-glycans with primary structures using GlySeeker. Rapid Commun. Mass Spectrom. 2018, 32, 142–148. [Google Scholar] [CrossRef]
  67. Pathan, M.; Keerthikumar, S.; Chisanga, D.; Alessandro, R.; Ang, C.S.; Askenase, P.; Batagov, A.O.; Benito-Martin, A.; Camussi, G.; Clayton, A.; et al. A novel community driven software for functional enrichment analysis of extracellular vesicles data. J. Extracell. Vesicles 2017, 6, 1321455. [Google Scholar] [CrossRef] [Green Version]
  68. Huang da, W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Workflow for the proteomic and glycoproteomic analysis of urinary EVs derived from patients with hepatocellular carcinoma and normal controls.
Figure 1. Workflow for the proteomic and glycoproteomic analysis of urinary EVs derived from patients with hepatocellular carcinoma and normal controls.
Molecules 28 01293 g001
Figure 2. Venn diagram depicted the overlap of proteins identified from HCC-derived urinary EVs with those published in the Vesiclepedia database (a) and top 100 EV proteins from Vesiclepedia (b). (c) GO cellular component annotations of the identified proteins in HCC-derived urinary EVs.
Figure 2. Venn diagram depicted the overlap of proteins identified from HCC-derived urinary EVs with those published in the Vesiclepedia database (a) and top 100 EV proteins from Vesiclepedia (b). (c) GO cellular component annotations of the identified proteins in HCC-derived urinary EVs.
Molecules 28 01293 g002
Figure 3. Identification and distribution of intact N-glycopeptides in urinary EVs from patients with HCC relative to NC. (a)The identified intact N-glycopeptide IDs (FDR ≤ 1%). (b) Sequence motifs for N-glycosites having either the N-XS or N-X-T sequon and their relative percentage in the unique Nglycosites identified. (c) The percentage of each N-glycan type in urinary EVs from patients with HCC and NC. (d) Distribution of glycosites (left) and glycans per glycosite (right) of the identified glycoproteins.
Figure 3. Identification and distribution of intact N-glycopeptides in urinary EVs from patients with HCC relative to NC. (a)The identified intact N-glycopeptide IDs (FDR ≤ 1%). (b) Sequence motifs for N-glycosites having either the N-XS or N-X-T sequon and their relative percentage in the unique Nglycosites identified. (c) The percentage of each N-glycan type in urinary EVs from patients with HCC and NC. (d) Distribution of glycosites (left) and glycans per glycosite (right) of the identified glycoproteins.
Molecules 28 01293 g003
Figure 4. Up-regulation of intact N-glycopeptide GLNLTEDTYKPR_N2H3F1S0 in urinary EVs from patients with HCC relative to NC; the N-glycosite is N398 on Galectin-3-binding protein (LG3BP_HUMAN, Q08380). (a) Paired precursor ions. (b,c) N-glycan and peptide backbone graphical fragmentation maps annotated with the matched fragment ions. (d) The MS/MS spectrum with the matched fragment ions.
Figure 4. Up-regulation of intact N-glycopeptide GLNLTEDTYKPR_N2H3F1S0 in urinary EVs from patients with HCC relative to NC; the N-glycosite is N398 on Galectin-3-binding protein (LG3BP_HUMAN, Q08380). (a) Paired precursor ions. (b,c) N-glycan and peptide backbone graphical fragmentation maps annotated with the matched fragment ions. (d) The MS/MS spectrum with the matched fragment ions.
Molecules 28 01293 g004
Figure 5. Down-regulation of intact N-glycopeptide KNQSSEDILR_N4H5F1S1 in urinary EVs from patients with HCC relative to NC; the N-glycosite is N478 on Apoptosis-stimulating of p53 protein 2 (ASPP2_HUMAN, Q13625). (a) Paired precursor ions. (b,c) N-glycan and peptide backbone graphical fragmentation maps annotated with the matched fragment ions. (d) The MS/MS spectrum with the matched fragment ions.
Figure 5. Down-regulation of intact N-glycopeptide KNQSSEDILR_N4H5F1S1 in urinary EVs from patients with HCC relative to NC; the N-glycosite is N478 on Apoptosis-stimulating of p53 protein 2 (ASPP2_HUMAN, Q13625). (a) Paired precursor ions. (b,c) N-glycan and peptide backbone graphical fragmentation maps annotated with the matched fragment ions. (d) The MS/MS spectrum with the matched fragment ions.
Molecules 28 01293 g005
Figure 6. Simultaneous up- and down-regulation on alpha-2-macroglobulin-like protein 1 (A2ML1_HUMAN, A8K2U0) and vasorin (VASN_HUMAN, Q6EMK4) with different intact N-glycopeptides.
Figure 6. Simultaneous up- and down-regulation on alpha-2-macroglobulin-like protein 1 (A2ML1_HUMAN, A8K2U0) and vasorin (VASN_HUMAN, Q6EMK4) with different intact N-glycopeptides.
Molecules 28 01293 g006
Figure 7. Gene ontology analysis of the intact N-glycoproteins corresponding to the differentially expressed intact N-glycopeptides in HCC relative to NC. (a) Cellular component. (b) Molecular function. (c) Biological process.
Figure 7. Gene ontology analysis of the intact N-glycoproteins corresponding to the differentially expressed intact N-glycopeptides in HCC relative to NC. (a) Cellular component. (b) Molecular function. (c) Biological process.
Molecules 28 01293 g007
Table 1. The intact N-glycopeptides identified from probable serine carboxypeptidase CPVL (Q9H3G5, CPVL_HUMAN) modified with N2H5F0S0 on N-glycosites N81, N307 and N346 and N2H4F0S0 and N2H6F0S0 on the N-glycosite N81 and N346 as well as N2H3F0S0 and N2H4F1S0 on N-glycosite 346.
Table 1. The intact N-glycopeptides identified from probable serine carboxypeptidase CPVL (Q9H3G5, CPVL_HUMAN) modified with N2H5F0S0 on N-glycosites N81, N307 and N346 and N2H4F0S0 and N2H6F0S0 on the N-glycosite N81 and N346 as well as N2H3F0S0 and N2H4F1S0 on N-glycosite 346.
GlycositePeptide SequenceCompositionGlycan-LinkageGF ScoreStructural Diagnostic Ions
81SYAGFLTVNKN2H4F0S001Y41Y41M(31M21M)61M9BII2,BI2,YI1,YI2,BI3,YII3,YI3,YII4,MH,YI1
N2H5F0S001Y41Y41M(31M)61M(31M)61M1124AI4,ZI1,YI1,YI2,YI3,BI4,
YI4,YII4,YIII3,MH,02XI1,ZI1,YI1
N2H6F0S001Y41Y41M(31M21M)61M(31M)61M3BI2,YI3,YI3
307LLDGDLTSDPSYFQNVTGCSNYYNFLRN2H5F0S001Y41Y41M(31M)61M(31M)61M3BI2,BI4,YI1
346QAIHVGNQTFNDGTIVEKN2H3F0S001Y41Y41M(31M)61M5BI3,YI1,YI2,ZI1,YI1
N2H4F0S001Y41Y41M(31M21M)61M12BII2,BI2,YI1,YI2,BI3,YII3,YI3,YII4,MH,02XI1,ZI1,YI1,YI2
N2H4F1S001Y(61F)41Y41M(31M)61M61M4BI4,YI3,YI4,YII3,YI1
N2H5F0S001Y41Y41M(31M)61M(31M)61M12ZI1,YI1,YI2,YI3,BI4,YI4,YII4,YIII3,02XI1,MH,ZI1,YI1,YI2,YI3
N2H6F0S001Y41Y41M(31M)61M(31M21M)61M2BI2,YI3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, D.; Jia, S.; Wang, S.; Hu, L. Glycoproteomic Analysis of Urinary Extracellular Vesicles for Biomarkers of Hepatocellular Carcinoma. Molecules 2023, 28, 1293. https://doi.org/10.3390/molecules28031293

AMA Style

Li D, Jia S, Wang S, Hu L. Glycoproteomic Analysis of Urinary Extracellular Vesicles for Biomarkers of Hepatocellular Carcinoma. Molecules. 2023; 28(3):1293. https://doi.org/10.3390/molecules28031293

Chicago/Turabian Style

Li, Dejun, Shengnan Jia, Shuyue Wang, and Lianghai Hu. 2023. "Glycoproteomic Analysis of Urinary Extracellular Vesicles for Biomarkers of Hepatocellular Carcinoma" Molecules 28, no. 3: 1293. https://doi.org/10.3390/molecules28031293

APA Style

Li, D., Jia, S., Wang, S., & Hu, L. (2023). Glycoproteomic Analysis of Urinary Extracellular Vesicles for Biomarkers of Hepatocellular Carcinoma. Molecules, 28(3), 1293. https://doi.org/10.3390/molecules28031293

Article Metrics

Back to TopTop