Next Article in Journal
Chondrocytes In Vitro Systems Allowing Study of OA
Next Article in Special Issue
PAUF Induces Migration of Human Pancreatic Cancer Cells Exclusively via the TLR4/MyD88/NF-κB Signaling Pathway
Previous Article in Journal
Glycyrrhizic Acid Derivatives Bearing Amino Acid Residues in the Carbohydrate Part as Dengue Virus E Protein Inhibitors: Synthesis and Antiviral Activity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 18
10.3390/ijms231810316
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prognostic Value, Immune Signature, and Molecular Mechanisms of the PHLDA Family in Pancreatic Adenocarcinoma

by
Yunjie Duan
1,†,
Yongxing Du
1,†,
Zongting Gu
2,
Xiaohao Zheng
1 and
Chengfeng Wang
1,*
1
State Key Lab of Molecular Oncology and Department of Pancreatic and Gastric Surgery, National Cancer Center/National Clinical Research Center for Cancer/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, 17 Panjiayuan Nanli, Chaoyang District, Beijing 100021, China
2
Department of Hepatobiliary and Pancreatic Surgery and Minimally Invasive Surgery, Zhejiang Provincial People’s Hospital, Hangzhou Medical College, Hangzhou 310000, China
*
Author to whom correspondence should be addressed.
†
These authors have contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(18), 10316; https://doi.org/10.3390/ijms231810316
Submission received: 5 August 2022 / Revised: 28 August 2022 / Accepted: 2 September 2022 / Published: 7 September 2022
(This article belongs to the Special Issue Molecular Targets in Preventing Pancreatic Cancer Progression)
Browse Figures
Versions Notes

Abstract

:
Background: Increasing evidence supports the belief that the pleckstrin homology domain family A (PHLDA) family is associated with the development of a variety of cancers. However, the function of the PHLDA family members in PAAD is still unclear. Methods: Comprehensive bioinformatic analyses using R (version 3.6.3), Cytoscape (version 3.9.1), UALCAN, etc., were performed to study the clinicopathological characteristics, prognostic value, immune features, and functional mechanisms of the PHLDA family members in PAAD. Results: The PHLDA family members showed significantly elevated expression in PAAD compared with paracancerous or normal tissues. Their high expression or amplification were significantly correlated with worse clinicopathological characteristics and prognosis in PAAD patients. In addition, the role of the PHLDA family members in the immune regulation is diverse and complex. Mechanistically, TP53 mutations were significantly associated with the promoter methylation and expression levels of the PHLDA family members, which were activated in multiple oncogenic pathways, including the EMT, RAS/MAPK, and TSC/mTOR pathways. Moreover, we found that their expression levels were significantly correlated with the sensitivity of multiple traditional chemotherapeutic drugs and novel targeted MEK1/2 inhibitors. Conclusion: The PHLDA family members play an oncogenic role in the development of PAAD and might serve as new biomarkers or therapeutic targets.

1. Introduction

PAAD is a fatal digestive malignancy with a 5-year overall survival rate of approximately 7% [1,2]. Among cancer-related causes of death, this malignancy ranks fourth in the United States and sixth in China [1,3]. Approximately 80% of patients are diagnosed at a late stage with locally advanced or widespread tumors, thus missing the opportunity for radical resection [1,4]. Despite advances in imaging diagnostic techniques, such as endosonography, computed tomography, magnetic resonance imaging, and endoscopic retrograde cholangiography, the early diagnosis and differential diagnosis of PAAD remain challenging for clinicians [5]. At present, there are no specific tumor markers for the diagnosis of PAAD. Mucin-associated markers such as CA 19-9 are currently the most widely used tumor markers for PAAD, but their value as screening markers is limited by their low specificity [6]. In addition, over the last ten years, despite advances in antimetabolite therapy, targeted therapy, and immunotherapy for PAAD, the overall survival rate of patients has not significantly improved due to the late pathological stage at diagnosis, its highly invasive phenotype, and its resistance to chemotherapy [7]. Therefore, the need to explore new therapeutic targets and prognostic biomarkers to improve the prognosis of PAAD patients is urgent.
The pleckstrin homology-like domain (PHLD) protein is a multifunctional protein containing the PH domain that binds to the phosphorylated head group of phosphatidylinositol 4,5-bisphosphate (PtdI(4,5)P2) lipids with high in vitro specificity [8]. The PHLDA family members express the PHLD protein and consist of three members: PHLDA1, PHLDA2, and PHLDA3 [9]. The PHLDA1 gene was originally discovered in T-cell hybridomas and encodes an evolutionarily conserved protein that is rich in proline and histidine residues [10,11]. It has been shown that the PH domain of the PHLDA1 protein can inhibit Akt function by binding to phosphatidylinositol (PIP), which can inhibit tumor progression [12]. However, a study by Sakthianandeswaren et al. found that the expression of PHLDA1 also contributes to the development of intestinal tumors and is associated with the high metastatic potential of osteosarcoma cells [13,14]. PHLDA2 is the first apoptosis-related gene to be imprinted and expressed exclusively from the maternal allele during normal development [15,16]. PHLDA2 is located in the important tumor suppressor gene region 11p15.5 and is believed to have an anticancer effect in lung cancer [17,18]. However, its function in cancer is complex; for example, Moon HG et al. found that PHLDA2 knockdown could inhibit the invasion and proliferation of breast cancer [19]. PHLDA3 is a p53-regulated Akt repressor that competes with the PH domain of Akt for binding to membrane lipids, thereby inhibiting Akt translocation to the cell membrane and activation [20].
However, the study of Lei et al. found that PHLDA3 could also promote the proliferation and invasion of lung adenocarcinoma cells by activating the Wnt signaling pathway [21]. Taken together, the above studies indicate that the PHLDA family members may have different functions in different tumor tissues and cell types, and their specific effects depend on the tumor cells and environment [9]. Most importantly, the role of the PHLDA family members in PAAD is still unclear, and relevant studies are very scarce. This emphasizes the importance of carrying out bioinformatics analyses of the function of the PHLDA family in PAAD.
In the present study, we comprehensively investigated the possible functions and mechanisms of PHLDA family members in the development and progression of PAAD using public databases and multiple bioinformatics analysis techniques. First, we used The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases to analyze the differences in the expression of the PHLDA family members in PAAD samples compared with paracancerous samples and normal pancreas samples. We then utilized the HPA database to explore the protein levels of the PHLDA family members in PAAD tissues. In addition, we analyzed the relationship between their expression level and clinicopathological features, overall survival, and drug sensitivity in patients with PAAD based on the UALCAN database and Kaplan-Meier Plotter database. Potential mechanisms of involvement of PHLDA family members in the development and progression of PAAD were explored by gene variant, immune infiltration, gene enrichment, PPI, and promoter methylation analyses. In summary, our study revealed that the PHLDA family members play an oncogenic role in the occurrence and development of PAAD and are expected to become novel biomarkers and therapeutic targets for patients with this dismal malignancy.

2. Results

2.1. Aberrant Expression of the PHLDA Family Members across Cancers

To investigate the transcript levels of the PHLDA family members in a variety of cancerous tissues, adjacent paracancerous tissues, and normal tissues, we performed an analysis based on the data from the TCGA database and the GTEx database. The results showed that compared with paracancerous tissues and normal tissues, the transcription levels of the PHLDA family members were significantly increased in 29 types of cancer tissue, including PAAD, glioblastoma multiforme (GBM), and thyroid carcinoma (THCA), but were significantly decreased in 15 types of cancer tissue, including adrenocortical carcinoma (ACC), breast invasive carcinoma (BRCA), and colon adenocarcinoma (COAD) (Figure 1A,B,C). The above findings indicate that the expression of the PHLDA family members was different across cancers and may play diverse roles in different cancers.

2.2. Abnormal Expression of PHLDA1/3 Proteins in PAAD Tissues

