BioMed Research International
Volume 2019 (2019), Article ID 2871289, 11 pages
https://doi.org/10.1155/2019/2871289

Research Article

Glycopatterns and Glycoproteins Changes in MCN and SCN: A Prospective Cohort Study

Ying Wang,¹ Yufa Sun,² Jia Feng,³ Zheng Li,⁴ Hanjie Yu,⁴ Xiang Ding,⁵ Fuquan Yang,⁵ and Enqiang Linghu⁶

¹Department of Gastroenterology, The Affiliated Fu Xing Hospital of Capital Medical University, Beijing 100038, China
²Department of Health Care, Central Guard Bureau, Beijing 100034, China
³Department of Gastroenterology, Bethune International Peace Hospital, Shi Jia Zhuang 050082, China
⁴Laboratory for Functional Glycomics, College of Life Sciences, Northwest University, Xi’an 710069, China
⁵Laboratory of Proteomics, Institute of Biophysics, Chinese Academy of Sciences, Beijing 100101, China
⁶Department of Gastroenterology, Chinese PLA General Hospital, Beijing 100853, China

Correspondence should be addressed to Ying Wang; 1073803334@qq.com and Enqiang Linghu; linghuenqiang@vip.sina.com

Received 23 March 2019; Revised 29 June 2019; Accepted 8 July 2019; Published 31 July 2019

Academic Editor: Toshimi Chiba

Copyright © 2019 Ying Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. Advances in imaging improve the detection of malignant pancreatic cystic including mucinous cystic neoplasm (MCN), intraductal papillary mucinous neoplasm (IPMN), and mucinous cystic adenocarcinoma (MCA), but the distinction between benign and malignant lesions remains a problem. In an effort to establish glycopatterns as potential biomarkers for differential diagnosis between MCN and SCN, we systematically investigated the alterations of glycopatterns in cystic fluids for both SCN and MCN. Methods. Among the 75 patients enrolled, 37 were diagnosed as MCN and 38 as SCN based on histology. Lectin microarray analysis was performed on each sample, and the fluorescence intensity was used to obtain the fold-change. Then, mixed cyst fluids of MCN group and SCN group were cross bonded with magnetic particles coupled by Lectin STL and WGA, respectively. Hydrophilic interaction liquid chromatography (HILIC) enrichment was performed, liquid chromatography (LC)/mass spectrometry (MS) analysis and bioinformatical analysis was conducted to find the differential glycoproteins between MCNs and SCNs. Results. Through analysis of lectin microarray between MCNs and SCNs, stronger lectin signal patterns were assigned to Lectin WFA, DBA, STL, WGA, and BPL; and weaker signal patterns were assigned to Lectin PTL-I, Con A, ACA, and MAL-I. The glycoproteins were enriched by STL or WGA-coupled magnetic particles. Furthermore, the 10 identified correspondding genes were found to be significantly elevated in the mucinous cystadenoma: CLU, A2M, FGA, FGB, FGG, PLG, SERPINA1, SERPING1, C5, C8A, and C9. Bioinformatics analysis revealed that the above genes may activate the KEGG pathway: immune complement system. Conclusion. This study shows changes in glycopatterns and glycoproteins are associated with MCNs and SCNs.

1. Introduction

Based on imaging studies, pancreatic cystic neoplasms (PCNs) are being diagnosed with increasing frequency, the prevalence of PCNs in the general population is estimated to be between 2.6% and 45% [1–4]. Clinical differential diagnosis between the potential malignant pancreatic mucinous cystic neoplasm (MCN) and the mostly benign serous cystic neoplasm (SCN) remains challenging [5–8]. According to the Cochrane evidence-based guideline on PCNs, there are no clinical biomarkers of DNA, RNA, or protein in the blood to differentiate PCNs, high-grade dysplasia, or adenocarcinoma [9–12]. According to the evidence, carcinoembryonic antigen (CEA) >192ng/L in cyst fluids could be reports on potential malignant PCNs. Mostly, gene mutations of KRAS at codon 12,13 are hallmark genetic alterations of pancreatic ductal adenocarcinoma (PDAC), and GNAS codon 201 mutations are hallmark genetic alterations of IPMNs [13, 14].

Back in 2004, Hägglund first applied hydrophilic interaction liquid chromatography (HILIC) enrichment method to glycopeptide separation [15]. Nowadays, HILIC is a separation and analysis method based on positive phase chromatography with water/organic phase as mobile phase and polar material as fixed phase [16, 17]. With its developments over years, HILIC has become an important method mainly used to enrich glycosylated peptides.

Aberrant protein glycosylation often occurs during malignant transformation and leads to the expression of specific tumour-associated glycans [18–20]. Therefore, changes in glycosylation of proteins have the tremendous potential to be used as biomarkers to help the diagnosis of malignant diseases [21–23]. We report 75 cases, with 37 in the MCN group and 38 in the SCN group, using lectin microarrays to quantitatively analyse the differential types of glycan patterns and using the magnetic particles of labelled lectins to enrich the glycoproteins. Then, trypsin was used to digest proteins to peptides, and HILIC enrichment was performed. Equal sample volumes were dissolved and loaded for 3 times, and then liquid chromatography (LC) and mass spectrometry (MS) were used to identify the protein. The aim of this paper was to study the differential expression profiles of glycan-related proteins in MCNs using lectin microarrays. The overall strategy is summarized in Figure 1.

Figure 1: The flow chart of differential lectins enriched the glycoproteins.

2. Materials and Methods

2.1. Study Design

This prospective study was approved by Fu Xing Hospital, Capital Medical University ethics committee (Approval Notice no. 2017FXHEC—KY044). The clinical records, EUS images, pathology, and surgical reports in this study were confirmed to be accurate.

