The integrated single-cell analysis developed a lactate metabolism-driven signature to improve outcomes and immunotherapy in lung adenocarcinoma

Machine learning algorithms integrate bulk and single-cell RNA data to unveil oxidative stress following intracerebral hemorrhage

BACKGROUND

Increased oxidative stress (OS) activity following intracerebral hemorrhage (ICH) had significantly impacting patient prognosis. Identifying optimal genes associated with OS could enhance the understanding of OS after ICH.

METHODS

We employed single-cell RNA sequencing (scRNA-seq) to investigate the heterogeneity of OS across various cellular tiers following ICH, aiming to acquire biological insights into ICH. We utilized AUCell, Ucell, singscore, ssgsea, and AddModuleScore algorithms, along with correlation analysis, to identify hub genes influencing high OS post-ICH. Furthermore, we employed four machine learning algorithms, eXtreme Gradient Boosting, Boruta, Random Forest, and Least Absolute Shrinkage and Selection Operator, to identify the optimal feature genes. To validate the accuracy of our analysis, we conducted validation in ICH animal experiments.

RESULTS

After analyzing the scRNA-seq dataset using various algorithms, we found that OS activity exhibited heterogeneity across different cellular layers following ICH, with particularly heightened activity observed in monocytes. Further integration of bulk data and machine learning algorithms revealed that ANXA2 and COTL1 were closely associated with high OS after ICH. Our animal experiments demonstrated an increase in OS expression post-ICH. Additionally, the protein expression of ANXA2 and COTL1 was significantly elevated and co-localized with microglia. Pearson correlation coefficient analysis revealed a significant correlation between ANXA2 and OS, indicating strong consistency (r = 0.84, p < 0.05). Similar results were observed for COTL1 and OS (r = 0.69, p < 0.05).

CONCLUSIONS

Following ICH, ANXA2 and COTL1 might penetrate the brain via monocytes, localize within microglia, and enhance OS activity. This might help us better understand OS after ICH.

The prognostic value of immune escape-related genes in lung adenocarcinoma

Background

Lung cancer is one of the most common cancers in humans, and lung adenocarcinoma (LUAD) has become the most common histological type of lung cancer. Immune escape promotes progression of LUAD from the early to metastatic late stages and is one of the main obstacles to improving clinical outcomes for immunotherapy targeting immune detection points. Our study aims to explore the immune escape related genes that are abnormally expressed in lung adenocarcinoma, providing assistance in predicting the prognosis of lung adenocarcinoma and targeted.

Methods

RNA data and related clinical details of patients with LUAD were obtained from The Cancer Genome Atlas (TCGA) database. Through weighted gene coexpression network analysis (WGCNA), 3112 key genes were screened and intersected with 182 immune escape genes obtained from a previous study to identify the immune escape-related genes (IERGs). The role of IERGs in LUAD was systematically explored through gene ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) analyses, which were used to enrich the relevant pathways of IERGs. The least absolute shrinkage and selection operator (LASSO) algorithm and multivariate Cox regression analysis were used to identify the key prognostic genes, and a prognostic risk model was constructed. Estimation of Stromal and Immune Cells in Malignant Tumor Tissues Using Expression Data (ESTIMATE) and microenvironment cell populations (MCP) counter methods (which can accurately assess the amount of eight immune cell populations and two stromal cell groups) were used to analyze the tumor immune status of the high and low risk subgroups. The protein expression level of the differentially expressed genes in lung cancer samples was determined by using the Human Protein Atlas (HPA) database. A nomogram was constructed, and the prognostic risk model was verified via the Gene Expression Omnibus (GEO) datasets GSE72094 and GSE30219.

Results

