HH

H. Holstege

info

Please Note

41 records found

Journal article (2026) - Matthijs W.J. de Waal, Sven J. van der Lee, Melanie Lunding, Lynn Boonkamp, Nolan Barrett, Jan Raska, Niccolo Tesi, Marc Hulsman, Henne Holstege, More authors...
INTRODUCTION: The sortilin-related receptor (SORL1) directs APP and Aβ trafficking within the retromer pathway. Cleavage at the cell surface releases soluble SORL1 (sSORL1) into cerebrospinal fluid (CSF). We examined whether CSF-sSORL1 can serve as an in vivo marker of genetically impaired SORL1. METHODS: CSF-sSORL1 was quantified by enzyme-linked immunosorbent assay (ELISA) in 218 participants: 90 carriers of SORL1 variants, 78 SORL1-wildtype (WT) AD patients, and 50 SORL1-WT controls. RESULTS: sSORL1 concentrations were significantly lower in carriers of protein-truncating and damaging missense variants. In SORL1-WT patients, CSF-sSORL1 correlated with pTau181 but not with Aβ42 among AD patients, and did not differ between patients and controls. DISCUSSION: These findings suggest that impaired SORL1 trafficking reduces receptor delivery to the cell surface and thereby decreases sSORL1 shedding, supporting its potential use as a pathway-specific biomarker. Highlights: Enzyme-linked immunosorbent assay (ELISA) enables quantitative measurement of soluble sortilin-related receptor (sSORL1) in cerebrospinal fluid (CSF). sSORL1 levels are reduced in CSF from carriers of a pathogenic SORL1 variant. CSF-sSORL1 levels correlate with tau pathology in Alzheimer's disease. sSORL1 levels represent an in vivo biomarker of SORL1 function. ...
Journal article (2025) - Marc Hulsman, Itziar de Rojas, Catarina B. Ferreira, Sven van der Lee, Iris Jansen, Gennady Roshchupkin, More Authors..., Sven van der Lee, Henne Holstege, Jeroen van Rooij, Jasper Van Dongen, Niccolo Tesí, Marcel J.T. Reinders, Marc Hulsman, Georgios Hadjigeorgiou
Traditional statistical approaches have advanced our understanding of the genetics of complex diseases, yet are limited to linear additive models. Here we applied machine learning (ML) to genome-wide data from 41,686 individuals in the largest European consortium on Alzheimer’s disease (AD) to investigate the effectiveness of various ML algorithms in replicating known findings, discovering novel loci, and predicting individuals at risk. We utilised Gradient Boosting Machines (GBMs), biological pathway-informed Neural Networks (NNs), and Model-based Multifactor Dimensionality Reduction (MB-MDR) models. ML approaches successfully captured all genome-wide significant genetic variants identified in the training set and 22% of associations from larger meta-analyses. They highlight 6 novel loci which replicate in an external dataset, including variants which map to ARHGAP25, LY6H, COG7, SOD1 and ZNF597. They further identify novel association in AP4E1, refining the genetic landscape of the known SPPL2A locus. Our results demonstrate that machine learning methods can achieve predictive performance comparable to classical approaches in genetic epidemiology and have the potential to uncover novel loci that remain undetected by traditional GWAS. These insights provide a complementary avenue for advancing the understanding of AD genetics. ...
Journal article (2025) - Harold Bae, Zeyuan Song, Amanat Ali, Niccolò Tesi, Marc Hulsman, Sven van der Lee, Natasja M. van Schoor, Marcel Reinders, Henne Holstege, More Authors...
We constructed a polygenic protective score specific to Alzheimer’s disease (AD PPS) based on the current literature among the participants enrolled in five studies of healthy aging and extreme longevity in the USA, Europe, and Asia. This AD PPS did not include variants on apolipoprotein E (APOE) gene. Comparisons of AD PPS in different data sets of healthy agers and centenarians showed that centenarians have stronger genetic protection against AD compared to individuals without familial longevity. The current study also shows evidence that this genetic protection increases with increasingly older ages in centenarians (centenarians who died before reaching age 105 years, semi-supercentenarians who reached age 105 to 109 years, and supercentenarians who reached age 110 years and older). However, the genetic protection was of modest size: the average increase in AD PPS was approximately one additional protective allele per 5 years of gained lifetime. Additionally, we show that the higher AD PPS was associated with better cognitive function and decreased mortality. Taken together, this analysis suggests that individuals who achieve the most extreme ages, on average, have the greatest protection against AD. This finding is robust to different genetic backgrounds with important implications for universal applicability of therapeutics that target this AD PPS. ...
Journal article (2025) - Sven J. Van Der Lee, Marc Hulsman, Rosalina Van Spaendonk, Jetske Van Der Schaar, Janna Dijkstra, Niccoló Tesi, Marcel Reinders, Itziar De Rojas, Henne Holstege, More Authors...
Background and ObjectivesIdentifying genetic causes of dementia in patients visiting memory clinics is important for patient care and family planning. Traditional clinical selection criteria for genetic testing may miss carriers of pathogenic variants in dementia-related genes. This study aimed identify how many carriers we are missing and to optimize criteria for selecting patients for genetic counseling in memory clinics.MethodsIn this clinical cohort study, we retrospectively genetically tested patients during 2.5 years (2010-2012) visiting the Alzheimer Center Amsterdam, a specialized memory clinic. Genetic tests consisted of a 54-gene dementia panel, focusing on Class IV/V variants per American College of Medical Genetics and Genomics guidelines, including APP duplications and the C9ORF72 repeat expansion. We determined the prevalence of pathogenic variants and propose new eligibility criteria for genetic testing in memory clinics. The eligibility criteria were prospectively applied for 1 year (2021-2022), and results were compared with the retrospective cohort.ResultsGenetic tests were retrospectively performed in in 1,022 of 1,138 patients (90%) who consecutively visited the memory clinic. Among these, 1,022 patients analyzed (mean age 62.1 ± 8.9 years; 40.4% were female), 34 pathogenic variant carriers were identified (3.3%), with 24 being symptomatic. Previous clinical criteria would have identified only 15 carriers (44% of all carriers, 65% of symptomatic carriers). The proposed criteria increased identification to 22 carriers (62.5% of all carriers, 91% of symptomatic carriers). In the prospective cohort, 148 (28.7%) of 515 patients were eligible for testing under the new criteria. Of the 90 eligible patients who consented to testing, 13 pathogenic carriers were identified, representing a 73% increase compared with the previous criteria.DiscussionWe found that patients who visit a memory clinic and carry a pathogenic genetic variant are often not eligible for genetic testing. The proposed new criteria improve the identification of patients with a genetic cause for their cognitive complaints. In systems without practical or financial barriers to genetic testing, the new criteria can enhance personalized care. In other countries where the health care systems differs and in other genetic ancestry groups, the performance of the criteria may be different. ...
Journal article (2025) - Henne Holstege, Alex N. Salazar, Lydian Knoop, Yolande A.L. Pijnenburg, Sven J. van der Lee, Sanduni Wijesekera, Jana Krizova, Mikko Hiltunen, Marcel JT Reinders, More Authors...
Background
Genome-wide association studies (GWAS) linked TMEM106B variants to susceptibility for neurodegenerative diseases, but the causal genetic elements remain unclear.

