NEBULA is a fast negative binomial mixed model for differential or co-expression analysis of large-scale multi-subject single-cell data

Single-nucleus RNA-seq identifies one galanin neuronal subtype in mouse preoptic hypothalamus activated during recovery from sleep deprivation

The preoptic area of the hypothalamus (POA) is essential for sleep regulation. However, the cellular makeup of the POA is heterogeneous, and the molecular identities of the sleep-promoting cells remain elusive. To address this question, this study compares mice during recovery sleep following sleep deprivation to mice allowed extended sleep. Single-nucleus RNA sequencing (single-nucleus RNA-seq) identifies one galanin inhibitory neuronal subtype that shows upregulation of rapid and delayed activity-regulated genes during recovery sleep. This cell type expresses higher levels of growth hormone receptor and lower levels of estrogen receptor compared to other galanin subtypes. single-nucleus RNA-seq also reveals cell-type-specific upregulation of purinergic receptor (P2ry14) and serotonin receptor (Htr2a) during recovery sleep in this neuronal subtype, suggesting possible mechanisms for sleep regulation. Studies with RNAscope validate the single-nucleus RNA-seq findings. Thus, the combined use of single-nucleus RNA-seq and activity-regulated genes identifies a neuronal subtype functionally involved in sleep regulation.

Association of 10 VEGF Family Genes with Alzheimer's Disease Endophenotypes at Single Cell Resolution

The cell-type specific role of the vascular endothelial growth factors (VEGFs) in the pathogenesis of Alzheimer's disease (AD) is not well characterized. In this study, we utilized a single-nucleus RNA sequencing dataset from Dorsolateral Prefrontal Cortex (DLFPC) of 424 donors from the Religious Orders Study and Memory and Aging Project (ROS/MAP) to investigate the effect of 10 VEGF genes ( VEGFA, VEGFB, VEGFC, VEGFD, PGF, FLT1, FLT4, KDR, NRP1 , and NRP2 ) on AD endophenotypes. Mean age of death was 89 years, among which 68% were females, and 52% has AD dementia. Negative binomial mixed models were used for differential expression analysis and for association analysis with β-amyloid load, PHF tau tangle density, and both cross-sectional and longitudinal global cognitive function. Intercellular VEGF-associated signaling was profiled using CellChat. We discovered prefrontal cortical FLT1 expression was upregulated in AD brains in both endothelial and microglial cells. Higher FLT1 expression was also associated with worse cross-sectional global cognitive function, longitudinal cognitive trajectories, and β-amyloid load. Similarly, higher endothelial FLT4 expression was associated with more β-amyloid load. In contrast to the receptors, VEGFB showed opposing effects on β-amyloid load whereby higher levels in oligodendrocytes was associated with high amyloid burden, while higher levels in inhibitory neurons was associated with lower amyloid burden. Finally, AD cells showed significant reduction in overall VEGF signaling comparing to those from cognitive normal participants. Our results highlight key changes in VEGF receptor expression in endothelial and microglial cells during AD, and the potential protective role of VEGFB in neurons.

Modeling type 1 diabetes progression using machine learning and single-cell transcriptomic measurements in human islets

Type 1 diabetes (T1D) is a chronic condition in which beta cells are destroyed by immune cells. Despite progress in immunotherapies that could delay T1D onset, early detection of autoimmunity remains challenging. Here, we evaluate the utility of machine learning for early prediction of T1D using single-cell analysis of islets. Using gradient-boosting algorithms, we model changes in gene expression of single cells from pancreatic tissues in T1D and non-diabetic organ donors. We assess if mathematical modeling could predict the likelihood of T1D development in non-diabetic autoantibody-positive donors. While most autoantibody-positive donors are predicted to be non-diabetic, select donors with unique gene signatures are classified as T1D. Our strategy also reveals a shared gene signature in distinct T1D-associated models across cell types, suggesting a common effect of the disease on transcriptional outputs of these cells. Our study establishes a precedent for using machine learning in early detection of T1D.

Transcriptomic and spatial dissection of human ex vivo right atrial tissue reveals proinflammatory microvascular changes in ischemic heart disease

Cardiovascular disease plays a central role in the electrical and structural remodeling of the right atrium, predisposing to arrhythmias, heart failure, and sudden death. Here, we dissect with single-nuclei RNA sequencing (snRNA-seq) and spatial transcriptomics the gene expression changes in the human ex vivo right atrial tissue and pericardial fluid in ischemic heart disease, myocardial infarction, and ischemic and non-ischemic heart failure using asymptomatic patients with valvular disease who undergo preventive surgery as the control group. We reveal substantial differences in disease-associated gene expression in all cell types, collectively suggesting inflammatory microvascular dysfunction and changes in the right atrial tissue composition as the valvular and vascular diseases progress into heart failure. The data collectively suggest that investigation of human cardiovascular disease should expand to all functionally important parts of the heart, which may help us to identify mechanisms promoting more severe types of the disease.

Recent advances in differential expression analysis for single-cell RNA-seq and spatially resolved transcriptomic studies

Differential expression (DE) analysis is a necessary step in the analysis of single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) data. Unlike traditional bulk RNA-seq, DE analysis for scRNA-seq or SRT data has unique characteristics that may contribute to the difficulty of detecting DE genes. However, the plethora of DE tools that work with various assumptions makes it difficult to choose an appropriate one. Furthermore, a comprehensive review on detecting DE genes for scRNA-seq data or SRT data from multi-condition, multi-sample experimental designs is lacking. To bridge such a gap, here, we first focus on the challenges of DE detection, then highlight potential opportunities that facilitate further progress in scRNA-seq or SRT analysis, and finally provide insights and guidance in selecting appropriate DE tools or developing new computational DE methods.

