3.1 The expression of PSEN1 in different human tissues and its differential expression pattern in pan-cancer
Based on the HPA data, we analyzed the expression of PSEN1 in various tissues and organs of the human body (Fig. 1A). PSEN1 expression is relatively high in gastrointestinal tissues and organs, especially in the jejunum and small intestine. And its expression is relatively low in skeletal muscle and epididymis. Furthermore, PSEN1 presents a differential expression pattern in the vast majority of human cancer types (Fig. 1B). For example, PSEN1 is highly expressed in tumor tissues of multiforming glioma (GBM), endometrial carcinoma (UCEC), breast invasive carcinoma (BRCA), cervical squamous carcinoma and adenocarcinoma (CESC), lung adenocarcinoma (LUAD), esophageal carcinoma (ESCA), gastric and esophageal carcinoma (STES), mixed renal carcinoma (KIRAN), colorectal adenocarcinoma (COADERAD), prostate cancer (PRAD), gastric adenocarcinoma (STAD), head and neck squamous cell carcinoma (HNSC), clear cell carcinoma, renal clear cell carcinoma (KIRC), lung squamous carcinoma (LUSC), hepatocellular carcinoma (LIHC), bladder urothelial carcinoma (B LCA), thyroid cancer (THCA), uterine sarcoma (UCS), pheochromocytomas, and paragangliomas (PCPG). Protein expression data from the CPTAC database also showed that PSEN1 shows high expression status in tumor tissues of many cancer types, such as endometrial cancer, lung cancer, liver cancer, etc. (Fig. 1C).
Moreover, we also tried to further reveal the expression status of PSEN1 in the cells of each tissue from the single-cell level. We found that PSEN1 was significantly enriched in oligodendrocytes and monocytes (Figure S1A). From the SINGLE CELL TISSUES OVERVIEW (Figure S1B), we found the level of single-cell RNA enrichment of PSEN1 in each tissue/organ varies by tissue cell specificity. For example, the brain and ocular tissues mainly include neuronal cells and glial cells, and PSEN1 is predominantly enriched in oligodendrocytes. Lung tissue included specialized epithelial cells, blood & immune cells, endothelial cells, mesenchymal cells, glandular epithelial cells, and muscle cells, while PSEN1 was mainly enriched in macrophages in blood & immune cells (Figure S1C).
3.2 Genomic variation landscape and its subcellular localization of PSEN1 in human cancers
Data from cBioPortal demonstrate the frequency and type of genomic variation of the PSEN1 in human cancers. We found that the frequency of PSEN1 gene alteration in human pan-cancer was not high, at about 1.2%, and it was dominated by missense mutation, followed by CNV amplification and deep deletion (Figure S2A). We also further explored the specific proportion of each mutation type of the PSEN1 in each cancer type by SangerBox, as shown in Figure S2B. Most genomic alterations of PSEN1 are dominated by missense mutation. The amino acid range of PSEN 1 is 467AA. Furthermore, we explored the genetic alteration status of the PSEN1 in different types of cancer in the TCGA cohort. As shown in Figure S2C, PSEN1 was most frequently altered in patients with mature B cell tumors with “mutation” and “amplification” as the major changes, over 6%. However, in patients with NSCLC, genomic alterations in PSEN1 are dominated by “mutation”, “amplification” and “deep deletion”. Figure S2D also shows the 3 D structure of PSEN1 and its specific mutation types and sites in cancer. Overall, PSEN1 is genetically altered with low frequency in human cancers, and the types of its genomic alterations vary by cancer type. This genomic alteration did not affect the survival of cancer patients (Figure S3).
Next, we also investigated the expression of PSEN1 in various cell lines (Fig. 2A). PSEN1 showed relatively high expression in PC-3 and HaCaT as well as THP-1. Among them, THP-1 is a monocyte, which is consistent with the previous single-cell RNA analysis results. The function of genes is closely related to their protein subcellular localization. Subsequently, we also investigated the protein localization of PSEN1 in the three human cancer cells based on the immunofluorescence data from HPA (Fig. 2B). We found that PSEN1 was mainly localized to golgi apparatus, nucleoplasm, and cell Junctions. Figure 2C shows this situation more intuitively.
3.3 PSEN1 was highly expressed in LUAD tissues and associated with poor prognosis
In the previous part, we systematically investigated the expression pattern and gene variation landscape of PSEN1 in pan-cancer. Next, we focused on the exploration of the role of PSEN1 in a specific cancer (LUAD). Based on the TCGA paired and unpaired LUAD samples, we observed a higher expression of PSEN1 in the tumor tissues (Fig. 3A-B). We also observed similar conditions in the three GEO datasets (Fig. 3C-D). Figure S4 illustrates the protein expression levels based on PSEN1 in the CPTAC dataset. Figure S5 also validated the above results at the IHC level. Subsequently, to validate the above results, we assessed the levels of PSEN1 protein expression in LUAD samples and adjacent normal tissues samples from the Renmin Hospital of Wuhan University. IHC showed that PSEN1 was overexpressed in LUAD tumor tissues compared with adjacent normal lung tissue (Fig. 3F), which confirms our preliminary observations from external datasets. To further confirm this expression pattern, we randomly detected PSEN1 protein expression in LUAD tumor tissues and paired adjacent normal tissues by Western blot (Fig. 3G), and the results also indicated that PSEN1 showed high expression status in tumor tissues.
Given that the potential role of PSEN1 in LUAD has not been fully elucidated, we also performed a Kaplan-Meier survival analysis based on the information of 535 LUAD samples from the TCGA-LUAD dataset. We found that high expression of PSEN1 was associated with worse overall survival (Fig. 3H). The same phenomenon was observed in subsequent subgroup analyses (Figure S6). In addition, we also performed survival analysis validation based on the information of 399 LUAD samples from the GSE72094 dataset. The result, as we expected, was that patients expressing high PSEN1 had a shorter overall survival (Fig. 3I).
3.4 PSEN1 promotes the malignant phenotype of LUAD cells in vitro
We further explored the effect of PSEN1 on the malignant biological behavior of LUAD cells. We transfected cells with the siRNA-PSEN1 and PSEN1 overexpression plasmids. The qRT-PCR results showed that the RNA expression level of PSEN1 in the si-PSEN1-#1 and si-PSEN1-#2 groups was significantly lower than that in the si-NC group (Fig. 4A), while the RNA expression level of PSEN1 in the OE-PSEN1 group was significantly higher than that in the control group (Vector group) (Fig. 4B). Subsequent Western blot results also confirmed the success of the PSEN 1 knockdown / overexpression treatment. In the CCK-8 assay, the proliferation rates in the si-PSEN1-#1 and si-PSEN1-#2 groups were markedly reduced compared to the si-NC group (Fig. 4C). Conversely, the cell proliferation rate was significantly higher in the OE-PSEN1 group than that in the Vector group (Fig. 4D). In the EdU assay, the number of EdU-positive cells in the si-PSEN1-#1 and si-PSEN1-#2 groups was significantly lower than that in the si-NC group (Fig. 4E-F). In contrast, a substantial increase in the proportion of EdU-positive cells was observed in the OE-PSEN1 group compared to the Vector group (Fig. 4E-F). These results suggest that PSEN1 can promote the proliferation of LUAD cells in vitro.
We subsequently investigated the effect of PSEN1 on the metastasis of LUAD cells. In the transwell migration assay, a significant reduction in the number of cells that migrated was noted in the si-PSEN1-#1 and si-PSEN1-#2 groups when compared to the si-NC group, as shown in Fig. 5A. Conversely, the migration of cells in the OE-PSEN1 group was considerably higher than that of the Vector group, illustrated in Fig. 5B. This result indicates that PSEN1 can promote the cell migration of both A549 and PC9 in vitro.
In the Transwell invasion assay, it was noted that the count of cells that successfully invaded through the membrane was considerably reduced in the si-PSEN1-#1 and si-PSEN1-#2 groups compared to the si-NC group, as shown in Fig. 5C. In contrast, the count of invasive cells in the OE-PSEN1 group was significantly elevated compared to the Vector group, as indicated in Fig. 5D. This result indicates that PSEN 1 promoted the cell invasion of A549 and PC9 in vitro.
Epithelial-mesenchymal transition (EMT) plays an important role in the development of LUAD. Therefore, we investigated the association of PSEN1 expression with EMT in LUAD cells. The expression levels of EMT-related proteins in LUAD cells (A549 and PC9) were determined by Western blot for 48 h after siRNA-based treatment. The results showed that the expression level of E-cadherin in si-PSEN1-#1 and si-PSEN1-#2 groups was significantly higher than that in si-NC group, while the expression level of N-cadherin and Vimentin was significantly lower than that in si-NC group (Fig. 5E-F). This initially indicates that knockdown of PSEN1 can inhibit EMT progression in human LUAD cells. Subsequently, we assessed the fluorescence intensity of proteins associated with EMT across various cell groups using immunofluorescence. The results showed that the fluorescence intensity of E-cadherin in si-PSEN1-#1 and si-PSEN1-#2 groups was significantly higher than that of si-NC group, conversely, the fluorescence intensity for N-cadherin and Vimentin was notably reduced compared to the si-NC group (Fig. 5G). This further suggested that knockdown of PSEN1 could inhibit EMT progression in human LUAD cells in vitro.
3.5 PSEN1 promotes the malignant phenotype of LUAD cells associated with the Notch1/EGFR pathway
To further investigate the mechanism of PSEN1 in LUAD, we compared differential genes between the PSEN1 high and low expression groups. Among them, 12780 genes were up-regulated in the PSEN1 high expression group, while 8064 genes were down-regulated (Fig. 6A). We also used the LinkedOmics database to conduct a pilot study on the genes coexpressed with PSEN1 in LUAD from the TCGA dataset. Figure 6B shows the genes showing a positive and negative correlation of expression with PSEN1. Subsequently, we crossed the differentially expressed genes in LUAD and the PSEN1-related genes, and finally obtained 3251 key genes (Fig. 6C). The pathway enrichment analysis revealed that the implicated genes predominantly participated in pathways such as the Notch signaling, ErbB signaling, and Ras signaling pathways (Fig. 6D). Considering the relevance of these pathways to cancer development and PSEN1, we speculate that the oncogenic role played by PSEN1 might also be relevant in the aforementioned pathways.
As the central catalytic component of the γ-secretase, variations and expression levels of PSEN1 can obviously directly influence the function and operational dynamics of the γ-secretase. Through the String analysis, we searched for proteins with existing interactions with PSEN1. We found that besides the other five subunits of γ-secretase, the proteins that directly interact with PSEN1 include NOTCH1, APP, CDH 1, TRAF6, CTNNB1, and CTNND2 (Fig. 7A). Notch, serving as a key transmembrane receptor, is among the primary substrates for the enzymatic action of γ-secretase. The Notch signaling pathway is also widely proven to be involved in cancer development. The functional enrichment analysis in this study also indicated that the oncogenic role of PSEN1 in LUAD may be associated with the Notch signaling pathway (Fig. 6D). To further clarify the interaction between NOTCH1 and PSEN1, we performed molecular docking on this topic. The results suggest that there is a high confidence interaction between NOTCH1 and PSEN 1 (Fig. 7B). In Western Blot, the expression level of Notch1 (NICD) in the si-PSEN1-#1 and si-PSEN1-#2 groups was significantly lower than that in the si-NC group (Fig. 7C-D). These results suggest that Notch1 is a direct downstream regulator of PSEN1. We also observed that the expression level of EGFR in the si-PSEN1-#1 and si-PSEN1-#2 groups was significantly lower than that in the si-NC group (Fig. 7C-D). Furthermore, at the transcriptional level, a pronounced positive association was observed between the expression levels of EGFR and PSEN1 (R = 0.393, p < 0.001) (Fig. 7E). The close relationship between Notch1 and EGFR was also fully demonstrated in several previous studies (15, 16). Therefore, we propose that PSEN1 can exert an oncogenic effect through regulation of the Notch1/EGFR pathway in LUAD.
3.6 YY1 enhances the malignant behavior of LUAD cells through the regulation of PSEN1 expression
Given that the transcription factors (TFs) play important roles in cancer, we therefore went on to explore the transcription factors that may regulate PSEN1 expression. First, we obtained the location of PSEN1 in the genome through the UCSC database, which is located in Chr 14:73136507–73223691. Their corresponding promoter sequences were also obtained. In the JASPAR database, by entering this promoter sequence and setting relative profile score threshold to 99.9%, we finally obtained 156 putative sites for a total of 66 TFs. The association of these 66 TFs with PSEN1 in LUAD is shown in Figure S7A. At the same time, considering the transcription factor results of PSEN1 predicted by ALGGEN-PROMO database (Figure S7B), we finally locked the transcription factors regulating PSEN1 to YY1, and Figure S7C shows the motif of YY1 binding to PSEN1promoter. We further confirmed the binding of YY1 to the PSEN1 promoter region by ChIP-qPCR (Fig. 7F). Subsequently, we performed YY1 knockdown in A549 cells, and the results of qRT-PCR and Western blot showed that the expression of PSEN1 was decreased after YY1 knockdown (Fig. 7G-H). In addition, based on the mRNA correlation analysis of the TCGA database, we also identified a distinct correlation linking YY1 with PSEN1 (R = 0.694, p < 0.001) (Fig. 7I). The above results indicate that YY1 is an upstream transcription factor of PSEN1 and can bind the PSEN1 promoter region to regulate the expression of PSEN1.
Previous studies (17, 18) showed that YY1 plays a crucial role in lung cancer, being overexpressed in tumor tissues and promoting the malignant cell phenotype. We next explored whether PSEN1 participates in the function of YY1 in LUAD cells. We noted that the suppression of YY1 substantially hindered the proliferation, migration, and invasion of LUAD cells. Conversely, elevating PSEN1 levels notably counteracted these inhibitory effects caused by YY1 knockdown (Fig. 8A-D). Since knockdown of PSEN1 inhibited EMT progression, we subsequently explored whether this process is regulated by YY1. We observed that the knockdown of YY1 enhanced the expression of E-cadherin and downregulated the expression of mesenchymal markers N-cadherin and Vimentin, and the overexpression of PSEN1 could largely reverse the EMT progression in LUAD cells caused by YY1 knockdown (Fig. 8E-F). Overall, these results indicate that PSEN1 is an important target of YY1 to promote LUAD progression.
3.7 PSEN1 promotes LUAD proliferation and EMT in vivo
To further investigate the role of PSEN1 in vivo, we injected the stable cell line PC9 as well as A549 cells into BALB/c nude mice to observe subcutaneous tumorigenesis in mice. Our findings indicated that the dimensions and mass of tumors in the PSEN1 knockdown group (KD) were considerably smaller than those in the control group (NC) (Fig. 9A-D). In addition, we analyzed the expression of EMT markers in tumor samples by IHC, and the results showed that the reduction in PSEN1 expression notably upregulated E-cadherin levels and concurrently reduced the levels of mesenchymal markers such as N-cadherin and Vimentin (Fig. 9E-G). These results suggest that PSEN1 can promote LUAD cell proliferation and EMT progression in vivo.