Single-nucleus and single-cell transcriptomes compared in matched cortical cell types

The use of single-cell RNA-seq to study heterogeneity at varying levels of virus-host interactions

The outcome of viral infection depends on the diversity of the infecting viral population and the heterogeneity of the cell population that is infected. Until almost a decade ago, the study of these dynamic processes during viral infection was challenging and limited to certain targeted measurements. Presently, with the use of single-cell sequencing technology, the complex interface defined by the interactions of cells with infecting virus can now be studied across the breadth of the transcriptome in thousands of individual cells simultaneously. In this review, we will describe the use of single-cell RNA sequencing (scRNA-seq) to study the heterogeneity of viral infections, ranging from individual virions to the immune response between infected individuals. In addition, we highlight certain key experimental limitations and methodological decisions that are critical to analyzing scRNA-seq data at each scale.

Single-cell RNA sequencing uncovers the cell type-dependent transcriptomic changes in the retrosplenial cortex after peripheral nerve injury

The retrosplenial cortex (RSC) is a vital area for storing remote memory and has recently been found to undergo broad changes after peripheral nerve injury. However, little is known about the role of RSC in pain regulation. Here, we examine the involvement of RSC in the pain of mice with nerve injury. Notably, reducing the activities of calcium-/calmodulin-dependent protein kinase type II-positive splenial neurons chemogenetically increases paw withdrawal threshold and extends thermal withdrawal latency in mice with nerve injury. The single-cell or single-nucleus RNA-sequencing results predict enhanced excitatory synaptic transmissions in RSC induced by nerve injury. Local infusion of 1-naphthyl acetyl spermine into RSC to decrease the excitatory synaptic transmissions relieves pain and induces conditioned place preference. Our data indicate that RSC is critical for regulating physiological and neuropathic pain. The cell type-dependent transcriptomic information would help understand the molecular basis of neuropathic pain.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single-cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding disease pathologies. However, several experimental and computational challenges impede transcriptomics-based deconvolution approaches using single-cell/nucleus RNA-seq reference atlases. Cells from the brain and blood have substantially different sizes, total mRNA, and transcriptional activities, and existing approaches may quantify total mRNA instead of cell type proportions. Further, standards are lacking for the use of cell reference atlases and integrative analyses of single-cell and spatial transcriptomics data. We discuss how to approach these key challenges with orthogonal "gold standard" datasets for evaluating deconvolution methods.

Single cell multi-omics reveal intra-cell-line heterogeneity across human cancer cell lines

Human cancer cell lines have long served as tools for cancer research and drug discovery, but the presence and the source of intra-cell-line heterogeneity remain elusive. Here, we perform single-cell RNA-sequencing and ATAC-sequencing on 42 and 39 human cell lines, respectively, to illustrate both transcriptomic and epigenetic heterogeneity within individual cell lines. Our data reveal that transcriptomic heterogeneity is frequently observed in cancer cell lines of different tissue origins, often driven by multiple common transcriptional programs. Copy number variation, as well as epigenetic variation and extrachromosomal DNA distribution all contribute to the detected intra-cell-line heterogeneity. Using hypoxia treatment as an example, we demonstrate that transcriptomic heterogeneity could be reshaped by environmental stress. Overall, our study performs single-cell multi-omics of commonly used human cancer cell lines and offers mechanistic insights into the intra-cell-line heterogeneity and its dynamics, which would serve as an important resource for future cancer cell line-based studies.

Single-cell analysis of innate spinal cord regeneration identifies intersecting modes of neuronal repair

Adult zebrafish have an innate ability to recover from severe spinal cord injury. Here, we report a comprehensive single nuclear RNA sequencing atlas that spans 6 weeks of regeneration. We identify cooperative roles for adult neurogenesis and neuronal plasticity during spinal cord repair. Neurogenesis of glutamatergic and GABAergic neurons restores the excitatory/inhibitory balance after injury. In addition, transient populations of injury-responsive neurons (iNeurons) show elevated plasticity between 1 and 3 weeks post-injury. Using cross-species transcriptomics and CRISPR/Cas9 mutagenesis, we found iNeurons are injury-surviving neurons that share transcriptional similarities with a rare population of spontaneously plastic mouse neurons. iNeurons are required for functional recovery and employ vesicular trafficking as an essential mechanism that underlies neuronal plasticity. This study provides a comprehensive resource of the cells and mechanisms that direct spinal cord regeneration and establishes zebrafish as a model of plasticity-driven neural repair.

Comparative Analysis of Single-cell and Single-nucleus RNA-sequencing in a Rabbit Model of Retinal Detachment-related Proliferative Vitreoretinopathy

Purpose

Proliferative vitreoretinopathy (PVR) is the most common cause of failure of retinal reattachment surgery, and the molecular changes leading to this aberrant wound healing process are currently unknown. Our ultimate goal is to study PVR pathogenesis by employing single-cell transcriptomics to dissect cellular heterogeneity.

Design

Here we aimed to compare single-cell RNA sequencing (scRNA-seq)  and single-nucleus RNA-sequencing (snRNA-seq) of retinal PVR samples in the rabbit model.

Participants

Unilateral induction of PVR lesions in rabbit eyes with contralateral eyes serving as controls.

Methods

Proliferative vitreoretinopathy was induced unilaterally in Dutch Belted rabbits. At different timepoints after PVR induction, retinas were dissociated into either cells or nuclei suspension and processed for scRNA-seq or snRNA-seq.

Main Outcome Measures

Single cell and nuclei transcriptomic profiles of retinas after PVR induction.

Results

Single-cell RNA sequencing and snRNA-seq were conducted on retinas at 4 hours and 14 days after disease induction. Although the capture rate of unique molecular identifiers and genes were greater in scRNA-seq samples, overall gene expression profiles of individual cell types were highly correlated between scRNA-seq and snRNA-seq. A major disparity between the 2 sequencing modalities was the cell type capture rate, however, with glial cell types overrepresented in scRNA-seq, and inner retinal neurons were enriched by snRNA-seq. Furthermore, fibrotic Müller glia were overrepresented in snRNA-seq samples, whereas reactive Müller glia were overrepresented in scRNA-seq samples. Trajectory analyses were similar between the 2 methods, allowing for the combined analysis of the scRNA-seq and snRNA-seq data sets.

Conclusions

These findings highlight limitations of both scRNA-seq and snRNA-seq analysis and imply that use of both techniques together can more accurately identify transcriptional networks critical for aberrant fibrogenesis in PVR than using either in isolation.

Financial Disclosures

Proprietary or commercial disclosure may be found after the references.

Cochlear transcriptome analysis of an outbred mouse population (CFW)

Age-related hearing loss (ARHL) is the most common cause of hearing loss and one of the most prevalent conditions affecting the elderly worldwide. Despite evidence from our lab and others about its polygenic nature, little is known about the specific genes, cell types, and pathways involved in ARHL, impeding the development of therapeutic interventions. In this manuscript, we describe, for the first time, the complete cell-type specific transcriptome of the aging mouse cochlea using snRNA-seq in an outbred mouse model in relation to auditory threshold variation. Cochlear cell types were identified using unsupervised clustering and annotated via a three-tiered approach-first by linking to expression of known marker genes, then using the NSForest algorithm to select minimum cluster-specific marker genes and reduce dimensional feature space for statistical comparison of our clusters with existing publicly-available data sets on the gEAR website, and finally, by validating and refining the annotations using Multiplexed Error Robust Fluorescence In Situ Hybridization (MERFISH) and the cluster-specific marker genes as probes. We report on 60 unique cell-types expanding the number of defined cochlear cell types by more than two times. Importantly, we show significant specific cell type increases and decreases associated with loss of hearing acuity implicating specific subsets of hair cell subtypes, ganglion cell subtypes, and cell subtypes within the stria vascularis in this model of ARHL. These results provide a view into the cellular and molecular mechanisms responsible for age-related hearing loss and pathways for therapeutic targeting.

