Purpose:

Intraductal papillary mucinous neoplasm (IPMN) is a precursor of pancreatic ductal adenocarcinoma. Low-grade dysplasia has a relatively good prognosis, whereas high-grade dysplasia and IPMN invasive carcinoma require surgical intervention. However, diagnostic distinction is difficult. We aimed to identify biomarkers in peripheral blood for accurate discrimination.

Experimental Design:

Sera were obtained from 302 patients with IPMNs and 88 healthy donors. For protein biomarkers, serum samples were analyzed on microarrays made of 2,977 antibodies. A support vector machine (SVM) algorithm was applied to define classifiers, which were validated on a separate sample set. For microRNA biomarkers, a PCR-based screen was performed for discovery. Biomarker candidates confirmed by quantitative PCR were used to train SVM classifiers, followed by validation in a different sample set. Finally, a combined SVM classifier was established entirely independent of the earlier analyses, again using different samples for training and validation.

Results:

Panels of 26 proteins or seven microRNAs could distinguish high- and low-risk IPMN with an AUC value of 95% and 94%, respectively. Upon combination, a panel of five proteins and three miRNAs yielded an AUC of 97%. These values were much better than those obtained in the same patient cohort by using the guideline criteria for discrimination. In addition, accurate discrimination was achieved between other patient subgroups.

Conclusions:

Protein and microRNA biomarkers in blood allow precise diagnosis and risk stratification of IPMN cases, which should improve patient management and thus the prognosis of IPMN patients.

See related commentary by Löhr and Pantel, p. 1387

Translational Relevance

IPMN usually develops from benign forms (low risk) to a premalignant stage and finally IPMN-associated carcinoma (high-risk). Patients with low-risk IPMN are kept under close surveillance, whereas the others have the tumor surgically resected. Currently, distinction between the two risk groups is difficult to achieve. Late intervention, however, results in poor prognosis. Inversely, if surgery is performed when not required, it poses an unnecessary risk to patient health. The definition of a highly accurate diagnostic panel of five proteins and three miRNAs in the serum of IPMN patients allows non-invasive, blood-based discrimination of the risk groups. In a cohort that embodies the actual representation of IPMN patients, for whom such a distinction is required, the biomarker signature markedly outperformed the discrimination accuracy obtained following the current clinical guidelines. The established process could substantially affect clinical disease management and thereby considerably improve patient prognosis.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive malignancies (1) and bound to become the second most frequent cause of cancer-related death in the Western world (2). There are two types of PDAC precursor lesions, namely pancreatic intraepithelial neoplasias (PanIN) and mucinous cysts. The latter are subdivided further into intraductal papillary mucinous neoplasms (IPMN) and mucinous cystic neoplasms (MCN). In particular IPMNs have a high potential to develop into malignant tumors with a bad prognosis and are therefore of high clinical relevance. They are frequently discovered through incidental cross-sectional imaging (3, 4). Diagnostically, three IPMN grades are described: low-grade and high-grade dysplasia (also called carcinoma in situ) as well as IPMN with an associated carcinoma (5). On the basis of clinical evidence, low-grade IPMN dysplasia are considered benign with a low risk of malignant progression. They have a 5-year disease-specific survival rate of 97%, as compared with 84% and 39% of patients with high-grade and invasive IPMN (6). They are kept under close surveillance to monitor further development of the disease. High-grade dysplasia is a premalignant form and IPMNs with an associated carcinoma have already become malignant. Because of the high risk, these IPMN patients require surgical tumor resection (7, 8). Distinguishing high-risk from low-risk IPMN patients accurately is critical to give patients the appropriate treatment (9). If the intervention is late, a highly malignant tumor could develop which results in poor prognosis. If surgery is performed when not required, it poses an unnecessary risk to the patient's health. There are quite a few indicators for surgery, such as a positive cytology for malignant, high-grade dysplasia, tumor-related jaundice, enhancing mural nodes, main pancreatic duct dilatation, tumor growth-rate, cyst diameter, increased carbohydrate antigen 19–9 (CA19–9) blood concentration, new-onset of diabetes mellitus, and diagnosis of pancreatitis (10, 11). Despite this large number of diagnostic features, however, diagnosis is not sufficiently accurate. Also imaging does not solve the problem; resected IPMN samples frequently reveal other grades of malignancy than predicted prior to resection (12).

Blood-based diagnostics of PDAC has been developing during the last few years. Blood is attractive for biomarker detection because of its stability and accessibility. Also, diagnosis would be minimally invasive, reducing the burden to the patients and cost. In addition, samples could be easily collected repeatedly from the same patient, allowing longitudinal monitoring. A large number of blood biomarkers or biomarker signatures for the diagnosis of PDAC has been published (for a comprehensive review see ref. 13). However, not much progress has been made toward actual translation to clinical routine. Even the utility of serum markers that are already in use in clinical diagnostics, such as CA19–9, is limited by a relatively high degree of false results.

There have been far fewer studies about blood-based detection of IPMN patients or the discrimination of IPMN grades than there are on PDAC and they are mostly lacking validation (14). In addition, performance is insufficient for application. Serum levels of CA19–9 are associated with malignant IPMN, but the diagnostic power is low (11, 15). Reported AUC values vary, ranging from 0.62 to 0.78. Another serum marker—carcinoembryonic antigen (CEA)—exhibits little utility in patient management. Also the blood type was proposed as a predictor of high-grade dysplasia but was of rather limited accuracy (16). At the RNA level, miRNA miR-4539 isolated from extracellular vesicles was reported to discriminate IPMN patients from healthy individuals with an AUC of 0.72; miR-6132 could distinguish benign IPMN and IPMN-derived carcinoma with an AUC of 0.77 (17). Five other miRNAs were found to correlate somehow with high-grade dysplasia or malignant carcinoma (AUC = 0.73; ref. 18) and a signature of 8 long noncoding RNAs yielded an AUC of 0.77 (19). With respect to proteins, MIC-1/GDF15 was considered, but proved ineffective (20). Also THBS2 showed very low accuracy with an AUC of 0.65 (21). Detection of anti-p53 antibody in combination with serum concentrations of CEA and CA19–9 could discriminate high-grade dysplasia and malignant IPMN from low-grade IPMN with a sensitivity of just 38.4% and specificity of 81.6% (22). The relatively strong anti-p53 antibody marker, however, could be detected in only 8.2% of the samples. Thus, for most patients, specificity was actually substantially lower than 81.6%. Recently, variations in the abundance of six serum proteins were described that distinguished low- and high-risk IPMNs with an accuracy of 83% (23). However, the analysis was based on a small sample number and not validated in an independent sample set.

We aimed at substantially better discrimination by performing a much more comprehensive study, analyzing a large number of serum samples and studying variations in two molecule classes, namely proteins and miRNAs. The samples were collected from patients with low-grade or high-grade dysplasia as well as invasive carcinoma. Serum from healthy blood donors was analyzed for reference. Besides producing separate protein and miRNA diagnostic signatures, we generated independently a classifier that combined both molecule classes as several studies suggested that a combination of diagnostic features of different molecule classes could yield better robustness (13). The main objective of our study was the discrimination between low- and high-risk patients, because of this is most relevant clinically. An accurate biomarker signature in serum could substantially simplify risk stratification and thus a decision about performing tumor resection or not.

Serum samples

Serum samples from 302 IPMN patients and 88 healthy donors were obtained from the Pancobank biorepository of the European Pancreas Center at the Surgery Department of University Hospital Heidelberg. Written informed consent had been obtained from all blood donors. Ethical approval was given by the Ethics Committee of Heidelberg University (ethics votes 159/2002 and 708/2019). Work was done in compliance with the provisions of the Declaration of Helsinki. All patients underwent IPMN resection because of the suspicion of malignancy or other indications for surgery. Only sera were included, for which pathologic diagnosis of the resected material confirmed IPMN with a specific grade of dysplasia. For further sample information, see Supplementary Table S1.

