Identification of gene-gene interactions for Alzheimer's disease using co-operative game theory

Thesis (Ph.D.)--Boston UniversityThe multifactorial nature of Alzheimer's Disease suggests that complex gene-gene interactions are present in AD pathways. Contemporary approaches to detect such interactions in genome-wide data are mathematically and computationally challenging. We investigated gene-gene interactions for AD using a novel algorithm based on cooperative game theory in 15 genome-wide association study (GWAS) datasets comprising of a total of 11,840 AD cases and 10,931 cognitively normal elderly controls from the Alzheimer Disease Genetics Consortium (ADGC). We adapted this approach, which was developed originally for solving multi-dimensional problems in economics and social sciences, to compute a Shapely value statistic to identify genetic markers that contribute most to coalitions of SNPs in predicting AD risk. Treating each GWAS dataset as independent discovery, markers were ranked according to their contribution to coalitions formed with other markers. Using a backward elimination strategy, markers with low Shapley values were eliminated and the statistic was recalculated iteratively. We tested all two-way interactions between top Shapley markers in regression models which included the two SNPs (main effects) and a term for their interaction. Models yielding a p-value<0.05 for the interaction term were evaluated in each of the other datasets and the results from all datasets were combined by meta-analysis. Statistically significant interactions were observed with multiple marker combinations in the APOE regions. My analyses also revealed statistically strong interactions between markers in 6 regions; CTNNA3-ATP11A (p=4.1E-07), CSMD1-PRKCQ (p=3.5E-08), DCC-UNC5CL (p=5.9e-8), CNTNAP2-RFC3 (p=1.16e-07), AACS-TSHZ3 (p=2.64e-07) and CAMK4-MMD (p=3.3e-07). The Shapley value algorithm outperformed Chi-Square and ReliefF in detecting known interactions between APOE and GAB2 in a previously published GWAS dataset. It was also more accurate than competing filtering methods in identifying simulated epistastic SNPs that are additive in nature, but its accuracy was low in identifying non-linear interactions. The game theory algorithm revealed strong interactions between markers in novel genes with weak main effects, which would have been overlooked if only markers with strong marginal association with AD were tested. This method will be a valuable tool for identifying gene-gene interactions for complex diseases and other traits

The APOE ε4 allele is associated with a reduction in FEV1/FVC in women: A cross-sectional analysis of the Long Life Family Study

Introduction Murine studies have shown that apolipoprotein E modulates pulmonary function during development, aging, and allergen-induced airway disease. It is not known whether the polymorphic human APOE gene influences pulmonary function. Objectives We assessed whether an association exists between the polymorphic human APOE ε2, ε3, and ε4 alleles and pulmonary function among participants in the Long Life Family Study. Methods Data from 4,468 Caucasian subjects who had genotyping performed for the APOE ε2, ε3, and ε4 alleles were analyzed, with and without stratification by sex. Statistical models were fitted considering the effects of the ε2 allele, defined as ε2/2 or ε2/3 genotypes, and the ε4 allele, defined as ε3/4 or ε4/4 genotypes, which were compared to the ε3/3 genotype. Results The mean FEV1/FVC ratio (the forced expiratory volume in one second divided by the forced vital capacity) was lower among women with the ε4 allele as compared to women with the ε3/3 genotype or the ε2 allele. Carriage of the APOE ε4 allele was associated with FEV1/FVC, which implied lower values. Further analysis showed that the association primarily reflected women without lung disease who were older than 70 years. The association was not mediated by lipid levels, smoking status, body mass index, or cardiovascular disease. Conclusions This study for the first time identifies that the APOE gene is associated with modified lung physiology in women. This suggests that a link may exist between the APOE ε4 allele, female sex, and a reduction in the FEV1/FVC ratio in older individuals

Association of Rare Coding Mutations With Alzheimer Disease and Other Dementias Among Adults of European Ancestry

