Next Article in Journal
Automatic Diabetic Foot Ulcer Recognition Using Multi-Level Thermographic Image Data
Next Article in Special Issue
Homeobox Gene Expression Dysregulation as Potential Diagnostic and Prognostic Biomarkers in Bladder Cancer
Previous Article in Journal
Skin Lesion Synthesis and Classification Using an Improved DCGAN Classifier
Previous Article in Special Issue
Artificial Intelligence in the Advanced Diagnosis of Bladder Cancer-Comprehensive Literature Review and Future Advancement
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diagnostics
Volume 13
Issue 16
10.3390/diagnostics13162636
diagnostics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differential Protein Expression Patterns of HOXA13 and HOXB13 Are Associated with Bladder Cancer Progression

by
Fee-Wai Chin
1,
Huzlinda Hussin
2,
De-Ming Chau
1,
Teng-Aik Ong
3,
Rosna Yunus
4,
Azad Hassan Abdul Razack
3,
Khatijah Yusoff
5,6,
Soon-Choy Chan
6,7,* and
Abhi Veerakumarasivam
1,6,8,*
1
Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, Universiti Putra Malaysia (UPM), Serdang 43400, Selangor, Malaysia
2
Department of Pathology, Faculty of Medicine and Health Sciences, Universiti Putra Malaysia (UPM), Serdang 43400, Selangor, Malaysia
3
Department of Surgery, Faculty of Medicine, Universiti Malaya, Kuala Lumpur 50603, Malaysia
4
Department of Pathology, Hospital Kuala Lumpur, Kuala Lumpur 50586, Malaysia
5
Department of Microbiology, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia (UPM), Serdang 43400, Selangor, Malaysia
6
Malaysia Genome and Vaccine Institute, National Institutes of Biotechnology Malaysia, Kajang 43000, Selangor, Malaysia
7
School of Liberal Arts, Science and Technology, Perdana University, Kuala Lumpur 50490, Malaysia
8
School of Medical and Life Sciences, Sunway University, Bandar Sunway 47500, Selangor, Malaysia
*
Authors to whom correspondence should be addressed.
Diagnostics 2023, 13(16), 2636; https://doi.org/10.3390/diagnostics13162636
Submission received: 30 May 2023 / Revised: 2 July 2023 / Accepted: 8 July 2023 / Published: 9 August 2023
(This article belongs to the Special Issue Advanced Diagnosis and Treatment Strategies for Bladder Cancer)
Browse Figures
Review Reports Versions Notes

Abstract

:
Bladder cancer is a common urological cancer and has the highest recurrence rate of any cancer. The aim of our study was to profile and characterize the protein expression of homeobox A13 (HOXA13) and homeobox B13 (HOXB13) genes in Malaysian bladder cancer patients. The protein expression of HOXA13 and HOXB13 in formalin-fixed paraffin-embedded (FFPE) bladder cancer tissues was determined by immunohistochemistry (IHC) analysis. The association between HOXA13/HOXB13 protein expression and demographic/clinicopathological characteristics of the bladder cancer patients was determined by chi-square analysis. Approximately 63.6% of the bladder cancer tissues harbored high HOXA13 expression. High HOXA13 expression was significantly associated with non-muscle invasive bladder cancer, lower tumor grade, higher number of lymph node metastases, and recurrence risk. In contrast, low HOXB13 expression (including those with negative expression) was observed in 71.6% of the bladder cancer tissues analyzed. Low HOXB13 expression was significantly associated with muscle-invasive bladder cancer, higher tumor stage, tumor grade, and metastatic risk. Both HOXA13 and HOXB13 protein expression were found to be associated with bladder tumorigenesis. The putative oncogenic and tumor suppressive roles of HOXA13 and HOXB13, respectively, suggest their potential utility as biomarkers in bladder cancer.

1. Introduction

Globally, more than half a million bladder cancer patients were diagnosed, and 212,536 deaths occurred in 2020 [1]. The predominant histologic type of bladder cancer is urothelial cell carcinoma (UCC), which accounts for more than 90% of the cases [2]. At diagnosis, approximately 70% of UCCs are non-muscle invasive bladder cancer (NMIBC; stages Ta–T1), while the remaining 30% are muscle invasive bladder cancer (MIBC; stages T2–T4) [3]. However, between 50 and 70% of NMIBC patients eventually develop tumor recurrence within five years of undergoing transurethral resection of the bladder tumor (TURBT). Of these, 10–20% of NMIBC patients will eventually progress to MIBC [4,5]. Due to the high risk of recurrence and disease progression, NMIBC patients require frequent and lifelong surveillance. Therefore, bladder cancer is associated with the highest medical cost per patient from first diagnosis to death as compared to other cancers [6,7].
Urine cytology and cystoscopy are routinely used for the diagnosis of bladder cancer. In addition, cystoscopy is also commonly used in surveillance follow-up. Urine cytology involves the microscopic examination of exfoliated cells/tumor cells from the voided urine of patients for the detection of bladder cancer. However, this test has low sensitivity for detecting low-grade tumors [8]. Cystoscopy is an invasive procedure that often causes significant discomfort and pain for patients. As a result, bladder cancer patients often have moderate to low compliance with surveillance follow-up [9,10]. Hence, the development of effective biomarkers is crucial so that cystoscopy is used only for confirmation upon the detection of a specific and sensitive biomarker(s) in the urine.
Bladder cancer is characterized by increasing genomic instability and the cohesive dysregulation of the transcriptome that influence cellular and clinical phenotypes, resulting in the identification of candidate genes that harbor significant biomarker potential [11,12,13]. The human homeobox gene family is a group of regulatory genes that encode embryogenesis-associated transcription factors. These genes include homeobox A13 (HOXA13) and homeobox B13 (HOXB13) [14,15]. HOXA13 is involved in the formation and development of the limbs, reproductive system, and urinary tract [16]. HOXB13 plays a crucial role in the development and maintenance of the skin [17].
Several studies have evaluated the utility of HOXA13 as a biomarker to detect bladder cancer using voided urine samples [18,19]. The detection of HOXA13 and CDC2 mRNA in voided urine had an accuracy of 80% in determining whether a particular UCC was of high stage or grade. Moreover, HOXA13 mRNA expression alone was able to detect low-stage or low-grade tumors with more than 90% specificity [20]. Their findings were supported by other studies [18,19] in which HOXA13 expression was elevated in low-stage tumors (pTa–pT1) of bladder cancer. In addition, a combination of HOXA13 and BLCA-4 mRNA expression had 80% sensitivity and specificity in distinguishing between low- and high-grade tumors [18,19]. Hence, HOXA13 has huge potential as a biomarker for the diagnosis and prognosis of bladder cancer.
HOXA13 overexpression has been reported in bladder, brain, esophageal, gastric, liver, and prostate cancers [21,22,23,24,25,26,27]. The overexpression of HOXA13 has been shown to promote tumor cell proliferation, migration, invasion, metastasis, and epithelial–mesenchymal transition, as well as inhibit apoptosis. Thus, HOXA13 may play a role as an oncogene in these cancers. In contrast, HOXB13 overexpression has been reported in liver, ovarian, and pancreatic cancers [28,29,30]. The overexpression of HOXB13 was shown to promote tumor cell proliferation and inhibit apoptosis, suggesting an oncogenic role for HOXB13 in these cancers. Nevertheless, some studies have reported that HOXB13 is downregulated in colorectal, kidney, and prostate cancers [31,32,33]. Induced expression of HOXB13 inhibited tumor cell proliferation and induced apoptosis, suggesting that HOXB13 acts as a tumor suppressor gene in these types of cancers. In this study, both HOXA13 and HOXB13 were selected after conducting a comprehensive review of the existing literature on the dysregulation of homeobox genes in bladder cancer.
Although previous studies have shown the association of HOXA13 and HOXB13 with bladder cancer as well as the potential utility of HOXA13 as a biomarker for bladder cancer [18,19,20], findings from those studies are far from conclusive and require further robust validation. As such, the role of HOXA13 and HOXB13 in bladder cancer is not well studied. Moreover, both genes have not been characterized in Malaysian bladder cancer patients. Therefore, our study aimed at profiling the protein expression of HOXA13 and HOXB13 in a cohort of formalin-fixed paraffin-embedded (FFPE) bladder cancer tissues from Malaysian bladder cancer patients. In addition, our study also aimed to determine the association of HOXA13/HOXB13 protein expression with various demographic and clinicopathological characteristics of bladder cancer patients.

