# A Genome Wide Association Study Links Glutamate Receptor Pathway to Sporadic Creutzfeldt-Jakob Disease Risk

Pascual Sanchez-Juan, Matthew T. Bishop, Gabor G Kovacs, Miguel Calero, Yurii S. Aulchenko, Anna Ladogana, Alison Boyd, Victoria Lewis, Claudia Ponto, Olga Calero, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Pascual Sanchez-Juan, Matthew T. Bishop, Gabor G Kovacs, Miguel Calero, Yurii S. Aulchenko, et al. A Genome Wide Association Study Links Glutamate Receptor Pathway to Sporadic Creutzfeldt-Jakob Disease Risk. PLoS ONE, 2015, 10 (4), pp.e0123654. 10.1371/journal.pone.0123654. hal-01224275

HAL Id: hal-01224275

https://hal.sorbonne-universite.fr/hal-01224275

Submitted on 4 Nov 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

RESEARCH ARTICLE

# A Genome Wide Association Study Links Glutamate Receptor Pathway to Sporadic Creutzfeldt-Jakob Disease Risk

Pascual Sanchez-Juan¹, Matthew T. Bishop², Gabor G. Kovacs³, Miguel Calero⁴,⁵, Yurii S. Aulchenko⁶,⁷,⁸, Anna Ladogana⁹, Alison Boyd¹⁰, Victoria Lewis¹⁰, Claudia Ponto¹¹, Olga Calero⁴, Anna Poleggi⁹, Ángel Carracedo¹²,¹³, Sven J. van der Lee⁶, Thomas Ströbel³, Fernando Rivadeneira⁶,¹⁴, Albert Hofman⁶, Stéphane Haïk¹⁵, Onofre Combarros¹, José Berciano¹, Andre G. Uitterlinden⁶,¹⁴, Steven J. Collins¹⁰, Herbert Budka³, Jean-Philippe Brandel¹⁵, Jean Louis Laplanche¹⁶, Maurizio Pocchiari⁹, Inga Zerr¹¹, Richard S. G. Knight², Robert G. Will², Cornelia M. van Duijn⁶*

1. Neurology Department, University Hospital "Marqués de Valdecilla". Instituto de Investigación "Marqués de Valdecilla" IDIVAL and Centro de Investigación Biomédica en Red sobre Enfermedades Neurodegenerativas (CIBERNED). Santander, Spain.
2. The National Creutzfeldt-Jakob disease Research and Surveillance Unit, University of Edinburgh, United Kingdom.
3. Institute of Neurology, Medical University Vienna, Vienna, Austria.
4. Chronic Disease Programme and CIBERNED. Carlos III Institute of Health. Madrid, Spain.
5. Alzheimer Disease Research Unit, CIEN Foundation, Carlos III Institute of Health, Alzheimer Center Reina Sofia Foundation, Madrid, Spain.
6. Department of Epidemiology, Erasmus Medical Centre, Rotterdam, the Netherlands.
7. Institute of Cytology and Genetics SB RAS, Novosibirsk, Russia.
8. Novosibirsk State University, Novosibirsk, Russia.
9. Department of Cell Biology and Neurosciences Instituto Superiore di Sanità, Roma, Italy.
10. Department of Pathology, The University of Melbourne, Parkville, 3010, Australia.
11. Department of Neurology, Clinical Dementia Centre, University Medical Center and German Center for Neurodegenerative Diseases (DZNE)—site Göttingen, Göttingen, Germany.
12. Fundación Pública Galega de Medicina Xenómica, CIBERER, Grupo de Medicina Xenómica-Universidad de Santiago de Compostela, Santiago de Compostela, Spain.
13. Center of Excellence in Genomic Medicine Research (CEGMR), King Abdulaziz University, Jeddah, KSA.
14. Department of Internal Medicine, Erasmus Medical Center, Rotterdam, the Netherlands.
15. Sorbonne Universités, UPMC Univ Paris 06 UMR S 1127, and Inserm, U 1127, and CNRS UMR 7225, and ICM, F-75013, Paris, France; AP-HP, Hôpital de la Pitié Salpêtrière, Cellule Nationale de Référence des maladies de Creutzfeldt-Jakob, F-75013, Paris, France.
16. Service de biochimie et biologie moleculaire, Laboratoire associé au CNR "ATNC", Hôpital Lariboisière, AP-HP, Paris, France

* c.vanduijn@erasmusmc.nl

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Sanchez-Juan P, Bishop MT, Kovacs GG, Calero M, Aulchenko YS, Ladogana A, et al. (2015) A Genome Wide Association Study Links Glutamate Receptor Pathway to Sporadic Creutzfeldt-Jakob Disease Risk. PLoS ONE 10(4): e0123654.
doi:10.1371/journal.pone.0123654

**Academic Editor:** Gianluigi Zanusso, University of Verona, ITALY

**Received:** January 12, 2015

**Accepted:** February 23, 2015

**Published:** April 28, 2015

