Abstract

This study aimed to examine HIF-2α, TWIST, and CXCR4 expression in papillary thyroid carcinoma (PTC) and assesses the association of their expression with clinicopathological indicators. HIF-2α, TWIST, and CXCR4 protein expression in 129 PTCs, 61 nodular hyperplasia, and 118 normal thyroid tissue specimens was analyzed using immunohistochemistry. The protein expression levels of these three molecules were upregulated in PTCs. High protein expression of HIF-2α, TWIST, and CXCR4 was significantly correlated with lymph node metastasis (LNM) ( ). Furthermore, HIF-2α, TWIST, and CXCR4 protein expression was correlated with one another. Concomitant high expression of these molecules had stronger correlation with LNM than did each alone ( for HIF-2α/TWIST, for HIF-2α/CXCR4, for TWIST/CXCR4, and for HIF-2α/TWIST/CXCR4). Additionally, HIF-2α, TWIST, and CXCR4 mRNA expression were assessed in 30 PTCs, 10 nodular hyperplasia, and 10 normal thyroid tissue specimens using real-time RT-PCR. TWIST and CXCR4 mRNA expression levels were up-regulated in PTCs, and high mRNA expression of TWIST and CXCR4 was significantly correlated with LNM ( and , resp.). These results demonstrated that the evaluation of HIF-2α, TWIST, and CXCR4 expression in PTC may be useful in predicting the risk of LNM.

1. Introduction

Papillary thyroid carcinoma (PTC) accounts for 80% of thyroid malignancy and is characterized by slow growth and an excellent prognosis. However, some cases show relatively early recurrence, severe invasion, lymph node metastasis (LNM), or distant metastasis [1]. It is important to identify the characteristics of PTC that have a high risk for invasion and metastasis.

Intratumoral hypoxia is an independent indicator of poor prognosis and contributes to a more aggressive tumor phenotype [2]. Tumor adaptation to hypoxia is predominantly regulated by two structurally related hypoxia inducible factors (HIFs), HIF-1α and HIF-2α [3, 4]. They activate the expression of genes involved in tumor growth, metabolism, angiogenesis, and metastasis. Although highly homologous, HIF-1α and HIF-2α play critical but nonoverlapping roles in tumor progression [5, 6]. Recent studies have shown that TWIST and CXCR4 are two direct target genes of HIF-2α, HIF-2α can activate the expression of TWIST and CXCR4, which subsequently promote invasion and metastasis of tumor cells [79].

TWIST is a class II member of basic helix-loop-helix (bHLH) transcription factors and is thought to regulate epithelial-mesenchymal transition (EMT) through the downregulation of key proteins that maintain epithelial cell characteristics and upregulation of proteins that confer a mesenchymal phenotype [10]. TWIST is overexpressed in many types of tumors including breast, gastric, hepatocellular, prostate, and bladder cancers. Its upregulation correlates with high cancer aggressiveness and poor patient survival rate [11, 12]. However, at present, studies barely assess the expression of TWIST and its role in PTC.

Chemokine receptors are transmembrane proteins that interact with specific chemokine ligands, resulting in G-protein-coupled signal transduction leading to chemotaxis or directional movement along a chemical gradient. Recently, functional chemokine receptors have shown to be expressed by a large number of human malignancies, leading to the hypothesis that chemokines may stimulate proliferation, chemotaxis, and site-directed metastasis of tumor cells [13, 14]. The CXC chemokine receptor 4 (CXCR4) interacts specifically with its chemokine ligand, stromal cell-derived factor 1 (SDF-1/CXCL12), to exert proliferative and chemotactic effects in CXCR4-expressing cancer cells, accounting for an association of CXCR4 with metastasis in numerous human tumors, including lung, breast, prostate, colon, and thyroid cancer [1519].

Given their association with tumor invasion and metastasis, the purpose of this study was to examine HIF-2α, TWIST, and CXCR4 expression in PTC, to assess association of HIF-2α, TWIST, and CXCR4 expression with several clinicopathological indicators, and to evaluate potential usefulness of the three molecules in prediction for invasion and metastasis of PTC.

2. Materials and Methods

2.1. Case Selection and Tissue Sample Preparation

