Hepatocyte Nuclear Factor 3β Plays a Suppressive Role in Colorectal Cancer Progression

Background and Objective: Hepatocyte nuclear factor 3β (HNF3β) is a key transcription factor in the development of the gastrointestinal tract. However, only few studies have examined its' expression, function and potential clinical significance in colorectal cancer tumorigenesis and progression.

Methods: HNF3β expression in colorectal cancer tissue samples of 174 patients was assessed by immunohistochemistry. The results were analyzed with respect to patients' clinicopathological characteristics and survival. Following the in vitro cell transfection, MTT, wound healing, and Transwell assays were used to test cell proliferation, migration, and invasion, respectively. Western blot was used to examine IL6, JAK1, and STAT3 protein expression. The potential for tumor formation was evaluated using a mouse xenograft model.

Results: HNF3β expression was lower in colon cancer tissue compared to normal tissue and correlated with UICC clinical stage (P = 0.001), depth of invasion (P = 0.004), regional lymph node metastasis (P = 0.007), distant metastasis (P = 0.048), and poor survival (P < 0.001) in patients with colorectal cancer. Furthermore, HNF3β overexpression impeded proliferation, migration and invasion of SW480 cells via JAK-STAT3 signaling in vitro. Moreso, HNF3β overexpression showed a significant growth inhibition of subcutaneous xenograft tumors in vivo.

Conclusions: The results show that HNF3β acts as a suppressor of colorectal cancer progression and decreased HNF3 β expression is closely related to the poor prognosis. Thus, HNF3β may be a potential molecular target for inhibition of colorectal cancer cells and development of new anti-tumor therapies.

Introduction

Although some progress has been made in the management of colorectal cancer based on early diagnosis, radical surgery, radiotherapy, and neoadjuvant chemotherapy, recurrence and metastases remain a major cancer management issue which has yet to be solved (1, 2). Therefore, there is an urgent need to explore and define potential molecular biomarkers for early diagnosis, as well as tumor recurrence and metastasis. However, regulatory mechanisms underlying colorectal cancer occurrence and recurrence or metastasis remain to be determined, as few prognostic biomarkers have been identified.

Hepatocyte nuclear factor 3β (HNF3β), is a key member of the hepatocyte nuclear factor (HNF) family with an important role in embryogenesis (3). Indeed, it has been shown that HNF3 transcription factors, especially HNF3α/FOXA1, HNF3β/FOXA2, and HNF3γ/FOXA3, are involved in the development and differentiation of the gastrointestinal tract (4). HNF3β mRNA is expressed during embryogenesis, and is detected during gastrulation in the anterior primitive streak and node with subsequent expression in the gut (5–7). It has been shown that HNF3β deficiency is lethal in the embryonic stage due to failed node and notochord development (8). Moreover, HNF3β is important for postnatal development due to its role in the regulation of normal bile acid and lipid homeostasis (9). Many studies have shown that HNF3β suppression is closely related to tumorigenesis, progression, and metastasis of several different types of cancer including hepatocellular carcinoma, lung cancer, pancreatic cancer, gastric cancer, and thyroid cancer (10–14). However, the role of HNF3β in colorectal cancer development and progression remains to be determined and thus we have decided to examine its role in colorectal cancer patients as well as in the colorectal cancer model in vivo and in vitro.

Materials and Methods

Patients and Tissue Samples

The use of human tissues was approved by Research Ethics Committee of the Fourth Military Medical University and each patient signed the informed consent for using their tissues in the study. All procedures involving human participants were performed in accordance with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards.

Eight inflammation tissues, 14 low-grade and 18 high-grade intraepithelial neoplasia tissues, 174 primary colon cancer tissues, as well as six normal colon tissues were obtained from surgical resections at the Department of Gastrointestinal Surgery of Xijing Hospital, Fourth Military Medical University from October 2010 to November 2011. None of these patients received neoadjuvant chemotherapy or radiotherapy before surgery and all received standard chemotherapy after surgery.

Tumor tissues were histologically examined by two independent experienced pathologists at the Department of Pathology, Xijing Hospital, Fourth Military Medical University. Tumor staging was performed according to the Union for International Cancer Control (UICC) staging system. The following clinicopathological parameters of patients and their tumors were obtained, including: age, gender, histological type and differentiation status, T status, lymph node and distant metastasis, and UICC stage (Table 1). Complete follow-up of patients for least 5 years for survival analysis was also available.

TABLE 1

Table 1. HNF3β protein expression and clinicopathological characteristics of 174 CRC patients.

Progression-free survival (PFS) was defined as the time from primary surgical resection to first local recurrence or distant metastasis, or the date of the last follow-up. Overall survival (OS) was considered the time from primary surgical resection to death or to the last date of follow-up.

Immunohistochemical (IHC) Analysis

