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

Untargeted Pixel-by-Pixel Imaging of Metabolite Ratio Pairs as a Novel Tool for Biomedical Discovery in Mass Spectrometry Imaging

Mass spectrometry imaging (MSI) is a powerful technology used to define the spatial distribution and relative abundance of structurally identified and yet-undefined metabolites across tissue cryosections. While numerous software packages enable pixel-by-pixel imaging of individual metabolites, the research community lacks a discovery tool that images all metabolite abundance ratio pairs. Importantly, recognition of correlated metabolite pairs informs discovery of unanticipated molecules contributing to shared metabolic pathways, uncovers hidden metabolic heterogeneity across cells and tissue subregions, and indicates single-timepoint flux through pathways of interest. Here, we describe the development and implementation of an untargeted R package workflow for pixel-by-pixel ratio imaging of all metabolites detected in an MSI experiment. Considering untargeted MSI studies of murine brain and embryogenesis, we demonstrate that ratio imaging minimizes systematic data variation introduced by sample handling and instrument drift, markedly enhances spatial image resolution, and reveals previously unrecognized metabotype-distinct tissue regions. Furthermore, ratio imaging facilitates identification of novel regional biomarkers and provides anatomical information regarding spatial distribution of metabolite-linked biochemical pathways. The algorithm described herein is applicable to any MSI dataset containing spatial information for metabolites, peptides or proteins, offering a potent tool to enhance knowledge obtained from current spatial metabolite profiling technologies.

Unravelling the metabolic landscape of cutaneous melanoma: Insights from single-cell sequencing analysis and machine learning for prognostic assessment of lactate metabolism

This manuscript presents a comprehensive investigation into the role of lactate metabolism-related genes as potential prognostic markers in skin cutaneous melanoma (SKCM). Bulk-transcriptome data from The Cancer Genome Atlas (TCGA) and GSE19234, GSE22153, and GSE65904 cohorts from GEO database were processed and harmonized to mitigate batch effects. Lactate metabolism scores were assigned to individual cells using the 'AUCell' package. Weighted Co-expression Network Analysis (WGCNA) was employed to identify gene modules correlated with lactate metabolism. Machine learning algorithms were applied to construct a prognostic model, and its performance was evaluated in multiple cohorts. Immune correlation, mutation analysis, and enrichment analysis were conducted to further characterize the prognostic model's biological implications. Finally, the function of key gene NDUFS7 was verified by cell experiments. Machine learning resulted in an optimal prognostic model, demonstrating significant prognostic value across various cohorts. In the different cohorts, the high-risk group showed a poor prognosis. Immune analysis indicated differences in immune cell infiltration and checkpoint gene expression between risk groups. Mutation analysis identified genes with high mutation loads in SKCM. Enrichment analysis unveiled enriched pathways and biological processes in high-risk SKCM patients. NDUFS7 was found to be a hub gene in the protein-protein interaction network. After the expression of NDUFS7 was reduced by siRNA knockdown, CCK-8, colony formation, transwell and wound healing tests showed that the activity, proliferation and migration of A375 and WM115 cell lines were significantly decreased. This study offers insights into the prognostic significance of lactate metabolism-related genes in SKCM.

Integrative Single-Cell Analysis of Cardiomyopathy Identifies Differences in Cell Stemness and Transcriptional Regulatory Networks among Fibroblast Subpopulations

Background

Cardiomyopathy encompasses a broad spectrum of diseases affecting myocardial tissue, characterized clinically by abnormalities in cardiac structure, heart failure, and/or arrhythmias. Clinically heterogeneous, major types include dilated cardiomyopathy (DCM), hypertrophic cardiomyopathy (HCM), restrictive cardiomyopathy (RM), ischemic cardiomyopathy (ICM), among which DCM is more prevalent, while ICM exhibits higher incidence and mortality rates. Myocardial injury during cardiomyopathy progression may lead to myocardial fibrosis. Failure to intervene early and inhibit the process of myocardial fibrosis may culminate in heart failure. Cardiac fibroblasts constitute crucial cellular components determining the extent and quality of myocardial fibrosis, with various subpopulations exerting diverse roles in cardiomyopathy progression. Despite this, understanding of the cellular plasticity and transcriptional regulatory networks of cardiac fibroblasts in cardiomyopathy remains limited. Therefore, in this study, we conducted comprehensive single-cell analysis of cardiac fibroblasts in cardiomyopathy to explore differences in cellular plasticity and transcriptional regulatory networks among fibroblast subpopulations, with the aim of providing as many useful references as possible for the diagnosis, prognosis, and treatment of cardiomyopathy.

Materials and Methods

Cells with mitochondrial gene expression comprising >20% of total expressed genes were excluded. Differential expression genes (DEGs) and stemness genes within cardiac fibroblast subpopulations were subjected to Gene Ontology (GO) analysis of biological processes (BP) and AUCell analysis. Monocle software was employed to analyze the pseudo-temporal trajectory of cardiac fibroblasts in cardiomyopathy. Additionally, the Python package SCENIC was utilized to assess enrichment of transcription factors and activity of regulators within cardiac fibroblast subpopulations in cardiomyopathy.

