Frequency of Variants in Mendelian Alzheimer's Disease Genes within the Alzheimer's Disease Sequencing Project (ADSP)

Background

Prior studies using the ADSP data examined variants within presenilin-2 ( PSEN2 ), presenilin-1 ( PSEN1 ), and amyloid precursor protein ( APP ) genes. However, previously-reported clinically-relevant variants and other predicted damaging missense (DM) variants have not been characterized in a newer release of the Alzheimer's Disease Sequencing Project (ADSP).

Objective

To characterize previously-reported clinically-relevant variants and DM variants in PSEN2, PSEN1, APP within the participants from the ADSP.

Methods

We identified rare variants (MAF <1%) previously-reported in PSEN2 , PSEN1, and APP in the available ADSP sample of 14,641 individuals with whole genome sequencing and 16,849 individuals with whole exome sequencing available for research-use (N total = 31,490). We additionally curated variants in these three genes from ClinVar, OMIM, and Alzforum and report carriers of variants in clinical databases as well as predicted DM variants in these genes.

Results

We detected 31 previously-reported clinically-relevant variants with alternate alleles observed within the ADSP: 4 variants in PSEN2 , 25 in PSEN1 , and 2 in APP . The overall variant carrier rate for the 31 clinically-relevant variants in the ADSP was 0.3%. We observed that 79.5% of the variant carriers were cases compared to 3.9% were controls. In those with AD, the mean age of onset of AD among carriers of these clinically-relevant variants was 19.6 ± 1.4 years earlier compared with non-carriers (p-value=7.8×10 -57 ).

Conclusion

A small proportion of individuals in the ADSP are carriers of a previously-reported clinically-relevant variant allele for AD and these participants have significantly earlier age of AD onset compared to non-carriers.

Key variants via the Alzheimer's Disease Sequencing Project whole genome sequence data

INTRODUCTION

Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci.

METHODS

We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases = 2184, N controls = 2383) and targeted analyses in subpopulations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants.

RESULTS

Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses.

DISCUSSION

This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Bone mineral density and the risk of incident dementia: A meta-analysis

BACKGROUND

It is not known whether bone mineral density (BMD) measured at baseline or as the rate of decline prior to baseline (prior bone loss) is a stronger predictor of incident dementia or Alzheimer's disease (AD).

METHODS

We performed a meta-analysis of three longitudinal studies, the Framingham Heart Study (FHS), the Rotterdam Study (RS), and the Rush Memory and Aging Project (MAP), modeling the time to diagnosis of dementia as a function of BMD measures accounting for covariates. We included individuals with one or two BMD assessments, aged ≥60 years, and free of dementia at baseline with follow-up available. BMD was measured at the hip femoral neck using dual-energy X-ray absorptiometry (DXA), or at the heel calcaneus using quantitative ultrasound to calculate estimated BMD (eBMD). BMD at study baseline ("baseline BMD") and annualized percentage change in BMD prior to baseline ("prior bone loss") were included as continuous measures. The primary outcome was incident dementia diagnosis within 10 years of baseline, and incident AD was a secondary outcome. Baseline covariates included age, sex, body mass index, ApoE4 genotype, and education.

RESULTS

The combined sample size across all three studies was 4431 with 606 incident dementia diagnoses, 498 of which were AD. A meta-analysis of baseline BMD across three studies showed higher BMD to have a significant protective association with incident dementia with a hazard ratio of 0.47 (95% CI: 0.23-0.96; p = 0.038) per increase in g/cm2 , or 0.91 (95% CI: 0.84-0.995) per standard deviation increase. We observed a significant association between prior bone loss and incident dementia with a hazard ratio of 1.30 (95% CI: 1.12-1.51; p < 0.001) per percent increase in prior bone loss only in the FHS cohort.

CONCLUSIONS

Baseline BMD but not prior bone loss was associated with incident dementia in a meta-analysis across three studies.

Key variants via Alzheimer's Disease Sequencing Project whole genome sequence data

INTRODUCTION

Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci.

METHODS

We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases=2,184, N controls=2,383) and targeted analyses in sub-populations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants.

RESULTS

Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses.

DISCUSSION

This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Exploiting family history in aggregation unit-based genetic association tests

The development of sequencing technology calls for new powerful methods to detect disease associations and lower the cost of sequencing studies. Family history (FH) contains information on disease status of relatives, adding valuable information about the probands' health problems and risk of diseases. Incorporating data from FH is a cost-effective way to improve statistical evidence in genetic studies, and moreover, overcomes limitations in study designs with insufficient cases or missing genotype information for association analysis. We proposed family history aggregation unit-based test (FHAT) and optimal FHAT (FHAT-O) to exploit available FH for rare variant association analysis. Moreover, we extended liability threshold model of case-control status and FH (LT-FH) method in aggregated unit-based methods and compared that with FHAT and FHAT-O. The computational efficiency and flexibility of the FHAT and FHAT-O were demonstrated through both simulations and applications. We showed that FHAT, FHAT-O, and LT-FH methods offer reasonable control of the type I error unless case/control ratio is unbalanced, in which case they result in smaller inflation than that observed with conventional methods excluding FH. We also demonstrated that FHAT and FHAT-O are more powerful than LT-FH and conventional methods in many scenarios. By applying FHAT and FHAT-O to the analysis of all cause dementia and hypertension using the exome sequencing data from the UK Biobank, we showed that our methods can improve significance for known regions. Furthermore, we replicated the previous associations in all cause dementia and hypertension and detected novel regions through the exome-wide analysis.

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.

Meta-analysis of genome-wide association studies identifies ancestry-specific associations underlying circulating total tau levels

Circulating total-tau levels can be used as an endophenotype to identify genetic risk factors for tauopathies and related neurological disorders. Here, we confirmed and better characterized the association of the 17q21 MAPT locus with circulating total-tau in 14,721 European participants and identified three novel loci in 953 African American participants (4q31, 5p13, and 6q25) at P < 5 × 10-8. We additionally detected 14 novel loci at P < 5 × 10-7, specific to either Europeans or African Americans. Using whole-exome sequence data in 2,279 European participants, we identified ten genes associated with circulating total-tau when aggregating rare variants. Our genetic study sheds light on genes reported to be associated with neurological diseases including stroke, Alzheimer's, and Parkinson's (F5, MAP1B, and BCAS3), with Alzheimer's pathological hallmarks (ADAMTS12, IL15, and FHIT), or with an important function in the brain (PARD3, ELFN2, UBASH3B, SLIT3, and NSD3), and suggests that the genetic architecture of circulating total-tau may differ according to ancestry.

