Li Wang1†
Guodong Song- 1Department of Gastrointestinal Surgery, The Second Hospital of Tianjin Medical University, Tianjin, China
- 2Department of Respiratory, The Second Hospital of Tianjin Medical University, Tianjin, China
- 3Radiology Department, The Second Hospital of Tianjin Medical University, Tianjin, China
Background: Gastric cancer (GC) remains the fifth most commonly diagnosed malignancy worldwide, with a poor prognosis. The lysyl oxidase (LOX) family, a type of secreted copper-dependent amine oxidases, is comprised of LOX and four LOX-like (LOXL) 1–4 isoforms and has been reported to be dysregulated in a number of different type cancers. However, the diverse expression patterns and prognostic values of LOX family in GC have yet to be systematically analyzed.
Methods: ONCOMINE, GEPIA, UALCAN, Kaplan–Meier Plotter, LOGpc, cBioPortal, GeneMANIA and Metascape databases were utilized in this study to analyze the expression, prognostic values, mutations and functional networks of LOX family in GC.
Results: The mRNA expression levels of LOX, LOXL1 and LOXL2 were significantly higher in GC, the expression level of LOXL3 was contrary in different databases, while the expression level of LOXL4 made no difference; the expression levels of LOX, LOXL1 and LOXL3 were higher in stages 2–4 than that of normal tissues and stage 1, while the mRNA level of LOXL2 in stage 1–4 was higher than normal tissues; patients with high expression of LOX and LOXL 2-4 had poor OS; the genes correlated with LOX and LOXL2 were enriched in extracellular matrix organization, vasculature development and skeletal system development.
Conclusion: Our results indicated that the LOX family, especially LOX and LOXL2, might play an important role in GC oncogenesis, and they may become biomarkers for predicting tumor prognosis and potential targets for tumor therapy.
Introduction
Even though declines in (GC)incidence and mortality rates have been observed consistently across world regions, GC remains the fifth most commonly diagnosed malignancy worldwide, with over 1 million estimated new cases in 2018 (Bray et al., 2018). However, due to its advanced-stage diagnosis, excess mortality from this cancer is high, making GC the third most common cause of cancer related death with 784,000 deaths globally (Bray et al., 2018). Despite major advances in understanding the epidemiology, pathology, and molecular mechanisms of GC and in implementing emerging therapeutic options such as targeted and immune-based therapies, not all patients respond to existing molecularly targeted agents developed for certain acknowledged biomarkers (Chau, 2017). As a result, it is of great importance to further identify novel diagnostic and prognostic biomarkers in terms of the biological complexity, poor prognosis and high reoccurrence of GC (Kang et al., 2018).
The extracellular matrix (ECM) is a complex network of secreted molecules, the function of which is to program cell behavior including supporting cell adhesion, survival and migration. The remodeling of the ECM in cancer plays an important role in controlling the progression of disease and influences cell growth, motility and survival (Yamauchi et al., 2018a). The lysyl oxidase (LOX) family, which are well known as ECM-modifying proteins, they participate in the crosslinking of collagens and elastin in the ECM, promoting its maturation (Yamauchi et al., 2018b; Chitty et al., 2019). The LOX family, a type of secreted copper-dependent amine oxidases, is comprised of five homologous members: LOX and lysyl oxidase-like proteins 1–4 (LOXL1, LOXL2, LOXL3 and LOXL4) (Molnar et al., 2003). Structurally, these members are all characterized by a highly conserved C-terminal domain and a variable N-terminal domain. The composition of C-terminal domain contains copper binding domain, amino acid residues forming lysine tryosylquinone (LTQ), cofactor formation, and a cytokine receptor-like (CRL) domain (Wang et al., 2016). The pro-domains are existed in the N-terminal region of LOX and LOXL1, whereas four scavenger-receptor cysteine-rich (SRCR) domains in the N-terminal are observed in LOXL2-4 (Xiao and Ge, 2012). Mature active forms of LOX and LOXL1 are obtained through a specific cleavage process induced by bone morphogenetic protein 1 (BMP-1), whereas LOXL2, LOXL3, and LOXL4 do not require this cleavage process to mature. In particular, a pre-pro-LOX protein is encoded by LOX mRNA and converted to the inactive LOX preprotein (pro-LOX) in the cytoplasm. The Pro-LOX protein is further cleaved by BMP-1 to form an active LOX with the LOX propeptide (LOX-PP) to perform its function (Wang et al., 2016). In healthy tissue, the synthesis of the LOX family is tightly regulated to control the amount of active LOX family members present. While the LOX family has been reported to be dysregulated in a number of different type cancers (Li et al., 2015; Salvador et al., 2017; Shao et al., 2019; Zeltz et al., 2019; Hu L. et al., 2020). The changes in LOX family member regulation, expression and subsequently enzymatic activity are therefore important factors in cancer progression (Setargew et al., 2021). The LOX family of enzymes may be favorable targets for anti-stromal therapeutics due to their importance in cancer development and progression when compared to healthy state ECM (Setargew et al., 2021). Additionally, highly LOX family expressing tumors have increased LOX family levels detectable in plasma (Rachman-Tzemah et al., 2017), and thus indicate the potential to be used as tumor serum markers.
The dysregulated expression level of LOX family and their relationship with clinicopathological features and prognosis have been partly reported in human GC. With the revolutionized development of microarray and bioinformatic technology, we conducted this systematical study using the data from The Cancer Genome Atlas (TCGA) and other versatile public databases to analyze the expression levels, mutations, functional networks and prognostic values of different LOX in GC, so as to reveal potential diagnostic, therapeutic, and prognostic targets for GC, and the results in different databases were verified with each other to make the results more convincible.
Materials and Methods
Oncomine Database
The mRNA expression levels of LOX family in various cancers and their normal tissue counterparts were analyzed using the Oncomine database (http://www.oncomine.org/) (Rhodes et al., 2007). A p-value of 0.001, a fold change of 2, and a gene rank in the top 10% were set as the significance thresholds. The p value was calculated using the Student’s t-test.
GEPIA Database
GEPIA (http://gepia2.cancer-pku.cn/) is a gene expression analysis web which contains 9,736 tumors and 8,587 normal samples from the TCGA and the Genotype-Tissue Expression (GTEx) project (Tang et al., 2019). Here we used GEPIA to compare the expression levels between TCGA cancer and matched TCGA normal and GTEx normal. The results were expressed as boxplots, and the cutoff criteria were set as p < 0.01 and |Log2FC| > 1.
UALCAN
UALCAN (http://ualcan.path.uab.edu/) is a comprehensive, user-friendly, and interactive web resource for analyzing cancer OMICS data (Chandrashekar et al., 2017). In this study, we used UALCAN to compare the expression levels of LOX family and their relationship with tumor stages. Student’s t-test was used to generate a p-value and the p-value cutoff was 0.05.
Survival Analysis
We used The Kaplan Meier plotter (http://kmplot.com/analysis/) and LOGpc (http://bioinfo.henu.edu.cn/DatabaseList.jsp) to evaluate the prognostic value of LOX family mRNA expression in which cancer patients were split into high and low expression group based on median values of mRNA expression and validated by K-M survival curves. The Kaplan Meier plotter is capable to assess the effect of 54 k genes (mRNA, miRNA, protein) on survival in 21 cancer types including breast, ovarian, lung and GC (Szász et al., 2016). LOGpc (Long-term Outcome and Gene Expression Profiling Database of pan-cancers) encompasses 209 expression datasets, provides 13 types of survival terms for 31,310 patients of 27 distinct malignancies (Liu et al., 2018). The log-rank test was used for computing p-value, with the hazard ratio (HR) and 95% confidence intervals (CI), and p < 0.05 was regarded as significant.
cBioPortal
The cBioPortal (https://www.cbioportal.org/) for Cancer Genomics provides visualization, analysis and download of large-scale cancer genomics data sets (Cerami et al., 2012). The stomach adenocarcinoma (TCGA, Firehose Legacy) dataset was selected to figure out the alterations of LOX family. We also estimated the correlations of each LOX family by analyzing their mRNA expression (RNA Seq V2 RSEM), and then the Spearman correlation coefficient was put into Microsoft Excel 2016 to draw the heat maps. Besides, genes with the highest expression correlation with each LOX protein were generated by cBioPortal, and the top 50 co-expressed genes with highest Spearman correlation score were included in the following functional enrichment analysis.
GeneMANIA Database
GeneMANIA (https://genemania.org/) is actively developed at the University of Toronto, in the Donnelly Centre for Cellular and Biomolecular Research, in the labs of Gary Bader and Quaid Morris, and it can find other genes that are related to a set of input genes, using a very large set of functional association data (Warde-Farley et al., 2010). GeneMANIA constructed protein–protein interaction (PPI) networks in terms of physical interaction, coexpression, predicted, colocalization, common pathway, genetic interaction, and shared protein domains. In this study, GeneMANIA was used to generate and analyze gene co-expression network.
Metascape
We used the Metascape web (http://metascape.org) (Zhou et al., 2019) to perform functional enrichment analysis by using the top 50 co-expressed genes of LOX family. The functional process and pathway, following the default, included Canonical Pathway (MSigDB), Hallmark Gene Sets (MSigDB), Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway and Gene Oncology (GO).
Results
Transcriptional Levels of LOX Family in GC and Other Cancers
The transcription level differences of LOX family between tumor and normal tissues were analyzed in multiple cancer types using the Oncomine database. As shown in Figure 1, the expressions of LOX, LOXL1 and LOXL2 have been up-regulated in most of the studied tumors, while LOXL3 and LOXL4 only have differences in expression in a small number of tumors.
FIGURE 1. Transcription levels of LOX family in different cancer types (Oncomine).
The mRNA expression level of LOX was significantly up-regulated in patients with GC in 4 analyses out of 23 in 2 datasets out of 7. Chen gastric statistics (Chen et al., 2003) indicated that LOX is overexpressed in gastric intestinal adenocarcinoma compared with gastric normal tissue with a fold change of 2.312, diffuse gastric adenocarcinoma with a fold change of 2.004, and gastric mixed adenocarcinoma with a fold change of 3.232 (Table 1). Wang gastric analysis (Wang et al., 2012) revealed that LOX is upregulated in GCwith a fold change of 2.287 (Table 1). The transcriptional level of LOXL1 was significantly up-regulated in GC in 3 analyses out of 23 in 3 datasets out of 7. In Chen’s (Chen et al., 2003) dataset, the expression of LOXL1 was 2.077 times higher in gastric mixed adenocarcinoma than normal tissues (Table 1). In Wang’s dataset (Wang et al., 2012), the expression of LOXL1 was 2.083 times higher in GC tissues than normal tissues (Table 1). In DErrico’s dataset (D'Errico et al., 2009), the expression of LOXL1 was 2.192 times higher in gastric mixed adenocarcinoma than normal tissues (Table 1). Upregulation of LOXL2 was observed in 4 analyses in GC tissues compared with normal tissues, with a fold change of 2.118 in Cui’s dataset (Cui et al., 2011), a fold change of 2.424 in Wang’s dataset, a fold change of 2.039 in Chen’s (Chen et al., 2003) dataset, and a fold change of 2.681 in DErrico’s dataset (D'Errico et al., 2009) respectively (Table 1). The analyses of Oncomine database showed no difference in transcriptional levels of LOXL3 and LOXL4 in GC (Figure 1).
TABLE1. Remarkable changes of LOX family expression in transcription level between GC and normal gastric tissues (ONCOMINE).
Furthermore, we used GEPIA to compare the mRNA expression of LOX family between 408 TCGA GC and 211 matched TCGA normal and GTEx normal, used UALCAN to compare the expression levels of LOX family between 415 TCGA GC and 34 TCGA normal. The GEPIA analyses showed that LOX/LOXL1/LOXL2 were higher in GC than in normal tissues (Figure 2A). The results of UALCAN indicated that LOX/LOXL1/LOXL2/LOXL3 were over-expressed in GC (Figure 2B).
FIGURE 2. Boxplot showing the expression levels of LOX family inGC. (A) GEPIA analysis. The number of normal samples are 211 (grey box), and number of primary tumor samples are 408 (red box), red star means p < 0.01; (B) UALCAN analysis. The number of normal samples are 34 (blue box), and number of primary tumor samples are 415 (red box) ***p < 0.001.
Relationship Between the mRNA Levels of LOX Family and the Cancer Stages of GC
Next, we analyzed the relationship between the mRNA expression of different LOX family members with patients’ individual cancer stages of GC patients by using UALCAN. LOX, LOXL1, LOXL2 and LOXL3 groups significantly varied, whereas LOXL4 groups did not significantly differ (Figure 3). According to clinical stages, the mRNA level of LOX in stage 2–4 was higher than normal tissues and stage 1, while there is no difference between stage 2–4 (Figure 3A). The similar result was found in expression of LOXL1 and LOXL3 (Figures 3B,D). The mRNA level of LOXL2 in stage 1–4 was higher than normal tissues and highest expression was found in stage 2 tissues (Figure 3C). There was no difference in expression of LOXL4 among different stages (Figure 3E).
FIGURE 3. Transcription levels of LOX family in different stage of patients with GC. LOX (A), LOXL1 (B), LOXL2 (C), and LOXL3 (D) groups significantly varied, whereas LOXL4 (E) groups did not significantly differ. (UALCAN, ***p < 0.001, **p < 0.01, *p < 0.05).
Prognostic Values of LOX Family in GC
We further explored the prognostic values of LOX family in patients with GC by using the Kaplan–Meier Plotter database and LOGpc database. We separated all GC patients into two groups (high vs. low) based on median expression values for each LOX protein across all GC samples and compared overall survival (OS) between the two groups. The Kaplan-Meier curve and log rank test analyses revealed that all of LOX family members were significantly associated with the OS (p < 0.05) in patients with GC (Figure 4A). Meanwhile, the results of LOGpc analyses indicated that the increased LOX, LOXL2-4 mRNA expression were associated with low OS (p < 0.05) in patients with GC, but the expression of LOXL1 had no correction with prognosis (p > 0.05) in patients with GC (Figure 4B).
FIGURE 4. Prognostic feature of mRNA expression of distinct LOX family members in GC patients. (A) OS of Kaplan–Meier plotter revealed that all of LOX family members were significantly associated with the OS (p < 0.05); (B) OS of LOGpc indicated that LOX, LOXL2-4 mRNA expression were associated with OS (p < 0.05).
Genetic Mutations and PPI Network of LOX Family
We analyzed the types and frequency of LOX Family alterations in a cohort of GC patients using cBioPortal. The LOX family were altered in 55 (14%) samples of 393 patients with stomach adenocarcinoma (Figure 5A). We also calculated the correlations of LOX family with each other by analyzing their mRNA expressions (RNA sequencing (RNA seq V2 RSEM)) via the cBioPortal online tool for stomach adenocarcinoma (TCGA, Firehose Legacy) and Pearson’s correction was included. The results indicated significant and positive correlations in the following pairs: LOX and LOXL2, LOX and LOXL3, LOXL1 and LOXL4 (Figure 5B). The mutation rates of LOX, LOXL1, LOXL2, LOXL3 and LOXL4 were 4, 2.5, 5, 2.5 and 2.8%, respectively (Figure 5C).
FIGURE 5. Genomic alterations (cBioPortal) and network (GeneMANIA) of LOX family inGC. (A) Distribution of LOX family genomic alterations according to cancer type. (B) Correlations of different LOX family members with each other inGC. (C) OncoPrint of LOX family alterations inGC. (D) Network of the 20 most frequently altered neighboring genes of LOX family.
Next, we used Gene-MANIA to construct a PPI network for LOX family, and the result is shown in Figure 5D. The most top 20 related genes are as follows: TLL1, BMP1, EGFL7, EFEMP2, SNAI1, COL1A2, COL1A1, MFAP4, FBLN5, COL3A1, ELN, JTB, DHDDS, SEPHS2, CD5, SCART1, CD163L1, CD163, SSC4D and SSC5D.
Functional Enrichment Analysis of Genes Co-expressed With LOX/LOXL2
Considering the expression level of LOX family in GC tumor tissues and their prognostic values in GC, LOX and LOXL2 were taken into next functional enrichment analysis. The top 50 genes, which had the most significant correlation with LOX/LOXL2 generated by cBioPortal, were included in the following functional enrichment analysis using Metascape.
The results shown in Figure 6 indicated that the most top 5 significant biological process with LOX and its co-expressed genes were GO: 0030198 (extracellular matrix organization), M5884 (NABA CORE MATRISOME), M18 (PID INTEGRINL PATHWAY), GP: 0001944 (vasculature development), and GO: 0001501 (skeletal system development) (Figure 6A); while the most top 5 with LOXL2 were GO: 0030198 (extracellular matrix organization), R-HAS-1474290 (Collagen formation), GO: 0035987 (endodermal cell differentiation), GP: 0001944 (vasculature development), and GO: 0001501 (skeletal system development) (Figure 6C). Figures 6B,D were networks that exhibited the interactions among cluster of genes enriched in biological processes and pathways mentioned above.
FIGURE 6. Functional enrichment analysis of genes co-expressed with LOX and LOXL2 using Metascape. (A,C) Heatmaps of the molecular functions, biological processes, or pathways enriched with LOX and LOXL2 co-expressed genes. The bar color shade was decided by the p value, the deeper the shade, the less the p value. (B,D) Networks exhibiting interactions among the clusters of genes enriched in the molecular functions, biological processes, or pathways presented in the heatmaps. The points in different colors represented clusters of genes enriched in different molecular functions, biological processes, and pathways.
Discussion
The LOX family has been reported to be dysregulated in a number of cancers (Li et al., 2015; Salvador et al., 2017; Shao et al., 2019; Zeltz et al., 2019; Hu L. et al., 2020). Although the role of LOX family in tumorigenesis and prognosis of several cancers has been partially confirmed, further bioinformatics analysis of GChas yet to be performed. In this study, we used multitalented public databases to reveal the dysregulated expression of the LOX family and their relations with tumor stage and prognosis. We mainly found that the mRNA expression levels of LOX, LOXL1 and LOXL2 were significantly higher in GC, the expression level of LOXL3 was contrary in different databases, while the expression level of LOXL4 made no difference; the expression levels of LOX, LOXL1 and LOXL3 were higher in stages 2–4 than that of normal tissues and stage 1, while the mRNA level of LOXL2 in stage 1–4 was higher than normal tissues; patients with high expression of LOX and LOXL 2-4 had poor OS; the genes correlated with LOXL2/4 were enriched in extracellular matrix organization, vasculature development and skeletal system development.
LOX is a secreted extracellular matrix protein that plays an important role in remodeling the extracellular matrix and promoting tumor progression. Higher LOX mRNA expression was detected in GC tissues than that in adjacent normal gastric tissues, and was significantly correlated with the invasion depth, tumor differentiation, lymph node metastasis, lymphatic invasion, venous invasion, and peritoneal metastasis in GC patients, predicting a poor prognosis of GC patients with high expression of LOX (Zhang et al., 2013). Besides, the mRNA and protein levels of LOX in GC cells and tissues were higher under hypoxia condition than that under normoxia condition (Kasashima et al., 2016). The number of migrating and invading GC cells in hypoxia was significantly decreased after knockdown of LOX (Kasashima et al., 2016). The mechanism involved in LOX-mediated proliferation facilitation in GC is enhancement of Warburg effect through regulation of HIF-1α and c-Myc (Li et al., 2019). The risk of macrophages high density, high microvessel density (MVD), low neomicrovessel maturation, MMP-9 expression and low type IV collagen was elevated after LOX overexpression, suggesting that LOX activated cancer stromal cells and facilitated the progression of GC (Peng et al., 2018). Combining LOX with CEA, CA724, CA199, and CA125 could increase the sensitivity of predicting lymph nodes metastasis and peritoneal metastasis in GC (Lai et al., 2014). Similar results have also been confirmed in our research, however, one study indicated that LOX expression was downregulated in GC, and LOX functioned as a tumor suppressor (Kaneda et al., 2004). Therefore, LOX function in GC needs to be further explored.
Relatively, few data are available on the role of LOXL1 in tumorigenesis. LOXL1 was overexpressed in GC cells, and high LOXL1 expression is a poor prognostic factor in GC patients (Kasashima et al., 2018). Moreover, LOXL1 is associated with peritoneal dissemination, potentially via promoting EMT in GC cells, and high LOXL1 expression was associated with poorly differentiated histological type, lymph node metastasis, and poor prognosis in GC (Hu Q. et al., 2020). Our study also revealed that LOXL1 is highly expressed in GC and may be related to the prognosis, although the results of the two survival databases are inconsistent.
It is first reported that LOXL2 promotes tumor progression and is associated with poor prognosis in breast cancer (Akiri et al., 2003). LOXL2 was overexpressed in GC versus normal tissues, and overexpression of LOXL2 was associated with depth of tumor invasion, lymph node metastasis and poorer overall survival (Peng et al., 2009). Furthermore, secreted LOXL2 promotes GC metastasis via Src kinase/Focal adhesion kinase (Src/FAK) pathway (Peng et al., 2009). LOXL2 expression in stromal cells was significantly associated with tumor invasion depth, lymph node metastasis, lymphatic invasion, venous invasion, peritoneal dissemination, and survival in GC patients (Kasashima et al., 2014). In our report, we illustrated that the expression of LOXL2 in GC tissues was higher than that in normal tissues, and this expression was markedly correlated with tumor stage and poor OS in patients with GC, which was consistent with reports above.
Even though LOXL3 expression was also detected in some kinds of tumors, studies have been conducted on LOXL3 were fewer. The expression of LOXL3 was detected mainly in the nucleus, and the expression of LOXL3 was correlated with tumor invasion, lymph node metastasis, and poorer prognosis of patients (Kasashima et al., 2018). Additionally, TGF-induced LOXL3 upregulation in GC cells, suggesting that LOXL3 was downstream from the TGF-signaling pathway (Kasashima et al., 2018). Our research showed that high expression of LOXL3 was confirmed in GC of TCGA data by using UALCAN database, and the expression was correlated with tumor stage, while there was no difference in the results of Oncomine and GEPIA. The survival analysis also verified that high expression was related to poor prognosis, we speculate that the prognosis of tumors are related to a variety of factors and this may be related to the target of certain drugs in the treatment of gastric cancer. Recently, more studies about LOXL3 have been published, the roles attributed to LOXL3 should be further determined.
LOXL4 was significantly up-regulated in gastric carcinoma tissues, and this over-expression is significantly correlated with tumor size, depth of tumor invasion, lymph node metastasis, higher TNM stages and poor prognosis (Li et al., 2015). LOXL4 may promote proliferation and metastasis via regulate FAK/Src pathway in GC cells (Li et al., 2015). The results about LOXL4 conducted by Kasashima et al. were similar with LOXL3 (Kasashima et al., 2018). In contrast, comprehensive bioinformatics analysis of multiple databases in our study did not yield positive results in expression of LOXL4. In other aspects, survival analysis showed that LOXL4 is associated with poor prognosis, suggesting that LOXL4 is implicated in the progression of GC. This prognostic-related reason may be the same as LOXL3.
In the functional enrichment analysis of genes co-expressed with LOX/LOXL2, the most significant biological process of LOX/LOXL2 and their co-expressed genes is ECM. ECM plays a key role in the occurrence and metastasis of gastric cancer. The destruction of the tightly coordinated ECM tissue will damage the structure and function of the gastric tissue, eventually leading to the progression of gastric cancer (Moreira et al., 2020). We speculate that LOX and LOXL2 affect the occurrence and development of gastric cancer by participating in the regulation of extracellular matrix, but further research is still needed.
Due to the secreted nature of the LOX family members, their detectable presence in the blood, and the well-established correlation between LOX family enzyme expression and prognosis in many cancers, the LOX family offers promise as an inexpensive and non-invasive companion biomarker for cancers (Setargew et al., 2021). The LOX family of enzymes are favorable targets for anti-stromal therapeutics because of their importance in cancer development and progression. A number of studies have examined the use of LOX family inhibitors in cancer therapy (Jiang et al., 2020; Smithen et al., 2020). CCT365623 is a LOX inhibitor based on methylaminothiophene. It has shown that its inhibitory effect can lead to delayed tumor development and reduced lung metastasis in mouse breast cancer models. But it has not yet been tested in a clinical (Smithen et al., 2020). XS-5382A, an oral LOXL2 inhibitor, has been shown to slow tumor growth and reduce collagen accumulation in LY2 oral cancer models and is currently being investigated in Phase 1 clinical trials in healthy adults (Clinical trial identifier: NCT04183517) (Mahjour et al., 2019). Although no inhibitors of the LOX family have currently been approved for routine clinical practice, the developing LOX family inhibitors have shown high specificity and low toxicity.
However, there are limitations in our research. Bioinformatics analysis alone cannot determine the specific mechanism of LOX family in GC. The role of the LOX family in GC might be complex, and more clinical studies and in-depth experiments are needed to verify the diagnostic value of these LOX family and explore the potential mechanism of LOX family affecting the development of GC.
Conclusion
In this study, through systematically analyzing the expression and prognostic value of LOX family in GC, we indicated that the LOX family, especially LOX and LOXL2, might play an important role in GC oncogenesis, and they may become biomarkers for predicting tumor prognosis and potential targets for tumor therapy.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author Contributions
LW and GS developed the idea, designed the research, and performed the data analysis work. SC and RZ reorganized the data and edited the manuscript and GS and YZ reviewed the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by National Natural Science Foundation of China (Grants 91959114, 81872106), Natural Science Foundation of Tianjin (19JCZDJC33900).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Akiri, G., Sabo, E., Dafni, H., Vadasz, Z., Kartvelishvily, Y., Gan, N., et al. (2003). Lysyl Oxidase-Related Protein-1 Promotes Tumor Fibrosis and Tumor Progression In Vivo. Cancer Res. 63 (7), 1657–1666.
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA: A Cancer J. Clinicians 68 (6), 394–424. doi:10.3322/caac.21492
Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data: Figure 1. Cancer Discov. 2 (5), 401–404. doi:10.1158/2159-8290.CD-12-0095
Chandrashekar, D. S., Bashel, B., Balasubramanya, S. A. H., Creighton, C. J., Ponce-Rodriguez, I., Chakravarthi, B. V. S. K., et al. (2017). UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia 19 (8), 649–658. doi:10.1016/j.neo.2017.05.002
Chau, I. (2017). Checkpoint Inhibition: an ATTRACTION in Advanced Gastric Cancer? The Lancet 390 (10111), 2418–2419. doi:10.1016/S0140-6736(17)32131-1
Chen, X., Leung, S. Y., Yuen, S. T., Chu, K.-M., Ji, J., Li, R., et al. (2003). Variation in Gene Expression Patterns in Human Gastric Cancers. MBoC 14 (8), 3208–3215. doi:10.1091/mbc.e02-12-0833
Chitty, J. L., Setargew, Y. F. I., and Cox, T. R. (2019). Targeting the Lysyl Oxidases in Tumour Desmoplasia. Biochem. Soc. Trans. 47 (6), 1661–1678. doi:10.1042/BST20190098
Cui, J., Chen, Y., Chou, W.-C., Sun, L., Chen, L., Suo, J., et al. (2011). An Integrated Transcriptomic and Computational Analysis for Biomarker Identification in Gastric Cancer. Nucleic Acids Res. 39 (4), 1197–1207. doi:10.1093/nar/gkq960
D’Errico, M., Rinaldis, E. d., Blasi, M. F., Viti, V., Falchetti, M., Calcagnile, A., et al. (2009). Genome-wide Expression Profile of Sporadic Gastric Cancers with Microsatellite Instability. Eur. J. Cancer 45 (3), 461–469. doi:10.1016/j.ejca.2008.10.032
Hu, L., Wang, J., Wang, Y., Wu, L., Wu, C., Mao, B., et al. (2020). LOXL1 Modulates the Malignant Progression of Colorectal Cancer by Inhibiting the Transcriptional Activity of YAP. Cell Commun Signal 18 (1), 148. doi:10.1186/s12964-020-00639-1
Hu, Q., Masuda, T., Kuramitsu, S., Tobo, T., Sato, K., Kidogami, S., et al. (2020). Potential Association of LOXL1 with Peritoneal Dissemination in Gastric Cancer Possibly via Promotion of EMT. PLoS One 15 (10), e0241140. doi:10.1371/journal.pone.0241140
Jiang, H., Torphy, R. J., Steiger, K., Hongo, H., Ritchie, A. J., Kriegsmann, M., et al. (2020). Pancreatic Ductal Adenocarcinoma Progression Is Restrained by Stromal Matrix. J. Clin. Invest. 130 (9), 4704–4709. doi:10.1172/JCI136760
Kaneda, A., Wakazono, K., Tsukamoto, T., Watanabe, N., Yagi, Y., Tatematsu, M., et al. (2004). Lysyl Oxidase Is a Tumor Suppressor Gene Inactivated by Methylation and Loss of Heterozygosity in Human Gastric Cancers. Cancer Res. 64 (18), 6410–6415. doi:10.1158/0008-5472.CAN-04-1543
Kang, M.-H., Choi, H., Oshima, M., Cheong, J.-H., Kim, S., Lee, J. H., et al. (2018). Estrogen-related Receptor Gamma Functions as a Tumor Suppressor in Gastric Cancer. Nat. Commun. 9 (1), 1920. doi:10.1038/s41467-018-04244-2
Kasashima, H., Yashiro, M., Kinoshita, H., Fukuoka, T., Morisaki, T., Masuda, G., et al. (2016). Lysyl Oxidase Is Associated with the Epithelial-Mesenchymal Transition of Gastric Cancer Cells in Hypoxia. Gastric Cancer 19 (2), 431–442. doi:10.1007/s10120-015-0510-3
Kasashima, H., Yashiro, M., Kinoshita, H., Fukuoka, T., Morisaki, T., Masuda, G., et al. (2014). Lysyl Oxidase-like 2 (LOXL2) from Stromal Fibroblasts Stimulates the Progression of Gastric Cancer. Cancer Lett. 354 (2), 438–446. doi:10.1016/j.canlet.201410.1016/j.canlet.2014.08.014
Kasashima, H., Yashiro, M., Okuno, T., Miki, Y., Kitayama, K., Masuda, G., et al. (2018). Significance of the Lysyl Oxidase Members Lysyl Oxidase like 1, 3, and 4 in Gastric Cancer. Digestion 98 (4), 238–248. doi:10.1159/000489558
Lai, H., Jin, Q., Lin, Y., Mo, X., Li, B., He, K., et al. (2014). Combined Use of Lysyl Oxidase, Carcino-Embryonic Antigen, and Carbohydrate Antigens Improves the Sensitivity of Biomarkers in Predicting Lymph Node Metastasis and Peritoneal Metastasis in Gastric Cancer. Tumor Biol. 35 (10), 10547–10554. doi:10.1007/s13277-014-2355-5
Li, Q., Zhu, C.-C., Ni, B., Zhang, Z.-Z., Jiang, S.-H., Hu, L.-P., et al. (2019). Lysyl Oxidase Promotes Liver Metastasis of Gastric Cancer via Facilitating the Reciprocal Interactions between Tumor Cells and Cancer Associated Fibroblasts. EBioMedicine 49, 157–171. doi:10.1016/j.ebiom.2019.10.037
Li, R.-k., Zhao, W.-y., Fang, F., Zhuang, C., Zhang, X.-x., Yang, X.-m., et al. (2015). Lysyl Oxidase-like 4 (LOXL4) Promotes Proliferation and Metastasis of Gastric Cancer via FAK/Src Pathway. J. Cancer Res. Clin. Oncol. 141 (2), 269–281. doi:10.1007/s00432-014-1823-z
Liu, J., Lichtenberg, T., Hoadley, K. A., Poisson, L. M., Lazar, A. J., Cherniack, A. D., et al. (2018). An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics. Cell 173 (2), 400–416.e11. doi:10.1016/j.cell.2018.02.052
Mahjour, F., Dambal, V., Shrestha, N., Singh, V., Noonan, V., Kantarci, A., et al. (2019). Mechanism for Oral Tumor Cell Lysyl Oxidase Like-2 in Cancer Development: Synergy with PDGF-AB. Oncogenesis 8 (5), 34. doi:10.1038/s41389-019-0144-0
Molnar, J., Fong, K. S. K., He, Q. P., Hayashi, K., Kim, Y., Fong, S. F. T., et al. (2003). Structural and Functional Diversity of Lysyl Oxidase and the LOX-like Proteins. Biochim. Biophys. Acta (Bba) - Proteins Proteomics 1647, 220–224. doi:10.1016/s1570-9639(03)00053-0
Moreira, A. M., Pereira, J., Melo, S., Fernandes, M. S., Carneiro, P., Seruca, R., et al. (2020). The Extracellular Matrix: An Accomplice in Gastric Cancer Development and Progression. Cells 9 (2), 394. doi:10.3390/cells9020394
Peng, C., Liu, J., Yang, G., and Li, Y. (2018). Lysyl Oxidase Activates Cancer Stromal Cells and Promotes Gastric Cancer Progression: Quantum Dot-Based Identification of Biomarkers in Cancer Stromal Cells. Int. J. Nanomedicine 13, 161–174. doi:10.2147/IJN.S143871
Peng, L., Ran, Y.-L., Hu, H., Yu, L., Liu, Q., Zhou, Z., et al. (2009). Secreted LOXL2 Is a Novel Therapeutic Target that Promotes Gastric Cancer Metastasis via the Src/FAK Pathway. Carcinogenesis 30 (10), 1660–1669. doi:10.1093/carcin/bgp178
Rachman-Tzemah, C., Zaffryar-Eilot, S., Grossman, M., Ribero, D., Timaner, M., Mäki, J. M., et al. (2017). Blocking Surgically Induced Lysyl Oxidase Activity Reduces the Risk of Lung Metastases. Cel. Rep. 19 (4), 774–784. doi:10.1016/j.celrep.2017.04.005
Rhodes, D. R., Kalyana-Sundaram, S., Mahavisno, V., Varambally, R., Yu, J., Briggs, B. B., et al. (2007). Oncomine 3.0: Genes, Pathways, and Networks in a Collection of 18,000 Cancer Gene Expression Profiles. Neoplasia 9 (2), 166–180. doi:10.1593/neo.07112
Salvador, F., Martin, A., López-Menéndez, C., Moreno-Bueno, G., Santos, V., Vázquez-Naharro, A., et al. (2017). Lysyl Oxidase-like Protein LOXL2 Promotes Lung Metastasis of Breast Cancer. Cancer Res. 77 (21), 5846–5859. doi:10.1158/0008-5472.CAN-16-3152
Setargew, Y. F. I., Wyllie, K., Grant, R. D., Chitty, J. L., and Cox, T. R. (2021). Targeting Lysyl Oxidase Family Meditated Matrix Cross-Linking as an Anti-stromal Therapy in Solid Tumours. Cancers 13 (3), 491. doi:10.3390/cancers13030491
Shao, J., Lu, J., Zhu, W., Yu, H., Jing, X., Wang, Y.-L., et al. (2019). Derepression of LOXL4 Inhibits Liver Cancer Growth by Reactivating Compromised P53. Cell Death Differ 26 (11), 2237–2252. doi:10.1038/s41418-019-0293-x
Smithen, D. A., Leung, L. M. H., Challinor, M., Lawrence, R., Tang, H., Niculescu-Duvaz, D., et al. (2020). 2-Aminomethylene-5-sulfonylthiazole Inhibitors of Lysyl Oxidase (LOX) and LOXL2 Show Significant Efficacy in Delaying Tumor Growth. J. Med. Chem. 63 (5), 2308–2324. doi:10.1021/acs.jmedchem.9b01112
Szász, A. M., Lánczky, A., Nagy, Á., Förster, S., Hark, K., Green, J. E., et al. (2016). Cross-validation of Survival Associated Biomarkers in Gastric Cancer Using Transcriptomic Data of 1,065 Patients. Oncotarget 7 (31), 49322–49333. doi:10.18632/oncotarget.10337
Tang, Z., Kang, B., Li, C., Chen, T., and Zhang, Z. (2019). GEPIA2: an Enhanced Web Server for Large-Scale Expression Profiling and Interactive Analysis. Nucleic Acids Res. 47, W556–W560. doi:10.1093/nar/gkz430
Wang, Q., Wen, Y.-G., Li, D.-P., Xia, J., Zhou, C.-Z., Yan, D.-W., et al. (2012). Upregulated INHBA Expression Is Associated with Poor Survival in Gastric Cancer. Med. Oncol. 29 (1), 77–83. doi:10.1007/s12032-010-9766-y
Wang, T.-H., Hsia, S.-M., and Shieh, T.-M. (2016). Lysyl Oxidase and the Tumor Microenvironment. Ijms 18 (1), 62. doi:10.3390/ijms18010062
Warde-Farley, D., Donaldson, S. L., Comes, O., Zuberi, K., Badrawi, R., Chao, P., et al. (2010). The GeneMANIA Prediction Server: Biological Network Integration for Gene Prioritization and Predicting Gene Function. Nucleic Acids Res. 38, W214–W220. doi:10.1093/nar/gkq537
Xiao, Q., and Ge, G. (2012). Lysyl Oxidase, Extracellular Matrix Remodeling and Cancer Metastasis. Cancer Microenvironment 5 (3), 261–273. doi:10.1007/s12307-012-0105-z
Yamauchi, M., Barker, T. H., Gibbons, D. L., and Kurie, J. M. (2018a). The Fibrotic Tumor Stroma. J. Clin. Invest. 128 (1), 16–25. doi:10.1172/JCI93554
Yamauchi, M., Taga, Y., Hattori, S., Shiiba, M., and Terajima, M. (2018b). Analysis of Collagen and Elastin Cross-Links. Methods Cel. Biol. 143, 115–132. doi:10.1016/bs.mcb.2017.08.006
Zeltz, C., Pasko, E., Cox, T. R., Navab, R., and Tsao, M.-S. (2019). LOXL1 Is Regulated by Integrin α11 and Promotes Non-small Cell Lung Cancer Tumorigenicity. Cancers 11 (5), 705. doi:10.3390/cancers11050705
Zhang, Q., Jin, X.-S., Yang, Z.-Y., Wei, M., Zhu, X.-C., Wang, P., et al. (2013). Upregulated Expression of LOX Is a Novel Independent Prognostic Marker of Worse Outcome in Gastric Cancer Patients after Curative Surgery. Oncol. Lett. 5 (3), 896–902. doi:10.3892/ol.2012.1092
Keywords: gastric cancer, LOX family, bioinformatics analysis, biomarker, prognosis
Citation: Wang L, Cao S, Zhai R, Zhao Y and Song G (2022) Systematic Analysis of Expression and Prognostic Values of Lysyl Oxidase Family in Gastric Cancer. Front. Genet. 12:760534. doi: 10.3389/fgene.2021.760534
Received: 18 August 2021; Accepted: 16 December 2021;
Published: 20 January 2022.
Edited by:
Frank Emmert-Streib, Tampere University, FinlandReviewed by:
Yuen Ting Lau, Queen Mary Hospital, Hong Kong SAR, ChinaChenyu Sun, AMITA Health St Joseph Hospital, United States
Shihori Tanabe, National Institute of Health Sciences (NIHS), Japan
Shan Yu, Fudan University, China
Shihai Liu, The Affiliated Hospital of Qingdao University, China
Copyright © 2022 Wang, Cao, Zhai, Zhao, Song. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guodong Song, tijdoctor@163.com
†These authors have contributed equally to this work and share first authorship