Results

Following batch effect correction, 179,927 cells were clustered into 32 clusters, designated as T_NK cells, endothelial cells, myeloid cells, fibroblasts, pericytes, SMCs, CMs, proliferating cells, EndoCs, and EPCs. Among them, 8148 fibroblasts were further subdivided into 4 subpopulations, namely C0 THBS4+ Fibroblasts, C1 LINC01133+ Fibroblasts, C2 FGF7+ Fibroblasts, and C3 AGT + Fibroblasts. Results from GO_BP and AUCell analyses suggest that C3 AGT + Fibroblasts may be associated with immune response activation, protein transport, and myocardial contractile function, correlating with disease progression in cardiomyopathy. Transcription factor enrichment analysis indicates that FOS is the most significant TF in C3 AGT + Fibroblasts, also associated with the M1 module, possibly implicated in protein hydrolysis, intracellular DNA replication, and cell proliferation. Moreover, correlation analysis of transcriptional regulatory activity between fibroblast subpopulations reveals a more pronounced heterogeneity within C3 AGT + Fibroblasts in cardiomyopathy.

Conclusion

C3 AGT + Fibroblasts exhibit increased sensitivity towards adverse outcomes in cardiomyopathy, such as myocardial fibrosis and impaired cardiac contractile function, compared to other cardiac fibroblast subpopulations. The differential cellular plasticity and transcriptional regulatory activity between C3 AGT + Fibroblasts and other subgroups offer new perspectives for targeting fibroblast subpopulation activity to treat cardiomyopathy. Additionally, stemness genes EPAS1 and MYC, along with the regulator FOS, may play roles in modulating the biological processes of cardiac fibroblasts in cardiomyopathy.

Identification of necroptosis-related gene signatures for predicting the prognosis of ovarian cancer

Ovarian cancer (OC) is one of the most prevalent and fatal malignant tumors of the female reproductive system. Our research aimed to develop a prognostic model to assist inclinical treatment decision-making.Utilizing data from The Cancer Genome Atlas (TCGA) and copy number variation (CNV) data from the University of California Santa Cruz (UCSC) database, we conducted analyses of differentially expressed genes (DEGs), gene function, and tumor microenvironment (TME) scores in various clusters of OC samples.Next, we classified participants into low-risk and high-risk groups based on the median risk score, thereby dividing both the training group and the entire group accordingly. Overall survival (OS) was significantly reduced in the high-risk group, and two independent prognostic factors were identified: age and risk score. Additionally, three genes-C-X-C Motif Chemokine Ligand 10 (CXCL10), RELB, and Caspase-3 (CASP3)-emerged as potential candidates for an independent prognostic signature with acceptable prognostic value. In Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, pathways related to immune responses and inflammatory cell chemotaxis were identified. Cellular experiments further validated the reliability and precision of our findings. In conclusion, necroptosis-related genes play critical roles in tumor immunity, and our model introduces a novel strategy for predicting the prognosis of OC patients.

A systematic pan-cancer analysis identifies LDHA as a novel predictor for immunological, prognostic, and immunotherapy resistance

Lactate dehydrogenase A (LDHA), a critical enzyme involved in glycolysis, is broadly involved multiple biological functions in human cancers. It is reported that LDHA can impact tumor immune surveillance and induce the transformation of tumor-associated macrophages, highlighting its unnoticed function of LDHA in immune system. However, in human cancers, the role of LDHA in prognosis and immunotherapy hasn't been investigated. In this study, we analyzed the expression pattern and prognostic value of LDHA in pan-cancer and explored its association between tumor microenvironment (TME), immune infiltration subtype, stemness scores, tumor mutation burden (TMB), and immunotherapy resistance. We found that LDHA expression is tumor heterogeneous and that its high expression is associated with poor prognosis in multiple human cancers. In addition, LDHA expression was positively correlated with the presence of mononuclear/macrophage cells, and also promoted the infiltration of a range of immune cells. Genomic alteration of LDHA was common in different types of cancer, while with prognostic value in pan-cancers. Pan-cancer analysis revealed that the significant correlations existed between LDHA expression and tumor microenvironment (including stromal cells and immune cells) as well as stemness scores (DNAss and RNAss) across cancer types. Drug sensitivity analysis also revealed that LDHA was able to predict response to chemotherapy and immunotherapy. Furthermore, it was confirmed that knockdown of LDHA reduced proliferation and migration ability of lung cancer cells. Taken together, LDHA could serve as a prognostic biomarker and a potential immunotherapy marker.

Targeting nucleotide metabolic pathways in colorectal cancer by integrating scRNA-seq, spatial transcriptome, and bulk RNA-seq data

BACKGROUND

