Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Oncol., 07 January 2026

Sec. Cancer Imaging and Image-directed Interventions

Volume 15 - 2025 | https://doi.org/10.3389/fonc.2025.1722238

This article is part of the Research TopicThe Contrast Enhanced Ultrasound and Ultrasound Cavitation in the Diagnostic and Therapeutic Application of Solid TumorView all 5 articles

Multimodal ultrasonography for predicting epidermal growth factor receptor mutation in subpleural non-small cell lung carcinoma

  • 1Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Ultrasound, Peking University Cancer Hospital & Institute, Beijing, China
  • 2Department of Ultrasonography, Shengli Oil Field Central Hospital, Dongying, Shandong, China

Objective: Currently precise target treatment based on gene status significantly improved the outcome for patients with non-small cell lung cancer (NSCLC) and epidermal growth factor receptor (EGFR) was the most important gene. We aimed to develop a multimodal ultrasound model for predicting EGFR mutation status in patients with subpleural NSCLC, to provide important information for precise target treatment.

Methods: 75 patients with pathologically confirmed NSCLC were included in this retrospective study. Patients were divided into two groups based on EGFR mutation status: wild-type (n=57) and mutant (n=18). The clinical characteristics (C), conventional ultrasound (US) features, contrast-enhanced ultrasound (CEUS) characteristics, and time-intensity curve (TIC) parameters of the lung lesions were analyzed and compared between the two groups. Univariate and multivariate logistic regression determined independent predictors of EGFR mutations. Two predictive models were constructed: a C+ US model and a FULL model. Both were presented using nomograms. Receiver operating characteristic and calibration curves evaluated predictive performance of two models, while decision curve analysis (DCA) assessed clinical utility.

Results: Multivariate analysis identified smoking status, lesion boundaries, and air bronchogram as predictors in the C + US model. The FULL model identified lesion boundaries and air bronchogram on US, enhancement intensity of lesions and internal necrosis on CEUS and RT (rise time) from TIC as predictors. The C + US model achieved an AUC of 0.843, and the FULL model achieved 0.939. DCA confirmed substantial net clinical benefits.

Conclusion: The models developed in this study enabled patients who are unable to undergo invasive procedures to predict EGFR mutation status noninvasively. These findings provided an ultrasound-based diagnostic reference to support clinician decision-making and personalized treatment planning.

Introduction

Lung cancer persists as the predominant contributor to global cancer-related morbidity and mortality, with non-small cell lung cancer (NSCLC) constituting approximately 85% of cases (1). Recent advances in molecular biology have identified several driver mutations, including epidermal growth factor receptor (EGFR), BRAF, ROS1, MET, and ALK, which are closely related to the pathogenesis and progression of lung cancer (2). Consequently, targeted therapies have transformed the treatment landscape for individuals with NSCLC (3, 4). Notably, epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) are widely administered to individuals harboring EGFR mutations and have been shown to significantly prolong survival (5, 6). The decision to initiate targeted therapy relies on EGFR mutation status evaluation (7) which achieved through pathological examination. However, this process can be limited by the general condition of the patients, biopsy operation ability of the doctors, sample quality and intratumoral heterogeneity (811). Therefore, the development of noninvasive methods to predict EGFR mutation status is essential to complement gene sequencing.

Although CT remains essential in lung cancer imaging, it involves exposure to ionizing radiation (12). In contrast, multimodal ultrasound imaging offers high sensitivity, real-time monitor, and nonionizing advantages in evaluating peripheral lung lesions. Its utility in distinguishing benign from malignant pulmonary lesions and in identifying histological subtypes of lung cancer has been well-proven by our center and other studies (13, 14). Beside the role of characterization of lung lesions with CEUS, we further explored the perspective of CEUS in prediction of gene status which related to target therapy. In this study, we aimed to identify ultrasound (US) and contrast-enhanced ultrasound (CEUS) features associated with EGFR mutation status and to construct predictive models for assessing EGFR genotype. To our knowledge, this work was the first study to apply ultrasonographic imaging for predicting EGFR mutation status in subpleural non-small cell lung cancer.

