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.

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.

Imaging genomics discovery of a new risk variant for Alzheimer's disease in the postsynaptic SHARPIN gene.

Molecular mechanisms underlying Alzheimer's disease (AD) are difficult to investigate, partly because diagnosis lags behind the insidious pathological processes. Therefore, identifying AD neuroimaging markers and their genetic modifiers may help study early mechanisms of neurodegeneration. We aimed to identify brain regions of the highest vulnerability to AD using a data-driven search in the AD Neuroimaging Initiative (ADNI, n = 1,100 subjects), and further explored genetic variants affecting this critical brain trait using both ADNI and the younger UK Biobank cohort (n = 8,428 subjects). Tensor-Based Morphometry (TBM) and Independent Component Analysis (ICA) identified the limbic system and its interconnecting white-matter as the most AD-vulnerable brain feature. Whole-genome analysis revealed a common variant in SHARPIN that was associated with this imaging feature (rs34173062, p = 2.1 × 10-10 ). This genetic association was validated in the UK Biobank, where it was correlated with entorhinal cortical thickness bilaterally (p = .002 left and p = 8.6 × 10-4 right), and with parental history of AD (p = 2.3 × 10-6 ). Our findings suggest that neuroanatomical variation in the limbic system and AD risk are associated with a novel variant in SHARPIN. The role of this postsynaptic density gene product in ß1-integrin adhesion is in line with the amyloid precursor protein (APP) intracellular signaling pathway and the recent genome-wide evidence.

Transcriptomic and genetic analyses reveal potential causal drivers for intractable partial epilepsy.

Mesial temporal lobe epilepsy with hippocampal sclerosis represents the most common epilepsy syndrome in adult patients with medically intractable partial epilepsy. Mesial temporal lobe epilepsy is usually regarded as a polygenic and complex disorder, still poorly understood but probably caused and perpetuated by dysregulation of numerous biological networks and cellular functions. The study of gene expression changes by single nucleotide polymorphisms in regulatory elements (expression quantitative trait loci, eQTLs) has been shown to be a powerful complementary approach to the detection and understanding of risk loci by genome-wide association studies. We performed a whole (gene and exon-level) transcriptome analysis on cortical tissue samples (Brodmann areas 20 and 21) from 86 patients with mesial temporal lobe epilepsy with hippocampal sclerosis and 75 neurologically healthy controls. Genome-wide genotyping data from the same individuals (patients and controls) were analysed and paired with the transcriptome data. We report potential epilepsy-risk eQTLs, some of which are specific to tissue from patients with mesial temporal lobe epilepsy with hippocampal sclerosis. We also found large transcriptional and splicing deregulation in mesial temporal lobe epilepsy with hippocampal sclerosis tissue as well as gene networks involving neuronal and glial mechanisms that provide new insights into the cause and maintenance of the seizures. These data (available via the 'Seizubraineac' web-tool resource, www.seizubraineac.org) will facilitate the identification of new therapeutic targets and biomarkers as well as genetic risk variants that could influence epilepsy and pharmacoresistance.

Increased brain expression of GPNMB is associated with genome wide significant risk for Parkinson's disease on chromosome 7p15.3.

Genome wide association studies (GWAS) for Parkinson's disease (PD) have previously revealed a significant association with a locus on chromosome 7p15.3, initially designated as the glycoprotein non-metastatic melanoma protein B (GPNMB) locus. In this study, the functional consequences of this association on expression were explored in depth by integrating different expression quantitative trait locus (eQTL) datasets (Braineac, CAGEseq, GTEx, and Phenotype-Genotype Integrator (PheGenI)). Top risk SNP rs199347 eQTLs demonstrated increased expressions of GPNMB, KLHL7, and NUPL2 with the major allele (AA) in brain, with most significant eQTLs in cortical regions, followed by putamen. In addition, decreased expression of the antisense RNA KLHL7-AS1 was observed in GTEx. Furthermore, rs199347 is an eQTL with long non-coding RNA (AC005082.12) in human tissues other than brain. Interestingly, transcript-specific eQTLs in immune-related tissues (spleen and lymphoblastoid cells) for NUPL2 and KLHL7-AS1 were observed, which suggests a complex functional role of this eQTL in specific tissues, cell types at specific time points. Significantly increased expression of GPNMB linked to rs199347 was consistent across all datasets, and taken in combination with the risk SNP being located within the GPNMB gene, these results suggest that increased expression of GPNMB is the causative link explaining the association of this locus with PD. However, other transcript eQTLs and subsequent functional roles cannot be excluded. This highlights the importance of further investigations to understand the functional interactions between the coding genes, antisense, and non-coding RNA species considering the tissue and cell-type specificity to understand the underlying biological mechanisms in PD.

