Next Article in Journal
A Single Dose of Novel PSMA-Targeting Radiopharmaceutical Agent [177Lu]Ludotadipep for Patients with Metastatic Castration-Resistant Prostate Cancer: Phase I Clinical Trial
Next Article in Special Issue
Identification and Validation of Potentially Clinically Relevant CpG Regions within the Class 2 Tumor Suppressor Gene SFRP1 in Pancreatic Cancer
Previous Article in Journal
Natural Compounds in Liposomal Nanoformulations of Potential Clinical Application in Glioblastoma
Previous Article in Special Issue
Risk Factors for Pancreatic Cancer in Patients with New-Onset Diabetes: A Systematic Review and Meta-Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CT Radiomics and Whole Genome Sequencing in Patients with Pancreatic Ductal Adenocarcinoma: Predictive Radiogenomics Modeling

1
Joint Department of Medical Imaging, Princess Margaret Hospital, University Health Network, University of Toronto, Toronto, ON M5G 2C1, Canada
2
Department of Biostatistics, Princess Margaret Cancer Centre, University Health Network, University of Toronto, Toronto, ON M5G 2C1, Canada
3
Division of Biostatistics, Dalla Lana School of Public Health, University of Toronto, Toronto, ON M5T 3M7, Canada
4
Princess Margaret Cancer Centre, University Health Network, Toronto, ON M5G 2C1, Canada
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(24), 6224; https://doi.org/10.3390/cancers14246224
Submission received: 26 October 2022 / Revised: 2 December 2022 / Accepted: 12 December 2022 / Published: 16 December 2022
(This article belongs to the Special Issue Advances in Diagnosis and Treatment of Pancreatic Cancer)

Abstract

:

Simple Summary

Linking imaging-derived radiomics features to underlying tumor biology and pathogenesis is a developing field of increasing interest, given the wide availability of imaging data in contrast to costs, expenses and logistical issues of molecular analyses. While invasive tissue sampling remains the gold standard for histologic characterization, the usage of noninvasive imaging techniques for diagnosis and detection of specific tumor characteristics could represent a potential additive or eventually an alternative, especially in patients with advanced, inoperable disease or when inaccessible to biopsy. Biomarkers are continuously expanding and several actionable targets have already been identified in pancreatic ductal adenocarcinoma (PDAC), in order to guide clinical decision making and help to develop novel treatment strategies. Our study indicates acceptable correlation of CT-derived radiomics features and driver gene mutations in PDAC.

Abstract

We investigate whether computed tomography (CT) derived radiomics may correlate with driver gene mutations in patients with pancreatic ductal adenocarcinoma (PDAC). In this retrospective study, 47 patients (mean age 64 ± 11 years; range: 42–86 years) with PDAC, who were treated surgically and who underwent preoperative CT imaging at our institution were included in the study. Image segmentation and feature extraction was performed semi-automatically with a commonly used open-source software platform. Genomic data from whole genome sequencing (WGS) were collected from our institution’s web-based resource. Two statistical models were then built, in order to evaluate the predictive ability of CT-derived radiomics feature for driver gene mutations in PDAC. 30/47 of all tumor samples harbored 2 or more gene mutations. Overall, 81% of tumor samples demonstrated mutations in KRAS, 68% of samples had alterations in TP53, 26% in SMAD4 and 19% in CDKN2A. Extended statistical analysis revealed acceptable predictive ability for KRAS and TP53 (Youden Index 0.56 and 0.67, respectively) and mild to acceptable predictive signal for SMAD4 and CDKN2A (Youden Index 0.5, respectively). Our study establishes acceptable correlation of radiomics features and driver gene mutations in PDAC, indicating an acceptable prognostication of genomic profiles using CT-derived radiomics. A larger and more homogenous cohort may further enhance the predictive ability.

1. Introduction

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive solid tumor and is associated with poor prognosis with a 5-year survival rate of 6–7% [1]. PDAC is considered the fourth most common cause of cancer-associated death in the United States, reflected by estimated 57,600 new cases and resulting in 47,050 deaths in 2020 [2]. Due to low detection rates at initial stages, the majority of patients are diagnosed with locally advanced or metastatic disease and are not candidates for potential curative surgery, resulting in median survival between 6.1 and 11 months [3].
In recent years, advanced genomic sequencing technologies and molecular profiling contributed to further characterization of recurrent genetic alterations in PDAC, yielding important insights in tumor biology and improved understanding of familial predisposition while delivering potentially valuable prognostic tools and targeted therapeutics approaches. Recent genomic studies of PDAC have identified high-frequency alterations in the function of many key oncogenes and tumor suppressor genes, including KRAS, TP53, SMAD4 and CDKN2A [4,5,6]. Moreover, RNA analyses of PDAC resection specimens have demonstrated gene expression signatures with prognostic and biological relevance [7,8].
Ongoing improvements in imaging modalities, techniques and postprocessing tools have intensified the radiological assessment and management of PDAC. Particularly recent developments of machine learning techniques and the huge growth of computational power has driven the field of radiomics [9]. The principles of radiomics include extraction of high-dimensional data from various sources of medical images, followed by an analysis of different classes of radiomic features, aiming to support clinical decision-making to overcome the limitations of solely visual image interpretation [10]. Previous studies have demonstrated the potential to describe underlying tumor biology of various cancer types using radiomics [11,12,13].
Therefore, decrypting pathogenesis and tumor biology in correlation with imaging-derived radiomics would be highly desirable, in order to guide clinical decision making and help to develop novel treatment strategies, including a targeted therapy approach. Thus, the primary aim of our study was to investigate whether CT-derived radiomics features may correlate with driver mutations from whole genome sequencing of PDAC resection specimens.

2. Materials and Methods