Materials and methods

Study design and patients

From November 2012 to July 2024, totally 186 consecutive patients with confirmed peripheral lung NSCLC underwent US and CEUS examination at our department were collected retrospectively. The inclusion criteria were as follows: (i) age ≥18 years; (ii) histologically confirmed NSCLC; (iii) documented pretreatment EGFR mutation status; and (iv) availability of pretreatment US and CEUS images. Exclusion criteria were as follows: (i) an interval >1 month between the ultrasound examination and EGFR testing; (ii) the presence of concurrent primary malignancies and (iii) unsatisfactory US, CEUS and TIC images. The final study population included 75 individuals (30 females, 45 males) with an average age of 64.08±11.01 years. All 75 enrolled patients including 47 patients tested by our internal pathology laboratory and 28 cases tested by external laboratories holding equivalent accreditation were tested for EGFR mutations using Amplification refractory mutation system polymerase chain reaction (ARMS-PCR) method. The commercial ADx-EGFR Mutation Detection Kit (AmoyDx) was utilized on the ABI 7500 Real-Time PCR system (Applied Biosystems). This assay is designed to detect 29 somatic mutations in the EGFR gene, including exon 19 deletions, the T790M mutation, and mutations in exons 18, 20, and 21 (such as L858R and L861Q). The patient selection flowchart was illustrated in Figure 1.

Figure 1
Flowchart detailing patient enrollment from November 2012 to July 2024 for a study on NSCLC. Starting with 186 patients, 93 were excluded for unavailable EGFR mutation status, 12 for CEUS and EGFR testing delays, and 1 for concurrent malignancies. From 80 patients with determined EGFR mutations, 5 were excluded for unsatisfactory images. Finally, 75 patients remained: 18 EGFR positive and 57 EGFR negative.

Figure 1. Flow diagram of the patient’s selection. EGFR, Epidermal Growth Factor Receptor; US, ultrasound; CEUS, contrast enhanced ultrasound; NSCLC, non-small cell lung cancer; TIC, time-intensity curve.

Ethical approval

We confirmed that this study adherences to the Declaration of Helsinki. This study was approved by our local ethics committee, and written informed consent was waved because of retrospective study.

CEUS examination methods

The contrast agent was SonoVue (Bracco SpA, Italy). The lyophilized powder was reconstituted with 5 mL of saline to produce a uniform microbubble suspension containing 8 µL/mL of sulfur hexafluoride encapsulated in a phospholipid shell. A 2.4 mL bolus of the suspension was injected into the antecubital vein over 2–3 s. CEUS images were acquired using the Logiq E9 ultrasound system. The frequency and depth settings used for CEUS were identical to those of conventional US. The mechanical index was adjusted to 0.11–0.13. The timer was activated simultaneously with contrast injection. Lesions were assessed using contrast-enhanced harmonic US with continuous recording for ≥ 180 seconds. Dynamic images were saved to the internal hard drive of the ultrasound system.

Enhancement process were classified as arterial and parenchymal phases, beginning at contrast administration. The arterial phase was defined as the interval from the initial appearance of the contrast agent in the lesion to the point of peak enhancement. Time-intensity curves were generated using the built-in postprocessing software of the Logiq E9 ultrasound equipment.

Evaluation of US, CEUS, and TIC

Two radiologists (W.Y. and J.B.), each with 5–10 years of experience in lung ultrasound and contrast-enhanced ultrasound, independently evaluated the US and CEUS images. They jointly reviewed the imaging data and formulated consensus reports. Both radiologists were blinded to the patients’ EGFR mutation status.

In the US assessment, lesion boundaries were classified as well defined or poor defined. Lesion shapes were categorized as suborbicular, irregular, or wedge-shaped. Using surrounding atelectasis normal lung tissue as a reference, echogenicity was defined as hypoechoic or isoechoic. An air bronchogram was defined as a tubular or branching air-filled structure within the tumor, resulting from ultrasound reflections of residual air-filled bronchi with posterior comet-tail artifacts. Pleural retraction was described as inward pleural displacement toward the tumor, attributed to fibrotic contraction within tumor tissue.