Anatomical Diversity of the Adult Corticospinal Tract Revealed by Single-Cell Transcriptional Profiling

The corticospinal tract (CST) forms a central part of the voluntary motor apparatus in all mammals. Thus, injury, disease, and subsequent degeneration within this pathway result in chronic irreversible functional deficits. Current strategies to repair the damaged CST are suboptimal in part because of underexplored molecular heterogeneity within the adult tract. Here, we combine spinal retrograde CST tracing with single-cell RNA sequencing (scRNAseq) in adult male and female mice to index corticospinal neuron (CSN) subtypes that differentially innervate the forelimb and hindlimb. We exploit publicly available datasets to confer anatomic specialization among CSNs and show that CSNs segregate not only along the forelimb and hindlimb axis but also by supraspinal axon collateralization. These anatomically defined transcriptional data allow us to use machine learning tools to build classifiers that discriminate between CSNs and cortical layer 2/3 and nonspinally terminating layer 5 neurons in M1 and separately identify limb-specific CSNs. Using these tools, CSN subtypes can be differentially identified to study postnatal patterning of the CST in vivo, leveraged to screen for novel limb-specific axon growth survival and growth activators in vitro, and ultimately exploited to repair the damaged CST after injury and disease.SIGNIFICANCE STATEMENT Therapeutic interventions designed to repair the damaged CST after spinal cord injury have remained functionally suboptimal in part because of an incomplete understanding of the molecular heterogeneity among subclasses of CSNs. Here, we combine spinal retrograde labeling with scRNAseq and annotate a CSN index by the termination pattern of their primary axon in the cervical or lumbar spinal cord and supraspinal collateral terminal fields. Using machine learning we have confirmed the veracity of our CSN gene lists to train classifiers to identify CSNs among all classes of neurons in primary motor cortex to study the development, patterning, homeostasis, and response to injury and disease, and ultimately target streamlined repair strategies to this critical motor pathway.

Cytokine enrichment in deep cerebellar nuclei is contributed by multiple glial populations and linked to reduced amyloid plaque pathology

Alzheimer's disease (AD) pathology and amyloid-beta (Aβ) plaque deposition progress slowly in the cerebellum compared to other brain regions, while the entorhinal cortex (EC) is one of the most vulnerable regions. Using a knock-in AD mouse model (App KI), we show that within the cerebellum, the deep cerebellar nuclei (DCN) has particularly low accumulation of Aβ plaques. To identify factors that might underlie differences in the progression of AD-associated neuropathology across regions, we profiled gene expression in single nuclei (snRNAseq) across all cell types in the DCN and EC of wild-type (WT) and App KI male mice at age 7 months. We found differences in expression of genes associated with inflammatory activation, PI3K-AKT signalling, and neuron support functions between both regions and genotypes. In WT mice, the expression of interferon-response genes in microglia is higher in the DCN than the EC and this enrichment is confirmed by RNA in situ hybridisation, and measurement of inflammatory cytokines by protein array. Our analyses also revealed that multiple glial populations are responsible for establishing this cytokine-enriched niche. Furthermore, homogenates derived from the DCN induced inflammatory gene expression in BV2 microglia. We also assessed the relationship between the DCN microenvironment and Aβ pathology by depleting microglia using a CSF1R inhibitor PLX5622 and saw that, surprisingly, the expression of a subset of inflammatory cytokines was increased while plaque abundance in the DCN was further reduced. Overall, our study revealed the presence of a cytokine-enriched microenvironment unique to the DCN that when modulated, can alter plaque deposition.

Single-Nuclei Multiome (ATAC + Gene Expression) Sequencing of a Primary Canine Osteosarcoma Elucidates Intra-Tumoral Heterogeneity and Characterizes the Tumor Microenvironment

Osteosarcoma (OSA) is a highly aggressive bone tumor primarily affecting pediatric or adolescent humans and large-breed dogs. Canine OSA shares striking similarities with its human counterpart, making it an invaluable translational model for uncovering the disease's complexities and developing novel therapeutic strategies. Tumor heterogeneity, a hallmark of OSA, poses significant challenges to effective treatment due to the evolution of diverse cell populations that influence tumor growth, metastasis, and resistance to therapies. In this study, we apply single-nuclei multiome sequencing, encompassing ATAC (Assay for Transposase-Accessible Chromatin) and GEX (Gene Expression, or RNA) sequencing, to a treatment-naïve primary canine osteosarcoma. This comprehensive approach reveals the complexity of the tumor microenvironment by simultaneously capturing the transcriptomic and epigenomic profiles within the same nucleus. Furthermore, these results are analyzed in conjunction with bulk RNA sequencing and differential analysis of the same tumor and patient-matched normal bone. By delving into the intricacies of OSA at this unprecedented level of detail, we aim to unravel the underlying mechanisms driving intra-tumoral heterogeneity, opening new avenues for therapeutic interventions in both human and canine patients. This study pioneers an approach that is broadly applicable, while demonstrating significant heterogeneity in the context of a single individual's tumor.

Enzalutamide Induces Cytotoxicity in Desmoplastic Small Round Cell Tumor Independent of the Androgen Receptor

Desmoplastic Small Round Cell Tumor (DSRCT) is a rare, pediatric cancer caused by the EWSR1::WT1 fusion protein. DSRCT predominantly occurs in males, which comprise 80-90% of the patient population. While the reason for this male predominance remains unknown, one hypothesis is that the androgen receptor (AR) plays a critical role in DSRCT and elevated testosterone levels in males help drive tumor growth. Here, we demonstrate that AR is highly expressed in DSRCT relative to other fusion-driven sarcomas and that the AR antagonists enzalutamide and flutamide reduce DSRCT growth. However, despite these findings, which suggest an important role for AR in DSRCT, we show that DSRCT cell lines form xenografts in female mice at the same rate as male mice and AR depletion does not significantly alter DSRCT growth in vitro. Further, we find that AR antagonists reduce DSRCT growth in cells depleted of AR, establishing an AR-independent mechanism of action. These findings suggest that AR dependence is not the reason for male predominance in DSRCT and that AR-targeted therapies may provide therapeutic benefit primarily through an AR-independent mechanism that requires further elucidation.

Functional genomics and systems biology in human neuroscience

Neuroscience research has entered a phase of key discoveries in the realm of neurogenomics owing to strong financial and intellectual support for resource building and tool development. The previous challenge of tissue heterogeneity has been met with the application of techniques that can profile individual cells at scale. Moreover, the ability to perturb genes, gene regulatory elements and neuronal activity in a cell-type-specific manner has been integrated with gene expression studies to uncover the functional underpinnings of the genome at a systems level. Although these insights have necessarily been grounded in model systems, we now have the opportunity to apply these approaches in humans and in human tissue, thanks to advances in human genetics, brain imaging and tissue collection. We acknowledge that there will probably always be limits to the extent to which we can apply the genomic tools developed in model systems to human neuroscience; however, as we describe in this Perspective, the neuroscience field is now primed with an optimal foundation for tackling this ambitious challenge. The application of systems-level network analyses to these datasets will facilitate a deeper appreciation of human neurogenomics that cannot otherwise be achieved from directly observable phenomena.

