Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC v1.0 BeadChip microarrays

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC v1.0 arrays. We conducted a comprehensive assessment of the EPIC v1.0 array probe reliability using 69 blood DNA samples, each measured twice, generated by the Alzheimer's Disease Neuroimaging Initiative study. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliability information for probes on the EPIC v1.0 array, will serve as a valuable resource for future DNAm studies.

Genome Wide Association Study of Neuropathic Ocular Pain

Purpose

To conduct a genome-wide association study (GWAS) of individuals with neuropathic ocular pain (NOP) symptoms to identify genomic variants that may predispose to NOP development.

Design

Prospective study of individuals with NOP.

Participants

Three hundred twenty-nine patients recruited from the Miami Veterans Affairs eye clinic.

Methods

The Neuropathic Pain Symptom Inventory modified for the eye (NPSI-Eye) was completed to calculate a NPSI-Eye-Sub-Score (summed ratings of burning and wind sensitivity) as an indicator of NOP severity. A GWAS was performed for the NPSI-Eye-Sub-Score with a significance threshold of P < 5 × 10-8. A gene-based analysis was performed using the multimarker analysis of genomic annotation software (in the functional mapping and annotation of GWAS online platform). The 13 865 778 single nucleotide polymorphisms (SNPs) from our GWAS analysis were mapped to 10 834 protein coding genes, and significant genes were run through gene set enrichment analysis.

Main Outcome Measures

Identification of SNPs and protein products that may be associated with the development of NOP.

Results

One hundred seventy-one SNPs reached a threshold of P < 10-5, of which 10 SNPs reached the suggestive level of significance of P < 5 × 10-7 and 1 SNP met our genome-wide significance threshold of P < 5 × 10-8. This lead SNP, rs140293404 (P = 1.23 × 10-8), is an intronic variant found within gene ENSG00000287251 coding for transcript ENST00000662732.1. Rs140293404 is in linkage disequilibrium with exon variant rs7926353 (r2 > 0.8) within ENSG00000279046 coding for transcript ENST00000624288.1. The most significant genes from gene-based tests were matrix metalloproteinase-19 (MMP19) (P = 1.12 × 10-5), zinc finger RNA-binding motif and serine/arginine rich-1 (ZRSR1) (P = 1.48 × 10-4), CTC-487M23.8 (P = 1.79 × 10-4), receptor expression-enhancing protein-5 (REEP5) (P = 2.36 × 10-4), and signal recognition particle-19 (SRP19) (P = 2.56 × 10-4). From gene set enrichment analysis, the sensory perception (false discovery rate = 6.57 × 10-3) and olfactory signaling (false discovery rate = 1.63 × 10-2) pathways were enriched with the most significant genes.

Conclusions

Our GWAS revealed genes with protein products that may impact sensory perception, lending biological plausibility to a role for SNPs identified by our GWAS in the development of NOP. A better understanding of the biological relevance of these genes and pathways in the pathophysiology associated with NOP may facilitate future novel mechanism-based treatments.

Financial Disclosures

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Sex-specific genetic architecture of late-life memory performance

BACKGROUND

Women demonstrate a memory advantage when cognitively healthy yet lose this advantage to men in Alzheimer's disease. However, the genetic underpinnings of this sex difference in memory performance remain unclear.

METHODS

We conducted the largest sex-aware genetic study on late-life memory to date (Nmales  = 11,942; Nfemales  = 15,641). Leveraging harmonized memory composite scores from four cohorts of cognitive aging and AD, we performed sex-stratified and sex-interaction genome-wide association studies in 24,216 non-Hispanic White and 3367 non-Hispanic Black participants.

RESULTS

We identified three sex-specific loci (rs67099044-CBLN2, rs719070-SCHIP1/IQCJ-SCHIP), including an X-chromosome locus (rs5935633-EGL6/TCEANC/OFD1), that associated with memory. Additionally, we identified heparan sulfate signaling as a sex-specific pathway and found sex-specific genetic correlations between memory and cardiovascular, immune, and education traits.

DISCUSSION

This study showed memory is highly and comparably heritable across sexes, as well as highlighted novel sex-specific genes, pathways, and genetic correlations that related to late-life memory.

HIGHLIGHTS

Demonstrated the heritable component of late-life memory is similar across sexes. Identified two genetic loci with a sex-interaction with baseline memory. Identified an X-chromosome locus associated with memory decline in females. Highlighted sex-specific candidate genes and pathways associated with memory. Revealed sex-specific shared genetic architecture between memory and complex traits.

Longitudinal change in memory performance as a strong endophenotype for Alzheimer's disease

INTRODUCTION

Although large-scale genome-wide association studies (GWAS) have been conducted on AD, few have been conducted on continuous measures of memory performance and memory decline.

METHODS

We conducted a cross-ancestry GWAS on memory performance (in 27,633 participants) and memory decline (in 22,365 participants; 129,201 observations) by leveraging harmonized cognitive data from four aging cohorts.

RESULTS

We found high heritability for two ancestry backgrounds. Further, we found a novel ancestry locus for memory decline on chromosome 4 (rs6848524) and three loci in the non-Hispanic Black ancestry group for memory performance on chromosomes 2 (rs111471504), 7 (rs4142249), and 15 (rs74381744). In our gene-level analysis, we found novel genes for memory decline on chromosomes 1 (SLC25A44), 11 (BSX), and 15 (DPP8). Memory performance and memory decline shared genetic architecture with AD-related traits, neuropsychiatric traits, and autoimmune traits.

DISCUSSION

We discovered several novel loci, genes, and genetic correlations associated with late-life memory performance and decline.

HIGHLIGHTS

Late-life memory has high heritability that is similar across ancestries. We discovered four novel variants associated with late-life memory. We identified four novel genes associated with late-life memory. Late-life memory shares genetic architecture with psychiatric/autoimmune traits.

