Rare genetic variation in fibronectin 1 (FN1) protects against APOEε4 in Alzheimer's disease

The risk of developing Alzheimer's disease (AD) significantly increases in individuals carrying the APOEε4 allele. Elderly cognitively healthy individuals with APOEε4 also exist, suggesting the presence of cellular mechanisms that counteract the pathological effects of APOEε4; however, these mechanisms are unknown. We hypothesized that APOEε4 carriers without dementia might carry genetic variations that could protect them from developing APOEε4-mediated AD pathology. To test this, we leveraged whole-genome sequencing (WGS) data in the National Institute on Aging Alzheimer's Disease Family Based Study (NIA-AD FBS), Washington Heights/Inwood Columbia Aging Project (WHICAP), and Estudio Familiar de Influencia Genetica en Alzheimer (EFIGA) cohorts and identified potentially protective variants segregating exclusively among unaffected APOEε4 carriers. In homozygous unaffected carriers above 70 years old, we identified 510 rare coding variants. Pathway analysis of the genes harboring these variants showed significant enrichment in extracellular matrix (ECM)-related processes, suggesting protective effects of functional modifications in ECM proteins. We prioritized two genes that were highly represented in the ECM-related gene ontology terms, (FN1) and collagen type VI alpha 2 chain (COL6A2) and are known to be expressed at the blood-brain barrier (BBB), for postmortem validation and in vivo functional studies. An independent analysis in a large cohort of 7185 APOEε4 homozygous carriers found that rs140926439 variant in FN1 was protective of AD (OR = 0.29; 95% CI [0.11, 0.78], P = 0.014) and delayed age at onset of disease by 3.37 years (95% CI [0.42, 6.32], P = 0.025). The FN1 and COL6A2 protein levels were increased at the BBB in APOEε4 carriers with AD. Brain expression of cognitively unaffected homozygous APOEε4 carriers had significantly lower FN1 deposition and less reactive gliosis compared to homozygous APOEε4 carriers with AD, suggesting that FN1 might be a downstream driver of APOEε4-mediated AD-related pathology and cognitive decline. To validate our findings, we used zebrafish models with loss-of-function (LOF) mutations in fn1b-the ortholog for human FN1. We found that fibronectin LOF reduced gliosis, enhanced gliovascular remodeling, and potentiated the microglial response, suggesting that pathological accumulation of FN1 could impair toxic protein clearance, which is ameliorated with FN1 LOF. Our study suggests that vascular deposition of FN1 is related to the pathogenicity of APOEε4, and LOF variants in FN1 may reduce APOEε4-related AD risk, providing novel clues to potential therapeutic interventions targeting the ECM to mitigate AD risk.

APOE loss-of-function variants: Compatible with longevity and associated with resistance to Alzheimer's disease pathology

The ε4 allele of apolipoprotein E (APOE) is the strongest genetic risk factor for sporadic Alzheimer's disease (AD). Knockdown of ε4 may provide a therapeutic strategy for AD, but the effect of APOE loss of function (LoF) on AD pathogenesis is unknown. We searched for APOE LoF variants in a large cohort of controls and patients with AD and identified seven heterozygote carriers of APOE LoF variants. Five carriers were controls (aged 71-90 years), one carrier was affected by progressive supranuclear palsy, and one carrier was affected by AD with an unremarkable age at onset of 75 years. Two APOE ε3/ε4 controls carried a stop-gain affecting ε4: one was cognitively normal at 90 years and had no neuritic plaques at autopsy; the other was cognitively healthy at 79 years, and lumbar puncture at 76 years showed normal levels of amyloid. These results suggest that ε4 drives AD risk through the gain of abnormal function and support ε4 knockdown as a viable therapeutic option.

In silico identification of putative causal genetic variants

Understanding the causal genetic architecture of complex phenotypes is essential for future research into disease mechanisms and potential therapies. Despite the widespread availability of genome-wide data, existing methods to analyze genetic data still primarily focus on marginal association models, which fall short of fully capturing the polygenic nature of complex traits and elucidating biological causal mechanisms. Here we present a computationally efficient causal inference framework for genome-wide detection of putative causal variants underlying genetic associations. Our approach utilizes summary statistics from potentially overlapping studies as input, constructs in silico knockoff copies of summary statistics as negative controls to attenuate confounding effects induced by linkage disequilibrium, and employs efficient ultrahigh-dimensional sparse regression to jointly model all genetic variants across the genome. Our method is computationally efficient, requiring less than 15 minutes on a single CPU to analyze genome-wide summary statistics. In applications to a meta-analysis of ten large-scale genetic studies of Alzheimer's disease (AD) we identified 82 loci associated with AD, including 37 additional loci missed by conventional GWAS pipeline via marginal association testing. The identified putative causal variants achieve state-of-the-art agreement with massively parallel reporter assays and CRISPR-Cas9 experiments. Additionally, we applied the method to a retrospective analysis of large-scale genome-wide association studies (GWAS) summary statistics from 2013 to 2022. Results reveal the method's capacity to robustly discover additional loci for polygenic traits beyond conventional GWAS and pinpoint potential causal variants underpinning each locus (on average, 22.7% more loci and 78.7% fewer proxy variants), contributing to a deeper understanding of complex genetic architectures in post-GWAS analyses. We are making the discoveries and software freely available to the community and anticipate that routine end-to-end in silico identification of putative causal genetic variants will become an important tool that will facilitate downstream functional experiments and future research into disease etiology, as well as the exploration of novel therapeutic avenues.

Post-translational modifications linked to preclinical Alzheimer's disease-related pathological and cognitive changes

INTRODUCTION

In this study, we leverage proteomic techniques to identify communities of proteins underlying Alzheimer's disease (AD) risk among clinically unimpaired (CU) older adults.

METHODS

