Multivariate profile and acute-phase correlates of cognitive deficits in a COVID-19 hospitalised cohort.

Background: Preliminary evidence has highlighted a possible association between severe COVID-19 and persistent cognitive deficits. Further research is required to confirm this association, determine whether cognitive deficits relate to clinical features from the acute phase or to mental health status at the point of assessment, and quantify rate of recovery. Methods: 46 individuals who received critical care for COVID-19 at Addenbrooke's hospital between 10th March 2020 and 31st July 2020 (16 mechanically ventilated) underwent detailed computerised cognitive assessment alongside scales measuring anxiety, depression and post-traumatic stress disorder under supervised conditions at a mean follow up of 6.0 (± 2.1) months following acute illness. Patient and matched control (N = 460) performances were transformed into standard deviation from expected scores, accounting for age and demographic factors using N = 66,008 normative datasets. Global accuracy and response time composites were calculated (G_SScore & G_RT). Linear modelling predicted composite score deficits from acute severity, mental-health status at assessment, and time from hospital admission. The pattern of deficits across tasks was qualitatively compared with normal age-related decline, and early-stage dementia. Findings: COVID-19 survivors were less accurate (G_SScore=-0.53SDs) and slower (G_RT=+0.89SDs) in their responses than expected compared to their matched controls. Acute illness, but not chronic mental health, significantly predicted cognitive deviation from expected scores (G_SScore (p=​​0.0037) and G_RT (p = 0.0366)). The most prominent task associations with COVID-19 were for higher cognition and processing speed, which was qualitatively distinct from the profiles of normal ageing and dementia and similar in magnitude to the effects of ageing between 50 and 70 years of age. A trend towards reduced deficits with time from illness (r∼=0.15) did not reach statistical significance. Interpretation: Cognitive deficits after severe COVID-19 relate most strongly to acute illness severity, persist long into the chronic phase, and recover slowly if at all, with a characteristic profile highlighting higher cognitive functions and processing speed. Funding: This work was funded by the National Institute for Health Research (NIHR) Cambridge Biomedical Research Centre (BRC), NIHR Cambridge Clinical Research Facility (BRC-1215-20014), the Addenbrooke's Charities Trust and NIHR COVID-19 BioResource RG9402. AH is funded by the UK Dementia Research Institute Care Research and Technology Centre and Imperial College London Biomedical Research Centre. ETB and DKM are supported by NIHR Senior Investigator awards. JBR is supported by the Wellcome Trust (220258) and Medical Research Council (SUAG/051 G101400). VFJN is funded by an Academy of Medical Sciences/ The Health Foundation Clinician Scientist Fellowship. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care.

Genetic associations at regulatory phenotypes improve fine-mapping of causal variants for 12 immune-mediated diseases.

The resolution of causal genetic variants informs understanding of disease biology. We used regulatory quantitative trait loci (QTLs) from the BLUEPRINT, GTEx and eQTLGen projects to fine-map putative causal variants for 12 immune-mediated diseases. We identify 340 unique loci that colocalize with high posterior probability (≥98%) with regulatory QTLs and apply Bayesian frameworks to fine-map associations at each locus. We show that fine-mapping credible sets derived from regulatory QTLs are smaller compared to disease summary statistics. Further, they are enriched for more functionally interpretable candidate causal variants and for putatively causal insertion/deletion (INDEL) polymorphisms. Finally, we use massively parallel reporter assays to evaluate candidate causal variants at the ITGA4 locus associated with inflammatory bowel disease. Overall, our findings suggest that fine-mapping applied to disease-colocalizing regulatory QTLs can enhance the discovery of putative causal disease variants and enhance insights into the underlying causal genes and molecular mechanisms.

Effect of age at puberty on risk of multiple sclerosis: A mendelian randomization study.

