Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 17 May 2021
Sec. Cancer Imaging and Image-directed Interventions
This article is part of the Research Topic Radiomics-Based Tumor Phenotyping in Precision Medicine View all 44 articles

Differentiating TP53 Mutation Status in Pancreatic Ductal Adenocarcinoma Using Multiparametric MRI-Derived Radiomics

Jing Gao&#x;Jing Gao1†Xiahan Chen&#x;Xiahan Chen2†Xudong Li,Xudong Li1,3Fei MiaoFei Miao4Weihuan FangWeihuan Fang5Biao LiBiao Li1Xiaohua Qian*Xiaohua Qian2*Xiaozhu Lin*Xiaozhu Lin1*
  • 1Department of Nuclear Medicine, Ruijin Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai, China
  • 2School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
  • 3Department of Nuclear Medicine, Qingdao Municipal Hospital, Qingdao, China
  • 4Department of Radiology, Ruijin Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai, China
  • 5Department of Radiology, Ruijin Hospital North, Shanghai, China

Objectives: This study assessed the preoperative prediction of TP53 status based on multiparametric magnetic resonance imaging (mpMRI) radiomics extracted from two-dimensional (2D) and 3D images.

Methods: 57 patients with pancreatic cancer who underwent preoperative MRI were included. The diagnosis and TP53 gene test were based on resections. Of the 57 patients included 37 mutated TP53 genes and the remaining 20 had wild-type TP53 genes. Two radiologists performed manual tumour segmentation on seven different MRI image acquisition sequences per patient, including multi-phase [pre-contrast, late arterial phase (ap), portal venous phase, and delayed phase] dynamic contrast enhanced (DCE) T1-weighted imaging, T2-weighted imaging (T2WI), Diffusion-weighted imaging (DWI), and apparent diffusion coefficient (ADC). PyRadiomics-package was used to generate 558 two-dimensional (2D) and 994 three-dimensional (3D) image features. Models were constructed by support vector machine (SVM) for differentiating TP53 status and DX score method were used for feature selection. The evaluation of the model performance included area under the curve (AUC), accuracy, calibration curves, and decision curve analysis.

Results: The 3D ADC-ap-DWI-T2WI model with 11 selected features yielded the best performance for differentiating TP53 status, with accuracy = 0.91 and AUC = 0.96. The model showed the good calibration. The decision curve analysis indicated that the radiomics model had clinical utility.

Conclusions: A non-invasive and quantitative mpMRI-based radiomics model can accurately predict TP53 mutation status in pancreatic cancer patients and contribute to the precision treatment.

Introduction

Pancreatic ductal adenocarcinoma (PDAC) is characterized by late diagnosis, high mortality rate, and low overall survival (1). The poor prognosis and inefficiency of current treatments are primarily caused by chemoresistance for PDAC. The absence of key genetic alterations is the main driver of chemoresistance, which can disorder the apoptotic (2). As one of the major genetic mutations,TP53 mutations were seen in 70% of pancreatic cancer (3). TP53 encodes the p53 protein, as a tumour suppressor gene, which restricts cell proliferation in many cellular processes, involved in DNA repair, cell-cycle arrest, and apoptosis (4). These mutations are linked to poor patient prognosis (5), although literature regarding its influence is controversial. In addition, resistance to some therapeutic modalities, such as gemcitabine and 5-fluorouracil-based therapies (6, 7). However, novel regimens that target pancreatic cancer cells are emerging. TP53 may be an attractive target for gene augmentation therapy in pancreatic cancer (8).

Unlike surgery resection or biopsy, radiomics as a non-invasive tool was used to detect TP53 mutations (9), which can assess tumour heterogeneity by evaluating the grey-level intensity of pixels and their position in a medical image (10). Currently, CT-based radiomics have been widely used for predicting gene expression and survival prediction (11, 12) in PDAC patients. Some studies use MRI texture recognizing the status of TP53 in many cancers (13, 14). Standard methods of texture analysis methods are 2D or 3D approaches. 3D imaging features may capture tissue properties of the entire tissue more accurately, which improves the predictive power of imaging biomarkers in pancreatic cancer. However, 3D whole-tumour analysis is complicated and time-consuming. Although 2D or 3D texture analysis (TA) has been previously used to extract CT/MRI image features to predict gene status (15), to our best knowledge, neither 2D or 3D TA of mpMRI has been aimed at predicting TP53 status before.

