Upregulated Ubiquitin D is a Favorable Prognostic Indicator for Rectal Cancer Patients Undergoing Preoperative Concurrent Chemoradiotherapy.

Purpose: For locally advanced rectal cancer, neoadjuvant concurrent chemoradiotherapy (CCRT) allows tumor downstaging and makes curative radical proctectomy possible. However, we lack a genetic biomarker to predict cancer prognosis or treatment response. We investigated the association between ubiquitin D (UBD) expression and clinical outcomes in rectal cancer patients receiving CCRT. Patients and
Methods: We analyzed the genes associated with the protein modification process (GO:0036211) and identified the UBD gene as the most relevant among the top 7 differentially expressed genes associated with CCRT resistance. We collected tissue specimens from 172 rectal cancer patients who had received CCRT followed by a curative proctectomy. We examine the relationship between UBD expression and patient characteristics, pathological findings, and patient survival, such as metastasis-free survival (MeFS) and disease-specific survival.
Results: Upregulated UBD expression was associated with lower pre-CCRT tumor T stage (P = 0.009), lower post-CCRT tumor T stage (P < 0.001), lower post-CCRT nodal stage (P < 0.001), less vascular invasion (P = 0.015), and better tumor regression (P < 0.001). Using univariate analysis, we found that high UBD expression was correlated with better disease-free survival (DFS) (P < 0.0001), local recurrence-free survival (LRFS) (P < 0.0001) and MeFS (P < 0.0001). Moreover, multivariate analysis demonstrated that high UBD expression was associated with superior DFS (P < 0.001), LRFS (P = 0.01), and MeFS (P = 0.004).
Conclusion: UBD upregulation was linked to better clinical prognosis, favorable pathological features, and good treatment response in rectal cancer patients undergoing CCRT. These results suggest UBD is a biomarker for rectal cancer.

Introduction

Colorectal cancer (CRC) is a common disease in the modern world. Although CRC mortality has been progressively declining since 1990,1 it still remains the third most common cause of cancer-related death in the United States in women and the second leading cause of death in men. Rectal cancer accounting for one-third of these cases. In contrast to these trends, the incidence of CRC in men and women under the age of 50 climbed at a rate of 2.1% per year from 1992 to 2012, and has continued to rise since.2 These increases are driven predominantly by left-sided cancers in general and rectal cancer in particular (3.9% per year).3 The Taiwan Cancer Registry shows that male young-onset rectal cancer incidence rates increased from 4.0 per 100,000 in 1995 to 8.3 per 100,000 in 2014. Female young-onset rectal cancer incidence also increased significantly from 3.8 per 100,000 to 6.4 per 100,000 over the same time period.4

Although surgical resection is the cornerstone of curative therapy for patients with potentially resectable rectal cancer, radiotherapy with concurrent chemotherapy (termed concurrent chemoradiotherapy [CCRT]) has emerged as an important component of curative therapy for transmural or node-positive rectal cancers because local recurrences are more common than with colon primaries. Neoadjuvant CCRT before surgery allows tumor downstaging and increases the rate of sphincter conservation. However, only approximately 20% of the patients can achieve pathological complete response (pCR).5 We had known that patients who achieve pCR have better outcomes and less recurrence of their cancer. Moreover, considering the declining age group of patients with rectal cancer, maintaining the long-term quality of life of rectal cancer patients after tumor resection has become very important and urgent. A biomarker predicting pCR could be useful in treatment decisions.

