Search

Huntington disease (HD) reflects the dominant consequences of a CAG-repeat expansion in HTT. Analysis of common SNP-based haplotypes has revealed that most European HD subjects have distinguishable HTT haplotypes on their normal and disease chromosomes and that ∼50% of the latter share the same major HD haplotype. We reasoned that sequence-level investigation of this founder haplotype could provide significant insights into the history of HD and valuable information for gene-targeting approaches. Consequently, we performed whole-genome sequencing of HD and control subjects from four independent families in whom the major European HD haplotype segregates with the disease. Analysis of the full-sequence-based HTT haplotype indicated that these four families share a common ancestor sufficiently distant to have permitted the accumulation of family-specific variants. Confirmation of new CAG-expansion mutations on this haplotype suggests that unlike most founders of human disease, the common ancestor of HD-affected families with the major haplotype most likely did not have HD. Further, availability of the full sequence data validated the use of SNP imputation to predict the optimal variants for capturing heterozygosity in personalized allele-specific gene-silencing approaches. As few as ten SNPs are capable of revealing heterozygosity in more than 97% of European HD subjects. Extension of allele-specific silencing strategies to the few remaining homozygous individuals is likely to be achievable through additional known SNPs and discovery of private variants by complete sequencing of HTT. These data suggest that the current development of gene-based targeting for HD could be extended to personalized allele-specific approaches in essentially all HD individuals of European ancestry.

The National Heart, Lung, and Blood Institute Family Heart Study (FHS) genome-wide linkage scan identified a region of chromosome 7q31-34 with a lod score of 4.9 for BMI at D7S1804 (131.9 Mb). We report the results of linkage and association to BMI in this region for two independent FHS samples. The first sample includes 225 FHS pedigrees with evidence of linkage to 7q31-34, using 1,132 single-nucleotide polymorphisms (SNPs) and 7 microsatellites. The second represents a case-control sample (318 cases; BMI >25 and 325 controls; BMI <25) derived from unrelated FHS participants who were not part of the genome scan. The latter set was genotyped for 606 SNPs, including 37 SNPs with prior evidence for association in the linked families. Although variance components linkage analysis using only SNPs generated a peak lod score that coincided with the original linkage scan at 131.9 Mb, a conditional linkage analysis showed evidence of a second quantitative trait locus (QTL) near 143 cM influencing BMI. Three SNPs (rs161339, rs12673281, and rs1993068) located near the three genes pleiotrophin (PTN), diacylglycerol (DAG) kinase iota (DGK iota), and cholinergic receptor, muscarinic 2 (CHRM2) demonstrated significant association in both linked families (P = 0.0005, 0.002, and 0.03, respectively) and the case-control sample (P = 0.01, 0.0003, and 0.03, respectively), regardless of the genetic model tested. These findings suggest that several genes may be associated with BMI in the 7q31-34 region. [ABSTRACT FROM AUTHOR]

Recent studies have demonstrated the importance of sex effects on the underlying genetic architecture of insulin-related traits. To explore sex-specific genetic effects on fasting insulin, we tested for genotype-by-sex interaction and conducted linkage analysis of fasting insulin in Hypertension Genetic Epidemiology Network families. Hypertensive siblings and their first-degree relatives were recruited from five field centers. We performed a genome scan for quantitative trait loci influencing fasting insulin among 1,505 European Americans and 1,616 African Americans without diabetes. Sex-stratified linear regression models, adjusted for race, center, and age, were explored. The Mammalian Genotyping Service typed 391 microsatellite markers, spaced roughly 9 cM. Variance component linkage analysis was performed in SOLAR using ethnic-specific marker allele frequencies and multipoint IBDs calculated in MERLIN. We detected a quantitative trait locus influencing fasting insulin in female subjects (logarithm of odds [LOD] = 3.4) on chromosome 2 at 95 cM (between GATA69E12 and GATA71G04) but not in male subjects (LOD = 0.0, P for interaction = 0.007). This sex-specific signal at 2p13.2 was detected in both European-American (LOD = 2.1) and African-American (LOD = 1.2) female subjects. Our findings overlap with several other linkage reports of insulin-related traits and demonstrate the importance of considering complex context-dependent interactions in the search for insulin-related genes. [ABSTRACT FROM AUTHOR]

Genome-wide linkage scans were carried out using a multipoint variance components method in white and black families of the NHLBI Hypertension Genetic Epidemiology Network (HyperGEN) study to identify quantitative trait loci (QTLs) for pancreatic beta-cell function and insulin sensitivity estimated through the newly released nonlinear computer version of homeostasis model assessment 2. Participants fasting <8 h, with diagnosed type 2 diabetes, or taking blood glucose or blood lipid-lowering medications were excluded. Both phenotypes were adjusted separately by race and sex for the effects of age, BMI, and field center before linkage scans using 370 microsatellite markers were performed. A total of 685 white families (1,180 sibpairs) and 773 black families (775 sibpairs) were evaluated as well as subsets including 267 obese white families (757 sibpairs) and 427 obese black families (599 sibpairs) identified through tree-linkage analyses using interacting covariates of age, sex, and BMI. For beta-cell function in the obese white families, significant (logarithm of odds [LOD] score >3.6) evidence supporting linkages was detected on chromosome 8q24 at D8S1179 (135 cM, LOD score 4.2, empirical P = 0.002) and at D8S1128 (140 cM, LOD score 3.7, empirical P = 0.003). In addition, two regions supported linkage for insulin sensitivity index in the obese black families on chromosome 7q11 at D7S3046 (79 cM, LOD score 3.0, empirical P = 0.018) and on chromosome 6q26 at D6S1277 (173 cM, LOD score 3.0, empirical P = 0.018). Reducing clinical heterogeneity using obesity data and improved estimates of beta-cell function and insulin sensitivity may have permitted identification of a QTL on chromosome 8q24 for beta-cell function in the presence of estimated insulin resistance and a QTL on chromosome 7q11 for insulin sensitivity. These regions replicate previous reports for type 2 diabetes-associated traits. [ABSTRACT FROM AUTHOR]

