Immune system‑associated genes increase malignant progression and can be used to predict clinical outcome in patients with hepatocellular carcinoma

Integrative single-cell analysis of LUAD: elucidating immune cell dynamics and prognostic modeling based on exhausted CD8+ T cells

Background

The tumor microenvironment (TME) plays a pivotal role in the progression and metastasis of lung adenocarcinoma (LUAD). However, the detailed characteristics of LUAD and its associated microenvironment are yet to be extensively explored. This study aims to delineate a comprehensive profile of the immune cells within the LUAD microenvironment, including CD8+ T cells, CD4+ T cells, and myeloid cells. Subsequently, based on marker genes of exhausted CD8+ T cells, we aim to establish a prognostic model for LUAD.

Method

Utilizing the Seurat and Scanpy packages, we successfully constructed an immune microenvironment atlas for LUAD. The Monocle3 and PAGA algorithms were employed for pseudotime analysis, pySCENIC for transcription factor analysis, and CellChat for analyzing intercellular communication. Following this, a prognostic model for LUAD was developed, based on the marker genes of exhausted CD8+ T cells, enabling effective risk stratification in LUAD patients. Our study included a thorough analysis to identify differences in TME, mutation landscape, and enrichment across varying risk groups. Moreover, by integrating risk scores with clinical features, we developed a new nomogram. The expression of model genes was validated via RT-PCR, and a series of cellular experiments were conducted, elucidating the potential oncogenic mechanisms of GALNT2.

Results

Our study developed a single-cell atlas for LUAD from scRNA-seq data of 19 patients, examining crucial immune cells in LUAD's microenvironment. We underscored pDCs' role in antigen processing and established a Cox regression model based on CD8_Tex-LAYN genes for risk assessment. Additionally, we contrasted prognosis and tumor environments across risk groups, constructed a new nomogram integrating clinical features, validated the expression of model genes via RT-PCR, and confirmed GALNT2's function in LUAD through cellular experiments, thereby enhancing our understanding and approach to LUAD treatment.

Conclusion

The creation of a LUAD single-cell atlas in our study offered new insights into its tumor microenvironment and immune cell interactions, highlighting the importance of key genes associated with exhausted CD8+ T cells. These discoveries have enabled the development of an effective prognostic model for LUAD and identified GALNT2 as a potential therapeutic target, significantly contributing to the improvement of LUAD diagnosis and treatment strategies.

The human factor H protein family - an update

Complement is an ancient and complex network of the immune system and, as such, it plays vital physiological roles, but it is also involved in numerous pathological processes. The proper regulation of the complement system is important to allow its sufficient and targeted activity without deleterious side-effects. Factor H is a major complement regulator, and together with its splice variant factor H-like protein 1 and the five human factor H-related (FHR) proteins, they have been linked to various diseases. The role of factor H in inhibiting complement activation is well studied, but the function of the FHRs is less characterized. Current evidence supports the main role of the FHRs as enhancers of complement activation and opsonization, i.e., counter-balancing the inhibitory effect of factor H. FHRs emerge as soluble pattern recognition molecules and positive regulators of the complement system. In addition, factor H and some of the FHR proteins were shown to modulate the activity of immune cells, a non-canonical function outside the complement cascade. Recent efforts have intensified to study factor H and the FHRs and develop new tools for the distinction, quantification and functional characterization of members of this protein family. Here, we provide an update and overview on the versatile roles of factor H family proteins, what we know about their biological functions in healthy conditions and in diseases.

Development and validation of web-based nomograms for predicting survival status in patients with intrahepatic cholangiocarcinoma depending on the surgical status: a SEER database analysis