Formalin-fixed paraffin-embedded tissues for the immunohistochemical analysis of protein expression were cut into 4-μm thick sections, which were subsequently deparaffinized in xylene and dehydrated in descending grades of ethanol. Sections were heated in 0.01 M sodium citrate (pH 6.0) for 15 min in the microwave for antigen retrieval. Endogenous peroxidase activity was blocked using methanol with 0.3% H₂O₂ for 10 min. The sections were blocked with 10% normal goat serum in phosphate-buffered saline (PBS) for 1 h at room temperature. The sections were incubated with rabbit monoclonal HNF3β antibody (1:250, Abcam 108422) in PBS at 4°C overnight in a humidified chamber. The next day sections were incubated with rabbit horseradish peroxidase-conjugated secondary antibody (1:500, Nichirei Japan) for 1 h at room temperature. The immunohistochemical reaction was visualized by incubating the sections with 0.1 mg/mL 3,3′-diaminobenzidine tetrahydrochloride (DAB) for 5 min at room temperature. All sections were evaluated by two independent trained pathologist.

Evaluation of IHC Staining

Results of the IHC staining were scored according to both the percentage and intensity of stained cells. Positive staining was scored 1, 2, 3, or 4 when the percentages of stained cells were 0, 1–25, 26–50, 51–75, and >75%, respectively. Staining intensity was scored as 0, 1, 2, or 3 when colorless, pallide-flavens, yellow, and brown. The final staining score was calculated by multiplying the scores for percentage and intensity. For product scores of nil, 1–4, 5–8, and 9–12, the final staining scores were: −, absent; +, weak; ++, moderate; and +++, strong. Tumors with absent or weak immunostaining were classified as having low target protein expression; tumors with moderate or strong immunostaining were classified as high target protein expression.

Cell Lines and Cell Culture

Cells of the human colon cancer cell line SW480 were purchased from Shanghai Cell Bank and cultured in McCoy's 5A growth medium (Invitrogen, Carlsbad, USA) supplemented with 10% heat inactivated fetal bovine serum (FBS, Gibco, USA) at 37°C in a humidified atmosphere with 5% CO₂. Cells of the human 293T cell line were purchased from Shanghai Cell Bank and cultured in Dulbecco's Modified Eagle's medium (DMEM; Invitrogen, Carlsbad, USA) at 37°C in a humidified atmosphere with 5% CO₂.

HNF3β Transfection

Full-length human HNF3β clone containing full-length human HNF3β cDNA was purchased from Abcam (UK). HNF3β cDNA was subcloned into the Sgf I and MluI digested fragment of pLVTHM/LENTI-BISVECTOR vector to obtain pLVTHM/LENTI-HNF3β. Recombinant lentivirus containing HNF3β was assembled in 293T cells by cotransfection of pLVTHM/LENTI-HNF3β, with pRsv-REV, pMDlg-pRRE, and pMD2G (Gemma, Germany). Lentiviral particles were harvested 48 h after transfection and purified by ultracentrifugation as lenti-HNF3β. As a control, lentiviruses without the HNF3β gene were generated and designated as lenti-green fluorescent protein (GFP). SW480 cells were transfected with lenti-HNF3β or lenti-GFP using Lipofectamine 2000 (Invitrogen, USA) in accordance with the manufacturer's instructions. Each transfection was performed in triplicates.

Cell Migration and Invasion Assays

To assess the effect of HNF3β on cell migration and invasion ability a transwell assay was performed using uncoated transwell chambers (Millipore, USA) or chambers coated with Matrigel (BD Bioscience, USA). In brief, the lower transwell chamber was filled with media containing 5% FBS, whereas the upper chamber was filled with 1 × 10⁶ cells suspended with medium containing 1% FBS. After incubation of the cells at 37°C for 24 h, the non-traversed cells on the upper surface of the filter (8 μm) were carefully wiped off with a cotton swab. Cells that migrated through the membrane on the bottom surface of the filter were fixed with methanol and stained with 0.1% crystal violet. Cell migration was assessed by counting cells that had migrated through the membrane at six randomly selected fields under the light microscope (200×, Olympus, Japan). For cell invasion assays, the same protocol was used except that the polycarbonate filter was pre-coated with Matrigel. Each experiment was performed in triplicate.

Wound Healing Assay

In addition to the transwell migration assay, a wound healing assay was performed to assess the effect of HNF3β on cell motility. Briefly, cells were seeded in 6-well plates and incubated at 37°C for 24 h to form a subconfluent cell monolayer (confluence of 90%). A straight linear wound (1 mm wide) was generated across the middle diameter of each well by gently dragging a sterile pipette tip through the monolayer. Cells were incubated at 37°C for 24 h to migrate across the open wound area. At the end of the incubation, the medium was removed and cells were washed and fixed in methanol. The wound healing process was photographed under a light microscope (Olympus, Japan) at 0, 12, and 24 h. Each experiment was performed in triplicate.

MTT Assay

