Abstract
Keywords
Introduction
Colorectal cancer (CRC) is one of the most common malignant tumors with an incidence rate ranking third and a mortality rate ranking second in the world. 1 Studies demonstrated that even in different patients at the same disease stage, the prognosis varies greatly, which may be related to changes in the tumor microenvironment. 2 High mobility group protein B-1 (HMGB-1) is a multifunctional non-histone protein located in the nucleus, involved in DNA repair, transcription, and genome stability. 3 HMGB-1 is highly expressed in CRC, being involved in lymph node metastasis and distant metastasis in patients with CRC, and can be used as a predictor of poor prognosis. 4 However, the mechanism by which HMGB-1 affects the progression of CRC remains unclear.
It has been reported that HMGB-1 not only participates in the alteration of the tumor microenvironment but also regulates tumor immune response. 5 The intestinal mucosal barrier is a part of the intestinal tumor microenvironment, and its destruction can accelerate the invasion and metastasis of CRC cells. 6 Previous studies confirmed that the destruction of the intestinal mucosal barrier involves intestinal flora and intestinal immune response. For example, the reduction of probiotic flora and the increase of pathogenic flora can lead to the destruction of the intestinal mucosal barrier, thereby promoting the occurrence and development of CRC.7,8 Tight junction protein are surface indicators of the intestinal mucosal barrier, and changes in their levels are related to the occurrence and development of CRC. 9 In addition, the increase of the expression of intestinal inflammatory mediators such as tumor necrosis factor-α (TNF-α) and interferon-γ (IFN-γ) can further affect the immune barrier of the intestinal mucosa, thereby promoting the initiation and progression of CRC. 10 Therefore, by examining the relationship between HMGB-1, intestinal flora, and the expression levels of Zonula occludens −1 (ZO-1), occludin, TNF-α and IFN-γ, this study explores the potential mechanism of HMGB-1 affecting the occurrence and development of CRC.
Materials and Methods
Clinical Information
We screened a total of 294 patients from the dataset GSE12945 (including 62 CRC patients), GSE17536 (including 177 CRC patients), and GSE17537 (including 55 CRC patients). Excluding incomplete data, 232 patients were finally included, 120 males and 112 females with the age range 35-82 years (median age 65.3 years), within which 118 patients were >65 and 114 patients were ≤65 years. According to the TNM staging (AJCC No. Version 11 ), 40 cases were diagnosed with stage I, 85 cases with stage II, 85 cases with stage III, and 22 cases with stage IV CRC. A retrospective analysis of 120 patients with CRC admitted to hospital from January 2016 to December 2016 was performed. The study inclusion criteria were as follows: (1) Concurrent chemoradiotherapy with radical surgery; (2) Pathology consistent with CRC; (3) Written informed consent has been obtained; (4) No previous malignant tumors or other accompanying malignant diseases; (5) Complete medical records, imaging data, and follow-up data. The study exclusion criteria were as follows: (1) CRC combined with other malignant tumors; (2) Surgical contraindications; (3) Postoperative survival ≤3 months; (4) Concurrent hematological diseases, acute injury, chronic inflammation, or intestinal perforation; (5) Incomplete patient information. According to the inclusion and exclusion criteria, a total of 110 CRC patients were included in this study, 71 males and 39 females with a median age of 59.61 years, among which 53 patients were ≤60 years, and 57 patients were >60 years. According to the TNM staging (AJCC 8th edition), 10 cases were diagnosed with stage I, 51 cases with stage II, and 49 cases with stage III CRC. Based on the degree of differentiation there were identified: 35 cases of poor differentiation, 62 cases of moderate differentiation, and 13 cases of high differentiation. There were 49 cases of lymphatic metastasis and 61 cases of no lymphatic metastasis. The median tumor diameter was 5.11 cm, 60 cases with a diameter of ≤5 cm and 50 cases with >5 cm. There were 17 cases positive for vascular tumor thrombus and 93 cases negative for vascular tumor thrombus. There were 23 cases with nerve invasion and 87 cases without nerve invasion.
Immunohistochemical Staining
All specimens were fixed in 10% formalin, dehydrated, transparent, soaked in wax, embedded in paraffin, dehydrated, and cut into 3 μm sections on a microtome. Next, antigen retrieval and inactivation of endogenous peroxidase were performed. Then, the samples were incubated in a blocking buffer (goat serum SAP-9100, Chian, Bei Jing, Beijing Zhong Shan-Golden Bridge Biological Technology Co., Ltd), followed by overnight incubation at 4 °C with primary antibody against HMGB-1 (EPR3507, Chian, Shanghai, Abcam; 1:350), ZO-1 (EPR19945-296, Chian, Shang Hai, Abcam; 1:250), occludin (EPR20992, Chian, Shang Hai, Abcam; 1:200), IFN-γ (EPR21704, Chian, Shang Hai, Abcam; 1:250), and TNF-α (P/T2, Chian, Shang Hai, Abcam; 1:100) added dropwise onto the tissue. The next day after the pre-warming of the slides, a secondary antibody (SE134, Chian, Shanghai, Abcam; 1:200) was added. The developing solution was applied to the sections on the slides. In the end, the sections were counterstained, mounted, and coverslips were sealed. The specimens were examined under a microscope by 2 experienced pathologists in a double-blind method. Three high-power fields of view were randomly selected for image collection. We adopted a semi-quantitative scoring method to investigate the expression of HMGB-1, ZO-1, occludin, IFN-γ, and TNF-α proteins using the following point scale: 0—no positive cells or <1% of positive cells, 1—1-10% of positive cells, 2—10-50% of positive cells, 3—50-75% of positive cells, 4— > 75% of positive cells. Color rendering evaluation was as follows: 0—no coloring, 1—light yellow, 2—brownish yellow, and 3—brown. When 2 points were multiplied together, 0-4 points were negative, and 5-12 points were positive.
Serum Specimen Enzyme-Linked Immunosorbent Assay (ELISA)
ELISA kit (E-EL-H1554c, Chian, Wu Han, Elabscience) was used to detect serum HMGB-1 levels in CRC patients. 100 mL/well of standards and samples were loaded into a 96-well plate, and after incubation with biotinylated antibodies, streptavidin-conjugated horseradish peroxidase (HRP) was added to each well and mixed with the HRP substrate solution. The optical density (OD) values were read by a microplate reader at 450 nm.
16S rDNA Sequencing
To detect the intestinal flora of the 2 groups of patients, CTAB or SDS method was used to extract the genomic DNA from patient stool samples. 16S V3-V4 region was amplified and PCR products were verified by agarose gel electrophoresis and purified using magnetic beads. Further, they were quantified by enzyme labeling, mixed in equal amounts according to the concentration of the PCR product, verified on 2% agarose gel, and recovered. TruSeq® DNA PCR-Free Sample Preparation Kit was used to prepare DNA library that was quantified with Qubit and qPCR. After quality control, whole-genome sequencing on NovaSeq 6000 system was performed. For the bioinformatic analysis, the original data was spliced and filtered. UPARSE algorithm was used to cluster the sequences at 97% similarity, and after chimera filtering and Operational Taxonomic Units (OTU) clustering the sequences were analyzed for species diversity (alpha and beta diversity) and abundance. Linear discriminant analysis effect size (LEfSe) analysis was performed using LEfSe software with a screening value of 4. The functional information of the whole genome of prokaryotic organisms in the KEGG database annotated by UProC and PAUDA methods was mapped to the SILVA database to realize the functional annotation with Tax4Fun.
Statistical Analysis
The data statistics were analyzed by SPSS26.0 software. The count data were described by percentages, and the χ2 test of 4 table data was used for comparison. The cutoff value of HMGB-1 was calculated by the ROC curve. The prognostic significance of HMGB1 expression on survival in CRC was carried out by Kaplan-Meier survival analysis. The log-rank test was used to conduct univariate analysis on the above indicators, and Cox regression was used to conduct multivariate analysis on the indicators with meaningful univariate analysis. The Pearson method was used to analyze the correlation. P < .05 was considered statistically significant.
Results
Elevated HMGB-1 Levels are Closely Associated With Poor Prognosis in CRC
To understand the impact of HMGB-1 level changes in CRC patients, ROC curve analysis was conducted to compare the predictive ability of HMGB-1 levels in 232 patients from the database for survival. The results showed that the area under the curve of HMGB-1 was 0.587, with a cutoff value of 10.24 ng/mL (Figure 1A). The CRC patients were divided into high (HMGB-1H, > 10.24 ng/mL) and low expression groups (HMGB-1L, ≤ 10.24 ng/mL) with 187 (80.6%) and 45 cases (19.4%), respectively. Compared to the HMGB-1L group, the TNM stage of the CRC patients in the HMGB-1H group was of a higher grade, and they were more likely to relapse, with more deaths and shorter overall survival (OS) and disease-free survival (DFS) (P = .010, P = .006, P < .001, P = .005, and P = .035). There was no significant difference between the 2 groups in gender and age (P > .05) (Table 1). Moreover, the results of Cox regression analysis showed that the high expression of HMGB-1 and recurrence were independent risk factors for the death of CRC patients (P = .013, P < .001), with the death risk in the HMGB-1H group being 6.204 times that of the HMGB-1L group (Table 2). Since December 2021, 170 CRC patients are still alive, with a 5-year OS rate of 73.3%, among which 67.9% and 95.6% for the high and low HMGB-1 expression group, respectively (P = .002) (Figure 1B).