Family history aggregation unit-based tests to detect rare genetic variant associations with application to the Framingham Heart Study

A challenge in standard genetic studies is maintaining good power to detect associations, especially for low prevalent diseases and rare variants. The traditional methods are most powerful when evaluating the association between variants in balanced study designs. Without accounting for family correlation and unbalanced case-control ratio, these analyses could result in inflated type I error. One cost-effective solution to increase statistical power is exploitation of available family history (FH) that contains valuable information about disease heritability. Here, we develop methods to address the aforementioned type I error issues while providing optimal power to analyze aggregates of rare variants by incorporating additional information from FH. With enhanced power in these methods exploiting FH and accounting for relatedness and unbalanced designs, we successfully detect genes with suggestive associations with Alzheimer disease, dementia, and type 2 diabetes by using the exome chip data from the Framingham Heart Study.

Common variants in Alzheimer's disease and risk stratification by polygenic risk scores

Genetic discoveries of Alzheimer's disease are the drivers of our understanding, and together with polygenetic risk stratification can contribute towards planning of feasible and efficient preventive and curative clinical trials. We first perform a large genetic association study by merging all available case-control datasets and by-proxy study results (discovery n = 409,435 and validation size n = 58,190). Here, we add six variants associated with Alzheimer's disease risk (near APP, CHRNE, PRKD3/NDUFAF7, PLCG2 and two exonic variants in the SHARPIN gene). Assessment of the polygenic risk score and stratifying by APOE reveal a 4 to 5.5 years difference in median age at onset of Alzheimer's disease patients in APOE ɛ4 carriers. Because of this study, the underlying mechanisms of APP can be studied to refine the amyloid cascade and the polygenic risk score provides a tool to select individuals at high risk of Alzheimer's disease.

Large‐scale sequencing studies expand the known genetic architecture of Alzheimer's disease

Introduction

Genes implicated by genome‐wide association studies and family‐based studies of Alzheimer's disease (AD) are largely discordant. We hypothesized that genes identified by sequencing studies like the Alzheimer's Disease Sequencing Project (ADSP) may bridge this gap and highlight shared biological mechanisms.

Methods

We performed structured literature review of genes prioritized by ADSP studies, genes underlying familial dementias, and genes nominated by genome‐wide association studies. Gene set enrichment analyses of each list identified enriched pathways.

Results

The genes prioritized by the ADSP, familial dementia studies, and genome‐wide association studies minimally overlapped. Each gene set identified dozens of enriched pathways, several of which were shared (e.g., regulation of amyloid beta clearance).

Discussion

Alternative study designs provide unique insights into AD genetics. Shared pathways enriched by different genes highlight their relevance to AD pathogenesis, while the patterns of pathway enrichment unique to each gene set provide additional targets for functional studies.

Evaluation of population stratification adjustment using genome-wide or exonic variants

Population stratification may cause an inflated type-I error and spurious association when assessing the association between genetic variations with an outcome. Many genetic association studies are now using exonic variants, which captures only 1% of the genome, however, population stratification adjustments have not been evaluated in the context of exonic variants. We compare the performance of two established approaches: principal components analysis (PCA) and mixed-effects models and assess the utility of genome-wide (GW) and exonic variants, by simulation and using a data set from the Framingham Heart Study. Our results illustrate that although the PCs and genetic relationship matrices computed by GW and exonic markers are different, the type-I error rate of association tests for common variants with additive effect appear to be properly controlled in the presence of population stratification. In addition, by considering single nucleotide variants (SNVs) that have different levels of confounding by population stratification, we also compare the power across multiple association approaches to account for population stratification such as PC-based corrections and mixed-effects models. We find that while these two methods achieve a similar power for SNVs that have a low or medium level of confounding by population stratification, mixed-effects model can reach a higher power for SNVs highly confounded by population stratification.

Analysis of Whole-Exome Sequencing Data for Alzheimer Disease Stratified by APOE Genotype

Importance

Previous genome-wide association studies of common variants identified associations for Alzheimer disease (AD) loci evident only among individuals with particular APOE alleles.

Objective

To identify APOE genotype-dependent associations with infrequent and rare variants using whole-exome sequencing.

Design, Setting, and Participants

The discovery stage included 10 441 non-Hispanic white participants in the Alzheimer Disease Sequencing Project. Replication was sought in 2 independent, whole-exome sequencing data sets (1766 patients with AD, 2906 without AD [controls]) and a chip-based genotype imputation data set (8728 patients with AD, 9808 controls). Bioinformatics and functional analyses were conducted using clinical, cognitive, neuropathologic, whole-exome sequencing, and gene expression data obtained from a longitudinal cohort sample including 402 patients with AD and 647 controls. Data were analyzed between March 2017 and September 2018.

Main Outcomes and Measures

Score, Firth, and sequence kernel association tests were used to test the association of AD risk with individual variants and genes in subgroups of APOE ε4 carriers and noncarriers. Results with P ≤ 1 × 10-5 were further evaluated in the replication data sets and combined by meta-analysis.

Results

Among 3145 patients with AD and 4213 controls lacking ε4 (mean [SD] age, 83.4 [7.6] years; 4363 [59.3.%] women), novel genome-wide significant associations were obtained in the discovery sample with rs536940594 in AC099552 (odds ratio [OR], 88.0; 95% CI, 9.08-852.0; P = 2.22 × 10-7) and rs138412600 in GPAA1 (OR, 1.78; 95% CI, 1.44-2.2; meta-P = 7.81 × 10-8). GPAA1 was also associated with expression in the brain of GPAA1 (β = -0.08; P = .03) and its repressive transcription factor, FOXG1 (β = 0.13; P = .003), and global cognition function (β = -0.53; P = .009). Significant gene-wide associations (threshold P ≤ 6.35 × 10-7) were observed for OR8G5 (P = 4.67 × 10-7), IGHV3-7 (P = 9.75 × 10-16), and SLC24A3 (P = 2.67 × 10-12) in 2377 patients with AD and 706 controls with ε4 (mean [SD] age, 75.2 [9.6] years; 1668 [54.1%] women).

Conclusions and Relevance

The study identified multiple possible novel associations for AD with individual and aggregated rare variants in groups of individuals with and without APOE ε4 alleles that reinforce known and suggest additional pathways leading to AD.

Exome Chip Analysis Identifies Low-Frequency and Rare Variants in MRPL38 for White Matter Hyperintensities on Brain Magnetic Resonance Imaging

