Genetic variants associated with Alzheimer's disease confer different cerebral cortex cell-type population structure

Genetic and multi-omic resources for Alzheimer disease and related dementia from the Knight Alzheimer Disease Research Center

The Knight-Alzheimer Disease Research Center (Knight-ADRC) at Washington University in St. Louis has pioneered and led worldwide seminal studies that have expanded our clinical, social, pathological, and molecular understanding of Alzheimer Disease. Over more than 40 years, research volunteers have been recruited to participate in cognitive, neuropsychologic, imaging, fluid biomarkers, genomic and multi-omic studies. Tissue and longitudinal data collected to foster, facilitate, and support research on dementia and aging. The Genetics and high throughput -omics core (GHTO) have collected of more than 26,000 biological samples from 6,625 Knight-ADRC participants. Samples available include longitudinal DNA, RNA, non-fasted plasma, cerebrospinal fluid pellets, and peripheral blood mononuclear cells. The GHTO has performed deep molecular profiling (genomic, transcriptomic, epigenomic, proteomic, and metabolomic) from large number of brain (n = 2,117), CSF (n = 2,012) and blood/plasma (n = 8,265) samples with the goal of identifying novel risk and protective variants, identify novel molecular biomarkers and causal and druggable targets. Overall, the resources available at GHTO support the increase of our understanding of Alzheimer Disease.

Brain Cell-based Genetic Subtyping and Drug Repositioning for Alzheimer Disease

Alzheimer's Disease (AD) is characterized by its complex and heterogeneous etiology and gradual progression, leading to high drug failure rates in late-stage clinical trials. In order to better stratify individuals at risk for AD and discern potential therapeutic targets we employed a novel procedure utilizing cell-based co-regulated gene networks and polygenic risk scores (cbPRSs). After defining genetic subtypes using extremes of cbPRS distributions, we evaluated correlations of the genetic subtypes with previously defined AD subtypes defined on the basis of domain-specific cognitive functioning and neuroimaging biomarkers. Employing a PageRank algorithm, we identified priority gene targets for the genetic subtypes. Pathway analysis of priority genes demonstrated associations with neurodegeneration and suggested candidate drugs currently utilized in diabetes, hypertension, and epilepsy for repositioning in AD. Experimental validation utilizing human induced pluripotent stem cell (hiPSC)-derived astrocytes demonstrated the modifying effects of estradiol, levetiracetam, and pioglitazone on expression of APOE and complement C4 genes, suggesting potential repositioning for AD.

Brain cell-type shifts in Alzheimer's disease, autism, and schizophrenia interrogated using methylomics and genetics

Few neuropsychiatric disorders have replicable biomarkers, prompting high-resolution and large-scale molecular studies. However, we still lack consensus on a more foundational question: whether quantitative shifts in cell types-the functional unit of life-contribute to neuropsychiatric disorders. Leveraging advances in human brain single-cell methylomics, we deconvolve seven major cell types using bulk DNA methylation profiling across 1270 postmortem brains, including from individuals diagnosed with Alzheimer's disease, schizophrenia, and autism. We observe and replicate cell-type compositional shifts for Alzheimer's disease (endothelial cell loss), autism (increased microglia), and schizophrenia (decreased oligodendrocytes), and find age- and sex-related changes. Multiple layers of evidence indicate that endothelial cell loss contributes to Alzheimer's disease, with comparable effect size to APOE genotype among older people. Genome-wide association identified five genetic loci related to cell-type composition, involving plausible genes for the neurovascular unit (P2RX5 and TRPV3) and excitatory neurons (DPY30 and MEMO1). These results implicate specific cell-type shifts in the pathophysiology of neuropsychiatric disorders.

Alzheimer's disease rewires gene coexpression networks coupling different brain regions

Connectome studies have shown how Alzheimer's disease (AD) disrupts functional and structural connectivity among brain regions. But the molecular basis of such disruptions is less studied, with most genomic/transcriptomic studies performing within-brain-region analyses. To inspect how AD rewires the correlation structure among genes in different brain regions, we performed an Inter-brain-region Differential Correlation (Inter-DC) analysis of RNA-seq data from Mount Sinai Brain Bank on four brain regions (frontal pole, superior temporal gyrus, parahippocampal gyrus and inferior frontal gyrus, comprising 264 AD and 372 control human post-mortem samples). An Inter-DC network was assembled from all pairs of genes across two brain regions that gained (or lost) correlation strength in the AD group relative to controls at FDR 1%. The differentially correlated (DC) genes in this network complemented known differentially expressed genes in AD, and likely reflects cell-intrinsic changes since we adjusted for cell compositional effects. Each brain region used a distinctive set of DC genes when coupling with other regions, with parahippocampal gyrus showing the most rewiring, consistent with its known vulnerability to AD. The Inter-DC network revealed master dysregulation hubs in AD (at genes ZKSCAN1, SLC5A3, RCC1, IL17RB, PLK4, etc.), inter-region gene modules enriched for known AD pathways (synaptic signaling, endocytosis, etc.), and candidate signaling molecules that could mediate region-region communication. The Inter-DC network generated in this study is a valuable resource of gene pairs, pathways and signaling molecules whose inter-brain-region functional coupling is disrupted in AD, thereby offering a new perspective of AD etiology.

Circulating blood circular RNA in Parkinson's Disease; a systematic study

We aimed to identify circRNAs associated with Parkinson's disease (PD) by leveraging 1,848 participants and 1,789 circRNA from two of the largest publicly available studies with longitudinal clinical and blood transcriptomic data. To comprehensively understand changes in circRNAs we performed a cross-sectional study utilizing the last visit of each participant, and a longitudinal (mix model) analysis that included 1,166 participants with at least two time points. We identified 192 circRNAs differentially expressed in PD participants compared to healthy controls, with effects that were consistent in the mixed models, mutation carriers, and diverse ancestry. Finally, we included the 149 circRNA in a model with a ROC AUC of 0.825, showing that have the potential to aid the diagnosis of PD. Overall, we demonstrated that circRNAs play an important role in PD and can be leveraged as biomarkers.

Computational saturation mutagenesis to explore the effect of pathogenic mutations on extra-cellular domains of TREM2 associated with Alzheimer's and Nasu-Hakola disease

CONTEXT

