Linkage disequilibrium and association of MAPT H1 in Parkinson disease

Critical reasoning on the co-expression module QTL in the dorsolateral prefrontal cortex

Expression quantitative trait locus (eQTL) analysis is a popular method of gaining insight into the function of regulatory variation. While cis-eQTL resources have been instrumental in linking genome-wide association study variants to gene function, complex trait heritability may be additionally mediated by other forms of gene regulation. Toward this end, novel eQTL methods leverage gene co-expression (module-QTL) to investigate joint regulation of gene modules by single genetic variants. Here we broadly define a "module-QTL" as the association of a genetic variant with a summary measure of gene co-expression. This approach aims to reduce the multiple testing burden of a trans-eQTL search through the consolidation of gene-based testing and provide biological context to eQTLs shared between genes. In this article we provide an in-depth examination of the co-expression module eQTL (module-QTL) through literature review, theoretical investigation, and real-data application of the module-QTL to three large prefrontal cortex genotype-RNA sequencing datasets. We find module-QTLs in our study that are disease associated and reproducible are not additionally informative beyond cis- or trans-eQTLs for module genes. Through comparison to prior studies, we highlight promises and limitations of the module-QTL across study designs and provide recommendations for further investigation of the module-QTL framework.

MAPT Locus in Parkinson's Disease Patients of Ashkenazi Origin: A Stratified Analysis

Introduction: MAPT locus is associated with Parkinson's disease (PD), which is located within a large inversion region of high linkage disequilibrium (LD). We aimed to determine whether the H2-haplotype protective effect and its effect size depends on the GBA1 or LRRK2 risk allele carrier status, and to further characterize genetic alterations that might contribute to its effect. Methods: LD analysis was performed using whole-genome sequencing data of 202 unrelated Ashkenazi Jewish (AJ) PDs. A haplotype-divergent variant was genotyped in a cohort of 1200 consecutively recruited AJ-PDs. The odd ratios were calculated using AJ-non-neuro cases from the gnomAD database as the controls in an un-stratified and a stratified manner according to the mutation carrier status, and the effect on the Age at Motor Symptom Onset (AMSO) was examined. Expression and splicing quantitative trait locus (eQTL and sQTL) analyses were carried out using brain tissues from a database. Results: The H2 haplotype exhibited significant association with PD protection, with a similar effect size in GBA1 carriers, LRRK2-G2019S carriers, and non-carriers (OR = 0.77, 0.69, and 0.82, respectively), and there was no effect on AMSO. The LD interval was narrowed to approximately 1.2 Mb. The H2 haplotype carried potential variants in candidate genes (MAPT and SPPL2C); structural deletions and segmental duplication (KANSL1); and variants affecting gene expression and intron excision ratio in brain tissues (LRRC37A/2). Conclusions: Our results demonstrate that H2 is associated with PD and its protective effect is not influenced by the GBA1/LRRK2 risk allele carrier status. This effect may be genetically complex, resulting from different levels of variations such as missense mutations in relevant genes, structural variations, epigenetic modifications, and RNA expression changes, which may operate independently or in synergy.

Genetic landscape of Parkinson's disease and related diseases in Luxembourg

Objectives

To explore the genetic architecture of PD in the Luxembourg Parkinson's Study including cohorts of healthy people and patients with Parkinson's disease (PD) and atypical parkinsonism (AP).

Methods

809 healthy controls, 680 PD and 103 AP were genotyped using the Neurochip array. We screened and validated rare single nucleotide variants (SNVs) and copy number variants (CNVs) within seven PD-causing genes (LRRK2, SNCA, VPS35, PRKN, PARK7, PINK1 and ATP13A2). Polygenic risk scores (PRSs) were generated using the latest genome-wide association study for PD. We then estimated the role of common variants in PD risk by applying gene-set-specific PRSs.

Results

We identified 60 rare SNVs in seven PD-causing genes, nine of which were pathogenic in LRRK2, PINK1 and PRKN. Eleven rare CNVs were detected in PRKN including seven duplications and four deletions. The majority of PRKN SNVs and CNVs carriers were heterozygous and not differentially distributed between cases and controls. The PRSs were significantly associated with PD and identified specific molecular pathways related to protein metabolism and signal transduction as drivers of PD risk.

Conclusion

We performed a comprehensive genetic characterization of the deep-phenotyped individuals of the Luxembourgish Parkinson's Study. Heterozygous SNVs and CNVs in PRKN were not associated with higher PD risk. In particular, we reported novel digenic variants in PD related genes and rare LRRK2 SNVs in AP patients. Our findings will help future studies to unravel the genetic complexity of PD.

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

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

Inversions on human chromosomes

Human chromosome inversions are types of balanced structural variations, making them difficult to analyze. Thanks to PEM (paired-end sequencing and mapping), there has been tremendous progress in studying inversions. Inversions play an important role as an evolutionary factor, contributing to the formation of gonosomes, speciation of chimpanzees and humans, and inv17q21.3 or inv8p23.1 exhibit the features of natural selection. Both inversions have been related to pathogenic phenotype by directly affecting a gene structure (e.g., inv5p15.1q14.1), regulating gene expression (e.g., inv7q21.3q35) and by predisposing to other secondary arrangements (e.g., inv7q11.23). A polymorphism of human inversions is documented by the InvFEST database (a database that stores information about clinical predictions, validations, frequency of inversions, etc.), but only a small fraction of these inversions is validated, and a detailed analysis is complicated by the frequent location of breakpoints within regions of repetitive sequences.

Association studies between chromosomal regions 1q21.3, 5q21.3, 14q21.2 and 17q21.31 and numbers of children in Poland

Number of children is an important human trait and studies have indicated associations with single-nucleotide polymorphisms (SNPs). Aim: to give further evidence for four associations using a large sample of Polish subjects. Data from the POPULOUS genetic database was provided from anonymous, healthy, unrelated, Polish volunteers of both sexes (N = 5760). SNPs (n = 173) studied: (a) 69 from the chromosome 17 H1/H2 inversion; (b) six from 1q21.3, 5q21.3 and 14q21.2; and (c) 98 random negative controls. Zero-inflated negative-binomial regression (z.i.) was performed (0-3 numbers of children per individual (NCI) set as non-events; adjustors: year of birth, sex). Significance level p = 0.05 with Bonferroni correction. Statistically-significant differences (with data from both sexes combined) were obtained from highly-linked inversion SNPs: representative rs12373123 gave means: homozygotes TT: 2.31 NCI (n = 1418); heterozygotes CT: 2.35 NCI (n = 554); homozygotes CC: 2.44 NCI (n = 43) (genotype p = 0.01; TTvs.CC p = 0.004; CTvs.CC p = 0.009). (Male data alone gave similar results.) Recessive modeling indicated that H2-homozygotes had 0.118 more children than H1-homozygotes + heterozygotes (z.i.-count estimates ± standard errors: CT, - 0.508 ± 0.194; TT, - 0.557 ± 0.191). The non-over-dispersed count model detected no interactions: of importance there was no significant interaction with age. No positive results were obtained from negative-control SNPs or (b). Conclusions: association between the H1/H2 inversion and numbers of children (previously reported in Iceland) has been confirmed, albeit using a different statistical model. One limitation is the small amount of data, despite initially ~ 6000 subjects. Causal studies require further investigation.

Cell-type-specific cis-eQTLs in eight human brain cell types identify novel risk genes for psychiatric and neurological disorders

