SpadaHC: a database to improve the classification of variants in hereditary cancer genes in the Spanish population.

Accurate classification of genetic variants is crucial for clinical decision-making in hereditary cancer. In Spain, genetic diagnostic laboratories have traditionally approached this task independently due to the lack of a dedicated resource. Here we present SpadaHC, a web-based database for sharing variants in hereditary cancer genes in the Spanish population. SpadaHC is implemented using a three-tier architecture consisting of a relational database, a web tool and a bioinformatics pipeline. Contributing laboratories can share variant classifications and variants from individuals in Variant Calling Format (VCF) format. The platform supports open and restricted access, flexible dataset submissions, automatic pseudo-anonymization, VCF quality control, variant normalization and liftover between genome builds. Users can flexibly explore and search data, receive automatic discrepancy notifications and access SpadaHC population frequencies based on many criteria. In February 2024, SpadaHC included 18 laboratory members, storing 1.17 million variants from 4306 patients and 16 343 laboratory classifications. In the first analysis of the shared data, we identified 84 genetic variants with clinically relevant discrepancies in their classifications and addressed them through a three-phase resolution strategy. This work highlights the importance of data sharing to promote consistency in variant classifications among laboratories, so patients and family members can benefit from more accurate clinical management. Database URL: https://spadahc.ciberisciii.es/.

Case report: Identification of a novel variant p.Gly215Arg in the CHN1 gene causing Moebius syndrome.

Background: Moebius Syndrome (MBS) is a rare congenital neurological disorder characterized by paralysis of facial nerves, impairment of ocular abduction and other variable abnormalities. MBS has been attributed to both environmental and genetic factors as potential causes. Until now only two genes, PLXND1 and REV3L have been identified to cause MBS. Results: We present a 9-year-old male clinically diagnosed with MBS, presenting facial palsy, altered ocular mobility, microglossia, dental anomalies and congenital torticollis. Radiologically, he lacks both abducens nerves and shows altered symmetry of both facial and vestibulocochlear nerves. Whole-exome sequence identified a de novo missense variant c.643G>A; p.Gly215Arg in CHN1, encoding the α2-chimaerin protein. The p.Gly215Arg variant is located in the C1 domain of CHN1 where other pathogenic gain of function variants have been reported. Bioinformatic analysis and molecular structural modelling predict a deleterious effect of the missense variant on the protein function. Conclusion: Our findings support that pathogenic variants in the CHN1 gene may be responsible for different cranial congenital dysinnervation syndromes, including Moebius and Duane retraction syndromes. We propose to include CHN1 in the genetic diagnoses of MBS.

High Performance of a Dominant/X-Linked Gene Panel in Patients with Neurodevelopmental Disorders.

Neurodevelopmental disorders (NDDs) affect 2-5% of the population and approximately 50% of cases are due to genetic factors. Since de novo pathogenic variants account for the majority of cases, a gene panel including 460 dominant and X-linked genes was designed and applied to 398 patients affected by intellectual disability (ID)/global developmental delay (GDD) and/or autism (ASD). Pathogenic variants were identified in 83 different genes showing the high genetic heterogeneity of NDDs. A molecular diagnosis was established in 28.6% of patients after high-depth sequencing and stringent variant filtering. Compared to other available gene panel solutions for NDD molecular diagnosis, our panel has a higher diagnostic yield for both ID/GDD and ASD. As reported previously, a significantly higher diagnostic yield was observed: (i) in patients affected by ID/GDD compared to those affected only by ASD, and (ii) in females despite the higher proportion of males among our patients. No differences in diagnostic rates were found between patients affected by different levels of ID severity. Interestingly, patients harboring pathogenic variants presented different phenotypic features, suggesting that deep phenotypic profiling may help in predicting the presence of a pathogenic variant. Despite the high performance of our panel, whole exome-sequencing (WES) approaches may represent a more robust solution. For this reason, we propose the list of genes included in our customized gene panel and the variant filtering procedure presented here as a first-tier approach for the molecular diagnosis of NDDs in WES studies.

Elimination of lipaemic interference by high-speed centrifugation.

