Identification and experimental validation of cuproptosis regulatory program in a sepsis immune microenvironment through a combination of single-cell and bulk RNA sequencing.

Background: In spite of its high mortality rate and poor prognosis, the pathogenesis of sepsis is still incompletely understood. This study established a cuproptosis-based risk model to diagnose and predict the risk of sepsis. In addition, the cuproptosis-related genes were identified for targeted therapy. Methods: Single-cell sequencing analyses were used to characterize the cuproptosis activity score (CuAS) and intercellular communications in sepsis. Differential cuproptosis-related genes (CRGs) were identified in conjunction with single-cell and bulk RNA sequencing. LASSO and Cox regression analyses were employed to develop a risk model. Three external cohorts were conducted to assess the model's accuracy. Differences in immune infiltration, immune cell subtypes, pathway enrichment, and the expression of immunomodulators were further evaluated in distinct groups. Finally, various in-vitro experiments, such as flow cytometry, Western blot, and ELISA, were used to explore the role of LST1 in sepsis. Results: ScRNA-seq analysis demonstrated that CuAS was highly enriched in monocytes and was closely related to the poor prognosis of sepsis patients. Patients with higher CuAS exhibited prominent strength and numbers of cell-cell interactions. A total of five CRGs were identified based on the LASSO and Cox regression analyses, and a CRG-based risk model was established. The lower riskScore cohort exhibited enhanced immune cell infiltration, elevated immune scores, and increased expression of immune modulators, indicating the activation of an antibacterial response. Ultimately, in-vitro experiments demonstrated that LST1, a key gene in the risk model, was enhanced in the macrophage in response to LPS, which was closely related to the decrease of macrophage survival rate, the enhancement of apoptosis and oxidative stress injury, and the imbalance of the M1/M2 phenotype. Conclusions: This study constructed a cuproptosis-related risk model to accurately predict the prognosis of sepsis. We further characterized the cuproptosis-related gene LST1 to provide a theoretical framework for sepsis therapy.

Construction of a new prognostic model for colorectal cancer based on bulk RNA-seq combined with The Cancer Genome Atlas data.

Background: Colorectal cancer (CRC) is one of the leading causes of cancer-related deaths, and improving the prognosis of CRC patients is an urgent concern. The aim of this study was to explore new immunotherapy targets to improve survival in CRC patients. Methods: We analyzed CRC-related single-cell data GSE201348 from the Gene Expression Omnibus (GEO) database, and identified differentially expressed genes (DEGs). Subsequently, we performed differential analysis on the rectum adenocarcinoma (READ) and colon adenocarcinoma (COAD) transcriptome sequencing data [The Cancer Genome Atlas (TCGA)-CRC queue] and clinical data downloaded from TCGA database. Subgroup analysis was performed using CIBERSORTx and cluster analysis. Finally, biomarkers were identified by one-way cox regression as well as least absolute shrinkage and selection operator (LASSO) analysis. Results: In this study, we analyzed CRC-related single-cell data GSE201348, and identified 5,210 DEGs. Subsequently, we performed differential analysis on the TCGA-CRC queue database, and obtained 4,408 DEGs. Then, we categorized the cancer samples in the sequencing data into three groups (k1, k2, and k3), with significant differences observed between the k1 and k2 groups via survival analysis. Further differential analysis on the samples in the k1 and k2 groups identified 1,899 DEGs. A total of 77 DEGs were selected among those DEGs obtained from three differential analyses. Through subsequent Cox univariate analysis and LASSO analysis, seven biomarkers (RETNLB, CLCA4, UGT2A3, SULT1B1, CCL24, BMP5, and ATOH1) were identified and selected to establish a risk score (RS). Conclusions: To sum up, this study demonstrates the potential of the seven-gene prognostic risk model as instrumental variables for predicting the prognosis of CRC.

Exploring the immune landscape and drug prediction of an M2 tumor-associated macrophage-related gene signature in EGFR-negative lung adenocarcinoma.