Tumor specimens for immunohistochemistry were obtained from 129 PTC patients who underwent thyroidectomy in the Department of Surgery, the First Affiliated Hospital, Chongqing Medical University, China between January 2010 and January 2013. There were 32 men and 97 women with a median age of 45 years. According to histopathologic diagnosis, there were 80 classic PTC, 21 follicular variant of PTC, 15 tall cell variant of PTC, and 13 oncocytic variant of PTC. The size of primary tumor ranged from 0.3 to 6.0 cm ( ). According to AJCC classification [20], there were 69 patients with stage I and stage II 60 with stage III and stage IV. Sixty-eight patients were confirmed to have LNM. Besides, benign thyroid disease specimens were obtained from 61 patients with nodular hyperplasia. One hundred and eighteen normal thyroid tissues were taken from the contralateral lobe of PTC specimens, which exhibit apparently normal morphology as a control. The study protocol was approved by the Research Ethics Committee of Chongqing Medical University and informed consent was obtained from all patients.

Tumor specimens for real-time RT-PCR were obtained from 30 PTC patients between January 2012 and January 2013, including 15 PTCs without LNM and 15 PTCs with LNM. The benign thyroid disease specimens were obtained from 10 patients with nodular hyperplasia. For controls, 10 normal thyroid tissue specimens were used. All specimens were immediately snap-frozen in liquid nitrogen and stored at −80°C up to subsequent real-time RT-PCR.

2.2. Tissue Microarray

Formalin-fixed, paraffin-embedded blocks were routinely prepared from surgical specimens of PTC, nodular hyperplasia, and normal thyroid tissue. Representative areas containing tumor, nodular hyperplasia, or normal thyroid tissue were identified by a pathologist. Duplicate tissue cores with a diameter of 0.6 mm were taken from each specimen (Beecher Instruments, Silver Springs, USA) and arrayed on a recipient paraffin block, using standard procedures [21]. Serial 5-μm-thick sections were cut with a Leica microtome (Leica Microsystems; Wetzlar, Germany) and mounted onto polylysine-coated slides.

2.3. Immunohistochemical Staining

Sections from TMA blocks were dewaxed and hydrated. Antigen retrieval was achieved by microwaving in 0.01 M citrate buffer (pH 6.0) for 10 min. After microwave treatment, the slides were treated with 3% hydrogen peroxide for 30 min to block the endogenous peroxidase and followed by blocking with 10% normal goat serum (50062Z, Invitrogen, USA) in PBS at room temperature for 1 h. The slides were then incubated overnight at 4°C in the primary rabbit polyclonal anti-HIF-2α antibody (1 : 50 dilution, ab20654; Abcam, USA), anti-TWIST antibody (1 : 50 dilution, ab49254; Abcam, USA), or anti-CXCR4 antibody (1 : 50 dilution, ab7199; Abcam, USA). For negative isotype controls, the sections were incubated in rabbit immunoglobulin G (1 : 1000, NI01-100UG; Merck Millipore, Germany). After defrosting at 37°C for 30 min, the slides were washed with PBS and incubated with a secondary biotinylated goat-anti-rabbit antibody (ZB-2010; Zhongshan Golden Bridge Biotechnology, China) for 30 min, peroxidase-labeled streptavidin (ZB-2404; Zhongshan Golden Bridge Biotechnology, China) for 20 min, and diaminobenzidine chromogen substrate (Sigma, USA) for 5 min. Slides were counterstained with hematoxylin, dehydrated, and mounted.

2.4. Immunohistochemical Scoring

A semiquantitative assessment of immunohistochemical (IHC) scoring was performed by two observers blinded to the diagnosis. The IHC score was assigned based on staining intensity and percentage of positive cells. The intensity score was assigned as 0 (no staining), 1 (weak staining), 2 (moderate staining), and 3 (strong staining). The proportion score was assigned as 0 (<5% positive cells), 1 (6–25% positive cells), 2 (26–50% positive cells), 3 (51–75% positive cells), and 4 (>75% positive cells). Multiplication of the intensity and percentage scores gave rise to the final staining score: 0 (negative), + (1–4), ++ (5–8), and +++ (9–12). For statistical analysis, a final staining score of negative or + was combined into the low expression group, and a final staining score of ++ or +++ was combined into the high expression group.

2.5. RNA Extraction, Reverse Transcription, and Real-Time PCR

