Next Article in Journal
Clinical Application of the New Prostate Imaging for Recurrence Reporting (PI-RR) Score Proposed to Evaluate the Local Recurrence of Prostate Cancer after Radical Prostatectomy
Next Article in Special Issue
p21-Activated Kinase: Role in Gastrointestinal Cancer and Beyond
Previous Article in Journal
Identification of Biomarkers Associated with Liver Metastasis Progression from Colorectal Cancer Using Exosomal RNA Profiling
Previous Article in Special Issue
Identification of Essential Tumor-Infiltrating Immune Cells and Relevant Genes in Left-Sided and Right-Sided Colon Cancers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A LETM2-Regulated PI3K-Akt Signaling Axis Reveals a Prognostic and Therapeutic Target in Pancreatic Cancer

1
Guangdong Provincial Key Laboratory of Malignant Tumor Epigenetics and Gene Regulation, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou 510120, China
2
Department of Gastroenterology, Sun Yat-sen Memorial Hospital, Zhongshan School of Medicine, Sun Yat-sen University, Guangzhou 510120, China
3
Department of Ophthalmology, No.903 Hospital of PLA Joint Logistic Support Force, Hangzhou 310013, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2022, 14(19), 4722; https://doi.org/10.3390/cancers14194722
Submission received: 1 September 2022 / Revised: 19 September 2022 / Accepted: 23 September 2022 / Published: 28 September 2022
(This article belongs to the Special Issue Signaling Pathway in Gastrointestinal Cancer)

Abstract

:

Simple Summary

LEMT2 was a newly discovered protein-encoding gene with little cancer research and an unclear mechanism. This study aimed to illustrate LETM2 as the crucial oncogene for tumor progression in pancreatic ductal adenocarcinoma (PDAC). We analyzed the expression level and prognostic value of LETM2 in multiple cancers using pan-cancer analysis and found that the LETM2 expression was the most significantly related to the dismal prognosis of PDAC. Immunohistochemical analyses showed that high LETM2 expression was correlated with poor outcomes of PDAC. In in vitro and in vivo experiments, LETM2 knockdown significantly inhibited tumor proliferation and metastasis, while LETM2 overexpression exerted the opposite effects. Then, we suggested that LETM2 may facilitate tumor progression by activating downstream PI3K-Akt signaling pathway in PDAC. In conclusion, the study enhanced our understanding of the LETM2-regulated PI3K-Akt signaling axis served as a prognostic and therapeutic target of pancreatic cancer.

Abstract

Pancreatic ductal adenocarcinoma (PDAC) is one of the highest mortalities malignant tumors, which is characterized by difficult diagnosis, rapid progression and high recurrence rate. Nevertheless, PDAC responds poorly to conventional therapies, which highlights the urgency to identify novel prognostic and therapeutic targets. LEMT2 was a newly discovered protein-encoding gene with little cancer research and an unclear mechanism. Thus, this study aimed to illustrate LETM2 as the crucial oncogene for tumor progression in PDAC. In this study, we analyzed the expression level and prognostic value of LETM2 in multiple cancers using pan-cancer analysis. The analyses based on the TCGA-GTEx dataset indicated that the LETM2 expression was obviously elevated in several cancers, and it was the most significantly related to the dismal prognosis of PDAC. Subsequently, we demonstrated the functional role and mechanism of LETM2 by clinical sample evaluation, and in in vitro and in vivo experiments. Immunohistochemical analyses showed that high expression of LETM2 was correlated with poor outcomes of PDAC. Moreover, we demonstrated that LETM2 knockdown significantly inhibited tumor proliferation and metastasis, and promoted cell apoptosis, while LETM2 overexpression exerted the opposite effects. Finally, the impairment caused by LETM2-knockdown could be recovered via excitation of the PI3k-Akt pathway in vitro and in vivo animal models, which suggested that LETM2 could activate the downstream PI3K-Akt pathway to participate in PDAC progression. In conclusion, the study enhanced our understanding of LETM2 as an oncogene hallmark of PDAC. LETM2 may facilitate tumor progression by activating the PI3K-Akt signaling pathway, which provides potential targets for the diagnosis, treatment, and prognosis of pancreatic cancer.

1. Introduction

Pancreatic ductal adenocarcinoma (PDAC), emanating from the exocrine compartment, represents over 90% of all pancreatic malignancies [1]. Being among the most aggressive and lethal malignant solid tumors, PDAC is distinguished by rapid progression and high mortality, as evidenced by it having the lowest overall 5-year survival rate of 8% among all types of cancer [2,3,4]. PDAC bears an extremely dismal prognosis mainly because of its difficult early detection and suboptimal efficacy of traditional therapies [5]. Therefore, the exploration of the crucial molecular pathways and functional mechanisms driving the development as well as the progression of PDAC is critically required to discover new diagnostic biomarkers and identify novel therapeutic strategies in clinical practice.
LETM2 (leucine zipper-EF-hand-containing transmembrane protein 2, LETM2) belongs to the LETM1 protein family, which comprises LETM2, localized on chromosome 8p11.2, and its important paralog LETM1, encoded on chromosome 4p16.3 [6,7,8]. LETM1 and LETM2 are initially identified as nucleus-encoded, mitochondrial inner-membrane proteins. Mitochondrial activity is extremely important for maintaining all sorts of normal physiological functions, the malfunctioning of which may induce dysregulation of biological function or occurrence of various disorders such as malignant tumors [9,10].
It was discovered that elevated expression levels of LETM1 were related to poor prognosis in multiple malignancies [11,12]. Furthermore, LETM1 performs a vital role in driving development and metastasis by PI3K-Akt signaling pathway activation in ovarian cancer, prostate cancer, gastric cancer, and so on [13,14,15]. However, studies involving the biological roles and functional mechanisms of LETM2 in malignant tumors are scarce. LETM2 is located in close proximity to the Wolf–Hirschhorn syndrome candidate-1-like gene-1 (WHSC1L1), which is associated with Wolf–Hirschhorn syndrome, a rare genetic condition characterized by intellectual disability, delayed growth, psychomotor retardation and seizures [16]. It has been suggested that there is LETM2 amplification in lung cancer [17]. Cheng et al. detected the high-level amplification of LETM2 in esophageal squamous cell carcinoma by whole-genome sequencing, and downregulation of LETM2 expression inhibited tumor proliferation in vitro [18]. To date, there are no relevant studies concerning the LETM2 expression and its functional roles in PDAC.
Regarded as one of the most fundamental intracellular signal transduction pathways, the PI3K-Akt pathway is considered the master regulator of cell growth, apoptosis, migration, metabolism, and so on [19,20]. PI3K (phosphoinositide 3-kinase) is an upstream effector of the PI3K/Akt signaling pathway. By promoting the phosphorylation at Ser473 as well as Thr308, PI3K may induce the activation of critical downstream effector Akt. This pathway is triggered under physiological conditions when it is responsive to insulin, cytokines, and growth hormones to govern important metabolic activities. On the contrary, aberrant activation of PI3K-Akt and its downstream pathway contributes to the tumorigenesis and metastasis of malignant tumors [21,22]. PDAC is closely connected with Akt activation, which has been detected in 60% of PDAC cases [23,24]. Moreover, previous studies demonstrated a tough correlation between high PI3K-Akt pathway activity and poor prognosis of PDAC patients [25]. In recent years, inhibitors targeting the PI3K-Akt signaling pathway rapidly became under focus for clinical investigation in anti-cancer therapy, and encouraging evidence indicates the potential of PI3K-Akt drugs in both hematological malignancies as well as solid tumors [26,27,28].
In this study, we screened critical genes and pathways for pancreatic cancer progression using pan-cancer analysis and bioinformatics techniques. Then, we observed differential expression of LETM2 in different malignancies in the TCGA-GTEx pan-cancer dataset, which was most significantly associated with a dismal prognosis of PDAC. Subsequently, immunohistochemical analyses showed that the obviously elevated expression of LETM2 was correlated with abysmal outcomes of pancreatic cancer patients. Moreover, we demonstrated that the knockdown of LETM2 significantly promoted cancer cell apoptosis and inhibited tumor migration and proliferation in both in vivo and in vitro experiments. Finally, we suggested that LETM2 could activate the downstream PI3K-Akt pathway to participate in PDAC development. In conclusion, the study proposed that LETM2 may facilitate tumor development by promoting the PI3K-Akt signaling pathway in PDAC, which provides prospective targets for the diagnosis, treatment, and prognosis of PDAC.

2. Materials and Methods

2.1. Evaluation of LETM2 Expression and Prognostic Significance in the Public Database

