Genetic association studies of Alzheimer disease using multi-phenotype tests and gene-based tests

The genome-wide association study (GWAS) approach has identified novel loci for a variety of complex diseases. However, for most of these disorder much of the heritability is not explained by this approach, which focuses on identifying common variants that are associated with disease risk. The unexplained heritability may be due to genetic or phenotypic heterogeneity or the influence of rare variants. The motivation behind this thesis was to uncover the unexplained heritability by applying joint analyses of sets of variants (gene-based association test) and multiple disease-related phenotypes (called multivariate gene-based association test). First, we evaluated multivariate gene-based methods for detecting association of common genetic variants with correlated phenotypes. An extensive simulation study showed that the method combining the MultiPhen and GATES software performed best for most tested scenarios especially when correlations among phenotypes are relatively low. We developed a new multivariate gene-based test using rare variants called VEMPHAS. A simulation study using VEMPHAS showed that this method correctly controls for type I error in all tested scenarios. We applied VEMPHAS to analysis of various phenotypes related to Alzheimer disease (AD) and found suggestive association (P < 4.15x10-6) with the gene TRIM22, which has been identified in a previous sequencing study of AD onset in PSEN1/2 mutation carriers. We also developed software with a graphical user interface which is designed for integrating information from different types of data sources including genetic data (from GWAS or sequencing), expression data (from RNA-Seq), and protein structures (from protein data banks). This software has several features including 1) testing associations between genetic variants and gene expressions; 2) locating amino acids, encoded by the variants, in a protein structure; and 3) retrieving genetic locations (chromosome and base pair positions) of amino acids of interest in the protein structure. The last feature can be applied for prioritizing coding variants for gene-based association testing. The methods and strategies developed for this dissertation project can effectively uncover a portion of the remaining heritability of complex diseases that is unexplained by traditional GWAS approaches

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

Clinical Outcomes of Single-level Posterior Percutaneous Endoscopic Cervical Foraminotomy for Patients with Less Cervical Lordosis

Objective Posterior percutaneous endoscopic cervical foraminotomy (P-PECF) is a minimally invasive surgical technique for treatment of cervical radiculopathies. Application of P-PECF to patients with preexisting loss of cervical curvature (<10°) is still controversial because violation to facet joint may lead to kyphotic change. Clinical outcomes of P-PECF was analyzed and compared according to preoperative cervical curvature. Methods In this retrospective nested case-control study, 71 patients who underwent P-PECF due to foraminal soft disc herniation or bony stenosis were reviewed. P-PECF was performed by a single senior surgeon, and surgical methods were as previously described. Visual analogue pain scale on arm (Arm-VAS) was assessed preoperatively and postoperatively (1, 3, 6, 12 month and yearly thereafter). All patients were clinically followed for 24.5±20.0 months. The minimal clinically important difference of the Arm-VAS was set at 2.5. Patients with preoperative cervical curvature ≥10° were included in group I (n=32) and cervical curvature <10° or kyphosis were included in the group II (n=39). Results At the last follow up, 68/71 (96%) patients showed significant reduction of arm pain (Pre-operation, 7.4±2.0; post-operation, 1.5±2.0) after 1.74±0.29 months (95% CI; 1.18-2.31). The preoperative cervical curvature did not influence the outcome (p=0.4, T-test) and time to reach the clinical endpoint (p=0.34, Cox-logistic regression analysis). Conclusion P-PECF effectively reduced radicular pain due to foraminal soft disc herniation or stenosis. Preexisting loss of lordosis is not a risk factor for outcomes of P-PECF

PLXNA4 is associated with Alzheimer disease and modulates tau phosphorylation

OBJECTIVE: Much of the genetic basis for Alzheimer disease (AD) is unexplained. We sought to identify novel AD loci using a unique family-based approach that can detect robust associations with infrequent variants (minor allele frequency < 0.10). METHODS: We conducted a genome-wide association study in the Framingham Heart Study (discovery) and NIA-LOAD (National Institute on Aging-Late-Onset Alzheimer Disease) Study (replication) family-based cohorts using an approach that accounts for family structure and calculates a risk score for AD as the outcome. Links between the most promising gene candidate and AD pathogenesis were explored in silico as well as experimentally in cell-based models and in human brain. RESULTS: Genome-wide significant association was identified with a PLXNA4 single nucleotide polymorphism (rs277470) located in a region encoding the semaphorin-3A (SEMA3A) binding domain (meta-analysis p value [meta-P] = 4.1 × 10(-8) ). A test for association with the entire region was also significant (meta-P = 3.2 × 10(-4) ). Transfection of SH-SY5Y cells or primary rat neurons with full-length PLXNA4 (TS1) increased tau phosphorylation with stimulated by SEMA3A. The opposite effect was observed when cells were transfected with shorter isoforms (TS2 and TS3). However, transfection of any isoform into HEK293 cells stably expressing amyloid β (Aβ) precursor protein (APP) did not result in differential effects on APP processing or Aβ production. Late stage AD cases (n = 9) compared to controls (n = 5) had 1.9-fold increased expression of TS1 in cortical brain tissue (p = 1.6 × 10(-4) ). Expression of TS1 was significantly correlated with the Clinical Dementia Rating score (ρ = 0.75, p = 2.2 × 10(-4) ), plaque density (ρ = 0.56, p = 0.01), and Braak stage (ρ = 0.54, p = 0.02). INTERPRETATION: Our results indicate that PLXNA4 has a role in AD pathogenesis through isoform-specific effects on tau phosphorylation

