Search

Abstract Mutations in GBA1 encoding the lysosomal enzyme β-glucocerebrosidase (GCase) are among the most prevalent genetic susceptibility factors for Parkinson’s disease (PD), with 10–30% of carriers developing the disease. To identify genetic modifiers contributing to the incomplete penetrance, we examined the effect of 1634 human transcription factors (TFs) on GCase activity in lysates of an engineered human glioblastoma line homozygous for the pathogenic GBA1 L444P variant. Using an arrayed CRISPR activation library, we uncovered 11 TFs as regulators of GCase activity. Among these, activation of MITF and TFEC increased lysosomal GCase activity in live cells, while activation of ONECUT2 and USF2 decreased it. While MITF, TFEC, and USF2 affected GBA1 transcription, ONECUT2 might control GCase trafficking. The effects of MITF, TFEC, and USF2 on lysosomal GCase activity were reproducible in iPSC-derived neurons from PD patients. Our study provides a systematic approach to identifying modulators of GCase activity and deepens our understanding of the mechanisms regulating GCase.

Abstract Understanding the molecular mechanisms underlying frontotemporal dementia (FTD) is essential for the development of successful therapies. Systematic studies on human post-mortem brain tissue of patients with genetic subtypes of FTD are currently lacking. The Risk and Modyfing Factors of Frontotemporal Dementia (RiMod-FTD) consortium therefore has generated a multi-omics dataset for genetic subtypes of FTD to identify common and distinct molecular mechanisms disturbed in disease. Here, we present multi-omics datasets generated from the frontal lobe of post-mortem human brain tissue from patients with mutations in MAPT, GRN and C9orf72 and healthy controls. This data resource consists of four datasets generated with different technologies to capture the transcriptome by RNA-seq, small RNA-seq, CAGE-seq, and methylation profiling. We show concrete examples on how to use the resulting data and confirm current knowledge about FTD and identify new processes for further investigation. This extensive multi-omics dataset holds great value to reveal new research avenues for this devastating disease.

Abstract The Monogenic Network of the Global Parkinson’s Genetics Program (GP2) aims to create an efficient infrastructure to accelerate the identification of novel genetic causes of Parkinson’s disease (PD) and to improve our understanding of already identified genetic causes, such as reduced penetrance and variable clinical expressivity of known disease-causing variants. We aim to perform short- and long-read whole-genome sequencing for up to 10,000 patients with parkinsonism. Important features of this project are global involvement and focusing on historically underrepresented populations.

Abstract Frontotemporal dementia is characterized by progressive atrophy of frontal and/or temporal cortices at an early age of onset. The disorder shows considerable clinical, pathological, and genetic heterogeneity. Here we investigated the proteomic signatures of frontal and temporal cortex from brains with frontotemporal dementia due to GRN and MAPT mutations to identify the key cell types and molecular pathways in their pathophysiology. We compared patients with mutations in the GRN gene (n = 9) or with mutations in the MAPT gene (n = 13) with non-demented controls (n = 11). Using quantitative proteomic analysis on laser-dissected tissues we identified brain region-specific protein signatures for both genetic subtypes. Using published single cell RNA expression data resources we deduced the involvement of major brain cell types in driving these different protein signatures. Subsequent gene ontology analysis identified distinct genetic subtype- and cell type-specific biological processes. For the GRN subtype, we observed a distinct role for immune processes related to endothelial cells and for mitochondrial dysregulation in neurons. For the MAPT subtype, we observed distinct involvement of dysregulated RNA processing, oligodendrocyte dysfunction, and axonal impairments. Comparison with an in-house protein signature of Alzheimer’s disease brains indicated that the observed alterations in RNA processing and oligodendrocyte function are distinct for the frontotemporal dementia MAPT subtype. Taken together, our results indicate the involvement of different brain cell types and biological mechanisms in genetic subtypes of frontotemporal dementia. Furthermore, we demonstrate that comparison of proteomic profiles of different disease entities can separate general neurodegenerative processes from disease-specific pathways, which may aid the development of disease subtype-specific treatment strategies.

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

