Skip to main content
Download PDF

Noninvasive KRAS mutation estimation in colorectal cancer using a deep learning method based on CT imaging

BMC Medical Imaging volume 20, Article number: 59 (2020) Cite this article

Abstract

Background

The detection of Kirsten rat sarcoma viral oncogene homolog (KRAS) gene mutations in colorectal cancer (CRC) is key to the optimal design of individualized therapeutic strategies. The noninvasive prediction of the KRAS status in CRC is challenging. Deep learning (DL) in medical imaging has shown its high performance in diagnosis, classification, and prediction in recent years. In this paper, we investigated predictive performance by using a DL method with a residual neural network (ResNet) to estimate the KRAS mutation status in CRC patients based on pre-treatment contrast-enhanced CT imaging.

Methods

We have collected a dataset consisting of 157 patients with pathology-confirmed CRC who were divided into a training cohort (n = 117) and a testing cohort (n = 40). We developed an ResNet model that used portal venous phase CT images to estimate KRAS mutations in the axial, coronal, and sagittal directions of the training cohort and evaluated the model in the testing cohort. Several groups of expended region of interest (ROI) patches were generated for the ResNet model, to explore whether tissues around the tumor can contribute to cancer assessment. We also explored a radiomics model with the random forest classifier (RFC) to predict KRAS mutations and compared it with the DL model.

Results

The ResNet model in the axial direction achieved the higher area under the curve (AUC) value (0.90) in the testing cohort and peaked at 0.93 with an input of ’ROI and 20-pixel’ surrounding area. AUC of radiomics model in testing cohorts were 0.818. In comparison, the ResNet model showed better predictive ability.

Conclusions

Our experiments reveal that the computerized assessment of the pre-treatment CT images of CRC patients using a DL model has the potential to precisely predict KRAS mutations. This new model has the potential to assist in noninvasive KRAS mutation estimation.

Peer Review reports

Background

CRC is the third most commonly diagnosed malignancy and the fourth leading cause of cancer-related deaths in the world, and its burden is expected to increase by 60% to more than 2.2 million new cases and 1.1 million cancer deaths by 2030 [1]. The distinct heterogeneity of prognosis and treatment response has been observed in clinical practice among various CRC patients, even in those who share similar pathological stages and treatment modalities. KRAS is a small G protein that plays a role in the epidermal growth factor receptor (EGFR) pathway. KRAS mutations are acquired early during colorectal tumorigenesis. In approximately 40% of all metastatic colorectal cancer tumors, one of several heterozygous KRAS codon 12 or 13 mutations is detected [2]. Targeted treatments with anti-EGFR monoclonal antibodies (cetuximab and panitumumab) are recommended by the National Comprehensive Cancer Network (NCCN) for metastatic CRC patients whose tumors do not harbor KRAS mutations [3, 4]. Therefore, the detection of KRAS gene mutations in CRC is key to the optimal design of individualized therapeutic strategies.

Invasive colonoscopy or surgery for the biopsy is the limitation of KRAS mutational testing. The samples from these invasive procedures may be limited by intratumoral heterogeneity and may not sufficiently represent the exact macroscopic status of the entire tumor. Moreover, the NCCN [3] recommends that fresh biopsies should not be obtained solely for the purpose of KRAS genotyping. Histological biopsies for KRAS testing are usually performed before the initial treatment and cannot be repeated in the subsequent treatment. Thus, it is very useful to explore a noninvasive method to assist in the diagnosis of KRAS mutations. The noninvasive prediction of the KRAS status in CRC is challenging. Medical imaging technology provides some exploratory methods for the detection of the KRAS status. Several previous studies have used positron emission tomography (PET)/computed tomography (CT) with 18-F fluorodeoxyglucose PET/CT to assess KRAS mutation information in CRC [5–7]. An alternative CT-based radiomics approach for noninvasive KRAS mutation estimation in CRC has also been applied [8].

