Next Article in Journal
WWOX Controls Cell Survival, Immune Response and Disease Progression by pY33 to pS14 Transition to Alternate Signaling Partners
Next Article in Special Issue
SALL4: An Intriguing Therapeutic Target in Cancer Treatment
Previous Article in Journal
Cortical Organoids to Model Microcephaly
Previous Article in Special Issue
Prognostic and Immunotherapeutic Roles of KRAS in Pan-Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gremlin-1 Promotes Colorectal Cancer Cell Metastasis by Activating ATF6 and Inhibiting ATF4 Pathways

1
Tomas Lindahl Nobel Laureate Laboratory, The Seventh Affiliated Hospital of Sun Yat-sen University (SYSU), No. 628, Zhenyuan Road, Guangming Dist., Shenzhen 518107, China
2
Center for Digestive Disease, The Seventh Affiliated Hospital of Sun Yat-sen University (SYSU), No. 628, Zhenyuan Road, Guangming Dist., Shenzhen 518107, China
3
Guangdong Provincial Key Laboratory of Digestive Cancer Research, The Seventh Affiliated Hospital of Sun Yat-sen University, No. 628 Zhenyuan Road, Shenzhen 518107, China
4
Department of Otolaryngology, The Seventh Affiliated Hospital of Sun Yat-sen University, No. 628, Zhenyuan Road, Guangming Dist., Shenzhen 518107, China
5
The Francis Crick Institute, 1 Midland Road, London NW1 1AT, UK
6
Ministry of Education Key Lab of Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430074, China
7
China-UK Institute for Frontier Science, Shenzhen 518107, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Cells 2022, 11(14), 2136; https://doi.org/10.3390/cells11142136
Submission received: 5 May 2022 / Revised: 17 June 2022 / Accepted: 29 June 2022 / Published: 7 July 2022
(This article belongs to the Collection Emerging Cancer Target Genes)

Abstract

:
Cancer cell survival, function and fate strongly depend on endoplasmic reticulum (ER) proteostasis. Although previous studies have implicated the ER stress signaling network in all stages of cancer development, its role in cancer metastasis remains to be elucidated. In this study, we investigated the role of Gremlin-1 (GREM1), a secreted protein, in the invasion and metastasis of colorectal cancer (CRC) cells in vitro and in vivo. Firstly, public datasets showed a positive correlation between high expression of GREM1 and a poor prognosis for CRC. Secondly, GREM1 enhanced motility and invasion of CRC cells by epithelial–mesenchymal transition (EMT). Thirdly, GREM1 upregulated expression of activating transcription factor 6 (ATF6) and downregulated that of ATF4, and modulation of the two key players of the unfolded protein response (UPR) was possibly through activation of PI3K/AKT/mTOR and antagonization of BMP2 signaling pathways, respectively. Taken together, our results demonstrate that GREM1 is an invasion-promoting factor via regulation of ATF6 and ATF4 expression in CRC cells, suggesting GREM1 may be a potential pharmacological target for colorectal cancer treatment.

1. Introduction

Colorectal cancer (CRC) is the third most common cancer worldwide and the second leading cause of cancer-related deaths [1]. CRC has accounted for 10% of malignancies and its mortality rate remains categorically high. More than half of CRC patients have metastases at the time they are diagnosed [2]. The current standard of care for CRC includes surgery, chemotherapy, radiation therapy and targeted therapy, with a 5-year relative survival rate of 65% for CRC patients, which devastatingly declines to 12% in patients with stage IV CRC [3].
Gremlin-1 (GREM1), a secreted glycoprotein, belongs to the DAN/Cerberus protein family, which is a member of the cysteine knot superfamily that includes, among others, transforming growth factor-β (TGF-β) and vascular endothelial growth factor (VEGF) [4]. It is also a BMP antagonist and known to be involved in embryogenesis [4,5], bone formation [6] and organ development [7]. In pathological states, GREM1 is involved in processes such as organ fibrosis [8], inflammation [9] and cancer [10,11]. Above all, GREM1 is known to induce fibrosis of organs, which requires the epithelial–mesenchymal transition (EMT) process [12,13]. Cancer-associated fibroblasts (CAFs) express GREM1, which inhibits BMP signaling and accelerates tumor cell proliferation [14]. In patients with breast cancer or CRC, high expression of GREM1 is usually associated with poor prognosis [15,16,17]. GREM1 secreted by glioma cancer stem cells maintains its own stemness and proliferation, whilst blocking differentiation of glioma cells [18]. Overexpression of GREM1 in intestinal epithelial cells supports the colorectal premalignant lesions [11,19]. In contrast, a clinical study showed that high expression of GREM1 was correlated with good prognosis in CRC [20]. Thus, the functional role of GREM1 in CRC remains elusive.
The endoplasmic reticulum (ER) in eukaryotic cells is involved in a variety of functions comprising lipid biosynthesis, calcium storage, protein transport and protein folding [21,22]. Stress stimuli disturb ER proteostasis [23,24], which leads to the activation of the unfolded protein response (UPR). The UPR signaling network is initiated by three ER transmembrane sensors, IRE1α, PERK and ATF6, and functions either to restore ER proteostasis or to confer apoptosis depending on the context [25,26,27]. In recent years, ER stress and UPR signaling have been extensively studied in various types of cancers, including both tumor cells per se, and the stromal cells within their microenvironment [28,29,30]. In the light of these efforts, targeting the key nodes of the UPR network has been proposed for novel therapeutic strategies in cancer treatment, and is under rapid translational development [31,32]. Despite this positive progress, relatively less is known about the implication of ER stress signaling in cancer metastasis. A recent study revealed that EMT was triggered by activation of the IRE1a and PERK pathways in lung adenocarcinoma cells A549 and H358 [33]. In pancreatic adenocarcinoma cells Capan-2 and SW1990, intracellular Ca2+ overload also triggered EMT through the IRE1a pathway [34]. However, its relevance to CRC metastasis has not been studied.
In this study, we performed immunohistochemistry (IHC) of GREM1 in CRC clinical specimens, and found that its expression was associated with poor prognosis. GREM1 activated EMT in CRC cells, which was mediated by the upregulation of ATF6 and downregulation of ATF4 pathways of the UPR. Further, our study preliminarily showed that GREM1 modulated the expression of ATF4 and ATF6 through BMP and VEGF signaling pathways, respectively. These results suggest that the crosstalk between GREM1 and the ER stress signaling network may provide a new theoretical basis for the treatment of advanced CRC.

2. Materials and Methods

2.1. Cell Lines and Culture

Human CRC cell lines (SW480 and HCT116), requested from the “Chinese Academy of Sciences”, were confirmed fingerprinting using “STR DNA” by Shanghai Biowing Applied Biotechnology Co., Ltd. (Shanghai, China). A Mycoplasma Stain Assay Kit (Beyotime, C0296) was used to identify the mycoplasma infection. The HCT116 cell line was cultivated with RPMI 1640 (Thermo Fisher Scientific, 31870082, Fremont, CA, USA) with 10% fetal bovine serum (FBS, Thermo Fisher Scientific), 100 μg/mL streptomycin and 100 U/mL penicillin (Strep/Pen, Thermo Fisher Scientific, 15140148) under the conditions of 5% CO2 at 37 °C. Similarly, SW480 cell line was cultivated with Dulbecco’s modified Eagle’s medium (DMEM, Thermo Fisher Scientific, 11965092) containing 10% FBS, 100 μg/mL streptomycin and 100 U/mL penicillin (Strep/Pen, Thermo Fisher Scientific, 15140148) in a humidified, 5% CO2 incubator at 37 °C.

2.2. Lentivirus Production and Infection

The shRNA was in the PTSB-SH-copGFP-2A-PURO lentiviral vector (shGREM1). The GREM1 gene was ligated into a pLV-EF1a-IRES-PURO lentiviral vector to construct the plasmid for GREM1 overexpression. To produce lentiviral particles, 1 × 107 HEK293T cells in a 55 cm2 dish were co-transfected with 10 μg lentiviral vector, 5 μg pCMV-dR8.91 (packaging-expressing plasmid) and 5 μg pMD2.G (envelope-expressing plasmid). The supernatant containing viral particles was harvested at 48 and 72 h post transfection, and was filtered through a Millex-GP Filter Unit (0.45 μm pore size, Millipore, SLHV033RB). To infect cancer cells with lentivirus, cells were infected twice at a 24 h interval with culture medium containing 2 mL lentivirus, 200 μL FBS and 5 mg/mL polybrene (Sigma, TR-1003, St. Louis, MO, USA) at 37 °C. To increase the overexpression and knockdown efficiency, infected cells underwent several days of puromycin (Sigma, P8833) selection. For lentiviral infection, a GREM1 expression vector (pLV-GREM1), a pLV-EF1a-IRES-PURO control lentiviral vector (pLV), the shRNA targeting GREM1 (shGREM1) and a PTSB-SH-copGFP-2A-PURO lentiviral vector (shNC) were designed and constructed by TranSheepBio Medical Biotechnology Co., Ltd. (Shanghai, China) The sequence of shRNA oligo-nucleotides against GREM1 (shGREM1) was: 5′-CCGGGCAGTGTCGTT-GCATATCCATCTCGAGATGGATATGCAACGACACTGCTTTTTT-3′. The mRNA and protein expression levels of GREM1 were then analyzed by qRT-PCR and Western blotting, respectively.

