Correlation Between PD-L2 Expression and Clinical Outcome in Solid Cancer Patients: A Meta-Analysis

Background: Immune checkpoint inhibitors targeting the programmed cell death 1 (PD-1)/programmed cell death ligand 1 (PD-L1) pathway are a paradigm-shifting cancer therapy. Programmed cell death ligand 2 (PD-L2) is another ligand of PD-1, but its prognostic significance in solid cancer patients after surgery remains controversial. In this study, we aimed to reveal the prognostic implication of PD-L2 in solid tumors through a meta-analysis.

Methods: We searched PubMed, Embase and the Cochrane library for studies reporting the relationship between PD-L2 expression and prognosis or clinicopathological features in solid cancer patients after surgery from inception to January 2018, with language restricted to English. Pooled hazard ratios (HRs) and 95% confidence intervals (CIs) were determined to explore the prognostic value of PD-L2 expression. Odds ratios (ORs) were also calculated to investigate the relationship between PD-L2 expression and clinicopathological parameters.

Results: Sixteen studies incorporating 3,533 patients were included in our meta-analysis. The pooled results revealed that PD-L2 overexpression was a weak negative predictor for overall survival (OS; HR = 1.38, 95% CI = 1.05–1.81, P = 0.021), as well as a strong predictor for poor disease-free survival (DFS)/progression-free survival (PFS) (HR = 1.44, 95% CI = 1.15–1.81, P = 0.001). In subgroup analyses, high PD-L2 expression revealed an unfavorable prognostic prediction for OS in hepatocellular carcinoma (HCC) (HR = 1.60, 95% CI = 1.12–2.29, P = 0.011) and for DFS/PFS in HCC (HR = 1.50, 95%CI = 1.04–2.16, P = 0.031) as well as clear cell renal cell carcinoma (HR = 1.45, 95% CI = 1.03–2.03, P = 0.033). Moreover, PD-L2 expression implied a weak trend toward the presence of lymphatic metastasis (presence vs. absence, OR = 1.61, 95% CI = 0.98–2.65, P = 0.061).

Conclusion: High PD-L2 expression may promote tumor metastasis and predict unfavorable prognosis in solid cancer patients after surgery, especially in HCC.

Introduction

Over the past decade, one of the most remarkable breakthroughs in cancer therapy has been tumor immunotherapy, in particular immune checkpoint inhibitory programmed cell death 1 (PD-1)/programmed cell death ligand 1 (PD-L1) blockade (1). Anti-PD-1/PD-L1 immunotherapy strengthens antitumor immunity by “releasing the brakes” for immune suppression in the tumor microenvironment and has exhibited inspiring efficacy in various cancer types, resulting in FDA approval and wide clinical application (2–4). PD-1 has been reported to possess two ligands, namely PD-L1 and PD-L2 (5), and the role of the PD-1/PD-L1 axis in cancer has been investigated in depth. It has been reported that PD-L1 predicts unfavorable prognosis in patients with various tumors (6). Nevertheless, PD-L2 has gained little attention, and its function in tumor immunity remains unclear.

Akin to PD-L1, PD-L2 interacts with PD-1 and suppresses T cell proliferation and cytokine release (5). Moreover, PD-L2 can bind to another exclusive partner, namely repulsive guidance molecule b, to promote immune tolerance, although the role of this pathway in cancer is obscure. Contrary to the extensive expression profile of PD-L1, PD-L2 is more confined to antigen-presenting cells, although it can be induced by microenvironmental stimuli (7). Additionally, PD-L2 expression has been discovered in a large number of tumor types (8, 9). One study that analyzed archival tumor samples across seven cancer types showed that in triple-negative breast cancer samples, PD-L2 expression was strongly consistent with PD-L1 expression, with no significantly discordant expression between PD-L1 and PD-L2 in any sample. Among other tumor types, although researchers observed discordant expression in a few samples, PD-L2 expression was also significantly associated with PD-L1 expression (10). While in esophageal adenocarcinoma, PD-L2 expression was identified in approximately half of tumor samples, together with only 2% PD-L1 positivity observed (11). Therefore, although PD-L1 may plays a more dominant role in immunologic modulation of the tumor microenvironment, PD-L2 may have been neglected as a potential target in tumor immunity.

Recently, studies have increasingly focused on the prognostic implications of PD-L2 in cancer patients after surgery. However, whether PD-L2 expression correlates with prognosis in solid cancer patients after surgery remains elusive. Some evidence supports that high PD-L2 expression in cancer specimens predicts impaired survival for various solid tumors (12–16). Nevertheless, several studies report negative results or even opposing findings (11, 17–24). Therefore, we assessed the consistency and magnitude of the prognostic effect of PD-L2 in solid cancer patients after surgery through a meta-analysis. We found that PD-L2 was a negative predictor for prognosis among solid cancer patients. To our knowledge, this study is the first to clarify the prognostic significance of PD-L2 expression in solid tumors by meta-analysis.

Materials and Methods

Literature Search Strategy