Type 2 diabetes susceptibility is determined by multiple genetic and environmental factors. Genome-wide linkage scans have localized common regions, possibly harboring susceptibility genes on chromosomes 1, 2, 12, and 20. Variability in linkage findings underscores the probable genetic heterogeneity of type 2 diabetes. Thus, we conducted a genome scan of diabetes status using maximum likelihood methods that model affection status by a liability threshold model. Hypertensive sibships and their offspring and/or parents in the Hypertension Genetic Epidemiology Network study were recruited from five field centers. The diabetes phenotype was derived using the World Health Organization criteria and adjusted for race/study center, age, age2, sex, and with and without percent body fat. In total, 567 diabetic participants were identified in 437 families. Variance component linkage analysis was performed among 1,545 Caucasians and 1,608 African Americans using race-specific marker allele frequencies. We detected a quantitative trait loci (QTLs) influencing diabetes variance (logarithm of odds = 3.4) on chromosome 22, which overlaps a positive type 2 diabetes finding among Canadian Oji-Cree Indians. We also observed suggestive evidence for linkage on chromosomes 1, 2, 5, 8, 14, 17, and 19. The identification and replication of type 2 diabetes QTLs will bring us closer to the detection of functional genes that influence diabetes susceptibility. [ABSTRACT FROM AUTHOR]

Transposable Elements (TEs) are implicated in aging and neurodegenerative disorders, but the impact of brain TE RNA dynamics on these phenomena is not fully understood. Therefore, we quantified TE RNA changes in aging post-mortem human and mouse brains and in the neurodegenerative disorders Huntington's Disease (HD) and Parkinson's Disease (PD). We tracked TE small RNAs (smRNAs) expression landscape to assess the relationship to the active processing from TE long RNAs (lnRNAs). Human brain transcriptomes from the BrainSpan Atlas displayed a significant shift of TE smRNA patterns at age 20 years, whereas aging mouse brains lacked any such marked change, despite clear shift in aging-associated mRNA levels. Human frontal cortex displayed pronounced sense TE smRNAs during aging with a negative relationship between the TE smRNAs and lnRNAs indicative of age associated regulatory effects. Our analysis revealed TE smRNAs dysregulation in HD, while PD showed a stronger impact on TE lnRNAs, potentially correlating with the early average age of death for HD relative to PD. Furthermore, TE-silencing factor TRIM28 was down-regulated only in aging human brains, possibly explaining the lack of substantial TE RNA changes in aging mouse brains. Our study suggests brain TE RNAs may serve as novel biomarkers of human brain aging and neurodegenerative disorders.

Huntington's disease (HD), due to expansion of a CAG repeat in HTT , is representative of a growing number of disorders involving somatically unstable short tandem repeats. We find that overlapping and distinct genetic modifiers of clinical landmarks and somatic expansion in blood DNA reveal an underlying complexity and cell-type specificity to the mismatch repair-related processes that influence disease timing. Differential capture of non-DNA-repair gene modifiers by multiple measures of cognitive and motor dysfunction argues additionally for cell-type specificity of pathogenic processes. Beyond trans modifiers, differential effects are also illustrated at HTT by a 5'-UTR variant that promotes somatic expansion in blood without influencing clinical HD, while, even after correcting for uninterrupted CAG length, a synonymous sequence change at the end of the CAG repeat dramatically hastens onset of motor signs without increasing somatic expansion. Our findings are directly relevant to therapeutic suppression of somatic expansion in HD and related disorders and provide a route to define the individual neuronal cell types that contribute to different HD clinical phenotypes.

Mitochondrial dysfunction is associated with neuronal damage in Huntington's disease (HD), but the precise mechanism of mitochondria-dependent pathogenesis is not understood yet. Herein, we found that colocalization of XIAP and p53 was prominent in the cytosolic compartments of normal subjects but reduced in HD patients and HD transgenic animal models. Overexpression of mutant Huntingtin (mHTT) reduced XIAP levels and elevated mitochondrial localization of p53 in striatal cells in vitro and in vivo. Interestingly, XIAP interacted directly with the C-terminal domain of p53 and decreased its stability via autophagy. Overexpression of XIAP prevented mitochondrially targeted-p53 (Mito-p53)-induced mitochondrial oxidative stress and striatal cell death, whereas, knockdown of XIAP exacerbated Mito-p53-induced neuronal damage in vitro. In vivo transduction of AAV-shRNA XIAP in the dorsal striatum induced rapid onset of disease and reduced the lifespan of HD transgenic (N171-82Q) mice compared to WT littermate mice. XIAP dysfunction led to ultrastructural changes of the mitochondrial cristae and nucleus morphology in striatal cells. Knockdown of XIAP exacerbated neuropathology and motor dysfunctions in N171-82Q mice. In contrast, XIAP overexpression improved neuropathology and motor behaviors in both AAV-mHTT-transduced mice and N171-82Q mice. Our data provides a molecular and pathological mechanism that deregulation of XIAP triggers mitochondria dysfunction and other neuropathological processes via the neurotoxic effect of p53 in HD. Together, the XIAP-p53 pathway is a novel pathological marker and can be a therapeutic target for improving the symptoms in HD., (Copyright © 2021. Published by Elsevier Ltd.)