Method
We used genotyping data from 5,792 Alzheimer disease cases and controls, and applied COJO to identify haplotypes in the TMEM106B locus that independently associated with AD. Then, we used long-read sequencing data from 513 individuals to annotate these haplotypes with structural variations that map into them.

Results
Analysis of the genotyping data revealed that the TMEM106B locus consists of four major haplotypes: HA/Ha (covering the coding region), and HB/Hb (covering the upstream regulatory region). These combine into four combinations with varying population-frequencies: HAB (57%), HaB (34%), Hab (9%), and HAb (<1%). Long-read sequencing of 513 individuals showed that HA haplotypes (marked by 185-Threonine) carry unique methylated CpG sites and an AluYb8-retrotransposon in the 3' UTR, while the Ha haplotypes are marked by the 185-Serine allele. Hb haplotypes carry several structural variants (SVs) in nearby distal enhancers, including a 19 Kbp rearrangement, absent in all other haplotypes. Joint association models revealed that the HAB combination (AluYb8+185-Threonine) is risk-increasing, while Hab (SVs+185-Serine) confers the protective effect. HaB (185-Serine only) is neutral, while HAb was too rare to assess. Relative to middle-aged non-demented controls, cognitively healthy centenarians were more enriched with Hab (OR=1.49, padj=2.18×10-2) than with HaB (OR=1.23, padj=5.06×10-2). Proteomic analysis of temporal cortex tissues (n = 182) indicated that relative to the neutral HaB combination, the protective Hab is associated with 1.1-fold lower TMEM106B C-terminal peptide abundance, while the risk-increasing HAB is associated with 1.16-fold higher abundance.

