LLGL2 mRNA expression and prognosis of breast cancer patients
We first investigated the association between LLGL2 mRNA expression level and prognosis of breast cancer patients with long-term follow up. A total of 624 breast cancer tissue samples were subjected to LLGL2 mRNA expression analysis. The associations between LLGL2 mRNA expression and clinicopathological characteristics are shown in Supplementary Table S1. Low LLGL2 mRNA levels were positively associated with larger tumor size (P = 0.047) and lymph node-negativity (P = 0.031). Low LLGL2 mRNA expression was positively associated with longer DFS in all breast cancer patients analyzed in this study (P = 0.023; Supplementary Figure S3). Furthermore, patients with tumors showing low LLGL2 mRNA expression showed a tendency towards longer OS (P = 0.072; Supplementary Figure S3).
LLGL2 was reported to be involved in prognosis only in ERα-positive breast cancer patients [8]. Therefore, we next investigated the association of LLGL2 mRNA expression with prognosis according to ERα status. As shown in Fig. 1A and B, positive associations were found between low LLGL2 mRNA expression and longer DFS and OS in ERα-positive breast cancer patients (n = 491; P = 0.0009 and P = 0.005, respectively). However, no association was observed between LLGL2 mRNA expression and prognosis in ERα-negative breast cancer patients (n = 132; Supplementary Figure S3). The clinicopathological characteristics of ERα-positive breast cancer patients are shown in Table 1. Low LLGL2 mRNA expression level was positively associated with lymph node negativity (P = 0.007).
LLGL2 was reported to be involved in resistance to tamoxifen in ERα-positive breast cancer patients [8]. Therefore, we investigated the association of LLGL2 mRNA expression with prognosis in ERα-positive breast cancer patients receiving adjuvant tamoxifen therapy (n = 272). As shown in Fig. 1C and D, positive associations were found between low LLGL2 mRNA expression and longer DFS and OS in ERα-positive breast cancer patients receiving adjuvant tamoxifen therapy (P = 0.016 and P = 0.018, respectively). Interestingly, no associations were identified between LLGL2 mRNA expression level and prognosis in ERα-positive breast cancer patients without adjuvant tamoxifen therapy (Supplementary Figure S3).
We next performed univariate and multivariate Cox regression analyses of clinicopathological factors associated with prognosis using stepwise linear regression in all breast cancer patients (Supplementary Table S2) and in ERα-positive breast cancer patients analyzed in this study (Table 1). Although low LLGL2 mRNA expression was not an independent favorable prognostic factor in all breast cancer patients, we showed that low LLGL2 was an independent favorable prognostic factor for both DFS and OS in ERα-positive breast cancer patients, as well as nodal status (P = 0.012 and P = 0.011, respectively; Table 2).
SLC7A5 mRNA expression and prognosis of breast cancer patients
Next, we investigated the association between SLC7A5 mRNA expression and prognosis of breast cancer patients. The characteristics of ERα-positive breast cancer patients according to SLC7A5 mRNA expression are shown in Supplementary Table S3. Low SLC7A5 mRNA expression was positively associated with favorable prognosis in both DFS and OS in all breast cancer patients analyzed (P = 0.002 and P = 0.0005, respectively). Low SLC7A5 mRNA expression was also positively associated with favorable prognosis in both DFS and OS in ERα-positive breast cancer patients (P = 0.004 and P = 0.004, respectively; Fig. 2A and B). However, no association was observed in ERα-negative breast cancer patients (Supplementary Figure S4). As shown in Fig. 2C, positive associations were identified between low SLC7A5 mRNA expression and longer DFS in ERα-positive breast cancer patients receiving adjuvant tamoxifen therapy (P = 0.014).
Combination of LLGL2 and SLC7A5 mRNA expression and prognosis of breast cancer patients
We then investigated the prognostic impact of the combination of LLGL2 and SLC7A5 mRNA expression. Supplementary Table S4 shows the characteristics of breast cancer patients classified by the combination of LLGL2 and SLC7A5 mRNA expression. Low LLGL2/SLC7A5 mRNA co-expression (LLGL2low/SLC7A5low) was positively associated with lower tumor grade, lymph node negativity, and ERα positivity. As shown in Fig. 3A and B, LLGL2low/SLC7A5low was associated with longer survival compared with other combination groups in all breast cancer patients analyzed in this study. The characteristics of ERα-positive breast cancer patients according to LLGL2/SLC7A5 mRNA co-expression are shown in Table 3. LLGL2low/SLC7A5low was associated with lower grade and lymph node negativity in ERα-positive breast cancer patients. As shown in Fig. 3C and D, LLGL2low/SLC7A5low showed longer survival compared with high LLGL2/SLC7A5 mRNA co-expression (LLGL2high/SLC7A5high) and a positive trend of longer survival compared with other combination groups in ERα-positive breast cancer patients. Then, we investigated the association of prognosis with the combination of LLGL2 and SLC7A5 mRNA expression in ERα-positive breast cancer patients receiving adjuvant tamoxifen therapy. As shown in Fig. 3E and F, LLGL2low/SLC7A5low showed longer survival than LLGL2high/SLC7A5high and a positive trend of longer survival compared with other combination groups in ERα-positive breast cancer patients receiving adjuvant tamoxifen therapy. However, no significant difference was observed between these four combination groups in ERα-positive breast cancer patients who had not received adjuvant tamoxifen therapy (Supplementary Figure S5).
We performed univariate and multivariate Cox regression analyses of clinicopathological factors associated with prognosis using stepwise linear regression in each group of breast cancer patients. Multivariate analyses demonstrated that LLGL2low/SLC7A5low was an independent favorable prognostic factor for DFS as well as lymph node negativity and ERα positivity in all breast cancer patients analyzed (Supplementary Table S5). Then, we performed univariate and multivariate analyses in ERα-positive breast cancer patients, which identified LLGL2low/SLC7A5low as an independent favorable prognostic factor for both DFS and OS, as well as lymph node negativity (Table 4).
Protein expression of LLGL2 and SLC7A5 in breast cancer patients
The expression levels of LLGL2 protein in breast cancer tissue samples were examined by IHC. LLGL2 protein expression was observed in the cytoplasm. Representative images of LLGL2 and SLC7A5 are as shown in Fig. 4A. Immunostaining results were evaluated using the Aperio scanscopeCS2 and eSlide manager application, and H-scores were calculated by this digital pathological system. In this study, a total of 285 consecutive breast cancer tissue samples for which mRNA expression data were available were analyzed for LLGL2 protein expression.
We also evaluated the expression levels of SLC7A5 protein in breast cancer tissue samples. The analysis was performed using the same tissue microarray as the protein expression analysis of LLGL2. SLC7A5 protein expression was observed in the cell membrane (Fig. 4A). Furthermore, SLC7A5 protein expression was observed in 10% of breast cancer tissues. We investigated the association between prognosis and the combination of LLGL2 and SLC7A5 protein expression. The LLGL2high/SLC7A5pos group seemed to show the worst prognosis among the four groups, although there was no statistically significant difference between them (Fig. 4B). The median H-score was used as the cutoff value for LLGL2, and SLC7A5 was divided into two groups based on the presence or absence of staining.