Chao Jiang
Meng Wang
Weikai Yao
Guoyue Lv
Xueyan Liu
Guangyi Wang1*- 1Department of Hepatobiliary Pancreatic Surgery I, The First Hospital of Jilin University, Changchun, China
- 2Multi-Organ Transplant Program, University Health Network, Toronto, ON, Canada
- 3Department of Pathology, The First Hospital of Jilin University, Changchun, China
- 4Cardiovascular Department, China-Japan Union Hospital of Jilin University, Changchun, China
Background: Pancreatic adenocarcinoma (PDAC) is the most aggressive among all solid malignancies with delayed disease detection and limited effective treatment. However, due to the intricate heterogeneity and exclusive tumor microenvironment (TME), the development of effective therapy has been facing enormous challenges. The lysyl oxidases (LOXs) underpin the shaping of the TME to promote cancer growth, metastasis and modulate response to treatment.
Materials and Methods: The mRNA expression, prognostic, and clinicopathological data for LOXs in PDAC from multiple open-access databases were summarized and analyzed. The protein expression was verified by immunohistochemistry (IHC). Co-expressed genes of LOXs were predicted and elaborated by LinkedOmics. Functional enrichment analysis of LOXs co-expressed genes was performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). TIMER and TISIDB were applied to analyze the relationship between LOXs expression and immune infiltration.
Results: The mRNA expression levels of LOX, LOXL1 and LOXL2 were significantly higher in PDAC, the expression levels of LOXL3 and LOXL4 were contrary in different databases. High mRNA levels of LOX and LOXL2 were associated with advanced PDAC stage, while elevated LOX and LOXL3 expression correlated with high tumor grade. The IHC staining showed higher expression levels of LOX, LOXL1 and LOXL2, lower expression level of LOXL3 in PDAC tissues, while the protein expression of LOXL4 made no difference. Functional enrichment analysis showed a close relationship with extracellular matrix (ECM) regulation, except that LOXL3 and its ligands were highly associated with immune-related functions. Further analysis suggested that LOX and LOXL3 strongly correlated with tumor-infiltrating lymphocytes (TILs), various immune signatures, and immune checkpoints. Finally, survival analysis revealed high LOX and LOXL2 expression predicted worse overall survival (OS), progression-free interval (PFI), and disease-specific survival (DSS).
Conclusion: These findings indicated that the LOX family, especially LOX and LOXL2, might have a prospective value in PDAC oncogenesis, and they may become prognostic biomarkers, revealing a promising field in targeted therapy.
Introduction
Pancreatic adenocarcinoma (PDAC), a highly fatal disease with a poor prognosis, ranks as the seventh leading cause of cancer mortality globally (Sung et al., 2021), is projected to become the second leading cancer-related death in America and Europe in 2030 (Rahib et al., 2014; Quante et al., 2016). However, over the past decades, albeit considerable advances have been made in the diagnosis, surgery, chemotherapy, radiotherapy and systemic therapy, however, due to non-specific symptoms, lacking reliable markers, low surgical resection rate and therapy resistance, only months of overall survival (OS) achieved (Conroy et al., 2011; Winter et al., 2012), with the 5-years survival remaining a grim 10% (Siegel et al., 2021). Therefore, further investigation of pathological features and the pathogenesis of PDAC is imperative for exceptional therapeutic strategies.
Such poor outcomes of PDAC have fueled ongoing efforts to exploit the tumor microenvironment (TME). The TME is a highly complex ecosystem composed of cancer cells, stromal cells, immune cells, and the extracellular matrix (ECM) interaction. Growing evidence has illustrated the pivotal roles of TME components in immune suppression and tumor progression. It is worth noting that the surrounding ECM has long been implicated in accelerating PDAC progression by directly promoting cellular transformation and metastasis. Furthermore, the ECM and its reorganization are of paramount importance in the evolution of both tumor and immune microenvironment and their interaction (Setargew et al., 2021). These alterations may greatly affect the biological behavior of tumor cells and stromal cells, and ultimately influence the prognosis of PDAC patients.
Lysyl oxidase (LOX) family includes five copper-dependent amino enzymes: LOX and LOX-like (LOXL) 1-4 (Molnar et al., 2003; Barker et al., 2012), of which the canonical function is catalyzing collagen and elastin cross-linking (Harris, 1976; Kagan and Li, 2003; Zhao et al., 2018). Lysyl Oxidases (LOXs) contribute to ECM stiffness and homeostasis, and their dysregulation is involved in several diseases, including tissue fibrosis (Al-U’datt et al., 2019; Li et al., 2020) and cancer (Wang et al., 2016; Tan et al., 2020). In reality, the activity of LOXs integrates a complex network within the TME linking the ECM and the immunological components. Intriguingly, LOXs play different roles in various tumor progression, exhibiting either enhancing (Wilgus et al., 2011; Kasashima et al., 2016)or suppressing (Johnston and Lopez, 2018; Zhao et al., 2020) function in different types of human malignancies.
It has been demonstrated that overexpressed LOXs are critically involved in the PDAC progression. Le Calve et al. (Le Calve et al., 2016) found higher mRNA expression of LOX, LOXL1 and LOXL2 in PDAC. Ma et al. (Ma et al., 2019) found that LOX expression might be an independent prognostic factor of progression-free survival (PFI) in PDAC patients. Higher LOXL2 expression elevated the invasiveness of PDAC cells and correlated with lower survival of PDAC patients (Kim et al., 2019). However, the potential roles of LOXs in the PDAC immune response have not been elucidated. Herein, various large-scale open-accessed databases were applied to explore the expression pattern, cancer-related functional aspects, prognostic value, and relationship with the immune infiltrating of LOXs in PDAC.
Materials and Methods
Oncomine
Oncomine gene expression array dataset (https://www.oncomine.org/resource/login.html) (Rhodes et al., 2004; Rhodes et al., 2007) was applied to analyze the transcripts levels of LOXs in different cancers. The data were compared by the t-test with a cut-off p-value of 0.001 and a fold change of 2.
GEPIA
Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia2.cancer-pku.cn/#index) (Tang et al., 2017) was used to analyze the relative transcriptional expression of LOXs between PDAC tissues and normal tissues. p < 0.01 was considered statistically significant.
Human Protein Atlas
The Human Protein Atlas (HPA, https://www.proteinatlas.org/) (Uhlen et al., 2015; Uhlen et al., 2017) was introduced to explore the direct contrast of protein expression of LOXL1-4 between normal tissues and PDAC tissues based on immunohistochemistry (IHC). The data of LOX were not available in HPA.
UALCAN
The UALCAN (http://ualcan.path.uab.edu/index.html) (Chandrashekar et al., 2017) was used to explore the protein expression of LOXs in PADC. The protein assembly data relative abundance was obtained from the Pancreatic Ductal Adenocarcinoma (PDAC) Discovery Study of Clinical Proteomics Tumor Analysis Consortium (CPTAC, https://proteomics.cancer.gov/programs/cptac). p < 0.01 was considered statistically significant.
Immunohistochemistry
As the feedback regulation loop of LOXs, the regulation of transcripts level and protein level are not necessarily the same. To further explore the LOXs protein level in PDAC patients, LOXs expression was evaluated through IHC in PDAC tissue and paired adjacent tissue(n = 6). The study was approved by the First Hospital of Jilin University Ethics Committee (2019-180). All patients provided written informed consent for the subsequent use of their resected tissues in this study. All clinical investigations were conducted following the principles of the Declaration of Helsinki. The PDAC tissues and paired adjacent tissues were prepared into 4 μm paraffin sections and incubated overnight at 4 °C with rabbit polyclonal antibodies against LOXs (1:150, Novus Biologicals, USA. LOX: NB100-2527; LOXL1: H00004016-D01P; LOXL2: NBP1-32954; LOXL3: NBP2-75964; LOXL4: NBP2-32692). After washing, the sections were conjugated with ImmPRESS anti-Rabbit kit for 30 min (LSBio, USA. LS-J1066). The negative controls demonstrated the specific signals of the primary antibody. Subsequently, all fields were observed under light microscopy. The average integrated option density (AOD) of LOXs was determined using ImageJ software (National Institutes of Health, USA).
Metascape
Metascape (https://metascape.org/gp/index.html#/main) (Zhou et al., 2019) was applied to explore the enrichment of the co-expressed genes with LOXs categorized by process and pathway. The most statistically significant 20 terms were shown in the diagram. Terms with p < 0.01, enrichment factor>1.5, and minimum count three were considered significant. The protein-protein interaction (PPI) was conducted using the Molecular Complex Detection (MCODE) algorithm.
LinkedOmics Database
The LinkedOmics database (http://www.linkedomics.org/login.php) (Vasaikar et al., 2018), a helpful third-party online tool containing The Cancer Genome Atlas (TCGA) data(Tomczak et al., 2015), was used to identify the co-expressed genes of LOXs. The correlation was expressed by the Pearson coefficient.
TIMER
Tumor Immune Estimation Resource (TIMER) database (http://timer.cistrome.org/) (Li et al., 2017; Li et al., 2020) was applied to assess the association of LOXs expression with tumor-infiltrating lymphocytes (TILs) via the Gene module and immune biomarker sets (Siemers et al., 2017) (including immune checkpoint gene sets) via Correlation module. The correlation was adjusted by tumor purity and expressed by the Spearman coefficient (Chen et al., 2020).
TISIDB
TISIDB (http://cis.hku.hk/TISIDB/) (Ru et al., 2019) was introduced to validate the correlation of LOXs expression with the abundance of TILs, and the relationship of LOXs expression with clinical characteristics in PDAC patients. The correlation was expressed by the Spearman coefficient.
Survival Analysis
The mRNA expression of LOXs in PDAC was obtained from c-BioPortal (www.cbioportal.org) (Cerami et al., 2012; Gao et al., 2013), and the clinical data of PDAC was extracted from integrated clinical data resources(Liu et al., 2018) to perform survival analysis. Patients were divided into two groups based on the median value of mRNA expression.
Data Visualization and Analysis
The Venn diagram was drawn via a friendly online tool (http://bioinformatics.psb.ugent.be/webtools/Venn/) to show the overlapping co-expressed genes of LOXs; the heatmap was conducted by TBtools (Chen et al., 2020); the KEGG Mapper (Shannon et al., 2003) was applied to conduct the cancer network using the overlapping genes retrieved from UniProt (https://www.uniprot.org/) (UniProt, 2021). The AOD of LOXs was expressed as the mean ± standard deviation. Data were analyzed with a two-tailed Student’s t-test for a parametric test. p < 0.05 was considered to indicate a statistically significant difference.
Results
High mRNA Expression of LOXs in PDAC Tissues
The mRNA expression of LOX, LOXL1, and LOXL2 was significantly higher in PDAC tissues than healthy tissues based on the Oncomine analysis (Figure 1 and Table 1). In Logsdon PDAC dataset(Logsdon et al., 2003), LOX and LOXL2 transcripts were upregulated in PDAC tissues by fold changes of 4.433 and 5.235 (p = 2.78E-4, 4.56E-6), respectively. Iacobuzio-Donahue (Iacobuzio-Donahue et al., 2003) found a 2.084-fold increase of LOXL1 and a 3.334-fold increase of LOXL2 in PDAC tissues. Moreover, other studies were added to support the upregulation of LOX, LOXL1, and LOXL2 in PDAC (Badea et al., 2008; Pei et al., 2009).
FIGURE 1. The transcriptional expression of distinct LOX family members in 20 types of cancers (Oncomine database). The numbers in each table represented the studies with statistically significant tumor tissue mRNA overexpression (red) or down-expression (blue), and LOX, LOXL1, and LOXL2 expression were upregulated in PDAC tissues.
TABLE 1. The significant changes of LOXs expression in transcriptional level between PAAD and normal pancreas tissues (Oncomine).
Furthermore, the mRNA expression patterns of LOXs were further measured by GEPIA based on the TCGA database. The GEPIA analyses showed all LOXs were higher in PDAC than in normal tissues (Figure 2).
FIGURE 2. Boxplot showing the expression levels of LOXs in PDAC. The mRNA expression of LOXs was significantly elevated in PDAC than normal pancreas tissues (GEPIA database). The number of normal samples is 171 (blue box), and PDAC samples are 179 (red box), red star means p < 0.01.
Protein Expression of LOXs in PDAC Tissues
To further explore the protein level in PDAC, IHC was performed in PDAC and adjacent pancreas. The IHC results demonstrated elevated LOX, LOXL1 and LOXL2 protein expression in the PDAC tissues (LOX: p < 0.01; LOXL1: p < 0.05; LOXL2: p < 0.05 Figures 3A–C). However, LOXL3 showed a remarkable decrease in PDAC tissues (p < 0.05, Figure 3D), and LOXL4 made no difference (Figure 3E).
FIGURE 3. Representative immunohistochemistry of LOXs in PDAC tissues and normal pancreas tissues. (A) IHC staining of LOX and quantification showed higher expression in PDAC than normal pancreas tissues. (B) IHC staining and quantification showed higher expression of LOXL1 in PDAC. (C) IHC staining and quantification showed higher expression of LOXL2 in PDAC. (D) IHC staining and quantification showed lower expression of LOXL3 in PDAC. (E) LOXL4 expression levels were similar in PDAC and normal tissues. *p < 0.01; NS, no significant difference. n = 6, repeated 5 times.
We also analyzed the protein expression patterns of LOXs in PDAC using HPA and CPTAC datasets. Consistent with our IHC results, HPA results indicated increased LOXL1 and decreased LOXL3 expression in PDAC tissues (Supplementary Figure S1 and Supplementary Table S1), while no significant difference was identified in LOXL2 and LOXL4 compared with exocrine glandular cells in normal tissue. However, the expression of LOLX3 was contrary in different databases. The protein expression level of LOXL3 was identified higher in PDAC tissues than normal tissues verified by CPTAC data (Figure 4).
FIGURE 4. Protein expression levels of LOXs in PDAC tissues and normal tissues (CPTAC database). Protein expression levels of LOX, LOXL1, 2, 3 were significantly up-regulated. ****p < 0.0001.
Relationship Between LOXs Expression and the Clinical Characteristics of PDAC Patients
We further investigated the association between the mRNA expression of LOXs and the clinical features of PDAC. As indicated in Figure 5, the higher PDAC stage was related to the higher mRNA levels of LOX and LOXL2. Meanwhile, tumor grade was significantly correlated with the expression of LOX and LOXL3.
FIGURE 5. Transcription levels of LOXs in different stages and grades of patients with PDAC (TISIDB). (A) The expression of LOX and LOXL2 was highly correlated to the stage of the PDAC. (B) The expression of LOX and LOXL3 was significantly associated with tumor grade in PDAC.
LOXs Correlated Genes and Functional Enrichment Analysis in PDAC
To gain insight into the potential role of LOXs in the PDAC process, we conducted interactome, gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the top 200 correlated genes of each LOX family member extracted from LinkedOmics. The results revealed that LOX, LOXL1, 2, 4 were mainly involved in ECM regulation, such as extracellular matrix structural constituent (GO: 0005201), collagen-containing extracellular matrix (GO: 0062023), extracellular matrix (GO: 0031012), external encapsulating structure (GO: 0030312), extracellular matrix organization (GO: 0030198), and extracellular structure organization (GO:0043062), while LOXL3 showed a close relationship with immune-related functions, such as myeloid leukocyte activation (GO: 0002274), leukocyte activation involved in immune response (GO: 0002366), myeloid leukocyte mediated immunity (GO: 0002444), myeloid cell activation involved in immune response (GO: 002275), regulation of leukocyte activation (GO: 0002694), leukocyte differentiation (GO: 0002521), and lymphocyte activation (GO: 0046649) (Figure 6 and Supplementary Table S2).
FIGURE 6. Top 20 statistically significant pathways associated with each LOXs identified by Metascape. Heatmap of the molecular function, biological processes, or pathways enriched with distinct LOX co-expressed genes. The bar color shade was decided by the p-value, the deeper the shade the less the p-value.
Besides the specific role of each LOX member plays in cancer, the involvement of several members of LOXs has also been illustrated. Thus, the 168 overlapping genes of each LOX member’s top 200 correlated genes were used to investigate the combined function of the LOX family in PDAC (Figure 7). As indicated in Figure 8A, we observed that functions of those genes were enriched in ECM remodeling, such as extracellular matrix organization (GO: 0030198), extracellular matrix structural constituent conferring tensile strength (GO: 0030020), blood vessel development (GO: 0001568), and collagen metabolic process (GO: 0032963). The networks of enrichment terms of LOXs according to cluster-ID were displayed in Figure 8B. Furthermore, the PPI network for LOXs and overlapping co-expressed genes was conducted, wherein the significant densely connected network constituents were assigned in a unique color (Figure 8C).
FIGURE 7. The Venn diagram represents the correlated genes of LOXs based on the Pearson correlation coefficient.
FIGURE 8. The enrichment analysis and visualized protein-protein interaction (PPI) enrichment analysis of co-expressed genes of LOXs (Metascape). (A) Heatmap of enriched terms regarding Gene Ontology analysis colored by p-values. (B) An interactive network of the top 20 enrichment terms. (C) PPI network and densely connected network components are identified by the Molecular Complex Detection (MCODE).
Similar functional enrichment results were achieved using DAVID (https://david.ncifcrf.gov/) (Huang da et al., 2009a; Huang da et al., 2009b) (Supplementary Figure S2). Considering the strong correlation between LOXs-related genes and pathways in cancer (hsa05200) as shown in the Metascape analysis, the KEGG pathway map was generated (Figure 9). The cancer network revealed that those correlated genes were essentially abundant in pathways in cancer (hsa05200), PI3K-AKT signaling pathway (has04151), focal adhesion (has04510), and ECM-receptor interaction (hsa04512), the details illustrated in Supplementary Figure S3.
FIGURE 9. Positional relationship of LOXs co-expressed genes in the cancer network. The overlapping genes were mainly enriched in cancer-related pathways, such as pathways in cancer (hsa05200), PI3K-AKT signaling pathway (has04151), focal adhesion (has04510), ECM-receptor interaction (hsa04512). The relevant positions of the corresponding genes in the pathway are marked with pink circles.
Association of LOXs With Immune Infiltration in PDAC
Considering the critical role of immune cells in the TME contributing to the progression of PDAC and a close relationship of LOXL3 with immune response demonstrated in the current study, we further investigated the correlation of LOXs transcriptional level with immune infiltration in PDAC using the TIMER and TISIDB database. The TIMER results showed that the LOXs expression was negatively related to tumor purity. The mRNA expression of LOX, LOXL1, 3, 4 was remarkably correlated with CD8+ T cells. The mRNA expression levels of LOX, LOXL1-3 had a significant correlation with infiltrating levels of macrophages and NK cells. In addition, the LOX family was highly associated with neutrophil and dendritic cell infiltration (Figures 10A–E). It was worth noting that LOX and LOXL3 remarkedly correlated with those TILs.
FIGURE 10. Correlation analysis of LOXs expression and immune infiltration in TME of PDAC. (A–E) (TIMER). The correlation between the abundance of immune cells and the expression of (A) LOX, (B) LOXL1, (C) LOXL2, (D) LOXL3, (E) LOXL4 in PDAC. (F) Correlation between the relative abundances of 28 tumor-infiltrating lymphocytes (TILs) and LOXs expression in PDAC (TISIDB).
Then, we explored the correlation between the abundances of TILs and LOXs expression via TISIDB. In line with the TIMER results, the TISIDB analysis demonstrated that high LOX and LOXL3 expression was strongly correlated with the Tcm CD8+ T cell, Tem CD8+ T cell, Tfh cell, Th1 cell, Treg cell, NK cell, MDSC cell, act DC cell, pDC cell, macrophage, and mast cell infiltration. Meanwhile, the expression of LOXL1, 2, 4 showed a relatively weak correlation with them (Figure 10F and Supplementary Table S3). These results suggested that LOXs, particularly LOX and LOXL3, might affect PDAC by influencing immune infiltration.
Correlation Between mRNA Expression of LOXs and Immune Biomarkers in PDAC
To better understand LOXs’ crosstalk with immune-related genes, we analyzed the association between the mRNA expression of LOXs and various immune biomarkers using TIMER. The results indicated that the expression of LOX and LOXL3 was highly related to biomarkers of CD8+ T cell, T cell, B cell, monocyte, TAM, M1 macrophage, M2 macrophage, neutrophil, DC cell, Th1, Th2, and Th17 in PDAC (Figure 11).
FIGURE 11. Correlation analysis of LOXs and immune signatures in PDAC (TIMER). The correlation was adjusted by purity.
We further explored the immune checkpoint expression level of PDAC cells. The results showed that LOX and LOXL3 were strongly correlated to PD-1, PD-L1, CLTA4, GZMB, TIM-3, and LAG3, while other members of LOX showed no relationship (Table 2). In brief, these results affirmed the role of LOX and LOXL3, especially LOXL3 in the immune response in PDAC.
TABLE 2. Correlation analysis of LOXs and immune biomarkers in PAAD (TIMER).
Prognostic Value of mRNA Expression of LOXs in PDAC Patients
As displayed in Figures 12A,B, C the mRNA levels of LOX and LOXL2 were significantly correlated with poor OS (p = 0.048 and 6e-04), progression-free interval (PFI) (p = 0.0029 and 0.018), and disease-specific survival (DSS) (p = 0.014 and 0.0031) in PDAC. However, the mRNA levels of LOXL1, 3, 4 showed no relationship with the overcomes in PDAC patients. The mRNA levels of LOX and LOXL2 highly correlated with the outcomes of PDAC patients, indicating LOX and LOXL2 may be served as biomarkers for predicting the prognosis.
FIGURE 12. The prognosis value of LOXs expression in PDAC. (A) Upregulated LOX and LOXL2 were highly correlated with poor overall survival (OS) in PDAC patients. (B) Upregulated LOX, LOXL2, and LOXL4 were significantly associated with inadequate progression-free interval (PFI). (C) LOX and LOXL2 expression was negatively related to disease-specific survival (DSS).
Discussion
As a kind of secreted protein family, LOXs play a pivotal role in PDAC pathogenesis by regulating the TME in two aspects, ECM remodeling and immune-related responses. It has been well known that the TME could significantly affect the behavior of tumor cells and ultimately influence the outcomes of patients (Setargew et al., 2021). The TME of PDAC is characterized by dense desmoplasia and extensive immunosuppression (Neesse et al., 2015; Ren et al., 2018). LOXs are major drivers in the biogenesis of the desmoplastic matrix in cancer. LOXs prevent collagen degradation and improve tissue stiffness through catalyzing the cross-linking process of collagen and elastin, and finally contribute to dense desmoplasia in PDAC. Dense desmoplasia plays a critical role in metastasis and limiting chemotherapy efficacy (Incio et al., 2016; Ren et al., 2018; Schizas et al., 2020). With the increasing interest in anti-stromal therapies, the LOX family has emerged as a potential target. Although the role of the LOX family in tumorigenesis and prognosis of several cancers has been partially confirmed, further bioinformatics analysis of PADC has yet to be performed. In this study, we used multitalented public databases to reveal the dysregulated expression of the LOXs and their relationship with tumor stage and prognosis. As for extensive immunosuppression in PDAC, we found that LOXs, particularly LOX and LOXL3 are closely related to TILs, various immune signatures, and immune checkpoints, however, the function in the immune microenvironment remains unclear which calls for further studies.
Elevated LOX expression has been noted in PDAC tissues compared to adjacent tissues and correlated with poor OS and DSS (Miller et al., 2015; Ma et al., 2019). A high expression of LOX was detected in metastatic pancreatic cancer in the mice model. Conversely, LOX inhibition increased immune cell recruitment, vascularization, and decreased fibrillar collagen (Miller et al., 2015). In line with this finding, Le Calve et al. (2016) found that LOX protein was upregulated in PDAC tissues. Compared to adjacent pancreas tissues, we found higher mRNA expression of LOX in PDAC tissues, which was correlated with tumor stage and grade of PDAC. The IHC showed a higher LOX protein level in PDAC tissues than normal pancreatic tissues. The functional enrichment analysis showed that LOX and its ligands were mainly involved in the ECM regulatory. The transcriptional level of LOX is strongly related to the infiltration of CD8+ T cell, macrophage, NK cell, neutrophil, dendritic cell infiltration, as well as various immune biomarkers and immune checkpoints. Furthermore, a high LOX expression was significantly associated with poor OS, PFI, and DSS in PDAC patients. These results together indicate the oncogenic role and prognostic value of LOX in PDAC, as well as the role in the immune response.
LOXL1 transcript was also detected upregulated in PDAC patients and associated with chemotherapy resistance (Le Calve et al., 2016). In our study, the mRNA and protein expression of LOXL1 in PDAC was higher than that in normal tissues. However, the mRNA expression of LOXL1 showed no relationship with tumor stage, grade, and outcomes of PDAC patients. We also revealed the role of LOXL1 in ECM shaping and a close relationship with infiltrating levels of CD8+ T cell, macrophage and NK cell, neutrophil, and dendritic cell. However, more studies should be carried out to identify the potential function of LOXL1 in PDAC progression.
Higher LOXL2 expression correlated with clinicopathological features of advanced disease and predicted a worse OS of PDAC patients (Park et al., 2016; Kim et al., 2019). Previous studies have revealed that LOXL2 induced epithelial-mesenchymal transition (EMT) through the PI3K-AKT signaling pathway (Wu and Zhu, 2015) to enhance the invasive and migratory capacity of PDAC cells (Park et al., 2016). In the current study, the expression of LOXL2 was higher in human pancreatic cancer tissues than in normal tissues, which highly correlated to the tumor stage. The functional enrichment analysis showed LOXL2 also had a close relationship with ECM regulation. In addition, LOXL2 expression showed a significant relationship with macrophage and NK cell, neutrophil, and dendritic cell infiltration. And Minici’s study (Minici et al., 2020) has demonstrated LOXL2 is expressed by infiltrating plasmablasts in PDAC. It is worth noting the high expressed LOXL2 predicted an inferior OS, PFI, and DSS in PDAC patients. Consequently, LOXL2 might be a representative biomarker for PDAC. However, phase II clinical trials failed to demonstrate the utility of LOXL2 inhibition in PDAC (Benson et al., 2017), which may be due to cancer-associated fibroblasts (CAFs). The discovery of the tumor-restraining functions of CAFs provides a potential explanation for the unsuccessful clinical trials of therapeutic agents targeting CAFs or stromal components (Van Cutsem et al., 2020). These observations suggest that future therapeutic strategies should avoid generic targeting of tumor-restraining CAF subpopulations in favor of precise reprogramming and normalization of tumor-promoting CAF subsets (Chen et al., 2021).
In this study, overexpressed mRNA of LOXL3 was identified in PDAC tissues and correlated with tumor grade, but showed no relationship with the outcomes of PDAC patients according to the GEPIA database and CPTAC database. Intriguingly, the LOXL3 protein level was found to be downregulated according to our IHC results and HPA database as well, which was controversial. As is known, the expression of LOXL3 was detected mainly in the nucleus correlated with tumor invasion, lymph node metastasis, and poorer prognosis of patients (Kasashima et al., 2018) through interacting with SNAIL and contributing to proliferation and metastasis by inducing epithelial-mesenchymal transition in pancreatic ductal adenocarcinoma cells (Eiseler et al., 2012). Additionally, TGF-induced LOXL3 upregulation in gastric cancer cells suggested that LOXL3 was downstream from the TGF-signaling pathway (Kasashima et al., 2018). However, studies on LOXL3 were mainly based on evidence from cell lines, limited evidence from in vivo studies was found. Meanwhile, LOXL3 expression in some tumors was contrary as in some cases LOXL3 was downregulated. Li Ma (Ma et al., 2017) found that LOXL3 was a dual-specificity enzyme involved in STAT3 deacetylation and deacetylimination modulation which reduced the activity of signal transducer and activator of transcription 3 (STAT3). STAT3 activation helps to promote tumor progression (Sun et al., 2021) and is associated with PDAC through PDAC patients’ biopsy study(Laklai et al., 2016). De-activation of STAT3 by LOXL3 would be a protective role in PDAC pathogenesis. Furthermore, unlike other LOX family members, the functional analysis of LOXL3 and its co-expressed genes showed enrichment in the immune response. LOXL3 expression showed a high correlation with CD8+ T cell, macrophage, NK cell, neutrophil, and dendritic cell infiltrating levels. In addition, LOXL3 was highly related to biomarkers of CD8+ T cell, T cell, B cell, monocyte, TAM, M1macrophage, M2 macrophage, neutrophil, DC cell, Th1, Th2, and Th17. In addition, the LOXL3 expression was positively associated with immune checkpoints. LOXL3 was also found to induce downregulation of Th17 and Treg activation as well (Ma et al., 2017) which play a role in cancer pathogenesis (Marques et al., 2021). Thus, we speculate LOXL3 plays a crucial role in the modulation and recruitment of immune cells and affects the expression of immune signatures in PDAC, and more studies should be carried out to validate.
Accumulating evidence suggest LOXL4 an oncoprotein in several tumors, such as hepatocellular carcinoma (Li et al., 2019), lung cancer (Zhang et al., 2019) and breast cancer (Choi et al., 2017). However, the expression and prognostic role of LOXL4 in PDAC have yet to be investigated. This report demonstrated a higher mRNA expression of LOXL4 in PDAC tissues with no correlation with tumor stage, grade, or outcomes in PDAC patients. The protein expression of LOLX4 showed no difference in PDAC and normal pancreas. The functional analysis of LOXL4 identified similar results to the canonical function of LOXs. In addition, the expression of LOXL4 showed a weak correlation with immune biomarkers. These results together suggest that LOXL4 function in pancreatic tumorigenesis needs to be further determined.
The involvement of one or more LOX members was verified in different cancers. Thus, the LOX family members are likely to present distinct functions specific for each member and partially overlapping activity. We further analyzed the synergy function of LOXs with overlapping co-expressed genes by GO enrichment analysis and KEGG pathway enrichment. GO analysis demonstrated the role of LOXs in ECM regulation. The cancer network showed that changes in LOXs expression mainly affect pathways in cancer (hsa05200), PI3K-AKT signaling pathway (has04151), focal adhesion (has04510), ECM-receptor interaction (hsa04512). These pathways are significantly associated with PDAC progression. Considering the major ECM component in PDAC, the predominant role of LOXs in the ECM remodeling, and involvement in the cancer-related pathways, we speculate that differentially expressed LOXs in PDAC are potential targets for drug therapy.
Current therapeutic strategies aim to deconstruct the surrounding desmoplastic stroma and target immunosuppressive pathways that have largely failed in PDAC. Even immunotherapeutic strategies, such as immune-checkpoint blockade, gained remarkable success in many malignancies that have not yet translated to PDAC (Johnson et al., 2017; Stromnes et al., 2017; Morrison et al., 2018). In addition, ablating the stromal barriers that restrict drug delivery has also demonstrated disappointing and contradictory responses. The key contributor is the multi-faceted functions of the complex TME components (Chronopoulos et al., 2016; Stromnes et al., 2017). Thus, combination strategies that target multiple features of the TME simultaneously might succeed. LOXs inhibitors reduce the fibrotic ECM and desmoplasia and overcome chemoresistance in triple-negative breast cancer by increasing the diffusion of chemotherapeutics into the tumor (Saatci et al., 2020). As demonstrated LOXs played a crucial role in modulating immune cells and affecting the expression of immune signatures in PDAC. Considering the close relationship of LOXs and the immune checkpoints, LOXs inhibitors might also play a role in promoting the anti-tumor effect of immune checkpoint inhibitors. A combination of LOXs inhibitor treatment may improve patient outcomes by priming the TME. However, further preclinical in vivo investigations and clinical trials need to be carried out to confirm our speculation.
In conclusion, our research identifies upregulated LOX, LOXL1 and LOXL2 mRNA and protein expression in PDAC tissue. LOX and LOXL2 could serve as novel biomarkers with prognostic significance in PDAC patients. The function of LOXs and their ligands is mainly involved in ECM regulation and pathways in cancer. Intriguing, LOXL3 and its co-expressed genes hold remarked relationship with immune response. Furthermore, a close association between the expression of LOX and LOXL3 with tumor immune infiltration was verified, uncovering the potential molecular mechanism underlying the carcinogenesis in PDAC. Our findings inspire new insight into prognostic biomarkers and new immunotherapeutic targets for PDAC. However, more experimental studies are needed to reveal the link between LOXs and PDAC, thereby promoting the clinical application of LOXs as prognostic indicators and combine therapy targets in PDAC.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committee of the First Hospital of Jilin University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
GW was responsible for the study concept and design; CJ and XL were involved in data collection, data screening, statistical analysis, and manuscript writing; WY was involved in the pathology. MW and GL modified and took charge of supervising the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding
This work was supported by the China Scholarship Council (201906175060) and the Jilin Provincial Department of Finance (JLSWSRCZX2021-016). However, the funders had no role in study design, data collection, analysis, decision to publish, or manuscript preparation.
Conflict of Interest
The authors declare that the study 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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.778857/full#supplementary-material
References
Al-U’datt, D. A., Allen, B. G., and Nattel, S. (2019). Role of the Lysyl Oxidase Enzyme Family in Cardiac Function and Disease. Cardiovasc. Res. 13, 1820–1837. doi:10.1093/cvr/cvz176
Badea, L., Herlea, V., Dima, S. O., Dumitrascu, T., and Popescu, I. (2008). Combined Gene Expression Analysis of Whole-Tissue and Microdissected Pancreatic Ductal Adenocarcinoma Identifies Genes Specifically Overexpressed in Tumor Epithelia. Hepatogastroenterology 55, 2016–2027.
Barker, H. E., Cox, T. R., and Erler, J. T. (2012). The Rationale for Targeting the LOX Family in Cancer. Nat. Rev. Cancer 12, 540–552. doi:10.1038/nrc3319
Benson, A. B., Wainberg, Z. A., Hecht, J. R., Vyushkov, D., Dong, H., Bendell, J., et al. (2017). A Phase II Randomized, Double-Blind, Placebo-Controlled Study of Simtuzumab or Placebo in Combination with Gemcitabine for the First-Line Treatment of Pancreatic Adenocarcinoma. The oncologist 22, 241–e15. doi:10.1634/theoncologist.2017-0024
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, 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, 649–658. doi:10.1016/j.neo.2017.05.002
Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020). TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 13, 1194–1202. doi:10.1016/j.molp.2020.06.009
Chen, Y., McAndrews, K. M., and Kalluri, R. (2021). Clinical and Therapeutic Relevance of Cancer-Associated Fibroblasts. Nat. Rev. Clin. Oncol. 18, 792–804. doi:10.1038/s41571-021-00546-5
Choi, S. K., Kim, H. S., Jin, T., and Moon, W. K. (2017). LOXL4 Knockdown Enhances Tumor Growth and Lung Metastasis through Collagen-dependent Extracellular Matrix Changes in Triple-Negative Breast Cancer. Oncotarget 8, 11977–11989. doi:10.18632/oncotarget.14450
Chronopoulos, A., Robinson, B., Sarper, M., Cortes, E., Auernheimer, V., Lachowski, D., et al. (2016). ATRA Mechanically Reprograms Pancreatic Stellate Cells to Suppress Matrix Remodelling and Inhibit Cancer Cell Invasion. Nat. Commun. 7, 12630. doi:10.1038/ncomms12630
Conroy, T., Desseigne, F., Ychou, M., Bouché, O., Guimbaud, R., Bécouarn, Y., et al. (2011). FOLFIRINOX versus Gemcitabine for Metastatic Pancreatic Cancer. N. Engl. J. Med. 364, 1817–1825. doi:10.1056/NEJMoa1011923
Eiseler, T., Köhler, C., Nimmagadda, S. C., Jamali, A., Funk, N., Joodi, G., et al. (2012). Protein Kinase D1 Mediates anchorage-dependent and -independent Growth of Tumor Cells via the Zinc finger Transcription Factor Snail1. J. Biol. Chem. 287, 32367–32380. doi:10.1074/jbc.M112.370999
Gao, J., Aksoy, B. A., Dogrusoz, U., Dresdner, G., Gross, B., Sumer, S. O., et al. (2013). Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal. Sci. Signal. 6, pl1. doi:10.1126/scisignal.2004088
Harris, E. D. (1976). Copper-induced Activation of Aortic Lysyl Oxidase In Vivo. Proc. Natl. Acad. Sci. 73, 371–374. doi:10.1073/pnas.73.2.371
Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009a). Bioinformatics Enrichment Tools: Paths toward the Comprehensive Functional Analysis of Large Gene Lists. Nucleic Acids Res. 37, 1–13. doi:10.1093/nar/gkn923
Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009b). Systematic and Integrative Analysis of Large Gene Lists Using DAVID Bioinformatics Resources. Nat. Protoc. 4, 44–57. doi:10.1038/nprot.2008.211
Iacobuzio-Donahue, C. A., Maitra, A., Olsen, M., Lowe, A. W., van Heek, N. T., Rosty, C., et al. (2003). Exploration of Global Gene Expression Patterns in Pancreatic Adenocarcinoma Using cDNA Microarrays. Am. J. Pathol. 162, 1151–1162. doi:10.1016/S0002-9440(10)63911-9
Incio, J., Liu, H., Suboj, P., Chin, S. M., Chen, I. X., Pinter, M., et al. (2016). Obesity-Induced Inflammation and Desmoplasia Promote Pancreatic Cancer Progression and Resistance to Chemotherapy. Cancer Discov. 6, 852–869. doi:10.1158/2159-8290.CD-15-1177
Johnson, B. A., Yarchoan, M., Lee, V., Laheru, D. A., and Jaffee, E. M. (2017). Strategies for Increasing Pancreatic Tumor Immunogenicity. Clin. Cancer Res. 23, 1656–1669. doi:10.1158/1078-0432.CCR-16-2318
Johnston, K. A., and Lopez, K. M. (2018). Lysyl Oxidase in Cancer Inhibition and Metastasis. Cancer Lett. 417, 174–181. doi:10.1016/j.canlet.2018.01.006
Kagan, H. M., and Li, W. (2003). Lysyl Oxidase: Properties, Specificity, and Biological Roles inside and outside of the Cell. J. Cel. Biochem. 88, 660–672. doi:10.1002/jcb.10413
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, 431–442. doi:10.1007/s10120-015-0510-3
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, 238–248. doi:10.1159/000489558
Kim, I.-k., Lee, Y. S., Kim, H. S., Dong, S. M., Park, J. S., and Yoon, D. S. (2019). Specific Protein 1(SP1) Regulates the Epithelial-Mesenchymal Transition via Lysyl Oxidase-like 2(LOXL2) in Pancreatic Ductal Adenocarcinoma. Sci. Rep. 9, 5933. doi:10.1038/s41598-019-42501-6
Laklai, H., Miroshnikova, Y. A., Pickup, M. W., Collisson, E. A., Kim, G. E., Barrett, A. S., et al. (2016). Genotype Tunes Pancreatic Ductal Adenocarcinoma Tissue Tension to Induce Matricellular Fibrosis and Tumor Progression. Nat. Med. 22, 497–505. doi:10.1038/nm.4082
Le Calvé, B., Griveau, A., Vindrieux, D., Maréchal, R., Wiel, C., Svrcek, M., et al. (2016). Lysyl Oxidase Family Activity Promotes Resistance of Pancreatic Ductal Adenocarcinoma to Chemotherapy by Limiting the Intratumoral Anticancer Drug Distribution. Oncotarget 7, 32100–32112. doi:10.18632/oncotarget.8527
Li, R., Wang, Y., Zhang, X., Feng, M., Ma, J., Li, J., et al. (2019). Exosome-mediated Secretion of LOXL4 Promotes Hepatocellular Carcinoma Cell Invasion and Metastasis. Mol. Cancer 18, 18. doi:10.1186/s12943-019-0948-8
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77, e108–e110. doi:10.1158/0008-5472.CAN-17-0307
Li, T., Fu, J., Zeng, Z., Cohen, D., Li, J., Chen, Q., et al. (2020). TIMER2.0 for Analysis of Tumor-Infiltrating Immune Cells. Nucleic Acids Res. 48, W509–W514. doi:10.1093/nar/gkaa407
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, 400–e11. doi:10.1016/j.cell.2018.02.052
Logsdon, C. D., Simeone, D. M., Binkley, C., Arumugam, T., Greenson, J. K., Giordano, T. J., et al. (2003). Molecular Profiling of Pancreatic Adenocarcinoma and Chronic Pancreatitis Identifies Multiple Genes Differentially Regulated in Pancreatic Cancer. Cancer Res. 63, 2649–2657.
Ma, L., Huang, C., Wang, X.-J., Xin, D. E., Wang, L.-s., Zou, Q. C., et al. (2017). Lysyl Oxidase 3 Is a Dual-Specificity Enzyme Involved in STAT3 Deacetylation and Deacetylimination Modulation. Mol. Cel. 65, 296–309. doi:10.1016/j.molcel.2016.12.002
Ma, W., Li, T., Wu, S., Li, J., Wang, X., and Li, H. (2019). LOX and ACSL5 as Potential Relapse Markers for Pancreatic Cancer Patients. Cancer Biol. Ther. 20, 787–798. doi:10.1080/15384047.2018.1564565
Marques, H. S., de Brito, B. B., da Silva, F. A. F., Santos, M. L. C., de Souza, J. C. B., Correia, T. M. L., et al. (2021). Relationship between Th17 Immune Response and Cancer. Wjco 12, 845–867. doi:10.5306/wjco.v12.i10.845
Miller, B. W., Morton, J. P., Pinese, M., Saturno, G., Jamieson, N. B., McGhee, E., et al. (2015). Targeting the LOX/Hypoxia axis Reverses many of the Features that Make Pancreatic Cancer Deadly: Inhibition of LOX Abrogates Metastasis and Enhances Drug Efficacy. EMBO Mol. Med. 7, 1063–1076. doi:10.15252/emmm.201404827
Minici, C., Rigamonti, E., Lanzillotta, M., Monno, A., Rovati, L., Maehara, T., et al. (2020). B Lymphocytes Contribute to Stromal Reaction in Pancreatic Ductal Adenocarcinoma. Oncoimmunology 9, 1794359. doi:10.1080/2162402X.2020.1794359
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
Morrison, A. H., Byrne, K. T., and Vonderheide, R. H. (2018). Immunotherapy and Prevention of Pancreatic Cancer. Trends Cancer 4, 418–428. doi:10.1016/j.trecan.2018.04.001
Neesse, A., Algül, H., Tuveson, D. A., and Gress, T. M. (2015). Stromal Biology and Therapy in Pancreatic Cancer: a Changing Paradigm. Gut 64, 1476–1484. doi:10.1136/gutjnl-2015-309304
Park, J. S., Lee, J.-h., Lee, Y. S., Kim, J. K., Dong, S. M., and Yoon, D. S. (2016). Emerging Role of LOXL2 in the Promotion of Pancreas Cancer Metastasis. Oncotarget 7, 42539–42552. doi:10.18632/oncotarget.9918
Pei, H., Li, L., Fridley, B. L., Jenkins, G. D., Kalari, K. R., Lingle, W., et al. (2009). FKBP51 Affects Cancer Cell Response to Chemotherapy by Negatively Regulating Akt. Cancer cell 16, 259–266. doi:10.1016/j.ccr.2009.07.016
Quante, A. S., Ming, C., Rottmann, M., Engel, J., Boeck, S., Heinemann, V., et al. (2016). Projections of Cancer Incidence and Cancer‐related Deaths in Germany by 2020 and 2030. Cancer Med. 5, 2649–2656. doi:10.1002/cam4.767
Rahib, L., Smith, B. D., Aizenberg, R., Rosenzweig, A. B., Fleshman, J. M., and Matrisian, L. M. (2014). Projecting Cancer Incidence and Deaths to 2030: the Unexpected burden of Thyroid, Liver, and Pancreas Cancers in the United States. Cancer Res. 74, 2913–2921. doi:10.1158/0008-5472.CAN-14-0155
Ren, B., Cui, M., Yang, G., Wang, H., Feng, M., You, L., et al. (2018). Tumor Microenvironment Participates in Metastasis of Pancreatic Cancer. Mol. Cancer 17, 108. doi:10.1186/s12943-018-0858-1
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, 166–180. doi:10.1593/neo.07112
Rhodes, D. R., Yu, J., Shanker, K., Deshpande, N., Varambally, R., Ghosh, D., et al. (2004). ONCOMINE: a Cancer Microarray Database and Integrated Data-Mining Platform. Neoplasia 6, 1–6. doi:10.1016/s1476-5586(04)80047-2
Ru, B., Wong, C. N., Tong, Y., Zhong, J. Y., Zhong, S. S. W., Wu, W. C., et al. (2019). TISIDB: an Integrated Repository portal for Tumor-Immune System Interactions. Bioinformatics 35, 4200–4202. doi:10.1093/bioinformatics/btz210
Saatci, O., Kaymak, A., Raza, U., Ersan, P. G., Akbulut, O., Banister, C. E., et al. (2020). Targeting Lysyl Oxidase (LOX) Overcomes Chemotherapy Resistance in Triple Negative Breast Cancer. Nat. Commun. 11, 2416. doi:10.1038/s41467-020-16199-4
Schizas, D., Charalampakis, N., Kole, C., Economopoulou, P., Koustas, E., Gkotsis, E., et al. (2020). Immunotherapy for Pancreatic Cancer: A 2020 Update. Cancer Treat. Rev. 86, 102016. doi:10.1016/j.ctrv.2020.102016
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, 491. doi:10.3390/cancers13030491
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13, 2498–2504. doi:10.1101/gr.1239303
Siegel, R. L., Miller, K. D., Fuchs, H. E., and Jemal, A. (2021). Cancer Statistics, 2021. CA A. Cancer J. Clin. 71, 7–33. doi:10.3322/caac.21654
Siemers, N. O., Holloway, J. L., Chang, H., Chasalow, S. D., Ross-MacDonald, P. B., Voliva, C. F., et al. (2017). Genome-wide Association Analysis Identifies Genetic Correlates of Immune Infiltrates in Solid Tumors. PloS one 12, e0179726. doi:10.1371/journal.pone.0179726
Stromnes, I. M., Hulbert, A., Pierce, R. H., Greenberg, P. D., and Hingorani, S. R. (2017). T-cell Localization, Activation, and Clonal Expansion in Human Pancreatic Ductal Adenocarcinoma. Cancer Immunol. Res. 5, 978–991. doi:10.1158/2326-6066.CIR-16-0322
Sun, L., Zhang, X., Song, Q., Liu, L., Forbes, E., Tian, W., et al. (2021). IGFBP2 Promotes Tumor Progression by Inducing Alternative Polarization of Macrophages in Pancreatic Ductal Adenocarcinoma through the STAT3 Pathway. Cancer Lett. 500, 132–146. doi:10.1016/j.canlet.2020.12.008
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A. Cancer J. Clin. 71, 209–249. doi:10.3322/caac.21660
Tan, H. Y., Wang, N., Zhang, C., Chan, Y. T., Yuen, M. F., and Feng, Y. (2021). Lysyl Oxidase‐Like 4 Fosters an Immunosuppressive Microenvironment during Hepatocarcinogenesis. Hepatology 73, 2326–2341. doi:10.1002/hep.31600
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a Web Server for Cancer and normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res. 45, W98–W102. doi:10.1093/nar/gkx247
Tomczak, K., Czerwińska, P., and Wiznerowicz, M. (2015). Review the Cancer Genome Atlas (TCGA): an Immeasurable Source of Knowledge. wo 1A, 68–77. doi:10.5114/wo.2014.47136
Uhlén, M., Fagerberg, L., Hallström, B. M., Lindskog, C., Oksvold, P., Mardinoglu, A., et al. (2015). Tissue-based Map of the Human Proteome. Science 347, 1260419. doi:10.1126/science.1260419
Uhlen, M., Zhang, C., Lee, S., Sjöstedt, E., Fagerberg, L., Bidkhori, G., et al. (2017). A Pathology Atlas of the Human Cancer Transcriptome. Science 357, eaan2507. doi:10.1126/science.aan2507
UniProt, C. (2021). UniProt: the Universal Protein Knowledgebase in 2021. Nucleic Acids Res. 49, D480–D9. doi:10.1093/nar/gkaa1100
Van Cutsem, E., Tempero, M. A., Sigal, D., Oh, D.-Y., Fazio, N., Macarulla, T., et al. (2020). Randomized Phase III Trial of Pegvorhyaluronidase Alfa with Nab-Paclitaxel Plus Gemcitabine for Patients with Hyaluronan-High Metastatic Pancreatic Adenocarcinoma. J. Clin. Oncol. 38, 3185–3194. doi:10.1200/JCO.20.00590
Vasaikar, S. V., Straub, P., Wang, J., and Zhang, B. (2018). LinkedOmics: Analyzing Multi-Omics Data within and across 32 Cancer Types. Nucleic Acids Res. 46, D956–D963. doi:10.1093/nar/gkx1090
Wang, T.-H., Hsia, S.-M., and Shieh, T.-M. (2016). Lysyl Oxidase and the Tumor Microenvironment. Int. J. Mol. Sci. 18, 62. doi:10.3390/ijms18010062
Wilgus, M.-L., Borczuk, A. C., Stoopler, M., Ginsburg, M., Gorenstein, L., Sonett, J. R., et al. (2011). Lysyl Oxidase: a Lung Adenocarcinoma Biomarker of Invasion and Survival. Cancer 117, 2186–2191. doi:10.1002/cncr.25768
Winter, J. M., Brennan, M. F., Tang, L. H., D’Angelica, M. I., Dematteo, R. P., Fong, Y., et al. (2012). Survival after Resection of Pancreatic Adenocarcinoma: Results from a Single Institution over Three Decades. Ann. Surg. Oncol. 19, 169–175. doi:10.1245/s10434-011-1900-3
Wu, L., and Zhu, Y. (2015). The Function and Mechanisms of Action of LOXL2 in Cancer (Review). Int. J. Mol. Med. 36, 1200–1204. doi:10.3892/ijmm.2015.2337
Zhang, Y., Jiang, W. L., Yang, J. Y., Huang, J., Kang, G., Hu, H. B., et al. (2019). Downregulation of Lysyl Oxidase‐like 4 LOXL4 by miR‐135a‐5p Promotes Lung Cancer Progression In Vitro and In Vivo. J. Cel Physiol 234, 18679–18687. doi:10.1002/jcp.28508
Zhao, R., Liu, W., Wang, M., Zhang, Y., Pan, L., Feng, F., et al. (2020). Lysyl Oxidase Inhibits TNF-α Induced Rat Nucleus Pulposus Cell Apoptosis via Regulating Fas/FasL Pathway and the P53 Pathways. Life Sci. 260, 118483. doi:10.1016/j.lfs.2020.118483
Zhao, W., Yang, A., Chen, W., Wang, P., Liu, T., Cong, M., et al. (2018). Inhibition of Lysyl Oxidase-like 1 (LOXL1) Expression Arrests Liver Fibrosis Progression in Cirrhosis by Reducing Elastin Crosslinking. Biochim. Biophys. Acta (Bba) - Mol. Basis Dis. 1864, 1129–1137. doi:10.1016/j.bbadis.2018.01.019
Keywords: pancreatic adenocarcinoma, lysyl oxidase, LOX family members, bioinformatics analysis, tumor microenvironment, immune infiltration, prognostic
Citation: Jiang C, Wang M, Yao W, Lv G, Liu X and Wang G (2022) Comprehensive Analysis on Prognosis and Immune Infiltration of Lysyl Oxidase Family Members in Pancreatic Adenocarcinoma With Experimental Verification. Front. Mol. Biosci. 9:778857. doi: 10.3389/fmolb.2022.778857
Received: 17 September 2021; Accepted: 07 March 2022;
Published: 01 April 2022.
Edited by:
Pedro Borralho, Novartis Oncology, PortugalReviewed by:
Emanuel Della Torre, Ospedale San Raffaele (IRCCS), ItalyMarta Bento Afonso, Research Institute for Medicines (iMed.ULisboa), Portugal
Copyright © 2022 Jiang, Wang, Yao, Lv, Liu and Wang. 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: Xueyan Liu, liuxueyan@jlu.edu.cn; Guangyi Wang, guangyi@jlu.edu.cn