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.

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.

Bi-allelic variants in HOPS complex subunit VPS41 cause cerebellar ataxia and abnormal membrane trafficking.

Membrane trafficking is a complex, essential process in eukaryotic cells responsible for protein transport and processing. Deficiencies in vacuolar protein sorting (VPS) proteins, key regulators of trafficking, cause abnormal intracellular segregation of macromolecules and organelles and are linked to human disease. VPS proteins function as part of complexes such as the homotypic fusion and vacuole protein sorting (HOPS) tethering complex, composed of VPS11, VPS16, VPS18, VPS33A, VPS39 and VPS41. The HOPS-specific subunit VPS41 has been reported to promote viability of dopaminergic neurons in Parkinson's disease but to date has not been linked to human disease. Here, we describe five unrelated families with nine affected individuals, all carrying homozygous variants in VPS41 that we show impact protein function. All affected individuals presented with a progressive neurodevelopmental disorder consisting of cognitive impairment, cerebellar atrophy/hypoplasia, motor dysfunction with ataxia and dystonia, and nystagmus. Zebrafish disease modelling supports the involvement of VPS41 dysfunction in the disorder, indicating lysosomal dysregulation throughout the brain and providing support for cerebellar and microglial abnormalities when vps41 was mutated. This provides the first example of human disease linked to the HOPS-specific subunit VPS41 and suggests the importance of HOPS complex activity for cerebellar function.

Recursive splicing in long vertebrate genes.

It is generally believed that splicing removes introns as single units from precursor messenger RNA transcripts. However, some long Drosophila melanogaster introns contain a cryptic site, known as a recursive splice site (RS-site), that enables a multi-step process of intron removal termed recursive splicing. The extent to which recursive splicing occurs in other species and its mechanistic basis have not been examined. Here we identify highly conserved RS-sites in genes expressed in the mammalian brain that encode proteins functioning in neuronal development. Moreover, the RS-sites are found in some of the longest introns across vertebrates. We find that vertebrate recursive splicing requires initial definition of an 'RS-exon' that follows the RS-site. The RS-exon is then excluded from the dominant mRNA isoform owing to competition with a reconstituted 5' splice site formed at the RS-site after the first splicing step. Conversely, the RS-exon is included when preceded by cryptic promoters or exons that fail to reconstitute an efficient 5' splice site. Most RS-exons contain a premature stop codon such that their inclusion can decrease mRNA stability. Thus, by establishing a binary splicing switch, RS-sites demarcate different mRNA isoforms emerging from long genes by coupling cryptic elements with inclusion of RS-exons.

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.

Mutations in HPCA cause autosomal-recessive primary isolated dystonia.

Reports of primary isolated dystonia inherited in an autosomal-recessive (AR) manner, often lumped together as "DYT2 dystonia," have appeared in the scientific literature for several decades, but no genetic cause has been identified to date. Using a combination of homozygosity mapping and whole-exome sequencing in a consanguineous kindred affected by AR isolated dystonia, we identified homozygous mutations in HPCA, a gene encoding a neuronal calcium sensor protein found almost exclusively in the brain and at particularly high levels in the striatum, as the cause of disease in this family. Subsequently, compound-heterozygous mutations in HPCA were also identified in a second independent kindred affected by AR isolated dystonia. Functional studies suggest that hippocalcin might play a role in regulating voltage-dependent calcium channels. The identification of mutations in HPCA as a cause of AR primary isolated dystonia paves the way for further studies to assess whether "DYT2 dystonia" is a genetically homogeneous condition or not.

A missense mutation in KCTD17 causes autosomal dominant myoclonus-dystonia.

Myoclonus-dystonia (M-D) is a rare movement disorder characterized by a combination of non-epileptic myoclonic jerks and dystonia. SGCE mutations represent a major cause for familial M-D being responsible for 30%-50% of cases. After excluding SGCE mutations, we identified through a combination of linkage analysis and whole-exome sequencing KCTD17 c.434 G>A p.(Arg145His) as the only segregating variant in a dominant British pedigree with seven subjects affected by M-D. A subsequent screening in a cohort of M-D cases without mutations in SGCE revealed the same KCTD17 variant in a German family. The clinical presentation of the KCTD17-mutated cases was distinct from the phenotype usually observed in M-D due to SGCE mutations. All cases initially presented with mild myoclonus affecting the upper limbs. Dystonia showed a progressive course, with increasing severity of symptoms and spreading from the cranio-cervical region to other sites. KCTD17 is abundantly expressed in all brain regions with the highest expression in the putamen. Weighted gene co-expression network analysis, based on mRNA expression profile of brain samples from neuropathologically healthy individuals, showed that KCTD17 is part of a putamen gene network, which is significantly enriched for dystonia genes. Functional annotation of the network showed an over-representation of genes involved in post-synaptic dopaminergic transmission. Functional studies in mutation bearing fibroblasts demonstrated abnormalities in endoplasmic reticulum-dependent calcium signaling. In conclusion, we demonstrate that the KCTD17 c.434 G>A p.(Arg145His) mutation causes autosomal dominant M-D. Further functional studies are warranted to further characterize the nature of KCTD17 contribution to the molecular pathogenesis of M-D.