2.2. Selection Flow Chart

Form April 2015 to May 2018, 211 patients with PCNs enrolled in this study. The selection strategy for the patients is described below. The first step was diagnosis according to clinical presentation and abdominal CT/MRI. The second step was whether it was feasible to perform the endoscopic ultrasonography-guided fine needle aspiration (EUS-FNA), and if it was, a 19G/22G needle was used (Cook Co., Boston, USA) along with a GF-UE260-AL5/GF-UM200 endoscope (Olympus Co., Tokyo, Japan) and an EU-ME2 endoscopic ultrasonography processor (Olympus Co., Tokyo, Japan). Needle puncture of the patient was used to obtain the fluid samples. The third step was the collection of cystic fluid for biochemical analysis, molecular diagnosis, and cytological examination. All patients with the following contraindications were excluded:

(i)pancreatic pseudocysts;(ii)severe acute pancreatitis;(iii)malignant tumours;(iv)severe cardiopulmonary circulatory system diseases;(v)blood coagulation disorders.

Every patient provided informed consent. For those who did not undergo surgery, follow-up was conducted every 3 months, during which the EUS was reviewed to determine whether the cyst had become too large for EUS-FNA. Patients with concern of neoplastic growth underwent surgery but did not receive additional intervention.

2.3. Tumour Puncture and Extraction of Cystic Fluid

The patients were given intravenous anaesthesia by aspiration with oropharyngeal intubation. The most suitable puncture route was chosen so that important structures, such as the abdominal viscera and blood vessels, would be avoided in EUS-FNA. The lesions were punctured and connected to the negative pressure syringe, after which the cystic fluid was extracted at constant negative pressure. The extracellular cystic fluid was divided into two 1.0 mL cryopreservation tubes by a specialist, and a protease inhibitor was added to each tube at a ratio of 1:100. After the solution was mixed gently by vortexing, it was stored in a freezer at −80°C at the biological specimen bank. The frozen tubes were removed from storage during each experiment, which was performed after dissolution.

2.4. Lectin Microarray Analysis

Lectin microarrays were constructed and analysed by 37 commercial lectins from Vector Laboratories (Burlingame, CA, USA), Sigma-Aldrich (St. Louis, MO, USA), or Merck (Darmstadt, Germany). Cystic fluid protein samples were labelled with the fluorescent dye Cy3 (GE Healthcare; Buckinghamshire, UK) and were applied to the lectin microarrays. The slides were incubated in a humidity-controlled incubator at 50% humidity overnight to allow lectin immobilization. After incubation, the slides were blocked with blocking buffer (50 mM ethanolamine, 2% (w/v) BSA and 500 mM glycine in 50 mM sodium borate buffer, pH 8) for one hour and were then rinsed three times with PBST (0.05% Tween 20 in 0.1 M phosphate buffer containing 0.15 M NaCl, pH 7.4). This was followed by a final rinse in PBS (0.1 M phosphate buffer containing 0.15 M NaCl, pH 7.4). Before use, the slides were dried by centrifugation at 600 rpm for 5 minutes. After incubation, the slides were scanned with a GenePix 4000B confocal scanner (Axon Instruments; Union City, CA, USA).

2.5. Isolation and Purification of Glycoproteins

The cystic fluids from MCN group and SCN group were mixed in two main types and diluted in binding solution to a final concentration of 1 mg/mL and were mixed with 300 μl of STL- or WGA-coupled magnetic particles, which were washed in binding solution. The mixtures were incubated for one hour at 25°C with shaking. Next, the magnetic particle conjugates were washed 5 times in washing solution to remove nonspecific proteins. Subsequently, the conjugates were incubated with elution buffer for one hour, and the glycoproteins enriched with STL or WGA were collected and quantitated by BCA assay.

2.6. Enzyme Cutting and Desalting

Before an appropriate amount of peptide sample was added, Lys-C hydrolytic enzyme was added in proportion to enzyme/protein at a ratio of 1:50; enzymolysis then occurred at 37°C for 6 hours. The sample was diluted with 50 mM NH4HCO3 solution so that the concentration of urea in the sample was lower than 1 M. Trypsin was added at an enzyme/protein ratio of 1:50, and the sample was incubated at 37°C in the enzyme solution overnight. Waters Oasis HLB solid-phase extraction (SPE) columns and an appropriate amount of 80% acetonitrile (ACN)/0.5% formic acid (FA) activated solid-phase extraction; this step was repeated twice. The FA balanced solid-phase extraction column step was also repeated twice. The enzymatic hydrolysis peptide sample flowed through the solid-phase extraction column, which was adsorbed on the filter. Moreover, the filtrate was collected, and passage through the column was repeated three times. Then, 0.5% FA was applied three times to clean the column. After washing the solid-phase extraction columns with an appropriate amount of 20% ACN/0.5% FA, 60% ACN/0.5% FA, and 80% ACN/0.5% FA, the three eluents that were collected were combined.

2.7. LC-MS/MS Analysis

