Circadian rhythm-related genes index: A predictor for HNSCC prognosis, immunotherapy efficacy, and chemosensitivity

Deciphering COPS5 influence on immune infiltration and prognosis in head and neck squamous cell carcinoma

Head and Neck Squamous Cell Carcinoma (HNSCC) is a widespread malignancy originating from the mucous epithelium of the oral cavity, pharynx, and larynx. Despite advances in diagnostic and therapeutic modalities, the prognosis of HNSCC remains challenging. This study investigates the intricate relationship among COPS5, immune infiltration patterns, and prognostic implications in HNSCC. Through comprehensive analyses of 519 HNSCC cases from TCGA and single-cell data from the GEO database, we utilize the CIBERSORT algorithm to discern immune cell dynamics influenced by COPS5 expression. Notably, Treg cells emerge as a central point in the interplay between COPS5 and immune modulation. Further analyses, encompassing differential gene expression, immune-related gene set enrichment, and protein-protein interaction networks, elucidate the molecular landscape associated with COPS5 in HNSCC. A prognostic risk model, incorporating CD27, TNFRSF4, FADD, and PSMD14, is formulated and validated across diverse datasets. The model demonstrates robust predictive power, underscoring its potential as a valuable prognostic tool. These genes, essential for immune regulation and cell cycle control, provide insights into the intricate mechanisms influencing HNSCC progression. In conclusion, this study not only reveals the impact of COPS5 on immune dynamics in HNSCC but also introduces a concise and effective prognostic model.

Construction of a circadian rhythm-relevant gene signature for hepatocellular carcinoma prognosis, immunotherapy and chemosensitivity prediction

Aims

This study explored the molecular and biologic mechanisms underlying the association between circadian rhythm disorders (CRD) and increased risk for hepatocellular carcinoma (HCC).

Background

CRD are linked to increased risk for HCC, but the molecular and biologic mechanisms underlying this association are limited.ObjectiveThe study constructed and validated a CRD related gene model as an independent prognostic factor for HCC, providing insight into the molecular mechanisms linking CRD to increased HCC risk and identifying potential indicators for the efficacy of immunotherapy and anticancer drugs. This helps provide important clues for personalized treatment strategies for HCC patients.

Methods

Gene sets correlated with circadian rhythm were obtained from the Molecular Signatures Database (MSigDB) to intersect with differentially expressed genes (DEGs) between tumor samples and control samples in The Cancer Genome Atlas (TCGA) and HCCDB18 from Hepatocellular Carcinoma Cell DataBase (HCCDB). The CRD related gene model was developed by univariate Cox and stepwise multivariate analysis. Immune checkpoint blockade (ICB) therapy and anticancer drugs were analyzed using the tumor immune dysfunction and exclusion (TIDE) and pRRophetic, respectively. Seurat determined the cell type of HCC by analyzing single-cell data, and malignant cells were identified using Copykat. To detect the mRNA levels of genes in the CRD related gene model, quantitative real-time polymerase chain reaction (qRT-PCR) was carried out.

Results

The activity of circadian rhythm in HCC tissue was significantly lower than that in control tissue. Subsequently, EZH2, IMPDH2, TYMS and SERPINE1 were selected to construct the CRD related gene model, which was an independent factor for HCC prognosis. Notably, low-risk patients had lower levels of immune cell infiltration and lower TIDE scores compared to high-risk patients with HCC, indicating that patients with a low risk may derive more benefit from immunotherapy. IMPDH2, TYMS and SERPINE1 expressed significantly higher in malignant cells than in benign epithelial cells.

Conclusions

This study presents a CRD related gene model to reveal the molecular perspective of the dependent mechanism of the association between CRD and cancer, which provides a potential indicator for understanding the preclinical efficacy of ICB and anticancer drugs.

Characterizing mitochondrial features in osteoarthritis through integrative multi-omics and machine learning analysis

Purpose

Osteoarthritis (OA) stands as the most prevalent joint disorder. Mitochondrial dysfunction has been linked to the pathogenesis of OA. The main goal of this study is to uncover the pivotal role of mitochondria in the mechanisms driving OA development.

Materials and methods

We acquired seven bulk RNA-seq datasets from the Gene Expression Omnibus (GEO) database and examined the expression levels of differentially expressed genes related to mitochondria in OA. We utilized single-sample gene set enrichment analysis (ssGSEA), gene set enrichment analysis (GSEA), and weighted gene co-expression network analysis (WGCNA) analyses to explore the functional mechanisms associated with these genes. Seven machine learning algorithms were utilized to identify hub mitochondria-related genes and develop a predictive model. Further analyses included pathway enrichment, immune infiltration, gene-disease relationships, and mRNA-miRNA network construction based on these hub mitochondria-related genes. genome-wide association studies (GWAS) analysis was performed using the Gene Atlas database. GSEA, gene set variation analysis (GSVA), protein pathway analysis, and WGCNA were employed to investigate relevant pathways in subtypes. The Harmonizome database was employed to analyze the expression of hub mitochondria-related genes across various human tissues. Single-cell data analysis was conducted to examine patterns of gene expression distribution and pseudo-temporal changes. Additionally, The real-time polymerase chain reaction (RT-PCR) was used to validate the expression of these hub mitochondria-related genes.

Results

In OA, the mitochondria-related pathway was significantly activated. Nine hub mitochondria-related genes (SIRT4, DNAJC15, NFS1, FKBP8, SLC25A37, CARS2, MTHFD2, ETFDH, and PDK4) were identified. They constructed predictive models with good ability to predict OA. These genes are primarily associated with macrophages. Unsupervised consensus clustering identified two mitochondria-associated isoforms that are primarily associated with metabolism. Single-cell analysis showed that they were all expressed in single cells and varied with cell differentiation. RT-PCR showed that they were all significantly expressed in OA.

Conclusion

SIRT4, DNAJC15, NFS1, FKBP8, SLC25A37, CARS2, MTHFD2, ETFDH, and PDK4 are potential mitochondrial target genes for studying OA. The classification of mitochondria-associated isoforms could help to personalize treatment for OA patients.

The expression of immune co-stimulators as a prognostic predictor of head and neck squamous cell carcinomas and oral squamous cell carcinomas

Background/purpose

T cells require second immune checkpoint molecules for activation and immune memory after antigen presentation. We found that inducible co-stimulator (ICOS) has been a favorable prognostic factor amongst B7 immune checkpoint co-stimulators (ICSs) families in head and neck squamous cell carcinoma (HNSCC) and oral SCC (OSCC).

Materials and methods

This study analyzed the expression of non-B7 tumor necrosis factor (TNF) superfamily ICSs in the Cancer Genome Atlas (TCGA) HNSCC cohort, our OSCC cohort, and TCGA pan-cancer datasets. The correlation in expression, prognosis, and immune status was assessed.

Results

The higher expression of CD27, CD30, CD40L, death domain 3 (DR3), and OX40, presumably on the T cell surface, defined better overall survival of HNSCC patients. Besides, CD27, CD30, CD40L, and OX40 were highly correlated with ICOS expression in tumors. CD27, CD40L, and DR3 expression are higher in HPV+ HNSCC tumors than in HPV- tumors. The combined expression level of CD27/OX40 or CD27/CD40L/OX40 enables the potent survival prediction of small, less nodal involvement, early stage, and HPV + tumor subsets. Tumors expressing high CD27, CD30, CD40L, ICOS, and OX40 exhibited enhanced immune cell infiltration. The high correlation in the expression of these ICSs was also noted in the vast majority of tumor types in TCGA datasets.

Conclusion

The findings of this study not only confirm the potential of the concordant stimulation of CD27, CD30, CD40L, ICOS, and OX40 as a crucial strategy in cancer immunotherapy but also inspire further exploration into the field, highlighting the promising future of cancer treatment.

Multi-omic validation of the cuproptosis-sphingolipid metabolism network: modulating the immune landscape in osteosarcoma

Background

The current understanding of the mechanisms by which metal ion metabolism promotes the progression and drug resistance of osteosarcoma remains incomplete. This study aims to elucidate the key roles and mechanisms of genes involved in cuproptosis-related sphingolipid metabolism (cuproptosis-SPGs) in regulating the immune landscape, tumor metastasis, and drug resistance in osteosarcoma cells.

Methods

This study employed multi-omics approaches to assess the impact of cuproptosis-SPGs on the prognosis of osteosarcoma patients. Lasso regression analysis was utilized to construct a prognostic model, while multivariate regression analysis was applied to identify key core genes and generate risk coefficients for these genes, thereby calculating a risk score for each osteosarcoma patient. Patients were then stratified into high-risk and low-risk groups based on their risk scores. The ESTIMATE and CIBERSORT algorithms were used to analyze the level of immune cell infiltration within these risk groups to construct the immune landscape. Single-cell analysis was conducted to provide a more precise depiction of the expression patterns of cuproptosis-SPGs among immune cell subtypes. Finally, experiments on osteosarcoma cells were performed to validate the role of the cuproptosis-sphingolipid signaling network in regulating cell migration and apoptosis.

Results

In this study, seven cuproptosis-SPGs were identified and used to construct a prognostic model for osteosarcoma patients. In addition to predicting survival, the model also demonstrated reliability in forecasting the response to chemotherapy drugs. The results showed that a high cuproptosis-sphingolipid metabolism score was closely associated with reduced CD8 T cell infiltration and indicated poor prognosis in osteosarcoma patients. Cellular functional assays revealed that cuproptosis-SPGs regulated the LC3B/ERK signaling pathway, thereby triggering cell death and impairing migration capabilities in osteosarcoma cells.

Conclusion

The impact of cuproptosis-related sphingolipid metabolism on the survival and migration of osteosarcoma cells, as well as on CD8 T cell infiltration, highlights the potential of targeting copper ion metabolism as a promising strategy for osteosarcoma patients.

Effect of time-of-day nivolumab and stereotactic body radiotherapy in metastatic head and neck squamous cell carcinoma: A secondary analysis of a prospective randomized trial

BACKGROUND

Prior work documented circadian rhythm impacts on efficacy and toxicity of cancer therapies.

METHODS

Secondary analysis of prospective, phase II trial of metastatic HNSCC randomized to nivolumab+/-SBRT. Used cutoffs of 1100 and 1630. Timing classified by first infusion or majority of SBRT (e.g., PM SBRT defined by two or three fractions after 1630).

RESULTS

Of 62 patients, there was no significant difference in median PFS between AM nivolumab (n = 7, 175 days), PM nivolumab (n = 21, 58 days), or Mid-Day nivolumab (n = 34, 67 days; p = 0.8). There was no significant difference in median PFS with AM SBRT (n = 4, 78 days), PM SBRT (n = 13, 111 days), or Mid-Day SBRT (n = 15, 63 days; p = 0.8). There was no significant difference in Grade 3-4 toxicity or ORR. Sensitivity analyses with other timepoints were negative.

CONCLUSIONS

Further work may elucidate circadian impacts on select patients, tumors, and therapies; however, we found no significant effect in this study.

KHDRBS1 as a novel prognostic signaling biomarker influencing hepatocellular carcinoma cell proliferation, migration, immune microenvironment, and drug sensitivity

Background

Human tumors pose significant challenges, with targeted therapy against specific molecular targets or signaling pathways being a mainstay alongside surgical resection. Previous studies have implicated KHDRBS1 in the oncogenesis of certain human tumors such as colorectal and prostate cancers, underscoring its potential as a therapeutic target. However, the comprehensive expression pattern of KHDRBS1 in hepatocellular carcinoma (HCC) warrants further exploration.

