Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “single cell RNA sequencing”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources

Next-generation sequencing in breast cancer: current clinical applications and future directions.

INTRODUCTION: Breast cancer is a heterogeneous disease that claims 670,000 lives by 2022. Omic technologies, particularly next generation sequencing (NGS) offers promising avenues for precision medicine. American Society of Clinical Oncology (ASCO) outlines genomic testing's utility, emphasizing prognostic and diagnostic potential. OBJECTIVES: This review succinctly explores NGS's evolution and clinical applications of NGS in breast cancer, thereby guiding future research to enhance patient care. METHODS: Comprehensive literature searches were conducted using databases such as PubMed, Google Scholar, and ResearchGate, focusing on keywords including breast cancer, HER-2 low breast cancer, circulating tumour DNA, single-cell RNA sequencing, and next-generation sequencing. Peer-reviewed, high-quality articles published in English were selected for inclusion. RESULTS: Previous studies have explored the evolution of NGS technology and its clinical applications in breast cancer, including genomic and transcriptomic characterization, treatment guidance, and resistance prediction. Molecular profiling of challenging entities such as early-onset breast cancer and HER-2 low tumours was summarized, with key findings highlighted. This review also discusses emerging technologies, including circulating DNA and single-cell sequencing, as promising avenues for discovery. CONCLUSION: NGS has revealed the genomic and transcriptomic diversity of breast cancer, identifying actionable alterations associated with chemotherapy response and resistance to therapies such as trastuzumab, TKIs, and CDK4/6 inhibitors. Circulating tumour DNA (ctDNA) shows potential for diagnosis, prediction, prognosis, and monitoring, despite tumour heterogeneity. Single-cell analysis enables exploration of individual cell transcriptomes, though high costs and low throughput remain barriers to widespread adoption. HER2-low tumours continue to pose significant research challenges.

Humans↗

A Functionally Constrained Immune Ecosystem in Microsatellite-stable Colorectal Cancer Resolved by Single-cell and Exome Profiling.

BACKGROUND/AIM: Microsatellite-stable (MSS) colorectal cancer (CRC) generally responds poorly to immune checkpoint blockade, but some MSS tumors are T-cell rich. We examined whether such infiltration reflected effective immunity or functional immune constraint. CASE REPORT: A 77-year-old woman underwent resection of a mismatch repair-proficient (pMMR), MSS, low-mutational-burden CRC with a synchronous adenoma. Whole-exome sequencing of tumor, adenoma and adjacent normal tissue detected no shared high-confidence somatic mutations between tumor and adenoma within the sensitivity of this WES analysis and identified tumor-specific APC, KRAS and TP53 alterations. Tumor single-cell RNA sequencing yielded 7,569 cells, with T-lineage populations comprising 83.5%. Cytotoxic T cells showed cytolytic and dysfunction-associated features, regulatory T cells (Tregs) showed suppressive remodeling, and Th17 cells showed inflammatory/profibrotic programs. CellChat nominated stromal MIF/FN1-CD74/CD44 and extracellular-matrix communication with T-cell compartments. CONCLUSION: This molecular case report shows that T-cell abundance and immune effectiveness can be uncoupled in MSS CRC.

Humans↗

Coordinated inflammatory macrophage and vascular smooth muscle cell remodeling signatures in human atherosclerosis: An integrative single-cell and bulk transcriptomic analysis.