Ubiquitination is a kind of posttranslational modification involved in the regulation of many biological processes, including the cell cycle, differentiation, transcription regulation, signal transmission, and damage repair.6 Ubiquitin D (UBD), also known as Human HLA-F adjacent transcript locus 10 (FAT10), is a protein containing 165 amino acid residues. It belongs to a class of ubiquitin-like proteins and is a member of the ubiquitin-like protein family, whose protein sequence and three-dimensional core structure resembles ubiquitin (UB).7 Binding of the UBD molecule to its target protein will result in the protein degrading due to the proteasome.8 Deregulation of UBD may result in aberrant modifications in apoptosis, cell division, or chromosomal instability, all of which are linked to cancer.9 Increased studies on UBD in malignancies have emerged in recent years, revealing a link between increased UBD expression and disease development in several tumors,10,11 such as hepatocellular carcinoma,12,13 colon cancer,14 gastric cancer,15 oral squamous cell carcinoma,16 breast cancer,17 and glioma.18 However, the role of UBD in rectal cancer had never been explored. This study aimed to investigate the associations between UBD expression and clinical outcomes in nonmetastatic rectal cancer patients receiving neoadjuvant CCRT.

Materials and Methods

Data Mining of a Public Transcriptome Dataset

A public rectal cancer dataset (GSE35452) with 46 patients receiving CCRT followed by curative resection was used for transcriptome profiling to determine the efficacy of preoperative CCRT. Biopsy specimens were taken during a colonoscopy prior to CCRT in this dataset. The raw CEL files from the Affymetrix Human Genome U133 Plus 2.0 microarray platform were loaded into the Nexus Expression 3 statistical program (BioDiscovery, El Segundo, CA, USA) to evaluate all probes without filtering in order to computerize expression levels. The tumor specimens were divided into two groups by tumor regression grade (TRG) based on their response to neoadjuvant CCRT: “responders” and “non-responders”. The statistical significance of each transcript was tested by comparing responders to non-responders. We selected genes with a p-value less than 0.0001 and absolute expression fold change >1 log2 ratio for further analysis among differentially expressed genes associated with the protein modification process (GO:0036211). From the same database, we profiled the expression level of all genes and probes and identified 200 top-ranking genes whose expression is positively and negatively correlated with UBD. These genes of selection were also submitted to perform enrichment ontology analysis on the GENEONTOLOGY website () and assess their molecular function.

Patient Eligibility and Enrollment

This study was approved by the Institutional Review Board of Chi Mei Medical Center (10302014). All patients had signed the informed consent. Our study complies with the Declaration of Helsinki. We collected formalin-fixed, paraffin-embedded Pre-Tx rectal tissue specimens from 172 rectal cancer patients. Imaging tests were used to determine the initial clinical stage, and those who had been diagnosed with distant metastases were ruled out. All patients received 5-fluorouracil-based chemotherapy as well as radiation (45–50 Gy) in 25 fractions over a 5-week period, followed by a curative proctectomy. Those with a nodal status higher than N1 or a pre- or post-CCRT tumor status greater than T3 received adjuvant chemotherapy. After diagnosis, all patients were continuously observed until death or the last follow-up.

Histopathological and Immunohistochemical Evaluations

Two independent pathologists, blinded to the patients’ clinical information, reviewed all tumor specimens to achieve more objective results. The seventh American Joint Committee on Cancer TNM staging system was used to identify the T and N stages. The tumor regression grade (TRG), which is predictive of the tumor response to CCRT, was evaluated in concordance with the description by Dworak et al.19 Immunohistochemical study was performed to assess the expression of UBD. In brief, tissue sections from Pre-Tx rectal tumor biopsies were cut from paraffin-embedded tissue blocks at 3 mm thickness onto precoated glass slides. Slides were then deparaffinized with xylene, rehydrated with ethanol and heated for 7 min by microwave for antigen retrieval in a 10-mM citrate buffer (pH 6). Endogenous peroxidase was blocked by using 3% H2O2. Slides were then washed with Tris-buffered saline for 15 min and then incubated with a primary poly-clonal antibody against UBD (Cat No. NBP2-81752, Novus Biologicals, dilution 1:100). Antibodies were then detected using a DAKO ChemMate EnVision Kit (K5001, Carpinteria, CA, USA). Cell blocks from cell lines known to express UBD were used as positive controls. Sections were processed without the primary anti-UBD antibody as negative controls. The H-score was applied to evaluate UBD immunoreactivity and was quantified with the following equation:

