PubMed HealthSearch

SEARCH · PubMed Health

Results for “Single-nucleus transcriptomics”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 19 recordsLinked to original sources

A single-nucleus transcriptomic atlas of human inner ear development.

Hearing and balance rely on coordinated activity of multiple inner ear cell types, yet the mechanisms governing their development and specification in humans remain unclear. Consequently, this limits our understanding of how disease genes affect cell type formation and function, limiting the development of targeted treatments, including gene therapies. Here we present the Human Inner Ear Development snRNA-seq Atlas (HIEDRA), a single-nucleus transcriptomic atlas of the human inner ear spanning the first and second trimesters. HIEDRA maps sensory and nonsensory epithelia, neurons and mesenchyme-associated populations, including undercharacterized secretory cells required for ion homeostasis. We identify selective vulnerability in sensory and secretory lineages to disease-associated genes, infer regulatory networks and show that Hedgehog signaling suppression is required for secretory cell specification. We validate this mechanism in human inner ear organoids, expanding the model to include all major cell types. Altogether, these findings provide insights into human inner ear cell type specification, improve in vitro models and establish HIEDRA as a resource for investigating human inner ear development.

Journal Article

Single-nucleus transcriptomics reveals cell type-specific remodeling and epilepsy-associated microglia.

Temporal lobe epilepsy (TLE) is the most common acquired epilepsy, causing refractory seizures and cognitive deficits. We performed single-nucleus RNA sequencing on hippocampal tissue from mice 3 and 6 weeks following pilocarpine-induced status epilepticus, a robust model of TLE. Epilepsy samples showed reductions in Cck and Lamp5-Lhx6 interneuron subclusters, alongside increases in Cajal-Retzius cells, dentate granule (DG) cell precursors, and a mature DG cell subcluster. Among glia, an astrocyte subcluster and a markedly expanded microglia sublcuster were increased. We term this microglia population epilepsy-associated microglia (EAM). The transcriptomic profile of EAM overlaps with microglia described in models of Alzheimer's disease and traumatic brain injury, including enrichment of Myo1e and Igf1. EAM display amoeboid morphology, can be found in clumps around pyramidal and granule cell body layers, and exhibit enlarged vesicles and mitochondria. Cell-cell interaction analysis predicts DG cells as their primary interaction partners. This dataset defines transcriptomic programs underlying key cellular alterations in TLE, enabling mechanistic dissection of epileptogenesis.

TLE

A single-nucleus transcriptome atlas of soybean anthers.

Anther development is crucial for plant sexual reproduction. However, a high-resolution, cell-type-specific transcriptomic atlas of this process is lacking for the legume crop soybean (Glycine max). Here, we construct a comprehensive transcriptional atlas of developing soybean anthers using single-nucleus RNA sequencing (snRNA-seq). We identify and characterize nine distinct cell types spanning both somatic and reproductive lineages. Our analysis reveals robust transcriptional continuity across anther developmental stages and dynamic reprogramming during key transitions. Notably, the shift from diploid meiocytes to haploid unicellular microspores is marked by the induction of previously inactive genes, despite an overall reduction in transcript abundance. Subsequently, within bicellular microspores, generative and vegetative cell lineages exhibit sharply divergent transcriptional programs: generative cells specialize in mRNA export and turnover, whereas vegetative cells up-regulate translational machinery. Evolutionary analysis further indicates that generative-cell-specific genes are subject to more relaxed purifying selection compared to those specific to vegetative cells. Functional validation using mutants generated by CRISPR/Cas9-mediated genome editing and EMS mutagenesis reveals the essential roles of OSD1A and PKSA in pollen development and fertility. This high-resolution atlas provides fundamental insights into the transcriptional regulation of soybean anther development and serves as a valuable resource for manipulating male fertility to advance hybrid breeding programs. The data are available at https://databases.genedenovo.com/pollen.

Glycine max

Precision-Based Filtering Facilitates Cross-Referencing of Conventional and Single-Nucleus Transcriptomes to Identify Time- and Temperature-Sensitive Cell Populations.