BACKGROUND: Improving immunotherapy efficacy for EGFR-negative lung adenocarcinoma (LUAD) patients remains a critical challenge, and the therapeutic effect of immunotherapy is largely determined by the tumor microenvironment (TME). Tumor-associated macrophages (TAMs) are the top-ranked immune infiltrating cells in the TME, and M2-TAMs exert potent roles in tumor promotion and chemotherapy resistance. An M2-TAM-based prognostic signature was constructed by integrative analysis of single-cell RNA-seq (scRNA-seq) and bulk RNA-seq data to reveal the immune landscape and select drugs in EGFR-negative LUAD. METHODS: M2-TAM-based biomarkers were obtained from the intersection of bulk RNA-seq data and scRNA-seq data. After consensus clustering of EGFR-negative LUAD into different clusters based on M2-TAM-based genes, we compared the prognosis, clinical features, estimate scores, immune infiltration, and checkpoint genes among the clusters. Next, we combined univariate Cox and LASSO regression analyses to establish an M2-TAM-based prognostic signature. RESULTS: CCL20, HLA-DMA, HLA-DRB5, KLF4, and TMSB4X were verified as prognostic M2-like TAM-related genes by univariate Cox and LASSO regression analyses. IPS and TMB analyses revealed that the high-risk group responded better to common immunotherapy. CONCLUSION: The study shows the potential of the M2-like TAM-related gene signature in EGFR-negative LUAD, explores the immune landscape based on M2-like TAM-related genes, and predict immunotherapy response of patients with EGFR-negative LUAD, providing a new insight for individualized treatment.

Circadian Rhythm Disruption in Hepatocellular Carcinoma Investigated by Integrated Analysis of Bulk and Single-Cell RNA Sequencing Data.

Circadian rhythms are essential regulators of a multitude of physiological and behavioral processes, such as the metabolism and function of the liver. Circadian rhythms are crucial to liver homeostasis, as the liver is a key metabolic organ accountable for the systemic equilibrium of the body. Circadian rhythm disruption alone is sufficient to cause liver cancer through the maintenance of hepatic metabolic disorder. Although there is evidence linking CRD to hepatocarcinogenesis, the precise cellular and molecular mechanisms that underlie the circadian crosstalk that leads to hepatocellular carcinoma remain unknown. The expression of CRD-related genes in HCC was investigated in this study via bulk RNA transcriptomic analysis and single-cell sequencing. Dysregulated CRD-related genes are predominantly found in hepatocytes and fibroblasts, according to the findings. By using a combination of single-cell RNA sequencing and bulk RNA sequencing analyses, the dysregulated CRD-related genes ADAMTS13, BIRC5, IGFBP3, MARCO, MT2A, NNMT, and PGLYRP2 were identified. The survival analysis using the Kaplan-Meier method revealed a significant correlation between the expression levels of BIRC5 and IGFBP3 and the survival of patients diagnosed with HCC.

Mitochondrial RNA modification-based signature to predict prognosis of lower grade glioma: a multi-omics exploration and verification study.

Mitochondrial RNA modification (MRM) plays a crucial role in regulating the expression of key mitochondrial genes and promoting tumor metastasis. Despite its significance, comprehensive studies on MRM in lower grade gliomas (LGGs) remain unknown. Single-cell RNA-seq data (GSE89567) was used to evaluate the distribution functional status, and correlation of MRM-related genes in different cell types of LGG microenvironment. We developed an MRM scoring system by selecting potential MRM-related genes using LASSO regression analysis and the Random Survival Forest algorithm, based on multiple bulk RNA-seq datasets from TCGA, CGGA, GSE16011, and E-MTAB-3892. Analysis was performed on prognostic and immunological features, signaling pathways, metabolism, somatic mutations and copy number variations (CNVs), treatment responses, and forecasting of potential small-molecule agents. A total of 35 MRM-related genes were selected from the literature. Differential expression analysis of 1120 normal brain tissues and 529 LGGs revealed that 22 and 10 genes were upregulated and downregulated, respectively. Most genes were associated with prognosis of LGG. METLL8, METLL2A, TRMT112, and METTL2B were extensively expressed in all cell types and different cell cycle of each cell type. Almost all cell types had clusters related to mitochondrial RNA processing, ribosome biogenesis, or oxidative phosphorylation. Cell-cell communication and Pearson correlation analyses indicated that MRM may promoting the development of microenvironment beneficial to malignant progression via modulating NCMA signaling pathway and ICP expression. A total of 11 and 9 MRM-related genes were observed by LASSO and the RSF algorithm, respectively, and finally 6 MRM-related genes were used to establish MRM scoring system (TRMT2B, TRMT11, METTL6, METTL8, TRMT6, and TRUB2). The six MRM-related genes were then validated by qPCR in glioma and normal tissues. MRM score can predict the malignant clinical characteristics, abundance of immune infiltration, gene variation, clinical outcome, the enrichment of signaling pathways and metabolism. In vitro experiments demonstrated that silencing METTL8 significantly curbs glioma cell proliferation and enhances apoptosis. Patients with a high MRM score showed a better response to immunotherapies and small-molecule agents such as arachidonyl trifluoromethyl ketone, MS.275, AH.6809, tacrolimus, and TTNPB. These novel insights into the biological impacts of MRM within the glioma microenvironment underscore its potential as a target for developing precise therapies, including immunotherapeutic approaches.