An X Chromosome Transcriptome Wide Association Study Implicates ARMCX6 in Alzheimer's Disease

Background

The X chromosome is often omitted in disease association studies despite containing thousands of genes that may provide insight into well-known sex differences in the risk of Alzheimer's disease (AD).

Objective

To model the expression of X chromosome genes and evaluate their impact on AD risk in a sex-stratified manner.

Methods

Using elastic net, we evaluated multiple modeling strategies in a set of 175 whole blood samples and 126 brain cortex samples, with whole genome sequencing and RNA-seq data. SNPs (MAF > 0.05) within the cis-regulatory window were used to train tissue-specific models of each gene. We apply the best models in both tissues to sex-stratified summary statistics from a meta-analysis of Alzheimer's Disease Genetics Consortium (ADGC) studies to identify AD-related genes on the X chromosome.

Results

Across different model parameters, sample sex, and tissue types, we modeled the expression of 217 genes (95 genes in blood and 135 genes in brain cortex). The average model R2 was 0.12 (range from 0.03 to 0.34). We also compared sex-stratified and sex-combined models on the X chromosome. We further investigated genes that escaped X chromosome inactivation (XCI) to determine if their genetic regulation patterns were distinct. We found ten genes associated with AD at p < 0.05, with only ARMCX6 in female brain cortex (p = 0.008) nearing the significance threshold after adjusting for multiple testing (α = 0.002).

Conclusions

We optimized the expression prediction of X chromosome genes, applied these models to sex-stratified AD GWAS summary statistics, and identified one putative AD risk gene, ARMCX6.

DNA from multiple viral species is associated with Alzheimer's disease risk

INTRODUCTION

Multiple infectious agents, including viruses, bacteria, fungi, and protozoa, have been linked to Alzheimer's disease (AD) risk by independent lines of evidence. We explored this association by comparing the frequencies of viral species identified in a large sample of AD cases and controls.

METHODS

DNA sequence reads that did not align to the human genome in sequences were mapped to viral reference sequences, quantified, and then were tested for association with AD in whole exome sequences (WES) and whole genome sequences (WGS) datasets.

RESULTS

Several viruses were significant predictors of AD according to the machine learning classifiers. Subsequent regression analyses showed that herpes simplex type 1 (HSV-1) (odds ratio [OR] = 3.71, p = 8.03 × 10-4) and human papillomavirus 71 (HPV-71; OR = 3.56, p = 0.02), were significantly associated with AD after Bonferroni correction. The phylogenetic-related cluster of Herpesviridae was significantly associated with AD in several strata of the data (p < 0.01).

DISCUSSION

Our results support the hypothesis that viral infection, especially HSV-1, is associated with AD risk.

MIAMI-AD (Methylation in Aging and Methylation in AD): an integrative knowledgebase that facilitates explorations of DNA methylation across sex, aging, and Alzheimer's disease

Alzheimer's disease (AD) is a common neurodegenerative disorder with a significant impact on aging populations. DNA methylation (DNAm) alterations have been implicated in both the aging processes and the development of AD. Given that AD affects more women than men, it is also important to explore DNAm changes that occur specifically in each sex. We created MIAMI-AD, a comprehensive knowledge base containing manually curated summary statistics from 97 published tables in 37 studies, all of which included at least 100 participants. MIAMI-AD enables easy browsing, querying, and downloading DNAm associations at multiple levels - at individual CpG, gene, genomic regions, or genome-wide, in one or multiple studies. Moreover, it also offers tools to perform integrative analyses, such as comparing DNAm associations across different phenotypes or tissues, as well as interactive visualizations. Using several use case examples, we demonstrated that MIAMI-AD facilitates our understanding of age-associated CpGs in AD and the sex-specific roles of DNAm in AD. This open-access resource is freely available to the research community, and all the underlying data can be downloaded. MIAMI-AD (https://miami-ad.org/) facilitates integrative explorations to better understand the interplay between DNAm across aging, sex, and AD.

Novel loci for Alzheimer's disease identified by a genome-wide association study in Ashkenazi Jews

INTRODUCTION

Most Alzheimer's disease (AD) loci have been discovered in individuals with European ancestry (EA).

METHODS

We applied principal component analysis using Gaussian mixture models and an Ashkenazi Jewish (AJ) reference genome-wide association study (GWAS) data set to identify Ashkenazi Jews ascertained in GWAS (n = 42,682), whole genome sequencing (WGS, n = 16,815), and whole exome sequencing (WES, n = 20,504) data sets. The association of AD was tested genome wide (GW) in the GWAS and WGS data sets and exome wide (EW) in all three data sets (EW). Gene-based analyses were performed using aggregated rare variants.

RESULTS

In addition to apolipoprotein E (APOE), GW analyses (1355 cases and 1661 controls) revealed associations with TREM2 R47H (p = 9.66 × 10-9 ), rs541586606 near RAB3B (p = 5.01 × 10-8 ), and rs760573036 between SPOCK3 and ANXA10 (p = 6.32 × 10-8 ). In EW analyses (1504 cases and 2047 controls), study-wide significant association was observed with rs1003710 near SMAP2 (p = 1.91 × 10-7 ). A significant gene-based association was identified with GIPR (p = 7.34 × 10-7 ).

DISCUSSION

Our results highlight the efficacy of founder populations for AD genetic studies.

A haptoglobin (HP) structural variant alters the effect of APOE alleles on Alzheimer's disease

BACKGROUND

Haptoglobin (HP) is an antioxidant of apolipoprotein E (APOE), and previous reports have shown HP binds with APOE and amyloid beta (Aβ) to aid its clearance. A common structural variant of the HP gene distinguishes it into two alleles: HP1 and HP2.

METHODS