The TCGA pan-cancer dataset and the GTEx pan-cancer dataset containing normalized data were contacted and required from the UCSC Xena database [29] to conduct the pan-cancer analysis. Additionally, we downloaded RNA-sequencing expression profiles and relevant clinical information for 178 samples of pancreatic cancer from the TCGA dataset. The abbreviations and full names of cancers mentioned are listed in Supplementary Materials Table S1. The limma package in R software was applied to make a comparison of the LETM2 expression between tumor and normal tissues of different malignancies, which may be validated by the GEPIA database [30]. We divided all patients into high and low groups on the basis of the median value of LETM2 expression in each cancer type to examine the impact of LETM2 expression on the overall survival (OS) as well as disease-free survival (DFS) in the survival analysis of PDAC and the pan-cancer analysis. Moreover, the GSVA package and ggstatsplot package in R software was used to present the correlations between LETM2 expression and signaling pathway as well as immune infiltration level. All the analysis methods and R packages were implemented by R version 4.0.3.

2.2. Tissue Specimens and Cell Culture

The paraffin-embedded PDAC specimens and matched non-tumor (paracancerous) tissue from 60 PDAC patients following the operation were collected from Sun Yat-sen Memorial Hospital between January 2016 and December 2021. Three kinds of PDAC cell lines (MIA PaCa-2, SW1990, and BxPc-3) and immortalized human pancreatic duct cells (hTERT-HPNE) were acquired from the American Type Culture Collection (ATCC, Manassas, VA, USA). The cells were cultivated in accordance with the manufacturers’ instructions.

2.3. Cell Counting Kit-8 (CCK-8) Assays

CCK8 assays were performed to evaluate cell proliferation in vitro using the CCK8 Kit (Beyotime Biotechnology, Shanghai, China) on the basis of the manufacturer’s manual. In brief, cells were seeded in 96-well plates with 2 × 103 cells per well. Then, 10 μL CCK8 solution was added to each well for incubation at 37 °C for 2 h. Finally, the absorbance was measured at 450 nm by a microplate reader.

2.4. Transwell Assays

With or without Matrigel (BD Biosciences, Franklin Lakes, NJ, USA), the PDAC cells were injected in the top transwell chamber (Corning, Corning, NY, USA) and incubated for 24 h or 48 h. Thereafter, cotton swabs were used to remove unmigrated and uninvaded PDAC cells. After being preserved in methanol for 10 to 15 min at room temperature, the PDAC cells should be stained with 0.5% crystal violet. Three triplicates of the tests were performed.

2.5. Quantitative Real-Time PCR (qRT–PCR)

Following the manufacturer’s directions, total RNA was extracted from tissues and cells using TRIzol Reagent (Invitrogen, Waltham, MA, USA). Using the PrimeScript RT reagent Kit and SYBR Premix Ex Taq (both from Takara, Japan), qRT–PCR was performed with the primer sequences shown in Supplementary Materials Table S2. At least three times each of the trials were conducted.

2.6. Western Blotting (WB)

After extraction of total protein by using RIPA lysate (Solarbio, Beijing, China), the BCA Protein Assay Kit (Thermo Fisher Scientific, Waltham, MA, USA) was utilized to measure the protein concentrations. Equal amounts of proteins (20 μg) were separated on a 10% SDS/PAGE gel and electrotransferred (100 V, 2 h) to nitrocellulose membranes. The membranes were blocked with 5% BSA for 1 h and were then incubated with primary antibody for LETM2 (Proteintech, Wuhan, China) and GAPDH (Abcam, Cambridge, UK) overnight at 4 °C. The corresponding secondary antibodies were applied on the following day. The blots containing target bands were exposed by enhanced chemiluminescence (ECL) reagent (Solarbio, China) on the Exposure meter. The densitometry readings of each band were detected by ImageJ software 1.8.0 in the gray value analysis, and the relative expression was calculated as intensity ratio = target protein gray density/GAPDH gray density

2.7. Plasmid Construction

Small interfering RNA (siRNA) of LETM2 (siLETM2) and small hairpin RNA (shRNA) targeting LETM2 (shLETM2) were designed as shown in Supplementary Materials Table S3 and both synthesized by IGE Biotech (Guangzhou, China). To generate stable LETM2-knockdown cell lines, the PDAC cells were infected with lentiviral via transduction and then screened under 8 μg/mL puromycin pressure. Full-length LETM2 cDNA was designed and synthesized by IGE Biotech (Guangzhou, China), which was cloned into pcDNA3.1.

2.8. Immunohistochemistry (IHC)

Tissue samples from PDAC patients as previously mentioned were fixed in formalin, embedded in paraffin, and cut into 4 μm sections. The paraffin-embedded sections were deparaffinized, hydrated, and blocked. Additionally, then the sections were incubated with primary antibodies overnight, followed by secondary antibodies for 2 h. Finally, DAB+ chromogen was used for color development and hematoxylin was applied to counterstain. IHC score for each case was calculated by the following formula: IHC score = positive percentage score × intensity score. Staining percentage score was graded as follows: 0; 0–25% of positively staining cells = 1; 25–50% of positively staining cells = 2; 50–75% of positively staining cells = 3; >75% of positively staining cells = 4. Staining intensity score was graded as follows: negative = 0; low positive = 1; positive = 2; high positive = 3. LETM2 expression was qualified as low (IHC score 0–6) and high (IHC score 8–12).

2.9. Animal Experiments

The Animal Research Committee of the Sun Yat-sen University Cancer Center gave its approval to all animal trials. Female Balb/c nude mice (Balb/c-nu, weighing ~15–20 g), were purchased from the Guangdong Medical Laboratory Animal Center. A total of 20 mice were maintained under controlled environmental conditions in a 12 h light/dark cycle with unrestricted access to sterile feed and filtered water. We injected 5 × 106 transfected BxPc-3 cells suspended in 0.2 ml of PBS into the flank of nude mice to construct the subcutaneous xenograft tumor model (five mice per group). Tumor size was surveyed once in 5 days (V = Length × Width2 × 0.5). The humane endpoints were defined as follows: (i) maximum tumor size > 1.5 cm3; (ii) weight loss >20%; (iii) tumor ulceration or necrosis. After four weeks, all of the nude mice in the study were sacrificed. The organs and tumors from the mice were dissected, photographed, weighed, and stained.

2.10. Statistical Analysis

The software involved in this study to perform statistical analysis contained SPSS 22.0 software, Stata 15.0 software, and R software 4.0.3. Student’s t-test was adopted for statistical analysis in the comparison of normally distributed continuous variables, while the Mann–Whitney U test was applied when our statistical data were not normally distributed. Additionally, categorical variables were processed by use of chi-square or Fisher exact test. Furthermore, we took advantage of Cox proportional hazards regression model for both univariate and multivariate analyses to calculate the hazard ratios (HR) along with associated 95% confidence intervals (CI). Variables with a p-value < 0.1 detected in univariate analyses were eligible for inclusion in the multivariate analyses. A value of p < 0.05 was supposed as an indication of a statistically significant difference.

3. Results

3.1. Pan-Cancer Analysis and Prognostic Value of LETM2

To assess the expression level and prognostic values of LETM2 in various cancers, pan-cancer analysis was implemented ground on the TCGA database and GTEx database. As indicated in Figure 1A,B, LETM2 was significantly elevated in 12 kinds of cancers compared to both paired adjacent tissues in TCGA and normal tissues in GTEx, including lung adenocarcinoma (LUAD, p = 5.7 × 10−4), colon adenocarcinoma (COAD, p = 7.0 × 10−15), colon adenocarcinoma/rectum adenocarcinoma (COADREAD, p = 1.4 × 10−15), breast invasive carcinoma (BRCA, p = 5.2 × 10−13), esophageal carcinoma (ESCA, p = 5.8 × 10−3), stomach and esophageal carcinoma (STES, p = 5.3 × 10−6), kidney renal papillary cell carcinoma (KIRP, p = 3.0 × 10−4), stomach adenocarcinoma(STAD, p = 8.4 × 10−4), head and neck squamous cell carcinoma (HNSC, p = 3.2 × 10−10), head and neck squamous cell carcinoma (LIHC, p = 3.8 × 10−7), rectum adenocarcinoma (READ, p = 0.03), pancreatic adenocarcinoma (PDAC, p = 0.02), pheochromocytoma/paraganglioma (PCPG, p = 3.4 × 10−3), cholangiocarcinoma (CHOL, p = 6.4 × 10−6), and cervical squamous cell carcinoma/endocervical adenocarcinoma (CESC, p = 0.05). Detail data of the statistical analysis are described in Supplementary Materials Tables S4 and S5. Then, we used the GEPIA database to further validate our results, which indeed revealed a significant increase in LETM2 expression in PDAC (Figure 1C). For exact values, see Supplementary Materials Table S6. Moreover, pan-cancer survival analysis on the ground of the TCGA database was implemented to evaluate the predictive value and clinical significance of LETM2. The analysis in Figure 1D represented that elevated expression of LETM2 was most significantly associated with dismal prognosis of pancreatic cancer, either poor OS or DFS.

