Single-Nucleus RNA-Seq Is Not Suitable for Detection of Microglial Activation Genes in Humans

Gene expression and chromatin conformation of microglia in virally suppressed people with HIV

The presence of HIV in sequestered reservoirs is a central impediment to a functional cure, allowing HIV to persist despite life-long antiretroviral therapy (ART), and driving a variety of comorbid conditions. Our understanding of the latent HIV reservoir in the central nervous system is incomplete, because of difficulties in accessing human central nervous system tissues. Microglia contribute to HIV reservoirs, but the molecular phenotype of HIV-infected microglia is poorly understood. We leveraged the unique "Last Gift" rapid autopsy program, in which people with HIV are closely followed until days or even hours before death. Microglial populations were heterogeneous regarding their gene expression profiles but showed similar chromatin accessibility landscapes. Despite ART, we detected occasional microglia containing cell-associated HIV RNA and HIV DNA integrated into open regions of the host's genome (∼0.005%). Microglia with detectable HIV RNA showed an inflammatory phenotype. These results demonstrate a distinct myeloid cell reservoir in the brains of people with HIV despite suppressive ART. Strategies for curing HIV and neurocognitive impairment will need to consider the myeloid compartment to be successful.

Profiling cell identity and tissue architecture with single-cell and spatial transcriptomics

Single-cell transcriptomics has broadened our understanding of cellular diversity and gene expression dynamics in healthy and diseased tissues. Recently, spatial transcriptomics has emerged as a tool to contextualize single cells in multicellular neighbourhoods and to identify spatially recurrent phenotypes, or ecotypes. These technologies have generated vast datasets with targeted-transcriptome and whole-transcriptome profiles of hundreds to millions of cells. Such data have provided new insights into developmental hierarchies, cellular plasticity and diverse tissue microenvironments, and spurred a burst of innovation in computational methods for single-cell analysis. In this Review, we discuss recent advancements, ongoing challenges and prospects in identifying and characterizing cell states and multicellular neighbourhoods. We discuss recent progress in sample processing, data integration, identification of subtle cell states, trajectory modelling, deconvolution and spatial analysis. Furthermore, we discuss the increasing application of deep learning, including foundation models, in analysing single-cell and spatial transcriptomics data. Finally, we discuss recent applications of these tools in the fields of stem cell biology, immunology, and tumour biology, and the future of single-cell and spatial transcriptomics in biological research and its translation to the clinic.

Purinergic exposure induces epigenomic and transcriptomic-mediated preconditioning resembling epilepsy-associated microglial states

Microglia play a crucial role in a range of neuropathologies through exacerbated activation. Microglial inflammatory responses can be influenced by prior exposures to noxious stimuli, like increased levels of extracellular adenosine and ATP. These are characteristic of brain insults like epileptic seizures and could potentially shape subsequent responses through epigenetic regulation. We investigated DNA methylation and expression changes in human microglia-like cells differentiated from monocytes following ATP-mediated preconditioning. We demonstrate that microglia-like cells display homeostatic microglial features, shown by surface markers, transcriptome, and DNA methylome. After exposure to ATP, TLR-mediated activation leads to an exacerbated pro-inflammatory response. These changes are accompanied by methylation and transcriptional reprogramming associated with enhanced immune-related functions. The reprogramming associated with ATP-mediated preconditioning leads to profiles found in microglial subsets linked to epilepsy. Purine-driven microglia immune preconditioning drives epigenetic and transcriptional changes that could contribute to altered functions of microglia during seizure development and progression.

Integrative multi-omic analysis reveals conserved cell-projection deficits in human Down syndrome brains

Down syndrome (DS) is the most common genetic cause of cognitive disability. However, it is largely unclear how triplication of a small gene subset may impinge on diverse aspects of DS brain physiopathology. Here, we took a multi-omic approach and simultaneously analyzed by RNA-seq and proteomics the expression signatures of two diverse regions of human postmortem DS brains. We found that the overexpression of triplicated genes triggered global expression dysregulation, differentially affecting transcripts, miRNAs, and proteins involved in both known and novel biological candidate pathways. Among the latter, we observed an alteration in RNA splicing, specifically modulating the expression of genes involved in cytoskeleton and axonal dynamics in DS brains. Accordingly, we found an alteration in axonal polarization in neurons from DS human iPSCs and mice. Thus, our study provides an integrated multilayer expression database capable of identifying new potential targets to aid in designing future clinical interventions for DS.

Single-cell multiregion dissection of Alzheimer's disease

Alzheimer's disease is the leading cause of dementia worldwide, but the cellular pathways that underlie its pathological progression across brain regions remain poorly understood1-3. Here we report a single-cell transcriptomic atlas of six different brain regions in the aged human brain, covering 1.3 million cells from 283 post-mortem human brain samples across 48 individuals with and without Alzheimer's disease. We identify 76 cell types, including region-specific subtypes of astrocytes and excitatory neurons and an inhibitory interneuron population unique to the thalamus and distinct from canonical inhibitory subclasses. We identify vulnerable populations of excitatory and inhibitory neurons that are depleted in specific brain regions in Alzheimer's disease, and provide evidence that the Reelin signalling pathway is involved in modulating the vulnerability of these neurons. We develop a scalable method for discovering gene modules, which we use to identify cell-type-specific and region-specific modules that are altered in Alzheimer's disease and to annotate transcriptomic differences associated with diverse pathological variables. We identify an astrocyte program that is associated with cognitive resilience to Alzheimer's disease pathology, tying choline metabolism and polyamine biosynthesis in astrocytes to preserved cognitive function late in life. Together, our study develops a regional atlas of the ageing human brain and provides insights into cellular vulnerability, response and resilience to Alzheimer's disease pathology.

Transcriptomic Approaches to Cardiomyocyte-Biomaterial Interactions: A Review

Biomaterials, essential for supporting, enhancing, and repairing damaged tissues, play a critical role in various medical applications. This Review focuses on the interaction of biomaterials and cardiomyocytes, emphasizing the unique significance of transcriptomic approaches in understanding their interactions, which are pivotal in cardiac bioengineering and regenerative medicine. Transcriptomic approaches serve as powerful tools to investigate how cardiomyocytes respond to biomaterials, shedding light on the gene expression patterns, regulatory pathways, and cellular processes involved in these interactions. Emerging technologies such as bulk RNA-seq, single-cell RNA-seq, single-nucleus RNA-seq, and spatial transcriptomics offer promising avenues for more precise and in-depth investigations. Longitudinal studies, pathway analyses, and machine learning techniques further improve the ability to explore the complex regulatory mechanisms involved. This review also discusses the challenges and opportunities of utilizing transcriptomic techniques in cardiomyocyte-biomaterial research. Although there are ongoing challenges such as costs, cell size limitation, sample differences, and complex analytical process, there exist exciting prospects in comprehensive gene expression analyses, biomaterial design, cardiac disease treatment, and drug testing. These multimodal methodologies have the capacity to deepen our understanding of the intricate interaction network between cardiomyocytes and biomaterials, potentially revolutionizing cardiac research with the aim of promoting heart health, and they are also promising for studying interactions between biomaterials and other cell types.

Proteomic features of gray matter layers and superficial white matter of the rhesus monkey neocortex: comparison of prefrontal area 46 and occipital area 17