Conclusion
Our data indicates that the genetic structure underlying the association of the TMEM106B locus with neurodegenerative diseases is driven by the effect of multiple haplotypes. ...
Journal article (2025) - Susan K. Rohde, Maruelle C. Luimes, Linda M.C. Lorenz, Patricia Fierro-Hernández, Annemieke J.M. Rozemuller, Marc Hulsman, Marieke J.I. Graat, Myke E. van der Hoorn, Henne Holstege, More authors...
IMPORTANCE Older individuals without dementia often have amyloid-beta (Aβ) Thal phases similar to patients with Alzheimer disease (AD), suggesting that Aβ pathology may be a benign consequence of aging. OBJECTIVE To explore whether Aβ pathology in centenarians is associated with cognitive performance. DESIGN, SETTING, AND PARTICIPANTS This longitudinal cohort study used cross-sectional data on antemortem cognitive performance and postmortem neuropathology of participants in the Dutch 100-plus Study. Cognitive performance was measured a median of 10 (IQR, 3-13) months before postmortem brain donation. From January 2013 to July 2022, 1187 centenarians who self-reported being cognitively healthy, confirmed by proxy, were approached: 406 were included and 95 donated their brain. Centenarians were compared with patients with clinicopathologically confirmed AD from the Netherlands Brain Bank. Data were analyzed from June 2022 to October 2024. MAIN OUTCOMES AND MEASURES Aβ pathology was assessed with the Thal phase for Aβ progression and by determining quantitative Aβ loads (percentage positive area) in the frontal, parietal, temporal, and occipital neocortices, 3 parahippocampal, and 5 hippocampal subregions. Aβ pathology was associated with performance on 13 neuropsychological tests assessing memory, fluency, attention/processing speed, and executive functioning, as well as 4 measures of global cognition. RESULTS This study evaluated Aβ pathology in 95 centenarians (median age at brain donation, 103.5 [IQR, 102.3-104.7] years; 71 female [75%] and 24 male [25%]) and 38 patients with AD (median age, 84 [IQR, 78-90] years; 18 female [47%] and 20 male [53%]). Global cognition parameters were available for all 95 centenarians and complete cognitive assessment for 72 centenarians (76%). A fraction of the centenarians had no Aβ load (9 of 95 [9%]), most had low Aβ load (53 of 95 [56%]) and, despite high Thal phases, about one-third (33 of 95 [35%]) had high Aβ load comparable with patients with AD. Centenarians with no or low Aβ load had significantly higher cognitive performance than centenarians with high Aβ loads. Higher Aβ loads across all 4 neocortical regions, cornu ammonis 3, cornu ammonis 1/subiculum, and the entorhinal cortex specifically affected executive functioning. Interestingly, 5 resilient centenarians maintained high cognitive performance despite having high Aβ loads; they had significantly less tau pathology compared with centenarians with high Aβ loads and low cognitive performance. CONCLUSIONS AND RELEVANCE These results indicate that Aβ pathology is not a benign consequence of aging. Even in the oldest individuals, Aβ and tau pathology interaction was consistent with the amyloid cascade hypothesis. ...
Journal article (2025) - Desiree D.S.H. de Bruin, Martin A. Haagmans, Kristiaan J. van der Gaag, Jerry Hoogenboom, Natalie E.C. Weiler, Niccoló Tesi, Henne Holstege, Marcel Reinders, Peter Henneman, More authors...
Introduction
The field of forensic DNA analysis has undergone rapid advancements in recent decades. The integration of massively parallel sequencing (MPS) has notably expanded the forensic toolkit, moving beyond identity matching to predicting phenotypic traits and biogeographical ancestry. This shift is of particular significance in cases where conventional DNA profiling fails to identify a single suspect. Supplementing forensic analyses with estimated biological age may be valuable but involves a complex and time-consuming DNA methylation analysis. This study explores and validates the performance of a comprehensive forensic third-generation sequencing assay utilizing Oxford Nanopore Technologies (ONT) in an adaptive and direct sequencing approach. We incorporated the most widely used forensic markers, i.e., STRs, SNPs, InDels, mitochondrial DNA (mtDNA), and two methylation-based clock classifiers, thereby combining forensic genetic and epigenetic analysis in one single workflow.

Methods and results
In our investigation, DNA from six anonymous individuals was sequenced using the ONT standard adaptive direct sequencing approach, reaching a mean percentage of on-target reads ranging from 6.6 % to 7.7 % per sample. ONT data was compared to standard MPS data and Illumina EPIC DNA methylation profiles. Basecalling employed recommended ONT software packages. TREAT was used for ONT-based analysis of autosomal and Y-chromosome STRs, achieving 90–92 % correct calls depending on allelic read depth thresholds. InDel analyses for two lower-quality samples proved challenging due to inadequate read depth, while the remaining four samples significantly contributed to the observed percentage markers (60.9 %) and correct calls (97.8 %). SNP analysis achieved a 98 % call rate, with only two mismatches and two missed alleles. ONT-generated DNA methylation data demonstrated Pearson’s correlation coefficients with EPIC data ranging from 0.67 to 0.97 for Horvath’s clock. Additional age-associated markers exhibited Pearson’s correlation coefficients with chronological age between 0.14 (ELOVL2) and 0.96 (FHL2) at read depths of <30 and <20, respectively. Despite excluding mtDNA from our targeted sequencing approach, adaptive proof-reading fragments covered the complete mtDNA with an average read depth of 21–72, showing 100 % concordance with reference data.