Methods

Integrating and analyzing multi-omics, multi-cohort data from public databases, coupled with clinical samples and molecular biology validation, we elucidate the oncogenic role of KHDRBS1 in HCC progression. Additionally, leveraging HCC single-cell sequencing data, we segregate malignant cells into KHDRBS1-positive and negative subsets, uncovering significant differences in their expression profiles and functional roles.

Results

Our study identifies KHDRBS1 as a tumor-promoting factor in HCC, with its positivity correlating with tumor progression. Furthermore, we highlight the clinical significance of KHDRBS1-positive malignant cells, aiming to further propel its clinical utility.

Conclusion

KHDRBS1 plays a key role in HCC development. This study provides crucial insights for further investigation into KHDRBS1 as a therapeutic target in HCC.

Experimentally validated oxidative stress -associated prognostic signatures describe the immune landscape and predict the drug response and prognosis of SKCM

Background

Skin Cutaneous Melanoma (SKCM) incidence is continually increasing, with chemotherapy and immunotherapy being among the most common cancer treatment modalities. This study aims to identify novel biomarkers for chemotherapy and immunotherapy response in SKCM and explore their association with oxidative stress.

Methods

Utilizing TCGA-SKCM RNA-seq data, we employed Weighted Gene Co-expression Network Analysis (WGCNA) and Protein-Protein Interaction (PPI) networks to identify six core genes. Gene co-expression analysis and immune-related analysis were conducted, and specific markers associated with oxidative stress were identified using Gene Set Variation Analysis (GSVA). Single-cell analysis revealed the expression patterns of Oxidative Stress-Associated Genes (OSAG) in the tumor microenvironment. TIDE analysis was employed to explore the association between immune therapy response and OSAG, while CIBERSORT was used to analyze the tumor immune microenvironment. The BEST database demonstrated the impact of the Oxidative Stress signaling pathway on chemotherapy drug resistance. Immunohistochemical staining and ROC curve evaluation were performed to assess the protein expression levels of core genes in SKCM and normal samples, with survival analysis utilized to determine their diagnostic value.

Results

We identified six central genes associated with SKCM metastasis, among which the expression of DSC2 and DSC3 involved in the oxidative stress pathway was closely related to immune cell infiltration. DSC2 influenced drug resistance in SKMC patients. Furthermore, downregulation of DSC2 and DSC3 expression enhanced the response of SKCM patients to immunotherapy.

Conclusion

This study identified two Oxidative Stress-Associated genes as novel biomarkers for SKCM. Additionally, targeting the oxidative stress pathway may serve as a new strategy in clinical practice to enhance SKCM chemotherapy and sensitivity.

Regulatory T cells subgroups in the tumor microenvironment cannot be overlooked: Their involvement in prognosis and treatment strategy in melanoma

BACKGROUND

Melanoma, the most lethal form of skin cancer, presents substantial challenges despite effective surgical interventions for in situ lesions. Regulatory T cells (Tregs) wield a pivotal immunomodulatory influence within the tumor microenvironment, yet their impact on melanoma prognosis and direct molecular interactions with melanoma cells remain elusive. This investigation employs single-cell analysis to unveil the intricate nature of Tregs in human melanoma.

METHODS

Single-cell RNA and bulk sequencing data, alongside clinical information, were obtained from public repositories. Initially, GO and GSEA analyses were employed to delineate functional disparities among distinct cell subsets. Pseudotime and cell-cell interconnection analyses were conducted, followed by an endeavor to construct a prognostic model grounded in Treg-associated risk scores. This model's efficacy was demonstrated via PCA and K-M analyses, with multivariate Cox regression affirming its independent prognostic value in melanoma patients. Furthermore, immune infiltration analysis, immune checkpoint gene expression scrutiny, and drug sensitivity assessments were performed to ascertain the clinical relevance of this prognostic model.

RESULTS

Following batch effect correction, 80 025 cells partitioned into 31 clusters, encompassing B cells, plasma cells, endothelial cells, fibroblasts, melanoma cells, monocytes, macrophages, and T_NK cells. Within these, 4240 CD4+ T cells were subclassified into seven distinct types. Functional analysis underscored the immunomodulatory function of Tregs within the melanoma tumor microenvironment, elucidating disparities among Treg subpopulations. Notably, the ITGB2 signaling pathway emerged as a plausible molecular nexus linking Tregs to melanoma cells. Our prognostic signature exhibited robust predictive capacities for melanoma prognosis and potential implications in evaluating immunotherapy response.

CONCLUSION

Tregs exert a critical role in immune suppression within the melanoma tumor microenvironment, revealing a potential molecular-level association with melanoma cells. Our innovative Treg-centered signature introduces a promising prognostic marker for melanoma, holding potential for future clinical prognostic assessments.

Unraveling the role of the circadian clock genes in cervical squamous cell carcinoma and endocervical adenocarcinoma: A prognostic indicator for prognostic, immunotherapy response, and chemotherapy sensitivity

Background: Cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) account for a significant proportion of gynecological malignancies and represent a major global health concern. Globally, CESC is ranked as the fourth most common cancer among women. Conventional treatment of this disease has a less favorable prognosis for most patients. However, the discovery of early molecular biomarkers is therefore important for the diagnosis of CESC, as well as for slowing down their progression process. Methods: To identify differentially expressed genes strongly associated with prognosis, univariate Cox proportional hazard analysis and least absolute shrinkage and selection operator (LASSO) regression analysis were used. Using multiple Cox proportional hazard regression, a multifactorial model for prognostic risk assessment was then created. Results: The expression of biological clock-related genes, which varied considerably among distinct subtypes and were associated with significantly diverse prognoses, was used to categorize CESC patients. These findings demonstrate how the nomogram developed based on the 7-CRGs signature may assist physicians in creating more precise, accurate, and successful treatment plans that can aid CESC patients at 1, 3, and 5 years. Conclusions: By using machine learning techniques, we thoroughly investigated the impact of CRGs on the prognosis of CESC patients in this study. By creating a unique nomogram, we were able to accurately predict patient prognosis. At the same time, we showed new perspectives on the development of CESC and its treatment by analyzing the associations of the prognostic model with immunity, enrichment pathways, chemotherapy sensitivity, and so on. This research provides a new direction for clinical treatment.

The comprehensive landscape of prognosis, immunity, and function of the GLI family by pan-cancer and single-cell analysis

The Hedgehog (Hh) signaling pathway has been implicated in the pathogenesis of various cancers. However, the roles of the downstream GLI family (GLI1, GLI2, and GLI3) in tumorigenesis remain elusive. This study aimed to unravel the genetic alterations of GLI1/2/3 in cancer and their association with the immune microenvironment and related signaling pathways. Firstly, we evaluated the expression profiles of GLI1/2/3 in different cancer types, analyzed their prognostic and predictive values, and assessed their correlation with tumor-infiltrating immune cells. Secondly, we explored the relationships between GLI1/2/3 and genetic mutations, epigenetic modifications, and clinically relevant drugs. Finally, we performed enrichment analysis to decipher the underlying mechanisms of GLI1/2/3 in cancer initiation and progression. Our results revealed that the expression levels of GLI1/2/3 were positively correlated in most cancer tissues, suggesting a cooperative role of these factors in tumorigenesis. We also identified tissue-specific expression patterns of GLI1/2/3, which may reflect the distinct functions of these factors in different cell types. Furthermore, GLI1/2/3 expression displayed significant associations with poor prognosis in several cancers, indicating their potential as prognostic biomarkers and therapeutic targets. Importantly, we found that GLI1/2/3 modulated the immune microenvironment by regulating the recruitment, activation, and polarization of cancer-associated fibroblasts, endothelial cells, and macrophages. Additionally, functional enrichment analyses indicated that GLI1/2/3 are involved in the regulation of epithelial-mesenchymal transition (EMT). Together, our findings shed new light on the roles of GLI1/2/3 in tumorigenesis and provide a potential basis for the development of novel therapeutic strategies targeting GLI-mediated signaling pathways in cancer.

Innovative prognostic modeling in ESCC: leveraging scRNA-seq and bulk-RNA for dendritic cell heterogeneity analysis

Background

Globally, esophageal squamous cell carcinoma (ESCC) stands out as a common cancer type, characterized by its notably high rates of occurrence and mortality. Recent advancements in treatment methods, including immunotherapy, have shown promise, yet the prognosis remains poor. In the context of tumor development and treatment outcomes, the tumor microenvironment (TME), especially the function of dendritic cells (DCs), is significantly influential. Our study aims to delve deeper into the heterogeneity of DCs in ESCC using single-cell RNA sequencing (scRNA-seq) and bulk RNA analysis.

Methods

In the scRNA-seq analysis, we utilized the SCP package for result visualization and functional enrichment analysis of cell subpopulations. CellChat was employed to identify potential oncogenic mechanisms in DCs, while Monocle 2 traced the evolutionary trajectory of the three DC subtypes. CopyKAT assessed the benign or malignant nature of cells, and SCENIC conducted transcription factor regulatory network analysis, offering a preliminary exploration of DC heterogeneity. In Bulk-RNA analysis, we constructed a prognostic model for ESCC prognosis and immunotherapy response, based on DC marker genes. This model was validated through quantitative PCR (qPCR) and immunohistochemistry (IHC), confirming the gene expression levels.

Results

In this study, through intercellular communication analysis, we identified GALECTIN and MHC-I signaling pathways as potential oncogenic mechanisms within dendritic cells. We categorized DCs into three subtypes: plasmacytoid (pDC), conventional (cDC), and tolerogenic (tDC). Our findings revealed that pDCs exhibited an increased proportion of cells in the G2/M and S phases, indicating enhanced cellular activity. Pseudotime trajectory analysis demonstrated that cDCs were in early stages of differentiation, whereas tDCs were in more advanced stages, with pDCs distributed across both early and late differentiation phases. Prognostic analysis highlighted a significant correlation between pDCs and tDCs with the prognosis of ESCC (P< 0.05), while no significant correlation was observed between cDCs and ESCC prognosis (P = 0.31). The analysis of cell malignancy showed the lowest proportion of malignant cells in cDCs (17%), followed by pDCs (29%), and the highest in tDCs (48%), with these results being statistically significant (P< 0.05). We developed a robust ESCC prognostic model based on marker genes of pDCs and tDCs in the GSE53624 cohort (n = 119), which was validated in the TCGA-ESCC cohort (n = 139) and the IMvigor210 immunotherapy cohort (n = 298) (P< 0.05). Additionally, we supplemented the study with a novel nomogram that integrates clinical features and risk assessments. Finally, the expression levels of genes involved in the model were validated using qPCR (n = 8) and IHC (n = 16), thereby confirming the accuracy of our analysis.

Conclusion

This study enhances the understanding of dendritic cell heterogeneity in ESCC and its impact on patient prognosis. The insights gained from scRNA-seq and Bulk-RNA analysis contribute to the development of novel biomarkers and therapeutic targets. Our prognostic models based on DC-related gene signatures hold promise for improving ESCC patient stratification and guiding treatment decisions.

Glycogen metabolism-mediated intercellular communication in the tumor microenvironment influences liver cancer prognosis