Atherosclerotic plaque progression is shaped by coordinated inflammatory and remodeling programs involving immune cells and vascular wall cells. Inflammatory macrophage activation and vascular smooth muscle cell (VSMC) phenotypic remodeling are central features of human atherosclerosis, but their transcriptomic relationships during plaque progression remain incompletely characterized. This study integrated single-cell and bulk transcriptomic datasets to examine highly inflammatory macrophage states, VSMC remodeling-related transcriptional programs, and candidate ligand-receptor expression patterns in human atherosclerotic plaques. Human atherosclerotic plaque single-cell RNA sequencing data from GSE260657 and bulk transcriptomic data from GSE28829 were analyzed. After quality control, 7628 cells were retained for single-cell analysis. Major cell types were annotated using canonical markers, followed by reclustering of macrophages and VSMC-related cells. Functional module scoring, differential expression analysis, Gene Ontology biological process enrichment, and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to characterize macrophage transcriptional states. Slingshot was applied to infer VSMC pseudotime ordering. CellChat and NicheNet were used to prioritize candidate ligand-receptor expression patterns and ligand-associated VSMC target gene programs. External bulk transcriptomic analysis was performed to examine whether single-cell-derived inflammatory and remodeling signatures were represented at the tissue-transcriptome level during plaque progression. Macrophage reclustering identified a highly inflammatory macrophage state characterized by prominent inflammatory activation, cytokine-response, and stress-response features. Genes upregulated in this population were enriched in pathways related to tumor necrosis factor (TNF) response, nuclear factor kappa B signaling, leukocyte activation, cytokine signaling, lipid and atherosclerosis, toll-like receptor signaling, and inflammasome-associated inflammation. VSMC reclustering revealed contractile VSMCs, PTHLH+ synthetic VSMCs, KRT7+ VSMC-like cells, interferon-responsive VSMCs, pericyte-like mural cells, and osteogenic/modulated VSMCs. Pseudotime analysis showed a broad contractile-to-osteogenic/modulated transcriptional continuum accompanied by increased expression of remodeling-associated genes and selected inflammatory or remodeling-associated receptor genes. CellChat and NicheNet analyses prioritized candidate ligand-receptor and ligand-associated target gene expression patterns involving SPP1-CD44, TNF-TNFRSF1A, IL1B-IL1R1/IL1RAP, MIF-ACKR3, PDGFB-PDGFRB, and FN1-SDC1/ITGB1. In GSE28829, inflammatory macrophage-, osteogenic/modulated VSMC-, candidate ligand-receptor expression-, SPP1-CD44 candidate axis-, and NicheNet-prioritized target program-related signatures were more prominent in advanced plaques and were positively correlated with each other. This integrative transcriptomic analysis identified a highly inflammatory macrophage state and a VSMC remodeling continuum in human atherosclerotic plaques. Candidate ligand-receptor and ligand-associated target gene expression patterns linked inflammatory macrophage activation with osteogenic/modulated VSMC remodeling at the computational level. External bulk data further showed coordinated enrichment of inflammatory and remodeling signatures in advanced plaques. These findings provide a descriptive and hypothesis-generating transcriptomic framework for understanding inflammatory macrophage activation and VSMC remodeling in human atherosclerosis.

atherosclerosis↗

ChemPerturb-seq screen identifies a small molecule cocktail enhancing human beta cell survival after subcutaneous transplantation.

Traditional chemical screens have focused on a single assay per screen, making them labor intensive and costly. Here, we combined a chemical screen with single-cell RNA sequencing (scRNA-seq) to perform Chemical Perturb-seq (ChemPerturb-seq), enabling a systematic analysis of the molecular changes of human beta cells upon individual small molecule treatments. Using this platform, we performed an in vivo barcoded screen and discovered a small molecule cocktail, including beta-lipotropin 61-91, insulin growth factor-1, and prostaglandin E2, with which preconditioning human beta cells and primary islets significantly enhanced function and survival when transplanted subcutaneously to female, but not to male, mice. We identified two additional molecules, serotonin and histamine, that promote islet function when transplanted subcutaneously to male mice using ChemPerturb-seq. Such small molecule cocktails could be applied to improve the current FDA-approved islet transplantation procedure. Finally, we developed an artificial intelligence (AI)-powered website, ChemPerturbDB, which provides user-friendly open access analysis of the extensive ChemPerturb-seq dataset.

Humans↗

Single-cell transcriptomics reveals distinct microglial state remodeling associated with the (R)-nicotine/diosmetin combination and galantamine in LPS-challenged BV2 cells.

