Non-alcoholic Fatty Liver Disease (NAFLD), noted for its widespread prevalence among adults, has become the leading chronic liver condition globally. Simultaneously, the annual disease burden, particularly liver cirrhosis caused by NAFLD, has increased significantly. Neutrophil Extracellular Traps (NETs) play a crucial role in the progression of this disease and are key to the pathogenesis of NAFLD. However, research into the specific roles of NETs-related genes in NAFLD is still a field requiring thorough investigation. Utilizing techniques like AddModuleScore, ssGSEA, and WGCNA, our team conducted gene screening to identify the genes linked to NETs in both single-cell and bulk transcriptomics. Using algorithms including Random Forest, Support Vector Machine, Least Absolute Shrinkage, and Selection Operator, we identified ZFP36L2 and PHLDA1 as key hub genes. The pivotal role of these genes in NAFLD diagnosis was confirmed using the training dataset GSE164760. This study identified 116 genes linked to NETs across single-cell and bulk transcriptomic analyses. These genes demonstrated enrichment in immune and metabolic pathways. Additionally, two NETs-related hub genes, PHLDA1 and ZFP36L2, were selected through machine learning for integration into a prognostic model. These hub genes play roles in inflammatory and metabolic processes. scRNA-seq results showed variations in cellular communication among cells with different expression patterns of these key genes. In conclusion, this study explored the molecular characteristics of NETs-associated genes in NAFLD. It identified two potential biomarkers and analyzed their roles in the hepatic microenvironment. These discoveries could aid in NAFLD diagnosis and management, with the ultimate goal of enhancing patient outcomes.

Background: Studies on immunogenic death (ICD) in lung adenocarcinoma are limited, and this study aimed to determine the function of ICD in LUAD and to construct a novel ICD-based prognostic model to improve immune efficacy in lung adenocarcinoma patients. Methods: The data for lung adenocarcinoma were obtained from the Cancer Genome Atlas (TCGA) database and the National Center for Biotechnology Information (GEO). The single-cell data were obtained from Bischoff P et al. To identify subpopulations, we performed descending clustering using TSNE. We collected sets of genes related to immunogenic death from the literature and identified ICD-related genes through gene set analysis of variance (GSVA) and weighted gene correlation network analysis (WGCNA). Lung adenocarcinoma patients were classified into two types using consistency clustering. The difference between the two types was analyzed to obtain differential genes. An immunogenic death model (ICDRS) was established using LASSO-Cox analysis and compared with lung adenocarcinoma models of other individuals. External validation was performed in the GSE31210 and GSE50081 cohorts. The efficacy of immunotherapy was assessed using the TIDE algorithm and the IMvigor210, GSE78220, and TCIA cohorts. Furthermore, differences in mutational profiles and immune microenvironment between different risk groups were investigated. Subsequently, ROC diagnostic curves and KM survival curves were used to screen ICDRS key regulatory genes. Finally, RT-qPCR was used to verify the differential expression of these genes. Results: Eight ICD genes were found to be highly predictive of LUAD prognosis and significantly correlated with it. Multivariate analysis showed that patients in the low-risk group had a higher overall survival rate than those in the high-risk group, indicating that the model was an independent predictor of LUAD. Additionally, ICDRS demonstrated better predictive ability compared to 11 previously published models. Furthermore, significant differences in biological function and immune cell infiltration were observed in the tumor microenvironment between the high-risk and low-risk groups. It is noteworthy that immunotherapy was also significant in both groups. These findings suggest that the model has good predictive efficacy. Conclusions: The ICD model demonstrated good predictive performance, revealing the tumor microenvironment and providing a new method for evaluating the efficacy of pre-immunization. This offers a new strategy for future treatment of lung adenocarcinoma.

