MIR-NATs repress MAPT translation and aid proteostasis in neurodegeneration.

The human genome expresses thousands of natural antisense transcripts (NAT) that can regulate epigenetic state, transcription, RNA stability or translation of their overlapping genes1,2. Here we describe MAPT-AS1, a brain-enriched NAT that is conserved in primates and contains an embedded mammalian-wide interspersed repeat (MIR), which represses tau translation by competing for ribosomal RNA pairing with the MAPT mRNA internal ribosome entry site3. MAPT encodes tau, a neuronal intrinsically disordered protein (IDP) that stabilizes axonal microtubules. Hyperphosphorylated, aggregation-prone tau forms the hallmark inclusions of tauopathies4. Mutations in MAPT cause familial frontotemporal dementia, and common variations forming the MAPT H1 haplotype are a significant risk factor in many tauopathies5 and Parkinson's disease. Notably, expression of MAPT-AS1 or minimal essential sequences from MAPT-AS1 (including MIR) reduces-whereas silencing MAPT-AS1 expression increases-neuronal tau levels, and correlate with tau pathology in human brain. Moreover, we identified many additional NATs with embedded MIRs (MIR-NATs), which are overrepresented at coding genes linked to neurodegeneration and/or encoding IDPs, and confirmed MIR-NAT-mediated translational control of one such gene, PLCG1. These results demonstrate a key role for MAPT-AS1 in tauopathies and reveal a potentially broad contribution of MIR-NATs to the tightly controlled translation of IDPs6, with particular relevance for proteostasis in neurodegeneration.

eQTL Catalogue 2023: New datasets, X chromosome QTLs, and improved detection and visualisation of transcript-level QTLs.

The eQTL Catalogue is an open database of uniformly processed human molecular quantitative trait loci (QTLs). We are continuously updating the resource to further increase its utility for interpreting genetic associations with complex traits. Over the past two years, we have increased the number of uniformly processed studies from 21 to 31 and added X chromosome QTLs for 19 compatible studies. We have also implemented Leafcutter to directly identify splice-junction usage QTLs in all RNA sequencing datasets. Finally, to improve the interpretability of transcript-level QTLs, we have developed static QTL coverage plots that visualise the association between the genotype and average RNA sequencing read coverage in the region for all 1.7 million fine mapped associations. To illustrate the utility of these updates to the eQTL Catalogue, we performed colocalisation analysis between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. Although most GWAS loci colocalised both with eQTLs and transcript-level QTLs, we found that visual inspection could sometimes be used to distinguish primary splicing QTLs from those that appear to be secondary consequences of large-effect gene expression QTLs. While these visually confirmed primary splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases.

IntroVerse: a comprehensive database of introns across human tissues.

Dysregulation of RNA splicing contributes to both rare and complex diseases. RNA-sequencing data from human tissues has shown that this process can be inaccurate, resulting in the presence of novel introns detected at low frequency across samples and within an individual. To enable the full spectrum of intron use to be explored, we have developed IntroVerse, which offers an extensive catalogue on the splicing of 332,571 annotated introns and a linked set of 4,679,474 novel junctions covering 32,669 different genes. This dataset has been generated through the analysis of 17,510 human control RNA samples from 54 tissues provided by the Genotype-Tissue Expression Consortium. IntroVerse has two unique features: (i) it provides a complete catalogue of novel junctions and (ii) each novel junction has been assigned to a specific annotated intron. This unique, hierarchical structure offers multiple uses, including the identification of novel transcripts from known genes and their tissue-specific usage, and the assessment of background splicing noise for introns thought to be mis-spliced in disease states. IntroVerse provides a user-friendly web interface and is freely available at https://rytenlab.com/browser/app/introverse.

Developmental Consequences of Defective ATG7-Mediated Autophagy in Humans.