In this novel large-scale multiplexed immunofluorescence study we comprehensively characterized and compared layer-specific proteomic features within regions of interest of the widely divergent dorsolateral prefrontal cortex (A46) and primary visual cortex (A17) of adult rhesus monkeys. Twenty-eight markers were imaged in rounds of sequential staining, and their spatial distribution precisely quantified within gray matter layers and superficial white matter. Cells were classified as neurons, astrocytes, oligodendrocytes, microglia, or endothelial cells. The distribution of fibers and blood vessels were assessed by quantification of staining intensity across regions of interest. This method revealed multivariate similarities and differences between layers and areas. Protein expression in neurons was the strongest determinant of both laminar and regional differences, whereas protein expression in glia was more important for intra-areal laminar distinctions. Among specific results, we observed a lower glia-to-neuron ratio in A17 than in A46 and the pan-neuronal markers HuD and NeuN were differentially distributed in both brain areas with a lower intensity of NeuN in layers 4 and 5 of A17 compared to A46 and other A17 layers. Astrocytes and oligodendrocytes exhibited distinct marker-specific laminar distributions that differed between regions; notably, there was a high proportion of ALDH1L1-expressing astrocytes and of oligodendrocyte markers in layer 4 of A17. The many nuanced differences in protein expression between layers and regions observed here highlight the need for direct assessment of proteins, in addition to RNA expression, and set the stage for future protein-focused studies of these and other brain regions in normal and pathological conditions.

SIMS: A deep-learning label transfer tool for single-cell RNA sequencing analysis

Cell atlases serve as vital references for automating cell labeling in new samples, yet existing classification algorithms struggle with accuracy. Here we introduce SIMS (scalable, interpretable machine learning for single cell), a low-code data-efficient pipeline for single-cell RNA classification. We benchmark SIMS against datasets from different tissues and species. We demonstrate SIMS's efficacy in classifying cells in the brain, achieving high accuracy even with small training sets (<3,500 cells) and across different samples. SIMS accurately predicts neuronal subtypes in the developing brain, shedding light on genetic changes during neuronal differentiation and postmitotic fate refinement. Finally, we apply SIMS to single-cell RNA datasets of cortical organoids to predict cell identities and uncover genetic variations between cell lines. SIMS identifies cell-line differences and misannotated cell lineages in human cortical organoids derived from different pluripotent stem cell lines. Altogether, we show that SIMS is a versatile and robust tool for cell-type classification from single-cell datasets.

A cell atlas of the adult female Aedes aegypti midgut revealed by single-cell RNA sequencing

Aedes aegypti is a primary vector for transmitting various arboviruses, including Yellow fever, dengue and Zika virus. The mosquito midgut is the principal organ for blood meal digestion, nutrient absorption and the initial site of arbovirus infection. Although a previous study delineated midgut's transcriptome of Ae. aegypti at the single-nucleus resolution, there still lacks an established protocol for isolating and RNA sequencing of single cells of Ae. aegypti midgut, which is required for investigating arbovirus-midgut interaction at the single-cell level. Here, we established an atlas of the midgut cells for Ae. aegypti by single-cell RNA sequencing. We annotated the cell clusters including intestinal stem cells/enteroblasts (ISC/EB), cardia cells (Cardia), enterocytes (EC, EC-like), enteroendocrine cells (EE), visceral muscle (VM), fat body cells (FBC) and hemocyte cells (HC). This study will provide a foundation for further studies of arbovirus infection in mosquito midgut at the single-cell level.

Inter- and trans-generational impacts of environmental exposures on the germline resolved at the single-cell level

Reproduction is a remarkably intricate process involving the interaction of multiple cell types and organ systems unfolding over long periods of time and that culminates with the production of gametes. The initiation of germ cell development takes place during embryogenesis but only completes decades later in humans. The complexity inherent to reproduction and its study has long hampered our ability to decipher how environmental agents disrupt this process. Single-cell approaches provide an opportunity for a deeper understanding of the action of toxicants on germline function and analyze how the response to their exposure is differentially distributed across tissues and cell types. In addition to single-cell RNA sequencing, other single-cell or nucleus level approaches such as ATAC-sequencing and multi-omics have expanded the strategies that can be implemented in reproductive toxicological studies to include epigenomic and the nuclear transcriptomic data. Here we will discuss the current state of single-cell technologies and how they can best be utilized to advance reproductive toxicological studies. We will then discuss case studies in two model organisms (Caenorhabditis elegans and mouse) studying different environmental exposures (alcohol and e-cigarettes respectively) to highlight the value of single-cell and single-nucleus approaches for reproductive biology and reproductive toxicology.

The impact of chronic pain on brain gene expression

Background

Chronic pain affects one fifth of American adults, contributing significant public health burden. Chronic pain mechanisms can be further understood through investigating brain gene expression.

Methods

We tested differentially expressed genes (DEGs) in chronic pain, migraine, lifetime fentanyl and oxymorphone use, and with chronic pain genetic risk in four brain regions (dACC, DLPFC, MeA, BLA) and imputed cell type expression data from 304 postmortem donors. We compared findings across traits and with independent transcriptomics resources, and performed gene-set enrichment.

Results

We identified two chronic pain DEGs: B4GALT and VEGFB in bulk dACC. We found over 2000 (primarily BLA microglia) chronic pain cell type DEGs. Findings were enriched for mouse microglia pain genes, and for hypoxia and immune response. Cross-trait DEG overlap was minimal.

Conclusions

Chronic pain-associated gene expression is heterogeneous across cell type, largely distinct from that in pain-related traits, and shows BLA microglia are a key cell type.

Cell type-specific gene expression dynamics during human brain maturation

The human brain undergoes protracted post-natal maturation, guided by dynamic changes in gene expression. Most studies exploring these processes have used bulk tissue analyses, which mask cell type-specific gene expression dynamics. Here, using single nucleus (sn)RNA-seq on temporal lobe tissue, including samples of African ancestry, we build a joint paediatric and adult atlas of 75 cell subtypes, which we verify with spatial transcriptomics. We explore the differences between paediatric and adult cell types, revealing the genes and pathways that change during brain maturation. Our results highlight excitatory neuron subtypes, including the LTK and FREM subtypes, that show elevated expression of genes associated with cognition and synaptic plasticity in paediatric tissue. The new resources we present here improve our understanding of the brain during its development and contribute to global efforts to build an inclusive brain cell map.

Xenografted human microglia display diverse transcriptomic states in response to Alzheimer's disease-related amyloid-β pathology

Microglia are central players in Alzheimer's disease pathology but analyzing microglial states in human brain samples is challenging due to genetic diversity, postmortem delay and admixture of pathologies. To circumvent these issues, here we generated 138,577 single-cell expression profiles of human stem cell-derived microglia xenotransplanted in the brain of the AppNL-G-F model of amyloid pathology and wild-type controls. Xenografted human microglia adopt a disease-associated profile similar to that seen in mouse microglia, but display a more pronounced human leukocyte antigen or HLA state, likely related to antigen presentation in response to amyloid plaques. The human microglial response also involves a pro-inflammatory cytokine/chemokine cytokine response microglia or CRM response to oligomeric Aβ oligomers. Genetic deletion of TREM2 or APOE as well as APOE polymorphisms and TREM2R47H expression in the transplanted microglia modulate these responses differentially. The expression of other Alzheimer's disease risk genes is differentially regulated across the distinct cell states elicited in response to amyloid pathology. Thus, we have identified multiple transcriptomic cell states adopted by human microglia in a multipronged response to Alzheimer's disease-related pathology, which should be taken into account in translational studies.

Epigenetic control of microglial immune responses

