Dataset Information

Identification of Candidate Parkinson Disease Genes by Integrating Genome-Wide Association Study, Expression, and Epigenetic Data Sets.

ABSTRACT:

Importance

Substantial genome-wide association study (GWAS) work in Parkinson disease (PD) has led to the discovery of an increasing number of loci shown reliably to be associated with increased risk of disease. Improved understanding of the underlying genes and mechanisms at these loci will be key to understanding the pathogenesis of PD.

Objective

To investigate what genes and genomic processes underlie the risk of sporadic PD.

Design and setting

This genetic association study used the bioinformatic tools Coloc and transcriptome-wide association study (TWAS) to integrate PD case-control GWAS data published in 2017 with expression data (from Braineac, the Genotype-Tissue Expression [GTEx], and CommonMind) and methylation data (derived from UK Parkinson brain samples) to uncover putative gene expression and splicing mechanisms associated with PD GWAS signals. Candidate genes were further characterized using cell-type specificity, weighted gene coexpression networks, and weighted protein-protein interaction networks.

Main outcomes and measures

It was hypothesized a priori that some genes underlying PD loci would alter PD risk through changes to expression, splicing, or methylation. Candidate genes are presented whose change in expression, splicing, or methylation are associated with risk of PD as well as the functional pathways and cell types in which these genes have an important role.

Results

Gene-level analysis of expression revealed 5 genes (WDR6 [OMIM 606031], CD38 [OMIM 107270], GPNMB [OMIM 604368], RAB29 [OMIM 603949], and TMEM163 [OMIM 618978]) that replicated using both Coloc and TWAS analyses in both the GTEx and Braineac expression data sets. A further 6 genes (ZRANB3 [OMIM 615655], PCGF3 [OMIM 617543], NEK1 [OMIM 604588], NUPL2 [NCBI 11097], GALC [OMIM 606890], and CTSB [OMIM 116810]) showed evidence of disease-associated splicing effects. Cell-type specificity analysis revealed that gene expression was overall more prevalent in glial cell types compared with neurons. The weighted gene coexpression performed on the GTEx data set showed that NUPL2 is a key gene in 3 modules implicated in catabolic processes associated with protein ubiquitination and in the ubiquitin-dependent protein catabolic process in the nucleus accumbens, caudate, and putamen. TMEM163 and ZRANB3 were both important in modules in the frontal cortex and caudate, respectively, indicating regulation of signaling and cell communication. Protein interactor analysis and simulations using random networks demonstrated that the candidate genes interact significantly more with known mendelian PD and parkinsonism proteins than would be expected by chance.

Conclusions and relevance

Together, these results suggest that several candidate genes and pathways are associated with the findings observed in PD GWAS studies.

SUBMITTER: Kia DA

PROVIDER: S-EPMC7851759 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Genome wide association studies (GWAS) are a population-scale approach to the identification of segments of the genome in which genetic variations may contribute to disease risk. Current methods focus on the discovery of single nucleotide polymorphisms (SNPs) associated with disease traits. As there are many SNPs within identified risk loci, and the majority of these are situated within non-coding regions, a key challenge is to identify and prioritize variants affecting regulatory sequences that are likely to contribute to the phenotype assessed.We focused investigation on SNPs within lung and breast cancer GWAS loci that reached genome-wide significance for potential roles in gene regulation with a specific focus on SNPs likely to disrupt transcription factor binding sites. Within risk loci, the regulatory potential of sub-regions was classified using relevant open chromatin and epigenetic high throughput sequencing data sets from the ENCODE project in available cancer and normal cell lines. Furthermore, transcription factor affinity altering variants were predicted by comparison of position weight matrix scores between disease and reference alleles. Lastly, ChIP-seq data of transcription associated factors and topological domains were included as binding evidence and potential gene target inference.The sets of SNPs, including both the disease-associated markers and those in high linkage disequilibrium with them, were significantly over-represented in regulatory sequences of cancer and/or normal cells; however, over-representation was generally not restricted to disease-relevant tissue specific regions. The calculated regulatory potential, allelic binding affinity scores and ChIP-seq binding evidence were the three criteria used to prioritize candidates. Fitting all three criteria, we highlighted breast cancer susceptibility SNPs and a borderline lung cancer relevant SNP located in cancer-specific enhancers overlapping multiple distinct transcription associated factor ChIP-seq binding sites.Incorporating high throughput sequencing epigenetic and transcription factor data sets from both cancer and normal cells into cancer genetic studies reveals potential functional SNPs and informs subsequent characterization efforts.

Project description:BackgroundOsteoarthritis (OA) is a common skeletal system disease that has been partially attributed to genetic factors. The hand is frequently affected, which seriously affects the patient's quality of life. However, the pathogenetic mechanism of hand osteoarthritis (hand OA) is still elusive.MethodsA genome-wide association study (GWAS) summary of hand OA was obtained from the UK Biobank dataset, which contains data from a total of 452,264 White British individuals, including 37,782 OA patients. The transcriptome-wide association study (TWAS) of hand OA was performed using FUnctional Summary-based ImputatiON (FUSION) with the skeletal muscle and blood as gene expression references. The significant genes identified by TWAS were further subjected to gene set enrichment analysis (GSEA) with the Database for Annotation, Visualization and Integrated Discovery (DAVID) tool. Furthermore, we compared the genes and gene sets identified by our TWAS with that of a knee OA mRNA expression profile to detect the genes and gene sets shared by TWAS and mRNA expression profiles in OA. The mRNA expression profiles of 18 normal knee cartilages and 20 OA knee cartilages were acquired from the Gene Expression Omnibus database (accession number: GSE114007).ResultsTWAS identified 177 genes with P < 0.05 for the skeletal muscle, including ANKRD44 (P = 0.0001), RIC3 (P = 0.0003), and AC005154.6 (P = 0.0004). TWAS identified 423 genes with P < 0.05 for the blood, including CRIM1 (P = 0.0002), ZNF880 (P = 0.0002), and NCKIPSD (P = 0.0003). After comparing the results of the TWAS to those of the mRNA expression profiling of OA, we identified 5 common genes, including DHRS3 (log2fold = - 1.85, P = 3.31 × 10- 9) and SKP2 (log2fold = 1.36, P = 1.62 × 10- 8). GSEA of TWAS identified 51 gene ontology (GO) terms for hand OA, for example, protein binding (P = 0.0003) and cytosol (P = 0.0020). We also detected 6 common GO terms shared by TWAS and mRNA expression profiling, including protein binding (PTWAS = 2.54 × 10- 4, PmRNA = 3.42 × 10- 8), extracellular exosome (PTWAS = 0.02, PmRNA = 1.18 × 10- 4), and cytoplasm (PTWAS = 0.0183, PmRNA = 0.0048).ConclusionIn this study, we identified 5 candidate genes and 6 GO terms related to hand OA, which may help to uncover the pathogenesis of hand OA. It should be noted that the possible difference in the gene expression profiles between hand OA and knee OA may affect our study results, which should be interpreted with caution.