The effect of HNF3β on cell proliferation was evaluated using the MTT assay. Briefly, cells were seeded in complete medium overnight in a 96-well plate at a density of 4 × 10³ cells per well. On the indicated day of culture, 20 μL of MTT (5.0 mg/mL) were added to each well and incubated for 4 h. MTT was incubated with 200 μL dimethyl sulfoxide per well to resolve the formed formazan crystals. The optical density (OD) was measured at a wavelength of 570 nm on a microplate reader (Bio-Rad, USA) every day for the next 5 days (days 2, 3, 4, 5, and 6). The proliferation rate (%) was calculated using the formula: OD_sample/OD_control × 100%. A cell growth curve was drawn according to the daily mean of OD values. Each experiment was performed in triplicate.

Western Blot Analysis

Cells were lysed in RIPA buffer (1% w/v Triton X-100, 0.1% w/v SDS, 10 mM Tris-HCl pH 7.5, 150 mM NaCl, 1% sodium deoxycholate) supplemented with protease and phosphatase inhibitors (Roche, Switzerland). Equal amounts of protein (40 μg/lane) were loaded and resolved via 12% SDS-PAGE (50 V, 3 h) and then transferred to a nitrocellulose membrane (Perkin Elmer, USA). Membranes were blocked with 5% non-fat milk in Tris-buffered saline for 1 h at room temperature and incubated overnight serially at 4°C with the following primary antibodies: GAPDH (loading control; ab8245; 1:500; Abcam), HNF3β (ab108422; 1:1,000; Abcam), IL6 (ab6672; 1:500; Abcam), JAK1 (ab47435; 1:500; Abcam), and STAT3 (ab68153; 1:1,000; Abcam). The membranes were incubated with the corresponding horseradish peroxidase-conjugated secondary goat-anti-rabbit (1:5,000; sc-2004, Sigma, USA) and goat-anti-mouse (sc-2005) antibodies (1:5,000; Sigma, USA) at room temperature for 1 h. Finally, the blots were visualized with enhanced chemiluminescence (Amersham Biosciences, USA) in according with the manufacturer's protocol and detected using an Image Quant LAS 4000 Mini Bimolecular Imager (GE Healthcare, USA). Each experiment was performed in triplicate.

Mouse Xenograft Model

The Animal Ethics Committee of the Fourth Military Medical University approved the tumor formation experiment with mice. All the mice were obtained from the Chinese Academy Shanghai Experimental Center, housed under specific pathogen-free conditions, and maintained in an animal facility in accordance with the guidelines of the Johns Hopkins University Animal Care and Use Committee. Eighteen 5-week-old male BALB-C/nude mice were randomized into three groups (each group consisted of six mice): lenti-HNF3β, lenti-GFP, and normal controls. Subcutaneous tumors were established by injecting SW480 cells (5 × 10⁷) per mouse to the dorsal right flank. The tumor diameter was measured every 2 days upon their appearance. Tumor volume was calculated using the length and width of the tumor. The mice were euthanized 6 weeks after the subcutaneously injection, and the tumors were excised and weighed.

Statistical Analysis

All statistical analyses were performed using IBM SPSS version 22.0 and GraphPad Prism 5 software. Continuous variables are shown as mean ± standard deviation. Categorical variables are stated as number (percentage). The chi-squared test or Fisher's exact test and unpaired Student's t-test, or one-way analysis of variance or Mann–Whitney U-tests were used to assess possible correlations between HNF3β in tumor samples and clinicopathological factors, depending on variable type and the data distribution. Survival analysis was performed by the Kaplan–Meier method and differences in survival were compared using the log-rank test. The Cox regression model was performed to test the independent prognostic value. P values were two-sided and P < 0.05 was considered statistically significant.

Results

HNF3β Protein Expression in Colonic Tissue Samples

Based on the staining immunoreactivity score, in 14.37% (25/174) of colon tumors no HNF3β protein expression was observed, that is, staining was negative. The remaining tumor samples were positive for HNF3β protein staining, and 41.95% (73/174), 37.93% (66/174), and 5.75% (107/174) were categorized as weak (+), moderate (++), and strong (+++). Taken together, weak or negative HNF3β protein expression was observed in 56.32% (98/174) of colon tumor samples (Figure 1Ae,f), while normal colorectal tissues examined in this study (6 samples) showed strong positive HNF3β staining (Figure 1Aa). Furthermore, all inflamed tissues (8 samples) were also positive for HNF3β staining (Figure 1Ab) and positive HNF3β protein expression was also detected in 92.86% (13/14) of low-grade intraepithelial neoplasias (Figure 1Ac). Furthermore, none of the examined high-grade intraepithelial neoplasias (18 samples) were positive for HNF3β expression (Figure 1Ad). The difference in HNF3β protein expression between the non-cancerous colon mucosal tissues and colon cancer tissues was statistically significant (P < 0.0001; Figure 1Ba). Furthermore, low HNF3β expression was detected in 0/5 stage 0, 10/24 (41.67%) stage I, 21/49 (42.86%) stage II, 58/87 (66.67%) stage III, and 9/9 (100%) stage IV patients (Figure 1Ba,b), showing a trend toward decreased HNF3β expression in advanced tumor stages (III and IV; P = 0.0012). This observation indicates that HNF3β low expression is associated with tumor progression and therefore HNF3β may have a tumor suppressor role in colon cancer.

