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.

Genetic variability and potential effects on clinical trial outcomes: perspectives in Parkinson's disease.

BACKGROUND: Classical randomisation of clinical trial patients creates a source of genetic variance that may be contributing to the high failure rate seen in neurodegenerative disease trials. Our objective was to quantify genetic difference between randomised trial arms and determine how imbalance can affect trial outcomes. METHODS: 5851 patients with Parkinson's disease of European ancestry data and two simulated virtual cohorts based on public data were used. Data were resampled at different sizes for 1000 iterations and randomly assigned to the two arms of a simulated trial. False-negative and false-positive rates were estimated using simulated clinical trials, and per cent difference in genetic risk score (GRS) and allele frequency was calculated to quantify variance between arms. RESULTS: 5851 patients with Parkinson's disease (mean (SD) age, 61.02 (12.61) years; 2095 women (35.81%)) as well as simulated patients from virtually created cohorts were used in the study. Approximately 90% of the iterations had at least one statistically significant difference in individual risk SNPs between each trial arm. Approximately 5%-6% of iterations had a statistically significant difference between trial arms in mean GRS. For significant iterations, the average per cent difference for mean GRS between trial arms was 130.87%, 95% CI 120.89 to 140.85 (n=200). Glucocerebrocidase (GBA) gene-only simulations see an average 18.86%, 95% CI 18.01 to 19.71 difference in GRS scores between trial arms (n=50). When adding a drug effect of -0.5 points in MDS-UPDRS per year at n=50, 33.9% of trials resulted in false negatives. CONCLUSIONS: Our data support the hypothesis that within genetically unmatched clinical trials, genetic heterogeneity could confound true therapeutic effects as expected. Clinical trials should undergo pretrial genetic adjustment or, at the minimum, post-trial adjustment and analysis for failed trials.

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.

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.

PINES: phenotype-informed tissue weighting improves prediction of pathogenic noncoding variants.

Functional characterization of the noncoding genome is essential for biological understanding of gene regulation and disease. Here, we introduce the computational framework PINES (Phenotype-Informed Noncoding Element Scoring), which predicts the functional impact of noncoding variants by integrating epigenetic annotations in a phenotype-dependent manner. PINES enables analyses to be customized towards genomic annotations from cell types of the highest relevance given the phenotype of interest. We illustrate that PINES identifies functional noncoding variation more accurately than methods that do not use phenotype-weighted knowledge, while at the same time being flexible and easy to use via a dedicated web portal.

A whole-genome sequence study identifies genetic risk factors for neuromyelitis optica.

Neuromyelitis optica (NMO) is a rare autoimmune disease that affects the optic nerve and spinal cord. Most NMO patients ( > 70%) are seropositive for circulating autoantibodies against aquaporin 4 (NMO-IgG+). Here, we meta-analyze whole-genome sequences from 86 NMO cases and 460 controls with genome-wide SNP array from 129 NMO cases and 784 controls to test for association with SNPs and copy number variation (total N = 215 NMO cases, 1244 controls). We identify two independent signals in the major histocompatibility complex (MHC) region associated with NMO-IgG+, one of which may be explained by structural variation in the complement component 4 genes. Mendelian Randomization analysis reveals a significant causal effect of known systemic lupus erythematosus (SLE), but not multiple sclerosis (MS), risk variants in NMO-IgG+. Our results suggest that genetic variants in the MHC region contribute to the etiology of NMO-IgG+ and that NMO-IgG+ is genetically more similar to SLE than MS.

Erratum: Sequence data and association statistics from 12,940 type 2 diabetes cases and controls.

This corrects the article DOI: 10.1038/sdata.2017.179.

Sequence data and association statistics from 12,940 type 2 diabetes cases and controls.

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

Radiographic endophenotyping in hip osteoarthritis improves the precision of genetic association analysis.