Microglia, the major population of brain-resident macrophages, are now recognized as a heterogeneous population comprising several cell subtypes with different (so far mostly supposed) functions in health and disease. A number of studies have performed molecular characterization of these different microglial activation states over the last years making use of "omics" technologies, that is transcriptomics, proteomics and, less frequently, epigenomics profiling. These approaches offer the possibility to identify disease mechanisms, discover novel diagnostic biomarkers, and develop new therapeutic strategies. Here, we focus on epigenetic profiling as a means to understand microglial immune responses beyond what other omics methods can offer, that is, revealing past and present molecular responses, gene regulatory networks and potential future response trajectories, and defining cell subtype-specific disease relevance through mapping non-coding genetic variants. We review the current knowledge in the field regarding epigenetic regulation of microglial identity and function, provide an exemplary analysis that demonstrates the advantages of performing joint transcriptomic and epigenomic profiling of single microglial cells and discuss how comprehensive epigenetic analyses may enhance our understanding of microglial pathophysiology.

Transcript Profiles of Microglia/Macrophage Cells at the Borders of Chronic Active and Subpial Gray Matter Lesions in Multiple Sclerosis

OBJECTIVE

Microglia/macrophages line the border of demyelinated lesions in both cerebral white matter and the cortex in the brains of multiple sclerosis patients. Microglia/macrophages associated with chronic white matter lesions are thought to be responsible for slow lesion expansion and disability progression in progressive multiple sclerosis, whereas those lining gray matter lesions are less studied. Profiling these microglia/macrophages could help to focus therapies on genes or pathways specific to lesion expansion and disease progression.

METHODS

We compared the morphology and transcript profiles of microglia/macrophages associated with borders of white matter (WM line) and subpial gray matter lesions (GM line) using laser capture microscopy. We performed RNA sequencing on isolated cells followed by immunocytochemistry to determine the distribution of translational products of transcripts increased in WM line microglia.

RESULTS

Cells in the WM line appear activated, with shorter processes and larger cell bodies, whereas those in the GM line appear more homeostatic, with smaller cell bodies and multiple thin processes. Transcript profiling revealed 176 genes in WM lines and 111 genes in GM lines as differentially expressed. Transcripts associated with immune activation and iron homeostasis were increased in WM line microglia, whereas genes belonging to the canonical Wnt signaling pathway were increased in GM line microglia.

INTERPRETATION

We propose that the mechanisms of demyelination and dynamics of lesion expansion are responsible for differential transcript expression in WM lines and GM lines, and posit that increased expression of the Fc epsilon receptor, spleen tyrosine kinase, and Bruton's tyrosine kinase, play a key role in regulating microglia/macrophage function at the border of chronic active white matter lesions. ANN NEUROL 2024;95:907-916.

Unraveling the enigma: housekeeping gene Ugt1a7c as a universal biomarker for microglia

Background

Microglia, brain resident macrophages, play multiple roles in maintaining homeostasis, including immunity, surveillance, and protecting the central nervous system through their distinct activation processes. Identifying all types of microglia-driven populations is crucial due to the presence of various phenotypes that differ based on developmental stages or activation states. During embryonic development, the E8.5 yolk sac contains erythromyeloid progenitors that go through different growth phases, eventually resulting in the formation of microglia. In addition, microglia are present in neurological diseases as a diverse population. So far, no individual biomarker for microglia has been discovered that can accurately identify and monitor their development and attributes.

Summary

Here, we highlight the newly defined biomarker of mouse microglia, UGT1A7C, which exhibits superior stability in expression during microglia development and activation compared to other known microglia biomarkers. The UGT1A7C sensing chemical probe labels all microglia in the 3xTG AD mouse model. The expression of Ugt1a7c is stable during development, with only a 4-fold variation, while other microglia biomarkers, such as Csf1r and Cx3cr1, exhibit at least a 10-fold difference. The UGT1A7C expression remains constant throughout its lifespan. In addition, the expression and activity of UGT1A7C are the same in response to different types of inflammatory activators' treatment in vitro.

Conclusion

We propose employing UGT1A7C as the representative biomarker for microglia, irrespective of their developmental state, age, or activation status. Using UGT1A7C can reduce the requirement for using multiple biomarkers, enhance the precision of microglia analysis, and even be utilized as a standard for gene/protein expression.

Single-cell dissection of the human motor and prefrontal cortices in ALS and FTLD

Amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) share many clinical, pathological, and genetic features, but a detailed understanding of their associated transcriptional alterations across vulnerable cortical cell types is lacking. Here, we report a high-resolution, comparative single-cell molecular atlas of the human primary motor and dorsolateral prefrontal cortices and their transcriptional alterations in sporadic and familial ALS and FTLD. By integrating transcriptional and genetic information, we identify known and previously unidentified vulnerable populations in cortical layer 5 and show that ALS- and FTLD-implicated motor and spindle neurons possess a virtually indistinguishable molecular identity. We implicate potential disease mechanisms affecting these cell types as well as non-neuronal drivers of pathogenesis. Finally, we show that neuron loss in cortical layer 5 tracks more closely with transcriptional identity rather than cellular morphology and extends beyond previously reported vulnerable cell types.

Omics Approaches Unveiling the Biology of Human Atherosclerotic Plaques

Atherosclerosis is a chronic inflammatory disease of the arterial wall, characterized by the buildup of plaques with the accumulation and transformation of lipids, immune cells, vascular smooth muscle cells, and necrotic cell debris. Plaques with collagen-poor thin fibrous caps infiltrated by macrophages and lymphocytes are considered unstable because they are at the greatest risk of rupture and clinical events. However, the current histologic definition of plaque types may not fully capture the complex molecular nature of atherosclerotic plaque biology and the underlying mechanisms contributing to plaque progression, rupture, and erosion. The advances in omics technologies have changed the understanding of atherosclerosis plaque biology, offering new possibilities to improve risk prediction and discover novel therapeutic targets. Genomic studies have shed light on the genetic predisposition to atherosclerosis, and integrative genomic analyses expedite the translation of genomic discoveries. Transcriptomic, proteomic, metabolomic, and lipidomic studies have refined the understanding of the molecular signature of atherosclerotic plaques, aiding in data-driven hypothesis generation for mechanistic studies and offering new prospects for biomarker discovery. Furthermore, advancements in single-cell technologies and emerging spatial analysis techniques have unveiled the heterogeneity and plasticity of plaque cells. This review discusses key omics-based discoveries that have advanced the understanding of human atherosclerotic plaque biology, focusing on insights derived from omics profiling of human atherosclerotic vascular specimens.

Differences in molecular sampling and data processing explain variation among single-cell and single-nucleus RNA-seq experiments

A mechanistic understanding of the biological and technical factors that impact transcript measurements is essential to designing and analyzing single-cell and single-nucleus RNA sequencing experiments. Nuclei contain the same pre-mRNA population as cells, but they contain a small subset of the mRNAs. Nonetheless, early studies argued that single-nucleus analysis yielded results comparable to cellular samples if pre-mRNA measurements were included. However, typical workflows do not distinguish between pre-mRNA and mRNA when estimating gene expression, and variation in their relative abundances across cell types has received limited attention. These gaps are especially important given that incorporating pre-mRNA has become commonplace for both assays, despite known gene length bias in pre-mRNA capture. Here, we reanalyze public data sets from mouse and human to describe the mechanisms and contrasting effects of mRNA and pre-mRNA sampling on gene expression and marker gene selection in single-cell and single-nucleus RNA-seq. We show that pre-mRNA levels vary considerably among cell types, which mediates the degree of gene length bias and limits the generalizability of a recently published normalization method intended to correct for this bias. As an alternative, we repurpose an existing post hoc gene length-based correction method from conventional RNA-seq gene set enrichment analysis. Finally, we show that inclusion of pre-mRNA in bioinformatic processing can impart a larger effect than assay choice itself, which is pivotal to the effective reuse of existing data. These analyses advance our understanding of the sources of variation in single-cell and single-nucleus RNA-seq experiments and provide useful guidance for future studies.