FIGURE 1

Figure 1. (A) Immunohistochemical analyses of HNF3β protein expression in colon cancer and non-cancerous colon tissues. (a) Normal colon tissue; (b) inflammation tissue; (c) low-grade intraepithelial neoplasia tissue; (d) high-grade intraepithelial neoplasia tissue; (e) colon cancer tissue; and (f) mucinous adenocarcinoma tissue. (B) HNF3β protein expression in colon cancer and non-cancerous colon tissues. (a) NO, normal colon tissue; IN, inflammation; LGIN, low-grade intraepithelial neoplasia; HGIN, high-grade intraepithelial neoplasia; CA, cancer. (b) HNF3β protein expression in different colon cancer stages. (C) Survival rates of patients with low (dashed lines) or high (bold lines) HNF3β protein expression. (a) OS; (b) PFS.

HNF3β Protein Expression and Clinicopathological Characteristics

HNF3β protein expression was associated with the progression of tumors from stage 0 to IV (P = 0.001) and with increasing depth of invasion (T classification) from Tis+T1 to T4 (P = 0.004; Table 1). Similar trend was observed in the mix population of TCGA colon cancer data sets (Supplementary Figure 1). In addition, low HNF3β expression was more frequently detected in tumors with regional lymph node and distant metastasis (P = 0.007 and P = 0.048, respectively). No correlation was observed between HNF3β protein expression and other clinicopathological characteristics, including patient's gender (P = 0.673), age (P = 0.495), tumor histological type (mucinous adenocarcinoma cf. adenocarcinoma, P = 0.601), and differentiation status (P = 0.107). Collectively these results again suggest that decrease in HNF3β protein expression may have a role in colorectal tumorigenesis and progression.

Decreased HNF3β Protein Expression Is Closely Related to Poor Prognosis of CRC Patients

During the entire follow-up period, 56.90% (99/174) of the colon cancer patients died and the postoperative median OS time of all recruited patients was 35 months (95% confidence interval [CI] 29.97–40.03). Nine patients had metastases at the time of the surgery. Of 165 remaining patients, 57.58% (95/165) had local recurrence or distant metastasis after surgery and the median PFS time of 30 months (95% CI 13.23–46.77). The 5-year OS rate of 174 patients was 43.10% and the 5-year PFS rate of 165 patients was 42.42%.

The Kaplan–Meier survival analysis and log-rank test revealed that colorectal cancer patients with low tumor HNF3β protein expression tended to have worse 5-year OS and PFS than patients with high tumor HNF3β protein expression (P < 0.001, both; Table 2). The postoperative median OS time of patients with low HNF3β expression was 30 months (95% CI 26.77–33.23) while that of patients with high HNF3β expression was 35 months (95% CI 29.97–40.03). The median PFS time of patients with low HNF3β expression in the tumor was 20 months (95% CI 16.45–23.55) while that of patients with high HNF3β tumor expression was 30 months (95% CI 13.23–46.77). The cumulative 5-year OS and PFS of patients with low HNF3β expression was lower than that of patients with high HNF3β expression (OS, 30.6% cf. 59.2%; PFS, 29.2% cf. 57.9%) (Figure 1C).

TABLE 2

Table 2. Univariate analysis of prognostic parameters affecting 5-year survival of HNF3β in CRC patients, n/N (%).

Univariate analyses revealed that a low level of HNF3β protein (P < 0.001), UICC clinical stages (P < 0.001), depth of invasion (T status; P < 0.001), positive lymph node invasion (P < 0.001), and distant metastasis (P < 0.001) were statistically significant predictors for poor OS. Furthermore, HNF3β protein expression (P < 0.001), UICC clinical stages (P < 0.001), depth of invasion (T status, P < 0.001), and lymph node invasion (P < 0.001) were associated with poor PFS (Table 2). However, there was no significant correlation between patients' gender, age, histological subtype, or differentiation status with either OS or PFS (Table 2).

In the multivariate analysis, the Cox proportional hazards model revealed that HNF3β expression [hazard ratio (HR) = 0.601, 95% CI, 0.387–0.933, P = 0.023], distant metastasis (HR = 5.89, 95% CI, 1.028–33.772, P = 0.046), and UICC clinical stage (HR = 3.737, 95% CI, 1.611–8.668, P = 0.002) were independent prognostic factors for OS. Finally, HNF3β expression (HR = 0.606, 95% CI, 0.392–0.937, P = 0.024), lymph node metastasis (HR = 0.171, 95% CI, 0.032–0.899, P = 0.037), and UICC clinical stage (HR = 9.935, 95% CI, 2.214–44.574, P = 0.003) were independent prognostic factors for PFS in colon cancer patients (Table 3).

TABLE 3

Table 3. Multivariate analysis of colorectal cancer patients' OS and PFS.

Collectively these findings suggest that low HNF3β protein expression in colorectal tumors could be closely related to poor prognosis of colorectal cancer patients.