Protein profiling

Microarrays containing 2,977 antibodies targeting 2,286 proteins were produced as described in detail earlier (24); a list of the antibodies’ target proteins is given in Supplementary Table S2. The antibodies were spotted onto epoxysilane-coated slides (Schott Nexterion) using contact printing (MicroGrid-2; BioRobotics). Positional marker molecules as well as negative and positive controls were included. The slides were stored in dry and dark conditions at 4°C until use. Microarray analyses were performed as described in detail elsewhere (25). The protein concentration of the serum samples was adjusted to 4 mg/mL. Proteins were labeled with the fluorescent dyes DY649 or DY549 (Dyomics), respectively, at a molar dye/protein ratio of 7.5 at 4°C for 2 hours. Unreacted dye was quenched by the addition of 10% glycine. As a common reference sample for normalization, we utilized the pooled cellular protein lysates of 24 PDAC cell lines; it was analyzed simultaneously to each serum sample labeled with the other respective dye. Incubation was at 40 μg/mL protein at 4°C overnight. After washing and drying, a PowerScanner system (Tecan) was used for image capture at constant laser power and photo-multiplier tube gain. The images were analyzed with the software package GenePix Pro 6.0 (Molecular Devices), generating numerical values of signal intensities.

miRNA profiling

Total RNA was purified from 100 μL serum with the miRNeasy Serum/Plasma Kit (Qiagen). cDNA was synthesized with the miRCURY LNA RT Kit and directly studied with the miRCURY LNA SYBR Green PCR Kit and the miRCURY LNA Human Serum/Plasma Focus miRNA PCR Panel (Qiagen) according to the manufacturer's protocols. It determined the abundance of 179 miRNAs (Supplementary Table S3). The cycle threshold (Ct) values were transformed to miRNA copy numbers by taking advantage of standard curves generated with the spike-in control cel-miR-39–3p.

qPCR

miRNA biomarker candidates were further studied by qPCR assays on more serum samples. The miRCURY LNA RT Kit was used for reverse transcription of RNA purified with the miRNeasy Serum/Plasma Kit. For each reaction, 5 μL of 2xmiRCURY SYBR Green Master mix and 1 μL PCR primer were added to 4 μL of cDNA, which had been diluted 1:20. After an initial heat activation for 2 minutes at 95°C, there were 45 PCR cycles at 95°C for 10 seconds and 56°C for 1 minute on a LightCycler 480 (Roche Diagnostics) according to the protocol of the miRCURY LNA miRNA SYBR Green PCR Kit. Simultaneously, a standard curve based on cel-miR-39–3p was produced to transform qPCR data into miRNA copy numbers.

Statistical analyses

Protein biomarkers

Analyses were conducted with the R software (version 4.0.4; RRID: SCR_001905). Data normalization and preprocessing were performed with the “limma” package (v. 3.50; ref. 26). The function “backgroundCorrect” was used for background correction (27). Variances within and between microarrays were corrected by applying the functions “normalizeWithinArrays,” “normalizeBetweenArrays,” and “removeBatchEffect” (28). For quality control, the function “boxplot” and cluster dendrograms from the “cluster” package (v. 2.1.3) were used (29). The samples were randomly divided into a discovery and validation cohort. In the discovery cohort, empirical Bayes test was used for comparing the data of two groups, such as healthy controls versus IPMN samples, using the “limma” package. To evaluate the diagnostic power of individual proteins, univariate logistic regression combined with ROC analysis (30) was performed; proteins with a P-value <0.05 were considered biomarker candidates. Least absolute shrinkage and selection operator (LASSO) regression based on the R package “glmnet” (v. 4.1–3; ref. 31) and recursive feature elimination (RFE) with 5-fold cross-validation of the “caret classification and regression” tool (v. 6.0–91) within the R software (RRID: SCR_001905) were applied for feature selection to enhance predictive accuracy (32). Identified protein biomarker candidates were used for the construction and training of a support vector machine (SVM) classifier for two-group classifications using the “e1071” tool (v. 1.7–9) within the R software (RRID: SCR_001905). During the training process, the three parameters cost (ranging from 0 to10), gamma (0 to 1), and epsilon (0 to 1) of the SVM classifier were optimized by an exhaustive search approach. The classifier was then applied with fixed parameters to assess its diagnostic performance with the data resulting from the validation cohort. Values for sensitivity, specificity and the AUC were calculated using the R package “pROC” (30).

miRNA biomarkers

For identifying miRNA biomarkers in serum, candidates were identified first by studying abundance differences in a discovery sample cohort detected by the miRCURY LNA PCR System with the R package “limma” (26). Next, qPCR was conducted for these miRNA candidates using serum samples of a training cohort. Differentially abundant miRNAs with P values <0.05 in two-tailed Student t test were used to define miRNA panels by RFE with 5-fold cross-validation. These panels were then applied to construct and train SVM classifiers. The SVM classifiers with fixed parameters were applied to the data generated from the validation sample cohort to assess the diagnostic performance. AUC values were calculated using the R package “pROC” (30).

Combining protein and miRNA biomarkers

Samples, for which both protein and qPCR miRNA profiling data were available, were randomly separated into a training and validation cohort. Empirical Bayes test and Student t test were used to confirm significant changes of protein and miRNA abundance, respectively. After RFE with 5-fold cross-validation, we used the combined miRNA and protein data resulting from the training cohort to construct and optimize a SVM classifier. It was then applied with fixed parameters to the data resulting from the validation samples to confirm its prediction and classification potential and determine the diagnostic value of the combined biomarker panel.

Data availability

The data generated in this study are available within the article and its supplementary data files. For additional data requests, please contact the corresponding author.

Clinical information about IPMN patients

Serum samples representing 302 IPMN patients and 88 healthy donors (Supplementary Table S1) were obtained from the European Pancreas Center at University Hospital Heidelberg. All patients underwent IPMN resection because of the suspicion of malignancy or other indications for surgery and represent a typical cohort. Blood had been collected prior to tumor resection. Only samples of patients were included, for whom pathologic analysis of the tissue confirmed IPMN with a specific grade of dysplasia. Demographic characteristics, preoperative symptoms, radiologic findings, laboratory data, and pathologic factors were collected. Preoperative cross-sectional imaging was reviewed for radiographic IPMN features, including cyst localization and size, presence of a solid component, thickened/enhancing cyst wall, lymphadenopathy, and main pancreatic duct dilatation. If multiple IPMN lesions were detected, the largest cyst size was recorded. IPMNs were classified as main-duct, branch-duct, or mixed-type based on preoperative imaging and postoperative histology according to international consensus guidelines (15). Preoperative serum levels of CA19–9 and bilirubin as well as the presence of symptoms for obstructive jaundice and pancreatitis were documented. Following surgical resection, an expert pancreatic pathologist assessed the specimens. The highest degree of dysplasia identified in each specimen was recorded as the grade of dysplasia. The study was done in accordance with the Standards for Reporting Diagnostic Accuracy Studies (33).

Protein signatures