In CEUS evaluation, enhancement patterns during the arterial phase were classified as: (i) centripetal perfusion, where contrast agent entered from the periphery toward the center; (ii) eccentric perfusion, where perfusion began on one side of the lesion and extended across; and (iii) disorganized enhancement, where perfusion occurred randomly throughout the lesion. Necrosis types were defined as: (i) type A, no necrosis; (ii) type B, single or multiple necrotic foci. Hydrothorax was recorded as present or absent. Homogeneity was classified as present or absent based on uniform echo intensity within the lesion, excluding necrotic areas at peak contrast enhancement. Enhancement intensity, relative to adjacent normal lung tissue, was classified as hyperenhancement or hypoenhancement.

The CEUS video was analyzed frame by frame to generate the TIC. The enhancement arrival time (AT) was defined as the interval from contrast injection to the initial enhancement of the lesion. The rise time (RT) was defined as the time from initial enhancement to peak enhancement. The time to peak (TTP) was calculated as the sum of AT and RT. The peak intensity (PI) was defined as the echo intensity at peak enhancement. The wash-in area under the curve (Wi-AUC) represented the area under the TIC during the wash-in phase of contrast enhancement.

Statistical analysis and model development

All statistical analyses were performed using R software (version 4.4.3). The 75 patients were stratified into EGFR mutant and EGFR wild-type groups based on mutation status. Categorical variables were described as frequencies and percentages, while continuous variables were reported as mean±standard deviation (SD) or median with range. Categorical variables were compared using the χ² test or Fisher’s exact test. Continuous variables were analyzed using Student’s t-test or the Mann–Whitney U-test, as appropriate. A p-value of <0.05 was considered statistically significant.

Variables showing statistical significance in clinical characteristics, US features, CEUS findings, and TIC parameters were further analyzed using univariate and stepwise multivariate logistic regression. Odds ratios (ORs) were calculated for all variables included in multivariate analysis. Two predictive models and their corresponding nomograms were constructed based on the multivariate logistic regression results. The discriminative performance of models was evaluated through Receiver operating characteristic (ROC). The area under the curve (AUC) was calculated, and the optimal cutoff value was determined using the Youden index. Calibration curves were used to assess model calibration, and decision curve analysis (DCA) evaluated clinical utility.

Results

Clinical characteristics

A total of 75 patients were enrolled in this study. Patients were stratified into an EGFR wild-type group (n = 57) and an EGFR-mutant group (n = 18) based on EGFR mutation status. The clinical characteristics and corresponding statistical comparisons between the two groups were presented in Table 1. Statistically significant differences between groups were observed in sex (p = 0.001) and smoking status (p < 0.001). Meanwhile, no significant differences were observed in age (p = 0.816), TNM stage (p = 0.312), histological subtype (p = 0.133), tumor location (p = 0.861), and tumor size (p = 0.664) between EGFR+ and EGFR- groups.

Table 1
www.frontiersin.org

Table 1. Comparison of clinical characteristics of patients with or without EGFR mutation.

Development of models

The US features, CEUS characteristics, and TIC parameters were summarized in Table 2. Significant differences were identified in US features, including lesion boundaries (p = 0.017) and air bronchogram (p = 0.012). In CEUS, enhancement pattern (p = 0.034), homogeneity (p < 0.001), enhancement intensity (p = 0.026), and internal necrosis (p = 0.011) differed significantly between the two groups. Additionally, there was significant difference in RT on TIC analysis (p = 0.015).

Table 2
www.frontiersin.org

Table 2. Comparsion of US, CEUS and TIC features of lung lesions with or without EGFR mutation.

