- 1Department of Orthopedics, The Fourth Affiliated Hospital of China Medical University, Shenyang, China
- 2Center for Translational Medicine, The Fourth Affiliated Hospital of China Medical University, Shenyang, China
Background: Programmed cell death 1 ligand-1 (PD-L1) is an immune checkpoint molecule that acts to protect cancer cells from immune surveillance and is considered as a prognostic biomarker in several cancers, but the prognostic value of PD-L1 in bone and soft tissue sarcomas remains inconclusive. In the present meta-analysis, the clinicopathological and prognostic value of PD-L1 in sarcomas was evaluated.
Method: We performed a systemic and comprehensive meta-analysis by searching the PubMed, Medline, Cochrane Library, EMBASE, and Web of Science databases up to October 31, 2019. Eligible articles were incorporated, and pooled hazard ratios (HRs) and odds ratios (ORs) with their 95% confidence intervals (CIs) were used to estimate the outcomes.
Results: Thirty-six articles containing 39 independent studies with 3,680 bone and soft tissue sarcoma patients were included in our meta-analysis. The pooled results showed that PD-L1 overexpression could predict poor overall survival (HR 1.45, 95% CI 1.11–1.90, P < 0.01), metastasis-free survival (HR 1.58, 95% CI 1.14–2.19, P < 0.01), and event-free survival (HR 2.82, 95% CI 1.69–4.71, P < 0.01) in sarcomas. Furthermore, PD-L1 overexpression was correlated with a higher rate of tumor metastasis (OR 2.95, 95% CI 1.32–6.60, P < 0.01), a more advanced tumor grade (OR 3.63, 95% CI 2.55–5.16, P < 0.01), and more T lymphocyte infiltration (OR 5.55, 95% CI 2.86–10.76, P < 0.01). No obvious publication bias was observed, and the sensitivity analysis showed that our results were robust.
Conclusion: The results of our meta-analysis indicate that high PD-L1 expression might serve as a valuable and predictive biomarker for adverse clinicopathological features and poor prognosis in patients with sarcoma.
Introduction
Bone and soft tissue sarcomas are a group of heterogeneous and rare neoplasms with more than 50 distinct histologic subtypes that originate from mesenchymal tissues, and these sarcomas account for <1% of adult cancers and 7–15% of pediatric malignancies (1). In recent years, due to the improvement of surgical techniques and the effective combining of chemotherapy or radiotherapy with surgery, a satisfactory but stagnant survival of over 60% has been achieved (2–4). Moreover, the response of sarcomas to conventional chemotherapy and radiotherapy varies, and for those patients with recurrent and metastatic disease, their survival remains unsatisfactory, with an average survival period of <1 year (5–7). Thus, novel treatment options for these patients are of utmost urgency.
Immunotherapy is considered as one of the foremost options for cancer treatment and has recently achieved tremendous success (8, 9). In the immune system, immune tolerance is an important process that inhibits the role of innate immunity in the cancer microenvironment. One of the immune tolerance mechanisms is the immune checkpoint mechanism. Immune checkpoint molecules, such as cytotoxic T lymphocyte antigen 4, lymphocyte activation gene 3, programmed cell death 1 (PD-1), and programmed cell death 1 ligand-1 (PD-L1), are expressed in the microenvironment of cancers (10), and suppressive signals are transmitted to T cells to prevent excessive immune responses. Thus, immune checkpoint inhibitors are emerging and becoming a prevailing form of systemic therapy regarded to “release the brakes” on the immune system. In malignant tumors, PD-L1 is often expressed on the surface of tumor cells and helps promote tumor evasion from the immune system by activating PD-1 and inhibiting the function of T cells. Immune checkpoint inhibitors targeting the PD-L1/PD-1 immune checkpoint axis have shown considerable success in cancers (11) but not in bone and soft tissue sarcomas (12–14). To investigate the underlying mechanism of PD-L1 in the progression of sarcomas and assess the clinical value of PD-L1, the correlation between PD-L1 expression and the clinical outcomes of sarcoma patients was assessed. However, more studies are emerging to show controversial results in the past several years (15–17). In addition, although several meta-analyses have been conducted to assess the prognostic value of PD-L1 in sarcomas, the results have been inconclusive (18, 19) due to the inconsistent results of PD-L1 assessment assays and the limited number of included studies.
In the present meta-analysis, the prognostic significance of PD-L1 expression was systemically and comprehensively evaluated with 39 independent studies of 3,680 bone and soft tissue sarcoma patients to investigate whether PD-L1 could serve as a prognostic biomarker in sarcomas.
Materials and Methods
Search Strategies
A systemic and comprehensive literature search was conducted in five databases with no language or publication date limitations: PubMed, Medline, Cochrane Library, EMBASE, and Web of Science. Articles published before October 2019 were included in the meta-analysis, and the last search was conducted on October 31, 2019. The following keywords were used to search the relevant literature: “sarcoma” or “bone sarcoma” or “soft tissue sarcoma” or “osteosarcoma” or “Ewing sarcoma” or “chondrosarcoma” or “myosarcoma” or “fibrosarcoma” or “synovial sarcoma” or “malignant fibrous histiocytoma” or “liposarcoma” or “epithelioid sarcoma” or “spindle cell sarcoma” or “carcinosarcoma” or “leiomyosarcoma” or “angiosarcoma” or “hemangiosarcoma” or “lymphangiosarcoma” or “malignant lymphangioendothelioma” and “PD-L1” or “B7-H1” or “CD274” or “programmed cell death 1 ligand 1 protein.” Moreover, we manually searched the references in each study to avoid missing potentially relevant data.
Inclusion and Exclusion Criteria
The eligible studies included in the meta-analysis had to be in accordance with the following criteria: (1) studies on patients with pathologically confirmed bone and/or soft tissue sarcomas, (2) studies focusing on the correlations between PD-L1 expression and the survival and/or clinicopathological outcome(s), and (3) studies utilizing defined cutoff values to classify the PD-L1 expression as “high” and “low” or “positive” and “negative.” Studies were excluded if they (1) were reviews, case reports, letters, or conference abstracts, (2) investigated cell experiments and/or animal experiments, (3) were not related to PD-L1 expression, or (4) comprised overlapping patients and/or insufficient data. The eligibility of the included studies was determined by two authors independently. Any discrepancy was resolved by a consensus after a discussion.
Data Extraction and Quality Assessment
The data of interest were extracted, including (1) basic characteristics (first author, publication year, tumor type, number and source of patients, presence of positive/high PD-L1 expression, cutoff value, assay method, biomarkers, and follow-up duration), (2) hazard ratios (HRs) with 95% confidence intervals (CIs) that either were given directly or could be calculated from published and available data by using the Cox regression method with a standard procedure (SPSS version 16.0) or Tierney's methods (20), and (3) data to estimate odds ratios (ORs) with 95% CIs for the correlations between PD-L1 and clinicopathological outcome(s). The HRs from the multivariate analysis were extracted if the HRs from both the univariate and the multivariate analysis were available.
The Newcastle–Ottawa Scale (NOS, https://www.ohri.ca/programs/clinical/epidemiology/oxford.asp) system was used to assess the methodological quality of each study, and the included studies were assessed for risk of bias by using the Quality and Prognosis Studies (QUIPS) method (21). The risk of bias was assessed for the following domains: study participation, study attrition, prognostic factor measurement, outcome measurement, study confounding, and statistical analysis. The score was classified and reported as low, moderate, or high. The assessment was performed by two authors independently. If there was any disagreement, it was resolved by a consensus after a discussion.
Statistical Analysis
Estimates of HRs, ORs, and their 95% CIs were used to assess the pooled survival related to and the clinicopathological significance of PD-L1 expression. Cochran's Q test and Higgins I2 statistic were employed to assess the statistical heterogeneity among studies. A fixed-effects model was appropriately used if there was no significant heterogeneity (P > 0.10 or I2 <50%); otherwise, a random-effects model was used, and a subgroup meta-analysis was conducted to identify the underlying heterogeneity when at least five studies were included. Publication bias was assessed by using a funnel plot analysis, Begg's tests, and Egger's tests, and a sensitivity analysis was used to assess the stability of the pooled results when at least five studies were included. P < 0.05 was considered as statistically significant. All analyses were conducted by using Stata version 14.0 (Stata Corporation, College Station, TX, USA).
Results
Study Selection
There were 734 articles identified initially from the systemic and comprehensive search. After the duplicates were excluded, 531 articles were subjected to further screening. After viewing the title and the abstract, 378 articles remained. Three hundred thirty articles were omitted due to irrelevance. After a full-text review, 12 articles were excluded, including one with overlapping data and 11 with insufficient data. In total, 36 articles containing 39 independent studies with 3,680 bone and soft tissue sarcoma patients were included (15–17, 22–54), among which 31 articles were used to evaluate the prognostic significance of PD-L1 expression and 27 articles were used to evaluate the clinicopathological significance. The selection process is shown in Figure 1.
Study Characteristics and Quality Assessment
The summarized study characteristics are shown in Table 1. Briefly, the publication year of the eligible studies, written in English, ranged from 2013 to 2019. The patients were collected internationally from Asian (n = 17) to non-Asian (n = 23) populations. The number of patients ranged from 11 to 492. A variety of sarcomas were investigated in the studies, including osteosarcoma (n = 11), soft tissue sarcoma (STS, n = 9), unclassified sarcoma (n = 4), angiosarcoma (n = 4), leiomyosarcoma (n = 3), liposarcoma (n = 2), chondrosarcoma (n = 2), Ewing's sarcoma (n = 1), myeloid sarcoma (n = 1), pleomorphic sarcoma (n = 1), and synovial sarcoma (n = 1). Immunohistochemistry (IHC), quantitative polymerase chain reaction (qPCR), or RNA sequencing methodologies were employed to evaluate the PD-L1 expression. The cutoff values were available in studies using IHC, and the positivity ranged from 5.90 to 75.90%. Moreover, data from studies were extracted to analyze the HRs of overall survival (OS, n = 30), as a primary endpoint, and other survival outcomes, including recurrence-free survival (RFS, n = 4), metastasis-free survival (MFS, n = 3), event-free survival (EFS, n = 3), disease-free survival (DFS, n = 5), and progression-free survival (PFS, n = 4). The follow-up duration ranged from 1 to 366 months based on the available data. The correlation between PD-L1 expression and clinicopathological outcomes was assessed with data from 29 studies, including patient age and sex, tumor size, depth and localization, tumor grade and stage, metastasis, recurrence, previous chemotherapy/radiotherapy, chemotherapy response, tumor-infiltrating lymphocytes (TILs), and programmed cell death 1 (PD1) positivity.
The included 39 studies were critically appraised for methodological quality and risk of bias by using the NOS system and the QUIPS criteria, respectively. The results showed that all of the studies were of high methodological quality, with scores that ranged from six to nine stars. Most of the studies were retrospective studies and were assessed as having a low or moderate risk of bias. Almost no study provided information about patients lost to follow-up. Due to a lack of sufficient data on survival outcome(s) or study confounding, studies with only clinicopathological features were scored with a high risk of bias. Some studies were scored with a high or moderate risk of study participation due to missing details and the limited number of patients. A few studies were scored with a high or moderate risk of study confounding because of insufficient survival analyses for some important factors. The results are shown in Supplementary Table 1.
Prognostic Value of PD-L1 Expression in Sarcomas
A random-effects model was used in the analysis for OS as obvious heterogeneity was observed (I2 = 64.10%, P = 0.00). The results showed that elevated PD-L1 expression was correlated with poor OS (HR 1.45, 95% CI 1.11–1.90, P < 0.01, Figure 2A). Moreover, to investigate the underlying source of heterogeneity, subgroup analyses by tumor type, publication date, sample size and source, analysis model of effects, PD-L1 positivity, assay method, and PD-L1 expression region (in the tumor or the tumor microenvironment) were conducted. The results showed that none of these factors were the underlying source of heterogeneity, and the overall value of PD-L1 did not change (Table 2). Moreover, as various tumor subtypes were observed in the present analysis, a stratified analysis by tumor subtype was conducted, in which no less than three studies were included. A random-effects model was used (I2 = 64.60%, P = 0.00), and the results showed that PD-L1 overexpression could predict poor OS specifically in osteosarcoma (P < 0.01) and leiomyosarcoma (P < 0.05, Figure 2B). Thus, it was indicated that PD-L1 was a predictor for poor OS in sarcoma patients.
Figure 2. Forest plots of pooled HR for pooled overall survival (OS) (A) and OS stratified by tumor subtypes (B), metastasis-free survival (C), and event-free survival (D).
Other extracted important survival outcomes in the present analysis, including RFS, MFS, PFS, EFS, and DFS, are survival parameters related to disease according to their definition, but these data had limited numbers and thus insufficient effect size. In our analysis, the results showed that an elevated PD-L1 expression predicted poor MFS (HR 1.58, 95% CI 1.14–2.19, P < 0.01, Figure 2C) and EFS (HR 2.82, 95% CI 1.69–4.71, P < 0.01, Figure 2D) but was not associated with other survival outcomes (Table 2). As obvious heterogeneity was observed (I2 = 74.50%, P = 0.003) in the DFS analysis, a random-effects model was used, and subgroup analyses were conducted. The results of the subgroup analyses showed that the heterogeneity was not reduced by stratification of the data by publication date, tumor type, sample size, PD-L1 positivity, cutoff value, and PD-L1 expression region. However, the pooled HR did not change, indicating that the result was stable. In addition, although stratified analyses by tumor type (bone sarcoma vs. non-bone sarcoma) were conducted for other survival outcomes, no conclusive and reliable results were observed due to the limited number of included studies, as shown in Supplementary Table 2.
Clinicopathological Value of PD-L1 Expression in Sarcomas
The clinicopathological value of PD-L1 expression is summarized in Table 3. Briefly, 29 studies were included, and comprehensive meta-analyses were conducted to determine the correlation between the PD-L1 expression and patient age (n = 4) and sex (n = 21), tumor size (n = 6), depth (n = 3) and localization (n = 8), tumor grade (n = 15) and stage (n = 8), metastasis (n = 11), recurrence (n = 7), previous chemotherapy/radiotherapy (n = 10 and n = 7, respectively), chemotherapy response (n = 5), and TILs (n = 12). A random-effects model or a fixed-effects model was used based on the results of the heterogeneity analysis for each parameter, as shown in Table 3. PD-L1 overexpression was positively correlated with a higher rate of metastasis (OR 2.95, 95% CI 1.32–6.60, P < 0.01), a higher tumor grade (OR 3.63, 95% CI 2.55–5.16, P < 0.01), and more TILs (OR 5.55, 95% CI 2.86–10.76, P < 0.01) but was not correlated with other clinicopathological outcomes. For the TIL subset analysis, the PD-L1 expression was not correlated with CD4-positive TILs (OR 2.45, 95% CI 0.69–8.73, P = 0.17). Moreover, a stratified analysis by bone sarcoma vs. non-bone sarcoma was conducted for each outcome to clarify whether the correlation was tumor type dependent. The results showed that PD-L1 was specifically correlated with metastasis and CD8+ TILs in bone sarcomas and CD3+ TILs in non-bone sarcomas. In addition, PD-L1 overexpression was correlated with a good response to chemotherapy in bone sarcomas but with a poor response in non-bone sarcomas (Supplementary Table 3). For the outcomes with obvious heterogeneity, subgroup analyses were conducted. The results showed that some factors might be the underlying source of heterogeneity in some of the clinicopathological outcomes, as shown in Supplementary Table 4.
Publication Bias and Sensitivity Analysis
Publication bias was assessed by using funnel plot analysis (Figure 3) and Begg's and Egger's tests (Table 4). For survival outcomes, no obvious publication bias was observed in the evaluation of OS (Figure 3A) or RFS, MFS, EFS, and DFS (Table 4), but publication bias was noted for PFS (P = 0.04 in Egger's tests). For clinicopathological outcomes, Egger's tests showed that there might be a publication bias in the analysis of PD1+ TILs (P = 0.03). No obvious publication bias was observed in any other analysis of clinicopathological outcomes.
A sensitivity analysis was performed to assess the effects of each individual study on the pooled HRs and ORs of our analysis. The results of the sensitivity analysis, which was performed by omitting one study at a time, showed that no unique study significantly affected the pooled HRs for OS (Figure 3B) and other outcomes (Supplementary Figures 1–14), indicating that the results of the present meta-analysis were stable and credible.
Discussion
PD-L1, one of the immune checkpoint molecules, can inhibit T cell activation and proliferation by binding with PD-1 and protect tumor cells from the host immunologic surveillance system. Immunotherapy, by administration of immune checkpoint inhibitors, has been considered as the foremost method in individualized medicine due to its tremendous success in cancer treatment. Clinical trials have shown satisfactory results by using anti-PD-1/PD-L1 blockade in some cancers (11). However, the results in bone and soft tissue sarcomas are unsatisfactory (14). Recent studies have indicated that the high expression of PD-L1 is associated with poor prognosis in some cancers (55, 56). Moreover, several meta-analyses have evaluated the clinical significance of PD-L1 expression in patients with sarcomas, but the results have been inconclusive (18, 19). In the present meta-analysis, 39 independent studies (including 3,680 patients) were included, and the pooled results showed that PD-L1 overexpression could predict poor OS, MFS, and EFS. Furthermore, an elevated PD-L1 expression was correlated with a higher tumor metastasis rate, a more advanced tumor grade, and more T lymphocyte infiltration.
In terms of survival outcomes, the results showed that PD-L1 overexpression could predict poor OS and EFS, which was similar to the findings of previous reports (18, 19). As obvious heterogeneity was observed for OS, a random-effects model was used to reduce heterogeneity, and subgroup analyses were conducted. Although the results did not identify the underlying source of heterogeneity, the pooled results did not change in each subgroup analysis. Moreover, a stratified analysis by tumor subtype showed that PD-L1 overexpression could specifically predict poor OS in osteosarcoma and leiomyosarcoma. There was no obvious publication bias, and the sensitivity analysis indicated that the results were stable. Therefore, our results were robust and reliable. Some other factors should be considered to explain the underlying heterogeneity. There are more than 50 different subtypes of sarcoma, and more than 10 subtypes were included in the present analysis, which introduced heterogeneity into the whole meta-analysis. Moreover, differences in the detection techniques, such as antibody selection and dilution, and the antigen retrieval methods may have affected the sensitivity of detection, and the positivity of the PD-L1 expression might have been affected by the different cutoff values. In addition, the region used to define the PD-L1 expression was different among the included studies. These might be other explanations for the heterogeneity.
It was initially reported that PD-L1 overexpression predicted poor MFS but was not associated with DFS, RFS, or PFS in patients with sarcoma. A previous meta-analysis showed a similar result that the PD1/PD-L1 expression was associated with lymph node metastasis in patients with osteosarcoma (57). PD-L1 is widely considered to function in the tumor immunologic surveillance system; however, the tumor-intrinsic roles of PD-L1 in regulating the ability of human cancers to disseminate and to metastasize are currently not understood. PD-L1 is suggested to be involved in cancer metastasis via several mechanisms. PD-L1 is shown to promote epithelial-to-mesenchymal transition and is involved in tumor proliferation, migration, and invasion (58, 59). Moreover, cancer stem cells (CSCs) are considered to be involved in tumor metastasis (60), and PD-L1 is suggested to sustain stemness by promoting OCT4 and Nanog expression in CSCs (61) and promoting tumor-initiating cell generation and virulence in cancer (62). Thus, it is suggested that PD-L1 is involved in the possible immune evasion mechanism employed by CSCs during metastasis (63). In addition, several signaling pathways have been suggested to modulate the role of PD-L1 in cancer initiation, metastasis, and chemoresistance, including the MAPK/ERK pathway, the PI3K/AKT pathway, and the RAS/MEK pathway (64, 65). Sarcomas are a group of malignancies with a high propensity for metastasis, and the prognosis of metastatic patients remains poor. Further studies are encouraged to investigate the mechanism by which PD-L1 is involved in sarcoma metastasis. However, as there were only three studies in the pooled analysis for MFS, the results should be interpreted with caution.
The results from the analyses of clinicopathological outcomes were similar to those from previous reports (18, 19). An elevated PD-L1 expression was correlated with a higher rate of metastasis, a higher tumor grade, and more T lymphocyte infiltration. In the analysis for tumor grade, the result was different from a previous report (19). In other cancers, the correlation between PD-L1 expression and tumor grade remains controversial (55, 66), which might arise from the inherent differences between cancers. Moreover, the prognosis of patients with high-grade sarcomas is not satisfactory. It was suggested that PD-L1 overexpression might predict the presence of advanced sarcoma. In addition, our results might be more conclusive and reliable than those from other meta-analyses because more studies were included (n = 15). Moreover, although the results showed that PD-L1 expression was not correlated with chemotherapy response, stratified analyses by tumor type showed that PD-L1 overexpression might predict good response to chemotherapy in bone sarcomas but poor response in non-bone sarcomas, and this difference might be due to the difference in chemotherapy sensitivity between these two types of sarcoma.
Another major result in our analysis was that high PD-L1 expression was correlated with more TILs, which was similar to the findings of previous analyses in sarcomas and some other cancers (19, 67, 68). However, in the analyses of T cell subtype, the results showed that elevated PD-L1 expression was mainly correlated with PD1+ T cells, CD3+ T cells, and CD8+ T cells but not with CD4+ T cells. In the process of tumor immune tolerance, the PD1/PD-L1 axis acts to suppress the T cell response. PD-L1 is expressed on tumor cells and binds with PD-1, which is expressed on tumor-infiltrating CD8+ T cells as well as CD4+ T cells and other inflammatory cells. Additionally, PD-L1 makes tumor cells less susceptible to specific CD8+ T cells (69). Therefore, the killing effect of T cells on tumor cells is inhibited (70, 71). Moreover, PD-L1 expression can be modulated by tumor cells. One feasible pathway is induced by IFNγ production and subsequent IFNGR/JAK/STAT signaling in tumor cells, which depends on TILs (72, 73). Thus, the correlation between PD-L1 overexpression and increased tumor-infiltrating T cells in sarcomas might indicate a higher capacity of sarcomas to evade immune surveillance. However, inhibitors targeting the PD1/PD-L1 axis have not resulted in satisfactory results in sarcomas. The combination of immunotherapy with other therapeutics is postulated to show improved responses over checkpoint inhibitor monotherapies (74, 75), and investigations to explore the underlying mechanisms by which PD-L1 is involved in sarcoma are encouraged.
Several limitations in our analysis should be acknowledged. First, there might be some potential sources of bias. The HRs and their corresponding 95% CIs were extracted indirectly or directly from the studies, including from multivariate and univariate analyses, which might lead to statistical bias. Moreover, most studies included in this meta-analysis were retrospective, investigating patients with osteosarcoma and STS. A limited number of studies of other rare sarcomas were included, which might lead to subject selection bias. In addition, some other underlying sources of publication bias should be considered, such as unpublished negative results (76) or studies that could not be included due to language limitations or insufficient data (77). Second, due to the limited number of studies included for some outcomes, the results should be interpreted with caution. Third, a lack of uniformity in the definitions of some clinicopathological outcomes was observed, such as tumor stage, size and depth, and patient age, which might be due to the authors' preferences or the characteristics of certain sarcomas, and the results should be interpreted with caution. In addition, although the onset of sarcomas is age dependent and the immune context might play a different role in young patients than in elderly patients, we could not obtain a conclusive result regarding whether PD-L1 expression was correlated with patient age due to the limited number of studies and the lack of uniformity in age classification. In addition, PD-L1 expression has been considered to correlate with tumor response to anti-PD-1/PD-L1 agents in melanoma patients (65), but no study investigated the correlation between PD-L1 expression and administration of PD-L1 inhibitors, which may be a very important factor for evaluating their efficacy. Because sufficient data could not be extracted, the correlation between PD-L1 expression and PD1-positive tumors was not assessed. Sarcomas are a group of rare malignancies with low morbidity, and there were not sufficient numbers of studies for some tumor subtypes; therefore, more investigations are encouraged to clarify the significance of PD-L1 for those rare tumor subtypes.
Moreover, the assessment of PD-L1 expression through IHC is advocated by many studies, and a variety of antibodies and cutoff values for positive PD-L1 labeling are used. Other assay methods, such as qPCR and RNA sequencing, are also used to assess PD-L1 expression. To avoid selection bias and missing important data, we included studies using IHC and other assay methods. The results of the subgroup analysis by assay method showed that high PD-L1 expression predicted poor OS at both the protein level (HR 1.36, 95% CI 1.00–1.86, P < 0.01) and the mRNA level (HR 2.06, 95% CI 1.42–2.98, P < 0.01), and the assay method was not an underlying source of heterogeneity. As different cutoff values for PD-L1 positivity might affect the reliability of the evaluation, a subgroup analysis by cutoff value for IHC data was conducted. The results showed that PD-L1 overexpression could not predict poor OS with a defined cutoff value of ≥1% of cells, ≥5% of cells, or ≥10% of cells (Supplementary Figure 15), suggesting that the cutoff value could impact the prognostic value of PD-L1. It has been reported that different IHC methodologies for PD-L1 assessment provide different results (78); therefore, other approaches should be employed to improve the evaluation of PD-L1 expression in clinical practice (79), including in studies of patients with sarcomas.
Conclusions
In this meta-analysis, the results showed that PD-L1 overexpression could predict poor survival and was correlated with adverse tumor status in sarcoma patients. The results suggest that PD-L1 is a valuable prognostic biomarker in bone and soft tissue sarcomas, although more well-designed prospective studies with appropriate multivariate analyses are needed to validate our results.
Data Availability Statement
All datasets analyzed for this study are included in the article/Supplementary Material.
Author Contributions
FW and ZZ conceived and designed the project. FW, TY, and HZ analyzed and interpreted the data. FW, TY, and ZZ wrote the paper. All the authors have read and approved the final version of the manuscript and participated in acquiring the data.
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.00749/full#supplementary-material
References
1. Burningham Z, Hashibe M, Spector L, Schiffman JD. The epidemiology of sarcoma. Clin Sarcoma Res. (2012) 2:14. doi: 10.1186/2045-3329-2-14
2. Liu CY, Yen CC, Chen WM, Chen TH, Chen PC, Wu HT, et al. Soft tissue sarcoma of extremities: the prognostic significance of adequate surgical margins in primary operation and reoperation after recurrence. Ann Surg Oncol. (2010) 17:2102–11. doi: 10.1245/s10434-010-0997-0
3. Friebele JC, Peck J, Pan X, Abdel-Rasoul M, Mayerson JL. Osteosarcoma: a meta-analysis and review of the literature. Am J Orthop. (2015) 44:547–53.
4. Lahat G, Tuvin D, Wei C, Anaya DA, Bekele BN, Lazar AJ, et al. New perspectives for staging and prognosis in soft tissue sarcoma. Ann Surg Oncol. (2008) 15:2739–48. doi: 10.1245/s10434-008-9970-6
5. Italiano A, Mathoulin-Pelissier S, Cesne AL, Terrier P, Bonvalot S, Collin F, et al. Trends in survival for patients with metastatic soft-tissue sarcoma. Cancer. (2011) 117:1049–54. doi: 10.1002/cncr.25538
6. Hagleitner MM, de Bont ES, Te Loo DM. Survival trends and long-term toxicity in pediatric patients with osteosarcoma. Sarcoma. (2012) 2012:636405. doi: 10.1155/2012/636405
7. Tsukushi S, Nishida Y, Urakawa H, Kozawa E, Ishiguro N. Prognostic significance of histological invasion in high grade soft tissue sarcomas. Springerplus. (2014) 3:544. doi: 10.1186/2193-1801-3-544
8. Maccio A, Madeddu C. Blocking inflammation to improve immunotherapy of advanced cancer. Immunol. (2020) 159:357–64. doi: 10.1111/imm.13164
9. Taefehshokr N, Baradaran B, Baghbanzadeh A, Taefehshokr S. Promising approaches in cancer immunotherapy. Immunobiology. (2019) 225:151875. doi: 10.1016/j.imbio.2019.11.010
10. Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. (2012) 12:252–64. doi: 10.1038/nrc3239
11. Romano E, Romero P. The therapeutic promise of disrupting the PD-1/PD-L1 immune checkpoint in cancer: unleashing the CD8 T cell mediated anti-tumor activity results in significant, unprecedented clinical efficacy in various solid tumors. J Immunother Cancer. (2015) 3:15. doi: 10.1186/s40425-015-0059-z
12. Tamura K, Hasegawa K, Katsumata N, Matsumoto K, Mukai H, Takahashi S, et al. Efficacy and safety of nivolumab in Japanese patients with uterine cervical cancer, uterine corpus cancer, or soft tissue sarcoma: multicenter, open-label phase 2 trial. Cancer Sci. (2019) 110:2894–904. doi: 10.1111/cas.14148
13. Groisberg R, Hong DS, Behrang A, Hess K, Janku F, Piha-Paul S, et al. Characteristics and outcomes of patients with advanced sarcoma enrolled in early phase immunotherapy trials. J Immunother Cancer. (2017) 5:100. doi: 10.1186/s40425-017-0301-y
14. Chawla SP, van Tine BA, Pollack S, Ganjoo KN, Elias AD, Riedel RF, et al. A phase II randomized study of CMB305 and atezolizumab versus atezolizumab in NY-ESO-1+ soft tissue sarcoma: analysis of immunogenicity, tumor control, and patient survival. J Clin Oncol. (2019) 37(Suppl. 15):11011. doi: 10.1200/JCO.2019.37.15_suppl.11011
15. Bagaria SP, Gatalica Z, Maney T, Serie D, Parasramka M, Attia S, et al. Association between programmed death-ligand 1 expression and the vascular endothelial growth factor pathway in Angiosarcoma. Front Oncol. (2018) 8:71. doi: 10.3389/fonc.2018.00071
16. Boxberg M, Steiger K, Lenze U, Rechl H, von Eisenhart-Rothe R, Wortler K, et al. PD-L1 and PD-1 and characterization of tumor-infiltrating lymphocytes in high grade sarcomas of soft tissue - prognostic implications and rationale for immunotherapy. Oncoimmunology. (2018) 7:e1389366. doi: 10.1080/2162402X.2017.1389366
17. Kawamura A, Kawamura T, Riddell M, Hikita T, Yanagi T. Regulation of programmed cell death ligand 1 expression by atypical protein kinase C lambda/iota in cutaneous angiosarcoma. Cancer Sci. (2019) 110:1780–9. doi: 10.1111/cas.13981
18. Zheng C, You W, Wan P, Jiang X, Chen J, Zheng Y, et al. Clinicopathological and prognostic significance of PD-L1 expression in sarcoma: a systematic review and meta-analysis. Medicine. (2018) 97:e11004. doi: 10.1097/MD.0000000000011004
19. Zhu Z, Jin Z, Zhang M, Tang Y, Yang G, Yuan X, et al. Prognostic value of programmed death-ligand 1 in sarcoma: a meta-analysis. Oncotarget. (2017) 8:59570–80. doi: 10.18632/oncotarget.19168
20. 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
21. Hayden JA, Van Der Windt DA, Cartwright JL, Côté P, Bombardier C. Assessing bias in studies of prognostic factors. Ann Intern Med. (2013) 158:280–6. doi: 10.7326/0003-4819-158-4-201302190-00009
22. Yoshida K, Okamoto M, Sasaki J, Kuroda C, Ishida H, Ueda K, et al. Clinical outcome of osteosarcoma and its correlation with programmed death-ligand 1 and T cell activation markers. Onco Targets Ther. (2019) 12:2513–8. doi: 10.2147/OTT.S198421
23. Yan L, Wang Z, Cui C, Guan X, Dong B, Zhao M, et al. Comprehensive immune characterization and T-cell receptor repertoire heterogeneity of retroperitoneal liposarcoma. Cancer Sci. (2019) 110:3038–48. doi: 10.1111/cas.14161
24. Zheng B, Ren T, Huang Y, Sun K, Wang S, Bao X, et al. PD-1 axis expression in musculoskeletal tumors and antitumor effect of nivolumab in osteosarcoma model of humanized mouse. J Hematol Oncol. (2018) 11:16. doi: 10.1186/s13045-018-0580-x
25. Yang Z, Zhu G, Yang X, Zeng K, Liu F, Sun J. Expression of PD-L1 associated with Ki-67 and chemotherapy response but not p53 in osteosarcoma. Int J Clin Exp Med. (2018) 11:7.
26. Yang X, Zhu G, Yang Z, Zeng K, Liu F, Sun J. Expression of PD-L1/PD-L2 is associated with high proliferation index of Ki-67 but not with TP53 overexpression in chondrosarcoma. Int J Biol Markers. (2018) 33:507–13. doi: 10.1177/1724600818774464
27. Patel KR, Martinez A, Stahl JM, Logan SJ, Perricone AJ, Ferris MJ, et al. Increase in PD-L1 expression after pre-operative radiotherapy for soft tissue sarcoma. Oncoimmunology. (2018) 7:e1442168. doi: 10.1080/2162402X.2018.1442168
28. Park HK, Kim M, Sung M, Lee SE, Kim YJ, Choi YL. Status of programmed death-ligand 1 expression in sarcomas. J Transl Med. (2018) 16:303. doi: 10.1186/s12967-018-1658-5
29. Oike N, Kawashima H, Ogose A, Hotta T, Hatano H, Ariizumi T, et al. Prognostic impact of the tumor immune microenvironment in synovial sarcoma. Cancer Sci. (2018) 109:3043–54. doi: 10.1111/cas.13769
30. Machado I, Lopez-Guerrero JA, Scotlandi K, Picci P, Llombart-Bosch A. Immunohistochemical analysis and prognostic significance of PD-L1, PD-1, and CD8+ tumor-infiltrating lymphocytes in Ewing's sarcoma family of tumors (ESFT). Virchows Arch. (2018) 472:815–24. doi: 10.1007/s00428-018-2316-2
31. Kostine M, Briaire-de Bruijn IH, Cleven AHG, Vervat C, Corver WE, Schilham MW, et al. Increased infiltration of M2-macrophages, T-cells and PD-L1 expression in high grade leiomyosarcomas supports immunotherapeutic strategies. Oncoimmunology. (2018) 7:e1386828. doi: 10.1080/2162402X.2017.1386828
32. Kawamoto K, Miyoshi H, Suzuki T, Kiyasu J, Yokoyama S, Sasaki Y, et al. Expression of programmed death ligand 1 is associated with poor prognosis in myeloid sarcoma patients. Hematol Oncol. (2018) 36:591–9. doi: 10.1002/hon.2506
33. Cohen JE, Eleyan F, Zick A, Peretz T, Katz D. Intratumoral immune-biomarkers and mismatch repair status in leiyomyosarcoma -potential predictive markers for adjuvant treatment: a pilot study. Oncotarget. (2018) 9:30847–54. doi: 10.18632/oncotarget.25747
34. Zheng B, Ren T, Huang Y, Guo W. Apatinib inhibits migration and invasion as well as PD-L1 expression in osteosarcoma by targeting STAT3. Biochem Biophys Res Commun. (2018) 495:1695–701. doi: 10.1016/j.bbrc.2017.12.032
35. van Erp AEM, Versleijen-Jonkers YMH, Hillebrandt-Roeffen MHS, van Houdt L, Gorris MAJ, van Dam LS, et al. Expression and clinical association of programmed cell death-1, programmed death-ligand-1 and CD8(+) lymphocytes in primary sarcomas is subtype dependent. Oncotarget. (2017) 8:71371–84. doi: 10.18632/oncotarget.19071
36. Sundara YT, Kostine M, Cleven AH, Bovee JV, Schilham MW, Cleton-Jansen AM. Increased PD-L1 and T-cell infiltration in the presence of HLA class I expression in metastatic high-grade osteosarcoma: a rationale for T-cell-based immunotherapy. Cancer Immunol Immunother. (2017) 66:119–28. doi: 10.1007/s00262-016-1925-3
37. Que Y, Xiao W, Guan YX, Liang Y, Yan SM, Chen HY, et al. PD-L1 Expression is associated with FOXP3+ regulatory T-cell infiltration of soft tissue sarcoma and poor patient prognosis. J Cancer. (2017) 8:2018–25. doi: 10.7150/jca.18683
38. Pollack SM, He Q, Yearley JH, Emerson R, Vignali M, Zhang Y, et al. T-cell infiltration and clonality correlate with programmed cell death protein 1 and programmed death-ligand 1 expression in patients with soft tissue sarcomas. Cancer. (2017) 123:3291–304. doi: 10.1002/cncr.30726
39. Palmerini E, Agostinelli C, Picci P, Pileri S, Marafioti T, Lollini PL, et al. Tumoral immune-infiltrate (IF), PD-L1 expression and role of CD8/TIA-1 lymphocytes in localized osteosarcoma patients treated within protocol ISG-OS1. Oncotarget. (2017) 8:111836–46. doi: 10.18632/oncotarget.22912
40. Liao Y, Chen L, Feng Y, Shen J, Gao Y, Cote G, et al. Targeting programmed cell death ligand 1 by CRISPR/Cas9 in osteosarcoma cells. Oncotarget. (2017) 8:30276–87. doi: 10.18632/oncotarget.16326
41. Costa Arantes DA, Goncalves AS, Jham BC, Duarte ECB, de Paula EC, de Paula HM, et al. Evaluation of HLA-G, HLA-E, and PD-L1 proteins in oral osteosarcomas. Oral Surg Oral Med Oral Pathol Oral Radiol. (2017) 123:e188–96. doi: 10.1016/j.oooo.2016.12.002
42. Budczies J, Mechtersheimer G, Denkert C, Klauschen F, Mughal SS. PD-L1 (CD274) copy number gain, expression, and immune cell infiltration as candidate predictors for response to immune checkpoint inhibitors in soft-tissue sarcoma. Oncoimmunology. (2017) 6:e1279777. doi: 10.1080/2162402X.2017.1279777
43. Botti G, Scognamiglio G, Marra L, Pizzolorusso A, Di Bonito M, De Cecio R, et al. Programmed death ligand 1 (PD-L1) expression in primary angiosarcoma. J Cancer. (2017) 8:3166–72. doi: 10.7150/jca.19060
44. Bertucci F, Finetti P, Perrot D, Leroux A, Collin F, Le Cesne A, et al. PDL1 expression is a poor-prognosis factor in soft-tissue sarcomas. Oncoimmunology. (2017) 6:e1278100. doi: 10.1080/2162402X.2016.1278100
45. Paydas S, Bagir EK, Deveci MA, Gonlusen G. Clinical and prognostic significance of PD-1 and PD-L1 expression in sarcomas. Med Oncol. (2016) 33:93. doi: 10.1007/s12032-016-0807-z
46. Kostine M, Cleven AH, de Miranda NF, Italiano A, Cleton-Jansen AM, Bovee JV. Analysis of PD-L1, T-cell infiltrate and HLA expression in chondrosarcoma indicates potential for response to immunotherapy specifically in the dedifferentiated subtype. Mod Pathol. (2016) 29:1028–37. doi: 10.1038/modpathol.2016.108
47. Koirala P, Roth ME, Gill J, Piperdi S, Chinai JM, Geller DS, et al. Immune infiltration and PD-L1 expression in the tumor microenvironment are prognostic in osteosarcoma. Sci Rep. (2016) 6:30093. doi: 10.1038/srep30093
48. Kim C, Kim EK, Jung H, Chon HJ, Han JW, Shin KH, et al. Prognostic implications of PD-L1 expression in patients with soft tissue sarcoma. BMC Cancer. (2016) 16:434. doi: 10.1186/s12885-016-2451-6
49. Honda Y, Otsuka A, Ono S, Yamamoto Y, Seidel JA, Morita S, et al. Infiltration of PD-1-positive cells in combination with tumor site PD-L1 expression is a positive prognostic factor in cutaneous angiosarcoma. Oncoimmunology. (2017) 6:e1253657. doi: 10.1080/2162402X.2016.1253657
50. Lussier DM, O'Neill L, Nieves LM, McAfee MS, Holechek SA, Collins AW, et al. Enhanced T-cell immunity to osteosarcoma through antibody blockade of PD-1/PD-L1 interactions. J Immunother. (2015) 38:96–106. doi: 10.1097/CJI.0000000000000065
51. D'Angelo SP, Shoushtari AN, Agaram NP, Kuk D, Qin LX, Carvajal RD, et al. Prevalence of tumor-infiltrating lymphocytes and PD-L1 expression in the soft tissue sarcoma microenvironment. Hum Pathol. (2015) 46:357–65. doi: 10.1016/j.humpath.2014.11.001
52. Chowdhury F, Dunn S, Mitchell S, Mellows T, Ashton-Key M, Gray J. PD-L1 and CD8(+)PD1(+) lymphocytes exist as targets in the pediatric tumor microenvironment for immunomodulatory therapy. OncoImmunology. (2015) 4:e1029701. doi: 10.1080/2162402X.2015.1029701
53. Shen JK, Cote GM, Choy E, Yang P, Harmon D, Schwab J, et al. Programmed cell death ligand 1 expression in osteosarcoma. Cancer Immunol Res. (2014) 2:690–8. doi: 10.1158/2326-6066.CIR-13-0224
54. Kim JR, Moon YJ, Kwon KS, Bae JS, Wagle S, Kim KM, et al. Tumor infiltrating PD1-positive lymphocytes and the expression of PD-L1 predict poor prognosis of soft tissue sarcomas. PLoS ONE. (2013) 8:e82870. doi: 10.1371/journal.pone.0082870
55. Zhu L, Sun J, Wang L, Li Z, Wang L, Li Z. Prognostic and clinicopathological significance of PD-L1 in patients with bladder cancer: a meta-analysis. Front Pharmacol. (2019) 10:962. doi: 10.3389/fphar.2019.00962
56. Li H, Xu Y, Wan B, Song Y, Zhan P, Hu Y, et al. The clinicopathological and prognostic significance of PD-L1 expression assessed by immunohistochemistry in lung cancer: a meta-analysis of 50 studies with 11,383 patients. Transl Lung Cancer Res. (2019) 8:429–49. doi: 10.21037/tlcr.2019.08.04
57. Huang X, Zhang W, Zhang Z, Shi D, Wu F, Zhong B, et al. Prognostic value of programmed cell death 1 ligand-1 (PD-L1) or PD-1 expression in patients with osteosarcoma: a meta-analysis. J Cancer. (2018) 9:2525–31. doi: 10.7150/jca.25011
58. Li J, Chen L, Xiong Y, Zheng X, Xie Q, Zhou Q, et al. Knockdown of PD-L1 in human gastric cancer cells inhibits tumor progression and improves the cytotoxic sensitivity to CIK therapy. Cell Physiol Biochem. (2017) 41:907–20. doi: 10.1159/000460504
59. Chen L, Xiong Y, Li J, Zheng X, Zhou Q, Turner A, et al. PD-L1 expression promotes epithelial to mesenchymal transition in human esophageal cancer. Cell Physiol Biochem. (2017) 42:2267–80. doi: 10.1159/000480000
60. Bighetti-Trevisan RL, Sousa LO. Cancer stem cells: powerful targets to improve current anticancer therapeutics. Stem Cells Int. (2019) 2019:9618065. doi: 10.1155/2019/9618065
61. Almozyan S, Colak D, Mansour F, Alaiya A, Al-Harazi O, Qattan A, et al. PD-L1 promotes OCT4 and Nanog expression in breast cancer stem cells by sustaining PI3K/AKT pathway activation. Int J Cancer. (2017) 141:1402–12. doi: 10.1002/ijc.30834
62. Gupta HB, Clark CA, Yuan B, Sareddy G, Pandeswara S, Padron AS, et al. Tumor cell-intrinsic PD-L1 promotes tumor-initiating cell generation and functions in melanoma and ovarian cancer. Signal Transduct Target Ther. (2016) 1:16030. doi: 10.1038/sigtrans.2016.30
63. Dong P, Xiong Y, Yue J, Hanley SJB, Watari H. Tumor-intrinsic PD-L1 signaling in cancer initiation, development and treatment: beyond immune evasion. Front Oncol. (2018) 8:386. doi: 10.3389/fonc.2018.00386
64. Liu S, Chen S, Yuan W, Wang H, Chen K, Li D, et al. PD-1/PD-L1 interaction up-regulates MDR1/P-gp expression in breast cancer cells via PI3K/AKT and MAPK/ERK pathways. Oncotarget. (2017) 8:99901–12. doi: 10.18632/oncotarget.21914
65. Abdel-Rahman O. PD-L1 expression and outcome of advanced melanoma patients treated with anti-PD-1/PD-L1 agents: a meta-analysis. Immunotherapy. (2016) 8:1081–9. doi: 10.2217/imt-2016-0025
66. Li S, Chen L, Jiang J. Role of programmed cell death ligand-1 expression on prognostic and overall survival of breast cancer: a systematic review and meta-analysis. Medicine. (2019) 98:e15201. doi: 10.1097/MD.0000000000015201
67. Huang W, Ran R, Shao B, Li H. Prognostic and clinicopathological value of PD-L1 expression in primary breast cancer: a meta-analysis. Breast Cancer Res Treat. (2019) 178:17–33. doi: 10.1007/s10549-019-05371-0
68. Kelly CM, Antonescu CR, Bowler T, Munhoz R, Chi P, Dickson MA, et al. Objective response rate among patients with locally advanced or metastatic sarcoma treated with talimogene laherparepvec in combination with pembrolizumab: a phase 2 clinical trial. JAMA Oncol. (2020) 6:402–8. doi: 10.1001/jamaoncol.2019.6152
69. Iwai Y, Ishida M, Tanaka Y, Okazaki T, Honjo T, Minato N. Involvement of PD-L1 on tumor cells in the escape from host immune system and tumor immunotherapy by PD-L1 blockade. Proc Natl Acad Sci USA. (2002) 99:12293–7. doi: 10.1073/pnas.192461099
70. Meng X, Huang Z, Teng F, Xing L, Yu J. Predictive biomarkers in PD-1/PD-L1 checkpoint blockade immunotherapy. Cancer Treat Rev. (2015) 41:868–76. doi: 10.1016/j.ctrv.2015.11.001
71. Schutz F, Stefanovic S, Mayer L, von Au A, Domschke C, Sohn C. PD-1/PD-L1 pathway in breast cancer. Oncol Res Treat. (2017) 40:294–7. doi: 10.1159/000464353
72. Schalper KA. PD-L1 expression and tumor-infiltrating lymphocytes: revisiting the antitumor immune response potential in breast cancer. Oncoimmunology. (2014) 3:e29288. doi: 10.4161/onci.29288
73. Lastwika KJ, Wilson W 3rd, Li QK, Norris J, Xu H, Ghazarian SR, et al. Control of PD-L1 expression by oncogenic activation of the AKT-mTOR pathway in non-small cell lung cancer. Cancer Res. (2016) 76:227–38. doi: 10.1158/0008-5472.CAN-14-3362
74. Hall F, Villalobos V, Wilky B. Future directions in soft tissue sarcoma treatment. Curr Probl Cancer. (2019) 43:300–7. doi: 10.1016/j.currproblcancer.2019.06.004
75. Thanindratarn P, Dean DC, Nelson SD, Hornicek FJ, Duan Z. Advances in immune checkpoint inhibitors for bone sarcoma therapy. J Bone Oncol. (2019) 15:100221. doi: 10.1016/j.jbo.2019.100221
76. Torabi A, Amaya CN, Wians FH Jr, Bryan BA. PD-1 and PD-L1 expression in bone and soft tissue sarcomas. Pathology. (2017) 49:506–13. doi: 10.1016/j.pathol.2017.05.003
77. Kim ST, Klempner SJ, Park SH, Park JO, Park YS, Lim HY, et al. Correlating programmed death ligand 1 (PD-L1) expression, mismatch repair deficiency, and outcomes across tumor types: implications for immunotherapy. Oncotarget. (2017) 8:77415–23. doi: 10.18632/oncotarget.20492
78. Sheffield BS, Fulton R, Kalloger SE, Milne K, Geller G, Jones M, et al. Investigation of PD-L1 biomarker testing methods for PD-1 axis inhibition in non-squamous non-small cell lung cancer. J Histochem Cytochem. (2016) 64:587–600. doi: 10.1369/0022155416665338
Keywords: programmed cell death 1 ligand-1, sarcoma, meta-analysis, prognosis, overall survival
Citation: Wang F, Yu T, Ma C, Yuan H, Zhang H and Zhang Z (2020) Prognostic Value of Programmed Cell Death 1 Ligand-1 in Patients With Bone and Soft Tissue Sarcomas: A Systemic and Comprehensive Meta-Analysis Based on 3,680 Patients. Front. Oncol. 10:749. doi: 10.3389/fonc.2020.00749
Received: 14 January 2020; Accepted: 20 April 2020;
Published: 02 June 2020.
Edited by:
Giovanna Schiavoni, Istituto Superiore di Sanità (ISS), ItalyReviewed by:
Patrizia Gasparini, Istituto Nazionale dei Tumori (IRCCS), ItalyRobin Lewis Jones, Royal Marsden Hospital, United Kingdom
Sheima Farag, Royal Marsden Hospital, United Kingdom, in collaboration with reviewer RLJ
Copyright © 2020 Wang, Yu, Ma, Yuan, Zhang and Zhang. 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: Feng Wang, d2FuZ2ZlbmcmI3gwMDA0MDtjbXUuZWR1LmNu; Zhiyu Zhang, enl6aGFuZyYjeDAwMDQwO2NtdS5lZHUuY24=