Introduction: In order to deliver high quality results, detection and elimination of possible analytical interferences, such as lipaemia, is crucial. The aim of this study is to evaluate the efficacy of high-speed centrifugation in eliminating lipaemic interference and to define own lipaemic index (LI) for the studied biochemical analytes. Materials and methods: Evaluated analytes were: albumin, alkaline phosphatase, alanine-aminotransferase (ALT), aspartate-aminotransferase (AST), calcium, creatinine, gamma-glutamyltransferase (GGT), glucose, phosphates, total proteins, urea and total bilirubin. Those analytes and LIs have been analysed in duplicate in the Roche Diagnostics-c8000 analyser in samples centrifuged at 3000 rpm/10 minutes in the SL16 (Thermo Scientific, Waltham, USA) centrifuge and according to an own high-speed centrifugation protocol (12,900 rpm/15 minutes) in the MicroCL17R (Thermo Scientific, Waltham, USA) centrifuge. Lipaemia has been measured in each sample. The efficiency of high-speed centrifugation is verified by the Wilcoxon test (P < 0.05). In cases where significant differences are observed, our own LI is calculated. For ALT and AST, it is verified by McNemar test (P < 0.05). For creatinine, both Wilcoxon and McNemar test were applied. Results: There were statistically significant differences in analyte concentration before and after high-speed centrifugation for: albumin, creatinine, GGT, glucose, phosphates, urea and total bilirrubin. Own LI is calculated. McNemar test shows statistically significant diferences in the proportion of delivered results before and after high-speed centrifugation in ALT, AST and creatinine. Conclusions: This study confirms the efficacy of high-speed centrifugation protocol for all the considered analytes, excepting calcium, alkaline phosphatase and total proteins.

De novo coding variants in the AGO1 gene cause a neurodevelopmental disorder with intellectual disability.

BACKGROUND: High-impact pathogenic variants in more than a thousand genes are involved in Mendelian forms of neurodevelopmental disorders (NDD). METHODS: This study describes the molecular and clinical characterisation of 28 probands with NDD harbouring heterozygous AGO1 coding variants, occurring de novo for all those whose transmission could have been verified (26/28). RESULTS: A total of 15 unique variants leading to amino acid changes or deletions were identified: 12 missense variants, two in-frame deletions of one codon, and one canonical splice variant leading to a deletion of two amino acid residues. Recurrently identified variants were present in several unrelated individuals: p.(Phe180del), p.(Leu190Pro), p.(Leu190Arg), p.(Gly199Ser), p.(Val254Ile) and p.(Glu376del). AGO1 encodes the Argonaute 1 protein, which functions in gene-silencing pathways mediated by small non-coding RNAs. Three-dimensional protein structure predictions suggest that these variants might alter the flexibility of the AGO1 linker domains, which likely would impair its function in mRNA processing. Affected individuals present with intellectual disability of varying severity, as well as speech and motor delay, autistic behaviour and additional behavioural manifestations. CONCLUSION: Our study establishes that de novo coding variants in AGO1 are involved in a novel monogenic form of NDD, highly similar to the recently reported AGO2-related NDD.

New genes involved in Angelman syndrome-like: Expanding the genetic spectrum.

Angelman syndrome (AS) is a neurogenetic disorder characterized by severe developmental delay with absence of speech, happy disposition, frequent laughter, hyperactivity, stereotypies, ataxia and seizures with specific EEG abnormalities. There is a 10-15% of patients with an AS phenotype whose genetic cause remains unknown (Angelman-like syndrome, AS-like). Whole-exome sequencing (WES) was performed on a cohort of 14 patients with clinical features of AS and no molecular diagnosis. As a result, we identified 10 de novo and 1 X-linked pathogenic/likely pathogenic variants in 10 neurodevelopmental genes (SYNGAP1, VAMP2, TBL1XR1, ASXL3, SATB2, SMARCE1, SPTAN1, KCNQ3, SLC6A1 and LAS1L) and one deleterious de novo variant in a candidate gene (HSF2). Our results highlight the wide genetic heterogeneity in AS-like patients and expands the differential diagnosis.

FBXO28 causes developmental and epileptic encephalopathy with profound intellectual disability.

