Single-nucleus RNA-sequencing of autosomal dominant Alzheimer disease and risk variant carriers

Genetic studies of Alzheimer disease (AD) have prioritized variants in genes related to the amyloid cascade, lipid metabolism, and neuroimmune modulation. However, the cell-specific effect of variants in these genes is not fully understood. Here, we perform single-nucleus RNA-sequencing (snRNA-seq) on nearly 300,000 nuclei from the parietal cortex of AD autosomal dominant (APP and PSEN1) and risk-modifying variant (APOE, TREM2 and MS4A) carriers. Within individual cell types, we capture genes commonly dysregulated across variant groups. However, specific transcriptional states are more prevalent within variant carriers. TREM2 oligodendrocytes show a dysregulated autophagy-lysosomal pathway, MS4A microglia have dysregulated complement cascade genes, and APOEε4 inhibitory neurons display signs of ferroptosis. All cell types have enriched states in autosomal dominant carriers. We leverage differential expression and single-nucleus ATAC-seq to map GWAS signals to effector cell types including the NCK2 signal to neurons in addition to the initially proposed microglia. Overall, our results provide insights into the transcriptional diversity resulting from AD genetic architecture and cellular heterogeneity. The data can be explored on the online browser ( http://web.hararilab.org/SNARE/ ).

Chi3l1/YKL-40 is controlled by the astrocyte circadian clock and regulates neuroinflammation and Alzheimer's disease pathogenesis

Regulation of glial activation and neuroinflammation are critical factors in the pathogenesis of Alzheimer's disease (AD). YKL-40, a primarily astrocytic protein encoded by the gene Chi3l1, is a widely studied cerebrospinal fluid biomarker that increases with aging and early in AD. However, the function of Chi3l1/YKL-40 in AD is unknown. In a cohort of patients with AD, we observed that a variant in the human CHI3L1 gene, which results in decreased CSF YKL-40 expression, was associated with slower AD progression. At baseline, Chi3l1 deletion in mice had no effect on astrocyte activation while modestly promoting microglial activation. In a mouse APP/PS1 model of AD, Chi3l1 deletion decreased amyloid plaque burden and increased periplaque expression of the microglial lysosomal marker CD68, suggesting that Chi3l1 may suppress glial phagocytic activation and promote amyloid accumulation. Accordingly, Chi3l1 knockdown increased phagocytosis of zymosan particles and of β-amyloid peptide in both astrocytes and microglia in vitro. We further observed that expression of Chi3l1 is regulated by the circadian clock, as deletion of the core clock proteins BMAL1 or CLOCK/NPAS2 strongly suppresses basal Chi3l1 expression, whereas deletion of the negative clock regulators PER1/PER2 increased Chi3l1 expression. Basal Chi3l1 mRNA was nonrhythmic because of a long mRNA half-life in astrocytes. However, inflammatory induction of Chi3l1 was gated by the clock. Our findings reveal Chi3l1/YKL-40 as a modulator of glial phagocytic activation and AD pathogenesis in both mice and humans and suggest that the astrocyte circadian clock regulates inflammatory Chi3l1 induction.

The TMEM106B FTLD-protective variant, rs1990621, is also associated with increased neuronal proportion

Apart from amyloid β deposition and tau neurofibrillary tangles, Alzheimer's disease (AD) is a neurodegenerative disorder characterized by neuronal loss and astrocytosis in the cerebral cortex. The goal of this study is to investigate genetic factors associated with the neuronal proportion in health and disease. To identify cell-autonomous genetic variants associated with neuronal proportion in cortical tissues, we inferred cellular population structure from bulk RNA-Seq derived from 1536 individuals. We identified the variant rs1990621 located in the TMEM106B gene region as significantly associated with neuronal proportion (p value = 6.40 × 10-07) and replicated this finding in an independent dataset (p value = 7.41 × 10-04) surpassing the genome-wide threshold in the meta-analysis (p value = 9.42 × 10-09). This variant is in high LD with the TMEM106B non-synonymous variant p.T185S (rs3173615; r2 = 0.98) which was previously identified as a protective variant for frontotemporal lobar degeneration (FTLD). We stratified the samples by disease status, and discovered that this variant modulates neuronal proportion not only in AD cases, but also several neurodegenerative diseases and in elderly cognitively healthy controls. Furthermore, we did not find a significant association in younger controls or schizophrenia patients, suggesting that this variant might increase neuronal survival or confer resilience to the neurodegenerative process. The single variant and gene-based analyses also identified an overall genetic association between neuronal proportion, AD and FTLD risk. These results suggest that common pathways are implicated in these neurodegenerative diseases, that implicate neuronal survival. In summary, we identified a protective variant in the TMEM106B gene that may have a neuronal protection effect against general aging, independent of disease status, which could help elucidate the relationship between aging and neuronal survival in the presence or absence of neurodegenerative disorders. Our findings suggest that TMEM106B could be a potential target for neuronal protection therapies to ameliorate cognitive and functional deficits.

An atlas of cortical circular RNA expression in Alzheimer disease brains demonstrates clinical and pathological associations

We generated parietal cortex RNA-seq data from individuals with and without Alzheimer disease (AD; n_control = 13; n_AD = 83) from the Knight-ADRC. Using this and an independent (MSBB) AD RNA-seq dataset, we quantified cortical circular RNA (circRNA) expression in the context of AD. We identified significant associations between circRNA expression and AD diagnosis, clinical dementia severity, and neuropathological severity. We demonstrated that a majority of circRNA AD-associations are independent from changes in cognate linear mRNA expression or brain cell-type proportions. We provided evidence for circRNA expression changes occurring early in pre-symptomatic AD, and in autosomal dominant AD. We also observed AD-associated circRNAs co-expressing with known AD genes. Finally, we identified potential microRNA binding sites in AD-associated circRNAs for microRNAs predicted to target AD genes. Together, these results highlight the importance of analyzing non-linear RNAs and support future studies exploring the potential roles of circRNAs in AD pathogenesis.

The MS4A gene cluster is a key modulator of soluble TREM2 and Alzheimer's disease risk

Soluble triggering receptor expressed on myeloid cells 2 (sTREM2) in cerebrospinal fluid (CSF) has been associated with Alzheimer's disease (AD). TREM2 plays a critical role in microglial activation, survival, and phagocytosis; however, the pathophysiological role of sTREM2 in AD is not well understood. Understanding the role of sTREM2 in AD may reveal new pathological mechanisms and lead to the identification of therapeutic targets. We performed a genome-wide association study (GWAS) to identify genetic modifiers of CSF sTREM2 obtained from the Alzheimer's Disease Neuroimaging Initiative. Common variants in the membrane-spanning 4-domains subfamily A (MS4A) gene region were associated with CSF sTREM2 concentrations (rs1582763; P = 1.15 × 10-15); this was replicated in independent datasets. The variants associated with increased CSF sTREM2 concentrations were associated with reduced AD risk and delayed age at onset of disease. The single-nucleotide polymorphism rs1582763 modified expression of the MS4A4A and MS4A6A genes in multiple tissues, suggesting that one or both of these genes are important for modulating sTREM2 production. Using human macrophages as a proxy for microglia, we found that MS4A4A and TREM2 colocalized on lipid rafts at the plasma membrane, that sTREM2 increased with MS4A4A overexpression, and that silencing of MS4A4A reduced sTREM2 production. These genetic, molecular, and cellular findings suggest that MS4A4A modulates sTREM2. These findings also provide a mechanistic explanation for the original GWAS signal in the MS4A locus for AD risk and indicate that TREM2 may be involved in AD pathogenesis not only in TREM2 risk-variant carriers but also in those with sporadic disease.

A single-nuclei RNA sequencing study of Mendelian and sporadic AD in the human brain

Background

Alzheimer’s disease (AD) is the most common form of dementia. This neurodegenerative disorder is associated with neuronal death and gliosis heavily impacting the cerebral cortex. AD has a substantial but heterogeneous genetic component, presenting both Mendelian and complex genetic architectures. Using bulk RNA-seq from the parietal lobes and deconvolution methods, we previously reported that brains exhibiting different AD genetic architecture exhibit different cellular proportions. Here, we sought to directly investigate AD brain changes in cell proportion and gene expression using single-cell resolution.

Methods

We generated unsorted single-nuclei RNA sequencing data from brain tissue. We leveraged the tissue donated from a carrier of a Mendelian genetic mutation, PSEN1 p.A79V, and two family members who suffer from sporadic AD, but do not carry any autosomal mutations. We evaluated alternative alignment approaches to maximize the titer of reads, genes, and cells with high quality. In addition, we employed distinct clustering strategies to determine the best approach to identify cell clusters that reveal neuronal and glial cell types and avoid artifacts such as sample and batch effects. We propose an approach to cluster cells that reduces biases and enable further analyses.

Results

We identified distinct types of neurons, both excitatory and inhibitory, and glial cells, including astrocytes, oligodendrocytes, and microglia, among others. In particular, we identified a reduced proportion of excitatory neurons in the Mendelian mutation carrier, but a similar distribution of inhibitory neurons. Furthermore, we investigated whether single-nuclei RNA-seq from the human brains recapitulate the expression profile of disease-associated microglia (DAM) discovered in mouse models. We also determined that when analyzing human single-nuclei data, it is critical to control for biases introduced by donor-specific expression profiles.

Conclusion

We propose a collection of best practices to generate a highly detailed molecular cell atlas of highly informative frozen tissue stored in brain banks. Importantly, we have developed a new web application to make this unique single-nuclei molecular atlas publicly available.

Electronic supplementary material

The online version of this article (10.1186/s13195-019-0524-x) contains supplementary material, which is available to authorized users.

TREM2 brain transcript-specific studies in AD and TREM2 mutation carriers

BACKGROUND

Low frequency coding variants in TREM2 are associated with Alzheimer disease (AD) risk and cerebrospinal fluid (CSF) TREM2 protein levels are different between AD cases and controls. Similarly, TREM2 risk variant carriers also exhibit differential CSF TREM2 levels. TREM2 has three different alternative transcripts, but most of the functional studies only model the longest transcript. No studies have analyzed TREM2 expression levels or alternative splicing in brains from AD and cognitively normal individuals. We wanted to determine whether there was differential expression of TREM2 in sporadic-AD cases versus AD-TREM2 carriers vs sex- and aged-matched normal controls; and if this differential expression was due to a particular TREM2 transcript.

METHODS

We analyzed RNA-Seq data from parietal lobe brain tissue from AD cases with TREM2 variants (n = 33), AD cases (n = 195) and healthy controls (n = 118), from three independent datasets using Kallisto and the R package tximport to determine the read count for each transcript and quantified transcript abundance as transcripts per million.

RESULTS

The three TREM2 transcripts were expressed in brain cortex in the three datasets. We demonstrate for the first time that the transcript that lacks the transmembrane domain and encodes a soluble form of TREM2 (sTREM2) has an expression level around 60% of the canonical transcript, suggesting that around 25% of the sTREM2 protein levels could be explained by this transcript. We did not observe a difference in the overall TREM2 expression level between cases and controls. However, the isoform which lacks the 5' exon, but includes the transmembrane domain, was significantly lower in TREM2- p.R62H carriers than in AD cases (p = 0.007).

CONCLUSION

Using bulk RNA-Seq data from three different cohorts, we were able to quantify the expression level of the three TREM2 transcripts, demonstrating: (1) all three transcripts of them are highly expressed in the human cortex, (2) that up to 25% of the sTREM2 may be due to the expression of a specific isoform and not TREM2 cleavage; and (3) that TREM2 risk variants do not affect expression levels, suggesting that the effect of the TREM2 variants on CSF levels occurs at post-transcriptional level.

Genetic variants associated with Alzheimer's disease confer different cerebral cortex cell-type population structure

Background

Alzheimer’s disease (AD) is characterized by neuronal loss and astrocytosis in the cerebral cortex. However, the specific effects that pathological mutations and coding variants associated with AD have on the cellular composition of the brain are often ignored.

Methods

We developed and optimized a cell-type-specific expression reference panel and employed digital deconvolution methods to determine brain cellular distribution in three independent transcriptomic studies.

Results

We found that neuronal and astrocyte relative proportions differ between healthy and diseased brains and also among AD cases that carry specific genetic risk variants. Brain carriers of pathogenic mutations in APP, PSEN1, or PSEN2 presented lower neuron and higher astrocyte relative proportions compared to sporadic AD. Similarly, the APOE ε4 allele also showed decreased neuronal and increased astrocyte relative proportions compared to AD non-carriers. In contrast, carriers of variants in TREM2 risk showed a lower degree of neuronal loss compared to matched AD cases in multiple independent studies.

Conclusions

These findings suggest that genetic risk factors associated with AD etiology have a specific imprinting in the cellular composition of AD brains. Our digital deconvolution reference panel provides an enhanced understanding of the fundamental molecular mechanisms underlying neurodegeneration, enabling the analysis of large bulk RNA-sequencing studies for cell composition and suggests that correcting for the cellular structure when performing transcriptomic analysis will lead to novel insights of AD.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0551-4) contains supplementary material, which is available to authorized users.