Artficial intelligence (AI) has become a hot topic in medical care support. Broadly speaking, machine learning can be divided into two major classes: radiomic analysis, which relies on multi-step pipelines, and the concept of the DL method, which simplifies this pipeline by learning predictive features on its own. Recent studies have shown the potential for DL in medical imaging due to its high performance in diagnosis, classification, and prediction [9–11]. Previous research has also shown the great potential of DL in the prediction of key molecular markers in gliomas [12, 13].A convolutional neural network (CNN) is a type of neural network developed specifically to learn hierarchical representations of imaging data. A ResNet is an effective exploration of a deep CNN, allowing the effective training of substantially deeper networks than those used previously while maintaining fast convergence times [14]. A ResNet has been used increasingly due to its utility and simplicity in clinical applications [15, 16].

To the best of our knowledge, there is no research on whether a CT-based DL model can predict the KRAS mutation status in CRC. Therefore, we developed and evaluated a ResNet model on contrast-enhanced CT before treatment to noninvasively predict the KRAS status in CRC. We also explored a radiomics model to compare with the DL model.

Methods

Data collection

This research was approved by the First Hospital of Jilin University medical ethics committee. All the patients involved in the study signed an informed consent. A total of 157 patients who had pathologically confirmed CRC and underwent KRAS mutation tests and contrast-enhanced CT pretreatment were identified retrospectively between December 2016 and January 2018. Inclusion criteria: (a) Archive data for patients with pathology-confirmed colorectal cancer from December 2016 to January 2018; (b) Patients before treatment who underwent a KRAS mutation test; (c) Pre-treatment contrast-enhanced CT available; and (d) Contrast-enhanced CT with a reconstruction slice thickness of 1.5 mm. Exclusion criteria: (a) CRC patients who underwent radiotherapy, chemotherapy or chemoradiotherapy before obtaining the pathological tissue sample; (b) Patients without enhanced CT before specimen collection. These patients were classified according to CT acquisition time and then assigned to either training cohort or testing cohort on a 3:1 ratio. There were 75 men and 42 women (mean age, 60 years; age range, 28 −88 years) assigned to the training cohort and 31 men and 9 women (mean age, 59 years; age range, 32 −73 years) assigned to the testing cohort. KRAS mutations (exons 2, 3, and 4) were analysed by a next-generation sequencing (NGS) method.

CT technique

All patients underwent contrast-enhanced CT examinations using one of two 64-detector row spiral CT scanners (Philips, Brilliance iCT). CT scan was performed after 65s delay following intravenous injection of 100ml Iopromide (Uitravist-300; Bayer Schering Pharma, Berlin Germany) at a rate of 3ml/s for enhancement. The scanning parameters were as follows: 120 kV; 150 mAs; rotation time, 0.5 s; and matrix size, 512 ×512. The section thickness was 1.5 mm, and the interval was 3 mm.

Data normalization

The portal venous phase images obtained from contrast-enhanced CT in the axial, coronal, and sagittal directions were first preprocessed by performing intensity normalization to reduce the noise and inconsistencies due to low-frequency non-uniformity or the inhomogeneity of intensities[17].The Z-score[18]method was applied to normalize the training and testing cohorts:

$$I'=\frac{I-\bar{I}}{\sigma_{I}} $$

where I indicates a slice from any cohort, \(\bar {I}\) indicates the average greyscale value of the two cohorts and σ1 indicates the greyscale standard deviation of the two cohorts.

ROI patch generation

Following normalization, the ROIs were manually delineated along the contour of tumor on the largest tumor cross section in axial direction, coronal direction and sagittal direction. Air area inside the tumor area was excluded from the contour. Two professional radiologists with 12 (reader 1, W.X.) and 7 (reader 2, X.X.Z.) years of experience delineated the ROI in a blinded fashion. The two-dimensional segmentation was completed using ITK-SNAP software (version 3.4; www.itk-snap.org). An independent samples t-test was used to evaluate the differences between the features generated by reader 1 and reader 2 (interobserver), and the differences between the twice-generated features by reader 1 (intraobserver). Inter- and intraclass correlation coefficients (ICCs) were used to assess the agreement of feature extraction. It was considered acceptable that the ICC is greater than 0.75.