Objective: The aim of this study was to search for genes/variants that modify the effect of LRRK2 mutations in terms of penetrance and age-at-onset of Parkinson's disease., Methods: We performed the first genomewide association study of penetrance and age-at-onset of Parkinson's disease in LRRK2 mutation carriers (776 cases and 1,103 non-cases at their last evaluation). Cox proportional hazard models and linear mixed models were used to identify modifiers of penetrance and age-at-onset of LRRK2 mutations, respectively. We also investigated whether a polygenic risk score derived from a published genomewide association study of Parkinson's disease was able to explain variability in penetrance and age-at-onset in LRRK2 mutation carriers., Results: A variant located in the intronic region of CORO1C on chromosome 12 (rs77395454; p value = 2.5E-08, beta = 1.27, SE = 0.23, risk allele: C) met genomewide significance for the penetrance model. Co-immunoprecipitation analyses of LRRK2 and CORO1C supported an interaction between these 2 proteins. A region on chromosome 3, within a previously reported linkage peak for Parkinson's disease susceptibility, showed suggestive associations in both models (penetrance top variant: p value = 1.1E-07; age-at-onset top variant: p value = 9.3E-07). A polygenic risk score derived from publicly available Parkinson's disease summary statistics was a significant predictor of penetrance, but not of age-at-onset., Interpretation: This study suggests that variants within or near CORO1C may modify the penetrance of LRRK2 mutations. In addition, common Parkinson's disease associated variants collectively increase the penetrance of LRRK2 mutations. ANN NEUROL 2021;90:82-94., (© 2021 The Authors. Annals of Neurology published by Wiley Periodicals LLC on behalf of American Neurological Association.)

Previous studies on Parkinson's disease mechanisms have shown dysregulated extracellular transport of α-synuclein and growth factors in the extracellular space. In the human brain these consist of perineuronal nets, interstitial matrices, and basement membranes, each composed of a set of collagens, non-collagenous glycoproteins, proteoglycans, and hyaluronan. The manner by which amyloidogenic proteins spread extracellularly, become seeded, oligomerize, and are taken up by cells, depends on intricate interactions with extracellular matrix molecules. We sought to assess the alterations to structure of glycosaminoglycans and proteins that occur in PD brain relative to controls of similar age. We found that PD differs markedly from normal brain in upregulation of extracellular matrix structural components including collagens, proteoglycans and glycosaminoglycan binding molecules. We also observed that levels of hemoglobin chains, possibly related to defects in iron metabolism, were enriched in PD brains. These findings shed important new light on disease processes that occur in association with PD.

Background: Huntington's disease (HD) is an inherited neurodegenerative disorder caused by an expanded CAG repeat in the HTT gene. It is diagnosed following a standardized examination of motor control and often presents with cognitive decline and psychiatric symptoms. Recent studies have detected genetic loci modifying the age at onset of motor symptoms in HD, but genetic factors influencing cognitive and psychiatric presentations are unknown., Methods: We tested the hypothesis that psychiatric and cognitive symptoms in HD are influenced by the same common genetic variation as in the general population by 1) constructing polygenic risk scores from large genome-wide association studies of psychiatric and neurodegenerative disorders and of intelligence and 2) testing for correlation with the presence of psychiatric and cognitive symptoms in a large sample (n = 5160) of patients with HD., Results: Polygenic risk score for major depression was associated specifically with increased risk of depression in HD, as was schizophrenia risk score with psychosis and irritability. Cognitive impairment and apathy were associated with reduced polygenic risk score for intelligence., Conclusions: Polygenic risk scores for psychiatric disorders, particularly depression and schizophrenia, are associated with increased risk of the corresponding psychiatric symptoms in HD, suggesting a common genetic liability. However, the genetic liability to cognitive impairment and apathy appears to be distinct from other psychiatric symptoms in HD. No associations were observed between HD symptoms and risk scores for other neurodegenerative disorders. These data provide a rationale for treatments effective in depression and schizophrenia to be used to treat depression and psychotic symptoms in HD., (Copyright © 2020 Society of Biological Psychiatry. Published by Elsevier Inc. All rights reserved.)

Background: The mechanisms underlying neurodegeneration in the striatum of Huntingon's Disease (HD) brain are currently unknown. While the striatum is massively degenerated in symptomatic individuals, which makes cellular characterization difficult, it is largely intact in asymptomatic HD gene positive (HD+) individuals. Unfortunately, as striatal tissue samples from HD+ individuals are exceedingly rare, recent focus has been on the Brodmann Area 9 (BA9), a relatively unaffected region, as a surrogate tissue. In this study, we analyze gene expression in caudate nucleus (CAU) from two HD+ individuals and compare the results with healthy and symptomatic HD brains., Methods: High-throughput mRNA sequencing (mRNA-Seq) datasets were generated from post-mortem CAU of 2 asymptomatic HD+ individuals and compared with 26 HD and 56 neurologically normal controls. Datasets were analyzed using a custom bioinformatic analysis pipeline to identify and interpret differentially expressed (DE) genes. Results were compared to publicly available brain mRNA-Seq datasets from the Genotype-Tissue Expression (GTEx) project. The analysis employed current state of the art bioinformatics tools and tailored statistical and machine learning methods., Results: The transcriptional profiles in HD+ CAU and HD BA9 samples are highly similar. Differentially expressed (DE) genes related to the heat shock response, particularly HSPA6 and HSPA1A, are common between regions. The most perturbed pathways show extensive agreement when comparing disease with control. A random forest classifier predicts that the two HD+ CAU samples strongly resemble HD BA9 and not control BA9. Nonetheless, when genes were prioritized by their specificity to HD+ CAU, pathways spanning many biological processes emerge. Comparison of HD+ BA9 with HD BA9 identified NPAS4 and REST1/2 as potential early responders to disease and reflect the active disease process., Conclusions: The caudate nucleus in HD brain is dramatically affected prior to symptom onset. Gene expression patterns observed in the HD BA9 are also present in the CAU, suggesting a common response to disease. Substantial caudate-specific differences implicate many different biological pathways including metabolism, protein folding, inflammation, and neurogenic processes. While these results are at best trends due to small sample sizes, these results nonetheless provide the most detailed insight to date into the primary HD disease process.

