Genome-wide determinants of mortality and motor progression in Parkinson's disease.

There are 90 independent genome-wide significant genetic risk variants for Parkinson's disease (PD) but currently only five nominated loci for PD progression. The biology of PD progression is likely to be of central importance in defining mechanisms that can be used to develop new treatments. We studied 6766 PD patients, over 15,340 visits with a mean follow-up of between 4.2 and 15.7 years and carried out genome-wide survival studies for time to a motor progression endpoint, defined by reaching Hoehn and Yahr stage 3 or greater, and death (mortality). There was a robust effect of the APOE ε4 allele on mortality in PD. We also identified a locus within the TBXAS1 gene encoding thromboxane A synthase 1 associated with mortality in PD. We also report 4 independent loci associated with motor progression in or near MORN1, ASNS, PDE5A, and XPO1. Only the non-Gaucher disease causing GBA1 PD risk variant E326K, of the known PD risk variants, was associated with mortality in PD. Further work is needed to understand the links between these genomic variants and the underlying disease biology. However, these may represent new candidates for disease modification in PD.

The annotation of GBA1 has been concealed by its protein-coding pseudogene GBAP1.

Mutations in GBA1 cause Gaucher disease and are the most important genetic risk factor for Parkinson's disease. However, analysis of transcription at this locus is complicated by its highly homologous pseudogene, GBAP1. We show that >50% of short RNA-sequencing reads mapping to GBA1 also map to GBAP1. Thus, we used long-read RNA sequencing in the human brain, which allowed us to accurately quantify expression from both GBA1 and GBAP1. We discovered significant differences in expression compared to short-read data and identify currently unannotated transcripts of both GBA1 and GBAP1. These included protein-coding transcripts from both genes that were translated in human brain, but without the known lysosomal function-yet accounting for almost a third of transcription. Analyzing brain-specific cell types using long-read and single-nucleus RNA sequencing revealed region-specific variations in transcript expression. Overall, these findings suggest nonlysosomal roles for GBA1 and GBAP1 with implications for our understanding of the role of GBA1 in health and disease.

Genome sequence analyses identify novel risk loci for multiple system atrophy.

Multiple system atrophy (MSA) is an adult-onset, sporadic synucleinopathy characterized by parkinsonism, cerebellar ataxia, and dysautonomia. The genetic architecture of MSA is poorly understood, and treatments are limited to supportive measures. Here, we performed a comprehensive analysis of whole genome sequence data from 888 European-ancestry MSA cases and 7,128 controls to systematically investigate the genetic underpinnings of this understudied neurodegenerative disease. We identified four significantly associated risk loci using a genome-wide association study approach. Transcriptome-wide association analyses prioritized USP38-DT, KCTD7, and lnc-KCTD7-2 as novel susceptibility genes for MSA within these loci, and single-nucleus RNA sequence analysis found that the associated variants acted as cis-expression quantitative trait loci for multiple genes across neuronal and glial cell types. In conclusion, this study highlights the role of genetic determinants in the pathogenesis of MSA, and the publicly available data from this study represent a valuable resource for investigating synucleinopathies.

Prenatally-diagnosed renal failure: an ethical framework for decision-making.

The Children's Hospital Working Group has developed an ethical framework to guide patient care and research for prenatally diagnosed severe renal anomalies. It identifies ethical challenges in communication, timing of decisions and scarce resources. Key elements include shared decision-making, establishing a trusting relationship, and managing disagreement. The ethical framework will be used to develop a clinical pathway that operationalizes the key values of trust, honesty, transparency, beneficence, nonmaleficence, respecting parental authority, professional integrity, and justice.

The Effect of Steroids on Prenatally Diagnosed Lung Lesions.