Transcriptome analysis via RNA sequencing (RNAseq) has become a ubiquitous method of molecular characterization from whole organisms, dissected tissues, and single cells. These experiments continue to provide an extraordinary volume of data describing molecular states and responses to many conditions. However, standard approaches to RNAseq analysis commonly use expression level filters that eliminate potentially useful data in the service of decreasing noise. Here we describe the implementation of a coefficient of variation-based filter for RNAseq gene expression data. This filter prioritizes consistent data across replicates, allowing lowly-expressed genes with low-variation measurements to be retained for downstream analysis. We show, using two independent Arabidopsis RNAseq datasets, that this filter allows for the inclusion of many more transcription factors than even a low-stringency expression level filter. This effect is independent of sequencing depth. We find that these lowly-expressed genes mark specific cell clusters in our single-nucleus (sn)RNAseq dataset and may facilitate future characterization of currently unknown cell types or states. We further characterize communities of co-expressed genes, sampled across the day at two growth temperatures, in relation to snRNAseq cell clusters, finding evidence for a highly photosynthetic cell population, and a cell state marked by high cell division and translation. These methods can be expanded to RNAseq analysis in many systems, facilitating the construction of more detailed models of tissue-specific gene regulatory networks.

Transcriptome analysis

Integrative TWAS and multi-omics analyses prioritize HSPE1 as a candidate risk gene for bipolar disorder with immune cell-specific regulatory evidence.

BACKGROUND: Bipolar disorder (BD) is a severe psychiatric disorder associated with substantial disability. Although genome-wide association studies have identified multiple BD-associated loci, the underlying genes and mechanisms remain incompletely understood. METHODS: We integrated a European-ancestry BD genome-wide association dataset with cross-tissue and tissue-specific transcriptome-wide association studies (TWAS) and complementary gene-based analysis. Candidate genes were further evaluated using differential expression analysis, consensus clustering, immune infiltration analysis, machine learning, summary-data-based Mendelian randomization, Mendelian randomization using single-cell expression quantitative trait locus data, single-nucleus transcriptomics, phenome-wide association analysis, and virtual screening. RESULTS: The integrative analyses prioritized 37 candidate genes. Peripheral-blood differential-expression analysis identified 14 genes that remained significant after FDR correction, and their expression profiles separated BD samples into two expression-defined clusters. Machine-learning analysis selected UNC50, LMAN2L, LYG2, HSPE1, and KANSL3 for an exploratory classification nomogram. SMR associated genetically predicted higher HSPE1 expression with increased BD risk in two blood eQTL datasets. Cell-type-specific analyses indicated HSPE1-related associations in T-cell and natural killer cell subsets, while single-nucleus analysis descriptively showed higher HSPE1 expression in medial thalamic T cells from BD samples. PheWAS identified no genome-wide significant associations for HSPE1, whereas virtual screening identified candidate compounds with favorable predicted docking scores against the HSPE1 structure. CONCLUSION: This integrative multi-omics study identified HSPE1 as a candidate BD risk gene with immune-cell-related regulatory evidence, providing insight into BD pathogenesis and supporting functional validation.

Humans

Distinct cellular phenotypes of language and executive decline in amyotrophic lateral sclerosis.

Cognitive manifestations, including impairments in language and executive functions, are seen in amyotrophic lateral sclerosis (ALS), but the underlying mechanisms remain unclear. We mapped prefrontal cortex regions from ALS patients by integrating spatial and single-nucleus transcriptomics in a cognitively stratified patient cohort. We uncover that cognitive impairment in ALS is associated with distinct patterns of neuronal dysfunction and glial-vascular dysregulation that vary by region and cognitive subtype. Executive dysfunction is linked to reduced mitochondrial and synaptic activity in deep-layer dorsolateral prefrontal cortex neurons, whereas language-related deficits track with a diffuse pan-regional response involving glial and vascular abnormalities. Our analyses, validated by multiplexed imaging, further identify signatures in the prefrontal cortex that span both motor and cognitive phenotypes, including a multicellular gliosis response. The findings reveal that clinical heterogeneity in ALS is driven by phenotype-specific cellular interactions in motor and non-motor regions of the brain.

Amyotrophic Lateral Sclerosis

Spatial profiling of the spleen in mouse and human myelofibrosis reveals complement-driven immune-stromal interactions as a therapeutic target.

