The bioinformatics analysis of TCGA STAD dataset indicated that a high expression of YY1is related to poor GC prognosis
First, we analyzed the expression of YY1in GC tumors and adjacent normal tissues and divided all enrolled patients into two groups based on their YY1 expression. The basic characteristics of patients are summarized in Table 1.
Table 1
Basic characteristics of patients
Characteristic | Low YY1(n = 187) | High YY1(n = 188) | p |
Age, meidan (IQR) | 66 (58, 73) | 68 (59, 73) | 0.522 |
Gender, n (%) | | | 0.718 |
Female | 69 (18.4%) | 65 (17.3%) | |
Male | 118 (31.5%) | 123 (32.8%) | |
T stage, n (%) | | | 0.066 |
T1 | 12 (3.3%) | 7 (1.9%) | |
T2 | 46 (12.5%) | 34 (9.3%) | |
T3 | 86 (23.4%) | 82 (22.3%) | |
T4 | 40 (10.9%) | 60 (16.3%) | |
N stage, n (%) | | | 0.159 |
N0 | 64 (17.9%) | 47 (13.2%) | |
N1 | 48 (13.4%) | 49 (13.7%) | |
N2 | 34 (9.5%) | 41 (11.5%) | |
N3 | 31 (8.7%) | 43 (12%) | |
M stage, n (%) | | | 0.398 |
M0 | 168 (47.3%) | 162 (45.6%) | |
M1 | 10 (2.8%) | 15 (4.2%) | |
Pathologic stage, n (%) | | | 0.189 |
Stage I | 30 (8.5%) | 23 (6.5%) | |
Stage II | 59 (16.8%) | 52 (14.8%) | |
Stage III | 69 (19.6%) | 81 (23%) | |
Stage IV | 14 (4%) | 24 (6.8%) | |
Histological type, n (%) | | | 0.387 |
Diffuse Type | 36 (9.6%) | 27 (7.2%) | |
Mucinous Type | 12 (3.2%) | 7 (1.9%) | |
Not Otherwise Specified | 102 (27.3%) | 105 (28.1%) | |
Papillary Type | 1 (0.3%) | 4 (1.1%) | |
Signet Ring Type | 5 (1.3%) | 6 (1.6%) | |
Tubular Type | 31 (8.3%) | 38 (10.2%) | |
Histologic grade, n (%) | | | 0.897 |
G1 | 6 (1.6%) | 4 (1.1%) | |
G2 | 69 (18.9%) | 68 (18.6%) | |
G3 | 109 (29.8%) | 110 (30.1%) | |
Anatomic neoplasm subdivision, n (%) | | | 0.543 |
Antrum/Distal | 70 (19.4%) | 68 (18.8%) | |
Cardia/Proximal | 28 (7.8%) | 20 (5.5%) | |
Fundus/Body | 58 (16.1%) | 72 (19.9%) | |
Gastroesophageal Junction | 22 (6.1%) | 19 (5.3%) | |
Furthermore, TCGA data showed that patients with higher levels of YY1 presented significant lower postoperative PFS [HR = 1.39 (0.98 − 1.98), p = 0.065] and DSS [HR = 1.99 (1.31 − 3.02), ** p < 0.01] rates (Fig. 1A-B). We also showed that the DSS was associated with YY1 expression and that it could be used as an independent DSS indicator [HR = 1.84 (1.16, 2.91), **p < 0.01]. Hence, these results indicated the prognostic value of YY1 as a GC biomarker (Fig. 1C-D). Additionally, compared to adjacent normal tissues, the tumor tissues exhibited significantly higher YY1 expression. (t = 4.042, ***p < 0.001, Fig. 2A). Next, we performed bioinformatics analysis on YY1 using TCGA STAD dataset to clarify the underlying promotive mechanisms of YY1 on the invasion, metastasis, and drug resistance of GC cells. We detected 3286 significant DEGs between patients with high and low YY1 expression. Additionally, 105 transcriptional targets and 11 transcription factors related to YY1 were predicted using the TRRUST v2 database. Then, 4741 YY1 co-expression genes were retrieved from the cbioportal (Fig. 2B-2D). Finally, 30 genes were identified in the intersection analysis between DEGs, co-expressed genes, regulators, and targets of YY1 (Venn diagram - Fig. 2E).
Transferrin is a potential key protein regulated by YY1
To identify key genes potentially regulated by, or significantly correlated to YY1 and that presented the most important interactions with YY1, we performed a LASSO Cox analysis based on the 30 hub genes and YY1. The LASSO results led to YY1 and 5 YY1-related hub genes (Fig. 3A-B). The receiver operating characteristic (ROC) curve was used to predict the DSS of 1, 3, and 5 years and showed that this gene group might be used to predict the DSS in the cohort from TCGA (AUC of 1, 3, 5 years were 0.746, 0.747, and 0.701, respectively; Fig. 3C). Moreover, the expressions of VEGFB, DNAJB4, CXCR4, and Transferrin were positively correlated with YY1, while COX7C was negatively correlated with YY1 (Fig. 3D). The Kaplan-Meier analysis of TCGA STAD indicated that patients with high Transferrin expression were the most correlated to adverse DSS [***p < 0.01, HR = 2.17 (1.40–3.37)] and PFS [**p < 0.01, HR = 1.84 (1.28–2.65)] prognoses after surgery (Fig. 3E). A previous study has also suggested that Transferrin participates in systemic iron homeostasis. Our current results indicated that Transferrin is a potential key protein regulated by YY1.
Overexpression of YY1 regulates the activity of the p53 signaling pathway and GC cell ferroptosis.
Next, we used GeneMANIA to predict the biological network integration for gene prioritization and performed a KEGG enrichment analysis using the 30 hub genes. The KEGG analysis identified that the YY1-regulated ferroptosis was highly enriched in GC tumors, and transferrin was involved in the ferroptosis process (Fig. 2F-G). Furthermore, the GSEA showed that the TP53_ACTIVITY_THROUGH_PHOSPHORYLATION, SIGNALING_BY_TGFB_FAMILY_MEMBERS, and GLYCEROPHOSPHOLIPID_BIOSYNTHESIS were significantly enriched pathways regulated by YY1(Fig. 2H). Altogether, these results indicated that the regulatory effect of YY1 on ferroptosis might be exerted through the p53 signaling pathway.
Overexpression of YY1 might suppress the infiltration of immune cells in GC
We also explored the significance of the relationship between YY1 expression and the tumor immune microenvironment via the Pearson's correlation coefficients between YY1 and immune infiltration. The expression of YY1 was negatively correlated with IMMUNE (r = -0.220, p < 0.001), STROMAL (r = -0.095, p = 0.065), and ESTIMATE (r = -0.175, p < 0.001) scores (Fig. 4A). The CIBERSORT analysis for TCGA STAD data showed that the expression of YY1 significantly suppressed immune cells infiltration. Finally, TIMER 2.0 verified that the infiltration of CD8 + T cells (Rho=-0.136, p < 0.001), B memory cells (Rho = -0.201, p < 0.001), active NK cells (Rho = -0.165, p < 0.001), and monocytes (Rho = -0.113, p < 0.001) was significantly reduced in tumor tissues, while the infiltration of NK resting cells (Rho = 0.15, p < 0.001) enhanced (Fig. 4B-C).
Overexpression of YY1 inhibits GC cell ferroptosis and mediates Apatinib-resistance via the p53 signaling pathway
The bioinformatics results showed that Transferrin was a potential target regulated by YY1. A previous study has also suggested that Transferrin plays a role in systemic iron homeostasis. Our current GSEA showed that TP53_ACTIVITY_THROUGH_PHOSPHORYLATION was a significantly enriched pathway regulated by YY1. Further, we verified the functionality of the YY1/Transferrin axis on GC cells. The Western blot showed that YY1 overexpression directly upregulated Transferrin, inhibited p53 expression thus upregulated SLC2A11 (Fig. 5A), which might constitute the mechanism of ferroptosis inhibition after YY1 overexpression.
To demonstrate the basis for YY1 in mediating ferroptosis and Apatinib resistance, YY1-overexpressed HGC-27 and MFC cells were treated with Fer-1, Erastin, Apatinib in vitro, cell death rate and relative MDA levels were also measured. First, YY1-overexpressed and YY1-overexpressed with SLC7A11 knock down HGC-27 and MFC cells were treated with Erastin in vitro. The cell death rate and relative MDA levels showed that YY1 overexpression could inhibit GC cell ferroptosis. Meanwhile, Erastin could reverse the inhibition effect (Fig. 5B-C). What’s more, loss of SLC7A11 could block the inhibition effect (Fig. 5B-5C, *p < 0.05, ANOVA). Considering the effect of Transferrin in maintaining systemic iron homeostasis ferroptosis[27–29], we tested relative intro-cellular levels of Fe2+ and GSH. Relative GSH level in YY1-overexpressed HGC-27 and MFC cells was elevated (Fig. 5D, *p < 0.05, ANOVA). However, Fe2+ maintains low level in YY1-overexpressed cell line and relative low level in YY1-overexpressed with SLC7A11 knock down cell line (Fig. 5E, *p < 0.05, ANOVA). Which shows Transferrin and p53 might both contribute to the inhibition effect of YY1. Hence, the overexpression of YY1 could induce Apatinib drug resistance (IC50 of HGC-27: 80.94 vs 20.74 ug/ml, IC50 of MFC: 27.09 vs 10.52 ug/ml, Nonlin-Fit). These results indicated that the inhibition of ferroptosis after YY1 overexpression via p53 signaling pathway and related to elevating Transferrin.
Compared to normal control GC cells, YY1-overexpressed GC cells presented significantly enhanced growth, migration, and invasion. Thus, these results indicated the promotive effect of YY1 on GC cell growth, invasion, and metastasis (Fig. 6A-6B).
Our previous study demonstrated that Interferon-α could remodel the hepatocellular microenvironment and potentiates anti-PD-1 efficacy[25]. These findings suggested that IFN-α might be an effective treatment for Apatinib and mPD-1-resistant GC cells.
First, to verify the promotive effect of YY1 on tumor growth in vivo, we injected YY1-overexpressed GC cells and measured the tumor progression. We found that overexpression of YY1 promoted tumor progression in vivo (Fig. 6D). Additionally, we analyzed the phenotypes of GC cells with YY1 overexpression after IFN-α treatment. The IFN-α treatment did not directly lead to decreased tumor growth in vivo and did not significantly reverse both Apatinib and mPD-1 antibody resistance in YY1-overexpressed subcutaneous tumors (Fig. 6D). Finally, the IHC analysis demonstrated that, after treatment with mIFN-α, both CD8 and CD27 were significantly upregulated in GC tissues, thereby indicating an improvement in the immune microenvironment of GC tissues. However, CD19 and PD-L1 has no significant different between the two groups (Fig. 6E).