Background and Purpose- White matter hyperintensities (WMH) on brain magnetic resonance imaging are typical signs of cerebral small vessel disease and may indicate various preclinical, age-related neurological disorders, such as stroke. Though WMH are highly heritable, known common variants explain a small proportion of the WMH variance. The contribution of low-frequency/rare coding variants to WMH burden has not been explored. Methods- In the discovery sample we recruited 20 719 stroke/dementia-free adults from 13 population-based cohort studies within the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, among which 17 790 were of European ancestry and 2929 of African ancestry. We genotyped these participants at ≈250 000 mostly exonic variants with Illumina HumanExome BeadChip arrays. We performed ethnicity-specific linear regression on rank-normalized WMH in each study separately, which were then combined in meta-analyses to test for association with single variants and genes aggregating the effects of putatively functional low-frequency/rare variants. We then sought replication of the top findings in 1192 adults (European ancestry) with whole exome/genome sequencing data from 2 independent studies. Results- At 17q25, we confirmed the association of multiple common variants in TRIM65, FBF1, and ACOX1 ( P<6×10^-7). We also identified a novel association with 2 low-frequency nonsynonymous variants in MRPL38 (lead, rs34136221; P_EA=4.5×10^-8) partially independent of known common signal ( P_{EA(conditional)}=1.4×10^-3). We further identified a locus at 2q33 containing common variants in NBEAL1, CARF, and WDR12 (lead, rs2351524; P_all=1.9×10^-10). Although our novel findings were not replicated because of limited power and possible differences in study design, meta-analysis of the discovery and replication samples yielded stronger association for the 2 low-frequency MRPL38 variants ( P_rs34136221=2.8×10^-8). Conclusions- Both common and low-frequency/rare functional variants influence WMH. Larger replication and experimental follow-up are essential to confirm our findings and uncover the biological causal mechanisms of age-related WMH.

Quality control and integration of genotypes from two calling pipelines for whole genome sequence data in the Alzheimer's disease sequencing project

The Alzheimer's Disease Sequencing Project (ADSP) performed whole genome sequencing (WGS) of 584 subjects from 111 multiplex families at three sequencing centers. Genotype calling of single nucleotide variants (SNVs) and insertion-deletion variants (indels) was performed centrally using GATK-HaplotypeCaller and Atlas V2. The ADSP Quality Control (QC) Working Group applied QC protocols to project-level variant call format files (VCFs) from each pipeline, and developed and implemented a novel protocol, termed "consensus calling," to combine genotype calls from both pipelines into a single high-quality set. QC was applied to autosomal bi-allelic SNVs and indels, and included pipeline-recommended QC filters, variant-level QC, and sample-level QC. Low-quality variants or genotypes were excluded, and sample outliers were noted. Quality was assessed by examining Mendelian inconsistencies (MIs) among 67 parent-offspring pairs, and MIs were used to establish additional genotype-specific filters for GATK calls. After QC, 578 subjects remained. Pipeline-specific QC excluded ~12.0% of GATK and 14.5% of Atlas SNVs. Between pipelines, ~91% of SNV genotypes across all QCed variants were concordant; 4.23% and 4.56% of genotypes were exclusive to Atlas or GATK, respectively; the remaining ~0.01% of discordant genotypes were excluded. For indels, variant-level QC excluded ~36.8% of GATK and 35.3% of Atlas indels. Between pipelines, ~55.6% of indel genotypes were concordant; while 10.3% and 28.3% were exclusive to Atlas or GATK, respectively; and ~0.29% of discordant genotypes were. The final WGS consensus dataset contains 27,896,774 SNVs and 3,133,926 indels and is publicly available.

Publisher Correction: Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genetic meta-analysis of diagnosed Alzheimer's disease identifies new risk loci and implicates Aβ, tau, immunity and lipid processing

Risk for late-onset Alzheimer's disease (LOAD), the most prevalent dementia, is partially driven by genetics. To identify LOAD risk loci, we performed a large genome-wide association meta-analysis of clinically diagnosed LOAD (94,437 individuals). We confirm 20 previous LOAD risk loci and identify five new genome-wide loci (IQCK, ACE, ADAM10, ADAMTS1, and WWOX), two of which (ADAM10, ACE) were identified in a recent genome-wide association (GWAS)-by-familial-proxy of Alzheimer's or dementia. Fine-mapping of the human leukocyte antigen (HLA) region confirms the neurological and immune-mediated disease haplotype HLA-DR15 as a risk factor for LOAD. Pathway analysis implicates immunity, lipid metabolism, tau binding proteins, and amyloid precursor protein (APP) metabolism, showing that genetic variants affecting APP and Aβ processing are associated not only with early-onset autosomal dominant Alzheimer's disease but also with LOAD. Analyses of risk genes and pathways show enrichment for rare variants (P = 1.32 × 10-7), indicating that additional rare variants remain to be identified. We also identify important genetic correlations between LOAD and traits such as family history of dementia and education.

Genetic and lifestyle risk factors for MRI-defined brain infarcts in a population-based setting

OBJECTIVE

To explore genetic and lifestyle risk factors of MRI-defined brain infarcts (BI) in large population-based cohorts.

METHODS

We performed meta-analyses of genome-wide association studies (GWAS) and examined associations of vascular risk factors and their genetic risk scores (GRS) with MRI-defined BI and a subset of BI, namely, small subcortical BI (SSBI), in 18 population-based cohorts (n = 20,949) from 5 ethnicities (3,726 with BI, 2,021 with SSBI). Top loci were followed up in 7 population-based cohorts (n = 6,862; 1,483 with BI, 630 with SBBI), and we tested associations with related phenotypes including ischemic stroke and pathologically defined BI.

RESULTS

The mean prevalence was 17.7% for BI and 10.5% for SSBI, steeply rising after age 65. Two loci showed genome-wide significant association with BI: FBN2, p = 1.77 × 10-8; and LINC00539/ZDHHC20, p = 5.82 × 10-9. Both have been associated with blood pressure (BP)-related phenotypes, but did not replicate in the smaller follow-up sample or show associations with related phenotypes. Age- and sex-adjusted associations with BI and SSBI were observed for BP traits (p value for BI, p [BI] = 9.38 × 10-25; p [SSBI] = 5.23 × 10-14 for hypertension), smoking (p [BI] = 4.4 × 10-10; p [SSBI] = 1.2 × 10-4), diabetes (p [BI] = 1.7 × 10-8; p [SSBI] = 2.8 × 10-3), previous cardiovascular disease (p [BI] = 1.0 × 10-18; p [SSBI] = 2.3 × 10-7), stroke (p [BI] = 3.9 × 10-69; p [SSBI] = 3.2 × 10-24), and MRI-defined white matter hyperintensity burden (p [BI] = 1.43 × 10-157; p [SSBI] = 3.16 × 10-106), but not with body mass index or cholesterol. GRS of BP traits were associated with BI and SSBI (p ≤ 0.0022), without indication of directional pleiotropy.

