Abstract

TXLNA (taxilin alpha), a binding partner of the syntaxin family, was identified as a key factor in the coordination of intracellular vesicle trafficking and highly expressed in various tumor cells. However, the accurate relation between TXLNA and tumorigenesis and progression of pancreatic adenocarcinoma (PAAD) is still unclear. The present study was designed to examine the expression profile of TXLNA and explore its prognostic significance in PAAD patients and the possible molecular regulatory mechanism by analyzing a series of data from databases, including GEPIA, LOGpc, STRING, and GeneMANIA. The results indicate that TXLNA mRNA and protein were remarkably increased in PAAD tissues compared with normal pancreatic tissues. The high TXLNA expression was significantly correlated with superior overall survival (OS), disease-free interval (DFI), disease specific survival (DSS), and progression-free interval (PFI) for PAAD patients. In summary, high TXLNA expression could predict favourable OS, DFI, DSS, and PFI for PAAD patients, and it might be as potential prognostic biomarkers and targets for PAAD.

1. Introduction

Pancreatic cancer (PC) is one of the leading causes of cancer mortality, resulting in a substantial global burden [1]. It has become the third leading cause of cancer death in the USA, with a 5-year survival rate of 8% [2]. In recent decades, the treatment of pancreatic cancer has been improved; however, the diagnostic and prognostic biomarkers for pancreatic cancer are still deficient due to its heterogeneous characteristics [3]. Up to now, carbohydrate antigen 19-9 (CA19-9) is the only available diagnostic biomarker for pancreatic adenocarcinoma (PAAD) approved by the Food and Drug Administration (FDA) in the U.S., but it had a limitation of poor sensitivity and specificity [4, 5]. Therefore, it is urgent to search for novel biomarkers.

Taxilin alpha (TXLNA) was identified as a binding partner of the syntaxin family [6]. The TXLNA was also named interleukin 14 (IL-14), and its alternative gene name was designated as TXLNA or IL14 in the NCBI. The function of TXLNA was identified as a key factor in the coordination of intracellular vesicle trafficking [6]. The syntaxin family was composed of three members, including alpha (α)-, beta (β)-, and gamma (γ)-taxilin. In addition, they share an extraordinarily long coiled-coil region homologous to that of Uso1, a yeast tethering factor homolog of p115 [7]. It was reported that the syntaxin family involved in the transport of the vesicles delivered to the plasma membrane [6, 7]. As we know, the intracellular vesicle tranfficking is a primary transportation system in eukaryotic cells, and vesicle tranfficking participates in the multiple processes of cell proliferation. The abnormal proliferation indicates a high probability of tumor or cancer [8].

Recent studies demonstrate that the TXLNA protein had a high expression in various tumor cells and its expression indicates a close relationship with histological grade and proliferative activity in hepatocellular carcinoma, as well as the metastatic and invasive potential of renal cell carcinoma [811]. Sakakibara et al. showed that a strong expression of TXLNA was detected in proliferating neuroepithelial cells during embryonic development of the rat brain [12]. Horii et al. have found that the expression of TXLNA was correlated with cell proliferation in the adult murine gastrointestinal tract and during postnatal development [9].

The relationship between TXLNA and pancreatic cancer prognosis remains to be determined. The present study was designed to evaluate the expression profile of TXLNA and to investigate its prognostic significance for PAAD patients.

2. Materials and Methods

2.1. TXLNA Gene Expression Analysis