Discussion
Our exploratory study using ONT adaptive sequencing for conventional forensic and age associated DNA methylation markers showed high sequencing accuracy for a significant number of markers, showcasing ONT as a promising (epi)genetic forensic method. Future studies must address three critical aspects: determining clear quantity and quality measures and detection thresholds for accuracy, optimizing input DNA quantity for forensic casework expectations, and addressing ethical considerations associated with phenotype and ancestry analysis to prevent ethnic biases. ...
Journal article (2025) - Tanya N. Weerakkody, Hanna Sabelström, Shan V. Andrews, Jean Paul Chadarevian, Marcus Y. Chin, Niccolò Tesi, Henne Holstege, Mathew Blurton-Jones, Kathryn M. Monroe, More authors...
The Alzheimer’s disease (AD) genetic landscape identified microglia as a key disease-modifying cell type. Paired immunoglobulin-like type 2 receptor alpha (PILRA) is an immunoreceptor tyrosine–based inhibitory motif domain–containing inhibitory receptor, expressed by myeloid cells such as microglia. The known protective PILRA G78R gene variant reduces AD risk in apolipoprotein E4 (APOE4) carriers and is enriched in a cohort of healthy centenarians. However, mechanisms underlying protective effects in microglia are undefined. Here, we identified biological functions of PILRA in human induced pluripotent stem cell–derived microglia (iMG) and chimeric AD mice. PILRA knockout (KO) in iMG rescued ApoE4-mediated immunometabolic deficits and prevented lipotoxicity through increased lipid storage, improved mitochondrial bioenergetics, and antioxidant activity. PILRA KO also enhanced microglial chemotaxis and attenuated inflammation. With pharmacological inhibitor studies, we showed that peroxisome proliferator–activated receptor and signal transducer and activator of transcription 1/3 mediated PILRA-dependent microglial functions. AD mice transplanted with human PILRA KO microglia exhibited reduced amyloid pathology and rescued synaptic markers. A high-affinity ligand blocking PILRA antibody phenocopied PILRA KO iMG. These findings suggest that PILRA is a pharmacologically tractable therapeutic target for AD. ...
Journal article (2025) - Susan K. Rohde, Maruelle C. Luimes, Annemieke J.M. Rozemuller, Marieke J.I. Graat, Myke E. van der Hoorn, Dominique A.H. Daatselaar, Marc Hulsman, Sietske A.M. Sikkes, Henne Holstege, More authors...
Background
The hippocampus is differentially affected in Alzheimer's disease neuropathologic change (ADNC) versus primary age-related tauopathy (PART), an amyloid-beta (Aβ)-independent tauopathy: the CA2/CA1 hyperphosphorylated tau (pTau)-ratio is higher in PART, which inversely correlates with Aβ-burden. However, as the aging brain often presents mixed rather than uniform pathologies, we questioned whether these distinct hippocampal pTau distributions persist into extreme ages and how hippocampal Aβ- and pTau-distributions correlate with cognition in centenarians.

Method
We quantified Aβ- (6F/3D) and pTau (AT8)-burdens across eight hippocampal and parahippocampal subregions in 112 centenarians (median age 104, IQR 102-105), alongside 11 AD (median age 84, IQR 72-86) and 7 PART cases for comparison (median age 88, IQR 78-92; Figure 1). We compared CA2/CA1-pTau-ratio in centenarians who met PART criteria (Thal phase ≤2, Braak stage I-IV; n = 49) with centenarians who met ADNC criteria (intermediate/high according to NIA-AA guidelines; Thal phase ≥3, Braak stage III-VI; n = 50). Cognitive performance was assessed using 13 neuropsychological tests shortly before brain donation (median 10 months, IQR5-14, n = 72). Robust linear regression models were used to associate subregional Aβ- and pTau-burdens with cognitive performance, while adjusting for age, sex, and education.

Result
In line with previous findings, CA2/CA1-pTau-ratios were higher in younger PART cases compared to AD patients (median 3.0, IQR 2.1-3.6, min-max 1.6-4.2 vs. median 1.2, IQR 0.9-1.4, min-max 0.8-1.4; p <0.001). Surprisingly, CA2/CA1-pTau-ratios in centenarians with PART were comparable to centenarians with ADNC (median 1.3, IQR 1.1-2.0, min-max 0.3-10.8 vs. median 1.2, IQR 1.0-1.8, min-max 0.2-6.2; p = 0.684). Accordingly, CA2/CA1-pTau-ratio in centenarians was unrelated to Aβ-burden, Thal phase or Braak stage. Higher Aβ- and pTau-burdens associated with lower cognition, though through different subregions: cognition associated with Aβ-burden in the hippocampus (CA4, CA3, CA2, CA1/subiculum), whereas pTau-burden in the parahippocampus (presubiculum, entorhinal cortex, fusiform gyrus) associated with cognition.