The implementation of this meta-analysis was conducted in accordance with PRISMA guidelines (25). We systematically reviewed PubMed, Embase, and the Cochrane library for literature published up to January 2018, with the language restricted to English. Keywords adopted to execute the retrieval are as follows: (programmed cell death 1 ligand 2 OR PD-L2 OR B7-DC) AND (cancer OR neoplasm OR malignancy OR carcinoma OR tumor) AND (prognostic OR prognosis). We also checked the references and citations of the retrieved papers. To ensure the reliability of the search results, three authors (HYY & XXZ & LJS) independently executed the retrieval according to the standardized process. Any disagreement between the three authors was resolved through discussion.

Inclusion and Exclusion Criteria

We adopted the following inclusion criterion: (1) studies published as original articles; (2) studies published in English; (3) studies reporting a correlation between PD-L2 expression in tumor specimens and solid tumor prognosis, such as overall survival (OS), disease-free survival (DFS), and progression-free survival (PFS) or clinicopathological characteristics; (4) available and clear prognostic data from which to directly extract hazard ratios (HR) and 95% confidence intervals (CI); (5) studies had a sample capacity of more than 50 individuals; (6) all enrolled patients received surgery with or without any adjuvant therapy. Studies were excluded if they failed to conform to any inclusion criteria. The analysis of PD-L2 expression in tumor infiltrating lymphocytes was also excluded. For duplicate publications, we included the most informative publication.

Data Collection and Quality Assessment

Three authors (HYY & XXZ & LJS) independently examined all selected publications and extracted the data in accordance with a standardized protocol. From every study, we pulled the following data: clinicopathological characteristics; first author's name; publication year; country of the patients; number of individuals; cancer type; trial design; cutoff value of high PD-L2; median follow-up time and the range of follow-up time; outcomes of patients; method adopted to probe PD-L2 expression; univariate or multivariate model. If there existed a multivariate HR and a univariate HR, the former one was chosen to avoid the influence of confounding factors.

Also, three authors (HYY & XXZ & LJS) independently assessed the quality of the selected publications on the basis of the Newcastle-Ottawa scale (NOS) criteria (25). Each study was graded on a scale of zero to nine according to the selection, comparability, and outcomes of the study cohorts.

Statistical Analyses

As DFS and PFS are very similar, we considered the two parameters together. We applied Stata 12.0 software (Stata Corporation, College Station, Texas, USA) to perform the statistical analysis. We applied the log HRs and 95% CIs for aggregation of the prognostic effects (26). We utilized the χ² test and I² statistic to assess the heterogeneity across studies (27). If P <0.10 for the χ² test or I² > 50%, significant heterogeneity was considered to exist and the random effects model was utilized (28); If not, a fixed-effects model was utilized (29). We also performed a sensitivity analysis in which one study was deleted every time to judge its impact on the results. We used Begg's funnel plot, Begg's test and Egger's tests to investigate the publication bias quantitatively (30, 31). We utilized the non-parametric “trim and fill” approach to evaluate the potential impact of publication bias, which considered hypothetical negative unpublished studies and recalculated a pooled estimate that comprised these hypothetical studies (32). For all analyses, two-sided P-values of <0.05 was regarded as statistically significant.

Results

Study Selection

We selected 142 studies by systematic literature retrieval. After screening the title and abstract, 112 studies were excluded for their irrelevance to the topic. Finally, we included 16 studies after reviewing the full texts (11–24, 33, 34). Studies were excluded for the following reasons: six had no information about OS/DFS/PFS or reported no clinicopathological data; six had insufficient data for quantitative analysis; and two were reviews. The details of our search process are presented in Figure 1.

FIGURE 1

Figure 1. Flow chart of the study screening process.

Characteristics of the Included Studies

We included data from 3,533 patients in this study. The main characteristics of the 16 studies are summed up in Table 1. One study was published in 2009, and the others were published between 2014 and 2017. Four studies were prospective cohort trials, and 12 were retrospective cohort trials. Patients in all studies were reported to receive resections. These studies were conducted in six countries. With regard to cancer type, studies investigating renal cell carcinoma (RCC; n = 3) and gastric cancer (n = 3) composed the two largest proportion among all included investigations, followed by hepatocellular carcinoma (HCC; n = 2), non-small cell lung cancer (NSCLC; n = 2), esophageal cancer (n = 2), breast cancer (n = 1), oral squamous cell cancer (OSCC; n = 1), neurological cancer (n = 1), and colorectal cancer (CRC; n = 1). More specifically, Two of RCCs were clear cell RCC (ccRCC) and one of RCCs was chromophobe RCC (chRCC); two of gastric cancers were gastric adenocarcinoma (GA) and one of gastric cancers was mixed with tubular adenocarcinoma (TA) and signet ring cell (SRC); one of NSCLCs was squamous cell carcinoma (SqCC) and another one was adenocarcinoma; one of esophageal cancers was SqCC and another was adenocarcinoma. Retrospective trials were designed in 12 studies, and prospective trials were designed in four studies. The cutoff values of high PD-L2 were discordant, while the most common criterion was the median score. With respect to spatial location, all the included studies focus PD-L2 expression on tumor cell. The median follow-up time ranged from 16 months to 7.18 years. The correlation between OS and PD-L2 expression was reported in 13 of the studies. The prognostic value of PD-L2 expression for DFS/PFS was reported in seven studies. In addition, the relationship between clinicopathological features and PD-L2 expression was presented in 13 studies. Seven studies calculated the HRs adjusted for PD-1 expression or PD-L1 expression, and nine studies didn't adjust for PD-1 expression or PD-L2 expression. All included studies used immunohistochemistry (IHC) to examine PD-L2 expression. The origins of PD-L2 antibodies utilized for IHC in the included studies varied, while 7 of the studies used the same antibody (clone 176611, R&D Systems). Nonetheless, only five studies have definitely checked the specificity of the PD-L2 antibody utilized for IHC on a positive control and none of the included studies has set negative control for PD-L2 antibody (Table S1). All included studies were allocated scores >5 on the Newcastle-Ottawa scale (NOS), suggesting that all possessed high methodological quality (Table 2).