LC-MS/MS analysis was performed on an EASY-nLC 1000 HPLC coupled with a Q Exactive™ mass spectrometer (Thermo Scientific, San Jose, CA, USA). Lyophilized peptide samples were resuspended in 0.1% FA and then separated on a reversed-phase C18 column in-house packed analytical column (75 μm ID×20 cm, 3 μm, ReproSil-Pur C18 AQ, Dr. Maisch GmbH, Germany) using mobile phase A (0.1% FA in water) and mobile phase B (0.1% FA in acetonitrile) with a gradient (5-8% B, 8 min; 8-22% B, 50 min; 22-32% B, 12 min; 32-95% B, 1 min; 95% B, 7 min) at a flow rate of 280 nL/min. The mass spectrometer was operated in positive ion mode and in the data-dependent acquisition mode. Full-scan MS spectra (from m/z 300 to 1600) were acquired in the Orbitrap at a high resolution of 70,000 (m/z 200) with an automatic gain control (AGC) of 3 × 10⁶ and a maximum injection time of 60 ms. The top 20 precursor ions were selected from each MS full scan with an isolation width of 2 m/z in the HCD collision cell at a normalized collision energy of 27%. Subsequently, MS/MS spectra were acquired in the Orbitrap with a resolution of 17,500 at m/z 200. The target value was 5 × 10⁴ with a maximum injection time of 80 milliseconds. Ions selected for MS/MS were dynamically excluded for a duration of 40 seconds. For nanoelectrospray ion source setting, the spray voltage was 2.0 kV, no sheath gas flow used, and the heated capillary temperature was 320°C. All raw data were viewed in Thermo Xcalibur version 2.2.

2.8. Data Analysis

The original data of microarrays need to be normalized by median normalization method for minimizing the possible systematic variation. The average background was subtracted, each sample was printed in triplicate, and data points with value outside of the average background ± standard deviation (SD) were removed from each data point. [24]. The normalized data of the MCN and SCN groups were compared with each other based on fold-change, according to the following criteria: p < 0.05 fold-change >1.45 or < 0.69; in pairs indicated upregulation or downregulation, respectively. The data were analysed by t-test using SPSS Version 22.0 (IBM, New York, USA).

Each enriched glycoprotein sample and the original mass spectrometry data were retrieved by the Proteome Discoverer™ Software. After the protein was identified, the intensity value was given, the protein intensity value was obtained for all three identifications, and the median value was extracted three times. The final protein intensity was expressed as the mean value of the three repetitions plus or minus the standard deviation value. The ratio of the normalized protein value in the MCN group to that in the SCN group was calculated to compare the relative changes in protein glycosylation. The fold-change > 1.21 as upregulated, and the fold-change <0.83 as downregulated, t-test and ratio criteria were selected as the final differential proteins. To better understand the WGA- and STL-binding glycoproteins in MCN and SCN, we applied the identified targets of interest to search tool for both UniProt database at https://www.uniprot.org and the Retrieval of Interacting proteins (STRING) 11.0 at http://string-db.org with a focus on the prediction with protein-protein interactions.

3. Result

3.1. Patient Characteristics

About 211 PCNs were evaluated by EUS-FNA from April 2015 to May 2018. Among them, 88 were men and 123 were women. The average age was 57.36 and the standard deviation of the ages was 14.82. 120 patients (56.9%) were pathologically diagnosed through surgical resection, including 50 patients (42.5%) diagnosed with MCNs, 2 patients with MCN-low grade dysplasia, 2 patients with MCN-intermediate grade dysplasia, and 3 patients with MCN-associated invasive adenocarcinoma. 13 cases (10.8%) were diagnosed with IPMN. 46 patients (38.3%) were diagnosed with SCNs, with 9 patients (7.5%) diagnosed with SPN and 2 patients (1.7%) diagnosed with pancreatic neuroendocrine neoplasm (PNEN). A total of 96 cases were diagnosed as MCN or SCN. Sufficient amount of cystic fluid was harvested from 75 patients for glycopattern analysis. The flow chart of patient selection is illustrated in Figure 2.

Figure 2: The flow chart of patients’ selection.

3.2. Glycopattern Analysis of MCN and SCN

The lectin microarray analysis, which included 37 lectins, 2 negative controls (BSA), and 1 positive control (Cy3-BSA), was performed to identify the glycopatterns of the two types of PCNs. The results revealed significant differences in the glycopatterns, as indicated by the white boxes (Figure 3(a)). Normalized relative signal intensities greater than 1.45 or lower than 0.69 were considered valid intensities. Based on intensity ratio of MCN to SCN, differential expressed lectins between two groups of patients were identified (Table 1). WFA, DBA, STL, and WGA were upregulated in MSN patients (Figure 3(b), MCN/SCN ratio > 1.45; p < 0.05), while PTL-1, Con A, and ACA were downregulated in MSN patients (Figure 3(c), MCN/SCN ratio < 0.69; p < 0.05).

Table 1: Differential Glycopattern Analysis of Glycoproteins in Cystic Fluid from MCN versus SCN Using Lectin Microarrays.

No.
Lectin
Binding structure
Normalized relative intensity
MCN/SCN 
 Ratio

MCN
SCN

1
WFA
GalNAcα/β1−3/6Gal
0.026 ± 0.021
0.012 ± 0.009
2.13

2
DBA
GalNAcα-Ser/Thr(Tn) GalNAcαl-3Gal
0.028 ± 0.023
0.016 ± 0.011
1.76

3
STL
GlcNAc oligomer
0.030 ± 0.020
0.021 ± 0.011
1.47

4
WGA
Multivalent Sia and (GlcNAc)n
0.050 ± 0.029
0.031 ± 0.014
1.6

5
BPL
Galβ1-3GalNAc
0.021 ± 0.009
0.013 ± 0.006
1.56

6
PTL-I
αGalNAc and Gal
0.008 ± 0.007
0.013 ± 0.006
0.59

7
ConA
Branched and Terminal Man terminal GlcNAc
0.031 ± 0.029
0.066 ± 0.058
0.47

8
ACA
Galβ1-3GalNAα-Ser/Thr
0.012 ± 0.011
0.021 ± 0.013
0.57

9
MAL-I
Galβ-1,4GlcNAc
0.012 ± 0.007
0.021 ± 0.007
0.57

