Nuclear extraction from endometrial tumors for single nuclei sequencing v1

Modified from: Slyper, M. et al. A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors. Nat. Med. 2020 26526, 792–802 (2020). A optimized protocol for nuclear extraction from endometrial tumors. Performed with endometrial adenocarcinoma, endometrioid type, FIGO grade 1.

Bulk and Single-nucleus Transcriptomics Highlight Intra-telencephalic and Somatostatin Neurons in Alzheimer's Disease

Background: Cortical neuron loss is a pathological hallmark of late-onset Alzheimer's disease (AD). However, it remains unclear which neuronal subtypes are most vulnerable to degeneration and contribute most to cognitive decline. Methods: We analyzed postmortem bulk brain RNA-sequencing (RNAseq) data collected from three studies of aging and AD comprising six neocortical regions (704 individuals; 1037 samples). We estimated relative cell type proportions from each brain sample using neuronal subclass-specific marker genes derived from ultra-high depth single-nucleus RNAseq data (snRNAseq). We associated cell type proportions with AD across all samples using mixed-effects mega-analyses. Bulk tissue analyses were complemented by analyses of three AD snRNAseq datasets using the same cell type definitions and diagnostic criteria (51 individuals). Lastly, we identified cell subtype associations with specific neuropathologies, cognitive decline, and residual cognition. Results: In our mega-analyses, we identified the strongest associations of AD with fewer somatostatin (SST) inhibitory neurons (β=-0.48, pbonf=8.98x10-9) and intra-telencephalic (IT) excitatory neurons (β=-0.45, pbonf =4.32x10-7). snRNAseq-based cell type proportion analyses especially supported the association of SST neurons. Analyses of cell type proportions with specific AD-related phenotypes in ROS/MAP consistently implicated fewer SST neurons with greater brain-wide postmortem tau and beta amyloid (β=-0.155, pFDR=3.1x10-4) deposition, as well as more severe cognitive decline prior to death (β=0.309, pFDR=3.9x10-6). Greater IT neuron proportions were associated strongly with improved cognition (β=0.173, pFDR=8.3x10-5) and residual cognition (β=0.175, pFDR=1.2x10-5), but not canonical AD neuropathology. Conclusions: Proportionally fewer SST and IT neurons were significantly associated with AD diagnosis across multiple studies and cortical regions. These findings support seminal work implicating somatostatin and pyramidal neurons in the pathogenesis of AD and improves our current understanding of neuronal vulnerability in AD.

Single-nucleus RNA-seq of normal-appearing brain regions in relapsing-remitting vs. secondary progressive multiple sclerosis

Multiple sclerosis (MS) is an immune-mediated demyelinating disease that alters central nervous system (CNS) functions. Relapsing-remitting MS (RRMS) is the most common form, which can transform into secondary-progressive MS (SPMS) that is associated with progressive neurodegeneration. Single-nucleus RNA sequencing (snRNA-seq) of MS lesions identified disease-related transcriptomic alterations; however, their relationship to non-lesioned MS brain regions has not been reported and which could identify prodromal or other disease susceptibility signatures. Here, snRNA-seq was used to generate high-quality RRMS vs. SPMS datasets of 33,197 nuclei from 8 normal-appearing MS brains, which revealed divergent cell type-specific changes. Notably, SPMS brains downregulated astrocytic sphingosine kinases (SPHK1/2), the enzymes required to phosphorylate and activate the MS drug, fingolimod. This reduction was modeled with astrocyte-specific Sphk1/2 null mice in which fingolimod lost activity, supporting functionality of observed transcriptomic changes. These data provide an initial resource for studies of single cells from non-lesioned RRMS and SPMS brains.

Neuronal Yin Yang1 in the prefrontal cortex regulates transcriptional and behavioral responses to chronic stress in mice

Although the synaptic alterations associated with the stress-related mood disorder major depression has been well-documented, the underlying transcriptional mechanisms remain poorly understood. Here, we perform complementary bulk nuclei- and single-nucleus transcriptome profiling and map locus-specific chromatin interactions in mouse neocortex to identify the cell type-specific transcriptional changes associated with stress-induced behavioral maladaptation. We find that cortical excitatory neurons, layer 2/3 neurons in particular, are vulnerable to chronic stress and acquire signatures of gene transcription and chromatin structure associated with reduced neuronal activity and expression of Yin Yang 1 (YY1). Selective ablation of YY1 in cortical excitatory neurons enhances stress sensitivity in both male and female mice and alters the expression of stress-associated genes following an abbreviated stress exposure. These findings demonstrate how chronic stress impacts transcription in cortical excitatory neurons and identify YY1 as a regulator of stress-induced maladaptive behavior in mice.

Single Nuclei Transcriptome Reveals Perturbed Brain Vascular Molecules in Alzheimer's Disease

Blood-brain barrier (BBB) dysfunction is well-known in Alzheimer's disease (AD), but the precise molecular changes contributing to its pathophysiology are unclear. To understand the transcriptional changes in brain vascular cells, we performed single nucleus RNA sequencing (snRNAseq) of temporal cortex tissue in 24 AD and control brains resulting in 79,751 nuclei, 4,604 of which formed three distinct vascular clusters characterized as activated pericytes, endothelia and resting pericytes. We identified differentially expressed genes (DEGs) and their enriched pathways in these clusters and detected the most transcriptional changes within activated pericytes. Using our data and a knowledge-based predictive algorithm, we discovered and prioritized molecular interactions between vascular and astrocyte clusters, the main cell types of the gliovascular unit (GVU) of the BBB. Vascular targets predicted to interact with astrocytic ligands have biological functions in signaling, angiogenesis, amyloid β metabolism and cytoskeletal structure. Top astrocytic and vascular interacting molecules include both novel and known AD risk genes such as APOE, APP and ECE1. Our findings provide information on transcriptional changes in predicted vascular-astrocytic partners at the GVU, bringing insights to the molecular mechanisms of BBB breakdown in AD.