BACKGROUND: Open fetal resection for large lung lesions has virtually been replaced by maternal steroid administration. Despite this paradigm shift, little is known about the effects steroids have on lung lesion growth in utero. METHODS: A 10-year retrospective review of all prenatally diagnosed lung lesions cared for at our fetal care center was performed. We evaluated the effects of prenatal steroids on congenital pulmonary airway malformation (CPAM)-volume-ratio (CVR), distinguishing change in CVR among CPAMs, bronchopulmonary sequestrations (BPS), and bronchial atresias. We also correlated fetal ultrasound and MRI findings with pathology to determine the accuracy of prenatal diagnosis. RESULTS: We evaluated 199 fetuses with a prenatal lung lesion. Fifty-four (27 %) were treated with prenatal steroids with a subsequent 21 % mean reduction in the CVR (2.1 ± 1.4 to 1.1 ± 0.4, p = 0.003). Fetuses with hydrops and mediastinal shift who were treated with steroids rarely had resolution of these radiographic findings. Postnatal pathology was available for 91/199 patients (45.7 %). The most common diagnosis was CPAM (42/91, 46 %), followed by BPS (30/91, 33 %), and bronchial atresia (14/91, 15 %). Fetuses who received steroids and had pathology consistent with CPAM were more likely to have a reduction in their CVR (p = 0.02). Fetal ultrasound correctly diagnosed the type of lung lesion in 75 % of cases and fetal MRI in 81 % of cases. CONCLUSIONS: Prenatally diagnosed CPAMs are more likely to respond to maternal steroids than BPS or bronchial atresias. Knowing the diagnosis in utero could aid to steward steroid usage, however, fetal imagining modalities are not perfect in distinguishing subtype. LEVEL OF EVIDENCE: III.

ensemblQueryR: fast, flexible and high-throughput querying of Ensembl LD API endpoints in R.

We present ensemblQueryR, an R package for querying Ensembl linkage disequilibrium (LD) endpoints. This package is flexible, fast and user-friendly, and optimised for high-throughput querying. ensemblQueryR uses functions that are intuitive and amenable to custom code integration, familiar R object types as inputs and outputs as well as providing parallelisation functionality. For each Ensembl LD endpoint, ensemblQueryR provides two functions, permitting both single- and multi-query modes of operation. The multi-query functions are optimised for large query sizes and provide optional parallelisation to leverage available computational resources and minimise processing time. We demonstrate improved computational performance of ensemblQueryR over an exisiting tool in terms of random access memory (RAM) usage and speed, delivering a 10-fold speed increase whilst using a third of the RAM. Finally, ensemblQueryR is near-agnostic to operating system and computational architecture through Docker and singularity images, making this tool widely accessible to the scientific community.

aws-s3-integrity-check: an open-source bash tool to verify the integrity of a dataset stored on Amazon S3.

Amazon Simple Storage Service (Amazon S3) is a widely used platform for storing large biomedical datasets. Unintended data alterations can occur during data writing and transmission, altering the original content and generating unexpected results. However, no open-source and easy-to-use tool exists to verify end-to-end data integrity. Here, we present aws-s3-integrity-check, a user-friendly, lightweight, and reliable bash tool to verify the integrity of a dataset stored in an Amazon S3 bucket. Using this tool, we only needed ∼114 min to verify the integrity of 1,045 records ranging between 5 bytes and 10 gigabytes and occupying ∼935 gigabytes of the Amazon S3 cloud. Our aws-s3-integrity-check tool also provides file-by-file on-screen and log-file-based information about the status of each integrity check. To our knowledge, this tool is the only open-source one that allows verifying the integrity of a dataset uploaded to the Amazon S3 Storage quickly, reliably, and efficiently. The tool is freely available for download and use at https://github.com/SoniaRuiz/aws-s3-integrity-check and https://hub.docker.com/r/soniaruiz/aws-s3-integrity-check.

Analysis of subcellular RNA fractions demonstrates significant genetic regulation of gene expression in human brain post-transcriptionally.

Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

Genome-wide Analysis of Motor Progression in Parkinson Disease.