Chromosome 1q41-q42 deletion syndrome is a rare cause of intellectual disability, seizures, dysmorphology, and multiple anomalies. Two genes in the 1q41-q42 microdeletion, WDR26 and FBXO28, have been implicated in monogenic disease. Patients with WDR26 encephalopathy overlap clinically with those with 1q41-q42 deletion syndrome, whereas only one patient with FBXO28 encephalopathy has been described. Seizures are a prominent feature of 1q41-q42 deletion syndrome; therefore, we hypothesized that pathogenic FBXO28 variants cause developmental and epileptic encephalopathies (DEEs). We describe nine new patients with FBXO28 pathogenic variants (four missense, including one recurrent, three nonsense, and one frameshift) and analyze all 10 known cases to delineate the phenotypic spectrum. All patients had epilepsy and 9 of 10 had DEE, including infantile spasms (3) and a progressive myoclonic epilepsy (1). Median age at seizure onset was 22.5 months (range 8 months to 5 years). Nine of 10 patients had intellectual disability, which was profound in six of nine and severe in three of nine. Movement disorders occurred in eight of 10 patients, six of 10 had hypotonia, four of 10 had acquired microcephaly, and five of 10 had dysmorphic features, albeit different to those typically seen in 1q41-q42 deletion syndrome and WDR26 encephalopathy. We distinguish FBXO28 encephalopathy from both of these disorders with more severe intellectual impairment, drug-resistant epilepsy, and hyperkinetic movement disorders.

Alcohol Induced Depression: Clinical, Biological and Genetic Features.

BACKGROUND: In clinical practice, there is the need to have clinical and biological markers to identify induced depression. The objective was to investigate clinical, biological and genetic differences between Primary Major Depression (Primary MD) and Alcohol Induced MD (AI-MD). METHODS: Patients, of both genders, were recruited from psychiatric hospitalisation units. The PRISM instrument was used to establish the diagnoses. Data on socio-demographic/family history, clinical scales for depression, anxiety, personality and stressful life events were recorded. A blood test was performed analysing biochemical parameters and a Genome Wide Association Study (GWAS) to identify genetic markers associated with AI-MD. RESULTS: A total of 80 patients were included (47 Primary MD and 33 AI-MD). The AI-MD group presented more medical comorbidities and less family history of depression. There were differences in traumatic life events, with higher scores in the AI-MD (14.21 ± 11.35 vs. 9.30 ± 7.38; p = 0.021). DSM-5 criteria were different between groups with higher prevalence of weight changes and less anhedonia, difficulties in concentration and suicidal thoughts in the AI-MD. None of the genetic variants reached significance beyond multiple testing thresholds; however, some suggestive variants were observed. CONCLUSIONS: This study has found clinical and biological features that may help physicians to identify AI-MD and improve its therapeutic approach.

A New Risk Variant for Multiple Sclerosis at 11q23.3 Locus Is Associated with Expansion of CXCR5+ Circulating Regulatory T Cells.

Genome-wide association studies and meta-analysis have contributed to the identification of more than 200 loci associated with multiple sclerosis (MS). However, a proportion of MS heritability remains unknown. We aimed to uncover new genetic variants associated with MS and determine their functional effects. For this, we resequenced the exons and regulatory sequences of 14 MS risk genes in a cohort of MS patients and healthy individuals (n = 1,070) and attempted to validate a selection of signals through genotyping in an independent cohort (n = 5,138). We identified three new MS-associated variants at C-X-C motif chemokine receptor 5 (CXCR5), Ts translation elongation factor, mitochondrial (TSFM) and cytochrome P450 family 24 subfamily A member 1 (CYP24A1). Rs10892307 resulted in a new signal at the CXCR5 region that explains one of the associations with MS within the locus. This polymorphism and three others in high linkage disequilibrium mapped within regulatory regions. Of them, rs11602393 showed allele-dependent enhancer activity in the forward orientation as determined by luciferase reporter assays. Immunophenotyping using peripheral blood mononuclear cells from MS patients associated the minor allele of rs10892307 with increased percentage of regulatory T cells expressing CXCR5. This work reports a new signal for the CXCR5 MS risk locus and points to rs11602393 as the causal variant. The expansion of CXCR5+ circulating regulatory T cells induced by this variant could cause its MS association.

Targeted resequencing reveals rare variants enrichment in multiple sclerosis susceptibility genes.