Total RNA was extracted from frozen thyroid tissues using TRIzol reagent (Invitrogen, Camarillo, CA, USA), and residual genomic DNA was eliminated by DNase I digestion (Ambion, USA). RNA purity was confirmed by spectrophotometry. Total RNA was reverse transcribed to cDNA by using SuperScript III Reverse Transcriptase (Invitrogen, USA) according to the manufacturer’s protocol. The final cDNA product was amounted to 25 μL and stored at −80°C.

Real-time PCR was performed by using SYBR-Green real-time PCR method on the ABI-Prism 7000 sequence detector (Applied Biosystems, USA). The primers are shown in Table 1. The predicated product size of the primers for HIF-2α, TWIST, and CXCR4 was 196 bp, 207 bp, and 109 bp, respectively. Quantities of gene specific mRNA expression were determined by the CT method. Samples were analyzed in triplicate. Average threshold cycle (CT) values for glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were used as an internal calibrator. The method was used for relative quantitation [22]. Results are presented as the mean ± standard deviation of three independent experiments. The real-time PCR mix was made on the basis of the prescription from the supplier: 6 μL sterile water, 1 μL sense and 1 μL antisense primers, 10 μL Platinum SYBR Green qPCR SuperMix-UDG w/ROX (Invitrogen, USA), and 2 μL target cDNA in a total volume of 20 μL. Run conditions were 50°C for 2 min, 95°C for 10 min, followed by 40 cycles at 95°C for 15 s and 60°C for 1 min.

2.6. Statistical Analysis

Statistical analysis was performed using SPSS 18.0 statistical software. Data are presented as percentages and mean and standard deviation, according to the distribution. Significance was assessed using chi-square, Spearman rank, and Mann-Whitney U tests as appropriate, to compare the groups. value <0.05 was considered statistically significant.

3. Results

3.1. Immunohistochemical Expression of HIF-2α, TWIST, and CXCR4 in PTCs, Nodular Hyperplasia, and Normal Thyroid Tissues

HIF-2α, TWIST, and CXCR4 protein expression was examined by immunohistochemistry and illustrated in Figure 1. The immunoreactivities of HIF-2α and TWIST were detected in the cytoplasm and nuclei and CXCR4 in the cytoplasm and cell membrane. In nodular hyperplasia tissues, there were only a few follicular cells with weak staining for HIF-2α (a), TWIST (b), and CXCR4 (c). However, in PTCs, some cases had quite a few tumor cells with moderate staining for these three molecules (d)–(f), and some cases had a lot of tumor cells with strong staining for the three molecules (g)–(i). As shown in Tables 2 and 3, like the normal thyroid tissues, the majority of nodular hyperplasia tissues have negative or 1 IHC score; no cases showed high expression (≥5) of these three molecules. However, in PTCs, the majority of cases have ≥3 IHC score; high expression (≥5) was present in 64 (49.6%), 56 (43.4%), and 71 (55.0%) of 129 cases for HIF-2α, TWIST, and CXCR4, respectively. The differences in HIF-2α, TWIST, and CXCR4 protein expression levels between PTCs and normal thyroid tissues as well nodular hyperplasia tissues were statistically significant ( ).

3.2. Correlation of HIF-2α, TWIST, and CXCR4 Protein Expression with Clinicopathological Features in PTCs

The correlation of HIF-2α, TWIST, and CXCR4 protein expression with clinicopathological data was assessed by chi-square test and summarized in Table 3. There were no significant differences in HIF-2α, TWIST, and CXCR4 protein expression between patients with different histologic subtypes of PTC ( , , and , resp.), between older (>45) and younger (≤45) patients ( , , and , resp.), between male and female patients ( , , and , resp.), between patients with larger (>2.3) and smaller (≤2.3) tumor size ( , , and , resp.), and between patients with higher (III-IV) and lower stage (I-II) ( , , and , resp.). However, HIF-2α, TWIST, and CXCR protein expression was significantly correlated with LNM; PTC patients with LNM showed higher protein expression of these three molecules than those without LNM ( for all the three molecules).

3.3. Correlation of HIF-2α, TWIST, and CXCR4 Protein Expression with One Another in PTCs