This study aimed to develop and validate prognostic nomograms that can estimate the probability of 1-, 3- and 5-year overall survival (OS) as well as cancer-specific survival (CSS) for Intrahepatic cholangiocarcinoma (ICCA) patients. Clinical data of 1446 patients diagnosed with ICCA between 2010 and 2017 from the Surveillance, Epidemiology, and End Results (SEER) database were analyzed. In both the OS and the CSS group, the training cohort and validation cohort were divided into a 7:3 ratio. Age, sex, AJCC T stage, AJCC N stage, AJCC M stage, surgical status, and tumor grade were selected as independent prognostic risk factors to build the nomograms. To compare the efficacy of predicting 1-, 3-, and 5-year OS and CSS rates of the nomogram with the 8th edition of the American Joint Committee on Cancer (AJCC) staging system, we evaluated the Harrell's index of concordance (C-index), area under the receiver operating characteristic curve (AUC) and decision curve analysis (DCA) in both cohorts. The results showed the nomogram for 1-, 3-, and 5-year OS and CSS prediction performed better than the AJCC staging system. In the subgroup analysis for patients could not receive surgery as the primary treatment. We developed two nomograms for predicting the 1-, and 2-year OS and CSS rates following the same analysis procedure. Results indicate that the performance of both nomograms, which contained sex, AJCC T stage, AJCC M stage, chemotherapy, and tumor grade and prognostic factors, was also superior to the AJCC staging system. Meanwhile, four dynamic network-based nomograms were published. The survival analysis showed the survival rate of patients classified as high-risk based on the nomogram score was significantly lower compared to those categorized as low-risk (P < 0.0001). Finally, accurate and convenient nomograms were established to assist clinicians in making more personalized prognosis predictions for ICCA patients.

Identification of a new gene signature for prognostic evaluation in cervical cancer: based on cuproptosis-associated angiogenesis and multi-omics analysis

Patients with recurrent or metastatic cervical cancer are in urgent need of novel prognosis assessment or treatment approaches. In this study, a novel prognostic gene signature was discovered by utilizing cuproptosis-related angiogenesis (CuRA) gene scores obtained through weighted gene co-expression network analysis (WGCNA) of The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets. To enhance its reliability, the gene signature was refined by integrating supplementary clinical variables and subjected to cross-validation. Meanwhile, the activation of the VEGF pathway was inferred from an analysis of cell-to-cell communication, based on the expression of ligands and receptors in cell transcriptomic datasets. High-CuRA patients had less infiltration of CD8 + T cells and reduced expression of most of immune checkpoint genes, which indicated greater difficulty in immunotherapy. Lower IC50 values of imatinib, pazopanib, and sorafenib in the high-CuRA group revealed the potential value of these drugs. Finally, we verified an independent prognostic gene SFT2D1 was highly expressed in cervical cancer and positively correlated with the microvascular density. Knockdown of SFT2D1 significantly inhibited ability of the proliferation, migration, and invasive in cervical cancer cells. CuRA gene signature provided valuable insights into the prediction of prognosis and immune microenvironment of cervical cancer, which could help develop new strategies for individualized precision therapy for cervical cancer patients.

Comprehensive analysis of necroptotic patterns and associated immune landscapes in individualized treatment of skin cutaneous melanoma

Skin cutaneous melanoma (SKCM) constitutes a malignant cutaneous neoplasm characterized by an exceedingly unfavorable prognosis. Over the past years, necroptosis, a manifestation of inflammatory programmed cell demise, has gained substantial traction in its application. However, a conclusive correlation between the expression of necroptosis-related genes (NRGs) and SKCM patient's prognosis remains elusive. In this endeavor, we have undertaken an integrative analysis of genomic data, aiming to provide an exhaustive evaluation of the intricate interplay between melanoma necroptosis and immune-infiltration nuances within the tumor microenvironment. Through meticulous scrutiny, we have endeavored to discern the prognostic potency harbored by individual necroptosis-associated genes. Our efforts culminated in the establishment of a risk stratification framework, allowing for the appraisal of necroptosis irregularities within each afflicted cutaneous melanoma patient. Notably, those SKCM patients classified within the low-risk cohort exhibited a markedly elevated survival quotient, in stark contrast to their high-risk counterparts (p < 0.001). Remarkably, the low-risk cohort not only displayed a more favorable survival rate but also exhibited an enhanced responsiveness to immunotherapeutic interventions, relative to their high-risk counterparts. The outcomes of this investigation proffer insights into a conceivable mechanistic underpinning linking necroptosis-related attributes to the intricacies of the tumor microenvironment. This prompts a conjecture regarding the plausible association between necroptosis characteristics and the broader tumor microenvironmental milieu. However, it is imperative to emphasize that the pursuit of discerning whether the expression profiles of NRG genes can indeed be regarded as viable therapeutic targets necessitates further comprehensive exploration and scrutiny. In conclusion, our study sheds light on the intricate interrelationship between necroptosis-related factors and the tumor microenvironment, potentially opening avenues for therapeutic interventions. However, the prospect of translating these findings into clinical applications mandates rigorous investigation.