For the identification of protein biomarkers, microarrays made of 2,977 antibodies were incubated with serum to acquire protein abundance information. The serum samples were divided randomly into a discovery and a validation cohort. The former consisted of 55 samples from healthy donors, 130 from patients of low-risk and 89 of high-risk IPMN; the overall design of the analysis is shown in Fig. 1. The normalized data from the microarray experiments is accessible in Supplementary Table S4. A relatively large number of proteins exhibited significant abundance variations in at least one of the pairwise comparisons of patient groups (Table 1). However, no protein was sufficiently informative on its own. Because of the molecular heterogeneity of pancreatic cancers, it can be expected that a biomarker panel will be required for an appropriately robust assay. Therefore, absolute shrinkage and selection operator (LASSO) regression and recursive feature elimination (RFE) were applied for further selection. The resulting biomarkers were used to train support vector machine (SVM) classifiers. The SVM algorithm determines the location of samples in a high-dimensional space based on the levels of the biomarker candidates in each sample; each axis in the space represents one protein. The SVM classifier draws a hyperplane, which is a decision boundary to best separate samples into two classes. During the training process, a combination of three parameters (cost, gamma, epsilon) was optimized to avoid over- and underfitting. The classifier with fixed parameter settings was subsequently assessed on data from separate validation samples from 33 healthy donors, 46 low-risk and 37 high-risk IPMN patients.

Figure 1.

Performance of serum protein panels to distinguish IPMN patients from healthy individuals or discriminate between different malignancy risks. A, Workflow of the process for the identification of serum protein biomarker classifiers. B, Performance results are presented as ROC curves and corresponding AUC values as determined in the discovery and validation cohorts, respectively. C, For some typical proteins, the difference is shown as a boxplots, indicating median, first and third quartile as well as maximum and minimum scores.

Figure 1.

Performance of serum protein panels to distinguish IPMN patients from healthy individuals or discriminate between different malignancy risks. A, Workflow of the process for the identification of serum protein biomarker classifiers. B, Performance results are presented as ROC curves and corresponding AUC values as determined in the discovery and validation cohorts, respectively. C, For some typical proteins, the difference is shown as a boxplots, indicating median, first and third quartile as well as maximum and minimum scores.

Close modal
Table 1.

List of serum-protein classifiers.

Number of biomarkersDiscovery cohortValidation cohort
Diagnostic objectiveDifferential abundance analysisUnivariate logistic regressionLASSO regressionRFE with 5-fold cross-validationSVM classifier parametersProtein biomarker panelAUCSenSpeAUCSenSpe
High-risk IPMN 841 728 43 Cost: 9.0 CD72, DOK5, ENO1, POU2F2, SUPT6H + NPY4R 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy     Gamma: 0.01        
     Epsilon: 0.2        
Low-risk IPMN 928 798 23 Cost: 0.2 CD72, DOK5, ENO1, POU2F2, SUPT6H 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy     Gamma: 0.1 + NPY4R + MMP26       
     Epsilon: 0.1        
IPMN 870 749 34 Cost: 2.0 CD72, DOK5, ENO1, POU2F2, SUPT6H + DNAJC17 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy     Gamma: 0.2        
     Epsilon: 0.1        
High-risk IPMN 371 362 97 26 Cost: 0.2 ABT1, BRAF, CCT4, EEF1A1, ETFDH, 0.98 0.99 0.90 0.95 0.84 0.96 
vs. low-risk IPMN     Gamma: 0.05 FGR, GSTM4, IL18, IMPDH1, IPO13,       
     Epsilon: 0.3 KLHL23, L1CAM, LRRC20, NCOR1, NETO2, OSR2, PHAX, PIM1, PVALB, RPH3AL, SERGEF, SYVN1, TECPR1, TJP2, TMEM161A, WDR54       
High-risk IPMN 371 362 97 17 Cost: 0.1 ABT1, EEF1A1, ETFDH, FGR, GSTM4, 0.96 0.94 0.89 0.91 0.84 0.85 
vs. low-risk IPMN     Gamma: 0.06 IPO13, KLHL23, L1CAM, LRRC20, NCOR1, NETO2, PHAX, PIM1, RPH3AL, SYVN1, TMEM161A, WDR54       
     Epsilon: 0.08        
Malignant IPMN 778 678 43 Cost: 1.0 CD72, DOK5, SUPT6H + DNAJC17 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy     Gamma: 0.05        
     Epsilon: 0.8        
Number of biomarkersDiscovery cohortValidation cohort
Diagnostic objectiveDifferential abundance analysisUnivariate logistic regressionLASSO regressionRFE with 5-fold cross-validationSVM classifier parametersProtein biomarker panelAUCSenSpeAUCSenSpe
High-risk IPMN 841 728 43 Cost: 9.0 CD72, DOK5, ENO1, POU2F2, SUPT6H + NPY4R 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy     Gamma: 0.01        
     Epsilon: 0.2        
Low-risk IPMN 928 798 23 Cost: 0.2 CD72, DOK5, ENO1, POU2F2, SUPT6H 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy     Gamma: 0.1 + NPY4R + MMP26       
     Epsilon: 0.1        
IPMN 870 749 34 Cost: 2.0 CD72, DOK5, ENO1, POU2F2, SUPT6H + DNAJC17 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy     Gamma: 0.2        
     Epsilon: 0.1        
High-risk IPMN 371 362 97 26 Cost: 0.2 ABT1, BRAF, CCT4, EEF1A1, ETFDH, 0.98 0.99 0.90 0.95 0.84 0.96 
vs. low-risk IPMN     Gamma: 0.05 FGR, GSTM4, IL18, IMPDH1, IPO13,       
     Epsilon: 0.3 KLHL23, L1CAM, LRRC20, NCOR1, NETO2, OSR2, PHAX, PIM1, PVALB, RPH3AL, SERGEF, SYVN1, TECPR1, TJP2, TMEM161A, WDR54       
High-risk IPMN 371 362 97 17 Cost: 0.1 ABT1, EEF1A1, ETFDH, FGR, GSTM4, 0.96 0.94 0.89 0.91 0.84 0.85 
vs. low-risk IPMN     Gamma: 0.06 IPO13, KLHL23, L1CAM, LRRC20, NCOR1, NETO2, PHAX, PIM1, RPH3AL, SYVN1, TMEM161A, WDR54       
     Epsilon: 0.08        
Malignant IPMN 778 678 43 Cost: 1.0 CD72, DOK5, SUPT6H + DNAJC17 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy     Gamma: 0.05        
     Epsilon: 0.8        

Abbreviations: Sen, sensitivity; Spe, specificity; vs., versus.

The resulting biomarker panels exhibited very high accuracy (Fig. 1). For the discrimination of IPMN and healthy individuals, a classifier was constructed from six proteins. It perfectly distinguished sera from IPMN patients and healthy individuals within the discovery cohort. Absolute accuracy was also achieved upon validation (Table 1). Five of the six biomarker molecules—CD72, DOK5, ENO1, POU2F2, and SUPT6H—are also part of the 7-protein signature that discriminated low-risk cases from healthy individuals and the 6-protein panel that separated high-risk patients from healthy individuals; in all these comparisons, the AUC value was 1.00 (Table 1).

Discrimination between patients with high- and low-risk IPMN was more difficult to achieve. A large group of 26 proteins was required (Table 1; Fig. 1), yielding the best performance with an AUC value of 0.98 in the discovery cohort. This value was slightly reduced to 0.95 upon validation. While discrimination power is very high with a specificity of 96%, the panel lacks some sensitivity at 84%. This is nevertheless substantially better than the performance of CA19–9, for example. Problematic in terms of translation to application is the large number of proteins. We therefore tried to define a smaller panel, setting an AUC threshold at 0.90. The smallest possible protein panel was made of 17 proteins achieving an AUC of 0.91 (Table 1; Fig. 1).

We also looked at distinguishing either high-risk IPMN from a combination of low-risk IPMN cases and healthy individuals or low-risk patients from high-risk IPMN and healthy donors. Again, mostly relatively large marker panels were required (Supplementary Table S5; Supplementary Fig. S1). The overall smallest number of biomarkers consists of four proteins, which discriminated samples from patients with malignant IPMN and those from healthy people with an accuracy of 100% (Table 1). These suggest that the number of required biomarkers gets smaller the wider the pathologic distance is between two patient groups.

