Next Article in Journal
Association between Angiotensin-Converting Enzyme Inhibitors and Lung Cancer—A Nationwide, Population-Based, Propensity Score-Matched Cohort Study
Next Article in Special Issue
Protein Expression in Metastatic Melanoma and the Link to Disease Presentation in a Range of Tumor Phenotypes
Previous Article in Journal
A Disintegrin and Metalloproteinase 9 (ADAM9) in Advanced Hepatocellular Carcinoma and Their Role as a Biomarker During Hepatocellular Carcinoma Immunotherapy
Previous Article in Special Issue
Proteome Profiling Uncovers an Autoimmune Response Signature That Reflects Ovarian Cancer Pathogenesis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Immune-Complexome Analysis Identifies Immunoglobulin-Bound Biomarkers That Predict the Response to Chemotherapy of Pancreatic Cancer Patients

1
Department of Molecular Biotechnology and Health Sciences, University of Turin, 10126 Torino, Italy
2
Center for Experimental Research and Medical Studies (CeRMS), University of Turin, 10126 Torino, Italy
3
Department of Computer Science, University of Turin, 10149 Torino, Italy
4
Department of Clinical and Biological Sciences, San Luigi Hospital, University of Turin, 10043 Orbassano, Turin, Italy
5
MD Anderson Cancer Center, University of Texas, Houston, TX 77030, USA
6
Centro Oncologico Ematologico Subalpino (COES), University of Turin, 10126 Torino, Italy
7
Cittá della salute e della scienza University Hospital of Turin, University of Turin, 10126 Torino, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 27 February 2020 / Revised: 14 March 2020 / Accepted: 18 March 2020 / Published: 21 March 2020
(This article belongs to the Special Issue The Cancer Proteome)

Abstract

:
Pancreatic Ductal Adenocarcinoma (PDA) is an aggressive malignancy with a very poor outcome. Although chemotherapy (CT) treatment has poor efficacy, it can enhance tumor immunogenicity. Tumor-Associated Antigens (TAA) are self-proteins that are overexpressed in tumors that may induce antibody production and can be PDA theranostic targets. However, the prognostic value of TAA-antibody association as Circulating Immune Complexes (CIC) has not yet been elucidated, mainly due to the lack of techniques that lead to their identification. In this study, we show a novel method to separate IgG, IgM, and IgA CIC from sera to use them as prognostic biomarkers of CT response. The PDA Immune-Complexome (IC) was identified using a LTQ-Orbitrap mass spectrometer followed by computational analysis. The analysis of the IC of 37 PDA patients before and after CT revealed differential associated antigens (DAA) for each immunoglobulin class. Our method identified different PDA-specific CIC in patients that were associated with poor prognosis patients. Finally, CIC levels were significantly modified by CT suggesting that they can be used as effective prognostic biomarkers to follow CT response in PDA patients.

1. Introduction

The role of the immune system in cancer has been extensively studied, due to its critical role in counteracting tumor development. The field of cancer immunology has elicited growing interest because of the success obtained in many cancers [1]. Even in tumors defined as immunologically “cold“ (i.e., with a low infiltrate of T cells), such as pancreatic cancer, the possibility to exploit the immune system to identify biomarkers or immunological targets has been demonstrated [2,3,4,5,6,7].
Immune complexes (IC) are those formed between the antibody and its antigen [8]. The self- or non-self-antigens are bound to immunoglobulins by non-covalent interactions. Circulating Immune Complexes (CIC) are present in healthy subjects and cooperate in the immune response [9]. In healthy subjects, IC are removed from the bloodstream, while, in pathological conditions, CIC can accumulate in the bloodstream or in different tissues, especially in the kidneys, inducing glomerulonephritis [8]. Pathological CIC were discovered as biomarkers of Raynaud’s syndrome [8]. Indeed, the accumulation of CIC could provide a novel source of biomarkers for the diagnosis and the follow-up of a pathological condition. The use of CIC as disease biomarkers was specifically studied for autoimmune diseases, such as systemic lupus erythematosus [9,10] and rheumatoid arthritis [8], but it was also observed in other pathologies, including infectious diseases and acquired thrombotic thrombocytopenic purpura [8,11,12]. Of note, the level of CIC increases in the presence of several malignancies. In the 1970s, CIC were found in melanoma, osteogenic sarcoma, lymphoma and carcinomas of the colon, lungs, breast, prostate, pancreas, stomach, esophagus, and uterine cervix [13,14]. However, a systematic analysis of CIC was provided only in 2013 in uveitis, using pull-down assays to separate IgG-IC [15]. One of the main reasons for this delay was the lack of a specific assay to detect and quantify CIC, which led to the missing identification of specific biomarkers in malignancies until now.
With the advent of high-throughput technologies, such as mass spectrometry, the study of CIC has become more accurate. Thus, several techniques have been used to study IC, including ultracentrifugation, ELISA, gel filtration, ultrafiltration, precipitation with polyethylene glycol, assays based on interactions with purified complement portions, quantitative dot-blot assay, and Raji cell radioimmunoassay. Furthermore, CIC were also studied through assays that measured their effects on other cells such as platelet aggregation tests, inhibition of antibody-dependent cytotoxicity assays, tests for the release of enzymes from eosinophils and mast cells, and macrophages inhibition assays [8,9,10,13]. However, most of these techniques allow the investigation of the presence of CIC bounded with already known antigens, without identifying novel CIC-associated antigens.
Immunoglobulin (Ig) G is the main Ig class in human blood, and it represents about 70% of the total Ig pool. For this reason, in studies performed on the humoral response, attention has mainly been focused on IgG class response, including our previous studies on antibody response in PDA patients [2,4,16]. However, the role of other Ig classes in cancer could also be relevant. Indeed, it was demonstrated that IgM and IgG antibodies were bound to allogeneic tumor cells [17]. Conversely, a subpopulation of IgA-secreting B cells that suppresses the activation of CD8 + T cells in prostate cancer, after oxaliplatin treatment, was identified [18]. Moreover, several studies identified IgM CIC containing squamous cell carcinoma antigen or carcinoembryonic antigen, which have diagnostic value in hepatocellular carcinoma and colorectal cancer, respectively [19,20].
Demonstrating the association between CIC levels of cofilin-1 in sera with cancer progression and poor prognosis has highlighted the importance of IC in PDA [21]. It has been well established that chemotherapy (CT) influences the humoral and cellular immune responses to the tumors [22], and, therefore, identifying CIC that can predict responses or resistance to CT could allow the maximization of the efficacy of treatments and could avoid useless toxic effects in non-responding patients [23].
Here, we propose a systematic analysis of IgG, IgM, and IgA CIC in PDA. For the first time, we have demonstrated that, in PDA, identifying novel Tumor-Associated Antigens (TAA) associated with different classes of Ig patients treated with CT can provide new prognostic markers.

2. Results

