The Parkinson's disease variant rs356182 regulates neuronal differentiation independently from alpha-synuclein.

One of the most significant risk variants for Parkinson's disease (PD), rs356182, is located at the PD-associated locus near the alpha-synuclein (α-syn) encoding gene, SNCA. SNCA-proximal variants, including rs356182, are thought to function in PD risk through enhancers via allele-specific regulatory effects on SNCA expression. However, this interpretation discounts the complex activity of genetic enhancers and possible non-conical functions of α-syn. Here we investigated a novel risk mechanism for rs356182. We use CRISPR-Cas9 in LUHMES cells, a model for dopaminergic midbrain neurons, to generate precise hemizygous lesions at rs356182. The PD-protective (A/-), PD-risk (G/-) and wild-type (A/G) clones were neuronally differentiated and then compared transcriptionally and morphologically. Among the affected genes was SNCA, whose expression was promoted by the PD-protective allele (A) and repressed in its absence. In addition to SNCA, hundreds of genes were differentially expressed and associated with neurogenesis and axonogenesis-an effect not typically ascribed to α-syn. We also found that the transcription factor FOXO3 specifically binds to the rs356182 A-allele in differentiated LUHMES cells. Finally, we compared the results from the rs356182-edited cells to our previously published knockouts of SNCA and found only minimal overlap between the sets of significant differentially expressed genes. Together, the data implicate a risk mechanism for rs356182 in which the risk-allele (G) is associated with abnormal neuron development, independent of SNCA expression. We speculate that these pathological effects manifest as a diminished population of dopaminergic neurons during development leading to the predisposition for PD later in life.

Gene expression asymmetry in Parkinson's Disease; variation of CCT and BEX gene expression levels are correlated with hemisphere specific severity.

Parkinson's Disease (PD) develops unilaterally, which may be related to brain hemispheric differences in gene expression. Here we measured bulk RNA-seq levels in neuronal nuclei obtained from prefrontal cortex postmortem brain samples from males and females with PD and from healthy controls. Left and right hemispheres from each brain were related the side of symptom onset and compared. We employed two a priori approaches; first we identified genes differentially expressed between PD and controls and between left vs right PD brain hemispheres. Second, we examined the presence of, and correlates to, variable asymmetry seen in candidate PD differentially expressed genes. We found large variation among individuals with PD, and PD stratification by gene expression similarity was required for patterns of genetic asymmetry to emerge. For a subset of PD brains, hemispherical variation of CCT and BEX gene levels correlated with the side of PD symptom onset.

Parkinson's disease risk enhancers in microglia.

Genome-wide association studies have identified thousands of single nucleotide polymorphisms that associate with increased risk for Parkinson's disease (PD), but the functions of most of them are unknown. Using assay for transposase-accessible chromatin (ATAC) and H3K27ac chromatin immunoprecipitation (ChIP) sequencing data, we identified 73 regulatory elements in microglia that overlap PD risk SNPs. To determine the target genes of a "risk enhancer" within intron two of SNCA, we used CRISPR-Cas9 to delete the open chromatin region where two PD risk SNPs reside. The loss of the enhancer led to reduced expression of multiple genes including SNCA and the adjacent gene MMRN1. It also led to expression changes of genes involved in glucose metabolism, a process that is known to be altered in PD patients. Our work expands the role of SNCA in PD and provides a connection between PD-associated genetic variants and underlying biology that points to a risk mechanism in microglia.

Post-GWAS knowledge gap: the how, where, and when.

Genetic risk for complex diseases very rarely reflects only Mendelian-inherited phenotypes where single-gene mutations can be followed in families by linkage analysis. More commonly, a large set of low-penetrance, small effect-size variants combine to confer risk; they are normally revealed in genome-wide association studies (GWAS), which compare large population groups. Whereas Mendelian inheritance points toward disease mechanisms arising from the mutated genes, in the case of GWAS signals, the effector proteins and even general risk mechanism are mostly unknown. Instead, the utility of GWAS currently lies primarily in predictive and diagnostic information. Although an amazing body of GWAS-based knowledge now exists, we advocate for more funding towards the exploration of the fundamental biology in post-GWAS studies; this research will bring us closer to causality and risk gene identification. Using Parkinson's Disease as an example, we ask, how, where, and when do risk loci contribute to disease?

Thermodynamic dissection of the polymerizing and editing modes of a DNA polymerase.