**Copyright:** © 2015 Sanchez-Juan et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This study was supported by: UK: the National CJD Research and Surveillance UK Unit is funded by the Department of Health and the Scottish Government Health Department. The National CJD Research and Surveillance Unit is funded by the Policy Research Programme in the Department of Health. Germany: This work was supported by a grant from the European Commission (Protecting the food chain from prions: shaping European priorities through basic and applied research (PRIORITY, N°

## Abstract

We performed a genome-wide association (GWA) study in 434 sporadic Creutzfeldt-Jakob disease (sCJD) patients and 1939 controls from the United Kingdom, Germany and The Netherlands. The findings were replicated in an independent sample of 1109 sCJD and 2264 controls provided by a multinational consortium. From the initial GWA analysis we selected 23 SNPs for further genotyping in 1109 sCJD cases from seven different countries. Five SNPs were significantly associated with sCJD after correction for multiple testing. Subsequently these five SNPs were genotyped in 2264 controls. The pooled analysis, including 1543 sCJD cases and 4203 controls, yielded two genome wide significant results: rs6107516 (p-value=7.62x10⁻⁹) a variant tagging the prion protein gene (PRNP); and rs6951643 (p-value=1.66x10⁻⁸) tagging the Glutamate Receptor Metabotropic 8 gene (GRM8). Next we analysed the data stratifying by country of origin combining samples from
222887) Project number: FP7-KBBE-2007-2A). The study was performed within the recently established Clinical Dementia Center at the University Medical Center Göttingen and was partly supported by grants from the EU Joint Programme – Neurodegenerative Disease Research (JPND - DEMTEST "Biomarker based diagnosis of rapid progressive dementias-optimization of diagnostic protocols", 01ED1201A). This study was funded by the Robert Koch Institute through funds from the Federal Ministry of Health (grant no. 1369-341). Italy: The Italian Registry of CJD and related disorders is funded by the Ministry of Health, National Centre for Disease Prevention and Control, Central Actions. This work was partly supported by grant from the EU Joint Programme – Neurodegenerative Disease Research (JPND - DEMTEST "Biomarker based diagnosis of rapid progressive dementias-optimization of diagnostic protocols", 01ED1201A). The Netherlands: the generation and management of genome-wide association study (GWAS) genotype data for the Rotterdam Study is supported by the Netherlands Organization of Scientific Research NWO Investments (nr. 175.010.2005.011, 911-03-012). This study is funded by the Research Institute for Diseases in the Elderly (014-93-015; RIDE2), the Netherlands Genomics Initiative (NGI)/Netherlands Organisation for Scientific Research (NWO) project nr. 050-060-810. The Rotterdam Study is funded by Erasmus Medical Center and Erasmus University, Rotterdam, Netherlands Organization for the Health Research and Development (ZonMw), the Research Institute for Diseases in the Elderly (RIDE), the Ministry of Education, Culture and Science, the Ministry for Health, Welfare and Sports, the European Commission (DG XII), and the Municipality of Rotterdam. YSA is supported by Russian Science Foundation (RSCF) grant 14-14-00313. Spain. PSJ was supported by a grant from FIS (PI12/02288) and JPND project DEMTEST (PI11/03028). AC is supported by PI13/01136 Acción Estratégica de Salud del Instituto de Salud Carlos III e INNOPHARMA. Australia: The Australian National Creutzfeldt-Jakob Disease Registry (ANCJDR) is funded by the Commonwealth Department of Health. SJC is supported by a NHMRC Practitioner Fellowship (#APP1005816). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

the pooled analysis with genotypes from the *1000 Genomes Project* and imputed genotypes from the *Rotterdam Study* (Total n=12967). The meta-analysis of the results showed that rs6107516 (p-value=3.00x10⁻⁸) and rs6951643 (p-value=3.91x10⁻⁵) remained as the two most significantly associated SNPs. Rs6951643 is located in an intronic region of *GRM8*, a gene that was additionally tagged by a cluster of 12 SNPs within our top100 ranked results. *GRM8* encodes for mGluR8, a protein which belongs to the metabotropic glutamate receptor family, recently shown to be involved in the transduction of cellular signals triggered by the prion protein. Pathway enrichment analyses performed with both Ingenuity Pathway Analysis and ALIGATOR postulates glutamate receptor signalling as one of the main pathways associated with sCJD. In summary, we have detected GRM8 as a novel, non-PRNP, genome-wide significant marker associated with heightened disease risk, providing additional evidence supporting a role of glutamate receptors in sCJD pathogenesis.

## Introduction

Sporadic Creutzfeldt-Jakob disease (sCJD), although rare, with a yearly incidence of one to two cases per million, is the most common form of human prion disease. This group of disorders is characterized by spongiform changes in the brain, as well as accumulation of misfolded, often protease-resistant, conformers ($PrP^{Sc}$) of the normal prion protein ($PrP^C$). The $PrP^C$ gene (*PRNP*) plays a central role in prion disease susceptibility. Expression of $PrP^C$ is indispensable for disease transmission [1] and the polymorphism coding for methionine (M) or valine (V) at codon 129 (*PRNP* M129V) has been linked to disease risk [2]. Homozygosity at *PRNP* M129V has been consistently associated to sCJD, being one of the strongest common genetic risk factors reported for neurodegenerative diseases. The remarkable disease-determining effect of this *PRNP* polymorphism is observed in variant CJD, a subtype acquired from dietary exposure to bovine spongiform encephalopathy [3], where all definite and probable clinical cases studied to date have been *PRNP*129MM [4].

Similar to other diseases, several genetic association studies with candidate genes have been performed on sCJD susceptibility [5]. Only one previous genome wide association study (GWAS) of sCJD risk has been published to date, showing that the *PRNP* locus was strongly associated with disease risk, specifically with rs1799990 (*PRNP* M129V) [6].

To further scrutinize genomic variations related to sCJD risk we have performed a three-stage GWAS encompassing a total of 1,543 sCJD cases and 4,203 controls, as well as a meta-analysis encompassing data from the *1000 Genomes Project* and imputations from the *Rotter-dam Study*.

## Results

Demographic and clinical features of the sCJD case populations are shown in S1 Table. The Q-Q plots for autosomal and X chromosome SNPs are given in S1 Fig; during discovery (stage one) the genomic inflation factor $\lambda$ was 1.053 for autosomal and 1.057 for X chromosome SNPs.

S2 Table shows the top 100 SNPs associated with sCJD at discovery stage, sorted by allelic differences p-value. A total of 23 SNPs were taken forward to replication. In stage two, we successfully genotyped 22 of the 23 SNPs with the Sequenom iPLEx GOLD platform in an independent population of 1,109 samples of sCJD. Table 1 shows that five SNPs remained significant
<table><caption>Table 1. SNPs genotyped in stage two.</caption><thead><tr><th rowspan="2">SNP</th><th rowspan="2">CJD risk allele in Discovery population</th><th rowspan="2">GENE tagged</th><th colspan="3">Stage one: Discovery</th><th colspan="2">Stage two: Replication in independent sCJD Cases</th></tr><tr><th>Controls Frequencies of sCJD risk allele (n = 1939)</th><th>Cases Frequencies of sCJD risk allele (n = 434)</th><th>P1df*</th><th>Cases Frequencies of sCJD risk allele (n = 1109)</th><th>P-value**</th></tr></thead><tbody><tr><td>rs10061929</td><td>A</td><td>FLJ43080</td><td>0.159</td><td>0.224</td><td>4.58E-05</td><td>0.200</td><td>5.89E-05</td></tr><tr><td>rs10915708</td><td>T</td><td>NA</td><td>0.140</td><td>0.207</td><td>1.35E-06</td><td>0.158</td><td>0.0626</td></tr><tr><td>rs11075924</td><td>C</td><td>NA</td><td>0.498</td><td>0.565</td><td>5.52E-05</td><td>0.500</td><td>0.8808</td></tr><tr><td>rs11245373</td><td>T</td><td>NA</td><td>0.037</td><td>0.069</td><td>2.48E-05</td><td>0.049</td><td>0.0238</td></tr><tr><td>rs12102156</td><td>T</td><td>NA</td><td>0.863</td><td>0.912</td><td>2.00E-05</td><td>0.849</td><td>0.1457</td></tr><tr><td>rs12188818</td><td>C</td><td>NA</td><td>0.796</td><td>0.859</td><td>6.44E-05</td><td>0.809</td><td>0.2544</td></tr><tr><td>rs12419710</td><td>A</td><td>NA</td><td>0.210</td><td>0.279</td><td>2.52E-05</td><td>0.219</td><td>0.3746</td></tr><tr><td>rs17060736</td><td>G</td><td>NA</td><td>0.112</td><td>0.160</td><td>7.13E-05</td><td>0.121</td><td>0.3134</td></tr><tr><td>rs17115017</td><td>A</td><td>GRIA1</td><td>0.050</td><td>0.082</td><td>2.23E-04</td><td>0.072</td><td>5.00E-04</td></tr><tr><td>rs17833759</td><td>C</td><td>GRIN2B</td><td>0.085</td><td>0.140</td><td>7.74E-07</td><td>0.095</td><td>0.1955</td></tr><tr><td>rs196940</td><td>C</td><td>ERN1</td><td>0.758</td><td>0.818</td><td>6.19E-05</td><td>0.758</td><td>0.9631</td></tr><tr><td>rs2240344</td><td>G</td><td>NA</td><td>0.432</td><td>0.600</td><td>4.47E-19</td><td>0.443</td><td>0.4194</td></tr><tr><td>rs2627829</td><td>A</td><td>INPP4B</td><td>0.939</td><td>0.970</td><td>1.83E-04</td><td>0.931</td><td>0.2188</td></tr><tr><td>rs392184</td><td>T</td><td>MACROD2</td><td>0.052</td><td>0.096</td><td>2.00E-06</td><td>0.07</td><td>0.003</td></tr><tr><td>rs565559</td><td>T</td><td>NA</td><td>0.553</td><td>0.618</td><td>2.02E-04</td><td>0.579</td><td>0.057</td></tr><tr><td>rs6027482</td><td>C</td><td>LOC100131710</td><td>0.187</td><td>0.247</td><td>6.69E-05</td><td>0.181</td><td>0.5476</td></tr><tr><td>rs6107516</td><td>C</td><td>PRNP</td><td>0.752</td><td>0.827</td><td>6.92E-06</td><td>0.804</td><td>3.27E-06</td></tr><tr><td>rs6463269</td><td>G</td><td>NA</td><td>0.084</td><td>0.134</td><td>3.53E-06</td><td>GENOTYPE FAILED</td><td>-</td></tr><tr><td>rs6496239</td><td>T</td><td>NA</td><td>0.805</td><td>0.862</td><td>4.96E-05</td><td>0.821</td><td>0.1195</td></tr><tr><td>rs6820644</td><td>T</td><td>NA</td><td>0.357</td><td>0.433</td><td>3.40E-06</td><td>0.647</td><td>0.7566</td></tr><tr><td>rs6951643</td><td>A</td><td>GRM8</td><td>0.507</td><td>0.592</td><td>8.09E-06</td><td>0.573</td><td>6.79E-07</td></tr><tr><td>rs9521699</td><td>A</td><td>COL4A2</td><td>0.137</td><td>0.194</td><td>3.84E-05</td><td>0.173</td><td>2.00E-04</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><th rowspan="2">SNP</th><th rowspan="2">CJD risk allele in Discovery population</th><th rowspan="2">GENE tagged</th><th colspan="3">Stage one: Discovery</th><th colspan="2">Stage two: Replication in independent sCJD Cases</th></tr><tr><th>Controls Frequencies of sCJD risk allele (n = 1939)</th><th>Cases Frequencies of sCJD risk allele (n = 434)</th><th>P1df*</th><th>Cases Frequencies of sCJD risk allele (n = 1109)</th><th>P-value**</th></tr></thead><tbody><tr><td>rs9830696</td><td>T</td><td>NA</td><td>0.120</td><td>0.173</td><td>1.02E-05</td><td>0.104</td><td>0.0711</td></tr></tbody></table>

sCJD: Sporadic Creutzfeldt-Jakob Disease

In bold SNPs significant after Bonferroni correction (Replication P-value <0.0023)

* P-value for allelic differences adjusted by Country of origin by PCA

**Chi Squared p-value from comparison between Allele frequencies of independent sCJD cases and stage one controls

NA: SNP in intergenic region.

doi:10.1371/journal.pone.0123654.t001

after Bonferroni correction. In stage three we genotyped those five SNPs in a population of 2,264 independent controls with Sequenom iPLEx GOLD. S2 Fig depicts neatly separated genotype clusters of the five SNPs studied in stage three. Only two SNPs were successfully replicated at stage three and reached genome wide significant p-values after meta-analysis of discovery and replication results: rs6107516 (p-value = 7.62x10⁻⁹) tagging PRNP and in linkage disequilibrium (LD) with PRNP M129V (rs1799990), and rs6951643 (p-value = 1.66x10⁻⁸) an intronic SNP within GRM8 (Glutamate Receptor Metabotropic 8) in chromosome 7. Table 2 shows differences in allelic frequencies in the discovery, replication and pooled population sets including 1,543 sCJD cases and 4,203 controls. The A allele of rs6951643 was associated with a 1.27 fold increased risk of sCJD (95%CI = 1.17–1.38), and this effect was consistently observed across all tested populations (Table 3).

Rs6107516 (p-value = 3.00x10⁻⁸) and rs6951643 (p-value = 3.91x10⁻⁵) remained as the two most statistically significant associations in a meta-analysis of the discovery data adjusted for PCAs and analysis of replication samples stratified by country. This meta-analysis included discovery and replication samples plus data from the *1000 Genomes Project* and imputed genotypes from the *Rotterdam Study* (Total n = 12967), (Table 4). Within the top 100 SNPs, there

<table><caption>**Table 2. SNPs studied in stage three.**</caption><thead><tr><th rowspan="2">ID</th><th rowspan="2">CJD risk allele in Discovery population</th><th rowspan="2">GENE tagged</th><th colspan="4">Discovery Population</th><th colspan="3">Replication Population</th><th colspan="3">Pooled Samples</th></tr><tr><th>Case, Control Frequencies of CJD risk allele (n = 434, n = 1939)</th><th>P-value</th><th>P-value*</th><th>OR (95%CI)</th><th>Case, Control Frequencies of CJD risk allele (n = 1109, n = 2264)</th><th>P-value</th><th>OR (95%CI)</th><th>Case, Control Frequencies of CJD risk allele (n = 1543, n = 4203)</th><th>P-value</th><th>OR (95%CI)</th></tr></thead><tbody><tr><td>rs10061929</td><td>A</td><td>FLJ43080</td><td>0.224, 0.159</td><td>4.30E-06</td><td>4.58E-05</td><td>1.53 (1.27–1.83)</td><td>0.200, 0.196</td><td>0.6882</td><td>1.03 (0.90–1.17)</td><td>0.207, 0.179</td><td>0.0007</td><td>1.20 (1.08–1.33)</td></tr><tr><td>rs17115017</td><td>A</td><td>GRIA1</td><td>0.082, 0.050</td><td>3.00E-04</td><td>2.23E-04</td><td>1.69 (1.27–2.24)</td><td>0.072, 0.077</td><td>0.4434</td><td>1.08 (0.89–1.31)</td><td>0.075, 0.065</td><td>0.0597</td><td>1.17 (0.99–1.37)</td></tr><tr><td>rs6107516</td><td>C</td><td>PRNP</td><td>0.827, 0.752</td><td>2.07E-06</td><td>6.92E-06</td><td>1.58 (1.31–1.91)</td><td>0.804, 0.766</td><td>5.00E-04</td><td>1.25 (1.10–1.42)</td><td>0.811, 0.759</td><td>7.62E-09</td><td>1.35 (1.22–1.50)</td></tr><tr><td>rs6951643</td><td>A</td><td>GRM8</td><td>0.592, 0.507</td><td>6.73E-06</td><td>8.09E-06</td><td>1.41 (1.21–1.64)</td><td>0.573, 0.529</td><td>7.00E-04</td><td>1.19 (1.08–1.32)</td><td>0.578, 0.519</td><td>1.66E-08</td><td>1.27 (1.17–1.38)</td></tr><tr><td>rs9521699</td><td>A</td><td>COL4A2</td><td>0.194, 0.137</td><td>2.41E-05</td><td>3.84E-05</td><td>1.51 (1.25–1.83)</td><td>0.173, 0.173</td><td>0.9723</td><td>1.00 (0.88–1.15)</td><td>0.179, 0.157</td><td>0.0045</td><td>1.17 (1.05–1.31)</td></tr></tbody></table>

* Adjusted by Country of origin by PCA.

doi:10.1371/journal.pone.0123654.t002
<table><caption>Table 3. - rs6951643 pooled genotypes by country.</caption><thead><tr><th>Population</th><th>GG n (%)</th><th>AG n (%)</th><th>AA n (%)</th><th>A allele n (%)</th><th>G allele n (%)</th></tr></thead><tbody><tr><td><strong>CONTOLS</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>UK</td><td>364 (24.6)</td><td>741 (50.0)</td><td>377 (25.4)</td><td>1495 (50.4)</td><td>1469 (49.6)</td></tr><tr><td>NL</td><td>118 (22.4)</td><td>265 (50.3)</td><td>144 (27.3)</td><td>553 (52.5)</td><td>501 (47.5)</td></tr><tr><td>Spain</td><td>488 (22.3)</td><td>1097 (50.0)</td><td>608 (27.7)</td><td>2313 (52.7)</td><td>2073 (47.3)</td></tr><tr><td><strong>CASES</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>UK</td><td>28 (10.5)</td><td>155 (58.3)</td><td>83 (31.2)</td><td>321 (60.3)</td><td>211 (39.7)</td></tr><tr><td>NL</td><td>23 (18.1)</td><td>62 (48.8)</td><td>42 (33.1)</td><td>146 (57.5)</td><td>108 (42.5)</td></tr><tr><td>Germany</td><td>66 (16.7)</td><td>217 (54.8)</td><td>113 (28.5)</td><td>443 (55.9)</td><td>349 (44.1)</td></tr><tr><td>Italy</td><td>43 (15.0)</td><td>140 (48.8)</td><td>104 (36.2)</td><td>348 (60.6)</td><td>226 (39.4)</td></tr><tr><td>Australia</td><td>10 (20.8)</td><td>26 (54.2)</td><td>12 (25.0)</td><td>50 (52.1)</td><td>46 (47.9)</td></tr><tr><td>France</td><td>23 (15.3)</td><td>77 (51.3)</td><td>50 (33.3)</td><td>177 (59.0)</td><td>123 (41.0)</td></tr><tr><td>Spain</td><td>35 (17.2)</td><td>110 (54.2)</td><td>58 (28.6)</td><td>226 (55.7)</td><td>180 (44.3)</td></tr><tr><td>Austria</td><td>11 (20.0)</td><td>27 (49.1)</td><td>17 (30.9)</td><td>61 (55.5)</td><td>49 (44.5)</td></tr><tr><td><strong>ALL CONTOLS</strong></td><td>970 (23.1)</td><td>2103 (50.0)</td><td>1129 (26.9)</td><td>4361 (51.9)</td><td>4043 (48.1)</td></tr><tr><td><strong>ALL CASES</strong></td><td>239 (15.6)</td><td>814 (53.1)</td><td>479 (31.3)</td><td>1772 (57.8)</td><td>1292 (42.2)</td></tr><tr><td><strong>HapMap Frequencies (CEU)</strong></td><td>27 (23.9)</td><td>59 (52.2)</td><td>27 (23.9)</td><td>113 (0.50)</td><td>113 (0.50)</td></tr></tbody></table>

doi:10.1371/journal.pone.0123654.t003

were 12 tagging GRM8 (S2 Table). Pathway analysis using IPA showed that glutamate receptor signalling was the most over-represented canonical pathway within our top results (p-value = 8.01x10⁻⁵) (Fig 1). Analysis with ALIGATOR software showed (S3 Table) that the GO category glutamate receptor activity was the 4th most represented amongst our top results. Three genes from the GO category "glutamate receptor pathway" (S2 Table) were within our top 100 SNPs (GRM8, GRIN2B and GRIA1). Fine mapping in search for functional variants linked to rs6951643 was attempted sequencing the 11 exons of GRM8 gene and the corresponding intronic flanking regions in 96 sCJD patients. We found 8 intronic SNPs (rs73231278, rs62468898, rs1008274, rs111546739, rs17685327, rs6951643, rs2074012, rs35648111) and one exonic non-coding SNP (rs34182595). Functional prediction analysis of these genetic variants by FuncPred or RESCUE-ESE showed no indication of regulatory implications. In the Gtex eQTL database we found no eQTL associated to any of the GRM8 SNPs. There is no co-expression of GRM8 and PRNP in blood eQTLs in Genenetwork. However, based on PITA prediction we found that rs34182595, which is an insertion/deletion SNP in

<table><caption>Table 4. Meta-analysis by country of origin.</caption><thead><tr><th rowspan="2">MarkerName</th><th colspan="3">Effect estimates of discovery analysis corrected for population substructure by PCA</th><th colspan="4">Meta-analysis of effect estimates of replication analysis stratified by country of origin</th><th colspan="3">Meta-analysis of discovery and replication</th></tr><tr><th>β</th><th>SE</th><th>P-value</th><th>β</th><th>SE</th><th>P-value</th><th>β Direction*</th><th>β</th><th>SE</th><th>P-value</th></tr></thead><tbody><tr><td>rs6107516</td><td>0.4243</td><td>0.0944</td><td>6.92E-06</td><td>0.2327</td><td>0.0637</td><td>0.000262</td><td>+++++</td><td>0.2927</td><td>0.0528</td><td>3.00E-08</td></tr><tr><td>rs17115017</td><td>-0.5363</td><td>0.1453</td><td>0.000223</td><td>-0.1420</td><td>0.1061</td><td>0.181</td><td>-----+</td><td>-0.2792</td><td>0.0857</td><td>0.001123</td></tr><tr><td>rs9521699</td><td>-0.4087</td><td>0.0993</td><td>3.84E-05</td><td>-0.0233</td><td>0.0723</td><td>0.7468</td><td>+---+</td><td>-0.157</td><td>0.0585</td><td>0.007242</td></tr><tr><td>rs6951643</td><td>-0.3472</td><td>0.0778</td><td>8.09E-06</td><td>-0.1033</td><td>0.0545</td><td>0.05821</td><td>++++</td><td>-0.1837</td><td>0.0447</td><td>3.91E-05</td></tr><tr><td>rs10061929</td><td>-0.3737</td><td>0.0917</td><td>4.58E-05</td><td>-0.0289</td><td>0.0701</td><td>0.6801</td><td>-----+</td><td>-0.1562</td><td>0.0557</td><td>0.005048</td></tr></tbody></table>

doi:10.1371/journal.pone.0123654.t004

<footnote>* Concordance between discovery population βs and each of the five replication populations; PCA, Principal Component Analysis; SE, standard Error.</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig 1. Canonical pathway analysis performed by IPA software including those genes tagged by SNPs with p-value <0.001 in our GWAs analysis.**
The horizontal axis represents the pathways identified. The ratio (vertical axis, right) is calculated by the numbers of genes in a given pathway that meet cutoff criteria, divided by total numbers of genes that make up that pathway. The orange line stands for the threshold above which there are statistically significantly (by default P<0.05). The vertical axis (left) shows the -log of the p-value calculated based on Fisher's exact test.</figcaption></figure>

doi:10.1371/journal.pone.0123654.g001

exon 11 in partial LD with rs6951643 ($D' = 0.41; r^2 = 0.16$), is located in a target site for microRNAs (miR-103 and miR-107), a finding that has been previously reported [7].

We then performed an immunohistochemical study in brain samples from 48 sCJD patients. Immunoreactivity for mGluR8 was observed in neurons and microglial cells (both in the white and grey matter) While in neurons we did not observe obvious differences we found a trend towards higher combined score of mGluR8 immunoreactivity in microglia related to the A allele at rs6951643 (Fig 2); however, these differences were not statistically significant (ordinal regression A-carriers versus A-non-carriers adjusting by c129 and levels of microglia in temporal cortex; p-value = 0.093). S3 Fig shows semi-quantitative assessments of mGluR8 expression in microglia (0 to 3) in the temporal region across the three rs6951643 genotypes. We found no association between rs6951643 genotypes and patient's age of disease onset or disease duration.

## Discussion

In this study we report a non-PRNP genetic risk variant for sCJD, GRM8. Moreover, building on findings from previous studies [8], pathway analyses yielded glutamate receptor signalling as one of the main pathways linked to sCJD pathogenesis.

A previous GWAS of prion diseases employed 1,259 sCJD samples, in addition to other disease subtypes and 6,015 shared controls [6]. In the sCJD sub-group and in a meta-analysis including all prion cases only variants in PRNP were found to be significantly associated with disease risk. Several SNPs outside PRNP were identified, but in contrast to our multi-national disease cohort the association with sCJD was not homogeneous across the different geographical groups. Our main novel finding (rs6951643) was not statistically significant in that analysis. One possible explanation for this discrepancy is the fact that both studies have limited power to detect small effects due to the overall restricted case numbers included for study and therefore a SNP with a relatively modest OR of 1.27 might not be detected by all analyses. The experience and insights drawn from GWAS performed in more prevalent disorders, such as Alzheimer's disease, shows that increasing sample size and performing meta-analysis might be essential in order to discover and confirm, relatively small-effect genetic risk factors. However, it is worth re-iterating that in our study, the association with rs6951643 was consistent across
<figure><img src="image_6.png" /><figcaption>Fig 2. Representative photomicrographs of a sCJD case with AA (a, c) and a sCJD case with GG (b, d) rs6951643 genotype showing differences in mGluR8 immunostaining of microglial cells in the temporal cortex (a, b) and temporal white matter (c, d). Bar represents 50 µm for all images.</figcaption></figure>

doi:10.1371/journal.pone.0123654.g002

the geographically diverse sCJD population tested, showing an over-representation of the A allele in cases, and across the different genotyping methods used.

A limitation of our primary analysis was the fact that we were not able to adjust by country of origin at the replication stage. In order to overcome this caveat we performed a meta-analysis including control data from the *1000 Genomes Project* and imputed data from the Rotterdam Study. In the meta-analysis of the stratified analysis by country of origin the association of both SNPs became less significant suggesting that some signal might be caused by population stratification. Despite the fact that both variants were nominally significant in the meta-analysis of the replication data, the *PRNP* variant was the only genome wide significantly associated to sCJD in the pooled meta-analysis. However, this meta-analysis is not perfect either, as ideally the cases and controls should originate from the same population and in our analysis we attempted to match on country of origin (see S4 Table) and this approach is statistically less powerful because of the uneven distribution of cases and controls, as might reflect the fact that the *PRNP* SNP p-value is also less significant in the meta-analysis than in the primary analysis. Still we acknowledge that the results are a call for caution and the association between *GRM8* genetic variants and sCJD should be confirmed in larger analysis with extensive control for population stratification.
GRM8 encodes for mGluR8, a protein that belongs to the metabotropic glutamate receptor family, a family that has recently been linked to the transduction of physiological and cytotoxic signals mediated by PrP<sup>C</sup>. MGluR1 and mGluR5 have been shown to interact with PrP<sup>C</sup>, with such associations appearing important for promoting neurite outgrowth [9, 10]. Additionally, a recent study demonstrated that mGluR5 coupled with PrP<sup>C</sup> mediated the cellular toxicity of soluble β-amyloid oligomers [10]. Further, in the APPswe/PS1dE9 Alzheimer's disease mouse model, altered PrP<sup>C</sup> processing and a selective increase in cortical mGluR1 expression has been reported. The authors hypothesized that complex processing of PrP<sup>C</sup> in connection with mGluR1 over-expression is triggered by β-amyloid peptides [11]. In the setting of accumulation of PrP<sup>Sc</sup>, our immunohistochemical assessment of a limited sample of sCJD patients found that carriers of the risk allele at rs6951643 tended to have higher mGluR8 expression in microglial cells compared to non-carriers. Although our functional prediction analysis for rs6951643 was uninformative, it is of interest that a nearby variant in partial LD (rs34182595) is located in a micro-RNA target site and could influence gene expression.

In conclusion, our study has detected a *GRM8* genetic variant as a suggestive marker for sCJD risk mapping outside the *PRNP* region. Our findings provide evidence supporting a role for glutamate receptor signalling pathways in sCJD susceptibility. The involvement of glutamate receptor pathway in sCJD should be addressed in future studies in order to provide new insights in its pathogenesis. Our results underscore the importance of increasing sample sizes in future studies in order to augment the likelihood of detecting additional non-*PRNP* genetic risk factors for sCJD.

# Methods

## Ethics statement

The present study was conducted according to the revised Declaration of Helsinki and Good Clinical Practice guidelines. A signed informed consent to participate in genetic research was obtained from all participants or patients' relatives. The study was approved by Comité de Ética de la Investigación y de Bienestar Animal (CEIyBA), National Health Institute Carlos III and Comité de Ética de Ensayos clínicos de Galicia and Comité de Ética de la Fundación Pública Gallega de Medicina Genómica (Servicio Gallego de Salud, SERGAS) (Spain); the Lothian Health Board, Lothian Research Ethics Committee (reference MCO/103/90) (UK); the Ethik-Kommission der Universitätsmedizin Göttingen (No. 9/6/0) (Germany); the Ethic Committee of the Istituto Superiore di Sanità (CE-ISS 09/266) (Italy); the Medical Ethics Committee of The Erasmus MC and the review board of The Netherlands Ministry of Health Welfare and Sports (The Netherlands); the Ethics Committee of the Medical University of Vienna (396/2011) (Austria); the Human Research Ethics Committee, based at The University of Melbourne (941450) (Australia).

## Populations and study design

All patients were ascertained by National CJD Surveillance Centers. Only definite or probable sCJD cases, according to accepted classification criteria, were included [12]. Table 5 summarizes the study design and populations included after quality control. All cases and controls were of Caucasian origin. Legal representatives gave written informed consent, and all samples were taken in accordance with the Helsinki declaration.
<table><caption>Table 5. Study Design and populations tested.</caption><thead><tr><th rowspan="2">STAGES</th><th colspan="2">POPULATION EFFECTIVELY ASSESSED*</th></tr><tr><th>sCJD CASES (N)</th><th>CONTROLS (N)</th></tr></thead><tbody><tr><td><strong>Stage 1: Discovery</strong></td><td>Germany (113)</td><td></td></tr><tr><td>Genotyping method: Affymetrix 500k array</td><td>UK (269)</td><td>UK WTCCC controls (1482)</td></tr><tr><td>Number of SNPs effectively assessed 279389</td><td>Netherlands (52)</td><td>Netherlands-RS controls (457)</td></tr><tr><td>Number of SNPs selected for replication: 23</td><td><strong>TOTAL sCJD (434)</strong></td><td><strong>TOTAL controls (1939)</strong></td></tr><tr><td><strong>Stage 2: Replication in independent sCJD Cases</strong></td><td>Germany (284)</td><td></td></tr><tr><td>Genotyping method: Sequenom iPLEX GOLD</td><td>Netherlands (76)</td><td></td></tr><tr><td>Number of SNPs effectively assessed: 22</td><td>Italy (292)</td><td></td></tr><tr><td>Number of SNPs replicated: 5</td><td>Australia (48)</td><td></td></tr><tr><td></td><td>France (150)</td><td></td></tr><tr><td></td><td>Spain (203)</td><td></td></tr><tr><td></td><td>Austria (56)</td><td></td></tr><tr><td></td><td><strong>TOTAL sCJD (1109)</strong></td><td></td></tr><tr><td><strong>Stage 3: Replication in independent Controls</strong></td><td></td><td>Spain-USC (2193)</td></tr><tr><td>Genotyping method: Sequenom iPLEX GOLD</td><td></td><td>Netherlands-RS (71)</td></tr><tr><td>Number of SNPs effectively assessed: 5</td><td></td><td><strong>TOTAL controls (2264)</strong></td></tr><tr><td>Number of SNPs replicated: 2</td><td></td><td><strong>TOTAL controls (4203)</strong></td></tr><tr><td><strong>POOLED ANALYSIS</strong></td><td><strong>TOTAL sCJD (1543)</strong></td><td></td></tr><tr><td></td><td></td><td>Netherlands-RS (6192)</td></tr><tr><td></td><td></td><td>British in England and Scotland-10KG (137)</td></tr><tr><td><strong>Controls added in analysis stratified by country</strong></td><td></td><td>Toscani in Italy-10KG (385)</td></tr><tr><td></td><td></td><td>Iberian populations in Spain-10KG (217)</td></tr><tr><td></td><td></td><td>Utah residents with Northern and Western European ancestry-10KG (290)</td></tr><tr><td></td><td><strong>TOTAL sCJD (1543)</strong></td><td><strong>TOTAL controls (11424)</strong></td></tr></tbody></table>

*Number of samples effectively genotyped that passed quality control
sCJD, sporadic Creutzfeldt-Jakob disease; WTCCC, Welcome Trust Case Control Consortium; RS, Rotterdam Study; USC, University of Santiago de Compostela; 10KG, 1000 Genome Project.

doi:10.1371/journal.pone.0123654.t005

Genotyping and quality control

At the discovery stage (stage one) patient's samples were characterized using the Affymetrix 500k array at the USC node of the Spanish National Genotyping Center. Out of the initially available 554 sCJD samples there was not enough DNA in 61 (36 samples from the UK and 25 from Germany), leaving 493 for genotyping. Genotypes were called first by the Dynamic Model (DM) and those samples with an overall call rate >93% were subsequently called by a Bayesian Robust linear model with Mahalanobis distance (BRLMM). Samples with BRLMM call rates <95% and with call and discrimination rates of the Modified Partitioning Around Medoids algorithm (MDR-MCR) >10% were excluded (31 patient samples out of 493 did not reach these thresholds, leaving in total 462 available for the analysis). Genotype and quality control of the WTCCC controls have been described elsewhere [13]. The other set of controls (Rotterdam Study) were genotyped at the Genetic Laboratory, Department of Internal Medicine at Erasmus Medical Center (Rotterdam) also using the Affymetrix 500K array following same protocols as those used for cases. After quality control, 28 patients and 25 controls were excluded resulting in a final discovery population of 434 sCJD patients and 1,939 controls. Replication was attempted in independent sCJD patients (n = 1,201) (stage two) and control series
(n = 2,264) (stage three) with the Sequenom iPLEx GOLD platform. SNPs for replication were selected as follows: after excluding outliers, we chose the 10 SNPs with lowest p-values (when a cluster of SNPs was in LD tagging the same gene we selected only the one with lowest p-value) and the remainder from the top 100 when a) they were in linkage disequilibrium with other top 100 SNPs or b) agreed in direction with the previous sCJD GWAS [6], or c) were connected to pathways of interest based on our previous study [14] (phosphatidylinositol) or based on the pathway analyses performed with the current data (glutamate receptors).