Summary: TSPEAR variants cause autosomal recessive ectodermal dysplasia (ARED) 14. The function of TSPEAR is unknown. The clinical features, the mutation spectrum, and the underlying mechanisms of ARED14 are poorly understood. Combining data from new and previously published individuals established that ARED14 is primarily characterized by dental anomalies such as conical tooth cusps and hypodontia, like those seen in individuals with WNT10A-related odontoonychodermal dysplasia. AlphaFold-predicted structure-based analysis showed that most of the pathogenic TSPEAR missense variants likely destabilize the β-propeller of the protein. Analysis of 100000 Genomes Project (100KGP) data revealed multiple founder TSPEAR variants across different populations. Mutational and recombination clock analyses demonstrated that non-Finnish European founder variants likely originated around the end of the last ice age, a period of major climatic transition. Analysis of gnomAD data showed that the non-Finnish European population TSPEAR gene-carrier rate is ∼1/140, making it one of the commonest AREDs. Phylogenetic and AlphaFold structural analyses showed that TSPEAR is an ortholog of drosophila Closca, an extracellular matrix-dependent signaling regulator. We, therefore, hypothesized that TSPEAR could have a role in enamel knot, a structure that coordinates patterning of developing tooth cusps. Analysis of mouse single-cell RNA sequencing (scRNA-seq) data revealed highly restricted expression of Tspear in clusters representing enamel knots. A tspeara−/−;tspearb−/− double-knockout zebrafish model recapitulated the clinical features of ARED14 and fin regeneration abnormalities of wnt10a knockout fish, thus suggesting interaction between tspear and wnt10a. In summary, we provide insights into the role of TSPEAR in ectodermal development and the evolutionary history, epidemiology, mechanisms, and consequences of its loss of function variants.

The transcriptional activity of Transposable Elements (TEs) has been involved in numerous pathological processes, including neurodegenerative diseases such as amyotrophic lateral sclerosis and frontotemporal lobar degeneration. The TE expression analysis from short-read sequencing technologies is, however, challenging due to the multitude of similar sequences derived from singular TEs subfamilies and the exaptation of TEs within longer coding or non-coding RNAs. Specialised tools have been developed to quantify the expression of TEs that either relies on probabilistic re-distribution of multimapper count fractions or allow for discarding multimappers altogether. Until now, the benchmarking across those tools was largely limited to aggregated expression estimates over whole TEs subfamilies. Here, we compared the performance of recently published tools (SQuIRE, TElocal, SalmonTE) with simplistic quantification strategies (featureCounts in unique, fraction and random modes) at the individual loci level. Using simulated datasets, we examined the false discovery rate and the primary driver of those false positive hits in the optimal quantification strategy. Our findings suggest a high false discovery number that exceeds the total number of correctly recovered active loci for all the quantification strategies, including the best performing tool TElocal. As a remedy, filtering based on the minimum number of read counts or baseMean expression improves the F1 score and decreases the number of false positives. Finally, we demonstrate that additional profiling of Transcription Start Site mapping statistics (using a k-means clustering approach) significantly improves the performance of TElocal while reporting a reliable set of detected and differentially expressed TEs in human simulated RNA-seq data.

Frontotemporal lobar degeneration (FTLD) includes a heterogeneous group of disorders pathologically characterized by the degeneration of the frontal and temporal lobes. In addition to major genetic contributors of FTLD such as mutations in MAPT, GRN, and C9orf72, recent work has identified several epigenetic modifications including significant differential DNA methylation in DLX1, and OTUD4 loci. As aging remains one of the major risk factors for FTLD, we investigated the presence of accelerated epigenetic aging in FTLD compared to controls. We calculated epigenetic age in both peripheral blood and brain tissues of multiple FTLD subtypes using several DNA methylation clocks, i.e., DNAmClockMulti, DNAmClockHannum, DNAmClockCortical, GrimAge, and PhenoAge, and determined age acceleration and its association with different cellular proportions and clinical traits. Significant epigenetic age acceleration was observed in the peripheral blood of both frontotemporal dementia (FTD) and progressive supranuclear palsy (PSP) patients compared to controls with DNAmClockHannum, even after accounting for confounding factors. A similar trend was observed with both DNAmClockMulti and DNAmClockCortical in post-mortem frontal cortex tissue of PSP patients and in FTLD cases harboring GRN mutations. Our findings support that increased epigenetic age acceleration in the peripheral blood could be an indicator for PSP and to a smaller extent, FTD.

