Overexpressing PLOD Family Genes Predict Poor Prognosis in Gastric Cancer Shan-Shan Li1*, Yi-Fan Lian2*, Yan-Lin Huang3, Yue-Hua Huang2,3, Jian Xiao1
Total Page:16
File Type:pdf, Size:1020Kb
Journal of Cancer 2020, Vol. 11 121 Ivyspring International Publisher Journal of Cancer 2020; 11(1): 121-131. doi: 10.7150/jca.35763 Research Paper Overexpressing PLOD family genes predict poor prognosis in gastric cancer Shan-Shan Li1*, Yi-Fan Lian2*, Yan-Lin Huang3, Yue-Hua Huang2,3, Jian Xiao1 1. Department of Medical Oncology, the Sixth Affiliated Hospital of Sun Yat-sen University, Guangzhou, China 2. Guangdong Provincial Key Laboratory of Liver Disease Research, the Third Affiliated Hospital of Sun Yat-sen University, Guangzhou, China 3. Department of Infectious Diseases, the Third Affiliated Hospital of Sun Yat-sen University, Guangzhou, China *These authors contributed equally to this work Corresponding authors: Jian Xiao, MD, PhD, Department of Medical Oncology, the Sixth Affiliated Hospital of Sun Yat-sen University, 26 Yuancun Erheng Road, Guangzhou, China, 510655, Tel: 8620-38285497, Email: [email protected]; Yue-Hua Huang, MD, PhD, Guangdong Provincial Key Laboratory of Liver Disease Research & Department of Infectious Diseases, the Third Affiliated Hospital of Sun Yat-sen University, 600 Tian He Rd. Guangzhou, China, 510630, Tel: 8620-85252702, Fax: 8620-85253305, Email: [email protected]. © The author(s). This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions. Received: 2019.04.15; Accepted: 2019.09.20; Published: 2020.01.01 Abstract Procollagen-lysine, 2-oxoglutarate 5-dioxygenases (PLODs) are a set of enzymes involved in the hydroxylation of lysine and stabilization of collagen by crosslinks. Previous studies have highlighted that overexpressing PLOD genes were related to the progression, migration and progression of different human cancers. However, the diverse expression patterns and prognostic values of PLOD genes remain to be elucidated in gastric cancer (GC). In this study, we mined the expression and survival data in GC patients through ONCOMINE, UALCAN and Kaplan-Meier Plotter database. STRING portal couple with DAVID was used to establish a functional protein interaction network of PLOD family genes and analyze the GO and KEGG enriched pathways. Differential gene expression correlated with PLOD family genes was identified with LinkedOmics. We found that PLOD1, 2 and 3 were up-regulated in GC patients compared with normal tissues. High expression levels of PLOD1 and PLOD3 were associated with shorter overall survival (OS), first progression (FP) and post progression survival (PPS) while high expression level of PLOD2 was only associated with shorter FP in all GC patients. Specifically, only high PLOD2 expression had significant correlation with shorter OS, FP and PPS in the diffuse type GC patients. Furthermore, combinatorial use of expressions of all PLOD genes was a superior prognostic indicator for GC patients. Pathway analysis confirmed that PLOD family genes mainly participate in regulating the collagen metabolism and extracellular matrix constitution, and the cellular adaptor protein SHC1, which helps to transduce an extracellular signal into an intracellular signal, could be the regulatory module mediating PLOD’s effect on GC. Therefore, we propose that individual PLOD genes or PLOD family genes as a whole could be potential prognostic biomarkers for GC. Key words: PLOD; gastric cancer; prognosis; bioinformatics Introduction As a major component of tumor micro- the balance of ECM signal will no doubt cause effect environment, extracellular matrix (ECM) plays an on cancer cell functions [3, 4]. Thus, identifying important role in cancer development and molecules correlated with ECM regulation would progression, which consists of structural proteins, help to develop novel diagnostic and therapeutic glycoproteins and proteoglycans. ECM proteins strategies for cancer treatment. deposition and crosslink provides necessary Among the various ECM components, collagens biochemical and biophysical supports for cancer cell are the most abundant proteins, which were reported proliferation, migration and invasion [1, 2]. Aberrant to modulate cancer cell migration, invasion [5], expression of ECM-related molecules and alteration of proliferation [6] and patient survival [7]. The collagen http://www.jcancer.org Journal of Cancer 2020, Vol. 11 122 family contains 28 members and can be divided into GC is generally poor due to the fact that many two groups: fibrillar collagen (type I, II, III, V, XI) and patients had metastasis by the time of diagnosis and non-fibrillar collagen (type IV, VIII, X, IX, XII, XIV, nearly inevitable tumor progression in spite of good XV, XVIII, XIX, XXI) [8]. After synthesis in the rough initial treatment response, resulting in a five-year endoplasmic reticulum, the collagen precursor went survival rate less than 10% [28]. Recently, increased through a process of modifications including lysine collagen deposition was found in GC compared to hydroxylation for maturation. Hydroxylysine residue non-neoplastic tissues, and collagen width was usually occurs in the Y position of the repeating demonstrated to be a powerful parameter for Gly-X-Y (X and Y represent proline or hydroxyl- predicting 5-year overall survival of GC [29]. Further proline) motif [9], and is critical for the formation of studies of collagen related oncogenes would collagen crosslink and glycosylation, stabilizing the contribute to identifying novel molecular markers of newly formed collagen fibers and enhancing the GC progression and help develop new diagnostic and stiffness of the matrix [10]. Abnormal lysyl therapeutic strategies. With respect to PLOD family hydroxylation contributes to the progression of many genes, the expression pattern and its relation with collagen-related diseases, such as fibrosis and cancer patient prognosis of individual PLODs in GC have [11]. Therefore, modification events of maturation and recently been reported. High PLOD1 mRNA the corresponding enzymes are important for proper expression level was observed in gastric cancer and function of collagen and related cancer cell biology. was associated with poor patient prognosis [18]. Procollagen-lysine, 2-oxoglutarate 5-dioxygen- PLOD2 was reported to promote cell invasiveness and ase (PLOD) is the lysyl hydroxylase responsible for migration in GC under hypoxia condition and lead to the lysyl hydroxylation of collagen [12]. In human peritoneal dissemination and poor prognosis [30]. genome, three PLODs, namely PLOD1, 2 and 3, are PLOD3 mRNA was also up-regulated in GC tissues identified and share high homology with an overall and cell lines, which was significantly correlated with identity being 47% in protein sequences [9]. PLOD1 larger tumor size and poor patient prognosis [31]. In and PLOD3 hydroxylate lysyl residues in the collagen the current study, we investigated on the PLOD triple helix, whereas PLOD2 hydroxylates lysyl family genes based on different online databases. residues in the telopeptides of collagen. In addition, Study of individual PLOD family genes would be PLOD3 has glycosylation activity towards mono- included, and we also highlight the prognostic values saccharide or disaccharide attaching to collagen of co-expression of three PLOD genes for GC patients, hydroxylysines [13]. Mounting evidences are suggest- and prediction of PLODs-related signaling pathways ing that increased collagen deposition and cross- in GC. linking can promote cancer progression by enhancing cancer cell proliferation, invasion and migration Methods [14-17], and overexpression of each PLOD gene has ONCOMINE analysis been reported in different human cancers. Higher expression level of PLOD1 was observed in The mRNA expressions of PLOD1, 2 and 3 in gastrointestinal carcinoma, and was related to different cancer types were analyzed using carcinogenesis and clinical outcome [18]. PLOD1 and ONCOMINE gene expression array database [32] PLOD2 were reportedly up-regulated in esophageal (www.oncomine.org), which is publicly accessible squamous-cell carcinoma and were associated with online. In this study, Student’s t test was used to cell apoptosis, cell cycle and metastasis regulation [19, generate the P-value for expression differences of 20]. Additionally, PLOD2 has been shown to promote PLOD family genes between cancer specimens and metastasis in lung cancer, renal cell carcinoma and normal controls. The fold change was defined as > 1 glioma [21-23]. Overexpression of PLOD3 was and the P-value was set up at 0.05. observed in lung cancer and glioma, and was also UALCAN analysis associated with tumor progression and metastasis [24, UALCAN [33] is an interactive web-portal to 25]. RNA interfering mediated down regulation of perform in-depth analyses of level 3 RNA-seq and PLOD2 expression inhibited proliferation and clinical data from 31 cancer type of The Cancer invasion of bladder cancer cells [26]. Therefore, the Genome Atlas (TCGA) database and it is publicly PLOD family genes could not only act as prognostic available at http://ualcan.path.uab.edu. We use it to signatures but also become druggable targets for analyze the relative expression of PLOD family genes cancer therapy. between GC and normal samples and to compare the Gastric cancer (GC) is the fifth leading expression differences among various sub-groups malignancy and the third leading cause of based on cancer stage. The expression level was cancer-related deaths globally [27]. The prognosis of normalized as transcripts