Background and Objectives: The genetic basis of Parkinson disease (PD) motor progression is largely unknown. Previous studies of the genetics of PD progression have included small cohorts and shown a limited overlap with genetic PD risk factors from case-control studies. Here, we have studied genomic variation associated with PD motor severity and early-stage progression in large longitudinal cohorts to help to define the biology of PD progression and potential new drug targets. Methods: We performed a GWAS meta-analysis of early PD motor severity and progression up to 3 years from study entry. We used linear mixed-effect models with additive effects, corrected for age at diagnosis, sex, and the first 5 genetic principal components to assess variability in axial, limb, and total Movement Disorder Society-Unified Parkinson's Disease Rating Scale (MDS-UPDRS) III scores. Results: We included 3,572 unrelated European ancestry patients with PD from 5 observational cohorts and 1 drug trial. The average AAO was 62.6 years (SD = 9.83), and 63% of participants were male. We found an average increase in the total MDS-UPDRS III score of 2.3 points/year. We identified an association between PD axial motor progression and variation at the GJA5 locus at 1q12 (ß = -0.25, SE = 0.04, p = 3.4e-10). Exploration of the regulation of gene expression in the region (cis-expression quantitative trait loci [eQTL] analysis) showed that the lead variant was associated with expression of ACP6, a lysophosphatidic acid phosphatase that regulates mitochondrial lipid biosynthesis (cis-eQTL p-values in blood and brain RNA expression data sets: <10-14 in eQTLGen and 10-7 in PsychEncode). Discussion: Our study highlights the potential role of mitochondrial lipid homeostasis in the progression of PD, which may be important in establishing new drug targets that might modify disease progression.

Regional genetic correlations highlight relationships between neurodegenerative disease loci and the immune system.

Neurodegenerative diseases, including Alzheimer's and Parkinson's disease, are devastating complex diseases resulting in physical and psychological burdens on patients and their families. There have been important efforts to understand their genetic basis leading to the identification of disease risk-associated loci involved in several molecular mechanisms, including immune-related pathways. Regional, in contrast to genome-wide, genetic correlations between pairs of immune and neurodegenerative traits have not been comprehensively explored, but could uncover additional immune-mediated risk-associated loci. Here, we systematically assess the role of the immune system in five neurodegenerative diseases by estimating regional genetic correlations between these diseases and immune-cell-derived single-cell expression quantitative trait loci (sc-eQTLs). We also investigate correlations between diseases and protein levels. We observe significant (FDR < 0.01) correlations between sc-eQTLs and neurodegenerative diseases across 151 unique genes, spanning both the innate and adaptive immune systems, across most diseases tested. With Parkinson's, for instance, RAB7L1 in CD4+ naïve T cells is positively correlated and KANSL1-AS1 is negatively correlated across all adaptive immune cell types. Follow-up colocalization highlight candidate causal risk genes. The outcomes of this study will improve our understanding of the immune component of neurodegeneration, which can warrant repurposing of existing immunotherapies to slow disease progression.

The non-specific lethal complex regulates genes and pathways genetically linked to Parkinson's disease.