OBJECTIVE: Osteoarthritis (OA) has a strong genetic component but the success of previous genome-wide association studies (GWAS) has been restricted due to insufficient sample sizes and phenotype heterogeneity. Our aim was to examine the effect of clinically relevant endophenotyping according to site of maximal joint space narrowing (maxJSN) and bone remodelling response on GWAS signal detection in hip OA. METHODS: A stratified GWAS meta-analysis was conducted in 2118 radiographically defined hip OA cases and 6500 population-based controls. Signals were followed up by analysing differential expression of proximal genes for bone remodelling endophenotypes in 33 pairs of macroscopically intact and OA-affected cartilage. RESULTS: We report suggestive evidence (p<5×10-6) of association at 6 variants with OA endophenotypes that would have been missed by using presence of hip OA as the disease end point. For example, in the analysis of hip OA cases with superior maxJSN versus cases with non-superior maxJSN we detected association with a variant in the LRCH1 gene (rs754106, p=1.49×10-7, OR (95% CIs) 0.70 (0.61 to 0.80)). In the comparison of hypertrophic with non-hypertrophic OA the most significant variant was located between STT3B and GADL1 (rs6766414, p=3.13×10-6, OR (95% CIs) 1.45 (1.24 to 1.69)). Both of these associations were fully attenuated in non-stratified analyses of all hip OA cases versus population controls (p>0.05). STT3B was significantly upregulated in OA-affected versus intact cartilage, particularly in the analysis of hypertrophic and normotrophic compared with atrophic bone remodelling pattern (p=4.2×10-4). CONCLUSIONS: Our findings demonstrate that stratification of OA cases into more homogeneous endophenotypes can identify genes of potential functional importance otherwise obscured by disease heterogeneity.

The genetic architecture of type 2 diabetes.

The genetic architecture of common traits, including the number, frequency, and effect sizes of inherited variants that contribute to individual risk, has been long debated. Genome-wide association studies have identified scores of common variants associated with type 2 diabetes, but in aggregate, these explain only a fraction of the heritability of this disease. Here, to test the hypothesis that lower-frequency variants explain much of the remainder, the GoT2D and T2D-GENES consortia performed whole-genome sequencing in 2,657 European individuals with and without diabetes, and exome sequencing in 12,940 individuals from five ancestry groups. To increase statistical power, we expanded the sample size via genotyping and imputation in a further 111,548 subjects. Variants associated with type 2 diabetes after sequencing were overwhelmingly common and most fell within regions previously identified by genome-wide association studies. Comprehensive enumeration of sequence variation is necessary to identify functional alleles that provide important clues to disease pathophysiology, but large-scale sequencing does not support the idea that lower-frequency variants have a major role in predisposition to type 2 diabetes.

Diagnosis of Parkinson's disease on the basis of clinical and genetic classification: a population-based modelling study.

BACKGROUND: Accurate diagnosis and early detection of complex diseases, such as Parkinson's disease, has the potential to be of great benefit for researchers and clinical practice. We aimed to create a non-invasive, accurate classification model for the diagnosis of Parkinson's disease, which could serve as a basis for future disease prediction studies in longitudinal cohorts. METHODS: We developed a model for disease classification using data from the Parkinson's Progression Marker Initiative (PPMI) study for 367 patients with Parkinson's disease and phenotypically typical imaging data and 165 controls without neurological disease. Olfactory function, genetic risk, family history of Parkinson's disease, age, and gender were algorithmically selected by stepwise logistic regression as significant contributors to our classifying model. We then tested the model with data from 825 patients with Parkinson's disease and 261 controls from five independent cohorts with varying recruitment strategies and designs: the Parkinson's Disease Biomarkers Program (PDBP), the Parkinson's Associated Risk Study (PARS), 23andMe, the Longitudinal and Biomarker Study in PD (LABS-PD), and the Morris K Udall Parkinson's Disease Research Center of Excellence cohort (Penn-Udall). Additionally, we used our model to investigate patients who had imaging scans without evidence of dopaminergic deficit (SWEDD). FINDINGS: In the population from PPMI, our initial model correctly distinguished patients with Parkinson's disease from controls at an area under the curve (AUC) of 0·923 (95% CI 0·900-0·946) with high sensitivity (0·834, 95% CI 0·711-0·883) and specificity (0·903, 95% CI 0·824-0·946) at its optimum AUC threshold (0·655). All Hosmer-Lemeshow simulations suggested that when parsed into random subgroups, the subgroup data matched that of the overall cohort. External validation showed good classification of Parkinson's disease, with AUCs of 0·894 (95% CI 0·867-0·921) in the PDBP cohort, 0·998 (0·992-1·000) in PARS, 0·955 (no 95% CI available) in 23andMe, 0·929 (0·896-0·962) in LABS-PD, and 0·939 (0·891-0·986) in the Penn-Udall cohort. Four of 17 SWEDD participants who our model classified as having Parkinson's disease converted to Parkinson's disease within 1 year, whereas only one of 38 SWEDD participants who were not classified as having Parkinson's disease underwent conversion (test of proportions, p=0·003). INTERPRETATION: Our model provides a potential new approach to distinguish participants with Parkinson's disease from controls. If the model can also identify individuals with prodromal or preclinical Parkinson's disease in prospective cohorts, it could facilitate identification of biomarkers and interventions. FUNDING: National Institute on Aging, National Institute of Neurological Disorders and Stroke, and the Michael J Fox Foundation.