miRNA signatures

We also investigated the serum samples’ miRNA content for its diagnostic utility. The overall workflow of defining miRNA biomarkers is displayed in Fig. 2. First a commercially available PCR-array assay that investigates the levels of 179 miRNAs was performed on a discovery cohort of 10 samples from healthy donors, 20 low-risk and 20 high-risk IPMN cases. We chose to focus on this subset of 179 miRNAs, because they represent molecules that are relatively abundant in serum and therefore likely to enable robust measurements. Identified biomarker candidates were then studied by qPCR in 24 healthy, 43 low-risk and 51 high-risk IPMN samples. The results were used to train SVM classifiers. For validation of their diagnostic performance, 16 normal, 19 low-risk, and 22 high-risk IPMN samples were analyzed at fixed classifier parameters. The normalized miRNA copy numbers resulting from all three analysis phases are shown in Supplementary Table S6. As before, diagnostic performance is presented as ROC curves based on the qPCR results generated on the training and validation sample cohorts (Fig. 2), respectively.

Figure 2.

Performance of serum miRNA panels to distinguish IPMN patients from healthy individuals or discriminate between different malignancy risks. A, Workflow of the process for the identification of serum miRNA biomarker classifiers. B, Performance results are presented as ROC curves and corresponding AUC values as determined in the training and validation cohorts, respectively. C, For some typical miRNA molecules, the variations in their abundance level across IPMN grades are shown as boxplots, indicating median, first and third quartile as well as maximum and minimum scores. A, adenoma; B, borderline; C, carcinoma in situ; I, invasive carcinoma.

Figure 2.

Performance of serum miRNA panels to distinguish IPMN patients from healthy individuals or discriminate between different malignancy risks. A, Workflow of the process for the identification of serum miRNA biomarker classifiers. B, Performance results are presented as ROC curves and corresponding AUC values as determined in the training and validation cohorts, respectively. C, For some typical miRNA molecules, the variations in their abundance level across IPMN grades are shown as boxplots, indicating median, first and third quartile as well as maximum and minimum scores. A, adenoma; B, borderline; C, carcinoma in situ; I, invasive carcinoma.

Close modal

In terms of the size of the biomarker panels, miRNA-based discrimination performed better than the proteins. Although the achieved accuracy was similar, the actual number of molecules was much smaller. As seen for the proteins, there was a common set of biomarkers—three in case of miRNA: miR-122–5p, miR-125b-5p, miR-146a-5p—that formed the core of the miRNA panels (Table 2). On their own, they allowed an accurate discrimination of high-risk IPMN from healthy individuals. By adding two more miRNAs—miR-365a-3p and miR-375—also low-risk IPMN could be split from healthy. Taking into account a sixth and seventh molecules, finally, permitted to distinguish between high- and low-risk IPMN with very high accuracy; the AUC upon validation was 0.94 (Table 2). Further comparisons of patient groups also yielded high AUC values (Supplementary Table S7; Supplementary Fig. S2). For none of them there was a panel larger than seven biomarkers.

Table 2.

List of serum-miRNA classifiers.

Training cohortValidation cohort
Diagnostic objectBiomarker panel sizeSVM classifier parametersmiRNA biomarker panelAUCSenSpeAUCSenSpe
High-risk IPMN Cost: 0.2 miR-122–5p, miR-125b-5p, miR-146a-5p 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy  Gamma: 0.2        
  Epsilon: 0.6        
Low-risk IPMN Cost: 9.0 miR-122–5p, miR-125b-5p, miR-146a-5p 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy  Gamma: 0.9 + miR-365a-3p, miR-375       
  Epsilon: 0.9        
IPMN Cost: 4.0 miR-122–5p, miR-125b-5p, miR-146a-5p 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy  Gamma: 0.9 + miR-365a-3p, miR-375       
  Epsilon: 0.9 + miR-130a-3p       
High-risk IPMN Cost: 0.7 miR-122–5p, miR-125b-5p, miR-146a-5p 0.95 0.94 0.84 0.94 0.86 0.95 
vs. low-risk IPMN  Gamma: 0.6 + miR-365a-3p, miR-375       
  Epsilon: 0.5 + miR-155–5p, miR-205–5p       
Training cohortValidation cohort
Diagnostic objectBiomarker panel sizeSVM classifier parametersmiRNA biomarker panelAUCSenSpeAUCSenSpe
High-risk IPMN Cost: 0.2 miR-122–5p, miR-125b-5p, miR-146a-5p 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy  Gamma: 0.2        
  Epsilon: 0.6        
Low-risk IPMN Cost: 9.0 miR-122–5p, miR-125b-5p, miR-146a-5p 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy  Gamma: 0.9 + miR-365a-3p, miR-375       
  Epsilon: 0.9        
IPMN Cost: 4.0 miR-122–5p, miR-125b-5p, miR-146a-5p 1.00 1.00 1.00 1.00 1.00 1.00 
vs. healthy  Gamma: 0.9 + miR-365a-3p, miR-375       
  Epsilon: 0.9 + miR-130a-3p       
High-risk IPMN Cost: 0.7 miR-122–5p, miR-125b-5p, miR-146a-5p 0.95 0.94 0.84 0.94 0.86 0.95 
vs. low-risk IPMN  Gamma: 0.6 + miR-365a-3p, miR-375       
  Epsilon: 0.5 + miR-155–5p, miR-205–5p       

Abbreviations: Sen, sensitivity; Spe, specificity; vs., versus.

Combination panel of serum proteins and miRNAs

Although IPMN cases could be distinguished from healthy individuals with an AUC of 1.00, no absolutely correct discrimination was achieved between high- and low-risk patients. Even accumulating all differentially abundant proteins or miRNAs did not result in an AUC value better than the values shown in Tables 1 and 2. Therefore, we looked into improving performance by generating a classifier that integrates differentially abundant proteins and miRNAs. Marker selection started from scratch. Serum samples were randomly separated in a training cohort of 118 and a validation cohort of 57 samples. Identification of differentially abundant molecules and subsequent RFE with 5-fold cross-validation on the data originating from the training cohort formed the basis for training and optimizing a classifier, which was then applied to the validation samples. A biomarker panel of five proteins (EEF1A1, RPH3AL, NCOR1, L1CAM, TMEM161A) and three miRNAs (miR-146a-5p, miR-155–5p, miR-375) was found that could distinguish high-risk from low-risk IPMN cases with AUC values of 0.99 in the training cohort and 0.97 upon validation (Fig. 3). For the results of other separating other patient groups, see Supplementary Table S8.

Figure 3.

Diagnostic performance of a combined panel of protein and miRNA biomarkers to discriminate high-risk from low-risk IPMN. The results are presented as ROC curves and corresponding AUC values as determined in the training and validation cohorts, respectively.

Figure 3.

Diagnostic performance of a combined panel of protein and miRNA biomarkers to discriminate high-risk from low-risk IPMN. The results are presented as ROC curves and corresponding AUC values as determined in the training and validation cohorts, respectively.

Close modal

Comparison with current clinical discrimination criteria