In this study, support vector machine (SVM) radiomics models were constructed using 2D and 3D texture features extracted from mpMRI for the assessment of pre-operative TP53 mutation in PDAC patients.

Materials and Methods

Study Population and Tissue Samples

Following local Institutional Review Board approval, this retrospective study was approved with a waiver to obtain written informed consent. Patients with pancreatic ductal adenocarcinoma (PDAC) treated with surgery-based strategy at Ruijin Hospital from January 2016 to December 2016 were included in this study. The inclusion criteria were as follows: 1) pathologic confirmation of PDAC; 2) available the next-generation sequencing(NGS)-based TP53 sequence analysis; 3) MR images contained all of the following sequences: multi-phase [pre-contrast, late arterial phase (ap), portal venous phase (pp), and delayed phase (dp)] dynamic contrast enhanced (DCE) T1-weighted imaging (T1WI, T1), T2-weighted imaging (T2WI, T2), DWI and ADC.

Tissue samples from surgical resections of all 57 PDAC patients were analyzed. Genomic DNA was extracted from formalin-fixed paraffin-embedded (FFPE) tissue, and TP53 mutations were examined using NGS approach.

Image Acquisition

MR images were acquired on 1.5 T MRI scanners (N=9) or 3.0 T MRI scanners (N=48). The MRI examination included different acquisition sequences, including axial turbo spin-echo T2 sequence with fat saturation, DWI using a single-shot echo-planar imaging pulse sequence with b-values (0, 600 or 800 s/mm2), pre- and post-contrast fat-suppressed T1-weighted gradient sequences with intravenous administration of gadopentetate dimeglumine-diethylenetriaminepentaacetic acid (Gd-DTPA) contrast. The apparent diffusion coefficient was calculated by using a monoexponential function with b-values of 0 s/mm2 and 600 or 800 s/mm2. Gd-DPTA dose of 0.1 mmol per kg and a flow rate of 2 mL per sec were achieved. Late arterial phase, portal vein phase and delayed phase were acquired approximately 35-, 60-, and 90- seconds after contrast injection. Scan parameters for the MRI sequences are summarized in the Supplementary Table 1.

Segmentation of Region of Interest (ROI)

ROIs of the tumour were manually segmented by a junior radiologist in ITK‐SNAP software [Version 3.6 (16)] and were validated by an experienced senior radiologist. ROIs were manually drawn along the margin of the tumor covering the largest possible region. ROIs were delineated on multi-phase (pre-contrast, ap, pp, dp) dynamic contrast enhanced T1WI, T2WI, DWI and ADC images.

MRI Image Feature Extraction

A total number of 558 2D image features from the largest cross-sectional area of a tumor and 994 3D image features were extracted from the entire tumor area for each image. In order to improve multiparametric MRI radiomic feature robustness, the image intensity of each sequence was normalized to the range of 0–1. MRI voxel was resampled to 1 mm× 1 mm× 1 mm to reduce the variability of different scanners. The features included first-order features, shape features, gray level co-occurrence matrix (GLCM) features, gray level dependence matrix (GLDM) features, gray level run length matrix (GLRLM) features, gray level size zone matrix (GLSZM) features and neighbouring gray-tone difference matrix (NGTDM) features of original images, wavelet transformed images and gradient images. The specific number of features are listed in Supplementary Table 2. The feature extraction procedure was implemented in the Pyradiomics package (python 3.6) (17).

Feature Selection and SVM Model Construction

In total, 378 radiomics models (i.e. 126×3 = 378) were formed based on 2D,3D, and 2D/3D combination from the seven different mpMRI datasets. For each model, feature selection and TP53 gene prediction were performed separately.

For 2D features, the number of models can be calculated by :

number=C71+C72+C73+C74+C75+C76+C77=126

Where

Cnm=n!m!(nm)!

Among the extracted features, some were highly correlated and some had poor ability to assess TP53 gene mutation. Besides, it’s unknown which model features play the main role. Therefore, in order to remove the most redundant and irrelevant features and choose the most important and typical features, we performed feature selection base on DX score method (18) before gene prediction. DX score is an effective method to measure the difference between positive and negative samples (19). The higher the score is, the stronger the discriminating ability to distinguish between two types of samples. It can be mathematically defined as:

