Forkhead Box q1 promotes invasion and metastasis in colorectal cancer by activating the epidermal growth factor receptor pathway.

BACKGROUND: Colorectal cancer (CRC) is an extremely malignant tumor with a high mortality rate. Little is known about the mechanism by which forkhead Box q1 (FOXQ1) causes CRC invasion and metastasis through the epidermal growth factor receptor (EGFR) pathway. AIM: To illuminate the mechanism by which FOXQ1 promotes the invasion and metastasis of CRC by activating the heparin binding epidermal growth factor (HB-EGF)/EGFR pathway.
METHODS: We investigated the differential expression and prognosis of FOXQ1 and HB-EGF in CRC using the Gene Expression Profiling Interactive Analysis (GEPIA) website (http://gepia.cancer-pku.cn/index.html). Quantitative real-time polymerase chain reaction (qRT-PCR) and Western blotting were used to detect the expression of FOXQ1 and HB-EGF in cell lines and tissues, and we constructed a stable low-expressing FOXQ1 cell line and verified it with the above method. The expression changes of membrane-bound HB-EGF (proHB-EGF) and soluble HB-EGF (sHB-EGF) in the low-expressing FOXQ1 cell line were detected by flow cytometry and ELISA. Western blotting was used to detect changes in the expression levels of HB-EGF and EGFR pathway-related downstream genes when exogenous recombinant human HB-EGF was added to FOXQ1 knockdown cells. Proliferation experiments, transwell migration experiments, and scratch experiments were carried out to determine the mechanism by which FOXQ1 activates the EGFR signaling pathway through HB-EGF, and then to evaluate the clinical relevance of FOXQ1 and HB-EGF.
RESULTS: GEPIA showed that the expression of FOXQ1 in CRC tissues was relatively high and was related to a lower overall survival rate. PCR array results showed that FOXQ1 is related to the HB-EGF and EGFR pathways. Knockdown of FOXQ1 suppressed the expression of HB-EGF, and led to a decrease in EGFR and its downstream genes AKT, RAF, KRAS expression levels. After knockdown of FOXQ1 in CRC cell lines, cell proliferation, migration and invasion were attenuated. Adding HB-EGF restored the migration and invasion ability of CRC, but not the cell proliferation ability. Kaplan-Meier survival analysis results showed that the combination of FOXQ1 and HB-EGF may serve to predict CRC survival.
CONCLUSION: Based on these collective data, we propose that FOXQ1 promotes the invasion and metastasis of CRC via the HB-EGF/EGFR pathway. ©The Author(s) 2022. Published by Baishideng Publishing Group Inc. All rights reserved.

Core Tip: Invasion and metastasis play important roles in tumorigenesis, resulting in the death of most colorectal cancer (CRC) patients. Forkhead Box q1 (FOXQ1) is a well-established oncogene in multiple tumors, including CRC. Our previous study suggested that FOXQ1 positively regulates the expression of heparin-binding epidermal growth factor (HB-EGF) and triggers the activation of the epidermal growth factor receptor (EGFR) pathway in CRC. However, the role and mechanism of how FOXQ1 promotes tumorigenesis in CRC by activating the HB-EGF/EGFR pathway remain unexplored. In the present study, our findings demonstrated that the essential role of FOXQ1-induced invasion and metastasis in CRC was related to activation of the HB-EGF/EGFR pathway.

INTRODUCTION

Colorectal cancer (CRC) is the third most common cancer in the world and the fourth most common cause of cancer death in the world[1]. Among all cancerous processes involved, local invasion and metastasis are the main factors related to cancer. The metastatic dissemination of primary tumors is directly related to the survival rate of patients, accounting for approximately 90% of all colon cancer deaths[2]. The median survival of metastatic CRC is less than 2 years[3]. Therefore, clarifying the mechanism of invasion and metastasis is the key to improving the survival rate of patients with CRC.

Forkhead Box q1 (FOXQ1) is a member of the forkhead transcription factor family[4], and it promotes tumorigenesis by activating cell proliferation, invasion and apoptosis[5]. Kaneda et al[6] found that compared to adjacent tissues, FOXQ1 is overexpressed in CRC. Overexpression of FOXQ1 reduces cell proliferation but increases cell tumorigenicity and tumor growth, and it inhibits apoptosis and promotes angiogenesis, thereby promoting CRC tumorigenesis. Liu et al[7] also demonstrated that the expression of FOXQ1 in either CRC tissue samples or cancer cell lines is higher than that in normal colorectal tissues and cell lines, and they reported that FOXQ1 promotes cancer metastasis by regulating PI3K/AKT signaling. Weng et al[8] verified that FOXQ1 can be used as an independent indicator of the prognosis of CRC patients.

The epidermal growth factor receptor (EGFR) signaling pathway plays an important role in physiological processes, such as cell growth, proliferation, and differentiation. The expression level of EGFR gradually increases from normal mucosa, adenomas with low-grade dysplasia, and adenomas with high-grade dysplasia to CRC, confirming that EGFR plays an important role in CRC[9]. Heparin-binding epidermal growth factor (HB-EGF) is one of the seven major ligands of EGFR. HB-EGF was originally identified as a secreted product from human macrophage U937 cells, and it induces cell proliferation and differentiation[10,11]. Soluble HB-EGF (sHB-EGF) is the main stimulator of cell proliferation. The affinity of sHB-EGF binding to target cells to promote proliferation of the target cells and activation of EGFR tyrosine kinase activity is 20-40 times higher than that of EGF. Thus, sHB-EGF is the most effective EGFR signaling pathway activator[12,13]. Membrane-bound HB-EGF (proHB-EGF) is affected by a variety of proteins, and it can be shed into active sHB-EGF to promote cell proliferation[14-16].

There is evidence that poor prognosis and low survival rates for CRC are associated with abnormally activated signaling pathways, including the EGFR signaling pathway[17]. In advanced CRC the most commonly used targeted therapies are the monoclonal antibodies cetuximab and panitumab, which block EGFR activation[18]. Activation of EGFR signaling leads to resistance to chemotherapy in CRC cells and promotes cell survival, while inhibition of EGFR signaling significantly reduces proliferation in CRC cells[19]. In nasopharyngeal carcinoma, Luo et al[20] reported that FOXQ1 induces vasculogenic mimicry through the EGFR signaling pathway, thereby promoting the metastasis of nasopharyngeal carcinoma cells. Our previous study suggested that FOXQ1 positively regulates the expression of HB-EGF and triggers the activation of the EGFR pathway in CRC[21]. However, the role and mechanism of how FOXQ1 promotes tumorigenesis in CRC by activating the HB-EGF/EGFR pathway remain largely unknown. Therefore, we analyzed the correlation between FOXQ1 and the HB-EGF/EGFR pathway by constructing FOXQ1 knockdown cells and using tissue microarrays, cell function experiments, quantitative real-time polymerase chain reaction (qRT-PCR), and western blotting. Our findings elucidated the critical role of the FOXQ1 and HB-EGF/EGFR pathways in CRC, providing theoretical support for the clinical application of targeted FOXQ1 in the treatment of CRC.

MATERIALS AND METHODS

Cell cultures

The human CRC cell lines, DLD1 and SW480, as well as the human embryonic kidney 293 (HEK293) cell line were purchased and authenticated from the Cell Bank of the Chinese Academy of Science in Shanghai, China. DLD1 cells were cultured in RPMI 1640. SW480 and HEK293 cells were cultured in DMEM. All media were supplemented with 10% heat-inactivated fetal bovine serum (Life Technologies BRL), and the cells were maintained in a 5% CO2-humidified atmosphere at 37 °C.

Plasmid construction and transfection

Three siRNAs targeting the human FOXQ1 sequence (NM_033260.3) were designed using siRNA Target Finder (InvivoGen, San Diego, CA, United States), and one scrambled siRNA was designed as a negative control, refer to our previous research[21]. The pSPAX2 packaging system from Addgene was used to construct a lentiviral PLKO.1 vector (PLKO.1-puro-shFOXQ1). After each group of recombinant plasmids was confirmed by sequencing, lentiviral vectors and packing vectors (pRSV-rev, pMDlg-pRRE and pCMV-VSV-G) were cotransfected into HEK293 cells using Lipofectamine® 2000 transfection reagent (Life Technologies). Lentivirus was collected to infect DLD1 and SW480 cells. Stable cells were generated after selection with puromycin (Solarbio, Beijing, China) (0.4 ng/μL for DLD1 cells and 0.1 ng/μL for SW480 cells) for 7-14 d after infection. The most effective knockdown cells were designated DLD1-shFOXQ1 and SW480-shFOXQ1, and the corresponding controls were named DLD1-shControl and SW480-shControl, respectively.

Flow cytometry

Cell surface HB-EGF was detected using APC-conjugated anti-human HB-EGF (eBioscience, San Diego, CA, United States). APC-conjugated mouse IgG2ak isotype was used as a control (eBioscience, San Diego, CA, United States) according to the manufacturer’s directions. Briefly, CRC cells were harvested and blocked with blocking buffer (PBS containing 2% BSA) for 10 min at 4 °C and then stained with 5 μL of monoclonal HB-EGF antibody (eBioscience, San Diego, CA, United States) for 30 min at 4 °C. After two washes, cells were resuspended in 100 μL of PBS. Samples were analyzed using a MoFlo flow cytometer (Beckman Coulter) and FlowJo software (Becton, Dickinson and Company).

Western blotting

Protein extracts were isolated from 106 treated cells using mammalian cell lysis reagent containing protease inhibitor cocktail (Sigma–Aldrich, St. Louis, MO, United States) and phosphatase inhibitors (Roche, United States) following the manufacturer’s directions. Equal amounts of protein (30 μg) were resolved on a 10% sodium dodecyl sulfate (SDS)-precast polyacrylamide gel (Bio–Rad Laboratories) and transferred to an Immobilon-polyvinylidene difluoride membrane (Millipore, Billerica, MA, United States). The membranes were blocked and incubated with the following primary antibodies: FOXQ1 (Abcam, Cambridge, MA, United States), HB-EGF (Abcam, Cambridge, MA, United States), phosphor-PI3K (Cell Signaling Technology, Cold Spring Harbor, NY, United States), PI3K (Proteintech, Wuhan, China), Akt (Proteintech, Wuhan, China), phosphor-Akt (Proteintech, Wuhan, China), phosphor-MAPK (Cell Signaling Technology, Cold Spring Harbor, NY, United States), MAPK (Proteintech, Wuhan, China), EGFR (Proteintech, Wuhan, China), KRAS (Proteintech, Wuhan, China), RAF (Proteintech, Wuhan, China), E-cadherin, N-cadherin, vimentin and β-actin (Proteintech, Wuhan, China). Blots were then incubated with the appropriate peroxidase-conjugated secondary antibody as follows: HRP-Rb-anti-goat (Cell Signaling Technology, Cold Spring Harbor, NY, United States) or HRP-goat-anti-mouse (Proteintech, Wuhan, China), respectively. The proteins were detected using an ECL system (Millipore, Braunschweig, Germany) and visualized with a ChemiDoc XRS system (Bio–Rad, Hercules, CA, United States).

RNA isolation and qRT–PCR

Total RNA was isolated using TRIzol (Invitrogen, Carlsbad, CA, United States) according to the manufacturer’s protocol. Total RNA (1 μg) was reverse transcribed into cDNA using PrimeScript RT reagent (TakaraBio, Japan), and qRT–PCR was performed using a LightCycler 480 (Roche, United States) with SYBR Premix Ex Taq (Takara, China). Each sample was analyzed in triplicate, and Glyceraldehyde-3-phosphate dehydrogenase(GAPDH) was used as the reference. Quantitative results were calculated using the 2-∆∆CT method. Primers used for qRT–PCR were designed and synthesized by Takara (Dalian, China) (Table 1).

Table 1

Oligonucleotide sequences used for real-time polymerase chain reaction

Gene name	Forward (5'→3')	Reverse (5'→3')
FOXQ1	TGACTTCAACAGCGACACCCA	CACCCTGTTGCTGTAGCCAAA
EGFR	AGACGCAGATAGTCGCCCAAAG	TCCATCAGGGCACGGTAGAAG
HB-EGF	CCATGTCTTCGGAAATACAAGGA	CCAGGATGGTTGTGTGGTCA
Akt	AGACGCAGATAGTCGCCCAAAG	TCCATCAGGGCACGGTAGAAG
RAF1	AGACGCAGATAGTCGCCCAAAG	TCCATCAGGGCACGGTAGAAG
KRAS	GTGGACGAATATGATCCAACAATAG	TGCTGTGTCGAGAATATCCAAGAG
GAPDH	GGCAACGGGCTACAGCTTTA	GGCACCCCACATACATAATCAA

Epidermal growth factor/platelet-derived growth factor pathway cDNA array assay

The human epidermal growth factor (EGF)/platelet-derived growth factor (PDGF) Signaling RT² Profiler™ PCR Array (SABiosciences), which profiles the expression of 84 genes related to EGF/PDGF-mediated signal transduction, five housekeeping genes and three controls, was used to analyze the effect of FOXQ1 on EGF/PDGF signaling-related gene expression (Table 2). Total RNA was extracted with TRIzol reagent according to the manufacturer’s manual. DNase treatment was performed by amplification grade I DNase I (Sigma–Aldrich, St. Louis, MO, United States) according to the manufacturer’s instructions. Each total RNA preparation (5 μg) was digested with 1 μL of DNase I (1 unit/μL) and 1 μL of 10 reaction buffer in a volume of 10 μL. After incubation and addition of Stop Solution, DNase I was denatured by incubation at 70 °C for 10 min. The RNA samples were kept on ice for another 5 min and then converted into cDNA with the RT2 PCR Array First Strand Kit (SuperArray) according to the manufacturer’s protocol. cDNA (20 ng) was combined with RT2 SYBR Green/Fluorescein PCR master mix (SuperArray), and equal amounts of this mixture (25 μL) were added to each well of the RT2 Profiler PCR plate containing the predispensed gene-specific primer sets. PCR cycles were performed according to the manufacturer's instructions. The relative level of mRNA expression for each gene in each sample was first normalized to the expression of GAPDH in that sample and then normalized to the level of mRNA expression in the DLD1-shControl.

Table 2

List of up- and down-regulated genes between DLDL-shFOXQ1 and DLD1-shControl with a fold-change ≥ 2.0 or ≤ -2.0 by using epidermal growth factor/platelet-derived growth factor SuperArray

Symbol	Gene information	Fold change (shFOXQ1/shControl)	P value	FOXQ1 binding sites (TSS: -2k_+1k)
PLAT	Plasminogen activator, tissue	-8.01	0.0375	3
COL1A1	Collagen, type I, alpha 1	-5.72	0.0138	1
ELK1	ELK1, member of ETS oncogene family	-5.45	0.0228	1
BCAR1	Breast cancer anti-estrogen resistance 1	-3.97	0.0479	2
PIK3R2	Phosphoinositide-3-kinase, regulatory subunit 2	-3.63	0.0197	3
NCK2	NCK adaptor protein 2	-3.51	0.0237	6
JUN	Jun proto-oncogene v-jun	-3.43	0.0077	12
EGR1	Early growth response 1	-3.36	0.0166	11
LTA	Lymphotoxin alpha (TNF superfamily, member 1)	-3.18	0.0796	3
PDGFB	Platelet-derived growth factor beta	-2.91	0.1896	5
STAT3	Signal transducer and activator of transcription 3 (acute-phase response factor)	-2.49	0.0038	6
HB-EGF	Heparin-binding EGF-like growth factor	-2.19	0.0278	3
FN1	Fibronectin 1	-2.10	0.0005	6
FASLG	Fas ligand (TNF superfamily, member 6)	-2.08	0.2026	17
HRAS	V-Ha-ras Harvey rat sarcoma viral oncogene homolog	-2.06	0.0058	2
MAP2K1	Mitogen-activated protein kinase kinase 1	-2.02	0.0243	5
AKT3	V-akt murine thymoma viral oncogene homolog 3	-2.02	0.1440	11
FOS	FBJ murine osteosarcoma viral oncogene homolog	-2.01	0.0470	5

ELISA

Exogenous recombinant HB-EGF protein at a final concentration of 50 ng/mL was added to the cell culture medium when the cell density reached 80%, and the culture medium was changed after incubation for 24 h[22]. The supernatant was collected after the cells were cultured for another 24 h. The protein concentrations of soluble HB-EGF, ADAM9, ADAM10, ADAM12 and MMP-7 in the cell culture medium were determined by ELISA detection kits against human HB-EGF, ADAM9, ADAM10, ADAM12 and MMP-7, respectively (R&D Systems, Minneapolis, United States).

Cell proliferation

DLD1-shFOXQ1 and SW480-shFOXQ1 cells (2,000 cells/well in a 96-well plate) were incubated with medium containing 10% FBS at 37 °C for 24, 48, 72 and 96 h. At the end of incubation, 10 μL of Cell Counting Kit-8 (CCK-8) solution (Beyotime Biotech, Shanghai, China) was added to each well with 100 μL of medium and incubated for another 4 h at 37 °C, and the OD at 450 nm was measured by a microplate reader (BioTek, Winooski, VT, United States). The effect of siRNA FOXQ1 on CRC cell viability was assessed as the percent of cell viability compared to vehicle-treated control cells, which were arbitrarily assigned as 100% viability.

Cell migration and wound-healing assay

Cell migration was analyzed using Transwell inserts with 8.0 μm membrane pores (BD, San Jose, CA, United States) according to the manufacturer’s protocol. Migration was additionally evaluated with the wound-healing assay. Briefly, DLD1-shFOXQ1/SW480-shFOXQ1 and DLD1-shControl/SW480-shControl cells were seeded in 6-well plates at a density that enabled a confluency of 80% to be attained 24 h after plating. A 10 μL filter tip was used to gently scratch the cell monolayer across the center of the well. Cells were then gently washed twice with PBS to remove the dislodged cells, and fresh medium was added. The first images of the scratch area were then acquired. Cells were cultured in serum-free medium for another 48 h, and a second set of images was then acquired to determine the extent of wound closure.

CRC tissue microarray

Tissue microarrays containing a total of 90 pairs of colorectal tumor tissues and matched adjacent normal tissues, together with pathological staging data in accordance with TNM classification of the American Joint Committee on Cancer (2010) and follow-up survival time after surgery, were obtained from Shanghai Biochip Co. Ltd., Shanghai, China (HCol-Ade180Sur-06). FOXQ1 (ab51340) and HB-EGF (ab192545) antibodies were purchased from Abcam (Cambridge, MA, United States). Tissue microarray analysis was performed using a standard immunohistochemistry (IHC) protocol. The median value of the immunoreactivity score (IRS) was selected as the cutoff for high and low protein expression levels based on a measure of heterogeneity according to the log-rank test with respect to disease-specific survival (DSS) as described previously. Cutoff values for the scoring system were assigned as follows: High expression of FOXQ1 and HB-EGF were defined as an IRS of ≥ 4 (4, 6, 8, 9 and 12); and low expression was defined as an IRS of < 4 (0, 1, 2 and 3)[23]. Immunostained sections were scanned using a microscope (Axiovert 200, Carl Zeiss, Gttingen, Germany). Data for 25 patients were excluded because the dots were off the chips during the experiment. In total, data for 65 patients with CRC were included in the final analysis (Tables 3 and 4).

Table 3

Correlation between heparin binding epidermal growth factor expression and clinicopathological characteristics of colorectal cancers in cohort of human colorectal cancer tissues

Clinicopathological variables		Cohort tumor HB-EGF expression		P value
		Negative (n = 26)	Positive (n = 39)
Age (yr)		65.72 (10.37)	69.43 (7.98)
Gender	Female	11	23	> 0.05
	Male	14	16
Size of tumor (mm)	≤ 10	5	7	> 0.05
	> 10	21	32
Lymphatic metastasis	Absent	19	19	> 0.05
	Present	7	20
Tumor differentiation	I-II	12	24	> 0.05
	III-IV	14	15
TNM stage	I-II	11	8	> 0.05
	III-IV	15	31
AJCC clinical stage according to 7th issue	1 and 2A	17	24	> 0.05
	3 and 3B	9	15

HB-EGF: Heparin binding epidermal growth factor; AJCC: American Journal of Critical Care.

Table 4

Correlation analysis of forkhead Box q1 and heparin binding epidermal growth factor expression in cohort (n = 65) colorectal cancer tissues

		FOXQ1, negative (n = 31)	FOXQ1, positive (n = 34)	χ2	P value
HB-EGF	Negative (n = 26)	20	6	4.116	< 0.05
	Positive (n = 39)	11	28

FOXQ1: Forkhead Box q1; HB-EGF: Heparin binding epidermal growth factor.

Statistical analysis

Statistical analyses were performed using GraphPad Prism 8 and SPSS v.19. An unpaired two-tailed Student’s t-test was performed for two-group comparisons, and one-way analysis of variance (ANOVA) was performed for multiple group comparisons. Survival curves were calculated using the Kaplan–Meier algorithm and log-rank test. P < 0.05 was considered to indicate a statistically significant difference.

RESULTS

Expression and prognosis of FOXQ1 and HB-EGF in CRC and normal colorectal tissues

To determine the expression levels of FOXQ1 and HB-EGF in CRC, we investigated the expression and prognosis of FOXQ1 and HB-EGF in CRC in the Gene Expression Profiling Interactive Analysis (GEPIA) online database (Figure 1). The results showed that FOXQ1 was upregulated in CRC compared to normal samples according to GEPIA (Figure 1A). Increased expression of FOXQ1 was also associated with worse overall survival (Figure 1C). The expression of HB-EGF in CRC was not significantly different from that in normal colorectal tissues (Figure 1B), and its expression had no significant effect on the overall survival of patients with CRC (Figure 1D).

Figure 1

Expression and prognosis of forkhead Box q1 and heparin binding epidermal growth factor in colorectal cancer and normal colorectal tissues. A and B: Forkhead Box q1 (FOXQ1) (left) and heparin binding epidermal growth factor (HB-EGF) (right) expression was analyzed in human colorectal cancer (CRC) tissues and normal tissues using the Gene Expression Profiling Interactive Analysis online database. aP < 0.05; C and D: Overall survival rate of FOXQ1 (left) and HB-EGF (right) in CRC patients. FOXQ1: Forkhead Box q1; HB-EGF: Heparin binding epidermal growth factor; GEPIA: Gene Expression Profiling Interactive Analysis.

Construction of FOXQ1 knockdown CRC cell lines

To elucidate the functional roles of FOXQ1 in CRC, we generated CRC lines with stable FOXQ1 knockdown. Among the CRC lines we tested, DLD1 and SW480 cells had relatively high endogenous FOXQ1 expression as described previously[21]. Because the high expression of FOXQ1 in these two cells lines has also been confirmed in several other independent studies[3,4,8,9,24], we selected these two cell lines for knockdown studies. qRT–PCR assays verified significant FOXQ1 knockdown in DLD1-shFOXQ1 and SW480-shFOXQ1 cells. FOXQ1 mRNA expression in DLD1-shFOXQ1 and SW480-shFOXQ1 cells was significantly reduced compared to that in the control DLD1-shControl and SW480-shControl cells (Figure 2A). Western blot analysis also verified that the expression of FOXQ1 was significantly downregulated (Figure 2B).

Figure 2

Forkhead Box q1 knockdown induces heparin binding epidermal growth factor suppression and epidermal growth factor/platelet-derived growth factor signaling pathway blockade A: Quantitative real-time polymerase chain reaction was performed to evaluate mRNA expression in the knockdown of forkhead Box q1 in DLD1 cells (DLD1-shFOXQ1) and SW480-shFOXQ1 cells (two-tailed, unpaired Student’s t test). Bars represent the mean ± SE of three independent experiments; B: Western blotting was performed to evaluate FOXQ1 protein expression in DLD1-shFOXQ1 and SW480-shFOXQ1 cells; C: Volcano plot illustrating differentially expressed genes using the epidermal growth factor (EGF)/platelet-derived growth factor (PDGF) array to compare DLD1-shFOXQ1 and DLD1-shControl cells. The vertical blue points in the plot represent 84 genes related to the EGF/PDGF signaling pathways. The vertical green line indicates a fold-change in gene expression of 1. The vertical pink lines indicate the threshold for fold-change (≥ 2.0 and ≤ -2.0). The horizontal blue line indicates the threshold for the P value of the t test. aP < 0.05; bP < 0.01; cP < 0.001. DLD1-shFOXQ1: The knockdown of forkhead Box q1 in DLD1 cells; SW480-shFOXQ1: The knockdown of forkhead Box q1 in SW480 cells.

Knockdown of FOXQ1 suppresses the expression of HB-EGF and blocks the EGF/PDGF signaling pathway in vitro

A panel of PCR arrays consisting of 84 representative genes related to the EGF/PDGF signaling pathways was used to detect the transcriptional signatures of DLD1-shFOXQ1 and DLD1-shControl cells. Differentially expressed genes with statistical significance were identified by volcano plot filtering. The results showed that 18 genes had expression changes with a fold-change ≥ 2.0 (P < 0.05)(Figure 2C). Among these 18 genes, 12 genes were associated with EGFR signaling pathways, including HB-EGF (FC = -2.19), and there were several novel genes of the three main EGF/PDGF downstream signaling pathways (MAPK/ERK1/2, PI3K/AKT and JAK/STAT3). These novel genes included RAS, PI3K and STAT3 as well as several intracellular transcription factors activated by these signaling pathways, such as c-JUN and c-FOS (Table 2).

Effect of FOXQ1 expression on HB-EGF expression and extracellular release ability in CRC cells

Western blot analysis confirmed that knockdown of FOXQ1 in CRC cells resulted in a significant decrease in the expression of HB-EGF in DLD1 cells (Figure 3A), but no significant changes in HB-EGF were observed in SW480 cells (Figure 3A). Flow cytometry analysis confirmed that knockdown of FOXQ1 expression in CRC cells did not affect the expression level of proHB-EGF (Figure 3B). ELISA results indicated that knockdown of FOXQ1 significantly reduced the shedding of soluble HB-EGF (sHB-EGF) in DLD1 cells but did not affect the shedding of HB-EGF in SW480 cells (Figure 3C). Among the four proteins that affect the release of the extracellular domain of proHB-EGF, three (ADAM9, ADAM12 and ADAM7) were decreased significantly in both DLD1-shFOXQ1 and SW480-shFOXQ1 cells, while ADAM10 was not changed. The feedback regulation of ADAM9, ADAM12 and MMP7 secretion by HB-EGF from CRC cells was studied by adding exogenous recombinant human HB-EGF (rhHB-EGF). The results showed that rhHB-EGF reversed the decline in the expression of ADAM9, ADAM12 and MMP7 in DLD-shFOXQ1 cells. However, only the decreased expression of ADAM9 was reversed by rhHB-EGF in SW480-shFOXQ1 cells (Figure 3D).

Figure 3

Effect of forkhead Box q1 expression on heparin binding epidermal growth factor expression and extracellular release ability in colorectal cancer cells. A: Quantitative analysis of western blots was performed to detect the expression of heparin binding epidermal growth factor (HB-EGF) with forkhead Box q1 (FOXQ1) inhibition; B: Flow cytometry was performed to detect the effect of changes in FOXQ1 expression on the expression of proHB-EGF in DLD1 and SW480 cells; C: ELISA verified the effect of FOXQ1 knockdown on sHB-EGF expression in DLD1 and SW480 cells; D: ELISA confirmed the effect of FOXQ1 on the expression of proHB-EGF extracellular release proteins (ADAM9, ADAM10, ADAM12 and MMP-7) in DLD1 and SW480 cells. aP < 0.05; bP < 0.01; cP < 0.001. FOXQ1: Forkhead Box q1; HB-EGF: Heparin binding epidermal growth factor; ProHB-EGF: Membrane-bound HB-EGF; SHB-EGF: Souble HB-EGF; ADAM9: A disintegrin and a metalloprotease 9; ADAM10: A disintegrin and a metalloprotease 10; ADAM12: A disintegrin and a metalloprotease 12; MMP7: Matrix metallopeptidase 7.

FOXQ1 regulates EGFR and downstream signaling pathways by regulating HB-EGF

To verify the role of FOXQ1 in activating the HB-EGF/EGFR signaling pathway, we performed qRT–PCR analysis in DLD1-shFOXQ1 and SW480-shFOXQ1 cells. The results confirmed that FOXQ1 knockdown resulted in a significant decrease in the mRNA expression of HB-EGF, EGFR and downstream genes (AKT, RAF and KRAS) in DLD1-shFOXQ1 and SW480-shFOXQ1 cells (Figure 4A and B). In addition, the Western blot analysis showed that knockdown of FOXQ1 resulted in a significant decrease in HB-EGF and EGFR expression as well as decreased AKT and MAPK phosphorylation in DLD1-shFOXQ1 cells compared to DLD1-shControl cells. In SW480-shFOXQ1 cells, EGFR and HB-EGF expression was significantly decreased, and AKT and PI3K phosphorylation was inhibited compared to that in SW480-shControl cells (Figure 4C). Furthermore, Western blot analysis confirmed that the decreased expression of important downstream genes was rescued by rhHB-EGF protein in either DLD1-shFOXQ or SW480-shFOXQ1 cells (Figure 4D).

Figure 4

Forkhead Box q1 regulates epidermal growth factor receptor and downstream signaling pathways by regulating heparin binding epidermal growth factor. A: Quantitative real-time polymerase chain reaction (qRT–PCR) was performed to analyze the levels of forkhead Box q1 (FOXQ1), heparin binding epidermal growth factor (HB-EGF), epidermal growth factor receptor (EGFR), Akt, kRas and Raf in DLD1-shFOXQ1 and DLD1-shControl cells; B: qRT–PCR was performed to analyze the levels of FOXQ1, HB-EGF, EGFR, AKT, RAF and KRAS in SW480-shFOXQ1 and SW480-shControl cells; C: Western blot analyses were performed to analyze the protein levels of FOXQ1, HB-EGF, EGFR, RAF and KRAS as well as the phosphorylation levels of PI3K, Akt, and MAPK in DLD1-shFOXQ1, SW480-shFOXQ1, DLD1-shControl and SW480-shControl cells. β-actin was used as the loading control; D: Exogenous rhHB-EGF (50 ng/mL) was added to the culture media of DLD1-shFOXQ1 and SW480-shFOXQ1 cells for 1.5 h and 2 h, respectively. Western blotting was performed to detect the expression of proteins in the HB-EGF/EGFR signaling pathway. aP < 0.05; bP < 0.01; cP < 0.001 indicate significant differences compared to the control group (two-tailed, unpaired Student’s t test). Bars represent the mean ± S.E. of three independent experiments. rhHB-EGF: Recombinant human HB-EGF; FOXQ1: Forkhead Box q1; HB-EGF: Heparin binding epidermal growth factor.

rhHB-EGF reverses the FOXQ1 knockdown-induced suppression of CRC cell proliferation and migration in vitro

CCK-8 results confirmed that FOXQ1 knockdown significantly inhibited the proliferation of DLD1-shFOXQ1 and SW480-shFOXQ1 cells, and the inhibitory effect was partially reversed by exogenous rhHB-EGF (Figure 5A). The results of the Transwell migration assay confirmed that FOXQ1 knockdown also reduced the migration of DLD1 and SW480 cells, and the inhibitory effect was also reversed to a large extent by rhHB-EGF (Figure 5B and C). Scratch experiment results confirmed that FOXQ1 knockdown reduced the wound-healing ability of DLD1 and SW480 cells, which was also reversed to a large extent by rhHB-EGF (Figure 5D and E). We next analyzed the protein expression during cell invasion and metastasis (Ecadherin, N-cadherin, Vimentin and Snail), Western blot analysis indicated that FOXQ1 knockdown reduced the expression in DLD1-shFOXQ1 and SW480-shFOXQ1 cells (Figure 5F).

Figure 5

Recombinant human heparin binding epidermal growth factor reverses the forkhead Box q1-induced suppression of colorectal cancer cell proliferation and migration A: A cell counting Kit-8 (CCK-8) assay was performed to detect the effects of recombinant human heparin binding epidermal growth factor (rhHB-EGF) on the proliferation of forkhead Box q1 (FOXQ1)-deficient DLD1 and SW480 cells; B: The Transwell migration assay was performed to evaluate the migration ability after knocking out FOXQ1 and after adding rhHB-EGF; C: A statistical analysis of the results of the Transwell migration assay was performed. bP < 0.01; D: After 0 h and 24 h, the scratch experiment evaluated the migration ability of DLD1 and SW480 cells with downregulated FOXQ1 and rhHB-EGF; E: Statistical analysis of the above scratch experiment results. aP < 0.05; bP < 0.01; cP < 0.001; F: Quantitative analysis of Western blot to detect the expression of E-cadherin, N-cadherin, vimentin and Snail protein levels in DLD-shFOXQ1 and SW480-shFOXQ1 cells compared to DLD1-shControl and SW480-shControl cells. β-actin was used as the loading control. CCK-8: Cell Counting Kit-8.

Prognostic value of the combination of FOXQ1 and HB-EGF

To verify the clinical relevance of our findings, we evaluated the expression of FOXQ1 and HB-EGF in human CRC tissue biopsies (cohort, n = 65) (Tables 3 and 4). IHC analysis showed that FOXQ1 was significantly upregulated in CRC tissues compared to adjacent nontumorous tissues and that HB-EGF was moderately upregulated (Figure 6A). Although the overexpression of HB-EGF had no significant correlation with any clinicopathological characteristics of CRCs (Table 3), further analysis verified the positive correlation between FOXQ1 and HB-EGF (Table 4). In a cohort of 65 CRC patients, Kaplan–Meier survival analysis results showed that CRC patients with positive expression of FOXQ1 had shorter overall survival than those with negative expression of FOXQ1. Furthermore, Kaplan–Meier survival analysis of CRC patients with positive coexpression of FOXQ1 and HB-EGF had the shortest overall survival times compared to the corresponding single-negative or double-negative groups in the cohort of 65 CRC patients (Figure 6B).

Figure 6

Prognostic value of the combination of forkhead Box q1 and heparin binding epidermal growth factor. A: Representative immunohistochemical images of forkhead Box q1 (FOXQ1) and heparin binding epidermal growth factor (HB-EGF) positive coexpression and negative coexpression in human colorectal cancer tissues and adjacent nontumorous tissues. Scale bars represent 200 μm; B: Kaplan-Meier survival analysis was used to evaluate the overall survival rate in the following patients: FOXQ1-positive and FOXQ1-negative patients; HB-EGF-positive and HB-EGF-negative patients; and FOXQ1-positive and FOXQ1-negative coexpressing HB-EGF-positive and HBEGF-negative patients. aP < 0.05; bP < 0.01; cP < 0.001 indicate a significant difference between the indicated groups (log-rank test). FOXQ1: Forkhead Box q1; HB-EGF: Heparin binding epidermal growth factor.

DISCUSSION

Studies have indicated that FOXQ1 is an oncogene in multiple tumors, including CRC[24], breast cancer[25], lung cancer[26], gastric cancer[27], liver cancer[28], pancreatic cancer[29], ovarian cancer[30], and neuroglioma[31]. In CRC, FOXQ1 promotes tumor invasion and metastasis through the Wnt signaling pathway, and it affects the prognosis of patients[24]. EGFR is overexpressed in a variety of cancers, including CRC. EGFR overexpression and activation have a positive effect on the cell growth and metastasis of a variety of solid tumors, including CRC[17]. Metastasis is a multistep process in which tumor cells spread from the primary site to distant sites and form secondary tumors[32], and is the leading cause of death in CRC patients. Many downstream targets of HB-EGF and EGFR include MAPK, p-MAPK, RAS, RAF, PI3K, p-PI3K, AKT, p-AKT, and AKT protein kinases, leading to many processes related to tumor progression, including cell growth[33], epithelial-mesenchymal transition (EMT)[34], metastasis[35], and angiogenesis[36]. FOXQ1 can regulate FAK, PI3K, AKT, and many other key proteins in the PI3K/AKT signaling pathway and promotes the phosphorylation of the above proteins to maintain the activation of PI3K/AKT signaling[37]. FOXQ1 can also combine with VEGFR2 and VE-cadherin to promote angiogenesis and endothelial cell migration and rearrangement[24].

In our previous studies, we found that the mRNA and protein expression levels of FOXQ1 gradually increase with the pathological development of colorectal adenoma to colorectal adenocarcinoma, and the increased expression of FOXQ is not only involved in the process of colorectal adenoma carcinogenesis, but is also closely related to invasion and metastasis of CRC[38]. The differential expression of genes that are involved in the process of colorectal adenoma carcinogenesis screened by the gene chip has been analyzed by signal pathway analysis, suggesting that the abnormally high expression of FOXQ1 is closely related to the activation of the EGFR signaling pathway[39]. Studies have shown that the abnormal activation of the EGFR pathway plays an important role in malignant growth, invasion, and metastasis of colorectal tumors. When sHB-EGF binds to EGFR, the tyrosine kinase activity of EGFR is activated, mainly by activating the MAPK/ERK1/2, PI3K/AKT, and JAK/STAT3 signaling pathways, leading to the proliferation, invasion, metastasis, and apoptosis of tumor cells[40]. The results of this study also proved that knockout of FOXQ1 caused a decrease in gene expression in these three signaling pathways.

Combined with the analysis in Figure 3B-D, these results suggested that FOXQ1 may regulate the EGFR pathway by promoting the separation of ProHB-EGF into sHB-EGF. AMAD7, AMAD9, and AMAD12 are also factors that affect the separation of ProHB-EGF[41]. In this study, decreased expression levels of AMAD7, AMAD9, and AMAD12 were also observed in DLD1-shFOXQ1 and SW480-shFOXQ1 cells, but this indirect regulation should be further verified in other CRC cell lines.

In this study, FOXQ1 was determined to be upregulated in CRC cell lines. The results showed that FOXQ1 knockdown inhibited the proliferation, migration and repair capabilities of CRC, which was consistent with our previous results[21,42]. When rhHB-EGF protein was added, the proliferation ability of the cells was not completely restored, but the migration and repair ability of CRC cells was partially restored. These results indicated that the effect of FOXQ1 in promoting the proliferation of CRC cells is not directly mediated by HB-EGF; however, the regulation of the invasion and metastasis of CRC cells by FOXQ1 is partly mediated by HB-EGF. HB-EGF is related to the abnormal proliferation of skin and mucosal cells, and high expression of HB-EGF is closely related to the occurrence and development of a variety of tumors; the expression of the HB-EGF gene is significantly increased in various human cancers and cancer-derived cell lines, indicating that HB-EGF plays an important role in tumor invasion and metastasis[43-46].

Studies have shown that FOXQ1 is related to the poor prognosis of CRC[47]. In this study, we conducted pathological and survival analyses on 65 CRC patients. These findings suggested that the expression of FOXQ1 and its coregulatory protein HB-EGF may have a prognostic correlation with CRC. Thus, FOXQ1 may serve as a therapeutic target for CRC treatment by blocking the HB-EGF/EGFR pathway. Our research suggests that FOXQ1 activates the EGFR signaling pathway by regulating the expression of HB-EGF, thereby affecting the invasion and metastasis of CRC. Regarding the limitations of this study, it is necessary to construct a FOXQ1 high-expressing cell line and combine it with the dual luciferase reporter gene system to further verify whether FOXQ1 is directly involved in the transcriptional regulation of HB-EGF. Next, we will further explore the regulation of FOXQ1 on EGFR and its downstream signaling pathways through in vivo and in vitro studies. We will conduct a more comprehensive study on the role of HB-EGF in the invasion and metastasis of CRC to provide more possibilities for the treatment of CRC.

CONCLUSION

In conclusion, we have demonstrated that decreased FOXQ1 expression was also associated with a lower ability to invade and metastasize in CRC. FOXQ1 promotes the invasion and metastasis of CRC by activating the HB-EGF/EGFR pathway. These data indicated that FOXQ1 and HB-EGF may be potential biomarkers to improve the accuracy of CRC diagnosis and treatment.

ARTICLE HIGHLIGHTS

Research background

Invasion and metastasis play important roles in tumorigenesis, resulting in the death of most colorectal cancer (CRC) patients. Forkhead Box q1 (FOXQ1) is a well-established oncogene in multiple tumors, including CRC. However, the role and mechanism of how FOXQ1 promotes tumorigenesis in CRC by activating the heparin-binding epidermal growth factor (HB-EGF)/epidermal growth factor receptor (EGFR) pathway remain largely unknown.

Research motivation

Our study aims to elucidate the critical role of FOXQ1 and HB-EGF/EGFR pathways in CRC, and to provide theoretical support for the clinical application of targeted FOXQ1 in the treatment of CRC.

Research objectives

To determine the role of FOXQ1-induced invasion and metastasis, which are related to activating the HB-EGF/EGFR pathway, and to explore the mechanism by which FOXQ1 promotes tumorigenesis by activating the HB-EGF/EGFR pathway in CRC.

Research methods

We analyzed the correlation between FOXQ1 and the HB-EGF/EGFR pathway by constructing FOXQ1 knockdown cells, tissue microarray, cell function experiments, quantitative real-time polymerase chain reaction (qRT–PCR), flow cytometry, ELISA, western blot, and the Gene Expression Profiling Interactive Analysis (GEPIA) website.

Research results

GEPIA showed that the expression of FOXQ1 in CRC tissues was relatively high and was related to a lower overall survival rate. PCR array results showed that FOXQ1 is related to the HB-EGF and EGFR pathways. Knockdown of FOXQ1 suppressed the expression of HB-EGF and led to a decrease in EGFR and its downstream genes AKT, RAF, KRAS expression levels. After knockdown of FOXQ1 in CRC cell lines, cell proliferation, migration, and invasion were attenuated. Adding HB-EGF restored the migration and invasion ability of CRC, but the cell proliferation ability was not restored. Kaplan–Meier survival analysis results showed that the combination of FOXQ1 and HB-EGF may serve to predict CRC survival.

Research conclusions

FOXQ1 promotes the invasion and metastasis of CRC by activating the HB-EGF/EGFR pathway.

Research perspectives

In this study, our results indicated that FOXQ1 and HB-EGF may be potential biomarkers to improve the accuracy of CRC diagnosis and treatment.