Search

Abstract Acetylcholine regulates various cognitive functions through broad cholinergic innervation. However, specific cholinergic subpopulations, circuits and molecular mechanisms underlying recognition memory remain largely unknown. Here we show that Ngfr + cholinergic neurons in the substantia innominate (SI)/nucleus basalis of Meynert (nBM)-medial prefrontal cortex (mPFC) circuit selectively underlies recency judgements. Loss of nerve growth factor receptor (Ngfr −/− mice) reduced the excitability of cholinergic neurons in the SI/nBM-mPFC circuit but not in the medial septum (MS)-hippocampus pathway, and impaired temporal order memory but not novel object and object location recognition. Expression of Ngfr in Ngfr −/− SI/nBM restored defected temporal order memory. Fiber photometry revealed that acetylcholine release in mPFC not only predicted object encounters but also mediated recency judgments of objects, and such acetylcholine release was absent in Ngfr −/− mPFC. Chemogenetic and optogenetic inhibition of SI/nBM projection to mPFC in ChAT-Cre mice diminished mPFC acetylcholine release and deteriorated temporal order recognition. Impaired cholinergic activity led to a depolarizing shift of GABAergic inputs to mPFC pyramidal neurons, due to disturbed KCC2-mediated chloride gradients. Finally, potentiation of acetylcholine signaling upregulated KCC2 levels, restored GABAergic driving force and rescued temporal order recognition deficits in Ngfr −/− mice. Thus, NGFR-dependent SI/nBM-mPFC cholinergic circuit underlies temporal order recognition memory.

Abstract Background The APOE gene is the strongest genetic risk factor for late-onset Alzheimer’s Disease (LOAD). However, the gene regulatory mechanisms at this locus remain incompletely characterized. Methods To identify novel AD-linked functional elements within the APOE locus, we integrated SNP variants with multi-omics data from human postmortem brains including 2,179 RNA-seq samples from 3 brain regions and two ancestries (European and African), 667 DNA methylation samples, and ChIP-seq samples. Additionally, we plotted the expression trajectory of APOE transcripts in human brains during development. Results We identified an AD-linked APOE transcript (jxn1.2.2) particularly observed in the dorsolateral prefrontal cortex (DLPFC). The APOE jxn1.2.2 transcript is associated with brain neuropathological features, cognitive impairment, and the presence of the APOE4 allele in DLPFC. We prioritized two independent functional SNPs (rs157580 and rs439401) significantly associated with jxn1.2.2 transcript abundance and DNA methylation levels. These SNPs are located within active chromatin regions and affect brain-related transcription factor-binding affinities. The two SNPs shared effects on the jxn1.2.2 transcript between European and African ethnic groups. Conclusion The novel APOE functional elements provide potential therapeutic targets with mechanistic insight into the disease etiology.

Abstract When somatic cells acquire complex karyotypes, they often are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, “CNV neurons”), but previous analyses of CNV neurons are limited by relatively small sample sizes. Here, we develop an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We apply this approach to 2,125 frontal cortical neurons from a neurotypical human brain. SCOVAL identifies 226 CNV neurons, which include a subclass of 65 CNV neurons with highly aberrant karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we find that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contain fewer, but longer, genes.

Abstract Schizophrenia is a complex neuropsychiatric disorder with sexually dimorphic features, including differential symptomatology, drug responsiveness, and male incidence rate. Prior large-scale transcriptome analyses for sex differences in schizophrenia have focused on the prefrontal cortex. Analyzing BrainSeq Consortium data (caudate nucleus: n = 399, dorsolateral prefrontal cortex: n = 377, and hippocampus: n = 394), we identified 831 unique genes that exhibit sex differences across brain regions, enriched for immune-related pathways. We observed X-chromosome dosage reduction in the hippocampus of male individuals with schizophrenia. Our sex interaction model revealed 148 junctions dysregulated in a sex-specific manner in schizophrenia. Sex-specific schizophrenia analysis identified dozens of differentially expressed genes, notably enriched in immune-related pathways. Finally, our sex-interacting expression quantitative trait loci analysis revealed 704 unique genes, nine associated with schizophrenia risk. These findings emphasize the importance of sex-informed analysis of sexually dimorphic traits, inform personalized therapeutic strategies in schizophrenia, and highlight the need for increased female samples for schizophrenia analyses.