TABLE 1

Table 1. Main characteristics of the eligible studies.

TABLE 2

Table 2. The Newcastle-Ottawa scale (NOS) quality assessment of the enrolled studies.

Prognostic Value of PD-L2 Expression for OS

Thirteen studies consisting of 2,845 patients reported OS. This study revealed that PD-L2 overexpression was a weak negative predictor for OS among patients with various solid tumors (HR = 1.38, 95% CI = 1.05–1.81, P = 0.021) (Figure 2A). This pooled meta-analysis was carried out using the random effects model on account of significant heterogeneity (I² = 64.2%, P = 0.001). To further explore the potential sources of heterogeneity, we utilized subgroup analyses, which are summarized in Figure 2B.

FIGURE 2

Figure 2. (A) Meta-analysis of the association between PD-L2 expression and OS among solid cancer patients after surgery; (B) Subgroup analyses of the correlation between PD-L2 and OS.

Subgroup analyses regarding cancer type clarified that high PD-L2 expression had a unfavorable prognostic value for OS in patients with HCC (HR = 1.60, 95% CI = 1.12–2.29, P = 0.011), while no significant association was observed in RCC, esophageal cancer or gastric cancer (OS for RCC: HR = 1.86, 95% CI = 0.55–6.27, P = 0.315; esophageal cancer: HR = 0.91, 95% CI = 0.60–1.38, P = 0.647; gastric cancer: HR = 1.33, 95% CI = 0.81–2.17). Notably, since there was only one study each of breast cancer, OSCC, neurological cancer, and CRC, we did not perform subgroup analyses in every cancer type. To further restricted the prognostic effect in homogeneous cancer histology, we performed subgroup analyses according to cancer histology. The pooled results showed that PD-L2 overexpression is not significantly associated with OS in gastric adenocarcinoma (HR = 1.63, 95% CI = 0.81–3.26, P = 0.168). By the same token, we didn't perform subgroup analyses in every cancer histology. Moreover, subgroup analyses regarding model and trial design was performed. The pooled results revealed that PD-L2 was a negative predictor for OS among studies adjusted for PD-1 expression or PD-L1 expression (HR = 1.42, 95% CI = 1.05–1.90, P = 0.021), whereas it was not significantly correlated with OS among studies not adjusted for PD-1/PD-L1 (HR = 1.39, 95% CI = 0.84–2.30, P = 0.201). Similarly, the pooled results were also not consistent between the subgroups divided by trial design (retrospective trial design: HR = 1.45, 95% CI = 1.02–2.07, P = 0.041; prospective trial design: HR = 1.22, 95% CI = 0.78–1.90, P = 0.387).

Prognostic Value of PD-L2 Expression for DFS/PFS

The HRs in DFS/PFS were available in seven studies comprising 1,030 patients. Similar to OS, PD-L2 overexpression in tumor samples showed unfavorable DFS/PFS (HR = 1.44, 95% CI = 1.15–1.81, P = 0.001) (Figure 3A). The fixed-effects model was applied as no significant heterogeneity was observed (I² = 21.4%, P = 0.267).

FIGURE 3

Figure 3. (A) Meta-analysis of the association between PD-L2 expression and DFS/PFS among solid cancer patients after surgery; (B) Subgroup analyses of the correlation between PD-L2 and DFS/PFS.

Because the grouping results of cancer type and cancer histology were identical, we directly performed subgroup analyses regarding cancer histology. The stratified analyses indicated that high PD-L2 expression predicted poor DFS/PFS among patients with HCC (HR = 1.50, 95% CI = 1.04–2.16, P = 0.031) and ccRCC (HR = 1.45, 95% CI = 1.03–2.03, P = 0.033). We did not carry out subgroup analyses for breast cancer, OSCC, or CRC because of the lack of more than one study reporting the associations between prognosis and each cancer histology. The results of subgroup analysis regarding the model for DFS/PFS was consistent with those of OS, namely, PD-L2 was an unfavorable predictor for DFS/PFS among studies adjusted for PD-1/PD-L1 (HR = 1.62, 95% CI = 1.24–2.13, P < 0.001), whereas it was not significantly correlated with DFS/PFS among studies not adjusted for PD-1/PD-L1 (HR = 1.12, 95% CI = 0.75-1.67, P = 0.591) (Figure 3B).

Correlation Between PD-L2 Expression and Clinicopathological Characteristics