A single nuclear transcriptomic characterisation of mechanisms responsible for impaired angiogenesis and blood-brain barrier function in Alzheimer's disease

Brain perfusion and blood-brain barrier (BBB) integrity are reduced early in Alzheimer's disease (AD). We performed single nucleus RNA sequencing of vascular cells isolated from AD and non-diseased control brains to characterise pathological transcriptional signatures responsible for this. We show that endothelial cells (EC) are enriched for expression of genes associated with susceptibility to AD. Increased β-amyloid is associated with BBB impairment and a dysfunctional angiogenic response related to a failure of increased pro-angiogenic HIF1A to increased VEGFA signalling to EC. This is associated with vascular inflammatory activation, EC senescence and apoptosis. Our genomic dissection of vascular cell risk gene enrichment provides evidence for a role of EC pathology in AD and suggests that reducing vascular inflammatory activation and restoring effective angiogenesis could reduce vascular dysfunction contributing to the genesis or progression of early AD.

BHLHE40/41 regulate microglia and peripheral macrophage responses associated with Alzheimer's disease and other disorders of lipid-rich tissues

Genetic and experimental evidence suggests that Alzheimer's disease (AD) risk alleles and genes may influence disease susceptibility by altering the transcriptional and cellular responses of macrophages, including microglia, to damage of lipid-rich tissues like the brain. Recently, sc/nRNA sequencing studies identified similar transcriptional activation states in subpopulations of macrophages in aging and degenerating brains and in other diseased lipid-rich tissues. We collectively refer to these subpopulations of microglia and peripheral macrophages as DLAMs. Using macrophage sc/nRNA-seq data from healthy and diseased human and mouse lipid-rich tissues, we reconstructed gene regulatory networks and identified 11 strong candidate transcriptional regulators of the DLAM response across species. Loss or reduction of two of these transcription factors, BHLHE40/41, in iPSC-derived microglia and human THP-1 macrophages as well as loss of Bhlhe40/41 in mouse microglia, resulted in increased expression of DLAM genes involved in cholesterol clearance and lysosomal processing, increased cholesterol efflux and storage, and increased lysosomal mass and degradative capacity. These findings provide targets for therapeutic modulation of macrophage/microglial function in AD and other disorders affecting lipid-rich tissues.

Mapping the cellular landscape of Atlantic salmon head kidney by single cell and single nucleus transcriptomics

Single-cell transcriptomics is the current gold standard for global gene expression profiling, not only in mammals and model species, but also in non-model fish species. This is a rapidly expanding field, creating a deeper understanding of tissue heterogeneity and the distinct functions of individual cells, making it possible to explore the complexities of immunology and gene expression on a highly resolved level. In this study, we compared two single cell transcriptomic approaches to investigate cellular heterogeneity within the head kidney of healthy farmed Atlantic salmon (Salmo salar). We compared 14,149 cell transcriptomes assayed by single cell RNA-seq (scRNA-seq) with 18,067 nuclei transcriptomes captured by single nucleus RNA-Seq (snRNA-seq). Both approaches detected eight major cell populations in common: granulocytes, heamatopoietic stem cells, erythrocytes, mononuclear phagocytes, thrombocytes, B cells, NK-like cells, and T cells. Four additional cell types, endothelial, epithelial, interrenal, and mesenchymal cells, were detected in the snRNA-seq dataset, but appeared to be lost during preparation of the single cell suspension submitted for scRNA-seq library generation. We identified additional heterogeneity and subpopulations within the B cells, T cells, and endothelial cells, and revealed developmental trajectories of heamatopoietic stem cells into differentiated granulocyte and mononuclear phagocyte populations. Gene expression profiles of B cell subtypes revealed distinct IgM and IgT-skewed resting B cell lineages and provided insights into the regulation of B cell lymphopoiesis. The analysis revealed eleven T cell sub-populations, displaying a level of T cell heterogeneity in salmon head kidney comparable to that observed in mammals, including distinct subsets of cd4/cd8-negative T cells, such as tcrγ positive, progenitor-like, and cytotoxic cells. Although snRNA-seq and scRNA-seq were both useful to resolve cell type-specific expression in the Atlantic salmon head kidney, the snRNA-seq pipeline was overall more robust in identifying several cell types and subpopulations. While scRNA-seq displayed higher levels of ribosomal and mitochondrial genes, snRNA-seq captured more transcription factor genes. However, only scRNA-seq-generated data was useful for cell trajectory inference within the myeloid lineage. In conclusion, this study systematically outlines the relative merits of scRNA-seq and snRNA-seq in Atlantic salmon, enhances understanding of teleost immune cell lineages, and provides a comprehensive list of markers for identifying major cell populations in the head kidney with significant immune relevance.

High-resolution single-cell transcriptomic survey of cardiomyocytes from patients with hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is a common inherited cardiovascular disease, which can cause heart failure and lead to death. In this study, we performed high-resolution single-cell RNA-sequencing of 2115 individual cardiomyocytes obtained from HCM patients and normal controls. Signature up- and down-regulated genes in HCM were identified by integrative analysis across 37 patients and 41 controls from our data and published human single-cell and single-nucleus RNA-seq datasets, which were further classified into gene modules by single-cell co-expression analysis. Using our high-resolution dataset, we also investigated the heterogeneity among individual cardiomyocytes and revealed five distinct clusters within HCM cardiomyocytes. Interestingly, we showed that some extracellular matrix (ECM) genes were up-regulated in the HCM cardiomyocytes, suggesting that they play a role in cardiac remodelling. Taken together, our study comprehensively profiled the transcriptomic programs of HCM cardiomyocytes and provided insights into molecular mechanisms underlying the pathogenesis of HCM.

Reduced progranulin increases tau and α-synuclein inclusions and alters mouse tauopathy phenotypes via glucocerebrosidase

Comorbid proteinopathies are observed in many neurodegenerative disorders including Alzheimer's disease (AD), increase with age, and influence clinical outcomes, yet the mechanisms remain ill-defined. Here, we show that reduction of progranulin (PGRN), a lysosomal protein associated with TDP-43 proteinopathy, also increases tau inclusions, causes concomitant accumulation of α-synuclein and worsens mortality and disinhibited behaviors in tauopathy mice. The increased inclusions paradoxically protect against spatial memory deficit and hippocampal neurodegeneration. PGRN reduction in male tauopathy attenuates activity of β-glucocerebrosidase (GCase), a protein previously associated with synucleinopathy, while increasing glucosylceramide (GlcCer)-positive tau inclusions. In neuronal culture, GCase inhibition enhances tau aggregation induced by AD-tau. Furthermore, purified GlcCer directly promotes tau aggregation in vitro. Neurofibrillary tangles in human tauopathies are also GlcCer-immunoreactive. Thus, in addition to TDP-43, PGRN regulates tau- and synucleinopathies via GCase and GlcCer. A lysosomal PGRN-GCase pathway may be a common therapeutic target for age-related comorbid proteinopathies.

