The Parkinson's disease-associated gene ITPKB protects against α-synuclein aggregation by regulating ER-to-mitochondria calcium release.

Inositol-1,4,5-triphosphate (IP3) kinase B (ITPKB) is a ubiquitously expressed lipid kinase that inactivates IP3, a secondary messenger that stimulates calcium release from the endoplasmic reticulum (ER). Genome-wide association studies have identified common variants in the ITPKB gene locus associated with reduced risk of sporadic Parkinson's disease (PD). Here, we investigate whether ITPKB activity or expression level impacts PD phenotypes in cellular and animal models. In primary neurons, knockdown or pharmacological inhibition of ITPKB increased levels of phosphorylated, insoluble α-synuclein pathology following treatment with α-synuclein preformed fibrils (PFFs). Conversely, ITPKB overexpression reduced PFF-induced α-synuclein aggregation. We also demonstrate that ITPKB inhibition or knockdown increases intracellular calcium levels in neurons, leading to an accumulation of calcium in mitochondria that increases respiration and inhibits the initiation of autophagy, suggesting that ITPKB regulates α-synuclein pathology by inhibiting ER-to-mitochondria calcium transport. Furthermore, the effects of ITPKB on mitochondrial calcium and respiration were prevented by pretreatment with pharmacological inhibitors of the mitochondrial calcium uniporter complex, which was also sufficient to reduce α-synuclein pathology in PFF-treated neurons. Taken together, these results identify ITPKB as a negative regulator of α-synuclein aggregation and highlight modulation of ER-to-mitochondria calcium flux as a therapeutic strategy for the treatment of sporadic PD.

Differences in the Presentation and Progression of Parkinson's Disease by Sex.

BACKGROUND: Previous studies reported various symptoms of Parkinson's disease (PD) associated with sex. Some were conflicting or confirmed in only one study. OBJECTIVES: We examined sex associations to PD phenotypes cross-sectionally and longitudinally in large-scale data. METHODS: We tested 40 clinical phenotypes, using longitudinal, clinic-based patient cohorts, consisting of 5946 patients, with a median follow-up of 3.1 years. For continuous outcomes, we used linear regressions at baseline to test sex-associated differences in presentation, and linear mixed-effects models to test sex-associated differences in progression. For binomial outcomes, we used logistic regression models at baseline and Cox regression models for survival analyses. We adjusted for age, disease duration, and medication use. In the secondary analyses, data from 17 719 PD patients and 7588 non-PD participants from an online-only, self-assessment PD cohort were cross-sectionally evaluated to determine whether the sex-associated differences identified in the primary analyses were consistent and unique to PD. RESULTS: Female PD patients had a higher risk of developing dyskinesia early during the follow-up period, with a slower progression in activities of daily living difficulties, and a lower risk of developing cognitive impairments compared with male patients. The findings in the longitudinal, clinic-based cohorts were mostly consistent with the results of the online-only cohort. CONCLUSIONS: We observed sex-associated contributions to PD heterogeneity. These results highlight the necessity of future research to determine the underlying mechanisms and importance of personalized clinical management. © 2020 International Parkinson and Movement Disorder Society.

Genomewide association study of Parkinson's disease clinical biomarkers in 12 longitudinal patients' cohorts.

BACKGROUND: Several reports have identified different patterns of Parkinson's disease progression in individuals carrying missense variants in GBA or LRRK2 genes. The overall contribution of genetic factors to the severity and progression of Parkinson's disease, however, has not been well studied. OBJECTIVES: To test the association between genetic variants and the clinical features of Parkinson's disease on a genomewide scale. METHODS: We accumulated individual data from 12 longitudinal cohorts in a total of 4093 patients with 22,307 observations for a median of 3.81 years. Genomewide associations were evaluated for 25 cross-sectional and longitudinal phenotypes. Specific variants of interest, including 90 recently identified disease-risk variants, were also investigated post hoc for candidate associations with these phenotypes. RESULTS: Two variants were genomewide significant. Rs382940(T>A), within the intron of SLC44A1, was associated with reaching Hoehn and Yahr stage 3 or higher faster (hazard ratio 2.04 [1.58-2.62]; P value = 3.46E-8). Rs61863020(G>A), an intergenic variant and expression quantitative trait loci for α-2A adrenergic receptor, was associated with a lower prevalence of insomnia at baseline (odds ratio 0.63 [0.52-0.75]; P value = 4.74E-8). In the targeted analysis, we found 9 associations between known Parkinson's risk variants and more severe motor/cognitive symptoms. Also, we replicated previous reports of GBA coding variants (rs2230288: p.E365K; rs75548401: p.T408M) being associated with greater motor and cognitive decline over time, and an APOE E4 tagging variant (rs429358) being associated with greater cognitive deficits in patients. CONCLUSIONS: We identified novel genetic factors associated with heterogeneity of Parkinson's disease. The results can be used for validation or hypothesis tests regarding Parkinson's disease. © 2019 International Parkinson and Movement Disorder Society.