To comprehensively dissect the role of PD-L2 expression as a biomarker in solid tumors, we explored the correlation between PD-L2 expression and clinicopathological characteristics. Thirteen studies comprising 3,282 patients reported the relationship between PD-L2 expression and clinicopathological parameters. The pooled results indicated that PD-L2 expression implied a weak trend toward the presence of lymphatic metastasis (presence vs. absence, OR = 1.61, 95% CI = 0.98–2.65, P = 0.061). Meanwhile, PD-L2 overexpression had no significant association with sex (male vs. female, OR = 1.08, 95% CI = 0.88–1.34, P = 0.467), depth of invasion (TIII+TIV vs. TI+TII, OR = 0.99, 95% CI = 0.47–2.01, P = 0.467), histopathological stage (III+IV vs. I+II, OR = 0.99, 95% CI = 0.51–1.92, P = 0.968), tumor metastasis (presence vs. absence, OR = 1.07, 95% CI = 0.70–1.63, P = 0.757), vascular invasion (presence vs. absence, OR = 1.27, 95% CI = 0.67–2.40, P = 0.459), recurrence (presence vs. absence, OR = 1.89, 95% CI = 0.72–5.01, P = 0.198), differentiation (poor vs. moderate or well, OR = 0.74, 95% CI = 0.48–1.15, P = 0.178) or tumor size (≥5 vs. <5 cm, OR = 1.04, 95% CI = 0.72–1.50, P = 0.851) (Table 3).

TABLE 3

Table 3. Meta-analysis of reported clinicopathological characteristics in the included studies.

Sensitivity Analysis and Publication Bias

All studies were successively omitted to judge the robustness of the pooled results. Furthermore, we combined Begg's funnel plot, Begg's test and Egger's test to evaluate whether a publication bias existed. As shown in Figure 4, the significance of the recalculated HRs did not change when any study was deleted except for Gao's study, which was published in 2017. When Gao's study was omitted, the HRs of the pooled OS became borderline to the statistical threshold (HR = 1.29, 95% CI = 0.99–1.69, P = 0.063), which indicated that the correlation between PD-L2 expression and OS was not robustly significant. Of note, Egger's test of OS indicated the existence of publication bias (P = 0.019), and visual estimation of Begg's funnel plot of OS revealed evident asymmetry, together with the results that Begg's test of OS was closed to the threshold of statistical significance (P = 0.059; Figure 5A). Hence, we performed a further sensitivity analysis with the trim and fill method to evaluate the effect of the potential publication bias of OS. As shown in Figure 5B, two hypothetical negative studies were added to diminish the asymmetry of Begg's funnel plot. However, the pooled results comprising the hypothetical resulting in losing the statistical significance of the correlation between PD-L2 expression and OS (HR = 1.274, 95% CI = 0.96–1.69). Therefore, we must be careful in drawing a conclusion regarding OS. As to DFS/PFS, the significance of the recalculated HRs did not change when any study was deleted (Figure 6A). In addition, the shape of the Begg's funnel plot did not display evident asymmetry (Figure 6B), indicating that no obvious publication bias existed, which was confirmed by Begg's test (P = 0.764) and Egger's test (P = 0.409).

FIGURE 4

Figure 4. Sensitivity analysis of the effect of individual studies on the pooled HRs for PD-L2 and OS.

FIGURE 5

Figure 5. Begg's funnel plots for publication bias of PD-L2 on OS. (A) Publication bias for OS without trim and fill; (B) publication bias for OS with trim and fill.

FIGURE 6

Figure 6. (A) Sensitivity analysis of the effect of individual studies on the pooled HRs for PD-L2 and DFS/PFS; (B) Funnel plots for publication bias of PD-L2 on DFS/PFS.

Discussion

Recently, a growing number of studies have investigated the prognostic implication of PD-L2 protein expression in tumors among patients with solid tumors, and the results remain controversial. A transcriptomic analysis performed by Danilova and colleagues used The Cancer Genome Atlas datasets and found that PD-L2 mRNA expression indicated good prognosis in skin cutaneous melanoma (SKCM) but had no prognostic role in other 8 tumors (35). However, mRNA expression may not accurately represent the protein expression. In this meta-analysis, we included studies reporting PD-L2 protein expression and revealed that PD-L2 overexpression was a weak negative predictor for OS and a strong negative predictor for DFS/PFS in solid cancer patients after surgery. Similarly, PD-L2 displayed inconsistent prognostic significance in different cancer types and histology. PD-L2 positive expression showed unfavorable prognostic prediction for OS in HCC and for DFS/PFS in HCC and ccRCC; however, the prognostic value remained not significant for other cancer types and histology. In addition, high PD-L2 expression implied a weak trend toward the presence of lymphatic metastasis. Taken together, high PD-L2 expression might promote tumor metastasis and predict unfavorable prognosis in solid cancer patients after surgery, especially in HCC.

