LRRK2 G2019S variant is associated with transcriptional changes in Parkinson's disease human myeloid cells under proinflammatory environment.

The G2019S mutation in the leucine-rich repeat kinase 2 (LRRK2) gene is a major risk factor for the development of Parkinson's disease (PD). LRRK2, although ubiquitously expressed, is highly abundant in cells of the innate immune system. Given the importance of central and peripheral immune cells in the development of PD, we sought to investigate the consequences of the G2019S mutation on microglial and monocyte transcriptome and function. We have generated large-scale transcriptomic profiles of isogenic human induced microglial cells (iMGLs) and patient derived monocytes carrying the G2019S mutation under baseline culture conditions and following exposure to the proinflammatory factors IFNγ and LPS. We demonstrate that the G2019S mutation exerts a profound impact on the transcriptomic profile of these myeloid cells, and describe corresponding functional differences in iMGLs. The G2019S mutation led to an upregulation in lipid metabolism and phagolysosomal pathway genes in untreated and LPS/IFNγ stimulated iMGLs, which was accompanied by an increased phagocytic capacity of myelin debris. We also identified dysregulation of cell cycle genes, with a downregulation of the E2F4 regulon. Transcriptomic characterization of human-derived monocytes carrying the G2019S mutation confirmed alteration in lipid metabolism associated genes. Altogether, these findings reveal the influence of G2019S on the dysregulation of the myeloid cell transcriptome under proinflammatory conditions.

Alzheimer's disease transcriptional landscape in ex-vivo human microglia.

Microglia are resident immune cells of the brain and are implicated in the etiology of Alzheimer's Disease (AD) and other diseases. Yet the cellular and molecular processes regulating their function throughout the course of the disease are poorly understood. Here, we present the transcriptional landscape of primary microglia from 189 human postmortem brains, including 58 healthy aging individuals and 131 with a range of disease phenotypes, including 63 patients representing the full spectrum of clinical and pathological severity of AD. We identified transcriptional changes associated with multiple AD phenotypes, capturing the severity of dementia and neuropathological lesions. Transcript-level analyses identified additional genes with heterogeneous isoform usage and AD phenotypes. We identified changes in gene-gene coordination in AD, dysregulation of co-expression modules, and disease subtypes with distinct gene expression. Taken together, these data further our understanding of the key role of microglia in AD biology and nominate candidates for therapeutic intervention.

TDP-43 loss induces extensive cryptic polyadenylation in ALS/FTD.