Colorectal cancer is a malignant tumor of the digestive system originating from abnormal cell proliferation in the colon or rectum, often leading to gastrointestinal symptoms and severe health issues. Nucleotide metabolism, which encompasses the synthesis of DNA and RNA, is a pivotal cellular biochemical process that significantly impacts both the progression and therapeutic strategies of colorectal cancer METHODS: For single-cell RNA sequencing (scRNA-seq), five functions were employed to calculate scores related to nucleotide metabolism. Cell developmental trajectory analysis and intercellular interaction analysis were utilized to explore the metabolic characteristics and communication patterns of different epithelial cells. These findings were further validated using spatial transcriptome RNA sequencing (stRNA-seq). A risk model was constructed using expression profile data from TCGA and GEO cohorts to optimize clinical decision-making. Key nucleotide metabolism-related genes (NMRGs) were functionally validated by further in vitro experiments.

RESULTS

In both scRNA-seq and stRNA-seq, colorectal cancer (CRC) exhibited unique cellular heterogeneity, with myeloid cells and epithelial cells in tumor samples displaying higher nucleotide metabolism scores. Analysis of intercellular communication revealed enhanced signaling pathways and ligand-receptor interactions between epithelial cells with high nucleotide metabolism and fibroblasts. Spatial transcriptome sequencing confirmed elevated nucleotide metabolism states in the core region of tumor tissue. After identifying differentially expressed NMRGs in epithelial cells, a risk prognostic model based on four genes effectively predicted overall survival and immunotherapy outcomes in patients. High-risk group patients exhibited an immunosuppressive microenvironment and relatively poorer prognosis and responses to chemotherapy and immunotherapy. Finally, based on data analysis and a series of cellular functional experiments, ACOX1 and CPT2 were identified as novel therapeutic targets for CRC.

CONCLUSION

In this study, a comprehensive analysis of NMRGs in CRC was conducted using a combination of single-cell sequencing, spatial transcriptome sequencing, and high-throughput data. The prognostic model constructed with NMRGs shows potential as a standalone prognostic marker for colorectal cancer patients and may significantly influence the development of personalized treatment approaches for CRC.

Deciphering cutaneous melanoma prognosis through LDL metabolism: Single-cell transcriptomics analysis via 101 machine learning algorithms

Cutaneous melanoma poses a formidable challenge within the field of oncology, marked by its aggressive nature and capacity for metastasis. Despite extensive research uncovering numerous genetic and molecular contributors to cutaneous melanoma development, there remains a critical knowledge gap concerning the role of lipids, notably low-density lipoprotein (LDL), in this lethal skin cancer. This article endeavours to bridge this knowledge gap by delving into the intricate interplay between LDL metabolism and cutaneous melanoma, shedding light on how lipids influence tumour progression, immune responses and potential therapeutic avenues. Genes associated with LDL metabolism were extracted from the GSEA database. We acquired and analysed single-cell sequencing data (GSE215120) and bulk-RNA sequencing data, including the TCGA data set, GSE19234, GSE22153 and GSE65904. Our analysis unveiled the heterogeneity of LDL across various cell types at the single-cell sequencing level. Additionally, we constructed an LDL-related signature (LRS) using machine learning algorithms, incorporating differentially expressed genes and highly correlated genes. The LRS serves as a valuable tool for assessing the prognosis, immunity and mutation status of patients with cutaneous melanoma. Furthermore, we conducted experiments on A375 and WM-115 cells to validate the function of PPP2R1A, a pivotal gene within the LRS. Our comprehensive approach, combining advanced bioinformatics analyses with an extensive review of current literature, presents compelling evidence regarding the significance of LDL within the cutaneous melanoma microenvironment.

Combining bulk and scRNA-seq to explore the molecular mechanisms governing the distinct efferocytosis activities of a macrophage subpopulation in PDAC

Pancreatic ductal adenocarcinoma (PDAC), a very aggressive tumour, is currently the third leading cause of cancer-related deaths. Unfortunately, many patients face the issue of inoperability at the diagnostic phase leading to a quite dismal prognosis. The onset of metastatic processes has a crucial role in the elevated mortality rates linked to PDAC. Individuals with metastatic advances receive only palliative therapy and have a grim prognosis. It is essential to carefully analyse the intricacies of the metastatic process to enhance the prognosis for individuals with PDAC. Malignancy development is greatly impacted by the process of macrophage efferocytosis. Our current knowledge about the complete range of macrophage efferocytosis activities in PDAC and their intricate interactions with tumour cells is still restricted. This work aims to resolve communication gaps and pinpoint the essential transcription factor that is vital in the immunological response of macrophage populations. We analysed eight PDAC tissue samples sourced from the gene expression omnibus. We utilized several software packages such as Seurat, DoubletFinder, Harmony, Pi, GSVA, CellChat and Monocle from R software together with pySCENIC from Python, to analyse the single-cell RNA sequencing (scRNA-seq) data collected from the PDAC samples. This study involved the analysis of a comprehensive sample of 22,124 cells, which were classified into distinct cell types. These cell types encompassed endothelial and epithelial cells, PDAC cells, as well as various immune cells, including CD4+ T cells, CD8+ T cells, NK cells, B cells, plasma cells, mast cells, monocytes, DC cells and different subtypes of macrophages, namely C0 macrophage TGM2+, C1 macrophage PFN1+, C2 macrophage GAS6+ and C3 macrophage APOC3+. The differentiation between tumour cells and epithelial cells was achieved by the implementation of CopyKat analysis, resulting in the detection and categorization of 1941 PDAC cells. The amplification/deletion patterns observed in PDAC cells on many chromosomes differ significantly from those observed in epithelial cells. The study of Pseudotime Trajectories demonstrated that the C0 macrophage subtype expressing TGM2+ had the lowest level of differentiation. Additionally, the examination of gene set scores related to efferocytosis suggested that this subtype displayed higher activity during the efferocytosis process compared to other subtypes. The most active transcription factors for each macrophage subtype were identified as BACH1, NFE2, TEAD4 and ARID3A. In conclusion, the examination of human PDAC tissue samples using immunofluorescence analysis demonstrated the co-localization of CD68 and CD11b within regions exhibiting the presence of keratin (KRT) and alpha-smooth muscle actin (α-SMA). This observation implies a spatial association between macrophages, fibroblasts, and epithelial cells. There is variation in the expression of efferocytosis-associated genes between C0 macrophage TGM2+ and other macrophage cell types. This observation implies that the diversity of macrophage cells might potentially influence the metastatic advancement of PDAC. Moreover, the central transcription factor of different macrophage subtypes offers a promising opportunity for targeted immunotherapy in the treatment of PDAC.

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.

