- Department of Radiation Oncology, The Affiliated Cancer Hospital of Zhengzhou University & Henan Cancer Hospital, Zhengzhou, China
Background: Lung cancer is a significant challenge to human health. Members of the high mobility group (HMG) superfamily (HMGB proteins) are implicated in a wide variety of physiological and pathophysiological processes, but the expression and prognostic value of HMGB family members in non-small cell lung cancer (NSCLC) have not been elucidated.
Methods: In this study, ONCOMINE, UALCAN, GEPIA, Kaplan–Meier Plotter, starBase, OncomiR databases, and GeneMANIA were utilized to evaluate the prognostic significance of HMGB family members in NSCLC.
Results: HMGB2/3 expression levels were higher in NSCLC patients. HMGB1 expression was higher in lung squamous cell carcinoma (LUSC) and was lower in lung adenocarcinoma (LUAD) tissue than in normal lung tissue. HMGB2 expression was related to cancer stage. Increased HMGB1 mRNA expression levels were associated with improved lung cancer prognosis, including overall survival (OS), first-progression survival (FP), and post-progression survival (PPS). There was no significant association between HMGB2 levels and prognostic indicators. HMGB3 expression was associated with poorer OS. GeneMANIA and GO/KEGG pathway analysis showed that HMGB family members mainly associated with chromosome condensation, regulation of chromatin organization, and nucleosome binding in NSCLC. HMGBs expression were closely correlated with infiltrating levels of specific types of immune cells in NSCLC, especially Th2 cells, Th17 cells, and mast cells. hsa-miR-25-3p, hsa-miR-374a-3p, and hsa-miR-93-5p were significantly positively correlated with HMGB1, HMGB2, and HMGB3, respectively. However, hsa-miR-30a-5p was predicted to significantly negatively regulate HMGB3 expression.
Conclusion: Our study revealed that HMGB1 is positively related to the improved prognosis in NSCLC, and demonstrate that HMGB3 might be a risk factor for poorer survival of NSCLC patients.
Introduction
Lung cancer is one of the most common malignant tumors and poses a severe threat to human health. Recent data from the International Agency for Research on Cancer (IARC) revealed that lung cancer accounted for the highest number of cancer-related deaths worldwide, and was responsible for nearly 1.8 million deaths in 2020, more than any other type of cancers (Baltayiannis et al., 2013; Bray et al., 2018). Small cell lung cancer (SCLC) and non-small cell lung cancer (NSCLC) are the two main histological types of lung cancer. NSCLC is the most common type of lung cancer, accounting for 85% of all lung cancer cases. NSCLC is further subtyped into adenocarcinomas (ADC; which comprise ∼40–50% of NSCLC), squamous cell carcinomas (SCC; which comprise ∼20–40% of NSCLC), and large cell carcinomas (LGC; which comprise ∼20% of NSCLC) (Wang Y et al., 2019; Relli et al., 2019). Despite the significant efforts devoted to the study of NSCLC, the exact underlying mechanisms that contribute to the development, progression, and recurrence of NSCLC remain to be clearly elucidated. Therefore, it is urgent to identify new highly specific and sensitive biomarkers and new molecular targets to improve care and outcomes for patients with NSCLC.
HMGB proteins are members of the high mobility group (HMG) superfamily. This family includes four members: HMGB1, HMGB2, HMGB3, and HMGB4 (Taniguchi et al., 2018; Niu et al., 2020; Wen et al., 2021). HMGB1/2/3 proteins are composed of two DNA-binding HMGB-boxes and an acidic tail. HMGB4 lacks the acidic tail and its nucleotide sequence is less well-conserved than the sequences of the other HMGB family members (Taniguchi et al., 2018; Niu et al., 2020; Wen et al., 2021). Since research into HMGB4 is limited, we mainly focused on the role of HMGB1/2/3 in NSCLC. HMGBs have been implicated in a wide variety of physiological and pathophysiological processes, involving infection, injury, and inflammation by regulating DNA replication, transcription, recombination, and repair (Niu et al., 2020). HMGBs also are known to play important role in tumorigenesis and cancer progression, including NSCLC (Taniguchi et al., 2018; Qiu et al., 2020; Wen et al., 2021).
Abnormal expression of some HMGB family members has been reported in NSCLC, and several HMGB genes have prognostic value in NSCLC. Qiu et al. (2020) demonstrated that HMGB1 is highly expressed in NSCLC tissue compared with normal tissue, and showed that HMGB1 facilitates NSCLC progression via activating the Wnt/β-catenin pathway. In addition, Wang et al. (2014) reported that HMGB1 was significantly and positively associated with tumor grade and stage in NSCLC, was positively associated with MMP-9, and may serve as a prognostic and predictive marker for NSCLC. Li S et al. (2018) showed a potential mechanism underlying the oncogenic role of HMGB2 in NSCLC development through the miR-181a-5p/HMGB2 axis. Similarly, Lou et al. (2021) found that HMGB2 protein expression was higher in tumor tissue compared with adjacent tissue, and correlated with lymph node (LYN) metastasis and advanced TNM stage in NSCLC, indicating a poor prognosis. HMGB3 is highly expressed in a variety of cancers, including breast cancer, gastric adenocarcinoma, lung cancer, bladder cancer, esophageal cancer, glioma, and prostate cancer (Niu et al., 2020; Wen et al., 2021). Song et al. (2019) and Xie et al. (2020) demonstrated that HMGB3 could accelerate the proliferation and colony formation of NSCLC cells, and reduced apoptosis of NSCLC cells regulated by circEPSTI1/miR-145, suggesting that HMGB3 may be a novel independent prognostic marker of worse outcome for NSCLC patients (Song et al., 2013).
These reports suggest that HMGB1/2/3 all act as oncogenes in NSCLC, but some other studies have contradicted these findings. Kang et al. (2013) and Luan et al. (2017) indicated that HMGB1 acts as tumor suppressor in pancreatic cancer and endometrial carcinoma. Furthermore, higher mRNA expression of HMGB1 led to a positive prognosis in lung cancer patients (Tang et al., 2010; Wu et al., 2018). In addition, Niu et al. (2020) reported that overexpression of HMGB2 associated with longer OS in cancers including invasive breast carcinoma (BRCA), cervical squamous cell carcinoma/endocervical adenocarcinoma (CESC), LUSC, rectal adenocarcinoma (READ), stomach adenocarcinoma (STAD), and thyroid carcinoma (THCA). In summary, the role of distinct HMGB family members remains unknown, especially regarding the development and progression of NSCLC. In the present study, bioinformatics analyses were performed to analyze the expression, prognosis, and mutations of different HMGB family members and to determine the relationship between HMGB family members and cancer stage in lung cancer patients. Additionally, we evaluated microRNA-HMGB interaction maps and explored the correlation of HMGB levels with tumor-infiltrating immune cells in NSCLC.
Materials and methods
Ethics statement
The study was approved by the Institutional Review Board of Zhengzhou University. All data were retrieved from the online databases, and written informed consent had already been obtained from patients prior to the deposition information or data in those databases.
ONCOMINE database
ONCOMINE is a cancer microarray database and web-based data-mining platform, which provides tools for genome-wide expression analysis (Rhodes et al., 2004). Using ONCOMINE, transcription levels of HMGB family members were compared between cancer tissues and corresponding adjacent normal tissues. A p-value of p < 0.05, a fold change of 2, and a gene rank in the top 10% were set as the significance thresholds.
UALCAN
UALCAN is an interactive web-portal used to perform to in-depth analyses of TCGA gene expression data. In our study, UALCAN was used to illustrate the distinct expression levels of HMBG family members in tumor and normal tissues. Differences were evaluated with Student’s t-tests, and a p-value of p < 0.05 was set as the significant level.
GEPIA
GEPIA is a web-based tool used deliver fast and customizable functionalities based on TCGA and GTEx data, and is freely available to all users. All the datasets on the server are computed by a standard pipeline and are compatible with each other. GEPIA covers a series of functions, including differential expression analysis, patient survival analysis, and correlation analysis (Tang et al., 2017). In this study, “Expression DIY-Stage Plot” was utilized to evaluate the relationship between HMGB family members and tumor stage. Differences were evaluated with Student’s t-test, and a p-value of p < 0.05 was set as the significant level.
Kaplan–Meier Plotter
The Kaplan-Meier method is the most common way to estimate survival times and probabilities. The Kaplan-Meier Plotter database can evaluate the effects of 54,000 genes on survival in 21 cancer types including breast, ovarian, lung, and gastric cancer (Nagy et al., 2018). Sources for the databases include GEO, EGA, and TCGA. In this study, KM Plotter was used to analyze the prognostic value of mRNA expression of HMGB family members in lung cancer in which cancer patients were split into high- and low- expression group, based on median expression values. Differences between groups were considered to be significantly different with a log-rank test p-value of p < 0.05.
Gene MANIA
Gene MANIA is a flexible, user-friendly web interface for searching genes that are related to a set of input genes, using a large set of functional association data. Association data include protein and genetic interactions, pathways, co-expression, co-localization, and protein domain similarity. Using Gene MANIA, we input “HMGB1; HMGB2; HMGB3” and analyzed by using the automatically selected weighting method in Homo sapiens.
Starbase
Starbase is a pan-cancer database used to identify interaction networks of lncRNAs, microRNAs, RNA-binding proteins, and mRNAs from large-scale CLIP-Seq data and tumor samples. In the “Target Gene” section, we obtained HMGB-related microRNAs by sequentially inputting “HMGB1, HMGB2, and HMGB3” genes into starBase. We then verified the correlation of each microRNA with HMGB genes in LUAD or LUSC under the “Pan-Cancer” section.
OncomiR
OncomiR is an online resource for exploring microRNA dysregulation in cancer. We retrieved a list of microRNAs closely associated with patient survival in LUAD or LUSC, with a log-rank p-value of p < 0.05 (univariate Cox analysis).
Results
Expression profiles of HMGB family members in NSCLC patients
To understand the relationship between HMGB family members and NSCLC, we first compared the definite genomic sites (Table 1) and protein structures of all members (Figure 1A). Details of the differences in HMGB gene sequences are shown in Supplementary Material S1. HMGB1/2/3 are highly similar in structure, containing two stable binding peptide motifs (A-box and B-box) and an acidic tail, while HMGB4 has no acidic tail. Subsequently, we focused on HMGB1/2/3 in patients with NSCLC. ONCOMINE and UALCAN were utilized to compare expression levels of HMGB family members in cancer vs. normal tissue samples. We identified dysregulated levels of HMGB gene transcripts in 20 different types of cancers, in comparison with normal samples, utilizing the ONCOMINE database (Figure 1B). The differences in HMGB gene expression in cancer and normal tissues indicated a potential role for HMGBs in tumorigenesis and cancer progression. In lung cancer samples, the mRNA expression levels of HMGB2 were slightly elevated, and HMGB3 was significantly up-regulated in cancer tissue vs. normal tissues in multiple datasets (Figure 1B). For example, in the Hou Lung dataset (Hou et al., 2010), a 5.974-fold increase in HMGB3 mRNA expression was identified in lung adenocarcinoma (LUAD) tissues vs. normal lung tissue (p = 3.38E-15; Table 2). The mRNA expression of HMGB2 in lung squamous cell carcinoma (LUSC) was higher by a fold change of 2.698 (p = 0.003) compared to normal lung tissue (Table 2; Bhattacharjee, et al., 2001). Furthermore, we compared the mRNA expression of HMGB family members between 515 LUAD and 59 normal tissues, and 503 LUSC and 52 normal tissues, using UALCAN. HMGB2/3 levels were significantly elevated in both LUAD and LUSC tumor tissues when compared to corresponding normal tissues (Figure 2); data from the Human Protein Atlas followed a similar pattern (Supplementary Figure S1). Of note, although HMGB1 was higher in LUSC tissue than in normal lung tissue, it was lower in LUAD tissue than in lung tissue; this was not the case with HMGB2/3. Based on these analyses, we concluded that HMGB family member genes were over-expressed in patients with NSCLC, with the exception of HMGB1 in LUAD.
FIGURE 1. (A) Schematic diagram of HMGB protein family members. HMGB1/2/3 have box-A/B and an acidic C-terminal tail, while HMGB4 has no acidic C-terminal tail. (B) The number of ONCOMINE analyses in which HMGB family members are significantly differentially expressed between cancer and normal tissue. Student’s t-test was used to compared the transcriptional differences. Cutoffs for p-value and fold change were as following: p-value: 0.05, fold change: 2, gene rank: 10%, data type: mRNA. Red: up-regulation; Blue: down-regulation.
TABLE 2. ONCOMINE analyses indicating significant differences in mRNA expression of HMGB family members between lung cancer and normal lung tissue.
FIGURE 2. The mRNA expression of HMGB family members in NSCLC tissues and normal lung tissues (UALCAN). mRNA expressions of HMGB1/2/3 were elevated in LUSC tissue compared to normal samples (A,C,E). HMGB2/3 was higher in LUAD tissue than in normal lung tissues (D,F), while the expression level of HMGB1 was lower in the LUAD tissue than in normal lung tissue (B). ***p < 0.001, *p < 0.05.
Correlation between mRNA expression of HMGB family members and NSCLC tumor stage
GEPIA was used analyze the relationships between mRNA expression of different HMGB family members and NSCLC tumor stage. The mRNA expression of HMGB2 was significantly related to cancer stages (Figure 3), while the mRNA expression of HMGB1 and HMGB3 were not significantly associated with NSCLC tumor stage. The molecular reason for these differences in association with tumor stage of HMGB2 vs. HMGB1/3 is not yet understood.
FIGURE 3. The relationship between mRNA expression of HMGB family members and tumor stage in patients with NSCLC (GEPIA). The mRNA expression of HMGB2 was significantly related to tumor stage (B), while the mRNA expression of HMGB1/3 was not associated with tumor stage (A,C).
Prognostic features of HMGB family members in lung cancer outcomes
To understand the prognostic features of HMGB family in patients with lung cancer, the Kaplan–Meier Plotter database (2015 version) was used to evaluate overall survival (OS), first-progression survival (FP), and post-progression survival (PPS). In each cohort, patients were divided into low and high-risk group based on the median cutoff value (see Figure 4; Table 3). Higher expression of HMGB1 mRNA were associated with improved OS (p = 0.0059), FP (p = 0.0331), and PPS (p = 5.6E-5). In contrast, higher HMGB3 mRNA expression was significantly associated with poorer OS (p = 0.0082). There was no significant correlation found between HMGB3 expression and FP or PPS in lung cancer. Similarly, there was no significant association observed between HMGB2 expression and OS, FP, or PPS in lung cancer.
FIGURE 4. OS, FP, and PPS survival curves of patients with lung cancer subset by expression of HMGB family members (Kaplan–Meier Plotter). The threshold of significance was a p-value of < 0.05. The association between prognostic value and (A) HMGB1, (B) HMGB2, (C) HMGB3. Red: high expression; Black: low expression. Left: OS; middle: FP; right: PPS.
TABLE 3. The prognostic values of HMGB family members in lung cancer patients (Kaplan–Meier plotter).
Protein and genetic interactions, pathways, and co-expression of HMGB family genes
After analyzing the prognostic value of HMGB family members in NSCLC patients, we used Gene MANIA to analyze protein and genetic interactions, pathways, and co-expression HMGB family members. HMGB family members were found to closely interact with histones, including H1-0, H1-8 and H1-10 (Figure 5A, Supplementary Material S2). In addition, we found that DFFB was the most frequent interacting partner with HMGB family members in this interactive network. DFFB is known to play a vital role in apoptosis (Takeshita et al., 2012; Han et al., 2020). Next, GO and KEGG analysis were performed in g:Profiler to analyze the potential role of HMGB family members and 20 related genes (Figures 5B,C). Chromosome condensation (GO: 0030261) and regulation of chromatin organization (GO: 001902275) were the top two GO Biological Pathways. The most enriched GO Cellular Components of the 16 genes were the nuclear lumen (GO: 003198) and nucleoplasm (GO: 0005654). The three most significantly enriched GO Molecular Functions were nucleosome binding (GO: 0031491), nucleosomal DNA binding (GO: 0031492), and protein-containing complex binding (GO: 0044877). In the KEGG analysis, no significant enrichment was found in relation to the function of HMGB family members in NSCLC. These results suggest that HMGB family members may participate in chromosome condensation, regulation of chromatin organization, and nucleosome binding in NSCLC.
FIGURE 5. Functional involvement of HMGB family members. GeneMANIA was used to predict the function of HMGB family and GO and KEGG analysis were performed to find related genes and pathways. (A) Schematic diagram of protein and genetic interactions, pathways, and co-expression of HMGB family members. (B,C) GO and KEGG analysis details.
HMGB expression correlates with cytotoxic immune cell infiltration in NSCLC
Immune infiltration in the tumor microenvironment is closely related to clinical outcomes in cancers, including NSCLC (Fridman et al., 2011; Chen et al., 2019; Guo et al., 2022). We used the ssGSEA algorithm on the XIANTAO platform (www.xiantao.love) to assess the correlation between expression of HMGB genes and immune infiltration profiles of 24 types of infiltrating immune cells in NSCLC samples. Spearman correlation analysis identified a significant positive correlation between HMGB1 expression levels and levels of infiltrating Th2 cells (R = 0.256, p < 0.001), and a significant negative correlation between HMGB1 levels and infiltrating eosinophils (R = −0.298, p < 0.001), Th17 cells (R = −0.286, p < 0.001), and mast cells (R = −0.262, p < 0.001) (Figure 6A; Table 4). Similarly, HMGB2 expression was positively correlation with infiltrating Th2 cells (R = 0.509, p < 0.001) and was negatively correlated with levels of infiltrating mast cells (R = −0.449, p < 0.001), eosinophils (R = −0.346, p < 0.001), and iDC (R = −0.336, p < 0.001) (Figure 6B; Table 4). In addition, HMGB3 expression was positively associated with infiltrating Th17 cells (R = 0.222, p < 0.001) and had a weak negative correlation with infiltrating neutrophils (R = −0.136, p < 0.001) and iDC (R = −0.135, p < 0.001; Figure 6C, Table 4). These findings suggest that HMGB genes play a role in mediating infiltration of specific types of immune cells in NSCLC, especially Th2 cells, Th17 cells, and mast cells.
FIGURE 6. Correlation of HMGB family members gene expression with levels of immune infiltration in NSCLC. Tcm, central memory T cells; Tem, effector memory T cells; TFH, T follicular helper cells; Tgd, gamma delta T cells; Treg, T regulatory cells.
TABLE 4. The Spearman correlation R value and p-value of HMGB family members with infiltrating immune cells in NSCLC samples.
microRNA–HMGB interactions in the prognosis of NSCLC
MicroRNAs play important roles in cancer initiation and development (Lages et al., 2012; Takasaki, 2015). First, we explored microRNA–HMGBs interaction maps in NSCLC samples from Argonaute CLIP-Seq and degradome-Seq data by using the starBase database (Supplementary Material S3). Then, we searched significant microRNAs in survival in LUAD and LUSC on the database of OncomiR (Supplementary Material S4). Finally, by taking the intersection of the above results (Figure 7), 7, 4, and 26 microRNAs, respectively, were found to potentially target HMGB1, HMGB2, and HMGB3 in NSCLC. We evaluated the correlation of HMGBs expression and microRNAs in NSCLC (Figure 7C, Supplementary Figure S2). We found that hsa-miR-25-3p, hsa-miR-374a-3p, and hsa-miR-93-5p were significantly positively correlated with HMGB1, HMGB2, and HMGB3, respectively. Furthermore, hsa-miR-30a-5p, which has been linked to favorable OS (Wei et al., 2016; Zhang et al., 2017), was predicted to negatively regulate HMGB3 expression, which is consistent with our aforementioned results.
FIGURE 7. Identification of microRNAs of prognostic value in NSCLC and their predicted interactions with HMGB family member genes. (A) Venn diagram depicting the intersection of HMGB-targeting microRNAs in LUAD from the starBase and OncomiR databases. (B) Venn diagram depicting the intersection of HMGB-targeting microRNAs in LUSC from the starBase and OncomiR databases. (C) Correlation of HMGB family members expression and expression of microRNAs that have prognostic value in NSCLC.
Discussion and conclusion
HMGB proteins are chromatin-associated architectural proteins and are fundamental for efficient genome organization and execution of DNA replication, DNA repair, and transcription (Niu et al., 2020; Lou et al., 2021; Wen et al., 2021). The HMGB family consists of four members, including HMGB1-4. HMGB family members are widely involved in the initiation and development of a variety of cancers, including colorectal cancer, hepatocellular carcinoma, and gastric cancer (Fang et al., 2020; Niu et al., 2020; Wen et al., 2021). In this study, we provide a comprehensive analysis of differential expression, immune infiltration and evaluation of the prognostic significance of HMGB1/2/3 in NSCLC.
HMGB1 is a non-histone nuclear protein and is present in almost all eukaryotic cells. HMGB1 works as a nuclear protein and can also be released or secreted into the extracellular interstices to perform a range of biological functions. HMGB1 plays a significant role in many diseases, including autoimmune diseases, infectious diseases, inflammatory diseases, and cancer (Kang et al., 2013; Tripathi et al., 2019). In cancer, HMGB1 plays dual roles as both a pro-cancer and an anti-cancer factor by regulating multiple cellular processes, including proliferation, angiogenesis, and invasion (Wu et al., 2018; Ren et al., 2021). HMGB1 has been shown to enhance tumor growth, angiogenesis, and metastasis of a variety of cancers, including breast cancer, prostate cancer, colorectal cancer, and lung cancer (Zhang J et al., 2014; Zhang L et al., 2014; Meng et al., 2014; Süren et al., 2014; Liu et al., 2015; Sun et al., 2015; Ren et al., 2021). In contrast, Kang et al. (2013) and Luan et al. (2017) reported that HMGB1 acts as a tumor suppressor in pancreatic cancer and endometrial carcinoma (Tang et al., 2010; Zuo et al., 2014). In terms of lung cancer, HMGB1 levels are higher in cancer tissue and in patient serum than in non-cancer tissue, and are correlated with TNM stages and tumor size in NSCLC (Liu et al., 2017; Wang M et al., 2019; Wang et al., 2020). In contrast, some studies have shown that HMGB1 expression is lower in NSCLC tissue compared with matched normal lung tissue, especially in patient with stage III-IV disease (Shen et al., 2009; Wu et al., 2018). In our present study, we found over-expression of HMGB1 only in LUSC vs. normal lung, but under-expressed in LUAD tissue vs. normal lung. These data demonstrate that the upregulation or downregulation of HMGB1 in NSCLC still needs to be further explored, and may be related to pathological tumor type. Previous studies have shown that HMGB1 can promote tumor growth and metastasis through the TAK1-NF-κB or RAGE/JNK/NF-κB pathways (van Beijnum et al., 2013; He et al., 2016; Nguyen et al., 2016; Nguyen et al., 2017; Wu et al., 2018). However, in this study we found that higher mRNA expression of HMGB1 associated with improved OS, FP, and PPS in lung cancer; this is similar to the findings of Dumitriu et al. (2007) and Suzuki et al. (2012). In terms of mechanism, we explored immune infiltration levels correlated with expression of HMGB genes in NSCLC. We identified a strong positive relationship between HMGB1 expression and levels of infiltrating Th2 cells, and a significantly negative correlation between HMGB1 expression and levels of infiltrating eosinophils, Th17 cells, and mast cells. Moreover, we found that HMGB1 was significantly positively correlated with expression of hsa-miR-25-3p, hsa-miR-142-5p, hsa-miR-93-5p, and hsa-let-7g-5p, and was negatively correlated with hsa-miR-375, hsa-miR-181c-5p and hsa-miR-29b-3p. These microRNAs all act as tumor suppressors in NSCLC, according to the OncomiR database. These results suggest that immune infiltration and microRNAs might serve as mechanisms by which HMGB1 expression influences the NSCLC prognosis. These discoveries raise the possibility that HMGB1 may play a significant role in the tumorigenesis of NSCLC, and may be able to be developed as a novel therapeutic target.
High-mobility group box 2 (HMGB2) is generally considered to be an oncogene (Niu et al., 2020; Lou et al., 2021). HMGB2 can promote the proliferation of different cancer cells, including prostate cancer, breast cancer, HCC, and cervical cancer, through the AKT and Wnt pathways (Niu et al., 2020; Wen et al., 2021). In addition, HMGB2 can also promote tumor invasion and migration by regulating noncoding RNAs and CENPU (Li H et al., 2018; Han et al., 2019; Li et al., 2019; Niu et al., 2020). In the context of NSCLC, Lou et al. (2021) reported that HMGB2 protein expression was higher in tumor tissue compared with adjacent tissue, which is consistent with our results in LUAD and LUSC. Lou et al. (2021) also revealed that high tumor mRNA expression of HMGB2 was associated with shorter disease-free survival (DFS) and shorter OS. However, Niu et al. (2020) showed that overexpression of HMGB2 was prognostic of longer OS in cancers including BRCA, CESC, LUSC, READ, STAD, and THCA. Unexpectedly, in our study, HMGB2 expression levels did not correlate with OS, FP, or PPS of lung cancer patients. We found that HMGB2 expression showed a strong positive correlation with levels of infiltrating Th2 cells and negative correlations with the levels of infiltrating mast cells, eosinophils, and iDC. Th2 cells are considered to have tumor-promoting activities, and Th17 cells promote tumor growth by inducing angiogenesis (via IL-17) and by exerting immunosuppressive functions. Mast cells and eosinophils play both anti- and pro-tumorigenic activities, depending on the cellular milieu of the tumor. In addition, we found that HMGB2 expression was significantly positively correlated with expression of hsa-miR-374a-3p, hsa-miR-320a and hsa-miR-362-3p, but negatively correlated with hsa-miR-181c-5p. And these microRNAs are all involved in the prognosis of NSCLC patients. Together with these reports, the role of HMGB2 in lung cancer patients requires the further exploration.
HMGB3 is mainly expressed in embryonic and bone marrow hematopoietic stem cells, and not expressed or is expressed at low levels in normal adult tissues (Wen et al., 2021). HMGB3 acts as an oncogene in the development of human cancers, including liver cancer, breast cancer, gastric cancer, leukemia, and colon cancer (Ren et al., 2020; Wen et al., 2021). HMGB3 can regulate the cell cycle to induce tumor development, promote the proliferation and invasion of cancer cells by regulating the Wnt/β-catenin, MAPK, Akt, HIF-1α, and other signaling pathways, and enhances the activity of tumor stem cell genes to promote the proliferation of cancer cells (Gu et al., 2019a; Gu et al., 2019b; Ren et al., 2020; Wen et al., 2021). However, there are few studies on HMGB3 in lung cancer. HMGB3 was reported to be overexpressed in NSCLC tumor tissue compared to adjacent or normal tissues (Wang J et al., 2019; Li et al., 2020; Ren et al., 2020). Similarly, in our investigation, we found higher mRNA expression of HMGB3 in NSCLC tissues, but mRNA expression of HMGB3 was not correlated with tumor stage of NSCLC. High mRNA expression of HMGB3 was positively related with tumor growth, proliferation, and invasion of NSCLC, suggesting that it may associate with poorer prognosis (Gu et al., 2019a; Gu et al., 2019b; Wang J et al., 2019; Li et al., 2020; Ren et al., 2020). We found that higher levels of HMGB3 were associated with worse OS in NSCLC patients. This might be related to the infiltration of Th17 cells and lower expression of has-miR-30a-5p, as suggested by our results. Based on these data, we suggest that HMGB3 may act as an oncogene in NSCLC and may be prognostic of poor outcomes for patients with NSCLC.
In summary, this study elucidated the prognostic significance, gene interactions, immune infiltration, and co-expression relationships related to the mRNA expression of HMGB family members in NSCLC. Our research revealed that HMGB1 may act as a tumor suppressor gene in NSCLC, which is in contrast with the results of previous studies. We found no significant association between HMGB2 expression and prognostic characteristics of lung cancer. We found that higher expression of HMGB3 associated with poorer outcomes, suggesting that HMGB3 may act as an oncogene in NSCLC. Furthermore, HMGB genes may affect patient survival by altering the immune microenvironment or by interacting with microRNAs. Together, our findings suggest that HMGB family members may act as novel potential therapeutic targets for NSCLC, but further research is needed to better understand the roles of HMGB family members and their potential clinical utility in NSCLC.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving human participants were reviewed and approved by The Zhengzhou University. The patients/participants provided their written informed consent to participate in this study. The animal study was reviewed and approved by Zhengzhou University.
Author contributions
Conceptualization, methodology: XZ, XW, and HG; software, XZ; validation: HG; formal analysis, investigation: XZ and YH; data curation: XZ and XW; writing—original draft preparation, writing—review and editing: XZ and YH; visualization: XZ; supervision, project administration: HG.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.844618/full#supplementary-material
References
Baltayiannis, N., Chandrinos, M., Anagnostopoulos, D., Zarogoulidis, P., Tsakiridis, K., Mpakas, A., et al. (2013). Lung cancer surgery: An up to date. J. Thorac. Dis. 5 (Suppl. 4), S425–S439. doi:10.3978/j.issn.2072-1439.2013.09.17
Bhattacharjee, A., Richards, W. G., Staunton, J., Li, C., Monti, S., Vasa, P., et al. (2001). Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses. Proc. Natl. Acad. Sci. USA 98, 13790–13795. doi:10.1073/pnas.191502998
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68, 394–424. doi:10.3322/caac.21492
Chen, F., Yang, Y., Zhao, Y., Pei, L., and Yan, H. (2019). Immune infiltration profiling in nonsmall cell lung cancer and their clinical significance: Study based on gene expression measurements. DNA Cell Biol. 38, 1387–1401. doi:10.1089/dna.2019.4899
Dumitriu, I. E., Bianchi, M. E., Bacci, M., Manfredi, A. A., and Rovere-Querini, P. (2007). The secretion of HMGB1 is required for the migration of maturing dendritic cells. J. Leukoc. Biol. 81, 84–91. doi:10.1189/jlb.0306171
Fang, J., Ge, X., Xu, W., Xie, J., Qin, Z., Shi, L., et al. (2020). Bioinformatics analysis of the prognosis and biological significance of HMGB1, HMGB2, and HMGB3 in gastric cancer. J. Cell Physiol. 235, 3438–3446. doi:10.1002/jcp.29233
Fridman, W. H., Galon, J., Dieu-Nosjean, M. C., Cremer, I., Fisson, S., Damotte, D., et al. (2011). Immune infiltration in human cancer: Prognostic significance and disease control. Curr. Top. Microbiol. Immunol. 344, 1–24. doi:10.1007/82_2010_46
Gu, J., Xu, T., Huang, Q.-H., Zhang, C.-M., and Chen, H.-Y. (2019). HMGB3 silence inhibits breast cancer cell proliferation and tumor growth by interacting with hypoxia-inducible factor 1α. Cancer Manag. Res. 11, 5075–5089. doi:10.2147/cmar.s204357
Gu, J., Xu, T., Zhang, C. M., Chen, H. Y., Huang, Q. H., and Zhang, Q. (2019). HMGB3 small interfere RNA suppresses mammosphere formation of MDA-MB-231 cells by down-regulating expression of HIF1α. Eur. Rev. Med. Pharmacol. Sci. 23, 9506–9516. doi:10.26355/eurrev_201911_19445
Guo, S., Li, T., Xu, D., Xu, J., Wang, H., Li, J., et al. (2022). Prognostic implications and immune infiltration characteristics of chromosomal instability-related dysregulated CeRNA in lung adenocarcinoma. Front. Mol. Biosci. 9, 843640. doi:10.3389/fmolb.2022.843640
Han, Q., Xu, L., Lin, W., Yao, X., Jiang, M., Zhou, R., et al. (2019). Long noncoding RNA CRCMSL suppresses tumor invasive and metastasis in colorectal carcinoma through nucleocytoplasmic shuttling of HMGB2. Oncogene 38, 3019–3032. doi:10.1038/s41388-018-0614-4
Han, D. S. C., Ni, M., Chan, R. W. Y., Chan, V. W. H., Lui, K. O., Chiu, R. W. K., et al. (2020). The biology of cell-free DNA fragmentation and the roles of DNASE1, DNASE1L3, and DFFB. Am. J. Hum. Genet. 106, 202–214. doi:10.1016/j.ajhg.2020.01.008
He, J., Zhou, M., Chen, X., Yue, D., Yang, L., Qin, G., et al. (2016). Inhibition of SALL4 reduces tumorigenicity involving epithelial-mesenchymal transition via Wnt/β-catenin pathway in esophageal squamous cell carcinoma. J. Exp. Clin. Cancer Res. 35, 98. doi:10.1186/s13046-016-0378-z
Hou, J., Aerts, J., den Hamer, B., van Ijcken, W., den Bakker, M., Riegman, P., et al. (2010). Gene expression based classific ation of non small cell lung carcinomas and survival prediction. PLoS One 5, e10312. doi:10.1371/journal.pone.0010312
Kang, R., Zhang, Q., Zeh, H. J., Lotze, M. T., and Tang, D. (2013). HMGB1 in cancer: Good, bad, or both? Clin. Cancer Res. 19, 4046–4057. doi:10.1158/1078-0432.ccr-13-0495
Lages, E., Ipas, H., Guttin, A., Nesr, H., Berger, F., and Issartel, J. P. (2012). MicroRNAs: Molecular features and role in cancer. Front. Biosci. 17, 2508–2540. doi:10.2741/4068
Landi, M. T., Dracheva, T., Rotunno, M., Figueroa, J. D., Liu, H., Dasgupta, A., et al. (2008). Gene expression signature of cigarette smoking and its role in lung adenocarcinoma development and survival. PLoS One 3, e1651. doi:10.1371/journal.pone.0001651
Li, W., Wang, Q., Feng, Q., Wang, F., Yan, Q., Gao, S.-J., et al. (2019). Oncogenic KSHV-encoded interferon regulatory factor upregulates HMGB2 and CMPK1 expression to promote cell invasion by disrupting a complex lncRNA-OIP5-AS1/miR-218-5p network. PLoS Pathog. 15, e1007578. doi:10.1371/journal.ppat.1007578
Li, Y., Ma, Y., Zhang, T., Feng, C., and Liu, Y. (2020). High-mobility group box 3 (HMGB3) silencing inhibits non-small cell lung cancer development through regulating Wnt/β-catenin pathway. Biol. Chem. 401, 1191–1198. doi:10.1515/hsz-2020-0144
Li H, H., Zhang, H., and Wang, Y. (2018). Centromere protein U facilitates metastasis of ovarian cancer cells by targeting high mobility group box 2 expression. Am. J. Cancer Res. 8, 835–851.
Li S, S., Yang, J., Xia, Y., Fan, Q., and Yang, K.-p. (2018). Long noncoding RNA NEAT1 promotes proliferation and invasion via targeting miR-181a-5p in non-small cell lung cancer. Oncol. Res. 26, 289–296. doi:10.3727/096504017x15009404458675
Liu, W., Zhang, Z., Zhang, Y., Chen, X., Guo, S., Lei, Y., et al. (2015). HMGB1-mediated autophagy modulates sensitivity of colorectal cancer cells to oxaliplatin via MEK/ERK signaling pathway. Cancer Biol. Ther. 16, 511–517. doi:10.1080/15384047.2015.1017691
Liu, P.-L., Liu, W.-L., Chang, J.-M., Chen, Y.-H., Liu, Y.-P., Kuo, H.-F., et al. (2017). MicroRNA-200c inhibits epithelial-mesenchymal transition, invasion, and migration of lung cancer by targeting HMGB1. PLoS One 12, e0180844. doi:10.1371/journal.pone.0180844
Lou, N., Zhu, T., Qin, D., Tian, J., and Liu, J. (2021). High-mobility group box 2 reflects exacerbated disease characteristics and poor prognosis in non-small cell lung cancer patients. Ir. J. Med. Sci. 191, 155. doi:10.1007/s11845-021-02549-8
Luan, X., Ma, C., Wang, P., and Lou, F. (2017). HMGB1 is negatively correlated with the development of endometrial carcinoma and prevents cancer cell invasion and metastasis by inhibiting the process of epithelial-to-mesenchymal transition. Onco Targets Ther. 10, 1389–1402. doi:10.2147/ott.s123085
Meng, Q., Zhao, J., Liu, H., Zhou, G., Zhang, W., Xu, X., et al. (2014). HMGB1 promotes cellular proliferation and invasion, suppresses cellular apoptosis in osteosarcoma. Tumor Biol. 35, 12265–12274. doi:10.1007/s13277-014-2535-3
Nagy, Á., Lánczky, A., Menyhárt, O., and Győrffy, B. (2018). Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets. Sci. Rep. 8, 9227. doi:10.1038/s41598-018-27521-y
Nguyen, A., Bhavsar, S., Riley, E., Caponetti, G., and Agrawal, D. (2016). Clinical value of high mobility group box 1 and the receptor for advanced glycation end-products in head and neck cancer: A systematic review. Int. Arch. Otorhinolaryngol. 20, 382–389. doi:10.1055/s-0036-1583168
Nguyen, A. H., Detty, S. Q., and Agrawal, D. K. (2017). Clinical implications of high-mobility group box-1 (HMGB1) and the receptor for advanced glycation end-products (RAGE) in cutaneous malignancy: A systematic review. Anticancer Res. 37, 1–8. doi:10.21873/anticanres.11282
Niu, L., Yang, W., Duan, L., Wang, X., Li, Y., Xu, C., et al. (2020). Biological functions and theranostic potential of HMGB family members in human cancers. Ther. Adv. Med. Oncol. 12, 1758835920970850. doi:10.1177/1758835920970850
Okayama, H., Kohno, T., Ishii, Y., Shimada, Y., Shiraishi, K., Iwakawa, R., et al. (2020). Identification of genes upregulated in ALK positive and EGFR/KRAS/ALK negative lung adenocarcinomas. Cancer Res. 72, 100–111. doi:10.1158/0008_5472.CAN111403
Qiu, Z., Pan, X. X., and You, D. y. (2020). LncRNA DSCAM-AS1 promotes non-small cell lung cancer progression via regulating miR-577/HMGB1 axis. Neoplasma 67, 871–879. doi:10.4149/neo_2020_190826n821
Relli, V., Trerotola, M., Guerra, E., and Alberti, S. (2019). Abandoning the notion of non-small cell lung cancer. Trends Mol. Med. 25, 585–594. doi:10.1016/j.molmed.2019.04.012
Ren, P., Xing, L., Hong, X., Chang, L., and Zhang, H. (2020). LncRNA PITPNA AS1 boosts the proliferation and migration of lung squamous cell carcinoma cells by recruiting TAF15 to stabilize HMGB3 mRNA. Cancer Med. 9, 7706–7716. doi:10.1002/cam4.3268
Ren, Y., Cao, L., Wang, L., Zheng, S., Zhang, Q., Guo, X., et al. (2021). Autophagic secretion of HMGB1 from cancer-associated fibroblasts promotes metastatic potential of non-small cell lung cancer cells via NFκB signaling. Cell Death Dis. 12, 858. doi:10.1038/s41419-021-04150-4
Rhodes, D. R., Yu, J., Shanker, K., Deshpande, N., Varambally, R., Ghosh, D., et al. (2004). ONCOMINE: A cancer microarray database and integrated data-mining platform. Neoplasia 6, 1–6. doi:10.1016/s1476-5586(04)80047-2
Selamat, S. A., Chung, B. S., Girard, L., Zhang, W., Zhang, Y., Campan, M., et al. (2012). Genome scale analysis of DNA methyla tion in lung adenocarcinoma and integration with mRNA expression. Genome Res. 22, 1197–1211. doi:10.1101/gr.132662.111
Shen, X., Hong, L., Sun, H., Shi, M., and Song, Y. (2009). The expression of high-mobility group protein box 1 correlates with the progression of non-small cell lung cancer. Oncol. Rep. 22, 535–539. doi:10.3892/or_00000468
Song, N., Liu, B., Wu, J.-L., Zhang, R.-F., Duan, L., He, W.-S., et al. (2013). Prognostic value of HMGB3 expression in patients with non-small cell lung cancer. Tumor Biol. 34, 2599–2603. doi:10.1007/s13277-013-0807-y
Song, N., Wang, B., Feng, G., Duan, L., Yuan, S., Jia, W., et al. (2019). Knockdown of high mobility group box 3 impairs cell viability and colony formation but increases apoptosis in A549 human non-small cell lung cancer cells. Oncol. Lett. 17, 2937–2945. doi:10.3892/ol.2019.9927
Su, Y., Bai, M., Zhu, L. P., Jin, Y., Zhang, X. J., Zhou, Q., et al. (2007). Pirh2 shRNA mediated by psiRNA hH1 vector plasmid effectively inhibits the proliferation of lung carcinoma cells: in vitro and in vivo experiments. Zhonghua Yi Xue Za Zhi 87, 1199–1203. doi:10.3760/j:issn:03762491.2007.17.013
Sun, S., Zhang, W., Cui, Z., Chen, Q., Xie, P., Zhou, C., et al. (2015). High mobility group box-1 and its clinical value in breast cancer. Onco Targets Ther. 8, 413–419. doi:10.2147/ott.s73366
Süren, D., Yıldırım, M., Demirpençe, Ö., Kaya, V., Alikanoğlu, A. S., Bülbüller, N., et al. (2014). The role of high mobility group box 1 (HMGB1) in colorectal cancer. Med. Sci. Monit. 20, 530–537. doi:10.12659/MSM.890531
Suzuki, Y., Mimura, K., Yoshimoto, Y., Watanabe, M., Ohkubo, Y., Izawa, S., et al. (2012). Immunogenic tumor cell death induced by chemoradiotherapy in patients with esophageal squamous cell carcinoma. Cancer Res. 72, 3967–3976. doi:10.1158/0008-5472.can-12-0851
Takasaki, S. (2015). Roles of microRNAs in cancers and development. Methods Mol. Biol. 1218, 375–413. doi:10.1007/978-1-4939-1538-5_24
Takeshita, H., Fujihara, J., Ueki, M., Iida, R., Koda, Y., Soejima, M., et al. (2012). Nonsynonymous single-nucleotide polymorphisms of the human apoptosis-related endonuclease - DNA fragmentation factor beta polypeptide, endonuclease G, and flap endonuclease-1 - genes show a low degree of genetic heterogeneity. DNA Cell Biol. 31, 36–42. doi:10.1089/dna.2011.1293
Tang, D., Loze, M. T., Zeh, H. J., and Kang, R. (2010). The redox protein HMGB1 regulates cell death and survival in cancer treatment. Autophagy 6, 1181–1183. doi:10.4161/auto.6.8.13367
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 45, W98–W102. doi:10.1093/nar/gkx247
Taniguchi, N., Kawakami, Y., Maruyama, I., and Lotz, M. (2018). HMGB proteins and arthritis. Hum. Cell 31, 1–9. doi:10.1007/s13577-017-0182-x
Tripathi, A., Shrinet, K., and Kumar, A. (2019). HMGB1 protein as a novel target for cancer. Toxicol. Rep. 6, 253–261. doi:10.1016/j.toxrep.2019.03.002
van Beijnum, J. R., Nowak-Sliwinska, P., van den Boezem, E., Hautvast, P., Buurman, W. A., and Griffioen, A. W. (2013). Tumor angiogenesis is enforced by autocrine regulation of high-mobility group box 1. Oncogene 32, 363–374. doi:10.1038/onc.2012.49
Wachi, S., Yoneda, K., and Wu, R. (2005). Interactome transcriptome analysis reveals the high centrality of genes differ entially expressed in lung cancer tissues. Bioinformatics 21, 4205–4208. doi:10.1093/bioinformatics/bti688
Wang, J.-L., Wu, D.-W., Cheng, Z.-Z., Han, W.-Z., Xu, S.-W., and Sun, N.-N. (2014). Expression of high mobility group box - B1 (HMGB-1) and matrix metalloproteinase-9 (MMP-9) in non-small cell lung cancer (NSCLC). Asian Pac. J. Cancer Prev. 15, 4865–4869. doi:10.7314/apjcp.2014.15.12.4865
Wang, X. H., Zhang, S. Y., Shi, M., and Xu, X. P. (2020). HMGB1 promotes the proliferation and metastasis of lung cancer by activating the wnt/β-catenin pathway. Technol. Cancer Res. Treat. 19, 1533033820948054. doi:10.1177/1533033820948054
Wang J, J., Sheng, Z., and Cai, Y. (2019). Effects of microRNA 513b on cell proliferation, apoptosis, invasion, and migration by targeting HMGB3 through regulation of mTOR signaling pathway in non small cell lung cancer. J. Cell. Physiology 234, 10934–10941. doi:10.1002/jcp.27921
Wang M, M., Gauthier, A., Daley, L., Dial, K., Wu, J., Woo, J., et al. (2019). The role of HMGB1, a nuclear damage-associated molecular pattern molecule, in the pathogenesis of lung diseases. Antioxidants Redox Signal. 31, 954–993. doi:10.1089/ars.2019.7818
Wang Y, Y., Xu, R., Zhang, D., Lu, T., Yu, W., Wo, Y., et al. (2019). Circ-ZKSCAN1 regulates FAM83A expression and inactivates MAPK signaling by targeting miR-330-5p to promote non-small cell lung cancer progression. Transl. Lung Cancer Res. 8, 862–875. doi:10.21037/tlcr.2019.11.04
Wei, W., Yang, Y., Cai, J., Cui, K., Li, R. x., Wang, H., et al. (2016). MiR-30a-5p suppresses tumor metastasis of human colorectal cancer by targeting ITGB3. Cell Physiol. Biochem. 39, 1165–1176. doi:10.1159/000447823
Wen, B., Wei, Y.-t., and Zhao, K. (2021). The role of high mobility group protein B3 (HMGB3) in tumor proliferation and drug resistance. Mol. Cell Biochem. 476, 1729–1739. doi:10.1007/s11010-020-04015-y
Wu, X. J., Chen, Y. Y., Gong, C. C., and Pei, D. S. (2018). The role of high mobility group protein box 1 in lung cancer. J. Cell. Biochem. 119, 6354–6365. doi:10.1002/jcb.26837
Xie, Y., Wang, L., and Yang, D. (2020). CircEPSTI1 promotes the progression of non-small cell lung cancer through miR-145/HMGB3 Axis. Cancer Manag. Res. 12, 6827–6836. doi:10.2147/cmar.s252893
Zhang, S., Liu, Q., Zhang, Q., and Liu, L. (2017). MicroRNA-30a-5p suppresses proliferation, invasion and tumor growth of hepatocellular cancer cells via targeting FOXA1. Oncol. Lett. 14, 5018–5026. doi:10.3892/ol.2017.6745
Zhang J, J., Kou, Y.-B., Zhu, J.-S., Chen, W.-X., and Li, S. (2014). Knockdown of HMGB1 inhibits growth and invasion of gastric cancer cells through the NF-κB pathway in vitro and in vivo. Int. J. Oncol. 44, 1268–1276. doi:10.3892/ijo.2014.2285
Zhang L, L., Han, J., Wu, H., Liang, X., Zhang, J., Li, J., et al. (2014). The association of HMGB1 expression with clinicopathological significance and prognosis in hepatocellular carcinoma: A meta-analysis and literature review. PLoS One 9, e110626. doi:10.1371/journal.pone.0110626
Keywords: non-small cell lung cancer, HMGB family, bioinformatics analysis, prognostic value, biomaker
Citation: Zheng X, Wang X, He Y and Ge H (2022) Systematic analysis of expression profiles of HMGB family members for prognostic application in non-small cell lung cancer. Front. Mol. Biosci. 9:844618. doi: 10.3389/fmolb.2022.844618
Received: 28 December 2021; Accepted: 27 June 2022;
Published: 18 July 2022.
Edited by:
Stella Liberman-Aronov, Ariel University, IsraelReviewed by:
Qing Lin, Johns Hopkins Medicine, United StatesLarance Ronsard, Ragon Institute, United States
Copyright © 2022 Zheng, Wang, He and Ge. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hong Ge, gehong_ro@163.com