To extract original patches, a minimum bounding rectangle was first drawn around each manual ROI. That ensured the entire tumor area, as well as the minimum peritumoral tissue, was captured. The region of the minimum bounding rectangle was cropped in this slice. Since the size of all the ROI patches was different, we counted the median height and median width from all the ROI patches for three groups. Then, nonlinear interpolation was used to resize the ROI patches for three groups. The ROI patch sizes after resizing were 60 ×60 pixels for the axial direction, 68 ×63 pixels for the sagittal direction and 72 ×67 pixels for the coronal direction.

Several groups of expended ROI patches were generated for the ResNet model, to explore whether tissues around the tumor can contribute to cancer assessment. The minimum bounding rectangle was enlarged in both height and width, with an interval of 10 pixels. We cropped and resized the patches in the same manner as in the original group ROI patches in three directions. For axial images, we generated 3 groups of patches as the subsequent experimental dataset. The sizes were 70 ×70 pixels for A1-set, 80 ×80 pixels for A2-set and 90 ×90 pixels for A3-set. Similarly, we generated 2 groups of patches for sagittal images (78 ×73-pixel patches of S1-set and 88 ×83-pixel patches of S2-set) and coronal images (82 ×87-pixel patches of C1-set and 92 ×97-pixel patches of C2-set). The CT images of all patients in the training cohort were constructed using the above steps. The detailed workflow of ROI patch generation and relative expansion is shown in Fig. 1.

Fig. 1
figure 1

The flow chart of ROI-patch generation. A minimum circumscribed rectangle is established around each irregular ROI. Then, the boundary of the minimum rectangle is equally expanded with an interval of 10 pixels. This procedure is repeated for the three orthogonal directions (axial, sagittal, and coronal)

Full size image

Data augmentation

DL networks are often data-driven architectures, and data augmentation is an effective way to improve the model and reduce overfitting. Therefore, we augmented the training dataset by introducing random rotations and translations, generating ’new’ training data. The augmentation technique allows us to further increase the size of our training cohort. All images were converted into 5 pixels, 10 pixels, 15 pixels, and 20 pixels and rotated by 3 ∘, 6 ∘, 9 ∘, 12 ∘, and 15 ∘. Then, the above 9 transformations and original images were separately rotated by 1 ∘, 2 ∘, 3 ∘, 4 ∘, and 5 ∘. Therefore, the number of images in the primary cohort expanded to 50 times its initial scale.

ResNet

A ResNet was applied to train the imaging data and build our neural network model. There are six residual learning blocks. As shown in Fig. 2, for input x, the residual block output y is defined as:

$$y=\sigma(F(x,W_{i})+x) $$

where the function is the i-th residual mapping and σ denotes the rectified linear unit (ReLU) [19] process. The kernel size of all the convolution layers is 5 ×5. Then,a max-pooling layer, a fully connected layer and a soft-max layer are implemented.

Fig. 2
figure 2

The structure of the employed residual neural network. There are six identity blocks, a pooling layer, a fully-connected layer and a softmax. Each identity block has three convolutional layers. The kernel size of all the convolution layers is 5 ×5. ReLu are adopted after every convolutional layer

Full size image

In the training phase, we trained the ResNet model with the constructed datasets in the portal venous phase images of the axial, coronal, and sagittal directions. The datasets in each direction contained both the original patch and the extended patch. We fed the ResNet with patches of different sizes and obtained 10 corresponding pre-trained models, which were models with original axial patches, A1-set, A2-set, and A3-set, original sagittal patches, S1-set and S2-set, and original coronal patches, C1-set and C2-set. In the testing phase, we evaluated the performances of the above 10 pre-trained models, respectively.