Glycogen metabolism plays a key role in the development of hepatocellular carcinoma (HCC), but the function of glycogen metabolism genes in the tumor microenvironment (TME) is still to be elucidated. Single-cell RNA-seq data were obtained from ten HCC tumor samples totaling 64,545 cells, and 65 glycogen metabolism genes were analyzed by a nonnegative matrix factorization (NMF). The prognosis and immune response of new glycogen TME cell clusters were predicted by using HCC and immunotherapy cohorts from public databases. HCC single-cell analysis was divided into fibroblasts, NT T cells, macrophages, endothelial cells, and B cells, which were separately divided into new cell clusters by glycogen metabolism gene annotation. Pseudo-temporal trajectory analysis demonstrated the temporal differentiation trajectory of different glycogen subtype cell clusters. Cellular communication analysis revealed extensive interactions between endothelial cells with glycogen metabolizing TME cell-related subtypes and different glycogen subtype cell clusters. SCENIC analysis of transcription factors upstream of TME cell clusters with different glycogen metabolism. In addition, TME cell clusters of glycogen metabolism were found to be enriched in expression in CAF subtypes, CD8 depleted, M1, and M2 types. Bulk-seq analysis showed the prognostic significance of glycogen metabolism-mediated TME cell clusters in HCC, while a significant immune response was found in the immunotherapy cohort in patients treated with immune checkpoint blockade (ICB), especially for CAFs, T cells, and macrophages. In summary, our study reveals for the first time that glycogen metabolism mediates intercellular communication in the hepatocellular carcinoma microenvironment while elucidating the anti-tumor mechanisms and immune prognostic responses of different subtypes of cell clusters.

Diagnostic plasma small extracellular vesicles miRNA signatures for pancreatic cancer using machine learning methods

BACKGROUND

Identifying biomarkers may lead to easier detection and a better understanding of pathogenesis of pancreatic ductal adenocarcinoma (PDAC).

METHODS

Plasma small extracellular vesicles (sEV) from 106 participants, including 20 healthy controls (HC), 12 chronic pancreatitis (CP) patients, 12 benign pancreatic tumour (BPT) patients, and 58 PDAC patients, were profiled for microRNA (miRNA) sequencing. Three machine learning methods were applied to establish and evaluate the diagnostic model.

RESULTS

The plasma sEV miRNA diagnostic signature (d-signature) selected using the three machine learning methods could distinguish PDAC patients from non-PDAC individuals, HC, and benign pancreatic disease (BPD, CP plus BPT) both in training and validation cohort. Combining the d-signature with carbohydrate antigen 19-9 (CA19-9) performed better than with each model alone. Plasma sEV miR-664a-3p was selected by all methods and used to predict PDAC diagnosis with high accuracy combined with CA19-9. Plasma sEV miR-664a-3p was significantly positively associated with the presence of vascular invasion, lower surgery ratio, and poor differentiation. MiR-664a-3p was mainly distributed in the PDAC cancer stroma, including fibers and vessels, and was accompanied by VEGFA expression. Overexpression of miR-664a-3p could promote the epithelial-mesenchymal transition (EMT) and angiogenesis.

CONCLUSION

In conclusion, our study demonstrated the potential utility of the sEV-miRNA d-signature in the diagnosis of PDAC via machine learning methods. A novel sEV biomarker, miR-664a-3p, was identified for the diagnosis of PDAC. It can also potentially promote angiogenesis and metastasis, provide insight into PDAC pathogenesis, and reveal novel regulators of this disease.

Comprehending the cuproptosis and cancer-immunity cycle network: delving into the immune landscape and its predictive role in breast cancer immunotherapy responses and clinical endpoints

Background

The role of cuproptosis, a phenomenon associated with tumor metabolism and immunological identification, remains underexplored, particularly in relation to the cancer-immunity cycle (CIC) network. This study aims to rigorously examine the impact of the cuproptosis-CIC nexus on immune reactions and prognostic outcomes in patients with breast cancer (BC), striving to establish a comprehensive prognostic model.

Methods

In the study, we segregated data obtained from TCGA, GEO, and ICGC using CICs retrieved from the TIP database. We constructed a genetic prognostic framework using the LASSO-Cox model, followed by its validation through Cox proportional hazards regression. This framework's validity was further confirmed with data from ICGC and GEO. Explorations of the tumor microenvironment were carried out through the application of ESTIMATE and CIBERSORT algorithms, as well as machine learning techniques, to identify potential treatment strategies. Single-cell sequencing methods were utilized to delineate the spatial distribution of key genes within the various cell types in the tumor milieu. To explore the critical role of the identified CICs, experiments were conducted focusing on cell survival and migration abilities.

Results

In our research, we identified a set of 4 crucial cuproptosis-CICs that have a profound impact on patient longevity and their response to immunotherapy. By leveraging these identified CICs, we constructed a predictive model that efficiently estimates patient prognoses. Detailed analyses at the single-cell level showed that the significance of CICs. Experimental approaches, including CCK-8, Transwell, and wound healing assays, revealed that the protein HSPA9 restricts the growth and movement of breast cancer cells. Furthermore, our studies using immunofluorescence techniques demonstrated that suppressing HSPA9 leads to a notable increase in ceramide levels.

Conclusion

This research outlines a network of cuproptosis-CICs and constructs a predictive nomogram. Our model holds great promise for healthcare professionals to personalize treatment approaches for individuals with breast cancer. The work provides insights into the complex relationship between the cuproptosis-CIC network and the cancer immune microenvironment, setting the stage for novel approaches to cancer immunotherapy. By focusing on the essential gene HSPA9 within the cancer-immunity cycle, this strategy has the potential to significantly improve the efficacy of treatments against breast cancer.

Systematic analysis of the prognostic value and immunological function of LTBR in human cancer

Lymphotoxin beta receptor (LTBR) is a positive T cell proliferation regulator gene. It is closely associated with the tumor immune microenvironment. However, its role in cancer and immunotherapy is unclear. Firstly, the expression level and prognostic value of LTBR were analyzed. Secondly, the expression of LTBR in clinical stages, immune subtypes, and molecular subtypes was analyzed. The correlation between LTBR and immune regulatory genes, immune checkpoint genes, and RNA modification genes was then analyzed. Correlations between LTBR and immune cells, scores, cancer-related functional status, tumor stemness index, mismatch repair (MMR) genes, and DNA methyltransferase were also analyzed. In addition, we analyzed the role of LTBR in DNA methylation, mutational status, tumor mutation burden (TMB), and microsatellite instability (MSI). Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA) were used to explore the role of LTBR in pan-cancer. Finally, the drugs associated with LTBR were analyzed. The expression of LTBR was confirmed using quantitative real-time PCR and Western blot. LTBR is significantly overexpressed in most cancers and is associated with low patient survival. In addition, LTBR expression was strongly correlated with immune cells, score, cancer-related functional status, tumor stemness index, MMR genes, DNA methyltransferase, DNA methylation, mutational status, TMB, and MSI. Enrichment analysis revealed that LTBR was associated with apoptosis, necroptosis, and immune-related pathways. Finally, multiple drugs targeting LTBR were identified. LTBR is overexpressed in several tumors and is associated with a poor prognosis. It is related to immune-related genes and immune cell infiltration.

Circadian rhythm-associated lncRNA RP11-414H17.5 as a key therapeutic target in osteosarcoma affects the tumor immune microenvironment and enhances malignancy

BACKGROUND

It has previously been proven that circadian rhythm disruption is associated with the incidence and deterioration of several tumors, which potentially leads to increased tumor susceptibility and a worse prognosis for tumor-bearing patients. However, their potential role in osteosarcoma has yet to be sufficiently investigated.

METHODS

Transcriptomic and clinical data of 84 osteosarcoma samples and 70 normal bone tissue samples were obtained from the TARGET and GTEx databases, circadian rhythm-related genes were obtained from Genecards, and circadian rhythm-related lncRNAs (CRLs) were obtained by Pearson correlation analysis, differential expression analysis, and protein-protein interaction (PPI) analysis. COX regression and LASSO regression were performed on the CRLs in order to construct a circadian rhythm-related prognostic prediction signature (CRPS). CRPS reliability was verified by Kaplan-Meier (KM), principal component analysis (PCA), nomogram, and receiver operating characteristic (ROC) curve. CRPS effects on the immune microenvironment of osteosarcoma were explored by enrichment analysis and immune infiltration analysis, and the effect of critical gene RP11-414H17.5 on osteosarcoma was experimentally verified.

RESULT

CRPS consisting of three CRLs was constructed and its area under the curve (AUC) values predicted that osteosarcoma prognosis reached 0.892 in the training group and 0.843 in the test group, with a p value of < 0.05 for the KM curve and stable performance across different clinical subgroups. PCA analysis found that CRPS could significantly distinguish between different risk subgroups, and exhibited excellent performance in the prediction of the immune microenvironment. The experiment verified that RP11-414H17.5 can promote metastasis and inhibit apoptosis of osteosarcoma cells.

CONCLUSION

The study revealed that circadian rhythm plays a crucial role in osteosarcoma progression and identified the impact of the key gene RP11-414H17.5 on osteosarcoma, which provides novel insights into osteosarcoma diagnosis and therapy.

Construction of a diagnostic model for hepatitis B-related hepatocellular carcinoma using machine learning and artificial neural networks and revealing the correlation by immunoassay

HBV infection profoundly escalates hepatocellular carcinoma (HCC) susceptibility, responsible for a majority of HCC cases. HBV-driven immune-mediated hepatocyte impairment significantly fuels HCC progression. Regrettably, inconspicuous early HCC symptoms often culminate in belated diagnoses. Nevertheless, surgically treated early-stage HCC patients relish augmented five-year survival rates. In contrast, advanced HCC exhibits feeble responses to conventional interventions like radiotherapy, chemotherapy, and surgery, leading to diminished survival rates. This investigation endeavors to unearth diagnostic hallmark genes for HBV-HCC leveraging a bioinformatics framework, thus refining early HBV-HCC detection. Candidate genes were sieved via differential analysis and Weighted Gene Co-Expression Network Analysis (WGCNA). Employing three distinct machine learning algorithms unearthed three feature genes (HHIP, CXCL14, and CDHR2). Melding these genes yielded an innovative Artificial Neural Network (ANN) diagnostic blueprint, portending to alleviate patient encumbrance and elevate life quality. Immunoassay scrutiny unveiled accentuated immune damage in HBV-HCC patients relative to solitary HCC. Through consensus clustering, HBV-HCC was stratified into two subtypes (C1 and C2), the latter potentially indicating milder immune impairment. The diagnostic model grounded in these feature genes showcased robust and transferrable prognostic potentialities, introducing a novel outlook for early HBV-HCC diagnosis. This exhaustive immunological odyssey stands poised to expedite immunotherapeutic curatives' emergence for HBV-HCC.

Identification of M5c regulator-medicated methylation modification patterns for prognosis and immune microenvironment in glioma

Glioma is a common intracranial tumor and is generally associated with poor prognosis. Recently, numerous studies illustrated the importance of 5-methylcytosine (m5C) RNA modification to tumorigenesis. However, the prognostic value and immune correlation of m5C in glioma remain unclear. We obtained RNA expression and clinical information from The Cancer Genome Atlas (TCGA) and The Chinese Glioma Genome Atlas (CGGA) datasets to analyze. Nonnegative matrix factorization (NMF) was used to classify patients into two subgroups and compare these patients in survival and clinicopathological characteristics. CIBERSORT and single-sample gene-set algorithm (ssGSEA) methods were used to investigate the relationship between m5C and the immune environment. The Weighted correlation network analysis (WGCNA) and univariate Cox proportional hazard model (CoxPH) were used to construct a m5C-related signature. Most of m5C RNA methylation regulators presented differential expression and prognostic values. There were obvious relationships between immune infiltration cells and m5C regulators, especially NSUN7. In the m5C-related module from WGCNA, we found SEPT3, CHI3L1, PLBD1, PHYHIPL, SAMD8, RAP1B, B3GNT5, RER1, PTPN7, SLC39A1, and MXI1 were prognostic factors for glioma, and they were used to construct the signature. The great significance of m5C-related signature in predicting the survival of patients with glioma was confirmed in the validation sets and CGGA cohort.

