frontiers in Pharmacology
Check for updates
New Insights Into PTBP3 in Human Cancers: Immune Cell Infiltration, TMB, MSI, PDCD1 and m6A Markers
Zhen Fang11, Peijuan Li21, Han Li3, Wei Chong4, Leping Li4, Liang Shang4* and Fei Li1*
1Department of General Surgery, Xuanwu Hospital, Capital Medical University, Beijing, China, 2Emergency Department, The First Affiliated Hospital of Dalian Medical University, Dalian, China, 3Department of General Surgery, Shandong Provincial Qianfoshan Hospital, The First Affiliated Hospital of Shandong First Medical University, Jinan, China, 4Department of Gastrointestinal Surgery, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, China
OPEN ACCESS
Edited by: Tzi Bun Ng, The Chinese University of Hong Kong, China
Reviewed by:
Mohd Ahmar Rauf, University of Miami Hospital, United States Angelique Nyinawabera, L.E.A.F. Pharmaceuticals, United States
*Correspondence:
Liang Shang docshang@163.com
Fei Li lfei21@126.com
+These authors have contributed equally to this work
Specialty section:
This article was submitted to Pharmacology of Anti-Cancer Drugs, a section of the journal Frontiers in Pharmacology Received: 08 November 2021 Accepted: 09 February 2022 Published: 10 March 2022
Citation:
Fang Z, Li P, Li H, Chong W, Li L, Shang L and Li F (2022) New Insights Into PTBP3 in Human Cancers: Immune Cell Infiltration, TMB, MSI, PDCD1 and m6A Markers. Front. Pharmacol. 13:811338. doi: 10.3389/fphar.2022.811338
Polypyrimidine tract binding protein 3 (PTBP3) plays a critical role in post-transcriptional regulation. The role of PTBP3 in various human tumours was explored and analysed in this study based on the Cancer Genome Atlas and Gene Expression Omnibus datasets. PTBP3 was highly expressed in most tumours, such as breast invasive carcinoma, colon adenocarcinoma and hepatocellular carcinoma. PTBP3 overexpression generally predicts poor overall survival and disease-free survival in patients with adrenocortical carcinoma, lung squamous cell carcinoma, and pancreatic adenocarcinoma. However, low PTBP3 expression predicts poor prognosis in kidney renal clear cell carcinoma. We also explored PTBP3 genetic alterations in different tumour tissues. The result found that the frequency of PTBP3 alteration (>4%) was the highest in uterine tumours with “mutation” as the primary type. Furthermore, we found a significant correlation between PTBP3 expression and tumour mutational burden and microsatellite instability in various human tumours, and found that PTBP3 expression was positively correlated with TMB in ACC, STAD, PAAD, LUAD, and SARC. Two enhanced phosphorylation levels of S30 and S426 in colon cancer, ovarian cancer, and uterine corpus endometrial carcinoma were found. Further analysis indicated that PTBP3 expression was positively correlated with the cancer-associated fibroblasts for most tumour types. This study also found a relationship between immune checkpoints and N6-methyladenosine-related markers and PTBP3 expression. Moreover, the “mRNA surveillance pathway” and “RNA degradation” were involved in the functional mechanisms of PTBP3. These results provide new insights for molecular studies, and integrative analysis provided a framework for determining the predictive, prognostic, and therapeutic relevance of PTBP3 in cancer patients.
Abbreviations: TCGA, The Cancer Genome Atlas; GTEx, The Genotype-Tissue Expression; TMB, Tumor Mutation Burden; MSI, Microsatellite Instability; ACC, Adrenocortical carcinoma; BLCA, Bladder Urothelial Carcinoma; BRCA, Breast invasive carcinoma; CESC, Cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOL, Cholangiocarcinoma; COAD, Colon adenocarcinoma; DLBC, Lymphoid Neoplasm Diffuse Large B-cell Lymphoma; ESCA, Esophageal carcinoma; GBM, Glioblastoma multiforme; HNSC, Head and Neck squamous cell carcinoma; KICH, Kidney Chromophobe; KIRC, Kidney renal clear cell carcinoma; KIRP, Kidney renal papillary cell carcinoma; LGG, Brain Lower Grade Glioma; LIHC, Liver hepatocellular carcinoma; LUAD, Lung adenocarcinoma; LUSC, Lung squamous cell carcinoma; MESO, Mesothelioma; OV, Ovarian serous cystadenocarcinoma; PAAD, Pancreatic adenocarcinoma; PCPG, Pheochromocytoma and Paraganglioma; PRAD, Prostate adenocarcinoma; READ, Rectum adenocarcinoma; SARC, Sarcoma; SKCM, Skin Cutaneous Melanoma; STAD, Stomach adenocarcinoma; TGCT, Testicular Germ Cell Tumours; THCA, Thyroid carcinoma; THYM, Thymoma; UCEC, Uterine Corpus Endometrial Carcinoma; UCS, Uterine Carcinosarcoma; UVM, Uveal Melanoma.
Keywords: PTBP3, pan-cancer, prognosis, immunotherapy, phosphorylation, m6A
INTRODUCTION
Polypyrimidine tract binding protein 3 (PTBP3), also known as the regulator of differentiation (ROD) 1, is located on human chromosome 9q32. PTBP3 is a member of the polypyrimidine tract binding (PTB) protein family, which also includes PTBP1 and PTBP2 (Hou et al., 2019). The PTB family is known to bind to CU-rich elements in exon and intron regions (Marinescu et al., 2007). The PTB family binding of CU-rich sequences can also influence splice site selection by preventing the transition from exons to introns, or by obstructing the definition of exons and introns (Sharma et al., 2008; Ma et al., 2020).
Previous studies have indicated that PTBP3 is crucial in the complex process of tumorigenesis (Hou et al., 2018; Liang et al., 2018; Yang et al., 2018; Hou et al., 2019; Ma et al., 2020; Wu et al., 2020; Chen et al., 2022). For example, PTBP3 regulates lung cancer cell proliferation through the cell cycle and may be a potential target for lung cancer molecular therapy (Chen et al., 2022). PTBP3 can promote colorectal cell proliferation, migration, and invasion in vitro and tumour growth and metastasis in vivo by binding to the 5’UTR HIF-1a mRNA to enhance HIF-la protein expression (Hou et al., 2019). PTBP3 overexpression also promotes pancreatic ductal adenocarcinoma proliferation in vitro and tumour growth in vivo (Ma et al., 2020). In addition, PTBP3 can regulate the expression of ZEB1 (epithelial-mesenchymal transition regulatory transcription factor) and promotes tumour cell invasive growth and metastasis in breast cancer (Hou et al., 2018). PTBP3 expression is positively correlated with lymph node metastasis and poor 5-year survival in patients with breast cancer (Hou et al., 2018). Therefore, these current findings imply that PTBP3 may be a new and promising tumour biomarker.
In this study, The Cancer Genome Atlas (TCGA) Project and Gene Expression Omnibus (GEO) databases were used to explore the expression profile of PTBP3 across various tumour types. The main investigated contents are as follows: gene expression profile, survival status, genetic alteration, protein phosphorylation, immune infiltration, N6-methyladenosine, and relevant cellular pathways. Comprehensive pan-cancer analysis aims to investigate the potential molecular mechanism of PTBP3 in the pathogenesis, clinical prognosis, and immunotherapy of various human tumours.
MATERIALS AND METHODS
Gene Expression Analysis
We used Tumour Immune Estimation Resource, Version 2 (TIMER2, http://timer.cistrome.org/) to analyse the differences in PTBP3 expression in different tumour types and adjacent normal tissues. However, lymphoid neoplasm diffuse large B-cell (DLBC) lymphoma, brain lower grade glioma (LGG), ovarian
serous cystadenocarcinoma (OV), testicular germ cell tumours (TGCT), and thymoma (THYM) had only tumour tissue but no adjacent normal tissue. The Gene Expression Profiling Interactive Analysis, version 2 (GEPIA2) website was used to analyse PTBP3 expression in DLBC, LGG, OV, TGCT, and THYM (Tang et al., 2019). UALCAN can provide protein expression analysis options using data from the Clinical Proteomic Tumour Analysis Consortium (CPTAC) dataset. We used the UALCAN website (http://ualcan.path.uab.edu/index.html) to analyse the protein expression of PTBP3 in breast cancer, colon cancer, uterine corpus endometrial carcinoma (UCEC), lung adenocarcinoma, and ovarian cancer. The GEPIA2 website was also used to analyse PTBP3 expression in different pathological stages of all TCGA cancers.
Survival Prognosis Analysis
The GEPIA2 tool was used to obtain the overall survival (OS) and disease-free survival (DFS) significance map data across all TCGA tumours (Tang et al., 2019). The median was used as the expression threshold for splitting the high- and low- expression groups. Month units were selected for the plotting. The Cox proportional hazard ratio and 95% confidence interval information were used in the survival plot. The log-rank test was used for the hypothesis testing.
Genetic Alteration Analysis
The cBioPortal tool was used to analyse the alteration frequency, mutation type, copy number alteration, and mutated site information across all TCGA tumours (Cerami et al., 2012; Gao et al., 2013). After logging in to the cBioPortal website, select “TCGA Pan Cancer Atlas Studies” to query the genetic alteration characteristics of PTBP3. We could observe results for alteration frequency, mutation type and copy number alterations (CNAs) in all TCGA tumours. The mutation site information of PTBP3 can be displayed through the “Mutations” module. At the same time, the 3D structure of the protein can also be obtained. The OS and DFS data for all TCGA cancer types were compared with or without PTBP3 genetic alterations.
Immune Infiltration and N6-Methyladenosine-Related Makers Analysis
We analysed the association between PTBP3 expression and immune infiltrates (e.g., cancer associated fibroblast and neutrophil) across all TCGA tumours and the relationship between PTBP3 expression and N6-methyladenosine (m6A)- related markers (e.g., METTL3, KIAA1429, FTO, ZC3H13, YTHDF1/2/3, RBM15) using the TIMER2 tool. The relationship between PTBP3 expression and immune checkpoints (e.g., PD1, PDL1, CTLA4, CD274, LAG3, SIGLEC15, TIGIT) was explored using the HOME FOR RESEARCHERS website tool (https://www.aclbi.com/).
A
TCGA dataset
PTBP3 Expression (log2 TPM)
8
**
**
*
*
6
LI
A
0
-
2
ACC.Tumor (n=79)
BLCA.Tumor (n=408)
BLCA.Normal (n=19)
BRCA. Tumor (n=1093)
BRCA.Normal (n=112)
BRCA-Basal. Tumor (n=190)
BRCA-Her2. Tumor (n=82)
BRCA-LumA. Tumor (n=564)
BRCA-LumB. Tumor (n=217)
CESC.Tumor (n=304)
CESC.Normal (n=3)
CHOL. Tumor (n=36)
CHOL.Normal (n=9)
COAD. Tumor (n=457)
COAD.Normal (n=41)
DLBC. Tumor (n=48)
ESCA. Tumor (n=184)
ESCA.Normal (n=11)
GBM. Tumor (n=153)
GBM.Normal (n=5)
HNSC.Tumor (n=520)
HNSC.Normal (n=44)
HNSC-HPV+. Tumor (n=97)
HNSC-HPV -. Tumor (n=421)
KICH. Tumor (n=66)
KICH.Normal (n=25)
KIRC. Tumor (n=533)
KIRC. Normal (n=72)
KIRP.Tumor (n=290)
KIRP.Normal (n=32)
LAML. Tumor (n=173)
LGG. Tumor (n=516)
LIHC. Tumor (n=371)
LIHC.Normal (n=50)
LUAD. Tumor (n=515)
LUAD.Normal (n=59)
LUSC. Tumor (n=501)
LUSC.Normal (n=51)
MESO. Tumor (n=87)
OV.Tumor (n=303)
PAAD. Tumor (n=178)
PAAD.Normal (n=4)
PCPG. Tumor (n=179)
PCPG.Normal (n=3)
PRAD.Tumor (n=497)
PRAD.Normal (n=52)
READ.Tumor (n=166)
READ.Normal (n=10)
SARC. Tumor (n=259)
SKCM.Tumor (n=103)
SKCM.Metastasis (n=368)
STAD. Tumor (n=415)
STAD.Normal (n=35)
TGCT.Tumor (n=150)
THCA. Tumor (n=501)
THCA.Normal (n=59)
THYM. Tumor (n=120)
UCEC. Tumor (n=545)
UCEC.Normal (n=35)
UCS. Tumor (n=57)
UVM. Tumor (n=80)
B TCGA + GTEx dataset
Expression -log2(TPM+1)
6
P
+
4
..:
1.
2
2
M
1
Y
C
.FR
N
9
1
0
DLBC
LGG
OV
PAAD
TGCT
THYM
(T=47; N=337) (T=518; N=207) (T=426; N=88) (T=179; N=171) (T=137; N=165) (T=118; N=339)
C CPTAC dataset
ACUNHO-NO
Breast cancer
Z-value
bb ANONA
Colon cancer
3 WON TONW
UCEC
3.
Lung adenocarcinoma
3
Ovarian cancer
Z-value
Z-value
WNO-NO
Z-value
Z-value
1.
2
Normal Primary tumor (n=18) (n=125)
Normal (n=100)
Primary tumor (n=97)
3
Normal Primary tumor (n=31) (n=100)
-4
Normal (n=111)
Primary tumor (n=111)
3
Normal Primary tumor (n=25) (n=100)
D
TCGA dataset
F value = 4.61
F value = 5.93
Pr(>F) = 0.000554
F value = 3.64
F value = 2.43
5
ACC
Pr(>F) = 0.00526
7
KIRC
6
LIHC
Pr(>F) = 0.0131
₪
SKCM
Pr(>F) = 0.0475
4
6
5
5.
5
3
+
+
+
2
0
3
0
2
2-
-
2
-
-
li
Ili
IV
li
III
IV
I
li
Iİl
IV
0
I
II
lil
IV
Stage
Stage
Stage
Stage
PTBP3-Related Gene Enrichment Analysis
The STRING tool (https://string-db.org/) was used to construct the protein-protein interaction network (Szklarczyk et al., 2021). The main parameters are as follow: minimum required interaction score [“Low confidence (0.150)”], max number of interactors to show (“no more than 50 interactors” in the 1st shell), meaning of network edges (“evidence”), and active interaction sources (“experiments”). The GEPIA2 tool was used to obtain the top 100 PTBP3-related targeting genes based on the database of all TCGA tissues. Furthermore, we performed pairwise gene correlation analysis of PTBP3 and the selected genes. Pearson’s correlation analysis was also used to calculate the correlation coefficient, while the heatmap of the selected genes contained the partial correlation and p-value in the purity-adjusted Spearman’s rank correlation. The R package of “tidyr” and “ggplot2” were used to visualise the enrichment pathways. R software (R-4.0.2, 64-bit) was used in this study (https://www.r-project.org/), and statistical significance was set at p < 0.05 (Cui et al., 2020; Huo et al., 2021).
RESULTS
Gene Expression Analysis Data
In this study, we have provided a comprehensive analysis of the role of PTBP3 (NM_001163788.4 for NP_001157260.1, Supplementary Figure S1) in humans. We analysed the expression of PTBP3 in different cell lines and normal tissues (Supplementary Figure S2). We found that PTBP3 expression in BeWo cells (human placental choriocarcinoma cells) was the highest in the Human Protein Atlas (HPA) datasets, followed by K-562 (chronic myelogenous leukaemia cells) and OE19 (human oesophageal cancer cells) (Supplementary Figure S2A). As shown in Supplementary Figure S2B, we collected data from the Consensus, HPA, Genotype-Tissue Expression (GTEx), and FANTOM5 (Function Annotation of the Mammalian Genome 5) datasets. This comparison showed high expression of PTBP3 in the oesophagus, tonsils, lymph nodes, and thymus. However, PTBP3 was expressed in all tissues (all consensus normalised expression values >1), and it showed low tissue specificity.
Next, we used TIMER2 to explore the expression of PTBP3 in different tumour types in the TCGA repository. As shown in Figure 1A, the expression levels of PTBP3 in breast invasive carcinoma (BRCA), cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), oesophageal carcinoma (ESCA), glioblastoma multiforme (GBM), head and neck squamous cell carcinoma (HNSC), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), rectum adenocarcinoma (READ), stomach adenocarcinoma (STAD), and uterine corpus endometrial carcinoma (UCEC) (all p < 0.05) were higher than those in the corresponding adjacent normal tissues. Only a few tumour types showed no differential expression (e.g., pancreatic adenocarcinoma [PAAD], pheochromocytoma and paraganglioma [PCPG], prostate adenocarcinoma [PRAD]). In contrast, PTBP3 showed lower
expression in the tumours of kidney chromophobe (KICH), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), and thyroid carcinoma (THCA) (all p < 0.05) relative to the corresponding adjacent normal tissues.
After including the GTEx dataset, we further analysed the differences in PTBP3 expression between the tumour and normal tissues. We obtained significant difference for six tumour types (e.g., DLBC, LGG, OV, PAAD, TGCT, and THYM) (all p < 0.05) (Figure 1B). There were no significant differences in other tumours (e.g., adrenocortical carcinoma (ACC), bladder urothelial carcinoma (BLCA), acute myeloid leukaemia, PCPG, PRAD, sarcoma (SARC), and uterine carcinosarcoma (UCS) (Supplementary Figure S3). In general, we found that the expression of PTBP3 is elevated in most human tumours.
We also assessed the PTBP3 protein levels. The results of the National Cancer Institute’s CPTAC dataset showed that the total protein expression of PTBP3 was significantly higher in BRCA, COAD, UCEC, LUAD, and OV (all p < 0.05) than in normal tissues (Figure 1C). We further analysed the relationship between PTBP3 expression and tumour pathological staging and found stage-specific expressional changes in PTBP3 expression (e.g., ACC, KIRC, LIHC, Skin Cutaneous Melanoma (SKCM)) (Figure 1D).
Survival Analysis Data
To understand the relationship between PTBP3 expression and prognosis and OS, we divided tumour cases into high- and low- expression groups according to the median expression of PTBP3. TCGA and GEO datasets were used to analyse the differences in survival rates between the two groups. As shown in Figure 2A, the OS rate of the high-expression group was worse in the ACC, LGG, LIHC, LUSC, and PAAD (all p < 0.05). DFS analysis data showed that high PTBP3 expression was associated with poor prognosis for ACC, LUSC, and PAAD (Figure 2B). However, the OS and DFS rates of the high-expression group were better in the KIRC group (all p < 0.01).
We also used the Kaplan-Meier plotter tool to analyse the relationship between PTBP3 expression and OS, progression-free survival (PFS), first progression and post-progression survival, DFS, disease-specific survival (DSS), distant metastasis-free survival, and relapse-free survival (RFS). As shown in Supplementary Figure S4, the OS rate of the high PTBP3 expression group was poor in LUSC, PAAD, LIHC, and BRCA (all p < 0.05). The RFS rate of the high PTBP3 expression group was also worse in the PAAD, LIHC, and BRCA groups (all p < 0.05). However, the OS rate of the high PTBP3 expression group was better in the KIRC and STAD groups (all p < 0.01).
Genetic Alteration Analysis Data
Genetic alterations play important roles in tumorigenesis and development. Therefore, we explored PTBP3 genetic alterations in different tumour tissues. This study indicated that the frequency of PTBP3 alteration (>4%) was the highest in uterine tumours with “mutation” as the primary type. The “amplification” type was the primary type in the ACC, with a
A
Overall Survival
log10(HR)
ENSG00000119314.15
☐
1.0
(PTBP3)
ACC BLCA
0.5
BRCA
CESC
☐
CHOL
COAD
DLBC
ESCA
GBM
HNSC
KICH
☐ KIRP LAML
KIRC
☐ ☐ ☐
LGG
LIHC
LUAD
LUSC
MESO
☐
PAAD
PCPG
PRAD
READ
SARC SKCM
OV
STAD
TGCT
THCA
THYM
UCEC
UCS
UVM
0.0
-0.5
-1.0
1.0
ACC
1.0
KIRC
1.0
LGG
1.0
LIHC
1.0
LUSC
1.0
PAAD
p=0.00092
Percent
0.6
0.6
0.6
0.6
0.6
0.6
0.2
0.2
2
0.2
0.2
0.2
p=1e-04
p=1.2e-05
0
p=9.2e-05
p=0.028
p=0.0025
0
50
100 150
0
50
100
150
0
100
200
0
40
80
120
0
50
100
150
0
40
80
Months
Months
Months
Low PTBP3 Group
Months
Months
Months
High PTBP3 Group
B Disease Free Survival
log10(HR)
ENSG00000119314.15
(PTBP3)
☐
ACC
BLCA
BRCA
CESC
CHOL
COAD
DLBC
ESCA
GBM
HNSC
KICH
☐
KIRC
KIRP
LGG
LIHC
LUAD
☐
0.4
LUSC
MESO
☐
OV
PAAD
PCPG
PRAD READ
SARC
SKCM
STAD TGCT
THCA
THYM
UCEC
UCS
UVM
0.0
-0.4
1.0
ACC
1.0
KIRC
1.0
LUSC
1,0
PAAD
Percent
0.6
0.6
0.6
0.6
N
0.2
0.2
0,2
0
p=1.5e+06
p=0.0015
p=0.0048
p=0.025’
0
50
100
150
0 40 80 120
0
50 100 150
0
40
Months
Months
Months
Months
80
Low PTBP3 Group
High PTBP3 Group
FIGURE 2 | Relationship between PTBP3 expression level and survival (TCGA dataset). (A) Overall survival (GEPIA2) and (B) disease-free survival (GEPIA2). The positive results of survival map and Kaplan-Meier curves are listed. GEPIA2: Gene Expression Profiling Interactive Analysis, version 2.
frequency of approximately 3% (Figure 3A). Mutations and locations within PTBP3 are shown in Figure 3B. We found that “Missense” was the main type of genetic alteration and that R122*/Q alterations were detected in four cases of UCEC (Figure 3B). Figure 3C shows the 3D structure of the PTBP3 protein. Furthermore, we analysed the association between R122*/Q alterations and prognosis and found that UCEC cases with R122*/Q alterations showed better OS (p = 0.0415) and DFS (p = 0.0069), but not PFS (p=0.0792) and DSS (p = 0.499) (Figure 3D).
We also analysed the relationship between PTBP3 expression and tumour mutational burden (TMB) and microsatellite instability (MSI) in some tumours. PTBP3 expression was positively correlated with TMB in ACC, STAD, PAAD, LUAD, and SARC, but negatively correlated with THYM, UVM, KICH, THCA, and KIRC. Moreover, we observed a positive correlation between PTBP3 expression and MSI for READ, UCEC, STAD, CHOL, and MESO, and a negative correlation for DLBC, HNSC, UCS, PRAD, and ESCA (Supplementary Figure S5).
Protein Phosphorylation Analysis Data
Phosphorylation-dephosphorylation cascades are ubiquitous in tumour development. Thus, we compared PTBP3 phosphorylation levels between normal and tumour tissues. As shown in Figure 4, we found that the S30 locus exhibited a higher phosphorylation level in the primary tumour tissues of colon and ovarian cancers. The S426 locus also exhibited a higher phosphorylation level in the primary tumour tissues of UCEC tissues.
Immune Infiltration Analysis Data
PTBP3 mediates alternative splicing regulation of pre-mRNA and plays a role in the regulation of cell proliferation, differentiation, and migration. Tumor-infiltrating immune cells are linked to the initiation, progression, or metastasis of cancer. Thus, the TIMER, CIBERSORT, CIBERSORT-ABS, TIDE, XCELL, MCPCOUNTER, QUANTISEQ, and EPIC algorithms were used to analyse the relationship between PTBP3 expression and tumour-infiltrating immune cells. This study found a positive correlation between PTBP3
A
· Mutation . Structural Variant · Amplification
· Deep Deletion . Multiple Alterations
C
Alteration Frequency
4%
3%
2%
1%
Structural variant data
Mutation data CNA data
UCEC
BLCA
ACC
STAD
SKCM
SARC
CESC
CRAD
PRAD
ESCA
LGG
HNSC
BRCA
LUAD
OV
LUSC
THCA
KIRP
GBM
LIHC
KIRC
PTBP3
RefSeq: NM_005156 Ensembl: ENST00000374255
B
# PTBP3 Mutations
CCDS: CCDS6784 UniProt: PTBP3_HUMAN
5
R122*/Q UCEC(n=4)
0
RRM_5
RRM_1
RRM_5
PF14259
0
100
200
300
400
500
552aa
67
Missense
0
Splice
4
Splice
9
Truncating
0
SV/Fusion
4
SV/Fusion
D
OS
PFS
DFS
DSS
100%
100%
100%
100%
80%
80%
80%
80%
Overall
Progression Free
Disease Free
Disease-specific_
60%
60%
60%
60%
40%
Logrank P-Value:0.0415
40%
Logrank P-Value:0.0792
40%
Logrank P-Value:6.933e-3
40%
Logrank P-Value:0.499
20%
Altered group
20%
Altered group
Altered group
Altered group
0%
Unaltered group
0%
Unaltered group
20%
Unaltered group
20%
Unaltered group
20 60 100 140 180 220
20 60 100 140 180 220 Months
0%
0%
Months
20 60 100 140 180 220
Months
20 60 100 140 180 220 Months
expression and neutrophil cell count in most tumour types, especially in BLCA, CESC, DLBC, KIRC, LIHC, PRAD, and STAD. (Figure 5). The relationship between PTBP3 expression and cancer-associated fibroblasts was analysed. We found a positive correlation in most tumour types, especially in ACC, GBM, LIHC, and SARC, while a negative correlation was found in the TGCT tumour (Supplementary Figure S6).
We also analysed the relationship between PTBP3 expression and immune checkpoints (PDCD1, PDCD1LG2, CTLA4, CD274, HAVCR2, LAG3, and TIGIT). In 29 of 33 tumours, the expression of PTBP3 was positively correlated with the expression of immune checkpoint CD274. In 27 of 33 tumours, the expression of PTBP3 was positively correlated
with the expression of immune checkpoints PDCD1LG2, while a negative correlation was found in the tumour of THYM. A positive relationship between PTBP3 expression and all immune checkpoints was found in the tumours of UVM, STAD, SKCM, OV, LIHC, and LGG (Supplementary Figure S7).
N6-Methyladenosine Analysis Data
We also analysed the relationship between PTBP3 expression and N6-methyladenosine (mºA) modification. The mºA-related genes include ZC3H13, YTHDF1/2/3, YTHDC1/2, RBM15, METTL3/14, KIAA1429, IGF2BP1/2/ 3, WTAP, HNRNPD, HNRNPC, HNRNPA2B1, FTO, FMR1, ALKBH5 and so on. As shown in Supplementary
A
S30
Colon cancer Ovarian cancer
S426
Uterine corpus endometrial carcinoma
RRM
RRM
RRM
RRM
B
Protein expression of PTBP3 in Colon cancer Protein expression of PTBP3 in Ovarian cancer Protein expression of PTBP3 in UCEC
4-
(NP_001157260.1:S30)
2-
(NP_001157260.1:S30)
3.
(NP_001157260.1:S426) P=1.82E-08
2.
P=2.68E-28
1-
P=2.505527E-03
2.
Z-value
Z-value
Z-value
0-
0-
1
-1-
0.
-2.
-2-
-1.
-4.
-3-
-2.
-6
Normal (n=100)
Primary tumor (n=97)
-4
Normal (n=19)
Primary tumor (n=84)
-3
Normal (n=31)
Primary tumor (n=100)
CPTAC samples
CPTAC samples
CPTAC samples
FIGURE 4 | Protein phosphorylation of PTBP3 in different tumours (CPTAC dataset). We analysed the level of PTBP3 phosphoproteins (S30 and S426 sites) between normal and primary tumour tissue in colon cancer, ovarian cancer, and UCEC. (A) Phosphoprotein sites of PTBP3 are depicted in a schematic diagram. (B) PTBP3 phosphoprotein levels in colon cancer, ovarian cancer, and UCEC. CPTAC: Clinical Proteomic Tumor Analysis Consortium; UCEC: Uterine corpus endometrial carcinoma.
Figure S8, we found that the expression of most genes was increased in 32 types of tumours, especially in DLBC and UVM.
Enrichment Analysis of PTBP3-Related Partners
To further study the molecular mechanism of PTBP3 in tumorigenesis and development, we screened PTBP3- interacting proteins and PTBP3 expression correlated genes for a series of pathway enrichment analyses. We obtained 44 experimentally detected TWF1-binding proteins using the STRING website and the top 100 genes which correlated with PTBP3 expression (Figure 6A). As shown in Figure 6B, the expression of PTBP3 was positively correlated with that of pre-mRNA processing factor 40 homolog A (PRPF40A, R = 0.72), MOB kinase activator 1A (MOB1A, R = 0.68), family with sequence similarity 120A (FAM120A, R = 0.7), capping actin protein of muscle Z-line subunit alpha 1 (CAPZA1, R = 0.64), N-alpha- acetyltransferase 15 (NAA15, R = 0.68) (all p < 0.001), and we used these five genes to compute a heatmap that showed their positive association with PTBP3 (Figure 6C).
In this study, we combined the two datasets to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The GO data showed that “single-stranded RNA binding” and “mRNA 3’-UTR binding” were among the top hits. This might be involved in the effect of PTBP3 on tumour pathogenesis (Figure 7A). Figure 7B shows the interaction network of these genes and their molecular functions. KEGG pathway enrichment analysis showed that most genes were related to the “mRNA surveillance pathway” and “RNA degradation” (Figure 7C).
DISCUSSION
PTBP3, a member of the PTB family, has been reported to be involved in tumorigenesis and progression, such as in lung cancer (Wu et al., 2020), pancreatic cancer (Ma et al., 2020), breast cancer (Hou et al., 2018; Zhou et al., 2018; Liang et al., 2020), gastric cancer (Chen et al., 2014; Liang et al., 2018; Chen et al., 2020), colorectal cancer (Hou et al., 2019), and hepatocellular carcinoma (Yang et al., 2018). However, the oncogenic role of PTBP3 in various human tumours remains unclear. A pan-cancer analysis can identify the difference between normal and tumour tissues, providing a comprehensive understanding of the molecular mechanisms of tumorigenesis and progression (Chakravorty et al., 2019; Chai et al., 2020; Chen et al., 2021; Fiala et al., 2021; Taghvaei et al., 2021). Thus, we performed a pan-cancer analysis of the gene expression profile, survival status, genetic alteration, protein phosphorylation, immune infiltration, N6- methyladenosine, and relevant cellular pathways of PTBP3 in various human tumours based on the TCGA, GEO, and CPTAC databases.
This study showed that PTBP3 is highly expressed in most tumours. The expression levels of PTBP3 in tumours of BRCA, CESC, CHOL, COAD, ESCA, GBM, HNSC, LIHC, LUAD, LUSC, READ, STAD, and UCEC were higher than those in the corresponding adjacent normal tissues, whereas low expression was observed in KICH, KIRC, KIRP, and THCA. Furthermore, we found that PTBP3 overexpression generally predicted poor OS in patients with ACC, LGG, LUSC, and PAAD, in which ACC and LGG were reported for the first time. However, the OS and DFS rates of the high-expression group were better only in the KIRC group. These results suggest that PTBP3 is a potential biomarker for predicting the prognosis of multiple cancers.
A
B
Neutrophil
BLCA
Punty
Neutrophil_TIMER
Purity
Neutrophil_CIBERSORT
Purity
Neutrophil_CIBERSORT-ABS
AHO 750.000
DE-35
p > 0.05
BLCA
BLCA
BLCA
p < 0.05
CIBERSORT
CIBERSORT-ABS
ROD1 Expression Level (log2 TPM)
2.
2.
MCPCOUNTER
ROD1 Expression Level (log2 TPM)
ROD1 Expression Level (log2 TPM)
.
7.
*
Partial_Cor
QUANTISEQ
0.25
0.50
0.75
1.00 0.0
Infiltration Level
02
as
0.25
0.50
0.75
1.00 0.00
005 010 015
0.20
025
0.50
0.75
004
Purity
Purity
Infiltration Level
Purity
1.00.0.00
0.02
Infiltration Level
1
TIMER
XCELL
ROD1 Expression Level [log2 TPM)
Purity
Neutrophil_MCPCOUNTER
ROD1 Expression Level [log2 TPM)
Punty
Neutrophil_QUANTISEQ
ROD1 Expression Level (log2 TPM)
Purity
Neutrogini_XCELL
Foder
HIRO 4 COLOR
TRUM
2459
CHOSE
0
BLCA
BLCA
BLCA
-1
..
2.
.
”,
ACC (n=79)
0.25
0.50
0.75
1.00
50
100
150
200
250
0.25
0.50
0.75
1.00 0.00
0.05 0.90 0.15
0.20
0.25
0.30
0.75
0 04
Purity
0.05
Infiltration Level
Infiltration Level
1.00.0.00
0.02
Purity
Purity
Infiltration Level
BLCA (n=408)
ROD1 Expression Level (log2 TPM)
Punty
Neutrophil_TIMER
Purity
Neutrophil_CIDERSORT
Neutrophil_CIBERSORT-ADS
KIRC
ROD1 Expression Level (log2 TPM)
ARO -9 112
ROD1 Expression Level (og2 TPM)
Punty
TRAN
TEH
BRCA (n=1100)
O
BRCA-Basal (n=191)
KIRC
KIRC
KIRC
BRCA-Her2 (n=82)
BRCA-LumA (n=568)
0.25
0.50
0.75
0.3
0.9
0.25
0.75
0.03
0.00
1.09
0.12
0.25
0.50
0.75
0.025
0.dso
0.075
Purity
1.00 0.0
Infiltration Level
Purity
1.000.00
Infiltration Level
Purity
1.000.000
Infiltration Level
ROD1 Expression Level (log2 TPM)
Purty
Neutrophil_MCPCOUNTER
ROD1 Expression Level (log2 TPM)
Purity
Neutrophil_QUANTISEQ
Neutroghi_XCELL
BRCA-LumB (n=219)
28
ROD1 Expression Level (log2 TPM)
Purity
NA8-2
9H
67
PARMA
OH
CESC (n=306)
KIRC
KIRC
KIRC
CHOL (n=36)
COAD (n=458)
0.25
0.50
0.75
1.00
100
200
300
0.25
0.60
0.75
1.00 0.0
02
0.3
0.50
0.75
0 02
DLBC (n=48)
Purity
Infiltration Level
0.25
Purity
Infiltration Level
Infiltration Level
0.DE
Purity
1.00.0.00
LIHC
ROD1 Expression Level (log2 TPM)
Purty
Neutrophil_TIMER
ROD1 Expression Level (log2 TPM)
Purity
Neutrophil_CIBERSORT
ROD1 Expression Level (log2 TPM)
Purity
Neutrophil_CIBERSORT-ABS
ESCA (n=185)
THEH
+
MEHR
THEH
THAM
BRAH
HA
GBM (n=153)
LIHC
LIHC
₹
HNSC (n=522)
HNSC-HPV- (n=422)
P
0.25
0.50
6.75
1.000.0
21
0.50
15
0.05
Purity
Infiltration Level
0.4
0.25
9.75
1.000.00
0.10
0 20
0.25
0.50
Purity
0.75
Purity
1.00-0.00
5.10
Infiltration Level
Infiltration Level
0.15
HNSC-HPV+ (n=98)
ROD1 Expression Level (log2 TPM)
Purity
Neutrophil_MOPCOUNTER
Neutrophil_QUANTISEQ
PROSE
ROD1 Expression Level (log2 TPM)
Punty
Dos
THEH
P + 1:560-65
KICH (n=66)
KIRC (n=533)
LIHC
LIHC
KIRP (n=290)
LGG (n=516)
0.25
0.50
0.75
1.000
100
200
300
0.25
0.50
0.06 0.10 0.15
0 20
Purity
0.25
Infiltration Level
Purity
0.75
1.000.00
Infiltration Level
LIHC (n=371)
PRAD
ROD1 Expression Level (log2 TPM)
Purity
Neutrophd_TIMER
ROD1 Expression Level (log2 TPM)
Purity
Neutrophil_CIBERSORT
Purity
Neutrophil_CIBERSORT-ABS
COTTA
HA
PROGRA
ROD1 Expression Level (log2 TPM)
19969 32
BILE
TEAM
LUAD (n=515)
LUSC (n=501)
PRAD
PRAD
RA
A
MESO (n=87)
OV (n=303)
0.25
0.50
0.75
1.00 0:00
025
0.75
LO
0.25
0.50
1.000.000 0.005 0.010 0.015 0.000
0.25
0.25
Purity
Infiltration Level
0.75
0.50
0.002
0.0050
Purity
Infiltration Level
0.0100
Purity
1.00.0000
Infiltration Level
PAAD (n=179)
ROD1 Expression Level (jog2 TPM)
Purity
Neutrophil_MCPCOUNTER
ROD1 Expression Level (log2 TPM)
Purity
Neutrophil_QUANTISEQ
NHO
- 1000-65
PCPG (n=181)
PRAD (n=498)
PRAD
PRAD
4
READ (n=166)
SARC (n=260)
0.25
0.50
0.75
1.00
800
025
0.50
0.75
1.00 0.00
0.05
0 10
0 15
Purity
Infiltration Level
Purity
Infiltration Level
ROD1 Expression Level (log2 TPM)
Purity
Purity
Purity
SKCM (n=471)
STAD
2
Neutrophi_TIMER
ROD1 Expression Level (log2 TPM)
Neutrophi_CIDERSORT
ROD1 Expression Level (log2 TPM)
Neutrophil_CIBERSORT-ADS
Doch GT
Josde
· · @ 176-23
SKCM Meta(n=368)
SKCM Pri (n=103)
STAD
STAD
STAD
1
STAD (n=415)
.
.
0.25
0.50
Purity
0.75
1.00 0.0
02
0.4
0.6
0.26
0.75
1.000.00
0:06
15
@20
0.50
0.75
0.05
TGCT (n=150)
Infiltration Level
Purity
0.25
Infiltration Level
Purity
1.00.0.00
0.15
Infiltration Level
THCA (n=509)
ROD1 Expression Level (log2 TPM)
Purity
Neutrophil_MCPCOUNTER
Purity
Neutrophil_QUANTISEQ
HiGRO
ROD1 Expression Level (log2 TPM)
THAT
THYM (n=120)
UCEC (n=545)
STAD
STAD
UCS (n=57)
-
UVM (n=80)
0.25
0.50
Purity
0.75
1.00
100
200
300
0.25
0.75
Infiltration Level
Purity
1.000.00
0.05
0.10
Infiltration Level
0.15
FIGURE 5 | The correlation between PTBP3 expression level and infiltration of neutrophil (TCGA dataset). TIMER, CIBERSORT, MCPCOUNTER, QUANTISEQ and XCELL algorithms were used for the analysis.
B
8
8
A
KHDRBS3
RBMXL2
log2(PRPF40A TPM)
p-value = 0
p-value = 0
RBM42
R = 0.72
log2(MOB1A TPM)
R = 0.68
KHDRBS2
6
6
ELAVL1
SF1
ROMX
CHEK2
ELAL4
QKI
+
4
PABPCI
CHEK1
2
2
ELAVL3
RBMXL3
KHORBS1
POUR2A
NCL
O
.
O
.
PA2G4
ELAVL2
LFNG
0
2
4
log2(PTBP3 TPM)
6
8
0
2
4
6
8
RBMXL1
log2(PTBP3 TPM)
ILF3
PTBP3
PABPC1L2.
log2(FAM120A TPM)
8
p-value = 0
PABPC4
TXN2
R = 0.7
log2(CAPZA1 TPM)
8
p-value =. 0
R = 0.64
ILF2
PABPC3
6
6
PABPC4L
PABPC5
MFNG
NENF
4
4
CALU
PRKČDBP
PABPCIL
2
2
ENSG00000257411
NRG3
NOVA2
PABPC1L28
0
.
O
.
CDH2
0
2
4
6
8
0
2
4
6
8
PPWD1
NOVA1
RFNG
log2(PTBP3 TPM)
log2(PTBP3 TPM)
MPP1
MOSPD2
6
p-value = 0
log2(NAA15 TPM)
R = 0.68
5
+
3
~
-
0
.
0
2
4
6
8
C
log2(PTBP3 TPM)
UCEC (n=545)
THCA (n=509)
TGCT (n=150)
STAD (n=415)
SKCM-Primary (n=103)
SKCM-Metastasis (n=368)
UVM (n=80)
UCS (n=57)
THYM (n=120)
SKCM (n=471)
SARC (n=260)
READ (n=166)
PRAD (n=498)
PCPG (n=181)
PAAD (n=179)
OV (n=303)
MESO (n=87)
LUSC (n=501)
LUAD (n=515)
LIHC (n=371)
LGG (n=516)
KIRP (n=290)
KIRC (n=533)
KICH (n=66)
HNSC-HPV+ (n=98)
HNSC-HPV- (n=422)
HNSC (n=522)
GBM (n=153)
ESCA (n=185)
DLBC (n=48)
COAD (n=458)
CHOL (n=36)
CESC (n=306)
BRCA-LumB (n=219)
BRCA-LumA (n=568)
BRCA-Her2 (n=82)
BRCA-Basal (n=191)
Partial_Cor
1
BRCA (n=1100)
BLCA (n=408)
0
ACC (n=79)
-1
p > 0.05
p < 0.05
CAPZA1
FAM120A
MOB1A
NAA15
PRPF40A
FIGURE 6 | PTBP3-related genes. (A) Protein network map of experimentally determined PTBP3-binding proteins (STRING tool). (B) Expression correlation between PTBP3 and top 5 TWF1-correlated genes (CAPZA1, FAM120A, MOB1A, NAA15, and PRPF40A) (GEPIA2). (C) The heatmap data of PTBP3 and top 5 TWF1- correlated genes in human tumour types. GEPIA2: Gene Expression Profiling Interactive Analysis, version 2.
Our analysis yielded interesting findings. Especially in N6- methyladenosine analysis, we found that PTBP3 expression is related to all m6A markers in almost all human tumours, such as DLBC, PRAD, PAAD, LIHC, KIRP, THYM, and UVM. We speculated that PTBP3 protein is involved in N6- methyladenosine modification, or that it is an important component of methylase. In the immune infiltration analysis, we found that high PTBP3 expression was positively correlated with the cancer-associated fibroblast
and neutrophil infiltration levels in most tumours, such as ACC, GBM, LIHC, and SARC. However, its expression was negatively correlated with cancer-associated fibroblasts in TGCT and STAD. In addition, the expression of PTBP3 was positively correlated with the expression of immune checkpoints PDCD1LG2 in most tumours, while a negative correlation was only found in THYM tumours. Our findings may provide a novel clinical biomarker for predicting response to immunotherapy.
Pathway enrichment
A
C
single-stranded RNA binding
mRNA 3’-UTR binding
Spliceosome
p.adjust
poly-purine tract binding
0.0005
0.0010
RNA transport
-log10(P.value)
poly-pyrimidine tract binding
0.0015
0.0020
6
poly(A) binding
4
Count
poly(U) RNA binding
4
RNA degradation
2
6
8
Gene number
protein C-terminus binding
10
4
double-stranded RNA binding
12
14
Other types of O-glycan
6
8
mRNA 3’-UTR AU-rich region binding
16
biosynthesis
10
AU-rich element binding
mRNA surveillance pathway
0.050 0.075 0.100
GeneRatio
5.0
10.0
15.0
Fold enrichment
B
category
-AU-rich element binding
poly-pyrimidine tract binding
- double-stranded RNA binding mRNA
poly-purine tract binding
-3’-UTR AU-rich region binding mRNA
single-stranded RNA binding
3’-UTR binding
poly-purine tract binding
- poly-pyrimidine tract binding
poly(A) binding
poly(A) binding
poly(U) RNA binding
-poly(U) RNA binding
PRRC2C
mRNA 3’-UTR binding
-protein C-terminus binding
- single-stranded RNA binding
CD2AP
TOPBP1
double-stranded RNA binding
SP1
size
MKI67
4
mRNA 3’-UTR AU-rich region binding
8
POLR2A
12
NCL
16
EIF2AK2
AU-rich element binding
MTDH
DDX21
protein C-terminus binding
RC3H2
KHDRBS1
ILF2
KHDRBS2
ELAVL4
PABPC1
ELAVL3
PABPC1L
ELAVL2
PABPC1L2A
ELAVL1
PABPC1L2B
HNRNPF
PABPC3
RBMX
PABPC4
ILF3
PTBP1 PATL1
PABPC4L
PABPC5
SYNCRIP
FIGURE 7 | PTBP3-related gene enrichment. (A) Dotplot and (B) cnetplot of GO analysis. (C) KEGG pathway analysis. GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes.
A recent study suggested that PTBP3 promotes migration of non-small cell lung cancer (NSCLC) cells by regulating E-cadherin in the epithelial-mesenchymal transition signalling pathway and that there was a significant association between high
PTBP3 expression in NSCLC tissues with poor OS (Wu et al., 2020). Moreover, PTBP3 expression is associated with differentiation, lymph node metastasis, and distant metastasis (Wu et al., 2020). However, in this study, we explored LUAD
(TCGA, n = 515) and LUSC (TCGA, n = 501) and found a correlation between the high PTBP3 expression and poor OS prognosis (p = 0.0025) and poor DFS (P-0.0048) specific for LUSC but not for LUAD. This was an important discovery, and further research is needed.
Based on the TCGA database, the survival analysis results indicated a correlation between high PTBP3 expression and poor OS and DFS for PAAD. Ma et al. (2020) reported that high PTBP3 expression results in increased resistance to gemcitabine, which is attributed to attenuated autophagy. PTBP3 expression is higher in pancreatic cancer tissues than in corresponding adjacent normal tissues (Ma et al., 2020). The results showed that PTBP3 plays an important role in PAAD.
CONCLUSION
This study is the first to report a comprehensive pan-cancer analysis of PTBP3 in human tumours. PTBP3 was highly expressed in most tumours, and predicts poor survival. PTBP3 expression was significant correlated with the immune cell infiltration, TMB, MSI, PDCD1 and m6A-related markers. PTBP3 may be a new and promising prognostic and immunotherapy biomarker in human tumors. We believe that these findings may lay the groundwork for prospective functional experiments and may eventually have a positive impact in clinical practice.
DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
REFERENCES
Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cBio Cancer Genomics portal: an Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discov. 2, 401-404. doi:10. 1158/2159-8290.CD-12-0095
Chai, P., Yu, J., Ge, S., Jia, R., and Fan, X. (2020). Genetic Alteration, RNA Expression, and DNA Methylation Profiling of Coronavirus Disease 2019 (COVID-19) Receptor ACE2 in Malignancies: a Pan-Cancer Analysis. J. Hematol. Oncol. 13, 43. doi:10.1186/s13045-020-00883-5
Chakravorty, S., Yan, B., Wang, C., Wang, L., Quaid, J. T., Lin, C. F., et al. (2019). Integrated Pan-Cancer Map of EBV-Associated Neoplasms Reveals Functional Host-Virus Interactions. Cancer Res. 79, 6010-6023. doi:10.1158/0008-5472. CAN-19-0615
Chen, B., Chen, W., Mu, X., Yang, L., Gu, X., Zhao, A., et al. (2020). PTBP3 Induced Inhibition of Differentiation of Gastric Cancer Cells through Alternative Splicing of Id1. Front. Oncol. 10, 1477. doi:10.3389/fonc.2020.01477
Chen, B., Zhao, A. G., Shao, J., Mu, X. Y., Jiang, L., and Liu, J. W. (2014). The Effects of PTBP3 Silencing on the Proliferation and Differentiation of MKN45 Human Gastric Cancer Cells. Life Sci. 114, 29-35. doi:10.1016/j.lfs.2014.07.038
Chen, X. J., Ren, A., Zheng, L., Zheng, E. D., and Jiang, T. (2021). Pan-Cancer Analysis Identifies Liver Metastases as Negative Predictive Factor for Immune Checkpoint Inhibitors Treatment Outcome. Front. Immunol. 12, 651086. doi:10.3389/fimmu.2021.651086
AUTHOR CONTRIBUTIONS
ZF: Conceptualization (lead); Methodology (lead); Project administration (lead); Supervision (equal); Writing-original draft (lead). PL: Visualization (lead); Writing-original draft (equal); Revising (lead). HL: Investigation (lead); Conceptualization (equal); Methodology (equal); Writing- original draft (equal). WC: Methodology (lead); Resources (equal). LL: Writing-review and editing (supporting). LS: Supervision (lead); Validation (lead). FL: Conceptualization (lead); Funding acquisition (lead); Methodology (lead); Project administration (lead); Validation (lead); Writing-review and editing (lead).
FUNDING
This work was supported by the Beijing Municipal Science and Technology Commission (No. Z171100001017077), Municipal Administration of Hospitals Clinical Medicine Development of Special Funding Support (No. XMLX201404), Key Research and Development Program of Shandong Province (No.2019JZZY010104); Academic promotion programme of Shandong First Medical University (No. 2019QL021).
SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.811338/ full#supplementary-material
Chen, Y., Ji, Y., Liu, S., Liu, Y., Feng, W., and Jin, L. (2022). PTBP3 Regulates Proliferation of Lung Squamous Cell Carcinoma Cells via CDC25A-mediated Cell Cycle Progression. Cancer Cel Int 22, 19. doi:10.1186/s12935-022-02448-7
Cui, X., Zhang, X., Liu, M., Zhao, C., Zhang, N., Ren, Y., et al. (2020). A Pan-Cancer Analysis of the Oncogenic Role of Staphylococcal Nuclease Domain- Containing Protein 1 (SND1) in Human Tumors. Genomics 112, 3958-3967. doi:10.1016/j.ygeno.2020.06.044
Fiala, E. M., Jayakumaran, G., Mauguen, A., Kennedy, J. A., Bouvier, N., Kemel, Y., et al. (2021). Prospective Pan-Cancer Germline Testing Using MSK-IMPACT Informs Clinical Translation in 751 Patients with Pediatric Solid Tumors. Nat. Cancer 2, 357-365. doi:10.1038/s43018-021-00172-1
Gao, J., Aksoy, B. A., Dogrusoz, U., Dresdner, G., Gross, B., Sumer, S. O., et al. (2013). Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal. Sci. Signal. 6, pl1. doi:10.1126/scisignal.2004088
Hou, P., Chen, F., Yong, H., Lin, T., Li, J., Pan, Y., et al. (2019). PTBP3 Contributes to Colorectal Cancer Growth and Metastasis via Translational Activation of HIF-1a. J. Exp. Clin. Cancer Res. 38, 301. doi:10.1186/s13046-019-1312-y
Hou, P., Li, L., Chen, F., Chen, Y., Liu, H., Li, J., et al. (2018). PTBP3-Mediated Regulation of ZEB1 mRNA Stability Promotes Epithelial-Mesenchymal Transition in Breast Cancer. Cancer Res. 78, 387-398. doi:10.1158/0008- 5472.CAN-17-0883
Huo, G., Wang, Y., Chen, J., Song, Y., Zhang, C., Guo, H., et al. (2021). A Pan- Cancer Analysis of the Oncogenic Role of Twinfilin Actin Binding Protein 1 in Human Tumors. Front. Oncol. 11, 692136. doi:10.3389/fonc.2021. 692136
Liang, X., Chen, W., Shi, H., Gu, X., Li, Y., Qi, Y., et al. (2018). PTBP3 Contributes to the Metastasis of Gastric Cancer by Mediating CAV1 Alternative Splicing. Cell Death Dis 9, 569. doi:10.1038/s41419-018-0608-8
Liang, Y., Song, X., Li, Y., Chen, B., Zhao, W., Wang, L., et al. (2020). LncRNA BCRT1 Promotes Breast Cancer Progression by Targeting miR-1303/PTBP3 axis. Mol. Cancer 19, 85. doi:10.1186/s12943-020-01206-5
Ma, J., Weng, L., Jia, Y., Liu, B., Wu, S., Xue, L., et al. (2020). PTBP3 Promotes Malignancy and Hypoxia-Induced Chemoresistance in Pancreatic Cancer Cells by ATG12 Up-Regulation. J. Cel Mol Med 24, 2917-2930. doi:10.1111/jcmm. 14896
Marinescu, V., Loomis, P. A., Ehmann, S., Beales, M., and Potashkin, J. A. (2007). Regulation of Retention of FosB Intron 4 by PTB. PLoS One 2, e828. doi:10. 1371/journal.pone.0000828
Sharma, S., Kohlstaedt, L. A., Damianov, A., Rio, D. C., and Black, D. L. (2008). Polypyrimidine Tract Binding Protein Controls the Transition from Exon Definition to an Intron Defined Spliceosome. Nat. Struct. Mol. Biol. 15, 183-191. doi:10.1038/nsmb.1375
Szklarczyk, D., Gable, A. L., Nastou, K. C., Lyon, D., Kirsch, R., Pyysalo, S., et al. (2021). The STRING Database in 2021: Customizable Protein-Protein Networks, and Functional Characterization of User-Uploaded Gene/ measurement Sets. Nucleic Acids Res. 49, D605-D612. doi:10.1093/nar/ gkaa1074
Taghvaei, S., Sabouni, F., and Minuchehr, Z. (2021). Evidence of Omics, Immune Infiltration, and Pharmacogenomic for SENP1 in the Pan-Cancer Cohort. Front. Pharmacol. 12, 700454. doi:10.3389/fphar.2021.700454
Tang, Z., Kang, B., Li, C., Chen, T., and Zhang, Z. (2019). GEPIA2: an Enhanced Web Server for Large-Scale Expression Profiling and Interactive Analysis. Nucleic Acids Res. 47, W556-W560. doi:10.1093/nar/gkz430
Wu, Q., Zhang, B., Li, B., Cao, X., Chen, X., and Xue, Q. (2020). PTBP3 Promotes Migration of Non-small Cell Lung Cancer through Regulating E-Cadherin in EMT Signaling Pathway. Cancer Cel Int 20, 172. doi:10.1186/s12935-020-01240-9
Yang, X., Qu, S., Wang, L., Zhang, H., Yang, Z., Wang, J., et al. (2018). PTBP3 Splicing Factor Promotes Hepatocellular Carcinoma by Destroying the Splicing Balance of NEAT1 and Pre-miR-612. Oncogene 37, 6399-6413. doi:10.1038/s41388-018-0416-8
Zhou, Y., Zou, H., Wu, E., Huang, L., Yin, R., Mei, Y., et al. (2018). Overexpression of ROD1 Inhibits Invasion of Breast Cancer Cells by Suppressing the Translocation of ß-catenin into the Nucleus. Oncol. Lett. 16, 2645-2653. doi:10.3892/ol.2018.8917
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.
Copyright @ 2022 Fang, Li, Li, Chong, Li, Shang and Li. 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.