eSVD-DE: cohort-wide differential expression in single-cell RNA-seq data using exponential-family embeddings

BACKGROUND

Single-cell RNA-sequencing (scRNA) datasets are becoming increasingly popular in clinical and cohort studies, but there is a lack of methods to investigate differentially expressed (DE) genes among such datasets with numerous individuals. While numerous methods exist to find DE genes for scRNA data from limited individuals, differential-expression testing for large cohorts of case and control individuals using scRNA data poses unique challenges due to substantial effects of human variation, i.e., individual-level confounding covariates that are difficult to account for in the presence of sparsely-observed genes.

RESULTS

We develop the eSVD-DE, a matrix factorization that pools information across genes and removes confounding covariate effects, followed by a novel two-sample test in mean expression between case and control individuals. In general, differential testing after dimension reduction yields an inflation of Type-1 errors. However, we overcome this by testing for differences between the case and control individuals' posterior mean distributions via a hierarchical model. In previously published datasets of various biological systems, eSVD-DE has more accuracy and power compared to other DE methods typically repurposed for analyzing cohort-wide differential expression.

CONCLUSIONS

eSVD-DE proposes a novel and powerful way to test for DE genes among cohorts after performing a dimension reduction. Accurate identification of differential expression on the individual level, instead of the cell level, is important for linking scRNA-seq studies to our understanding of the human population.

eSVD-DE: Cohort-wide differential expression in single-cell RNA-seq data using exponential-family embeddings

Background

Single-cell RNA-sequencing (scRNA) datasets are becoming increasingly popular in clinical and cohort studies, but there is a lack of methods to investigate differentially expressed (DE) genes among such datasets with numerous individuals. While numerous methods exist to find DE genes for scRNA data from limited individuals, differential-expression testing for large cohorts of case and control individuals using scRNA data poses unique challenges due to substantial effects of human variation, i.e., individual-level confounding covariates that are difficult to account for in the presence of sparsely-observed genes.

Results

We develop the eSVD-DE, a matrix factorization that pools information across genes and removes confounding covariate effects, followed by a novel two-sample test in mean expression between case and control individuals. In general, differential testing after dimension reduction yields an inflation of Type-1 errors. However, we overcome this by testing for differences between the case and control individuals' posterior mean distributions via a hierarchical model. In previously published datasets of various biological systems, eSVD-DE has more accuracy and power compared to other DE methods typically repurposed for analyzing cohort-wide differential expression.

Conclusions

eSVD-DE proposes a novel and powerful way to test for DE genes among cohorts after performing a dimension reduction. Accurate identification of differential expression on the individual level, instead of the cell level, is important for linking scRNA-seq studies to our understanding of the human population.

Reciprocal activity of AgRP and POMC neurons governs coordinated control of feeding and metabolism

Agouti-related peptide (AgRP)-expressing and proopiomelanocortin (POMC)-expressing neurons reciprocally regulate food intake. Here, we combine non-interacting recombinases to simultaneously express functionally opposing chemogenetic receptors in AgRP and POMC neurons for comparing metabolic responses in male and female mice with simultaneous activation of AgRP and inhibition of POMC neurons with isolated activation of AgRP neurons or isolated inhibition of POMC neurons. We show that food intake is regulated by the additive effect of AgRP neuron activation and POMC neuron inhibition, while systemic insulin sensitivity and gluconeogenesis are differentially modulated by isolated-versus-simultaneous regulation of AgRP and POMC neurons. We identify a neurocircuit engaging Npy1R-expressing neurons in the paraventricular nucleus of the hypothalamus, where activated AgRP neurons and inhibited POMC neurons cooperate to promote food consumption and activate Th+ neurons in the nucleus tractus solitarii. Collectively, these results unveil how food intake is precisely regulated by the simultaneous bidirectional interplay between AgRP and POMC neurocircuits.

Advances in the Application of Single-Cell Transcriptomics in Plant Systems and Synthetic Biology

Plants are complex systems hierarchically organized and composed of various cell types. To understand the molecular underpinnings of complex plant systems, single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for revealing high resolution of gene expression patterns at the cellular level and investigating the cell-type heterogeneity. Furthermore, scRNA-seq analysis of plant biosystems has great potential for generating new knowledge to inform plant biosystems design and synthetic biology, which aims to modify plants genetically/epigenetically through genome editing, engineering, or re-writing based on rational design for increasing crop yield and quality, promoting the bioeconomy and enhancing environmental sustainability. In particular, data from scRNA-seq studies can be utilized to facilitate the development of high-precision Build-Design-Test-Learn capabilities for maximizing the targeted performance of engineered plant biosystems while minimizing unintended side effects. To date, scRNA-seq has been demonstrated in a limited number of plant species, including model plants (e.g., Arabidopsis thaliana), agricultural crops (e.g., Oryza sativa), and bioenergy crops (e.g., Populus spp.). It is expected that future technical advancements will reduce the cost of scRNA-seq and consequently accelerate the application of this emerging technology in plants. In this review, we summarize current technical advancements in plant scRNA-seq, including sample preparation, sequencing, and data analysis, to provide guidance on how to choose the appropriate scRNA-seq methods for different types of plant samples. We then highlight various applications of scRNA-seq in both plant systems biology and plant synthetic biology research. Finally, we discuss the challenges and opportunities for the application of scRNA-seq in plants.

Single-nucleus RNA sequencing demonstrates an autosomal dominant Alzheimer's disease profile and possible mechanisms of disease protection

