Genome-wide association studies (GWAS) have been highly informative in discovering disease-associated loci but are not designed to capture all structural variations in the human genome. Using long-read sequencing data, we discovered widespread structural variation within SINE-VNTR-Alu (SVA) elements, a class of great ape-specific transposable elements with gene-regulatory roles, which represents a major source of structural variability in the human population. We highlight the presence of structurally variable SVAs (SV-SVAs) in neurological disease-associated loci, and we further associate SV-SVAs to disease-associated SNPs and differential gene expression using luciferase assays and expression quantitative trait loci data. Finally, we genetically deleted SV-SVAs in the BIN1 and CD2AP Alzheimer's disease-associated risk loci and in the BCKDK Parkinson's disease-associated risk locus and assessed multiple aspects of their gene-regulatory influence in a human neuronal context. Together, this study reveals a novel layer of genetic variation in transposable elements that may contribute to identification of the structural variants that are the actual drivers of disease associations of GWAS loci.

@en

Studying the genome of centenarians may give insights into the molecular mechanisms underlying extreme human longevity and the escape of age-related diseases. Here, we set out to construct polygenic risk scores (PRSs) for longevity and to investigate the functions of longevity-associated variants. Using a cohort of centenarians with maintained cognitive health (N = 343), a population-matched cohort of older adults from 5 cohorts (N = 2905), and summary statistics data from genome-wide association studies on parental longevity, we constructed a PRS including 330 variants that significantly discriminated between centenarians and older adults. This PRS was also associated with longer survival in an independent sample of younger individuals (p =. 02), leading up to a 4-year difference in survival based on common genetic factors only. We show that this PRS was, in part, able to compensate for the deleterious effect of the APOE-ϵ4 allele. Using an integrative framework, we annotated the 330 variants included in this PRS by the genes they associate with. We find that they are enriched with genes associated with cellular differentiation, developmental processes, and cellular response to stress. Together, our results indicate that an extended human life span is, in part, the result of a constellation of variants each exerting small advantageous effects on aging-related biological mechanisms that maintain overall health and decrease the risk of age-related diseases.

@en

Human longevity is influenced by the genetic risk of age-related diseases. As Alzheimer’s disease (AD) represents a common condition at old age, an interplay between genetic factors affecting AD and longevity is expected. We explored this interplay by studying the prevalence of AD-associated single-nucleotide-polymorphisms (SNPs) in cognitively healthy centenarians, and replicated findings in a parental-longevity GWAS. We found that 28/38 SNPs that increased AD-risk also associated with lower odds of longevity. For each SNP, we express the imbalance between AD- and longevity-risk as an effect-size distribution. Based on these distributions, we grouped the SNPs in three groups: 17 SNPs increased AD-risk more than they decreased longevity-risk, and were enriched for β-amyloid metabolism and immune signaling; 11 variants reported a larger longevity-effect compared to their AD-effect, were enriched for endocytosis/immune-signaling, and were previously associated with other age-related diseases. Unexpectedly, 10 variants associated with an increased risk of AD and higher odds of longevity. Altogether, we show that different AD-associated SNPs have different effects on longevity, including SNPs that may confer general neuro-protective functions against AD and other age-related diseases.

@en

BACKGROUND: Several collaborative genome-wide-association studies (GWAS) have characterized the genetic landscape of Alzheimer's disease (AD), which now counts >70 single-nucleotide polymorphisms (SNPs) associated with AD-risk. METHOD: We linked these SNPs to their affected biological pathways, and combined the effect of multiple SNPs into pathway-specific polygenic-risk-scores (PRS). Using a genetically homogeneous dataset of 2,458 AD-patients (age=70.2±10.4), 3,848 healthy controls (age=61.1±14.8), 343 cognitively healthy centenarians (age=101.1±1.8) and 705 individuals with mild-cognitive-impairment (MCI, age=69.1±8.9), we studied the association between pathway-PRS and AD-risk, MCI-to-AD conversion, and resilience against AD. RESULT: With an integrative strategy, we linked 72 AD-associated SNPs to 5 clusters of biological pathways: beta-amyloid cluster, immune cluster, vascular cluster, endocytosis cluster and a trafficking cluster (Figure 1). A PRS comprising all 72 SNPs (excluding APOE-SNPs) was significantly associated with AD-risk (OR=1.53, p<2e-16), progression of MCI to AD dementia (OR=1.45, p=5.5e-4), and, albeit in the opposite direction, resilience against AD in healthy centenarians (OR=0.82, p=6.4e-4) (Figure 2, Figure 3). At the pathway level, all pathway-PRS significantly associated with increased AD-risk while specifically the vascular-, endocytosis- and trafficking-PRS associated with resilience against AD (p<0.05) (Figure 4). Interestingly, the beta-amyloid-, endocytosis and trafficking-PRSs were significantly associated with MCI-to-AD progression (p<0.05), with a higher-PRS leading up to 1.5-year shorter conversion time. CONCLUSION: Our results add on the efficacy of pathway-specific PRSs to identify individuals at highest or the lowest risk for the development of AD. Further, our results suggest that specific biological pathways are more important at different stages of the disease or in fact, associate with the escape of disease. AD diagnostic procedures may benefit from the identification of individual, pathway-specific vulnerabilities. In the future, pathway-specific PRSs may contribute to the selection of patients who may benefit most from pathway-specific treatment strategies.@en