HNF3β Overexpression Suppresses the Growth, Migration, and Invasion Ability of SW480 Cells

Since we identified a correlation between HNF3β low protein expression in tumor samples and clinicopathological parameters of colon cancer invasion and metastasis, we hypothesized that HNF3β may play a tumor suppressive role in colon cancer progression. By quantitative PCR, we identified HNF3β expression in several colon cell lines and found that it was lowest in SW480 cells (Supplementary Figures 2, 3). Therefore, MTT, transwell migration and invasion, and wound healing assays were performed in SW480 cells to explore the effect of HNF3β on cell growth, migration, and invasion (Figure 2A). Compared with the lenti-GFP transfected and non-transfected cells, the proliferation of the lenti-HNF3β-transfected SW480 cells was significantly decreased from day 3 to day 6 (P < 0.05). Furthermore, the migration and invasion abilities were significantly suppressed by HNF3β overexpression in transfected cells, compared with cells transfected with lenti-GFP the non-transfected cells (migration, 12.67 ± 1.37 cf. 23.33 ± 1.37 cf. 23.83 ± 1.47, P < 0.0001; invasion, 9.67 ± 1.21 cf. 19.00 ± 1.41 cf. 20.17 ± 2.14, P < 0.0001; Figures 2B–G).

FIGURE 2

Figure 2. (A) MTT assay assessing the proliferative effects of HNF3β overexpression in SW480 cells P < 0.05; ^#control cf. HNF3β; *GFP cf. HNF3β. (B–G) Effect of HNF3β overexpression on migration, invasion, and wound healing ability of SW480 cells.

The wound-healing assay also showed that cell migration was restrained by HNF3β overexpression compared with lenti-GFP and non-transfected SW480 cells at 12 and 24 h (mean relative migration distance: at 12 h, 5.33 ± 1.63 cf. 20.50 ± 2.35 cf. 22.67 ± 2.34, P < 0.0001; at 24 h: 20.00 ± 3.16 cf. 50.83 ± 3.76 cf. 55.00 ± 4.47, P < 0.0001; Figures 2D–F). Consistently, these results showed that HNF3β overexpression significantly impeded the proliferation, migration, and invasion ability of SW480 cells.

HNF3β Overexpression Inhibits the Growth of Subcutaneous Xenograft Tumors in Nude Mice

To examine the potential inhibitory effect of HNF3β overexpression on colon cancer growth in vivo, we set up a subcutaneous xenograft tumor model by subcutaneously injecting lenti-HNF3β SW480, lenti-GFP SW480, or SW480 cells in nude mice. The average final tumor volume and weight in the lenti-HNF3β group (0.55 ± 0.05 cm³ and 0.32 ± 0.08 g) was significantly suppressed compared with the lenti-GFP group (1.37 ± 0.08 cm³ and 0.68 ± 0.13 g) or control group (1.38 ± 0.07 cm³ and 0.73 ± 0.12 g; P < 0.0001 and P < 0.0001; Figure 3). Collectively, these results indicate that HNF3β overexpression in tumor cells may have a suppressive role in the progression of colon cancer.

FIGURE 3

Figure 3. Effect of HNF3β overexpression on xenograft tumor growth. (A) The tumors removed from mice. (B) Tumor weight. (C) The average tumor volume over the time of the experiment. P < 0.05; #control cf. HNF3β; *GFP cf. HNF3β.

HNF3β Overexpression Inhibited IL6, JAK1, and STAT3 Expression in SW480 Cells

To understand the underlying molecular mechanism through which HNF3β inhibits the growth, migration, and invasion of colon cancer cells, we screened the pathways that HNF3β was probably involved in using western blot. We found that HNF3β overexpression accompanied with decrease in expression of several key genes in JAK-STAT signaling. These included the following: IL6 (control cf. GFP cf. HNF3β, 0.36 ± 0.05 cf. 0.35 ± 0.04 cf. 0.21 ± 0.03, P = 0.0237); JAK1 (control cf. GFP cf. HNF3β, 0.37 ± 0.02 cf. 0.32 ± 0.03 cf. 0.13 ± 0.02, P = 0.0004); and STAT3 (control cf. GFP cf. HNF3β, 0.37 ± 0.02 cf. 0.32 ± 0.04 cf. 0.19 ± 0.02, P = 0.0044; Figure 4). However, no direct interaction between HNF3β and STAT3 was detected by co-immunoprecipitation (Supplementary Figure 5). To explore if HNF3β was involved in interferon γ (IFNγ)-STAT1 signaling so as to affect the STAT3 signal pathway, we have examined the IFNγ expression by quantitative PCR and ELISA. Non-significant change was found in HNF3β transfected cells comparing to the control (Supplementary Figure 4). Analysis of HNF3β involved pathways in STRING database showed that various of molecules correlated to STAT3 signaling, such as AKT1, mTOR, GSK3β, etc., and were likely to interact with HNF3β (Supplementary Figure 6).

FIGURE 4