H-score = SPi (i + 1)

where Pi is the percentage of stained tumor cells for each intensity, which ranges from 0% to 100%, and i is the staining intensity (0 to 3+). The H-score, which ranged from 100 to 400, was calculated using a combination of the intensity and proportion of positively stained tumor cells. High UBD expression was defined as H-scores more than or equal to the median of all scored instances.

Statistical Analysis

The associations between clinicopathological characteristics and UBD expression were measured using the chi-square (X2) test. The Kaplan–Meier method was used to plot survival curves, and the Log rank test was employed to quantify and compare the time from surgery to death (or last seen alive) or recurrence (or last seen relapse-free). The Cox proportional hazards model for multivariate analysis includes those parameters clinically significant in the univariate analysis. SPSS software version 22.0 (IBM Corporation, Armonk, NY, USA) was used for all statistical analyses, and two-tailed tests with p < 0.05 were considered statistically significant.

Results

UBD is Recognized as the Most Significant Differentially Expressed Gene Connected with the Protein Modification Process

A published rectal cancer transcriptome dataset (GSE35452) was used for data mining to investigate the potential biomarkers of rectal cancer cells sensitive to preoperative CCRT. The dataset included 46 patients who had neoadjuvant CCRT followed by standardized curative resection. Twenty-four patients (52.2%) were classed as responders, whereas 22 patients (47.8%) were classified as non-responders, based on their response to CCRT. We found seven transcripts focusing on the protein modification process (GO:0036211) (Table 1 and Figure 1). UBD was chosen for further investigation because its expression was significantly lower in CCRT non-responders (p < 0.0001). This discovery prompted us to look into the clinical significance of UBD expression in rectal cancer.

Table 1

Summary of Differentially Expressed Genes Associated with Protein Modification Process (GO:0036211) in Relation to Response to CCRT in Rectal Carcinoma

Probe	Comparison Log Ratio	Comparison p-value	Gene Symbol	Gene Name	Biological Process	Molecular Function
202954_at	−0.4053	0.0085	UBE2C	Ubiquitin-conjugating enzyme E2C	Cell cycle, cell division, cyclin catabolic process, mitosis, phosphoinositide-mediated signaling, positive regulation of exit from mitosis, protein modification process, protein ubiquitination, spindle organization and biogenesis, ubiquitin cycle, ubiquitin-dependent protein catabolic process	Ligase activity, small conjugating protein ligase activity, ubiquitin-protein ligase activity
205890_s_at	−1.0343	<0.0001	UBD	Ubiquitin D	G-protein-coupled receptor protein signaling pathway, protein modification process, proteolysis, signal transduction, ubiquitin cycle	G-protein-coupled receptor activity, GABA-B receptor activity, metabotropic glutamate; GABA-B-like receptor activity, receptor activity, signal transducer activity
212445_s_at	0.6247	0.0004	NEDD4L	Neural precursor cell expressed; developmentally down-regulated 4-like	Cellular sodium ion homeostasis, excretion, positive regulation of endocytosis, protein modification process, protein ubiquitination, regulation of protein catabolic process, response to metal ion, sodium ion transport, ubiquitin cycle, water homeostasis	Ligase activity, protein binding, sodium channel regulator activity, ubiquitin-protein ligase activity
214971_s_at	−0.5482	0.0088	ST6GAL1	ST6 beta-galactoside alpha-2;6-sialyltranferase 1	Humoral immune response, oligosaccharide metabolic process, protein amino acid glycosylation, protein modification process	Beta-galactoside alpha-2,6-sialyltransferase activity, sialyltransferase activity, transferase activity, transferase activity; transferring glycosyl groups
222446_s_at	0.4044	0.0064	BACE2	Beta-site APP-cleaving enzyme 2	Membrane protein ectodomain proteolysis, negative regulation of amyloid precursor protein biosynthetic process, peptide hormone processing, protein modification process, protein secretion, proteolysis	Aspartic-type endopeptidase activity, aspartic-type signal peptidase activity, hydrolase activity, pepsin A activity, peptidase activity
227914_s_at	0.1259	0.0068	UBE3B	Ubiquitin protein ligase E3B	Protein modification process, ubiquitin cycle	Ligase activity, ubiquitin-protein ligase activity
236294_at	−0.2411	0.0004	HUWE1	HECT; UBA and WWE domain containing 1	Histone ubiquitination, mRNA transport, protein modification process, protein polyubiquitination, ubiquitin cycle	DNA binding, ligase activity, protein binding, ubiquitin-protein ligase activity

