PubMed HealthSearch

SEARCH · PubMed Health

Results for “single-nucleus sequencing”

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.

48 records · Page 3Linked to original sources

Evolutionary Process Underlying Receptor Gene Expansion and Cellular Divergence of Olfactory Sensory Neurons in Honeybees.

Olfaction is crucial for animals' survival and adaptation. Unlike the strict singular expression of odorant receptor (OR) genes in vertebrate olfactory sensory neurons (OSNs), insects exhibit complex OR gene expression patterns. In honeybees (Apis mellifera), a significant expansion of OR genes implies a selection preference for the olfactory demands of social insects. However, the mechanisms underlying receptor expression specificity and their contribution to OSN divergence remain unclear. In this study, we used single-nucleus multiomics profiling to investigate the transcriptional regulation of OR genes and the cellular identity of OSNs in A. mellifera. We identified three distinct OR expression patterns, singular OR expression, co-expression of multiple OR genes with a single active promoter, and co-expression of multiple OR genes with multiple active promoters. Notably, ∼50% of OSNs co-expressed multiple OR genes, driven by polycistronic transcription of tandemly duplicated OR genes via a single active promoter. In these OSNs, their identity was determined by the first transcribed receptor. The divergent activation of the promoter for duplicated OR genes ensures the coordinated increased divergence of OSN population. By integrating multiomics data with genomic architecture, we illustrate how fundamental genetic mechanisms drive OR gene expansion and influence flanking regulatory elements, ultimately contributing to the cellular divergence of OSNs. Our findings highlight the interplay between gene duplication and regulatory evolution in shaping OSN diversity, providing new insights into the evolution and adaptation of olfaction in social insects. This study also sheds light on how genetic innovations contribute to the evolution of complex traits.

Animals

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

Single-nucleus profiling of postmortem diffuse midline gliomas identifies mitochondrial biogenesis as a resistance mechanism to imipridone therapy.

BACKGROUND: Imipridone ONC201 is the first FDA-approved therapy for H3K27-altered diffuse midline glioma; however, clinical responses remain limited. Defining tumor-intrinsic determinants and microenvironmental, extrinsic factors that shape sensitivity or resistance to imipridones will identify actionable therapeutic opportunities and inform improved clinical strategies. METHODS: To identify mechanisms of imipridone resistance, we obtained postmortem brain tissue from DMG patients who had received imipridones and/or standard care. Single-nucleus RNA and open-chromatin sequencing were performed on N = 22 cases. Immunofluorescence-based myeloid phenotyping was performed on N = 46 cases. Mitochondrial copy-number analysis was performed on N = 19 cases. Validation of imipridone sensitivity, its effect on mitochondrial density, and its synergy with inhibition of mitochondrial biogenesis were assessed in DMG primary cells. RESULTS: We established a single-cell RNA/open-chromatin atlas from postmortem DMG cases and found imipridone treatment resulting in regressed mesenchymal transition, reduced myeloid-derived suppressive cells, and reversed aberrant H3K27-altered enhancer activity. Resistant tumors showed increased mitochondrial density, turnover, and membrane potential. Mitochondrial biogenesis and PPARGC1A emerged as resistance biomarkers and actionable targets. CONCLUSIONS: These studies implicate mitochondrial biogenesis as a biomarker of imipridone resistance and a focus for the development of combinatorial strategies to provide effective therapeutic options for a challenging pediatric brain tumor.

Humans

Alignment-free integration of single-nucleus ATAC-seq across species with sPYce.

Changes in gene regulation largely contribute to differences in cellular identities and phenotypes between species. Single-nucleus assays for transposase-accessible chromatin with sequencing (snATAC-seq) are an efficient strategy to identify putative gene regulatory elements and provide new insight into evolutionary divergence of regulatory programmes. However, no dedicated framework exists to integrate and compare snATAC-seq data across species, while methods designed for single-cell gene expression data have serious limitations. Here we present sPYce, a cross-species snATAC-seq integration method that relies on sequence composition similarities through k-mer histograms of regulatory regions, removing the need for genome alignments to anchor data from different species. sPYce can embed datasets from multiple species into the same mathematical space and permits further downstream analysis steps. We benchmarked sPYce against existing approaches on two publicly available datasets spanning more than 160 myr of evolution, showing that it successfully uncovers conserved cellular programmes while preserving biologically relevant species-specific differences. By comparing cerebellar development in mice and opossums, sPYce identifies regulatory divergence in granule cell differentiation programmes, particularly driven by nuclear factor 1. As an easy-to-use, alignment-free cross-species snATAC-seq integration approach, sPYce opens new perspectives to compare gene regulatory evolution across species.

Animals

A complete and near-perfect rhesus macaque reference genome: lessons from subtelomeric repeats and sequencing bias.