Although genome-wide association studies have identified a number of common variants associated with multiple sclerosis (MS) susceptibility, little is known about the relevance of rare variants. Here, we aimed to explore the role of rare variants in 14 MS risk genes (FCRL1, RGS1, TIMMDC1, HHEX, CXCR5, LTBR, TSFM, GALC, TRAF3, STAT3, TNFSF14, IFI30, CD40, and CYP24A1) by targeted resequencing in an Iberian population of 524 MS cases and 546 healthy controls. Four rare variants-enriched regions within CYP24A1, FCRL1, RGS1, and TRAF3 were identified as significantly associated with MS. Functional studies revealed significantly decreased regulator of G protein signaling 1 (RGS1) gene expression levels in peripheral blood mononuclear cells from MS patients with RGS1 rare variants compared to noncarriers, whereas no significant differences in gene expression were observed for CYP24A1, FCRL1, and TRAF3 between rare variants carriers and noncarriers. Immunophenotyping showed significant decrease in RGS1 expression in peripheral blood B lymphocytes from MS patients with RGS1 rare variants relative to noncarriers. Lastly, peripheral blood mononuclear cell from MS patients carrying RGS1 rare variants showed significantly lower induction of RGS1 gene expression by interferon-ß compared to MS patients lacking RGS1 variants. The presence of rare variants in RGS1 reinforce the ideas of high genetic heterogeneity and a role of rare variants in MS pathogenesis.

Genome-phenome explorer (GePhEx): a tool for the visualization and interpretation of phenotypic relationships supported by genetic evidence.

MOTIVATION: Association studies based on SNP arrays and Next Generation Sequencing technologies have enabled the discovery of thousands of genetic loci related to human diseases. Nevertheless, their biological interpretation is still elusive, and their medical applications limited. Recently, various tools have been developed to help bridging the gap between genomes and phenomes. To our knowledge, however none of these tools allows users to retrieve the phenotype-wide list of genetic variants that may be linked to a given disease or to visually explore the joint genetic architecture of different pathologies. RESULTS: We present the Genome-Phenome Explorer (GePhEx), a web-tool easing the visual exploration of phenotypic relationships supported by genetic evidences. GePhEx is primarily based on the thorough analysis of linkage disequilibrium between disease-associated variants and also considers relationships based on genes, pathways or drug-targets, leveraging on publicly available variant-disease associations to detect potential relationships between diseases. We demonstrate that GePhEx does retrieve well-known relationships as well as novel ones, and that, thus, it might help shedding light on the patho-physiological mechanisms underlying complex diseases. To this end, we investigate the potential relationship between schizophrenia and lung cancer, first detected using GePhEx and provide further evidence supporting a functional link between them. AVAILABILITY AND IMPLEMENTATION: GePhEx is available at: https://gephex.ega-archive.org/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A pharmacogenetic study implicates NINJ2 in the response to Interferon-ß in multiple sclerosis.

BACKGROUND: Multiple sclerosis (MS) is a disease in which biomarker identification is fundamental to predict response to treatments and to deliver the optimal drug to patients. We previously found an association between rs7298096, a polymorphism upstream to the NINJ2 gene, and the 4-year response to interferon-ß (IFNß) treatment in MS patients. OBJECTIVES: To analyse the association between rs7298096 and time to first relapse (TTFR) during IFNß therapy in MS patients and to better investigate its functional role. METHODS: Survival analysis was applied in three MS cohorts from different countries (n = 1004). We also studied the role of the polymorphism on gene expression using GTEx portal and a luciferase assay. We interrogated GEO datasets to explore the relationship between NINJ2 expression, IFNß and TTFR. RESULTS: Rs7298096AA patients show a shorter TTFR than rs7298096G-carriers (Pmeta-analysis = 3 × 10-4, hazard ratio = 1.41). Moreover, rs7298096AA is associated with a higher NINJ2 expression in blood (p = 7.0 × 10-6), which was confirmed in vitro (p = 0.009). Finally, NINJ2 expression is downregulated by IFNß treatment and related to TTFR. CONCLUSIONS: Rs7298096 could influence MS disease activity during IFNß treatment by modulating NINJ2 expression in blood. The gene encodes for an adhesion molecule involved in inflammation and endothelial cells activation, supporting its role in MS.

NLRP3 polymorphisms and response to interferon-beta in multiple sclerosis patients.