Genetic variants conferring risks for Parkinson's disease have been highlighted through genome-wide association studies, yet exploration of their specific disease mechanisms is lacking. Two Parkinson's disease candidate genes, KAT8 and KANSL1, identified through genome-wide studies and a PINK1-mitophagy screen, encode part of the histone acetylating non-specific lethal complex. This complex localizes to the nucleus, where it plays a role in transcriptional activation, and to mitochondria, where it has been suggested to have a role in mitochondrial transcription. In this study, we sought to identify whether the non-specific lethal complex has potential regulatory relationships with other genes associated with Parkinson's disease in human brain. Correlation in the expression of non-specific lethal genes and Parkinson's disease-associated genes was investigated in primary gene co-expression networks using publicly-available transcriptomic data from multiple brain regions (provided by the Genotype-Tissue Expression Consortium and UK Brain Expression Consortium), whilst secondary networks were used to examine cell type specificity. Reverse engineering of gene regulatory networks generated regulons of the complex, which were tested for heritability using stratified linkage disequilibrium score regression. Prioritized gene targets were then validated in vitro using a QuantiGene multiplex assay and publicly-available chromatin immunoprecipitation-sequencing data. Significant clustering of non-specific lethal genes was revealed alongside Parkinson's disease-associated genes in frontal cortex primary co-expression modules, amongst other brain regions. Both primary and secondary co-expression modules containing these genes were enriched for mainly neuronal cell types. Regulons of the complex contained Parkinson's disease-associated genes and were enriched for biological pathways genetically linked to disease. When examined in a neuroblastoma cell line, 41% of prioritized gene targets showed significant changes in mRNA expression following KANSL1 or KAT8 perturbation. KANSL1 and H4K8 chromatin immunoprecipitation-sequencing data demonstrated non-specific lethal complex activity at many of these genes. In conclusion, genes encoding the non-specific lethal complex are highly correlated with and regulate genes associated with Parkinson's disease. Overall, these findings reveal a potentially wider role for this protein complex in regulating genes and pathways implicated in Parkinson's disease.

Local genetic correlations exist among neurodegenerative and neuropsychiatric diseases.

Genetic correlation ([Formula: see text]) between traits can offer valuable insight into underlying shared biological mechanisms. Neurodegenerative diseases overlap neuropathologically and often manifest comorbid neuropsychiatric symptoms. However, global [Formula: see text] analyses show minimal [Formula: see text] among neurodegenerative and neuropsychiatric diseases. Importantly, local [Formula: see text] s can exist in the absence of global relationships. To investigate this possibility, we applied LAVA, a tool for local [Formula: see text] analysis, to genome-wide association studies of 3 neurodegenerative diseases (Alzheimer's disease, Lewy body dementia and Parkinson's disease) and 3 neuropsychiatric disorders (bipolar disorder, major depressive disorder and schizophrenia). We identified several local [Formula: see text] s missed in global analyses, including between (i) all 3 neurodegenerative diseases and schizophrenia and (ii) Alzheimer's and Parkinson's disease. For those local [Formula: see text] s identified in genomic regions containing disease-implicated genes, such as SNCA, CLU and APOE, incorporation of expression quantitative trait loci identified genes that may drive genetic overlaps between diseases. Collectively, we demonstrate that complex genetic relationships exist among neurodegenerative and neuropsychiatric diseases, highlighting putative pleiotropic genomic regions and genes. These findings imply sharing of pathogenic processes and the potential existence of common therapeutic targets.

The IPDGC/GP2 Hackathon - an open science event for training in data science, genomics, and collaboration using Parkinson's disease data.

Open science and collaboration are necessary to facilitate the advancement of Parkinson's disease (PD) research. Hackathons are collaborative events that bring together people with different skill sets and backgrounds to generate resources and creative solutions to problems. These events can be used as training and networking opportunities, thus we coordinated a virtual 3-day hackathon event, during which 49 early-career scientists from 12 countries built tools and pipelines with a focus on PD. Resources were created with the goal of helping scientists accelerate their own research by having access to the necessary code and tools. Each team was allocated one of nine different projects, each with a different goal. These included developing post-genome-wide association studies (GWAS) analysis pipelines, downstream analysis of genetic variation pipelines, and various visualization tools. Hackathons are a valuable approach to inspire creative thinking, supplement training in data science, and foster collaborative scientific relationships, which are foundational practices for early-career researchers. The resources generated can be used to accelerate research on the genetics of PD.

The contribution of Neanderthal introgression and natural selection to neurodegenerative diseases.

