Screening and identification of angiogenesis-related genes as potential novel prognostic biomarkers of hepatocellular carcinoma through bioinformatics analysis

Prognostic Hypoxia-Angiogenesis-Related Gene Signature in Hepatocellular Carcinoma, in Which HILPDA Contributes to Tumor Progression

Objective

Hepatocellular carcinoma (HCC) is the predominant form of liver cancer. Hypoxia can be involved in HCC tumor growth, invasion and metastasis through inducing angiogenesis. Nevertheless, the assessment of the impact of hypoxia and angiogenesis on the prognosis of HCC remains inadequate.

Methods

According to hypoxia-angiogenesis-related genes (HARGs) expression information and clinical data from patients within the Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA-LIHC) cohort, we constructed a prognostic model (HARG-score) using bioinformatic tools. In addition to assessing the predictive ability of this prognostic model in both Liver Cancer-Riken-Japan (LIRI-JP) and GSE14520 cohorts, we analyzed the correlation between HARG-score and clinical characteristics, immune infiltration and immunotherapy efficacy. Moreover, we investigated the exact role and underlying mechanism of key HARGs through molecular experiments.

Results

We constructed a 5-gene prognostic model HARG-score consisting of hypoxia-inducible lipid droplet-associated (HILPDA), erythropoietin (EPO), solute carrier family 2 member 1 (SLC2A1), proteasome subunit alpha type 7 (PSMA7) and cAMP responsive element-binding protein 1 (CREB1) through differentially expressed HARGs. The findings demonstrated that HARG-score was a good predictor of the prognosis of HCC patients from distinct cohorts and was correlated with clinical characteristics and immune infiltration. Furthermore, the HARG-score was identified as an independent prognostic factor. Lower HARG-score implied greater immunotherapy efficacy and better response. The expression and prognostic significance of these 5 genes were additionally validated in clinical data. In addition, experimental data revealed that the key gene HILPDA contributes to the progression of HCC through facilitating angiogenesis and affecting the expression of cytotoxic T-lymphocyte-associated protein 4 (CTLA4).

Conclusion

HARG-score has promising applications in prognosis prediction of HCC patients, in which HILPDA may be a latent prognostic biomarker and therapeutic target, providing a foundation for further research and treatment of HCC.

MRI Radiomics Combined with Clinicopathological Factors for Predicting 3-Year Overall Survival of Hepatocellular Carcinoma After Hepatectomy

Background

A limited number of studies have examined the use of radiomics to predict 3-year overall survival (OS) after hepatectomy in patients with hepatocellular carcinoma (HCC). This study develops 3-year OS prediction models for HCC patients after liver resection using MRI radiomics and clinicopathological factors.

Materials and Methods

A retrospective analysis of 141 patients who underwent surgical resection of HCC was performed. Patients were randomized into two set: the training set (n=98) and the validation set (n=43) including the survival groups (n=111) and non-survival groups (n=30) based on 3-year survival after hepatectomy. Furthermore, x2 or Fisher's exact test, univariate and multivariate logistic regression analyses were conducted to determine independent clinicopathological risk factors associated with 3-year OS. 1688 quantitative imaging features were extracted from preoperative T2-weighted imaging (T2WI) and contrast-enhanced magnetic resonance imaging (CE-MRI) of arterial phase (AP), portal venous phases (PVP)and delay period (DP). The features were selected using the variance threshold method, the select K best method and the least absolute shrinkage and selection operator (LASSO) algorithm. By using Bernoulli Naive Bayes (BernoulliNB) and Multinomial Naive Bayes (MultinomialNB) classifiers, we constructed models based on the independent clinicopathological factors and Rad-scores. To determine the best model, receiver operating characteristics (ROC) and Delong's test were used. Moreover, calibration curves were used to determine the calibration ability of the model, while decision curve analysis (DCA) was implemented to evaluate its clinical benefit.

Results

The fusion model showed excellent prediction precision with AUC of 0.910 and 0.846 in training and validation set and revealed significant diagnostic accuracy and value in the calibration curve and DCA analysis.

Conclusion

Nomograms based on MRI radiomics and clinicopathological factors have significant predictive value for 3-year OS after hepatectomy and can be used for risk classification.

Identification of VEGFs-related gene signature for predicting microangiogenesis and hepatocellular carcinoma prognosis