Targeting BMP2 for therapeutic strategies against hepatocellular carcinoma.

OBJECTIVES: This study aimed to investigate the role of BMP2 in hepatocellular carcinoma (HCC) growth and metastasis using a dual approach combining single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq. METHODS: scRNA-seq data from the GEO database and bulk RNA-seq data from the TCGA database were analyzed. Differentially expressed marker genes of endothelial cells were identified and analyzed using enrichment analysis, PPI analysis, correlation analysis, and GSEA. In vitro, experiments were conducted using the Huh-7 HCC cell line, and in vivo, models of HCC growth and metastasis were established by knocking down BMP2. RESULTS: The scRNA-seq analysis identified BMP2 as a key marker gene in endothelial cells of HCC samples. Elevated BMP2 expression correlated with poor prognosis in HCC. In vitro experiments showed that silencing BMP2 inhibited the proliferation, migration, and invasion of liver cancer cells. In vivo studies confirmed increased BMP2 expression in HCC tissues, promoting angiogenesis and HCC growth. CONCLUSION: This study highlights the role of BMP2 in tumor angiogenesis and HCC progression. Targeting BMP2 could be a promising therapeutic strategy against HCC.

Based on scRNA-seq and bulk RNA-seq to establish tumor immune microenvironment-associated signature of skin melanoma and predict immunotherapy response.

Skin cutaneous melanoma (SKCM), a form of skin cancer, ranks among the most formidable and lethal malignancies. Exploring tumor microenvironment (TME)-based prognostic indicators would help improve the efficacy of immunotherapy for SKCM patients. This study analyzed SKCM scRNA-seq data to cluster non-malignant cells that could be used to explore the TME into nine immune/stromal cell types, including B cells, CD4 T cells, CD8 T cells, dendritic cells, endothelial cells, Fibroblasts, macrophages, neurons, and natural killer (NK) cells. Using data from The Cancer Genome Atlas (TCGA), we employed SKCM expression profiling to identify differentially expressed immune-associated genes (DEIAGs), which were then incorporated into weighted gene co-expression network analysis (WGCNA) to investigate TME-associated hub genes. Discover candidate small molecule drugs based on pivotal genes. Tumor immune microenvironment-associated genes (TIMAGs) for constructing TIMAS were identified and validated. Finally, the characteristics of TIAMS subgroups and the ability of TIMAS to predict immunotherapy outcomes were analyzed. We identified five TIMAGs (CD86, CD80, SEMA4D, C1QA, and IRF1) and used them to construct TIMAS. In addition, five potential SKCM drugs were identified. The results showed that TIMAS-low patients were associated with immune-related signaling pathways, high MUC16 mutation frequency, high T cell infiltration, and M1 macrophages, and were more favorable for immunotherapy. Collectively, TIMAS constructed by comprehensive analysis of scRNA-seq and bulk RNA-seq data is a promising marker for predicting ICI treatment outcomes and improving individualized therapy for SKCM patients.

Meningioma transcriptomic landscape demonstrates novel subtypes with regional associated biology and patient outcome.