Abstract The polygenic architecture of schizophrenia implicates several molecular pathways involved in synaptic function. However, it is unclear how polygenic risk funnels through these pathways to translate into syndromic illness. Using tensor decomposition, we analyze gene co-expression in the caudate nucleus, hippocampus, and dorsolateral prefrontal cortex of post-mortem brain samples from 358 individuals. We identify a set of genes predominantly expressed in the caudate nucleus and associated with both clinical state and genetic risk for schizophrenia that shows dopaminergic selectivity. A higher polygenic risk score for schizophrenia parsed by this set of genes predicts greater dopamine synthesis in the striatum and greater striatal activation during reward anticipation. These results translate dopamine-linked genetic risk variation into in vivo neurochemical and hemodynamic phenotypes in the striatum that have long been implicated in the pathophysiology of schizophrenia.

Abstract Primary human trophoblast stem cells (TSCs) and TSCs derived from human pluripotent stem cells (hPSCs) can potentially model placental processes in vitro. Yet, the pluripotent states and factors involved in the differentiation of hPSCs to TSCs remain poorly understood. In this study, we demonstrate that the primed pluripotent state can generate TSCs by activating pathways such as Epidermal Growth Factor (EGF) and Wingless-related integration site (WNT), and by suppressing tumor growth factor beta (TGFβ), histone deacetylases (HDAC), and Rho-associated protein kinase (ROCK) signaling pathways, all without the addition of exogenous Bone morphogenetic protein 4 (BMP4)—a condition we refer to as the TS condition. We characterized this process using temporal single-cell RNA sequencing to compare TS conditions with differentiation protocols involving BMP4 activation alone or BMP4 activation in conjunction with WNT inhibition. The TS condition consistently produced a stable, proliferative cell type that closely mimics first-trimester placental cytotrophoblasts, marked by the activation of endogenous retroviral genes and the absence of amnion expression. This was observed across multiple cell lines, including various primed induced pluripotent stem cell (iPSC) and embryonic stem cell (ESC) lines. Primed-derived TSCs can proliferate for over 30 passages and further specify into multinucleated syncytiotrophoblasts and extravillous trophoblast cells. Our research establishes that the differentiation of primed hPSCs to TSC under TS conditions triggers the induction of TMSB4X, BMP5/7, GATA3, and TFAP2A without progressing through a naive state. These findings propose that the primed hPSC state is part of a continuum of potency with the capacity to differentiate into TSCs through multiple routes.

Abstract Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

Abstract Our earlier work has shown that genomic risk for schizophrenia converges with early life complications in affecting risk for the disorder and sex-biased neurodevelopmental trajectories. Here, we identify specific genes and potential mechanisms that, in placenta, may mediate such outcomes. We performed TWAS in healthy term placentae (N = 147) to derive candidate placental causal genes that we confirmed with SMR; to search for placenta and schizophrenia-specific associations, we performed an analogous analysis in fetal brain (N = 166) and additional placenta TWAS for other disorders/traits. The analyses in the whole sample and stratifying by sex ultimately highlight 139 placenta and schizophrenia-specific risk genes, many being sex-biased; the candidate molecular mechanisms converge on the nutrient-sensing capabilities of placenta and trophoblast invasiveness. These genes also implicate the Coronavirus-pathogenesis pathway and showed increased expression in placentae from a small sample of SARS-CoV-2-positive pregnancies. Investigating placental risk genes for schizophrenia and candidate mechanisms may lead to opportunities for prevention that would not be suggested by study of the brain alone.