(a)

(b)

(c)

Figure 3: Glycan-binding protein targeting SCN and MCN by lectin microarray. (a) Lectin microarrays: representative MCN and SCN lectin microarrays: words in red indicate higher signal lectins in MCN compared with SCN, while words in white indicate lower signal lectins in MCN compared with SCN. (b) The recognition glycopattern of lectins was significantly higher in MCN samples than in SCN samples; box plots of lectins (R value > 1.45): WFA, DBA, STL, WGA, and BPL, p < 0.05, p<0.01. (c) The recognition glycopattern of lectins was significantly lower in the MCN samples than in the SCN samples; box plots of lectins (R value < 0.69): PTL-I, ConA, ACA, and MAL-I, p < 0.05, p<0.01.

3.3. Construction of Diagnostic Models Based on Glycopattern Abundances

To test whether these lectins can be used to predict the types of patients, we built a logistic regression model using all lectins in the lectin microarray. The diagnostic accuracy of Model referred to four lectins: Erythrina Cristagalli Letin (ECA), Hippeastrum Hybrid Lectin (HHL), Triticum vulgaris Agglutinin (WGA), and Bauhinia Purpurea Lectin (BPL) were fed into logistic regression. The receiver operating characteristic (ROC) curve indicated that Model could distinguish MCN from SCN with high sensitivity and specificity (AUC: 0.907, sensitivity 0.971, and specificity 0.806) (Figure 4).

Figure 4: The diagnosis accuracy of the diagnostic models and selected lectins analysed by ROC analysis.

3.4. LC-MS/MS Analysis of MCN/SCN Glycoproteins by Differential Lectin Binding

For the lectin affinity enrichment method, this study used the lectins STL and WGA to enrich glycosylated proteins and conducted LC-MS/MS analysis. In the three replications of each sample, the proteins identified in at least one of the replications were included in downstream analysis. 175 glycoproteins were enriched by WGA-coupled magnetic particles, of which 106 proteins and 150 proteins were identified in MSN group and SCN group, respectively. 167 glycoproteins were enriched by STL-coupled magnetic particles, of which 120 proteins and 134 proteins were in MSN group and SCN group, respectively (Table 2). According to the enriched lectins STL and WGA using peptide-spectrum matches (PSM) as the quantitative basis, a four-way WGA-MCN, WGA-SCN, STL-MCN, and STL-SCN Venn diagram shows the distribution of unique and shared proteins (Figure 5(a)).

Table 2: LC-MS/MS analysis of Glycoproteins in MCN versus SCN.

WGA-SCN
WGA-MCN
WGA
STL-SCN
STL-MCN
STL

150
106
175
134
120
167

(a)

(b)

(c)

(d)

(e)

Figure 5: Glycoproteins by differential lectin binding in MCN and SCN. (a) Four-way Venn diagram: WGA- MCN, WGA-SCN, STL-MCN, and STL-SCN in a cross-statistic chart. (b) The Heat map of WGA-binding glycoproteins in MCN and SCN. The samples were listed in columns, and the proteins were listed in rows. The colour and intensity of each square indicated expression levels relative to the other data in the row. Red, high; green, low; black, medium. (c) The Heat map of STL-binding glycoproteins in MCN and SCN. The colour and intensity of each square indicated expression levels relative to the other data in the row. Red, high; green, low; black, medium. (d) String network interaction of APOA1, VTN, C5, and S100A9 genes significantly enriched in WGA-binding glycoproteomics. (e) String network interaction of APOE, ECM1, MUC5B, CP, KLKB1, and PLG genes significantly enriched in STL-binding glycoproteomics.

The heat map in Figures 5(b) and 5(c) shows the expression of WGA- and STL-binding glycoproteins in MCN and SCN accordingly. In the three replications of each sample, these glycoproteins were based on the quantitative LC-MS/MS results, because the glycoproteins are mostly low-abundance proteins [25]. In this study, the label free quantification (LFQ) intensity was found generally in the order of . Sixteen differentially expressed proteins were detected in WGA-binding glycoproteins (p < 0.05), with eight upregulated proteins (MCN/SCN ratio > 1.21) and seven downregulated proteins (MCN/SCN ratio < 0.83) (Table 3). Sixteen differentially expressed proteins were identified in STL-binding glycoproteins (p < 0.05), and ten expressed upregulated glycoproteins and six downregulated glycoproteins (Table 4). Among these thirty-two different glycoproteins, the expressions of four glycoproteins (Apolipoprotein A-I, Actin gamma 1, Haptoglobin, and Clusterin) were increased in both lectin enrichment groups, and the expression of four glycoproteins (Immunoglobulin heavy constant alpha 1, Hemoglobin subunit beta, Apolipoprotein A-II, Alpha-2-macroglobulin) were decreased in both lectin enrichment groups.

Table 3: The protein enrichment by WGA in MCN versus SCN.

Protein name
Gene name
Majority protein IDs
Normalized relative intensity
MCN/SCN 
 Ratio

MCN
SCN

(A) Increased expression

Ceruloplasmin
CP
P00450
3.58 ± 0.45
2.38 ± 0.15
1.50

Haptoglobin
HP
P00738-2
21.43 ± 1.85
11.86 ± 1.31
1.81

Apolipoprotein A-I
APOA1
P02647
203.23 ± 0.80
163.91 ± 9.36
1.24

Fibrinogen beta chain
FGB
P02675
55.94 ± 5.12
32.06 ± 1.62
1.74

Serotransferrin
TF
P02787
47.99 ± 0.83
28.92 ± 4.17
1.66