Nearly 100 loci have been identified for pulmonary function, almost exclusively in studies of European ancestry populations. We extend previous research by meta-analyzing genome-wide association studies of 1000 Genomes imputed variants in relation to pulmonary function in a multiethnic population of 90,715 individuals of European (N = 60,552), African (N = 8429), Asian (N = 9959), and Hispanic/Latino (N = 11,775) ethnicities. We identify over 50 additional loci at genome-wide significance in ancestry-specific or multiethnic meta-analyses. Using recent fine-mapping methods incorporating functional annotation, gene expression, and differences in linkage disequilibrium between ethnicities, we further shed light on potential causal variants and genes at known and newly identified loci. Several of the novel genes encode proteins with predicted or established drug targets, including KCNK2 and CDK12. Our study highlights the utility of multiethnic and integrative genomics approaches to extend existing knowledge of the genetics of lung function and clinical relevance of implicated loci.

Objective: To investigate the feasibility of microRNA (miRNA) levels in CSF as biomarkers for prodromal Huntington disease (HD)., Methods: miRNA levels were measured in CSF from 60 PREDICT-HD study participants using the HTG protocol. Using a CAG-Age Product score, 30 prodromal HD participants were selected based on estimated probability of imminent clinical diagnosis of HD (i.e., low, medium, high; n = 10/group). For comparison, participants already diagnosed (n = 15) and healthy controls (n = 15) were also selected., Results: A total of 2,081 miRNAs were detected and 6 were significantly increased in the prodromal HD gene expansion carriers vs controls at false discovery rate q < 0.05 (miR-520f-3p, miR-135b-3p, miR-4317, miR-3928-5p, miR-8082, miR-140-5p). Evaluating the miRNA levels in each of the HD risk categories, all 6 revealed a pattern of increasing abundance from control to low risk, and from low risk to medium risk, which then leveled off from the medium to high risk and HD diagnosed groups., Conclusions: This study reports miRNAs as CSF biomarkers of prodromal and diagnosed HD. Importantly, miRNAs were detected in the prodromal HD groups furthest from diagnosis where treatments are likely to be most consequential and meaningful. The identification of potential biomarkers in the disease prodrome may prove useful in evaluating treatments that may postpone disease onset., Clinicaltrialsgov Identifier: NCT00051324., (Copyright © 2017 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology.)

Huntington's and Parkinson's Diseases (HD and PD) are neurodegenerative disorders that share some pathological features but are disparate in others. For example, while both diseases are marked by aberrant protein aggregation in the brain, the specific proteins that aggregate and types of neurons affected differ. A better understanding of the molecular similarities and differences between these two diseases may lead to a more complete mechanistic picture of both the individual diseases and the neurodegenerative process in general. We sought to characterize the common transcriptional signature of HD and PD as well as genes uniquely implicated in each of these diseases using mRNA-Seq data from post mortem human brains in comparison to neuropathologically normal controls. The enriched biological pathways implicated by HD differentially expressed genes show remarkable consistency with those for PD differentially expressed genes and implicate the common biological processes of neuroinflammation, apoptosis, transcriptional dysregulation, and neuron-associated functions. Comparison of the differentially expressed (DE) genes highlights a set of consistently altered genes that span both diseases. In particular, processes involving nuclear factor kappa-light-chain-enhancer of activated B cells (NFkB) and transcription factor cAMP response element-binding protein (CREB) are the most prominent among the genes common to HD and PD. When the combined HD and PD data are compared to controls, relatively few additional biological processes emerge as significantly enriched, suggesting that most pathways are independently seen within each disorder. Despite showing comparable numbers of DE genes, DE genes unique to HD are enriched in far more coherent biological processes than the DE genes unique to PD, suggesting that PD may represent a more heterogeneous disorder. The complexity of the biological processes implicated by this analysis provides impetus for the development of better experimental models to validate the results.