Single-Nucleus RNA-Sequencing in Brain Tissue

Neuroscience research greatly benefits from single-cell sequencing technologies, which can reveal transcriptional alterations on a cellular level. However, preparing single-cell suspensions is technically challenging, requires experience, and has several limitations that can influence the transcriptional readout. Performing sequencing of single nuclei instead of single cells alleviates several of the challenges of sample preparation and highlights acute nuclear transcription. Here, we provide a protocol to prepare a nuclei suspension for single-nucleus RNA-sequencing for cell type-specific transcriptional profiling of brain tissue using the 10x Genomics single-cell gene expression assay. Furthermore, we highlight important aspects to consider during experimental design and data analysis. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Preparation of single-nucleus suspension Basic Protocol 2: Preparation and sequencing of single-nucleus libraries for RNA-seq.

Single-nucleus transcriptome inventory of giant panda reveals cellular basis for fitness optimization under low metabolism

BACKGROUND

Energy homeostasis is essential for the adaptation of animals to their environment and some wild animals keep low metabolism adaptive to their low-nutrient dietary supply. Giant panda is such a typical low-metabolic mammal exhibiting species specialization of extremely low daily energy expenditure. It has low levels of basal metabolic rate, thyroid hormone, and physical activities, whereas the cellular bases of its low metabolic adaptation remain rarely explored.

RESULTS

In this study, we generate a single-nucleus transcriptome atlas of 21 organs/tissues from a female giant panda. We focused on the central metabolic organ (liver) and dissected cellular metabolic status by cross-species comparison. Adaptive expression mode (i.e., AMPK related) was prominently displayed in the hepatocyte of giant panda. In the highest energy-consuming organ, the heart, we found a possibly optimized utilization of fatty acid. Detailed cell subtype annotation of endothelial cells showed the uterine-specific deficiency of blood vascular subclasses, indicating a potential adaptation for a low reproductive energy expenditure.

CONCLUSIONS

Our findings shed light on the possible cellular basis and transcriptomic regulatory clues for the low metabolism in giant pandas and helped to understand physiological adaptation response to nutrient stress.

Bibliometric analysis of research on neurodegenerative diseases and single-cell RNA sequencing: Opportunities and challenges

Neurodegeneration, characterized by the progressive deterioration in neuronal structure or function, presents an elusive mechanism. The use of single-cell RNA sequencing (scRNA-seq) technology in the clinic is becoming increasingly prevalent in recent decades. This technology offers unparalleled cell-level insights into neurodegenerative diseases, establishing itself as a potent tool for elucidating these diseases underlying mechanisms. Here, we made a deep investigation for scRNA-seq research in neurodegenerative diseases using bibliometric analysis from 2009 to 2022. We observed a robust upward trajectory in the number of publications on this subject. The United States stood out as the principal contributor to this expanding field. Specifically, the University of California System exhibited notable research prowess in this field. Alzheimer disease and Parkinson disease were the diseases most frequently investigated. Key research hotspots include the creation of a molecular brain atlas and identification of vulnerable neuronal subpopulations and potential therapeutic targets at the transcriptomic level.

Interindividual variation in human cortical cell type abundance and expression

Single-cell transcriptomic studies have identified a conserved set of neocortical cell types from small postmortem cohorts. We extended these efforts by assessing cell type variation across 75 adult individuals undergoing epilepsy and tumor surgeries. Nearly all nuclei map to one of 125 robust cell types identified in the middle temporal gyrus. However, we found interindividual variance in abundances and gene expression signatures, particularly in deep-layer glutamatergic neurons and microglia. A minority of donor variance is explainable by age, sex, ancestry, disease state, and cell state. Genomic variation was associated with expression of 150 to 250 genes for most cell types. This characterization of cellular variation provides a baseline for cell typing in health and disease.

Neurodevelopmental deficits and cell-type-specific transcriptomic perturbations in a mouse model of HNRNPU haploinsufficiency

Heterozygous de novo loss-of-function mutations in the gene expression regulator HNRNPU cause an early-onset developmental and epileptic encephalopathy. To gain insight into pathological mechanisms and lay the potential groundwork for developing targeted therapies, we characterized the neurophysiologic and cell-type-specific transcriptomic consequences of a mouse model of HNRNPU haploinsufficiency. Heterozygous mutants demonstrated global developmental delay, impaired ultrasonic vocalizations, cognitive dysfunction and increased seizure susceptibility, thus modeling aspects of the human disease. Single-cell RNA-sequencing of hippocampal and neocortical cells revealed widespread, yet modest, dysregulation of gene expression across mutant neuronal subtypes. We observed an increased burden of differentially-expressed genes in mutant excitatory neurons of the subiculum-a region of the hippocampus implicated in temporal lobe epilepsy. Evaluation of transcriptomic signature reversal as a therapeutic strategy highlights the potential importance of generating cell-type-specific signatures. Overall, this work provides insight into HNRNPU-mediated disease mechanisms and provides a framework for using single-cell RNA-sequencing to study transcriptional regulators implicated in disease.

Glyoxal as an alternative fixative for single-cell RNA sequencing

Single-cell RNA sequencing has become an important method to identify cell types, delineate the trajectories of cell differentiation in whole organisms, and understand the heterogeneity in cellular responses. Nevertheless, sample collection and processing remain a severe bottleneck for single-cell RNA sequencing experiments. Cell isolation protocols often lead to significant changes in the transcriptomes of cells, requiring novel methods to preserve cell states. Here, we developed and benchmarked protocols using glyoxal as a fixative for single-cell RNA sequencing applications. Using Drop-seq methodology, we detected a large number of transcripts and genes from glyoxal-fixed Drosophila cells after single-cell RNA sequencing. The effective glyoxal fixation of transcriptomes in Drosophila and human cells was further supported by a high correlation of gene expression data between glyoxal-fixed and unfixed samples. Accordingly, we also found highly expressed genes overlapping to a large extent between experimental conditions. These results indicated that our fixation protocol did not induce considerable changes in gene expression and conserved the transcriptome for subsequent single-cell isolation procedures. In conclusion, we present glyoxal as a suitable fixative for Drosophila cells and potentially cells of other species that allow high-quality single-cell RNA sequencing applications.

Celf4 controls mRNA translation underlying synaptic development in the prenatal mammalian neocortex

Abnormalities in neocortical and synaptic development are linked to neurodevelopmental disorders. However, the molecular and cellular mechanisms governing initial synapse formation in the prenatal neocortex remain poorly understood. Using polysome profiling coupled with snRNAseq on human cortical samples at various fetal phases, we identify human mRNAs, including those encoding synaptic proteins, with finely controlled translation in distinct cell populations of developing frontal neocortices. Examination of murine and human neocortex reveals that the RNA binding protein and translational regulator, CELF4, is expressed in compartments enriched in initial synaptogenesis: the marginal zone and the subplate. We also find that Celf4/CELF4-target mRNAs are encoded by risk genes for adverse neurodevelopmental outcomes translating into synaptic proteins. Surprisingly, deleting Celf4 in the forebrain disrupts the balance of subplate synapses in a sex-specific fashion. This highlights the significance of RNA binding proteins and mRNA translation in evolutionarily advanced synaptic development, potentially contributing to sex differences.