Meningiomas, although mostly benign, can be recurrent and fatal. World Health Organization (WHO) grading of the tumor does not always identify high-risk meningioma, and better characterizations of their aggressive biology are needed. To approach this problem, we combined 13 bulk RNA sequencing (RNA-seq) datasets to create a dimension-reduced reference landscape of 1,298 meningiomas. The clinical and genomic metadata effectively correlated with landscape regions, which led to the identification of meningioma subtypes with specific biological signatures. The time to recurrence also correlated with the map location. Further, we developed an algorithm that maps new patients onto this landscape, where the nearest neighbors predict outcome. This study highlights the utility of combining bulk transcriptomic datasets to visualize the complexity of tumor populations. Further, we provide an interactive tool for understanding the disease and predicting patient outcomes. This resource is accessible via the online tool Oncoscape, where the scientific community can explore the meningioma landscape.

RNA-Seq of an LPS-Induced Inflammation Model Reveals Transcriptional Profile Patterns of Inflammatory Processes.

The LPS-induced inflammation model is widely used for studying inflammatory processes due to its cost-effectiveness, reproducibility, and faithful representation of key hallmarks. While researchers often validate this model using clinical cytokine markers, a comprehensive understanding of gene regulatory mechanisms requires extending investigation beyond these hallmarks. Our study leveraged multiple whole-blood bulk RNA-seq datasets to rigorously compare the transcriptional profiles of the well-established LPS-induced inflammation model with those of several human diseases characterized by systemic inflammation. Beyond conventional inflammation-associated systems, we explored additional systems indirectly associated with inflammatory responses (i.e., ISR, RAAS, and UPR) using a customized core inflammatory gene list. Our cross-condition-validation approach spanned four distinct conditions: systemic lupus erythematosus (SLE) patients, dengue infection, candidemia infection, and staphylococcus aureus exposure. This analysis approach, utilizing the core gene list aimed to assess the model's suitability for understanding the gene regulatory mechanisms underlying inflammatory processes triggered by diverse factors. Our analysis resulted in elevated expressions of innate immune-associated genes, coinciding with suppressed expressions of adaptive immune-associated genes. Also, upregulation of genes associated with cellular stresses and mitochondrial innate immune responses underscored oxidative stress as a central driver of the corresponding inflammatory processes in both the LPS-induced and other inflammatory contexts.

Error modelled gene expression analysis (EMOGEA) provides a superior overview of time course RNA-seq measurements and low count gene expression.

Temporal RNA-sequencing (RNA-seq) studies of bulk samples provide an opportunity for improved understanding of gene regulation during dynamic phenomena such as development, tumor progression or response to an incremental dose of a pharmacotherapeutic. Moreover, single-cell RNA-seq (scRNA-seq) data implicitly exhibit temporal characteristics because gene expression values recapitulate dynamic processes such as cellular transitions. Unfortunately, temporal RNA-seq data continue to be analyzed by methods that ignore this ordinal structure and yield results that are often difficult to interpret. Here, we present Error Modelled Gene Expression Analysis (EMOGEA), a framework for analyzing RNA-seq data that incorporates measurement uncertainty, while introducing a special formulation for those acquired to monitor dynamic phenomena. This method is specifically suited for RNA-seq studies in which low-count transcripts with small-fold changes lead to significant biological effects. Such transcripts include genes involved in signaling and non-coding RNAs that inherently exhibit low levels of expression. Using simulation studies, we show that this framework down-weights samples that exhibit extreme responses such as batch effects allowing them to be modeled with the rest of the samples and maintain the degrees of freedom originally envisioned for a study. Using temporal experimental data, we demonstrate the framework by extracting a cascade of gene expression waves from a well-designed RNA-seq study of zebrafish embryogenesis and an scRNA-seq study of mouse pre-implantation and provide unique biological insights into the regulation of genes in each wave. For non-ordinal measurements, we show that EMOGEA has a much higher rate of true positive calls and a vanishingly small rate of false negative discoveries compared to common approaches. Finally, we provide two packages in Python and R that are self-contained and easy to use, including test data.

Heterogeneity of immune cells and their communications unveiled by transcriptome profiling in acute inflammatory lung injury.