Figure 4. Influence of HNF3β overexpression on protein expression of key JAK/STAT signaling in SW480 cells. (A) western blot analysis of IL6, JAK1, and STAT3 protein expression. (B–D) relative protein expression of IL6, JAK1, and STAT3 in non-transfected, empty vector transfected, and lenti-HNF3β tranfected SW480 cells P < 0.05; ^#control cf. HNF3β; *GFP cf. HNF3β.

Taken together, these findings suggest that HNF3β suppresses the growth, migration, and invasion of colon cancer cells by JAK-STAT3 signaling.

Discussion

HNF3β has been described in the scientific literature as one of the most important transcription factors involved in gastrointestinal tract differentiation. It has also recently been implicated in tumorigenesis and thus we decided to examine its status in colorectal cancer.

The HNF3 gene family, including HNF3α, HNF3β, and HNF3γ (also known as FOXA1, FOXA2, and FOXA3, respectively), has a highly conserved DNA binding domain, the winged helix domain, which binds to the DNA sequence of the target gene in a monomeric form (15). HNF3β is the founding member of the HNF3 family of transcription factors (16), and its low expression was reported previously in tumor tissues and cell lines. Indeed, over the years accumulating evidence from a number of studies has supported the view that HNF3β genetic alterations or suppression were closely related to the occurrence, progression, or metastasis of several kinds of tumors, including gastric, thyroid, and lung cancer (10, 11, 14). For instance, in the study by Basseres et al. (17), 24.2% (42/173) of non-small cell lung cancers had a HNF3β deletion and in 55.4% (96/173) of these tumors low HNF3β expression was observed. Halmos et al. (10) further confirmed that HNF3β expression in lung cancer cell lines was low compared with normal lung cells. Furthermore, HNF3β expression negatively correlated with the migration ability of lung cancer cells (17, 18), while HNF3β overexpression could induce growth reduction, promote proliferation arrest and apoptosis, and inhibit metastasis of lung cancers (10). Other studies have also shown that HNF3β expression was decreased in papillary thyroid carcinoma and melanoma cell lines (11, 19). HNF3β putatively has a role in the progression of Barrett's esophagus to esophageal adenocarcinoma (20). More recently, inactivation of HNF3β has been implicated in inflammation-mediated liver cancer carcinogenesis (21, 22). However, HNF3β expression and its function in colon cancer has not been systematically investigated.

We analyzed HNF3β expression in the public TCGA database, which consisted of a mixed population including only 12 Asian cases. It showed a non-significant trend of decreasing HNF3β expression with increasing tumor stages, suggesting that more data on HNF3β in East Asian population were required and of great value. Our results show that HNF3β expression was significantly decreased in colon cancer tissues, compared with normal mucous tissue. Furthermore, low HNF3β expression significantly correlated with advanced tumor stage, cancer invasion, regional lymph nodes, and distant metastasis. Colon cancer patients with HNF3β low expression had shorter 5-year OS and PFS compared with those with higher HNF3β expression. Moreover, HNF3β expression was an independent prognostic factor for colon cancer patients. Besides the findings obtained from the IHC analysis of colorectal tumors, we have also shown that forced expression of HNF3β significantly impeded the proliferation, migration, and invasion of colon cancer cells in vitro and in vivo.

The JAK /STAT signaling pathway is one of the most important that is involved in the control of inflammatory response, cell proliferation, apoptosis, telomerase activity, angiogenesis, tumor invasion, and metastasis (23). Some studies reported a regulatory interaction between STAT3 and HNF3β (24, 25). However, the correlation between HNF3β and STAT3 signal pathway has not been well documented and no related study was conducted in colon cancer. In the present study we have demonstrated that HNF3β suppression of colorectal cancer growth, migration, and invasion may be involved in this pathway and that HNF3β upregulation inhibited the expression of its molecules. Although no direct interaction was found between HNF3β and STAT3, analysis by STRING database showed that HNF3β was likely to regulate STAT3 signaling by a complex, dynamic network via such as AKT1, mTOR and/or GSK3β. In addition, aberrant Akt was found to contribute to maintaining stemness through a NF-kB/IL-6/STAT3 pathway in lung cancer cells (26). The expression of STAT3 can be reduced by mTOR blocking, which is closely related to the invasion and metastasis of liver cancer cells (27). STAT3 is also a target of GSK3β and GSK3β inhibition reduced STAT3 phosphorylation while higher GSK3β expression promoted esophageal squamous cell carcinoma progression through STAT3 in vitro and in vivo (28). Therefore, HNF3β probably regulates STAT3 signaling via the network including those and more other molecules, and inhibition of this inflammatory signaling pathway may be the mechanism through which HNF3β inhibits the malignant progression of colon cancer. STAT1 signaling that was mainly activated by IFNγ was known to regulate STAT3 signaling (29). However, our data showed that HNF3β was not involved in IFNγ-STAT1 pathway, which implied that the underlying mechanism of how HNF3β is involved in STAT3 signaling remains to be defined in future studies. Finally, gene expression in cell lines is not exactly the same as in tissues due to their different histological origins, microenvironment and tumor heterogeneity, nevertheless, observations obtained from clinical samples are usually more representative of the real conditions in vivo. In addition, examining expression of target genes in larger sets of normal tissue specimens would be of great use, however obtaining them is often difficult.