Single-nucleus RNA sequencing of developing superior colliculus identifies neuronal diversity and candidate mediators of circuit assembly

The superior colliculus (SC) is a sensorimotor structure in the midbrain that integrates input from multiple sensory modalities to initiate motor commands. It undergoes well-characterized steps of circuit assembly during development, rendering the mouse SC a popular model to study establishment of neural connectivity. Here we perform single-nucleus RNA-sequencing analysis of the mouse SC isolated at various developmental time points. Our study provides a transcriptomic landscape of the cell types that comprise the SC across murine development with particular emphasis on neuronal heterogeneity. We report a repertoire of genes differentially expressed across the different postnatal ages, many of which are known to regulate axon guidance and synapse formation. Using these data, we find that Pax7 expression is restricted to a subset of GABAergic neurons. Our data provide a valuable resource for interrogating the mechanisms of circuit development and identifying markers for manipulating specific SC neuronal populations and circuits.

Application of single-cell multi-omics approaches in horticulture research

Cell heterogeneity shapes the morphology and function of various tissues and organs in multicellular organisms. Elucidation of the differences among cells and the mechanism of intercellular regulation is essential for an in-depth understanding of the developmental process. In recent years, the rapid development of high-throughput single-cell transcriptome sequencing technologies has influenced the study of plant developmental biology. Additionally, the accuracy and sensitivity of tools used to study the epigenome and metabolome have significantly increased, thus enabling multi-omics analysis at single-cell resolution. Here, we summarize the currently available single-cell multi-omics approaches and their recent applications in plant research, review the single-cell based studies in fruit, vegetable, and ornamental crops, and discuss the potential of such approaches in future horticulture research.

A scATAC-seq atlas of chromatin accessibility in axolotl brain regions

Axolotl (Ambystoma mexicanum) is an excellent model for investigating regeneration, the interaction between regenerative and developmental processes, comparative genomics, and evolution. The brain, which serves as the material basis of consciousness, learning, memory, and behavior, is the most complex and advanced organ in axolotl. The modulation of transcription factors is a crucial aspect in determining the function of diverse regions within the brain. There is, however, no comprehensive understanding of the gene regulatory network of axolotl brain regions. Here, we utilized single-cell ATAC sequencing to generate the chromatin accessibility landscapes of 81,199 cells from the olfactory bulb, telencephalon, diencephalon and mesencephalon, hypothalamus and pituitary, and the rhombencephalon. Based on these data, we identified key transcription factors specific to distinct cell types and compared cell type functions across brain regions. Our results provide a foundation for comprehensive analysis of gene regulatory programs, which are valuable for future studies of axolotl brain development, regeneration, and evolution, as well as on the mechanisms underlying cell-type diversity in vertebrate brains.

Single-Nucleus RNA-Seq: Open the Era of Great Navigation for FFPE Tissue

Single-cell sequencing (scRNA-seq) has revolutionized our ability to explore heterogeneity and genetic variations at the single-cell level, opening up new avenues for understanding disease mechanisms and cell-cell interactions. Single-nucleus RNA-sequencing (snRNA-seq) is emerging as a promising solution to scRNA-seq due to its reduced ionized transcription bias and compatibility with richer samples. This approach will provide an exciting opportunity for in-depth exploration of billions of formalin-fixed paraffin-embedded (FFPE) tissues. Recent advancements in single-cell/nucleus gene expression workflows tailored for FFPE tissues have demonstrated their feasibility and provided crucial guidance for future studies utilizing FFPE specimens. In this review, we provide a broad overview of the nuclear preparation strategies, the latest technologies of snRNA-seq applicable to FFPE samples. Finally, the limitations and potential technical developments of snRNA-seq in FFPE samples are summarized. The development of snRNA-seq technologies for FFPE samples will lay a foundation for transcriptomic studies of valuable samples in clinical medicine and human sample banks.

Cortical somatostatin interneuron subtypes form cell-type-specific circuits

The cardinal classes are a useful simplification of cortical interneuron diversity, but such broad subgroupings gloss over the molecular, morphological, and circuit specificity of interneuron subtypes, most notably among the somatostatin interneuron class. Although there is evidence that this diversity is functionally relevant, the circuit implications of this diversity are unknown. To address this knowledge gap, we designed a series of genetic strategies to target the breadth of somatostatin interneuron subtypes and found that each subtype possesses a unique laminar organization and stereotyped axonal projection pattern. Using these strategies, we examined the afferent and efferent connectivity of three subtypes (two Martinotti and one non-Martinotti) and demonstrated that they possess selective connectivity with intratelecephalic or pyramidal tract neurons. Even when two subtypes targeted the same pyramidal cell type, their synaptic targeting proved selective for particular dendritic compartments. We thus provide evidence that subtypes of somatostatin interneurons form cell-type-specific cortical circuits.

Heterogeneous tumor microenvironment in pancreatic ductal adenocarcinoma: An emerging role of single-cell analysis

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is one of the deadliest malignancies in the world, for which the mortality is almost as high as the disease incidence and is predicted to be the second-highest cause of cancer-related deaths by 2030. These cancerous tumors consist of diversified gene expressions within the different cellular subpopulations that include neoplastic ductal cells, cancer-associated fibroblasts, and immune cells, all of which collectively facilitate cellular heterogeneity in the PDAC tumor microenvironment (TME). Active intratumoral interaction within the cell populations in TME induces the proliferation of cancerous cells, accounting for tumorigenesis and rapid metastasis.

METHODS

This review will focus on novel findings uncovering PDAC heterogeneity in different cellular subpopulations using single-cell RNA-sequencing (scRNA-seq) and other single-cell analysis technologies. It will further explore the emerging role of single-cell technologies in assessing the role of different subpopulations of neoplastic ductal cells, cancer-associated fibroblasts, and immune cells in PDAC progression.

RESULTS AND CONCLUSION

The application of scRNA-seq in PDAC has started to unveil associations between disease progression and heterogeneity in pancreatic TME and could influence future PDAC treatment. Recent advances in scRNA-seq have uncovered comprehensive analyses of heterogeneous ecosystems present within the TME. These emerging findings underpins further need for a more in-depth understanding of intratumoral heterogeneity in the PDAC microenvironment.

Simultaneous profiling of chromatin architecture and transcription in single cells

The three-dimensional structure of chromatin plays a crucial role in development and disease, both of which are associated with transcriptional changes. However, given the heterogeneity in single-cell chromatin architecture and transcription, the regulatory relationship between the three-dimensional chromatin structure and gene expression is difficult to explain based on bulk cell populations. Here we develop a single-cell, multimodal, omics method allowing the simultaneous detection of chromatin architecture and messenger RNA expression by sequencing (single-cell transcriptome sequencing (scCARE-seq)). Applying scCARE-seq to examine chromatin architecture and transcription from 2i to serum single mouse embryonic stem cells, we observe improved separation of cell clusters compared with single-cell chromatin conformation capture. In addition, after defining the cell-cycle phase of each cell through chromatin architecture extracted by scCARE-seq, we find that periodic changes in chromatin architecture occur in parallel with transcription during the cell cycle. These findings highlight the potential of scCARE-seq to facilitate comprehensive analyses that may boost our understanding of chromatin architecture and transcription in the same single cell.