Machine learning-enhanced insights into sphingolipid-based prognostication: revealing the immunological landscape and predictive proficiency for immunomotherapy and chemotherapy responses in pancreatic carcinoma

Background: With a poor prognosis for affected individuals, pancreatic adenocarcinoma (PAAD) is known as a complicated and diverse illness. Immunocytes have become essential elements in the development of PAAD. Notably, sphingolipid metabolism has a dual function in the development of tumors and the invasion of the immune system. Despite these implications, research on the predictive ability of sphingolipid variables for PAAD prognosis is strikingly lacking, and it is yet unclear how they can affect PAAD immunotherapy and targeted pharmacotherapy. Methods: The investigation process included SPG detection while also being pertinent to the prognosis for PAAD. Both the analytical capability of CIBERSORT and the prognostic capability of the pRRophetic R package were used to evaluate the immunological environments of the various HCC subtypes. In addition, CCK-8 experiments on PAAD cell lines were carried out to confirm the accuracy of drug sensitivity estimates. The results of these trials, which also evaluated cell survival and migratory patterns, confirmed the usefulness of sphingolipid-associated genes (SPGs). Results: As a result of this thorough investigation, 32 SPGs were identified, each of which had a measurable influence on the dynamics of overall survival. This collection of genes served as the conceptual framework for the development of a prognostic model, which was carefully assembled from 10 chosen genes. It should be noted that this grouping of patients into cohorts with high and low risk was a sign of different immune profiles and therapy responses. The increased abundance of SPGs was identified as a possible sign of inadequate responses to immune-based treatment approaches. The careful CCK-8 testing carried out on PAAD cell lines was of the highest importance for providing clear confirmation of drug sensitivity estimates. Conclusion: The significance of Sphingolipid metabolism in the complex web of PAAD development is brought home by this study. The novel risk model, built on the complexity of sphingolipid-associated genes, advances our understanding of PAAD and offers doctors a powerful tool for developing personalised treatment plans that are specifically suited to the unique characteristics of each patient.

Identifying endoplasmic reticulum stress-related molecular subtypes and prognostic model for predicting the immune landscape and therapy response in pancreatic cancer

Endoplasmic reticulum stress (ERS) is caused by the accumulation of intracellular misfolded or unfolded proteins and is associated with cancer development. In this study, pan-cancer analysis revealed complex genetic variations, including copy number variation, methylation, and somatic mutations for ERS-related genes (ERGs) in 33 kinds of cancer. Consensus clustering divided pancreatic cancer (PC) patients from TCGA and GEO databases into two ERS-related subtypes: ERGcluster A and B. Compared with ERGcluster A, ERGcluster B had a more active ERS state and worse prognosis. Subsequently, the ERS-related prognostic model was established to quantify the ERS score for a single sample. The patient with a low ERS score had remarkably longer survival times. ssGSEA and CIBERSORT algorithms revealed that activated B cells and CD8+ T cells had higher infiltration in the low ERS score group, but higher infiltration of activated CD4+ T cells, activated dendritic cells, macrophages, and neutrophils in the high ERS score group. Drug sensitivity analysis indicated the low ERS score group had a better response to gemcitabine, paclitaxel, 5-fluorouracil, oxaliplatin, and irinotecan. RT-qPCR validated that MET, MUC16, and KRT7 in the model had higher expression levels in pancreatic tumour tissues. Single-cell analysis further revealed that MET, MUC16, and KRT7 were mainly expressed in cancer cells in PC tumour microenvironment. In all, we first constructed the ERS-related molecular subtypes and prognostic model in PC. Our research highlighted the vital role of ERS in PC and contributed to further research on molecular mechanisms and novel therapeutic strategies for PC in the future.

Deciphering the role of HPV-mediated metabolic regulation in shaping the tumor microenvironment and its implications for immunotherapy in HNSCC

Head and neck squamous cell carcinoma (HNSCC), as a complex and variable malignancy, poses a significant threat to human health. Since the intricate association between HPV and HNSCC emerged, its role within the TME has garnered extensive attention. HPV+HNSCC exhibits distinct immunological characteristics within the TME, intricately intertwined with mechanisms of immune evasion. HPV employs multifaceted pathways to intervene in metabolic regulation within the TME, exerting influence over immune cell functionality and neoplastic cell genesis. Furthermore, the heightened immune reactivity exhibited by HPV+HNSCC within the TME augments responses to immune interventions such as immune checkpoint inhibitors. Therefore, amidst the current limitations of therapeutic approaches, immunotherapy stands as a promising strategy to overcome the conventional confines of treating HNSCC. This article comprehensively outlines the impact of HPV on the inception and progression of HNSCC while discussing the amalgamation of metabolic regulation within the TME and immunotherapeutic strategies. By intervening in the reciprocal interactions between HPV and HNSCC within the TME, the potential to modulate the efficacy of immune-based treatments becomes evident. Concurrently, a synthesis of pertinent biomarker development is summarized. Such endeavors hold paramount significance for personalized therapeutic approaches and the more effective management of HNSCC patients.

Constructing a prognostic model for head and neck squamous cell carcinoma based on glucose metabolism related genes

Background

Glucose metabolism (GM) plays a crucial role in cancer cell proliferation, tumor growth, and survival. However, the identification of glucose metabolism-related genes (GMRGs) for effective prediction of prognosis in head and neck squamous cell carcinoma (HNSC) is still lacking.

Methods

We conducted differential analysis between HNSC and Normal groups to identify differentially expressed genes (DEGs). Key module genes were obtained using weighted gene co-expression network analysis (WGCNA). Intersection analysis of DEGs, GMRGs, and key module genes identified GMRG-DEGs. Univariate and multivariate Cox regression analyses were performed to screen prognostic-associated genes. Independent prognostic analysis of clinical traits and risk scores was implemented using Cox regression. Gene set enrichment analysis (GSEA) was used to explore functional pathways and genes between high- and low-risk groups. Immune infiltration analysis compared immune cells between the two groups in HNSC samples. Drug prediction was performed using the Genomics of Drug Sensitivity in Cancer (GDSC) database. Quantitative real-time fluorescence PCR (qRT-PCR) validated the expression levels of prognosis-related genes in HNSC patients.

Results

We identified 4973 DEGs between HNSC and Normal samples. Key gene modules, represented by black and brown module genes, were identified. Intersection analysis revealed 76 GMRG-DEGs. Five prognosis-related genes (MTHFD2, CDKN2A, TPM2, MPZ, and DNMT1) were identified. A nomogram incorporating age, lymph node status (N), and risk score was constructed for survival prediction in HNSC patients. Immune infiltration analysis showed significant differences in five immune cell types (Macrophages M0, memory B cells, Monocytes, Macrophages M2, and Dendritic resting cells) between the high- and low-risk groups. GDSC database analysis identified 53 drugs with remarkable differences between the groups, including A.443654 and AG.014699. DNMT1 and MTHFD2 were up-regulated, while MPZ was down-regulated in HNSC.

Conclusion

Our study highlights the significant association of five prognosis-related genes (MTHFD2, CDKN2A, TPM2, MPZ, and DNMT1) with HNSC. These findings provide further evidence of the crucial role of GMRGs in HNSC.

The integrated single-cell analysis developed an immunogenic cell death signature to predict lung adenocarcinoma prognosis and immunotherapy

BACKGROUND

Research on immunogenic cell death (ICD) in lung adenocarcinoma (LUAD) has been relatively limited. This study aims to create ICD-related signatures for accurate survival prognosis prediction in LUAD patients, addressing the challenge of lacking reliable early prognostic indicators for this type of cancer.

METHODS

Using single-cell RNA sequencing (scRNA-seq) analysis, ICD activity in cells was calculated by AUCell algorithm, divided into high- and low-ICD groups according to median values, and key ICD regulatory genes were identified through differential analysis, and these genes were integrated into TCGA data to construct prognostic signatures using LASSO and COX regression analysis, and multi-dimensional analysis of ICD-related signatures in terms of prognosis, immunotherapy, tumor microenvironment (TME), and mutational landscape.

RESULTS

The constructed signature reveals a pronounced disparity in prognosis between the high- and low-risk groups of LUAD patients. The statistical discrepancies in survival times among LUAD patients from both the TCGA and GEO databases further corroborate this observation. Additionally, heightened levels of immune cell infiltration expression are evidenced in the low-risk group, suggesting a potential benefit from immunotherapeutic interventions for these patients. The expression levels of pivotal risk-associated genes in tissue samples were assessed utilizing qRT-PCR, thereby unveiling PITX3 as a plausible therapeutic target in the context of LUAD.

CONCLUSIONS

Our constructed ICD-related signatures provide help in predicting the prognosis and immunotherapy of LUAD patients, and to some extent guide the clinical treatment of LUAD patients.

Integrated single-cell and bulk RNA sequencing revealed the molecular characteristics and prognostic roles of neutrophils in pancreatic cancer

Pancreatic cancer, one of the most prevalent tumors of the digestive system, has a dismal prognosis. Cancer of the pancreas is distinguished by an inflammatory tumor microenvironment rich in fibroblasts and different immune cells. Neutrophils are important immune cells that infiltrate the microenvironment of pancreatic cancer tumors. The purpose of this work was to examine the complex mechanism by which neutrophils influence the carcinogenesis and development of pancreatic cancer and to construct a survival prediction model based on neutrophil marker genes. We incorporated the GSE111672 dataset, comprising RNA expression data from 27,000 cells obtained from 3 patients with PC, and conducted single-cell data analysis. Thorough investigation of pancreatic cancer single-cell RNA sequencing data found 350 neutrophil marker genes. Using The Cancer Genome Atlas (TCGA), GSE28735, GSE62452, GSE57495, and GSE85916 datasets to gather pancreatic cancer tissue transcriptome data, and consistent clustering was used to identify two categories for analyzing the influence of neutrophils on pancreatic cancer. Using the Random Forest algorithm and Cox regression analysis, a survival prediction model for pancreatic cancer was developed, the model showed independent performance for survival prognosis, clinic pathological features, immune infiltration, and drug sensitivity. Multivariate Cox analysis findings revealed that the risk scores derived from predictive models is independent prognostic markers for pancreatic patients. In conclusion, based on neutrophil marker genes, this research created a molecular typing and prognostic grading system for pancreatic cancer, this system was very accurate in predicting the prognosis, tumor immune microenvironment status, and pharmacological treatment responsiveness of pancreatic cancer patients.

Reprogramming of Treg cells in the inflammatory microenvironment during immunotherapy: a literature review

Regulatory T cells (Treg), as members of CD4+ T cells, have garnered extensive attention in the research of tumor progression. Treg cells have the function of inhibiting the immune effector cells, preventing tissue damage, and suppressing inflammation. Under the stimulation of the tumor inflammatory microenvironment (IM), the reprogramming of Treg cells enhances their suppression of immune responses, ultimately promoting tumor immune escape or tumor progression. Reducing the number of Treg cells in the IM or lowering the activity of Treg cells while preventing their reprogramming, can help promote the body's anti-tumor immune responses. This review introduces a reprogramming mechanism of Treg cells in the IM; and discusses the regulation of Treg cells on tumor progression. The control of Treg cells and the response to Treg inflammatory reprogramming in tumor immunotherapy are analyzed and countermeasures are proposed. This work will provide a foundation for downregulating the immunosuppressive role of Treg in the inflammatory environment in future tumor immunotherapy.