Highly penetrant autosomal dominant Alzheimer's disease (ADAD) comprises a distinct disease entity as compared to the far more prevalent form of AD in which common variants collectively contribute to risk. The downstream pathways that distinguish these AD forms in specific cell types have not been deeply explored. We compared single-nucleus transcriptomes among a set of 27 cases divided among PSEN1-E280A ADAD carriers, sporadic AD, and controls. Autophagy genes and chaperones clearly defined the PSEN1-E280A cases compared to sporadic AD. Spatial transcriptomics validated the activation of chaperone-mediated autophagy genes in PSEN1-E280A. The PSEN1-E280A case in which much of the brain was spared neurofibrillary pathology and harbored a homozygous APOE3-Christchurch variant revealed possible explanations for protection from AD pathology including overexpression of LRP1 in astrocytes, increased expression of FKBP1B, and decreased PSEN1 expression in neurons. The unique cellular responses in ADAD and sporadic AD require consideration when designing clinical trials.

Single-cell RNA-seq reveals the metabolic status of immune cells response to immunotherapy in triple-negative breast cancer

Immune checkpoint blockade (ICB) therapy offers promise in the treatment of triple-negative breast cancer (TNBC); however, its limited efficacy in certain TNBC patients poses a challenge. In this study, we elucidated the metabolic mechanism at 'sub-subtype' resolution underlying the non-response to ICB therapy in TNBC. Here, an analytic pipeline was developed to reveal the metabolic heterogeneity, which is correlated with the ICB outcomes, within each immune cell subtype. First, we identified metabolic 'sub-subtypes' within certain cell subtypes, predominantly T cell subsets, which are enriched in ICB non-responders and named as non-responder-enriched (NR-E) clusters. Notably, most of NR-E T metabolic cells exhibit globally higher metabolic activities compared to other cells within the same individual subtype. Further, we investigated the extra-cellular signals that trigger the metabolic status of NR-E T cells. In detail, the prediction of cell-to-cell communication indicated that NR-E T cells are regulated by plasmatic dendritic cells (pDCs) through TNFSF9, as well as by macrophages expressing SIGLEC9. In addition, we also validate the communication between TNFSF9+ pDCs and NR-E T cells utilizing deconvolution of spatial transcriptomics analysis. In summary, our research identified specific metabolic 'sub-subtypes' associated with ICB non-response and uncovered the mechanisms of their regulation in TNBC. And the proposed analytical pipeline can be used to examine metabolic heterogeneity within cell types that correlate with diverse phenotypes.

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction

Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1, which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.

A rat liver cell atlas reveals intrahepatic myeloid heterogeneity

The large size and vascular accessibility of the laboratory rat (Rattus norvegicus) make it an ideal hepatic animal model for diseases that require surgical manipulation. Often, the disease susceptibility and outcomes of inflammatory pathologies vary significantly between strains. This study uses single-cell transcriptomics to better understand the complex cellular network of the rat liver, as well as to unravel the cellular and molecular sources of inter-strain hepatic variation. We generated single-cell and single-nucleus transcriptomic maps of the livers of healthy Dark Agouti and Lewis rat strains and developed a factor analysis-based bioinformatics analysis pipeline to study data covariates, such as strain and batch. Using this approach, we discovered transcriptomic variation within the hepatocyte and myeloid populations that underlie distinct cell states between rat strains. This finding will help provide a reference for future investigations on strain-dependent outcomes of surgical experiment models.

Computational single cell oncology: state of the art

Single cell computational analysis has emerged as a powerful tool in the field of oncology, enabling researchers to decipher the complex cellular heterogeneity that characterizes cancer. By leveraging computational algorithms and bioinformatics approaches, this methodology provides insights into the underlying genetic, epigenetic and transcriptomic variations among individual cancer cells. In this paper, we present a comprehensive overview of single cell computational analysis in oncology, discussing the key computational techniques employed for data processing, analysis, and interpretation. We explore the challenges associated with single cell data, including data quality control, normalization, dimensionality reduction, clustering, and trajectory inference. Furthermore, we highlight the applications of single cell computational analysis, including the identification of novel cell states, the characterization of tumor subtypes, the discovery of biomarkers, and the prediction of therapy response. Finally, we address the future directions and potential advancements in the field, including the development of machine learning and deep learning approaches for single cell analysis. Overall, this paper aims to provide a roadmap for researchers interested in leveraging computational methods to unlock the full potential of single cell analysis in understanding cancer biology with the goal of advancing precision oncology. For this purpose, we also include a notebook that instructs on how to apply the recommended tools in the Preprocessing and Quality Control section.

Dissecting the human leptomeninges at single-cell resolution

Emerging evidence shows that the meninges conduct essential immune surveillance and immune defense at the brain border, and the dysfunction of meningeal immunity contributes to aging and neurodegeneration. However, no study exists on the molecular properties of cell types within human leptomeninges. Here, we provide single nuclei profiling of dissected postmortem leptomeninges from aged individuals. We detect diverse cell types, including unique meningeal endothelial, mural, and fibroblast subtypes. For immune cells, we show that most T cells express CD8 and bear characteristics of tissue-resident memory T cells. We also identify distinct subtypes of border-associated macrophages (BAMs) that display differential gene expressions from microglia and express risk genes for Alzheimer's Disease (AD), as nominated by genome-wide association studies (GWAS). We discover cell-type-specific differentially expressed genes in individuals with Alzheimer's dementia, particularly in fibroblasts and BAMs. Indeed, when cultured, leptomeningeal cells display the signature of ex vivo AD fibroblasts upon amyloid-β treatment. We further explore ligand-receptor interactions within the leptomeningeal niche and computationally infer intercellular communications in AD. Thus, our study establishes a molecular map of human leptomeningeal cell types, providing significant insight into the border immune and fibrotic responses in AD.