To date, most expression quantitative trait loci (eQTL) studies, which investigate how genetic variants contribute to gene expression, have been performed in heterogeneous brain tissues rather than specific cell types. In this study, we performed an eQTL analysis using single-nuclei RNA sequencing from 192 individuals in eight brain cell types derived from the prefrontal cortex, temporal cortex and white matter. We identified 7,607 eGenes, a substantial fraction (46%, 3,537/7,607) of which show cell-type-specific effects, with strongest effects in microglia. Cell-type-level eQTLs affected more constrained genes and had larger effect sizes than tissue-level eQTLs. Integration of brain cell type eQTLs with genome-wide association studies (GWAS) revealed novel relationships between expression and disease risk for neuropsychiatric and neurodegenerative diseases. For most GWAS loci, a single gene co-localized in a single cell type, providing new clues into disease etiology. Our findings demonstrate substantial contrast in genetic regulation of gene expression among brain cell types and reveal potential mechanisms by which disease risk genes influence brain disorders.

Identification of recurrent genetic patterns from targeted sequencing panels with advanced data science: a case-study on sporadic and genetic neurodegenerative diseases

Background

Targeted Next Generation Sequencing is a common and powerful approach used in both clinical and research settings. However, at present, a large fraction of the acquired genetic information is not used since pathogenicity cannot be assessed for most variants. Further complicating this scenario is the increasingly frequent description of a poli/oligogenic pattern of inheritance showing the contribution of multiple variants in increasing disease risk. We present an approach in which the entire genetic information provided by target sequencing is transformed into binary data on which we performed statistical, machine learning, and network analyses to extract all valuable information from the entire genetic profile. To test this approach and unbiasedly explore the presence of recurrent genetic patterns, we studied a cohort of 112 patients affected either by genetic Creutzfeldt–Jakob (CJD) disease caused by two mutations in the PRNP gene (p.E200K and p.V210I) with different penetrance or by sporadic Alzheimer disease (sAD).

Results

Unsupervised methods can identify functionally relevant sources of variation in the data, like haplogroups and polymorphisms that do not follow Hardy–Weinberg equilibrium, such as the NOTCH3 rs11670823 (c.3837 + 21 T > A). Supervised classifiers can recognize clinical phenotypes with high accuracy based on the mutational profile of patients. In addition, we found a similar alteration of allele frequencies compared the European population in sporadic patients and in V210I-CJD, a poorly penetrant PRNP mutation, and sAD, suggesting shared oligogenic patterns in different types of dementia. Pathway enrichment and protein–protein interaction network revealed different altered pathways between the two PRNP mutations.

Conclusions

We propose this workflow as a possible approach to gain deeper insights into the genetic information derived from target sequencing, to identify recurrent genetic patterns and improve the understanding of complex diseases. This work could also represent a possible starting point of a predictive tool for personalized medicine and advanced diagnostic applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01173-4.

SARS-CoV-2 susceptibility and COVID-19 disease severity are associated with genetic variants affecting gene expression in a variety of tissues

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we identify 4 genomic loci with suggestive associations for SARS-CoV-2 susceptibility and 19 for COVID-19 disease severity. Four of these 23 loci likely have an ethnicity-specific component. Genome-wide association study (GWAS) signals in 11 loci colocalize with expression quantitative trait loci (eQTLs) associated with the expression of 20 genes in 62 tissues/cell types (range: 1:43 tissues/gene), including lung, brain, heart, muscle, and skin as well as the digestive system and immune system. We perform genetic fine mapping to compute 99% credible SNP sets, which identify 10 GWAS loci that have eight or fewer SNPs in the credible set, including three loci with one single likely causal SNP. Our study suggests that the diverse symptoms and disease severity of COVID-19 observed between individuals is associated with variants across the genome, affecting gene expression levels in a wide variety of tissue types.

Tau and MAPT genetics in tauopathies and synucleinopathies

MAPT encodes the microtubule-associated protein tau, which is the main component of neurofibrillary tangles (NFTs) and found in other protein aggregates. These aggregates are among the pathological hallmarks of primary tauopathies such as frontotemporal dementia (FTD). Abnormal tau can also be observed in secondary tauopathies such as Alzheimer's disease (AD) and synucleinopathies such as Parkinson's disease (PD). On top of pathological findings, genetic data also links MAPT to these disorders. MAPT variations are a cause or risk factors for many tauopathies and synucleinopathies and are associated with certain clinical and pathological features in affected individuals. In addition to clinical, pathological, and genetic overlap, evidence also suggests that tau and alpha-synuclein may interact on the molecular level, and thus might collaborate in the neurodegenerative process. Understanding the role of MAPT variations in tauopathies and synucleinopathies is therefore essential to elucidate the role of tau in the pathogenesis and phenotype of those disorders, and ultimately to develop targeted therapies. In this review, we describe the role of MAPT genetic variations in tauopathies and synucleinopathies, several genotype-phenotype and pathological features, and discuss their implications for the classification and treatment of those disorders.

[Reconstruction of gene network associated with Parkinson disease for gene targets search]

Accumulation of genetic data in the field of Parkinson's disease research culminated in identifying risk factors and confident prediction of the disease occurrence. To find new gene-targets for diagnostics and therapy we have to reconstruct gene network of the disease, to cluster genes in the network, to reveal key (hub) genes with largest number of interactions in the network. Using the on-line bioinformatics tools OMIM, PANTHER, g:Profiler, GeneMANIA, and STRING-DB, we have analyzed the current array of data related to Parkinson's disease, calculated the categories of gene ontologies for a large list of genes, visualized them, and built gene networks containing the identified key objects and their relationships. However, translating the results into biological understanding is still a promising major challenge. The analysis of the genes associated with the disease, the assessment of their place in the gene network (connectivity) allows us to evaluate them as target genes for medicinal effects.

An integrative systems-based analysis of substance use: eQTL-informed gene-based tests, gene networks, and biological mechanisms

Genome-wide association studies have identified multiple genetic risk factors underlying susceptibility to substance use, however, the functional genes and biological mechanisms remain poorly understood. The discovery and characterization of risk genes can be facilitated by the integration of genome-wide association data and gene expression data across biologically relevant tissues and/or cell types to identify genes whose expression is altered by DNA sequence variation (expression quantitative trait loci; eQTLs). The integration of gene expression data can be extended to the study of genetic co-expression, under the biologically valid assumption that genes form co-expression networks to influence the manifestation of a disease or trait. Here, we integrate genome-wide association data with gene expression data from 13 brain tissues to identify candidate risk genes for 8 substance use phenotypes. We then test for the enrichment of candidate risk genes within tissue-specific gene co-expression networks to identify modules (or groups) of functionally related genes whose dysregulation is associated with variation in substance use. We identified eight gene modules in brain that were enriched with gene-based association signals for substance use phenotypes. For example, a single module of 40 co-expressed genes was enriched with gene-based associations for drinks per week and biological pathways involved in GABA synthesis, release, reuptake and degradation. Our study demonstrates the utility of eQTL and gene co-expression analysis to uncover novel biological mechanisms for substance use traits.

Transcript Variants of Genes Involved in Neurodegeneration Are Differentially Regulated by the APOE and MAPT Haplotypes

Genetic variations at the Apolipoprotein E (ApoE) and microtubule-associated protein tau (MAPT) loci have been implicated in multiple neurogenerative diseases, but their exact molecular mechanisms are unclear. In this study, we performed transcript level linear modelling using the blood whole transcriptome data and genotypes of the 570 subjects in the Parkinson’s Progression Markers Initiative (PPMI) cohort. ApoE, MAPT haplotypes and two SNPs at the SNCA locus (rs356181, rs3910105) were used to detect expression quantitative trait loci eQTLs associated with the transcriptome and differential usage of transcript isoforms. As a result, we identified 151 genes associated with the genotypic variations, 29 cis and 122 trans eQTL positions. Profound effect with genome-wide significance of ApoE e4 haplotype on the expression of TOMM40 transcripts was identified. This finding potentially explains in part the frequently established genetic association with the APOE e4 haplotypes in neurodegenerative diseases. Moreover, MAPT haplotypes had significant differential impact on 23 transcripts from the 17q21.31 and 17q24.1 loci. MAPT haplotypes had also the largest up-regulating (256) and the largest down-regulating (−178) effect sizes measured as β values on two different transcripts from the same gene (LRRC37A2). Intronic SNP in the SNCA gene, rs3910105, differentially induced expression of three SNCA isoforms. In conclusion, this study established clear association between well-known haplotypic variance and transcript specific regulation in the blood. APOE e4 and MAPT H1/H2 haplotypic variants are associated with the expression of several genes related to the neurodegeneration.