2. Materials and Methods

2.1. Patients and Tissue Samples

Bladder cancer patients admitted to Hospital Kuala Lumpur, Malaysia, from 2002–2013 were involved in this study. Ethical approval was obtained from the Medical Research and Ethics Committee, Ministry of Health Malaysia (approval number NMRR-12-970-13096). In addition, written informed consent was obtained from the patients. Retrieval of 573 bladder cancer patients’ medical records was conducted from April 2013 to August 2016. The collected information includes patients’ demographic (gender, ethnicity, and age) and clinicopathological characteristics (tumor invasiveness, tumor stage, tumor grade, lymph node metastasis, metastasis, and tumor recurrence).
The FFPE bladder cancer tissue blocks of 164 patients were retrospectively collected from the Department of Pathology, Hospital Kuala Lumpur, and stored at room temperature. Histological examination was conducted to select tissue blocks with at least 80% tumor cellularity, whereas tissue blocks with large areas of tumor necrosis or an insufficient amount of tissue were excluded. A total of 110 bladder cancer patients were selected for this study, in which their tissue blocks were obtained from cystoprostatectomy, cystectomy, TURBT, cystectomy plus total abdominal hysterectomy, bilateral salpingo-oophorectomy (TAHBSO), and biopsy.

2.2. Tissue Sectioning

Sectioning of FFPE tissues was performed using a LEICA RM2135 microtome (Leica Biosystems, Wetzlar, Germany). The tissue blocks were sectioned into 4 µm thick sections. The tissue sections were placed in a Leica HI1210 water bath (Leica Biosystems, Germany) set at 42 °C and then fixed onto SuperFrost Plus® glass slides (Thermo Fisher Scientific, Waltham, MA, USA). The slides were dried overnight at room temperature and then stored at 4 °C until further analysis.

2.3. Hematoxylin and Eosin (H&E) Staining

H&E staining was carried out using a Tissue-Tek Prisma Automated Slide Stainer (Sakura Finetek, Torrance, CA, USA). The H&E-stained slides were examined by a pathologist to ensure that the selected tissue blocks contained sufficient tumor cellularity for immunohistochemistry (IHC) analysis.

2.4. IHC Analysis

A total of 110 and 102 FFPE bladder cancer tissues were used for the IHC analysis of HOXA13 and HOXB13, respectively. There was insufficient tissue material in eight tissue blocks after the completion of the HOXA13 IHC analysis. Hence, HOXB13 IHC analysis was not conducted on these samples. IHC analysis was performed using the Dako REALTM EnvisionTM Detection System, Peroxidase/DAB+, Rabbit/Mouse kit (Agilent Technologies, Santa Clara, CA, USA). FFPE human skin tissue was used as the positive control for HOXA13 IHC analysis, while FFPE human prostate cancer tissue was used as the positive control for HOXB13 IHC analysis. In addition, matched tissues stained without the primary antibody were used as negative controls. In every batch of IHC analysis, FFPE bladder cancer tissues were analyzed together with a positive and a negative control.
Tissue slides were deparaffinized by incubating at 60 °C for 1 h in an oven (Venticell, Munich, Germany), followed by immersing them twice in xylene for 5 min. The slides were rehydrated by immersing them in a series of decreasing concentrations of ethanol: 100%, 100%, 95%, 80%, and 70% for 5 min each. The slides were then rinsed under slow-running tap water for 5 min. Antigen retrieval was performed using a microwave oven (National, Kadoma, Japan), where the slides were immersed in 1× Tris-EDTA buffer (pH 9) and heated at high mode for 5 min until boiling, followed by defrosting mode for 10 min. The slides were then cooled at room temperature for 35 min. Next, a hydrophobic circle was drawn around the tissue section onto the glass slides using a Dako Pen (Agilent Technologies, Santa Clara, CA, USA). The slides were then washed thrice by immersing them in 1× TBS with 0.05% Tween 20 (TBST) buffer for 5 min each time. The slides were then incubated with 200 µL of Dako REALTM Peroxidase-Blocking Solution (Agilent Technologies, Santa Clara, CA, USA) for 10 min. The slides were then washed thrice with 1× TBST buffer for 5 min each time.
For HOXA13 IHC analysis, 200 µL of the rabbit polyclonal anti-HOXA13 primary antibody (ab106503; Abcam, Cambridge, UK) at an optimal 1:400 dilution (working concentration of 2.5 µg/mL) was added onto each slide and incubated at room temperature for 30 min. For HOXB13 IHC analysis, 200 µL of the rabbit monoclonal anti-HOXB13 primary antibody (ab201682; Abcam, Cambridge, UK) at an optimal 1:200 dilution (working concentration of 1.725 µg/mL) was added onto each slide and incubated at room temperature for 1 h. The same optimal dilution of the primary antibody was used for all the tissue samples analyzed. A negative control was prepared, whereby the primary antibody was replaced with Dako Antibody Diluent (Agilent Technologies, Santa Clara, CA, USA). The slides were then washed thrice with 1× TBST buffer for 5 min each time. Subsequently, 200 µL of the Dako REALTM EnvisionTM/HRP, Rabbit/Mouse (Agilent Technologies, Santa Clara, CA, USA), which is the secondary antibody, was added to each slide and incubated at room temperature for 30 min. The slides were then washed with 1× TBST buffer for 5 min each time.
Next, 200 µL of Dako REALTM DAB+ Chromogen (Agilent Technologies, Santa Clara, CA, USA), which was diluted with Dako REALTM Substrate Buffer (Agilent Technologies, Santa Clara, CA, USA) at 1:50 dilution, was added to each slide and incubated for 10 min. The slides were rinsed twice under slow-running tap water for 5 min each time. After that, the slides were counterstained with Dako Mayer’s Hematoxylin (Agilent Technologies, Santa Clara, CA, USA) at a 1:2 dilution for 30 s. The slides were then rinsed twice under slow-running tap water for 5 min each time. After the slides had been stained with the primary and secondary antibodies and counterstained, they were dehydrated by immersing them in a series of increasing concentrations of ethanol: 70%, 80%, 95%, 100%, and 100% for 3 min each time. The slides were then immersed in xylene twice for 3 min each time. Finally, the slides were mounted with DePex mounting medium (VWR, Lutterworth, UK), covered with glass coverslips (Marienfeld, Lauda-Konigshofen, Germany), and dried overnight at room temperature.

2.5. Scoring of Immunohistochemical Staining

The IHC-stained slides were examined under a light microscope (Olympus, Center Valley, PA, USA) and scored by a pathologist who was blinded to the patients’ demographic and clinicopathological characteristics. The slides were scored based on the percentage of tumor cells expressing positive staining and the intensity of staining. The final IHC score was calculated by multiplying the percentage of positive cells by the staining intensity. The IHC scoring system for HOXA13 was based on Hu et al. [24], as shown in Table 1. Low HOXA13 expression was determined for tissues with final IHC scores of 0, 1, 2, and 3, while high HOXA13 expression was determined for tissues with final IHC scores of 4, 6, and 9.
Meanwhile, the IHC scoring system for HOXB13 was based on Kristiansen et al. [34], as shown in Table 2. Low HOXB13 expression was determined for tissues with final IHC scores of 0, 1, 2, 3, and 4, while tissues with final IHC scores of 6, 8, 9, and 12 were determined to have high HOXB13 expression.

2.6. Statistical Analysis

