- 1Department of General Surgery, Tianjin Medical University General Hospital, Tianjin, China
- 2Department of Gastrointestinal Surgery, Wuhan Central Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Background: Piezo2 is a transmembrane-spanning ion channel protein implicated in multiple physiological processes, including cell proliferation and angiogenesis in many cell types. However, Piezo2 was recognized as representing a double-edged sword in terms of tumor growth. The prognostic and immunotherapeutic roles of Piezo2 in pan-cancer have not been reported.
Methods: In this study, several databases available including the UCSC Xena database, HPA, TIDE, GSEA, and cBioportal were used to investigate the expression, alterations, associations with immune indicators, and prognostic roles of Piezo2 across pan-cancer. R software and Perl scripts were used to process the raw data acquired from the UCSC Xena database.
Results: Based on processed data, our results suggested that Piezo2 expression levels were tissue-dependent in different tumor tissues. Meanwhile, the survival analysis reflected that patients suffering from KIRC, LUAD, and USC with high Piezo2 expression had good OS, while those suffering from KIRP and SARC with high Piezo2 expression had poor OS. In addition, our results showed that Piezo2 expression was associated with the infiltration of CD4+ T memory cells, mast cells, and dendritic cells. These results suggested that Piezo2 may involve tumor progression by influencing immune infiltration or regulating immune cell function. Further analysis indicated that Piezo2 could influence TME by regulating T-cell dysfunction. We also found that gene mutation was the most common genetic alteration of Piezo2. The GSEA analysis revealed that Piezo2 was associated with calcium ion transport, the activation of the immune response, antigen processing and presentation pathways.
Conclusion: Our study showed the expression and prognostic features of Piezo2 and highlighted its associations with genetic alterations and immune signatures in pan-cancer. Moreover, we provided several novel insights for further research on the therapeutic potential of Piezo2.
Introduction
As a common and frequently-occurring disease at present, cancer seriously threatens human health, reduces the quality of life, and brings a heavy economic burden to patients (Sung et al., 2021). With in-depth studies of cancer, we know that cancer initiation and progression are accompanied by multiple accumulations of genetic alterations, epigenetic alterations (Dawson and Kouzarides, 2012), cell metabolism (Justus et al., 2015), and cancer microenvironment changes (Wang et al., 2018). A growing body of evidence demonstrated that genetic alterations could drive genome instability and mutation, a hallmark of cancer. Hence, it is necessary to figure out the associations between genetic alterations and genome variation.
The tumor microenvironment (TME) mainly consists of tumor cells and various tissue components, including cancer-associated fibroblasts (CAF), microvessels, immune cells, cytokines, and chemokines (Wu and Dai, 2017)et al. Among them, immune cells, an essential component of TME, are recruited to the ecological niche in response to signals released by tumor cells. The anti-tumor effector functions of immune cells in TME were hampered or reversed into pro-tumor effector functions (Whiteside, 2008). Cancer cells implemented proliferation, migration, and invasion by escaping from immune surveillance. Thus, current cancer therapy is mainly focused on maintaining immune homeostasis and restoring immune surveillance by reactivating the capacity of immune cells (Galluzzi et al., 2020). Recently, emerging therapeutic targets and potential pathologic mechanisms have gradually been revealed. Many drugs targeting immune checkpoints, including PD-1, CTLA-4, TIM3, etc., have been approved for clinical application (Rotte, 2019). A study also has shown that angiopoietin-2 in melanoma was a potential predictive and prognostic biomarker for immune checkpoint therapy (Wu et al., 2017). Unfortunately, immune checkpoint inhibitors were effective in improving patient survival rates with multiple cancers, but they might bring adverse events to patients. And many patients still respond poorly to these therapies currently available. It is still essential to find novel targets for immune therapy and overcoming immune therapy resistance.
Piezo2, a member of the piezo family, is a cation-selective mechanical channel and mediates cation influx in which calcium ions (Ca2+) predominate (Coste et al., 2010). Piezo1 and Piezo2 played significant roles in regulating diverse physiological processes including sensing tactile sensations (Wang et al., 2019), regulating blood pressure (Zeng et al., 2018), prompting vascular development (Yang et al., 2016), and conducting mechanical force (Coste et al., 2010) et al. In addition to these, increasing shreds of evidence demonstrated that piezo protein could influence cancer development, progression, and invasion. During the past years, Piezo1 has been widely studied in cancers. Some studies have confirmed that Piezo1 could promote malignant biological behavior involving viability, migration, and metastasis of colon cancer cells (Sun et al., 2020; Jiang et al., 2021). Moreover, Piezo1 has been identified as a cancer-promoting factor in a variety of cancers, including hepatocellular carcinoma (Liu et al., 2021), gastric cancer (Wang et al., 2021), oral squamous cell carcinoma (Hasegawa et al., 2021), and breast cancer (Luo et al., 2022). However, the functional role of Piezo2 has not received much attention in cancers (Yang et al., 2016). Suggested that Piezo2 expression was elevated in glioma tumor cells, and Piezo2 downregulation inhibited the proliferative capacity of tumor cells and reduced tumor angiogenesis. As well, it has also been revealed that Piezo2 was associated with the progression of laryngeal squamous cell carcinoma and was a promising biomarker used in diagnosis and prognosis (Cheng et al., 2017). In contrast, a study has shown that Piezo2 expression levels were downregulated in non-small cell lung cancer (NSCLC) compared with normal tissues. Meanwhile, high Piezo2 expression was associated with better overall survival (OS) for patients (Huang et al., 2019). Likewise, a previous study reported that decreased expression of Piezo2 predicted poor prognosis of patients with breast cancer (Lou et al., 2019). Thus, Piezo2 might exert cancer-promoting or cancer‐suppressing effects in various cancers. It is, therefore, essential to explore the specific roles of Piezo2 in pan-cancer.
In the present study, we performed a comprehensive analysis to investigate the expression, genetic alteration, associations with immune indicators, and prognostic roles of Piezo2 across 33 cancers using the UCSC Xena, HPA, TIDE, GSEA, and cBioportal databases. Our study may further broaden the cancer landscape of Piezo2 and provide novel insights into cancer therapy.
Materials and methods
Data collection and processing
The gene expression data, clinical phenotype data, and mutation data of 33 cancers were downloaded from the UCSC Xena database (https://xena.ucsc.edu/). The database contained transcriptomic data of tumor samples and normal samples. The raw data were processed for ID conversion and other format conversion steps by R software (https://www.R-project.org) and Strawberry Perl. RNAseq data in FPKM format was converted into TPM format and log2 conversion was performed. The workflow of our study was shown in Figure 1.
The human protein atlas database
HPA (http://www.proteinatlas.org) is an open-access resource used for analyzing protein expression in tissues, single-cell types, immune cells, cell lines, and other subcellular levels. The protein subcellular localization, expression levels in cell lines, and expression levels in tissues of Piezo2 were obtained from the HPA database. The subcellular localization of Piezo2 was analyzed in the Hela cell line, Hep G2 cell line, and U-2 OS cell line. The relationship between the expression of Piezo2, endoplasmic reticulum (ER), microtubules, and nucleus was explored in the HPA database.
Relationships with pathological staging and prognostic value of Piezo2 in pan-cancer
Piezo2 mRNA expression differences and gene activity differences in 33 cancers were analyzed using “lima” and “GSVA” packages and visualized by “ggpubr” packages in R software. The gene activity scores of Piezo2 were measured by ssGSEA. The Kaplan-Meier method and univariate Cox regression analysis were used for survival analyses including OS, disease-free survival (DFS), disease-specific survival (DSS), and progression-free survival (PFS) across various cancers. Data analysis and graphical output were performed by “survival”, “survminer”, and “forestplot” packages.
To comprehensively evaluate the relationship between Piezo2 and clinical staging, we used the TCGA database and made a pathological stage plot of different cancer types. And the relationships between Piezo2 and T stage, N stage, and M stage were further analyzed separately. The data was processed and visualized by the “ggpubr” package in R software.
Analysis of the correlations between Piezo2 expression and the immune microenvironment in pan-cancer
The immune scores and stromal scores of samples were obtained using the ESTIMATE algorithm (Yoshihara et al., 2013) in R software. Calculation of immune scores and stromal scores were based on the different expression levels of specific genes in different immune and stromal cells. The correlation analysis between Piezo2 expression and immune/stromal scores was performed by “ggplot2”, “ggpubr”, and “ggExtra” packages. The relationships between Piezo2 expression and immune-related indicators including lymphocytes, immunoinhibitory, chemokines, and chemokine receptors were analyzed using the TISIDB database (http://cis.hku.hk/TISIDB/index.php).
Gene ontology and kyoto encyclopedia of genes and genomes
The GO and KEGG gene annotation files were downloaded from the GSEA website (http://www.gsea-msigdb.org/). According to the expression levels of Piezo2, the samples were divided into two groups with high and low expression and then matched with GO and KEGG annotation files using “limma”, “clusterProfiler”, “enrichplot”, and “ggplot2” packages in R software. The top 5 pathways and functions were displayed.
Analysis of epigenetic alterations
The data of tumor mutational burden (TMB) and microsatellite instability (MSI) in pan-cancer were downloaded as mentioned above and then integrated by Strawberry Perl. The correlations between Piezo2 expression and TMB/MSI were analyzed using the Spearman method, then visualized by the “fmsb” package in R software. The Piezo2 copy number alterations and genetic alterations in different cancers were analyzed by the cBioportal database (http://www.cbioportal.org/).
The tumor immune dysfunction and exclusion database
Tumor Immune Dysfunction and Exclusion (TIDE) database (Fu et al., 2020), which contained 33,197 samples and 189 cancer studies, was a predictor of immune checkpoint blockade therapy (http://tide.dfci.harvard.edu/). The Piezo2 predictive performance was calculated in TIDE compared with some published predictive markers in immune checkpoint blockade (ICB) cohorts. The associations between Piezo2 expression with ICB response outcomes and T cell dysfunction levels were evaluated in diverse cohorts. And then, the roles of Piezo2 methylation in T cell dysfunction and its risk scores were explored in various cancers.
Statistical analyses
All statistical analyses were performed using R (4.1.0) software and the abovementioned R packages. The statistical analyses were analyzed by Student’s t-test or ANOVA when appropriate. The Cox regression model and Kaplan-Meier method were used to performing survival analysis. p < 0.05 was considered statistically significant if not particularly specified.
Results
Piezo2 localization, composition, and expression of in normal, neoplastic cell lines and tissues
The essential functions of protein were always associated with its localization, structure, and sequence (Whisstock and Lesk, 2003). Most methods used in predicting protein function always rely on contrasting amino acid sequences and spatial structures with known proteins. To explore the intracellular protein localization of Piezo2, we retrieved its protein expression and distribution in the HPA database. Piezo2 was mainly detected in plasma membrane and cytoplasm, as shown in the simulation diagram (Figure 2A). As shown in Figure 2B, the distribution of Piezo2 was co-localized with the ER, microtubules, and nucleus in Hela, Hep G2, and U-2 OS cells. The protein topology structure of Piezo2 suggested that Piezo2 had 34 transmembrane domains which were the same as described in the literature (Figure 2C). Furthermore, we assessed the mRNA expression levels of Piezo2 in various normal tissues, including the digestive system, nervous system, cardiovascular system, etc., and multiple human tumor cell lines. We found that Piezo2 was widely expressed in a wide variety of normal human tissues and tumor cell lines (Figures 2D,E).
FIGURE 2. The localization, composition, and expression levels of Piezo2 in normal cells and tissues. (A) The expression of Piezo2 was mainly located at the plasma membrane and cytoplasm. (B) The expression and localization analyses of Piezo2, endoplasmic reticulum (ER), microtubules, and nucleus in Hela, Hep G2, and U-2 OS by immunofluorescence as adopted from the HPA database. (C) The protein topology structure of Piezo2. The structure showed the association between the protein structure of Piezo2 and the cell membrane. (D) Piezo2 mRNA expression levels in various normal tissues. The mRNA of Piezo2 was expressed in similar levels in most tissues. (E) Piezo2 mRNA expression levels were different in various tumor cell lines.
The mRNA and protein expression of Piezo2 in pan-cancer
To identify the associations between Piezo2 and cancers, we explored the expression differences of Piezo2 in tumor tissues. Cancer type abbreviations were shown in Supplementary Table S1. The mRNA expression levels of Piezo2 in tumor tissues were higher than in normal tissues in CHOL, HNSC, KIRC, LIHC, PCPG, STAD, and THCA (p < 0.05). On the contrary, Piezo2 expression levels in tumor tissues were lower than in normal tissues in BLCA, CESC, COAD, LUAD, etc. (Figure 3A). We sorted the expression levels of Piezo2 in tumor tissues and found that Piezo2 expression levels were the highest in MESO and lowest in CESC (Figure 3B). Furthermore, we evaluated the activity of Piezo2, which suggested that the gene activity of Piezo2 was grossly consistent with the trend of expression level (Figure 3C). Interestingly, the analysis results of Piezo2 gene activity were not the same as gene expression levels in different tumors (Figure 3D). Piezo2 gene activity was highest in KIRC and lowest in LAML. These results suggested that Piezo2 may have different co-expressed genes in different cancers. Further analysis of protein expression of Piezo2 in tumors showed that Piezo2 protein was also expressed in various tumor tissues, including the cerebral cortex, kidney, stomach, etc (Supplementary Figure S1).
FIGURE 3. Piezo2 expression levels and gene activity in 33 cancers. (A) Expression differences of Piezo2 in diverse cancers compared with normal tissues. Among them, Piezo2 was upregulated in CHOL, HNSC, KIRC, LIHC, PCPG, STAD, and THCA compared with normal tissues. In BLCA, COAD, GBM, KIRP, LUAD, LUSC, PRAD, and READ, Piezo2 was down-regulated in tumor tissues. (B) Relative expression differences of Piezo2 in diverse cancers. (C) Gene activity of Piezo2 in diverse cancers compared with normal tissues. Piezo2 had higher activity in tumor tissues of GBM, and KIRC. In BLCA, BRCA, CESC, CHOL, COAD, KICH, KIRP, LIHC, LUAD, LUSC, PRAD, READ, and UCEC, the gene activity of Piezo2 was higher in normal tissues. (D) Gene activity differences of Piezo2 in diverse cancers. *p < 0.05, **p < 0.01, ***p < 0.001.
Relationship between Piezo2 and clinical staging in pan-cancer
Owing to the abnormal expression of Piezo2 in tumor tissues, we accessed the relationship between Piezo2 expression and the clinical staging of tumors (Figure 4). As shown in Figure 4A, Piezo2 expression levels were associated with different clinical staging in BRCA, KIRC, KIRP, and THCA. We also analyzed the associations between Piezo2 expression and T, N, and M stages in depth. There was no apparent association between Piezo2 expression and tumor size in different tumors (Figures 4B–D). Meanwhile, we noticed that Piezo2 was upregulated in tumors with lymph node metastasis (Figures 4E–G). However, Piezo2 was down-regulated in tumors that occurred with distant metastases (Figures 4H,I).
FIGURE 4. The relationship between Piezo2 expression levels and tumor stage. (A) The relationship between Piezo2 mRNA expression levels and tumor stage in pan-cancer. (B–D) The relationship between Piezo2 mRNA expression levels and tumor size in STAD, LUADLUSC, and KIRC. (E–G) The relationship between Piezo2 mRNA expression levels and lymph node metastasis in KIRP, PAAD, and CRC. (H,I) The relationship between Piezo2 mRNA expression levels and distant metastasis in HNSC and KIRC. *p < 0.05, **p < 0.01, ***p < 0.001.
Correlations between Piezo2 and prognosis in pan-cancer
We further accessed the impact of Piezo2 expression on the prognosis of patients with multiple cancers. The survival analysis was verified by Cox regression and Kaplan-Meier analysis, respectively. The two methods were combined to explore the impact of Piezo2 expression on patient prognosis. As shown in Figures 5, 6, the relationship between the expression levels of Piezo2 and prognosis was not consistent in different cancers. The results of OS showed that Piezo2 was a protective factor in KIRC, LUAD, and UCS and an indicator of poor prognosis in KIRP and SARC (Figures 5A–F). The analyses of DFS revealed that Piezo2 acted as an unfavorable factor in KIRP, PAAD, and SARC (Figures 6A–E). Furthermore, the DSS analyses revealed that Piezo2 acted as a protective factor in KIRC, and LUAD and a risk factor in KIRP and PAAD (Figures 6F–J). We also explored whether Piezo2 expression affected the PFS of cancer patients. The results showed that increased Piezo2 expression was associated with a good prognosis for patients with KIRC and UCS (Figures 6K–O). While the PFS of patients from KIRP and PAAD with higher Piezo2 expression was shorter than those with lower expression. Thus, Piezo2 may play both cancer-promoting or cancer-suppressing functions in different tumor types.
FIGURE 5. Piezo2 expression levels were associated with OS in multiple cancers. (A) Cox regression of Piezo2 and OS in pan-cancer. (B–F) Kaplan-Meier analysis of Piezo2 expression and OS in KIRC, LUAD, KIRP, SARC, and UCS.
FIGURE 6. Piezo2 expression levels were associated with DFS, DSS, and PFS. Cox regression and Kaplan-Meier analysis of Piezo2 with DFS (A–E), DSS (F–J), and PFS (K–O).
Associations between Piezo2 and immune-related indicators
Growing evidence indicated that the tumor immune microenvironment may play an important role in the maintenance and progression of tumors. Hence, in the current study, we explored the associations between Piezo2 expression and immune signatures by scoring immune-related genes and stromal-related genes (Figure 7, Supplementary Figure S2). Overall, the associations between Piezo2 expression and immune/stromal scores showed a positive correlation. Among them, the strongest observed correlation between immune scores and Piezo2 was in PAAD (r = 0.45, p = 5.4e-10). And the strongest observed correlation between stromal scores and Piezo2 was in PRAD (r = 0.69, p < 2.2e-16). To further analyze whether Piezo2 can affect the immune environment, we explored the associations with other immune indicators, including lymphocytes, immunoinhibitory factors, chemokines, and chemokine receptors in TIBSD (Figure 8). The strongest indicators of positive and negative correlations were exhibited. Surprisingly, Piezo2 showed a positive correlation with memory B cells in PRAD (r = 0.629) and a negative correlation with activated CD8+ cytotoxic cells in TGCT (r = −0.516) (Figures 8A–C). In an analysis of immunoinhibitory factors, Piezo2 showed a strong positive correlation with KDR in numerous immunosuppressive factors (Figures 8D–F). And in the correlation analysis between Piezo2 expression and chemokine, Piezo2 showed the positively strongest correlation with CXCL12 (r = 0.573) in COAD and the negatively strongest correlation with CCL26 (r = −0.49) in TGCG (Figures 8G–I). In the observation of chemokine receptors, Piezo2 had the strongest positive correlation with CCR8 (r = 0.567) in PAAD and the strongest positive correlation with CCR10 (r = −0.389) in UVM (Figures 10J–L).
FIGURE 7. The associations between Piezo2 expression and immune scores in multiple cancers. Only results with p < 0.05 and |R| > 0.3 were considered for analysis.
FIGURE 8. Heatmap and scatter plot showed the correlations of Piezo2 expression with immune-related indicators, including lymphocytes (A–C), immunoinhibitory factors (D–F), chemokines (G–I), and chemokine receptors (J–L). The strongest representative of the positive and negative correlations of the indicators was shown in scatter plots.
Epigenetic alterations of Piezo2 and its associations with immunotherapy of cancers
We accessed the biomarker relevance of Piezo2 by comparing it with well-recognized biomarkers according to their predictive power of response outcomes of immune checkpoint blockade (ICB) therapy. The results showed that Piezo2 alone had an area under the receiver operating characteristic curve (AUC) of >0.5 in 6 of the 23 ICB cohorts (Figure 9A). Although Piezo2 was not as predictive as most “star molecules” such as CD8, CD247, MSI, etc., the predictive value of Piezo2 had the same effect as B. Clonality.
FIGURE 9. Piezo2 genetic alterations and associations with tumor immunity responses. (A) The bar chart of the receiver operating characteristic curve (AUC) showed the relevance of Piezo2 compared with published tumor immune evasion biomarkers. The size of the AUC value reflected predictive performances of the garget markers on the immune checkpoint blockade. “Custom” represented our target gene “Piezo2”. (B) The analyses of Piezo2 associations with T-cell dysfunction, outcomes using ICB treatment, and lymphocyte-mediated tumor killing. (C) Piezo2 alteration frequencies within various cancers. The mutation was the most common genetic variant of Piezo2. (D) The analysis of Piezo2 copy number alteration frequencies. Gene gain was the most frequent CNV of Piezo2. (E,F) The correlations of Piezo2 with TMB and MSI.
We also evaluated the relationship between Piezo2 expression levels and immunotherapy in multiple cancers (Figure 9B). Our results revealed that high Piezo2 expression correlated with worse outcomes of anti-CTLA4 treatment in melanoma (ICB_Vanllen2015_CTLA-4), anti-PD-1 treatment in melanoma (ICB_Gide2019 PD-1), anti-PD-L1 treatment in bladder cancer (ICB_Mariathasan2018_PD-L1) but achieved good outcomes in kidney cancer (ICB_Miao2018_ICB). Meanwhile, analysis of the T dysfunction value in the core dataset suggested that Piezo2 played an important role in breast cancer (METABRIC) and endometrial cancer (TCGA). At the same time, we found that Piezo2-knockout was an important factor of tumor-killing mediated by NK cells in melanoma (Freeman 2019 NK) and (Kearney 2018 NK_20).
Subsequently, we queried the epigenetic alterations of Piezo2, including genetic alterations frequencies, copy-number alterations, MSI, and TMB. We found that gene mutation was the most common genetic variant of Piezo2, followed by gene amplification and deep deletion, while a structural variant of Piezo2 was less frequent (Figure 9C). Moreover, the results of copy-number alterations (CNV) revealed that gene gain was the most frequent CNV of Piezo2, followed by gene shallow deletion, while deep deletion and amplification were less frequent (Figure 9D). In addition, we found that Piezo2 was related to MSI in HNSC, PAAD, STAD, and UCEC. And the results revealed that Piezo2 could impact TMB in BLCA, CESC, HNSC, PAAD, etc. (Figures 9E–11F). Furthermore, we assessed the effects of Piezo2 methylation on T-cell dysfunction in multiple cancers (Figure 10A). The results suggested that hypomethylation of Piezo2 was correlated with T-cell phenotypic dysfunction. After the discovery of correlations between Piezo2 hypomethylation and T-cell dysfunction, we evaluated whether hypomethylation of Piezo2 could affect the prognosis of cancer patients. Surprisingly, hypomethylation of Piezo2 reduced survival time among patients with lung cancer, uveal cancer, stomach cancer, and brain cancer. In contrast, hypomethylation of Piezo2 was associated with a good prognosis for patients with kidney cancer, bladder cancer, and breast cancer. (Figure 10B).
FIGURE 10. Epigenetic modification of Piezo2 influenced T cell phenotypes dysfunction, cytotoxic T-cell levels, and prognoses of different cancers. (A) The effects of Piezo2 methylation on T cell dysfunction, cytotoxic T-cell levels in cancers. (B) Kaplan-Meier curves of OS in patients with different methylation levels of Piezo2.
GSEA of Piezo2
Due to the ion channel properties of Piezo2, we considered whether Piezo2 could affect downstream signaling and thereby involve in the progression and metastasis of the tumor. As presented in Supplementary Figure S3, GO and KEGG analysis showed that Piezo2 was mainly enriched in cancer calcium signaling pathways. This was also in line with the calcium-gated properties of Piezo2. Surprisingly, we found that Piezo2 was also involved in the activation of immune response in OV, ACC, ESCA, LUAD, etc (Supplementary Figure S3A). Moreover, Piezo2 was enriched in antigen processing and presentation pathways (Supplementary Figure S3B). These indicated that Piezo2 may play some roles in immune processes.
Discussion
In the occurrence of malignant tumors, especially transformation and invasion, tumor cells were subject to mechanical stresses from the surrounding environment when they crossed the basement membrane and migrated to distant sites (Mierke, 2019). Piezo2, a member of the mechanosensitive ion channel piezo family, senses cell-matrix environment and mechanical changes, then regulates cation influx to influence cell biological behavior (Coste et al., 2010). Since it was discovered in 2010 (Coste et al., 2010), the piezo family which won the Nobel Prize in Physiology or Medicine in 2021 has gradually been studied and explored. The associations between Piezo1 and cancers have been extensively studied, such as colon cancer (Sun et al., 2020), breast cancer (Luo et al., 2022), and bladder cancer (Etem et al., 2018) et al., whereases the role of Piezo2 in pan-cancer has not been well studied. Therefore, in this article, we mainly analyzed the implications of Piezo2 in the prognosis and immunotherapy of multiple cancers using the UCSC Xena database and some online analysis tools.
In this study, we interrogated that the expression levels and epigenetic alterations of Piezo2 were associated with tumor staging, the prognosis of patients, and immunotherapy for various cancers. Piezo2 was mainly expressed in the cytoplasm, located on the cell membrane and co-localized with microtubules and endoplasmic reticulum in cells. To gain more knowledge about Piezo2 function, we explored Piezo2 mRNA and protein expression differences in multiple cell lines and tissues. Piezo2 presented a similar expression level in normal tissues. In subsequent detection of tumor cell lines, we found that Piezo2 was expressed in the kidney and urinary bladder, lung, mesenchymal, and lymphoid tumor cell et al. Through the analysis of TCGA data, the results showed that Piezo2 was up-regulated in tumor tissues from CHOL, HNSC, KIRC, LIHC, PCPG, STAD, and THCA compared with normal tissues. In contrast, Piezo2 was down-regulated in tumor tissues from BLCA, COAD, GBM, KIRP, LUAD, LUSC, PRAD and READ. Huang et al. (Huang et al., 2019) also found that Piezo2 mRNA expression levels were decreased in tumor tissues of LUAD and LUSC. Subsequent functional experiments showed that loss of Piezo2 could enhance the biological behavior of tumor cells including migration and proliferation ability, which suggested that Piezo2 may play an inhibitory function in tumor progression of non-small cell lung cancer (NSCLC). Although there was no statistical difference in the down-regulation of Piezo2 in tumor tissues of BRCA in our study, its gene activity was significantly reduced in tumor tissues. Owing to the correlation between low Piezo2 expression and poor prognosis in breast cancer (Lou et al., 2019), believed that Piezo2 could be used as a prognostic indicator in breast cancer. In contrast to our study, some scholars found that Piezo2 was highly expressed in GBM (Yang et al., 2016) and BLCA (Etem et al., 2018). Additional experiments were needed to prove the expression levels of Piezo2 in different cancers. These findings indicated that Piezo2 might play different functions in different types of tumors.
Although there were correlations between Piezo2 expression and clinical stage in multiple cancers such as BRCA, KIRC, KIRP, and THCA, we found no clear trend between different stages. For patients suffering from KIRP, PAAD, and CRC, the Piezo2 expression levels were lower in patients without lymph node metastases than those with lymph node metastases. In contrast, Piezo2 expression levels were relatively low for patients with distant metastases of HNSC and KIRC. Due to the inconsistency between the correlation of Piezo2 expression and various clinical stages in different cancers, we next examined the relationship between Piezo2 expression levels and the prognosis of patients with different cancers. Cox regression analysis and Kaplan-Meier analysis indicated that high Piezo2 expression suggested a good prognosis of OS for patients suffering from KIRC, LUAD, and USC and a poor prognosis for those with KIRP and SARC. With regard to the exploration of DFS, DSS, and PFS, we found that Piezo2 seemed to be a bad indicator for patients with KIRP and PAAD, whereas a good indicator for patients with KIRC. These findings were also consistent with the previous results of Piezo2 expression levels in tumor tissues.
In a variety of tumors such as lung cancer and bladder cancer (Cao et al., 2021; Wu et al., 2021), TME was used to assess the prognosis and efficacy of immune immunotherapy. We investigated immune scores and stromal scores to determine whether there was a correlation between Piezo2 and TME. In multiple cancers, the expression levels of Piezo2 were positively correlated with the immune scores and were also positively associated with the stromal scores. The immune cells in TME played either tumor-promoting or tumor-suppressing functions in angiogenesis, cell invasion, metastasis, etc. While cytotoxic T cells played a major role in the tumor microenvironment, resting CD4+ T memory cells also were potential response markers for precision detection (Zhang et al., 2020). In the construction and validation of lung cancer-related prognostic genes (Zhang et al., 2021), found that resting mast cells had the potential ability to predict individualized regimens of immunotherapy. And dendritic cells were also a promising target of immunotherapy, which were recruited to the extracellular matrix, diminished T cell functions, and promoted tumor progression (Wylie et al., 2019). Our results showed associations between Piezo2 and CD4+ T memory cells, mast cells, and dendritic cells, suggesting that Piezo2 may involve in tumor progression by influencing immune infiltration or regulating immune cell function. Interestingly, in our findings, Piezo2 showed a positive correlation with M2 macrophages and a negative correlation with regulatory T (Treg) cells in UVM. M2 macrophages and Tregs are both immunosuppressive cells, and play inhibitory effects in tumor progression. So, there were complex interactions between Piezo2 and the anti-tumor or pro-tumor response of immune cells. Furthermore, our study explored the correlations between Piezo2 expression with immunosuppressive genes, chemokines, and chemokine receptors. These results revealed that Piezo2 was closely related to indictors of immune infiltration, and provided novel ideas for immunotherapy. Recently, DNA methylation was demonstrated to induce T-cell depletion and dysfunction and was a barrier to restoring T-cell function through ICB therapy (Ghoneim et al., 2017). Therefore, we postulated that methylation statuses of Piezo2 might drive the transcriptomic alterations of cancer cells. It was intriguing to note that the methylation status of Piezo2 could mediate T-cell dysfunction and cytotoxic T-cell levels. And hypomethylation of Piezo2 mediated a longer life span for breast cancer, bladder cancer, and kidney cancer patients and a shorter life span for brain cancer, stomach cancer, uveal cancer, and lung cancer patients. These further confirmed that Piezo2 showed a tissue-dependent manner of regulating the immune environment and affecting prognosis in patients with cancer.
Immune checkpoint inhibitor treatments, which target the PD-1/PD-L1 axis, have drawn increasing attention in tumor treatment in recent years (Makuku et al., 2021). Our results revealed that kidney cancer with high expression of Piezo2 exhibited better clinical benefits of therapy against the PD-1 and PD-L1 axis. By contrast, higher Piezo2 expression was associated with worse outcomes of anti-CTAL4 or anti-PD-L1 treatment in melanoma and bladder cancer. These suggested that Piezo2 may exhibit opposite results targeting different immunotherapy checkpoints in different tumor environments. In addition, based on the reactions of Piezo2 to ICB, we evaluated the predictive power of Piezo2 compared with several well-recognized cancer immune evasion biomarkers. Surprisingly, while the predictive power of Piezo2 was not as good as most of the classic biomarkers, it was similar to B. Clonality. Recent studies have identified an essential role for B lymphocytes against anti-tumor immunotherapy across various cancers, including CRC (Berntsson et al., 2016), breast cancer (Garaud et al., 2019), and hepatocellular carcinoma (Shi et al., 2013). Considering the associations between T cell functions and calcium influx (Freitas et al., 2018), we speculated that Piezo2 might play certain roles in orchestrating TME through regulating calcium signaling. Ca2+ channels had important implications in the occurrence and progression of cancers, and considerable results have confirmed the potential roles of Ca2+ influx in tumor growth (Cantonero et al., 2019). In many types of cancer, calcium influx has been proven to exert effects on cancer cell behavior including cell proliferation (Capiod, 2011), apoptosis (Orrenius et al., 2003), and cell death (Trump and Berezesky, 1995) et al. We found that Piezo2 was significantly associated with calcium ion transport and calcium signaling pathway which was consistent with its ion channel properties. Furthermore, analysis of KEEG and GO in cancers suggested that Piezo2 could influence the activation of the immune response, antigen processing, and presentation pathways. These findings also verified Piezo2’s considerable impact on immune responses.
The changes in gene expression patterns caused by genetic alterations were thought to transform normal cells into malignant cells through hyperplastic and immortalization (Garnis et al., 2004). The analyses of gene alteration types that were already known for tumor suppressors or oncogenes should provide more insights into the functional roles in tumor advancement and metastasis (Hanahan and Weinberg, 2000). We explored the genetic alteration types and frequencies of Piezo2 in different tumors. We found that gene mutation was the most common genetic alteration of Piezo2, while structural variant was less frequent. Consequently, analysis of CNV suggested that gene gain and shallow deletion variants had the highest frequency of Piezo2. Most copy number variants were polymorphisms and were common in healthy people, although some copy number variants were pathogenic (Zack et al., 2013). In most instances, Piezo2 maintained its original diploid morphology and occurred with mild loss of copy number and amplification to become monoploid or triploid. The gene dosage imbalances caused by CNV led to a sequence of genetic abnormalities, finally promoting cancer tumorigenesis and progression (Wong et al., 2007; Wee et al., 2018). As mentioned above, more precise predictive biomarkers were needed for immunotherapy. And several reports have suggested that TMB could predict the therapeutic effect in multiple cancers (Rizvi et al., 2015; Van Allen et al., 2015; Samstein et al., 2019). Moreover, correlations between high TMB and the ICB therapeutic benefit were observed in some cohorts, including non-small cell lung cancer treated with anti-PD-1/L1 therapy (Rizvi et al., 2015) and melanoma treated with CTLA4 inhibitors (Van Allen et al., 2015) et al. Our study confirmed that Piezo2 expression was negatively correlated with TMB in BRCA, BLCA, STAD, KIRP, and other 11 cancers. Moreover, our study also revealed correlations between Piezo2 expression and MSI. These results confirmed that Piezo2 could influence TMB and MSI in cancers and impact immunotherapy response in patients, which provided new insights to improve therapy outcomes.
Although we have extensively explored the prognostic and immunological roles of Piezo2 expression in pan-cancer, there were some limitations to our analysis. First, there may be bias and heterogeneity in multiple analyses of different databases. Therefore, the contradictory results needed further verification. Second, external validation studies were warranted to confirm the associations between Piezo2 expression and tumor prognosis. Third, it would be better to confirm the relationship between Piezo2 and immune cells through simple experiments in vivo or in vitro.
Conclusion
To our knowledge, our study was the first to explore the roles of Piezo2 in pan-cancer. In brief, our analysis suggested that Piezo2 expression had cell-type-dependent and tissue-dependent in different types of cancer. Piezo2 expression levels were associated with the prognosis of patients with cancer. Piezo2 played underlying roles in the TME by regulating T-cell dysfunction and affecting immune-related factors. So, Piezo2 could be a novel target used for cancer treatment in the future.
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 author.
Author contributions
XL designed the study and performed the data analysis. YJ and ZZ wrote the manuscript. ZW and WF helped the revision. All authors read and approved the final manuscript.
Funding
This study was funded by Tianjin Key Medical Discipline (Specialty) Construction Project (TJSYXZDZK015).
Acknowledgments
We acknowledge the contributions from UCSC, TCGA, HPA, TIDE, TISIDB, and cBioPortal databases.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.918977/full#supplementary-material
References
Berntsson, J., Nodin, B., Eberhard, J., Micke, P., and Jirström, K. (2016). Prognostic impact of tumour-infiltrating B cells and plasma cells in colorectal cancer. Int. J. Cancer 139 (5), 1129–1139. doi:10.1002/ijc.30138
Cantonero, C., Sanchez-Collado, J., Gonzalez-Nuñez, M. A., Salido, G. M., Lopez, J. J., Jardin, I., et al. (2019). Store-independent Orai1-mediated Ca(2+) entry and cancer. Cell Calcium 80, 1–7. doi:10.1016/j.ceca.2019.02.012
Cao, R., Yuan, L., Ma, B., Wang, G., and Tian, Y. (2021). Tumour microenvironment (TME) characterization identified prognosis and immunotherapy response in muscle-invasive bladder cancer (MIBC). Cancer Immunol. Immunother. 70 (1), 1–18. doi:10.1007/s00262-020-02649-x
Capiod, T. (2011). Cell proliferation, calcium influx and calcium channels. Biochimie 93 (12), 2075–2079. doi:10.1016/j.biochi.2011.07.015
Cheng, L., Shen, Z., and Zhou, C. (2017). Promoter hypermethylation of PIEZO2 is a risk factor and potential clinical biomarker for laryngeal squamous cell carcinoma. Int. J. Clin. Exp. Pathol. 10 (12), 11635–11643.
Coste, B., Mathur, J., Schmidt, M., Earley, T. J., Ranade, S., Petrus, M. J., et al. (2010). Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels. Science 330 (6000), 55–60. doi:10.1126/science.1193270
Dawson, M. A., and Kouzarides, T. (2012). Cancer epigenetics: from mechanism to therapy. Cell 150 (1), 12–27. doi:10.1016/j.cell.2012.06.013
Etem, E., Ceylan, G. G., Özaydın, S., Ceylan, C., Özercan, I., and Kuloğlu, T. (2018). The increased expression of Piezo1 and Piezo2 ion channels in human and mouse bladder carcinoma. Adv. Clin. Exp. Med. 27 (8), 1025–1031. doi:10.17219/acem/71080
Freitas, C. M. T., Johnson, D. K., and Weber, K. S. (2018). T cell calcium signaling regulation by the Co-receptor CD5. Int. J. Mol. Sci. 19 (5), E1295. doi:10.3390/ijms19051295
Fu, J., Li, K., Zhang, W., Wan, C., Zhang, J., Jiang, P., et al. (2020). Large-scale public data reuse to model immunotherapy response and resistance. Genome Med. 12 (1), 21. doi:10.1186/s13073-020-0721-z
Galluzzi, L., Humeau, J., Buqué, A., Zitvogel, L., and Kroemer, G. (2020). Immunostimulation with chemotherapy in the era of immune checkpoint inhibitors. Nat. Rev. Clin. Oncol. 17 (12), 725–741. doi:10.1038/s41571-020-0413-z
Garaud, S., Buisseret, L., Solinas, C., Gu-Trantien, C., de Wind, A., Van den Eynden, G., et al. (2019). Tumor infiltrating B-cells signal functional humoral immune responses in breast cancer. JCI insight 5 (18), 129641. doi:10.1172/jci.insight.129641
Garnis, C., Buys, T. P., and Lam, W. L. (2004). Genetic alteration and gene expression modulation during cancer progression. Mol. Cancer 3, 9. doi:10.1186/1476-4598-3-9
Ghoneim, H. E., Fan, Y., Moustaki, A., Abdelsamed, H. A., Dash, P., Dogra, P., et al. (2017). De novo epigenetic programs inhibit PD-1 blockade-mediated T cell rejuvenation. Cell 170 (1), 142–157. e19. doi:10.1016/j.cell.2017.06.007
Hanahan, D., and Weinberg, R. A. (2000). The hallmarks of cancer. Cell 100 (1), 57–70. doi:10.1016/s0092-8674(00)81683-9
Hasegawa, K., Fujii, S., Matsumoto, S., Tajiri, Y., Kikuchi, A., and Kiyoshima, T. (2021). YAP signaling induces PIEZO1 to promote oral squamous cell carcinoma cell proliferation. J. Pathol. 253 (1), 80–93. doi:10.1002/path.5553
Huang, Z., Sun, Z., Zhang, X., Niu, K., Wang, Y., Zheng, J., et al. (2019). Loss of stretch-activated channels, PIEZOs, accelerates non-small cell lung cancer progression and cell migration. Biosci. Rep. 39 (3), BSR20181679. doi:10.1042/BSR20181679
Jiang, Y., Song, J., Xu, Y., Liu, C., Qian, W., Bai, T., et al. (2021). Piezo1 regulates intestinal epithelial function by affecting the tight junction protein claudin-1 via the ROCK pathway. Life Sci. 275, 119254. doi:10.1016/j.lfs.2021.119254
Justus, C. R., Sanderlin, E. J., and Yang, L. V. (2015). Molecular connections between cancer cell metabolism and the tumor microenvironment. Int. J. Mol. Sci. 16 (5), 11055–11086. doi:10.3390/ijms160511055
Liu, S., Xu, X., Fang, Z., Ning, Y., Deng, B., Pan, X., et al. (2021). Piezo1 impairs hepatocellular tumor growth via deregulation of the MAPK-mediated YAP signaling pathway. Cell Calcium 95, 102367. doi:10.1016/j.ceca.2021.102367
Lou, W., Liu, J., Ding, B., Jin, L., Xu, L., Li, X., et al. (2019). Five miRNAs-mediated PIEZO2 downregulation, accompanied with activation of Hedgehog signaling pathway, predicts poor prognosis of breast cancer. Aging 11 (9), 2628–2652. doi:10.18632/aging.101934
Luo, M., Cai, G., Ho, K. K. Y., Wen, K., Tong, Z., Deng, L., et al. (2022). Compression enhances invasive phenotype and matrix degradation of breast Cancer cells via Piezo1 activation. BMC Mol. Cell Biol. 23 (1), 1. doi:10.1186/s12860-021-00401-6
Makuku, R., Khalili, N., Razi, S., Keshavarz-Fathi, M., and Rezaei, N. (2021). Current and future perspectives of PD-1/PDL-1 blockade in cancer immunotherapy. J. Immunol. Res. 2021, 6661406. doi:10.1155/2021/6661406
Mierke, C. T. (2019). The matrix environmental and cell mechanical properties regulate cell migration and contribute to the invasive phenotype of cancer cells. Rep. Prog. Phys. 82 (6), 064602. doi:10.1088/1361-6633/ab1628
Orrenius, S., Zhivotovsky, B., and Nicotera, P. (2003). Regulation of cell death: the calcium-apoptosis link. Nat. Rev. Mol. Cell Biol. 4 (7), 552–565. doi:10.1038/nrm1150
Rizvi, N. A., Hellmann, M. D., Snyder, A., Kvistborg, P., Makarov, V., Havel, J. J., et al. (2015). Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science 348 (6230), 124–128. doi:10.1126/science.aaa1348
Rotte, A. (2019). Combination of CTLA-4 and PD-1 blockers for treatment of cancer. J. Exp. Clin. Cancer Res. 38 (1), 255. doi:10.1186/s13046-019-1259-z
Samstein, R. M., Lee, C. H., Shoushtari, A. N., Hellmann, M. D., Shen, R., Janjigian, Y. Y., et al. (2019). Tumor mutational load predicts survival after immunotherapy across multiple cancer types. Nat. Genet. 51 (2), 202–206. doi:10.1038/s41588-018-0312-8
Shi, J. Y., Gao, Q., Wang, Z. C., Zhou, J., Wang, X. Y., Min, Z. H., et al. (2013). Margin-infiltrating CD20(+) B cells display an atypical memory phenotype and correlate with favorable prognosis in hepatocellular carcinoma. Clin. Cancer Res. 19 (21), 5994–6005. doi:10.1158/1078-0432.CCR-12-3497
Sun, Y., Li, M., Liu, G., Zhang, X., Zhi, L., Zhao, J., et al. (2020). The function of Piezo1 in colon cancer metastasis and its potential regulatory mechanism. J. Cancer Res. Clin. Oncol. 146 (5), 1139–1152. doi:10.1007/s00432-020-03179-w
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA. Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Trump, B. F., and Berezesky, I. K. (1995). Calcium-mediated cell injury and cell death. FASEB J. 9 (2), 219–228. doi:10.1096/fasebj.9.2.7781924
Van Allen, E. M., Miao, D., Schilling, B., Shukla, S. A., Blank, C., Zimmer, L., et al. (2015). Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science 350 (6257), 207–211. doi:10.1126/science.aad0095
Wang, J. J., Lei, K. F., and Han, F. (2018). Tumor microenvironment: recent advances in various cancer treatments. Eur. Rev. Med. Pharmacol. Sci. 22 (12), 3855–3864. doi:10.26355/eurrev_201806_15270
Wang, L., Zhou, H., Zhang, M., Liu, W., Deng, T., Zhao, Q., et al. (2019). Structure and mechanogating of the mammalian tactile channel PIEZO2. Nature 573 (7773), 225–229. doi:10.1038/s41586-019-1505-8
Wang, X., Cheng, G., Miao, Y., Qiu, F., Bai, L., Gao, Z., et al. (2021). Piezo type mechanosensitive ion channel component 1 facilitates gastric cancer omentum metastasis. J. Cell. Mol. Med. 25 (4), 2238–2253. doi:10.1111/jcmm.16217
Wee, Y., Wang, T., Liu, Y., Li, X., and Zhao, M. (2018). A pan-cancer study of copy number gain and up-regulation in human oncogenes. Life Sci. 211, 206–214. doi:10.1016/j.lfs.2018.09.032
Whisstock, J. C., and Lesk, A. M. (2003). Prediction of protein function from protein sequence and structure. Q. Rev. Biophys. 36 (3), 307–340. doi:10.1017/s0033583503003901
Whiteside, T. L. (2008). The tumor microenvironment and its role in promoting tumor growth. Oncogene 27 (45), 5904–5912. doi:10.1038/onc.2008.271
Wong, K. K., deLeeuw, R. J., Dosanjh, N. S., Kimm, L. R., Cheng, Z., Horsman, D. E., et al. (2007). A comprehensive analysis of common copy-number variations in the human genome. Am. J. Hum. Genet. 80 (1), 91–104. doi:10.1086/510560
Wu, J., Li, L., Zhang, H., Zhao, Y., Zhang, H., Wu, S., et al. (2021). A risk model developed based on tumor microenvironment predicts overall survival and associates with tumor immunity of patients with lung adenocarcinoma. Oncogene 40 (26), 4413–4424. doi:10.1038/s41388-021-01853-y
Wu, T., and Dai, Y. (2017). Tumor microenvironment and therapeutic response. Cancer Lett. 387, 61–68. doi:10.1016/j.canlet.2016.01.043
Wu, X., Giobbie-Hurder, A., Liao, X., Connelly, C., Connolly, E. M., Li, J., et al. (2017). Angiopoietin-2 as a biomarker and target for immune checkpoint therapy. Cancer Immunol. Res. 5 (1), 17–28. doi:10.1158/2326-6066.CIR-16-0206
Wylie, B., Macri, C., Mintern, J. D., and Waithman, J. (2019). Dendritic cells and cancer: from biology to therapeutic intervention. Cancers 11 (4), E521. doi:10.3390/cancers11040521
Yang, H., Liu, C., Zhou, R. M., Yao, J., Li, X. M., Shen, Y., et al. (2016). Piezo2 protein: a novel regulator of tumor angiogenesis and hyperpermeability. Oncotarget 7 (28), 44630–44643. doi:10.18632/oncotarget.10134
Yoshihara, K., Shahmoradgoli, M., Martínez, E., Vegesna, R., Kim, H., Torres-Garcia, W., et al. (2013). Inferring tumour purity and stromal and immune cell admixture from expression data. Nat. Commun. 4, 2612. doi:10.1038/ncomms3612
Zack, T. I., Schumacher, S. E., Carter, S. L., Cherniack, A. D., Saksena, G., Tabak, B., et al. (2013). Pan-cancer patterns of somatic copy number alteration. Nat. Genet. 45 (10), 1134–1140. doi:10.1038/ng.2760
Zeng, W. Z., Marshall, K. L., Min, S., Daou, I., Chapleau, M. W., Abboud, F. M., et al. (2018). PIEZOs mediate neuronal sensing of blood pressure and the baroreceptor reflex. Science 362 (6413), 464–467. doi:10.1126/science.aau6324
Zhang, A., Yang, J., Ma, C., Li, F., and Luo, H. (2021). Development and validation of a robust ferroptosis-related prognostic signature in lung adenocarcinoma. Front. Cell Dev. Biol. 9, 616271. doi:10.3389/fcell.2021.616271
Keywords: Piezo2, pan-cancer, prognostic, genetic alterations, tumor microenvironment
Citation: Liu X, Jia Y, Wang Z, Zhang Z and Fu W (2022) A pan-cancer analysis reveals the genetic alterations and immunotherapy of Piezo2 in human cancer. Front. Genet. 13:918977. doi: 10.3389/fgene.2022.918977
Received: 10 June 2022; Accepted: 12 July 2022;
Published: 04 August 2022.
Edited by:
Jinhui Liu, Nanjing Medical University, ChinaReviewed by:
Tianxiang Lei, The First Affiliated Hospital of Sun Yat-sen University, ChinaQinfei Zhao, First Affiliated Hospital of Gannan Medical University, China
Taotao Liang, Henan Provincial Cancer Hospital, China
Copyright © 2022 Liu, Jia, Wang, Zhang and Fu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Weihua Fu, tjmughgs_fwh@163.com
†These authors have contributed equally to this work