Ionic current mechanisms generating vertebrate primary cardiac pacemaker activity at the single cell level: an integrative view.
Explore the source record for details and available documents.
Amphibians↗
SEARCH · Search PubMed
Results for “single cell integration”
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.
Stage-specific ROMO1 in rheumatoid arthritis: predictive immune insights into the MIF pathway and HLA-DR/IL2RA axis via integrated GWAS, transcriptomic, single-cell, and spatial profiling.
Emerging evidence links reactive oxygen species modulator 1 (ROMO1), a key mitochondrial ROS regulator, to rheumatoid arthritis (RA) pathogenesis. However, its exact mechanism remains elusive given the conflicting evidence about its specific function. We used a four-level integrative framework combining multi-omics data and literature‑supported mechanistic inference. At the genetic level, Mendelian randomization (MR) was performed to explore potential causal relationships between ROMO1, IL2RA, HLA-DR, MIF, and RA risk, followed by differential expression analysis and machine learning-based feature selection to identify key mROS genes. The temporal expression dynamics of ROMO1 were assessed in RA progression. At the cellular and tissue levels, we integrated single-cell RNA sequencing and spatial transcriptomics to map cell-type-specific expression and synovial localization of ROMO1-related immune cells and pathways. Finally, our multi-omics findings were contextualized with literature-supported mechanistic inference. (1) MR results were consistent with a potential protective effect of ROMO1 on RA (OR = 0.52) and its potential regulation of risk factors IL2RA (OR = 0.46) and HLA-DR (OR = 0.40). Conversely, IL2RA (OR = 1.42), HLA-DR (OR = 1.88), and MIF (OR = 1.17) were positively associated with RA risk. Additionally, ROMO1 was identified as a top candidate diagnostic predictor with stage-specific dynamics: downregulated in the early but upregulated in the late/remission stages. (2) Single-cell RNA sequencing showed ROMO1's cell-specific expression in CD14+ HLA-DR+ CD74+ monocytes and CD4+ IL2RA+ T cells. Cell communication analysis further suggested that these cells may participate in MIF pathway regulation. Spatial transcriptomics subsequently identified that ROMO1-related cells localized to synovial pathological regions, with MIF pathway changes correlated with RA progression. (3) Finally, literature-supported mechanistic inference suggests that ROMO1 may modulate mROS levels to promote anti-inflammatory M2 macrophage polarization, which could theoretically contribute to reduced systemic inflammation and the alleviation of multi-organ decline in RA. This integrated multi-omics investigation, supported by literature-based mechanistic inference, suggests ROMO1 as a stage-dependent biomarker candidate and potential immune regulator in RA.
Humans↗
Integrated method for single-cell DNA extraction, PCR amplification, and sequencing of ribosomal DNA from harmful dinoflagellates Cochlodinium polykrikoides and Alexandrium catenella.
A simplified technique was developed for DNA sequence-based diagnosis of harmful dinoflagellate species. This protocol integrates procedures for DNA extraction and polymerase chain reaction (PCR) amplification into a single tube. DNA sequencing reactions were performed directly, using unpurified PCR products as the DNA template for subsequent sequencing reactions. PCR reactions using DNA extracted from single cells of Cocodinium polykrikoides and Alexandrium catenella successfully amplified the target ribosomal DNA regions. DNA sequencing of the unpurified PCR products showed that DNA sequences corresponded to the expected locus of ribosomal DNA regions of both A. catenella and C. polykrikoides (each zero genetic distance and 100% sequence similarity). Using the protocol described in this article, there was little DNA loss during the purification step, and the technique was found to be rapid and inexpensive. This protocol clearly resolves the taxonomic ambiguities of closely related algal species (such as Alexandrium and Cochlodinium), and it constitutes a significant breakthrough for the molecular analysis of nonculturable dinoflagellate species.
Animals↗
Basal ganglia precursors found in aggregates following embryonic transplantation adopt a striatal phenotype in heterotopic locations.
Transplantation of immature CNS-derived cells into the developing brain is a powerful approach to investigate the factors that regulate neuronal position and phenotype. CNS progenitor cells dissociated from the embryonic striatum and implanted into the brain of embryos of the same species generate cells that reaggregate to form easily recognizable structures that we previously called clusters and cells that disperse and integrate as single cells into the host brain. We sought to determine if the neurons in the clusters differentiate according to their final location or acquire a striatal phenotype in heterotopic positions. We transplanted dissociated cells from the E14 rat medial and lateral ganglionic eminences, either combined or in isolation, into the E16 embryonic rat brain. At all time points, we found clusters of BrdU- and DiI-labelled donor cells located in the forebrain and hindbrain, without any apparent preference for striatum. Immunocytochemical analyses revealed that cells in the clusters expressed DARPP-32 and ARPP-21, two antigens typically co-expressed in striatal medium-sized spiny neurons. In agreement with observations previously noted by several groups, isolated cells integrated into heterologous host areas do not express basal ganglia phenotypes. These data imply that immature striatal neuronal progenitors exert a community effect on each other that is permissive and/or instructive for development of a striatal phenotype in heterotopic locations.
Animals↗
Immune-Like Malignant Epithelial Programs Shape Tumor-Immune Interactions and Inform Prognostic Stratification in Lung Adenocarcinoma.
Lung adenocarcinoma (LUAD) is characterized by marked cellular heterogeneity, yet how malignant epithelial states contribute to immune regulation and clinical outcomes remains incompletely defined. We integrated single-cell RNA-sequencing data to map the cellular landscape of LUAD and identify malignant epithelial cells based on inferred copy-number alterations. Epithelial states were further examined through trajectory inference, transcription factor analysis, and cell-cell communication profiling. Single-cell-derived genes were subsequently integrated with TCGA and independent GEO cohorts to construct and validate a machine learning-based prognostic signature. Malignant epithelial cells displayed distinct functional programs, including an immune-like state associated with genomic instability, immune-related transcriptional activity, tumor-immune communication, and patient outcomes. The resulting immune-like malignant epithelial cell signature (IMEC-Sig) consistently stratified survival across multiple cohorts. Low IMEC-Sig scores were accompanied by greater immune infiltration, higher immune checkpoint expression, and increased immunophenoscore, whereas high scores were linked to a comparatively immunosuppressive phenotype. Pan-cancer analyses further identified KRT8 as a gene associated with unfavorable prognosis, and functional experiments showed that KRT8 silencing suppressed proliferation, migration, invasion, and colony formation in LUAD cells. Together, these findings connect malignant epithelial heterogeneity with the immune context and clinical outcomes, support IMEC-Sig as a biologically informed prognostic tool, and nominate KRT8 as a potential therapeutic target in LUAD.
Humans↗
Network methods for diagonal integration of unpaired single-cell multiomics data: a review.
MOTIVATION: Advances in single-cell sequencing have enabled multiomics profiling at unprecedented resolution; however, mass spectrometry-based single-cell proteomics (scMS) remains inherently destructive, precluding simultaneous transcriptomic capture. Unlike antibody-based methods such as CITE-seq, which permit paired profiling but are restricted to targeted protein panels, scMS provides unbiased, genome-scale coverage of the intracellular proteome yet necessitates post hoc integration of unpaired datasets. This diagonal integration challenge, where transcriptomes and proteomes are measured in separate cells lacking shared anchors, remains underserved by existing reviews, which focus predominantly on vertical integration strategies enabled by non-destructive assays. RESULTS: We survey the complete computational pipeline for constructing mechanistic proteogenomic networks from unpaired single-cell data, covering: (i) unimodal network inference such as knowledge-based approaches, probabilistic graphical models, temporal directionality inference, and generative and foundation model strategies that establish the transcriptomic scaffold; (ii) cross-modal integration architectures such as network propagation, graph neural networks (scMRDR, scmFormer, scCotag), and consensus frameworks designed explicitly for the unpaired proteomics setting; and (iii) benchmarking paradigms spanning network reconstruction (BEELINE, GRETA, CausalBench) and multi-task integration evaluation (scMultiBench, SCMMIB), with guidance on metric selection under network sparsity and class imbalance. We identify three principal axes of future development: generative proteomic translation from transcriptomic precursors, inductive prior embedding in next-generation architectures, and perturbation-based causal benchmarking. AVAILABILITY AND IMPLEMENTATION: This is a review article; no novel software is distributed. A curated benchmark resource table, methods starter guide, and per-method bottleneck annotations are provided in the Supplementary Material.
Multiomics↗
Integrated Bulk and Single-Cell RNA-Seq Analysis Reveals Transcriptional Activation of PTGS2 by FOS in Progression From T2DM to T2DM-Associated NAFLD.
Type 2 diabetes mellitus (T2DM) and nonalcoholic fatty liver disease (NAFLD) frequently coexist, exacerbating disease burden. However, the molecular mechanisms underlying the progression from T2DM to T2DM-associated NAFLD remain unclear. This study investigated the regulatory function of FOS-mediated PTGS2 activation in this transition. We integrated bulk RNA-seq data from GEO, single-cell transcriptomic data and transcriptomes from patients with T2DM-associated NAFLD. Differentially expressed genes were identified using the limma package, and T2DM-related gene modules were defined by weighted gene co-expression network analysis. LASSO regression and random forest identified 14 candidate genes, with PTGS2 and FOS prioritised. Single-cell analysis showed increased FOS and PTGS2 expression in monocytes, CD8+ T cells and Kupffer cells. Transcription factor prediction and dual-luciferase assays confirmed that FOS directly binds the PTGS2 promoter and drives its transcription. In vitro, FOS silencing decreased PTGS2 expression, cytokine secretion and apoptosis under high-glucose and free fatty acid conditions, whereas PTGS2 overexpression exacerbated inflammation and apoptosis independently of FOS expression. These findings demonstrate that FOS transcriptionally activates PTGS2, contributing to hepatic inflammation and apoptosis during the progression from T2DM to NAFLD. PTGS2 may serve as a promising biomarker and therapeutic target for T2DM-associated NAFLD.
Single-Cell Gene Expression Analysis↗
Single-cell multiomics reveals exosome-mediated reprogramming and clonotypic remodeling of T cells in triple-negative breast cancer.
Triple-negative breast cancer (TNBC) is an aggressive and immunogenic subtype lacking targeted therapies. While tumor-derived exosomes are known to modulate immune function, their direct impact on human T cell plasticity and antigen specificity remains poorly defined. Here, we conducted a comprehensive single-cell multiomic analysis of primary human T cells exposed to exosomes derived from 17 genomically diverse TNBC cell lines and 35 patient samples. Integrating single-cell RNA-seq, V(D)J sequencing, non-coding RNA profiling, bulk and single-cell cytokine analyses, we uncovered conserved and subtype-specific immunomodulatory programs induced by TNBC exosomes. Exosome-treated T cells displayed skewing toward regulatory and dysfunctional phenotypes, including Th17-like, Treg, and PD-1⁺/PD-L1⁺ Tfh cells. Functional profiling revealed suppression of early activation markers and cytokine responses, alongside selective preservation of cytotoxic features in γδ T and NKT subsets. Transcriptomic and miRNA network analyses demonstrated widespread downregulation of immune effector genes (e.g., HBEGF and TNFSF9) mediated by exosome-delivered regulatory miRNAs (has-miR-98-5p). Notably, exosome-stimulated T cells displayed distinct clonotypic expansions, characterized by the emergence of five tumor-specific γδ TCR clonotypes and 30 unique αβ TCR CDR3 sequences that were absent in mock-treated controls, underscoring the role of exosomes in shaping TCR repertoire dynamics.
Humans↗
transFusion: a novel comprehensive platform for integration analysis of single-cell and spatial transcriptomics.
MOTIVATION: Understanding spatial organization, intercellular interactions, and regulatory networks within the spatial context of tissues is crucial for uncovering complex biological processes and disease mechanisms. Spatial transcriptomics technologies have revolutionized this field by enabling the spatially resolved profiling of gene expression. 10× Visium has emerged as the predominant spatial technology, but its low resolution and the complexity of integrating multimodal datasets present significant analytical challenges, particularly for researchers with limited computational and statistical expertise. Current spatial transcriptomics analysis platforms generally fall short of effectively integrating multimodal data and maximizing the utility of spatial information-such as uncovering complex cellular spatial dependencies, multimodal gradient patterns, and spatial coexpression of ligand-receptor pairs and regulatory networks related to disease or biological states-thereby limiting their ability to provide comprehensive end-to-end analytical workflows when analyzing 10× Visium data. RESULTS: To address these limitations, we developed transFusion, a novel, advanced web-based platform specializing in the most comprehensive and effective integration analysis of scRNA-seq and 10× Visium spatial transcriptomics data. transFusion offers 12 key functions, from basic visualization to advanced analyses, including intercellular dependency analysis, ligand-receptor coexpression identification and visualization, and spatial multimodal gradient variation patterns. Two case studies were used to demonstrate transFusion's capabilities in exploring tissue architecture, intercellular communication, dependency networks, and multimodal gradient variation patterns with minimal computational skills and statistical expertise. transFusion provides a flexible and powerful framework for multimodal data integration analysis. AVAILABILITY AND IMPLEMENTATION: transFusion is freely available at https://github.com/WQLin8/transFusion.
Spatial Transcriptomics↗
Integrated Pan-Cancer, Single-Cell, and Spatial Transcriptomic Analyses Identify ZDHHC12 as a Biomarker Associated with Macrophage Infiltration and the Immune Landscape in Glioma.
BACKGROUND: The tumor immune microenvironment (TME) critically influences cancer progression and therapeutic response. However, the pan-cancer expression landscape, prognostic relevance, and spatial distribution of ZDHHC12 remain incompletely characterized. This study investigated the prognostic value of ZDHHC12 and its associations with immune microenvironmental features and drug sensitivity. METHODS: Data from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) datasets were used to evaluate ZDHHC12 expression and prognosis across cancer types. Immune infiltration analyses, single-cell RNA sequencing, and spatial transcriptomics were integrated to characterize the associations of ZDHHC12 with the cancer immunity cycle and the spatial architecture of glioma. Drug sensitivity and immunotherapy-related metrics were assessed using pharmacogenomic databases and computational prediction models. RESULTS: ZDHHC12 was aberrantly expressed across multiple tumors and was associated with patient prognosis. Its expression was broadly correlated with immune cell recruitment- and activation-related signatures. In glioma, single-cell and spatial transcriptomic analyses showed enrichment of ZDHHC12 in monocyte/macrophage populations and spatial co-localization with BAK1, CD68, and CD163. ZDHHC12 expression was also associated with predicted drug sensitivity and immunotherapy-related metrics. CONCLUSION: ZDHHC12 may serve as a candidate pan-cancer prognostic biomarker. In glioma, its expression is associated with macrophage-enriched and immunosuppressive microenvironmental features. Functional studies are required to establish causality and determine its therapeutic relevance.
GBM↗
A pan-cancer single-cell atlas uncovers the role of sex hormones and chromosomes in sex-divergent reprogramming of the tumor microenvironment.
BACKGROUND: Sex bias is pervasive in tumors; however, how sex chromosomes and hormone-responsive signaling shape the tumor microenvironment (TME) remains insufficiently characterized. Considering the critical impact of the TME on tumor progression and response to immunotherapy, a pan-cancer investigation of sex-specific and cancer-context-dependent TME features is warranted. METHOD: Based on stringent inclusion criteria, we constructed a high-resolution pan-cancer single-cell sequencing atlas by integrating 31 publicly available single-cell RNA-seq datasets, comprising a total of 1,831,436 cells by integrating 468 samples from eight types of non-sex-specific solid tumors (282 males and 186 females). After correcting for batch effects, we identified major and minor cellular subsets. Multiple computational approaches were applied to investigate sex-associated differences in cellular composition, gene expression, pathway activity, malignant cell states and intercellular communication. RESULTS: We systematically compared sex-specific TME features across eight common solid malignancies. Male-biased CD8+ T cell exhaustion emerged as a recurrent but non-uniform feature, with its magnitude varying across cancer types and being modified by tissue-specific contexts. This pattern was associated with androgen-response signature scores and expression-based loss of the Y chromosome (LOY) scores. M2-like macrophage polarization showed a more cancer-type-dependent pattern; although female-biased enrichment was observed in selected malignancies, it did not represent a uniform pan-cancer feature. Expression-based X chromosome inactivation (XCI)/XCI escape-related programs, estrogen-response signature scores and stromal components, including fibroblasts and endothelial cells, were associated with macrophage and immune-regulatory states in specific tumor contexts. Tumor cells of male origin displayed higher genomic instability and more aggressive phenotypes, with androgen-response signatures and LOY contributing to the development of a male biased malignant state. Furthermore, expression-based LOY scores in malignant cells were associated with CD8+ T cell exhaustion based on transcriptomic proxies. CONCLUSION: Our study uncovers extensive but heterogeneous sex-specific differences in the TME across multiple cancer types. We propose a regulatory framework linking sex chromosomes, hormone-responsive signaling and TME interactions, which is consistent with recurrent male-biased CD8⁺ T cell exhaustion and context-dependent M2-like macrophage polarization. Importantly, the magnitude and, in some cancers, the direction of these sex-biased features are modified by tissue-specific contexts. These findings underscore the need to include sex chromosome and hormone status as essential biological variables in studies of the tumor microenvironment and the design of immunotherapies.
Tumor Microenvironment↗
Large-scale neural model for visual attention: integration of experimental single-cell and fMRI data.
A computational neuroscience framework is proposed to better understand the role and the neuronal correlate of spatial attention modulation in visual perception. The model consists of several interconnected modules that can be related to the different areas of the dorsal and ventral paths of the visual cortex. Competitive neural interactions are implemented at both microscopic and interareal levels, according to the biased competition hypothesis. This hypothesis has been experimentally confirmed in studies in humans using functional magnetic resonance imaging (fMRI) techniques and also in single-cell recording studies in monkeys. Within this neuro-dynamical approach, numerical simulations are carried out that describe both the fMRI and the electrophysiological data. The proposed model draws together data of different spatial and temporal resolution, as are the above-mentioned imaging and single-cell results.
Algorithms↗
SSB deficiency-induced R-loop accumulation triggers podocyte inflammation in DKD.
INTRODUCTION: Diabetic kidney disease (DKD) is fundamentally a podocytopathy in which sterile inflammation plays a central pathogenic role, yet the upstream triggers that initiate inflammatory cascades in podocytes remain elusive. R-loops are critical regulators of genomic stability, and their pathological accumulation triggers DNA damage and innate immune activation. Whether R-loop dysregulation contributes to podocyte-driven inflammation in DKD is unknown. METHODS: We integrated single-cell transcriptomic profiling, dual machine learning algorithms, and functional experiments to dissect the R-loop regulatory network in the diabetic kidney. RESULTS: Integrated analysis of human diabetic kidney single-cell RNA-seq data revealed a globally compromised R-loop regulatory network selectively within podocytes. Intersection of podocyte-specific transcriptomic shifts with validated R-loop regulators identified 93 candidate genes, from which dual machine learning algorithms pinpointed SSB (Sjögren syndrome antigen B) as the principal podocyte-selective R-loop resolver and a superior diagnostic biomarker (AUC = 0.983). SSB expression was selectively downregulated in diabetic podocytes and showed the strongest positive correlation with the R-loop resolution module. Mechanistically, SSB loss impaired RNA splicing and stability pathways, leading to aberrant R-loop accumulation that activated the cGAS-dependent inflammatory signaling in podocytes. In two murine DKD models and high glucose-challenged podocytes, SSB was markedly reduced. Remarkably, SSB knockdown in podocytes alone sufficed to trigger R-loop accumulation and pro-inflammatory cytokine expression, whereas both RNase H1-mediated R-loop removal and cGAS co-depletion blunted this response. DISCUSSION: These findings suggest that an SSB-governed R-loop -cGAS -inflammatory signaling axis may link genomic instability to podocyte inflammation and contribute to DKD progression, nominating R-loop homeostasis as a previously unrecognized potential therapeutic target.
Podocytes↗
Single-cell transcriptome revealed the aberrant keratinocytes activation in antigen presentation in atopic dermatitis.
BACKGROUND: Atopic dermatitis (AD), a common chronic inflammatory skin disease, has been extensively studied using single-cell genomics. However, keratinocytes, as key effector cells in AD, have underlying mechanisms remain incompletely understood and require further investigation. METHODS: We integrated single-cell transcriptomic data from skin tissues of healthy controls, chronic active AD patients, spontaneously healed AD (SHAD) patients, and an ovalbumin-induced AD mouse model. The study particularly emphasized the gene expression and cellular dynamics of keratinocytes across the different groups, as well as their interactions with immune cells. RESULTS: Compared to healthy controls, we observed significant changes in the keratinocyte transcriptome, cellular state, and keratinocyte-immune cell ligand-receptor interactions in AD skin, particularly the marked activation of genes involved in antigen processing and presentation. Interestingly, such gene activation was not observed in keratinocytes from the ovalbumin-induced AD mouse model, despite its phenotype closely resembling human AD. Furthermore, in SHAD, we identified a recovery of both the ligand-receptor interaction patterns and antigen processing and presentation genes, accompanied by a notable shift in the transcriptome. This involved a significant downregulation of genes related to cytoplasmic transcription and oxidative phosphorylation. Notably, this pattern was not observed in the self-healing mouse model following the removal of ovalbumin stimulation. CONCLUSION: Our results suggest that the persistent activation of antigen processing and presentation pathways in keratinocytes may be a key driver of chronic inflammation in AD. Therefore, redirecting anti-allergic therapeutic strategies from solely targeting immune cells to targeting of keratinocyte-mediated antigen presentation may offer a more effective approach. Furthermore, we raise concerns about the use of ovalbumin-induced mouse models to recapitulate human chronic AD, as the underlying mechanisms may differ significantly.
Dermatitis, Atopic↗
Immune cell-specific genetic architecture of Alzheimer's disease revealed by multi-omics analysis for therapeutic target discovery and prioritization.
Alzheimer's disease (AD) is a multifactorial neurodegenerative condition in which accumulating genetic and molecular evidence implicates dysregulation of peripheral immune processes in disease pathogenesis. Nevertheless, the contribution of distinct peripheral immune cell subsets and associated gene regulatory landscapes to AD risk remains incompletely defined. To address this gap, we integrated single-cell expression quantitative trait loci (sc‑eQTL) data from the OneK1K cohort with AD GWAS summary statistics. We systematically interrogated immune cell-specific genes for their contributions to AD risk by integrating genetic causal inference with Bayesian colocalization analyses, and identified 24 eGenes that passed both the MR significance threshold (P < 0.05) and the criterion for strong shared genetic signals (PP.H4 > 0.8). Notable candidates included GATS, HLA-DOB, HLA-DQA1, PM20D1, and others, with each gene demonstrating a cell-type-specific association restricted to its corresponding immune cell type, such as monocytes, CD8 + T cells, or B cells. Independent peripheral blood single-cell transcriptomic data further supported disease-associated shifts in cell-type-specific expression patterns in AD. Phenome-wide association studies (PheWAS) indicated limited associations with off-target traits, indicating a favorable safety profile for therapeutic intervention, with the exceptions of B4GALNT3, PM20D1, and CNN2. Integration of immune gene targets with pharmacological databases yielded three candidate compound, including NSC321521 (targeting HLA-DQA1), phenoxybenzamine (targeting GSTP1), and rimexolone (targeting BIN1). Among these compounds, Predicted blood-brain barrier permeability was observed only for phenoxybenzamine and rimexolone, with docking studies indicating stable interactions, such as those between NSC321521 and HLA-DQA1, phenoxybenzamine and GSTP1, and rimexolone and BIN1. This integrative approach highlights key immune‑cell‑specific genes involved in AD and proposes repurposable drugs with central nervous system potential, paving the way for more targeted immunomodulatory strategies in AD.
Humans↗
Methylation-Associated Differentiation Features Define Biological and Prognostic Heterogeneity in CMS4 Colorectal Cancer.
Consensus molecular subtype 4 (CMS4) colorectal cancer (CRC) is associated with an aggressive clinical course and poor survival, yet the biological basis of heterogeneity within this subtype remains incompletely understood. DNA methylation is an epigenetic mechanism involved in transcriptional regulation, cellular differentiation, and colorectal tumorigenesis. Here, we integrated single-cell RNA sequencing (scRNA-seq), bulk data, and promoter DNA methylation data to characterize CMS4-associated cancer cell states and methylation-related features. Using the scAB algorithm, we integrated scRNA-seq with bulk CMS4 data and identified CMS4-related cells distributed across multiple patients. Single-cell analyses of cell-cell communication and transcriptional regulation revealed a CMS4-related cancer cell population characterized by macrophage migration inhibitory factor (MIF)-centered intercellular communication, enhanced caudal type homeobox 1 (CDX1) and Kruppel-like factor 5 (KLF5) regulon activity, and gene modules enriched in differentiation-related pathways. CytoTRACE analysis further stratified CMS4 cancer cells into poorly and well-differentiated states, yielding 802 differentially expressed genes (DEGs). Linking these differentiation-associated DEGs with bulk expression and promoter methylation data identified 218 methylation-associated DEGs showing significant inverse methylation expression correlations, suggesting a link between differentiation-related heterogeneity and promoter methylation. Univariable Cox regression followed by LASSO regression further prioritized eight genes for construction of the methylation and differentiation-related prognostic model (MeDiff-PM). MeDiff-PM consistently stratified overall survival in the TCGA CMS4 cohort and two independent validation cohorts, with cutoff-independent continuous Cox analyses further supporting its prognostic association across cohorts. And MeDiff-PM remained prognostically significant after adjustment for available clinical variables. High MeDiff-PM risk scores were associated with activation of P53, WNT, and ubiquitin-mediated proteolysis pathways and with consistent predicted drug response differences for compounds across three CMS4 cohorts. While individual in silico knockout analysis suggested links between MeDiff-PM genes and metallothionein-related and immune-associated transcriptional responses. Collectively, these findings indicate that methylation-associated differentiation features represent a molecular dimension of intra-CMS4 heterogeneity and provide a biologically informed framework for prognostic stratification within CMS4 CRC.
Humans↗
Mapping the Immune cell-specific gene regulatory network in bipolar disorder: A framework from scTWMR to exploratory drug-target annotation.
BACKGROUND: Although the involvement of the immune system in the genetic susceptibility of bipolar disorder (BD) is widely acknowledged, the causal relationship between gene expression in specific immune cell subtypes and BD requires systematic elucidation. METHODS: We implemented an analytical framework integrating single-cell transcriptome-wide Mendelian randomization (scTWMR) with colocalization analysis. This approach utilized cis-expression quantitative trait loci (cis-eQTLs) derived from 14 distinct immune cell types as instrumental variables to interrogate BD genome-wide association study (GWAS) summary statistics (comprising 41,917 cases and 371,549 controls). Subsequent investigations encompassed functional enrichment analysis, protein-protein interaction (PPI) network construction, phenome-wide association study (PheWAS), and performed an exploratory drug-target annotation. RESULTS: Our analysis identified 33 gene-immune cell associations. Colocalization analysis provided robust evidence (PPH4 > 90%) for shared causal variants implicating the MAD1L1, APOM, and NFKBIL1 loci. Significantly enriched biological pathways included cell cycle regulation, circadian rhythm entrainment, and neuroinflammation. The PPI network revealed a core regulatory module centered on histone-encoding and immune-related genes. Exploratory drug-target annotation nominated compounds for further investigation for compounds targeting APOM, TMEM258, and NFKBIL1. CONCLUSION: This study systematically delineates a genetically supported regulatory network of immune cell-specific gene expression in BD, predominantly implicating CD8⁺ effector T cells, plasma cells, and B cells. The findings corroborate established pathological pathways while uncovering novel cell type-specific therapeutic targets, thereby providing a genetic framework for prioritizing candidate targets for future investigation.
Bipolar disorder↗
Protocol to perform integrative analysis of high-dimensional single-cell multimodal data using an interpretable deep learning technique.
The advent of single-cell multi-omics sequencing technology makes it possible for researchers to leverage multiple modalities for individual cells. Here, we present a protocol to perform integrative analysis of high-dimensional single-cell multimodal data using an interpretable deep learning technique called moETM. We describe steps for data preprocessing, multi-omics integration, inclusion of prior pathway knowledge, and cross-omics imputation. As a demonstration, we used the single-cell multi-omics data collected from bone marrow mononuclear cells (GSE194122) as in our original study. For complete details on the use and execution of this protocol, please refer to Zhou et al.1.
Deep Learning↗