Background There is mounting evidence for a connection between the gut and Parkinson’s disease (PD). Dysbiosis of gut microbiota could explain several features of PD. Objective To determine if PD involves dysbiosis of gut microbiome, disentangle effects of confounders, and identify candidate taxa and functional pathways to guide research. Methods 197 PD cases and 130 controls were studied. Microbial composition was determined by 16S rRNA gene sequencing of DNA extracted from stool. Metadata were collected on 39 potential confounders including medications, diet, gastrointestinal symptoms, and demographics. Statistical analyses were conducted while controlling for potential confounders and correcting for multiple testing. We tested differences in the overall microbial composition, taxa abundance, and functional pathways. Results Independent microbial signatures were detected for PD (P=4E-5), subjects’ region of residence within the United States (P=3E-3), age (P=0.03), sex (P=1E-3) and dietary fruits/vegetables (P=0.01). Among patients, independent signals were detected for catechol-O-methyltransferase-inhibitors (P=4E-4), anticholinergics (P=5E-3), and possibly carbidopa/levodopa (P=0.05). We found significantly altered abundance of Bifidobacteriaceae, Christensenellaceae, [Tissierellaceae], Lachnospiraceae, Lactobacillaceae, Pasteurellaceae and Verrucomicrobiaceae families. Functional predictions revealed changes in numerous pathways including metabolism of plant-derived compounds and xenobiotics degradation. Conclusion PD is accompanied by dysbiosis of gut microbiome. Results coalesce divergent findings of prior studies, reveal altered abundance of several taxa, nominate functional pathways, and demonstrate independent effects of PD medications on the microbiome. The findings provide new leads and testable hypotheses on the pathophysiology and treatment of PD.

Our aim was to identify genes that influence the inverse association of coffee with the risk of developing Parkinson's disease (PD). We used genome-wide genotype data and lifetime caffeinated-coffee-consumption data on 1,458 persons with PD and 931 without PD from the NeuroGenetics Research Consortium (NGRC), and we performed a genome-wide association and interaction study (GWAIS), testing each SNP's main-effect plus its interaction with coffee, adjusting for sex, age, and two principal components. We then stratified subjects as heavy or light coffee-drinkers and performed genome-wide association study (GWAS) in each group. We replicated the most significant SNP. Finally, we imputed the NGRC dataset, increasing genomic coverage to examine the region of interest in detail. The primary analyses (GWAIS, GWAS, Replication) were performed using genotyped data. In GWAIS, the most significant signal came from rs4998386 and the neighboring SNPs in GRIN2A. GRIN2A encodes an NMDA-glutamate-receptor subunit and regulates excitatory neurotransmission in the brain. Achieving P2df = 10−6, GRIN2A surpassed all known PD susceptibility genes in significance in the GWAIS. In stratified GWAS, the GRIN2A signal was present in heavy coffee-drinkers (OR = 0.43; P = 6×10−7) but not in light coffee-drinkers. The a priori Replication hypothesis that “Among heavy coffee-drinkers, rs4998386_T carriers have lower PD risk than rs4998386_CC carriers” was confirmed: ORReplication = 0.59, PReplication = 10−3; ORPooled = 0.51, PPooled = 7×10−8. Compared to light coffee-drinkers with rs4998386_CC genotype, heavy coffee-drinkers with rs4998386_CC genotype had 18% lower risk (P = 3×10−3), whereas heavy coffee-drinkers with rs4998386_TC genotype had 59% lower risk (P = 6×10−13). Imputation revealed a block of SNPs that achieved P2df<5×10−8 in GWAIS, and OR = 0.41, P = 3×10−8 in heavy coffee-drinkers. This study is proof of concept that inclusion of environmental factors can help identify genes that are missed in GWAS. Both adenosine antagonists (caffeine-like) and glutamate antagonists (GRIN2A-related) are being tested in clinical trials for treatment of PD. GRIN2A may be a useful pharmacogenetic marker for subdividing individuals in clinical trials to determine which medications might work best for which patients.