For the IPMN patients, information was available concerning the aspects of preoperative CA19–9 levels, diagnosis of pancreatitis, cyst diameter, thickened/enhancing cyst wall, detection of a solid component, main pancreatic duct dilatation, lymphadenopathy and obstructive jaundice (Supplementary Table S9), which are diagnostically relevant parameters according to international IPMN guidelines (10, 15, 34). We used the same approach for the training of a SVM classifier that had been applied to the protein and miRNA data. The training cohort consisted of 145 sera (92 low-risk, 53 high-risk), whereas the validation was done in 65 samples (40 low-risk, 25 high-risk). Upon validation, an AUC value of 0.81 was achieved (Fig. 4). Because this was higher than expected, we checked further to avoid potential overfitting. Besides re-running the analysis independently several times with different training and validation sets, we also reversed the analysis, using the 65 validation sera for an entirely independent creation of an SVM classifier and the 145 samples originally used for training now utilized for validation. In all analyses, results were obtained that are very similar to the initial finding; in the reverse analysis, validation produced an AUC value of 0.82 (Fig. 4). Furthermore, we imputed missing values by running the “mice” R package. Training and validation cohorts were randomly selected four times. The AUC values in the respective validation cohorts ranged from 0.80 to 0.82 (Supplementary Fig. S3). The results document that the blood-based analysis is substantially more accurate in comparison to the criteria recommended in the guidelines, with an AUC of 0.97 compared with 0.82.

Figure 4.

Diagnostic performance of clinical parameters to discriminate high-risk from low-risk IPMN. The 8 patient characteristics of preoperative CA19–9 levels, diagnosis of pancreatitis, cyst diameter, thickened/enhancing cyst wall, detection of a solid component, main pancreatic duct dilatation, lymphadenopathy, and obstructive jaundice were used for the training of an SVM classifier. The training and validation cohorts consisted of 145 and 65 samples, respectively. A, The discriminative performance of the classifier is shown as ROC curves and corresponding AUC values. B, As a means to check for overfitting, also the reverse analysis was performed. The 65 samples were used for independent training, whereas validation was done on the 145 samples. Again, the resulting ROC curves and corresponding AUC values are shown.

Figure 4.

Diagnostic performance of clinical parameters to discriminate high-risk from low-risk IPMN. The 8 patient characteristics of preoperative CA19–9 levels, diagnosis of pancreatitis, cyst diameter, thickened/enhancing cyst wall, detection of a solid component, main pancreatic duct dilatation, lymphadenopathy, and obstructive jaundice were used for the training of an SVM classifier. The training and validation cohorts consisted of 145 and 65 samples, respectively. A, The discriminative performance of the classifier is shown as ROC curves and corresponding AUC values. B, As a means to check for overfitting, also the reverse analysis was performed. The 65 samples were used for independent training, whereas validation was done on the 145 samples. Again, the resulting ROC curves and corresponding AUC values are shown.

Close modal

In this study, protein and miRNA biomarker panels were identified that discriminate at high accuracy different IPMN progression stages through SVM classifiers that were made of serum proteins, circulating miRNAs or a mixture of both. This includes particularly the discrimination of healthy individuals from IPMN patients and high-risk from low-risk IPMN cases. The former is an important finding toward detection of IPMN; it is a matter of further research to determine, how early during tumor development this might work. The latter was the key objective of our study and allows a risk-stratification once IPMN has been diagnosed and should help improving clinical decision-making. At low risk of developing a malignant tumor, an operation could be avoided, keeping patients under close surveillance instead. At high risk, tumor resection would be performed, improving the patients’ prognosis substantially. The detection of high risk may even be possible earlier during the tumor pathogenesis process, because the test could be easily repeated frequently during the course of the disease. As a third result of clinical potential, we found biomarkers that effectively distinguished between healthy individuals and patients with malignant IPMN. They may become useful for the detection of early-stage asymptomatic PDAC.

The analyses were performed on a sample cohort collected at a single medical institution. An advantage of this is that there should be minimal effects of differences in sample handling. Consequently, the detected variations are likely to be the result of the biological differences of IPMN grades rather than biases caused by serum processing. For an application in a wider clinical setting, however, the performance of the defined biomarker panels has to be re-evaluated in a multicenter and prospective study so as to compensate for such potential handling effects, although common standard procedures were applied. Furthermore, the parameters of the SVM classifier probably require refinement for it to become commonly applicable at many different sites for a diagnosis of individual patients. The more pathologically confirmed samples the classifier would be based on, the more its performance will be fine-tuned and optimized, thereby continuously improving accuracy and robustness, although eventually only by subtle increments. Because IPMN cases are relatively rare and the time span will be long, it will be a challenge to perform such a study. Concerning eventual clinical testing, qPCR-based assays could be used for miRNA biomarkers, whereas multiplex ELISAs could be applied for protein biomarkers.

Already now, however, the number of analyzed samples was large enough to warrant a solid evaluation. We are confident that the identified marker panels are robust. This is supported by the repeated and thus reproducible identification of particular marker molecules in different, independent analyses. In all comparisons, the AUC value achieved with the validation cohort was identical or slightly less than the one calculated for the training cohort. In none of them, the difference was significant, however. This suggests that there was no major overfitting or underfitting generated during classifier calculation and optimization. This is further supported by the fact that repetition with changed training and validation sets or reversal of the analysis did not make a difference. The quality of the identified biomarkers is indicated by the fact that CA19–9 did not contribute to the creation of any of the protein classifiers, although it was also assayed on the antibody microarray. Although CA19–9 exhibited differential abundance, other factors were found to be superior in diagnostic performance.

For our initial screening for miRNA biomarker candidates, we used the miRCURY LNA miRNA PCR system. The studied miRNAs are present in serum at high concentrations. This was done on purpose to focus on molecules, which could act as robust biomarkers. These miRNAs are likely to be detectable and also relatively small changes in abundance may be big enough in absolute value to be significant. This focus on abundant miRNAs may be at the sacrifice of excluding informative miRNA molecules, which are present in serum at only low concentrations but may nevertheless be specific. However, they could be particularly prone to lacking sensitivity, thereby negatively affecting overall performance.

Thirty-seven of the overall 67 protein biomarkers identified in the pairwise comparisons (Table 1; Supplementary Table S1) exhibited correlation with IPMN progression. Of these molecules, only VEGFA had been reported in connection with cystic pancreatic tumors. It was found at elevated concentrations in the cyst fluid of serous cystic neoplasms and could distinguish it from other cystic lesions including IPMN (35). Eighteen proteins—ANXA1, AQP1, CAPN2, CD59, EEF1A1, ELMO2, ENO1, GSTM4, HMMR, IGF2BP3, L1CAM, NCOR1, NETO2, PRTN3, TJP2, TSLP, TTK, and VEGFA—have been shown to be involved in promoting PDAC proliferation or invasion (36–53). Their correlation with IPMN progression suggests that they may be factors that are critical for tumor functioning. However, meaningful conclusions on functional aspects are difficult to be drawn, because most functional information about proteins is based on cellular studies. Variations in serum may have different causes or consequences and represent the status of not just the tumor but the entire organism. Therefore, potential explanations are rather speculative.

Recently, extracellular vesicle-associated MUC5AC was found as a possible marker for invasive IPMN (54), lacking any assessment of non-invasive high-grade samples, however. The antibody microarray used in our study did not include MUC5AC, unfortunately, although assaying 2,286 proteins. Therefore, we were not able to analyze its performance in our study. Besides analyzing blood, there has been interest in developing protein biomarkers from pancreatic cyst fluid for IPMN classification. Protein IL1B was found to be upregulated in both serum samples and pancreatic cyst fluid from high-risk IPMN compared with low-risk IPMN (55). However, IL1B was not included in our SVM classifier for IPMN risk stratification because of its comparably small discriminative power, although more abundant in serum samples of high- than low-risk IPMN patients.