Abstract Patient-derived cellular models become an increasingly powerful tool to model human diseases for precision medicine approaches. The identification of robust cellular disease phenotypes in these models paved the way towards high throughput screenings (HTS) including the implementation of laboratory advanced automation. However, maintenance and expansion of cells for HTS remains largely manual work. Here, we describe an integrated, complex automated platform for HTS in a translational research setting also designed for maintenance and expansion of different cell types. The comprehensive design allows automation of all cultivation steps and is flexible for development of methods for variable cell types. We demonstrate protocols for controlled cell seeding, splitting and expansion of human fibroblasts, induced pluripotent stem cells (iPSC), and neural progenitor cells (NPC) that allow for subsequent differentiation into different cell types and image-based multiparametric screening. Furthermore, we provide automated protocols for neuronal differentiation of NPC in 2D culture and 3D midbrain organoids for HTS. The flexibility of this multitask platform makes it an ideal solution for translational research settings involving experiments on different patient-derived cellular models for precision medicine.

Machado-Joseph disease (MJD/SCA3) is a neurodegenerative polyglutamine disorder exhibiting a wide spectrum of phenotypes. The abnormal size of the (CAG)n at ATXN3 explains ~55% of the age at onset variance, suggesting the involvement of other factors, namely genetic modifiers, whose identification remains limited. Our aim was to find novel genetic modifiers, analyse their epistatic effects and identify disease-modifying pathways contributing to MJD variable expressivity. We performed whole-exome sequencing in a discovery sample of four age at onset concordant and four discordant first-degree relative pairs of Azorean patients, to identify candidate variants which genotypes differed for each discordant pair but were shared in each concordant pair. Variants identified by this approach were then tested in an independent multi-origin cohort of 282 MJD patients. Whole-exome sequencing identified 233 candidate variants, from which 82 variants in 53 genes were prioritized for downstream analysis. Eighteen disease-modifying pathways were identified; two of the most enriched pathways were relevant for the nervous system, namely the neuregulin signaling and the agrin interactions at neuromuscular junction. Variants at PARD3, NFKB1, CHD5, ACTG1, CFAP57, DLGAP2, ITGB1, DIDO1 and CERS4 modulate age at onset in MJD, with those identified in CFAP57, ACTG1 and DIDO1 showing consistent effects across cohorts of different geographical origins. Network analyses of the nine novel MJD modifiers highlighted several important molecular interactions, including genes/proteins previously related with MJD pathogenesis, namely between ACTG1/APOE and VCP/ITGB1. We describe novel pathways, modifiers, and their interaction partners, providing a broad molecular portrait of age at onset modulation to be further exploited as new disease-modifying targets for MJD and related diseases.

The cellular prion protein PrPC is necessary for prion replication, and its reduction greatly increases life expectancy in animal models of prion infection. Hence the factors controlling the levels of PrPC may represent therapeutic targets against human prion diseases. Here we performed an arrayed whole-transcriptome RNA interference screen to identify modulators of PrPC expression. We cultured human U251-MG glioblastoma cells in the presence of 64’752 unique siRNAs targeting 21’584 annotated human genes, and measured PrPC using a one-pot fluorescence-resonance energy transfer immunoassay in 51’128 individual microplate wells. This screen yielded 743 candidate regulators of PrPC. When downregulated, 563 of these candidates reduced and 180 enhanced PrPC expression. Recursive candidate attrition through multiple secondary screens yielded 54 novel regulators of PrPC, 9 of which were confirmed by CRISPR interference as robust regulators of PrPC biosynthesis and degradation. The phenotypes of 6 of the 9 candidates were inverted in response to transcriptional activation using CRISPRa. The RNA-binding post-transcriptional repressor Pumilio-1 was identified as a potent limiter of PrPC expression through the degradation of PRNP mRNA. Because of its hypothesis-free design, this comprehensive genetic-perturbation screen delivers an unbiased landscape of the genes regulating PrPC levels in cells, most of which were unanticipated, and some of which may be amenable to pharmacological targeting in the context of antiprion therapies. Author summary The cellular prion protein (PrPC) acts as both, the substrate for prion formation and mediator of prion toxicity during the progression of all prion diseases. Suppressing the levels of PrPC is a viable therapeutic strategy as PRNP null animals are resistant to prion disease and the knockout of PRNP is not associated with any severe phenotypes. Motivated by the scarcity of knowledge regarding the molecular regulators of PrPC biosynthesis and degradation, which might serve as valuable targets to control its expression, here, we present a cell-based genome wide RNAi screen in arrayed format. The screening effort led to the identification of 54 regulators, nine of which were confirmed by an independent CRISPR-based method. Among the final nine targets, we identified PUM1 as a regulator of PRNP mRNA by acting on the 3’UTR promoting its degradation. The newly identified factors involved in the life cycle of PrPC provided by our study may also represent themselves as therapeutic targets for the intervention of prion diseases.