To investigate the expression of the PHLDA family proteins in PAAD tissues, we used the data of the PAAD project from the HPA database to analyze the immunohistochemical staining of PHLDA1/3 proteins (the PHLDA2 data were missing). The cancer tissues of 10 PAAD patients were stained with the CAB016160 antibody, and the results showed that the PHLDA1 indicator was positive in 7 of the 10 tissue samples. The ratio of high staining–moderate staining–low staining–unstained was 1:4:2:3 (Figure 1D). In addition, the cancer tissues of 12 PAAD patients were stained with the HPA027601 antibody, and the results showed that the PHLDA3 indicators were positive in all of the 12 tissue samples. The ratio of moderate staining to low staining was 5:7 (Figure 1E). Based on the above results, we found that the PHLDA1 protein was positively stained in 70% of the samples, mainly with moderate and high staining, while the PHLDA3 protein was positively stained in all the samples. Thus, PHLDA1/3 proteins were both highly expressed in PAAD tissues, which was consistent with their expression at the transcriptional level.

2.3. The Relationship between PHLDA Family Member Expression and the Clinicopathological Characteristics and Prognosis of PAAD Patients

First, we analyzed the relationship between the expression of PHLDA family members and the clinicopathological characteristics of patients with PAAD using R (version 3.6.3) based on the PAAD project data in the UALCAN database. Our results showed that high PHLDA 2/3 expression was significantly associated with a history of pancreatitis (p = 1.838 × 10−2/p = 8.663 × 10−3) (Figure 2A), and the high expression of PHLDA1/2 was significantly correlated with a worse initial treatment effect (p = 0.003/p = 0.045) (Figure 2B). In addition, it is very interesting that PHLDA 1/3 is expressed at significantly higher levels in TP53 mutant PAAD tissues than in TP53 wild-type tumors (p = 1.642 × 10−2/p = 7.435 × 10−4) (Figure 2C). However, the expression levels of the PHLDA family members were not significantly correlated with any of these clinicopathological features, such as tumor grade, tumor stage, N stage, or diabetes history (Supplementary Figure S1).
Next, we analyzed the correlation of the expression of PHLDA family members and the prognosis of patients with PAAD using the Kaplan-Meier Plotter database. The results showed that the increased expression levels of PHLDA1 and PHLDA2 were significantly correlated with shorter OS and RFS (p < 0.05) (Figure 2D,E). In conclusion, the present study found that high expression of the PHLDA family members in PAAD patients was related to poor clinicopathological features and prognosis. This suggests that the PHLDA family members might play oncogenic roles and could serve as novel biomarkers in PAAD.

2.4. Mutations in the PHLDA Family Member Genes Were Associated with Worse Clinicopathological Characteristics and a Worse Prognosis for PAAD Patients

To further explore the related mechanisms by which PHLDA family members are differentially expressed in PAAD, we analyzed gene variants of PHLDA family members using the cBioPortal online tool. Genetic mutations occurred in the PHLDA family member genes in 26 samples (17%) from PAAD patients. Among them, the gene with the highest mutation frequency was PHLDA3 (9%), and the main types of mutations were amplification and mRNA high (Figure 3A). On this basis, we analyzed the clinicopathological characteristics and prognosis of PAAD patients with and without mutations in the PHLDA family member genes. The results showed that there were significant correlations between a high frequency of PHLDA2 mRNA mutations and two clinicopathological features of PAAD: the invasion of surrounding tissues and tumor sites (p < 0.05) (Figure 3B,C). Moreover, there was a significant correlation between the amplification mutations of PHLDA3 and a history of drinking in PAAD patients (p < 0.05) (Figure 3D). In terms of prognosis, PAAD patients with a high frequency of PHLDA1 mRNA mutations had shorter OS and DFS (p < 0.05) (Figure 3E,F), but PHLDA2/3 amplification mutations did not significantly affect the OS and DFS of PAAD patients (Supplementary Figure S2).

2.5. Correlations between the Methylation Levels of the PHLDA Family Member Promoters and the Clinicopathological Characteristics of PAAD Patients

We next explored the association of epigenetic variations in PHLDA family members with clinicopathological features of patients with PAAD. Based on the data of the PAAD project from the UALCAN database, we analyzed the relationship between the methylation levels of the PHLDA family member promoters and the clinicopathological characteristics of PAAD patients. As shown in Figure 3G, the promoter methylation levels of PHLDA1/3 in the tumor tissues of patients were significantly higher than those in the paracancerous tissues (p = 4.506 × 10−2/p = 6.353 × 10−5). However, the methylation level of the PHLDA2 promoter in paracancerous tissues was significantly higher than that in tumor tissues (p = 3.170 × 10−2). In addition, we found that the methylation level of the PHLDA2 promoter was significantly decreased in patients with TP53 mutations (p = 1.210 × 10−4) (Figure 3H), while the methylation level of the PHLDA3 promoter was significantly increased (p = 5.025 × 10−4) (Figure 3H). Therefore, the effect of TP53 mutations on the epigenetic modification of PHLDA2/3 promoter methylation may be an important factor for their diverse functions in the occurrence and development of PAAD. However, the methylation levels of the PHLDA family member promoters were not associated with the clinicopathological characteristics of tumor stage, tumor grade, or N stage (Supplementary Figure S3).

2.6. The Immune Landscape of PHLDA Family Member Expression and Variations in PAAD Patients

Next, using the data of the PAAD project from the TISIDB database, we explored the relationship between the expression of the PHLDA family members and the level of tumor-infiltrating cells and multiple immunomodulators based on various immunological markers in PAAD (Supplementary Material S1). Our results showed that the role of the PHLDA family members in tumor immune regulation was not consistent. First, PHLDA1 was positively correlated with the infiltration levels of a variety of immunocompetent tumor-infiltrating cells, including Act CD4, CD56dim, and Act DC (Table 1, Supplementary Figure S4A), and 23 types of immune promoters, including IL6, NT5E, and TNFSF9 (Figure 4A). Therefore, high expression levels of PHLDA1 may mainly enhance the immune response to PAAD. In contrast, PHLDA2 was negatively correlated with the infiltration levels of several immunocompetent tumor-infiltrating cells, including TemCD8, Th1, and NKT cells (Table 1, Supplementary Figure S4B), and PHLDA2 was negatively correlated with most immunomodulators (immune promoters, MHC molecules, chemokines, and chemokine receptors) in PAAD (Figure 4B–E). The above findings indicate that PHLDA2 mainly plays an immunosuppressive role in the PAAD tumor microenvironment, which might be associated with the immune escape of PAAD. However, PHLDA3 was only associated with the infiltration level of a small number of tumor-infiltrating cells (Table 1, Supplementary Figure S4C), indicating its limited immunomodulatory role in PAAD. In addition, we analyzed published single-cell sequencing datasets for pancreatic cancer (GSE148073 and GSE154778) and the final results show that the expression level of the PHLDA family is high in tumor tissues such as the TFF1 ductal cell and MET epithelial cell, but very low in immune infiltrating cells (Supplementary Figure S5). We speculate that the high expression of the PHLDA family on pancreatic cancer tumor cells may recruit relevant immune cells.
Then, we further explored the relationship between the copy number alterations of the PHLDA family members and the infiltration level of tumor-infiltrating cells based on the data from the PAAD project in the TIMER database. The results showed that copy number alterations of PHLDA1/3 may affect the infiltration level of six dominant immune cells, especially arm-level gain and high amplification (Figure 4F,G,H) (Supplementary Material S2).
Finally, we investigated the relationship between the expression levels of the PHLDA family members and immune subtypes. As shown in Figure 4I, the expression levels of the PHLDA family members were significantly different among the five types of immune subtypes. Among them, the highest expression level of the PHLDA family members was found in the PAAD tissues of the C1 subtype, while the lowest expression level of the PHLDA family members was found in the PAAD tissues of the C3 subtype.
In summary, the role of the PHLDA family members in the immune regulation of the PAAD tumor microenvironment is complex and diverse. It is necessary to pay special attention to the immune escape of PAAD cells caused by the high expression of PHLDA2, which may be a potential target for the treatment of PAAD.

2.7. Enrichment Analysis of the PHLDA Family Members and the 600 Co-Expressed Genes

