Alternative splicing (AS) is an important source of proteome diversity in eukaryotes. However, how this affects protein repertoires at a single-cell level remains an open question. Here, we show that many 3'-terminal exons are persistently co-expressed with their alternatives in mammalian neurons. In an important example of this scenario, cell polarity gene Cdc42, a combination of polypyrimidine tract-binding, protein-dependent, and constitutive splicing mechanisms ensures a halfway switch from the general (E7) to the neuron-specific (E6) alternative 3'-terminal exon during neuronal differentiation. Perturbing the nearly equimolar E6/E7 ratio in neurons results in defects in both axonal and dendritic compartments and suggests that Cdc42E7 is involved in axonogenesis, whereas Cdc42E6 is required for normal development of dendritic spines. Thus, co-expression of a precise blend of functionally distinct splice isoforms rather than a complete switch from one isoform to another underlies proper structural and functional polarization of neurons.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease resulting in severe muscle weakness and eventual death by respiratory failure. Although little is known about its pathogenesis, mutations in fused in sarcoma/translated in liposarcoma (FUS) are causative for familial ALS. FUS is a multifunctional protein that is involved in many aspects of RNA processing. To elucidate the role of FUS in ALS, we overexpressed wild-type and two mutant forms of FUS in HEK-293T cells, as well as knocked-down FUS expression. This was followed by RNA-Seq to identify genes which displayed differential expression or altered splicing patterns. Pathway analysis revealed that overexpression of wild-type FUS regulates ribosomal genes, whereas knock-down of FUS additionally affects expression of spliceosome related genes. Furthermore, cells expressing mutant FUS displayed global transcription patterns more similar to cells overexpressing wild-type FUS than to the knock-down condition. This observation suggests that FUS mutants do not contribute to the pathogenesis of ALS through a loss-of-function. Finally, our results demonstrate that the R521G and R522G mutations display differences in their influence on transcription and splicing. Taken together, these results provide additional insights into the function of FUS and how mutations contribute to the development of ALS.

Many subtypes of cortical interneurons (CINs) are found in adult mouse cortices, but the mechanism generating their diversity remains elusive. We performed single-cell RNA sequencing on the mouse embryonic medial ganglionic eminence (MGE), the major birthplace for CINs, and on MGE-like cells differentiated from embryonic stem cells. Two distinct cell types were identified as proliferating neural progenitors and immature neurons, both of which comprised sub-populations. Although lineage development of MGE progenitors was reconstructed and immature neurons were characterized as GABAergic, cells that might correspond to precursors of different CINs were not identified. A few non-neuronal cell types were detected, including microglia. In vitro MGE-like cells resembled bona fide MGE cells but expressed lower levels of Foxg1 and Epha4. Together, our data provide detailed understanding of the embryonic MGE developmental program and suggest how CINs are specified.

In multiple sclerosis (MS) and other neurological diseases, the failure to repair demyelinated lesions contributes to axonal damage and clinical disability. Here, we provide evidence that Mertk, a gene highly expressed by microglia that alters MS risk, is required for efficient remyelination. Compared to wild-type (WT) mice, Mertk-knockout (KO) mice show impaired clearance of myelin debris and remyelination following demyelination. Using single-cell RNA sequencing, we characterize Mertk-influenced responses to cuprizone-mediated demyelination and remyelination across different cell types. Mertk-KO brains show an attenuated microglial response to demyelination but an elevated proportion of interferon (IFN)-responsive microglia. In addition, we identify a transcriptionally distinct subtype of surviving oligodendrocytes specific to demyelinated lesions. The inhibitory effect of myelin debris on remyelination is mediated in part by IFNγ, which further impedes microglial clearance of myelin debris and inhibits oligodendrocyte differentiation. Together, our work establishes a role for Mertk in microglia activation, phagocytosis, and migration during remyelination.