## Statistical analysis of genetic data

The statistical analyses were conducted using GenABEL [15]. For the individual SNP analysis, we excluded those SNPs: 1) with call rates <98% within either group (n = 219,182); 2) with minor allele frequency<0.01 (n = 1); 3) with controls not in Hardy-Weinberg equilibrium (False Discovery Rate for unacceptably high individual heterozygosity<1%) (n = 1,300). The quality control further included a check of the sex chromosomes against the reported sex and unexpected sample duplicates (Identity by State, IBS>95%). All quality control was performed with the 'check.marker' function of GenABEL. After quality control, 279,389 out of 499,872 SNPs were included in the analyses. In the discovery analyses (stage one), we conducted a one degree of freedom additive score test between cases and controls coding the presence of the minor genotype 0 for non-carriers, 1 for heterozygous and 2 for homozygous carriers. We controlled subpopulation structure using principal components analysis (PCA). In brief, we selected a random set of 10,000 SNPs and calculated a genomic kinship matrix using pair wise identity by state statistic (function 'ibs'). We then performed PCA analysis (function 'cmdscale'), and adjusted our analysis by the three main IBS matrix principal components (function 'mlreg'). We calculated the genomic inflation factor lambda (λ) for both autosomal and X chromosome SNPs. In the replication phase (stages two and three), we compared allele frequencies in cases and controls by a Chi-squared test implemented in Haploview [16]. We used a Bonferroni correction to adjust for multiple testing, setting the threshold for significance to a p-value of 0.0023 (0.05/22 SNPs successfully genotyped). We combined the gene discovery and replication series of patients for all validated SNPs and used as the criterion for genome wide significance a p-value of 5x10⁻⁸. We attempted to adjust for country of origin after stage three. We analysed the data from stage two and three with genotypes of the *1000 Genomes Project* [17] and imputed genotypes from the *Rotterdam Study* [18]. We extracted the SNPs from the *1000 Genomes* reference set excluding the children and other family members (Version: phase I v3). Additionally we imputed the same SNPs from the *Rotterdam Study* (imputations of the same phase I v3 of the *1000 Genomes*). Imputed genotypes had a very high quality score (R² > = 0.99). Next we paired our cases and controls by country of origin (S4 Table) and calculated effect estimates of the SNPs per country. The effect estimates of the SNPs from the discovery stage (adjusted by PCAs) and the effect estimates from the replication were meta-analysed using an inverse variance weighted meta-analysis (METAL version released 2011-03-25). [19]. We performed pathway analyses using ALIGATOR, a method for studying groups of genes by testing for over-representation of members of those groups within lists of genes containing significantly associated SNPs from GWA studies [20]. This analysis used 410 SNPs with a p-value for significance of <0.001 out of 115,565 within-gene SNPs from a total of 279,389 available. We found these SNPs were associated with 94 genes from a total of 13,092 genes with GO annotation covered in this study. Additionally we used Ingenuity Pathway Analysis (IPA) www.ingenuity.com) to determine the functional pathways in the genes tagged by our top ranked SNPs. We again selected those genes tagged by SNPs with p-value <0.001 and selected the canonical pathway analysis implemented in IPA software.
## Sequencing