BACKGROUND: Microglial neuroinflammation contributes to the progression of neurodegenerative diseases, yet it remains challenging to attenuate inflammatory responses while preserving cellular function. The effects of (R)-nicotine, diosmetin, their combined administration, and galantamine on heterogeneous BV2 transcriptional states have not been compared at single-cell resolution. METHODS: LPS-stimulated BV2 microglial-like cells were treated with (R)-nicotine, diosmetin, their combination (DR), or galantamine. Single-cell RNA sequencing was performed with three biological replicates per group and integrated with RNA velocity and SCENIC regulon inference to characterize treatment-associated state redistribution, inferred local transcriptional directionality and regulon-activity patterns. Functional validation included CCK-8 metabolic activity assays, multiplex cytokine ELISA, BDNF/GDNF quantification, qPCR, and high-content immunofluorescence analysis of iNOS and Arg1 at single-cell resolution. RESULTS: LPS decreased the relative abundance of the Itgae+/Plk4+ cluster while increasing the Nmur1+/Limk2+ cluster and inflammatory effector programs. DR treatment suppressed pro-inflammatory cytokine release without significantly reducing CCK-8-assessed metabolic activity, increased the Itgae+/Plk4+ cluster proportion, and increased BDNF/GDNF relative to LPS. RNA velocity and SCENIC analyses suggested that DR and galantamine showed distinct transcriptional and regulatory patterns: DR attenuated Batf-associated inflammatory regulons and was associated with increased Atf3-linked stress-response activity, whereas galantamine preferentially engaged DNA repair and genome-maintenance programs. CONCLUSION: These findings indicate that the DR condition was associated with remodeling of LPS-challenged BV2 microglial-like states, attenuation of inflammatory programs, and increased neurotrophic outputs relative to LPS, without significantly reducing CCK-8-assessed metabolic activity. This study provides a single-cell characterization of distinct treatment-associated responses to (R)-nicotine, diosmetin, their combined administration, and galantamine.

Microglia↗

Single-cell and spatial transcriptomic technologies for lung cancer tumor microenvironment analysis.

Lung cancer remains one of the leading causes of cancer-related mortality worldwide; beyond its rising incidence, its marked molecular heterogeneity and complex tumor microenvironment (TME) hinder treatment response and drive resistance, contributing directly to its high mortality rate. Single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) provide complementary approaches for dissecting these features. scRNA-seq enables high-resolution analysis of cellular diversity and transcriptional states but requires tissue dissociation and therefore loses spatial context. In contrast, ST preserves tissue architecture and provides insights into how gene-expression programs within the TME are organized, although no currently available spatial platform combines whole-transcriptome coverage with true single-cell resolution over large tissue areas. Together, these technologies have enabled detailed mapping of tumor, immune and stromal populations and of their spatial interactions, revealing functionally distinct cellular niches that contribute to immune evasion, metastasis and response to therapy. In this narrative review we organize the primary literature around a single question, how spatially structured cellular ecosystems, rather than individual cell types, determine therapeutic response and resistance in lung cancer - and we explicitly separate observations that are reproducible across independent cohorts and platforms from those that remain confined to single studies. We further summarize the technical, analytical and logistic barriers that currently prevent spatially resolved signatures from entering routine diagnostic pathology. Understanding dysregulated pathways and spatially constrained intercellular communication within the TME helps identify candidate biomarkers and may support the identification of therapeutic approaches directed at tumor-intrinsic programs as well as at microenvironment-driven resistance mechanisms.

Cell-cell communication↗

MitoScribe single-cell molecular recorder logs graded signaling dynamics into mitochondrial DNA.

Genetically encoded DNA recorders convert transient biological events into stable genomic mutations, offering a means to reconstruct past cellular states. However, current approaches to log historical events by modifying genomic DNA have limited capacity to record the magnitude of biological signals within individual cells. Here, we introduce MitoScribe, a mitochondrial DNA (mtDNA)-based recording platform that uses mtDNA base editors (DdCBEs) to write graded biological signals into mtDNA as neutral, single-nucleotide substitutions at a defined site. Taking advantage of the hundreds to thousands of mitochondrial genome copies per cell, we demonstrate MitoScribe enables reproducible, highly sensitive, non-destructive, durable, and high-throughput measurements of molecular signals, including hypoxia, NF-κB activity, BMP and Wnt signaling. We show multiple modes of operation, including multiplexed recordings of two independent signals, and coincidence detection of temporally overlapping signals. Coupling MitoScribe with single-cell RNA sequencing and mitochondrial transcript enrichment, we further reconstruct signaling dynamics at the single-cell transcriptome level. Applying this approach during the directed differentiation of human induced pluripotent stem cells (iPSCs) toward mesoderm, we show that early heterogeneity in response to a differentiation cue predicts the later cell state. Together, MitoScribe provides a scalable platform for high-resolution molecular recording in complex cellular contexts.