OBJECTIVE: To investigate the potential for a causal effect of age at puberty on multiple sclerosis (MS) susceptibility using a mendelian randomization (MR) approach. METHODS: We used 372 genetic variants strongly associated with age at menarche in a genome-wide association study (GWAS) involving 329,245 women. The genetic architecture of pubertal timing across both sexes is highly correlated (genetic correlation [r g] = 0.75, p = 1.2 × 10-79), allowing these variants to provide reliable insight into pubertal timing in males as well. The effect of pubertal timing on risk of MS was measured with summary statistics from a GWAS of 14,802 cases with MS and 26,703 controls from the International Multiple Sclerosis Genetics Consortium. Multivariable MR controlling for effects of body mass index (BMI) using genetic data from additional consortia investigated whether pubertal effects on MS were dependent on weight status. RESULTS: A 1-year increase in genetically predicted age at puberty decreased odds of MS by 8% (odds ratio [OR] 0.92, 95% confidence interval [CI] 0.86-0.99, p = 0.03). However, multivariable MR analysis showed that after accounting for effects on adult BMI, the association of age at puberty with MS susceptibility attenuated (OR 0.96, 95% CI 0.88-1.04, p = 0.36). Similar results were obtained when childhood BMI was incorporated. Sensitivity analyses provided no evidence of major bias from genetic pleiotropy. CONCLUSIONS: We found support for an association between higher age at puberty and decreased risk of MS with a magnitude comparable to that reported in observational studies. This effect appears to be largely mediated by the strong association between age at puberty and obesity. A large causal effect of pubertal timing independent of BMI is unlikely.

Population- and individual-specific regulatory variation in Sardinia.

Genetic studies of complex traits have mainly identified associations with noncoding variants. To further determine the contribution of regulatory variation, we combined whole-genome and transcriptome data for 624 individuals from Sardinia to identify common and rare variants that influence gene expression and splicing. We identified 21,183 expression quantitative trait loci (eQTLs) and 6,768 splicing quantitative trait loci (sQTLs), including 619 new QTLs. We identified high-frequency QTLs and found evidence of selection near genes involved in malarial resistance and increased multiple sclerosis risk, reflecting the epidemiological history of Sardinia. Using family relationships, we identified 809 segregating expression outliers (median z score of 2.97), averaging 13.3 genes per individual. Outlier genes were enriched for proximal rare variants, providing a new approach to study large-effect regulatory variants and their relevance to traits. Our results provide insight into the effects of regulatory variants and their relationship to population history and individual genetic risk.

Immunoseq: the identification of functionally relevant variants through targeted capture and sequencing of active regulatory regions in human immune cells.

BACKGROUND: The observation that the genetic variants identified in genome-wide association studies (GWAS) frequently lie in non-coding regions of the genome that contain cis-regulatory elements suggests that altered gene expression underlies the development of many complex traits. In order to efficiently make a comprehensive assessment of the impact of non-coding genetic variation in immune related diseases we emulated the whole-exome sequencing paradigm and developed a custom capture panel for the known DNase I hypersensitive site (DHS) in immune cells - "Immunoseq". RESULTS: We performed Immunoseq in 30 healthy individuals where we had existing transcriptome data from T cells. We identified a large number of novel non-coding variants in these samples. Relying on allele specific expression measurements, we also showed that our selected capture regions are enriched for functional variants that have an impact on differential allelic gene expression. The results from a replication set with 180 samples confirmed our observations. CONCLUSIONS: We show that Immunoseq is a powerful approach to detect novel rare variants in regulatory regions. We also demonstrate that these novel variants have a potential functional role in immune cells.

Polymorphism in a lincRNA Associates with a Doubled Risk of Pneumococcal Bacteremia in Kenyan Children.

Bacteremia (bacterial bloodstream infection) is a major cause of illness and death in sub-Saharan Africa but little is known about the role of human genetics in susceptibility. We conducted a genome-wide association study of bacteremia susceptibility in more than 5,000 Kenyan children as part of the Wellcome Trust Case Control Consortium 2 (WTCCC2). Both the blood-culture-proven bacteremia case subjects and healthy infants as controls were recruited from Kilifi, on the east coast of Kenya. Streptococcus pneumoniae is the most common cause of bacteremia in Kilifi and was thus the focus of this study. We identified an association between polymorphisms in a long intergenic non-coding RNA (lincRNA) gene (AC011288.2) and pneumococcal bacteremia and replicated the results in the same population (p combined = 1.69 × 10(-9); OR = 2.47, 95% CI = 1.84-3.31). The susceptibility allele is African specific, derived rather than ancestral, and occurs at low frequency (2.7% in control subjects and 6.4% in case subjects). Our further studies showed AC011288.2 expression only in neutrophils, a cell type that is known to play a major role in pneumococcal clearance. Identification of this novel association will further focus research on the role of lincRNAs in human infectious disease.