3.2. Prognostic Value and Clinical Significance of LETM2 for PDAC in TCGA Database

Therefore, we hereafter focused on the expression level and prognostic values of LETM2 in pancreatic cancer. The impact of LETM2 on OS and DFS of PDAC patients in the TCGA database was investigated. As presented in Figure 2A–C, the outcome showed that LETM2 overexpression negatively correlated with clinical outcome as OS (HR = 2.1, p = 0.00046) and DFS (HR = 1.7, p = 0.016) in PDAC patients. Subsequently, multivariate Cox regression analysis in TCGA database showed that LETM2 expression was an independent risk factor for OS of PDAC (Table 1). Additionally, we observed that LETM2 expression correlated with tumor inflammation signature (R = -0.15, p = 0.047), PI3K-AKT-mTOR pathway (R = −0.21, p = 0.006), and P53 pathway (R = 0.25, p = 0.001) (Figure 2D–F, Figure S1). On account of the relation between immune infiltration level and prognosis of malignancies, we further evaluated the correlation between LETM2 expression and immune infiltration level. As shown in Figure 2G,H and Figure S2, LETM2 expression correlated with immune infiltration level of B cell expression (R = −0.34, p = 3.95 × 10−6), and T cell CD8+ expression (R = 0.15, p = 0.045). The Sankey diagram in Figure 2I demonstrated the association between LETM2 expression and the histological grade of PDAC. The above findings suggested that LETM2 could bear the responsibility of a novel prognostic target for PDAC.

3.3. Elevated LETM2 Expression Correlates with Clinicopathological Features and Dismal Prognosis in PDAC Patients

To clarify the prognostic significance of LETM2 in PDAC patients, we first examined the LETM2 protein expression in six paired samples of PDAC tumor and non-tumor (paracarcinoma) tissues. The results showed that the LETM2 protein expression was remarkably increased in PDAC tumor tissues (Figure 3A,B). The densitometry readings/intensity ratio of each band and the original whole blot could be found in the Supplemental Materials Table S7 and Figure S3. The corresponding gene expression results of qRT–PCR also corroborated the WB findings, as indicated in Figure 3C. Moreover, IHC staining in specimens of 60 PDAC patients enrolled in our study indicated that LETM2 expression was significantly elevated in tumor tissue in comparison with the matched non-tumor tissue (n = 60, p < 0.001; Figure 3D,E). Then, we revealed significant positive correlations between IHC score and clinicopathological characteristics such as T stage, distant metastasis and pathologic stage, while there was no statistical significance between IHC score and lymph node metastasis (Figure 3F–H, Table S8). Additionally, Kaplan–Meier survival analysis was put into practice to investigate the prognostic significance of LETM2 in PDAC. The results as showcased in Figure 3I,J found that elevated LETM2 expression was in connection with worse OS and DFS. In brief, these findings implied elevated LETM2 expression was significantly linked to poor prognosis in PDAC.

3.4. LETM2 Drives PDAC Tumorigenesis and Metastasis In Vitro

We evaluated the LETM2 expression in PDAC cell lines to investigate the functional roles of LETM2 in PDAC. In contrast to human pancreatic duct cells (hTERT-HPNE), the findings showed a considerable elevation of LETM2 mRNA as well as protein levels in PDAC cell lines (MIA PaCa-2, SW1990, and BxPc-3) (Figure 4A–C). In addition, consistent data was obtained in the analysis of the CCLE (Cancer Cell Line Encyclopedia) database (Figure S4).
Therefore, we selected MIA PaCa-2 and SW1990 cells for subsequent experiments of LETM2 overexpression, and BxPc-3 and SW1990 cells for LETM2 knockdown experimental validation. We designed the overexpression lentivirus to construct overexpression cell lines named oeLETM2, while the names of the matching control cell lines were Ctrl. Moreover, BxPc-3 and SW1990 cell lines were transfected with siRNA of negative control (siNC) or targeting LETM2 (siLETM2-1, siLETM2-2 and siLETM2-3). The LETM2 expression of treated cells was confirmed with qRT–PCR and Western blotting. The findings suggested that the expression of LETM2 protein and mRNA was markedly increased in overexpression cell lines (Figure 4D–F), while both LETM2 protein and mRNA expression were significantly downregulated in knockdown cells (Figure 5A–C).
Firstly, we performed cell apoptosis analyses by flow cytometry. The results revealed that upregulation of LETM2 might lightly decrease cell apoptosis as demonstrated in Figure 4G, while the percentage of apoptotic cells was significantly higher in LETM2-knockdown PDAC cells (Figure 5D). Sequentially, we measured the growth of the above cells with LETM2 overexpression or knockdown by use of cell counting kit-8 (CCK-8) assays. Overexpression of LETM2 significantly promoted cell proliferation, and conversely (Figure 4I), LETM2 knockdown markedly led to a reduced trend in cell proliferation (Figure 5E). In addition, transwell assays were conducted to assess the migratory and invasive ability. As shown in Figure 4H, upregulation of LETM2 displayed a facilitation effect on not only migration but also invasion of PDAC cells. Additionally, the migratory and invasive ability of PDAC cells was prominently impeded after LETM2 knockdown (Figure 5F).

3.5. LETM2 Activates Downstream PI3K-Akt Signaling Axis in PDAC

To identify the target genes and key downstream pathways of LETM2, we performed a gene correlation analysis of LETM2 in the GEPIA database. Furthermore, the KEGG pathway and GO functional enrichment analyses were implemented on the top 100 co-expressed genes. According to the KEGG pathway enrichment analysis, several significantly enriched pathways (FDR < 0.1 and p < 0.05) were identified, including PI3K-Akt signaling pathway, axon guidance, small cell lung cancer, NF-kappa B signaling pathway, and so on (Figure 6A). Meanwhile, the molecular functions (MF) of GO functional enrichment analyses indicated signaling receptor binding, cadherin binding, extracellular matrix structural constituent, chemorepellent activity, etc. (Figure 6B). The biological processes (BP) of GO functional enrichment analyses indicated response to external stimulus, locomotion, cell adhesion, biological adhesion, etc. (Figure 6C). The cellular components (CC) of GO functional enrichment analyses indicated extracellular region, extracellular region part, vesicle, extracellular region space, cell junction, etc. (Figure 6D).
Based on our experimental analysis and previous reports in the literature, we thought the PI3K-Akt signaling pathway served as the critical downstream pathway of LETM2 in PDAC. Then, Western blotting was conducted to validate the regulatory relationship between LETM2 and the PI3K-Akt pathway. As showcased in Figure 6E,F, the key proteins of the PI3K-Akt pathway including p-PI3K and p-Akt significantly downregulated in LETM2-knockdown cells, while there observed no discernible differences in PI3K and Akt. These results implied that LETM2 may activate downstream PI3K-Akt signaling pathway through enhancing the level of PI3K phosphorylation as well as Akt phosphorylation.

3.6. LETM2 Accelerates PDAC Malignant Progression by Activating the PI3K-Akt Signaling Pathway

In order to further demonstrate the biological functions and molecular mechanisms of LETM2-regulated PI3K-Akt signaling axis in PDAC, we used the PI3K-Akt pathway inhibitor and activator to conduct rescue experiments in vitro and in vivo. As indicated in Figure 7A, the flow cytometry results implied that inhibiting the PI3K-Akt pathway can also increase PDAC cell apoptosis, which was consistent with the trend of LETM2 downregulation. Additionally, meanwhile, the recovery of PI3K-Akt pathway activities via activator would partially rescue the cell apoptosis after LETM2 knockdown. Similarly, the PI3K-Akt pathway activator can abolish or decrease the effect of LETM2 knockdown including inhibiting PDAC cell proliferation, migration, and invasion during the rescue experiments (Figure 7B,C).
Consistent with the above observations, it was evident that the downregulation of LETM2 and inhibition of the PI3k-Akt pathway had a negative influence on tumor proliferation in the PDAC subcutaneous tumor model of nude mice. The temporal scheme of in vivo experiments was shown in Figure 8A. However, the impairment caused by LETM2-knockdown could be recovered via excitation of the PI3k-Akt pathway (Figure 8B–D). Furthermore, IHC staining of the nude mice tissues (Figure 8E) also verified that the Ki67 was downregulated in the LETM2 downregulation group and it was also reversed in the tumor tissues with the PI3k-Akt pathway activator stimulation.
In conclusion, the findings of the present investigation confirmed that LETM2 accelerates PDAC malignant progression via activating the PI3K-Akt pathway.

