HECTD1 mRNA was downregulated in breast cancer tissues
We first investigated if HECTD1 mRNA expression levels are altered in breast cancer tissues by analyzing 114 matched tumor and adjacent normal tissues in TCGA cohort. We found that HECTD1 mRNA expression levels were significantly lower in the breast tumor tissues compared with the adjacent normal tissues (P<0.001) (Fig. 1).
HECTD1 expression differed among breast cancer subtypes
We examined if HECTD1 mRNA expression levels are different among breast cancer subtypes in the NCU cohort. HECTD1 mRNA expression levels were significantly lower in ER-negative tumors (P<0.001), PgR-negative tumors (P<0.001), and HER2-positive tumors (P<0.001) compared with the respective counterparts (Fig. 2). We further stratified the patients into four subtypes based on IHC findings: luminal (ER-positive and/or PgR-positive, HER2-negative), luminal-HER2 (ER-positive and/or PgR-positive, HER2-positive), HER2 (ER-negative, PgR-negative, HER2-positive), and triple-negative (ER-negative, PgR-negative, HER2-negative). HECTD1 mRNA expression levels were significantly lower in the HER2 and triple-negative subtypes compared with the luminal subtype (P=0.001 and P<0.001, respectively) (Fig. 2).
HECTD1 mRNA level was negatively associated with aggressive tumor characteristics
Based on the downregulation of HECTD1 mRNA in tumor tissues (Fig. 1) and the previous reports showing that HECTD1 suppressed EMT10,11, we investigated the association between HECTD1 mRNA expression and tumor aggressiveness in the NCU cohort. HECTD1 mRNA expression was significantly lower in larger tumors (P=0.009) and higher-grade tumors (grade 3 vs. grade 1, P<0.001; grade 3 vs. grade 2, P<0.001; grade 2 vs. grade 1, P=0.042) (Fig. 3). However, HECTD1 mRNA expression was not associated with lymph node metastasis.
We also investigated the associations between HECTD1 mRNA expression levels and other patient demographics in the NCU cohort. High HECTD1 mRNA expression was associated with a higher proportion of invasive lobular carcinomas compared with tumors with low HECTD1 expression (P=0.009) (Supplementary Table S1). However, because of the small number of lobular carcinomas, the results could not be determined to be statistically significant. There was no significant association between HECTD1 mRNA expression and other examined factors, including age, sex, or menopausal status.
HECTD1 was inversely associated with genes involved in mitochondrial cellular respiratory function
HECTD1 was shown to negatively regulate EMT (Duhamel et al. 2018). We therefore investigated the relationship between HECTD1 mRNA expression and EMT using TCGA cohorts. In contrast to the previous report (Duhamel et al. 2018), GSEA revealed no correlation between HECTD1 mRNA expression and the EMT gene set (Fig. 4a). In contrast, we found that the mitochondrial gene set, mitochondrion (P=0.009), and especially the respiratory-related gene sets for OXPHOS (P=0.011), the respiratory chain complex (P=0.006), and the ROS pathway (P=0.009) were inversely correlated with HECTD1 mRNA expression in TCGA cohort (Fig. 4b). These findings thus indicated that HECTD1 might be directly or indirectly involved in the negative regulation of mitochondrial cellular respiratory function in breast cancer.
Low HECTD1 mRNA expression was associated with poor prognosis in breast cancer
We analyzed the impact of HECTD1 mRNA expression on the survival of breast cancer patients. Patients with HECTD1 mRNA low-expressing tumors had significantly shorter DFS and OS among all breast cancer patients (P=0.013 and P=0.038, respectively) in the NCU cohort (Fig. 5a). Univariate analyses showed that larger tumor size (DFS, P=0.032; OS, P=0.002), lymph node metastasis (DFS, P<0.001; OS, P<0.001), higher histological grade (DFS, P=0.003; OS, P=0.002), ER-negativity (DFS, P=0.004; OS, P=0.005), HER2-positivity (DFS, P=0.015; OS, P=0.011), and low HECTD1 mRNA expression level (DFS, P=0.014; OS, P=0.039) were all significantly associated with shorter DFS and OS (Table 1). Multivariate analyses identified low HECTD1 mRNA expression level as an independent factor associated with lower DFS (hazard ratio (HR): 1.54, 95% confidence interval (CI): 1.11–2.13, P=0.009), together with lymph node metastasis (P<0.001) and ER status (P=0.009) (Table 1). Low HECTD1 mRNA expression was also an independent prognostic factor for poor OS (HR: 1.50, 95% CI: 1.01–2.24, P=0.046), together with lymph node metastasis (P<0.001) and ER status (P=0.037) (Table 1).
Given that HECTD1 mRNA expression levels differed according to subtype, we further investigated the impact of HECTD1 mRNA expression on patient survival in each subtype. There was no significant difference in either DFS or OS between the HECTD1 high- and low-expressing tumors in ER-positive (DFS, P=0.273; OS, P=0.388) (Fig. 5b) and HER2-positive patients (DFS, P=0.617; OS, P=0.972) (Supplementary Fig. S1 online). However, low HECTD1 expression tended to be associated with worse prognosis in ER-negative (DFS, P=0.079; OS, P=0.139) (Fig. 5c) and triple-negative patients (DFS, P=0.123; OS, P=0.106) (Supplementary Fig. S1 online).
Lack of association of HECTD1 protein expression with mRNA expression or patient survival in breast cancer
We investigated if HECTD1 mRNA expression was correlated with its protein expression levels in breast cancer. Representative images of HECTD1 protein expression detected by IHC are shown in Fig. 6a. There was no significant correlation between HECTD1 mRNA expression levels and HECTD1 H-score in the NCU cohort (r=−0.089, P=0.134) (Fig. 6b) or protein expression by mass spectrometry in TCGA cohort (r=0.205, P=0.079) (Fig. 6c). In addition, there was no significant difference in DFS or OS between patients with HECTD1 high- and low-expressing tumors in the NCU cohort (P=0.698 and P=0.819, respectively) (Fig. 6d).