Crohn's disease (CD), a form of inflammatory bowel disease, has a higher prevalence in Ashkenazi Jewish than in non-Jewish European populations. To define the role of nonsynonymous mutations, we performed exome sequencing of Ashkenazi Jewish patients with CD, followed by array-based genotyping and association analysis in 2066 CD cases and 3633 healthy controls. We detected association signals in the LRRK2 gene that conferred risk for CD (N2081D variant, P = 9.5 × 10 -10 ) or protection from CD (N551K variant, tagging R1398H-associated haplotype, P = 3.3 × 10 -8 ). These variants affected CD age of onset, disease location, LRRK2 activity, and autophagy. Bayesian network analysis of CD patient intestinal tissue further implicated LRRK2 in CD pathogenesis. Analysis of the extended LRRK2 locus in 24,570 CD cases, patients with Parkinson's disease (PD), and healthy controls revealed extensive pleiotropy, with shared genetic effects between CD and PD in both Ashkenazi Jewish and non-Jewish cohorts. The LRRK2 N2081D CD risk allele is located in the same kinase domain as G2019S, a mutation that is the major genetic cause of familial and sporadic PD. Like the G2019S mutation, the N2081D variant was associated with increased kinase activity, whereas neither N551K nor R1398H variants on the protective haplotype altered kinase activity. We also confirmed that R1398H, but not N551K, increased guanosine triphosphate binding and hydrolyzing enzyme (GTPase) activity, thereby deactivating LRRK2. The presence of shared LRRK2 alleles in CD and PD provides refined insight into disease mechanisms and may have major implications for the treatment of these two seemingly unrelated diseases., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Huntington's disease (HD) is an autosomal dominant neurodegenerative disease caused by expansion of a CAG trinucleotide repeat in HTT, resulting in an extended polyglutamine tract in huntingtin. We and others have previously determined that the HD-causing expansion occurs on multiple different haplotype backbones, reflecting more than one ancestral origin of the same type of mutation. In view of the therapeutic potential of mutant allele-specific gene silencing, we have compared and integrated two major systems of HTT haplotype definition, combining data from 74 sequence variants to identify the most frequent disease-associated and control chromosome backbones and revealing that there is potential for additional resolution of HD haplotypes. We have used the large collection of 4078 heterozygous HD subjects analyzed in our recent genome-wide association study of HD age at onset to estimate the frequency of these haplotypes in European subjects, finding that common genetic variation at HTT can distinguish the normal and CAG-expanded chromosomes for more than 95% of European HD individuals. As a resource for the HD research community, we have also determined the haplotypes present in a series of publicly available HD subject-derived fibroblasts, induced pluripotent cells, and embryonic stem cells in order to facilitate efforts to develop inclusive methods of allele-specific HTT silencing applicable to most HD patients. Our data providing genetic guidance for therapeutic gene-based targeting will significantly contribute to the developments of rational treatments and implementation of precision medicine in HD.

Huntington's disease (HD) is a dominantly inherited neurodegenerative disease caused by an expanded CAG repeat in HTT. Many clinical characteristics of HD such as age at motor onset are determined largely by the size of HTT CAG repeat. However, emerging evidence strongly supports a role for other genetic factors in modifying the disease pathogenesis driven by mutant huntingtin. A recent genome-wide association analysis to discover genetic modifiers of HD onset age provided initial evidence for modifier loci on chromosomes 8 and 15 and suggestive evidence for a locus on chromosome 3. Here, genotyping of candidate single nucleotide polymorphisms in a cohort of 3,314 additional HD subjects yields independent confirmation of the former two loci and moves the third to genome-wide significance at MLH1, a locus whose mouse orthologue modifies CAG length-dependent phenotypes in a Htt-knock-in mouse model of HD. Both quantitative and dichotomous association analyses implicate a functional variant on ∼32% of chromosomes with the beneficial modifier effect that delays HD motor onset by 0.7 years/allele. Genomic DNA capture and sequencing of a modifier haplotype localize the functional variation to a 78 kb region spanning the 3'end of MLH1 and the 5'end of the neighboring LRRFIP2, and marked by an isoleucine-valine missense variant in MLH1. Analysis of expression Quantitative Trait Loci (eQTLs) provides modest support for altered regulation of MLH1 and LRRFIP2, raising the possibility that the modifier affects regulation of both genes. Finally, polygenic modification score and heritability analyses suggest the existence of additional genetic modifiers, supporting expanded, comprehensive genetic analysis of larger HD datasets., (© The Author 2017. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

We have previously shown that exon 1 of the huntingtin gene does not always splice to exon 2 resulting in the production of a small polyadenylated mRNA (HTTexon1) that encodes the highly pathogenic exon 1 HTT protein. The level of this read-through product is proportional to CAG repeat length and is present in all knock-in mouse models of Huntington's disease (HD) with CAG lengths of 50 and above and in the YAC128 and BACHD mouse models, both of which express a copy of the human HTT gene. We have now developed specific protocols for the quantitative analysis of the transcript levels of HTTexon1 in human tissue and applied these to a series of fibroblast lines and post-mortem brain samples from individuals with either adult-onset or juvenile-onset HD. We found that the HTTexon1 mRNA is present in fibroblasts from juvenile HD patients and can also be readily detected in the sensory motor cortex, hippocampus and cerebellum of post-mortem brains from HD individuals, particularly in those with early onset disease. This finding will have important implications for strategies to lower mutant HTT levels in patients and the design of future therapeutics.

Background: Next generation sequencing provides a count of RNA molecules in the form of short reads, yielding discrete, often highly non-normally distributed gene expression measurements. Although Negative Binomial (NB) regression has been generally accepted in the analysis of RNA sequencing (RNA-Seq) data, its appropriateness has not been exhaustively evaluated. We explore logistic regression as an alternative method for RNA-Seq studies designed to compare cases and controls, where disease status is modeled as a function of RNA-Seq reads using simulated and Huntington disease data. We evaluate the effect of adjusting for covariates that have an unknown relationship with gene expression. Finally, we incorporate the data adaptive method in order to compare false positive rates., Results: When the sample size is small or the expression levels of a gene are highly dispersed, the NB regression shows inflated Type-I error rates but the Classical logistic and Bayes logistic (BL) regressions are conservative. Firth's logistic (FL) regression performs well or is slightly conservative. Large sample size and low dispersion generally make Type-I error rates of all methods close to nominal alpha levels of 0.05 and 0.01. However, Type-I error rates are controlled after applying the data adaptive method. The NB, BL, and FL regressions gain increased power with large sample size, large log2 fold-change, and low dispersion. The FL regression has comparable power to NB regression., Conclusions: We conclude that implementing the data adaptive method appropriately controls Type-I error rates in RNA-Seq analysis. Firth's logistic regression provides a concise statistical inference process and reduces spurious associations from inaccurately estimated dispersion parameters in the negative binomial framework.