4. Discussion

With a high mortality rate and rapid progression, PDAC remains one of the most treatment-refractory malignancies [4]. PDAC presents an extremely poor prognosis, and the incidence has been rising year by year. The 5-year OS rate of PDAC inferior to 3-10% is dismal in various studies [31,32,33]. However, the progress of diagnosis and treatment in PDAC is slow on account of its complexity and difficulty [34,35]. Meanwhile, PDAC treatment is still a significant challenge due to the poor efficacy of the current therapies. Surgical resection continues to be the sole means for a cure in patients with PDAC to date, but the majority are diagnosed too late for surgery due to the lack of specific diagnostic indicators. Many patients still experience relapse even after surgery and chemotherapy [36]. Therefore, the mechanism of pancreatic cancer should be further explored and it is of great significance to discover effective diagnostic and therapeutic targets.
In the current research, we firstly identified LETM2 as an oncogenic hallmark of PDAC by pan-cancer analysis of the TCGA + GTEx database and the results of GEPIA analysis also supported this finding. Pan-cancer survival analysis ground on the TCGA database was implemented to evaluate the predictive value and clinical significance of LETM2, which represented that elevated expression of LETM2 was most significantly associated with dismal prognosis of pancreatic cancer, either poor OS or DFS. Moreover, we hereafter focused on the prognostic values of LETM2 in pancreatic cancer. PDAC patients with elevated LETM2 expression had substantially shorter OS times than those with low LETM2 expression. Additionally, the ROC curve of LETM2 indicated a high predictive power. Furthermore, IHC staining in specimens of 60 PDAC patients enrolled in our study verified that LETM2 expression was significantly elevated in tumor tissue in comparison with the matched non-tumor tissue. Additionally, there revealed significant positive correlations between IHC score and clinicopathological characteristics such as stage and distant metastasis. Additionally, elevated LETM2 expression was significantly linked to poor prognosis in PDAC. The above findings suggested that LETM2 could bear the responsibility of a novel prognostic target for PDAC.
LETM2 was first found within a region on chromosome 8p11.2. LETM2 is a kind of protein-coding gene exhibiting high sequence similarity to LETM1. The 8p11.2 region comprises a cluster of genes containing LETM2, FGFR1, and WHSC1L1 [16]. Rearrangements of chromosome 8p are continually reported in various types of cancers, which were revealed to be tightly connected with malignant tumor development and progression [37,38]. There have been many reports about the close relationship between FGFR1 and cancer development [39,40,41]. In addition, several studies focused on the oncogenic role of WHSC1L1 in breast cancer, pancreatic adenocarcinoma, and lung tumor in recent years [42,43]. However, although there is differential expression of LETM2 in multiple cancers, research on the function and mechanism of LETM2 in tumor development and progression is rare. The prognostic value and clinical significance of LETM2 need to be further explored.
On the basis of this perspective, we next demonstrated the functional role and mechanism of LETM2 by in vitro and in vivo trials. In comparison to normal pancreatic duct cells, LETM2 mRNA and protein expression were observed to be markedly elevated in PDAC cells. LETM2 knockdown could not only suppress PDAC cell proliferation, migration, and invasion, but also promote PDAC cell apoptosis in vitro, while LETM2 overexpression may speed up tumor development and progression. Additionally, the tumorigenicity of PDAC cells was decreased by the knockdown of the LETM2 in vivo xenograft tumor model of nude mice. The above-mentioned results presented the important role of LETM2 indicating the prognosis and promoting tumor progression.
To identify the target genes and key downstream pathways of LETM2, we performed gene correlation analysis of LETM2 in the GEPIA database. Furthermore, the KEGG pathway and GO functional analyses were implemented on the top 100 co-expressed genes. According to the KEGG analysis, we considered the PI3K-Akt pathway as the critical downstream pathway of LETM2 in PDAC. Then, Western blotting was conducted to validate the regulatory relationship between LETM2 and the PI3K-Akt pathway. The key proteins of the PI3K-Akt pathway including p-PI3K and p-Akt significantly downregulated in LETM2-knockdown cells, while no discernible differences in PI3K and Akt were observed. These results implied that LETM2 may activate the downstream PI3K-Akt signaling pathway through enhancing the level of PI3K phosphorylation as well as Akt phosphorylation.
Dysregulation of the PI3K-Akt signaling pathway appears as one of the most frequent oncogenic events related to tumor progression in all kinds of malignancies [44,45]. The main downstream effectors of PI3K are AKT and mTOR (mammalian targets of rapamycin), and the pathway bears a fundamental oncogenic role in driving tumor evolution and progression in pancreatic cancer. Inhibitors aimed at the pathway could provide novel and important strategies for the treatment of pancreatic cancer [28,46].
Therefore, we used PI3K inhibitor GDC-0941 [47] and AKT activator SC-79 [48] to conduct rescue experiments in vitro as well as in vivo, in which we noticed that the PI3K-Akt pathway activator can abolish or decrease the effect of inhibiting tumor progression after LETM2 knockdown. Consistent with the above observations, it was evident that the downregulation of LETM2 and inhibition of the PI3k-Akt pathway had a negative influence on tumorigenesis in the PDAC subcutaneous xenograft tumor model of nude mice. The impairment caused by LETM2-knockdown could be recovered via activation of the PI3k-Akt pathway.
However, we noticed that knockdown LETM2 can suppress PDAC malignant progression more efficiently than Akt inhibitors. Additionally, after the recovery of the PI3K-Akt pathway by the activator, both the proliferation and metastasis ability of PDAC cells were incapable to return to the original level completely, which dropped a hint that LETM2 may adjust tumor progression by means of diverse signaling pathways besides the PI3K-Akt pathway in PDAC. Therefore, LETM2 may act as not only a prognostic but also a therapeutic target for early diagnosis and novel treatment of PDAC. We will be pursuing the feasibility of inhibitors targeting LETM2 such as delivering siLETM2s by smart nanoparticles in subsequent studies.

5. Conclusions

To summarize, we observed differential expression of LETM2 in various malignancies in the pan-cancer analysis, which was most prominently correlated with the dismal prognosis of PDAC. Results of in vivo as well as in vitro trials confirmed that LETM2 accelerates PDAC malignant progression via activating the downstream PI3K-Akt pathway. This is the first report on the functional role and molecular mechanism of the newly protein-coding gene, LETM2, in pancreatic cancer. We identified the LETM2-regulated PI3K-Akt signaling axis as the crucial element of PDAC development and progression, which provides effective targets for the diagnosis, treatment, and prognosis of PDAC.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers14194722/s1, Figure S1: Correlation analysis between LETM2 expression and signaling pathway. Figure S2: Correlation analysis between LETM2 expression and immune infiltration. Figure S3: Source Western blotting images for Figure 3A. Figure S4: The bioinformatic analysis of relative expression of LETM2 in PAAD cell lines via CCLE database. Figure S5: Source Western blotting images for Figure 4A. Figure S6: Source Western blotting images for Figure 4D. Figure S7: Source Western blotting images for Figure 5A. Figure S8: Source Western blotting images for Figure 6E. Table S1: Abbreviation and full name of cancer in TCGA + GTEx dataset. Table S2: Primers used in qRT–PCR analysis. Table S3: LETM2 siRNA sequences. Table S4: Pan-cancer analysis of LETM2 expression in TCGA dataset. Table S5: Pan-cancer analysis of LETM2 expression in TCGA + GTEx dataset. Table S6: Pan-cancer analysis of LETM2 expression in GEPIA database. Table S7: Source data for Figure 3A. Table S8: The correlation between IHC score of LETM2 and the clinical characteristic in PDAC patients. Table S9: Source data for Figure 4A. Table S10: Source data for Figure 4D. Table S11: Source data for Figure 5A. Table S12: Source data for Figure 6E.

Author Contributions

