Abstract
Genome-wide association study (GWAS) has seen great strides in revealing initial insights into the genetic architecture of Parkinson’s disease (PD). Since GWAS signals often reside in non-coding regions, relatively few of the associations have implicated specific biological mechanisms. Here, we aimed to integrate the GWAS results with large-scale expression quantitative trait loci (eQTL) in 13 brain tissues to identify candidate causal genes for PD. We conducted a transcriptome-wide association study (TWAS) for PD using the summary statistics of over 480,000 individuals from the most recent PD GWAS. We identified 18 genes significantly associated with PD after Bonferroni corrections. The most significant gene, LRRC37A2, was associated with PD in all 13 brain tissues, such as in the hypothalamus (P = 6.12 × 10−22) and nucleus accumbens basal ganglia (P = 5.62 × 10−21). We also identified eight conditionally independent genes, including four new genes at known PD loci: CD38, LRRC37A2, RNF40, and ZSWIM7. Through conditional analyses, we demonstrated that several of the GWAS significant signals on PD could be driven by genetically regulated gene expression. The most significant TWAS gene LRRC37A2 accounts for 0.855 of the GWAS signal at its loci, and ZSWIM7 accounts for all the GWAS signals at its loci. We further identified several phenotypes previously associated with PD by querying the single nucleotide polymorphisms (SNPs) in the final model of the identified genes in phenome databases. In conclusion, we prioritized genes that are likely to affect PD by using a TWAS approach and identified phenotypes associated with PD.
Similar content being viewed by others
Introduction
Parkinson’s disease (PD) is the second most common age-related neurodegenerative disorder after Alzheimer’s disease, characterized by the loss of nigrostriatal neurons in the substantia nigra1. PD is more common in the elderly, and the prevalence increased from 1% in people over 60 years to 3–4% in those over 80 years2. Based on twin and family studies, the heritability of PD has been estimated to be at least 27% and up to 60%3,4, suggesting a substantial involvement of genetic factors drives the phenotypic variance. Genome-wide association study (GWAS) has seen great strides and invaluable utilities in revealing initial insights into PD’s genetic architecture. The largest-to-date GWAS for PD identified 90 independent genome-wide significant association signals, including 38 novel risk variants5. Together, associated loci represent only a small fraction of PD’s genetic etiology, leaving a substantial proportion of genetic risk factors uncharacterized6.
Most disease- and trait-associated variants mapped through GWAS lie within non-coding regions and are not in linkage disequilibrium with any nonsynonymous coding single nucleotide polymorphisms (SNPs)7,8. Thus, the causal variants and target susceptible genes of most GWAS risk loci have not been identified9. Some studies, including the most recently published PD GWAS, have linked the GWAS risk loci to the nearest gene, which will inevitably lead to biases against longer genes and may not accurately depict the locus’s real effect. Many of these genetic variants are located in the gene regulatory elements10,11 and have hypothesized to influence complex traits by modulating gene regulatory programs, an intermediate between genetic variation and complex disease. Expression quantitative trait locus (eQTL) analysis seeks to identify genetic variants that affect the gene expression; several studies have successfully used this approach to identify putative susceptibility genes at GWAS risk loci for PD12. The enrichment of eQTLs of trait-associated variants also showed the importance of gene expression regulation13,14.
Recent studies have reported that regulatory variants may account for a large proportion of disease heritability that has not yet been identified through GWAS15. Many of these variants may have modest effect sizes, and it is difficult to identify them via a typical individual SNP-based GWAS study, even with a considerable sample size. Transcriptome-wide association studies (TWAS) that systematically investigate the relationship between genetically predicted gene expression and disease risk, providing a powerful approach to identify disease risk genes and uncovering possible causal genes at loci identified previously by GWAS16,17,18,19,20. Considering the essential role of epigenetic features in predicting gene expression, we developed an epigenetic element-based TWAS21. Briefly, for a given gene, we used eQTL data to impute the total expression across a large cohort of genotyped individuals with the epigenetic features as prior, followed by a test of association with disease risk. Compared with the single SNP-based GWAS study, TWAS can increase power in identifying disease-related genes, either by reducing the burden of multiple comparisons or by aggregating multiple expression-altering variants into a single test.
In the current study, we carried out a TWAS to prioritize candidate PD genes and to better understand the primary mechanisms that underlie PD genetic risk factors using the largest PD cohort currently available with 15,056 PD cases, 18,618 proxies, and 449,056 controls. We prioritized 18 genes whose predicted expression is significantly associated with PD risk, comprising 67 associations, after Bonferroni corrections. Through conditional analyses, we demonstrated that several of the GWAS significant signals on PD could be driven by genetically regulated gene expression. We identified eight conditionally independent genes through conditional analyses, including four new genes at known PD loci. We identified several phenotypes previously associated with PD by querying the SNPs in the final models of identified genes in phenome databases. Overall, our study suggests that a TWAS approach considering both genetic and epigenetic effects on gene expression is a powerful method to identify specific genes and mechanisms at each GWAS locus as determinants of PD risk.
Results
Transcriptome-wide association study
To identify genes associated with PD, we conducted a TWAS using summary data comprising 15,056 PD cases, 18,618 UK Biobank proxies, and 449,056 controls (Fig. 1). We used genotyping and transcriptome data from 13 brain tissues to build epigenetic-based gene expression prediction models. The tissue abbreviations and sample sizes are listed in Supplementary Table 1. We identified a total of 18 genes significantly associated with PD after Bonferroni corrections (P < 2.55 × 10−6, Fig. 2 and Supplementary Table 2), comprising 67 associations. The most significant gene, LRRC37A2, was associated with PD in all 13 tissues tested, such as the hypothalamus (P = 6.12 × 10−22) and the nucleus accumbens basal ganglia (P = 5.62 × 10−21). Amongst the signals, we found six genes only significant in one tissue, including RNF40 (P = 4.42 × 10−10) and CPLX1 (P = 6.36 × 10−7) in the cerebellar hemisphere, VKORC1 in the cortex (P = 2.66 × 10−9), MAP3K14 (P = 3.52 × 10−8) and GAK (P = 1.85 × 10−7) in the cerebellum, as well as CENPV in the frontal cortex BA9 (P = 6.07 × 10−7). Most of the genes (14/18) remained significant after permutation, suggesting that their signals were genuine and not due to chance.
GWAS, genome-wide association study; TFBS, transcription factor binding site; DHS, DNase I hypersensitive sites; HMM, chromatin state segmentation by hidden Markov model; TWAS, transcriptome-wide association study; PD, Parkinson’s disease; PheWAS, phenome-wide association study.
The blue line represents the Bonferroni-corrected significant thresholds, P = 2.55 × 10−6. Conditionally independent genes are listed in blue letters.
Conditionally testing GWAS signals
Since most of the TWAS-identified genes overlapped with GWAS PD loci, we performed the conditional and joint analyses to check whether these signals were due to multiple associated features and how much GWAS signal remains after the gene expression is removed. We identified eight conditionally independent genes, including LRRC37A2, MMRN1, CD38, RNF40, GPNMB, ZSWIM7, GAK, and CPLX1. We observed that ZSWIM7 accounts for all the signals at its loci (rs4566208 lead SNPGWAS P = 3.90 × 10−8, conditioned on ZSWIM7 lead SNPGWAS P = 1) (Fig. 3a). Similarly, conditioning on LRRC37A2 accounts for most of the signal at its loci (rs199452 lead SNPGWAS P = 4.80 × 10−21, conditioned on LRRC37A2 lead SNPGWAS P = 0.17, accounting for 0.855 of the variances) (Fig. 3b). Conditioning on RNF40 accounts for most of the loci’s variance (rs8050588 lead SNPGWAS P = 1.60 × 10−10, conditioned on RNF40 lead SNPGWAS P = 0.066, accounting for 0.713 of the variances) (Fig. 3c). We also found that conditioned on the expression GPNMB accounts for 0.8 of the variances (rs466225 lead SNPGWAS P = 4.30 × 10−9, conditioned on GPNMB lead SNPGWAS P = 0.241) (Fig. 3d).
a Chromosome 17 regional association plot (part 1). b Chromosome 17 regional association plot (part 2). c Chromosome 16 regional association plot. d Chromosome 7 regional association plot. The marginally associated TWAS genes are shown in blue, and the conditionally significant genes are shown in green. The bottom panel shows a regional Manhattan plot of the GWAS data before (gray) and after (blue) conditioning on the green genes’ predicted expression.
Replication analyses
To assess our results’ robustness, we performed replication analyses to validate the conditionally independent genes. We found 46 reference models for these eight genes and five genes (LRRC37A2, MMRN1, GPNMB, ZSWIM7, CPLX1) significantly associated with PD in at least one brain tissue (P = 1.09 × 10−3, Supplementary Table 3). Similarly, LRRC37A2 was associated with PD in all 13 tissues in the replication data, such as the hippocampus (P = 5.86 × 10−11) and the hypothalamus (P = 1.41 × 10−10).
Partitioned heritability of TWAS-identified genes
We partitioned the heritability explained by SNPs around TWAS-identified genes and found that the genes identified by ten out of 13 tissues significantly contribute to PD’s heritability comparing to the percentage of SNPs (Table 1). For example, the genes identified in the brain cerebellum panel explained 5.52% (1.35%) of the estimated heritability, a 127.7× enrichment compared to the percentage of SNPs (P = 2.50 × 10−5). Partitioned heritability of TWAS genes in all the 13 tissues explained 22.17% of the estimated heritability, a 32.3× enrichment compared to the percentage of SNPs (P = 1.68 × 10−5).
Gene set enrichment analyses
The Genetic Association Database (GAD) disease enrichment analyses detected six diseases (Fig. 4a), including PD (P = 0.011), tobacco use disorder (P = 0.012), cholesterol (P = 0.023), LDL cholesterol (P = 0.025), bone mineral density (P = 0.038), and cleft lip and cleft palate (P = 0.044). Three KEGG pathways were significantly enriched (Fig. 4b), including mTOR signaling pathway (P = 9.90 × 10−4), selenocompound metabolism (P = 0.031), as well as PPAR signaling pathway (P = 0.038). We also identified several GO terms enriched in TWAS-identified PD genes (Fig. 4c–e), such as regulation of bone resorption (biological processes, P = 2.87 × 10−4), immunological synapse (cellular component, P = 1.97 × 10−4), and syntaxin-1 binding (molecular function, P = 5.75 × 10−3).
a GAD disease enrichment analyses of TWAS genes. b KEGG pathway enrichment analyses of TWAS genes. c–e Top ten GO terms enriched in TWAS genes, including biological processes (c), cellular component (d), and molecular function (e). f GWAS Catalog enrichment analyses of TWAS genes. The histogram shows the expected number of genes with P < 0.01 based on 10,000 random permutations. The large red point shows the observed number of previously known PD genes that fall below this threshold.
GWAS Catalog enrichment analyses
We found that nine TWAS-identified genes had been reported in the GWAS Catalog. The null expectation based on 10,000 randomly drawn gene sets proved that PD had a significant enrichment (P = 1.0 × 10−4) of GWAS Catalog-reported genes (Fig. 4f) in the TWAS results, which suggested that there are likely to be genuine disease associations among the genes that fail to meet strict genome-wide significance.
Phenome-wide association studies
We identified 166 phenotypes associated with the SNPs in the final model of the TWAS-identified conditionally independent genes, including activities such as neurological, psychiatric, and cognitive (Supplementary Fig. 1). We further conducted the genetic correlation between PD and 122 identified traits with currently available GWAS data to determine their relationship. The latest GWAS summary statistics are listed in Supplementary Table 4. We found that the age at first sexual intercourse positively correlated with PD (Fig. 5). We also reported five phenotypes negatively correlated with PD, including three impedance measures, heel bone mineral density, and current tobacco smoking (Fig. 5). Importantly, tobacco use disorder and bone mineral density have been identified through GAD disease enrichment analyses.
*P < 0.05/122, **P < 0.01/122, error bars indicate standard error of the genetic correlations.
Discussion
PD is a long-term movement disorder that affects approximately seven million people globally. Although recent GWAS has seen great strides in identifying risk loci associated with PD, the functional significance of these associations remains elusive. We conducted a PD TWAS using the summary statistics of over 480,000 individuals from the most recent PD GWAS. This approach creates genotype-expression reference panels using public consortia through Lasso and Elastic Net with epigenetic annotations as prior, allowing for imputation and association testing of independent large-scale data. We identified 18 genes associated with PD after Bonferroni corrections, comprising 67 associations localized to seven different regions in the genome.
Conditional and joint analyses identified eight independent genes and demonstrated that the TWAS expression signals were driving the significance for several previously implicated PD loci. For example, the most significant TWAS gene LRRC37A2 accounts for 0.855 of the GWAS signal at its loci, and ZSWIM7 accounts for all the GWAS signals at its loci. These results imply a limited residual association signal from the genetic variant in the GWAS locus after considering these predicted expression signals. Our identifications provide further support for three genes previously implicated by GWAS, whose expression was significantly associated with a possible causal change in PD risk by summary-based Mendelian randomization, including MMRN1, GPNMB, and GAK5. Another gene, CPLX1, was possibly associated with at least one QTL in public reference datasets but did not pass the Bonferroni corrections.
Additionally, our TWAS implicates four new genes at known PD loci, including CD38, LRRC37A2, RNF40, and ZSWIM7. For instance, CD38 encodes the cluster of differentiation 38, also known as cyclic ADP ribose hydrolase, a glycoprotein on many immune cells’ surfaces. CD38 strongly expressed in brain cells, including neurons, astrocytes, and microglial cells22. Of note, several data tend to indicate that CD38 expression increase in the brain as a consequence of aging23,24, the primary risk factor associated with the vast majority of neurodegenerative diseases, including PD. Moreover, several experimental data demonstrated that CD38 knockout mice are protected from neurodegenerative and neuroinflammatory insults25,26. Future studies could interrogate whether expression differences of other candidate genes are consistent with our findings.
In parallel to ours, Kia and colleagues27 conducted TWAS using expression weights from the CommonMind Consortium (CMC) dorsolateral prefrontal cortex to identify genes associated with PD. While considering the essential role of epigenetic features in predicting gene expression, our approach integrated the epigenetic information in the original TWAS to find gene–trait associations. Among the 18 genes we identified, eight genes (MMRN1, CD38, NUPL2, GPNMB, RNF40, VKORC1, ZSWIM7, and CENPV) were significantly heritable in CMC dorsolateral prefrontal cortex and qualified for the TWAS analyses. According to the published results, seven of the above eight genes associated with PD risk at an FDR level of 0.05 (Supplementary Table 5), including two genes (CD38, and GPNMB) with solid evidence for colocalization (PPH4 > 0.75)27. At the same time, numerous previously reported PD risk loci did not implicate in our TWAS. For example, Li and colleagues28 found the predicted gene expression of SNCA was associated with PD in peripheral monocytes, and published evidence suggested the association between MAPT expression and PD29,30. The most significant gene we identified, LRRC37A2, is located at the end of the MAPT locus, and the second most significant finding is MMRN1 at the SNCA locus. However, MAPT and SNCA were not included in this project since they were not significantly heritable in any brain tissues at current sample sizes. Considering the previous evidence and the fact that we did not perform these associations in our project, we cannot determine the driver genes of PD at these loci, which can be regarded as one of the limitations of TWAS. Further studies with a larger sample size of reference data are needed to identify the driver genes of PB at these loci.
The GAD disease enrichment analyses of TWAS-identified genes detected six GAD diseases, including PD itself. Other enriched phenotypes, such as tobacco use disorder, cholesterol, and bone mineral density, also have been reported to be associated with PD. For example, tobacco use disorder, also known as nicotine dependence, is a chronic, relapsing disease defined as a compulsive craving to use it, despite harmful social consequences31. Epidemiological studies show that smoking is associated with a lower incidence of PD32. Moreover, nicotine stimulates striatal dopamine neurons that are damaged in PD and protect against neuronal insults in experimental models. Besides, previous researches have suggested that higher total- and LDL-cholesterol levels may be associated with lower risk and beneficial outcomes in PD33,34. Bone mineral density is the most widely used predictor for osteoporosis, and increasing evidence suggests that neurological conditions, including PD, are associated with an excess rate of osteoporosis and fracture risk35.
KEGG pathway enrichment analyses identified three pathways, including the mTOR signaling pathway, PPAR signaling pathway, and selenocompound metabolism. mTOR is a serine/threonine kinase that is the central component of mTORC1 and mTORC2 multiprotein complexes. mTOR regulates many integrated physiological functions of the nerve system, including neuronal development, synaptic plasticity, memory storage, and cognition36. The deregulation of mTOR signaling appears to be a common hallmark of human neurological disorders, including PD37, and mTORC1-induced transcripts enriched in a cluster of genes related to PD38. We also identified multiple GO terms that had been reported to be associated with PD. For example, above we have discussed the correlation between PD and osteoporosis, and three GO terms associated with osteoporosis were detected, such as regulation of bone resorption, regulation of bone remodeling, and bone resorption.
In summary, by using the TWAS method that we generated recently, we identified 18 genes associated with PD after Bonferroni corrections, comprising 67 associations localized to seven different regions in the genome. We identified eight conditionally independent genes, and we demonstrated that several of the GWAS significant signals on PD could be driven by genetically regulated gene expression. Our TWAS implicates four new genes at known PD loci: CD38, LRRC37A2, RNF40, and ZSWIM7. We further identified several phenotypes associated with PD by querying the SNPs in the final model of identified genes in phenome databases. In conclusion, we prioritized genes that are likely to affect PD by using a TWAS approach and identified phenotypes associated with PD.
Methods
GWAS summary statistics
Discovery GWAS data
We used the most recent GWAS summary statistics for PD. Details on participant ascertainment and quality control were previously reported by Nalls et al.5. The summary statistics comprising 15,056 PD cases, 18,618 UK Biobank proxy-cases (individuals who do not have PD but have a first-degree relative that does), and 449,056 controls. Since our approach depends on having dense summary-level data to overlap with the expression weights closely, we did not prune SNPs in the summary data.
Replicate GWAS data
Another GWAS summary statistic of PD, comprising 4238 PD cases and 4239 controls, was used as the replication data of this study. Details about ascertainment and quality control were previously reported by Pankratz et al.39. Similarly, all datasets employed standard UK Brain Bank criteria40 for the diagnosis of PD, with a modification to allow the inclusion of cases that had a family history of PD, since familial PD cases may have a more substantial genetic contribution than sporadic PD, making them potentially more informative for genetic studies.
Transcriptome-wide association study
We used transcriptome and high-density genotyping data of European decedent from the Genotype-Tissue Expression (GTEx) study Pilot Project V8 (dbGap accession: phs000424.v8.p2) to establish gene expression prediction models41. In this project, we used genotyping and transcriptome data from 13 brain tissues to build epigenetic-based gene expression prediction models.
As described in our previously reported method21, we performed a TWAS using reference panels derived from tissue-specific gene expression coupled with genotypic data with the epigenetic features. For each gene x, we trained and evaluated the models for gene expression prediction in each round y of ten-fold cross-validation using the following steps. (a) We performed eQTL analyses with SNPs located within 1 Mb of the transcription start/end sites of the gene using the training data; (b) We then annotated the SNPs with epigenetic annotations. For each SNP, an epigenomic feature was labeled if the SNP overlapped with the feature. (c) We obtained multiple SNP sets according to the eQTL P-value threshold and epigenetic annotation. (d) For each SNP set z, we built an expression prediction model in the training dataset by using the Lasso and the Elastic Net (α = 0.5) methods as implemented in the glmnet R package. For each model, we evaluated its prediction performance by the coefficient of determination R2 between the predicted gene expression and the observed gene expression of the testing data, and averaged all the cross-validation data. For each gene x, the model with the highest mean R2 in the testing data was selected as the best model. Based on the parameters of the best model, we performed the eQTL analyses again using all the samples in the reference data and constructed each gene’s final prediction model. We estimated the associations between predicted expressions and PD with the combination of SNP-trait effect sizes while accounting for linkage disequilibrium among SNPs via Functional Summary-based Imputation (FUSION, http://gusevlab.org/projects/fusion/).
Since highly heritable genes were significantly enriched in trait associations18, we only focused on genes whose heritability did not overlap with zero with 95% confidence interval. We used a strict Bonferroni-corrected study-wise threshold with P = 2.55 × 10−6 (0.05/19,632, the total number of highly heritable genes across tissues). We applied the 1000 Genomes v3 LD panel for the TWAS. To assess the possibility of inflated association statistics from TWAS, we performed a permutation test. We shuffled the eQTL weights (n = 10,000) and recomputed an empirical association statistics conditional on the GWAS effects at the locus via FUSION.
Conditional analyses
We performed conditional and joint analyses at each genome-wide Bonferroni-corrected TWAS genes to determine how much GWAS signal remains after the expression association from TWAS is removed. Moreover, for regions where TWAS identified multiple associated features, we jointly modeled these genes to determine the independent signals. Each PD GWAS association was conditioned on the joint gene model, one SNP at a time. We set the overlapping genes in the range of 1 Mb around each SNP, and the defined regions included only the transcribed region of the genes. We used the FUSION tool to perform the conditional and joint analyses with the cis-genetic component of expression we generated.
Partitioned heritability estimation
The partitioned analysis is to quantify the heritability directly explained by SNPs in each functional category using the summary statistics compared with the null expectation, equal to the percentage of SNPs in the gene set. We estimated the partitioned heritability of PD by SNPs around TWAS-identified genes to see whether the identified genes significantly contribute to the PD heritability. Partitioned heritability of PD was estimated using LD score regression (LDSC) following the previously described methodology42. We partitioned the heritability explained by TWAS-identified genes with a less stringent threshold (P < 0.01) in each tissue by SNPs within 2 kb of the genes. We generated the LD score files using the open-source software available at https://github.com/bulik/ldsc/wiki/Partitioned-Heritability.
Gene set enrichment analyses
To demonstrate TWAS’s ability to identify PD-related genes, we performed the GAD disease enrichment analyses via the Database for Annotation, Visualization, and Integrated Discovery (DAVID) tool (https://david.ncifcrf.gov/home.jsp) using default settings. A relaxed threshold of 0.01, rather than Bonferroni-correction, was used for GAD disease enrichment analyses since Bonferroni-correction assumes independence while genes tend to correlate due to co-expression. More genes will allow for better recapitulation and prioritization of appropriate pathways. We performed the GO and pathway enrichment analyses using the clusterProfiler R package. The GO terms (including biological processes, cellular components, and molecular functions) and pathways from KEGG (https://www.genome.jp/kegg/) were analyzed in this study.
GWAS Catalog enrichment analyses
We tested whether the TWAS-identified genes enriched known PD-associated genes. The reported PD genes derived from the NHGRI GWAS Catalog43 identified using GWAS were regarded as the set of known PD-associated genes. We excluded studies that included the discovery dataset to make sure our known gene list was independent of the current analysis. We then counted the number of known disease-associated genes that had a TWAS P-value below 0.01. We compared this count to the null expectation based on 10,000 randomly drawn gene sets of similar size to the known disease gene set to derive an enrichment P-value.
Phenome-wide association studies
To identify phenotypes that may be associated with PD, we conducted a phenome-wide association study (pheWAS) for each SNP in the final model of the identified genes. We reported the first five phenotypes (excluding PD) through public data provided by GWASAtlas (https://atlas.ctglab.nl) according to the P-values. We further conducted the genetic correlation between PD and identified traits with currently available GWAS data to determine their relationship. The analyses were carried out using LDSC at https://github.com/bulik/ldsc. The latest GWAS summary statistics were used for the correlation.
Standard protocol approvals and participant consents
This study protocol was approved by the Ethics Committee of Hunan Brain Hospital. The study based on GWAS summary statistics does not require informed consent from all study participants. The methods were carried out in accordance with the approved guidelines.
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
The discovery GWAS summary statistic of PD was obtained from the link (https://bit.ly/2ofzGrk) shared by Nalls and colleagues5. The replicate GWAS summary statistic was obtained from the Genome-Wide Repository of Associations Between SNPs and Phenotypes database (https://grasp.nhlbi.nih.gov/FullResults.aspx). Other data that support the findings of this study are available from the corresponding author upon reasonable request.
Code availability
The source code is available at https://github.com/studentyaoshi/etwas.
References
Dauer, W. & Przedborski, S. Parkinson’s disease: mechanisms and models. Neuron 39, 889–909 (2003).
Ascherio, A. & Schwarzschild, M. A. The epidemiology of Parkinson’s disease: risk factors and prevention. Lancet Neurol. 15, 1257–1272 (2016).
Do, C. B. et al. Web-based genome-wide association study identifies two novel loci and a substantial genetic component for Parkinson’s disease. PLoS Genet. 7, e1002141 (2011).
Hamza, T. H. & Payami, H. The heritability of risk and age at onset of Parkinson’s disease after accounting for known genetic risk factors. J. Hum. Genet. 55, 241–243 (2010).
Nalls, M. A. et al. Identification of novel risk loci, causal insights, and heritable risk for Parkinson’s disease: a meta-analysis of genome-wide association studies. Lancet Neurol. 18, 1091–1102 (2019).
Konig, I. R. Challenges in disentangling the genetic background of Parkinson’s disease. Lancet Neurol. 18, 1069–1070 (2019).
Maurano, M. T. et al. Systematic localization of common disease-associated variation in regulatory DNA. Science 337, 1190–1195 (2012).
Hindorff, L. A. et al. Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc. Natl Acad. Sci. USA 106, 9362–9367 (2009).
Spisak, S. et al. CAUSEL: an epigenome- and genome-editing pipeline for establishing function of noncoding GWAS variants. Nat. Med. 21, 1357–1363 (2015).
Dunham, I. et al. An integrated encyclopedia of DNA elements in the human genome. Nature 489, 57–74 (2012).
Kundaje, A. et al. Integrative analysis of 111 reference human epigenomes. Nature 518, 317–330 (2015).
Chang, D. et al. A meta-analysis of genome-wide association studies identifies 17 new Parkinson’s disease risk loci. Nat. Genet. 49, 1511–1516 (2017).
Nicolae, D. L. et al. Trait-associated SNPs are more likely to be eQTLs: annotation to enhance discovery from GWAS. PLoS Genet. 6, e1000888 (2010).
Cookson, W. et al. Mapping complex disease traits with global gene expression. Nat. Rev. Genet. 10, 184–194 (2009).
Gusev, A. et al. Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases. Am. J. Hum. Genet. 95, 535–552 (2014).
Gamazon, E. R. et al. A gene-based association method for mapping traits using reference transcriptome data. Nat. Genet. 47, 1091 (2015).
Gusev, A. et al. Transcriptome-wide association study of schizophrenia and chromatin activity yields mechanistic disease insights. Nat. Genet. 50, 538 (2018).
Gusev, A. et al. Integrative approaches for large-scale transcriptome-wide association studies. Nat. Genet. 48, 245–252 (2016).
Barbeira, A. N. et al. Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics. Nat. Commun. 9, 1825 (2018).
Manor, O. & Segal, E. Robust prediction of expression differences among human individuals using only genotype information. PLoS Genet. 9, e1003396 (2013).
Yao, S. et al. Epigenetic element-based transcriptome-wide association study identifies novel genes for bipolar disorder. Schizophr. Bull. 47, sbab023 (2021).
Guerreiro, S., Privat, A. L., Bressac, L. & Toulorge, D. CD38 in neurodegeneration and neuroinflammation. Cells 9, 471 (2020).
Camacho-Pereira, J. et al. CD38 dictates age-related NAD decline and mitochondrial dysfunction through an SIRT3-dependent mechanism. Cell Metab. 23, 1127–1139 (2016).
Braidy, N. et al. Mapping NAD+ metabolism in the brain of ageing Wistar rats: potential targets for influencing brain senescence. Biogerontology 15, 177–198 (2014).
Ma, Y., Wu, D., Ding, X. & Ying, W. CD38 plays key roles in both antioxidation and cell survival of H2O2-treated primary rodent astrocytes. Int. J. Physiol. Pathophysiol. Pharmacol. 6, 102–108 (2014).
Kou, W. et al. CD38 regulation in activated astrocytes: implications for neuroinflammation and HIV-1 brain infection. J. Neurosci. Res. 87, 2326–2339 (2009).
Kia, D. A. et al. Identification of candidate Parkinson disease genes by integrating genome-wide association study, expression, and epigenetic data sets. JAMA Neurol. 78, 464–472 (2021).
Li, Y., Wong, G., Humphrey, J. & Raj, T. Prioritizing Parkinson’s disease genes using population-scale transcriptomic data. Nat. Commun. 10, 994 (2019).
Tobin, J. E. et al. Haplotypes and gene expression implicate the MAPT region for Parkinson disease: the GenePD Study. Neurology 71, 28–34 (2008).
Rittman, T. et al. Regional expression of the MAPT gene is associated with loss of hubs in brain networks and cognitive impairment in Parkinson disease and progressive supranuclear palsy. Neurobiol. Aging 48, 153–160 (2016).
Falcone, M., Lee, B., Lerman, C. & Blendy, J. A. Translational research on nicotine dependence. Curr. Top. Behav. Neurosci. 28, 121–150 (2016).
Quik, M. Smoking nicotine and Parkinson’s disease. Trends Neurosci. 27, 561–568 (2004).
Huang, X. et al. Lower low-density lipoprotein cholesterol levels are associated with Parkinson’s disease. Mov. Disord. 22, 377–381 (2007).
Scigliano, G. et al. Reduced risk factors for vascular disorders in Parkinson disease patients: a case-control study. Stroke 37, 1184–1188 (2006).
Torsney, K. M. et al. Bone health in Parkinson’s disease: a systematic review and meta-analysis. Mov. Disord. 29, S551 (2014).
Bockaert, J. & Marin, P. mTOR in brain physiology and pathologies. Physiol. Rev. 95, 1157–1187 (2015).
Laplante, M. & Sabatini, D. M. mTOR signaling in growth control and disease. Cell 149, 274–293 (2012).
Duvel, K. et al. Activation of a metabolic gene regulatory network downstream of mTOR complex 1. Mol. Cell 39, 171–183 (2010).
Pankratz, N. et al. Meta-analysis of Parkinson’s disease: identification of a novel locus, RIT2. Ann. Neurol. 71, 370–384 (2012).
Hughes, A. J., Daniel, S. E., Kilford, L. & Lees, A. J. Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: a clinico-pathological study of 100 cases. J. Neurol. Neurosurg. Psychiatry 55, 181–184 (1992).
Consortium, G. T. Human genomics. The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science 348, 648–660 (2015).
Finucane, H. K. et al. Partitioning heritability by functional annotation using genome-wide association summary statistics. Nat. Genet. 47, 1228–1235 (2015).
Buniello, A. et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 47, D1005–D1012 (2019).
Acknowledgements
This study is supported by Science and Technology Major Project of Hunan Provincial Science and Technology Department (2020SK2123, 2018SK1030), National Natural Science Foundation of China (82101601), Natural Science Foundation of Hunan Province (2019JJ80024, 2020JJ8081), Scientific Research Fund of Hunan Provincial Health and Family Planning Commission (C2019043), and Hunan Science and Technology Innovation Project (2018SK50605).
Author information
Authors and Affiliations
Contributions
S.Y. and X.Z. contributed equally and thus share first authorship. Study concept and design: S.Y., X.Z., X.L., L.L., and X.W. Collection, analysis, and interpretation of data: S.Y., X.Z., S.Z., Y.Z., B.L., W.K., and Y.G. Drafting of the manuscript: S.Y. Critical revision of the manuscript: S.Y., X.Z., S.Z., Y.Z., B.L., W.K., Y.G., X.L., L.L., and X.W.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yao, S., Zhang, X., Zou, SC. et al. A transcriptome-wide association study identifies susceptibility genes for Parkinson’s disease. npj Parkinsons Dis. 7, 79 (2021). https://doi.org/10.1038/s41531-021-00221-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41531-021-00221-7
This article is cited by
-
Novel and experimental therapeutics for the management of motor and non-motor Parkinsonian symptoms
Neurological Sciences (2024)
-
Expanding causal genes for Parkinson’s disease via multi-omics analysis
npj Parkinson's Disease (2023)
-
Analysis of subcellular RNA fractions demonstrates significant genetic regulation of gene expression in human brain post-transcriptionally
Scientific Reports (2023)