Genetic risk of Parkinson disease and progression:: An analysis of 13 longitudinal cohorts.

OBJECTIVE: To determine if any association between previously identified alleles that confer risk for Parkinson disease and variables measuring disease progression. METHODS: We evaluated the association between 31 risk variants and variables measuring disease progression. A total of 23,423 visits by 4,307 patients of European ancestry from 13 longitudinal cohorts in Europe, North America, and Australia were analyzed. RESULTS: We confirmed the importance of GBA on phenotypes. GBA variants were associated with the development of daytime sleepiness (p.N370S: hazard ratio [HR] 3.28 [1.69-6.34]) and possible REM sleep behavior (p.T408M: odds ratio 6.48 [2.04-20.60]). We also replicated previously reported associations of GBA variants with motor/cognitive declines. The other genotype-phenotype associations include an intergenic variant near LRRK2 and the faster development of motor symptom (Hoehn and Yahr scale 3.0 HR 1.33 [1.16-1.52] for the C allele of rs76904798) and an intronic variant in PMVK and the development of wearing-off effects (HR 1.66 [1.19-2.31] for the C allele of rs114138760). Age at onset was associated with TMEM175 variant p.M393T (-0.72 [-1.21 to -0.23] in years), the C allele of rs199347 (intronic region of GPNMB, 0.70 [0.27-1.14]), and G allele of rs1106180 (intronic region of CCDC62, 0.62 [0.21-1.03]). CONCLUSIONS: This study provides evidence that alleles associated with Parkinson disease risk, in particular GBA variants, also contribute to the heterogeneity of multiple motor and nonmotor aspects. Accounting for genetic variability will be a useful factor in understanding disease course and in minimizing heterogeneity in clinical trials.

Molecular Genetic Anatomy and Risk Profile of Hirschsprung's Disease.

BACKGROUND: Hirschsprung's disease, or congenital aganglionosis, is a developmental disorder of the enteric nervous system and is the most common cause of intestinal obstruction in neonates and infants. The disease has more than 80% heritability, including significant associations with rare and common sequence variants in genes related to the enteric nervous system, as well as with monogenic and chromosomal syndromes. METHODS: We genotyped and exome-sequenced samples from 190 patients with Hirschsprung's disease to quantify the genetic burden in patients with this condition. DNA sequence variants, large copy-number variants, and karyotype variants in probands were considered to be pathogenic when they were significantly associated with Hirschsprung's disease or another neurodevelopmental disorder. Novel genes were confirmed by functional studies in the mouse and human embryonic gut and in zebrafish embryos. RESULTS: The presence of five or more variants in four noncoding elements defined a widespread risk of Hirschsprung's disease (48.4% of patients and 17.1% of controls; odds ratio, 4.54; 95% confidence interval [CI], 3.19 to 6.46). Rare coding variants in 24 genes that play roles in enteric neural-crest cell fate, 7 of which were novel, were also common (34.7% of patients and 5.0% of controls) and conferred a much greater risk than noncoding variants (odds ratio, 10.02; 95% CI, 6.45 to 15.58). Large copy-number variants, which were present in fewer patients (11.4%, as compared with 0.2% of controls), conferred the highest risk (odds ratio, 63.07; 95% CI, 36.75 to 108.25). At least one identifiable genetic risk factor was found in 72.1% of the patients, and at least 48.4% of patients had a structural or regulatory deficiency in the gene encoding receptor tyrosine kinase (RET). For individual patients, the estimated risk of Hirschsprung's disease ranged from 5.33 cases per 100,000 live births (approximately 1 per 18,800) to 8.38 per 1000 live births (approximately 1 per 120). CONCLUSIONS: Among the patients in our study, Hirschsprung's disease arose from common noncoding variants, rare coding variants, and copy-number variants affecting genes involved in enteric neural-crest cell fate that exacerbate the widespread genetic susceptibility associated with RET. For individual patients, the genotype-specific odds ratios varied by a factor of approximately 67, which provides a basis for risk stratification and genetic counseling. (Funded by the National Institutes of Health.).

