Skip to main content

LMOD1, an oncogene associated with Lauren classification, regulates the metastasis of gastric cancer cells through the FAK-AKT/mTOR pathway

Abstract

Background

The Lauren classification of gastric tumors strongly correlates with prognosis. The purpose of this study was to explore the specific molecular mechanism of Lauren classification of gastric cancer and provide a possible theoretical basis for the treatment of gastric cancer.

Methods

We standardized the gene expression data of five Gene Expression Omnibus gastric cancer databases and constructed a Weighted Co-expression Network Analysis (WGCNA) model based on clinicopathological information. The overall survival (OS) and disease-free survival (DFS) curves were extracted from the Cancer Genome Atlas (TCGA) and GSE62254 databases. Western blotting was used to measure protein expression in cells and tissues. Scratch and transwell experiments were used to test the migration ability of tumor cells. Immunohistochemistry was used to measure tissue protein expression in clinical tissue samples to correlate to survival data.

Results

The WGCNA model demonstrated that blue cyan was highly correlated with the Lauren classification of the tumor (r = 0.24, P = 7 × 1016). A protein-protein interaction network was used to visualize the genes in the blue cyan module. The OS and PFS TCGA analysis revealed that LMOD1 was a gene of interest. The proportion of diffuse gastric cancer patients with high expression of LMOD1 was significantly higher than that of intestinal type patients. LMOD1 promoted the migration of gastric cancer cells by regulating the FAK-Akt/mTOR pathway in vitro. Additionally, a Gene Set Enrichment Analysis using the TCGA and GSE62254 databases, and western blot data, showed that LMOD1 could promote an epithelial-mesenchymal transition (EMT), thus potentially affecting the occurrence of peritoneal metastasis of gastric cancer. Immunohistochemistry showed that LMOD1 was highly expressed in cancer tissues, and the prognosis of patients with high LMOD1 expression was poor.

Conclusion

LMOD1 is an oncogene associated with diffuse gastric cancer and can affect the occurrence and development of EMT by regulating the FAK-Akt/mTOR pathway. LMOD1 can therefore promote peritoneal metastasis of gastric cancer cells and can be used as a novel therapeutic target for gastric cancer.

Peer Review reports

Introduction

Gastric cancer is the fifth most common type of cancer, and fourth leading cause of cancer-related death [1]. The occurrence and development of gastric cancer is a multi-gene and multi-stage process, including the proliferation and metastasis of tumor cells [2, 3]. In recent years, there has been great progress in the treatment of gastric cancer. The application of targeted therapy and immunotherapy has greatly improved the survival rate of advanced gastric cancer, but not all patients respond to treatment, especially for patients with peritoneal metastasis [4, 5]. Peritoneal metastasis is one of the main factors affecting the prognosis and quality of life of gastric cancer patients. Therefore, it is of great significance to identify the key genes affecting the progression of peritoneal metastasis of gastric cancer and target them in clinical treatment to improve the prognosis of patients.

The classification systems used with gastric cancer include Lauren classification, Borrmann classification, and the serosal classification. The Lauren classification divides tumors into an intestinal type, diffuse type and mixed type [6]. The survival time of patients with diffuse gastric cancer is significantly shorter than that of patients with intestinal gastric cancer, where the 5-year survival rate for the intestinal type is 57.7%, and the 5-year survival rate for the diffuse type is only 45.6% [7]. Peritoneal metastasis is reported in 80% of diffuse gastric cancer patients, which is much higher than for intestinal gastric cancer patients [8]. Gene therapy for Lauren typing may greatly reduce peritoneal metastasis of gastric cancer and could be a potential approach for the treatment of gastric cancer.

Weighted Gene Co-expression Network Analysis (WGCNA) is an efficient and accurate data mining method [9]. By inducing genes with similar expression changes, we can identify co-expression modules and explore the relationship between gene network and phenotype. An increasing number of scholars use WGCNA in gene mining of various tumors, which greatly improves the efficiency of basic research [10, 11]. If the expression of genes in different patient tissues is similar, these genes can be thought of as functionally related. Genomes with similar expression profiles form a module. Leiomodin 1 (LMOD1) was initially identified as 64 kDa autoantigen in a serum expression screen in Hashimoto’s thyroiditis. The gene transcript was found to be expressed in all tissues tested, with the highest levels being in the thyroid, eye muscle, and skeletal muscle [12]. In recent years, immunohistochemistry and western blotting for LMOD1 have demonstrated its highly restricted expression pattern in smooth muscle cell lines, such as adult aorta and bladder, suggesting that LMOD1 may be related to the phenotype of smooth muscle differentiation.

In this study, we constructed a WGCNA analysis model by combing the information on gastric cancer patients from Gene Expression Omnibus (GEO) databases, and selected LMOD1 as our research gene target for the modules highly related to the Lauren classification data. Western blotting results and functional tests showed that LMOD1 can induce gastric cancer cell migration by regulating the FAK-Akt/mTOR pathway and can promote the occurrence and development of gastric cancer peritoneal metastasis by affecting the epithelial-mesenchymal transition (EMT) of gastric cancer cells.