Single cell genomics as a transformative approach for aquaculture research and innovation

Single cell genomics encompasses a suite of rapidly maturing technologies that measure the molecular profiles of individual cells within target samples. These approaches provide a large up-step in biological information compared to long-established 'bulk' methods that profile the average molecular profiles of all cells in a sample, and have led to transformative advances in understanding of cellular biology, particularly in humans and model organisms. The application of single cell genomics is fast expanding to non-model taxa, including aquaculture species, where numerous research applications are underway with many more envisaged. In this review, we highlight the potential transformative applications of single cell genomics in aquaculture research, considering barriers and potential solutions to the broad uptake of these technologies. Focusing on single cell transcriptomics, we outline considerations for experimental design, including the essential requirement to obtain high quality cells/nuclei for sequencing in ectothermic aquatic species. We further outline data analysis and bioinformatics considerations, tailored to studies with the under-characterized genomes of aquaculture species, where our knowledge of cellular heterogeneity and cell marker genes is immature. Overall, this review offers a useful source of knowledge for researchers aiming to apply single cell genomics to address biological challenges faced by the global aquaculture sector though an improved understanding of cell biology.

Cell-specific MAPT gene expression is preserved in neuronal and glial tau cytopathologies in progressive supranuclear palsy

Microtubule-associated protein tau (MAPT) aggregates in neurons, astrocytes and oligodendrocytes in a number of neurodegenerative diseases, including progressive supranuclear palsy (PSP). Tau is a target of therapy and the strategy includes either the elimination of pathological tau aggregates or reducing MAPT expression, and thus the amount of tau protein made to prevent its aggregation. Disease-associated tau affects brain regions in a sequential manner that includes cell-to-cell spreading. Involvement of glial cells that show tau aggregates is interpreted as glial cells taking up misfolded tau assuming that glial cells do not express enough MAPT. Although studies have evaluated MAPT expression in human brain tissue homogenates, it is not clear whether MAPT expression is compromised in cells accumulating pathological tau. To address these perplexing aspects of disease pathogenesis, this study used RNAscope combined with immunofluorescence (AT8), and single-nuclear(sn) RNAseq to systematically map and quantify MAPT expression dynamics across different cell types and brain regions in controls (n = 3) and evaluated whether tau cytopathology affects MAPT expression in PSP (n = 3). MAPT transcripts were detected in neurons, astrocytes and oligodendrocytes, and varied between brain regions and within each cell type, and were preserved in all cell types with tau aggregates in PSP. These results propose a complex scenario in all cell types, where, in addition to the ingested misfolded tau, the preserved cellular MAPT expression provides a pool for local protein production that can (1) be phosphorylated and aggregated, or (2) feed the seeding of ingested misfolded tau by providing physiological tau, both accentuating the pathological process. Since tau cytopathology does not compromise MAPT gene expression in PSP, a complete loss of tau protein expression as an early pathogenic component is less likely. These observations provide rationale for a dual approach to therapy by decreasing cellular MAPT expression and targeting removal of misfolded tau.

Methods and Insights from Single-Cell Expression Quantitative Trait Loci

Recent advancements in single-cell technologies have enabled expression quantitative trait locus (eQTL) analysis across many individuals at single-cell resolution. Compared with bulk RNA sequencing, which averages gene expression across cell types and cell states, single-cell assays capture the transcriptional states of individual cells, including fine-grained, transient, and difficult-to-isolate populations at unprecedented scale and resolution. Single-cell eQTL (sc-eQTL) mapping can identify context-dependent eQTLs that vary with cell states, including some that colocalize with disease variants identified in genome-wide association studies. By uncovering the precise contexts in which these eQTLs act, single-cell approaches can unveil previously hidden regulatory effects and pinpoint important cell states underlying molecular mechanisms of disease. Here, we present an overview of recently deployed experimental designs in sc-eQTL studies. In the process, we consider the influence of study design choices such as cohort, cell states, and ex vivo perturbations. We then discuss current methodologies, modeling approaches, and technical challenges as well as future opportunities and applications.

Single-cell spatial transcriptome reveals cell-type organization in the macaque cortex

Elucidating the cellular organization of the cerebral cortex is critical for understanding brain structure and function. Using large-scale single-nucleus RNA sequencing and spatial transcriptomic analysis of 143 macaque cortical regions, we obtained a comprehensive atlas of 264 transcriptome-defined cortical cell types and mapped their spatial distribution across the entire cortex. We characterized the cortical layer and region preferences of glutamatergic, GABAergic, and non-neuronal cell types, as well as regional differences in cell-type composition and neighborhood complexity. Notably, we discovered a relationship between the regional distribution of various cell types and the region's hierarchical level in the visual and somatosensory systems. Cross-species comparison of transcriptomic data from human, macaque, and mouse cortices further revealed primate-specific cell types that are enriched in layer 4, with their marker genes expressed in a region-dependent manner. Our data provide a cellular and molecular basis for understanding the evolution, development, aging, and pathogenesis of the primate brain.

An integrated organoid omics map extends modeling potential of kidney disease

Kidney organoids are a promising model to study kidney disease, but their use is constrained by limited knowledge of their functional protein expression profile. Here, we define the organoid proteome and transcriptome trajectories over culture duration and upon exposure to TNFα, a cytokine stressor. Older organoids increase deposition of extracellular matrix but decrease expression of glomerular proteins. Single cell transcriptome integration reveals that most proteome changes localize to podocytes, tubular and stromal cells. TNFα treatment of organoids results in 322 differentially expressed proteins, including cytokines and complement components. Transcript expression of these 322 proteins is significantly higher in individuals with poorer clinical outcomes in proteinuric kidney disease. Key TNFα-associated protein (C3 and VCAM1) expression is increased in both human tubular and organoid kidney cell populations, highlighting the potential for organoids to advance biomarker development. By integrating kidney organoid omic layers, incorporating a disease-relevant cytokine stressor and comparing with human data, we provide crucial evidence for the functional relevance of the kidney organoid model to human kidney disease.

Opportunities and challenges in the application of single-cell and spatial transcriptomics in plants

Single-cell and spatial transcriptomics have diverted researchers' attention from the multicellular level to the single-cell level and spatial information. Single-cell transcriptomes provide insights into the transcriptome at the single-cell level, whereas spatial transcriptomes help preserve spatial information. Although these two omics technologies are helpful and mature, further research is needed to ensure their widespread applicability in plant studies. Reviewing recent research on plant single-cell or spatial transcriptomics, we compared the different experimental methods used in various plants. The limitations and challenges are clear for both single-cell and spatial transcriptomic analyses, such as the lack of applicability, spatial information, or high resolution. Subsequently, we put forth further applications, such as cross-species analysis of roots at the single-cell level and the idea that single-cell transcriptome analysis needs to be combined with other omics analyses to achieve superiority over individual omics analyses. Overall, the results of this review suggest that combining single-cell transcriptomics, spatial transcriptomics, and spatial element distribution can provide a promising research direction, particularly for plant research.

Human Alzheimer's disease reactive astrocytes exhibit a loss of homeostastic gene expression