Experience from two decades of the Cambridge Rapid Access Neurology Clinic.

We report on the evolution of the rapid access neurology clinic (established in 1995) at Addenbrooke's Hospital, Cambridge. Annualised attendance data demonstrate an ever increasing demand, with primary headache disorders now accounting for more than 40% of referrals. Secondary causes of headache (including intracranial tumours, idiopathic intracranial hypertension, carotid or vertebral artery dissection and subdural haematomas) remain infrequent. In all such cases, there were additional diagnostic clues. The number of patients referred with problems related to chronic neurological diseases has fallen considerably, reflecting the roles of specialist nurses and clinics. Imaging investigation of choice shifted from computerised tomography scan (45 to 16%) towards magnetic resonance imaging (17 to 47%). Management is increasingly on an outpatient basis, often without the need for a follow-up appointment. The experience presented here should inform further development of rapid access neurology clinics across the UK and suggests the need for acute headache services, in line with those for transient ischaemic attack and first seizure.

Dissection of a Complex Disease Susceptibility Region Using a Bayesian Stochastic Search Approach to Fine Mapping.

Identification of candidate causal variants in regions associated with risk of common diseases is complicated by linkage disequilibrium (LD) and multiple association signals. Nonetheless, accurate maps of these variants are needed, both to fully exploit detailed cell specific chromatin annotation data to highlight disease causal mechanisms and cells, and for design of the functional studies that will ultimately be required to confirm causal mechanisms. We adapted a Bayesian evolutionary stochastic search algorithm to the fine mapping problem, and demonstrated its improved performance over conventional stepwise and regularised regression through simulation studies. We then applied it to fine map the established multiple sclerosis (MS) and type 1 diabetes (T1D) associations in the IL-2RA (CD25) gene region. For T1D, both stepwise and stochastic search approaches identified four T1D association signals, with the major effect tagged by the single nucleotide polymorphism, rs12722496. In contrast, for MS, the stochastic search found two distinct competing models: a single candidate causal variant, tagged by rs2104286 and reported previously using stepwise analysis; and a more complex model with two association signals, one of which was tagged by the major T1D associated rs12722496 and the other by rs56382813. There is low to moderate LD between rs2104286 and both rs12722496 and rs56382813 (r2 ≃ 0:3) and our two SNP model could not be recovered through a forward stepwise search after conditioning on rs2104286. Both signals in the two variant model for MS affect CD25 expression on distinct subpopulations of CD4+ T cells, which are key cells in the autoimmune process. The results support a shared causal variant for T1D and MS. Our study illustrates the benefit of using a purposely designed model search strategy for fine mapping and the advantage of combining disease and protein expression data.

Statistical colocalization of genetic risk variants for related autoimmune diseases in the context of common controls.

Determining whether potential causal variants for related diseases are shared can identify overlapping etiologies of multifactorial disorders. Colocalization methods disentangle shared and distinct causal variants. However, existing approaches require independent data sets. Here we extend two colocalization methods to allow for the shared-control design commonly used in comparison of genome-wide association study results across diseases. Our analysis of four autoimmune diseases--type 1 diabetes (T1D), rheumatoid arthritis, celiac disease and multiple sclerosis--identified 90 regions that were associated with at least one disease, 33 (37%) of which were associated with 2 or more disorders. Nevertheless, for 14 of these 33 shared regions, there was evidence that the causal variants differed. We identified new disease associations in 11 regions previously associated with one or more of the other 3 disorders. Four of eight T1D-specific regions contained known type 2 diabetes (T2D) candidate genes (COBL, GLIS3, RNLS and BCAR1), suggesting a shared cellular etiology.

Genetic variants are major determinants of CSF antibody levels in multiple sclerosis.