Negative Binomial Mixture Model for Identification of Noise in Antigen-Specificity Predictions by LIBRA-seq

Motivation

LIBRA-seq (linking B cell receptor to antigen specificity by sequencing) provides a powerful tool for interrogating the antigen-specific B cell compartment and identifying antibodies against antigen targets of interest. Identification of noise in LIBRA-seq antigen count data is critical for improving antigen binding predictions for downstream applications including antibody discovery and machine learning technologies.

Results

In this study, we present a method for denoising LIBRA-seq data by clustering antigen counts into signal and noise components with a negative binomial mixture model. This approach leverages the VRC01 negative control cells included in a recent LIBRA-seq study(Abu-Shmais et al.) to provide a data-driven means for identification of technical noise. We apply this method to a dataset of nine donors representing separate LIBRA-seq experiments and show that our approach provides improved predictions for in vitro antibody-antigen binding when compared to the standard scoring method used in LIBRA-seq, despite variance in data size and noise structure across samples. This development will improve the ability of LIBRA-seq to identify antigen-specific B cells and contribute to providing more reliable datasets for future machine learning based approaches to predicting antibody-antigen binding as the corpus of LIBRA-seq data continues to grow.

4R-Tau seeding activity unravels molecular subtypes in patients with Progressive Supranuclear Palsy

Progressive Supranuclear palsy (PSP) is a 4-repeat (4-R) tauopathy. We hypothesized that the molecular diversity of tau could explain the heterogeneity seen in PSP disease progression. To test this hypothesis, we performed an extensive biochemical characterisation of the high molecular weight tau species (HMW-Tau) in 20 different brain regions of 25 PSP patients. We found a correlation between the HMW-Tau species and tau seeding capacity in the primary motor cortex, where we confirmed that an elevated 4R-Tau seeding activity correlates with a shorter disease duration. To identify factors that contribute to these differences, we performed proteomic and spatial transcriptomic analysis that revealed key mechanistic pathways, in particular those involving the immune system, that defined patients demonstrating high and low tau seeding capacity. These observations suggest that differences in the tau seeding activity may contribute to the considerable heterogeneity seen in disease progression of patients suffering from PSP.

Epigenomic dissection of Alzheimer's disease pinpoints causal variants and reveals epigenome erosion

Recent work has identified dozens of non-coding loci for Alzheimer's disease (AD) risk, but their mechanisms and AD transcriptional regulatory circuitry are poorly understood. Here, we profile epigenomic and transcriptomic landscapes of 850,000 nuclei from prefrontal cortexes of 92 individuals with and without AD to build a map of the brain regulome, including epigenomic profiles, transcriptional regulators, co-accessibility modules, and peak-to-gene links in a cell-type-specific manner. We develop methods for multimodal integration and detecting regulatory modules using peak-to-gene linking. We show AD risk loci are enriched in microglial enhancers and for specific TFs including SPI1, ELF2, and RUNX1. We detect 9,628 cell-type-specific ATAC-QTL loci, which we integrate alongside peak-to-gene links to prioritize AD variant regulatory circuits. We report differential accessibility of regulatory modules in late AD in glia and in early AD in neurons. Strikingly, late-stage AD brains show global epigenome dysregulation indicative of epigenome erosion and cell identity loss.

In vivo screening characterizes chromatin factor functions during normal and malignant hematopoiesis

Cellular differentiation requires extensive alterations in chromatin structure and function, which is elicited by the coordinated action of chromatin and transcription factors. By contrast with transcription factors, the roles of chromatin factors in differentiation have not been systematically characterized. Here, we combine bulk ex vivo and single-cell in vivo CRISPR screens to characterize the role of chromatin factor families in hematopoiesis. We uncover marked lineage specificities for 142 chromatin factors, revealing functional diversity among related chromatin factors (i.e. barrier-to-autointegration factor subcomplexes) as well as shared roles for unrelated repressive complexes that restrain excessive myeloid differentiation. Using epigenetic profiling, we identify functional interactions between lineage-determining transcription factors and several chromatin factors that explain their lineage dependencies. Studying chromatin factor functions in leukemia, we show that leukemia cells engage homeostatic chromatin factor functions to block differentiation, generating specific chromatin factor-transcription factor interactions that might be therapeutically targeted. Together, our work elucidates the lineage-determining properties of chromatin factors across normal and malignant hematopoiesis.

Single-cell transcriptomics and epigenomics unravel the role of monocytes in neuroblastoma bone marrow metastasis

Metastasis is the major cause of cancer-related deaths. Neuroblastoma (NB), a childhood tumor has been molecularly defined at the primary cancer site, however, the bone marrow (BM) as the metastatic niche of NB is poorly characterized. Here we perform single-cell transcriptomic and epigenomic profiling of BM aspirates from 11 subjects spanning three major NB subtypes and compare these to five age-matched and metastasis-free BM, followed by in-depth single cell analyses of tissue diversity and cell-cell interactions, as well as functional validation. We show that cellular plasticity of NB tumor cells is conserved upon metastasis and tumor cell type composition is NB subtype-dependent. NB cells signal to the BM microenvironment, rewiring via macrophage mgration inhibitory factor and midkine signaling specifically monocytes, which exhibit M1 and M2 features, are marked by activation of pro- and anti-inflammatory programs, and express tumor-promoting factors, reminiscent of tumor-associated macrophages. The interactions and pathways characterized in our study provide the basis for therapeutic approaches that target tumor-to-microenvironment interactions.

scViewer: An Interactive Single-Cell Gene Expression Visualization Tool

