Parkinson-associated risk variant in enhancer element produces subtle effect on target gene expression

Genome-wide association studies (GWAS) have identified numerous genetic variants associated with complex diseases but mechanistic insights are impeded by the lack of understanding of how specific risk variants functionally contribute to the underlying pathogenesis(1). It has been proposed that cis-acting effects of non-coding risk variants on gene expression are a major factor for phenotypic variation of complex traits and disease susceptibility. Recent genome-scale chromatin mapping studies have highlighted the enrichment of GWAS variants in regulatory DNA elements of disease-relevant cell types(2–6). Furthermore, single nucleotide polymorphism (SNP)-specific changes in transcription factor (TF) binding are correlated with heritable alterations in chromatin state and considered a major mediator of sequence-dependent regulation of gene expression(7–10). Here we describe a novel strategy to functionally dissect the cis-acting effect of genetic risk variants in regulatory elements on gene expression by combining genome-wide epigenetic information with clustered regularly-interspaced short palindromic repeats (CRISPR)/Cas9 genome editing in human pluripotent stem cells (hPSCs). By generating a genetically precisely controlled experimental system we identify a common Parkinson’s disease (PD)-associated risk variant in a non-coding distal enhancer element that regulates the expression of alpha-synuclein (SNCA), a key gene implicated in the pathogenesis of PD. Our data suggest that the transcriptional deregulation of SNCA is associated with sequence-dependent binding of the brain-specific TFs EMX2 and NKX6-1. This work establishes an experimental paradigm to functionally connect genetic variation with disease relevant phenotypes.