We constructed a protein co-expression network using 3869 cerebrospinal fluid (CSF) proteins quantified by SomaLogic, Inc., in a cohort of participants along the AD clinical spectrum. We then replicated this network in an independent cohort of CU older adults and related these modules to clinically-relevant outcomes.

RESULTS

We discovered modules enriched for phosphorylation and ubiquitination that were associated with abnormal amyloid status, as well as p-tau181 (M4: β = 2.44, p < 0.001, M7: β = 2.57, p < 0.001) and executive function performance (M4: β = -2.00, p = 0.005, M7: β = -2.39, p < 0.001).

DISCUSSION

In leveraging CSF proteomic data from individuals spanning the clinical spectrum of AD, we highlight the importance of post-translational modifications for early cognitive and pathological changes.

APOE Genotype and Alzheimer Disease Risk Across Age, Sex, and Population Ancestry

Importance

Apolipoprotein E (APOE)*2 and APOE*4 are, respectively, the strongest protective and risk-increasing, common genetic variants for late-onset Alzheimer disease (AD), making APOE status highly relevant toward clinical trial design and AD research broadly. The associations of APOE genotypes with AD are modulated by age, sex, race and ethnicity, and ancestry, but these associations remain unclear, particularly among racial and ethnic groups understudied in the AD and genetics research fields.

Objective

To assess the stratified associations of APOE genotypes with AD risk across sex, age, race and ethnicity, and global population ancestry.

Design, Setting, Participants

This genetic association study included case-control, family-based, population-based, and longitudinal AD-related cohorts that recruited referred and volunteer participants. Data were analyzed between March 2022 and April 2023. Genetic data were available from high-density, single-nucleotide variant microarrays, exome microarrays, and whole-exome and whole-genome sequencing. Summary statistics were ascertained from published AD genetic studies.

Main Outcomes and Measures

The main outcomes were risk for AD (odds ratios [ORs]) and risk of conversion to AD (hazard ratios [HRs]), with 95% CIs. Risk for AD was evaluated through case-control logistic regression analyses. Risk of conversion to AD was evaluated through Cox proportional hazards regression survival analyses.

Results

Among 68 756 unique individuals, analyses included 21 852 East Asian (demographic data not available), 5738 Hispanic (68.2% female; mean [SD] age, 75.4 [8.8] years), 7145 non-Hispanic Black (hereafter referred to as Black) (70.8% female; mean [SD] age, 78.4 [8.2] years), and 34 021 non-Hispanic White (hereafter referred to as White) (59.3% female; mean [SD] age, 77.0 [9.1] years) individuals. There was a general, stepwise pattern of ORs for APOE*4 genotypes and AD risk across race and ethnicity groups. Odds ratios for APOE*34 and AD risk attenuated following East Asian (OR, 4.54; 95% CI, 3.99-5.17),White (OR, 3.46; 95% CI, 3.27-3.65), Black (OR, 2.18; 95% CI, 1.90-2.49) and Hispanic (OR, 1.90; 95% CI, 1.65-2.18) individuals. Similarly, ORs for APOE*22+23 and AD risk attenuated following White (OR, 0.53, 95% CI, 0.48-0.58), Black (OR, 0.69, 95% CI, 0.57-0.84), and Hispanic (OR, 0.89; 95% CI, 0.72-1.10) individuals, with no association for Hispanic individuals. Deviating from the global pattern of ORs, APOE*22+23 was not associated with AD risk in East Asian individuals (OR, 0.97; 95% CI, 0.77-1.23). Global population ancestry could not explain why Hispanic individuals showed APOE associations with less pronounced AD risk compared with Black and White individuals. Within Black individuals, decreased global African ancestry or increased global European ancestry showed a pattern of APOE*4 dosage associated with increasing AD risk, but no such pattern was apparent for APOE*2 dosage with AD risk. The sex-by-age-specific interaction effect of APOE*34 among White individuals (higher risk in women) was reproduced but shifted to ages 60 to 70 years (OR, 1.48; 95% CI, 1.10-2.01) and was additionally replicated in a meta-analysis of Black individuals and Hispanic individuals (OR, 1.72; 95% CI, 1.01-2.94).

Conclusion and Relevance

Through recent advances in AD-related genetic cohorts, this study provided the largest-to-date overview of the association of APOE with AD risk across age, sex, race and ethnicity, and population ancestry. These novel insights are critical to guide AD clinical trial design and research.

A 3'UTR Insertion Is a Candidate Causal Variant at the TMEM106B Locus Associated with Increased Risk for FTLD-TDP

Background and Objectives

Single nucleotide variants near TMEM106B associate with risk of frontotemporal lobar dementia with TDP-43 inclusions (FTLD-TDP) and Alzheimer's disease (AD) in genome-wide association studies (GWAS), but the causal variant at this locus remains unclear. Here we asked whether a novel structural variant on TMEM106B is the causal variant.

Methods

An exploratory analysis identified structural variants on neurodegeneration-related genes. Subsequent analyses focused on an Alu element insertion on the 3'UTR of TMEM106B. This study included data from longitudinal aging and neurogenerative disease cohorts at Stanford University, case-control cohorts in the Alzheimer's Disease Sequencing Project (ADSP), and expression and proteomics data from Washington University in St. Louis (WUSTL). 432 individuals from two Stanford aging cohorts were whole-genome long-read and short-read sequenced. 16,906 samples from ADSP were short-read sequenced. Genotypes, transcriptomics, and proteomics data were available in 1,979 participants from an aging and dementia cohort at WUSTL. Selection criteria were specific to each cohort. In primary analyses, the linkage disequilibrium between the TMEM106B locus variants in the FTLD-TDP GWAS and the 3'UTR insertion was estimated. We then estimated linkage by ancestry in the ADSP and evaluated the effect of the TMEM106B lead variant on mRNA and protein levels.

Results