The results of univariate logistic regression for all statistically significant parameters were displayed in Table 3. The clinical characteristics (sex, smoking status) and US features (boundaries, air bronchogram) were entered into a stepwise multivariate logistic regression. Smoking status (p = 0.003), lesion boundaries (p = 0.090) and air bronchogram (p = 0.007) were identified as independent predictors. A combined C + US model was constructed using these three variables (Table 4), and its corresponding nomogram was illustrated in Figure 2.

Table 3
www.frontiersin.org

Table 3. Univariate logistic regression analysis of clinical, US and CEUS features.

Table 4
www.frontiersin.org

Table 4. Stepwise multivariate logistic regression analysis for C+US model.

Figure 2
A nomogram for predicting EGFR mutation probability using factors such as smoking status, boundaries, and airbronchogram presence. Points range from 0 to 100 for each factor, and total points convert to mutation probability from 0.1 to 0.7. Smoking status ranges from former to never smoker, boundaries from poorly to well defined, and airbronchogram absent to present. Total points align from 0 to 260.

Figure 2. Nomogram for C+US model. C, Clinical characteristics; US, ultrasound; EGFR, Epidermal Growth Factor Receptor.

Subsequently, clinical characteristics (sex, smoking status), US features (boundaries, air bronchogram), CEUS characteristics (homogeneity, enhancement intensity, necrosis), and the TIC parameter RT were analyzed via stepwise multivariate logistic regression. Lesion boundaries (p = 0.002), air bronchogram (p = 0.013), enhancement intensity (p = 0.020), internal necrosis (p = 0.002), and RT (p = 0.008) were identified as independent predictors for EGFR status. Based on these five variables, the FULL model was developed (Table 5), with a corresponding nomogram shown in Figure 3. Figures 4A–C, 5A–C demonstrated ultrasound images and pathological image of EGFR-wild type and EGFR-mutant lesions, respectively.

Table 5
www.frontiersin.org

Table 5. Stepwise multivariate logistic regression analysis for FULL model.

Figure 3
Nomogram with horizontal lines representing factors influencing EGFR mutation probability, including boundaries, air bronchogram, enhancement intensity,necrosis,and RT. Scales for total points and mutation probability range from zero to two hundred twenty and zero to 0.9, respectively.

Figure 3. Nomogram for the FULL model. RT, rise time; EGFR, Epidermal Growth Factor Receptor.

Figure 4
Figure A shows a grayscale ultrasound image, Figure B shows a contrast-enhanced ultrasound (CEUS) image, and Figure C shows the corresponding pathological image.

Figure 4. The figure 4A–4C showed a case of 70-year-old male NSCLC with EGFR-wild type. (A) The gray-scale ultrasound shown a hypo-echoic lesion with poor defined boundary. (B) On the contrast-enhanced ultrasound, the lesion exhibited hypo-enhancement with multiple necrotic focis. (C) The pathological image confirmed the lesion as a NSCLC.

Figure 5
Figure A shows a grayscale ultrasound image, Figure B shows a contrast-enhanced ultrasound (CEUS) image, and Figure C shows the corresponding pathological image.

Figure 5. The 5A–C showed a case of 50-year-old female NSCLC with EGFR-mutant. (A) The gray-scale ultrasound shown a hypo-echoic lesion with well defined boundary and air bronchograms. (B) On the contrast-enhanced ultrasound, the lesion exhibited homogeneous hyperenhancement. (C) The pathological image confirmed the lesion as a NSCLC.

Performance and validation of two models

The ROC curves for both predictive models were presented in Figure 6. The FULL model demonstrated superior diagnostic performance compared to the C + US model (AUC: 0.939 vs. 0.843, p < 0.05). The predictive metrics of both models, including specificity, sensitivity, negative predictive value, positive predictive value, and accuracy were summarized in Table 6. Calibration curves (Figures 7, 8) demonstrated high concordance of predictions with observations in both models. DCA further confirmed the clinical utility of both models. As shown in Figure 9, the FULL model provided greater net clinical benefit than the C + US model across most threshold probabilities. The FULL model demonstrated a substantial improvement in predictive performance compared to the C+US model, attributable to the integration of CEUS features and TIC parameters.

Figure 6
ROC curve comparison graph showing two models: C+US Model in red with an AUC of 0.843 and FULL Model in blue with an AUC of 0.939. The y-axis represents sensitivity, and the x-axis represents one minus specificity. The FULL Model demonstrates a superior performance curve. Points include coordinates 0.222 (0.860, 0.944) for the FULL Model and 0.224 (0.702, 0.833) for the C+US Model.

Figure 6. Receiver operating characteristic (ROC) curves of the two models with AUC of 0.843 and 0.939. C, Clinical characteristics; US, ultrasound.

Table 6
www.frontiersin.org

Table 6. Predictive performance of two models.

Figure 7
Calibration curve showing observed proportion versus predicted probability. Includes three lines: dotted for apparent, solid for bias-corrected, and dashed red for ideal. Mean absolute error is 0.052, based on seventy-five samples with one thousand bootstrap repetitions.

Figure 7. Calibration curves of C+US model. C, Clinical characteristics; US, ultrasound.

Figure 8
Calibration curve graph showing observed proportion versus predicted probability. The graph includes three lines: dotted for apparent, solid for bias-corrected, and dashed red for ideal. Key details include B equals 1000 repetitions and mean absolute error of 0.03 with n equal to 75.

Figure 8. Calibration curves of FULL model. C, Clinical characteristics; US, ultrasound; CEUS, contrast enhanced ultrasound; NSCLC, non-small cell lung cancer; TIC, time-intensity curve.

Figure 9
Chart showing a decision curve analysis with net benefit on the y-axis and high risk threshold on the x-axis. The red line represents the C+US model, the blue line represents the FULL model, and two gray lines represent “All” and “None.” The FULL model maintains a higher net benefit compared to the C+US model across most thresholds.

Figure 9. Decision curve analysis showing the net clinical benefit of the two models. C, Clinical characteristics; US, ultrasound; NSCLC, non-small cell lung cancer.

Discussion

Since targeted therapy provides meaningful survival benefits and quality-of-life improvements in advanced lung cancer, accurate detection of EGFR status before treatment is crucial throughout clinical management.

Multiple studies support the predictive value of radiomic features for EGFR genotype in NSCLC (10, 11, 15, 16). Several studies have suggested that female individuals and never-smokers are more likely to harbor EGFR mutations (1719). Additionally, prior research using computed tomography (CT) and PET/CT imaging has demonstrated that imaging features such as ground-glass opacity, convergence of surrounding vessels, pleural retraction, and adjacent invasion are correlated with EGFR mutation status (2022). Compared to CT, US, particularly CEUS, offers real-time monitoring, absence of radiation exposure, accessibility, and cost-effectiveness, making it a viable supplemental imaging tool for the differential diagnosis of subpleural pulmonary lesions. Existing studies have confirmed the high sensitivity of US and CEUS in diagnosing pulmonary and pleural lesions and in distinguishing benign from malignant pathologies (13). To our knowledge, this is the first study to predict EGFR mutation status in NSCLC with ultrasound. Our findings demonstrated that US enables noninvasive, radiation-free prediction of EGFR status with robust diagnostic performance.

In our study, univariate analysis revealed a predominance of female and non-smoking patients in the EGFR-mutant group, consistent with previous studies (10, 11). Prior research has demonstrated that adenocarcinoma is more prone to EGFR mutations. Although our cohort included squamous cell carcinoma cases, no significant association between histologic type and EGFR status was observed, possibly due to the relatively small sample size (23). In the US results, well defined boundaries and air bronchograms were exhibited predominantly in the EGFR-mutant group. This imaging phenotype may relate to the lepidic and expansile growth patterns linked to EGFR mutations. Preserved bronchial structure and air-filled spaces enable bronchogram formation, while the expansile growth pattern likely contributes to well-defined tumor boundaries (20). Regarding CEUS findings, lesions with EGFR mutations were more likely to exhibit homogeneous hyperenhancement and lack necrosis. This imaging phenomenon may reflect EGFR downstream signaling. EGFR mutations activate PI3K–AKT–mTOR pathways, upregulating hypoxia-inducible factor-1 (HIF-1), which promotes vascular endothelial growth factor (VEGF) production. VEGF induces neovascularization by influencing endothelial cells in existing vessels and promoting their migration. The development of abnormal vasculature promotes tumor growth and metastasis. Enhanced vascularization may explain the observed homogeneous hyperenhancement and reduced necrosis (24, 25).