Dysfunction of the endolysosomal-autophagy network is emerging as an important pathogenic process in Alzheimer's disease. Mutations in the sorting receptor-encoding gene SORL1 cause autosomal-dominant Alzheimer's disease, and SORL1 variants increase risk for late-onset AD. To understand the contribution of SORL1 mutations to AD pathogenesis, we analyze the effects of a SORL1 truncating mutation on SORL1 protein levels and endolysosome function in human neurons. We find that truncating mutation results in SORL1 haploinsufficiency and enlarged endosomes in human neurons. Analysis of isogenic SORL1 wild-type, heterozygous, and homozygous null neurons demonstrates that, whereas SORL1 haploinsufficiency results in endosome dysfunction, complete loss of SORL1 leads to additional defects in lysosome function and autophagy. Neuronal endolysosomal dysfunction caused by loss of SORL1 is relieved by extracellular antisense oligonucleotide-mediated reduction of APP protein, demonstrating that PSEN1, APP, and SORL1 act in a common pathway regulating the endolysosome system, which becomes dysfunctional in AD.

@en

Genetic association studies are frequently used to study the genetic basis of numerous human phenotypes. However, the rapid interrogation of how well a certain genomic region associates across traits as well as the interpretation of genetic associations is often complex and requires the integration of multiple sources of annotation, which involves advanced bioinformatic skills. We developed snpXplorer, an easy-to-use web-server application for exploring Single Nucleotide Polymorphisms (SNP) association statistics and to functionally annotate sets of SNPs. snpXplorer can superimpose association statistics from multiple studies, and displays regional information including SNP associations, structural variations, recombination rates, eQTL, linkage disequilibrium patterns, genes and gene-expressions per tissue. By overlaying multiple GWAS studies, snpXplorer can be used to compare levels of association across different traits, which may help the interpretation of variant consequences. Given a list of SNPs, snpXplorer can also be used to perform variant-to-gene mapping and gene-set enrichment analysis to identify molecular pathways that are overrepresented in the list of input SNPs. snpXplorer is freely available at https://snpxplorer.net. Source code, documentation, example files and tutorial videos are available within the Help section of snpXplorer and at https://github.com/TesiNicco/snpXplorer.

@en

Genetic discoveries of Alzheimer’s disease are the drivers of our understanding, and together with polygenetic risk stratification can contribute towards planning of feasible and efficient preventive and curative clinical trials. We first perform a large genetic association study by merging all available case-control datasets and by-proxy study results (discovery n = 409,435 and validation size n = 58,190). Here, we add six variants associated with Alzheimer’s disease risk (near APP, CHRNE, PRKD3/NDUFAF7, PLCG2 and two exonic variants in the SHARPIN gene). Assessment of the polygenic risk score and stratifying by APOE reveal a 4 to 5.5 years difference in median age at onset of Alzheimer’s disease patients in APOE ɛ4 carriers. Because of this study, the underlying mechanisms of APP can be studied to refine the amyloid cascade and the polygenic risk score provides a tool to select individuals at high risk of Alzheimer’s disease.

@en