To the best of our knowledge, our study is the first to identify downregulated HNF3β as an independent poor prognostic factor in colon cancer. In addition, in this study we are the first to demonstrate its inhibitory role on malignant behavior of colon cancer in vitro and in vivo and identified the possible underlying molecular mechanism of its action through modulation of JAK-STAT signaling. Based on these data, we speculate that HNF3β plays a suppressive role in colon cancer tumorigenesis, and that it can be used as a potential new biomarker and therapeutic target in the control of colon cancer progression.

Data Availability Statement

The datasets generated and analyzed during the present study are available from the corresponding author on reasonable request.

Ethics Statement

The studies involving human participants were reviewed and approved by Research Ethics Committee of the Fourth Military Medical University. The patients/participants provided their written informed consent to participate in this study. The animal study was reviewed and approved by Animal Ethics Committee of the Fourth Military Medical University.

Author Contributions

ZH, GJ, and JW conceived and designed the study. JW, HL, and WW performed the experiments. JW wrote the paper. NZ and YW participated in the statistical analysis. WW, ZH, and JW reviewed and edited the manuscript. All authors read and approved the manuscript and agree to be accountable for all aspects of the research in ensuring that the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Funding

This work was supported by the grant from the National Natural Science Foundation of China (Key Program 81502401).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.01096/full#supplementary-material

References

1. Center MM, Jemal A, Smith RA, Ward E. Worldwide variations in colorectal cancer. CA Cancer J Clin. (2009) 59:366–78. doi: 10.3322/caac.20038

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2015. CA Cancer J Clin. (2015) 65:5–29. doi: 10.3322/caac.21254

CrossRef Full Text | Google Scholar

3. Schrem H, Klempnauer J, Borlak J. Liver-enriched transcription factors in liver function and development. Part I: the hepatocyte nuclear factor network and liver-specific gene expression. Pharmacol Rev. (2002) 54:129–58. doi: 10.1124/pr.54.1.129

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Richmond CA, Breault DT. Regulation of gene expression in the intestinal epithelium. Prog Mol Biol Transl Sci. (2010) 96:207–29. doi: 10.1016/B978-0-12-381280-3.00009-9

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Ang SL, Wierda A, Wong D, Stevens KA, Cascio S, Rossant J, et al. The formation and maintenance of the definitive endoderm lineage in the mouse: involvement of HNF3/forkhead proteins. Development. (1993) 119:1301–15.

PubMed Abstract | Google Scholar

6. Monaghan AP, Kaestner KH, Grau E, Schutz G. Postimplantation expression patterns indicate a role for the mouse forkhead/HNF-3 alpha, beta and gamma genes in determination of the definitive endoderm, chordamesoderm and neuroectoderm. Development. (1993) 119:567–78.

PubMed Abstract | Google Scholar

7. Sasaki H, Hogan BL. Differential expression of multiple fork head related genes during gastrulation and axial pattern formation in the mouse embryo. Development. (1993) 118:47–59.

PubMed Abstract | Google Scholar