Genetic architecture of sporadic frontotemporal dementia and overlap with Alzheimer's and Parkinson's diseases.

BACKGROUND: Clinical, pathological and genetic overlap between sporadic frontotemporal dementia (FTD), Alzheimer's disease (AD) and Parkinson's disease (PD) has been suggested; however, the relationship between these disorders is still not well understood. Here we evaluated genetic overlap between FTD, AD and PD to assess shared pathobiology and identify novel genetic variants associated with increased risk for FTD. METHODS: Summary statistics were obtained from the International FTD Genomics Consortium, International PD Genetics Consortium and International Genomics of AD Project (n>75 000 cases and controls). We used conjunction false discovery rate (FDR) to evaluate genetic pleiotropy and conditional FDR to identify novel FTD-associated SNPs. Relevant variants were further evaluated for expression quantitative loci. RESULTS: We observed SNPs within the HLA, MAPT and APOE regions jointly contributing to increased risk for FTD and AD or PD. By conditioning on polymorphisms associated with PD and AD, we found 11 loci associated with increased risk for FTD. Meta-analysis across two independent FTD cohorts revealed a genome-wide signal within the APOE region (rs6857, 3'-UTR=PVRL2, p=2.21×10-12), and a suggestive signal for rs1358071 within the MAPT region (intronic=CRHR1, p=4.91×10-7) with the effect allele tagging the H1 haplotype. Pleiotropic SNPs at the HLA and MAPT loci associated with expression changes in cis-genes supporting involvement of intracellular vesicular trafficking, immune response and endo/lysosomal processes. CONCLUSIONS: Our findings demonstrate genetic pleiotropy in these neurodegenerative diseases and indicate that sporadic FTD is a polygenic disorder where multiple pleiotropic loci with small effects contribute to increased disease risk.

Unravelling cell type-specific responses to Parkinson's Disease at single cell resolution.

Parkinson's Disease (PD) is the second most common neurodegenerative disorder. The pathological hallmark of PD is loss of dopaminergic neurons and the presence of aggregated α-synuclein, primarily in the substantia nigra pars compacta (SNpc) of the midbrain. However, the molecular mechanisms that underlie the pathology in different cell types is not currently understood. Here, we present a single nucleus transcriptome analysis of human post-mortem SNpc obtained from 15 sporadic Parkinson's Disease (PD) cases and 14 Controls. Our dataset comprises ∼84K nuclei, representing all major cell types of the brain, allowing us to obtain a transcriptome-level characterization of these cell types. Importantly, we identify multiple subpopulations for each cell type and describe specific gene sets that provide insights into the differing roles of these subpopulations. Our findings reveal a significant decrease in neuronal cells in PD samples, accompanied by an increase in glial cells and T cells. Subpopulation analyses demonstrate a significant depletion of tyrosine hydroxylase (TH) enriched astrocyte, microglia and oligodendrocyte populations in PD samples, as well as TH enriched neurons, which are also depleted. Moreover, marker gene analysis of the depleted subpopulations identified 28 overlapping genes, including those associated with dopamine metabolism (e.g., ALDH1A1, SLC6A3 & SLC18A2). Overall, our study provides a valuable resource for understanding the molecular mechanisms involved in dopaminergic neuron degeneration and glial responses in PD, highlighting the existence of novel subpopulations and cell type-specific gene sets.

Analysis of subcellular RNA fractions demonstrates significant genetic regulation of gene expression in human brain post-transcriptionally.

Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

Systematic visualisation of molecular QTLs reveals variant mechanisms at GWAS loci.

Splicing quantitative trait loci (QTLs) have been implicated as a common mechanism underlying complex trait associations. However, utilising splicing QTLs in target discovery and prioritisation has been challenging due to extensive data normalisation which often renders the direction of the genetic effect as well as its magnitude difficult to interpret. This is further complicated by the fact that strong expression QTLs often manifest as weak splicing QTLs and vice versa, making it difficult to uniquely identify the underlying molecular mechanism at each locus. We find that these ambiguities can be mitigated by visualising the association between the genotype and average RNA sequencing read coverage in the region. Here, we generate these QTL coverage plots for 1.7 million molecular QTL associations in the eQTL Catalogue identified with five quantification methods. We illustrate the utility of these QTL coverage plots by performing colocalisation between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. We find that while visually confirmed 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. All our association summary statistics and QTL coverage plots are freely available at https://www.ebi.ac.uk/eqtl/.

Leveraging omic features with F3UTER enables identification of unannotated 3'UTRs for synaptic genes.

There is growing evidence for the importance of 3' untranslated region (3'UTR) dependent regulatory processes. However, our current human 3'UTR catalogue is incomplete. Here, we develop a machine learning-based framework, leveraging both genomic and tissue-specific transcriptomic features to predict previously unannotated 3'UTRs. We identify unannotated 3'UTRs associated with 1,563 genes across 39 human tissues, with the greatest abundance found in the brain. These unannotated 3'UTRs are significantly enriched for RNA binding protein (RBP) motifs and exhibit high human lineage-specificity. We find that brain-specific unannotated 3'UTRs are enriched for the binding motifs of important neuronal RBPs such as TARDBP and RBFOX1, and their associated genes are involved in synaptic function. Our data is shared through an online resource F3UTER ( https://astx.shinyapps.io/F3UTER/ ). Overall, our data improves 3'UTR annotation and provides additional insights into the mRNA-RBP interactome in the human brain, with implications for our understanding of neurological and neurodevelopmental diseases.