CONCLUSION

In this multiethnic GWAS meta-analysis, including over 20,000 population-based participants, we identified genetic risk loci for BI requiring validation once additional large datasets become available. High BP, including genetically determined, was the most significant modifiable, causal risk factor for BI.

Network analysis of drug effect on triglyceride-associated DNA methylation

BACKGROUND

DNA methylation, an epigenetic modification, can be affected by environmental factors and thus regulate gene expression levels that can lead to alterations of certain phenotypes. Network analysis has been used successfully to discover gene sets that are expressed differently across multiple disease states and suggest possible pathways of disease progression. We applied this framework to compare DNA methylation levels before and after lipid-lowering medication and to identify modules that differ topologically between the two time points, revealing the association between lipid medication and these triglyceride-related methylation sites.

METHODS

We performed quality control using beta-mixture quantile normalization on 463,995 cytosine-phosphate-guanine (CpG) sites and deleted problematic sites, resulting in 423,004 probes. We identified 14,850 probes that were nominally associated with triglycerides prior to treatment and performed weighted gene correlation network analysis (WGCNA) to construct pre- and posttreatment methylation networks of these probes. We then applied both WGCNA module preservation and generalized Hamming distance (GHD) to identify modules with topological differences between the pre- and posttreatment. For modules with structural changes between 2 time points, we performed pathway-enrichment analysis to gain further insight into the biological function of the genes from these modules.

RESULTS

Six triglyceride-associated modules were identified using pretreatment methylation probes. The same 3 modules were not preserved in posttreatment data using both the module-preservation and the GHD methods. Top-enriched pathways for the 3 differentially methylated modules are sphingolipid signaling pathway, proteoglycans in cancer, and metabolic pathways (p values < 0.005). One module in particular included an enrichment of lipid-related pathways among the top results.

CONCLUSIONS

The same 3 modules, which were differentially methylated between pre- and posttreatment, were identified using both WGCNA module-preservation and GHD methods. Pathway analysis revealed that triglyceride-associated modules contain groups of genes that are involved in lipid signaling and metabolism. These 3 modules may provide insight into the effect of fenofibrate on changes in triglyceride levels and these methylation sites.

Integrative methylation score to identify epigenetic modifications associated with lipid changes resulting from fenofibrate treatment in families

Epigenome-wide association studies (EWAS) have traditionally focused on the association test of single epigenetic markers with complex traits. However, it is possible that multiple cytosine-phosphate-guanine (CpG) sites at the same locus could jointly exert their effects on human traits. Therefore, a region-based test that combines multiple markers could be more powerful. We used 2 different region-based tests to investigate the association between changes in DNA methylation and drug response, including the median methylation level test (MMLT) and sequence kernel association test (SKAT). No genes were found to be significantly associated with the drug response (for triglycerides, the false discovery rate ranged from 0.855 to 0.999; for high-density lipoprotein cholesterol, and the false discovery rate ranged from 0.584 to 0.915). Further evidence is needed to explore potential application of gene-level methylation association analysis.

Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes

Stroke has multiple etiologies, but the underlying genes and pathways are largely unknown. We conducted a multiancestry genome-wide-association meta-analysis in 521,612 individuals (67,162 cases and 454,450 controls) and discovered 22 new stroke risk loci, bringing the total to 32. We further found shared genetic variation with related vascular traits, including blood pressure, cardiac traits, and venous thromboembolism, at individual loci (n = 18), and using genetic risk scores and linkage-disequilibrium-score regression. Several loci exhibited distinct association and pleiotropy patterns for etiological stroke subtypes. Eleven new susceptibility loci indicate mechanisms not previously implicated in stroke pathophysiology, with prioritization of risk variants and genes accomplished through bioinformatics analyses using extensive functional datasets. Stroke risk loci were significantly enriched in drug targets for antithrombotic therapy.

Evaluation of a Genetic Risk Score to Improve Risk Prediction for Alzheimer's Disease

Effective prevention of Alzheimer's disease (AD) requires the development of risk prediction tools permitting preclinical intervention. We constructed a genetic risk score (GRS) comprising common genetic variants associated with AD, evaluated its association with incident AD and assessed its capacity to improve risk prediction over traditional models based on age, sex, education, and APOEɛ4. In eight prospective cohorts included in the International Genomics of Alzheimer's Project (IGAP), we derived weighted sum of risk alleles from the 19 top SNPs reported by the IGAP GWAS in participants aged 65 and older without prevalent dementia. Hazard ratios (HR) of incident AD were estimated in Cox models. Improvement in risk prediction was measured by the difference in C-index (Δ-C), the integrated discrimination improvement (IDI) and continuous net reclassification improvement (NRI>0). Overall, 19,687 participants at risk were included, of whom 2,782 developed AD. The GRS was associated with a 17% increase in AD risk (pooled HR = 1.17; 95% CI =   [1.13-1.21] per standard deviation increase in GRS; p-value =  2.86×10-16). This association was stronger among persons with at least one APOEɛ4 allele (HRGRS = 1.24; 95% CI =   [1.15-1.34]) than in others (HRGRS = 1.13; 95% CI =   [1.08-1.18]; pinteraction = 3.45×10-2). Risk prediction after seven years of follow-up showed a small improvement when adding the GRS to age, sex, APOEɛ4, and education (Δ-Cindex =  0.0043 [0.0019-0.0067]). Similar patterns were observed for IDI and NRI>0. In conclusion, a risk score incorporating common genetic variation outside the APOEɛ4 locus improved AD risk prediction and may facilitate risk stratification for prevention trials.

Whole genome sequence analyses of brain imaging measures in the Framingham Study

Objective

We sought to identify rare variants influencing brain imaging phenotypes in the Framingham Heart Study by performing whole genome sequence association analyses within the Trans-Omics for Precision Medicine Program.

Methods