Astrocytes are one of the brain's major cell types and are responsible for maintaining neuronal homeostasis via regulating the extracellular environment, providing metabolic support, and modulating synaptic activity. In neurodegenerative diseases, such as Alzheimer's disease, astrocytes can take on a hypertrophic appearance. These reactive astrocytes are canonically associated with increases in cytoskeletal proteins, such as glial fibrillary acidic protein and vimentin. However, the molecular alterations that characterize astrocytes in human disease tissues have not been extensively studied with single cell resolution. Using single nucleus RNA sequencing data from normal, pathologic aging, and Alzheimer's disease brains, we identified the transcriptomic changes associated with reactive astrocytes. Deep learning-based clustering algorithms denoised expression data for 17,012 genes and clustered 15,529 astrocyte nuclei, identifying protoplasmic, gray matter and fibrous, white matter astrocyte clusters. RNA trajectory analyses revealed a spectrum of reactivity within protoplasmic astrocytes characterized by a modest increase of reactive genes and a marked decrease in homeostatic genes. Amyloid but not tau pathology correlated with astrocyte reactivity. To identify reactivity-associated genes, linear regressions of gene expression versus reactivity were used to identify the top 52 upregulated and 144 downregulated genes. Gene Ontology analysis revealed that upregulated genes were associated with cellular growth, responses to metal ions, inflammation, and proteostasis. Downregulated genes were involved in cellular interactions, neuronal development, ERBB signaling, and synapse regulation. Transcription factors were significantly enriched among the downregulated genes. Using co-immunofluorescence staining of Alzheimer's disease brain tissues, we confirmed pathologic downregulation of ERBB4 and transcription factor NFIA in reactive astrocytes. Our findings reveal that protoplasmic, gray matter astrocytes in Alzheimer's disease exist within a spectrum of reactivity that is marked by a strong loss of normal function.

Cell-type specific transcriptional adaptations of nucleus accumbens interneurons to amphetamine

Parvalbumin-expressing (PV+) interneurons of the nucleus accumbens (NAc) play an essential role in the addictive-like behaviors induced by psychostimulant exposure. To identify molecular mechanisms of PV+ neuron plasticity, we isolated interneuron nuclei from the NAc of male and female mice following acute or repeated exposure to amphetamine (AMPH) and sequenced for cell type-specific RNA expression and chromatin accessibility. AMPH regulated the transcription of hundreds of genes in PV+ interneurons, and this program was largely distinct from that regulated in other NAc GABAergic neurons. Chromatin accessibility at enhancers predicted cell-type specific gene regulation, identifying transcriptional mechanisms of differential AMPH responses. Finally, we assessed expression of PV-enriched, AMPH-regulated genes in an Mecp2 mutant mouse strain that shows heightened behavioral sensitivity to psychostimulants to explore the functional importance of this transcriptional program. Together these data provide novel insight into the cell-type specific programs of transcriptional plasticity in NAc neurons that underlie addictive-like behaviors.

Isolation of Projection-Specific and Behavior-Relevant Amygdala Circuit for RNA Sequencing

One of the most sought-after topics in neuroscience is to understand how the environment regulates the activity and function of neural circuitry and subsequently influences relevant behaviors. In response to alterations in the environment, the neural circuits undergo adaptive changes ranging from gene expression changes to altered cellular function. Performing sequencing of the transcriptome involved in these behavior-related circuits will provide clues to accurately dissect the detailed mechanisms of related behavior. Here, we describe methods for marking and collecting the ventral hippocampus-projecting basolateral amygdala neurons, which have been repeatedly implicated in regulation of anxiety-like behavior, and subsequently constructing a library ready for sequencing. Specifically, the reported approaches include adeno-associated virus injection, acute brain slice isolation, cell suspension preparation, cell extraction, and cDNA library construction. By utilizing the techniques described here, researchers can comprehensively investigate the transcriptional levels of neural clusters embedded in particular circuits and discover potential pathogenic and therapeutic targets for behavior-relevant disorders. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Tagging of behavior-related neural circuits Basic Protocol 2: Isolation and capture of fluorescent-positive cells Basic Protocol 3: Foundation of sequencing library.

Harmonized cross-species cell atlases of trigeminal and dorsal root ganglia

Peripheral sensory neurons in the dorsal root ganglion (DRG) and trigeminal ganglion (TG) are specialized to detect and transduce diverse environmental stimuli including touch, temperature, and pain to the central nervous system. Recent advances in single-cell RNA-sequencing (scRNA-seq) have provided new insights into the diversity of sensory ganglia cell types in rodents, non-human primates, and humans, but it remains difficult to compare transcriptomically defined cell types across studies and species. Here, we built cross-species harmonized atlases of DRG and TG cell types that describe 18 neuronal and 11 non-neuronal cell types across 6 species and 19 studies. We then demonstrate the utility of this harmonized reference atlas by using it to annotate newly profiled DRG nuclei/cells from both human and the highly regenerative axolotl. We observe that the transcriptomic profiles of sensory neuron subtypes are broadly similar across vertebrates, but the expression of functionally important neuropeptides and channels can vary notably. The new resources and data presented here can guide future studies in comparative transcriptomics, simplify cell type nomenclature differences across studies, and help prioritize targets for future pain therapy development.

Single-cell transcriptome atlas of Drosophila gastrula 2.0

During development, positional information directs cells to specific fates, leading them to differentiate with their own transcriptomes and express specific behaviors and functions. However, the mechanisms underlying these processes in a genome-wide view remain ambiguous, partly because the single-cell transcriptomic data of early developing embryos containing accurate spatial and lineage information are still lacking. Here, we report a single-cell transcriptome atlas of Drosophila gastrulae, divided into 77 transcriptomically distinct clusters. We find that the expression profiles of plasma-membrane-related genes, but not those of transcription-factor genes, represent each germ layer, supporting the nonequivalent contribution of each transcription-factor mRNA level to effector gene expression profiles at the transcriptome level. We also reconstruct the spatial expression patterns of all genes at the single-cell stripe level as the smallest unit. This atlas is an important resource for the genome-wide understanding of the mechanisms by which genes cooperatively orchestrate Drosophila gastrulation.

New technologies to study helminth development and host-parasite interactions

How parasites develop and survive, and how they stimulate or modulate host immune responses are important in understanding disease pathology and for the design of new control strategies. Microarray analysis and bulk RNA sequencing have provided a wealth of data on gene expression as parasites develop through different life-cycle stages and on host cell responses to infection. These techniques have enabled gene expression in the whole organism or host tissue to be detailed, but do not take account of the heterogeneity between cells of different types or developmental stages, nor the spatial organisation of these cells. Single-cell RNA-seq (scRNA-seq) adds a new dimension to studying parasite biology and host immunity by enabling gene profiling at the individual cell level. Here we review the application of scRNA-seq to establish gene expression cell atlases for multicellular helminths and to explore the expansion and molecular profile of individual host cell types involved in parasite immunity and tissue repair. Studying host-parasite interactions in vivo is challenging and we conclude this review by briefly discussing the applications of organoids (stem-cell derived mini-tissues) to examine host-parasite interactions at the local level, and as a potential system to study parasite development in vitro. Organoid technology and its applications have developed rapidly, and the elegant studies performed to date support the use of organoids as an alternative in vitro system for research on helminth parasites.