Conceptualization, S.Z. and Z.Z.; methodology, Y.L. (Yanzong Lu); software, H.Y.; validation, Y.Z., T.G. and K.Y.; writing—original draft preparation, Y.L. (Yunlong Li); writing—review and editing, Y.L. (Yaqing Li) and S.C.; supervision, K.H.; project administration, G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (grant no. 81874057, 82103142, and 82002505), and Medical Science and Technology Foundation of Guangdong Province (grant no. A2022371).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki. All enrolled PDAC patients obtained written consent and the collection of the clinical sample were approved by Institutional Review Board of Sun Yat-sen Memorial Hospital, Sun Yat-sen University at the date 24 December 2021 under the grant number: SYSEC-KY-KS-2021-337. The in vivo study was approved by the Institutional Animal Care and Use Committee (IACUC), Sun Yat-sen University at the date 26 February 2021 under the grant number: SYSU-IACUC-2021-B0269. Ethical review and approval were not applicable in bioinformatics and in vitro analysis performed in this study.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patient(s) to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank the patients who participated in our study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vincent, A.; Herman, J.; Schulick, R.; Hruban, R.H.; Goggins, M. Pancreatic cancer. Lancet 2011, 378, 607–620. [Google Scholar] [CrossRef]
  2. Raimondi, S.; Maisonneuve, P.; Lowenfels, A.B. Epidemiology of pancreatic cancer: An overview. Nat. Rev. Gastroenterol. Hepatol. 2009, 6, 699–708. [Google Scholar] [CrossRef] [PubMed]
  3. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2018. CA Cancer J. Clin. 2018, 68, 7–30. [Google Scholar] [CrossRef] [PubMed]
  4. Rahib, L.; Smith, B.D.; Aizenberg, R.; Rosenzweig, A.B.; Fleshman, J.M.; Matrisian, L.M. Projecting cancer incidence and deaths to 2030: The unexpected burden of thyroid, liver, and pancreas cancers in the United States. Cancer Res. 2014, 74, 2913–2921. [Google Scholar] [CrossRef]
  5. Neoptolemos, J.P.; Kleeff, J.; Michl, P.; Costello, E.; Greenhalf, W.; Palmer, D.H. Therapeutic developments in pancreatic cancer: Current and future perspectives. Nat. Rev. Gastroenterol. Hepatol. 2018, 15, 333–348. [Google Scholar] [CrossRef]
  6. Lin, Q.T.; Stathopulos, P.B. Molecular Mechanisms of Leucine Zipper EF-Hand Containing Transmembrane Protein-1 Function in Health and Disease. Int. J. Mol. Sci. 2019, 20, 286. [Google Scholar] [CrossRef] [PubMed]
  7. Endele, S.; Fuhry, M.; Pak, S.J.; Zabel, B.U.; Winterpacht, A. LETM1, a novel gene encoding a putative EF-hand Ca(2+)-binding protein, flanks the Wolf-Hirschhorn syndrome (WHS) critical region and is deleted in most WHS patients. Genomics 1999, 60, 218–225. [Google Scholar] [CrossRef] [PubMed]
  8. Tamai, S.; Iida, H.; Yokota, S.; Sayano, T.; Kiguchiya, S.; Ishihara, N.; Hayashi, J.; Mihara, K.; Oka, T. Characterization of the mitochondrial protein LETM1, which maintains the mitochondrial tubular shapes and interacts with the AAA-ATPase BCS1L. J. Cell Sci. 2008, 121 (Pt 15), 2588–2600. [Google Scholar] [CrossRef]
  9. Bock, F.J.; Tait, S.W.G. Mitochondria as multifaceted regulators of cell death. Nat. Rev. Mol. Cell. Biol. 2020, 21, 85–100. [Google Scholar] [CrossRef]
  10. Wallace, D.C. Mitochondria and cancer. Nat. Rev. Cancer 2012, 12, 685–698. [Google Scholar] [CrossRef] [Green Version]
  11. Austin, S.; Nowikovsky, K. LETM1: Essential for Mitochondrial Biology and Cation Homeostasis? Trends Biochem. Sci. 2019, 44, 648–658. [Google Scholar] [CrossRef] [PubMed]
  12. Tran, Q.; Lee, H.; Jung, J.H.; Chang, S.; Shrestha, R.; Kong, G.; Park, J.; Kim, S.; Park, K.; Rhee, H.; et al. Emerging role of LETM1/GRP78 axis in lung cancer. Cell Death Dis. 2022, 13, 543. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, J.; Ding, W.; Xu, Y.; Tao, E.; Mo, M.; Xu, W.; Cai, X.; Chen, X.; Yuan, J.; Wu, X. Long non-coding RNA RHPN1-AS1 promotes tumorigenesis and metastasis of ovarian cancer by acting as a ceRNA against miR-596 and upregulating LETM1. Aging 2020, 12, 4558–4572. [Google Scholar] [CrossRef] [PubMed]
  14. Piao, L.; Li, H.; Feng, Y.; Li, X.; Cui, Y.; Xuan, Y. Leucine Zipper-EF-Hand Containing Transmembrane Protein 1 Is a Potential Prognostic Biomarker and Promotes Cell Progression in Prostate Cancer. Cancer Manag. Res. 2020, 12, 1649–1660. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Chen, L.; Cao, Y.; Chen, S.; Xu, C.; Xing, J.; Zhang, K. LETM1 Promotes Gastric Cancer Cell Proliferation, Migration, and Invasion via the PI3K/Akt Signaling Pathway. J. Gastric Cancer 2020, 20, 139–151. [Google Scholar] [CrossRef]
  16. Stec, I.; van Ommen, G.J.; den Dunnen, J.T. WHSC1L1, on human chromosome 8p11.2, closely resembles WHSC1 and maps to a duplicated region shared with 4p16.3. Genomics 2001, 76, 5–8. [Google Scholar] [CrossRef]
  17. Dutt, A.; Ramos, A.H.; Hammerman, P.S.; Mermel, C.; Cho, J.; Sharifnia, T.; Chande, A.; Tanaka, K.E.; Stransky, N.; Greulich, H.; et al. Inhibitor-sensitive FGFR1 amplification in human non-small cell lung cancer. PLoS ONE 2011, 6, e20351. [Google Scholar] [CrossRef]
  18. Cheng, C.; Zhou, Y.; Li, H.; Xiong, T.; Li, S.; Bi, Y.; Kong, P.; Wang, F.; Cui, H.; Li, Y.; et al. Whole-Genome Sequencing Reveals Diverse Models of Structural Variations in Esophageal Squamous Cell Carcinoma. Am. J. Hum. Genet 2016, 98, 256–274. [Google Scholar] [CrossRef]
  19. Engelman, J.A.; Luo, J.; Cantley, L.C. The evolution of phosphatidylinositol 3-kinases as regulators of growth and metabolism. Nat. Rev. Genet. 2006, 7, 606–619. [Google Scholar] [CrossRef]
  20. Fruman, D.A.; Chiu, H.; Hopkins, B.D.; Bagrodia, S.; Cantley, L.C.; Abraham, R.T. The PI3K Pathway in Human Disease. Cell 2017, 170, 605–635. [Google Scholar] [CrossRef] [Green Version]
  21. Hoxhaj, G.; Manning, B.D. The PI3K-AKT network at the interface of oncogenic signalling and cancer metabolism. Nat. Rev. Cancer 2020, 20, 74–88. [Google Scholar] [CrossRef] [PubMed]
  22. Yang, J.; Nie, J.; Ma, X.; Wei, Y.; Peng, Y.; Wei, X. Targeting PI3K in cancer: Mechanisms and advances in clinical trials. Mol. Cancer 2019, 18, 26. [Google Scholar] [CrossRef] [PubMed]
  23. Conway, J.R.; Herrmann, D.; Evans, T.J.; Morton, J.P.; Timpson, P. Combating pancreatic cancer with PI3K pathway inhibitors in the era of personalised medicine. Gut 2019, 68, 742–758. [Google Scholar] [CrossRef] [PubMed]
  24. Schlieman, M.G.; Fahy, B.N.; Ramsamooj, R.; Beckett, L.; Bold, R.J. Incidence, mechanism and prognostic value of activated AKT in pancreas cancer. Br. J. Cancer 2003, 89, 2110–2115. [Google Scholar] [CrossRef]
  25. Kennedy, A.L.; Morton, J.P.; Manoharan, I.; Nelson, D.M.; Jamieson, N.B.; Pawlikowski, J.S.; McBryan, T.; Doyle, B.; McKay, C.; Oien, K.A.; et al. Activation of the PIK3CA/AKT pathway suppresses senescence induced by an activated RAS oncogene to promote tumorigenesis. Mol. Cell 2011, 42, 36–49. [Google Scholar] [CrossRef]
  26. Vanhaesebroeck, B.; Perry, M.W.D.; Brown, J.R.; Andre, F.; Okkenhaug, K. PI3K inhibitors are finally coming of age. Nat. Rev. Drug Discov. 2021, 20, 741–769. [Google Scholar] [CrossRef]
  27. Alzahrani, A.S. PI3K/Akt/mTOR inhibitors in cancer: At the bench and bedside. Semin. Cancer Biol. 2019, 59, 125–132. [Google Scholar] [CrossRef]
  28. Song, M.; Bode, A.M.; Dong, Z.; Lee, M.H. AKT as a Therapeutic Target for Cancer. Cancer Res. 2019, 79, 1019–1031. [Google Scholar] [CrossRef]
  29. Kent, W.J.; Sugnet, C.W.; Furey, T.S.; Roskin, K.M.; Pringle, T.H.; Zahler, A.M.; Haussler, D. The human genome browser at UCSC. Genome Res. 2002, 12, 996–1006. [Google Scholar]
  30. Li, C.; Tang, Z.; Zhang, W.; Ye, Z.; Liu, F. GEPIA2021: Integrating multiple deconvolution-based analysis into, G.E.P.I.A. Nucleic Acids Res. 2021, 49, W242–W246. [Google Scholar] [CrossRef]
  31. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [Green Version]
  32. Chen, S.; Chen, S.; Lian, G.; Li, Y.; Ye, X.; Zou, J.; Li, R.; Tan, Y.; Li, X.; Zhang, M.; et al. Development and validation of a novel nomogram for pretreatment prediction of liver metastasis in pancreatic cancer. Cancer Med. 2020, 9, 2971–2980. [Google Scholar] [CrossRef] [PubMed]
  33. Arnold, M.; Rutherford, M.J.; Bardot, A.; Ferlay, J.; Andersson, T.M.; Myklebust, T.A.; Tervonen, H.; Thursfield, V.; Ransom, D.; Shack, L.; et al. Progress in cancer survival, mortality, and incidence in seven high-income countries 1995-2014 (ICBP SURVMARK-2): A population-based study. Lancet Oncol. 2019, 20, 1493–1505. [Google Scholar] [CrossRef]
  34. Singhi, A.D.; Koay, E.J.; Chari, S.T.; Maitra, A. Early Detection of Pancreatic Cancer: Opportunities and Challenges. Gastroenterology 2019, 156, 2024–2040. [Google Scholar] [CrossRef]
  35. Pereira, S.P.; Oldfield, L.; Ney, A.; Hart, P.A.; Keane, M.G.; Pandol, S.J.; Li, D.; Greenhalf, W.; Jeon, C.Y.; Koay, E.J.; et al. Early detection of pancreatic cancer. Lancet Gastroenterol. Hepatol. 2020, 5, 698–710. [Google Scholar] [CrossRef]
  36. Fan, J.Q.; Wang, M.F.; Chen, H.L.; Shang, D.; Das, J.K.; Song, J. Current advances and outlooks in immunotherapy for pancreatic ductal adenocarcinoma. Mol. Cancer 2020, 19, 32. [Google Scholar] [CrossRef] [PubMed]
  37. Pole, J.C.; Courtay-Cahen, C.; Garcia, M.J.; Blood, K.A.; Cooke, S.L.; Alsop, A.E.; Tse, D.M.; Caldas, C.; Edwards, P.A. High-resolution analysis of chromosome rearrangements on 8p in breast, colon and pancreatic cancer reveals a complex pattern of loss, gain and translocation. Oncogene 2006, 25, 5693–5706. [Google Scholar] [CrossRef] [PubMed]
  38. Tabares-Seisdedos, R.; Rubenstein, J.L. Chromosome 8p as a potential hub for developmental neuropsychiatric disorders: Implications for schizophrenia, autism and cancer. Mol. Psychiatry 2009, 14, 563–589. [Google Scholar] [CrossRef]
  39. Servetto, A.; Kollipara, R.; Formisano, L.; Lin, C.C.; Lee, K.M.; Sudhan, D.R.; Gonzalez-Ericsson, P.I.; Chatterjee, S.; Guerrero-Zotano, A.; Mendiratta, S.; et al. Nuclear FGFR1 Regulates Gene Transcription and Promotes Antiestrogen Resistance in ER(+) Breast Cancer. Clin. Cancer Res. 2021, 27, 4379–4396. [Google Scholar] [CrossRef]
  40. Wang, K.; Ji, W.; Yu, Y.; Li, Z.; Niu, X.; Xia, W.; Lu, S. FGFR1-ERK1/2-SOX2 axis promotes cell proliferation, epithelial-mesenchymal transition, and metastasis in FGFR1-amplified lung cancer. Oncogene 2018, 37, 5340–5354. [Google Scholar] [CrossRef]
  41. Helsten, T.; Elkin, S.; Arthur, E.; Tomson, B.N.; Carter, J.; Kurzrock, R. The FGFR Landscape in Cancer: Analysis of 4,853 Tumors by Next-Generation Sequencing. Clin. Cancer Res. 2016, 22, 259–267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Irish, J.C.; Mills, J.N.; Turner-Ivey, B.; Wilson, R.C.; Guest, S.T.; Rutkovsky, A.; Dombkowski, A.; Kappler, C.S.; Hardiman, G.; Ethier, S.P. Amplification of WHSC1L1 regulates expression and estrogen-independent activation of ERalpha in SUM-44 breast cancer cells and is associated with ERalpha over-expression in breast cancer. Mol. Oncol. 2016, 10, 850–865. [Google Scholar] [CrossRef] [PubMed]
  43. Mahmood, S.F.; Gruel, N.; Nicolle, R.; Chapeaublanc, E.; Delattre, O.; Radvanyi, F.; Bernard-Pierrot, I. PPAPDC1B and WHSC1L1 are common drivers of the 8p11-12 amplicon, not only in breast tumors but also in pancreatic adenocarcinomas and lung tumors. Am. J. Pathol. 2013, 183, 1634–1644. [Google Scholar] [CrossRef]
  44. Toren, P.; Zoubeidi, A. Targeting the PI3K/Akt pathway in prostate cancer: Challenges and opportunities (review). Int. J. Oncol. 2014, 45, 1793–1801. [Google Scholar] [CrossRef] [PubMed]
  45. Duan, Y.; Haybaeck, J.; Yang, Z. Therapeutic Potential of PI3K/AKT/mTOR Pathway in Gastrointestinal Stromal Tumors: Rationale and Progress. Cancers 2020, 12, 2972. [Google Scholar] [CrossRef] [PubMed]
  46. Mehra, S.; Deshpande, N.; Nagathihalli, N. Targeting PI3K Pathway in Pancreatic Ductal Adenocarcinoma: Rationale and Progress. Cancers 2021, 13, 4434. [Google Scholar] [CrossRef]
  47. Folkes, A.J.; Ahmadi, K.; Alderton, W.K.; Alix, S.; Baker, S.J.; Box, G.; Chuckowree, I.S.; Clarke, P.A.; Depledge, P.; Eccles, S.A.; et al. The identification of 2-(1H-indazol-4-yl)-6-(4-methanesulfonyl-piperazin-1-ylmethyl)-4-morpholin-4-yl-t hieno [3,2-d]pyrimidine (GDC-0941) as a potent, selective, orally bioavailable inhibitor of class I PI3 kinase for the treatment of cancer. J. Med. Chem. 2008, 51, 5522–5532. [Google Scholar] [CrossRef]
  48. Jo, H.; Mondal, S.; Tan, D.; Nagata, E.; Takizawa, S.; Sharma, A.K.; Hou, Q.; Shanmugasundaram, K.; Prasad, A.; Tung, J.K.; et al. Small molecule-induced cytosolic activation of protein kinase Akt rescues ischemia-elicited neuronal death. Proc. Natl. Acad. Sci. USA 2012, 109, 10581–10586. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Expression of LETM2 and its prognostic value in pan-cancer analysis and PDAC. (A) Comparison of LETM2 expression in tumor tissue and adjacent normal tissue in the TCGA database. Cancers indicated in red or blue denotes have significantly different expressions of LETM2 in comparison with tumor tissues and normal tissues. (B) Comparison of LETM2 expression in tumor tissue and normal tissue in samples combined TCGA database and GTEx database. The abbreviations and full name of cancer mentioned in (A) and (B) are listed in Table S1. (C) Comparison of LETM2 expression in tumor tissue and normal tissue in the GEPIA database. Cancers indicated in red or blue denotes have significantly different expressions of LETM2 in comparison with tumor tissues and normal tissues. (D) Pan-cancer survival analysis of LETM2 based on the TCGA database. Statistical analysis was compared to the control or normal groups: * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001.