Figure 1

Expression profiling of genes related to the protein modification process (GO:0036211) and their relationship to CCRT response. We identified UBD as the most significantly upregulated gene among CCRT responders.

Clinicopathological Features of Patients with Rectal Carcinoma in Our Cohort

We obtained medical records from a total of 172 rectal cancer patients (Table 2). The median age was 63. Most patients were male (n = 108, 62.8%) and less than 70 years old (n = 106, 61.6%). Fifty-three percent of the patients (n = 91) had T3-T4 lesions; 27% of the patients (n = 47) had lymph node metastases before CCRT. After CCRT, 50% of the patients (n = 86) had T3-T4 lesions; 28% of the patients (n = 49) had lymph node metastases. The surgical specimens revealed vascular invasion in 8% of the patients (n =15) and perineural invasion in 3% of the patients (n=5). Tumor regression grade, which was used to evaluate treatment response to CCRT, varied from grade 0 to 1 (n = 37, 22%), grade 2 to 3 (n = 118, 69%) and grade 4 (n = 17, 10%).

Table 2

Associations and Comparisons Between UBD Expression and Clinicopathological Factors in 172 Rectal Cancer Patients Receiving Neoadjuvant CCRT

Parameter		No.	UBD Expression		p-value
			High Exp.	Low Exp.
Gender	Male	108	51	57	0.999
	Female	64	35	29
Age	<70	106	51	55	0.117
	≥70	66	35	31
Pre-Tx tumor status (Pre-T)	T1-T2	81	49	32	0.009*
	T3-T4	91	37	54
Pre-Tx nodal status (Pre-N)	N0	125	67	58	0.124
	N1-N2	47	19	28
Post-Tx tumor status (Post-T)	T1-T2	86	59	27	<0.001*
	T3-T4	86	27	59
Post-Tx nodal status (Post-N)	N0	123	74	49	<0.001*
	N1-N2	49	12	37
Vascular invasion	Absent	157	83	74	0.015*
	Present	15	3	12
Perineural invasion	Absent	167	85	82	0.173
	Present	5	1	4
Tumor regression grade	Grade 0–1	37	5	32	<0.001*
	Grade 2–3	118	66	52
	Grade 4	17	15	2

Note: *Statistically significant.

Correlation of UBD Immunoexpression with Clinicopathological Parameters

The immunohistochemical stain showed UBD immunoactivity was significantly higher in CCRT-responders than in CCRT non-responders (Figure 2). Table 2 shows the correlation between UBD immunoexpression and clinicopathological parameters. Upregulated UBD expression was associated with lower pre-CCRT tumor T stage (P = 0.009), lower post-CCRT tumor T stage (P < 0.001), lower post-CCRT nodal stage (P < 0.001), less vascular invasion (P = 0.015), and better tumor regression (P < 0.001).

Figure 2

Immunohistochemical expression of UBD. Representative images of rectal cancer exhibiting high UBD expression among CCRT responders. (A) Normal colon mucosa with no UBD expression. (B) CCRT non-responder with low UBD expression. (C) CCRT responder with high UBD expression.

Survival Analyses and Clinical Implications of UBD Expression