Based on a database search in cBioPortal, 55 patients with PDAC who underwent Whipple’s procedure or distal pancreatectomy between June 2008 and November 2015 were identified. Patients were excluded due to the lack of preoperative CT imaging (n = 5) and/or missing whole genome sequencing data (n = 3).
Demographic patient data is provided in Table 1. Our study received institutional review board and local ethics committee approval (CAPCR/UHN REB #: 20-6105). All patients provided a written informed consent prior to the study inclusion for evaluation of genome sequencing.

2.1. Image Acquisition

All CT examinations were performed on a 64-slice single source CT scanner (SOMATOM Definition, Siemens Healthineers, Erlangen, Germany). Acquisition parameters were as follows: tube voltage, 120 kVp, maximum allowable tube current set at 200 mAs using automated exposure control (CAREDose 4D), gantry rotation time, 500ms, pitch 1, beam collimation 32  ×  0.6 mm. A bolus of 140 mL of iopromide (Ultravist 370; Bayer Schering Pharma AG, Berlin, Germany) was administered at a rate of 3 mL/s through a 20-gauge angiographic catheter inserted into an antecubital vein. Portalvenous phase images were performed using a bolus tracking technique, 75 s delay after the aortic attenuation at the level of the diaphragm had reached 100 Hounsfield units (HU). Additional arterial phase was performed in 19/47 patients. Images with a transverse pixel size of 1.00 and slice thickness of 5 mm were reconstructed in the axial plane using a soft tissue kernel.

2.2. Image Segmentation, Radiomic Feature Extraction and Genomic Data

Image segmentation and radiomics features extraction was performed with a commonly used open-source software platform (LIFEx, Version 6.30; [14]) by one independent reader with 7 years of experience in oncologic radiology. The contours of the primary pancreatic tumor were manually delineated and segmented in venous phase of the preoperative CT in slice-by-slice fashion (Figure 1).
Next generation sequencing for whole genome sequencing was performed using Illumina HiSeq 2000/2500 reagents and instruments on paired-end library at the Ontario Cancer Research Institute (Toronto, ON, Canada). Data analysis protocol was based on previously described methods [15]. Illumina’s CASAVA software (version 1.8.2) converted the sequencing base calls to fastq format reads. Reads were aligned to the human reference genome (hg19_random) using Burrows-Wheeler Aligner (version 0.6.2) [16]. Somatic single nucleotide mutations were called using both Strelka (version 1.0.7) [17], and MuTect (version 1.1.4) [18], while indels were called using only Strelka (version 1.0.7). ANNOVAR [19] was used to annotate all the final mutation calls. Somatic copy number variation was assessed using CELLULOID [20].
Multiple genes, that have been previously reported as altered in PDCA were collected from a web-based resource (cBioPortal.ca for Cancer Genomics [21,22]) and were categorized as follows: significantly mutated genes, oncogenes, DNA damage repair genes and chromatin modification genes [7,8,23]. The 4 most prevalent genes were then included for statistical analysis (KRAS, TP53, SMAD4, CDKN2A) (Table 2).

2.3. Statistical Analysis and Modelling

2.3.1. Feature Selection

To reduce the number of radiomic features fed into the model building, a complete clustering algorithm was used to group correlated features to form a feature dendogram. Prior to clustering, features were normalized using a data-dependent normalization algorithm (R package bestNormalize, Version 1.8.3 [24]).
Features were selected from seven clusters, based on visual inspection of the cluster trees to strike a balance between cluster size and number of clusters. These clusters were independent of mutation status and clustering was performed once for the entire feature set. A reduced set of features was then selected for each mutation by choosing the feature from each cluster with the largest squared Spearman correlation with mutation status. Because our goal was to build a verifiable prediction model independent of our sample, we used raw values for the modelling of normally distributed features and a log transformation for those features with a mean/median ratio of >1.5. This step resulted in a reduced data set of seven features for each mutation.

2.3.2. Model Building

Prior to model building, 1000 cross-validation (CV) sets were selected for internal validation. Each set was comprised of 37 training participants and 10 testing participants, randomly allocated.
From the reduced set of features identified for each mutation all possible subset logistic regression (APS-LR) was performed to determine which combination of variables maximized the Youden Index (Sensitivity + Specificity − 1). APS-LR is computationally intensive, with processing time increasing exponentially with the number of variables under consideration and is possible because of the pre-selection of influential features. Two approaches were used for model building.
Method 1 involved performing APS-LR for each of the 1000 CV training sets and evaluating the model on the corresponding testing set. This approach allowed us to determine the most important features, defined as those appearing most frequently in the list of model predictors.
Method 2 involved performing the APS-LR step once, on the complete dataset to identify the model most predictive of each mutation, fixing the feature set for each mutation. The regression parameters from this model were then re-estimated on each of the 1000 CV training sets and evaluated on the testing sets.
Both methods allow us to inspect the distribution of the statistics associated with the model and therefore to estimate the distribution of the performance statistics, which cannot be achieved with a single testing/training set.

3. Results

3.1. Baseline Characteristics of the Study Cohort

Overall, 47 patients (21 females, 26 males; mean age 64 ± 11 years, range: 42–86 years) who received preoperative CT imaging and whole genome sequencing of the pancreatic resection specimen were included in our study. No systemic treatment was administered prior to surgery.
At the time of diagnosis, the majority of patients presented with advanced local disease with T3 or T4 PDAC in 42/47 (89%) patients, lymph node metastasis (N1) in 37/47 (78%) patients and M1 disease in 12/47 (26%) patients.
All patients were treated surgically, either with Whipple’s procedure (41/47; 87%) or distal pancreatectomy (6/47; 13%). The majority of patients underwent adjuvant chemotherapy (42/47; 89%) (Table 1).

3.2. Whole Genome Sequencing

More than half (30/47) of the analyzed tumor samples harbored 2 or more gene mutations. Overall, 81% of tumor samples demonstrated a missense mutation in KRAS, 68% of samples had alterations in TP53, 26% in SMAD4 and 19% in CDKN2A (Table 2).

3.3. Method 1

Method 1 revealed a set of important radiomics features for KRAS, TP53, SMAD4 and CDKN2, respectively (Table 3). Feature importance was calculated as the number of cross-validation models a feature appeared in. The most important features for each gene were as follows: HU_Skewness for both KRAS (number of cross-validation 918/1000) and TP53 (number of cross-validation 981/1000), GLZLM_SZLGE for SMAD 4 (number of cross-validation 743/1000) and NGLDM_Coarseness for CDKN2A (number of cross-validation 584/1000) (Table 3). The Youden Index for model 1 were as follows: 0.50 for KRAS, 0.62 for TP53, 0.44 for SMAD4 and 0.50 for CDKN2A (Table 4 and Figure 2).

3.4. Method 2

As compared to the model with varying radiomics features (method 1), the model performance slightly improved with a fixed set of the most important radiomics features (method 2). The most important features for each gene are displayed in Table 5. The Youden Index for method 2 were as follows: 0.56 for KRAS, 0.67 for TP53 and 0.50 for both SMAD4 and CDKN2A (Table 4 and Figure 2).

4. Discussion

In our study we assessed the association of CT-derived radiomics features with genetic driver mutations in surgical PDAC resection specimens. The main finding of our study demonstrates an acceptable correlation of the most important radiomics features with TP53 and KRAS gene mutation (Youden index of 0.67 and 0.56, respectively), indicating a reasonable predictive ability of genomic profiles using CT-derived radiomics.
Linking radiomics features to underlying tumor biology and pathogenesis is a developing field of increasing interest, given the wide availability of imaging data in contrast to costs, expenses and logistical issues of molecular analyses. While invasive tissue sampling remains the gold standard for histologic characterization, the usage of noninvasive imaging techniques for diagnosis and detection of specific tumor characteristics could represent a potential additive or eventually an alternative, especially in patients with advanced, inoperable disease or when inaccessible to biopsy. Biomarkers are continuously expanding and several actionable targets have already been identified in PDAC [25,26], partly already guiding clinical decision making and helping to develop novel targeted treatment strategies. Several genomic sequencing studies have also reported poor outcome in the presence of driver gene mutations in the PDAC resection specimen [27,28,29]. Although early stage precancerous lesions already appear to harbor driver genes mutations like KRAS and TP53, suggesting that these gene alterations play an important role in tumor onset and progression, genomic analysis of resection specimens represents only a snapshot of existing mutations [30,31].
Chemotherapy and combined chemo-radiotherapy is used both in adjuvant and neoadjuvant treatment of PDAC, however many patients are either intrinsically resistant or will acquire resistance to these treatments [32]. Therefore, it is crucial to understand related molecular mechanisms and pathogenesis, since specific gene mutations have been identified to contribute to reduced drug sensitivity and effectiveness of alternative therapeutic strategies. For example, a recent study showed that KRAS mutation plays an important role in the metabolic reprogramming of cancer cells, shifting them towards an anabolic metabolism and support unconstrained proliferation, leading to reduced activity of both gemcitabine and paclitaxel [33,34]. Wang et al. [35] also demonstrated resistance to radiotherapy in SMAD4-mutated PDAC through the induction of reactive oxygen species production and increased level of radiation-induced autophagy. This may indicate that the assessment of SMAD4 status represents a vital screening test, not only prior to radiotherapy but also for overall disease estimation, since multiple clinical studies showed that SMAD4 mutated PDAC commonly had distant disease progression as compared to those with intact SMAD4, who more commonly had local dominant disease pattern [36,37]. Several studies demonstrated that functional CDKN2A inactivation caused by mutations and deep deletions were associated with poor prognosis in PDAC patients [38,39,40]. Oshima et al. [41] investigated 106 surgically treated patients with PDAC and discovered that the loss of CDKN2A was significantly associated with lymphatic invasion and postoperative metastases. Another study, however, investigating chemosensitivity profiling of PDAC cell lines and patient-derived organoids, found that CDKN2A inactivation correlated with increased sensitivity to paclitaxel and an active metabolite of irinotecan [38].
In general, radiomics features can be classified into three main categories [42]. Shape features describe semantic and/or geometric properties (e.g., volume, maximum diameter). First order features describe the distribution of individual voxel values without any concern for their spatial relationships (e.g., skewness, kurtosis). Second order features describe the statistical relationship between neighboring voxels, providing a measure of the spatial arrangement of the voxel intensities and hence of intra-lesion heterogeneity (e.g., grey-level co-occurrence matrix, grey-level run length matrix).
Several studies so far investigated the correlation of radiomics feature and gene expression profile in various malignancies with a main focus on Non-Small Cell Lung Cancer (NSCLC) [43,44,45]. For example, Zhang et al. [44] constructed a clinical-radiological-radiomics model (C-R-R), based on the combination of CT radiomics feature signature with clinical and radiological features, aiming to predict epidermal growth factor receptor (EGFR) status among 420 patients with lung adenocarcinoma. The model demonstrated excellent diagnostic performance and high sensitivity in predicting EGFR mutation status. Another study by Digumarthy et al. [46] concluded that CT-derived radiomics features of NSCLC can help distinguish between EGFR positive and wild-type adenocarcinomas.
Further studies focused on the utility of CT radiomics features in PDAC, applied over a wide range of study purposes. Chu et al. [47] determine the utility of CT radiomics features in differentiating PDAC from normal pancreatic tissue in a cohort of 190 patients. The authors concluded that radiomics can be used to differentiate between CT cases from patients with PDAC and healthy control subjects with normal pancreas, using the maximally relevant-features, including gray-level-co-occurrence matrix, summation of the entropy and shape features of the whole pancreas. A similar work by Chen et al. [48] assessed the ability of CT radiomics features to distinguish between patients with and without PDAC in a large retrospective study, including 436 patients. The authors identified several features in the cancerous patch with lower values reflecting image intensity (first order: median, mean and 90th percentile) and higher values for features reflecting heterogeneity (NGTDM busyness, GLDM gray-level nonuniformity and GLDM dependence nonuniformity), consistent with the notion that PDACs typically manifest with heterogenous hypoenhancement on CT. Multiple other studies investigated survival prediction in PDAC using radiomics [49,50,51,52]. For example, Khalvati et al. [51] conducted a multicenter study, including 98 patients, aiming to assess the reproducibility and prognostic value of CT-derived radiomics features for resectable PDAC. The authors identified 2 significant second order features, among those, entropy-related features. Entropy measures the degree of randomness or non-uniformity in the image and it has been hypothesized that it can act as a surrogate for tumor heterogeneity. Notably, the majority of significant radiomics features, which were included in our prediction models consistent of either grey-level zone length matrix (GLZLM) or grey-level run length matrix (GLRLM), indicating spatial voxel distribution with estimation of tumor heterogeneity.
To the best of our knowledge there is only very limited literature investigating the correlation of CT radiomics features and gene expression profile in patients with PDAC. Attiyeh et al. [53] demonstrated a strong correlation of CT derived radiomics features with PDAC genetic alterations. In a cohort of 34 patients, 46% were found to have a mutation in SMAD4 and 83% showed alterations in TP53 expression. The authors found a good discriminatory power between patients with normal and abnormal SMAD4 status using CT-derived radiomics, whereas the radiomics feature model for TP53 mutations did not clearly discriminate between normal and those with altered TP53 status. The results of our study, which included a larger set of patients with PDAC, demonstrated superior predictive ability of TP53 mutation status (Youden Index 0.67). SMAD4 mutation status demonstrated only acceptable predictive signal (Youden Index 0.5), which may be explained by the relative low prevalence of SMAD4 mutations in our cohort (26% vs. 46%). Further, we followed a more holistic approach searching for more than twenty important gene mutations in PDAC, resulting in the inclusion of 4 main driver gene mutations (TP53, KRAS, SMAD4, CDKN2A) in our analysis. This may be of particular importance, since the interaction of synchronous gene mutations plays a crucial role in the development of drug resistances [54].
A further study by Iwatate et al. [55] evaluated in 107 patients with PDAC whether the expression of TP53 could be predicted using machine learning. The results of their study revealed the expression of TP53 to be predictable from CT texture features with a sensitivity of 0.667 and a specificity of 0.813 (AUC 0.795). The results of our analysis demonstrated higher sensitivity for predicting TP53 mutation. Our models revealed a maximal Youden Index of 0.67, with a sensitivity of 0.83 (95%CI 0.50–1.00) and a specificity of 0.50 (95%CI 0.00–1.00). However, the authors only investigated TP53, whereas we pursued a more holistic approach, including 4 main driver gene mutation in PDAC, which adds value to the literature.
As compared to the mentioned prior studies [53,55], we believe that we applied a more extended and more rigorous statistical analysis, reflecting a possible explanation for certain differences. It is common in high-dimensional studies, with many features relative to cases to perform split-validation, where the sample is divided into a training set (typically 80% of cases) and a distinct testing set. This strategy is used to ensure that the built model does not over-fit the data and to provide information about the reproducibility of the findings. Our strategy followed that of Papp et al. [56] whereby the sample was split into 1000 cross validation samples, each containing distinct training and testing sub-samples. We feel this approach is more robust because it enables us to obtain a distribution of the validation statistics and thereby shows the distribution of predictive accuracies, as opposed to a single feature estimate. While a single validation sample may indicate excellent discrimination the most clinically useful results arise when a large proportion of validation samples display high sensitivity and specificity. Additionally, the ability to place some confidence bounds around the reproducibility of the results is gained when multiple validation sets are used.
The following study limitations must be acknowledged. First, there are inherent drawbacks, due to the retrospective nature of the study and the relatively small sample size. Second, we investigated a relative inhomogeneous set of patients consisting of various tumor stages and grades, thus depicting different timepoints in the course of the disease. Third, due the relatively small sample size of our cohort, we did not discriminate between different molecular subtypes of PDAC, although we are aware of clinical and prognostic implications. Fourth, the Youden Index for the included mutations was 0.5–0.67, reflecting acceptable predictive ability, however given applied extended and rigorous statistical analysis, the respective values may demonstrate a reasonable result in this setting.

5. Conclusions

In conclusion, our study indicates acceptable correlation of CT-derived radiomics features and driver gene mutations in PDAC. Given small sample size and relative heterogeneity of the study cohort, we believe that a larger and more homogenous cohort may enhance the predictive ability.

Author Contributions

Conceptualization, R.H. and P.V.-H.; methodology, R.H. and P.V.-H.; validation, R.H. and P.V.-H.; formal analysis, L.A., R.H. and P.V.-H.; investigation, R.H.; resources, P.V.-H.; data curation, L.A., R.H. and P.V.-H.; writing—original draft preparation, R.H.; writing—review and editing, R.K., A.K., C.O., U.M., N.-A.P., M.H. and P.V.-H.; visualization, P.V.-H.; supervision, P.V.-H.; project administration, P.V.-H.; funding acquisition, P.V.-H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of the UNIVERSITY OF TORONTO (protocol code CAPCR/UHN REB #: 20-6105, 11/01/2021).

Informed Consent Statement

Patient consent was waived due to the retrospective nature of the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Elbanna, K.Y.; Jang, H.-J.; Kim, T.K. Imaging diagnosis and staging of pancreatic ductal adenocarcinoma: A comprehensive review. Insights Imaging 2020, 11, 58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2020. CA Cancer J. Clin. 2020, 70, 7–30. [Google Scholar] [CrossRef] [PubMed]
  3. Campen, C.J.; Dragovich, T.; Baker, A.F. Management strategies in pancreatic cancer. Am. J. Health Syst. Pharm. 2011, 68, 573–584. [Google Scholar] [CrossRef] [PubMed]
  4. Waddell, N.; Pajic, M.; Patch, A.-M.; Chang, D.K.; Kassahn, K.S.; Bailey, P.; Johns, A.L.; Miller, D.; Nones, K.; Quek, K. Whole genomes redefine the mutational landscape of pancreatic cancer. Nature 2015, 518, 495–501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Biankin, A.V.; Waddell, N.; Kassahn, K.S.; Gingras, M.-C.; Muthuswamy, L.B.; Johns, A.L.; Miller, D.K.; Wilson, P.J.; Patch, A.-M.; Wu, J. Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes. Nature 2012, 491, 399–405. [Google Scholar] [CrossRef] [Green Version]
  6. Jones, S.; Zhang, X.; Parsons, D.W.; Lin, J.C.-H.; Leary, R.J.; Angenendt, P.; Mankoo, P.; Carter, H.; Kamiyama, H.; Jimeno, A. Core signaling pathways in human pancreatic cancers revealed by global genomic analyses. Science 2008, 321, 1801–1806. [Google Scholar] [CrossRef] [Green Version]
  7. Bailey, P.; Chang, D.K.; Nones, K.; Johns, A.L.; Patch, A.-M.; Gingras, M.-C.; Miller, D.K.; Christ, A.N.; Bruxner, T.J.; Quinn, M.C. Genomic analyses identify molecular subtypes of pancreatic cancer. Nature 2016, 531, 47–52. [Google Scholar] [CrossRef]
  8. Raphael, B.J.; Hruban, R.H.; Aguirre, A.J.; Moffitt, R.A.; Yeh, J.J.; Stewart, C.; Robertson, A.G.; Cherniack, A.D.; Gupta, M.; Getz, G. Integrated genomic characterization of pancreatic ductal adenocarcinoma. Cancer Cell 2017, 32, 185–203.e113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Obermeyer, Z.; Emanuel, E.J. Predicting the Future—Big Data, Machine Learning, and Clinical Medicine. N. Engl. J. Med. 2016, 375, 1216–1219. [Google Scholar] [CrossRef] [Green Version]
  10. Baessler, B.; Nestler, T.; Dos Santos, D.P.; Paffenholz, P.; Zeuch, V.; Pfister, D.; Maintz, D.; Heidenreich, A. Radiomics allows for detection of benign and malignant histopathology in patients with metastatic testicular germ cell tumors prior to post-chemotherapy retroperitoneal lymph node dissection. Eur. Radiol. 2020, 30, 2334–2345. [Google Scholar] [CrossRef]
  11. GGillies, R.; Kinahan, P.; Hricak, H. Radiomics: Images Are More than Pictures, They Are Data. Radiology 2016, 278, 563–577. [Google Scholar] [CrossRef] [Green Version]
  12. Pinker, K.; Shitano, F.; Sala, E.; Do, R.K.; Young, R.J.; Wibmer, A.G.; Hricak, H.; Sutton, E.J.; Morris, E.A. Background, current role, and potential applications of radiogenomics. J. Magn. Reson. Imaging 2018, 47, 604–620. [Google Scholar] [CrossRef]
  13. Salinas-Miranda, E.; Healy, G.M.; Grünwald, B.; Jain, R.; Deniffel, D.; O’Kane, G.M.; Grant, R.; Wilson, J.; Knox, J.; Gallinger, S.; et al. Correlation of transcriptional subtypes with a validated CT radiomics score in resectable pancreatic ductal adenocarcinoma. Eur. Radiol. 2022, 32, 6712–6722. [Google Scholar] [CrossRef]
  14. Nioche, C.; Orlhac, F.; Boughdad, S.; Reuzé, S.; Goya-Outi, J.; Robert, C.; Pellot-Barakat, C.; Soussan, M.; Frouin, F.; Buvat, I. LIFEx: A Freeware for Radiomic Feature Calculation in Multimodality Imaging to Accelerate Advances in the Characterization of Tumor Heterogeneity. Cancer Res. 2018, 78, 4786–4789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Lathrop, M.; Gut, I.; Heath, S.; Tost, J.; Gress, T.; Hudson, T. International network of cancer genome projects (The International Cancer Genome Consortium). Nature 2010, 464, 993–998. [Google Scholar]
  16. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Saunders, C.T.; Wong, W.S.; Swamy, S.; Becq, J.; Murray, L.J.; Cheetham, R.K. Strelka: Accurate somatic small-variant calling from sequenced tumor–normal sample pairs. Bioinformatics 2012, 28, 1811–1817. [Google Scholar] [CrossRef] [Green Version]
  18. Cibulskis, K.; Lawrence, M.S.; Carter, S.L.; Sivachenko, A.; Jaffe, D.; Sougnez, C.; Gabriel, S.; Meyerson, M.; Lander, E.S.; Getz, G. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat. Biotechnol. 2013, 31, 213–219. [Google Scholar] [CrossRef]
  19. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef] [PubMed]
  20. Notta, F.; Chan-Seng-Yue, M.; Lemire, M.; Li, Y.; Wilson, G.W.; Connor, A.A.; Denroche, R.E.; Liang, S.-B.; Brown, A.M.; Kim, J.C. A renewed model of pancreatic cancer evolution based on genomic rearrangement patterns. Nature 2016, 538, 378–382. [Google Scholar] [CrossRef] [Green Version]
  21. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E. The cBio cancer genomics portal: An open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [PubMed]
  22. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci. Signal. 2013, 6, 269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Pham, N.-A.; Radulovich, N.; Ibrahimov, E.; Martins-Filho, S.N.; Li, Q.; Pintilie, M.; Weiss, J.; Raghavan, V.; Cabanero, M.; Denroche, R.E. Patient-derived tumor xenograft and organoid models established from resected pancreatic, duodenal and biliary cancers. Sci. Rep. 2021, 11, 10619. [Google Scholar] [CrossRef] [PubMed]
  24. Peterson, R.A. Finding Optimal Normalizing Transformations via best Normalize. R J. 2021, 13, 310–329. [Google Scholar] [CrossRef]
  25. Golan, T.; Hammel, P.; Reni, M.; Van Cutsem, E.; Macarulla, T.; Hall, M.J.; Park, J.-O.; Hochhauser, D.; Arnold, D.; Oh, D.-Y. Maintenance olaparib for germline BRCA-mutated metastatic pancreatic cancer. N. Engl. J. Med. 2019, 381, 317–327. [Google Scholar] [CrossRef]
  26. O’Reilly, E.M.; Lee, J.W.; Lowery, M.A.; Capanu, M.; Stadler, Z.K.; Moore, M.J.; Dhani, N.; Kindler, H.L.; Estrella, H.; Maynard, H. Phase 1 trial evaluating cisplatin, gemcitabine, and veliparib in 2 patient cohorts: Germline BRCA mutation carriers and wild-type BRCA pancreatic ductal adenocarcinoma. Cancer 2018, 124, 1374–1382. [Google Scholar] [CrossRef] [Green Version]
  27. Qian, Z.R.; Rubinson, D.A.; Nowak, J.A.; Morales-Oyarvide, V.; Dunne, R.F.; Kozak, M.M.; Welch, M.W.; Brais, L.K.; Da Silva, A.; Li, T. Association of alterations in main driver genes with outcomes of patients with resected pancreatic ductal adenocarcinoma. JAMA Oncol. 2018, 4, e173420. [Google Scholar] [CrossRef]
  28. McIntyre, C.A.; Lawrence, S.A.; Richards, A.L.; Chou, J.F.; Wong, W.; Capanu, M.; Berger, M.F.; Donoghue, M.T.; Yu, K.H.; Varghese, A.M. Alterations in driver genes are predictive of survival in patients with resected pancreatic ductal adenocarcinoma. Cancer 2020, 126, 3939–3949. [Google Scholar] [CrossRef]
  29. Miyabayashi, K.; Nakagawa, H.; Koike, K. Molecular and phenotypic profiling for precision medicine in pancreatic cancer: Current advances and future perspectives. Front. Oncol. 2021, 11, 2373. [Google Scholar] [CrossRef]
  30. Fischer, C.G.; Wood, L.D. From somatic mutation to early detection: Insights from molecular characterization of pancreatic cancer precursor lesions. J. Pathol. 2018, 246, 395–404. [Google Scholar] [CrossRef] [Green Version]
  31. Bernard, V.; Semaan, A.; Huang, J.; San Lucas, F.A.; Mulu, F.C.; Stephens, B.M.; Guerrero, P.A.; Huang, Y.; Zhao, J.; Kamyabi, N. Single-Cell Transcriptomics of Pancreatic Cancer Precursors Demonstrates Epithelial and Microenvironmental Heterogeneity as an Early Event in Neoplastic ProgressionSingle-Cell RNA Sequencing of Pancreatic Cancer. Clin. Cancer Res. 2019, 25, 2194–2205. [Google Scholar] [CrossRef]
  32. Yoneyama, H.; Takizawa-Hashimoto, A.; Takeuchi, O.; Watanabe, Y.; Atsuda, K.; Asanuma, F.; Yamada, Y.; Suzuki, Y. Acquired resistance to gemcitabine and cross-resistance in human pancreatic cancer clones. Anti Cancer Drugs 2015, 26, 90–100. [Google Scholar] [CrossRef] [PubMed]
  33. 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]
  34. Kimmelman, A.C. Metabolic dependencies in RAS-driven cancers. Clin. Cancer Res. 2015, 21, 1828–1834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wang, F.; Xia, X.; Yang, C.; Shen, J.; Mai, J.; Kim, H.-C.; Kirui, D.; Kang, Y.A.; Fleming, J.B.; Koay, E.J. SMAD4 gene mutation renders pancreatic cancer resistance to radiotherapy through promotion of autophagy. Clin. Cancer Res. 2018, 24, 3176–3185. [Google Scholar] [CrossRef] [Green Version]
  36. Tascilar, M.; Skinner, H.G.; Rosty, C.; Sohn, T.; Wilentz, R.E.; Offerhaus, G.J.A.; Adsay, V.; Abrams, R.A.; Cameron, J.L.; Kern, S.E. The SMAD4 protein and prognosis of pancreatic ductal adenocarcinoma. Clin. Cancer Res. 2001, 7, 4115–4121. [Google Scholar]
  37. Crane, C.H.; Varadhachary, G.R.; Yordy, J.S.; Staerkel, G.A.; Javle, M.M.; Safran, H.; Haque, W.; Hobbs, B.D.; Krishnan, S.; Fleming, J.B. Phase II trial of cetuximab, gemcitabine, and oxaliplatin followed by chemoradiation with cetuximab for locally advanced (T4) pancreatic adenocarcinoma: Correlation of Smad4 (Dpc4) immunostaining with pattern of disease progression. J. Clin. Oncol. 2011, 29, 3037. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Lin, J.-C.; Liu, T.-P.; Yang, P.-M. CDKN2A-inactivated pancreatic ductal adenocarcinoma exhibits therapeutic sensitivity to paclitaxel: A bioinformatics study. J. Clin. Med. 2020, 9, 4019. [Google Scholar] [CrossRef] [PubMed]
  39. Chen, J.; Li, D.; Killary, A.M.; Sen, S.; Amos, C.I.; Evans, D.B.; Abbruzzese, J.L.; Frazier, M.L. Polymorphisms of p16, p27, p73, and MDM2 modulate response and survival of pancreatic cancer patients treated with preoperative chemoradiation. Ann. Surg. Oncol. 2009, 16, 431–439. [Google Scholar] [CrossRef]
  40. Luo, Y.; Tian, L.; Feng, Y.; Yi, M.; Chen, X.; Huang, Q. The predictive role of p16 deletion, p53 deletion, and polysomy 9 and 17 in pancreatic ductal adenocarcinoma. Pathol. Oncol. Res. 2013, 19, 35–40. [Google Scholar] [CrossRef]
  41. Oshima, M.; Okano, K.; Muraki, S.; Haba, R.; Maeba, T.; Suzuki, Y.; Yachida, S. Immunohistochemically detected expression of 3 major genes (CDKN2A/p16, TP53, and SMAD4/DPC4) strongly predicts survival in patients with resectable pancreatic cancer. Ann. Surg. 2013, 258, 336–346. [Google Scholar] [CrossRef]
  42. Marti-Bonmati, L.; Cerdá-Alberich, L.; Pérez-Girbés, A.; Díaz Beveridge, R.; Montalvá Orón, E.; Pérez Rojas, J.; Alberich-Bayarri, A. Pancreatic cancer, radiomics and artificial intelligence. Br. J. Radiol. 2022, 95, 20220072. [Google Scholar] [CrossRef] [PubMed]
  43. Avesani, G.; Tran, H.E.; Cammarata, G.; Botta, F.; Raimondi, S.; Russo, L.; Persiani, S.; Bonatti, M.; Tagliaferri, T.; Dolciami, M. CT-Based Radiomics and Deep Learning for BRCA Mutation and Progression-Free Survival Prediction in Ovarian Cancer Using a Multicentric Dataset. Cancers 2022, 14, 2739. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, G.; Cao, Y.; Zhang, J.; Ren, J.; Zhao, Z.; Zhang, X.; Li, S.; Deng, L.; Zhou, J. Predicting EGFR mutation status in lung adenocarcinoma: Development and validation of a computed tomography-based radiomics signature. Am. J. Cancer Res. 2021, 11, 546–560. [Google Scholar] [PubMed]
  45. Veeraraghavan, H.; Friedman, C.F.; DeLair, D.F.; Ninčević, J.; Himoto, Y.; Bruni, S.G.; Cappello, G.; Petkovska, I.; Nougaret, S.; Nikolovski, I. Machine learning-based prediction of microsatellite instability and high tumor mutation burden from contrast-enhanced computed tomography in endometrial cancers. Sci. Rep. 2020, 10, 17769. [Google Scholar] [CrossRef] [PubMed]
  46. Digumarthy, S.R.; Padole, A.M.; Gullo, R.L.; Sequist, L.V.; Kalra, M.K. Can CT radiomic analysis in NSCLC predict histology and EGFR mutation status? Medicine 2019, 98, e13963. [Google Scholar] [CrossRef]
  47. Chu, L.C.; Park, S.; Kawamoto, S.; Fouladi, D.F.; Shayesteh, S.; Zinreich, E.S.; Graves, J.S.; Horton, K.M.; Hruban, R.H.; Yuille, A.L. Utility of CT radiomics features in differentiation of pancreatic ductal adenocarcinoma from normal pancreatic tissue. Am. J. Roentgenol. 2019, 213, 349–357. [Google Scholar] [CrossRef]
  48. Chen, P.-T.; Chang, D.; Yen, H.; Liu, K.-L.; Huang, S.-Y.; Roth, H.; Wu, M.-S.; Liao, W.-C.; Wang, W. Radiomic features at CT can distinguish pancreatic cancer from noncancerous pancreas. Radiol. Imaging Cancer 2021, 3, e210010. [Google Scholar] [CrossRef] [PubMed]
  49. Park, S.; Sham, J.G.; Kawamoto, S.; Blair, A.B.; Rozich, N.; Fouladi, D.F.; Shayesteh, S.; Hruban, R.H.; He, J.; Wolfgang, C.L. CT Radiomics–Based Preoperative Survival Prediction in Patients with Pancreatic Ductal Adenocarcinoma. Am. J. Roentgenol. 2021, 217, 1104–1112. [Google Scholar] [CrossRef]
  50. Gao, Y.; Cheng, S.; Zhu, L.; Wang, Q.; Deng, W.; Sun, Z.; Wang, S.; Xue, H. A systematic review of prognosis predictive role of radiomics in pancreatic cancer: Heterogeneity markers or statistical tricks? Eur. Radiol. 2022, 32, 8443–8452. [Google Scholar] [CrossRef]
  51. Khalvati, F.; Zhang, Y.; Baig, S.; Lobo-Mueller, E.M.; Karanicolas, P.; Gallinger, S.; Haider, M.A. Prognostic value of CT radiomic features in resectable pancreatic ductal adenocarcinoma. Sci. Rep. 2019, 9, 5449. [Google Scholar] [CrossRef] [Green Version]
  52. Xie, T.; Wang, X.; Li, M.; Tong, T.; Yu, X.; Zhou, Z. Pancreatic ductal adenocarcinoma: A radiomics nomogram outperforms clinical model and TNM staging for survival estimation after curative resection. Eur. Radiol. 2020, 30, 2513–2524. [Google Scholar] [CrossRef] [PubMed]
  53. Attiyeh, M.A.; Chakraborty, J.; McIntyre, C.A.; Kappagantula, R.; Chou, Y.; Askan, G.; Seier, K.; Gonen, M.; Basturk, O.; Balachandran, V.P. CT radiomics associations with genotype and stromal content in pancreatic ductal adenocarcinoma. Abdom. Radiol. 2019, 44, 3148–3157. [Google Scholar] [CrossRef] [PubMed]
  54. Pao, W.; Wang, T.Y.; Riely, G.J.; Miller, V.A.; Pan, Q.; Ladanyi, M.; Zakowski, M.F.; Heelan, R.T.; Kris, M.G.; Varmus, H.E. KRAS mutations and primary resistance of lung adenocarcinomas to gefitinib or erlotinib. PLoS Med. 2005, 2, e17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Iwatate, Y.; Hoshino, I.; Yokota, H.; Ishige, F.; Itami, M.; Mori, Y.; Chiba, S.; Arimitsu, H.; Yanagibashi, H.; Nagase, H.; et al. Radiogenomics for predicting p53 status, PD-L1 expression, and prognosis with machine learning in pancreatic cancer. Br. J. Cancer 2020, 123, 1253–1261. [Google Scholar] [CrossRef] [PubMed]
  56. Papp, L.; Spielvogel, C.; Grubmüller, B.; Grahovac, M.; Krajnc, D.; Ecsedi, B.; Sareshgi, R.; Mohamad, D.; Hamboeck, M.; Rausch, I. Supervised machine learning enables non-invasive lesion characterization in primary prostate cancer with [68Ga] Ga-PSMA-11 PET/MRI. Eur. J. Nucl. Med. Mol. Imaging 2021, 48, 1795–1805. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Axial CT image (a) and triplanar CT images with segmentation (displayed in pink; b, c, d) in a representative 62-year-old woman with PDAC in the pancreatic head.