All experiments were performed on the workstation of a Windows 10 64-bit operating system with a 64-GB memory and an NVIDIA GeForce GTX 1080 GPU. Data normalization and ROI generation were performed in MATLAB 2016b. Data augmentation, training and testing for all the ResNet models were developed on the Keras library with a TensorFlow backend. When training the ResNet, the Adam optimization function was used with a batch size of 40 and a learning rate of 0.001.

Radiomics model

We also explored a radiomics model with RFC to predict KRAS mutations and compared it with the DL model. Random forest classifier has been a prevalent data mining and statistical tool because of its transparency and great success in classification and regression task [20, 21]. A total of 1025 features, including tumor intensity, shape and size, texture, and wavelet characteristics, were extracted from the primary tumors based on the manually delineated ROI. Detailed descriptions of these features are shown in Supplementary Information 4.1. Feature selection and modelling were based on the training cohort. A univariate analysis was performed for each feature. Features with P values <0.05 were considered associated with KRAS mutations and were incorporated into the least absolute shrinkage and selection operator (LASSO) logistic regression model with 10-fold cross-validation. We established a radiomics model with an RFC according to the low-dimensional radiomics feature signature. The RFC consists of multiple classification and regression trees (CARTs), which are highly accurate and tolerant to exception values and noise without being prone to overfitting. Detailed descriptions of the radiomics method are shown in Additional file 1.

Results

Patient demographics

The demographic and tumor characteristics in the training and testing cohorts are listed in Table 1. Based on the results of KRAS status, the patients were classified into two groups: the mutated group and the wild-type group. Regarding gender, age and TNM stage, no demographic differences were observed between the two groups in both cohorts (p >0.05). However, for tumor size, statistical significance was found in both cohorts (P values <0.05). There were significant differences in tumor location between the two groups in the testing cohort, but these differences were not confirmed in the training cohort.

Table 1 Demographic differences in the training and testing cohorts No, number; m, median; SD standard deviation
Full size table

Predictive performance of the ResNet classification model

We investigated the predictive model for KRAS gene mutation for CRC patients based on CT images in three positions: axial direction, coronal direction and sagittal direction. Metrics of the AUC, sensitivity and specificity were used to evaluate the performances of the networks. In addition, several groups with different input sizes were carried out to explore the effect of the surrounding tissue on classification accuracy. All the results for the testing cohort are displayed in Table 2. Figure 3a shows the AUC in the axial direction peaked at 0.93 when the input was ’ROI and 20-pixel’, and the increase in the AUC value was limited by the expansion of the input size. Figure 3b and c show that the inclusion of surrounding tissues did not contribute to the ResNet model either in the coronal or sagittal direction. Figure 4 shows that the ResNet model in the axial direction reached the higher AUC value compared with the coronal and sagittal positions.

Fig. 3
figure 3

ROC curves for four KRAS mutations predicted by the residual neural network and radiomics models in testing cohort. a ResNet and raidiomics predictions on different input in axial direction. b ResNet predictions on different input in coronal direction. c ResNet predictions on different input in sagittal direction

Full size image
Fig. 4
figure 4

The line chart of AUC values for each CNN with different inputs in testing cohort. ResNet model in the axial direction reached the higher AUC value compared with the coronal and sagittal

Full size image
Table 2 Performance of Models of the KRAS mutation prediction in testing cohorts
Full size table

Predictive performance of the radiomics model

For the feature selection of the high-dimensional radiomic features of the training cohort, 6-dimensional features of P <0.05 are obtained by t-test. Then, the 4-dimensional radiomics feature signatures are selected by LASSO regularization with 10-fold cross-validation. The feature names, values and coefficients of the radiomics feature signatures are shown in Additional file 1. We use RFC as the basic model. The model parameters are n_estimators = 200, max_depth = 4; that is, the number of CART is 200, and the maximum tree depth is 4. In the training cohort, the AUC was 0.945 (sensitivity: 0.75; specificity: 0.94), and in the testing cohort, the AUC was0.818 (sensitivity: 0.70; specificity: 0.85).In comparison, the ResNet model showed better predictive ability. The ROC curves are shown in Fig. 3a.