Bone marrow-derived mesenchymal stem cells (BMSCs) exhibit multi-lineage differentiation potential and robust proliferative capacity. The late stage of differentiation signifies the functional maturation and characterization of specific cell lineages, which is crucial for studying lineage-specific differentiation mechanisms. However, the molecular processes governing late-stage BMSC differentiation remain poorly understood. This study aimed to elucidate the key biological processes involved in late-stage BMSC differentiation. Publicly available transcriptomic data from human BMSCs were analyzed after approximately 14 days of osteogenic, adipogenic, and chondrogenic differentiation. Thirty-one differentially expressed genes (DEGs) associated with differentiation were identified. Pathway enrichment analysis indicated that the DEGs were involved in extracellular matrix (ECM)-receptor interactions, focal adhesion, and glycolipid biosynthesis, a ganglion series process. Subsequently, the target genes were validated using publicly available single-cell RNA-seq data from mouse BMSCs. Lamc1 exhibited predominant distribution in adipocytes and osteoblasts, primarily during the G2/M phase. Tln2 and Hexb were expressed in chondroblasts, osteoblasts, and adipocytes, while St3gal5 was abundantly distributed in stem cells. Cell communication analysis identified two receptors that interact with LAMCI. q-PCR results confirmed the upregulation of Lamc1, Tln2, Hexb, and St3gal5 during osteogenic differentiation and their downregulation during adipogenic differentiation. Knockdown of Lamc1 inhibited adipogenic and osteogenic differentiation. In conclusion, this study identified four genes, Lamc1, Tln2, Hexb, and St3gal5, that may play important roles in the late-stage differentiation of BMSCs. It elucidated their interactions and the pathways they influence, providing a foundation for further research on BMSC differentiation.

OBJECTIVE: Ovarian metastasis of gastric cancer (GC), commonly referred to as Krukenberg tumors, leads to a poor prognosis. However, the cause of metastasis remains unknown. Here, we present an integrated single-cell RNA sequencing (scRNA-Seq) analysis of the immunological microenvironment of two paired clinical specimens with ovarian metastasis of GC.
METHODS: scRNA-Seq was performed to determine the immunological microenvironment in ovarian metastasis of gastric cancer. CellChat was employed to analyze cell-cell communications across different cell types. Functional enrichment analysis was done by enrichKEGG in clusterProfiler. GEPIA2 was used to assess the influence of certain genes and gene signatures on prognosis.
RESULTS: The ovarian metastasis tissues exhibit a heterogenous immunological microenvironment compared to the primary tumors. Exhaustion of T and B cells is observed in the ovarian metastasis tissues. Compared to the paired adjacent non-tumoral and primary tumors, the ratio of endothelial cells and fibroblasts is high in the ovarian metastasis tissues. Compared to primary ovarian cancers, we identify a specific group of tumor-associated fibroblasts with MFAP4 and CAPNS1 expression in the ovarian metastatic tissues of GC. We further define metastasis-related-endothelial and metastasis-related-fibroblast signatures and indicate that patients with these high signature scores have a poor prognosis. In addition, the ovarian metastasis tissue has a lower level of intercellular communications compared to the primary tumor.
CONCLUSIONS: Our findings reveal the immunological microenvironment of ovarian metastasis of gastric cancer and will promote the discovery of new therapeutic strategies for ovarian metastasis in gastric cancer.

Metastatic clear cell renal cell carcinoma has heterogenous tumor microenvironment (TME). Among the metastatic lesions, pancreas metastasis is rare and controversy in treatment approaches. Here, extensive primary and metastatic lesion samples were included by single-cell RNA-seq to decipher the distinct metastasis TME. The hypoxic and inflammatory TME of pancreas metastasis was decoded in this study, and the activation of PAX8-myc signaling, and metabolic reprogramming were observed. The active components including endothelial cells, fibroblasts and T cells were profiled. Meanwhile, we also evaluated the effect of anti-angiogenesis treatment in the pancreas metastasis patient. The potential mechanisms of pancreatic tropism, instability of genome, and the response of immunotherapy were also discussed in this work. Taken together, our findings suggest a clue to the heterogeneity in metastasis TME and provide evidence for the treatment of pancreas metastasis in renal cell carcinoma patients.