The primary analysis included 432 participants (52.5% females, age range 45-92 years old). We identified a 316 bp Alu insertion overlapping the TMEM106B 3'UTR tightly linked with top GWAS variants rs3173615(C) and rs1990622(A). In ADSP European-ancestry participants, this insertion is in equivalent linkage with rs1990622(A) (R2=0.962, D'=0.998) and rs3173615(C) (R2=0.960, D'=0.996). In African-ancestry participants, the insertion is in stronger linkage with rs1990622(A) (R2=0.992, D'=0.998) than with rs3173615(C) (R2=0.811, D'=0.994). In public datasets, rs1990622 was consistently associated with TMEM106B protein levels but not with mRNA expression. In the WUSTL dataset, rs1990622 is associated with TMEM106B protein levels in plasma and cerebrospinal fluid, but not with TMEM106B mRNA expression.

Discussion

We identified a novel Alu element insertion in the 3'UTR of TMEM106B in tight linkage with the lead FTLD-TDP risk variant. The lead variant is associated with TMEM106B protein levels, but not expression. The 3'UTR insertion is a lead candidate for the causal variant at this complex locus, pending confirmation with functional studies.

Improving genetic risk prediction across diverse population by disentangling ancestry representations

Risk prediction models using genetic data have seen increasing traction in genomics. However, most of the polygenic risk models were developed using data from participants with similar (mostly European) ancestry. This can lead to biases in the risk predictors resulting in poor generalization when applied to minority populations and admixed individuals such as African Americans. To address this issue, largely due to the prediction models being biased by the underlying population structure, we propose a deep-learning framework that leverages data from diverse population and disentangles ancestry from the phenotype-relevant information in its representation. The ancestry disentangled representation can be used to build risk predictors that perform better across minority populations. We applied the proposed method to the analysis of Alzheimer's disease genetics. Comparing with standard linear and nonlinear risk prediction methods, the proposed method substantially improves risk prediction in minority populations, including admixed individuals, without needing self-reported ancestry information.

APOE-ε4 and BIN1 increase risk of Alzheimer's disease pathology but not specifically of Lewy body pathology

Lewy body (LB) pathology commonly occurs in individuals with Alzheimer's disease (AD) pathology. However, it remains unclear which genetic risk factors underlie AD pathology, LB pathology, or AD-LB co-pathology. Notably, whether APOE-ε4 affects risk of LB pathology independently from AD pathology is controversial. We adapted criteria from the literature to classify 4,985 subjects from the National Alzheimer's Coordinating Center (NACC) and the Rush University Medical Center as AD-LB co-pathology (AD+LB+), sole AD pathology (AD+LB-), sole LB pathology (AD-LB+), or no pathology (AD-LB-). We performed a meta-analysis of a genome-wide association study (GWAS) per subpopulation (NACC/Rush) for each disease phenotype compared to the control group (AD-LB-), and compared the AD+LB+ to AD+LB- groups. APOE-ε4 was significantly associated with risk of AD+LB- and AD+LB+ compared to AD-LB-. However, APOE-ε4 was not associated with risk of AD-LB+ compared to AD-LB- or risk of AD+LB+ compared to AD+LB-. Associations at the BIN1 locus exhibited qualitatively similar results. These results suggest that APOE-ε4 is a risk factor for AD pathology, but not for LB pathology when decoupled from AD pathology. The same holds for BIN1 risk variants. These findings, in the largest AD-LB neuropathology GWAS to date, distinguish the genetic risk factors for sole and dual AD-LB pathology phenotypes. Our GWAS meta-analysis summary statistics, derived from phenotypes based on postmortem pathologic evaluation, may provide more accurate disease-specific polygenic risk scores compared to GWAS based on clinical diagnoses, which are likely confounded by undetected dual pathology and clinical misdiagnoses of dementia type.

Multiancestry analysis of the HLA locus in Alzheimer's and Parkinson's diseases uncovers a shared adaptive immune response mediated by HLA-DRB1*04 subtypes

Across multiancestry groups, we analyzed Human Leukocyte Antigen (HLA) associations in over 176,000 individuals with Parkinson's disease (PD) and Alzheimer's disease (AD) versus controls. We demonstrate that the two diseases share the same protective association at the HLA locus. HLA-specific fine-mapping showed that hierarchical protective effects of HLA-DRB1*04 subtypes best accounted for the association, strongest with HLA-DRB1*04:04 and HLA-DRB1*04:07, and intermediary with HLA-DRB1*04:01 and HLA-DRB1*04:03. The same signal was associated with decreased neurofibrillary tangles in postmortem brains and was associated with reduced tau levels in cerebrospinal fluid and to a lower extent with increased Aβ42. Protective HLA-DRB1*04 subtypes strongly bound the aggregation-prone tau PHF6 sequence, however only when acetylated at a lysine (K311), a common posttranslational modification central to tau aggregation. An HLA-DRB1*04-mediated adaptive immune response decreases PD and AD risks, potentially by acting against tau, offering the possibility of therapeutic avenues.

Sex- and APOE-specific genetic risk factors for late-onset Alzheimer's disease: Evidence from gene-gene interaction of longevity-related loci

Advanced age is the largest risk factor for late-onset Alzheimer's disease (LOAD), a disease in which susceptibility correlates to almost all hallmarks of aging. Shared genetic signatures between LOAD and longevity were frequently hypothesized, likely characterized by distinctive epistatic and pleiotropic interactions. Here, we applied a multidimensional reduction approach to detect gene-gene interactions affecting LOAD in a large dataset of genomic variants harbored by genes in the insulin/IGF1 signaling, DNA repair, and oxidative stress pathways, previously investigated in human longevity. The dataset was generated from a collection of publicly available Genome Wide Association Studies, comprising a total of 2,469 gene variants genotyped in 20,766 subjects of Northwestern European ancestry (11,038 LOAD cases and 9,728 controls). The stratified analysis according to APOE*4 status and sex corroborated evidence that pathways leading to longevity also contribute to LOAD. Among the significantly interacting genes, PTPN1, TXNRD1, and IGF1R were already found enriched in gene-gene interactions affecting survival to old age. Furthermore, interacting variants associated with LOAD in a sex- and APOE-specific way. Indeed, while in APOE*4 female carriers we found several inter-pathway interactions, no significant epistasis was found in APOE*4 negative females; conversely, in males, significant intra- and inter-pathways epistasis emerged according to APOE*4 status. These findings suggest that interactions of risk factors may drive different trajectories of cognitive aging. Beyond helping to disentangle the genetic architecture of LOAD, such knowledge may improve precision in predicting the risk of dementia and enable effective sex- and APOE-stratified preventive and therapeutic interventions for LOAD.

APOE - ε 4 and BIN1 increase risk of Alzheimer's disease pathology but not specifically of Lewy body pathology

Lewy body (LB) pathology commonly occurs in individuals with Alzheimer's disease (AD) pathology. However, it remains unclear which genetic risk factors underlie AD pathology, LB pathology, or AD-LB co-pathology. Notably, whether APOE - ε 4 affects risk of LB pathology independently from AD pathology is controversial. We adapted criteria from the literature to classify 4,985 subjects from the National Alzheimer's Coordinating Center (NACC) and the Rush University Medical Center as AD-LB co-pathology (AD + LB + ), sole AD pathology (AD + LB - ), sole LB pathology (AD - LB + ), or no pathology (AD - LB - ). We performed a meta-analysis of a genome-wide association study (GWAS) per subpopulation (NACC/Rush) for each disease phenotype compared to the control group (AD - LB - ), and compared the AD + LB + to AD + LB - groups. APOE - ε 4 was significantly associated with risk of AD + LB - and AD + LB + compared to AD - LB - . However, APOE - ε 4 was not associated with risk of AD - LB + compared to AD - LB - or risk of AD + LB + compared to AD + LB - . Associations at the BIN1 locus exhibited qualitatively similar results. These results suggest that APOE - ε 4 is a risk factor for AD pathology, but not for LB pathology when decoupled from AD pathology. The same holds for BIN1 risk variants. These findings, in the largest AD-LB neuropathology GWAS to date, distinguish the genetic risk factors for sole and dual AD-LB pathology phenotypes. Our GWAS meta-analysis summary statistics, derived from phenotypes based on postmortem pathologic evaluation, may provide more accurate disease-specific polygenic risk scores compared to GWAS based on clinical diagnoses, which are likely confounded by undetected dual pathology and clinical misdiagnoses of dementia type.

Decoding the heterogeneity of Alzheimer's disease diagnosis and progression using multilayer networks

Alzheimer's disease (AD) is a multifactorial and heterogeneous disorder, which makes early detection a challenge. Studies have attempted to combine biomarkers to improve AD detection and predict progression. However, most of the existing work reports results in parallel or compares normalized findings but does not analyze data simultaneously. We tested a multi-dimensional network framework, applied to 490 subjects (cognitively normal [CN] = 147; mild cognitive impairment [MCI] = 287; AD = 56) from ADNI, to create a single model capable of capturing the heterogeneity and progression of AD. First, we constructed subject similarity networks for structural magnetic resonance imaging, amyloid-β positron emission tomography, cerebrospinal fluid, cognition, and genetics data and then applied multilayer community detection to find groups with shared similarities across modalities. Individuals were also followed-up longitudinally, with AD subjects having, on average, 4.5 years of follow-up. Our findings show that multilayer community detection allows for accurate identification of present and future AD (≈90%) and is also able to identify cases that were misdiagnosed clinically. From all MCI participants who developed AD or reverted to CN, the multilayer model correctly identified 90.8% and 88.5% of cases respectively. We observed similar subtypes across the full sample and when examining multimodal data from subjects with no AD pathology (i.e., amyloid negative). Finally, these results were also validated using an independent testing set. In summary, the multilayer framework is successful in detecting AD and provides unique insight into the heterogeneity of the disease by identifying subtypes that share similar multidisciplinary profiles of neurological, cognitive, pathological, and genetics information.

Association of African Ancestry-Specific APOE Missense Variant R145C With Risk of Alzheimer Disease

Importance

Numerous studies have established the association of the common APOE ε2 and APOE ε4 alleles with Alzheimer disease (AD) risk across ancestries. Studies of the interaction of these alleles with other amino acid changes on APOE in non-European ancestries are lacking and may improve ancestry-specific risk prediction.

Objective

To determine whether APOE amino acid changes specific to individuals of African ancestry modulate AD risk.

Design, Setting, and Participants

Case-control study including 31 929 participants and using a sequenced discovery sample (Alzheimer Disease Sequencing Project; stage 1) followed by 2 microarray imputed data sets derived from the Alzheimer Disease Genetic Consortium (stage 2, internal replication) and the Million Veteran Program (stage 3, external validation). This study combined case-control, family-based, population-based, and longitudinal AD cohorts, which recruited participants (1991-2022) in primarily US-based studies with 1 US/Nigerian study. Across all stages, individuals included in this study were of African ancestry.

Exposures

Two APOE missense variants (R145C and R150H) were assessed, stratified by APOE genotype.

Main Outcomes and Measures

The primary outcome was AD case-control status, and secondary outcomes included age at AD onset.

Results

Stage 1 included 2888 cases (median age, 77 [IQR, 71-83] years; 31.3% male) and 4957 controls (median age, 77 [IQR, 71-83] years; 28.0% male). In stage 2, across multiple cohorts, 1201 cases (median age, 75 [IQR, 69-81] years; 30.8% male) and 2744 controls (median age, 80 [IQR, 75-84] years; 31.4% male) were included. In stage 3, 733 cases (median age, 79.4 [IQR, 73.8-86.5] years; 97.0% male) and 19 406 controls (median age, 71.9 [IQR, 68.4-75.8] years; 94.5% male) were included. In ε3/ε4-stratified analyses of stage 1, R145C was present in 52 individuals with AD (4.8%) and 19 controls (1.5%); R145C was associated with an increased risk of AD (odds ratio [OR], 3.01; 95% CI, 1.87-4.85; P = 6.0 × 10-6) and was associated with a reported younger age at AD onset (β, -5.87 years; 95% CI, -8.35 to -3.4 years; P = 3.4 × 10-6). Association with increased AD risk was replicated in stage 2 (R145C was present in 23 individuals with AD [4.7%] and 21 controls [2.7%]; OR, 2.20; 95% CI, 1.04-4.65; P = .04) and was concordant in stage 3 (R145C was present in 11 individuals with AD [3.8%] and 149 controls [2.7%]; OR, 1.90; 95% CI, 0.99-3.64; P = .051). Association with earlier AD onset was replicated in stage 2 (β, -5.23 years; 95% CI, -9.58 to -0.87 years; P = .02) and stage 3 (β, -10.15 years; 95% CI, -15.66 to -4.64 years; P = 4.0 × 10-4). No significant associations were observed in other APOE strata for R145C or in any APOE strata for R150H.

Conclusions and Relevance

In this exploratory analysis, the APOE ε3[R145C] missense variant was associated with an increased risk of AD among individuals of African ancestry with the ε3/ε4 genotype. With additional external validation, these findings may inform AD genetic risk assessment in individuals of African ancestry.

APOE effects on regional tau in preclinical Alzheimer's disease

BACKGROUND

APOE variants are strongly associated with abnormal amyloid aggregation and additional direct effects of APOE on tau aggregation are reported in animal and human cell models. The degree to which these effects are present in humans when individuals are clinically unimpaired (CU) but have abnormal amyloid (Aβ+) remains unclear.

METHODS

We analyzed data from CU individuals in the Anti-Amyloid Treatment in Asymptomatic AD (A4) and Longitudinal Evaluation of Amyloid Risk and Neurodegeneration (LEARN) studies. Amyloid PET data were available for 4486 participants (3163 Aβ-, 1323 Aβ+) and tau PET data were available for a subset of 447 participants (55 Aβ-, 392 Aβ+). Linear models examined APOE (number of e2 and e4 alleles) associations with global amyloid and regional tau burden in medial temporal lobe (entorhinal, amygdala) and early neocortical regions (inferior temporal, inferior parietal, precuneus). Consistency of APOE4 effects on regional tau were examined in 220 Aβ + CU and mild cognitive impairment (MCI) participants from the Alzheimer's Disease Neuroimaging Initiative (ADNI).

RESULTS

APOE2 and APOE4 were associated with lower and higher amyloid positivity rates, respectively. Among Aβ+ CU, e2 and e4 were associated with reduced (-12 centiloids per allele) and greater (+15 centiloids per allele) continuous amyloid burden, respectively. APOE2 was associated with reduced regional tau in all regions (-0.05 to -0.09 SUVR per allele), whereas APOE4 was associated with greater regional tau (+0.02 to +0.07 SUVR per allele). APOE differences were confirmed by contrasting e3/e3 with e2/e3 and e3/e4. Mediation analyses among Aβ+ s showed that direct effects of e2 on regional tau were present in medial temporal lobe and early neocortical regions, beyond an indirect pathway mediated by continuous amyloid burden. For e4, direct effects on regional tau were only significant in medial temporal lobe. The magnitude of protective e2 effects on regional tau was consistent across brain regions, whereas detrimental e4 effects were greatest in medial temporal lobe. APOE4 patterns were confirmed in Aβ+ ADNI participants.

CONCLUSIONS

APOE influences early regional tau PET burden, above and beyond effects related to cross-sectional amyloid PET burden. Therapeutic strategies targeting underlying mechanisms related to APOE may modify tau accumulation among Aβ+ individuals.

GhostKnockoff inference empowers identification of putative causal variants in genome-wide association studies

Recent advances in genome sequencing and imputation technologies provide an exciting opportunity to comprehensively study the contribution of genetic variants to complex phenotypes. However, our ability to translate genetic discoveries into mechanistic insights remains limited at this point. In this paper, we propose an efficient knockoff-based method, GhostKnockoff, for genome-wide association studies (GWAS) that leads to improved power and ability to prioritize putative causal variants relative to conventional GWAS approaches. The method requires only Z-scores from conventional GWAS and hence can be easily applied to enhance existing and future studies. The method can also be applied to meta-analysis of multiple GWAS allowing for arbitrary sample overlap. We demonstrate its performance using empirical simulations and two applications: (1) a meta-analysis for Alzheimer's disease comprising nine overlapping large-scale GWAS, whole-exome and whole-genome sequencing studies and (2) analysis of 1403 binary phenotypes from the UK Biobank data in 408,961 samples of European ancestry. Our results demonstrate that GhostKnockoff can identify putatively functional variants with weaker statistical effects that are missed by conventional association tests.

Deep neural networks with controlled variable selection for the identification of putative causal genetic variants

Deep neural networks (DNNs) have been successfully utilized in many scientific problems for their high prediction accuracy, but their application to genetic studies remains challenging due to their poor interpretability. Here we consider the problem of scalable, robust variable selection in DNNs for the identification of putative causal genetic variants in genome sequencing studies. We identified a pronounced randomness in feature selection in DNNs due to its stochastic nature, which may hinder interpretability and give rise to misleading results. We propose an interpretable neural network model, stabilized using ensembling, with controlled variable selection for genetic studies. The merit of the proposed method includes: flexible modelling of the nonlinear effect of genetic variants to improve statistical power; multiple knockoffs in the input layer to rigorously control the false discovery rate; hierarchical layers to substantially reduce the number of weight parameters and activations, and improve computational efficiency; and stabilized feature selection to reduce the randomness in identified signals. We evaluate the proposed method in extensive simulation studies and apply it to the analysis of Alzheimer's disease genetics. We show that the proposed method, when compared with conventional linear and nonlinear methods, can lead to substantially more discoveries.

A Fast and Robust Strategy to Remove Variant-Level Artifacts in Alzheimer Disease Sequencing Project Data

Background and Objectives

Exome sequencing (ES) and genome sequencing (GS) are expected to be critical to further elucidate the missing genetic heritability of Alzheimer disease (AD) risk by identifying rare coding and/or noncoding variants that contribute to AD pathogenesis. In the United States, the Alzheimer Disease Sequencing Project (ADSP) has taken a leading role in sequencing AD-related samples at scale, with the resultant data being made publicly available to researchers to generate new insights into the genetic etiology of AD. To achieve sufficient power, the ADSP has adapted a study design where subsets of larger AD cohorts are collected and sequenced across multiple centers, using a variety of sequencing platforms. This approach may lead to variable variant quality across sequencing centers and/or platforms. In this study, we sought to implement and evaluate filters that can be applied fast to robustly remove variant-level artifacts in the ADSP data.

Methods

We implemented a robust quality control procedure to handle ADSP data. We evaluated this procedure while performing exome-wide and genome-wide association analyses on AD risk using the latest ADSP whole ES (WES) and whole GS (WGS) data releases (NG00067.v5).

Results

We observed that many variants displayed large variation in allele frequencies across sequencing centers/platforms and contributed to spurious association signals with AD risk. We also observed that sequencing platform/center adjustment in association models could not fully account for these spurious signals. To address this issue, we designed and implemented variant filters that could capture and remove these center-specific/platform-specific artifactual variants.

Discussion

We derived a fast and robust approach to filter variants that represent sequencing center-related or platform-related artifacts underlying spurious associations with AD risk in ADSP WES and WGS data. This approach will be important to support future robust genetic association studies on ADSP data, as well as other studies with similar designs.

Association of Rare APOE Missense Variants V236E and R251G With Risk of Alzheimer Disease

Importance

The APOE ε2 and APOE ε4 alleles are the strongest protective and risk-increasing, respectively, genetic variants for late-onset Alzheimer disease (AD). However, the mechanisms linking APOE to AD-particularly the apoE protein's role in AD pathogenesis and how this is affected by APOE variants-remain poorly understood. Identifying missense variants in addition to APOE ε2 and APOE ε4 could provide critical new insights, but given the low frequency of additional missense variants, AD genetic cohorts have previously been too small to interrogate this question robustly.

Objective

To determine whether rare missense variants on APOE are associated with AD risk.

Design, Setting, and Participants

Association with case-control status was tested in a sequenced discovery sample (stage 1) and followed up in several microarray imputed cohorts as well as the UK Biobank whole-exome sequencing resource using a proxy-AD phenotype (stages 2 and 3). This study combined case-control, family-based, population-based, and longitudinal AD-related cohorts that recruited referred and volunteer participants. Stage 1 included 37 409 nonunique participants of European or admixed European ancestry, with 11 868 individuals with AD and 11 934 controls passing analysis inclusion criteria. In stages 2 and 3, 475 473 participants were considered across 8 cohorts, of which 84 513 individuals with AD and proxy-AD and 328 372 controls passed inclusion criteria. Selection criteria were cohort specific, and this study was performed a posteriori on individuals who were genotyped. Among the available genotypes, 76 195 were excluded. All data were retrieved between September 2015 and November 2021 and analyzed between April and November 2021.

Main Outcomes and Measures

In primary analyses, the AD risk associated with each missense variant was estimated, as appropriate, with either linear mixed-model regression or logistic regression. In secondary analyses, associations were estimated with age at onset using linear mixed-model regression and risk of conversion to AD using competing-risk regression.

Results

A total of 544 384 participants were analyzed in the primary case-control analysis; 312 476 (57.4%) were female, and the mean (SD; range) age was 64.9 (15.2; 40-110) years. Two missense variants were associated with a 2-fold to 3-fold decreased AD risk: APOE ε4 (R251G) (odds ratio, 0.44; 95% CI, 0.33-0.59; P = 4.7 × 10-8) and APOE ε3 (V236E) (odds ratio, 0.37; 95% CI, 0.25-0.56; P = 1.9 × 10-6). Additionally, the cumulative incidence of AD in carriers of these variants was found to grow more slowly with age compared with noncarriers.

Conclusions and Relevance

In this genetic association study, a novel variant associated with AD was identified: R251G always coinherited with ε4 on the APOE gene, which mitigates the ε4-associated AD risk. The protective effect of the V236E variant, which is always coinherited with ε3 on the APOE gene, was also confirmed. The location of these variants confirms that the carboxyl-terminal portion of apoE plays an important role in AD pathogenesis. The large risk reductions reported here suggest that protein chemistry and functional assays of these variants should be pursued, as they have the potential to guide drug development targeting APOE.

Challenges at the APOE locus: a robust quality control approach for accurate APOE genotyping

BACKGROUND

Genetic variants within the APOE locus may modulate Alzheimer's disease (AD) risk independently or in conjunction with APOE*2/3/4 genotypes. Identifying such variants and mechanisms would importantly advance our understanding of APOE pathophysiology and provide critical guidance for AD therapies aimed at APOE. The APOE locus however remains relatively poorly understood in AD, owing to multiple challenges that include its complex linkage structure and uncertainty in APOE*2/3/4 genotype quality. Here, we present a novel APOE*2/3/4 filtering approach and showcase its relevance on AD risk association analyses for the rs439401 variant, which is located 1801 base pairs downstream of APOE and has been associated with a potential regulatory effect on APOE.

METHODS

We used thirty-two AD-related cohorts, with genetic data from various high-density single-nucleotide polymorphism microarrays, whole-genome sequencing, and whole-exome sequencing. Study participants were filtered to be ages 60 and older, non-Hispanic, of European ancestry, and diagnosed as cognitively normal or AD (n = 65,701). Primary analyses investigated AD risk in APOE*4/4 carriers. Additional supporting analyses were performed in APOE*3/4 and 3/3 strata. Outcomes were compared under two different APOE*2/3/4 filtering approaches.

RESULTS

Using more conventional APOE*2/3/4 filtering criteria (approach 1), we showed that, when in-phase with APOE*4, rs439401 was variably associated with protective effects on AD case-control status. However, when applying a novel filter that increases the certainty of the APOE*2/3/4 genotypes by applying more stringent criteria for concordance between the provided APOE genotype and imputed APOE genotype (approach 2), we observed that all significant effects were lost.

CONCLUSIONS

We showed that careful consideration of APOE genotype and appropriate sample filtering were crucial to robustly interrogate the role of the APOE locus on AD risk. Our study presents a novel APOE filtering approach and provides important guidelines for research into the APOE locus, as well as for elucidating genetic interaction effects with APOE*2/3/4.

Genome-wide analysis of common and rare variants via multiple knockoffs at biobank scale, with an application to Alzheimer disease genetics

Knockoff-based methods have become increasingly popular due to their enhanced power for locus discovery and their ability to prioritize putative causal variants in a genome-wide analysis. However, because of the substantial computational cost for generating knockoffs, existing knockoff approaches cannot analyze millions of rare genetic variants in biobank-scale whole-genome sequencing and whole-genome imputed datasets. We propose a scalable knockoff-based method for the analysis of common and rare variants across the genome, KnockoffScreen-AL, that is applicable to biobank-scale studies with hundreds of thousands of samples and millions of genetic variants. The application of KnockoffScreen-AL to the analysis of Alzheimer disease (AD) in 388,051 WG-imputed samples from the UK Biobank resulted in 31 significant loci, including 14 loci that are missed by conventional association tests on these data. We perform replication studies in an independent meta-analysis of clinically diagnosed AD with 94,437 samples, and additionally leverage single-cell RNA-sequencing data with 143,793 single-nucleus transcriptomes from 17 control subjects and AD-affected individuals, and proteomics data from 735 control subjects and affected indviduals with AD and related disorders to validate the genes at these significant loci. These multi-omics analyses show that 79.1% of the proximal genes at these loci and 76.2% of the genes at loci identified only by KnockoffScreen-AL exhibit at least suggestive signal (p < 0.05) in the scRNA-seq or proteomics analyses. We highlight a potentially causal gene in AD progression, EGFR, that shows significant differences in expression and protein levels between AD-affected individuals and healthy control subjects.

Common X-Chromosome Variants Are Associated with Parkinson Disease Risk

OBJECTIVE

The objective of this study was to identify genetic variants on the X-chromosome associated with Parkinson disease (PD) risk.

METHODS

We performed an X-chromosome-wide association study (XWAS) of PD risk by meta-analyzing results from sex-stratified analyses. To avoid spurious associations, we designed a specific harmonization pipeline for the X-chromosome and focused on a European ancestry sample. We included 11,142 cases, 280,164 controls, and 5,379 proxy cases, based on parental history of PD. Additionally, we tested the association of significant variants with (1) PD risk in an independent replication with 1,561 cases and 2,465 controls and (2) putamen volume in 33,360 individuals from the UK Biobank.

RESULTS

In the discovery meta-analysis, we identified rs7066890 (odds ratio [OR] = 1.10, 95% confidence interval [CI] = 1.06-1.14, p = 2.2 × 10-9 ), intron of GPM6B, and rs28602900 (OR = 1.10, 95% CI = 1.07-1.14, p = 1.6 × 10-8 ) in a high gene density region including RPL10, ATP6A1, FAM50A, and PLXNA3. The rs28602900 association with PD was replicated (OR = 1.16, 95% CI = 1.03-1.30, p = 0.016) and shown to colocalize with a significant expression quantitative locus (eQTL) regulating RPL10 expression in the putamen and other brain tissues in the Genotype-Tissue Expression Project. Additionally, the rs28602900 locus was found to be associated with reduced brain putamen volume. No results reached genome-wide significance in the sex-stratified analyses.

INTERPRETATION

We report the first XWAS of PD and identify 2 genome-wide significant loci. The rs28602900 association was replicated in an independent PD dataset and showed concordant effects in its association with putamen volume. Critically, rs26802900 is a significant eQTL of RPL10. These results support a role for ribosomal proteins in PD pathogenesis and show that the X-chromosome contributes to PD genetic risk. ANN NEUROL 2021;90:22-34.

Long-term ovarian hormone deprivation alters functional connectivity, brain neurochemical profile and white matter integrity in the Tg2576 amyloid mouse model of Alzheimer's disease

Premenopausal bilateral ovariectomy is considered to be one of the risk factors of Alzheimer's disease (AD). However, the underlying mechanisms remain unclear. Here, we aimed to investigate long-term neurological consequences of ovariectomy in a rodent AD model, TG2576 (TG), and wild-type mice (WT) that underwent an ovariectomy or sham-operation, using in vivo MRI biomarkers. An increase in osmoregulation and energy metabolism biomarkers in the hypothalamus, a decrease in white matter integrity, and a decrease in the resting-state functional connectivity was observed in ovariectomized TG mice compared to sham-operated TG mice. In addition, we observed an increase in functional connectivity in ovariectomized WT mice compared to sham-operated WT mice. Furthermore, genotype (TG vs. WT) effects on imaging markers and GFAP immunoreactivity levels were observed, but there was no effect of interaction (Genotype × Surgery) on amyloid-beta-and GFAP immunoreactivity levels. Taken together, our results indicated that both genotype and ovariectomy alters imaging biomarkers associated with AD.

A novel age-informed approach for genetic association analysis in Alzheimer's disease

BACKGROUND

Many Alzheimer's disease (AD) genetic association studies disregard age or incorrectly account for it, hampering variant discovery.

METHODS

Using simulated data, we compared the statistical power of several models: logistic regression on AD diagnosis adjusted and not adjusted for age; linear regression on a score integrating case-control status and age; and multivariate Cox regression on age-at-onset. We applied these models to real exome-wide data of 11,127 sequenced individuals (54% cases) and replicated suggestive associations in 21,631 genotype-imputed individuals (51% cases).

RESULTS

Modeling variable AD risk across age results in 5-10% statistical power gain compared to logistic regression without age adjustment, while incorrect age adjustment leads to critical power loss. Applying our novel AD-age score and/or Cox regression, we discovered and replicated novel variants associated with AD on KIF21B, USH2A, RAB10, RIN3, and TAOK2 genes.

CONCLUSION

Our AD-age score provides a simple means for statistical power gain and is recommended for future AD studies.

Association of Klotho-VS Heterozygosity With Risk of Alzheimer Disease in Individuals Who Carry APOE4

Importance

Identification of genetic factors that interact with the apolipoprotein e4 (APOE4) allele to reduce risk for Alzheimer disease (AD) would accelerate the search for new AD drug targets. Klotho-VS heterozygosity (KL-VSHET+ status) protects against aging-associated phenotypes and cognitive decline, but whether it protects individuals who carry APOE4 from AD remains unclear.

Objectives

To determine if KL-VSHET+ status is associated with reduced AD risk and β-amyloid (Aβ) pathology in individuals who carry APOE4.

Design, Setting, and Participants

This study combined 25 independent case-control, family-based, and longitudinal AD cohorts that recruited referred and volunteer participants and made data available through public repositories. Analyses were stratified by APOE4 status. Three cohorts were used to evaluate conversion risk, 1 provided longitudinal measures of Aβ CSF and PET, and 3 provided cross-sectional measures of Aβ CSF. Genetic data were available from high-density single-nucleotide variant microarrays. All data were collected between September 2015 and September 2019 and analyzed between April 2019 and December 2019.

Main Outcomes and Measures

The risk of AD was evaluated through logistic regression analyses under a case-control design. The risk of conversion to mild cognitive impairment (MCI) or AD was evaluated through competing risks regression. Associations with Aβ, measured from cerebrospinal fluid (CSF) or brain positron emission tomography (PET), were evaluated using linear regression and mixed-effects modeling.

Results

Of 36 530 eligible participants, 13 782 were excluded for analysis exclusion criteria or refusal to participate. Participants were men and women aged 60 years and older who were non-Hispanic and of Northwestern European ancestry and had been diagnosed as being cognitively normal or having MCI or AD. The sample included 20 928 participants in case-control studies, 3008 in conversion studies, 556 in Aβ CSF regression analyses, and 251 in PET regression analyses. The genotype KL-VSHET+ was associated with reduced risk for AD in individuals carrying APOE4 who were 60 years or older (odds ratio, 0.75 [95% CI, 0.67-0.84]; P = 7.4 × 10-7), and this was more prominent at ages 60 to 80 years (odds ratio, 0.69 [95% CI, 0.61-0.79]; P = 3.6 × 10-8). Additionally, control participants carrying APOE4 with KL-VS heterozygosity were at reduced risk of converting to MCI or AD (hazard ratio, 0.64 [95% CI, 0.44-0.94]; P = .02). Finally, in control participants who carried APOE4 and were aged 60 to 80 years, KL-VS heterozygosity was associated with higher Aβ in CSF (β, 0.06 [95% CI, 0.01-0.10]; P = .03) and lower Aβ on PET scans (β, -0.04 [95% CI, -0.07 to -0.00]; P = .04).

Conclusions and Relevance

The genotype KL-VSHET+ is associated with reduced AD risk and Aβ burden in individuals who are aged 60 to 80 years, cognitively normal, and carrying APOE4. Molecular pathways associated with KL merit exploration for novel AD drug targets. The KL-VS genotype should be considered in conjunction with the APOE genotype to refine AD prediction models used in clinical trial enrichment and personalized genetic counseling.

KL∗VS heterozygosity reduces brain amyloid in asymptomatic at-risk APOE∗4 carriers

KLOTHO∗VS heterozygosity (KL∗VS^HET+) was recently shown to be associated with reduced risk of Alzheimer's disease (AD) in APOE∗4 carriers. Additional studies suggest that KL∗VS^HET+ protects against amyloid burden in cognitively normal older subjects, but sample sizes were too small to draw definitive conclusions. We performed a well-powered meta-analysis across 5 independent studies, comprising 3581 pre-clinical participants ages 60-80, to investigate whether KL∗VS^HET+ reduces the risk of having an amyloid-positive positron emission tomography scan. Analyses were stratified by APOE∗4 status. KL∗VS^HET+ reduced the risk of amyloid positivity in APOE∗4 carriers (odds ratio = 0.67 [0.52-0.88]; p = 3.5 × 10^-3), but not in APOE∗4 non-carriers (odds ratio = 0.94 [0.73-1.21]; p = 0.63). The combination of APOE∗4 and KL∗VS genotypes should help enrich AD clinical trials for pre-symptomatic subjects at increased risk of developing amyloid aggregation and AD. KL-related pathways may help elucidate protective mechanisms against amyloid accumulation and merit exploration for novel AD drug targets. Future investigation of the biological mechanisms by which KL interacts with APOE∗4 and AD are warranted.