Single-cell transcriptomics analysis of cellular heterogeneity and immune mechanisms in neurodegenerative diseases

Single-cell transcriptomics analysis is an advanced technology that can describe the intracellular transcriptome in complex tissues. It profiles and analyses datasets by single-cell RNA sequencing. Neurodegenerative diseases are identified by the abnormal apoptosis of neurons in the brain with few or no effective therapy strategies at present, which has been a growing healthcare concern and brought a great burden to society. The transcriptome of individual cells provides deep insights into previously unforeseen cellular heterogeneity and gene expression differences in neurodegenerative disorders. It detects multiple cell subsets and functional changes during pathological progression, which deepens the understanding of the molecular underpinnings and cellular basis of neurodegenerative diseases. Furthermore, the transcriptome analysis of immune cells shows the regulation of immune response. Different subtypes of immune cells and their interaction are found to contribute to disease progression. This finding enables the discovery of novel targets and biomarkers for early diagnosis. In this review, we emphasize the principles of the technology, and its recent progress in the study of cellular heterogeneity and immune mechanisms in neurodegenerative diseases. The application of single-cell transcriptomics analysis in neurodegenerative disorders would help explore the pathogenesis of these diseases and develop novel therapeutic methods.

Selective Vulnerability to Neurodegenerative Disease: Insights from Cell Type-Specific Translatome Studies

Neurodegenerative diseases (NDs) manifest a wide variety of clinical symptoms depending on the affected brain regions. Gaining insights into why certain regions are resistant while others are susceptible is vital for advancing therapeutic strategies. While gene expression changes offer clues about disease responses across brain regions, the mixture of cell types therein obscures experimental results. In recent years, methods that analyze the transcriptomes of individual cells (e.g., single-cell RNA sequencing or scRNAseq) have been widely used and have provided invaluable insights into specific cell types. Concurrently, transgene-based techniques that dissect cell type-specific translatomes (CSTs) in model systems, like RiboTag and bacTRAP, offer unique advantages but have received less attention. This review juxtaposes the merits and drawbacks of both methodologies, focusing on the use of CSTs in understanding conditions like amyotrophic lateral sclerosis (ALS), Huntington's disease (HD), Alzheimer's disease (AD), and specific prion diseases like fatal familial insomnia (FFI), genetic Creutzfeldt-Jakob disease (gCJD), and acquired prion disease. We conclude by discussing the emerging trends observed across multiple diseases and emerging methods.

Opioid-driven disruption of the septal complex reveals a role for neurotensin-expressing neurons in withdrawal

Because opioid withdrawal is an intensely aversive experience, persons with opioid use disorder (OUD) often relapse to avoid it. The lateral septum (LS) is a forebrain structure that is important in aversion processing, and previous studies have linked the lateral septum (LS) to substance use disorders. It is unclear, however, which precise LS cell types might contribute to the maladaptive state of withdrawal. To address this, we used single-nucleus RNA-sequencing to interrogate cell type specific gene expression changes induced by chronic morphine and withdrawal. We discovered that morphine globally disrupted the transcriptional profile of LS cell types, but Neurotensin-expressing neurons (Nts; LS-Nts neurons) were selectively activated by naloxone. Using two-photon calcium imaging and ex vivo electrophysiology, we next demonstrate that LS-Nts neurons receive enhanced glutamatergic drive in morphine-dependent mice and remain hyperactivated during opioid withdrawal. Finally, we showed that activating and silencing LS-Nts neurons during opioid withdrawal regulates pain coping behaviors and sociability. Together, these results suggest that LS-Nts neurons are a key neural substrate involved in opioid withdrawal and establish the LS as a crucial regulator of adaptive behaviors, specifically pertaining to OUD.

Transcriptomic signatures of individual cell types in cerebral cavernous malformation

Cerebral cavernous malformation (CCM) is a hemorrhagic neurovascular disease with no currently available therapeutics. Prior evidence suggests that different cell types may play a role in CCM pathogenesis. The contribution of each cell type to the dysfunctional cellular crosstalk remains unclear. Herein, RNA-seq was performed on fluorescence-activated cell sorted endothelial cells (ECs), pericytes, and neuroglia from CCM lesions and non-lesional brain tissue controls. Differentially Expressed Gene (DEG), pathway and Ligand-Receptor (LR) analyses were performed to characterize the dysfunctional genes of respective cell types within CCMs. Common DEGs among all three cell types were related to inflammation and endothelial-to-mesenchymal transition (EndMT). DEG and pathway analyses supported a role of lesional ECs in dysregulated angiogenesis and increased permeability. VEGFA was particularly upregulated in pericytes. Further pathway and LR analyses identified vascular endothelial growth factor A/ vascular endothelial growth factor receptor 2 signaling in lesional ECs and pericytes that would result in increased angiogenesis. Moreover, lesional pericytes and neuroglia predominantly showed DEGs and pathways mediating the immune response. Further analyses of cell specific gene alterations in CCM endorsed potential contribution to EndMT, coagulation, and a hypoxic microenvironment. Taken together, these findings motivate mechanistic hypotheses regarding non-endothelial contributions to lesion pathobiology and may lead to novel therapeutic targets. Video Abstract.

Single Cell Sequencing Technology and Its Application in Alzheimer's Disease

Alzheimer's disease (AD) involves degeneration of cells in the brain. Due to insidious onset and slow progression, AD is often not diagnosed until it gets progressed to a more severe stage. The diagnosis and treatment of AD has been a challenge. In recent years, high-throughput sequencing technologies have exhibited advantages in exploring the pathogenesis of diseases. However, the types of cells of the central nervous system are complex and traditional bulk sequencing cannot reflect their heterogeneity. Single-cell sequencing technology enables study at the individual cell level and has an irreplaceable advantage in the study of complex diseases. In recent years, this field has expanded rapidly and several types of single-cell sequencing technologies have emerged, including transcriptomics, epigenomics, genomics and proteomics. This review article provides an overview of these single-cell sequencing technologies and their application in AD.

Associations of cortical SPP1 and ITGAX with cognition and common neuropathologies in older adults

INTRODUCTION

The secreted phosphoprotein 1 (SPP1) gene expressed by CD11c+ cells is known to be associated with microglia activation and neuroinflammatory diseases. As most studies rely on mouse models, we investigated these genes and proteins in the cortical brain tissue of older adults and their role in Alzheimer's disease (AD) and related disorders.

METHODS

We leveraged protein measurements, single-nuclei, and RNASeq data from the Religious Orders Study and Rush Memory and Aging Project (ROSMAP) of over 1200 samples for association analysis.

RESULTS

Expression of SPP1 and its encoded protein osteopontin were associated with faster cognitive decline and greater odds of common neuropathologies. At single-cell resolution,  integrin subunit alpha X (ITGAX) was highly expressed in microglia, where specific subpopulations were associated with AD and cerebral amyloid angiopathy.

DISCUSSION

The study provides evidence of SPP1 and ITGAX association with cognitive decline and common neuropathologies identifying a microglial subset associated with disease.

Single-cell spatial multi-omics and deep learning dissect enhancer-driven gene regulatory networks in liver zonation