2.3. RNA-Seq and Gene Set Enrichment Analysis (GSEA)

Total RNA was isolated using the Trizol reagent (Invitrogen, 15596026, Carlsbad, CA, USA). A Bioanalyzer 2100 (Agilent Technologies, Berlin, Germany) was used to monitor the integrity of the RNA samples. RNA-seq was performed with the Illumina-PE150 sequencer (Illumina, San Diego, CA, USA). The NOISeq method was employed to identify differentially expressed genes with a 1.5-fold change between two groups (HCT116-pLV-GREM1 and HCT116-pLV cells) [35]. To clarify the molecular mechanisms underlying the effect of GREM1 on the biological behavior of colorectal cancer cells, gene set enrichment analysis (GSEA) was performed as described previously [36], and was used to explore the enrichment terminology in the pathway of the gene signature Kyoto Encyclopedia of Genes and Genomes (KEGG) related to EMT, UPR, TGFβ, VEGF and PI3K/AKT/mTOR. A heatmap was generated from this data using the R heatmap function. A p-value of <0.05 was considered statistically significant.

2.4. qRT-PCR

The total RNA of both cell lines was extracted with protocols provided by Trizol (Invitrogen, 15596026). According to the manufacturer’s instructions, cDNA was generated using the PrimeScript RT Reagent Kit with gDNA Eraser (Accurate Biology, AG11706, Hunan, China). The SYBR Green Premix Pro Taq HS qPCR Kit (Accurate Biology, AG11701) was then used to quantify mRNA expression according to the manual instructions. For these experiments, gene-specific and GAPDH primers were employed (Table A1). Data were analyzed using BioRad CFX Manager Software (BioRad, Hercules, CA, USA). Each experiment was prepared in triplicate, and data are represented as mean ± SD of at least three independent experiments. To normalize sample variation, expression of GAPDH was determined as the internal control.

2.5. Western Blot Analysis