(A) The ROC curve of HMGB-1 levels. (B) Analysis of the overall survival in colorectal cancer and HMGB-1 levels. (C) HMGB-1 positive group. (D) HMGB-1 negative group. (E) The relationship between HMGB-1 expression and the overall survival rate of patients with colorectal cancer.
The Association of HMGB-1 With Clinical Features and Prognosis of Colorectal Cancer.
Abbreviations: HMGB-1, high mobility group box B-1; OS, overall survival; DFS, disease-free survival.
P < .05 was considered statistically significant.
Multivariate Cox Regression Analysis of Death From Colorectal Cancer.
Abbreviations: HMGB-1, high mobility group box B-1; OS, overall survival; DFS, disease-free survival.
P < .05 was considered statistically significant.
To further examine the influence of HMGB-1 expression on CRC pathophysiology, 110 tissue samples from CRC patients from our hospital were selected for immunohistochemical detection. The results demonstrated that HMGB-1 was mainly expressed in the nuclei of CRC cells, with a brown or brownish-yellow color. The positive expression rate was 66.4% (73/110) (Figure 1C), and the negative expression rate was 33.6% (37/110) (Figure 1D). Furthermore, the HMGB-1 expression was related to the degree of differentiation, lymph node metastasis, vascular tumor thrombus, and TNM stage in CRC patients (P < .05), but not to gender, age, tumor diameter, and nerve invasion (P > .05) (Table 3). The results of Cox regression analysis revealed that HMGB-1 expression and vascular tumor thrombus were independent risk factors for poor prognosis in CRC patients (P < .05) with the death risk in the HMGB- 1H group being 0.437 times that of the HMGB-1L group (P < .036). The death risk of CRC patients with vascular tumor thrombus was 0.461 times that of patients without vascular tumor thrombus (P < .026) (Table 4). Based on the follow-up until December 2021, 62 of the 110 patients survived, and the 5-year OS was 56.4% with 45.2% and 78.4% for the HMGB-1-positive and negative patients, respectively (P = .003) (Figure 1E). Therefore, our findings demonstrate that the increase in HMGB-1 level is closely related to the poor prognosis of CRC.
Comparison of HMGB-1 Expression With Clinical Data and Pathological Characteristics.
Abbreviations: HMGB-1, high mobility group box B-1.
P < .05 was considered statistically significant.
Multivariate Cox Regression Analysis of Prognostic Factors for Colorectal Cancer.
Abbreviations: HMGB-1, high mobility group box B-1.
P < .05 was considered statistically significant.
Increase in HMGB-1 Expression Affects the Redistribution and Function of Intestinal Flora in Patients With CRC
To explore whether the change in HMGB-1 level is associated with the alteration of intestinal flora, we performed ELISA on serum samples from randomly selected 88 of 110 CRC patients. We found that the average value of HMGB-1 expression was 10.33 ng/mL, which was similar to the database cutoff value (10.24 ng/mL). The 88 samples were divided into high expression group (HMGB-1H) and low expression group (HMGB-1L) with 13.09 ± 2.456 ng/mL and 7.381 ± 1.634 ng/mL, respectively (P < .001) (Figure 2A).