CGB5, INHBA and TRAJ19 Hold Prognostic Potential as Immune Genes for Patients with Gastric Cancer

BACKGROUND

Gastric cancer (GC) seriously threatens people's health and life quality worldwide.

AIM

The current study sought to explore prognostic immune genes and their regulatory network in GC.

METHODS

First, expression data in GC and normal samples were analyzed based on bioinformatics analysis. Immune-related genes were identified and confirmed with univariate/multivariate Cox analysis and receiver-operating characteristic curve. The upstream transcription factors of immune genes were subsequently predicted, and their regulatory network was constructed. GC and adjacent normal tissues were obtained from 76 patients with GC to determine the expression patterns of immune genes and their correlation with overall prognosis. CD8⁺ T-cell infiltration of patients with high or low risk was detected by means of immunohistochemistry.

RESULTS

Bioinformatics analysis highlighted 3689 differentially expressed genes in GC, including 87 immune genes, 8 of which were significantly associated with patient survival. CGB5 and INHBA were high-risk genes, while TRAJ19 was identified as a low-risk gene, all of which were found to be regulated by 11 different transcription factors. Furthermore, CGB5 and INHBA exhibited negative correlation with the prognosis of GC patients; however, TRAJ19 was positively correlated with GC patient prognosis. The incidence of lymph node metastasis was higher, the pathological stage was advanced and the infiltrated CD8⁺ T cells were fewer in the high-risk GC group.

CONCLUSIONS

Overall, our findings identified the key roles of CGB5, INHBA and TRAJ19 in prognosis GC patients, serving as an important gene set for prognostic prediction.

Identification of a Solute Carrier Family-Based Signature for Predicting Overall Survival in Osteosarcoma

Given the important role of SLC family in essential physiological processes including nutrient uptake, ion transport, and waste removal, and that their dysregulation was found in distinct forms of cancer, here we identified a novel gene signature of SLC family for patient risk stratification in osteosarcoma. Gene expression data and relevant clinical materials of osteosarcoma samples were retrieved from The Cancer Genome Atlas (TCGA) database. Prognosis-related SLC genes were identified by performing univariate Cox regression analysis and were utilized to construct a four-SLC gene signature in osteosarcoma. It allowed patients to be classified into high- and low-risk groups, and Kaplan-Meier survival analysis in the training, testing, entire, and external GSE21257 cohorts suggested that the overall survival of patients in high-risk group was consistently worse than that in low-risk group, suggesting the promising accuracy and generalizability of the SLC-based signature in predicting the prognosis of patients with osteosarcoma. Moreover, univariate and multivariate Cox regression analyses indicated that the derived risk score was the only independent prognostic factor for osteosarcoma patients in TCGA and GSE21257 cohorts. Besides, a prognostic nomogram comprising the derived risk score and clinical features including gender and age was developed for clinical decision-making. Functional enrichment analyses of the differentially expressed genes between high- and low-risk group revealed that immune-related biological processes and pathways were significantly enriched. Estimation of tumor immune microenvironment using ESTIMATE algorithm revealed that patients with lower risk score had higher stromal, immune, and ESTIMATE score, and lower tumor purity. ssGSEA analyses indicated that the scores of various immune subpopulations including CD8+ T cells, DCs, and TIL were lower in high-risk group than these in low-risk group in both cohorts. As for the related immune functions, the scores of APC co-inhibition, CCR, check-point, T cell co-stimulation, and Type II IFN response were lower in high-risk group than these in low-risk group in both cohorts. In all, we identified a novel prognostic signature based on four SLC family genes that accurately predicted overall survival in osteosarcoma patients. Furthermore, the signature is linked to differences in immunological status and immune cell infiltrations in the tumor microenvironment.