We performed association analyses of cerebral and hippocampal volumes and white matter hyperintensity (WMH) in up to 2,180 individuals by testing the association of rank-normalized residuals from mixed-effect linear regression models adjusted for sex, age, and total intracranial volume with individual variants while accounting for familial relatedness. We conducted gene-based tests for rare variants using (1) a sliding-window approach, (2) a selection of functional exonic variants, or (3) all variants.

Results

We detected new loci in 1p21 for cerebral volume (minor allele frequency [MAF] 0.005, p = 10⁻⁸) and in 16q23 for hippocampal volume (MAF 0.05, p = 2.7 × 10⁻⁸). Previously identified associations in 12q24 for hippocampal volume (rs7294919, p = 4.4 × 10⁻⁴) and in 17q25 for WMH (rs7214628, p = 2.0 × 10⁻³) were confirmed. Gene-based tests detected associations (p ≤ 2.3 × 10⁻⁶) in new loci for cerebral (5q13, 8p12, 9q31, 13q12-q13, 15q24, 17q12, 19q13) and hippocampal volumes (2p12) and WMH (3q13, 4p15) including Alzheimer disease– (UNC5D) and Parkinson disease–associated genes (GBA). Pathway analyses evidenced enrichment of associated genes in immunity, inflammation, and Alzheimer disease and Parkinson disease pathways.

Conclusions

Whole genome sequence–wide search reveals intriguing new loci associated with brain measures. Replication of novel loci is needed to confirm these findings.

Rare coding variants in PLCG2, ABI3, and TREM2 implicate microglial-mediated innate immunity in Alzheimer's disease

We identified rare coding variants associated with Alzheimer's disease in a three-stage case-control study of 85,133 subjects. In stage 1, we genotyped 34,174 samples using a whole-exome microarray. In stage 2, we tested associated variants (P < 1 × 10-4) in 35,962 independent samples using de novo genotyping and imputed genotypes. In stage 3, we used an additional 14,997 samples to test the most significant stage 2 associations (P < 5 × 10-8) using imputed genotypes. We observed three new genome-wide significant nonsynonymous variants associated with Alzheimer's disease: a protective variant in PLCG2 (rs72824905: p.Pro522Arg, P = 5.38 × 10-10, odds ratio (OR) = 0.68, minor allele frequency (MAF)cases = 0.0059, MAFcontrols = 0.0093), a risk variant in ABI3 (rs616338: p.Ser209Phe, P = 4.56 × 10-10, OR = 1.43, MAFcases = 0.011, MAFcontrols = 0.008), and a new genome-wide significant variant in TREM2 (rs143332484: p.Arg62His, P = 1.55 × 10-14, OR = 1.67, MAFcases = 0.0143, MAFcontrols = 0.0089), a known susceptibility gene for Alzheimer's disease. These protein-altering changes are in genes highly expressed in microglia and highlight an immune-related protein-protein interaction network enriched for previously identified risk genes in Alzheimer's disease. These genetic findings provide additional evidence that the microglia-mediated innate immune response contributes directly to the development of Alzheimer's disease.

Evaluation of logistic regression models and effect of covariates for case-control study in RNA-Seq analysis

Background

Next generation sequencing provides a count of RNA molecules in the form of short reads, yielding discrete, often highly non-normally distributed gene expression measurements. Although Negative Binomial (NB) regression has been generally accepted in the analysis of RNA sequencing (RNA-Seq) data, its appropriateness has not been exhaustively evaluated. We explore logistic regression as an alternative method for RNA-Seq studies designed to compare cases and controls, where disease status is modeled as a function of RNA-Seq reads using simulated and Huntington disease data. We evaluate the effect of adjusting for covariates that have an unknown relationship with gene expression. Finally, we incorporate the data adaptive method in order to compare false positive rates.

Results

When the sample size is small or the expression levels of a gene are highly dispersed, the NB regression shows inflated Type-I error rates but the Classical logistic and Bayes logistic (BL) regressions are conservative. Firth’s logistic (FL) regression performs well or is slightly conservative. Large sample size and low dispersion generally make Type-I error rates of all methods close to nominal alpha levels of 0.05 and 0.01. However, Type-I error rates are controlled after applying the data adaptive method. The NB, BL, and FL regressions gain increased power with large sample size, large log2 fold-change, and low dispersion. The FL regression has comparable power to NB regression.

Conclusions

We conclude that implementing the data adaptive method appropriately controls Type-I error rates in RNA-Seq analysis. Firth’s logistic regression provides a concise statistical inference process and reduces spurious associations from inaccurately estimated dispersion parameters in the negative binomial framework.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1498-y) contains supplementary material, which is available to authorized users.

Novel genetic loci underlying human intracranial volume identified through genome-wide association

Intracranial volume reflects the maximally attained brain size during development, and remains stable with loss of tissue in late life. It is highly heritable, but the underlying genes remain largely undetermined. In a genome-wide association study of 32,438 adults, we discovered five previously unknown loci for intracranial volume and confirmed two known signals. Four of the loci were also associated with adult human stature, but these remained associated with intracranial volume after adjusting for height. We found a high genetic correlation with child head circumference (ρgenetic = 0.748), which indicates a similar genetic background and allowed us to identify four additional loci through meta-analysis (Ncombined = 37,345). Variants for intracranial volume were also related to childhood and adult cognitive function, and Parkinson's disease, and were enriched near genes involved in growth pathways, including PI3K-AKT signaling. These findings identify the biological underpinnings of intracranial volume and their link to physiological and pathological traits.

Novel microRNA discovery using small RNA sequencing in post-mortem human brain

BACKGROUND

MicroRNAs (miRNAs) are short, non-coding RNAs that regulate gene expression mainly through translational repression of target mRNA molecules. More than 2700 human miRNAs have been identified and some are known to be associated with disease phenotypes and to display tissue-specific patterns of expression.

METHODS

We used high-throughput small RNA sequencing to discover novel miRNAs in 93 human post-mortem prefrontal cortex samples from individuals with Huntington's disease (n = 28) or Parkinson's disease (n = 29) and controls without neurological impairment (n = 36). A custom miRNA identification analysis pipeline was built, which utilizes miRDeep* miRNA identification and result filtering based on false positive rate estimates.

RESULTS

Ninety-nine novel miRNA candidates with a false positive rate of less than 5 % were identified. Thirty-four of the candidate miRNAs show sequence similarity with known mature miRNA sequences and may be novel members of known miRNA families, while the remaining 65 may constitute previously undiscovered families of miRNAs. Nineteen of the 99 candidate miRNAs were replicated using independent, publicly-available human brain RNA-sequencing samples, and seven were experimentally validated using qPCR.

CONCLUSIONS