(A) Comparison of HMGB-1 levels between 2 groups. (B) Venn diagram of intestinal flora among 2 groups. (C) Histogram of relative abundance of genus level in 2 groups. (D) Alpha diversity analysis (Chao 1, Observed species, Shannon, and Simpson) of differential bacteria between groups. (E) Unweighted UniFrac distance analysis in 2 groups. (F) Weighted UniFrac distance analysis in 2 groups. (G) Linear discriminant analysis effect size (LEfSe) for differentially abundant taxa from the 2 groups. (H) Taxonomic cladogram from the LEfSe showing differences in fecal taxa. (I) Tax4fun prediction of the differential microbiome.
Based on the ELISA results, we selected 10 samples from the HMGB-1H group and HMGB-1L group for 16S rDNA analysis. Statistical analysis of OTUs at a 97% similarity level revealed that apart from 468 common species there were 259 and 166 unique species of bacterial flora in the HMGB-1H and HMGB-1L groups, respectively, indicating that the 2 groups have their own unique flora abundance (Figure 2B). The species diversity and richness suggested by alpha and beta diversity analysis were higher in the HMGB-1H group than those of the HMGB-1L group (P < .05) (Figure 2D-F). The species composition was different between the 2 groups. The HMGB-1H group was enriched with Staphylococcus, Faecalibacterium, and Enterococcus, while the HMGB-1L group was enriched with Proteus, Ekmansella, Escherichia coli, and Shigella bacilli (Figure 2C). The results of LEfSe analysis demonstrated that there was no significant difference between the HMGB-1H group and the HMGB-1L group; only Providencia content was found to be significantly increased in the HMGB-1L group. Elevation of HMGB-1 level had little effect on gut microbiota (Figure 2G-H). To explore how the HMGB-1 level affects the gut microbiota, we analyzed the functional predictions of the HMGB-1H and HMGB-1L group differential flora. The results showed that the functional enrichment of the differential flora between the groups was mostly concentrated in DNA damage repair, metabolic changes such as glycolysis and amino acid metabolism, and cell motility (P < .05) (Figure 2I). The above results show that the increase of HMGB-1 levels affects the redistribution and function of intestinal flora.
Increased HMGB-1 Levels Affect the Expression of Proteins Related to the Intestinal Barrier in CRC Patients
To examine the impact of HMGB-1 level changes on the intestinal barrier, we selected ZO-1, occludin, TNF-α, and IFN-γ for immunohistochemical detection. Our results showed that ZO-1, occludin, TNF-α, and IFN-γ were diffusely distributed in the membrane and cytoplasm of CRC cells (brown or brownish yellow), and a small amount was expressed in the nucleus (Figure 3). Among HMGB-1-positive patients, 14 cases (19.2%) were positive for ZO-1 protein, and 59 cases (80.8%) were negative (P = .034). The positive and negative expression rates of occludin in HMGB-1-positive patients were 27.4% and 72.6%, respectively (P = .003). There were 56 cases of TNF-α positive patients, accounting for 76.7% and 17 cases of TNF-α negative patients, accounting for 23.3% (P = .015). In HMGB-1 positive patients, 48 cases were IFN-γ positive, accounting for 65.8%, and 25 patients were IFN-γ negative, accounting for 34.2% (P = .001). Our findings demonstrated that the expressions of ZO-1 and occludin in HMGB-1 positive patients were lower than those in HMGB-1 negative patients (R = 0.035, R = 0.003, P < .05), while the expressions of TNF-α and IFN-γ were higher than those in negative patients (R = 0.016, R = 0.001, P < .05) (Table 5), indicating that the increase in HMGB-1 levels may affect the intestinal mucosal barrier and immune barrier in CRC patients.

