- Open Access
HOXA10 mediates epithelial-mesenchymal transition to promote gastric cancer metastasis partly via modulation of TGFB2/Smad/METTL3 signaling axis
Journal of Experimental & Clinical Cancer Research volume 40, Article number: 62 (2021)
Homeobox A10 (HOXA10) belongs to the HOX gene family, which plays an essential role in embryonic development and tumor progression. We previously demonstrated that HOXA10 was significantly upregulated in gastric cancer (GC) and promoted GC cell proliferation. This study was designed to investigate the role of HOXA10 in GC metastasis and explore the underlying mechanism.
Immunohistochemistry (IHC) was used to evaluate the expression of HOXA10 in GC. In vitro cell migration and invasion assays as well as in vivo mice metastatic models were utilized to investigate the effects of HOXA10 on GC metastasis. GSEA, western blot, qRT-PCR and confocal immunofluorescence experiments preliminarily analyzed the relationship between HOXA10 and EMT. ChIP-qPCR, dual-luciferase reporter (DLR), co-immunoprecipitation (CoIP), colorimetric m6A assay and mice lung metastasis rescue models were performed to explore the mechanism by which HOXA10 accelerated the EMT process in GC.
In this study, we demonstrated HOXA10 was upregulated in GC patients and the difference was even more pronounced in patients with lymph node metastasis (LNM) than without. Functionally, HOXA10 promoted migration and invasion of GC cells in vitro and accelerated lung metastasis in vivo. EMT was an important mechanism responsible for HOXA10-involved metastasis. Mechanistically, we revealed HOXA10 enriched in the TGFB2 promoter region, promoted transcription, increased secretion, thus triggered the activation of TGFβ/Smad signaling with subsequent enhancement of Smad2/3 nuclear expression. Moreover, HOXA10 upregulation elevated m6A level and METTL3 expression in GC cells possible by regulating the TGFB2/Smad pathway. CoIP and ChIP-qPCR experiments demonstrated that Smad proteins played an important role in mediating METTL3 expression. Furthermore, we found HOXA10 and METTL3 were clinically relevant, and METTL3 was responsible for the HOXA10-mediated EMT process by performing rescue experiments with western blot and in vivo mice lung metastatic models.
Our findings indicated the essential role of the HOXA10/TGFB2/Smad/METTL3 signaling axis in GC progression and metastasis.
Gastric cancer (GC) is the fifth most common cancer and the third leading cause of cancer-related mortality worldwide due to aggressive properties and high metastasis potentials [1,2,3]. Most GC patients are first diagnosed in the metastasized stage, making the overall survival rate still low [4,5,6,7]. Therefore, it is imperative and urgent to elucidate the metastasis mechanism of GC and develop more efficient treatment strategies.
Homeobox (HOX) genes play pivotal roles in cell proliferation, apoptosis, metabolism, and migration . HOXA10 is a member of the HOX gene family. Previous studies have shown that upregulated HOXA10 promoted proliferation, inhibited GC cells’ apoptosis and indicated poor prognosis for GC patients [9, 10]. However, the underlying molecular mechanism of HOXA10-mediated invasion and metastasis in GC has not been fully elucidated.
Epithelial-mesenchymal transition (EMT) is mainly executed by EMT transcription factors (EMT-TFs) and it has been widely recognized as a crucial regulator of tumor cell invasion and metastasis to distant organs . Several EMT-TFs, such as Snail and Slug, could repress the E-cadherin expression, which is the fundamental event in the EMT process [12,13,14].
Notably, multiple studies have recognized that the transforming growth factor beta (TGFβ) signaling pathway plays a crucial role in promoting the EMT process in GC [12, 13, 15,16,17,18]. TGFβ pathway could be activated by TGF-βs (TGFβ1, TGFβ2 and TGFβ3), which are secreted by cells and have similar but not identical functions . Interestingly, recent research has revealed that TGFB2 is upregulated in GC and predicted poor survival . Besides, a previous study in myeloid progenitor cells showed HOXA10 could activate TGFB2 transcription and increase TGFβ2 production . So far, however, the relationship between HOXA10 and TGFβ signaling in GC has not been clearly understood, which prompted us to perform further exploration.
TGF-βs could activate canonical TGFβ signaling and subsequently phosphorylate Smad2/3, which are the main intracellular effectors of the TGFβ pathway [13, 21, 22]. It is noteworthy that Alessandro Bertero et al. gave us new insights that Smad2/3 could interact with the METTL3-METTL14-WTAP complex in hESCs and hiPSCs, thus promoting N6-methyladenosine (m6A) modification on specific Smad2/3 transcriptional targets .
M6A is the most common mRNA modification in eukaryotes and regulates mRNA stability, splicing, translation and degradation [24, 25]. METTL3, METTL14 and WTAP are the main m6A “writers,” of which METTL3 is the primary catalytic subunit and dysregulated METTL3 could result in changes of the total m6A methylation [26,27,28].
Recent studies have shown elevated METTL3 promotes EMT in cancers [29,30,31], and several EMT-TFs (such as Snail and Slug) are key molecules of EMT [12, 13]. It should also be noted that Lin et al. demonstrated m6A modification on Snail, which acted as a Smad2/3 target, promoted its translation in cancer cells and accelerated the EMT progression [32, 33]. Therefore, we mainly focused on the relationship between HOXA10 and METTL3 to explore whether HOXA10 might promote the EMT process and accelerate GC metastasis via activating TGFβ/Smad signaling and regulating m6A modification.
In this study, we demonstrated that HOXA10 promoted metastasis of GC by inducing EMT. The relevant mechanism could be explained as follows: elevated HOXA10 enhanced the transcription and secretion of TGFB2, activated TGFβ/Smad signaling, promoted the interaction of Smad2/3 with METTL3 in the nucleus, and upregulated the expression of Snail and Slug through m6A modification.
Materials and methods
Patients and specimens
A total of 127 paired GC tissue specimens collected for immunohistochemistry and other 80 paired GC tissue samples selected for qRT-PCR were obtained from Shanghai General Hospital for this research. This study was approved by the Ethics Committee of Shanghai General Hospital and all patients provided informed consent.
Immunohistochemistry (IHC) staining
The immunohistochemistry (IHC) analysis was performed using a Dako K5007 kit (Dako, Denmark) according to the manual instructions. The staining intensity was scored as (0 = negative, 1 = weak, 2 = moderate, 3 = strong) according to the expression level and the staining area was scored as (0 = 0%, 1 = 1–25%, 2 = 26–50%, 3 = 51–100%) based on the percentage of positive cells. The immunostaining score (IS) was calculated by adding the score of staining intensity and area according to the method of Yuan et al. . The positive specimen for HOXA10 and E-cadherin individually was defined as an IS≥3, while negative staining with an IS < 3. The primary antibodies used in immunohistochemistry included: HOXA10 (1:100, Santa Cruz), E-cadherin (1:400, CST), N-cadherin (1:125, CST), Vimentin (1:200, CST), METTL3 (1:500, Abcam).
Cell culture, lentivirus construction and transfection
Stable HOXA10 knockdown (BGC823-shHOXA10) or overexpression (AGS-HOXA10) GC cells have been successfully established . The cells were divided into the categories below: AGS-Mock, without any treatment; AGS-Ctrl, infected with the lentivirus containing the control fragment; AGS-HOXA10, infected with the lentivirus containing the HOXA10 fragment; BGC823-Mock, without any treatment; BGC823-Ctrl, infected with the control Lenti-shRNA;BGC823-shHOXA10, infected with the HOXA10 Lenti-shRNA. The cells were cultured using an RPMI1640 medium containing 10% fetal bovine serum (Gibco, USA).
The lentivirus containing METTL3 shRNA (sh-METTL3) or negative control shRNA (sh-NC) and the METTL3 overexpression lentivirus (METTL3) or control vector lentivirus (Vector) were obtained from Obio Technology (Shanghai). Stable METTL3 overexpression or knockdown in previously selected HOXA10 differential expressing GC cells was achieved by transfecting the lentivirus.
Cell scratch wound-healing assay
The cells were plated in six-well plates and cultured until full confluence. Then suitable micropipette tips were used to generate uniformed scratches from the center of each well. After washed with phosphate-buffered saline (PBS), the cells were then cultured in fetal bovine serum (FBS)-free medium and taken photos at 0 h, 48 h. The cell migration was measured by comparing the gap change with the initial gap at 0 h. The experiment was performed independently in triplicate.
Cell transwell migration and invasion assays
The transwell 24-well chambers with 8.0 μm pore size (Corning, USA) were used for cell migration and invasion assays. The cells dispersed in 200 μL serum-free medium (2 × 104/ well) were seeded into the upper chamber, and 600 μL medium containing 10% FBS was added to the lower section. The upper compartment was coated with Matrigel (BD Bioscience, USA) for invasion assay or without for migration assay. The cells migrated or invaded to the lower chamber at time point 24 h or 36 h respectively were stained with 0.1% crystal violet and then counted. Each experiment was conducted independently in triplicate.
Nude mice metastasis models
4-week male athymic BALB/c nude mice were used to establish in vivo metastasis mice models. For the lung metastasis models, 5 × 106 different stable transfected cells were injected through the mice tail vein (n = 5 for each group). Luciferase imaging was performed to monitor metastatic progression with an IVIS-100 system (Caliper Life Sciences, USA) after injection of D-Luciferin potassium salt (Synchem, Germany) intraperitoneally.
For the liver metastasis models, 5 × 106 cells (BGC823-Ctrl or BGC823-shHOXA10) resuspended in sterile PBS were injected into the spleen of nude mice.
All the nude mice were killed after 6 weeks, and all the lungs and livers were collected for HE staining or immunohistochemistry analysis. All the animal studies were performed strictly following the Animal Care Guidelines approved by the Animal Care Committee of Shanghai General Hospital.
RNA-sequencing and GSEA
For RNA-Sequencing, total RNA separated from BGC823-Ctrl and BGC-shHOXA10 were extracted using Trizol reagent (Takara, Japan). Majorbio Bio-pharm Technology Corporation performed the construction of the transcriptome libraries and Illumina sequencing, and each group contained three biological replicates. Gene Set Enrichment Analysis (GSEA) was conducted using the GSEA software downloaded from the Broad Institute website (https://www.gsea-msigdb.org/gsea/index.jsp).
Western blot analysis
Total cellular proteins or nuclear proteins were separately extracted using RIPA buffer (Beyotime Biotechnology, China) or Minute™ Cytoplasmic and Nuclear Extraction Kit for Cells (Invent Biotechnologies, USA). In brief, the protein samples were separated using SDS-PAGE gel (7.5–12%) electrophoresis and then transferred to PVDF membranes. After blockade in 5% skimmed milk, the membranes were incubated with primary antibodies at 4 °C overnight followed by incubation of specific second antibodies for 1 h at room temperature. The target proteins were then detected using ECL reagent kits (Millipore, USA). The primary antibodies used in western blot were listed in Table 1.
Quantitative real-time PCR (qRT-PCR)
For qRT-PCR, total RNA extraction from all cell lines and GC tissues was performed using Eastep® Super Total RNA Extraction Kit (Promega, USA) according to the instructions. RNA was reversely transcribed to cDNA using PrimeScript™ RT Master Mix (Takara, Japan). qRT-PCR reactions were carried out using the Hieff™ qPCR SYBR Green Master Mix (Yeasen, China).
The following primer sequences used for qRT-PCR were listed in Table 2:
HOXA10; E-cadherin; N-cadherin; Vimentin; Snail; Slug; TGFB1; TGFB2; TGFB3; METTL3; GAPDH.
The relative mRNA quantification was calculated using the 2−ΔΔCt method and normalized by GAPDH. All the experiments were repeated in triplicate.
Confocal immunofluorescence staining
Cells were seeded on confocal dishes (Cellvis, USA) for 24 h before fixed with paraformaldehyde and permeabilized by Triton X-100. After blocking with 5% normal donkey serum, the cells were incubated at 4 °C overnight with primary antibodies specific against E-cadherin (1:200, CST), N-cadherin (1:200, CST), Vimentin (1:100, CST), Smad2/3 (1:400, CST), followed by fluorescence-conjugated secondary antibodies for 1 h at room temperature. The cells were then counterstained with DAPI (Sigma, Germany) and were photographed using a confocal laser scanning microscope (Leica, Germany).
Chromatin immunoprecipitation (ChIP)
ChIP assay was performed using the SimpleChIP® Plus Enzymatic Chromatin IP Kit (Magnetic Beads) (CST, USA) following the same protocol as previously described . For the ChIP-qPCR assay, primers explicitly designed to target the TGFB2 or METTL3 promoter region were listed in Table 3. The products of ChIP-qPCR were qualitatively analyzed by 3% agarose gel electrophoresis.
Dual-luciferase reporter gene assay
For the dual-luciferase reporter assay, pGL3-TGFB2 (TGFB2-WT) and pGL3-mutated TGFB2 (TGFB2-mut) plasmids were constructed and synthesized by Zorin (Zorin, China). Besides, Flag-tagged HOXA10 overexpressing plasmid and control plasmid were reconstructed and acquired from GeneChem (GeneChem, China). 293T, AGS and BGC823 cells were co-transfected with the luciferase reporter plasmid and other indicated expression constructs using Lipofectamine 2000 (Invitrogen, USA).
The luciferase assays were performed after transfection for 48 h by using a Dual-Luciferase Reporter Assay System kit (Promega, USA). Then the relative luciferase activity was detected and normalized to Renilla activity.
TCA protein precipitation assay
The TCA Protein Precipitation Kit (Sangon Biotech, China) was used to separate protein pellets from cellular supernatant in different HOXA10 expression cells according to the manufacturers’ instructions. The protein extracts from TCA assays were then used for western blot detection.
TGFβ2 concentration in cellular supernatant from altered HOXA10 expression cells was quantified with a human TGFβ2 ELISA kit (MultiSciences, China) according to the manual guides. All the samples were repeated in triplicate.
For each immunoprecipitation assay (IP), nuclear proteins extracted from 2 × 15 cm dishes BGC823 cells with full confluence were used according to the manufacturer’s protocol (Invent Biotechnologies, USA). The nuclear protein lysates were immunoprecipitated with specific antibodies targeting against Smad2/3 (CST, USA), METTL3 (Abcam, USA) or IgG negative control (CST, USA) at 4 °C overnight in the rotator, then added into 50 μL Protein A/G Magnetic Beads (MedChemExpress, USA) pre-cleaned with PBST to rotate for 2 h at 4 °C. Subsequently, the magnetic beads were washed 4 times with PBST buffer, and the immunoprecipitation complexes were separated from the beads for further western blot analysis.
Colorimetric m6A assay
The relative N6-methyladenosine (m6A) levels in GC cells with different HOXA10 expression were measured via the EpiQuik™ m6A RNA Methylation Quantification Kit (Colorimetric) (Epigentek, USA) according to the manufacture’s guidance.
All data were analyzed with IBM SPSS Statistics version 23.0 software (SPSS, USA). The statistical significance between HOXA10 or E-cadherin protein levels with the clinicopathologic characteristics levels was measured using the Chi-square test.
The correlation between HOXA10 and E-cadherin protein was calculated with the Spearman correlation coefficient test, while the correlation between HOXA10 and METTL3 mRNA expression was analyzed using Pearson correlation analysis.
Quantitative results for in vitro and in vivo experiments were shown as mean ± SD or box plot, and the differences were determined using the two-tailed Student’s t-test. For all studies, P-values < 0.05 indicates statistically significant differences.
HOXA10 overexpression is associated with metastasis in GC
We previously demonstrated that HOXA10 was upregulated in GC tissues, and the Kaplan-Meier Plotter database showed higher expression of HOXA10 predicted poor prognosis in GC patients .
Herein to further investigate the role of HOXA10 in GC progression, we employed 127 pairs of GC samples to analyze the expression levels of HOXA10. The results showed that HOXA10 was significantly upregulated in the GC tumor samples than the adjacent normal mucosae (mean IS = 3.98 vs. 1.50, P < 0.001, Fig.1a, b). More importantly, we found HOXA10 had higher expression in GC tumor samples with LNM than those of non-LNM (mean IS = 4.45 vs. 3.19, P < 0.001; Table 4 and Fig. 1c). Besides, we observed that HOXA10 showed positive staining in metastatic lymph nodes (Fig. 1d). Taken together, these results suggested HOXA10 upregulation might be relevant to metastasis of GC.
HOXA10 promotes GC cells migratory and invasive capacities in vitro
Given that HOXA10 was upregulated in GC tissues and was associated with metastasis, we then investigated whether HOXA10 could influence the GC cell metastatic capacity in vitro.
We have established the HOXA10 overexpressing tool cell (AGS-HOXA10) and HOXA10 knockdown tool cell (BGC823-shHOXA10) . Wound-healing assay verified that AGS-HOXA10 cells closed the wound area significantly faster while BGC823-shHOXA10 cells relatively slower than respective control cells. (Fig. 2a-b). Moreover, the transwell assays indicated that overexpression of HOXA10 significantly promoted the migratory and invasive abilities of AGS cells. On the contrary, HOXA10 downregulation obviously reduced the abilities in BGC823 cells as well (Fig. 2c-d). These results suggested that HOXA10 promoted GC cells metastatic abilities in vitro.
HOXA10 enhances GC cells metastasis in vivo
To investigate the in vivo effect of HOXA10 on GC cells metastasis, we performed two metastatic assays. After AGS-HOXA10 and control cells were injected into nude mice via the tail vein, we collected images at different time points to dynamically observe the location and enlargement of lung metastases. The result showed that HOXA10 overexpression remarkably increased the number of metastatic lung nodules (Fig. 3a-b). Also, the nude mice liver metastasis model was established via injecting modified BGC823 cells into the spleens. Notably, HOXA10 silencing obviously reduced the number of metastatic liver tumors (Fig. 3c-d). The above experiments demonstrated that HOXA10 enhanced GC cells metastasis in vivo.
HOXA10 promotes GC cells metastasis by modulating EMT
To better understand the role of HOXA10 in GC metastasis, transcriptome sequencing analysis (RNA-Seq) was used to compare the relative changes of gene expression in BGC823-Ctrl and BGC823-shHOXA10 cells. Then Gene Set Enrichment Analysis (GSEA) demonstrated that HOXA10 knockdown was negatively correlated with EMT, which meant downregulated HOXA10 inhibited EMT pathway in GC (Fig. S1).
EMT is a crucial mechanism that leads to cancer metastasis, and the loss of E-cadherin is commonly recognized as a molecular marker in EMT progression . Then qRT-PCR, Western blot and confocal immunofluorescence analysis were performed to validate further whether EMT was responsible for HOXA10 induced metastatic phenotype in GC.
As shown in Fig. 4a, c, compared with the AGS-Ctrl cells, AGS-HOXA10 cells expressed a lower level of epithelial marker E-cadherin and a higher level of the mesenchymal markers N-cadherin, Vimentin as well as EMT-TFs Snail and Slug. In contrast, knockdown of HOXA10 in BGC823 cells elevated the expression of E-cadherin, which was accompanied by the suppression of N-cadherin, Vimentin, Snail and Slug (Fig. 4b, d). Moreover, the changes of characterized EMT markers in different HOXA10-expressing cells were further validated by confocal immunofluorescence (Fig. 4e).
Subsequently, we performed IHC in typical GC serial sections with or without LNM and found that HOXA10 expression was consistent with N-cadherin and Vimentin but seemed to be opposite to E-cadherin (Fig. 4f). Furthermore, we analyzed the expression correlation between HOXA10 and E-cadherin in 127 paired GC tumor samples and found that HOXA10 was inversely correlated with E-cadherin (r = − 0.469, P < 0.001; Table 5).
Besides, the relationships between clinicopathologic features and the expression of HOXA10 or E-cadherin in GC patients were summarized in Table 4. Interestingly, both proteins were significantly related to LNM (P < 0.01). In summary, these findings suggested that HOXA10 might enhance GC cells metastasis through EMT.
HOXA10 modulates TGFβ signaling pathway possibly by regulating TGFB2 in GC
TGFβ signaling has been recognized as essential in inducing GC EMT progression . To further study whether HOXA10 promoted EMT and GC metastasis through the canonical TGFβ pathway, we first treated AGS-HOXA10 cells with TGFβ/Smad signaling inhibitor HY-N0439 (MedChemExpress, USA). As shown in Fig. 5a, HY-N0439 upregulated the expression of E-cadherin while downregulated Snail and Slug. This result suggested that HOXA10-induced EMT progress might be suppressed by blocking TGFβ/Smad signaling.
Then we investigated which isoform of TGFβ played a role in HOXA10 mediating EMT. It should be noted that qPCR and western blot results showed only TGFB2 expression, but not TGFB1 or TGFB3, was significantly increased by upregulating HOXA10 (Fig. 5b, c).
Moreover, AGS-HOXA10 cells were further processed with selective TGFβ inhibitor (HY-N0439), TGFβ2 inhibitor (HY-B0673) (MedChemExpress) or TGFβ1 inhibitor (HY-P0118) (MedChemExpress). The results showed cell migration/invasion induced by HOXA10 overexpression was significantly inhibited by 47.6%/31.9% or 48.2%/30.8% after being treated with HY-N0439 or HY-B0673, while only 13.4%/7.7% with HY-P0118 (Fig. 5d; Fig. S2). Similarly, cell migration/invasion ability was partly restored in BGC823-shHOXA10 cells after cytokine hTGFβ2 treatment (Fig. 5e; Fig. S2).
Additionally, western blot analysis demonstrated the decrease of Snail and Slug accompanied by the elevation of E-cadherin when TGFβ2 activity was inhibited using HY-B0673 in AGS-HOXA10 cells. In contrast, the opposite trend was observed in BGC823-shHOXA10 cells when stimulated with hTGFβ2 (Fig. 5f).
The above results indicated that HOXA10 was involved in modulating TGFβ signaling possibly by regulating TGFB2 expression.
HOXA10 activates TGFB2 expression through transcriptional regulation
To further verify whether HOXA10 took part in the transcriptional regulation of TGFB2, we conducted the ChIP-qPCR analysis. According to the prediction of the database JASPAR (http://jaspar.genereg.net/), we selected four putative HOXA10 binding sites (P1-P4) in the TGFB2 promoter, designed and synthesized the corresponding primers. ChIP-qPCR results showed that HOXA10 had a high enrichment in the P3 area (Fig. 6a), which was qualitatively confirmed by 3% agarose gel electrophoresis (Fig. S3).
Subsequently, we constructed wild-type (WT) and mutant type (Mut) TGFB2 promoter-luciferase reporter plasmids and performed dual-luciferase reporter (DLR) gene assays (Fig. S3). As shown in Fig. 6b-d, HOXA10 overexpression significantly elevated TGFB2 promoter activity after the co-transfection of wild-type TGFB2. Nevertheless, luciferase activity remained unchanged when co-transfected with mutant TGFB2 plasmid.
Following up, we performed TCA precipitation and ELISA assays with cell supernatants to investigate whether HOXA10 could increase the secretion of TGFβ2 into extracellular space. The results showed that TGFβ2 content in the cell supernatants increased with HOXA10 overexpression, which indicated that HOXA10 could promote the secretion of TGFβ2 (Fig. 6e, f).
These results above confirmed TGFB2 as a target gene transcriptionally regulated by HOXA10.
HOXA10 increases Smad2/3 expression in the cell nucleus and promotes the deposition of METTL3
Elevated phosphorylated Smad2/3 protein and increased Smad2/3 expression in the nucleus are considered a symbol of the TGFβ signaling activation.
Therefore, to investigate whether HOXA10 overexpression activated the TGFβ/Smad pathway, we detected Smad2, Smad3, P-Smad2, P-Smad3, and nuclear Smad2/3 by western blot. As illustrated in Fig. 7a, increased expression of P-Smad2, P-Smad3 and nuclear Smad2/3 were detected in AGS-HOXA10 cells. Besides, confocal immunofluorescence showed that Smad2/3 markedly elevated in the nucleus with HOXA10 overexpression (Fig. 7b). These findings suggested that overexpression of HOXA10 activated the TGFβ/Smad signaling pathway in GC.
Moreover, as shown in Fig. 7a, the protein expressions of P-Smad2, P-Smad3 and nuclear Smad2/3 were reduced when AGS-HOXA10 cells were treated with TGFB2 inhibitor (HY-B0673), which indicated that HOXA10 could not activate the Smad pathway after TGFB2 depletion.
Notably, Alessandro Bertero et al. developed an optimized CoIP protocol and identified 89 putative interacting proteins of Smad2/3 in given cell types. The study focused on the relationship between Smad2/3 and the m6A methyltransferase complex and uncovered that Smad2/3 could interact with the METTL3-METTL14-WTAP complex .
Then we decided to explore whether HOXA10 was involved in m6A modification. The colorimetric m6A assay results demonstrated that the m6A levels were remarkably higher in AGS-HOXA10 cells, which meant HOXA10 overexpression elevated the m6A level in GC cells (Fig. 7c).
According to the previously reported researches, METTL3 was the main m6A catalytic subunit and influenced the total m6A methylation . Therefore, we intended to evaluate whether the METTL3 expression in GC cells would vary with the change of HOXA10 expression. Western blot showed that METTL3 expression increased in AGS-HOXA10 cells but decreased in BGC823-shHOXA10 cells. Besides, we also found METTL3 expression could be inhibited by TGFB2 inhibitor or enhanced by hTGFβ2, which suggested TGFB2/Smad pathway exerted important functions in HOXA10-mediated METTL3 expression (Fig. 7d).
Subsequently, we performed CoIP and ChIP-qPCR experiments to explore the relationship between Smad2/3 and METTL3. The CoIP results demonstrated that Smad2/3 could directly interact with METTL3 in the nucleus of BGC823 cells (Fig. 7e).
Moreover, to further investigate whether Smad2/3 could also participate in the transcriptional regulation of METTL3, we selected three putative Smad2/3 binding sites (P1-P3) in the METTL3 promoter according to JASPAR and performed the ChIP experiment. ChIP-qPCR results showed that Smad2/3 had a relatively high enrichment in the P1 area of the METTL3 promoter region (Fig. 7f).
Furthermore, we used ITD-1 (Selleck, USA), which blocked P-Smad2/3 induced by TGFβ2, to investigate whether suppressing P-Smad2/3 could reduce the expression of METTL3. As shown in Fig. 7g, the protein expressions of P-Smad2/3 and METTL3 in AGS-HOXA10 cells were both inhibited after being treated with ITD-1. The above two experiments showed that Smad proteins played an important role in mediating METTL3 expression.
These results gave us valuable hints that HOXA10 might further promote METTL3 expression and elevate m6A level in GC cells following the activation of the TGFβ/Smad signaling pathway.
HOXA10 regulates EMT progression partly through modulating METTL3
To identify the clinical relevance of HOXA10 and METTL3, the other 80 pairs of GC samples were used for qPCR experiments. qPCR results showed HOXA10 and METTL3 mRNA expression were both upregulated in the GC samples (HOXA10, 66/80, 82.5%; METTL3, 64/80, 80%) and the Pearson correlation coefficient of r = 0.433 (P < 0.001) (Fig. 8a-b).
Furthermore, we determined to explore whether METTL3 could participate in the HOXA10-mediated EMT process in GC cells. Firstly, AGS-Ctrl and AGS-HOXA10 cells were transfected with lentiviruses carrying sh-METTL3 or sh-NC. Besides, BGC823-Ctrl and BGC823-shHOXA10 cells were transfected with lentiviruses carrying METTL3 or Vector. As illustrated in Fig. 8c, knockdown of METTL3 partly elevated E-cadherin levels but reduced the expression of N-cadherin, Vimentin, Snail and Slug in AGS-HOXA10 cells. In contrast, overexpression of METTL3, to a certain extent, weakened the inhibition of the EMT process caused by HOXA10 knockdown (Fig. 8d).
More importantly, we demonstrated the rescue effects of METTL3 for HOXA10-induced EMT in nude mice lung metastasis models. The increased metastatic potential induced by HOXA10 overexpression could be partially weakened by METTL3 inhibition, whereas the suppressive effect of HOXA10 knockdown on lung metastasis could be in part reversed by METTL3 upregulation (Fig. S4).
Moreover, immunohistochemistry results showed that the protein level of METTL3 was upregulated with increased lung metastasis burden or decreased with fewer lung nodules caused by altered HOXA10 expression (Fig. 8e-f). Taken together, these findings demonstrated that the upregulation of METTL3 was partly responsible for HOXA10-mediated EMT progression in GC (Fig. 9).
In recent years, researches concerning the relation between HOXA10 and tumor progression have been gradually increasing [35,36,37]. However, the underlying metastatic mechanism of HOXA10 in GC remains poorly elucidated. In the present study, we have provided experimental evidence that HOXA10 might accelerate the EMT progression in GC by activating TGFB2/Smad/METTL3 signaling, thus promoting the invasion and metastasis of GC cells.
Immunohistochemistry revealed that the expression of HOXA10 was elevated in GC, and it was noteworthy that the IS score was even higher in samples with LNM. By performing in vitro cell assays and in vivo mice metastasis models, we further observed HOXA10 upregulation significantly promoted migratory, invasive and metastatic abilities of GC cells. These findings suggested that HOXA10 might play an essential role in accelerating GC progression and metastasis.
To further elucidate the potential metastatic mechanism, RNA-Seq and GSEA were performed and the results indicated that HOXA10 knockdown inhibited the EMT pathway. EMT was thought to be a crucial contributor to promote invasion and metastasis in cancer cells. Snail (also known as SNAI1) and Slug (also known as SNAI2) are important EMT-TFs. They are both members of the zinc-fingerE-box binding homeobox factor and could repress the E-cadherin expression, which is a crucial event in triggering the EMT process [12, 13, 38].
Following the GSEA indication, cell assays showed that HOXA10 knockdown increased the E-cadherin levels but reduced the expression of Snail and Slug. Moreover, in GC tumor samples, HOXA10 was negatively correlated with E-cadherin (r = − 0.469, P < 0.001). These results above supported that HOXA10 promoted the EMT process in GC and triggered metastasis.
EMT progression could be activated by multiple intracellular signaling pathways, including the TGFβ, Wnt, Notch, etc. . Accumulating research has recognized that TGFβ signaling plays a vital role in inducing GC EMT progression [15,16,17, 19].
To explore whether HOXA10 acted in promoting EMT through TGFβ signaling, we treated AGS-HOXA10 cells with TGFβ/Smad inhibitor (HY-N0439). Then we observed an elevation of E-cadherin, suggesting that HOXA10 might regulate the TGFβ pathway to promote EMT.
Three TGF-βs (TGFβ1, TGFβ2 and TGFβ3) have been identified, which could activate TGFβ signaling by interacting with specific receptors in the cell-surface . Previous studies demonstrated that higher expression of TGFβs indicated poor survival in cancer patients [19, 39]. To determine which subtype of TGF-βs was the main downstream effector of HOXA10, qRT-PCR and western blot were performed. The results showed that only TGFB2 expression was significantly increased or reduced by upregulating or inhibiting HOXA10.
Then we observed that the change of cell migration or invasion abilities induced by HOXA10 were partly negatively and positively affected with TGFβ2 inhibitors and hTGFβ2 cytokine. Therefore, we further convinced that TGFB2, among three TGF-β subtypes, might be an essential effector by which HOXA10 exert functions.
To further understand the mechanism of HOXA10-regulated TGFB2 gene expression, we conducted ChIP experiments and DLR assays. ChIP-qPCR and DLR assays results indicated that HOXA10 was capable of binding to the promoter region of TGFB2 and could enhance TGFB2 transcription activity, suggesting that HOXA10 could transcriptionally upregulate TGFB2 expression.
However, in order to activate TGFβ signaling pathway, TGFβ2 needed to be secreted into the extracellular space and then bonded to the TGFβ receptor on the cell surface . Consequently, we adopted TCA precipitation and ELISA experiments to detect extracellular TGFβ2 levels. The results indicated that HOXA10 upregulation could increase the secretion of TGFβ2 in extracellular supernatants.
The TGFβ signaling pathway exerted its effects through phosphorylating Smad2/3 and then translocating them into the cell nucleus . As illustrated in Fig. 7a-b, the protein expression of P-Smad2, P-Smad3 and nuclear Smad2/3 elevated with HOXA10 overexpression, which suggested that HOXA10 activated the TGFβ pathway.
Notably, suppressed expression of P-Smad2, P-Smad3 and nuclear Smad2/3 were observed after AGS-HOXA10 cells were treated with TGFB2 inhibitor (HY-B0673), which indicated that TGFB2 was a key molecule for HOXA10 in activating the Smad pathway.
A recent research showed that Smad2/3 could interact with 89 potential interactomes in given cell types. Most notably, the study validated the interaction of Smad2/3 with METTL3-METTL14-WTAP complex and Smad2/3 facilitated m6A modification onto nascent transcripts .
m6A modification is the most common mRNA modification mainly mediated by m6A “writers,” “erasers,” and “readers.” It has been associated with tumor progression in glioblastoma, liver cancer, breast cancer, cervical cancer, etc. [25, 33, 41]. METTL3, METTL14 and WTAP are the main partners of “writers” belonging to m6A modulators, of which METTL3 is the primary catalytic subunit and has proved essential in promoting metastasis in several cancers [26,27,28, 32].
Most importantly, recent studies have demonstrated that METTL3 enhanced the EMT process by affecting Snail expression in HCC [32, 33]. As critical members of zinc finger transcription factors belonging to the Snail family, Snail and Slug are critical downstream effectors of Smad2/3 and they accelerate EMT by repressing E-cadherin expression [12,13,14, 21, 22].
Therefore, we performed the colorimetric m6A assay and found the m6A levels of AGS-HOXA10 cells were obviously higher than that in AGS-Ctrl cells. Then we investigated and observed that the changes of METTL3 expression were consistent with the altered expression of HOXA10. Besides, the upregulation of METTL3 in AGS-HOXA10 could be suppressed by HY-B0673, the TGFβ2 inhibitor, suggesting the TGFB2/Smad pathway played an important role in mediating METTL3 expression.
Afterward, we performed CoIP experiments and validated that Smad2/3 and METTL3 could be coprecipitated with each other, which was in accord with the previous results by Alessandro Bertero et al. Interestingly, we also found that Smad2/3 might bound to the promoter of METTL3 by conducting the ChIP-qPCR experiment. Moreover, after AGS-HOXA10 cells were treated with ITD-1, the expressions of METTL3 and P-Smad2/3 were both reduced. In brief, the above experiments showed that Smad proteins exerted an important effect in modulating METTL3 expression.
Subsequently, we tested the correlation between HOXA10 and METTL3 in other 80 paired GC samples and the results showed evident correlation (r = 0.433, P < 0.001). More importantly, western blot showed that the expressions of important EMT-TFs, Snail and Slug, were reduced, and E-cadherin was upregulated when METTL3 was inhibited in AGS-HOXA10 cells. Furthermore, the in vivo lung metastatic models demonstrated that METTL3 knockdown weakened the metastatic potentials induced by HOXA10. These experiments suggested that METTL3 was an important molecule involved in HOXA10-mediated EMT.
In summary, our study preliminarily revealed that elevated HOXA10 expression in GC was significantly correlated with LNM. Overexpression of HOXA10 played a pivotal role in promoting EMT and GC metastasis by modulating TGFB2/Smad/METTL3 signaling axis, thus might help provide a novel candidate for GC treatment and develop new therapeutic strategies for GC metastasis.
Availability of data and materials
All data generated and analyzed in this study are included in this article and its supplementary information files.
Transforming growth factor beta
Gene Set Enrichment Analysis
Lymph node metastasis
Fetal bovine serum
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Huang L-Y, Zhao J, Chen H, Wan L, Inuzuka H, Guo J, et al. SCFFBW7-mediated degradation of Brg1 suppresses gastric cancer metastasis. Nat Commun. 2018;9(1):3569.
Schuhmacher C, Gretschel S, Lordick F, Reichardt P, Hohenberger W, Eisenberger CF, et al. Neoadjuvant chemotherapy compared with surgery alone for locally advanced cancer of the stomach and cardia: European Organisation for Research and Treatment of Cancer randomized trial 40954. J Clin Oncol. 2010;28(35):5210–8.
Coburn N, Cosby R, Klein L, Knight G, Malthaner R, Mamazza J, et al. Staging and surgical approaches in gastric cancer: a systematic review. Cancer Treat Rev. 2018;63:104–15.
Van Cutsem E, Sagaert X, Topal B, Haustermans K, Prenen H. Gastric cancer. Lancet. 2016;388(10060):2654–64.
Fontana E, Smyth EC. Novel targets in the treatment of advanced gastric cancer: a perspective review. Ther Adv Med Oncol. 2016;8(2):113–25.
Bernards N, Creemers GJ, Nieuwenhuijzen GAP, Bosscha K, Pruijt JFM, Lemmens VEPP. No improvement in median survival for patients with metastatic gastric cancer despite increased use of chemotherapy. Ann Oncol. 2013;24(12):3056–60.
Li B, Huang Q, Wei G-H. The role of HOX transcription factors in Cancer predisposition and progression. Cancers (Basel). 2019;11(4):528.
Song C, Han Y, Luo H, Qin Z, Chen Z, Liu Y, et al. HOXA10 induces BCL2 expression, inhibits apoptosis, and promotes cell proliferation in gastric cancer. Cancer Med. 2019;8(12):5651–61.
Lim JY, Yoon SO, Seol S-Y, Hong SW, Kim JW, Choi SH, et al. Overexpression of miR-196b and HOXA10 characterize a poor-prognosis gastric cancer subtype. World J Gastroenterol. 2013;19(41):7078–88.
Santamaria PG, Moreno-Bueno G, Portillo F, Cano A. EMT: Present and future in clinical oncology. Mol Oncol. 2017;11(7):718–38.
Dongre A, Weinberg RA. New insights into the mechanisms of epithelial-mesenchymal transition and implications for cancer. Nat Rev Mol Cell Biol. 2019;20(2):69–84.
Colak S, Ten Dijke P. Targeting TGF-β signaling in Cancer. Trends Cancer. 2017;3(1):56–71.
Thiery JP, Acloque H, Huang RYJ, Nieto MA. Epithelial-mesenchymal transitions in development and disease. Cell. 2009;139(5):871–90.
Xu G, Chen Y, Fu M, Zang X, Cang M, Niu Y, et al. Circular RNA CCDC66 promotes gastric cancer progression by regulating c-Myc and TGF-β signaling pathways. J Cancer. 2020;11(10):2759–68.
Xiong R, Gao J-L, Yin T. G3BP1 activates the TGF-β/Smad signaling pathway to promote gastric cancer. Onco Targets Ther. 2019;12:7149–56.
Jeong M-H, Park S-Y, Lee S-H, Seo JS, Yoo J-Y, Park S-H, et al. EPB41L5 mediates TGFβ-induced metastasis of gastric Cancer. Clin Cancer Res. 2019;25(12):3617–29.
Sang Cheul O, Sohn BH, Cheong J-H, Kim S-B, Lee JE, Park KC, et al. Clinical and genomic landscape of gastric cancer with a mesenchymal phenotype. Nat Commun. 2018;9(1):1777.
Yang BW, Bai J, Shi R, Shao X, Yang Y, Jin Y, et al. TGFB2 serves as a link between epithelial-mesenchymal transition and tumor mutation burden in gastric cancer. Int Immunopharmacol. 2020;84:106532.
Shah CA, Wang H, Bei L, Platanias LC, Eklund EA. HoxA10 regulates transcription of the gene encoding transforming growth factor beta2 (TGFbeta2) in myeloid cells. J Biol Chem. 2011;286(4):3161–76.
Lin Y-T, Kou-Juey W. Epigenetic regulation of epithelial-mesenchymal transition: focusing on hypoxia and TGF-β signaling. J Biomed Sci. 2020;27(1):39.
Angela Nieto M, Huang RY-J, Jackson RA, Thiery JP. EMT: 2016. Cell. 2016;166(1):21–45.
Bertero A, Brown S, Madrigal P, Osnato A, Ortmann D, Yiangou L, et al. The SMAD2/3 interactome reveals that TGFβ controls m6A mRNA methylation in pluripotency. Nature. 2018;555(7695):256–9.
Chen M, Wei L, Law C-T, Tsang FH-C, Shen J, Cheng CL-H, et al. RNA N6-methyladenosine methyltransferase-like 3 promotes liver cancer progression through YTHDF2-dependent posttranscriptional silencing of SOCS2. Hepatology. 2018;67(6):2254–70.
Lin S, Choe J, Peng D, Triboulet R, Gregory RI. The m(6) A methyltransferase METTL3 promotes translation in human Cancer cells. Mol Cell. 2016;62(3):335–45.
Wang Q, Geng W, Guo H, Wang Z, Xu K, Chen C, et al. Emerging role of RNA methyltransferase METTL3 in gastrointestinal cancer. J Hematol Oncol. 2020;13(1):57.
Wang H, Hu X, Huang M, Liu J, Gu Y, Ma L, et al. Mettl3-mediated mRNA m6A methylation promotes dendritic cell activation. Nat Commun. 2019;10(1):1898.
Chen Y, Peng C, Chen J, Chen D, Yang B, He B, et al. WTAP facilitates progression of hepatocellular carcinoma via m6A-HuR-dependent epigenetic silencing of ETS1. Mol Cancer. 2019;18(1):127.
Yu X, Zhao H, Cao Z. The m6A methyltransferase METTL3 aggravates the progression of nasopharyngeal carcinoma through inducing EMT by m6A-modified snail mRNA. Minerva Med. 2020. https://doi.org/10.23736/S0026-4806.20.06653-7. Online ahead of print.
Hua W, Zhao Y, Jin X, Yu D, He J, Xie D, et al. METTL3 promotes ovarian carcinoma growth and invasion through the regulation of AXL translation and epithelial to mesenchymal transition. Gynecol Oncol. 2018;151(2):356–65.
Liu T, Yang S, Sui J, Xu S-Y, Cheng Y-P, Shen B, et al. Dysregulated N6-methyladenosine methylation writer METTL3 contributes to the proliferation and migration of gastric cancer. J Cell Physiol. 2020;235(1):548–62.
Xu H, Wang H, Zhao W, Sirui F, Li Y, Ni W, et al. SUMO1 modification of methyltransferase-like 3 promotes tumor progression via regulating snail mRNA homeostasis in hepatocellular carcinoma. Theranostics. 2020;10(13):5671–86.
Lin X, Chai G, Wu Y, Li J, Chen F, Liu J, et al. RNA m6A methylation regulates the epithelial mesenchymal transition of cancer cells and translation of snail. Nat Commun. 2019;10(1):2065.
Yuan C, Zhu X, Han Y, Song C, Liu C, Lu S, et al. Elevated HOXA1 expression correlates with accelerated tumor cell proliferation and poor prognosis in gastric cancer partly via cyclin D1. J Exp Clin Cancer Res. 2016;35:15.
Cui Y-P, Xie M, Pan W-X, Zhang Z-Y, Li W-F. HOXA10 promotes the development of bladder cancer through regulating FOSL1. Eur Rev Med Pharmacol Sci. 2020;24(6):2945–54.
Zhang Y, Chen J, Wu S-S, Lv M-J, Yu Y-S, Tang Z-H, et al. HOXA10 knockdown inhibits proliferation, induces cell cycle arrest and apoptosis in hepatocellular carcinoma cells through HDAC1. Cancer Manag Res. 2019;11:7065–76.
Chen W, Wu G, Zhu Y, Zhang W, Zhang H, Zhou Y, et al. HOXA10 deteriorates gastric cancer through activating JAK1/STAT3 signaling pathway. Cancer Manag Res. 2019;11:6625–35.
Bure IV, Nemtsova MV, Zaletaev DV. Roles of E-cadherin and Noncoding RNAs in the Epithelial-mesenchymal Transition and Progression in Gastric Cancer. Int J Mol Sci. 2019;20(12):2870.
Jackstadt R, van Hooff SR, Leach JD, Cortes-Lavaud X, Lohuis JO, Ridgway RA, et al. Epithelial NOTCH Signaling Rewires the Tumor Microenvironment of Colorectal Cancer to Drive Poor-Prognosis Subtypes and Metastasis. Cancer Cell. 2019;36(3):319–36.e7.
Zhang Z, Deng Y, Zheng G, Jia X, Xiong Y, Luo K, et al. SRGN-TGFβ2 regulatory loop confers invasion and metastasis in triple-negative breast cancer. Oncogenesis. 2017;6(7):e360.
Lan Q, Liu PY, Haase J, Bell JL, Hüttelmaier S, Liu T. The critical role of RNA m6A methylation in Cancer. Cancer Res. 2019;79(7):1285–92.
This study was funded by grants from the National Natural Science Foundation of China (No.81672346).
Ethics approval and consent to participate
Every human sample was obtained with informed consent and was approved by the Ethics Committee of Shanghai General Hospital. All the animal studies were performed strictly in accordance with the Animal Care Guidelines approved by the Animal Care Committee of Shanghai General Hospital.
Consent for publication
The authors declare that they have no conflict of interest.
GSEA comparing BGC823-Ctrl cells with BGC823-shHOXA10 cells: the level of HOXA10 mRNA was positively correlated with Hallmark Epithelial-Mesenchymal Transition.
The statistical graph of the migration and invasion cells in Fig. 5d and e. **P < 0.01, ***P < 0.001.
Qualitative analysis of ChIP-qPCR results with 3% agarose gel electrophoresis and a schematic diagram of the TGFB2-luciferase mutation sites.
In vivo metastasis models of AGS or BGC823 cells by suppressing or elevating METTL3 expression, respectively. **P < 0.01, ***P < 0.001.
About this article
Cite this article
Song, C., Zhou, C. HOXA10 mediates epithelial-mesenchymal transition to promote gastric cancer metastasis partly via modulation of TGFB2/Smad/METTL3 signaling axis. J Exp Clin Cancer Res 40, 62 (2021). https://doi.org/10.1186/s13046-021-01859-0
- Gastric cancer