We performed a systematic analysis of the PDA IC. The CIC-containing different Ig subclasses were sequentially separated and then analyzed by high resolution mass spectrometry (further details in Section 4). Sera from PDA before and after gemcitabine-based CT (namely BCT and ACT, respectively) were collected and analyzed. The IC bound to IgG, IgM, and IgA classes obtained from 37 PDA patients were subdivided into five groups characterized by different disease prognosis. Group one consisted of patients with the best prognosis, while in group five were patients with the worst prognosis (Supplementary Table S1A and Figure 1A). In each group, sera were pooled and the pooled IC content was analyzed using an LTQ-Orbitrap mass spectrometer (Figure 1A). For detailed information, see the Materials and Methods section). In order to identify CIC markers that could be informative for PDA prognosis and CT response, we analyzed the pooled sera of each BCT and ACT group. In this analysis, we identified 3547 proteins bound to IgG, 1043 proteins bound to IgM, and 1143 proteins bound to IgA before and after CT (Supplementary Table S1B–D). For each class of immunoglobulins, we performed a three-step computational analysis. Firstly, to investigate the effect of CT on the CIC content, the presence of each antigen before and after CT in each group was analyzed. Then, the different levels of CIC measured in the five BCT and ACT groups were analyzed to identify possible trends (i.e., a CIC that is found to increase or decrease progressively from group one to group five). Finally, the ACT/BCT ratio for each CIC in each group of patients was considered.

2.1. Analysis of IgG-Circulating Immune Complexes

The analysis of IgG-bound proteins BCT and ACT revealed that only the TTR (thyroid hormone-binding protein) CIC was significantly down-modulated by gemcitabine treatment (Figure 1B, p = 9.64 × 10 4 , Table 1 and Supplementary Table S2A). Considering BCT and ACT results separately, the analysis of different CIC levels in the five cohorts showed a statistically significant trend for IGKV6D-21 (Immunoglobulin kappa variable 6D-21) and IGLV7-46 (Immunoglobulin lambda variable 7-46) CIC in BCT data (Table 2). Conversely, IgG CIC of APOC2 (Apolipoprotein C-II), C4B (Complement C4-B), and S100A9 (Protein S100-A9) increased significantly from group one to group five in ACT data (p = 0.0274 and p = 0.043, respectively) (Table 3). Interestingly, results from the trend analysis demonstrated that IGLV7-46 (BCT) and S100A9 (ACT) showed the most significant trend (p = 0.027) as they were highly abundant in the low survival patient group and absent in the high survival patient blood (Figure 1C, Supplementary Table S3A). Trend analysis of CIC level BCT w.r.t. ACT (ACT/BCT ratio) in the five groups revealed 79 differentially recognized proteins (Supplementary Table S4A and Table S5A), with IGHV4-34 (immunoglobulin heavy variable 4-34) displaying a progressively increasing abundant trend (ratio ACT/BCT: 1, 0.57, 0.07, 60.2, and 119.3 for groups 1, 2, 3, 4, and 5, respectively).

2.2. Analysis of IgM Circulating Immune Complexes

The analysis of IgM-bound proteins BCT and ACT revealed the statistically significant variation of eight different protein levels, namely ADIPOQ (adiponectin), CDH5 (Cadherin-5), FCGBP (IgGFc-binding protein), SERPINA3 (alpha-1-antichymotrypsin), NEO1 (Neogenin) and SERPIND1 (Heparin cofactor 2), which were increased ACT, and HP (Haptoglobin) and KRT16 (Keratin, type I cytoskeletal 16), which were decreased ACT (Figure 1D, Table 1 and Supplementary Table S2B). Analyzing the IgM CIC levels in the five cohorts, considering BCT and ACT separately, revealed that seven proteins were significantly modulated w.r.t. patient survival: AOC3 (membrane primary amine oxidase), APCS (serum amyloid P-component), APOB (apolipoprotein B-100) and IGLV3-25 (immunoglobulin lambda variable 3-25) in the BCT data (Table 2). GAPDH (glyceraldehyde-3-phosphate dehydrogenase), IGLV1-47 (immunoglobulin lambda variable 1-47), and OLFM2 (noelin-2) were increased significantly in the ACT data (Table 3). CIC containing AOC3, APOB, IGLV3-25, IGLV1-47, and OLFM2 decreased from group 1 to group 5 while those with APCS and GAPDH increased in abundance (Figure 1E, Table 2 and Table 3 and Supplementary Table S3B). Comparing the ACT/BCT ratios between the five groups, 30 proteins were identified as being significantly modulated (Supplementary Tables S4B and S5B). Of such proteins, ADIPOQ, AOC3, and SERPINA1 showed the most significant trends (p-value < 0.0001), and their level progressively increased from group 1 to group 5 (Supplementary Table S4B).

2.3. Analysis of IgA Circulating Immune Complexes

The analysis of IgA-bound proteins BCT and ACT revealed the significant modulation of seven proteins: ABCA13 (ATP Binding Cassette Subfamily A Member 13), IGLV1-40 (Immunoglobulin Lambda Variable 1-40) and XIRP2 (Xin Actin Binding Repeat Containing 2) which increased ACT and ECM1 (Extracellular Matrix Protein 1), IGHG3 (Immunoglobulin Heavy Constant Gamma 3), IGKV3-20 (Immunoglobulin Kappa Variable 3-20) and IGLV7-43 (Immunoglobulin Lambda Variable 7-43) which decreased ACT (Figure 1F, Table 3 and Supplementary Table S2C). Considering results BCT and ACT separately, the analysis of CIC levels in the five cohorts showed a statistically significant trend for 10 proteins: APOL1 (apolipoprotein L1), IGHV3OR16-9 (immunoglobulin heavy variable 3/OR16-9), MMRN2 (multimerin-2), SERPINC1 (antithrombin-III), and SERPIND1 (heparin cofactor 2) in the BCT condition (Table 2), while C9 (complement component C9), FREM2 (FRAS1-related extracellular matrix protein 2), IGHD (immunoglobulin heavy constant delta), KCNQ2 (potassium voltage-gated channel subfamily KQT member 2), and PTGDS (prostaglandin-H2 D-isomerase) were detected in the ACT condition (Table 3 and Supplementary Table S3C). Interestingly, all these CIC, increased from group 1 to group 5, except for SERPINC1, which decreased in abundance (Figure 1G and Table 2). In ACT CIC containing PTGDS progressively decreased from group 1 to group 5 (Figure 1G and Table 3). Comparing the modulation of CT in the trends of CIC BCT and ACT in the five groups, 47 proteins were identified as statistically significantly related to patient survival (Supplementary Tables S4C and S5C). PTGDS, KCNQ2, PDE4A, MAN1C1, IGHV3-23, OBSNC, and BRSK2 were characterized as having the most significant trend that increased from group 1 to group 5 while MRPS26 decreased significantly from group 1 to group 5.

2.4. Enrichment Analysis of IgA, IgM and IgG CIC

To predict candidate biological processes and signaling pathways that are enriched in the set of antigens bound to the CIC, we performed an enrichment analysis of combined IgA, IgM, and IgG DAA. This revealed different enriched processes and phenotypes including “PDA genesis“, “type II diabetes mellitus onset“, and “extracellular vesicles“ (Supplementary Table S6). PDA genesis, insulin secretion, and diabetes pathways confirm the accuracy of our analysis.

3. Discussion

The relationship between CIC and tumors has been known for many years, but a full and complete comprehension of the role of CIC in cancer is still lacking. The observation of elevated levels of CIC in the presence of malignancies dates back to the 1970s [13,14]. However, the study of CIC remained a challenge until the development of mass spectrometry. In fact, in the past, CIC were studied with methods that separated them only roughly, and assumed for the analysis the prior knowledge of the antigen under study, thus not allowing the discovery of new antigens [8,9,10,13]. Instead, the combination of a proper CIC separation and mass spectrometry as presented here leads to the identification of the proteins involved in the disease, opening new perspectives to investigate the role of these CIC in cancer.
In this article, we have proposed a novel approach to investigate CIC that has led us to identifying important TAA suitable to be follow-up biomarkers. The analysis was performed on the sera of patients with PDA before and after gemcitabine-based CT. Currently, we are recruiting a large number of PDA patients in a multi-center clinical study in order to confirm and validate the prognostic role of CIC-bound TAA identified in this study.
In contrast to previously used methods that allowed the separation of all immunoglobulin repertoire (ultracentrifugation, gel filtration, ultrafiltration, precipitation with polyethylene glycol) [8] or the enrichment of CIC containing IgG only [15,16], we presented a systematic analysis of CIC containing IgG, IgM, and IgA isotypes that we define as the “immune-complexome“. This represents a technical and conceptual improvement that was found of particular interest. In fact, we demonstrated that the proteins bound to different classes of immunoglobulins, and their changes in relation to CT, defined different patterns of response to CT, underlying the need to consider them separately. Of interest, proteins bounded to immunoglobulins also showed marked differences with respect to serum-free circulating proteins [24].
The most interesting outcome of this study is that the clinically relevant CIC trends were from the IgM fraction. It is known that IgM CIC are more frequent w.r.t. IgG ones because the clearance of IgG CIC is maintained by (i) the activation of complement by a classical pathway and (ii) the Fc gamma receptors (Fc γ R) system which is present on the surface of different cells of the immune system. This second modality of CIC clearance is not employed for IgM CIC. For this reason, it is possible that proteins bound to IgM reflect the earliest humoral response to PDA after CT. Of note, the IgG fraction is not the most important class associated with protein markers in our analysis; although sera IgG were the most present quantitatively, they showed less informative cues when the purpose was to identify of IgG-associated biomarkers following CT treatment.
Enrichment analysis strongly reinforced the robustness of our study, highlighting many processes related to pancreatic cancer (PDA genesis, type II diabetes mellitus onset, and insulin secretion). Notably, the enrichment in protein localized in the extracellular vesicle was particularly interesting, since a significant enrichment of exosome in the Ig-bound fraction was demonstrated in the sera of PDA patients compared to healthy subjects [16]. Moreover, several identified antigens have been previously correlated to PDA, although not in association with immunoglobulins, namely TTR (thyroid hormone-binding protein) and S100A9 (protein S100-A9) were found in IgG CIC [25,26], ADIPOQ (adiponectin), CDH5 (cadherin-5), FCGBP (IgGFc-binding protein), HP (haptoglobin), GAPDH (glyceraldehyde-3-phosphate dehydrogenase), and SERPINA1 (serine proteinase inhibitor A1) in IgM CIC [24,27,28,29], BRSK2 (BR Serine/Threonine-Protein Kinase 2) in IgA CIC.
TTR and HP bound to fucose have been already proposed as a biomarkers of PDA [30,31]; CDH5 in pancreatic neuroendocrine tumors developed in Von Hippel–Lindau disease patients [32] and FCGBP in mucinous cysts in the pancreas [33] and PDA [34]. S100A9 is a calcium- and zinc-binding protein with a key function in regulating inflammation and immune responses, and has been previously associated with PDA inflammation status, invasiveness, and metastasis [25,35]; in our experimental setting, S100A9 levels increased in case of worst prognosis. ADIPOQ levels were correlated to various kinds of cancers [36], but its role has not been well understood, yet. Indeed, even though the preponderance of evidence suggests an inverse correlation of serum ADIPOQ levels with cancer, other studies correlated the increase of ADIPOQ levels with cancer progression [36]. In the case of PDA, the scenario is particularly disputable. Adiponectin inhibits PDA proliferation of Pan02 cells both in vitro and in vivo when cells were transplanted in a murine model, by induction of caspase-3 and caspase-7 and apoptosis [37]. Conversely, in the same year, it was showed that ADIPOQ inhibited apoptosis of pancreatic cancer cells in the same murine model via activation of AMPK-Sirt1-PGC1 α signaling [27]. In human studies, the situation is controversial; adiponectin in pancreatic cancer patients has been separately observed to be elevated [38,39,40], decreased [41,42], and unchanged [43,44]. A possible explanation of these conflicting results could be that, as shown in this study, ADIPOQ is not important by itself, but the ADIPOQ bound to IgM. In fact, we showed a very interesting trend, with IgM-ADIPOQ levels significantly increasing with worse prognoses.
GAPDH, which we observed to increase in being bound to IgM from the group one to the group five ACT, has already been correlated to PDA with a worse prognosis [24]. SERPINA1 is up-regulated in PDA tissue and serum [45] and has been proposed as a prognostic factor in pancreatic cancer [46]. BRSK2, which we saw to increase from group one to group five, was demonstrated to provide a survival advantage to PDA cells [47].
Other identified antigens have been correlated to other types of cancer, namely ABCA13 in ovarian cancer [48], XIRP2 in gastric cancer [49], and APOC2 expression in Non-Small Cell Lung Cancer [50]. Typical serum proteins were among the identified ones, despite the immune depletion before the mass spectrometry analysis. We found complement proteins, immunoglobulins, apolipoproteins, and inflammation proteins, some of which were significantly modulated by CT, namely APOC2 (apolipoprotein C-II), C4B (complement C4-B), IGLV7-46, IGHV4-34 (immunoglobulin heavy variable 4-34) in IgG CIC, APOB (apolipoprotein B-100), IGLV3-25 (immunoglobulin lambda variable 3-25), IGLV1-47 (immunoglobulin lambda variable 1-47) and SERPINA1 in IgM CIC, APOL1 (apolipoprotein L1), C9 (complement component C9), IGHG3, IGKV3-20, IGLV1-40, IGLV7-43, IGHV3OR16-9 (immunoglobulin heavy variable 3/OR16-9), IGHD (immunoglobulin heavy constant delta), IGHV3-23, SERPINC1 (antithrombin-III) and SERPIND1 (heparin cofactor 2) in IgA CIC. Concerning the role of CIC in cancer, very little data are available. It has been reported that proteins bound to Ig could have a decoy function, inhibiting antibody-dependent cytotoxicity [16]. This aspect could be relevant to explain the increase of the majority of the identified CIC after CT in parallel with prognosis getting worst; however, it is certainly not the only role. It is likely that CIC binding could enhance the immune response against a protein, but it could also sequester a protein rendering it ineffective. Depending on the role of the antigen, this behavior could lead to a good or bad prognosis.

4. Materials and Methods

4.1. Patients