Twenty differentially expressed IERGs were obtained. GO analysis of these 20 IERGs revealed that they were mainly associated with the regulation of immune system processes, immune responses, and interferon-γ enrichment in mediating signaling pathways and apoptotic signaling pathways; meanwhile, KEGG analysis revealed that IERGs were associated with necroptosis, antigen processing and presentation, programmed cell death ligand 1 (PD-L1) expression and programmed cell death 1 (PD-1) pathway in tumors, cytokine-cytokine receptor interactions, T helper cell 1 (Th1) and Th2 differentiation, and tumor necrosis factor signaling pathways. Using LASSO and Cox regression analysis, we constructed a four-gene model that could predict the prognosis of patients with LUAD, and the model was validated with a validation cohort. The immunohistochemical results of the HPA database showed that AHSA1 and CEP55 had low expression in normal lung tissue but high expression in lung cancer tissue.

Conclusions

We constructed an IERG-based model for predicting the prognosis of LUAD. Among the genes identified, CEP55 and AHSA1 may be potential prognostic and therapeutic targets, and reducing their expression may represent a novel approach in the treatment of LUAD.

Prognostic and immunotherapeutic potential of regulatory T cell-associated signature in ovarian cancer

Tumour-induced immunosuppressive microenvironments facilitate oncogenesis, with regulatory T cells (Tregs) serving as a crucial component. The significance of Treg-associated genes within the context of ovarian cancer (OC) remains elucidated insufficiently. Utilizing single-cell RNA sequencing (scRNA-Seq) for the identification of Treg-specific biomarkers, this investigation employed single-sample gene set enrichment analysis (ssGSEA) for the derivation of a Treg signature score. Weighted gene co-expression network analysis (WGCNA) facilitated the identification of Treg-correlated genes. Machine learning algorithms were employed to determine an optimal prognostic model, subsequently exploring disparities across risk strata in terms of survival outcomes, immunological infiltration, pathway activation and responsiveness to immunotherapy. Through WGCNA, a cohort of 365 Treg-associated genes was discerned, with 70 implicated in the prognostication of OC. A Tregs-associated signature (TAS), synthesized from random survival forest (RSF) and Least Absolute Shrinkage and Selection Operator (LASSO) algorithms, exhibited robust predictive validity across both internal and external cohorts. Low TAS OC patients demonstrated superior survival outcomes, augmented by increased immunological cell infiltration, upregulated immune checkpoint expression, distinct pathway enrichment and differential response to immunotherapeutic interventions. The devised TAS proficiently prognosticates patient outcomes and delineates the immunological milieu within OC, offering a strategic instrument for the clinical stratification and selection of patients.

Integrated single-cell sequencing, spatial transcriptome sequencing and bulk RNA sequencing highlights the molecular characteristics of parthanatos in gastric cancer

BACKGROUND

Parthanatos is a novel programmatic form of cell death based on DNA damage and PARP-1 dependency. Nevertheless, its specific role in the context of gastric cancer (GC) remains uncertain.

METHODS

In this study, we integrated multi-omics algorithms to investigate the molecular characteristics of parthanatos in GC. A series of bioinformatics algorithms were utilized to explore clinical heterogeneity of GC and further predict the clinical outcomes.

RESULTS

Firstly, we conducted a comprehensive analysis of the omics features of parthanatos in various human tumors, including genomic mutations, transcriptome expression, and prognostic relevance. We successfully identified 7 cell types within the GC microenvironment: myeloid cell, epithelial cell, T cell, stromal cell, proliferative cell, B cell, and NK cell. When compared to adjacent non-tumor tissues, single-cell sequencing results from GC tissues revealed elevated scores for the parthanatos pathway across multiple cell types. Spatial transcriptomics, for the first time, unveiled the spatial distribution characteristics of parthanatos signaling. GC patients with different parthanatos signals often exhibited distinct immune microenvironment and metabolic reprogramming features, leading to different clinical outcomes. The integration of parthanatos signaling and clinical indicators enabled the creation of novel survival curves that accurately assess patients' survival times and statuses.

CONCLUSIONS