Contrary to PD-L1, PD-L2 is less well-studied, and the underlying mechanism involved in the correlation between PD-L2 positive expression and decreased survival remains unclear. Previous studies showed that PD-L2 exerts its main physiological and pathological function in immune tolerance via dampening and modulating T helper type 2 (Th2) response (36–38). However, Th1 immuno-responses are the most dominant in the context of antitumor immunity. In esophageal squamous cell carcinoma (ESCC), PD-L2 was shown to negatively associate with the number of PD-1+ tumor infiltrating lymphocytes (TILs), indicating that PD-L2 expressed on ESCC cells might restrain PD-1+ TILs activity, ultimately facilitating immune escape (39). This inverse correlation between PD-L2 expression and abundance of TILs has also been reported elsewhere (40). Recently, a multivariate analysis revealed that in CRC, high tumor PD-L2 expression was correlated with a weak Crohn's-like lymphoid reaction, which was deemed to play a significant role in adaptive immune responses against tumors, suggesting that PD-L2 expression might reduce Crohn's-like lymphoid reactions to suppress antitumor immunity (41). Moreover, Pinato et al. revealed that PD-L2 upregulation was associated with tumor hypoxia (14). As tumor hypoxia is generally recognized as a contributor to immune resistance in the tumor microenvironment (42), it is conceivable that PD-L2 may potentially mediate immune resistance via promoting tumor hypoxia. Furthermore, high PD-L2 expression in tumor tissue was independently associated with better clinical response to pembrolizumab in patients with head and neck squamous cell carcinoma, suggesting that PD-L2 expression in tumor tissue may exert a vital role in response to PD-1 blockade therapy, although this study was not included in our meta-analysis (10). However, whether PD-L2 has a causal relationship with immune suppression in tumor immunity as well as the underlying mechanism requires further investigation.

In recent years, immune checkpoint inhibitors, particularly those targeting the PD-1 pathway, have become a paradigm-shifting therapy in cancer treatment. Blockade of PD-L2 in a pancreatic murine model also displayed evident anti-tumor effects with decreased tumor outgrowth rates (43). In contrast, a PD-L2 knockout mouse bearing a CT26 tumor exhibited a weakened tumor-specific cytotoxic T lymphocyte response and more rapid tumor growth compared with a WT mouse bearing CT26 (44). These inconsistent outcomes leave the targeting potential of PD-L2 unresolved. Recently, Ahmad and colleagues described CD4+ and CD8+ PD-L2-specific T cells, which can directly exert cytotoxic activity against PD-L2-expressing target cells and indirectly release pro-inflammatory cytokines into the tumor microenvironment, indicating the possibility of novel PD-L2-based vaccines complementary to checkpoint inhibitors. The results of our meta-analysis further suggest that PD-L2 is a pejorative prognostic biomarker for solid tumors, especially in HCC. Further investigations studying the mechanism of the role of PD-L2 in solid tumor immunity are required to understand the targeting potential of PD-L2.

Regarding the tumor types and tumor histology, the unfavorable prognostic effect of PD-L2 was consistently significant both for OS and DFS/PFS in HCC, which has startling heterogeneity and lacks an efficient therapeutic approach (45, 46), suggesting that PD-L2 may serve as prognostic marker especially in HCC. However, in esophageal cancer, high PD-L2 expression implied a favorable prognosis trend, although with no statistical significance. This discordant result suggested that perhaps PD-L2 has different effects on immune suppression among different cancer types. More studies reporting the prognostic roles of PD-L2 in specific cancer types and histology are thus required. Given that PD-L2 expression is widely reported to be correlated with PD-L1 expression and PD-L1 has been reported to predict pejorative prognosis in solid tumors (6, 10), it is crucial to exclude the confounding effect of PD-L1 before we draw conclusion for PD-L2. Indeed, subgroup analyses regarding model revealed that the prognostic role of PD-L2 for both OS and DFS/PFS remained significant in studies adjusted for PD-1 expression or PD-L1 expression, indicating that PD-L2 overexpression predicts poor prognosis in solid tumors independent of PD-1/PD-L1 expression and other confounding factors. Of note, all studies included in our study utilized IHC to determine PD-L2 expression in specimens. IHC is the major approach by which to probe the expression level of protein in resected tumor samples and has been widely applied. Therefore, our findings are very feasible for clinical practice.

Our study had several limitations. First, sensitivity analyses revealed that the correlation between PD-L2 expression and OS was not robustly stable, which might be explained by the publication bias that existed in this meta-analysis. To validate our hypothesis, we utilized the trim and fill method to include two hypothetical negative results, revealing that the correlation between PD-L2 expression and OS lost significance. It seemed that because positive findings tend to be published compared with negative findings, the association between PD-L2 expression and OS might potentially be exaggerated. Nevertheless, the association between PD-L2 expression and DFS/PFS was robust and stable, with no publication bias observed. Considering that PD-L2 expression has an unstable weak significant prognostic effect for OS and stable significant prognostic effect for DPS/PFS, we believe that our study reveals meaningful statistical evidence endorsing the prognostic value of PD-L2 to predict unfavorable outcome in solid cancer patients after surgery. Second, the OS results displayed heterogeneity in our meta-analysis, which may be because of the different analysis strategy of IHC and heterogenous cut-off values of PD-L2 expression. However, we performed subgroup analyses in an attempt to minimize the impact of heterogeneity. We found that PD-L2 was an unfavorable predictor of OS in a homogeneous HCC population. In addition, because of the extensively heterogenous cut-off values of PD-L2 expression, we were unable to stratify the studies according to cut-off values. Furthermore, most of the included studies didn't check the specificity of PD-L2 antibody utilized for IHC, which might weaken the reliability of the result. Thus, to remove the heterogeneity in current clinical protocols, a unified approach for protein expression assessment and a rigorous check for antibody specificity are urgently required. Last, all studies selected for this meta-analysis focused on PD-L2 expressed by tumor cells. However, whether PD-L2 expressed by tumor cells or other cells in the tumor microenvironment plays the dominant role is unclear. PD-L2 has also been found to be expressed by stromal cells and appears to be functional (47, 48). It should be noted that perhaps not only PD-L2 expression in tumor cells but also its expression in stromal cells plays a significant role in immune suppression and affects prognosis. We highly recommend further studies to investigate the prognostic roles of PD-L2 expressed by stroma cells in addition to tumor cells.

