{'fraction_non_alphanumeric': 0.029726977353183647, 'fraction_numerical': 0.019058241598370593, 'mean_word_length': 5.310281517747858, 'pattern_counts': {'":': 0, '<': 6, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Identification of Critical Pathways and Hub Genes in LanCL1-Overexpressed Prostate Cancer Cells

Background Prostate cancer is one of the most common malignancies in urology, especially in developed countries. Our previous studies showed that Lanthionine synthase C-like protein 1 (LanCL1) can promote the proliferation of prostate cancer cells and protect cells from oxidative stress. Also, LanCL1 protects cells by inhibiting the JNK signaling pathway after H2O2 treatment. Materials and Methods In our study, we analyzed the data of RNA-seq to identify the DEGs after LanCL1 overexpression. We performed a functional enrichment analysis with gene set enrichment analysis (GSEA) and a database for annotation, visualization, and integrated discovery (DAVID). We also identified the critical hub gene correlated with disease prognosis by Cox regression analysis. Results A total of 8928 DEGs were identified. Through the analysis of GO and KEGG, we found that DEGs are significantly enriched in categories related to metabolism, cancer-related signaling pathways, and inflammation. The top 15 hub genes were then identified and ranked by degree from the protein–protein interaction network. Survival analysis showed 4 hub genes related to disease prognosis and ICAM1 expression is an independent risk factor for the prognosis. Conclusion Our results suggest the critical genes and pathways that might play key roles after LanCL1 overexpression in prostate cancer. We also provide candidate gene targets that might play important roles in prostate cancer development.


Introduction
Prostate cancer remains the most common cause of cancer-related death in men, particularly in developed countries. 1 Up to now, many studies have found out multiple critical genes, cellular processes, and signaling pathways that contribute to disease progression. In previous study, we focused on LanCL1 and explored whether it contributes to prostate cancer disease progression. As a result, we determined that LanCL1 expression correlated with the progression of human prostate cancer and that LanCL1 could also promote prostate cancer proliferation. Moreover, LanCL1 protects prostate cancer cells from reactive oxygen species (ROS) via suppression of JNK pathway. 2 Lanthionine synthase C-like protein 1 (also known as P40 or GRP69A) is a mammalian member of the LanC-like protein superfamily. It contains highly differentiated peptide-modifying enzymes found in bacteria and plants (LanCs). 3 Multiple studies reveal the role of LanCL1 in relieving oxidative stress in neurons during normal development by binding zinc ion and GSH. In addition, LanCL1 can catalyze the formation of thioether products in neurons, thereby protecting them from oxidative stress. [4][5][6] In breast cancer, LanCL1 expression is related to breast cancer survival, making it a marker of aging. 7 In prostate cancer, we have preliminarily explored the role played by LanCL1 in the disease process. However, further mechanisms are still unclear. As we all know, a gene or protein usually exerts its function by affecting the signaling pathways downstream. In this study, we attempt to further study and explore the mechanism underlying to discover new potential targets for disease diagnosis and treatment.

Cell Lines and Cell Culture
LNCaP and PC-3 cells (ATCC) were maintained in RPMI1640 supplemented with 10% FBS. LNCaP and PC-3 cells were co-transfected with pPB-CAG-LanCL1 and pCMVPBase to obtain control and LanCL1 stable overexpressing LNCaP/PC-3 cells. We used 2 μg/mL puromycin (Amresco) to screen for 2 weeks for selecting the stable cell lines. Stable cell lines were identified by Western blotting. All cells were added with antibioticantifungal solution (100 units/mL penicillin, 0.1 mg/mL streptomycin, and 0.25 mg/mL amphotericin B), and under standard cell culture conditions (5% CO2, 95% humidity) grow at 37°C.
Library Preparation, and Illumina Hiseq4000 Sequencing RNA-seq was performed according to the manufacturer's guidelines (C-10,365; Life Technologies) and previous research. 9 RNA-seq transcriptome library was prepared following TruSeq TM RNA sample preparation Kit from Illumina (San Diego, CA) using 5μg of total RNA. Libraries were size selected for cDNA target fragments of 200-300 bp on 2% Low Range Ultra Agarose followed by PCR amplified using Phusion DNA polymerase (NEB) for 15 PCR cycles. After quantified by TBS380, pairedend RNA-seq sequencing library was sequenced with the Illumina HiSeq 4000 (2 × 150 bp read length). RNA-seq was performed on three biological replicates.

Analysis of RNA-Seq and Differentially Expressed Genes (DEG)
RNA-seq was performed according to the manufacturer's guidelines (C-10,365; Life Technologies) and previous research. 9 We then used EdgeR, an R package to examine differential expressed genes as previously described. 10,11 DEG was determined using the following criteria: | Fold Change (FC) | ≥ 2 and P value < 0.05. DEG is used for further bioinformatics analysis.
Gene Set Enrichment Analysis (GSEA) GSEA v3.0 was used to analyze the mRNA expression level of biological function annotations and pathways after LanCL1 overexpression. The results indicate the effects of LanCL1 overexpression on various biological functional gene sets. The number of permutations is set to 100. In the analysis, we used the following values as cut-off values: the P-value cut-off value is 0.05 and the false discovery rate (FDR q-val) is <0.25.

Functional Annotation and Pathway Enrichment Analysis
We uploaded the top 2000 DEGs and the top 2000 downgraded DEGs to the database, respectively, for annotation, visualization, and integrated discovery (DAVID) website for GO analysis and KEGG analysis, 12 while specifying a P-value <0.05 for statistical significance.

Integration of Protein-Protein Interaction (PPI) Networks and Module Analysis
We used STRING database, which provides a critical assessment and integration of PPI, to analyze the protein interaction. 13 To evaluate the interactive relationships among the DEGs above, we uploaded the top 1000 upregulated and top 1000 downregulated DEGs together to STRING. Then, only experimentally validated interactions (total score > 0.4) were selected as meaningful. The PPI network was then used for module screening by Molecular Complex Detection (MCODE) in Cytoscape (scores > 3, nodes > 4). 14

Publicly Available Gene Expression and Clinical Datasets for Prognosis Analysis
We download the gene expression data with corresponding prognostic profiles of 140 prostate cancer patients from online dataset GSE21032 for prognosis analysis. 15 The Kaplan-Meier method using Log-rank test was used to calculate the clinical results between different gene expression groups by Graphpad. Values of P <0.05 were considered statistically significant.

Statistical Analysis
The Kaplan-Meier method using Log-rank test was used to calculate the clinical results by Graphpad. Multivariate Cox regression analyses were used to investigate the association between prognosis and hub genes expression/clinical factors. The FDR in GSEA was tested several times using the Benjamini-Hochberg program to control the FDR. 16,17 Values of P <0.05 were considered statistically significant. All statistical analyses were performed using SPSS v19, Graphpad, and R 3.3.0.

RNA-Seq and Identification of DEGs
To further explore the critical signaling pathways or the regulatory network downstream, we obtained the LanCL1 stably overexpressed LNCaP cell line. Western blotting showed upregulation of LanCL1 in stable cell lines ( Figure 1A). We then performed RNA-seq on a stable cell line with LanCL1 overexpression or an empty vector as a control. RNA-seq was performed on three biological replicates.
Based on computer analysis, we used P < 0.05 and Fold Control (FC)≥ 2.0 or ≤0.5 as criteria, and identified 8928 genes as DEGs, of which 4755 were upregulated and 4173 were downregulated. The volcano map of DEGs is shown in Figure 1B, while the heat map of DEGs expression for the first 100 genes is shown in Figure 1C.

GO and KEGG Analyses of DEGs
To further analyze the DEGs at the functional level, we submitted the top 2000 upregulated and top 2000 downregulated DEGs, respectively, online for further analyses using DAVID. In the upregulation group, the GO-BP (Biological Process) analysis of DEGs ( Figure 1D) suggested significant enrichment in extracellular matrix organization, collagen fibril organization, proximal/distal pattern formation, cell adhesion, peptidyl-tyrosine phosphorylation, axon guidance, anterior/posterior pattern specification, chemical synaptic transmission, regulation of ion transmembrane transport, ureteric bud morphogenesis, angiogenesis, negative regulation of neuronal apoptosis, ureteric bud development, calcium ion-regulated exocytosis of neurotransmitter, and positive regulation of synapse assembly. We also showed the results of GO-CC (Cellular Component) and GO-MF (Molecular Function) in Supplementary Figures 1A and B. While KEGG pathway analysis showed enrichment mainly in focal adhesion, ECM-receptor interaction, PI3K-AKT signaling pathway, arrhythmogenic right ventricular cardiomyopathy (ARVC), axon guidance, pathways in cancer, Ras signaling pathway, cell adhesion molecules (CAMs), amoebiasis, cAMP signaling pathway, hypertrophic cardiomyopathy (HCM), protein digestion and absorption, proteoglycans in cancer, platelet activation, and Rap1 signaling pathway ( Figure 1E).
In downregulation group, in GO-BP (Biological Process) analysis, DEGs were significantly enriched in oxidation-reduction process, homophilic cell adhesion via plasma membrane adhesion molecules, inflammatory response, establishment of skin barrier, nucleosome assembly, phosphatidylserine acyl-chain remodeling, acute-phase response, angiogenesis, phosphatidylinositol acyl-chain remodeling, positive regulation of cell adhesion, phosphatidylethanolamine acyl-chain remodeling, phospholipid metabolic process, negative regulation of viral genome replication, positive regulation of inflammatory response, and transmembrane transport ( Figure 1F). The results of GO-CC (Cellular Component) and GO-MF (Molecular Function) were shown in Supplementary Figures 1C and D. KEGG pathway analysis showed enrichment mainly in arachidonic acid metabolism, alcoholism, systemic lupus erythematosus, linoleic acid metabolism, alpha-Linolenic acid metabolism, drug metabolism -cytochrome P450, chemical carcinogenesis, influenza A, glutathione metabolism, mineral absorption, ether lipid metabolism, transcriptional misregulation in cancer, Aldosterone-regulated sodium reabsorption, Hepatitis C, and Bladder cancer ( Figure 1G).

GSEA
To determine the effect of LanCL1 overexpression on prostate cancer progression, we then analyzed various

PPI Network Analysis
To investigate the protein-protein interaction and find out the hub genes, we then uploaded the top 1000 upregulated DEGs and 1000 downregulated DEGs to the   (Figure 3). Functional annotations in module 1 genes are richer in G protein-coupled receptor signaling pathways, chemotaxis, positive regulation of cytosolic calcium ion concentration, chemokine signaling pathways, neuroactive ligand-receptor interactions, and Cytokine-cytokine receptor interaction. Enrichment analysis showed that the genes in module 2 were mainly associated with interferon-gamma-mediated signaling pathway, type I interferon signaling pathway, collagen catabolic process, protein digestion and absorption, ECM-receptor interaction, and herpes simplex

Dovepress
Tang et al

Clinical Outcome and Hub Genes
To select the hub genes that are closely related to clinical prognosis, we use the data and clinical information of GSE21032. 15 Figure 4 and Table 2 show that 4 of the 15 hub genes strongly associated with clinical outcomes, including APOE, FGF2, ICAM1, and KIT. Lower APOE expression level correlated with better outcome, while higher FGF2, ICAM1, and KIT correlated with better outcome. Then, we reorganized the clinical data of GSE21032 and found that it contains two clinical factors, Gleason score and Tumor stage (Table 3). To investigate the association between prognosis and these clinical two factors, and to evaluate whether the expressions of the 4 hub genes above are independently associated with disease recurrence, we first did Kaplan-Meier analysis. The results suggest that both Gleason score and Tumor stage significantly correlated with the prognosis of the disease (Table 2). Then, we conducted a multivariate COX regression analysis. Results showed that only the expression of ICAM1 in the four genes, as well as Gleason score and T stage, is an independent risk factor for the prognosis (Table 4). At last, to confirm our results, we investigated the expression of the 4 critical hub genes after LanCL1 overexpression in another prostate cancer cell line, PC-3, which revealed the same results as in LNCaP cells ( Figure 5). These

DovePress
OncoTargets and Therapy 2020:13 results suggested that these 4 genes, especially ICAM1, might be potential key genes for further research.

Discussion
Prostate cancer is still one of the most common malignant tumors in male urinary system all around the world, especially in developed countries. In-depth research on the downstream pathways and regulatory networks of key genes involved in disease progression is of great significance for understanding the disease and finding new therapeutic targets. In our previous study, we revealed the role of LanCL1 as an oncogene in prostate cancer to promote cell proliferation and protect cell from attack of ROS via JNK pathway. However, important downstream genes and regulatory signaling pathways are still not clear.
Here in this study, we overexpressed LanCL1 gene in LNCaP cells and did RNA-Seq analysis. By selecting DEGs and performing GO and KEGG analysis on the

DovePress
top 4000 DEGs, we found that LanCL1 overexpression was significantly associated with the metabolism, axon and synaptic function, and several signaling targets or pathways, such as cAMP signaling pathway, Ras signaling pathway, PI3K-AKT-mTOR signaling pathway, and Rap1 signaling pathway. Previous researches have shown that Rap1 (ras-related protein 1) signaling pathway and Ras signaling pathway are related to JNK and participate in cellular oxidative stress and the production and elimination of ROS in a variety of cellular backgrounds. [18][19][20][21] In addition, the cAMP signaling pathway is also involved in the production, metabolism, and clearance of ROS in the body, which has the same important role as the JNKrelated pathway. 22,23 These results suggest that LanCL1 plays an important role in assisting cells in processing ROS, which is consistent with our previous findings. Besides, because the oxidative stress state and related signaling pathways in the cells will change after H 2 O 2 treatment, there may be significant differences in the signaling pathways between treated and untreated cells. This might be the main reason why JNK pathway was not enriched in the RNA-seq results of untreated cells in this study. JNK pathway is the main mechanism by which LanCL1 protects prostate cancer cells from oxidate stress in H 2 O 2 treated cells. Our results here also show that the mechanism of LanCL1 in dealing with ROS is more complicated. It also validates our previous experiments, that is, changing the activity of JNK pathway in rescue experiments only partially restored the ROS-related phenotype. In addition to oxidative stress such as ROS, we also found that many pathways related to tumorigenesis and progression appeared enriched after LanCL1 overexpression, mainly including WNT beta-catenin signaling, PI3K-AKT signaling pathway, and KRAS signaling. These pathways play important roles in the pathogenesis and progression of various tumors, including prostate cancer.
In prostate cancer, WNT/β-catenin and AR signals inhibit each other. β-catenin activation and PTEN deletion synergistically drive non-AR independent castrationresistant prostate cancer. [24][25][26] In most cancer types, the PI3K signaling pathway is activated and promotes

DovePress
OncoTargets and Therapy 2020:13 tumorigenesis by regulating nutrient metabolism, cell proliferation, survival, migration, and angiogenesis. 26 The potential mechanism of PI3K/AKT activation is mainly due to the deletion or mutation of its key negative regulatory gene, PTEN. 27 In fact, it is known that the mutual interference between PI3K/AKT and other pathways (such as AR, WNT, and ERK signaling pathways) plays a vital role in the disease progression and drug resistance of prostate cancer, making co-targeting the PI3K/AKT signaling pathway and its synergistic pathways critical to improve the anticancer efficacy of PI3K/AKT inhibitors. 28,29 Our results indicate that LanCL1 may contribute to prostate cancer by influencing these important tumor regulatory signal pathways or signal networks, especially its promotion effect on prostate cancer cell proliferation. LanCL1 also may affect one of these pathways. As for the specific regulatory mechanism, further molecular biology experiments are needed to verify it.
In the PPI network analysis, we identified the top 15 central genes with the highest degree of interaction. TNF (tumor necrosis factor) is the top hub gene, which is mainly secreted by macrophages. TNF can bind to its receptors TNFRSF1A/TNFR1 and TNFRSF1B/TNFBR to exert its function. TNF could involve in the development of tumors by regulating inflammatory and immune responses. Among the 15 hub genes, we also identified 4 hub genes that strongly associated with clinical outcomes in prostate cancer, including APOE, FGF2, ICAM1, and KIT. We used data from a classic and authoritative prostate cancer RNA-seq database GSE21032 to analyze. These 4 genes might play critical roles in the development of prostate cancer and are closely related to the prognosis of the disease. Then, we conducted a multivariate COX regression analysis and found that only the expression of ICAM1 is an independent risk factor for the prognosis. This result emphasizes the important potential role of ICAM1 in prostate cancer initiation and development. ICAM1 (Intercellular adhesion molecule 1) is a cell surface glycoprotein expressed on endothelial cells and immune system cells. It binds to CD11a/CD18 or CD11b/CD18 type integrins and is also used by rhinovirus as a receptor. 30,31 Previous studies have shown that ICAM1 may facilitate breast and lung. 32-34 Researchers

DovePress
reported that CCN3-ICAM1 and MiRNA-296-3p-ICAM1 axes could influence cell migration in prostate cancer cells. 35,36 However, the role and mechanism of ICAM1 in prostate cancer are still largely unknown. ICAM1 may be further studied as key a gene to find more clues about the initiation and progression of prostate cancer disease. As for APOE, FGF2, and KIT, many previous studies have investigated their roles in prostate cancer, and our results further confirmed these findings.
Our study contained one limitation. In our study, we analyzed the data of RNA-seq to identify the DEGs after LanCL1 overexpression and figured out the critical genes and pathways that might play key roles after LanCL1 overexpression in prostate cancer using bioinformatics analysis. The mechanism and validation of LanCL1 overexpression in prostate cancer still need further research in clinical and molecular biology experiments.

Conclusions
In conclusion, our results revealed the critical genes and pathways that might play key roles after LanCL1 overexpression in prostate cancer. We also provide 4 critical hub genes that could play critical roles in the development of prostate cancer and are closely related to the prognosis of the disease,

DovePress
OncoTargets and Therapy 2020:13 which may help in developing prostate cancer treatment strategies.

Publish your work in this journal
OncoTargets and Therapy is an international, peer-reviewed, open access journal focusing on the pathological basis of all cancers, potential targets for therapy and treatment protocols employed to improve the management of cancer patients. The journal also focuses on the impact of management programs and new therapeutic agents and protocols on patient perspectives such as quality of life, adherence and satisfaction. The manuscript management system is completely online and includes a very quick and fair peer-review system, which is all easy to use. Visit http://www.dovepress.com/ testimonials.php to read real quotes from published authors.