BACKGROUND: Dementia in families can be caused by one genetic variant. Identifying these so-called monogenic causes of dementia is important, because it explains the origin of dementia in families and raises the possibility of predictive testing for relatives. Still, we do not know how frequent these monogenic causes are, due to strict selection criteria for DNA testing based on age at onset and a positive family history. We aimed to identify the mutation frequency of monogenetic dementia by screening a large, unselected cohort of patients from a memory clinic specialized in early onset dementia. METHOD: Between 1-1-2010 and 1-7-2012, 1,138 patients visited the Alzheimer Center Amsterdam and all patients underwent the same diagnostic trajectory (van der Flier et al 2018). Of these, 1,093 (96%) consented to research and donated blood. An additional 73 patients (7%) were excluded because suboptimal DNA quality. Whole exome sequencing and C9orf72-repeat length analysis was performed in the remaining 1,020 patients (90% of all patients). All variants in 54 dementia related genes were analysed and classified. C9orf72 repeat length >30 repeat units was considered pathogenic. Variants classified as likely pathogenic and pathogenic were considered a monogenic cause of dementia. RESULT: The average age at presentation was 62(±6SD) years and 419(41%) were female. We found a monogenic cause of dementia in 34 (3.4%) patients. The most frequently occurring (likely) pathogenic variants were found in the genes C9orf72, PSEN1, NOTCH3 and MAPT. These genes explained 65% of the (likely) pathogenic variants. The clinical characteristics of the patients screened positive; 50% were diagnosed with dementia and 50% had a first degree relative with dementia; and 40% had an age at presentation of <60 years and 44% was between 60 to 70 years, the remainder was older. Surprisingly, only 47% of the patients with a (likely) pathogenic variant fulfilled our current criteria (based on diagnosis, age and familial history) for offering diagnostic genetic testing (21% of patients). CONCLUSION: In our large unselected cohort of patients from a memory clinic specialized in early onset dementia, the majority of the patients with a monogenetic predisposition for dementia did not fulfill our current criteria for genetic testing.@en

The practical use of graph-based reference genomes depends on the ability to align reads to them. Performing substring queries to paths through these graphs lies at the core of this task. The combination of increasing pattern length and encoded variations inevitably leads to a combinatorial explosion of the search space. Instead of heuristic filtering or pruning steps to reduce the complexity, we propose CHOP, a method that constrains the search space by exploiting haplotype information, bounding the search space to the number of haplotypes so that a combinatorial explosion is prevented. We show that CHOP can be applied to large and complex datasets, by applying it on a graph-based representation of the human genome encoding all 80 million variants reported by the 1000 Genomes Project.

@en

The human genome harbors numerous structural variants (SVs) which, due to their repetitive nature, are currently underexplored in short-read whole-genome sequencing approaches. Using single-molecule, real-time (SMRT) long-read sequencing technology in combination with FALCON-Unzip, we generated a de novo assembly of the diploid genome of a 115-year-old Dutch cognitively healthy woman. We combined this assembly with two previously published haploid assemblies (CHM1 and CHM13) and the GRCh38 reference genome to create a compendium of SVs that occur across five independent human haplotypes using the graph-based multi-genome aligner REVEAL. Across these five haplotypes, we detected 31,680 euchromatic SVs (>50 bp). Of these, ~62% were comprised of repetitive sequences with ‘variable number tandem repeats’ (VNTRs), ~10% were mobile elements (Alu, L1, and SVA), while the remaining variants were inversions and indels. We observed that VNTRs with GC-content >60% and repeat patterns longer than 15 bp were 21-fold enriched in the subtelomeric regions (within 5 Mb of the ends of chromosome arms). VNTR lengths can expand to exceed a critical length which is associated with impaired gene transcription. The genes that contained most VNTRs, of which PTPRN2 and DLGAP2 are the most prominent examples, were found to be predominantly expressed in the brain and associated with a wide variety of neurological disorders. Repeat-induced variation represents a sizeable fraction of the genetic variation in human genomes and should be included in investigations of genetic factors associated with phenotypic traits, specifically those associated with neurological disorders. We make available the long and short-read sequence data of the supercentenarian genome, and a compendium of SVs as identified across 5 human haplotypes.

@en

Objectives: The objectives of this study were to investigate the effect of genetic and social factors on depressive symptoms and depression over time and to test whether social factors moderate the relationship between depressive symptoms and its underlying genetics in later life. Methods: The study included 2,279 participants with a mean follow-up of 15 years from the Longitudinal Aging Study Amsterdam with genotyping data. The personal genetic loading for depression was estimated for each participant by calculating a polygenic risk scores (PRS-D), based on 23,032 single nucleotide polymorphisms associated with major depression in a large genome-wide association study. Partner status, network size, received and given emotional support were assessed via questionnaires and depressive symptoms were assessed using the CES-D Scale. A CES-D Scale of 16 and higher was considered as clinically relevant depression. Results: Higher PRS-D was associated with more depressive symptoms whereas having a partner and having a larger network size were independently associated with less depressive symptoms. After extra adjustment for education, cognitive function and functional limitations, giving more emotional support was also associated with less depressive symptoms. No evidence for gene-environment interaction between PRS-D and social factors was found. Similar results were found for clinically relevant depression. Conclusion: Genetic and social factors are independently associated with depressive symptoms over time in older adults. Strategies that boost social functioning should be encouraged in the general population of older adults regardless of the genetic liability for depression.