Immunological hallmarks of multiple sclerosis include the production of antibodies in the central nervous system, expressed as presence of oligoclonal bands and/or an increased immunoglobulin G index-the level of immunoglobulin G in the cerebrospinal fluid compared to serum. However, the underlying differences between oligoclonal band-positive and -negative patients with multiple sclerosis and reasons for variability in immunoglobulin G index are not known. To identify genetic factors influencing the variation in the antibody levels in the cerebrospinal fluid in multiple sclerosis, we have performed a genome-wide association screen in patients collected from nine countries for two traits, presence or absence of oligoclonal bands (n = 3026) and immunoglobulin G index levels (n = 938), followed by a replication in 3891 additional patients. We replicate previously suggested association signals for oligoclonal band status in the major histocompatibility complex region for the rs9271640*A-rs6457617*G haplotype, correlated with HLA-DRB1*1501, and rs34083746*G, correlated with HLA-DQA1*0301 (P comparing two haplotypes = 8.88 × 10(-16)). Furthermore, we identify a novel association signal of rs9807334, near the ELAC1/SMAD4 genes, for oligoclonal band status (P = 8.45 × 10(-7)). The previously reported association of the immunoglobulin heavy chain locus with immunoglobulin G index reaches strong evidence for association in this data set (P = 3.79 × 10(-37)). We identify two novel associations in the major histocompatibility complex region with immunoglobulin G index: the rs9271640*A-rs6457617*G haplotype (P = 1.59 × 10(-22)), shared with oligoclonal band status, and an additional independent effect of rs6457617*G (P = 3.68 × 10(-6)). Variants identified in this study account for up to 2-fold differences in the odds of being oligoclonal band positive and 7.75% of the variation in immunoglobulin G index. Both traits are associated with clinical features of disease such as female gender, age at onset and severity. This is the largest study population so far investigated for the genetic influence on antibody levels in the cerebrospinal fluid in multiple sclerosis, including 6950 patients. We confirm that genetic factors underlie these antibody levels and identify both the major histocompatibility complex and immunoglobulin heavy chain region as major determinants.

The correlation between reading and mathematics ability at age twelve has a substantial genetic component.

Dissecting how genetic and environmental influences impact on learning is helpful for maximizing numeracy and literacy. Here we show, using twin and genome-wide analysis, that there is a substantial genetic component to children's ability in reading and mathematics, and estimate that around one half of the observed correlation in these traits is due to shared genetic effects (so-called Generalist Genes). Thus, our results highlight the potential role of the learning environment in contributing to differences in a child's cognitive abilities at age twelve.

No evidence for shared genetic basis of common variants in multiple sclerosis and amyotrophic lateral sclerosis.

Genome-wide association studies have been successful in identifying common variants that influence the susceptibility to complex diseases. From these studies, it has emerged that there is substantial overlap in susceptibility loci between diseases. In line with those findings, we hypothesized that shared genetic pathways may exist between multiple sclerosis (MS) and amyotrophic lateral sclerosis (ALS). While both diseases may have inflammatory and neurodegenerative features, epidemiological studies have indicated an increased co-occurrence within individuals and families. To this purpose, we combined genome-wide data from 4088 MS patients, 3762 ALS patients and 12 030 healthy control individuals in whom 5 440 446 single-nucleotide polymorphisms (SNPs) were successfully genotyped or imputed. We tested these SNPs for the excess association shared between MS and ALS and also explored whether polygenic models of SNPs below genome-wide significance could explain some of the observed trait variance between diseases. Genome-wide association meta-analysis of SNPs as well as polygenic analyses fails to provide evidence in favor of an overlap in genetic susceptibility between MS and ALS. Hence, our findings do not support a shared genetic background of common risk variants in MS and ALS.

Glucocerebrosidase mutations influence the natural history of Parkinson's disease in a community-based incident cohort.

Carriers of mutations in the glucocerebrosidase gene (GBA) are at increased risk of developing Parkinson's disease. The frequency of GBA mutations in unselected Parkinson's disease populations has not been established. Furthermore, no previous studies have investigated the influence of GBA mutations on the natural history of Parkinson's disease using prospective follow-up. We studied DNA from 262 cases who had been recruited at diagnosis into one of two independent community-based incidence studies of Parkinson's disease. In 121 cases, longitudinal data regarding progression of motor disability and cognitive function were derived from follow-up assessments conducted every 18 months for a median of 71 months. Sequencing of the GBA was performed after two-stage polymerase chain reaction amplification. The carrier frequency of genetic variants in GBA was determined. Baseline demographic and clinical variables were compared between cases who were either GBA mutation carriers, polymorphism carriers or wild-type homozygotes. Cox regression analysis was used to model progression to major motor (Hoehn and Yahr stage 3), and cognitive (dementia) end-points in cases followed longitudinally. We show that in a representative, unselected UK Parkinson's disease population, GBA mutations are present at a frequency of 3.5%. This is higher than the prevalence of other genetic mutations currently associated with Parkinson's disease and indicates that GBA mutations make an important contribution to Parkinson's disease encountered in the community setting. Baseline clinical characteristics did not differ significantly between cases with and without GBA sequence variants. However, the hazard ratio for progression both to dementia (5.7, P = 0.003) and Hoehn and Yahr stage 3 (4.2, P = 0.003) were significantly greater in GBA mutation carriers. We also show that carriers of polymorphisms in GBA which are not generally considered to increase Parkinson's disease risk are at significantly increased risk of progression to Hoehn and Yahr stage 3 (3.2, P = 0.004). Our results indicate that genetic variation in GBA has an important impact on the natural history of Parkinson's disease. To our knowledge, this is the first time a genetic locus has been shown to influence motor progression in Parkinson's disease. If confirmed in further studies, this may indicate that GBA mutation status could be used as a prognostic marker in Parkinson's disease. Elucidation of the molecular mechanisms that underlie this effect will further our understanding of the pathogenesis of the disease and may in turn suggest novel therapeutic strategies.