Materials and methods

Cell culture

The GES-1, SGC-7901, MGC-803, HGC-27 and AGS cell lines were purchased from the Cell Culture Collection of the Chinese Academy of Sciences (Shanghai, China). The Sv5 cell line (human peritoneal mesothelial cell, HPMC) was kindly provided by You-Ming Peng (Second Hospital of Zhongnan University, China). DMEM and RPMI-1640 medium were used for cell culture (Biological Industries, Israel), and the proportion of serum in cell culture medium was 10% (Biological Industries, Israel).

Data sources and data preprocessing

We searched “stomach neoplasms” [MeSH terms] or gastric cancer [All Fields]) and Lauren [All Fields] in the public database, and identified five datasets (GSE15460, GSE29272, GSE62254, GSE13861 and GSE26253) containing a total of 1001 gastric cancer patients with complete clinical information. The data from GSE29272 and GSE62254 were normalized by robust multichip average (RMA). The Illumina data (GSE13861, GSE26253) were normalized by quantity normalization and log2 transformation.

Construction of the co-expression network

We extracted and screened gene expression data from 8135 genes from five GEO databases, analyzed the expression of each gene in 1001 patient samples by ANOVA, and took the top 25% of genes (2035) for subsequent WGCNA analysis. We used the WGCNA package for the WGCNA analysis, and the cluster analysis was performed on the expression data, and clinical information of 1001 patients to measure their distribution, and eliminate abnormal values. After the modules were obtained, the modules with similarity greater than 0.75 were merged to obtain 12 co-expression modules.

Gene set enrichment analysis

For the Gene Set Enrichment Analysis (GSEA), gene expression in the TCGA and GSE62254 datasets were divided into two groups according to the high and low expression of LOMD1. Gene pathway enrichment was analyzed by GSEA software (4.1.0) [13].

Vector construction and transfection

The silenced GFP-LMOD1 Lentivirus was purchased from Beijing Syngenesis (China). The type of GFP-LMOD1 virus was PHS-ASR-0921(pLV-hU6-shRNA-hef1a-mNeogreen-P2A-Puro) and lentivirus target was GGAGATGTCCATGGATGAAAG. Puromycin (5 μg/ml) was used to screen viral-transduced stable cell lines. The LMOD1 overexpression plasmid was purchased from Genechem (Shanghai, China). The nucleotide sequence for the LMOD1 overexpression plasmid is shown in the Supplement word. Higene was used as the transfection reagent for plasmid transfection, and G418 was used for stable transformation screening.

Western blot

Western blotting was performed as previously described [14, 15]. Briefly, cells were collected, lysed, and centrifuged. SDS-PAGE gel was used to separate equal amounts of proteins. Proteins were then transferred added to PVDF membranes (Millipore). The membranes were blocked in 5% non-fat milk and incubated with primary antibodies at 4 °C overnight. After 1 h incubation with secondary antibodies, the membranes were developed with a chemiluminescence solution. The primary antibodies used were: LMOD1 (Accept), FAK (Proteintech), AKT (CST), p-AKT (CST), mTOR (Proteintech), p-mTOR (Proteintech), E-cadherin (CST), N-cadherin (CST), vimentin (CST), and actin (Proteintech). We cut the membrane to facilitate the better hybridization of antibodies and avoid the interaction of antibodies. The original pictures were attached to the Supplementary figure_ Original Western blots.

Real-time quantitative PCR analysis

After total RNA was extracted from the cells, 1 μg of RNA was added into collection tubes, and an appropriate amount of reverse transcription reagent (Takara, Japan) was added to obtain the cDNA product for quantitative analysis. SYBR dye (Takara, Japan) was added to the reactions, and the number of cycles on the thermal cycler was set to 45. GAPDH was used as a reference gene. The primer sequences were: LMOD1 5′- GAAGAACTCCCGTGACCAGCTA-3′ (sense) and 5′-AGCCTGGTCCTACTGAAGCAGT-3′ (antisense), GAPDH 5′- GTCTCCTCTGACTTCAACAGCG-3′ (sense) and 5′- ACCACCCTGTTGCTGTAGCCAA-3′ (antisense).

Wound-healing assay

Cells were trypsinized, resuspended, and 5 × 105 cells were evenly seeded in a 6-well plate. After the cells reached confluency a 200 μL pipette tip was used to make a scratch in the cell layer. Wells were washed twice with PBS and 2 mL of serum-free medium was added. An image was acquired at time zero. After 24 h, wells were washed twice with PBS and imaged.

Migration assay

Cells were trypsinized and resuspended in serum-free medium. A total of 3 × 105 / 200 μL medium were added into a Transwell chamber (Corning, USA), and 600 μL serum-containing medium was added into the lower chamber. After 18-22 h, the chambers were moved and washed twice with PBS. Chambers were immersed in 95% alcohol for 20 min, washed in PBS, and 0.4% trypan blue solution was added for 15 min. After washing with PBS, the inner surface of the chamber membrane was wiped with a cotton swab. Images were acquired.