Deciphering the molecular landscape: integrating single-cell transcriptomics to unravel myofibroblast dynamics and therapeutic targets in clear cell renal cell carcinomas

Background

Clear cell renal cell carcinomas (ccRCCs) epitomize the most formidable clinical subtype among renal neoplasms. While the impact of tumor-associated fibroblasts on ccRCC progression is duly acknowledged, a paucity of literature exists elucidating the intricate mechanisms and signaling pathways operative at the individual cellular level.

Methods

Employing single-cell transcriptomic analysis, we meticulously curated UMAP profiles spanning substantial ccRCC populations, delving into the composition and intrinsic signaling pathways of these cohorts. Additionally, Myofibroblasts were fastidiously categorized into discrete subpopulations, with a thorough elucidation of the temporal trajectory relationships between these subpopulations. We further probed the cellular interaction pathways connecting pivotal subpopulations with tumors. Our endeavor also encompassed the identification of prognostic genes associated with these subpopulations through Bulk RNA-seq, subsequently validated through empirical experimentation.

Results

A notable escalation in the nFeature and nCount of Myofibroblasts and EPCs within ccRCCs was observed, notably enriched in oxidation-related pathways. This phenomenon is postulated to be closely associated with the heightened metabolic activities of Myofibroblasts and EPCs. The Myofibroblasts subpopulation, denoted as C3 HMGA1+ Myofibroblasts, emerges as a pivotal subset, displaying low differentiation and positioning itself at the terminal point of the temporal trajectory. Intriguingly, these cells exhibit a high degree of interaction with tumor cells through the MPZ signaling pathway network, suggesting that Myofibroblasts may facilitate tumor progression via this pathway. Prognostic genes associated with C3 were identified, among which TUBB3 is implicated in potential resistance to tumor recurrence. Finally, experimental validation revealed that the knockout of the key gene within the MPZ pathway, MPZL1, can inhibit tumor activity, proliferation, invasion, and migration capabilities.

Conclusion

This investigation delves into the intricate mechanisms and interaction pathways between Myofibroblasts and ccRCCs at the single-cell level. We propose that targeting MPZL1 and the oxidative phosphorylation pathway could serve as potential key targets for treating the progression and recurrence of ccRCC. This discovery paves the way for new directions in the treatment and prognosis diagnosis of ccRCC in the future.

Investigating the mechanisms of drug resistance and prognosis in ovarian cancer using single-cell RNA sequencing and bulk RNA sequencing

Ovarian cancer stands as a prevalent malignancy within the realm of gynecology, and the emergence of resistance to chemotherapeutic agents remains a pivotal impediment to both prognosis and treatment. Through a single-cell level investigation, we scrutinize the drug resistance and mitotic activity of the core tumor cells in ovarian cancer. Our study revisits the interrelationships and temporal trajectories of distinct epithelial cells (EPCs) subpopulations, while identifying genes associated with ovarian cancer prognosis. Notably, our findings establish a strong association between the drug resistance of EPCs and oxidative phosphorylation pathways. Subsequently, through subpopulation and temporal trajectory analysis, we confirm the intermediate position of EPCs subpopulation C0. Furthermore, we delve into the immunological functions and differentially expressed genes associated with the prognosis of C0, shedding light on the potential for constructing novel ovarian cancer prognosis models and identifying new therapeutic targets.

Unraveling the intricacies of glioblastoma progression and recurrence: insights into the role of NFYB and oxidative phosphorylation at the single-cell level

Background

Glioblastoma (GBM), with its high recurrence and mortality rates, makes it the deadliest neurological malignancy. Oxidative phosphorylation is a highly active cellular pathway in GBM, and NFYB is a tumor-associated transcription factor. Both are related to mitochondrial function, but studies on their relationship with GBM at the single-cell level are still scarce.

Methods