Discussion

Targeted treatments with anti-EGFR monoclonal antibodies (cetuximab and panitumumab) are recommended by the NCCN for metastatic CRC patients whose tumors do not harbour KRAS mutations [3, 4].The determination of the KRAS mutational status in CRC patients is essential for the management of CRC. However, several issues with KRAS testing limit its utility, as described previously. The aim of our study was to employ a ResNet model to provide a noninvasive preoperative prediction of the KRAS mutational status. To the best of our knowledge, this is the first study concerning the predictive value of DL based on CT images in CRC. Previous studies have attempted to investigate the relationship between image characteristics and genetic mutations [6, 22, 23] or pathological subtype [24, 25]. Compared with other medical imaging technology, one of the challenges of DL is the requirement of a large training dataset. A limited number of suitable patients might lead to insufficient training data. Therefore, in the study, data augmentation was used to increase the size of the training dataset and prevent overfitting. Besides, the ResNet with simple architecture and short time consumption was enough to learn the predictive features according to the satisfactory results in our study and the other researchers [26, 27].

To explore the lesion-based prediction in three positions: axial, coronal, and sagittal, and the influence of surrounding tissues on the classified results, we cropped several groups of patches around the manually annotated ROI. When the input was ‘ROI only’, the ResNet model in the axial direction reached the higher AUC value compared with the coronal and sagittal direction. That shows the reconstructed images based on 2D axial images do not contribute to the improvement of model performance. One interesting finding is that the inclusion of surrounding tissues contributed to the ResNet just in the axial direction, and the continuous expansion of the input size limited the improvement of AUC.The AUC in the axial direction peaked at 0.93 when the input was ‘ROI and 20-pixel’, which indicated that the surrounding tissue could contribute to the model’s performance. However, when further expand the surrounding edge information, more noise may be introduced, so the classification ability of the network decreased. The inclusion of surrounding tissues for KRAS mutational classification in CRC which has biological importance should be explored further.

Our study also compared the performances of the DL method and the CT-based radiomics method. Some studies found that slice thickness affected radiomic feature values and the performance of diagnostic models [28, 29]. To reduce the effects of reconstruction slice thickness on the performance of the radiomic model, we only included images of the same thickness. However, the diagnostic performance of the ResNet model with 2D input data was unaffected by slice thickness. Thus, the proposed deep learning model is suitable for routine CT imaging with other slice thickness. In our study, the proposed CT-based radiomics signature incorporated four radiomics features: two high-dimensional features obtained by the LoG and two high-dimensional features obtained by the wavelet transform, which are all conducive to classification, shown in Additional file 1. We chose the RFC as the machine learning classifier. The RFC consists of multiple CARTs, and each CART trains the sub-classifier by bootstrapping. Finally, the RFC predicts results through the voting of each CART, so it has good generalization [30]. In the training cohort, the AUC was 0.945 (sensitivity: 0.75; specificity: 0.94), and in the testing cohort, the AUC was 0.818 (sensitivity: 0.70; specificity: 0.85). These results are consistent with the results of other investigators [8]. The clinical background and tumor stage were not used in the radiomics model or the ResNet model to ensure the consistency of the input data type. The comparison of the two models indicates that the CT-based DL model could reach the predictive level of the radiomics model in predicting the KRAS mutation status of CRC patients, while the DL can simplify the multi-step pipeline of the conventional radiomics method with little pre-processing and relatively greater reproducibility [31].

There are several limitations to this preliminary study. First, the number of study subjects was limited, and although data augmentation was performed to increase the size of the training cohort, it was still difficult to observe a robust outcome. However, the deep neural network has achieved good performance in many clinical medical applications, both with large and relatively small sample sizes [10, 32, 33]. Second, the input data type was 2D and not 3D in the DL model. 3D deep learning requires significantly higher computation power than sequential 2D image analyses. A further study with the differences between 2D and 3D models of DL in predicting KRAS mutations in CRC will need to be undertaken. Third, our study analyzed the KRAS gene status of only the RAS-RAF-MAPK pathway. The V-raf murine sarcoma viral oncogene homolog B (BRAF) mutation is another marker of anti-EGFR resistance in non-first-line treatment of metastatic CRC. We did not accumulate enough BRAF mutation cases because of its low prevalence in CRC. Fourth, we included only a single-center cohort with an internal testing set. In the future, large, multi-center cohorts should be recruited for evaluation.