IMPORTANCE Some of the unexplained heritability of Alzheimer disease (AD) may be due to rare variants whose effects are not captured in genome-wide association studies because very large samples are needed to observe statistically significant associations. OBJECTIVE To identify genetic variants associated with AD risk using a nonstatistical approach. DESIGN, SETTING, AND PARTICIPANTS Genetic association study in which rare variants were identified by whole-exome sequencing in unrelated individuals of European ancestry from the Alzheimer’s Disease Sequencing Project (ADSP). Data were analyzed between March 2017 and September 2018. MAIN OUTCOMES AND MEASURES Minor alleles genome-wide and in 95 genes previously associated with AD, AD-related traits, or other dementias were tabulated and filtered for predicted functional impact and occurrence in participants with AD but not controls. Support for several findings was sought in a whole-exome sequencing data set comprising 19 affected relative pairs from Utah high-risk pedigrees and whole-genome sequencing data sets from the ADSP and Alzheimer’s Disease Neuroimaging Initiative. RESULTS Among 5617 participants with AD (3202 [57.0%] women; mean [SD] age, 76.4 [9.3] years) and 4594 controls (2719 [59.0%] women; mean [SD] age, 86.5 [4.5] years), a total of 24 variants with moderate or high functional impact from 19 genes were observed in 10 or more participants with AD but not in controls. These variants included a missense mutation (rs149307620 [p.A284T], n = 10) in NOTCH3, a gene in which coding mutations are associated with cerebral autosomal-dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL), that was also identified in 1 participant with AD and 1 participant with mild cognitive impairment in the whole genome sequencing data sets. Four participants with AD carried the TREM2 rs104894002 (p.Q33X) high-impact mutation that, in homozygous form, causes Nasu-Hakola disease, a rare disorder characterized by early-onset dementia and multifocal bone cysts, suggesting an intermediate inheritance model for the mutation. Compared with controls, participants with AD had a significantly higher burden of deleterious rare coding variants in dementia-associated genes (2314 vs 3354 cumulative variants, respectively; P = .006). CONCLUSIONS AND RELEVANCE Different mutations in the same gene or variable dose of a mutation may be associated with result in distinct dementias. These findings suggest that minor differences in the structure or amount of protein may be associated with in different clinical outcomes. Understanding these genotype-phenotype associations may provide further insight into the pathogenic nature of the mutations, as well as offer clues for developing new therapeutic targets

Ultra-rare mutations in SRCAP segregate in Caribbean Hispanic families with Alzheimer disease.

OBJECTIVE: To identify rare coding variants segregating with late-onset Alzheimer disease (LOAD) in Caribbean Hispanic families. METHODS: Whole-exome sequencing (WES) was completed in 110 individuals from 31 Caribbean Hispanic families without APOE ε4 homozygous carriers. Rare coding mutations segregating in families were subsequently genotyped in additional families and in an independent cohort of Caribbean Hispanic patients and controls. SRCAP messenger RNA (mRNA) expression was assessed in whole blood from mutation carriers with LOAD, noncarriers with LOAD, and healthy elderly controls, and also from autopsied brains in 2 clinical neuropathologic cohort studies of aging and dementia. RESULTS: Ten ultra-rare missense mutations in the Snf2-related CREBBP, activator protein (SRCAP), were found in 12 unrelated families. Compared with the frequency in Caribbean Hispanic controls and the Latino population in the Exome Aggregation Consortium, the frequency of SRCAP mutations among Caribbean Hispanic patients with LOAD was significantly enriched (p = 1.19e-16). mRNA expression of SRCAP in whole blood was significantly lower in mutation carriers with LOAD, while the expression in whole blood and in the brain was significantly higher in nonmutation carriers with LOAD. Brain expression also correlated with clinical and neuropathologic endophenotypes. CONCLUSIONS: WES in Caribbean Hispanic families with LOAD revealed ultra-rare missense mutations in SRCAP, a gene expressed in the brain and mutated in Floating-Harbor syndrome. SRCAP is a potent coactivator of the CREB-binding protein and a regulator of DNA damage response involving ATP-dependent chromatin remodeling. We hypothesize that increased expression in LOAD suggests a compensatory mechanism altered in mutation carriers

Rare Variants Imputation in Admixed Populations: Comparison Across Reference Panels and Bioinformatics Tools