Overall, 10 miRNA molecules were part of the biomarker panels that discriminated best between the different IPMN stages (Table 2; Supplementary Table S2). Five of them have already been studied in IPMN. In our analysis, miR-125b-5p showed negative correlation with IPMN progression, consistent with changes observed in pancreatic cystic fluid (56). Decreases in serum levels of miR-130a-3p from benign to malignant IPMN match findings of such a reduction in IPMN tissues (57). Variations in miR-375 levels were observed in tissues of serous cyst-adenoma, mucinous cystic lesions, and PDAC (58). For miR-146a-5p and miR-155–5p, there was less in IPMN sera compared with samples from healthy individuals; the opposite has been found in tissues (59–61). Such asymmetric distribution has been observed in cancer before (e.g., 62) and suggests an actively regulated secretion of miRNA molecules from cells. Four of the other five molecules—miR-122–5p, miR-205–5p, miR-326, miR-365a-3p—are known to be associated with PDAC (63–66). No such link has been reported for miR-154–5p, but it has been found to act as a tumor suppressor in several tumor entities (67). As for proteins, a functional interpretation of variations in miRNA levels in serum is rather speculative. Also, the entire picture of all changes should be taken into account to such an end rather than focusing on individual molecules.

In conclusion, a small set of serum-based protein and miRNA biomarkers was defined that permits discrimination between patients with high- and low-risk IPMN. This minimally invasive and accurate process could contribute substantially to patient management and improve the prognosis of patients with cystic pancreatic tumors.

S. Roth reports grants from German Cancer Aid during the conduct of the study. No disclosures were reported by the other authors.

C. Zhang: Conceptualization, software, formal analysis, validation, investigation, visualization, methodology, writing–original draft, writing–review and editing. F.N. Al-Shaheri: Formal analysis, validation, investigation, methodology, writing–review and editing. M.S.S. Alhamdani: Software, investigation, methodology. A.S. Bauer: Data curation, software, formal analysis. J.D. Hoheisel: Conceptualization, formal analysis, supervision, investigation, methodology, project administration, writing–review and editing. M. Schenk: Investigation. U. Hinz: Data curation, investigation. P. Goedecke: Data curation, investigation. K. Al-Halabi: Data curation, formal analysis, investigation. M.W. Büchler: Resources, supervision, project administration, writing–review and editing. N.A. Giese: Resources, data curation, formal analysis. T. Hackert: Conceptualization, resources, formal analysis, supervision, validation, project administration, writing–review and editing. S. Roth: Conceptualization, resources, data curation, formal analysis, supervision, validation, project administration, writing–review and editing.

Chaoyang Zhang was supported by a scholarship of the China Scholarship Council.

The publication costs of this article were defrayed in part by the payment of publication fees. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).