A truly complete, telomere-to-telomere (T2T), and error-free reference genome remains a foundational resource-and long-standing goal-for unbiased comparative and functional genomics. While recent T2T assemblies of humans and other primates have made substantial progress, most still contain thousands of base-level errors, particularly within highly repetitive regions. Here, we present T2T-MMU8v2.0, a near-perfect T2T assembly of the rhesus macaque (Macaca mulatta), representing the highest base-level accuracy reported in a primate genome to date. By employing an optimized ONT-only assembly strategy, we identify subtelomeric satellite-rich regions as the principal bottleneck to improving assembly quality, owing to technological biases in long-read platforms and limitations in current hybrid assembly frameworks. We discover 268 previously unannotated repeat families and resolve ~8 Mbp of SATR satellite arrays, with over 99-fold enrichment in historically misassembled subtelomeric regions. These satellites form four distinct genomic architectures, each with unique SATR satellite composition, segmental duplication organization, and epigenetic signatures, distinct from the subtelomeric architectures observed in hominid genomes. Notably, in contrast to the largely gene-poor subtelomeric regions in African hominids, the SATR architectures in macaques harbor 58 actively transcribed genes, supported by open chromatin and expression data, suggesting gene innovation within these repetitive regions. Functionally, T2T-MMU8v2.0 improves read mappability and accuracy across sequencing platforms, and results in a 19% improvement of transcription start site enrichment scores and 5,821 additional chromatin accessibility peaks on average, thereby enhancing variant detection, regulatory annotation, and transcriptomic resolution in population genetics or single-nucleus studies. Together, this work establishes a new benchmark for genomics, offers a roadmap for resolving complex repetitive regions, and reveals previously unrecognized features of subtelomeric genome structure and evolution.

Journal Article

Multiomic single-nucleus profiling reveals cell-type-specific epigenetic and transcriptional dysregulation in major depressive disorder brain.

OBJECTIVE: Major depressive disorder (MDD) is a leading global cause of disability, marked by persistent mood disturbances, cognitive deficits, and changes in prefrontal cortex neural circuitry. In this study, we aimed to define cell-type-specific molecular and regulatory mechanisms underlying MDD by mapping gene-expression and chromatin-accessibility changes in the dorsolateral prefrontal cortex (PFC) (dlPFC). METHODS: Postmortem dlPFC (BA9) tissue from 7 MDD and 8 well-matched controls was analyzed using 10× Genomics snRNA-seq and paired ATAC+RNA multiome sequencing. Sequencing data were processed with Cell Ranger pipelines, nuclei were filtered for quality and doublets/debris, and datasets were integrated and clustered using Seurat/Signac packages. Differential gene expression, chromatin accessibility, and transcription factor motif activity were tested between MDD and controls within each cell type, followed by peak-to-gene linkage and Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and PsyGeNET enrichment to interpret dysregulated regulatory mechanisms. RESULTS: A total of 20 distinct clusters encompassing major neuronal and non-neuronal populations were identified. Differential analyses uncovered extensive cell type-specific changes in chromatin accessibility and gene expression, particularly within excitatory layer 5/6 and inhibitory Pvalb neurons, as well as glial and vascular populations. Functional enrichment indicated dysregulation of synaptic organization, neurotransmission, myelination, stress-response, and immune-regulatory pathways across neuronal and non-neuronal cells. Notably, glucocorticoid-responsive transcription factors NR3C1/NR3C2 exhibited conserved regulatory networks implicating stress signaling in MDD pathophysiology. CONCLUSIONS: Together, these findings provide a comprehensive single-nucleus atlas of gene regulation in the MDD PFC, highlighting coordinated dysfunction across neurons, glia, and vascular cells.

Major Depressive Disorder

Drought recovery in plants triggers a cell-state-specific immune activation.

All organisms experience stress as an inevitable part of life, from single-celled microorganisms to complex multicellular beings. The ability to recover from stress is a fundamental trait that determines the overall resilience of an organism, yet stress recovery is understudied. To investigate how plants recover from drought, we examine a fine-scale time series of RNA sequencing starting 15 min after rehydration following moderate drought. We reveal that drought recovery is a rapid process involving the activation of thousands of recovery-specific genes. To capture these rapid recovery responses in different Arabidopsis thaliana (A. thaliana) leaf cell types, we perform a single-nucleus transcriptome analysis at the onset of drought recovery, identifying a cell type-specific transcriptional state developing independently across cell types. To further validate the cell-type specific transcriptional changes observed during drought recovery, we employ spatial transcriptomics using multiplexed error-robust fluorescence in situ hybridization (MERFISH), revealing anatomical localization of recovery-induced gene expression programs across Arabidopsis leaf tissues. Furthermore, we reveal a recovery-induced activation of the immune system that occurs autonomously, and which enhances pathogen resistance in vivo in A. thaliana, wild tomato (Solanum pennellii) and domesticated tomato (Solanum lycopersicum cv. M82). Since rehydration promotes microbial proliferation and thereby increases the risk of infection, the activation of drought recovery-induced immunity may be crucial for plant survival in natural environments. These findings indicate that drought recovery coincides with a preventive defense response, unraveling the complex regulatory mechanisms that facilitate stress recovery in different plant cell types.

Arabidopsis

De Novo TRIO Missense Variants Disrupt Ras-GEF Domains and Cause Congenital Ventriculomegaly and Hydrocephalus.