Single-cell RNA sequencing (scRNA-seq) is an attractive technology for researchers to gain valuable insights into the cellular processes and cell type diversity present in all tissues. The data generated by the scRNA-seq experiment are high-dimensional and complex in nature. Several tools are now available to analyze the raw scRNA-seq data from public databases; however, simple and easy-to-explore single-cell gene expression visualization tools focusing on differential expression and co-expression are lacking. Here, we present scViewer, an interactive graphical user interface (GUI) R/Shiny application designed to facilitate the visualization of scRNA-seq gene expression data. With the processed Seurat RDS object as input, scViewer utilizes several statistical approaches to provide detailed information on the loaded scRNA-seq experiment and generates publication-ready plots. The major functionalities of scViewer include exploring cell-type-specific gene expression, co-expression analysis of two genes, and differential expression analysis with different biological conditions considering both cell-level and subject-level variations using negative binomial mixed modeling. We utilized a publicly available dataset (brain cells from a study of Alzheimer's disease to demonstrate the utility of our tool. scViewer can be downloaded from GitHub as a Shiny app with local installation. Overall, scViewer is a user-friendly application that will allow researchers to visualize and interpret the scRNA-seq data efficiently for multi-condition comparison by performing gene-level differential expression and co-expression analysis on the fly. Considering the functionalities of this Shiny app, scViewer can be a great resource for collaboration between bioinformaticians and wet lab scientists for faster data visualizations.

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.

scRNASequest: an ecosystem of scRNA-seq analysis, visualization, and publishing

BACKGROUND

Single-cell RNA sequencing is a state-of-the-art technology to understand gene expression in complex tissues. With the growing amount of data being generated, the standardization and automation of data analysis are critical to generating hypotheses and discovering biological insights.

RESULTS

Here, we present scRNASequest, a semi-automated single-cell RNA-seq (scRNA-seq) data analysis workflow which allows (1) preprocessing from raw UMI count data, (2) harmonization by one or multiple methods, (3) reference-dataset-based cell type label transfer and embedding projection, (4) multi-sample, multi-condition single-cell level differential gene expression analysis, and (5) seamless integration with cellxgene VIP for visualization and with CellDepot for data hosting and sharing by generating compatible h5ad files.

CONCLUSIONS

We developed scRNASequest, an end-to-end pipeline for single-cell RNA-seq data analysis, visualization, and publishing. The source code under MIT open-source license is provided at https://github.com/interactivereport/scRNASequest . We also prepared a bookdown tutorial for the installation and detailed usage of the pipeline: https://interactivereport.github.io/scRNAsequest/tutorial/docs/ . Users have the option to run it on a local computer with a Linux/Unix system including MacOS, or interact with SGE/Slurm schedulers on high-performance computing (HPC) clusters.

Single-nucleus RNA-sequencing of autosomal dominant Alzheimer disease and risk variant carriers

Genetic studies of Alzheimer disease (AD) have prioritized variants in genes related to the amyloid cascade, lipid metabolism, and neuroimmune modulation. However, the cell-specific effect of variants in these genes is not fully understood. Here, we perform single-nucleus RNA-sequencing (snRNA-seq) on nearly 300,000 nuclei from the parietal cortex of AD autosomal dominant (APP and PSEN1) and risk-modifying variant (APOE, TREM2 and MS4A) carriers. Within individual cell types, we capture genes commonly dysregulated across variant groups. However, specific transcriptional states are more prevalent within variant carriers. TREM2 oligodendrocytes show a dysregulated autophagy-lysosomal pathway, MS4A microglia have dysregulated complement cascade genes, and APOEε4 inhibitory neurons display signs of ferroptosis. All cell types have enriched states in autosomal dominant carriers. We leverage differential expression and single-nucleus ATAC-seq to map GWAS signals to effector cell types including the NCK2 signal to neurons in addition to the initially proposed microglia. Overall, our results provide insights into the transcriptional diversity resulting from AD genetic architecture and cellular heterogeneity. The data can be explored on the online browser ( http://web.hararilab.org/SNARE/ ).

Discovery of Targets for Immune-Metabolic Antitumor Drugs Identifies Estrogen-Related Receptor Alpha

Drugs that kill tumors through multiple mechanisms have the potential for broad clinical benefits. Here, we first developed an in silico multiomics approach (BipotentR) to find cancer cell-specific regulators that simultaneously modulate tumor immunity and another oncogenic pathway and then used it to identify 38 candidate immune-metabolic regulators. We show the tumor activities of these regulators stratify patients with melanoma by their response to anti-PD-1 using machine learning and deep neural approaches, which improve the predictive power of current biomarkers. The topmost identified regulator, ESRRA, is activated in immunotherapy-resistant tumors. Its inhibition killed tumors by suppressing energy metabolism and activating two immune mechanisms: (i) cytokine induction, causing proinflammatory macrophage polarization, and (ii) antigen-presentation stimulation, recruiting CD8+ T cells into tumors. We also demonstrate a wide utility of BipotentR by applying it to angiogenesis and growth suppressor evasion pathways. BipotentR (http://bipotentr.dfci.harvard.edu/) provides a resource for evaluating patient response and discovering drug targets that act simultaneously through multiple mechanisms.

SIGNIFICANCE

BipotentR presents resources for evaluating patient response and identifying targets for drugs that can kill tumors through multiple mechanisms concurrently. Inhibition of the topmost candidate target killed tumors by suppressing energy metabolism and effects on two immune mechanisms. This article is highlighted in the In This Issue feature, p. 517.

Mostly natural sequencing-by-synthesis for scRNA-seq using Ultima sequencing

Here we introduce a mostly natural sequencing-by-synthesis (mnSBS) method for single-cell RNA sequencing (scRNA-seq), adapted to the Ultima genomics platform, and systematically benchmark it against current scRNA-seq technology. mnSBS uses mostly natural, unmodified nucleotides and only a low fraction of fluorescently labeled nucleotides, which allows for high polymerase processivity and lower costs. We demonstrate successful application in four scRNA-seq case studies of different technical and biological types, including 5' and 3' scRNA-seq, human peripheral blood mononuclear cells from a single individual and in multiplex, as well as Perturb-Seq. Benchmarking shows that results from mnSBS-based scRNA-seq are very similar to those using Illumina sequencing, with minor differences in results related to the position of reads relative to annotated gene boundaries, owing to single-end reads of Ultima being closer to gene ends than reads from Illumina. The method is thus compatible with state-of-the-art scRNA-seq libraries independent of the sequencing technology. We expect mnSBS to be of particular utility for cost-effective large-scale scRNA-seq projects.

TREM2-independent microgliosis promotes tau-mediated neurodegeneration in the presence of ApoE4

In addition to tau and Aβ pathologies, inflammation plays an important role in Alzheimer's disease (AD). Variants in APOE and TREM2 increase AD risk. ApoE4 exacerbates tau-linked neurodegeneration and inflammation in P301S tau mice and removal of microglia blocks tau-dependent neurodegeneration. Microglia adopt a heterogeneous population of transcriptomic states in response to pathology, at least some of which are dependent on TREM2. Previously, we reported that knockout (KO) of TREM2 attenuated neurodegeneration in P301S mice that express mouse Apoe. Because of the possible common pathway of ApoE and TREM2 in AD, we tested whether TREM2 KO (T2KO) would block neurodegeneration in P301S Tau mice expressing ApoE4 (TE4), similar to that observed with microglial depletion. Surprisingly, we observed exacerbated neurodegeneration and tau pathology in TE4-T2KO versus TE4 mice, despite decreased TREM2-dependent microgliosis. Our results suggest that tau pathology-dependent microgliosis, that is, TREM2-independent microgliosis, facilitates tau-mediated neurodegeneration in the presence of ApoE4.

Integrated single-cell transcriptomics of cerebrospinal fluid cells in treatment-naïve multiple sclerosis

Multiple sclerosis (MS) is a chronic and often disabling autoimmune disease of the central nervous system (CNS). Cerebrospinal fluid (CSF) surrounds and protects the CNS. Analysis of CSF can aid the diagnosis of CNS diseases, help to identify the prognosis, and underlying mechanisms of diseases. Several recent studies have leveraged single-cell RNA-sequencing (scRNA-seq) to identify MS-associated changes in CSF cells that are considerably more altered than blood cells in MS. However, not all alterations were replicated across all studies. We therefore integrated multiple available scRNA-seq datasets of CSF cells from MS patients with early relapsing-remitting (RRMS) disease. We provide a searchable and interactive resource of this integrated analysis ( https://CSFinMS.bxgenomics.com ) facilitating diverse visualization and analysis methods without requiring computational skills. In the present joint analysis, we replicated the known expansion of B lineage and the recently described expansion of natural killer (NK) cells and some cytotoxic T cells and decrease of monocytes in the CSF in MS. The previous observation of the abundance of Th1-like Th17 effector memory cells in the CSF was not replicated. Expanded CSF B lineage cells resembled class-switched plasmablasts/-cells (e.g., SDC1/CD138, MZB1) as expected. Our integrative analysis thus validates increased cell type diversity and B cell maturation in the CSF in MS and improves accessibility of available data.

Cell-type-specific co-expression inference from single cell RNA-sequencing data

The inference of gene co-expressions from microarray and RNA-sequencing data has led to rich insights on biological processes and disease mechanisms. However, the bulk samples analyzed in most studies are a mixture of different cell types. As a result, the inferred co-expressions are confounded by varying cell type compositions across samples and only offer an aggregated view of gene regulations that may be distinct across different cell types. The advancement of single cell RNA-sequencing (scRNA-seq) technology has enabled the direct inference of co-expressions in specific cell types, facilitating our understanding of cell-type-specific biological functions. However, the high sequencing depth variations and measurement errors in scRNA-seq data present significant challenges in inferring cell-type-specific gene co-expressions, and these issues have not been adequately addressed in the existing methods. We propose a statistical approach, CS-CORE, for estimating and testing cell-type-specific co-expressions, built on a general expression-measurement model that explicitly accounts for sequencing depth variations and measurement errors in the observed single cell data. Systematic evaluations show that most existing methods suffer from inflated false positives and biased co-expression estimates and clustering analysis, whereas CS-CORE has appropriate false positive control, unbiased co-expression estimates, good statistical power and satisfactory performance in downstream co-expression analysis. When applied to analyze scRNA-seq data from postmortem brain samples from Alzheimer’s disease patients and controls and blood samples from COVID-19 patients and controls, CS-CORE identified cell-type-specific co-expressions and differential co-expressions that were more reproducible and/or more enriched for relevant biological pathways than those inferred from other methods.

ApoE in Alzheimer’s disease: pathophysiology and therapeutic strategies

Alzheimer’s disease (AD) is the most common cause of dementia worldwide, and its prevalence is rapidly increasing due to extended lifespans. Among the increasing number of genetic risk factors identified, the apolipoprotein E (APOE) gene remains the strongest and most prevalent, impacting more than half of all AD cases. While the ε4 allele of the APOE gene significantly increases AD risk, the ε2 allele is protective relative to the common ε3 allele. These gene alleles encode three apoE protein isoforms that differ at two amino acid positions. The primary physiological function of apoE is to mediate lipid transport in the brain and periphery; however, additional functions of apoE in diverse biological functions have been recognized. Pathogenically, apoE seeds amyloid-β (Aβ) plaques in the brain with apoE4 driving earlier and more abundant amyloids. ApoE isoforms also have differential effects on multiple Aβ-related or Aβ-independent pathways. The complexity of apoE biology and pathobiology presents challenges to designing effective apoE-targeted therapeutic strategies. This review examines the key pathobiological pathways of apoE and related targeting strategies with a specific focus on the latest technological advances and tools.

Broad transcriptomic dysregulation occurs across the cerebral cortex in ASD

Neuropsychiatric disorders classically lack defining brain pathologies, but recent work has demonstrated dysregulation at the molecular level, characterized by transcriptomic and epigenetic alterations1-3. In autism spectrum disorder (ASD), this molecular pathology involves the upregulation of microglial, astrocyte and neural-immune genes, the downregulation of synaptic genes, and attenuation of gene-expression gradients in cortex1,2,4-6. However, whether these changes are limited to cortical association regions or are more widespread remains unknown. To address this issue, we performed RNA-sequencing analysis of 725 brain samples spanning 11 cortical areas from 112 post-mortem samples from individuals with ASD and neurotypical controls. We find widespread transcriptomic changes across the cortex in ASD, exhibiting an anterior-to-posterior gradient, with the greatest differences in primary visual cortex, coincident with an attenuation of the typical transcriptomic differences between cortical regions. Single-nucleus RNA-sequencing and methylation profiling demonstrate that this robust molecular signature reflects changes in cell-type-specific gene expression, particularly affecting excitatory neurons and glia. Both rare and common ASD-associated genetic variation converge within a downregulated co-expression module involving synaptic signalling, and common variation alone is enriched within a module of upregulated protein chaperone genes. These results highlight widespread molecular changes across the cerebral cortex in ASD, extending beyond association cortex to broadly involve primary sensory regions.

APOE4 impairs myelination via cholesterol dysregulation in oligodendrocytes

APOE4 is the strongest genetic risk factor for Alzheimer's disease1-3. However, the effects of APOE4 on the human brain are not fully understood, limiting opportunities to develop targeted therapeutics for individuals carrying APOE4 and other risk factors for Alzheimer's disease4-8. Here, to gain more comprehensive insights into the impact of APOE4 on the human brain, we performed single-cell transcriptomics profiling of post-mortem human brains from APOE4 carriers compared with non-carriers. This revealed that APOE4 is associated with widespread gene expression changes across all cell types of the human brain. Consistent with the biological function of APOE2-6, APOE4 significantly altered signalling pathways associated with cholesterol homeostasis and transport. Confirming these findings with histological and lipidomic analysis of the post-mortem human brain, induced pluripotent stem-cell-derived cells and targeted-replacement mice, we show that cholesterol is aberrantly deposited in oligodendrocytes-myelinating cells that are responsible for insulating and promoting the electrical activity of neurons. We show that altered cholesterol localization in the APOE4 brain coincides with reduced myelination. Pharmacologically facilitating cholesterol transport increases axonal myelination and improves learning and memory in APOE4 mice. We provide a single-cell atlas describing the transcriptional effects of APOE4 on the aging human brain and establish a functional link between APOE4, cholesterol, myelination and memory, offering therapeutic opportunities for Alzheimer's disease.

Benchmarking methods for detecting differential states between conditions from multi-subject single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) enables researchers to quantify transcriptomes of thousands of cells simultaneously and study transcriptomic changes between cells. scRNA-seq datasets increasingly include multisubject, multicondition experiments to investigate cell-type-specific differential states (DS) between conditions. This can be performed by first identifying the cell types in all the subjects and then by performing a DS analysis between the conditions within each cell type. Naïve single-cell DS analysis methods that treat cells statistically independent are subject to false positives in the presence of variation between biological replicates, an issue known as the pseudoreplicate bias. While several methods have already been introduced to carry out the statistical testing in multisubject scRNA-seq analysis, comparisons that include all these methods are currently lacking. Here, we performed a comprehensive comparison of 18 methods for the identification of DS changes between conditions from multisubject scRNA-seq data. Our results suggest that the pseudobulk methods performed generally best. Both pseudobulks and mixed models that model the subjects as a random effect were superior compared with the naïve single-cell methods that do not model the subjects in any way. While the naïve models achieved higher sensitivity than the pseudobulk methods and the mixed models, they were subject to a high number of false positives. In addition, accounting for subjects through latent variable modeling did not improve the performance of the naïve methods.

Detecting Fear-Memory-Related Genes from Neuronal scRNA-seq Data by Diverse Distributions and Bhattacharyya Distance

The detection of differentially expressed genes (DEGs) is one of most important computational challenges in the analysis of single-cell RNA sequencing (scRNA-seq) data. However, due to the high heterogeneity and dropout noise inherent in scRNAseq data, challenges in detecting DEGs exist when using a single distribution of gene expression levels, leaving much room to improve the precision and robustness of current DEG detection methods. Here, we propose the use of a new method, DEGman, which utilizes several possible diverse distributions in combination with Bhattacharyya distance. DEGman can automatically select the best-fitting distributions of gene expression levels, and then detect DEGs by permutation testing of Bhattacharyya distances of the selected distributions from two cell groups. Compared with several popular DEG analysis tools on both large-scale simulation data and real scRNA-seq data, DEGman shows an overall improvement in the balance of sensitivity and precision. We applied DEGman to scRNA-seq data of TRAP; Ai14 mouse neurons to detect fear-memory-related genes that are significantly differentially expressed in neurons with and without fear memory. DEGman detected well-known fear-memory-related genes and many novel candidates. Interestingly, we found 25 DEGs in common in five neuron clusters that are functionally enriched for synaptic vesicles, indicating that the coupled dynamics of synaptic vesicles across in neurons plays a critical role in remote memory formation. The proposed method leverages the advantage of the use of diverse distributions in DEG analysis, exhibiting better performance in analyzing composite scRNA-seq datasets in real applications.

Pyramidal neuron subtype diversity governs microglia states in the neocortex

Microglia are specialized macrophages in the brain parenchyma that exist in multiple transcriptional states and reside within a wide range of neuronal environments1-4. However, how and where these states are generated remains poorly understood. Here, using the mouse somatosensory cortex, we demonstrate that microglia density and molecular state acquisition are determined by the local composition of pyramidal neuron classes. Using single-cell and spatial transcriptomic profiling, we unveil the molecular signatures and spatial distributions of diverse microglia populations and show that certain states are enriched in specific cortical layers, whereas others are broadly distributed throughout the cortex. Notably, conversion of deep-layer pyramidal neurons to an alternate class identity reconfigures the distribution of local, layer-enriched homeostatic microglia to match the new neuronal niche. Leveraging the transcriptional diversity of pyramidal neurons in the neocortex, we construct a ligand-receptor atlas describing interactions between individual pyramidal neuron subtypes and microglia states, revealing rules of neuron-microglia communication. Our findings uncover a fundamental role for neuronal diversity in instructing the acquisition of microglia states as a potential mechanism for fine-tuning neuroimmune interactions within the cortical local circuitry.

Recommendations of scRNA-seq Differential Gene Expression Analysis Based on Comprehensive Benchmarking

To guide analysts to select the right tool and parameters in differential gene expression analyses of single-cell RNA sequencing (scRNA-seq) data, we developed a novel simulator that recapitulates the data characteristics of real scRNA-seq datasets while accounting for all the relevant sources of variation in a multi-subject, multi-condition scRNA-seq experiment: the cell-to-cell variation within a subject, the variation across subjects, the variability across cell types, the mean/variance relationship of gene expression across genes, library size effects, group effects, and covariate effects. By applying it to benchmark 12 differential gene expression analysis methods (including cell-level and pseudo-bulk methods) on simulated multi-condition, multi-subject data of the 10x Genomics platform, we demonstrated that methods originating from the negative binomial mixed model such as glmmTMB and NEBULA-HL outperformed other methods. Utilizing NEBULA-HL in a statistical analysis pipeline for single-cell analysis will enable scientists to better understand the cell-type-specific transcriptomic response to disease or treatment effects and to discover new drug targets. Further, application to two real datasets showed the outperformance of our differential expression (DE) pipeline, with unified findings of differentially expressed genes (DEG) and a pseudo-time trajectory transcriptomic result. In the end, we made recommendations for filtering strategies of cells and genes based on simulation results to achieve optimal experimental goals.

Allele-specific analysis reveals exon- and cell-type-specific regulatory effects of Alzheimer's disease-associated genetic variants

Elucidating regulatory effects of Alzheimer's disease (AD)-associated genetic variants is critical for unraveling their causal pathways and understanding the pathology. However, their cell-type-specific regulatory mechanisms in the brain remain largely unclear. Here, we conducted an analysis of allele-specific expression quantitative trait loci (aseQTLs) for 33 AD-associated variants in four brain regions and seven cell types using ~3000 bulk RNA-seq samples and >0.25 million single nuclei. We first develop a flexible hierarchical Poisson mixed model (HPMM) and demonstrate its superior statistical power to a beta-binomial model achieved by unifying samples in both allelic and genotype-level expression data. Using the HPMM, we identified 24 (~73%) aseQTLs in at least one brain region, including three new eQTLs associated with CA12, CHRNE, and CASS4. Notably, the APOE ε4 variant reduces APOE expression across all regions, even in AD-unaffected controls. Our results reveal region-dependent and exon-specific effects of multiple aseQTLs, such as rs2093760 with CR1, rs7982 with CLU, and rs3865444 with CD33. In an attempt to pinpoint the cell types responsible for the observed tissue-level aseQTLs using the snRNA-seq data, we detected many aseQTLs in microglia or monocytes associated with immune-related genes, including HLA-DQB1, HLA-DQA2, CD33, FCER1G, MS4A6A, SPI1, and BIN1, highlighting the regulatory role of AD-associated variants in the immune response. These findings provide further insights into potential causal pathways and cell types mediating the effects of the AD-associated variants.

A systematic comparison of novel and existing differential analysis methods for CyTOF data

Cytometry techniques are widely used to discover cellular characteristics at single-cell resolution. Many data analysis methods for cytometry data focus solely on identifying subpopulations via clustering and testing for differential cell abundance. For differential expression analysis of markers between conditions, only few tools exist. These tools either reduce the data distribution to medians, discarding valuable information, or have underlying assumptions that may not hold for all expression patterns. Here, we systematically evaluated existing and novel approaches for differential expression analysis on real and simulated CyTOF data. We found that methods using median marker expressions compute fast and reliable results when the data are not strongly zero-inflated. Methods using all data detect changes in strongly zero-inflated markers, but partially suffer from overprediction or cannot handle big datasets. We present a new method, CyEMD, based on calculating the earth mover's distance between expression distributions that can handle strong zero-inflation without being too sensitive. Additionally, we developed CYANUS - CYtometry ANalysis Using Shiny - a user-friendly R Shiny App allowing the user to analyze cytometry data with state-of-the-art tools, including well-performing methods from our comparison. A public web interface is available at https://exbio.wzw.tum.de/cyanus/.