Adhesion and invasion assay

Adhesion and invasion experiments relating to peritoneal metastasis were carried out according to previous experiments [15]. For the adhesion experiment, HPMCs cocultured with gastric cancer cells were plated in 12 well plates. 2 × 105 GC cells incubated in 5 μmol /L Calcein-AM (Sigma) for 1 h, then added to 12 well plates and cultured for 2 h. The plates were washed with PBS to remove suspended GC cells. The number of GC cells was counted under a fluorescence microscope. For the invasion experiment, HPMCs cocultured with GC cells were plated in a chamber (8-μm pore size; Corning) containing matrigel for 4-6 h. 5 × 105 GC cells resuspended in 100 μl serum-free medium were added to the HPMC layer. 600 μl of medium containing 10% FBS was added to the lower chamber. After 28-32 h, we removed the upper chamber for fixation and staining, and counted the number of cells under the microscope.

Samples and patients

We collected 102 cases of gastric cancer tissues and 40 cases of paracancerous tissues in our department from 2011 to 2012. These tissues were embedded in paraffin for storage. Clinicopathological data were also collected. In addition, we extracted proteins from 48 pairs of fresh gastric cancer tissues and stored them at − 80 °C for western blotting. All patients signed informed consents, and all procedures were in accordance with the Standards of Ethics Committee of China Medical University.

Immunohistochemistry

Immunohistochemistry was carried out according to previously published methods [14, 15]. The primary antibody against LMOD1 was purchased from Accepta (China). The percent positive cells and staining intensity were used to determine the degree of expression. The percentage of positive cells was divided into five groups: 0: < 5%; 1:5 – 25%; 2:25 – 50%; 3:50 – 75%; 4: > 75%. The staining intensity was divided into four groups: 0: colorless; 1: light yellow; 2: brown; 3: dark brown. Scores 0-4 were considered low expression, and scores 6-12 were considered high expression.

In vivo gastric cancer peritoneal dissemination assay

Gastric cancer peritoneal dissemination experiments were carried out according to previous experimental methods [15]. 10 Blab / C nude mice were randomly divided into two groups with five mice in each group. A mixed suspension of 2× 106 SGC-7901 shCtrl/ shLMOD1 cells and 5 × 105 HPMCs and SGC-7901 and HPMCs were injected intraperitoneally. Mice were euthanized 6 weeks after tumor inoculation. The location and number of abdominal macroscopic lesions were recorded, and the tumors were weighed. All methods were in accordance with the requirements of the Animal Committee of China Medical University.

Statistical analysis

SPSS 20.0 and GraphPad Prism 7.0 were used for data statistics and analysis. A Student’s t-test was used to compare the differences between the two means, and the correlation with clinicopathological factors was analyzed using a chi-square test. The Kaplan – Meier method was used to calculate the difference in overall survival, and X-tile was used to identify Optimal cut-off value. P < 0.05 was considered statistically significant.

Results

Construction and analysis of WGCNA from GEO databases

We extracted gene expression information from 8135 genes from five gastric cancer GEO datasets, analyzed the expression of each gene in 1001 patient samples using ANOVA, and took the top 25% of genes (2035) for subsequent WGCNA analysis. The expression and clinical information from 1001 patients were analyzed using cluster analysis (Fig. 1A), followed by calculation of the soft threshold. When the soft threshold was set at three, the scale-free topology fitting index (R2) was 0.853 (Fig. 1B and C). Finally, 12 co-expression modules were obtained (Fig. 1D and E). The correlation coefficient between the blue module and the Lauren classification was 0.24 (Fig. 2A). The genes in the blue module were selected for analysis, and the module membership and the significance of genes were calculated. The correlation between the two was 0.81, and the P value was statistically significant (Fig. 2B and C). The relationship between genes in the blue module was visualized using Cytoscape, and the relationship between genes with weighted values in the top 500 were selected for input into Cytoscape. The Maximal Clique Centrality (MCC) algorithm in the plug-in was used to calculate. The higher the weight of the network in the module, the greater the representative effect, and six key genes (MFAP4, MYL9, MYH11, LMOD1, FXYD6 and VIP) were selected from the network (Fig. 2D).

Fig. 1
figure 1

Data processing and construction of WGCNA. A Clustering dendrogram of the clinical traits. B, C Calculation of the soft threshold by scale independence and mean connectivity. When the soft threshold was 3, the scale free topology fitting index R2 was 0.853. D Clustering dendrogram of genes based on topological overlap

Fig. 2
figure 2

Analysis of the WGCNA module. A Heatmap of modules for gene-traits relationship. Red represents a positive correlation, and green represents a negative correlation. B The correlation between blue module membership and Lauren classification was significant. The correlation coefficient was 0.81. C Network heatmap plot. D Utilization of Cytoscape to visualize the relationship between genes in the blue module. Six key genes were selected from the network using the Maximal Clique Centrality algorithm