We have used small RNA sequencing to identify 99 putative novel miRNAs that are present in human brain samples.

Rare Functional Variant in TM2D3 is Associated with Late-Onset Alzheimer's Disease

We performed an exome-wide association analysis in 1393 late-onset Alzheimer's disease (LOAD) cases and 8141 controls from the CHARGE consortium. We found that a rare variant (P155L) in TM2D3 was enriched in Icelanders (~0.5% versus <0.05% in other European populations). In 433 LOAD cases and 3903 controls from the Icelandic AGES sub-study, P155L was associated with increased risk and earlier onset of LOAD [odds ratio (95% CI) = 7.5 (3.5-15.9), p = 6.6x10-9]. Mutation in the Drosophila TM2D3 homolog, almondex, causes a phenotype similar to loss of Notch/Presenilin signaling. Human TM2D3 is capable of rescuing these phenotypes, but this activity is abolished by P155L, establishing it as a functionally damaging allele. Our results establish a rare TM2D3 variant in association with LOAD susceptibility, and together with prior work suggests possible links to the β-amyloid cascade.

Six Novel Loci Associated with Circulating VEGF Levels Identified by a Meta-analysis of Genome-Wide Association Studies

Vascular endothelial growth factor (VEGF) is an angiogenic and neurotrophic factor, secreted by endothelial cells, known to impact various physiological and disease processes from cancer to cardiovascular disease and to be pharmacologically modifiable. We sought to identify novel loci associated with circulating VEGF levels through a genome-wide association meta-analysis combining data from European-ancestry individuals and using a dense variant map from 1000 genomes imputation panel. Six discovery cohorts including 13,312 samples were analyzed, followed by in-silico and de-novo replication studies including an additional 2,800 individuals. A total of 10 genome-wide significant variants were identified at 7 loci. Four were novel loci (5q14.3, 10q21.3, 16q24.2 and 18q22.3) and the leading variants at these loci were rs114694170 (MEF2C, P = 6.79x10^-13), rs74506613 (JMJD1C, P = 1.17x10^-19), rs4782371 (ZFPM1, P = 1.59x10^-9) and rs2639990 (ZADH2, P = 1.72x10^-8), respectively. We also identified two new independent variants (rs34528081, VEGFA, P = 1.52x10^-18; rs7043199, VLDLR-AS1, P = 5.12x10^-14) at the 3 previously identified loci and strengthened the evidence for the four previously identified SNPs (rs6921438, LOC100132354, P = 7.39x10^-1467; rs1740073, C6orf223, P = 2.34x10^-17; rs6993770, ZFPM2, P = 2.44x10^-60; rs2375981, KCNV2, P = 1.48x10^-100). These variants collectively explained up to 52% of the VEGF phenotypic variance. We explored biological links between genes in the associated loci using Ingenuity Pathway Analysis that emphasized their roles in embryonic development and function. Gene set enrichment analysis identified the ERK5 pathway as enriched in genes containing VEGF associated variants. eQTL analysis showed, in three of the identified regions, variants acting as both cis and trans eQTLs for multiple genes. Most of these genes, as well as some of those in the associated loci, were involved in platelet biogenesis and functionality, suggesting the importance of this process in regulation of VEGF levels. This work also provided new insights into the involvement of genes implicated in various angiogenesis related pathologies in determining circulating VEGF levels. The understanding of the molecular mechanisms by which the identified genes affect circulating VEGF levels could be important in the development of novel VEGF-related therapies for such diseases.

Vascular Endothelial Growth Factor (VEGF) is a protein with a fundamental role in development of vascular system. The protein, produced by many types of cells, is released in the blood. High levels of VEGF have been observed in different pathological conditions especially in cancer, cardiovascular, and inflammatory diseases. Therefore, identifying the genetic factors influencing VEGF levels is important for predicting and treating such pathologies. The number of genetic variants associated with VEGF levels has been limited. To identify new loci, we have performed a Genome Wide Association Study meta-analysis on a sample of more than 16,000 individuals from 10 cohorts, using a high-density genetic map. This analysis revealed 10 variants associated with VEGF circulating levels, 6 of these being novel associations. The 10 variants cumulatively explain more than 50% of the variability of VEGF serum levels. Our analyses have identified genes known to be involved in angiogenesis related diseases and genes implicated in platelet metabolism, suggesting the importance of links between this process and VEGF regulation. Overall, these data have improved our understanding of the genetic variation underlying circulating VEGF levels. This in turn could guide our response to the challenge posed by various VEGF-related pathologies.

GWAS for executive function and processing speed suggests involvement of the CADM2 gene

To identify common variants contributing to normal variation in two specific domains of cognitive functioning, we conducted a genome-wide association study (GWAS) of executive functioning and information processing speed in non-demented older adults from the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) consortium. Neuropsychological testing was available for 5429-32,070 subjects of European ancestry aged 45 years or older, free of dementia and clinical stroke at the time of cognitive testing from 20 cohorts in the discovery phase. We analyzed performance on the Trail Making Test parts A and B, the Letter Digit Substitution Test (LDST), the Digit Symbol Substitution Task (DSST), semantic and phonemic fluency tests, and the Stroop Color and Word Test. Replication was sought in 1311-21860 subjects from 20 independent cohorts. A significant association was observed in the discovery cohorts for the single-nucleotide polymorphism (SNP) rs17518584 (discovery P-value=3.12 × 10(-8)) and in the joint discovery and replication meta-analysis (P-value=3.28 × 10(-9) after adjustment for age, gender and education) in an intron of the gene cell adhesion molecule 2 (CADM2) for performance on the LDST/DSST. Rs17518584 is located about 170 kb upstream of the transcription start site of the major transcript for the CADM2 gene, but is within an intron of a variant transcript that includes an alternative first exon. The variant is associated with expression of CADM2 in the cingulate cortex (P-value=4 × 10(-4)). The protein encoded by CADM2 is involved in glutamate signaling (P-value=7.22 × 10(-15)), gamma-aminobutyric acid (GABA) transport (P-value=1.36 × 10(-11)) and neuron cell-cell adhesion (P-value=1.48 × 10(-13)). Our findings suggest that genetic variation in the CADM2 gene is associated with individual differences in information processing speed.

Genome-wide linkage analyses of non-Hispanic white families identify novel loci for familial late-onset Alzheimer's disease

INTRODUCTION

Few high penetrance variants that explain risk in late-onset Alzheimer's disease (LOAD) families have been found.

METHODS