Historically, association of disease with the major histocompatibility complex (HLA) genes has been tested with HLA alleles that encode antigen-binding affinity. The association with Parkinson disease (PD), however, was discovered with noncoding SNPs in a genome-wide association study (GWAS). We show here that several HLA-region SNPs that have since been associated with PD form two blocks tagged by rs3129882 (p = 9 × 10(-11)) and by rs9268515 and/or rs2395163 (p = 3 × 10(-11)). We investigated whether these SNP-associations were driven by HLA-alleles at adjacent loci. We imputed class I and class II HLA-alleles for 2000 PD cases and 1986 controls from the NeuroGenetics Research Consortium GWAS and sequenced a subset of 194 cases and 204 controls. We were therefore able to assess accuracy of two imputation algorithms against next-generation-sequencing while taking advantage of the larger imputed data sets for disease study. Additionally, we imputed HLA alleles for 843 cases and 856 controls from another GWAS for replication. PD risk was positively associated with the B(∗)07:02_C(∗)07:02_DRB5(∗)01_DRB1(∗)15:01_DQA1(∗)01:02_DQB1(∗)06:02 haplotype and negatively associated with the C(∗)03:04, DRB1(∗)04:04 and DQA1(∗)03:01 alleles. The risk haplotype and DQA1(∗)03:01 lost significance when conditioned on the SNPs, but C(∗)03:04 (OR = 0.72, p = 8 × 10(-6)) and DRB1(∗)04:04 (OR = 0.65, p = 4 × 10(-5)) remained significant. Similarly, rs3129882 and the closely linked rs9268515 and rs2395163 remained significant irrespective of HLA alleles. rs3129882 and rs2395163 are expression quantitative trait loci (eQTLs) for HLA-DR and HLA-DQ (9 × 10(-5) ≥ PeQTL ≥ 2 × 10(-79)), suggesting that HLA gene expression might influence PD. Our data suggest that PD is associated with both structural and regulatory elements in HLA. Furthermore, our study demonstrates that noncoding SNPs in the HLA region can be associated with disease irrespective of HLA alleles, and that observed associations with HLA alleles can sometimes be secondary to a noncoding variant.

Levels of stool immune factors indicate that intestinal inflammation is present in patients with Parkinson's disease. © 2018 International Parkinson and Movement Disorder Society.

Prior studies have established an inverse association between cigarette-smoking and the risk of developing Parkinson's disease (PD), and currently, disease-modifying potential of the nicotine-patch is being tested in clinical trials. To identify genes that interact with the effect of smoking/nicotine, we conducted genome-wide interaction studies in humans and in Drosophila. We identified SV2C which encodes a synaptic-vesicle protein in PD-vulnerable substantia-nigra (P=1×10-7 for gene-smoking interaction on PD risk), and CG14691 which is predicted to encode a synaptic-vesicle protein in Drosophila (P=2×10-11 for nicotine-paraquat interaction on gene-expression). SV2C is biologically plausible because nicotine enhances release of dopamine through synaptic vesicles, and PD is caused by depletion of dopamine. Effect of smoking on PD varied by SV2C genotype from protective to neutral to harmful (P=5×10-10). Taken together, cross-validating evidence from humans and Drosophila suggest SV2C is involved in PD pathogenesis and it might be a useful marker for pharmacogenomics studies involving nicotine.

This study quantifies the effects of naturally occurring X-linked variation on immune response in Drosophila melanogaster to assess associations between immunity genotypes and innate immune response. We constructed a set of 168 X-chromosomal extraction lines, incorporating X chromosomes from a natural population into co-isogenic autosomal backgrounds, and genotyped the lines at 88 SNPs in 20 X-linked immune genes. We find that genetic variation in many of the genes is associated with immune response phenotypes, including bacterial load and immune gene expression. Many of the associations act in a sex-specific or sexually antagonistic manner, supporting the theory that with the selective pressures facing genes on the X chromosome, sexually antagonistic variation may be more easily maintained.