LMOD1 is identified as an oncogene and is associated with Lauren classification

To identify the genes with a high correlation between expression and Lauren classification, we used the TCGA database to measure the OS and PFS of hub genes in the blue module. Patients with high MFAP4 and LMOD1 expression had a worse prognosis, while the other four genes were not statistically correlated with prognosis (Fig. 3A-L). It has previously been reported that MFAP4 is a cancer-promoting gene with a positive effect on the development of gastric cancer. Since, LMOD1 had not been previously studied in relation to gastric cancer, this gene was chosen for further study. By analyzing the relationship between LMOD1 expression and Lauren classification in GEO databases, we found that patients with high LMOD1 expression had diffuse gastric cancer, while patients with low LMOD1 expression were more likely to have intestinal gastric cancer (Fig. 3M). The expression of LMOD1 was negatively correlated with HER-2 in GEO and TCGA databases (Fig. 3N). We chose the genes in the blue module that were highly related to LMOD1 to construct a PPI network (Fig. 3O).

Fig. 3
figure 3

Identification of the key gene related to Lauren classification. A-L Overall survival and disease-free survival curves for six key genes in the TCGA and GSE62254 databases (MFAP4, MYL9, MYH11, LMOD1, FXYD6, VIP). M Relationship between gene expression of LMOD1 and Lauren classification of patients. N The correlation between LMOD1 and HER-2 (ERBB2) expression in TCGA and GEO databases. M Analysis with the MCOD algorithm yielding related sub-modules of LMOD1 in blue cyan

LMOD1 regulates the migration of gastric cancer cells through the FAK-Akt/mTOR pathway in vitro

The expression of LMOD1 protein and mRNA was detected in cell lines. The expression of LMOD1 in GES-1 cells was lower than that in gastric cancer cell lines, further supporting that LMOD1 is a cancer-promoting gene (Fig. 4A and B). Moreover, LMOD1 was highly expressed in SGC-7901 and MGC-803 cells, and lowly expressed in HGC-27 and AGS cells. To explore the role of LMOD1 in the occurrence and development of gastric cancer, we used TCGA and GSE62254 databases to analyze LMOD1 with GSEA. LMOD1 was enriched in Ca + signaling pathways and the mTOR signaling pathway (Fig. 4C-F) in both databases. We next generated LMOD1-silenced (shLMOD1) SGC-7901 and MGC-803 cell lines and LMOD1-overexpressing (LMOD1-OE) HGC-27 cell lines. The expression of FAK, p-Akt and p-mTOR in shLMOD1 cells was decreased (Fig. 4G and H), while the expression of FAK, p-Akt and p-mTOR in LMOD1 OE cells increased (Fig. 4I), indicating that LMOD1 can further affect the phosphorylation of Akt and mTOR by altering the expression of FAK. Meanwhile, compared with the control group, the migration ability of shLMOD1 cells was significantly decreased (Fig. 4J, K and M), while the migration ability of LMOD1 OE cells was significantly increased (Fig. 4L and N).

Fig. 4
figure 4

LMOD1 enhanced the migration ability of gastric cancer cells through the FAK-Akt/mTOR pathway. A The expression of LMOD1 in cell lines was measured using western blots. B The expression of LMOD1 in cell lines was measured by RT-qPCR. C-F A GSEA analysis of LMOD1 enriched pathways in TCGA and GSE62254 databases. G, H Western blotting was used to measure protein markers of the FAK-Akt/mTOR pathway in Control or LMOD1- silenced SGC-7901 and MGC-803 cell lines. I Western blotting was used to measure protein markers of the FAK-Akt/mTOR pathway in the Control or LMOD1-overexpressed HGC-27 cell line. J, K The migration ability of cells was detected using a scratch test. M, N The migration ability of cells was measured using a transwell test. The membrane was cut out due to the different secondary antibodies incubated. *, P < 0.05, **, P < 0.01, ***, P < 0.001, ****, P < 0.0001

LMOD1 promotes peritoneal metastasis of gastric cancer cells by regulating EMT

Patients with diffuse gastric cancer often suffer from peritoneal metastasis, which greatly affects the prognosis of gastric cancer patients. Studies have shown that mTOR phosphorylation can affect downstream EMT and further promote the invasion and metastasis of gastric cancer cells [16, 17]. To answer whether LMOD1 could promote peritoneal metastasis of gastric cancer by affecting EMT, we predicted that LMOD1 could be enriched in the EMT pathway in the TCGA and GSE62254 databases (Fig. 5A and B). We then measured the expression of EMT-related protein markers in LMOD1-silenced and LMOD1-OE cells. The results showed that compared with the control group, the levels of E-cadherin in LMOD1-silenced cells were increased, and the levels of N-cadherin and vimentin were decreased (Fig. 5C-E). In contrast, the levels of E-cadherin in LMOD1-OE cells were decreased, and the levels of N-cadherin and vimentin were increased. Furthermore, in adhesion and invasion tests, compared with the control group, the adhesion and penetration of LMOD1-silenced cells were decreased, while the adhesion and penetration of LMOD1-OE cells were increased (Fig. 5F-K). In the model of peritoneal metastasis in nude mice, the number and weight of tumor nodules from SGC-7901 shLMOD1 cells were significantly lower than that of the control group (Fig. 5L, M; Supplement figure 1).