Conclusion
In the oldest-old, PART and ADNC are less distinguishable by determinants observed in younger individuals: centenarians with ADNC may show age-related Aβ accumulation alongside PART-like pTau patterns, while centenarians meeting PART criteria do not always show PART-like pTau patterns. However, hippocampal Aβ-burden and parahippocampal pTau-burden associate with cognitive decline, highlighting subregional-specific vulnerability to pathology-driven cognitive decline. ...
Journal article (2025) - Olav M. Andersen, Matthijs W. J. de Waal, Giulia Monti, Niccolo Tesi, Anne Mette G. Jensen, Christa de Geus, Marcel J. T. Reinders, Marc Hulsman, Henne Holstege, More authors...
Background
Protein truncating variants (PTVs) in SORL1 are observed almost exclusively in Alzheimer’s Disease (AD) cases, but the effect of rare SORL1 missense variants is unclear.

Methods
To identify high-priority missense variants (HPVs), we applied ‘domain mapping of disease mutations’ for the 637 unique coding SORL1 variants detected in 18,959 AD-cases and 21,893 non-demented controls.

Results
In this sample, PTVs and HPVs associated with respectively a 35- and 10-fold increased risk of early onset AD and 17- and 6-fold increased risk of overall AD. The median age at onset (AAO) of PTV- and HPV-carriers was 62 and 64 years, and APOE-genotype contributed to AAO-variability. The median AAO of PTV- and HPV-carriers is ~8–10 years earlier than wild-type SORL1 carriers, matched for APOE-genotype. Specific HPVs are highly penetrant and lead to earlier AAOs than PTVs, suggesting possible dominant negative effects.