Our understanding of human fetal cerebellum development during the late second trimester, a critical period for the generation of astrocytes, oligodendrocytes, and unipolar brush cells (UBCs), remains limited. Here, we performed single-cell RNA sequencing (scRNA-seq) in human fetal cerebellum samples from gestational weeks (GWs) 18-25. We find that proliferating UBC progenitors distribute in the subventricular zone of the rhombic lip (RLSVZ) near white matter (WM), forming a layer structure. We also delineate two trajectories from astrogenic radial glia (ARGs) to Bergmann glial progenitors (BGPs) and recognize oligodendrogenic radial glia (ORGs) as one source of primitive oligodendrocyte progenitor cells (PriOPCs). Additionally, our scRNA-seq analysis of the trisomy 21 fetal cerebellum at this stage reveals abnormal upregulated genes in pathways such as the cell adhesion pathway and focal adhesion pathway, which potentially promote neuronal differentiation. Overall, our research provides valuable insights into normal and abnormal development of the human fetal cerebellum.

OBJECTIVE: Osteosarcoma, a highly malignant primary bone tumor primarily affecting adolescents, frequently develops resistance to initial chemotherapy, leading to metastasis and limited treatment options. Our study aims to uncover novel therapeutic targets for metastatic and recurrent osteosarcoma.
METHODS: In this study, we proved the potential of modulating the YAP1-regulated glutamine metabolic pathway to augment the response of OS to DFMO. We initially employed single-cell transcriptomic data to gauge the activation level of polyamine metabolism in MTAP-deleted OS patients. This was further substantiated by transcriptome sequencing data from recurrent and non-recurrent patient tissues, confirming the activation of polyamine metabolism in progressive OS. Through high-throughput drug screening, we pinpointed CIL56, a YAP1 inhibitor, as a promising candidate for a combined therapeutic strategy with DFMO. In vivo, we utilized PDX and CDX models to validate the therapeutic efficacy of this drug combination. In vitro, we conducted western blot analysis, qPCR analysis, immunofluorescence staining, and PuMA experiments to monitor alterations in molecular expression, distribution, and tumor metastasis capability. We employed CCK-8 and colony formation assays to assess the proliferative capacity of cells in the experimental group. We used flow cytometry and reactive oxygen probes to observe changes in ROS and glutamine metabolism within the cells. Finally, we applied RNA-seq in tandem with metabolomics to identify metabolic alterations in OS cells treated with a DFMO and CIL56 combination. This enabled us to intervene and validate the role of the YAP1-mediated glutamine metabolic pathway in DFMO resistance.
RESULTS: Through single-cell RNA-seq data analysis, we pinpointed a subset of late-stage OS cells with significantly upregulated polyamine metabolism. This upregulation was further substantiated by transcriptomic profiling of recurrent and non-recurrent OS tissues. High-throughput drug screening revealed a promising combination strategy involving DFMO and CIL56. DFMO treatment curbs the phosphorylation of YAP1 protein in OS cells, promoting nuclear entry and initiating the YAP1-mediated glutamine metabolic pathway. This reduces intracellular ROS levels, countering DFMO\'s anticancer effect. The therapeutic efficacy of DFMO can be amplified both in vivo and in vitro by combining it with the YAP1 inhibitor CIL56 or the glutaminase inhibitor CB-839. This underscores the significant potential of targeting the YAP1-mediated glutamine metabolic pathway to enhance efficacy of DFMO.
CONCLUSIONS: Our findings elucidate YAP1-mediated glutamine metabolism as a crucial bypass mechanism against DFMO, following the inhibition of polyamine metabolism. Our study provides valuable insights into the potential role of DFMO in an \"One-two Punch\" therapy of metastatic and recurrent osteosarcoma.