BackgroundImputation has become a standard approach in genome-wide association studies (GWAS) to infer in silico untyped markers. Although feasibility for common variants imputation is well established, we aimed to assess rare and ultra-rare variants’ imputation in an admixed Caribbean Hispanic population (CH).MethodsWe evaluated imputation accuracy in CH (N = 1,000), focusing on rare (0.1% ≤ minor allele frequency (MAF) ≤ 1%) and ultra-rare (MAF &lt; 0.1%) variants. We used two reference panels, the Haplotype Reference Consortium (HRC; N = 27,165) and 1000 Genome Project (1000G phase 3; N = 2,504) and multiple phasing (SHAPEIT, Eagle2) and imputation algorithms (IMPUTE2, MACH-Admix). To assess imputation quality, we reported: (a) high-quality variant counts according to imputation tools’ internal indexes (e.g., IMPUTE2 “Info” ≥ 80%). (b) Wilcoxon Signed-Rank Test comparing imputation quality for genotyped variants that were masked and imputed; (c) Cohen’s kappa coefficient to test agreement between imputed and whole-exome sequencing (WES) variants; (d) imputation of G206A mutation in the PSEN1 (ultra-rare in the general population an more frequent in CH) followed by confirmation genotyping. We also tested ancestry proportion (European, African and Native American) against WES-imputation mismatches in a Poisson regression fashion.ResultsSHAPEIT2 retrieved higher percentage of imputed high-quality variants than Eagle2 (rare: 51.02% vs. 48.60%; ultra-rare 0.66% vs. 0.65%, Wilcoxon p-value &lt; 0.001). SHAPEIT-IMPUTE2 employing HRC outperformed 1000G (64.50% vs. 59.17%; 1.69% vs. 0.75% for high-quality rare and ultra-rare variants, respectively, Wilcoxon p-value &lt; 0.001). SHAPEIT-IMPUTE2 outperformed MaCH-Admix. Compared to 1000G, HRC-imputation retrieved a higher number of high-quality rare and ultra-rare variants, despite showing lower agreement between imputed and WES variants (e.g., rare: 98.86% for HRC vs. 99.02% for 1000G). High Kappa (K = 0.99) was observed for both reference panels. Twelve G206A mutation carriers were imputed and all validated by confirmation genotyping. African ancestry was associated with higher imputation errors for uncommon and rare variants (p-value &lt; 1e-05).ConclusionReference panels with larger numbers of haplotypes can improve imputation quality for rare and ultra-rare variants in admixed populations such as CH. Ethnic composition is an important predictor of imputation accuracy, with higher African ancestry associated with poorer imputation accuracy

SORL1 mutations in early- and late-onset Alzheimer disease.

OBJECTIVE: To characterize the clinical and molecular effect of mutations in the sortilin-related receptor (SORL1) gene. METHODS: We performed whole-exome sequencing in early-onset Alzheimer disease (EOAD) and late-onset Alzheimer disease (LOAD) families followed by functional studies of select variants. The phenotypic consequences associated with SORL1 mutations were characterized based on clinical reviews of medical records. Functional studies were completed to evaluate β-amyloid (Aβ) production and amyloid precursor protein (APP) trafficking associated with SORL1 mutations. RESULTS: SORL1 alterations were present in 2 EOAD families. In one, a SORL1 T588I change was identified in 4 individuals with AD, 2 of whom had parkinsonian features. In the second, an SORL1 T2134 alteration was found in 3 of 4 AD cases, one of whom had postmortem Lewy bodies. Among LOAD cases, 4 individuals with either SORL1 A528T or T947M alterations had parkinsonian features. Functionally, the variants weaken the interaction of the SORL1 protein with full-length APP, altering levels of Aβ and interfering with APP trafficking. CONCLUSIONS: The findings from this study support an important role for SORL1 mutations in AD pathogenesis by way of altering Aβ levels and interfering with APP trafficking. In addition, the presence of parkinsonian features among select individuals with AD and SORL1 mutations merits further investigation

Epigenetic and genetic risk of Alzheimer disease from autopsied brains in two ethnic groups.