Splenomegaly is a defining feature of myelofibrosis, yet the contribution of splenic mesenchymal stroma to disease progression remains unclear. We combined spatial and single-nucleus transcriptomics of patient spleens with spatial and single-cell transcriptomics, as well as imaging analyses, of murine spleens to map extramedullary hematopoiesis niches. Activated red pulp reticular cells localize near hematopoietic stem and progenitor cells, and early disease is characterized by marginal zone disruption with lymphoid depletion preceding stromal remodeling. Trajectory analyses reveal a shift in reticular cells from hematopoiesis-supportive to inflammatory and pro-fibrotic states, driven by macrophage- and megakaryocyte-derived signals that activate complement and induce tumor necrosis factor α (TNF-α), transforming growth factor β (TGF-β), extracellular matrix, and Thbs1 programs. Non-hematopoietic complement component C3 deficiency or pharmacological C3 inhibition suppresses these pathways, restores splenic architecture, and reduces splenomegaly and bone marrow fibrosis. These findings identify complement-dependent stromal reprogramming as a mechanism governing hematopoietic niches and as a targetable axis in myelofibrosis.

Animals

A single-nucleus and spatial transcriptomic atlas of poplar leaves reveals the regulation of leaf polarity and cuticle deposition.

Leaf adaxial-abaxial polarity is fundamental for plant morphogenesis and environmental adaptation through asymmetric cell differentiation. Emerging evidence reveals dorsoventral metabolic gradients act downstream of transcriptional networks to fine-tune cellular specialization. While conserved transcription factors (e.g., HD-ZIP III and KANADI) establish initial polarity, the molecular networks driving position-specific cellular differentiation and their integration with metabolic adaptation remain unclear. Leveraging single-nucleus and spatial transcriptomics, we resolve major cell classes (mesophyll, epidermal, and vascular-associated) and their adaxial-abaxial subtypes, revealing dorsoventral polarity in transcriptional profiles and metabolic pathways. Adaxial cells are enriched in phenylpropanoid/flavonoid biosynthesis, while abaxial cells show preferential activation of stress and hormone signaling. Notably, we identify MYC2 as a key regulator of adaxial cuticle biosynthesis, binding to promoters of lipid biosynthetic and transport genes (e.g., CER10 and LTPG1) and promoting cuticle thickening. Our study uncovers how positional identity shapes transcriptional and metabolic polarity in leaves, with MYC2 emerging as a central regulator coordinating organ-specific adaptations. These findings provide insights into the spatial regulation of plant development and stress resilience, offering potential strategies for engineering stress-tolerant woody crops.

Plant Leaves

Decoding the landscape of cell-type-specific co-expressed transcription factors in soybean.

Soybean (Glycine max) is an essential source of protein and oil with high nutritional value for human and animal consumption. To enhance our understanding of soybean biology, it is essential to have accurate information regarding the expression of each of its protein-coding genes. Here, we present Tabula Glycine max, a soybean single-cell resolution transcriptome atlas. This atlas comprises single-nucleus RNA-sequencing data from ten different G. max organs and morphological structures constituting the entire soybean plant. These nuclei are grouped into 156 different clusters based on their transcriptomic profiles. The breadth of various organs, tissues and cell types represented in Tabula Glycine max reveals that the pattern of co-expressed transcription factor genes is sufficient to define most cell types based on their function and organ of origin. Defining cell-type-specific co-expressed transcription factor genes offers a new perspective to engineer cell-type-specific programmes and enhance the biology of unique soybean cell types. This cellular resolution and breadth make the Tabula Glycine max an exceptional resource for the plant and soybean communities.

Journal Article

A Molecularly Anchored Spatial Transcriptomic Framework for Precise CA1-Subiculum Parcellation and Region-Resolved Analysis in Alzheimer's Disease.