BACKGROUND: Autophagy is the major intracellular degradation route in mammalian cells. Systemic ablation of core autophagy-related (ATG) genes in mice leads to embryonic or perinatal lethality, and conditional models show neurodegeneration. Impaired autophagy has been associated with a range of complex human diseases, yet congenital autophagy disorders are rare. METHODS: We performed a genetic, clinical, and neuroimaging analysis involving five families. Mechanistic investigations were conducted with the use of patient-derived fibroblasts, skeletal muscle-biopsy specimens, mouse embryonic fibroblasts, and yeast. RESULTS: We found deleterious, recessive variants in human ATG7, a core autophagy-related gene encoding a protein that is indispensable to classical degradative autophagy. Twelve patients from five families with distinct ATG7 variants had complex neurodevelopmental disorders with brain, muscle, and endocrine involvement. Patients had abnormalities of the cerebellum and corpus callosum and various degrees of facial dysmorphism. These patients have survived with impaired autophagic flux arising from a diminishment or absence of ATG7 protein. Although autophagic sequestration was markedly reduced, evidence of basal autophagy was readily identified in fibroblasts and skeletal muscle with loss of ATG7. Complementation of different model systems by deleterious ATG7 variants resulted in poor or absent autophagic function as compared with the reintroduction of wild-type ATG7. CONCLUSIONS: We identified several patients with a neurodevelopmental disorder who have survived with a severe loss or complete absence of ATG7, an essential effector enzyme for autophagy without a known functional paralogue. (Funded by the Wellcome Centre for Mitochondrial Research and others.).

Genome-Wide Analysis of Structural Variants in Parkinson Disease.

OBJECTIVE: Identification of genetic risk factors for Parkinson disease (PD) has to date been primarily limited to the study of single nucleotide variants, which only represent a small fraction of the genetic variation in the human genome. Consequently, causal variants for most PD risk are not known. Here we focused on structural variants (SVs), which represent a major source of genetic variation in the human genome. We aimed to discover SVs associated with PD risk by performing the first large-scale characterization of SVs in PD. METHODS: We leveraged a recently developed computational pipeline to detect and genotype SVs from 7,772 Illumina short-read whole genome sequencing samples. Using this set of SV variants, we performed a genome-wide association study using 2,585 cases and 2,779 controls and identified SVs associated with PD risk. Furthermore, to validate the presence of these variants, we generated a subset of matched whole-genome long-read sequencing data. RESULTS: We genotyped and tested 3,154 common SVs, representing over 412 million nucleotides of previously uncatalogued genetic variation. Using long-read sequencing data, we validated the presence of three novel deletion SVs that are associated with risk of PD from our initial association analysis, including a 2 kb intronic deletion within the gene LRRN4. INTERPRETATION: We identified three SVs associated with genetic risk of PD. This study represents the most comprehensive assessment of the contribution of SVs to the genetic risk of PD to date. ANN NEUROL 2023;93:1012-1022.

Adaptive Long-Read Sequencing Reveals GGC Repeat Expansion in ZFHX3 Associated with Spinocerebellar Ataxia Type 4.

BACKGROUND: Spinocerebellar ataxia type 4 (SCA4) is an autosomal dominant ataxia with invariable sensory neuropathy originally described in a family with Swedish ancestry residing in Utah more than 25 years ago. Despite tight linkage to the 16q22 region, the molecular diagnosis has since remained elusive. OBJECTIVES: Inspired by pathogenic structural variation implicated in other 16q-ataxias with linkage to the same locus, we revisited the index SCA4 cases from the Utah family using novel technologies to investigate structural variation within the candidate region. METHODS: We adopted a targeted long-read sequencing approach with adaptive sampling on the Oxford Nanopore Technologies (ONT) platform that enables the detection of segregating structural variants within a genomic region without a priori assumptions about any variant features. RESULTS: Using this approach, we found a heterozygous (GGC)n repeat expansion in the last coding exon of the zinc finger homeobox 3 (ZFHX3) gene that segregates with disease, ranging between 48 and 57 GGC repeats in affected probands. This finding was replicated in a separate family with SCA4. Furthermore, the estimation of this GGC repeat size in short-read whole genome sequencing (WGS) data of 21,836 individuals recruited to the 100,000 Genomes Project in the UK and our in-house dataset of 11,258 exomes did not reveal any pathogenic repeats, indicating that the variant is ultrarare. CONCLUSIONS: These findings support the utility of adaptive long-read sequencing as a powerful tool to decipher causative structural variation in unsolved cases of inherited neurological disease. © 2024 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