DNA Methylation-Driven Genes for Developing Survival Nomogram for Low-Grade Glioma

Low-grade gliomas (LGG) are heterogeneous, and the current predictive models for LGG are either unsatisfactory or not user-friendly. The objective of this study was to establish a nomogram based on methylation-driven genes, combined with clinicopathological parameters for predicting prognosis in LGG. Differential expression, methylation correlation, and survival analysis were performed in 516 LGG patients using RNA and methylation sequencing data, with accompanying clinicopathological parameters from The Cancer Genome Atlas. LASSO regression was further applied to select optimal prognosis-related genes. The final prognostic nomogram was implemented together with prognostic clinicopathological parameters. The predictive efficiency of the nomogram was internally validated in training and testing groups, and externally validated in the Chinese Glioma Genome Atlas database. Three DNA methylation-driven genes, ARL9, CMYA5, and STEAP3, were identified as independent prognostic factors. Together with IDH1 mutation status, age, and sex, the final prognostic nomogram achieved the highest AUC value of 0.930, and demonstrated stable consistency in both internal and external validations. The prognostic nomogram could predict personal survival probabilities for patients with LGG, and serve as a user-friendly tool for prognostic evaluation, optimizing therapeutic regimes, and managing LGG patients.

Screening and Prognostic Analysis of Immune-Related Genes in Pancreatic Cancer

Pancreatic cancer remains to have a high mortality, which is partly due to the lack of effective treatment strategies. In this study, genes with potential associations with immunophenotyping of pancreatic cancer were screened through bioinformatics analysis and the correlation between immune-related genes and the prognosis of pancreatic cancer patients was assessed. Firstly, differentially expressed immune genes were extracted from the pancreatic cancer-related datasets obtained for purposes of this study. The samples were processed by the "Consensus Cluster Plus" R package to determine the number of immune subtypes. Then, the pancreatic cancer immunophenotyping-related gene modules were determined. Differential analysis of immune gene modules was performed, and the function of genes related to pancreatic cancer immune subtypes was identified. The number of immune cells in the samples was calculated, followed by the differential expression analysis of immune cell numbers in each immune subtype of pancreatic cancer. The immune infiltration score was also estimated, and the correlation between the immune infiltration score and the patient prognosis with different immune subtypes was determined. Gene differences between each immune subtype were identified by differential expression analysis, and key immune genes affecting immunophenotyping were obtained. Following the analysis, 426 immune-related genes were identified to have potential involvement in the occurrence and development of pancreatic cancer, of which CD19 may be the most critical gene affecting the immunophenotyping of pancreatic cancer. CD19 played a significant role in the occurrence and development of IS2 and IS3 immune subtypes of pancreatic cancer through its action on B cells and T cells. Moreover, the expression of CD19 was increased in the collected pancreatic cancer tissues. Overall, our findings uncovered the critical role of CD19 in the prognosis of pancreatic cancer patients.

Establishment of an Immune-Related Gene Signature for Risk Stratification for Patients with Glioma

Glioma is a frequently seen primary malignant intracranial tumor, characterized by poor prognosis. The study is aimed at constructing a prognostic model for risk stratification in patients suffering from glioma. Weighted gene coexpression network analysis (WGCNA), integrated transcriptome analysis, and combining immune-related genes (IRGs) were used to identify core differentially expressed IRGs (DE IRGs). Subsequently, univariate and multivariate Cox regression analyses were utilized to establish an immune-related risk score (IRRS) model for risk stratification for glioma patients. Furthermore, a nomogram was developed for predicting glioma patients' overall survival (OS). The turquoise module (cor = 0.67; P < 0.001) and its genes (n = 1092) were significantly pertinent to glioma progression. Ultimately, multivariate Cox regression analysis constructed an IRRS model based on VEGFA, SOCS3, SPP1, and TGFB2 core DE IRGs, with a C-index of 0.811 (95% CI: 0.786-0.836). Then, Kaplan-Meier (KM) survival curves revealed that patients presenting high risk had a dismal outcome (P < 0.0001). Also, this IRRS model was found to be an independent prognostic indicator of gliomas' survival prediction, with HR of 1.89 (95% CI: 1.252-2.85) and 2.17 (95% CI: 1.493-3.14) in the Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) datasets, respectively. We established the IRRS prognostic model, capable of effectively stratifying glioma population, convenient for decision-making in clinical practice.

