Skip to main content

Tissue proteome analysis for profiling proteins associated with lymph node metastasis in gallbladder cancer

Abstract

Lymph node (LN) metastasis is the earliest sign of metastatic spread and an established predictor of poor outcome in gallbladder cancer (GBC). Patients with LN positive GBC have a significantly worse survival (median survival- 7 months) than patients with LN negative disease (median survival- ~ 23 months) in spite of standard treatment which includes extended surgery followed by chemotherapy, radiotherapy and targeted therapy. This study aims at understanding the underlying molecular processes associated with LN metastasis in GBC. Here, we used iTRAQ-based quantitative proteomic analysis using tissue cohort comprising of primary tumor of LN negative GBC (n = 3), LN positive GBC (n = 4) and non-tumor controls (Gallstone disease, n = 4), to identify proteins associated with LN metastasis. A total of 58 differentially expressed proteins (DEPs) were found to be specifically associated with LN positive GBC based on the criteria of p value ≤ 0.05, fold change ≥ 2 and unique peptides ≥ 2. These include the cytoskeleton and associated proteins such as keratin, type II cytoskeletal 7 (KRT7), keratin type I cytoskeletal 19 (KRT19), vimentin (VIM), sorcin (SRI) and nuclear proteins such as nucleophosmin Isoform 1 (NPM1), heterogeneous nuclear ribonucleoproteins A2/B1 isoform X1 (HNRNPA2B1). Some of them are reported to be involved in promoting cell invasion and metastasis. Bioinformatic analysis of the deregulated proteins in LN positive GBC using STRING database identified ‘neutrophil degranulation’ and ‘HIF1 activation’ to be among the top deregulated pathways. Western blot and IHC analysis showed a significant overexpression of KRT7 and SRI in LN positive GBC in comparison to LN negative GBC. KRT7, SRI and other proteins may be further explored for their diagnostics and therapeutic applications in LN positive GBC.

Peer Review reports

Introduction

Lymph node (LN) metastasis is the initial sign of metastasis and is one of the established predictors of prognosis in gallbladder cancer (GBC) patients. As per AJCC, 8th Edition [1], GBC stage III is catergorized into Stage IIIA (LN negative) and IIIB (LN positive) based on LN status. In LN negative GBC, tumor is spread to liver and/or one nearby extrahepatic organ but not to nearby lymph nodes, while in LN positive GBC, tumor cells disseminates to nearby lymph nodes. After the spread of tumor cells to lymph nodes, the survival of the patients is significantly worse [2]. Five-year survival rates of GBC cases with node negative pN0 was reported to be ~ 80%, while for node positive pN1 and pN2 was ~ 57% and ~ 23% respectively [3, 4]. Negi et al. reported that lymph node ratio (LNR- ratio of metastatic and total number of examined nodes) independently may predict survival after resection [5]. Tumor recurrence is common in LN positive GBC inspite of standard treatment as there is no effective targeted therapy available for these patients. Molecular insights into the mechanisms associated with survival and growth of tumor cells in LNs may be useful for improved therapy to combat LN metastases.

Recent findings provide understanding of the mechanisms causing growth of tumor cells in LNs. These tumor cells as well as the molecules secreted by them decrease anti-tumor immunity and promote tumor growth via engagement of stromal, lymphoid and myeloid cells present in the organ-restricted tumor as well as in the lymphatic system. Vascular endothelial growth factors (VEGFs) e.g. VEGF-C, VEGF-D and several chemokines such as CCL21, CXCL12 and CCL1 facilitate LN metastasis in several cancers [6]. The levels of tissue and serum VEGF-C and D are increased in GBC patients and promote LN metastasis [7-9]. Yao et al.detected the expression of nuclear chemokine receptor type 4 (CXCR4) in the nucleus and/or cytoplasm of GB cancer cells and is reported to be associated with LN metastases [10]. CCR7 is found to be overexpressed in GBC and modulates TNF-α-induced LN metastasis of GBC [11]. Shah et al.evaluated cytokeratin CK7 as marker for detecting micrometastatic disease in lymph nodes (LN with 0.2–2 mm tumor size) using IHC analysis and found low incidence rates (7.5%) of LN micrometastasis in GBC [12].

Application of high throughput proteomic approaches have been successful to understand the complex biological processes and has been applied to understand the molecular processes involved in LN metastasis in various cancers [13-15]. However, there is no high throughput proteomic analysis carried out to study LN metastasis in GBC till date. The present study applied iTRAQ-based quantitative proteomics to identify differentially expressed proteins (DEPs) in primary tumor of LN positive GBC in comparison to LN negative GBC and gallstone disease (GSD) as non-tumor controls. Further, bioinformatic analysis was performed to identify deregulated processes, pathways and networks. The expression of some of the functionally relevant proteins was verified by Western blot analysis.

Methodology

Clinical samples

Adult patients with age ≥ 20 years diagnosed with GBC or GSD cases (non-tumor control) visiting Govind Ballabh Pant Institute of Postgraduate Medical Education and Research (GIPMER), New Delhi, were recruited for the study after approval from the Institutional Human Ethics Committee [Maulana Azad Medical College- Institutional Ethics Committee and ICMR-National Institute of Pathology- Institutional Ethics Committee, New Delhi]. Inclusion criteria includes (1) Patients ≥ 20 year age (2) GBC cases with adenocarcinomas. Exclusion criteria includes (1) Patients < 20 year age (2) GBC cases other than adenocarcinomas (3) presence of biliary tract or any other malignancy other than gallbladder carcinoma (4) Patients who are morbidly ill or having other malignancies (5) Those who have already taken the treatment. Tumor Staging was done on the basis of clinical data of patients, histopathological evaluation and imaging tools, as per AJCC, 8th edition staging system [1]. The tissue sample from GBC cases had ~ 70% tumor cells. Fresh frozen tissue samples from LN negative GBC cases (n = 4), LN positive GBC cases (n = 4) and GSD cases with no dysplasia as ‘non-tumor controls’ (n = 5) were used in this study. Tissue samples were collected immediately after surgical resection from patients with GBC or GSD and stored at -80° C until used for further analysis. Formalin-fixed tissue samples from LN negative GBC cases (n = 15), LN positive GBC cases (n = 15) and GSD cases with no dysplasia as ‘non-tumor controls’ (n = 15) were used for immunohistochemistry (IHC) analysis. Clinico-pathological data of these subjects are detailed in Table 1. Clinical parameters such as TNM, Stage, grade, white cell count, liver enzymes (AST/ ALT/ ALP), bilirubin and co-morbidities (jaundice, diabetes melitus, hypertension, loss of appetite and loss of weight) for the GBC patients and control groups as available (~ 68%) are provided in Supplementary Table S1.

Table 1 Clinico-pathological parameters of the patients used for the study