Fig. 5
figure 5

LMOD1 promoted peritoneal metastasis of gastric cancer cells by regulating the epithelial-mesenchymal transition. A, B GSEA analysis of LMOD1 enriched pathways in TCGA and GSE62254 databases. C, D Western blotting was used to verify the regulation of epithelial-mesenchymal transition (EMT)-related proteins in Control or shLMOD1 SGC-7901 and MGC-803 cell lines. E Western blotting was used to verify the regulation of EMT-related proteins in Control/\ or overexpressed-LMOD1 HGC-27 cell lines. F-H The adhesion ability of cells was verified using an adhesion test. I-K An invasion assay was used to verify the ability of peritoneal infiltration. L, M The effect of LMOD1 silencing on invasion and migration of gastric cancer cells in SGC-7901 cells was verified by intraperitoneal metastasis model in nude mice. We cut the membrane to facilitate the better hybridization of antibodies and avoid the interaction of antibodies. *, P < 0.05, **, P < 0.01, ***, P < 0.001, ****, P < 0.0001

Expression and survival analyses of LMOD1 in clinical tissue samples

To explore the expression and prognostic significance of LMOD1 in clinical tissues of patients with gastric cancer, we measured the protein expression of LMOD1 in 48 pairs of gastric cancer and adjacent tissue samples. The expression of LMOD1 in cancer tissues was significantly higher than that in adjacent tissues, and the difference was statistically significant (P = 0.014; Fig. 6A and B). Next, we examined the expression of LMOD1 in 102 gastric cancer tissues and 40 paracancerous tissues using immunohistochemistry (Fig. 6C), which also showed that the expression of LMOD1 in gastric cancer tissues was higher than that in paracancerous tissues (Fig. 6D). According to the LMOD1 immunohistochemical scores, 102 patients with gastric cancer were divided into a low-expression group (38/102) or a high-expression group (64/102). Survival analysis showed that patients with high expression of LMOD1 had worse prognosis than patients with low LMOD1 expression (Fig. 6E), Meanwhile, in diffuse gastric cancer, the survival time of patients with high expression of LOMD1 was significantly shorter than that of patients with low expression (Fig. 6F). The clinicopathological analysis showed that the expression of LMOD1 was positively correlated with Borrmann type, Lauren type and T stage (Table 1). The high expression of LMOD1 accounted for 59.4% in patients with diffused gastric cancer and 28.9% in patients with intestinal type.

Fig. 6
figure 6

Expression and survival significance of LMOD1 in clinical tissue samples. A, B Western blotting was used to measure the protein expression of LMOD1 in 48 pairs of gastric cancer tissues. C, D LMOD1 protein expression was measured by immunohistochemistry in 102 cases of gastric cancer tissues and 40 cases of adjacent tissues. E Overall survival curve of LMOD1in gastric cancer. F Overall survival curve of LMOD1 in diffuse gastric cancer. G Schematic diagram showed that the regulatory effect of LMOD1 on peritoneal metastasis in diffused gastric cancer. The membrane was cut out due to the different secondary antibodies incubated. ****, P < 0.0001

Table 1 Correlation between LMOD1 expression and clinicopathologic characteristics in GC patients

Discussion

The Lauren classification designates special types of gastric cancer. Compared with the diffuse type, the intestinal type is less prone to peritoneal metastasis and has a better prognosis. Therefore, it is especially important to identify the key genes that determine Lauren classification typing for targeted therapy of gastric cancer. In this paper, we constructed a WGCNA model by analyzing the gene expression and clinical information of patients in GEO gastric cancer databases, and analyzed the blue cyan modules related to Lauren classification. Eventually, LMOD1 was selected as the target of our analysis. In vitro and in vivo experiments verified that LMOD1 could affect the migration ability of gastric cancer cells by regulating the FAK-Akt/mTOR pathway, which can alter the EMT of cells to promote the occurrence of peritoneal metastasis in gastric cancer (Fig. 6G).

Peritoneal metastasis is an independent prognostic factor in gastric cancer [18, 19]. About 20% of patients are diagnosed with peritoneal metastasis preoperatively or intraoperatively, and the median survival time is only 6-7.5 months [20]. Diffuse gastric cancer is most likely to develop peritoneal metastasis, accounting for about 37%, while intestinal gastric cancer most likely to develop distant metastasis [21]. The adhesion and invasion ability of gastric cancer cells can be increased through EMT, which permeates the barrier of mesothelial cells, allowing for entry into the abdominal cavity to produce peritoneal metastasis. However, due to the limitations of technology and surgical methods, researchers can rarely obtain tissue samples from patients with peritoneal metastasis of gastric cancer. In this study, we extracted and analyzed gene modules highly related to Lauren classification by using clinicopathological factors and gene expression in the GEO database and WGCNA algorithm, to further predict and explore the mechanism of peritoneal metastasis of gastric cancer.