Background: Acute Respiratory Distress Syndrome (ARDS) or its earlier stage Acute lung injury (ALI), is a worldwide health concern that jeopardizes human well-being. Currently, the treatment strategies to mitigate the incidence and mortality of ARDS are severely restricted. This limitation can be attributed, at least in part, to the substantial variations in immunity observed in individuals with this syndrome. Methods: Bulk and single cell RNA sequencing from ALI mice and single cell RNA sequencing from ARDS patients were analyzed. We utilized the Seurat program package in R and cellmarker 2.0 to cluster and annotate the data. The differential, enrichment, protein interaction, and cell-cell communication analysis were conducted. Results: The mice with ALI caused by pulmonary and extrapulmonary factors demonstrated differential expression including Clec4e, Retnlg, S100a9, Coro1a, and Lars2. We have determined that inflammatory factors have a greater significance in extrapulmonary ALI, while multiple pathways collaborate in the development of pulmonary ALI. Clustering analysis revealed significant heterogeneity in the relative abundance of immune cells in different ALI models. The autocrine action of neutrophils plays a crucial role in pulmonary ALI. Additionally, there was a significant increase in signaling intensity between B cells and M1 macrophages, NKT cells and M1 macrophages in extrapulmonary ALI. The CXCL, CSF3 and MIF, TGFß signaling pathways play a vital role in pulmonary and extrapulmonary ALI, respectively. Moreover, the analysis of human single-cell revealed DCs signaling to monocytes and neutrophils in COVID-19-associated ARDS is stronger compared to sepsis-related ARDS. In sepsis-related ARDS, CD8+ T and Th cells exhibit more prominent signaling to B-cell nucleated DCs. Meanwhile, both MIF and CXCL signaling pathways are specific to sepsis-related ARDS. Conclusion: This study has identified specific gene signatures and signaling pathways in animal models and human samples that facilitate the interaction between immune cells, which could be targeted therapeutically in ARDS patients of various etiologies.

imply: improving cell-type deconvolution accuracy using personalized reference profiles.

Using computational tools, bulk transcriptomics can be deconvoluted to estimate the abundance of constituent cell types. However, existing deconvolution methods are conditioned on the assumption that the whole study population is served by a single reference panel, ignoring person-to-person heterogeneity. Here, we present imply, a novel algorithm to deconvolute cell type proportions using personalized reference panels. Simulation studies demonstrate reduced bias compared with existing methods. Real data analyses on longitudinal consortia show disparities in cell type proportions are associated with several disease phenotypes in Type 1 diabetes and Parkinson's disease. imply is available through the R/Bioconductor package ISLET at https://bioconductor.org/packages/ISLET/ .

scINRB: single-cell gene expression imputation with network regularization and bulk RNA-seq data.

Single-cell RNA sequencing (scRNA-seq) facilitates the study of cell type heterogeneity and the construction of cell atlas. However, due to its limitations, many genes may be detected to have zero expressions, i.e. dropout events, leading to bias in downstream analyses and hindering the identification and characterization of cell types and cell functions. Although many imputation methods have been developed, their performances are generally lower than expected across different kinds and dimensions of data and application scenarios. Therefore, developing an accurate and robust single-cell gene expression data imputation method is still essential. Considering to maintain the original cell-cell and gene-gene correlations and leverage bulk RNA sequencing (bulk RNA-seq) data information, we propose scINRB, a single-cell gene expression imputation method with network regularization and bulk RNA-seq data. scINRB adopts network-regularized non-negative matrix factorization to ensure that the imputed data maintains the cell-cell and gene-gene similarities and also approaches the gene average expression calculated from bulk RNA-seq data. To evaluate the performance, we test scINRB on simulated and experimental datasets and compare it with other commonly used imputation methods. The results show that scINRB recovers gene expression accurately even in the case of high dropout rates and dimensions, preserves cell-cell and gene-gene similarities and improves various downstream analyses including visualization, clustering and trajectory inference.

Integrated analysis of scRNA-seq and bulk RNA-seq identifies FBXO2 as a candidate biomarker associated with chemoresistance in HGSOC.