Phenome-wide association studies across large population cohorts support drug target validation.

Phenome-wide association studies (PheWAS) have been proposed as a possible aid in drug development through elucidating mechanisms of action, identifying alternative indications, or predicting adverse drug events (ADEs). Here, we select 25 single nucleotide polymorphisms (SNPs) linked through genome-wide association studies (GWAS) to 19 candidate drug targets for common disease indications. We interrogate these SNPs by PheWAS in four large cohorts with extensive health information (23andMe, UK Biobank, FINRISK, CHOP) for association with 1683 binary endpoints in up to 697,815 individuals and conduct meta-analyses for 145 mapped disease endpoints. Our analyses replicate 75% of known GWAS associations (P < 0.05) and identify nine study-wide significant novel associations (of 71 with FDR < 0.1). We describe associations that may predict ADEs, e.g., acne, high cholesterol, gout, and gallstones with rs738409 (p.I148M) in PNPLA3 and asthma with rs1990760 (p.T946A) in IFIH1. Our results demonstrate PheWAS as a powerful addition to the toolkit for drug discovery.

LRRK2 p.Ile1371Val Mutation in a Case with Neuropathologically Confirmed Multi-System Atrophy.

BACKGROUND: Mutations in the leucine rich repeat kinase 2 (LRRK2) gene are among the most common genetic causes of Lewy body Parkinson's disease (PD). However, LRRK2 mutations can also lead to a variety of pathological phenotypes other than typical PD, including relatively pure nigrostriatal cell loss without alpha-synuclein-positive Lewy bodies or Lewy neurites, progressive supranuclear palsy (PSP), and multiple system atrophy (MSA). The mechanisms behind this remarkable pleomorphic pathology are currently unclear. OBJECTIVE: To genetically and pathologically characterize a case with a LRRK2, p.Ile1371Val rare variant and pathologically proven MSA. METHODS: From the brain donation program at the Parkinson's Institute and Clinical Center, we selected 26 brains with family history and a with clinicopathological diagnosis of PD (n = 20), MSA (n = 4), or PSP (n = 2). We performed neuropathological evaluation, including alpha-synuclein and tau immunohistochemistry and sequenced 188 genes that have been reported as causative for or associated with neurodegenerative diseases. RESULTS: We identified a known LRRK2, p.Ile1371Val genetic variant in a case with clinically diagnosed and pathologically proven MSA. Neuropathology revealed that the olivopontocerebellar system was more affected than the striatonigral system. CONCLUSIONS: Our data suggest that genetic variants in the LRRK2 gene can present clinically and neuropathologically as MSA. One other LRRK2 genetic variant (LRRK2, p.Ile2020Thr) has been reported with a neuropathological diagnosis of MSA. Interestingly, the LRRK2 variant (LRRK2, p.Ile1371Val) identified here has been reported previously in a postmortem case with Lewy body PD.Future studies are critical to discover the mechanisms leading to different neurodegenerative trajectories both in neuronal and glial cell populations.

Parkinson's disease associated with pure ATXN10 repeat expansion.