@en

Developing Alzheimer’s disease (AD) is influenced by multiple genetic variants that are involved in five major AD-pathways. Per individual, these pathways may differentially contribute to the modification of the AD-risk. The pathways involved in the resilience against AD have thus far been poorly addressed. Here, we investigated to what extent each molecular mechanism associates with (i) the increased risk of AD and (ii) the resilience against AD until extreme old age, by comparing pathway-specific polygenic risk scores (pathway-PRS). We used 29 genetic variants associated with AD to develop pathway-PRS for five major pathways involved in AD. We developed an integrative framework that allows multiple genes to associate with a variant, and multiple pathways to associate with a gene. We studied pathway-PRS in the Amsterdam Dementia Cohort of well-phenotyped AD patients (N = 1895), Dutch population controls from the Longitudinal Aging Study Amsterdam (N = 1654) and our unique 100-plus Study cohort of cognitively healthy centenarians who avoided AD (N = 293). Last, we estimated the contribution of each pathway to the genetic risk of AD in the general population. All pathway-PRS significantly associated with increased AD-risk and (in the opposite direction) with resilience against AD (except for angiogenesis, p < 0.05). The pathway that contributed most to the overall modulation of AD-risk was β-amyloid metabolism (29.6%), which was driven mainly by APOE-variants. After excluding APOE variants, all pathway-PRS associated with increased AD-risk (except for angiogenesis, p < 0.05), while specifically immune response (p = 0.003) and endocytosis (p = 0.0003) associated with resilience against AD. Indeed, the variants in these latter two pathways became the main contributors to the overall modulation of genetic risk of AD (45.5% and 19.2%, respectively). The genetic variants associated with the resilience against AD indicate which pathways are involved with maintained cognitive functioning until extreme ages. Our work suggests that a favorable immune response and a maintained endocytosis pathway might be involved in general neuro-protection, which highlight the need to investigate these pathways, next to β-amyloid metabolism.

@en

To determine the small vessel disease spectrum associated with cysteine-altering NOTCH3 variants in community-dwelling individuals by analyzing the clinical and neuroimaging features of UK Biobank participants harboring such variants. The exome and genome sequencing datasets of the UK Biobank (n = 50,000) and cohorts of cognitively healthy elderly (n = 751) were queried for cysteine-altering NOTCH3 variants. Brain MRIs of individuals harboring such variants were scored according to Standards for Reporting Vascular Changes on Neuroimaging criteria, and clinical information was extracted with ICD-10 codes. Clinical and neuroimaging data were compared to age- and sex-matched UK Biobank controls and clinically diagnosed patients from the Dutch cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) registry. We identified 108 individuals harboring a cysteine-altering NOTCH3 variant (2.2 of 1,000), of whom 75% have a variant that has previously been reported in CADASIL pedigrees. Almost all variants were located in 1 of the NOTCH3 protein epidermal growth factor-like repeat domains 7 to 34. White matter hyperintensity lesion load was higher in individuals with NOTCH3 variants than in controls (p = 0.006) but lower than in patients with CADASIL with the same variants (p < 0.001). Almost half of the 24 individuals with brain MRI had a Fazekas score of 0 or 1 up to age 70 years. There was no increased risk of stroke. Although community-dwelling individuals harboring a cysteine-altering NOTCH3 variant have a higher small vessel disease MRI burden than controls, almost half have no MRI abnormalities up to age 70 years. This shows that NOTCH3 cysteine altering variants are associated with an extremely broad phenotypic spectrum, ranging from CADASIL to nonpenetrance.

@en

Human longevity is heritable, but genome-wide association (GWA) studies have had limitedsuccess. Here, we perform two meta-analyses of GWA studies of a rigorous longevityphenotype definition including 11,262/3484 cases surviving at or beyond the age corre-sponding to the 90th/99th survival percentile, respectively, and 25,483 controls whose ageat death or at last contact was at or below the age corresponding to the 60th survivalpercentile. Consistent with previous reports, rs429358 (apolipoprotein E (ApoE)ε4) isassociated with lower odds of surviving to the 90th and 99th percentile age, while rs7412(ApoEε2) shows the opposite. Moreover, rs7676745, located nearGPR78, associates withlower odds of surviving to the 90th percentile age. Gene-level association analysis reveals arole for tissue-specific expression of multiple genes in longevity. Finally, genetic correlation ofthe longevity GWA results with that of several disease-related phenotypes points to a sharedgenetic architecture between health and longevity@en