We re-analyzed the single-cell profiles of GBM from patients with different subtypes by single-cell transcriptomic analysis and further subdivided the large population of Glioma cells into different subpopulations, explored the interrelationships and active pathways among cell stages and clinical subtypes of the populations, and investigated the relationship between the transcription factor NFYB of the key subpopulations and GBM, searching for the prognostic genes of GBM related to NFYB, and verified by experiments.

Results

Glioma cells and their C5 subpopulation had the highest percentage of G2M staging and rGBM, which we hypothesized might be related to the higher dividing and proliferating ability of both Glioma and C5 subpopulations. Oxidative phosphorylation pathway activity is elevated in both the Glioma and C5 subgroup, and NFYB is a key transcription factor for the C5 subgroup, suggesting its possible involvement in GBM proliferation and recurrence, and its close association with mitochondrial function. We also identified 13 prognostic genes associated with NFYB, of which MEM60 may cause GBM patients to have a poor prognosis by promoting GBM proliferation and drug resistance. Knockdown of the NFYB was found to contribute to the inhibition of proliferation, invasion, and migration of GBM cells.

Conclusion

These findings help to elucidate the key mechanisms of mitochondrial function in GBM progression and recurrence, and to establish a new prognostic model and therapeutic target based on NFYB.

Decoding the tumor microenvironment and molecular mechanism: unraveling cervical cancer subpopulations and prognostic signatures through scRNA-Seq and bulk RNA-seq analyses

Background

Cervical carcinoma (CC) represents a prevalent gynecological neoplasm, with a discernible rise in prevalence among younger cohorts observed in recent years. Nonetheless, the intrinsic cellular heterogeneity of CC remains inadequately investigated.

Methods

We utilized single-cell RNA sequencing (scRNA-seq) transcriptomic analysis to scrutinize the tumor epithelial cells derived from four specimens of cervical carcinoma (CC) patients. This method enabled the identification of pivotal subpopulations of tumor epithelial cells and elucidation of their contributions to CC progression. Subsequently, we assessed the influence of associated molecules in bulk RNA sequencing (Bulk RNA-seq) cohorts and performed cellular experiments for validation purposes.

Results

Through our analysis, we have discerned C3 PLP2+ Tumor Epithelial Progenitor Cells as a noteworthy subpopulation in cervical carcinoma (CC), exerting a pivotal influence on the differentiation and progression of CC. We have established an independent prognostic indicator-the PLP2+ Tumor EPCs score. By stratifying patients into high and low score groups based on the median score, we have observed that the high-score group exhibits diminished survival rates compared to the low-score group. The correlations observed between these groups and immune infiltration, enriched pathways, single-nucleotide polymorphisms (SNPs), drug sensitivity, among other factors, further underscore their impact on CC prognosis. Cellular experiments have validated the significant impact of ATF6 on the proliferation and migration of CC cell lines.

Conclusion

This study enriches our comprehension of the determinants shaping the progression of CC, elevates cognizance of the tumor microenvironment in CC, and offers valuable insights for prospective CC therapies. These discoveries contribute to the refinement of CC diagnostics and the formulation of optimal therapeutic approaches.

Single cell sequencing revealed the mechanism of CRYAB in glioma and its diagnostic and prognostic value

Background

We explored the characteristics of single-cell differentiation data in glioblastoma and established prognostic markers based on CRYAB to predict the prognosis of glioblastoma patients. Aberrant expression of CRYAB is associated with invasive behavior in various tumors, including glioblastoma. However, the specific role and mechanisms of CRYAB in glioblastoma are still unclear.

Methods

We assessed RNA-seq and microarray data from TCGA and GEO databases, combined with scRNA-seq data on glioma patients from GEO. Utilizing the Seurat R package, we identified distinct survival-related gene clusters in the scRNA-seq data. Prognostic pivotal genes were discovered through single-factor Cox analysis, and a prognostic model was established using LASSO and stepwise regression algorithms. Moreover, we investigated the predictive potential of these genes in the immune microenvironment and their applicability in immunotherapy. Finally, in vitro experiments confirmed the functional significance of the high-risk gene CRYAB.

Results

By analyzing the ScRNA-seq data, we identified 28 cell clusters representing seven cell types. After dimensionality reduction and clustering analysis, we obtained four subpopulations within the oligodendrocyte lineage based on their differentiation trajectory. Using CRYAB as a marker gene for the terminal-stage subpopulation, we found that its expression was associated with poor prognosis. In vitro experiments demonstrated that knocking out CRYAB in U87 and LN229 cells reduced cell viability, proliferation, and invasiveness.

Conclusion

The risk model based on CRYAB holds promise in accurately predicting glioblastoma. A comprehensive study of the specific mechanisms of CRYAB in glioblastoma would contribute to understanding its response to immunotherapy. Targeting the CRYAB gene may be beneficial for glioblastoma patients.

Interplay of Sphingolipid Metabolism in Predicting Prognosis of GBM Patients: Towards Precision Immunotherapy