Insights into the Influence of Specific Splicing Events on the Structural Organization of LRRK2.

Leucine-rich repeat kinase 2 (LRRK2) is a large protein of unclear function. Rare mutations in the LRRK2 gene cause familial Parkinson's disease (PD) and inflammatory bowel disease. Genome-wide association studies (GWAS) have revealed significant association of the abovementioned diseases at the LRRK2 locus. Cell and systems biology research has led to potential roles that LRRK2 may have in PD pathogenesis, especially the kinase domain (KIN). Previous human expression studies showed evidence of mRNA expression and splicing patterns that may contribute to our understanding of the function of LRRK2. In this work, we investigate and identified significant regional differences in LRRK2 expression at the mRNA level, including a number of splicing events in the Ras of complex protein (Roc) and C-terminal of Roc domain (COR) of LRRK2, in the substantia nigra (SN) and occipital cortex (OCTX). Our findings indicate that the predominant form of LRRK2 mRNA is full length, with shorter isoforms present at a lower copy number. Our molecular modelling study suggests that splicing events in the ROC/COR domains will have major consequences on the enzymatic function and dimer formation of LRRK2. The implications of these are highly relevant to the broader effort to understand the biology and physiological functions of LRRK2, and to better characterize the role(s) of LRRK2 in the underlying mechanism leading to PD.

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.

Mutations in ANO3 cause dominant craniocervical dystonia: ion channel implicated in pathogenesis.

In this study, we combined linkage analysis with whole-exome sequencing of two individuals to identify candidate causal variants in a moderately-sized UK kindred exhibiting autosomal-dominant inheritance of craniocervical dystonia. Subsequent screening of these candidate causal variants in a large number of familial and sporadic cases of cervical dystonia led to the identification of a total of six putatively pathogenic mutations in ANO3, a gene encoding a predicted Ca(2+)-gated chloride channel that we show to be highly expressed in the striatum. Functional studies using Ca(2+) imaging in case and control fibroblasts demonstrated clear abnormalities in endoplasmic-reticulum-dependent Ca(2+) signaling. We conclude that mutations in ANO3 are a cause of autosomal-dominant craniocervical dystonia. The locus DYT23 has been reserved as a synonym for this gene. The implication of an ion channel in the pathogenesis of dystonia provides insights into an alternative mechanism that opens fresh avenues for further research.

Analysis of gene expression data using a linear mixed model/finite mixture model approach: application to regional differences in the human brain.

MOTIVATION: Gene expression data exhibit common information over the genome. This article shows how data can be analysed from an efficient whole-genome perspective. Further, the methods have been developed so that users with limited expertise in bioinformatics and statistical computing techniques could use and modify this procedure to their own needs. The method outlined first uses a large-scale linear mixed model for the expression data genome-wide, and then uses finite mixture models to separate differentially expressed (DE) from non-DE transcripts. These methods are illustrated through application to an exceptional UK Brain Expression Consortium involving 12 human frozen post-mortem brain regions. RESULTS: Fitting linear mixed models has allowed variation in gene expression between different biological states (e.g. brain regions, gender, age) to be investigated. The model can be extended to allow for differing levels of variation between different biological states. Predicted values of the random effects show the effects of each transcript in a particular biological state. Using the UK Brain Expression Consortium data, this approach yielded striking patterns of co-regional gene expression. Fitting the finite mixture model to the effects within each state provides a convenient method to filter transcripts that are DE: these DE transcripts can then be extracted for advanced functional analysis. AVAILABILITY: The data for all regions except HYPO and SPCO are available at the Gene Expression Omnibus (GEO) site, accession number GSE46706. R code for the analysis is available in the Supplementary file.

Resolving the polymorphism-in-probe problem is critical for correct interpretation of expression QTL studies.

