Skip to main content
Download PDF

Genomewide association study for onset age in Parkinson disease

BMC Medical Genetics volume 10, Article number: 98 (2009) Cite this article

Abstract

Background

Age at onset in Parkinson disease (PD) is a highly heritable quantitative trait for which a significant genetic influence is supported by multiple segregation analyses. Because genes associated with onset age may represent invaluable therapeutic targets to delay the disease, we sought to identify such genetic modifiers using a genomewide association study in familial PD. There have been previous genomewide association studies (GWAS) to identify genes influencing PD susceptibility, but this is the first to identify genes contributing to the variation in onset age.

Methods

Initial analyses were performed using genotypes generated with the Illumina HumanCNV370Duo array in a sample of 857 unrelated, familial PD cases. Subsequently, a meta-analysis of imputed SNPs was performed combining the familial PD data with that from a previous GWAS of 440 idiopathic PD cases. The SNPs from the meta-analysis with the lowest p-values and consistency in the direction of effect for onset age were then genotyped in a replication sample of 747 idiopathic PD cases from the Parkinson Institute Biobank of Milan, Italy.

Results

Meta-analysis across the three studies detected consistent association (p < 1 × 10-5) with five SNPs, none of which reached genomewide significance. On chromosome 11, the SNP with the lowest p-value (rs10767971; p = 5.4 × 10-7) lies between the genes QSER1 and PRRG4. Near the PARK3 linkage region on chromosome 2p13, association was observed with a SNP (rs7577851; p = 8.7 × 10-6) which lies in an intron of the AAK1 gene. This gene is closely related to GAK, identified as a possible PD susceptibility gene in the GWAS of the familial PD cases.

Conclusion

Taken together, these results suggest an influence of genes involved in endocytosis and lysosomal sorting in PD pathogenesis.

Peer Review reports

Background

Parkinson disease (PD), the second most common neurodegenerative disorder, is characterized by debilitating symptoms of tremor, rigidity, and bradykinesia, usually occurring late in life. PD incidence increases with age from 1.7/10,000 person-years between ages 50 to 59 to 9.3/10,000 person-years between ages 70 to 79 and has a prevalence of approximately 1.8% among people over the age of 65 [1]. While the average age of onset of PD is approximately 60 years, there is wide variation, with some individuals having onset before age 20 and others not until after age 90 [2, 3].

Five monogenic forms of PD have been identified and characterized by mode of transmission, onset age and protein affected by mutation. These include α-synuclein (SNCA or PARK1) [4], parkin (PARK2) [5], PTEN-induced putative kinase 1 (PINK1 or PARK6) [6], DJ-1 (PARK7) [7], and leucine-rich repeat kinase 2 (LRRK2 or PARK8) [8]. Onset for PARK1 is younger than that seen for idiopathic PD [4]. PARK2 (parkin) is a recessive form with young onset, commonly before age 40. Heterozygous mutations in parkin are also associated with earlier onset of PD, typically in the early to mid sixth decade [9, 10]. By contrast, PD associated with LRRK2 mutations presents an onset distribution very similar to that seen in idiopathic PD, as well as clear age-dependent penetrance [11–13].

Onset of PD has been shown to be correlated between siblings with PD [14] suggesting that genetic modifiers influence onset age. Segregation analyses in three independent studies showed evidence of a genetic effect influencing age of onset of PD [15–17]. Notably, all three of these segregation analyses showed stronger evidence for the presence of "major genes" influencing onset age or penetrance, than for genes influencing susceptibility. Furthermore, age is one of the strongest risk factors for PD, suggesting that age related penetrance is strongly associated with disease expression. By identifying genes related to onset age, it may be possible to identify pathogenic mechanisms and therapeutic targets capable of delaying onset of disease symptoms. Effectively postponing disease onset will reduce disease prevalence and ease the burden of PD in our aging population.

All prior PD genome wide association studies (GWAS) have focused exclusively on the detection of susceptibility genes and none has investigated association to genes influencing onset age [18–20]. In this study, we describe the first GWAS of onset age. This GWAS included 857 PD cases with a positive family history of PD. In addition, we performed a second GWAS with onset age as the phenotype using publicly available data from 440 randomly ascertained PD cases [19]. We conducted a meta-analysis of the two studies comprising approximately 2 million SNPs imputed using HapMap data. Finally, a replication study of the top findings from the meta-analysis was performed in an independent sample of 747 randomly ascertained PD cases from Milan, Italy.

Methods

PD Cases

One PD case (n = 935) from each family recruited from two ongoing studies of familial PD, the GenePD study and the PROGENI study, was selected for the GWAS. Both studies recruited families consisting of at least two members meeting diagnostic criteria for PD. PD cases underwent a uniform neurological evaluation that employed PD diagnostic criteria based on a modified version of the United Kingdom PD Society Brain Bank Criteria [21]. Detailed descriptions of the inclusion and exclusion criteria are described elsewhere for the PROGENI [22] and GenePD [14] studies. The analyzed sample was exclusively white, non-Hispanic.

No subject known to have a disease-producing mutation was included in this analysis. All cases were known to be negative for the LRRK2 G2019S mutation, and many, but not all, were also screened for PARK1(SNCA)(N screened = 702), PARK2 (parkin)(N = 593), PARK7 (DJ1)(N = 328), and NR4A2 (N = 550) [9, 10, 12, 23–27]. PD onset age was determined by interview and reflected the age of first symptom of PD, which commonly preceded age of physician diagnosis. The reliability of ascertaining age of onset though interview compared to medical records has been estimated as 0.94 [28].

In addition, control samples (n = 895) obtained from the NINDS Human Genetics Resource Center DNA and Cell Line Repository (Camden, NJ) were intermixed on the same plates as the cases when genotyped. All selected control samples were reported to be white, non-Hispanic. While not used in the analysis of onset age of PD, these samples were used for SNP and sample quality assessment and genotype imputation as described below.

Appropriate written informed consent was obtained for all samples included in this study.

Microarray Genotyping and Quality Assessment

Genotyping was performed by the Center for Inherited Disease Research (CIDR) using the Illumina HumanCNV370 version1_C BeadChips (Illumina, San Diego, CA, USA) and the Illumina Infinium II assay protocol [29]. As previously described [20], 78 cases and 28 controls were removed due to a low genotype call rate (<98% of SNPs), cryptic relatives in the sample, or population stratification. The final sample included 857 PD cases (all derived from whole blood) and 867 control DNA samples (all derived from lymphoblastoid cell lines). Genotype calls and quality scores were determined from allele cluster definitions for each SNP as determined by the Illumina BeadStudio Genotyping Module version 3.1.14 and the combined intensity data from 96% of study samples as previously described [20]. Genotype calls with a quality score (Gencall value) of 0.25 or higher were considered acceptable. Blind duplicate reproducibility was 99.98%. Although we performed additional SNP filtering described below, the CIDR released dataset contained 344,301 SNPs and is available on dbGaP (http://www.ncbi.nlm.nih.gov/sites/entrez?db=gap; Accession number: phs000126.v1.p1).

SNPs were removed if: 1) the call rate of the SNPs was lower than 98% (n = 7,764), 2) the minor allele frequency was less than 0.01 in the combined case and control dataset (n = 7,667), 3) there were differential rates of missing genotypes in the cases and controls (n = 75) or males and females (n = 271), or 4) significant deviation from Hardy Weinberg equilibrium was observed in the control sample (n = 906). Many markers failed multiple tests. The final dataset consisted of 328,189 SNPs that passed all quality control measures (94.6% of all attempted SNPs).

Mayo-Perlegen LEAPS sample

