Gene‐based polygenic score analysis identifies novel immune genes with deleterious variants that associate with Alzheimer’s disease

Abstract

AbstractBackgroundIn the Alzheimer’s disease sequencing project (ADSP) whole exome dataset, only the well‐known TREM2.pR47H variant shows genome‐wide significance by conventional single variant analysis. To identify novel AD genes related to TREM2, we used gene‐based polygenic score analysis (PGSA) to evaluate an immune, temporal cortical co‐expression network with 229 genes including TREM2.MethodThe most recent whole exome variant call files from the ADSP (NG00067.v7) were subjected to stringent quality control. Post‐QC samples were divided into discovery (4174 samples sequenced at Broad), test (3299 samples sequenced at WashU), and follow‐up (2273 samples sequenced at Baylor) sets. The pruned variants (r2 < 0.2) with a minor allele count over 20 in each gene were analyzed by Broad/WashU PGSA thereby obtaining P‐values (gPbr.wa) for each gene in WashU samples. To assess significance in independent Baylor samples, cumulative Broad/Baylor PGSA was performed using a novel approach in which variants were stratified by both gene (gPbr.wa) and variant (Broad.P).ResultAfter excluding TREM2, the five variants with gPbr.wa ≤ 0.08 and Broad.P ≤ 0.03 showed significant association in Broad/Baylor PGSA. One of these variants showed no evidence of association with AD in the combined dataset. Polygenic scores for the other four variants, which were in three genes, showed significant association in Baylor samples, and significant, essentially identical association in WashU samples. Polygenic scores for nine additional variants in these genes (0.03> Broad.P ≤ 0.13) also showed significant, essentially identical association in Baylor and WashU samples. The association of these 13 variants with AD was driven primarily by six deleterious, protein‐altering variants with minor allele frequencies less than 1%, but the other 7 variants also contributed.ConclusionUsing gene‐based polygenic score analysis to analyze an immune co‐expression network, we identify three novel immune genes with actionable nonsynonymous and stopgain variants that show significant, replicable association with AD in the most recent release of the ADSP whole exome dataset.