Background: In spite of numerous existing bio-surveillance systems for predicting glioma (GBM) prognosis, enhancing the efficacy of immunotherapy remains an ongoing conundrum. The continual scrutiny of the dynamic interplay between the sphingolipid metabolic pathway and tumor immunophenotypes has unveiled potential implications. However, the intricate orchestration of functional and regulatory mechanisms by long non-coding RNAs (lncRNAs) in GBM, particularly in the context of sphingolipid metabolism, remains cryptic. Methods: We harnessed established R packages to intersect gene expression profiles of GBM patients within the The Cancer Genome Atlas (TCGA) database with the compilation of sphingolipid metabolism genes from GeneCards. This enabled us to discern markedly distinct lncRNAs, which were subsequently deployed to construct a robust prognostic model utilizing Lasso-Cox regression analysis. We then scrutinized the immune microenvironment across various risk strata using the ssGSEA and CIBERSORT algorithms. To evaluate mutation patterns and drug resistance profiles within patient subgroups, we devised the "Prophytic" and "Maftools" packages, respectively. Results: Our investigation scrutinized lncRNAs linked to sphingolipid metabolism, utilizing glioma specimens from TCGA. We meticulously curated 1224 sphingolipid-associated genes gleaned from GeneCards and pinpointed 272 differentially expressed mRNAs via transcriptomic analysis. Enrichment analyses underscored their significance in sphingolipid processes. A prognostic model founded on 17 meticulously selected lncRNAs was systematically constructed and validated. This model adeptly stratified GBM patients into high- and low-risk categories, yielding highly precise prognostic insights. We also discerned correlations between immune cell infiltration and genetic mutation discrepancies, along with distinct therapeutic responses through drug sensitivity analysis. Notably, computational findings were corroborated through experimental validation by RT-PCR. Conclusion: In summation, our exhaustive inquiry underscores the multifaceted utility of the sphingolipid metabolic pathway as an autonomous diagnostic and prognostic indicator for glioma patients. Furthermore, we amalgamate a profusion of substantiated evidence concerning immune infiltration and gene mutations, thereby reinforcing the proposition that sphingolipid metabolism may function as a pivotal determinant in the panorama of immunotherapeutic interventions.

Exploring the immune microenvironment of osteosarcoma through T cell exhaustion-associated gene expression: a study on prognosis prediction

Background

Osteosarcoma is a highly aggressive type of bone cancer with a poor prognosis. In the tumor immune microenvironment, T-cell exhaustion can occur due to various factors, leading to reduced tumor-killing ability. The purpose of this study was to construct a prognostic model based on T-cell exhaustion-associated genes in osteosarcoma.

Methods

Patient data for osteosarcoma were retrieved from the TARGET and GEO databases. Consensus clustering was employed to identify two novel molecular subgroups. The dissimilarities in the tumor immune microenvironment between these subgroups were evaluated using the "xCell" algorithm. GO and KEGG analyses were conducted to elucidate the underlying mechanisms of gene expression. Predictive risk models were constructed using the least absolute shrinkage and selection operator algorithm and Cox regression analysis. To validate the prognostic significance of the risk gene expression model at the protein level, immunohistochemistry assays were performed on osteosarcoma patient samples. Subsequently, functional analysis of the key risk gene was carried out through in vitro experimentation.

Results

Four gene expression signatures (PLEKHO2, GBP2, MPP1, and VSIG4) linked to osteosarcoma prognosis were identified within the TARGET-osteosarcoma cohort, categorizing patients into two subgroups. The resulting prognostic model showed strong predictive capability, with area under the receiver operating characteristic curve (AUC) values of 0.728/0.740, 0.781/0.658, and 0.788/0.642 for 1, 3, and 5-year survival in both training and validation datasets. Notably, patients in the low-risk group had significantly higher stromal, immune, and ESTIMATE scores compared to high-risk counterparts. Additionally, a nomogram was developed, exhibiting high accuracy in predicting the survival outcome of osteosarcoma patients. Immunohistochemistry, Kaplan-Meier, and time-dependent AUC analyses consistently supported the prognostic value of the risk model within our osteosarcoma patient cohort. In vitro experiments provided additional validation by demonstrating that the downregulation of GBP2 promoted the proliferation, migration, and invasion of osteosarcoma cells while inhibiting apoptosis.

Conclusion

The current study established a prognostic signature associated with TEX-related genes and elucidated the impact of the pivotal gene GBP2 on osteosarcoma cells via in vitro experiments. Consequently, it introduces a fresh outlook for clinical prognosis prediction and sets the groundwork for targeted therapy investigations in osteosarcoma.

Combined signature of G protein-coupled receptors and tumor microenvironment provides a prognostic and therapeutic biomarker for skin cutaneous melanoma

BACKGROUND

G protein-coupled receptors (GPCRs) have been shown to have an important role in tumor development and metastasis, and abnormal expression of GPCRs is significantly associated with poor prognosis of tumor patients. In this study, we analyzed the GPCRs-related gene (GPRGs) and tumor microenvironment (TME) in skin cutaneous melanoma (SKCM) to construct a prognostic model to help SKCM patients obtain accurate clinical treatment strategies.

METHODS