The non-specific lethal complex regulates genes and pathways genetically linked to Parkinson's disease.

Genetic variants conferring risks for Parkinson's disease have been highlighted through genome-wide association studies, yet exploration of their specific disease mechanisms is lacking. Two Parkinson's disease candidate genes, KAT8 and KANSL1, identified through genome-wide studies and a PINK1-mitophagy screen, encode part of the histone acetylating non-specific lethal complex. This complex localizes to the nucleus, where it plays a role in transcriptional activation, and to mitochondria, where it has been suggested to have a role in mitochondrial transcription. In this study, we sought to identify whether the non-specific lethal complex has potential regulatory relationships with other genes associated with Parkinson's disease in human brain. Correlation in the expression of non-specific lethal genes and Parkinson's disease-associated genes was investigated in primary gene co-expression networks using publicly-available transcriptomic data from multiple brain regions (provided by the Genotype-Tissue Expression Consortium and UK Brain Expression Consortium), whilst secondary networks were used to examine cell type specificity. Reverse engineering of gene regulatory networks generated regulons of the complex, which were tested for heritability using stratified linkage disequilibrium score regression. Prioritized gene targets were then validated in vitro using a QuantiGene multiplex assay and publicly-available chromatin immunoprecipitation-sequencing data. Significant clustering of non-specific lethal genes was revealed alongside Parkinson's disease-associated genes in frontal cortex primary co-expression modules, amongst other brain regions. Both primary and secondary co-expression modules containing these genes were enriched for mainly neuronal cell types. Regulons of the complex contained Parkinson's disease-associated genes and were enriched for biological pathways genetically linked to disease. When examined in a neuroblastoma cell line, 41% of prioritized gene targets showed significant changes in mRNA expression following KANSL1 or KAT8 perturbation. KANSL1 and H4K8 chromatin immunoprecipitation-sequencing data demonstrated non-specific lethal complex activity at many of these genes. In conclusion, genes encoding the non-specific lethal complex are highly correlated with and regulate genes associated with Parkinson's disease. Overall, these findings reveal a potentially wider role for this protein complex in regulating genes and pathways implicated in Parkinson's disease.

Functional genomics provide key insights to improve the diagnostic yield of hereditary ataxia.

Improvements in functional genomic annotation have led to a critical mass of neurogenetic discoveries. This is exemplified in hereditary ataxia, a heterogeneous group of disorders characterised by incoordination from cerebellar dysfunction. Associated pathogenic variants in more than 300 genes have been described, leading to a detailed genetic classification partitioned by age-of-onset. Despite these advances, up to 75% of patients with ataxia remain molecularly undiagnosed even following whole genome sequencing, as exemplified in the 100 000 Genomes Project. This study aimed to understand whether we can improve our knowledge of the genetic architecture of hereditary ataxia by leveraging functional genomic annotations, and as a result, generate insights and strategies that raise the diagnostic yield. To achieve these aims, we used publicly-available multi-omics data to generate 294 genic features, capturing information relating to a gene's structure, genetic variation, tissue-specific, cell-type-specific and temporal expression, as well as protein products of a gene. We studied these features across genes typically causing childhood-onset, adult-onset or both types of disease first individually, then collectively. This led to the generation of testable hypotheses which we investigated using whole genome sequencing data from up to 2182 individuals presenting with ataxia and 6658 non-neurological probands recruited in the 100 000 Genomes Project. Using this approach, we demonstrated a high short tandem repeat (STR) density within childhood-onset genes suggesting that we may be missing pathogenic repeat expansions within this cohort. This was verified in both childhood- and adult-onset ataxia patients from the 100 000 Genomes Project who were unexpectedly found to have a trend for higher repeat sizes even at naturally-occurring STRs within known ataxia genes, implying a role for STRs in pathogenesis. Using unsupervised analysis, we found significant similarities in genomic annotation across the gene panels, which suggested adult- and childhood-onset patients should be screened using a common diagnostic gene set. We tested this within the 100 000 Genomes Project by assessing the burden of pathogenic variants among childhood-onset genes in adult-onset patients and vice versa. This demonstrated a significantly higher burden of rare, potentially pathogenic variants in conventional childhood-onset genes among individuals with adult-onset ataxia. Our analysis has implications for the current clinical practice in genetic testing for hereditary ataxia. We suggest that the diagnostic rate for hereditary ataxia could be increased by removing the age-of-onset partition, and through a modified screening for repeat expansions in naturally-occurring STRs within known ataxia-associated genes, in effect treating these regions as candidate pathogenic loci.