We sequenced the 11 exons and intronic flanking areas of the *GRM8* gene in 96 sCJD patients using specific primers (S5 Table). The amplification reactions were carried out with 25 ng of genomic DNA and 0.5 units of Taq DNA Polymerase in a volume of 25 µl. The final concentrations of other reactants were: 1x Taq DNA Polymerase Buffer, 0.1 mM dNTPs and 0.1 µM of each primer. The final concentration of MgCl₂ was 2 mM for the amplification of exons 1–7 and 10; 4 mM for exons 8 and 9a; 1.5 mM for exon 9b and 1 mM for exon 11. Additionally, the amplification reactions of exons 1 and 8 were supplied with DMSO 10% and 5%, respectively. The PCR cycling conditions were as follows: initial denaturation at 96°C for 3 min followed by 30 cycles of 96°C for 30 s, annealing temperature (see S5 Table) for 30 s and extension temperature of 72°C for 1 min with a final extension at 72°C for 10 min. A 2 µl aliquot of the amplification reaction was sequenced using 0.1µM of the above primers.

## Functional prediction analysis

Functional prediction analysis of genetic variants was performed by the use of FuncPred online software (http://snpinfo.niehs.nih.gov/snpinfo/snpfunc.htm) [21]. Exonic splicing enhancers were analyzed by RESCUE-ESE Web Server (http://genes.mit.edu/burgelab/rescue-ese/) [22] and PITA online software (http://genie.weizmann.ac.il/pubs/mir07/index.html) [23] was used to assess potential micro-RNA target sites. We searched for correlations between our genetic variants and eQTL in GTEX [24] and Genenetwork (http://genenetwork.nl:8080/GeneNetwork/) [25].

## Phenotypic correlations

We performed an immunohistochemical study in brain samples from 48 sCJD patients. The rs6951643 genotype distribution was as follows: 13 AA, 25 AG, and 9 GG. Sections from the hippocampal region CA1 sub-region, temporal cortex and white matter were immunostained for mGluR8. The intensity of mGluR8 (1:50; polyclonal rabbit antibody, Novus Biologicals, Cambridge, UK) immunostaining was evaluated using a scale of 0–3 (0: no; 1: weak; 2: moderate; 3: strong staining). The frequency of mGluR8 positive cells was scored semi-quantitatively using three categories: 1, <10%; 2, 10–50%; 3, >50% and was evaluated to provide information about the relative number of mGluR8 positive cells within the tissue. The product of these two values (intensity and frequency scores) was used to give the overall scores (total scores). In adjacent sections microglial activation using immunostaining for HLA-DP, DQ, DR (clone CR3/43, 1:00, monoclonal mouse), and spongiosis were also estimated semi-quantitatively. The analysis was performed (GGK) blinded to rs6951643 genotype. Ordinal regression was employed to test for association between rs6951643 genotypes and brain semi-quantitative expression of mGluR8. We adjusted by disease duration, PRNP M129V genotype, gliosis and spongiosis. We correlated clinical variables (age at death and patient's disease onset) with the presence of 0, 1 or 2 sCJD risk alleles. In order to assess rs6951643 influence on age at onset and disease duration we performed a time to event analysis using Cox regression adjusting by PRNPM129V genotype and other relevant factors like sex and country of origin.

## Supporting Information

S1 Fig. Q-Q plots for autosomal and X chromosome SNPs.
(TIF)
S2 Fig. Genotype clusters of the five SNPs studied in stage three with Sequenom iPLEx GOLD.

(TIF)

S3 Fig. Semi-quantitative assessments of mGluR8 expression in microglia (0 to 3) in the temporal region across the three rs6951643 genotypes in 48 sCJD patients.

(TIF)

S1 Table. Demographic and clinical features of cases.

(DOCX)

S2 Table. Top 100 SNPs after GWA analysis.

(XLS)

S3 Table. Pathway analysis with ALIGATOR.

(XLS)

S4 Table. Meta-analysis Sample population matching.

(DOCX)

S5 Table. Pairs of primers used for sequencing of the 11 exons and intronic flanking areas of GRM8 gene.

(DOCX)

S6 Table. GWA analysis results X chromosome SNPs.

(ZIP)

S7 Table. GWA analysis results autosomic SNPs.

(ZIP)

## Acknowledgments

The author thanks Simon Mead for sharing his data in order to select SNPs for replication. We also thank the Wellcome Trust Case Control Consortium (WTCCC) for the use of their control data. **UK.** The national UK CJD Surveillance Unit is grateful to clinicians, patients, and family members throughout the UK for a remarkable level of cooperation with the CJD surveillance program. **France.** We thank all members of the French National Surveillance Network for Creutzfeldt-Jakob disease and all physicians for case notification. **Italy.** We are very grateful to neurologists and neuropathologists, patients, and family members throughout Italy for their collaboration. **The Netherlands.** We thank Pascal Arp, Mila Jhamai, Marijn Verkerk, Lizbeth Herrera, and Marjolein Peters for their help in creating the GWAS database, and Karol Estrada and Maksim V. Struchalin for their support in creation and analysis of imputed data. The authors are grateful to the study participants, the staff from the Rotterdam Study, and the participating general practitioners and pharmacists. **Spain.** The authors thank Biobanco Valdecilla and IDIVAL for its support throughout this project thanks to J.M. Polo for their critical advice, thanks to Beatriz Sobrino and María Torres (CEGEN, Spanish National Genotyping Center) for all their work in the sample genotyping, and Weyma Notel for her help in editing the manuscript. **Australia.** The ANCJDR thanks all patients, families, clinicians and allied health personnel for their support and cooperation to facilitate this study.

## Author Contributions

Conceived and designed the experiments: CMvD RGW RSGK MTB YSA PSJ. Performed the experiments: AC GGK MC O. Calero. Analyzed the data: PSJ YSA CMvD SJvdL. Contributed
reagents/materials/analysis tools: MTB GGK MC AL JPB JLL SJC AB VL AP AC SJvdL FR AH O. Combarros JB AGU SH HB MP IZ RSGK RGW CMvD CP TS. Wrote the paper: PSJ MTB GGK MC O. Combarros MP SJC IZ RSGK RGW CMvD.

## References

1. Bueler H, Aguzzi A, Sailer A, Greiner RA, Autenried P, Aguet M, et al. Mice devoid of PrP are resistant to scrapie. Cell. 1993; 73:1339-47. PMID: 8100741

2. Alperovitch A, Zerr I, Pocchiari M, Mitrova E, de Pedro Cuesta J, Hegyi I, et al. Codon 129 prion protein genotype and sporadic Creutzfeldt-Jakob disease. Lancet. 1999; 353:1673-1674. PMID: 10335789

3. Bruce ME, Will RG, Ironside JW, McConnell I, Drummond D, Suttie A, et al. Transmissions to mice indicate that 'new variant' CJD is caused by the BSE agent. Nature. 1997; 389:498-501. PMID: 9333239

4. Bishop MT, Pennington C, Heath CA, Will RG, Knight RS. PRNP variation in UK sporadic and variant Creutzfeldt Jakob disease highlights genetic risk factors and a novel non-synonymous polymorphism. BMC Med Genet. 2009; 10:146. doi: 10.1186/1471-2350-10-146 PMID: 20035629

5. Lloyd SE, Mead S, Collinge J. Genetics of prion diseases. Curr Opin Genet Dev. 2013; 23:345-51. doi: 10.1016/j.gde.2013.02.012 PMID: 23518043

6. Mead S, Uphill J, Beck J, Poulter M, Campbell T, Lowe J, et al. Genome-wide association study in multiple human prion diseases suggests genetic risk factors additional to PRNP. Hum Mol Genet. 2012; 21:1897-906. doi: 10.1093/hmg/ddr607 PMID: 22210626

7. Duan S, Mi S, Zhang W, Dolanet ME. Comprehensive analysis of the impact of SNPs and CNVs on human microRNAs and their regulatory genes. RNA Biol. 2009; 6: 412-425. PMID: 19458495

8. Muller WE, Ushijima H, Schroder HC, Forrest JM, Schatton WF, Rytik PG, et al. Cytoprotective effect of NMDA receptor antagonists on prion protein (PrionSc)-induced toxicity in rat cortical cell cultures. Eur J Pharmacol. 1993; 15:246-261-7.

9. Beraldo FH, Arantes CP, Santos TG, Machado CF, Roffe M, Hajj GN, et al. Metabotropic glutamate receptors transduce signals for neurite outgrowth after binding of the prion protein to laminin gamma1 chain. FASEB J. 2011; 25, 265-279. doi: 10.1096/fj.10-161653 PMID: 20876210

10. Um JW, Kaufman AC, Kostylev M, Heiss JK, Stagi M, Takahashi H, et al. Metabotropic glutamate receptor 5 is a coreceptor for Alzheimer aβ oligomer bound to cellular prion protein. Neuron. 2013; 79:887-902. doi: 10.1016/j.neuron.2013.06.036 PMID: 24012003

11. Ostapchenko VG, Beraldo FH, Guimarães AL, Mishra S, Guzman M, Fan J, et al. Increased prion protein processing and expression of metabotropic glutamate receptor 1 in a mouse model of Alzheimer's disease. J Neurochem. 2013; 127:415-25. doi: 10.1111/jnc.12296 PMID: 23651058

12. WHO. Human transmissible spongiform encephalopathies. Wkly Epidemiol Rec. 1998; 47:361-365. PMID: 9844549

13. The Wellcome Trust, Case Control Consortium. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007; 447: 661-678. PMID: 17554300

14. Sanchez-Juan P, Bishop MT, Aulchenko YS, Brandel JP, Rivadeneira F, Struchalin M, et al. Genome-wide study links MTMR7 gene to variant Creutzfeldt-Jakob risk. Neurobiol Aging. 2012; 33:1487. doi: 10.1016/j.neurobiolaging.2011.10.011 PMID: 22137330

15. Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007; 23: 1294-1296. PMID: 17384015

16. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005; 21:263-5. PMID: 15297300

17. The 1000 Genomes Project Consortium. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012; 491: 56-65. doi: 10.1038/nature11632 PMID: 23128226

18. Hofman A, Darwish Murad S, van Duijn CM, Franco OH, Goedegebure A, Ikram MA, et al. The Rotterdam Study: 2014 objectives and design update. Eur J Epidemiol. 2013; 28: 889-926. doi: 10.1007/s10654-013-9866-z PMID: 24258680

19. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010; 26: 2190-2191. doi: 10.1093/bioinformatics/btq340 PMID: 20616382

20. Holmans P, Green EK, Pahwa JS, Ferreira MA, Purcell SM, Sklar P, et al. Gene ontology analysis of GWA study data sets provides insights into the biology of bipolar disorder. Am J Hum Genet. 2009; 85:13-24. doi: 10.1016/j.ajhg.2009.05.011 PMID: 19539887

21. Xu Z, Taylor JA. SNPinfo: integrating GWAS and candidate gene information into functional SNP selection for genetic association studies. Nucleic Acids Res. 2009; 37: W600-5. doi: 10.1093/nar/gkp290 PMID: 19417063
22. Fairbrother WG, Yeh RF, Sharp PA, Burge CB. Predictive identification of exonic splicing enhancers in human genes. Science. 2002; 297: 1007-13. PMID: 12114529

23. Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E. The role of site accessibility in microRNA target recognition. Nat Genet. 2007; 39:1278-84. PMID: 17893677

24. Lonsdale J, Thomas J, Salvatore M, Phillips R, Lo E, Shad S, et al. The genotype-tissue expression (GTEx) project. Nat Genet. 2013; 45:580-585. doi: 10.1038/ng.2653 PMID: 23715323

25. Williams RW, Mulligan MK. Genetic and molecular network analysis of behavior. Int Rev Neurobiol. 2012; 104:135-57. doi: 10.1016/B978-0-12-398323-7.00006-9 PMID: 23195314