In this study, the molecular characteristics of parthanatos' unicellular and spatial transcriptomics in GC were revealed for the first time. Our model based on parthanatos signals can be used to distinguish individual heterogeneity and predict clinical outcomes in patients with GC.

Establishment of an immunogenic cell death-related model for prognostic prediction and identification of therapeutic targets in endometrial carcinoma

OBJECTIVE

Studies have firmly established the pivotal role of the immunogenic cell death (ICD) in the development of tumors. This study seeks to develop a risk model related to ICD to predict the prognosis of patients with endometrial carcinoma (EC).

MATERIALS AND METHODS

RNA-seq data of EC retrieved from TCGA database were analyzed using R software. We determined clusters based on ICD-related genes (ICDRGs) expression levels. Cox and LASSO analyses were further used to build the prediction model, and its accuracy was evaluated in the train and validation sets. Finally, in vitro and in vivo experiments were conducted to confirm the impact of the high-risk gene IFNA2 on EC.

RESULTS

Patients were sorted into two ICD clusters, with survival analysis revealing divergent prognoses between the clusters. The Cox regression analysis identified prognostic risk genes, and the LASSO analysis constructed a model based on 9 of these genes. Notably, this model displayed excellent predictive accuracy when validated. Finally, increased IFNA2 levels led to decreased vitality, proliferation, and invasiveness in vitro. IFNA2 also has significant tumor inhibiting effect in vivo.

CONCLUSIONS

The ICD-related model can accurately predict the prognosis of patients with EC, and IFNA2 may be a potential treatment target.

Single-cell transcriptomics reveals heterogeneity in esophageal squamous epithelial cells and constructs models for predicting patient prognosis and immunotherapy

Background

Esophageal squamous cell carcinoma (ESCC), characterized by its high invasiveness and malignant potential, has long been a formidable challenge in terms of treatment.

Methods

A variety of advanced analytical techniques are employed, including single-cell RNA sequencing (scRNA-seq), cell trajectory inference, transcription factor regulatory network analysis, GSVA enrichment analysis, mutation profile construction, and the inference of potential immunotherapeutic drugs. The purpose is to conduct a more comprehensive exploration of the heterogeneity among malignant squamous epithelial cell subgroups within the ESCC microenvironment and establish a model for predicting the prognosis and immunotherapy outcomes of ESCC patients.

Results

An analysis was conducted through scRNA-seq, and three Cluster of malignant epithelial cells were identified using the infer CNV method. Cluster 0 was found to exhibit high invasiveness, whereas Cluster 1 displayed prominent characteristics associated with epithelial-mesenchymal transition. Confirmation of these findings was provided through cell trajectory analysis, which positioned Cluster 0 at the initiation stage of development and Cluster 1 at the final developmental stage. The abundance of Cluster 0-2 groups in TCGA-LUAD samples was assessed using ssGSEA and subsequently categorized into high and low-expression groups. Notably, it was observed that Cluster 0-1 had a significant impact on survival (p<0.05). Furthermore, GSVA enrichment analysis demonstrated heightened activity in hallmark pathways for Cluster 0, whereas Cluster 1 exhibited notable enrichment in pathways related to cell proliferation. It is noteworthy that a prognostic model was established utilizing feature genes from Cluster 0-1, employing the Lasso and stepwise regression methods. The results revealed that in TCGA and GSE53624 cohorts, the low-risk group demonstrated significantly higher overall survival and increased levels of immune infiltration. An examination of four external immunotherapy cohorts unveiled that the low-risk group exhibited improved immunotherapeutic efficacy. Additionally, more meaningful treatment options were identified for the low-risk group.

Conclusion

The findings revealed distinct interactions between malignant epithelial cells of ESCC and subgroups within the tumor microenvironment. Two cell clusters, strongly linked to survival, were pinpointed, and a signature was formulated. This signature is expected to play a crucial role in identifying and advancing precision medicine approaches for the treatment of ESCC.

Single-cell analysis reveals exosome-associated biomarkers for prognostic prediction and immunotherapy in lung adenocarcinoma