To further explore the functional mechanism of the PHLDA family members in the development of PAAD, we performed GO and KEGG enrichment analyses (Supplementary Material S4) on the PHLDA family members and their 600 co-expressed genes (Supplementary Material S3), which were based on the data from the PAAD project in the LinkedOmics database and Metascape database. GO analysis showed that the PHLDA family members and their 600 co-expressed genes were mainly enriched in the BP terms “actin cytoskeleton organization”, “positive regulation of locomotion”, and “regulation of actin filament-based process” (Figure 5A), the CC terms “cell-substrate junction”, “cell-cell junction”, and “actin cytoskeleton” (Supplementary Figure S6A), and the MF terms “cadherin binding”, “protein domain specific binding”, and “kinase binding” (Supplementary Figure S6B). KEGG analysis showed that the PHLDA family members and their 600 co-expressed genes may play important roles in various pathways related to the occurrence and development of cancer, such as the “Hippo signaling pathway”, “apoptosis”, the “TNF signaling pathway”, and the “MAPK signaling pathway” (Figure 5B). On this basis, we used the data of the PAAD project from the GSCA database to further analyze the mechanism of action of the PHLDA family members in a variety of oncogenic pathways. As shown in Figure 5C, the PHLDA family played an active role in the EMT signaling pathway, and PHLDA1 was an important activator of the “RAS/MAPK signaling pathway”, “RTK signaling pathway”, and “TSC/mTOR signaling pathway”. In addition, the PHLDA family was mainly involved in the “DNA damage response signaling pathway”, “hormone AR signaling pathway”, and “PI3K/AKT signaling pathway”. The above findings reveal a possible mechanism by which the PHLDA family members may promote the development of pancreatic cancer, namely, by activating multiple oncogenic pathways and enhancing the survival and motility of PAAD cells.

2.8. Construction and Analysis of the PPI Network Associated with the PHLDA Family Members

To construct and analyze the PPI network of the PHLDA family members in PAAD patients, we used the STRING database to identify the 23 types of genes that have the strongest PPIs with the PHLDA family members (Figure 5D) (Supplementary Material S5) and then mapped the related PPI network using Cytoscape software. The larger circles and darker colors indicate that the number of PPIs associated with the gene was higher. The results showed that TP53 demonstrated the higher score in the PPI network of the PHLDA family members (Figure 5E) (Supplementary Material S6), suggesting that TP53 plays an important role in the PPI network closely related to the PHLDA family members. In addition, an interesting finding was that there was a direct PPI interaction between PHLDA3 and TP53 (Figure 5D), which suggests that PHLDA3 may be involved in the biological function of TP53 in PAAD.

2.9. The Expression Levels of the PHLDA Family Members Affect the Treatment Sensitivity of Multiple Drugs