Based on univariate analysis, high UBD expression was correlated with better disease-free survival (DFS) (P < 0.0001), local recurrence-free survival (LRFS) (P < 0.0001), and metastases-free survival (MeFS) (P < 0.0001). (Table 3 and Figure 3) Lower post-CCRT tumor T stage and better tumor regression grade were associated with better DFS, LRFS, and MeFS (all P < 0.009). Lower pre-CCRT tumor T stage was associated with better DFS (P = 0.086). Pre-CCRT nodal metastases correlated with worse DFS (P = 0.001) and worse LRFS (P = 0.007). Vascular invasion was correlated with worse DFS (P = 0.0029) and LRFS (P = 0.0028). Based on multivariate analysis (Table 4), high UBD expression was associated with superior DFS (P < 0.001), LRFS (P = 0.01), and MeFS (P = 0.004). Worse tumor regression grade and pre-CCRT nodal metastases were correlated with inferior DFS (P = 0.078 and P = 0.023, respectively).

Table 3

Univariate Log Rank Analysis for Important Clinicopathological Variables and UBD Expression

Parameter		No. of Cases	DFS		LRFS		MeFS
			No. of Events	p-value	No. of Events	p-value	No. of Events	p-value
Gender	Male	108	34	0.6344	20	0.2250	17	0.3520
	Female	64	20		7		14
Age	<70	106	35	0.6999	18	0.6615	20	0.7427
	≥70	66	19		9		11
Pre-Tx tumor status (Pre-T)	T1-T2	81	20	0.0486*	10	0.2261	11	0.1745
	T3-T4	91	34		17		20
Pre-Tx nodal status (Pre-N)	N0	125	33	0.0010*	15	0.0070*	19	0.0973
	N1-N2	47	21		12		12
Post-Tx tumor status (Post-T)	T1-T2	86	15	0.0002*	7	0.0040*	8	0.0033*
	T3-T4	86	39		20		23
Post-Tx nodal status (Post-N)	N0	123	35	0.2338	16	0.1320	20	0.4634
	N1-N2	49	19		11		11
Vascular invasion	Absent	157	45	0.0029*	21	0.0028*	27	0.4470
	Present	15	9		6		4
Perineural invasion	Absent	167	51	0.0647	25	0.0940	30	0.9083
	Present	5	3		2		1
Tumor regression grade	Grade 0–1	37	23	<0.0001*	10	0.0090*	14	0.0006*
	Grade 2~3	118	30		17		16
	Grade 4	17	1		0		1
UBD expression	High Exp.	86	8	<0.0001*	4	<0.0001*	4	<0.0001*
	Low Exp.	86	46		23		27

Note: *Statistically significant.

Abbreviations: DFS, disease-free survival; LRFS, local (pelvic) recurrence-free survival; MeFS, metastasis-free survival.

Figure 3

Kaplan–Meier analysis showed UBD overexpression was associated with superior disease-specific survival (A), local recurrence-free survival (B) and metastasis-free survival (C).

Table 4

Multivariate Analysis for Important Clinicopathological Variables and UBD Expression

Parameter	DFS			LRFS			MeFS
	HR	95% CI	p-value	HR	95% CI	p-value	HR	95% CI	p-value
Tumor regression grade	0.604	0.345–1.059	0.078	0.616	0.280–1.358	0.229	0.607	0.297–1.239	0.170
UBD expression	5.019	2.167–11.627	<0.001*	4.752	1.454–15.524	0.010*	5.276	1.706–16.313	0.004*
Vascular invasion	1.218	0.553–2.681	0.624	1.764	0.643–4.838	0.270	-	-	-
Post-Tx tumor status (Post-T)	1.535	0.802–2936	0.196	1.518	0.610–3.779	0.370	1.576	0.672–3.694	0.295
Pre-Tx nodal status (Pre-N)	2.209	1.113–4.386	0.023*	2.162	0.923–5.066	0.076	-	-	-
Pre-Tx tumor status (Pre-T)	0.863	0.446–1.670	0.663	-	-	-	-	-	-

