- 1Department of Radiology, Sichuan Provincial People’s Hospital, University of Electronic Science and Technology of China, Chengdu, China
- 2Department of Radiology, Chinese Academy of Sciences Sichuan Translational Medicine Research Hospital, Chengdu, China
- 3Department of Radiology, Lanzhou University Second Hospital, Lanzhou, China
- 4Department of Radiology, Fifth Affiliated Hospital of Zunyi Medical University, Zhuhai, China
- 5Department of Radiology, Affiliated Hospital of Qinghai University, Xining, China
- 6Department of Pharmaceuticals Diagnosis, GE Healthcare, Beijing, China
- 7Clinical Science Department, Philips (China) Investment Co., Ltd., Chengdu, China
- 8Department of Radiology and BME, University of California Irvine, Irvine, CA, United States
Background: This study aimed to noninvasively predict the mutation status of epidermal growth factor receptor (EGFR) molecular subtype in lung adenocarcinoma based on CT radiomics features.
Methods: In total, 728 patients with lung adenocarcinoma were included, and divided into three groups according to EGFR mutation subtypes. 1727 radiomics features were extracted from the three-dimensional images of each patient. Wilcoxon test, least absolute shrinkage and selection operator regression, and multiple logistic regression were used for feature selection. ROC curve was used to evaluate the predictive performance of the model. Nomogram was constructed by combining radiomics features and clinical risk factors. Calibration curve was used to evaluate the goodness of fit of the model. Decision curve analysis was used to evaluate the clinical applicability of the model.
Results: There were three, two, and one clinical factor and fourteen, thirteen, and four radiomics features, respectively, which were significantly related to each EGFR molecular subtype. Compared with the clinical and radiomics models, the combined model had the highest predictive performance in predicting EGFR molecular subtypes [Del-19 mutation vs. wild-type, AUC=0.838 (95% CI, 0.799-0.877); L858R mutation vs. wild-type, AUC=0.855 (95% CI, 0.817-0.894); and Del-19 mutation vs. L858R mutation, AUC=0.906 (95% CI, 0.869-0.943), respectively], and it has a stable performance in the validation set [AUC was 0.813 (95% CI, 0.740-0.886), 0.852 (95% CI, 0.790-0.913), and 0.875 (95% CI, 0.781-0.929), respectively].
Conclusion: Our combined model showed good performance in predicting EGFR molecular subtypes in patients with lung adenocarcinoma. This model can be applied to patients with lung adenocarcinoma.
Introduction
Targeted therapy has brought based on recognizing the importance of acquired gene driver mutations, such as epidermal growth factor receptor (EGFR) mutations, kristen rat sarcoma (KRAS) mutations and anaplastic lymphoma kinase (ALK) rearrangements, in non-small cell lung cancer (NSCLC) new hope to patients with these gene mutations. In the Asian population, about 50% of lung adenocarcinoma patients have known carcinogenic driver genes (1, 2). There are currently targeted drugs used in clinical practice for these mutations, such as gefitinib and osimertinib for EGFR mutations. In contrast, patients without these mutations are not candidates for targeted therapy (3). Furthermore, there are molecular differences between each molecular mutation and molecular subtype, and these differences lead to different therapeutic effects after using other targeted drugs (4). EGFR mutations mainly include exon 18-21 mutations. Among them, exon 19 deletion (Del-19) mutation and 21 L858R point (L858R) mutation are the two most common activating mutations, and they are also the two most sensitive mutation sites for tyrosine kinase inhibitors (TKI) treatment (5). In a single targeted therapy, patients with Del-19 mutation benefited more from osimertinib (6), while patients with L858R mutation benefited significantly from dacomitinib (7); in addition, combination therapy and immunotherapy brought patients with L858R for more potential benefits (8, 9). Therefore, the detection of specific EGFR mutation subtypes can make targeted therapies more precise and allow patients receiving these treatments to benefit the most.
Currently, the detection of EGFR mutation status from histological specimens is the most common detection method. However, in clinical practice, these detection techniques also have some limitations. For example, tissue samples are obtained through invasive methods such as biopsy or surgery; sometimes the amount of tissue samples obtained due to operational errors is insufficient; biopsy can increase the risk of tumor metastasis; in addition, a small part of the tissue obtained does not represent the heterogeneity of the entire tumor, etc. (10–12). In addition, another noninvasive detection strategy for EGFR mutations is ‘liquid biopsy’, which is a biological detection method on the blood. For patients with advanced NSCLC, ‘liquid biopsy’ is a promising method to isolate circulating tumor DNA from blood samples (13). However, ‘liquid biopsy’ has a high risk of false-negative results (30%) (14). Therefore, until this defect is effectively resolved, ‘liquid biopsy’ is far from substitute for histological testing. Because of this, there is an urgent need for a simple and noninvasive method to detect EGFR mutation subtypes before targeted drug therapy.
The radiological features have been shown to reflect EGFR mutation status in lung adenocarcinoma (12, 15–17). However, the clinical applicability of these studies needs to be confirmed by further research. Compared with traditional CT, radiomics converts medical images into mineable data and extracts a large number of features that cannot be observed by the human naked eye system, thereby reflecting more characteristics of tumors (18). To our knowledge, some studies have used radiomics to predict EGFR mutation status (19–22). Although the prediction performance of these studies is different, this shows that it is feasible to predict EGFR mutations noninvasively through radiomics. However, only a few studies have used radiomics methods to predict the mutation status of EGFR molecular subtypes (23–26). Unfortunately, the sample size included in these studies is limited, and the accuracy of the obtained prediction model was only 65.5-79.0%. In this study, we retrospectively collected a relatively large data set and constructed a model based on CT radiomics signature to noninvasively predict the mutation status of EGFR molecular subtype in lung adenocarcinoma.
Materials and Methods
Patient Population
This retrospective study was ethically approved by the Institutional Review Board of the Sichuan Provincial People’s Hospital and Lanzhou University Second Hospital, and the need for patient informed consent was waived. Clinical data and chest CT images of these patients were obtained from the picture archiving and communication system (PACS). The inclusion criteria were as follows: (1) patients with the histologic type of lung adenocarcinoma; (2) patients with complete CT thin-slice images (1.25 mm) and clinical data; (3) patients who did not receive lung cancer-related treatment before CT scan; (4) patients who underwent biopsy or surgery within one month after CT scan; (5) patients with EGFR exon Del-19 mutation, exon L858R mutation, and wild-type. The exclusion criteria were as follows: (1) patients whose tumor boundary is difficult to be recognized by the naked eye on CT images; (2) patients younger than 18 years old.
According to the above inclusion and exclusion criteria, 728 patients (median age, 57.0 years, age range, 21-82 years, 370 males and 358 females) were finally selected from 2,557 patients in the two hospitals. Among them, a total of 540 patients from Sichuan Provincial People’s Hospital were used as the training set from January 2018 to March 2021, and 188 patients from Lanzhou University Second Hospital were identified as the external validation set from January 2019 to September 2020. The patient recruitment flowchart is shown in Figure 1.
Demographic and clinical data include the patient’s sex, age, smoking history [including non-smoking (never smoked) and smoking (former or current smoking)], carcinoembryonic antigen (CEA) level, and tumor lobe location of the tumor (including right upper, right middle, right lower, left upper and left lower lobes). If the tumor crosses the fissure, the lobe location is defined as the lobe in which the tumor predominates.
EGFR Mutation Status Detection
The polymerase chain reaction-amplified refractory mutation system (PCR-ARMS) detected EGFR mutation status. The human EGFR gene detection kit (Beijing SinoMD Gene Detection Technology Co., Ltd., China; Amoy Diagnostics, Xiamen, China) detected EGFR exon 18 to 21 mutation status.
CT Image Acquisition
CT scans ranged from the thoracic inlet to the level of the lower edge of the 12th rib were completed by three spiral CT scanners (Discovery CT750 HD, GE Healthcare; Philips iCT 256, Koninklijke Philips N.V.; Somatom Sensation 64, Siemens Healthineers). Scanning parameters were as follows: (1) tube voltage 120 kVp, tube current adjusted automatically for the Sensation 64 scanner, and (2) tube voltage 120 kVp, tube current 150 to 200 mA for the other two scanners. For all scanners, 0.5-1.0 second tube rotation time, and field of view (FOV): 350 mm; matrix, 512 × 512; the layer thickness and spacing were both 5 mm; the reconstruction layer thickness and spacing were both 1.25 mm. All images were exported in DICOM format to facilitate feature extraction.
Tumor Segmentation and Radiomics Feature Extraction
Radiomics feature extraction and analysis workflow are shown in Figure 2. To ensure the accuracy and consistency of the data, two readers (radiologists with 6 and 4 years of experience in chest CT diagnosis, respectively) independently used the open-source software ITK-SNAP 3.8.0 (http://www.itksnap.org) to segment the tumor on the thin-slice CT lung window (window width: 1500HU; window level: -500HU). Training cases were segmented by reader 1 (G.J.Z), and validation cases were segmented by reader 2 (L.N.D). Both readers were blinded to all patients’ clinical data, pathological records, and EGFR status results. When the two readers were unsure, a consultant radiologist (J.Z) confirmed the segmentation with 17 years of experience. The region of interest (ROI) was manually segmented on CT axial images with tumor tissue and confirmed on sagittal and coronal images.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/889293/fonc-12-889293-HTML/image_m/fonc-12-889293-g002.jpg)
Figure 2 Flowchart of the process of radiomics. (A) The tumours were segmented on CT images to form the region of interest (ROI). (B) Radiomics feature extraction from the ROI. (C) Radiomics feature dimensionality reduction process. (D) Construct a radiomics model.
In order to evaluate the robustness and repeatability of the radiomics feature extraction process, one month later, 40 patients were randomly selected from the training set and segmented again by readers 1 and 2 to construct a re-segmentation set, and 40 patients were randomly selected from each CT scanner to construct different CT scanner sets for calculating the intra-/interclass correlation coefficients (ICC), respectively. ICC values > 0.8 reflected good consistency (26).
The open-source Python software package PyRadiomics 3.0.1 automatically extracted radiomics features from the three-dimensional (3D) tumor volume after segmentation. Radiomics features were divided into three main categories: 16 shape features, 324 first-order features and 1387 texture features. Details of radiomics features were included in the Supplementary Material (Methods). According to the recommendations of International Symposium on Biomedical Imaging (ISBI), we only resampled the image and set the bin width of gray discretization to 25. We performed z-score preprocessing on the extracted radiomics features.
Radiomics Feature Selection
To avoid model overfitting and improve accuracy, we used three steps for feature selection to obtain the optimal feature subset. First, Wilcoxon test was used to retain the features with P-value less than 0.05. Secondly, the least absolute shrinkage and selection operator (LASSO) regression with 10-fold cross-validation was used to eliminate the collinearity features. LASSO is a recognized algorithm that has been used for feature selection of higher-dimensional variables (27). Finally, multiple logistic regression was used to select the features, and the minimum features of akaike information criterion (AIC) were retained.
For clinical factors, the Chi-square and Student’s t-tests were first used in the training set to screen for clinical characteristics that were correlated between each group. P-value was set to 0.05. Clinical factors with a P-value greater than 0.05 were excluded. Next, use logistic regression to further analyze the most relevant variables.
Radiomics Model Establishment
Logistic regression was used in the training set to build a model for predicting Del-19 or L858R mutations, and its performance was evaluated in the external validation set. Logistic regression is a classic method in radiomics research. It is easy to understand, explain, and combine discrete and continuous variables (28, 29). To identify Del 19 and L858R mutations, we used logistic regression, support vector machine (SVM), and random forest (RF) to construct prediction models. The clinical and radiomics models were constructed based on clinical factors and radiomics features, respectively, while the combined model was constructed based on clinical and radiomics models. Additionally, clinical models were constructed using logistic regression.
Statistical analysis
All statistical analyses were performed using R 3.6.0 (http://www.R-project.org). Two-sided P-values less than 0.05 were considered to be statistically different. Kolmogorov-Smirnov test was used to evaluate the normal distribution of the data. Categorical variables were expressed as percentiles, and the Chi-square test or Fisher’s exact test was used to analyze groups. Continuous variables were expressed as mean ± standard deviation (SD), and Student’s t-test or Mann-Whitney U test was used for analysis between groups. Receiver operating characteristic (ROC) curve was used to evaluate the performance of the model, and the area under the curve (AUC), sensitivity, specificity, accuracy, positive predictive value (PPV), and negative predictive value (NPV) were calculated. Delong test was used to compare the performance differences of the prediction models.
Based on the above-screened clinical factors and radiomics features, we constructed a personalized nomogram to predict the mutation status of the EGFR molecular subtype. Calibration curve and Hosmer-Lemeshow (H-L) test were used to evaluate the model’s goodness of fit. Decision curve analysis (DCA) was used to assess the clinical applicability of the model.
Results
Clinical Characteristics of Patients
There were no significant differences in clinical factors (including sex, smoking history, and CEA), the mutation rate of Del-19 or L858R, and tumor location in each EGFR mutant subtype group (all P > 0.05), while there were significant differences in age between the two EGFR mutant subtype groups (Del-19 vs. wild type, Del-19 vs. L858R) (Supplementary Tables S1–S3).
Univariate analysis revealed that age, sex and smoking history were significantly different between the Del-19 mutation and wild-type groups (P < 0.05), sex and smoking history were significantly different between the L858R mutation and wild-type groups (P < 0.05), and age was significantly different between the Del-19 mutation and L858R mutation groups (P < 0.05). Multivariate analysis revealed that age (OR, 0.972; 95% CI, 0.948-0.996; P = 0.021) and sex (OR, 3.193; 95% CI, 1.836-5.565; P < 0.001) were correlated independently with the task of Del-19 vs. wild-type (Table 1), sex (OR, 2.612; 95% CI, 1.548-4.457; P < 0.001) and smoking history (OR, 0.427; 95% CI, 0.238-0.761; P = 0.001) were correlated independently with the task of L858R vs. wild-type (Table 2), and age (OR, 1.050; 95% CI, 1.022-1.081; P < 0.001) was correlated independently with the task of Del-19 vs. L858R (Table 3). Based on multivariate analysis results, clinical factors with P < 0.05 in each task were incorporated in the clinical model.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/889293/fonc-12-889293-HTML/image_m/fonc-12-889293-t001.jpg)
Table 1 The relationship between clinical variables of patients and EGFR molecular subtypes (Del-19 mutation vs. Wild-type) in the training set.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/889293/fonc-12-889293-HTML/image_m/fonc-12-889293-t002.jpg)
Table 2 The relationship between clinical variables of patients and EGFR molecular subtypes (L858R mutation vs. Wild-type) in the training set.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/889293/fonc-12-889293-HTML/image_m/fonc-12-889293-t003.jpg)
Table 3 The relationship between clinical variables of patients and EGFR molecular subtypes (Del-19 mutation vs. L858R mutation) in the training set.
Radiomics Feature Selection and Model Establishment
In total, 1727 radiomics features were extracted from the 3D images of each ROI. The ICC values of the radiomics features extracted from two readers and different CT scanners were all greater than 0.80, reflecting good consistency. Fourteen radiomics features were highly correlated with Del-19 mutation (Table S4; Figure S1), thirteen radiomics features were highly correlated with L858R mutation (Table S4; Figure S2). For Del-19 mutation vs. L858R mutation, only four radiomics features were screened after using the Wilcoxon test (Table S4; Figure S3). Therefore, we retained these four features to construct the prediction model.
Correlation analysis showed that the correlation between each feature is weak and independent in the training and validation sets (Figures S4–S6).
Based on the above-screened radiomics features and clinical factors, the clinical, radiomics, and combined models were established in the training set, respectively, to predict the EGFR molecular subtype mutation status.
Predictive performance and Validation Based on Clinical, Radiomics, and Combined Models
The predictive performance of different models in the training and validation sets is shown in Figure 3 and Table 4. The predictive performance of the combined model was higher than that of other single models. In the training set, the AUC of the combined model was 0.838 (95% CI, 0.799-0.877), 0.855 (95% CI, 0.817-0.894), and 0.906 (95% CI, 0.869-0.943), respectively. In addition, we used an external validation set to verify the accuracy of the combined model, and the AUC was 0.813 (95% CI, 0.740-0.886), 0.852 (95% CI, 0.790-0.913), and 0.875 (95% CI, 0.781-0.929), respectively. In addition, when distinguishing between Del-19 and L858R mutations, the prediction model’s performance constructed using random forest was higher than that of other single models. The AUC of the training and validation sets were 0.881 (95% CI, 0.840-0.921) and 0.871 (95% CI, 0.802-0.941), respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/889293/fonc-12-889293-HTML/image_m/fonc-12-889293-g003.jpg)
Figure 3 Receiver operating characteristic (ROC) curves of the three models were used to predict the mutant status of EGFR molecular subtypes. (A, B) Del-19 mutation vs. wild-type. (C, D) L858R mutation vs. wild-type. (E, F) Del-19 mutation vs. L858R mutation. (A, C, E) Training set. (B, D, F) Validation set.
Delong test showed that there were significant differences in AUC values of the three models in the training set between EGFR Del-19 mutation or L858R mutation and wild-type groups (all P < 0.05); However, only the AUC value of combined model and clinical model was significantly different in the validation set (P < 0.05), and the AUC values between other models were not statistically significant (P > 0.05) (Figures S7A–D). There were significant differences in AUC values of the combined model and clinical, SVM or Logistic models in the training and validation sets between EGFR Del-19 mutation and L858R mutation groups (all P < 0.05). However, the AUC value between combined model and RF model was not statistically significant in the both sets (P > 0.05) (Figures S7E, F).
Clinical Application of the Combined Model
Based on radiomics score and clinical risk factors, we constructed two user-friendly nomograms to predict the mutation status of EGFR molecular subtypes (Figures 4A and 5A). The detailed formula for calculating the radiomics score is shown in the Supplementary Material (Result). The calibration curve analysis showed that the probability of Del-19 mutation or L858R mutation predicted by the combined model was highly consistent with the actual possibility, indicating that the model had the best discriminant ability (Figures 4B, C, and 5B, C). Decision curve analysis showed that the combined model threshold in range of 0.18-0.77 have higher net benefit for Del-19 vs. wild type and the cutoff value was 0.440 fall in this rang; the combined model threshold in range of 0.16-0.715 have higher net benefit for L858R vs. wild type and the cutoff value was 0.389 fall in this range (Figures 4D, E and 5D, E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/889293/fonc-12-889293-HTML/image_m/fonc-12-889293-g004.jpg)
Figure 4 Nomogram was used to identify Del-19 mutation and wild-type. (A) Construct a nomogram in the training set based on the combined model. (B, C) Calibration curve of the combined model in the training (B) and validation (C) sets. The x-axis represents the use of the combined model to predict the risk of Del-19 mutation. The y-axis represents the actual Del-19 mutation rate. The green, red, and blue lines represent the distinguishing ability of the clinical, radiomics, and combined models, respectively, while the gray diagonal line represents the ideal evaluation of the ideal model. The closer the fit to the diagonal line indicates the better discrimination ability. (D, E) Decision curve analysis for the combined model in the training (D) and validation (E) sets. The x-axis shows the threshold probability, and the y-axis measures the net benefit. The gray line represents all patients with Del-19 mutation, and the black line represents all patients without Del-19 mutation. The green, red, and blue lines represent the clinical, radiomics, and combined models, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/889293/fonc-12-889293-HTML/image_m/fonc-12-889293-g005.jpg)
Figure 5 Nomogram was used to identify L858R mutation and wild-type. (A) Construct a nomogram in the training set based on the combined model. (B, C) Calibration curve of the combined model in the training (B) and validation (C) sets. (D, E) Decision curve analysis for the combined model in the training (D) and validation (E) sets.
The precision-recall curves showed that the combined model constructed by the RF model combined with clinical factors had better performance than other single models in predicting Del-19 and L858R mutations (Figures 6A, B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/889293/fonc-12-889293-HTML/image_m/fonc-12-889293-g006.jpg)
Figure 6 Precision-recall (PR) curves of the different models in the training (A) and validation sets (B). PR represents the relationship between precision and recall.The larger the area under the PR curve, the better the model performance.
Discussion
Preoperative noninvasive prediction of EGFR mutant subtypes is a new field that attracts researchers’ attention. It can well overcome some shortcomings of molecular mutation detection based on histology and provide critical information for the rational formulation of targeted therapy in clinical practice. This study established different models based on clinical factors and radiomics features to predict EGFR mutation subtypes. Among them, the combined model showed good predictive performance in the training set. It also had good stability when evaluating the model’s performance in the external validation set, which reveals the feasibility of predicting EGFR molecular subtypes through radiomics features.
Previous studies have found that some clinical variables such as female, non-smokers, patients with histological type of lung adenocarcinoma, and East Asian populations are significantly associated with EGFR mutations (16, 17, 20, 30, 31). Our previous research has also confirmed this (12, 15). However, these studies did not further analyze the correlation between EGFR mutation subtypes and clinical variables. In our study, sex, smoking history, and age were significantly different in the EGFR mutation subtypes group. Compared with EGFR wild-type patients, Del-19 mutation patients were more common in females and non-smokers, and L858R mutation patients were more common in females and non-smokers. Compared with patients with Del-19 mutation, patients with L858R mutation were relatively older. Only patients with Del-19 mutation and L858R mutation were selected because they are the most common mutations in EGFR mutation subtypes. The two mutation sites most related to the sensitivity of EGFR TKI treatment.
Some previous studies have predicted the mutation status of EGFR molecular subtypes based on CT radiomics features and achieved promising results. For example, Li and colleagues (26) retrospectively collected 312 patients with NSCLC, and 580 radiomics features were extracted from each patient’s CT images to construct a model to predict EGFR mutation subtypes (Del 19 and L858R). The test set’s AUC for predicting Del 19 and L858R mutations were 79.3% and 77.5%, respectively. Similarly, Zhao et al. (25) included 637 patients with lung adenocarcinoma in their study to predict EGFR mutation subtypes, and extracted 475 radiomics features to construct a model. The results showed that the AUC in the training and validation datasets were 68.9% and 75.7%, respectively. However, these studies did not distinguish between the Del-19 mutation and the L858R mutation, limiting the clinical applicability of these studies. In this study, we not only distinguished between EGFR Del-19 mutation or L858R mutation and EGFR wild-type. More importantly, we further distinguished the Del-19 mutation and the L858R mutation, and achieved good prediction performance. The training and validation sets’ AUC was 90.6% and 87.5%, respectively. Therefore, our research may be more in line with actual clinical needs.
In this study, whether in the training or validation sets, the combined model established by clinical factors combined with radiomics features can improve the diagnostic performance of identifying EGFR molecular subtypes. Liu et al. (24) included 263 patients with lung adenocarcinoma in their study to detect EGFR mutation status and its molecular subtypes. Among the 6 models established, the combined model had better distinguishing ability than the model that only uses radiomics features or clinical factors. Tu et al. (20) included 404 NSCLC patients in their study to predict EGFR mutation status, and the comprehensive model showed higher predictive performance than any other single model. Jia et al. (32) also showed that a comprehensive model with radiomics features combined with clinical factors had better diagnostic performance than a single model. It shows that adding clinical factors to the radiomics model can improve the diagnostic performance of the model.
In recent years, the study of radiomics in predicting tumor gene mutations has attracted extensive attention from researchers (20, 22, 32, 33). The intrinsic relationship between the radiomics features and EGFR mutation status in patients with lung adenocarcinoma can be further explored through data mining to guide clinical decision-making, predict prognosis and evaluate efficacy (19, 22, 32). This study investigated the relationship between radiomics features and EGFR molecular subtypes. Among these features, most of them were texture features, indicating that texture features were more closely related to EGFR molecular subtypes. The human visual system cannot recognize these features, nor can they be interpreted as specific meanings (34, 35). We observed that radiomics features, including logarithm_glcm_Correlation, wavelet.LLL_glszm_Zone Entropy (ZE), and gradient_glszm_Gray Level Non-Uniformity Normalized (GLNN), etc. were associated significantly with Del-19 mutation. Among them, logarithm_glcm_Correlation and wavelet.LLL_glszm_ZE reflected the image texture of the tumor area, and gradient_glszm_GLNN reflected the gray image value of the tumor area (36). Compared with the EGFR wild-type group, the values of these features were higher in Del-19 mutation, indicating that the image texture and gray image values were related to Del-19 mutation. Lbp.3D.m1_firstorder_10Percentile and Lbp.3D.m1_firstorder_Skewness, etc. were associated significantly with L858R mutation. They reflected the voxel intensity of the image (36). Compared with the EGFR wild-type group, the values of these features were higher in L858R mutation, indicating that the image voxel intensity was related to L858R mutation. Therefore, radiomics features as a new auxiliary tool can predict EGFR molecular subtypes.
Compared with the radiomics model based on only containing radiomics features, incorporating preoperative clinical factors of the nomogram showed the best predictive performance. This user-friendly nomogram will help clinicians easily predict EGFR molecular subtypes in clinical practice. The results were more practical than a single model and can be used for clinical applications in patients with lung adenocarcinoma undergoing CT scans. The task of Del-19 vs. wild-type and L858R vs. wild-type build with linear model (logistic regression) could obtain a satisfactory result, and the linear model is easy for application. Such we didn’t applied nonlinear model. The task of Del-19 vs. L858R was hard, the performance of linear model was not satisfactory, so we add nonlinear model for comparison and selected best model for radiomics score construction.
Our study had several limitations. First, although we collected data from two large medical centers, this was a retrospective study and there may be inevitable selection bias. The conclusions of this study need to be prospectively verified in more centers to improve the clinical applicability of our model. Second, although our study included 728 patients, increasing the sample size will further improve the accuracy of the results of this study. Finally, two radiologists spent a lot of time manually segmenting ROI. Therefore, ROI can be segmented automatically and effectively in future research.
Conclusion
In conclusion, demonstrated the feasibility of identifying EGFR molecular subtypes through the radiomics features of patients with lung adenocarcinoma, making the formulation of clinically individualized targeted therapy programs more precise and more in line with actual clinical needs, so as to benefit the patients with candidate targeted therapy the most.
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 the medical ethics committees of Sichuan Provincial People’s Hospital and Lanzhou University Second Hospital. The ethics committee waived the requirement to participate in written informed consent.
Author Contributions
GZ, LD, WK, JZ and HP contributed to conception and design of the study. GZ, LD, SL, RQ, SP and XZ organized the database. GZ, YC, JZ and JR performed the statistical analysis. GZ, LD, JZ and YC wrote the first draft of the manuscript. YC, JZ and ZZ wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This study received funding from the Sichuan Provincial Cadre Health Research Project (No. Chuan Gan Yan 2022-208), Qinghai Province "Kunlun Talents High-end Innovation and Entrepreneurial Talents" Top Talent Cultivation Project, Medical Science and Technology Research Fund Project of Guangdong Province (B2022144), Science and Technology Plan Fund of Guizhou Provincial (Qiankehe Foundation-ZK [2022] General 634), and Doctoral research start-up fund project of Zunyi Medical University (BS2021-03).
Conflict of Interest
Author JR was employed by GE Healthcare. Author XZ was employed by Philips (China) Investment Co., Ltd.
The remaining 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.
Publisher’s Note
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.889293/full#supplementary-material
References
1. Harrison PT, Vyse S, Huang PH. Rare Epidermal Growth Factor Receptor (EGFR) Mutations in non-Small Cell Lung Cancer. Semin Cancer Biol (2020) 61:167–79. doi: 10.1016/j.semcancer.2019.09.015
2. D’Angelo SP, Pietanza MC, Johnson ML, Riely GJ, Miller VA, Sima CS, et al. Incidence of EGFR Exon 19 Deletions and L858R in Tumor Specimens From Men and Cigarette Smokers With Lung Adenocarcinomas. J Clin Oncol (2011) 29(15):2066–70. doi: 10.1200/jco.2010.32.6181
3. Li WQ, Cui JW. Non-Small Cell Lung Cancer Patients With Ex19del or Exon 21 L858R Mutation: Distinct Mechanisms, Different Efficacies to Treatments. J Cancer Res Clin Oncol (2020) 146(9):2329–38. doi: 10.1007/s00432-020-03296-6
4. Jänne PA, Yang JC, Kim DW, Planchard D, Ohe Y, Ramalingam SS, et al. AZD9291 in EGFR Inhibitor-Resistant non-Small-Cell Lung Cancer. N Engl J Med (2015) 372(18):1689–99. doi: 10.1056/NEJMoa1411817
5. Lee DH. Treatments for EGFR-Mutant non-Small Cell Lung Cancer (NSCLC): The Road to a Success, Paved With Failures. Pharmacol Ther (2017) 174:1–21. doi: 10.1016/j.pharmthera.2017.02.001
6. Ramalingam SS, Vansteenkiste J, Planchard D, Cho BC, Gray JE, Ohe Y, et al. Overall Survival With Osimertinib in Untreated, EGFR-Mutated Advanced NSCLC. N Engl J Med (2020) 382(1):41–50. doi: 10.1056/NEJMoa1913662
7. Mok TS, Cheng Y, Zhou X, Lee KH, Nakagawa K, Niho S, et al. Improvement in Overall Survival in a Randomized Study That Compared Dacomitinib With Gefitinib in Patients With Advanced Non-Small-Cell Lung Cancer and EGFR-Activating Mutations. J Clin Oncol (2018) 36(22):2244–50. doi: 10.1200/jco.2018.78.7994
8. Mok TS, Wu YL, Ahn MJ, Garassino MC, Kim HR, Ramalingam SS, et al. Osimertinib or Platinum-Pemetrexed in EGFR T790M-Positive Lung Cancer. N Engl J Med (2017) 376(7):629–40. doi: 10.1056/NEJMoa1612674
9. Soria JC, Ohe Y, Vansteenkiste J, Reungwetwattana T, Chewaskulyong B, Lee KH, et al. Osimertinib in Untreated EGFR-Mutated Advanced Non-Small-Cell Lung Cancer. N Engl J Med (2018) 378(2):113–25. doi: 10.1056/NEJMoa1713137
10. Yang X, Dong X, Wang J, Li W, Gu Z, Gao D, et al. Computed Tomography-Based Radiomics Signature: A Potential Indicator of Epidermal Growth Factor Receptor Mutation in Pulmonary Adenocarcinoma Appearing as a Subsolid Nodule. Oncologist (2019) 24(11):e1156–e64. doi: 10.1634/theoncologist.2018-0706
11. Zhang G, Cao Y, Zhang J, Zhao Z, Zhang W, Zhou J. Epidermal Growth Factor Receptor Mutations in Lung Adenocarcinoma: Associations Between Dual-Energy Spectral CT Measurements and Histologic Results. J Cancer Res Clin Oncol (2021) 147(4):1169–78. doi: 10.1007/s00432-020-03402-8
12. Zhang G, Zhang J, Cao Y, Zhao Z, Li S, Deng L, et al. Nomogram Based on Preoperative CT Imaging Predicts the EGFR Mutation Status in Lung Adenocarcinoma. Transl Oncol (2021) 14(1):100954. doi: 10.1016/j.tranon.2020.100954
13. Rolfo C, Mack PC, Scagliotti GV, Baas P, Barlesi F, Bivona TG, et al. Liquid Biopsy for Advanced Non-Small Cell Lung Cancer (NSCLC): A Statement Paper From the IASLC. J Thorac Oncol (2018) 13(9):1248–68. doi: 10.1016/j.jtho.2018.05.030
14. Goldman JW, Noor ZS, Remon J, Besse B, Rosenfeld N. Are Liquid Biopsies a Surrogate for Tissue EGFR Testing? Ann Oncol (2018) 29(suppl_1):i38–46. doi: 10.1093/annonc/mdx706
15. Zhang G, Zhao Z, Cao Y, Zhang J, Li S, Deng L, et al. Relationship Between Epidermal Growth Factor Receptor Mutations and CT Features in Patients With Lung Adenocarcinoma. Clin Radiol (2021) 76(6):473.e17–24. doi: 10.1016/j.crad.2021.02.012
16. Liu Y, Kim J, Qu F, Liu S, Wang H, Balagurunathan Y, et al. CT Features Associated With Epidermal Growth Factor Receptor Mutation Status in Patients With Lung Adenocarcinoma. Radiology (2016) 280(1):271–80. doi: 10.1148/radiol.2016151455
17. Lee HJ, Kim YT, Kang CH, Zhao B, Tan Y, Schwartz LH, et al. Epidermal Growth Factor Receptor Mutation in Lung Adenocarcinomas: Relationship With CT Characteristics and Histologic Subtypes. Radiology (2013) 268(1):254–64. doi: 10.1148/radiol.13112553
18. Gillies RJ, Kinahan PE, Hricak H. Radiomics: Images Are More Than Pictures, They Are Data. Radiology (2016) 278(2):563–77. doi: 10.1148/radiol.2015151169
19. Zhang G, Cao Y, Zhang J, Ren J, Zhao Z, Zhang X, et al. Predicting EGFR Mutation Status in Lung Adenocarcinoma: Development and Validation of a Computed Tomography-Based Radiomics Signature. Am J Cancer Res (2021) 11(2):546–60.
20. Tu W, Sun G, Fan L, Wang Y, Xia Y, Guan Y, et al. Radiomics Signature: A Potential and Incremental Predictor for EGFR Mutation Status in NSCLC Patients, Comparison With CT Morphology. Lung Cancer (2019) 132:28–35. doi: 10.1016/j.lungcan.2019.03.025
21. Hong D, Xu K, Zhang L, Wan X, Guo Y. Radiomics Signature as a Predictive Factor for EGFR Mutations in Advanced Lung Adenocarcinoma. Front Oncol (2020) 10:28. doi: 10.3389/fonc.2020.00028
22. 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(5):441–48.e6. doi: 10.1016/j.cllc.2016.02.001
23. Mei D, Luo Y, Wang Y, Gong J. CT Texture Analysis of Lung Adenocarcinoma: Can Radiomic Features be Surrogate Biomarkers for EGFR Mutation Statuses. Cancer Imaging (2018) 18(1):52. doi: 10.1186/s40644-018-0184-2
24. Liu G, Xu Z, Ge Y, Jiang B, Groen H, Vliegenthart R, et al. 3D Radiomics Predicts EGFR Mutation, Exon-19 Deletion and Exon-21 L858R Mutation in Lung Adenocarcinoma. Transl Lung Cancer Res (2020) 9(4):1212–24. doi: 10.21037/tlcr-20-122
25. Zhao W, Wu Y, Yn X, Sun Y, Gao P, Tan M, et al. The Potential of Radiomics Nomogram in Non-Invasively Prediction of Epidermal Growth Factor Receptor Mutation Status and Subtypes in Lung Adenocarcinoma. Front Oncol (2020) 9:1485(1485). doi: 10.3389/fonc.2019.01485
26. Li S, Ding C, Zhang H, Song J, Wu L. Radiomics for the Prediction of EGFR Mutation Subtypes in Non-Small Cell Lung Cancer. Med Phys (2019) 46(10):4545–52. doi: 10.1002/mp.13747
27. Zhang L, Chen B, Liu X, Song J, Fang M, Hu C, et al. Quantitative Biomarkers for Prediction of Epidermal Growth Factor Receptor Mutation in Non-Small Cell Lung Cancer. Transl Oncol (2018) 11(1):94–101. doi: 10.1016/j.tranon.2017.10.012
28. Ji GW, Zhang YD, Zhang H, Zhu FP, Wang K, Xia YX, et al. Biliary Tract Cancer at CT: A Radiomics-Based Model to Predict Lymph Node Metastasis and Survival Outcomes. Radiology (2019) 290(1):90–8. doi: 10.1148/radiol.2018181408
29. Huang Y, Liu Z, He L, Chen X, Pan D, Ma Z, et al. Radiomics Signature: A Potential Biomarker for the Prediction of Disease-Free Survival in Early-Stage (I or II) Non-Small Cell Lung Cancer. Radiology (2016) 281(3):947–57. doi: 10.1148/radiol.2016152234
30. Zhou JY, Zheng J, Yu ZF, Xiao WB, Zhao J, Sun K, et al. Comparative Analysis of Clinicoradiologic Characteristics of Lung Adenocarcinomas With ALK Rearrangements or EGFR Mutations. Eur Radiol (2015) 25(5):1257–66. doi: 10.1007/s00330-014-3516-z
31. Dearden S, Stevens J, Wu YL, Blowers D. Mutation Incidence and Coincidence in Non Small-Cell Lung Cancer: Meta-Analyses by Ethnicity and Histology (Mutmap). Ann Oncol (2013) 24(9):2371–6. doi: 10.1093/annonc/mdt205
32. Jia TY, Xiong JF, Li XY, Yu W, Xu ZY, Cai XW, et al. Identifying EGFR Mutations in Lung Adenocarcinoma by Noninvasive Imaging Using Radiomics Features and Random Forest Modeling. Eur Radiol (2019) 29(9):4742–50. doi: 10.1007/s00330-019-06024-y
33. Park JE, Park SY, Kim HJ, Kim HS. Reproducibility and Generalizability in Radiomics Modeling: Possible Strategies in Radiologic and Statistical Perspectives. Korean J Radiol (2019) 20(7):1124–37. doi: 10.3348/kjr.2018.0070
34. Wu Q, Yao K, Liu Z, Li L, Zhao X, Wang S, et al. Radiomics Analysis of Placenta on T2WI Facilitates Prediction of Postpartum Haemorrhage: A Multicentre Study. EBioMedicine (2019) 50:355–65. doi: 10.1016/j.ebiom.2019.11.010
35. Zhang J, Yao K, Liu P, Liu Z, Han T, Zhao Z, et al. A Radiomics Model for Preoperative Prediction of Brain Invasion in Meningioma non-Invasively Based on MRI: A Multicentre Study. EBioMedicine (2020) 58:102933. doi: 10.1016/j.ebiom.2020.102933
Keywords: NSCLC, lung adenocarcinoma, EGFR, computed tomography, radiomics
Citation: Zhang G, Deng L, Zhang J, Cao Y, Li S, Ren J, Qian R, Peng S, Zhang X, Zhou J, Zhang Z, Kong W and Pu H (2022) Development of a Nomogram Based on 3D CT Radiomics Signature to Predict the Mutation Status of EGFR Molecular Subtypes in Lung Adenocarcinoma: A Multicenter Study. Front. Oncol. 12:889293. doi: 10.3389/fonc.2022.889293
Received: 04 March 2022; Accepted: 04 April 2022;
Published: 29 April 2022.
Edited by:
Yiyan Liu, University of Louisville, United StatesReviewed by:
Wenjuan Ma, Tianjin Medical University Cancer Institute and Hospital, ChinaCheng Dong, The Affiliated Hospital of Qingdao University, China
Copyright © 2022 Zhang, Deng, Zhang, Cao, Li, Ren, Qian, Peng, Zhang, Zhou, Zhang, Kong and Pu. 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: Guojin Zhang, zhanggj15@lzu.edu.cn; Hong Pu, pqs19951211@126.com; Weifang Kong, 13936034@qq.com
†These authors have contributed equally to this work