In the mammalian liver, hepatocytes exhibit diverse metabolic and functional profiles based on their location within the liver lobule. However, it is unclear whether this spatial variation, called zonation, is governed by a well-defined gene regulatory code. Here, using a combination of single-cell multiomics, spatial omics, massively parallel reporter assays and deep learning, we mapped enhancer-gene regulatory networks across mouse liver cell types. We found that zonation affects gene expression and chromatin accessibility in hepatocytes, among other cell types. These states are driven by the repressors TCF7L1 and TBX3, alongside other core hepatocyte transcription factors, such as HNF4A, CEBPA, FOXA1 and ONECUT1. To examine the architecture of the enhancers driving these cell states, we trained a hierarchical deep learning model called DeepLiver. Our study provides a multimodal understanding of the regulatory code underlying hepatocyte identity and their zonation state that can be used to engineer enhancers with specific activity levels and zonation patterns.

Investigative needle core biopsies for multi-omics in Glioblastoma

Glioblastoma (GBM) is a primary brain cancer with an abysmal prognosis and few effective therapies. The ability to investigate the tumor microenvironment before and during treatment would greatly enhance both understanding of disease response and progression, as well as the delivery and impact of therapeutics. Stereotactic biopsies are a routine surgical procedure performed primarily for diagnostic histopathologic purposes. The role of investigative biopsies - tissue sampling for the purpose of understanding tumor microenvironmental responses to treatment using integrated multi-modal molecular analyses ('Multi-omics") has yet to be defined. Secondly, it is unknown whether comparatively small tissue samples from brain biopsies can yield sufficient information with such methods. Here we adapt stereotactic needle core biopsy tissue in two separate patients. In the first patient with recurrent GBM we performed highly resolved multi-omics analysis methods including single cell RNA sequencing, spatial-transcriptomics, metabolomics, proteomics, phosphoproteomics, T-cell clonotype analysis, and MHC Class I immunopeptidomics from biopsy tissue that was obtained from a single procedure. In a second patient we analyzed multi-regional core biopsies to decipher spatial and genomic variance. We also investigated the utility of stereotactic biopsies as a method for generating patient derived xenograft models in a separate patient cohort. Dataset integration across modalities showed good correspondence between spatial modalities, highlighted immune cell associated metabolic pathways and revealed poor correlation between RNA expression and the tumor MHC Class I immunopeptidome. In conclusion, stereotactic needle biopsy cores are of sufficient quality to generate multi-omics data, provide data rich insight into a patient's disease process and tumor immune microenvironment and can be of value in evaluating treatment responses.

One sentence summary

Integrative multi-omics analysis of stereotactic needle core biopsies in glioblastoma.

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.

Spatiotemporal characterization of glial cell activation in an Alzheimer's disease model by spatially resolved transcriptomics

The molecular changes that occur with the progression of Alzheimer's disease (AD) are well known, but an understanding of the spatiotemporal heterogeneity of changes in the brain is lacking. Here, we investigated the spatially resolved transcriptome in a 5XFAD AD model at different ages to understand regional changes at the molecular level. Spatially resolved transcriptomic data were obtained from 5XFAD AD models and age-matched control mice. Differentially expressed genes were identified using spots clustered by anatomical structures. Gene signatures of activation of microglia and astrocytes were calculated and mapped on the spatially resolved transcriptomic data. We identified early alterations in the white matter (WM) of the AD model before the definite accumulation of amyloid plaques in the gray matter (GM). Changes in the early stage of the disease involved primarily glial cell activation in the WM, whereas the changes in the later stage of pathology were prominent in the GM. We confirmed that disease-associated microglia (DAM) and astrocyte (DAA) signatures also showed initial changes in WM and that activation spreads to GM. Trajectory inference using microglial gene sets revealed the subdivision of DAMs with different spatial patterns. Taken together, these results help to understand the spatiotemporal changes associated with reactive glial cells as a major pathophysiological characteristic of AD. The heterogeneous spatial molecular changes apply to identifying diagnostic and therapeutic targets caused by amyloid accumulation in AD.

Universal method for the isolation of microvessels from frozen brain tissue: A proof-of-concept multiomic investigation of the neurovasculature

The neurovascular unit, comprised of vascular cell types that collectively regulate cerebral blood flow to meet the needs of coupled neurons, is paramount for the proper function of the central nervous system. The neurovascular unit gatekeeps blood-brain barrier properties, which experiences impairment in several central nervous system diseases associated with neuroinflammation and contributes to pathogenesis. To better understand function and dysfunction at the neurovascular unit and how it may confer inflammatory processes within the brain, isolation and characterization of the neurovascular unit is needed. Here, we describe a singular, standardized protocol to enrich and isolate microvessels from archived snap-frozen human and frozen mouse cerebral cortex using mechanical homogenization and centrifugation-separation that preserves the structural integrity and multicellular composition of microvessel fragments. For the first time, microvessels are isolated from postmortem ventromedial prefrontal cortex tissue and are comprehensively investigated as a structural unit using both RNA sequencing and Liquid Chromatography with tandem mass spectrometry (LC-MS/MS). Both the transcriptome and proteome are obtained and compared, demonstrating that the isolated brain microvessel is a robust model for the NVU and can be used to generate highly informative datasets in both physiological and disease contexts.

Microglial Transcriptional Signatures in the Central Nervous System: Toward A Future of Unraveling Their Function in Health and Disease

Microglia, the resident immune cells of the central nervous system (CNS), are primarily derived from the embryonic yolk sac and make their way to the CNS during early development. They play key physiological and immunological roles across the life span, throughout health, injury, and disease. Recent transcriptomic studies have identified gene transcript signatures expressed by microglia that may provide the foundation for unprecedented insights into their functions. Microglial gene expression signatures can help distinguish them from macrophage cell types to a reasonable degree of certainty, depending on the context. Microglial expression patterns further suggest a heterogeneous population comprised of many states that vary according to the spatiotemporal context. Microglial diversity is most pronounced during development, when extensive CNS remodeling takes place, and following disease or injury. A next step of importance for the field will be to identify the functional roles performed by these various microglial states, with the perspective of targeting them therapeutically.

Unraveling Neuronal Identities Using SIMS: A Deep Learning Label Transfer Tool for Single-Cell RNA Sequencing Analysis

Large single-cell RNA datasets have contributed to unprecedented biological insight. Often, these take the form of cell atlases and serve as a reference for automating cell labeling of newly sequenced samples. Yet, classification algorithms have lacked the capacity to accurately annotate cells, particularly in complex datasets. Here we present SIMS (Scalable, Interpretable Machine Learning for Single-Cell), an end-to-end data-efficient machine learning pipeline for discrete classification of single-cell data that can be applied to new datasets with minimal coding. We benchmarked SIMS against common single-cell label transfer tools and demonstrated that it performs as well or better than state of the art algorithms. We then use SIMS to classify cells in one of the most complex tissues: the brain. We show that SIMS classifies cells of the adult cerebral cortex and hippocampus at a remarkably high accuracy. This accuracy is maintained in trans-sample label transfers of the adult human cerebral cortex. We then apply SIMS to classify cells in the developing brain and demonstrate a high level of accuracy at predicting neuronal subtypes, even in periods of fate refinement, shedding light on genetic changes affecting specific cell types across development. Finally, we apply SIMS to single cell datasets of cortical organoids to predict cell identities and unveil genetic variations between cell lines. SIMS identifies cell-line differences and misannotated cell lineages in human cortical organoids derived from different pluripotent stem cell lines. When cell types are obscured by stress signals, label transfer from primary tissue improves the accuracy of cortical organoid annotations, serving as a reliable ground truth. Altogether, we show that SIMS is a versatile and robust tool for cell-type classification from single-cell datasets.

A public resource of single cell transcriptomes and multiscale networks from persons with and without Alzheimer's disease

The emergence of technologies that can support high-throughput profiling of single cell transcriptomes offers to revolutionize the study of brain tissue from persons with and without Alzheimer's disease (AD). Integration of these data with additional complementary multiomics data such as genetics, proteomics and clinical data provides powerful opportunities to link observed cell subpopulations and molecular network features within a broader disease-relevant context. We report here single nucleus RNA sequencing (snRNA-seq) profiles generated from superior frontal gyrus cortical tissue samples from 101 exceptionally well characterized, aged subjects from the Banner Brain and Body Donation Program in combination with whole genome sequences. We report findings that link common AD risk variants with CR1 expression in oligodendrocytes as well as alterations in peripheral hematological lab parameters, with these observations replicated in an independent, prospective cohort study of ageing and dementia. We also observed an AD-associated CD83(+) microglial subtype with unique molecular networks that encompass many known regulators of AD-relevant microglial biology, and which are associated with immunoglobulin IgG4 production in the transverse colon. These findings illustrate the power of multi-tissue molecular profiling to contextualize snRNA-seq brain transcriptomics and reveal novel disease biology. The transcriptomic, genetic, phenotypic, and network data resources described within this study are available for access and utilization by the scientific community.

Assessing Markovian and Delay Models for Single-Nucleus RNA Sequencing

The serial nature of reactions involved in the RNA life-cycle motivates the incorporation of delays in models of transcriptional dynamics. The models couple a transcriptional process to a fairly general set of delayed monomolecular reactions with no feedback. We provide numerical strategies for calculating the RNA copy number distributions induced by these models, and solve several systems with splicing, degradation, and catalysis. An analysis of single-cell and single-nucleus RNA sequencing data using these models reveals that the kinetics of nuclear export do not appear to require invocation of a non-Markovian waiting time.

A facile method to generate cerebral organoids from human pluripotent stem cells

Human cerebral organoids (COs) are self-organizing three-dimensional (3D) neural structures that provide a human-specific platform to study the cellular and molecular processes that underlie different neurological events. The first step of CO generation from human pluripotent stem cells (hPSCs) is neural induction, which is an in vitro simulation of neural ectoderm development. Several signaling pathways cooperate during neural ectoderm development and in vitro differentiation of hPSCs toward neural cell lineages is also affected by them. In this study, we considered some of the known sources of these variable signaling cues arising from cell culture media components and sought to modulate their effects by applying a comprehensive combination of small molecules and growth factors for CO generation. Histological analysis demonstrated that these COs recapitulate the neural progenitor zone and early cortical layer organization, containing different types of neuronal and glial cells which was in accordance with single-nucleus transcriptome profiling results. Moreover, patch clamp and intracellular Ca2+ dynamic studies demonstrated that the COs behave as a functional neural network. Thus, this method serves as a facile protocol for generating hPSC-derived COs that faithfully mimic the features of their in vivo counterparts in the developing human brain. See also Figure 1(Fig. 1).

Glioblastoma initiation, migration, and cell types are regulated by core bHLH transcription factors ASCL1 and OLIG2

Glioblastomas (GBMs) are highly aggressive, infiltrative, and heterogeneous brain tumors driven by complex driver mutations and glioma stem cells (GSCs). The neurodevelopmental transcription factors ASCL1 and OLIG2 are co-expressed in GBMs, but their role in regulating the heterogeneity and hierarchy of GBM tumor cells is unclear. Here, we show that oncogenic driver mutations lead to dysregulation of ASCL1 and OLIG2, which function redundantly to initiate brain tumor formation in a mouse model of GBM. Subsequently, the dynamic levels and reciprocal binding of ASCL1 and OLIG2 to each other and to downstream target genes then determine the cell types and degree of migration of tumor cells. Single-cell RNA sequencing (scRNA-seq) reveals that a high level of ASCL1 is key in defining GSCs by upregulating a collection of ribosomal protein, mitochondrial, neural stem cell (NSC), and cancer metastasis genes - all essential for sustaining the high proliferation, migration, and therapeutic resistance of GSCs.

Single-nucleus transcriptomes reveal spatiotemporal symbiotic perception and early response in Medicago

Establishing legume-rhizobial symbiosis requires precise coordination of complex responses in a time- and cell type-specific manner. Encountering Rhizobium, rapid changes of gene expression levels in host plants occur in the first few hours, which prepare the plants to turn off defence and form a symbiotic relationship with the microbes. Here, we applied single-nucleus RNA sequencing to characterize the roots of Medicago truncatula at 30 min, 6 h and 24 h after nod factor treatment. We found drastic global gene expression reprogramming at 30 min in the epidermis and cortex and most of these changes were restored at 6 h. Moreover, plant defence response genes are activated at 30 min and subsequently suppressed at 6 h in non-meristem cells. Only in the cortical cells but not in other cell types, we found the flavonoid synthase genes required to recruit rhizobia are highly expressed 30 min after inoculation with nod factors. A gene module enriched for symbiotic nitrogen fixation genes showed that MtFER (MtFERONIA) and LYK3 (LysM domain receptor-like kinase 3) share similar responses to symbiotic signals. We further found that MtFER can be phosphorylated by LYK3 and it participates in rhizobial symbiosis. Our results expand our understanding of dynamic spatiotemporal symbiotic responses at the single-cell level.

Single-Nucleus Transcriptome Profiling of Dorsolateral Prefrontal Cortex: Mechanistic Roles for Neuronal Gene Expression, Including the 17q21.31 Locus, in PTSD Stress Response

OBJECTIVE

Multidisciplinary studies of posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) implicate the dorsolateral prefrontal cortex (DLPFC) in disease risk and pathophysiology. Postmortem brain studies have relied on bulk-tissue RNA sequencing (RNA-seq), but single-cell RNA-seq is needed to dissect cell-type-specific mechanisms. The authors conducted the first single-nucleus RNA-seq postmortem brain study in PTSD to elucidate disease transcriptomic pathology with cell-type-specific resolution.