Protein S100-A9
S100A9
P06702
4.82 ± 0.38
2.08 ± 0.28
2.32

Clusterin
CLU
P10909-4
9.19 ± 0.84
6.99 ± 0.42
1.32

Actin gamma 1
ACTG1
P63261
3.69 ± 0.30
2.59 ± 0.27
1.43

(B) Decreased expression

Serpin peptidase inhibitor clade A member 1
SERPINA1
P01009
41.05 ± 2.37
53.74 ± 2.02
0.76

Alpha-1-antichymotrypsin
AACT
P01011
3.88 ± 0.22
7.91 ± 0.52
0.49

Alpha-2-macroglobulin
A2M
P01023
15.58 ± 2.81
33.14 ± 3.32
0.47

Immunoglobulin heavy constant alpha 1
IGHA1
P01876
27.16 ± 3.34
43.11 ± 7.02
0.63

Apolipoprotein A-II
APOA2
P02652
44.02 ± 1.52
49.89 ± 2.11
0.88

Complement component 9
C9
P02748
3.21 ± 0.14
3.89 ± 0.24
0.83

Plasma protease C1 inhibitor
SERPING1
P05155
9.05 ± 3.04
26.13 ± 5.88
0.35

Hemoglobin subunit beta
HBB
P68871
86.97 ± 13.09
320.97 ± 14.96
0.27

Table 4: The protein enrichment by STL in MCN versus SCN.

Protein name
Gene name
Majority protein IDs
Normalized relative intensity
MCN/SCN 
 Ratio

MCN
SCN

(A) Increased expression

    Ig lambda-7 chain C region
IGLL1
A0M8Q6
12.35 ± 10.23
1.38 ± 0.08
8.95

    Haptoglobin
HP
P00738
15.42 ± 1.67
10.13 ± 0.88
1.52

    Plasminogen
PLG
P00747
14.19 ± 1.29
10.48 ± 0.50
1.35

    Apolipoprotein A-I
APOA1
P02647
195.25 ± 8.29
160.44 ± 14.22
1.22

    Apolipoprotein E
APOE
P02649
48.10 ± 2.26
34.68 ± 0.53
1.39

    Fibrinogen Alpha A
FGA
P02671
44.56 ± 3.31
29.25 ± 1.01
1.52

    Fibrinogen gamma chain
FGG
P02679
35.71 ± 4.92
21.61 ± 1.32
1.65

    Clusterin
CLU
P10909
8.70 ± 1.07
5.36 ± 0.50
1.62

    Inter-alpha-trypsin inhibitor heavy chain H2
ITIH2
P19823
10.27 ± 0.70
6.62 ± 0.96
1.55

    Actin, gamma 1
ACTG1
P63261
4.21 ± 0.32
2.34 ± 0.18
1.8

(B) Decreased expression

    Alpha-2-macroglobulin
A2M
P01023
13.27 ± 0.43
35.82 ± 4.42
0.37

    Ig alpha-1 chain C region
IGHA1
P01876
18.50 ± 1.83
29.41 ± 0.56
0.63

    Apolipoprotein A-II
APOA2
P02652
39.36 ± 7.77
62.11 ± 7.96
0.63

    Apolipoprotein A-IV
APOA4
P06727
2.45 ± 0.09
5.36 ± 1.04
0.46

    Complement component C8 alpha chain
C8A
P07357
2.17 ± 0.12
2.73 ± 0.16
0.8

    Hemoglobin subunit beta
HBB
P68871
97.91 ± 6.10
308.30 ± 10.39
0.32

3.5. STRING Aggregates and KEGG Pathways of DEGs

Using differentially expressed proteins identified, STRING aggregates and KEGG pathway were used to further annotate their functions. As one of the largest databases for protein-protein interaction, STRING aggregates the most available information derived from proteomic context, high throughput experiments, coexpression, and previous knowledge [26]. Using the STRING aggregates http://string-db.org, WGA-binding glycoproteins were mapped into an interaction network, and upregulated glycoproteins were in the red circle (Figure 5(d)), STL-binding proteins were mapped into an interaction network, and 10 increased expression were in the red circle (Figure 5(e)). Certain core genes have been reported to play important roles in pancreatic cystic neoplasms via KEGG pathway analysis. With the KEGG pathway enrichment analysis, the top three enriched pathways were complement and coagulation cascades, coagulation and fibrinolytic activation, and platelet activation (Table 5).

Table 5: KEGG pathway analysis of the DEGs in MCN versus SCN.

Pathway name
Genes number
Genes ID

hsa 04610:Complement and coagulation cascades
10
CLU, A2M, FGA, FGB, FGG, PLG, SERPINA1, SERPING1, C8A, C9

hsa 04611:Platelet activation
4
FGA, FGB, FGG, ACTG1

hsa 04979:Cholesterol metabolism
3
APOA1, APOE, APOA2, APOA4

4. Discussion

Mammalian plasma membranes are known to contain 2% to 10% saccharides, most of which are oligosaccharides and glycoproteins. Linear or branched side chains of glycoproteins may contain two or more monosaccharide residues. At the end of the monosaccharide unit, there is often a negatively charged residue of N-acetylneuraminic acid and a siallic acid [27, 28]. Lectin is currently recognized as a monosaccharide or polysaccharide protein animals and plants and can be used to identify specific saccharide structure [29–31]. Different lectins have different binding affinity to glycan epitopes, which can be divided into mannose type, fucose type, sialic acid type/N-acetylglucosamine type, and Galactose type/N-acetylgalactose type [32, 33]. In order to improve the diagnostic accuracy of lectin microarray for MCN, we chose four lectins—Hippeastrum Hybrid Lectin (HHL) as the representative of mannose type, Triticum vulgaris Agglutinin (WGA) as the representative of sialic acid type/N-acetylglucosamine type, and Erythrina Cristagalli Letin (ECA) and Bauhinia Purpurea Lectin (BPL) as the representative of Galactose type/N-acetylgalactose type. The mathematical model was constructed to distinguish MCN from SCN with highest sensitivity and specificity (Figure 4).