Journal Article↗

Atherosclerotic plaque fibroblasts derive from adventitial and medial Pdgfra-lineage-positive cells and predominantly maintain fibroblast identity.

AIMS: Fibroblasts are mesenchymal cells in the healthy vascular adventitia. In atherosclerosis, single-cell sequencing datasets suggest fibroblasts are abundant in plaques. However, their identity, origin, and fate during plaque progression remain unclear, which we aim to unravel here. APPROACH AND RESULTS: To robustly define fibroblast identity, origin, and fate, we employed meta-analyses of 54 single-cell RNA sequencing libraries, including murine smooth muscle cell (Myh11) and endothelial cell (EC) (Cdh5) lineage reporter mice with and without atherosclerosis; human control and atherosclerotic arteries; and murine adventitia and atherosclerotic plaques processed separately from low-density lipoprotein (LDL) receptor knockout (Ldlr-/-) mice. These meta-analyses showed that murine and human plaque fibroblast identity was robustly defined by Pdgfra, Pi16, Cygb, and Serpinf1 mRNA. Ninety-five percent of plaque fibroblasts do not derive from the Myh11 lineage, while no Cdh5-lineage-positive cells were present in the fibroblast cluster. We identified five murine arterial fibroblast subsets in atherosclerotic murine aorta: progenitor fibroblasts, matrix fibroblasts, inflammatory fibroblasts, an EC-like fibroblast subset, detected in both adventitia and plaques, and Col5a3+ fibroblasts, unique to the adventitia. We next studied fibroblast identity, origin, and fate using pseudotime analysis and Pdgfra-CreERT2/tdTomato lineage reporter mice (Pdgfra Lin+). Healthy Pdgfra Lin+ reporter mice showed predominant adventitial tdTomato expression, and infrequent medial and intimal Pdgfra Lin+ cells co-expressing MYH11 and PECAM1, respectively. The Pdgfra Lin+ plaque area increased with diet duration. Pdgfra Lin+ cells largely maintain fibroblast identity in the plaque, while <10% co-express SMC markers (MYH11, SM22&#x3b1;), or contribute to ACTA2+ cap cells. ECs gaining mesenchymal markers are transcriptionally distinct from Cdh5-lineage-negative fibroblasts gaining EC markers. Plaque-resident EC-like fibroblasts displayed a mesenchymal-to-endothelial transition transcriptome, which was induced in human primary fibroblasts in vitro by starvation, and dampened or reversed by IL1B, TGFB1, TGFB3, and oxidized LDL. Cross-species integration showed that all murine plaque fibroblasts were conserved in human atherosclerosis, with one additional subset partially resembling murine subsets, and three human-specific subsets. Importantly, human fibroblast subsets differentially correlated to human plaque traits, with EC-like fibroblasts correlating to plaque instability. CONCLUSION: Our results indicate that 95% of plaque-residing fibroblasts are Myh11 Lin- Plaque fibroblasts have a dual origin, predominantly adventitial Pdgfra Lin+ progenitor fibroblasts, with a minor contribution from medial Pdgfra Lin+ &#xa0;Myh11+ SMCs. Most plaque fibroblasts maintain fibroblast identity. Murine plaque fibroblast subsets were conserved in human atherosclerosis. EC-like fibroblasts are linked to human plaque instability. Intervening in progenitor-to-specific fibroblast transitions could present a new avenue to promote plaque stability in atherosclerosis.

Atherosclerosis↗

A Standardized Protocol for Generating iPSC-Derived Human Microglia for Functional Genomic Assays.

Human induced pluripotent stem cell (iPSC)-derived microglia (iMG) provide an in vitro experimental system for studying human microglial biology, neuroinflammation, and genetic risk mechanisms associated with neurological disease. This chapter describes a standardized, scalable, and reproducible protocol for the differentiation of human iPSCs into functional microglia-like cells, with particular emphasis on applications in transcriptional and epigenomic network analysis. The protocol supports high-viability floating iMG production, compatibility with pooled CRISPR perturbation approaches, and downstream multiomic profiling, including single-cell RNA sequencing, chromatin accessibility assays, and proteomics. Detailed procedures are provided for iPSC maintenance, hematopoietic progenitor cell generation, microglial maturation, functional genomics integration, and quality control.

Humans↗

An isoform-resolution transcriptomic atlas of colorectal cancer from long-read single-cell sequencing.

Colorectal cancer (CRC) ranks as the second leading cause of cancer deaths globally. In recent years, short-read single-cell RNA sequencing (scRNA-seq) has been instrumental in deciphering tumor heterogeneities. However, these studies only enable gene-level quantification but neglect alterations in transcript structures arising from alternative end processing or splicing. In this study, we integrated short- and long-read scRNA-seq of CRC samples to build an isoform-resolution CRC transcriptomic atlas. We identified 394 dysregulated transcript structures in tumor epithelial cells, including 299 resulting from various combinations of splicing events. Second, we characterized genes and isoforms associated with epithelial lineages and subpopulations exhibiting distinct prognoses. Among 31,935 isoforms with novel junctions, 330 were supported by The Cancer Genome Atlas RNA-seq and mass spectrometry data. Finally, we built an algorithm that integrated novel peptides derived from open reading frames of recurrent tumor-specific transcripts with mass spectrometry data and identified recurring neoepitopes that may aid the development of cancer vaccines.

Colorectal Neoplasms↗

Multi-level Transcriptomic and Machine-learning Analyses Identify MZT1 as a Proliferation-associated Prognostic Marker in Lung Adenocarcinoma.

BACKGROUND/AIM: Lung adenocarcinoma (LUAD) exhibits substantial molecular heterogeneity and variable clinical outcomes, highlighting the need for biomarkers that reflect core tumor biological processes. Centrosome-associated proteins regulate mitotic fidelity and genome stability, yet their roles in LUAD remain incompletely defined. In this study, we systematically characterized mitotic spindle organizing protein 1 (MOZART1; MZT1) and related family members in LUAD. MATERIALS AND METHODS: We performed integrated analyses combining bulk transcriptomic datasets, survival modeling, gene set enrichment, immune deconvolution, machine-learning based prognostic modeling, and single-cell RNA sequencing. Expression patterns and clinical associations of MZT family genes were evaluated across pan-cancer and LUAD cohorts. RESULTS: MZT family genes were consistently upregulated in tumor tissues, with MZT1 showing the most robust expression pattern. Elevated MZT1 expression was significantly associated with reduced overall survival. Functional analyses revealed coordinated activation of proliferative and genome maintenance pathways, including G2/M checkpoint regulation, E2F and MYC signaling, and DNA repair. A multivariable analysis indicated that the prognostic association of MZT1 was reduced after adjusting for canonical proliferation markers, suggesting partial overlap with established proliferation signals. The LASSO-based Cox model demonstrated stable time-dependent predictive performance at 1-, 3-, and 5-year survival. Immune analyses indicated associations between MZT1 expression and tumor microenvironmental features. Single-cell analysis showed that MZT1 expression was predominantly enriched in malignant epithelial cells and associated with proliferative cellular states. Protein-level validation supported concordance with transcriptomic findings. CONCLUSION: MZT1 is a proliferation-associated marker that integrates clinical risk, transcriptional programs, cellular heterogeneity, and predictive modeling in LUAD, providing a potential framework for biomarker development and risk stratification.

Humans↗

Uncovering heterogeneous effects via localized feature selection.

Identifying features that interact to trigger disease, while accounting for heterogeneity across diverse populations, is essential for the development of precision and targeted medicine. Despite the availability of vast and complex health-related datasets, most existing works focus on identifying disease-associated features at the population level or within a few subpopulations, often overlooking individual-level heterogeneity within these groups. To address this limitation, we propose a framework that utilizes localized test statistics to identify disease-associated features tailored to individual profiles. Our method leverages the recently developed knockoffs methodology to control the noise level of the selection set so that the results are replicable. Moreover, it allows for the discovery of hidden heterogeneous effects within the data, as demonstrated in an application to single-cell RNA sequencing data for Alzheimer's disease. By aggregating localized feature selection results, our framework also enables powerful population-level feature selection. Our framework provides a powerful tool for exploratory studies of precision medicine, offering the potential to generate novel hypotheses for confirmatory biological experiments.

Alzheimer Disease↗

PreDigs: A Database of Context-specific Cell Type Markers and Precise Cell Subtypes for Digestive Cell Annotation.

Research on cell type markers helps investigators explore the diverse cellular composition of gastrointestinal tumors, thereby enhancing our understanding of tumor heterogeneity and its impact on disease progression and treatment response. However,&#xa0;the integration of large-scale datasets and the standardization of cell type identification remain challenging. Here, we developed PreDigs, a user-friendly database of predicted signatures for the digestive system, which offers 124 curated single-cell RNA sequencing datasets, covering over 3.4 million cells, all available for download. After unsupervised clustering, we unified the identification and nomenclature of cell subtype labels, constructing a cell ontology tree with 142 cell types across 8 hierarchical levels. Meanwhile, we calculated three different context-specific cell type markers, including "Cell Markers", "Subtype Markers", and "TPN Markers", based on various application requirements within or across tissues. Through the integrated analysis of PreDigs data, we identified distinct cell subpopulations exclusive to tumors, one of which corresponds to tumor-specific endothelial cells. Additionally, PreDigs offers online cell annotation tools, allowing users to classify single cells with greater flexibility. PreDigs is accessible at https://www.biosino.org/predigs/.

Humans↗

scRNA-seq and bulk RNA-seq reveal the characteristics of macrophage copper metabolism and establish a risk signature in hepatocellular carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) is a prevalent malignancy with an urgent need for improved prognostic stratification and treatment-response prediction. This study aimed to explore a macrophage copper metabolism-associated prognostic model and to investigate the relationship between this risk model and the tumor immune microenvironment. METHODS: The FindClusters function was used to analyze cell clusters, and CellChat and CellPhoneDB/LIANA were employed for cell-cell communication analysis. Copper metabolism-related genes were sourced from the MSigDB database. A prognostic risk model was established using least absolute shrinkage and selection operator (LASSO) analysis and multivariate Cox regression analysis, and a nomogram was constructed by integrating the prognostic model with clinicopathological factors. Additional analyses were performed to map the seven model genes in single-cell data, assess model uncertainty and robustness, evaluate macrophage/copper/cuproptosis-related transcriptional programs, and examine the correlations between risk score, immune infiltration and predicted drug sensitivity. RESULTS: Using single-cell RNA sequencing (scRNA-seq) data, we identified four macrophage subpopulations. Macrophages with high SPP1 expression showed close interaction with T cell populations and were associated with copper ion metabolism. By incorporating 141 copper metabolism-related genes and using The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) cohort, we constructed a seven-gene risk prediction model. Additional single-cell mapping showed that the model genes were detectable in the HCC single-cell dataset and showed a macrophage-associated expression pattern. The model showed moderate prognostic discrimination in TCGA-LIHC, whereas its external performance was heterogeneous and remained evaluable across external cohorts, with performance varying among datasets. Immune and mechanism-related analyses suggested that the risk signature was associated with macrophage-related infiltration, copper metabolism and cuproptosis-related transcriptional programs. Drug sensitivity analysis nominated Daporinad as a computationally predicted candidate compound, supporting Daporinad as a pharmacogenomic candidate for follow-up investigation. CONCLUSIONS: By integrating scRNA-seq and bulk RNA sequencing (RNA-seq) data, we constructed a macrophage copper metabolism-associated prognostic signature for HCC. The risk score was associated with survival, immune microenvironment features and predicted drug response, providing a transcriptomic framework for risk stratification and therapeutic hypothesis generation.