BACKGROUND

Exosomes play a crucial role in tumor initiation and progression, yet the precise involvement of exosome-related genes (ERGs) in lung adenocarcinoma (LUAD) remains unclear.

METHODS

We conducted a comprehensive investigation of ERGs within the tumor microenvironment (TME) of LUAD using single-cell RNA sequencing (scRNA-seq) analysis. Multiple scoring methods were employed to assess exosome activity (EA). Differences in cell communication were examined between high and low EA groups, utilizing the "CellChat" R package. Subsequently, we leveraged multiple bulk RNA-seq datasets to develop and validate exosome-associated signatures (EAS), enabling a multifaceted exploration of prognosis and immunotherapy outcomes between high- and low-risk groups.

RESULTS

In the LUAD TME, epithelial cells demonstrated the highest EA, with even more elevated levels observed in advanced LUAD epithelial cells. The high-EA group exhibited enhanced intercellular interactions. EAS were established through the analysis of multiple bulk RNA-seq datasets. Patients in the high-risk group exhibited poorer overall survival (OS), reduced immune infiltration, and decreased expression of immune checkpoint genes. Finally, we experimentally validated the high expression of SEC61G in LUAD cell lines and demonstrated that knockdown of SEC61G reduced the proliferative capacity of LUAD cells using colony formation assays.

CONCLUSION

The integration of single-cell and bulk RNA-seq analyses culminated in the development of the profound and significant EAS, which imparts invaluable insights for the clinical diagnosis and therapeutic management of LUAD patients.

Unraveling the role of disulfidptosis-related LncRNAs in colon cancer: a prognostic indicator for immunotherapy response, chemotherapy sensitivity, and insights into cell death mechanisms

Background: Colon cancer, a prevalent and deadly malignancy worldwide, ranks as the third leading cause of cancer-related mortality. Disulfidptosis stress triggers a unique form of programmed cell death known as disulfidoptosis, characterized by excessive intracellular cystine accumulation. This study aimed to establish reliable bioindicators based on long non-coding RNAs (LncRNAs) associated with disulfidptosis-induced cell death, providing novel insights into immunotherapeutic response and prognostic assessment in patients with colon adenocarcinoma (COAD). Methods: Univariate Cox proportional hazard analysis and Lasso regression analysis were performed to identify differentially expressed genes strongly associated with prognosis. Subsequently, a multifactorial model for prognostic risk assessment was developed using multiple Cox proportional hazard regression. Furthermore, we conducted comprehensive evaluations of the characteristics of disulfidptosis response-related LncRNAs, considering clinicopathological features, tumor microenvironment, and chemotherapy sensitivity. The expression levels of prognosis-related genes in COAD patients were validated using quantitative real-time fluorescence PCR (qRT-PCR). Additionally, the role of ZEB1-SA1 in colon cancer was investigated through CCK8 assays, wound healing experiment and transwell experiments. Results: disulfidptosis response-related LncRNAs were identified as robust predictors of COAD prognosis. Multifactorial analysis revealed that the risk score derived from these LncRNAs served as an independent prognostic factor for COAD. Patients in the low-risk group exhibited superior overall survival (OS) compared to those in the high-risk group. Accordingly, our developed Nomogram prediction model, integrating clinical characteristics and risk scores, demonstrated excellent prognostic efficacy. In vitro experiments demonstrated that ZEB1-SA1 promoted the proliferation and migration of COAD cells. Conclusion: Leveraging medical big data and artificial intelligence, we constructed a prediction model for disulfidptosis response-related LncRNAs based on the TCGA-COAD cohort, enabling accurate prognostic prediction in colon cancer patients. The implementation of this model in clinical practice can facilitate precise classification of COAD patients, identification of specific subgroups more likely to respond favorably to immunotherapy and chemotherapy, and inform the development of personalized treatment strategies for COAD patients based on scientific evidence.

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 entanglement of extracellular matrix molecules and immune checkpoint inhibitors in cancer: a systematic review of the literature