We aimed to investigate whether NLR family, pyrin domain containing 3 (NLRP3) polymorphisms are associated with the response to interferon-beta (IFNß) in multiple sclerosis (MS) patients. A total of 14 NLRP3 polymorphisms were genotyped in a cohort of 665 relapsing-remitting MS patients recruited across 5 centers and classified into responders and non-responders according to clinical-radiological criteria after 1 year of IFNß treatment. A meta-analysis failed to demonstrate significant associations between the response to IFNß and NLRP3 polymorphisms. These findings do not support a role of polymorphisms located in the NLRP3 gene and the response to IFNß in MS patients.

Analysis of known amyotrophic lateral sclerosis and frontotemporal dementia genes reveals a substantial genetic burden in patients manifesting both diseases not carrying the C9orf72 expansion mutation.

Amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) are part of a clinical, pathological and genetic continuum. OBJECTIVES: The purpose of the present study was to assess the mutation burden that is present in patients with concurrent ALS and FTD (ALS/FTD) not carrying the chromosome 9 open reading frame 72 (C9orf72) hexanucleotide repeat expansion, the most important genetic cause in both diseases. METHODS: From an initial group of 973 patients with ALS, we retrospectively selected those patients fulfilling diagnostic criteria of concomitant ALS and FTD lacking the repeat expansion mutation in C9orf72. Our final study group consisted of 54 patients clinically diagnosed with ALS/FTD (16 with available postmortem neuropathological diagnosis). Data from whole exome sequencing were used to screen for mutations in known ALS and/or FTD genes. RESULTS: We identified 11 patients carrying a probable pathogenic mutation, representing an overall mutation frequency of 20.4%. TBK1 was the most important genetic cause of ALS/FTD (n=5; 9.3%). The second most common mutated gene was SQSTM1, with three mutation carriers (one of them also harboured a TBK1 mutation). We also detected probable pathogenic genetic alterations in TAF15, VCP and TARDBP and possible pathogenic mutations in FIG4 and ERBB4. CONCLUSION: Our results indicate a high genetic burden underlying the co-occurrence of ALS and FTD and expand the phenotype associated with TAF15, FIG4 and ERBB4 to FTD. A systematic screening of ALS and FTD genes could be indicated in patients manifesting both diseases without the C9orf72 expansion mutation, regardless of family history of disease.

Antagonistic pleiotropy and mutation accumulation influence human senescence and disease.

Senescence has long been a public health challenge as well as a fascinating evolutionary problem. There is neither a universally accepted theory for its ultimate causes, nor a consensus about what may be its impact on human health. Here we test the predictions of two evolutionary explanations of senescence-mutation accumulation and antagonistic pleiotropy-which postulate that genetic variants with harmful effects in old ages can be tolerated, or even favoured, by natural selection at early ages. Using data from genome-wide association studies (GWAS), we study the effects of genetic variants associated with diseases appearing at different periods in life, when they are expected to have different impacts on fitness. Data fit theoretical expectations. Namely, we observe higher risk allele frequencies combined with large effect sizes for late-onset diseases, and detect a significant excess of early-late antagonistically pleiotropic variants that, strikingly, tend to be harboured by genes related to ageing. Beyond providing systematic, genome-wide evidence for evolutionary theories of senescence in our species and contributing to the long-standing question of whether senescence is the result of adaptation, our approach reveals relationships between previously unrelated pathologies, potentially contributing to tackling the problem of an ageing population.

Detection of genomic rearrangements from targeted resequencing data in Parkinson's disease patients.

BACKGROUND: The analysis of coverage depth in next-generation sequencing data allows the detection of gene dose alterations. We explore the frequency of such structural events in a Spanish cohort of sporadic PD cases. METHODS: Gene dose alterations were detected with the eXome-Hidden Markov Model (XHMM) software from depth of coverage in resequencing data available for 38 Mendelian and other risk PD loci in 394 individuals (249 cases and 145 controls) and subsequently validated by quantitative PCR. RESULTS: We identified 10 PD patients with exon dosage alterations in PARK2, GBA-GBAP1, and DJ1. Additional functional variants, including 2 novel nonsense mutations (p.Arg1552Ter in LRRK2 and p.Trp90Ter in PINK1), were confirmed by Sanger sequencing. This combined approach disclosed the genetic cause of 12 PD cases. CONCLUSIONS: Gene dose alterations related to PD can be correctly identified from targeting resequencing data. This approach substantially improves the detection rate of cases with causal genetic alterations. © 2016 The Authors. Movement Disorders published by Wiley Periodicals, Inc. on behalf of International Parkinson and Movement Disorder Society.