DNA polymerases with intrinsic proofreading activity interact with DNA primer/templates in two distinct modes, corresponding to the complexes formed during the 5'-3' polymerization or 3'-5' editing of a nascent DNA chain. Thermodynamic measurements designed to quantify the energetic contributions of individual DNA-protein contacts in either the polymerizing or editing complexes are complicated by the fact that both species exist in solution and are not resolved in conventional DNA-protein binding assays. To overcome this problem, we have developed a new binding analysis that combines information from steady-state and time-resolved fluorescence experiments and uses the Klenow fragment of Escherichia coli DNA polymerase I (KF) and fluorescently labeled primer/template oligonucleotides as a model polymerase-DNA system. Steady-state fluorescence titrations are used to evaluate the overall affinity of KF for the primer/template, while time-resolved fluorescence anisotropy is used to quantify the equilibrium fractions of the primer/template bound in the polymerizing and editing modes. From a combined analysis of both data, the equilibrium constant and hence standard free energy change associated with each binding mode can be obtained unequivocally. This method is initially used to determine the equilibrium constants describing binding of a correctly base-paired primer/template to the 5'-3' polymerase and 3'-5' exonuclease sites of KF. It is then extended to quantify the extent to which these parameters are affected by the introduction of mismatches into the primer/template, and by rearrangement of specific side-chains in the exonuclease domain of the protein. While these perturbants were originally designed to demonstrate the utility of our new approach, they are also relevant in their own right since they have helped identify some hitherto unknown determinants of polymerase fidelity.

Dimerization of the Klenow fragment of Escherichia coli DNA polymerase I is linked to its mode of DNA binding.

Upon associating with a proofreading polymerase, the nascent 3' end of a DNA primer/template has two possible fates. Depending upon its suitability as a substrate for template-directed extension or postsynthetic repair, it will bind either to the 5'-3' polymerase active site, yielding a polymerizing complex, or to the 3'-5' exonuclease site, yielding an editing complex. In this investigation, we use a combination of biochemical and biophysical techniques to probe the stoichiometry, thermodynamic, and kinetic stability of the polymerizing and editing complexes. We use the Klenow fragment of Escherichia coli DNA polymerase I (KF) as a model proofreading polymerase and oligodeoxyribonucleotide primer/templates as model DNA substrates. Polymerizing complexes are produced by mixing KF with correctly base paired (matched) primer/templates, whereas editing complexes are produced by mixing KF with multiply mismatched primer/templates. Electrophoretic mobility shift titrations carried out with matched and multiply mismatched primer/templates give rise to markedly different electrophoretic patterns. In the case of the matched primer/template, the KF.DNA complex is represented by a slow moving band. However, in the case of the multiply mismatched primer/template, the complex is predominantly represented by a fast moving band. Analytical ultracentrifugation measurements indicate that the fast and slow moving bands correspond to 1:1 and 2:1 KF.DNA complexes, respectively. Fluorescence anisotropy titrations reveal that KF binds with a higher degree of cooperativity to the matched primer/template. Taken together, these results indicate that KF is able to dimerize on a DNA primer/template and that dimerization is favored when the first molecule is bound in the polymerizing mode, but disfavored when it is bound in the editing mode. We suggest that self-association of the polymerase may play an important and as yet unexplored role in coordinating high-fidelity DNA replication.

Determinants of DNA mismatch recognition within the polymerase domain of the Klenow fragment.

The Klenow fragment of Escherichia coli DNA polymerase I catalyzes template-directed synthesis of DNA and uses a separate 3'-5' exonuclease activity to edit misincorporated bases. The polymerase and exonuclease activities are contained in separate structural domains. In this study, nine Klenow fragment derivatives containing mutations within the polymerase domain were examined for their interaction with model primer-template duplexes. The partitioning of the DNA primer terminus between the polymerase and 3'-5' exonuclease active sites of the mutant proteins was assessed by time-resolved fluorescence anisotropy, utilizing a dansyl fluorophore attached to the DNA. Mutation of N845 or R668 disrupted favorable interactions between the Klenow fragment and a duplex containing a matched terminal base pair but had little effect when the terminus was mismatched. Thus, N845 and R668 are required for recognition of correct terminal base pairs in the DNA substrate. Mutation of N675, R835, R836, or R841 resulted in tighter polymerase site binding of DNA, suggesting that the side chains of these residues induce strain in the DNA and/or protein backbone. A double mutant (N675A/R841A) showed an even greater polymerase site partitioning than was displayed by either single mutation, indicating that such strain is additive. In both groups of mutant proteins, the ability to discriminate between duplexes containing matched or mismatched base pairs was impaired. In contrast, mutation of K758 or Q849 had no effect on partitioning relative to wild type, regardless of DNA mismatch character. These results demonstrate that DNA mismatch recognition is dependent on specific amino acid residues within the polymerase domain and is not governed solely by thermodynamic differences between correct and mismatched base pairs. Moreover, this study suggests a mechanism whereby the Klenow fragment is able to recognize polymerase errors following a misincorporation event, leading to their eventual removal by the 3'-5' exonuclease activity.