BACKGROUND: The precise molecular delineation of the interface between the Subiculum (Sub) and cornu ammonis 1 (CA1) is a challenge in hippocampal research, as conventional cytoarchitectural boundaries are often ambiguous and limit reproducible regional annotation. Here, we developed a molecularly anchored spatial transcriptomic framework to define CA1-Sub regional identities using high-definition spatial transcriptomics (Stereo-seq) and single-nucleus RNA sequencing (snRNA-seq) references. FINDINGS: Using a human hippocampal Stereo-seq dataset from 12 donors, we established a data-driven parcellation framework that defines reproducible molecular features distinguishing CA1 and Sub while capturing the transition between these regions. FN1 was identified as a Sub-enriched marker in a subset of EX_Sub and, together with ETV1 and additional regional markers, enabled molecular assignment of CA1 and Sub identities across datasets. The Sub association of FN1 and ETV1 was further supported by human 10X Genomics spatial transcriptomics, mouse in situ hybridization data, and a mouse spatial transcriptomic dataset. Applying this framework to Alzheimer's disease (AD) tissues revealed region-specific transcriptional alterations across CA1 and Sub, including enrichment of mitochondrial energy metabolism-related transcripts in the Sub, suggesting exploratory transcriptional associations of altered metabolic function. CONCLUSIONS: This study provides a molecularly anchored framework for human CA1-Sub parcellation that complements conventional annotation. By defining regional molecular states while preserving the biological continuum across CA1-Sub interface, this approach enables more consistent regional analysis of human hippocampus tissue across donors, datasets, and disease conditions.

Journal Article

Combined somatic mutation and transcriptome analysis reveals region-specific differences in clonal architecture in human cortex.

The human cerebral cortex is specialized into regions, but little is known about how human cellular lineages shape cortical regional variation and neuronal cell-type distribution during development. Here, we map single-cell lineages of human cortical regions and neuronal subtypes using >1,000 somatic single-nucleotide variants (sSNVs) identified from deep bulk whole-genome sequencing and analyzed over 25 regions and >72,000 single cells. In the fronto-parietal cortex, sSNVs are rarely restricted, marking neuron-generating clones that disperse into neighboring regions. In contrast, the primary visual cortex harbors 30%-70% more sSNVs than the neighboring secondary visual cortex. Clones at this border exhibit more restricted dispersion, suggesting late developmental lineage segregation. Single-nucleus sSNV and whole-transcriptome analysis reveal glutamatergic neuron clones with modest regional restrictions that share low-mosaic sSNVs with some GABAergic neurons, suggesting a recent dorsal cortical progenitor. Our analysis reveals human-specific cortical lineage patterns, regional differences in clonal patterns, and late divergence of some glutamatergic/GABAergic lineages.

Humans

Genetics-Informed Mapping Identifies a CRIM1-Associated Endocardial Inflammatory Remodeling State in Acute Myocardial Infarction.

BACKGROUND Acute myocardial infarction (AMI) reflects inherited susceptibility and inflammatory remodeling, but the cellular contexts linking genetic risk to disease remain unclear. MATERIAL AND METHODS We integrated a meta-transcriptome-wide association study (TWAS) with a human cardiac single-nucleus RNA-sequencing atlas contained 11 individuals (5 AMI and 6 donor) to identify genetics-informed cellular programs. Composite program states were defined by global score quartiles. A fixed 5-gene panel was evaluated for nucleus-level endocardial low-transcriptional-state (Endo_LTS) vs endocardial high-transcriptional-state (Endo_HTS) discrimination within the AMI endocardium using 5-fold leave-1-patient-out cross-validation. Functional follow-up used CRIM1 silencing in hypoxia-treated human induced pluripotent stem cell (hiPSC)-derived endocardial endothelial-like cells and complementary peripheral blood analyses. RESULTS The endocardium exhibited the most prominent infarction-associated increase in TWAS-anchored program activity, with expansion of program-high states and higher CytoTRACE scores. A consensus 5-gene panel (RPS8, PLEC, CFDP1, CRIM1, TNS2) was identified. Among 2163 AMI endocardial nuclei from 5 patients, the state classifier included 364 Endo_LTS and 751 Endo_HTS nuclei; 1048 Endo_MTS nuclei were excluded. Pooled out-of-fold ROC-AUCs ranged from 0.665 to 0.831. The panel also showed discriminatory value in an independent peripheral-blood AMI-vs-control cohort. CRIM1 was prioritized as a candidate linked to the remodeling program. CRIM1 silencing attenuated ACTA2/alpha-SMA, vimentin, LDHA, CCL2, and VEGFA and partially restored CD31, whereas TGF-ß remained elevated. CONCLUSIONS These findings identify a genetics-informed endocardial inflammatory remodeling state in AMI and define a 5-gene surrogate of its activated state. CRIM1 is prioritized as a candidate linked to selected inflammatory, metabolic, and structural outputs. Persistent TGF-b elevation after CRIM1 silencing argues against a simple linear regulatory model and indicates that further mechanistic validation is required.