Unveiling the role of regulatory T cells in the tumor microenvironment of pancreatic cancer through single-cell transcriptomics and in vitro experiments

Background

In order to investigate the impact of Treg cell infiltration on the immune response against pancreatic cancer within the tumor microenvironment (TME), and identify crucial mRNA markers associated with Treg cells in pancreatic cancer, our study aims to delve into the role of Treg cells in the anti-tumor immune response of pancreatic cancer.

Methods

The ordinary transcriptome data for this study was sourced from the GEO and TCGA databases. It was analyzed using single-cell sequencing analysis and machine learning. To assess the infiltration level of Treg cells in pancreatic cancer tissues, we employed the CIBERSORT method. The identification of genes most closely associated with Treg cells was accomplished through the implementation of weighted gene co-expression network analysis (WGCNA). Our analysis of single-cell sequencing data involved various quality control methods, followed by annotation and advanced analyses such as cell trajectory analysis and cell communication analysis to elucidate the role of Treg cells within the pancreatic cancer microenvironment. Additionally, we categorized the Treg cells into two subsets: Treg1 associated with favorable prognosis, and Treg2 associated with poor prognosis, based on the enrichment scores of the key genes. Employing the hdWGCNA method, we analyzed these two subsets to identify the critical signaling pathways governing their mutual transformation. Finally, we conducted PCR and immunofluorescence staining in vitro to validate the identified key genes.

Results

Based on the results of immune infiltration analysis, we observed significant infiltration of Treg cells in the pancreatic cancer microenvironment. Subsequently, utilizing the WGCNA and machine learning algorithms, we ultimately identified four Treg cell-related genes (TRGs), among which four genes exhibited significant correlations with the occurrence and progression of pancreatic cancer. Among them, CASP4, TOB1, and CLEC2B were associated with poorer prognosis in pancreatic cancer patients, while FYN showed a correlation with better prognosis. Notably, significant differences were found in the HIF-1 signaling pathway between Treg1 and Treg2 cells identified by the four genes. These conclusions were further validated through in vitro experiments.

Conclusion

Treg cells played a crucial role in the pancreatic cancer microenvironment, and their presence held a dual significance. Recognizing this characteristic was vital for understanding the limitations of Treg cell-targeted therapies. CASP4, FYN, TOB1, and CLEC2B exhibited close associations with infiltrating Treg cells in pancreatic cancer, suggesting their involvement in Treg cell functions. Further investigation was warranted to uncover the mechanisms underlying these associations. Notably, the HIF-1 signaling pathway emerged as a significant pathway contributing to the duality of Treg cells. Targeting this pathway could potentially revolutionize the existing treatment approaches for pancreatic cancer.

Integrated multi-omics identified the novel intratumor microbiome-derived subtypes and signature to predict the outcome, tumor microenvironment heterogeneity, and immunotherapy response for pancreatic cancer patients

Background: The extremely malignant tumour known as pancreatic cancer (PC) lacks efficient prognostic markers and treatment strategies. The microbiome is crucial to how cancer develops and responds to treatment. Our study was conducted in order to better understand how PC patients' microbiomes influence their outcome, tumour microenvironment, and responsiveness to immunotherapy. Methods: We integrated transcriptome and microbiome data of PC and used univariable Cox regression and Kaplan-Meier method for screening the prognostic microbes. Then intratumor microbiome-derived subtypes were identified using consensus clustering. We utilized LASSO and Cox regression to build the microbe-related model for predicting the prognosis of PC, and utilized eight algorithms to assess the immune microenvironment feature. The OncoPredict package was utilized to predict drug treatment response. We utilized qRT-PCR to verify gene expression and single-cell analysis to reveal the composition of PC tumour microenvironment. Results: We obtained a total of 26 prognostic genera in PC. And PC samples were divided into two microbiome-related subtypes: Mcluster A and B. Compared with Mcluster A, patients in Mcluster B had a worse prognosis and higher TNM stage and pathological grade. Immune analysis revealed that neutrophils, regulatory T cell, CD8+ T cell, macrophages M1 and M2, cancer associated fibroblasts, myeloid dendritic cell, and activated mast cell had remarkably higher infiltrated levels within the tumour microenvironment of Mcluster B. Patients in Mcluster A were more likely to benefit from CTLA-4 blockers and were highly sensitive to 5-fluorouracil, cisplatin, gemcitabine, irinotecan, oxaliplatin, and epirubicin. Moreover, we built a microbe-derived model to assess the outcome. The ROC curves showed that the microbe-related model has good predictive performance. The expression of LAMA3 and LIPH was markedly increased within pancreatic tumour tissues and was linked to advanced stage and poor prognosis. Single-cell analysis indicated that besides cancer cells, the tumour microenvironment of PC was also rich in monocytes/macrophages, endothelial cells, and fibroblasts. LIPH and LAMA3 exhibited relatively higher expression in cancer cells and neutrophils. Conclusion: The intratumor microbiome-derived subtypes and signature in PC were first established, and our study provided novel perspectives on PC prognostic indicators and treatment options.

The core role of macrophages in hepatocellular carcinoma: the definition of molecular subtypes and the prognostic risk system

Background: In patients with hepatocellular carcinoma (HCC), the tumor microenvironment (TME) is resistant to immunotherapy because of its specificity. It is meaningful to explore the role of macrophage, which is one of the most abundant immune cells in the TME, in cellular communication and its effect on the prognosis and immunotherapy of HCC. Methods: Dimensionality reduction and clustering of the single-cell RNA-seq data from the GSE149614 dataset were carried out to identify the cellular composition of HCC. CellChat was used to analyze the communication between different cells. The specifically highly expressed genes of macrophages were extracted for univariate Cox regression analysis to obtain prognostic genes for HCC cluster analysis, and the risk system of macrophage-specifically highly expressed genes was developed by random forest analysis and multivariate Cox regression analysis. Prognosis, TME infiltration, potential responses to immunotherapy, and antineoplastic drugs were compared among molecular subtypes and between risk groups. Results: We found that HCC included nine identifiable cell types, of which macrophages had the highest communication intensity with each of the other eight cell types. Of the 179 specifically highly expressed genes of macrophage, 56 were significantly correlated with the prognosis of HCC, which classified HCC into three subtypes, which were reproducible and produced different survival outcomes, TME infiltration, and immunotherapy responses among the subtypes. In the integration of four macrophage-specifically highly expressed genes for the development of a risk system, the risk score was significantly involved in higher immune cell infiltration, poor prognosis, immunotherapy response rate, and sensitivity of six drugs. Conclusion: In this study, through single-cell RNA-seq data, we identified nine cell types, among which macrophage had the highest communication intensity with the rest of the cell types. Based on specifically highly expressed genes of macrophage, we successfully divided HCC patients into three clusters with distinct prognosis, TME, and therapeutic response. Additionally, a risk system was constructed, which provided a potential reference index for the prognostic target and preclinical individualized treatment of HCC.

Multi-omics analyses based on genes associated with oxidative stress and phospholipid metabolism revealed the intrinsic molecular characteristics of pancreatic cancer

Oxidative stress (OS), which impacts lipid metabolic reprogramming, can affect the biological activities of cancer cells. How oxidative stress and phospholipid metabolism (OSPM) influence the prognosis of pancreatic cancer (PC) needs to be elucidated. The metabolic data of 35 pancreatic tumor samples, 34 para-carcinoma samples, and 31 normal pancreatic tissues were obtained from the previously published literature. Pan-cancer samples were obtained from The Cancer Genome Atlas (TCGA). And the Gene Expression Omnibus (GEO), International Cancer Genome Consortium (ICGC), ArrayExpress, and the Genotype-Tissue Expression (GTEx) databases were searched for more PC and normal pancreatic samples. The metabolites in PC were compared with normal and para-carcinoma tissues. The characteristics of the key OSPM genes were summarized in pan-cancer. The random survival forest analysis and multivariate Cox regression analysis were utilized to construct an OSPM-related signature. Based on this signature, PC samples were divided into high- and low-risk subgroups. The dysregulations of the tumor immune microenvironment were further investigated. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was conducted to investigate the expression of genes in the signature in PC and normal tissues. The protein levels of these genes were further demonstrated. In PC, metabolomic studies revealed the alteration of PM, while transcriptomic studies showed different expressions of OSPM-related genes. Then 930 PC samples were divided into three subtypes with different prognoses, and an OSPM-related signature including eight OSPM-related genes (i.e., SLC2A1, MMP14, TOP2A, MBOAT2, ANLN, ECT2, SLC22A3, and FGD6) was developed. High- and low-risk subgroups divided by the signature showed different prognoses, expression levels of immune checkpoint genes, immune cell infiltration, and tumor microenvironment. The risk score was negatively correlated with the proportion of TIL, pDC, Mast cell, and T cell co-stimulation. The expression levels of genes in the signature were verified in PC and normal samples. The protein levels of SLC2A1, MMP14, TOP2A, MBOAT2, ANLN, and SLC22A3 showed up-regulation in PC samples compared with normal tissues. After integrating metabolomics and transcriptomics data, the alterations in OSPM in PC were investigated, and an OSPM-related signature was developed, which was helpful for the prognostic assessment and individualized treatment for PC.

Revealing the role of regulatory T cells in the tumor microenvironment of lung adenocarcinoma: a novel prognostic and immunotherapeutic signature

Background

Regulatory T cells (Tregs), are a key class of cell types in the immune system. In the tumor microenvironment (TME), the presence of Tregs has important implications for immune response and tumor development. Relatively little is known about the role of Tregs in lung adenocarcinoma (LUAD).

Methods

Tregs were identified using but single-cell RNA sequencing (scRNA-seq) analysis and interactions between Tregs and other cells in the TME were investigated. Next, we used multiple bulk RNA-seq datasets to construct risk models based on marker genes of Tregs and explored differences in prognosis, mutational landscape, immune cell infiltration and immunotherapy between high- and low-risk groups, and finally, qRT-PCR and cell function experiments were performed to validate the model genes.

Results

The cellchat analysis showed that MIF-(CD74+CXCR4) pairs play a key role in the interaction of Tregs with other cell subpopulations, and the Tregs-associated signatures (TRAS) could well classify multiple LUAD cohorts into high- and low-risk groups. Immunotherapy may offer greater potential benefits to the low-risk group, as indicated by their superior survival, increased infiltration of immune cells, and heightened expression of immune checkpoints. Finally, the experiment verified that the model genes LTB and PTTG1 were relatively highly expressed in cancer tissues, while PTPRC was relatively highly expressed in paracancerous tissues. Colony Formation assay confirmed that knockdown of PTTG1 reduced the proliferation ability of LUAD cells.

Conclusion

TRAS were constructed using scRNA-seq and bulk RNA-seq to distinguish patient risk subgroups, which may provide assistance in the clinical management of LUAD patients.

The combined signatures of telomere and immune cell landscape provide a prognostic and therapeutic biomarker in glioma

Background

Gliomas, the most prevalent primary malignant tumors of the central nervous system in adults, exhibit slow growth in lower-grade gliomas (LGG). However, the majority of LGG cases progress to high-grade gliomas, posing challenges for prognostication. The tumor microenvironment (TME), characterized by telomere-related genes and immune cell infiltration, strongly influences glioma growth and therapeutic response. Therefore, our objective was to develop a Telomere-TME (TM-TME) classifier that integrates telomere-related genes and immune cell landscape to assess prognosis and therapeutic response in glioma.