Investigation of association between hip osteoarthritis susceptibility loci and radiographic proximal femur shape.

OBJECTIVE: To test whether previously reported hip morphology or osteoarthritis (OA) susceptibility loci are associated with proximal femur shape as represented by statistical shape model (SSM) modes and as univariate or multivariate quantitative traits. METHODS: We used pelvic radiographs and genotype data from 929 subjects with unilateral hip OA who had been recruited previously for the Arthritis Research UK Osteoarthritis Genetics Consortium genome-wide association study. We built 3 SSMs capturing the shape variation of the OA-unaffected proximal femur in the entire mixed-sex cohort and for male/female-stratified cohorts. We selected 41 candidate single-nucleotide polymorphisms (SNPs) previously reported as being associated with hip morphology (for replication analysis) or OA (for discovery analysis) and for which genotype data were available. We performed 2 types of analysis for genotype-phenotype associations between these SNPs and the modes of the SSMs: 1) a univariate analysis using individual SSM modes and 2) a multivariate analysis using combinations of SSM modes. RESULTS: The univariate analysis identified association between rs4836732 (within the ASTN2 gene) and mode 5 of the female SSM (P = 0.0016) and between rs6976 (within the GLT8D1 gene) and mode 7 of the mixed-sex SSM (P = 0.0003). The multivariate analysis identified association between rs5009270 (near the IFRD1 gene) and a combination of modes 3, 4, and 9 of the mixed-sex SSM (P = 0.0004). Evidence of associations remained significant following adjustment for multiple testing. All 3 SNPs had previously been associated with hip OA. CONCLUSION: These de novo findings suggest that rs4836732, rs6976, and rs5009270 may contribute to hip OA susceptibility by altering proximal femur shape.

Exome sequencing in amyotrophic lateral sclerosis identifies risk genes and pathways.

Amyotrophic lateral sclerosis (ALS) is a devastating neurological disease with no effective treatment. We report the results of a moderate-scale sequencing study aimed at increasing the number of genes known to contribute to predisposition for ALS. We performed whole-exome sequencing of 2869 ALS patients and 6405 controls. Several known ALS genes were found to be associated, and TBK1 (the gene encoding TANK-binding kinase 1) was identified as an ALS gene. TBK1 is known to bind to and phosphorylate a number of proteins involved in innate immunity and autophagy, including optineurin (OPTN) and p62 (SQSTM1/sequestosome), both of which have also been implicated in ALS. These observations reveal a key role of the autophagic pathway in ALS and suggest specific targets for therapeutic intervention.

Whole Genome Sequencing Reveals a Chromosome 9p Deletion Causing DOCK8 Deficiency in an Adult Diagnosed with Hyper IgE Syndrome Who Developed Progressive Multifocal Leukoencephalopathy.