Humans are thought to be more susceptible to neurodegeneration than equivalently-aged primates. It is not known whether this vulnerability is specific to anatomically-modern humans or shared with other hominids. The contribution of introgressed Neanderthal DNA to neurodegenerative disorders remains uncertain. It is also unclear how common variants associated with neurodegenerative disease risk are maintained by natural selection in the population despite their deleterious effects. In this study, we aimed to quantify the genome-wide contribution of Neanderthal introgression and positive selection to the heritability of complex neurodegenerative disorders to address these questions. We used stratified-linkage disequilibrium score regression to investigate the relationship between five SNP-based signatures of natural selection, reflecting different timepoints of evolution, and genome-wide associated variants of the three most prevalent neurodegenerative disorders: Alzheimer's disease, amyotrophic lateral sclerosis and Parkinson's disease. We found no evidence for enrichment of positively-selected SNPs in the heritability of Alzheimer's disease, amyotrophic lateral sclerosis and Parkinson's disease, suggesting that common deleterious disease variants are unlikely to be maintained by positive selection. There was no enrichment of Neanderthal introgression in the SNP-heritability of these disorders, suggesting that Neanderthal admixture is unlikely to have contributed to disease risk. These findings provide insight into the origins of neurodegenerative disorders within the evolution of Homo sapiens and addresses a long-standing debate, showing that Neanderthal admixture is unlikely to have contributed to common genetic risk of neurodegeneration in anatomically-modern humans.

Piloting a multidisciplinary approach to improve outcomes of fetal whole exome sequencing: An overview of workflow and case example.

INTRODUCTION: Whole exome sequencing (WES) has increasingly become integrated into prenatal care and genetic testing pathways. Current studies of prenatal WES have focused on diagnostic yield. The possibility of obtaining a variant of uncertain significance and lack of provider expertise are frequently described as common barriers to clinical integration of prenatal WES. We describe the implementation and workflow for a multidisciplinary approach to effectively integrate prenatal WES into maternal-fetal care to overcome these barriers. METHODS: A multidisciplinary team reviews and approves potential cases for WES. This team reviews WES results, reclassifying variants as appropriate and provides recommendations for postnatal care. A detailed description of this workflow is provided, and a case example is included to demonstrate effectiveness of this approach. Our team has approved 62 cases for WES with 45 patients ultimately pursuing WES. We have achieved a diagnostic yield of 40% and the multidisciplinary team has played a role in variant interpretation in 50% of the reported variants of uncertain significance. CONCLUSIONS: This approach facilitates communication between prenatal and postnatal care teams and provides accurate interpretation and recommendations for identified fetal variants. This model can be replicated to ensure appropriate patient care and effective integration of novel genomic technologies into prenatal settings.

Functional genomics provide key insights to improve the diagnostic yield of hereditary ataxia.

Improvements in functional genomic annotation have led to a critical mass of neurogenetic discoveries. This is exemplified in hereditary ataxia, a heterogeneous group of disorders characterised by incoordination from cerebellar dysfunction. Associated pathogenic variants in more than 300 genes have been described, leading to a detailed genetic classification partitioned by age-of-onset. Despite these advances, up to 75% of patients with ataxia remain molecularly undiagnosed even following whole genome sequencing, as exemplified in the 100 000 Genomes Project. This study aimed to understand whether we can improve our knowledge of the genetic architecture of hereditary ataxia by leveraging functional genomic annotations, and as a result, generate insights and strategies that raise the diagnostic yield. To achieve these aims, we used publicly-available multi-omics data to generate 294 genic features, capturing information relating to a gene's structure, genetic variation, tissue-specific, cell-type-specific and temporal expression, as well as protein products of a gene. We studied these features across genes typically causing childhood-onset, adult-onset or both types of disease first individually, then collectively. This led to the generation of testable hypotheses which we investigated using whole genome sequencing data from up to 2182 individuals presenting with ataxia and 6658 non-neurological probands recruited in the 100 000 Genomes Project. Using this approach, we demonstrated a high short tandem repeat (STR) density within childhood-onset genes suggesting that we may be missing pathogenic repeat expansions within this cohort. This was verified in both childhood- and adult-onset ataxia patients from the 100 000 Genomes Project who were unexpectedly found to have a trend for higher repeat sizes even at naturally-occurring STRs within known ataxia genes, implying a role for STRs in pathogenesis. Using unsupervised analysis, we found significant similarities in genomic annotation across the gene panels, which suggested adult- and childhood-onset patients should be screened using a common diagnostic gene set. We tested this within the 100 000 Genomes Project by assessing the burden of pathogenic variants among childhood-onset genes in adult-onset patients and vice versa. This demonstrated a significantly higher burden of rare, potentially pathogenic variants in conventional childhood-onset genes among individuals with adult-onset ataxia. Our analysis has implications for the current clinical practice in genetic testing for hereditary ataxia. We suggest that the diagnostic rate for hereditary ataxia could be increased by removing the age-of-onset partition, and through a modified screening for repeat expansions in naturally-occurring STRs within known ataxia-associated genes, in effect treating these regions as candidate pathogenic loci.