Parkinson's disease (PD) was recently found to be associated with HLA in a genome-wide association study (GWAS). Follow-up GWAS's replicated the PD-HLA association but their top hits differ. Do the different hits tag the same locus or is there more than one PD-associated variant within HLA? We show that the top GWAS hits are not correlated with each other (0.00≤r2≤0.15). Using our GWAS (2000 cases, 1986 controls) we conducted step-wise conditional analysis on 107 SNPs with P<10−3 for PD-association; 103 dropped-out, four remained significant. Each SNP, when conditioned on the other three, yielded PSNP1 = 5×10−4, PSNP2 = 5×10−4, PSNP3 = 4×10−3 and PSNP4 = 0.025. The four SNPs were not correlated (0.01≤r2≤0.20). Haplotype analysis (excluding rare SNP2) revealed increasing PD risk with increasing risk alleles from OR = 1.27, P = 5×10−3 for one risk allele to OR = 1.65, P = 4×10−8 for three. Using additional 843 cases and 856 controls we replicated the independent effects of SNP1 (Pconditioned-on-SNP4 = 0.04) and SNP4 (Pconditioned-on-SNP1 = 0.04); SNP2 and SNP3 could not be replicated. In pooled GWAS and replication, SNP1 had ORconditioned-on-SNP4 = 1.23, Pconditioned-on-SNP4 = 6×10−7; SNP4 had ORconditioned-on-SNP1 = 1.18, Pconditioned-on-SNP1 = 3×10−3; and the haplotype with both risk alleles had OR = 1.48, P = 2×10−12. Genotypic OR increased with the number of risk alleles an individual possessed up to OR = 1.94, P = 2×10−11 for individuals who were homozygous for the risk allele at both SNP1 and SNP4. SNP1 is a variant in HLA-DRA and is associated with HLA-DRA, DRB5 and DQA2 gene expression. SNP4 is correlated (r2 = 0.95) with variants that are associated with HLA-DQA2 expression, and with the top HLA SNP from the IPDGC GWAS (r2 = 0.60). Our findings suggest more than one PD-HLA association; either different alleles of the same gene, or separate loci.

Parkinson’s disease (PD) is the most common cause of neurodegenerative movement disorder and the second most common cause of dementia. Genes are thought to have a stronger effect on age-at-onset of PD than on risk, yet there has been a phenomenal success in identifying risk loci but not age-at-onset modifiers. We conducted a genome-wide study for age-at-onset. We analysed familial and non-familial PD separately, per prior evidence for strong genetic effect on age-at-onset in familial PD. GWAS was conducted in 431 unrelated PD individuals with at least one affected relative (familial PD) and 1544 non-familial PD from the NeuroGenetics Research Consortium (NGRC); an additional 737 familial PD and 2363 non-familial PD were used for replication. In familial PD, two signals were detected and replicated robustly: one mapped to LHFPL2 on 5q14.1 (PNGRC = 3E-8, PReplication = 2E-5, PNGRC + Replication = 1E-11), the second mapped to TPM1 on 15q22.2 (PNGRC = 8E-9, PReplication = 2E-4, PNGRC + Replication = 9E-11). The variants that were associated with accelerated onset had low frequencies (<0.02). The LHFPL2 variant was associated with earlier onset by 12.33 [95% CI: 6.2; 18.45] years in NGRC, 8.03 [2.95; 13.11] years in replication, and 9.79 [5.88; 13.70] years in the combined data. The TPM1 variant was associated with earlier onset by 15.30 [8.10; 22.49] years in NGRC, 9.29 [1.79; 16.79] years in replication, and 12.42 [7.23; 17.61] years in the combined data. Neither LHFPL2 nor TPM1 was associated with age-at-onset in non-familial PD. LHFPL2 (function unknown) is overexpressed in brain tumours. TPM1 encodes a highly conserved protein that regulates muscle contraction, and is a tumour-suppressor gene.