The H1 haplotype of the microtubule-associated protein tau (MAPT) gene is a common genetic risk factor for some neurodegenerative diseases such as progressive supranuclear palsy, corticobasal degeneration, and Parkinson’s disease. The molecular mechanism causing the increased risk for the named diseases, however, remains unclear. In this paper, we present a valuable tool of eight small molecule neural precursor cell lines (smNPC) homozygous for the MAPT haplotypes (four H1/H1 and four H2/H2 cell lines), which can be used to identify MAPT-dependent phenotypes. The employed differentiation protocol is fast due to overexpression of NEUROGENIN-2 and therefore suitable for high-throughput approaches. A basic characterization of all human cell lines was performed, and their TAU and α-SYNUCLEIN profiles were compared during a differentiation time of 30 days. We could identify higher levels of conformationally altered TAU in cell lines carrying the H2 haplotype. Additionally, we found increased expression levels of α-SYNUCLEIN in H1/H1 cells. With this resource, we aim to fill a gap in neurodegenerative disease modeling with induced pluripotent stem cells (iPSC) for sporadic tauopathies.

Genetic variants underlying complex traits, including disease susceptibility, are enriched within the transcriptional regulatory elements, promoters and enhancers. There is emerging evidence that regulatory elements associated with particular traits or diseases share similar patterns of transcriptional activity. Accordingly, shared transcriptional activity (coexpression) may help prioritise loci associated with a given trait, and help to identify underlying biological processes. Using cap analysis of gene expression (CAGE) profiles of promoter- and enhancer-derived RNAs across 1824 human samples, we have analysed coexpression of RNAs originating from trait-associated regulatory regions using a novel quantitative method (network density analysis; NDA). For most traits studied, phenotype-associated variants in regulatory regions were linked to tightly-coexpressed networks that are likely to share important functional characteristics. Coexpression provides a new signal, independent of phenotype association, to enable fine mapping of causative variants. The NDA coexpression approach identifies new genetic variants associated with specific traits, including an association between the regulation of the OCT1 cation transporter and genetic variants underlying circulating cholesterol levels. NDA strongly implicates particular cell types and tissues in disease pathogenesis. For example, distinct groupings of disease-associated regulatory regions implicate two distinct biological processes in the pathogenesis of ulcerative colitis; a further two separate processes are implicated in Crohn's disease. Thus, our functional analysis of genetic predisposition to disease defines new distinct disease endotypes. We predict that patients with a preponderance of susceptibility variants in each group are likely to respond differently to pharmacological therapy. Together, these findings enable a deeper biological understanding of the causal basis of complex traits.

Background: PLA2G6-associated neurodegeneration (PLAN) is a recessive neurodegenerative disorder characterized by three distinct phenotypes: infantile neuroaxonal dystrophy (INAD), atypical neuroaxonal dystrophy (atypical NAD), and PLA2G6-related dystonia–parkinsonism. Methods: A consanguineous index case from Turkey was diagnosed with early-onset Parkinsonism at the Istanbul Faculty of Medicine. She and her unaffected brother were subjected to whole-genome sequencing. Results: In this report, we describe a 33-year-old index case with parental consanguinity and early-onset Parkinsonism. Whole-genome sequencing of this individual revealed that a homozygous p.R747W mutation in PLA2G6 segregates with the disease in this family Discussion: This result supports the importance of prioritizing this gene in mutational analysis of autosomal recessive Parkinsonism, and confirms the clinical heterogeneity of PLAN.

'Microtubule-associated protein tau' (MAPT), 'granulin' (GRN) and 'chromosome 9 open reading frame72' (C9ORF72) gene mutations are the major known genetic causes of frontotemporal dementia (FTD). Recent studies suggest that mutations in these genes may also be associated with other forms of dementia. Therefore we investigated whether MAPT, GRN and C9ORF72 gene mutations are major contributors to dementia in a random, unselected Turkish cohort of dementia patients. A combination of whole-exome sequencing, Sanger sequencing and fragment analysis/Southern blot was performed in order to identify pathogenic mutations and novel variants in these genes as well as other FTD-related genes such as the 'charged multivesicular body protein 2B' (CHMP2B), the 'FUS RNA binding protein' (FUS), the 'TAR DNA binding protein' (TARDBP), the 'sequestosome1' (SQSTM1), and the 'valosin containing protein' (VCP). We determined one pathogenic MAPT mutation (c.1906C>T, p.P636L) and one novel missense variant (c.38A>G, p.D13G). In GRN we identified a probably pathogenic TGAG deletion in the splice donor site of exon 6. Three patients were found to carry the GGGGCC expansions in the non-coding region of the C9ORF72 gene. In summary, a complete screening for mutations in MAPT, GRN and C9ORF72 genes revealed a frequency of 5.4% of pathogenic mutations in a random cohort of 93 Turkish index patients with dementia.