HP genotypes were imputed in 29 cohorts from the Alzheimer's Disease Genetics Consortium (N = 20,512). Associations between the HP polymorphism and Alzheimer's disease (AD) risk and age of onset through APOE interactions were investigated using regression models.

RESULTS

The HP polymorphism significantly impacts AD risk in European-descent individuals (and in meta-analysis with African-descent individuals) by modifying both the protective effect of APOE ε2 and the detrimental effect of APOE ε4. The effect is particularly significant among APOE ε4 carriers.

DISCUSSION

The effect modification of APOE by HP suggests adjustment and/or stratification by HP genotype is warranted when APOE risk is considered. Our findings also provided directions for further investigations on potential mechanisms behind this association.

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC BeadChip microarrays

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC arrays. We conducted a comprehensive assessment of the EPIC array probe reliability using 138 duplicated blood DNAm samples generated by the Alzheimer's Disease Neuroimaging Initiative study. We introduced a novel statistical measure, the modified intraclass correlation, to better account for the disagreement in duplicate measurements. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliable information for probes on the EPIC array, will serve as a valuable resource for future DNAm studies.

Critical evaluation of the reliability of DNA methylation probes on the Illumina MethylationEPIC BeadChip microarrays

DNA methylation (DNAm) plays a crucial role in a number of complex diseases. However, the reliability of DNAm levels measured using Illumina arrays varies across different probes. Previous research primarily assessed probe reliability by comparing duplicate samples between the 450k-450k or 450k-EPIC platforms, with limited investigations on Illumina EPIC arrays. We conducted a comprehensive assessment of the EPIC array probe reliability using 138 duplicated blood DNAm samples generated by the Alzheimer's Disease Neuroimaging Initiative study. We introduced a novel statistical measure, the modified intraclass correlation, to better account for the disagreement in duplicate measurements. We observed higher reliability in probes with average methylation beta values of 0.2 to 0.8, and lower reliability in type I probes or those within the promoter and CpG island regions. Importantly, we found that probe reliability has significant implications in the analyses of Epigenome-wide Association Studies (EWAS). Higher reliability is associated with more consistent effect sizes in different studies, the identification of differentially methylated regions (DMRs) and methylation quantitative trait locus (mQTLs), and significant correlations with downstream gene expression. Moreover, blood DNAm measurements obtained from probes with higher reliability are more likely to show concordance with brain DNAm measurements. Our findings, which provide crucial reliable information for probes on the EPIC array, will serve as a valuable resource for future DNAm studies.

Associations of Sex, Race, and Apolipoprotein E Alleles With Multiple Domains of Cognition Among Older Adults

Importance

Sex differences are established in associations between apolipoprotein E (APOE) ε4 and cognitive impairment in Alzheimer disease (AD). However, it is unclear whether sex-specific cognitive consequences of APOE are consistent across races and extend to the APOE ε2 allele.

Objective

To investigate whether sex and race modify APOE ε4 and ε2 associations with cognition.

Design, Setting, and Participants

This genetic association study included longitudinal cognitive data from 4 AD and cognitive aging cohorts. Participants were older than 60 years and self-identified as non-Hispanic White or non-Hispanic Black (hereafter, White and Black). Data were previously collected across multiple US locations from 1994 to 2018. Secondary analyses began December 2021 and ended September 2022.

Main Outcomes and Measures

Harmonized composite scores for memory, executive function, and language were generated using psychometric approaches. Linear regression assessed interactions between APOE ε4 or APOE ε2 and sex on baseline cognitive scores, while linear mixed-effect models assessed interactions on cognitive trajectories. The intersectional effect of race was modeled using an APOE × sex × race interaction term, assessing whether APOE × sex interactions differed by race. Models were adjusted for age at baseline and corrected for multiple comparisons.

Results

Of 32 427 participants who met inclusion criteria, there were 19 007 females (59%), 4453 Black individuals (14%), and 27 974 White individuals (86%); the mean (SD) age at baseline was 74 years (7.9). At baseline, 6048 individuals (19%) had AD, 4398 (14%) were APOE ε2 carriers, and 12 538 (38%) were APOE ε4 carriers. Participants missing APOE status were excluded (n = 9266). For APOE ε4, a robust sex interaction was observed on baseline memory (β = -0.071, SE = 0.014; P = 9.6 × 10-7), whereby the APOE ε4 negative effect was stronger in females compared with males and did not significantly differ among races. Contrastingly, despite the large sample size, no APOE ε2 × sex interactions on cognition were observed among all participants. When testing for intersectional effects of sex, APOE ε2, and race, an interaction was revealed on baseline executive function among individuals who were cognitively unimpaired (β = -0.165, SE = 0.066; P = .01), whereby the APOE ε2 protective effect was female-specific among White individuals but male-specific among Black individuals.

Conclusions and Relevance

In this study, while race did not modify sex differences in APOE ε4, the APOE ε2 protective effect could vary by race and sex. Although female sex enhanced ε4-associated risk, there was no comparable sex difference in ε2, suggesting biological pathways underlying ε4-associated risk are distinct from ε2 and likely intersect with age-related changes in sex biology.

Extended genome-wide association study employing the African Genome Resources Panel identifies novel susceptibility loci for Alzheimer's Disease in individuals of African ancestry

INTRODUCTION

Despite a two-fold increased risk, individuals of African ancestry have been significantly underrepresented in Alzheimer's Disease (AD) genomics efforts.

METHODS

GWAS of 2,903 AD cases and 6,265 cognitive controls of African ancestry. Within-dataset results were meta-analyzed, followed by gene-based and pathway analyses, and analysis of RNAseq and whole-genome sequencing data.

RESULTS