Genome-wide association studies (GWASs) have revealed numerous loci for areal bone mineral density (aBMD). We completed the first GWAS meta-analysis (n = 15,275) of lumbar spine volumetric BMD (vBMD) measured by quantitative computed tomography (QCT), allowing for examination of the trabecular bone compartment. SNPs that were significantly associated with vBMD were also examined in two GWAS meta-analyses to determine associations with morphometric vertebral fracture (n = 21,701) and clinical vertebral fracture (n = 5893). Expression quantitative trait locus (eQTL) analyses of iliac crest biopsies were performed in 84 postmenopausal women, and murine osteoblast expression of genes implicated by eQTL or by proximity to vBMD-associated SNPs was examined. We identified significant vBMD associations with five loci, including: 1p36.12, containing WNT4 and ZBTB40; 8q24, containing TNFRSF11B; and 13q14, containing AKAP11 and TNFSF11. Two loci (5p13 and 1p36.12) also contained associations with radiographic and clinical vertebral fracture, respectively. In 5p13, rs2468531 (minor allele frequency [MAF] = 3%) was associated with higher vBMD (β = 0.22, p = 1.9 × 10 -8 ) and decreased risk of radiographic vertebral fracture (odds ratio [OR] = 0.75; false discovery rate [FDR] p = 0.01). In 1p36.12, rs12742784 (MAF = 21%) was associated with higher vBMD (β = 0.09, p = 1.2 × 10 -10 ) and decreased risk of clinical vertebral fracture (OR = 0.82; FDR p = 7.4 × 10 -4 ). Both SNPs are noncoding and were associated with increased mRNA expression levels in human bone biopsies: rs2468531 with SLC1A3 (β = 0.28, FDR p = 0.01, involved in glutamate signaling and osteogenic response to mechanical loading) and rs12742784 with EPHB2 (β = 0.12, FDR p = 1.7 × 10 -3 , functions in bone-related ephrin signaling). Both genes are expressed in murine osteoblasts. This is the first study to link SLC1A3 and EPHB2 to clinically relevant vertebral osteoporosis phenotypes. These results may help elucidate vertebral bone biology and novel approaches to reducing vertebral fracture incidence. © 2016 American Society for Bone and Mineral Research., (© 2016 American Society for Bone and Mineral Research.)

Background: Leucine-rich repeat kinase 2 (LRRK2) mutation 6055G→A (Gly2019Ser) accounts for roughly 1% of patients with Parkinson's disease in white populations, 13-30% in Ashkenazi Jewish populations, and 30-40% in North African Arab-Berber populations, although age of onset is variable. Some carriers have early-onset parkinsonism, whereas others remain asymptomatic despite advanced age. We aimed to use a genome-wide approach to identify genetic variability that directly affects LRRK2 Gly2019Ser penetrance., Methods: Between 2006 and 2012, we recruited Arab-Berber patients with Parkinson's disease and their family members (aged 18 years or older) at the Mongi Ben Hamida National Institute of Neurology (Tunis, Tunisia). Patients with Parkinson's disease were diagnosed by movement disorder specialists in accordance with the UK Parkinson's Disease Society Brain Bank criteria, without exclusion of familial parkinsonism. LRRK2 carrier status was confirmed by Sanger sequencing or TaqMan SNP assays-on-demand. We did genome-wide linkage analysis using data from multi-incident Arab-Berber families with Parkinson's disease and LRRK2 Gly2019Ser (with both affected and unaffected family members). We assessed Parkinson's disease age of onset both as a categorical variable (dichotomised by median onset) and as a quantitative trait. We used data from another cohort of unrelated Tunisian LRRK2 Gly2019Ser carriers for subsequent locus-specific genotyping and association analyses. Whole-genome sequencing in a subset of 14 unrelated Arab-Berber individuals who were LRRK2 Gly2019Ser carriers (seven with early-onset disease and seven elderly unaffected individuals) subsequently informed imputation and haplotype analyses. We replicated the findings in separate series of LRRK2 Gly2019Ser carriers originating from Algeria, France, Norway, and North America. We also investigated associations between genotype, gene, and protein expression in human striatal tissues and murine LRRK2 Gly2019Ser cortical neurons., Findings: Using data from 41 multi-incident Arab-Berber families with Parkinson's disease and LRRK2 Gly2019Ser (150 patients and 103 unaffected family members), we identified significant linkage on chromosome 1q23.3 to 1q24.3 (non-parametric logarithm of odds score 2·9, model-based logarithm of odds score 4·99, θ=0 at D1S2768). In a cohort of unrelated Arab-Berber LRRK2 Gly2019Ser carriers, subsequent association mapping within the linkage region suggested genetic variability within DNM3 as an age-of-onset modifier of disease (n=232; rs2421947; haplotype p=1·07 × 10 -7 ). We found that DNM3 rs2421947 was a haplotype tag for which the median onset of LRRK2 parkinsonism in GG carriers was 12·5 years younger than that of CC carriers (Arab-Berber cohort, hazard ratio [HR] 1·89, 95% CI 1·20-2·98). Replication analyses in separate series from Algeria, France, Norway, and North America (n=263) supported this finding (meta-analysis HR 1·61, 95% CI 1·15-2·27, p=0·02). In human striatum, DNM3 expression varied as a function of rs2421947 genotype, and dynamin-3 localisation was perturbed in murine LRRK2 Gly2019Ser cortical neurons., Interpretation: Genetic variability in DNM3 modifies age of onset for LRRK2 Gly2019Ser parkinsonism and informs disease-relevant translational neuroscience. Our results could be useful in genetic counselling for carriers of this mutation and in clinical trial design., Funding: The Canada Excellence Research Chairs (CERC), Leading Edge Endowment Fund (LEEF), Don Rix BC Leadership Chair in Genetic Medicine, National Institute on Aging, National Institute of Neurological Disorders and Stroke, the Michael J Fox Foundation, Mayo Foundation, the Roger de Spoelberch Foundation, and GlaxoSmithKline., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Background: MicroRNAs (miRNAs) are short, non-coding RNAs that regulate gene expression mainly through translational repression of target mRNA molecules. More than 2700 human miRNAs have been identified and some are known to be associated with disease phenotypes and to display tissue-specific patterns of expression., Methods: We used high-throughput small RNA sequencing to discover novel miRNAs in 93 human post-mortem prefrontal cortex samples from individuals with Huntington's disease (n = 28) or Parkinson's disease (n = 29) and controls without neurological impairment (n = 36). A custom miRNA identification analysis pipeline was built, which utilizes miRDeep* miRNA identification and result filtering based on false positive rate estimates., Results: Ninety-nine novel miRNA candidates with a false positive rate of less than 5 % were identified. Thirty-four of the candidate miRNAs show sequence similarity with known mature miRNA sequences and may be novel members of known miRNA families, while the remaining 65 may constitute previously undiscovered families of miRNAs. Nineteen of the 99 candidate miRNAs were replicated using independent, publicly-available human brain RNA-sequencing samples, and seven were experimentally validated using qPCR., Conclusions: We have used small RNA sequencing to identify 99 putative novel miRNAs that are present in human brain samples.