Large, non-coding pentanucleotide repeat expansions of ATTCT in intron 9 of the ATXN10 gene typically cause progressive spinocerebellar ataxia with or without seizures and present neuropathologically with Purkinje cell loss resulting in symmetrical cerebellar atrophy. These ATXN10 repeat expansions can be interrupted by sequence motifs which have been attributed to seizures and are likely to act as genetic modifiers. We identified a Mexican kindred with multiple affected family members with ATXN10 expansions. Four affected family members showed clinical features of spinocerebellar ataxia type 10 (SCA10). However, one affected individual presented with early-onset levodopa-responsive parkinsonism, and one family member carried a large repeat ATXN10 expansion, but was clinically unaffected. To characterize the ATXN10 repeat, we used a novel technology of single-molecule real-time (SMRT) sequencing and CRISPR/Cas9-based capture. We sequenced the entire span of ~5.3-7.0 kb repeat expansions. The Parkinson's patient carried an ATXN10 expansion with no repeat interruption motifs as well as an unaffected sister. In the siblings with typical SCA10, we found a repeat pattern of ATTCC repeat motifs that have not been associated with seizures previously. Our data suggest that the absence of repeat interruptions is likely a genetic modifier for the clinical presentation of l-Dopa responsive parkinsonism, whereas repeat interruption motifs contribute clinically to epilepsy. Repeat interruptions are important genetic modifiers of the clinical phenotype in SCA10. Advanced sequencing techniques now allow to better characterize the underlying genetic architecture for determining accurate phenotype-genotype correlations.

Trans-ethnic fine-mapping of genetic loci for body mass index in the diverse ancestral populations of the Population Architecture using Genomics and Epidemiology (PAGE) Study reveals evidence for multiple signals at established loci.

Most body mass index (BMI) genetic loci have been identified in studies of primarily European ancestries. The effect of these loci in other racial/ethnic groups is less clear. Thus, we aimed to characterize the generalizability of 170 established BMI variants, or their proxies, to diverse US populations and trans-ethnically fine-map 36 BMI loci using a sample of >102,000 adults of African, Hispanic/Latino, Asian, European and American Indian/Alaskan Native descent from the Population Architecture using Genomics and Epidemiology Study. We performed linear regression of the natural log of BMI (18.5-70 kg/m2) on the additive single nucleotide polymorphisms (SNPs) at BMI loci on the MetaboChip (Illumina, Inc.), adjusting for age, sex, population stratification, study site, or relatedness. We then performed fixed-effect meta-analyses and a Bayesian trans-ethnic meta-analysis to empirically cluster by allele frequency differences. Finally, we approximated conditional and joint associations to test for the presence of secondary signals. We noted directional consistency with the previously reported risk alleles beyond what would have been expected by chance (binomial p < 0.05). Nearly, a quarter of the previously described BMI index SNPs and 29 of 36 densely-genotyped BMI loci on the MetaboChip replicated/generalized in trans-ethnic analyses. We observed multiple signals at nine loci, including the description of seven loci with novel multiple signals. This study supports the generalization of most common genetic loci to diverse ancestral populations and emphasizes the importance of dense multiethnic genomic data in refining the functional variation at genetic loci of interest and describing several loci with multiple underlying genetic variants.

The role of rare variants in systolic blood pressure: analysis of ExomeChip data in HyperGEN African Americans.

Cardiovascular diseases are among the most significant health problems in the United States today, with their major risk factor, hypertension, disproportionately affecting African Americans (AAs). Although GWAS have identified dozens of common variants associated with blood pressure (BP) and hypertension in European Americans, these variants collectively explain <2.5% of BP variance, and most of the genetic variants remain yet to be identified. Here, we report the results from rare-variant analysis of systolic BP using 94,595 rare and low-frequency variants (minor allele frequency, MAF, <5%) from the Illumina exome array genotyped in 2,045 HyperGEN AAs. In addition to single-variant analysis, 4 gene-level association tests were used for analysis: burden and family-based SKAT tests using MAF cutoffs of 1 and 5%. The gene-based methods often provided lower p values than the single-variant approach. Some consistency was observed across these 4 gene-based analysis options. While neither the gene-based analyses nor the single-variant analysis produced genome-wide significant results, the top signals, which had supporting evidence from multiple gene-based methods, were of borderline significance. Though additional molecular validations are required, 6 of the 16 most promising genes are biologically plausible with physiological connections to BP regulation.

Expression quantitative trait Loci acting across multiple tissues are enriched in inherited risk for coronary artery disease.