Statistical analysis was performed using the IBM SPSS Statistics version 20 software (SPSS Inc., Chicago, IL, USA). The association between the expression of both proteins and the demographic or clinicopathological characteristics of the bladder cancer patients was determined using chi-square analysis. A p-value of less than 0.05 was considered statistically significant. Patients were excluded from the chi-square analysis if missing data were present for a particular clinicopathological characteristic that was being analyzed.

3. Results

3.1. Protein Expression of HOXA13 and HOXB13 in Bladder Cancer Tissues

The majority of the FFPE bladder cancer tissues (94.5%, 104/110) positively expressed HOXA13, while only 5.5% (6/110) of these tissues were negative for HOXA13. The HOXA13 protein was predominantly localized in the nucleus of the bladder cancer cells, but in some tissues, it was also localized in both the nucleus and cytoplasm (Figure 1).
Both H&E and IHC staining of HOXA13 are shown in Figure 2. For HOXA13 IHC analysis, 40.0% (44/110) of the tissues showed strong staining, 31.8% (35/110) showed moderate staining, and 22.7% (25/110) showed weak staining (Figure 2). It was found that 63.6% (70/110) of the tissues had high HOXA13 expression, while 36.4% (40/110) had low HOXA13 expression (including those with negative expression).
Positive HOXB13 expression was found in 57.8% (59/102) of the FFPE bladder cancer tissues. The remaining tissues (42.2%, 43/102) did not express HOXB13; there was no observable staining. When expressed, the HOXB13 protein was mainly localized in the nucleus of the bladder cancer cells, followed by nucleus and cytoplasm co-localization. In a small fraction of tissues, expression was limited to the cytoplasm only (Figure 3).
Both H&E and IHC staining of HOXB13 are shown in Figure 4. For HOXB13 IHC analysis, 24.5% (25/102) showed weak staining, 21.6% (22/102) showed moderate staining, and 11.8% (12/102) showed strong staining (Figure 4). Most of the tissues (71.6%, 73/102) had low HOXB13 expression (including those with negative expression), whereas the remaining (28.4%, 29/102) had high HOXB13 expression.

3.2. Protein Expression of HOXA13 and HOXB13 across Different Stages and Grades of Bladder Cancer

HOXA13 protein expression in 110 FFPE bladder cancer tissues across different tumor stages and grades is depicted in Figure 5A,B, respectively. High HOXA13 expression was observed in 78.3% (18/23) of pTa tumors, 73.1% (19/26) of pT1 tumors, 66.7% (14/21) of pT3 tumors, and 58.8% (10/17) of pT4 tumors, as compared to 39.1% (9/23) of pT2 tumors. Correspondingly, high HOXA13 expression was also observed in 85.0% (17/20) of grade 1 tumors and 67.2% (43/64) of grade 3 tumors, as compared to 38.5% (10/26) of grade 2 tumors.
HOXB13 protein expression in 102 FFPE bladder cancer tissues across different tumor stages and grades is depicted in Figure 6A,B, respectively. Low HOXB13 expression was observed in 87.0% (20/23) pT2 tumors, 94.4% (17/18) pT3 tumors, and 93.3% (14/15) pT4 tumors as compared to 42.9% (9/21) pTa tumors and 52% (13/25) pT1 tumors. Correspondingly, low HOXB13 expression was also observed in 68% (17/25) grade 2 tumors, and 83.1% (49/59) grade 3 tumors, as compared to 38.9% (7/18) grade 1 tumors.

3.3. Association of HOXA13 and HOXB13 Protein Expression with Demographic and Clinicopathological Characteristics of Bladder Cancer Patients

The association between HOXA13 protein expression and various patient demographics as well as clinicopathological characteristics is summarized in Table 3. HOXA13 protein expression was significantly associated with tumor invasiveness (p = 0.020), tumor grade (p = 0.003), lymph node metastasis (p = 0.027), and recurrence (p = 0.043). Based on the strength of the association, a strong association was found for tumor grade, while a moderate association was found for tumor invasiveness and lymph node metastasis. However, a weak association was found for recurrence. Notably, high HOXA13 expression was significantly associated with NMIBC, lower tumor grade, higher number of lymph node metastases, and recurrence risk (Table 3). In addition, HOXA13 protein localization was significantly associated with age at diagnosis (p = 0.001) and recurrence (p = 0.023) (Table 4).
Meanwhile, the association between HOXB13 protein expression and various patient demographics as well as clinicopathological characteristics is summarized in Table 5. HOXB13 protein expression was significantly associated with tumor invasiveness (p ≤ 0.001), tumor stage (p ≤ 0.001), tumor grade (p = 0.001), and metastasis (p = 0.005). There was a strong association for tumor invasiveness, tumor stage, and tumor grade, while there was a moderate association for metastasis. Notably, low HOXB13 expression was significantly associated with MIBC, higher tumor stage, higher tumor grade, and metastatic risk (Table 5). In addition, HOXB13 protein localization was significantly associated with tumor invasiveness (p ≤ 0.001), tumor stage (p ≤ 0.001), tumor grade (p = 0.032), and metastasis (p = 0.001). Moreover, nuclear HOXB13 expression was significantly associated with NMIBC, as a high frequency of nuclear localization was observed in NMIBC, indicating nuclear HOXB13 expression is able to distinguish NMIBC from MIBC (Table 6).

4. Discussion