D(X)=(mpositivemnegative)2dpositive2+dnegative2(1)

where mpositive and mnegative are the mean value respectively, dpositive and dnegative are the standard deviation of the feature X.

m=1Ni=1Nxi(2)
d=1N1i=1N(xim)2(3)

where x is the value of the feature X with respect to positive (or negative) samples, N is the total number of positive (or negative) samples.

Then, features were ranked from the most important to the least important. The top n(1≤n≤N) features were input to the SVM classifier model of the initial parameter in turn. We assessed its classification performance using the accuracy rate, so the number of features in the feature set is the abscissa (X-axis), and the classification accuracy of each feature set is plotted on the ordinate (Y-axis). The feature set with the highest accuracy was selected.

Next, the SVM model was constructed using the features selected above. A grid search was conducted on the trade-off coefficient (C) and the kernel function parameter (gamma). Then, features were selected again using the optimal SVM model. The SVM package LIBSVM (20) was used for SVM model construction due to its well-known performance. Five-fold cross-validation was applied to evaluate model performance during the experiment to avoid over-fitting.

Statistical Analysis

Independent sample t tests and chi-squared test were used to compare continuous variables and categorical variables, respectively. Receiver operating characteristic (ROC) curve, calibration curve, and decision curve analysis (DCA) were used to evaluate and choose between the generated models. The corresponding values of the area under the ROC curve (AUC) and the Brier score (BS) were also calculated. The higher the AUC value, the better the model performance. On the other hand, the Brier score is the measure of the calibration curve. For a set of prediction values, the lower the Brier score, the better the prediction calibration. BS score is mathematically calculated as:

BS=14t=1N(ptot)2(4)

where pt is the probability of prediction, ot is the real probability of sample t and N is the number of samples. The calibration curve and Brier score were obtained with Sklearn (python 3.6). A t-test was conducted to evaluate significant differences between the two models. Decision curve analysis was used to quantify the net benefit of models to guide subsequent actions (21). The ROC and AUC were used to evaluate the diagnostic performance of radiomics models. The ROCs of radiomics models were compared using the DeLong test (22) implemented in MATLAB R2018a(Mathworks, Natick, MA, USA). A p value < 0.05 was considered statistically different.

Results

Patients

Out of 57 patients, 37 had mutated TP53 genes, and 20 had wild-type TP53 genes. The characteristics of patients are summarized in Table 1. There were no statistical differences in age (p = 0.770), gender (p = 0.397), histologic stage (p = 0.402) between the TP53 mutation and TP53 wild-type groups.

TABLE 1
www.frontiersin.org

Table 1 Clinical and pathological analysis of patients with or without TP53 mutation.

Radiomics Model Performance and Feature Selection

In total, 378 radiomics models were formed based on the combination of 7 sequences in the 2D, 3D, and 2D/3D combination. The model including ADC, ap (DCE T1WI late arterial phase), dp (DCE T1WI delayed phase), DWI, pp (DCE T1WI portal venous phase), T2 (T2 weighted imaging) sequence (ADC_ap_dp_DWI_pp_T2) generated the best performance with an AUC of 98.02%, including 41 features. Models with fewer predictors and similar predictive performance from all combinations of modalities were also selected, to avoid over-fitting and redundancy. Comparing the AUCs from the above-mentioned model presenting the highest value of AUC revealed that there were seven models with no statistically significant difference (P>0.05) in their abilities to differentiate (Table 2). The models with the top two AUC values and the least number of features were selected and named as the model I, II, and III. Details of extracted radiomic features in the three model was listed in the Table 3. The overall ability of model I, model II, and model III for identifying TP53 status are represented by their AUC values (96.15,94.81,87.85%, respectively; Figure 1A). Furthermore, Figure 1B shows the Brier score, measuring overall model performance via mean squared error between predicted probabilities and expected values.

TABLE 2
www.frontiersin.org

Table 2 Comparison of Performance of the models between the highest one versus the other models with fewer predictors.

TABLE 3
www.frontiersin.org