Constructing co-expression modules using the WGCNA algorithm has been widely used to identify potential key or prognostic genes. In a study by Wan et al. data from melanoma patients from the TCGA was analyzed using WGCNA, yielding 21 modules, and four key genes related to melanoma recurrence [22]. Zhai et al. used data from the E-GEOD-39582 colorectal cancer microarray to obtain 434 differentially genes. After performing a WGCNA analysis, they obtained 5 differential genes related to recurrence and prognosis [23]. In the current study, we used the WGCNA algorithm to find gene enrichment modules associated with diffuse and intestinal gastric cancer and selected the hub gene, LMOD1, to determine the molecular mechanism of Lauren classification. The research and development of the WGCNA method was originally applied to the development of space technology. However, WGCNA can now play an important role in exploring the specific mechanism of tumor occurrence and development through interdisciplinary research.

The LMOD1 protein has a putative transmembrane domain and two kinds of repeated blockades. It was initially identified as a 64 kDa autoantigen in a serum screening study of Hashimoto’s thyroiditis patients. It is reported that LMOD1 is a new target gene of smooth muscle cell-restricted serum response factor/myocardial protein [12]. In the cancer setting, Zhao et al. demonstrated that LMOD1 could be used as a key gene to determine the diagnosis and prognosis of early-onset colorectal cancer [24]. Kawahara et al. showed that the expression of LMOD1 was significantly different between low-grade and high-grade of prostate cancer, which could be used as a biological marker for the diagnosis of prostate cancer [25]. In this study, a GSEA analysis and western blotting demonstrated that LMOD1 could alter the migration ability of gastric cancer cells by regulating FAK-Akt / mTOR pathway. In vivo and in vitro experiments further confirmed that LMOD1 could promote peritoneal metastasis by regulating EMT of gastric cancer cells.

The Lauren classification is one of the most important phenotypes of gastric cancer. HER-2 as a therapeutic target for breast cancer has been widely confirmed [26, 27]. In gastric cancer patients, the HER-2 positive rate is about 10 - 30%. ToGA study showed that for HER-2-positive gastric cancer patients, trastuzumab combined with chemotherapy resulted in the longest overall survival and enhanced efficacy for intestinal gastric cancer, with the worst prognosis being for diffuse gastric cancer patients. Meanwhile, the positive rate of HER-2 in intestinal patients was significantly higher than that of diffuse type patients [4]. Qiu et al. reported that the prognosis of patients with HER-2-negative intestinal gastric cancer was better than that of patients with HER-2-positive diffuse gastric cancer [28]. Our results showed that LMOD1 was highly expressed in diffuse gastric cancer and negatively correlated with HER-2 expression. However, the relationship between HER-2 and prognosis in gastric cancer is not clear. Based on these studies, it is therefore critical to explore the biological functional genes that determine development of the intestinal type and diffuse type of gastric cancer, to provide a promising target for the future treatment of gastric cancer.

“Seed and soil theory” is a highly recognized theoretical mechanism of peritoneal metastasis of gastric cancer [29, 30]. The abscission of cancer cells and the change of extracellular matrix components promotes the implantation of cancer cells in the peritoneum. EMT is the main mechanism for the peritoneal metastasis of cancer cells. It is reported that the exosomes in malignant ascites can promote the EMT of gastric cancer cells and increase the diffusion degree of gastric cancer cells in the abdominal cavity [31]. Chen et al. showed that FNDC1 can be closely related to the peritoneal metastasis of gastric cancer by regulating the EMT pathway of cells [32]. Moreover, studies have shown that FAK activation is an important upstream protein regulating EMT [33, 34]. Our results suggest that LMOD1 can upregulate the expression of FAK, activate the phosphorylation level of the Akt / mTOR pathway, promote EMT in gastric cancer cells, increase cell invasion and migration, and promote the occurrence and development of peritoneal metastasis of gastric cancer.

To summarize, we used clinical patient information and gene expression data from GEO databases of gastric cancer to build a WGCNA model. We then selected the modules relating to Lauren classification, and identified the oncogene LMOD1 in the module for further study. LMOD1 was found to regulate the FAK-Akt/mTOR pathway to induce EMT of gastric cancer cells, to further increase invasion and adhesion of gastric cancer cells and promote peritoneal metastasis of gastric cancer cells. These results indicate that LMOD1 could sever as a novel biomarker and therapeutic target for gastric cancer.

Availability of data and materials

