- Department of Urology, The First Affiliated Hospital of Harbin Medical University, Harbin, China
Background: Pre-mRNA processing factor 19 (PRPF19) is an E3 ligase that plays a crucial role in repairing tumor-damaged cells and promoting cell survival. However, the predictive value and biological function of PRPF19 in bladder urothelial carcinoma (BLCA) require further investigation.
Methods: In this study, we utilized transcriptomic data and bladder cancer tissue microarrays to identify the high expression of PRPF19 in BLCA, suggesting its potential as a prognostic biomarker. To gain a better understanding of the role of PRPF19 in the immune microenvironment of BLCA, we performed single cell analysis and employed the LASSO method. Additionally, we examined the methylation profiles of PRPF19 using the SMART website. Our investigation confirmed the correlation between PRPF19 and BLCA cell senescence and stemness. Furthermore, we constructed a PRPF19-miR-125a-5p-LINC02693-MIR4435-2HG ceRNA network using the ENCORI and miRWALK databases.
Results: Our comprehensive analysis reveals that PRPF19 can serve as a prognostic marker for BLCA and is significantly associated with various immune-infiltrating cells in BLCA. Moreover, our findings suggest that PRPF19 influences cellular senescence through the regulation of stemness. Finally, we developed a ceRNA network that has the potential to predict the prognosis of BLCA patients.
Conclusion: We confirmed the prognostic value and multiple biological functions of PRPF19 in BLCA. Furthermore, the specific ceRNA network can be used as a potential therapeutic target for BLCA.
1 Introduction
Bladder cancer, positioned among the top ten most prevalent cancers globally, encompasses muscle-invasive bladder cancer (MIBC) and non-muscle-invasive bladder cancer (NMIBC). Despite the use of standard treatments such as radical cystectomy or concurrent radiotherapy (1), the prognosis for patients diagnosed with bladder cancer remains unfavorable. Therefore, it is crucial to conduct a comprehensive investigation into the underlying mechanisms that lead to the development and progression of bladder cancer. Such an analysis is of utmost importance in facilitating screening and preventive measures.
PRPF19, also referred to as pre-mRNA processing factor 19, functions as an E3 ligase, which regulates cellular damage repair (2). Furthermore, it helps to avoid cell cycle arrest and prolong cell survival (3, 4). E3 ubiquitin ligases play a vital part in various cellular procedures. They are responsible for protein degradation, controlling protein-protein interactions and substrate activation or inactivation. Dysfunctional E3 frequently leads to cancer development (5, 6). This potential link to cancer immunotherapy shows promise for treating cancer by targeting E3 ubiquitin ligases. PRPF19, an E3 ubiquitin ligase, promotes colorectal cancer liver metastasis through activation of the Src-YAP1 pathway via K63-linked MYL9 ubiquitination (7). PRPF19 functions as an E3 ligase, serving as a fresh therapeutic target gene for pathological cardiac hypertrophy to inhibit cardiomyocyte hypertrophy (8). The heightened invasiveness of neuroblastoma cells may be linked to PRPF19’s regulation of the Hippo-YAP pathway (9). Additionally, PRPF19 governs MDM4 mRNA splicing in order to alter p53-dependent cellular senescence (10). Moreover, PRPF19 has been discovered to advance the development of tongue cancer and alter resistance to radiotherapy (11). Nevertheless, there exists a deficiency of research on the correlation between PRPF19 and bladder cancer, leading to the need for further study.
In this work, we examined bladder cancer data from TCGA and GEO to gather clinical data and RNA-seq gene expression information (12). We identified the expression pattern and functional network of PRPF19 in BLCA using multidimensional analysis. In addition, infiltrating immune cells linked to PRPF19 that may be relevant for prognosis were found using a machine learning technique. We investigated the relevance of PRPF19 to senescence-associated secretory phenotype (SASP)-related genes in bladder cancer. We also looked into how methylation of PRPF19 affects bladder cancer. Finally, we anticipated the PRPF19-connected ceRNA network using PRPF19 as the primary molecule. This offers a foundation for figuring out potential PRPF19-related molecular pathways.
2 Materials and methods
2.1 Samples, datasets and antibody
RNA-sequencing expression profiles (level 3) and associated clinical information for BLCA were acquired from TCGA (13), and GEO databases. The TCGA database consists of 406 bladder cancer samples and 19 normal bladder tissues, while the GSE13507 data set includes 165 bladder cancer samples and 9 normal bladder tissues. Genetic alterations data were obtained from the TCGA database. Bladder cancer tissue microarrays were purchased from Shanghai Outdo Biotech Company, and a total of 59 cases were included. Samples with lost visits or incomplete information were deleted (patients with unclear clinical staging and lost to follow-up), and finally 56 bladder cancer samples (including 32 paired normal bladder cancer tissues) were included in the study. All patients were operated from May 2007 to November 2011, and the follow-up was in March 2014, lasting 2.3-7 years. The study was ethically approved and informed consent was obtained from the patients. PRPF19 antibody (ab126776) was purchased from Abcam.
2.2 Prognostic analysis
A log-rank test was used to examine the variations in survival rates between these groups. The findings of each variable (P-value, HR, and 95% CI) were presented in the form of a forest plot using the “forestplot” program after a multivariate Cox regression analysis.
2.3 Immunohistochemical analysis of PRPF19 expression in bladder cancer tissue microarrays
After dewaxing the microarray, endogenous peroxidase blocker was added dropwise to the tissue and incubated at room temperature for 10 minutes for antigen repair, then after sealing was completed PRPF19 antibody (Dilution ratio is 1:60) was added and incubated overnight at 4°C, washed with PRPF19 antibody and then incubated with the secondary antibody for 30 minutes, and then finally the color was developed, re-stained, and sealed. Finally, the staining results were independently scored by 2 pathologists and the intensity of staining was categorized as low, medium, and high (assigned a score of 1,2,3 in that order) and the range of staining was categorized as 0-25%, 26%-50%, 51%-75%, and 76%-100% (assigned a score of 1,2,3, and 4 in that order).Finally, the two results were multiplied to obtain a total score, with ≤6 being low expression and >6 being high expression.
2.4 Gene enrichment analysis
The R software’s limma package was used to study the mRNA with differential expression. The threshold for differential expression of mRNAs was set at “adjusted P 0.05.” ClusterProfiler is a R utility that evaluates the GO function of candidate mRNAs and enriches the KEGG pathway, giving information on target genes’ oncogenic involvement. Another use for the ClusterProfiler tool is gene set enrichment analysis (GSEA) (14).
2.5 Immune infiltration analysis and machine learning
The present study aimed to investigate the potential association between PRPF19 expression and immune cells by analyzing the Tumor Immunization Single Cell Center (TISCH) database (15). This analysis aimed to shed light on the potential effects of PRPF19 at the individual cell level. Furthermore, the XCELL approach was used to select crucial immune cells from a group of 38 immune cells utilizing the Least Absolute Shrinkage and Selection Operator (LASSO) (16). The Kaplan-Meier survival analysis was used to analyze the disparity in survival rates between the aforementioned groups. Additionally, a time-dependent receiver operating characteristic (timeROC) study was performed to assess the predictive accuracy of the model.
2.6 DNA methylation analysis
The UALCAN database was used to get data on the variance in PRPF19 methylation levels between bladder cancer and normal bladder tissues (17). In addition, the SMARP database offered valuable insights into the variations in expression and prognostic implications of methylation probes targeting PRPF19 in the context of bladder cancer (18).
2.7 Stemness score of tumor cells
The determination of mRNAsi, as developed by Malta et al., may be achieved by the use of the OCLR technique. The creation of the gene expression profile is informed by the mRNA expression signature, which encompasses a total of 11,774 genes. The study used Spearman correlations to assess the similarity of RNA expression patterns. The dryness index is bounded inside the interval [0,1] due to the computation of subtracting the minimum value from the highest value and then dividing the result by the maximum value (19).
2.8 CeRNA network analysis
The MiRWALK platform (20) was used in our study to examine the microRNAs (miRNAs) that have the ability to target PRPF19.The probable interaction between screened miRNAs and long non-coding RNAs (lncRNAs) was predicted using ENCORI (21). In addition, the subcellular localization of ceRNAs was investigated with the Genecards program. The building of the ceRNA network was carried out based on the ceRNA hypothesis, which posits a negative association between mRNAs and lncRNAs and miRNAs. This was supported by RNA expression analysis and total survival analysis of the BLCA cohort.
2.9 Statistical analysis
All statistical analyses were conducted using version 4.0.3 of the R statistical software. The Wilcox test was utilized to compare the statistical disparity of two groups. The Kaplan-Meier curves for survival analysis underwent log-rank test. A p-value less than 0.05 was deemed to be statistically significant. *p < 0.05, **p < 0.01 and ***p < 0.001.
3 Result
3.1 PRPF expression and mutation analysis in BLCA
Initially, a comparative analysis of PRPF19 mRNA expression was conducted in bladder cancer and normal bladder tissues using the TCGA and GEO databases. The results showed a significant increase in PRPF19 mRNA expression in bladder cancer tissues compared to normal bladder tissues. This observation held true for both bladder cancer and GSE13507 primitive bladder tumor samples inside the TCGA database (Figures 1A, B). Notable alterations were seen in the expression of PRPF19 mRNA when comparing low-grade and high-grade bladder cancer samples (Figures 1C, D). Our analysis of the correlation between PRPF19 expression and the survival outcomes of bladder cancer patients, using the TCGA-BLCA dataset and the GSE13507 dataset, revealed that patients with high PRPF19 expression exhibited a significantly unfavorable prognosis in terms of both overall survival and disease-specific survival (Figures 1E–H). A multifactorial COX regression analysis was conducted, including age, gender, T-stage, N-stage, and M-stage as variables. The results of this analysis provided confirmation that PRPF19 serves as a valuable predictive biomarker for patients with bladder cancer (Figures 1I, J). PRPF19 mRNA levels were significantly higher in bladder cancer samples compared to paired paraneoplastic tissues in the TCGA bladder cancer dataset (Figure 1K). The present study used Oncoplot to visually represent the somatic landscape of the TCGA bladder cancer cohort, whereby samples were classified into high and low PRPF19 expression categories. The genes in question exhibited a higher frequency of mutations in the group with high PRPF19 expression, as shown by their ranking based on the number of mutations they underwent (Figure 1L). Therefore, it is plausible that PRPF19 might potentially have an impact on other genetic modifications and function as a prognostic indicator for persons afflicted with bladder cancer.
Figure 1 PRPF19 is significantly expressed in bladder cancer and is linked to high mutation frequency and a poor prognosis. (A, B) PRPF19 mRNA expression level in TCGA-BLCA databases and GSE13507 databases. (C, D) Expression levels of PRPF19 mRNA in low-grade and high-grade samples of TCGA-BLCA databases and GSE13507 databases. (E, F) Correlation between PRPF19 and overall survival in bladder cancer patients. (G, H) Correlation between PRPF19 and disease-free survival in bladder cancer patients. (I, J) Multifactorial COX regression analysis of the prognostic value of PRPF19 in bladder cancer samples. (K) Expression levels of PRPF19 in paired samples of bladder cancer in TCGA database. (L) Top 10 mutated genes in PRPF19 high expression and PRPF19 low expression groups from the TCGA database, shown as a heat map. **p < 0.01 and ***p < 0.001.
3.2 Validation of PRPF19 prognosis and expression in bladder cancer
To confirm the previously mentioned findings, we performed immunohistochemistry tests to investigate the expression and prognostic relevance of PRPF19 in bladder cancer. A set of 56 bladder cancer tissue samples and 32 normal bladder tissue samples were analyzed to further assess the expression and prognostic significance of PRPF19. The outcomes of our research are in line with bioinformatics-based investigations already conducted. Our observations reveal a noteworthy elevation in PRPF19 expression in bladder cancer tissues as against normal bladder tissues. (Figures 2A–D). In the context of paired bladder cancer tissue samples, it was shown that the expression of PRPF19 was markedly elevated compared to normal bladder tissue (Figure 2E). The difference in PRPF19 expression between bladder cancer and normal bladder cancer samples was shown by the use of box plots (Figure 2F). In this study, we conducted an analysis using Table 1 to examine the relationship between PRPF19 expression and several clinical parameters in patients diagnosed with bladder cancer. The findings of our analysis provided confirmation that PRPF19 expression exhibited a significant association with both pTNM-stage and status. Following a comprehensive examination of the survival duration among patients categorized into the PRPF19 high-expression and PRPF19 low-expression groups, alongside an assessment of the ratio of survival to mortality, it was observed that patients in the PRPF19 low-expression group exhibited significantly prolonged survival periods compared to those with high PRPF19 expression (Figure 2G). Furthermore, it is worth noting that the mortality rate observed in the PRPF19 low expression group was much lower compared to the PRPF19 high expression group, as seen in Figure 2H. Figure 2I presents a Sankey diagram illustrating the distribution of high and low expression levels of PRPF19 in bladder cancer samples with respect to several factors such as age, gender, tumor size, stage, lymph node metastasis, and patient survival. Table 2 presents the results of univariate COX regression analysis, which revealed significant associations between the prognosis of bladder cancer patients and variables such as pTNM-Stage, external aggression, and PRPF19. Furthermore, the findings from multivariate COX regression analysis indicated that PRPF19 holds potential as a prognostic biomarker for patients with bladder cancer. Patients exhibiting elevated levels of PRPF19 expression have a markedly worse prognosis, as seen in Figure 2J. In addition, it was ascertained that the expression of PRPF19 exhibited predictive capabilities for 1-, 3-, and 5-year survival in patients with bladder cancer, as shown by the utilization of nomograms and the assessment of area under the curve (AUC) (Figures 2K–M). In summary, our findings have substantiated the elevated expression of PRPF19 in bladder cancer and established its potential use as a predictive biomarker for individuals diagnosed with this malignancy.
Figure 2 Tissue microarray analysis confirms high expression of PRPF19 in bladder cancer and its prognostic value. (A–D) Expression of PRPF19 in bladder cancer and normal tissues analyzed by immunohistochemical staining. (E) Expression of PRPF19 in paired bladder cancer tissues. (F) Expression of PRPF19 in 56 bladder cancer tissues and 32 normal bladder tissues. (G) Analysis of the difference in survival time between patients in the high PRPF19 expression group and low PRPF19 expression group. (H) Analysis of the ratio of survival/death of patients in the high PRPF19 expression group and low PRPF19 expression group. (I) Trends in the distribution of bladder cancer patients by age, sex, and other clinical factors. (J) Kaplan-Meier survival analysis of PRPF19 expression. (K–M) Analysis of the predictive ability of PRPF19 expression on the 1,3,5-year prognosis of bladder cancer patients. **p < 0.01 and ***p < 0.001.
3.3 Enrichment analysis of PRPF19 in bladder cancer
In the TCGA-BLCA dataset, samples from patients with bladder cancer were classified into high-expression and low-expression groups of PRPF19, based on the level of PRPF19 expression. The purpose of this study was to investigate the potential involvement of PRPF19 in bladder cancer (Figures 3A, B). The results of the GO enrichment analysis indicated that PRPF19 governs the regulation of vital biological processes, including nuclear division and the cell cycle (Figures 3C, D). The findings from the KEGG analysis demonstrate that PRPF19 has a regulatory role in several aspects of bladder cancer cell biology, including senescence, cell cycle progression, and resistance to platinum-based drugs (Figures 3E, F).
Figure 3 Analysis of the functions of PRPF19 in bladder cancer. (A) High-expression and low-expression subsets of PRPF19 were used in enrichment analysis and volcano mapping, respectively. (B) With the PRPF19 high expression and PRPF19 low expression groups, an enrichment analysis and a heat map were done. (C, D) Results of differential gene GO term enrichment. (E, F) Differential gene KEGG pathway enrichment results.
3.4 Gene set enrichment analysis of PRPF19 in BLCA
Gene set enrichment analysis (GSEA) offered more elucidation into the role of PRPF19 in the context of bladder cancer. The results obtained from GSEA demonstrated a statistically significant correlation between PRPF19 and Sarscov2 innate immunity evasion and cellspecific immune response, additionally, PRPF19 was shown to be associated with Mhc pathway (Figures 4A, B).T lymphocytes get activated and launch a targeted immune response upon the presentation of antigenic peptides by MHC molecules, resulting in the identification of both MHC molecules and antigenic peptides by T cells (22).Furthermore, recent research has shown evidence that the major histocompatibility complex (MHC) may function as a biomarker for immune checkpoint inhibitors, as indicated by many studies (23).Therefore, we posit that PRPF19 may have a substantial role in the immunological microenvironment of bladder cancer via the MHC route. In the interim, our examination of gene set enrichment further confirmed the noteworthy correlation between PRPF19 and various biological processes, including DNA methylation, senescence-associated secretory phenotype, senescence induced by oncogenes, pluripotency pathways in embryonic stem cells, senescence-related glycolysis, cell cycle, and the E2F pathway in BLCA (Figures 4C–I).
Figure 4 GSEA results confirm the involvement of PRPF19 in multiple biological functions of bladder cancer. (A) Sarscov2 innate immunity evasion and cellspecific immune response. (B) Mhc pathway. (C) DNA methylation. (D) Senescence associated secretory phenotype. (E) Oncogene induced senescence. (F) Embryonic stem cell pluripotency pathway. (G) Glycolysis in senescence. (H) cell cycle. (I) E2f pathway.
3.5 Correlation between PRPF19 expression and immune infiltration of bladder cancer
In order to examine the correlation between the distribution of immune cells and the expression levels of PRPF19 at the individual cell level, we obtained two distinct datasets related to BLCA from the scRNA-seq TISCH database (Figure 5A). Elevated levels of PRPF19 expression were seen in actively proliferating CD4 Tconv, CD8T, and NK cells within the BLCA GSE145281 aPDL1 and BLCA GSE149652 datasets, as shown in Figures 5A, B. The location and expression of PRPF19 in different immune cells were determined using clustered plots of single-cell RNA sequencing data (Figure 5C). The results of this study indicate a positive correlation between the expression levels of PRPF19 and the different kinds of immune cells seen in BLCA.
Figure 5 Investigating PRPF19 expression in BLCA’s various immune cells using scRNA-seq. (A) Independent scRNA-seq database heatmap showing PRPF19’s association with immune cell infiltration levels. (B) Violin plot of PRPF19 and immune cell infiltration. (C) Single-cell atlas including all of the cells that are included in the GSE145281 and GSE149652 datasets, as well as the expression and dispersion of PRPF19.
3.6 Analysis of key immune cells associated with PRPF19 in bladder cancer
The presence of immune cells inside tumors, known as immune infiltration, has a substantial influence on tumor development. The composition of immune cells residing in tumors has been acknowledged as a pivotal factor in determining the efficacy of cancer therapy (24). The XCELL deconvolution methodology was used to determine the relative abundance of 38 distinct immune cell types present within each BLCA sample. The examination of PRPF19 expression using subgroup analysis has shown a statistically significant correlation between PRPF19 and 16 discrete subcategories of infiltrating immune cells (Figure 6A). The LASSO technique was then used to develop prognostic models using a total of 8 of these 16 immune cells (Figures 6B, C). The risk score = (35.811)*Mast cell+(0.191)*Monocyte+(-7.033)*T cell NK+(-3.223)* Class-switched memory B cell+(1.052)* T cell CD4+Th1+(2.962)* stroma score+(-10.547)* T cell CD4+ central memory+(-7.270)* T cell CD8+. The BLCA patients were categorized into two groups based on their risk ratings. Figure 6D depicts the expression patterns of the eight genes over a range of risk scores, survival rates, and patient outcomes. Based on the Kaplan-Meier survival curves, it can be shown that patients belonging to the high-risk group had a less favorable outcome compared to those in the low-risk group (Figure 6E). The AUC values for the durations of 1, 3, and 5 years were seen to be 0.679, 0.694, and 0.640, respectively (Figure 6F).
Figure 6 Identification of the key infiltrated immune cells by LASSO. (A) Immune infiltration varies across groups with high and low PRPF19 expression. (B) Diagram depicting penalties for 38 immune cells using the LASSO model. (C) Selective feature coefficients are shown as a function of the lambda parameter. (D) Eight prognostic immune infiltration-related genes and their expression patterns in BLCA patients with varying risk scores and survival outcomes. (E, F) ROC curves and overall survival curves for high- and low-risk BLCA patients, respectively. *p < 0.05, **p < 0.01 and ***p < 0.001.
3.7 Relationship between prognostic PRPF19 expression and critical immune infiltrate cells
The results of our univariate prognostic analysis, utilizing the major immune infiltrating cells collected from the previous section (Figure 7A), demonstrated a significant association between mast cell, stroma score, T cell CD4+ central memory, T cell CD8+, and prognosis in patients with BLCA. Subsequently, an examination was conducted to assess the connection among the 16 immune infiltrating cells. The results indicated that a significant proportion of these cells had a positive relationship with one another. This finding was substantiated by our examination of the association among the 16 immune invading cells, as seen in Figure 7B. Violin plots were used to demonstrate the notable dissimilarity in expression levels of the four prognosis-related immune infiltration cells between the PRPF19 high expression group and the PRPF19 low expression group (Figure 7C). In conclusion, we conducted a comparative analysis of the expression levels of four immune infiltrating cells and PRPF19. Our findings indicate a positive association between Mast cell and PRPF19 expression, whereas a negative correlation was seen between the expression levels of the other three immune infiltrating cells (Figure 7D).
Figure 7 Validation and analysis of the persistence of the main infiltrating immune cells. (A) Analysis using univariate COX regression of the prognostic value of sixteen immune cells. (B) 16 immune cell correlation analysis. (C) Prognosis-related immune cell differences between the PRPF19 high and PRPF19 low expression groups. (D) Correlation of prognosis-related immune cell expression with PRPF19 expression. *p < 0.05 and ***p < 0.001.
3.8 An investigation of the relationship between the expression of PRPF19 and the methylation of DNA
The UALCAN database was used to assess the presence of a statistically significant difference in PRPF19 methylation levels between bladder cancers and normal bladder tissues. It was observed that a difference existed between the methylation levels of PRPF19 in bladder tumors and normal bladder tissues. Specifically, the methylation levels of PRPF19 were shown to be much lower in bladder tumors compared to normal bladder tissues. This disparity became evident when analyzing bladder tumor patients at stage 2 of the disease (Figures 8A, B). Additionally, it was shown that the methylation levels of PRPF19 in bladder tumors of stages N0-3 were significantly lower compared to those in normal bladder tissue (Figure 8C). The SMART database was employed for analysis in order to investigate the involvement of PRPF19 methylation in bladder cancer. This analysis involved examining the chromosomal distribution of methylation probes associated with PRPF19, as well as exploring comprehensive genomic information related to PRPF19 (Figures 8D, E). The statistical analysis revealed a significant difference in the expression levels of 10 out of the 14 methylation probes between bladder cancer and normal bladder tissue, as seen in Figure 8F. The predictive relevance of the 10 differentially methylated probes in bladder cancer was then assessed, revealing a high association between cg07509094 and cg26375147 and the prognosis of individuals with bladder cancer (Figure 8G). Subsequently, the expression patterns of the aforementioned methylation probes in bladder cancer were successfully delineated (Figure 8H), followed by an investigation into their correlation with PRPF19 expression. The present study revealed a significant association between PRPF19 expression and the DNA methylation levels of cg26375147, cg00778920, and cg10733145.
Figure 8 Analysis of PRPF19 methylation levels in bladder cancer. (A) Differences in PRPF19 methylation levels in bladder cancer and normal tissues. (B) Differential expression of PRPF19 methylation levels in different stages of bladder cancer. (C) Differential expression of PRPF19 methylation levels in different N stages of bladder cancer. (D) Methylation probes linked to PRPF19 and their chromosomal distribution. (E) Detailed genomic information of PRPF19. (F) Expression of different methylation probes in bladder cancer tissues and normal samples. (G) Prognostic value of different methylation probes in bladder cancer. (H) Expression distribution of different methylation probes in bladder cancer samples and their correlation with PRPF19 expression. **p < 0.01, ***p < 0.001 and ns not significant.
3.9 Correlation analysis of PRPF19 with senescence associated secretory phenotypes
Cellular senescence has historically been seen as a protective process that safeguards tissues against the development of cancer by permanently halting the cell cycle, restricting proliferation, and minimizing the transfer of unfavorable genetic alterations (25). Senescent cells refer to cancer cells that have ceased proliferating but retain their metabolic activity. These cells may lead to the adverse consequence of a modified cytosolic secretome, which is linked to localized inflammation, alterations in the extracellular environment, and heightened growth factor activity (26). In the first phase of our study, we conducted an investigation on the TCGA BLCA dataset to examine the potential correlation between PRPF19 expression levels and the expression levels of genes associated with SASP. The findings presented in Figure 9A revealed a negative correlation between PRPF19 and CCL16 as well as IGFBP3. The present study aimed to investigate potential disparities in the expression of genes associated with SASP across two distinct groups categorized based on their levels of PRPF19 expression. In addition to EREG, MMP13, CXCL2, IL1B, and PLAU, the results indicated significant differences in the expression levels of the remaining genes between the PRPF19 high expression group and the PRPF19 low expression group, as seen in Figure 9B. The genes associated with SASP were shown to be positively connected with a tumor stemness score. This observation suggests that senescent tumor cells could undergo reprogramming into tumor stem cells through this particular method (27). A robust association was seen between tumor stemness scores and a significant proportion of genes implicated in the secretory phenotype often linked to senescence (Figure 9C). Additionally, an investigation was conducted to examine the impact of PRPF19 on bladder cancer stem cells. The results revealed that cells with elevated levels of PRPF19 expression exhibited a greater stemness score compared to cells with lower levels of PRPF19 expression, as seen in Figure 9D. A single differential gene was identified to have a high correlation with the tumor stemness score and a negative association with PRPF19, as seen in Figure 9E via the use of Venn diagrams. The TCGA bladder cancer dataset revealed a significant correlation between the expression levels of CCL16 and PRPF19 (r=0.166) (Figure 9F). Cancer stem cells (CSCs) are a distinct population of cells inside tumors that possess the ability to undergo unlimited division and generate new tumors. These cells play a critical role in the progression, dissemination, and persistence of malignant tumors, as well as in the acquisition of resistance to cytotoxic treatment (28). The process of inducing a stem cell-like state in cancer cells is associated with an elevated level of resistance to therapy in tumor cells. The establishment of this state resembling stem cells is facilitated by cell-intrinsic proteins that are activated by senescence-associated secretory phenotypes (29). Subsequently, we proceeded to show the distribution of stemness scores and the distribution of bladder cancer samples. The former exhibited a range spanning from low to high values, while the latter depicted the samples in a sorted manner (Figures 9G, H). In summary, the findings of our study demonstrate a significant correlation between the presence of PPRF19 in bladder cancer specimens and the manifestation of a secretory phenotype often associated with senescence. There exists a potential regulatory relationship between PRPF19 and CCL16, which may contribute to the promotion of stemness in bladder cancer cells.
Figure 9 PRPF19 is significantly associated with senescence-associated secretory phenotypes. (A) Correlation analysis of PRPF19 with genes related to senescence-associated secretory phenotypes. (B) PRPF19 expression levels correlate differently with senescence-related secretory phenotypes. (C) Differential senescence-associated secretory phenotype genes positively associated with PRPF19 in the TCGA BLCA cohort and GSE13507 cohort. (D) Kaplan-Meier overall survival curves of CCL16. (E) Scatter plot of PRPF19 and CCL16 correlation. (F) Differences in stemness scores in bladder cancer samples between high and low PRPF19 expressing groups. (G, H) Distribution of stemness scores from low to high in bladder cancer samples.
3.10 Construction of PRPF19-related ceRNA network
Multiple studies have shown the role of ceRNA in the process of tumor formation (30). This study conducted an analysis of the ceRNA network related with PRPF19 in BLCA. A total of 1610 target microRNAs for PRPF19 were first acquired from the MiRWALK website. The ceRNA hypothesis posits a negative link between the expression of microRNAs (miRNAs) and ceRNAs, such as messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs). Among the microRNAs examined, it was shown that only hsa-miR-125a-5p exhibited a substantial association with the prognosis of patients with bladder cancer. Additionally, this particular microRNA had a significant differential expression between bladder cancer tissues and normal tissues, as seen in Figures 10A–C. By using the Genecards online platform, we conducted a comprehensive analysis to determine the subcellular localization of hsa-miR-125a-5p. Our findings indicate that this microRNA mostly localizes inside the nucleus of cells, as seen in Figure 10D. Furthermore, the interaction between lncRNA and miRNA was shown using a Venn diagram. This was achieved by extracting miRNAs that were differently expressed from the ENCORI database. The elevated expression levels of LINC02693 and MIR4435-2HG were shown to be strongly linked to worse prognosis among individuals diagnosed with bladder cancer. Furthermore, their increased expression levels exhibited a negative correlation with hsa-miR-125a-5p. Consequently, these two LncRNAs were selected for further investigation. (Figures 10E–G, I–K). In addition, it was observed that both LINC02693 and MIR4435-2HG exhibited comparable localization inside the nucleus, as seen in Figures 10H, L. The improvement seen may be attributed to the PRPR19-related lncRNA, miRNA, and mRNA triple regulatory networks. These networks include the interaction of several RNAs, including LINC02693 and MIR4435-2HG, which act as ceRNAs to upregulate the expression of PRPF19 by sequestering hsa-miR-125a-5p.
Figure 10 Construction of ceRNA network of PRPF19 in BLCA. (A) Venn diagram of miRNA interacting with PRPF19. (B) Scatter plot of PRPF19 and hsa-miR-125a-5p correlation. (C) Overall survival curves determined using Kaplan and Meier methods for hsa-miR-125a-5p. (D) The subcellular localization of miRNAs that were screened. (E) Venn diagram of LncRNA and hsa-miR-125a-5p interactions. (F) Scatter plot of LINC02693 and hsa-miR-125a-5p correlation. (G) Kaplan-Meier overall survival curves of LINC02693. (H) The subcellular localization of screened LncRNA. (I) Targeted binding sites of hsa-miR-125a-5p to LINC02693 and MIR4435-2HG. (J) Scatter plot of MIR4435-2HG and hsa-miR-125a-5p correlation. (K) MIR4435-2HG’s Kaplan-Meier overall survival curves. (L) The subcellular localization of LncRNA that were screened.
4 Discussion
PRPF19, sometimes referred to as Pso4, SNEV, and NMP200, is situated on chromosome 11q12.2 (31). It is crucial in responding to DNA damage and processing mRNA. Various physiological processes that exert its influence include mRNA splicing and genome maintenance (32). In the context of cellular biology, it has been shown that PRPF19 has the ability to act as a stress response protein, hence enhancing the cellular capacity to withstand oxidative stress (33, 34). Furthermore, it has been found to have a role in mitigating apoptosis and DNA damage (35), as well as impeding cellular senescence (10). PRPF19, apart from its involvement in cellular biological mechanisms, exerts a significant regulatory function in many pathological conditions, such as cancers, by virtue of its multipotency. A research investigation revealed an upregulation of PRPF19 expression in hepatocellular carcinoma, which exhibited a strong association with worse disease prognosis (36). This work aimed to investigate the variations in PRPF19 expression in bladder cancer by using the BLCA dataset from the TCGA database, as well as the GSE13507 and GSE120736 datasets from the GEO database. Our analysis revealed that PRPF19 exhibited elevated expression levels in bladder cancer, with a notable increase seen in higher-grade tumors. This suggests that PRPF19 exerts an oncogenic function in the context of bladder cancer. In a similar vein, the findings of the prognostic analysis substantiated this assertion, as individuals exhibiting elevated levels of PRPF19 expression had a very unfavorable prognosis.
In order to investigate the probable mechanism behind the oncogenic involvement of PRPF19 in bladder cancer, a gene enrichment analysis was conducted. The MHC pathway was discovered by our research team as a potentially significant channel by which PRPF19 may exert regulatory control on the immune microenvironment in cases of bladder cancer. The complex biological interplay between immune cells inside the tumor microenvironment and tumor cells may give rise to diverse immunotherapeutic responses. Our investigation using single-cell RNA sequencing of BLCA tissues shown that PRPF19 exhibited predominant distribution in CD4 Tconv, CD8T, and NK cells. This finding suggests a potential association between PRPF19 and invading immune cells. The prognostic model used in our study utilized the LASSO algorithm and included eight essential immune infiltrating cells. Notably, individuals classified in the high-risk group had significantly worse prognoses compared to those in the low-risk category. The prognostic AUC values at 1, 3, and 5 years were all found to be higher than 0.6 in our model for patients with bladder cancer, as seen by the ROC curves.
There is an increasing body of data indicating that DNA methylation has an influence on the metabolic activities of cancer cells (37). DNA methylation plays a significant role in several metabolic processes, including glycolysis, the methionine cycle, and the tricarboxylic acid (TCA) cycle (38–40). In recent times, there has been a growing awareness of the potential use of DNA methylation status and its associated genes in the identification, management, and prediction of many illnesses (41, 42). In the present investigation, it was shown that there were significant variations in the methylation levels of PRPF19 between bladder cancer and normal bladder tissue. Hence, we conducted a more comprehensive examination of PRPF19 expression in bladder cancer by the utilization of methylation probes that possess predictive significance in this particular malignancy. In our investigation, it was shown that there exists a negative correlation between cg26375147 and PRPF19 in the context of bladder cancer. Additionally, a significant correlation was seen between cg26375147 and the prognosis of individuals diagnosed with bladder cancer. It is plausible that PRPF19 methylation may exert its influence on bladder cancer via the involvement of cg26375147.
The phenomenon of cellular senescence has been extensively acknowledged in the scientific community as a significant mechanism for the prevention of tumor growth (43). In previous studies, bladder cancer cells (44), colorectal cancer cells (45), and breast cancer cells (46) were subjected to cellular senescence, resulting in the induction of persistent cell cycle arrest and subsequent inhibition of growth in these tumor cells. While this strategy first shown efficacy in combating malignancies, it subsequently resulted in the development of tumor resistance as an unintended consequence. Senescent tumor cells do not die outright, but are transformed into tumor stem cells (47). The signaling pathways including p16, p21, and p53 are implicated in both senescence and stemness processes. This suggests that senescence could play a role in gene reprogramming and activation of stemness, hence contributing to tumor development driven by cancer stem cells (48). In this work, an investigation was conducted to examine the correlation between PRPF19 and genes linked to senescence-associated secretory characteristics across several bladder cancer datasets. The gene CCL16 was shown to have the strongest association with PRPF19 in the context of the senescence-associated secretory phenotype. According to recent research findings, it has been shown that the regulation of stemness in breast cancer cells is influenced by CCL16 (49). Furthermore, an examination was conducted to assess the significance of PRPF19 in relation to bladder tumor stem cells. The findings revealed that elevated expression of PRPF19 was associated with a greater stemness score. Hence, it is postulated that PRPF19 potentially modulates the stemness of bladder tumor cells by influencing the expression of the senescence-associated secretory phenotype-associated gene CCL16. The ceRNA network refers to a kind of post-transcriptional control that is facilitated by miRNAs. The task at hand is establishing the connection between the functions of coding and non-coding RNAs. The ceRNA network regulates the production of mRNA, potentially impacting several biological processes and contributing to the development of numerous diseases. This regulation occurs via the competitive binding of lncRNAs or circRNAs to miRNAs (50). In order to ascertain the presence of has-miR-125a-5p, an initial screening was conducted to find miRNAs that exhibited a connection or underwent survival analysis in conjunction with PRPF19. Two other LncRNAs were subjected to the same experimental methodology, leading to the identification of PRPF19, hsa-miR-125a-5p, LINC02693, and MIR4435-2HG as nuclear-localized entities. In summary, the construction of the ceRNA network aimed to elucidate the interaction profile of the BLCA gene. This was achieved by investigating the interactions involving PRPF19, hsa-miR-125a-5p, LINC02693, and MIR4435-2HG. This study analysed the role of PRPF19 in bladder cancer prognosis as well as in the immune microenvironment from several perspectives, but threw some limitations. We only confirmed the expression and prognostic value of PRPF19 in bladder cancer by bladder cancer tissue microarrays, and further experimental analysis of the role of PRPF19 and the immune microenvironment of bladder cancer is necessary, which will be further explored in our subsequent studies.
5 Conclusion
In summary, our study provides evidence that PRPF19 has a substantial impact on the immunological environment of bladder cancer and has the potential to serve as a prognostic indicator in patients afflicted with this condition. Additionally, it was shown that PRPF19 had a significant correlation with both the secretory phenotype often seen in senescence and the stemness characteristics of bladder cancer cells. Finally, it is worth noting that the ceRNA network under investigation has potential therapeutic implications in the context of bladder cancer.
Data availability statement
The original contributions presented in the study are included in the article/supplementary materials. Further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by Shanghai outdo Ethical Committee. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
YW: Conceptualization, Formal Analysis, Writing – original draft. JFW: Data curation, Writing – original draft. JH: Writing – original draft. BJ: Investigation, Writing – review & editing. ZP: Formal Analysis, Investigation, Writing – review & editing. JSW: Methodology, Writing – review & editing. YL: Supervision, Writing – review & editing. MR: Funding acquisition, Resources, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Horizontal Project of Heilongjiang Renxin Medical Assistance Foundation, Heilongjiang Charity Fund General Association, Natural Science Foundation of Heilongjiang Province (No. LH2019H030).
Acknowledgments
The findings that were shown here are wholly or partially based on information that was produced by the TCGA and GEO Network.
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Yang CJ, Tan ZL, Yang JD, Hsu FT, Chiang CH. Fluoxetine inactivates STAT3/NF-κB signaling and promotes sensitivity to cisplatin in bladder cancer. BioMed Pharmacother (2023) 164:114962. doi: 10.1016/j.biopha.2023.114962
2. Zhang GC, Yu XN, Guo HY, Sun JL, Liu ZY, Zhu JM, et al. PRP19 enhances esophageal squamous cell carcinoma progression by reprogramming SREBF1-dependent fatty acid metabolism. Cancer Res (2023) 83(4):521–37. doi: 10.1158/0008-5472.CAN-22-2156
3. Maréchal A, Li JM, Ji XY, Wu CS, Yazinski SA, Nguyen HD, et al. PRP19 transforms into a sensor of RPA-ssDNA after DNA damage and drives ATR activation via a ubiquitin-mediated circuitry. Mol Cell (2014) 53(2):235–46. doi: 10.1016/j.molcel.2013.11.002
4. Mahajan K. hPso4/hPrp19: a critical component of DNA repair and DNA damage checkpoint complexes. Oncogene. (2016) 35(18):2279–86. doi: 10.1038/onc.2015.321
5. Shen J, Li P, Shao X, Yang Y, Liu X, Feng M, et al. The E3 ligase RING1 targets p53 for degradation and promotes cancer cell proliferation and survival. Cancer Res (2018) 78(2):359–71. doi: 10.1158/0008-5472.CAN-17-1805
6. Rape M. Ubiquitylation at the crossroads of development and disease. Nat Rev Mol Cell Biol (2018) 19(1):59–70. doi: 10.1038/nrm.2017.83
7. Zhou R, Chen J, Xu Y, Ye Y, Zhong G, Chen T, et al. PRPF19 facilitates colorectal cancer liver metastasis through activation of the Src-YAP1 pathway via K63-linked ubiquitination of MYL9. Cell Death Dis (2023) 14(4):258. doi: 10.1038/s41419-023-05776-2
8. Cheng L, Xu C, Chen Y, Li J, Wang J, Xiao J, et al. DEC1 represses cardiomyocyte hypertrophy by recruiting PRP19 as an E3 ligase to promote ubiquitination-proteasome-mediated degradation of GATA4. J Mol Cell Cardiol (2022) 169:96–110. doi: 10.1016/j.yjmcc.2022.05.005
9. Cai Y, Chen K, Cheng C, Xu Y, Cheng Q, Xu G, et al. Prp19 is an independent prognostic marker and promotes neuroblastoma metastasis by regulating the hippo-YAP signaling pathway. Front Oncol (2020) 10:575366. doi: 10.3389/fonc.2020.575366
10. Yano K, Takahashi RU, Shiotani B, Abe J, Shidooka T, Sudo Y, et al. PRPF19 regulates p53-dependent cellular senescence by modulating alternative splicing of MDM4 mRNA. J Biol Chem (2021) 297(1):100882. doi: 10.1016/j.jbc.2021.100882
11. He Y, Huang C, Cai K, Liu P, Chen X, Xu YI, et al. PRPF19 promotes tongue cancer growth and chemoradiotherapy resistance. Acta Biochim Biophys Sin (Shanghai) (2021) 53(7):893–902. doi: 10.1093/abbs/gmab059
12. Xu Z, Pei C, Cheng H, Song K, Yang J, Li Y, et al. Comprehensive analysis of FOXM1 immune infiltrates, m6a, glycolysis and ceRNA network in human hepatocellular carcinoma. Front Immunol (2023) 14:1138524. doi: 10.3389/fimmu.2023.1138524
13. Yang H. Comprehensive analysis of the expression and clinical significance of a ferroptosis-related genome in ovarian serous cystadenocarcinoma: A study based on TCGA data. Oncologie 24(4):835–63. doi: 10.32604/oncologie.2022.026447
14. Wang Y, Zhu H, Wang X. Prognosis and immune infiltration analysis of endoplasmic reticulum stress-related genes in bladder urothelial carcinoma. Front Genet (2022) 13:965100. doi: 10.3389/fgene.2022.965100
15. Han Y, Wang Y, Dong X, Sun D, Liu Z, Yue J, et al. TISCH2: expanded datasets and new tools for single-cell transcriptome analyses of the tumor microenvironment. Nucleic Acids Res (2023) 51(D1):D1425–31. doi: 10.1093/nar/gkac959
16. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol (2017) 18(1):220. doi: 10.1186/s13059-017-1349-1
17. Chandrashekar DS, Karthikeyan SK, Korla PK, Patel H, Shovon AR, Athar M, et al. UALCAN: An update to the integrated cancer data analysis platform. Neoplasia. (2022) 25:18–27. doi: 10.1016/j.neo.2022.01.001
18. Li Y, Ge D, Lu C. The SMART App: an interactive web application for comprehensive DNA methylation analysis and visualization. Epigenet Chromatin (2019) 12(1):71. doi: 10.1186/s13072-019-0316-3
19. Malta TM, Sokolov A, Gentles AJ, Burzykowski T, Poisson L, Weinstein JN, et al. Machine learning identifies stemness features associated with oncogenic dedifferentiation. Cell. (2018) 173(2):338–354.e15. doi: 10.1016/j.cell.2018.03.034
20. Sticht C, de la Torre C, Parveen A, Gretz N. miRWalk: An online resource for prediction of microRNA binding sites. PloS One (2018) 13(10):e0206239. doi: 10.1371/journal.pone.0206239
21. Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res (2014) 42(Database issue):D92–7. doi: 10.1093/nar/gkt1248
22. Cruz FM, Chan A, Rock KL. Pathways of MHC I cross-presentation of exogenous antigens. Semin Immunol (2023) 66:101729. doi: 10.1016/j.smim.2023.101729
23. Cai W, Zhou D, Wu W, Tan WL, Wang J, Zhou C, et al. MHC class II restricted neoantigen peptides predicted by clonal mutation analysis in lung adenocarcinoma patients: implications on prognostic immunological biomarker and vaccine design. BMC Genomics (2018) 19(1):582. doi: 10.1186/s12864-018-4958-5
24. Li Y, McGrail DJ, Xu J, Mills GB, Sahni N, Yi S. Gene regulatory network perturbation by genetic and epigenetic variation. Trends Biochem Sci (2018) 43(8):576–92. doi: 10.1016/j.tibs.2018.05.002
25. Lucas V, Cavadas C, Aveleira CA. Cellular senescence: from mechanisms to current biomarkers and senotherapies. Pharmacol Rev (2023) 75(4):675–713. doi: 10.1124/pharmrev.122.000622
26. Han X, Lei Q, Xie J, Liu H, Li J, Zhang X, et al. Potential regulators of the senescence-associated secretory phenotype during senescence and aging. J Gerontol A Biol Sci Med Sci (2022) 77(11):2207–18. doi: 10.1093/gerona/glac097
27. Lagunas AM, Francis M, Maniar NB, Nikolova G, Wu J, Crowe DL. Paracrine interaction of cancer stem cell populations is regulated by the senescence-associated secretory phenotype (SASP). Mol Cancer Res (2019) 17(7):1480–92. doi: 10.1158/1541-7786.MCR-18-1356
28. Pan Y, Gu Z, Lyu Y, Yang Y, Chung M, Pan X, et al. Link between senescence and cell fate: senescence-associated secretory phenotype and its effects on stem cell fate transition. Rejuvenation Res (2022) 25(4):160–72. doi: 10.1089/rej.2022.0021
29. Heft Neal ME, Brenner JC, Prince MEP, Chinn SB. Advancement in cancer stem cell biology and precision medicine-review article head and neck cancer stem cell plasticity and the tumor microenvironment. Front Cell Dev Biol (2022) 9:660210. doi: 10.3389/fcell.2021.660210
30. Wang L, Ma X, Chen Y, Zhang J, Zhang J, Wang W, et al. MiR-145-5p suppresses hepatocellular carcinoma progression by targeting ABHD17C. Oncologie 24(4):897–912. doi: 10.32604/oncologie.2022.025693
31. Benjamin AB, Zhou X, Isaac O, Zhao H, Song Y, Chi X, et al. PRP19 upregulation inhibits cell proliferation in lung adenocarcinomas by p21-mediated induction of cell cycle arrest. BioMed Pharmacother (2014) 68(4):463–70. doi: 10.1016/j.biopha.2014.03.006
32. de Moura TR, Mozaffari-Jovin S, Szabó CZK, Schmitzová J, Dybkov O, Cretu C, et al. Prp19/pso4 is an autoinhibited ubiquitin ligase activated by stepwise assembly of three splicing factors. Mol Cell (2018) 69(6):979–992.e6. doi: 10.1016/j.molcel.2018.02.022
33. Chanarat S, Seizl M, Strässer K. The Prp19 complex is a novel transcription elongation factor required for TREX occupancy at transcribed genes. Genes Dev (2011) 25(11):1147–58. doi: 10.1101/gad.623411
34. Dellago H, Khan A, Nussbacher M, Gstraunthaler A, Lämmermann I, Schosserer M, et al. ATM-dependent phosphorylation of SNEVhPrp19/hPso4 is involved in extending cellular life span and suppression of apoptosis. Aging (Albany NY) (2012) 4(4):290–304. doi: 10.18632/aging.100452
35. Ohi MD, Vander Kooi CW, Rosenberg JA, Ren L, Hirsch JP, Chazin WJ, et al. Structural and functional analysis of essential pre-mRNA splicing factor Prp19p. Mol Cell Biol (2005) 25(1):451–60. doi: 10.1128/MCB.25.1.451-460.2005
36. Yang M, Qiu Y, Yang Y, Wang W. An integrated analysis of the identified PRPF19 as an onco-immunological biomarker encompassing the tumor microenvironment, disease progression, and prognoses in hepatocellular carcinoma. Front Cell Dev Biol (2022) 10:840010. doi: 10.3389/fcell.2022.840010
37. Chen C, Wang Z, Qin Y. Connections between metabolism and epigenetics: mechanisms and novel anti-cancer strategy. Front Pharmacol (2022) 13:935536. doi: 10.3389/fphar.2022.935536
38. Zeng JD, Wu WKK, Wang HY, Li XX. Serine and one-carbon metabolism, a bridge that links mTOR signaling and DNA methylation in cancer. Pharmacol Res (2019) 149:104352. doi: 10.1016/j.phrs.2019.104352
39. Rashkovan M, Ferrando A. Metabolic dependencies and vulnerabilities in leukemia. Genes Dev (2019) 33(21-22):1460–74. doi: 10.1101/gad.326470.119
40. Shanmugam M, McBrayer SK, Rosen ST. Targeting the Warburg effect in hematological Malignancies: from PET to therapy. Curr Opin Oncol (2009) 21(6):531–6. doi: 10.1097/CCO.0b013e32832f57ec
41. Papanicolau-Sengos A, Aldape K. DNA methylation profiling: an emerging paradigm for cancer diagnosis. Annu Rev Pathol (2022) 17:295–321. doi: 10.1146/annurev-pathol-042220-022304
42. Cai G, Cai M, Feng Z, Liu R, Liang L, Zhou P, et al. A multilocus blood-based assay targeting circulating tumor DNA methylation enables early detection and early relapse prediction of colorectal cancer. Gastroenterology. (2021) 161(6):2053–2056.e2. doi: 10.1053/j.gastro.2021.08.054
43. Bai Z, Zhou Y, Ye X, Li Y, Peng Y, Guan Q, et al. Survivin suppression heightens BZML-induced mitotic catastrophe to overcome multidrug resistance by removing therapy-induced senescent A549/Taxol cells. Biochim Biophys Acta Mol Cell Res (2022) 1869(2):119174. doi: 10.1016/j.bbamcr.2021.119174
44. Jia M, Su B, Mo L, Qiu W, Ying J, Lin P, et al. Circadian clock protein CRY1 prevents paclitaxel−induced senescence of bladder cancer cells by promoting p53 degradation. Oncol Rep (2021) 45(3):1033–43. doi: 10.3892/or.2020.7914
45. Wang R, Sun L, Xia S, Wu H, Ma Y, Zhan S, et al. B7-H3 suppresses doxorubicin-induced senescence-like growth arrest in colorectal cancer through the AKT/TM4SF1/SIRT1 pathway. Cell Death Dis (2021) 12(5):453. doi: 10.1038/s41419-021-03736-2
46. Bojko A, Staniak K, Czarnecka-Herok J, Sunderland P, Dudkowska M, Śliwińska MA, et al. Improved autophagic flux in escapers from doxorubicin-induced senescence/polyploidy of breast cancer cells. Int J Mol Sci (2020) 21(17):6084. doi: 10.3390/ijms21176084
47. Jiang X, Li W, Ge L, Lu M. Mesenchymal stem cell senescence during aging:From mechanisms to rejuvenation strategies. Aging Dis (2023) 14(5):1651–76. doi: 10.14336/AD.2023.0208
48. Zhang DY, Monteiro MJ, Liu JP, Gu WY. Mechanisms of cancer stem cell senescence: Current understanding and future perspectives. Clin Exp Pharmacol Physiol (2021) 48(9):1185–202. doi: 10.1111/1440-1681.13528
49. Shen W, Zhang X, Tang J, Zhang Z, Du R, Luo D, et al. CCL16 maintains stem cell-like properties in breast cancer by activating CCR2/GSK3β/β-catenin/OCT4 axis. Theranostics. (2021) 11(5):2297–317. doi: 10.7150/thno.51000
Keywords: prognosis, immune infiltration, senescence, DNA methylation, ceRNA
Citation: Wang Y, Wang J, He J, Ji B, Pang Z, Wang J, Liu Y and Ren M (2023) Comprehensive analysis of PRPF19 immune infiltrates, DNA methylation, senescence-associated secretory phenotype and ceRNA network in bladder cancer. Front. Immunol. 14:1289198. doi: 10.3389/fimmu.2023.1289198
Received: 05 September 2023; Accepted: 25 October 2023;
Published: 06 November 2023.
Edited by:
Hongda Liu, Nanjing Medical University, ChinaReviewed by:
Ziheng Wang, University of Macau, ChinaYunfei Liu, Ludwig Maximilian University of Munich, Germany
Copyright © 2023 Wang, Wang, He, Ji, Pang, Wang, Liu and Ren. 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: MingHua Ren, cmVubWluZ2h1YTE5NzJAMTYzLmNvbQ==