Introduction

Immune-checkpoint inhibitors (ICIs) have emerged as a core pillar of cancer therapy as single agents or in combination regimens both in adults and children. Unfortunately, ICIs provide a long-lasting therapeutic effect in only one third of the patients. Thus, the search for predictive biomarkers of responsiveness to ICIs remains an urgent clinical need. The efficacy of ICIs treatments is strongly affected not only by the specific characteristics of cancer cells and the levels of immune checkpoint ligands, but also by other components of the tumor microenvironment, among which the extracellular matrix (ECM) is emerging as key player. With the aim to comprehensively describe the relation between ECM and ICIs' efficacy in cancer patients, the present review systematically evaluated the current literature regarding ECM remodeling in association with immunotherapeutic approaches.

Methods

This review followed the Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) guidelines and was registered at the International Prospective Register of Systematic Reviews (PROSPERO, CRD42022351180). PubMed, Web of Science, and Scopus databases were comprehensively searched from inception to January 2023. Titles, abstracts and full text screening was performed to exclude non eligible articles. The risk of bias was assessed using the QUADAS-2 tool.

Results

After employing relevant MeSH and key terms, we identified a total of 5070 studies. Among them, 2540 duplicates, 1521 reviews or commentaries were found and excluded. Following title and abstract screening, the full text was analyzed, and 47 studies meeting the eligibility criteria were retained. The studies included in this systematic review comprehensively recapitulate the latest observations associating changes of the ECM composition following remodeling with the traits of the tumor immune cell infiltration. The present study provides for the first time a broad view of the tight association between ECM molecules and ICIs efficacy in different tumor types, highlighting the importance of ECM-derived proteolytic products as promising liquid biopsy-based biomarkers to predict the efficacy of ICIs.

Conclusion

ECM remodeling has an important impact on the immune traits of different tumor types. Increasing evidence pinpoint at ECM-derived molecules as putative biomarkers to identify the patients that would most likely benefit from ICIs treatments.

Systematic review registration

https://www.crd.york.ac.uk/prospero/display_record.php?ID=CRD42022351180, identifier CRD42022351180.

Deciphering Treg cell roles in esophageal squamous cell carcinoma: a comprehensive prognostic and immunotherapeutic analysis

Background: Esophageal squamous cell carcinoma (ESCC) is a prevalent and aggressive form of cancer that poses significant challenges in terms of prognosis and treatment. Regulatory T cells (Treg cells) have gained attention due to their influential role in immune modulation within the tumor microenvironment (TME). Understanding the intricate interactions between Treg cells and the tumor microenvironment is essential for unraveling the mechanisms underlying ESCC progression and for developing effective prognostic models and immunotherapeutic strategies. Methods: A combination of single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq analysis was utilized to explore the role of Treg cells within the TME of ESCC. The accuracy and applicability of the prognostic model were assessed through multi-dimensional evaluations, encompassing an examination of the model's performance across various dimensions, such as the mutation landscape, clinical relevance, enrichment analysis, and potential implications for immunotherapy strategies. Results: The pivotal role of the macrophage migration inhibitory factor (MIF) signaling pathway within the ESCC TME was investigated, with a focus on its impact on Treg cells and other subpopulations. Through comprehensive integration of bulk sequencing data, a Treg-associated signature (TAS) was constructed, revealing that ESCC patients with elevated TAS (referred to as high-TAS individuals) experienced significantly improved prognoses. Heightened immune infiltration and increased expression of immune checkpoint markers were observed in high-TAS specimens. The model's validity was established through the IMvigor210 dataset, demonstrating its robustness in predicting prognosis and responsiveness to immunotherapy. Heightened therapeutic benefits were observed in immune-based interventions for high-TAS ESCC patients. Noteworthy differences in pathway enrichment patterns emerged between high and low-TAS cohorts, highlighting potential avenues for therapeutic exploration. Furthermore, the clinical relevance of key model genes was substantiated by analyzing clinical samples from ten paired tumor and adjacent tissues, revealing differential expression levels. Conclusion: The study established a TAS that enables accurate prediction of patient prognosis and responsiveness to immunotherapy. This achievement holds significant implications for the clinical management of ESCC, offering valuable insights for informed therapeutic interventions.