Nuclear depletion and cytoplasmic aggregation of the RNA-binding protein TDP-43 is the hallmark of ALS, occurring in over 97% of cases. A key consequence of TDP-43 nuclear loss is the de-repression of cryptic exons. Whilst TDP-43 regulated cryptic splicing is increasingly well catalogued, cryptic alternative polyadenylation (APA) events, which define the 3' end of last exons, have been largely overlooked, especially when not associated with novel upstream splice junctions. We developed a novel bioinformatic approach to reliably identify distinct APA event types: alternative last exons (ALE), 3'UTR extensions (3'Ext) and intronic polyadenylation (IPA) events. We identified novel neuronal cryptic APA sites induced by TDP-43 loss of function by systematically applying our pipeline to a compendium of publicly available and in house datasets. We find that TDP-43 binding sites and target motifs are enriched at these cryptic events and that TDP-43 can have both repressive and enhancing action on APA. Importantly, all categories of cryptic APA can also be identified in ALS and FTD post mortem brain regions with TDP-43 proteinopathy underlining their potential disease relevance. RNA-seq and Ribo-seq analyses indicate that distinct cryptic APA categories have different downstream effects on transcript and translation. Intriguingly, cryptic 3'Exts occur in multiple transcription factors, such as ELK1, SIX3, and TLX1, and lead to an increase in wild-type protein levels and function. Finally, we show that an increase in RNA stability leading to a higher cytoplasmic localisation underlies these observations. In summary, we demonstrate that TDP-43 nuclear depletion induces a novel category of cryptic RNA processing events and we expand the palette of TDP-43 loss consequences by showing this can also lead to an increase in normal protein translation.

Long-read RNA-seq atlas of novel microglia isoforms elucidates disease-associated genetic regulation of splicing.

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. We previously mapped the genetic regulation of gene expression and mRNA splicing in human microglia, identifying several loci where common genetic variants in microglia-specific regulatory elements explain disease risk loci identified by GWAS. However, identifying genetic effects on splicing has been challenging due to the use of short sequencing reads to identify causal isoforms. Here we present the isoform-centric microglia genomic atlas (isoMiGA) which leverages the power of long-read RNA-seq to identify 35,879 novel microglia isoforms. We show that the novel microglia isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ethnic meta-analysis of 555 human microglia short-read RNA-seq samples from 391 donors, the largest to date, and found associations with genetic risk loci in Alzheimer's disease and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice site usage.

The human microglia responsome: a resource to better understand microglia states in health and disease.

Microglia, the immune cells of the brain, are increasingly implicated in neurodegenerative disorders through genetic studies. However, how genetic risk factors for these diseases are related to microglial gene expression, microglial function, and ultimately disease, is still largely unknown. Microglia change rapidly in response to alterations in their cellular environment, which is regulated through changes in transcriptional programs, which are as yet poorly understood. Here, we compared the effects of a set of inflammatory and restorative stimuli (lipopolysaccharide, interferon-gamma, resiquimod, tumor necrosis factor-alpha, adenosine triphosphate, dexamethasone, and interleukin-4) on human microglial cells from 67 different donors (N = 398 samples) at the gene and transcript level. We show that microglia from different anatomical brain regions show distinct responses to inflammatory stimuli. We observed a greater overlap between human stimulated microglia and human monocytes than with mouse microglia. We define specific microglial signatures across conditions which are highly relevant for a wide range of biological functions and complex human diseases. Finally, we used our stimulation signatures to interpret associations from Alzheimer's disease (AD) genetic studies and microglia by integrating our inflammatory gene expression profiles with common genetic variants to map cis -expression QTLs (eQTLs). Together, we provide the most comprehensive transcriptomic database of the human microglia responsome. Highlights: RNA-sequencing of 398 human microglial samples exposed to six different triggers.Microglia from different anatomical regions show distinct stimulation responses.Responses in human microglia show a greater overlap with human monocytes than murine microglia.Mapping of response Quantitative Trait Loci identifies interactions between genotype and effect of stimulation on gene expression.Our atlas provides a reference map for interpreting microglia signatures in health and disease.

Brain DNA methylomic analysis of frontotemporal lobar degeneration reveals OTUD4 in shared dysregulated signatures across pathological subtypes.

Frontotemporal lobar degeneration (FTLD) is an umbrella term describing the neuropathology of a clinically, genetically and pathologically heterogeneous group of diseases, including frontotemporal dementia (FTD) and progressive supranuclear palsy (PSP). Among the major FTLD pathological subgroups, FTLD with TDP-43 positive inclusions (FTLD-TDP) and FTLD with tau-positive inclusions (FTLD-tau) are the most common, representing about 90% of the cases. Although alterations in DNA methylation have been consistently associated with neurodegenerative diseases, including Alzheimer's disease and Parkinson's disease, little is known for FTLD and its heterogeneous subgroups and subtypes. The main goal of this study was to investigate DNA methylation variation in FTLD-TDP and FTLD-tau. We used frontal cortex genome-wide DNA methylation profiles from three FTLD cohorts (142 FTLD cases and 92 controls), generated using the Illumina 450K or EPIC microarrays. We performed epigenome-wide association studies (EWAS) for each cohort followed by meta-analysis to identify shared differentially methylated loci across FTLD subgroups/subtypes. In addition, we used weighted gene correlation network analysis to identify co-methylation signatures associated with FTLD and other disease-related traits. Wherever possible, we also incorporated relevant gene/protein expression data. After accounting for a conservative Bonferroni multiple testing correction, the EWAS meta-analysis revealed two differentially methylated loci in FTLD, one annotated to OTUD4 (5'UTR-shore) and the other to NFATC1 (gene body-island). Of these loci, OTUD4 showed consistent upregulation of mRNA and protein expression in FTLD. In addition, in the three independent co-methylation networks, OTUD4-containing modules were enriched for EWAS meta-analysis top loci and were strongly associated with the FTLD status. These co-methylation modules were enriched for genes implicated in the ubiquitin system, RNA/stress granule formation and glutamatergic synaptic signalling. Altogether, our findings identified novel FTLD-associated loci, and support a role for DNA methylation as a mechanism involved in the dysregulation of biological processes relevant to FTLD, highlighting novel potential avenues for therapeutic development.

Integrative transcriptomic analysis of the amyotrophic lateral sclerosis spinal cord implicates glial activation and suggests new risk genes.

Amyotrophic lateral sclerosis (ALS) is a progressively fatal neurodegenerative disease affecting motor neurons in the brain and spinal cord. In this study, we investigated gene expression changes in ALS via RNA sequencing in 380 postmortem samples from cervical, thoracic and lumbar spinal cord segments from 154 individuals with ALS and 49 control individuals. We observed an increase in microglia and astrocyte gene expression, accompanied by a decrease in oligodendrocyte gene expression. By creating a gene co-expression network in the ALS samples, we identified several activated microglia modules that negatively correlate with retrospective disease duration. We mapped molecular quantitative trait loci and found several potential ALS risk loci that may act through gene expression or splicing in the spinal cord and assign putative cell types for FNBP1, ACSL5, SH3RF1 and NFASC. Finally, we outline how common genetic variants associated with splicing of C9orf72 act as proxies for the well-known repeat expansion, and we use the same mechanism to suggest ATXN3 as a putative risk gene.

Transcriptome deregulation of peripheral monocytes and whole blood in GBA-related Parkinson's disease.

BACKGROUND: Genetic mutations in beta-glucocerebrosidase (GBA) represent the major genetic risk factor for Parkinson's disease (PD). GBA participates in both the endo-lysosomal pathway and the immune response, two important mechanisms involved in the pathogenesis of PD. However, modifiers of GBA penetrance have not yet been fully elucidated. METHODS: We characterized the transcriptomic profiles of circulating monocytes in a population of patients with PD and healthy controls (CTRL) with and without GBA variants (n = 23 PD/GBA, 13 CTRL/GBA, 56 PD, 66 CTRL) and whole blood (n = 616 PD, 362 CTRL, 127 PD/GBA, 165 CTRL/GBA). Differential expression analysis, pathway enrichment analysis, and outlier detection were performed. Ultrastructural characterization of isolated CD14+ monocytes in the four groups was also performed through electron microscopy. RESULTS: We observed hundreds of differentially expressed genes and dysregulated pathways when comparing manifesting and non-manifesting GBA mutation carriers. Specifically, when compared to idiopathic PD, PD/GBA showed dysregulation in genes involved in alpha-synuclein degradation, aging and amyloid processing. Gene-based outlier analysis confirmed the involvement of lysosomal, membrane trafficking, and mitochondrial processing in manifesting compared to non-manifesting GBA-carriers, as also observed at the ultrastructural levels. Transcriptomic results were only partially replicated in an independent cohort of whole blood samples, suggesting cell-type specific changes. CONCLUSIONS: Overall, our transcriptomic analysis of primary monocytes identified gene targets and biological processes that can help in understanding the pathogenic mechanisms associated with GBA mutations in the context of PD.

Genetics of the human microglia regulome refines Alzheimer's disease risk loci.

Microglia are brain myeloid cells that play a critical role in neuroimmunity and the etiology of Alzheimer's disease (AD), yet our understanding of how the genetic regulatory landscape controls microglial function and contributes to AD is limited. Here, we performed transcriptome and chromatin accessibility profiling in primary human microglia from 150 donors to identify genetically driven variation and cell-specific enhancer-promoter (E-P) interactions. Integrative fine-mapping analysis identified putative regulatory mechanisms for 21 AD risk loci, of which 18 were refined to a single gene, including 3 new candidate risk genes (KCNN4, FIBP and LRRC25). Transcription factor regulatory networks captured AD risk variation and identified SPI1 as a key putative regulator of microglia expression and AD risk. This comprehensive resource capturing variation in the human microglia regulome provides insights into the etiology of neurodegenerative disease.

17q21.31 sub-haplotypes underlying H1-associated risk for Parkinson's disease are associated with LRRC37A/2 expression in astrocytes.

BACKGROUND: Parkinson's disease (PD) is genetically associated with the H1 haplotype of the MAPT 17q.21.31 locus, although the causal gene and variants underlying this association have not been identified. METHODS: To better understand the genetic contribution of this region to PD and to identify novel mechanisms conferring risk for the disease, we fine-mapped the 17q21.31 locus by constructing discrete haplotype blocks from genetic data. We used digital PCR to assess copy number variation associated with PD-associated blocks, and used human brain postmortem RNA-seq data to identify candidate genes that were then further investigated using in vitro models and human brain tissue. RESULTS: We identified three novel H1 sub-haplotype blocks across the 17q21.31 locus associated with PD risk. Protective sub-haplotypes were associated with increased LRRC37A/2 copy number and expression in human brain tissue. We found that LRRC37A/2 is a membrane-associated protein that plays a role in cellular migration, chemotaxis and astroglial inflammation. In human substantia nigra, LRRC37A/2 was primarily expressed in astrocytes, interacted directly with soluble α-synuclein, and co-localized with Lewy bodies in PD brain tissue. CONCLUSION: These data indicate that a novel candidate gene, LRRC37A/2, contributes to the association between the 17q21.31 locus and PD via its interaction with α-synuclein and its effects on astrocytic function and inflammatory response. These data are the first to associate the genetic association at the 17q21.31 locus with PD pathology, and highlight the importance of variation at the 17q21.31 locus in the regulation of multiple genes other than MAPT and KANSL1, as well as its relevance to non-neuronal cell types.

Integrating whole-genome sequencing with multi-omic data reveals the impact of structural variants on gene regulation in the human brain.

Structural variants (SVs), which are genomic rearrangements of more than 50 base pairs, are an important source of genetic diversity and have been linked to many diseases. However, it remains unclear how they modulate human brain function and disease risk. Here we report 170,996 SVs discovered using 1,760 short-read whole genomes from aged adults and individuals with Alzheimer's disease. By applying quantitative trait locus (SV-xQTL) analyses, we quantified the impact of cis-acting SVs on histone modifications, gene expression, splicing and protein abundance in postmortem brain tissues. More than 3,200 SVs were associated with at least one molecular phenotype. We found reproducibility of 65-99% SV-eQTLs across cohorts and brain regions. SV associations with mRNA and proteins shared the same direction of effect in more than 87% of SV-gene pairs. Mediation analysis showed ~8% of SV-eQTLs mediated by histone acetylation and ~11% by splicing. Additionally, associations of SVs with progressive supranuclear palsy identified previously known and novel SVs.

TDP-43 loss and ALS-risk SNPs drive mis-splicing and depletion of UNC13A.

Variants of UNC13A, a critical gene for synapse function, increase the risk of amyotrophic lateral sclerosis and frontotemporal dementia1-3, two related neurodegenerative diseases defined by mislocalization of the RNA-binding protein TDP-434,5. Here we show that TDP-43 depletion induces robust inclusion of a cryptic exon in UNC13A, resulting in nonsense-mediated decay and loss of UNC13A protein. Two common intronic UNC13A polymorphisms strongly associated with amyotrophic lateral sclerosis and frontotemporal dementia risk overlap with TDP-43 binding sites. These polymorphisms potentiate cryptic exon inclusion, both in cultured cells and in brains and spinal cords from patients with these conditions. Our findings, which demonstrate a genetic link between loss of nuclear TDP-43 function and disease, reveal the mechanism by which UNC13A variants exacerbate the effects of decreased TDP-43 function. They further provide a promising therapeutic target for TDP-43 proteinopathies.

Genetic analysis of the human microglial transcriptome across brain regions, aging and disease pathologies.

Microglia have emerged as important players in brain aging and pathology. To understand how genetic risk for neurological and psychiatric disorders is related to microglial function, large transcriptome studies are essential. Here we describe the transcriptome analysis of 255 primary human microglial samples isolated at autopsy from multiple brain regions of 100 individuals. We performed systematic analyses to investigate various aspects of microglial heterogeneities, including brain region and aging. We mapped expression and splicing quantitative trait loci and showed that many neurological disease susceptibility loci are mediated through gene expression or splicing in microglia. Fine-mapping of these loci nominated candidate causal variants that are within microglia-specific enhancers, finding associations with microglial expression of USP6NL for Alzheimer's disease and P2RY12 for Parkinson's disease. We have built the most comprehensive catalog to date of genetic effects on the microglial transcriptome and propose candidate functional variants in neurological and psychiatric disorders.

Transcriptomic analysis of frontotemporal lobar degeneration with TDP-43 pathology reveals cellular alterations across multiple brain regions.

Frontotemporal lobar degeneration (FTLD) is a group of heterogeneous neurodegenerative disorders affecting the frontal and temporal lobes of the brain. Nuclear loss and cytoplasmic aggregation of the RNA-binding protein TDP-43 represents the major FTLD pathology, known as FTLD-TDP. To date, there is no effective treatment for FTLD-TDP due to an incomplete understanding of the molecular mechanisms underlying disease development. Here we compared postmortem tissue RNA-seq transcriptomes from the frontal cortex, temporal cortex, and cerebellum between 28 controls and 30 FTLD-TDP patients to profile changes in cell-type composition, gene expression and transcript usage. We observed downregulation of neuronal markers in all three regions of the brain, accompanied by upregulation of microglia, astrocytes, and oligodendrocytes, as well as endothelial cells and pericytes, suggesting shifts in both immune activation and within the vasculature. We validate our estimates of neuronal loss using neuropathological atrophy scores and show that neuronal loss in the cortex can be mainly attributed to excitatory neurons, and that increases in microglial and endothelial cell expression are highly correlated with neuronal loss. All our analyses identified a strong involvement of the cerebellum in the neurodegenerative process of FTLD-TDP. Altogether, our data provides a detailed landscape of gene expression alterations to help unravel relevant disease mechanisms in FTLD.

Multi-omic insights into Parkinson's Disease: From genetic associations to functional mechanisms.

Genome-Wide Association Studies (GWAS) have elucidated the genetic components of Parkinson's Disease (PD). However, because the vast majority of GWAS association signals fall within non-coding regions, translating these results into an interpretable, mechanistic understanding of the disease etiology remains a major challenge in the field. In this review, we provide an overview of the approaches to prioritize putative causal variants and genes as well as summarise the primary findings of previous studies. We then discuss recent efforts to integrate multi-omics data to identify likely pathogenic cell types and biological pathways implicated in PD pathogenesis. We have compiled full summary statistics of cell-type, tissue, and phentoype enrichment analyses from multiple studies of PD GWAS and provided them in a standardized format as a resource for the research community (https://github.com/RajLabMSSM/PD_omics_review). Finally, we discuss the experimental, computational, and conceptual advances that will be necessary to fully elucidate the effects of functional variants and genes on cellular dysregulation and disease risk.

Dysregulated coordination of MAPT exon 2 and exon 10 splicing underlies different tau pathologies in PSP and AD.

Understanding regulation of MAPT splicing is important to the etiology of many nerurodegenerative diseases, including Alzheimer disease (AD) and progressive supranuclear palsy (PSP), in which different tau isoforms accumulate in pathologic inclusions. MAPT, the gene encoding the tau protein, undergoes complex alternative pre-mRNA splicing to generate six isoforms. Tauopathies can be categorized by the presence of tau aggregates containing either 3 (3R) or 4 (4R) microtubule-binding domain repeats (determined by inclusion/exclusion of exon 10), but the role of the N-terminal domain of the protein, determined by inclusion/exclusion of exons 2 and 3 has been less well studied. Using a correlational screen in human brain tissue, we observed coordination of MAPT exons 2 and 10 splicing. Expressions of exon 2 splicing regulators and subsequently exon 2 inclusion are differentially disrupted in PSP and AD brain, resulting in the accumulation of 1N4R isoforms in PSP and 0N isoforms in AD temporal cortex. Furthermore, we identified different N-terminal isoforms of tau present in neurofibrillary tangles, dystrophic neurites and tufted astrocytes, indicating a role for differential N-terminal splicing in the development of disparate tau neuropathologies. We conclude that N-terminal splicing and combinatorial regulation with exon 10 inclusion/exclusion is likely to be important to our understanding of tauopathies.

Fine-mapping of Parkinson's disease susceptibility loci identifies putative causal variants.

Recent genome-wide association studies have identified 78 loci associated with Parkinson's disease susceptibility but the underlying mechanisms remain largely unclear. To identify likely causal variants for disease risk, we fine-mapped these Parkinson's-associated loci using four different fine-mapping methods. We then integrated multi-assay cell type-specific epigenomic profiles to pinpoint the likely mechanism of action of each variant, allowing us to identify Consensus single nucleotide polymorphism (SNPs) that disrupt LRRK2 and FCGR2A regulatory elements in microglia, an MBNL2 enhancer in oligodendrocytes, and a DYRK1A enhancer in neurons. This genome-wide functional fine-mapping investigation of Parkinson's disease substantially advances our understanding of the causal mechanisms underlying this complex disease while avoiding focus on spurious, non-causal mechanisms. Together, these results provide a robust, comprehensive list of the likely causal variants, genes and cell-types underlying Parkinson's disease risk as demonstrated by consistently greater enrichment of our fine-mapped SNPs relative to lead GWAS SNPs across independent functional impact annotations. In addition, our approach prioritized an average of 3/85 variants per locus as putatively causal, making downstream experimental studies both more tractable and more likely to yield disease-relevant, actionable results. Large-scale studies comparing individuals with Parkinson's disease to age-matched controls have identified many regions of the genome associated with the disease. However, there is widespread correlation between different parts of the genome, making it difficult to tell which genetic variants cause Parkinson's and which are simply co-inherited with causal variants. We therefore applied a suite of statistical models to identify the most likely causal genetic variants (i.e. fine-mapping). We then linked these genetic variants with epigenomic and gene expression signatures across a wide variety of tissues and cell types to identify how these variants cause disease. Therefore, this study provides a comprehensive and robust list of cellular and molecular mechanisms that may serve as targets in the development of more effective Parkinson's therapeutics.

echolocatoR: an automated end-to-end statistical and functional genomic fine-mapping pipeline.

SUMMARY: echolocatoR integrates a diverse suite of statistical and functional fine-mapping tools to identify, test enrichment in, and visualize high-confidence causal consensus variants in any phenotype. It requires minimal input from users (a summary statistics file), can be run in a single R function, and provides extensive access to relevant datasets (e.g. reference linkage disequilibrium panels, quantitative trait loci, genome-wide annotations, cell-type-specific epigenomics), thereby enabling rapid, robust and scalable end-to-end fine-mapping investigations. AVAILABILITY AND IMPLEMENTATION: echolocatoR is an open-source R package available through GitHub under the GNU General Public License (Version 3) license: https://github.com/RajLabMSSM/echolocatoR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

HnRNP K mislocalisation is a novel protein pathology of frontotemporal lobar degeneration and ageing and leads to cryptic splicing.

Heterogeneous nuclear ribonucleoproteins (HnRNPs) are a group of ubiquitously expressed RNA-binding proteins implicated in the regulation of all aspects of nucleic acid metabolism. HnRNP K is a member of this highly versatile hnRNP family. Pathological redistribution of hnRNP K to the cytoplasm has been linked to the pathogenesis of several malignancies but, until now, has been underexplored in the context of neurodegenerative disease. Here we show hnRNP K mislocalisation in pyramidal neurons of the frontal cortex to be a novel neuropathological feature that is associated with both frontotemporal lobar degeneration and ageing. HnRNP K mislocalisation is mutually exclusive to TDP-43 and tau pathological inclusions in neurons and was not observed to colocalise with mitochondrial, autophagosomal or stress granule markers. De-repression of cryptic exons in RNA targets following TDP-43 nuclear depletion is an emerging mechanism of potential neurotoxicity in frontotemporal lobar degeneration and the mechanistically overlapping disorder amyotrophic lateral sclerosis. We silenced hnRNP K in neuronal cells to identify the transcriptomic consequences of hnRNP K nuclear depletion. Intriguingly, by performing RNA-seq analysis we find that depletion of hnRNP K induces 101 novel cryptic exon events. We validated cryptic exon inclusion in an SH-SY5Y hnRNP K knockdown and in FTLD brain exhibiting hnRNP K nuclear depletion. We, therefore, present evidence for hnRNP K mislocalisation to be associated with FTLD and for this to induce widespread changes in splicing.