Evaluation of Gene-Based Family-Based Methods to Detect Novel Genes Associated With Familial Late Onset Alzheimer Disease

Gene-based tests to study the combined effect of rare variants on a particular phenotype have been widely developed for case-control studies, but their evolution and adaptation for family-based studies, especially studies of complex incomplete families, has been slower. In this study, we have performed a practical examination of all the latest gene-based methods available for family-based study designs using both simulated and real datasets. We examined the performance of several collapsing, variance-component, and transmission disequilibrium tests across eight different software packages and 22 models utilizing a cohort of 285 families (N = 1,235) with late-onset Alzheimer disease (LOAD). After a thorough examination of each of these tests, we propose a methodological approach to identify, with high confidence, genes associated with the tested phenotype and we provide recommendations to select the best software and model for family-based gene-based analyses. Additionally, in our dataset, we identified PTK2B, a GWAS candidate gene for sporadic AD, along with six novel genes (CHRD, CLCN2, HDLBP, CPAMD8, NLRP9, and MAS1L) as candidate genes for familial LOAD.

Polygenic risk score of sporadic late-onset Alzheimer's disease reveals a shared architecture with the familial and early-onset forms

OBJECTIVE

To determine whether the extent of overlap of the genetic architecture among the sporadic late-onset Alzheimer's Disease (sLOAD), familial late-onset AD (fLOAD), sporadic early-onset AD (sEOAD), and autosomal dominant early-onset AD (eADAD).