Figure 1. Expression of LETM2 and its prognostic value in pan-cancer analysis and PDAC. (A) Comparison of LETM2 expression in tumor tissue and adjacent normal tissue in the TCGA database. Cancers indicated in red or blue denotes have significantly different expressions of LETM2 in comparison with tumor tissues and normal tissues. (B) Comparison of LETM2 expression in tumor tissue and normal tissue in samples combined TCGA database and GTEx database. The abbreviations and full name of cancer mentioned in (A) and (B) are listed in Table S1. (C) Comparison of LETM2 expression in tumor tissue and normal tissue in the GEPIA database. Cancers indicated in red or blue denotes have significantly different expressions of LETM2 in comparison with tumor tissues and normal tissues. (D) Pan-cancer survival analysis of LETM2 based on the TCGA database. Statistical analysis was compared to the control or normal groups: * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001.
Cancers 14 04722 g001
Figure 2. Validation of LETM2 prognostic value and clinical significance for pancreatic cancer in TCGA database. (A) TCGA + GTEx database analysis revealed that LETM2 expression was significantly elevated in PDAC tumor tissue in comparison with normal tissue. Red represents tumor tissue and blue represents normal tissue. (B,C) High levels of LETM2 expression correlated with significantly poor (B) OS and (C) DFS of PDAC patients based on survival analysis of the TCGA database. The dotted lines represented 95% confidence intervals. (D) LETM2 expression correlated with tumor inflammation signature. (E) LETM2 expression correlated with PI3K-AKT-mTOR pathway. (F) LETM2 expression correlated with the P53 pathway. (G) LETM2 expression correlated with immune infiltration level of B cell expression. (H) LETM2 expression correlated with immune infiltration level of T cell CD8+ expression. (I) Sankey diagram demonstrated the association between LETM2 expression and histological grade and living status. * p < 0.05.
Figure 2. Validation of LETM2 prognostic value and clinical significance for pancreatic cancer in TCGA database. (A) TCGA + GTEx database analysis revealed that LETM2 expression was significantly elevated in PDAC tumor tissue in comparison with normal tissue. Red represents tumor tissue and blue represents normal tissue. (B,C) High levels of LETM2 expression correlated with significantly poor (B) OS and (C) DFS of PDAC patients based on survival analysis of the TCGA database. The dotted lines represented 95% confidence intervals. (D) LETM2 expression correlated with tumor inflammation signature. (E) LETM2 expression correlated with PI3K-AKT-mTOR pathway. (F) LETM2 expression correlated with the P53 pathway. (G) LETM2 expression correlated with immune infiltration level of B cell expression. (H) LETM2 expression correlated with immune infiltration level of T cell CD8+ expression. (I) Sankey diagram demonstrated the association between LETM2 expression and histological grade and living status. * p < 0.05.
Cancers 14 04722 g002
Figure 3. Elevated LETM2 expression correlates with clinicopathological features and poor prognosis in PDAC patients. (A) The protein expression level of LETM2 was examined by Western blotting in paired samples of PDAC tumor and non-tumor tissues. The densitometry readings/intensity ratio of each band and the original whole blot could be found in Table S7 and Figure S3. (B) Quantitative analysis of Western blotting and (C) qRT-PCR analysis revealed LETM2 protein expression was remarkably increased in PDAC tumor tissues. The two-tailed t-test in paired samples. (D) Representative images of IHC staining for LETM2 protein of PDAC tumor tissues and their adjacent non-tumor tissues. Scale bars, left: 600 μm, right: 100 μm. (E) Quantification of IHC staining score. (FH) Relationship between IHC score of LETM2 and clinicopathological features of PDAC patients, such as (F) T stage, (G) distant metastasis and (H) pathological stages. Overall survival (I) and disease-free survival (J) analyses in PDAC patients in the high LETM2 expression group and low LETM2 expression group were conducted by Kaplan–Meier method with the two-tailed log-rank test. The data are represented as mean ± SEM, and Student’s t-test was performed. ns, not significant; * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001.
Figure 3. Elevated LETM2 expression correlates with clinicopathological features and poor prognosis in PDAC patients. (A) The protein expression level of LETM2 was examined by Western blotting in paired samples of PDAC tumor and non-tumor tissues. The densitometry readings/intensity ratio of each band and the original whole blot could be found in Table S7 and Figure S3. (B) Quantitative analysis of Western blotting and (C) qRT-PCR analysis revealed LETM2 protein expression was remarkably increased in PDAC tumor tissues. The two-tailed t-test in paired samples. (D) Representative images of IHC staining for LETM2 protein of PDAC tumor tissues and their adjacent non-tumor tissues. Scale bars, left: 600 μm, right: 100 μm. (E) Quantification of IHC staining score. (FH) Relationship between IHC score of LETM2 and clinicopathological features of PDAC patients, such as (F) T stage, (G) distant metastasis and (H) pathological stages. Overall survival (I) and disease-free survival (J) analyses in PDAC patients in the high LETM2 expression group and low LETM2 expression group were conducted by Kaplan–Meier method with the two-tailed log-rank test. The data are represented as mean ± SEM, and Student’s t-test was performed. ns, not significant; * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001.
Cancers 14 04722 g003
Figure 4. LETM2 promotes PDAC cell proliferation, migration, and invasion in vitro. (A) Western blotting analysis of LETM2 protein expression in three types PDAC cell lines (MIA PaCa-2, SW1990, and BxPc-3) compared with human pancreatic duct cells (hTERT-HPNE). The densitometry readings/intensity ratio of each band and the original whole blot could be found in Table S9 and Figure S5. (B) Quantitative analysis of Western blotting and (C) qRT-PCR analysis revealed LETM2 protein and mRNA expression were significantly higher in PDAC cell lines. (D) Western blotting analysis for verifying the efficiency of LETM2 overexpression. The densitometry readings/intensity ratio of each band and the original whole blot could be found in Table S10 and Figure S6. (E) Quantitative analysis of Western blotting and (F) qRT-PCR analysis revealed that LETM2 protein and mRNA expression were significantly elevated in overexpression cells. oeLETM2 indicates MIA PaCa-2 cells and SW1990 cells transfected with an overexpressing plasmid. The relative quantification was calculated using the 2−ΔΔCt method and normalized based on GAPDH. (G) Cell apoptosis of LETM2 overexpression cell lines was measured using flow cytometry. (H) Cell migratory and invasive ability of LETM2 overexpression cell lines were measured using transwell assay. Representative photographs (left) and quantification (right) were shown. Scale bars, 50 μm. (I) Cell growth of LETM2 overexpression cell lines was measured using CCK-8 assays. The data are represented as mean ± SEM (n = 3 biologically independent experiments), and Student’s t-test was performed. * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 4. LETM2 promotes PDAC cell proliferation, migration, and invasion in vitro. (A) Western blotting analysis of LETM2 protein expression in three types PDAC cell lines (MIA PaCa-2, SW1990, and BxPc-3) compared with human pancreatic duct cells (hTERT-HPNE). The densitometry readings/intensity ratio of each band and the original whole blot could be found in Table S9 and Figure S5. (B) Quantitative analysis of Western blotting and (C) qRT-PCR analysis revealed LETM2 protein and mRNA expression were significantly higher in PDAC cell lines. (D) Western blotting analysis for verifying the efficiency of LETM2 overexpression. The densitometry readings/intensity ratio of each band and the original whole blot could be found in Table S10 and Figure S6. (E) Quantitative analysis of Western blotting and (F) qRT-PCR analysis revealed that LETM2 protein and mRNA expression were significantly elevated in overexpression cells. oeLETM2 indicates MIA PaCa-2 cells and SW1990 cells transfected with an overexpressing plasmid. The relative quantification was calculated using the 2−ΔΔCt method and normalized based on GAPDH. (G) Cell apoptosis of LETM2 overexpression cell lines was measured using flow cytometry. (H) Cell migratory and invasive ability of LETM2 overexpression cell lines were measured using transwell assay. Representative photographs (left) and quantification (right) were shown. Scale bars, 50 μm. (I) Cell growth of LETM2 overexpression cell lines was measured using CCK-8 assays. The data are represented as mean ± SEM (n = 3 biologically independent experiments), and Student’s t-test was performed. * p < 0.05; ** p < 0.01; *** p < 0.001.
Cancers 14 04722 g004
Figure 5. Downregulation of LETM2 inhibits PDAC cell proliferation and metastasis and promoted cell apoptosis in vitro. (A) Western blotting analysis for verifying the efficiency of LETM2 knockdown. The densitometry readings/intensity ratio of each band and the original whole blot could be found in Table S11 and Figure S7. (B) Quantitative analysis of Western blotting and (C) qRT-PCR analysis revealed LETM2 protein and mRNA expression were significantly downregulated in knockdown cells. siLETM2 indicates BxPc-3 cells and SW1990 cells transfected with siRNA targeting LETM2. The relative quantification was calculated using the 2−ΔΔCt method and normalized based on GAPDH. (D) Cell apoptosis of LETM2 knockdown cell lines was measured using flow cytometry. (E) Cell growth of LETM2 knockdown cell lines was measured using CCK-8 assays. (F) Cell migratory and invasive ability of LETM2 knockdown cell lines were measured using transwell assay. Representative photographs (left) and quantification (right) were shown. Scale bars, 50 μm. The data are represented as mean ± SEM (n = 3 biologically independent experiments), and Student’s t-test was performed. * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 5. Downregulation of LETM2 inhibits PDAC cell proliferation and metastasis and promoted cell apoptosis in vitro. (A) Western blotting analysis for verifying the efficiency of LETM2 knockdown. The densitometry readings/intensity ratio of each band and the original whole blot could be found in Table S11 and Figure S7. (B) Quantitative analysis of Western blotting and (C) qRT-PCR analysis revealed LETM2 protein and mRNA expression were significantly downregulated in knockdown cells. siLETM2 indicates BxPc-3 cells and SW1990 cells transfected with siRNA targeting LETM2. The relative quantification was calculated using the 2−ΔΔCt method and normalized based on GAPDH. (D) Cell apoptosis of LETM2 knockdown cell lines was measured using flow cytometry. (E) Cell growth of LETM2 knockdown cell lines was measured using CCK-8 assays. (F) Cell migratory and invasive ability of LETM2 knockdown cell lines were measured using transwell assay. Representative photographs (left) and quantification (right) were shown. Scale bars, 50 μm. The data are represented as mean ± SEM (n = 3 biologically independent experiments), and Student’s t-test was performed. * p < 0.05; ** p < 0.01; *** p < 0.001.
Cancers 14 04722 g005
Figure 6. LETM2 may activate the downstream PI3K-Akt signaling pathway by increasing the phosphorylation of PI3K and Akt. (A) KEGG pathway enrichment analyses were performed on LETM2 and co-expressed genes in PDAC. (B) The molecular functions (MF) of GO functional enrichment analyses were performed on LETM2 and co-expressed genes in PDAC. (C) The biological processes (BP) of GO functional enrichment analyses were performed on LETM2 and co-expressed genes in PDAC. (D) The cellular components (CC) of GO functional enrichment analyses were performed on LETM2 and co-expressed genes in PDAC. (E) Western blotting analysis for evaluating the key protein level of PI3K-Akt signaling pathway after LETM2 knockdown. The densitometry readings/intensity ratio of each band and the original whole blot could be found in Table S12 and Figure S8. (F) Quantitation of Western blotting implied that p-PI3K and p-Akt were significantly inhibited after LETM2 knockdown. Student’s t-test was performed. ns, not significant; ** p < 0.01.
Figure 6. LETM2 may activate the downstream PI3K-Akt signaling pathway by increasing the phosphorylation of PI3K and Akt. (A) KEGG pathway enrichment analyses were performed on LETM2 and co-expressed genes in PDAC. (B) The molecular functions (MF) of GO functional enrichment analyses were performed on LETM2 and co-expressed genes in PDAC. (C) The biological processes (BP) of GO functional enrichment analyses were performed on LETM2 and co-expressed genes in PDAC. (D) The cellular components (CC) of GO functional enrichment analyses were performed on LETM2 and co-expressed genes in PDAC. (E) Western blotting analysis for evaluating the key protein level of PI3K-Akt signaling pathway after LETM2 knockdown. The densitometry readings/intensity ratio of each band and the original whole blot could be found in Table S12 and Figure S8. (F) Quantitation of Western blotting implied that p-PI3K and p-Akt were significantly inhibited after LETM2 knockdown. Student’s t-test was performed. ns, not significant; ** p < 0.01.
Cancers 14 04722 g006
Figure 7. LETM2 accelerates PDAC malignant progression by activating the PI3K-Akt pathway in vitro. (A) Cell apoptosis was measured using flow cytometry. (B) Cell migratory and invasive ability were measured using transwell assay. Representative photographs (left) and quantification (right) were shown. Scale bars, 50 μm. (C) Cell growth was measured using CCK-8 assays. The data are represented as mean ± SEM (n = 3 biologically independent experiments), and Student’s t-test was performed. Ns, not significant; * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 7. LETM2 accelerates PDAC malignant progression by activating the PI3K-Akt pathway in vitro. (A) Cell apoptosis was measured using flow cytometry. (B) Cell migratory and invasive ability were measured using transwell assay. Representative photographs (left) and quantification (right) were shown. Scale bars, 50 μm. (C) Cell growth was measured using CCK-8 assays. The data are represented as mean ± SEM (n = 3 biologically independent experiments), and Student’s t-test was performed. Ns, not significant; * p < 0.05; ** p < 0.01; *** p < 0.001.
Cancers 14 04722 g007
Figure 8. LETM2 accelerates PDAC malignant progression by activating the PI3K-Akt pathway in vivo. (A) Temporal scheme of in vivo experiments. (B) Gross images of PDAC subcutaneous tumor. (C) The tumors were extracted and weighed after 4 weeks. (D) The tumor volume was monitored every five days, and tumor growth curves were generated. (E) H&E and IHC staining with an antibody specific for Ki67 in sections of tumors. The data shown were representative of three experiments. Scale bars, 200 μm. Student’s t-test was performed. ns, not significant; ** p < 0.01; *** p < 0.001.
Figure 8. LETM2 accelerates PDAC malignant progression by activating the PI3K-Akt pathway in vivo. (A) Temporal scheme of in vivo experiments. (B) Gross images of PDAC subcutaneous tumor. (C) The tumors were extracted and weighed after 4 weeks. (D) The tumor volume was monitored every five days, and tumor growth curves were generated. (E) H&E and IHC staining with an antibody specific for Ki67 in sections of tumors. The data shown were representative of three experiments. Scale bars, 200 μm. Student’s t-test was performed. ns, not significant; ** p < 0.01; *** p < 0.001.
Cancers 14 04722 g008
Table 1. Univariate and multivariate Cox regression analyses of overall survival of PDAC in TCGA database.
Table 1. Univariate and multivariate Cox regression analyses of overall survival of PDAC in TCGA database.
CharacteristicsnUnivariate AnalysisMultivariate Analysis
Hazard Ratio (95% CI)p ValueHazard Ratio (95% CI)p Value
Age178
≤6593Reference
>65851.290 (0.854–1.948)0.227
Gender178
Female80Reference
Male980.809 (0.537–1.219)0.311
T stage176
T1 + T231Reference
T3 + T41452.023 (1.072–3.816)0.030 *1.322 (0.658–2.654)0.433
N stage173
N050Reference
N11232.154 (1.282–3.618)0.004 *2.071 (1.179–3.637)0.011 *
M stage84
M079Reference
M150.756 (0.181–3.157)0.701
Pathologic stage175
Stage I + II167Reference
Stage III + IV80.673 (0.212–2.135)0.501
Histologic grade176
G1 + G2126Reference
G3 + G4501.538 (0.996–2.376)0.0521.196 (0.764–1.871)0.434
LETM2 expression1781.524 (0.989–2.350)0.0561.658 (1.052–2.613)0.030 *
CI, confidence interval; variables with a p-value < 0.1 in univariate analyses were included in the multivariate analyses. * p < 0.05 indicates significance.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, S.; Zhong, Z.; Lu, Y.; Li, Y.; Yao, H.; Zhao, Y.; Guo, T.; Yang, K.; Li, Y.; Chen, S.; et al. A LETM2-Regulated PI3K-Akt Signaling Axis Reveals a Prognostic and Therapeutic Target in Pancreatic Cancer. Cancers 2022, 14, 4722. https://doi.org/10.3390/cancers14194722

