Cell-type-specific expression quantitative trait loci associated with Alzheimer disease in blood and brain tissue

Because regulation of gene expression is heritable and context-dependent, we investigated AD-related gene expression patterns in cell types in blood and brain. Cis-expression quantitative trait locus (eQTL) mapping was performed genome-wide in blood from 5257 Framingham Heart Study (FHS) participants and in brain donated by 475 Religious Orders Study/Memory & Aging Project (ROSMAP) participants. The association of gene expression with genotypes for all cis SNPs within 1 Mb of genes was evaluated using linear regression models for unrelated subjects and linear-mixed models for related subjects. Cell-type-specific eQTL (ct-eQTL) models included an interaction term for the expression of "proxy" genes that discriminate particular cell type. Ct-eQTL analysis identified 11,649 and 2533 additional significant gene-SNP eQTL pairs in brain and blood, respectively, that were not detected in generic eQTL analysis. Of note, 386 unique target eGenes of significant eQTLs shared between blood and brain were enriched in apoptosis and Wnt signaling pathways. Five of these shared genes are established AD loci. The potential importance and relevance to AD of significant results in myeloid cell types is supported by the observation that a large portion of GWS ct-eQTLs map within 1 Mb of established AD loci and 58% (23/40) of the most significant eGenes in these eQTLs have previously been implicated in AD. This study identified cell-type-specific expression patterns for established and potentially novel AD genes, found additional evidence for the role of myeloid cells in AD risk, and discovered potential novel blood and brain AD biomarkers that highlight the importance of cell-type-specific analysis.

Set-Based Rare Variant Expression Quantitative Trait Loci in Blood and Brain from Alzheimer Disease Study Participants

Because studies of rare variant effects on gene expression have limited power, we investigated set-based methods to identify rare expression quantitative trait loci (eQTL) related to Alzheimer disease (AD). Gene-level and pathway-level cis rare-eQTL mapping was performed genome-wide using gene expression data derived from blood donated by 713 Alzheimer's Disease Neuroimaging Initiative participants and from brain tissues donated by 475 Religious Orders Study/Memory and Aging Project participants. The association of gene or pathway expression with a set of all cis potentially regulatory low-frequency and rare variants within 1 Mb of genes was evaluated using SKAT-O. A total of 65 genes expressed in the brain were significant targets for rare expression single nucleotide polymorphisms (eSNPs) among which 17% (11/65) included established AD genes HLA-DRB1 and HLA-DRB5. In the blood, 307 genes were significant targets for rare eSNPs. In the blood and the brain, GNMT, LDHC, RBPMS2, DUS2, and HP were targets for significant eSNPs. Pathway enrichment analysis revealed significant pathways in the brain (n = 9) and blood (n = 16). Pathways for apoptosis signaling, cholecystokinin receptor (CCKR) signaling, and inflammation mediated by chemokine and cytokine signaling were common to both tissues. Significant rare eQTLs in inflammation pathways included five genes in the blood (ALOX5AP, CXCR2, FPR2, GRB2, IFNAR1) that were previously linked to AD. This study identified several significant gene- and pathway-level rare eQTLs, which further confirmed the importance of the immune system and inflammation in AD and highlighted the advantages of using a set-based eQTL approach for evaluating the effect of low-frequency and rare variants on gene expression.

Genome-wide association study of phenotypes measuring progression from first cocaine or opioid use to dependence reveals novel risk genes