Publicly available SNP genotype data for 443 PD cases was accessed from dbGaP for the GWAS "Mayo-Perlegen LEAPS (Linked Efforts to Accelerate Parkinson's Solutions) Collaboration" (http://www.ncbi.nlm.nih.gov/sites/entrez?db=gap; Accession number: phs000048.v1.p1; [19]. This data set included genotyping for the 198,345 SNPs meeting the quality control standards described by Maraganore et al. [19]. We examined these data for population stratification and for cryptic relatedness. While no population outliers were observed, three individuals were removed due to apparent second-degree relationships identified during our review of the data, leaving a sample of 440 PD cases included in these analyses. While data are also publicly available for 267 PD cases from the Fung et al. GWAS [18], that study only included cases whose age of PD onset was 55 years or greater resulting in an age distribution that was significantly different from the other two samples (p < 0.0001). Therefore, we did not include the cases from that study in our meta-analysis due to limited variability in the onset distribution.

Population stratification

Both the GenePD-PROGENI and Mayo-Perlegen LEAPS sample sets were screened for population outliers during initial QC. In addition, after the final study sample sets were determined, principle components were recalculated independently for each study using only the samples included in the final analyses. Association to onset age was then tested for the first six principle components for each study. No association (p < 0.05) between the first six principle components and onset age was seen in either study sample.

Imputation

Imputation was performed to increase the power of the meta-analysis and to facilitate the joint analysis of results generated from Mayo-Perlegen LEAPS and our GWAS, which were genotyped on different platforms. The program MACH 1.0 (compiled using Intel's optimized compiler) was used to impute genotypes for 2,543,887 autosomal SNPs characterized in the HapMap project [30–32].

The GenePD-PROGENI sample and the Mayo-Perlegen LEAPS cohort were imputed separately using phased haplotype data downloaded from the HapMap project website [33].

In the GenePD-PROGENI sample, imputation was performed using both cases and controls. A subset of 200 cases and controls with high call rates was selected to perform the initial model parameter calculation. Next, imputation was performed on all participants using all autosomal SNPs where the strand was not ambiguous (i.e. not an A/T or G/C SNP) and that passed all other quality control measures described previously. In the Mayo-Perlegen LEAPS cohort, all PD cases were used to compute the initial model parameters. All unambiguous autosomal SNPs passing quality review, as described in the original Mayo-Perlegen LEAPS GWAS [19], were used for this imputation. In addition, because the genotyping platform used in the Mayo-Perlegen LEAPS study included a high percentage of ambiguous SNPs, those ambiguous SNPs which could be confidently matched to HapMap strands through the comparison of minor allele frequencies were also used for imputation. For each study, initial model parameters were calculated using 100 iterations. SNP quality was assessed using the Rsq metric, which estimates the squared correlation between imputed and actual genotypes with SNPs having an Rsq <0.3 excluded from further study (N = 61,271).

At the time of this study, only autosomal imputation is supported by MACH; therefore, the program IMPUTE was used to impute SNPs on the X chromosome [34, 35]http://www.stats.ox.ac.uk/~marchini/software/gwas/impute.html. Haplotype, legend, and recombination rate files based upon HapMap (rel#21-NCBI build 35) were downloaded from IMPUTE's website for use in imputing 64,621 SNPs in the non-pseudoautosomal region of the X chromosome. Imputation was run in all participants in the GenePD-PROGENI study and separately in all PD cases in the Mayo-Perlegen LEAPS GWAS, using all X-chromosome SNPs that passed quality assessment as described above and were found in the legend file provided by IMPUTE.

Statistical Analyses

In order to meta-analyze results from both studies, association to onset age was performed using exclusively imputed, not genotyped data, for 1,861,750 SNPs passing imputation QC and with minor allele frequencies greater than 10%. Imputed SNPs with minor allele frequencies less than 10% in the GenePD-PROGENI sample (N = 620,866) were excluded from the association analyses to avoid false positives that can occur in the analysis of low minor allele frequency SNPs. To evaluate SNP association under an additive mode of inheritance, the predicted allele dosage for each genotype estimated by MACH for autosomes or by IMPUTE for the X chromosome was used. To model the recessive and dominant modes of inheritance, the genotype probabilities calculated by MACH were used. While a combination of additive and recessive models can have high power to detect dominant effects, this power decreases as the minor allele frequency of the SNP nears 0.5 [36]. Therefore, we studied additive, dominant and recessive models. The total probability of having either one copy or two copies of the minor allele was used for the dominant model and the probability of having two copies of the minor allele was used for the recessive model. Recessive and dominant models were not studied for the X chromosome. Linear regression analyses were performed using SAS v9.1.

Meta-Analysis

Since the genotypic data was generated on different arrays (Illumina and Perlegen) with few SNPs in common, we employed a conservative meta-analytic approach to combine results from the two studies. Meta-analysis of the results of the linear regression of the imputed data for the GenePD-PROGENI sample and the Mayo-Perlegen LEAPS cohort was performed using METAL http://www.sph.umich.edu/csg/abecasis/metal/. As is common in GWAS meta-analyses, a fixed effects model with standard error weighting was used, as random-effects models may be too conservative in GWA studies with a small number of studies[37].

Replication Study

In order to validate the top findings from the GWAS meta-analysis, an additional sample of 896 PD cases with reported ages of PD onset was provided by the Parkinson Institute - Istituti Clinici di Perfezionamento, Milan, Italy from the "Human Genetic Bank of Patients Affected by PD and Parkinsonisms". These cases were recruited irrespective of family history or onset age, and similar to the GenePD-PROGENI, used the UK Parkinson's Disease Society Brain Bank criteria to confirm idiopathic PD [38] and defined onset as the age of first symptom of PD.

Twenty-four SNPs were selected based upon the following criteria (1) a p-value less than 0.00001 in the meta-analysis of the two GWAS, with (2) a consistent direction of effect in both studies. Nine SNPs meeting these criteria were identified from the additive inheritance model, thirteen from the dominant model and eleven from the recessive model. Four SNPs were identified in both the additive and dominant model, one SNP was identified in both the additive and recessive model and one SNP was identified by all three models. Three gene regions were identified under two different genetic models with different SNPs. For each multiply-nominated gene region, the SNP from the model with the smaller p-value was selected for replication. These SNPs were genotyped using TaqMan technology implemented on the ABI PRISM® 7900HT Sequence Detection system (Applied Biosystems: Foster City, CA) at Boston University School of Medicine. Individual samples (149) with genotyping call rates of less than 95% were excluded from further analysis.

Association of the 24 SNPs to onset age was evaluated in the 747 remaining Italian cases using linear regression performed with the software Plink v1.01 http://pngu.mgh.harvard.edu/purcell/plink/[39] using the corresponding genetic model (additive, recessive or dominant) by which each SNP was originally identified. A final fixed effects meta-analysis of all three studies was performed using METAL.

To distinguish whether associations observed were to age in general, as opposed to age at onset of PD, linear regression to censoring age was performed in the 867 NINDS control samples genotyped with the GenePD-PROGENI cases.

Results

Demographic characteristics of the three samples studied are shown in Table 1. All three studies have a similar percentage of male participants. The GenePD-PROGENI and Mayo-Perlegen LEAPS samples have similar mean ages of PD onset, while the Italian sample has a somewhat younger average age at onset. The GenePD-PROGENI sample has the widest range of onset ages from 19 to 90 years while the Mayo-Perlegen LEAPS has no participants under 30 years of age and the Italian sample has no participants over 81 years of age. No significant differences in onset age are seen between men and women in any of the studies.

Table 1 Age and sex distribution of the GenePD- PROGENI, Mayo-Perlegen LEAPS, and Milan, Italian study samples.
Full size table

Supplementary Tables S1 (additive model), S2 (dominant model) and S3 (recessive model) (see Additional file 1) present the top SNPs from each region with a meta-analysis p < 0.0001 for imputed SNP data in the GenePD-PROGENI and Mayo-Perlegen LEAPS studies. Twenty-four SNPs with meta-analysis p-values of less than 0.00001 and with a consistent direction of effect for both GWAS were genotyped in the Italian replication sample of 747 PD cases (Table 2). SNPs genotyped in either the GenePD-PROGENI or Mayo-Perlegen LEAPS GWAS platforms are distinguished by notation from those SNPs imputed by both studies. The results of the association analysis in the Italian sample, as well as the combined meta-analysis of all three samples, are shown in Tables 3 (additive), 4 (dominant), and 5 (recessive) and in Figure 1.

Table 2 Location and minor allele frequency (MAF) information for 24 SNPs selected for replication study in Italian sample
Full size table
Table 3 Additive model study-specific and meta-analysis results
Full size table
Table 4 Dominant model study-specific and meta-analysis results
Full size table
Table 5 Recessive model study-specific and meta-analysis results
Full size table
Figure 1
figure 1

Forest plots showing study-specific and pooled effects for top six results in final meta-analysis. The effect (in years) on age at PD onset and 95% confidence intervals for the indicated genetic model are shown for each of the three independent samples and the combined analysis. A. Additive model of the A allele in SNP rs17565841 B. Dominant model of the A allele in SNP rs17565841 C. Dominant model of the C allele in SNP rs1941184 D. Dominant model of the A allele in SNP rs10918270 E. Recessive model of the T allele in SNP rs10767971 F. Recessive model of the T allele in SNP rs7577851.

Full size image

Although ten of the replication SNPs showed a consistent direction of effect across all three studies, only two SNPs (both of which were genotyped on the Mayo-Perlegen LEAPS GWAS platforms) resulted in increased statistical evidence of association when combining the three studies, and both of these were in a recessive model. The most strongly associated SNP, rs10767971, located on chromosome 11 between the genes QSER1 and PRRG4, was associated with a 3.2 year older PD onset in individuals with 2 copies of the minor allele (p = 5.4 × 10-7 in the 3 sample meta-analysis, compared to 4.3 × 10-6 in the 2 sample meta-analysis). Conversely, an estimated 6.9 year earlier age of onset (p = 8.7 × 10-6 in the 3 sample meta-analysis compared to 9.8 × 10-6 in the 2 sample meta-analysis) was observed for the SNP rs7577851, located in the 16th intron of the gene AAK1 on chromosome 2.

The most highly associated SNP in both the additive and dominant three sample meta-analyses, rs17565841, is located approximately 3 kb from the 3' end of the gene OCA2 on chromosome 15. This SNP was associated with an average 2.8 years younger onset age (p = 2.6 × 10-6) under an additive model and a 3.3 years younger onset age (p = 2.1 × 10-6) under a dominant model. Inclusion of the Italian cases in the meta-analysis did not strengthen the evidence of association as compared to the results seen in the two sample meta-analyses (9.1 × 10-7 for additive and 1.9 × 10-6 for dominant). However, in the Italian replication sample, this SNP did provide modest statistical association to onset age (p = 0.05 for additive and p = 0.04 for dominant) with the same direction of effect seen in the two other studies.

Also showing consistent directions of effect and p-values in the three-sample meta-analysis at the level of p < 1 × 10-5 were two SNPs located in the genes DSG3 and ATF6. The SNP rs1941184, had the second best p-value identified under the dominant model and is located in the third intron of the gene DSG3 on chromosome 18. This SNP was associated with an average 2.3 year younger age of onset of PD across the three studies (p = 4.3 × 10-6). The SNP rs10918270, located in the 15th intron of the gene ATF6 on chromosome 1, was identified under both additive and dominant modes of inheritance, but showed stronger association in the three sample meta-analysis under the dominant model (p = 7.5 × 10-6) with an average 2.3 year younger onset of PD. No association to age with direction of effect consistent with that observed for onset age of PD was seen in the control sample for any of the 24 SNPs at a significance level equal to 0.05.

Discussion

We present results from the first GWAS for age at onset of PD, including a meta-analysis with the publicly available Mayo-Perlegen LEAPS GWAS data (dbGaP Study Accession: phs000048.v1.p1) and a follow-up replication study in an independent PD sample recruited in Milan, Italy. Differences were observed in the age distributions in the three populations used in this study (Table 1). However, there were no imposed age restrictions in any of the studies and a wide distribution of ages was represented in all three populations, each with a range of greater than 60 years.

No SNP reached the commonly accepted criterion for genome-wide significance of p < 5 × 10-8 [40], which is based on recent estimates of independent genomewide sequence variation to maintain 5% genomewide type I error rate [41, 42]. While this criterion provides an appropriate cutoff for determining significance for the large number of SNPs provided by imputation, this measure does not account for the testing of multiple genetic models, as was performed in this study. Despite the lack of genomewide significance, the meta-analysis in this study showed evidence of several interesting loci with consistent effects on onset age of PD across the three independent populations studied.

The SNP with the strongest evidence for association to onset age, rs10767971, is associated with a later age of onset under a recessive model. This SNP is nearly equidistant between the genes PRRG4 and QSER1 on chromosome 11, at just under 20 kb from each. The association results for nearby SNPs studied in the meta-analysis of the two GWAS as well as the LD structure and recombination rates for the region are shown in Figure 2A and 2B under recessive and dominant models. While the strongest evidence of association appeared with a recessive model, many nearby SNPs in strong LD showed similar association under the dominant model. The pattern suggests a single association, spanning several SNPs that may be more evident around the PRRG4 gene and may extend toward the CCDC73 and EIF3M genes.

Figure 2
figure 2

Evidence of association in Chromosome 11 region. Along the x axis is the physical position in the region (in kB) with known genes shown in their orientation The left Y-axis denotes the association test result as -log(p-value) corresponding to diamonds in the figure. The blue diamond identifies the primary SNP result labeled with an rs# and p-value. The color of additional diamonds depicts the pairwise linkage disequilibrium with the primary SNP: red indicates r2 > 0.8, orange 0.5 < r2 < 0.8, yellow 0.2 < r2 < 0.5, white r2 < 0.2, grey unknown LD. R2 values were obtained from the CEPH HapMap data. The right Y-axis indicates the recombination rate, also obtained from CEPH HapMap data, and shown within the figure by the solid light-blue line. A. Recessive model B. Dominant model.

Full size image

The SNP showing the strongest association to PD onset across the three populations under both the additive and dominant models lies 3 kb from the 3' end of the gene OCA2 (oculocutaneous albinism II). Mutations in OCA2 are responsible for the most common form of albinism, oculocutaneous albinism type 2 [43] and several studies have shown this gene to be associated with common skin, hair and eye color variation found in European populations [43–48]. Variations in OCA2 have been associated with susceptibility to melanoma [49], which has also been reported to occur with increased frequency among PD cases [50–52]. Because of the targeted degradation of pigmented neurons in PD brains, these associations have lead to a hypothesized link between genes involved in pigmentation, such as OCA2, and PD, likely mediated through common elements in melanin and neuromelanin synthesis[53]. The protein encoded by OCA2, 'P' protein, is involved in the transport of tyrosine, a precursor to melanin, as well as in the regulation of melanosomal pH which may be key to the initiation of the enzyme controlling melanin synthesis in melanocytes [54]. It is not clear whether synthesis of neuromelanin is also regulated through melanosomal pH or is affected by the 'P' protein in a similar way [55]. Nevertheless, the association between OCA2 and younger age of onset are suggestive of a neuromelanin-related mechanism of effect.

The association of an intronic SNP in the gene DSG3 (desmoglein 3 (pemphigus vulgaris antigen)) may also be indicative of a neuromelanin related effect on onset age of PD. The protein encoded by DSG3 is the autoantigen for the autoimmune skin disease pemphigus vulgaris and this gene is expressed primarily in skin. Interestingly, some reports have demonstrated increased expression of DSG3 in melanocytes compared to keratinocytes (the most common cell type in the epidermis) [56].

An intronic SNP in the gene ATF6 also showed strong association to earlier age of onset of PD. ATF6 (activating transcription factor 6) transcribes a transcription factor localized to the endoplasmic reticulum (ER). The ATF6 protein is a critical regulator of the unfolded protein response (UPR), a highly conserved pathway activated in response to ER stress, and is a protective cellular response to the accumulation of misfolded proteins [57–59]. The UPR has been implicated in neurotoxin based cellular models of PD [60] and has also been shown to be activated by the over-expression of α-synuclein in yeast cells [61]. More recently, postmortem studies of PD case and control brains have shown activation of the UPR in cases, but not controls and that this activation is associated with the aggregation of α-synuclein [62].

The finding of an intronic SNP in AAK1 associated with PD onset age is intriguing because of the gene's genomic location, its function, and its close relation to a gene identified for PD susceptibility. The AAK1 (AP2 associated kinase 1) gene is located on chromosome 2p14, near 2p13, which has previously been implicated as the PARK3 locus[63], and for which linkage to onset age was demonstrated in both the GenePD [64] and PROGENI [65] studies. The AAK1 gene itself has not been previously implicated by positional mapping, but a microarray study of PD brain compared to controls demonstrated differential expression of AAK1 [66]. In our past GWAS of PD susceptibility in the GenePD-PROGENI cohort [20], the region containing the gene GAK (cyclin G-associated kinase) had the strongest evidence for association (chromosome 4p). The AAK1 and GAK genes both function at multiple steps in clathrin-mediated vesicular transport and the two kinases likely have some redundant functions [67] related to their homologous serine/threonine-kinase domain [68]. Recently, cathepsin D was implicated as the main lysosomal enzyme involved in α-synuclein degredation [69], and depletion of GAK was shown to impair the lysosomal sorting of cathepsin D [68]. Thus, the finding that AAK1 influences onset age and GAK influences risk in familial PD suggest that pathways involving lysosomal activity influence PD risk.

Several previous genome scans have provided evidence of loci linked to onset age, including 2p13 seen in both the GenePD and PROGENI studies, as mentioned above. Evidence of linkage has also been reported on chromosomes 1p, 1q, 8q, 9q, 10q, 20 and 21[64, 65, 70], but there is little overlap between these linkage regions and our top meta-analysis association results. Aside from the identification of the AAK1 gene near 2p13, the strongest association result observed under one of the previously reported linkage regions occurs near the gene sortilin-related VPS10 domain containing receptor 3 (SORCS3) located near the LOD score peak at 10q, originally identified in a combined linkage scan of onset age in PD and Alzheimer's disease [70]. This association is seen most strongly under a dominant model (p = 2.5 × 10-5) with minor allele carriers having an estimated older onset age by 3.3 years (see Additional file 1: Table S2).

A final region of interest is 15q26.2 that includes the gene MCTP2. Several SNPs in this region showed association to earlier onset age of PD in the meta-analysis of the two GWAS under the recessive model (see Additional file 1: Table S3). These SNPs did not reach the criteria for inclusion in the replication study, as the strongest p-value seen was 2.2 × 10-5 with rs17504636. This SNP was associated with a 9.2 year earlier average PD onset. This region overlaps with a SNP reported in the susceptibility GWAS including these cases [20]. In that study the SNP rs4476132 was associated with PD susceptibility with an odds ratio of 1.3 (p = 7.7 × 10-5- meta-analysis with additive model). This overlap is consistent with a locus associated with a risk for younger-onset PD or with an effect modifying age dependent penetrance. The gene in this region, MCTP2 (multiple C2 domains, transmembrane 2), is expressed in the brain and has been implicated in linkage and association studies of abdominal fat [71] and major depression [72].

Important distinctions can be made between those genes that influence susceptibility for developing disease, and the genetic modifiers that influence penetrance or, as studied here, onset age. Perhaps the best examples for genetic modifiers are seen for Huntington's disease (HD) where an expanded CAG trinucleotide repeat on chromosome 4p16.3 causes the disease, but wide variation in onset age is evident for individuals with identical repeat lengths. The identification of the genes that presumably interact with huntingtin to produce relatively younger or older onset for a given repeat size provide insight into the pathogenic mechanisms for HD, as well as therapeutic targets for intervention [73–75]. Similarly, identifying those genes and their products that are associated with older onset in PD may provide insight into the disease mechanisms and processes for delaying onset with implications for novel treatments.

Conclusion

The identification of the 15q26.2 region as well as the related genes AAK1 and GAK in the studies of PD onset age and susceptibility highlights the importance of the continued study of both of these traits, both separately and in combination. The direct overlap in affection and onset age association results in the 15q26.2 region shows this to be a candidate region that would benefit from further examination with consideration of important age-related effects, for example in studies correlating expression of genes in this region with onset age. The identification of association to onset age with the gene AAK1, in the same pathway as a previously identified susceptibility-associated gene GAK highlights the importance of genetic pathways in PD etiology, showing that the genes along the same pathway may have redundant effects or may modify disease pathology different ways, observed by differences in disease onset and progression. Studying PD in the context of onset age provides fundamental insight into the disease process and is essential to understanding mechanisms that modify disease penetrance and therefore may be key in identifying therapeutic targets.

References

  1. Mayeux R: Epidemiology of neurodegeneration. Annual review of neuroscience. 2003, 26: 81-104. 10.1146/annurev.neuro.26.043002.094919.

    Article  CAS  PubMed  Google Scholar 

  2. Schrag A, Ben-Shlomo Y, Brown R, Marsden CD, Quinn N: Young-onset Parkinson's disease revisited--clinical features, natural history, and mortality. Mov Disord. 1998, 13 (6): 885-894. 10.1002/mds.870130605.

    Article  CAS  PubMed  Google Scholar 

  3. Bower JH, Maraganore DM, McDonnell SK, Rocca WA: Incidence and distribution of parkinsonism in Olmsted County, Minnesota, 1976-1990. Neurology. 1999, 52 (6): 1214-1220.

    Article  CAS  PubMed  Google Scholar 

  4. Polymeropoulos MH, Lavedan C, Leroy E, Ide SE, Dehejia A, Dutra A, Pike B, Root H, Rubenstein J, Boyer R, Stenroos ES, Chandrasekharappa S, Athanassiadou A, Papapetropoulos T, Johnson WG, Lazzarini AM, Duvoisin RC, Di Iorio G, Golbe LI, Nussbaum RL: Mutation in the alpha-synuclein gene identified in families with Parkinson's disease. Science (New York, NY). 1997, 276 (5321): 2045-2047.

    Article  CAS  Google Scholar 

  5. Kitada T, Asakawa S, Hattori N, Matsumine H, Yamamura Y, Minoshima S, Yokochi M, Mizuno Y, Shimizu N: Mutations in the parkin gene cause autosomal recessive juvenile parkinsonism. Nature. 1998, 392 (6676): 605-608. 10.1038/33416.

    Article  CAS  PubMed  Google Scholar 

  6. Valente EM, Bentivoglio AR, Dixon PH, Ferraris A, Ialongo T, Frontali M, Albanese A, Wood NW: Localization of a novel locus for autosomal recessive early-onset parkinsonism, PARK6, on human chromosome 1p35-p36. Am J Hum Genet. 2001, 68 (4): 895-900. 10.1086/319522.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. van Duijn CM, Dekker MC, Bonifati V, Galjaard RJ, Houwing-Duistermaat JJ, Snijders PJ, Testers L, Breedveld GJ, Horstink M, Sandkuijl LA, van Swieten JC, Oostra BA, Heutink P: Park7, a novel locus for autosomal recessive early-onset parkinsonism, on chromosome 1p36. Am J Hum Genet. 2001, 69 (3): 629-634. 10.1086/322996.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Paisan-Ruiz C, Jain S, Evans EW, Gilks WP, Simon J, Brug van der M, Lopez de Munain A, Aparicio S, Gil AM, Khan N, Johnson J, Martinez JR, Nicholl D, Carrera IM, Pena AS, de Silva R, Lees A, Marti-Masso JF, Perez-Tur J, Wood NW, Singleton AB: Cloning of the gene containing mutations that cause PARK8-linked Parkinson's disease. Neuron. 2004, 44 (4): 595-600. 10.1016/j.neuron.2004.10.023.

    Article  CAS  PubMed  Google Scholar 

  9. Sun M, Latourelle JC, Wooten GF, Lew MF, Klein C, Shill HA, Golbe LI, Mark MH, Racette BA, Perlmutter JS, Parsian A, Guttman M, Nicholson G, Xu G, Wilk JB, Saint-Hilaire MH, DeStefano AL, Prakash R, Williamson S, Suchowersky O, Labelle N, Growdon JH, Singer C, Watts RL, Goldwurm S, Pezzoli G, Baker KB, Pramstaller PP, Burn DJ, Chinnery PF, Sherman S, Vieregge P, Litvan I, Gillis T, MacDonald ME, Myers RH, Gusella JF: Influence of heterozygosity for parkin mutation on onset age in familial Parkinson disease: the GenePD study. Arch Neurol. 2006, 63 (6): 826-832. 10.1001/archneur.63.6.826.

    Article  PubMed  Google Scholar 

  10. Foroud T, Uniacke SK, Liu L, Pankratz N, Rudolph A, Halter C, Shults C, Marder K, Conneally PM, Nichols WC: Heterozygosity for a mutation in the parkin gene leads to later onset Parkinson disease. Neurology. 2003, 60 (5): 796-801.

    Article  CAS  PubMed  Google Scholar 

  11. Nichols WC, Pankratz N, Hernandez D, Paisan-Ruiz C, Jain S, Halter CA, Michaels VE, Reed T, Rudolph A, Shults CW, Singleton A, Foroud T: Genetic screening for a single common LRRK2 mutation in familial Parkinson's disease. Lancet. 2005, 365 (9457): 410-412.

    CAS  PubMed  Google Scholar 

  12. Latourelle JC, Sun M, Lew MF, Suchowersky O, Klein C, Golbe LI, Mark MH, Growdon JH, Wooten GF, Watts RL, Guttman M, Racette BA, Perlmutter JS, Ahmed A, Shill HA, Singer C, Goldwurm S, Pezzoli G, Zini M, Saint-Hilaire MH, Hendricks AE, Williamson S, Nagle MW, Wilk JB, Massood T, Huskey KW, Laramie JM, Destefano AL, Baker KB, Itin I, Litvan I, Nicholson G, Corbett A, Nance M, Drasby E, Isaacson S, Burn DJ, Chinnery PF, Pramstaller PP, Al-Hinti J, Moller AT, Ostergaard K, Sherman SJ, Roxburgh R, Snow B, Slevin JT, Cambi F, Gusella JF, Myers RH: The Gly2019Ser mutation in LRRK2 is not fully penetrant in familial Parkinson's Disease: the GenePD study. BMC Med. 2008, 6 (1): 32-10.1186/1741-7015-6-32.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Goldwurm S, Zini M, Mariani L, Tesei S, Miceli R, Sironi F, Clementi M, Bonifati V, Pezzoli G: Evaluation of LRRK2 G2019S penetrance: relevance for genetic counseling in Parkinson disease. Neurology. 2007, 68 (14): 1141-1143. 10.1212/01.wnl.0000254483.19854.ef.

    Article  CAS  PubMed  Google Scholar 

  14. Maher NE, Golbe LI, Lazzarini AM, Mark MH, Currie LJ, Wooten GF, Saint-Hilaire M, Wilk JB, Volcjak J, Maher JE, Feldman RG, Guttman M, Lew M, Waters CH, Schuman S, Suchowersky O, Lafontaine AL, Labelle N, Vieregge P, Pramstaller PP, Klein C, Hubble J, Reider C, Growdon J, Watts R, Montgomery E, Baker K, Singer C, Stacy M, Myers RH: Epidemiologic study of 203 sibling pairs with Parkinson's disease: the GenePD study. Neurology. 2002, 58 (1): 79-84.

    Article  CAS  PubMed  Google Scholar 

  15. Maher NE, Currie LJ, Lazzarini AM, Wilk JB, Taylor CA, Saint-Hilaire MH, Feldman RG, Golbe LI, Wooten GF, Myers RH: Segregation analysis of Parkinson disease revealing evidence for a major causative gene. American journal of medical genetics. 2002, 109 (3): 191-197. 10.1002/ajmg.10335.

    Article  CAS  PubMed  Google Scholar 

  16. McDonnell SK, Schaid DJ, Elbaz A, Strain KJ, Bower JH, Ahlskog JE, Maraganore DM, Rocca WA: Complex segregation analysis of Parkinson's disease: The Mayo Clinic Family Study. Annals of neurology. 2006, 59 (5): 788-795. 10.1002/ana.20844.

    Article  CAS  PubMed  Google Scholar 

  17. Zareparsi S, Taylor TD, Harris EL, Payami H: Segregation analysis of Parkinson disease. American journal of medical genetics. 1998, 80 (4): 410-417. 10.1002/(SICI)1096-8628(19981204)80:4<410::AID-AJMG21>3.0.CO;2-2.

    Article  CAS  PubMed  Google Scholar 

  18. Fung HC, Scholz S, Matarin M, Simon-Sanchez J, Hernandez D, Britton A, Gibbs JR, Langefeld C, Stiegert ML, Schymick J, Okun MS, Mandel RJ, Fernandez HH, Foote KD, Rodriguez RL, Peckham E, De Vrieze FW, Gwinn-Hardy K, Hardy JA, Singleton A: Genome-wide genotyping in Parkinson's disease and neurologically normal controls: first stage analysis and public release of data. Lancet neurology. 2006, 5 (11): 911-916. 10.1016/S1474-4422(06)70578-6.

    Article  CAS  PubMed  Google Scholar 

  19. Maraganore DM, de Andrade M, Lesnick TG, Strain KJ, Farrer MJ, Rocca WA, Pant PV, Frazer KA, Cox DR, Ballinger DG: High-resolution whole-genome association study of Parkinson disease. American journal of human genetics. 2005, 77 (5): 685-693. 10.1086/496902.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Pankratz N, Wilk JB, Latourelle JC, DeStefano AL, Halter C, Pugh EW, Doheny KF, Gusella JF, Nichols WC, Foroud T, Myers RH: Genomewide association study for susceptibility genes contributing to familial Parkinson disease. Human genetics. 2009, 124 (6): 593-605. 10.1007/s00439-008-0582-9.

    Article  CAS  PubMed  Google Scholar 

  21. Gibb WR, Lees AJ: The relevance of the Lewy body to the pathogenesis of idiopathic Parkinson's disease. Journal of neurology, neurosurgery, and psychiatry. 1988, 51 (6): 745-752. 10.1136/jnnp.51.6.745.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Pankratz N, Nichols WC, Uniacke SK, Halter C, Rudolph A, Shults C, Conneally PM, Foroud T: Genome screen to identify susceptibility genes for Parkinson disease in a sample without parkin mutations. American journal of human genetics. 2002, 71 (1): 124-135. 10.1086/341282.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Pankratz N, Pauciulo MW, Elsaesser VE, Marek DK, Halter CA, Rudolph A, Shults CW, Foroud T, Nichols WC: Mutations in LRRK2 other than G2019S are rare in a north American-based sample of familial Parkinson's disease. Mov Disord. 2006, 21 (12): 2257-2260. 10.1002/mds.21162.

    Article  PubMed  Google Scholar 

  24. Nichols WC, Elsaesser VE, Pankratz N, Pauciulo MW, Marek DK, Halter CA, Rudolph A, Shults CW, Foroud T: LRRK2 mutation analysis in Parkinson disease families with evidence of linkage to PARK8. Neurology. 2007, 69 (18): 1737-1744. 10.1212/01.wnl.0000278115.50741.4e.

    Article  CAS  PubMed  Google Scholar 

  25. Nichols WC, Uniacke SK, Pankratz N, Reed T, Simon DK, Halter C, Rudolph A, Shults CW, Conneally PM, Foroud T: Evaluation of the role of Nurr1 in a large sample of familial Parkinson's disease. Mov Disord. 2004, 19 (6): 649-655. 10.1002/mds.20097.

    Article  PubMed  Google Scholar 

  26. Karamohamed S, Golbe LI, Mark MH, Lazzarini AM, Suchowersky O, Labelle N, Guttman M, Currie LJ, Wooten GF, Stacy M, Saint-Hilaire M, Feldman RG, Liu J, Shoemaker CM, Wilk JB, DeStefano AL, Latourelle JC, Xu G, Watts R, Growdon J, Lew M, Waters C, Vieregge P, Pramstaller PP, Klein C, Racette BA, Perlmutter JS, Parsian A, Singer C, Montgomery E, Baker K, Gusella JF, Herbert A, Myers RH: Absence of previously reported variants in the SCNA (G88C and G209A), NR4A2 (T291D and T245G) and the DJ-1 (T497C) genes in familial Parkinson's disease from the GenePD study. Mov Disord. 2005, 20 (9): 1188-1191. 10.1002/mds.20515.

    Article  PubMed  Google Scholar 

  27. Pankratz N, Nichols WC, Elsaesser VE, Pauciulo MW, Marek DK, Halter CA, Wojcieszek J, Rudolph A, Pfeiffer RF, Foroud T: Alpha-synuclein and familial Parkinson's disease. Mov Disord. 2009, 24 (8): 1125-1131. 10.1002/mds.22524.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Reider CR, Halter CA, Castelluccio PF, Oakes D, Nichols WC, Foroud T: Reliability of reported age at onset for Parkinson's disease. Mov Disord. 2003, 18 (3): 275-279. 10.1002/mds.10391.

    Article  PubMed  Google Scholar 

  29. Gunderson KL, Steemers FJ, Ren H, Ng P, Zhou L, Tsan C, Chang W, Bullis D, Musmacker J, King C, Lebruska LL, Barker D, Oliphant A, Kuhn KM, Shen R: Whole-genome genotyping. Methods in enzymology. 2006, 410: 359-376. 10.1016/S0076-6879(06)10017-8.

    Article  CAS  PubMed  Google Scholar 

  30. Sanna S, Jackson AU, Nagaraja R, Willer CJ, Chen WM, Bonnycastle LL, Shen H, Timpson N, Lettre G, Usala G, Chines PS, Stringham HM, Scott LJ, Dei M, Lai S, Albai G, Crisponi L, Naitza S, Doheny KF, Pugh EW, Ben-Shlomo Y, Ebrahim S, Lawlor DA, Bergman RN, Watanabe RM, Uda M, Tuomilehto J, Coresh J, Hirschhorn JN, Shuldiner AR, Schlessinger D, Collins FS, Davey Smith G, Boerwinkle E, Cao A, Boehnke M, Abecasis GR, Mohlke KL: Common variants in the GDF5-UQCC region are associated with variation in human height. Nat Genet. 2008, 40 (2): 198-203. 10.1038/ng.74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, Clarke R, Heath SC, Timpson NJ, Najjar SS, Stringham HM, Strait J, Duren WL, Maschio A, Busonero F, Mulas A, Albai G, Swift AJ, Morken MA, Narisu N, Bennett D, Parish S, Shen H, Galan P, Meneton P, Hercberg S, Zelenika D, Chen WM, Li Y, Scott LJ, Scheet PA, Sundvall J, Watanabe RM, Nagaraja R, Ebrahim S, Lawlor DA, Ben-Shlomo Y, Davey-Smith G, Shuldiner AR, Collins R, Bergman RN, Uda M, Tuomilehto J, Cao A, Collins FS, Lakatta E, Lathrop GM, Boehnke M, Schlessinger D, Mohlke KL, Abecasis GR: Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet. 2008, 40 (2): 161-169. 10.1038/ng.76.

    Article  CAS  PubMed  Google Scholar 

  32. Li Y, Abecasis G: Mach 1.0: Rapid Haplotype Reconstruction and Missing Genotype Inference. American journal of human genetics. 2006, S79: 2290-

    Google Scholar 

  33. The International HapMap Project. Nature. 2003, 426 (6968): 789-796. 10.1038/nature02168.

  34. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007, 447 (7145): 661-678. 10.1038/nature05911.

  35. Marchini J, Howie B, Myers S, McVean G, Donnelly P: A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet. 2007, 39 (7): 906-913. 10.1038/ng2088.

    Article  CAS  PubMed  Google Scholar 

  36. Lettre G, Lange C, Hirschhorn JN: Genetic model testing and statistical power in population-based association studies of quantitative traits. Genetic epidemiology. 2007, 31 (4): 358-362. 10.1002/gepi.20217.

    Article  PubMed  Google Scholar 

  37. de Bakker PI, Ferreira MA, Jia X, Neale BM, Raychaudhuri S, Voight BF: Practical aspects of imputation-driven meta-analysis of genome-wide association studies. Hum Mol Genet. 2008, 17 (R2): R122-128. 10.1093/hmg/ddn288.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Biasiotto G, Goldwurm S, Finazzi D, Tunesi S, Zecchinelli A, Sironi F, Pezzoli G, Arosio P: HFE gene mutations in a population of Italian Parkinson's disease patients. Parkinsonism Relat Disord. 2008, 14 (5): 426-430. 10.1016/j.parkreldis.2007.10.011.

    Article  PubMed  Google Scholar 

  39. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics. 2007, 81 (3): 559-575. 10.1086/519795.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Risch N, Merikangas K: The future of genetic studies of complex human diseases. Science. 1996, 273 (5281): 1516-1517. 10.1126/science.273.5281.1516.

    Article  CAS  PubMed  Google Scholar 

  41. Pe'er I, Yelensky R, Altshuler D, Daly MJ: Estimation of the multiple testing burden for genomewide association studies of nearly all common variants. Genetic epidemiology. 2008, 32 (4): 381-385. 10.1002/gepi.20303.

    Article  PubMed  Google Scholar 

  42. Dudbridge F, Gusnanto A: Estimation of significance thresholds for genomewide association scans. Genetic epidemiology. 2008, 32 (3): 227-234. 10.1002/gepi.20297.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Lee ST, Nicholls RD, Bundey S, Laxova R, Musarella M, Spritz RA: Mutations of the P gene in oculocutaneous albinism, ocular albinism, and Prader-Willi syndrome plus albinism. N Engl J Med. 1994, 330 (8): 529-534. 10.1056/NEJM199402243300803.

    Article  CAS  PubMed  Google Scholar 

  44. Rebbeck TR, Kanetsky PA, Walker AH, Holmes R, Halpern AC, Schuchter LM, Elder DE, Guerry D: P gene as an inherited biomarker of human eye color. Cancer Epidemiol Biomarkers Prev. 2002, 11 (8): 782-784.

    CAS  PubMed  Google Scholar 

  45. Sulem P, Gudbjartsson DF, Stacey SN, Helgason A, Rafnar T, Magnusson KP, Manolescu A, Karason A, Palsson A, Thorleifsson G, Jakobsdottir M, Steinberg S, Palsson S, Jonasson F, Sigurgeirsson B, Thorisdottir K, Ragnarsson R, Benediktsdottir KR, Aben KK, Kiemeney LA, Olafsson JH, Gulcher J, Kong A, Thorsteinsdottir U, Stefansson K: Genetic determinants of hair, eye and skin pigmentation in Europeans. Nat Genet. 2007, 39 (12): 1443-1452. 10.1038/ng.2007.13.

    Article  CAS  PubMed  Google Scholar 

  46. Zhu G, Evans DM, Duffy DL, Montgomery GW, Medland SE, Gillespie NA, Ewen KR, Jewell M, Liew YW, Hayward NK, Sturm RA, Trent JM, Martin NG: A genome scan for eye color in 502 twin families: most variation is due to a QTL on chromosome 15q. Twin Res. 2004, 7 (2): 197-210. 10.1375/136905204323016186.

    Article  PubMed  Google Scholar 

  47. Duffy DL, Montgomery GW, Chen W, Zhao ZZ, Le L, James MR, Hayward NK, Martin NG, Sturm RA: A three-single-nucleotide polymorphism haplotype in intron 1 of OCA2 explains most human eye-color variation. American journal of human genetics. 2007, 80 (2): 241-252. 10.1086/510885.

    Article  CAS  PubMed  Google Scholar 

  48. Frudakis T, Terravainen T, Thomas M: Multilocus OCA2 genotypes specify human iris colors. Human genetics. 2007, 122 (3-4): 311-326. 10.1007/s00439-007-0401-8.

    Article  PubMed  Google Scholar 

  49. Jannot AS, Meziani R, Bertrand G, Gerard B, Descamps V, Archimbaud A, Picard C, Ollivaud L, Basset-Seguin N, Kerob D, Lanternier G, Lebbe C, Saiag P, Crickx B, Clerget-Darpoux F, Grandchamp B, Soufir N, Melan C: Allele variations in the OCA2 gene (pink-eyed-dilution locus) are associated with genetic susceptibility to melanoma. Eur J Hum Genet. 2005, 13 (8): 913-920. 10.1038/sj.ejhg.5201415.

    Article  CAS  PubMed  Google Scholar 

  50. Driver JA, Logroscino G, Buring JE, Gaziano JM, Kurth T: A prospective cohort study of cancer incidence following the diagnosis of Parkinson's disease. Cancer Epidemiol Biomarkers Prev. 2007, 16 (6): 1260-1265. 10.1158/1055-9965.EPI-07-0038.

    Article  PubMed  Google Scholar 

  51. Olsen JH, Friis S, Frederiksen K: Malignant melanoma and other types of cancer preceding Parkinson disease. Epidemiology. 2006, 17 (5): 582-587. 10.1097/01.ede.0000229445.90471.5e.

    Article  PubMed  Google Scholar 

  52. Olsen JH, Friis S, Frederiksen K, McLaughlin JK, Mellemkjaer L, Moller H: Atypical cancer pattern in patients with Parkinson's disease. Br J Cancer. 2005, 92 (1): 201-205. 10.1038/sj.bjc.6602279.

    Article  CAS  PubMed  Google Scholar 

  53. Herrero Hernandez E: Pigmentation genes link Parkinson's disease to melanoma, opening a window on both etiologies. Med Hypotheses. 2009, 72 (3): 280-284. 10.1016/j.mehy.2008.10.011.

    Article  CAS  PubMed  Google Scholar 

  54. Sturm RA, Teasdale RD, Box NF: Human pigmentation genes: identification, structure and consequences of polymorphic variation. Gene. 2001, 277 (1-2): 49-62. 10.1016/S0378-1119(01)00694-1.

    Article  CAS  PubMed  Google Scholar 

  55. Fedorow H, Tribl F, Halliday G, Gerlach M, Riederer P, Double KL: Neuromelanin in human dopamine neurons: comparison with peripheral melanins and relevance to Parkinson's disease. Prog Neurobiol. 2005, 75 (2): 109-124. 10.1016/j.pneurobio.2005.02.001.

    Article  CAS  PubMed  Google Scholar 

  56. Smith AP, Hoek K, Becker D: Whole-genome expression profiling of the melanoma progression pathway reveals marked molecular differences between nevi/melanoma in situ and advanced-stage melanomas. Cancer Biol Ther. 2005, 4 (9): 1018-1029.

    Article  CAS  PubMed  Google Scholar 

  57. Marciniak SJ, Ron D: Endoplasmic reticulum stress signaling in disease. Physiol Rev. 2006, 86 (4): 1133-1149. 10.1152/physrev.00015.2006.

    Article  CAS  PubMed  Google Scholar 

  58. Lindholm D, Wootz H, Korhonen L: ER stress and neurodegenerative diseases. Cell Death Differ. 2006, 13 (3): 385-392. 10.1038/sj.cdd.4401778.

    Article  CAS  PubMed  Google Scholar 

  59. Rao RV, Bredesen DE: Misfolded proteins, endoplasmic reticulum stress and neurodegeneration. Curr Opin Cell Biol. 2004, 16 (6): 653-662. 10.1016/j.ceb.2004.09.012.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Ryu EJ, Harding HP, Angelastro JM, Vitolo OV, Ron D, Greene LA: Endoplasmic reticulum stress and the unfolded protein response in cellular models of Parkinson's disease. J Neurosci. 2002, 22 (24): 10690-10698.

    CAS  PubMed  Google Scholar 

  61. Cooper AA, Gitler AD, Cashikar A, Haynes CM, Hill KJ, Bhullar B, Liu K, Xu K, Strathearn KE, Liu F, Cao S, Caldwell KA, Caldwell GA, Marsischky G, Kolodner RD, Labaer J, Rochet JC, Bonini NM, Lindquist S: Alpha-synuclein blocks ER-Golgi traffic and Rab1 rescues neuron loss in Parkinson's models. Science. 2006, 313 (5785): 324-328. 10.1126/science.1129462.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Hoozemans JJ, van Haastert ES, Eikelenboom P, de Vos RA, Rozemuller JM, Scheper W: Activation of the unfolded protein response in Parkinson's disease. Biochem Biophys Res Commun. 2007, 354 (3): 707-711. 10.1016/j.bbrc.2007.01.043.

    Article  CAS  PubMed  Google Scholar 

  63. Gasser T: Genetics of Parkinson's disease. Annals of neurology. 1998, 44 (3 Suppl 1): S53-57.

    Article  CAS  PubMed  Google Scholar 

  64. DeStefano AL, Lew MF, Golbe LI, Mark MH, Lazzarini AM, Guttman M, Montgomery E, Waters CH, Singer C, Watts RL, Currie LJ, Wooten GF, Maher NE, Wilk JB, Sullivan KM, Slater KM, Saint-Hilaire MH, Feldman RG, Suchowersky O, Lafontaine AL, Labelle N, Growdon JH, Vieregge P, Pramstaller PP, Klein C, Hubble JP, Reider CR, Stacy M, MacDonald ME, Gusella JF, Myers RH: PARK3 influences age at onset in Parkinson disease: a genome scan in the GenePD study. American journal of human genetics. 2002, 70 (5): 1089-1095. 10.1086/339814.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Pankratz N, Uniacke SK, Halter CA, Rudolph A, Shults CW, Conneally PM, Foroud T, Nichols WC: Genes influencing Parkinson disease onset: replication of PARK3 and identification of novel loci. Neurology. 2004, 62 (9): 1616-1618.

    Article  CAS  PubMed  Google Scholar 

  66. Zhang Y, James M, Middleton FA, Davis RL: Transcriptional analysis of multiple brain regions in Parkinson's disease supports the involvement of specific protein processing, energy metabolism, and signaling pathways, and suggests novel disease mechanisms. Am J Med Genet B Neuropsychiatr Genet. 2005, 137B (1): 5-16. 10.1002/ajmg.b.30195.

    Article  PubMed  Google Scholar 

  67. Henderson DM, Conner SD: A novel AAK1 splice variant functions at multiple steps of the endocytic pathway. Mol Biol Cell. 2007, 18 (7): 2698-2706. 10.1091/mbc.E06-09-0831.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Kametaka S, Moriyama K, Burgos PV, Eisenberg E, Greene LE, Mattera R, Bonifacino JS: Canonical interaction of cyclin G associated kinase with adaptor protein 1 regulates lysosomal enzyme sorting. Mol Biol Cell. 2007, 18 (8): 2991-3001. 10.1091/mbc.E06-12-1162.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Sevlever D, Jiang P, Yen SH: Cathepsin D is the main lysosomal enzyme involved in the degradation of alpha-synuclein and generation of its carboxy-terminally truncated species. Biochemistry. 2008, 47 (36): 9678-9687. 10.1021/bi800699v.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Li YJ, Scott WK, Hedges DJ, Zhang F, Gaskell PC, Nance MA, Watts RL, Hubble JP, Koller WC, Pahwa R, Stern MB, Hiner BC, Jankovic J, Allen FA, Goetz CG, Mastaglia F, Stajich JM, Gibson RA, Middleton LT, Saunders AM, Scott BL, Small GW, Nicodemus KK, Reed AD, Schmechel DE, Welsh-Bohmer KA, Conneally PM, Roses AD, Gilbert JR, Vance JM, Haines JL, Pericak-Vance MA: Age at onset in two common neurodegenerative diseases is genetically controlled. American journal of human genetics. 2002, 70 (4): 985-993. 10.1086/339815.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Bouchard L, Bouchard C, Chagnon YC, Perusse L: Evidence of linkage and association with body fatness and abdominal fat on chromosome 15q26. Obesity (Silver Spring). 2007, 15 (8): 2061-2070. 10.1038/oby.2007.245.

    Article  CAS  Google Scholar 

  72. Verma R, Holmans P, Knowles JA, Grover D, Evgrafov OV, Crowe RR, Scheftner WA, Weissman MM, DePaulo JR, Potash JB, Levinson DF: Linkage disequilibrium mapping of a chromosome 15q25-26 major depression linkage region and sequencing of NTRK3. Biol Psychiatry. 2008, 63 (12): 1185-1189. 10.1016/j.biopsych.2008.02.005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Djousse L, Knowlton B, Hayden MR, Almqvist EW, Brinkman RR, Ross CA, Margolis RL, Rosenblatt A, Durr A, Dode C, Morrison PJ, Novelletto A, Frontali M, Trent RJ, McCusker E, Gomez-Tortosa E, Mayo Cabrero D, Jones R, Zanko A, Nance M, Abramson RK, Suchowersky O, Paulsen JS, Harrison MB, Yang Q, Cupples LA, Mysore J, Gusella JF, MacDonald ME, Myers RH: Evidence for a modifier of onset age in Huntington disease linked to the HD gene in 4p16. Neurogenetics. 2004, 5 (2): 109-114. 10.1007/s10048-004-0175-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Wexler NS, Lorimer J, Porter J, Gomez F, Moskowitz C, Shackell E, Marder K, Penchaszadeh G, Roberts SA, Gayan J, Brocklebank D, Cherny SS, Cardon LR, Gray J, Dlouhy SR, Wiktorski S, Hodes ME, Conneally PM, Penney JB, Gusella J, Cha JH, Irizarry M, Rosas D, Hersch S, Hollingsworth Z, MacDonald M, Young AB, Andresen JM, Housman DE, De Young MM, Bonilla E, Stillings T, Negrette A, Snodgrass SR, Martinez-Jaurrieta MD, Ramos-Arroyo MA, Bickham J, Ramos JS, Marshall F, Shoulson I, Rey GJ, Feigin A, Arnheim N, Acevedo-Cruz A, Acosta L, Alvir J, Fischbeck K, Thompson LM, Young A, Dure L, O'Brien CJ, Paulsen J, Brickman A, Krch D, Peery S, Hogarth P, Higgins DS, Landwehrmeyer B: Venezuelan kindreds reveal that genetic and environmental factors modulate Huntington's disease age of onset. Proceedings of the National Academy of Sciences of the United States of America. 2004, 101 (10): 3498-3503. 10.1073/pnas.0308679101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Andresen JM, Gayan J, Cherny SS, Brocklebank D, Alkorta-Aranburu G, Addis EA, Cardon LR, Housman DE, Wexler NS: Replication of twelve association studies for Huntington's disease residual age of onset in large Venezuelan kindreds. J Med Genet. 2007, 44 (1): 44-50. 10.1136/jmg.2006.045153.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

This project was supported by R01 NS37167, R01 NS036711, the Robert P. & Judith N. Goldberg Foundation, the Bumpus Foundation and the Harvard NeuroDiscovery Center. This study used samples from the NINDS Human Genetics Resource Center DNA and Cell Line Repository http://ccr.coriell.org/ninds, as well as clinical data.

DNA samples contributed by the Parkinson Institute - Istituti Clinici di Perfezionamento, Milan, Italy were from the "Human genetic bank of patients affected by PD and parkinsonisms", supported by Italian Telethon grant n. GTB07001 and by the "Fondazione Grigioni per il Morbo di Parkinson".

A portion of this research was conducted using the Linux Cluster for Genetic Analysis (LinGA-II) funded by the Robert Dawson Evans Endowment of the Department of Medicine at Boston University School of Medicine and Boston Medical Center.

Genotyping services were provided by the Center for Inherited Disease Research (CIDR). CIDR is fully funded through a federal contract from the National Institutes of Health to The Johns Hopkins University, contract number HHSN268200782096C. We particularly thank Drs. Kimberly Doheny and Elizabeth Pugh from CIDR and Mr. Justin Paschall from NCBI for their assistance.

The data used in this study is available at http://www.ncbi.nlm.nih.gov/sites/entrez?db=gap through dbGaP accession number: phs000126.v1.p1.

Funding support for the Mayo-Perlegen LEAPS (Linked Efforts to Accelerate Parkinson's Solutions) Collaboration Study was provided by the Michael J. Fox Foundation. Additional support was provided by the NIH grants ES10751 and NS33978. The datasets used for the analyses described in this manuscript were obtained from dbGaP found at http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000048.v1.p1 through dbGaP study accession numbers phs000048. Data for the Mayo-Perlegen LEAPS Collaboration Study were provided by Demetrius M. Maraganore, MD on behalf of his collaborators at the Mayo Clinic and at Perlegen Sciences, Inc.

The following are members of the PROGENI Steering Committee. University of Tennessee Health Science Center: R. F. Pfeiffer; University of Rochester: F. Marshall, D. Oakes, A. Rudolph, A. Shinaman; Columbia University Medical Center: K. Marder; Indiana University School of Medicine: P.M. Conneally, T. Foroud, C. Halter; University of Kansas Medical Center: K. Lyons; Eli Lilly & Company: E. Siemers; Medical College of Ohio: L. Elmers; University of California, Irvine: N. Hermanowicz.

The following are members of the GenePD Steering Committee. University of Virginia Health System: G.F. Wooten; UMDNJ-Robert Wood Johnson Medical School: L. Golbe; Center for Human Genetic Research, Massachusetts General Hospital, Harvard Medical School: J.F. Gusella; Boston University School of Medicine: R.H. Myers.

We thank the subjects for their participation in this research study.

PSG-PROGENI Investigators and Coordinators: Albany Medical College: S Factor, D Higgins, S Evans; Barrow Neurological Institute: H Shill, M Stacy, J Danielson, L Marlor, K Williamson; Baylor College of Medicine: J Jankovic, C Hunter; Beth Israel Deaconess Medical Center: D Simon, P Ryan, L Scollins; Beth Israel Medical Center: R Saunders-Pullman, K Boyar, C Costan-Toth, E Ohmann; Brigham & Women's Hospital: L Sudarsky, C Joubert; Brown University (Memorial Hospital of RI): J Friedman, K Chou, H Fernandez, M Lannon; Cleveland Clinic Florida-Weston: N Galvez-Jimenez, A Podichetty, K Thompson; Clinical Neuroscience Center: P Lewitt, M DeAngelis; Colorado Neurological Institute: C O'Brien, L Seeberger, C Dingmann, D Judd; Columbia University Medical Center: K Marder, J Fraser, J Harris; Creighton University: J Bertoni, C Peterson; Evanston Northwestern Healthcare: M Rezak, G Medalle; Hotel-Dieu Hospital-Chum: S Chouinard, M Panisset, J Hall, H Poiffaut; Hunter Homes McGuire Veterans Medical Center: V Calabrese, P Roberge; Indiana University School of Medicine: J Wojcieszek, J Belden; Institute For Neurodegenerative Disorders: D Jennings, K Marek, S Mendick; Johns Hopkins University: S Reich, B Dunlop; London Health Sciences Centre: M Jog, C Horn; Mayo Clinic Jacksonville: R Uitti, M Turk; McFarland Neurosciences: T Ajax, J Mannetter; Medical College of Georgia: K Sethi, J Carpenter, B Dill, L Hatch, K Ligon, S Narayan; Medical College of Wisconsin: K Blindauer, K Abou-Samra, J Petit; Medical University of Ohio: L Elmer, E Aiken, K Davis, C Schell, S Wilson; Mount Sinai School of Medicine: M Velickovic, W Koller (deceased), S Phipps; North Shore-LIJ Health System: A Feigin, M Gordon, J Hamann, E Licari, M Marotta-Kollarus, B Shannon, R Winnick; Northwestern University: T Simuni, A Videnovic, A Kaczmarek, K Williams, M Wolff; Ochsner Clinic Foundation: J Rao, M Cook; Ohio State University: M Fernandez, S Kostyk, J Hubble, A Campbell, C Reider, A Seward; Oregon Health & Science University: R Camicioli, J Carter, J Nutt, P Andrews, S Morehouse, C Stone; Ottawa Hospital Civic Site: T Mendis, D Grimes, C Alcorn-Costa, P Gray, K Haas, J Vendette; Pacific Neuroscience Medical Group: J Sutton, B Hutchinson, J Young; Saskatoon Dist Health Board Royal Univ Hosp: A Rajput, A Rajput, L Klassen, T Shirley; Scott & White Hospital/Texas A&M University: B Manyam, P Simpson, J Whetteckey, B Wulbrecht; The Parkinson's & Movement Disorder Institute: D Truong, M Pathak, K Frei, N Luong, T Tra, A Tran, J Vo; Toronto Western Hospital, University Health: A Lang, G Kleiner-Fisman, A Nieves, L Johnston, J So; UMDNJ-School of Osteopathic Medicine: G Podskalny, L Giffin; University of Alabama at Birmingham: P Atchison, C Allen; University of Alberta: W Martin, M Wieler; University of Calgary: O Suchowersky, M Klimek; University of California Irvine: N Hermanowicz, S Niswonger; University of California San Diego: C Shults (deceased), D Fontaine; University of California San Francisco: M Aminoff, C Christine, M Diminno, J Hevezi; University of Chicago: A Dalvi, U Kang, J Richman, S Uy, J Young; University of Cincinnati: A Dalvi, A Sahay, M Gartner, D Schwieterman; University of Colorado Health Sciences Center: D Hall, M Leehey, S Culver, T Derian; University of Connecticut: T Demarcaida, S Thurlow; University of Iowa: R Rodnitzky, J Dobson; University of Kansas Medical Center: K Lyons, R Pahwa, T Gales, S Thomas; University of Maryland School of Medicine: L Shulman, S Reich, W Weiner, K Dustin; University of Miami: K Lyons, C Singer, W Koller (deceased), W Weiner, L Zelaya; University of Minnesota: P Tuite, V Hagen, S Rolandelli, R Schacherer, J Kosowicz; University of New Mexico: P Gordon, J Werner; University of Puerto Rico School of Medicine: C Serrano, S Roque; University of Rochester: R Kurlan, D Berry, I Gardiner; University of South Florida: R Hauser, J Sanchez-Ramos, T Zesiewicz, H Delgado, K Price, P Rodriguez, S Wolfrath; University of Tennessee Health Science Center: R Pfeiffer, L Davis, B Pfeiffer; University of Texas Southwestern Medical Center: R Dewey, B Hayward, A Johnson, M Meacham, B Estes; Wake Forest University School of Medicine: F Walker, V Hunt, C O'Neill; Washington University: B Racette, L Good, M Rundle

PROGENI Molecular Genetic Laboratory: Division of Human Genetics, Cincinnati Children's Hospital Medical Center: Michael W. Pauciulo, Diane K. Marek, Veronika E. Elsaesser

GenePD Investigators and Coordinators: University Southern California School of Medicine: M. Lew; University of Calgary: O. Suchowersky, S. Furtado; University of Lübeck, Germany: C. Klein; UMDNJ-Robert Wood Johnson Medical School: L. Golbe, M.H. Mark; Massachusetts General Hospital, Harvard Medical School: J. Growdon, N. Huggins; University of Virginia Health System: G.F. Wooten; University of Alabama at Birmingham: R. Watts; University of Toronto: M. Guttman; Washington University School of Medicine: B. Racette, J. Perlmutter; Barrow Neurological Institute: L. Marlor, Sun Health Research Institute: H. Shill; University of Miami: C. Singer; Boston University School of Medicine: M.H. Saint-Hilaire, T. Massood; Cleveland Clinic Foundation: K. Baker, I. Itin, A. Ahmed; University of Louisville School of Medicine: I. Litvan; University of Sydney ANZAC Research Institute, Concord Hospital, Sydney, Australia: G. Nicholson, A. Corbett; Struthers Parkinson's Center, Minneapolis: M. Nance; Port City Neurology, Scarborough, ME: E. Drasby; Parkinson's Disease and Movement Disorder Center of Boca Raton: S. Isaacson; Newcastle University, Newcastle upon Tyne, UK: D. Burn, P. Chinnery; General Regional Hospital Bolzano, Bolzano, Italy: P. Pramstaller; University of Arkansas for Medical Sciences: J. Al-hinti; Aarhus University Hospital, Aarhus, Denmark: A. Moller, K. Ostergaard; University of Arizona: S. Sherman; Auckland City Hospital, Auckland, New Zealand: R. Roxburgh, B. Snow; University of Kentucky College of Medicine: J. Slevin, F. Cambi.

GenePD Molecular Genetics Laboratories: Center for Human Genetic Research, Massachusetts General Hospital, Harvard Medical School: M.E. McDonald, M. Sun, L. Mysore, M.A, Anderson, D. Lucente; Neurogenetics Laboratory, Boston University School of Medicine: S. Williamson, M.W. Nagle, B. Brandler.

Author information

Author notes

    Authors and Affiliations

    1. Boston University School of Medicine, Boston, MA, USA

      Jeanne C Latourelle, Alexandra Dumitriu, Jemma B Wilk, Anita L DeStefano & Richard H Myers

    2. Indiana University School of Medicine, Indianapolis, IN, USA

      Nathan Pankratz, Cheryl Halter & Tatiana Foroud

    3. Parkinson Institute, Istituti Clinici di Perfezionamento, Milano, Italy

      Stefano Goldwurm, Gianni Pezzoli & Claudio B Mariani

    4. Boston University School of Public Health, Boston, MA, USA

      Anita L DeStefano

    5. Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA

      James F Gusella

    6. Cincinnati Children's Hospital Medical Center and University of Cincinnati College of Medicine, Cincinnati, Ohio, USA

      William C Nichols

    Authors
    1. Jeanne C Latourelle
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Nathan Pankratz
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Alexandra Dumitriu
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Jemma B Wilk
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Stefano Goldwurm
      View author publications

      You can also search for this author in PubMed Google Scholar

    6. Gianni Pezzoli
      View author publications

      You can also search for this author in PubMed Google Scholar

    7. Claudio B Mariani
      View author publications

      You can also search for this author in PubMed Google Scholar

    8. Anita L DeStefano
      View author publications

      You can also search for this author in PubMed Google Scholar

    9. Cheryl Halter
      View author publications

      You can also search for this author in PubMed Google Scholar

    10. James F Gusella
      View author publications

      You can also search for this author in PubMed Google Scholar

    11. William C Nichols
      View author publications

      You can also search for this author in PubMed Google Scholar

    12. Richard H Myers
      View author publications

      You can also search for this author in PubMed Google Scholar

    13. Tatiana Foroud
      View author publications

      You can also search for this author in PubMed Google Scholar

    Consortia

    the PSG -PROGENI and GenePD Investigators, Coordinators and Molecular Genetic Laboratories

    Corresponding authors

    Correspondence to Richard H Myers or Tatiana Foroud.

    Additional information

    Competing interests

    The authors declare that they have no competing interests.

    Authors' contributions

    JCL participated in conception and design of the study, conducted statistical analyses, participated in the interpretation of data, and drafted the manuscript. NP, AD and JBW participated in the conception and design of the study, conducted statistical analyses, participated in the interpretation of data, and revised the article critically for important intellectual content. SG, GP, CBM, ALD, CH, JFG, WCN, RHM, and TF participated in the conception and design of the study, participated in the interpretation of data, and revised article critically for important intellectual content. All authors read and approved the final manuscript.

    Jeanne C Latourelle, Nathan Pankratz contributed equally to this work.

    Electronic supplementary material

    12881_2009_512_MOESM1_ESM.XLS

    Additional file 1: Supplementary Tables. The excel file includes three supplementary tables. Supplementary Table S1 includes regions showing associations at p < 0.0001 for the additive model GWAS meta-analysis and study-specific results. Supplementary Table S2 includes regions showing associations at p < 0.0001 for the dominant model GWAS meta-analysis and study-specific results. Supplementary Table S3 includes Regions showing associations at p < 0.0001 for the recessive model GWAS meta-analysis and study-specific results (XLS 74 KB)

    Authors’ original submitted files for images

    Below are the links to the authors’ original submitted files for images.

    Authors’ original file for figure 1

    Authors’ original file for figure 2

    Rights and permissions

    Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

    Reprints and Permissions

    About this article

    Cite this article

    Latourelle, J.C., Pankratz, N., Dumitriu, A. et al. Genomewide association study for onset age in Parkinson disease. BMC Med Genet 10, 98 (2009). https://doi.org/10.1186/1471-2350-10-98

    Download citation

    • Received:

    • Accepted:

    • Published:

    • DOI: https://doi.org/10.1186/1471-2350-10-98

    Keywords

    • Genetic Modifier
    • Genome Wide Association Study
    • Parkinson Disease
    • Unfold Protein Response
    • Recessive Model

    BMC Medical Genetics

    ISSN: 1471-2350

    Contact us