Background: High-grade serous ovarian carcinoma (HGSOC) is the most prevalent and aggressive histological subtype of epithelial ovarian cancer. Around 80% of individuals will experience a recurrence within five years because of resistance to chemotherapy, despite initially responding well to platinum-based treatment. Biomarkers associated with chemoresistance are desperately needed in clinical practice. Methods: We jointly analyzed the transcriptomic profiles of single-cell and bulk datasets of HGSOC to identify cell types associated with chemoresistance. Copy number variation (CNV) inference was performed to identify malignant cells. We subsequently analyzed the expression of candidate biomarkers and their relationship with patients' prognosis. The enrichment analysis and potential biological function of candidate biomarkers were explored. Then, we validated the candidate biomarker using in vitro experiments. Results: We identified 8871 malignant epithelial cells in a single-cell RNA sequencing dataset, of which 861 cells were associated with chemoresistance. Among these malignant epithelial cells, FBXO2 (F-box protein 2) is highly expressed in cells related to chemoresistance. Moreover, FBXO2 expression was found to be higher in epithelial cells from chemoresistance samples compared to those from chemosensitivity samples in a separate single-cell RNA sequencing dataset. Patients exhibiting elevated levels of FBXO2 experienced poorer outcomes in terms of both overall survival (OS) and progression-free survival (PFS). FBXO2 could impact chemoresistance by influencing the PI3K-Akt signaling pathway, focal adhesion, and ECM-receptor interactions and regulating tumorigenesis. The 50% maximum inhibitory concentration (IC50) of cisplatin decreased in A2780 and SKOV3 ovarian carcinoma cell lines with silenced FBXO2 during an in vitro experiment. Conclusions: We determined that FBXO2 is a potential biomarker linked to chemoresistance in HGSOC by combining single-cell RNA-seq and bulk RNA-seq dataset. Our results suggest that FBXO2 could serve as a valuable prognostic marker and potential target for drug development in HGSOC.

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.

A novel pyroptosis-related gene signature exhibits distinct immune cells infiltration landscape in Wilms' tumor.

BACKGROUND: Wilms' tumor (WT) is the most common renal tumor in childhood. Pyroptosis, a type of inflammation-characterized and immune-related programmed cell death, has been extensively studied in multiple tumors. In the current study, we aim to construct a pyroptosis-related gene signature for predicting the prognosis of Wilms' tumor. METHODS: We acquired RNA-seq data from TARGET kidney tumor projects for constructing a gene signature, and snRNA-seq data from GEO database for validating signature-constructing genes. Pyroptosis-related genes (PRGs) were collected from three online databases. We constructed the gene signature by Lasso Cox regression and then established a nomogram. Underlying mechanisms by which gene signature is related to overall survival states of patients were explored by immune cell infiltration analysis, differential expression analysis, and functional enrichment analysis. RESULTS: A pyroptosis-related gene signature was constructed with 14 PRGs, which has a moderate to high predicting capacity with 1-, 3-, and 5-year area under the curve (AUC) values of 0.78, 0.80, and 0.83, respectively. A prognosis-predicting nomogram was established by gender, stage, and risk score. Tumor-infiltrating immune cells were quantified by seven algorithms, and the expression of CD8( +) T cells, B cells, Th2 cells, dendritic cells, and type 2 macrophages are positively or negatively correlated with risk score. Two single nuclear RNA-seq samples of different histology were harnessed for validation. The distribution of signature genes was identified in various cell types. CONCLUSIONS: We have established a pyroptosis-related 14-gene signature in WT. Moreover, the inherent roles of immune cells (CD8( +) T cells, B cells, Th2 cells, dendritic cells, and type 2 macrophages), functions of differentially expressed genes (tissue/organ development and intercellular communication), and status of signaling pathways (proteoglycans in cancer, signaling pathways regulating pluripotent of stem cells, and Wnt signaling pathway) have been elucidated, which might be employed as therapeutic targets in the future.

Multi-omics analysis reveals the association between specific solute carrier proteins gene expression patterns and the immune suppressive microenvironment in glioma.