A novel AD risk locus was identified in MPDZ on chromosome 9p23 (rs141610415, MAF=.002, P =3.68×10 -9 ). Two additional novel common and nine novel rare loci approached genome-wide significance at P <9×10 -7 . Comparison of association and LD patterns between datasets with higher and lower degrees of African ancestry showed differential association patterns at chr12q23.2 ( ASCL1 ), suggesting that the association is modulated by regional origin of local African ancestry.

DISCUSSION

Increased sample sizes and sample sets from Africa covering as much African genetic diversity as possible will be critical to identify additional disease-associated loci and improve deconvolution of local genetic ancestry effects.

Identification of novel genes for age-at-onset of Alzheimer's disease by combining quantitative and survival trait analyses

INTRODUCTION

Our understanding of the genetic predisposition for age-at-onset (AAO) of Alzheimer's disease (AD) is limited. Here, we sought to identify genes modifying AAO and examined whether any have sex-specific effects.

METHODS

Genome-wide association analysis were performed on imputed genetic data of 9219 AD cases and 10,345 controls from 20 cohorts of the Alzheimer's Disease Genetics Consortium. AAO was modeled from cases directly and as a survival outcome.

RESULTS

We identified 11 genome-wide significant loci (P < 5 × 10-8 ), including six known AD-risk genes and five novel loci, UMAD1, LUZP2, ARFGEF2, DSCAM, and 4q25, affecting AAO of AD. Additionally, 39 suggestive loci showed strong association. Twelve loci showed sex-specific effects on AAO including CD300LG and MLX/TUBG2 for females and MIR4445 for males.

DISCUSSION

Genes that influence AAO of AD are excellent therapeutic targets for delaying onset of AD. Several loci identified include genes with promising functional implications for AD.

Admixture mapping identifies novel Alzheimer's disease risk regions in African Americans

BACKGROUND

This study used admixture mapping to prioritize the genetic regions associated with Alzheimer's disease (AD) in African American (AA) individuals, followed by ancestry-aware regression analysis to fine-map the prioritized regions.

METHODS

We analyzed 10,271 individuals from 17 different AA datasets. We performed admixture mapping and meta-analyzed the results. We then used regression analysis, adjusting for local ancestry main effects and interactions with genotype, to refine the regions identified from admixture mapping. Finally, we leveraged in silico annotation and differential gene expression data to prioritize AD-related variants and genes.

RESULTS

Admixture mapping identified two genome-wide significant loci on chromosomes 17p13.2 (p = 2.2 × 10-5 ) and 18q21.33 (p = 1.2 × 10-5 ). Our fine mapping of the chromosome 17p13.2 and 18q21.33 regions revealed several interesting genes such as the MINK1, KIF1C, and BCL2.

DISCUSSION

Our ancestry-aware regression approach showed that AA individuals have a lower risk of AD if they inherited African ancestry admixture block at the 17p13.2 locus.

HIGHLIGHTS

We identified two genome-wide significant admixture mapping signals: on chromosomes 17p13.2 and 18q21.33, which are novel in African American (AA) populations. Our ancestry-aware regression approach showed that AA individuals have a lower risk of Alzheimer's disease (AD) if they inherited African ancestry admixture block at the 17p13.2 locus. We found that the overall proportion of African ancestry does not differ between the cases and controls that suggest African genetic ancestry alone is not likely to explain the AD prevalence difference between AA and non-Hispanic White populations.

Distinct CSF biomarker-associated DNA methylation in Alzheimer's disease and cognitively normal subjects

BACKGROUND

Growing evidence has demonstrated that DNA methylation (DNAm) plays an important role in Alzheimer's disease (AD) and that DNAm differences can be detected in the blood of AD subjects. Most studies have correlated blood DNAm with the clinical diagnosis of AD in living individuals. However, as the pathophysiological process of AD can begin many years before the onset of clinical symptoms, there is often disagreement between neuropathology in the brain and clinical phenotypes. Therefore, blood DNAm associated with AD neuropathology, rather than with clinical data, would provide more relevant information on AD pathogenesis.

METHODS

We performed a comprehensive analysis to identify blood DNAm associated with cerebrospinal fluid (CSF) pathological biomarkers for AD. Our study included matched samples of whole blood DNA methylation, CSF Aβ₄₂, phosphorylated tau₁₈₁ (pTau₁₈₁), and total tau (tTau) biomarkers data, measured on the same subjects and at the same clinical visits from a total of 202 subjects (123 CN or cognitively normal, 79 AD) in the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. To validate our findings, we also examined the association between premortem blood DNAm and postmortem brain neuropathology measured on a group of 69 subjects in the London dataset.

RESULTS

We identified a number of novel associations between blood DNAm and CSF biomarkers, demonstrating that changes in pathological processes in the CSF are reflected in the blood epigenome. Overall, the CSF biomarker-associated DNAm is relatively distinct in CN and AD subjects, highlighting the importance of analyzing omics data measured on cognitively normal subjects (which includes preclinical AD subjects) to identify diagnostic biomarkers, and considering disease stages in the development and testing of AD treatment strategies. Moreover, our analysis revealed biological processes associated with early brain impairment relevant to AD are marked by DNAm in the blood, and blood DNAm at several CpGs in the DMR on HOXA5 gene are associated with pTau₁₈₁ in the CSF, as well as tau-pathology and DNAm in the brain, nominating DNAm at this locus as a promising candidate AD biomarker.

CONCLUSIONS

Our study provides a valuable resource for future mechanistic and biomarker studies of DNAm in AD.

Distinct CSF biomarker-associated DNA methylation in Alzheimer's disease and cognitively normal subjects