PURPOSE: A 30 year-old man with a history of recurrent skin infections as well as elevated serum IgE and eosinophils developed neurological symptoms and had T2-hyperintense lesions observed in cerebral MRI. The immune symptoms were attributed to Hyper IgE syndrome (HIES) and the neurological symptoms with presence of JC virus in cerebrospinal fluid were diagnosed as Progressive Multifocal Leukoencephalopathy (PML). The patient was negative for STAT3 mutations. To determine if other mutations explain HIES and/or PML in this subject, his DNA was analyzed by whole genome sequencing. METHODS: Whole genome sequencing was completed to 30X coverage, and whole genome SNP typing was used to complement these data. The methods revealed single nucleotide variants, structural variants, and copy number variants across the genome. Genome-wide data were analyzed for homozygous or compound heterozygous null mutations for all protein coding genes. Mutations were confirmed by PCR and/or Sanger sequencing. RESULTS: Whole genome analysis revealed deletions near the telomere of both copies of chromosome 9p. Several genes, including DOCK8, were impacted by the deletions but it was unclear whether each chromosome had identical or distinct deletions. PCR across the impacted region combined with Sanger sequencing of selected fragments confirmed a homozygous deletion from position 10,211 to 586,751. CONCLUSION: While several genes are impacted by the deletion, DOCK8 deficiency is the most probable cause of HIES in this patient. DOCK8 deficiency may have also predisposed the patient to develop PML.

A meta-analysis of genome-wide association studies identifies novel variants associated with osteoarthritis of the hip.

OBJECTIVES: Osteoarthritis (OA) is the most common form of arthritis with a clear genetic component. To identify novel loci associated with hip OA we performed a meta-analysis of genome-wide association studies (GWAS) on European subjects. METHODS: We performed a two-stage meta-analysis on more than 78,000 participants. In stage 1, we synthesised data from eight GWAS whereas data from 10 centres were used for 'in silico' or 'de novo' replication. Besides the main analysis, a stratified by sex analysis was performed to detect possible sex-specific signals. Meta-analysis was performed using inverse-variance fixed effects models. A random effects approach was also used. RESULTS: We accumulated 11,277 cases of radiographic and symptomatic hip OA. We prioritised eight single nucleotide polymorphism (SNPs) for follow-up in the discovery stage (4349 OA cases); five from the combined analysis, two male specific and one female specific. One locus, at 20q13, represented by rs6094710 (minor allele frequency (MAF) 4%) near the NCOA3 (nuclear receptor coactivator 3) gene, reached genome-wide significance level with p=7.9×10(-9) and OR=1.28 (95% CI 1.18 to 1.39) in the combined analysis of discovery (p=5.6×10(-8)) and follow-up studies (p=7.3×10(-4)). We showed that this gene is expressed in articular cartilage and its expression was significantly reduced in OA-affected cartilage. Moreover, two loci remained suggestive associated; rs5009270 at 7q31 (MAF 30%, p=9.9×10(-7), OR=1.10) and rs3757837 at 7p13 (MAF 6%, p=2.2×10(-6), OR=1.27 in male specific analysis). CONCLUSIONS: Novel genetic loci for hip OA were found in this meta-analysis of GWAS.

Evaluation of the genetic overlap between osteoarthritis with body mass index and height using genome-wide association scan data.