Methods

This study encompassed LGG patients from the TCGA and CCGA databases. TM score and TME score were derived from the expression signatures of telomere-related genes and the presence of immune cells in LGG, respectively. The TM-TME classifier was established by combining TM and TME scores to effectively predict prognosis. Subsequently, we conducted Kaplan-Meier survival estimation, univariate Cox regression analysis, and receiver operating characteristic curves to validate the prognostic prediction capacity of the TM-TME classifier across multiple cohorts. Gene Ontology (GO) analysis, biological processes, and proteomaps were performed to annotate the functional aspects of each subgroup and visualize the cellular signaling pathways.

Results

The TM_low+TME_high subgroup exhibited superior prognosis and therapeutic response compared to other subgroups (P<0.001). This finding could be attributed to distinct tumor somatic mutations and cancer cellular signaling pathways. GO analysis indicated that the TM_low+TME_high subgroup is associated with the neuronal system and modulation of chemical synaptic transmission. Conversely, the TM_high+TME_low subgroup showed a strong association with cell cycle and DNA metabolic processes. Furthermore, the classifier significantly differentiated overall survival in the TCGA LGG cohort and served as an independent prognostic factor for LGG patients in both the TCGA cohort (P<0.001) and the CGGA cohort (P<0.001).

Conclusion

Overall, our findings underscore the significance of the TM-TME classifier in predicting prognosis and immune therapeutic response in glioma, shedding light on the complex immune landscape within each subgroup. Additionally, our results suggest the potential of integrating risk stratification with precision therapy for LGG.

Single-cell and bulk RNA sequencing reveal cancer-associated fibroblast heterogeneity and a prognostic signature in prostate cancer

Cancer-associated fibroblasts (CAFs), the central players in the tumor microenvironment (TME), can promote tumor progression and metastasis via various functions. However, the properties of CAFs in prostate cancer (PCa) have not been fully assessed. Therefore, we aimed to examine the CAF characteristics in PCa and construct a CAF-derived signature to predict PCa prognosis. CAFs were identified using single-cell RNA sequencing (scRNA-seq) data from 3 studies. We performed the FindAllMarkers function to extract CAF marker genes and constructed a signature to predict the biochemical relapse-free survival (bRFS) of PCa in the Cancer Genome Atlas (TCGA) cohort. Subsequently, different algorithms were applied to reveal the differences of the TME, immune infiltration, treatment responses in the high- and low-risk groups. Additionally, the CAF heterogeneity was assessed in PCa, which were confirmed by the functional enrichment analysis, gene set enrichment analysis (GSEA), and AUCell method. The scRNA-seq analysis identified a CAF cluster with 783 cells and determined 183 CAF marker genes. Cell-cell communication revealed extensive interactions between fibroblasts and immune cells. A CAF-related prognostic model, containing 7 genes (ASPN, AEBP1, ALDH1A1, BGN, COL1A1, PAGE4 and RASD1), was developed to predict bRFS and validated by 4 independent bulk RNA-seq cohorts. Moreover, the high-risk group of the signature score connected with an immunosuppressive TME, such as a higher level of M2 macrophages and lower levels of plasma cells and CD8+ T cells, and a reduced reaction rate for immunotherapy compared with low-risk group. After re-clustering CAFs via unsupervised clustering, we revealed 3 biologically distinct CAF subsets, namely myofibroblast-like CAFs (myCAFs), immune and inflammatory CAFs (iCAFs) and antigen-presenting CAFs (apCAFs). In conclusion, the CAF-derived signature, the first of its kind, can effectively predict PCa prognosis and serve as an indicator for immunotherapy. Furthermore, our study identified 3 CAF subpopulations with distinct functions in PCa.

GADD45B regulates the carcinogenesis process of chronic atrophic gastritis and the metabolic pathways of gastric cancer

Background

Gastric cancer continues to be a significant global healthcare challenge, and its burden remains substantial. The development of gastric cancer (GC) is closely linked to chronic atrophic gastritis (CAG), yet there is a scarcity of research exploring the underlying mechanisms of CAG-induced carcinogenesis.

Methods

In this study, we conducted a comprehensive investigation into the oncogenes involved in CAG using both bulk transcriptome and single-cell transcriptome data. Our approach employed hdWGCNA to identify pathogenic genes specific to CAG, with non-atrophic gastritis (NAG) serving as the control group. Additionally, we compared CAG with GC, using normal gastric tissue as the control group in the single-cell transcriptome analysis. By intersecting the identified pathogenic genes, we pinpointed key network molecules through protein interaction network analysis. To further refine the gene selection, we applied LASSO, SVM-RFE, and RF techniques, which resulted in a set of cancer-related genes (CRGs) associated with CAG. To identify CRGs potentially linked to gastric cancer progression, we performed a univariate COX regression analysis on the gene set. Subsequently, we explored the relationship between CRGs and immune infiltration, drug sensitivity, and clinical characteristics in gastric cancer patients. We employed GSVA to investigate how CRGs regulated signaling pathways in gastric cancer cells, while an analysis of cell communication shed light on the impact of CRGs on signal transmission within the gastric cancer tumor microenvironment. Lastly, we analyzed changes in metabolic pathways throughout the progression of gastric cancer.

Results

Using hdWGCNA, we have identified a total of 143 pathogenic genes that were shared by CAG and GC. To further investigate the underlying mechanisms, we conducted protein interaction network analysis and employed machine learning screening techniques. As a result, we have identified 15 oncogenes that are specifically associated with chronic atrophic gastritis. By performing ROC reanalysis and prognostic analysis, we have determined that GADD45B is the most significant gene involved in the carcinogenesis of CAG. Immunohistochemical staining and differential analysis have revealed that GADD45B expression was low in GC tissues while high in normal gastric tissues. Moreover, based on prognostic analysis, high expression of GADD45B has been correlated with poor prognosis in GC patients. Additionally, an analysis of immune infiltration has shown a relationship between GADD45B and the infiltration of various immune cells. By correlating GADD45B with clinical characteristics, we have found that it primarily affects the depth of invasion in GC. Through cell communication analysis, we have discovered that the CD99 signaling pathway network and the CDH signaling pathway network are the main communication pathways that significantly alter the microenvironment of gastric tissue during the development of chronic atrophic gastritis. Specifically, GADD45B-low GC cells were predominantly involved in the network communication of the CDH signaling pathway, while GADD45B-high GC cells played a crucial role in both signaling pathways. Furthermore, we have identified several metabolic pathways, including D-Glutamine and D-glutamate metabolism and N-Glycan biosynthesis, among others, that played important roles in the occurrence and progression of GC, in addition to the six other metabolic pathways. In summary, our study highlighted the discovery of 143 pathogenic genes shared by CAG and GC, with a specific focus on 15 oncogenes associated with CAG. We have identified GADD45B as the most important gene in the carcinogenesis of CAG, which exhibited differential expression in GC tissues compared to normal gastric tissues. Moreover, GADD45B expression was correlated with patient prognosis and is associated with immune cell infiltration. Our findings also emphasized the impact of the CD99 and CDH signaling pathway networks on the microenvironment of gastric tissue during the development of CAG. Additionally, we have identified key metabolic pathways involved in GC progression.

Conclusion

GADD45B, an oncogene implicated in chronic atrophic gastritis, played a critical role in GC development. Decreased expression of GADD45B was associated with the onset of GC. Moreover, GADD45B expression levels were closely tied to poor prognosis in GC patients, influencing the infiltration patterns of various cells within the tumor microenvironment, as well as impacting the metabolic pathways involved in GC progression.

Proposing new early detection indicators for pancreatic cancer: Combining machine learning and neural networks for serum miRNA-based diagnostic model

Background

Pancreatic cancer (PC) is a lethal malignancy that ranks seventh in terms of global cancer-related mortality. Despite advancements in treatment, the five-year survival rate remains low, emphasizing the urgent need for reliable early detection methods. MicroRNAs (miRNAs), a group of non-coding RNAs involved in critical gene regulatory mechanisms, have garnered significant attention as potential diagnostic and prognostic biomarkers for pancreatic cancer (PC). Their suitability stems from their accessibility and stability in blood, making them particularly appealing for clinical applications.

Methods

In this study, we analyzed serum miRNA expression profiles from three independent PC datasets obtained from the Gene Expression Omnibus (GEO) database. To identify serum miRNAs associated with PC incidence, we employed three machine learning algorithms: Support Vector Machine-Recursive Feature Elimination (SVM-RFE), Least Absolute Shrinkage and Selection Operator (LASSO), and Random Forest. We developed an artificial neural network model to assess the accuracy of the identified PC-related serum miRNAs (PCRSMs) and create a nomogram. These findings were further validated through qPCR experiments. Additionally, patient samples with PC were classified using the consensus clustering method.

Results

Our analysis revealed three PCRSMs, namely hsa-miR-4648, hsa-miR-125b-1-3p, and hsa-miR-3201, using the three machine learning algorithms. The artificial neural network model demonstrated high accuracy in distinguishing between normal and pancreatic cancer samples, with verification and training groups exhibiting AUC values of 0.935 and 0.926, respectively. We also utilized the consensus clustering method to classify PC samples into two optimal subtypes. Furthermore, our investigation into the expression of PCRSMs unveiled a significant negative correlation between the expression of hsa-miR-125b-1-3p and age.

Conclusion

Our study introduces a novel artificial neural network model for early diagnosis of pancreatic cancer, carrying significant clinical implications. Furthermore, our findings provide valuable insights into the pathogenesis of pancreatic cancer and offer potential avenues for drug screening, personalized treatment, and immunotherapy against this lethal disease.

Unveiling efferocytosis-related signatures through the integration of single-cell analysis and machine learning: a predictive framework for prognosis and immunotherapy response in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) represents a prominent gastrointestinal malignancy with a grim clinical outlook. In this regard, the discovery of novel early biomarkers holds substantial promise for ameliorating HCC-associated mortality. Efferocytosis, a vital immunological process, assumes a central position in the elimination of apoptotic cells. However, comprehensive investigations exploring the role of efferocytosis-related genes (EFRGs) in HCC are sparse, and their regulatory influence on HCC immunotherapy and targeted drug interventions remain poorly understood.

Methods

RNA sequencing data and clinical characteristics of HCC patients were acquired from the TCGA database. To identify prognostically significant genes in HCC, we performed the limma package and conducted univariate Cox regression analysis. Subsequently, machine learning algorithms were employed to identify hub genes. To assess the immunological landscape of different HCC subtypes, we employed the CIBERSORT algorithm. Furthermore, single-cell RNA sequencing (scRNA-seq) was utilized to investigate the expression levels of ERFGs in immune cells and to explore intercellular communication within HCC tissues. The migratory capacity of HCC cells was evaluated using CCK-8 assays, while drug sensitivity prediction reliability was determined through wound-healing assays.

Results

We have successfully identified a set of nine genes, termed EFRGs, that hold significant potential for the establishment of a hepatocellular carcinoma-specific prognostic model. Furthermore, leveraging the individual risk scores derived from this model, we were able to stratify patients into two distinct risk groups, unveiling notable disparities in terms of immune infiltration patterns and response to immunotherapy. Notably, the model's capacity to accurately predict drug responses was substantiated through comprehensive experimental investigations, encompassing wound-healing assay, and CCK8 experiments conducted on the HepG2 and Huh7 cell lines.

Conclusions