Construction of a Promising Tumor-Infiltrating CD8+ T Cells Gene Signature to Improve Prediction of the Prognosis and Immune Response of Uveal Melanoma

Background

CD8+ T cells work as a key effector of adaptive immunity and are closely associated with immune response for killing tumor cells. It is crucial to understand the role of tumor-infiltrating CD8+ T cells in uveal melanoma (UM) to predict the prognosis and response to immunotherapy.

Materials and Methods

Single-cell transcriptomes of UM with immune-related genes were combined to screen the CD8+ T-cell-associated immune-related genes (CDIRGs) for subsequent analysis. Next, a prognostic gene signature referred to tumor-infiltrating CD8+ T cells was constructed and validated in several UM bulk RNA sequencing datasets. The risk score of UM patients was calculated and classified into high- or low-risk subgroup. The prognostic value of risk score was estimated by using multivariate Cox analysis and Kaplan–Meier survival analysis. Moreover, the potential ability of gene signature for predicting immunotherapy response was further explored.

Results

In total, 202 CDIRGs were screened out from the single-cell RNA sequencing of GSE139829. Next, a gene signature containing three CDIRGs (IFNGR1, ANXA6, and TANK) was identified, which was considered as an independent prognostic indicator to robustly predict overall survival (OS) and metastasis-free survival (MFS) of UM. In addition, the UM patients were classified into high- and low-risk subgroups with different clinical characteristics, distinct CD8+ T-cell immune infiltration, and immunotherapy response. Gene set enrichment analysis (GSEA) showed that immune pathways such as allograft rejection, inflammatory response, interferon alpha and gamma response, and antigen processing and presentation were all positively activated in low-risk phenotype.

Conclusion

Our work gives an inspiration to explain the limited response for the current immune checkpoint inhibitors to UM. Besides, we constructed a novel gene signature to predict prognosis and immunotherapy responses, which may be regarded as a promising therapeutic target.

An Effective Hypoxia-Related Long Non-Coding RNAs Assessment Model for Prognosis of Clear Cell Renal Carcinoma

Hypoxia is a significant clinical feature and regulates various tumor processes in clear cell renal carcinoma (ccRCC). Increasing evidence has demonstrated that long non-coding RNAs (lncRNAs) are closely associated with the survival outcomes of ccRCC patients and regulates hypoxia-induced tumor processes. Thus, this study aimed to develop a hypoxia-related lncRNA (HRL) prognostic model for predicting the survival outcomes in ccRCC. LncRNAs in ccRCC samples were extracted from The Cancer Genome Atlas database. Hypoxia-related genes were downloaded from the Molecular Signatures Database. A co-expression analysis between differentially expressed lncRNAs and hypoxia-related genes in ccRCC samples was performed to identify HRLs. Univariate and multivariate Cox regression analyses were performed to select nine optimal lncRNAs for developing the HRL model. The prognostic model showed good performance in predicting prognosis among patients with ccRCC, and the validation sets reached consistent results. The model was also found to be related to the clinicopathologic parameters of tumor grade and tumor stage and to tumor immune infiltration. In conclusion, our findings indicate that the hypoxia-lncRNA assessment model may be useful for prognostication in ccRCC cases. Furthermore, the nine HRLs included in the model might be useful targets for investigating the tumorigenesis of ccRCC and designing individualized treatment strategies.