METHODS

Polygenic risk scores (PRSs) were constructed using previously identified 21 genome-wide significant loci for LOAD risk.

RESULTS

We found that there is an overlap in the genetic architecture among sEOAD, fLOAD, and sLOAD. The highest association of the PRS and risk (odds ratio [OR] = 2.27; P = 1.29 × 10-7) was observed in sEOAD, followed by fLOAD (OR = 1.75; P = 1.12 × 10-7) and sLOAD (OR = 1.40; P = 1.21 × 10-3). The PRS was associated with cerebrospinal fluid ptau181-Aβ42 on eADAD (P = 4.36 × 10-2).

CONCLUSION

Our analysis confirms that the genetic factors identified for LOAD modulate risk in sLOAD and fLOAD and also sEOAD cohorts. Specifically, our results suggest that the burden of these risk variants is associated with familial clustering and earlier onset of AD. Although these variants are not associated with risk in the eADAD, they may be modulating age at onset.

Assessment of the Genetic Architecture of Alzheimer's Disease Risk in Rate of Memory Decline

Many genetic studies for Alzheimer's disease (AD) have been focused on the identification of common genetic variants associated with AD risk and not on other aspects of the disease, such as age at onset or rate of dementia progression. There are multiple approaches to untangling the genetic architecture of these phenotypes. We hypothesized that the genetic architecture of rate of progression is different than the risk for developing AD dementia. To test this hypothesis, we used longitudinal clinical data from ADNI and the Knight-ADRC at Washington University, and we calculated PRS (polygenic risk score) based on the IGAP study to compare the genetic architecture of AD risk and dementia progression. Dementia progression was measured by the change of Clinical Dementia Rating Sum of Boxes (CDR)-SB per year. Out of the 21 loci for AD risk, no association with the rate of dementia progression was found. The PRS rate was significantly associated with the rate of dementia progression (β= 0.146, p = 0.03). In the case of rare variants, TREM2 (β= 0.309, p = 0.02) was also associated with the rate of dementia progression. TREM2 variant carriers showed a 23% faster rate of dementia compared with non-variant carriers. In conclusion, our results indicate that the recently identified common and rare variants for AD susceptibility have a limited impact on the rate of dementia progression in AD patients.