Hepatocellular carcinoma (HCC)↗

Multimodal Analysis Reveals Aberrant Expression of SUMO2 and Its Significant Association With Key Mechanisms of Metabolic Pathways in Hepatocellular Carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related deaths worldwide. However, the role of small ubiquitin-like modifier 2 (SUMO2), a core member of the small ubiquitin-like modifier (SUMO) family, regarding its expression patterns and metabolism-related functions in HCC remains inadequately understood. METHODS: A multidimensional analytical framework was applied, integrating immunohistochemistry (153 HCC vs. 21 non-HCC samples), proteomics (159 paired samples), bulk transcriptomics (3240 HCC vs. 2267 non-HCC samples), single-cell RNA sequencing (RNA-seq) (10 HCC vs. 8 non-HCC samples), spatial transcriptomics, and external CRISPR/Cas9 functional genomics data. Systematic analyses included standardized mean difference (SMD), pathway enrichment, pseudotime trajectory inference, in silico knockout, cell-cell communication, metabolic flux scoring, immune infiltration, clinical correlation, drug sensitivity prediction, and molecular docking. RESULTS: At the protein level, immunohistochemistry (nuclear positivity) and external proteomic data collectively demonstrated consistent SUMO2 overexpression in HCC. Consistent upregulation was also observed at the mRNA level across large-scale cohorts. Single-cell RNA-seq and spatial transcriptomics localized SUMO2 enrichment to malignant hepatocytes and tumor-dominant regions. CRISPR-mediated SUMO2 knockout suppressed proliferation in multiple HCC cell lines. Mechanistically, high SUMO2 expression was significantly associated with metabolic reprogramming involving glycolysis/gluconeogenesis, pyruvate metabolism, and the tricarboxylic acid cycle. SUMO2-high malignant hepatocyte subpopulations exhibited enhanced activity of the macrophage migration inhibitory factor signaling axis and enhanced iron-sensor interactions. Further, the immune infiltration analysis revealed a negative correlation between SUMO2 expression and M1 macrophages and a positive correlation between follicular helper T cells and regulatory T cells. Clinically, elevated SUMO2 levels were found to be associated with adverse prognostic features. Furthermore, high SUMO2 expression was associated with increased sensitivity to dasatinib, and molecular docking simulations predicted potential binding between SUMO2 and dasatinib, with a Vina score of -8.5 kcal/mol. CONCLUSIONS: SUMO2 is aberrantly expressed at the protein, mRNA, single-cell, and spatial transcriptomic levels in HCC and is significantly associated with metabolic reprogramming and altered migration inhibitory factor (MIF)-mediated intercellular communication, suggesting its potential as a novel biomarker for diagnosis and treatment.

