Downregulation of PTCD1 in Bladder Urothelial Carcinoma Predicts Poor Prognosis and Levels of Immune Infiltration

Zhou, Zhongbao; Zhou, Yongjian; Zhou, Xin; Huang, Yongjin; Cui, Yuanshan; Zhang, Yong

doi:https://doi.org/10.1155/2022/1146186

Journal of Oncology

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 1146186 | https://doi.org/10.1155/2022/1146186

Downregulation of PTCD1 in Bladder Urothelial Carcinoma Predicts Poor Prognosis and Levels of Immune Infiltration

Zhongbao Zhou,¹Yongjian Zhou,¹Xin Zhou,¹Yongjin Huang,¹Yuanshan Cui,²and Yong Zhang¹

Academic Editor: Zhen-Jian Zhuo

Received31 Jan 2022

Revised27 May 2022

Accepted06 Jun 2022

Published28 Jun 2022

Abstract

Pentatricopeptide repeat domain 1 (PTCD1) was reported to regulate mitochondrial metabolism and oxidative phosphorylation. However, the effect and mechanism of PTCD1 in the development of bladder urothelial carcinoma (BLCA) remain unclear. The databases from The Cancer Genome Atlas (TCGA) and Human Protein Atlas (HPA) were used to analyze the expression changes, clinical features, and prognostic values of PTCD1. A nomogram was built to predict the prognostic outcomes of BLCA cases. The potential genes interacting with PTCD1 were explored by Weighted Gene Coexpression Network Analysis (WGCNA). The estimation of associations between PTCD1 and tumor mutations, tumor immunities, and m6A methylations was performed. The study found that the gradual decrease of PTCD1 expression was observed with the increase of stage and grade. Low PTCD1 expression was greatly correlated with higher pathological stage, N stage, and poor prognosis in TCGA cohorts; interestingly, low-grade BLCA cases all exhibited high expression of PTCD1. HPA database analysis implied that the expression of PTCD1 protein in BLCA was lower than that in normal bladder tissue, and the protein expression of PTCD1 in high-grade BLCA was lower than that in low-grade BLCA. Multivariate Cox regression analysis indicated that PTCD1 may serve as an independent factor influencing prognosis of BLCA. Mechanistically, PTCD1 played a regulatory role in BLCA progression through multiple tumor-related pathways containing PI3K-Akt signaling, ECM-receptor interaction, oxidative phosphorylation, and extracellular matrix organization. WGCNA reported that PTCD1 had a strong positive correlation with POLR2J, ZNHT1, ATP5MF, PDAP1, BUD31, and COPS6. Besides, the mRNA expression of PTCD1 was negatively associated with immune cells’ infiltrations, immune functions, and checkpoints, especially with some m6A methylation regulators in BLCA. In sum, downregulation of PTCD1 expression may be involved in the development of BLCA and remarkably correlated with poor prognosis. Meantime, it showed an influence in immune cell infiltration and may serve as an agreeable prognostic indicator in BLCA.

1. Introduction

Bladder urothelial carcinoma (BLCA) is considered as one of the most invasive cancers in the world [1]. It is reported that the five-year overall survival (OS) in BLCA is approximately 80% [2]. Furthermore, it is estimated that more than 80 thousand new cases and 17 thousand deaths annually occur in the United States alone [3]. Bladder tumors are assorted into high and low grade according to histomorphological features following the description from the World Health Organization [4]. About eighty precent of BLCA cases express non-muscle-invasive bladder cancer (NMIBC) in the diagnosis, while the rest cases show muscle-invasive bladder cancer (MIBC) or even distant metastasis [5]. NMIBC generally does not pose threats to the survival of patients and has better prognosis for effective treatment options [6]. However, T1 tumors in NMIBC act as an invasive subtype with forty percent relapse and fifteen percent progression to MIBC within five years [7]. MIBC has clinical aggression with a rapid development of metastases to bone, brain, liver, lungs, and lymph nodes, which usually poses a fatal threat to patients [8]. Despite significant progress in revealing the molecular mechanisms and surgery techniques of BLCA in recent years, there is no significant decrease for the mortality rate [9–11]. Therefore, exploring molecular biomarker is helpful to the early disease diagnosis, progression, and prediction of prognostic and individualized therapy for specific patient [12].

Pentatricopeptide repeat (PPR) proteins are a large class of RNA-binding proteins characterized by a typical thirty-five residue repeat motif [13]. Only seven PPR proteins are found in human beings, all of which are located in mitochondria and involved in various cell procedures, notably in RNA supersession [14, 15]. The discovery of essential role of PPR proteins in mitochondria gene expression and energy metabolism underscores its biological importance [15]. For example, mitochondrial RNA polymerase regulates the transcription of mitochondrial mRNA (mt-mRNA) [16]. PPR domain 1 (PTCD1), PPR domain 2 (PTCD2), and only protein RNase P catalytic subunit are related to the process of mt-mRNA [17, 18]. Leucine-rich pentatricopeptide repeat motif-containing protein (LRPPRC) is involved in the evolution and stability regulations of mt-mRNA, and PPR domain 3 (PTCD3) is related in protein composite in mitochondria [19, 20]. In recent years, an increasing number of PPRs have been identified with important regulation effects of BLCA carcinogenesis. For instance, LRPPRC promotes BLCA tumorigenesis by regulating redox homeostasis through the circANKHD1/FOXM1 [15]. Black et al. reported that FOXM1 translocated into mitochondria and inhibited mitochondria respiration and oxidative phosphorylation by enhancing PTCD1 and may be a promising molecular target [21]. As above shown, PTCD1 is very likely involved in formation and development of cancer. However, the effect of PTCD1 in BLCA is still not well clarified.

In this research, RNA-seq data obtained from The Cancer Genome Atlas (TCGA) database was employed to comprehensively assess the diversities in PTCD1 expression and their relationships with patient outcomes. The Human Protein Atlas (HPA) database was used to explore the protein expression of PTCD1 in normal bladder and tumor tissues. Moreover, Gene Set Enrichment Analysis (GSEA) was conducted on the high- and low-PTCD1 expression groups, respectively, to show the latent mechanisms. Next, we explored the potential genes interacting with PTCD1 by Weighted Gene Coexpression Network Analysis (WGCNA). Finally, the association between the expression of PTCD1 and extent of immune cell infiltration was investigated to examine the probable mechanism of PTCD1-induced tumorigenesis and development.