Genetic variants and epigenetic features&nbsp;both contribute to the risk of Alzheimer's disease (AD). We studied the AD association of CpG-related single nucleotide polymorphisms (CGS), which act as a hub of both the genetic and epigenetic effects, in Caribbean Hispanics (CH) and generalized the findings to Non-Hispanic Whites (NHW). First, we conducted a genome-wide, sliding-window-based association with AD, in 7,155 CH and 1,283 NHW&nbsp;participants. Next, using&nbsp;data from the dorsolateral prefrontal cortex in 179 CH brains, we tested the cis- and trans-effects of AD-associated CGS on brain DNA methylation to mRNA expression. For the genes with significant cis- and trans-effects, we investigated their enriched pathways. We identified six genetic loci in CH with CGS dosage associated with AD at genome-wide significance levels: ADAM20 (Score = 55.19, P = 4.06 × 10-8), the intergenic region between VRTN and SYNDIG1L (Score = - 37.67, P = 2.25 × 10-9), SPG7 (16q24.3) (Score = 40.51, P = 2.23 × 10-8), PVRL2 (Score = 125.86, P = 1.64 × 10-9), TOMM40 (Score = - 18.58, P = 4.61 × 10-8), and APOE (Score = 75.12, P = 7.26 × 10-26). CGSes in PVRL2 and APOE were also significant in NHW. Except for ADAM20, CGSes in the other five loci were associated with CH brain methylation levels (mQTLs) and CGSes in SPG7, PVRL2, and APOE were also mQTLs in NHW. Except for SYNDIG1L (P = 0.08), brain methylation levels in the other five loci affected downstream mRNA expression in CH (P &lt; 0.05), and methylation at VRTN and TOMM40 were also associated with mRNA expression in NHW. Gene expression in these six loci were also regulated by CpG sites in genes that were enriched in the neuron projection and glutamatergic synapse pathways&nbsp;(FDR &lt; 0.05). DNA methylation at all six loci and mRNA expression of SYNDIG1 and TOMM40 were significantly associated with Braak Stage in CH. In summary,&nbsp;we identified six CpG-related genetic loci associated with AD in CH, harboring both genetic and epigenetic risks. However, their downstream effects on mRNA expression maybe ethnic specific and different from NHW

ABCA7 frameshift deletion associated with Alzheimer disease in African Americans

Objective: To identify a causative variant(s) that may contribute to Alzheimer disease (AD) in African Americans (AA) in the ATP-binding cassette, subfamily A (ABC1), member 7 (ABCA7) gene, a known risk factor for late-onset AD. Methods: Custom capture sequencing was performed on ∼150 kb encompassing ABCA7 in 40 AA cases and 37 AA controls carrying the AA risk allele (rs115550680). Association testing was performed for an ABCA7 deletion identified in large AA data sets (discovery n = 1,068; replication n = 1,749) and whole exome sequencing of Caribbean Hispanic (CH) AD families. Results: A 44-base pair deletion (rs142076058) was identified in all 77 risk genotype carriers, which shows that the deletion is in high linkage disequilibrium with the risk allele. The deletion was assessed in a large data set (531 cases and 527 controls) and, after adjustments for age, sex, and APOE status, was significantly associated with disease (p = 0.0002, odds ratio [OR] = 2.13 [95% confidence interval (CI): 1.42–3.20]). An independent data set replicated the association (447 cases and 880 controls, p = 0.0117, OR = 1.65 [95% CI: 1.12–2.44]), and joint analysis increased the significance (p = 1.414 × 10−5, OR = 1.81 [95% CI: 1.38–2.37]). The deletion is common in AA cases (15.2%) and AA controls (9.74%), but in only 0.12% of our non-Hispanic white cohort. Whole exome sequencing of multiplex, CH families identified the deletion cosegregating with disease in a large sibship. The deleted allele produces a stable, detectable RNA strand and is predicted to result in a frameshift mutation (p.Arg578Alafs) that could interfere with protein function. Conclusions: This common ABCA7 deletion could represent an ethnic-specific pathogenic alteration in AD

Episodic memory performance in a multi-ethnic longitudinal study of 13,037 elderly