The specialised family of triggering receptors expressed on myeloid cells (TREMs) plays a pivotal role in causing neurodegenerative disorders and activating microglial anti-inflammatory responses. Nasu-Hakola disease (NHD), a rare autosomal recessive disorder, has been associated with mutations in TREM2, which is also responsible for raising the risk of Alzheimer's disease (AD). Herein, we have made an endeavour to differentiate the confirmed pathogenic variants in TREM2 extra-cellular domain (ECD) linked with NHD and AD using mutation-induced fold stability change (∆∆G), with the computation of 12distinct structure-based methods through saturation mutagenesis. Correlation analysis between relative solvent accessibility (RSA) and ∆∆G expresses the discrete distributive behaviour of mutants associated with TREM2 in AD (R2 = 0.061) and NHD (R2 = 0.601). Our findings put an emphasis on W50 and V126 as major players in maintaining V-like domain in TREM2. Interestingly, we discern that both of them interact with a common residue Y108, which is dissolved upon mutation. This Y108 could have structural or functional role for TREM2 which can be an ideal candidate for further study. Furthermore, the residual interaction network highlights the importance of R47 and R62 in maintaining the CDR loops that are crucial for ligand binding. Future studies using biophysical characterisation of ligand interactions in TREM2-ECD would be helpful for the development of novel therapeutics for AD and NHD.

METHODS

ConSurf algorithm and ENDscript were used to determine the position and conservation of each residue in the wild-type ECD of TREM2. The mutation-induced fold stability change (∆∆G) of confirmed pathogenic mutants associated with NHD and AD was estimated using 12 state-of-the-art structure-based protein stability tools. Furthermore, we also computed the effect of random mutation on these sites using computational saturation mutagenesis. Linear regression analysis was performed using mutants ∆∆G and RSA through GraphPad software. In addition, a comprehensive non-bonded residual interaction network (RIN) of wild type and its mutants of TREM2-ECD was enumerated using RING3.0.

Brain expression quantitative trait locus and network analyses reveal downstream effects and putative drivers for brain-related diseases

Identification of therapeutic targets from genome-wide association studies (GWAS) requires insights into downstream functional consequences. We harmonized 8,613 RNA-sequencing samples from 14 brain datasets to create the MetaBrain resource and performed cis- and trans-expression quantitative trait locus (eQTL) meta-analyses in multiple brain region- and ancestry-specific datasets (n ≤ 2,759). Many of the 16,169 cortex cis-eQTLs were tissue-dependent when compared with blood cis-eQTLs. We inferred brain cell types for 3,549 cis-eQTLs by interaction analysis. We prioritized 186 cis-eQTLs for 31 brain-related traits using Mendelian randomization and co-localization including 40 cis-eQTLs with an inferred cell type, such as a neuron-specific cis-eQTL (CYP24A1) for multiple sclerosis. We further describe 737 trans-eQTLs for 526 unique variants and 108 unique genes. We used brain-specific gene-co-regulation networks to link GWAS loci and prioritize additional genes for five central nervous system diseases. This study represents a valuable resource for post-GWAS research on central nervous system diseases.

Conserved gene signatures shared among MAPT mutations reveal defects in calcium signaling

Introduction: More than 50 mutations in the MAPT gene result in heterogeneous forms of frontotemporal lobar dementia with tau inclusions (FTLD-Tau). However, early pathogenic events that lead to disease and the degree to which they are common across MAPT mutations remain poorly understood. The goal of this study is to determine whether there is a common molecular signature of FTLD-Tau. Methods: We analyzed genes differentially expressed in induced pluripotent stem cell-derived neurons (iPSC-neurons) that represent the three major categories of MAPT mutations: splicing (IVS10 + 16), exon 10 (p.P301L), and C-terminal (p.R406W) compared with isogenic controls. The genes that were commonly differentially expressed in MAPT IVS10 + 16, p.P301L, and p.R406W neurons were enriched in trans-synaptic signaling, neuronal processes, and lysosomal function. Many of these pathways are sensitive to disruptions in calcium homeostasis. One gene, CALB1, was significantly reduced across the three MAPT mutant iPSC-neurons and in a mouse model of tau accumulation. We observed a significant reduction in calcium levels in MAPT mutant neurons compared with isogenic controls, pointing to a functional consequence of this disrupted gene expression. Finally, a subset of genes commonly differentially expressed across MAPT mutations were also dysregulated in brains from MAPT mutation carriers and to a lesser extent in brains from sporadic Alzheimer disease and progressive supranuclear palsy, suggesting that molecular signatures relevant to genetic and sporadic forms of tauopathy are captured in a dish. The results from this study demonstrate that iPSC-neurons capture molecular processes that occur in human brains and can be used to pinpoint common molecular pathways involving synaptic and lysosomal function and neuronal development, which may be regulated by disruptions in calcium homeostasis.

MS4A4A modifies the risk of Alzheimer disease by regulating lipid metabolism and immune response in a unique microglia state

Genome-wide association studies (GWAS) have identified many modifiers of Alzheimer disease (AD) risk enriched in microglia. Two of these modifiers are common variants in the MS4A locus (rs1582763: protective and rs6591561: risk) and serve as major regulators of CSF sTREM2 levels. To understand their functional impact on AD, we used single nucleus transcriptomics to profile brains from carriers of these variants. We discovered a "chemokine" microglial subpopulation that is altered in MS4A variant carriers and for which MS4A4A is the major regulator. The protective variant increases MS4A4A expression and shifts the chemokine microglia subpopulation to an interferon state, while the risk variant suppresses MS4A4A expression and reduces this subpopulation of microglia. Our findings provide a mechanistic explanation for the AD variants in the MS4A locus. Further, they pave the way for future mechanistic studies of AD variants and potential therapeutic strategies for enhancing microglia resilience in AD pathogenesis.

The Cholesteryl Ester Transfer Protein (CETP) raises Cholesterol Levels in the Brain

The cholesteryl ester transfer protein (CETP) is a lipid transfer protein responsible for the exchange of cholesteryl esters and triglycerides between lipoproteins. Decreased CETP activity is associated with longevity, cardiovascular health, and maintenance of good cognitive performance. Interestingly, mice lack the CETP-encoding gene and have very low levels of LDL particles compared with humans. Currently, the molecular mechanisms induced because of CETP activity are not clear. To understand how CETP activity affects the brain, we utilized CETP transgenic (CETPtg) mice that show elevated LDL levels upon induction of CETP expression through a high-cholesterol diet. CETPtg mice on a high-cholesterol diet showed up to 22% higher cholesterol levels in the brain. Using a microarray on mostly astrocyte-derived mRNA, we found that this cholesterol increase is likely not because of elevated de novo synthesis of cholesterol. However, cholesterol efflux is decreased in CETPtg mice along with an upregulation of the complement factor C1Q, which plays a role in neuronal cholesterol clearance. Our data suggest that CETP activity affects brain health through modulating cholesterol distribution and clearance. Therefore, we propose that CETPtg mice constitute a valuable research tool to investigate the impact of cholesterol metabolism on brain function.