BACKGROUND: Despite recent discoveries of new genetic risk factors, the majority of risk for coronary artery disease (CAD) remains elusive. As the most proximal sensor of DNA variation, RNA abundance can help identify subpopulations of genetic variants active in and across tissues mediating CAD risk through gene expression. METHODS AND RESULTS: By generating new genomic data on DNA and RNA samples from the Stockholm Atherosclerosis Gene Expression (STAGE) study, 8156 cis-acting expression quantitative trait loci (eQTLs) for 6450 genes across 7 CAD-relevant tissues were detected. The inherited risk enrichments of tissue-defined sets of these eQTLs were assessed using 2 independent genome-wide association data sets. eQTLs acting across increasing numbers of tissues were found increasingly enriched for CAD risk and resided at regulatory hot spots. The risk enrichment of 42 eQTLs acting across 5 to 6 tissues was particularly high (≤7.3-fold) and confirmed in the combined genome-wide association data from Coronary Artery Disease Genome Wide Replication And Meta-Analysis Consortium. Sixteen of the 42 eQTLs associated with 19 master regulatory genes and 29 downstream gene sets (n>30) were further risk enriched comparable to that of the 153 genome-wide association risk single-nucleotide polymorphisms established for CAD (8.4-fold versus 10-fold). Three gene sets, governed by the master regulators FLYWCH1, PSORSIC3, and G3BP1, segregated the STAGE patients according to extent of CAD, and small interfering RNA targeting of these master regulators affected cholesterol-ester accumulation in foam cells of the THP1 monocytic cell line. CONCLUSIONS: eQTLs acting across multiple tissues are significant carriers of inherited risk for CAD. FLYWCH1, PSORSIC3, and G3BP1 are novel master regulatory genes in CAD that may be suitable targets.

Fine Mapping and Identification of BMI Loci in African Americans.

Genome-wide association studies (GWASs) primarily performed in European-ancestry (EA) populations have identified numerous loci associated with body mass index (BMI). However, it is still unclear whether these GWAS loci can be generalized to other ethnic groups, such as African Americans (AAs). Furthermore, the putative functional variant or variants in these loci mostly remain under investigation. The overall lower linkage disequilibrium in AA compared to EA populations provides the opportunity to narrow in or fine-map these BMI-related loci. Therefore, we used the Metabochip to densely genotype and evaluate 21 BMI GWAS loci identified in EA studies in 29,151 AAs from the Population Architecture using Genomics and Epidemiology (PAGE) study. Eight of the 21 loci (SEC16B, TMEM18, ETV5, GNPDA2, TFAP2B, BDNF, FTO, and MC4R) were found to be associated with BMI in AAs at 5.8 × 10(-5). Within seven out of these eight loci, we found that, on average, a substantially smaller number of variants was correlated (r(2) > 0.5) with the most significant SNP in AA than in EA populations (16 versus 55). Conditional analyses revealed GNPDA2 harboring a potential additional independent signal. Moreover, Metabochip-wide discovery analyses revealed two BMI-related loci, BRE (rs116612809, p = 3.6 × 10(-8)) and DHX34 (rs4802349, p = 1.2 × 10(-7)), which were significant when adjustment was made for the total number of SNPs tested across the chip. These results demonstrate that fine mapping in AAs is a powerful approach for both narrowing in on the underlying causal variants in known loci and discovering BMI-related loci.

A systematic mapping approach of 16q12.2/FTO and BMI in more than 20,000 African Americans narrows in on the underlying functional variation: results from the Population Architecture using Genomics and Epidemiology (PAGE) study.