Multicellular communities are perturbed in the aging human brain and Alzheimer's disease

The role of different cell types and their interactions in Alzheimer's disease (AD) is a complex and open question. Here, we pursued this question by assembling a high-resolution cellular map of the aging frontal cortex using single-nucleus RNA sequencing of 24 individuals with a range of clinicopathologic characteristics. We used this map to infer the neocortical cellular architecture of 638 individuals profiled by bulk RNA sequencing, providing the sample size necessary for identifying statistically robust associations. We uncovered diverse cell populations associated with AD, including a somatostatin inhibitory neuronal subtype and oligodendroglial states. We further identified a network of multicellular communities, each composed of coordinated subpopulations of neuronal, glial and endothelial cells, and we found that two of these communities are altered in AD. Finally, we used mediation analyses to prioritize cellular changes that might contribute to cognitive decline. Thus, our deconstruction of the aging neocortex provides a roadmap for evaluating the cellular microenvironments underlying AD and dementia.

A single-cell level comparison of human inner ear organoids with the human cochlea and vestibular organs

Inner ear disorders are among the most common congenital abnormalities; however, current tissue culture models lack the cell type diversity to study these disorders and normal otic development. Here, we demonstrate the robustness of human pluripotent stem cell-derived inner ear organoids (IEOs) and evaluate cell type heterogeneity by single-cell transcriptomics. To validate our findings, we construct a single-cell atlas of human fetal and adult inner ear tissue. Our study identifies various cell types in the IEOs including periotic mesenchyme, type I and type II vestibular hair cells, and developing vestibular and cochlear epithelium. Many genes linked to congenital inner ear dysfunction are confirmed to be expressed in these cell types. Additional cell-cell communication analysis within IEOs and fetal tissue highlights the role of endothelial cells on the developing sensory epithelium. These findings provide insights into this organoid model and its potential applications in studying inner ear development and disorders.

Ambient RNAs removal of cortex-specific snRNA-seq reveals Apoe+ microglia/macrophage after deeper cerebral hypoperfusion in mice

BACKGROUND

Ambient RNAs contamination in single-nuclei RNA sequencing (snRNA-seq) is a challenging problem, but the consequences of ambient RNAs contamination of damaged and/or diseased tissues are poorly understood. Cognitive impairments and white/gray matter injuries are characteristic of deeper cerebral hypoperfusion mouse models induced by bilateral carotid artery stenosis (BCAS), but the molecular mechanisms still need to be further explored. More importantly, the BCAS mice can also offer an excellent model to examine the signatures of ambient RNAs contamination in damaged tissues when performing snRNA-seq.

METHODS

After the sham and BCAS mice were established, cortex-specific single-nuclei libraries were constructed. Single-nuclei transcriptomes were described informatically by the R package Seurat, and ambient RNA markers of were identified in each library. Then, after removing ambient RNAs in each sample using the in silico approaches, the combination of CellBender and subcluster cleaning, single-nuclei transcriptomes were reconstructed. Next, the comparison of ambient RNA contamination was performed using irGSEA analysis before and after the in silico approaches. Finally, further bioinformatic analyses were performed.

RESULTS

The ambient RNAs are more predominant in the BCAS group than the sham group. The contamination mainly originated from damaged neuronal nuclei, but could be reduced largely using the in silico approaches. The integrative analysis of cortex-specific snRNA-seq data and the published bulk transcriptome revealed that microglia and other immune cells were the primary effectors. In the sequential microglia/immune subgroups analysis, the subgroup of Apoe⁺ MG/Mac (microglia/macrophages) was identified. Interestingly, this subgroup mainly participated in the pathways of lipid metabolism, associated with the phagocytosis of cell debris.

CONCLUSIONS

Taken together, our current study unravels the features of ambient RNAs in snRNA-seq datasets under diseased conditions, and the in silico approaches can effectively eliminate the incorrected cell annotation and following misleading analysis. In the future, snRNA-seq data analysis should be carefully revisited, and ambient RNAs removal needs to be taken into consideration, especially for those diseased tissues. To our best knowledge, our study also offers the first cortex-specific snRNA-seq data of deeper cerebral hypoperfusion, which provides with novel therapeutic targets.

Mapping visual functions onto molecular cell types in the mouse superior colliculus

The superficial superior colliculus (sSC) carries out diverse roles in visual processing and behaviors, but how these functions are delegated among collicular neurons remains unclear. Here, using single-cell transcriptomics, we identified 28 neuron subtypes and subtype-enriched marker genes from tens of thousands of adult mouse sSC neurons. We then asked whether the sSC's molecular subtypes are tuned to different visual stimuli. Specifically, we imaged calcium dynamics in single sSC neurons in vivo during visual stimulation and then mapped marker gene transcripts onto the same neurons ex vivo. Our results identify a molecular subtype of inhibitory neuron accounting for ∼50% of the sSC's direction-selective cells, suggesting a genetic logic for the functional organization of the sSC. In addition, our studies provide a comprehensive molecular atlas of sSC neuron subtypes and a multimodal mapping method that will facilitate investigation of their respective functions, connectivity, and development.

Cochlear transcriptome analysis of an outbred mouse population (CFW)