Conclusions

The aim of the present research was to develop and evaluate a DL model on contrast-enhanced CT before treatment to noninvasively predict the KRAS status in CRC. We also explored a radiomics model with RFC to compare with the DL model (ResNet). The notable findings of our study can be summarized as follows:

(1) Computerized assessment of the pre-treatment CT images of CRC patients using DL has the potential to precisely predict KRAS mutations. (2) The 2D ResNet model in the axial direction reached the higher AUC value compared with the coronal and sagittal direction. (3) The continuous expansion of the input size with surrounding tissues limited the performance of the ResNet model. (4) CT-based DP model may reach the predictive level of the radiomics model in predicting the KRAS mutation status of CRC patients, while the DL can simplify the multi-step pipeline of the conventional radiomics method.However, there are still several limitations in the DL research field, including the insufficient number of patients, and further studies should be performed to optimize and verify its utility.

Availability of data and materials

The datasets generated and/or analysed during the study are not publicly available due to patient privacy protection, but are available from the corresponding author on reasonable request.

Abbreviations

RFC:

Random forest classifier

AUC:

Area under the curve

EGFR:

Epidermal growth factor receptor

NCCN:

National comprehensive cancer network

PET:

Positron emission tomography

CT:

Computed tomography

AI:

Artificial intelligence

CNN:

Convolutional neural network

LASSO:

Least absolute shrinkage and selection operator

CARTs:

Classification and regression trees

BRAF:

V-raf murine sarcoma viral oncogene homolog B