Note: *Statistically significant.

Abbreviations: DFS, disease-free survival; LRFS, local (pelvic) recurrence-free survival; MeFS, metastasis-free survival.

UBD Positively and Negatively Correlated Genes

The top 200 genes whose expression is most positively correlated () or negatively correlated () with UBD in rectal adenocarcinoma from the same transcriptome database (GSE35452) were examined to further predict the biological activities of UBD-interacting networks. The top 10 co-upregulated genes are CXCL11, COL12A1, CTHRC1, COL13A1, PHTF2, EGFL6, CD47, PPT1, BCAT1, and LARP6. The top 10 co-downregulated genes are MOGAT2, STARD5, CD300LG, PTPN18, RASAL1, FAM22A, DNAJC4, GNG13, PHKB, and FAM98C. These genes of selection were also submitted to perform enrichment ontology analysis of their molecular function. We identified the top 200 positively correlated genes are enriched to several pathways affecting molecular functions. Platelet-derived growth factor binding (GO:0048407, fold enrichment:76.4), extracellular matrix (ECM) structural constituent conferring tensile strength (GO:0030020, fold enrichment: 38.4), and IgG binding (GO:0019864, fold enrichment: 38.2) are the top three GO terms associated with UBD upregulation (). Meanwhile, enoyl-CoA hydratase activity (GO:0004300, fold enrichment: 41.3) is the only GO term for those molecular functions negatively associated with UBD ().

Discussion

UB has been established as a versatile molecule that is involved in both proteolytic and non-proteolytic activities. Ubiquitination is a multi-step process that involves an E1 ubiquitin-activating enzyme, E2 ubiquitin-conjugating enzymes, and E3 ubiquitin ligases to attach ubiquitin to a target protein.20 At least four ubiquitin moieties linked by lysine target a wide range of cytosolic and nucleoplasmic proteins for proteasomal destruction. Polyubiquitylated substrates are transported to the proteasome by various ubiquitin receptors, where ubiquitins are recycled and broken into tiny peptides by proteasomic proteases. UB governs several physiological activities indirectly through proteolysis, including cell cycle progression, differentiation, and development, cell responses to extracellular effectors and stress, induction of inflammatory and immunological responses, and DNA damage repair or tolerance.21,22 The UBD, also called FAT10, is a ubiquitin-like protein and was first discovered in 1996.23 While ubiquitin is recycled from degraded target proteins, UBD degrades with its target, resulting in a relatively short half-life.24,25 It is a protein belonging to the immune system, which can be strongly upregulated by pro-inflammatory cytokines.10 The description of the role of UBD in cancer in the published literature is contradictory. Several clinical studies have shown that UBD could confer malignant characteristics to non-tumorigenic cells and enhance the malignant-related characteristics of cancer cells.26,27 Tumor UBD expression shows some tissue specificity, with transcriptional upregulation observed in liver, uterine cervix, ovarian, pancreatic, gastric and small intestine adenocarcinomas, but not in thyroid, prostate or kidney cancers.28,29 Ji et al demonstrated that overexpression of UBD in gastric cancer is correlated with metastasis and tumor staging, and both UBD mRNA and protein levels were identified as independent prognostic factors for the disease.15 Increased UBD has been also linked to mutant p53 expression, which could activate UBD expression and aid gastric cancer growth indirectly.30 There are also several studies exploring the role of UBD in colon cancer. Yan et al31 demonstrated that increased cytoplasmic UBD immunohistochemical staining was significantly associated with depth of cancer invasion, lymph node metastasis, distant metastasis, tumor histologic grade, advanced clinical stage, and Ki-67 proliferative index. After radical surgery for stage II and III colon cancer, individuals with UBD-positive tumors had a considerably higher disease recurrence rate and shorter survival than patients with UBD-negative tumors. These findings suggest that UBD may contribute to colon carcinogenesis and can be used as a new prognostic indicator for predicting recurrence in stage II and III patients following curative surgery. Another study revealed that UBD-positive expression was a significant, independent predictive high-risk factor for overall survival and DFS in stage IIB–IIC colon cancer patients treated with 5-fluorouracil-based adjuvant chemotherapy.32 However, our study yielded different results. Our study revealed that high UBD expression correlated to better DFS, LRFS, and MeFS in rectal cancer patients receiving CCRT followed by curative resection. The differences between this study and others are that our patients have received radiation therapy while others have not. This is the first study to explore the relationship between UBD and chemoradiotherapy. In a study presented by Yu Wang et al,33 they found that UBD expression was significantly positively correlated with all analyzed immune cells in skin cutaneous melanoma (B cell, CD8+T cell, CD4+T cell, macrophage, neutrophil, and dendritic cell). In this study, high UBD expression predicted better outcomes. Since several studies proposed that tumor-infiltrating immune cells in the tumor microenvironment may improve the response to neoadjuvant CCRT and improve the recurrence-free survival in rectal cancer patients,34–37 the relationships between UBD, tumor-infiltration immune cells and response to neoadjuavant CCRT are worth further exploration. In fact, the conclusion that UBD promotes carcinogenesis is still open to question. Lukasiak et al38 suggested that overexpression of UBD in colon cancer is a consequence of a proinflammatory immune response and so there is no reason to propose that UBD contributes to carcinogenesis. In fact, UBD failed to exert transforming activity in their study. Therefore, it is difficult to conclude whether UBD overexpression in more advanced cancer is the cause or the effect.