Association between the LRP1B and APOE loci and the development of Parkinson's disease dementia.

Parkinson's disease is one of the most common age-related neurodegenerative disorders. Although predominantly a motor disorder, cognitive impairment and dementia are important features of Parkinson's disease, particularly in the later stages of the disease. However, the rate of cognitive decline varies among Parkinson's disease patients, and the genetic basis for this heterogeneity is incompletely understood. To explore the genetic factors associated with rate of progression to Parkinson's disease dementia, we performed a genome-wide survival meta-analysis of 3923 clinically diagnosed Parkinson's disease cases of European ancestry from four longitudinal cohorts. In total, 6.7% of individuals with Parkinson's disease developed dementia during study follow-up, on average 4.4 ± 2.4 years from disease diagnosis. We have identified the APOE ε4 allele as a major risk factor for the conversion to Parkinson's disease dementia [hazard ratio = 2.41 (1.94-3.00), P = 2.32 × 10-15], as well as a new locus within the ApoE and APP receptor LRP1B gene [hazard ratio = 3.23 (2.17-4.81), P = 7.07 × 10-09]. In a candidate gene analysis, GBA variants were also identified to be associated with higher risk of progression to dementia [hazard ratio = 2.02 (1.21-3.32), P = 0.007]. CSF biomarker analysis also implicated the amyloid pathway in Parkinson's disease dementia, with significantly reduced levels of amyloid ß42 (P = 0.0012) in Parkinson's disease dementia compared to Parkinson's disease without dementia. These results identify a new candidate gene associated with faster conversion to dementia in Parkinson's disease and suggest that amyloid-targeting therapy may have a role in preventing Parkinson's disease dementia.

Large-scale rare variant burden testing in Parkinson's disease.

Parkinson's disease has a large heritable component and genome-wide association studies have identified over 90 variants with disease-associated common variants, providing deeper insights into the disease biology. However, there have not been large-scale rare variant analyses for Parkinson's disease. To address this gap, we investigated the rare genetic component of Parkinson's disease at minor allele frequencies <1%, using whole genome and whole exome sequencing data from 7184 Parkinson's disease cases, 6701 proxy cases and 51 650 healthy controls from the Accelerating Medicines Partnership Parkinson's disease (AMP-PD) initiative, the National Institutes of Health, the UK Biobank and Genentech. We performed burden tests meta-analyses on small indels and single nucleotide protein-altering variants, prioritized based on their predicted functional impact. Our work identified several genes reaching exome-wide significance. Two of these genes, GBA1 and LRRK2, have variants that have been previously implicated as risk factors for Parkinson's disease, with some variants in LRRK2 resulting in monogenic forms of the disease. We identify potential novel risk associations for variants in B3GNT3, AUNIP, ADH5, TUBA1B, OR1G1, CAPN10 and TREML1 but were unable to replicate the observed associations across independent datasets. Of these, B3GNT3 and TREML1 could provide new evidence for the role of neuroinflammation in Parkinson's disease. To date, this is the largest analysis of rare genetic variants in Parkinson's disease.

Repeat expansion disorders.

An increasing number of repeat expansion disorders have been found to cause both rare and common neurological disease. This is exemplified in recent discoveries of novel repeat expansions underlying a significant proportion of several late-onset neurodegenerative disorders, such as CANVAS (cerebellar ataxia, neuropathy and vestibular areflexia syndrome) and spinocerebellar ataxia type 27B. Most of the 60 described repeat expansion disorders to date are associated with neurological disease, providing substantial challenges for diagnosis, but also opportunities for management in a clinical neurology setting. Commonalities in clinical presentation, overarching diagnostic features and similarities in the approach to genetic testing justify considering these disorders collectively based on their unifying causative mechanism. In this review, we discuss the characteristics and diagnostic challenges of repeat expansion disorders for the neurologist and provide examples to highlight their clinical heterogeneity. With the ready availability of clinical-grade whole-genome sequencing for molecular diagnosis, we discuss the current approaches to testing for repeat expansion disorders and application in clinical practice.