Humans↗

Amaranth: enhanced single-cell transcript assembly via discriminative modelling of UMI reads and internal reads.

MOTIVATION: Single-cell RNA sequencing (scRNA-seq) has transformed transcriptome profiling at cellular resolution, yet accurate reconstruction of full-length transcripts for individual cells remains a central challenge. Emerging scRNA-seq protocols can produce reads that span entire transcripts, enabling isoform-level expression analysis. For example, Smart-seq protocols combine unique molecular identifier (UMI)-linked reads that index and stitch together multiple reads from the same molecule, with internal reads filling coverage gaps. We demonstrate that these read types exhibit markedly different biological and statistical properties in strandness, 5'/3' coverage bias, and genomic locality. Existing assemblers fail to leverage these distinctions, yielding suboptimal assembly. RESULTS: We developed Amaranth, a novel single-cell assembler that discriminatively models UMI and internal reads. Amaranth implements heuristics specifically designed to address the distinct biases of UMI-linked and internal reads, enabling accurate strandness assignment for internal reads, reliable splicing graph refinement, and precise transcript start site determination. We also developed Amaranth-meta, which integrates information across cells to enhance individual cell assemblies. Benchmarked on Smart-seq3 datasets from human HEK293T and mouse fibroblast cells, Amaranth outperformed other state-of-the-art assemblers in assembling individual cells and in meta-assembly. Amaranth advances isoform-level analysis in single-cell transcriptomics, facilitating detailed studies at cellular resolution. AVAILABILITY AND IMPLEMENTATION: Amaranth is implemented in C++ and is freely available at https://github.com/Shao-Group/amaranth under the BSD-3-Clause license. Scripts, documentation, and data for reproducing experiments in this manuscript are available at https://github.com/Shao-Group/amaranth-test.