In our previous study, we use the same dataset to evaluate the genes associated with maintenance of gastrointestinal epithelium (GO: 0030277) and validate MUC2 as a significant survival biomarker by the same patient cohort.39 Based on this, we believe the pathogenesis of rectal cancer is complex and all biomarkers are different in nature and biology. Searching other possible molecular biomarkers is mandatory. In this study, we use the same transcriptome dataset, where genes related to the protein modification process (GO:0036211) were assessed to better understand the role of UBD in rectal cancer.

To further predict the molecular functions of UBD, the top 200 genes co-upregulated () or co-downregulated () with UBD in colorectal adenocarcinoma from the transcriptomic database were evaluated. We found that the genes negatively associated with UBD expression are significantly enriched to enoyl-CoA hydratase activity (GO:0004300, Fold enrichment:41.3) (). Moreover, among the top 10 pathways to which the positively correlated genes are enriched (), ECM-related pathways accounted for four places: ECM structural constituent conferring tensile strength (GO:0030020, fold enrichment: 38.4), ECM structural constituent (GO:0005201, fold enrichment: 19.4), ECM structural constituent conferring compression resistance (GO:0030021, fold enrichment: 19.1), and ECM binding (GO:0050840, fold enrichment: 13.1). The ECM proteins confer distinct biochemical and biophysical properties that influence cell phenotype. The ECM’s composition and organization are spatiotemporally regulated to control cell behavior and differentiation, but dysregulation of ECM dynamics may lead to the development of cancer.40 ECM does not only influence malignancy and growth of the tumor but also its response to therapy.41 The exact mechanisms to link UBD and the ECM-related pathway need further exploration and analysis. There are still several limitations to the existing studies. First, this is a single institution’s retrospective analysis of rectal cancer patients who received preoperative CCRT. Second, this study lacks a validation cohort. Third, the role of UBD function in cancer progression is generally indirect. Fourth, if the ideal cutoff for high UBD expression is not established, clinical assessment will be useless. And finally, there were not enough samples. Therefore, more multicenter investigations should be conducted to confirm the significance of UBD expression.

Conclusion

In conclusion, our study demonstrated that high UBD expression has been linked to less advanced rectal cancer characteristics and is a unique predictive indicator of better patient outcomes for rectal cancer patients undergoing neoadjuvant CCRT.