Identification of Candidate Parkinson Disease Genes by Integrating Genome-Wide Association Study, Expression, and Epigenetic Data Sets.

Importance: Substantial genome-wide association study (GWAS) work in Parkinson disease (PD) has led to the discovery of an increasing number of loci shown reliably to be associated with increased risk of disease. Improved understanding of the underlying genes and mechanisms at these loci will be key to understanding the pathogenesis of PD. Objective: To investigate what genes and genomic processes underlie the risk of sporadic PD. Design and Setting: This genetic association study used the bioinformatic tools Coloc and transcriptome-wide association study (TWAS) to integrate PD case-control GWAS data published in 2017 with expression data (from Braineac, the Genotype-Tissue Expression [GTEx], and CommonMind) and methylation data (derived from UK Parkinson brain samples) to uncover putative gene expression and splicing mechanisms associated with PD GWAS signals. Candidate genes were further characterized using cell-type specificity, weighted gene coexpression networks, and weighted protein-protein interaction networks. Main Outcomes and Measures: It was hypothesized a priori that some genes underlying PD loci would alter PD risk through changes to expression, splicing, or methylation. Candidate genes are presented whose change in expression, splicing, or methylation are associated with risk of PD as well as the functional pathways and cell types in which these genes have an important role. Results: Gene-level analysis of expression revealed 5 genes (WDR6 [OMIM 606031], CD38 [OMIM 107270], GPNMB [OMIM 604368], RAB29 [OMIM 603949], and TMEM163 [OMIM 618978]) that replicated using both Coloc and TWAS analyses in both the GTEx and Braineac expression data sets. A further 6 genes (ZRANB3 [OMIM 615655], PCGF3 [OMIM 617543], NEK1 [OMIM 604588], NUPL2 [NCBI 11097], GALC [OMIM 606890], and CTSB [OMIM 116810]) showed evidence of disease-associated splicing effects. Cell-type specificity analysis revealed that gene expression was overall more prevalent in glial cell types compared with neurons. The weighted gene coexpression performed on the GTEx data set showed that NUPL2 is a key gene in 3 modules implicated in catabolic processes associated with protein ubiquitination and in the ubiquitin-dependent protein catabolic process in the nucleus accumbens, caudate, and putamen. TMEM163 and ZRANB3 were both important in modules in the frontal cortex and caudate, respectively, indicating regulation of signaling and cell communication. Protein interactor analysis and simulations using random networks demonstrated that the candidate genes interact significantly more with known mendelian PD and parkinsonism proteins than would be expected by chance. Conclusions and Relevance: Together, these results suggest that several candidate genes and pathways are associated with the findings observed in PD GWAS studies.

Incomplete annotation has a disproportionate impact on our understanding of Mendelian and complex neurogenetic disorders.