Conclusion
Our results justify a debate on whether HPV carriers should be considered for clinical counseling. ...
Journal article (2025) - Yaran Zhang, Marc Hulsman, Alex Salazar, Niccolò Tesi, Lydian Knoop, Sven van der Lee, Sanduni Wijesekera, Jana Krizova, Erik Jan Kamsteeg, Henne Holstege
Tandem repeats (TRs) occupy a significant portion of the human genome and are a source of polymorphisms due to variations in sizes and motif compositions. Some of these variations have been associated with various neuropathological disorders, highlighting the clinical importance of assessing the motif structure of TRs. Moreover, assessing the TR motif variation can offer valuable insights into evolutionary dynamics and population structure. Previously, characterizations of TRs were limited by short-read sequencing technology, which lacks the ability to accurately capture the full TR sequences. As long-read sequencing becomes more accessible and can capture the full complexity of TRs, there is now also a need for tools to characterize and analyze TRs using long-read data across multiple samples. In this study, we present MotifScope, a novel algorithm for the characterization and visualization of TRs based on a de novo k-mer approach for motif discovery. Comparative analysis against established tools reveals that MotifScope can identify a greater number of motifs and more accurately represent the underlying repeat sequences. Moreover, MotifScope has been specifically designed to enable motif composition comparisons across assemblies of different individuals, as well as across long-read sequencing reads within an individual, through combined motif discovery and sequence alignment. We showcase potential applications of MotifScope in diverse fields, including population genetics, clinical settings, and forensic analyses. ...
Journal article (2024) - Niccolo’ Tesi, Sven van der Lee, Marc Hulsman, Natasja M. van Schoor, Martijn Huisman, Yolande Pijnenburg, Wiesje M. van der Flier, Marcel Reinders, Henne Holstege
BACKGROUND: Alzheimer's disease (AD) prevalence increases with age, yet a small fraction of the population reaches ages > 100 years without cognitive decline. We studied the genetic factors associated with such resilience against AD. METHODS: Genome-wide association studies identified 86 single nucleotide polymorphisms (SNPs) associated with AD risk. We estimated SNP frequency in 2281 AD cases, 3165 age-matched controls, and 346 cognitively healthy centenarians. We calculated a polygenic risk score (PRS) for each individual and investigated the functional properties of SNPs enriched/depleted in centenarians. RESULTS: Cognitively healthy centenarians were enriched with the protective alleles of the SNPs associated with AD risk. The protective effect concentrated on the alleles in/near ANKH, GRN, TMEM106B, SORT1, PLCG2, RIN3, and APOE genes. This translated to >5-fold lower PRS in centenarians compared to AD cases (P = 7.69 × 10−71), and 2-fold lower compared to age-matched controls (P = 5.83 × 10−17). DISCUSSION: Maintaining cognitive health until extreme ages requires complex genetic protection against AD, which concentrates on the genes associated with the endolysosomal and immune systems. Highlights: Cognitively healthy cent enarians are enriched with the protective alleles of genetic variants associated with Alzheimer's disease (AD). The protective effect is concentrated on variants involved in the immune and endolysosomal systems. Combining variants into a polygenic risk score (PRS) translated to > 5-fold lower PRS in centenarians compared to AD cases, and ≈ 2-fold lower compared to middle-aged healthy controls. ...
Journal article (2024) - Niccoló Tesi, Alex Salazar, Yaran Zhang, Sven van der Lee, Marc Hulsman, Lydian Knoop, Sanduni Wijesekera, Marcel Reinders, Henne Holstege, More authors...
Tandem repeats (TRs) play important roles in genomic variation and disease risk in humans. Long-read sequencing allows for the accurate characterization of TRs; however, the underlying bioinformatics perspectives remain challenging. We present otter and TREAT: otter is a fast targeted local assembler, cross-compatible across different sequencing platforms. It is integrated in TREAT, an end-to-end workflow for TR characterization, visualization, and analysis across multiple genomes. In a comparison with existing tools based on long-read sequencing data from both Oxford Nanopore Technology (ONT, Simplex and Duplex) and Pacific Bioscience (PacBio, Sequel II and Revio), otter and TREAT achieve state-of-the-art genotyping and motif characterization accuracy. Applied to clinically relevant TRs, TREAT/otter significantly identify individuals with pathogenic TR expansions. When applied to a case-control setting, we replicate previously reported associations of TRs with Alzheimer's disease, including those near or within APOC1 (P = 2.63 × 10−9), SPI1 (P = 6.5 × 10−3), and ABCA7 (P = 0.04) genes. Finally, we use TREAT/otter to systematically evaluate potential biases when genotyping TRs using diverse ONT and PacBio long-read sequencing data sets. We show that, in rare cases (0.06%), long-read sequencing from coverage drops in TRs, including the disease-associated TRs in ABCA7 and RFC1 genes. Such coverage drops can lead to TR misgenotyping, hampering the accurate characterization of TR alleles. Taken together, our tools can accurately genotype TRs across different sequencing technologies and with minimal requirements, allowing end-to-end analysis and comparisons of TRs in human genomes, with broad applications in research and clinical fields. ...
Journal article (2024) - Floor C. Loonstra, Daniel Álvarez Sirvent, Niccoló Tesi, Henne Holstege, Eva M.M. Strijbis, Alex N. Salazar, Marc Hulsman, Sven J. Van Der Lee, Bernard Uitdehaag
Background and Objectives
More than 200 genetic variants have been associated with multiple sclerosis (MS) susceptibility. However, it is unclear to what extent genetic factors influence lifetime risk of MS. Using a population-based birth-year cohort, we investigate the effect of genetics on lifetime risk of MS.

Methods
In the Project Y study, we tracked down almost all persons with MS (pwMS) from birth year 1966 in the Netherlands. As control participants, we included non-MS participants from the Project Y cohort (born 1965–1967 in the Netherlands) and non-MS participants from the Amsterdam Dementia Cohort born between 1963 and 1969. Genetic variants associated with MS were determined in pwMS and control participants using genotyping or imputation methods. Polygenic risk scores (PRSs) based on variants and weights from the largest genetic study in MS were calculated for each participant and assigned into deciles based on the PRS distribution in the control participants. We examined the lifetime risk for each decile and the association between PRS and MS disease variables, including age at onset and time to secondary progression.

Results
MS-PRS was calculated for 285 pwMS (mean age 53.0 ± 0.9 years, 72.3% female) and 267 control participants (mean age 51.8 ± 3.2 years, 58.1% female). Based on the lifetime risk estimation, we observed that 1:2,739 of the women with the lowest 30% genetic risk developed MS, whereas 1:92 of the women with the top 10% highest risk developed MS. For men, only 1:7,900 developed MS in the lowest 30% genetic risk group, compared with 1:293 men with the top 10% genetic risk. The PRS was not significantly associated with age at onset and time to secondary progression in both sexes.