Uncovering the immune microenvironment and molecular subtypes of hepatitis B-related liver cirrhosis and developing stable a diagnostic differential model by machine learning and artificial neural networks

Background: Hepatitis B-related liver cirrhosis (HBV-LC) is a common clinical disease that evolves from chronic hepatitis B (CHB). The development of cirrhosis can be suppressed by pharmacological treatment. When CHB progresses to HBV-LC, the patient's quality of life decreases dramatically and drug therapy is ineffective. Liver transplantation is the most effective treatment, but the lack of donor required for transplantation, the high cost of the procedure and post-transplant rejection make this method unsuitable for most patients. Methods: The aim of this study was to find potential diagnostic biomarkers associated with HBV-LC by bioinformatics analysis and to classify HBV-LC into specific subtypes by consensus clustering. This will provide a new perspective for early diagnosis, clinical treatment and prevention of HCC in HBV-LC patients. Two study-relevant datasets, GSE114783 and GSE84044, were retrieved from the GEO database. We screened HBV-LC for feature genes using differential analysis, weighted gene co-expression network analysis (WGCNA), and three machine learning algorithms including least absolute shrinkage and selection operator (LASSO), support vector machine recursive feature elimination (SVM-RFE), and random forest (RF) for a total of five methods. After that, we constructed an artificial neural network (ANN) model. A cohort consisting of GSE123932, GSE121248 and GSE119322 was used for external validation. To better predict the risk of HBV-LC development, we also built a nomogram model. And multiple enrichment analyses of genes and samples were performed to understand the biological processes in which they were significantly enriched. And the different subtypes of HBV-LC were analyzed using the Immune infiltration approach. Results: Using the data downloaded from GEO, we developed an ANN model and nomogram based on six feature genes. And consensus clustering of HBV-LC classified them into two subtypes, C1 and C2, and it was hypothesized that patients with subtype C2 might have milder clinical symptoms by immune infiltration analysis. Conclusion: The ANN model and column line graphs constructed with six feature genes showed excellent predictive power, providing a new perspective for early diagnosis and possible treatment of HBV-LC. The delineation of HBV-LC subtypes will facilitate the development of future clinical treatment of HBV-LC.

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.

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 tumor immune evasion mechanisms: abnormal expression of transporters on immune cells in the tumor microenvironment

The tumor microenvironment (TME) is a crucial driving factor for tumor progression and it can hinder the body's immune response by altering the metabolic activity of immune cells. Both tumor and immune cells maintain their proliferative characteristics and physiological functions through transporter-mediated regulation of nutrient acquisition and metabolite efflux. Transporters also play an important role in modulating immune responses in the TME. In this review, we outline the metabolic characteristics of the TME and systematically elaborate on the effects of abundant metabolites on immune cell function and transporter expression. We also discuss the mechanism of tumor immune escape due to transporter dysfunction. Finally, we introduce some transporter-targeted antitumor therapeutic strategies, with the aim of providing new insights into the development of antitumor drugs and rational drug usage for clinical cancer therapy.

Identification of peripheral blood immune infiltration signatures and construction of monocyte-associated signatures in ovarian cancer and Alzheimer's disease using single-cell sequencing

Background

Ovarian cancer (OC) is a common tumor of the female reproductive system, while Alzheimer's disease (AD) is a prevalent neurodegenerative disease that primarily affects cognitive function in the elderly. Monocytes are immune cells in the blood that can enter tissues and transform into macrophages, thus participating in immune and inflammatory responses. Overall, monocytes may play an important role in Alzheimer's disease and ovarian cancer.

