Eur J Neurol. 2025;32:e16499. | 1 of 9 https://doi.org/10.1111/ene.16499 wileyonlinelibrary.com/journal/ene Received:7May2024 | Accepted:16September2024 DOI: 10.1111/ene.16499 ORIGINAL ARTICLE Understanding Parkinson disease in Spain: Genetic and clinical insights Pilar GómezGarre1,2 | Miguel MartínBórnez1 | Laura MuñozDelgado1,2 | Rafael DíazBelloso1,2,3 | María Teresa Periñán1,2,4 | Marta BonillaToribio1,2 | Dolores BuizaRueda1,2 | Daniel MacíasGarcía1,2 | Silvia Jesús1,2 | Astrid AdarmesGómez1,2 | Elena Ojeda1,2 | Antonio LuqueAmbrosiani1 | Sergio GarcíaDíaz1 | Rocío Pineda Sánchez1,2 | Fátima Carrillo1,2,3 | Pablo Mir1,2,3 1UnidaddeTrastornosdelMovimiento,ServiciodeNeurología,InstitutodeBiomedicinadeSevilla,IBiS/HospitalUniversitarioVirgendelRocío/CSIC/ UniversidaddeSevilla,Seville,Spain 2CentrodeInvestigaciónBiomédicaenRedsobreEnfermedadesNeurodegenerativas,InstitutodeSaludCarlosIII,Madrid,Spain 3DepartamentodeMedicina,FacultaddeMedicina,UniversidaddeSevilla,Seville,Spain 4CentreforPreventiveNeurologyUnit,WolfsonInstituteofPopulationHealth,QueenMaryUniversityofLondon,London,UK ThisisanopenaccessarticleunderthetermsoftheCreativeCommonsAttribution-NonCommercialLicense,whichpermitsuse,distributionandreproduction inanymedium,providedtheoriginalworkisproperlycitedandisnotusedforcommercialpurposes. ©2024TheAuthor(s).European Journal of NeurologypublishedbyJohnWiley&SonsLtdonbehalfofEuropeanAcademyofNeurology. Correspondence PilarGómez-GarreandPabloMir,Unidad de Trastornos del Movimiento, Instituto deBiomedicinadeSevilla(IBiS),Hospital UniversitarioVirgendelRocío/CSIC/ UniversidaddeSevilla,AvanzadaManuel Siurots/n,Seville41013,Spain. Email:mgomez-[email protected] and
[email protected] Funding information ConsejeríadeEconomía,Innovación, CienciayEmpleodelaJuntadeAndalucía, Grant/AwardNumber:CTS-7685,CVI- 02526andPY20_00896;Consejeríade SaludyBienestarSocialdelaJuntade Andalucía,Grant/AwardNumber:PE- 0 1 8 6 - 2 0 1 9 , P E - 0 2 1 0 - 2 0 1 8 , P I - 0 4 5 9 - 2 0 1 8 andPI-0471-2013;Grant/AwardNumber: RTC2019-007150-1,InstitutodeSalud CarlosIII(ISCIII)andco-fundedbythe EuropeanUnion,Grant/AwardNumber: PI14/01823,PI16/01575,PI18/01898, PI19/01576andPI21/01875 Abstract Background and purpose: Parkinsondisease(PD)isacomplexandheterogeneousneurodegenerativedisorderwithabroadspectrumofclinicalmanifestations,determinedbya complexinterplayofenvironmentalandgeneticfactors.Thisstudyaimedtoinvestigate geneticvariantsassociatedwithPDandassesstheirimpactonthediseasephenotype through genotype–phenotype correlations. Methods: Weemployedatargetedresequencingpaneltoanalyze27geneslinkedto PD in a cohort of 1185 PD patients from southern Spain. Variants were categorized basedontheAmericanCollegeofMedicalGeneticsandGenomicspathogenicitycriteria. Demographic and clinical data were also collected. Results: Amongthepatientsanalyzed,13.5%carriedpotentialdisease-causingpathogenicorlikelypathogenicvariantsin12differentgenes,indicatingsignificantgenetic heterogeneity.ThemostfrequentlyaffectedgeneswereLRRK2, PRKN, and GBA1(accountingfor 72.1%of positive cases).Sex-specific differenceswere observed, with a higherproportionoffemalepatientscarryingLRRK2variants.Differencesinageatonset andclinicalfeatureswerealsoobservedamongthedifferentmutatedgenes.Notably, variants in genes associated with atypical parkinsonism presented distinct clinical presentations,highlightingtheimportanceofgeneticfactorsinthedifferentialdiagnosis. Conclusions: OurstudyprovidesvaluableinformationonthegeneticlandscapeofPD anditsclinicalmanifestations.Theobservedgenotype–phenotypecorrelations, along withsex-specificdifferences,emphasizethecomplexityofPDpathogenesis,underliningtheimportanceofpersonalizedapproachestoPDdiagnosisandtreatment.Further
2 of 9 | GÓMEZ-GARRE et al. INTRODUCTION Parkinsondisease(PD)isacommonneurodegenerativedisease, clinically characterized by three cardinal motor symptoms: bradykinesia,resting tremor,andmusclerigidity.However,PDisa complex,heterogeneousdisorderwithabroadspectrumofclinical manifestations. Beyond motor symptoms, PD also presents numerousnonmotorsymptoms,oftenappearinglongbeforethe onsetofclassicmotor symptomsandsignificantly affecting the quality of life. Dementia is particularly prevalent, occurring in 83%ofpatientswithPD,withadiseasedurationofapproximately 20 years[1]. The risk of developing PD is determined by a complex interplayofenvironmentalandgeneticriskfactors[2], with age being themostimportantriskfactor.TheincidenceandprevalenceofPD increase almost exponentially with age, peaking at >80 years [3]. Likewise,sexisanestablishedriskfactor,withamale/femaleratio ofapproximately3:2[1]. Withinitsmultifactorialetiopathogenesis,PDisrecognizedas ahighlycomplexgeneticdisease.NumerousstudieshaveinvestigatedriskandprotectivefactorsofPD,withsignificantprogressin geneticresearchduetorapidadvancesinsequencingmethods.The geneticcontributiontoPDisexplainedbyawidespectrumofgeneticvariants,rangingfromcommonvariantsthatconfersusceptibilitytodevelopingPD,withamoderatetoweakeffectsize,torare variantforms,highlypenetrant,wherethepresenceofthevariant issufficienttocausethedisease[4]. To date, pathogenic variants associatedwithmonogenicPDhavebeendescribedin>20 genes, mostofwhicharehighlypenetrantandoftencauseearlyonsetor atypicalPD[4].Althoughithasbeenreportedthatonly5%–10%of PDcasesaremonogenic[5], these data are mainly based on individualsofEuropean/Whiteancestry.Thefrequency,penetrance,and clinicalimpactofgeneticvariantscanvaryacrossethnicities.Spain has a complex demographic history, characterized by continuous occupation since ancient times by various populations with diverse origins,includinganeight-centuryperiodofArabdomination(from MiddleEastandNorthAfrica),withacleargeneticimpact[6]. This study aimed to assess the presence and contribution of variantsinPD-relatedgenesinacohortofPDpatientsfromsouthernSpain,andtoestablishgenotype–phenotypecorrelations.We usedatargetedresequencingpanel,focusingon27genesassociatedwithPD,andanalyzedthedatausinganin-housepipeline.Our findingssignificantlycontributetounderstandingofPD,shedding lightontheeffectofgeneticvariantsondiseasedevelopmentand offeringimportantinsightsintodiagnosisandtreatmentofPDinour population. PATIENTS AND METHODS SeeAppendixS9. Study participants and clinical assessments Atotalof1185PDpatientswereincludedfrom2008.Thediagnoses wereperformedfollowingtheUnitedKingdomParkinson'sDisease SocietyBrainBankcriteria[7]forpatientsrecruitedupto2018and MovementDisorderSocietyclinicaldiagnosticcriteria[8]from2019 onward. Genetic analysis Genomic DNA was isolated from peripheral blood samples from eachparticipant.Acombinationofhigh-resolutionmeltinganalysis, multiplexligation-dependentprobeamplificationanalysis,andnext generation sequencing (NGS)-based targeted resequencing was applied(Figure S1).VariantswerecategorizedaccordingtotheinternationalguidelinesoftheAmericanCollegeofMedicalGenetics andGenomics(ACMG)[9](Figure S2).Onlyvariantsconformingto thegene'sinheritancepatternwereconsideredpotentialdisease- causing variants, and the genetic diagnosis was considered positive. Bioinformatics analysis for NGS results All NGS data were analyzed employing an in-house pipeline (AppendixS1). TobetterdefinetheroleofGBA1variantsinPD,variantsinthis gene were classified using the GBA1-PD browser (ht tps:// pdgen e t i c s . s h i n y a p p s . i o / G B A 1 B r o w s e r / )[10]. Inaddition,copynumbervariations(CNVs)werealsoanalyzed. Lossesweredefinedbyamaximumlog2ratioof−0.55andaminimumlog2ratioof0.4forgains[11]. Statistical analysis Associationanalysesof the variants and PD risk were performed withPLINKsoftware[12].Ararevariantburdenanalysiswasconductedwithconsistentsummarycounts-basedrarevariantburden test[13].Ascontrols,weusedtwosubcohorts:non-neuroandcontrolsintheNon-FinnishEuropean-SouthernEuropean(NFE-SEU) investigationsintogeneticinteractionsandpopulation-specificeffectsarewarrantedto enhanceourunderstandingofPDetiologyandimprovepatientcare. KEYWORDS clinicalinsights,geneticinsights,Parkinsondisease,targetedsequencing 14681331, 2025, 1, Downloaded from https://onlinelibrary.wiley.com/doi/10.1111/ene.16499 by Readcube (Labtiva Inc.), Wiley Online Library on [14/04/2025]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
| 3 of 9 PARKINSON DISEASE IN SPAIN cohort.WealsouseddatafromtheCIBERERSpanishVariantServer (CSVS;h t t p s : // c s v s . c l i n b i o i n f o s s p a . e s / )[14]. RESULTS Demographic characteristics Weanalyzed27genesassociatedwithPDinacohortof1185unrelatedPDpatients(Table 1).Amongthem,162(13.7%)reporteda familyhistoryofPD,and387(32.7%)hadearlyonsetPD(EOPD).The ageatonset(AAO)wasunavailableforthreepatients,buttheywere classifiedashavingEOPDbasedontheircurrentage(<50 years). Genetic spectrum Weidentifiedatotalof1014variantsinthecodingsequenceand spliceregionsofthe27analyzedgenes,allwithaminorallelefrequency(MAF)< 5%inthegeneralpopulation.Thesevariantswere categorizedaccordingtoACMGcriteriaasbenign(B),likelybenign (LB),variantofunknownsignificance(VUS),likelypathogenic(LP), andpathogenic(P).Amongtheidentifiedvariants,898werelocated incodingregionsand116inspliceregions(Table S1A and Figure S3A). Specifically,71variantswereclassifiedasP/LP,156asVUS,and788 asB/LB(Table S1B and Figure S3B).TheP/LPvariantswerefoundin 207patients(17.5%),affecting17genes,includingATP13A2, DCTN1, DNAJC6, FBXO7, LRP10, LRRK2, PRKN, PARK7, PINK1, PLA2G6, POLG, SMPD1, SPR, SYNJ1, VPS13A, VPS13C, and GBA1. CNVs were detected in the PRKN and SNCAgenes.Consideringtheinheritance patternofeachaffectedgene,potentiallydisease-causingP/LPvariantswereidentifiedin160patients(13.5%ofthecohort),involving12distinctgenes(DCTN1, DNAJC6, LRP10, LRRK2, PRKN, PINK1, POLG, SMPD1, SNCA, SPR, VPS13A, and GBA1).Amongthese160 patientswithpositivegeneticdiagnosis(PGD),sevencarriedP/LP variantsinmorethanonegene.Nodisease-causingvariantswere identifiedin960patients(80.9%),whowereconsideredtohavea negativegenetic diagnosis (NGD).Additionally, 65patients (5.5%; Figure 1a and Table S2),wereclassifiedashavinganunknowngeneticdiagnosis(UGD),ThisgroupincludedpatientscarryingVUS, those lacking a second variant in recessive genes, and those with GBA1variantsconsideredPDriskfactorsaccordingtotheGBA1-PD browser.Geneticanalysisresults foreachpatientareprovided in Table S3. Variantsinthreegenes(LRRK2, PRKN, and GBA1)accountedfor 72.1%ofthepositivecases.ThemostprevalentgenewasLRRK2, with61patients(37.9%)carryingP/LPvariants(including49with theG2019Svariant),followedbyPRKN(18%,including12patients withtwovariantsbutwithoutsegregationanalysis)andGBA1(18%; Figure 1b). We identified 29 patients (2.45%) carrying P/LP variants in GBA1. Severe variants were present in 10 patients (0.84%), mild variantsin14patients(1.18%),and unknown variantsinfour patients. Moreover, 21 patients (1.76%) carried risk variants. Three patients carried two variants in GBA1butdidnotexhibitsymptoms of Gaucher disease (GD). The most frequent P/LP variants were p.N409Sandp.L483P,andthemostfrequentvariantwastherisk variantp.E365K. ThedistributionofpatientswithPGD,NGD,orUGDacrossdifferentgroupsisshowninFigure 1c. Genotype–phenotype correlations Patients with PGD had a younger mean AAO (50.6 ± 12.4 years) compared to those with NGD (56.7 ± 11.6 years) or UGD (55.0 ± 12.3 years).TheAAOvariedamongdifferentmutatedgenes, withheterogeneityevenamongcarriersofvariantsinthesamegene (Figure 2).TheproportionofpatientswithafamilyhistoryofPDwas similarbetweenthosewithPGDandUGD(18.8%and18.5%,respectively),andhigherthaninthosewithNGD(12.5%;Figure 3a).A higherproportionofwomen(44.4%)wasnotedamongpatientswith PGD,particularlyamongthosecarryingP/LPvariantsintheLRRK2 gene(Figure 3b). Wealsoexploredtheassociationsbetweenthegenotypesand clinical characteristics in PGD patients (Table S4). Despite some missing clinical data, GBA1PDpatientsexhibitedahigherproportionofrapideyemovementsleepbehaviordisorder(RBD)andolfactorydysfunction(88%and68.4%,respectively),whereasalower frequencyofRBDwasobservedinLRRK2PDpatients(39.6%)and a lower frequency of olfactory dysfunction in PRKN PD patients (6.7%; Table S4). Cognitive impairment (CI) was more frequent amongpatientswithPGDrelatedtoGBA1, DCTN1, POLG, and SPR, withahigherprevalenceinmen(Table 2). Cohort Participants, nSex, men, n (%) Age, years, mean ± SD AAO, years, mean ± SD Family PD, n (%) Total 1185 707(59.7) 64.6 ± 11 55.8 ± 11.9 162(13.7) EOPD 387 242(62.5) 53.3 ± 9.2 42.5 ± 7.4 53(13.7) LOPD 798 465(58.3) 69.6 ± 7.8 62.3 ± 7.5 109(13.7) FamilyPD 162 89(54.9) 65.8 ± 10.2 56.0 ± 12.3 162(100) SporadicPD 1023 618(60.4) 64.4 ± 11.1 55.8 ± 11.9 0(0) Abbreviations:AAO,ageatonset;EOPD,earlyonsetPD(ageatonset≤ 50 years);LOPD,lateonset PD;PD,Parkinsondisease. TABLE 1 DemographicdataofthePD cohort. 14681331, 2025, 1, Downloaded from https://onlinelibrary.wiley.com/doi/10.1111/ene.16499 by Readcube (Labtiva Inc.), Wiley Online Library on [14/04/2025]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
4 of 9 | GÓMEZ-GARRE et al. Genotype–phenotype correlations in wellestablished PDrelated genes LRRK2 Dystoniawasobservedin68.4%ofLRRK2PDpatients,andresting tremorwasobservedin59.3%.Asubstantialmajorityexperienced motorfluctuation(91.8%)anddyskinesias(83.3%),withCIandRBD presentin31.0%and39.6%,respectively. A familyhistory of PD wasnotedonlyamongcarriersofthepathogenicG2019Svariant (Table S4). In our cohort, LRRK2 PD patients were predominantly women (60.7%).Thishigherproportionwasonlyobservedincarriersofthe G2019S variant, among whom 60% were women. No significant differencesinAAOwereobservedbasedonthespecificvariantcarried among LRRK2PDpatients,although23patients(37.7%)presented EOPD.Amongthese,fivepatientscarriedP/LPvariantsinothergenes. Onepatient(EP-40)carriedP/LPvariantsinboththeLRRK2 and DCTN1genes,withanAAOof51 yearsandabenigndiseasecourse without notable nonmotor symptoms. Twopatients(EP-10andEP-16)werefoundtocarrytheG2019S variantandhomozygousPRKNvariants,exhibitingclinicalfeatures similar to those associated with PRKNPD. Lastly, two patients (EP-215 and EP-223) carried VUSs and were therefore included in the UGD group. Their AAOs were 62 and54 years,respectively,andtheypresentedwithrestingtremor. Despitetheirrelativelybriefclinicalfollow-up,theirclinicaldisease course has been benign. FIGURE 1 Thefrequencyofcausative genes in patients initially diagnosed withParkinsondisease(PD).(a)Overall diagnosticyield.(b)Distributionof disease-causativegenesin160probands. (c)FeaturesoftheanalyzedPDcohort. EOPD,earlyonsetPD(ageatonset ≤ 50 years);LOPD,lateonsetPD(age at onset > 50 years);NGD,patients withanegativegeneticdiagnosis;PGD, patients with a positive genetic diagnosis; UGD,patientswithanunknowngenetic diagnosis. FIGURE 2 Boxplotsofageatonset underspecificgenes.Dotsrepresent the age at onset in each patient carrying putative causative variants in each gene. Horizontallinesrepresentthemedian onset.M,mild;NGD,Parkinsondisease (PD)patientswithanegativegenetic diagnosis;PGD,PDpatientswitha positivegeneticdiagnosis;RF,riskfactor; S,severe;T,totalcohortwithpathogenic or likely pathogenic variants in GBA1 and those with risk variants in GBA1;UGD, PDpatientswithanunknowngenetic diagnosis. 14681331, 2025, 1, Downloaded from https://onlinelibrary.wiley.com/doi/10.1111/ene.16499 by Readcube (Labtiva Inc.), Wiley Online Library on [14/04/2025]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
| 5 of 9 PARKINSON DISEASE IN SPAIN SNCA Twopatients(EP-114andEP-135)carriedduplicationsoftheentire SNCAgeneinheterozygosis.However,thecutoffpointscouldnot bedetermined.OnepatientpresentedwithalateAAO,whereasthe otherhadanearlydiseaseonset.ThepatientwiththeearlierAAO also carried a risk variant in the GBA1 gene. PRKN AlthoughmostPRKNPDpatientsexhibitedanearlyAAOof PD,14.3%presentedlateAAO(rangingfrom55to59 years), motor fluctuations and dyskinesias were common (82.1% and77.8%,respectively),andCIandhallucinations/illusions wererelativelyrare(25.9%and18.9%,respectively),with FIGURE 3 Bargraphsshowingfamily historyofParkinsondisease(PD)and genderdistributionacrossspecific genes.NGD,PDpatientswithanegative geneticdiagnosis;PGD,PDpatients withapositivegeneticdiagnosis;UGD, PDpatientswithanunknowngenetic diagnosis. TABLE 2 Cognitiveimpairmentinthepositivegeneticdiagnosiscohort. LRRK2 PRKN GBA1 GBA1, severe GBA1, mild DCTN1 POLG SMPD1 SPR n42 27 25 912 3 9 8 8 CI(%) 31.0 25.9 64 55.6 58.3 66.7 66.7 25 62.5 CI-F(%) 33.3 8.3 28.6 25 0 0 66.7 050 CI-M(%) 20 40.0 77.8 80 70 100 66.7 28.6 66.7 D P - C I 15 ± 7.9 16.6 ± 7.8 7.8 ± 6.0 4.8 ± 1.8 11.1 ± 7.7 7.5 ± 3.5 10.0 ± 5.0 7.5 ± 3.5 9.4 ± 2.9 D P - C I - F 15.8 ± 6.9 15 4.5 ± 3.5 2 0 0 11.5 ± 6.4 010 D P - C I - M 16.5 ± 13.4 16.8 ± 8.5 8.2 ± 6.3 5.1 ± 1.0 11.1 ± 7.7 7.5 ± 3.5 9.3 ± 5.1 7.5 ± 3.5 9.3 ± 3.3 Note:ThedevelopmentofCIwasconsideredthroughoutthefollow-up.Themeanfollow-uptimewas15.1 ± 8.5 years. Abbreviations:CI,cognitiveimpairment;CI(%),percentageofPGD-patientswithCI;CI-F(%),percentageofwomenwithPGDandCI;CI-M(%),percentage ofmenwithPGDandCI;DP-CI,diseaseprogressiontoCI(meanofyears ± SD);DP-CI-F,diseaseprogressiontoCI(meanofyears ± SD)inwomen;DP- CI-M,diseaseprogressiontoCI(meanofyears ± SD)inmen;n,sampleswithpositivegeneticdiagnosis(PGD)andCIdata;PGD,positivegeneticdiagnosis. 14681331, 2025, 1, Downloaded from https://onlinelibrary.wiley.com/doi/10.1111/ene.16499 by Readcube (Labtiva Inc.), Wiley Online Library on [14/04/2025]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
6 of 9 | GÓMEZ-GARRE et al. CI more prevalent among men (40%) than women (8.3%; Table 2). PINK1 Onepatient(EP-3)carriedLPmissensevariantsinputativecompound heterozygosityinthePINK1gene.ThispatientpresentedwithearlyAAO (49 years),motorfluctuations,dyskinesias,and dystonia.Additionally, althoughCIwasnotpresent,hallucinations/illusionsoccurred. GBA1 Wedetected13differentvariantsintheGBA1gene.Accordingto theireffectonGDandPD,fiveweredescribedassevere,threeas mild,andtwoasriskvariants,theirfrequencybeing4.9%inEOPD and3.1%inlateonsetPD.TheAAOofthosecarryingtheriskvariantsp.E365Kandp.T408Mwaslaterthantherest. Motorfluctuationsanddyskinesiasweremorefrequentinpatients with mild variants compared with those with severe variants (61.5%vs.50%,and69.2%vs.66.7%,respectively).Although1patient(EP-105)presentedwithtwoseverevariantsinGBA1, no signs of GD have been reported in his clinical history, suggesting that those variants are probably in cis. RBD and olfactory dysfunction were observed with higher frequency in patients carrying severe variants (87.5% and 71.4%, respectively) compared with those with mild variants (84.6% and 66.7%,respectively). CIwaspresentinapproximatelyhalfoftheGBA1patients(affectingahigherproportionofmalepatients),anddiseaseprogressiontoCIwaslowerinthepatientscarryingGBA1 severe variants compared with those with GBA1mildvariants(Table 2). Genotype–phenotype in genes linked with parkinsonism with atypical signs DCTN1 Sixpatients(0.5%ofthePDcohort)harboredanLPvariantinthe DCTN1 gene: c.414 + 1G > A (rs576198476). This variant displayed acombined annotation dependent depletion(CADD) valueof 33 andwasdescribedasLPaccordingtoACMGcriteria.Inourpopulation, DCTN1PDpatientspresentedameanAAOof59.5 ± 6.9 years (Figure 2),goodresponsetolevodopa,andmotorfluctuationsand dyskinesia.Mostofthesepatientsexperiencedparticularlysevere neuropsychiatricsymptoms,andthereforetheydidnottoleratehigh dosesofdopaminergicmedication.Otherfeatures,suchasrespiratory symptoms and weight loss, were absent in their clinical history. Additionally,onepatientcarriedboththeDCTN1c.414 + 1G > Avariant and the LRRK2G2019Svariant,showinganearlierAAOanda lackofneuropsychiatricsymptoms. DNAJC6 Onepatient(EP-47)withalateAAOharboredtwoLPvariantsin DNAJC6, presenting with a typical PD phenotype. Both variants were frameshift changes likely in compound heterozygosity but withoutbiallelicconfirmation. LRP10 Onepatient(EP-29)withnofamilyhistoryofPDcarriedanonsense variant in the LRP10gene.ThispatientpresentedalateAAO anddevelopedCIandhallucinations/illusions5 yearsafterdisease onset. POLG Seven different P/LP variants were identified in 11 patients with a wide range of AAOs. It is known that two ofthesevariants(p.Thr251Ileandp.Pro587Leu)arelinked andin cis. Therefore,all POLG patients were considered monoallelic. Motorfluctuationswerepresentin81.8%ofthePOLGPDpatients. Dyskinesias, CI, RBD, and olfactory dysfunction were also frequentlyobserved(70%,80%,70%,and71.4%,respectively).No patientwasspecificallydescribedashavingataxiaorperipheralneuropathy.Onepatient(EP-54)presentedprogressiveexternalophthalmoplegia(PEO). SMPD1 ElevendifferentP/LPvariantsinSMPD1weredetectedin15PDpatients,whofrequentlyexhibitedmotorfluctuationsandRBD(71.4% and75%,respectively).Dyskinesias(61.5%)andhyposmia(62.5%) were also common, although hallucinations/illusions were rarely reported(23.1%). SPR AnLPvariant(p.Asp69Glu;c.207C > G)intheSPR gene was detected ineightpatientswithPD.Motorfluctuationsanddyskinesiaswere highly prevalent in SPRPDpatients(87.5%and75%,respectively); however,dystoniawasobservedin50%.OnlyoneSPRPDpatient reportedafamilyhistoryofPD. VPS13A OnepatientcarriedanLPvariantinVPS13A;however,clinicalinformationforthisindividualisunavailable. 14681331, 2025, 1, Downloaded from https://onlinelibrary.wiley.com/doi/10.1111/ene.16499 by Readcube (Labtiva Inc.), Wiley Online Library on [14/04/2025]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
| 7 of 9 PARKINSON DISEASE IN SPAIN Gene burden analysis Asignificantassociationwithveryrarevariantswasfoundforthe DNAJC13 gene, and PRKNshowedasimilarpatterninthecontextof loss-of-functionvariantswithaCADD> 25.Furthermore,PLA2G6 and POLG showed associations withnoncommonloss-of-function variants(Table S5). DISCUSSION This studyaimedtodefinetheclinicalspectrum and implications ofgeneticvariantsacross27genesinaPDpopulationfromsouthernSpain.Ourfindingsrevealedasignificantgeneticheterogeneity,withpotentiallycausativeP/LPvariantsidentifiedin12genes (DCTN1, DNAJC6, LRP10, LRRK2, PRKN, PINK1, POLG, SMPD1, SNCA, SPR, VPS13A, and GBA1). The analysis of typical and atypical forms of monogenic PD revealed several notable differences. Patients with atypical formsexhibitedahigherAAOandareducedfamilyhistoryofPD. Furthermore,weobservedsex-specificdifferencesintheprevalence ofP/LPvariants,particularlyintheLRRK2gene.Ahigherproportion ofwomenwerefoundtocarryLRRK2 variants, driven particularly bytheG2019Svariant(with60%ofthecarriersbeingwomen)but notamongcarriersofotherLRRK2variants.Thisfindingalignswith previousstudiesindicatingafemalepredominanceamongG2019S variantcarriers[15].Althoughotherstudieshavenotreportedthis gendereffect[16],alargemeta-analysisin2018ofLRRK2-related clinical features reported a similar female predominance among G2019Scarriers[17]. LRRK2 emerged as the most prevalent gene, with a notable proportionofpatientscarryingtheG2019Svariant(4%ofcases), similartootherEuropeans[18] but slightly higher than that reportedinthebroaderSpanishpopulation[19].Clinically,LRRK2- associatedPDpatientsexhibitedfeaturessimilartothosewith sporadicPD.However,therelativelylowerprevalenceofCIand RBDsuggestspotentialdifferencesindiseaseprogressionand underlyingneuropathology.Theidentificationofafamilyhistory ofPDexclusivelyamongcarriersoftheG2019Svariantunderscorestherelevanceofspecificgeneticvariantsinfamilialforms ofPDandemphasizestheimportanceofgeneticcounseling[20, 21]. Inlinewithotherstudies,ourcohortshowsahigherfrequency of the deleterious p.Asn52fs variant in PRKN compared to other populations. This variant has been related with an intronic signal in PRKN,describedasanAAOmodifier,highlightingtheimportanceof studying diverse populations to better understand the genetic contributionstoPD[19]. Interestingly,nosignificantdifferencesinAAOwereobserved among LRRK2 variant carriers, although some patients with LRRK2 and PRKNvariantspresentedwithAAOsoutsidethetypicalrange. Thisreinforcesthenotionthatadditionalgeneticandenvironmental factorslikelymodulatediseasepresentation. Moreover,CIwaslesscommonandtendedtoappearlaterinpatientswithPGDassociatedwithLRRK2 or PRKN compared to those withPGDassociatedwithothergenes,suchasGBA1, DCTN1, POLG, and SPR. AlthoughPINK1-associatedPDhasbeenreportedasthesecond mostcommoncauseofautosomal-recessivePD,itwasveryrarein our cohort. TwopatientswerefoundtocarryduplicationsintheSNCA gene; however,theprecisebreakpointsoftheseduplicationscouldnotbe determinated.Interestingly,oneofthesepatientsalsocarriedthe p.T408MriskvariantinGBA1andhadamuchearlierAAO(38vs. 61 years),suggestingthatgeneticinteractionsbetweenthesegenes mayinfluencediseasecourse. VariantsintheGBA1genehavebeenclassifiedaccordingtotheir effectonGDandPD[10, 22].Contrarytopreviousfindings[23], our studyfoundthatpatientscarryingmildGBA1 variants had an earlier AAOcomparedtothosewithseverevariants.Additionally,twopatientscarryingbothaseverevariantandariskvariant(withoutevidenceofGD)presentedwithalaterAAO,suggestinganinteraction between these two variants. Variantsingenesassociatedwithatypicalparkinsonism,suchas LRP10, DCTN1, DNAJC6, POLG, SPR, and SMPD1, were linked to distinct clinical presentations and disease courses, highlighting the importanceofconsideringgeneticfactorsinthedifferentialdiagnosis ofparkinsoniansyndromes. Therefore, although mutations in POLG gene, inherited in an autosomaldominantorrecessivemanner,havebeenlinkedtoPEO [24],inourcohortonlyonepatientexhibitedit,carryingthep.Tyr955Cysvariation,whichhasbeenpreviouslyassociatedwithPEO [24]. This suggests a genotype–phenotype correlation. The LRP10genehasbeenproposedascausativeinautosomal- dominantPD,withorwithoutdementia[25, 26]. One patient in our cohort presented clinical characteristics consistent with LRP10- related PD,althoughthe absenceofa familyhistoryofPD could suggestlowpenetranceforthisvariant. DCTN1genevariantshavebeenlinkedtoPerrysyndrome(PS) [27], arareautosomal-dominant adultonset neurodegenerative disordercharacterizedbyatypicalparkinsonism,depression/apathy,weightloss,andrespiratorysymptoms[28].PathogenicmutationsinPSaretypicallyfoundwithinthecytoskeleton-associated protein glycine-rich(CAP-Gly) domain[29].However,ourstudy identifiedthec.414 + 1G > Avariantoutsidethisdomain.Although carriersofthisvariantwereclassifiedashavingaPGD,thereremains uncertainty about its actual involvement in PD development.Thisvariantwasalsopresentinhealthycontrolsfromthe CSVSpopulation(MAF = 0.0019),althoughwithslightlylowerfrequencycomparedwithourPDpopulation(MAF = 0.0025).Thus, moreevidenceisneededtoconfirmitspathogenicity.Ourfindings suggest that variants in the DCTN1geneoutsidetheCAP-GlydomaincouldhaveuniqueimplicationsinPD,possiblypresentinga riskfactor. HomozygousvariantsinDNAJC6havebeenlinkedtoautosomal- recessiveearlyonsetPD.Inourcohort,onepatientcarriedtwoLP 14681331, 2025, 1, Downloaded from https://onlinelibrary.wiley.com/doi/10.1111/ene.16499 by Readcube (Labtiva Inc.), Wiley Online Library on [14/04/2025]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
8 of 9 | GÓMEZ-GARRE et al. variants,withalaterAAOthanpreviouslyreported.Additionaltests arerequiredtoconfirmtheircausality. EightindividualscarriedanLPvariant(p.Asp69Glu;c.207C > G) in the SPR gene, suggesting a founder effect in our cohort. A population-specificroleforthisgenehasbeensuggested[30]. The SPR gene encodes sepiapterin reductase, a key protein in tetrahydrobiopterin(BH4)synthesis,acofactorcriticalinthesynthesisand secretionofneurotransmitterssuchasdopamine.BH4 deficiency hasbeenassociatedwithconditionscharacterizedbymovementdisordersandmentalretardation,amongothers[31].However,there isconflictingevidenceregardingtheroleofSPRinPDpathogenesis [31].DescribedasaVUSinClinVar,thevariantc.207C > Ghasbeen previously associated with autosomal-dominant dopa-responsive dystonia[32].Inourcohort,62%ofSPRpatientspresentedwithCI, andonly50%presentedwithdystonia,suggestingthatthisvariant maybeaPDriskfactorinourpopulation. HeterozygousvariantsinSMPD1(encodingacidsphingomyelinase)havebeenidentifiedasriskfactorsforPD[33]. In our cohort, 1.26%ofpatientscarriedP/LPvariantsinSMPD1, supporting its role inPD.Notably,halfofthesepatientsexhibitedRBD,underscoring theclinicalimportanceofthisgene. Actually,itiscriticaltoincreaseourknowledgeofgeneticrisks for PD in individuals of varied genetic backgrounds and increase theavailabilityofdetailedclinaldataassociatedwiththem.Inthis context,studyingourPDpopulationfromaregionwithhistorically uniquedemographiccharacteristicsinEuropeaddsvaluabledatato ourunderstandingofPDpathogenesis. Toourknowledge,thisisthelargeststudyofPDtodateconductedinsouthernSpain,providinginformationintothegenetic variants and their associated clinical features. However, we acknowledge several limitations, including the inability to establish a cleargenotype–phenotypecorrelationforsomegenesduetothe limitedsizeofthesampleswithPGD,affectingthosegenesand challengesrelatedtomissingclinicaldata.Inaddition,theuseofa targetedgenepanelrestrictedourabilitytoanalyzerecentlydiscovered genes such as RAB32[34]. Finally, it is important to note thattherelativelylowfrequencyofGBA1 pathogenic variants in our population could be due to technical challenges related to the existenceofahighlysimilarpseudogene.Thismayhaveledtoan underestimation ofthetrue prevalenceofGBA1variantsinPD patients. Despite these limitations, our study provides important insights intothegeneticbasisofPDanditsclinicalmanifestations.Theobserved genotype–phenotype correlations and sex-specific differenceshighlighttheimportanceofpersonalizedapproachestoPD diagnosis and treatment. CONCLUSIONS OurresultshighlightthecomplexityofgeneticcontributionstoPD phenotype and prognosis. Comprehensive genetic profiling and personalizedriskstratificationarecrucialforguidingprognosisand therapeutic decision-making in PD patients. There are numerous mechanismsassociatedwithPDpathogenesisandprogressionthat leadtoPD.Furtherinvestigationsintogeneticinteractionsandthe inclusionoffamilyhistoryandsex-specificdatawillbekeytobetter understandingPDpathogenesis. AUTHOR CONTRIBUTIONS Pilar GómezGarre:Conceptualization;supervision;formalanalysis; fundingacquisition;projectadministration;writing–originaldraft; writing – review and editing. Miguel MartínBórnez: Data curation; formalanalysis;methodology;software;writing–reviewandediting. Laura MuñozDelgado: Data curation; investigation; resources; writing – review and editing. Rafael DíazBelloso: Formal analysis; investigation; writing – review and editing. María Teresa Periñán: Formal analysis; investigation; writing – review and editing. Marta BonillaToribio: Methodology; writing – review and editing. Dolores BuizaRueda: Methodology; writing – review and editing. Daniel MacíasGarcía: Resources; writing – review and editing. Silvia Jesús: Resources; writing – review and editing. Astrid AdarmesGómez: Resources; writing – review and editing. Elena Ojeda: Resources; writing – review and editing. Antonio LuqueAmbrosiani: Resources; writing – review and editing. Sergio GarcíaDíaz:Writing–review and editing. Rocío Pineda Sánchez: Methodology. Fátima Carrillo: Resources. Pablo Mir: Conceptualization; supervision; funding acquisition; project administration; writing – review and editing; resources. ACKNOWLEDGMENTS Theauthorswouldliketothankthedonors,theUniversityHospital Virgen del Rocio, and the Biomedical Institute of Seville Biobank (Andalusian Public Health System Biobank and ISCIII-Red de BiobancosPT20/00069)forthehumanspecimensusedinthisstudy. FUNDING INFORMATION Thisworkwassupportedby theSpanish MinistryofScienceand Innovation (RTC2019-007150-1); the Instituto de Salud Carlos III (ISCIII)cofundedbytheEuropeanUnion(PI14/01823,PI16/01575, PI18/01898,PI19/01576,PI21/01875);theConsejeríadeEconomía, Innovación,CienciayEmpleodelaJuntadeAndalucía(CVI-02526, CTS-7685,PY20_00896);andtheConsejeríadeSaludyBienestar SocialdelaJuntadeAndalucía(PI-0471-2013,PE-0210-2018,PI- 0459-2018,PE-0186-2019).P.G.-G.wassupportedbytheNicolás Monardesprogram(C-0048-2017)ofAndalusianRegionalMinistry of Health. D.M.-G. was supported by the Juan Rodés program (JR22/00073) of ISCIII cofunded by the European Union FSE+ (ISCIII-FEDER).L.M.-D.wassupportedbytheRíoHortegaprogram (CM21/00051) of ISCII cofunded by the European Union FSE+ (ISCIII-FEDER).Thefundershadnoroleinthestudydesign,data collection and analysis, decision to publish, or preparation of the manuscript. CONFLICT OF INTEREST STATEMENT Theauthorsdeclarenoconflictofinterest. 14681331, 2025, 1, Downloaded from https://onlinelibrary.wiley.com/doi/10.1111/ene.16499 by Readcube (Labtiva Inc.), Wiley Online Library on [14/04/2025]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License
| 9 of 9 PARKINSON DISEASE IN SPAIN DATA AVAILABILITY STATEMENT Thedatathatsupportthefindingsofthisstudyareavailableinthe supplementarymaterialofthisarticle. ORCID Pilar GómezGarre https://orcid.org/0000-0002-0437-6182 Miguel MartínBórnez https://orcid.org/0000-0001-8221-3816 Rafael DíazBelloso https://orcid.org/0000-0003-4713-9248 Daniel MacíasGarcía https://orcid.org/0000-0002-4822-1529 Pablo Mir https://orcid.org/0000-0003-1656-302X REFERENCES 1. Kalia LV, Lang AE. Parkinson's disease. Lancet. 2015;386(9996):896-912. 2. PangSY,HoPW,LiuHF,etal.Theinterplayofaging,geneticsand environmentalfactorsinthepathogenesisofParkinson'sdisease. Transl Neurodegener.2019;8:23. 3. Driver JA, Logroscino G, Gaziano JM, Kurth T. Incidence and remaining lifetime risk of Parkinson disease in advanced age. Neurology.2009;72(5):432-438. 4. BlauwendraatC,NallsMA,SingletonAB.Thegeneticarchitecture ofParkinson'sdisease.Lancet Neurol.2020;19(2):170-178. 5. LunatiA,LesageS,BriceA.ThegeneticlandscapeofParkinson's disease. Rev Neurol (Paris).2018;174(9):628-643. 6. BycroftC,Fernandez-RozadillaC,Ruiz-PonteC,etal.Patternsof geneticdifferentiationandthefootprintsofhistoricalmigrationsin theIberianPeninsula.Nat Commun.2019;10(1):551. 7. GibbWR,LeesAJ.TherelevanceoftheLewybodytothepathogenesis of idiopathic Parkinson's disease. J Neurol Neurosurg Psychiatry.1988;51(6):745-752. 8. PostumaRB,BergD,SternM,etal.MDSclinicaldiagnosticcriteria forParkinson'sdisease.Mov Disord.2015;30(12):1591-1601. 9. RichardsS,AzizN,BaleS,etal.Standardsandguidelinesforthe interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and GenomicsandtheAssociationforMolecularPathology.Genet Med. 2015;17(5):405-424. 10. Parlar SC, Grenn FP, Kim JJ, Baluwendraat C, Gan-Or Z. ClassificationofGBA1variantsinParkinson'sdisease:theGBA1- P D Browser.Mov Disord.2023;38(3):489-495. 11. PughTJ,AmrSS,BowserMJ,etal.VisCap:inferenceandvisualization of germ-line copy-number variants from targeted clinical sequencingdata.Genet Med.2016;18(7):712-719. 12. PurcellS,NealeB,Todd-BrownK,etal.PLINK:atoolsetforwhole- genomeassociationandpopulation-basedlinkageanalyses.Am J Hum Genet.2007;81(3):559-575. 13. ChenW,WangS,TithiSS,EllisonDW,SchaidDJ,WuG.Ararevariant analysisframeworkusingpublicgenotypesummarycountstoprioritizedisease-predispositiongenes.Nat Commun.2022;13(1):2592. 14. Pena-Chilet M, Roldan G, Perez-Florido J, et al. CSVS, a crowdsourcing database of the Spanish population genetic variability. Nucleic Acids Res.2021;49(D1):D1130-D1137. 15. Orr-UrtregerA,ShifrinC,RozovskiU,etal.TheLRRK2G2019SmutationinAshkenaziJewswithParkinsondisease:isthereagender effect?Neurology.2007;69(16):1595-1602. 16. Gan-Or Z, Leblond CS, Mallett V, Orr-Urtreger A, Dion PA, Rouleau GA. LRRK2mutationsinParkinsondisease;asexeffect or lack thereof? A meta-analysis. Parkinsonism Relat Disord. 2015;21(7):778-782. 17. ShuL,ZhangY,PanH,etal.ClinicalheterogeneityamongLRRK2 variants in Parkinson's disease: a meta-analysis. Front Aging Neurosci.2018;10:283. 18. Healy DG, Falchi M, O'Sullivan SS, et al. Phenotype, genotype,andworldwidegeneticpenetranceofLRRK2-associated Parkinson's disease: a case–control study. Lancet Neurol. 2008;7(7):583-590. 19. Bandres-CigaS,AhmedS,SabirMS,etal.Thegeneticarchitecture of Parkinson disease in Spain: characterizing population-specific risk,differentialhaplotypestructures,andprovidingetiologicinsight. Mov Disord.2019;34(12):1851-1863. 20. TaymansJM,FellM,GreenamyreT,etal.PerspectiveonthecurrentstateoftheLRRK2field.NPJ Parkinsons Dis.2023;9(1):104. 21. RochaEM,KeeneyMT,DiMaioR,DeMirandaBR,Greenamyre JT. LRRK2 and idiopathic Parkinson's disease. Trends Neurosci. 2022;45(3):224-236. 22. StirnemannJ,BelmatougN,CamouF,etal.AreviewofGaucher disease pathophysiology, clinical presentation and treatments. Int J Mol Sci.2017;18(2):441. 23. Gan-Or Z, Amshalom I, Kilarski LL, et al. Differential effects of severe vs. mild GBA mutations on Parkinson disease. Neurology. 2015;84(9):880-887. 24. DaPozzoP,CardaioliE,RubegniA,etal.NovelPOLG mutations and variable clinical phenotypes in 13 Italian patients. Neurol Sci. 2017;38(4):563-570. 25. Quadri M, Mandemakers W, Grochowska MM, et al. LRP10 genetic variants in familial Parkinson's disease and dementia with Lewybodies:agenome-widelinkageandsequencingstudy.Lancet Neurol.2018;17(7):597-608. 26. ManiniA,StranieroL,MonfriniE,etal.ScreeningofLRP10 mutationsinParkinson'sdiseasepatientsfromItaly.Parkinsonism Relat Disord.2021;89:17-21. 27. FarrerMJ,HulihanMM,KachergusJM,etal.DCTN1 mutations in Perrysyndrome.Nat Genet.2009;41(2):163-165. 28. Mishima T, Fujioka S, Tomiyama H, et al. Establishing diagnostic criteria for Perry syndrome. J Neurol Neurosurg Psychiatry. 2018;89(5):482-487. 29. DulskiJ,KogaS,Liberski PP,etal.Perrydisease:expanding the genetic basis. Mov Disord Clin Pract.2023;10(7):1136-1142. 30. SharmaM,MaraganoreDM,IoannidisJP,etal.Roleofsepiapterin reductasegeneatthePARK3locusinParkinson'sdisease.Neurobiol Aging.2011;32(11):2108. 31. WuY,ChenP,SunL,etal.Sepiapterinreductase:characteristics and role in diseases. J Cell Mol Med.2020;24(17):9495-9506. 32. ShalashAS,RoslerTW,MullerSH,etal.c.207C>Gmutationin sepiapterinreductasecausesautosomaldominantdopa-responsive dystonia. Neurol Genet.2017;3(6):e197. 33. Alcalay RN, Mallett V, Vanderperre B, et al. SMPD1 mutations, activity,andalpha-synucleinaccumulationinParkinson'sdisease. Mov Disord.2019;34(4):526-535. 34. GustavssonEK,FollettJ,TrinhJ,etal.RAB32Ser71ArginautosomaldominantParkinson'sdisease:linkage,association,andfunctional analyses. Lancet Neurol.2024;23:603-614. SUPPORTING INFORMATION Additional supporting information can be found online in the SupportingInformationsectionattheendofthisarticle. How to cite this article: Gómez-GarreP,Martín-BórnezM, Muñoz-DelgadoL,etal.UnderstandingParkinsondiseasein Spain:Geneticandclinicalinsights.Eur J Neurol. 2025;32:e16499. doi:10.1111/ene.16499 14681331, 2025, 1, Downloaded from https://onlinelibrary.wiley.com/doi/10.1111/ene.16499 by Readcube (Labtiva Inc.), Wiley Online Library on [14/04/2025]. See the Terms and Conditions (https://onlinelibrary.wiley.com/terms-and-conditions) on Wiley Online Library for rules of use; OA articles are governed by the applicable Creative Commons License