METHOD

Profiling of 32 DLPFC samples from 11 individuals with PTSD, 10 with MDD, and 11 control subjects was conducted (∼415K nuclei; >13K cells per sample). A replication sample included 15 DLPFC samples (∼160K nuclei; >11K cells per sample).

RESULTS

Differential gene expression analyses identified significant single-nucleus RNA-seq differentially expressed genes (snDEGs) in excitatory (EX) and inhibitory (IN) neurons and astrocytes, but not in other cell types or bulk tissue. MDD samples had more false discovery rate-corrected significant snDEGs, and PTSD samples had a greater replication rate. In EX and IN neurons, biological pathways that were differentially enriched in PTSD compared with MDD included glucocorticoid signaling. Furthermore, glucocorticoid signaling in induced pluripotent stem cell (iPSC)-derived cortical neurons demonstrated greater relevance in PTSD and opposite direction of regulation compared with MDD, especially in EX neurons. Many snDEGs were from the 17q21.31 locus and are particularly interesting given causal roles in disease pathogenesis and DLPFC-based neuroimaging (PTSD: ARL17B, LINC02210-CRHR1, and LRRC37A2; MDD: LRRC37A and LRP4), while others were regulated by glucocorticoids in iPSC-derived neurons (PTSD: SLC16A6, TAF1C; MDD: CDH3).

CONCLUSIONS

The study findings point to cell-type-specific mechanisms of brain stress response in PTSD and MDD, highlighting the importance of examining cell-type-specific gene expression and indicating promising novel biomarkers and therapeutic targets.