Cortical circuit activity is shaped by the parvalbumin (PV) and somatostatin (SST) interneurons that inhibit principal excitatory (EXC) neurons and the vasoactive intestinal peptide (VIP) interneurons that suppress activation of other interneurons. To understand the molecular-genetic basis of functional specialization and identify potential drug targets specific to each neuron subtype, we performed a genome wide assessment of both gene expression and splicing across EXC, PV, SST and VIP neurons from male and female mouse brains. These results reveal numerous examples where neuron subtype-specific gene expression, as well as splice-isoform usage, can explain functional differences between neuron subtypes, including in presynaptic plasticity, postsynaptic receptor function, and synaptic connectivity specification. We provide a searchable web resource for exploring differential mRNA expression and splice form usage between excitatory, PV, SST, and VIP neurons (http://research-pub.gene.com/NeuronSubtypeTranscriptomes). This resource, combining a unique new dataset and novel application of analysis methods to multiple relevant datasets, identifies numerous potential drug targets for manipulating circuit function, reveals neuron subtype-specific roles for disease-linked genes, and is useful for understanding gene expression changes observed in human patient brains.SIGNIFICANCE STATEMENT Understanding the basis of functional specialization of neuron subtypes and identifying drug targets for manipulating circuit function requires comprehensive information on cell-type-specific transcriptional profiles. We sorted excitatory neurons and key inhibitory neuron subtypes from mouse brains and assessed differential mRNA expression. We used a genome-wide analysis which not only examined differential gene expression levels but could also detect differences in splice isoform usage. This analysis reveals numerous examples of neuron subtype-specific isoform usage with functional importance, identifies potential drug targets, and provides insight into the neuron subtypes involved in psychiatric disease. We also apply our analysis to two other relevant datasets for comparison, and provide a searchable website for convenient access to the resource.

Synapse loss and Tau pathology are hallmarks of Alzheimer's disease (AD) and other tauopathies, but how Tau pathology causes synapse loss is unclear. We used unbiased proteomic analysis of postsynaptic densities (PSDs) in Tau-P301S transgenic mice to identify Tau-dependent alterations in synapses prior to overt neurodegeneration. Multiple proteins and pathways were altered in Tau-P301S PSDs, including depletion of a set of GTPase-regulatory proteins that leads to actin cytoskeletal defects and loss of dendritic spines. Furthermore, we found striking accumulation of complement C1q in the PSDs of Tau-P301S mice and AD patients. At synapses, C1q decorated perisynaptic membranes, accumulated in correlation with phospho-Tau, and was associated with augmented microglial engulfment of synapses and decline of synapse density. A C1q-blocking antibody inhibited microglial synapse removal in cultured neurons and in Tau-P301S mice, rescuing synapse density. Thus, inhibiting complement-mediated synapse removal by microglia could be a potential therapeutic target for Tau-associated neurodegeneration.

RNA-Seq and microarray platforms have emerged as important tools for detecting changes in gene expression and RNA processing in biological samples. We present ExpressionPlot, a software package consisting of a default back end, which prepares raw sequencing or Affymetrix microarray data, and a web-based front end, which offers a biologically centered interface to browse, visualize, and compare different data sets. Download and installation instructions, a user's manual, discussion group, and a prototype are available at http://expressionplot.com/.

Damage-associated microglia (DAM) profiles observed in Alzheimer's disease (AD)-related mouse models reflect an activation state that could modulate AD risk or progression. To learn whether human AD microglia (HAM) display a similar profile, we develop a method for purifying cell types from frozen cerebrocortical tissues for RNA-seq analysis, allowing better transcriptome coverage than typical single-nucleus RNA-seq approaches. The HAM profile we observe bears little resemblance to the DAM profile. Instead, HAM display an enhanced human aging profile, in addition to other disease-related changes such as APOE upregulation. Analyses of whole-tissue RNA-seq and single-cell/nucleus RNA-seq datasets corroborate our findings and suggest that the lack of DAM response in human microglia occurs specifically in AD tissues, not other neurodegenerative settings. These results, which can be browsed at http://research-pub.gene.com/BrainMyeloidLandscape, provide a genome-wide picture of microglial activation in human AD and highlight considerable differences between mouse models and human disease.

Microglia and complement can mediate neurodegeneration in Alzheimer's disease (AD). By integrative multi-omics analysis, here we show that astrocytic and microglial proteins are increased in TauP301S synapse fractions with age and in a C1q-dependent manner. In addition to microglia, we identified that astrocytes contribute substantially to synapse elimination in TauP301S hippocampi. Notably, we found relatively more excitatory synapse marker proteins in astrocytic lysosomes, whereas microglial lysosomes contained more inhibitory synapse material. C1q deletion reduced astrocyte-synapse association and decreased astrocytic and microglial synapses engulfment in TauP301S mice and rescued synapse density. Finally, in an AD mouse model that combines β-amyloid and Tau pathologies, deletion of the AD risk gene Trem2 impaired microglial phagocytosis of synapses, whereas astrocytes engulfed more inhibitory synapses around plaques. Together, our data reveal that astrocytes contact and eliminate synapses in a C1q-dependent manner and thereby contribute to pathological synapse loss and that astrocytic phagocytosis can compensate for microglial dysfunction.

Heterozygous mutations in the granulin (GRN) gene are a leading cause of frontotemporal lobar degeneration with TDP-43 aggregates (FTLD-TDP). Polymorphisms in TMEM106B have been associated with disease risk in GRN mutation carriers and protective TMEM106B variants associated with reduced levels of TMEM106B, suggesting that lowering TMEM106B might be therapeutic in the context of FTLD. Here, we tested the impact of full deletion and partial reduction of TMEM106B in mouse and iPSC-derived human cell models of GRN deficiency. TMEM106B deletion did not reverse transcriptomic or proteomic profiles in GRN-deficient microglia, with a few exceptions in immune signaling markers. Neither homozygous nor heterozygous Tmem106b deletion normalized disease-associated phenotypes in Grn -/-mice. Furthermore, Tmem106b reduction by antisense oligonucleotide (ASO) was poorly tolerated in Grn -/-mice. These data provide novel insight into TMEM106B and GRN function in microglia cells but do not support lowering TMEM106B levels as a viable therapeutic strategy for treating FTD-GRN.

A common approach to understanding neurodegenerative disease is comparing gene expression in diseased versus healthy tissues. We illustrate that expression profiles derived from whole tissue RNA highly reflect the degenerating tissues' altered cellular composition, not necessarily transcriptional regulation. To accurately understand transcriptional changes that accompany neuropathology, we acutely purify neurons, astrocytes and microglia from single adult mouse brains and analyse their transcriptomes by RNA sequencing. Using peripheral endotoxemia to establish the method, we reveal highly specific transcriptional responses and altered RNA processing in each cell type, with Tnfr1 required for the astrocytic response. Extending the method to an Alzheimer's disease model, we confirm that transcriptomic changes observed in whole tissue are driven primarily by cell type composition, not transcriptional regulation, and identify hundreds of cell type-specific changes undetected in whole tissue RNA. Applying similar methods to additional models and patient tissues will transform our understanding of aberrant gene expression in neurological disease.

Loss-of-function mutations in GRN cause frontotemporal dementia (FTD) with transactive response DNA-binding protein of 43 kD (TDP-43)-positive inclusions and neuronal ceroid lipofuscinosis (NCL). There are no disease-modifying therapies for either FTD or NCL, in part because of a poor understanding of how mutations in genes such as GRN contribute to disease pathogenesis and neurodegeneration. By studying mice lacking progranulin (PGRN), the protein encoded by GRN, we discovered multiple lines of evidence that PGRN deficiency results in impairment of autophagy, a key cellular degradation pathway. PGRN-deficient mice are sensitive to Listeria monocytogenes because of deficits in xenophagy, a specialized form of autophagy that mediates clearance of intracellular pathogens. Cells lacking PGRN display reduced autophagic flux, and pathological forms of TDP-43 typically cleared by autophagy accumulate more rapidly in PGRN-deficient neurons. Our findings implicate autophagy as a novel therapeutic target for GRN-associated NCL and FTD and highlight the emerging theme of defective autophagy in the broader FTD/amyotrophic lateral sclerosis spectrum of neurodegenerative disease.

Antisense oligonucleotide (ASO) therapeutics are being investigated for a broad range of neurological diseases. While ASOs have been effective in the clinic, improving productive ASO internalization into target cells remains a key area of focus in the field. Here, we investigated how the delivery of ASO-loaded lipid nanoparticles (LNPs) affects ASO activity, subcellular trafficking, and distribution in the brain. We show that ASO-LNPs increase ASO activity up to 100-fold in cultured primary brain cells as compared to non-encapsulated ASO. However, in contrast to the widespread ASO uptake and activity observed following free ASO delivery in vivo, LNP-delivered ASOs did not downregulate mRNA levels throughout the brain after intracerebroventricular injection. This lack of activity was likely due to ASO accumulation in cells lining the ventricles and blood vessels. Furthermore, we reveal a formulation-dependent activation of the immune system post dosing, suggesting that LNP encapsulation cannot mask cellular ASO backbone-mediated toxicities. Together, these data provide insights into how LNP encapsulation affects ASO distribution as well as activity in the brain, and a foundation that enables future optimization of brain-targeting ASO-LNPs.

Loss-of-function TREM2 mutations strongly increase Alzheimer's disease (AD) risk. Trem2 deletion has revealed protective Trem2 functions in preclinical models of β-amyloidosis, a prominent feature of pre-diagnosis AD stages. How TREM2 influences later AD stages characterized by tau-mediated neurodegeneration is unclear. To understand Trem2 function in the context of both β-amyloid and tau pathologies, we examined Trem2 deficiency in the pR5-183 mouse model expressing mutant tau alone or in TauPS2APP mice, in which β-amyloid pathology exacerbates tau pathology and neurodegeneration. Single-cell RNA sequencing in these models revealed robust disease-associated microglia (DAM) activation in TauPS2APP mice that was amyloid-dependent and Trem2-dependent. In the presence of β-amyloid pathology, Trem2 deletion further exacerbated tau accumulation and spreading and promoted brain atrophy. Without β-amyloid pathology, Trem2 deletion did not affect these processes. Therefore, TREM2 may slow AD progression and reduce tau-driven neurodegeneration by restricting the degree to which β-amyloid facilitates the spreading of pathogenic tau.

A propensity for rewiring genetic and epigenetic regulatory networks, thus enabling sustained cell proliferation, suppression of apoptosis, and the ability to evade the immune system, is vital to cancer cell propagation. An increased understanding of how this is achieved is critical for identifying or improving therapeutic interventions. In this study, using acute myeloid leukemia (AML) human cell lines and a custom CRISPR/Cas9 screening platform, we identify the H3K9 methyltransferase SETDB1 as a novel, negative regulator of innate immunity. SETDB1 is overexpressed in many cancers, and loss of this gene in AML cells triggers desilencing of retrotransposable elements that leads to the production of double-stranded RNAs (dsRNAs). This is coincident with induction of a type I interferon response and apoptosis through the dsRNA-sensing pathway. Collectively, our findings establish a unique gene regulatory axis that cancer cells can exploit to circumvent the immune system.

The aggregation of intracellular tau protein is a major hallmark of Alzheimer's disease (AD). The extent and the stereotypical spread of tau pathology in the AD brain are correlated with cognitive decline during disease progression. Here we present an in-depth analysis of endogenous tau fragmentation in a well-characterized cohort of AD and age-matched control subjects. Using protein mass spectrometry and Edman degradation to interrogate endogenous tau fragments in the human brain, we identified two novel proteolytic sites, G323 and G326, as major tau cleavage events in both normal and AD cortex. These sites are located within the sequence recently identified as the structural core of tau protofilaments, suggesting an inhibitory mechanism of fibril formation. In contrast, a different set of novel cleavages showed a distinct increase in late stage AD. These disease-associated sites are located outside of the protofilament core sequence. We demonstrate that calpain 1 specifically cleaves at both the normal and diseased sites in vitro, and the site selection is conformation-dependent. Monomeric tau is predominantly cleaved at G323/G326 (normal sites), whereas oligomerization increases cleavages at the late-AD-associated sites. The fragmentation patterns specific to disease and healthy states suggest novel regulatory mechanisms of tau aggregation in the human brain.

TREM2 is an Alzheimer's disease (AD) risk gene expressed in microglia. To study the role of Trem2 in a mouse model of β-amyloidosis, we compared PS2APP transgenic mice versus PS2APP mice lacking Trem2 (PS2APP;Trem2ko) at ages ranging from 4 to 22 months. Microgliosis was impaired in PS2APP;Trem2ko mice, with Trem2-deficient microglia showing compromised expression of proliferation/Wnt-related genes and marked accumulation of ApoE. Plaque abundance was elevated in PS2APP;Trem2ko females at 6-7 months; but by 12 or 19-22 months of age, it was notably diminished in female and male PS2APP;Trem2ko mice, respectively. Across all ages, plaque morphology was more diffuse in PS2APP;Trem2ko brains, and the Aβ42:Aβ40 ratio was elevated. The amount of soluble, fibrillar Aβ oligomers also increased in PS2APP;Trem2ko hippocampi. Associated with these changes, axonal dystrophy was exacerbated from 6 to 7 months onward in PS2APP;Trem2ko mice, notwithstanding the reduced plaque load at later ages. PS2APP;Trem2ko mice also exhibited more dendritic spine loss around plaque and more neurofilament light chain in CSF. Thus, aggravated neuritic dystrophy is a more consistent outcome of Trem2 deficiency than amyloid plaque load, suggesting that the microglial packing of Aβ into dense plaque is an important neuroprotective activity.SIGNIFICANCE STATEMENT Genetic studies indicate that TREM2 gene mutations confer increased Alzheimer's disease (AD) risk. We studied the effects of Trem2 deletion in the PS2APP mouse AD model, in which overproduction of Aβ peptide leads to amyloid plaque formation and associated neuritic dystrophy. Interestingly, neuritic dystrophies were intensified in the brains of Trem2-deficient mice, despite these mice displaying reduced plaque accumulation at later ages (12-22 months). Microglial clustering around plaques was impaired, plaques were more diffuse, and the Aβ42:Aβ40 ratio and amount of soluble, fibrillar Aβ oligomers were elevated in Trem2-deficient brains. These results suggest that the Trem2-dependent compaction of Aβ into dense plaques is a protective microglial activity, limiting the exposure of neurons to toxic Aβ species.

Non-neuronal responses in neurodegenerative disease have received increasing attention as important contributors to disease pathogenesis and progression. Here we utilize single-cell RNA sequencing to broadly profile 13 cell types in three different mouse models of Alzheimer disease (AD), capturing the effects of tau-only, amyloid-only, or combined tau-amyloid pathology. We highlight microglia, oligodendrocyte, astrocyte, and T cell responses and compare them across these models. Notably, we identify two distinct transcriptional states for oligodendrocytes emerging differentially across disease models, and we determine their spatial distribution. Furthermore, we explore the impact of Trem2 deletion in the context of combined pathology. Trem2 knockout mice exhibit severely blunted microglial responses to combined tau and amyloid pathology, but responses from non-microglial cell types (oligodendrocytes, astrocytes, and T cells) are relatively unchanged. These results delineate core transcriptional states that are engaged in response to AD pathology, and how they are influenced by a key AD risk gene, Trem2.

The gene GPNMB is known to play roles in phagocytosis and tissue repair, and is upregulated in microglia in many mouse models of neurodegenerative disease as well as in human patients. Nearby genomic variants are associated with both elevated Parkinson's disease (PD) risk and higher expression of this gene, suggesting that inhibiting GPNMB activity might be protective in Parkinson's disease. We tested this hypothesis in three different mouse models of neurological diseases: a remyelination model and two models of alpha-synuclein pathology. We found that Gpnmb deletion had no effect on histological, cellular, behavioral, neurochemical or gene expression phenotypes in any of these models. These data suggest that Gpnmb does not play a major role in the development of pathology or functional defects in these models and that further work is necessary to study its role in the development or progression of Parkinson's disease.

Tumor progression locus 2 (TPL2) (MAP3K8) is a central signaling node in the inflammatory response of peripheral immune cells. We find that TPL2 kinase activity modulates microglial cytokine release and is required for microglia-mediated neuron death in vitro. In acute in vivo neuroinflammation settings, TPL2 kinase activity regulates microglia activation states and brain cytokine levels. In a tauopathy model of chronic neurodegeneration, loss of TPL2 kinase activity reduces neuroinflammation and rescues synapse loss, brain volume loss, and behavioral deficits. Single-cell RNA sequencing analysis indicates that protection in the tauopathy model was associated with reductions in activated microglia subpopulations as well as infiltrating peripheral immune cells. Overall, using various models, we find that TPL2 kinase activity can promote multiple harmful consequences of microglial activation in the brain including cytokine release, iNOS (inducible nitric oxide synthase) induction, astrocyte activation, and immune cell infiltration. Consequently, inhibiting TPL2 kinase activity could represent a potential therapeutic strategy in neurodegenerative conditions.