LRRK2; a dynamic regulator of cellular trafficking

Parkinson's disease (PD) represents the second most common neurodegenerative disorder, characterized clinically by bradykinesia, resting tremor, rigidity and postural instability, and a variety of non-motor features. The etiology of PD is unknown, however genetic, environmental and inflammatory factors may influence disease onset and progression. Genetic variability in leucine-rich repeat kinase 2 confers significant genotypic and population-attributable risk for LRRK2-parkinsonism that is clinically indistinguishable from idiopathic PD. Nevertheless, the age-associated midbrain pathology observed post-mortem in LRRK2-parkinsonism may involve the abnormal accumulation of either α-synuclein or tau, or just the loss of dopaminergic neurons and gliosis. While diverse biological functions have been described for this multi-domain protein in many cell types, evidence suggests LRRK2 may sense endosomal trafficking to orchestrate dynamic changes in vesicular flux and cytoskeletal architecture. This review posits the long-held belief that synaptic-axonal dysfunction and terminal degeneration may precede dopaminergic cell loss, and provocatively questions how facets of LRRK2 biology may influence this molecular pathogenesis.

Screening non-MAPT genes of the Chr17q21 H1 haplotype in Parkinson's disease

INTRODUCTION

The microtubule-associated protein tau (MAPT) gene is considered a strong genetic risk factor for Parkinson's disease (PD) in Caucasians. MAPT is located within an inversion region of high linkage disequilibrium designated as H1 and H2 haplotype, and contains eight other genes which have been implicated in neurodegeneration. The aim of the current study was to identify common coding variants in strong linkage disequilibrium (LD) within the associated loci on chr17q21 harboring MAPT.

METHODS

Sanger sequencing of coding exons in 90 Caucasian late-onset PD (LOPD) patients was performed. Specific gene sequencing for LRRC37A, LRRC37A2, ARL17A and ARL17B was not possible given the high homology, presence of pseudogenes and copy number variants that are in the region, and therefore four genes (NSF, KANSL1, SPPL2C, and CRHR1) were included in the analysis. Coding variants from these four genes that did not perfectly tag (r² = 1) the MAPT H1/H2 haplotype were genotyped in an independent replication series of Caucasian PD cases (N = 851) and controls (N = 730).

RESULTS

In the 90 LOPD cases we identified 30 coding variants. Eleven non-synonymous variants tagged the MAPT H1/H2 haplotype, including two SPPL2C variants (rs12185233 and rs12373123) that had high pathogenic combined annotation dependent depletion (CADD) scores of >20. In the replication series, the non-synonymous KANSL1 rs17585974 variant was in very strong LD with MAPT H1/H2 and had a high CADD score of 24.7.

CONCLUSION

We have identified several non-synonymous variants across neighboring genes of MAPT that may warrant further genetic and functional investigation within the biological etiology of PD.

Association of MAPT subhaplotypes with clinical and demographic features in Parkinson's disease

Objective

To determine whether distinct microtubule‐associated protein tau MAPT H1 subhaplotypes are associated with clinical and demographic features in Parkinson’s disease.

Methods

A retrospective cohort study included 855 unrelated Caucasian patients with Parkinson’s disease who were seen by Movement Disorder specialists at the Mayo Clinic Florida between 1998 and 2016. The primary outcome measures were specific demographic and clinical features of Parkinson’s disease, including age at onset, disease progression, survival, motor signs, dementia, dystonia, dyskinesia, autonomic dysfunction, impulse control disorder, psychiatric features, REM sleep behavior disorder, restless legs syndrome, and Parkinson’s disease subtype. Specific clinical features were measured at the initial visit and most recent visit. These outcomes were assessed for association with MAPT H1 subhaplotypes, which were defined by six haplotype tagging variants.

Results

Median onset age was 64 years (range: 22‐94 years); 548 (64%) of patients were male. Significant associations (P < 0.0029) were observed between MAPT H1b and orthostatic hypotension (OR = 1.72, P = 0.001); between H1j and rest tremor (OR = 0.15; P < 0.001) as well as REM sleep behavior disorder (OR = 3.87, P < 0.001); between H1r and bradykinesia (OR = 0.11; P < 0.001); and between H1v and restless legs syndrome (OR = 4.02, P = 0.002).

Interpretation

Four MAPT H1 subhaplotypes, but not the H2 haplotype, were significantly associated with specific clinical features in Parkinson’s disease. MAPT haplotypic structure may explain some of the phenotypic variability in disease. Replication of our findings will be critical to fully resolve the Parkinson’s disease risk association signal at Chr17q21.

Epigenetics in Lewy Body Diseases: Impact on Gene Expression, Utility as a Biomarker, and Possibilities for Therapy

Lewy body disorders (LBD) include Parkinson's disease (PD) and dementia with Lewy bodies (DLB). They are synucleinopathies with a heterogeneous clinical manifestation. As a cause of neuropathological overlap with other neurodegenerative diseases, the establishment of a correct clinical diagnosis is still challenging, and clinical management may be difficult. The combination of genetic variation and epigenetic changes comprising gene expression-modulating DNA methylation and histone alterations modifies the phenotype, disease course, and susceptibility to disease. In this review, we summarize the results achieved in the deciphering of the LBD epigenome. To provide an appropriate context, first LBD genetics is briefly outlined. Afterwards, a detailed review of epigenetic modifications identified for LBD in human cells, postmortem, and peripheral tissues is provided. We also focus on the difficulty of identifying epigenome-related biomarker candidates and discuss the results obtained so far. Additionally, epigenetic changes as therapeutic targets, as well as different epigenome-based treatments, are revised. The number of studies focusing on PD is relatively limited and practically inexistent for DLB. There is a lack of replication studies, and some results are even contradictory, probably due to differences in sample collection and analytical techniques. In summary, we show the current achievements and directions for future research.

Early cognitive decline after bilateral subthalamic deep brain stimulation in Parkinson's disease patients with GBA mutations

BACKGROUND

Subthalamic nucleus deep brain stimulation (STN-DBS) has demonstrated its efficacy on motor complications in advanced Parkinson's disease (PD) but does not modify disease progression. Genetic forms of PD have been associated with different cognitive progression profiles.

OBJECTIVE

To assess the effect of PD-related genetic mutations on cognitive outcome after STN-DBS.

METHODS

Patients with STN-DBS were screened for LRRK2, GBA, and PRKN mutations at the Pitié-Salpêtrière Hospital between 1997 and 2009. Patients with known monogenetic forms of PD from six other centers were also included. The Mattis Dementia Rating Scale (MDRS) was used to evaluate cognition at baseline and one-year post-surgery. The standardized Unified PD Rating Scale (UPDRS) evaluation On and Off medication/DBS was also administered. A generalized linear model adjusted for sex, ethnicity, age at onset, and disease duration was used to evaluate the effect of genetic factors on MDRS changes.

RESULTS

We analyzed 208 patients (131 males, 77 females, 54.3 ± 8.8 years) including 25 GBA, 18 LRRK2, 22 PRKN, and 143 PD patients without mutations. PRKN patients were younger and had a longer disease duration at baseline. A GBA mutation was the only significant genetic factor associated with MDRS change (β = -2.51, p = 0.009). GBA mutation carriers had a more pronounced post-operative MDRS decline (3.2 ± 5.1) than patients with LRRK2 (0.9 ± 4.8), PRKN (0.5 ± 2.7) or controls (1.4 ± 4.4). The motor response to DBS was similar between groups.