Our study had several limitations. First, the relatively small sample size and limited number of EGFR‐mutant events may compromise the stability of the multivariable model and restrict its generalizability to broader populations. Therefore, the model is exploratory and we will further confirm the results with large samples in multiple centers in the next study. Second, the long study period (2012–2024) introduces the potential of temporal bias, as clinical management, EGFR testing techniques, and ultrasound/CEUS equipment and protocols may have evolved over time. Because the EGFR testing became popular for lung cancer patients in the recent years, we believe we can collect more data during short term period in the future. Third, beyond EGFR, other driver mutations are also relevant to targeted therapy, which need further investigate in the future.

Conclusions

The models developed in this study enable patients who are unable to undergo invasive procedures to predict EGFR mutation status noninvasively. These findings provided an ultrasound-based diagnostic reference to support clinician decision-making and personalized treatment planning.

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 humans were approved by Institutional Review Board of the Peking University School of Oncology. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

JB: Writing – original draft, Writing – review & editing. QZ: Writing – original draft, Writing – review & editing. SW: Software, Writing – original draft. HW: Investigation, Conceptualization, Writing – original draft. KY: Data curation, Methodology, Supervision, Writing – original draft. WZ: Writing – original draft, Project administration, Formal analysis. LD: Writing – original draft, Validation, Project administration. WY: Writing – original draft, Methodology, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by national natural science foundation of China (82472002), Capital Medical Development Program (2024-2-2157) and Dongying City Major Special Project for Scientific and Technological Innovation (2024ZDJH88).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

References

1. Bray F, Laversanne M, Sung H, Ferlay J, Siegel RL, Soerjomataram I, et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2024) 74:229–63. doi: 10.3322/caac.21834

PubMed Abstract | Crossref Full Text | Google Scholar

2. Li T, Kung HJ, Mack PC, and Gandara DR. Genotyping and genomic profiling of non–small-cell lung cancer: Implications for current and future therapies. J Clin Oncol. (2013) 31:1039–49. doi: 10.1200/JCO.2012.45.3753

PubMed Abstract | Crossref Full Text | Google Scholar

3. Kris MG, Johnson BE, Berry LD, Kwiatkowski DJ, Iafrate AJ, Wistuba II, et al. Using multiplexed assays of oncogenic drivers in lung cancers to select targeted drugs. JAMA. (2014) 311:1998–2006. doi: 10.1001/jama.2014.3741

PubMed Abstract | Crossref Full Text | Google Scholar

4. Lim SM, Kim EY, Kim HR, Ali SM, Greenbowe JR, Shim HS, et al. Genomic profiling of lung adenocarcinoma patients reveals therapeutic targets and confers clinical benefit when standard molecular testing is negative. Oncotarget. (2016) 7:24172–8. doi: 10.18632/oncotarget.8138

PubMed Abstract | Crossref Full Text | Google Scholar

5. An N, Zhang Y, Niu H, Li Z, Cai J, Zhao Q, et al. EGFR-TKIs versus taxanes agents in therapy for nonsmall-cell lung cancer patients. Med (Baltimore). (2016) 95:e5353. doi: 10.1097/MD.0000000000005601

PubMed Abstract | Crossref Full Text | Google Scholar

6. Hou J, Li H, Ma S, He Z, Yang S, Hao L, et al. EGFR exon 20 insertion mutations in advanced non-small-cell lung cancer: current status and perspectives. biomark Res. (2022) 10:21. doi: 10.1186/s40364-022-00372-6

PubMed Abstract | Crossref Full Text | Google Scholar