Genome-wide association studies (GWAS) have identified the GAK/DGKQ/IDUA region on 4p16.3 among the top three risk loci for Parkinson's disease (PD), but the specific gene and risk mechanism are unclear. Here, we report transcripts containing the 3' clathrin-binding domain of GAK identified by RNA deep-sequencing in post-mortem human brain tissue as having increased expression in PD. Furthermore, carriers of 4p16.3 PD GWAS risk SNPs show decreased expression of one of these transcripts, GAK25 (Gencode Transcript 009), which correlates with the expression of genes functioning in the synaptic vesicle membrane. Together, these findings provide strong evidence for GAK clathrin-binding- and J-domain transcripts' influence on PD pathogenicity, and for a role for GAK in regulating synaptic function in PD.

Genome-wide association studies (GWAS) have identified numerous genetic variants associated with complex diseases, but mechanistic insights are impeded by a lack of understanding of how specific risk variants functionally contribute to the underlying pathogenesis. It has been proposed that cis-acting effects of non-coding risk variants on gene expression are a major factor for phenotypic variation of complex traits and disease susceptibility. Recent genome-scale epigenetic studies have highlighted the enrichment of GWAS-identified variants in regulatory DNA elements of disease-relevant cell types. Furthermore, single nucleotide polymorphism (SNP)-specific changes in transcription factor binding are correlated with heritable alterations in chromatin state and considered a major mediator of sequence-dependent regulation of gene expression. Here we describe a novel strategy to functionally dissect the cis-acting effect of genetic risk variants in regulatory elements on gene expression by combining genome-wide epigenetic information with clustered regularly-interspaced short palindromic repeats (CRISPR)/Cas9 genome editing in human pluripotent stem cells. By generating a genetically precisely controlled experimental system, we identify a common Parkinson's disease associated risk variant in a non-coding distal enhancer element that regulates the expression of α-synuclein (SNCA), a key gene implicated in the pathogenesis of Parkinson's disease. Our data suggest that the transcriptional deregulation of SNCA is associated with sequence-dependent binding of the brain-specific transcription factors EMX2 and NKX6-1. This work establishes an experimental paradigm to functionally connect genetic variation with disease-relevant phenotypes.

Objective: The goal of this study was to compare the microRNA (miRNA) profile of Parkinson's disease (PD) frontal cortex with normal control brain, allowing for the identification of PD specific signatures as well as study the disease-related phenotypes of onset age and dementia., Methods: Small RNA sequence analysis was performed from prefrontal cortex for 29 PD samples and 33 control samples. After sample QC, normalization and batch correction, linear regression was employed to identify miRNAs altered in PD, and a PD classifier was developed using weighted voting class prediction. The relationship of miRNA levels to onset age and PD with dementia (PDD) was also characterized in case-only analyses., Results: One twenty five miRNAs were differentially expressed in PD at a genome-wide level of significance (FDR q < 0.05). A set of 29 miRNAs classified PD from non-diseased brain (93.9% specificity, 96.6% sensitivity). The majority of differentially expressed miRNAs (105/125) showed an ordinal relationship from control, to PD without dementia (PDN), to PDD. Among PD brains, 36 miRNAs classified PDD from PDN (sensitivity = 81.2%, specificity = 88.9%). Among differentially expressed miRNAs, miR-10b-5p had a positive association with onset age (q = 4.7e-2)., Conclusions: Based on cortical miRNA levels, PD brains were accurately classified from non-diseased brains. Additionally, the PDD miRNA profile exhibited a more severe pattern of alteration among those differentially expressed in PD. To evaluate the clinical utility of miRNAs as potential clinical biomarkers, further characterization and testing of brain-related miRNA alterations in peripheral biofluids is warranted.

Leptin is an adipocyte-secreted hormone, the circulating levels of which correlate closely with overall adiposity. Although rare mutations in the leptin (LEP) gene are well known to cause leptin deficiency and severe obesity, no common loci regulating circulating leptin levels have been uncovered. Therefore, we performed a genome-wide association study (GWAS) of circulating leptin levels from 32,161 individuals and followed up loci reaching P<10(-6) in 19,979 additional individuals. We identify five loci robustly associated (P<5 × 10(-8)) with leptin levels in/near LEP, SLC32A1, GCKR, CCNL1 and FTO. Although the association of the FTO obesity locus with leptin levels is abolished by adjustment for BMI, associations of the four other loci are independent of adiposity. The GCKR locus was found associated with multiple metabolic traits in previous GWAS and the CCNL1 locus with birth weight. Knockdown experiments in mouse adipose tissue explants show convincing evidence for adipogenin, a regulator of adipocyte differentiation, as the novel causal gene in the SLC32A1 locus influencing leptin levels. Our findings provide novel insights into the regulation of leptin production by adipose tissue and open new avenues for examining the influence of variation in leptin levels on adiposity and metabolic health.

