Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular biology at an unprecedented resolution, enabling the characterization of cellular heterogeneity, identification of rare but significant cell types, and exploration of cell–cell communications and interactions. Its broad applications span both basic and clinical research domains. In this comprehensive review, we survey the current landscape of scRNA-seq analysis methods and tools, focusing on count modeling, cell-type annotation, data integration, including spatial transcriptomics, and the inference of cell–cell communication. We review the challenges encountered in scRNA-seq analysis, including issues of sparsity or low expression, reliability of cell annotation, and assumptions in data integration, and discuss the potential impact of suboptimal clustering and differential expression analysis tools on downstream analyses, particularly in identifying cell subpopulations. Finally, we discuss recent advancements and future directions for enhancing scRNA-seq analysis. Specifically, we highlight the development of novel tools for annotating single-cell data, integrating and interpreting multimodal datasets covering transcriptomics, epigenomics, and proteomics, and inferring cellular communication networks. By elucidating the latest progress and innovation, we provide a comprehensive overview of the rapidly advancing field of scRNA-seq analysis.

Acute myeloid leukemia (AML) is a hematological malignancy that is characterized by an expansion of immature myeloid precursors. Despite therapeutic advances, the prognosis of AML patients remains poor and there is a need for the evaluation of promising therapeutic candidates to treat the disease. The objective of this study was to evaluate the efficacy of duocarmycin Stable A (DSA) in AML cells in vitro. We hypothesized that DSA would induce DNA damage in the form of DNA double-strand breaks (DSBs) and exert cytotoxic effects on AML cells within the picomolar range. Human AML cell lines Molm-14 and HL-60 were used to perform 3-(4,5-dimethylthiazolyl-2)-2,5-diphenyltetrazolium bromide (MTT), DNA DSBs, cell cycle, 5-ethynyl-2-deoxyuridine (EdU), colony formation unit (CFU), Annexin V, RNA sequencing and other assays described in this study. Our results showed that DSA induced DNA DSBs, induced cell cycle arrest at the G2M phase, reduced proliferation and increased apoptosis in AML cells. Additionally, RNA sequencing results showed that DSA regulates genes that are associated with cellular processes such as DNA repair, G2M checkpoint and apoptosis. These results suggest that DSA is efficacious in AML cells and is therefore a promising potential therapeutic candidate that can be further evaluated for the treatment of AML.

The H3K27me2/me3 histone demethylase KDM6B is over-expressed in neuroblastoma and is essential to neuroblastoma cell survival. While the KDM6B inhibitor, GSK-J4, has shown activity in in vitro and in vivo preclinical models, the mechanism of action remains poorly defined. We demonstrate that genetic and pharmacologic inhibition of KDM6B downregulates the pRB-E2F transcriptome and MYCN expression. Chemical genetic analyses show that high expression of the E2F transcriptome is positively correlated with sensitivity of cancer cells to GSK-J4. Mechanistically, inhibition of KDM6B activity 1) reduces the chromatin accessibility of E2F target genes and MYCN, 2) selectively leads to an increase of H3K27me3 but a decrease of the enhancer mark H3K4me1 at the CTCF and BORIS binding sites, which may, consequently, disrupt the long-range chromatin interaction of MYCN and E2F target genes, and 3) phenocopies the transcriptome induced by the specific CDK4/6 inhibitor palbociclib. Overexpression of CDK4/6 or Rb1 knockout confers neuroblastoma cell resistance to both palbociclib and GSK-J4. A gene signature targeted by KDM6B inhibition is associated with poor survival of patients with neuroblastoma regardless of the MYCN status. These data indicate that KDM6B activity promotes an oncogenic CDK4/6-pRB-E2F pathway in neuroblastoma cells via H3K27me3-dependent enhancer-promoter interactions, providing a rationale to target KDM6B for high-risk neuroblastoma.

Neuroblastoma accounts for approximately 15% of all pediatric cancer-related deaths, largely due to disease relapse following intensive multimodal therapy. A critical barrier to cure neuroblastoma is the emergence of therapy-resistant tumor cells. Neuroblastoma comprises two major cell states, adrenergic (ADRN) and mesenchymal (MES), which are believed to interconvert and contribute to therapeutic resistance through lineage plasticity. To investigate the mechanisms underlying this plasticity, we subjected human and murine neuroblastoma models to repeated treatment with indisulam, a molecular glue compound that selectively degrades the splicing factor RBM39, until full drug resistance emerged. We then generated datasets from these models, including bulk transcriptomic data, ATAC-seq, and H3K27ac CUT&Tag. These data comprehensively characterize transcriptomic and epigenetic landscapes of resistant ADRN and MES neuroblastoma cell states. We present this resource to facilitate reuse by the scientific community. These datasets may support efforts to decipher lineage switching, identify regulators of therapy resistance, and discover potential therapeutic vulnerabilities in resistant neuroblastoma.