CONCLUSION

GBA mutations are associated with early cognitive decline following STN-DBS. Neuropsychological assessment and discussions on the benefit/risk ratio of DBS are particularly important for this population.

Bayesian meta-analysis across genome-wide association studies of diverse phenotypes

Genome-wide association studies (GWAS) are a powerful tool for understanding the genetic basis of diseases and traits, but most studies have been conducted in isolation, with a focus on either a single or a set of closely related phenotypes. We describe MetABF, a simple Bayesian framework for performing integrative meta-analysis across multiple GWAS using summary statistics. The approach is applicable across a wide range of study designs and can increase the power by 50% compared with standard frequentist tests when only a subset of studies have a true effect. We demonstrate its utility in a meta-analysis of 20 diverse GWAS which were part of the Wellcome Trust Case Control Consortium 2. The novelty of the approach is its ability to explore, and assess the evidence for a range of possible true patterns of association across studies in a computationally efficient framework.

[The genetic predictors of cognitive impairment in Parkinson's disease]

Parkinson's disease (PD) is a common neurodegenerative disorder that can be both sporadic and familial. A number of studies are devoted to the study of non-motor symptoms in PD today. Cognitive deficits, and especially dementia, are one of the most severe and disabling non-motor symptoms of PD. More than a quarter of patients in the early stages of PD have a moderate cognitive impairment, more than half of patients with PD develop dementia within 10 years from the date of diagnosis. Using genome-wide association studies (GWAS), a number of genes associated with cognitive impairment have been identified based on a comparison of genetic and clinical phenotypes. These genes can be divided into three groups: genes that lead to the development of PD and are inherited according to the laws of Mendel (SNCA), genes that are risk factors for PD development (GBA, MAPT) and genes associated with the development of cognitive impairment, but not with PD (COMT, APOE, BDNF). This review examines the effect of genetic variants in the above-mentioned genes on cognitive functions in patients with PD. The elucidation of the genetic basis of cognitive deficits in PD could help in choice of treatment tactics and in development of new therapeutic strategies.

Advances in the discovery of genetic risk factors for complex forms of neurodegenerative disorders: contemporary approaches, success, challenges and prospects

Neurodegenerative diseases constitute a large proportion of disorders in elderly, majority being sporadic in occurrence with ∼5-10% familial. A strong genetic component underlies the Mendelian forms but nongenetic factors together with genetic vulnerability contributes to the complex sporadic forms. Several gene discoveries in the familial forms have provided novel insights into the pathogenesis of neurodegeneration with implications for treatment. Conversely, findings from genetic dissection of the sporadic forms, despite large genomewide association studies and more recently whole exome and whole genome sequencing, have been limited. This review provides a concise account of the genetics that we know, the pathways that they implicate, the challenges that are faced and the prospects that are envisaged for the sporadic, complex forms of neurodegenerative diseases, taking four most common conditions, namely Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis and Huntington disease as examples. Poor replication across studies, inability to establish genotype-phenotype correlations and the overall failure to predict risk and/or prevent disease in this group poses a continuing challenge. Among others, clinical heterogeneity emerges as the most important impediment warranting newer approaches. Advanced computational and system biology tools to analyse the big data are being generated and the alternate strategy such as subgrouping of case-control cohorts based on deep phenotyping using the principles of Ayurveda to overcome current limitation of phenotype heterogeneity seem to hold promise. However, at this point, with advances in discovery genomics and functional analysis of putative determinants with translation potential for the complex forms being minimal, stem cell therapies are being attempted as potential interventions. In this context, the possibility to generate patient derived induced pluripotent stem cells, mutant/gene/genome correction through CRISPR/Cas9 technology and repopulating the specific brain regions with corrected neurons, which may fulfil the dream of personalized medicine have been mentioned briefly. Understanding disease pathways/biology using this technology, with implications for development of novel therapeutics are optimistic expectations in the near future.

Genetic risk factors in Parkinson's disease

Over the last two decades, we have witnessed a revolution in the field of Parkinson's disease (PD) genetics. Great advances have been made in identifying many loci that confer a risk for PD, which has subsequently led to an improved understanding of the molecular pathways involved in disease pathogenesis. Despite this success, it is predicted that only a relatively small proportion of the phenotypic variability has been explained by genetics. Therefore, it is clear that common heritable components of disease are still to be identified. Dissecting the genetic architecture of PD constitutes a critical effort in identifying therapeutic targets and although such substantial progress has helped us to better understand disease mechanism, the route to PD disease-modifying drugs is a lengthy one. In this review, we give an overview of the known genetic risk factors in PD, focusing not on individual variants but the larger networks that have been implicated following comprehensive pathway analysis. We outline the challenges faced in the translation of risk loci to pathobiological relevance and illustrate the need for integrating big-data by noting success in recent work which adopts a broad-scale screening approach. Lastly, with PD genetics now progressing from identifying risk to predicting disease, we review how these models will likely have a significant impact in the future.

Inversion variants in human and primate genomes

For many years, inversions have been proposed to be a direct driving force in speciation since they suppress recombination when heterozygous. Inversions are the most common large-scale differences among humans and great apes. Nevertheless, they represent large events easily distinguishable by classical cytogenetics, whose resolution, however, is limited. Here, we performed a genome-wide comparison between human, great ape, and macaque genomes using the net alignments for the most recent releases of genome assemblies. We identified a total of 156 putative inversions, between 103 kb and 91 Mb, corresponding to 136 human loci. Combining literature, sequence, and experimental analyses, we analyzed 109 of these loci and found 67 regions inverted in one or multiple primates, including 28 newly identified inversions. These events overlap with 81 human genes at their breakpoints, and seven correspond to sites of recurrent rearrangements associated with human disease. This work doubles the number of validated primate inversions larger than 100 kb, beyond what was previously documented. We identified 74 sites of errors, where the sequence has been assembled in the wrong orientation, in the reference genomes analyzed. Our data serve two purposes: First, we generated a map of evolutionary inversions in these genomes representing a resource for interrogating differences among these species at a functional level; second, we provide a list of misassembled regions in these primate genomes, involving over 300 Mb of DNA and 1978 human genes. Accurately annotating these regions in the genome references has immediate applications for evolutionary and biomedical studies on primates.

Meta-analysis of the association between variants in MAPT and neurodegenerative diseases

Microtubule-associated protein tau (MAPT) gene is compelling among the susceptibility genes of neurodegenerative diseases which include Alzheimer’s disease (AD), Parkinson’s disease (PD), progressive supranuclear palsy (PSP), corticobasal degeneration (CBD), frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS). Our meta-analysis aimed to find the association between MAPT and the risk of these diseases. Published literatures were retrieved from MEDLINE and other databases, and 82 case-control studies were recruited. Six haplotype tagging single-nucleotide polymorphisms (rs1467967, rs242557, rs3785883, rs2471738, del-In9 and rs7521) and haplotypes (H2 and H1c) were significantly associated with the above diseases. The odds ratios (ORs) and 95 % confidence intervals (CIs) were evaluated by comparison in minor and major allele frequency using the R software. This study demonstrated that different variants in MAPT were associated with AD (rs2471738: OR= 1.04, 95%CI = 1.00 - 1.09; H2: OR = 0.94, 95% CI = 0.91 - 0.97), PD (H2: OR = 0.76, 95% CI = 0.74 - 0.79), PSP (rs242557: OR = 1. 96, 95% CI = 1. 71 - 2.25; rs2471738: OR = 1. 85, 95% CI = 1. 48 - 2.31; H2: OR = 0.20, 95% CI = 0.18 - 0.23), CBD (rs242557: OR = 2.51, 95%CI = 1. 66 -3.78; rs2471738: OR = 2.07, 95%CI = 1. 32 -3.23; H2: OR = OR = 0.30, 95% CI = 0.23 - 0.41) and ALS (H2: OR = 0.92, 95% CI = 0.86 - 0.98) instead of FTD (H2: OR = 1.02, 95% CI = 0.78 - 1.32). In conclusion, MAPT is associated with risk of neurodegenerative diseases, suggesting crucial roles of tau in neurodegenerative processes.