Single nucleus RNA-sequencing reveals altered intercellular communication and dendritic cell activation in nonobstructive hypertrophic cardiomyopathy

End stage, nonobstructive hypertrophic cardiomyopathy (HCM) is an intractable condition with no disease-specific therapies. To gain insights into the pathogenesis of nonobstructive HCM, we performed single nucleus RNA-sequencing (snRNA-seq) on human HCM hearts explanted at the time of cardiac transplantation and organ donor hearts serving as controls. Differential gene expression analysis revealed 64 differentially expressed genes linked to specific cell types and molecular functions. Analysis of ligand-receptor pair gene expression to delineate potential intercellular communication revealed significant reductions in expressed ligand-receptor pairs affecting the extracellular matrix, growth factor binding, peptidase regulator activity, platelet-derived growth factor binding and protease binding in the HCM tissue. Changes in Integrin-beta1 receptor expression were responsible for many changes related to extracellular matrix interactions, by increasing in dendritic, smooth muscle and pericyte cells while decreasing in endothelial and fibroblast cells, suggesting potential mechanisms for fibrosis and microvascular disease in HCM and a potential role for dendritic cells. In contrast, there was an increase in ligand-receptor pair expression associated with adenylate cyclase binding, calcium channel molecular functions, channel inhibitor activity, ion channel inhibitor activity, phosphatase activator activity, protein kinase activator activity and titin binding, suggesting important shifts in various signaling cascades in nonobstructive, end stage HCM.

Transcriptomic Profiling Reveals Cancer-Associated Fibroblasts as Potential Targets for the Prognosis and Treatment of Cervical Squamous Cell Carcinoma

To understand the etiological, structural, and immunological characteristics of cervical squamous cell carcinoma (CSCC), we conducted single nucleus RNA sequencing (snRNA-seq) and spatial transcriptomics (ST) experiments for cervical samples from 20 individuals. When exploring the possible factors shaping the intra-individual immune heterogeneity in CSCC, we identified a cluster of cancer-associated fibroblasts (CAFs) enriched around some tumors, which highly expressed ACTA2, POSTN, ITGB4, and FAP. Results showed that the CAFs might support the growth and metastasis of tumors by inhibiting lymphocyte infiltration and remodeling the tumor extracellular matrix. Moreover, high CAF signals predicted poorer clinical outcomes for CSCC patients. Our data also revealed the infection profiles of HPV in tumors, the critical factors involved in the progression of cervical cancerous lesions, and the association between tumor metabolism and immune response intensity. Collectively, our findings may improve the prognostic and therapeutic methods for CSCC.

Cellular development and evolution of the mammalian cerebellum

The expansion of the neocortex, one of the hallmarks of mammalian evolution, was accompanied by an increase in the number of cerebellar neurons. However, little is known about the evolution of the cellular programs underlying cerebellum development in mammals. In this study, we generated single-nucleus RNA-sequencing data for ~400,000 cells to trace the development of the cerebellum from early neurogenesis to adulthood in human, mouse, and the marsupial opossum. Our cross-species analyses revealed that the cellular composition and differentiation dynamics throughout cerebellum development are largely conserved, except for human Purkinje cells. Global transcriptome profiles, conserved cell state markers, and gene expression trajectories across neuronal differentiation show that the cerebellar cell type-defining programs have been overall preserved for at least 160 million years. However, we also discovered differences. We identified 3,586 genes that either gained or lost expression in cerebellar cells in one of the species, and 541 genes that evolved new expression trajectories during neuronal differentiation. The potential functional relevance of these cross-species differences is highlighted by the diverged expression patterns of several human disease-associated genes. Altogether, our study reveals shared and lineage-specific programs governing the cellular development of the mammalian cerebellum, and expands our understanding of the evolution of mammalian organ development.

Mapping pediatric brain tumors to their origins in the developing cerebellum

Understanding the cellular origins of childhood brain tumors is key for discovering novel tumor-specific therapeutic targets. Previous strategies mapping cellular origins typically involved comparing human tumors to murine embryonal tissues, a potentially imperfect approach due to spatio-temporal gene expression differences between species. Here we use an unprecedented single-nucleus atlas of the developing human cerebellum (Sepp, Leiss, et al) and extensive bulk and single-cell transcriptome tumor data to map their cellular origins with focus on three most common pediatric brain tumors - pilocytic astrocytoma, ependymoma, and medulloblastoma. Using custom bioinformatics approaches, we postulate the astroglial and glial lineages as the origins for posterior fossa ependymomas and radiation-induced gliomas (secondary tumors after medulloblastoma treatment), respectively. Moreover, we confirm that SHH, Group3 and Group4 medulloblastomas stem from granule cell/unipolar brush cell lineages, whereas we propose pilocytic astrocytoma to originate from the oligodendrocyte lineage. We also identify genes shared between the cerebellar lineage of origin and corresponding tumors, and genes that are tumor specific; both gene sets represent promising therapeutic targets. As a common feature among most cerebellar tumors, we observed compositional heterogeneity in terms of similarity to normal cells, suggesting that tumors arise from or differentiate into multiple points along the cerebellar "lineage of origin".