Genome-wide pleiotropy analysis of neuropathological traits related to Alzheimer’s disease

Background Simultaneous consideration of two neuropathological traits related to Alzheimer’s disease (AD) has not been attempted in a genome-wide association study. Methods We conducted genome-wide pleiotropy analyses using association summary statistics from the Beecham et al. study (PLoS Genet 10:e1004606, 2014) for AD-related neuropathological traits, including neuritic plaque (NP), neurofibrillary tangle (NFT), and cerebral amyloid angiopathy (CAA). Significant findings were further examined by expression quantitative trait locus and differentially expressed gene analyses in AD vs. control brains using gene expression data. Results Genome-wide significant pleiotropic associations were observed for the joint model of NP and NFT (NP + NFT) with the single-nucleotide polymorphism (SNP) rs34487851 upstream of C2orf40 (alias ECRG4, P = 2.4 × 10−8) and for the joint model of NFT and CAA (NFT + CAA) with the HDAC9 SNP rs79524815 (P = 1.1 × 10−8). Gene-based testing revealed study-wide significant associations (P ≤ 2.0 × 10−6) for the NFT + CAA outcome with adjacent genes TRAPPC12, TRAPPC12-AS1, and ADI1. Risk alleles of proxy SNPs for rs79524815 were associated with significantly lower expression of HDAC9 in the brain (P = 3.0 × 10−3), and HDAC9 was significantly downregulated in subjects with AD compared with control subjects in the prefrontal (P = 7.9 × 10−3) and visual (P = 5.6 × 10−4) cortices. Conclusions Our findings suggest that pleiotropy analysis is a useful approach to identifying novel genetic associations with complex diseases and their endophenotypes. Functional studies are needed to determine whether ECRG4 or HDAC9 is plausible as a therapeutic target

Influence of Genetic Variation in PDE3A on Endothelial Function and Stroke

We aimed to characterize the genetics of endothelial function and how this influences risk for cardiovascular diseases such as ischemic stroke. We integrated genetic data from a study of ultrasound flow-mediated dilatation of brachial artery in adolescents from ALSPAC (Avon Longitudinal Study of Parents and Children; n=5214) with a study of ischemic stroke (MEGASTROKE: n=60 341 cases and 452 969 controls) to identify variants that confer risk of ischemic stroke through altered endothelial function. We identified a variant in PDE3A (Phosphodiesterase 3A), encoding phosphodiesterase 3A, which was associated with flow-mediated dilatation in adolescents (9-12 years of age; beta[SE], 0.38 [0.070]; P=3.8x10(-8)) and confers risk of ischemic stroke (odds ratio, 1.04 [95% CI, 1.02-1.06]; P=5.2x10(-6)). Bayesian colocalization analyses showed the same underlying variation is likely to lead to both associations (posterior probability, 97%). The same variant was associated with flow-mediated dilatation in a second study in young adults (age, 24-27 years; beta[SE], 0.47 [0.23]; P=0.047) but not in older adults (beta[SE], -0.012 [0.13]; P=0.89). We conclude that a genetic variant in PDE3A influences endothelial function in early life and leads to increased risk of ischemic stroke. Subtle, measurable changes to the vasculature that are influenced by genetics also influence risk of ischemic stroke

Genome-wide association study of cerebral small vessel disease reveals established and novel loci