Our study demonstrates that the HOXA13 protein is heterogeneously expressed in bladder cancer tissues with a predominantly nuclear localization. This finding is in contrast to that of Hu et al. [24], who reported that HOXA13 was mainly localized in the cytoplasm of bladder cancer cells. The cytoplasmic localization of HOXA13 might be due to the modulation of nuclear localization signals, the nuclear export sequence, and its interaction with cytoplasmic anchoring factors [35,36]. It is worth noting that the HOXA13 protein is a transcription factor that is synthesized in the endoplasmic reticulum, which is located in the cytoplasm. However, most transcription factors contain nuclear localization signals that mediate the transport of transcription factors into the nucleus to regulate gene transcription [37]. This explains the observation in our study that the HOXA13 protein is mainly localized in the nucleus. A similar finding was reported by Cantile et al. [38], in which their study found that HOXA13 was predominantly expressed in the nucleus of thyroid cancer cells rather than the cytoplasm.
Some of the findings of our study are similar to those of Hu et al. [24], in which high HOXA13 protein expression in bladder cancer was significantly associated with positive lymph node metastasis. High HOXA13 expression was more prevalent in positive lymph node metastasis, whereas low HOXA13 expression was more frequently observed in negative lymph node metastasis. In addition, our study found a strong association between HOXA13 protein expression and tumor grade. A meta-analysis study on various human cancers (bladder, esophageal, gastric, liver, oral, ovarian, pancreatic, and prostate cancers) that was conducted by Wen et al. [39] also showed that high HOXA13 expression was significantly associated with positive lymph node metastasis. In addition, high HOXA13 expression was also significantly associated with higher tumor grade and TNM stage as well as poor overall survival.
Interestingly, our study revealed a pattern of high HOXA13 expression in 73–78% of pTa and pT1 tumors and 58–67% in pT3 and pT4 tumors, but only in 39% of pT2 tumors. Correspondingly, high HOXA13 expression was found in 67–85% of grade 1 and grade 3 tumors but only 39% of grade 2 tumors. Our study partially contradicts that of Hu et al. [24], who reported that high HOXA13 expression was associated with high stage (pT2–pT4) and high grade (grades 2 and 3) bladder cancer. The discrepancy between the findings of these studies may be due to population heterogeneity, different types of primary antibodies used, and variations in the IHC protocol [40].
Our study suggests that HOXA13 may be involved in different stages of bladder cancer progression. The roles of HOXA13 in the differentiation and morphogenesis of the genitourinary tracts [41] may explain the higher HOXA13 expression in the well-differentiated pTa and pT1 tumors. When bladder cancer cells progress to a more advanced stage, the HOXA13 protein is overexpressed and potentially plays a different role. It has been suggested that the HOXA13 protein may be involved in the morphogenesis of cancer cells by causing abnormalities in cell adhesion, which affect the interactions between cancer cells and their extracellular matrix [42], thus enabling tumor invasion and metastasis. Taken together, HOXA13 potentially plays different roles at the early and late stages of bladder cancer progression. These findings are significant, and further studies are needed to elucidate the precise stage-dependent role of HOXA13 in bladder cancer.
Findings from our study suggest that HOXA13 is a putative oncogene in bladder cancer and that it has the potential to be developed into a clinically applicable biomarker for the detection of bladder cancer. Thus, HOXA13 can potentially complement the use of cystoscopy, whereby cystoscopy can be used for confirmation upon the detection of high HOXA13 expression. After collecting voided urine samples from bladder cancer patients, RNA or protein can be extracted from the exfoliated tumor cells that are present in the urine. The expression of HOXA13 can be measured either at RNA or protein levels using RT-qPCR or ELISA, respectively. The use of HOXA13 as a potential clinically applicable biomarker will help to avoid the need to impose mandatory cystoscopies, which often inhibit patients from coming and getting treatment or screening for recurrence. Perhaps the likelihood of compliance with surveillance follow-up will be higher if the HOXA13 test is performed first instead of a cystoscopy.
Our study demonstrated that positive HOXB13 expression was found in 57.8% of the tissues analyzed, whereas negative HOXB13 expression was found in 42.2%, with the majority of the tissues being MIBC. In addition, low HOXB13 expression was significantly associated with MIBC, higher tumor stage (pT2–pT4) and grade (grades 2 and 3), and metastatic risk. These findings corroborate with those studies that were conducted on gastric and kidney cancers [33,43]. HOXB13 expression was significantly lower in gastric cancer patients with lymph node metastasis, higher stage, and poorly differentiated (grade 3) tumors [43]. In addition, HOXB13 methylation (loss of expression) in renal cell carcinoma was positively correlated with tumor grade and invasion [33].
Findings from our study showed that the HOXB13 protein was heterogeneously expressed in bladder cancer tissues and was mainly localized in the nucleus. Nuclear HOXB13 expression was significantly associated with NMIBC, allowing it to distinguish NMIBC from MIBC. Notably, the loss of HOXB13 protein expression was associated with progression to MIBC, suggesting HOXB13 may have a tumor suppressive role in bladder cancer. Incidentally, the putative role of HOXB13 as a tumor suppressor gene has been reported in other cancers, including colorectal [31], gastric [43], kidney [33], melanoma [44], and prostate cancer [45]. In contrast, Marra et al. [46] reported that HOXB13 is mainly localized in the cytoplasm of bladder cancer cells. Their study found that cytoplasmic HOXB13 protein expression was significantly higher in MIBC compared to NMIBC, indicating HOXB13 promotes bladder cancer invasiveness. Thus, further investigation is needed to clarify the contradictory role of HOXB13 in bladder cancer.
Our study demonstrated the association between both HOXA13 and HOXB13 expression with bladder tumorigenesis. However, only 110 patients fulfilled the inclusion criteria for this study, which is considered a small sample size for biomarker studies involving retrospectively collected specimens [47]. Small sample size studies may increase the chance of making inaccurate conclusions [48]. Hence, the findings of this study should be further validated by future studies with a larger sample size. This is crucial to thoroughly evaluating the potential of both HOXA13 and HOXB13 as biomarkers for bladder cancer.

5. Conclusions

Our study supports the oncogenic role of HOXA13 in bladder cancer, as this gene is involved in the early and later stages of bladder tumorigenesis, while HOXB13 may play a tumor suppressive role in bladder cancer. Thus, both the HOXA13 and HOXB13 genes have biomarker potential. Future work will focus on validating the clinical significance of HOXA13 and HOXB13 as biomarkers for bladder cancer. Moreover, it is necessary to conduct functional mechanistic candidate-gene approaches to clarify the emerging role of both genes in bladder cancer.

Author Contributions

Conceptualization, A.V., S.-C.C., K.Y., T.-A.O., R.Y. and A.H.A.R.; methodology, F.-W.C. and H.H.; formal analysis, F.-W.C., S.-C.C. and H.H.; investigation, F.-W.C.; writing—original draft, F.-W.C.; writing—review and editing, A.V., S.-C.C., H.H. and D.-M.C.; visualization, F.-W.C.; supervision, A.V., S.-C.C., D.-M.C. and H.H.; funding acquisition, A.V., S.-C.C., K.Y., T.-A.O., R.Y. and A.H.A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Higher Education, Malaysia, under the Fundamental Research Grant Scheme (FRGS), grant number FRGS/2/2014/SKK01/PERDANA/03/01; the Ministry of Energy, Science, Technology, Environmental and Climate Change, Malaysia, under the DSTIN Flagship Programme, grant number FP0514B0021-2 (DSTIN); and the Ministry of Higher Education, Malaysia, under the University of Malaya High Impact Research (HIR) Grant Scheme, grant number UM.C/625/1/HIR/130. F.-W.C. was sponsored by the MyBrain15 scholarship awarded by the Ministry of Higher Education, Malaysia.

Institutional Review Board Statement