Bulk and Single-Nucleus Transcriptomics Highlight Intra-Telencephalic and Somatostatin Neurons in Alzheimer's Disease

Cortical neuron loss is a pathological hallmark of late-onset Alzheimer's disease (AD). However, it remains unclear which neuronal subtypes beyond broad excitatory and inhibitory classes are most vulnerable. Here, we analyzed cell subtype proportion differences in AD compared to non-AD controls using 1037 post-mortem brain samples from six neocortical regions. We identified the strongest associations of AD with fewer somatostatin (SST) inhibitory neurons (β = -0.48, p bonf = 8.98 × 10-9) and intra-telencephalic (IT) excitatory neurons (β = -0.45, p bonf = 4.32 × 10-7). Replication in three AD case-control single-nucleus RNAseq datasets most strongly supported the bulk tissue association of fewer SST neurons in AD. In depth analyses of cell type proportions with specific AD-related neuropathological and cognitive phenotypes revealed fewer SST neurons with greater brain-wide post-mortem tau and beta amyloid, as well as a faster rate of antemortem cognitive decline. In contrast, greater IT neuron proportions were associated with a slower rate of cognitive decline as well as greater residual cognition-a measure of cognitive resilience-but not canonical AD neuropathology. Our findings implicate somatostatin inhibitory and intra-telencephalic excitatory neuron subclasses in the pathogenesis of AD and in cognitive resilience to AD pathology, respectively.

Comprehensive evaluation of deconvolution methods for human brain gene expression

Transcriptome deconvolution aims to estimate the cellular composition of an RNA sample from its gene expression data, which in turn can be used to correct for composition differences across samples. The human brain is unique in its transcriptomic diversity, and comprises a complex mixture of cell-types, including transcriptionally similar subtypes of neurons. Here, we carry out a comprehensive evaluation of deconvolution methods for human brain transcriptome data, and assess the tissue-specificity of our key observations by comparison with human pancreas and heart. We evaluate eight transcriptome deconvolution approaches and nine cell-type signatures, testing the accuracy of deconvolution using in silico mixtures of single-cell RNA-seq data, RNA mixtures, as well as nearly 2000 human brain samples. Our results identify the main factors that drive deconvolution accuracy for brain data, and highlight the importance of biological factors influencing cell-type signatures, such as brain region and in vitro cell culturing.

Transcriptome deconvolution aims to estimate cellular composition based on gene expression data. Here the authors evaluate deconvolution methods for human brain transcriptome and conclude that partial deconvolution algorithms work best, but that appropriate cell-type signatures are also important.

Circular RNA detection identifies circPSEN1 alterations in brain specific to autosomal dominant Alzheimer's disease

BACKGROUND

Autosomal-dominant Alzheimer's disease (ADAD) is caused by pathogenic mutations in APP, PSEN1, and PSEN2, which usually lead to an early age at onset (< 65). Circular RNAs are a family of non-coding RNAs highly expressed in the nervous system and especially in synapses. We aimed to investigate differences in brain gene expression of linear and circular transcripts from the three ADAD genes in controls, sporadic AD, and ADAD.

METHODS

We obtained and sequenced RNA from brain cortex using standard protocols. Linear counts were obtained using the TOPMed pipeline; circular counts, using python package DCC. After stringent quality control (QC), we obtained the counts for PSEN1, PSEN2 and APP genes. Only circPSEN1 passed QC. We used DESeq2 to compare the counts across groups, correcting for biological and technical variables. Finally, we performed in-silico functional analyses using the Circular RNA interactome website and DIANA mirPath software.

RESULTS

Our results show significant differences in gene counts of circPSEN1 in ADAD individuals, when compared to sporadic AD and controls (ADAD = 21, AD = 253, Controls = 23-ADADvsCO: log2FC = 0.794, p = 1.63 × 10-04, ADADvsAD: log2FC = 0.602, p = 8.22 × 10-04). The high gene counts are contributed by two circPSEN1 species (hsa_circ_0008521 and hsa_circ_0003848). No significant differences were observed in linear PSEN1 gene expression between cases and controls, indicating that this finding is specific to the circular forms. In addition, the high circPSEN1 levels do not seem to be specific to PSEN1 mutation carriers; the counts are also elevated in APP and PSEN2 mutation carriers. In-silico functional analyses suggest that circPSEN1 is involved in several pathways such as axon guidance (p = 3.39 × 10-07), hippo signaling pathway (p = 7.38 × 10-07), lysine degradation (p = 2.48 × 10-05) or Wnt signaling pathway (p = 5.58 × 10-04) among other KEGG pathways. Additionally, circPSEN1 counts were able to discriminate ADAD from sporadic AD and controls with an AUC above 0.70.

CONCLUSIONS

Our findings show the differential expression of circPSEN1 is increased in ADAD. Given the biological function previously ascribed to circular RNAs and the results of our in-silico analyses, we hypothesize that this finding might be related to neuroinflammatory events that lead or that are caused by the accumulation of amyloid-beta.

Multi-ancestry GWAS reveals excitotoxicity associated with outcome after ischaemic stroke