Protein extraction

Tissue from individual cases [tumor tissue from GBC patients] or controls [GB tissue from GSD cases] was grinded in liquid nitrogen followed by addition of modified RIPA buffer with 2% protease inhibitor cocktail. The tissue homogenate was then sonicated and centrifuged at 13,000 g for 20 min at 4 °C. The supernatant was collected and protein estimation was done using Bradford assay. SDS-PAGE was performed to analyze the protein profile of the tissue lysate from different groups and normalized the protein concentration based on total density [16].

iTRAQ Labeling

For iTRAQ experiments, a pool of GSD tissue lysate (n = 4) was used as control while individual tissue lysate from GBC cases (n = 3 for stage IIIA and n = 4 for stage IIIB) was used for the quantitative proteomic analysis. The experimental design is shown in Supplementary Fig. S1.

In brief, proteins (100 µg) from control (n = 4, pooled sample) and LN Negative GBC (Stage IIIA, n = 3, individual samples) and LN Positive GBC (Stage IIIB, n = 4, individual samples) was reduced, alkylated and digested with trypsin followed by labelling of peptides with 8-plex iTRAQ reagents separately with specific iTRAQ labels (Reagent 113, 114, 115, 116, 117, 118, 119 and 121) as per the manufacturer’s instructions (iTRAQ Reagents Multiplex kit; Applied Biosystems) [16]. The labeled samples were pooled vacuum-dried and subjected to strong cation exchange (SCX) clean up (Cation exchange cartridge, Sciex, US), and desalted using C18 column (Zorbax 300SB-C18, Agilent Technologies, US) as per the manufacturer’s instructions. The samples were then vacuum-dried and used for mass spectrometric analysis (nano-LC MS/MS analysis) [16].

LC–MS/MS analysis

Nanoflow electrospray ionization tandem mass spectrometric analysis was carried out using Orbitrap Fusion (Thermo Scientific, Bremen, Germany) interfaced with Easy-nLC 1000 nanoflow LC system [16]. Peptides from each sample were enriched using a C18 trap column (75 μm × 2 cm) at a flow rate of 3 μl/min and fractionated on an analytical column (75 μm × 50 cm) at a flow rate of 280 nl/min using a linear gradient of 8–60% acetonitrile (ACN) over 46 min. Mass spectrometric analysis was performed in a data dependent manner with a cycle time of 3 s using the Orbitrap mass analyzer at a mass resolution of 120,000 at m/z 200. For each MS cycle, top most intense precursor ions were selected and subjected to MS/MS fragmentation and detected at a mass resolution of 50,000 at m/z 200. The fragmentation was carried out using higher-energy collision dissociation (HCD) mode. Normalized collision energy (CE) of 30% was used to obtain release of reporter ions from all peptides detected in the full scan. The ions selected for fragmentation were excluded for next 30 s. The automatic gain control for full FT MS and FT MS/MS was set to 3e6 ions and 1e5 ions respectively with a maximum time of accumulation of 50 ms for MS and 75 ms for MS/MS. The lock mass with 10 ppm error window option was enabled for accurate mass measurements [16]. The LC–MS/MS analysis was performed in four replicates.

Identification and quantification of proteins

Protein identification, quantification and annotations of DEPs were carried out as described earlier by Priya et al. [17]. The MS/MS data was analyzed using Proteome Discoverer (Thermo Fisher Scientific, version 2.2) with Mascot and Sequest HT search engine nodes using NCBI RefSeq database (release 89). Search parameters included trypsin as the enzyme with 2 missed cleavage allowed; precursor and fragment mass tolerance were set to 10 ppm and 0.1 Da, respectively; Methionine oxidation and deamidation of asparagines and glutamine amino acids was set as a dynamic modification while methylation modification at cysteine and iTRAQ modification at N-terminus of the peptide and lysines were set as static modifications. The peptide and protein information were extracted using high peptide confidence and top one peptide rank filters. The FDR was calculated using percolator node in proteome discoverer 2.2. High confidence peptide identifications were obtained by setting a target FDR threshold of 1% at the peptide level.

The iTRAQ intensity of proteins from each of the four replicates was used for the PCA plot analysis using MetaboAnalyst 5.0 [18] to determine the correlation among the four replicates dataset as well as the correlation of GSD, LN Negative and LN Positive proteome dataset.

Relative quantitation of proteins was carried out based on the intensities of reporter ions released during MS/MS fragmentation of peptides. The proteins identified in all the four replicates were used for the analysis. The average relative intensities of the two reporter ions for each of the unique peptide identifiers for a protein were used to determine relative quantity of a protein and percentage variability. Proteins identified with ≥ 2 unique peptides, with two-fold change or above and FDR adjusted p value < 0.05 were considered significant and used for further analysis [16]. The volcano maps for ‘LN negative GBC vs GSD’ and ‘LN positive GBC vs GSD’ were prepared by using log2 fold change and -log10 (p-value) as the co-ordinates and significant fold change ≥ 2.0 and p-value < 0.05 were considered to screen the proteins.

The data was analyzed for DEPs in individual patient with LN metastasis or without LN metastasis vs GSD (non-tumor control) and represented as Venn diagram. Further, the list of DEPs in LN positive GBC was derived and used for bioinformatics analysis.

Bioinformatic analysis