Microangiogenesis is an important prognostic factor in various cancers, including hepatocellular carcinoma (HCC). The Vascular Endothelial Growth Factor (VEGF) has been shown to contribute to tumor angiogenesis. Recently, several studies have investigated the regulation of VEGF production by a single gene, with few researchers exploring all genes that affect VEGF production. In this study, we comprehensively analyzed all genes affecting VEGF production in HCC and developed a risk model and gene-based risk score based on VEGF production. Moreover, the model's predictive capacity on prognosis of HCCs was verified using training and validation datasets. The developed model showed good prediction of the overall survival rate. Patients with a higher risk score experienced poor outcomes compared to those with a lower risk score. Furthermore, we identified the immunological causes of the poor prognosis of patients with high-risk scores comparing with those with low-risk scores.

A novel angiogenesis-associated risk score predicts prognosis and characterizes the tumor microenvironment in colon cancer

Background

Angiogenesis of the tumor microenvironment (TME) can promote the proliferation and metastases of colon cancer (CC). However, there is a lack of bioinformatics analysis to comprehensively clarify the molecular characteristics, immune interaction characteristics and predictive values of angiogenesis characteristics in CC patients. This study aimed to perform a comprehensive elucidation of the correlation between angiogenesis and CC for the purpose of improving the clinical management of CC.

Methods

Angiogenesis-associated genes (AAGs) were evaluated in the population of CC patients from the Cancer Genome Atlas database and Gene Expression Omnibus dataset. The expression, prognostic role, and immune cell infiltration of AAGs were assessed first. And then we established the AAGs score to further explore the prognosis and treatment response of angiogenesis characteristics in individual patient.

Results

Totally, we identified two different molecular subtypes of angiogenesis, and there was a significant difference in the background of genome, expression profiles, prognosis, and characteristics of TME between two subtypes. And the AAGs score was independently associated with over survival in CC patients, the prognostic value was significant and confirmed in the entire cohort. And we also constructed a nomogram based on the risk score and clinical parameters to maximize the predictive ability of the risk score. Additionally, the AAGs score was significantly correlated with the tumor mutation burden score, cancer stem cell score and drug sensitivity.

Conclusions

Our study elucidated the role of angiogenesis characteristics in CC and the AAGs score could help clinicians plan for individual management with chemotherapy agents and promote the development of immunotherapy in CC. Prospective studies need to be conducted to further confirm our findings.

Crosstalk of angiogenesis-related subtypes, establishment of a prognostic signature and immune infiltration characteristics in colorectal adenocarcinoma

Background

Colorectal adenocarcinoma (COAD) is one of the most common malignancies and angiogenesis is vital to the development of cancer. Here, we explored the roles of angiogenesis-related genes (ARGs) that affect the prognosis of COAD and constructed risk models to assess patient prognosis, immune characteristics, and treatment outcomes.

Methods

We comprehensively characterized the transcriptional and genetic modifications of 48 ARGs in COAD and evaluated the expression patterns. We identified two ARG subgroups using the consensus clustering algorithm. Based on the differentially expressed genes (DEGs) of two ARG subtypes, we calculated risk score, namely ARG_scores, and calssified COAD patients into different risk groups. To investigate the expression of ARG_score-related genes, qRT-PCR was performed. Subsequently, we mapped the nomogram to visually and accurately describe the value of the application of ARG_score. Finally, the correlation between ARG_score and clinical features, immune infiltration along with drug sensitivity were explored.

Results

We identified two ARG related subgroups and there were great differences in overall survival (OS) and tumor microenvironment. Then, we created an ARG_score for predicting overall survival based on eight DEGs and confirmed its reliable predictive power in COAD patients, with higher ARG_score associated with worse prognosis. Furthermore, eight ARG_score-related genes expression was investigated by qRT-PCR. To make the ARG_score clinically feasible, we created a highly reliable nomogram. We also found a higher proportion of microsatellite instability-high (MSI-H) and higher tumor mutational burden (TMB) in the high-risk group. In addition, ARG_score was notably correlated with cancer stem cell indices and drug sensitivity.

Conclusion

This scoring model has potential clinical application value in the prognosis, immune microenvironment and therapeutic drug sensitivity of COAD, which provides new insights for personalized treatment.