SKCM expression data and clinical information were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Differential expression analysis, LASSO algorithm, and univariate and multivariate cox regression analysis were used to screen prognosis-related genes (GPR19, GPR146, S1PR2, PTH1R, ADGRE5, CXCR3, GPR143, and OR2I1P) and multiple prognosis-good immune cells; the data set was analyzed according to above results and build up a GPR-TME classifier. The model was further subjected to immune infiltration, functional enrichment, tumor mutational load, immunotherapy prediction, and scRNA-seq data analysis. Finally, cellular experiments were conducted to validate the functionality of the key gene GPR19 in the model.

RESULTS

The findings indicate that high expression of GPRGs is associated with a poor prognosis in patients with SKCM, highlighting the significant role of GPRGs and the tumor microenvironment (TME) in SKCM development. Notably, the group characterized by low GPR expression and a high TME exhibited the most favorable prognosis and immunotherapeutic efficacy. Furthermore, cellular assays demonstrated that knockdown of GPR19 significantly reduced the proliferation, migration, and invasive capabilities of melanoma cells in A375 and A2058 cell lines.

CONCLUSION

This study provides novel insights for the prognosis evaluation and treatment of melanoma, along with the identification of a new biomarker, GPR19.

Integrative lactylation and tumor microenvironment signature as prognostic and therapeutic biomarkers in skin cutaneous melanoma

BACKGROUND

The incidence of skin cutaneous melanoma (SKCM), one of the most aggressive and lethal skin tumors, is increasing worldwide. However, for advanced SKCM, we still lack an accurate and valid way to predict its prognosis, as well as novel theories to guide the planning of treatment options for SKCM patients. Lactylation (LAC), a novel post-translational modification of histones, has been shown to promote tumor growth and inhibit the antitumor response of the tumor microenvironment (TME) in a variety of ways. We hope that this study will provide new ideas for treatment options for SKCM patients, as well as research on the molecular mechanisms of SKCM pathogenesis and development.

METHODS

At the level of the RNA sequencing set (TCGA, GTEx), we used differential expression analysis, LASSO regression analysis, and multifactor Cox regression analysis to screen for prognosis-related genes and calculate the corresponding LAC scores. The content of TME cells in the tumor tissue was calculated using the CIBERSORT algorithm, and the TME score was calculated based on its results. Finally, the LAC-TME classifier was established and further analyzed based on the two scores, including the construction of a prognostic model, analysis of clinicopathological characteristics, and correlation analysis of tumor mutation burden (TMB) and immunotherapy. Based on single-cell RNA sequencing data, this study analyzed the cellular composition in SKCM tissues and explored the role of LAC scores in intercellular communication. To validate the functionality of the pivotal gene CLPB in the model, cellular experiments were ultimately executed.

RESULTS

We screened a total of six prognosis-related genes (NDUFA10, NDUFA13, CLPB, RRM2B, HPDL, NARS2) and 7 TME cells with good prognosis. According to Kaplan-Meier survival analysis, we found that the LAClow/TMEhigh group had the highest overall survival (OS) and the LAChigh/TMElow group had the lowest OS (p value < 0.05). In further analysis of immune infiltration, tumor microenvironment (TME), functional enrichment, tumor mutational load and immunotherapy, we found that immunotherapy was more appropriate in the LAClow/TMEhigh group. Moreover, the cellular assays exhibited substantial reductions in proliferation, migration, and invasive potentials of melanoma cells in both A375 and A2058 cell lines upon CLPB knockdown.

CONCLUSIONS

The prognostic model using the combined LAC score and TME score was able to predict the prognosis of SKCM patients more consistently, and the LAC-TME classifier was able to significantly differentiate the prognosis of SKCM patients across multiple clinicopathological features. The LAC-TME classifier has an important role in the development of immunotherapy regimens for SKCM patients.

Single-cell sequencing analysis related to sphingolipid metabolism guides immunotherapy and prognosis of skin cutaneous melanoma

Background

We explore sphingolipid-related genes (SRGs) in skin melanoma (SKCM) to develop a prognostic indicator for patient outcomes. Dysregulated lipid metabolism is linked to aggressive behavior in various cancers, including SKCM. However, the exact role and mechanism of sphingolipid metabolism in melanoma remain partially understood.

Methods

We integrated scRNA-seq data from melanoma patients sourced from the GEO database. Through the utilization of the Seurat R package, we successfully identified distinct gene clusters associated with patient survival in the scRNA-seq data. Key prognostic genes were identified through single-factor Cox analysis and used to develop a prognostic model using LASSO and stepwise regression algorithms. Additionally, we evaluated the predictive potential of these genes within the immune microenvironment and their relevance to immunotherapy. Finally, we validated the functional significance of the high-risk gene IRX3 through in vitro experiments.

Results

Analysis of scRNA-seq data identified distinct expression patterns of 4 specific genes (SRGs) in diverse cell subpopulations. Re-clustering cells based on increased SRG expression revealed 7 subgroups with significant prognostic implications. Using marker genes, lasso, and Cox regression, we selected 11 genes to construct a risk signature. This signature demonstrated a strong correlation with immune cell infiltration and stromal scores, highlighting its relevance in the tumor microenvironment. Functional studies involving IRX3 knockdown in A375 and WM-115 cells showed significant reductions in cell viability, proliferation, and invasiveness.

Conclusion