Subjects with incidental Lewy body disease (iLBD) may represent the premotor stage of Parkinson's disease (PD). To elucidate molecular mechanisms underlying neuronal dysfunction and alpha-synuclein pathology in the premotor phase of PD, we investigated the transcriptome of the substantia nigra (SN) of well-characterized iLBD, PD donors and age-matched controls with Braak alpha-synuclein stage ranging from 0-6. In Braak alpha-synuclein stages 1 and 2, we observed deregulation of pathways linked to axonal degeneration, immune response and endocytosis, including axonal guidance signaling, mTOR signaling, EIF2 signaling and clathrin-mediated endocytosis in the SN. In Braak stages 3 and 4, we observed deregulation of pathways involved in protein translation and cell survival, including mTOR and EIF2 signaling. In Braak stages 5 and 6, we observed deregulation of dopaminergic signaling, axonal guidance signaling and thrombin signaling. Throughout the progression of PD pathology, we observed a deregulation of mTOR, EIF2 and regulation of eIF4 and p70S6K signaling in the SN. Our results indicate that molecular mechanisms related to axonal dysfunction, endocytosis and immune response are an early event in PD pathology, whereas mTOR and EIF2 signaling are impaired throughout disease progression. These pathways may hold the key to altering the disease progression in PD.

A Swiss frontotemporal dementia (FTD) kindred with extrapyramidal-like features and without motor neuron disease shows a brain pathology with ubiquitin-positive but tau-negative inclusions. Tau and neurofilament modifications are now studied here in three recently deceased family members. No major and specific decrease of tau was observed as described by others in, e.g., sporadic cases of FTD with absence of tau-positive inclusions. However, a slight decrease of tau, neurofilament, and synaptic proteins, resulting from frontal atrophy was detected. In parallel, polymorphic markers on chromosome 17q21–22, the centromeric region of chromosome 3 and chromosome 9, were tested. Haplotype analysis showed several recombination events for chromosomes 3 and 17, but patients shared a haplotype on chromosome 9q21–22. However as one of the patients exhibited Alzheimer and vascular dementia pathology with uncertain concomitant FTD, this locus is questionable. Altogether, these data indicate principally that the Swiss kindred is unlinked to locus 17q21–22, and that tau is not at the origin of FTD in this family.

In this work we investigate the role of CHIP in a new CHIP-mutation related ataxia and the therapeutic potential of trehalose. The patient's fibroblasts with a new form of hereditary ataxia, related to STUB1 gene (CHIP) mutations, and three age and sex-matched controls were treated with epoxomicin and trehalose. The effects on cell death, protein misfolding and proteostasis were evaluated. Recent studies have revealed that mutations in STUB-1 gene lead to a growing list of molecular defects as deregulation of protein quality, inhibition of proteasome, cell death, decreased autophagy and alteration in CHIP and HSP70 levels. In this CHIP-mutant patient fibroblasts the inhibition of proteasome with epoxomicin induced severe pathophysiological age-associated changes, cell death and protein ubiquitination. Additionally, treatment with epoxomicin produced a dose-dependent increase in the number of cleaved caspase-3 positive cells. However, co-treatment with trehalose, a disaccharide of glucose present in a wide variety of organisms and known as a autophagy enhancer, reduced these pathological events. Trehalose application also increased CHIP and HSP70 expression and GSH free radical levels. Furthermore, trehalose augmented macro and chaperone mediated autophagy (CMA), rising the levels of LC3, LAMP2, CD63 and increasing the expression of Beclin-1 and Atg5-Atg12. Trehalose treatment in addition increased the percentage of immunoreactive cells to HSC70 and LAMP2 and reduced the autophagic substrate, p62. Although this is an individual case based on only one patient and the statistical comparisons are not valid between controls and patient, the low variability among controls and the obvious differences with this patient allow us to conclude that trehalose, through its autophagy activation capacity, anti-aggregation properties, anti-oxidative effects and lack of toxicity, could be very promising for the treatment of CHIP-mutation related ataxia, and possibly a wide spectrum of neurodegenerative disorders related to protein disconformation.