Figure 1. Axial CT image (a) and triplanar CT images with segmentation (displayed in pink; b, c, d) in a representative 62-year-old woman with PDAC in the pancreatic head.
Cancers 14 06224 g001
Figure 2. Boxplots demonstrating the Youden Index comparing method 1 (upper row) and method 2 (lower row).
Figure 2. Boxplots demonstrating the Youden Index comparing method 1 (upper row) and method 2 (lower row).
Cancers 14 06224 g002
Table 1. Demographic data of the study cohort.
Table 1. Demographic data of the study cohort.
Characteristicsn = 47
Age (mean ± SD; range)63.7 ± 10.7 (42–86)
Sex
Females21 (45%)
Male26 (55%)
Current or former smoker21 (45%)
Race
Asian6 (13%)
Non-Asian41 (87%)
Tumor Grade
G18 (17%)
G226 (55%)
G312 (26%)
G41 (2%)
TNM
T25 (11%)
T3-442 (89%)
N137 (78%)
M18 (17%)
Treatment
Surgery47 (100%)
Whipple’s procedure41 (87%)
Distal Pancreatectomy6 (13%)
Adjuvant Chemotherapy42 (89%)
Table 2. Prevalence of genomic alterations in PDAC resection specimens. The 4 most prevalent gene mutations are presented in bold.
Table 2. Prevalence of genomic alterations in PDAC resection specimens. The 4 most prevalent gene mutations are presented in bold.
GenesPrevalence of Mutations N (%)
Significantly mutated genes
KRAS38 (81)
TP5332 (68)
SMAD412 (26)
CDKN2A9 (19)
ARID1A2 (4)
TGFBR21 (2)
NF11 (2)
Oncogenes
BRAF1 (2)
EGFR1 (2)
ERBB21 (2)
FGFR11 (2)
GNAS1 (2)
MET1 (2)
PAK41 (2)
DNA damage repair genes
BRCA21 (2)
ATM1 (2)
Tumor suppressor genes
RNF432 (4)
ACVR2A2 (4)
SMAD31 (2)
TSC21 (2)
Chromatin modification genes
KDM6A1 (2)
KMT2C3 (6)
PBRM11 (2)
SMARCA21 (2)
SMARCA41 (2)
Table 3. Feature importance of method 1 for each gene mutation with a varying set of features. Feature importance was calculated as the number of cross-validation models a feature appeared in 1000 sets.
Table 3. Feature importance of method 1 for each gene mutation with a varying set of features. Feature importance was calculated as the number of cross-validation models a feature appeared in 1000 sets.
FeatureKRASTP53SMAD4CDKN2A
HU_Skewness918981
HU_Q3860
GLZLM_LZLGE690
GLRLM_LRLGE664
GLRLM_SRHGE578
GLRLM_LGRE514
NGLDM_Contrast470
GLZLM_SZHGE 782
NGLDM_Coarseness 752423584
GLCM_Energy 617
GLZLM_ZLNU 502
PARAMS_ZSpatialResampling 353
HU_max 351
GLZLM_SZLGE 743397
HU_peakSphere1mL 633
HUmin 541465
GLZLM_GLNU 528561
HU_Q1 491580
SHAPE_Volume 396
GLCM_Correlation 526
SHAPE_Sphericity 461
Table 4. Statistical parameters (median and 95% CI) for method 1 and method 2 across the cross-validation samples.
Table 4. Statistical parameters (median and 95% CI) for method 1 and method 2 across the cross-validation samples.
KRASTP53SMAD4CDKN2A
Model 1
NPV0.00 (0.00, 1.00)0.50 (0.00, 1.00)0.75 (0.44, 1.00)0.80 (0.56, 1.00)
PPV0.86 (0.56, 1.00)0.75 (0.43, 1.00)0.00 (0.00, 0.50)0.00 (0.00, 0.50)
Sensitivity0.88 (0.57, 1.00)0.80 (0.44, 1.00)0.00 (0.00, 0.67)0.00 (0.00, 0.50)
Specificity0.00 (0.00, 1.00)0.50 (0.00, 1.00)0.86 (0.50, 1.00)0.89 (0.56, 1.00)
Youden Index0.50 (0.31, 0.94)0.62 (0.31, 0.94)0.44 (0.28, 0.72)0.50 (0.28, 0.69)
Model 2
NPV0.33 (0.00, 1.00)0.50 (0.00, 1.00)0.75 (0.50, 1.00)0.80 (0.56, 1.00)
PPV0.88 (0.62, 1.00)0.78 (0.43, 1.00)0.00 (0.00, 1.00)0.00 (0.00, 1.00)
Sensitivity0.88 (0.62, 1.00)0.83 (0.50, 1.00)0.00 (0.00, 0.50)0.00 (0.00, 0.50)
Specificity0.25 (0.00, 1.00)0.50 (0.00, 1.00)0.89 (0.57, 1.00)0.89 (0.56, 1.00)
Youden Index0.56 (0.33, 0.97)0.67 (0.33, 0.94)0.50 (0.31, 0.69)0.50 (0.31, 0.75)
Table 5. Feature importance of method 2 for each gene mutation with a fixed set of features based on the complete sample.
Table 5. Feature importance of method 2 for each gene mutation with a fixed set of features based on the complete sample.
FeatureKRAS Mutations N (%)TP53SMAD4CDKN2A
HUSkewnessxx
HU_Q3x
GLZLM_LZLGEx
GLRLM_LRLGEx
NGLDM_Contrastx
GLZLM_SZHGE x
NGLDM_Coarseness xxx
GLCM_Energy x
GLZLM_ZLNU x
GLZLM_SZLGE xx
HUpeakSphere1mL xx
GLZLM_GLNU xx
HU_Q1 x
GLCM_Correlation x
SHAPE_Sphericity x
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hinzpeter, R.; Kulanthaivelu, R.; Kohan, A.; Avery, L.; Pham, N.-A.; Ortega, C.; Metser, U.; Haider, M.; Veit-Haibach, P. CT Radiomics and Whole Genome Sequencing in Patients with Pancreatic Ductal Adenocarcinoma: Predictive Radiogenomics Modeling. Cancers 2022, 14, 6224. https://doi.org/10.3390/cancers14246224