We constructed an EFRGs model that serves as valuable tools for prognostic assessment and decision-making support in the context of immunotherapy and chemotherapy.

PANoptosis-related molecular subtype and prognostic model associated with the immune microenvironment and individualized therapy in pancreatic cancer

Background

PANoptosis is an inflammatory type of programmed cell death regulated by PANopotosome. Mounting evidence has shown that PANoptosis could be involved in cancer pathogenesis and the tumor immune microenvironment. Nevertheless, there have been no studies on the mechanism of PANoptosis on pancreatic cancer (PC) pathogenesis.

Methods

We downloaded the data on transcriptomic and clinical features of PC patients from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus databases. Additionally, the data on copy number variation (CNV), methylation and somatic mutations of genes in 33 types of cancers were obtained from TCGA. Next, we identified the PANoptosis-related molecular subtype using the consensus clustering analysis, and constructed and validated the PANoptosis-related prognostic model using LASSO and Cox regression analyses. Moreover, RT-qPCR was performed to determine the expression of genes involved in the model.

Results

We obtained 66 PANoptosis-related genes (PANRGs) from published studies. Of these, 24 PC-specific prognosis-related genes were identified. Pan-cancer analysis revealed complex genetic changes, including CNV, methylation, and mutation in PANRGs were identified in various cancers. By consensus clustering analysis, PC patients were classified into two PANoptosis-related patterns: PANcluster A and B. In PANcluster A, the patient prognosis was significantly worse compared to PANcluster B. The CIBERSORT algorithm showed a significant increase in the infiltration of CD8⁺ T cells, monocytes, and naïve B cells, in patients in PANcluster B. Additionally, the infiltration of macrophages, activated mast cells, and dendritic cells were higher in patients in PANcluster A. Patients in PANcluster A were more sensitive to erlotinib, selumetinib and trametinib, whereas patients in PANcluster B were highly sensitive to irinotecan, oxaliplatin and sorafenib. Moreover, we constructed and validated the PANoptosis-related prognostic model to predict the patient's survival. Finally, the GEPIA and Human Protein Atlas databases were analyzed, and RT-qPCR was performed. Compared to normal tissues, a significant increase in CXCL10 and ITGB6 (associated with the model) expression was observed in PC tissues.

Conclusion

We first identified the PANoptosis-related molecular subtypes and established a PANoptosis-related prognostic model for predicting the survival of patients with PC. These results would aid in exploring the mechanisms of PANoptosis in PC pathogenesis.

Prognostic model for hepatocellular carcinoma based on anoikis-related genes: immune landscape analysis and prediction of drug sensitivity

Background

Hepatocellular carcinoma (HCC) represents a complex ailment characterized by an unfavorable prognosis in advanced stages. The involvement of immune cells in HCC progression is of significant importance. Moreover, metastasis poses a substantial impediment to enhanced prognostication for HCC patients, with anoikis playing an indispensable role in facilitating the distant metastasis of tumor cells. Nevertheless, limited investigations have been conducted regarding the utilization of anoikis factors for predicting HCC prognosis and assessing immune infiltration. This present study aims to identify hepatocellular carcinoma-associated anoikis-related genes (ANRGs), establish a robust prognostic model for HCC, and delineate distinct immune characteristics based on the anoikis signature. Cell migration and cytotoxicity experiments were performed to validate the accuracy of the ANRGs model.

Methods

Consensus clustering based on ANRGs was employed in this investigation to categorize HCC samples obtained from both TCGA and Gene Expression Omnibus (GEO) cohorts. To assess the differentially expressed genes, Cox regression analysis was conducted, and subsequently, prognostic gene signatures were constructed using LASSO-Cox methodology. External validation was performed at the International Cancer Genome Conference. The tumor microenvironment (TME) was characterized utilizing ESTIMATE and CIBERSORT algorithms, while machine learning techniques facilitated the identification of potential target drugs. The wound healing assay and CCK-8 assay were employed to evaluate the migratory capacity and drug sensitivity of HCC cell lines, respectively.

Results

Utilizing the TCGA-LIHC dataset, we devised a nomogram integrating a ten-gene signature with diverse clinicopathological features. Furthermore, the discriminative potential and clinical utility of the ten-gene signature and nomogram were substantiated through ROC analysis and DCA. Subsequently, we devised a prognostic framework leveraging gene expression data from distinct risk cohorts to predict the drug responsiveness of HCC subtypes.

Conclusion

In this study, we have established a promising HCC prognostic ANRGs model, which can serve as a valuable tool for clinicians in selecting targeted therapeutic drugs, thereby improving overall patient survival rates. Additionally, this model has also revealed a strong connection between anoikis and immune cells, providing a potential avenue for elucidating the mechanisms underlying immune cell infiltration regulated by anoikis.

Advances in HPV-associated tumor management: Therapeutic strategies and emerging insights

With the rapid increase in the incidence of cervical cancer, anal cancer and other cancers, human papillomavirus (HPV) infection has become a growing concern. Persistent infection with high-risk HPV is a major cause of malignant tumors. In addition, microbiota and viruses such as human immunodeficiency virus, herpes simplex virus, and Epstein-Barr virus are closely associated with HPV infection. The limited effectiveness of existing treatments for HPV-associated tumors and the high rates of recurrence and metastasis in patients create an urgent need for novel and effective approaches. In recent years, HPV vaccine coverage has increased and can reduce the incidence of serious adverse events. Overall, this article provides a comprehensive overview of HPV biology, microbiome, and other viral interactions in cancer development, highlighting the need for a more comprehensive approach to cancer prevention and treatment. Current and emerging HPV-related cancer control and treatment strategies are also further explored.

Innovative breakthroughs facilitated by single-cell multi-omics: manipulating natural killer cell functionality correlates with a novel subcategory of melanoma cells

Background

Melanoma is typically regarded as the most dangerous form of skin cancer. Although surgical removal of in situ lesions can be used to effectively treat metastatic disease, this condition is still difficult to cure. Melanoma cells are removed in great part due to the action of natural killer (NK) and T cells on the immune system. Still, not much is known about how the activity of NK cell-related pathways changes in melanoma tissue. Thus, we performed a single-cell multi-omics analysis on human melanoma cells in this study to explore the modulation of NK cell activity.

Materials and methods

Cells in which mitochondrial genes comprised > 20% of the total number of expressed genes were removed. Gene ontology (GO), gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), and AUCcell analysis of differentially expressed genes (DEGs) in melanoma subtypes were performed. The CellChat package was used to predict cell-cell contact between NK cell and melanoma cell subtypes. Monocle program analyzed the pseudotime trajectories of melanoma cells. In addition, CytoTRACE was used to determine the recommended time order of melanoma cells. InferCNV was utilized to calculate the CNV level of melanoma cell subtypes. Python package pySCENIC was used to assess the enrichment of transcription factors and the activity of regulons in melanoma cell subtypes. Furthermore, the cell function experiment was used to confirm the function of TBX21 in both A375 and WM-115 melanoma cell lines.

Results

Following batch effect correction, 26,161 cells were separated into 28 clusters and designated as melanoma cells, neural cells, fibroblasts, endothelial cells, NK cells, CD4+ T cells, CD8+ T cells, B cells, plasma cells, monocytes and macrophages, and dendritic cells. A total of 10137 melanoma cells were further grouped into seven subtypes, i.e., C0 Melanoma BIRC7, C1 Melanoma CDH19, C2 Melanoma EDNRB, C3 Melanoma BIRC5, C4 Melanoma CORO1A, C5 Melanoma MAGEA4, and C6 Melanoma GJB2. The results of AUCell, GSEA, and GSVA suggested that C4 Melanoma CORO1A may be more sensitive to NK and T cells through positive regulation of NK and T cell-mediated immunity, while other subtypes of melanoma may be more resistant to NK cells. This suggests that the intratumor heterogeneity (ITH) of melanoma-induced activity and the difference in NK cell-mediated cytotoxicity may have caused NK cell defects. Transcription factor enrichment analysis indicated that TBX21 was the most important TF in C4 Melanoma CORO1A and was also associated with M1 modules. In vitro experiments further showed that TBX21 knockdown dramatically decreases melanoma cells' proliferation, invasion, and migration.

Conclusion

The differences in NK and T cell-mediated immunity and cytotoxicity between C4 Melanoma CORO1A and other melanoma cell subtypes may offer a new perspective on the ITH of melanoma-induced metastatic activity. In addition, the protective factors of skin melanoma, STAT1, IRF1, and FLI1, may modulate melanoma cell responses to NK or T cells.

Circadian rhythms and cancer

Circadian rhythms are autonomous oscillators developed by the molecular circadian clock, essential for coordinating internal time with the external environment in a 24-h daily cycle. In mammals, this circadian clock system plays a major role in all physiological processes and severely affects human health. The regulation of the circadian clock extends beyond the clock genes to involve several clock-controlled genes. Hence, the aberrant expression of these clock genes leads to the downregulation of important targets that control the cell cycle and the ability to undergo apoptosis. This may lead to genomic instability and promotes carcinogenesis. Alteration in the clock genes and their modulation is recognized as a new approach for the development of effective treatment against several diseases, including cancer. Until now, there has been a lack of understanding of circadian rhythms and cancer disease. For that, this chapter aims to represent the core components of circadian rhythms and their function in cancer pathogenesis and progression. In addition, the clinical impacts, current clock drugs, and potential therapeutic targets have been discussed.

Efferocytosis signatures as prognostic markers for revealing immune landscape and predicting immunotherapy response in hepatocellular carcinoma

Background: Hepatocellular carcinoma (HCC) is a highly lethal liver cancer with late diagnosis; therefore, the identification of new early biomarkers could help reduce mortality. Efferocytosis, a process in which one cell engulfs another cell, including macrophages, dendritic cells, NK cells, etc., plays a complex role in tumorigenesis, sometimes promoting and sometimes inhibiting tumor development. However, the role of efferocytosis-related genes (ERGs) in HCC progression has been poorly studied, and their regulatory effects in HCC immunotherapy and drug targeting have not been reported. Methods: We downloaded efferocytosis-related genes from the Genecards database and screened for ERGs that showed significant expression changes between HCC and normal tissues and were associated with HCC prognosis. Machine learning algorithms were used to study prognostic gene features. CIBERSORT and pRRophetic R packages were used to evaluate the immune environment of HCC subtypes and predict treatment response. CCK-8 experiments conducted on HCC cells were used to assess the reliability of drug sensitivity prediction. Results: We constructed a prognostic prediction model composed of six genes, and the ROC curve showed good predictive accuracy of the risk model. In addition, two ERG-related subgroups in HCC showed significant differences in tumor immune landscape, immune response, and prognostic stratification. The CCK-8 experiment conducted on HCC cells confirmed the reliability of drug sensitivity prediction. Conclusion: Our study emphasizes the importance of efferocytosis in HCC progression. The risk model based on efferocytosis-related genes developed in our study provides a novel precision medicine approach for HCC patients, allowing clinicians to customize treatment plans based on unique patient characteristics. The results of our investigation carry noteworthy implications for the development of individualized treatment approaches involving immunotherapy and chemotherapy, thereby potentially facilitating the realization of personalized and more efficacious therapeutic interventions for HCC.

Multi-omic analysis of glycolytic signatures: exploring the predictive significance of heterogeneity and stemness in immunotherapy response and outcomes in hepatocellular carcinoma