Summary: Genetics have nominated many schizophrenia risk genes and identified convergent signals between schizophrenia and neurodevelopmental disorders. However, functional interpretation of the nominated genes in the relevant brain cell types is often lacking. We executed interaction proteomics for six schizophrenia risk genes that have also been implicated in neurodevelopment in human induced cortical neurons. The resulting protein network is enriched for common variant risk of schizophrenia in Europeans and East Asians, is down-regulated in layer 5/6 cortical neurons of individuals affected by schizophrenia, and can complement fine-mapping and eQTL data to prioritize additional genes in GWAS loci. A sub-network centered on HCN1 is enriched for common variant risk and contains proteins (HCN4 and AKAP11) enriched for rare protein-truncating mutations in individuals with schizophrenia and bipolar disorder. Our findings showcase brain cell-type-specific interactomes as an organizing framework to facilitate interpretation of genetic and transcriptomic data in schizophrenia and its related disorders.

Abstract Urbanization is increasing globally, and is associated with stress and increased mental health risks, including for depression. However, it remains unclear, especially at the level of brain function, how urbanicity, social threat stressors, and psychiatric risk may be linked. Here, we aim to define the structural and functional MRI neural correlates of social stress, childhood urbanicity, and their putative mechanistic relevance to depressive illness risk, in terms of behavioral traits and genetics. We studied a sample of healthy adults with divergent urban and rural childhoods. We examined childhood urbanicity effects on brain structure as suggested by MRI, and its functional relevance to depression risk, through interactions between urbanicity and trait anxiety-depression, as well as between urbanicity and polygenic risk for depression, during stress-related medial prefrontal cortex (mPFC) engagement. Subjects with divergent rural and urban childhoods were similar in adult socioeconomic status and were genetically homogeneous. Urban childhood was associated with relatively reduced mPFC gray matter volumes as suggested by MRI. MPFC engagement under social status threat correlated with the higher trait anxiety-depression in subjects with urban childhoods, but not in their rural counterparts, implicating an exaggerated physiological response to the threat context with urbanicity, in association with behavioral risk for depression. Stress-associated mPFC engagement also interacted with polygenic risk for depression, significantly predicting a differential mPFC response in individuals with urban but not rural childhoods. Developmental urbanicity, therefore, appears to interact with genetic and behavioral risk for depression on the mPFC neural response to a threat context.

The authors provide a comprehensive, single base resolution view of association between genetic variation and DNA methylation in human brain. They also show that heritability attributed to schizophrenia GWAS-associated variants reflects the epigenetic plasticity of the brain.

Abstract Background Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells. Results Here, the Brain Somatic Mosaicism Network conducts a coordinated, multi-institutional study to examine the ability of existing methods to detect simulated somatic single-nucleotide variants (SNVs) in DNA mixing experiments, generate multiple replicates of whole-genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual, devise strategies to discover somatic SNVs, and apply various approaches to validate somatic SNVs. These efforts lead to the identification of 43 bona fide somatic SNVs that range in variant allele fractions from ~ 0.005 to ~ 0.28. Guided by these results, we devise best practices for calling mosaic SNVs from 250× whole-genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrate that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees. Conclusions This study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

DNA methylation (DNAm) is a key epigenetic regulator of gene expression across development. The developing prenatal brain is a highly dynamic tissue, but our understanding of key drivers of epigenetic variability across development is limited. We, therefore, assessed genomic methylation at over 39 million sites in the prenatal cortex using whole-genome bisulfite sequencing and found loci and regions in which methylation levels are dynamic across development. We saw that DNAm at these loci was associated with nearby gene expression and enriched for enhancer chromatin states in prenatal brain tissue. Additionally, these loci were enriched for genes associated with neuropsychiatric disorders and genes involved with neurogenesis. We also found autosomal differences in DNAm between the sexes during prenatal development, though these have less clear functional consequences. We lastly confirmed that the dynamic methylation at this critical period is specifically CpG methylation, with generally low levels of CpH methylation. Our findings provide detailed insight into prenatal brain development as well as clues to the pathogenesis of psychiatric traits seen later in life.

Here, authors present results of a hiPSC transcriptomics study on corticogenesis from multiple donors across four transitions in differentiation. They present a bulk data deconvolution method and show that co-culturing human NPCs with rodent astrocytes results in mutually synergistic maturation.