Background: Parkinson disease (PD) is a neurodegenerative disease characterized by the accumulation of alpha-synuclein (SNCA) and other proteins in aggregates termed "Lewy Bodies" within neurons. PD has both genetic and environmental risk factors, and while processes leading to aberrant protein aggregation are unknown, past work points to abnormal levels of SNCA and other proteins. Although several genome-wide studies have been performed for PD, these have focused on DNA sequence variants by genome-wide association studies (GWAS) and on RNA levels (microarray transcriptomics), while genome-wide proteomics analysis has been lacking., Methods: This study employed two state-of-the-art technologies, three-stage Mass Spectrometry Tandem Mass Tag Proteomics (12 PD, 12 controls) and RNA-sequencing transcriptomics (29 PD, 44 controls), evaluated in the context of PD GWAS implicated loci and microarray transcriptomics (19 PD, 24 controls). The technologies applied for this study were performed in a set of overlapping prefrontal cortex (Brodmann area 9) samples obtained from PD patients and sex and age similar neurologically healthy controls., Results: After appropriate filters, proteomics robustly identified 3558 unique proteins, with 283 of these (7.9 %) significantly different between PD and controls (q-value < 0.05). RNA-sequencing identified 17,580 protein-coding genes, with 1095 of these (6.2 %) significantly different (FDR p-value < 0.05); only 166 of the FDR significant protein-coding genes (0.94 %) were present among the 3558 proteins characterized. Of these 166, eight genes (4.8 %) were significant in both studies, with the same direction of effect. Functional enrichment analysis of the proteomics results strongly supports mitochondrial-related pathways, while comparable analysis of the RNA-sequencing results implicates protein folding pathways and metallothioneins. Ten of the implicated genes or proteins co-localized to GWAS loci. Evidence implicating SNCA was stronger in proteomics than in RNA-sequencing analyses., Conclusions: We report the largest analysis of proteomics in PD to date, and the first to combine this technology with RNA-sequencing to investigate GWAS implicated loci. Notably, differentially expressed protein-coding genes were more likely to not be characterized in the proteomics analysis, which lessens the ability to compare across platforms. Combining multiple genome-wide platforms offers novel insights into the pathological processes responsible for this disease by identifying pathways implicated across methodologies.

Huntington's disease (HD) is an autosomal-dominant neurodegenerative disorder resulting from expansion of CAG repeats in the Huntingtin (HTT) gene. Previous studies have shown mutant HTT can alter expression of genes associated with dysregulated epigenetic modifications. One of the most widely studied chromatin modifications is trimethylated lysine 4 of histone 3 (H3K4me3). Here, we conducted the first comprehensive study of H3K4me3 ChIP-sequencing in neuronal chromatin from the prefrontal cortex of six HD cases and six non-neurologic controls, and its association with gene expression measured by RNA-sequencing. We detected 2,830 differentially enriched H3K4me3 peaks between HD and controls, with 55% of them down-regulated in HD. Although H3K4me3 signals are expected to be associated with mRNA levels, we found an unexpected discordance between altered H3K4me3 peaks and mRNA levels. Gene ontology (GO) term enrichment analysis of the genes with differential H3K4me3 peaks, revealed statistically significantly enriched GO terms only in the genes with down-regulated signals in HD. The most frequently implicated biological process terms are organ morphogenesis and positive regulation of gene expression. More than 9,000 H3K4me3 peaks were located not near any recognized transcription start sites and approximately 36% of these "distal" peaks co-localized to known enhancer sites. Six transcription factors and chromatin remodelers are differentially enriched in HD H3K4me3 distal peaks, including EZH2 and SUZ12, two core subunits of the polycomb repressive complex 2 (PRC2). Moreover, PRC2 repressive state was significantly depleted in HD-enriched peaks, suggesting the epigenetic role of PRC2 inhibition associated with up-regulated H3K4me3 in Huntington's disease. In summary, our study provides new insights into transcriptional dysregulation of Huntington's disease by analyzing the differentiation of H3K4me3 enrichment.

Huntington's Disease (HD) is a devastating neurodegenerative disorder that is caused by an expanded CAG trinucleotide repeat in the Huntingtin (HTT) gene. Transcriptional dysregulation in the human HD brain has been documented but is incompletely understood. Here we present a genome-wide analysis of mRNA expression in human prefrontal cortex from 20 HD and 49 neuropathologically normal controls using next generation high-throughput sequencing. Surprisingly, 19% (5,480) of the 28,087 confidently detected genes are differentially expressed (FDR<0.05) and are predominantly up-regulated. A novel hypothesis-free geneset enrichment method that dissects large gene lists into functionally and transcriptionally related groups discovers that the differentially expressed genes are enriched for immune response, neuroinflammation, and developmental genes. Markers for all major brain cell types are observed, suggesting that HD invokes a systemic response in the brain area studied. Unexpectedly, the most strongly differentially expressed genes are a homeotic gene set (represented by Hox and other homeobox genes), that are almost exclusively expressed in HD, a profile not widely implicated in HD pathogenesis. The significance of transcriptional changes of developmental processes in the HD brain is poorly understood and warrants further investigation. The role of inflammation and the significance of non-neuronal involvement in HD pathogenesis suggest anti-inflammatory therapeutics may offer important opportunities in treating HD.