The contribution of Neanderthal introgression and natural selection to neurodegenerative diseases.

Humans are thought to be more susceptible to neurodegeneration than equivalently-aged primates. It is not known whether this vulnerability is specific to anatomically-modern humans or shared with other hominids. The contribution of introgressed Neanderthal DNA to neurodegenerative disorders remains uncertain. It is also unclear how common variants associated with neurodegenerative disease risk are maintained by natural selection in the population despite their deleterious effects. In this study, we aimed to quantify the genome-wide contribution of Neanderthal introgression and positive selection to the heritability of complex neurodegenerative disorders to address these questions. We used stratified-linkage disequilibrium score regression to investigate the relationship between five SNP-based signatures of natural selection, reflecting different timepoints of evolution, and genome-wide associated variants of the three most prevalent neurodegenerative disorders: Alzheimer's disease, amyotrophic lateral sclerosis and Parkinson's disease. We found no evidence for enrichment of positively-selected SNPs in the heritability of Alzheimer's disease, amyotrophic lateral sclerosis and Parkinson's disease, suggesting that common deleterious disease variants are unlikely to be maintained by positive selection. There was no enrichment of Neanderthal introgression in the SNP-heritability of these disorders, suggesting that Neanderthal admixture is unlikely to have contributed to disease risk. These findings provide insight into the origins of neurodegenerative disorders within the evolution of Homo sapiens and addresses a long-standing debate, showing that Neanderthal admixture is unlikely to have contributed to common genetic risk of neurodegeneration in anatomically-modern humans.

ggtranscript: an R package for the visualization and interpretation of transcript isoforms using ggplot2.

MOTIVATION: The advent of long-read sequencing technologies has increased demand for the visualization and interpretation of transcripts. However, tools that perform such visualizations remain inflexible and lack the ability to easily identify differences between transcript structures. Here, we introduce ggtranscript, an R package that provides a fast and flexible method to visualize and compare transcripts. As a ggplot2 extension, ggtranscript inherits the functionality and familiarity of ggplot2 making it easy to use. AVAILABILITY AND IMPLEMENTATION: ggtranscript is an R package available at https://github.com/dzhang32/ggtranscript (DOI: https://doi.org/10.5281/zenodo.6374061) via an open-source MIT licence. Further documentation is available at https://dzhang32.github.io/ggtranscript/.

Regulation of mitophagy by the NSL complex underlies genetic risk for Parkinson's disease at 16q11.2 and MAPT H1 loci.

Parkinson's disease is a common incurable neurodegenerative disease. The identification of genetic variants via genome-wide association studies has considerably advanced our understanding of the Parkinson's disease genetic risk. Understanding the functional significance of the risk loci is now a critical step towards translating these genetic advances into an enhanced biological understanding of the disease. Impaired mitophagy is a key causative pathway in familial Parkinson's disease, but its relevance to idiopathic Parkinson's disease is unclear. We used a mitophagy screening assay to evaluate the functional significance of risk genes identified through genome-wide association studies. We identified two new regulators of PINK1-dependent mitophagy initiation, KAT8 and KANSL1, previously shown to modulate lysine acetylation. These findings suggest PINK1-mitophagy is a contributing factor to idiopathic Parkinson's disease. KANSL1 is located on chromosome 17q21 where the risk associated gene has long been considered to be MAPT. While our data do not exclude a possible association between the MAPT gene and Parkinson's disease, they provide strong evidence that KANSL1 plays a crucial role in the disease. Finally, these results enrich our understanding of physiological events regulating mitophagy and establish a novel pathway for drug targeting in neurodegeneration.

PhenoExam: gene set analyses through integration of different phenotype databases.