The expression of ZO-1, occludin, TNF-α, and IFN-γ in colorectal cancer tissues (HE, 400 ×). (A) ZO-1 positive group. (B) ZO-1 negative group. (C) Occludin positive group. (D) Occludin negative group. (E) TNF-α positive group. (F) TNF-α negative group. (G) IFN-γ positive group. (H) IFN-γ negative group.
The Association of HMGB-1 With ZO-1, Occludin, TNF-α, and IFN-γ in Colorectal Cancer.
Abbreviations: HMGB-1, high mobility group box B-1; ZO-1, Zonula occludens-1; TNF-α, tumor necrosis factor-α; IFN-γ, interferon-γ.
P < .05 was considered statistically significant.
Discussion
HMGB-1 plays a vital role in the occurrence and development of CRC. Previous studies confirmed that the high expression of HMGB-1 is a poor prognostic factor in CRC. In this study, we reported that the investigated patients with high expression of HMGB-1 had a higher TNM stage. Moreover, their OS and DFS shortened, and they were more likely to relapse and metastasize. Based on the HMGB-1 cutoff value, they were divided into 2 groups, HMGB-1H and HMGB-1L. 16SrDNA sequencing demonstrated that the intestinal flora between these 2 groups had different diversity and richness, among which Providencia was substantially increased in the HMGB-1L group. In addition, this study also found that the increase in HMGB-1 expression may be related to the damage of the intestinal mucosal physical barrier and the change of the local immune barrier.
The analysis of the patient clinical phenotype revealed that the expression of HMGB-1 was related to the death, recurrence, OS, DFS, and TNM stage of CRC patients, and HMGB-1 expression is an independent risk factor for the death and recurrence of CRC patients. The 5-year OS of patients with low HMGB-1 expression was higher than that of CRC patients with high HMGB-1 expression. We further verified that the HMGB-1 expression was related to the degree of differentiation, lymph node metastasis, vascular tumor embolism, and TNM stage of CRC patients. These results correspond to those in the database. The results of this study are also in line with many previous studies, and they confirmed that HMGB-1 can be used as a predictive factor for the CRC prognosis.12–15 Our research confirms that HMGB-1 is closely related to CRC stage, lymph node metastasis and differentiation, and may be an important marker of its occurrence and development, which will help to improve the accuracy of CRC diagnosis.
Furthermore, we performed 16S rDNA sequencing on the feces of CRC patients according to the HMGB-1 expression. The Venn diagram (Venn) showed that the intestinal flora is a complex micro-ecological environment, and there were many similar species between low and high HMGB-1 expression groups. However, the abundance of specific bacterial flora in the HMGB-1H group was still higher than that in the HMGB-1L group. It has been shown previously that the reduction or depletion of intestinal flora could significantly reduce CRC liver metastasis and primary liver cancer, suggesting that the reduction of the abundance of intestinal flora may be beneficial to the development of CRC. 16 In addition, we also know from the alpha and beta diversity that the species diversity and richness of the HMGB-1H group were greater than those of the HMGB-1L group. It suggests that the reduction of flora species was beneficial to the development of CRC. We used the LEfSe histogram and evolutionary map to perform a high-dimensional classification analysis of the intestinal flora in the 2 groups. At the genus level, only Providencia was enriched in the HMGB-1L group, which suggested that the increase of Providencia content may inhibit the development of CRC. It has been previously reported that the increase of Providencia could lead to massive secretion of lethal swelling toxin, causing cell swelling, and blocking eukaryotic cell proliferation in the G2/M phase, resulting in cell death. 17 In vivo studies have also found that Providencia could translocate from the intestinal lumen to extraintestinal organs, confirming its potential role in disseminating infection. 18 Although previous studies on Providencia considered it pathogenic, the research did not determine a causal link between Providencia species and cancer. Thus, it needs further investigation whether Providencia has a tumor suppressor effect, its content is reduced due to the increase in HMGB-1 level, or produces specific metabolites limiting the CRC progression.
An intact mucosal barrier and mucosal immune system are critical for maintaining intestinal homeostasis. 19 ZO-1 and occludin play an important role in the intestinal inflammatory response and immune response by maintaining the physical barrier function of the intestinal mucosa. 20 This study found that the ZO-1 and occludin proteins were reduced in HMGB-1-positive patients, which is consistent with previous research results.21–24 The possible mechanism may be that HMGB1 stimulates Toll-like receptor-4 (TLR4) leading to a decreased expression of ZO-1 and occludin, thereby damaging endothelial cells and changing their permeability. 25 Alternatively, HMGB-1 may stimulate the production of trimethylamine nitrogen oxides, which decreases the expression of ZO-2 and occludin and destroys the tight junction between cells, increasing the physical barrier damage of intestinal mucosa and finally leading to CRC. 26 It has been demonstrated that HMGB-1 destroys intercellular connections in a paracrine or autocrine manner, and when the receptor for advanced glycation end products (RAGE) is silenced, the expressions of ZO-1 and occludin increased significantly, thereby weakening the disruption of the intestinal mucosa. 27
Imbalances in the structure of the gut barrier may unleash an uncontrolled immune response in the gut microenvironment. This study showed that HMGB-1 could increase TNF-α and INF-γ expression in CRC. Similar to the previous studies, it was found that TNF-α level in the serum of CRC patients was substantially increased, and thus anti-TNF-α agents could be applied in CRC treatment,28–31 potentially due to the activation of EGFR/ERK1/2 and Wnt/β-catenin signaling. 32 HMGB-1 can increase IFN-γ production by promoting the activation of CD8[ + ] T cells induced by dendritic cells, and IFN-γ is positively correlated with HMGB1 at the mRNA level. 33 Therefore, HMGB-1 can participate in inflammatory and immune responses, affecting CRC microenvironment and promoting its initiation and progression.
However, this study also has some limitations. It did not explore the relationship between HMGB-1 and CRC molecular typing, and its influence on the local recurrence rate and distant metastasis rate. Moreover, the association of HMGB-1 with inflammatory markers (TNF-α, and IFN-γ) and the specific mechanism of ZO-1 and occludin should be studied more in-depth.
Conclusion
In summary, we demonstrated that the increase in HMGB-1 levels is closely related to the poor prognosis of CRC and may lead to the decrease of Providencia content, ZO-1, and occludin expression, and the increase of TNF-α, and IFN-γ expression. This can cause damage to the intestinal mucosal barrier ultimately leading to the occurrence and development of CRC.
Footnotes
Abbreviation
Data Availability
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by Training Project of Excellent Clinical Medical Talents in Hebei Provincial Government Funding (No.4).
Statement of Human Rights
This study was approved by the Ethics Committee of the First Affiliated Hospital of Hebei Northern University (K2021125), and all patients signed an informed consent form.