Tau haplotypes support the Asian ancestry of the Roma population settled in the Basque Country

We examined tau haplotype frequencies in two different ethnical groups from the Basque Country (BC): Roma people and residents of European ancestry (general population). In addition, we analyzed the spatial distribution of tau haplotypes in Eurasian populations to explore the genetic affinities of the Romani groups living in Europe in a broader scope. The 17q21.31 genomic region was characterized through the genotyping of two diagnostic single nucleotide polymorphisms, SNPs (rs10514879 and rs199451), which allow the identification of H1 and H2 haplotypes. A significant heterozygous deficit was detected in the Romani for rs10514879. The H2 haplotype frequency proved to be more than twice in the BC general population (0.283) than in the Roma people (0.127). In contrast, H2 frequency proved to be very similar between Basque and Hungarian Romani, and similar to the H2 frequencies found in northwestern India and Pakistan as well. Several statistical analyses unveiled genetic structuring for the MAPT diversity, mirrored in a significant association between geography and genetic distances, with an upward trend of H2 haplotype frequencies from Asia to Europe. Yet, Roma samples did not fit into this general spatial patterning because of their discrepancy between geographical position and H2 frequency. Despite the long spatial coexistence in the Basque region between the residents of European ancestry and the Roma, the latter have preserved their Asian genetic ancestry. Bearing in mind the lack of geographical barriers between both ethnical groups, these findings support the notion that sociocultural mores might promote assortative matings in human populations.

Leucine-Rich Repeat Kinase (LRRK2) Genetics and Parkinson's Disease

The discovery of LRRK2 mutations as a cause of Parkinson's disease (PD), including the sporadic late-onset form, established the decisive role of genetics in the field of PD research. Among LRRK2 mutations, the G2019S, mostly lying in a haplotype originating from a common Middle Eastern ancestor, has been identified in different populations worldwide. The G2385R and R1628P variants represent validated risk factors for PD in Asian populations. Here, we describe in detail the origin, the present worldwide epidemiology, and the penetrance of LRRK2 mutations. Furthermore, this chapter aims to characterize other definitely/probably pathogenic mutations and risk variants of LRRK2. Finally, we provide some general guidelines for a LRRK2 genetic testing and counseling. In summary, LRRK2 discovery revolutionized the understanding of PD etiology and laid the foundation for a promising future of genetics in PD research.

The Role of MAPT Haplotype H2 and Isoform 1N/4R in Parkinsonism of Older Adults

Background and Objective

Recently, we have shown that the Parkinson’s disease (PD) susceptibility locus MAPT (microtubule associated protein tau) is associated with parkinsonism in older adults without a clinical diagnosis of PD. In this study, we investigated the relationship between parkinsonian signs and MAPT transcripts by assessing the effect of MAPT haplotypes on alternative splicing and expression levels of the most common isoforms in two prospective clinicopathologic studies of aging.

Materials and Methods

using regression analysis, controlling for age, sex, study and neuropathology, we evaluated 976 subjects with clinical, genotyping and brain pathology data for haplotype analysis. For transcript analysis, we obtained MAPT gene and isoform-level expression from the dorsolateral prefrontal cortex for 505 of these subjects.

Results

The MAPT H2 haplotype was associated with lower total MAPT expression (p = 1.2x10-14) and global parkinsonism at both study entry (p = 0.001) and proximate to death (p = 0.050). Specifically, haplotype H2 was primarily associated with bradykinesia in both assessments (p<0.001 and p = 0.008). MAPT total expression was associated with age and decreases linearly with advancing age (p<0.001). Analysing MAPT alternative splicing, the expression of 1N/4R isoform was inversely associated with global parkinsonism (p = 0.008) and bradykinesia (p = 0.008). Diminished 1N/4R isoform expression was also associated with H2 (p = 0.001).

Conclusions

Overall, our results suggest that age and H2 are associated with higher parkinsonism score and decreased total MAPT RNA expression. Additionally, we found that H2 and parkinsonism are associated with altered expression levels of specific isoforms. These findings may contribute to the understanding of the association between MAPT locus and parkinsonism in elderly subjects and in some extent to age-related neurodegenerative diseases.

Exploring the Effects of Genetic Variants on Clinical Profiles of Parkinson's Disease Assessed by the Unified Parkinson's Disease Rating Scale and the Hoehn-Yahr Stage

Many genetic variants have been linked to familial or sporadic Parkinson’s disease (PD), among which those identified in PARK16, BST1, SNCA, LRRK2, GBA and MAPT genes have been demonstrated to be the most common risk factors worldwide. Moreover, complex gene-gene and gene-environment interactions have been highlighted in PD pathogenesis. Compared to studies focusing on the predisposing effects of genes, there is a relative lack of research investigating how these genes and their interactions influence the clinical profiles of PD. In a cohort consisting of 2,011 Chinese Han PD patients, we selected 9 representative variants from the 6 above-mentioned common PD genes to analyze their main and epistatic effects on the Unified Parkinson’s Disease Rating Scale (UPDRS) and the Hoehn and Yahr (H-Y) stage of PD. With multiple linear regression models adjusting for medication status, disease duration, gender and age at onset, none of the variants displayed significant main effects on UPDRS or the H-Y scores. However, for gene-gene interaction analyses, 7 out of 37 pairs of variants showed significant or marginally significant associations with these scores. Among these, the GBA rs421016 (L444P)×LRRK2 rs33949390 (R1628P) interaction was consistently significant in relation to UPDRS III and UPDRS total (I+II+III), even after controlling for the family-wise error rate using False Discovery Rate (FDR-corrected p values are 0.0481 and 0.0070, respectively). Although the effects of the remaining pairs of variants did not survive the FDR correction, they showed marginally significant associations with either UPDRS or the H-Y stage (raw p<0.05). Our results highlight the importance of epistatic effects of multiple genes on the determination of PD clinical profiles and may have implications for molecular classification and personalized intervention of the disease.

Copy number variation analysis of the 17q21.31 region and its role in neurodegenerative diseases