8. Weinstein DC, Ruiz i Altaba A, Chen WS, Hoodless P, Prezioso VR, Jessell TM, et al. The winged-helix transcription factor HNF-3 beta is required for notochord development in the mouse embryo. Cell. (1994) 78:575–88. doi: 10.1016/0092-8674(94)90523-1

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Montgomery RK, Mulberg AE, Grand RJ. Development of the human gastrointestinal tract: twenty years of progress. Gastroenterology. (1999) 116:702–31. doi: 10.1016/S0016-5085(99)70193-9

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Halmos B, Basseres DS, Monti S, D'Alo F, Dayaram T, Ferenczi K, et al. A transcriptional profiling study of CCAAT/enhancer binding protein targets identifies hepatocyte nuclear factor 3 beta as a novel tumor suppressor in lung cancer. Cancer Res. (2004) 64:4137–47. doi: 10.1158/0008-5472.CAN-03-4052

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Akagi T, Luong QT, Gui D, Said J, Selektar J, Yung A, et al. Induction of sodium iodide symporter gene and molecular characterisation of HNF3 beta/FoxA2, TTF-1 and C/EBP beta in thyroid carcinoma cells. Br J Cancer. (2008) 99:781–8. doi: 10.1038/sj.bjc.6604544

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Song Y, Washington MK, Crawford HC. Loss of FOXA1/2 is essential for the epithelial-to-mesenchymal transition in pancreatic cancer. Cancer Res. (2010) 70:2115–25. doi: 10.1158/0008-5472.CAN-09-2979

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wang J, Zhu CP, Hu PF, Qian H, Ning BF, Zhang Q, et al. FOXA2 suppresses the metastasis of hepatocellular carcinoma partially through matrix metalloproteinase-9 inhibition. Carcinogenesis. (2014) 35:2576–83. doi: 10.1093/carcin/bgu180

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhu CP, Wang J, Shi B, Hu PF, Ning BF, Zhang Q, et al. The transcription factor FOXA2 suppresses gastric tumorigenesis in vitro and in vivo. Dig Dis Sci. (2015) 60:109–17. doi: 10.1007/s10620-014-3290-4

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Kaestner KH, Knochel W, Martinez DE. Unified nomenclature for the winged helix/forkhead transcription factors. Genes Dev. (2000) 14:142–6. doi: 10.1101/gad.14.2.142

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Hannenhalli S, Kaestner KH. The evolution of Fox genes and their role in development and disease. Nat Rev Genet. (2009) 10:233–40. doi: 10.1038/nrg2523

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Basseres DS, D'Alo F, Yeap BY, Lowenberg EC, Gonzalez DA, Yasuda H, et al. Frequent downregulation of the transcription factor Foxa2 in lung cancer through epigenetic silencing. Lung Cancer. (2012) 77:31–7. doi: 10.1016/j.lungcan.2012.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Tang Y, Shu G, Yuan X, Jing N, Song J. FOXA2 functions as a suppressor of tumor metastasis by inhibition of epithelial-to-mesenchymal transition in human lung cancers. Cell Res. (2011) 21:316–26. doi: 10.1038/cr.2010.126

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yu KS, Jo JY, Kim SJ, Lee Y, Bae JH, Chung YH, et al. Epigenetic regulation of the transcription factor Foxa2 directs differential elafin expression in melanocytes and melanoma cells. Biochem Biophys Res Commun. (2011) 408:160–6. doi: 10.1016/j.bbrc.2011.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Watts JA, Zhang C, Klein-Szanto AJ, Kormish JD, Fu J, Zhang MQ, et al. Study of FoxA pioneer factor at silent genes reveals Rfx-repressed enhancer at Cdx2 and a potential indicator of esophageal adenocarcinoma development. PLoS Genet. (2011) 7:e1002277. doi: 10.1371/journal.pgen.1002277

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Xu L, Hui L, Wang S, Gong J, Jin Y, Wang Y, et al. Expression profiling suggested a regulatory role of liver-enriched transcription factors in human hepatocellular carcinoma. Cancer Res. (2001) 61:3176–81. doi: 10.1097/00002820-200104000-00012

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Liu M, Lee DF, Chen CT, Yen CJ, Li LY, Lee HJ, et al. IKKalpha activation of NOTCH links tumorigenesis via FOXA2 suppression. Mol Cell. (2012) 45:171–84. doi: 10.1016/j.molcel.2011.11.018

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Thomas SJ, Snowden JA, Zeidler MP, Danson SJ. The role of JAK/STAT signalling in the pathogenesis, prognosis and treatment of solid tumours. Br J Cancer. (2015) 113:365–71. doi: 10.1038/bjc.2015.233

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Ibrahim S, Dayoub R, Krautbauer S, Liebisch G, Wege AK, Melter M, et al. Bile acid-induced apoptosis and bile acid synthesis are reduced by over-expression of Augmenter of Liver Regeneration (ALR) in a STAT3-dependent mechanism. Exp Cell Res. (2019) 374:189–97. doi: 10.1016/j.yexcr.2018.11.023

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Hao Y, Kuang Z, Jing J, Miao J, Mei LY, Lee RJ, et al. Mycoplasma pneumoniae modulates STAT3-STAT6/EGFR-FOXA2 signaling to induce overexpression of airway mucins. Infect Immun. (2014) 82:5246–55. doi: 10.1128/IAI.01989-14

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Malanga D, De Marco C, Guerriero I, Colelli F, Rinaldo N, Scrima M, et al. The Akt1/IL-6/STAT3 pathway regulates growth of lung tumor initiating cells. Oncotarget. (2015) 6:42667–86. doi: 10.18632/oncotarget.5626

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Pu X, Guo QX, Long HA, Yang CW. Effects of mTOR-STAT3 on the migration and invasion abilities of hepatoma cell and mTOR-STAT3 expression in liver cancer. Asian Pac J Trop Med. (2014) 7:368–72. doi: 10.1016/S1995-7645(14)60057-2

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Gao S, Li S, Duan X, Gu Z, Ma Z, Yuan X, et al. Inhibition of glycogen synthase kinase 3 beta (GSK3β) suppresses the progression of esophageal squamous cell carcinoma by modifying STAT3 activity. Mol Carcinog. (2017) 56:2301–16. doi: 10.1002/mc.22685

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Costa-Pereira AP, Tininini S, Strobl B, Alonzi T, Schlaak JF, Is'harc H, et al. Mutational switch of an IL-6 response to an interferon- gamma-like response. Proc Natl Acad Sci USA. (2002) 99:8043–7. doi: 10.1073/pnas.122236099

PubMed Abstract | CrossRef Full Text | Google Scholar