Methods

The CIBERSORT algorithm results indicate a potential crucial role of monocytes/macrophages in OC and AD. To identify monocyte marker genes, single-cell RNA-seq data of peripheral blood mononuclear cells (PBMCs) from OC and AD patients were analyzed. Enrichment analysis of various cell subpopulations was performed using the "irGSEA" R package. The estimation of cell cycle was conducted with the "tricycle" R package, and intercellular communication networks were analyzed using "CellChat". For 134 monocyte-associated genes (MRGs), bulk RNA-seq data from two diseased tissues were obtained. Cox regression analysis was employed to develop risk models, categorizing patients into high-risk (HR) and low-risk (LR) groups. The model's accuracy was validated using an external GEO cohort. The different risk groups were evaluated in terms of immune cell infiltration, mutational status, signaling pathways, immune checkpoint expression, and immunotherapy. To identify characteristic MRGs in AD, two machine learning algorithms, namely random forest and support vector machine (SVM), were utilized.

Results

Based on Cox regression analysis, a risk model consisting of seven genes was developed in OC, indicating a better prognosis for patients in the LR group. The LR group had a higher tumor mutation burden, immune cell infiltration abundance, and immune checkpoint expression. The results of the TIDE algorithm and the IMvigor210 cohort showed that the LR group was more likely to benefit from immunotherapy. Finally, ZFP36L1 and AP1S2 were identified as characteristic MRGs affecting OC and AD progression.

Conclusion

The risk profile containing seven genes identified in this study may help further guide clinical management and targeted therapy for OC. ZFP36L1 and AP1S2 may serve as biomarkers and new therapeutic targets for patients with OC and AD.

A fibroblast-associated signature predicts prognosis and immunotherapy in esophageal squamous cell cancer

Background

Current paradigms of anti-tumor therapies are not qualified to evacuate the malignancy ascribing to cancer stroma's functions in accelerating tumor relapse and therapeutic resistance. Cancer-associated fibroblasts (CAFs) has been identified significantly correlated with tumor progression and therapy resistance. Thus, we aimed to probe into the CAFs characteristics in esophageal squamous cancer (ESCC) and construct a risk signature based on CAFs to predict the prognosis of ESCC patients.

Methods

The GEO database provided the single-cell RNA sequencing (scRNA-seq) data. The GEO and TCGA databases were used to obtain bulk RNA-seq data and microarray data of ESCC, respectively. CAF clusters were identified from the scRNA-seq data using the Seurat R package. CAF-related prognostic genes were subsequently identified using univariate Cox regression analysis. A risk signature based on CAF-related prognostic genes was constructed using Lasso regression. Then, a nomogram model based on clinicopathological characteristics and the risk signature was developed. Consensus clustering was conducted to explore the heterogeneity of ESCC. Finally, PCR was utilized to validate the functions that hub genes play on ESCC.

Results

Six CAF clusters were identified in ESCC based on scRNA-seq data, three of which had prognostic associations. A total of 642 genes were found to be significantly correlated with CAF clusters from a pool of 17080 DEGs, and 9 genes were selected to generate a risk signature, which were mainly involved in 10 pathways such as NRF1, MYC, and TGF-Beta. The risk signature was significantly correlated with stromal and immune scores, as well as some immune cells. Multivariate analysis demonstrated that the risk signature was an independent prognostic factor for ESCC, and its potential in predicting immunotherapeutic outcomes was confirmed. A novel nomogram integrating the CAF-based risk signature and clinical stage was developed, which exhibited favorable predictability and reliability for ESCC prognosis prediction. The consensus clustering analysis further confirmed the heterogeneity of ESCC.

Conclusion

The prognosis of ESCC can be effectively predicted by CAF-based risk signatures, and a comprehensive characterization of the CAF signature of ESCC may aid in interpreting the response of ESCC to immunotherapy and offer new strategies for cancer treatment.

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.