Background Growing evidence has demonstrated that DNA methylation (DNAm) plays an important role in Alzheimer's disease (AD) and that DNAm differences can be detected in the blood of AD subjects. Most studies have correlated blood DNAm with the clinical diagnosis of AD in living individuals. However, as the pathophysiological process of AD can begin many years before the onset of clinical symptoms, there is often disagreement between neuropathology in the brain and clinical phenotypes. Therefore, blood DNAm associated with AD neuropathology, rather than with clinical data, would provide more relevant information on AD pathogenesis. Methods We performed a comprehensive analysis to identify blood DNAm associated with cerebrospinal fluid (CSF) pathological biomarkers for AD. Our study included matched samples of whole blood DNA methylation, CSF Aβ ₄₂ , phosphorylated tau ₁₈₁ (pTau ₁₈₁ ), and total tau (tTau) biomarkers data, measured on the same subjects and at the same clinical visits from a total of 202 subjects (123 CN or cognitively normal, 79 AD) in the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. To validate our findings, we also examined the association between premortem blood DNAm and postmortem brain neuropathology measured on a group of 69 subjects in the London dataset. Results We identified a number of novel associations between blood DNAm and CSF biomarkers, demonstrating that changes in pathological processes in the CSF are reflected in the blood epigenome. Overall, the CSF biomarker-associated DNAm is relatively distinct in CN and AD subjects, highlighting the importance of analyzing omics data measured on cognitively normal subjects (which includes preclinical AD subjects) to identify diagnostic biomarkers, and considering disease stages in the development and testing of AD treatment strategies. Moreover, our analysis revealed biological processes associated with early brain impairment relevant to AD are marked by DNAm in the blood, and blood DNAm at several CpGs in the DMR on HOXA5 gene are associated with pTau ₁₈₁ in the CSF, as well as tau-pathology and DNAm in the brain, nominating DNAm at this locus as a promising candidate AD biomarker. Conclusions Our study provides a valuable resource for future mechanistic and biomarker studies of DNAm in AD.

Explainable machine learning aggregates polygenic risk scores and electronic health records for Alzheimer's disease prediction

Alzheimer's disease (AD) is the most common late-onset neurodegenerative disorder. Identifying individuals at increased risk of developing AD is important for early intervention. Using data from the Alzheimer Disease Genetics Consortium, we constructed polygenic risk scores (PRSs) for AD and age-at-onset (AAO) of AD for the UK Biobank participants. We then built machine learning (ML) models for predicting development of AD, and explored feature importance among PRSs, conventional risk factors, and ICD-10 codes from electronic health records, a total of > 11,000 features using the UK Biobank dataset. We used eXtreme Gradient Boosting (XGBoost) and SHapley Additive exPlanations (SHAP), which provided superior ML performance as well as aided ML model explanation. For participants age 40 and older, the area under the curve for AD was 0.88. For subjects of age 65 and older (late-onset AD), PRSs were the most important predictors. This is the first observation that PRSs constructed from the AD risk and AAO play more important roles than age in predicting AD. The ML model also identified important predictors from EHR, including urinary tract infection, syncope and collapse, chest pain, disorientation and hypercholesterolemia, for developing AD. Our ML model improved the accuracy of AD risk prediction by efficiently exploring numerous predictors and identified novel feature patterns.

Neuropathological lesions and their contribution to dementia and cognitive impairment in a heterogeneous clinical population

INTRODUCTION

Alzheimer disease (AD) and related dementias are characterized by damage caused by neuropathological lesions in the brain. These include AD lesions (plaques and tangles) and non-AD lesions such as vascular injury or Lewy bodies. We report here an assessment of lesion association to dementia in a large clinic-based population.

METHODS

We identified 5272 individuals with neuropathological data from the National Alzheimer's Coordinating Center. Individual lesions, as well as a neuropathological composite score (NPCS) were tested for association with dementia, and both functional and neurocognitive impairment using regression models.

RESULTS

Most individuals exhibited mixed pathologies, especially AD lesions in combination with non-AD lesions. All lesion types were associated with one or more clinical outcomes; most even while controlling for AD pathology. The NPCS was also associated with clinical outcomes.

DISCUSSION

These data suggest mixed-type pathologies are extremely common in a clinic-based population and may contribute to dementia and cognitive impairment.

Exome sequencing identifies rare damaging variants in ATP8B4 and ABCA1 as risk factors for Alzheimer's disease

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.

Distinct sex-specific DNA methylation differences in Alzheimer's disease

BACKGROUND

Sex is increasingly recognized as a significant factor contributing to the biological and clinical heterogeneity in AD. There is also growing evidence for the prominent role of DNA methylation (DNAm) in Alzheimer's disease (AD).

METHODS

We studied sex-specific DNA methylation differences in the blood samples of AD subjects compared to cognitively normal subjects, by performing sex-specific meta-analyses of two large blood-based epigenome-wide association studies (ADNI and AIBL), which included DNA methylation data for a total of 1284 whole blood samples (632 females and 652 males). Within each dataset, we used two complementary analytical strategies, a sex-stratified analysis that examined methylation to AD associations in male and female samples separately, and a methylation-by-sex interaction analysis that compared the magnitude of these associations between different sexes. After adjusting for age, estimated immune cell type proportions, batch effects, and correcting for inflation, the inverse-variance fixed-effects meta-analysis model was used to identify the most consistent DNAm differences across datasets. In addition, we also evaluated the performance of the sex-specific methylation-based risk prediction models for AD diagnosis using an independent external dataset.

RESULTS

In the sex-stratified analysis, we identified 2 CpGs, mapped to the PRRC2A and RPS8 genes, significantly associated with AD in females at a 5% false discovery rate, and an additional 25 significant CpGs (21 in females, 4 in males) at P-value < 1×10^-5. In methylation-by-sex interaction analysis, we identified 5 significant CpGs at P-value < 10^-5. Out-of-sample validations using the AddNeuroMed dataset showed in females, the best logistic prediction model included age, estimated immune cell-type proportions, and methylation risk scores (MRS) computed from 9 of the 23 CpGs identified in AD vs. CN analysis that are also available in AddNeuroMed dataset (AUC = 0.74, 95% CI: 0.65-0.83). In males, the best logistic prediction model included only age and MRS computed from 2 of the 5 CpGs identified in methylation-by-sex interaction analysis that are also available in the AddNeuroMed dataset (AUC = 0.70, 95% CI: 0.56-0.82).

