Copy Number Variation and Haplotype Analysis of 17q21.31 Reveals Increased Risk Associated with Progressive Supranuclear Palsy and Gene Expression Changes in Neuronal Cells
AbstractBackgroundProgressive supranuclear palsy (PSP) is a neurodegenerative disorder involving pathological deposition of tau that includes glial inclusions and specific regional vulnerability patterns. Therapeutic developments are hampered by incomplete understanding of disease mechanisms. Few studies have examined its cell type‐specific effects. We performed RNA‐sequencing in a large cohort of PSP patients and controls to understand molecular pathways disrupted across prefrontal cortex and caudate nucleus in PSP. We also performed multi‐omic single‐nucleus sequencing to deconvolute cell type‐specific changes.MethodRNA‐sequencing was conducted in 250 PSP patients and age‐ and sex‐matched controls, with a subset (N = 23) undergoing multi‐omic single‐nucleus RNA‐sequencing and ATAC‐sequencing (10x Genomics). Sequencing data was aligned to GRCh38 and quantified using STAR and RSEM (RNA‐seq) and CellRanger (single‐nucleus) and corrected for technical and biological covariates. Differential expression was calculated with Dream (RNA‐seq) and Dreamlet (single‐nucleus sequencing), and co‐expression network analysis (WGCNA) used to extract disease‐relevant gene modules.ResultRNA‐sequencing analysis identified ∼4000 differentially expressed genes in both brain regions. Despite differences in cellular composition and degree of disease impact, both regions share patterns of dysregulation related to protein folding and processing, indicating common pathological mechanisms. Differences between the regions includes disruption of cellular homeostasis and intracellular trafficking in caudate, while prefrontal cortex shows suppression of metabolic processes and immune response. Preliminary findings from the single cell genomic analyses suggest the greatest changes occur in neurons of both regions, where there is a shared pattern of protein processing dysregulation and a distinct enrichment in metabolic processing within prefrontal cortex interneurons.ConclusionWe identified thousands of dysregulated genes in two impacted brain regions in PSP and show that this dysregulation is reflected in specific neuronal populations in these regions. Identifying regional and specific cell type‐specific mechanisms related to PSP will further our understanding of causative disease mechanisms and provide new avenues for therapeutic developments.
Rare Genetic Risk in Progressive Supranuclear Palsy
AbstractBackgroundGenetic risk in common variants have been identified for Progressive Supranuclear Palsy (PSP), but common variants only account for a portion of disease heritability. A proportion of the remaining disease heritability is likely contributed by rare variants identifiable via sequencing. Leveraging whole‐genome sequencing of 1719 PSP cases and 2940 controls from the Alzheimer’s Disease Sequencing Project, we detected single nucleotide variants (SNVs), insertion/deletions (INDELs), and copy number variants (CNVs). Note that all samples are non‐Hispanic whites with identity by descent (IBD) >0.25, and no outlier is >6 standard deviation in the principal component analysis. We conducted association analysis for variants, genes and gene sets to identify rare genetic risk for PSP.MethodAfter quality control (variants and genotypes notated PASS by the calling algorithms), we detected 85,228,234 rare SNVs/INDELs and 163,987 rare CNVs (minor allele frequency (MAF) < 1%). For rare SNVs/INDELs, the Sequence Kernel Association Test Optimized was performed to identify rare genetic PSP risk from protein truncating variants (PTV) and damaging missense variants (MAF<0.01%) in both protein coding genes and gene sets that were previously associated with PSP based on protein co‐expression networks. For rare CNVs, a permutation‐based association test implemented in Plink was performed to identify CNVs associated with PSP.ResultCombining protein truncating variants (PTV) and damaging missense variants, we identified one genome‐wide significant gene, FBXO38 (FDR = 0.0003, OR = 1.78). We also found enrichment of PTVs and damaging missense variants in genes comprising the C2 module, a neuronal module previously identified in post mortem PSP brain (FDR = 0.0016). For CNV analysis, there were 17 significant deletions and 2 significant duplications that were associated with PSP (Padjust < 0.05). The deletion (18:58488230‐58488489) in ALPK2, which has been identified as an Alzheimer’s Disease risk gene, was the only CNV that conferred a higher‐risk for PSP. Particularly, one deletion (17:46009357‐46009595) was in linkage disequilibrium with MAPT H2 haplotype which has been associated various neurodegenerative diseases including PSP.ConclusionWe identified several new genes and a protein co‐expression network in which risk variants were associated with PSP. These genes and pathways provide potential candidate drug targets for future investigation.