Genetic variants in intron 1 of the fat mass- and obesity-associated (FTO) gene have been consistently associated with body mass index (BMI) in Europeans. However, follow-up studies in African Americans (AA) have shown no support for some of the most consistently BMI-associated FTO index single nucleotide polymorphisms (SNPs). This is most likely explained by different race-specific linkage disequilibrium (LD) patterns and lower correlation overall in AA, which provides the opportunity to fine-map this region and narrow in on the functional variant. To comprehensively explore the 16q12.2/FTO locus and to search for second independent signals in the broader region, we fine-mapped a 646-kb region, encompassing the large FTO gene and the flanking gene RPGRIP1L by investigating a total of 3,756 variants (1,529 genotyped and 2,227 imputed variants) in 20,488 AAs across five studies. We observed associations between BMI and variants in the known FTO intron 1 locus: the SNP with the most significant p-value, rs56137030 (8.3 × 10(-6)) had not been highlighted in previous studies. While rs56137030was correlated at r(2)>0.5 with 103 SNPs in Europeans (including the GWAS index SNPs), this number was reduced to 28 SNPs in AA. Among rs56137030 and the 28 correlated SNPs, six were located within candidate intronic regulatory elements, including rs1421085, for which we predicted allele-specific binding affinity for the transcription factor CUX1, which has recently been implicated in the regulation of FTO. We did not find strong evidence for a second independent signal in the broader region. In summary, this large fine-mapping study in AA has substantially reduced the number of common alleles that are likely to be functional candidates of the known FTO locus. Importantly our study demonstrated that comprehensive fine-mapping in AA provides a powerful approach to narrow in on the functional candidate(s) underlying the initial GWAS findings in European populations.

Effects of rare and common blood pressure gene variants on essential hypertension: results from the Family Blood Pressure Program, CLUE, and Atherosclerosis Risk in Communities studies.

RATIONALE: Hypertension affects ≈30% of adults in industrialized countries and is the major risk factor for cardiovascular disease. OBJECTIVE: We sought to study the genetic effect of coding and conserved noncoding variants in syndromic hypertension genes on systolic blood pressure (BP) and diastolic BP to assess their overall impact on essential hypertension. METHODS AND RESULTS: We resequenced 11 genes (AGT, CYP11B1, CYP17A1, HSD11B2, NR3C1, NR3C2, SCNN1A, SCNN1B, SCNN1G, WNK1, and WNK4) in 560 European American (EA) and African American ancestry GenNet participants with extreme systolic BP. We investigated genetic associations of 2535 variants with BP in 19997 EAs and in 6069 African Americans in 3 types of analyses. First, we studied the combined effects of all variants in GenNet. Second, we studied 1000 Genomes imputed polymorphic variants in 9747 EA and 3207 African American Atherosclerosis Risk in Communities subjects. Finally, we genotyped 37 missense and common noncoding variants in 6591 EAs and in 6521 individuals (3659 EA/2862 African American) from the CLUE and Family Blood Pressure Program studies, respectively. None of the variants individually reached significant false-discovery rates ≤0.05 for systolic BP and diastolic BP. However, on pooling all coding and noncoding variants, we identified at least 5 loci (AGT, CYP11B1, NR3C2, SCNN1G, and WNK1) with higher association at evolutionary conserved sites. CONCLUSIONS: Both rare and common variants at these genes affect BP in the general population with modest effects sizes (<0.05 standard deviation units), and much larger sample sizes are required to assess the impact of individual genes. Collectively, conserved noncoding variants affect BP to a greater extent than missense mutations.

Upper temperature limits of tropical marine ectotherms: global warming implications.

Animal physiology, ecology and evolution are affected by temperature and it is expected that community structure will be strongly influenced by global warming. This is particularly relevant in the tropics, where organisms are already living close to their upper temperature limits and hence are highly vulnerable to rising temperature. Here we present data on upper temperature limits of 34 tropical marine ectotherm species from seven phyla living in intertidal and subtidal habitats. Short term thermal tolerances and vertical distributions were correlated, i.e., upper shore animals have higher thermal tolerance than lower shore and subtidal animals; however, animals, despite their respective tidal height, were susceptible to the same temperature in the long term. When temperatures were raised by 1°C hour(-1), the upper lethal temperature range of intertidal ectotherms was 41-52°C, but this range was narrower and reduced to 37-41°C in subtidal animals. The rate of temperature change, however, affected intertidal and subtidal animals differently. In chronic heating experiments when temperature was raised weekly or monthly instead of every hour, upper temperature limits of subtidal species decreased from 40°C to 35.4°C, while the decrease was more than 10°C in high shore organisms. Hence in the long term, activity and survival of tropical marine organisms could be compromised just 2-3°C above present seawater temperatures. Differences between animals from environments that experience different levels of temperature variability suggest that the physiological mechanisms underlying thermal sensitivity may vary at different rates of warming.