Single-Cell Gene Expression Analysis↗

Colorectal Liver Metastasis Pathomics Model: Integrating Single-Cell and Spatial Transcriptome Analysis With Pathomics for Predicting Liver Metastasis in Colorectal Cancer.

The liver is the primary target organ for hematologic metastasis of colorectal cancer (CRC), and CRC liver metastasis (CRLM) often precludes radical resection, making it the leading cause of death in patients with CRC. To improve the identification and prediction of liver metastasis risk, we identified a cell type of liver metastasis--triggering malignant cells (LMTMCs) through integrating single-cell RNA sequencing and spatial transcriptome analysis. Multiomics cell communication analysis indicated that the interaction between fibroblasts and LMTMCs through the COL1A1-CD44/SDC4 and LAMA4-CD44 signaling axes could promote CRLM. By applying the one-class logistic regression algorithm, we developed a CRLM scoring system in the bulk RNA-sequencing data according to the abundance of LMTMCs in each individual. Using the grouping labels derived from the CRLM scoring system in the bulk data and the corresponding whole-slide images without any manual annotations at the region or pixel level, processed via slide-level weakly supervised learning, a deep-learning model based on the ResNet18 architecture, called Colorectal Liver Metastasis Pathomics Model, was developed to predict the risk of liver metastasis in patients with CRC. The Colorectal Liver Metastasis Pathomics Model achieved an area under the curve of 0.84 at the internal test set of The Cancer Genome Atlas-CRC histology images. In the external independent validation sets, namely the Affiliated Hospital of Southwest Medical University and the Affiliated Traditional Chinese Medicine Hospital of Southwest Medical University cohorts, the areas under the curve were 0.89 and 0.72, respectively, indicating effective classification performances. This study provided new insights and tools for the early identification of CRLM and demonstrated the potential of combining multiomics with deep learning-based pathomics in cancer research.