Congenital hydrocephalus (CH), characterized by congenital ventriculomegaly (CV), affects approximately 0.5-1 per 1000 live births and is a common cause of pediatric neurosurgical intervention, yet its genetic architecture remains incompletely defined. We report a child with syndromic CH requiring cerebrospinal fluid diversion who harbored a pathogenic de novo missense variant in TRIO (c.3232C > T; p.(Arg1078Trp)), a gene previously associated with autosomal dominant neurodevelopmental disorders featuring variable head circumference. This case prompted systematic evaluation of TRIO variation in our CV/CH cohort (2,697 patient-parent trios) using exome sequencing. We identified five additional unrelated probands with de novo TRIO variants, including two novel substitutions affecting the same residue within the Ras-GEF1 domain (p.(Glu1299Lys) and p.(Glu1299Gly)), yielding significant gene-level enrichment for protein-damaging de novo variants (adjusted p = 6.12 × 10-5). All affected individuals exhibited CV, frequently accompanied by developmental delay and additional structural brain abnormalities. In silico structural modeling predicted that associated variants destabilize critical TRIO Ras-GEF domains required for Rho GTPase activation. Analysis of single-nucleus transcriptomic data from the developing human neocortex revealed enrichment of TRIO expression in multipotent progenitor populations. A systematic literature review identified six additional individuals with TRIO de novo variants and reported CV or CH, including an unrelated patient with the same p.(Arg1078Trp) substitution. Together, these findings expand the phenotypic spectrum associated with pathogenic TRIO variation to include CV/CH and support TRIO as a clinically relevant gene in the genetic evaluation of syndromic CV/CH patients.

Humans

Positional grammar of transcription factor binding partitions developmental and stress-response regulation in plants.

Understanding how transcription factor binding site (TFBS) position influences gene regulation remains a fundamental challenge in plants. Here, we integrate conserved multiDAP TFBS maps for 244 transcription factors (TFs) with single-nucleus chromatin accessibility, cell type-resolved gene expression, and hormone-response datasets across Brassicaceae species to determine how TFBS position relates to regulatory function. Although conserved TFBSs are enriched near transcription start sites (TSSs), TSS-proximal accessibility poorly predicts cell type-specific expression. Instead, cell type-specific expression correlates best with conserved TFBSs embedded in cell type-restricted chromatin, with TF family-specific distributions across distal promoters and introns. In contrast, TSS-proximal TFBSs in broadly accessible chromatin are associated with rapid transcriptional responses to abiotic and biotic stress hormones. Coding sequence TFBSs mark a distinct regulatory context in which the same DNA sequence encodes both amino acid sequence and TF motifs, including evidence that CDS-localized ABR1 binding may contribute to repression during hormone response. Finally, distal upstream regions contain conserved multi-family TF clusters with enhancer-like features overlapping rare cell type-specific accessible chromatin and enriched near genes controlling embryonic, meristematic, and hormone-dependent developmental patterning. Together, these results support a positional grammar in which TFBS position and chromatin context jointly partition developmental, stress-responsive, and repressive regulatory output in plants.

Transcription Factors

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

Decoding ALS from the tail end of RNA.

In this issue of Cell Genomics, McKeever et al.1 generate a single-nucleus transcriptomic atlas of ALS/FTLD brain and reveal widespread alternative polyadenylation changes. Their findings highlight 3' end RNA processing as a central integrator of stress responses, cell-type specificity, and disease susceptibility, offering new mechanistic insight and potential therapeutic directions.

Cell Nucleus

VINE-seq and MultiVINE-seq for single-nucleus and multiome profiling of the brain vasculature.

The human cerebrovasculature is a critical yet historically understudied component of neurological health. Dysfunction of the diverse endothelial, mural, and perivascular cells that comprise cerebral vessels is central to diseases ranging from stroke to Alzheimer's disease. However, characterizing these cell populations at a molecular level has proven exceptionally challenging. Encased within a robust basement membrane, vascular cells resist standard dissociation methods, leading to their systematic depletion and underrepresentation in existing single-nucleus genomic atlases. This has created a major blind spot in neuroscience. To overcome this barrier, we developed vessel isolation and nucleus extraction for sequencing (VINE-seq) and its advanced iteration, MultiVINE-seq. The protocol provides a robust, reproducible workflow for the enrichment and high-resolution profiling of vascular, perivascular, and immune cells from fresh or frozen human and mouse brain tissue. First, intact vessels (predominantly capillaries and small arterioles/venules, 100 µm in diameter) are isolated from homogenized brain tissue via dextran-based density-gradient centrifugation, separating the vascular pellet from myelin and the parenchymal fraction. Second, the collected vessels are rigorously washed over a cell strainer to remove trapped contaminants. A critical innovation lies in the third stage: the optimized extraction of nuclei from purified vessels using enzymatic digestion. After extraction, the protocol uses fluorescence-activated cell sorting (FACS) to ensure collection of high-purity nuclei suitable for widely used droplet-based sequencing platforms (e.g., 10x Genomics single cell 3' or multiome). This protocol requires 4-5 h to complete and can be carried out by researchers with single-cell and flow cytometry training.

Journal Article