We performed genome-wide linkage and identity-by-descent (IBD) analyses on 41 non-Hispanic white families exhibiting likely dominant inheritance of LOAD, and having no mutations at known familial Alzheimer's disease (AD) loci, and a low burden of APOE ε4 alleles.

RESULTS

Two-point parametric linkage analysis identified 14 significantly linked regions, including three novel linkage regions for LOAD (5q32, 11q12.2-11q14.1, and 14q13.3), one of which replicates a genome-wide association LOAD locus, the MS4A6A-MS4A4E gene cluster at 11q12.2. Five of the 14 regions (3q25.31, 4q34.1, 8q22.3, 11q12.2-14.1, and 19q13.41) are supported by strong multipoint results (logarithm of odds [LOD*] ≥1.5). Nonparametric multipoint analyses produced an additional significant locus at 14q32.2 (LOD* = 4.18). The 1-LOD confidence interval for this region contains one gene, C14orf177, and the microRNA Mir_320, whereas IBD analyses implicates an additional gene BCL11B, a regulator of brain-derived neurotrophic signaling, a pathway associated with pathogenesis of several neurodegenerative diseases.

DISCUSSION

Examination of these regions after whole-genome sequencing may identify highly penetrant variants for familial LOAD.

Evaluation of power of the Illumina HumanOmni5M-4v1 BeadChip to detect risk variants for human complex diseases

Although emerging sequencing technologies can characterize all genetic variants, the cost is still high. Illumina released the HumanOmni5M-4v1 (Omni5) genotype array with ~4.3M assayed SNPs, a much denser array compared with other available arrays. The Omni5 balances both cost and array density. In this article, we illustrate the power of Omni5 to detect genetic associations. The Omni5 includes variants with a wide range of minor allele frequencies down to <1%. The theoretical power calculation examples indicate the increased power of the Omni5 array compared with other arrays with lower density when evaluating associations with some known loci, although there are exceptions. We further evaluate the genetic associations between known loci and several quantitative traits in the Framingham Heart Study: femoral neck bone mineral density, lumbar spine bone mineral density and hippocampal volume. Finally, we search genome wide for novel associations using the Omni5 genotypes. We compare our association results from Affymetrix 500K+MIPS 50K arrays and two imputed data sets: (1) HapMap Phase II and (2) 1000 Genomes reference panel. We observed increased evidence for genotype-phenotype associations with smaller P-values for selected known loci using the Omni5 genotypes. With limited sample sizes, we identify novel variants with genome-wide significant P-values. Our observations support the notion that dense genotyping using the Omni5 can be powerful in detecting novel associated variants. Comparison with imputed data with higher density also suggests that imputation helps but cannot replace genotyping, especially when imputation quality is low.

Rare and Coding Region Genetic Variants Associated With Risk of Ischemic Stroke: The NHLBI Exome Sequence Project

IMPORTANCE

Stroke is the second leading cause of death and the third leading cause of years of life lost. Genetic factors contribute to stroke prevalence, and candidate gene and genome-wide association studies (GWAS) have identified variants associated with ischemic stroke risk. These variants often have small effects without obvious biological significance. Exome sequencing may discover predicted protein-altering variants with a potentially large effect on ischemic stroke risk.

OBJECTIVE

To investigate the contribution of rare and common genetic variants to ischemic stroke risk by targeting the protein-coding regions of the human genome.

DESIGN, SETTING, AND PARTICIPANTS

The National Heart, Lung, and Blood Institute (NHLBI) Exome Sequencing Project (ESP) analyzed approximately 6000 participants from numerous cohorts of European and African ancestry. For discovery, 365 cases of ischemic stroke (small-vessel and large-vessel subtypes) and 809 European ancestry controls were sequenced; for replication, 47 affected sibpairs concordant for stroke subtype and an African American case-control series were sequenced, with 1672 cases and 4509 European ancestry controls genotyped. The ESP's exome sequencing and genotyping started on January 1, 2010, and continued through June 30, 2012. Analyses were conducted on the full data set between July 12, 2012, and July 13, 2013.

MAIN OUTCOMES AND MEASURES

Discovery of new variants or genes contributing to ischemic stroke risk and subtype (primary analysis) and determination of support for protein-coding variants contributing to risk in previously published candidate genes (secondary analysis).

RESULTS

We identified 2 novel genes associated with an increased risk of ischemic stroke: a protein-coding variant in PDE4DIP (rs1778155; odds ratio, 2.15; P = 2.63 × 10(-8)) with an intracellular signal transduction mechanism and in ACOT4 (rs35724886; odds ratio, 2.04; P = 1.24 × 10(-7)) with a fatty acid metabolism; confirmation of PDE4DIP was observed in affected sibpair families with large-vessel stroke subtype and in African Americans. Replication of protein-coding variants in candidate genes was observed for 2 previously reported GWAS associations: ZFHX3 (cardioembolic stroke) and ABCA1 (large-vessel stroke).

CONCLUSIONS AND RELEVANCE

Exome sequencing discovered 2 novel genes and mechanisms, PDE4DIP and ACOT4, associated with increased risk for ischemic stroke. In addition, ZFHX3 and ABCA1 were discovered to have protein-coding variants associated with ischemic stroke. These results suggest that genetic variation in novel pathways contributes to ischemic stroke risk and serves as a target for prediction, prevention, and therapy.

Polygenic Overlap Between C-Reactive Protein, Plasma Lipids, and Alzheimer Disease

BACKGROUND

Epidemiological findings suggest a relationship between Alzheimer disease (AD), inflammation, and dyslipidemia, although the nature of this relationship is not well understood. We investigated whether this phenotypic association arises from a shared genetic basis.

METHODS AND RESULTS

Using summary statistics (P values and odds ratios) from genome-wide association studies of >200 000 individuals, we investigated overlap in single-nucleotide polymorphisms associated with clinically diagnosed AD and C-reactive protein (CRP), triglycerides, and high- and low-density lipoprotein levels. We found up to 50-fold enrichment of AD single-nucleotide polymorphisms for different levels of association with C-reactive protein, low-density lipoprotein, high-density lipoprotein, and triglyceride single-nucleotide polymorphisms using a false discovery rate threshold <0.05. By conditioning on polymorphisms associated with the 4 phenotypes, we identified 55 loci associated with increased AD risk. We then conducted a meta-analysis of these 55 variants across 4 independent AD cohorts (total: n=29 054 AD cases and 114 824 healthy controls) and discovered 2 genome-wide significant variants on chromosome 4 (rs13113697; closest gene, HS3ST1; odds ratio=1.07; 95% confidence interval=1.05-1.11; P=2.86×10(-8)) and chromosome 10 (rs7920721; closest gene, ECHDC3; odds ratio=1.07; 95% confidence interval=1.04-1.11; P=3.38×10(-8)). We also found that gene expression of HS3ST1 and ECHDC3 was altered in AD brains compared with control brains.