Parkinson disease polygenic risk score is associated with Parkinson disease status and age at onset but not with alpha-synuclein cerebrospinal fluid levels

BACKGROUND

The genetic architecture of Parkinson's Disease (PD) is complex and not completely understood. Multiple genetic studies to date have identified multiple causal genes and risk loci. Nevertheless, most of the expected genetic heritability remains unexplained. Polygenic risk scores (PRS) may provide greater statistical power and inform about the genetic architecture of multiple phenotypes. The aim of this study was to test the association between PRS and PD risk, age at onset and cerebrospinal fluid (CSF) biomarkers (α-synuclein, Aβ_1-42, t-tau and p-tau).

METHODS

The weighted PRS was created using the genome-wide loci from Nalls et al., 2014 PD GWAs meta-analysis. The PRS was tested for association with PD status, age at onset and CSF biomarker levels in 829 cases and 432 controls of European ancestry.

RESULTS

The PRS was associated with PD status (p = 5.83×10^-08) and age at onset (p = 5.70×10^-07). The CSF t-tau levels showed a nominal association with the PRS (p = 0.02). However, CSF α-synuclein, amyloid beta and phosphorylated tau were not found to be associated with the PRS.

CONCLUSION

Our study suggests that there is an overlap in the genetic architecture of PD risk and onset, although the different loci present different weights for those phenotypes. In our dataset we found a marginal association of the PRS with CSF t-tau but not with α-synuclein CSF levels, suggesting that the genetic architecture for the CSF biomarker levels is different from that of PD risk.