Polymorphisms in the target mRNA sequence can greatly affect the binding affinity of microarray probe sequences, leading to false-positive and false-negative expression quantitative trait locus (QTL) signals with any other polymorphisms in linkage disequilibrium. We provide the most complete solution to this problem, by using the latest genome and exome sequence reference data to identify almost all common polymorphisms (frequency >1% in Europeans) in probe sequences for two commonly used microarray panels (the gene-based Illumina Human HT12 array, which uses 50-mer probes, and exon-based Affymetrix Human Exon 1.0 ST array, which uses 25-mer probes). We demonstrate the impact of this problem using cerebellum and frontal cortex tissues from 438 neuropathologically normal individuals. We find that although only a small proportion of the probes contain polymorphisms, they account for a large proportion of apparent expression QTL signals, and therefore result in many false signals being declared as real. We find that the polymorphism-in-probe problem is insufficiently controlled by previous protocols, and illustrate this using some notable false-positive and false-negative examples in MAPT and PRICKLE1 that can be found in many eQTL databases. We recommend that both new and existing eQTL data sets should be carefully checked in order to adequately address this issue.

MAPT expression and splicing is differentially regulated by brain region: relation to genotype and implication for tauopathies.

The MAPT (microtubule-associated protein tau) locus is one of the most remarkable in neurogenetics due not only to its involvement in multiple neurodegenerative disorders, including progressive supranuclear palsy, corticobasal degeneration, Parksinson's disease and possibly Alzheimer's disease, but also due its genetic evolution and complex alternative splicing features which are, to some extent, linked and so all the more intriguing. Therefore, obtaining robust information regarding the expression, splicing and genetic regulation of this gene within the human brain is of immense importance. In this study, we used 2011 brain samples originating from 439 individuals to provide the most reliable and coherent information on the regional expression, splicing and regulation of MAPT available to date. We found significant regional variation in mRNA expression and splicing of MAPT within the human brain. Furthermore, at the gene level, the regional distribution of mRNA expression and total tau protein expression levels were largely in agreement, appearing to be highly correlated. Finally and most importantly, we show that while the reported H1/H2 association with gene level expression is likely to be due to a technical artefact, this polymorphism is associated with the expression of exon 3-containing isoforms in human brain. These findings would suggest that contrary to the prevailing view, genetic risk factors for neurodegenerative diseases at the MAPT locus are likely to operate by changing mRNA splicing in different brain regions, as opposed to the overall expression of the MAPT gene.

Mutations in the autoregulatory domain of ß-tubulin 4a cause hereditary dystonia.

Dystonia type 4 (DYT4) was first described in a large family from Heacham in Norfolk with an autosomal dominantly inherited whispering dysphonia, generalized dystonia, and a characteristic hobby horse ataxic gait. We carried out a genetic linkage analysis in the extended DYT4 family that spanned 7 generations from England and Australia, revealing a single LOD score peak of 6.33 on chromosome 19p13.12-13. Exome sequencing in 2 cousins identified a single cosegregating mutation (p.R2G) in the ß-tubulin 4a (TUBB4a) gene that was absent in a large number of controls. The mutation is highly conserved in the ß-tubulin autoregulatory MREI (methionine-arginine-glutamic acid-isoleucine) domain, highly expressed in the central nervous system, and extensive in vitro work has previously demonstrated that substitutions at residue 2, specifically R2G, disrupt the autoregulatory capability of the wild-type ß-tubulin peptide, affirming the role of the cytoskeleton in dystonia pathogenesis.

Novel C12orf65 mutations in patients with axonal neuropathy and optic atrophy.

OBJECTIVE: Charcot-Marie Tooth disease (CMT) forms a clinically and genetically heterogeneous group of disorders. Although a number of disease genes have been identified for CMT, the gene discovery for some complex form of CMT has lagged behind. The association of neuropathy and optic atrophy (also known as CMT type 6) has been described with autosomaldominant, recessive and X-linked modes of inheritance. Mutations in Mitofusin 2 have been found to cause dominant forms of CMT6. Phosphoribosylpyrophosphate synthetase-I mutations cause X-linked CMT6, but until now, mutations in the recessive forms of disease have never been identified. METHODS: We here describe a family with three affected individuals who inherited in an autosomal recessive fashion a childhood onset neuropathy and optic atrophy. Using homozygosity mapping in the family and exome sequencing in two affected individuals we identified a novel protein-truncating mutation in the C12orf65 gene, which encodes for a protein involved in mitochondrial translation. Using a variety of methods we investigated the possibility of mitochondrial impairment in the patients cell lines. RESULTS: We described a large consanguineous family with neuropathy and optic atrophy carrying a loss of function mutation in the C12orf65 gene. We report mitochondrial impairment in the patients cell lines, followed by multiple lines of evidence which include decrease of complex V activity and stability (blue native gel assay), decrease in mitochondrial respiration rate and reduction of mitochondrial membrane potential. CONCLUSIONS: This work describes a mutation in the C12orf65 gene that causes recessive form of CMT6 and confirms the role of mitochondrial dysfunction in this complex axonal neuropathy.