AMA Style

Hinzpeter R, Kulanthaivelu R, Kohan A, Avery L, Pham N-A, Ortega C, Metser U, Haider M, Veit-Haibach P. CT Radiomics and Whole Genome Sequencing in Patients with Pancreatic Ductal Adenocarcinoma: Predictive Radiogenomics Modeling. Cancers. 2022; 14(24):6224. https://doi.org/10.3390/cancers14246224

Chicago/Turabian Style

Hinzpeter, Ricarda, Roshini Kulanthaivelu, Andres Kohan, Lisa Avery, Nhu-An Pham, Claudia Ortega, Ur Metser, Masoom Haider, and Patrick Veit-Haibach. 2022. "CT Radiomics and Whole Genome Sequencing in Patients with Pancreatic Ductal Adenocarcinoma: Predictive Radiogenomics Modeling" Cancers 14, no. 24: 6224. https://doi.org/10.3390/cancers14246224

APA Style

Hinzpeter, R., Kulanthaivelu, R., Kohan, A., Avery, L., Pham, N. -A., Ortega, C., Metser, U., Haider, M., & Veit-Haibach, P. (2022). CT Radiomics and Whole Genome Sequencing in Patients with Pancreatic Ductal Adenocarcinoma: Predictive Radiogenomics Modeling. Cancers, 14(24), 6224. https://doi.org/10.3390/cancers14246224

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