ApoE4 markedly exacerbates tau-mediated neurodegeneration in a mouse model of tauopathy

APOE4 is the strongest genetic risk factor for late-onset Alzheimer disease. ApoE4 increases brain amyloid-β pathology relative to other ApoE isoforms. However, whether APOE independently influences tau pathology, the other major proteinopathy of Alzheimer disease and other tauopathies, or tau-mediated neurodegeneration, is not clear. By generating P301S tau transgenic mice on either a human ApoE knock-in (KI) or ApoE knockout (KO) background, here we show that P301S/E4 mice have significantly higher tau levels in the brain and a greater extent of somatodendritic tau redistribution by three months of age compared with P301S/E2, P301S/E3, and P301S/EKO mice. By nine months of age, P301S mice with different ApoE genotypes display distinct phosphorylated tau protein (p-tau) staining patterns. P301S/E4 mice develop markedly more brain atrophy and neuroinflammation than P301S/E2 and P301S/E3 mice, whereas P301S/EKO mice are largely protected from these changes. In vitro, E4-expressing microglia exhibit higher innate immune reactivity after lipopolysaccharide treatment. Co-culturing P301S tau-expressing neurons with E4-expressing mixed glia results in a significantly higher level of tumour-necrosis factor-α (TNF-α) secretion and markedly reduced neuronal viability compared with neuron/E2 and neuron/E3 co-cultures. Neurons co-cultured with EKO glia showed the greatest viability with the lowest level of secreted TNF-α. Treatment of P301S neurons with recombinant ApoE (E2, E3, E4) also leads to some neuronal damage and death compared with the absence of ApoE, with ApoE4 exacerbating the effect. In individuals with a sporadic primary tauopathy, the presence of an ε4 allele is associated with more severe regional neurodegeneration. In individuals who are positive for amyloid-β pathology with symptomatic Alzheimer disease who usually have tau pathology, ε4-carriers demonstrate greater rates of disease progression. Our results demonstrate that ApoE affects tau pathogenesis, neuroinflammation, and tau-mediated neurodegeneration independently of amyloid-β pathology. ApoE4 exerts a 'toxic' gain of function whereas the absence of ApoE is protective.