The detection of genetic loci associated with Alzheimer’s disease (AD) requires large numbers of cases and controls because variant effect sizes are mostly small. We hypothesized that variant effect sizes should increase when individuals who represent the extreme ends of a disease spectrum are considered, as their genomes are assumed to be maximally enriched or depleted with disease-associated genetic variants. We used 1,073 extensively phenotyped AD cases with relatively young age at onset as extreme cases (66.3 ± 7.9 years), 1,664 age-matched controls (66.0 ± 6.5 years) and 255 cognitively healthy centenarians as extreme controls (101.4 ± 1.3 years). We estimated the effect size of 29 variants that were previously associated with AD in genome-wide association studies. Comparing extreme AD cases with centenarian controls increased the variant effect size relative to published effect sizes by on average 1.90-fold (SE = 0.29, p = 9.0 × 10⁻⁴). The effect size increase was largest for the rare high-impact TREM2 (R74H) variant (6.5-fold), and significant for variants in/near ECHDC3 (4.6-fold), SLC24A4-RIN3 (4.5-fold), NME8 (3.8-fold), PLCG2 (3.3-fold), APOE-ε2 (2.2-fold), and APOE-ε4 (twofold). Comparing extreme phenotypes enabled us to replicate the AD association for 10 variants (p < 0.05) in relatively small samples. The increase in effect sizes depended mainly on using centenarians as extreme controls: the average variant effect size was not increased in a comparison of extreme AD cases and age-matched controls (0.94-fold, p = 6.8 × 10⁻¹), suggesting that on average the tested genetic variants did not explain the extremity of the AD cases. Concluding, using centenarians as extreme controls in AD case–control studies boosts the variant effect size by on average twofold, allowing the replication of disease-association in relatively small samples.

@en

Next-generation sequencing has contributed to our understanding of the genetics of Alzheimer's disease (AD) and has explained a substantial part of the missing heritability of familial AD. We sequenced 19 exomes from 8 Dutch families with a high AD burden and identified EIF2AK3, encoding for protein kinase RNA-like endoplasmic reticulum kinase (PERK), as a candidate gene. Gene-based burden analysis in a Dutch AD exome cohort containing 547 cases and 1070 controls showed a significant association of EIF2AK3 with AD (OR 1.84 [95% CI 1.07–3.17], p-value 0.03), mainly driven by the variant p.R240H. Genotyping of this variant in an additional cohort from the Rotterdam Study showed a trend toward association with AD (p-value 0.1). Immunohistochemical staining with pPERK and peIF2α of 3 EIF2AK3 AD carriers showed an increase in hippocampal neuronal cells expressing these proteins compared with nondemented controls, but no difference was observed in AD noncarriers. This study suggests that rare variants in EIF2AK3 may be associated with disease risk in AD.

@en

Accumulating evidence suggests that genetic variants in the SORL1 gene are associated with Alzheimer disease (AD), but a strategy to identify which variants are pathogenic is lacking. In a discovery sample of 115 SORL1 variants detected in 1908 Dutch AD cases and controls, we identified the variant characteristics associated with SORL1 variant pathogenicity. Findings were replicated in an independent sample of 103 SORL1 variants detected in 3193 AD cases and controls. In a combined sample of the discovery and replication samples, comprising 181 unique SORL1 variants, we developed a strategy to classify SORL1 variants into five subtypes ranging from pathogenic to benign. We tested this pathogenicity screen in SORL1 variants reported in two independent published studies. SORL1 variant pathogenicity is defined by the Combined Annotation Dependent Depletion (CADD) score and the minor allele frequency (MAF) reported by the Exome Aggregation Consortium (ExAC) database. Variants predicted strongly damaging (CADD score >30), which are extremely rare (ExAC-MAF <1 × 10 '5) increased AD risk by 12-fold (95% CI 4.2-34.3; P=5 × 10 '9). Protein-truncating SORL1 mutations were all unknown to ExAC and occurred exclusively in AD cases. More common SORL1 variants (ExAC-MAF≥1 × 10 '5) were not associated with increased AD risk, even when predicted strongly damaging. Findings were independent of gender and the APOE-I 4 allele. High-risk SORL1 variants were observed in a substantial proportion of the AD cases analyzed (2%). Based on their effect size, we propose to consider high-risk SORL1 variants next to variants in APOE, PSEN1, PSEN2 and APP for personalized risk assessments in clinical practice.

@en