The study was approved by the Medical Research and Ethics Committee, Ministry of Health Malaysia (approval number NMRR-12-970-13096).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors thank Juita Chupri and Zamzarina Ahmad Bajari from the Department of Pathology, Faculty of Medicine and Health Sciences, Universiti Putra Malaysia, for their assistance in the H&E staining and H&E image capturing, respectively.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  2. Kirkali, Z.; Chan, T.; Manoharan, M.; Algaba, F.; Busch, C.; Cheng, L.; Kiemeney, L.; Kriegmair, M.; Montironi, R.; Murphy, W.M.; et al. Bladder cancer: Epidemiology, staging and grading, and diagnosis. Urology 2005, 66 (Suppl. S1), 4–34. [Google Scholar] [CrossRef] [PubMed]
  3. Epstein, J.I.; Amin, M.B.; Reuter, V.R.; Mostofi, F.K. The World Health Organization/International Society of Urological Pathology consensus classification of urothelial (transitional cell) neoplasms of the urinary bladder. Bladder Consensus Conference Committee. Am. J. Surg. Pathol. 1998, 22, 1435–1448. [Google Scholar] [CrossRef] [PubMed]
  4. Fernandez-Gomez, J.; Madero, R.; Solsona, E.; Unda, M.; Martinez-Piñeiro, L.; Gonzalez, M.; Portillo, J.; Ojea, A.; Pertusa, C.; Rodriguez-Molina, J.; et al. Predicting nonmuscle invasive bladder cancer recurrence and progression in patients treated with bacillus Calmette-Guerin: The CUETO scoring model. J. Urol. 2009, 182, 2195–2203. [Google Scholar] [CrossRef] [PubMed]
  5. Sylvester, R.J.; van der Meijden, A.P.M.; Oosterlinck, W.; Witjes, J.A.; Bouffioux, C.; Denis, L.; Newling, D.W.W.; Kurth, K. Predicting recurrence and progression in individual patients with stage Ta T1 bladder cancer using EORTC risk tables: A combined analysis of 2596 patients from seven EORTC trials. Eur. Urol. 2006, 49, 466–477. [Google Scholar] [CrossRef]
  6. Sievert, K.D.; Amend, B.; Nagele, U.; Schilling, D.; Bedke, J.; Horstmann, M.; Hennenlotter, J.; Kruck, S.; Stenzl, A. Economic aspects of bladder cancer: What are the benefits and costs? World J. Urol. 2009, 27, 295–300. [Google Scholar] [CrossRef]
  7. Stenzl, A.; Hennenlotter, J.; Schilling, D. Can we still afford bladder cancer? Curr. Opin. Urol. 2008, 18, 488–492. [Google Scholar] [CrossRef]
  8. Liu, S.; Hou, J.; Zhang, H.; Wu, Y.; Hu, M.; Zhang, L.; Xu, J.; Na, R.; Jiang, H.; Ding, Q. The evaluation of the risk factors for non-muscle invasive bladder cancer (NMIBC) recurrence after transurethral resection (TURBt) in Chinese population. PLoS ONE 2015, 10, e0123617. [Google Scholar] [CrossRef]
  9. Zhu, C.Z.; Ting, H.N.; Ng, K.H.; Ong, T.A. A review on the accuracy of bladder cancer detection methods. J. Cancer 2019, 10, 4038–4044. [Google Scholar] [CrossRef] [PubMed]
  10. Chamie, K.; Saigal, C.S.; Lai, J.; Hanley, J.M.; Setodji, C.M.; Konety, B.R.; Litwin, M.S. Urologic Diseases in America Project Compliance with guidelines for patients with bladder cancer: Variation in the delivery of care. Cancer 2012, 117, 5392–5401. [Google Scholar] [CrossRef]
  11. Aldoghachi, A.F.; Baharudin, A.; Ahmad, U.; Chan, S.C.; Ong, T.A.; Yunus, R.; Razack, A.H.; Yusoff, K.; Veerakumarasivam, A. Evaluation of CERS2 gene as a potential biomarker for bladder cancer. Dis. Markers 2019, 2019, 3875147. [Google Scholar] [CrossRef]
  12. Gurung, P.M.S.; Veerakumarasivam, A.; Williamson, M.; Counsell, N.; Douglas, J.; Tan, W.S.; Feber, A.; Crabb, S.J.; Short, S.C.; Freeman, A.; et al. Loss of expression of the tumour suppressor gene AIMP3 predicts survival following radiotherapy in muscle-invasive bladder cancer. Int. J. Cancer 2015, 136, 709–720. [Google Scholar] [CrossRef]
  13. Veerakumarasivam, A.; Goldstein, L.D.; Saeb-Parsy, K.; Scott, H.E.; Warren, A.; Thorne, N.P.; Mills, I.G.; Venkitaraman, A.; Neal, D.E.; Kelly, J.D. AURKA overexpression accompanies dysregulation of DNA-damage response genes in invasive urothelial cell carcinoma. Cell Cycle 2008, 7, 3525–3533. [Google Scholar] [CrossRef] [PubMed]
  14. Gehring, W.J.; Hiromi, Y. Homeotic genes and the homeobox. Annu. Rev. Genet. 1986, 20, 147–173. [Google Scholar] [CrossRef] [PubMed]
  15. Krumlauf, R. Hox genes in vertebrate development. Cell 1994, 78, 191–201. [Google Scholar] [CrossRef] [PubMed]
  16. HOXA13 Gene. Available online: https://ghr.nlm.nih.gov/gene/HOXA13 (accessed on 4 July 2019).
  17. HOXB13 Gene. Available online: https://ghr.nlm.nih.gov/gene/HOXB13 (accessed on 4 July 2019).
  18. Guo, B.; Che, T.; Shi, B.; Guo, L.; Yin, Y.; Li, L.; Wang, J.; Yan, D.; Chen, Y. Screening and identification of specific markers for bladder transitional cell carcinoma from urine urothelial cells with suppressive subtractive hybridization and cDNA microarray. Can. Urol. Assoc. J. 2011, 5, E129–E137. [Google Scholar] [CrossRef]
  19. Guo, B.; Che, T.; Shi, B.; Guo, L.; Zhang, Z.; Li, L.; Cai, C.; Chen, Y. Interaction network analysis of differentially expressed genes and screening of cancer marker in the urine of patients with invasive bladder cancer. Int. J. Clin. Exp. Med. 2015, 8, 3619–3628. [Google Scholar]
  20. Holyoake, A.; O’Sullivan, P.; Pollock, R.; Best, T.; Watanabe, J.; Kajita, Y.; Matsui, Y.; Ito, M.; Nishiyama, H.; Kerr, N.; et al. Development of a multiplex RNA urine test for the detection and stratification of transitional cell carcinoma of the bladder. Clin. Cancer Res. 2008, 14, 742–749. [Google Scholar] [CrossRef]
  21. Duan, R.; Han, L.; Wang, Q.; Wei, J.; Chen, L.; Zhang, J.; Kang, C.; Wang, L. HOXA13 is a potential GBM diagnostic marker and promotes glioma invasion by activating the Wnt and TGF-β pathways. Oncotarget 2015, 6, 27778–27793. [Google Scholar] [CrossRef]
  22. Gu, Z.D.; Shen, L.Y.; Wang, H.; Chen, X.M.; Li, Y.; Ning, T.; Chen, K.N. HOXA13 promotes cancer cell growth and predicts poor survival of patients with esophageal squamous cell carcinoma. Cancer Res. 2009, 69, 4969–4973. [Google Scholar] [CrossRef]
  23. He, Y.X.; Song, X.H.; Zhao, Z.Y.; Zhao, H. HOXA13 upregulation in gastric cancer is associated with enhanced cancer cell invasion and epithelial-to-mesenchymal transition. Eur. Rev. Med. Pharmacol. Sci. 2017, 21, 258–265. [Google Scholar] [CrossRef]
  24. Hu, H.Y.; Chen, Y.B.; Cheng, S.; Li, G.H.; Zhang, Z.G. Dysregulated expression of homebox gene HOXA13 is correlated with the poor prognosis in bladder cancer. Wien. Klin. Wochenschr. 2017, 129, 391–397. [Google Scholar] [CrossRef]
  25. Quagliata, L.; Quintavalle, C.; Lanzafame, M.; Matter, M.S.; Novello, C.; di Tommaso, L.; Pressiani, T.; Rimassa, L.; Tornillo, L.; Roncalli, M.; et al. High expression of HOXA13 correlates with poorly differentiated hepatocellular carcinomas and modulates sorafenib response in in vitro models. Lab. Investig. 2018, 98, 95–105. [Google Scholar] [CrossRef]
  26. Dong, Y.; Cai, Y.; Liu, B.; Jiao, X.; Li, Z.; Guo, D.; Li, X.; Wang, Y.; Yang, D. HOXA13 is associated with unfavorable survival and acts as a novel oncogene in prostate carcinoma. Futur. Oncol. 2017, 13, 1505–1516. [Google Scholar] [CrossRef] [PubMed]
  27. Qu, L.; Zhong, Y.; Zheng, Z.; Zhao, R. CDH17 is a downstream effector of HOXA13 in modulating the Wnt/β-catenin signaling pathway in gastric cancer. Eur. Rev. Med. Pharmacol. Sci. 2017, 21, 1234–1241. [Google Scholar]
  28. Miao, J.Y.; Wang, Z.C.; Provencher, H.; Muir, B.; Dahiya, S.; Carney, E.; Leong, C.O.; Sgroi, D.C.; Orsulic, S. HOXB13 promotes ovarian cancer progression. Proc. Natl. Acad. Sci. USA 2007, 104, 17093–17098. [Google Scholar] [CrossRef]
  29. Zhai, L.L.; Wu, Y.; Cai, C.Y.; Tang, Z.G. Overexpression of homeobox B-13 correlates with angiogenesis, aberrant expression of EMT markers, aggressive characteristics and poor prognosis in pancreatic carcinoma. Int. J. Clin. Exp. Pathol. 2015, 8, 6919–6927. [Google Scholar]
  30. Zhu, J.Y.; Sun, Q.K.; Wang, W.; Jia, W.D. High-level expression of HOXB13 is closely associated with tumor angiogenesis and poor prognosis of hepatocellular carcinoma. Int. J. Clin. Exp. Pathol. 2014, 7, 2925–2933. [Google Scholar] [PubMed]
  31. Jung, C.; Kim, R.S.; Zhang, H.; Lee, S.J.; Sheng, H.; Loehrer, P.J.; Gardner, T.A.; Jeng, M.H.; Kao, C. HOXB13 is downregulated in colorectal cancer to confer TCF4-mediated transactivation. Br. J. Cancer 2005, 92, 2233–2239. [Google Scholar] [CrossRef]
  32. Jung, C.Y.; Kim, R.S.; Zhang, H.J.; Lee, S.J.; Jeng, M.H. HOXB13 induces growth suppression of prostate cancer cells as a repressor of hormone-activated androgen receptor signaling. Cancer Res. 2004, 64, 9185–9192. [Google Scholar] [CrossRef] [PubMed]
  33. Okuda, H.; Toyota, M.; Ishida, W.; Furihata, M.; Tsuchiya, M.; Kamada, M.; Tokino, T.; Shuin, T. Epigenetic inactivation of the candidate tumor suppressor gene HOXB13 in human renal cell carcinoma. Oncogene 2006, 25, 1733–1742. [Google Scholar] [CrossRef] [PubMed]
  34. Kristiansen, I.; Stephan, C.; Jung, K.; Dietel, M.; Rieger, A.; Tolkach, Y.; Kristiansen, G. Sensitivity of HOXB13 as a diagnostic immunohistochemical marker of prostatic origin in prostate cancer metastases: Comparison to PSA, prostein, androgen receptor, ERG, NKX3.1, PSAP, and PSMA. Int. J. Mol. Sci. 2017, 18, 1151. [Google Scholar] [CrossRef] [PubMed]
  35. Haller, K.; Rambaldi, I.; Daniels, E.; Featherstone, M. Subcellular localization of multiple PREP2 isoforms is regulated by actin, tubulin, and nuclear export. J. Biol. Chem. 2004, 279, 49384–49394. [Google Scholar] [CrossRef] [PubMed]
  36. Stevens, K.E.; Mann, R.S. A balance between two nuclear localization sequences and a nuclear export sequence governs extradenticle subcellular localization. Genetics 2007, 175, 1625–1636. [Google Scholar] [CrossRef] [PubMed]
  37. Liu, Y.; Li, P.; Fan, L.; Wu, M. The nuclear transportation routes of membrane-bound transcription factors. Cell Commun. Signal. 2018, 16, 12. [Google Scholar] [CrossRef]
  38. Cantile, M.; Scognamiglio, G.; La Sala, L.; La Mantia, E.; Scaramuzza, V.; Valentino, E.; Tatangelo, F.; Losito, S.; Pezzullo, L.; Chiofalo, M.G.; et al. Aberrant expression of posterior HOX genes in well differentiated histotypes of thyroid cancers. Int. J. Mol. Sci. 2013, 14, 21727–21740. [Google Scholar] [CrossRef]
  39. Wen, Y.; Shu, F.; Chen, Y.; Chen, Y.; Lan, Y.; Duan, X.; Zhao, S.C.; Zeng, G. The prognostic value of HOXA13 in solid tumors: A meta-analysis. Clin. Chim. Acta 2018, 483, 64–68. [Google Scholar] [CrossRef]
  40. Heidari, Z.; Mahmoudzadeh-Sagheb, H.; Jahantigh, M.; Gorgich, E.A.C. Immunohistochemical expression of KI67 and HER2 in colorectal cancer compared to adenomatous and normal samples. Int. J. Cancer Manag. 2017, 10, e12252. [Google Scholar] [CrossRef]
  41. Scott, V.; Morgan, E.A.; Stadler, H.S. Genitourinary functions of Hoxa13 and Hoxd13. J. Biochem. 2005, 137, 671–676. [Google Scholar] [CrossRef]
  42. Hirohashi, S.; Kanai, Y. Cell adhesion system and human cancer morphogenesis. Cancer Sci. 2003, 94, 575–581. [Google Scholar] [CrossRef]
  43. Sui, B.Q.; Zhang, C.D.; Liu, J.C.; Wang, L.; Dai, D.Q. HOXB13 expression and promoter methylation as a candidate biomarker in gastric cancer. Oncol. Lett. 2018, 15, 8833–8840. [Google Scholar] [CrossRef]
  44. Muthusamy, V.; Duraisamy, S.; Bradbury, C.M.; Hobbs, C.; Curley, D.P.; Nelson, B.; Bosenberg, M. Epigenetic silencing of novel tumor suppressors in malignant melanoma. Cancer Res. 2006, 66, 11187–11193. [Google Scholar] [CrossRef]
  45. Jung, C.Y.; Kim, R.S.; Lee, S.J.; Wang, C.H.; Jeng, M.H. HOXB13 homeodomain protein suppresses the growth of prostate cancer cells by the negative regulation of T-cell factor 4. Cancer Res. 2004, 64, 3046–3051. [Google Scholar] [CrossRef] [PubMed]
  46. Marra, L.; Cantile, M.; Scognamiglio, G.; Marra, L.; Perdona, S.; La Mantia, E.; Cerrone, M.; Gigantino, V.; Cillo, C.; Caraglia, M.; et al. Deregulation of HOX B13 expression in urinary bladder cancer progression. Curr. Med. Chem. 2013, 20, 833–839. [Google Scholar] [CrossRef]
  47. Popat, S.; Chen, Z.; Zhao, D.; Pan, H.; Hearle, N.; Chandler, I.; Shao, Y.; Aherne, W.; Houlston, R.S. A prospective, blinded analysis of thymidylate synthase and p53 expression as prognostic markers in the adjuvant treatment of colorectal cancer. Ann. Oncol. 2006, 17, 1810–1817. [Google Scholar] [CrossRef] [PubMed]
  48. Faber, J.; Fonseca, L.M. How sample size influences research outcomes. Dental Press J. Orthod. 2014, 19, 27–29. [Google Scholar] [CrossRef] [PubMed]