Humans↗

Single-cell transcriptomics reveals heterogeneous stress responses and Mg2+-mediated survival mechanisms in Lactobacillus delbrueckii subsp. bulgaricus during freeze-drying and storage.

Maintaining the viability of lactic acid bacteria during dehydration and subsequent storage remains a significant challenge. Here, we employed single-cell RNA sequencing to reveal the heterogeneous stress responses of Lactobacillus delbrueckii subsp. bulgaricus, identifying seven distinct transcriptional clusters across the liquid culture, freeze-drying, and storage phases. The dominant clusters in the freeze-drying and storage were not completely consistent, showing significant functional differentiation. Genomic stability may be important for survival during freeze-drying and storage, while intracellular energy homeostasis appears important for viability during storage. The magnesium transporter mgtB was highly expressed in clusters tolerant to freeze-drying and storage, suggesting a critical role for Mg2+ homeostasis. Further experimental validation confirmed that Mg2+ treatment significantly bolstered stress resistance, increasing immediate post-freeze-drying survival by over 2-fold (up to 92.90%) and post-storage survival by over 5-fold (up to 5.98%). Proteomic data indicated that Mg2+ supplementation correlated with the maintenance of several biological functions potentially relevant to bacterial survival during freeze-drying and storage, including DNA repair, translation, and central carbon metabolism. These findings provide a map of microbial stress resistance through population heterogeneity and offer a potential strategy that may be adapted for enhancing the stability of other industrial lactic acid bacteria products.

Freeze Drying↗