Common variants in the HLA-DRB1-HLA-DQA1 HLA class II region are associated with susceptibility to visceral leishmaniasis.

To identify susceptibility loci for visceral leishmaniasis, we undertook genome-wide association studies in two populations: 989 cases and 1,089 controls from India and 357 cases in 308 Brazilian families (1,970 individuals). The HLA-DRB1-HLA-DQA1 locus was the only region to show strong evidence of association in both populations. Replication at this region was undertaken in a second Indian population comprising 941 cases and 990 controls, and combined analysis across the three cohorts for rs9271858 at this locus showed P(combined) = 2.76 × 10(-17) and odds ratio (OR) = 1.41, 95% confidence interval (CI) = 1.30-1.52. A conditional analysis provided evidence for multiple associations within the HLA-DRB1-HLA-DQA1 region, and a model in which risk differed between three groups of haplotypes better explained the signal and was significant in the Indian discovery and replication cohorts. In conclusion, the HLA-DRB1-HLA-DQA1 HLA class II region contributes to visceral leishmaniasis susceptibility in India and Brazil, suggesting shared genetic risk factors for visceral leishmaniasis that cross the epidemiological divides of geography and parasite species.

Systematic review and UK-based study of PARK2 (parkin), PINK1, PARK7 (DJ-1) and LRRK2 in early-onset Parkinson's disease.

Approximately 3.6% of patients with Parkinson's disease develop symptoms before age 45. Early-onset Parkinson's disease (EOPD) patients have a higher familial recurrence risk than late-onset patients, and 3 main recessive EOPD genes have been described. We aimed to establish the prevalence of mutations in these genes in a UK cohort and in previous studies. We screened 136 EOPD probands from a high-ascertainment regional and community-based prevalence study for pathogenic mutations in PARK2 (parkin), PINK1, PARK7 (DJ-1), and exon 41 of LRRK2. We also carried out a systematic review, calculating the proportion of cases with pathogenic mutations in previously reported studies. We identified 5 patients with pathogenic PARK2, 1 patient with PINK1, and 1 with LRRK2 mutations. The rate of mutations overall was 5.1%. Mutations were more common in patients with age at onset (AAO) < 40 (9.5%), an affected first-degree relative (6.9%), an affected sibling (28.6%), or parental consanguinity (50%). In our study EOPD mutation carriers were more likely to present with rigidity and dystonia, and 6 of 7 mutation carriers had lower limb symptoms at onset. Our systematic review included information from >5800 unique cases. Overall, the weighted mean proportion of cases with PARK2 (parkin), PINK1, and PARK7 (DJ-1) mutations was 8.6%, 3.7%, and 0.4%, respectively. PINK1 mutations were more common in Asian subjects. The overall frequency of mutations in known EOPD genes was lower than previously estimated. Our study shows an increased likelihood of mutations in patients with lower AAO, family history, or parental consanguinity.

Genetic and pathological links between Parkinson's disease and the lysosomal disorder Sanfilippo syndrome.