Aim: Substance use disorders (SUD) result in substantial morbidity and mortality worldwide. Opioids, and to a lesser extent cocaine, contribute to a large percentage of this health burden. Despite their high heritability, few genetic risk loci have been identified for either opioid or cocaine dependence (OD or CD, respectively). A genome-wide association study of OD and CD related phenotypes reflecting the time between first self-reported use of these substances and a first DSM-IV dependence diagnosis was conducted. Methods: Cox proportional hazards regression in a discovery sample of 6,188 African-Americans (AAs) and 6,835 European-Americans (EAs) participants in a genetic study of multiple substance dependence phenotypes were used to test for association between genetic variants and these outcomes. The top findings were tested for replication in two independent cohorts. Results: In the discovery sample, three independent regions containing variants associated with time to dependence at P < 5 x 10-8 were identified, one (rs61835088 = 1.03 x 10-8) for cocaine in the combined EA-AA meta-analysis in the gene FAM78B on chromosome 1, and two for opioids in the AA portion of the sample in intergenic regions of chromosomes 4 (rs4860439, P = 1.37 x 10-8) and 9 (rs7032521, P = 3.30 x 10-8). After meta-analysis with data from the replication cohorts, the signal at rs61835088 improved (HR = 0.87, P = 3.71 x 10-9 and an intergenic SNP on chromosome 21 (rs2825295, HR = 1.14, P = 2.57 x 10-8) that missed the significance threshold in the AA discovery sample became genome-wide significant (GWS) for CD. Conclusions: Although the two GWS variants are not in genes with obvious links to SUD biology and have modest effect sizes, they are statistically robust and show evidence for association in independent samples. These results may point to novel pathways contributing to disease progression and highlight the utility of related phenotypes to better understand the genetics of SUDs.

Genome-wide association study of phenotypes measuring progression from first cocaine or opioid use to dependence reveals novel risk genes

Aim

Substance use disorders (SUD) result in substantial morbidity and mortality worldwide. Opioids, and to a lesser extent cocaine, contribute to a large percentage of this health burden. Despite their high heritability, few genetic risk loci have been identified for either opioid or cocaine dependence (OD or CD, respectively). A genome-wide association study of OD and CD related phenotypes reflecting the time between first self-reported use of these substances and a first DSM-IV dependence diagnosis was conducted.

Methods

Cox proportional hazards regression in a discovery sample of 6,188 African-Americans (AAs) and 6,835 European-Americans (EAs) participants in a genetic study of multiple substance dependence phenotypes were used to test for association between genetic variants and these outcomes. The top findings were tested for replication in two independent cohorts.

Results

In the discovery sample, three independent regions containing variants associated with time to dependence at P < 5 x 10^-8 were identified, one (rs61835088 = 1.03 x 10^-8) for cocaine in the combined EA-AA meta-analysis in the gene FAM78B on chromosome 1, and two for opioids in the AA portion of the sample in intergenic regions of chromosomes 4 (rs4860439, P = 1.37 x 10^-8) and 9 (rs7032521, P = 3.30 x 10^-8). After meta-analysis with data from the replication cohorts, the signal at rs61835088 improved (HR = 0.87, P = 3.71 x 10^-9 and an intergenic SNP on chromosome 21 (rs2825295, HR = 1.14, P = 2.57 x 10^-8) that missed the significance threshold in the AA discovery sample became genome-wide significant (GWS) for CD.

Conclusions

Although the two GWS variants are not in genes with obvious links to SUD biology and have modest effect sizes, they are statistically robust and show evidence for association in independent samples. These results may point to novel pathways contributing to disease progression and highlight the utility of related phenotypes to better understand the genetics of SUDs.

Identifying factors associated with opioid cessation in a biracial sample using machine learning

Aim

Racial disparities in opioid use disorder (OUD) management exist, however, and there is limited research on factors that influence opioid cessation in different population groups.

Methods

We employed multiple machine learning prediction algorithms least absolute shrinkage and selection operator, random forest, deep neural network, and support vector machine to assess factors associated with ceasing opioid use in a sample of 1,192 African Americans (AAs) and 2,557 individuals of European ancestry (EAs) who met Diagnostic and Statistical Manual of Mental Disorders, 5th Edition criteria for OUD. Values for nearly 4,000 variables reflecting demographics, alcohol and other drug use, general health, non-drug use behaviors, and diagnoses for other psychiatric disorders, were obtained for each participant from the Semi-Structured Assessment for Drug Dependence and Alcoholism, a detailed semi-structured interview.

Results

Support vector machine models performed marginally better on average than other machine learning methods with maximum prediction accuracies of 75.4% in AAs and 79.4% in EAs. Subsequent stepwise regression considered the 83 most highly ranked variables across all methods and models and identified less recent cocaine use (AAs: odds ratio (OR) = 1.82, P = 9.19 × 10^-5; EAs: OR = 1.91, P = 3.30 × 10^-15), shorter duration of opioid use (AAs: OR = 0.55, P = 5.78 × 10^-6; EAs: OR = 0.69, P = 3.01 × 10^-7), and older age (AAs: OR = 2.44, P = 1.41 × 10^-12; EAs: OR = 2.00, P = 5.74 × 10^-9) as the strongest independent predictors of opioid cessation in both AAs and EAs. Attending self-help groups for OUD was also an independent predictor (P < 0.05) in both population groups, while less gambling severity (OR = 0.80, P = 3.32 × 10^-2) was specific to AAs and post-traumatic stress disorder recovery (OR = 1.93, P = 7.88 × 10^-5), recent antisocial behaviors (OR = 0.64, P = 2.69 × 10^-3), and atheism (OR = 1.45, P = 1.34 × 10^-2) were specific to EAs. Factors related to drug use comprised about half of the significant independent predictors in both AAs and EAs, with other predictors related to non-drug use behaviors, psychiatric disorders, overall health, and demographics.

Conclusions

These proof-of-concept findings provide avenues for hypothesis-driven analysis, and will lead to further research on strategies to improve OUD management in EAs and AAs.

Novel Alzheimer Disease Risk Loci and Pathways in African American Individuals Using the African Genome Resources Panel: A Meta-analysis

Importance

Compared with non-Hispanic White individuals, African American individuals from the same community are approximately twice as likely to develop Alzheimer disease. Despite this disparity, the largest Alzheimer disease genome-wide association studies to date have been conducted in non-Hispanic White individuals. In the largest association analyses of Alzheimer disease in African American individuals, ABCA7, TREM2, and an intergenic locus at 5q35 were previously implicated.

Objective

To identify additional risk loci in African American individuals by increasing the sample size and using the African Genome Resource panel.

Design, Setting, and Participants

This genome-wide association meta-analysis used case-control and family-based data sets from the Alzheimer Disease Genetics Consortium. There were multiple recruitment sites throughout the United States that included individuals with Alzheimer disease and controls of African American ancestry. Analysis began October 2018 and ended September 2019.

Main Outcomes and Measures

Diagnosis of Alzheimer disease.

Results

A total of 2784 individuals with Alzheimer disease (1944 female [69.8%]) and 5222 controls (3743 female [71.7%]) were analyzed (mean [SD] age at last evaluation, 74.2 [13.6] years). Associations with 4 novel common loci centered near the intracellular glycoprotein trafficking gene EDEM1 (3p26; P = 8.9 × 10-7), near the immune response gene ALCAM (3q13; P = 9.3 × 10-7), within GPC6 (13q31; P = 4.1 × 10-7), a gene critical for recruitment of glutamatergic receptors to the neuronal membrane, and within VRK3 (19q13.33; P = 3.5 × 10-7), a gene involved in glutamate neurotoxicity, were identified. In addition, several loci associated with rare variants, including a genome-wide significant intergenic locus near IGF1R at 15q26 (P = 1.7 × 10-9) and 6 additional loci with suggestive significance (P ≤ 5 × 10-7) such as API5 at 11p12 (P = 8.8 × 10-8) and RBFOX1 at 16p13 (P = 5.4 × 10-7) were identified. Gene expression data from brain tissue demonstrate association of ALCAM, ARAP1, GPC6, and RBFOX1 with brain β-amyloid load. Of 25 known loci associated with Alzheimer disease in non-Hispanic White individuals, only APOE, ABCA7, TREM2, BIN1, CD2AP, FERMT2, and WWOX were implicated at a nominal significance level or stronger in African American individuals. Pathway analyses strongly support the notion that immunity, lipid processing, and intracellular trafficking pathways underlying Alzheimer disease in African American individuals overlap with those observed in non-Hispanic White individuals. A new pathway emerging from these analyses is the kidney system, suggesting a novel mechanism for Alzheimer disease that needs further exploration.

Conclusions and Relevance

While the major pathways involved in Alzheimer disease etiology in African American individuals are similar to those in non-Hispanic White individuals, the disease-associated loci within these pathways differ.

A large-scale genome-wide association study meta-analysis of cannabis use disorder

BACKGROUND

Variation in liability to cannabis use disorder has a strong genetic component (estimated twin and family heritability about 50-70%) and is associated with negative outcomes, including increased risk of psychopathology. The aim of the study was to conduct a large genome-wide association study (GWAS) to identify novel genetic variants associated with cannabis use disorder.

METHODS

To conduct this GWAS meta-analysis of cannabis use disorder and identify associations with genetic loci, we used samples from the Psychiatric Genomics Consortium Substance Use Disorders working group, iPSYCH, and deCODE (20 916 case samples, 363 116 control samples in total), contrasting cannabis use disorder cases with controls. To examine the genetic overlap between cannabis use disorder and 22 traits of interest (chosen because of previously published phenotypic correlations [eg, psychiatric disorders] or hypothesised associations [eg, chronotype] with cannabis use disorder), we used linkage disequilibrium score regression to calculate genetic correlations.

FINDINGS

We identified two genome-wide significant loci: a novel chromosome 7 locus (FOXP2, lead single-nucleotide polymorphism [SNP] rs7783012; odds ratio [OR] 1·11, 95% CI 1·07-1·15, p=1·84 × 10-9) and the previously identified chromosome 8 locus (near CHRNA2 and EPHX2, lead SNP rs4732724; OR 0·89, 95% CI 0·86-0·93, p=6·46 × 10-9). Cannabis use disorder and cannabis use were genetically correlated (rg 0·50, p=1·50 × 10-21), but they showed significantly different genetic correlations with 12 of the 22 traits we tested, suggesting at least partially different genetic underpinnings of cannabis use and cannabis use disorder. Cannabis use disorder was positively genetically correlated with other psychopathology, including ADHD, major depression, and schizophrenia.

INTERPRETATION

These findings support the theory that cannabis use disorder has shared genetic liability with other psychopathology, and there is a distinction between genetic liability to cannabis use and cannabis use disorder.

FUNDING

National Institute of Mental Health; National Institute on Alcohol Abuse and Alcoholism; National Institute on Drug Abuse; Center for Genomics and Personalized Medicine and the Centre for Integrative Sequencing; The European Commission, Horizon 2020; National Institute of Child Health and Human Development; Health Research Council of New Zealand; National Institute on Aging; Wellcome Trust Case Control Consortium; UK Research and Innovation Medical Research Council (UKRI MRC); The Brain & Behavior Research Foundation; National Institute on Deafness and Other Communication Disorders; Substance Abuse and Mental Health Services Administration (SAMHSA); National Institute of Biomedical Imaging and Bioengineering; National Health and Medical Research Council (NHMRC) Australia; Tobacco-Related Disease Research Program of the University of California; Families for Borderline Personality Disorder Research (Beth and Rob Elliott) 2018 NARSAD Young Investigator Grant; The National Child Health Research Foundation (Cure Kids); The Canterbury Medical Research Foundation; The New Zealand Lottery Grants Board; The University of Otago; The Carney Centre for Pharmacogenomics; The James Hume Bequest Fund; National Institutes of Health: Genes, Environment and Health Initiative; National Institutes of Health; National Cancer Institute; The William T Grant Foundation; Australian Research Council; The Virginia Tobacco Settlement Foundation; The VISN 1 and VISN 4 Mental Illness Research, Education, and Clinical Centers of the US Department of Veterans Affairs; The 5th Framework Programme (FP-5) GenomEUtwin Project; The Lundbeck Foundation; NIH-funded Shared Instrumentation Grant S10RR025141; Clinical Translational Sciences Award grants; National Institute of Neurological Disorders and Stroke; National Heart, Lung, and Blood Institute; National Institute of General Medical Sciences.

Expanding the genetic architecture of nicotine dependence and its shared genetics with multiple traits

Cigarette smoking is the leading cause of preventable morbidity and mortality. Genetic variation contributes to initiation, regular smoking, nicotine dependence, and cessation. We present a Fagerström Test for Nicotine Dependence (FTND)-based genome-wide association study in 58,000 European or African ancestry smokers. We observe five genome-wide significant loci, including previously unreported loci MAGI2/GNAI1 (rs2714700) and TENM2 (rs1862416), and extend loci reported for other smoking traits to nicotine dependence. Using the heaviness of smoking index from UK Biobank (N = 33,791), rs2714700 is consistently associated; rs1862416 is not associated, likely reflecting nicotine dependence features not captured by the heaviness of smoking index. Both variants influence nearby gene expression (rs2714700/MAGI2-AS3 in hippocampus; rs1862416/TENM2 in lung), and expression of genes spanning nicotine dependence-associated variants is enriched in cerebellum. Nicotine dependence (SNP-based heritability = 8.6%) is genetically correlated with 18 other smoking traits (rg = 0.40-1.09) and co-morbidities. Our results highlight nicotine dependence-specific loci, emphasizing the FTND as a composite phenotype that expands genetic knowledge of smoking.

Association of OPRM1 Functional Coding Variant With Opioid Use Disorder: A Genome-Wide Association Study

IMPORTANCE

With the current opioid crisis, it is important to improve understanding of the biological mechanisms of opioid use disorder (OUD).

OBJECTIVES

To detect genetic risk variants for OUD and determine genetic correlations and causal association with OUD and other traits.

DESIGN, SETTING, AND PARTICIPANTS

A genome-wide association study of electronic health record-defined OUD in the Million Veteran Program sample was conducted, comprising 8529 affected European American individuals and 71 200 opioid-exposed European American controls (defined by electronic health record trajectory analysis) and 4032 affected African American individuals and 26 029 opioid-exposed African American controls. Participants were enrolled from January 10, 2011, to May 21, 2018, with electronic health record data for OUD diagnosis from October 1, 1999, to February 7, 2018. Million Veteran Program results and additional OUD case-control genome-wide association study results from the Yale-Penn and Study of Addiction: Genetics and Environment samples were meta-analyzed (total numbers: European American individuals, 10 544 OUD cases and 72 163 opioid-exposed controls; African American individuals, 5212 cases and 26 876 controls). Data on Yale-Penn participants were collected from February 14, 1999, to April 1, 2017, and data on Study of Addiction: Genetics and Environment participants were collected from 1990 to 2007. The key result was replicated in 2 independent cohorts: proxy-phenotype buprenorphine treatment in the UK Biobank and newly genotyped Yale-Penn participants. Genetic correlations between OUD and other traits were tested, and mendelian randomization analysis was conducted to identify potential causal associations.

MAIN OUTCOMES AND MEASURES

Main outcomes were International Classification of Diseases, Ninth Revision-diagnosed OUD or International Statistical Classification of Diseases and Related Health Problems, Tenth Revision-diagnosed OUD (Million Veteran Program), and DSM-IV-defined opioid dependence (Yale-Penn and Study of Addiction: Genetics and Environment).

RESULTS

A total of 114 759 individuals (101 016 men [88%]; mean [SD] age, 60.1 [12.8] years) were included. In 82 707 European American individuals, a functional coding variant (rs1799971, encoding Asn40Asp) in OPRM1 (μ-opioid receptor gene, the main biological target for opioid drugs; OMIM 600018) reached genome-wide significance (G allele: β = -0.066 [SE = 0.012]; P = 1.51 × 10-8). The finding was replicated in 2 independent samples. Single-nucleotide polymorphism-based heritability of OUD was 11.3% (SE = 1.8%). Opioid use disorder was genetically correlated with 83 traits, including multiple substance use traits, psychiatric illnesses, cognitive performance, and others. Mendelian randomization analysis revealed the following associations with OUD: risk of tobacco smoking, depression, neuroticism, worry neuroticism subcluster, and cognitive performance. No genome-wide significant association was detected for African American individuals or in transpopulation meta-analysis.

CONCLUSIONS AND RELEVANCE

This genome-wide meta-analysis identified a significant association of OUD with an OPRM1 variant, which was replicated in 2 independent samples. Post-genome-wide association study analysis revealed associated pleiotropic characteristics. Recruitment of additional individuals with OUD for future studies-especially those of non-European ancestry-is a crucial next step in identifying additional significant risk loci.

Analysis of brain region-specific co-expression networks reveals clustering of established and novel genes associated with Alzheimer disease

BACKGROUND

Identifying and understanding the functional role of genetic risk factors for Alzheimer disease (AD) has been complicated by the variability of genetic influences across brain regions and confounding with age-related neurodegeneration.

METHODS

A gene co-expression network was constructed using data obtained from the Allen Brain Atlas for multiple brain regions (cerebral cortex, cerebellum, and brain stem) in six individuals. Gene network analyses were seeded with 52 reproducible (i.e., established) AD (RAD) genes. Genome-wide association study summary data were integrated with the gene co-expression results and phenotypic information (i.e., memory and aging-related outcomes) from gene knockout studies in Drosophila to generate rankings for other genes that may have a role in AD.

RESULTS

We found that co-expression of the RAD genes is strongest in the cortical regions where neurodegeneration due to AD is most severe. There was significant evidence for two novel AD-related genes including EPS8 (FDR p = 8.77 × 10-3) and HSPA2 (FDR p = 0.245).

CONCLUSIONS

Our findings indicate that AD-related risk factors are potentially associated with brain region-specific effects on gene expression that can be detected using a gene network approach.

Genetic variants and functional pathways associated with resilience to Alzheimer's disease

Approximately 30% of older adults exhibit the neuropathological features of Alzheimer's disease without signs of cognitive impairment. Yet, little is known about the genetic factors that allow these potentially resilient individuals to remain cognitively unimpaired in the face of substantial neuropathology. We performed a large, genome-wide association study (GWAS) of two previously validated metrics of cognitive resilience quantified using a latent variable modelling approach and representing better-than-predicted cognitive performance for a given level of neuropathology. Data were harmonized across 5108 participants from a clinical trial of Alzheimer's disease and three longitudinal cohort studies of cognitive ageing. All analyses were run across all participants and repeated restricting the sample to individuals with unimpaired cognition to identify variants at the earliest stages of disease. As expected, all resilience metrics were genetically correlated with cognitive performance and education attainment traits (P-values < 2.5 × 10-20), and we observed novel correlations with neuropsychiatric conditions (P-values < 7.9 × 10-4). Notably, neither resilience metric was genetically correlated with clinical Alzheimer's disease (P-values > 0.42) nor associated with APOE (P-values > 0.13). In single variant analyses, we observed a genome-wide significant locus among participants with unimpaired cognition on chromosome 18 upstream of ATP8B1 (index single nucleotide polymorphism rs2571244, minor allele frequency = 0.08, P = 2.3 × 10-8). The top variant at this locus (rs2571244) was significantly associated with methylation in prefrontal cortex tissue at multiple CpG sites, including one just upstream of ATPB81 (cg19596477; P = 2 × 10-13). Overall, this comprehensive genetic analysis of resilience implicates a putative role of vascular risk, metabolism, and mental health in protection from the cognitive consequences of neuropathology, while also providing evidence for a novel resilience gene along the bile acid metabolism pathway. Furthermore, the genetic architecture of resilience appears to be distinct from that of clinical Alzheimer's disease, suggesting that a shift in focus to molecular contributors to resilience may identify novel pathways for therapeutic targets.