In some researches, Shan Li has reported that oligosaccharide containing GalNAc and GlcNAc structural changes in Huh7 cells surface was associated with Epithelial Mesenchymal Transition (EMT) [34]. Gpc-1, an exosome-associated polysaccharide protein identified in the blood, has been shown to be significantly different among healthy individuals, pancreatic cancer patients, patients with early pancreatic cancer, and those with advanced pancreatic cancer [35]. Moreover, Gbormittah analysed pancreatic cystic tumours and found an increase in lectin WGA, which specifically recognized GlcNAc polymers and polyvalent sialic acid structure, accompanied by increasing level of MUC5AC and CA19-9 [36]. Due to the defect of lectin that one lectin can only enrich one class of glycoproteins, the method of multilectin binding is usually adopted to enrich glycoproteins. In this study, STL/WGA was selected to enrich glycoproteins to compare the differential lectin in MCN/SCN. The significance of the identified proteins was further shown by pathway analysis, protein-protein interactions, and chromosomal location investigations.

Seventy-five cyst fluids were collected by fine needle aspiration directly from the pancreatic cysts. In the cystic fluid of PCLs, Fibrinogen Alpha chain (FGA), fibrinogen beta chain (FGB), and fibrinogen gamma chain (FGG) were markedly elevated in MCN. They have been reported to elevate in prostate cancer, lung cancer, hepatocellular carcinoma, pancreatic cancer, infectious diseases, and myocardial infarction [37, 38]. Ghezzi reported that plasma fibrinogen level might be a potential index to predict prognosis, improving the early diagnosis of endometrial cancer and optimize the treatment schedule [39]. RT-PCR analysis showed that the expression level of FGG mRNA was significantly increased in adjacent tissues of liver cancer tissues, especially in liver cancer metastases [40].

Gene APOA1 is located on a well-known gene cluster at chromosome 11 and the accumulating evidence has shown that single nucleotide polymorphisms (SNPs) in the APOA1/A5/A4/C3 cluster are associated with high triglyceride (TG) levels [41]. Shavva found that the expression of endogenous APOA1 in human monocytes and macrophages was moderate and showed that proinflammatory cytokine tumour necrosis factor α (TNFα) increased APOA1 mRNA and stimulated APOA1 protein secretion by human monocytes and macrophages [42]. APOE was identified to be essential for TG-modulating, and three genes APOA1, APOA2, and APOE were reported to be associated with hypertriglyceridemia, which may increase the risk of acute pancreatitis [43].

S100 calcium binding protein A9 (S100A9) is a calcium- and zinc-binding protein which plays a prominent role in the regulation of inflammatory processes and immune response. It was reported to be dysregulated in many types of cancer, and it induced the overexpression of PD-1/PD-L1, which contributes to ineffective haematopoiesis in myelodysplastic syndromes [44]. Elevated serum S100A9 indicated poor prognosis in hepatocellular carcinoma after curative resection [45].

5. Conclusions

In summary, the study of alterations in glycoproteins associated with pancreatic cystic neoplasms can help in the search for tumour markers to improve cancer diagnosis. It may provide insights into the glycobiology mechanism underlying the progression of pancreatic cystic neoplasms and provide valuable information for the development of new therapies.

Data Availability

The data used to support the findings of this study are included within the article.

Disclosure

The funding sources had no role in writing, data collection, analysis, or interpretation, or any aspect pertinent to this study. The decision to submit this manuscript has been made by the corresponding authors.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Ying Wang, Yufa Sun, and Hanjie Yu conceived and designed the experiments; Jia Feng, Zheng Li, Fuquan Yang, and Xiang Ding performed the experiments; Ying Wang, Yufa Sun, and Enqiang Linghu analysed the data; Ying Wang wrote the manuscript. All authors read, reviewed, and approved the final manuscript. Ying Wang and Yufa Sun have contributed equally to this work.

Acknowledgments

The authors would like to thank all patients involved in this study and thank AJE (American Journal Experts) for language editing. This work was supported by the Beijing Natural Science Foundation (no. 7182066).

References

J. Broughton, J. Lipschitz, M. Cantor, D. Moffatt, A. Abdoh, and A. McKay, “Determining the natural history of pancreatic cystic neoplasms: a manitoban cohort study,” HPB, vol. 18, no. 4, pp. 383–388, 2016.

W. R. Brugge, “Pancreatic cyst surveillance: threat or opportunity?” Gastrointestinal Endoscopy, vol. 83, no. 6, pp. 1118–1120, 2016.

R. Talar-Wojnarowska, M. Pazurek, L. Durko et al., “Pancreatic cyst fluid analysis for differential diagnosis between benign and malignant lesions,” Oncology Letters, vol. 5, no. 2, pp. 613–616, 2013.

J. M. Scheiman, J. H. Hwang, and P. Moayyedi, “American gastroenterological association technical review on the diagnosis and management of asymptomatic neoplastic pancreatic cysts,” Gastroenterology, vol. 148, no. 4, pp. 824–848.e22, 2015.

B. B. Haab, A. Porter, T. Yue et al., “Glycosylation variants of mucins and CEACAMs as candidate biomarkers for the diagnosis of pancreatic cystic neoplasms,” Annals of Surgery, vol. 251, no. 5, pp. 937–945, 2010.