Humans

Big data analytics for CLEC5A dynamics based on single cell genomics and proteomics reveal its diverse functions in human diseases.

BACKGROUND: CLEC5A (C-type lectin domain family 5 member A) is an innate immune receptor implicated in inflammatory signaling, contributing to hyperinflammatory responses in infections and sterile inflammation. However, CLEC5A dynamics in human diseases remain to be identified. Here, we systematically characterized CLEC5A dynamics in humans across cells, tissues, and disease states, and to explore the functional significance of CLEC5A in macrophage activation based on single-cell genomics. METHODS: With multi-omics (scRNA-seq, proteomics and big data analytics), we analyzed extensive human transcriptomic datasets (>42,000 samples) to profile CLEC5A expression by cell type, tissue, and disease. Single-nucleus RNA-seq (snRNA-seq) from pediatric congenital heart disease and a virtual CLEC5A gene knockout were also performed to characterize CLEC5A dynamics in humans. RESULTS: CLEC5A is highly enriched in innate immune cells, particularly in macrophages and neutrophils. Baseline CLEC5A in most tissues is low, but it is markedly upregulated in inflammatory and infectious diseases. CLEC5A expression has sex-specific differences in certain organs. Single-cell analysis showed that CLEC5A can be considered novel marker of proinflammatory macrophages with elevated cytokine production, antigen presentation, and impaired phagocytosis. Virtual CLEC5A knockout analysis identified coordinated perturbation of immune-regulatory pathways and overlapping genes linking CLEC5A to macrophage activation networks. CONCLUSION: CLEC5A is predominantly expressed in myeloid cells and acts as a key amplifier of inflammation in human diseases. Our findings highlight CLEC5A as a potential biomarker and therapeutic target in myeloid-driven hyperinflammatory conditions, warranting further experimental and translational validation.

Humans

From wild to domestic: Single-cell transcriptomic perspectives on hippocampal regulation and evolution.

How domestication shapes brain evolution remains an open question. In this study, we integrated single-nucleus RNA sequencing (snRNA-seq), population genomics, and machine learning to investigate the hippocampal evolution under domestication. Across-species comparisons revealed that hippocampal cell type profiles are largely conserved across vertebrate species, while supporting the presence of adult hippocampal neurogenesis in birds. We further found that domestication and selective breeding likely influence the cellular composition and molecular regulation of the hippocampus. Our findings provide cellular evidence supporting the hypothesis that domestication affects adult hippocampal neurogenesis. Additionally, we showed that genes associated with neural progenitor cells (NPC) states and cell-marker programs are enriched for signatures of selection. Many of these genes function as regulators of neurogenesis and pathways mediating stress and fear reduction. Specifically, we identified selection at the FKBP5 promoter that may influence its expression in the NPC lineage, potentially contributing to stress-response regulation during domestication. Collectively, these results suggest that domestication is associated with hippocampal remodeling as part of an adaptive response to human-managed environments. This study provides a cellular and genetic perspective on how domestication reshapes the brain and offers a basis for further investigation into the mechanisms of neural evolution within the context of microevolution.

Animals

Spatiotemporal mapping of tertiary lymphoid structure heterogeneity shapes immune niches and clinical outcomes in intrahepatic cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (iCCA) is a highly lethal malignancy with limited therapeutic options. The spatial architecture and functional diversity of tertiary lymphoid structures (TLSs) in iCCA remain unclear. Here, we present a multimodal spatial atlas of TLSs and identified intratumoral TLSs (iTLSs) as independent prognostic markers. Bulk proteomic profiling of 214 discovery and 155 validation cases identified a four-tier TLS-based tumor microenvironment classification system and supported development of a TLS-predictive random forest classifier. Imaging mass cytometry revealed that iTLS+ tumors harbor structured immune architectures, where M1-like tissue-resident macrophages (RTMs), dendritic cells, and CXCL13+ CD4+ T cells colocalize to form antigen-presenting neighborhoods (apc-CNs) spatially coupled to TLS core regions (TLScore-CNs). Single-cell spatial transcriptomics further resolved 61 TLSs into 14 spatial niches and defined a pseudotemporal maturation continuum: aggregated, activated, and postactivated. Intraniche communication, primarily mediated by ifnCAFs, iCAFs, and CXCL12+ macrophages, evolved dynamically with maturation. Single-nucleus RNA sequencing combined with Tangram-based spatial mapping revealed CXCL12+ macrophages and iCAFs forming a peripheral band in aggregated TLSs, whereas ifnCAFs infiltrated TLS interiors during activation. These findings define TLS heterogeneity and provide insights for stroma-directed immunotherapy.