Growing evidence suggests that human gene annotation remains incomplete; however, it is unclear how this affects different tissues and our understanding of different disorders. Here, we detect previously unannotated transcription from Genotype-Tissue Expression RNA sequencing data across 41 human tissues. We connect this unannotated transcription to known genes, confirming that human gene annotation remains incomplete, even among well-studied genes including 63% of the Online Mendelian Inheritance in Man-morbid catalog and 317 neurodegeneration-associated genes. We find the greatest abundance of unannotated transcription in brain and genes highly expressed in brain are more likely to be reannotated. We explore examples of reannotated disease genes, such as SNCA, for which we experimentally validate a previously unidentified, brain-specific, potentially protein-coding exon. We release all tissue-specific transcriptomes through vizER: http://rytenlab.com/browser/app/vizER We anticipate that this resource will facilitate more accurate genetic analysis, with the greatest impact on our understanding of Mendelian and complex neurogenetic disorders.

Regulatory sites for splicing in human basal ganglia are enriched for disease-relevant information.

Genome-wide association studies have generated an increasing number of common genetic variants associated with neurological and psychiatric disease risk. An improved understanding of the genetic control of gene expression in human brain is vital considering this is the likely modus operandum for many causal variants. However, human brain sampling complexities limit the explanatory power of brain-related expression quantitative trait loci (eQTL) and allele-specific expression (ASE) signals. We address this, using paired genomic and transcriptomic data from putamen and substantia nigra from 117 human brains, interrogating regulation at different RNA processing stages and uncovering novel transcripts. We identify disease-relevant regulatory loci, find that splicing eQTLs are enriched for regulatory information of neuron-specific genes, that ASEs provide cell-specific regulatory information with evidence for cellular specificity, and that incomplete annotation of the brain transcriptome limits interpretation of risk loci for neuropsychiatric disease. This resource of regulatory data is accessible through our web server, http://braineacv2.inf.um.es/.

Genetic variability in response to amyloid beta deposition influences Alzheimer's disease risk.

Genome-wide association studies of late-onset Alzheimer's disease risk have previously identified genes primarily expressed in microglia that form a transcriptional network. Using transgenic mouse models of amyloid deposition, we previously showed that many of the mouse orthologues of these risk genes are co-expressed and associated with amyloid pathology. In this new study, we generate an improved RNA-seq-derived network that is expressed in amyloid-responsive mouse microglia and we statistically compare this with gene-level variation in previous human Alzheimer's disease genome-wide association studies to predict at least four new risk genes for the disease (OAS1, LAPTM5, ITGAM/CD11b and LILRB4). Of the mouse orthologues of these genes Oas1a is likely to respond directly to amyloid at the transcriptional level, similarly to established risk gene Trem2, because the increase in Oas1a and Trem2 transcripts in response to amyloid deposition in transgenic mice is significantly higher than both the increase of the average microglial transcript and the increase in microglial number. In contrast, the mouse orthologues of LAPTM5, ITGAM/CD11b and LILRB4 (Laptm5, Itgam/CD11b and Lilra5) show increased transcripts in the presence of amyloid plaques similar in magnitude to the increase of the average microglial transcript and the increase in microglia number, except that Laptm5 and Lilra5 transcripts increase significantly quicker than the average microglial transcript as the plaque load becomes dense. This work suggests that genetic variability in the microglial response to amyloid deposition is a major determinant for Alzheimer's disease risk, and identification of these genes may help to predict the risk of developing Alzheimer's disease. These findings also provide further insights into the mechanisms underlying Alzheimer's disease for potential drug discovery.

An additional k-means clustering step improves the biological features of WGCNA gene co-expression networks.