AMA Style

Zhou S, Zhong Z, Lu Y, Li Y, Yao H, Zhao Y, Guo T, Yang K, Li Y, Chen S, et al. A LETM2-Regulated PI3K-Akt Signaling Axis Reveals a Prognostic and Therapeutic Target in Pancreatic Cancer. Cancers. 2022; 14(19):4722. https://doi.org/10.3390/cancers14194722

Chicago/Turabian Style

Zhou, Shurui, Ziyi Zhong, Yanzong Lu, Yunlong Li, Hanming Yao, Yue Zhao, Tairan Guo, Kege Yang, Yaqing Li, Shaojie Chen, and et al. 2022. "A LETM2-Regulated PI3K-Akt Signaling Axis Reveals a Prognostic and Therapeutic Target in Pancreatic Cancer" Cancers 14, no. 19: 4722. https://doi.org/10.3390/cancers14194722

APA Style

Zhou, S., Zhong, Z., Lu, Y., Li, Y., Yao, H., Zhao, Y., Guo, T., Yang, K., Li, Y., Chen, S., Huang, K., & Lian, G. (2022). A LETM2-Regulated PI3K-Akt Signaling Axis Reveals a Prognostic and Therapeutic Target in Pancreatic Cancer. Cancers, 14(19), 4722. https://doi.org/10.3390/cancers14194722

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