Table 3 Details of extracted radiomics features in the three models.

FIGURE 1
www.frontiersin.org

Figure 1 The ROC curves and calibration curves of three classification models. (A) the 3D-ADC-ap-DWI-T2 model (best one, AUC=0.9615), 3D-ADC-DWI-pp-T2 model (AUC=0.9481) and 3D-ADC-ap-pp-T2 model including the fewest features model (AUC=0.8786). (B) Observed (y-axis) versus the predicted probability frequency (x-axis). The closer the points appear along the main diagonal, the better calibrated. 3D-ADC-ap-DWI-T2 is the closet to the diagonal dotted line, which represents perfect calibration.

In our study, the 3D-ADC-ap-DWI-T2 model had the highest AUC and the lowest Brier score, indicating superior model performance. The optimal feature subset selected by DX score based on SVM classifier accuracy (using 5-fold CV accuracy), included the top-11 feature set for predicting TP53. Figure 2 illustrates the feature selection process for the best model (3D-ADC-ap-DWI-T2).

FIGURE 2
www.frontiersin.org

Figure 2 The accuracy of 5-fold CV by adding features sequentially. The best performance was achieved using the Top-11 feature set in the 3D-ADC-ap-DWI-T2 multiparametric model.

Clinical Usefulness

Decision curve analysis was used to assess the clinical usefulness of the radiomics models to guide identifying TP53 status (Figure 3). Three radiomics models had clinical utility due to the net benefit of models were greater than treating all and none patients (Figure 4). Figure 4 showed two cases of multiparametric MRI images from pancreatic ductal adenocarcinoma patients with wild-TP53(A-D) and TP53 mutation(E-H). MRI findings were similar, and gene mutation could not be distinguished. However, gene mutations can be accurately classified by the model. Net benefit was maximized with threshold probabilities of 0%-19% by the “3D-ADC-DWI-pp-T2” model. If the threshold probability was more than 19%, the “3D-ADC-ap-DWI-T2” model would add more net benefit compared with the other radiomics models across the majority of the range of threshold probabilities.

FIGURE 3
www.frontiersin.org

Figure 3 Decision curve analysis for 3 radiomics models. The y-axis measures the net benefit. The x-axis represented the threshold probability. The dashed line represents the assumption that all patients underwent model I, model II and model III test; the horizontal black line represents the assumption that no patients underwent MRI test; The blue line represents the 3D-ADC-ap-DWI-T2 model; the orange line represents the 3D-ADC-DWI-pp-T2 model; the green line represents the 3D-ADC-ap-pp-T2 model.

FIGURE 4
www.frontiersin.org

Figure 4 Two cases of multiparametric MRI images from pancreatic ductal adenocarcinoma patients with wild-TP53 (A–D) from a 78 year-old women and TP53 mutation (E–H) from a 53 year-old men. ADC map showed hypointense lesion (A, E). Slightly hypovascular lesion on late arterial phase (B, F). DWI sequence showing hyperintense lesion (C, G). T2WI showed slight hyperintensity of the pancreatic head mass (D, H).

Discussion

The results of our study suggested that radiomics models using different MRI-based multisequence TA has the potential for identifying the TP53 mutation status in PDAC.

We extracted several quantitative 2D and 3D radiomics features from seven different MRI sequences and constructed 378 SVM classifiers, including 2D,3D and the combined 2D and 3D modalities. Subsequently, three models with better predictive performance were selected to predict TP53 status in pancreatic cancer. The candidate features were reduced to 11 potential variables using the DX score, based on the accuracy of the SVM classifiers. SVM and DX scores were used at the same time, achieving a good classification accuracy of the TP53 mutation status in PDAC patients. All features extracted from 3D in consist with the previous study,3D analysis may provide a more representative evaluation of tumor heterogeneity (23). Although no independent validation cohort was available, the model was validated using five-fold cross-validation, which is a common procedure for model validation with a limited number of samples. Here, the combined model proved valid in the cross-validation yielding high diagnostic accuracies above 91.0%.