T. Ingkakul, Y. Sadakari, J. Ienaga, N. Satoh, S. Takahata, and M. Tanaka, “Predictors of the presence of concomitant invasive ductal carcinoma in intraductal papillary mucinous neoplasm of the pancreas,” Annals of Surgery, vol. 251, no. 1, pp. 70–75, 2010.

Y. Wang, N. Chai, J. Feng, and E. Linghu, “A prospective study of endoscopic ultrasonography features, cyst fluid carcinoembryonic antigen, and fluid cytology for the differentiation of small pancreatic cystic neoplasms,” Endoscopic Ultrasound, vol. 7, no. 5, pp. 335–342, 2018.

M. Skaro, N. Nanda, C. Gauthier et al., “Prevalence of germline mutations associated with cancer risk in patients with intraductal papillary mucinous neoplasms,” Gastroenterology, vol. S0016-5085, no. 19, p. 30353, 2019.

M. Jones, Z. Zheng, J. Wang et al., “Impact of next-generation sequencing on the clinical diagnosis of pancreatic cysts,” Gastrointestinal Endoscopy, vol. 83, no. 1, pp. 140–148, 2016.

G. K. Ma, D. S. Goldberg, N. Thiruvengadam et al., “Comparing american gastroenterological association pancreatic cyst management guidelines with fukuoka consensus guidelines as predictors of advanced neoplasia in patients with suspected pancreatic cystic neoplasms,” Journal of the American College of Surgeons, vol. 223, no. 5, pp. 729–737, 2016.

M. Kearns and N. A. Ahmad, “Diagnosis and management of pancreatic cystic neoplasms,” Current Treatment Options in Gastroenterology, vol. 15, no. 4, pp. 587–602, 2017.

W. Zhang, E. Linghu, N. Chai, and H. Li, “New criteria to differentiate between mucinous cystic neoplasm and serous cystic neoplasm in pancreas by endoscopic ultrasound: a preliminarily confirmed outcome of 41 patients,” Endoscopic Ultrasound, vol. 6, no. 2, pp. 116–122, 2017.

A. E. Frampton, J. Stebbing, T. M. Gall, B. Silver, L. R. Jiao, and J. Krell, “Activating mutations of GNAS and KRAS in cystic fluid can help detect intraductal papillary mucinous neoplasms of the pancreas,” Expert Review of Molecular Diagnostics, vol. 15, no. 3, pp. 325–328, 2015.

H. Matthaei, J. Wu, M. D. Molin et al., “GNAS sequencing identifies IPMN-specific mutations in a subgroup of diminutive pancreatic cysts referred to as "Incipient IPMNs",” The American Journal of Surgical Pathology, vol. 38, no. 3, pp. 360–363, 2014.

P. Hägglund, J. Bunkenborg, F. Elortza, O. N. Jensen, and P. Roepstorff, “A new strategy for identification of N-glycosylated proteins and unambiguous assignment of their glycosylation sites using HILIC enrichment and partial deglycosylation,” Journal of Proteome Research, vol. 3, no. 3, pp. 556–566, 2004.

K. Alagesan, S. K. Khilji, and D. Kolarich, “It is all about the solvent: on the importance of the mobile phase for ZIC-HILIC glycopeptide enrichment,” Analytical and Bioanalytical Chemistry, vol. 409, no. 2, pp. 529–538, 2017.

M. Yuan, D. M. Kremer, H. Huang et al., “Ex vivo and in vivo stable isotope labelling of central carbon metabolism and related pathways with analysis by LC–MS/MS,” Nature Protocols, vol. 14, no. 2, pp. 313–330, 2019.

T. Yue, I. J. Goldstein, M. A. Hollingsworth, K. Kaul, R. E. Brand, and B. B. Haab, “The prevalence and nature of glycan alterations on specific proteins in pancreatic cancer patients revealed using antibody-lectin sandwich arrays,” Molecular & Cellular Proteomics, vol. 8, no. 7, pp. 1697–1707, 2009.

S. Pan, T. A. Brentnall, and R. Chen, “Proteomics analysis of bodily fluids in pancreatic cancer,” Proteomics, vol. 15, no. 15, pp. 2705–2715, 2015.

M. Lipinski, M. Degowska, and G. Rydzewska, “Cystic fluid neutrophil gelatinase-associated lipocalin (NGAL) concentration in differential diagnosis of pancreatic cystic lesions: a new factor enters the scene?” Przegląd Gastroenterologiczny, vol. 13, no. 2, pp. 132–136, 2018.

Y. Nakai, T. Iwashita, S. Shinoura et al., “Role of serial EUS-guided FNA on pancreatic cystic neoplasms: a retrospective analysis of repeat carcinoembryonic antigen measurements,” Gastrointestinal Endoscopy, vol. 84, no. 5, pp. 780–784, 2016.

Y. Qin, Y. Zhong, G. Yang et al., “Profiling of concanavalin a-binding glycoproteins in human hepatic stellate cells activated with transforming growth factor-β1,” Molecules, vol. 19, no. 12, pp. 19845–19867, 2014.

X. Guo, X. Zhan, and Z. Li, “Molecular analyses of aspirated cystic fluid for the differential diagnosis of cystic lesions of the pancreas: a systematic review and meta-analysis,” Gastroenterology Research and Practice, vol. 2016, Article ID 3546085, 8 pages, 2016.

Y. Qin, Y. Zhong, L. Dang et al., “Alteration of protein glycosylation in human hepatic stellate cells activated with transforming growth factor-β1,” Journal of Proteomics, vol. 75, no. 13, pp. 4114–4123, 2012.