From 2006 to 2012, 37 patients with PDA, enrolled in ENOAPA project (https://www.epiclin.it/enoapa), provided serum samples. The protocol was approved by the local research ethical committee (Azienda Ospedaliera Città della Salute e della Scienza di Torino, Turin) ethical code #0114454, November 23th, 2016.) and was performed according to the Helsinki Declaration principles. All participants in the protocol signed a declaration of informed consent. No patients were subjected to surgery and all were treated with gemcitabine-based CT (gemcitabine with oxaliplatin or alone). Serum samples were isolated from venous blood before and after CT and stored at −80 C until use. Patients’ characteristics are listed in Supplementary Table S1A. Patients were grouped into five groups on the basis of response to therapy and survival, as indicated in Supplementary Table S2. An equal amount of serum from each patient in each group was taken to reach a total volume of 300 μ L.

4.2. IgG-IC Purification

HiTrap Protein G HP (1 mL) was used (GE HealthCare Life Sciences, Milan, Italy). The column was operated with a peristaltic pump. All buffers were filtered before use. Serum was filtered with a SPARTAN Syringe Filter, RC Filter Media, 13 mm diameter, 0.2 μ m pore size (Whatman, Maidstone, UK) and then diluted at 1:3 in binding buffer (20 mM sodium phosphate, pH 7.0). The column was washed with 5 mL of milliQ water and with 5 mL of binding buffer before sample application. The sample was then passed through the column along with 9 mL of binding buffer, collecting the “unbound fraction“ exiting from the column. IgG-IC were eluted from the column with 10 mL of elution buffer (0.1 M glycine-HCl, pH 2.7) and collected in a tube containing 1 mL of neutralizing buffer (1 M Tris-HCl, pH 9.0). All timing and volumes were optimized evaluating IgG presence by Western blotting.

4.3. IgM-IC Purification

After the evaluation of their presence by Western blotting, IgM-IC were separated from the IgG unbound fraction by HiTrap IgM purification HP 1 mL (GE HealthCare Life Sciences, Milan, Italy). The column was operated with a peristaltic pump. The sample should have the same concentration of ammonium sulfate as the binding buffer, so small amounts of solid ammonium sulfate were gradually added to the sample until the final concentration was 1M, stirring slowly and continuously, to avoid IgM precipitation. Samples were then filtered with SPARTAN Syringe Filter, RC Filter Media, 13 mm diameter, 0.2 μ m pore size (Whatman, Maidstone, UK). All buffers were filtered before use. Columns were washed with 5 mL of binding buffer (20 mM sodium phosphate, 1 M ammonium, pH 7.5), 5 mL of elution buffer (20 mM sodium phosphate, pH 7.5) and 5 mL of regeneration buffer (20 mM sodium phosphate, pH 7.5 with 30% isopropanol) and then equilibrated with 5 mL of binding buffer. This was then applied to the sample. The unbound fraction was washed out with 15 mL of binding buffer. Elution was performed with 12 mL of elution buffer. Columns were then regenerated with regeneration buffer and equilibrated with a binding buffer for the subsequent separation. All timings and volumes were optimized evaluating IgM presence by Western blotting.

4.4. IgA-IC Separation

IgA-IC were present with the eluted IgM-IC fraction, as verified by Western blotting. Indeed, starting from the elution containing IgM-IC, IgA-IC were separated by immunoprecipitation with anti-human IgA (a chain specific)-agarose antibody (Sigma Aldrich, St. Louis, MO, USA). After adding 1.9 mL of antibody, samples were incubated overnight at 4 C on a rocking platform. After three washes with PBS, IgA-IC were eluted with two elutions of 4 mL each of elution buffer (0.1 M glycine-HCl, pH 2.7). Elutions were collected in a tube containing 0.8 mL of neutralizing buffer (1 M Tris-HCl, pH 9.0).

4.5. LC-MS Analysis

All samples were concentrated to a volume of 500 μ L using a 3KDa cut-off concentrator (Merck Millipore, Burlington, MA, USA) at 4000 g for 60 min. They were then reduced with TCEP (35 μ L in each sample; incubation at 37 C for 30 min) and alkylated with acrylamide (20 mg in each sample; incubation at room temperature for 30 min in the dark). After buffer exchange with a solution containing 20 mM Tris, 4 M urea, 3% isopropanol, pH 8.0, samples were subjected to reverse phase separation as previously described [51]. Samples were then fractionated in protein level using the Shimadzu HPLC system. The protein sample was loaded onto the RPGS reversed-phase column (4.6 mm I.D. ×150 mm, 15 μ m, 1000Å, Column Technology Inc., #RP5D46150) and desalted for 5 min with 95% mobile-phase A (0.1% TFA in 95% H 2 O) at a flow rate of 3 mL/min. Proteins were eluted from the column at a flow rate of 2.1 mL/min with a gradient elution including: (i) an increase from 5 to 70% mobile phase B (0.1% TFA in 95% ACN) for more than 25 min, (ii) an increase from 70% to 95% mobile phase B for 3 min, (iii) a wash step to hold at 95% mobile phase B for 2 min, (iv) a re-equilibration step at 95% mobile phase A for 5 min. The collected fractions were pooled into 15 fractions. All fractions obtained in this way were lyophilized and proteins were digested with trypsin. The analysis was performed using Orbitrap Elite. The tryptic peptides were separated by reversed-phase chromatography using an EASYnano HPLC system (Thermo Scientific, Waltham, MA, USA) coupled online with a LTQ-Orbitrap ELITE mass spectrometer (Thermo Scientific, Waltham, MA, USA). Mass spectrometer parameters were: spray voltage 2.5 kV, capillary temperature 280 C, Furier Transform (FT) resolution 60,000, FT target value 1 × 10 6 , LTQ target value 3 × 10 4 , 1 FT microscan with 500 ms injection time, and 1 LTQ microscan with 10 ms injection time. Mass spectra were acquired in a data-dependent mode with the m/z range of 400–2000. The full mass spectrum (MS scan) was acquired by the FT, while the tandem mass spectrum (MS/MS scan) was acquired by the LTQ with 35% normalized collision energy. Acquisition of each full mass spectrum was followed by the acquisition of MS/MS spectra for the five most intense +2 or +3 ions within a one second duty cycle. The minimum signal threshold (counts) for a precursor occurring during a MS scan was set at 1000 for triggering an MS/MS scan.
The acquired LC-MS/MS data were processed by the Proteome Discoverer 1.4 (Thermo Scientific, Waltham, MA, USA). Sequest HT was used as a search engine with input parameters including cysteine (Cys) alkylated with acrylamide (71.03714@C) as a fixed modification and methionine (Met) oxidation (15.99491@M) as a variable modification. Data were searched against the Uniprot Human database 2017 and further filtered with a False Discovery Rate (FDR) < = 5%. Total mass spectra counts for each protein were used as a measure of protein abundance.

4.6. Protein Identification and Bioinformatics Aanalysis

A new R computational pipeline to pre-process, normalize, and analyze obtained proteomic data was applied. Our pipeline consisted of data characterization, individuation of differentially associated antigen (DAA), trend analysis, and functional enrichment analysis. Using our approach, we identified DAA and their trends upon CT. We also characterized the enriched biological processes and phenotypes of the antigens identified.

4.7. Data Pre-Processing and Normalization

Spectra obtained from LTQ-Orbitrap were analyzed using SEQUEST software. SEQUEST identified peptide sequences from tandem mass spectra. Peptides-Spectrum Matches (PSMs) obtained from this analysis were filtered in order to take into consideration only IC with PSM > 1 or IC with PSM equal to 1 occurring in at least two groups. Filtered PSMs were normalized using the quantile normalization method in R.

4.8. Differential Associated Antigens and Antigen Recognition Trends

Comparisons of IC abundance (before and after CT) were computed using the DESeq2 R package [52]. Analysis of IC abundances’ trends among the five sample pools was performed using MaSigPro R package [53] with default parameters.

4.9. Gene Set Enrichment Analysis

Gene list coding for identified proteins was retrieved using UniProt and BioDBnet and analyzed using an Enrichr tool [54].

5. Conclusions

In conclusion, our method leads: i) to the identification of novel PDA specific CIC before and after gemcitabine treatment, ii) to the characterization of CT effects on CIC levels. In this study, we developed a new method that combines the capture of IC and mass spectrometry, allowing the discovery that CIC-bound TAAs could be markers of CT response. Some of the identified candidate CIC could be suitable for use as novel non-invasive prognostic markers.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6694/12/3/746/s1, Table S1: Clinical features of patients and list of identified proteins. A. Table containing the clinical information of patients enrolled in the study B. List of proteins bound to IgG identified in the study C. List of proteins bound to IgM identified in the study D. List of proteins bound to IgA identified in the study, Table S2: Differential recognized antigens. A. List of IgG antigens significantly modulated by gemcitabine B. List of IgM antigens significantly modulated by gemcitabine C. List of IgA antigens significantly modulated by gemcitabine, Table S3: Significant Trends before and after CT. A. List of IgG CICs trends BCT and ACT with p < 0.05 B. List of IgM CICs trends BCT and ACT with p < 0.05 C. List of IgA CICs trends BCT and ACT with p < 0.05, Table S4: MaSig Pro Results of significant trends upon CT treatment. A. IgG significant trends computed by MaSig Pro. Only antigens recognized by at least three out of five groups of patients are represented. B. IgM significant trends computed by MaSig Pro. Only antigens recognized by at least three out of five groups of patients are represented. C. IgA significant trends computed by MaSig Pro. Only antigens recognized by at least three out of five groups of patients are represented, Table S5: ACT/BCT ratios of the different groups of patients. A. IgG ACT/BCT ratio of antigens with significant trends computed with MaSigPro. B. IgM ACT/BCT ratio of antigens with significant trends computed with MaSigPro. C. IgA ACT/BCT ratio of antigens with significant trends computed with MaSigPro, Table S6: EnrichR results. A. Analysis of combined IgA, IgM and IgG differential recognized antigens.