1.
Siegel
RL
,
Miller
KD
,
Fuchs
HE
,
Jemal
A
.
Cancer statistics, 2022
.
CA Cancer J Clin
2022
;
72
:
7
33
.
2.
Rahib
L
,
Smith
BD
,
Aizenberg
R
,
Rosenzweig
AB
,
Fleshman
JM
,
Matrisian
LM
.
Projecting cancer incidence and deaths to 2030: the unexpected burden of thyroid, liver, and pancreas cancers in the United States
.
Cancer Res
2014
;
74
:
2913
21
.
3.
Singhi
AD
,
Koay
EJ
,
Chari
ST
,
Maitra
A
.
Early detection of pancreatic cancer: opportunities and challenges
.
Gastroenterology
2019
;
156
:
2024
40
.
4.
Assarzadegan
N
,
Babaniamansour
S
,
Shi
J
.
Updates in the diagnosis of intraductal neoplasms of the pancreas
.
Front. Physiol
2022
;
13
:
856803
.
5.
Distler
M
,
Aust
D
,
Weitz
J
,
Pilarsky
C
,
Grützmann
R
.
Precursor lesions for sporadic pancreatic cancer: PanIN, IPMN, and MCN
.
BioMed Res Int
2014
;
2014
:
474905
.
6.
Rezaee
N
,
Barbon
C
,
Zaki
A
,
He
J
,
Salman
B
,
Hruban
RH
, et al
.
Intraductal papillary mucinous neoplasm (IPMN) with high-grade dysplasia is a risk factor for the subsequent development of pancreatic ductal adenocarcinoma
.
HPB
2016
;
18
:
236
46
.
7.
Carmicheal
J
,
Patel
A
,
Dalal
V
,
Atri
P
,
Dhaliwal
AS
,
Wittel
UA
, et al
.
Elevating pancreatic cystic lesion stratification: Current and future pancreatic cancer biomarker(s)
.
Biochim Biophys Acta Rev Cancer
2020
;
1873
:
188318
.
8.
Basturk
O
,
Hong
SM
,
Wood
LD
,
Adsay
NV
,
Albores-Saavedra
J
,
Biankin
AV
, et al
.
A revised classification system and recommendations from the Baltimore Consensus Meeting for neoplastic precursor lesions in the pancreas
.
Am J Surg Pathol
2015
;
39
:
1730
41
.
9.
Kaiser J Scheifele
C
,
Hinz
U
,
Leonhardt
CS
,
Hank
T
,
Koenig
AK
, et al
.
IPMN-associated pancreatic cancer: survival, prognostic staging and impact of adjuvant chemotherapy
.
Eur J of Surg Oncol
2022
;
48
:
1309
20
.
10.
European Study Group on Cystic Tumours of the Pancreas
.
European evidence-based guidelines on pancreatic cystic neoplasms
.
Gut
2018
;
67
:
789
804
.
11.
Tanaka
M
.
Clinical management and surgical decision-making of IPMN of the pancreas
.
Methods Mol Biol
2019
;
1882
:
9
22
.
12.
Khoury
RE
,
Kabir
C
,
Maker
VK
,
Banulescu
M
,
Wasserman
M
,
Maker
AV
.
What is the incidence of malignancy in resected IPMN? An analysis of over 100 US institutions in a single year
.
Ann Surg Oncol
2018
;
25
:
1746
51
.
13.
Al-Shaheri
FN
,
Alhamdani
MSS
,
Bauer
AS
,
Giese
N
,
Büchler
MW
,
Hackert
T
, et al
.
Blood biomarkers for differential diagnosis and early detection of pancreatic cancer
.
Cancer Treat Rev
2021
;
96
:
102193
.
14.
Midlöv
F
,
Ansari
D
,
Aronsson
L
,
Andersson
R
.
Blood-based markers for predicting malignancy in intraductal papillary mucinous neoplasms of the pancreas: what's beyond the horizon?
Scand J Gastroenterol
2020
;
55
:
515
6
.
15.
Tanaka
M
,
Fernández-Del Castillo
C
,
Kamisawa
T
,
Jang
JY
,
Levy
P
,
Ohtsuka
T
, et al
.
Revisions of international consensus Fukuoka guidelines for the management of IPMN of the pancreas
.
Pancreatology
2017
;
17
:
738
53
.
16.
Poruk
KE
,
Griffin
J
,
Makary
MA
,
He
J
,
Cameron
JL
,
Weiss
MJ
, et al
.
Blood type as a predictor of high-grade dysplasia and associated malignancy in patients with intraductal papillary mucinous neoplasms
.
J Gastrointest Surg
2019
;
23
:
477
83
.
17.
Sato
Y
,
Suzuki
R
,
Takagi
T
,
Sugimoto
M
,
Ohira
H
.
Circulating extracellular vesicle-encapsulated microRNA as screening biomarkers for intraductal papillary mucinous neoplasm
.
Oncol Lett
2020
;
20
:
315
.
18.
Permuth-Wey
J
,
Chen
DT
,
Fulp
WJ
,
Yoder
SJ
,
Zhang
Y
,
Georgeades
C
, et al
.
Plasma microRNAs as novel biomarkers for patients with intraductal papillary mucinous neoplasms of the pancreas
.
Cancer Prev Res
2015
;
8
:
826
34
.
19.
Permuth
JB
,
Chen
DT
,
Yoder
SJ
,
Li
J
,
Smith
AT
,
Choi
JW
, et al
.
Linc-ing circulating long non-coding RNAs to the diagnosis and malignant prediction of intraductal papillary mucinous neoplasms of the pancreas
.
Sci Rep
2017
;
7
:
10484
.
20.
O'Neill
RS
,
Emmanuel
S
,
Williams
D
,
Stoita
A
.
Macrophage inhibitory cytokine-1/growth differentiation factor-15 in premalignant and neoplastic tumours in a high-risk pancreatic cancer cohort
.
World J Gastroenterol
2020
;
26
:
1660
73
.
21.
Simpson
RE
,
Yip-Schneider
MT
,
Wu
H
,
Fan
H
,
Liu
Z
,
Korc
M
, et al
.
Circulating Thrombospondin-2 enhances prediction of malignant intraductal papillary mucinous neoplasm
.
Am J Surg
2019
;
217
:
425
8
.
22.
Hata
T
,
Mizuma
M
,
Motoi
F
,
Iseki
M
,
Omori
Y
,
Hayashi
H
, et al
.
Serum anti-p53 antibody can serve as a predictive marker for histological grade of intraductal papillary mucinous neoplasms of the pancreas
.
Pancreas
2020
;
49
:
768
73
.
23.
Roth
S
,
Bose
P
,
Alhamdani
MSS
,
Mustafa
SA
,
Tjaden
C
,
Zamzow
K
, et al
.
Noninvasive discrimination of low and high-risk pancreatic intraductal papillary mucinous neoplasms
.
Ann Surg
2021
;
273
:
e273
5
.
24.
Marzoq
AJ
,
Mustafa
SA
,
Heidrich
L
,
Hoheisel
JD
,
Alhamdani
ASS
.
Impact of the secretome of activated pancreatic stellate cells on growth and differentiation of pancreatic tumor cells
.
Sci Rep
2019
;
9
:
5303
.
25.
Mustafa
S
,
Pan
L
,
Marzoq
A
,
Fawaz
M
,
Sander
L
,
Rückert
F
, et al
.
Comparison of the tumor cell secretome and patient sera for an accurate serum-based diagnosis of pancreatic ductal adenocarcinoma
.
Oncotarget
2017
;
8
:
11963
76
.
26.
Ritchie
ME
,
Phipson
B
,
Wu
D
,
Hu
Y
,
Law
CW
,
Shi
W
, et al
.
limma powers differential expression analyses for RNA-sequencing and microarray studies
.
Nucleic Acids Res
2015
;
43
:
e47
.
27.
Ritchie
ME
,
Silver
J
,
Oshlack
A
,
Holmes
M
,
Diyagama
D
,
Holloway
A
, et al
.
A comparison of background correction methods for two-colour microarrays
.
Bioinformatics
2007
;
23
:
2700
7
.
28.
Smyth
GK
,
Speed
T
.
Normalization of cDNA microarray data
.
Methods
2003
;
31
:
265
73
.
29.
Maechler
M
,
Rousseeuw
P
,
Struyf
A
,
Hubert
M
,
Hornik
K
.
Cluster: cluster analysis basics and extensions
.
R packages
2012
;
1
.
30.
Robin
X
,
Turck
N
,
Hainard
A
,
Tiberti
N
,
Lisacek
F
,
Sanchez
JC
, et al
.
pROC: an open-source package for R and S+ to analyze and compare ROC curves
.
BMC Bioinf
2011
;
12
:
77
.
31.
Friedman
J
,
Hastie
T
,
Tibshirani
R
.
Regularization paths for generalized linear models via coordinate descent
.
J Stat Softw
2010
;
33
:
1
22
.
32.
Tibshirani
R
.
Regression shrinkage and selection via the lasso: a retrospective
.
J. R. Statist. Soc.B
2011
;
73
:
273
82
.
33.
Cohen
JF
,
Korevaar
DA
,
Altman
DG
,
Bruns
DE
,
Gatsonis
CA
,
Hooft
L
, et al
.
STARD 2015 guidelines for reporting diagnostic accuracy studies: explanation and elaboration
.
BMJ Open
2016
;
6
:
e012799
.
34.
Vege
SS
,
Ziring
B
,
Jain
R
,
Moayyedi
P
.
Clinical Guidelines Committee; American Gastroenterology Association. American gastroenterological association institute guideline on the diagnosis and management of asymptomatic neoplastic pancreatic cysts
.
Gastroenterology
2015
;
148
:
819
22
.
35.
Carr
RA
,
Yip-Schneider
MT
,
Dolejs
S
,
Hancock
BA
,
Wu
H
,
Radovich
M
, et al
.
Pancreatic cyst fluid vascular endothelial growth factor A and carcinoembryonic antigen: a highly accurate test for the diagnosis of serous cystic neoplasm
.
J Am Coll Surg
2017
;
S1072–7515(17)30452–0
.
36.
Novizio
N
,
Belvedere
R
,
Pessolano
E
,
Morello
S
,
Tosco
A
,
Campiglia
P
, et al
.
ANXA1 contained in EVs regulates macrophage polarization in tumor microenvironment and promotes pancreatic cancer progression and metastasis
.
Int J Mol Sci
2021
;
22
:
11018
.
37.
Zou
W
,
Yang
Z
,
Li
D
,
Liu
Z
,
Zou
Q
,
Yuan
Y
.
AQP1 and AQP3 expression are associated with severe symptoms and poor-prognosis of the pancreatic ductal adenocarcinoma
.
Appl Immunohistochem Mol Morphol
2019
;
27
:
40
47
.
38.
Szymczak-Pajor
I
,
Fleszar
K
,
Kasznicki
J
,
Gralewska
P
,
Śliwińska
A
.
A potential role of calpains in sulfonylureas (SUs)-mediated death of human pancreatic cancer cells (1.2B4)
.
Toxicol In Vitro
2021
;
73
:
105128
.
39.
Zhang
R
,
Liu
Q
,
Peng
J
,
Wang
M
,
Gao
X
,
Liao
Q
, et al
.
Pancreatic cancer-educated macrophages protect cancer cells from complement-dependent cytotoxicity by up-regulation of CD59
.
Cell Death Dis
2019
;
10
:
836
.
40.
Selga
E
,
Oleaga
C
,
Ramírez
S
,
de Almagro
MC
,
Noé
V
,
Ciudad
CJ
.
Networking of differentially expressed genes in human cancer cells resistant to methotrexate
.
Genome Med
2009
;
1
:
83
.
41.
Wang
Y
,
Li
H
,
Li
F
.
ELMO2 association with Gαi2 regulates pancreatic cancer cell chemotaxis and metastasis
.
PeerJ
2020
;
8
:
e8910
.
42.
Principe
M
,
Borgoni
S
,
Cascione
M
,
Chattaragada
MS
,
Ferri-Borgogno
S
,
Capello
M
, et al
.
Alpha-enolase (ENO1) controls alpha v/beta 3 integrin expression and regulates pancreatic cancer adhesion, invasion, and metastasis
.
J Hematol Oncol
2017
;
10
:
16
.
43.
Wei
L
,
Wen
JY
,
Chen
J
,
Ma
XK
,
Wu
DH
,
Chen
ZH
, et al
.
Oncogenic ADAM28 induces gemcitabine resistance and predicts a poor prognosis in pancreatic cancer
.
World J Gastroenterol
2019
;
25
:
5590
603
.
44.
Lin
A
,
Feng
J
,
Chen
X
,
Wang
D
,
Wong
M
,
Zhang
G
, et al
.
High levels of truncated RHAMM cooperate with dysfunctional p53 to accelerate the progression of pancreatic cancer
.
Cancer Lett
2021
;
514
:
79
89
.
45.
Ennajdaoui
H
,
Howard
JM
,
Sterne-Weiler
T
,
Jahanbani
F
,
Coyne
DJ
,
Uren
PJ
, et al
.
IGF2BP3 modulates the interaction of invasion-associated transcripts with RISC
.
Cell Rep
2016
;
15
:
1876
83
.
46.
Na'ara
S
,
Amit
M
,
Gil
Z
.
L1CAM induces perineural invasion of pancreas cancer cells by upregulation of metalloproteinase expression
.
Oncogene
2019
;
38
:
596
608
.
47.
Aaquist
T
,
Dembic
M
,
Thomassen
M
,
de Stricker
K
,
Bertelsen
M
,
Christensen
LG
, et al
.
Synchronous detection of pancreatic adenocarcinoma and paraganglioma in a Whipple resection specimen
.
Pathol Res Pract
2021
;
226
:
153590
.
48.
Li
Y
,
Zhang
Y
,
Liu
J
.
NETO2 promotes pancreatic cancer cell proliferation, invasion and migration via activation of the STAT3 signaling pathway
.
Cancer Manag Res
2019
;
11
:
5147
56
.
49.
Hu
D
,
Ansari
D
,
Zhou
Q
,
Sasor
A
,
Said Hilmersson
K
,
Andersson
R
.
Low P4HA2 and high PRTN3 expression predicts poor survival in patients with pancreatic cancer
.
Scand J Gastroenterol
2019
;
54
:
246
51
.
50.
Zhou
L
,
Tan
X
,
Wang
W
,
Wang
B
,
Dai
X
,
Liu
J
.
Analysis of invasion-metastasis in pancreatic cancer: correlation between the expression and arrangement of tight junction protein-2 and cell dissociation in pancreatic cancer cells
.
Mol Med Rep
2010
;
3
:
149
53
.
51.
Brunetto
E
,
De Monte
L
,
Balzano
G
,
Camisa
B
,
Laino
V
,
Riba
M
, et al
.
The IL-1/IL-1 receptor axis and tumor cell released inflammasome adaptor ASC are key regulators of TSLP secretion by cancer associated fibroblasts in pancreatic cancer
.
J Immunother Cancer
2019
;
7
:
45
.
52.
Kaistha
BP
,
Honstein
T
,
Müller
V
,
Bielak
S
,
Sauer
M
,
Kreider
R
, et al
.
Key role of dual specificity kinase TTK in proliferation and survival of pancreatic cancer cells
.
Br J Cancer
2014
;
111
:
1780
7
.
53.
Quan
L
,
Ohgaki
R
,
Hara
S
,
Okuda
S
,
Wei
L
,
Okanishi
H
, et al
.
Amino acid transporter LAT1 in tumor-associated vascular endothelium promotes angiogenesis by regulating cell proliferation and VEGF-A-dependent mTORC1 activation
.
J Exp Clin Cancer Res
2020
;
39
:
266
.
54.
Yang
KS
,
Ciprani
D
,
O'Shea
A
,
Liss
AS
,
Yang
R
,
Fletcher-Mercaldo
S
, et al
.
Extracellular vesicle analysis allows for identification of invasive IPMN
.
Gastroenterology
2021
;
160
:
1345
58
.
55.
Maker
AV
,
Katabi
N
,
Qin
LX
,
Klimstra
DS
,
Schattner
M
,
Brennan
MF
, et al
.
Cyst fluid interleukin-1beta (IL1beta) levels predict the risk of carcinoma in intraductal papillary mucinous neoplasms of the pancreas
.
Clin Cancer Res
2011
;
17
:
1502
8
.
56.
Henry
JC
,
Bassi
C
,
Giovinazzo
F
,
Bloomston
M
.
MicroRNA from pancreatic duct aspirate differentiates cystic lesions of the pancreas
.
Ann Surg Oncol
2013
;
20
:
S661
6
.
57.
Permuth-Wey
J
,
Chen
YA
,
Fisher
K
,
McCarthy
S
,
Qu
X
,
Lloyd
MC
, et al
.
A genome-wide investigation of microRNA expression identifies biologically-meaningful microRNAs that distinguish between high-risk and low-risk intraductal papillary mucinous neoplasms of the pancreas
.
PLoS One
2015
;
10
:
e0116869
.
58.
Lee
LS
,
Szafranska-Schwarzbach
AE
,
Wylie
D
,
Doyle
LA
,
Bellizzi
AM
,
Kadiyala
V
, et al
.
Investigating microRNA expression profiles in pancreatic cystic neoplasms
.
Clin Transl Gastroenterol
2014
;
5
:
e47
.
59.
Habbe
N
,
Koorstra
JB
,
Mendell
JT
,
Offerhaus
GJ
,
Ryu
JK
,
Feldmann
G
, et al
.
MicroRNA miR-155 is a biomarker of early pancreatic neoplasia
.
Cancer Biol Ther
2009
;
8
:
340
6
.
60.
Vila-Navarro
E
,
Vila-Casadesús
M
,
Moreira
L
,
Duran-Sanchon
S
,
Sinha
R
,
Ginés
À
, et al
MicroRNAs for detection of pancreatic neoplasia: Biomarker discovery by next-generation sequencing and validation in 2 independent cohorts
.
Ann Surg
2017
;
265
:
1226
34
.
61.
Lubezky
N
,
Loewenstein
S
,
Ben-Haim
M
,
Brazowski
E
,
Marmor
S
,
Pasmanik-Chor
M
, et al
.
MicroRNA expression signatures in intraductal papillary mucinous neoplasm of the pancreas
.
Surgery
2013
;
153
:
663
72
.
62.
Ragusa
M
,
Maugeri
M
,
Statello
L
,
Barbagallo
C
,
Passanisi
R
,
Alhamdani
MSS
, et al
.
Variations in the miRNA and protein content between colorectal cancer cells and their exosomes following Cetuximab treatment; biomolecular and translational implications
.
Oncoscience
2014
;
1
:
132
57
.
63.
Schultz
NA
,
Dehlendorff
C
,
Jensen
BV
,
Bjerregaard
JK
,
Nielsen
KR
,
Bojesen
SE
, et al
.
MicroRNA biomarkers in whole blood for detection of pancreatic cancer
.
JAMA
2014
;
311
:
392
404
.
64.
Diaz-Riascos
ZV
,
Ginesta
MM
,
Fabregat
J
,
Serrano
T
,
Busquets
J
,
Buscail
L
, et al
.
Expression and role of microRNAs from the miR-200 family in the tumor formation and metastatic propensity of pancreatic cancer
.
Mol Ther Nucleic Acids
2019
;
17
:
491
503
.
65.
Zhang
ZL
,
Bai
ZH
,
Wang
XB
,
Bai
L
,
Miao
F
,
Pei
HH
.
miR-186 and 326 predict the prognosis of pancreatic ductal adenocarcinoma and affect the proliferation and migration of cancer cells
.
PLoS One
2015
;
10
:
e0118814
.
66.
Yin
L
,
Xiao
X
,
Georgikou
C
,
Yin
Y
,
Liu
L
,
Karakhanova
S
, et al
.
MicroRNA-365a-3p inhibits c-Rel-mediated NF-κB signaling and the progression of pancreatic cancer
.
Cancer Lett
2019
;
452
:
203
12
.
67.
Zhang
J
,
Mao
S
,
Wang
L
,
Zhang
W
,
Zhang
Z
,
Guo
Y
, et al
.
MicroRNA-154 functions as a tumor suppressor in bladder cancer by directly targeting ATG7
.
Oncol Rep
2019
;
41
:
819
28
.
This open access article is distributed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license.