Age-related changes in memory are not uniform, even in the absence of dementia. Characterization of non-disease associated cognitive changes is crucial to gain a more complete understanding of brain aging. Episodic memory was investigated in 13,037 ethnically diverse elderly (ages 72 to 85 years) with two to 15 years of follow-up, and with known dementia status, age, sex, education, and APOE genotypes. Adjusted trajectories of episodic memory performance over time were estimated using Latent Class Mixed Models. Analysis was conducted using two samples at baseline evaluation: i) non-cognitively impaired individuals, and ii) all individuals regardless of dementia status. We calculated the age-specific annual incidence rates of dementia in the non-demented elderly (n = 10,220). Two major episodic memory trajectories were estimated: 1) Stable—consisting of individuals exhibiting a constant or improved memory function, and 2) Decliner—consisting of individuals whose memory function declined. The majority of the study participants maintain their memory performance over time. Compared to those with Stable trajectory, individuals characterized as Decliners were more likely to have non-white ethnic background, fewer years of education, a higher frequency of ε4 allele at APOE gene and five times more likely to develop dementia. The steepest decline in episodic memory was observed in Caribbean-Hispanics compared to non-Hispanic whites (p = 4.3 x 10−15). The highest incident rates of dementia were observed in the oldest age group, among those of Caribbean-Hispanics ancestry and among Decliners who exhibited rates five times higher than those with Stable trajectories (11 per 100 person-years versus 3 per 100 person-years. Age, education, ethnic background and APOE genotype influence the maintenance of episodic memory. Declining memory is one of the strongest predictors of incident dementia

Sex differences in the genetic architecture of cognitive resilience to Alzheimer\u27s disease.

Approximately 30% of elderly adults are cognitively unimpaired at time of death despite the presence of Alzheimer\u27s disease neuropathology at autopsy. Studying individuals who are resilient to the cognitive consequences of Alzheimer\u27s disease neuropathology may uncover novel therapeutic targets to treat Alzheimer\u27s disease. It is well established that there are sex differences in response to Alzheimer\u27s disease pathology, and growing evidence suggests that genetic factors may contribute to these differences. Taken together, we sought to elucidate sex-specific genetic drivers of resilience. We extended our recent large scale genomic analysis of resilience in which we harmonized cognitive data across four cohorts of cognitive ageing, in vivo amyloid PET across two cohorts, and autopsy measures of amyloid neuritic plaque burden across two cohorts. These data were leveraged to build robust, continuous resilience phenotypes. With these phenotypes, we performed sex-stratified [n (males) = 2093, n (females) = 2931] and sex-interaction [n (both sexes) = 5024] genome-wide association studies (GWAS), gene and pathway-based tests, and genetic correlation analyses to clarify the variants, genes and molecular pathways that relate to resilience in a sex-specific manner. Estimated among cognitively normal individuals of both sexes, resilience was 20-25% heritable, and when estimated in either sex among cognitively normal individuals, resilience was 15-44% heritable. In our GWAS, we identified a female-specific locus on chromosome 10 [rs827389, β (females) = 0.08, P (females) = 5.76 × 10-09, β (males) = -0.01, P(males) = 0.70, β (interaction) = 0.09, P (interaction) = 1.01 × 10-04] in which the minor allele was associated with higher resilience scores among females. This locus is located within chromatin loops that interact with promoters of genes involved in RNA processing, including GATA3. Finally, our genetic correlation analyses revealed shared genetic architecture between resilience phenotypes and other complex traits, including a female-specific association with frontotemporal dementia and male-specific associations with heart rate variability traits. We also observed opposing associations between sexes for multiple sclerosis, such that more resilient females had a lower genetic susceptibility to multiple sclerosis, and more resilient males had a higher genetic susceptibility to multiple sclerosis. Overall, we identified sex differences in the genetic architecture of resilience, identified a female-specific resilience locus and highlighted numerous sex-specific molecular pathways that may underly resilience to Alzheimer\u27s disease pathology. This study illustrates the need to conduct sex-aware genomic analyses to identify novel targets that are unidentified in sex-agnostic models. Our findings support the theory that the most successful treatment for an individual with Alzheimer\u27s disease may be personalized based on their biological sex and genetic context