Age-related hearing loss (ARHL) is the most common cause of hearing loss and one of the most prevalent conditions affecting the elderly worldwide. Despite evidence from our lab and others about its polygenic nature, little is known about the specific genes, cell types and pathways involved in ARHL, impeding the development of therapeutic interventions. In this manuscript, we describe, for the first time, the complete cell-type specific transcriptome of the aging mouse cochlea using snRNA-seq in an outbred mouse model in relation to auditory threshold variation. Cochlear cell types were identified using unsupervised clustering and annotated via a three-tiered approach - first by linking to expression of known marker genes, then using the NS-Forest algorithm to select minimum cluster-specific marker genes and reduce dimensional feature space for statistical comparison of our clusters with existing publicly-available data sets on the gEAR website (https://umgear.org/), and finally, by validating and refining the annotations using Multiplexed Error Robust Fluorescence In Situ Hybridization (MERFISH) and the cluster-specific marker genes as probes. We report on 60 unique cell-types expanding the number of defined cochlear cell types by more than two times. Importantly, we show significant specific cell type increases and decreases associated with loss of hearing acuity implicating specific subsets of hair cell subtypes, ganglion cell subtypes, and cell subtypes withing the stria vascularis in this model of ARHL. These results provide a view into the cellular and molecular mechanisms responsible for age-related hearing loss and pathways for therapeutic targeting.

Mapping human adult hippocampal neurogenesis with single-cell transcriptomics: Reconciling controversy or fueling the debate?

The notion of exploiting the regenerative potential of the human brain in physiological aging or neurological diseases represents a particularly attractive alternative to conventional strategies for enhancing or restoring brain function. However, a major first question to address is whether the human brain does possess the ability to regenerate. The existence of human adult hippocampal neurogenesis (AHN) has been at the center of a fierce scientific debate for many years. The advent of single-cell transcriptomic technologies was initially viewed as a panacea to resolving this controversy. However, recent single-cell RNA sequencing studies in the human hippocampus yielded conflicting results. Here, we critically discuss and re-analyze previously published AHN-related single-cell transcriptomic datasets. We argue that, although promising, the single-cell transcriptomic profiling of AHN in the human brain can be confounded by methodological, conceptual, and biological factors that need to be consistently addressed across studies and openly discussed within the scientific community.

Single-Nucleus RNA Sequencing of Developing and Mature Superior Colliculus Identifies Neuronal Diversity and Candidate Mediators of Circuit Assembly

The superior colliculus (SC) is a sensorimotor structure in the midbrain that integrates input from multiple sensory modalities to initiate motor commands. It undergoes well-characterized steps of circuit assembly during development, rendering the mouse SC a popular model to study establishment and refinement of neural connectivity. Here we performed single nucleus RNA-sequencing analysis of the mouse SC isolated at various developmental time points. Our study provides a transcriptomic landscape of the cell types that comprise the SC across murine development with particular emphasis on neuronal heterogeneity. We used these data to identify Pax7 as a marker for an anatomically homogeneous population of GABAergic neurons. Lastly, we report a repertoire of genes differentially expressed across the different postnatal ages, many of which are known to regulate axon guidance and synapse formation. Our data provide a valuable resource for interrogating the mechanisms of circuit development, and identifying markers for manipulating specific SC neuronal populations and circuits.

Single-cell and spatial transcriptomics: deciphering brain complexity in health and disease

In the past decade, single-cell technologies have proliferated and improved from their technically challenging beginnings to become common laboratory methods capable of determining the expression of thousands of genes in thousands of cells simultaneously. The field has progressed by taking the CNS as a primary research subject - the cellular complexity and multiplicity of neuronal cell types provide fertile ground for the increasing power of single-cell methods. Current single-cell RNA sequencing methods can quantify gene expression with sufficient accuracy to finely resolve even subtle differences between cell types and states, thus providing a great tool for studying the molecular and cellular repertoire of the CNS and its disorders. However, single-cell RNA sequencing requires the dissociation of tissue samples, which means that the interrelationships between cells are lost. Spatial transcriptomic methods bypass tissue dissociation and retain this spatial information, thereby allowing gene expression to be assessed across thousands of cells within the context of tissue structural organization. Here, we discuss how single-cell and spatially resolved transcriptomics have been contributing to unravelling the pathomechanisms underlying brain disorders. We focus on three areas where we feel these new technologies have provided particularly useful insights: selective neuronal vulnerability, neuroimmune dysfunction and cell-type-specific treatment response. We also discuss the limitations and future directions of single-cell and spatial RNA sequencing technologies.

Neural cell diversity in the light of single-cell transcriptomics

The development of highly parallel and affordable high-throughput single-cell transcriptomics technologies has revolutionized our understanding of brain complexity. These methods have been used to build cellular maps of the brain, its different regions, and catalog the diversity of cells in each of them during development, aging and even in disease. Now we know that cellular diversity is way beyond what was previously thought. Single-cell transcriptomics analyses have revealed that cell types previously considered homogeneous based on imaging techniques differ depending on several factors including sex, age and location within the brain. The expression profiles of these cells have also been exploited to understand which are the regulatory programs behind cellular diversity and decipher the transcriptional pathways driving them. In this review, we summarize how single-cell transcriptomics have changed our view on the cellular diversity in the human brain, and how it could impact the way we study neurodegenerative diseases. Moreover, we describe the new computational approaches that can be used to study cellular differentiation and gain insight into the functions of individual cell populations under different conditions and their alterations in disease.

Approaches for single-cell RNA sequencing across tissues and cell types

Single-cell sequencing of RNA (scRNA-seq) has advanced our understanding of cellular heterogeneity and signaling in developmental biology and disease. A large number of complementary assays have been developed to profile transcriptomes of individual cells, also in combination with other readouts, such as chromatin accessibility or antibody-based analysis of protein surface markers. As scRNA-seq technologies are advancing fast, it is challenging to establish robust workflows and up-to-date protocols that are best suited to address the large range of research questions. Here, we review scRNA-seq techniques from mRNA end-counting to total RNA in relation to their specific features and outline the necessary sample preparation steps and quality control measures. Based on our experience in dealing with the continuously growing portfolio from the perspective of a central single-cell facility, we aim to provide guidance on how workflows can be best automatized and share our experience in coping with the continuous expansion of scRNA-seq techniques.

Dissociation protocols used for sarcoma tissues bias the transcriptome observed in single-cell and single-nucleus RNA sequencing

BACKGROUND

Single-cell RNA-seq has emerged as an innovative technology used to study complex tissues and characterize cell types, states, and lineages at a single-cell level. Classification of bulk tumors by their individual cellular constituents has also created new opportunities to generate single-cell atlases for many organs, cancers, and developmental models. Despite the tremendous promise of this technology, recent evidence studying epithelial tissues and diverse carcinomas suggests the methods used for tissue processing, cell disaggregation, and preservation can significantly bias gene expression and alter the observed cell types. To determine whether sarcomas - tumors of mesenchymal origin - are subject to the same technical artifacts, we profiled patient-derived tumor explants (PDXs) propagated from three aggressive subtypes: osteosarcoma (OS), Ewing sarcoma (ES), desmoplastic small round cell tumor (DSRCT). Given the rarity of these sarcoma subtypes, we explored whether single-nuclei RNA-seq from more widely available archival frozen specimens could accurately be identified by gene expression signatures linked to tissue phenotype or pathognomonic fusion proteins.

RESULTS

We systematically assessed dissociation methods across different sarcoma subtypes. We compared gene expression from single-cell and single-nucleus RNA-sequencing of 125,831 whole-cells and nuclei from ES, DSRCT, and OS PDXs. We detected warm dissociation artifacts in single-cell samples and gene length bias in single-nucleus samples. Classic sarcoma gene signatures were observed regardless of the dissociation method. In addition, we showed that dissociation method biases could be computationally corrected.

CONCLUSIONS

We highlighted transcriptional biases, including warm dissociation and gene-length biases, introduced by the dissociation method for various sarcoma subtypes. This work is the first to characterize how the dissociation methods used for sc/snRNA-seq may affect the interpretation of the molecular features in sarcoma PDXs.

Integrated multimodal cell atlas of Alzheimer's disease

Alzheimer's disease (AD) is the most common cause of dementia in older adults. Neuropathological and imaging studies have demonstrated a progressive and stereotyped accumulation of protein aggregates, but the underlying molecular and cellular mechanisms driving AD progression and vulnerable cell populations affected by disease remain coarsely understood. The current study harnesses single cell and spatial genomics tools and knowledge from the BRAIN Initiative Cell Census Network to understand the impact of disease progression on middle temporal gyrus cell types. We used image-based quantitative neuropathology to place 84 donors spanning the spectrum of AD pathology along a continuous disease pseudoprogression score and multiomic technologies to profile single nuclei from each donor, mapping their transcriptomes, epigenomes, and spatial coordinates to a common cell type reference with unprecedented resolution. Temporal analysis of cell-type proportions indicated an early reduction of Somatostatin-expressing neuronal subtypes and a late decrease of supragranular intratelencephalic-projecting excitatory and Parvalbumin-expressing neurons, with increases in disease-associated microglial and astrocytic states. We found complex gene expression differences, ranging from global to cell type-specific effects. These effects showed different temporal patterns indicating diverse cellular perturbations as a function of disease progression. A subset of donors showed a particularly severe cellular and molecular phenotype, which correlated with steeper cognitive decline. We have created a freely available public resource to explore these data and to accelerate progress in AD research at SEA-AD.org.