A common haplotype lowers PU.1 expression in myeloid cells and delays onset of Alzheimer's disease

A genome-wide survival analysis of 14,406 Alzheimer's disease (AD) cases and 25,849 controls identified eight previously reported AD risk loci and 14 novel loci associated with age at onset. Linkage disequilibrium score regression of 220 cell types implicated the regulation of myeloid gene expression in AD risk. The minor allele of rs1057233 (G), within the previously reported CELF1 AD risk locus, showed association with delayed AD onset and lower expression of SPI1 in monocytes and macrophages. SPI1 encodes PU.1, a transcription factor critical for myeloid cell development and function. AD heritability was enriched within the PU.1 cistrome, implicating a myeloid PU.1 target gene network in AD. Finally, experimentally altered PU.1 levels affected the expression of mouse orthologs of many AD risk genes and the phagocytic activity of mouse microglial cells. Our results suggest that lower SPI1 expression reduces AD risk by regulating myeloid gene expression and cell function.

Genome-wide association study identifies four novel loci associated with Alzheimer's endophenotypes and disease modifiers

More than 20 genetic loci have been associated with risk for Alzheimer's disease (AD), but reported genome-wide significant loci do not account for all the estimated heritability and provide little information about underlying biological mechanisms. Genetic studies using intermediate quantitative traits such as biomarkers, or endophenotypes, benefit from increased statistical power to identify variants that may not pass the stringent multiple test correction in case-control studies. Endophenotypes also contain additional information helpful for identifying variants and genes associated with other aspects of disease, such as rate of progression or onset, and provide context to interpret the results from genome-wide association studies (GWAS). We conducted GWAS of amyloid beta (Aβ42), tau, and phosphorylated tau (ptau181) levels in cerebrospinal fluid (CSF) from 3146 participants across nine studies to identify novel variants associated with AD. Five genome-wide significant loci (two novel) were associated with ptau181, including loci that have also been associated with AD risk or brain-related phenotypes. Two novel loci associated with Aβ42 near GLIS1 on 1p32.3 (β = -0.059, P = 2.08 × 10-8) and within SERPINB1 on 6p25 (β = -0.025, P = 1.72 × 10-8) were also associated with AD risk (GLIS1: OR = 1.105, P = 3.43 × 10-2), disease progression (GLIS1: β = 0.277, P = 1.92 × 10-2), and age at onset (SERPINB1: β = 0.043, P = 4.62 × 10-3). Bioinformatics indicate that the intronic SERPINB1 variant (rs316341) affects expression of SERPINB1 in various tissues, including the hippocampus, suggesting that SERPINB1 influences AD through an Aβ-associated mechanism. Analyses of known AD risk loci suggest CLU and FERMT2 may influence CSF Aβ42 (P = 0.001 and P = 0.009, respectively) and the INPP5D locus may affect ptau181 levels (P = 0.009); larger studies are necessary to verify these results. Together the findings from this study can be used to inform future AD studies.