Background: Hepatocellular carcinoma (HCC) is a global health challenge with complex pathophysiology, characterized by high mortality rates and poor early detection due to significant tumor heterogeneity. Stemness significantly contributes to the heterogeneity of HCC tumors, and glycolysis is crucial for maintaining stemness. However, the predictive significance of glycolysis-related metabolic genes (GMGs) in HCC remains unknown. Therefore, this study aimed to identify critical GMGs and establish a reliable model for HCC prognosis. Methods: GMGs associated with prognosis were identified by evaluating genes with notable expression changes between HCC and normal tissues retrieved from the MsigDB database. Prognostic gene characteristics were established using univariate and multivariate Cox regression studies for prognosis prediction and risk stratification. The "CIBERSORT" and "pRRophetic" R packages were respectively used to evaluate the immunological environment and predict treatment response in HCC subtypes. The HCC stemness score was obtained using the OCLR technique. The precision of drug sensitivity prediction was evaluated using CCK-8 experiments performed on HCC cells. The miagration and invasion ability of HCC cell lines with different riskscores were assessed using Transwell and wound healing assays. Results: The risk model based on 10 gene characteristics showed high prediction accuracy as indicated by the receiver operating characteristic (ROC) curves. Moreover, the two GMG-related subgroups showed considerable variation in the risk of HCC with respect to tumor stemness, immune landscape, and prognostic stratification. The in vitro validation of the model's ability to predict medication response further demonstrated its reliability. Conclusion: Our study highlights the importance of stemness variability and inter-individual variation in determining the HCC risk landscape. The risk model we developed provides HCC patients with a novel method for precision medicine that enables clinical doctors to customize treatment plans based on unique patient characteristics. Our findings have significant implications for tailored immunotherapy and chemotherapy methods, and may pave the way for more personalized and effective treatment strategies for HCC.

An amino acid metabolism-based seventeen-gene signature correlates with the clinical outcome and immune features in pancreatic cancer

Background: Pancreatic cancer is an aggressive tumor with a low 5-year survival rate and primary resistance to most therapy. Amino acid (AA) metabolism is highly correlated with tumor growth, crucial to the aggressive biological behavior of pancreatic cancer; nevertheless, the comprehensive predictive significance of genes that regulate AA metabolism in pancreatic cancer remains unknown. Methods: The mRNA expression data downloaded from The Cancer Genome Atlas (TCGA) were derived as the training cohort, and the GSE57495 cohort from Gene Expression Omnibus (GEO) database was applied as the validation cohort. Random survival forest (RSF) and the least absolute shrinkage and selection operator (LASSO) regression analysis were employed to screen genes and construct an AA metabolism-related risk signature (AMRS). Kaplan-Meier analysis and receiver operating characteristic (ROC) curve were performed to assess the prognostic value of AMRS. We performed genomic alteration analysis and explored the difference in tumor microenvironment (TME) landscape associated with KRAS and TP53 mutation in both high- and low-AMRS groups. Subsequently, the relationships between AMRS and immunotherapy and chemotherapy sensitivity were evaluated. Results: A 17-gene AA metabolism-related risk model in the TCGA cohort was constructed according to RSF and LASSO. After stratifying patients into high- and low-AMRS groups based on the optimal cut-off value, we found that high-AMRS patients had worse overall survival (OS) in the training cohort (a median OS: 13.1 months vs. 50.1 months, p < 0.0001) and validation cohort (a median OS: 16.2 vs. 30.5 months, p = 1e-04). Genetic mutation analysis revealed that KRAS and TP53 were significantly more mutated in high-AMRS group, and patients with KRAS and TP53 alterations had significantly higher risk scores than those without. Based on the analysis of TME, low-AMRS group displayed significantly higher immune score and more enrichment of T Cell CD8⁺ cells. In addition, high-AMRS-group exhibited higher TMB and significantly lower tumor immune dysfunction and exclusion (TIDE) score and T Cells dysfunction score, which suggested a higher sensitive to immunotherapy. Moreover, high-AMRS group was also more sensitive to paclitaxel, cisplatin, and docetaxel. Conclusion: Overall, we constructed an AA-metabolism prognostic model, which provided a powerful prognostic predictor for the clinical treatment of pancreatic cancer.

A novel signature predicts prognosis and immunotherapy in lung adenocarcinoma based on cancer-associated fibroblasts

Background

Extensive research has established the significant correlations between cancer-associated fibroblasts (CAFs) and various stages of cancer development, including initiation, angiogenesis, progression, and resistance to therapy. In this study, we aimed to investigate the characteristics of CAFs in lung adenocarcinoma (LUAD) and develop a risk signature to predict the prognosis of patients with LUAD.

Methods

We obtained single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data from the public database. The Seurat R package was used to process the scRNA-seq data and identify CAF clusters based on several biomarkers. CAF-related prognostic genes were further identified using univariate Cox regression analysis. To reduce the number of genes, Lasso regression was performed, and a risk signature was established. A novel nomogram that incorporated the risk signature and clinicopathological features was developed to predict the clinical applicability of the model. Additionally, we conducted immune landscape and immunotherapy responsiveness analyses. Finally, we performed in vitro experiments to verify the functions of EXO1 in LUAD.

Results

We identified 5 CAF clusters in LUAD using scRNA-seq data, of which 3 clusters were significantly associated with prognosis in LUAD. A total of 492 genes were found to be significantly linked to CAF clusters from 1731 DEGs and were used to construct a risk signature. Moreover, our immune landscape exploration revealed that the risk signature was significantly related to immune scores, and its ability to predict responsiveness to immunotherapy was confirmed. Furthermore, a novel nomogram incorporating the risk signature and clinicopathological features showed excellent clinical applicability. Finally, we verified the functions of EXP1 in LUAD through in vitro experiments.

Conclusions

The risk signature has proven to be an excellent predictor of LUAD prognosis, stratifying patients more appropriately and precisely predicting immunotherapy responsiveness. The comprehensive characterization of LUAD based on the CAF signature can predict the response of LUAD to immunotherapy, thus offering fresh perspectives into the management of LUAD patients. Our study ultimately confirms the role of EXP1 in facilitating the invasion and growth of tumor cells in LUAD. Nevertheless, further validation can be achieved by conducting in vivo experiments.

FAM family gene prediction model reveals heterogeneity, stemness and immune microenvironment of UCEC

Background: Endometrial cancer (UCEC) is a highly heterogeneous gynecologic malignancy that exhibits variable prognostic outcomes and responses to immunotherapy. The Familial sequence similarity (FAM) gene family is known to contribute to the pathogenesis of various malignancies, but the extent of their involvement in UCEC has not been systematically studied. This investigation aimed to develop a robust risk profile based on FAM family genes (FFGs) to predict the prognosis and suitability for immunotherapy in UCEC patients. Methods: Using the TCGA-UCEC cohort from The Cancer Genome Atlas (TCGA) database, we obtained expression profiles of FFGs from 552 UCEC and 35 normal samples, and analyzed the expression patterns and prognostic relevance of 363 FAM family genes. The UCEC samples were randomly divided into training and test sets (1:1), and univariate Cox regression analysis and Lasso Cox regression analysis were conducted to identify the differentially expressed genes (FAM13C, FAM110B, and FAM72A) that were significantly associated with prognosis. A prognostic risk scoring system was constructed based on these three gene characteristics using multivariate Cox proportional risk regression. The clinical potential and immune status of FFGs were analyzed using CiberSort, SSGSEA, and tumor immune dysfunction and rejection (TIDE) algorithms. qRT-PCR and IHC for detecting the expression levels of 3-FFGs. Results: Three FFGs, namely, FAM13C, FAM110B, and FAM72A, were identified as strongly associated with the prognosis of UCEC and effective predictors of UCEC prognosis. Multivariate analysis demonstrated that the developed model was an independent predictor of UCEC, and that patients in the low-risk group had better overall survival than those in the high-risk group. The nomogram constructed from clinical characteristics and risk scores exhibited good prognostic power. Patients in the low-risk group exhibited a higher tumor mutational load (TMB) and were more likely to benefit from immunotherapy. Conclusion: This study successfully developed and validated novel biomarkers based on FFGs for predicting the prognosis and immune status of UCEC patients. The identified FFGs can accurately assess the prognosis of UCEC patients and facilitate the identification of specific subgroups of patients who may benefit from personalized treatment with immunotherapy and chemotherapy.

Integrating multiple machine learning methods to construct glutamine metabolism-related signatures in lung adenocarcinoma

Background

Glutamine metabolism (GM) is known to play a critical role in cancer development, including in lung adenocarcinoma (LUAD), although the exact contribution of GM to LUAD remains incompletely understood. In this study, we aimed to discover new targets for the treatment of LUAD patients by using machine learning algorithms to establish prognostic models based on GM-related genes (GMRGs).

Methods

We used the AUCell and WGCNA algorithms, along with single-cell and bulk RNA-seq data, to identify the most prominent GMRGs associated with LUAD. Multiple machine learning algorithms were employed to develop risk models with optimal predictive performance. We validated our models using multiple external datasets and investigated disparities in the tumor microenvironment (TME), mutation landscape, enriched pathways, and response to immunotherapy across various risk groups. Additionally, we conducted in vitro and in vivo experiments to confirm the role of LGALS3 in LUAD.

Results

We identified 173 GMRGs strongly associated with GM activity and selected the Random Survival Forest (RSF) and Supervised Principal Components (SuperPC) methods to develop a prognostic model. Our model's performance was validated using multiple external datasets. Our analysis revealed that the low-risk group had higher immune cell infiltration and increased expression of immune checkpoints, indicating that this group may be more receptive to immunotherapy. Moreover, our experimental results confirmed that LGALS3 promoted the proliferation, invasion, and migration of LUAD cells.

Conclusion

Our study established a prognostic model based on GMRGs that can predict the effectiveness of immunotherapy and provide novel approaches for the treatment of LUAD. Our findings also suggest that LGALS3 may be a potential therapeutic target for LUAD.

Identification and validation of TME-related signatures to predict prognosis and response to anti-tumor therapies in skin cutaneous melanoma

The tumor microenvironment (TME) dynamically regulates cancer progression and affects clinical outcomes. This study aimed to identify molecular subtypes and construct a prognostic risk model based on TME-related signatures in skin cutaneous melanoma (SKCM) patients. We categorized SKCM patients based on transcriptome data of SKCM from The Cancer Genome Atlas (TCGA) database and 29 TME-related gene signatures. Differentially expressed genes were identified using univariate Cox regression and Lasso regression analysis, which were used for risk model construction. The robustness of this model was validated in independent external cohorts. Genetic landscape alterations, immune characteristics, and responsiveness to immunotherapy/chemotherapy were evaluated. Three TME-related subtypes were identified, and subtype C3 exhibited the most favorable prognosis, had enriched immune-related pathways, and possessed more infiltration of T_cells_CD8, T_cells_CD4_memory_activated, and Macrophages_M1 but a lower TumorPurity, whereas Macrophages_M2 were increased in subtype C1 and subtype C2. Subtype C1 was more sensitive to Cisplatin, subtype C2 was more sensitive to Temozolomide, and subtype C3 was more sensitive to Paclitaxel; 8 TME-related genes (NOTCH3, HEYL, ZNF703, ABCC2, PAEP, CCL8, HAPLN3, and HPDL) were screened for risk model construction. High-risk patients had dismal prognosis with good prediction performance. Moreover, low-risk patients were more sensitive to Paclitaxel and Temozolomide, whereas high-risk patients were more sensitive to Cisplatin. This risk model had robustness in predicting prognosis in SKCM patients. The results facilitate the understanding of TME-related genes in SKCM and provide a TME-related genes-based predictive model in prognosis and direction of personalized options for SKCM patients.