CONCLUSIONS

Overall, our results show that the DNA methylation differences in AD are largely distinct between males and females. Our best-performing sex-specific methylation-based prediction model in females performed better than that for males and additionally included estimated cell-type proportions. The significant discriminatory classification of AD samples with our methylation-based prediction models demonstrates that sex-specific DNA methylation could be a predictive biomarker for AD. As sex is a strong factor underlying phenotypic variability in AD, the results of our study are particularly relevant for a better understanding of the epigenetic architecture that underlie AD and for promoting precision medicine in AD.

Genetic architecture of RNA editing regulation in Alzheimer's disease across diverse ancestral populations

Most Alzheimer's disease (AD)-associated genetic variants do not change protein coding sequence and thus likely exert their effects through regulatory mechanisms. RNA editing, the post-transcriptional modification of RNA bases, is a regulatory feature that is altered in AD patients that differs across ancestral backgrounds. Editing QTLs (edQTLs) are DNA variants that influence the level of RNA editing at a specific site. To study the relationship of DNA variants genome-wide, and particularly in AD-associated loci, with RNA editing, we performed edQTL analyses in self-reported individuals of African American (AF) or White (EU) race with corresponding global genetic ancestry averaging 82.2% African ancestry (AF) and 96.8% European global ancestry (EU) in the two groups, respectively. We used whole-genome genotyping array and RNA sequencing data from peripheral blood of 216 AD cases and 212 age-matched, cognitively intact controls. We identified 2144 edQTLs in AF and 3579 in EU, of which 1236 were found in both groups. Among these, edQTLs in linkage disequilibrium (r2 > 0.5) with AD-associated genetic variants in the SORL1, SPI1 and HLA-DRB1 loci were associated with sites that were differentially edited between AD cases and controls. While there is some shared RNA editing regulatory architecture, most edQTLs had distinct effects on the rate of RNA editing in different ancestral populations suggesting a complex architecture of RNA editing regulation. Altered RNA editing may be one possible mechanism for the functional effect of AD-associated variants and may contribute to observed differences in the genetic etiology of AD between ancestries.

Cross-tissue analysis of blood and brain epigenome-wide association studies in Alzheimer's disease

To better understand DNA methylation in Alzheimer's disease (AD) from both mechanistic and biomarker perspectives, we performed an epigenome-wide meta-analysis of blood DNA methylation in two large independent blood-based studies in AD, the ADNI and AIBL studies, and identified 5 CpGs, mapped to the SPIDR, CDH6 genes, and intergenic regions, that are significantly associated with AD diagnosis. A cross-tissue analysis that combined these blood DNA methylation datasets with four brain methylation datasets prioritized 97 CpGs and 10 genomic regions that are significantly associated with both AD neuropathology and AD diagnosis. An out-of-sample validation using the AddNeuroMed dataset showed the best performing logistic regression model includes age, sex, immune cell type proportions, and methylation risk score based on prioritized CpGs in cross-tissue analysis (AUC = 0.696, 95% CI: 0.616 - 0.770, P-value = 2.78 × 10^-5). Our study offers new insights into epigenetics in AD and provides a valuable resource for future AD biomarker discovery.

Genome-wide association and multi-omics studies identify MGMT as a novel risk gene for Alzheimer's disease among women

INTRODUCTION

Variants in the tau gene (MAPT) region are associated with breast cancer in women and Alzheimer's disease (AD) among persons lacking apolipoprotein E ε4 (ε4-).

METHODS

To identify novel genes associated with tau-related pathology, we conducted two genome-wide association studies (GWAS) for AD, one among 10,340 ε4- women in the Alzheimer's Disease Genetics Consortium (ADGC) and another in 31 members (22 women) of a consanguineous Hutterite kindred.

RESULTS

We identified novel associations of AD with MGMT variants in the ADGC (rs12775171, odds ratio [OR] = 1.4, P = 4.9 × 10-8 ) and Hutterite (rs12256016 and rs2803456, OR = 2.0, P = 1.9 × 10-14 ) datasets. Multi-omics analyses showed that the most significant and largest number of associations among the single nucleotide polymorphisms (SNPs), DNA-methylated CpGs, MGMT expression, and AD-related neuropathological traits were observed among women. Furthermore, promoter capture Hi-C analyses revealed long-range interactions of the MGMT promoter with MGMT SNPs and CpG sites.

DISCUSSION

These findings suggest that epigenetically regulated MGMT expression is involved in AD pathogenesis, especially in women.

MethReg: estimating the regulatory potential of DNA methylation in gene transcription

Epigenome-wide association studies often detect many differentially methylated sites, and many are located in distal regulatory regions. To further prioritize these significant sites, there is a critical need to better understand the functional impact of CpG methylation. Recent studies demonstrated that CpG methylation-dependent transcriptional regulation is a widespread phenomenon. Here, we present MethReg, an R/Bioconductor package that analyzes matched DNA methylation and gene expression data, along with external transcription factor (TF) binding information, to evaluate, prioritize and annotate CpG sites with high regulatory potential. At these CpG sites, TF-target gene associations are often only present in a subset of samples with high (or low) methylation levels, so they can be missed by analyses that use all samples. Using colorectal cancer and Alzheimer's disease datasets, we show MethReg significantly enhances our understanding of the regulatory roles of DNA methylation in complex diseases.

Manifestations of Alzheimer's disease genetic risk in the blood are evident in a multiomic analysis in healthy adults aged 18 to 90