BACKGROUND: Weighted Gene Co-expression Network Analysis (WGCNA) is a widely used R software package for the generation of gene co-expression networks (GCN). WGCNA generates both a GCN and a derived partitioning of clusters of genes (modules). We propose k-means clustering as an additional processing step to conventional WGCNA, which we have implemented in the R package km2gcn (k-means to gene co-expression network, https://github.com/juanbot/km2gcn ). RESULTS: We assessed our method on networks created from UKBEC data (10 different human brain tissues), on networks created from GTEx data (42 human tissues, including 13 brain tissues), and on simulated networks derived from GTEx data. We observed substantially improved module properties, including: (1) few or zero misplaced genes; (2) increased counts of replicable clusters in alternate tissues (x3.1 on average); (3) improved enrichment of Gene Ontology terms (seen in 48/52 GCNs) (4) improved cell type enrichment signals (seen in 21/23 brain GCNs); and (5) more accurate partitions in simulated data according to a range of similarity indices. CONCLUSIONS: The results obtained from our investigations indicate that our k-means method, applied as an adjunct to standard WGCNA, results in better network partitions. These improved partitions enable more fruitful downstream analyses, as gene modules are more biologically meaningful.

Frontotemporal dementia: insights into the biological underpinnings of disease through gene co-expression network analysis.

BACKGROUND: In frontotemporal dementia (FTD) there is a critical lack in the understanding of biological and molecular mechanisms involved in disease pathogenesis. The heterogeneous genetic features associated with FTD suggest that multiple disease-mechanisms are likely to contribute to the development of this neurodegenerative condition. We here present a systems biology approach with the scope of i) shedding light on the biological processes potentially implicated in the pathogenesis of FTD and ii) identifying novel potential risk factors for FTD. We performed a gene co-expression network analysis of microarray expression data from 101 individuals without neurodegenerative diseases to explore regional-specific co-expression patterns in the frontal and temporal cortices for 12 genes (MAPT, GRN, CHMP2B, CTSC, HLA-DRA, TMEM106B, C9orf72, VCP, UBQLN2, OPTN, TARDBP and FUS) associated with FTD and we then carried out gene set enrichment and pathway analyses, and investigated known protein-protein interactors (PPIs) of FTD-genes products. RESULTS: Gene co-expression networks revealed that several FTD-genes (such as MAPT and GRN, CTSC and HLA-DRA, TMEM106B, and C9orf72, VCP, UBQLN2 and OPTN) were clustering in modules of relevance in the frontal and temporal cortices. Functional annotation and pathway analyses of such modules indicated enrichment for: i) DNA metabolism, i.e. transcription regulation, DNA protection and chromatin remodelling (MAPT and GRN modules); ii) immune and lysosomal processes (CTSC and HLA-DRA modules), and; iii) protein meta/catabolism (C9orf72, VCP, UBQLN2 and OPTN, and TMEM106B modules). PPI analysis supported the results of the functional annotation and pathway analyses. CONCLUSIONS: This work further characterizes known FTD-genes and elaborates on their biological relevance to disease: not only do we indicate likely impacted regional-specific biological processes driven by FTD-genes containing modules, but also do we suggest novel potential risk factors among the FTD-genes interactors as targets for further mechanistic characterization in hypothesis driven cell biology work.

A genome-wide gene-expression analysis and database in transgenic mice during development of amyloid or tau pathology.

We provide microarray data comparing genome-wide differential expression and pathology throughout life in four lines of "amyloid" transgenic mice (mutant human APP, PSEN1, or APP/PSEN1) and "TAU" transgenic mice (mutant human MAPT gene). Microarray data were validated by qPCR and by comparison to human studies, including genome-wide association study (GWAS) hits. Immune gene expression correlated tightly with plaques whereas synaptic genes correlated negatively with neurofibrillary tangles. Network analysis of immune gene modules revealed six hub genes in hippocampus of amyloid mice, four in common with cortex. The hippocampal network in TAU mice was similar except that Trem2 had hub status only in amyloid mice. The cortical network of TAU mice was entirely different with more hub genes and few in common with the other networks, suggesting reasons for specificity of cortical dysfunction in FTDP17. This Resource opens up many areas for investigation. All data are available and searchable at http://www.mouseac.org.

Genetic variability in the regulation of gene expression in ten regions of the human brain.

Germ-line genetic control of gene expression occurs via expression quantitative trait loci (eQTLs). We present a large, exon-specific eQTL data set covering ten human brain regions. We found that cis-eQTL signals (within 1 Mb of their target gene) were numerous, and many acted heterogeneously among regions and exons. Co-regulation analysis of shared eQTL signals produced well-defined modules of region-specific co-regulated genes, in contrast to standard coexpression analysis of the same samples. We report cis-eQTL signals for 23.1% of catalogued genome-wide association study hits for adult-onset neurological disorders. The data set is publicly available via public data repositories and via http://www.braineac.org/. Our study increases our understanding of the regulation of gene expression in the human brain and will be of value to others pursuing functional follow-up of disease-associated variants.