Discussion
Our results show that the lifetime risk of MS is strongly influenced by genetic factors. Our findings have the potential to support diagnostic certainty in individuals with suspected MS: a high PRS could strengthen a diagnosis, but especially a PRS from the lowest tail of the PRS distribution should be considered a red flag and could prevent misdiagnosing conditions that mimic MS. ...
Journal article (2023) - M. Zhang, G.A. Bouland, H. Holstege, M.J.T. Reinders
Background and Objectives With age, somatic mutations accumulated in human brain cells can lead to various neurologic disorders and brain tumors. Because the incidence rate of Alzheimer disease (AD) increases exponentially with age, investigating the association between AD and the accumulation of somatic mutation can help understand the etiology of AD. Methods We designed a somatic mutation detection workflow by contrasting genotypes derived from whole-genome sequencing (WGS) data with genotypes derived from scRNA-seq data and applied this workflow to 76 participants from the Religious Order Study and the Rush Memory and Aging Project (ROSMAP) cohort. We focused only on excitatory neurons, the dominant cell type in the scRNA-seq data. Results We identified 196 sites that harbored at least 1 individual with an excitatory neuron–specific somatic mutation (ENSM), and these 196 sites were mapped to 127 genes. The single base substitution (SBS) pattern of the putative ENSMs was best explained by signature SBS5 from the Catalogue of Somatic Mutations in Cancer (COSMIC) mutational signatures, a clock-like pattern correlating with the age of the individual. The count of ENSMs per individual also showed an increasing trend with age. Among the mutated sites, we found 2 sites tend to have more mutations in older individuals (16:6899517 [RBFOX1], p = 0.04; 4:21788463 [KCNIP4], p < 0.05). In addition, 2 sites were found to have a higher odds ratio to detect a somatic mutation in AD samples (6:73374221 [KCNQ5], p = 0.01 and 13:36667102 [DCLK1], p = 0.02). Thirty-two genes that harbor somatic mutations unique to AD and the KCNQ5 and DCLK1 genes were used for gene ontology (GO)–term enrichment analysis. We found the AD-specific ENSMs enriched in the GO-term “vocalization behavior” and “intraspecies interaction between organisms.” Of interest we observed both age-specific and AD-specific ENSMs enriched in the K + channel–associated genes. Discussion Our results show that combining scRNA-seq and WGS data can successfully detect putative somatic mutations. The putative somatic mutations detected from ROSMAP data set have provided new insights into the association of AD and aging with brain somatic mutagenesis. ...
Journal article (2023) - Niccolo' Tesi, Sven van der Lee, Marc Hulsman, Henne Holstege, Marcel Reinders
The authors of the manuscript Identification of five potential predictive biomarkers for Alzheimer’s disease by integrating the unified test for molecular signatures and weighted gene co-expression network analysis in this issue of the Medical Sciences Section of the Journals of Gerontology Series A have exploited a number of different bioinformatic tools to analyze the vast amount of data that they were confronted with (1). This is increasingly happening as it is becoming easier and cheaper to generate comprehensive molecular and phenotypic data with which biological hypotheses can be sharpened. Here, we will give a basic understanding of the methods used by Zhou et al., which are becoming standard practices when analyzing high-throughput biological data (1). ...
Journal article (2023) - Emily Baker, Ganna Leonenko, Karl Michael Schmidt, Matthew Hill, Amanda J. Myers, Maryam Shoai, Itziar de Rojas, Niccoló Tesi, Henne Holstege, More Authors...
Introduction Both late-onset Alzheimer’s disease (AD) and ageing have a strong genetic component. In each case, many associated variants have been discovered, but how much missing heritability remains to be discovered is debated. Variability in the estimation of SNP-based heritability could explain the differences in reported heritability. Methods We compute heritability in five large independent cohorts (N = 7,396, 1,566, 803, 12,528 and 3,963) to determine whether a consensus for the AD heritability estimate can be reached. These cohorts vary by sample size, age of cases and controls and phenotype definition. We compute heritability a) for all SNPs, b) excluding APOE region, c) excluding both APOE and genome-wide association study hit regions, and d) SNPs overlapping a microglia gene-set. Results SNP-based heritability of late onset Alzheimer’s disease is between 38 and 66% when age and genetic disease architecture are correctly accounted for. The heritability estimates decrease by 12% [SD = 8%] on average when the APOE region is excluded and an additional 1% [SD = 3%] when genome-wide significant regions were removed. A microglia gene-set explains 69–84% of our estimates of SNP-based heritability using only 3% of total SNPs in all cohorts. Conclusion The heritability of neurodegenerative disorders cannot be represented as a single number, because it is dependent on the ages of cases and controls. Genome-wide association studies pick up a large proportion of total AD heritability when age and genetic architecture are correctly accounted for. Around 13% of SNP-based heritability can be explained by known genetic loci and the remaining heritability likely resides around microglial related genes. ...
Journal article (2023) - Meng Zhang, Andrea B. Ganz, Susan Rohde, Linda Lorenz, Annemieke J.M. Rozemuller, Marcel J.T. Reinders, Marc Hulsman, Jeroen J.M. Hoozemans, Henne Holstege, More authors...
INTRODUCTION: Neuropathological substrates associated with neurodegeneration occur in brains of the oldest old. How does this affect cognitive performance?. METHODS: The 100-plus Study is an ongoing longitudinal cohort study of centenarians who self-report to be cognitively healthy; post mortem brain donation is optional. In 85 centenarian brains, we explored the correlations between the levels of 11 neuropathological substrates with ante mortem performance on 12 neuropsychological tests. RESULTS: Levels of neuropathological substrates varied: we observed levels up to Thal-amyloid beta phase 5, Braak-neurofibrillary tangle (NFT) stage V, Consortium to Establish a Registry for Alzheimer's Disease (CERAD)-neuritic plaque score 3, Thal-cerebral amyloid angiopathy stage 3, Tar-DNA binding protein 43 (TDP-43) stage 3, hippocampal sclerosis stage 1, Braak-Lewy bodies stage 6, atherosclerosis stage 3, cerebral infarcts stage 1, and cerebral atrophy stage 2. Granulovacuolar degeneration occurred in all centenarians. Some high performers had the highest neuropathology scores. DISCUSSION: Only Braak-NFT stage and limbic-predominant age-related TDP-43 encephalopathy (LATE) pathology associated significantly with performance across multiple cognitive domains. Of all cognitive tests, the clock-drawing test was particularly sensitive to levels of multiple neuropathologies. ...
Journal article (2022) - H. Holstege, M.J. Hulsman, Camille Charbonnier, Banjamin Grenier-Boley, Olivier Quenez, Detelina Grozeva, Jeroen G.J. van Rooij, M.J.T. Reinders, S.J. van der Lee, More authors...
Alzheimer’s disease (AD), the leading cause of dementia, has an estimated
heritability of approximately 70%1. The genetic component of AD has been mainly assessed using genome-wide association studies, which do not capture the risk contributed by rare variants2. Here, we compared the gene-based burden of rare damaging variants in exome sequencing data from 32,558 individuals—16,036 AD cases and 16,522 controls. Next to variants in TREM2, SORL1 and ABCA7, we observed a significant association of rare, predicted damaging variants in ATP8B4 and ABCA1 with AD risk, and a suggestive signal in ADAM10. Additionally, the rare-variant burden in RIN3, CLU, ZCWPW1 and ACE highlighted these genes as potential drivers of respective AD-genome-wide association study loci. Variants associated with the strongest effect on AD risk, in particular loss-of-function variants, are enriched in early-onset AD cases. Our results provide additional evidence for a major role for amyloid-β precursor protein processing, amyloid-β aggregation, lipid metabolism and microglial function in AD. ...
Journal article (2022) - Merel O. Mol, Jeroen G.J. van Rooij, Sven J. van der Lee, Marc Hulsman, Yolande A.L. Pijnenburg, Phillip Scheltens, Harro Seelaar, John C. van Swieten, Laura Donker Kaat, Henne Holstege
Background: Many families with clinical early-onset Alzheimer’s disease (EOAD) remain genetically unexplained. A combination of genetic factors is not standardly investigated. In addition to monogenic causes, we evaluated the possible polygenic architecture in a large series of families, to assess if genetic testing of familial EOAD could be expanded. Methods: Thirty-six pedigrees (77 patients) were ascertained from a larger cohort of patients, with relationships determined by genetic data (exome sequencing data and/or SNP arrays). All families included at least one AD patient with symptom onset <70 years. We evaluated segregating rare variants in known dementia-related genes, and other genes or variants if shared by multiple families. APOE was genotyped and duplications in APP were assessed by targeted test or using SNP array data. We computed polygenic risk scores (PRS) compared with a reference population-based dataset, by imputing SNP arrays or exome sequencing data. Results: In eight families, we identified a pathogenic variant, including the genes APP, PSEN1, SORL1, and an unexpected GRN frameshift variant. APOE-ε4 homozygosity was present in eighteen families, showing full segregation with disease in seven families. Eight families harbored a variant of uncertain significance (VUS), of which six included APOE-ε4 homozygous carriers. PRS was not higher in the families combined compared with the population mean (beta 0.05, P = 0.21), with a maximum increase of 0.61 (OR = 1.84) in the GRN family. Subgroup analyses indicated lower PRS in six APP/PSEN1 families compared with the rest (beta −0.22 vs. 0.10; P = 0.009) and lower APOE burden in all eight families with monogenic cause (beta 0.29 vs. 1.15, P = 0.010). Nine families remained without a genetic cause or risk factor identified. Conclusion: Besides monogenic causes, we suspect a polygenic disease architecture in multiple families based on APOE and rare VUS. The risk conveyed by PRS is modest across the studied families. Families without any identified risk factor render suitable candidates for further in-depth genetic evaluation. ...