C. Schmelter, S. Funke, J. Treml et al., “Comparison of two solid-phase extraction (SPE) methods for the identification and quantification of porcine retinal protein markers by LC-MS/MS,” International Journal of Molecular Sciences, vol. 19, no. 12, Article ID 3847, 2018.

Y. Zhong, J. Zhang, H. Yu et al., “Characterization and sub-cellular localization of GalNAc-binding proteins isolated from human hepatic stellate cells,” Biochemical and Biophysical Research Communications, vol. 468, no. 4, pp. 906–912, 2015.

Y. Zhong, X. X. Sun, P. Zhang et al., “Identification and localization of xylose-binding proteins as potential biomarkers for liver fibrosis/cirrhosis,” Molecular BioSystems, vol. 12, no. 2, pp. 598–605, 2016.

Y. Qin, Y. Zhong, M. J. Zhu et al., “Age- and sex-associated differences in the glycopatterns of human salivary glycoproteins and their roles against influenza A virus,” Journal of Proteome Research, vol. 12, no. 6, pp. 2742–2754, 2013.

K. W. Oppong, M. F. Dawwas, R. M. Charnley et al., “EUS and EUS–FNA diagnosis of suspected pancreatic cystic neoplasms: Is the sum of the parts greater than the CEA?” Pancreatology, vol. 15, no. 5, pp. 531–537, 2015.

K. Nagata, M. Horinouchi, M. Saitou et al., “Mucin expression profile in pancreatic cancer and the precursor lesions,” Journal of Hepato-Biliary-Pancreatic Surgery, vol. 14, no. 3, pp. 243–254, 2007.

A. Sultana, R. Jackson, G. Tim et al., “What is the best way to identify malignant transformation within pancreatic IPMN: a systematic review and meta-analyses,” Clinical and Translational Gastroenterology, vol. 6, p. e130, 2015.

Y. Ji, S. Wei, J. Hou et al., “Integrated proteomic and N-glycoproteomic analyses of doxorubicin sensitive and resistant ovarian cancer cells reveal glycoprotein alteration in protein abundance and glycosylation,” Oncotarget, vol. 8, no. 8, pp. 13413–13427, 2017.

L. Zou and J. Qian, “Decline of serum CA724 as a probable predictive factor for tumor response during chemotherapy of advanced gastric carcinoma,” Chinese Journal of Cancer Research, vol. 26, no. 4, pp. 404–409, 2014.

S. Li, C. Mo, Q. Peng et al., “Cell surface glycan alterations in epithelial mesenchymal transition process of Huh7 hepatocellular carcinoma cell,” Plos One, vol. 8, no. 8, Article ID e71273, 2013.

S. A. Melo, L. B. Luecke, C. Kahlert et al., “Glypican-1 identifies cancer exosomes and detects early pancreatic cancer,” Nature, vol. 523, no. 7559, pp. 177–182, 2015.

F. O. Gbormittah, B. B. Haab, K. Partyka, C. Garcia-Ott, M. Hancapie, and W. S. Hancock, “Characterization of glycoproteins in pancreatic cyst fluid using a high-performance multiple lectin affinity chromatography platform,” Journal of Proteome Research, vol. 13, no. 1, pp. 289–299, 2014.

C. Chen, L. G. Zhang, J. Liu et al., “Bioinformatics analysis of differentially expressed proteins in prostate cancer based on proteomics data,” OncoTargets and Therapy, vol. 9, Article ID 15451557, 2016.

S. Duan, B. Gong, P. Wang, H. Huang, L. Luo, and F. Liu, “Novel prognostic biomarkers of gastric cancer based on gene expression microarray: COL12A1, GSTA3, FGA and FGG,” Molecular Medicine Reports, vol. 18, no. 4, pp. 3727–3736, 2018.

F. Ghezzi, A. Cromi, G. Siesto et al., “Prognostic significance of preoperative plasma fibrinogen in endometrial cancer,” Gynecologic Oncology, vol. 119, no. 2, pp. 309–313, 2010.

B. L. Fan, W. L. Zhu, G. L. Zou, G. S. Luo, C. L. Xu, and W. X. Zhao, “Cloning and identification of fibrinogen gamma polypeptide (FGG) gene differentially expressed in human hepatocellular carcinoma,” Chinese Journal of Cancer, vol. 23, no. 3, Article ID 249253, 2004.

Q. Lin, Y. Cao, and J. Gao, “Decreased expression of the APOA1-APOC3-APOA4 gene cluster is associated with risk of alzheimer’s disease,” Drug Design, Development and Therapy, vol. 9, pp. 5421–5431, 2015.

V. S. Shavva, D. A. Mogilenko, E. V. Nekrasova et al., “Tumor necrosis factor α stimulates endogenous apolipoprotein A-I expression and secretion by human monocytes and macrophages: role of MAP-kinases, NF-κB, and nuclear receptors PPARα and LXRs,” Molecular and Cellular Biochemistry, vol. 448, no. 1-2, pp. 211–223, 2018.

J. L. Jin, D. Sun, Y. X. Cao et al., “Intensive genetic analysis for chinese patients with very high triglyceride levels: relations of mutations to triglyceride levels and acute pancreatitis,” EBioMedicine, vol. 38, pp. 171–177, 2018.

P. Cheng, E. A. Eksioglu, X. Chen et al., “S100A9-induced overexpression of PD-1/PD-L1 contributes to ineffective hematopoiesis in myelodysplastic syndromes,” Leukemia, 2019.

J. Meng, F. Gu, H. Fang, and B. Qu, “Elevated serum S100A9 indicated poor prognosis in hepatocellular carcinoma after curative resection,” Journal of Cancer, vol. 10, no. 2, pp. 408–415, 2019.