References

  1. Ferlay J, Soerjomataram I, Dikshit R. Cancer incidence and mortality worldwide: sources, methods and major patterns in globocan 2012. Int J Cancer. 2015; 136(5):359–86.

    Article  Google Scholar 

  2. Roth AD, Tejpar S, Delorenzi M. Prognostic role of kras and braf in stage ii and iii resected colon cancer: results of the translational study on the petacc-3, eortc 40993, sakk 60-00 trial. J Clin Oncol. 2010; 28(3):466–74.

    Article  CAS  PubMed  Google Scholar 

  3. Van Cutsem E, Lenz HJ, Kohne CH. Fluorouracil, leucovorin, and irinotecan plus cetuximab treatment and ras mutations in colorectal cancer. J Clin Oncol. 2015; 33(7):692–700.

    Article  CAS  PubMed  Google Scholar 

  4. Peeters M, Oliner KS, Price TJ. Analysis of kras/nras mutations in a phase iii study of panitumumab with folfri compared with folfri alone as second-line treatment for metastatic colorectal cancer. Clin Cancer Res. 2015; 21(24):5469–79.

    Article  CAS  PubMed  Google Scholar 

  5. Kawada K, Nakamoto Y, Kawada M. Relationship between 18f-uorodeoxyglucose accumulation and kras/braf mutations in colorectal cancer. Clin Cancer Res. 2012; 18(6):1696–703.

    Article  CAS  PubMed  Google Scholar 

  6. Chen SW, Chiang HC, Chen WT. Correlation between pet/ct parameters and kras expression in colorectal cancer. Clin Nucl Med. 2014; 39(8):685–9.

    Article  PubMed  Google Scholar 

  7. Lovinfosse P, Koopmansch B, Lambert F. (18)f-fdg pet/ct imaging in rectal cancer: relationship with the ras mutational status. Br J Radiol. 2016; 89(1063):20160212.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Yang L, Dong D, Fang M. Can ct-based radiomics signature predict kras/nras/braf mutations in colorectal cancer?Eur Radiol. 2018; 28(5):2058–67.

    Article  PubMed  Google Scholar 

  9. Cha KH, Hadjiiski L, Samala RK. Urinary bladder segmentation in ct urography using deep-learning convolutional neural network and level sets. Med Phys. 2016; 43(4):1882.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Esteva A, Kuprel B, Novoa RA. Corrigendum: Dermatologist-level classification of skin cancer with deep neural networks. Nature. 2017; 546(7660):686.

    Article  CAS  PubMed  Google Scholar 

  11. Oakden-Rayner L, Carneiro G, Bessen T. Precision radiology: Predicting longevity using feature engineering and deep learning methods in a radiomics framework. Sci Rep. 2017; 7(1):1684.

    Article  Google Scholar 

  12. Akkus Z, Ali I, Sedlar J. Deletion of chromosomal arms 1p/19q in low-grade gliomas from mr images using machine intelligence. J Digit Imaging. 2017; 30(4):469–76.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Korfiatis P, Kline TL, Lachance DH. Residual deep convolutional neural network predicts mgmt methylation status. J Digit Imaging. 2017; 30(5):622–8.

    Article  PubMed  PubMed Central  Google Scholar 

  14. He KM, Zhang XY, Ren SQ. Deep residual learning for image recognition. In: 2016 Ieee Conference on Computer Vision and Pattern Recognition. IEEE Conference on Computer Vision and Pattern Recognition: 2016. p. 770–8. https://doi.org/10.1109/cvpr.2016.90.

  15. Chen T, Liu S, Li Y. Developed and validated a prognostic nomogram for recurrence-free survival after complete surgical resection of local primary gastrointestinal stromal tumors based on deep learning. EBioMedicine. 2019; 39:272–9.

    Article  PubMed  Google Scholar 

  16. Chang K, Bai HX, Zhou H. Residual convolutional neural network for the determination of idh status in low- and high-grade gliomas from mr imaging. Clin Cancer Res. 2018; 24(5):1073–81.

    Article  CAS  PubMed  Google Scholar 

  17. Tustison NJ, Avants BB, Cook PA. N4itk: improved n3 bias correction. IEEE Trans Med Imaging. 2010; 29(6):1310–20.

    Article  PubMed  PubMed Central  Google Scholar 

  18. ZZhang ZH, Cheng Y, Liu NC. Comparison of the effect of mean-based method and z-score for field normalization of citations at the level of web of science subject categories. Scientometrics. 2014; 101(3):1679–93.

    Article  Google Scholar 

  19. Hara K, Saito D, Shouno H. Analysis of function of rectified linear unit used in deep learning. In: International Joint Conference on Neural Networks. IEEE International Joint Conference on Neural Networks (IJCNN): 2015.

  20. Shi T, Horvath S. Unsupervised learning with random forest predictors. J Comput Graph Stat. 2005; 15(1):118–38.

    Article  Google Scholar 

  21. Fernandez-Delgado M, Cernadas E, Barro S. Do we need hundreds of classifiers to solve real world classification problems?J Mach Learn Res. 2014; 15:3133–81.

    Google Scholar 

  22. Sun R, Limkin EJ, Vakalopoulou M. A radiomics approach to assess tumour-infiltrating cd8 cells and response to anti-pd-1 or anti-pd-l1 immunotherapy: an imaging biomarker, retrospective multicohort study. Lancet Oncol. 2018; 19(9):1180–91.

    Article  CAS  PubMed  Google Scholar 

  23. Zhou M, Leung A, Echegaray S. Non-small cell lung cancer radiogenomics map identifies relationships between molecular and imaging phenotypes with prognostic implications. Radiology. 2018; 286(1):307–15.

    Article  PubMed  Google Scholar 

  24. Zhang S, Song G, Zang Y. Non-invasive radiomics approach potentially predicts non-functioning pituitary adenomas subtypes before surgery. Eur Radiol. 2018; 28(9):3692–701.

    Article  PubMed  Google Scholar 

  25. Zhu X, Dong D, Chen Z. Radiomic signature as a diagnostic factor for histologic subtype classification of non-small cell lung cancer. Eur Radiol. 2018; 28(7):2772–8.

    Article  PubMed  Google Scholar 

  26. Chen T, Liu SQ, Li Y. Developed and validated a prognostic nomogram for recurrence-free survival after complete surgical resection of local primary gastrointestinal stromal tumors based on deep learning. Ebiomedicine. 2019; 39:272–9.

    Article  PubMed  Google Scholar 

  27. Kitamura G, Chung CY, Moore BE. 2nd.ankle fracture detection utilizing a convolutional neural network ensemble implemented with a small sample, de novo training, and multiview incorporation. J Digit Imaging. 2019; 32(4):672–7.

    Article  PubMed  PubMed Central  Google Scholar 

  28. He L, Huang Y, Ma Z. Effects of contrast-enhancement, reconstruction slice thickness and convolution kernel on the diagnostic performance of radiomics signature in solitary pulmonary nodule. Sci Rep. 2016; 6:34921.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Larue R, van Timmeren JE, de Jong EEC. Influence of gray level discretization on radiomic feature stability for different CT scanners, tube currents and slice thicknesses: a comprehensive phantom study. Acta Oncol. 2017; 56(11):1544–53.

    Article  PubMed  Google Scholar 

  30. van Griethuysen JJM, Fedorov A, Parmar C. Computational radiomics system to decode the radiographic phenotype. Cancer Res. 2017; 77(21):104–7.

    Article  Google Scholar 

  31. Lecun Y, Bottou L, Bengio Y. Gradient-based learning applied to document recognition. Proc IEEE. 1998; 86(11):2278–324.

    Article  Google Scholar 

  32. Bibault JE, Giraud P, Housset M. Author correction: Deep learning and radiomics predict complete response after neo-adjuvant chemoradiation for locally advanced rectal cancer. Sci Rep. 2018; 8(1):16914.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Cha KH, Hadjiiski L, Chan HP. Bladder cancer treatment response assessment in ct using radiomics with deep-learning. Sci Rep. 2017; 7(1):8738.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