In conclusion, our study has clinical significance because it clarifies the correlation between high PD-L2 expression and unfavorable prognosis for solid cancer patients after surgery, especially in HCC.

Author Contributions

HY designed the study, performed the literature search and screening, performed the data analyses and wrote the manuscript. XZ designed the study, retrieved the literature and data, analyzed the retrieved data and participated in the writing of manuscript. LS assisted in the designing of the study, performed the literature search and screening, assisted in the data analyses and participated in the writing of manuscript. YM designed the study and supervised the study. Every author approved the final version of this study.

Funding

This study was supported by grants from the CAMS Innovation Fund for Medical Sciences (CIFMS) (No.2016-I2M-1-001) and the National High-tech Research and Development Projects (863) (No.2015AA020303).

Conflict of Interest Statement

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.

Acknowledgments

We thank H. Nikki March, Ph.D, from Liwen Bianji, Edanz Editing China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.00047/full#supplementary-material

References

1. Postow MA, Callahan MK, Wolchok JD. Immune checkpoint blockade in cancer therapy. J Clin Oncol. (2015) 33:1974–82. doi: 10.1200/JCO.2014.59.4358

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Herzberg B, Campo MJ, Gainor JF. Immune checkpoint inhibitors in non-small cell lung cancer. Oncologist. (2017) 22:81–88. doi: 10.1634/theoncologist.2016-0189

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Rittmeyer A, Barlesi F, Waterkamp D, Park K, Ciardiello F, von Pawel J, et al. Atezolizumab versus docetaxel in patients with previously treated non-small-cell lung cancer (OAK): a phase 3, open-label, multicentre randomised controlled trial. Lancet (2017) 389:255–65. doi: 10.1016/S0140-6736(16)32517-X

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Balar AV, Galsky MD, Rosenberg JE, Powles T, Petrylak DP, Bellmunt J, et al. Atezolizumab as first-line treatment in cisplatin-ineligible patients with locally advanced and metastatic urothelial carcinoma: a single-arm, multicentre, phase 2 trial. Lancet (2017) 389:67–76. doi: 10.1016/S0140-6736(16)32455-2

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Latchman Y, Wood CR, Chernova T, Chaudhary D, Borde M, Chernova I, et al. PD-L2 is a second ligand for PD-1 and inhibits T cell activation. Nat Immunol. (2001) 2:261–8. doi: 10.1038/85330

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Wang Q, Liu F, Liu L. Prognostic significance of PD-L1 in solid tumor: An updated meta-analysis. Medicine (2017) 96:e6369. doi: 10.1097/MD.0000000000006369

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Lesterhuis WJ, Punt CJ, Hato SV, Eleveld-Trancikova D, Jansen BJ, Nierkens S, et al. Platinum-based drugs disrupt STAT6-mediated suppression of immune responses against cancer in humans and mice. J Clin Invest. (2011) 121:3100–8. doi: 10.1172/JCI43656

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Ahmad SM, Martinenaite E, Holmstrom M, Jorgensen MA, Met O, Nastasi C, et al. The inhibitory checkpoint, PD-L2, is a target for effector T cells: Novel possibilities for immune therapy. Oncoimmunology (2018) 7:e1390641. doi: 10.1080/2162402X.2017.1390641

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Panjwani PK, Charu V, DeLisser M, Molina-Kirsch H, Natkunam Y, Zhao S. Programmed death-1 ligands PD-L1 and PD-L2 show distinctive and restricted patterns of expression in lymphoma subtypes. Hum Pathol. (2018) 71:91–9. doi: 10.1016/j.humpath.2017.10.029

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Yearley JH, Gibson C, Yu N, Moon C, Murphy E, Juco J, et al. PD-L2 Expression in Human Tumors: Relevance to Anti-PD-1 Therapy in Cancer. Clin Cancer Res. (2017) 23:3158–67. doi: 10.1158/1078-0432.CCR-16-1761

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Derks S, Nason KS, Liao X, Stachler MD, Liu KX, Liu JB, et al. Epithelial PD-L2 expression marks Barrett's esophagus and esophageal adenocarcinoma. Cancer Immunol Res. (2015) 3:1123–9. doi: 10.1158/2326-6066.CIR-15-0046

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Jung HI, Jeong D, Ji S, Ahn TS, Bae SH, Chin S, et al. Overexpression of PD-L1 and PD-L2 is associated with poor prognosis in patients with hepatocellular carcinoma. Cancer Res Treat. (2017) 49:246–54. doi: 10.4143/crt.2016.066

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Gao Y, Li S, Xu D, Chen S, Cai Y, Jiang W, et al. Prognostic value of programmed death-1, programmed death-ligand 1, programmed death-ligand 2 expression, and CD8(+) T cell density in primary tumors and metastatic lymph nodes from patients with stage T1-4N+M0 gastric adenocarcinoma. Chin J Cancer (2017) 36:61. doi: 10.1186/s40880-017-0226-3

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Pinato DJ, Black JR, Trousil S, Dina RE, Trivedi P, Mauri FA, et al. Programmed cell death ligands expression in phaeochromocytomas and paragangliomas: Relationship with the hypoxic response, immune evasion and malignant behavior. Oncoimmunology (2017) 6:e1358332. doi: 10.1080/2162402X.2017.1358332

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wang H, Yao H, Li C, Liang L, Zhang Y, Shi H, et al. PD-L2 expression in colorectal cancer: Independent prognostic effect and targetability by deglycosylation. Oncoimmunology (2017) 6:e1327494. doi: 10.1080/2162402X.2017.1327494

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Shin SJ, Jeon YK, Kim PJ, Cho YM, Koh J, Chung DH, et al. Clinicopathologic analysis of PD-L1 and PD-L2 expression in renal cell carcinoma: association with oncogenic proteins status. Ann Surg Oncol. (2016) 23:694–702. doi: 10.1245/s10434-015-4903-7

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wu Y, Cao D, Qu L, Cao X, Jia Z, Zhao T, et al. PD-1 and PD-L1 co-expression predicts favorable prognosis in gastric cancer. Oncotarget (2017) 8:64066–82. doi: 10.18632/oncotarget.19318

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Tanaka K, Miyata H, Sugimura K, Kanemura T, Hamada-Uematsu M, Mizote Y, et al. Negative influence of programmed death-1-ligands on the survival of esophageal cancer patients treated with chemotherapy. Cancer Sci. (2016) 107:726–33. doi: 10.1111/cas.12938

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Shin SJ, Jeon YK, Cho YM, Lee JL, Chung DH, Park JY, et al. The association between PD-L1 expression and the clinical outcomes to vascular endothelial growth factor-targeted therapy in patients with metastatic clear cell renal cell carcinoma. Oncologist (2015) 20:1253–60. doi: 10.1634/theoncologist.2015-0151

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Baptista MZ, Sarian LO, Derchain SF, Pinto GA, Vassallo J. Prognostic significance of PD-L1 and PD-L2 in breast cancer. Hum Pathol. (2016) 47:78–84. doi: 10.1016/j.humpath.2015.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Dong M, Wang HY, Zhao XX, Chen JN, Zhang YW, Huang Y, et al. Expression and prognostic roles of PIK3CA, JAK2, PD-L1, and PD-L2 in Epstein-Barr virus-associated gastric carcinoma. Hum Pathol. (2016) 53:25–34. doi: 10.1016/j.humpath.2016.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Erlmeier F, Weichert W, Autenrieth M, Wiedemann M, Schrader AJ, Hartmann A, et al. PD-L2: A prognostic marker in chromophobe renal cell carcinoma? Med Oncol. (2017) 34:71. doi: 10.1007/s12032-017-0926-1

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Gao Q, Wang XY, Qiu SJ, Yamato I, Sho M, Nakajima Y, et al. Overexpression of PD-L1 significantly associates with tumor aggressiveness and postoperative recurrence in human hepatocellular carcinoma. Clin Cancer Res. (2009) 15:971–9. doi: 10.1158/1078-0432.CCR-08-1608

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Kogashiwa Y, Yasuda M, Sakurai H, Nakahira M, Sano Y, Gonda K, et al. PD-L1 expression confers better prognosis in locally advanced oral squamous cell carcinoma. Anticancer Res. (2017) 37:1417–24. doi: 10.21873/anticanres.11465

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Wells GA SB, O'Connell D, Peterson J, Welch V, Losos M, et al. The Newcastle-Ottawa Scale (NOS) for Assessing the Quality if Nonrandomized Studies in Meta-Analyses (2009). Available online at: http://www.ohri.ca/programs/clinical_epidemiology/oxford.htm

26. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med. (2002) 21:1539–58. doi: 10.1002/sim.1186

PubMed Abstract | CrossRef Full Text | Google Scholar

28. DerSimonian R, Kacker R. Random-effects model for meta-analysis of clinical trials: an update. Contemp Clin Trials (2007) 28:105–14. doi: 10.1016/j.cct.2006.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst. (1959) 22:719–48.

PubMed Abstract | Google Scholar

30. Begg CB, Mazumdar M. Operating characteristics of a rank correlation test for publication bias. Biometrics (1994) 50:1088–101. doi: 10.2307/2533446

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Egger M, Davey Smith G, 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

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Duval S, Tweedie R. Trim and fill: A simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics (2000) 56:455–63. doi: 10.1111/j.0006-341X.2000.00455.x

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Kim S, Kim MY, Koh J, Go H, Lee DS, Jeon YK, et al. Programmed death-1 ligand 1 and 2 are highly expressed in pleomorphic carcinomas of the lung: Comparison of sarcomatous and carcinomatous areas. Eur J Cancer (2015) 51:2698–707. doi: 10.1016/j.ejca.2015.08.013

CrossRef Full Text | Google Scholar

34. Zhang Y, Wang L, Li Y, Pan Y, Wang R, Hu H, et al. Protein expression of programmed death 1 ligand 1 and ligand 2 independently predict poor prognosis in surgically resected lung adenocarcinoma. Onco Targets Ther. (2014) 7:567–73. doi: 10.2147/OTT.S59959

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Danilova L, Wang H, Sunshine J, Kaunitz GJ, Cottrell TR, Xu H, et al. Association of PD-1/PD-L axis expression with cytolytic activity, mutational load, and prognosis in melanoma and other solid tumors. Proc Natl Acad Sci USA (2016) 113:E7769–77. doi: 10.1073/pnas.1607836113

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Ansari MJ, Salama AD, Chitnis T, Smith RN, Yagita H, Akiba H, et al. The programmed death-1 (PD-1) pathway regulates autoimmune diabetes in nonobese diabetic (NOD) mice. J Exp Med. (2003) 198:63–9. doi: 10.1084/jem.20022125

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Keir ME, Liang SC, Guleria I, Latchman YE, Qipo A, Albacker LA, et al. Tissue expression of PD-L1 mediates peripheral T cell tolerance. J Exp Med. (2006) 203:883–95. doi: 10.1084/jem.20051776

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Ishiwata K, Watanabe N, Guo M, Tomihara K, Brumlik MJ, Yagita H, et al. Costimulator B7-DC attenuates strong Th2 responses induced by Nippostrongylus brasiliensis. J Immunol. (2010) 184:2086–94. doi: 10.4049/jimmunol.0804051

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Leng C, Li Y, Qin J, Ma J, Liu X, Cui Y, et al. Relationship between expression of PD-L1 and PD-L2 on esophageal squamous cell carcinoma and the antitumor effects of CD8+ T cells. Oncol Rep. (2016) 35:699–708. doi: 10.3892/or.2015.4435

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Sridharan V, Gjini E, Liao X, Chau NG, Haddad RI, Severgnini M, et al. Immune profiling of adenoid cystic carcinoma: PD-L2 expression and associations with tumor-infiltrating lymphocytes. Cancer Immunol Res. (2016) 4:679–87. doi: 10.1158/2326-6066.CIR-16-0031

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Masugi Y, Nishihara R, Hamada T, Song M, da Silva A, Kosumi K, et al. Tumor PDCD1LG2 (PD-L2) expression and the lymphocytic reaction to colorectal cancer. Cancer Immunol Res. (2017) 5:1046–55. doi: 10.1158/2326-6066.CIR-17-0122

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Noman MZ, Hasmim M, Messai Y, Terry S, Kieda C, Janji B, et al. Hypoxia: a key player in antitumor immune response. A review in the theme: cellular responses to hypoxia. Am J Physiol Cell Physiol. (2015) 309:C569–79. doi: 10.1152/ajpcell.00207.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Okudaira K, Hokari R, Tsuzuki Y, Okada Y, Komoto S, Watanabe C, et al. Blockade of B7-H1 or B7-DC induces an anti-tumor effect in a mouse pancreatic cancer model. Int J Oncol. (2009) 35:741–9. doi: 10.3892/ijo_00000387

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Shin T, Yoshimura K, Shin T, Crafton EB, Tsuchiya H, Housseau F, et al. In vivo costimulatory role of B7-DC in tuning T helper cell 1 and cytotoxic T lymphocyte responses. J Exp Med. (2005) 201:1531–41. doi: 10.1084/jem.20050072

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Vijgen S, Terris B, Rubbia-Brandt L. Pathology of intrahepatic cholangiocarcinoma. Hepatobiliary Surg Nutr. (2017) 6:22–34. doi: 10.21037/hbsn.2016.11.04

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Weiskirchen R. Intratumor heterogeneity, variability and plasticity: questioning the current concepts in classification and treatment of hepatocellular carcinoma. Hepatobiliary Surg Nutr (2016) 5:183–7. doi: 10.3978/j.issn.2304-3881.2016.02.04

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Nazareth MR, Broderick L, Simpson-Abelson MR, Kelleher RJ Jr, Yokota SJ, Bankert RB. Characterization of human lung tumor-associated fibroblasts and their ability to modulate the activation of tumor-associated T cells. J Immunol. (2007) 178:5552–62. doi: 10.4049/jimmunol.178.9.5552

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Yang H, Bueso-Ramos C, DiNardo C, Estecio MR, Davanlou M, Geng QR, et al. Expression of PD-L1, PD-L2, PD-1 and CTLA4 in myelodysplastic syndromes is enhanced by treatment with hypomethylating agents. Leukemia (2014) 28:1280–8. doi: 10.1038/leu.2013.355

PubMed Abstract | CrossRef Full Text | Google Scholar