The H1 haplotype of the 17q21.31 inversion polymorphism has been consistently associated with progressive supranuclear palsy, corticobasal degeneration, and Parkinson's disease in Caucasians. Recently, large polymorphic segmental duplications resulting into complex rearrangements at this locus with a high diversity range in human populations have been revealed. We sought to explore whether the two multi-allelic copy number variants that are present in the H1 clade (with segmental duplications of 300 and 218 kilobases in length) could be responsible for the known H1-related risk of developing these neurodegenerative disorders. A total of 857 Spanish subjects including 330 patients with Parkinson's disease, 96 with progressive supranuclear palsy, 55 with corticobasal degeneration, 51 dementia with Lewy bodies, and 325 neurologically healthy controls, were genotyped for the H1/H2 haplotype. Subsequently, the two copy number variants that are characteristic of the H1 haplotype were evaluated through a high-resolution approach based on droplet digital polymerase chain reaction, in all H1 homozygous subjects. The H1 allele was significantly overrepresented in all diagnostic groups compared with controls (Parkinson's disease, P = 0.0001; progressive supranuclear palsy, P = 1.22 × 10(-6) ; corticobasal degeneration, P = 0.0002; and dementia with Lewy bodies, P = 0.032). However, no dosage differences were found for any of the two copy number variants analyzed. The H1 haplotype is associated with the risk of several neurodegenerative disorders, including dementia with Lewy bodies. However, common structural diversity within the 17q21.31-H1 clade does not explain this genetic association.

Functional Impact and Evolution of a Novel Human Polymorphic Inversion That Disrupts a Gene and Creates a Fusion Transcript

Despite many years of study into inversions, very little is known about their functional consequences, especially in humans. A common hypothesis is that the selective value of inversions stems in part from their effects on nearby genes, although evidence of this in natural populations is almost nonexistent. Here we present a global analysis of a new 415-kb polymorphic inversion that is among the longest ones found in humans and is the first with clear position effects. This inversion is located in chromosome 19 and has been generated by non-homologous end joining between blocks of transposable elements with low identity. PCR genotyping in 541 individuals from eight different human populations allowed the detection of tag SNPs and inversion genotyping in multiple populations worldwide, showing that the inverted allele is mainly found in East Asia with an average frequency of 4.7%. Interestingly, one of the breakpoints disrupts the transcription factor gene ZNF257, causing a significant reduction in the total expression level of this gene in lymphoblastoid cell lines. RNA-Seq analysis of the effects of this expression change in standard homozygotes and inversion heterozygotes revealed distinct expression patterns that were validated by quantitative RT-PCR. Moreover, we have found a new fusion transcript that is generated exclusively from inverted chromosomes around one of the breakpoints. Finally, by the analysis of the associated nucleotide variation, we have estimated that the inversion was generated ~40,000–50,000 years ago and, while a neutral evolution cannot be ruled out, its current frequencies are more consistent with those expected for a deleterious variant, although no significant association with phenotypic traits has been found so far.

Since the discovery of chromosomal inversions almost 100 years ago, how they are maintained in natural populations has been a highly debated issue. One of the hypotheses is that inversion breakpoints could affect genes and modify gene expression levels, although evidence of this came only from laboratory mutants. In humans, a few inversions have been shown to associate with expression differences, but in all cases the molecular causes have remained elusive. Here, we have carried out a complete characterization of a new human polymorphic inversion and determined that it is specific to East Asian populations. In addition, we demonstrate that it disrupts the ZNF257 gene and, through the translocation of the first exon and regulatory sequences, creates a previously nonexistent fusion transcript, which together are associated to expression changes in several other genes. Finally, we investigate the potential evolutionary and phenotypic consequences of the inversion, and suggest that it is probably deleterious. This is therefore the first example of a natural polymorphic inversion that has position effects and creates a new chimeric gene, contributing to answer an old question in evolutionary biology.

Role for the microtubule-associated protein tau variant p.A152T in risk of α-synucleinopathies

OBJECTIVE

To assess the importance of MAPT variant p.A152T in the risk of synucleinopathies.

METHODS

In this case-control study, we screened a large global series of patients and controls, and assessed associations between p.A152T and disease risk. We included 3,229 patients with clinical Parkinson disease (PD), 442 with clinical dementia with Lewy bodies (DLB), 181 with multiple system atrophy (MSA), 832 with pathologically confirmed Lewy body disease (LBD), and 2,456 healthy controls.

RESULTS

The minor allele frequencies (MAF) in clinical PD cases (0.28%) and in controls (0.2%) were not found to be significantly different (odds ratio [OR] 1.37, 95% confidence interval [CI] 0.63-2.98, p = 0.42). However, a significant association was observed with clinical DLB (MAF 0.68%, OR 5.76, 95% CI 1.62-20.51, p = 0.007) and LBD (MAF 0.42%, OR 3.55, 95% CI 1.04-12.17, p = 0.04). Additionally, p.A152T was more common in patients with MSA compared to controls (MAF 0.55%, OR 4.68, 95% CI 0.85-25.72, p = 0.08) but this was not statistically significant and therefore should be interpreted with caution.

CONCLUSIONS

Overall, our findings suggest that MAPT p.A152T is a rare low penetrance variant likely associated with DLB that may be influenced by coexisting LBD and AD pathology. Given the rare nature of the variant, further studies with greater sample size are warranted and will help to fully explain the role of p.A152T in the pathogenesis of the synucleinopathies.

Human inversions and their functional consequences

Polymorphic inversions are a type of structural variants that are difficult to analyze owing to their balanced nature and the location of breakpoints within complex repeated regions. So far, only a handful of inversions have been studied in detail in humans and current knowledge about their possible functional effects is still limited. However, inversions have been related to phenotypic changes and adaptation in multiple species. In this review, we summarize the evidences of the functional impact of inversions in the human genome. First, given that inversions have been shown to inhibit recombination in heterokaryotes, chromosomes displaying different orientation are expected to evolve independently and this may lead to distinct gene-expression patterns. Second, inversions have a role as disease-causing mutations both by directly affecting gene structure or regulation in different ways, and by predisposing to other secondary arrangements in the offspring of inversion carriers. Finally, several inversions show signals of being selected during human evolution. These findings illustrate the potential of inversions to have phenotypic consequences also in humans and emphasize the importance of their inclusion in genome-wide association studies.

Complex and multi-allelic copy number variation in human disease

Hundreds of copy number variants are complex and multi-allelic, in that they have many structural alleles and have rearranged multiple times in the ancestors who contributed chromosomes to current humans. Not only are the relationships of these multi-allelic CNVs (mCNVs) to phenotypes generally unknown, but many mCNVs have not yet been described at the basic levels—alleles, allele frequencies, structural features—that support genetic investigation. To date, most reported disease associations to these variants have been ascertained through candidate gene studies. However, only a few associations have reached the level of acceptance defined by durable replications in many cohorts. This likely stems from longstanding challenges in making precise molecular measurements of the alleles individuals have at these loci. However, approaches for mCNV analysis are improving quickly, and some of the unique characteristics of mCNVs may assist future association studies. Their various structural alleles are likely to have different magnitudes of effect, creating a natural allelic series of growing phenotypic impact and giving investigators a set of natural predictions and testable hypotheses about the extent to which each allele of an mCNV predisposes to a phenotype. Also, mCNVs’ low-to-modest correlation to individual single-nucleotide polymorphisms (SNPs) may make it easier to distinguish between mCNVs and nearby SNPs as the drivers of an association signal, and perhaps, make it possible to preliminarily screen candidate loci, or the entire genome, for the many mCNV–disease relationships that remain to be discovered.

Interaction between Neuromelanin and Alpha-Synuclein in Parkinson's Disease

Parkinson's disease (PD) is a very common neurodegenerative disorder characterized by the accumulation of α-synuclein (α-syn) into Lewy body (LB) inclusions and the loss of neuronmelanin (NM) containing dopamine (DA) neurons in the substantia nigra (SN). Pathological α-syn and NM are two prominent hallmarks in this selective and progressive neurodegenerative disease. Pathological α-syn can induce dopaminergic neuron death by various mechanisms, such as inducing oxidative stress and inhibiting protein degradation systems. Therefore, to explore the factors that trigger α-syn to convert from a non-toxic protein to toxic one is a pivotal question to clarify the mechanisms of PD pathogenesis. Many triggers for pathological α-syn aggregation have been identified, including missense mutations in the α-syn gene, higher concentration, and posttranslational modifications of α-Syn. Recently, the role of NM in inducing α-syn expression and aggregation has been suggested as a mechanism for this pigment to modulate neuronal vulnerability in PD. NM may be responsible for PD and age-associated increase and aggregation in α-syn. Here, we reviewed our previous study and other recent findings in the area of interaction between NM and α-syn.

Lack of exacerbation of neurodegeneration in a double transgenic mouse model of mutant LRRK2 and tau

LRRK2 (leucine-rich repeat kinase) mutations constitute the most common cause of familial Parkinson's disease (PD). Microtubule-associated protein tau mutations cause a group of neurodegenerative diseases termed tauopathies. Genome-wide association studies show that, after α-synuclein, polymorphisms in the tau gene have the second strongest genetic association with PD. In a proportion of PD patients with LRRK2 mutations, and in several transgenic animal models of LRRK2, tau hyperphosphorylation and aggregation, rather than α-synuclein aggregation, are the most prominent neuropathologic findings. To further examine the relationship between LRRK2 and tau, we crossed LRRK2 R1441G BAC transgenic mice (Mus musculus) with tau P301S mutant transgenic mice and characterized their behavioral, neuropathological and biochemical phenotypes. We found that the combination of the two mutations does not increase tau hyperphosphorylation or aggregation nor does it exacerbate the behavioral and pathological deficits previously described in the tau P301S mice. The double-mutant mice had no shortening of lifespan and no worsening of motor or memory deficits. There was no increase in the aggregation of tau or α-synuclein. Dopaminergic neuron cell counts and striatal levels of dopamine and its metabolites were unaltered. There was no exacerbation of cell loss, microgliosis or astrogliosis in multiple brain regions. These results suggest that LRRK2 and tau do not interact to exacerbate behavioral, biochemical or pathological abnormalities in neurodegeneration and that LRRK2 and tau exert their pathogenic effects through independent mechanisms.

Polygenic determinants of Parkinson's disease in a Chinese population

It has been reported that some single-nucleotide polymorphisms (SNPs) are associated with the risk of Parkinson's disease (PD), but whether a combination of these SNPs would have a stronger association with PD than any individual SNP is unknown. Sixteen SNPs located in the 8 genes and/or loci (SNCA, LRRK2, MAPT, GBA, HLA-DR, BST1, PARK16, and PARK17) were analyzed in a Chinese cohort consisting of 1061 well-characterized PD patients and 1066 control subjects from Central South of Mainland China. We found that Rep1, rs356165, and rs11931074 in SNCA gene; G2385R in LRRK2 gene; rs4698412 in BST1 gene; rs1564282 in PARK17; and L444P in GBA gene were associated with PD with adjustment of sex and age (p < 0.05) in the analysis of 16 variants. PD risk increased when Rep1 and rs11931074, G2385R, rs1564282, rs4698412; rs11931074 and G2385R, rs1564282, rs4698412; G2385R and rs1564282, rs4698412; and rs1564282 and rs4698412 were combined for the association analysis. In addition, PD risk increased cumulatively with the increasing number of variants (odds ratio for carrying 3 variants, 3.494). In summary, we confirmed that Rep1, rs356165, and rs11931074 in SNCA gene, G2385R in LRRK2 gene, rs4698412 in BST1 gene, rs1564282 in PARK17, and L444P in GBA gene have an independent and combined significant association with PD. SNPs in these 4 genes have a cumulative effect with PD.

A graphical weighted power improving multiplicity correction approach for SNP selections

Controlling for the multiplicity effect is an essential part of determining statistical significance in large-scale single-locus association genome scans on Single Nucleotide Polymorphisms (SNPs). Bonferroni adjustment is a commonly used approach due to its simplicity, but is conservative and has low power for large-scale tests. The permutation test, which is a powerful and popular tool, is computationally expensive and may mislead in the presence of family structure. We propose a computationally efficient and powerful multiple testing correction approach for Linkage Disequilibrium (LD) based Quantitative Trait Loci (QTL) mapping on the basis of graphical weighted-Bonferroni methods. The proposed multiplicity adjustment method synthesizes weighted Bonferroni-based closed testing procedures into a powerful and versatile graphical approach. By tailoring different priorities for the two hypothesis tests involved in LD based QTL mapping, we are able to increase power and maintain computational efficiency and conceptual simplicity. The proposed approach enables strong control of the familywise error rate (FWER). The performance of the proposed approach as compared to the standard Bonferroni correction is illustrated by simulation and real data. We observe a consistent and moderate increase in power under all simulated circumstances, among different sample sizes, heritabilities, and number of SNPs. We also applied the proposed method to a real outbred mouse HDL cholesterol QTL mapping project where we detected the significant QTLs that were highlighted in the literature, while still ensuring strong control of the FWER.

Genetics and genomics of Parkinson's disease

Parkinson’s disease (PD) is a progressively debilitating neurodegenerative syndrome. Although best described as a movement disorder, the condition has prominent autonomic, cognitive, psychiatric, sensory and sleep components. Striatal dopaminergic innervation and nigral neurons are progressively lost, with associated Lewy pathology readily apparent on autopsy. Nevertheless, knowledge of the molecular events leading to this pathophysiology is limited. Current therapies offer symptomatic benefit but they fail to slow progression and patients continue to deteriorate. Recent discoveries in sporadic, Mendelian and more complex forms of parkinsonism provide novel insight into disease etiology; 28 genes, including those encoding alpha-synuclein (SNCA), leucine-rich repeat kinase 2 (LRRK2) and microtubule-associated protein tau (MAPT), have been linked and/or associated with PD. A consensus regarding the affected biological pathways and molecular processes has also started to emerge. In early-onset and more a typical PD, deficits in mitophagy pathways and lysosomal function appear to be prominent. By contrast, in more typical late-onset PD, chronic, albeit subtle, dysfunction in synaptic transmission, early endosomal trafficking and receptor recycling, as well as chaperone-mediated autophagy, provide a unifying synthesis of the molecular pathways involved. Disease-modification (neuroprotection) is no longer such an elusive goal given the unparalleled opportunity for diagnosis, translational neuroscience and therapeutic development provided by genetic discovery.

Holm multiple correction for large-scale gene-shape association mapping

Background

Linkage Disequilibrium (LD) is a powerful approach for the identification and characterization of morphological shape, which usually involves multiple genetic markers. However, multiple testing corrections substantially reduce the power of the associated tests. In addition, the principle component analysis (PCA), used to quantify the shape variations into several principal phenotypes, further increases the number of tests. As a result, a powerful multiple testing correction for simultaneous large-scale gene-shape association tests is an essential part of determining statistical significance. Bonferroni adjustments and permutation tests are the most popular approaches to correcting for multiple tests within LD based Quantitative Trait Loci (QTL) models. However, permutations are extremely computationally expensive and may mislead in the presence of family structure. The Bonferroni correction, though simple and fast, is conservative and has low power for large-scale testing.

Results

We propose a new multiple testing approach, constructed by combining an Intersection Union Test (IUT) with the Holm correction, which strongly controls the family-wise error rate (FWER) without any additional assumptions on the joint distribution of the test statistics or dependence structure of the markers. The power improvement for the Holm correction, as compared to the standard Bonferroni correction, is examined through a simulation study. A consistent and moderate increase in power is found under the majority of simulated circumstances, including various sample sizes, Heritabilities, and numbers of markers. The power gains are further demonstrated on real leaf shape data from a natural population of poplar, Populus szechuanica var tietica, where more significant QTL associated with morphological shape are detected than under the previously applied Bonferroni adjustment.

Conclusion

The Holm correction is a valid and powerful method for assessing gene-shape association involving multiple markers, which not only controls the FWER in the strong sense but also improves statistical power.

MAPT rs242562 and GSK3B rs334558 are associated with Parkinson's Disease in central China

Background

Microtubule-associated protein tau (MAPT) is a neuronal protein involved in the pathogenesis of several neurodegenerative diseases including Parkinson’s Disease (PD). Glycogen synthase kinase 3 beta (GSK3B) catalyzes phosphorylation in multiple sites of tau protein. However, whether or not there is any association between the GSK3B gene alteration, MAPT haplotype and PD remains unexplored in Chinese population, especially in central Chinese population.

Results

Here, we aimed at studying the effect of MAPT rs242562 and GSK3B rs334558 on the risk of PD by performing a case-control association study in central China. Our data showed that all PD patients and controls were of MAPT H1/H1 diplotype in our study, thus confirming that the distribution of the MAPT H1 haplotype is common in China. GG genotype of MAPT rs242562 serves protection effect on PD risk in central Chinese population, while genotype of GSK3B rs334558 showed no difference between PD patients and controls.

Conclusions

We conclude that the MAPT rs242562 is associated with PD in central China in the background of MAPT H1/H1 diplotype. The GG genotype of rs242562 displays protection against PD in subgroup with GSK3B rs334558 T carrier.

DNA methylation of the MAPT gene in Parkinson's disease cohorts and modulation by vitamin E in vitro

Parkinson's disease (PD) is a neurodegenerative disorder for which environmental factors influence disease risk and may act via an epigenetic mechanism. The microtubule-associated protein tau (MAPT) is a susceptibility gene for idiopathic PD. Methylation levels were determined by pyrosequencing of bisulfite-treated DNA in a leukocyte cohort (358 PD patients and 1084 controls) and in two brain cohorts (Brain1, comprising 69 cerebellum controls; and Brain2, comprising 3 brain regions from 28 PD patients and 12 controls). In vitro assays involved the transfection of methylated promoter-luciferase constructs or treatment with an exogenous micronutrient. In normal leukocytes, the MAPT H1/H2 diplotype and sex were predictors of MAPT methylation. Haplotype-specific pyrosequencing confirmed that the H1 haplotype had higher methylation than the H2 haplotype in normal leukocytes and brain tissues. MAPT methylation was negatively associated with MAPT expression in the Brain1 cohort and in transfected cells. Methylation levels differed between three normal brain regions (Brain2 cohort, putamen < cerebellum < anterior cingulate cortex). In PD samples, age at onset was positively associated with MAPT methylation in leukocytes. Moreover, there was hypermethylation in the cerebellum and hypomethylation in the putamen of PD patients compared with controls (Brain2 cohort). Finally, leukocyte methylation status was positively associated with blood vitamin E levels, and the effect was more significant in H2 haplotype carriers; this result was confirmed in cells that were exposed to 100 μM vitamin E. The significant effects of sex, diplotype, and brain region suggest that hypermethylation of the MAPT gene is neuroprotective by reducing MAPT expression. The effect of vitamin E on MAPT represents a possible gene-environment interaction.

Key role of CRF in the skin stress response system

The discovery of corticotropin-releasing factor (CRF) or CRH defining the upper regulatory arm of the hypothalamic-pituitary-adrenal (HPA) axis, along with the identification of the corresponding receptors (CRFRs 1 and 2), represents a milestone in our understanding of central mechanisms regulating body and local homeostasis. We focused on the CRF-led signaling systems in the skin and offer a model for regulation of peripheral homeostasis based on the interaction of CRF and the structurally related urocortins with corresponding receptors and the resulting direct or indirect phenotypic effects that include regulation of epidermal barrier function, skin immune, pigmentary, adnexal, and dermal functions necessary to maintain local and systemic homeostasis. The regulatory modes of action include the classical CRF-led cutaneous equivalent of the central HPA axis, the expression and function of CRF and related peptides, and the stimulation of pro-opiomelanocortin peptides or cytokines. The key regulatory role is assigned to the CRFR-1α receptor, with other isoforms having modulatory effects. CRF can be released from sensory nerves and immune cells in response to emotional and environmental stressors. The expression sequence of peptides includes urocortin/CRF→pro-opiomelanocortin→ACTH, MSH, and β-endorphin. Expression of these peptides and of CRFR-1α is environmentally regulated, and their dysfunction can lead to skin and systemic diseases. Environmentally stressed skin can activate both the central and local HPA axis through either sensory nerves or humoral factors to turn on homeostatic responses counteracting cutaneous and systemic environmental damage. CRF and CRFR-1 may constitute novel targets through the use of specific agonists or antagonists, especially for therapy of skin diseases that worsen with stress, such as atopic dermatitis and psoriasis.

FUS-mediated alternative splicing in the nervous system: consequences for ALS and FTLD

Mutations in fused in sarcoma (FUS) in a subset of patients with amyotrophic lateral sclerosis (ALS) linked this DNA/RNA-binding protein to neurodegeneration. Most of the mutations disrupt the nuclear localization signal which strongly suggests a loss-of-function pathomechanism, supported by cytoplasmic inclusions. FUS-positive neuronal cytoplasmic inclusions are also found in a subset of patients with frontotemporal lobar degeneration (FTLD). Here, we discuss recent data on the role of alternative splicing in FUS-mediated pathology in the central nervous system. Several groups have shown that FUS binds broadly to many transcripts in the brain and have also identified a plethora of putative splice targets; however, only ABLIM1, BRAF, Ewing sarcoma protein R1 (EWSR1), microtubule-associated protein tau (MAPT), NgCAM cell adhesion molecule (NRCAM), and netrin G1 (NTNG1) have been identified in at least three of four studies. Gene ontology analysis of all putative targets unanimously suggests a role in axon growth and cytoskeletal organization, consistent with the altered morphology of dendritic spines and axonal growth cones reported upon loss of FUS. Among the axonal targets, MAPT/tau and NTNG1 have been further validated in biochemical studies. The next challenge will be to confirm changes of FUS-mediated alternative splicing in patients and define their precise role in the pathophysiology of ALS and FTLD.

Protective effect of LRRK2 p.R1398H on risk of Parkinson's disease is independent of MAPT and SNCA variants

The best validated susceptibility variants for Parkinson's disease are located in the α-synuclein (SNCA) and microtubule-associated protein tau (MAPT) genes. Recently, a protective p.N551K-R1398H-K1423K haplotype in the leucine-rich repeat kinase 2 (LRRK2) gene was identified, with p.R1398H appearing to be the most likely functional variant. To date, the consistency of the protective effect of LRRK2 p.R1398H across MAPT and SNCA variant genotypes has not been assessed. To address this, we examined 4 SNCA variants (rs181489, rs356219, rs11931074, and rs2583988), the MAPT H1-haplotype-defining variant rs1052553, and LRRK2 p.R1398H (rs7133914) in Caucasian (n = 10,322) and Asian (n = 2289) series. There was no evidence of an interaction of LRRK2 p.R1398H with MAPT or SNCA variants (all p ≥ 0.10); the protective effect of p.R1398H was observed at similar magnitude across MAPT and SNCA genotypes, and the risk effects of MAPT and SNCA variants were observed consistently for LRRK2 p.R1398H genotypes. Our results indicate that the association of LRRK2 p.R1398H with Parkinson's disease is independent of SNCA and MAPT variants, and vice versa, in Caucasian and Asian populations.

Advances in the genetics of Parkinson disease

Parkinson disease (PD) is a multifactorial neurodegenerative disease that was long considered the result of environmental factors. In the past 15 years, however, a genetic aetiology for PD has begun to emerge. Here, we review results from linkage and next-generation sequencing studies of familial parkinsonism, as well as candidate gene and genome-wide association findings in sporadic PD. In these studies, many of the genetic findings overlap, despite different designs and study populations, highlighting novel therapeutic targets. The molecular results delineate a sequence of pathological events whereby deficits in synaptic exocytosis and endocytosis, endosomal trafficking, lysosome-mediated autophagy and mitochondrial maintenance increase susceptibility to PD. These discoveries provide the rationale, molecular insight and research tools to develop neuroprotective and disease-modifying therapies.

A prognostic view on the application of individualized genomics in Parkinson's disease

Parkinson's disease is a common age-related progressive neurodegenerative disorder. Over the last 15 years advances have been made in our understanding of the etiology of the disease, with the greatest insights perhaps coming from genetic studies. The identification of a number of genes that harbor pathogenic mutations causing Parkinson's disease have on the whole driven the development of disease model systems and nominated a number of therapeutic targets. As we move towards an era of personalized medicine, genetic determinants will become even more crucial for accurate diagnosis, and assessing prognosis and outcomes. The individual genomic profile and risk assessments will in the long-term determine clinical trial participation, treatment plans and therapeutic dosing. Herein we discuss the status of genetics in Parkinson's disease and how these factors may affect the patient care in the future.