This research is supported by the Jilin provincial health and family planning commission (Grant Number: 2017J073), Jilin provincial finance department (Grant Number: 2018SCZWSZX-026) and National natural science foundation of China (Grant Number: 81571737). Jilin provincial health and family planning commission had a role in the design of the study, collection, analysis, and interpretation of data, and in writing the manuscript. Jilin provincial finance department and National natural science of foundation of China provided funding for study design, analysis of data.

Author information

Authors and Affiliations

  1. The First Hospital of Jilin University, Department of Radiology, Changchun, China

    Kan He & Huimao Zhang

  2. Jilin University, College of Electronic Science and Engineering, State Key Laboratory on Integrated Optoelectronics, No.2699 Qianjin Street, Changchun, China

    Xiaoming Liu, Mingyang Li, Xueyan Li & Hualin Yang

Authors
  1. Kan He
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xiaoming Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Mingyang Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Xueyan Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Hualin Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Huimao Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

KH conceived the project. KH performed and participated in data collection. XML and XYL designed the experiments. ML and HY analyzed the data. Kan He, XML and HZ wrote the paper. All authors revised the manuscript and approved the final version prior to submission.

Corresponding author

Correspondence to Huimao Zhang.

Ethics declarations

Ethics approval and consent to participate

This research was approved by the First Hospital of Jilin University medical ethics committee and carried out in accordance with the guidelines and regulations of NCCN and Helsinki, and the need for informed consent was waived by the retrospective nature of the study.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1

Supplementary file.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

He, K., Liu, X., Li, M. et al. Noninvasive KRAS mutation estimation in colorectal cancer using a deep learning method based on CT imaging. BMC Med Imaging 20, 59 (2020). https://doi.org/10.1186/s12880-020-00457-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12880-020-00457-4

Keywords

BMC Medical Imaging

ISSN: 1471-2342

Contact us