Finally, we analyzed the associations between the expression levels of the PHLDA family members and the therapeutic sensitivity to various chemotherapeutic and targeted drugs using the GSCALite database. The results are shown in Figure 5F, and the expression levels of the PHLDA family members were negatively correlated with the sensitivity to various targeted or chemotherapeutic drugs for PAAD, including gemcitabine, docetaxel, and fluorouracil. Therefore, there is a good prospect for the PHLDA family members to serve as a therapeutic target for PAAD in drug development. In addition, an interesting finding was that the expression levels of PHLDA1/2 were positively correlated with the sensitivity to three targeted drugs, namely, selumetinib, trametinib, and PD318088, all of which are MEK1/2 inhibitors. We searched the ClinicalTrials.gov database (https://clinicaltrials.gov/) (accessed on 1 July 2022) for clinical studies related to the above three drugs and found that clinical trials of selumetinib and trametinib in the treatment of PAAD are in progress, but the clinical trial of PD318088 in the treatment of PAAD has not yet been carried out. As a promising biomarker for the prediction of the sensitivity to multiple drugs for PAAD, further experimental studies on the PHLDA family members should be carried out as soon as possible. The results of these studies are highly anticipated.

3. Discussion

In recent years, people have paid increasing attention to the function of the PHLDA family members in tumors. Relevant studies have found that PHLDA family members have potential tumor suppressor or cancer-promoting roles in various malignancies, such as breast, rectal, and lung cancer [13,21,22]. However, functional studies of PHLDA family members in PAAD are currently lacking, and no specific bioinformatics analysis has been performed. This study is the first bioinformatics analysis of the function of the PHLDA family members in PAAD and comprehensively reveals the functions and possible mechanisms of PHLDA family members involved in the development and progression of PAAD from the aspects of gene expression, gene variation, promoter methylation, immune cell infiltration, gene function, pathway enrichment, PPI and drug sensitivity.
As members of the PHLD protein family, the expression products of the PHLDA family members all contain the PH domain [9]. This domain can compete with the PH domain of Akt to bind membrane lipids to inhibit the activation of the Akt signaling pathway, thus playing an inhibitory role in a variety of tumors [9]. However, related studies have found that the PH domain can also promote the proliferation and migration of tumor cells [23]. Therefore, it is necessary to explore the function and molecular mechanism of the PH domain in PAAD. Interestingly, our study found that members of the PHLDA family can activate the EMT signaling pathway. The EMT signaling pathway is a classic cancer-promoting pathway that leads to the formation of secondary metastatic lesions by activating the movement and invasive ability of tumor cells in many tumors, including pancreatic cancer, prostate cancer, and breast cancer [24,25,26,27]. Therefore, we speculate that the PH domain of the PHLDA family expression products may be an important factor in activating the EMT signaling pathway, and this activation may be the potential mechanism by which the PHLDA family acts in promoting the occurrence and development of PAAD, which needs further experimental verification. In addition, PHLDA1 mainly played a stimulating role in the RAS/MAPK, RTK, and TSC/mTOR signaling pathways. The RAS/MAPK and RTK signaling pathways can promote the growth and proliferation of tumor cells [28,29], and the TSC/mTOR signaling pathway can promote tumor angiogenesis [30]. In short, the PHLDA family members may promote the occurrence and development of PAAD by activating a variety of carcinogenic pathways.
Next, it is worth noting that we discovered a possible target for the occurrence and development of PAAD driven by TP53 mutations. TP53, located on the short arm of human chromosome 17, is a tumor suppressor gene and the most common mutant gene in the human cancer cell [31]. In addition, TP53 mutations can be detected in 50–60% of human cancers, and some mutated proteins play a cancer-promoting role in tumors such as colon cancer, head and neck squamous cell carcinoma, and lung adenocarcinoma [32,33,34]. Thus, it is an important prognostic and predictive marker in cancer. In this study, we found that TP53 plays an important role in the PPI network of the PHLDA family members, and its expression products interact with a variety of proteins. Based on this rich interaction, the functional changes in proteins caused by TP53 mutations will inevitably affect the function of many interacting proteins and will further lead to the functional disorder of many interaction pathways [31]. Interestingly, this study also found that TP53 mutations were significantly correlated with the increased expression of PHLDA1 and PHLDA3 in PAAD, and the high expression of PHLDA1/3 was associated with poor prognosis in patients with PAAD. These findings suggest that PHLDA1/3 may be potential downstream targets of mutant TP53 that promote the occurrence and development of PAAD. This finding is helpful for researchers to further understand the mechanism by which TP53 mutations cause PAAD. Further study of this potential mechanism is expected to find a new target that can prevent the occurrence and development of PAAD driven by TP53 mutations, thus improving the prognosis of patients. Another concern is that we found that TP53 mutations may affect the promoter methylation levels of most members of the PHLDA family, and whether this effect can promote the occurrence and development of PAAD is worthy of future in-depth study. In addition, it is traditionally believed that CpG methylation in the promoter region is usually related to the inactivation of gene expression, so it is generally regarded as an inhibitory marker. However, recent studies have found that the function of DNA methylation in varied tissues and different gene regions is diverse [35,36,37]. Our study found that compared with TP53 wild-type PAAD tissues, the promoter methylation level and expression level of PHLDA3 in TP53 mutant PAAD tissues were both increased. These findings indicate that the methylation of the PHLDA3 promoter might increase its expression in PAAD.
Finally, we analyzed the role of the PHLDA family in the immunity of pancreatic cancer, and explored the relationship between the expression level of the PHLDA family and pancreatic cancer targeting or chemotherapeutic drug sensitivity. The final analysis results show that the PHLDA family is closely related to the clinical treatment of pancreatic cancer. Firstly, the immune system is complex. In addition to acting as a first line of defense against various pathogens, immune cells can provide surveillance functions by recognizing and destroying latent cancer cells. Immunocompetent tumor-infiltrating cells are a unique type of immune cells [38] that infiltrate the tumor microenvironment by detecting tumor antigens and releasing proinflammatory and immune molecules that regulate immune function. In addition, although cancer immunotherapy has been shown to improve survival in a variety of cancers, remission rates in patients with PAAD remain low [39]. Therefore, a comprehensive study of the interaction between tumors and immune cells will help to clarify the pathogenesis of cancer and develop a new immunotherapy strategy. We studied the relationship between the infiltration levels of immunocompetent tumor-infiltrating cells and the expression levels of the PHLDA family members in PAAD. The results show that the roles of PHLDA1 and PHLDA2 in PAAD immune regulation were opposite, which indicated that the roles of the PHLDA family members in PAAD tumor immune regulation are diverse and complex. Similarly, Bolandi et al. have also found that members of the B7 family may play different roles in tumor immune regulation [40]. This finding indicates that our incidence of members of the same family playing different roles in tumor immune regulation is not an isolated case, which may be very important to maintaining the balance between immune effectiveness and autoimmune suppression [40]. However, regardless of whether their regulatory functions are consistent, the functions of PHLDA1/2 in the tumor immunoregulation of PAAD are very active. For example, PHLDA1 was positively correlated with the infiltration levels of a variety of immunocompetent tumor-infiltrating cells, including Act CD4, CD56dim, and Act DC and 23 types of immune promoters, including IL6, NT5E, and TNFSF9. Therefore, high expression levels of PHLDA1 may mainly enhance the immune response to PAAD. In contrast, PHLDA2 was negatively correlated with the infiltration levels of several immunocompetent tumor-infiltrating cells, including TemCD8, Th1, and NKT cells, and PHLDA2 was negatively correlated with most immunomodulators (immune promoters, MHC molecules, chemokines, and chemokine receptors) in PAAD, so they can be used as potential targets of PAAD immunotherapy and molecular indicators to predict the efficacy of immunotherapy. In addition, our study revealed that the expression of the PHLDA family members was negatively correlated with the sensitivity to a variety of PAAD-targeting or chemotherapeutic drugs, including gemcitabine, docetaxel, and fluorouracil. Therefore, the research and development of PAAD-targeting drugs for the PHLDA family members will have good prospects. In addition, an interesting finding was that the expression levels of PHLDA1/2 were positively correlated with the sensitivities to three MEK 1/2 inhibitors, including selumetinib, trametinib, and PD318088. Previous studies have found that the activation of the MEK/ERK/P21 signaling pathway can promote the development of PAAD [41], and MEK inhibitors can activate the immune recognition of PAAD to inhibit the progression of PAAD [42]. In addition, clinical trials have confirmed that trametinib inhibits PAAD tumors [43]. Therefore, for PAAD patients with high expression of PHLDA1/2, it may be better to choose MEK1/2 inhibitors for treatment.
This research is the first bioinformatics analysis of the functions of the PHLDA family members in PAAD. We found obvious differences in the expression levels of the PHLDA family members across cancers, and these proteins may play diverse roles in different cancers. The high levels of expression and amplification mutations of the PHLDA family members are significantly related to worse clinicopathological features and a worse prognosis for PAAD patients. In addition, the PHLDA family members may promote the occurrence and development of PAAD by activating a variety of oncogenic pathways and regulating PAAD tumor immunity. Another finding is that TP53 mutations have a significant effect on the promoter methylation levels of most members of the PHLDA family, and PHLDA1/3 may be potential downstream targets to drive the occurrence and development of PAAD with TP53 mutations. Moreover, we found that the PHLDA family members are associated with the sensitivity to a variety of PAAD-targeted or chemotherapeutic drug treatments. Therefore, the PHLDA family members could be used as targets for the development of therapeutic drugs for PAAD or as molecular markers for predicting drug sensitivity. However, the present study has certain limitations. For example, the number of databases included in this study is somewhat inadequate. Moreover, the present study is just a bioinformatics analysis of the functions and mechanisms of the PHLDA family involved in the development and progression of PAAD. Future experimental studies could further confirm the tumor-promoting role of PHLDA family members in PAAD.

4. Methods and Materials

Ethics statement: This study was approved by the Ethics Committee of the Cancer Hospital, Chinese Academy of Medical Sciences and was conducted strictly according to the principles of the Declaration of Helsinki. All data in this study were retrieved from online databases. No human or animal tests were involved.

4.1. Expression Analysis

In this study, R (version 3.6.3) was used to analyze the expression levels of the PHLDA family member genes in cancerous and paracancerous tissues in the PAAD project of the TCGA database and normal pancreas of the GTEx database (https://commonfund.nih.gov/GTEx) (accessed on 2 July 2022). The statistical significance of the difference in expression was evaluated by the Wilcoxon test. Then, the UALCAN database (http://ualcan.path.uab.edu) (accessed on 4 July 2022) was used to obtain the expression levels of the PHLDA family member genes in PAAD patients with different clinicopathological characteristics, and the statistical analysis was performed by Welch’s t test. p < 0.05 was considered statistically significant [44,45].

4.2. HPA Analysis

The immunohistochemical staining of PHLDA1/2 in PAAD tissues was analyzed using the HPA database (https://www.proteinatlas.org/) (accessed on 5 July 2022). According to dyeing intensity and quantification, the protein expression levels were classified into four categories. The classification criteria for the protein expression levels were as follows: negative, not detected; weak and <25%, not detected; weak combined with either 25–75% or 75%, low; moderate, and <25%, low; moderate combined with either 25–75% or 75%, medium; strong and <25%, medium; and strong combined with either 25–75% or 75%, high [46].

4.3. Survival Analysis

The Kaplan-Meier Plotter database (http://www.kmplot.com/) (accessed on 6 July 2022) was used to analyze the correlations between PHLDA family member expression levels and overall survival (OS) and recurrence-free survival (RFS) in PAAD [47]. The “survival” R package (version 2.38) was utilized to calculate log-rank P values, and p < 0.05 was considered statistically significant [48].

4.4. Gene Variation Analysis

The cBioPortal database (http://www.cbioportal.org/) (accessed on 3 July 2022) was used to analyze the genetic variations of the PHLDA family members in PAAD, and the correlations between the variations and some clinicopathological characteristics were further determined [49]. Statistical significance was assessed by the chi-squared test, and p < 0.05 was considered statistically significant. In addition, the log-rank test was used to evaluate the statistical significance between PHLDA family member gene mutations and OS and disease-free survival (DFS) in PAAD patients, and p < 0.05 was considered statistically significant [50].

4.5. Methylation Analysis

To evaluate the associations between the expression levels of the PHLDA family members and the levels of promoter methylation, promoter methylation levels in PAAD under different conditions were analyzed in the UALCAN database. The statistical significance of the difference was evaluated by Welch’s t test, and p < 0.05 was considered statistically significant [51,52].

4.6. Immune Infiltration Analysis

The correlations between the expression levels of the PHLDA family members and the infiltration levels of infiltrating immune cells and the expression levels of immune molecules in PAAD were analyzed using the TISIDB database (http://cis.hku.hk/TISIDB/) (accessed on 7 July 2022). Statistical significance was assessed by the Spearman test, and p < 0.05 was considered statistically significant [53]. The database was then used to analyze the differences in the expression levels of the PHLDA family members in different immune subtypes of PAAD. The statistical significance of the difference in expression levels was evaluated by the Kruskal-Wallis test, and p < 0.05 was considered statistically significant [53]. In addition, the TIMER database (http://timer.cistrome.org/) (accessed on 8 July 2022) was used to analyze the correlations between the copy number alterations of the PHLDA family members and the infiltration levels of the six types of immune cells in PAAD. The statistical significance of the difference was evaluated by the Wilcoxon rank-sum test, and p < 0.05 was considered statistically significant [54,55].

4.7. Gene Enrichment Analysis

The LinkedOmics database (http://www.linkedomics.org/) (accessed on 10 July 2022) was used to select the 600 genes most closely related to the co-expression of the PHLDA family members [56]. The Metascape database (https://metascape.org) (accessed on 14 July 2022) was used to visualize the biological process (BP), cellular components (CC), molecular function (MF), and Kyoto Encyclopedia of Genes and Genomes (KEGG) of the PHLDA family members and their 600 co-expressed genes [57]. In addition, pathway enrichment of the PHLDA family members was performed using the GSCALite database (http://bioinfo.life.hust.edu.cn/web/GSCALite/) (accessed on 16 July 2022) [58].

4.8. Construction of the Functional PPI Network

The genes with the strongest PPIs with the PHLDA family members were obtained by using the STRING database (https://string-db.org/) (accessed on 18 July 2022), and the functional PPI network of the PHLDA family members was established [59]. Cytoscape (version 3.9.1) was used to score the effect of different genes in the PPI network [60,61].

4.9. Drug Sensitivity Analysis

The GSCALite database was used to analyze the correlation between the expression levels of the PHLDA family members and a variety of chemotherapy or targeted therapy drugs for PAAD. The statistical significance of the differences was assessed by the Spearman test. p < 0.05 was considered statistically significant [58].

5. Conclusions

In summary, this comprehensive study on the functions of the PHLDA family members reveals that they might play an oncogenic role in the development and progression of PAAD and serve as new biomarkers or therapeutic targets for patients with PAAD. Further functional verification and mechanism exploration studies should be promptly carried out.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms231810316/s1.

Author Contributions

Y.D. (Yunjie Duan) and Y.D. (Yongxing Du) conceived and designed the study. Y.D. (Yunjie Duan), Z.G. and X.Z. performed the analyses. Y.D. (Yunjie Duan) and Y.D. (Yongxing Du) wrote the manuscript. C.W. supervised the study. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Natural Science Foundation of China (81972314 and 81802463) and Clinical Research Fund of Jieping Wu Medical Foundation (320.6750.2020-19-4).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials; further inquiries can be directed to the corresponding author.

Acknowledgments

We thank Y.D., Z.G., X.Z. and C.W. from the department of Pancreatic and Gastric Surgery, National Cancer Center/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College for their kind technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

PAAD: pancreatic adenocarcinoma; PHLDA: pleckstrin homology domain family A; PPI: protein–protein interaction; PHLD: pleckstrin homology-like domain; PIP: phosphatidylinositol; TCGA: The Cancer Genome Atlas; GTEx: Genotype-Tissue Expression; OS: overall survival; RFS: recurrence-free survival; DFS: disease-free survival; BP: biological process; CC: cellular components; MF: molecular function; KEGG: Kyoto Encyclopedia of Genes and Genomes; GBM: glioblastoma multiforme; THCA: thyroid carcinoma; ACC: adrenocortical carcinoma; BRCA: breast invasive carcinoma; COAD: colon adenocarcinoma; PtdI(4,5)P2: phosphatidylinositol 4,5-bisphosphate.

References

  1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2019. CA Cancer J. Clin. 2019, 69, 7–34. [Google Scholar] [CrossRef] [PubMed]
  2. Moore, A.; Donahue, T. Pancreatic Cancer. JAMA 2019, 322, 1426. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, W.; Zheng, R.; Baade, P.D.; Zhang, S.; Zeng, H.; Bray, F.; Jemal, A.; Yu, X.Q.; He, J. Cancer statistics in China, 2015. CA Cancer J. Clin. 2016, 66, 115–132. [Google Scholar] [CrossRef] [PubMed]
  4. Allemani, C.; Matsuda, T.; Di Carlo, V.; Harewood, R.; Matz, M.; Nikšić, M.; Bonaventure, A.; Valkov, M.; Johnson, C.J.; Estève, J.; et al. Global surveillance of trends in cancer survival 2000–2014 (CONCORD-3): Analysis of individual records for 37,513,025 patients diagnosed with one of 18 cancers from 322 population-based registries in 71 countries. Lancet 2018, 391, 1023–1075. [Google Scholar] [CrossRef]
  5. Adamek, H.E.; Albert, J.; Breer, H.; Weitz, M.; Schilling, D.; Riemann, J.F. Pancreatic cancer detection with magnetic resonance cholangiopancreatography and endoscopic retrograde cholangiopancreatography: A prospective controlled study. Lancet 2000, 356, 190–193. [Google Scholar] [CrossRef]
  6. Ringel, J.; Löhr, M. The MUC gene family: Their role in diagnosis and early detection of pancreatic cancer. Mol. Cancer. 2003, 2, 9. [Google Scholar] [CrossRef]
  7. Cui, X.H.; Hu, S.Y.; Zhu, C.F.; Qin, X.H. Expression and prognostic analyses of the insulin-like growth factor 2 mRNA binding protein family in human pancreatic cancer. BMC Cancer 2020, 20, 1160. [Google Scholar] [CrossRef]
  8. Haslam, R.J.; Koide, H.B.; Hemmings, B.A. Pleckstrin domain homology. Nature 1993, 363, 309–310. [Google Scholar] [CrossRef]
  9. Fuselier, T.T.; Lu, H. PHLD Class Proteins: A Family of New Players in the p53 Network. Int. J. Mol. Sci. 2020, 21, 3543. [Google Scholar] [CrossRef] [PubMed]
  10. Park, C.G.; Lee, S.Y.; Kandala, G.; Lee, S.Y.; Choi, Y. A novel gene product that couples TCR signaling to Fas(CD95) expression in activation-induced cell death. Immunity 1996, 4, 583–591. [Google Scholar] [CrossRef] [Green Version]
  11. Bonatto, N.; Carlini, M.J.; de Bessa Garcia, S.A.; Nagai, M.A. PHLDA1 (pleckstrin homology-like domain, family A, member 1) knockdown promotes migration and invasion of MCF10A breast epithelial cells. Cell Adhes. Migr. 2018, 12, 37–46. [Google Scholar] [CrossRef] [PubMed]
  12. Chen, Y.; Takikawa, M.; Tsutsumi, S.; Yamaguchi, Y.; Okabe, A.; Shimada, M.; Kawase, T.; Sada, A.; Ezawa, I.; Takano, Y.; et al. PHLDA1, another PHLDA family protein that inhibits Akt. Cancer Sci. 2018, 109, 3532–3542. [Google Scholar] [CrossRef] [PubMed]
  13. Sakthianandeswaren, A.; Christie, M.; D’Andreti, C.; Tsui, C.; Jorissen, R.N.; Li, S.; Fleming, N.I.; Gibbs, P.; Lipton, L.; Malaterre, J.; et al. PHLDA1 expression marks the putative epithelial stem cells and contributes to intestinal tumorigenesis. Cancer Res. 2011, 71, 3709–3719. [Google Scholar] [CrossRef]
  14. Nagai, M.A. Pleckstrin homology-like domain, family A, member 1 (PHLDA1) and cancer. Biomed. Rep. 2016, 4, 275–281. [Google Scholar] [CrossRef]
  15. Müller, S.; van den Boom, D.; Zirkel, D.; Köster, H.; Berthold, F.; Schwab, M.; Westphal, M.; Zumkeller, W. Retention of imprinting of the human apoptosis-related gene TSSC3 in human brain tumors. Hum. Mol. Genet. 2000, 9, 757–763. [Google Scholar] [CrossRef]
  16. McMinn, J.; Wei, M.; Schupf, N.; Cusmai, J.; Johnson, E.B.; Smith, A.C.; Weksberg, R.; Thaker, H.M.; Tycko, B. Unbalanced placental expression of imprinted genes in human intrauterine growth restriction. Placenta 2006, 27, 540–549. [Google Scholar] [CrossRef]
  17. Anderson, J.; Gordon, A.; McManus, A.; Shipley, J.; Pritchard-Jones, K. Disruption of imprinted genes at chromosome region 11p15.5 in paediatric rhabdomyosarcoma. Neoplasia 1999, 1, 340–348. [Google Scholar] [CrossRef]
  18. Xu, X.L.; Wu, L.C.; Du, F.; Davis, A.; Peyton, M.; Tomizawa, Y.; Maitra, A.; Tomlinson, G.; Gazdar, A.F.; Weissman, B.E.; et al. Inactivation of human SRBC, located within the 11p15.5-p15.4 tumor suppressor region, in breast and lung cancers. Cancer Res. 2001, 61, 7943–7949. [Google Scholar] [PubMed]
  19. Moon, H.G.; Oh, K.; Lee, J.; Lee, M.; Kim, J.Y.; Yoo, T.K.; Seo, M.W.; Park, A.K.; Ryu, H.S.; Jung, E.J.; et al. Prognostic and functional importance of the engraftment-associated genes in the patient-derived xenograft models of triple-negative breast cancers. Breast Cancer Res. Treat. 2015, 154, 13–22. [Google Scholar] [CrossRef]
  20. Kawase, T.; Ohki, R.; Shibata, T.; Tsutsumi, S.; Kamimura, N.; Inazawa, J.; Ohta, T.; Ichikawa, H.; Aburatani, H.; Tashiro, F.; et al. PH domain-only protein PHLDA3 is a p53-regulated repressor of Akt. Cell 2009, 136, 535–550. [Google Scholar] [CrossRef] [Green Version]
  21. Lei, L.; Wang, Y.; Li, Z.H.; Fei, L.R.; Huang, W.J.; Zheng, Y.W.; Liu, C.C.; Yang, M.Q.; Wang, Z.; Zou, Z.F.; et al. PHLDA3 promotes lung adenocarcinoma cell proliferation and invasion via activation of the Wnt signaling pathway. Lab. Investig. 2021, 101, 1130–1141. [Google Scholar] [CrossRef]
  22. Mangone, F.R.; Valoyes, M.A.; do Nascimento, R.G.; Conceição, M.P.F.; Bastos, D.R.; Pavanelli, A.C.; Soares, I.C.; de Mello, E.S.; Nonogaki, S.; de T Osório, C.A.B.; et al. Prognostic and predictive value of Pleckstrin homology-like domain, family A family members in breast cancer. Biomark. Med. 2020, 14, 1537–1552. [Google Scholar] [CrossRef] [PubMed]
  23. Indarte, M.; Puentes, R.; Maruggi, M.; Ihle, N.T.; Grandjean, G.; Scott, M.; Ahmed, Z.; Meuillet, E.J.; Zang, S.; Lemos, R.; et al. An Inhibitor of the Pleckstrin Homology Domain of CNK1 Selectively Blocks the Growth of Mutant KRAS Cells and Tumors. Cancer Res. 2019, 79, 3100–3111, Erratum in Cancer Res. 2019, 79, 5457. [Google Scholar] [CrossRef]
  24. Nowak, E.; Bednarek, I. Aspects of the Epigenetic Regulation of EMT Related to Cancer Metastasis. Cells 2021, 10, 3435. [Google Scholar] [CrossRef]
  25. Peng, X.; Yang, R.; Song, J.; Wang, X.; Dong, W. Calpain2 Upregulation Regulates EMT-Mediated Pancreatic Cancer Metastasis via the Wnt/β-Catenin Signaling Pathway. Front. Med. 2022, 9, 783592. [Google Scholar] [CrossRef]
  26. Zhao, Y.; Cai, C.; Zhang, M.; Shi, L.; Wang, J.; Zhang, H.; Ma, P.; Li, S. Ephrin-A2 promotes prostate cancer metastasis by enhancing angiogenesis and promoting EMT. J. Cancer Res. Clin. Oncol. 2021, 147, 2013–2023. [Google Scholar] [CrossRef]
  27. Kong, D.; Zhou, H.; Neelakantan, D.; Hughes, C.J.; Hsu, J.Y.; Srinivasan, R.R.; Lewis, M.T.; Ford, H.L. VEGF-C mediates tumor growth and metastasis through promoting EMT-epithelial breast cancer cell crosstalk. Oncogene 2021, 40, 964–979. [Google Scholar] [CrossRef]
  28. Rezatabar, S.; Karimian, A.; Rameshknia, V.; Parsian, H.; Majidinia, M.; Kopi, T.A.; Bishayee, A.; Sadeghinia, A.; Yousefi, M.; Monirialamdari, M.; et al. RAS/MAPK signaling functions in oxidative stress, DNA damage response and cancer progression. J. Cell Physiol 2019. online ahead of print. [Google Scholar] [CrossRef] [PubMed]
  29. Ledda, F.; Paratcha, G. Negative Regulation of Receptor Tyrosine Kinase (RTK) Signaling: A Developing Field. Biomark. Insights 2007, 2, 45–58. [Google Scholar] [CrossRef]
  30. Jham, B.C.; Ma, T.; Hu, J.; Chaisuparat, R.; Friedman, E.R.; Pandolfi, P.P.; Schneider, A.; Sodhi, A.; Montaner, S. Amplification of the angiogenic signal through the activation of the TSC/mTOR/HIF axis by the KSHV vGPCR in Kaposi’s sarcoma. PLoS ONE 2011, 6, e19103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Monti, P.; Menichini, P.; Speciale, A.; Cutrona, G.; Fais, F.; Taiana, E.; Neri, A.; Bomben, R.; Gentile, M.; Gattei, V.; et al. Heterogeneity of TP53 Mutations and P53 Protein Residual Function in Cancer: Does It Matter? Front. Oncol. 2020, 10, 593383. [Google Scholar] [CrossRef] [PubMed]
  32. Marbun, V.M.G.; Erlina, L.; Lalisang, T.J.M. Genomic landscape of pathogenic mutation of APC, KRAS, TP53, PIK3CA, and MLH1 in Indonesian colorectal cancer. PLoS ONE 2022, 17, e0267090. [Google Scholar] [CrossRef]
  33. Selvaraj, J.; Yasothkumar, D.; Vishnu Priya, V.; Raj, A.T.; Babu, S.D.; Patil, S. Development and tumorigenic potential of TP53: A therapeutic target for head and neck squamous cell carcinoma. Oral Oncol. 2022, 130, 105922. [Google Scholar] [CrossRef]
  34. Vokes, N.I.; Chambers, E.; Nguyen, T.; Coolidge, A.; Lydon, C.A.; Le, X.; Sholl, L.; Heymach, J.V.; Nishino, M.; Van Allen, E.M.; et al. Concurrent TP53 Mutations Facilitate Resistance Evolution in EGFR-Mutant Lung Adenocarcinoma. J. Thorac. Oncol. 2022, 17, 779–792. [Google Scholar] [CrossRef] [PubMed]
  35. Jones, P.A. Functions of DNA methylation: Islands, start sites, gene bodies and beyond. Nat. Rev. Genet. 2012, 13, 484–492. [Google Scholar] [CrossRef]
  36. Ceccarelli, V.; Racanicchi, S.; Martelli, M.P.; Nocentini, G.; Fettucciari, K.; Riccardi, C.; Marconi, P.; Di Nardo, P.; Grignani, F.; Binaglia, L.; et al. Eicosapentaenoic acid demethylates a single CpG that mediates expression of tumor suppressor CCAAT/enhancer-binding protein delta in U937 leukemia cells. J. Biol. Chem. 2011, 286, 27092–27102. [Google Scholar] [CrossRef]
  37. Barcena-Varela, M.; Colyn, L.; Fernandez-Barrena, M.G. Epigenetic Mechanisms in Hepatic Stellate Cell Activation During Liver Fibrosis and Carcinogenesis. Int. J. Mol. Sci. 2019, 20, 2507. [Google Scholar] [CrossRef] [PubMed]
  38. Liu, L.; Wang, C.; Li, S.; Qu, Y.; Xue, P.; Ma, Z.; Zhang, X.; Bai, H.; Wang, J. ERO1L Is a Novel and Potential Biomarker in Lung Adenocarcinoma and Shapes the Immune-Suppressive Tumor Microenvironment. Front. Immunol. 2021, 12, 677169. [Google Scholar] [CrossRef] [PubMed]
  39. Ribas, A.; Wolchok, J.D. Cancer immunotherapy using checkpoint blockade. Science 2018, 359, 1350–1355. [Google Scholar] [CrossRef]
  40. Bolandi, N.; Derakhshani, A.; Hemmat, N.; Baghbanzadeh, A.; Asadzadeh, Z.; Afrashteh Nour, M.; Brunetti, O.; Bernardini, R.; Silvestris, N.; Baradaran, B. The Positive and Negative Immunoregulatory Role of B7 Family: Promising Novel Targets in Gastric Cancer Treatment. Int. J. Mol. Sci. 2021, 22, 10719. [Google Scholar] [CrossRef]
  41. Zhang, R.; Ma, L.; Wei, Y.; Wei, K.; Song, T.; Du, Z.; Feng, Z. KIF22 Promotes Development of Pancreatic Cancer by Regulating the MEK/ERK/P21 Signaling Axis. Biomed. Res. Int. 2022, 2022, 6000925. [Google Scholar] [CrossRef] [PubMed]
  42. Jiang, H.; Courau, T.; Borison, J.; Ritchie, A.J.; Mayer, A.T.; Krummel, M.F.; Collisson, E.A. Activating Immune Recognition in Pancreatic Ductal Adenocarcinoma via Autophagy Inhibition, MEK Blockade, and CD40 Agonism. Gastroenterology 2022, 162, 590–603.e14, Erratum in Gastroenterology 2022, 162, 1785. [Google Scholar] [CrossRef]
  43. Kawaguchi, K.; Igarashi, K.; Miyake, K.; Lwin, T.M.; Miyake, M.; Kiyuna, T.; Hwang, H.K.; Murakami, T.; Delong, J.C.; Singh, S.R.; et al. MEK inhibitor trametinib in combination with gemcitabine regresses a patient-derived orthotopic xenograft (PDOX) pancreatic cancer nude mouse model. Tissue Cell 2018, 52, 124–128. [Google Scholar] [CrossRef] [PubMed]
  44. Chandrashekar, D.S.; Bashel, B.; Balasubramanya, S.A.H.; Creighton, C.J.; Ponce-Rodriguez, I.; Chakravarthi, B.; Varambally, S. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia 2017, 19, 649–658. [Google Scholar] [CrossRef]
  45. Chandrashekar, D.S.; Karthikeyan, S.K.; Korla, P.K.; Patel, H.; Shovon, A.R.; Athar, M.; Netto, G.J.; Qin, Z.S.; Kumar, S.; Manne, U.; et al. UALCAN: An update to the integrated cancer data analysis platform. Neoplasia 2022, 25, 18–27. [Google Scholar] [CrossRef]
  46. Uhlén, M.; Fagerberg, L.; Hallström, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, Å.; Kampf, C.; Sjöstedt, E.; Asplund, A.; et al. Proteomics. Tissue-based map of the human proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef]
  47. Lánczky, A.; Győrffy, B. Web-Based Survival Analysis Tool Tailored for Medical Research (KMplot): Development and Implementation. J. Med. Internet Res. 2021, 23, e27633. [Google Scholar] [CrossRef] [PubMed]
  48. Nagy, Á.; Munkácsy, G.; Győrffy, B. Pancancer survival analysis of cancer hallmark genes. Sci Rep. 2021, 11, 6047. [Google Scholar] [CrossRef]
  49. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci. Signal. 2013, 6, pl1. [Google Scholar] [CrossRef]
  50. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The cBio cancer genomics portal: An open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012, 2, 401–404, Erratum in Cancer Discov. 2012, 2, 960. [Google Scholar] [CrossRef] [Green Version]
  51. Men, C.; Chai, H.; Song, X.; Li, Y.; Du, H.; Ren, Q. Identification of DNA methylation associated gene signatures in endometrial cancer via integrated analysis of DNA methylation and gene expression systematically. J. Gynecol. Oncol. 2017, 28, e83. [Google Scholar] [CrossRef]
  52. Shinawi, T.; Hill, V.K.; Krex, D.; Schackert, G.; Gentle, D.; Morris, M.R.; Wei, W.; Cruickshank, G.; Maher, E.R.; Latif, F. DNA methylation profiles of long- and short-term glioblastoma survivors. Epigenetics 2013, 8, 149–156. [Google Scholar] [CrossRef] [PubMed]
  53. Ru, B.; Wong, C.N.; Tong, Y.; Zhong, J.Y.; Zhong, S.; Wu, W.C.; Chu, K.C.; Wong, C.Y.; Lau, C.Y.; Chen, I.; et al. TISIDB: An integrated repository portal for tumor-immune system interactions. Bioinformatics 2019, 35, 4200–4202. [Google Scholar] [CrossRef] [PubMed]
  54. Li, T.; Fan, J.; Wang, B.; Traugh, N.; Chen, Q.; Liu, J.S.; Li, B.; Liu, X.S. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 2017, 77, e108–e110. [Google Scholar] [CrossRef]
  55. Mermel, C.H.; Schumacher, S.E.; Hill, B.; Meyerson, M.L.; Beroukhim, R.; Getz, G. GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol. 2011, 12, R41. [Google Scholar] [CrossRef] [PubMed]
  56. Vasaikar, S.V.; Straub, P.; Wang, J.; Zhang, B. LinkedOmics: Analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018, 46, D956–D963. [Google Scholar] [CrossRef]
  57. Zhou, Y.; Zhou, B.; Pache, L.; Chang, M.; Khodabakhshi, A.H.; Tanaseichuk, O.; Benner, C.; Chanda, S.K. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 2019, 10, 1523. [Google Scholar] [CrossRef]
  58. Liu, C.J.; Hu, F.F.; Xia, M.X.; Han, L.; Zhang, Q.; Guo, A.Y. GSCALite: A web server for gene set cancer analysis. Bioinformatics 2018, 34, 3771–3772. [Google Scholar] [CrossRef] [PubMed]
  59. Von Mering, C.; Huynen, M.; Jaeggi, D.; Schmidt, S.; Bork, P.; Snel, B. STRING: A database of predicted functional associations between proteins. Nucleic Acids Res. 2003, 31, 258–261. [Google Scholar] [CrossRef]
  60. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  61. Doncheva, N.T.; Morris, J.H.; Gorodkin, J.; Jensen, L.J. Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J. Proteome Res. 2019, 18, 623–632. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 10316 g001 550
Figure 1. Expression analysis of the PHLDA family in pan-cancer and HPA analysis of PHLDA1/3 in PAAD. (A) The expression level of PHLDA1; (B) the expression level of PHLDA2; (C) the expression level of PHLDA3. ns, p ≥ 0.05; * p < 0.05; ** p < 0.01; *** p < 0.001. (D) The protein level of PHLDA1 (40×); (E) the protein level of PHLDA3 (40×).
Figure 1. Expression analysis of the PHLDA family in pan-cancer and HPA analysis of PHLDA1/3 in PAAD. (A) The expression level of PHLDA1; (B) the expression level of PHLDA2; (C) the expression level of PHLDA3. ns, p ≥ 0.05; * p < 0.05; ** p < 0.01; *** p < 0.001. (D) The protein level of PHLDA1 (40×); (E) the protein level of PHLDA3 (40×).
Ijms 23 10316 g001
Ijms 23 10316 g002 550
Figure 2. The relationship between expression levels of PHLDA family members and clinicopathologic features, TP53 mutation status, and prognosis of PAAD patients. (A) High expression of PHLDA2/3 was significantly correlated with a previous history of pancreatitis; (B) high expression of PHLDA1/2 was significantly correlated with a worse initial treatment effect; (C) the expression levels of PHLDA1/3 in the TP53 mutant PAAD tissues were significantly higher than those in the TP53 wild−type PAAD tissues; (D,E) the increased expression levels of PHLDA1 and PHLDA2 were significantly correlated with shorter OS and RFS. ** p < 0.01; *** p < 0.001.
Figure 2. The relationship between expression levels of PHLDA family members and clinicopathologic features, TP53 mutation status, and prognosis of PAAD patients. (A) High expression of PHLDA2/3 was significantly correlated with a previous history of pancreatitis; (B) high expression of PHLDA1/2 was significantly correlated with a worse initial treatment effect; (C) the expression levels of PHLDA1/3 in the TP53 mutant PAAD tissues were significantly higher than those in the TP53 wild−type PAAD tissues; (D,E) the increased expression levels of PHLDA1 and PHLDA2 were significantly correlated with shorter OS and RFS. ** p < 0.01; *** p < 0.001.
Ijms 23 10316 g002
Ijms 23 10316 g003 550
Figure 3. Genetic variation analysis and promoter methylation analysis of the PHLDA family in PAAD. (A) Gene variation characteristics of the PHLDA family in PAAD; (B,C) the mRNA high mutation of PHLDA2 was significantly correlated with PAAD invasion of surrounding tissues and tumor sites; (D) the amplification variation of PHLDA3 was significantly correlated with previous drinking history of PAAD patients; (E,F) PAAD patients with PHLDA1 mRNA high mutation had shorter OS and DFS; (G) the promoter methylation levels of PHLDA1/3 in the tumor tissues of patients were significantly higher than those in the paracancerous tissues; (H) the promoter methylation level of PHLDA2 was significantly lower in TP53 mutant PAAD tissues than that in TP53 wild-type PAAD tissues, but PHLDA3 was the opposite.
Figure 3. Genetic variation analysis and promoter methylation analysis of the PHLDA family in PAAD. (A) Gene variation characteristics of the PHLDA family in PAAD; (B,C) the mRNA high mutation of PHLDA2 was significantly correlated with PAAD invasion of surrounding tissues and tumor sites; (D) the amplification variation of PHLDA3 was significantly correlated with previous drinking history of PAAD patients; (E,F) PAAD patients with PHLDA1 mRNA high mutation had shorter OS and DFS; (G) the promoter methylation levels of PHLDA1/3 in the tumor tissues of patients were significantly higher than those in the paracancerous tissues; (H) the promoter methylation level of PHLDA2 was significantly lower in TP53 mutant PAAD tissues than that in TP53 wild-type PAAD tissues, but PHLDA3 was the opposite.
Ijms 23 10316 g003
Ijms 23 10316 g004 550
Figure 4. Analysis of the immune infiltration of the PHLDA family in PAAD. (A–E) PHLDA1 was positively correlated with 23 types of immune promoters, including IL6, NT5E, and TNFSF9, while PHLDA2 was negatively correlated with most immunomodulators (immune promoters, MHC molecules, chemokines, and chemokine receptors); (F–H) the copy number alteration of PHLDA1/3 may affect the infiltration levels of six types of dominant immune cells; (I) expression of PHLDA family was significantly different among the five immune subtypes. ns, p ≥ 0.05; * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 4. Analysis of the immune infiltration of the PHLDA family in PAAD. (A–E) PHLDA1 was positively correlated with 23 types of immune promoters, including IL6, NT5E, and TNFSF9, while PHLDA2 was negatively correlated with most immunomodulators (immune promoters, MHC molecules, chemokines, and chemokine receptors); (F–H) the copy number alteration of PHLDA1/3 may affect the infiltration levels of six types of dominant immune cells; (I) expression of PHLDA family was significantly different among the five immune subtypes. ns, p ≥ 0.05; * p < 0.05; ** p < 0.01; *** p < 0.001.
Ijms 23 10316 g004
Ijms 23 10316 g005 550
Figure 5. Gene enrichment analysis, PPI network construction, and drug sensitivity analysis of PHLDA Family in PAAD. (A) The GO enrichment of the BP terms of the PHLDA family and its 600 co-expressed genes; (B) the KEGG enrichment of the PHLDA family and its 600 co-expressed genes; (C) the PHLDA family played an activating role in a variety of oncogenic pathways; (D,E) TP53 played an important role in PPI network, which was closely related to the PHLDA family; (F) the expression of PHLDA family was negatively correlated with the sensitivity of various PAAD-targeting or chemotherapeutic drugs, including gemcitabine, docetaxel, and fluorouracil.
Figure 5. Gene enrichment analysis, PPI network construction, and drug sensitivity analysis of PHLDA Family in PAAD. (A) The GO enrichment of the BP terms of the PHLDA family and its 600 co-expressed genes; (B) the KEGG enrichment of the PHLDA family and its 600 co-expressed genes; (C) the PHLDA family played an activating role in a variety of oncogenic pathways; (D,E) TP53 played an important role in PPI network, which was closely related to the PHLDA family; (F) the expression of PHLDA family was negatively correlated with the sensitivity of various PAAD-targeting or chemotherapeutic drugs, including gemcitabine, docetaxel, and fluorouracil.
Ijms 23 10316 g005
Table
Table 1. Relationship between PHLDA family expression and the abundance of tumor-infiltrating immune cells in PAAD.
Table 1. Relationship between PHLDA family expression and the abundance of tumor-infiltrating immune cells in PAAD.
DescriptionPHLDA1 PHLDA2 PHLDA3
Purity Purity Purity
CorpCorpCorp
Act CD8−0.0730.334−0.363***−0.0130.858
Tcm CD80.1460.05140.0080.9110.0610.415
Tem CD80.0210.776−0.244**−0.0280.715
Act CD40.249***−0.0420.576−0.0870.247
Tcm CD40.341***0.209**0.298***
Tem CD4−0.1060.158−0.546***−0.313***
Tfh−0.050.509−0.294***−0.0770.305
Tgd0.0730.329−0.1410.05890.1350.0707
Th1−0.0010.993−0.251***−0.0240.755
Th170.0270.7220.363***−0.0210.782
Th20.238**−0.157*−0.174*
Treg0.0360.627−0.304***−0.0510.495
Act B−0.218**−0.363***−0.1420.0587
Imm B−0.0880.242−0.326***−0.0710.345
Mem B−0.0420.578−0.429***−0.198**
NK0.0440.557−0.272***−0.0410.587
CD56bright0.279***0.328***0.405***
CD56dim0.34***0.459***0.412***
MDSC0.0180.811−0.171*0.0340.655
NKT0.0510.499−0.243**−0.0870.248
Act DC0.199**0.1360.07010.178*
pDC0.257***−0.0530.4810.0490.51
iDC−0.130.084−0.192*−0.194***
Macrophage0.0060.932−0.251***0.0290.696
Eosinophil−0.1060.159−0.379***−0.253***
Mast0.0180.812−0.391***−0.1190.111
Monocyte0.0740.3250.163*0.0070.925
Neutrophil0.0320.668−0.0780.299−0.0970.196
Notes: * p < 0.05; ** p < 0.01; *** p < 0.001.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Duan, Y.; Du, Y.; Gu, Z.; Zheng, X.; Wang, C. Prognostic Value, Immune Signature, and Molecular Mechanisms of the PHLDA Family in Pancreatic Adenocarcinoma. Int. J. Mol. Sci. 2022, 23, 10316. https://doi.org/10.3390/ijms231810316

AMA Style

Duan Y, Du Y, Gu Z, Zheng X, Wang C. Prognostic Value, Immune Signature, and Molecular Mechanisms of the PHLDA Family in Pancreatic Adenocarcinoma. International Journal of Molecular Sciences. 2022; 23(18):10316. https://doi.org/10.3390/ijms231810316

Chicago/Turabian Style

Duan, Yunjie, Yongxing Du, Zongting Gu, Xiaohao Zheng, and Chengfeng Wang. 2022. "Prognostic Value, Immune Signature, and Molecular Mechanisms of the PHLDA Family in Pancreatic Adenocarcinoma" International Journal of Molecular Sciences 23, no. 18: 10316. https://doi.org/10.3390/ijms231810316

APA Style

Duan, Y., Du, Y., Gu, Z., Zheng, X., & Wang, C. (2022). Prognostic Value, Immune Signature, and Molecular Mechanisms of the PHLDA Family in Pancreatic Adenocarcinoma. International Journal of Molecular Sciences, 23(18), 10316. https://doi.org/10.3390/ijms231810316

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop