Abstract
Introduction
Primary lung cancer is a malignant tumor that commonly originates from the epithelial cells of the airways, especially in the trachea and bronchus. 1 Based on the most recent data released by the World Health Organization in 2023, lung cancer continues to be the primary contributor to cancer-related fatalities, with approximately 350 individuals succumbing to the disease daily, nearly 2.5 times the mortality rate of the second leading cause of cancer-related death. 2
Thoracic surgery has historically been a cornerstone of the treatment of lung neoplasms. When considering surgical interventions, evaluating several factors, including the patient’s lung function, performance status, and tumor characteristics, such as size, number, and involvement of the trachea and bronchus, is important. Unfortunately, approximately 70% of lung cancer cases were diagnosed at advanced stages, often with distant metastasis.3,4 Moreover, over 80% of cases were identified as non–small-cell lung cancer (NSCLC), which is unsuitable for surgical intervention. 5
For locally advanced NSCLC, which is commonly treated in radiotherapy departments, great progress has been made in terms of clinical efficacy in recent years. The 5-year survival rate of radical concurrent chemoradiotherapy plus immune consolidation was 42.9%, which is encouraging. For patients who cannot tolerate concurrent chemoradiotherapy, radiotherapy alone or sequential chemoradiotherapy followed by consolidation immunotherapy was also the recommended treatment mode. 6
Epidermal growth factor receptor (EGFR) is one of the most mature targets in NSCLC research. Over the past decade, targeted therapeutic strategies based on driver genes have proven effective in the treatment of NSCLC, significantly improving the prognosis of patients with driver gene-positive advanced NSCLC. In the meantime, several researchers have verified that the EGFR signaling pathway can promote vascular endothelial growth factor (VEGF) expression in tumors, 7 which may increase the sensitivity of patients with EGFR mutations to bevacizumab. 8 Increasing research indicates that the combination of EGFR-tyrosine kinase inhibitors (TKIs) with bevacizumab achieves sound therapeutic effects. 9 Although in a series of clinical studies, these 2 drugs have significantly improved the overall therapeutic effect, progression-free survival (PFS), and quality of life of cancer patients, long-term effectiveness is often difficult to achieve when either drug is used alone because of the development of resistance. 10 Consequently, the objective of this meta-analysis was to assess the effectiveness and safety of using EGFR-TKIs in conjunction with bevacizumab in treating patients with newly identified or recurring EGFR mutation-positive advanced NSCLC who did not receive previous chemotherapy, providing an evidence-based reference for the selection of first-line treatment options for this patient group in clinical settings.
Materials and Methods
This protocol was registered with the International Platform of Registered Systematic Review and Meta-Analysis Protocols (INPLASY) on December 14, 2023, and was last updated on December 14, 2023 (registration number INPLASY2023120059).The reporting of this study adheres to Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA; Supplemental Material) guidelines. 11
Study Inclusion Criteria
1. Clinical trials of EGFR-TKIs plus Bevacizumab in the treatment of NSCLC harboring EGFR-sensitizing mutations.
2. Individuals identified with EGFR-mutant NSCLC through pathology or clinical assessments, who had not undergone surgery or other treatments prior to joining the clinical study.
3. After statistical analysis, there was no significant difference between the 2 groups in terms of age, gender, stage, pathological type, gene mutation status, and other aspects (
4. The experimental group received EGFR-TKIs plus Bevacizumab, while the control group received EGFR-TKIs alone.
5. Outcome indicators: complete response (CR), partial response (PR), objective response rate (ORR), and disease control rate (DCR) which was according to the internationally accepted criteria for evaluating the efficacy of solid tumors (RECIST 1.1). And overall survival (OS) of 1, 2, or 3 years was used as primary efficacy outcomes. Also includes PFS of 1, 2, or 3 years used to explore the efficacy and safety of combined medications. The secondary outcome measure was the number of adverse reactions, including diarrhea, rash, anorexia, proteinuria, pneumonitis, and hypertension.
The PFS is the time from randomization to tumor progression or death; OS is the time from randomization to death, which is considered the best treatment endpoint in cancer clinical trials; ORR is the proportion of patients whose tumors shrank to a certain degree and remained unchanged for a while, including cases of CR, PR, progressive disease (PD), and stable disease (SD).
6. Published RCTs (randomized clinical trials).
7. Medical records are complete.
8. These articles conform to ethical norms.
Study Exclusion Criteria
Non-interventional clinical trials or observational or retrospective studies.
Studies which enrolled patients had had other treatments before enrollment.
Studies with inconsistent intervention measures. In addition, single-arm studies and studies with inconsistent design and outcome indicators should be excluded.
Studies like review, animal or basic experiment, and so on should be excluded.
Non-randomized controlled trials.
Studies still being carried.
Search Strategy and Study Selection
We comprehensively searched for clinical studies of EGFR-TKIs combined with Bevacizumab for EGFR-mutant NSCLC in the PubMed, Cochrane Library, Web of Science, Embase, Wanfang Data, Chinese National Knowledge Infrastructure (CNKI), Chinese Biological Medicine (CBM) Database, and VIP Database until October 2023. We retrieved the following terms from the Medical Subject Headings (MeSH): Non small cell lung cancer; EGFR mutation; Bevacizumab; and VEGF inhibitor. Independently, a pair of researchers (XQ and ZW) meticulously examined the literature, adhering strictly to the set criteria for inclusion and exclusion. Should there be a discrepancy, the decision to incorporate the study rests with the third researcher.
Data Extraction and Quality Assessment
Relevant data were extracted independently by 2 authors (XQ and ZW) according to predesigned data tables, including authors, year of publication, country, sample size, age, sex ratio, interventions, CR/PR/ORR/DCR rates, 1-, 2-, or 3-year OS, and 1-,2-, or 3-year PFS after treatment in the experimental and control groups. Two independent authors strictly assessed the quality of RCT studies according to the tools provided by the Cochrane Collaboration Network, using MINORS entries to evaluate the quality of non-randomized clinical trials.
Statistical Analysis
The meta-analysis of data from the studies was conducted using Review Manager (RevMan) versions 5.3.5 and Stata MP versions 14.0 and 16.0 (Stata Corporation, College Station, Texas). Results from binary data will undergo analysis in the form of odds ratios (ORs), while continuous data will be examined as mean difference (MD) or standardized MD, with 95% confidence intervals (CIs) provided for every estimate. The diversity among studies was assessed through the application of the Q statistic (Chi-square test) and the Higgins I2 index. Heterogeneity was deemed negligible when
Results
Characteristics of studies
We initially retrieved a total of 547 publications, including 214 duplicates, 178 studies with irrelevant content, 115 reviews and case reports. On reviewing the entire text, 26 studies were omitted due to their irregular experimental designs, including 1 study with data entry errors, 7 studies still being carried, 1 study focusing on designing RCT trials, and 17 studies with wrong disease stage included. Fourteen studies were ultimately selected for meta-analysis, and the literature screening process is shown in Figure 1. In total, 1873 patients were included, 934 in the experimental group and 939 in the control group, and 1 patient in the control group quitted before treating. Fourteen articles were all based on randomized controlled clinical trials. The experimental group’s patients were treated with both EGFR-TKI and Bevacizumab, in contrast to the control group’s patients, who solely received EGFR-TKI therapy. In the included literature, 2 studies choose Gefitinib among EGFR-TKIs, 1 study chooses Afatinib, 8 studies choose Erlotinib, 3 studies choose Osimertinib, and the dosage use of EGFR-TKIs is as directed by the instructions. Most studies’ dosage of Bevacizumab is 15 mg/kg, q3w. A summary of the detailed characteristics for each included study is provided in Table 1, while Table 2 systematically displays the baseline patient characteristics.

Flow chart of study screening.
Study design and efficacy outcomes of included studies.
Abbreviations: 1L: first line; 2L: second line; Bev + TKI: Bevacizumab + EGFR-TKI (combination therapy); TKI: EGFR-TKI (monotherapy); NR: not reached.
Indicates mean value.
Baseline patient characteristics.
Quality assessment
The results of the quality evaluation of the 14 RCT studies are shown in Figure 2. All of the studies were RCTs. Two studies assigned simple random sampling, and 4 studies assigned stratified sampling. Although none of the remaining studies described their methods, they do declare that patients were randomly assigned to 2 groups. Most of the included studies did not provide sufficient information to assess whether the allocation concealment was adequate, except Yosuke Kawashima, 2021. Each of the included studies provided comprehensive examination data with no evidence of selective reporting or bias. Minors (methodological index for non-randomized studies) was used to assess the quality of the 14 RCTs, quality analysis showed that all the included studies were of high or moderate quality.

Risk of bias summary (A) and bias graph (B).
Efficiency
CR
Seven studies12,15,17,20,22,24,25 reported CR rates. There was no significant heterogeneity between studies (

Forest plot for CR of EGFR-TKIs combined with Bevacizumab group and EGFR-TKIs alone group.
PR
Seven studies12,15,17,20,22,24,25 reported PR rates. There was no significant heterogeneity between studies (

Forest plot for PR of EGFR-TKIs combined with Bevacizumab group and EGFR-TKIs alone group.
SD
Seven studies12,15,17,20,22,24,25 reported SD rates. There was no significant heterogeneity between studies (

Forest plot for SD of EGFR-TKIs combined with Bevacizumab group and EGFR-TKIs alone group.
PD
Seven studies12,15,17,20,22,24,25 reported PD rates. Significant heterogeneity was observed across these studies (

Forest plot for PD of EGFR-TKIs combined with Bevacizumab group and EGFR-TKIs alone group.
ORR
Seven studies12,15,17,20,22,24,25 reported ORR rates. There was no significant heterogeneity between studies (

Forest plot for ORR of EGFR-TKIs combined with Bevacizumab group and EGFR-TKIs alone group.
DCR
Seven studies12,15,17,20,22,24,25 reported DCRs. Significant heterogeneity was observed among the included studies (

Forest plot for DCR of EGFR-TKIs combined with Bevacizumab group and EGFR-TKIs alone group.
1-year, 2-year, and 3-year PFS
Eleven studies12-14,16-23 reported 1-year PFS rates; 9 studies12-14,16,18-22 reported 2-year PFS rates; and 5 studies12-14,16,19 reported 3-year PFS rates. As the Hiroaki Akamatsus study showed a substantial difference from the others, it was excluded from the analysis. After exclusion, the remaining studies exhibited no significant heterogeneity for 2- and 3-year PFS, and thus, a fixed-effects model was applied for these outcomes (2-year PFS:

Forest plot for 1-year (A), 2-year (B), and 3-year (C) progression-free survival rates of EGFR-TKIs combined with Bevacizumab group and EGFR-TKIs alone group.
1-year, 2-year, and 3-year OS
Ten studies12-14,16-21,23 reported 1-year OS rates; 10 studies12-14,16-21,23 reported 2-year OS rates; and 8 studies12-14,16,18-21 reported 3-year OS rates. These studies did not show heterogeneity, so a fixed-effects model was employed (1-year OS:

Forest plot for 1-year (A), 2-year (B), and 3-year (C) survival rates of EGFR-TKIs combined with Bevacizumab group and EGFR-TKIs alone group.
Adverse effects
In the analysis of adverse events, we uniformly applied the fixed-effect model to all groups that demonstrated no significant heterogeneity. Eight studies reported cases of diarrhea. The meta-analysis showed a significantly higher incidence of diarrhea in the experimental group compared to the control group (OR = 1.34, 95% CI = 1.05-1.72,

Forest plot for adverse reactions of EGFR-TKIs combined with Bevacizumab group and EGFR-TKIs alone group. Anorexia (B): After subgroup analysis excluding the study by Takashi Ninomiya, heterogeneity was resolved, and results are derived from a fixed-effect model.
Evaluation of the sensitivity and publication bias
We performed a sensitivity analysis by excluding the included studies sequentially from the pooled effect and estimated the effect of each study on the overall outcome by observing the change in outcome after deletion. The results showed that the significances of OR did not change after the exclusion of any of the studies, suggesting the reliability and reasonableness of our meta-analysis results. We plotted the funnel diagram by STATA (Figure 12) and performed Begg and Egger funnel plots to detect publication bias. The results showed that the funnel plot of PR and ORR showed publication bias; this may lead to an undue overestimation of the efficacy of the combination therapy, and we plotted the funnel plot of the cut-and-patch method (Figures 13 and 14), which showed that in the future we could eliminate publication bias by adding 3 more in PR and 1 more in ORR. Other funnel plots revealed an absence of notable bias in publication with

Funnel plots.

Begg funnel plot (A) and sensitivity analysis (B) of PR. Begg test (

Begg funnel plot (A) and sensitivity analysis (B) of ORR. Begg test (
Result of Egger regression tests.
Sensitivity analysis of pooled odds ratios using the trim and fill method.
Discussion
The EGFR is a glycoprotein belonging to the receptor tyrosine kinase family. They also play crucial roles in cell growth and proliferation. It is activated upon binding to its corresponding ligand, forming a dimer, and participating in signal transduction. Abnormal EGFR function, including gene mutations or overexpression, leads to the sustained activation of EGFR and uncontrolled cell growth and proliferation, thereby inducing carcinogenesis. 26 Clinically, patients with EGFR mutations account for 40% to 60% of cases in Asia, with 90% of these mutations being exon 19 deletions (EX19 DEL) or exon 21 leucine to arginine substitutions (EX21 L858R).27,28 For patients with EGFR mutations, EGFR-TKIs have become first-line treatment options. Tyrosine kinase inhibitors act as competitive inhibitors for the binding of adenosine triphosphate (ATP) to tyrosine kinase or as analogs of tyrosine, blocking the activity of tyrosine kinase and inhibiting cell proliferation while avoiding the adverse effects of conventional chemotherapy. 29 Currently, the clinically used EGFR-TKIs include erlotinib, gefitinib, afatinib, and second-generation osimertinib. Among these, gefitinib, which can effectively inhibit angiogenesis, hinder tumor growth and metastasis, and promote tumor cell apoptosis, is the most commonly used. However, long-term clinical practice has proven that long-term single use of EGFR-TKIs can easily lead to acquired resistance. Common resistance mechanisms include MET amplification, HER-2 overexpression, and histological transformation into small-cell lung cancer. This has prompted researchers to search for novel pathways.
Angiogenic factors, such as VEGF, play a central role in tumor angiogenesis, so a large amount of research has focused on developing VEGF inhibitors, among which bevacizumab is the most famous. Bevacizumab is a recombinant humanized monoclonal antibody targeting VEGF, which specifically binds to VEGF, blocking the binding of VEGF to vascular endothelial growth factor receptor (VEGFR), interrupting the signal transduction pathway of angiogenesis, and inhibiting the formation of new blood vessels in tumors, thereby inhibiting tumor cell growth and exerting an anti-tumor effect.30,31 The ECOG4599 study was the first Phase III randomized controlled trial (RCT) to include bevacizumab as a first-line treatment regimen for advanced NSCLC. The experimental group received maintenance treatment with bevacizumab in addition to the chemotherapy control group. Significant differences were in the improvements in OS, PFS, ORR, and other indicators between the bevacizumab and chemotherapy groups. This is the first clinical study to extend the median OS of patients with advanced NSCLC to >12 months. The incidence of hypertension, proteinuria, bleeding, neutropenia, febrile neutropenia, thrombocytopenia, hyponatremia, rash, and headache in the bevacizumab group was significantly higher than that in the chemotherapy group (
Based on these questions, we conceived this study to explore the efficacy of combination therapy, which was a systematic review and meta-analysis of RCTs to compare the combination of Bev and EGFR-TKIs to EGFR-TKIs alone in the treatment of NSCLC harboring EGFR-sensitizing mutations. In this meta-analysis, we employed the trim-and-fill method to address potential publication bias and heterogeneity. This approach indicated that certain studies (eg, the removal of Haruhiro Saito et al from the PD analysis, and Su DW et al from the DCR analysis) might represent outliers. The discrepancies between these studies and the others primarily stemmed from their distinct patient populations or study designs. For instance, the study by Haruhiro Saito et al enrolled a higher proportion of patient subgroups with poor prognosis (such as those with a greater incidence of brain metastases or inferior performance status), which may have contributed to its significantly higher progressive disease (PD) rate compared to other studies. Similarly, the atypical DCR observed in the study by Su DW et al could be related to the specific response evaluation criteria applied or distinct mechanisms of treatment response within their cohort. Regarding the studies by Hiroaki Akamatsu et al (excluded from PFS analysis) and Takashi Ninomiya et al (excluded from the anorexia analysis), the observed heterogeneity may originate from differences in follow-up duration and progression assessment criteria, as well as variations in the recording and reporting standards for adverse events, respectively. These subtle methodological and population differences, objectively identified through the trim-and-fill method, explain the sources of the observed heterogeneity. After excluding these studies, the pooled results of the remaining studies demonstrated good consistency, thereby strengthening the robustness of our primary conclusions. Our study demonstrated that the combination treatment significantly prolonged 1- and 2-year PFS and 1-year OS, also elevated the PR, PD, and ORR in the treatment of NSCLC harboring EGFR-sensitizing mutations. Bev combined with EGFR-TKIs did not prolong 3-year PFS and 2- and 3-year OS in the treatment of NSCLC. Our study has similar results with 1 published network meta-analysis of Bev and Erlotinib for patients with advanced NSCLC. In short time, the therapeutic schedule of EGFR-TKIs plus Bev apparently prolongs patients’ survival, but in the long run, the experimental group does not show better efficacy than the control group. There are many reasons for it. First, most of our patients were found with advanced NSCLC, some even also with brain or bone metastases, their scores of ECOG (Eastern Cooperative Oncology Group) were low, and it is hard to reverse the death trend completely. Second, there is basic research showing that the expression levels of VEGF undergo dynamic changes during targeted therapy, and resistance to EGFR-TKIs is often accompanied by an increase in VEGF levels. 34 Third, significant differences in prognosis exist among the different mutation types, which may also contribute to these results. The EGFR mutations are also divided into several types: among the 3 types with higher mutation frequency, exon 19 deletion, exon 20 T790M, and Exon 21 L858, there was a significant difference between different mutation types and metastatic sites, with a higher percentage of bone metastases (68.18%) occurring in Exon19 deletion and a higher percentage of brain metastases (56.00%) in exon 20; T790M had a higher percentage of brain metastases (56.00%). 35 Given that all selected studies were real-world RCTs, certain biases inherently exist in their implementation—such as inadequate allocation concealment and difficulties in blinding procedures. Concurrently, RCTs reporting positive results (indicating intervention efficacy) are more likely to be published than those with negative or null results, potentially leading to overly optimistic estimates of treatment effects in the literature. Furthermore, our meta-analysis relied on aggregate published data rather than individual patient data (IPD), which may introduce additional bias. Notably, the specific EGFR-TKIs investigated varied across studies, with only a few RCTs available for each agent. Many included trials also featured immature follow-up data.
With regard to adverse reactions, the occurrences of diarrhea, anorexia, and proteinuria were more frequent in the group treated with EGFR-TKIs combined with Bev than in the control group. To our surprise, the pneumonitis and hypertension rates showed significant heterogeneity between the 2 groups, which may be caused by the low sample size, but there is still a lot of potential to cause these symptoms due to combination treatment. To our surprise, the rate of rash was almost the same in the experimental and control groups, and all the adverse effects are common in clinical practice so that we can deal with them easily.
Our meta-analysis confirms and extends the findings of prior systematic reviews, such as the study by Yang et al, 36 which established the efficacy of combining bevacizumab with EGFR-TKIs for improving PFS in advanced EGFR-mutant NSCLC. Our study adds to the existing literature by providing a substantially more comprehensive and updated synthesis of the evidence. By incorporating 14 randomized trials and nearly doubling the patient population, our analysis offers greater statistical power and robustness. This expanded dataset allows us to demonstrate a significant improvement in 1-year OS, a key survival benefit not firmly established in previous meta-analyses.37-39 Furthermore, our more granular safety analysis refines the risk profile of the combination therapy, indicating that it does not significantly increase the incidence of rash or pneumonitis. Therefore, our work strengthens the evidence base for the efficacy of this combination and provides a more precise and contemporary assessment of its survival benefits and safety profile for clinical decision-making.
A randomized, open-controlled, single-center study was initiated by Zhao et al to compare Gefitinib + Bev vs gefitinib alone as the first-line treatment for advanced EGFR L858R NSCLC. This study is currently in progress in enlisting participants, and the results are highly anticipated. Further studies are needed to explore the potential of EGFR-TKIs plus Bev in the treatment of NSCLC harboring EGFR-sensitizing mutations. 40
Conclusions
This meta-analysis demonstrates that, as a first-line treatment for advanced EGFR-mutant NSCLC, the combination of bevacizumab and EGFR-TKIs significantly improves PR rate, ORR, 1-year PFS, and 1-year OS compared to TKI monotherapy while effectively delaying disease progression. Although the combination increases the risks of diarrhea, hypertension, and proteinuria, it does not significantly raise the incidence of rash or pneumonitis, presenting a manageable safety profile. Thus, this combination represents a valuable first-line therapeutic option that provides substantial short-to-medium-term survival benefits.
Limitations
Some questions remain regarding our study. Subgroup analysis of these studies is difficult due to differences in patient enrollment, so the issues requiring further study include the identification of the best EGFR-TKI OR Bevacizumab regimen, different effects in smokers and non-smokers, and so on. Moreover, most of the included studies did not include long-term follow-up. Additional studies are required to evaluate the relationship between EGFR-TKI plus bevacizumab therapy and long-term survival.
Supplemental Material
sj-docx-1-onc-10.1177_11795549251414657 – Supplemental material for Bevacizumab Plus EGFR-TKIs vs EGFR-TKIs Alone for Advanced EGFR-Mutant NSCLC: A Meta-Analysis
Supplemental material, sj-docx-1-onc-10.1177_11795549251414657 for Bevacizumab Plus EGFR-TKIs vs EGFR-TKIs Alone for Advanced EGFR-Mutant NSCLC: A Meta-Analysis by Zexian Wang, Yaru Guo, Xiaohan Qin, Zhiling Wan, Xiaojin Wu and Chen Liu in Clinical Medicine Insights: Oncology
Footnotes
Ethical Considerations
Consent to participate
Author contributions
Funding
Declaration of conflicting interests
Data availability statement
Supplemental material
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