We propose a novel optimization framework that integrates imaging and genetics data for simultaneous biomarker identification and disease classification. The generative component of our model uses a dictionary learning framework to project the imaging and genetic data into a shared low dimensional space. We have coupled both the data modalities by tying the linear projection coefficients to the same latent space. The discriminative component of our model uses logistic regression on the projection vectors for disease diagnosis. This prediction task implicitly guides our framework to find interpretable biomarkers that are substantially different between a healthy and disease population. We exploit the interconnectedness of different brain regions by incorporating a graph regularization penalty into the joint objective function. We also use a group sparsity penalty to find a representative set of genetic basis vectors that span a low dimensional space where subjects are easily separable between patients and controls. We have evaluated our model on a population study of schizophrenia that includes two task fMRI paradigms and single nucleotide polymorphism (SNP) data. Using ten-fold cross validation, we compare our generative-discriminative framework with canonical correlation analysis (CCA) of imaging and genetics data, parallel independent component analysis (pICA) of imaging and genetics data, random forest (RF) classification, and a linear support vector machine (SVM). We also quantify the reproducibility of the imaging and genetics biomarkers via subsampling. Our framework achieves higher class prediction accuracy and identifies robust biomarkers. Moreover, the implicated brain regions and genetic variants underlie the well documented deficits in schizophrenia.

Abstract Background DNA methylation (DNAm) is a critical regulator of both development and cellular identity and shows unique patterns in neurons. To better characterize maturational changes in DNAm patterns in these cells, we profile the DNAm landscape at single-base resolution across the first two decades of human neocortical development in NeuN+ neurons using whole-genome bisulfite sequencing and compare them to non-neurons (primarily glia) and prenatal homogenate cortex. Results We show that DNAm changes more dramatically during the first 5 years of postnatal life than during the entire remaining period. We further refine global patterns of increasingly divergent neuronal CpG and CpH methylation (mCpG and mCpH) into six developmental trajectories and find that in contrast to genome-wide patterns, neighboring mCpG and mCpH levels within these regions are highly correlated. We integrate paired RNA-seq data and identify putative regulation of hundreds of transcripts and their splicing events exclusively by mCpH levels, independently from mCpG levels, across this period. We finally explore the relationship between DNAm patterns and development of brain-related phenotypes and find enriched heritability for many phenotypes within identified DNAm features. Conclusions By profiling DNAm changes in NeuN-sorted neurons over the span of human cortical development, we identify novel, dynamic regions of DNAm that would be masked in homogenate DNAm data; expand on the relationship between CpG methylation, CpH methylation, and gene expression; and find enrichment particularly for neuropsychiatric diseases in genomic regions with cell type-specific, developmentally dynamic DNAm patterns.

Abstract Background RNA sequencing offers advantages over other quantification methods for microRNA (miRNA), yet numerous biases make reliable quantification challenging. Previous evaluations of these biases have focused on adapter ligation bias with limited evaluation of reverse transcription bias or amplification bias. Furthermore, evaluations of the quantification of isomiRs (miRNA isoforms) or the influence of starting amount on performance have been very limited. No study had yet evaluated the quantification of isomiRs of altered length or compared the consistency of results derived from multiple moderate starting inputs. We therefore evaluated quantifications of miRNA and isomiRs using four library preparation kits, with various starting amounts, as well as quantifications following removal of duplicate reads using unique molecular identifiers (UMIs) to mitigate reverse transcription and amplification biases. Results All methods resulted in false isomiR detection; however, the adapter-free method tested was especially prone to false isomiR detection. We demonstrate that using UMIs improves accuracy and we provide a guide for input amounts to improve consistency. Conclusions Our data show differences and limitations of current methods, thus raising concerns about the validity of quantification of miRNA and isomiRs across studies. We advocate for the use of UMIs to improve accuracy and reliability of miRNA quantifications.

Christina M. Lill, who contributed to analysis of data, was inadvertently omitted from the author list in the originally published version of this article. This has now been corrected in both the PDF and HTML versions of the article.