7. Lindeman NI, Cagle PT, Aisner DL, Arcila ME, Beasley MB, Bernicker EH, et al. Updated molecular testing guideline for the selection of lung cancer patients for treatment with targeted tyrosine kinase inhibitors. J Mol Diagn. (2018) 20:129–59. doi: 10.1016/j.jmoldx.2017.11.004

PubMed Abstract | Crossref Full Text | Google Scholar

8. Guo L, Chen Z, Xu C, Zhang X, Yan H, Su J, et al. Intratumoral heterogeneity of EGFR-activating mutations in advanced NSCLC patients at the single-cell level. BMC Cancer. (2019) 19:847. doi: 10.1186/s12885-019-5555-y

PubMed Abstract | Crossref Full Text | Google Scholar

9. Mu W, Jiang L, Zhang J, Shi Y, Gray JE, Tunali I, et al. Non-invasive decision support for NSCLC treatment using PET/CT radiomics. Nat Commun. (2020) 11:5228. doi: 10.1038/s41467-020-19116-x

PubMed Abstract | Crossref Full Text | Google Scholar

10. Wang S, Shi J, Ye Z, Dong D, Yu D, Zhou M, et al. Predicting EGFR mutation status in lung adenocarcinoma on computed tomography image using deep learning. Eur Respir J. (2019) 53:1800986. doi: 10.1183/13993003.00986-2018

PubMed Abstract | Crossref Full Text | Google Scholar

11. Hasegawa M, Sakai F, Ishikawa R, Kimura F, Ishida H, and Kobayashi K. CT features of EGFR-mutated adenocarcinoma of the lung: Comparison with nonmutated adenocarcinoma. J Thorac Oncol. (2016) 11:819–26. doi: 10.1016/j.jtho.2016.02.010

PubMed Abstract | Crossref Full Text | Google Scholar

12. Smith-Bindman R, Chu PW, Azman FH, Stewart C, Malekhedayat M, Alber S, et al. Projected lifetime cancer risks from current computed tomography imaging. JAMA Intern Med. (2025) 185:1–8. doi: 10.1001/jamainternmed.2025.0505

PubMed Abstract | Crossref Full Text | Google Scholar

13. Bai J, Du YQ, Yang W, Bai XM, Wang S, Wu W, et al. The role of contrast-enhanced ultrasound plus color parametric imaging in the differential diagnosis of subpleural pulmonary lesions. J Ultrasound Med. (2023) 42:2777–89. doi: 10.1002/jum.16316

PubMed Abstract | Crossref Full Text | Google Scholar

14. Yu QD, Jing B, Bing W, Song W, Fei ZQ, Kun Y, et al. Differential diagnosis of pathological type of peripheral lung cancer with multimodal contrast-enhanced ultrasound. Ultrasound Med Biol. (2024) 50:1485–93. doi: 10.1016/j.ultrasmedbio.2024.05.017

PubMed Abstract | Crossref Full Text | Google Scholar

15. Mahajan A, Kania V, Agarwal U, Ashtekar R, Shukla S, Patil VM, et al. Deep-learning-based predictive imaging biomarker model for EGFR mutation status in non-small cell lung cancer from CT imaging. Cancers (Basel). (2024) 16:1123. doi: 10.3390/cancers16061130

PubMed Abstract | Crossref Full Text | Google Scholar

16. Zhang J, Zhao X, Zhao Y, Zhang J, Zhang Z, Wang J, et al. Value of pre-therapy 18F-FDG PET/CT radiomics in predicting EGFR mutation status in patients with non-small cell lung cancer. Eur J Nucl Med Mol Imaging. (2019) 47:1137–46. doi: 10.1007/s00259-019-04592-1

PubMed Abstract | Crossref Full Text | Google Scholar

17. Dogan S, Shen R, Ang DC, Johnson ML, D’Angelo SP, Paik PK, et al. Molecular epidemiology of EGFR and KRAS mutations in 3,026 lung adenocarcinomas: Higher susceptibility of women to smoking-related KRAS-mutant cancers. Clin Cancer Res. (2012) 18:6169–77. doi: 10.1158/1078-0432.CCR-11-3265