The data that support the findings of this study are openly available in GEO databases and The Cancer Genome Atlas database (https://www.ncbi.nlm.nih.gov/gds, https://www.cancer.gov/aboutnci/organization/ccg/research/structural-genomics/tcga).

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  Google Scholar 

  2. Ganesh K, Massagué J. Targeting metastatic cancer. Nat Med. 2021;27(1):34–44.

    Article  CAS  Google Scholar 

  3. Klein CA. Cancer progression and the invisible phase of metastatic colonization. Nat Rev Cancer. 2020;20(11):681–94.

    Article  CAS  Google Scholar 

  4. Bang YJ, Van Cutsem E, Feyereislova A, Chung HC, Shen L, Sawaki A, et al. Trastuzumab in combination with chemotherapy versus chemotherapy alone for treatment of HER2-positive advanced gastric or gastro-oesophageal junction cancer (ToGA): a phase 3, open-label, randomised controlled trial. Lancet. 2010;376(9742):687–97.

    Article  CAS  Google Scholar 

  5. Fukuoka S, Hara H, Takahashi N, Kojima T, Kawazoe A, Asayama M, et al. Regorafenib plus Nivolumab in patients with advanced gastric or colorectal Cancer: An open-label, dose-escalation, and dose-expansion phase Ib trial (REGONIVO, EPOC1603). J Clin Oncol. 2020;38(18):2053–61.

    Article  CAS  Google Scholar 

  6. Smyth EC, Nilsson M, Grabsch HI, van Grieken NC, Lordick F. Gastric cancer. Lancet. 2020;396(10251):635–48.

    Article  CAS  Google Scholar 

  7. Chen YC, Fang WL, Wang RF, Liu CA, Yang MH, Lo SS, et al. Clinicopathological variation of Lauren classification in gastric Cancer. Pathol Oncol Res. 2016;22(1):197–202.

    Article  CAS  Google Scholar 

  8. Zhao J, Qin R, Chen H, Yang Y, Qin W, Han J, et al. A nomogram based on glycomic biomarkers in serum and clinicopathological characteristics for evaluating the risk of peritoneal metastasis in gastric cancer. Clin Proteomics. 2020;17:34.

    Article  CAS  Google Scholar 

  9. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.

    Article  Google Scholar 

  10. Niemira M, Collin F, Szalkowska A, Bielska A, Chwialkowska K, Reszec J, et al. Molecular signature of subtypes of non-small-cell lung Cancer by large-scale transcriptional profiling: identification of key modules and genes by weighted gene co-expression network analysis (WGCNA). Cancers (Basel). 2019;12(1):37.

    Article  Google Scholar 

  11. Ding M, Li F, Wang B, Chi G, Liu H. A comprehensive analysis of WGCNA and serum metabolomics manifests the lung cancer-associated disordered glucose metabolism. J Cell Biochem. 2019;120(6):10855–63.

    Article  CAS  Google Scholar 

  12. Nanda V, Miano JM. Leiomodin 1, a new serum response factor-dependent target gene expressed preferentially in differentiated smooth muscle cells. J Biol Chem. 2012;287(4):2459–67.

    Article  CAS  Google Scholar 

  13. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.

    Article  CAS  Google Scholar 

  14. Tan YE, Xing Y, Ran BL, Zhang C, Pan SW, An W, et al. LINC01235-TWIST2 feedback loop facilitates epithelial-mesenchymal transition in gastric cancer by inhibiting THBS2. Aging (Albany NY). 2020;12(24):25060–75.

    Article  CAS  Google Scholar 

  15. Yin S, Miao Z, Tan Y, Wang P, Xu X, Zhang C, et al. SPHK1-induced autophagy in peritoneal mesothelial cell enhances gastric cancer peritoneal dissemination. Cancer medicine. 2019;8(4):1731–43.

    Article  CAS  Google Scholar 

  16. Karimi Roshan M, Soltani A, Soleimani A, Rezaie Kahkhaie K, Afshari AR, Soukhtanloo M. Role of AKT and mTOR signaling pathways in the induction of epithelial-mesenchymal transition (EMT) process. Biochimie. 2019;165:229–34.

    Article  CAS  Google Scholar 

  17. Wu YJ, Lin SH, Din ZH, Su JH, Liu CI. Sinulariolide inhibits gastric Cancer cell migration and invasion through downregulation of the EMT process and suppression of FAK/PI3K/AKT/mTOR and MAPKs signaling pathways. Mar Drugs. 2019;17(12):668.

    Article  CAS  Google Scholar 

  18. Chen D, Liu Z, Liu W, Fu M, Jiang W, Xu S, et al. Predicting postoperative peritoneal metastasis in gastric cancer with serosal invasion using a collagen nomogram. Nat Commun. 2021;12(1):179.

    Article  CAS  Google Scholar 

  19. Ishigami H, Yamaguchi H, Yamashita H, Asakage M, Kitayama J. Surgery after intraperitoneal and systemic chemotherapy for gastric cancer with peritoneal metastasis or positive peritoneal cytology findings. Gastric Cancer. 2017;20(Suppl 1):128–34.

    Article  Google Scholar 

  20. Nashimoto A, Akazawa K, Isobe Y, Miyashiro I, Katai H, Kodera Y, et al. Gastric cancer treated in 2002 in Japan: 2009 annual report of the JGCA nationwide registry. Gastric Cancer. 2013;16(1):1–27.

    Article  Google Scholar 

  21. Lee JH, Chang KK, Yoon C, Tang LH, Strong VE, Yoon SS. Lauren histologic type is the Most important factor associated with pattern of recurrence following resection of gastric adenocarcinoma. Ann Surg. 2018;267(1):105–13.

    Article  Google Scholar 

  22. Wan Q, Tang J, Han Y, Wang D. Co-expression modules construction by WGCNA and identify potential prognostic markers of uveal melanoma. Exp Eye Res. 2018;166:13–20.

    Article  CAS  Google Scholar 

  23. Zhai X, Xue Q, Liu Q, Guo Y, Chen Z. Colon cancer recurrence-associated genes revealed by WGCNA co-expression network analysis. Mol Med Rep. 2017;16(5):6499–505.

    Article  CAS  Google Scholar 

  24. Zhao B, Baloch Z, Ma Y, Wan Z, Huo Y, Li F, et al. Identification of potential key genes and pathways in early-onset colorectal Cancer through bioinformatics analysis. Cancer Control. 2019;26(1):1073274819831260.

    Article  Google Scholar 

  25. Kawahara R, Recuero S, Nogueira FCS, Domont GB, Leite KRM, Srougi M, et al. Tissue proteome signatures associated with five grades of prostate Cancer and benign prostatic hyperplasia. Proteomics. 2019;19(21-22):e1900174.

    Article  Google Scholar 

  26. von Minckwitz G, Huang CS, Mano MS, Loibl S, Mamounas EP, Untch M, et al. Trastuzumab Emtansine for residual invasive HER2-positive breast Cancer. N Engl J Med. 2019;380(7):617–28.

    Article  Google Scholar 

  27. von Minckwitz G, Procter M, de Azambuja E, Zardavas D, Benyunes M, Viale G, et al. Adjuvant Pertuzumab and Trastuzumab in early HER2-positive breast Cancer. N Engl J Med. 2017;377(2):122–31.

    Article  Google Scholar 

  28. Qiu M, Zhou Y, Zhang X, Wang Z, Wang F, Shao J, et al. Lauren classification combined with HER2 status is a better prognostic factor in Chinese gastric cancer patients. BMC Cancer. 2014;14:823.

    Article  Google Scholar 

  29. Zhao Y, Li J, Li D, Wang Z, Zhao J, Wu X, et al. Tumor biology and multidisciplinary strategies of oligometastasis in gastrointestinal cancers. Semin Cancer Biol. 2020;60:334–43.

    Article  CAS  Google Scholar 

  30. Greig RG, Trainer DL. Shaping future strategies for the pharmacological control of tumor cell metastases. Cancer Metastasis Rev. 1986;5(1):3–14.

    Article  CAS  Google Scholar 

  31. Hu Y, Qi C, Liu X, Zhang C, Gao J, Wu Y, et al. Malignant ascites-derived exosomes promote peritoneal tumor cell dissemination and reveal a distinct miRNA signature in advanced gastric cancer. Cancer Lett. 2019;457:142–50.

    Article  CAS  Google Scholar 

  32. Jiang T, Gao W, Lin S, Chen H, Du B, Liu Q, et al. FNDC1 promotes the invasiveness of gastric Cancer Wnt/β-catenin signaling pathway and correlates with peritoneal metastasis and prognosis. Front Oncol. 2020;10:590492.

    Article  Google Scholar 

  33. Taliaferro-Smith L, Oberlick E, Liu T, McGlothen T, Alcaide T, Tobin R, et al. FAK activation is required for IGF1R-mediated regulation of EMT, migration, and invasion in mesenchymal triple negative breast cancer cells. Oncotarget. 2015;6(7):4757–72.

    Article  Google Scholar 

  34. Zhang PF, Li KS, Shen YH, Gao PT, Dong ZR, Cai JB, et al. Galectin-1 induces hepatocellular carcinoma EMT and sorafenib resistance by activating FAK/PI3K/AKT signaling. Cell Death Dis. 2016;7:e2201.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by National Natural Science Foundation of China (No 81772549, 81572334).

Author information

Authors and Affiliations

Authors

Contributions

Tan and Xing work at paper writing, and Chen and Pan work at data analysis. An is responsible for data download. Zhang and Xu provide ideas of research. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Huimian Xu, Yao Xing or Jianjun Zhang.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Medical Ethics Research Association of the first affiliated Hospital of China Medical University, and each GC patient signed a written informed consent form. All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors have declared that no competing interest exists.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tan, Y., Chen, Q., Pan, S. et al. LMOD1, an oncogene associated with Lauren classification, regulates the metastasis of gastric cancer cells through the FAK-AKT/mTOR pathway. BMC Cancer 22, 474 (2022). https://doi.org/10.1186/s12885-022-09541-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-09541-0

Keywords