Glioma is the most prevalent malignant brain tumour. Currently, reshaping its tumour microenvironment has emerged as an appealing strategy to enhance therapeutic efficacy. As the largest group of transmembrane transport proteins, solute carrier proteins (SLCs) are responsible for the transmembrane transport of various metabolites and ions. They play a crucial role in regulating the metabolism and functions of malignant cells and immune cells within the tumour microenvironment, making them a promising target in cancer therapy. Through multidimensional data analysis and experimental validation, we investigated the genetic landscape of SLCs in glioma. We established a classification system comprising 7-SLCs to predict the prognosis of glioma patients and their potential responses to immunotherapy and chemotherapy. Our findings unveiled specific SLC expression patterns and their correlation with the immune-suppressive microenvironment and metabolic status. The 7-SLC classification system was validated in distinguishing subgroups within the microenvironment, specifically identifying subsets involving malignant cells and tumour-associated macrophages. Furthermore, the orphan protein SLC43A3, a core member of the 7-SLC classification system, was identified as a key facilitator of tumour cell proliferation and migration, suggesting its potential as a novel target for cancer therapy.

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.

Comprehensive analysis of scRNA-seq and bulk RNA-seq reveal the characteristics of disulfidptosis and a prognostic signature in BLCA.

Disulfidptosis is a newly discovered mode of cell death. However, its biological mechanism in bladder cancer (BLCA) is still uncharacterized. In this investigation, we firstly examined the expression and mutation of disulfidptosis-related genes (DRGs) in BLCA. Two disulfidptosis phenotypes associated with DRGs expression patterns and immune cell infiltration were built. A disulfidptosis risk score signature was constructed based on ten differentially expressed genes (DEGs) between the disulfidptosis subtypes, which allowed patients to be stratified into high- and low-risk groups. We further confirmed that the disulfidptosis risk score signature has great power to predict prognosis, immune cell infiltration, and immunotherapy efficacy in BLCA. Additionally, we analyzed the differences in therapeutic sensitivities between high- and low-risk groups concerning targeted inhibitor therapy and immunotherapy. Analysis of single-cell RNA sequencing was conducted of the ten hub DRGs. Of the ten genes, we found that DUSP2 and SLCO1B3 were differentially expressed in BLCA tissues and adjacent normal tissues, and were markedly associated with patients' prognosis. Functional experiments revealed that overexpression of DUSP2 or knockdown of SLCO1B3 significantly inhibited cell proliferation, migration, and invasion in BLCA cells. In all, we present a fresh disulfidptosis-related prognostic signature, which has a remarkable capacity to characterize the immunological landscape and prognosis of BLCA patients.

Identification and validation of a novel signature based on macrophage marker genes for predicting prognosis and drug response in kidney renal clear cell carcinoma by integrated analysis of single cell and bulk RNA sequencing.

Macrophages are found in a variety of tumors and play a critical role in shaping the tumor microenvironment, affecting tumor progression, metastasis, and drug resistance. However, the clinical relevance of marker genes associated with macrophage in kidney renal clear cell carcinoma (KIRC) has yet to be documented. In this study, we initiated a thorough examination of single-cell RNA sequencing (scRNA-seq) data for KIRC retrieved from the Gene Expression Omnibus (GEO) database and determined 244 macrophage marker genes (MMGs). Univariate analysis, LASSO regression, and multivariate regression analysis were performed to develop a five-gene prognostic signature in The Cancer Genome Atlas (TCGA) database, which could divide KIRC patients into low-risk (L-R) and high-risk (H-R) groups. Then, a nomogram was constructed to predict the survival rate of KIRC patients at 1, 3, and 5 years, which was well assessed by receiver operating characteristic curve (ROC), calibration curve, and decision curve analyses (DCA). Functional enrichment analysis showed that immune-related pathways (such as immunoglobulin complex, immunoglobulin receptor binding, and cytokine-cytokine receptor interaction) were mainly enriched in the H-R group. Additionally, in comparison to the L-R cohort, patients belonging to the H-R cohort exhibited increased immune cell infiltration, elevated expression of immune checkpoint genes (ICGs), and a higher tumor immune dysfunction and exclusion (TIDE) score. This means that patients in the H-R group may be less sensitive to immunotherapy than those in the L-R group. Finally, IFI30 was validated to increase the ability of KIRC cells to proliferate, invade and migrate in vitro. In summary, our team has for the first time developed and validated a predictive model based on macrophage marker genes to accurately predict overall survival (OS), immune characteristics, and treatment benefit in KIRC patients.