Intracerebral haemorrhage and small vessel ischaemic stroke (SVS) are the most acute manifestations of cerebral small vessel disease, with no established preventive approaches beyond hypertension management. Combined genome-wide association study (GWAS) of these two correlated diseases may improve statistical power to detect novel genetic factors for cerebral small vessel disease, elucidating underlying disease mechanisms that may form the basis for future treatments. Because intracerebral haemorrhage location is an adequate surrogate for distinct histopathological variants of cerebral small vessel disease (lobar for cerebral amyloid angiopathy and non-lobar for arteriolosclerosis), we performed GWAS of intracerebral haemorrhage by location in 1813 subjects (755 lobar and 1005 non-lobar) and 1711 stroke-free control subjects. Intracerebral haemorrhage GWAS results by location were meta-analysed with GWAS results for SVS from MEGASTROKE, using 'Multi-Trait Analysis of GWAS' (MTAG) to integrate summary data across traits and generate combined effect estimates. After combining intracerebral haemorrhage and SVS datasets, our sample size included 241 024 participants (6255 intracerebral haemorrhage or SVS cases and 233 058 control subjects). Genome-wide significant associations were observed for non-lobar intracerebral haemorrhage enhanced by SVS with rs2758605 [MTAG P-value (P) = 2.6 x 10(-8)] at 1q22;rs72932727 (P = 1.7 x 10(-8)) at 2q33;and rs9515201 (P = 5.3 x 10(-10)) at 13q34. In the GTEx gene expression library, rs2758605 (1q22), rs72932727 (2q33) and rs9515201 (13q34) are significant cis-eQTLs for PMF1 (P = 1 x 10(-4) in tibial nerve), NBEAL1, FAM117B and CARF (P<2.1 x 10(-7) in arteries) and COL4A2 and COL4A1 (P<0.01 in brain putamen), respectively. Leveraging S-PrediXcan for gene-based association testing with the predicted expression models in tissues related with nerve, artery, and non-lobar brain, we found that experiment-wide significant (P<8.5 x 10(-7)) associations at three genes at 2q33 including NBEAL1, FAM117B and WDR12 and genome-wide significant associations at two genes including ICA1L at 2q33 and ZCCHC14 at 16q24. Brain cell-type specific expression profiling libraries reveal that SEMA4A, SLC25A44 and PMF1 at 1q22 and COL4A1 and COL4A2 at 13q34 were mainly expressed in endothelial cells, while the genes at 2q33 (FAM117B, CARF and NBEAL1) were expressed in various cell types including astrocytes, oligodendrocytes and neurons. Our cross-phenotype genetic study of intracerebral haemorrhage and SVS demonstrates novel genome-wide associations for non-lobar intracerebral haemorrhage at 2q33 and 13q34. Our replication of the 1q22 locus previous seen in traditional GWAS of intracerebral haemorrhage, as well as the rediscovery of 13q34, which had previously been reported in candidate gene studies with other cerebral small vessel disease-related traits strengthens the credibility of applying this novel genome-wide approach across intracerebral haemorrhage and SVS

Two novel loci, COBL and SLC10A2, for Alzheimer's disease in African Americans

INTRODUCTION: African Americans' (AAs) late-onset Alzheimer's disease (LOAD) genetic risk profile is incompletely understood. Including clinical covariates in genetic analyses using informed conditioning might improve study power. METHODS: We conducted a genome-wide association study (GWAS) in AAs employing informed conditioning in 1825 LOAD cases and 3784 cognitively normal controls. We derived a posterior liability conditioned on age, sex, diabetes status, current smoking status, educational attainment, and affection status, with parameters informed by external prevalence information. We assessed association between the posterior liability and a genome-wide set of single-nucleotide polymorphisms (SNPs), controlling for APOE and ABCA7, identified previously in a LOAD GWAS of AAs. RESULTS: Two SNPs at novel loci, rs112404845 (P = 3.8 × 10-8), upstream of COBL, and rs16961023 (P = 4.6 × 10-8), downstream of SLC10A2, obtained genome-wide significant evidence of association with the posterior liability. DISCUSSION: An informed conditioning approach can detect LOAD genetic associations in AAs not identified by traditional GWAS

Rare missense functional variants at COL4A1 and COL4A2 in sporadic intracerebral Hhmorrhage

Objective: To test the genetic contribution of rare missense variants in COL4A1 and COL4A2 in which common variants are genetically associated with sporadic intracerebral hemorrhage (ICH), we performed rare variant analysis in multiple sequencing data for the risk for sporadic ICH. Methods: We performed sequencing across 559Kbp at 13q34 including COL4A1 and COL4A2 among 2,133 individuals (1,055 ICH cases; 1,078 controls) in US-based and 1,492 individuals (192 ICH cases; 1,189 controls) from Scotland-based cohorts, followed by sequence annotation, functional impact prediction, genetic association testing, and in silico thermodynamic modeling. Results: We identified 107 rare nonsynonymous variants in sporadic ICH, of which two missense variants, rs138269346 (COL4A1I110T) and rs201716258 (COL4A2H203L), were predicted to be highly functional and occurred in multiple ICH cases but not in controls from the US-based cohort. The minor allele of rs201716258 was also present in Scottish ICH patients, and rs138269346 was observed in two ICH-free controls with a history of hypertension and myocardial infarction. Rs138269346 was nominally associated with non-lobar ICH risk (P=0.05), but not with lobar ICH (P=0.08), while associations between rs201716258 and ICH subtypes were non-significant (P&gt;0.12). Both variants were considered pathogenic based on minor allele frequency (&lt;0.00035 in EUR), predicted functional impact (deleterious or probably damaging), and in silico modeling studies (substantially altered physical length and thermal stability of collagen). Conclusions: We identified rare missense variants in COL4A1/A2 in association with sporadic ICH. Our annotation and simulation studies suggest that these variants are highly functional and may represent targets for translational follow-up