STRING (http://www.string-db.org) is a database to visualise protein–protein interaction including physical and functional interactions. Mapping of DEPs in LN positive GBC for localization, associated molecular functions, pathways and protein–protein interaction analysis was performed using the STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database by setting up the parameters as Homo Sapiens and combined confidence score greater than 0.4 [19].

Western blot analysis

Proteins were separated by SDS-PAGE and electrotransferred to polyvinyl difluoride (PVDF) membrane [20]. The blots were blocked with 5% skimmed milk powder in TBST [1 × tris buffered saline (10 mM Tris–Cl, pH 7.4 and 30 mM NaCl) with 0.05% Tween 20 and 0.005% Triton-X-100] at RT for 1 h. The blots were then incubated with primary antibodies against Keratin, type II cytoskeletal 7 (KRT7) (1:1000, Cat. No. MA5-11,986, Thermo Scientific, USA) and keratin type I cytoskeletal 19 (KRT19) (1:500, Cat No. MA5-12,663, Thermo Scientific, USA), Sorcin (SRI) (1:1000, Cat. no. PA5-23,143, Thermo Scientific, USA) and Nucleophosmin 1 (NPM1) (1:1000, Cat No. MA5-17,141, Thermo Scientific, USA) at 4 °C overnight. Secondary antibodies against rabbit IgG (1:20,000, Thermo Scientific, USA) was used for SRI while mouse IgG (1:20,000, Thermo Scientific, USA) was used for KRT7, KRT19 and NPM1 respectively at RT for 1 h. The blots were developed using the enhanced chemiluminescent (ECL) Kit (Millipore, USA). The images were acquired using Chemidoc MP imager and immunoblots were analyzed using Image Lab 4.1 software (Bio-Rad, USA). Densitometric analysis of the specific band showing reactivity was performed to get relative expression of target proteins in LN positive (n = 4) vs LN negative (n = 4) and LN positive (n = 4) vs GSD (non-tumor control) (n = 4). The maximum density of GSD cases was used to calculate the fold change in expression in LN positive GBC. Statistical analysis was performed using GraphPad Prism 5 [21]. Differences in expression of target proteins between LN positive vs LN negative, LN positive GBC vs GSD and LN negative GBC vs GSD were tested with unpaired t-test (two-tailed) with confidence intervals of 95%. The p-values less than 0.05 was used to indicate statistical significance.

Immunohistochemistry analysis

IHC was performed on FFPE tissues using individual tissue sections from non-tumor controls, GSD cases (n = 15), LN negative GBC (n = 15) and LN positive GBC cases (n = 15) (Supplementary Table S1) to analyze the expression of KRT7, KRT19 and SRI protein. IHC analysis was performed as described earlier by Akhtar et al. [20]. In brief, after deparaffinization and rehydration of FFPE tissue sections, antigen retrieval was performed by immersing the slide in antigen retrieval buffer (20 mM Tris buffer, pH 9.0) at 90 °C for 20 min. Endogenous peroxidases were blocked with 0.03% hydrogen peroxide, and nonspecific binding was blocked with protein blocking reagent. Sections were then incubated for 1 h at RT with primary antibody against KRT7 (dilution 1:250, catalogue no. 307 M-96, Merck, USA), KRT19 (ready-to-use, catalogue no. PR138, Pathnsitu Biotechnologies, USA) and SRI (dilution 1:3000, catalogue no. PA5-23,143, Thermo Scientific, USA) followed by incubation with PolyExcel PolyHRP for 40 min at RT. Tissue sections were then incubated with Stunn DAB working solution for 5 min at RT (PathnSitu Biotechnologies, USA). Sections were counter stained with Mayer’s hematoxylin, dehydrated and images were taken under the microscope. The distribution of staining and staining intensity across the section was observed under the microscope. The scoring criteria were based on both staining intensity and distribution. For KRT7, KRT19 and SRI, the 2 + or higher intensity, with ≥ 25% distribution was considered as ‘Positive’, while 1 + positivity or < 25% distribution was considered as ‘Negative’. The data was analyzed for nuclear expression only as well as both nuclear and cytoplasmic expression for SRI. IHC data analysis was done by two independent pathologists. The statistical analysis (Chi square test, two-tailed) was performed using GraphPad Prism 5 [21] to study the correlation of KRT7, KRT17 and SRI expression among LN positive vs LN negative GBC, LN positive vs GSD and LN negative vs GSD. The p-value less than 0.05 indicated statistical significance.

Results

In the present study, we performed the differential tissue proteome profiling to identify the proteins and associated molecular processes and pathways in LN positive GBC followed by verification of selected proteins based on their significant fold change in LN metastatic GBC and their functional association with LN metastasis using Western blot and/or IHC analysis. The overall work plan of the study is shown in Fig. 1.

Fig. 1
figure 1

Overall workflow of the study. GSD- Gallstone disease; GBC- Gallbladder carcinoma, LN- Lymph node

Identification of differentially expressed proteins in LN positive GBC

Tissue lysate was prepared using individual samples followed by protein estimation by Bradford assay. SDS-PAGE analysis was performed to ensure equal loading using densitometric analysis (Supplementary Fig. S2). Further, iTRAQ based quantitative proteomic analysis to identify proteins associated with LN metastasis in GBC.

The mass spectrometric analysis led to the identification of a total of 1076 proteins. A comparison of four technical replicates of mass spectrometric runs was performed and the data is presented as Venn diagram (Supplementary Fig. S3). The analysis revealed around 700 proteins in each replicate run and a total of 468 proteins were found to be common among the four replicates. PCA plot analysis for the assessment of replicate data showed a significant correlation among the four replicates for each sample (Fig. 2). We observed the non-tumor control group and LN negative GBC to be more similar to each other, while patients with LN positive GBC clustered as a distinct group. We found 129 DEPs in LN Negative GBC and a total of 132 DEPs in LN Positive GBC with ≥ two-fold change and adjusted p-value ≤ 0.05.

Fig. 2
figure 2

Principal component analysis for the control, LN negative and LN positive cases. We observed a significant correlation among the four technical replicates for each sample and found the majority of the LN negative cases group together. PCA was performed based on the log2-transformed intensity of proteins identified in the four replicates. Green area is depicting the LN negative and blue area indicates LN positive GBC. Red dots indicate controls. The PCA plot is derived using the iTRAQ reporter intensity from the quantitative proteomics data using metaboanalyst

The comparison of LN positive and LN negative group showed 53 proteins specific to LN positive GBC while 79 proteins were common between the two stages (Fig. 3). Among the 79 common proteins, expression of 5 proteins showed opposite correlation. Overall, we found 58 proteins (53 + 5) to be differentially expressed in LN positive GBC and the list is provided in Table 2. This includes functionally relevant proteins associated with LN metastasis in cancer such as KRT7, KRT19, SRI, NPM1, Annexin A2 (ANXA2), Annexin A5 (ANXA5). The detailed protein and peptide list is provided in Supplementary Table S2.

Fig. 3
figure 3

Venn Diagram showing the comparison of DEPs in LN negative and LN positive GBC. A total of 79 proteins are common to both stage LN negative and LN positive GBC, while 50 proteins were specific to LN negative GBC and 53 proteins are specific to LN positive. A total of 5 proteins showed opposite correlation in expression among the proteins common to both the stages. n signifies the number of non redundant proteins in both the stages

Table 2 List of 58 proteins differentially expressed in LN positive GBC. This includes 53 proteins specific to LN positive GBC (A) and 5 proteins showed opposite correlation in expression among the proteins common to both the stages (B). The detailed protein and peptide list is provided in Supplementary Table S1

In order to screen the most relevant proteins, we applied a criteria of significant fold change (≥ two-fold, p value < 0.05) in ≥ 50% LN positive GBC cases and no change in expression level (≤ 1.3 fold change) in LN negative GBC cases and found five proteins (KRT7, KRT17, SRI, NPM1 and HNRNPA2B1) showed overexpression (≥ two-fold, p value < 0.05) in ≥ 50% of the LN positive patients and ≤ 1.3 fold change in LN negative patients.

Bioinformatic analysis

A gene ontology analysis for the localization of 58 DEPs showed that 31% of them belong to cytosol, 25% belong to cytoskeleton, 17% are associated to endoplasmic reticulum, 12% are from extracellular region, 7% are associated to mitochondria and 5% are associated with other localization (Fig. 4A). The top molecular functions include RNA binding, unfolded protein binding, oxidoreductase activity and MHC class II protein complex binding (Fig. 4B, Supplementary Table S3A) and the significantly altered pathways include ‘neutrophil degranulation’ and ‘HIF1 activation’ (Fig. 4C, Supplementary Table S3B). ‘Tumor associated neutrophils’ has been shown to be positively associated with LN metastasis in early gastric cancer [22] ‘neutrophil degranulation’ is reported to be a contributing factor to LN metastasis in breast cancer [23]. Pezutto et al. reported the role of HIF1 activation in cancer progression via regulation of VEGFs and endothelial mesenchymal transition (EMT) related genes [24].

Fig. 4
figure 4

Gene ontology of 58 proteins DEPs in LN positive GBC. (A) Subcellular localisation according to UniProt database. (B) Molecular functions and (C) Reactome pathways according to STRING database

We further analyzed the DEPs within individual patients and the data is represented as Volcano plots in Supplementary Fig. S4. Protein–Protein Interaction Networks using 58 DEPs showed hub molecules including HSP90AA1, HSPA8, HSPA9, MDH2, PKM, HNRNPA1, HNRNPK, NPM1 and ANXA2. Of these, HSPA8, HSPA9, PKM, HNRNPA1, NPM1, ANXA2 are already reported to be associated with LN metastasis in GBC or other cancers (Fig. 5).

Fig. 5
figure 5

Protein–protein interaction analysis of 58 DEPs in LN positive GBC using STRING database. We observed HSP90AA1, HSPA8, HSPA9, MDH2, PKM, HNRNPA1, HNRNPK, NPM1 and ANXA2 as the hub molecules. The proteins HSPA8, HSPA9, PKM, HNRNPA1, NPM1, ANXA2 are already reported to be associated with LN metastasis in GBC or other cancers

Clinical verification using Western blot and IHC analysis

Based on the significant fold change in LN positive GBC in mass spectrometric data and literature survey for their functional relevance and association with LN metastasis in cancer, we selected four proteins namely KRT7, KRT19, SRI and NPM1 for verification using Western blot analysis. The expression level of these proteins as per the quantitative proteomics data is shown in Fig. 6.

Fig. 6
figure 6

Altered levels of functionally relevant proteins in individual GBC patients as observed in quantitative proteomics data. The bar diagram showing the upregulated levels of KRT7, KRT19, SRI, NPM1 GBC LN positive GBC. LN Negative- GBC 1–3, LN Positive- GBC 4–7

We found a significant overexpression of three proteins, KRT7, KRT19 and SRI, with p value 0.03, 0.04 and 0.04 respectively while NPM1 did not show any significant difference in expression between the two groups. Western blot results are shown in Fig. 7 and the full-length blot images are presented in Supplementary Fig. S5.

Fig. 7
figure 7

KRT7, KRT19, SRI and NPM1 protein expression in the LN positive GBC, LN negative GBC and GSD cases. (A) Western blot images showing expression of KRT7, KRT19, SRI and NPM1 in the individual tissue samples from three groups. (B) A significant overexpression (p value < 0.05) of KRT7, KRT19 and SRI was found in LN positive GBC cases in comparison to LN negative GBC cases or GSD cases (non-tumor control). No significant change in expression was observed for NPM1 in LN positive cases

Based on the Western blot results, we performed IHC analysis to study the expression of three proteins, KRT7, KRT19 and SRI, in GSD (non-tumor control), LN negative GBC and LN positive GBC. The expression of KRT7 was found to be ‘positive’ in 100% GSD cases, 33% LN negative GBC and 86% LN positive GBC. KRT19 expression was found to be ‘positive’ in 100% GSD cases, 93% LN negative GBC and 100% LN positive GBC. The expression of SRI was found to be ‘positive’ in 6% GSD cases, 20% LN negative GBC and 60% LN positive GBC. Figure 8A shows the representative IHC images of controls, LN negative GBC and LN positive GBC. The statistical analysis between LN positive vs LN negative GBC showed a significant difference in KRT7 and SRI expression (p value 0.002 and 0.02 respectively) while KRT19 did not show any significant difference in expression in LN positive GBC (Fig. 8B).

Fig. 8
figure 8

IHC analysis to study the expression of KRT7, KRT19 and SRI in the GBC cases and control group. (A) Representative IHC images showing the expression of KRT7, KRT19 and SRI in controls, LN Negative GBC and LN Positive GBC cases. IHC was performed on formalin-fixed paraffin-embedded (FFPE) individual tissue sections of 15 controls (GSD cases with no dysplasia), 15 LN Negative GBC cases and 15 LN Positive GBC cases. The IHC results showed that the expression of KRT7 was found to be ‘positive’ in 100% GSD cases, 33% LN negative GBC and 86% LN positive GBC. The expression of KRT19 was found to be ‘positive’ in 100% GSD cases, 93% LN negative GBC and 100% LN positive GBC. The expression of SRI was found to be ‘positive’ in 6% GSD cases, 20% LN negative GBC and 60% LN positive GBC. (B) The statistical analysis between LN positive vs LN negative GBC showed a significant difference in KRT7 and SRI expression (p value 0.002 and 0.02 respectively) while KRT19 did not show any significant difference in expression. The p values ≤ 0.05, ≤ 0.01 are marked with ‘*’ and ‘**’ respectively. The scale bar is shown as white line. IHC scoring is shown in Methodology section 2.8. GSD- Gallstone disease, LN- Lymph node

Discussion

LN metastasis is a major prognostic factor for patients with GBC as cancer cells from LNs may spread to other metastatic sites by lymphatic or blood vessels. It would be important to carry out an indepth analysis to understand the molecular processes associated with LN metastasis. Here, we have used iTRAQ-based quantitative proteomics approach and analyzed the differential proteome in LN metastatic GBC which led to the identification of 58 DEPs. Of these, 6 proteins are earlier reported to be associated with LN metastasis in GBC (KRT7, KRT19, VIM, PKM, LASP-1, TYMP), 33 proteins are associated with LN metastasis in other cancers (NPM1, SRI, FABP5, PPIA, HSPA9, AOC3, ANXA2, ANXA5, ANXA10, TAGLN2, HNRNPK, CAPG) while 19 proteins (DLST, MDH2, CKMT1A, APOC3, RDH16, VCAN) are novel to LN metastasis in cancer (Supplementary Table S4). Out of 58 proteins, five proteins (KRT7, KRT17, SRI, NPM1 and HNRNPA2B1) showed overexpression (≥ two-fold, p value < 0.05) in ≥ 50% of the LN positive patients and ≤ 1.3 fold change in LN negative patients. Three of these proteins (KRT7, KRT17 and SRI) are cytoskeleton proteins or cytoskeletal associated proteins and two proteins (NPM1 and HNRNPA2B1) are nuclear proteins. We further analyzed the expression of four proteins, KRT7, KRT19, SRI and NPM1 by Western blot and/or analysis.

The various components of ‘cytoskeleton and their associated proteins’ are well integrated in normal cells, however, remodeling of these proteins is well-established in cancer cells. Cytoskeleton proteins not only facilitate the invasion and migration of tumor cells but also regulate their intercellular signaling to facilitate tumor progression [25]. Some of the cytoskeleton proteins such as actin, vimentin, cytokeratins, transgelins are already reported to be involved in LN metastasis in various cancers [26-29]. In the present study, we observed cytoskeleton and associated proteins (KRT7, KRT17, KRT19, KRT6A, SRI) to be significantly overexpressed in LN positive GBC (Fig. 4A).

Cytokeratins (CKs or KRTs), the major group of proteins identified in this study, belongs to the intermediate filament (IF) protein family and epithelial cell markers that are routinely used in cancer diagnostics. It has been reported that cytokeratins release from proliferating or apoptotic cells, distinctly reflecting ongoing cell activity which makes them useful markers for epithelial malignancies [30]. Increased tissue expression levels of KRT7 is reported to be associated with LN metastasis in gallbladder cancer, lung cancer and colorectal cancer patients [28, 31, 32], KRT17 in gastric cancer, esophageal carcinoma and papillary thyroid carcinoma [33-35], KRT6A in lung adenocarcinoma [36], KRT19 in adenocarcinoma of gallbladder [37], thoracic tumors [38]. KRT7 and KRT19 staining of lymph nodes is reported to have a potential in detecting micrometastatic foci in regional lymph nodes of patients with GBC [12, 32] and cervical cancer [39] respectively. KRT19 expression was found to be the independent prognostic factor for hepatocellular carcinoma [40] and pancreatic ductal adenocarcinoma with LN metastasis [10]. Wang et al. reported KRT17 knockdown inhibited the invasion and proliferation of lung cancer cells. Moreover, its overexpression upregulated the activity of β-catenin and expression levels of Wnt target genes, such as cyclin D1, c-Myc, and MMP7 and also promoted EMT by increasing the levels of Vimentin [41]. Liu et al., 2019 explored the potential of KRT17 in proliferation and metastasis in esophageal squamous carcinoma. The volume and weight of KRT17 knockout tumors were smaller than the controls in vivo [35]. Here, we analyzed the expression of KRT7 and KRT19 by Western blot and IHC analysis. We found a significant overexpression of both the proteins in LN positive GBC in comparison to LN negative or GSD cases by Western blot analysis. However, IHC analysis confirmed the overexpression of KRT7 in LN positive GBC in comparison to LN negative GBC while KRT19 did not show any difference between the two groups. The previous studies on CK7 expression have been carried out in lymph node tissue (metastatic vs non-metastatic lymph nodes). Here, we confirm its overexpression in the primary tumor from LN positive GBC cases. All the GSD cases showed high expression of KRT7 and KRT19. The ‘low’ expression of KRT7 and KRT19 in GSD cases observed in proteomics data and Western blot data could be due to the presence of mixed tissue (epithelial mucosa, muscle tissue, connective tissue).

We also found sorcin, a cytoskeleton-associated protein, to be differentially expressed in LN positive GBC in the proteomics dataset. Sorcin is a 22 kDa soluble resistance related calcium-dependent protein reported to be involved in multidrug resistance in cancer [42]. This protein is associated with vimentin, a cytoskeletal protein already reported to be linked to LN metastasis in GBC. Deng et al.reported the association of sorcin with LN metastasis in gastric cancer patients by immunohistochemistry analysis [43]. It is well reported to play a role in EMT via downregulating the levels of epithelial marker ‘E-cadherin’ and upregulating the levels of mesenchymal marker ‘vimentin’ [44]. Sorcin exert its oncogenic effect by regulating key molecules such as VEGFs, MMPs, NF-κB, ERK1/ 2, CTSZ, Akt, STAT3, caspases involved in carcinogenesis and invasion of tumor cells and modulating signaling pathways including ERK, MAPK/ERK, and PI3K/Akt in various cancers. Supporting this fact, Hu et al.reported that sorcin silencing led to acquisition of epithelial-like morphology, attenuation of EMT and suppression of breast cancer metastasis in vivo [33]. Tong et al.reported the regulatory role of sorcin in EMT in colorectal cancer [45]. In the present study, the overexpression of SRI in LN positive GBC was confirmed by Western blot and IHC analysis. Sorcin is reported to be localized in nucleus and cytoplasm in various cancers as per HPA data (https://www.proteinatlas.org/ENSG00000075142-SRI/pathology). We observed both nuclear and cytoplasmic expression of sorcin in GBC. We found a significant difference in nuclear SRI expression in LN positive GBC in comparison to LN negative GBC, however, no significant difference in SRI expression was observed among these two groups when both nuclear and cytoplamic expression were considered. Deng et al.reported its cytoplasmic expression in gastric cancer and was associated with LN metastasis [43], whereas, we find its nuclear expression to be linked to LN metastasis in GBC. The role of nuclear SRI in LN metastasis is not clear and needs to be explored further.

We also observed an overexpression of two nuclear proteins, HNRNPA2B1 and NPM1 in our quantitative proteomics data. HNRNPA2B1 is an RNA binding protein involved in the transcription, splicing processing, transport, stability, telomere maintenance and DNA repair and its overexpression is reported to mediate EMT in different cancers [46, 47]. Higher HNRNPA2B1 gene expression was found to be associated with LN metastasis in esophageal carcinoma [48]. NPM1 is a nucleocytoplasmic shuttling protein having an essential role in cellular processes such as maintenance of genomic stability, DNA repair, centrosome duplication, ribosome biogenesis, cell cycle progression and regulation of activity of tumor suppressor genes p53 and ARF [49]. Correlation of higher NPM1 expression with LN metastasis was also found in patients with oral squamous carcinoma and colon cancer [50, 51]. Liu et al.reported an increased expression of NPM1 and its association with LN metastasis poor survival of colon cancer patients. Small interfering RNA (siRNA) of NPM1 inhibited migration and invasiveness of metastatic colon cancer in the HCT116 cell line [51]. Zhang et al.also reported inhibition of proliferation of HepG2 cells on NSC348884 treatment (a small molecular inhibitor of NPM1) [52]. Qi et al.explored the potential of NSC348884 as an anticancer drug and reported upregulation of p53 and apoptosis in a dose-dependent manner after treatment of several different cancer cell types with the inhibitor. NPM1 inhibitor also has a synergistic effect with chemotherapeutic drugs [53]. The expression of NPM1 was analyzed using Western blot analysis and we observed the overexpression of NPM1 in one of the LN positive GBC in comparison to LN negative GBC (25% positivity), however, based on the statistical analysis there was no significant difference in expression between LN positive GBC vs LN negative GBC.

Over all, the study identified a total of 58 proteins associated with LN metastasis in GBC which includes the cytoskeleton proteins and nuclear proteins. The overexpression of KRT7 and SRI in LN positive GBC in comparison to LN negative GBC was confirmed by Western blot and IHC analysis. The limitations include the low sample size used in the study. Also, the verification was restricted to the primary tumor and was not performed in the metastatic lymph nodes. KRT7 has been previously reported for its potential use in detection of micrometastasis in GBC and our study confirmed high positivity rate in primary tumor (gallbladder) from LN metastatic GBC. Further, the functional role of SRI using knockdown studies may be investigated in GBC cell lines for their potential as therapeutic targets for LN positive GBC cases.

Availability of data and materials

All data generated or analyzed for this study is included in the main article and supplementary information files and is also available with the corresponding author.

References

  1. Giannis D, Cerullo M, Moris D, Shah KN, Herbert G, Zani S, et al. Validation of the 8th Edition American Joint Commission on Cancer (AJCC) Gallbladder Cancer Staging System: Prognostic Discrimination and Identification of Key Predictive Factors. Cancers (Basel). 2021;13(3):547.

    Article  PubMed  Google Scholar 

  2. Behari A, Sikora SS, Wagholikar GD, Kumar A, Saxena R, Kapoor VK. Long term survival after extended resections in patients with gallbladder cancer. J Am Coll Surg. 2003;196(1):82–8.

    Article  PubMed  Google Scholar 

  3. Shirai Y, Sakata J, Wakai T, Ohashi T, Ajioka Y, Hatakeyama K. Assessment of lymph node status in gallbladder cancer: location, number, or ratio of positive nodes. World J SurgOncol. 2012;10:87.

    Google Scholar 

  4. Choi BG, Kim CY, Cho SH, Kim HJ, Koh YS, Kim JC, et al. Impact of lymph node ratio as a valuable prognostic factor in gallbladder carcinoma, focusing on stage IIIB gallbladder carcinoma. J Korean Surg Soc. 2013;84(3):168–77.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Negi SS, Singh A, Chaudhary A. Lymph nodal involvement as prognostic factor in gallbladder cancer: location, count or ratio? J Gastrointest Surg. 2011;15(6):1017–25.

    Article  PubMed  Google Scholar 

  6. Jones D, Pereira ER, Padera TP. Growth and Immune Evasion of Lymph Node Metastasis. Front Oncol. 2018;8:36.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Jiang L, Chen YL, She FF, Tang NH, Li XJ, Wang XX. Expressions of VEGF-C and VEGF-D and their correlation with lymphangiogenesis and angiogenesis in gallbladder carcinoma. ZhonghuaZhong Liu ZaZhi. 2010;32(3):190–5.

    CAS  Google Scholar 

  8. Liu MC, Jiang L, Hong HJ, Meng ZW, Du Q, Zhou LY, et al. Serum vascular endothelial growth factors C and D as forecast tools for patients with gallbladder carcinoma. Tumour Biol. 2015;36(8):6305–12.

    Article  CAS  PubMed  Google Scholar 

  9. Jiang L, Liu M, Cai X, Xie L, She F, Chen Y. Serum vascular endothelial growth factor-C levels predict lymph node metastasis and prognosis of patients with gallbladder cancer. OncolLett. 2018;16(5):6065–70.

    CAS  Google Scholar 

  10. Yao H, Yang Z, Liu Z, Miao X, Yang L, Li D, et al. Glypican-3 and KRT19 are markers associating with metastasis and poor prognosis of pancreatic ductal adenocarcinoma. Cancer Biomark. 2016;17(4):397–404.

    Article  CAS  PubMed  Google Scholar 

  11. Hong H, He C, Zhu S, Zhang Y, Wang X, She F, et al. CCR7 mediates the TNF-α-induced lymphatic metastasis of gallbladder cancer through the “ERK1/2 - AP-1” and “JNK - AP-1” pathways. J Exp Clin Cancer Res. 2016;35:51.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Shah SH, Gupta N, Gupta G, Mehta A, Singh S. Lymph node micrometastasis in gallbladder cancer. Indian J Gastroenterol. 2017;36(4):268–74.

    Article  PubMed  Google Scholar 

  13. Mori K, Toiyama Y, Otake K, Ide S, Imaoka H, Okigami M, et al. Successful identification of a predictive biomarker for lymph node metastasis in colorectal cancer using a proteomic approach. Oncotarget. 2017;8(63):106935–47.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Ma Y, Li YF, Wang T, Pang R, Xue YW, Zhao SP. Identification of proteins associated with lymph node metastasis of gastric cancer. J Cancer Res Clin Oncol. 2014;140(10):1739–49.

    Article  CAS  PubMed  Google Scholar 

  15. Suzuki Y, Takadate T, Mizuma M, Shima H, Suzuki T, Tachibana T, et al. Stromal expression of hemopexin is associated with lymph-node metastasis in pancreatic ductal adenocarcinoma. PLoS ONE. 2020;15(7): e0235904.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Akhtar J, Jain V, Kansal R, Priya R, Sakhuja P, Goyal S, et al. Quantitative tissue proteome profile reveals neutrophil degranulation and remodeling of extracellular matrix proteins in early stage gallbladder cancer. Front Oncol. 2023;12:1046974.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Priya R, Jain V, Akhtar J, Chauhan G, Sakhuja P, Goyal S, et al. Plasma-derived candidate biomarkers for detection of gallbladder carcinoma. Sci Rep. 2021;11(1):23554.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Pang Z, Zhou G, Ewald J, et al. Using MetaboAnalyst 5.0 for LC–HRMS spectra processing, multi-omics integration and covariate adjustment of global metabolomics data. Nat Protoc. 2022;17:1735–61.

    Article  CAS  PubMed  Google Scholar 

  19. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.

    Article  CAS  PubMed  Google Scholar 

  20. Akhtar J, Priya R, Jain V, Sakhuja P, Agarwal AK, Goyal S, et al. Immunoproteomics approach revealed elevated autoantibody levels against ANXA1 in early stage gallbladder carcinoma. BMC Cancer. 2020;20(1):1175.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Unpaired t-test and Chi Square test were done using GraphPad Prism software version 5.0 www.graphpad.com. Accessed February 2022.

  22. Wang Y, Zhai J, Zhang T, Han S, Zhang Y, Yao X, et al. Tumor-Associated Neutrophils Can Predict Lymph Node Metastasis in Early Gastric Cancer. Front Oncol. 2020;10: 570113.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Stevenson J, Barrow-McGee R, Yu L, Paul A, Mansfield D, Owen J, et al. Proteomics of REPLICANT perfusate detects changes in the metastatic lymph node microenvironment. NPJ Breast Cancer. 2021;7(1):24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Pezzuto A, Carico E. Role of HIF-1 in Cancer Progression: Novel Insights. A Review Curr Mol Med. 2018;18(6):343–51.

    Article  CAS  PubMed  Google Scholar 

  25. Fife CM, McCarroll JA, Kavallaris M. Movers and shakers: cell cytoskeleton in cancer metastasis. Br J Pharmacol. 2014;171(24):5507–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Lee DY, Kang Y, Im NR, Kim B, Kwon TK, Jung KY, et al. Actin-Associated Gene Expression is Associated with Early Regional Metastasis of Tongue Cancer. Laryngoscope. 2021;131(4):813–9.

    Article  CAS  PubMed  Google Scholar 

  27. Liu S, Liu L, Ye W, Ye D, Wang T, Guo W, et al. High Vimentin Expression Associated with Lymph Node Metastasis and Predicated a Poor Prognosis in Oral Squamous Cell Carcinoma. Sci Rep. 2016;6:38834.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Luo HT, Liang CX, Luo RC, Gu WG. Identification of relevant prognostic values of cytokeratin 20 and cytokeratin 7 expressions in lung cancer. Biosci Rep. 2017;37(6):BSR20171086.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Dvorakova M, Nenutil R, Bouchal P. Transgelins, cytoskeletal proteins implicated in different aspects of cancer development. Expert Rev Proteomics. 2014;11(2):149–65.

    Article  CAS  PubMed  Google Scholar 

  30. Barak V, Goike H, Panaretakis KW, Einarsson R. Clinical utility of cytokeratins as tumor markers. Clin Biochem. 2004;37(7):529–40.

    Article  CAS  PubMed  Google Scholar 

  31. Fei F, Li C, Cao Y, Liu K, Du J, Gu Y, et al. CK7 expression associates with the location, differentiation, lymph node metastasis, and the Dukes’ stage of primary colorectal cancers. J Cancer. 2019;10(11):2510–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Tajima Y, Tomioka T, Ikematsu Y, Ichinose K, Inoue K, Kanematsu T. Immunohistochemical demonstration of cytokeratin is useful for detecting micrometastatic foci from gallbladder carcinoma in regional lymph nodes. Jpn J Clin Oncol. 1999;29(9):425–8.

    Article  CAS  PubMed  Google Scholar 

  33. Hu Y, Li S, Yang M, Yan C, Fan D, Zhou Y, et al. Sorcin silencing inhibits epithelial-to-mesenchymal transition and suppresses breast cancer metastasis in vivo. Breast Cancer Res Treat. 2014;143(2):287–99.

    Article  CAS  PubMed  Google Scholar 

  34. Kim HS, Lee JJ, Do SI, Kim K, Do IG, Kim DH, et al. Overexpression of cytokeratin 17 is associated with the development of papillary thyroid carcinoma and the presence of lymph node metastasis. Int J Clin Exp Pathol. 2015;8(5):5695–701.

    PubMed  PubMed Central  Google Scholar 

  35. Liu Z, Yu S, Ye S, Shen Z, Gao L, Han Z, et al. Keratin 17 activates AKT signalling and induces epithelial-mesenchymal transition in oesophageal squamous cell carcinoma. J Proteomics. 2020;211: 103557.

    Article  CAS  PubMed  Google Scholar 

  36. Yang B, Zhang W, Zhang M, Wang X, Peng S, Zhang R. KRT6A Promotes EMT and Cancer Stem Cell Transformation in Lung Adenocarcinoma. Technol Cancer Res Treat. 2020;19:1533033820921248.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Li J, Yang ZL, Zou Q, Yuan Y, Li J, Liang L, et al. Squamous cell/adenosquamous carcinomas and adenocarcinomas of the gallbladder: an immunohistochemistry study of prognostic markers. Pathol Oncol Res. 2014;20(2):285–92.

    Article  CAS  PubMed  Google Scholar 

  38. Masai K, Nakagawa K, Yoshida A, Sakurai H, Watanabe SI, et al. Cytokeratin 19 expression in primary thoracic tumors and lymph node metastases. Lung Cancer. 2014;86(3):318–23.

    Article  PubMed  Google Scholar 

  39. Wang HY, Sun JM, Lu HF, Shi DR, Ou ZL, Ren YL, et al. Micrometastases detected by cytokeratin 19 expression in sentinel lymph nodes of patients with early-stage cervical cancer. Int J Gynecol Cancer. 2006;16(2):643–8.

    Article  CAS  PubMed  Google Scholar 

  40. Zhuang PY, Zhang JB, Zhu XD, Zhang W, Wu WZ, Tan YS, et al. Two pathologic types of hepatocellular carcinoma with lymph node metastasis with distinct prognosis on the basis of CK19 expression in tumor. Cancer. 2008;112(12):2740–8.

    Article  PubMed  Google Scholar 

  41. Lee J, Park HL, Jeong CW, Ha J, Jo K, Kim MH, et al. CYFRA 21–1 in Lymph Node Fine Needle Aspiration Washout Improves Diagnostic Accuracy for Metastatic Lymph Nodes of Differentiated Thyroid Cancer. Cancers (Basel). 2019;11(4):487.

    Article  PubMed  Google Scholar 

  42. Colotti G, Poser E, Fiorillo A, Genovese I, Chiarini V, Ilari A. Sorcin, a calcium binding protein involved in the multidrug resistance mechanisms in cancer cells. Molecules. 2014;19(9):13976–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Deng L, Su T, Leng A, Zhang X, Xu M, Yan L, et al. Upregulation of soluble resistance-related calcium-binding protein (sorcin) in gastric cancer. Med Oncol. 2010;27(4):1102–8.

    Article  CAS  PubMed  Google Scholar 

  44. Shabnam B, Padmavathi G, Banik K, Girisa S, Monisha J, Sethi G, et al. Sorcin a Potential Molecular Target for Cancer Therapy. Transl Oncol. 2018;11(6):1379–89.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Tong W, Sun D, Wang Q, Suo J. Sorcin Enhances Metastasis and Promotes Epithelial-to-Mesenchymal Transition of Colorectal Cancer. Cell Biochem Biophys. 2015;72(2):453–9.

    Article  CAS  PubMed  Google Scholar 

  46. Liu Y, Li H, Liu F, Gao LB, Han R, Chen C, Ding X, Li S, Lu K, Yang L, Tian HM, Chen BB, Li X, Xu DH, Deng XL, Shi SL. Heterogeneous nuclear ribonucleoprotein A2/B1 is a negative regulator of human breast cancer metastasis by maintaining the balance of multiple genes and pathways. EBioMedicine. 2020;51: 102583.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Zhu F, Yang T, Yao M, Shen T, Fang C. HNRNPA2B1, as a m6A Reader, Promotes Tumorigenesis and Metastasis of Oral Squamous Cell Carcinoma. Front Oncol. 2021;11: 716921.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Li K, Chen J, Lou X, Li Y, Qian B, Xu D, et al. HNRNPA2B1 Affects the Prognosis of Esophageal Cancer by Regulating the miR-17-92 Cluster. Front Cell Dev Biol. 2021;9: 658642.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Karimi Dermani F, Gholamzadeh Khoei S, Afshar S, Amini R. The potential role of nucleophosmin (NPM1) in the development of cancer. J Cell Physiol. 2021;236(11):7832–52.

  50. Peng HH, Ko HH, Chi NC, Wang YP, Lee HC, Pan PY, et al. Upregulated NPM1 is an independent biomarker to predict progression and prognosis of oral squamous cell carcinomas in Taiwan. Head Neck. 2020;42(1):5–13.

    Article  PubMed  Google Scholar 

  51. Liu Y, Zhang F, Zhang XF, Qi LS, Yang L, Guo H, et al. Expression of nucleophosmin/NPM1 correlates with migration and invasiveness of colon cancer cells. J Biomed Sci. 2012;19(1):53.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Zhang J, Zhao HL, He JF, Li HY. Inhibitory effect of NSC348884, a small molecular inhibitor of nucleophosmin, on the growth of hepatocellular carcinoma cell line hepG2. Zhongguo Yi Xue Ke Xue Yuan Xue Bao. 2012;34(1):58–61.

    PubMed  Google Scholar 

  53. Qi W, Shakalya K, Stejskal A, Goldman A, Beeck S, Cooke L, et al. NSC348884, a nucleophosmin inhibitor disrupts oligomer formation and induces apoptosis in human cancer cells. Oncogene. 2008;27(30):4210–20.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

The work reported here was financially supported by Science and Engineering Research Board (SERB) (Project ID- CRG/2020/002100) and Indian Council of Medical Research (ICMR), Govt. of India, New Delhi. Vaishali Jain is a Ph.D. student registered at recipient of Senior Research Fellowship (SRF) from the CSIR, Govt. of India. Javed Akhtar is a Ph.D. student registered at Jamia Hamdard, New Delhi and a recipient of Senior Research Fellowship (SRF) from the Council of Scientific and Industrial Research (CSIR) and ICMR- National Institute of Pathology, Govt. of India. Ratna Priya is a Ph.D. student registered at Jamia Hamdard, New Delhi and DST- Innovation in Science Pursuit for Inspired Research (INSPIRE)- Senior Research Fellow (SRF). Vivek Ghose is SRF from Department of Biotechnology (DBT), Govt. of India.

Author information

Authors and Affiliations

Authors

Contributions

PG, PS designed the experiment. PS, SG, AKA, VJ, JA, and RP contributed to clinical sample collection and clinical data management. VJ and JA did experimental work and data compilation. Mass spectrometric data acquisition and analysis- VG, VJ, JA. Literature search and interpretation of overall data was done by PG, RS, RVP, PS, AKA, FS, VJ, JA, RP. Drafting and editing of the manuscript was done by PG, VJ, PS, SG, AKA, RS, RVP, PS, VG. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Puja Sakhuja or Poonam Gautam.

Ethics declarations

Ethics approval and consent to participate

All the methods were performed in accordance with the Declaration of Helsinki guidelines and regulations. Clinical samples from participants visiting GIPMER, Delhi, were collected for the study after approval from the Institutional Human Ethics Committee [F.1/IEC/MAMC/80/08/2020/No. 314) and NIP-IEC/10–12-19/06]. All the participants provided informed consent to participate in the study and written informed consent was obtained.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table S1.

Clinical parameters of patients such as age, gender, TNM, stage, grade, white cell count, liver enzymes (AST/ ALT/ ALP), bilirubin and co-morbidities (jaundice, pulmonary tuberculosis, asthma, diabetes melitus, hypertension (HTN), loss of appetite(LOA) and LOW- Loss of weight; Foot note: NA- Not applicable, (-) Not available. Supplementary Table S2. List of 58 DEPs in LN Positive GBC. The proteins with ≥2.0 fold change and adjusted p value <0.05 were used as DEPs. The details are provided in the Methods section. The Table describes the total number of proteins and peptides identified in mass spectrometry runs along with their quantity values. Supplementary Table S3. Bioinformatic analysis of 58 DEPs in LN positive GBC using STRING databse. (A) Top 10 Molecular Functions (B) Top 10 reactome pathways. Supplementary Table S4. Literature survey of 58 DEPs for their association with LN metastasis in cancer. A total of 6 proteins are reported in GBC, 33 proteins in other cancers and 19 were found to be novel to cancer with respect to LN metastasis.

Additional file 2: Supplementary Figure S1.

Experimental design of iTRAQ-based quantitative proteomic analysis for identification of differentially expressed proteins in LN Metastatic GBC. Supplementary Figure S2. SDS-PAGE profile of tissue lysate GBC Stage IIIA and IIIB and control (GSD). A total of 15 µg protein was resolved on 12% gel and stained with Coomassie Brilliant Blue R250 to visualize the protein bands. The protein load for different samples was normalized based on the total density of proteins in each lane. Supplementary Figure S3. Venn diagram showing proteins identified in four technical replicate runs. We found a total of ~700 proteins in each replicate, of which 468 proteins were identified by all four technical replicates. Venn diagram was prepared using the BioVenn software. Supplementary Figure S4. Volcano plot showing DEPS in individual GBC patients. The volcano map was prepared by using log2 fold change and -log10 (p-value) as the co-ordinates and significant fold change ≥ 2.0 and p-value <0.05 were considered to screen the proteins. Dots in orange, blue and grey represents proteins that are upregulated, downregulated and unchanged respectively. GBC- Gallbladder cancer. Supplementary Figure S5. Full-length blot images of Fig. 7 for the expression of KRT7, KRT19, SRI and NPM1 in the individual tissue samples from LN positive GBC, LN negative GBC and GSD cases. For NPM1, the main image includes a) the lane with MW marker from the ‘low exposure’ image and b) other lanes showing NPM1 expression from the blot with ‘high exposure’. The pooled tissue lysate from GSD or LN negative or LN positive GBC was used for Negative control (the blot without primary antibody). Negative control data is not presented in the main image. The cropping of the blot images is indicated with red dashed line.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jain, V., Akhtar, J., Priya, R. et al. Tissue proteome analysis for profiling proteins associated with lymph node metastasis in gallbladder cancer. BMC Cancer 23, 402 (2023). https://doi.org/10.1186/s12885-023-10840-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-10840-3

Keywords