Genomewide association study in cervical dystonia demonstrates possible association with sodium leak channel.

Dystonia is a common movement disorder. A number of monogenic causes have been identified. However, the majority of dystonia cases are not explained by single gene defects. Cervical dystonia is one of the commonest forms without genetic causes identified. This pilot study aimed to identify large effect-size risk loci in cervical dystonia. A genomewide association study (GWAS) was performed. British resident cervical dystonia patients of European descent were genotyped using the Illumina-610-Quad. Comparison was made with controls of European descent from the Wellcome Trust Case Control Consortium using logistic regression algorithm from PLINK. SNPs not genotyped by the array were imputed with 1000 Genomes Project data using the MaCH algorithm and minimac. Postimputation analysis was done with the mach2dat algorithm using a logistic regression model. After quality control measures, 212 cases were compared with 5173 controls. No single SNP passed the genomewide significant level of 5 × 10(-8) in the analysis of genotyped SNP in PLINK. Postimputation, there were 5 clusters of SNPs that had P value <5 × 10(-6) , and the best cluster of SNPs was found near exon 1 of NALCN, (sodium leak channel) with P = 9.76 × 10(-7) . Several potential regions were found in the GWAS and imputation analysis. The lowest P value was found in NALCN. Dysfunction of this ion channel is a plausible cause for dystonia. Further replication in another cohort is needed to confirm this finding. We make this data publicly available to encourage further analyses of this disorder.

Epilepsy, hippocampal sclerosis and febrile seizures linked by common genetic variation around SCN1A.

Epilepsy comprises several syndromes, amongst the most common being mesial temporal lobe epilepsy with hippocampal sclerosis. Seizures in mesial temporal lobe epilepsy with hippocampal sclerosis are typically drug-resistant, and mesial temporal lobe epilepsy with hippocampal sclerosis is frequently associated with important co-morbidities, mandating the search for better understanding and treatment. The cause of mesial temporal lobe epilepsy with hippocampal sclerosis is unknown, but there is an association with childhood febrile seizures. Several rarer epilepsies featuring febrile seizures are caused by mutations in SCN1A, which encodes a brain-expressed sodium channel subunit targeted by many anti-epileptic drugs. We undertook a genome-wide association study in 1018 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 7552 control subjects, with validation in an independent sample set comprising 959 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 3591 control subjects. To dissect out variants related to a history of febrile seizures, we tested cases with mesial temporal lobe epilepsy with hippocampal sclerosis with (overall n = 757) and without (overall n = 803) a history of febrile seizures. Meta-analysis revealed a genome-wide significant association for mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures at the sodium channel gene cluster on chromosome 2q24.3 [rs7587026, within an intron of the SCN1A gene, P = 3.36 × 10(-9), odds ratio (A) = 1.42, 95% confidence interval: 1.26-1.59]. In a cohort of 172 individuals with febrile seizures, who did not develop epilepsy during prospective follow-up to age 13 years, and 6456 controls, no association was found for rs7587026 and febrile seizures. These findings suggest SCN1A involvement in a common epilepsy syndrome, give new direction to biological understanding of mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures, and open avenues for investigation of prognostic factors and possible prevention of epilepsy in some children with febrile seizures.

Initial assessment of the pathogenic mechanisms of the recently identified Alzheimer risk Loci.

Recent genome wide association studies have identified CLU, CR1, ABCA7 BIN1, PICALM and MS4A6A/MS4A6E in addition to the long established APOE, as loci for Alzheimer's disease. We have systematically examined each of these loci to assess whether common coding variability contributes to the risk of disease. We have also assessed the regional expression of all the genes in the brain and whether there is evidence of an eQTL explaining the risk. In agreement with other studies we find that coding variability may explain the ABCA7 association, but common coding variability does not explain any of the other loci. We were not able to show that any of the loci had eQTLs within the power of this study. Furthermore the regional expression of each of the loci did not match the pattern of brain regional distribution in Alzheimer pathology. Although these results are mainly negative, they allow us to start defining more realistic alternative approaches to determine the role of all the genetic loci involved in Alzheimer's disease.

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/.