Diagnostics 13 02636 g001
Figure 1. HOXA13 protein localization in FFPE bladder cancer tissues. HOXA13 protein was predominantly localized in the nucleus. In a small proportion of tissues, HOXA13 protein was also localized in both the nucleus and cytoplasm.
Figure 1. HOXA13 protein localization in FFPE bladder cancer tissues. HOXA13 protein was predominantly localized in the nucleus. In a small proportion of tissues, HOXA13 protein was also localized in both the nucleus and cytoplasm.
Diagnostics 13 02636 g001
Diagnostics 13 02636 g002
Figure 2. Hematoxylin and eosin (H&E) and immunohistochemical (IHC) staining of HOXA13 in FFPE bladder cancer tissues with nuclear localization. (1A–6A) are H&E staining, where (1A,2A) are FFPE human skin tissues and (3A–6A) are FFPE bladder cancer tissues. (1B–6B) are IHC staining of HOXA13. FFPE human skin tissues were used as (1B) positive and (2B) negative controls. FFPE bladder cancer tissues showed (3B) no staining, (4B) weak staining, (5B) moderate staining, and (6B) strong staining. All images were taken at 400× magnification.
Figure 2. Hematoxylin and eosin (H&E) and immunohistochemical (IHC) staining of HOXA13 in FFPE bladder cancer tissues with nuclear localization. (1A–6A) are H&E staining, where (1A,2A) are FFPE human skin tissues and (3A–6A) are FFPE bladder cancer tissues. (1B–6B) are IHC staining of HOXA13. FFPE human skin tissues were used as (1B) positive and (2B) negative controls. FFPE bladder cancer tissues showed (3B) no staining, (4B) weak staining, (5B) moderate staining, and (6B) strong staining. All images were taken at 400× magnification.
Diagnostics 13 02636 g002
Diagnostics 13 02636 g003
Figure 3. HOXB13 protein localization in FFPE bladder cancer tissues. HOXB13 protein was predominantly localized in the nucleus. In a smaller proportion of tissues, HOXB13 protein was localized in either both the nucleus and cytoplasm, or the cytoplasm only.
Figure 3. HOXB13 protein localization in FFPE bladder cancer tissues. HOXB13 protein was predominantly localized in the nucleus. In a smaller proportion of tissues, HOXB13 protein was localized in either both the nucleus and cytoplasm, or the cytoplasm only.
Diagnostics 13 02636 g003
Diagnostics 13 02636 g004
Figure 4. Hematoxylin and eosin (H&E) and immunohistochemical (IHC) staining of HOXB13 in FFPE bladder cancer tissues with nuclear localization. (1A–6A) are H&E staining where (1A,2A) are FFPE human prostate cancer tissues and (3A–6A) are FFPE bladder cancer tissues. (1B–6B) are IHC staining of HOXB13. FFPE human prostate cancer tissues were used as (1B) positive and (2B) negative controls. FFPE bladder cancer tissues showed (3B) no staining, (4B) weak staining, (5B) moderate staining, and (6B) strong staining. All images were taken at 400× magnification.
Figure 4. Hematoxylin and eosin (H&E) and immunohistochemical (IHC) staining of HOXB13 in FFPE bladder cancer tissues with nuclear localization. (1A–6A) are H&E staining where (1A,2A) are FFPE human prostate cancer tissues and (3A–6A) are FFPE bladder cancer tissues. (1B–6B) are IHC staining of HOXB13. FFPE human prostate cancer tissues were used as (1B) positive and (2B) negative controls. FFPE bladder cancer tissues showed (3B) no staining, (4B) weak staining, (5B) moderate staining, and (6B) strong staining. All images were taken at 400× magnification.
Diagnostics 13 02636 g004
Diagnostics 13 02636 g005
Figure 5. Protein expression of HOXA13 in FFPE bladder cancer tissues. The distribution of HOXA13 protein expression was stratified according to different (A) tumor stages and (B) tumor grades. HOXA13 protein expression was categorized into low and high HOXA13 expression.
Figure 5. Protein expression of HOXA13 in FFPE bladder cancer tissues. The distribution of HOXA13 protein expression was stratified according to different (A) tumor stages and (B) tumor grades. HOXA13 protein expression was categorized into low and high HOXA13 expression.
Diagnostics 13 02636 g005
Diagnostics 13 02636 g006
Figure 6. Protein expression of HOXB13 in FFPE bladder cancer tissues. The distribution of HOXB13 protein expression was stratified according to different (A) tumor stages and (B) tumor grades. HOXB13 protein expression was categorized into low and high HOXB13 expression.
Figure 6. Protein expression of HOXB13 in FFPE bladder cancer tissues. The distribution of HOXB13 protein expression was stratified according to different (A) tumor stages and (B) tumor grades. HOXB13 protein expression was categorized into low and high HOXB13 expression.
Diagnostics 13 02636 g006
Table 1. IHC scoring system for HOXA13 protein expression.
Table 1. IHC scoring system for HOXA13 protein expression.
Percentage of Positive CellsScoreStaining IntensityScore
Less than 1% positive tumor cells0No staining0
1 to 30% positive tumor cells1Weak staining1
30 to 70% positive tumor cells2Moderate staining2
More than 70% positive tumor cells3Strong staining3
Final IHC score = score of percentage of positive cells × score of staining intensity; final IHC score 0, 1, 2, and 3 = low HOXA13 expression; final IHC score 4, 6, and 9 = high HOXA13 expression.
Table 2. IHC scoring system for HOXB13 protein expression.
Table 2. IHC scoring system for HOXB13 protein expression.
Percentage of Positive CellsScoreStaining IntensityScore
0% positive tumor cells0No staining0
1 to 10% positive tumor cells1Weak staining1
11 to 50% positive tumor cells2Moderate staining2
50 to 80% positive tumor cells3Strong staining3
More than 80% positive tumor cells4
Final IHC score = score of percentage of positive cells × score of staining intensity; final IHC score 0, 1, 2, 3, and 4 = low HOXB13 expression; final IHC score 6, 8, 9, and 12 = high HOXB13 expression.
Table 3. Association between HOXA13 protein expression and the demographic as well as clinicopathological characteristics of bladder cancer patients.
Table 3. Association between HOXA13 protein expression and the demographic as well as clinicopathological characteristics of bladder cancer patients.
CharacteristicNumber of Cases (%)HOXA13 Expressionχ2
Value		p
Value		Cramer’s
V Value		Strength of Association
Low 2 (%)High (%)
Gender 1.2650.261--
Male105(95.5)37(35.2)68(64.8)
Female5(4.5)3(60.0)2(40.0)
Ethnicity 3.1320.372--
Malay65(59.1)25(38.5)40(61.5)
Chinese28(25.5)7(25.0)21(75.0)
Indian14(12.7)6(42.9)8(57.1)
Others3(2.7)2(66.7)1(33.3)
Age 1.9470.163--
<5013(11.8)7(53.8)6(46.2)
≥5097(88.2)33(34.0)64(66.0)
Tumor invasiveness 5.3840.020 *0.221Moderate
NMIBC49(44.5)12(24.5)37(75.5)
MIBC61(55.5)28(45.9)33(54.1)
Tumor stage 9.3500.053--
pTa23(20.9)5(21.7)18(78.3)
pT126(23.6)7(26.9)19(73.1)
pT223(20.9)14(60.9)9(39.1)
pT321(19.1)7(33.3)14(66.7)
pT417(15.5)7(41.2)10(58.8)
Tumor grade 1 11.4140.003 *0.322Strong
G120(18.2)3(15.0)17(85.0)
G226(23.6)16(61.5)10(38.5)
G364(58.2)21(32.8)43(67.2)
Lymph node metastasis 4.9010.027 *0.288Moderate
Yes26(23.6)9(34.6)17(65.4)
No33(30.0)21(63.6)12(36.4)
N/A51(46.4)----
Metastasis 1.8920.169--
Yes30(27.3)14(46.7)16(53.3)
No80(72.7)26(32.5)54(67.5)
Recurrence 4.0920.043 *0.193Weak
Yes66(60.0)19(28.8)47(71.2)
No44(40.0)21(47.7)23(52.3)
Abbreviations: MIBC = muscle invasive bladder cancer; NMIBC = non-muscle invasive bladder cancer; 1 WHO 1973 classification system; 2 low HOXA13 expression (including those with negative expression); * statistically significant at p < 0.05.
Table 4. Association between HOXA13 protein localization and the demographic as well as clinicopathological characteristics of bladder cancer patients.
Table 4. Association between HOXA13 protein localization and the demographic as well as clinicopathological characteristics of bladder cancer patients.
CharacteristicNumber of Cases (%)HOXA13 Localizationχ2
Value		p-ValueCramer’s
V Value		Strength of Association
NoneNucleusNucleus and Cytoplasm
Gender 0.5060.777--
Male105(95.5)6(5.7)86(81.9)13(12.4)
Female5(4.5)0(0.0)4(80.0)1(20.0)
Ethnicity 11.7280.068--
Malay65(59.1)3(4.6)53(81.5)9(13.8)
Chinese28(25.5)0(0.0)24(85.7)4(14.3)
Indian14(12.7)2(14.3)12(85.7)0(0.0)
Others3(2.7)1(33.3)1(33.3)1(33.3)
Age 15.1560.001 *0.371Strong
<5013(11.8)0(0.0)7(53.8)6(46.2)
≥5097(88.2)6(6.2)83(85.6)8(8.2)
Tumor invasiveness 2.0680.356--
NMIBC49(44.5)1(2.0)42(85.7)6(12.2)
MIBC61(55.5)5(8.2)48(78.7)8(13.1)
Tumor stage 11.2570.188--
pTa23(20.9)1(4.3)20(87.0)2(8.7)
pT126(23.6)0(0.0)22(84.6)4(15.4)
pT223(20.9)4(17.4)17(73.9)2(8.7)
pT321(19.1)1(4.8)18(85.7)2(9.5)
pT417(15.5)0(0.0)13(76.5)4(23.5)
Tumor grade 1 3.4380.487--
G120(18.2)1(5.0)15(75.0)4(20.0)
G226(23.6)1(3.8)20(76.9)5(19.2)
G364(58.2)4(6.2)55(85.9)5(7.8)
Lymph node metastasis 1.3360.513--
Yes26(23.6)1(3.8)21(80.8)4(15.4)
No33(30.0)4(12.1)25(75.8)4(12.1)
N/A51(46.4)------
Metastasis 0.1400.933--
Yes30(27.3)2(6.7)24(80.0)4(13.3)
No80(72.7)4(5.0)66(82.5)10(12.5)
Recurrence 7.5660.023 *0.262Moderate
Yes66(60.0)1(1.5)59(89.4)6(9.1)
No44(40.0)5(11.4)31(70.5)8(18.2)
Abbreviations: MIBC = muscle invasive bladder cancer; NMIBC = non-muscle invasive bladder cancer; 1 WHO 1973 classification system; * statistically significant at p < 0.05.
Table 5. Association between HOXB13 protein expression and the demographic as well as clinicopathological characteristics of bladder cancer patients.
Table 5. Association between HOXB13 protein expression and the demographic as well as clinicopathological characteristics of bladder cancer patients.
CharacteristicNumber of Cases (%)HOXB13 Expressionχ2
Value		p-ValueCramer’s
V Value		Strength of Association
Low 2 (%)High (%)
Gender 0.3460.556--
Male97(95.1)70(72.2)27(27.8)
Female5(4.9)3(60.0)2(40.0)
Ethnicity 3.0000.392--
Malay58(56.9)43(74.1)15(25.9)
Chinese27(26.5)18(66.7)9(33.3)
Indian14(13.7)11(78.6)3(21.4)
Others3(2.9)1(33.3)2(66.7)
Age 0.0400.841--
<5013(12.7)9(69.2)4(30.8)
≥5089(87.3)64(71.9)25(28.1)
Tumor invasiveness 23.212≤0.001 *0.477Strong
NMIBC46(45.1)22(47.8)24(52.2)
MIBC56(54.9)51(91.1)5(8.9)
Tumor stage 24.010≤0.0010.485Strong
pTa21(20.6)9(42.9)12(57.1)
pT125(24.5)13(52.0)12(48.0)
pT223(22.5)20(87.0)3(13.0)
pT318(17.6)17(94.4)1(5.6)
pT415(14.7)14(93.3)1(6.7)
Tumor grade 1 13.4270.001 *0.363Strong
G118(17.6)7(38.9)11(61.1)
G225(24.5)17(68.0)8(32.0)
G359(57.8)49(83.1)10(16.9)
Lymph node metastasis 0.7890.375--
Yes24(23.5)24(100.0)0(0.0)
No31(30.4)30(96.8)1(3.2)
N/A47(46.1)----
Metastasis 7.9760.005 *0.280Moderate
Yes27(26.5)25(92.6)2(7.4)
No75(73.5)48(64.0)27(36.0)
Recurrence 3.4100.065--
Yes63(61.8)41(65.1)22(34.9)
No39(38.2)32(82.1)7(17.9)
Abbreviations: MIBC = muscle invasive bladder cancer; NMIBC = non-muscle invasive bladder cancer; 1 WHO 1973 classification system; 2 low HOXB13 expression (including those with negative expression); * statistically significant at p < 0.05.
Table 6. Association between HOXB13 protein localization and the demographic as well as clinicopathological characteristics of bladder cancer patients.
Table 6. Association between HOXB13 protein localization and the demographic as well as clinicopathological characteristics of bladder cancer patients.
CharacteristicNumber of Cases (%)HOXB13 Localizationχ2
Value		p-ValueCramer’s
V Value		Strength of Association
NoneNucleusNucleus and CytoplasmCytoplasm
Gender 1.4370.697--
Male97(95.1)41(42.3)34(35.1)19(19.6)3(3.1)
Female5(4.9)2(40.0)1(20.0)2(40.0)0(0.0)
Ethnicity 2.2160.988--
Malay58(56.9)25(43.1)18(31.0)13(22.4)2(3.4)
Chinese27(26.5)10(37.0)11(40.7)5(18.5)1(3.7)
Indian14(13.7)7(50.0)5(35.7)2(14.3)0(0.0)
Others3(2.9)1(33.3)1(33.3)1(33.3)0(0.0)
Age 3.3460.341--
<5013(12.7)5(38.5)3(23.1)5(38.5)0(0.0)
≥5089(87.3)38(42.7)32(36.0)16(18.0)3(3.4)
Tumor invasiveness 26.927≤0.001 *0.514Strong
NMIBC46(45.1)8(17.4)27(58.7)10(21.7)1(2.2)
MIBC56(54.9)35(62.5)8(14.3)11(19.6)2(3.6)
Tumor stage 43.160≤0.001 *0.376Strong
pTa21(20.6)2(9.5)16(76.2)2(9.5)1(4.8)
pT125(24.5)6(24.0)11(44.0)8(32.0)0(0.0)
pT223(22.5)13(56.5)3(13.0)7(30.4)0(0.0)
pT318(17.6)12(66.7)3(16.7)1(5.6)2(11.1)
pT415(14.7)10(66.7)2(13.3)3(20.0)0(0.0)
Tumor grade 1 13.7690.032 *0.260Moderate
G118(17.6)1(5.6)10(55.6)6(33.3)1(5.6)
G225(24.5)13(52.0)9(36.0)3(12.0)0(0.0)
G359(57.8)29(49.2)16(27.1)12(20.3)2(3.4)
Lymph node metastasis 0.1540.985--
Yes24(23.5)17(70.8)3(12.5)3(12.5)1(4.2)
No31(30.4)23(74.2)4(12.9)3(9.7)1(3.2)
N/A47(46.1)--------
Metastasis 16.3090.001 *0.400Strong
Yes27(26.5)20(74.1)3(11.1)3(11.1)1(3.7)
No75(73.5)23(30.7)32(42.7)18(24.0)2(2.7)
Recurrence 5.6110.132--
Yes63(61.8)21(33.3)24(38.1)16(25.4)2(3.2)
No39(38.2)22(56.4)11(28.2)5(12.8)1(2.6)
Abbreviations: MIBC = muscle invasive bladder cancer; NMIBC = non-muscle invasive bladder cancer; 1 WHO 1973 classification system; * statistically significant at p < 0.05.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chin, F.-W.; Hussin, H.; Chau, D.-M.; Ong, T.-A.; Yunus, R.; Abdul Razack, A.H.; Yusoff, K.; Chan, S.-C.; Veerakumarasivam, A. Differential Protein Expression Patterns of HOXA13 and HOXB13 Are Associated with Bladder Cancer Progression. Diagnostics 2023, 13, 2636. https://doi.org/10.3390/diagnostics13162636