Chitinase-3-like 1 protein (CHI3L1) locus influences cerebrospinal fluid levels of YKL-40

Background

Alzheimer’s disease (AD) pathology appears several years before clinical symptoms, so identifying ways to detect individuals in the preclinical stage is imperative. The cerebrospinal fluid (CSF) Tau/Aβ₄₂ ratio is currently the best known predictor of AD status and cognitive decline, and the ratio of CSF levels of chitinase-3-like 1 protein (CHI3L1, YKL-40) and amyloid beta (Aβ₄₂) were reported as predictive, but individual variability and group overlap inhibits their utility for individual diagnosis making it necessary to find ways to improve sensitivity of these biomarkers.

Methods

We used linear regression to identify genetic loci associated with CSF YKL-40 levels in 379 individuals (80 cognitively impaired and 299 cognitively normal) from the Charles F and Joanne Knight Alzheimer’s Disease Research Center. We tested correlations between YKL-40 and CSF Tau/Aβ₄₂ ratio, Aβ₄₂, tau, and phosphorylated tau (ptau₁₈₁). We used studentized residuals from a linear regression model of the log-transformed, standardized protein levels and the additive reference allele counts from the most significant locus to adjust YKL-40 values and tested the differences in correlations with CSF Tau/Aβ₄₂ ratio, Aβ₄₂, tau, and ptau₁₈₁.

Results

We found that genetic variants on the CH13L1 locus were significantly associated with CSF YKL-40 levels, but not AD risk, age at onset, or disease progression. The most significant variant is a reported expression quantitative trait locus for CHI3L1, the gene which encodes YKL-40, and explained 12.74 % of the variance in CSF YKL-40 in our study. YKL-40 was positively correlated with ptau₁₈₁ (r = 0.521) and the strength of the correlation significantly increased with the addition of genetic information (r = 0.573, p = 0.006).

Conclusions

CSF YKL-40 levels are likely a biomarker for AD, but we found no evidence that they are an AD endophenotype. YKL-40 levels are highly regulated by genetic variation, and by including genetic information the strength of the correlation between YKL-40 and ptau₁₈₁ levels is significantly improved. Our results suggest that studies of potential biomarkers may benefit from including genetic information.

Electronic supplementary material

The online version of this article (doi:10.1186/s12883-016-0742-9) contains supplementary material, which is available to authorized users.

Genetic studies of plasma analytes identify novel potential biomarkers for several complex traits

Genome-wide association studies of 146 plasma protein levels in 818 individuals revealed 56 genome-wide significant associations (28 novel) with 47 analytes. Loci associated with plasma levels of 39 proteins tested have been previously associated with various complex traits such as heart disease, inflammatory bowel disease, Type 2 diabetes and multiple sclerosis. These data suggest that these plasma protein levels may constitute informative endophenotypes for these complex traits. We found three potential pleiotropic genes: ABO for plasma SELE and ACE levels, FUT2 for CA19-9 and CEA plasma levels and APOE for ApoE and CRP levels. We also found multiple independent signals in loci associated with plasma levels of ApoH, CA19-9, FetuinA, IL6r and LPa. Our study highlights the power of biological traits for genetic studies to identify genetic variants influencing clinically relevant traits, potential pleiotropic effects and complex disease associations in the same locus.

Role of ABCA7 loss-of-function variant in Alzheimer's disease: a replication study in European-Americans

INTRODUCTION

A recent study found a significant increase of ABCA7 loss-of-function variants in Alzheimer's disease (AD) cases compared to controls. Some variants were located on noncoding regions, but it was demonstrated that they affect splicing. Here, we try to replicate the association between AD risk and ABCA7 loss-of-function variants at both the single-variant and gene level in a large and well-characterized European American dataset.

METHODS

We genotyped the GWAS common variant and four rare variants previously reported for ABCA7 in 3476 European-Americans.