Construction and validation of an angiogenesis-related scoring model to predict prognosis, tumor immune microenvironment and therapeutic response in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors in the world with high morbidity and mortality. Identifying an effective marker for predicting the prognosis and therapeutic response is extremely meaningful. Angiogenesis-related genes (ARGs) play important roles in the tumor progression and immune-suppressive microenvironment formation.

Methods

The differential expressed ARGs associated with the prognosis of HCC were identified in the TCGA dataset. Univariate Cox and least absolute shrinkage selection operator (LASSO) regression were applied to construct a ARGs Scoring model. The prognostic value of the ARGs Scoring model was assessed by Cox regression, Kaplan-Meier (KM) and ROC curve analyses. Then the model was further validated in an external dataset, ICGC dataset. The patients were split into two groups based on the ARGs Score and the clinical features were compared. TIMER, CIBERSORT and xCell algorithms were utilized to analyze the correlation between the ARGs Score and tumor immune microenvironment (TIME). Furthermore, we analyzed the efficacy of the model in predicting the therapeutic response for immunotherapy, targeted therapy and TACE treatment in different cohorts.

Results

A total of 97 differential expressed ARGs were identified relating to the prognosis of HCC patients from the TCGA dataset. Then the ARGs Scoring model based on a 9-gene signature was constructed using the Cox and LASSO regression analyses. Higher ARGs Score had a poor clinical outcome and was considered to be an independent prognostic predictor for HCC in the multivariate Cox analysis. The ARGs Score was related to the enrichment of various immune cells, such as CD4+ T cells, Treg, macrophage, neutrophil and dendritic cells, exhibiting a more immunosuppressive phenotype. Higher ARGs Score was correlated with higher expression of immune checkpoint genes and poor response to immunotherapy. Furthermore, higher ARGs Score indicated poor therapeutic response in the sorafenib and TACE treatment cohorts, individually.

Conclusions

The ARGs Scoring model exhibited robust predictive value for the prognosis and TIME for HCC patients. Higher ARGs Score indicated poor therapeutic response of the immunotherapy, sorafenib and TACE treatment. The ARGs Scoring model could be used as a biomarker to help physicians to develop more individualized treatment for HCC patients.

Dissecting a hypoxia-related angiogenic gene signature for predicting prognosis and immune status in hepatocellular carcinoma

Background

Hypoxia and angiogenesis, as prominent characteristics of malignant tumors, are implicated in the progression of hepatocellular carcinoma (HCC). However, the role of hypoxia in the angiogenesis of liver cancer is unclear. Therefore, we explored the regulatory mechanisms of hypoxia-related angiogenic genes (HRAGs) and the relationship between these genes and the prognosis of HCC.

Methods

The transcriptomic and clinical data of HCC samples were downloaded from public datasets, followed by identification of hypoxia- and angiogenesis-related genes in the database. A gene signature model was constructed based on univariate and multivariate Cox regression analyses, and validated in independent cohorts. Kaplan-Meier survival and time-dependent receiver operating characteristic (ROC) curves were generated to evaluate the model’s predictive capability. Gene set enrichment analysis (GSEA) was performed to explore signaling pathways regulated by the gene signature. Furthermore, the relationships among gene signature, immune status, and response to anti-angiogenesis agents and immune checkpoint blockade (ICB) were analyzed.

Results

The prognostic model was based on three HRAGs (ANGPT2, SERPINE1 and SPP1). The model accurately predicted that low-risk patients would have longer overall survival than high-risk patients, consistent with findings in other cohorts. GSEA indicated that high-risk group membership was significantly associated with hypoxia, angiogenesis, the epithelial-mesenchymal transition, and activity in immune-related pathways. The high-risk group also had more immunosuppressive cells and higher expression of immune checkpoints such as PD-1 and PD-L1. Conversely, the low-risk group had a better response to anti-angiogenesis and ICB therapy.

Conclusions

The gene signature based on HRAGs was predictive of prognosis and provided an immunological perspective that will facilitate the development of personalized therapies.

Dissecting the Role of Immune Checkpoint Regulation Patterns in Tumor Microenvironment and Prognosis of Gastric Cancer