Identification of Hub Genes Associated With Immune Infiltration and Predict Prognosis in Hepatocellular Carcinoma via Bioinformatics Approaches

Aims

In the cancer-related research field, there is currently a major need for a greater number of valuable biomarkers to predict the prognosis of hepatocellular carcinoma (HCC). In this study, we aimed to screen hub genes related to immune cell infiltration and explore their prognostic value for HCC.

Methods

We analyzed five datasets (GSE46408, GSE57957, GSE74656, GSE76427, and GSE87630) from the Gene Expression Omnibus database to screen the differentially expressed genes (DEGs). A protein-protein interaction network of the DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes; then, the hub genes were identified. Functional enrichment of the genes was performed on the Metascape website. Next, the expression of these hub genes was validated in several databases, including Oncomine, Gene Expression Profiling Interactive Analysis 2 (GEPIA2), and Human Protein Atlas. We explored the correlations between the hub genes and infiltrated immune cells in the TIMER2.0 database. The survival curves were generated in GEPIA2, and the univariate and multivariate Cox regression analyses were performed using TIMER2.0.

Results

The top ten hub genes [DNA topoisomerase II alpha (TOP2A), cyclin B2 (CCNB2), protein regulator of cytokinesis 1 (PRC1), Rac GTPase-activating protein 1 (RACGAP1), aurora kinase A (AURKA), cyclin-dependent kinase inhibitor 3 (CDKN3), nucleolar and spindle-associated protein 1 (NUSAP1), cell division cycle-associated 5 (CDCA5), abnormal spindle microtubule assembly (ASPM), and non-SMC condensin I complex subunit G (NCAPG)] were identified in subsequent analysis. These genes are most markedly enriched in cell division, suggesting their close association with tumorigenesis. Multi-database analyses validated that the hub genes were upregulated in HCC tissues. All hub genes positively correlated with several types of immune infiltration, including B cells, CD4⁺ T cells, macrophages, and dendritic cells. Furthermore, these hub genes served as independent prognostic factors, and the expression of these hub genes combing with the macrophage levels could help predict an unfavorable prognosis of HCC.

Conclusion

In sum, these hub genes (TOP2A, CCNB2, PRC1, RACGAP1, AURKA, CDKN3, NUSAP1, CDCA5, ASPM, and NCAPG) may be pivotal markers for prognostic prediction as well as potentially work as targets for immune-based intervention strategies in HCC.

Identification of a Five-Autophagy-Related-lncRNA Signature as a Novel Prognostic Biomarker for Hepatocellular Carcinoma

Although great progresses have been made in the diagnosis and treatment of hepatocellular carcinoma (HCC), its prognostic marker remains controversial. In this current study, weighted correlation network analysis and Cox regression analysis showed significant prognostic value of five autophagy-related long non-coding RNAs (AR-lncRNAs) (including TMCC1-AS1, PLBD1-AS1, MKLN1-AS, LINC01063, and CYTOR) for HCC patients from data in The Cancer Genome Atlas. By using them, we constructed a five-AR-lncRNA prognostic signature, which accurately distinguished the high- and low-risk groups of HCC patients. All of the five AR lncRNAs were highly expressed in the high-risk group of HCC patients. This five-AR-lncRNA prognostic signature showed good area under the curve (AUC) value (AUC = 0.751) for the overall survival (OS) prediction in either all HCC patients or HCC patients stratified according to several clinical traits. A prognostic nomogram with this five-AR-lncRNA signature predicted the 3- and 5-year OS outcomes of HCC patients intuitively and accurately (concordance index = 0.745). By parallel comparison, this five-AR-lncRNA signature has better prognosis accuracy than the other three recently published signatures. Furthermore, we discovered the prediction ability of the signature on therapeutic outcomes of HCC patients, including chemotherapy and immunotherapeutic responses. Gene set enrichment analysis and gene mutation analysis revealed that dysregulated cell cycle pathway, purine metabolism, and TP53 mutation may play an important role in determining the OS outcomes of HCC patients in the high-risk group. Collectively, our study suggests a new five-AR-lncRNA prognostic signature for HCC patients.