Association between the LRP1B and APOE loci and the development of Parkinson's disease dementia.

Parkinson's disease is one of the most common age-related neurodegenerative disorders. Although predominantly a motor disorder, cognitive impairment and dementia are important features of Parkinson's disease, particularly in the later stages of the disease. However, the rate of cognitive decline varies among Parkinson's disease patients, and the genetic basis for this heterogeneity is incompletely understood. To explore the genetic factors associated with rate of progression to Parkinson's disease dementia, we performed a genome-wide survival meta-analysis of 3923 clinically diagnosed Parkinson's disease cases of European ancestry from four longitudinal cohorts. In total, 6.7% of individuals with Parkinson's disease developed dementia during study follow-up, on average 4.4 ± 2.4 years from disease diagnosis. We have identified the APOE ε4 allele as a major risk factor for the conversion to Parkinson's disease dementia [hazard ratio = 2.41 (1.94-3.00), P = 2.32 × 10-15], as well as a new locus within the ApoE and APP receptor LRP1B gene [hazard ratio = 3.23 (2.17-4.81), P = 7.07 × 10-09]. In a candidate gene analysis, GBA variants were also identified to be associated with higher risk of progression to dementia [hazard ratio = 2.02 (1.21-3.32), P = 0.007]. CSF biomarker analysis also implicated the amyloid pathway in Parkinson's disease dementia, with significantly reduced levels of amyloid ß42 (P = 0.0012) in Parkinson's disease dementia compared to Parkinson's disease without dementia. These results identify a new candidate gene associated with faster conversion to dementia in Parkinson's disease and suggest that amyloid-targeting therapy may have a role in preventing Parkinson's disease dementia.

IntroVerse: a comprehensive database of introns across human tissues.

Dysregulation of RNA splicing contributes to both rare and complex diseases. RNA-sequencing data from human tissues has shown that this process can be inaccurate, resulting in the presence of novel introns detected at low frequency across samples and within an individual. To enable the full spectrum of intron use to be explored, we have developed IntroVerse, which offers an extensive catalogue on the splicing of 332,571 annotated introns and a linked set of 4,679,474 novel junctions covering 32,669 different genes. This dataset has been generated through the analysis of 17,510 human control RNA samples from 54 tissues provided by the Genotype-Tissue Expression Consortium. IntroVerse has two unique features: (i) it provides a complete catalogue of novel junctions and (ii) each novel junction has been assigned to a specific annotated intron. This unique, hierarchical structure offers multiple uses, including the identification of novel transcripts from known genes and their tissue-specific usage, and the assessment of background splicing noise for introns thought to be mis-spliced in disease states. IntroVerse provides a user-friendly web interface and is freely available at https://rytenlab.com/browser/app/introverse.