Radiomics have been previously correlated with a specific genotype or molecular phenotype in different cancer types such as lung (24), brain (25), and rectal cancer (26). However, radiomics had previously seen limited use in PDAC characterization. Marc A. Attiyeh et al. (27) found an association between resectable PDAC imaging features and SMAD4 status using CT texture analysis, but the model did not predict TP53 status. Si Shi et al. (28) found a correlation of PET-imaging features with TP53 status in terms of metabolic tumour burden. Yosuke et al. had reported a model for predicting TP53 mutations in pancreatic cancer from CT images using machine learning (29), and its AUC value was 0.795. But CT had no high contrast resolution that can reflect indistinguishable lesions in PDAC. Our results found that 3D texture features from T2WI, ADC, DWI and DCE T1WI in ap were influential in the analysis to find a classifier for TP53 status characterization.

SGT-53 is a gene therapy anti-cancer therapeutic agent comprised of a cationic liposome encapsulating a plasmid encoding wild-type p53. A phase I trial showed that SGT-53 is well tolerated, exhibits anticancer activity and reaches metastatic lesions in patients with different solid tumors (30). And a Phase II clinical trial of SGT-53 plus gemcitabine and nab-paclitaxel (NCT02340117) was used for metastatic pancreatic cancer. In the future, TP53 status in PDAC may play a greater role in treatment selection.

There are several limitations in our study. First, this was a retrospective study from a single center and population was small. More data from multiple institutions were needed to validate our results and acquired the optimized model. Moreover, the effect of field strength (1.5 T and 3 T) on radiomics had been unclear due to the limitation of the retrospective data. Thus, we used normalization of the values extracted from images to improve the repeatability of features. Second, the method of including patients who underwent surgery would also lead to selection bias. Finally, we used only one software for the texture analysis and one method. Therefore, the applicability of our model to other software and algorithms is uncertain.

Conclusion

In conclusion, the radiomics model derived from mp-MRI provided a non-invasive, quantitative method to predict TP53 mutation status in PDAC. Therefore, this radiomics model may help clinicians to select optimal therapies in patients with PDAC.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Board of Ruijin Hospital affiliated to Shanghai Jiao Tong University School of Medicine. The ethics committee waived the requirement of written informed consent for participation.

Author Contributions

Study concept and design: XzL and XhQ. Acquisition and segmentation of image: JG, XdL, and WhF. Analysis and interpretation of data: XhC. Clinical studies: BL and FM. Draft of the initial manuscript: JG and XhC. Further review and editing of the manuscript: XzL and XhQ. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by research grants from the National Natural Science Foundation of China (62073218), Shanghai Science and Technology Committee of Shanghai Municipality (20Y11912300), and Shanghai Municipal Key Clinical Specialty (shslczdzk03403).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.632130/full#supplementary-material

References