CONCLUSIONS

We demonstrate genetic overlap between AD, C-reactive protein, and plasma lipids. By conditioning on the genetic association with the cardiovascular phenotypes, we identify novel AD susceptibility loci, including 2 genome-wide significant variants conferring increased risk for AD.

Strategies to design and analyze targeted sequencing data: cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study

BACKGROUND

Genome-wide association studies have identified thousands of genetic variants that influence a variety of diseases and health-related quantitative traits. However, the causal variants underlying the majority of genetic associations remain unknown. Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study aims to follow up genome-wide association study signals and identify novel associations of the allelic spectrum of identified variants with cardiovascular-related traits.

METHODS AND RESULTS

The study included 4231 participants from 3 CHARGE cohorts: the Atherosclerosis Risk in Communities Study, the Cardiovascular Health Study, and the Framingham Heart Study. We used a case-cohort design in which we selected both a random sample of participants and participants with extreme phenotypes for each of 14 traits. We sequenced and analyzed 77 genomic loci, which had previously been associated with ≥1 of 14 phenotypes. A total of 52 736 variants were characterized by sequencing and passed our stringent quality control criteria. For common variants (minor allele frequency ≥1%), we performed unweighted regression analyses to obtain P values for associations and weighted regression analyses to obtain effect estimates that accounted for the sampling design. For rare variants, we applied 2 approaches: collapsed aggregate statistics and joint analysis of variants using the sequence kernel association test.

CONCLUSIONS

We sequenced 77 genomic loci in participants from 3 cohorts. We established a set of filters to identify high-quality variants and implemented statistical and bioinformatics strategies to analyze the sequence data and identify potentially functional variants within genome-wide association study loci.

Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson's disease

We conducted a meta-analysis of Parkinson's disease genome-wide association studies using a common set of 7,893,274 variants across 13,708 cases and 95,282 controls. Twenty-six loci were identified as having genome-wide significant association; these and 6 additional previously reported loci were then tested in an independent set of 5,353 cases and 5,551 controls. Of the 32 tested SNPs, 24 replicated, including 6 newly identified loci. Conditional analyses within loci showed that four loci, including GBA, GAK-DGKQ, SNCA and the HLA region, contain a secondary independent risk variant. In total, we identified and replicated 28 independent risk variants for Parkinson's disease across 24 loci. Although the effect of each individual locus was small, risk profile analysis showed substantial cumulative risk in a comparison of the highest and lowest quintiles of genetic risk (odds ratio (OR) = 3.31, 95% confidence interval (CI) = 2.55-4.30; P = 2 × 10(-16)). We also show six risk loci associated with proximal gene expression or DNA methylation.

Gene-wide analysis detects two new susceptibility genes for Alzheimer's disease

BACKGROUND

Alzheimer's disease is a common debilitating dementia with known heritability, for which 20 late onset susceptibility loci have been identified, but more remain to be discovered. This study sought to identify new susceptibility genes, using an alternative gene-wide analytical approach which tests for patterns of association within genes, in the powerful genome-wide association dataset of the International Genomics of Alzheimer's Project Consortium, comprising over 7 m genotypes from 25,580 Alzheimer's cases and 48,466 controls.

PRINCIPAL FINDINGS

In addition to earlier reported genes, we detected genome-wide significant loci on chromosomes 8 (TP53INP1, p = 1.4×10-6) and 14 (IGHV1-67 p = 7.9×10-8) which indexed novel susceptibility loci.

SIGNIFICANCE

The additional genes identified in this study, have an array of functions previously implicated in Alzheimer's disease, including aspects of energy metabolism, protein degradation and the immune system and add further weight to these pathways as potential therapeutic targets in Alzheimer's disease.

Serum brain-derived neurotrophic factor and the risk for dementia: the Framingham Heart Study

IMPORTANCE

In animal studies, brain-derived neurotrophic factor (BDNF) has been shown to impact neuronal survival and function and improve synaptic plasticity and long-term memory. Circulating BDNF levels increase with physical activity and caloric restriction, thus BDNF may mediate some of the observed associations between lifestyle and the risk for dementia. Some prior studies showed lower circulating BDNF in persons with Alzheimer disease (AD) compared with control participants; however, it remains uncertain whether reduced levels precede dementia onset.

OBJECTIVE

To examine whether higher serum BDNF levels in cognitively healthy adults protect against the future risk for dementia and AD and to identify potential modifiers of this association.

DESIGN, SETTING, AND PARTICIPANTS

Framingham Study original and offspring participants were followed up from 1992 and 1998, respectively, for up to 10 years. We used Cox models to relate BDNF levels to the risk for dementia and AD and adjusted for potential confounders. We also ran sensitivity analyses stratified by sex, age, and education, as well as related BDNF genetic variants to AD risk. This community-based, prospective cohort study involved 2131 dementia-free participants aged 60 years and older (mean [SD] age, 72 [7] years; 56% women).

MAIN OUTCOMES AND MEASURES

Ten-year incidence of dementia and AD.

RESULTS

During follow-up, 140 participants developed dementia, 117 of whom had AD. Controlling for age and sex, each standard-deviation increment in BDNF was associated with a 33% lower risk for dementia and AD (P = .006 and P = .01, respectively) and these associations persisted after additional adjustments. Compared with the bottom quintile, BDNF levels in the top quintile were associated with less than half the risk for dementia and AD (hazard ratio, 0.49; 95% CI, 0.28-0.85; P = .01; and hazard ratio, 0.46; 95% CI, 0.24-0.86; P = .02, respectively). These associations were apparent only among women, persons aged 80 years and older, and those with college degrees (hazard ratios for AD: 0.65, [95% CI, 0.50-0.85], P = .001; 0.63 [95% CI, 0.47-0.85], P = .002; and 0.27 [95% CI, 0.11-0.65], P = .003, respectively). Brain-derived neurotrophic factor genetic variants were not associated with AD risk.

CONCLUSIONS AND RELEVANCE

Higher serum BDNF levels may protect against future occurrence of dementia and AD. Our findings suggest a role for BDNF in the biology and possibly in the prevention of dementia and AD, especially in select subgroups of women and older and more highly educated persons.