RESULTS

We were not able to replicate the association at the single-variant level, likely due to a lower effect size on the European American population which led to limited statistical power. However, we did replicate the association at the gene level; we found a significant enrichment of ABCA7 loss-of-function variants in AD cases compared to controls (P = 0.0388; odds ratio =1.54). We also confirmed that the association of the loss-of-function variants is independent of the previously reported genome-wide association study signal.

CONCLUSIONS

Although the effect size for the association of ABCA7 loss-of-function variants with AD risk is lower in our study (odds ratio = 1.54) compared to the original report (odds ratio = 2.2), the replication of the findings of the original report provides a stronger foundation for future functional applications. The data indicate that different independent signals that modify risk for complex traits may exist on the same locus. Additionally, our results suggest that replication of rare-variant studies should be performed at the gene level rather than focusing on a single variant.

Alzheimer's disease: rare variants with large effect sizes

Recent advances in sequencing technology and novel genotyping arrays (focused on low-frequency and coding variants) have made it possible to identify novel coding variants with large effect sizes and also novel genes (TREM2, PLD3, UNC5C, and AKAP9) associated with Alzheimer's disease (AD) risk. The major advantages of these studies over the classic genome-wide association studies (GWAS) include the identification of the functional variant and the gene-driven association. In addition to the large effect size, these studies make it possible to model these variants and genes using cell and animal systems. On the other hand, the underlying population-variability of these very low allele frequency variants poses a great challenge to replicating results. Studies that include very large datasets (>10,000 cases and controls) and combine sequencing and genotyping approaches will lead to the identification of novel genes for Alzheimer's disease.

Transethnic meta-analysis suggests genetic variation in the HEME pathway influences potassium response in patients treated with hydrochlorothiazide

Hypokalemia is a recognized adverse effect of thiazide diuretic treatment. This phenomenon, which may impair insulin secretion, has been suggested to be a reason for the adverse effects on glucose metabolism associated with thiazide diuretic treatment of hypertension. However, the mechanisms underlying thiazide diuretic-induced hypokalemia are not well understood. In an effort to identify genes or genomic regions associated with potassium response to hydrochlorothiazide, without a priori knowledge of biologic effects, we performed a genome-wide association study and a multiethnic meta-analysis in 718 European- and African-American hypertensive participants from two different pharmacogenetic studies. Single-nucleotide polymorphisms rs10845697 (Bayes factor=5.560) on chromosome 12, near to the HEME binding protein 1 gene, and rs11135740 (Bayes factor=5.258) on chromosome 8, near to the Mitoferrin-1 gene, reached genome-wide association study significance (Bayes factor >5). These results, if replicated, suggest a novel mechanism involving effects of genes in the HEME pathway influencing hydrochlorothiazide-induced renal potassium loss.

Genome-wide association analyses suggest NELL1 influences adverse metabolic response to HCTZ in African Americans

Hydrochlorothiazide (HCTZ) is one of the most widely prescribed antihypertensive medications. Although it is well known that HCTZ is associated with hyperglycemia and hypertriglyceridemia, the mechanisms underlying these adverse effects are not well understood. We performed a genome-wide association study and meta-analysis of the change in fasting plasma glucose and triglycerides in response to HCTZ from two different clinical trials: the Pharmacogenomic Evaluation of Antihypertensive Responses and the Genetic Epidemiology of Responses to Antihypertensive studies. Two single-nucleotide polymorphisms (rs12279250 and rs4319515 (r²=0.73)), located at 11p15.1 in the NELL1 gene, achieved genome-wide significance for association with change in fasting plasma triglycerides in African Americans, whereby each variant allele was associated with a 28 mg dl⁻¹ increase in the change in triglycerides. NELL1 encodes a cytoplasmic protein that contains epidermal growth factor-like repeats and has been shown to represses adipogenic differentiation. These findings may represent a novel mechanism underlying HCTZ-induced adverse metabolic effects.