During the first hours after stroke onset, neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between the National Institutes of Health Stroke Scale (NIHSS) within 6 h of stroke onset and NIHSS at 24 h. A total of 5876 individuals from seven countries (Spain, Finland, Poland, USA, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of NIHSS at 24 h of variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture from that of stroke risk. Eight loci (1p21.1, 1q42.2, 2p25.1, 2q31.2, 2q33.3, 5q33.2, 7p21.2 and 13q31.1) were genome-wide significant and explained 1.8% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each locus. Expression quantitative trait loci mapping and summary data-based Mendelian randomization indicate that ADAM23 (log Bayes factor = 5.41) was driving the association for 2q33.3. Gene-based analyses suggested that GRIA1 (log Bayes factor = 5.19), which is predominantly expressed in the brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated GNPAT (log Bayes factor = 7.64) ABCB5 (log Bayes factor = 5.97) for the 1p21.1 and 7p21.1 loci. Human brain single-nuclei RNA-sequencing indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23, a presynaptic protein and GRIA1, a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provide the first genetic evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischaemic stroke.

Murine roseolovirus does not accelerate amyloid-β pathology and human roseoloviruses are not over-represented in Alzheimer disease brains

BACKGROUND

The role of viral infection in Alzheimer Disease (AD) pathogenesis is an area of great interest in recent years. Several studies have suggested an association between the human roseoloviruses, HHV-6 and HHV-7, and AD. Amyloid-β (Aβ) plaques are a hallmark neuropathological finding of AD and were recently proposed to have an antimicrobial function in response to infection. Identifying a causative and mechanistic role of human roseoloviruses in AD has been confounded by limitations in performing in vivo studies. Recent -omics based approaches have demonstrated conflicting associations between human roseoloviruses and AD. Murine roseolovirus (MRV) is a natural murine pathogen that is highly-related to the human roseoloviruses, providing an opportunity to perform well-controlled studies of the impact of roseolovirus on Aβ deposition.

METHODS

We utilized the 5XFAD mouse model to test whether MRV induces Aβ deposition in vivo. We also evaluated viral load and neuropathogenesis of MRV infection. To evaluate Aβ interaction with MRV, we performed electron microscopy. RNA-sequencing of a cohort of AD brains compared to control was used to investigate the association between human roseolovirus and AD.

RESULTS

We found that 5XFAD mice were susceptible to MRV infection and developed neuroinflammation. Moreover, we demonstrated that Aβ interacts with viral particles in vitro and, subsequent to this interaction, can disrupt infection. Despite this, neither peripheral nor brain infection with MRV increased or accelerated Aβ plaque formation. Moreover, -omics based approaches have demonstrated conflicting associations between human roseoloviruses and AD. Our RNA-sequencing analysis of a cohort of AD brains compared to controls did not show an association between roseolovirus infection and AD.

CONCLUSION

Although MRV does infect the brain and cause transient neuroinflammation, our data do not support a role for murine or human roseoloviruses in the development of Aβ plaque formation and AD.

A diet rich in docosahexaenoic acid enhances reactive astrogliosis and ramified microglia morphology in apolipoprotein E epsilon 4-targeted replacement mice

Docosahexaenoic acid (DHA) consumption reduces spatial memory impairment in mice carrying the human apolipoprotein E ε4 (APOE4) allele. The current study evaluated whether astrocyte and microglia morphology contribute to the mechanism of this result. APOE3 and APOE4 mice were fed either a DHA-enriched diet or a control diet from 4 to 12 months of age. Coronal brain sections were immunostained for GFAP, Iba1, and NeuN. Astrocytes from APOE4 mice exhibited signs of reactive astrogliosis compared to APOE3 mice. Consumption of DHA exacerbated reactive astrocyte morphology in APOE4 carriers. Microglia from APOE4-control mice exhibited characteristics of amoeboid morphology and other characteristics of ramified morphology (more processes, greater process complexity, and greater distance between neighboring microglia). DHA enhanced ramified microglia morphology in APOE4 mice. In addition, APOE4 mice fed the DHA diet had lower hippocampal concentrations of interleukin-7, lipopolysaccharide-induced CXC chemokine and monocyte chemoattractant protein 1, and higher concentration of interferon-gamma compared to APOE4-control mice. Our results indicate that a diet rich in DHA enhances reactive astrogliosis and ramified microglia morphology in APOE4 mice.

Alzheimer's disease alters oligodendrocytic glycolytic and ketolytic gene expression

INTRODUCTION

Sporadic Alzheimer's disease (AD) is strongly correlated with impaired brain glucose metabolism, which may affect AD onset and progression. Ketolysis has been suggested as an alternative pathway to fuel the brain.

METHODS

RNA-seq profiles of post mortem AD brains were used to determine whether dysfunctional AD brain metabolism can be determined by impairments in glycolytic and ketolytic gene expression. Data were obtained from the Knight Alzheimer's Disease Research Center (62 cases; 13 controls), Mount Sinai Brain Bank (110 cases; 44 controls), and the Mayo Clinic Brain Bank (80 cases; 76 controls), and were normalized to cell type: astrocytes, microglia, neurons, oligodendrocytes.

RESULTS

In oligodendrocytes, both glycolytic and ketolytic pathways were significantly impaired in AD brains. Ketolytic gene expression was not significantly altered in neurons, astrocytes, and microglia.

DISCUSSION

Oligodendrocytes may contribute to brain hypometabolism observed in AD. These results are suggestive of a potential link between hypometabolism and dysmyelination in disease physiology. Additionally, ketones may be therapeutic in AD due to their ability to fuel neurons despite impaired glycolytic metabolism.

Advances in Genetic and Molecular Understanding of Alzheimer's Disease

Alzheimer's disease (AD) has become a common disease of the elderly for which no cure currently exists. After over 30 years of intensive research, we have gained extensive knowledge of the genetic and molecular factors involved and their interplay in disease. These findings suggest that different subgroups of AD may exist. Not only are we starting to treat autosomal dominant cases differently from sporadic cases, but we could be observing different underlying pathological mechanisms related to the amyloid cascade hypothesis, immune dysfunction, and a tau-dependent pathology. Genetic, molecular, and, more recently, multi-omic evidence support each of these scenarios, which are highly interconnected but can also point to the different subgroups of AD. The identification of the pathologic triggers and order of events in the disease processes are key to the design of treatments and therapies. Prevention and treatment of AD cannot be attempted using a single approach; different therapeutic strategies at specific disease stages may be appropriate. For successful prevention and treatment, biomarker assays must be designed so that patients can be more accurately monitored at specific points during the course of the disease and potential treatment. In addition, to advance the development of therapeutic drugs, models that better mimic the complexity of the human brain are needed; there have been several advances in this arena. Here, we review significant, recent developments in genetics, omics, and molecular studies that have contributed to the understanding of this disease. We also discuss the implications that these contributions have on medicine.

Emerging genetic complexity and rare genetic variants in neurodegenerative brain diseases

Knowledge of the molecular etiology of neurodegenerative brain diseases (NBD) has substantially increased over the past three decades. Early genetic studies of NBD families identified rare and highly penetrant deleterious mutations in causal genes that segregate with disease. Large genome-wide association studies uncovered common genetic variants that influenced disease risk. Major developments in next-generation sequencing (NGS) technologies accelerated gene discoveries at an unprecedented rate and revealed novel pathways underlying NBD pathogenesis. NGS technology exposed large numbers of rare genetic variants of uncertain significance (VUS) in coding regions, highlighting the genetic complexity of NBD. Since experimental studies of these coding rare VUS are largely lacking, the potential contributions of VUS to NBD etiology remain unknown. In this review, we summarize novel findings in NBD genetic etiology driven by NGS and the impact of rare VUS on NBD etiology. We consider different mechanisms by which rare VUS can act and influence NBD pathophysiology and discuss why a better understanding of rare VUS is instrumental for deriving novel insights into the molecular complexity and heterogeneity of NBD. New knowledge might open avenues for effective personalized therapies.

Comparison of machine learning approaches for enhancing Alzheimer's disease classification

Alzheimer's disease (AD) is a progressive neurodegenerative disorder, accounting for nearly 60% of all dementia cases. The occurrence of the disease has been increasing rapidly in recent years. Presently about 46.8 million individuals suffer from AD worldwide. The current absence of effective treatment to reverse or stop AD progression highlights the importance of disease prevention and early diagnosis. Brain structural Magnetic Resonance Imaging (MRI) has been widely used for AD detection as it can display morphometric differences and cerebral structural changes. In this study, we built three machine learning-based MRI data classifiers to predict AD and infer the brain regions that contribute to disease development and progression. We then systematically compared the three distinct classifiers, which were constructed based on Support Vector Machine (SVM), 3D Very Deep Convolutional Network (VGGNet) and 3D Deep Residual Network (ResNet), respectively. To improve the performance of the deep learning classifiers, we applied a transfer learning strategy. The weights of a pre-trained model were transferred and adopted as the initial weights of our models. Transferring the learned features significantly reduced training time and increased network efficiency. The classification accuracy for AD subjects from elderly control subjects was 90%, 95%, and 95% for the SVM, VGGNet and ResNet classifiers, respectively. Gradient-weighted Class Activation Mapping (Grad-CAM) was employed to show discriminative regions that contributed most to the AD classification by utilizing the learned spatial information of the 3D-VGGNet and 3D-ResNet models. The resulted maps consistently highlighted several disease-associated brain regions, particularly the cerebellum which is a relatively neglected brain region in the present AD study. Overall, our comparisons suggested that the ResNet model provided the best classification performance as well as more accurate localization of disease-associated regions in the brain compared to the other two approaches.

Quantitative endophenotypes as an alternative approach to understanding genetic risk in neurodegenerative diseases

Endophenotypes, as measurable intermediate features of human diseases, reflect underlying molecular mechanisms. The use of quantitative endophenotypes in genetic studies has improved our understanding of pathophysiological changes associated with diseases. The main advantage of the quantitative endophenotypes approach to study human diseases over a classic case-control study design is the inferred biological context that can enable the development of effective disease-modifying treatments. Here, we summarize recent progress on biomarkers for neurodegenerative diseases, including cerebrospinal fluid and blood-based, neuroimaging, neuropathological, and clinical studies. This review focuses on how endophenotypic studies have successfully linked genetic modifiers to disease risk, disease onset, or progression rate and provided biological context to genes identified in genome-wide association studies. Finally, we review critical methodological considerations for implementing this approach and future directions.

Commonalities between Copper Neurotoxicity and Alzheimer's Disease

Alzheimer's disease, a highly prevalent form of dementia, targets neuron function beginning from the hippocampal region and expanding outwards. Alzheimer's disease is caused by elevated levels of heavy metals, such as lead, zinc, and copper. Copper is found in many areas of daily life, raising a concern as to how this metal and Alzheimer's disease are related. Previous studies have not identified the common pathways between excess copper and Alzheimer's disease etiology. Our review corroborates that both copper and Alzheimer's disease target the hippocampus, cerebral cortex, cerebellum, and brainstem, affecting motor skills and critical thinking. Additionally, Aβ plaque formation was analyzed beginning from synthesis at the APP parent protein site until Aβ plaque formation was completed. Structural changes were also noted. Further analysis revealed a relationship between amyloid-beta plaques and copper ion concentration. As copper ion levels increased, it bound to the Aβ monomer, expediting the plaque formation process, and furthering neurodegeneration. These conclusions can be utilized in the medical community to further research on the etiology of Alzheimer's disease and its relationships to copper and other metal-induced neurotoxicity.

Multi-ancestry genetic study in 5,876 patients identifies an association between excitotoxic genes and early outcomes after acute ischemic stroke

During the first hours after stroke onset neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between NIH stroke scale (NIHSS) within six hours of stroke onset and NIHSS at 24h (ΔNIHSS). A total of 5,876 individuals from seven countries (Spain, Finland, Poland, United States, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of ΔNIHSS variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture than that of stroke risk. Seven loci (2p25.1, 2q31.2, 2q33.3, 4q34.3, 5q33.2, 6q26 and 7p21.1) were genome-wide significant and explained 2.1% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each loci. eQTL mapping and SMR indicate that ADAM23 (log Bayes Factor (LBF)=6.34) was driving the association for 2q33.3. Gene based analyses suggested that GRIA1 (LBF=5.26), which is predominantly expressed in brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated PARK2 (LBF=5.30) and ABCB5 (LBF=5.70) for the 6q26 and 7p21.1 loci. Human brain single nuclei RNA-seq indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23 , a pre-synaptic protein, and GRIA1 , a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provides the first evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischemic stroke.

RESEARCH INTO CONTEXT

Evidence before this study: No previous genome-wide association studies have investigated the genetic architecture of early outcomes after ischemic stroke.Added Value of this study: This is the first study that investigated genetic influences on early outcomes after ischemic stroke using a genome-wide approach, revealing seven genome-wide significant loci. A unique aspect of this genetic study is the inclusion of all of the major ethnicities by recruiting from participants throughout the world. Most genetic studies to date have been limited to populations of European ancestry.Implications of all available evidence: The findings provide the first evidence that genes implicating excitotoxicity contribute to human acute ischemic stroke, and demonstrates proof of principle that GWAS of acute ischemic stroke patients can reveal mechanisms involved in ischemic brain injury.

Deciphering cellular transcriptional alterations in Alzheimer's disease brains

Large-scale brain bulk-RNAseq studies identified molecular pathways implicated in Alzheimer's disease (AD), however these findings can be confounded by cellular composition changes in bulk-tissue. To identify cell intrinsic gene expression alterations of individual cell types, we designed a bioinformatics pipeline and analyzed three AD and control bulk-RNAseq datasets of temporal and dorsolateral prefrontal cortex from 685 brain samples. We detected cell-proportion changes in AD brains that are robustly replicable across the three independently assessed cohorts. We applied three different algorithms including our in-house algorithm to identify cell intrinsic differentially expressed genes in individual cell types (CI-DEGs). We assessed the performance of all algorithms by comparison to single nucleus RNAseq data. We identified consensus CI-DEGs that are common to multiple brain regions. Despite significant overlap between consensus CI-DEGs and bulk-DEGs, many CI-DEGs were absent from bulk-DEGs. Consensus CI-DEGs and their enriched GO terms include genes and pathways previously implicated in AD or neurodegeneration, as well as novel ones. We demonstrated that the detection of CI-DEGs through computational deconvolution methods is promising and highlight remaining challenges. These findings provide novel insights into cell-intrinsic transcriptional changes of individual cell types in AD and may refine discovery and modeling of molecular targets that drive this complex disease.

Silencing of aquaporin1 activates the Wnt signaling pathway to improve cognitive function in a mouse model of Alzheimer's disease

BACKGROUND

Aquaporins (AQPs) are water channel proteins robustly presenting in the central nervous system (CNS). Increasing evidence suggests the crucial role of AQP1 in the pathogenesis of CNS injury but scarce data are provided for the potential role of AQP1 in Alzheimer's disease (AD). Thus, the present study aimed to investigate the effects of AQP1 on cognitive function in a mouse model of AD.

METHODS

A mouse model of AD was established by using the β-amyloid isoform Aβ_1-42, and then assessed by the step-through test and water maze experiment. The expression of AQP1 was quantified in the AD model. The effects of AQP1 on the cognitive function of AD mice and the Wnt signaling pathway were elucidated using gain- and loss-of-function approaches. Furthermore, hippocampal neurons were isolated and treated with Aβ_1-42 for in vitro experiments and the effects of the Wnt signaling pathway on hippocampal neuron apoptosis were analyzed with the use of inhibitor or activator of this pathway.

RESULTS

AQP1 was highly-expressed in the AD mouse model while AQP1 silencing improved cognitive function in AD mice. Besides, silencing of AQP1 exhibited protective effects on hippocampal neurons in AD mice. Furthermore, AQP1 inhibited the Wnt signaling pathway while AQP1 promoted neuronal apoptosis by inhibiting the Wnt signaling pathway, thereby damaging the cognitive function.

CONCLUSIONS

AQP1 silencing attenuates the cognitive impairment in AD through activation of the Wnt signaling pathway, highlighting a novel therapeutic target against AD.

Overlapping genetic architecture between Parkinson disease and melanoma

Epidemiologic studies have reported inconsistent results regarding an association between Parkinson disease (PD) and cutaneous melanoma (melanoma). Identifying shared genetic architecture between these diseases can support epidemiologic findings and identify common risk genes and biological pathways. Here, we apply polygenic, linkage disequilibrium-informed methods to the largest available case-control, genome-wide association study summary statistic data for melanoma and PD. We identify positive and significant genetic correlation (correlation: 0.17, 95% CI 0.10-0.24; P = 4.09 × 10-06) between melanoma and PD. We further demonstrate melanoma and PD-inferred gene expression to overlap across tissues (correlation: 0.14, 95% CI 0.06 to 0.22; P = 7.87 × 10-04) and highlight seven genes including PIEZO1, TRAPPC2L, and SOX6 as potential mediators of the genetic correlation between melanoma and PD. These findings demonstrate specific, shared genetic architecture between PD and melanoma that manifests at the level of gene expression.

The TMEM106B FTLD-protective variant, rs1990621, is also associated with increased neuronal proportion

Apart from amyloid β deposition and tau neurofibrillary tangles, Alzheimer's disease (AD) is a neurodegenerative disorder characterized by neuronal loss and astrocytosis in the cerebral cortex. The goal of this study is to investigate genetic factors associated with the neuronal proportion in health and disease. To identify cell-autonomous genetic variants associated with neuronal proportion in cortical tissues, we inferred cellular population structure from bulk RNA-Seq derived from 1536 individuals. We identified the variant rs1990621 located in the TMEM106B gene region as significantly associated with neuronal proportion (p value = 6.40 × 10-07) and replicated this finding in an independent dataset (p value = 7.41 × 10-04) surpassing the genome-wide threshold in the meta-analysis (p value = 9.42 × 10-09). This variant is in high LD with the TMEM106B non-synonymous variant p.T185S (rs3173615; r2 = 0.98) which was previously identified as a protective variant for frontotemporal lobar degeneration (FTLD). We stratified the samples by disease status, and discovered that this variant modulates neuronal proportion not only in AD cases, but also several neurodegenerative diseases and in elderly cognitively healthy controls. Furthermore, we did not find a significant association in younger controls or schizophrenia patients, suggesting that this variant might increase neuronal survival or confer resilience to the neurodegenerative process. The single variant and gene-based analyses also identified an overall genetic association between neuronal proportion, AD and FTLD risk. These results suggest that common pathways are implicated in these neurodegenerative diseases, that implicate neuronal survival. In summary, we identified a protective variant in the TMEM106B gene that may have a neuronal protection effect against general aging, independent of disease status, which could help elucidate the relationship between aging and neuronal survival in the presence or absence of neurodegenerative disorders. Our findings suggest that TMEM106B could be a potential target for neuronal protection therapies to ameliorate cognitive and functional deficits.

An atlas of cortical circular RNA expression in Alzheimer disease brains demonstrates clinical and pathological associations

We generated parietal cortex RNA-seq data from individuals with and without Alzheimer disease (AD; n_control = 13; n_AD = 83) from the Knight-ADRC. Using this and an independent (MSBB) AD RNA-seq dataset, we quantified cortical circular RNA (circRNA) expression in the context of AD. We identified significant associations between circRNA expression and AD diagnosis, clinical dementia severity, and neuropathological severity. We demonstrated that a majority of circRNA AD-associations are independent from changes in cognate linear mRNA expression or brain cell-type proportions. We provided evidence for circRNA expression changes occurring early in pre-symptomatic AD, and in autosomal dominant AD. We also observed AD-associated circRNAs co-expressing with known AD genes. Finally, we identified potential microRNA binding sites in AD-associated circRNAs for microRNAs predicted to target AD genes. Together, these results highlight the importance of analyzing non-linear RNAs and support future studies exploring the potential roles of circRNAs in AD pathogenesis.

The MS4A gene cluster is a key modulator of soluble TREM2 and Alzheimer's disease risk

Soluble triggering receptor expressed on myeloid cells 2 (sTREM2) in cerebrospinal fluid (CSF) has been associated with Alzheimer's disease (AD). TREM2 plays a critical role in microglial activation, survival, and phagocytosis; however, the pathophysiological role of sTREM2 in AD is not well understood. Understanding the role of sTREM2 in AD may reveal new pathological mechanisms and lead to the identification of therapeutic targets. We performed a genome-wide association study (GWAS) to identify genetic modifiers of CSF sTREM2 obtained from the Alzheimer's Disease Neuroimaging Initiative. Common variants in the membrane-spanning 4-domains subfamily A (MS4A) gene region were associated with CSF sTREM2 concentrations (rs1582763; P = 1.15 × 10-15); this was replicated in independent datasets. The variants associated with increased CSF sTREM2 concentrations were associated with reduced AD risk and delayed age at onset of disease. The single-nucleotide polymorphism rs1582763 modified expression of the MS4A4A and MS4A6A genes in multiple tissues, suggesting that one or both of these genes are important for modulating sTREM2 production. Using human macrophages as a proxy for microglia, we found that MS4A4A and TREM2 colocalized on lipid rafts at the plasma membrane, that sTREM2 increased with MS4A4A overexpression, and that silencing of MS4A4A reduced sTREM2 production. These genetic, molecular, and cellular findings suggest that MS4A4A modulates sTREM2. These findings also provide a mechanistic explanation for the original GWAS signal in the MS4A locus for AD risk and indicate that TREM2 may be involved in AD pathogenesis not only in TREM2 risk-variant carriers but also in those with sporadic disease.

A single-nuclei RNA sequencing study of Mendelian and sporadic AD in the human brain

Background

Alzheimer’s disease (AD) is the most common form of dementia. This neurodegenerative disorder is associated with neuronal death and gliosis heavily impacting the cerebral cortex. AD has a substantial but heterogeneous genetic component, presenting both Mendelian and complex genetic architectures. Using bulk RNA-seq from the parietal lobes and deconvolution methods, we previously reported that brains exhibiting different AD genetic architecture exhibit different cellular proportions. Here, we sought to directly investigate AD brain changes in cell proportion and gene expression using single-cell resolution.

Methods

We generated unsorted single-nuclei RNA sequencing data from brain tissue. We leveraged the tissue donated from a carrier of a Mendelian genetic mutation, PSEN1 p.A79V, and two family members who suffer from sporadic AD, but do not carry any autosomal mutations. We evaluated alternative alignment approaches to maximize the titer of reads, genes, and cells with high quality. In addition, we employed distinct clustering strategies to determine the best approach to identify cell clusters that reveal neuronal and glial cell types and avoid artifacts such as sample and batch effects. We propose an approach to cluster cells that reduces biases and enable further analyses.

Results

We identified distinct types of neurons, both excitatory and inhibitory, and glial cells, including astrocytes, oligodendrocytes, and microglia, among others. In particular, we identified a reduced proportion of excitatory neurons in the Mendelian mutation carrier, but a similar distribution of inhibitory neurons. Furthermore, we investigated whether single-nuclei RNA-seq from the human brains recapitulate the expression profile of disease-associated microglia (DAM) discovered in mouse models. We also determined that when analyzing human single-nuclei data, it is critical to control for biases introduced by donor-specific expression profiles.

Conclusion

We propose a collection of best practices to generate a highly detailed molecular cell atlas of highly informative frozen tissue stored in brain banks. Importantly, we have developed a new web application to make this unique single-nuclei molecular atlas publicly available.

Electronic supplementary material

The online version of this article (10.1186/s13195-019-0524-x) contains supplementary material, which is available to authorized users.

TREM2 brain transcript-specific studies in AD and TREM2 mutation carriers

BACKGROUND

Low frequency coding variants in TREM2 are associated with Alzheimer disease (AD) risk and cerebrospinal fluid (CSF) TREM2 protein levels are different between AD cases and controls. Similarly, TREM2 risk variant carriers also exhibit differential CSF TREM2 levels. TREM2 has three different alternative transcripts, but most of the functional studies only model the longest transcript. No studies have analyzed TREM2 expression levels or alternative splicing in brains from AD and cognitively normal individuals. We wanted to determine whether there was differential expression of TREM2 in sporadic-AD cases versus AD-TREM2 carriers vs sex- and aged-matched normal controls; and if this differential expression was due to a particular TREM2 transcript.

METHODS

We analyzed RNA-Seq data from parietal lobe brain tissue from AD cases with TREM2 variants (n = 33), AD cases (n = 195) and healthy controls (n = 118), from three independent datasets using Kallisto and the R package tximport to determine the read count for each transcript and quantified transcript abundance as transcripts per million.

RESULTS

The three TREM2 transcripts were expressed in brain cortex in the three datasets. We demonstrate for the first time that the transcript that lacks the transmembrane domain and encodes a soluble form of TREM2 (sTREM2) has an expression level around 60% of the canonical transcript, suggesting that around 25% of the sTREM2 protein levels could be explained by this transcript. We did not observe a difference in the overall TREM2 expression level between cases and controls. However, the isoform which lacks the 5' exon, but includes the transmembrane domain, was significantly lower in TREM2- p.R62H carriers than in AD cases (p = 0.007).

CONCLUSION

Using bulk RNA-Seq data from three different cohorts, we were able to quantify the expression level of the three TREM2 transcripts, demonstrating: (1) all three transcripts of them are highly expressed in the human cortex, (2) that up to 25% of the sTREM2 may be due to the expression of a specific isoform and not TREM2 cleavage; and (3) that TREM2 risk variants do not affect expression levels, suggesting that the effect of the TREM2 variants on CSF levels occurs at post-transcriptional level.

Association Between Aldehyde dehydrogenase-2 Polymorphisms and Risk of Alzheimer's Disease and Parkinson's Disease: A Meta-Analysis Based on 5,315 Individuals

Objective: A number of studies have reported that aldehyde dehydrogenase-2 (ALDH2) polymorphisms maybe associated with the risk of Alzheimer's disease (AD) and Parkinson's disease (PD). However, the results of such studies are inconsistent. We therefore conducted a meta-analysis to clarify the association between ALDH2 polymorphisms and the risk of AD and PD.

Methods: Five online databases were searched and the relevant studies were reviewed from inception through May 10, 2018. Odds ratios (ORs) and corresponding 95% confidence intervals (CIs) were calculated in each genetic model of the general population and various subgroups. Furthermore, we simultaneously performed heterogeneity, cumulative, sensitivity, and publication bias analyses.

Results: Overall, nine case-control studies involving 5,315 subjects were included in this meta-analysis. Potential associations were found between the ALDH2 rs671 G>A polymorphism and the risk of AD (A vs. G: OR = 1.46, 95%CI = 1.01–2.11, P = 0.05, I² = 84.2%; AA vs. GG: OR = 2.22, 95%CI = 1.03–4.77, P = 0.04, I² = 79.2%; AA vs. GG+GA: OR = 1.94, 95%CI = 1.03–3.64, P =0.04, I² = 71.1%). In addition, some similar results were observed in other subgroups. Moreover, no significant association between ALDH2 polymorphisms and PD risk.

Conclusions: In conclusion, our meta-analysis indicated that the ALDH2 rs671 G>A polymorphism plays an important role in AD development.

Analysis of whole genome-transcriptomic organization in brain to identify genes associated with alcoholism

Alcohol exposure triggers changes in gene expression and biological pathways in human brain. We explored alterations in gene expression in the Pre-Frontal Cortex (PFC) of 65 alcoholics and 73 controls of European descent, and identified 129 genes that showed altered expression (FDR < 0.05) in subjects with alcohol dependence. Differentially expressed genes were enriched for pathways related to interferon signaling and Growth Arrest and DNA Damage-inducible 45 (GADD45) signaling. A coexpression module (thistle2) identified by weighted gene co-expression network analysis (WGCNA) was significantly correlated with alcohol dependence, alcohol consumption, and AUDIT scores. Genes in the thistle2 module were enriched with genes related to calcium signaling pathways and showed significant downregulation of these pathways, as well as enrichment for biological processes related to nicotine response and opioid signaling. A second module (brown4) showed significant upregulation of pathways related to immune signaling. Expression quantitative trait loci (eQTLs) for genes in the brown4 module were also enriched for genetic associations with alcohol dependence and alcohol consumption in large genome-wide studies included in the Psychiatric Genetic Consortium and the UK Biobank's alcohol consumption dataset. By leveraging multi-omics data, this transcriptome analysis has identified genes and biological pathways that could provide insight for identifying therapeutic targets for alcohol dependence.

Effects of APOE Genotype on Brain Proteomic Network and Cell Type Changes in Alzheimer's Disease

Polymorphic alleles in the apolipoprotein E (APOE) gene are the main genetic determinants of late-onset Alzheimer's disease (AD) risk. Individuals carrying the APOE E4 allele are at increased risk to develop AD compared to those carrying the more common E3 allele, whereas those carrying the E2 allele are at decreased risk for developing AD. How ApoE isoforms influence risk for AD remains unclear. To help fill this gap in knowledge, we performed a comparative unbiased mass spectrometry-based proteomic analysis of post-mortem brain cortical tissues from pathologically-defined AD or control cases of different APOE genotypes. Control cases (n = 10) were homozygous for the common E3 allele, whereas AD cases (n = 24) were equally distributed among E2/3, E3/3, and E4/4 genotypes. We used differential protein expression and co-expression analytical approaches to assess how changes in the brain proteome are related to APOE genotype. We observed similar levels of amyloid-β, but reduced levels of neurofibrillary tau, in E2/3 brains compared to E3/3 and E4/4 AD brains. Weighted co-expression network analysis revealed 33 modules of co-expressed proteins, 12 of which were significantly different by APOE genotype in AD. The modules that were significantly different by APOE genotype were associated with synaptic transmission and inflammation, among other biological processes. Deconvolution and analysis of brain cell type changes revealed that the E2 allele suppressed homeostatic and disease-associated cell type changes in astrocytes, microglia, oligodendroglia, and endothelia. The E2 allele-specific effect on brain cell type changes was validated in a separate cohort of 130 brains. Our systems-level proteomic analyses of AD brain reveal alterations in the brain proteome and brain cell types associated with allelic variants in APOE, and suggest further areas for investigation into the upstream mechanisms that drive ApoE-associated risk for AD.

Integrative system biology analyses of CRISPR-edited iPSC-derived neurons and human brains reveal deficiencies of presynaptic signaling in FTLD and PSP

Mutations in the microtubule-associated protein tau (MAPT) gene cause autosomal dominant frontotemporal lobar degeneration with tau inclusions (FTLD-tau). MAPT p.R406W carriers present clinically with progressive memory loss and neuropathologically with neuronal and glial tauopathy. However, the pathogenic events triggered by the expression of the mutant tau protein remain poorly understood. To identify the genes and pathways that are dysregulated in FTLD-tau, we performed transcriptomic analyses in induced pluripotent stem cell (iPSC)-derived neurons carrying MAPT p.R406W and CRISPR/Cas9-corrected isogenic controls. We found that the expression of the MAPT p.R406W mutation was sufficient to create a significantly different transcriptomic profile compared with that of the isogeneic controls and to cause the differential expression of 328 genes. Sixty-one of these genes were also differentially expressed in the same direction between MAPT p.R406W carriers and pathology-free human control brains. We found that genes differentially expressed in the stem cell models and human brains were enriched for pathways involving gamma-aminobutyric acid (GABA) receptors and pre-synaptic function. The expression of GABA receptor genes, including GABRB2 and GABRG2, were consistently reduced in iPSC-derived neurons and brains from MAPT p.R406W carriers. Interestingly, we found that GABA receptor genes, including GABRB2 and GABRG2, are significantly lower in symptomatic mouse models of tauopathy, as well as in brains with progressive supranuclear palsy. Genome wide association analyses reveal that common variants within GABRB2 are associated with increased risk for frontotemporal dementia (P < 1 × 10^-3). Thus, our systems biology approach, which leverages molecular data from stem cells, animal models, and human brain tissue can reveal novel disease mechanisms. Here, we demonstrate that MAPT p.R406W is sufficient to induce changes in GABA-mediated signaling and synaptic function, which may contribute to the pathogenesis of FTLD-tau and other primary tauopathies.