Bulk brain tissue cell-type deconvolution with bias correction for single-nuclei RNA sequencing data using DeTREM

BACKGROUND

Quantifying cell-type abundance in bulk tissue RNA-sequencing enables researchers to better understand complex systems. Newer deconvolution methodologies, such as MuSiC, use cell-type signatures derived from single-cell RNA-sequencing (scRNA-seq) data to make these calculations. Single-nuclei RNA-sequencing (snRNA-seq) reference data can be used instead of scRNA-seq data for tissues such as human brain where single-cell data are difficult to obtain, but accuracy suffers due to sequencing differences between the technologies.

RESULTS

We propose a modification to MuSiC entitled 'DeTREM' which compensates for sequencing differences between the cell-type signature and bulk RNA-seq datasets in order to better predict cell-type fractions. We show DeTREM to be more accurate than MuSiC in simulated and real human brain bulk RNA-sequencing datasets with various cell-type abundance estimates. We also compare DeTREM to SCDC and CIBERSORTx, two recent deconvolution methods that use scRNA-seq cell-type signatures. We find that they perform well in simulated data but produce less accurate results than DeTREM when used to deconvolute human brain data.

CONCLUSION

DeTREM improves the deconvolution accuracy of MuSiC and outperforms other deconvolution methods when applied to snRNA-seq data. DeTREM enables accurate cell-type deconvolution in situations where scRNA-seq data are not available. This modification improves characterization cell-type specific effects in brain tissue and identification of cell-type abundance differences under various conditions.

Human astrocytes and microglia show augmented ingestion of synapses in Alzheimer's disease via MFG-E8

Synapse loss correlates with cognitive decline in Alzheimer's disease (AD). Data from mouse models suggests microglia are important for synapse degeneration, but direct human evidence for any glial involvement in synapse removal in human AD remains to be established. Here we observe astrocytes and microglia from human brains contain greater amounts of synaptic protein in AD compared with non-disease controls, and that proximity to amyloid-β plaques and the APOE4 risk gene exacerbate this effect. In culture, mouse and human astrocytes and primary mouse and human microglia phagocytose AD patient-derived synapses more than synapses from controls. Inhibiting interactions of MFG-E8 rescues the elevated engulfment of AD synapses by astrocytes and microglia without affecting control synapse uptake. Thus, AD promotes increased synapse ingestion by human glial cells at least in part via an MFG-E8 opsonophagocytic mechanism with potential for targeted therapeutic manipulation.

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.

Molecular consequences of peripheral Influenza A infection on cell populations in the murine hypothalamus

Infection with Influenza A virus (IAV) causes the well-known symptoms of the flu, including fever, loss of appetite, and excessive sleepiness. These responses, mediated by the brain, will normally disappear once the virus is cleared from the system, but a severe respiratory virus infection may cause long-lasting neurological disturbances. These include encephalitis lethargica and narcolepsy. The mechanisms behind such long lasting changes are unknown. The hypothalamus is a central regulator of the homeostatic response during a viral challenge. To gain insight into the neuronal and non-neuronal molecular changes during an IAV infection, we intranasally infected mice with an H1N1 virus and extracted the brain at different time points. Using single-nucleus RNA sequencing (snRNA-seq) of the hypothalamus, we identify transcriptional effects in all identified cell populations. The snRNA-seq data showed the most pronounced transcriptional response at 3 days past infection, with a strong downregulation of genes across all cell types. General immune processes were mainly impacted in microglia, the brain resident immune cells, where we found increased numbers of cells expressing pro-inflammatory gene networks. In addition, we found that most neuronal cell populations downregulated genes contributing to the energy homeostasis in mitochondria and protein translation in the cytosol, indicating potential reduced cellular and neuronal activity. This might be a preventive mechanism in neuronal cells to avoid intracellular viral replication and attack by phagocytosing cells. The change of microglia gene activity suggest that this is complemented by a shift in microglia activity to provide increased surveillance of their surroundings.

Multi-OMIC analysis of Huntington disease reveals a neuroprotective astrocyte state

Huntington disease (HD) is an incurable neurodegenerative disease characterized by neuronal loss and astrogliosis. One hallmark of HD is the selective neuronal vulnerability of striatal medium spiny neurons. To date, the underlying mechanisms of this selective vulnerability have not been fully defined. Here, we employed a multi-omic approach including single nucleus RNAseq (snRNAseq), bulk RNAseq, lipidomics, HTT gene CAG repeat length measurements, and multiplexed immunofluorescence on post-mortem brain tissue from multiple brain regions of HD and control donors. We defined a signature of genes that is driven by CAG repeat length and found it enriched in astrocytic and microglial genes. Moreover, weighted gene correlation network analysis showed loss of connectivity of astrocytic and microglial modules in HD and identified modules that correlated with CAG-repeat length which further implicated inflammatory pathways and metabolism. We performed lipidomic analysis of HD and control brains and identified several lipid species that correlate with HD grade, including ceramides and very long chain fatty acids. Integration of lipidomics and bulk transcriptomics identified a consensus gene signature that correlates with HD grade and HD lipidomic abnormalities and implicated the unfolded protein response pathway. Because astrocytes are critical for brain lipid metabolism and play important roles in regulating inflammation, we analyzed our snRNAseq dataset with an emphasis on astrocyte pathology. We found two main astrocyte types that spanned multiple brain regions; these types correspond to protoplasmic astrocytes, and fibrous-like - CD44-positive, astrocytes. HD pathology was differentially associated with these cell types in a region-specific manner. One protoplasmic astrocyte cluster showed high expression of metallothionein genes, the depletion of this cluster positively correlated with the depletion of vulnerable medium spiny neurons in the caudate nucleus. We confirmed that metallothioneins were increased in cingulate HD astrocytes but were unchanged or even decreased in caudate astrocytes. We combined existing genome-wide association studies (GWAS) with a GWA study conducted on HD patients from the original Venezuelan cohort and identified a single-nucleotide polymorphism in the metallothionein gene locus associated with delayed age of onset. Functional studies found that metallothionein overexpressing astrocytes are better able to buffer glutamate and were neuroprotective of patient-derived directly reprogrammed HD MSNs as well as against rotenone-induced neuronal death in vitro. Finally, we found that metallothionein-overexpressing astrocytes increased the phagocytic activity of microglia in vitro and increased the expression of genes involved in fatty acid binding. Together, we identified an astrocytic phenotype that is regionally-enriched in less vulnerable brain regions that can be leveraged to protect neurons in HD.

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.

Investigating microglia-neuron crosstalk by characterizing microglial contamination in human and mouse patch-seq datasets

Microglia are cells with diverse roles, including the regulation of neuronal excitability. We leveraged Patch-seq to assess the presence and effects of microglia in the local microenvironment of recorded neurons. We first quantified the amounts of microglial transcripts in three Patch-seq datasets of human and mouse neocortical neurons, observing extensive contamination. Variation in microglial contamination was explained foremost by donor identity, particularly in human samples, and additionally by neuronal cell type identity in mice. Gene set enrichment analysis suggests that microglial contamination is reflective of activated microglia, and that these transcriptional signatures are distinct from those captured via single-nucleus RNA-seq. Finally, neurons with greater microglial contamination differed markedly in their electrophysiological characteristics, including lowered input resistances and more depolarized action potential thresholds. Our results generalize beyond Patch-seq to suggest that activated microglia may be widely present across brain slice preparations and contribute to neuron- and donor-related electrophysiological variability in vitro.