The total lysate was prepared according to the designated processing procedure. CEB lysis buffer (Invitrogen, FNN0011) was used to harvest the cells. Quantification of proteins was performed according to the manufacturer’s procedure using Pierce’s BCA Protein Assay Reagent Kit (Pierce Biotechnology, 23227, Rockford, IL, USA). Protein bands were visualized using ECL Western blotting substrate (Thermo Fisher, 32132, Rockford, IL, USA) and the ChemiDoc Imager System (Bio-Rad), and analyzed using ImageJ software. Membranes were blocked with bovine serum albumin (BSA, Sigma–Aldrich, A1933, St. Louis, MO, USA) and incubated with one of the following primary antibodies purchased from Cell Signaling Technology (CST, MA, USA): ZEB1 (Cat. #3195, 1:1000), Snail (Cat. #3879, 1:1000), Slug (Cat. #9585, 1:1000), E-cadherin (Cat. #3195, 1:1000), Vimentin (Cat. #5741, 1:1000), ZO-1 (Cat. #8193, 1:1000), β-catenin (Cat. #8480, 1:1000), ATF4 (Cat. #11815, 1:1000), IRE1α (Cat. #3294, 1:1000), ATF6 (Cat. #65880, 1:1000), p-PI3K (Cat. #4228, 1:1000), PI3K (Cat. #4249, 1:1000), p-Akt (Cat. #4060, 1:1000), Akt (Cat. #4691, 1:1000), mTOR (Cat. #2983, 1:1000), p-mTOR (Cat. #5536, 1:1000), Smad1 (Cat. #6944, 1:1000), p-Smad1/Smad5/Smad9 (Cat. #13820, 1:1000), Bip (Cat. #3177, 1:1000) and GAPDH (Cat. #51332, 1:5000). Secondary antibodies were as follows: anti-rabbit antibody (Proteintech, SA00001-2, 1:5000, Wuhan, China) and anti-mouse antibody (Proteintech, SA00001-1, 1:5000). Pre-stained protein molecular weight marker (Thermo Scientific, 26616) was used for protein size calculation. Each experiment was performed in triplicate and the most representative images were selected. The relative band intensity was assessed by densitometric analysis of digitalized autographic images using ImageJ software. To normalize sample variation, expression of GAPDH was determined as the internal control.

2.6. Invasion Assay

Cells (1 × 105) were plated in serum-free media in the 24-well “Corning Transwell plate” that was coated with Matrigel (Corning Inc., 354480, NY, USA) for invasion experiments. The lower chamber was filled with DMEM containing 10% FBS. Movements of SW480 and HCT116 cells were measured with five random visual fields quantified by microscopy after a 36 h incubation followed by staining with 0.1% crystal violet. Invasive capacity of cells was assessed with ImageJ software.

2.7. Wound Healing Assay

Cells were inoculated into 6-well plates after centrifugation and digestion with 0.25% trypsin. When the cell density reached up to 90% or more, 3 vertical lines were scratched in each well with a 10 µL pipette tip and the floating cells were gently washed away with 1× PBS. Complete medium was added, and a photo was taken of the scratches at 0 h. Three different fields of view were selected for each well for photography. After photography, the medium was changed to a serum-free medium. The healing of the wound was photographed at the same location at the corresponding time of 24 h, 48 h and 72 h after incubation. The experiment was repeated 3 times. To observe the healing effect of drugs on the cell wound, the corresponding drugs (Tunicamycin 0.5 μg/mL, Sigma, T7765; GSK621 30 μM, MCE, HY-100548; CeapinA7 16 μM, Sigma, 2323027-38-7) were added into the serum-free medium. Afterwards, the wound healing was observed under the microscope at 0 h, 24 h, 48 h and 72 h and analyzed using ImageJ software.

2.8. Animal Experiments

Animal experimental methods were approved by the Ethical Committee of Sun Yat-sen University (Institutional Animal Care and Use KY-2021-096-02). Six-week-old nude mice were purchased from Zhejiang Vital River Laboratory Animal Technology Co., Ltd. (Hangzhou, China). Mice were kept in specific pathogen-free conditions: 20–24 °C, 12/12 h of dark/light cycle, 60 ± 5% humidity and plastic cage (four mice/cage). HCT116 cells were stably transduced with a lentiviral overexpression vector encoding luciferase according to the manufacturer’s protocol (Shanghai Genechem Co., Ltd., Shanghai, China). Here, HCT116 pLV-GREM1-luc and HCT116 pLV-luc cells (2 × 106) were intravenously injected via mouse tails. GSK621 (30 mg/kg) and CeapinA7 (470 μg/kg) were injected intraperitoneally into the mice once daily, respectively, until day 12. Lung tumor colonization was observed by bioluminescence imaging on the 5th day and 12th day.

2.9. Immunohistochemistry

Informed consent was obtained for all specimens collected. The experimental protocol of this study was approved by the Ethics Committee of the Seventh Affiliated Hospital, Sun Yat-sen University. This study abided by the Declaration of Helsinki principles. Immunohistochemical staining for GREM1 (1:100 dilution, Biorbyt, orb10741, Cambridge, UK) was performed on formalin-fixed, paraffin-embedded samples from 55 patients with clinical stage IV CRC. After dewaxing, hydration and antigen repair, the rest of the experimental procedures were performed according to the instructions of the SP Immunohistochemistry Kit (ZSBIO, PV9000, Beijing, China). Finally, after DAB staining, hematoxylin re-staining and neutral resin sealing, the sections were observed under the microscope for the degree of staining.

2.10. Statistical Analysis

All experiments were performed in triplicate. GraphPad Prism 8.0 was used for statistical analysis. Student’s t-test was used to compare the differences between two groups. Statistical differences of multiple groups were determined by one-way analysis of variance (ANOVA) followed by Tukey’s post hoc tests. The mean ± standard error of the mean (SEM) were used to express the data. A p-value of less than 0.05 was considered statistically significant. *, p < 0.05.

3. Results

3.1. The GREM1 Level Is Associated with Poor Prognosis of CRC

To investigate the clinical relevance of GREM1 in CRC, we first stained patient specimens to determine GREM1 expression via IHC. Compared with normal tissues, GREM1 was significantly upregulated in tumor tissues (Figure 1A). Moreover, RNA-seq of the CRC dataset in the Human Protein ATLAS (https://www.proteinatlas.org/ accessed on 30 April 2022) showed that GREM1 overexpression was significantly associated with poor overall survival (OS). The Kaplan–Meier plots were generated as per the best cutoff value of 5.31, suggesting that a level of GREM1 expression (FPKM) lower than 5.31 was defined as “Low”, otherwise as “High”(p = 0.046; Figure 1B). Thus, by distinguishing high vs. low expression of GREM1, we found that the corresponding survival curves between the two groups were separated significantly. These findings suggest that GREM1 may be associated with human CRC progression.

3.2. GREM1 Promotes Invasion, Migration and ER Stress of CRC Cells

To evaluate the functional role of GREM1 in CRC cells, we stably overexpressed or knocked down GREM1 in human CRC SW480 and HCT116 cells using a lentivirus-based system including a GREM1 gene vector and its control (pLV-GREM1 vs. pLV), and the shRNA targeting GREM1 and its control (shGREM1 vs. shNC). The protein expression levels of GREM1 were detected by Western blotting that showed that GREM1 protein expression levels were substantially elevated in HCT116 and SW480 cells expressing pLV-GREM1 compared to GREM1 negative control, whereas GREM1 protein expression levels were significantly reduced in shRNA-expressing lentivirus-infected cells (Figure 2A).
Firstly, to identify candidate genes that are sensitive to GREM1 overexpression, we performed RNA-seq and the subsequent GSEA analysis revealed remarkable cellular phenotypes upon GREM1 overexpression, such as enriched gene sets associated with EMT activation and unfolded protein response (UPR) (NES = 2.0395, NP = 0.0000; NES = 1.7594, NP = 0.0000, respectively) (Figure 2B). In parallel, the pertinent heatmap derived from the RNA-seq showed that expression of EMT-inducing genes, such as VEGFA [37], CXCL1 [38] and JUN [39], were significantly increased, whereas a significant decrease was found in EMT suppressor genes, such as MCM7 [40], FBLN1 [41,42] and DPYSL3 [43]. Interestingly, the majority of UPR downstream genes were significantly upregulated, except for the PERK/ATF4 arm of UPR signaling-associated genes, such as EIF2S1 (EIF-2a) and EIF2S2 (EIF-2b) [44] that were significantly downregulated in HCT116 pLV-GREM1 vs. in HCT116 pLV cells (Figure 2C).
Further, to investigate the changes in migration and invasion capacity upon GREM1 overexpression or underexpression, we performed wound healing assays which showed that GREM1-overexpressing CRC cells migrated markedly faster than control cells, in keeping with the results of the invasion assays. Meanwhile, knockdown of GREM1 by GREM1-specific shRNA reduced the invasion and migration of HCT116 and SW480 cells (Figure 2D,E). All these results indicated that GREM1 promoted migration and invasion of CRC cells. We postulated that these cellular phenotypes might be attributed to EMT and ER stress presented in our aforementioned transriptomic signature. Consistently, expression of canonical EMT factors ZEB1, Vimentin, Snail and ZO-1 were significantly upregulated, and expression of E-cadherin was markedly downregulated at both the mRNA (Figure 2F,G) and the protein levels (Figure 2H) in GREM1-overexpressing CRC cells. Previous studies have shown that EMT in a variety of tumor cells is inextricably linked to the activation of ERS, which could be seen in our study that the canonical ER stress chaperone and sensors Bip, ATF6 and IRE1a were upregulated whilst ATF4 was downregulated at the mRNA and protein expression levels in GREM1-overexpressing HCT116 and SW480 cells (Figure 2I–K). Taken together, these results indicate that GREM1 promotes EMT and regulates ER stress signaling in CRC cells.

3.3. ER Stress Activator Tunicamycin Inhibits Invasion and Migration of CRC

To further determine the sequential role of ER stress and EMT in CRC progression, we assessed the influence of ER stress on the motility of SW480 and HCT116 cells that were treated with a general ER stresser, tunicamycin (TM), and subjected to the wound healing assay. We found that the treatment group showed slower migration, which filled only half of the wound region 48 h after injury (Figure 3A). In keeping, qRT-PCR (Figure 3B,C) and Western blotting (Figure 3D) confirmed that TM treatment reduced the expression of ZEB1 and Vimentin and increased the expression of E-cadherin. As expected, TM treatment significantly increased the expression of ATF4 and Bip in these cells, although ATF6 expression was not affected to a substantial extent (Figure 3D). These results suggest that activated ATF4 signaling may be involved in inhibition of CRC motility and EMT.

3.4. GREM1 Promotes CRC Invasion and Migration through Activating ATF6 but Inhibiting the ATF4 Signaling Pathways

To further delineate the explicit arms of the UPR signaling that might be regulated by GREM1, we used our RNA-seq data to perform GSEA that presented negative regulation of PERK (NES = 1.6542, NP = 0.0000) but positive regulation of ATF6 signaling pathways (NES = 1.3553, NP = 0.0000) (Figure 4A). GSK621, a known AMPK activator, also induces PERK phosphorylation and activates the downstream eIF2α/ATF4 signaling pathway [45]. CeapinA7 is a selective blocker of ATF6α signaling in response to ER stress [46,47]. To gain insight into the underlying molecular mechanisms, we used the small molecules GSK621 and CeapinA7 to modulate the activity of these UPR pathways. GSK621 and CeapinA7 caused a significant reduction in cell invasion and migration detected by Transwell assay (Figure 4B) and wound healing assay (Figure 4C,D). Consistent with this result, the EMT-related genes (ZEB1, ZO-1 and Snail) were significantly downregulated, and E-cadherin was significantly upregulated in the GSK621 and CeapinA7 treatment group at both the mRNA (Figure 5A–D) and protein levels (Figure 5E,F), respectively, suggesting an EMT suppression role of these two compounds. These results collectively consolidate that ATF6 and ATF4 of the UPR pathways mediate the GREM1-induced CRC invasion and migration in vitro.

3.5. GREM1 Modulates ATF4 and ATF6 via Inhibiting BMP and Activating VEGF-VEGFR2-PI3K-AKT Signaling Pathways

GREM1 is known as an antagonist of BMP2 to inhibit TGFβ/BMP signaling [48]. In contrast, activation of the BMP2 signaling pathway promotes the expression of ATF4 [49]. Hence, we postulated that overexpression of GREM1 in CRC cell lines would inhibit ATF4 expression by canonically suppressing the BMP2 signaling pathway. On the other hand, GREM1 was reported to non-canonically bind VEGFR2 and activate various intracellular effectors downstream [50]. The PI3K-AKT-mTOR axis, a classical downstream pathway of VEGF [51], could promote ATF6 expression [52]. We thus hypothesized that GREM1 might activate VEGFR2 and its downstream PI3K-AKT-mTOR axis to elevate ATF6 expression. In support, GSEA of our in-house RNA-seq dataset revealed the enrichment of the PI3K/AKT/mTOR, VEGF-VEGFR2 and TGFβ signaling pathways in HCT116 pLV-GREM1 cells (Figure 6A,B). In keeping, the heatmaps derived from the RNA-seq displayed that genes associated with VEGF/VEGFR2 and PI3K/AKT/mTOR signaling pathways, such as MAP2K1 (MEK1), MAPK14 (p38), CDC42, [51,53], FOXO4 [54], LDLR [55] and HPRT1 [56], were upregulated in response to GREM1 overexpression, whilst genes of the TGFβ signaling pathway were suppressed, such as ID1, ID2 and TGFB1 [57] (Figure 6C). Next, to validate these findings, we assessed the changes in protein expression levels of the key downstream molecules involved in these pathways by overexpression or knockdown of GREM1 in CRC cells. The results showed that overexpression of GREM1 promoted the activation of the PI3K/AKT/mTOR axis, downstream of the VEGF-VEGFR2 pathway, but impeded BMP-Smad1/5/9 signaling (Figure 6D,E). On the contrary, knockdown of GREM1 resulted in declined expression of phosphorylated PI3K, phosphorylated AKT and phosphorylated mTOR but in elevated expression of phosphorylated Smad1/5/9 (Figure 6D,E).

3.6. Effects of GREM1, GSK621 and CeapinA7 on CRC Metastasis In Vivo

To further evaluate our findings in vivo, luciferase-encoding lentiviral overexpression vectors (luc) were successfully transduced into HCT116 cells. HCT116-luc cells (pLV-luc or pLV-GREM1-luc) were injected into the tail veins of nude mice and tumor progression was then monitored with bioluminescence imaging. GREM1 overexpression significantly increased tumor metastasis to the lung. By contrast, treatments by GSK621 which could activate ATF4 or by CeapinA7, an ATF6 inhibitor, significantly suppressed tumor metastasis rates (Figure 7A–C). These results further corroborate that GREM1 could promote CRC metastasis through activating the ATF6 pathway but inhibiting the ATF4 pathway (Figure 8).

4. Discussion

GREM1, a BMP antagonist, is thought to play a role in organogenesis, body patterning and tissue differentiation [48,58,59]. Canonically, Gremlin1 directly binds and inhibits BMP2 and BMP4 [48]. Furthermore, GREM1 also participates in the EMT through regulating the STAT3-MMP3 signaling, in addition to the BMP signaling pathway [60]. Non-canonically, previous studies have demonstrated that GREM1 binds to VEGFR2 and promotes angiogenesis [50]. The effect of GREM1 on the prognosis of malignant tumors remains elusive. Some studies have linked high GREM1 expression to poor prognosis in several malignancies, including ER-negative breast cancer [16], cervical cancer [61], ex-trahepatic cholangiocarcinoma [62], basal cell carcinoma [63] and renal cell carcinoma [64]. In CRC, high GREM1 protein expression has been associated with low tumor stage and extended survival [20]. However, our study showed that high GREM1 expression was significantly correlated with CRC progression.
GREM1 was found to be overexpressed in an array of common neoplasms of cervix [65], ovary [66], lung [67], stomach [68], breast [15] and kidney [64], and to interact with YWHAH protein to exert a pro-carcinogenic effect [69]. GREM1 was also overexpressed in malignant mesothelioma [70], pancreatic neuroendocrine tumors [71] and hepatocellular carcinoma associated with hepatitis C infection [72]. It was reported that high GREM1 expression in gliomas and cervical cancers played an important role in maintaining the stemness of cancer stem cells [65,73]. In our study, we found that GREM1 was highly expressed in stage IV CRC tissues and was strongly associated with poor prognosis. We also found that high GREM1 expression promoted EMT in CRC cells. There are sporadic reports of detailed insights into the molecular mechanisms of GREM1-mediated tumor-igenesis. GREM1 is associated with AKT/mTOR signaling in lung malignant mesothelioma [18]. GREM1 in breast cancer cells promote ERK activation [74]. GREM1 in Caco2 colon cancer cells inhibits differentiation by suppressing p21/CKDN1A expression [75]. Several reports have identified that GREM1 amplifies TGFβ1 signaling to drive EMT in CRC [76] and esophageal squamous cell carcinoma [77]. Our study illustrated that overexpression of GREM1 in CRC cells modulated ER stress, through which GREM1 promoted EMT via activation of the ATF6 and inhibition of the ATF4 pathways.
Protein handling, modification and folding in the endoplasmic reticulum (ER) are tightly regulated processes that determine cell function, fate and survival [30]. Studies have shown that the signaling molecules of ER stress affected tumor metastasis by reg-ulating EMT-transcription factors [78,79,80]. As Tunicamycin (TM) is widely used as an ER stress inducer in experimental settings due to its inhibitory effect on N-linked glycosylation [81,82], we used TM to treat CRC cells to observe the effect of ER stress on EMT. The UPR induced by ER stress is one of the most important mechanisms regulating cellular adaptation to an adverse microenvironment [83]. TM treatment in CRC cells profoundly suppressed the EMT, which is largely due to activation of the PERK/ATF4 arm of UPR. Since we did not observe a drastic change of ATF6 thereafter, it promoted us to investigate further how these two signature pathways orchestrate the EMT. We used small molecules to modulate the activity of these UPR pathways, namely GSK621 to induce PERK-eIF2α-ATF4 pathway [84] and CeapinA7 to block the activation of ATF6α [85]. By those approaches, our study demonstrated that GREM1 could regulate ER stress and EMT by activating ATF6 and inhibiting ATF4 signaling pathways. To further explore the intrinsic molecular mechanism underlying GREM1 regulation of ATF4 and ATF6 signaling, we conducted transcriptomic analysis and integrated it with previous studies. We observed that GREM1 inhibited TGFβ/BMP signaling, which was possibly through antagonizing BMP2 [48] and thereby suppressing ATF4 expression [49]. It was reported that GREM1 integrates with VEGFR2 to activate the PI3K/AKT/mTOR signaling pathway thereby promoting ATF6 expression [51,52]. We found that GREM1 activated PI3K/AKT/mTOR, a classical downstream of the VEGF-VEGFR2 signaling pathway that was likely to promote ATF6 expression. Taken together, these results implied that GREM1 induced EMT in CRC cells by inhibiting ATF4 expression and promoting ATF6 expression through the regulation of TGFβ/BMP and VEGF/VEGFR2, respectively. Detailed mechanisms as to how ATF6 and ATF4 divergently regulate CRC cell metastasis are worth further investigation in the future.
In conclusion, our study demonstrated that GREM1 was associated with CRC pro-gression, and it may play an important pro-metastatic role. This could be mediated, at least in part, by the divergent regulation of the ATF6 and ATF4 of the UPR pathways through activation of PI3K/AKT/mTOR and antagonization of BMP2 signaling pathways.

Author Contributions

R.L.: Investigation (equal); Writing-original draft (equal); Visualization. H.Z.: Investigation (equal); Writing—original draft (equal). M.L.: Resources. Q.M.: Data Curation. Z.F.: Methodology. Y.J.: Methodology. C.L.: Methodology. Y.G.: Methodology. Y.F.: Resources; Funding acquisition; Writing—review and editing (equal). K.W.: Resources. C.D.C.: Writing—review and editing; Supervision. X.S.: Conceptualization; Writing—review and editing (equal); Supervision. Y.H.: Methodology; Supervision. N.L.: Conceptualization; Writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (81874176 and 82072766, N.L.); Technology Planning Project of Guangdong Province, China (No. 2019B030316031); Shenzhen Science, Technology and Innovation Commission (SZSTI) Basic Research Program (JCYJ20190809154411427, N.L.); Shenzhen Sanming Project of Medicine (SZSM201911003, N.L.; SZSM202111005, Y.F.); the Hundred Talents Program of Sun Yat-sen University (392007, N.L.); Preclinical Development Program of the Seventh Affiliated Hospital of Sun Yat-sen University.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of the Seventh Affiliated Hospital of Sun Yat-sen University (protocol code: KY-2021-096-02 and date of approval: 25 November 2021).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

All data described are in this article.

Acknowledgments

The overall survival curve shown in Figure 1B was created and modified by analyzing the data from http://www.proteinatlas.org/ (last accessed on 30 April 2022). We thank Tomas Lindahl, William Richardson, Axel Behrens and Xianan Fu for their advice and guidance.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

This article contains supporting information. Table A1: Sequences of the primers for qRT-PCR.
Table A1. Sequences of the primers for qRT-PCR.
Table A1. Sequences of the primers for qRT-PCR.
GeneqRT-PCR forward Primer (5′-3′)qRT-PCR Reversed Primer (5′-3′)
GREM1CTGCTGAAGGGAAAA AGAAGATGGATATGCAACGACACT
E-cadherinATTTTTCCCTCGACACCCGATTCCCAGGCGTAGACCAAGA
ZEB-1CGAGTCAGATGCAGAAAATGAGCAAACCCAGACTGCGTCACATGTCTT
ZO-1AGGCGGATGGTGCTACAAGTGAAGAGGACCGTGTAATGGCAGACT
SnailTCGGAAGCCTAACTACAGCGAAGATGAGCATTGGCAGCGAG
VimentinAGTCCACTGAGTACCGGAGACCATTTCACGCATCTGGCGTTC
β-CateninGAGTGCTGAAGGTGCTATCTGTCTGTTCTGAACAAGACGTTGACTTGGA
Hspa5 (BiP)GACGGGCAAAGATGTCAGGAGCCCGTTTGGCCTTTTCTAC
ATF4ATGACCGAAATGAGCTTCCTGGCTGGAGAACCCATGAGGT
ATF6AGCAGCACCCAAGACTCAAACGCATAAGCGTTGGTACTGTCTGA
IRE1aAGTATGTGGAGCAGAAGGACGTTGTGTGGCTTTAGGTCTC
GAPDHGATTTGGTCGTATTGGGCGTGGAAGATGGTGATGGGAT

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  2. Poturnajova, M.; Furielova, T.; Balintova, S.; Schmidtova, S.; Kucerova, L.; Matuskova, M. Molecular features and gene expression signature of metastatic colorectal cancer. Oncol. Rep. 2021, 45, 10. [Google Scholar] [CrossRef]
  3. Miller, K.D.; Nogueira, L.; Mariotto, A.B.; Rowland, J.H.; Yabroff, K.R.; Alfano, C.M.; Jemal, A.; Kramer, J.L.; Siegel, R.L. Cancer treatment and survivorship statistics, 2019. CA Cancer J. Clin. 2019, 69, 363–385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Xiang, Q.; Hong, D.; Liao, Y.; Cao, Y.; Liu, M.; Pang, J.; Zhou, J.; Wang, G.; Yang, R.; Wang, M.; et al. Overexpression of Gremlin1 in Mesenchymal Stem Cells Improves Hindlimb Ischemia in Mice by Enhancing Cell Survival. J. Cell. Physiol. 2017, 232, 996–1007. [Google Scholar] [CrossRef] [PubMed]
  5. Zuniga, A.; Laurent, F.; Lopez-Rios, J.; Klasen, C.; Matt, N.; Zeller, R. Conserved cis-regulatory regions in a large genomic landscape control SHH and BMP-regulated Gremlin1 expression in mouse limb buds. BMC Dev. Biol. 2012, 12, 23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Gazzerro, E.; Smerdel-Ramoya, A.; Zanotti, S.; Stadmeyer, L.; Durant, D.; Economides, A.N.; Canalis, E. Conditional deletion of gremlin causes a transient increase in bone formation and bone mass. J. Biol. Chem. 2007, 282, 31549–31557. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Church, R.H.; Ali, I.; Tate, M.; Lavin, D.; Krishnakumar, A.; Kok, H.M.; Hombrebueno, J.R.; Dunne, P.D.; Bingham, V.; Goldschmeding, R.; et al. Gremlin1 plays a key role in kidney development and renal fibrosis. Am. J. Physiol. Ren. Physiol. 2017, 312, F1141–F1157. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, Y.Q.; Wan, L.Y.; He, X.M.; Ni, Y.R.; Wang, C.; Liu, C.B.; Wu, J.F. Gremlin1 Accelerates Hepatic Stellate Cell Activation Through Upregulation of TGF-Beta Expression. DNA Cell Biol. 2017, 36, 603–610. [Google Scholar] [CrossRef]
  9. Lavoz, C.; Alique, M.; Rodrigues-Diez, R.; Pato, J.; Keri, G.; Mezzano, S.; Egido, J.; Ruiz-Ortega, M. Gremlin regulates renal inflammation via the vascular endothelial growth factor receptor 2 pathway. J. Pathol. 2015, 236, 407–420. [Google Scholar] [CrossRef]
  10. Park, S.-A.; Sung, N.J.; Choi, B.-J.; Kim, W.; Kim, S.H.; Surh, Y.-J. Gremlin-1 augments the oestrogen-related receptor α signalling through EGFR activation: Implications for the progression of breast cancer. Br. J. Cancer 2020, 123, 988–999. [Google Scholar] [CrossRef]
  11. Davis, H.; Irshad, S.; Bansal, M.; Rafferty, H.; Boitsova, T.; Bardella, C.; Jaeger, E.; Lewis, A.; Freeman-Mills, L.; Giner, F.C.; et al. Aberrant epithelial GREM1 expression initiates colonic tumorigenesis from cells outside the stem cell niche. Nat. Med. 2015, 21, 62–70. [Google Scholar] [CrossRef] [PubMed]
  12. Rodrigues-Diez, R.; Rodrigues-Diez, R.R.; Lavoz, C.; Carvajal, G.; Droguett, A.; Garcia-Redondo, A.B.; Rodriguez, I.; Ortiz, A.; Egido, J.; Mezzano, S.; et al. Gremlin activates the Smad pathway linked to epithelial mesenchymal transdifferentiation in cultured tubular epithelial cells. Biomed. Res. Int. 2014, 2014, 802841. [Google Scholar] [CrossRef] [PubMed]
  13. Kalluri, R.; Neilson, E.G. Epithelial-mesenchymal transition and its implications for fibrosis. J. Clin. Investig. 2003, 112, 1776–1784. [Google Scholar] [CrossRef] [PubMed]
  14. Kobayashi, H.; Gieniec, K.A.; Wright, J.A.; Wang, T.; Asai, N.; Mizutani, Y.; Lida, T.; Ando, R.; Suzuki, N.; Lannagan, T.R.M.; et al. The Balance of Stromal BMP Signaling Mediated by GREM1 and ISLR Drives Colorectal Carcinogenesis. Gastroenterology 2021, 160, 1224–1239.e30. [Google Scholar] [CrossRef]
  15. Ren, J.; Smid, M.; Iaria, J.; Salvatori, D.C.F.; van Dam, H.; Zhu, H.J.; Martens, J.W.M.; Ten Dijke, P. Cancer-associated fibroblast-derived Gremlin 1 promotes breast cancer progression. Breast Cancer Res. 2019, 21, 109. [Google Scholar] [CrossRef] [Green Version]
  16. Neckmann, U.; Wolowczyk, C.; Hall, M.; Almaas, E.; Ren, J.; Zhao, S.; Johannessen, B.; Skotheim, R.I.; Bjorkoy, G.; Ten Dijke, P.; et al. GREM1 is associated with metastasis and predicts poor prognosis in ER-negative breast cancer patients. Cell Commun. Signal. 2019, 17, 140. [Google Scholar] [CrossRef] [Green Version]
  17. Karagiannis, G.S.; Berk, A.; Dimitromanolakis, A.; Diamandis, E.P. Enrichment map profiling of the cancer invasion front suggests regulation of colorectal cancer progression by the bone morphogenetic protein antagonist, gremlin-1. Mol. Oncol. 2013, 7, 826–839. [Google Scholar] [CrossRef] [PubMed]
  18. Yan, K.; Wu, Q.; Yan, D.H.; Lee, C.H.; Rahim, N.; Tritschler, I.; DeVecchio, J.; Kalady, M.F.; Hjelmeland, A.B.; Rich, J.N. Glioma cancer stem cells secrete Gremlin1 to promote their maintenance within the tumor hierarchy. Genes Dev. 2014, 28, 1085–1100. [Google Scholar] [CrossRef] [Green Version]
  19. Kawasaki, K.; Fujii, M.; Sugimoto, S.; Ishikawa, K.; Matano, M.; Ohta, Y.; Toshimitsu, K.; Takahashi, S.; Hosoe, N.; Sekine, S.; et al. Chromosome Engineering of Human Colon-Derived Organoids to Develop a Model of Traditional Serrated Adenoma. Gastroenterology 2020, 158, 638–651.e8. [Google Scholar] [CrossRef] [Green Version]
  20. Pelli, A.; Vayrynen, J.P.; Klintrup, K.; Makela, J.; Makinen, M.J.; Tuomisto, A.; Karttunen, T.J. Gremlin1 expression associates with serrated pathway and favourable prognosis in colorectal cancer. Histopathology 2016, 69, 831–838. [Google Scholar] [CrossRef]
  21. Wang, M.; Kaufman, R.J. The impact of the endoplasmic reticulum protein-folding environment on cancer development. Nat. Rev. Cancer 2014, 14, 581–597. [Google Scholar] [CrossRef] [PubMed]
  22. Song, M.J.; Malhi, H. The unfolded protein response and hepatic lipid metabolism in non alcoholic fatty liver disease. Pharmacol. Ther. 2019, 203, 107401. [Google Scholar] [CrossRef] [PubMed]
  23. Zhang, K.; Wang, S.; Malhotra, J.; Hassler, J.R.; Back, S.H.; Wang, G.; Chang, L.; Xu, W.; Miao, H.; Leonardi, R.; et al. The unfolded protein response transducer IRE1alpha prevents ER stress-induced hepatic steatosis. EMBO J. 2011, 30, 1357–1375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wu, J.; Chen, Y.J.; Dobbs, N.; Sakai, T.; Liou, J.; Miner, J.J.; Yan, N. STING-mediated disruption of calcium homeostasis chronically activates ER stress and primes T cell death. J. Exp. Med. 2019, 216, 867–883. [Google Scholar] [CrossRef] [PubMed]
  25. Hetz, C. The unfolded protein response: Controlling cell fate decisions under ER stress and beyond. Nat. Rev. Mol. Cell Biol. 2012, 13, 89–102. [Google Scholar] [CrossRef]
  26. Zhou, K.; Zheng, Z.; Li, Y.; Han, W.; Zhang, J.; Mao, Y.; Chen, H.; Zhang, W.; Liu, M.; Xie, L.; et al. TFE3, a potential therapeutic target for Spinal Cord Injury via augmenting autophagy flux and alleviating ER stress. Theranostics 2020, 10, 9280–9302. [Google Scholar] [CrossRef]
  27. Bhardwaj, M.; Leli, N.M.; Koumenis, C.; Amaravadi, R.K. Regulation of autophagy by canonical and non-canonical ER stress responses. Semin. Cancer Biol. 2020, 66, 116–128. [Google Scholar] [CrossRef]
  28. Urra, H.; Dufey, E.; Avril, T.; Chevet, E.; Hetz, C. Endoplasmic Reticulum Stress and the Hallmarks of Cancer. Trends Cancer 2016, 2, 252–262. [Google Scholar] [CrossRef]
  29. Wang, Y.; Wang, K.; Jin, Y.; Sheng, X. Endoplasmic reticulum proteostasis control and gastric cancer. Cancer Lett. 2019, 449, 263–271. [Google Scholar] [CrossRef]
  30. Chen, X.; Cubillos-Ruiz, J.R. Endoplasmic reticulum stress signals in the tumour and its microenvironment. Nat. Rev. Cancer 2021, 21, 71–88. [Google Scholar] [CrossRef]
  31. Zhang, T.; Li, N.; Sun, C.; Jin, Y.; Sheng, X. MYC and the unfolded protein response in cancer: Synthetic lethal partners in crime? EMBO Mol. Med. 2020, 12, e11845. [Google Scholar] [CrossRef] [PubMed]
  32. Hetz, C.; Axten, J.M.; Patterson, J.B. Pharmacological targeting of the unfolded protein response for disease intervention. Nat. Chem. Biol. 2019, 15, 764–775. [Google Scholar] [CrossRef]
  33. Shah, P.P.; Beverly, L.J. Regulation of VCP/p97 demonstrates the critical balance between cell death and epithelial-mesenchymal transition (EMT) downstream of ER stress. Oncotarget 2015, 6, 17725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Sheng, W.; Wang, G.; Tang, J.; Shi, X.; Cao, R.; Sun, J.; Lin, Y.H.; Jia, C.; Chen, C.; Zhou, J. Calreticulin promotes EMT in pancreatic cancer via mediating Ca 2+ dependent acute and chronic endoplasmic reticulum stress. J. Exp. Clin. Cancer Res. 2020, 39, 209. [Google Scholar] [CrossRef] [PubMed]
  35. Tarazona, S.; García-Alcalde, F.; Dopazo, J.; Ferrer, A.; Conesa, A. Differential expression in RNA-seq: A matter of depth. Genome Res. 2011, 21, 2213–2223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Yu, J.; Cao, Q.; Mehra, R.; Laxman, B.; Yu, J.; Tomlins, S.A.; Creighton, C.J.; Dhanasekaran, S.M.; Shen, R.; Chen, G. Integrative genomics analysis reveals silencing of β-adrenergic signaling by polycomb in prostate cancer. Cancer Cell 2007, 12, 419–431. [Google Scholar] [CrossRef] [Green Version]
  37. Kim, M.; Jang, K.; Miller, P.; Picon-Ruiz, M.; Yeasky, T.M.; El-Ashry, D.; Slingerland, J.M. VEGFA links self-renewal and metastasis by inducing Sox2 to repress miR-452, driving Slug. Oncogene 2017, 36, 5199–5211. [Google Scholar] [CrossRef]
  38. Park, Y.L.; Kim, H.P.; Ock, C.Y.; Min, D.W.; Kang, J.K.; Lim, Y.J.; Song, S.H.; Han, S.W.; Kim, T.Y. EMT-mediated regulation of CXCL1/5 for resistance to anti-EGFR therapy in colorectal cancer. Oncogene 2022, 41, 2026–2038. [Google Scholar] [CrossRef]
  39. Raychaudhuri, K.; Chaudhary, N.; Gurjar, M.; D’Souza, R.; Limzerwala, J.; Maddika, S.; Dalal, S.N. 14-3-3σ Gene Loss Leads to Activation of the Epithelial to Mesenchymal Transition Due to the Stabilization of c-Jun Protein. J. Biol. Chem. 2016, 291, 16068–16081. [Google Scholar] [CrossRef] [Green Version]
  40. Zhang, Y.; Xu, L.; Wang, P.; Jian, H.; Shi, X.; Jia, M.; Mo, L.; Hu, Z.; Li, H.; Li, J. Phenotypic transition of tumor cells between epithelial- and mesenchymal-like state during adaptation to acidosis. Cell Cycle 2019, 18, 1938–1947. [Google Scholar] [CrossRef]
  41. Harikrishnan, K.; Cooley, M.A.; Sugi, Y.; Barth, J.L.; Rasmussen, L.M.; Kern, C.B.; Argraves, K.M.; Argraves, W.S. Fibulin-1 suppresses endothelial to mesenchymal transition in the proximal outflow tract. Mech. Dev. 2015, 136, 123–132. [Google Scholar] [CrossRef] [PubMed]
  42. Harikrishnan, K.; Joshi, O.; Madangirikar, S.; Balasubramanian, N. Cell Derived Matrix Fibulin-1 Associates With Epidermal Growth Factor Receptor to Inhibit Its Activation, Localization and Function in Lung Cancer Calu-1 Cells. Front. Cell Dev. Biol. 2020, 8, 522. [Google Scholar] [CrossRef] [PubMed]
  43. Tan, F.; Wahdan-Alaswad, R.; Yan, S.; Thiele, C.J.; Li, Z. Dihydropyrimidinase-like protein 3 expression is negatively regulated by MYCN and associated with clinical outcome in neuroblastoma. Cancer Sci. 2013, 104, 1586–1592. [Google Scholar] [CrossRef] [PubMed]
  44. Artero-Castro, A.; Perez-Alea, M.; Feliciano, A.; Leal, J.A.; Genestar, M.; Castellvi, J.; Peg, V.; Ramón, Y.C.S.; Lleonart, M.E. Disruption of the ribosomal P complex leads to stress-induced autophagy. Autophagy 2015, 11, 1499–1519. [Google Scholar] [CrossRef] [Green Version]
  45. Sujobert, P.; Poulain, L.; Paubelle, E.; Zylbersztejn, F.; Grenier, A.; Lambert, M.; Townsend, E.C.; Brusq, J.-M.; Nicodeme, E.; Decrooqc, J. Co-activation of AMPK and mTORC1 induces cytotoxicity in acute myeloid leukemia. Cell Rep. 2015, 11, 1446–1457. [Google Scholar] [CrossRef]
  46. Gallagher, C.M.; Garri, C.; Cain, E.L.; Ang, K.K.-H.; Wilson, C.G.; Chen, S.; Hearn, B.R.; Jaishankar, P.; Aranda-Diaz, A.; Arkin, M.R. Ceapins are a new class of unfolded protein response inhibitors, selectively targeting the ATF6α branch. eLife 2016, 5, e11878. [Google Scholar] [CrossRef]
  47. Gallagher, C.M.; Walter, P. Ceapins inhibit ATF6α signaling by selectively preventing transport of ATF6α to the Golgi apparatus during ER stress. eLife 2016, 5, e11880. [Google Scholar] [CrossRef]
  48. Church, R.H.; Krishnakumar, A.; Urbanek, A.; Geschwindner, S.; Meneely, J.; Bianchi, A.; Basta, B.; Monaghan, S.; Elliot, C.; Strömstedt, M.; et al. Gremlin1 preferentially binds to bone morphogenetic protein-2 (BMP-2) and BMP-4 over BMP-7. Biochem. J. 2015, 466, 55–68. [Google Scholar] [CrossRef] [Green Version]
  49. Hughes, A.; Oxford, A.E.; Tawara, K.; Jorcyk, C.L.; Oxford, J.T. Endoplasmic reticulum stress and unfolded protein response in cartilage pathophysiology; contributing factors to apoptosis and osteoarthritis. Int. J. Mol. Sci. 2017, 18, 665. [Google Scholar] [CrossRef] [Green Version]
  50. Wellbrock, J.; Sheikhzadeh, S.; Bonk, V.; Oliveira-Ferrer, L.; Klaetschke, K.; Streichert, T.; Bokemeyer, C.; von Kodolitsch, Y.; Fiedler, W. Gremlin-1 Is Overexpressed in Endothelial Cells of Patients with Loeys-Dietz Syndrome Due to Dysregulation of TGF-β Signalling. Blood 2011, 118, 3269. [Google Scholar] [CrossRef]
  51. Claesson-Welsh, L.; Welsh, M. VEGFA and tumour angiogenesis. J. Intern. Med. 2013, 273, 114–127. [Google Scholar] [CrossRef] [PubMed]
  52. Song, Q.; Han, C.; Xiong, X.; He, F.; Gan, W.; Wei, S.; Liu, H.; Li, L.; Xu, H. PI3K-Akt-mTOR signal inhibition affects expression of genes related to endoplasmic reticulum stress. Genet. Mol. Res. 2016, 15, gmr.15037868. [Google Scholar] [CrossRef] [PubMed]
  53. Lamalice, L.; Houle, F.; Jourdan, G.; Huot, J. Phosphorylation of tyrosine 1214 on VEGFR2 is required for VEGF-induced activation of Cdc42 upstream of SAPK2/p38. Oncogene 2004, 23, 434–445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Zhang, Y.; Gan, B.; Liu, D.; Paik, J.H. FoxO family members in cancer. Cancer Biol. 2011, 12, 253–259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Deng, C.F.; Zhu, N.; Zhao, T.J.; Li, H.F.; Gu, J.; Liao, D.F.; Qin, L. Involvement of LDL and ox-LDL in Cancer Development and Its Therapeutical Potential. Front. Oncol. 2022, 12, 803473. [Google Scholar] [CrossRef]
  56. Wu, T.; Jiao, Z.; Li, Y.; Su, X.; Yao, F.; Peng, J.; Chen, W.; Yang, A. HPRT1 Promotes Chemoresistance in Oral Squamous Cell Carcinoma via Activating MMP1/PI3K/Akt Signaling Pathway. Cancers 2022, 14, 855. [Google Scholar] [CrossRef] [PubMed]
  57. Miyazono, K.; Miyazawa, K. Id: A target of BMP signaling. Sci. STKE 2002, 2002, pe40. [Google Scholar] [CrossRef]
  58. O’Reilly, S. Gremlin: A complex molecule regulating wound healing and fibrosis. Cell. Mol. Life Sci. 2021, 78, 7917–7923. [Google Scholar] [CrossRef]
  59. Koppens, M.A.J.; Davis, H.; Valbuena, G.N.; Mulholland, E.J.; Nasreddin, N.; Colombe, M.; Antanaviciute, A.; Biswas, S.; Friedrich, M.; Lee, L.; et al. Bone Morp.phogenetic Protein Pathway Antagonism by Grem1 Regulates Epithelial Cell Fate in Intestinal Regeneration. Gastroenterology 2021, 161, 239–254.e9. [Google Scholar] [CrossRef]
  60. Sung, N.J.; Kim, N.H.; Surh, Y.-J.; Park, S.-A. Gremlin-1 promotes metastasis of breast cancer cells by activating STAT3-MMP13 signaling pathway. Int. J. Mol. Sci. 2020, 21, 9227. [Google Scholar] [CrossRef]
  61. Sun, Q.; Qi, X.; Zhang, W.; Li, X. Knockdown of circRNA_0007534 suppresses the tumorigenesis of cervical cancer via miR-206/GREM1 axis. Cancer Cell Int. 2021, 21, 54. [Google Scholar] [CrossRef] [PubMed]
  62. Wu, Z.; Liu, R.; Miao, X.; Li, D.; Zou, Q.; Yuan, Y.; Yang, Z. Prognostic and clinicopathological significance of Hapto and Gremlin1 expression in extrahepatic cholangiocarcinoma. J. Cancer 2020, 11, 199–207. [Google Scholar] [CrossRef] [PubMed]
  63. Kim, H.S.; Shin, M.S.; Cheon, M.S.; Kim, J.W.; Lee, C.; Kim, W.H.; Kim, Y.S.; Jang, B.G. GREM1 is expressed in the cancer-associated myofibroblasts of basal cell carcinomas. PLoS ONE 2017, 12, e0174565. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. van Vlodrop, I.J.; Baldewijns, M.M.; Smits, K.M.; Schouten, L.J.; van Neste, L.; van Criekinge, W.; van Poppel, H.; Lerut, E.; Schuebel, K.E.; Ahuja, N.; et al. Prognostic significance of Gremlin1 (GREM1) promoter CpG island hypermethylation in clear cell renal cell carcinoma. Am. J. Pathol. 2010, 176, 575–584. [Google Scholar] [CrossRef] [Green Version]
  65. Sato, M.; Kawana, K.; Fujimoto, A.; Yoshida, M.; Nakamura, H.; Nishida, H.; Inoue, T.; Taguchi, A.; Takahashi, J.; Adachi, K.; et al. Clinical significance of Gremlin 1 in cervical cancer and its effects on cancer stem cell maintenance. Oncol. Rep. 2016, 35, 391–397. [Google Scholar] [CrossRef] [Green Version]
  66. Sneddon, J.B.; Zhen, H.H.; Montgomery, K.; van de Rijn, M.; Tward, A.D.; West, R.; Gladstone, H.; Chang, H.Y.; Morganroth, G.S.; Oro, A.E.; et al. Bone morphogenetic protein antagonist gremlin 1 is widely expressed by cancer-associated stromal cells and can promote tumor cell proliferation. Proc. Natl. Acad. Sci. USA 2006, 103, 14842–14847. [Google Scholar] [CrossRef] [Green Version]
  67. Mulvihill, M.S.; Kwon, Y.W.; Lee, S.; Fang, L.T.; Choi, H.; Ray, R.; Kang, H.C.; Mao, J.H.; Jablons, D.; Kim, I.J. Gremlin is overexpressed in lung adenocarcinoma and increases cell growth and proliferation in normal lung cells. PLoS ONE 2012, 7, e42264. [Google Scholar] [CrossRef]
  68. Yamasaki, Y.; Ishigami, S.; Arigami, T.; Kita, Y.; Uchikado, Y.; Kurahara, H.; Kijima, Y.; Maemura, K.; Natsugoe, S. Expression of gremlin1 in gastric cancer and its clinical significance. Med. Oncol. 2018, 35, 30. [Google Scholar] [CrossRef] [Green Version]
  69. Namkoong, H.; Shin, S.M.; Kim, H.K.; Ha, S.A.; Cho, G.W.; Hur, S.Y.; Kim, T.E.; Kim, J.W. The bone morphogenetic protein antagonist gremlin 1 is overexpressed in human cancers and interacts with YWHAH protein. BMC Cancer 2006, 6, 74. [Google Scholar] [CrossRef] [Green Version]
  70. Wang, D.J.; Zhi, X.Y.; Zhang, S.C.; Jiang, M.A.; Liu, P.; Han, X.P.; Li, J.; Chen, Z.; Wang, C.L. The bone morphogenetic protein antagonist Gremlin is overexpressed in human malignant mesothelioma. Oncol. Rep. 2012, 27, 58–64. [Google Scholar] [CrossRef]
  71. Chen, M.H.; Yeh, Y.C.; Shyr, Y.M.; Jan, Y.H.; Chao, Y.; Li, C.P.; Wang, S.E.; Tzeng, C.H.; Chang, P.M.; Liu, C.Y.; et al. Expression of gremlin 1 correlates with increased angiogenesis and progression-free survival in patients with pancreatic neuroendocrine tumors. J. Gastroenterol. 2013, 48, 101–108. [Google Scholar] [CrossRef]
  72. Guimei, M.; Baddour, N.; Elkaffash, D.; Abdou, L.; Taher, Y. Gremlin in the pathogenesis of hepatocellular carcinoma complicating chronic hepatitis C: An immunohistochemical and PCR study of human liver biopsies. BMC Res. Notes 2012, 5, 390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Li, C.; Rezov, V.; Joensuu, E.; Vartiainen, V.; Ronty, M.; Yin, M.; Myllarniemi, M.; Koli, K. Pirfenidone decreases mesothelioma cell proliferation and migration via inhibition of ERK and AKT and regulates mesothelioma tumor microenvironment in vivo. Sci. Rep. 2018, 8, 10070. [Google Scholar] [CrossRef] [PubMed]
  74. Sung, N.J.; Kim, N.H.; Bae, N.Y.; Jo, H.S.; Park, S.A. DHA inhibits Gremlin-1-induced epithelial-to-mesenchymal transition via ERK suppression in human breast cancer cells. Biosci. Rep. 2020, 40, BSR20200164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Kosinski, C.; Li, V.S.W.; Chan, A.S.Y.; Zhang, J.; Ho, C.; Tsui, W.Y.; Chan, T.L.; Mifflin, R.C.; Powell, D.W.; Yuen, S.T.; et al. Gene expression patterns of human colon tops and basal crypts and BMP antagonists as intestinal stem cell niche factors. Proc. Natl. Acad. Sci. USA 2007, 104, 15418–15423. [Google Scholar] [CrossRef] [Green Version]
  76. Karagiannis, G.S.; Afaloniati, H.; Karamanavi, E.; Poutahidis, T.; Angelopoulou, K. BMP pathway suppression is an early event in inflammation-driven colon neoplasmatogenesis of uPA-deficient mice. Tumor Biol. 2016, 37, 2243–2255. [Google Scholar] [CrossRef]
  77. Hong, D.; Liu, T.; Huang, W.; Liao, Y.; Wang, L.; Zhang, Z.; Chen, H.; Zhang, X.; Xiang, Q. Gremlin1 Delivered by Mesenchymal Stromal Cells Promoted Epithelial-Mesenchymal Transition in Human Esophageal Squamous Cell Carcinoma. Cell. Physiol. Biochem. 2018, 47, 1785–1799. [Google Scholar] [CrossRef]
  78. Santamaria, P.G.; Mazon, M.J.; Eraso, P.; Portillo, F. UPR: An Upstream Signal to EMT Induction in Cancer. J. Clin. Med. 2019, 8, 624. [Google Scholar] [CrossRef] [Green Version]
  79. Han, C.C.; Wan, F.S. New Insights into the Role of Endoplasmic Reticulum Stress in Breast Cancer Metastasis. J. Breast Cancer 2018, 21, 354–362. [Google Scholar] [CrossRef]
  80. Li, H.; Chen, X.; Gao, Y.; Wu, J.; Zeng, F.; Song, F. XBP1 induces snail expression to promote epithelial- to-mesenchymal transition and invasion of breast cancer cells. Cell Signal. 2015, 27, 82–89. [Google Scholar] [CrossRef]
  81. Azim, M.; Surani, H.J.C. Glycoprotein synthesis and inhibition of glycosylation by tunicamycin in preimplantation mouse embryos: Compaction and trophoblast adhesion. Cell 1979, 18, 217–227. [Google Scholar] [CrossRef]
  82. Wu, J.; Chen, S.; Liu, H.; Zhang, Z.; Ni, Z.; Chen, J.; Yang, Z.; Nie, Y.; Fan, D. Tunicamycin specifically aggravates ER stress and overcomes chemoresistance in multidrug-resistant gastric cancer cells by inhibiting N-glycosylation. J. Exp. Clin. Cancer Res. 2018, 37, 272. [Google Scholar] [CrossRef] [PubMed]
  83. Hetz, C.; Zhang, K.; Kaufman, R.J. Mechanisms, regulation and functions of the unfolded protein response. Nat. Rev. Mol. Cell Biol. 2020, 21, 421–438. [Google Scholar] [CrossRef] [PubMed]
  84. Grenier, A.; Poulain, L.; Mondesir, J.; Jacquel, A.; Bosc, C.; Stuani, L.; Mouche, S.; Larrue, C.; Sahal, A.; Birsen, R.; et al. AMPK-PERK axis represses oxidative metabolism and enhances apoptotic priming of mitochondria in acute myeloid leukemia. Cell Rep. 2022, 38, 110197. [Google Scholar] [CrossRef] [PubMed]
  85. Ma, M.; Li, H.; Wang, P.; Yang, W.; Mi, R.; Zhuang, J.; Jiang, Y.; Lu, Y.; Shen, X.; Wu, Y.; et al. ATF6 aggravates angiogenesis-osteogenesis coupling during ankylosing spondylitis by mediating FGF2 expression in chondrocytes. iScience 2021, 24, 102791. [Google Scholar] [CrossRef] [PubMed]
Figure 1. GREM1 expression is correlated with poor survival in CRC patients. (A) Representative images of GREM1 in para-carcinoma and carcinoma tissues of patients with stage IV CRC. Scale bar, 250 µm. (B) Kaplan–Meier survival analysis of 597 colon cancer patients stratified by GREM1 expression levels. Kaplan–Meier survival analysis and GREM1 mRNA expression according to FPKM value 5.31 as a cutoff value in the whole cohort. p = 0.046 by log-rank (Mantel–Cox) test.
Figure 1. GREM1 expression is correlated with poor survival in CRC patients. (A) Representative images of GREM1 in para-carcinoma and carcinoma tissues of patients with stage IV CRC. Scale bar, 250 µm. (B) Kaplan–Meier survival analysis of 597 colon cancer patients stratified by GREM1 expression levels. Kaplan–Meier survival analysis and GREM1 mRNA expression according to FPKM value 5.31 as a cutoff value in the whole cohort. p = 0.046 by log-rank (Mantel–Cox) test.
Cells 11 02136 g001
Figure 2. GREM1 expression promotes cell invasion and migration of HCT116 and SW480 cells and ER stress. (A) Overexpression or underexpression of GREM1 was detected by Western blotting. (B) GSEA was performed using our RNA-seq dataset to identify sets of hallmark genes that were positively associated with GREM1 expression. Enrichment plots showed EMT (NES = 2.0395, NP = 0.0000) and UPR (NES = 1.7594, NP = 0.0000). (C) The heatmap showed signature genes involving EMT and UPR. (D,E) Wound healing test and Boyden chambers with Matrigel revealed that GREM1 overexpression promoted HCT116 and SW480 cell migration and invasion. (FH) Cells were analyzed after harvesting, for protein and mRNA expression. qRT-PCR and Western blotting analysis showed expression of mRNA and proteins involved in EMT. (IK) qRT-PCR and Western blotting results showed the expression of proteins and mRNA that were involved in UPR and ER stress.
Figure 2. GREM1 expression promotes cell invasion and migration of HCT116 and SW480 cells and ER stress. (A) Overexpression or underexpression of GREM1 was detected by Western blotting. (B) GSEA was performed using our RNA-seq dataset to identify sets of hallmark genes that were positively associated with GREM1 expression. Enrichment plots showed EMT (NES = 2.0395, NP = 0.0000) and UPR (NES = 1.7594, NP = 0.0000). (C) The heatmap showed signature genes involving EMT and UPR. (D,E) Wound healing test and Boyden chambers with Matrigel revealed that GREM1 overexpression promoted HCT116 and SW480 cell migration and invasion. (FH) Cells were analyzed after harvesting, for protein and mRNA expression. qRT-PCR and Western blotting analysis showed expression of mRNA and proteins involved in EMT. (IK) qRT-PCR and Western blotting results showed the expression of proteins and mRNA that were involved in UPR and ER stress.
Cells 11 02136 g002
Figure 3. Tunicamycin inhibits epithelial-to-mesenchymal transition in human CRC cell lines. SW480 and HCT116 cells were exposed to tunicamycin (0.1 μg/mL, 0.5 μg/mL). (A) Wound healing test revealed that tunicamycin inhibited HCT116 and SW480 cell migration. (BD) SW480 and HCT116 cells were exposed to tunicamycin (0.1 μg/mL, 0.5 μg/mL) for 6 h, and then analyzed for protein and mRNA expression. qRT-PCR and Western blotting results showed the expression of proteins and mRNA involved in EMT.
Figure 3. Tunicamycin inhibits epithelial-to-mesenchymal transition in human CRC cell lines. SW480 and HCT116 cells were exposed to tunicamycin (0.1 μg/mL, 0.5 μg/mL). (A) Wound healing test revealed that tunicamycin inhibited HCT116 and SW480 cell migration. (BD) SW480 and HCT116 cells were exposed to tunicamycin (0.1 μg/mL, 0.5 μg/mL) for 6 h, and then analyzed for protein and mRNA expression. qRT-PCR and Western blotting results showed the expression of proteins and mRNA involved in EMT.
Cells 11 02136 g003
Figure 4. CeapinA7 and GSK621 inhibit GREM1-induced CRC cellline invasion and migration. (A) The GSEA was performed using an in-house RNA-seq dataset to identify gene sets in the KEGG collections that have positive correlations with GREM1 expression. The negative regulating PERK-mediated UPR (NES = 1.6542 NP = 0.0000) and ATF6-mediated UPR signaling pathway (NES = 1.3553, NP = 0.0000) with GREM1 overexpression in HCT116 cells. (B) HCT116 pLV-GREM1 and SW480 pLV-GREM1 cells were seeded into Matrigel-coated inserts with medium containing CeapinA7 (16 μM) or GSK621 (30 μM) and incubated for 36 h, followed by the invasion assay. (C,D) Wound healing test revealed that CeapinA7 and GSK621 inhibited cell migration.
Figure 4. CeapinA7 and GSK621 inhibit GREM1-induced CRC cellline invasion and migration. (A) The GSEA was performed using an in-house RNA-seq dataset to identify gene sets in the KEGG collections that have positive correlations with GREM1 expression. The negative regulating PERK-mediated UPR (NES = 1.6542 NP = 0.0000) and ATF6-mediated UPR signaling pathway (NES = 1.3553, NP = 0.0000) with GREM1 overexpression in HCT116 cells. (B) HCT116 pLV-GREM1 and SW480 pLV-GREM1 cells were seeded into Matrigel-coated inserts with medium containing CeapinA7 (16 μM) or GSK621 (30 μM) and incubated for 36 h, followed by the invasion assay. (C,D) Wound healing test revealed that CeapinA7 and GSK621 inhibited cell migration.
Cells 11 02136 g004
Figure 5. CeapinA7 and GSK621 reverse EMT in GREM1-overexpressing CRC cell lines. HCT116 pLV-GREM1 and SW480 pLV-GREM1 cells were exposed to CeapinA7 (6 μM,16 μM) or GSK621 (30 μM) for 24 h, respectively. After harvesting, the cells were analyzed for protein and mRNA expression. (AD) qRT-PCR analysis revealed the expression of mRNA involved in EMT. (E,F) Western blotting analysis revealed the expression of proteins involved in the EMT.
Figure 5. CeapinA7 and GSK621 reverse EMT in GREM1-overexpressing CRC cell lines. HCT116 pLV-GREM1 and SW480 pLV-GREM1 cells were exposed to CeapinA7 (6 μM,16 μM) or GSK621 (30 μM) for 24 h, respectively. After harvesting, the cells were analyzed for protein and mRNA expression. (AD) qRT-PCR analysis revealed the expression of mRNA involved in EMT. (E,F) Western blotting analysis revealed the expression of proteins involved in the EMT.
Cells 11 02136 g005
Figure 6. GREM1 negatively regulates ATF4 and but positively regulates ATF6 via suppressing BMP and activating the VEGF-PI3K-AKT signaling pathways. (A,B) GSEA of RNA-seq data showed enriched gene sets associated with VEGF/VEGFR2, PI3K/Akt/mTOR and TGFβ signaling pathways annotated in the KEGG by comparing GREM1 overexpression versus negative control CRC cells. (C) Heatmap showing gene expression differences in the pertinent signaling pathways. (D,E) Western blotting analysis revealed the expression of proteins involved in the pertinent signaling pathways.
Figure 6. GREM1 negatively regulates ATF4 and but positively regulates ATF6 via suppressing BMP and activating the VEGF-PI3K-AKT signaling pathways. (A,B) GSEA of RNA-seq data showed enriched gene sets associated with VEGF/VEGFR2, PI3K/Akt/mTOR and TGFβ signaling pathways annotated in the KEGG by comparing GREM1 overexpression versus negative control CRC cells. (C) Heatmap showing gene expression differences in the pertinent signaling pathways. (D,E) Western blotting analysis revealed the expression of proteins involved in the pertinent signaling pathways.
Cells 11 02136 g006
Figure 7. Anti-tumor activity of CeapinA7 and GSK621 in vivo. (A) Animal experimental grouping and plans. (B) After the injection of different cells (2 × 10 6 cells/nude mouse) via tail vein of nude mice, the distribution of CRC cells in these mice was detected via a bioluminescence imaging (BLI) system. The scale to the right of the BLI images describes the color map for the luminescent signal. (C) Fluorescence intensity of lung metastasis areas in different groups of the CRC cells. n = 21. * p < 0.05, ** p < 0.01.
Figure 7. Anti-tumor activity of CeapinA7 and GSK621 in vivo. (A) Animal experimental grouping and plans. (B) After the injection of different cells (2 × 10 6 cells/nude mouse) via tail vein of nude mice, the distribution of CRC cells in these mice was detected via a bioluminescence imaging (BLI) system. The scale to the right of the BLI images describes the color map for the luminescent signal. (C) Fluorescence intensity of lung metastasis areas in different groups of the CRC cells. n = 21. * p < 0.05, ** p < 0.01.
Cells 11 02136 g007
Figure 8. Mechanistic diagram of the role of GREM1 in ER stress, EMT and metastasis in CRC cells.
Figure 8. Mechanistic diagram of the role of GREM1 in ER stress, EMT and metastasis in CRC cells.
Cells 11 02136 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, R.; Zhou, H.; Li, M.; Mai, Q.; Fu, Z.; Jiang, Y.; Li, C.; Gao, Y.; Fan, Y.; Wu, K.; et al. Gremlin-1 Promotes Colorectal Cancer Cell Metastasis by Activating ATF6 and Inhibiting ATF4 Pathways. Cells 2022, 11, 2136. https://doi.org/10.3390/cells11142136

AMA Style

Li R, Zhou H, Li M, Mai Q, Fu Z, Jiang Y, Li C, Gao Y, Fan Y, Wu K, et al. Gremlin-1 Promotes Colorectal Cancer Cell Metastasis by Activating ATF6 and Inhibiting ATF4 Pathways. Cells. 2022; 11(14):2136. https://doi.org/10.3390/cells11142136

Chicago/Turabian Style

Li, Ruohan, Huaixiang Zhou, Mingzhe Li, Qiuyan Mai, Zhang Fu, Youheng Jiang, Changxue Li, Yunfei Gao, Yunping Fan, Kaiming Wu, and et al. 2022. "Gremlin-1 Promotes Colorectal Cancer Cell Metastasis by Activating ATF6 and Inhibiting ATF4 Pathways" Cells 11, no. 14: 2136. https://doi.org/10.3390/cells11142136

APA Style

Li, R., Zhou, H., Li, M., Mai, Q., Fu, Z., Jiang, Y., Li, C., Gao, Y., Fan, Y., Wu, K., Costa, C. D., Sheng, X., He, Y., & Li, N. (2022). Gremlin-1 Promotes Colorectal Cancer Cell Metastasis by Activating ATF6 and Inhibiting ATF4 Pathways. Cells, 11(14), 2136. https://doi.org/10.3390/cells11142136

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop