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/ ).

Metabolomic and lipidomic signatures in autosomal dominant and late-onset Alzheimer's disease brains.

INTRODUCTION: The identification of multiple genetic risk factors for Alzheimer's disease (AD) suggests that many pathways contribute to AD onset and progression. However, the metabolomic and lipidomic profiles in carriers of distinct genetic risk factors are not fully understood. The metabolome can provide a direct image of dysregulated pathways in the brain. METHODS: We interrogated metabolomic signatures in the AD brain, including carriers of pathogenic variants in APP, PSEN1, and PSEN2 (autosomal dominant AD; ADAD), APOE ɛ4, and TREM2 risk variant carriers, and sporadic AD (sAD). RESULTS: We identified 133 unique and shared metabolites associated with ADAD, TREM2, and sAD. We identified a signature of 16 metabolites significantly altered between groups and associated with AD duration. DISCUSSION: AD genetic variants show distinct metabolic perturbations. Investigation of these metabolites may provide greater insight into the etiology of AD and its impact on clinical presentation. HIGHLIGHTS: APP/PSEN1/PSEN2 and TREM2 variant carriers show distinct metabolic changes. A total of 133 metabolites were differentially abundant in AD genetic groups. ß-citrylglutamate is differentially abundant in autosomal dominant, TREM2, and sporadic AD. A 16-metabolite profile shows differences between Alzheimer's disease (AD) genetic groups. The identified metabolic profile is associated with duration of disease.

Circular RNA detection identifies circPSEN1 alterations in brain specific to autosomal dominant Alzheimer's disease.

BACKGROUND: Autosomal-dominant Alzheimer's disease (ADAD) is caused by pathogenic mutations in APP, PSEN1, and PSEN2, which usually lead to an early age at onset (< 65). Circular RNAs are a family of non-coding RNAs highly expressed in the nervous system and especially in synapses. We aimed to investigate differences in brain gene expression of linear and circular transcripts from the three ADAD genes in controls, sporadic AD, and ADAD. METHODS: We obtained and sequenced RNA from brain cortex using standard protocols. Linear counts were obtained using the TOPMed pipeline; circular counts, using python package DCC. After stringent quality control (QC), we obtained the counts for PSEN1, PSEN2 and APP genes. Only circPSEN1 passed QC. We used DESeq2 to compare the counts across groups, correcting for biological and technical variables. Finally, we performed in-silico functional analyses using the Circular RNA interactome website and DIANA mirPath software. RESULTS: Our results show significant differences in gene counts of circPSEN1 in ADAD individuals, when compared to sporadic AD and controls (ADAD = 21, AD = 253, Controls = 23-ADADvsCO: log2FC = 0.794, p = 1.63 × 10-04, ADADvsAD: log2FC = 0.602, p = 8.22 × 10-04). The high gene counts are contributed by two circPSEN1 species (hsa_circ_0008521 and hsa_circ_0003848). No significant differences were observed in linear PSEN1 gene expression between cases and controls, indicating that this finding is specific to the circular forms. In addition, the high circPSEN1 levels do not seem to be specific to PSEN1 mutation carriers; the counts are also elevated in APP and PSEN2 mutation carriers. In-silico functional analyses suggest that circPSEN1 is involved in several pathways such as axon guidance (p = 3.39 × 10-07), hippo signaling pathway (p = 7.38 × 10-07), lysine degradation (p = 2.48 × 10-05) or Wnt signaling pathway (p = 5.58 × 10-04) among other KEGG pathways. Additionally, circPSEN1 counts were able to discriminate ADAD from sporadic AD and controls with an AUC above 0.70. CONCLUSIONS: Our findings show the differential expression of circPSEN1 is increased in ADAD. Given the biological function previously ascribed to circular RNAs and the results of our in-silico analyses, we hypothesize that this finding might be related to neuroinflammatory events that lead or that are caused by the accumulation of amyloid-beta.

Genomic atlas of the proteome from brain, CSF and plasma prioritizes proteins implicated in neurological disorders.

Understanding the tissue-specific genetic controls of protein levels is essential to uncover mechanisms of post-transcriptional gene regulation. In this study, we generated a genomic atlas of protein levels in three tissues relevant to neurological disorders (brain, cerebrospinal fluid and plasma) by profiling thousands of proteins from participants with and without Alzheimer's disease. We identified 274, 127 and 32 protein quantitative trait loci (pQTLs) for cerebrospinal fluid, plasma and brain, respectively. cis-pQTLs were more likely to be tissue shared, but trans-pQTLs tended to be tissue specific. Between 48.0% and 76.6% of pQTLs did not co-localize with expression, splicing, DNA methylation or histone acetylation QTLs. Using Mendelian randomization, we nominated proteins implicated in neurological diseases, including Alzheimer's disease, Parkinson's disease and stroke. This first multi-tissue study will be instrumental to map signals from genome-wide association studies onto functional genes, to discover pathways and to identify drug targets for neurological diseases.

Long runs of homozygosity are associated with Alzheimer's disease.

Long runs of homozygosity (ROH) are contiguous stretches of homozygous genotypes, which are a footprint of inbreeding and recessive inheritance. The presence of recessive loci is suggested for Alzheimer's disease (AD); however, their search has been poorly assessed to date. To investigate homozygosity in AD, here we performed a fine-scale ROH analysis using 10 independent cohorts of European ancestry (11,919 AD cases and 9181 controls.) We detected an increase of homozygosity in AD cases compared to controls [ßAVROH (CI 95%) = 0.070 (0.037-0.104); P = 3.91 × 10-5; ßFROH (CI95%) = 0.043 (0.009-0.076); P = 0.013]. ROHs increasing the risk of AD (OR > 1) were significantly overrepresented compared to ROHs increasing protection (p < 2.20 × 10-16). A significant ROH association with AD risk was detected upstream the HS3ST1 locus (chr4:11,189,482‒11,305,456), (ß (CI 95%) = 1.09 (0.48 ‒ 1.48), p value = 9.03 × 10-4), previously related to AD. Next, to search for recessive candidate variants in ROHs, we constructed a homozygosity map of inbred AD cases extracted from an outbred population and explored ROH regions in whole-exome sequencing data (N = 1449). We detected a candidate marker, rs117458494, mapped in the SPON1 locus, which has been previously associated with amyloid metabolism. Here, we provide a research framework to look for recessive variants in AD using outbred populations. Our results showed that AD cases have enriched homozygosity, suggesting that recessive effects may explain a proportion of AD heritability.

Exome-wide rare variant analysis in familial essential tremor.

INTRODUCTION: Essential tremor (ET) is one of the most common movement disorders. Despite its high prevalence and heritability, its genetic etiology remains elusive with only a few susceptibility genes identified and poorly replicated. Our aim was to find novel candidate genes involved in ET predisposition through whole exome sequencing. METHODS: We studied eight multigenerational families (N = 40 individuals) with an autosomal-dominant inheritance using a comprehensive strategy combining whole exome sequencing followed by case-control association testing of prioritized variants in a separate cohort comprising 521 ET cases and 596 controls. We further performed gene-based burden analyses in an additional dataset comprising 789 ET patients and 770 healthy individuals to investigate whether there was an enrichment of rare deleterious variants within our candidate genes. RESULTS: Fifteen variants co-segregated with disease status in at least one of the families, among which rs749875462 in CCDC183, rs535864157 in MMP10 and rs114285050 in GPR151 showed a nominal association with ET. However, we found no significant enrichment of rare variants within these genes in cases compared with controls. Interestingly, MMP10 protein is involved in the inflammatory response to neuronal damage and has been previously associated with other neurological disorders. CONCLUSIONS: We prioritized a set of promising genes, especially MMP10, for further genetic and functional studies in ET. Our study suggests that rare deleterious coding variants that markedly increase susceptibility to ET are likely to be found in many genes. Future studies are needed to replicate and further infer biological mechanisms and potential disease causality for our identified genes.

Functional genomic analyses uncover APOE-mediated regulation of brain and cerebrospinal fluid beta-amyloid levels in Parkinson disease.

Alpha-synuclein is the main protein component of Lewy bodies, the pathological hallmark of Parkinson's disease. However, genetic modifiers of cerebrospinal fluid (CSF) alpha-synuclein levels remain unknown. The use of CSF levels of amyloid beta1-42, total tau, and phosphorylated tau181 as quantitative traits in genetic studies have provided novel insights into Alzheimer's disease pathophysiology. A systematic study of the genomic architecture of CSF biomarkers in Parkinson's disease has not yet been conducted. Here, genome-wide association studies of CSF biomarker levels in a cohort of individuals with Parkinson's disease and controls (N = 1960) were performed. PD cases exhibited significantly lower CSF biomarker levels compared to controls. A SNP, proxy for APOE ε4, was associated with CSF amyloid beta1-42 levels (effect = - 0.5, p = 9.2 × 10-19). No genome-wide loci associated with CSF alpha-synuclein, total tau, or phosphorylated tau181 levels were identified in PD cohorts. Polygenic risk score constructed using the latest Parkinson's disease risk meta-analysis were associated with Parkinson's disease status (p = 0.035) and the genomic architecture of CSF amyloid beta1-42 (R2 = 2.29%; p = 2.5 × 10-11). Individuals with higher polygenic risk scores for PD risk presented with lower CSF amyloid beta1-42 levels (p = 7.3 × 10-04). Two-sample Mendelian Randomization revealed that CSF amyloid beta1-42 plays a role in Parkinson's disease (p = 1.4 × 10-05) and age at onset (p = 7.6 × 10-06), an effect mainly mediated by variants in the APOE locus. In a subset of PD samples, the APOE ε4 allele was associated with significantly lower levels of CSF amyloid beta1-42 (p = 3.8 × 10-06), higher mean cortical binding potentials (p = 5.8 × 10-08), and higher Braak amyloid beta score (p = 4.4 × 10-04). Together these results from high-throughput and hypothesis-free approaches converge on a genetic link between Parkinson's disease, CSF amyloid beta1-42, and APOE.

Examination of the Effect of Rare Variants in TREM2, ABI3, and PLCG2 in LOAD Through Multiple Phenotypes.

BACKGROUND: Rare variants in PLCG2 (p.P522R), ABI3 (p.S209F), and TREM2 (p.R47H, p.R62H) have been associated with late onset Alzheimer's disease (LOAD) risk in Caucasians. After the initial report, several studies have found positive results in cohorts of different ethnic background and with different phenotype. OBJECTIVE: In this study, we aim to evaluate the association of rare coding variants in PLCG2, ABI3, and TREM2 with LOAD risk and their effect at different time points of the disease. METHODS: We used a European American cohort to assess the association of the variants prior onset (using CSF Aß42, tau, and pTau levels, and amyloid imaging as endophenotypes) and after onset (measured as rate of memory decline). RESULTS: We confirm the association with LOAD risk of TREM2 p.R47H, p.R62H and ABI3 p.S209F variants, and the protective effect of PLCG2 p.P522R. In addition, ABI3 and TREM2 gene-sets showed significant association with LOAD risk. TREM2 p.R47H and PLCG2 p.P522R variants were also statistically associated with increase of amyloid imaging and AD progression, respectively. We did not observe any association of ABI3 p.S209F with any of the other AD endophenotypes. CONCLUSION: The results of this study highlight the importance of including biomarkers and alternative phenotypes to better understand the role of novel candidate genes with the disease.

Correction to: Overlapping genetic architecture between Parkinson disease and melanoma.

The original version of this article unfortunately contained a mistake. Supplementary Tables 3 and 4 are not available with the rest of the supplementary material available online.

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.

Overlapping genetic architecture between Parkinson disease and melanoma.

Epidemiologic studies have reported inconsistent results regarding an association between Parkinson disease (PD) and cutaneous melanoma (melanoma). Identifying shared genetic architecture between these diseases can support epidemiologic findings and identify common risk genes and biological pathways. Here, we apply polygenic, linkage disequilibrium-informed methods to the largest available case-control, genome-wide association study summary statistic data for melanoma and PD. We identify positive and significant genetic correlation (correlation: 0.17, 95% CI 0.10-0.24; P = 4.09 × 10-06) between melanoma and PD. We further demonstrate melanoma and PD-inferred gene expression to overlap across tissues (correlation: 0.14, 95% CI 0.06 to 0.22; P = 7.87 × 10-04) and highlight seven genes including PIEZO1, TRAPPC2L, and SOX6 as potential mediators of the genetic correlation between melanoma and PD. These findings demonstrate specific, shared genetic architecture between PD and melanoma that manifests at the level of gene expression.

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

Parietal cortex RNA-sequencing (RNA-seq) data were generated from individuals with and without Alzheimer disease (AD; ncontrol = 13; nAD = 83) from the Knight Alzheimer Disease Research Center (Knight ADRC). Using this and an independent (Mount Sinai Brain Bank (MSBB)) AD RNA-seq dataset, cortical circular RNA (circRNA) expression was quantified in the context of AD. Significant associations were identified between circRNA expression and AD diagnosis, clinical dementia severity and neuropathological severity. It was demonstrated that most circRNA-AD associations are independent of changes in cognate linear messenger RNA expression or estimated brain cell-type proportions. Evidence was provided for circRNA expression changes occurring early in presymptomatic AD and in autosomal dominant AD. It was also observed that AD-associated circRNAs co-expressed with known AD genes. Finally, potential microRNA-binding sites were identified in AD-associated circRNAs for miRNAs 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.

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.

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.

Genome-wide association study for variants that modulate relationships between cerebrospinal fluid amyloid-beta 42, tau, and p-tau levels.

BACKGROUND: A relationship quantitative trait locus exists when the correlation between multiple traits varies by genotype for that locus. Relationship quantitative trait loci (rQTL) are often involved in gene-by-gene (G×G) interactions or gene-by-environmental interactions, making them a powerful tool for detecting G×G. METHODS: We performed genome-wide association studies to identify rQTL between tau and Aß42 and ptau and Aß42 with over 3000 individuals using age, gender, series, APOE ε2, APOE ε4, and two principal components for population structure as covariates. Each significant rQTL was separately screened for interactions with other loci for each trait in the rQTL model. Parametric bootstrapping was used to assess significance. RESULTS: We found four significant tau/Aß42 rQTL from three unique locations and six ptau/Aß42 rQTL from five unique locations. G×G screens with these rQTL produced four significant G×G interactions (one Aß42, two ptau, and one tau) with four rQTL where each second locus was from a unique location. On follow-up, rs1036819 and rs74025622 were associated with Alzheimer's disease (AD) case/control status; rs15205 and rs79099429 were associated with rate of decline. CONCLUSIONS: The two most significant rQTL (rs8027714 and rs1036819) for ptau/Aß42 are on different chromosomes and both are strong hits for pelvic organ prolapse. While diseases of the nervous system can cause pelvic organ prolapse, it is unlikely related to the ptau/Aß42 relationship but may suggest that these two loci share a pathway. In addition to a ptau/Aß42 rQTL and association with AD case/control status, rs1036819 is a strong rQTL for case/control status/Aß42 and for tau/Aß42. It resides in the ZFAT gene, which is related to autoimmune thyroid disease. For tau, rs9817620 interacts with the tau/Aß42 rQTL rs74025622. It is in the CHL1 gene, which is a neural cell adhesion molecule and may be involved in signal transduction pathways. CHL1 is related to BACE1, which is a ß-secretase enzyme that initiates production of the ß-amyloid peptide involved in AD and is a primary drug target. Overall, there are numerous loci that affect the relationship between these important AD endophenotypes and some are due to interactions with other loci. Some affect the risk of AD and/or rate of progression.

CYP2A6 metabolism in the development of smoking behaviors in young adults.