BACKGROUND: Breast cancer (BC) is a heterogeneous disease, with the ductal subtype exhibiting significant cellular diversity that influences prognosis and response to treatment. Single-cell RNA sequencing data from the GEO database were utilized in this study to investigate the underlying mechanisms of cellular heterogeneity and to identify potential prognostic markers and therapeutic targets.
METHODS: Bioinformatics analysis was conducted using R packages to analyze the single-cell sequencing data. The presence of highly variable genes and differences in malignant potency within the same BC samples were examined. Differential gene expression and biological function between Type 1 and Type 2 ductal epithelial cells were identified. Lasso regression and Cox proportional hazards regression analyses were employed to identify genes associated with patient prognosis. Experimental validation was performed in vitro and in vivo to confirm the functional relevance of the identified genes.
RESULTS: The analysis revealed notable heterogeneity among BC cells, with the presence of highly variable genes and differences in malignant behavior within the same samples. Significant disparities in gene expression and biological function were identified between Type 1 and Type 2 ductal epithelial cells. Through regression analyses, CYP24A1 and TFPI2 were identified as pivotal genes associated with patient prognosis. Kaplan-Meier curves demonstrated their prognostic significance, and experimental validation confirmed their inhibitory effects on malignant behaviors of ductal BC cells.
CONCLUSIONS: This study highlights the cellular heterogeneity in ductal subtype breast cancer and delineates the differential gene expressions and biological functions between Type 1 and Type 2 ductal epithelial cells. The genes CYP24A1 and TFPI2 emerged as promising prognostic markers and therapeutic targets, exhibiting inhibitory effects on BC cell malignancy in vitro and in vivo. These findings offer the potential for improved BC management and the development of targeted treatment strategies.

Single-cell RNA sequencing (scRNA-seq) is the sequencing technology of a single cell whose expression reflects the overall characteristics of the individual cell, facilitating the research of problems at the cellular level. However, the problems of scRNA-seq such as dimensionality reduction processing of massive data, technical noise in data, and visualization of single-cell type clustering cause great difficulties for analyzing and processing scRNA-seq data. In this paper, we propose a new single-cell data analysis model using denoising autoencoder and multi-type graph neural networks (scDMG), which learns cell-cell topology information and latent representation of scRNA-seq data. scDMG introduces the zero-inflated negative binomial (ZINB) model into a denoising autoencoder (DAE) to perform dimensionality reduction and denoising on the raw data. scDMG integrates multiple-type graph neural networks as the encoder to further train the preprocessed data, which better deals with various types of scRNA-seq datasets, resolves dropout events in scRNA-seq data, and enables preliminary classification of scRNA-seq data. By employing TSNE and PCA algorithms for the trained data and invoking Louvain algorithm, scDMG has better dimensionality reduction and clustering optimization. Compared with other mainstream scRNA-seq clustering algorithms, scDMG outperforms other state-of-the-art methods in various clustering performance metrics and shows better scalability, shorter runtime, and great clustering results.

Gene set scoring (GSS) has been routinely conducted for gene expression analysis of bulk or single-cell RNA sequencing (RNA-seq) data, which helps to decipher single-cell heterogeneity and cell type-specific variability by incorporating prior knowledge from functional gene sets. Single-cell assay for transposase accessible chromatin using sequencing (scATAC-seq) is a powerful technique for interrogating single-cell chromatin-based gene regulation, and genes or gene sets with dynamic regulatory potentials can be regarded as cell type-specific markers as if in single-cell RNA-seq (scRNA-seq). However, there are few GSS tools specifically designed for scATAC-seq, and the applicability and performance of RNA-seq GSS tools on scATAC-seq data remain to be investigated. Here, we systematically benchmarked ten GSS tools, including four bulk RNA-seq tools, five scRNA-seq tools, and one scATAC-seq method. First, using matched scATAC-seq and scRNA-seq datasets, we found that the performance of GSS tools on scATAC-seq data was comparable to that on scRNA-seq, suggesting their applicability to scATAC-seq. Then, the performance of different GSS tools was extensively evaluated using up to ten scATAC-seq datasets. Moreover, we evaluated the impact of gene activity conversion, dropout imputation, and gene set collections on the results of GSS. Results show that dropout imputation can significantly promote the performance of almost all GSS tools, while the impact of gene activity conversion methods or gene set collections on GSS performance is more dependent on GSS tools or datasets. Finally, we provided practical guidelines for choosing appropriate preprocessing methods and GSS tools in different application scenarios.