- 1Department of Radiation Oncology, Affiliated Cancer Hospital and Institute of Guangzhou Medical University, Guangzhou, China
- 2Health Ward, Affiliated Cancer Hospital and Institute of Guangzhou Medical University, Guangzhou, China
- 3Department of Radiation Oncology, Guangzhou Concord Cancer Center, Guangzhou, China
Background: The optimal second-line systemic treatment model for recurrent and/or metastatic head and neck squamous cell carcinoma (R/M HNSCC) remains controversial. A Bayesian network meta-analysis (NMA) was performed to address this issue with regard to efficacy and toxicity.
Methods: By searching MEDLINE (via PubMed), Embase, the Cochrane Central Register of Controlled Trials and Web of Science, we extracted eligible studies. Efficacy, represented as overall survival (OS) and progression-free survival (PFS), and overall toxicity, represented as ≥ grade 3 severe acute events (sAE), were assessed to compare the following 7 treatment models through an NMA: standard-of-care therapy (SoC), single targeted therapy different from SoC (ST), double targeted therapy (DT), targeted therapy combined with chemotherapy (T+C), single immune checkpoint inhibitor therapy (SI), double immune checkpoint inhibitor therapy (DI) and single chemotherapy different from SoC (SC). Rank probabilities according to the values of the surface under the cumulative ranking curve (SUCRA) were separately determined for efficacy and toxicity.
Results: In total, 5285 patients from 24 eligible studies were ultimately screened, with 5184, 4532 and 4026 involved in the NMA of OS, PFS and sAE, respectively. All qualifying studies were absent from first-line immune checkpoint inhibitor therapy. In terms of OS, SI was superior to the other treatments, followed by DI, ST, T+C, SoC, DT and SC. Other than SI and SC, all treatments tended to be consistent, with hazard ratios (HRs) close to 1 between groups. For PFS, ST ranked first, while DT ranked last. For the toxicity profiles, compared with the other models, SI resulted in the lowest incidences of sAE, with statistical significance over SoC (odds ratio [OR] 0.31, 95% credible interval [CrI] 0.11 to 0.90), ST (OR 0.23, 95% CrI 0.06 to 0.86) and DT (OR 0.11, 95% CrI 0.02 to 0.53), while DT was the worst. When the SUCRA values of OS and sAE were combined, a cluster plot illustrated the superiority of SI, which demonstrated the best OS and tolerability toward sAE.
Conclusion: For R/M HNSCC patients without immune checkpoint inhibitors in the first-line setting, SI may serve as the optimal second-line systemic treatment model, demonstrating the best OS and least sAE.
Highlights
This network meta-analysis, encompassing 5285 individuals from 24 trials which were absent from first-line immune checkpoint inhibitor therapy, revealed that single immune checkpoint inhibitor therapy achieved the best overall survival as well as the least severe acute events in the second-line treatment of R/M HNSCC patients who were unable or did not have access to immunosuppressants in the first-line setting. It provides recommendations for future clinical decision-making and trial designs.
Introduction
Head and neck squamous cell carcinoma (HNSCC), originating from the epithelial tissue of the head and neck region, includes mainly oral cancer, oropharyngeal cancer, hypopharyngeal cancer and laryngeal cancer. As estimated by the International Agency for Research of Cancer, there were approximately 830,000 patients newly diagnosed with HNSCC worldwide in 2018, and approximately 430,000 patients died of this type of tumor throughout the year (1).
HNSCC is a biologically and clinically heterogeneous disease, and its prognosis is quite less than satisfactory (2). Approximately 4% ~ 26% of patients present with metastatic disease at the first diagnosis, and among patients who have not been diagnosed with metastasis, despite radical treatment, more than 30% will eventually develop recurrent and/or metastatic (R/M) diseases (3–8). Treatment options for R/M HNSCC are limited, especially in regard to second-line regimens, mainly referring to chemotherapy, targeted therapy or a combination of the two. Adoption of these regimens leads to median progression-free survival (PFS) and overall survival (OS) times of only 2 ~ 3 and 6 ~ 8 months, respectively (9–11). Encouragingly, the advent of immunotherapy has led to survival hope for patients with metastatic tumors in recent years (12). However, its application is controversial in the administration of different immune checkpoint inhibitors (ICIs) in HNSCC (13–16). In the CheckMate 141 and KEYNOTE 040 trials, the results showed that ICIs significantly increased OS over standard therapy as second-line treatments (13, 14). Therefore, both nivolumab and pembrolizumab have been recommended by the National Comprehensive Cancer Network (NCCN) guidelines. Notably, in the EAGLE trial, which explored the use of single and double ICIs as second-line regimens for R/M HNSCC, no statistically superior efficacy was observed with respect to standard therapy (16). This opposing outcome poses a challenge to the earlier research. In addition, dual ICIs exhibited no survival benefit over single applications in the EAGLE trial, which is consistent with the results of the CONDOR trial (15). Hence, controversy remains in regard to whether ICIs can be used as the optimal treatment and to how it should be used in combinations.
KEYNOTE 048 is a notable study evaluating the efficacy of pembrolizumab monotherapy or with chemotherapy in comparison to the Extreme regimen, which marked a remarkable survival superiority when pembrolizumab was administered (17). Its advent facilitates the widespread application of pembrolizumab in the first-line treatment of R/M HNSCC, however, pembrolizumab is presently not so affordable and there are many developing countries where ICIs are not yet approved in the front line setting and only in the second-line setting owing to the fact that treatment with pembrolizumab-related regimens is a cost-effective strategy in the developed country like the USA, whereas in the developing country like China, it is the opposite and the Extreme regimen still stands out as a cost-effective choice (18, 19). In patients who have experienced ICIs in the first-line setting, the subsequence choices will definitely catch public attention, however, this would be a longer-than-expected event in certain patient populations since preferring pembrolizumab in the first-line is not a cost-effective strategy in certain developing countries (18, 19). Under these circumstances, in patients who are absent from first-line immunotherapy, the second-line choice concerning ICIs is still worthy of explorations. Further determination of the exact role immunotherapy plays in second-line treatment is imperative.
Since few trials have directly compared immunotherapy with other treatments, this further exploration of the best second-line systemic treatment model for R/M HNSCC can be addressed by conducting a network meta-analysis (NMA), which can directly and/or indirectly compare different groups from the current literature and, to a certain extent, make up for the deficiency in the limited treatment comparisons in clinical trials to date in order to determine the best treatment group (20). Therefore, by performing an NMA, we conducted this work to determine the optimal second-line systemic treatment model for R/M HNSCC in terms of efficacy and toxicity.
Materials and Methods
Abiding by the Preferred Items for Systematic Reviews and Meta-Analyses (PRISMA) reporting guidelines as well as the extension statement of NMA (21, 22), this work was undertaken, and the relevant protocol was registered in INPLASY (No. 2020110041). The Institutional Review Board of Affiliated Cancer Hospital and Institute of Guangzhou Medical University granted this review exemption due to the innocuousness of this research, and no ethical approval was required. All necessary data were extracted from online reported studies on MEDLINE (via PubMed), Embase, the Cochrane Central Register of Controlled Trials and Web of Science. Two reviewers (Zhan ZJ, Zhang F) independently executed the subsequent search process. Whenever disagreements arose, a consultation referring to a senior investigator (Qiu WZ) was carried out to achieve a consensus.
Literature Search
Supplementary Table 1 describes in detail the search strategy for this NMA. First, up to October 2020, which was finally updated on June 2021, screening was performed via the relevant research titles and abstracts, and duplicated publications were removed. Then, for the assessment for final inclusion, we thoroughly browsed the full text of the articles as well as the corresponding reference lists. Those publications that did not strictly meet the following inclusion criteria were excluded.
Study Selection
The inclusion criteria were as follows:
a. Randomized controlled trials (RCTs);
b. Patients with pathologically confirmed R/M HNSCC;
c. Use of treatments including chemotherapy, targeted therapy, immunotherapy or a combination thereof;
d. Survival data and/or toxicity profiles available in the study or could be calculated with the relevant reported data.
Exclusion criteria were as follows:
a. Pathologically confirmed nasopharyngeal carcinoma as the only type of tumor included in the whole cohort;
b. Use of local treatments such as surgery and radiotherapy in a second-line setting;
c. A second primary HNSCC;
d. Studies in which toxicity profiles were not presented as ≥ grade 3 toxicity overall based on National Cancer Institute-Common Terminology Criteria for Adverse Events (NCI-AECTC) (available at http://ctep.cancer.gov);
e. Studies in which full-text reports were unavailable.
Data Extraction and Risk of Bias Assessment
We collated the relevant data from all eligible studies onto an electronic spreadsheet, which mainly included first author, year of publication, sample sizes, institutions, trial design, interventions and outcome measures. Hazard ratios (HRs) were determined for efficacy outcome assessments of OS/PFS, and odds ratios (ORs) were used for safety profile measures of treatment-related severe acute events (sAE), along with their 95% credible intervals (CrI). OS and PFS were defined as the interval from randomization to last follow-up or death and the interval from randomization to tumor progression or death, respectively. The sAE was defined as ≥ grade 3 overall toxicity. For studies lacking HRs, the Engauge Digitizer (version 4.1) and a calculation spreadsheet were utilized to estimate HRs through Kaplan-Meier curves along with their 95% CrI according to the methods described by Tierney and colleagues (23). If the related outcomes were unavailable in the study, an attempt was made to consult the study authors to obtain individual patient data of interest.
Adhering to the Cochrane Risk of Bias Tool (24), the risk of bias of the included RCTs was assessed. It includes 7 domains, namely, random sequence generation, allocation concealment, blinding of participants and personnel, blinding of outcome assessment, incomplete outcome data, selective outcome reporting and other sources of bias.
Statistical Analysis
The NMA is a mixed-treatment comparison meta-analysis that simultaneously integrates both direct and indirect comparisons. For example, if a direct comparison between treatments A and C is lacking, a closed loop (namely, a network) consisting of A-B-C-(A) can be established under a Bayesian framework if treatments A and C were commonly compared with treatment B across trials, enabling the identification of the relative effect between any two treatments (22, 25). In the present work, the primary endpoint was OS in the available patient cohort, and the secondary endpoints included PFS and sAE. For outcomes for analysis, HRs and ORs along with their 95% CrI were collected or calculated for OS/PFS and sAE, respectively.
By adopting the model of the lower deviance information criterion (DIC), which indicates better feasibility (26), a Bayesian NMA of OS and PFS was performed with WinBUGS (version 1.4.3); for treatment-related sAE, it was performed with the gemtc package (version 0.8-2) (27), which was implemented in R (X64 version 3.5.1). By comparing the surface under the cumulative ranking curve (SUCRA) generated by these tools, we can obtain evidence to determine the rank probability of the competing arms to rank them as first, second, third and so on (25, 28, 29). The highest SUCRA value denotes the highest probability of ranking the best, while the lowest value denotes the highest probability of ranking worst. To fit the model, parameters were set as follows for the analyses in accordance with noninformative uniform and normal prior distributions (30): 3 chains with dispersed initial values; 100,000 iterations for each chain for a total of 300,000 in the three chains for the posterior distributions; 5000 burn-ins; and an interval of 10 iterations per chain.
Network and funnel plots were generated with STATA (version 14.0) to visually illustrate the relationships of each treatment model of the available trials and evaluate the studies for publication bias (31), respectively. The inconsistency statistic (I2) (32) was calculated to quantify the global inconsistency across studies (i.e., between-study heterogeneity). An I2 < 25% indicated low heterogeneity, whereas 25% ≤ I2 ≤ 50% and I2 > 50% were deemed to be moderate and high heterogeneity, respectively (32). Local inconsistency, which may be generated in the network, was tested by the “node-splitting” technique (33) and served as a tool in the comparison between direct and indirect evidence among the entire network (34). Statistical significance was considered at p < 0.05.
Results
Eligible Studies
Figure 1 demonstrates the screening of eligible studies. After searching the online databases, 3803 records were identified, with 1784, 1381, 597 and 41 recognized in the MEDLINE, Embase, Cochrane Central Register of Controlled Trials and Web of Science databases, respectively. 688 duplicate records were removed. Then, 3115 records were extensively screened via their titles and abstracts as extracted from the online databases. Twenty-four RCTs with 5285 participants conformed to the inclusion and exclusion criteria and thus were finally screened (9–11, 13–16, 35–51). Among the eligible studies, which consisted of treatment models of standard-of-care therapy (SoC), single targeted therapy different from SoC (ST), double targeted therapy (DT), targeted therapy combined with chemotherapy (T+C), single immune checkpoint inhibitor therapy (SI), double immune checkpoint inhibitor therapy (DI) and single chemotherapy different from SoC (SC), SoC was the most common controlled arm with the largest number participants, followed by ST, DT, T+C, SI, DI and SC, as indicated in the network plot in Figure 2. Of note, none of the eligible RCTs applied immunosuppressant therapy in the first-line setting and the treatment model of T+C in the second-line setting represents single targeted therapy combined with single chemotherapy. Supplementary Table 2 summarizes the characteristics of all eligible RCTs; the HRs of the OS and PFS data and the overall sAE profiles could not be extracted from 1, 2, and 6 of these trials, respectively. As a result, 5184, 4532 and 4026 participants were finally selected for OS, PFS and sAE analyses, respectively. The network plots of the survival and sAE analyses are separately shown in Supplementary Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719650/fimmu-12-719650-HTML-r1/image_m/fimmu-12-719650-g002.jpg)
Figure 2 A network plot of all eligible trials. SoC, standard-of-care therapy; ST, single targeted therapy different from SoC; DT, double targeted therapy; T+C, targeted therapy combined with chemotherapy; SI, single immune checkpoint therapy; DI, double immune checkpoint therapy; SC, single chemotherapy different from SoC.
Quality of Evidence
The assessments of all eligible studies for risk of bias are presented in Supplementary Table 4 and Supplementary Figure 2; two trials were detected as having a high risk of bias. In the inconsistency analyses, all global analyses concerning OS, PFS and sAE showed low heterogeneity (I² = 14%, I² = 8%, I² = 5%, respectively). Further exploration of local inconsistency revealed that except for the p value between SoC and DI in the OS analysis, all other p values between direct and indirect evidence were greater than 0.05, suggesting that satisfactory consistency was present among the trials, as exhibited in detail in Figure 3. The relevant funnel plots are shown in Figure 4, which demonstrated the satisfactory symmetry for publication bias in the pairwise comparisons. In terms of the transitivity analysis of the present NMA, only RCTs complying with the selection criterion were strictly included, ensuring the balance for cross-study transitivity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719650/fimmu-12-719650-HTML-r1/image_m/fimmu-12-719650-g003.jpg)
Figure 3 Local inconsistency of OS (A), PFS (B) and sAE (C) by “node-splitting” analyses. OS, overall survival; PFS, progression-free survival; sAE, severe acute events; SoC, standard-of-care therapy; ST, single targeted therapy different from SoC; DT, double targeted therapy; T+C, targeted therapy combined with chemotherapy; SI, single immune checkpoint therapy; DI, double immune checkpoint therapy; SC, single chemotherapy different from SoC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719650/fimmu-12-719650-HTML-r1/image_m/fimmu-12-719650-g004.jpg)
Figure 4 Funnel plots of OS (A), PFS (B) and sAE (C). OS, overall survival; PFS, progression-free survival; sAE, severe acute events; SoC, standard-of-care therapy; ST, single targeted therapy different from SoC; DT, double targeted therapy; T+C, targeted therapy combined with chemotherapy; SI, single immune checkpoint therapy; DI, double immune checkpoint therapy; SC, single chemotherapy different from SoC.
NMA for Efficacy and Safety
Considering treatment efficacy, 5184 and 4532 individuals involved in 23 and 22 studies were included in the analyses of OS and PFS, respectively; their network plots are displayed in Supplementary Figure 1. For OS, as shown in Table 1, significant superiority was observed for SI when compared with ST (HR 0.88, 95% CrI 0.66 to 0.99), SoC (HR 0.76, 95% CrI 0.64 to 0.89), DT (HR 0.75, 95% CrI 0.57 to 0.96) and SC (HR 0.52, 95% CrI 0.29 to 0.91). Except for SI and SC, all treatments (DI, ST, T+C, SoC, DT) tended to be consistent, with HRs close to 1 between groups. For PFS, as described in Table 1, ST tended to be superior to the other treatments, although only one pairwise comparison showed a significant difference (SoC vs ST: HR 1.20, 95% CrI 1.08 to 1.33).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719650/fimmu-12-719650-HTML-r1/image_m/fimmu-12-719650-t001.jpg)
Table 1 Estimate results according to the network meta-analysis on OS (lower left) and PFS (upper right), along with the ranking distribution by SUCRA values (in the below arrow shape).
In terms of safety profiles, 16 trials with 4026 participants were available for an analysis of overall sAE. The relevant network plot for sAE is shown in Supplementary Figure 1. Table 2 shows that there were significantly lower incidences of sAE with SI than with SoC (OR 0.31, 95% CrI 0.11 to 0.90), ST (OR 0.23, 95% CrI 0.06 to 0.86) and DT (OR 0.11, 95% CrI 0.02 to 0.53). However, pairwise comparisons between the other treatment models demonstrated no statistically significant differences.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719650/fimmu-12-719650-HTML-r1/image_m/fimmu-12-719650-t002.jpg)
Table 2 Estimate results according to the network meta-analysis on treatment-related sAE, along with the ranking distribution by SUCRA values (in the below arrow shape).
Rank Probabilities
The arrows below Tables 1 and 2 indicate the ranking distributions of safety and efficacy, respectively. Among all competing therapy models involved, as shown in Table 1, SI (SUCRA = 94.91%) was deemed to have the greatest probability of ranking first in terms of OS, followed by DI (SUCRA = 61.98%), ST (SUCRA = 58.08%), T+C (SUCRA = 55.69%), SoC (SUCRA = 36.34%), DT (SUCRA = 35.37%) and SC (SUCRA = 7.62%). The PFS ranking probability was not quite the same, even demonstrating outcomes opposite those of OS. The ranking sequence from best to worst for PFS were ST (SUCRA = 76.63%), T+C (SUCRA = 72.85%), SI (SUCRA = 56.37%), DI (SUCRA = 43.94%), SC (SUCRA = 37.20%), SoC (SUCRA = 36.71%) and DT (SUCRA = 26.30%). In Table 2, as sorted by the SUCRA values, the incidences of sAE ranked from low to high were SI (SUCRA = 89.61%), DI (SUCRA = 75.54%), SoC (SUCRA = 55.47%), SC (SUCRA = 42.41%), ST (SUCRA = 38.69%), T+C (SUCRA = 36.96%) and DT (SUCRA = 11.33%). SI performed best while DT was denoted worst.
Additionally, a cluster plot (Supplementary Figure 3) presenting both the efficacy and safety profiles that were defined by the SUCRA values for OS and sAE was plotted. It demonstrates the overall distribution in this setting and indicates that SI presented with the best OS and the optimal management of sAE.
Sensitivity Analysis
For trials with small sample sizes, the relevant conclusions are not as convincing as those for larger trials. Therefore, the sensitivity analysis in this section was performed by excluding studies in which any treatment arm included no more than 50 individuals (9, 11, 35, 38, 40, 43–45, 48, 50, 51). Then, the SUCRA values for OS, PFS and sAE in the remaining trials were generated, and the rank probabilities defined by the SUCRA values were subsequently obtained. Supplementary Table 3 denotes that SI still performed best in terms of OS and sAE, and ST was considered the optimal treatment model regarding PFS, which was consistent with the result from the initial NMA. These similar rank probabilities indicate the intrinsic robustness of the NMA to some extent, thus confirming the final results.
Discussion
This study is the first seeking to determine the optimal second-line systemic treatment model for R/M HNSCC using an NMA and revealed that SI had been shown to be the best second-line systemic treatment model in terms of OS and sAE compared with SoC, ST, DT, T+C, DI and SC based on literature retrieval up to June 2021, without any of the eligible RCTs applying immunosuppressant therapy in the first-line setting. Furthermore, the targeted agent-involved therapeutic modalities, including ST, DT and T+C, more frequently suffered from sAE, despite the superiority of ST and T+C with regard to PFS. These findings, to some extent, could substantially guide clinical options toward favorable regimens regarding both efficacy and toxicity in the second-line treatment of R/M HNSCC patients who are unable or do not have access to immunosuppressants in the first-line setting and pave the way for future research.
For patients who have progressive disease after first-line treatment, the prognosis is poor. SoC has long been established as a cornerstone treatment (35–37, 39, 41) and was used as a control arm in many more clinical trials than other treatments, including targeted therapy (9, 10), immunotherapy (13, 14, 16), and a combination of targeted therapy and chemotherapy (45, 49). This NMA revealed that SI was the best choice for OS, while ST was best for PFS. It was not until the emergence of the CheckMate 141 and KEYNOTE 040 trials (13, 14), which compared the administration of SI (nivolumab, pembrolizumab, respectively) with SoC as a second-line treatment for R/M HNSCC, that the OS could be significantly prolonged (7.5 months vs 5.1 months and 8.4 months vs 6.9 months, respectively; both p < 0.05). The present NMA further confirmed the superiority of SI on OS as a second-line treatment modality.
However, the EAGLE trial, which explored the application of durvalumab and tremelimumab, proved no significant OS advantage over SoC (durvalumab vs SoC: HR 0.88, 95% CrI 0.72 to 1.08, p = 0.20; durvalumab plus tremelimumab vs SoC: HR 1.04, 95% CrI 0.85 to 1.26, p = 0.76) (16). The underlying reasons why SoC provided a longer-than-expected OS may be attributed to the differences in the choice of SoC regimen, the subsequent regimen and the imbalanced baseline clinical characteristics between groups, confounding the finding that immunotherapy resulted in non-superior survival over SoC. Nevertheless, an OS benefit was indeed observed in the durvalumab arm, with a median OS of 9.8 months in the group with programmed cell death ligand 1 (PD-L1) expression in tumor cells (TCs) > 25% versus 7.6 months for PD-L1 expression in TCs < 25%. Similarly, for pembrolizumab in KEYNOTE 040, the patients were stratified by PD-L1 combined positive score (CPS) and tumor proportion score (TPS), and compared with SoC, a significant advantage was evident in patients with PD-L1 CPS ≥ 1% (8.7 months vs 7.1 months, p = 0.0049) and in those with PD-L1 TPS ≥ 50% (11.6 months vs 6.6 months, p = 0.0014) (14). In the application of nivolumab in the Check 141 trial (13), however, the subsequent follow-up and subgroup analyses revealed that a superior OS was continuously observed in the nivolumab group regardless of PD-L1 expression, human papilloma virus (HPV) status and prior cetuximab exposure (52, 53). Hence, the level of PD-L1 expression should be measured when certain ICIs, such as durvalumab and pembrolizumab, are prescribed in the second-line treatment of R/M HNSCC, whereas the use of nivolumab in this setting may be ignored as the present evidence indicates.
In terms of PFS, we found ST to be the best choice, while DT was inferior to the other treatment models. Of note, the targeted agent-based treatment regimens, namely, ST and T+C, were the top two choices in the prolonging of PFS, but disappointingly, PFS superiority did not convert into an OS benefit but notably increased the incidence of sAE. In regard to DT, the acute toxicities were too aggressive to be tolerated, with the incidence of sAE reaching more than 50% (42), which may account for its ranking worst in terms of PFS in the present work. Notably, SI and DI rank the top two in terms of OS while it’s not the case of PFS. The potential reason why PFS is not well correlated with OS when ICIs are introduced may be ascribed to the unique mechanisms to rehabilitate or activate self-immunity toward tumors, leading to delayed clinical effects and long-term survival benefits (54). The delayed effects of ICIs also arise owing to the phenomenon of disease progression (PD) followed by either tumor shrinkage (pseudoprogression) (55, 56), or a long post-PD survival. As such, PFS exhibited only moderate-to-poor correlation with OS for ICI trials (57, 58).
For the toxicity profiles, SI was revealed to be the safest choice in this NMA, while the targeted agent-based treatment regimens were correlated with higher incidences of sAE. Attention should be paid in regard to the quality of life of patients with R/M HNSCC, especially in the setting of second-line treatments. In clinical practice, first-line treatment tends to be intense and aggressive for patients with R/M HNSCC. In both the Extreme (59) and TPExtreme trials (60), cetuximab in combination with cytotoxic agents platinum plus 5-fluorouracil with or without docetaxel was the basis utilization, which indeed caused conspicuous side effects. In the Extreme trial (59), six cycles of platinum plus 5-fluorouracil with or without weekly cetuximab both resulted in profound ≥ grade 3 overall toxicity (82% and 76%, respectively), and the maintenance of cetuximab resulted in a certain proportion of ≥ grade 3 infusion-related adverse effects, such as septic shock (7%), hypomagnesemia (9%) and skin reactions (9%). Although TPExtreme (60) made efforts to improve patients’ treatment tolerance, the incidences of sAE were still far greater than envisaged, rendering subsequent second-line treatments somewhat troublesome. As recommended by the NCCN guidelines, the second-line choice for patients with R/M HNSCC includes mainly monotherapy, such as SI, SoC, SC and ST. However, attempts have been made to introduce combination methods such as T+C or DT (11, 38–40, 43–45, 49–51). The results demonstrated that neither T+C nor DT had an OS advantage over ST or SoC; instead, sAE were much more frequent in the combination arms. This NMA again emphasizes that SI showed overwhelming superiority in terms of tolerance over the other models and provides persuasive evidence from RCTs. The other monotherapy regimens, SC and ST, were similar in terms of sAE. The evidence presented here may guide directions for strategies when making second-line treatment decisions, among which SI monotherapy should be prioritized, while combination options, such as T+C or DT, should be considered with caution.
However, inevitably, some limitations should be acknowledged for this NMA. First, not all eligible trials simultaneously reported survival data and overall sAE, resulting in an outflow of participants in some trials. At the same time, except for SoC, indirect comparisons were made for most treatment models only, so the direct evidence among all models was insufficient, and it is difficult to form sufficient closed loops in the analyses of survival and sAE. In addition, although this NMA was strictly limited to RCTs, the randomization procedures were not thoroughly illustrated, making it difficult to accurately assess the risk of bias from the included studies. In view of the intrinsic limitations of meta analyses, this work was strictly abidied by the Preferred Items for Systematic Reviews and Meta-Analyses (PRISMA) reporting guidelines as well as the extension statement of NMA (21, 22), thus making all analyses more convincing.
In summary, SI was shown to be the best second-line systemic treatment model in terms of OS and sAE. Owing to the promising survival superiority that KEYNOTE 048 demonstrates (17), pembrolizumab was approved in the first-line setting of R/M HNSCC. Thus, the subsequence choices following first-line utilization of pembrolizumab will definitely catch public attention since preferring SI in the second-line setting will encounter new problems. There will be new clinical data challenging immunosuppressant therapy in the second-line again, which is an issue worthy of further exploration and being looked forward to in the near future.
Conclusion
This NMA, encompassing 5285 individuals from 24 trials, revealed that compared with SoC, ST, DT, T+C, DI and SC, SI achieved the best OS as well as the least sAE. SI may thus serve as the optimal second-line systemic treatment model for R/M HNSCC patients who are unable or do not have access to ICIs in the first-line setting.
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.
Author Contributions
Y-WY, R-HZ and T-ZY were guarantors of the entire study. Z-JZ, Y-WY, R-HZ and T-ZY built up the study concepts and designed this study. Z-JZ, FZ and W-YY collected the relevant data all needed. W-ZQ, K-L, J-HF, J-YT, and HL rechecked data and helped submit this manuscript. Z-JZ, W-YY and FZ performed the statistical analyses and wrote this paper. Y-WY, R-HZ and T-ZY revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from the National Natural Science Foundation of China (Grant number: 81773354) and the Key Clinical Technology of Guangzhou (Grant number: 2019ZD17).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
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/fimmu.2021.719650/full#supplementary-material
References
1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2018) 68:394–424. doi: 10.3322/caac.21492
2. Cancer Genome Atlas N. Comprehensive Genomic Characterization of Head and Neck Squamous Cell Carcinomas. Nature (2015) 517:576–82. doi: 10.1038/nature14129
3. Leemans CR, Tiwari R, Nauta JJP, Waal I, Snow GB. Regional Lymph Node Involvement and Its Significance in the Development of Distant Metastases in Head and Neck Carcinoma. Cancer (1993) 71:452–6. doi: 10.1002/1097-0142(19930115)71:2<452::AID-CNCR2820710228>3.0.CO;2-B
4. Ferlito A, Shaha AR, Silver CE, Rinaldo A, Mondin V. Incidence and Sites of Distant Metastases From Head and Neck Cancer. ORL (2001) 63:202–7. doi: 10.1159/000055740
5. León X, Quer M, Orús C, Venegas MP, López M. Distant Metastases in Head and Neck Cancer Patients Who Achieved Loco-Regional Control. Head Neck (2000) 22:680–6. doi: 10.1002/1097-0347(200010)22:7<680::AID-HED7>3.0.CO;2-J
6. Leoncini E, Ricciardi W, Cadoni G, Arzani D, Petrelli L, Paludetti G, et al. Adult Height and Head and Neck Cancer: A Pooled Analysis Within the INHANCE Consortium. Eur J Epidemiol (2014) 29:35–48. doi: 10.1007/s10654-013-9863-2
7. Seiwert TY, Cohen EEW. State-Of-the-Art Management of Locally Advanced Head and Neck Cancer. Br J Cancer (2005) 92:1341–8. doi: 10.1038/sj.bjc.6602510
8. Marur S, Forastiere AA. Head and Neck Squamous Cell Carcinoma: Update on Epidemiology, Diagnosis, and Treatment. Mayo Clin Proc (2016) 91:386–96. doi: 10.1016/j.mayocp.2015.12.017
9. Machiels JP, Haddad RI, Fayette J, Licitra LF, Tahara M, Vermorken JB, et al. Afatinib Versus Methotrexate as Second-Line Treatment in Patients With Recurrent or Metastatic Squamous-Cell Carcinoma of the Head and Neck Progressing on or After Platinum-Based Therapy (LUX-Head & Neck 1): An Open-Label, Randomised Phase 3 Trial. Lancet Oncol (2015) 16:583–94. doi: 10.1016/S1470-2045(15)70124-5
10. Guo Y, Ahn MJ, Chan A, Wang CH, Kang JH, Kim SB, et al. Afatinib Versus Methotrexate as Second-Line Treatment in Asian Patients With Recurrent or Metastatic Squamous Cell Carcinoma of the Head and Neck Progressing on or After Platinumbased Therapy (LUX-Head & Neck 3): An Open-Label, Randomised Phase III Trial. Ann Oncol (2019) 30:1831–9. doi: 10.1093/annonc/mdz388
11. Kochanny SE, Worden FP, Adkins DR, Lim DW, Bauman JE, Wagner SA, et al. A Randomized Phase 2 Network Trial of Tivantinib Plus Cetuximab Versus Cetuximab in Patients With Recurrent/Metastatic Head and Neck Squamous Cell Carcinoma. Cancer (2020) 126:2146–52. doi: 10.1002/cncr.32762
12. Farkona S, Diamandis EP, Blasutig IM. Cancer Immunotherapy: The Beginning of the End of Cancer? BMC Med (2016) 14:1–18. doi: 10.1186/s12916-016-0623-5
13. Ferris RL, Blumenschein G, Fayette J, Guigay J, Colevas AD, Licitra L, et al. Nivolumab for Recurrent Squamous-Cell Carcinoma of the Head and Neck. N Engl J Med (2016) 375:1856–67. doi: 10.1056/NEJMoa1602252
14. Cohen EEW, Soulières D, Le Tourneau C, Dinis J, Licitra L, Ahn M-J, et al. Pembrolizumab Versus Methotrexate, Docetaxel, or Cetuximab for Recurrent or Metastatic Head-and-Neck Squamous Cell Carcinoma (KEYNOTE-040): A Randomised, Open-Label, Phase 3 Study. Lancet (2019) 393:156–67. doi: 10.1016/S0140-6736(18)31999-8
15. Siu LL, Even C, Mesía R, Remenar E, Daste A, Delord JP, et al. Safety and Efficacy of Durvalumab With or Without Tremelimumab in Patients With PD-L1-Low/Negative Recurrent or Metastatic HNSCC: The Phase 2 CONDOR Randomized Clinical Trial. JAMA Oncol (2019) 5:195–203. doi: 10.1001/jamaoncol.2018.4628
16. Ferris RL, Haddad R, Even C, Tahara M, Dvorkin M, Ciuleanu TE, et al. Durvalumab With or Without Tremelimumab in Patients With Recurrent or Metastatic Head and Neck Squamous Cell Carcinoma: EAGLE, a Randomized, Open-Label Phase III Study. Ann Oncol (2020) 31:942–50. doi: 10.1016/j.annonc.2020.04.001
17. Burtness B, Harrington KJ, Greil R, Soulières D, Tahara M, de Castro G, et al. Pembrolizumab Alone or With Chemotherapy Versus Cetuximab With Chemotherapy for Recurrent or Metastatic Squamous Cell Carcinoma of the Head and Neck (KEYNOTE-048): A Randomised, Open-Label, Phase 3 Study. Lancet (2019) 394:1915–28. doi: 10.1016/S0140-6736(19)32591-7
18. Lang Y, Dong D, Wu B. Pembrolizumab vs the EXTREME Regimen in Recurrent or Metastatic Head and Neck Squamous Cell Carcinoma: A Cost-Effectiveness Analysis. Clin Drug Invest (2020) 40:1137–46. doi: 10.1007/s40261-020-00973-9
19. Zhou K, Li Y, Liao W, Zhang M, Bai L, Li Q. Pembrolizumab Alone or With Chemotherapy for Squamous Cell Carcinoma of the Head and Neck: A Cost-Effectiveness Analysis From Chinese Perspective. Oral Oncol (2020) 107:104754. doi: 10.1016/j.oraloncology.2020.104754
20. Ter Veer E, Van Oijen MGH, Van Laarhoven HWM. The Use of (Network) Meta-Analysis in Clinical Oncology. Front Oncol (2019) 9:822. doi: 10.3389/fonc.2019.00822
21. Moher D, Liberati A, Tetzlaff J, Altman DG, Group P. Preferred Reporting Items for Systematic Reviews and Meta-Analyses: The PRISMA Statement. BMJ (2009) 339:b2535. doi: 10.1136/bmj.b2535
22. Hutton B, Salanti G, Caldwell DM, Chaimani A, Schmid CH, Cameron C, et al. The PRISMA Extension Statement for Reporting of Systematic Reviews Incorporating Network Meta-Analyses of Health Care Interventions: Checklist and Explanations. Ann Intern Med (2015) 162:777–84. doi: 10.7326/M14-2385
23. Tierney JF, Stewart LA, Ghersi D, Burdett S, Sydes MR. Practical Methods for Incorporating Summary Time-to-Event Data Into Meta-Analysis. Trials (2007) 8:16. doi: 10.1186/1745-6215-8-16
24. Higgins JP, Altman DG, Gotzsche PC, Juni P, Moher D, Oxman AD, et al. The Cochrane Collaboration’s Tool for Assessing Risk of Bias in Randomised Trials. BMJ (2011) 343:d5928. doi: 10.1136/bmj.d5928
25. Tonin FS, Rotta I, Mendes AM, Pontarolo R. Network Meta-Analysis: A Technique to Gather Evidence From Direct and Indirect Comparisons. Pharm Pract (Granada) (2017) 15:943. doi: 10.18549/PharmPract.2017.01.943
26. Oravecz Z, Muth C. Fitting Growth Curve Models in the Bayesian Framework. Psychon Bull Rev (2018) 25:235–55. doi: 10.3758/s13423-017-1281-0
27. Neupane B, Richer D, Bonner AJ, Kibret T, Beyene J. Network Meta-Analysis Using R: A Review of Currently Available Automated Packages. PloS One (2014) 9:e115065. doi: 10.1371/journal.pone.0115065
28. Jansen JP, Fleurence R, Devine B, Itzler R, Barrett A, Hawkins N, et al. Interpreting Indirect Treatment Comparisons and Network Meta-Analysis for Health-Care Decision Making: Report of the ISPOR Task Force on Indirect Treatment Comparisons Good Research Practices: Part 1. Value Health (2011) 14:417–28. doi: 10.1016/j.jval.2011.04.002
29. Salanti G, Ades AE, Ioannidis JP. Graphical Methods and Numerical Summaries for Presenting Results From Multiple-Treatment Meta-Analysis: An Overview and Tutorial. J Clin Epidemiol (2011) 64:163–71. doi: 10.1016/j.jclinepi.2010.03.016
30. Sutton A, Ades AE, Cooper N, Abrams K. Use of Indirect and Mixed Treatment Comparisons for Technology Assessment. Pharmacoeconomics (2008) 26:753–67. doi: 10.2165/00019053-200826090-00006
31. Egger M, Smith GD, Schneider M, Minder C. Bias in Meta-Analysis Detected by a Simple, Graphical Test. BMJ (1997) 315:629–34. doi: 10.1136/bmj.315.7109.629
32. Higgins JPT, Thompson SG, Deeks JJ, Altman DG. Measuring Inconsistency in Meta-Analyses. BMJ (2003) 327:557–60. doi: 10.1136/bmj.327.7414.557
33. Van Valkenhoef G, Dias S, Ades AE, Welton NJ. Automated Generation of Node-Splitting Models for Assessment of Inconsistency in Network Meta-Analysis. Res Synth Methods (2016) 7:80–93. doi: 10.1002/jrsm.1167
34. Stang A. Critical Evaluation of the Newcastle-Ottawa Scale for the Assessment of the Quality of Nonrandomized Studies in Meta-Analyses. Eur J Epidemiol (2010) 25:603–5. doi: 10.1007/s10654-010-9491-z
35. Pivot X, Wadler S, Kelly C, Ruxer R, Tortochaux J, Stern J, et al. Result of Two Randomized Trials Comparing Nolatrexed (Thymitaq™) Versus Methotrexate in Patients With Recurrent Head and Neck Cancer. Ann Oncol (2001) 12:1595–9. doi: 10.1023/A:1013185402896
36. Stewart JSW, Cohen EEW, Licitra L, Van Herpen CML, Khorprasert C, Soulieres D, et al. Phase III Study of Gefitinib Compared With Intravenous Methotrexate for Recurrent Squamous Cell Carcinoma of the Head and Neck. J Clin Oncol (2009) 27:1864–71. doi: 10.1200/JCO.2008.17.0530
37. Machiels JP, Subramanian S, Ruzsa A, Repassy G, Lifirenko I, Flygare A, et al. Zalutumumab Plus Best Supportive Care Versus Best Supportive Care Alone in Patients With Recurrent or Metastatic Squamous-Cell Carcinoma of the Head and Neck After Failure of Platinum-Based Chemotherapy: An Open-Label, Randomised Phase 3 Trial. Lancet Oncol (2011) 12:333–43. doi: 10.1016/S1470-2045(11)70034-1
38. Gilbert J, Lee JW, Argiris A, Haigentz M, Feldman LE, Jang M, et al. Phase II 2-Arm Trial of the Proteasome Inhibitor, PS-341 (Bortezomib) in Combination With Irinotecan or PS-341 Alone Followed by the Addition of Irinotecan at Time of Progression in Patients With Locally Recurrent or Metastatic Squamous Cell Carcinoma of the Head and Neck (E1304): A Trial of the Eastern Cooperative Oncology Group. Head Neck (2013) 35:942–8. doi: 10.1002/hed.23046
39. Argiris A, Ghebremichael M, Gilbert J, Lee JW, Sachidanandam K, Kolesar JM, et al. Phase III Randomized, Placebo-Controlled Trial of Docetaxel With or Without Gefitinib in Recurrent or Metastatic Head and Neck Cancer: An Eastern Cooperative Oncology Group Trial. J Clin Oncol (2013) 31:1405–14. doi: 10.1200/JCO.2012.45.4272
40. Limaye S, Riley S, Zhao S, O’Neill A, Posner M, Adkins D, et al. A Randomized Phase II Study of Docetaxel With or Without Vandetanib in Recurrent or Metastatic Squamous Cell Carcinoma of Head and Neck (SCCHN). Oral Oncol (2013) 49:835–41. doi: 10.1016/j.oraloncology.2013.04.010
41. Seiwert TY, Fayette J, Cupissol D, del Campo JM, Clement PM, Hitt R, et al. A Randomized, Phase II Study of Afatinib Versus Cetuximab in Metastatic or Recurrent Squamous Cell Carcinoma of the Head and Neck. Ann Oncol (2014) 25:1813–20. doi: 10.1093/annonc/mdu216
42. Ruzsa A, Sen M, Evans M, Lee LW, Hideghety K, Rottey S, et al. Phase 2, Open-Label, 1:1 Randomized Controlled Trial Exploring the Efficacy of EMD 1201081 in Combination With Cetuximab in Second-Line Cetuximab-Naïve Patients With Recurrent or Metastatic Squamous Cell Carcinoma of the Head and Neck (R/M SCCHN). Invest New Drugs (2014) 32:1278–84. doi: 10.1007/s10637-014-0117-2
43. Jimeno A, Shirai K, Choi M, Laskin J, Kochenderfer M, Spira A, et al. A Randomized, Phase II Trial of Cetuximab With or Without PX-866, an Irreversible Oral Phosphatidylinositol 3-Kinase Inhibitor, in Patients With Relapsed or Metastatic Head and Neck Squamous Cell Cancer. Ann Oncol (2015) 26:556–61. doi: 10.1093/annonc/mdu574
44. Gilbert J, Schell MJ, Zhao X, Murphy B, Tanvetyanon T, Leon ME, et al. A Randomized Phase II Efficacy and Correlative Studies of Cetuximab With or Without Sorafenib in Recurrent and/or Metastatic Head and Neck Squamous Cell Carcinoma. Oral Oncol (2015) 51:376–82. doi: 10.1016/j.oraloncology.2014.12.011
45. Jimeno A, Bauman JE, Weissman C, Adkins D, Schnadig I, Beauregard P, et al. A Randomized, Phase 2 Trial of Docetaxel With or Without PX-866, an Irreversible Oral Phosphatidylinositol 3-Kinase Inhibitor, in Patients With Relapsed or Metastatic Head and Neck Squamous Cell Cancer. Oral Oncol (2015) 51:383–8. doi: 10.1016/j.oraloncology.2014.12.013
46. Machiels JP, Van Maanen A, Vandenbulcke JM, Filleul B, Seront E, Henry S, et al. Randomized Phase II Study of Cabazitaxel Versus Methotrexate in Patients With Recurrent and/or Metastatic Squamous Cell Carcinoma of the Head and Neck Previously Treated With Platinum-Based Therapy. Oncologist (2016) 21(12):1416–e17. doi: 10.1634/theoncologist.2016-0296
47. Fayette J, Wirth L, Oprean C, Udrea A, Jimeno A, Rischin D, et al. Randomized Phase II Study of Duligotuzumab (MEHD7945A) vs. Cetuximab in Squamous Cell Carcinoma of the Head and Neck (MEHGAN Study). Front Oncol (2016) 6:232. doi: 10.3389/fonc.2016.00232
48. Joshi A, Patil V, Noronha V, Dhumal S, Pande N, Chandrasekharan A, et al. Results of a Phase II Randomized Controlled Clinical Trial Comparing Efficacy of Cabazitaxel Versus Docetaxel as Second Line or Above Therapy in Recurrent Head and Neck Cancer. Oral Oncol (2017) 75:54–60. doi: 10.1016/j.oraloncology.2017.10.018
49. Soulières D, Faivre S, Mesía R, Remenár É, Li S-H, Karpenko A, et al. Buparlisib and Paclitaxel in Patients With Platinum-Pretreated Recurrent or Metastatic Squamous Cell Carcinoma of the Head and Neck (BERIL-1): A Randomised, Double-Blind, Placebo-Controlled Phase 2 Trial. Lancet Oncol (2017) 18:323–35. doi: 10.1016/S1470-2045(17)30064-5
50. Ferrarotto R, William WN, Tseng JE, Marur S, Shin DM, Murphy B, et al. Randomized Phase II Trial of Cixutumumab Alone or With Cetuximab for Refractory Recurrent/Metastatic Head and Neck Squamous Cell Carcinoma. Oral Oncol (2018) 82:83–90. doi: 10.1016/j.oraloncology.2018.05.014
51. Seiwert TY, Kochanny S, Wood K, Worden FP, Adkins D, Wade JL, et al. A Randomized Phase 2 Study of Temsirolimus and Cetuximab Versus Temsirolimus Alone in Recurrent/Metastatic, Cetuximab-Resistant Head and Neck Cancer: The MAESTRO Study. Cancer (2020) 126:3237–43. doi: 10.1002/cncr.32929
52. Ferris RL, Blumenschein G, Fayette J, Guigay J, Colevas AD, Licitra L, et al. Nivolumab vs Investigator’s Choice in Recurrent or Metastatic Squamous Cell Carcinoma of the Head and Neck: 2-Year Long-Term Survival Update of CheckMate 141 With Analyses by Tumor PD-L1 Expression. Oral Oncol (2018) 81:45–51. doi: 10.1016/j.oraloncology.2018.04.008
53. Ferris RL, Licitra L, Fayette J, Even C, Blumenschein G, Harrington KJ, et al. Nivolumab in Patients With Recurrent or Metastatic Squamous Cell Carcinoma of the Head and Neck: Efficacy and Safety in CheckMate 141 by Prior Cetuximab Use. Clin Cancer Res (2019) 25:5221–30. doi: 10.1158/1078-0432.CCR-18-3944
54. Constantinidou A, Alifieris C, Trafalis DT. Targeting Programmed Cell Death -1 (PD-1) and Ligand (PD-L1): A New Era in Cancer Active Immunotherapy. Pharmacol Ther (2019) 194:84–106. doi: 10.1016/j.pharmthera.2018.09.008
55. Abdel-Rahman O. Nonconventional Patterns of Benefit of Solid Tumors Treated With PD-(L)1 Inhibitors: A Systematic Review. Immunotherapy (2017) 9:995–1004. doi: 10.2217/imt-2017-0074
56. Long GV, Weber JS, Larkin J, Atkinson V, Grob JJ, Schadendorf D, et al. Nivolumab for Patients With Advanced Melanoma Treated Beyond Progression: Analysis of 2 Phase 3 Clinical Trials. JAMA Oncol (2017) 3:1511–9. doi: 10.1001/jamaoncol.2017.1588
57. Mushti SL, Mulkey F, Sridhara R. Evaluation of Overall Response Rate and Progression-Free Survival as Potential Surrogate Endpoints for Overall Survival in Immunotherapy Trials. Clin Cancer Res (2018) 24:2268–75. doi: 10.1158/1078-0432.CCR-17-1902
58. Nie RC, Chen FP, Yuan SQ, Luo YS, Chen S, Chen YM, et al. Evaluation of Objective Response, Disease Control and Progression-Free Survival as Surrogate End-Points for Overall Survival in Anti-Programmed Death-1 and Anti-Programmed Death Ligand 1 Trials. Eur J Cancer (2019) 106:1–11. doi: 10.1016/j.ejca.2018.10.011
59. Vermorken JB, Mesia R, Rivera F, Remenar E, Kawecki A, Rottey S, et al. Platinum-Based Chemotherapy Plus Cetuximab in Head and Neck Cancer. N Engl J Med (2008) 359:1116–27. doi: 10.1056/NEJMoa0802656
60. Guigay J, Fayette J, Dillies AF, Sire C, Kerger JN, Tennevet I, et al. Cetuximab, Docetaxel and Cisplatin as First-Line Treatment in Patients With Recurrent or Metastatic Head and Neck Squamous Cell Carcinoma: A Multicentre, Phase II GORTEC Study. Ann Oncol (2015) 26:1941–7. doi: 10.1093/annonc/mdv268
Keywords: recurrent and/or metastatic head and neck squamous cell carcinoma, treatment, efficacy, toxicity, network meta-analysis
Citation: Zhan Z-J, Yao W-Y, Zhang F, Qiu W-Z, Liao K-, Feng J-H, Tan J-Y, Liu H, Yuan T-Z, Zheng R-H and Yuan Y-W (2021) The Optimal Second-Line Systemic Treatment Model for Recurrent and/or Metastatic Head and Neck Squamous Cell Carcinoma: A Bayesian Network Meta-Analysis. Front. Immunol. 12:719650. doi: 10.3389/fimmu.2021.719650
Received: 02 June 2021; Accepted: 20 July 2021;
Published: 02 August 2021.
Edited by:
Linlang Go, Southern Medical University, ChinaReviewed by:
Siddharth Sheth, University of North Carolina at Chapel Hill, United StatesBenhua Xu, Fujian Medical University Union Hospital, China
Copyright © 2021 Zhan, Yao, Zhang, Qiu, Liao, Feng, Tan, Liu, Yuan, Zheng and Yuan. 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: Ya-Wei Yuan, yuanyawei@gzhmu.edu.cn; Rong-Hui Zheng, zhengronghui@gzhmu.edu.cn; Tai-Ze Yuan, taize.yuan@ccm.cn
†These authors have contributed equally to this work and share first authorship