The direct estimation of heritability from genome-wide common variant data as implemented in the program Genome-wide Complex Trait Analysis (GCTA) has provided a means to quantify heritability attributable to all interrogated variants. We have quantified the variance in liability to disease explained by all SNPs for two phenotypically-related neurobehavioral disorders, obsessive-compulsive disorder (OCD) and Tourette Syndrome (TS), using GCTA. Our analysis yielded a heritability point estimate of 0.58 (se = 0.09, p = 5.64e-12) for TS, and 0.37 (se = 0.07, p = 1.5e-07) for OCD. In addition, we conducted multiple genomic partitioning analyses to identify genomic elements that concentrate this heritability. We examined genomic architectures of TS and OCD by chromosome, MAF bin, and functional annotations. In addition, we assessed heritability for early onset and adult onset OCD. Among other notable results, we found that SNPs with a minor allele frequency of less than 5% accounted for 21% of the TS heritability and 0% of the OCD heritability. Additionally, we identified a significant contribution to TS and OCD heritability by variants significantly associated with gene expression in two regions of the brain (parietal cortex and cerebellum) for which we had available expression quantitative trait loci (eQTLs). Finally we analyzed the genetic correlation between TS and OCD, revealing a genetic correlation of 0.41 (se = 0.15, p = 0.002). These results are very close to previous heritability estimates for TS and OCD based on twin and family studies, suggesting that very little, if any, heritability is truly missing (i.e., unassayed) from TS and OCD GWAS studies of common variation. The results also indicate that there is some genetic overlap between these two phenotypically-related neuropsychiatric disorders, but suggest that the two disorders have distinct genetic architectures.

Association studies have identified several signals at the LRRK2 locus for Parkinson's disease (PD), Crohn's disease (CD) and leprosy. However, little is known about the molecular mechanisms mediating these effects. To further characterize this locus, we fine-mapped the risk association in 5,802 PD and 5,556 controls using a dense genotyping array (ImmunoChip). Using samples from 134 post-mortem control adult human brains (UK Human Brain Expression Consortium), where up to ten brain regions were available per individual, we studied the regional variation, splicing and regulation of LRRK2. We found convincing evidence for a common variant PD association located outside of the LRRK2 protein coding region (rs117762348, A>G, P = 2.56×10(-8), case/control MAF 0.083/0.074, odds ratio 0.86 for the minor allele with 95% confidence interval [0.80-0.91]). We show that mRNA expression levels are highest in cortical regions and lowest in cerebellum. We find an exon quantitative trait locus (QTL) in brain samples that localizes to exons 32-33 and investigate the molecular basis of this eQTL using RNA-Seq data in n = 8 brain samples. The genotype underlying this eQTL is in strong linkage disequilibrium with the CD associated non-synonymous SNP rs3761863 (M2397T). We found two additional QTLs in liver and monocyte samples but none of these explained the common variant PD association at rs117762348. Our results characterize the LRRK2 locus, and highlight the importance and difficulties of fine-mapping and integration of multiple datasets to delineate pathogenic variants and thus develop an understanding of disease mechanisms.