1. Rawla P, Sunkara T, Gaduputi V. Epidemiology of Pancreatic Cancer: Global Trends, Etiology and Risk Factors. World J Oncol (2019) 10:10–27. doi: 10.14740/wjon1166

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Adamska A, Elaskalani O, Emmanouilidi A, Kim M, Abdol Razak NB, Metharom P, et al. Molecular and Cellular Mechanisms of Chemoresistance in Pancreatic Cancer. Adv Biol Regul (2018) 68:77–87. doi: 10.1016/j.jbior.2017.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Hwang RF, Gordon EM, Anderson WF, Parekh D. Gene Therapy for Primary and Metastatic Pancreatic Cancer With Intraperitoneal Retroviral Vector Bearing the Wild-Type p53 Gene. Surgery (1998) 124(2):143–50. doi: 10.1016/S0039-6060(98)70114-X

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Morton JP, Timpson P, Karim SA, Ridgway RA, Athineos D, Doyle B, et al. Mutant p53 Drives Metastasis and Overcomes Growth Arrest/Senescence in Pancreatic Cancer. Proc Natl Acad Sci USA (2010) 107:246–51. doi: 10.1073/pnas.0908428107

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Temraz S, Shamseddine A, Mukherji D, Charafeddine M, Tfayli A, Assi H, et al. Ki67 and P53 in Relation to Disease Progression in Metastatic Pancreatic Cancer: A Single Institution Analysis. Pathol Oncol Res (2019) 25:1059–66. doi: 10.1007/s12253-018-0464-y

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Fiorini C, Cordani M, Padroni C, Blandino G, Di Agostino S, Donadelli M. Mutant p53 Stimulates Chemoresistance of Pancreatic Adenocarcinoma Cells to Gemcitabine. Biochim Biophys Acta (2015) 1853:89–100. doi: 10.1016/j.bbamcr.2014.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Weissmueller S, Manchado E, Saborowski M, Morris JPT, Wagenblast E, Davis CA, et al. Mutant p53 Drives Pancreatic Cancer Metastasis Through Cell-Autonomous PDGF Receptor Beta Signaling. Cell (2014) 157:382–94. doi: 10.1016/j.cell.2014.01.066

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kurtanich T, Roos N, Wang G, Yang J, Wang A, Chung EJ. Pancreatic Cancer Gene Therapy Delivered by Nanoparticles. SLAS Technol (2019) 24:151–60. doi: 10.1177/2472630318811108

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zhang X, Tian Q, Wang L, Liu Y, Li B, Liang Z, et al. Radiomics Strategy for Molecular Subtype Stratification of Lower-Grade Glioma: Detecting IDH and TP53 Mutations Based on Multimodal MRI. J Magn Reson Imaging (2018) 48:916–26. doi: 10.1002/jmri.25960

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chicklore S, Goh V, Siddique M, Roy A, Marsden PK, Cook GJ. Quantifying Tumour Heterogeneity in 18F-FDG PET/CT Imaging by Texture Analysis. Eur J Nucl Med Mol Imaging (2013) 40:133–40. doi: 10.1007/s00259-012-2247-0

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Li K, Xiao J, Yang J, Li M, Xiong X, Nian Y, et al. Association of Radiomic Imaging Features and Gene Expression Profile as Prognostic Factors in Pancreatic Ductal Adenocarcinoma. Am J Trans Res (2019) 11:4491–9.

Google Scholar

12. Attiyeh MA, Chakraborty J, Doussot A, Langdon-Embry L, Mainarich S, Gonen M, et al. Survival Prediction in Pancreatic Ductal Adenocarcinoma by Quantitative Computed Tomography Image Analysis. Ann Surg Oncol (2018) 25:1034–42. doi: 10.1245/s10434-017-6323-3

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Moon WK, Chen HH, Shin SU, Han W, Chang RF. Evaluation of TP53/PIK3CA Mutations Using Texture and Morphology Analysis on Breast MRI. Magn Reson Imaging (2019) 63:60–9. doi: 10.1016/j.mri.2019.08.026

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Dang M, Lysack JT, Wu T, Matthews TW, Chandarana SP, Brockton NT, et al. MRI Texture Analysis Predicts p53 Status in Head and Neck Squamous Cell Carcinoma. AJNR Am J Neuroradiol (2015) 36:166–70. doi: 10.3174/ajnr.A4110

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Chu LC, Park S, Kawamoto S, Fouladi DF, Shayesteh S, Zinreich ES, et al. Utility of CT Radiomics Features in Differentiation of Pancreatic Ductal Adenocarcinoma From Normal Pancreatic Tissue. AJR Am J Roentgenol (2019) 213:349–57. doi: 10.2214/ajr.18.20901

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-Guided 3D Active Contour Segmentation of Anatomical Structures: Significantly Improved Efficiency and Reliability. Neuroimage (2006) 31:1116–28. doi: 10.1016/j.neuroimage.2006.01.015

PubMed Abstract | CrossRef Full Text | Google Scholar

17. van Griethuysen JJM, Fedorov A, Parmar C, Hosny A, Aucoin N, Narayan V, et al. Computational Radiomics System to Decode the Radiographic Phenotype. Cancer Res (2017) 77:e104–7. doi: 10.1158/0008-5472.Can-17-0339

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Tan H, Bao J, Zhou X. A Novel Missense-Mutation-Related Feature Extraction Scheme for ‘Driver’ Mutation Identification. Bioinformatics (2012) 28:2948–55. doi: 10.1093/bioinformatics/bts558

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Qian X, Tan H, Zhang J, Zhuang X, Branch L, Sanger C, et al. Objective Classification System for Sagittal Craniosynostosis Based on Suture Segmentation. Med Phys (2015) 42:5545–58. doi: 10.1118/1.4928708

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Chang CC, Lin CJ. LIBSVM: A Library for Support Vector Machines. ACM Trans Intell Syst Technol (2011) 2(3):1–27. doi: 10.1145/1961189.1961199