Cholangiocarcinoma

Wound healing in Atlantic spiny dogfish sharks.

Field observations and limited experimental studies indicate that elasmobranchs can repair substantial skin injuries, but the temporal course and cellular composition of wound healing in Atlantic spiny dogfish remain poorly characterized. We conducted an exploratory laboratory study in 20 female Atlantic spiny dogfish (Squalus acanthias) using standardized full-thickness skin wounds monitored by serial photography for 35 days, histological analysis at defined post-injury time points, and pooled single-nucleus RNA sequencing of intact and wounded skin. A continuous neoepithelial layer covered all examined wound beds by Day 1, whereas macroscopic wound area decreased progressively over 35 days and dermal denticles remained absent from the repaired surface. Histological examination showed progressive neoepidermal maturation, basement-membrane reformation, collagen deposition, and granulation-tissue organization, indicating that epithelial coverage preceded restoration of normal skin architecture. Single-nucleus RNA sequencing identified epithelial, stromal, vascular, pigment, neural, and immune-cell populations. T and B cells were detected in intact skin, and their relative abundance, together with that of several other leukocyte populations, increased at Day 1 and generally declined by Day 14. Because samples were pooled by time point, these transcriptomic changes are descriptive. These findings characterize rapid early reepithelialization followed by slower tissue remodeling in Atlantic spiny dogfish and provide a foundation for future comparative studies of elasmobranch skin repair.

Animals

Complementation testing identifies genes mediating effects at quantitative trait loci underlying fear-related behavior.

Knowing the genes involved in quantitative traits provides an entry point to understanding the biological bases of behavior, but there are very few examples where the pathway from genetic locus to behavioral change is known. To explore the role of specific genes in fear behavior, we mapped three fear-related traits, tested fourteen genes at six quantitative trait loci (QTLs) by quantitative complementation, and identified six genes. Four genes, Lamp, Ptprd, Nptx2, and Sh3gl, have known roles in synapse function; the fifth, Psip1, was not previously implicated in behavior; and the sixth is a long non-coding RNA, 4933413L06Rik, of unknown function. Variation in transcriptome and epigenetic modalities occurred preferentially in excitatory neurons, suggesting that genetic variation is more permissible in excitatory than inhibitory neuronal circuits. Our results relieve a bottleneck in using genetic mapping of QTLs to uncover biology underlying behavior and prompt a reconsideration of expected relationships between genetic and functional variation.

Animals

A single-cell meta-analysis evidences transposable element dysregulation in sex-based differences in Parkinson's disease.

Transposable elements (TEs) (mobile genetic elements comprising ∼45% of the human genome) have recently emerged as potential contributors to Parkinson's disease (PD); however their role and sex-specific impact remain poorly understood. Here, we present the first integrative meta-analysis of TE expression across 4 substantia nigra single-nucleus RNA-seq datasets, comprising a total of 66 donors, generating a cell-type-resolved atlas of TE dysregulation in PD. We identified widespread TE activation across major brain cell types (i.e. neurons, astrocytes, oligodendrocytes and microglia), with marked upregulation of L1s in neurons and HERVs in oligodendrocytes. Sex-stratified analyses revealed distinct male- and female-biased TE signatures, indicating regulatory programs uniquely affected in each sex, including MIR elements in microglia and Alu subfamilies in neurons. Correlation and genomic proximity analyses also uncovered TE-gene associations linked to important PD pathways such as neuroinflammation or myelination. Collectively, our study positions TEs as potential sex-modulated contributors to PD pathology and also provides a public web resource (PATOSS) to explore PD-associated TE transcriptional deregulation.

Parkinson's disease