SRG-based risk signature holds promise for precise melanoma prognosis. An in-depth exploration of SRG characteristics offers insights into immunotherapy response. Therapeutic targeting of the IRX3 gene may benefit melanoma patients.

Machine learning and single-cell transcriptome profiling reveal regulation of fibroblast activation through THBS2/TGFβ1/P-Smad2/3 signalling pathway in hypertrophic scar

Hypertrophic scar (HS) is a chronic inflammatory skin disorder characterized by excessive deposition of extracellular matrix, and the mechanisms underlying their formation remain poorly understood. We analysed scRNA-seq data from samples of normal skin and HS. Using the hdWGCNA method, key gene modules of fibroblasts in HS were identified. Non-negative matrix factorization was employed to perform subtype analysis of HS patients using these gene modules. Multiple machine learning algorithms were applied to screen and validate accurate gene signatures for identifying and predicting HS, and a convolutional neural network (CNN) based on deep learning was established and validated. Quantitative reverse transcription-polymerase chain reaction and western blotting were performed to measure mRNA and protein expression. Immunofluorescence was used for gene localization analysis, and biological features were assessed through CCK8 and wound healing assay. Single-cell sequencing revealed distinct subpopulations of fibroblasts in HS. HdWGCNA identified key gene characteristics of this population, and pseudotime analysis was conducted to investigate gene variation during fibroblast differentiation. By employing various machine learning algorithms, the gene range was narrowed down to three key genes. A CNN was trained using the expression of these key genes and immune cell infiltration, enabling diagnosis and prediction of HS. Functional experiments demonstrated that THBS2 is associated with fibroblast proliferation and migration in HS and affects the formation and development of HS through the TGFβ1/P-Smad2/3 pathway. Our study identifies unique fibroblast subpopulations closely associated with HS and provides biomarkers for the diagnosis and treatment of HS.

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.

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.

Unraveling molecular networks in thymic epithelial tumors: deciphering the unique signatures

Thymic epithelial tumors (TETs) are a rare and diverse group of neoplasms characterized by distinct molecular signatures. This review delves into the complex molecular networks of TETs, highlighting key aspects such as chromosomal abnormalities, molecular subtypes, aberrant gene mutations and expressions, structural gene rearrangements, and epigenetic changes. Additionally, the influence of the dynamic tumor microenvironment on TET behavior and therapeutic responses is examined. A thorough understanding of these facets elucidates TET pathogenesis, offering avenues for enhancing diagnostic accuracy, refining prognostic assessments, and tailoring targeted therapeutic strategies. Our review underscores the importance of deciphering TETs' unique molecular signatures to advance personalized treatment paradigms and improve patient outcomes. We also discuss future research directions and anticipated challenges in this intriguing field.

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.

Combining single-cell sequencing and spatial transcriptome sequencing to identify exosome-related features of glioblastoma and constructing a prognostic model to identify BARD1 as a potential therapeutic target for GBM patients

Background

Glioblastoma (GBM) is a malignant primary brain tumor. This study focused on exploring the exosome-related features of glioblastoma to better understand its cellular composition and molecular characteristics.

Methods

Single-cell RNA sequencing (scRNA-seq) and spatial transcriptome RNA sequencing (stRNA-seq) were used to analyze the heterogeneity of glioblastomas. After data integration, cell clustering, and annotation, five algorithms were used to calculate scores for exosome-related genes(ERGs). Cell trajectory analysis and intercellular communication analysis were performed to explore exosome-related communication patterns. Spatial transcriptome sequencing data were analyzed to validate the findings. To further utilize exosome-related features to aid in clinical decision-making, a prognostic model was constructed using GBM's bulk RNA-seq.

Results

Different cell subpopulations were observed in GBM, with Monocytes/macrophages and malignant cells in tumor samples showing higher exosome-related scores. After identifying differentially expressed ERGs in malignant cells, pseudotime analysis revealed the cellular status of malignant cells during development. Intercellular communication analysis highlighted signaling pathways and ligand-receptor interactions. Spatial transcriptome sequencing confirmed the high expression of exosome-related gene features in the tumor core region. A prognostic model based on six ERGs was shown to be predictive of overall survival and immunotherapy outcome in GBM patients. Finally, based on the results of scRNA-seq and prognostic modeling as well as a series of cell function experiments, BARD1 was identified as a novel target for the treatment of GBM.

Conclusion

This study provides a comprehensive understanding of the exosome-related features of GBM in both scRNA-seq and stRNA-seq, with malignant cells with higher exosome-related scores exhibiting stronger communication with Monocytes/macrophages. In terms of spatial data, highly scored malignant cells were also concentrated in the tumor core region. In bulk RNA-seq, patients with a high exosome-related index exhibited an immunosuppressive microenvironment, which was accompanied by a worse prognosis as well as immunotherapy outcomes. Prognostic models constructed using ERGs are expected to be independent prognostic indicators for GBM patients, with potential implications for personalized treatment strategies for GBM. Knockdown of BARD1 in GBM cell lines reduces the invasive and value-added capacity of tumor cells, and thus BARD1-positively expressing malignant cells are a risk factor for GBM patients.