Author Contributions

Project design, G.M. and L.F.; proteomic analyses, G.M.; methodology, G.M. and L.F.; bioinformatic analyses, L.F.; writing—original draft preparation, G.M., L.F., and G.F.; writing—review and editing, G.M., L.F., G.F., and F.N.; proteomic analyses support H.K., W.H., and A.A.M.; statistical analyses of patients survival, A.E.; Sera collection, D.G., R.S., and M.A.S.; supervision of the proteomic data analysis, M.C. and S.M.H.; supervision of the bioinformatic analyses, F.C.; funding acquisition, F.N.; project administration, F.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the Associazione Italiana Ricerca sul Cancro (5 x mille no. 12182, IG no. 15257); the Italian Ministry of Health-Progetti Ricerca Finalizzata (RF-2013-02354892), the Fondazione Ricerca Molinette Onlus, the Fondazione Nadia Valsecchi awarded to FN; and by the National Institutes of Health of USA (Grant No. R01GM088342) awarded to SMH.

Conflicts of Interest

The authors have declared that there are no competing interests.

Abbreviations

The following abbreviations are used in this manuscript:
PDAPancreatic Ductal Adenocarcinoma
CTChemoTherapy
ICImmune-Complexome
CICCirculating Immune Complexes
DAADifferential Associated Antigens
TAATumor Associated Antigens