Genetics play an important role in late-onset Alzheimer's Disease (AD) etiology and dozens of genetic variants have been implicated in AD risk through large-scale GWAS meta-analyses. However, the precise mechanistic effects of most of these variants have yet to be determined. Deeply phenotyped cohort data can reveal physiological changes associated with genetic risk for AD across an age spectrum that may provide clues to the biology of the disease. We utilized over 2000 high-quality quantitative measurements obtained from blood of 2831 cognitively normal adult clients of a consumer-based scientific wellness company, each with CLIA-certified whole-genome sequencing data. Measurements included: clinical laboratory blood tests, targeted chip-based proteomics, and metabolomics. We performed a phenome-wide association study utilizing this diverse blood marker data and 25 known AD genetic variants and an AD-specific polygenic risk score (PGRS), adjusting for sex, age, vendor (for clinical labs), and the first four genetic principal components; sex-SNP interactions were also assessed. We observed statistically significant SNP-analyte associations for five genetic variants after correction for multiple testing (for SNPs in or near NYAP1, ABCA7, INPP5D, and APOE), with effects detectable from early adulthood. The ABCA7 SNP and the APOE2 and APOE4 encoding alleles were associated with lipid variability, as seen in previous studies; in addition, six novel proteins were associated with the e2 allele. The most statistically significant finding was between the NYAP1 variant and PILRA and PILRB protein levels, supporting previous functional genomic studies in the identification of a putative causal variant within the PILRA gene. We did not observe associations between the PGRS and any analyte. Sex modified the effects of four genetic variants, with multiple interrelated immune-modulating effects associated with the PICALM variant. In post-hoc analysis, sex-stratified GWAS results from an independent AD case-control meta-analysis supported sex-specific disease effects of the PICALM variant, highlighting the importance of sex as a biological variable. Known AD genetic variation influenced lipid metabolism and immune response systems in a population of non-AD individuals, with associations observed from early adulthood onward. Further research is needed to determine whether and how these effects are implicated in early-stage biological pathways to AD. These analyses aim to complement ongoing work on the functional interpretation of AD-associated genetic variants.

New insights into the genetic etiology of Alzheimer's disease and related dementias

Characterization of the genetic landscape of Alzheimer's disease (AD) and related dementias (ADD) provides a unique opportunity for a better understanding of the associated pathophysiological processes. We performed a two-stage genome-wide association study totaling 111,326 clinically diagnosed/'proxy' AD cases and 677,663 controls. We found 75 risk loci, of which 42 were new at the time of analysis. Pathway enrichment analyses confirmed the involvement of amyloid/tau pathways and highlighted microglia implication. Gene prioritization in the new loci identified 31 genes that were suggestive of new genetically associated processes, including the tumor necrosis factor alpha pathway through the linear ubiquitin chain assembly complex. We also built a new genetic risk score associated with the risk of future AD/dementia or progression from mild cognitive impairment to AD/dementia. The improvement in prediction led to a 1.6- to 1.9-fold increase in AD risk from the lowest to the highest decile, in addition to effects of age and the APOE ε4 allele.

An association test of the spatial distribution of rare missense variants within protein structures identifies Alzheimer's disease-related patterns

More than 90% of genetic variants are rare in most modern sequencing studies, such as the Alzheimer's Disease Sequencing Project (ADSP) whole-exome sequencing (WES) data. Furthermore, 54% of the rare variants in ADSP WES are singletons. However, both single variant and unit-based tests are limited in their statistical power to detect an association between rare variants and phenotypes. To best use missense rare variants and investigate their biological effect, we examine their association with phenotypes in the context of protein structures. We developed a protein structure-based approach, protein optimized kernel evaluation of missense nucleotides (POKEMON), which evaluates rare missense variants based on their spatial distribution within a protein rather than their allele frequency. The hypothesis behind this test is that the three-dimensional spatial distribution of variants within a protein structure provides functional context to power an association test. POKEMON identified three candidate genes (TREM2, SORL1, and EXOC3L4) and another suggestive gene from the ADSP WES data. For TREM2 and SORL1, two known Alzheimer's disease (AD) genes, the signal from the spatial cluster is stable even if we exclude known AD risk variants, indicating the presence of additional low-frequency risk variants within these genes. EXOC3L4 is a novel AD risk gene that has a cluster of variants primarily shared by case subjects around the Sec6 domain. This cluster is also validated in an independent replication data set and a validation data set with a larger sample size.

Genome-Wide Linkage Study Meta-Analysis of Male Sexual Orientation

Male sexual orientation is a scientifically and socially important trait shown by family and twin studies to be influenced by environmental and complex genetic factors. Individual genome-wide linkage studies (GWLS) have been conducted, but not jointly analyzed. Two main datasets account for > 90% of the published GWLS concordant sibling pairs on the trait and are jointly analyzed here: MGSOSO (Molecular Genetic Study of Sexual Orientation; 409 concordant sibling pairs in 384 families, Sanders et al. (2015)) and Hamer (155 concordant sibling pairs in 145 families, Mustanski et al. (2005)). We conducted multipoint linkage analyses with Merlin on the datasets separately since they were genotyped differently, integrated genetic marker positions, and combined the resultant LOD (logarithm of the odds) scores at each 1 cM grid position. We continue to find the strongest linkage support at pericentromeric chromosome 8 and chromosome Xq28. We also incorporated the remaining published GWLS dataset (on 55 families) by using meta-analytic approaches on published summary statistics. The meta-analysis has maximized the positional information from GWLS of currently available family resources and can help prioritize findings from genome-wide association studies (GWAS) and other approaches. Although increasing evidence highlights genetic contributions to male sexual orientation, our current understanding of contributory loci is still limited, consistent with the complexity of the trait. Further increasing genetic knowledge about male sexual orientation, especially via large GWAS, should help advance our understanding of the biology of this important trait.