AMA Style

Chin F-W, Hussin H, Chau D-M, Ong T-A, Yunus R, Abdul Razack AH, Yusoff K, Chan S-C, Veerakumarasivam A. Differential Protein Expression Patterns of HOXA13 and HOXB13 Are Associated with Bladder Cancer Progression. Diagnostics. 2023; 13(16):2636. https://doi.org/10.3390/diagnostics13162636

Chicago/Turabian Style

Chin, Fee-Wai, Huzlinda Hussin, De-Ming Chau, Teng-Aik Ong, Rosna Yunus, Azad Hassan Abdul Razack, Khatijah Yusoff, Soon-Choy Chan, and Abhi Veerakumarasivam. 2023. "Differential Protein Expression Patterns of HOXA13 and HOXB13 Are Associated with Bladder Cancer Progression" Diagnostics 13, no. 16: 2636. https://doi.org/10.3390/diagnostics13162636

APA Style

Chin, F. -W., Hussin, H., Chau, D. -M., Ong, T. -A., Yunus, R., Abdul Razack, A. H., Yusoff, K., Chan, S. -C., & Veerakumarasivam, A. (2023). Differential Protein Expression Patterns of HOXA13 and HOXB13 Are Associated with Bladder Cancer Progression. Diagnostics, 13(16), 2636. https://doi.org/10.3390/diagnostics13162636

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

Article Metrics

Back to TopTop