References

  1. Soares, K.C.; Zheng, L.; Edil, B.; Jaffee, E.M. Vaccines for pancreatic cancer. Cancer J. Sudbury Mass. 2012, 18, 642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Tomaino, B.; Cappello, P.; Capello, M.; Fredolini, C.; Ponzetto, A.; Novarino, A.; Ciuffreda, L.; Bertetto, O.; De Angelis, C.; Gaia, E.; et al. Autoantibody signature in human ductal pancreatic adenocarcinoma. J. Proteome Res. 2007, 6, 4025–4031. [Google Scholar] [CrossRef] [PubMed]
  3. Cappello, P.; Tomaino, B.; Chiarle, R.; Ceruti, P.; Novarino, A.; Castagnoli, C.; Migliorini, P.; Perconti, G.; Giallongo, A.; Milella, M.; et al. An integrated humoral and cellular response is elicited in pancreatic cancer by α-enolase, a novel pancreatic ductal adenocarcinoma-associated antigen. Int. J. Cancer 2009, 125, 639–648. [Google Scholar] [CrossRef] [PubMed]
  4. Tomaino, B.; Cappello, P.; Capello, M.; Fredolini, C.; Sperduti, I.; Migliorini, P.; Salacone, P.; Novarino, A.; Giacobino, A.; Ciuffreda, L.; et al. Circulating autoantibodies to phosphorylated α-enolase are a hallmark of pancreatic cancer. J. Proteome Res. 2010, 10, 105–112. [Google Scholar] [CrossRef] [PubMed]
  5. Cappello, P.; Rolla, S.; Chiarle, R.; Principe, M.; Cavallo, F.; Perconti, G.; Feo, S.; Giovarelli, M.; Novelli, F. Vaccination with ENO1 DNA prolongs survival of genetically engineered mice with pancreatic cancer. Gastroenterology 2013, 144, 1098–1106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Capello, M.; Cappello, P.; Linty, F.C.; Chiarle, R.; Sperduti, I.; Novarino, A.; Salacone, P.; Mandili, G.; Naccarati, A.; Sacerdote, C.; et al. Autoantibodies to Ezrin are an early sign of pancreatic cancer in humans and in genetically engineered mouse models. J. Hematol. Oncol. 2013, 6, 67. [Google Scholar] [CrossRef] [Green Version]
  7. Cappello, P.; Curcio, C.; Mandili, G.; Roux, C.; Bulfamante, S.; Novelli, F. Next, Generation Immunotherapy for Pancreatic Cancer: DNA Vaccination is Seeking New Combo Partners. Cancers 2018, 10, 51. [Google Scholar] [CrossRef] [Green Version]
  8. Ritzmann, S.E.; Daniels, J.C. Immune complexes: Characteristics, clinical correlations, and interpretive approaches in the clinical laboratory. Clin. Chem. 1982, 28, 1259–1271. [Google Scholar] [CrossRef]
  9. Nezlin, R. A quantitative approach to the determination of antigen in immune complexes. J. Immunol. Methods 2000, 237, 1–17. [Google Scholar] [CrossRef]
  10. Wener, M.H. Tests for circulating immune complexes. In Systemic Lupus Erythematosus; Springer: Berlin/Heidelberg, Germany, 2014; pp. 47–57. [Google Scholar]
  11. Ferrari, S.; Palavra, K.; Gruber, B.; Hovinga, J.A.K.; Knöbl, P.; Caron, C.; Cromwell, C.; Aledort, L.; Plaimauer, B.; Turecek, P.L.; et al. Persistence of circulating ADAMTS13-specific immune complexes in patients with acquired thrombotic thrombocytopenic purpura. Haematologica 2014, 99, 779–787. [Google Scholar] [CrossRef] [Green Version]
  12. Korolevskaya, L.; Shmagel, K.; Saidakova, E.; Shmagel, N.; Chereshnev, V. Immune complexes that contain HIV antigens activate peripheral blood T cells. In Doklady Biological Sciences; Springer: Berlin/Heidelberg, Germany, 2016; Volume 469, pp. 196–198. [Google Scholar]
  13. Theofilopoulos, A.N.; Andrews, B.S.; Urist, M.M.; Morton, D.L.; Dixon, F.J. The nature of immune complexes in human cancer sera. J. Immunol. 1977, 119, 657–663. [Google Scholar] [PubMed]
  14. Seth, P.; Balachandran, N.; Malaviya, A.; Kumar, R. Circulating immune complexes in carcinoma of uterine cervix. Clin. Exp. Immunol. 1979, 38, 77. [Google Scholar] [PubMed]
  15. Merl, J.; Deeg, C.A.; Swadzba, M.E.; Ueffing, M.; Hauck, S.M. Identification of autoantigens in body fluids by combining pull-downs and organic precipitations of intact immune complexes with quantitative label-free mass spectrometry. J. Proteome Res. 2013, 12, 5656–5665. [Google Scholar] [CrossRef] [PubMed]
  16. Capello, M.; Vykoukal, J.V.; Katayama, H.; Bantis, L.E.; Wang, H.; Kundnani, D.L.; Aguilar-Bonavides, C.; Aguilar, M.; Tripathi, S.C.; Dhillon, D.S.; et al. Exosomes harbor B cell targets in pancreatic adenocarcinoma and exert decoy function against complement-mediated cytotoxicity. Nat. Commun. 2019, 10, 254. [Google Scholar] [CrossRef]
  17. Carmi, Y.; Spitzer, M.H.; Linde, I.L.; Burt, B.M.; Prestwood, T.R.; Perlman, N.; Davidson, M.G.; Kenkel, J.A.; Segal, E.; Pusapati, G.V.; et al. Allogeneic IgG combined with dendritic cell stimuli induce antitumour T-cell immunity. Nature 2015, 521, 99. [Google Scholar] [CrossRef]
  18. Shalapour, S.; Karin, M. Immunity, inflammation, and cancer: An eternal fight between good and evil. J. Clin. Invest. 2015, 125, 3347–3355. [Google Scholar] [CrossRef] [Green Version]
  19. Castaldi, F.; Marino, M.; Beneduce, L.; Belluco, C.; De Marchi, F.; Mammano, E.; Nitti, D.; Lise, M.; Fassina, G. Detection of circulating CEA-IgM complexes in early stage colorectal cancer. Int. J. Biol. Markers 2005, 20, 204–208. [Google Scholar] [CrossRef]
  20. Beneduce, L.; Prayer-Galetti, T.; Giustinian, A.M.G.; Gallotta, A.; Betto, G.; Pagano, F.; Fassina, G. Detection of prostate-specific antigen coupled to immunoglobulin M in prostate cancer patients. Cancer Detect. Prev. 2007, 31, 402–407. [Google Scholar] [CrossRef]
  21. Satoh, M.; Takano, S.; Sogawa, K.; Noda, K.; Yoshitomi, H.; Ishibashi, M.; Mogushi, K.; Takizawa, H.; Otsuka, M.; Shimizu, H.; et al. Immune-complex level of cofilin-1 in sera is associated with cancer progression and poor prognosis in pancreatic cancer. Cancer Sci. 2017, 108, 795–803. [Google Scholar] [CrossRef]
  22. Bracci, L.; Schiavoni, G.; Sistigu, A.; Belardelli, F. Immune-based mechanisms of cytotoxic chemotherapy: Implications for the design of novel and rationale-based combined treatments against cancer. Cell Death Differ. 2014, 21, 15. [Google Scholar] [CrossRef] [Green Version]
  23. Grasso, C.; Jansen, G.; Giovannetti, E. Drug resistance in pancreatic cancer: Impact of altered energy metabolism. Crit. Rev. Oncol. Hematol. 2017, 114, 139–152. [Google Scholar] [CrossRef] [PubMed]
  24. Follia, L.; Ferrero, G.; Mandili, G.; Beccuti, M.; Giordano, D.; Spadi, R.; Satolli, M.A.; Evangelista, A.; Katayama, H.; Wang, H.; et al. Integrative analysis of novel metabolic subtypes in pancreatic cancer fosters new prognostic biomarkers. Front. Oncol. 2019, 9, 115. [Google Scholar] [CrossRef] [PubMed]
  25. Basso, D.; Bozzato, D.; Padoan, A.; Moz, S.; Zambon, C.F.; Fogar, P.; Greco, E.; Scorzeto, M.; Simonato, F.; Navaglia, F.; et al. Inflammation and pancreatic cancer: Molecular and functional interactions between S100A8, S100A9, NT-S100A8 and TGFβ1. Cell Cmmun. Signal. 2014, 12, 20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Lv, S.; Gao, J.; Zhu, F.; Li, Z.; Gong, Y.; Xu, G.; Ma, L. Transthyretin, identified by proteomics, is overabundant in pancreatic juice from pancreatic carcinoma and originates from pancreatic islets. Diagn. Cytopathol. 2011, 39, 875–881. [Google Scholar] [CrossRef]
  27. Huang, B.; Cheng, X.; Wang, D.; Peng, M.; Xue, Z.; Da, Y.; Zhang, N.; Yao, Z.; Li, M.; Xu, A.; et al. Adiponectin promotes pancreatic cancer progression by inhibiting apoptosis via the activation of AMPK/Sirt1/PGC-1α signaling. Oncotarget 2014, 5, 4732. [Google Scholar] [CrossRef] [Green Version]
  28. Murphy, S.J.; Hart, S.N.; Lima, J.F.; Kipp, B.R.; Klebig, M.; Winters, J.L.; Szabo, C.; Zhang, L.; Eckloff, B.W.; Petersen, G.M.; et al. Genetic alterations associated with progression from pancreatic intraepithelial neoplasia to invasive pancreatic tumor. Gastroenterology 2013, 145, 1098–1109. [Google Scholar] [CrossRef] [Green Version]
  29. Morishita, K.; Ito, N.; Koda, S.; Maeda, M.; Nakayama, K.; Yoshida, K.; Takamatsu, S.; Yamada, M.; Eguchi, H.; Kamada, Y.; et al. Haptoglobin phenotype is a critical factor in the use of fucosylated haptoglobin for pancreatic cancer diagnosis. Clin. Chim. Acta 2018, 487, 84–89. [Google Scholar] [CrossRef]
  30. Chen, J.; Chen, L.J.; Xia, Y.L.; Zhou, H.C.; Yang, R.B.; Wu, W.; Lu, Y.; Hu, L.W.; Zhao, Y. Identification and verification of transthyretin as a potential biomarker for pancreatic ductal adenocarcinoma. J. Cancer Res. Clin. Oncol. 2013, 139, 1117–1127. [Google Scholar] [CrossRef]
  31. Miyoshi, E.; Kamada, Y. Application of glycoscience to the early detection of pancreatic cancer. Cancer Sci. 2016, 107, 1357–1362. [Google Scholar] [CrossRef] [Green Version]
  32. Speisky, D.; Duces, A.; Bièche, I.; Rebours, V.; Hammel, P.; Sauvanet, A.; Richard, S.; Bedossa, P.; Vidaud, M.; Murat, A.; et al. Molecular profiling of pancreatic neuroendocrine tumors in sporadic and Von Hippel-Lindau patients. Clin. Cancer Res. 2012, 18, 2838–2849. [Google Scholar] [CrossRef] [Green Version]
  33. Park, J.; Yun, H.S.; Lee, K.H.; Lee, K.T.; Lee, J.K.; Lee, S.Y. Discovery and validation of biomarkers that distinguish mucinous and nonmucinous pancreatic cysts. Cancer Res. 2015, 75, 3227–3235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Stamp, L.A.; Braxton, D.R.; Wu, J.; Akopian, V.; Hasegawa, K.; Chandrasoma, P.T.; Hawes, S.M.; McLean, C.; Petrovic, L.M.; Wang, K.; et al. The GCTM-5 epitope associated with the mucin-like glycoprotein FCGBP marks progenitor cells in tissues of endodermal origin. Stem Cells 2012, 30, 1999–2009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Nedjadi, T.; Evans, A.; Sheikh, A.; Barerra, L.; Al-Ghamdi, S.; Oldfield, L.; Greenhalf, W.; Neoptolemos, J.P.; Costello, E. S100A8 and S100A9 proteins form part of a paracrine feedback loop between pancreatic cancer cells and monocytes. BMC Cancer 2018, 18, 1255. [Google Scholar] [CrossRef] [PubMed]
  36. Hebbard, L.; Ranscht, B. Multifaceted roles of adiponectin in cancer. Best Pract. Res. Clin. Eendocrinol. Metab. 2014, 28, 59–69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kato, M.; Watabe, K.; Tsujii, M.; Funahashi, T.; Shimomura, I.; Takehara, T. Adiponectin inhibits murine pancreatic cancer growth. Dig. Dis. Sci. 2014, 59, 1192–1196. [Google Scholar] [CrossRef] [PubMed]
  38. Chang, M.C.; Chang, Y.T.; Su, T.C.; Yang, W.S.; Chen, C.L.; Tien, Y.W.; Liang, P.C.; Wei, S.C.; Wong, J.M. Adiponectin as a potential differential marker to distinguish pancreatic cancer and chronic pancreatitis. Pancreas 2007, 35, 16–21. [Google Scholar] [CrossRef] [Green Version]
  39. Dalamaga, M.; Migdalis, I.; Fargnoli, J.L.; Papadavid, E.; Bloom, E.; Mitsiades, N.; Karmaniolas, K.; Pelecanos, N.; Tseleni-Balafouta, S.; Dionyssiou-Asteriou, A.; et al. Pancreatic cancer expresses adiponectin receptors and is associated with hypoleptinemia and hyperadiponectinemia: A case–control study. Cancer Caus. Control 2009, 20, 625–633. [Google Scholar] [CrossRef] [Green Version]
  40. Dranka-Bojarowska, D.; Lekstan, A.; Olakowski, M.; Jablonska, B.; Lewinski, A.; Musialski, P.; Sobczyk, W.; Kapalka, A.; Lampe, P. The assessment of serum concentration of adiponectin, leptin and serum carbohydrate antigen-19.9 in patients with pancreatic cancer and chronic pancreatitis. J. Physiol. Pharmacol. 2015, 66, 653–663. [Google Scholar]
  41. Stolzenberg-Solomon, R.Z.; Weinstein, S.; Pollak, M.; Tao, Y.; Taylor, P.R.; Virtamo, J.; Albanes, D. Prediagnostic adiponectin concentrations and pancreatic cancer risk in male smokers. Am. J. Epidemiol. 2008, 168, 1047–1055. [Google Scholar] [CrossRef]
  42. Bao, Y.; Giovannucci, E.L.; Kraft, P.; Stampfer, M.J.; Ogino, S.; Ma, J.; Buring, J.E.; Sesso, H.D.; Lee, I.M.; Gaziano, J.M.; et al. A prospective study of plasma adiponectin and pancreatic cancer risk in five US cohorts. J. Natl. Cancer Ins. 2012, 105, 95–103. [Google Scholar] [CrossRef]
  43. Pezzilli, R.; Barassi, A.; Corsi, M.M.; Morselli-Labate, A.M.; Campana, D.; Casadei, R.; Santini, D.; Corinaldesi, R.; D’Eril, G.M. Serum leptin, but not adiponectin and receptor for advanced glycation end products, is able to distinguish autoimmune pancreatitis from both chronic pancreatitis and pancreatic neoplasms. Scand. J. Gastroenterol. 2010, 45, 93–99. [Google Scholar] [CrossRef]
  44. Kadri, M.C.; Bostanci, E.; Ozdemir, Y.; Dalgic, T.; Aksoy, E.; Ozer, I.; Ozogul, Y.; Oter, V. Roles of adiponectin and leptin as diagnostic markers in pancreatic cancer. Bratislavske Lekarske Listy 2017, 118, 394–398. [Google Scholar]
  45. Trachte, A.L.; Suthers, S.E.; Lerner, M.R.; Hanas, J.S.; Jupe, E.R.; Sienko, A.E.; Adesina, A.M.; Lightfoot, S.A.; Brackett, D.J.; Postier, R.G. Increased expression of alpha-1-antitrypsin, glutathione S-transferase π and vascular endothelial growth factor in human pancreatic adenocarcinoma. Am. J. Surg. 2002, 184, 642–647. [Google Scholar] [CrossRef]
  46. Trichopoulos, D.; Tzonou, A.; Kalapothaki, V.; Sparos, L.; Kremastinou, T.; Skoutari, M. Alpha1-antitrypsin and survival in pancreatic cancer. Int. J. Cancers 1990, 45, 685–686. [Google Scholar] [CrossRef]
  47. Hexige Saiyin, N.N.; Han, X.; Fang, Y.; Wu, Y.; Lou, W.; Yang, X. BRSK2 induced by nutrient deprivation promotes Akt activity in pancreatic cancer via downregulation of mTOR activity. Oncotarget 2017, 8, 44669. [Google Scholar]
  48. Nymoen, D.A.; Holth, A.; Falkenthal, T.E.H.; Tropé, C.G.; Davidson, B. CIAPIN1 and ABCA13 are markers of poor survival in metastatic ovarian serous carcinoma. Mol. Cancer 2015, 14, 44. [Google Scholar] [CrossRef] [Green Version]
  49. Li, X.; Wu, W.K.; Xing, R.; Wong, S.H.; Liu, Y.; Fang, X.; Zhang, Y.; Wang, M.; Wang, J.; Li, L.; et al. Distinct subtypes of gastric cancer defined by molecular characterization include novel mutational signatures with prognostic capability. Cancer Res. 2016, 76, 1724–1732. [Google Scholar] [CrossRef] [Green Version]
  50. Zeng, X.; Hood, B.L.; Zhao, T.; Conrads, T.P.; Sun, M.; Gopalakrishnan, V.; Grover, H.; Day, R.S.; Weissfeld, J.L.; Wilson, D.O.; et al. Lung cancer serum biomarker discovery using label-free liquid chromatography-tandem mass spectrometry. J. Thor. Oncol. 2011, 6, 725–734. [Google Scholar] [CrossRef] [Green Version]
  51. Katayama, H.; Paczesny, S.; Prentice, R.; Aragaki, A.; Faca, V.M.; Pitteri, S.J.; Zhang, Q.; Wang, H.; Silva, M.; Kennedy, J.; et al. Application of serum proteomics to the Women’s Health Initiative conjugated equine estrogens trial reveals a multitude of effects relevant to clinical findings. Genome Med. 2009, 1, 47. [Google Scholar] [CrossRef]
  52. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  53. Nueda, M.J.; Tarazona, S.; Conesa, A. Next, maSigPro: Updating maSigPro bioconductor package for RNA-seq time series. Bioinformatics 2014, 30, 2598–2602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016, 44, W90–W97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Sample Availability: All relevant data are available within the Article and the Supplementary files, or available from the authors upon request.
