Upregulation of DARS2 Expression in BLCA
We analyzed TCGA data to investigate the expression of DARS2 in tumor and normal tissues. Our findings revealed that DARS2 expression was significantly higher than that in normal tissues across 28 types of tumors, including BLCA, ACC, BRCA, CESC, CHOL, COAD, DLBC, ESCA, GBM, HNSC, KICH, KIRC, LAML, LGG, LIHC, LUAD, LUSC, OV, PAAD, PRAD, READ, SKCM, STAD, TGCT, THCA, THYM, UCEC, and UCS (Fig. 1A). Furthermore, our analysis of TCGA data demonstrated a significant upregulation of DARS2 expression in BLCA (P < 0.001, Fig. 1B). Additionally, analysis of paired mRNA expression data from TCGA showed a significant increase in DARS2 expression in tumor tissues (P < 0.001, Fig. 1C).
Subsequently, qPCR was performed to validate the mRNA expression levels of DARS2 in BLCA using paired samples from five patients. The results showed an increase in DARS2 mRNA expression in four pairs of cancer tissues compared to adjacent normal tissues (Fig. 1D). Protein immunoblotting results further demonstrated elevated expression of DARS2 in both cancer cell lines and tumor tissues (Fig. 1E-F). IHC analysis of ten pairs of bladder cancer tissues revealed increased expression of DARS2 in bladder cancer compared to normal bladder mucosal epithelial tissues (P < 0.05, Fig. 1G-H).
Subsequently, patients were grouped based on clinical characteristics to determine the correlation between DARS2 expression levels and clinical features. TCGA analysis results revealed significant differences in higher DARS2 expression across histological grade (P < 0.001), pathological stage (P < 0.05), T stage (P < 0.05), and M stage (P < 0.01) (Supplement Table 1). Furthermore, to further confirm the relationship between clinical features of tumors and DARS2 protein expression, we conducted clinical baseline data analysis based on the results of IHC analysis in 37 patient samples. The patients were divided into high and low DARS2 expression groups based on the median DARS2 expression. The analysis showed significant differences in higher DARS2 expression concerning N stage (P < 0.05), M stage (P < 0.05), and histological grade (P < 0.05) (Supplement Table 2). Notably, the results for M stage and histological grade from IHC analysis were consistent with the TCGA results.
These findings suggest that DARS2 expression is upregulated in BLCA at both the transcriptional and translational levels.
DARS2 as an Independent Prognostic Factor in BLCA
To determine the prognostic value of DARS2 in BLCA, we divided patients in the TCGA dataset into DARS2 low-expression and DARS2 high-expression groups based on the median DARS2 expression for survival analysis. High DARS2 expression was associated with poorer overall survival (OS) in BLCA (HR = 1.48, P < 0.01) (Fig. 2A).
To further validate the prognostic role of DARS2, we analyzed follow-up data from collected patients. Based on IHC analysis results, patients were categorized into high-expression and low-expression groups. The analysis revealed that patients in the high-expression group had worse OS (HR = 3.02, P < 0.01) (Fig. 2B).
Next, we performed univariate and multivariate Cox regression analyses on clinical data from TCGA to determine the correlation between overall survival in BLCA and multiple factors. Univariate analysis showed that six clinical features, including T3 stage (HR = 1.970, P < 0.001), T4 stage (HR = 2.987, P < 0.001), N stage (HR = 2.250, P < 0.001), M stage (HR = 3.112, P < 0.005), age > 70 (HR = 1.424, P < 0.05), and DARS2 expression (HR = 1.480, P < 0.005), were significantly associated with patient overall survival (Table 3; Fig. 2C). Multivariate analysis data revealed that DARS2 expression (HR = 1.953, P = 0.022) is an independent prognostic factor (Supplement Table 3).
Identification and Enrichment Analysis of Differentially Expressed Genes
Through gene differential analysis, a total of 466 genes were identified as differentially expressed genes (DEGs) between the high DARS2 group and the low DARS2 group (Fig. 3A). To identify signaling pathways regulated by abnormal DARS2 expression, we compared the DARS2 high-expression and low-expression groups using a signature gene set based on the TCGA dataset.
The analysis results indicated the top five upregulated pathways as follows: Formation of the cornified envelope; WP retinoblastoma gene in cancer; G2 mDNA damage checkpoint; Processing of DNA double-strand break ends; Meiosis (Fig. 3B)
Conversely, the top five downregulated pathways were identified as Initial triggering of complement; Complement cascade; Scavenging of heme from plasma; Creation of C4 and C2 activators; and CD22 mediated by regulation (Fig. 3C)
These findings provide insights into the potential pathways and processes influenced by DARS2 dysregulation in BLCA, offering valuable information for further mechanistic investigations and therapeutic targeting.
DARS2 Expression Correlates with Immune Infiltration and PD-L1 Expression
Our analysis using TIMER2 showed that DARS2 expression was negatively correlated with immune-active cells, including CD4+ T cells (R= -0.251, P < 0.001) and NK cells (R= -0.067, P < 0.001). Conversely, DARS2 expression was positively correlated with immunosuppressive cells, such as MDSCs (R = 0.372, P < 0.001) and macrophages (R = 0.196, P < 0.01). Interestingly, DARS2 expression was positively correlated with CD8+T cells (R = 0.203, P < 0.001) (Fig. 4A).
It is known that blocking immune checkpoint receptors such as PD-1/PD-L1 and CTLA-4 can alleviate CD8+T cell exhaustion and reactivate immune cell cytotoxicity to eliminate antigen-expressing tumor cells (Bibliography Farhood, et al. 2019). Therefore, we analyzed the correlation between DARS2 expression and the expression of PD-1, PD-L1 (CD274), and CTLA-4 using TIMER2. We found that DARS2 expression was positively correlated with PD-L1 expression (R = 0.202, P < 0.001), but showed no significant correlation with PD-1 and CTLA-4 (Fig. 4B). Additionally, protein immunoblotting from five pairs of bladder cancer and adjacent normal tissues showed increased expression of both PD-L1 and DARS2 in cancer tissues (Fig. 1F).
Therefore, DARS2 expression plays a significant role in immune infiltration and may serve as a potential biomarker for immune therapy response in BLCA patients.
Knockdown of DARS2 Inhibits Cell Proliferation, Migration, and Invasion, and downgraded PD-L1 Expression
The above results indicate that DARS2 is upregulated in BLCA, suggesting its crucial role in BLCA tumorigenesis. Therefore, we used siRNA to knock down DARS2 expression (Fig. 5A) to observe its biological effects on T24 and EJ cells. We assessed the impact of DARS2 on cell proliferation using CCK8 and EDU assays. Our results showed that knocking down DARS2 inhibited cell proliferation (Fig. 5B-C). Furthermore, we found that migration and invasion abilities were reduced in T24 and EJ cells with reduced DARS2 expression (Fig. 5D).
In addition, since our immune infiltration analysis revealed a positive correlation between DARS2 expression and PD-L1, we examined PD-L1 expression in cells after DARS2 knockdown. Interestingly, we observed a decrease in PD-L1 expression in both bladder cancer cell lines after knocking down DARS2 (Fig. 5E). However, the exact mechanisms underlying this decrease are currently unclear.
Finally, we conducted animal experiments and found that T24 cells with DARS2 knockdown formed smaller subcutaneous tumors in nude mice compared to the NC group (Fig. 5F). IHC images of tumor tissues showed a decrease in both DARS2 and the immune marker PD-L1 after DARS2 knockdown (Fig. 5G).
In summary, these results indicate that besides affecting cell proliferation, migration, and invasion, DARS2 can also participate in the BLCA process by influencing the cellular immune microenvironment.