BACKGROUND: Gene set enrichment analysis (detecting phenotypic terms that emerge as significant in a set of genes) plays an important role in bioinformatics focused on diseases of genetic basis. To facilitate phenotype-oriented gene set analysis, we developed PhenoExam, a freely available R package for tool developers and a web interface for users, which performs: (1) phenotype and disease enrichment analysis on a gene set; (2) measures statistically significant phenotype similarities between gene sets and (3) detects significant differential phenotypes or disease terms across different databases. RESULTS: PhenoExam generates sensitive and accurate phenotype enrichment analyses. It is also effective in segregating gene sets or Mendelian diseases with very similar phenotypes. We tested the tool with two similar diseases (Parkinson and dystonia), to show phenotype-level similarities but also potentially interesting differences. Moreover, we used PhenoExam to validate computationally predicted new genes potentially associated with epilepsy. CONCLUSIONS: We developed PhenoExam, a freely available R package and Web application, which performs phenotype enrichment and disease enrichment analysis on gene set G, measures statistically significant phenotype similarities between pairs of gene sets G and G' and detects statistically significant exclusive phenotypes or disease terms, across different databases. We proved with simulations and real cases that it is useful to distinguish between gene sets or diseases with very similar phenotypes. Github R package URL is https://github.com/alexcis95/PhenoExam . Shiny App URL is https://alejandrocisterna.shinyapps.io/phenoexamweb/ .

Modeling multifunctionality of genes with secondary gene co-expression networks in human brain provides novel disease insights.

MOTIVATION: Co-expression networks are a powerful gene expression analysis method to study how genes co-express together in clusters with functional coherence that usually resemble specific cell type behavior for the genes involved. They can be applied to bulk-tissue gene expression profiling and assign function, and usually cell type specificity, to a high percentage of the gene pool used to construct the network. One of the limitations of this method is that each gene is predicted to play a role in a specific set of coherent functions in a single cell type (i.e. at most we get a single for each gene). We present here GMSCA (Gene Multifunctionality Secondary Co-expression Analysis), a software tool that exploits the co-expression paradigm to increase the number of functions and cell types ascribed to a gene in bulk-tissue co-expression networks. RESULTS: We applied GMSCA to 27 co-expression networks derived from bulk-tissue gene expression profiling of a variety of brain tissues. Neurons and glial cells (microglia, astrocytes and oligodendrocytes) were considered the main cell types. Applying this approach, we increase the overall number of predicted triplets by 46.73%. Moreover, GMSCA predicts that the SNCA gene, traditionally associated to work mainly in neurons, also plays a relevant function in oligodendrocytes. AVAILABILITYAND IMPLEMENTATION: The tool is available at GitHub, https://github.com/drlaguna/GMSCA as open-source software. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A systems-level analysis highlights microglial activation as a modifying factor in common epilepsies.

AIMS: The causes of distinct patterns of reduced cortical thickness in the common human epilepsies, detectable on neuroimaging and with important clinical consequences, are unknown. We investigated the underlying mechanisms of cortical thinning using a systems-level analysis. METHODS: Imaging-based cortical structural maps from a large-scale epilepsy neuroimaging study were overlaid with highly spatially resolved human brain gene expression data from the Allen Human Brain Atlas. Cell-type deconvolution, differential expression analysis and cell-type enrichment analyses were used to identify differences in cell-type distribution. These differences were followed up in post-mortem brain tissue from humans with epilepsy using Iba1 immunolabelling. Furthermore, to investigate a causal effect in cortical thinning, cell-type-specific depletion was used in a murine model of acquired epilepsy. RESULTS: We identified elevated fractions of microglia and endothelial cells in regions of reduced cortical thickness. Differentially expressed genes showed enrichment for microglial markers and, in particular, activated microglial states. Analysis of post-mortem brain tissue from humans with epilepsy confirmed excess activated microglia. In the murine model, transient depletion of activated microglia during the early phase of the disease development prevented cortical thinning and neuronal cell loss in the temporal cortex. Although the development of chronic seizures was unaffected, the epileptic mice with early depletion of activated microglia did not develop deficits in a non-spatial memory test seen in epileptic mice not depleted of microglia. CONCLUSIONS: These convergent data strongly implicate activated microglia in cortical thinning, representing a new dimension for concern and disease modification in the epilepsies, potentially distinct from seizure control.

Heritability Enrichment Implicates Microglia in Parkinson's Disease Pathogenesis.