OBJECTIVES: Obesity as measured by body mass index (BMI) is one of the major risk factors for osteoarthritis. In addition, genetic overlap has been reported between osteoarthritis and normal adult height variation. We investigated whether this relationship is due to a shared genetic aetiology on a genome-wide scale. METHODS: We compared genetic association summary statistics (effect size, p value) for BMI and height from the GIANT consortium genome-wide association study (GWAS) with genetic association summary statistics from the arcOGEN consortium osteoarthritis GWAS. Significance was evaluated by permutation. Replication of osteoarthritis association of the highlighted signals was investigated in an independent dataset. Phenotypic information of height and BMI was accounted for in a separate analysis using osteoarthritis-free controls. RESULTS: We found significant overlap between osteoarthritis and height (p=3.3×10(-5) for signals with p≤0.05) when the GIANT and arcOGEN GWAS were compared. For signals with p≤0.001 we found 17 shared signals between osteoarthritis and height and four between osteoarthritis and BMI. However, only one of the height or BMI signals that had shown evidence of association with osteoarthritis in the arcOGEN GWAS was also associated with osteoarthritis in the independent dataset: rs12149832, within the FTO gene (combined p=2.3×10(-5)). As expected, this signal was attenuated when we adjusted for BMI. CONCLUSIONS: We found a significant excess of shared signals between both osteoarthritis and height and osteoarthritis and BMI, suggestive of a common genetic aetiology. However, only one signal showed association with osteoarthritis when followed up in a new dataset.

Genome-wide association analysis of imputed rare variants: application to seven common complex diseases.

Genome-wide association studies have been successful in identifying loci contributing effects to a range of complex human traits. The majority of reproducible associations within these loci are with common variants, each of modest effect, which together explain only a small proportion of heritability. It has been suggested that much of the unexplained genetic component of complex traits can thus be attributed to rare variation. However, genome-wide association study genotyping chips have been designed primarily to capture common variation, and thus are underpowered to detect the effects of rare variants. Nevertheless, we demonstrate here, by simulation, that imputation from an existing scaffold of genome-wide genotype data up to high-density reference panels has the potential to identify rare variant associations with complex traits, without the need for costly re-sequencing experiments. By application of this approach to genome-wide association studies of seven common complex diseases, imputed up to publicly available reference panels, we identify genome-wide significant evidence of rare variant association in PRDM10 with coronary artery disease and multiple genes in the major histocompatibility complex (MHC) with type 1 diabetes. The results of our analyses highlight that genome-wide association studies have the potential to offer an exciting opportunity for gene discovery through association with rare variants, conceivably leading to substantial advancements in our understanding of the genetic architecture underlying complex human traits.

Identification of new susceptibility loci for osteoarthritis (arcOGEN): a genome-wide association study.

BACKGROUND: Osteoarthritis is the most common form of arthritis worldwide and is a major cause of pain and disability in elderly people. The health economic burden of osteoarthritis is increasing commensurate with obesity prevalence and longevity. Osteoarthritis has a strong genetic component but the success of previous genetic studies has been restricted due to insufficient sample sizes and phenotype heterogeneity. METHODS: We undertook a large genome-wide association study (GWAS) in 7410 unrelated and retrospectively and prospectively selected patients with severe osteoarthritis in the arcOGEN study, 80% of whom had undergone total joint replacement, and 11,009 unrelated controls from the UK. We replicated the most promising signals in an independent set of up to 7473 cases and 42,938 controls, from studies in Iceland, Estonia, the Netherlands, and the UK. All patients and controls were of European descent. FINDINGS: We identified five genome-wide significant loci (binomial test p≤5·0×10(-8)) for association with osteoarthritis and three loci just below this threshold. The strongest association was on chromosome 3 with rs6976 (odds ratio 1·12 [95% CI 1·08-1·16]; p=7·24×10(-11)), which is in perfect linkage disequilibrium with rs11177. This SNP encodes a missense polymorphism within the nucleostemin-encoding gene GNL3. Levels of nucleostemin were raised in chondrocytes from patients with osteoarthritis in functional studies. Other significant loci were on chromosome 9 close to ASTN2, chromosome 6 between FILIP1 and SENP6, chromosome 12 close to KLHDC5 and PTHLH, and in another region of chromosome 12 close to CHST11. One of the signals close to genome-wide significance was within the FTO gene, which is involved in regulation of bodyweight-a strong risk factor for osteoarthritis. All risk variants were common in frequency and exerted small effects. INTERPRETATION: Our findings provide insight into the genetics of arthritis and identify new pathways that might be amenable to future therapeutic intervention. FUNDING: arcOGEN was funded by a special purpose grant from Arthritis Research UK.