The correlation of HIF-2α, TWIST, and CXCR4 protein expression with one another was assessed by Spearman rank test. As shown in Table 4, 40/129 PTCs showed high expression and 49/129 displayed low expression for both HIF-2α and TWIST. The correlation between HIF-2α and TWIST expression was statistically significant ( , ). Similarly, there was a statistically significant correlation between expression of HIF-2α and CXCR4 ( , ). For both HIF-2α and CXCR4, 48/129 PTCs showed high expression. In addition, high expression for both TWIST and CXCR4 was present in 44/129 PTCs. A significantly positive correlation ( , ) was also present between expression of TWIST and CXCR4.

3.4. Association of Concomitant HIF-2α, TWIST, and CXCR4 High Expression with LNM in PTCs

Given that HIF-2α, TWIST, and CXCR4 protein expression were correlated with one another and statistical analysis showed that the incidence of LNM tends to be higher in PTCs with high protein expression of HIF-2α, TWIST, and CXCR4, we further evaluated the correlation of LNM with concomitant high expression of HIF-2α/TWIST, HIF-2α/CXCR4, or TWIST/CXCR4. As shown in Table 5, the incidence of LNM is significantly higher in patients (87.5%) with high expression of HIF-2α/TWIST than in those patients (67.5%) with high expression of only one of these two molecules or in those patients (12.2%) without high expression for either of these two molecules. Similar results were observed in PTCs with high expression of HIF-2α/CXCR4 and TWIST/CXCR4. There were statistically significant differences in the incidence of LNM between patients with high expression of only one and any two of the three molecules ( for HIF-2α/TWIST, for HIF-2α/CXCR4, and for TWIST/CXCR4). In addition, statistical analysis showed that concomitant high expression of all the three molecules is significantly associated with LNM as compared with cases not showing such expression ( ). As demonstrated in Figures 2(a)2(c) are representative of PTC without LNM showing only HIF-2α high expression and low expression of the other two molecules, TWIST and CXCR4; (d)-(e) are representative of PTC with LNM showing high expression of all the three molecules, HIF-2α, TWIST, and CXCR4, respectively.

3.5. Correlation of HIF-2α, TWIST, and CXCR4 mRNA Expression with LNM in PTCs

To compare gene expression of HIF-2α, TWIST, and CXCR4 in PTCs without and with LNM, fifteen PTCs without LNM and another fifteen PTCs with LNM were collected to analyze HIF-2α, TWIST, and CXCR4 mRNA levels using real-time RT-PCR. Ten nodular hyperplasia and ten normal thyroid tissues were used for comparison and as a control. As shown in Table 6, mean HIF-2α mRNA levels were -fold in normal thyroid tissues, -fold in nodular hyperplasia tissues, -fold in PTCs without LNM, and -fold in PTCs with LNM. There was no statistically significant difference in HIF-2α mRNA levels between nodular hyperplasia and normal thyroid tissues ( ), between PTCs and normal thyroid tissues ( ), and between PTCs with and without LNM ( ). In contrast, TWIST and CXCR4 mRNA levels were significantly higher in PTCs compared with normal thyroid tissues ( for both), while there were no statistically significant differences in TWIST and CXCR4 mRNA levels between nodular hyperplasia and normal thyroid tissues ( and , resp.). Moreover, PTCs with LNM showed increased mRNA levels of TWIST and CXCR4 compared with those without LNM. The differences in TWIST and CXCR4 mRNA levels between PTCs with and without LNM were statistically significant ( and , resp.).

4. Discussion