BACKGROUND: Parkinson's disease (PD) is a common neurodegenerative disorder of unknown etiology. The characteristic α-synuclein aggregation of PD is also a feature of Sanfilippo syndrome, a storage disorder caused by α-N-acetylglucosaminidase (NAGLU) gene mutations. We explored genetic links between these disorders and studied the pathology of Sanfilippo syndrome to investigate a common pathway toward α-synuclein aggregation. METHODS: We typed the 2 single-nucleotide polymorphisms that tag the common haplotypes of NAGLU in 926 PD patients and 2308 controls and also stained cortical tissue from 2 cases of Sanfilippo A syndrome using the anti-α-synuclein antibody, Per7. RESULTS: Allelic analysis showed an association between rs2071046 and risk for PD (P 1.3 × 10(-3) ). Intracellular α-synuclein accumulation was observed in the cortical tissue of both Sanfilippo A syndrome cases. CONCLUSIONS: This study suggests a possible role of NAGLU in susceptibility to PD while extending evidence for α-synuclein aggregation in the brain in lysosomal storage disorders. Our findings support a mechanism involving lysosomal dysfunction more generally in the pathogenesis of PD.

Interrogating the complex role of chromosome 16p13.13 in multiple sclerosis susceptibility: independent genetic signals in the CIITA-CLEC16A-SOCS1 gene complex.

Multiple sclerosis (MS) is a neurodegenerative, autoimmune disease of the central nervous system, and numerous studies have shown that MS has a strong genetic component. Independent studies to identify MS-associated genes have often indicated multiple signals in physically close genomic regions, although by their proximity it is not always clear if these data indicate redundant or truly independent genetic signals. Recently, three MS study samples were genotyped in parallel using an Illumina Custom BeadChip. These revealed multiple significantly associated single-nucleotide polymorphisms within a 600 kb stretch on chromosome 16p13. Here we present a detailed analysis of variants in this region that clarifies the independent nature of these signals. The linkage disequilibrium patterns in the region and logistic regression analysis of the associations suggest that this region likely harbors three independent MS disease loci. Further, we examined cis-expression QTLs, histone modifications and CCCTC-binding factor (CTCF) binding data in the region. We also tested for correlated expression of the genes from the region using whole-genome expression array data from lymphoblastoid cell lines. Three of the genes show expression correlations across loci. Furthermore, in the GM12878 lymphoblastoid cell line, these three genes are in a continuous region devoid of H3K27 methylation, suggesting an open chromatin configuration. This region likely only contributes minimal risk to MS; however, investigation of this region will undoubtedly provide insight into the functional mechanisms of these genes. These data highlight the importance of taking a closer look at the expression and function of chromosome 16p13 in the pathogenesis of MS.

Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes.

Metformin is the most commonly used pharmacological therapy for type 2 diabetes. We report a genome-wide association study for glycemic response to metformin in 1,024 Scottish individuals with type 2 diabetes with replication in two cohorts including 1,783 Scottish individuals and 1,113 individuals from the UK Prospective Diabetes Study. In a combined meta-analysis, we identified a SNP, rs11212617, associated with treatment success (n = 3,920, P = 2.9 × 10(-9), odds ratio = 1.35, 95% CI 1.22-1.49) at a locus containing ATM, the ataxia telangiectasia mutated gene. In a rat hepatoma cell line, inhibition of ATM with KU-55933 attenuated the phosphorylation and activation of AMP-activated protein kinase in response to metformin. We conclude that ATM, a gene known to be involved in DNA repair and cell cycle control, plays a role in the effect of metformin upstream of AMP-activated protein kinase, and variation in this gene alters glycemic response to metformin.

Dissection of the genetics of Parkinson's disease identifies an additional association 5' of SNCA and multiple associated haplotypes at 17q21.

We performed a genome-wide association study (GWAS) in 1705 Parkinson's disease (PD) UK patients and 5175 UK controls, the largest sample size so far for a PD GWAS. Replication was attempted in an additional cohort of 1039 French PD cases and 1984 controls for the 27 regions showing the strongest evidence of association (P< 10(-4)). We replicated published associations in the 4q22/SNCA and 17q21/MAPT chromosome regions (P< 10(-10)) and found evidence for an additional independent association in 4q22/SNCA. A detailed analysis of the haplotype structure at 17q21 showed that there are three separate risk groups within this region. We found weak but consistent evidence of association for common variants located in three previously published associated regions (4p15/BST1, 4p16/GAK and 1q32/PARK16). We found no support for the previously reported SNP association in 12q12/LRRK2. We also found an association of the two SNPs in 4q22/SNCA with the age of onset of the disease.