OBJECTIVE: Understanding how different parts of the immune system contribute to pathogenesis in Parkinson's disease is a burning challenge with important therapeutic implications. We studied enrichment of common variant heritability for Parkinson's disease stratified by immune and brain cell types. METHODS: We used summary statistics from the most recent meta-analysis of genomewide association studies in Parkinson's disease and partitioned heritability using linkage disequilibrium score regression, stratified for specific cell types, as defined by open chromatin regions. We also validated enrichment results using a polygenic risk score approach and intersected disease-associated variants with epigenetic data and expression quantitative loci to nominate and explore a putative microglial locus. RESULTS: We found significant enrichment of Parkinson's disease risk heritability in open chromatin regions of microglia and monocytes. Genomic annotations overlapped substantially between these 2 cell types, and only the enrichment signal for microglia remained significant in a joint model. We present evidence suggesting P2RY12, a key microglial gene and target for the antithrombotic agent clopidogrel, as the likely driver of a significant Parkinson's disease association signal on chromosome 3. INTERPRETATION: Our results provide further support for the importance of immune mechanisms in Parkinson's disease pathogenesis, highlight microglial dysregulation as a contributing etiological factor, and nominate a targetable microglial gene candidate as a pathogenic player. Immune processes can be modulated by therapy, with potentially important clinical implications for future treatment in Parkinson's disease. ANN NEUROL 2021;89:942-951.

Regional brain iron and gene expression provide insights into neurodegeneration in Parkinson's disease.

The mechanisms responsible for the selective vulnerability of specific neuronal populations in Parkinson's disease are poorly understood. Oxidative stress secondary to brain iron accumulation is one postulated mechanism. We measured iron deposition in 180 cortical regions of 96 patients with Parkinson's disease and 35 control subjects using quantitative susceptibility mapping. We estimated the expression of 15 745 genes in the same regions using transcriptomic data from the Allen Human Brain Atlas. Using partial least squares regression, we then identified the profile of gene transcription in the healthy brain that underlies increased cortical iron in patients with Parkinson's disease relative to controls. Applying gene ontological tools, we investigated the biological processes and cell types associated with this transcriptomic profile and identified the sets of genes with spatial expression profiles in control brains that correlated significantly with the spatial pattern of cortical iron deposition in Parkinson's disease. Gene ontological analyses revealed that these genes were enriched for biological processes relating to heavy metal detoxification, synaptic function and nervous system development and were predominantly expressed in astrocytes and glutamatergic neurons. Furthermore, we demonstrated that the genes differentially expressed in Parkinson's disease are associated with the pattern of cortical expression identified in this study. Our findings provide mechanistic insights into regional selective vulnerabilities in Parkinson's disease, particularly the processes involving iron accumulation.

A genetic link between risk for Alzheimer's disease and severe COVID-19 outcomes via the OAS1 gene.

Recently, we reported oligoadenylate synthetase 1 (OAS1) contributed to the risk of Alzheimer's disease, by its enrichment in transcriptional networks expressed by microglia. However, the function of OAS1 within microglia was not known. Using genotyping from 1313 individuals with sporadic Alzheimer's disease and 1234 control individuals, we confirm the OAS1 variant, rs1131454, is associated with increased risk for Alzheimer's disease. The same OAS1 locus has been recently associated with severe coronavirus disease 2019 (COVID-19) outcomes, linking risk for both diseases. The single nucleotide polymorphisms rs1131454(A) and rs4766676(T) are associated with Alzheimer's disease, and rs10735079(A) and rs6489867(T) are associated with severe COVID-19, where the risk alleles are linked with decreased OAS1 expression. Analysing single-cell RNA-sequencing data of myeloid cells from Alzheimer's disease and COVID-19 patients, we identify co-expression networks containing interferon (IFN)-responsive genes, including OAS1, which are significantly upregulated with age and both diseases. In human induced pluripotent stem cell-derived microglia with lowered OAS1 expression, we show exaggerated production of TNF-α with IFN-γ stimulation, indicating OAS1 is required to limit the pro-inflammatory response of myeloid cells. Collectively, our data support a link between genetic risk for Alzheimer's disease and susceptibility to critical illness with COVID-19 centred on OAS1, a finding with potential implications for future treatments of Alzheimer's disease and COVID-19, and development of biomarkers to track disease progression.