Hypoxia is a common condition found in a wide range of solid tumors and has been increasingly recognized to play a central role in different stages of tumor progression [2]. Tumor adaptation to hypoxia is predominantly regulated by HIF-1α and HIF-2α, which activate the expression of genes involved in proliferation, metabolism, angiogenesis, and metastasis. While highly homologous, HIF-1α and HIF-2α have unique tissue distributions and play critical but nonoverlapping roles in tumor progression [5, 6]. HIF-2α promotes metastasis through regulation of critical factors controlling tumor cell metastatic potential, such as CXCR4 and TWIST [79]. To date, studies have shown that HIF-2α, TWIST, and CXCR4 are overexpressed in several human tumors such as clear cell renal carcinoma (ccRCC) [23], nonsmall cell lung cancer (NSCLC) [24], neuroblastoma [25], and breast, prostate, gastric, hepatocellular, colon, and bladder cancers [1118]. CXCR4 has been observed to be upregulated in thyroid cancer [19]. However, studies barely examined HIF-2α and TWIST expression in PTC; moreover, no study investigated simultaneously the expression of HIF-2α, TWIST, and CXCR4 and assessed correlation of their expression with clinicopathological features in PTC. In our present study, we examined HIF-2α, TWIST, and CXCR4 protein expression in PTCs, nodular hyperplasia and normal thyroid tissues using immunohistochemistry. The results demonstrated that no cases of normal thyroid tissue and nodular hyperplasia show high protein expression of HIF-2α, TWIST, and CXCR4. However, in PTCs, high protein expression was present in 49.6%, 43.4%, and 55.0% of cases for HIF-2α, TWIST and CXCR4, respectively. The differences in HIF-2α, TWIST, and CXCR4 protein expression between PTCs and normal thyroid tissues as well nodular hyperplasia were statistically significant ( ). Then we assessed the correlation of HIF-2α, TWIST, and CXCR4 protein expression with several clinicopathological indicators. We found that HIF-2α, TWIST, and CXCR4 protein expression were not associated with histologic subtype, gender, age, tumor size, and TNM stage. However, there was a significant correlation between LNM and single HIF-2α, TWIST, and CXCR4 protein expression. High protein expression of HIF-2α, TWIST, and CXCR4 was associated with positive LNM. These results suggested that HIF-2α, TWIST, and CXCR4 may play an important role in invasion and metastasis of PTC.

Additionally, our study, for the first time, demonstrated a significantly positive correlation between HIF-2α, TWIST, and CXCR4 expression in PTCs. HIF-2α expression is positively correlated with TWIST expression ( , ) and CXCR4 expression ( , ). This finding could be supported by several recent studies. Gort et al. reported that TWIST expression in human cancer cells is enhanced by hypoxia in a HIF-2α-dependent manner. TWIST is a direct target of HIF-2α, and HIF-2α can promote the invasion and migration of cancer cells by upregulation of TWIST [7]. Meanwhile, studies have also showed that HIF-2α can promote invasion and metastasis through upregulation of CXCR4 expression [8, 9]. In addition, a significantly positive correlation ( , ) was present between expression of TWIST and CXCR4. To date, there is no report to explain this positive correlation. It is necessary to further explore mechanisms underlying this correlation.

Given that HIF-2α, TWIST, and CXCR4 protein expression were positively correlated with one another and the expression of these single molecules was related to LNM, we subsequently evaluated the association of concomitant expression of HIF-2α, TWIST, and CXCR4 with LNM in PTCs. The results showed that concomitant expression of any two of these three molecules had stronger correlation with LNM than did each alone. Concomitant expression of all three molecules strongly correlates with LNM.

Lastly, we analyzed HIF-2α, TWIST, and CXCR4 mRNA expression levels in PTCs using real-time RT-PCR. The results demonstrated that TWIST and CXCR4 mRNA levels were significantly higher in PTCs than in normal thyroid tissues. Moreover, as their protein expression, TWIST and CXCR4 mRNA expression levels were also correlated with LNM; PTCs with LNM showed to have higher mRNA levels of TWIST and CXCR4 than those without LNM. However, there was no statistically significant difference in HIF-2α mRNA expression levels between PTCs, nodular hyperplasia, and normal thyroid tissues. Possibly, HIF-2α, like HIF-1α [3, 4], is predominantly regulated through hypoxia-dependent protein stabilization and genetic alterations, and under both normoxia and hypoxia, HIF-2α mRNA is also expressed constitutively in many cell types, including cancer cells.

5. Conclusions

In summary, our results, for the first time, demonstrated a positive correlation of HIF-2α, TWIST, and CXCR4 expression in PTCs. High expression of HIF-2α, TWIST, and CXCR4 was associated with LNM. Concomitant high expression of any two or all of the three molecules had stronger correlation with LNM than did each alone. Consequently, our results provide a possible basis for prediction of LNM in PTC. Future studies in larger sets of patients will be necessary to determine the utility of these molecules as biomarkers of tumor diagnosis and prognosis in PTC.

Conflict of Interests

The authors declare that they have no conflict of interests.

Acknowledgments

This work was supported by the Chongqing Natural Science Foundation of China (CSTC, 2011BA5038), the National Natural Science Foundation of China (81072183), the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Personnel Ministry [2011], and the Scientific Research Foundation of Chongqing Medical University (XBZD201002).