Genome-Wide Linkage and Association Study of Childhood Gender Nonconformity in Males

Male sexual orientation is influenced by environmental and complex genetic factors. Childhood gender nonconformity (CGN) is one of the strongest correlates of homosexuality with substantial familiality. We studied brothers in families with two or more homosexual brothers (409 concordant sibling pairs in 384 families, as well as their heterosexual brothers), who self-recalled their CGN. To map loci for CGN, we conducted a genome-wide linkage scan (GWLS) using SNP genotypes. The strongest linkage peaks, each with significant or suggestive two-point LOD scores and multipoint LOD score support, were on chromosomes 5q31 (maximum two-point LOD = 4.45), 6q12 (maximum two-point LOD = 3.64), 7q33 (maximum two-point LOD = 3.09), and 8q24 (maximum two-point LOD = 3.67), with the latter not overlapping with previously reported strongest linkage region for male sexual orientation on pericentromeric chromosome 8. Family-based association analyses were used to identify associated variants in the linkage regions, with a cluster of SNPs (minimum association p = 1.3 × 10^-8) found at the 5q31 linkage peak. Genome-wide, clusters of multiple SNPs in the 10^-6 to 10^-8 p-value range were found at chromosomes 5p13, 5q31, 7q32, 8p22, and 10q23, highlighting glutamate-related genes. This is the first reported GWLS and genome-wide association study on CGN. Further increasing genetic knowledge about CGN and its relationships to male sexual orientation should help advance our understanding of the biology of these associated traits.

Genomic evidence consistent with antagonistic pleiotropy may help explain the evolutionary maintenance of same-sex sexual behaviour in humans

Human same-sex sexual behaviour (SSB) is heritable, confers no immediately obvious direct reproductive or survival benefit and can divert mating effort from reproductive opportunities. This presents a Darwinian paradox: why has SSB been maintained despite apparent selection against it? We show that genetic effects associated with SSB may, in individuals who only engage in opposite-sex sexual behaviour (OSB individuals), confer a mating advantage. Using results from a recent genome-wide association study of SSB and a new genome-wide association study on number of opposite-sex sexual partners in 358,426 individuals, we show that, among OSB individuals, genetic effects associated with SSB are associated with having more opposite-sex sexual partners. Computer simulations suggest that such a mating advantage for alleles associated with SSB could help explain how it has been evolutionarily maintained. Caveats include the cultural specificity of our UK and US samples, the societal regulation of sexual behaviour in these populations, the difficulty of measuring mating success and the fact that measured variants capture a minority of the total genetic variation in the traits.

Genomewide Association Studies of LRRK2 Modifiers of Parkinson's Disease

Objective

The aim of this study was to search for genes/variants that modify the effect of LRRK2 mutations in terms of penetrance and age‐at‐onset of Parkinson's disease.

Methods

We performed the first genomewide association study of penetrance and age‐at‐onset of Parkinson's disease in LRRK2 mutation carriers (776 cases and 1,103 non‐cases at their last evaluation). Cox proportional hazard models and linear mixed models were used to identify modifiers of penetrance and age‐at‐onset of LRRK2 mutations, respectively. We also investigated whether a polygenic risk score derived from a published genomewide association study of Parkinson's disease was able to explain variability in penetrance and age‐at‐onset in LRRK2 mutation carriers.

Results

A variant located in the intronic region of CORO1C on chromosome 12 (rs77395454; p value = 2.5E‐08, beta = 1.27, SE = 0.23, risk allele: C) met genomewide significance for the penetrance model. Co‐immunoprecipitation analyses of LRRK2 and CORO1C supported an interaction between these 2 proteins. A region on chromosome 3, within a previously reported linkage peak for Parkinson's disease susceptibility, showed suggestive associations in both models (penetrance top variant: p value = 1.1E‐07; age‐at‐onset top variant: p value = 9.3E‐07). A polygenic risk score derived from publicly available Parkinson's disease summary statistics was a significant predictor of penetrance, but not of age‐at‐onset.

Interpretation

This study suggests that variants within or near CORO1C may modify the penetrance of LRRK2 mutations. In addition, common Parkinson's disease associated variants collectively increase the penetrance of LRRK2 mutations. ANN NEUROL 2021;90:82–94

Sex-specific DNA methylation differences in Alzheimer's disease pathology

Sex is an important factor that contributes to the clinical and biological heterogeneities in Alzheimer's disease (AD), but the regulatory mechanisms underlying sex disparity in AD are still not well understood. DNA methylation is an important epigenetic modification that regulates gene transcription and is known to be involved in AD. We performed the first large-scale sex-specific meta-analysis of DNA methylation differences in AD neuropathology, by re-analyzing four recent epigenome-wide association studies totaling more than 1000 postmortem prefrontal cortex brain samples using a uniform analytical pipeline. For each cohort, we employed two complementary analytical strategies, a sex-stratified analysis that examined methylation-Braak stage associations in male and female samples separately, and a sex-by-Braak stage interaction analysis that compared the magnitude of these associations between different sexes. Our analysis uncovered 14 novel CpGs, mapped to genes such as TMEM39A and TNXB that are associated with the AD Braak stage in a sex-specific manner. TMEM39A is known to be involved in inflammation, dysregulated type I interferon responses, and other immune processes. TNXB encodes tenascin proteins, which are extracellular matrix glycoproteins demonstrated to modulate synaptic plasticity in the brain. Moreover, for many previously implicated genes in AD neuropathology, such as MBP and AZU1, our analysis provided the new insights that they were predominately driven by effects in only one sex. These sex-specific DNA methylation differences were enriched in divergent biological processes such as integrin activation in females and complement activation in males. Our study implicated multiple new loci and biological processes that affected AD neuropathology in a sex-specific manner.