The TXLNA mRNA level in PAAD tissue and normal pancreatic tissue was examined by on-line analysis GEPIA (Gene Expression Profiling Interactive Analysis) (http://gepia.cancer-pku.cn). GEPIA is a web server which analyzes the RNA expression data of 9,736 tumors and 8,587 normal samples from TCGA and GTEx [13].

2.2. Prognosis Analysis

The relationship between the expression of TXLNA and survival terms, including overall survival (OS), disease-free interval (DFI), disease specific survival (DSS), and progression-free interval (PFI), in 542 PAAD patients was determined by the prognosis analysis web server LOGpc [14] (Long-term Outcome and Gene Expression Profiling Database of pan-cancers, http://bioinfo.henu.edu.cn/DatabaseList.jsp) using TCGA and ICGA_Seq as we previously reported [15]. This LOGpc tool encompasses 171 expression datasets and provides 13 types of survival terms for 27929 patients of 26 distinct malignancies.

2.3. Gene Interaction Network Construction

Search Tool for the Retrieval of Interacting Genes (STRING, version 11.0, available online: http://www.string-db.org/) is an open-source database, which was applied to predict protein-protein interaction (PPI) networks of differentially expressed genes (DEGs) [16]. In addition, only the interactions with a combined score >0.4 were considered to have a statistically significant difference. The present new version contains 5090 organisms, 24584628 proteins, and 3123056667 total interactions. GeneMANIA (http://genemania.org/) was utilized to conduct analysis of TXLNA correlated genes. GeneMANIA is a prediction server to search for other genes related to the targeted genes using analysis of functional association, including protein and genetic interactions, pathways, coexpression, colocalization, and protein domain similarity [17].

2.4. Statistical Analysis

Univariate and multivariate Cox regression models were utilized to assess the prognostic role of TXLNA in terms of OS, DFI, PFI, and DSS for PAAD patients using SPSS. Risk factors () analyzed by univariate analysis were selected for multivariate Cox regression analysis. A value of was considered statistically significant.

3. Results

3.1. TXLNA Expression in PAAD Tissue Was Increased Compared with Normal Tissue

As shown in Figure 1(a), the GEPIA analysis based on TCGA and the GTEx data indicates that PAAD tissues exhibited a significantly increased TXLNA expression, compared with the normal tissue (). There was no significant difference of TXLNA expression between the TNM stages (, Figure 1(b)) and different neoplasm histologic grades (, Figure S1).

3.2. High TXLNA Expression Was Correlated with Favourable OS, DFI, DSS, and PFI in PAAD Patients

The prognostic value of TXLNA was tested by KaplanMeier Plotter Analysis of LOGpc [18, 19]. Analysis based on TCGA data in LOGpc showed that the high TXLNA expression group had favourable OS (; HR: 0.5519; 95% CI: 0.338–0.9013; Figure 2(a)), DFI (; HR: 0.4971; 95% CI: 0.2832–0.8725; Figure 2(b)), DSS (; HR: 0.5549; 95% CI: 0.3079–1.0001; Figure 2(c)), and PFI (; HR: 0.5515; 95% CI: 0.3409–0.8921; Figure 2(d)) for PAAD patients. Meanwhile, analysis results from ICGA_Seq data demonstrated that the increased TXLNA expression group had favourable OS (; HR: 0.5786; 95% CI: 0.4128–0.8109; Figure 2(e)) and favourable DFI (; HR: 0.6577; 95% CI: 0.4367–0.9904; Figure 2(f)) for PAAD patients. These results indicate that high TXLNA expression was significantly correlated with favourable OS, DFI, DSS, and PFI in PAAD patients.

3.3. Significant Correlation between the Increased TXLNA Expression and Favourable OS in PAAD Patients Was Restricted to Female Patients or Patients with Lymph Nodes, No-Smoking History, or Alcohol History

To examine whether the correlation was influenced by the different clinical characteristics of PAAD patients, such as gender, lymph nodes, and smoking or alcohol, each character was used to divide PAAD patients into two subgroups (male and female, lymph nodes and no-lymph nodes, smoking and no-smoking, and alcohol and no-alcohol) for analysis. The high TXLNA expression group had favourable OS in female PAAD patients (; HR: 0.3509; 95% CI: 0.1698–0.7248; Figure 3(a)), but not in male PAAD patients (; HR: 0.7952; 95% CI: 0.4053–1.5602; Figure 3(b)). The increased TXLNA expression group had favourable OS in PAAD patients with lymph nodes (; HR: 0.5773; 95% CI: 0.3531–0.9437; Figure 3(c)), but not in PAAD patients with no-lymph nodes (; Figure 3(d)). The high TXLNA expression group had favourable OS in PAAD patients with no-smoking history (; HR: 0.366; 95% CI: 0.1701–0.7872; Figure 3(e)), but not in PAAD patients with smoking history (; HR: 0.6551; 95% CI: 0.2828–1.5174; Figure 3(f)). The increased TXLNA expression group had favourable OS in PAAD patients with alcohol history (; HR: 0.5088; 95% CI: 0.2625–0.9861; Figure 3(g)), but not in PAAD patients with no-alcohol history (; HR: 0.5048; 95% CI: 0.2211–1.1525; Figure 3(h)). These results demonstrated that high TXLNA expression indicated a favourable OS in female patients or patients with lymph nodes, no-smoking history, or alcohol history.

3.4. TXLNA Expression Was an Independent Prognostic Biomarker for Favourable OS, DFI, and PFI in PAAD

Univariate analysis demonstrated that the lymph nodes and increased TXLNA expression correlated with favourable OS, DFI, PFI, and DSS for PAAD (Table 1). Multivariate analysis indicated that high TXLNA expression was an independent prognostic indicator in terms of OS (HR: 1.641; 95% CI: 1.003–2.685; ), DFI (HR: 1.848; 95% CI: 1.051–3.250; ), and PFI (HR: 1.668; 95% CI: 1.029–2.702; ; Table 1).

3.5. Construction of Protein-Protein and Gene-Gene Interaction Network

To explore the potentially interacted proteins with TXLNA, the PPI network was constructed by STRING. After removing unconnected nodes, the PPI network containing 41 nodes and 310 edges was established (Figure 4(a)). To further investigate the hub proteins in the network, the top 10 hub proteins were retrieved: STX3, STX4, STX1A, STXBP2, VAMP2, SNAP25, TXLNG, DDX17, TXLNB, and QARS (Figure 4(b)). Gene-gene interaction analysis using GeneMANIA showed 20 correlated genes with TXLNA, such as TXLNB, TXLNG, NACA, TTC27, IKBKG, TXNDC12, WIPI2, BRCA1, VAMP2, and KIF3C (Figure 4(c)).

4. Discussion

Pancreatic cancer is a highly heterogeneous disease with a low overall survival rate [2, 3]. Therefore, it is essential to discover effective tumor biomarkers to improve the clinical diagnosis, prognosis, and treatment. Recently, a variety of prognostic markers for pancreatic cancer are reported. Li et al. had found that the combination of SRPX2 and RAB31 were independent prognostic factors that associated with OS and DFS of pancreatic cancer [20]. Chen et al. reported that three hypomethylated genes (SULT1E1, IGF2BP3, and MAP4K4) were associated with poor overall survival of pancreatic cancer patients [21]. Pang et al. have indicated that IL22RA1, BCL2L1, STAT1, MYC, and STAT2 involved in the Jak-STAT signaling pathway may be significantly associated with prognosis of pancreatic cancer [22].

Recent studies indicate that TXLNA was involved in inflammation, endocrine, and respiratory diseases. IL-14, a cytokine, was identified as high molecular weight B-cell growth factor [23]. Its mRNA transcript level was increased in vitro in human T cells responding to alloantigen stimulation [24]. The IL-14, IL-17, and TNF-α level in whole lung homogenates were upregulated in RAGE transgenic mice [25]. Helen Kemp et al. had demonstrated that the thyroid follicular cells in the thyroid had the ability to express IL-14 mRNA [26]. Nevertheless, the role of TXLNA on pancreatic cancer remains unknown.

For PAAD patients, prognostic factors could help to guide the personalized treatments. The correlation of TXLNA expression and prognostic value was analyzed using LOGpc. We discovered that increased TXLNA expression was significantly correlated with favourable OS, DFI, PFI, and DSS in PAAD patients in TCGA data. Analysis from ICGA_Seq data also indicated that the upregulated TXLNA expression was significantly correlated with favourable OS and DFI. In addition, further analysis suggested that the correlation of high TXLNA expression and favourable OS was influenced by gender, smoking, alcohol, and lymph history in patients with PAAD. Gender was a key factor for prognostic analysis. Previous studies indicated that gender had a significant influence on the outcome of prognostic analysis, such as disease-free survival, and distant metastasis [2729]. Our study showed that the high TXLNA expression group exhibited favourable OS in female PAAD patients, but not in male. The gender differences in the prognosis of PAAD could be explained by genetic variants affected by the hormonal environment. It suggests that high TXLNA expression was suitable for prognosis of female (not male) PAAD patients. Multivariate analysis showed that high TXLNA expression was an independent prognostic indicator in terms of OS, DFI, and PFI in PAAD patient.

In conclusion, TXLNA has a high expression in PAAD tissue compared with normal tissue, and its upregulated expression might serve as an independent prognostic indicator of OS, DFI, DSS, and PFI in PAAD patients. Further experimental studies are needed to fully understand the association between TXLNA expression and prognosis of PAAD patients.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

XG and YY developed the idea and designed the research. SL, GZ, and LX collected and analyzed the data. SL, ZY, and QW wrote the draft of the manuscript. YL, LZ, YH, HL, and YD contributed to revise the writing. All authors read and approved the submitted version. Shuangyu Lv and Guosen Zhang authors have contributed equally to this work.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 81602362, 81600974, and 81971280), the Key Science and Technology Program of Henan Province in China (no. 162102310391 and 192102310080), the supporting grants of Henan University (no. 2015YBZR048 and B2015151), the Research Program for Young Talent of Henan University School of Medicine (Grant no. 2019018), the Program for Innovative Talents of Science and Technology in Henan Province (no. 18HASTIT048), the Yellow River Scholar Program (no. H2016012), the Kaifeng Science and Technology Major Project (no. 18ZD008), and Key Science and Technology Program of Kaifeng City (no. 1803034 and 1903019).

Supplementary Materials

Supplementary Figure 1: the TXLNA gene expression in PAAD patients with different grades. The histogram indicates the expression level of the TXLNA gene at different neoplasm histologic grades (G1, G2, G3, G4, and GX) in patients with PAAD in the TCGA. The values were expressed as the mean ± S.E.M. A one-way analysis of variance (ANOVA) followed by the Tukey’s test for post hoc comparisons was used for the statistical evaluation. (Supplementary Materials)