Summary: Better understanding of the progression of neural stem cells (NSCs) in the developing cerebral cortex is important for modeling neurogenesis and defining the pathogenesis of neuropsychiatric disorders. Here, we use RNA sequencing, cell imaging, and lineage tracing of mouse and human in vitro NSCs and monkey brain sections to model the generation of cortical neuronal fates. We show that conserved signaling mechanisms regulate the acute transition from proliferative NSCs to committed glutamatergic excitatory neurons. As human telencephalic NSCs develop from pluripotency in vitro, they transition through organizer states that spatially pattern the cortex before generating glutamatergic precursor fates. NSCs derived from multiple human pluripotent lines vary in these early patterning states, leading differentially to dorsal or ventral telencephalic fates. This work furthers systematic analyses of the earliest patterning events that generate the major neuronal trajectories of the human telencephalon.

Patients with schizophrenia show varied response to antipsychotics. Here, the authors demonstrate in patients under antipsychotics treatment that a haplotype associated with lower dysbindin-1 expression correlated with better executive functions, providing further mechanistic support from mouse models.

Cognitive function is associated with health and important life outcomes. Here, the authors perform a genome-wide association study for general cognitive function in 300,486 individuals and identify genetic loci that implicate neural and cell developmental pathways in this trait.

Here, we present a large (n = 107,207) genome-wide association study (GWAS) of general cognitive ability (“g”), further enhanced by combining results with a large-scale GWAS of educational attainment. We identified 70 independent genomic loci associated with general cognitive ability. Results showed significant enrichment for genes causing Mendelian disorders with an intellectual disability phenotype. Competitive pathway analysis implicated the biological processes of neurogenesis and synaptic regulation, as well as the gene targets of two pharmacologic agents: cinnarizine, a T-type calcium channel blocker, and LY97241, a potassium channel inhibitor. Transcriptome-wide and epigenome-wide analysis revealed that the implicated loci were enriched for genes expressed across all brain regions (most strongly in the cerebellum). Enrichment was exclusive to genes expressed in neurons but not oligodendrocytes or astrocytes. Finally, we report genetic correlations between cognitive ability and disparate phenotypes including psychiatric disorders, several autoimmune disorders, longevity, and maternal age at first birth.

Abstract Before their disappearance from the fossil record approximately 40,000 years ago, Neanderthals, the ancient hominin lineage most closely related to modern humans, interbred with ancestors of present-day humans. The legacy of this gene flow persists through Neanderthal-derived variants that survive in modern human DNA; however, the neural implications of this inheritance are uncertain. Here, using MRI in a large cohort of healthy individuals of European-descent, we show that the amount of Neanderthal-originating polymorphism carried in living humans is related to cranial and brain morphology. First, as a validation of our approach, we demonstrate that a greater load of Neanderthal-derived genetic variants (higher “NeanderScore”) is associated with skull shapes resembling those of known Neanderthal cranial remains, particularly in occipital and parietal bones. Next, we demonstrate convergent NeanderScore-related findings in the brain (measured by gray- and white-matter volume, sulcal depth, and gyrification index) that localize to the visual cortex and intraparietal sulcus. This work provides insights into ancestral human neurobiology and suggests that Neanderthal-derived genetic variation is neurologically functional in the contemporary population.

The hippocampus in mammalian brain varies in size across individuals. Here, Hibar and colleagues perform a genome-wide association meta-analysis to find six genetic loci with significant association to hippocampus volume.

Human genetic studies have recently suggested that the postsynaptic activity-regulated cytoskeleton-associated protein (Arc) complex is a convergence signal for several genes implicated in schizophrenia. However, the functional significance of Arc in schizophrenia-related neurobehavioral phenotypes and brain circuits is unclear. Here, we find that, consistent with schizophrenia-related phenotypes, disruption of Arc in mice produces deficits in sensorimotor gating, cognitive functions, social behaviors, and amphetamine-induced psychomotor responses. Furthermore, genetic disruption of Arc leads to concomitant hypoactive mesocortical and hyperactive mesostriatal dopamine pathways. Application of a D1 agonist to the prefrontal cortex or a D2 antagonist in the ventral striatum rescues Arc-dependent cognitive or psychomotor abnormalities, respectively. Our findings demonstrate a role for Arc in the regulation of dopaminergic neurotransmission and related behaviors. The results also provide initial biological support implicating Arc in dopaminergic and behavioral abnormalities related to schizophrenia.