Multiple lines of evidence suggest that fatty acids (FA) play an important role in cognitive function. However, little is known about the functional genetic pathways involved in cognition. The main goals of this study were to replicate previously reported interaction effects between breast feeding (BF) and FA desaturase (FADS) genetic variation on IQ and to investigate the possible mechanisms by which these variants might moderate BF effect, focusing on brain expression. Using a sample of 534 twins, we observed a trend in the moderation of BF effects on IQ by FADS2 variation. In addition, we made use of publicly available gene expression databases from both humans (193) and mice (93) and showed that FADS2 variants also correlate with FADS1 brain expression (P-value

Major depressive disorder (MDD) is a psychiatric disorder, characterized by periods of low mood of more than two weeks, loss of interest in normally enjoyable activities and behavioral changes. MDD is a complex disorder and does not have a single genetic cause. In 2009 a genome wide association study (GWAS) was performed on the Dutch GAIN-MDD cohort. Many of the top signals of this GWAS mapped to a region spanning the gene PCLO, and the non-synonymous coding single nucleotide polymorphism (SNP) rs2522833 in the PCLO gene became genome wide significant after post-hoc analysis. We performed resequencing of PCLO, GRM7, and SLC6A4 in 50 control samples from the GAIN-MDD cohort, to detect new genomic variants. Subsequently, we genotyped these variants in the entire GAIN-MDD cohort and performed association analysis to investigate if rs2522833 is the causal variant or simply in linkage disequilibrium with a more associated variant. GRM7 and SLC6A4 are both candidate genes for MDD from literature. We aimed to gather more evidence that rs2522833 is indeed the causal variant in the GAIN-MDD cohort or to find a previously undetected common variant in either PCLO, GRM7, or SLC6A4 with a higher association in this cohort. After next generation sequencing and association analysis we excluded the possibility of an undetected common variant to be more associated. For neither PCLO nor GRM7 we found a more associated variant. For SLC6A4, we found a new SNP that showed a lower P-value (P = 0.07) than in the GAIN-MDD GWAS (P = 0.09). However, no evidence for genome-wide significance was found. Although we did not take into account rare variants, we conclude that our results provide further support for the hypothesis that the non-synonymous coding SNP rs2522833 in the PCLO gene is indeed likely to be the causal variant in the GAIN-MDD cohort.

More than 800 published genetic association studies have implicated dozens of potential risk loci in Parkinson's disease (PD). To facilitate the interpretation of these findings, we have created a dedicated online resource, PDGene, that comprehensively collects and meta-analyzes all published studies in the field. A systematic literature screen of -27,000 articles yielded 828 eligible articles from which relevant data were extracted. In addition, individual-level data from three publicly available genome-wide association studies (GWAS) were obtained and subjected to genotype imputation and analysis. Overall, we performed meta-analyses on more than seven million polymorphisms originating either from GWAS datasets and/or from smaller scale PD association studies. Meta-analyses on 147 SNPs were supplemented by unpublished GWAS data from up to 16,452 PD cases and 48,810 controls. Eleven loci showed genome-wide significant (P < 5 × 10(-8)) association with disease risk: BST1, CCDC62/HIP1R, DGKQ/GAK, GBA, LRRK2, MAPT, MCCC1/LAMP3, PARK16, SNCA, STK39, and SYT11/RAB25. In addition, we identified novel evidence for genome-wide significant association with a polymorphism in ITGA8 (rs7077361, OR 0.88, P = 1.3 × 10(-8)). All meta-analysis results are freely available on a dedicated online database (www.pdgene.org), which is cross-linked with a customized track on the UCSC Genome Browser. Our study provides an exhaustive and up-to-date summary of the status of PD genetics research that can be readily scaled to include the results of future large-scale genetics projects, including next-generation sequencing studies.

Parkinson's disease (PD) occurs in both familial and sporadic forms, and both monogenic and complex genetic factors have been identified. Early onset PD (EOPD) is particularly associated with autosomal recessive (AR) mutations, and three genes, PARK2, PARK7 and PINK1, have been found to carry mutations leading to AR disease. Since mutations in these genes account for less than 10% of EOPD patients, we hypothesized that further recessive genetic factors are involved in this disorder, which may appear in extended runs of homozygosity.We carried out genome wide SNP genotyping to look for extended runs of homozygosity (ROHs) in 1,445 EOPD cases and 6,987 controls. Logistic regression analyses showed an increased level of genomic homozygosity in EOPD cases compared to controls. These differences are larger for ROH of 9 Mb and above, where there is a more than three-fold increase in the proportion of cases carrying a ROH. These differences are not explained by occult recessive mutations at existing loci. Controlling for genome wide homozygosity in logistic regression analyses increased the differences between cases and controls, indicating that in EOPD cases ROHs do not simply relate to genome wide measures of inbreeding. Homozygosity at a locus on chromosome19p13.3 was identified as being more common in EOPD cases as compared to controls. Sequencing analysis of genes and predicted transcripts within this locus failed to identify a novel mutation causing EOPD in our cohort.There is an increased rate of genome wide homozygosity in EOPD, as measured by an increase in ROHs. These ROHs are a signature of inbreeding and do not necessarily harbour disease-causing genetic variants. Although there might be other regions of interest apart from chromosome 19p13.3, we lack the power to detect them with this analysis.

Major depressive disorder (MDD) is a psychiatric disorder that is characterized--amongst others--by persistent depressed mood, loss of interest and pleasure and psychomotor retardation. Environmental circumstances have proven to influence the aetiology of the disease, but MDD also has an estimated 40% heritability, probably with a polygenic background. In 2009, a genome wide association study (GWAS) was performed on the Dutch GAIN-MDD cohort. A non-synonymous coding single nucleotide polymorphism (SNP) rs2522833 in the PCLO gene became only nominally significant after post-hoc analysis with an Australian cohort which used similar ascertainment. The absence of genome-wide significance may be caused by low SNP coverage of genes. To increase SNP coverage to 100% for common variants (m.a.f.>0.1, r(2)>0.8), we selected seven genes from the GAIN-MDD GWAS: PCLO, GZMK, ANPEP, AFAP1L1, ST3GAL6, FGF14 and PTK2B. We genotyped 349 SNPs and obtained the lowest P-value for rs2715147 in PCLO at P = 6.8E-7. We imputed, filling in missing genotypes, after which rs2715147 and rs2715148 showed the lowest P-value at P = 1.2E-6. When we created a haplotype of these SNPs together with the non-synonymous coding SNP rs2522833, the P-value decreased to P = 9.9E-7 but was not genome wide significant. Although our study did not identify a more strongly associated variant, the results for PCLO suggest that the causal variant is in high LD with rs2715147, rs2715148 and rs2522833.

BACKGROUND:Understanding the aetiologies of neurodegenerative diseases such as Alzheimer's disease (AD), Pick's disease (PiD), Progressive Supranuclear Palsy (PSP) and Frontotemporal dementia (FTD) is often hampered by the considerable clinical and molecular overlap between these diseases and normal ageing. The development of high throughput genomic technologies such as microarrays provide a new molecular tool to gain insight in the complexity and relationships between diseases, as they provide data on the simultaneous activity of multiple genes, gene networks and cellular pathways. METHODOLOGY/PRINCIPAL FINDINGS:We have constructed genome wide expression profiles from snap frozen post-mortem tissue from the medial temporal lobe of patients with four neurodegenerative disorders (5 AD, 5 PSP, 5 PiD and 5 FTD patients) and 5 control subjects. All patients were matched for age, gender, ApoE-epsilon and MAPT (tau) haplotype. From all groups a total of 790 probes were shown to be differently expressed when compared to control individuals. The results from these experiments were then used to investigate the correlations between clinical, pathological and molecular findings. From the 790 identified probes we extracted a gene set of 166 probes whose expression could discriminate between these disorders and normal ageing. CONCLUSIONS/SIGNIFICANCE:From genome wide expression profiles we extracted a gene set of 166 probes whose expression could discriminate between neurological disorders and normal ageing. This gene set can be further developed into an accurate microarray-based classification test. Furthermore, from this dataset we extracted a disease specific set of genes and identified two aging related transcription factors (FOXO1A and FOXO3A) as possible drug targets related to neurodegenerative disease.

A single nucleotide polymorphism (rs5848) located in the 3'- untranslated region of GRN has recently been associated with a risk of frontotemporal lobar degeneration (FTLD) in North American population particularly in pathologically confirmed cases with neural inclusions immunoreactive for ubiquitin and TAR DNA-binding protein 43 (TDP-43), but negative for tau and alpha-synuclein (FTLD-TDP).In an effort to replicate these results in a different population, rs5848 was genotyped in 256 FTLD cases and 1695 controls from the Netherlands. Single SNP gender-adjusted logistic regression analysis revealed no significant association between variation at rs5848 and FTLD. Fisher's exact test, failed to find any significant association between rs5848 and a subset of 23 pathology confirmed FTLD-TDP cases.The evidence presented here suggests that variation at rs5848 does not contribute to the etiology of FTLD in the Dutch population.

The human genome encodes a limited number of genes yet contributes to individual differences in a vast array of heritable traits. A possible explanation for the capacity our genome to generate this virtually unlimited range of phenotypic variation in complex traits is to assume functional interactions between genes. Therefore we searched two mammalian genomes to identify potential epistatic interactions by looking for co-adapted genes marked by excess two-locus genetic differentiation between populations/lineages using publicly available SNP genotype data. The practical motivation for this effort is to reduce the number of pair-wise tests that need to be performed in genome-wide association studies aimed at detecting GxG interactions, by focusing on pairs predicted to be more likely to jointly affect variation in complex traits. Hence, this approach generates a list of candidate interactions that can be empirically tested. In both the mouse and human data we observed two-locus genetic differentiation in excess of what can be expected from chance alone based on simulations. In an attempt to validate our hypothesis that pairs of genes showing excess genetic divergence represent potential functional interactions, we selected a small set of gene combinations postulated to be interacting based on our analyses and looked for a combined effect of the selected genes on variation in complex traits in both mice and man. In both cases the individual effect of the genes were not significant, instead we observed marginally significant interaction effects. These results show that genome wide searches for gene-gene interactions based on population genetic data are feasible and can generate interesting candidate gene pairs to be further tested for their contribution to phenotypic variation in complex traits.