2. Materials and Methods

2.1. Gathering of PTCD1 Expression Data from TCGA and HPA Databases

We followed the methods of Dr. Sun et al. [22]. The expression of PTCD1 combined clinical data (containing 414 BLCA cases) were gathered from TCGA database (https://portal.gdc.cancer.gov/). The clinical data mainly included age, sex, grade, pathological stage, TNM stage, survival condition, and survival time. The exclusion criteria were as follows: (a) The pathological diagnosis did not meet BLCA; (b) The PTCD1 expression profile and clinical data were incomplete. We used the HPA database (https://www.proteinatlas.org/) to explore the protein level of PTCD1 in tissues. Normalization was performed by log 2 (FPKM + 1) to transform expression data for data analysis.

2.2. Survival Analysis

Survival data were analyzed using R software and R package “survival” and “survminer.” The threshold was identified by the best cut-off expression of PTCD1 in tumor samples and the queues were classified into low- and high-expression groups. The correlations between PTCD1 expression and patient outcomes, including OS, grade, and stage, were examined. The ENCORI database (https://starbase.sysu.edu.cn/) was applied to verify the prognosis of PTCD1 in BLCA.

2.3. Diagnostic Value and Construction of Prognostic Nomogram

The univariate and multivariate Cox regression analyses were used to assess whether PTCD1 expression served as an independent prognostic factor in BLCA using “glmnet” package in R software. Based on the results of multivariate Cox regression analysis, we performed a prognostic nomogram to quantify the prediction of prognosis for BLCA patients. The concordance index (C-index) and calibration curve were employed to assess the reliability and accurateness of nomogram. In addition, we evaluated the effects of PTCD1 expression on the clinical prognosis in BLCA patients.

2.4. Differentially Expressed Genes and Functional Analysis

Background correction and normalization were performed on public microarray data, and R package “limma” was applied to identify differential expression genes (DEGs) between low- and high-PTCD1 expression groups. The values were corrected for multiple test correction by false discovery rate (FDR). DEGs were included in the Gene Ontology (GO) term enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using the R package “clusterProfiler,” which were recruited when |log2FC| value is more than 1 and FDR less than 0.001. FDR q-value less than 0.05 and P adjustment less than 0.05 were greatly enriched for gene sets.

2.5. WGCNA and Module Identification

WGCNA was employed to determine the gene module closely related with PTCD1 expression. To identify gene clusters that were likely to be highly coexpressed, a coexpression network was constructed using profiles of DEGs dataset with R package “WGCNA.” First, cluster analysis was performed on tumor samples using the “hclust” function to check and remove outliers. Second, Pearson’s correlation between every pair of extracted gene was analyzed to generate an adjacency matrix. Third, a soft-threshold parameter performance value (β) was built, which can accentuate the strong correlation of genes while penalizing the low correlation to guarantee the construction of a scale-free network. Furthermore, based on TOM-based dissimilarity (1-TOM) with a dendrogram of more than 30 genes, hierarchical clustering analysis was conducted to examine modules with genes of similar expression profiles. A cutoff (less than 0.25) was then chosen to incorporate the similar modules to ensure the outcomes trustier. “DynamicTreeCut” algorithm was employed to construct networks and detect consensus modules. The relationships between modules and BLCA were computed by the module-trait associations with WGCNA. The module eigengene (ME, representing the gene expression profile of a module) was considered as the first major component of a given module. Finally, clinically significant modules were identified by computing the association between clinical features and MEs. By selecting modules where PTCD1 was located, the network construction and core gene screening were carried out.

2.6. Tumor Mutation Profile Analysis

We used the somatic mutation data obtained from the TCGA database and processed in R software using “maftools” package to identify the somatic variants and visualize somatic landscape. Mutant signatures from BLCA tissues were described with “Somatic Signatures” package in R software.

2.7. Relationships between PTCD1 Expressions and Immune Features

Immune cell infiltration data were gathered from TIMER website (https://timer.cistrome.org/). “GSVA,” “limma,” and “GSEABase” packages were employed to analyze and visualize the data. The relationships between different PTCD1 expression groups and immune infiltrations were examined by Wilcoxon rank-sum test. Immune function analysis of BLCA was conducted using the single-sample GSEA function with the R package “GSVA.” Moreover, the association between PTCD1 expression and immune checkpoint biomarkers was shown by Wilcoxon rank-sum test.

2.8. Correlation between PTCD1 Expression and m6A RNA Methylation Regulators

The matrixes of expressions were taken from transcriptome profiles datasets, containing regulators on writers (VIRMA, METTL3, METTL14, WTAP, RBM15, RBM15B, METTL16, ZC3H13, and PCIF1), readers (TRMT112, ZCCHC4, NUDT21, CPSF6, CBLL1, SETD2, HNRNPC, HNRNPG, HNRNPA2B1, IGF2BP1, IGF2BP2, IGF2BP3, YTHDC1, YTHDF1, YTHDF2, YTHDF3, YTHDC2, SRSF3, SRSF10, XRN1, FMR1, NXF1, and PRRC2A), and erasers (FTO, ALKBH5, and ALKBH3). The analysis of m6A RNA methylation regulators was conducted between the two differential PTCD1 expression groups using “limma” package with the Wilcoxon rank-sum test. P less than 0.05 was used as the significant significance.

2.9. Statistics Analysis

All statistical analyses were performed using R software (version 4.0.2). Wilcoxon rank-sum test was used to evaluate the difference of PTCD1 expression among multiple stratified clinical indicators. Categorical variables were presented as proportions, and continuous variables were expressed as mean ± standard deviation (SD). The association of clinicopathological variables in BLCA patients between high- and low-PTCD1 cohorts was subjected to a chi-square test. The survival curves were generated via the log-rank test for Kaplan–Meier analysis. The univariate cox regression model was used to analyze the effects of individual variables on survival, and the multivariate cox regression model was used to confirm independent factors associated with survival. Wilcoxon rank-sum test was applied to analyze the relationships between PTCD1 expression and molecular functions. The difference was considered significant when .

3. Results

3.1. The Expression Characteristics of PTCD1 in BLCA

The stratified analysis according to clinicopathological characteristics was conducted to evaluate the difference of PTCD1 expression, including age (less than 70 years old and more than 70 years old), sex (male and female), stage (stage I, stage II, stage III, and stage IV), and grade (high grade and low grade), which showed that the expression of PTCD1 gradually reduced with the growth of stage (), and the expression of PTCD1 in low-grade bladder cancer was importantly higher than that in high-grade bladder cancer (), but age () or gender () had no effect on the expression of PTCD1 (Figures 1(a)–1(d)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

Figure 1

The PTCD1 expression level and survival analysis in BLCA. (a) The PTCD1 expressions of patients with BLCA according to different clinical characteristics including age (a), gender (b), stage (c), and grade (d). The PTCD1 expression (e), OS (f), PFS (g), and DFS (h) of stages I-II versus stages III-IV in the TCGA-BLCA dataset. The K-M curve between low-PTCD1 group and high-PTCD1 group in the TCGA database (i) and AUC curve related to OS (j). The survival curve based on the median expression of PTCD1 as the threshold in ENCORI database (k). Protein expression of PTCD1 in normal bladder tissue (l), low-grade (m) and high-grade (n) bladder cancer. ; ; and .

3.2. Relationship between PTCD1 Expression and Prognosis

We divided the patients in TCGA-BLCA data set into G1 (stages I-II) and G2 (stages III-IV) and found that the expression of PTCD1 in G2 group was importantly lower than that in G1 group (Figure 1(e)). Great differences were shown in OS (), progression-free survival (PFS; ), and disease-free survival (DFS; ) between the two groups (Figures 1(f)–1(h)). Next, TCGA-BLCA cohorts were classified into low- or high-PTCD1 groups based on the best cut-off expression of PTCD1 (FPKM = 1.128054) in tumor samples as the threshold. The features of BLCA patients were available in Table 1, and we found that the expression of PTCD1 was correlated with pathologic stage (), N stage (), and survival status () but was not associated with age (), sex (), grade (), T stage (), and M stage (). Taken together, these results elucidated that the low expression of PTCD1 was greatly associated with advanced N stage and pathologic stage in BLCA; impressively, low-grade BLCA cases all exhibited high expression of PTCD1. The survival curve revealed that the survival ability of cases in the high-PTCD1 group was remarkably higher than that of the low-PTCD1 group () (Figure 1(i)). The area under curve (AUC) for OS reached 0.591, implying a good predictive value (Figure 1(j)). The ENCORI database also further confirmed the role of PTCD1 in the survival and prognosis of BLCA (; Figure 1(k)). HPA database analysis implied that the expression of PTCD1 protein in BLCA was lower than that in normal bladder tissue, and the protein expression of PTCD1 in high-grade cases was lower than that in low-grade cases (Figures 1(l)–1(n)). Sankey diagram showed the interrelation between TNM stage, PTCD1 expression, and survival status, which showed that most patients with high expression of PTCD1 were in alive status (Figure 2).

3.3. Construction and Evaluation of Nomogram

The univariate and multivariate Cox regression analysis showed that PTCD1 expression [1; hazard ratio (HR) = 0.715] and age (; HR = 1.035) were independent risk factors for survival prognosis of BLCA (Figures 3(a)-3(b)). The PTCD1 expression and patient’s age were applied to construct a nomogram using the “rms” package in R software to predict 1-, 3-, and 5-year OS of BLCA patients (Figure 3(c)). The higher the PTCD1 expression of tissues, the better the prognosis of patients. The calibration curves were used to observe whether the actual prognostic value was consistent with the predicted value of nomogram, and it was found that the calibration curve of 1-, 3-, 5-year OS was almost consistent with the nomogram (Figure 3(d)).

(a)

(b)

(c)

(d)

3.4. Difference Analysis and Functional Analysis Based on Low- and High-PTCD1 Expression

Finally, 1,445 DEGs between low-PTCD1 group and high-PTCD1 group were determined, which were expressed in the Supplementary File 1. Next, we used the GSEA analysis to investigate the molecular mechanisms linked to the DEGs (Supplementary File 2). The highly abundant GO term in the biological process (BP) category was “extracellular matrix organization,” “extracellular structure organization,” “neutrophil mediated immunity,” and “cell-substrate adhesion.” Remarkably abundant GO terms associated with the cellular component (CC) category contained “cell-substrate junction,” “focal adhesion,” and “collagen-containing extracellular matrix.” In the molecular function (MF) category, DEGs were most abundant in the terms “actin binding,” “cadherin binding,” and “extracellular matrix structural constituent” (Figures 4(a) and 4(b)). Furthermore, KEGG pathway analysis indicated that these genes were especially correlated with PI3K-Akt signaling, focal adhesion, proteoglycans in cancer, ECM-receptor interaction, and oxidative phosphorylation (Figures 4(c) and 4(d)).

(a)

(b)

(c)

(d)

3.5. WGCNA and Module Identification

To identify the module associated with PTCD1 expression, coexpression analysis was conducted to build the gene coexpression networks from DEGs dataset. Pearson’s correlation and average linkage methods were used to cluster tumor samples of DEGs dataset (Figure 5(a)). No abnormal samples were detected or rejected. Optimal β = 6 (scale-free R₂ = 0.9) was chosen to make sure to construct scale-free networks in DEGs dataset (Figure 5(b)). A cutoff of 0.25 and minimal module size of 30 and six modules in DEGs dataset (Figure 5(c)) were used to retain for consecutive analyses (gray modules indicate no assignment to any cluster). The gene interaction network of all modules was shown in Figure 5(d). We then examined the relationship between modules and traits by the heatmap (Figure 5(e)). Among them, the turquoise module including 290 DEGs was where PTCD1 was located, which were positively associated with the blue module (Figure 5(f)).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5

Identification of the module related with PTCD1 in DEGs dataset. (a) Clustering dendrograms of samples as well as traits; (b) the left panel shows the scale-free fitting indices for various soft-thresholding powers (β); (c) cluster dendrogram of coexpression network modules based on the 1-TOM matrix; (d) the gene interaction network of all modules; (e) heatmap of the correlation between module eigengenes and traits of BLCA; and (f) correlation heatmap between each module. Each module represents a cluster of corelated genes and was assigned a unique color.

3.6. Coexpression Network Construction and Linear Correlation Analysis

We selected DEGs directly related to PTCD1 in the turquoise module to construct the expression network by Cytoscape v3.6.0 software (Figure 6(a)). The dot heatmap showed the strength of interaction of each gene (Figure 6(b)). We determined the number of adjacent nodes for each gene and visualized the associated role of the top six core genes with PTCD1 (Figure 6(c)), which identified that PTCD1 had a higher coexpression relationship with POLR2J (R = 0.59), ZNHT1 (R = 0.54), ATP5MF (R = 0.68), PDAP1(R = 0.55), BUD31 (R = 0.6), and COPS6 (R = 0.5), suggesting their potential interaction.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

3.7. The Correlation between PTCD1 Expression and Somatic Mutation Load

By analyzing somatic mutation load data, we found that only 7 BLCA patients had PTCD1 mutation. The further analysis of relationship between somatic mutation load and PTCD1 expression in 7 samples was performed, indicating dominant type of variation classification was missense mutation (Figure 7(a)). A detailed comparison showed that single-nucleotide variant (SNV) was more common than insertions or deletions (Figure 7(b)). Furthermore, main SNV class showed the results of C > T and C > G (Figure 7(b)). Moreover, we ranked the top ten mutated genes in seven cases with calculated percentages (Figure 7(b)).

(a)

(b)

Figure 7

Mutation feature of PTCD1 in BLCA. The alteration frequency with mutation type and mutation site of PTCD1 (a). Distribution of frequently mutated genes in Patients with PTCD1 mutation. The upper bar plot shows the tumor mutation burden (TMB) for each patient, whereas the left bar plot indicates the gene mutation frequency in different groups (a). The top 10 BLCA-correlated mutation genes in patients with PTCD1 mutation and the mutation frequency, variant classification, variant type, and SNV class of the mutated genes in high- and low-PTCD1 subgroups (b).

3.8. The Relationship between PTCD1 Expression and Immune Cell Infiltration

We summed up the results of 405 BLCA cases computed with various algorithm and compared all the immune cell subtypes in two groups. The proportion of partial infiltrating cell subtypes was obviously different between the two groups, among which mostly myeloid dendritic cell, T cell CD8⁺, neutrophil, macrophage M1, cancer-associated fibroblast, monocyte, myeloid dendritic cell, and macrophage M2 had a higher infiltration scale in the low-PTCD1 group, while T cell regulatory (Tregs), T cell follicular helper, and other types had a lower proportion (Figure 8(a), all ).

(a)

(b)

(c)

The immune functions of high-PTCD1 group and low-PTCD1 group were analyzed, respectively, using “GSVA” and “GSEABase” package in R software, and it was found that almost immune-related functions (APC coinhibition, APC costimulation, CCR, checkpoint, cytolytic activity, HLA, inflammation promoting, MHC class I, para-inflammation, T cell coinhibition, T cell costimulation, and Type I IFN response) in the low-PTCD1 group were markedly activated (Figure 8(b), all ), indicating a significant change of immunophenotype in the low-PTCD1 group. We further explored the expressions of immune checkpoint-related biomarkers in two groups and found that some markers (LAIR1, IDO1, PDCD1, TNFRSF8, PDCD1LG2, CD86, CD28, CD44, CD276, TNFSF14, TIGIT, HAVCR2, NRP1, LAG3, TNFSF4, BTLA, CD80, CD274, ICOS, IDO2, CTLA4, CD70, CD200, TNFRSF9, and CD48) in the low-PTCD1 group were upregulated, and some markers (TMIGD2, TNFRSF14, and TNFRSF25) were downregulated, presenting the presence of immunosuppressive and exhausted phenotypes in the low-PTCD1 group (Figure 8(c), all ). According to the above analysis, we found that two groups had a remarkable different pattern of immune infiltration, which might lead to distinct survival benefits.

3.9. The Association between PTCD1 Expression and m6A Regulatory Genes

The expressions of m6A-related genes in two groups were analyzed and found that IGF2BP2, XRN1, ALKBH5, ZC3H13, IGF2BP3, VIRMA, YTHDF3, FMR1, and FTO in the low-PTCD1 group were upregulated, and YTHDC1, RBM15 B, YTHDF1, METTL3, TRMT112, CBLL1, and NXF1 were downregulated, indicating some genetic epigenetic changes in the low-PTCD1 group (Figure 9, all less than 0.05). Based on the above results, we found that patients had significantly different m6A patterns in two groups, which may lead to different survival outcomes.

4. Discussion

PTCD1 is a mitochondrial matrix protein containing eight PPR domains [17]. Different from other PPR domain proteins, PTCD1 is a low-abundance protein correlated with leucine tRNAs and precursor RNAs that contain leucine tRNAs [17]. Previous studies showed that PTCD1 is straightly linked to and exhausted leucine mitochondrial tRNAs, restraining the translation of mitochondrial-encoded proteins and inhibiting oxidative phosphorylation in the mitochondria [17, 23]. Its abnormal expression may lead to many diseases, such as metabolic dysfunction, obesity, and neoplasms [14, 21, 24–26]. However, the implication of its expression in prognosis and diagnosis of patients with BLCA remains unclear. In this research, we attempted to address the role of PTCD1 in BLCA for the first time.

First, we analyzed the association between the transcriptional level of PTCD1 in BLCA and clinicopathological characteristics through TCGA database. The results revealed that PTCD1 expression in bladder cancer tissues gradually reduced with the growth of stage and grade. HPA database analysis found that the expression of PTCD1 protein in BLCA was lower than that in normal bladder tissue, which in high-grade tumors was lower than that in low-grade tumors. Kaplan–Meier curves indicated that downregulation of PTCD1 led to poor prognosis of BLCA, which was reflected in the late clinical features of cancer pathological staging and lymph node metastasis. Multivariate Cox regression analysis further verified that low-PTCD1 expression may serve as an independent prognostic factor for OS in patients with BLCA. These results suggest that PTCD1 exerts a tumor-suppressive role in BLCA. In this research, we noticed that M stage indicated the distant metastasis, but there was no effect on the expression of PTCD1. Combined with follow-up analyses, PTCD1 may influence cancer progression mainly by regulation of the immune microenvironment, but this might not enhance distant metastasis of tumor cells. When analyzing the basic characteristics of patients, all patients in the low-expression group belong to high-grade tumors, which may indicate that PTCD1 expression is a factor in predicting tumor grade.

The GO analysis indicated that PTCD1-related DEGs was highly abundant in some terms including “extracellular matrix organization,” “neutrophil degranulation,” “cell-matrix adhesion,” and “cell-substrate junction.” Genetic and epigenetic alterations might cause changed molecular pathways concerned tumor procession and, afterwards, force bladder cancer cells towards epithelial-to-mesenchymal transition (EMT) [27, 28]. EMT confers complex reprogramming of sessile, nonmotile urothelial cells with loss of apical polarity and increased migratory capacity [29–31]. Molecular regulators regulate the transition from epithelium to mesenchyme, enabling some urothelial cancer cells to acquire mesenchymal features with the ability to self-renew, evade immune mechanisms, and infiltrate the surrounding basement membrane [29, 32]. Researchers reported that EMT coordinately regulated maintenance of cancer stemness, drug resistance, angiogenesis, and muscle invasion/metastasis, and had been proofed to be a major characteristic [33–35]. Successful elimination of subpopulation of UroCSCs and its differentiated progenies through targeted treatment of EMT and its core partners can clinically reverse EMT, prevent tumor recurrence, and prolong patient survival [36, 37].

Additionally, these genes uncovered by KEGG pathway analysis are remarkably related to PI3K-Akt signaling pathway, focal adhesion, proteoglycans in cancer, ECM-receptor interaction, and oxidative phosphorylation. Regarding bladder cancer, constitutive activation of PI3K/AKT pathway was observed in up to 40% of tumors [38], which pathway combined a serial of external signals to regulate downstream signaling involved in the cell growth, differentiation, angiogenesis, and protein synthesis, ultimately resulting in bladder carcinogenesis, metastatic potential, and therapy resistance [39–43]. The extracellular matrix (ECM) comprises noncellular components playing an important part in the cell behaviors of regulation [44], which dynamically changes in tumor microenvironment and has a crucial role in tumor progression [45]. Tumor cell growth, migration, and apoptosis via regulation of receptor interaction were reported to be related to the extracellular matrix protein 1 (ECM1) expression, which was considered as a predictive parameter in the carcinogenesis and postoperative recurrence of BLCA [46].

One of the factors affecting the occurrence and progression of tumor is tumor mutation burden (TMB) [47]. Some genes involved in TMB have been reported to be of predictive accuracy in OS among BLCA patients [47–49]. In the present work, PTCD1 expression level was highly impacted by the mutation rate of same cancer-related genes, which were of high mutation probability in BLCA. The relationship between mutant genes and PTCD1 expression implies that PTCD1 may play a regulatory part in BLCA, which is needed to be validated by further studies.

The tumor immune cell infiltration is of essential importance for the prognosis of bladder cancer and interferes the response to immunotherapy [50–52]. Several studies have demonstrated the importance of tumor-infiltrating immune cells and other immune molecules (including myeloid dendritic cell, T cell, neutrophil, and macrophage) in the prognosis of bladder cancer [53–56]. Accordingly, infiltration rates of immune cells with the low- and high-PTCD1 expression levels were further analyzed, suggesting that PTCD1 expression was negatively associated with immune cell infiltration, especially for myeloid dendritic cell, T cell, neutrophil, and macrophage, which was consistent with previous work. The level of immune checkpoint molecules was determined to obtain a deeper understanding into the immune landscape of BLCA. Heterogeneous expression of immune checkpoint proteins in the immune microenvironment of BLCA was discovered by the correlation between its expression and PTCD1 expression. Many antigens presentation-related functions were found by the enrichment score analysis of immune-related pathways. For instance, in the low-PTCD1 group, functions including coinhibition and costimulation of APC, CCR, checkpoint, cytolytic activity, HLA, inflammation promoting, MHC class I, para-inflammation, coinhibition, costimulation of T cell, and Type I IFN response were observed. Taken together, we assume that PTCD1 level may indirectly affect the occurrence and development of BLCA via immune cell.

RNA N6-methyladenosine (m6A) modification as the most abundant mRNA modification has attracted much attention recently [57]. The modification of m6A plays a role in the regulation of occurrence and development of tumors by modulating expressions of key genes, becoming the hot spot of research recently [58]. It was reported by previous research that the development of BLCA may be affected by the abnormal modification of m6A methylation [59–62]. The expression levels of 16 regulators related to m6A found in this work greatly depend on the PTCD1 expression level of BLCA. Reasonably, downregulation of PTCD1 in bladder cancer likely causes the abnormal regulation of m6A-related regulators, where more studies need to be done to prove the hypothesis.

In conclusion, PTCD1 was downregulated in BLCA, which was greatly associated with the poor clinical features and prognosis. Further, our data suggested that PTCD1 may participate in the immunological functioning and immune cell infiltration and might play certain roles in tumor progression and metastasis, which could be acted as a potential prognostic target of BLCA. Likewise, there are few limitations existing in this work. First, more long-term survival data of clinical samples are required to support our findings because of the retrospective data and the limitations of clinical information from TCGA. Second, the absence of PTCD1 protein expression level in tumor specimens should be further detected. Thirdly, potential mechanisms underlying the effects of PTCD1 on the clinical outcomes in patients with BLCA were not fully explored. Last but not the least, more studies are needed to investigate the interplay between PTCD1 expression and immune infiltration.

5. Conclusions

In this work, the role of PTCD1 in BLCA was explored from the perspectives of clinical values and potential mechanisms of PTCD1. Downregulation of PTCD1 is likely related to BLCA progression. PTCD1 may play some certain parts in cell invasion and metabolic and immune microenvironment of BLCA. All in all, PTCD1 could act as a potential predictor for diagnosis and prognosis and novel therapeutic target in BLCA.

Data Availability

Publicly available datasets were analyzed in this study. HPA database and TCGA database belong to public databases. Users can download relevant data for free for research and publish relevant articles.

Conflicts of Interest

The authors have no conflicts of interest to declare regarding the publication of this paper.

Authors’ Contributions

All authors constructed this study. ZZB, CYS, ZYJ, ZX, and HYJ performed the data analysis, figures plotting, and writing. ZY was responsible for the critical reading of the manuscript. All authors contributed to the article and approved the submitted version. Yong Zhang was the first corresponding author of this article and Yuanshan Cui was the co corresponding author of this article. Yong Zhang was the first corresponding author of this article and Yuanshan Cui was the co corresponding author of this article.

Acknowledgments

The authors would like to thank 51runse for the English language editing service. This work was supported by Beijing Municipal Administration of Hospitals’ Ascent Plan (Code: DFL20190502); Beijing Municipal Administration of Hospitals Clinical Medicine Development of Special Funding Support (Code: ZYLX201820); and Joint Fund of Shandong Natural Science Foundation(Code: ZR2021LSW019).

Supplementary Materials

Supplementary File 1. The differential expression genes between low-PTCD1 group and high-PTCD1 group. Supplementary File 2. The results of GO enrichment analysis and KEGG pathway analysis. (Supplementary Materials)

References

S. Antoni, J. Ferlay, I. Soerjomataram, A. Znaor, A. Jemal, and F. Bray, “Bladder cancer incidence and mortality: a global overview and recent trends,” European Urology, vol. 71, no. 1, pp. 96–108, 2017.
View at: Publisher Site | Google Scholar
A. Lopez-Beltran, A. Cimadamore, A. Blanca et al., “Immune checkpoint inhibitors for the treatment of bladder cancer,” Cancers, vol. 13, no. 1, p. 131, 2021.
View at: Publisher Site | Google Scholar
A. T. Lenis, P. M. Lec, K. Chamie, and M. D. Mshs, “Bladder cancer: a review,” JAMA, vol. 324, no. 19, p. 1980, 2020.
View at: Publisher Site | Google Scholar
H. Kaseb and N. R. Aeddula, Bladder Cancer. StatPearls. Treasure Island (FL), StatPearls Publishing LLC, St. Petersburg, FL, USA, 2022.
S. K. Bhanvadia, “Bladder cancer survivorship,” Current Urology Reports, vol. 19, no. 12, p. 111, 2018.
View at: Publisher Site | Google Scholar
M. Babjuk, M. Burger, O. Capoun et al., “European association of Urology guidelines on non-muscle-invasive bladder cancer (Ta, T1, and carcinoma in situ),” European Urology, vol. 81, no. 1, pp. 75–94, 2022.
View at: Publisher Site | Google Scholar
S. van den Bosch and J. Alfred Witjes, “Long-term cancer-specific survival in patients with high-risk, non-muscle-invasive bladder cancer and tumour progression: a systematic review,” European Urology, vol. 60, no. 3, pp. 493–500, 2011.
View at: Publisher Site | Google Scholar
D. J. McConkey and W. Choi, “Molecular subtypes of bladder cancer,” Current Oncology Reports, vol. 20, no. 10, p. 77, 2018.
View at: Publisher Site | Google Scholar
S. S. Sugeeta, A. Sharma, K. Ng, A. Nayak, and N. Vasdev, “Biomarkers in bladder cancer surveillance,” Front Surg, vol. 8, Article ID 735868, 2021.
View at: Publisher Site | Google Scholar
C. Zhang, J. Hu, H. Li et al., “Emerging biomarkers for predicting bladder cancer lymph node metastasis,” Frontiers Oncology, vol. 11, Article ID 648968, 2021.
View at: Publisher Site | Google Scholar
J. Dobruch and M. Oszczudłowski, “Bladder cancer: current challenges and future directions,” Medicina, vol. 57, no. 8, p. 749, 2021.
View at: Publisher Site | Google Scholar
I. Faiena, C. J. Rosser, K. Chamie, and H. Furuya, “Diagnostic biomarkers in non-muscle invasive bladder cancer,” World Journal of Urology, vol. 37, no. 10, pp. 2009–2016, 2019.
View at: Publisher Site | Google Scholar
R. N. Lightowlers and Z. M. Chrzanowska-Lightowlers, “Human pentatricopeptide proteins: only a few and what do they do?” RNA Biology, vol. 10, no. 9, pp. 1433–1438, 2013.
View at: Publisher Site | Google Scholar
K. L. Perks, N. Ferreira, T. R. Richman et al., “Adult-onset obesity is triggered by impaired mitochondrial gene expression,” Science Advances, vol. 3, no. 8, Article ID e1700677, 2017.
View at: Publisher Site | Google Scholar
W. S. Wei, N. Wang, M. H. Deng et al., “LRPPRC regulates redox homeostasis via the circANKHD1/FOXM1 axis to enhance bladder urothelial carcinoma tumorigenesis,” Redox Biology, vol. 48, Article ID 102201, 2021.
View at: Publisher Site | Google Scholar
I. Kühl, M. Miranda, V. Posse et al., “POLRMT regulates the switch between replication primer formation and gene expression of mammalian mtDNA,” Science Advances, vol. 2, no. 8, Article ID e1600963, 2016.
View at: Publisher Site | Google Scholar
O. Rackham, S. M. K. Davies, A. M. J. Shearwood, K. L. Hamilton, J. Whelan, and A. Filipovska, “Pentatricopeptide repeat domain protein 1 lowers the levels of mitochondrial leucine tRNAs in cells,” Nucleic Acids Research, vol. 37, no. 17, pp. 5859–5867, 2009.
View at: Publisher Site | Google Scholar
J. Holzmann, P. Frank, E. Löffler, K. L. Bennett, C. Gerner, and W. Rossmanith, “RNase P without RNA: identification and functional reconstitution of the human mitochondrial tRNA processing enzyme,” Cell, vol. 135, no. 3, pp. 462–474, 2008.
View at: Publisher Site | Google Scholar
F. Sasarman, C. Brunel-Guitton, H. Antonicka, T. Wai, and E. A. Shoubridge, “LRPPRC and SLIRP interact in a ribonucleoprotein complex that regulates posttranscriptional gene expression in mitochondria,” Molecular Biology of the Cell, vol. 21, no. 8, pp. 1315–1323, 2010.
View at: Publisher Site | Google Scholar
S. M. Davies, O. Rackham, A. M. J. Shearwood et al., “Pentatricopeptide repeat domain protein 3 associates with the mitochondrial small ribosomal subunit and regulates translation,” FEBS Letters, vol. 583, no. 12, pp. 1853–1858, 2009.
View at: Publisher Site | Google Scholar
M. Black, P. Arumugam, S. Shukla et al., “FOXM1 nuclear transcription factor translocates into mitochondria and inhibits oxidative phosphorylation,” Molecular Biology of the Cell, vol. 31, no. 13, pp. 1411–1424, 2020.
View at: Publisher Site | Google Scholar
Z. Sun, Y. Mao, X. Zhang et al., “Identification of ARHGEF38, NETO2, GOLM1, and SAPCD2 associated with prostate cancer progression by bioinformatic analysis and experimental validation,” Frontiers in Cell and Developmental Biology, vol. 9, Article ID 718638, 2021.
View at: Publisher Site | Google Scholar
C. Schild, D. Hahn, A. Schaller et al., “Mitochondrial leucine tRNA level and PTCD1 are regulated in response to leucine starvation,” Amino Acids, vol. 46, no. 7, pp. 1775–1783, 2014.
View at: Publisher Site | Google Scholar
D. Fleck, L. Phu, E. Verschueren et al., “PTCD1 is required for mitochondrial oxidative-phosphorylation: possible genetic association with alzheimer’s disease,” Journal of Neuroscience, vol. 39, no. 24, pp. 4636–4656, 2019.
View at: Publisher Site | Google Scholar
K. L. Perks, N. Ferreira, J. A. Ermer et al., “Reduced mitochondrial translation prevents diet-induced metabolic dysfunction but not inflammation,” Aging (Albany NY), vol. 12, no. 19, pp. 19677–19700, 2020.
View at: Publisher Site | Google Scholar
J. Pa, S. J. Andrews, and R. H. Swerdlow, “Mitochondria and alzheimer’s: is PTCD1 the smoking gun?” Trends in Neurosciences, vol. 42, no. 11, pp. 759–762, 2019.
View at: Publisher Site | Google Scholar
M. Sánchez-Carbayo, “Hypermethylation in bladder cancer: biological pathways and translational applications,” Tumor Biology, vol. 33, no. 2, pp. 347–361, 2012.
View at: Publisher Site | Google Scholar
T. Reinert, C. Modin, F. M. Castano et al., “Comprehensive genome methylation analysis in bladder cancer: identification and validation of novel methylated genes and application of these as urinary tumor markers,” Clinical Cancer Research, vol. 17, no. 17, pp. 5582–5592, 2011.
View at: Publisher Site | Google Scholar
M. Garg and R. Singh, “Epithelial-to-mesenchymal transition: event and core associates in bladder cancer,” Frontiers in Bioscience, vol. 11, no. 1, pp. 853–865, 2019.
View at: Publisher Site | Google Scholar
N. Dillon, “Gene regulation and large-scale chromatin organization in the nucleus,” Chromosome Research, vol. 14, no. 1, pp. 117–126, 2006.
View at: Publisher Site | Google Scholar
F. Yang, L. Sun, Q. Li et al., “SET8 promotes epithelial-mesenchymal transition and confers TWIST dual transcriptional activities,” The EMBO Journal, vol. 31, no. 1, pp. 110–123, 2012.
View at: Publisher Site | Google Scholar
T. Müller, G. Bain, X. Wang, and J. Papkoff, “Regulation of epithelial cell migration and tumor formation by beta-catenin signaling,” Experimental Cell Research, vol. 280, no. 1, pp. 119–133, 2002.
View at: Publisher Site | Google Scholar
N. V. Jordan, G. L. Johnson, and A. N. Abell, “Tracking the intermediate stages of epithelial-mesenchymal transition in epithelial stem cells and cancer,” Cell Cycle, vol. 10, no. 17, pp. 2865–2873, 2011.
View at: Publisher Site | Google Scholar
E. M. Peek, D. R. Li, H. Zhang et al., “Stromal modulation of bladder cancer-initiating cells in a subcutaneous tumor model,” American Journal of Cancer Research, vol. 2, no. 6, pp. 745–751, 2012.
View at: Google Scholar
C. A. Franzen, R. H. Blackwell, V. Todorovic et al., “Urothelial cells undergo epithelial-to-mesenchymal transition after exposure to muscle invasive bladder cancer exosomes,” Oncogenesis, vol. 4, no. 8, 2015.
View at: Publisher Site | Google Scholar
M. A. Skowron, G. Niegisch, G. Fritz et al., “Phenotype plasticity rather than repopulation from CD90/CK14⁺ cancer stem cells leads to cisplatin resistance of urothelial carcinoma cell lines,” Journal of Experimental and Clinical Cancer Research, vol. 34, no. 1, p. 144, 2015.
View at: Publisher Site | Google Scholar
Y. Wang, G. Wang, X. Zhang et al., “γ-Secretase inhibitor inhibits bladder cancer cell drug resistance and invasion by reducing epithelial-mesenchymal transition,” Molecular Medicine Reports, vol. 12, no. 2, pp. 2821–2827, 2015.
View at: Publisher Site | Google Scholar
S. T. Liu, G. Hui, C. Mathis, K. Chamie, A. J. Pantuck, and A. Drakaki, “The current status and future role of the phosphoinositide 3 kinase/AKT signaling pathway in urothelial cancer: an old pathway in the new immunotherapy era,” Clinical Genitourinary Cancer, vol. 16, no. 2, pp. e269–e276, 2018.
View at: Publisher Site | Google Scholar
X. Gao, J. Liu, D. Fan et al., “Berberine enhances gemcitabine-induced cytotoxicity in bladder cancer by downregulating Rad51 expression through inactivating the PI3K/Akt pathway,” Oncology Reports, vol. 47, no. 2, 2021.
View at: Publisher Site | Google Scholar
G. Yang, Z. Li, L. Dong, and F. Zhou, “lncRNA ADAMTS9-AS1 promotes bladder cancer cell invasion, migration, and inhibits apoptosis and autophagy through PI3K/AKT/mTOR signaling pathway,” The International Journal of Biochemistry and Cell Biology, vol. 140, Article ID 106069, 2021.
View at: Publisher Site | Google Scholar
Y. Xue, Y. Zhong, T. Wu et al., “Anti-proliferative and apoptosis-promoting effect of microRNA-125b on pancreatic cancer by targeting NEDD9 via PI3K/AKT signaling,” Cancer Management and Research, vol. 12, pp. 7363–7373, 2020.
View at: Publisher Site | Google Scholar
Y. L. Su, H. L. Luo, C. C. Huang et al., “Galectin-1 overexpression activates the FAK/PI3K/AKT/mTOR pathway and is correlated with upper urinary urothelial carcinoma progression and survival,” Cells, vol. 9, no. 4, p. 806, 2020.
View at: Publisher Site | Google Scholar
A. Sathe and R. Nawroth, “Targeting the PI3K/AKT/mTOR pathway in bladder cancer,” Methods in Molecular Biology, vol. 1655, pp. 335–350, 2018.
View at: Publisher Site | Google Scholar
S. Appunni, V. Anand, M. Khandelwal, N. Gupta, M. Rubens, and A. Sharma, “Small Leucine Rich Proteoglycans (decorin, biglycan and lumican) in cancer,” Clinica Chimica Acta, vol. 491, pp. 1–7, 2019.
View at: Publisher Site | Google Scholar
P. Lu, V. M. Weaver, and Z. Werb, “The extracellular matrix: a dynamic niche in cancer progression,” Journal of Cell Biology, vol. 196, no. 4, pp. 395–406, 2012.
View at: Publisher Site | Google Scholar
Z. Wang, Q. Zhou, A. Li, W. Huang, Z. Cai, and W. Chen, “Extracellular matrix protein 1 (ECM1) is associated with carcinogenesis potential of human bladder cancer,” OncoTargets and Therapy, vol. 12, pp. 1423–1432, 2019.
View at: Publisher Site | Google Scholar
M. G. McNamara, T. Jacobs, A. Lamarca, R. A. Hubner, J. W. Valle, and E. Amir, “Impact of high tumor mutational burden in solid tumors and challenges for biomarker application,” Cancer Treatment Reviews, vol. 89, Article ID 102084, 2020.
View at: Publisher Site | Google Scholar
G. Zhu, L. Pei, Y. Li, and X. Gou, “EP300 mutation is associated with tumor mutation burden and promotes antitumor immunity in bladder cancer patients,” Aging (Albany NY), vol. 12, no. 3, pp. 2132–2141, 2020.
View at: Publisher Site | Google Scholar
D. Liu, P. Abbosh, D. Keliher et al., “Mutational patterns in chemotherapy resistant muscle-invasive bladder cancer,” Nature Communications, vol. 8, no. 1, p. 2193, 2017.
View at: Publisher Site | Google Scholar
Z. Yang, Y. Xu, Y. Bi et al., “Immune escape mechanisms and immunotherapy of urothelial bladder cancer,” Journal of Clinical Translational Research, vol. 7, no. 4, pp. 485–500, 2021.
View at: Google Scholar
H. Gil-Julio, F. Perea, A. Rodriguez-Nicolas et al., “Tumor escape phenotype in bladder cancer is associated with loss of HLA class I expression, T-cell exclusion and stromal changes,” International Journal of Molecular Sciences, vol. 22, no. 14, 2021.
View at: Publisher Site | Google Scholar
Q. Zhong, J. Shou, J. Ying et al., “High PD-L1 expression on immune cells, but not on tumor cells, is a favorable prognostic factor in urothelial carcinoma,” Future Oncology, vol. 17, no. 22, pp. 2893–2905, 2021.
View at: Publisher Site | Google Scholar
R. Rossi, M. Lichtner, F. Iori et al., “Dendritic cells in blood and urine samples from bladder cancer patients undergoing BCG immunotherapy,” Italian Archives of Urology, Andrology, vol. 85, no. 4, 2013.
View at: Publisher Site | Google Scholar
T. Anami, Y. Komohara, Y. Miura et al., “High T-cell infiltration in tumor tissue and younger age predict the response to pembrolizumab in recurrent urothelial cancer,” Medical Molecular Morphology, vol. 54, no. 4, pp. 316–323, 2021.
View at: Publisher Site | Google Scholar
S. Shinde-Jadhav, J. J. Mansure, R. F. Rayes et al., “Role of neutrophil extracellular traps in radiation resistance of invasive bladder cancer,” Nature Communications, vol. 12, no. 1, 2021.
View at: Publisher Site | Google Scholar
S. Jin, H. Zeng, Z. Liu et al., “Stromal tumor-associated macrophage infiltration predicts poor clinical outcomes in muscle-invasive bladder cancer patients,” Annals of Surgical Oncology, vol. 29, no. 4, pp. 2495–2503, 2022.
View at: Publisher Site | Google Scholar
Y. An and H. Duan, “The role of m6A RNA methylation in cancer metabolism,” Molecular Cancer, vol. 21, no. 1, p. 14, 2022.
View at: Publisher Site | Google Scholar
N. Zhang, Y. Zuo, Y. Peng, and L. Zuo, “Function of N6-methyladenosine modification in tumors,” Journal of Oncology, vol. 2021, Article ID 6461552, 10 pages, 2021.
View at: Publisher Site | Google Scholar
M. Chen, Z. Y. Nie, X. H. Wen, Y. H. Gao, H. Cao, and S. F. Zhang, “m6A RNA methylation regulators can contribute to malignant progression and impact the prognosis of bladder cancer,” Bioscience Reports, vol. 39, no. 12, Article ID BSR20192892, 2019.
View at: Publisher Site | Google Scholar
M. Cheng, L. Sheng, Q. Gao et al., “The m (6)A methyltransferase METTL3 promotes bladder cancer progression via AFF4/NF-κB/MYC signaling network,” Oncogene, vol. 38, no. 19, pp. 3667–3680, 2019.
View at: Publisher Site | Google Scholar
H. Jin, X. Ying, B. Que et al., “N (6)-methyladenosine modification of ITGA6 mRNA promotes the development and progression of bladder cancer,” EBioMedicine, vol. 47, pp. 195–207, 2019.
View at: Publisher Site | Google Scholar
H. Xie, J. Li, Y. Ying et al., “METTL3/YTHDF2 m (6) A axis promotes tumorigenesis by degrading SETD7 and KLF4 mRNAs in bladder cancer,” Journal of Cellular and Molecular Medicine, vol. 24, no. 7, pp. 4092–4104, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Zhongbao Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

250

Downloads

467

Citations