Many studies suggest that immune checkpoint molecules play a vital role in tumor progression and immune responses. However, the impact of the comprehensive regulation pattern of immune checkpoint molecules on immune responses, tumor microenvironment (TME) formation, and patient prognosis is poorly understood. In this study, we evaluated immune checkpoint regulation patterns in 1,174 gastric cancer (GC) samples based on 31 immune checkpoint genes (ICGs). Three distinct immune checkpoint regulation patterns with significant prognostic differences were ultimately identified. Moreover, GC patients were divided into two subgroups according to immune checkpoint score (ICscore). Patients with lower ICscore were characterized by a favorable prognosis and enhanced immune infiltration as well as an increased tumor mutation burden, non-recurrence, and microsatellite instability-high. Collectively, this study indicated that immune checkpoint regulation patterns were essential to forming the diversity of TME and a better understanding of that will contribute to assessing the characteristics of TME in GC, which intends to improve the development of immunotherapy.

Downregulation of Low-density lipoprotein receptor-related protein 1B (LRP1B) inhibits the progression of hepatocellular carcinoma cells by activating the endoplasmic reticulum stress signaling pathway

Hepatocellular carcinoma (HCC) has a high recurrence rate and mortality rate even after surgery. Low-density lipoprotein receptor-related protein 1B (LRP1B) has proven to be involved in tumor development and progression of multiple malignancies. However, the function of LRP1B in HCC progression has not been fully elucidated. Thus, we conducted this study to explore the relationship between LRP1B and HCC. Bioinformatic analyses implied that LRP1B was highly expressed in HCC tissues. High LRP1B expression was shown to be related to poor outcomes and the determination of HCC patients’ tumor stage. LRP1B deletion impeded the proliferation, migration, and invasion of HCC cells. Further investigation demonstrated that silencing LRP1B expression enhanced the sensitivity of HCC cells to doxorubicin. LRP1B deletion inhibited HCC progression by regulating the PERK-ATF4-CHOP signaling pathway. Additionally, we probed the genomic alterations of LRP1B in HCC and the impact on the prognosis of patients. Collectively, our results suggest that LRP1B plays an essential role in the promotion of HCC progression by regulating the PERK-ATF4-CHOP signaling pathway, which is a potential prognostic biomarker and a promising therapeutic target of HCC.

Pyroptosis-Related Risk Signature Exhibits Distinct Prognostic, Immune, and Therapeutic Landscapes in Hepatocellular Carcinoma

Background: Hepatocellular carcinoma (HCC) is a common abdominal cancer. The existing therapeutic approaches often fail to achieve satisfactory results. Pyroptosis, an inflammatory form of programmed cell death, provides new ideas for anticancer treatment. However, the roles of pyroptosis-related (PR) genes (PRGs) in HCC remain elusive. Methods: Differentially expressed genes (DEGs) (n = 22) were screened out using TCGA and GTEx databases. A novel PR risk signature was constructed through Lasso regression analysis. Its prognostic value was evaluated through a series of survival analyses and was tested in ICGC and GSE14520 cohorts. CIBERSORT, ssGSEA, and ESTIMATE methods were employed to determine the effects of the PR risk score on the tumor immune microenvironment (TIM). The TIDE scoring system, IMvigor210 cohort, GSE109211 dataset, and GSDC database were applied to explore the associations of the PR risk score with therapeutic effects. The biofunctions of WNK1 in hepatocellular cancer (HC) cells were confirmed through qPCR, colony formation, and Transwell assays. Results: Overall, 22 of 45 PRGs (48.9%) were abnormally expressed in HCC samples. Then, a PR risk signature consisting of eight PRGs was constructed. A high PR risk score led to an unfavorable prognosis. The PR risk score was identified as an independent prognostic factor of HCC and could increase the decision-making benefit of the traditional TNM model. In addition, we established a nomogram containing the clinical stage and PR risk score to predict the survival rates of HCC patients. The prognostic value of the PR model was successfully validated in ICGC and GSE14520 cohorts. Moreover, high PR risk conferred the decreased infiltration level of CD8⁺ T cells and weakened the activities of "cytolytic activity" pathways. As for therapeutic correlation, a high PR risk score seemed to imply a poor efficacy of PD-1/L1 inhibitors and sorafenib. Finally, the overexpression of WNK1 could promote the proliferation, migration, and invasion of HC cells. Conclusions: The PR risk score was closely related to the prognosis, antitumor immune process, therapeutic outcomes, and malignant progression of HCC. WNK1, the core regulator of pyroptosis, possesses pro-oncogenic abilities, showing promise as a novel treatment target.