PubMed Abstract | Crossref Full Text | Google Scholar

18. Shi Y, Au JS, Thongprasert S, Srinivasan S, Tsai CM, Khoa MT, et al. A prospective, molecular epidemiology study of EGFR mutations in Asian patients with advanced non-small-cell lung cancer of adenocarcinoma histology (PIONEER). J Thorac Oncol. (2014) 9:154–62. doi: 10.1097/JTO.0000000000000033

PubMed Abstract | Crossref Full Text | Google Scholar

19. 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:113–25. doi: 10.1056/NEJMoa1713137

PubMed Abstract | Crossref Full Text | Google Scholar

20. Zhang H, Cai W, Wang Y, Liao M, and Tian S. CT and clinical characteristics that predict risk of EGFR mutation in non-small cell lung cancer: a systematic review and meta-analysis. Int J Clin Oncol. (2019) 24:649–59. doi: 10.1007/s10147-019-01403-3

PubMed Abstract | Crossref Full Text | Google Scholar

21. Zhu P, Xu XJ, Zhang MM, and Fan SF. High-resolution computed tomography findings independently predict EGFR mutation status in ground-glass nodular lung adenocarcinoma. World J Clin Cases. (2021) 9:9792–803. doi: 10.12998/wjcc.v9.i32.9792

PubMed Abstract | Crossref Full Text | Google Scholar

22. Wang S, Yu H, Gan Y, Wu Z, Li E, Li X, et al. Mining whole-lung information by artificial intelligence for predicting EGFR genotype and targeted therapy response in lung cancer: a multicohort study. Lancet Digit Health. (2022) 4:e309–19. doi: 10.1016/S2589-7500(22)00024-3

PubMed Abstract | Crossref Full Text | Google Scholar

23. Warth A, Penzel R, Lindenmaier H, Brandt R, Stenzinger A, Herpel E, et al. EGFR, KRAS, BRAF and ALK gene alterations in lung adenocarcinomas: patient outcome, interplay with morphology and immunophenotype. Eur Respir J. (2014) 43:872–83. doi: 10.1183/09031936.00018013

PubMed Abstract | Crossref Full Text | Google Scholar

24. Sharma M, Basu DA, Nathany S, Amrith BP, and Batra U. A narrative review of the role of common EGFR mutations in pathogenesis and treatment of non-small-cell lung carcinoma. Cancer Res Stat Treat. (2022) 5:507–18. doi: 10.4103/crst.crst_222_22

Crossref Full Text | Google Scholar

25. Wang Q, Zeng A, Zhu M, and Song L. Dual inhibition of EGFR–VEGF: an effective approach to the treatment of advanced non-small cell lung cancer with EGFR mutation. Int J Oncol. (2023) 62:45–53. doi: 10.3892/ijo.2023.5474

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: contrast-enhanced ultrasound, EGFR mutation, non-small cell lung cancer, predictive model, ultrasonography

Citation: Bai J, Zhang Q, Wang S, Wang H, Yan K, Zhou W, Dong L and Yang W (2026) Multimodal ultrasonography for predicting epidermal growth factor receptor mutation in subpleural non-small cell lung carcinoma. Front. Oncol. 15:1722238. doi: 10.3389/fonc.2025.1722238

Received: 10 October 2025; Accepted: 11 December 2025; Revised: 11 December 2025;
Published: 07 January 2026.

Edited by:

Hang Zhou, Zhejiang University, China

Reviewed by:

Yao Luo, Sichuan University, China
Abir Swaidan, University of California, Los Angeles, United States

Copyright © 2026 Bai, Zhang, Wang, Wang, Yan, Zhou, Dong and Yang. 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: Wei Yang, MTM2ODE0MDgxODNAMTYzLmNvbQ==; Liang Dong, ZG9uZ2xpYW5nMTk4MzYxMEAxNjMuY29t

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.