CrossRef Full Text | Google Scholar

21. Vickers AJ, Elkin EB. Decision Curve Analysis: A Novel Method for Evaluating Prediction Models. Med Decis Making (2006) 26:565–74. doi: 10.1177/0272989x06295361

PubMed Abstract | CrossRef Full Text | Google Scholar

22. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the Areas Under Two or More Correlated Receiver Operating Characteristic Curves: A Nonparametric Approach. Biometrics (1988) 44:837–45. doi: 10.2307/2531595

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ng F, Kozarski R, Ganeshan B, Goh V. Assessment of Tumor Heterogeneity by CT Texture Analysis: can the Largest Cross-Sectional Area be Used as an Alternative to Whole Tumor Analysis? Eur J Radiol (2013) 82:342–8. doi: 10.1016/j.ejrad.2012.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Liu Y, Kim J, Balagurunathan Y, Li Q, Garcia AL, Stringfield O, et al. Radiomic Features are Associated With EGFR Mutation Status in Lung Adenocarcinomas. Clin Lung Cancer (2016) 17:441–8.e6. doi: 10.1016/j.cllc.2016.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Li Y, Qian Z, Xu K, Wang K, Fan X, Li S, et al. MRI Features Predict p53 Status in Lower-Grade Gliomas Via a Machine-Learning Approach. NeuroImage Clin (2018) 17:306–11. doi: 10.1016/j.nicl.2017.10.030

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Oh JE, Kim MJ, Lee J, Hur BY, Kim B, Kim DY, et al. Magnetic Resonance-Based Texture Analysis Differentiating KRAS Mutation Status in Rectal Cancer. Cancer Res Treat (2020) 52:51–9. doi: 10.4143/crt.2019.050

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Attiyeh MA, Chakraborty J, McIntyre CA, Kappagantula R, Chou Y, Askan G, et al. CT Radiomics Associations With Genotype and Stromal Content in Pancreatic Ductal Adenocarcinoma. Abdom Radiol (NY) (2019) 44:3148–57. doi: 10.1007/s00261-019-02112-1

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Shi S, Ji S, Qin Y, Xu J, Zhang B, Xu W, et al. Metabolic Tumor Burden is Associated With Major Oncogenomic Alterations and Serum Tumor Markers in Patients With Resected Pancreatic Cancer. Cancer Lett (2015) 360:227–33. doi: 10.1016/j.canlet.2015.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Iwatate Y, Hoshino I, Yokota H, Ishige F, Itami M, Mori Y, et al. Radiogenomics for Predicting p53 Status, PD-L1 Expression, and Prognosis With Machine Learning in Pancreatic Cancer. Br J Cancer (2020) 123:1253–61. doi: 10.1038/s41416-020-0997-1

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Senzer N, Nemunaitis J, Nemunaitis D, Bedell C, Edelman G, Barve M, et al. Phase I Study of a Systemically Delivered p53 Nanoparticle in Advanced Solid Tumors. Mol Ther (2013) 21:1096–103. doi: 10.1038/mt.2013.32

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: pancreatic ductal adenocarcinoma, TP53, radiomics, support vector machine, multiparametric MRI

Citation: Gao J, Chen X, Li X, Miao F, Fang W, Li B, Qian X and Lin X (2021) Differentiating TP53 Mutation Status in Pancreatic Ductal Adenocarcinoma Using Multiparametric MRI-Derived Radiomics. Front. Oncol. 11:632130. doi: 10.3389/fonc.2021.632130

Received: 22 November 2020; Accepted: 27 April 2021;
Published: 17 May 2021.

Edited by:

Bo Gao, Affiliated Hospital of Guizhou Medical University, China

Reviewed by:

Zhibing Ruan, Affiliated Hospital of Guizhou Medical University, China
Xiaohu Li, First Affiliated Hospital of Anhui Medical University, China

Copyright © 2021 Gao, Chen, Li, Miao, Fang, Li, Qian and Lin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xiaohua Qian, xiaohua.qian@sjtu.edu.cn; Xiaozhu Lin, lxz11357@rjh.com.cn

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.