Cytochrome P450 2A6 (CYP2A6) encodes the enzyme responsible for the majority of nicotine metabolism. Previous studies support that slow metabolizers smoke fewer cigarettes once nicotine dependent but provide conflicting results on the role of CYP2A6 in the development of dependence. By focusing on the critical period of young adulthood, this study examines the relationship of CYP2A6 variation and smoking milestones. A total of 1209 European American young adults enrolled in the Collaborative Study on the Genetics of Alcoholism were genotyped for CYP2A6 variants to calculate a previously well-validated metric that estimates nicotine metabolism. This metric was not associated with the transition from never smoking to smoking initiation nor with the transition from initiation to daily smoking (P > 0.4). But among young adults who had become daily smokers (n = 506), decreased metabolism was associated with increased risk of nicotine dependence (P = 0.03) (defined as Fagerström Test for Nicotine Dependence score ≥4). This finding was replicated in the Collaborative Genetic Study of Nicotine Dependence with 335 young adult daily smokers (P = 0.02). Secondary meta-analysis indicated that slow metabolizers had a 53 percent increased odds (OR = 1.53, 95 percent CI 1.11-2.11, P = 0.009) of developing nicotine dependence compared with normal metabolizers. Furthermore, secondary analyses examining four-level response of time to first cigarette after waking (>60, 31-60, 6-30, ≤5 minutes) demonstrated a robust effect of the metabolism metric in Collaborative Study on the Genetics of Alcoholism (P = 0.03) and Collaborative Genetic Study of Nicotine Dependence (P = 0.004), illustrating the important role of this measure of dependence. These findings highlight the complex role of CYP2A6 variation across different developmental stages of smoking behaviors.

Analysis of neurodegenerative Mendelian genes in clinically diagnosed Alzheimer Disease.

Alzheimer disease (AD), Frontotemporal lobar degeneration (FTD), Amyotrophic lateral sclerosis (ALS) and Parkinson disease (PD) have a certain degree of clinical, pathological and molecular overlap. Previous studies indicate that causative mutations in AD and FTD/ALS genes can be found in clinical familial AD. We examined the presence of causative and low frequency coding variants in the AD, FTD, ALS and PD Mendelian genes, in over 450 families with clinical history of AD and over 11,710 sporadic cases and cognitive normal participants from North America. Known pathogenic mutations were found in 1.05% of the sporadic cases, in 0.69% of the cognitively normal participants and in 4.22% of the families. A trend towards enrichment, albeit non-significant, was observed for most AD, FTD and PD genes. Only PSEN1 and PINK1 showed consistent association with AD cases when we used ExAC as the control population. These results suggest that current study designs may contain heterogeneity and contamination of the control population, and that current statistical methods for the discovery of novel genes with real pathogenic variants in complex late onset diseases may be inadequate or underpowered to identify genes carrying pathogenic mutations.

Cis-regulatory variants affect CHRNA5 mRNA expression in populations of African and European ancestry.

Variants within the gene cluster encoding α3, α5, and ß4 nicotinic receptor subunits are major risk factors for substance dependence. The strongest impact on risk is associated with variation in the CHRNA5 gene, where at least two mechanisms are at work: amino acid variation and altered mRNA expression levels. The risk allele of the non-synonymous variant (rs16969968; D398N) primarily occurs on the haplotype containing the low mRNA expression allele. In populations of European ancestry, there are approximately 50 highly correlated variants in the CHRNA5-CHRNA3-CHRNB4 gene cluster and the adjacent PSMA4 gene region that are associated with CHRNA5 mRNA levels. It is not clear which of these variants contribute to the changes in CHRNA5 transcript level. Because populations of African ancestry have reduced linkage disequilibrium among variants spanning this gene cluster, eQTL mapping in subjects of African ancestry could potentially aid in defining the functional variants that affect CHRNA5 mRNA levels. We performed quantitative allele specific gene expression using frontal cortices derived from 49 subjects of African ancestry and 111 subjects of European ancestry. This method measures allele-specific transcript levels in the same individual, which eliminates other biological variation that occurs when comparing expression levels between different samples. This analysis confirmed that substance dependence associated variants have a direct cis-regulatory effect on CHRNA5 transcript levels in human frontal cortices of African and European ancestry and identified 10 highly correlated variants, located in a 9 kb region, that are potential functional variants modifying CHRNA5 mRNA expression levels.