Properties of human disease genes and the role of genes linked to Mendelian disorders in complex disease aetiology.

Do genes presenting variation that has been linked to human disease have different biological properties than genes that have never been related to disease? What is the relationship between disease and fitness? Are the evolutionary pressures that affect genes linked to Mendelian diseases the same to those acting on genes whose variation contributes to complex disorders? The answers to these questions could shed light on the architecture of human genetic disorders and may have relevant implications when designing mapping strategies in future genetic studies. Here we show that, relative to non-disease genes, human disease (HD) genes have specific evolutionary profiles and protein network properties. Additionally, our results indicate that the mutation-selection balance renders an insufficient account of the evolutionary history of some HD genes and that adaptive selection could also contribute to shape their genetic architecture. Notably, several biological features of HD genes depend on the type of pathology (complex or Mendelian) with which they are related. For example, genes harbouring both causal variants for Mendelian disorders and risk factors for complex disease traits (Complex-Mendelian genes), tend to present higher functional relevance in the protein network and higher expression levels than genes associated only with complex disorders. Moreover, risk variants in Complex-Mendelian genes tend to present higher odds ratios than those on genes associated with the same complex disorders but with no link to Mendelian diseases. Taken together, our results suggest that genetic variation at genes linked to Mendelian disorders plays an important role in driving susceptibility to complex disease.

Signatures of Evolutionary Adaptation in Quantitative Trait Loci Influencing Trace Element Homeostasis in Liver.

Essential trace elements possess vital functions at molecular, cellular, and physiological levels in health and disease, and they are tightly regulated in the human body. In order to assess variability and potential adaptive evolution of trace element homeostasis, we quantified 18 trace elements in 150 liver samples, together with the expression levels of 90 genes and abundances of 40 proteins involved in their homeostasis. Additionally, we genotyped 169 single nucleotide polymorphism (SNPs) in the same sample set. We detected significant associations for 8 protein quantitative trait loci (pQTL), 10 expression quantitative trait loci (eQTLs), and 15 micronutrient quantitative trait loci (nutriQTL). Six of these exceeded the false discovery rate cutoff and were related to essential trace elements: 1) one pQTL for GPX2 (rs10133290); 2) two previously described eQTLs for HFE (rs12346) and SELO (rs4838862) expression; and 3) three nutriQTLs: The pathogenic C282Y mutation at HFE affecting iron (rs1800562), and two SNPs within several clustered metallothionein genes determining selenium concentration (rs1811322 and rs904773). Within the complete set of significant QTLs (which involved 30 SNPs and 20 gene regions), we identified 12 SNPs with extreme patterns of population differentiation (FST values in the top 5% percentile in at least one HapMap population pair) and significant evidence for selective sweeps involving QTLs at GPX1, SELENBP1, GPX3, SLC30A9, and SLC39A8. Overall, this detailed study of various molecular phenotypes illustrates the role of regulatory variants in explaining differences in trace element homeostasis among populations and in the human adaptive response to environmental pressures related to micronutrients.

Mendelian genes for Parkinson's disease contribute to the sporadic forms of the disease.

Parkinson's disease (PD) can be divided into familial (Mendelian) and sporadic forms. A number of causal genes have been discovered for the Mendelian form, which constitutes 10-20% of the total cases. Genome-wide association studies have successfully uncovered a number of susceptibility loci for sporadic cases but those only explain a small fraction (6-7%) of PD heritability. It has been observed that some genes that confer susceptibility to PD through common risk variants also contain rare causing mutations for the Mendelian forms of the disease. These results suggest a possible functional link between Mendelian and sporadic PD and led us to investigate the role that rare and low-frequency variants could have on the sporadic form. Through a targeting approach, we have resequenced at 49× coverage the exons and regulatory regions of 38 genes (including Mendelian and susceptibility PD genes) in 249 sporadic PD patients and 145 unrelated controls of European origin. Unlike susceptibility genes, Mendelian genes show a clear general enrichment of rare functional variants in PD cases, observed directly as well as with Tajima's D statistic and several collapsing methods. Our findings suggest that rare variation on PD Mendelian genes may have a role in the sporadic forms of the disease.