Figure 1. Analyses of IgG, IgM, and IgA immuno-complexes upon chemotherapy (CT) in sera of PDA patients. (A) scheme of the experimental design of the immune-complexome proteomic analysis using pulled sera from PDA patients collected before and after CT (BCT and ACT). Samples were divided into pools based on the months of disease-free survival; (B) Bar plot representing the circulating level of TTR (IgG differential associated antigen (DAA)) BCT and ACT; (C) bar plots representing the most significant trends of IgG Circulating Immune Complexes (CIC) considering BCT and ACT results separately. Both IgG-IGLV7-46 (left panel) and IgG-S100A9 (right panel) increased from good to bad prognosis in the five groups; (D) heat map showing IgM DAA in the five groups of patients; (E) bar plots representing the most significant trends of IgM CIC considering BCT and ACT results separately. Both IgM-APCS (left panel) and IgM-GAPDH (right panel) increased from good to bad prognosis in the five groups; (F) heat map showing IgA DAA in the five groups of patients; (G) bar plots representing the most significant trends of IgA CIC considering BCT and ACT results separately. IgA-SERPINC1 (left panel) decreased from good to bad prognosis and IgA-PTGDS (right panel) increased from good to bad prognosis in the five groups. All CIC are represented as number of Peptide-Spectrum Matches (PSM).
Figure 1. Analyses of IgG, IgM, and IgA immuno-complexes upon chemotherapy (CT) in sera of PDA patients. (A) scheme of the experimental design of the immune-complexome proteomic analysis using pulled sera from PDA patients collected before and after CT (BCT and ACT). Samples were divided into pools based on the months of disease-free survival; (B) Bar plot representing the circulating level of TTR (IgG differential associated antigen (DAA)) BCT and ACT; (C) bar plots representing the most significant trends of IgG Circulating Immune Complexes (CIC) considering BCT and ACT results separately. Both IgG-IGLV7-46 (left panel) and IgG-S100A9 (right panel) increased from good to bad prognosis in the five groups; (D) heat map showing IgM DAA in the five groups of patients; (E) bar plots representing the most significant trends of IgM CIC considering BCT and ACT results separately. Both IgM-APCS (left panel) and IgM-GAPDH (right panel) increased from good to bad prognosis in the five groups; (F) heat map showing IgA DAA in the five groups of patients; (G) bar plots representing the most significant trends of IgA CIC considering BCT and ACT results separately. IgA-SERPINC1 (left panel) decreased from good to bad prognosis and IgA-PTGDS (right panel) increased from good to bad prognosis in the five groups. All CIC are represented as number of Peptide-Spectrum Matches (PSM).
Cancers 12 00746 g001
Table 1. List of differentially recognized Circulating Immune Complexes.
Table 1. List of differentially recognized Circulating Immune Complexes.
Gene NameIg SubtypeVariation Following Chemotherapy
TTRIgGdecreased
ADIPOQIgMincreased
CDH5IgMincreased
FCGBPIgMincreased
HPIgMdecreased
KRT16IgMdecreased
NEO1IgMincreased
SERPINA3IgMincreased
SERPIND1IgMincreased
ABCA13IgAincreased
ECM1IgAdecreased
IGHG3IgAdecreased
IGKV3-20IgAdecreased
IGLV1-40IgAincreased
IGLV7-43IgAdecreased
XIRP2IgAincreased
Table 2. List of significant Circulating Immune Complexes trends before chemotherapy.
Table 2. List of significant Circulating Immune Complexes trends before chemotherapy.
Gene NameIg SubtypeTrend (from Group 1 to 5)
IGKV6D-21IgGincreasing
IGLV7-46IgGincreasing
AOC3IgMdecreasing
APCSIgMincreasing
APOBIgMdecreasing
IGLV3-25IgMdecreasing
APOL1IgAincreasing
IGHV3OR16-9IgAincreasing
MMRN2IgAincreasing
SERPINC1IgAdecreasing
SERPIND1IgAincreasing
Table 3. List of significant Circulating Immune Complexes trends after chemotherapy.
Table 3. List of significant Circulating Immune Complexes trends after chemotherapy.
Gene NameIg SubtypeTrend (from Groups 1 to 5)
APOC2IgGdecreasing
C4BIgGdecreasing
S100A9IgGincreasing
GAPDHIgMincreasing
IGLV1-47IgMdecreasing
OLFM2IgMdecreasing
C9IgAincreasing
FREM2IgAincreasing
IGHDIgAincreasing
KCNQ2IgAincreasing
PTGDSIgAincreasing