Pharmacological studies indicate that dopamine D1-like receptors (D1 and D5) are critically involved in cognitive function. However, the lack of pharmacological ligands selective for either the D1 or D5 receptors has made it difficult to determine the unique contributions of the D1-like family members. To circumvent these pharmacological limitations, we used D5 receptor homozygous (-/-) and heterozygous (+/-) knockout mice, to identify the specific role of this receptor in higher order cognitive functions. We identified a novel role for D5 receptors in the regulation of spatial working memory and temporal order memory function. The D5 mutant mice acquired a discrete paired-trial variable-delay T-maze task at normal rates. However, both D5+/- and D5-/- mice exhibited impaired performance compared to D5+/+ littermates when a higher burden on working memory faculties was imposed. In a temporal order object recognition task, D5+/- exhibited significant memory deficits. No D5-dependent differences in locomotor functions and interest in exploring objects were evident. Molecular biomarkers of dopaminergic functions within the prefrontal cortex (PFC) revealed a selective gene-dose effect on Akt phosphorylation at Ser473 with increased levels in D5-/- knockout mice. A trend toward reduced levels in CaMKKbeta brain-specific band (64 kDa) in D5-/- compared to D5+/+ was also evident. These findings highlight a previously unidentified role for D5 receptors in working memory function and associated molecular signatures within the PFC.

In the original version of this Article, references in the Methods section incorrectly referred to references in the Supplementary References section. The relevant references (now numbered 20, 27, 42, 47, 69–80) have been removed from the Supplementary References section of the Supplementary Information file and added to the References section of the main manuscript, in both the PDF and HTML versions of the Article.

During the anticipation of task demands frontal control is involved in the assembly of stimulus-response mappings based on current goals. It is not clear whether prefrontal modulations occur in higher-order cortical regions, likely reflecting cognitive anticipation processes. The goal of this paper was to investigate prefrontal modulation during anticipation of upcoming working memory demands as revealed by magnetoencephalography (MEG). Twenty healthy volunteers underwent MEG while they performed a variation of the Sternberg Working Memory (WM) task. Beta band (14–30 Hz) SAM (Synthetic Aperture Magnetometry) analysis was performed. During the preparatory periods there was an increase in beta power (event-related synchronization) in dorsolateral prefrontal cortex (DLPFC) bilaterally, left inferior prefrontal gyrus, left parietal, and temporal areas. Our results provide support for the hypothesis that, during preparatory states, the prefrontal cortex is important for biasing higher order brain regions that are going to be engaged in the upcoming task.

Background: Statistical deconvolution strategies have emerged over the past decade to estimate the proportion of various cell populations in homogenate tissue sources like brain using gene expression data. However, no study has been undertaken to assess the extent to which expression-based and DNAm-based cell type composition estimates agree. Results: Using estimated neuronal fractions from DNAm data, from the same brain region (i.e., matched) as our bulk RNA-Seq dataset, as proxies for the true unobserved cell-type fractions (i.e., as the gold standard), we assessed the accuracy (RMSE) and concordance (R 2) of four reference-based deconvolution algorithms: Houseman, CIBERSORT, non-negative least squares (NNLS)/MIND, and MuSiC. We did this for two cell-type populations - neurons and non-neurons/glia - using matched single nuclei RNA-Seq and mismatched single cell RNA-Seq reference datasets. With the mismatched single cell RNA-Seq reference dataset, Houseman, MuSiC, and NNLS produced concordant (high correlation; Houseman R 2 = 0.51, 95% CI [0.39, 0.65]; MuSiC R 2 = 0.56, 95% CI [0.43, 0.69]; NNLS R 2 = 0.54, 95% CI [0.32, 0.68]) but biased (high RMSE, >0.35) neuronal fraction estimates. CIBERSORT produced more discordant (moderate correlation; R 2 = 0.25, 95% CI [0.15, 0.38]) neuronal fraction estimates, but with less bias (low RSME, 0.09). Using the matched single nuclei RNA-Seq reference dataset did not eliminate bias (MuSiC RMSE = 0.17). Conclusions: Our results together suggest that many existing RNA deconvolution algorithms estimate the RNA composition of homogenate tissue, e.g. the amount of RNA attributable to each cell type, and not the cellular composition, which relates to the underlying fraction of cells.