Share and Cite

MDPI and ACS Style

Mandili, G.; Follia, L.; Ferrero, G.; Katayama, H.; Hong, W.; Momin, A.A.; Capello, M.; Giordano, D.; Spadi, R.; Satolli, M.A.; et al. Immune-Complexome Analysis Identifies Immunoglobulin-Bound Biomarkers That Predict the Response to Chemotherapy of Pancreatic Cancer Patients. Cancers 2020, 12, 746. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12030746

AMA Style

Mandili G, Follia L, Ferrero G, Katayama H, Hong W, Momin AA, Capello M, Giordano D, Spadi R, Satolli MA, et al. Immune-Complexome Analysis Identifies Immunoglobulin-Bound Biomarkers That Predict the Response to Chemotherapy of Pancreatic Cancer Patients. Cancers. 2020; 12(3):746. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12030746

Chicago/Turabian Style

Mandili, Giorgia, Laura Follia, Giulio Ferrero, Hiroyuki Katayama, Wang Hong, Amin A. Momin, Michela Capello, Daniele Giordano, Rosella Spadi, Maria Antonietta Satolli, and et al. 2020. "Immune-Complexome Analysis Identifies Immunoglobulin-Bound Biomarkers That Predict the Response to Chemotherapy of Pancreatic Cancer Patients" Cancers 12, no. 3: 746. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12030746

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

Article Metrics

Back to TopTop