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 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 stimulant dependence

Stimulant dependence is heritable, but specific genetic factors underlying the trait have not been identified. A genome-wide association study for stimulant dependence was performed in a discovery cohort of African- (AA) and European-ancestry (EA) subjects ascertained for genetic studies of alcohol, opioid, and cocaine use disorders. The sample comprised individuals with DSM-IV stimulant dependence (393 EA cases, 5288 EA controls; 155 AA cases, 5603 AA controls). An independent cohort from the family-based Collaborative Study on the Genetics of Alcoholism (532 EA cases, 7635 EA controls; 53 AA cases, AA 3352 controls) was used for replication. One variant in SLC25A16 (rs2394476, p = 3.42 × 10-10, odds ratio [OR] = 3.70) was GWS in AAs. Four other loci showed suggestive evidence, including KCNA4 in AAs (rs11500237, p = 2.99 × 10-7, OR = 2.31) which encodes one of the potassium voltage-gated channel protein that has been linked to several other substance use disorders, and CPVL in the combined population groups (rs1176440, p = 3.05 × 10-7, OR = 1.35), whose expression was previously shown to be upregulated in the prefrontal cortex from users of cocaine, cannabis, and phencyclidine. Analysis of the top GWAS signals revealed a significant enrichment with nicotinic acetylcholine receptor genes (adjusted p = 0.04) and significant pleiotropy between stimulant dependence and alcohol dependence in EAs (padj = 3.6 × 10-3), an anxiety disorder in EAs (padj = 2.1 × 10-4), and ADHD in both AAs (padj = 3.0 × 10-33) and EAs (padj = 6.7 × 10-35). Our results implicate novel genes and pathways as having roles in the etiology of stimulant dependence

Genome-wide association study of opioid cessation

The United States is experiencing an epidemic of opioid use disorder (OUD) and overdose-related deaths. However, the genetic basis for the ability to discontinue opioid use has not been investigated. We performed a genome-wide association study (GWAS) of opioid cessation (defined as abstinence from illicit opioids for \u3e1 year or \u3c6 months before the interview date) in 1130 African American (AA) and 2919 European ancestry (EA) participants recruited for genetic studies of substance use disorders and who met lifetime Diagnostic and Statistical Manual of Mental Disorders, 5th Edition (DSM-5) criteria for OUD. Association tests performed separately within each ethnic group were combined by meta-analysis with results obtained from the Comorbidity and Trauma Study. Although there were no genome-wide significant associations, we found suggestive associations with nine independent loci, including three which are biologically relevant: rs4740988 i

Genome-wide association study identifies loci associated with liability to alcohol and drug dependence that is associated with variability in reward-related ventral striatum activity in African- and European-Americans.

Genetic influences on alcohol and drug dependence partially overlap, however, specific loci underlying this overlap remain unclear. We conducted a genome-wide association study (GWAS) of a phenotype representing alcohol or illicit drug dependence (ANYDEP) among 7291 European-Americans (EA; 2927 cases) and 3132 African-Americans (AA: 1315 cases) participating in the family-based Collaborative Study on the Genetics of Alcoholism. ANYDEP was heritable (h 2 in EA = 0.60, AA = 0.37). The AA GWAS identified three regions with genome-wide significant (GWS; P < 5E-08) single nucleotide polymorphisms (SNPs) on chromosomes 3 (rs34066662, rs58801820) and 13 (rs75168521, rs78886294), and an insertion-deletion on chromosome 5 (chr5:141988181). No polymorphisms reached GWS in the EA. One GWS region (chromosome 1: rs1890881) emerged from a trans-ancestral meta-analysis (EA + AA) of ANYDEP, and was attributable to alcohol dependence in both samples. Four genes (AA: CRKL, DZIP3, SBK3; EA: P2RX6) and four sets of genes were significantly enriched within biological pathways for hemostasis and signal transduction. GWS signals did not replicate in two independent samples but there was weak evidence for association between rs1890881 and alcohol intake in the UK Biobank. Among 118 AA and 481 EA individuals from the Duke Neurogenetics Study, rs75168521 and rs1890881 genotypes were associated with variability in reward-related ventral striatum activation. This study identified novel loci for substance dependence and provides preliminary evidence that these variants are also associated with individual differences in neural reward reactivity. Gene discovery efforts in non-European samples with distinct patterns of substance use may lead to the identification of novel ancestry-specific genetic markers of risk

Genetic risk prediction and neurobiological understanding of alcoholism

We have used a translational Convergent Functional Genomics (CFG) approach to discover genes involved in alcoholism, by gene-level integration of genome-wide association study (GWAS) data from a German alcohol dependence cohort with other genetic and gene expression data, from human and animal model studies, similar to our previous work in bipolar disorder and schizophrenia. A panel of all the nominally significant P-value single-nucleotide length polymorphisms (SNPs) in the top candidate genes discovered by CFG (n = 135 genes, 713 SNPs) was used to generate a genetic risk prediction score (GRPS), which showed a trend towards significance (P = 0.053) in separating alcohol dependent individuals from controls in an independent German test cohort. We then validated and prioritized our top findings from this discovery work, and subsequently tested them in three independent cohorts, from two continents. In order to validate and prioritize the key genes that drive behavior without some of the pleiotropic environmental confounds present in humans, we used a stress-reactive animal model of alcoholism developed by our group, the D-box binding protein (DBP) knockout mouse, consistent with the surfeit of stress theory of addiction proposed by Koob and colleagues. A much smaller panel (n = 11 genes, 66 SNPs) of the top CFG-discovered genes for alcoholism, cross-validated and prioritized by this stress-reactive animal model showed better predictive ability in the independent German test cohort (P = 0.041). The top CFG scoring gene for alcoholism from the initial discovery step, synuclein alpha (SNCA) remained the top gene after the stress-reactive animal model cross-validation. We also tested this small panel of genes in two other independent test cohorts from the United States, one with alcohol dependence (P = 0.00012) and one with alcohol abuse (a less severe form of alcoholism; P = 0.0094). SNCA by itself was able to separate alcoholics from controls in the alcohol-dependent cohort (P = 0.000013) and the alcohol abuse cohort (P = 0.023). So did eight other genes from the panel of 11 genes taken individually, albeit to a lesser extent and/or less broadly across cohorts. SNCA, GRM3 and MBP survived strict Bonferroni correction for multiple comparisons. Taken together, these results suggest that our stress-reactive DBP animal model helped to validate and prioritize from the CFG-discovered genes some of the key behaviorally relevant genes for alcoholism. These genes fall into a series of biological pathways involved in signal transduction, transmission of nerve impulse (including myelination) and cocaine addiction. Overall, our work provides leads towards a better understanding of illness, diagnostics and therapeutics, including treatment with omega-3 fatty acids. We also examined the overlap between the top candidate genes for alcoholism from this work and the top candidate genes for bipolar disorder, schizophrenia, anxiety from previous CFG analyses conducted by us, as well as cross-tested genetic risk predictions. This revealed the significant genetic overlap with other major psychiatric disorder domains, providing a basis for comorbidity and dual diagnosis, and placing alcohol use in the broader context of modulating the mental landscape

ACSL6 Is Associated with the Number of Cigarettes Smoked and Its Expression Is Altered by Chronic Nicotine Exposure

Individuals with schizophrenia tend to be heavy smokers and are at high risk for tobacco dependence. However, the nature of the comorbidity is not entirely clear. We previously reported evidence for association of schizophrenia with SNPs and SNP haplotypes in a region of chromosome 5q containing the SPEC2, PDZ-GEF2 and ACSL6 genes. In this current study, analysis of the control subjects of the Molecular Genetics of Schizophrenia (MGS) sample showed similar pattern of association with number of cigarettes smoked per day (numCIG) for the same region. To further test if this locus is associated with tobacco smoking as measured by numCIG and FTND, we conducted replication and meta-analysis in 12 independent samples (n\u3e16,000) for two markers in ACSL6 reported in our previous schizophrenia study. In the meta-analysis of the replication samples, we found that rs667437 and rs477084 were significantly associated with numCIG (p = 0.00038 and 0.00136 respectively) but not with FTND scores. We then used in vitro and in vivo techniques to test if nicotine exposure influences the expression of ACSL6 in brain. Primary cortical culture studies showed that chronic (5-day) exposure to nicotine stimulated ACSL6 mRNA expression. Fourteen days of nicotine administration via osmotic mini pump also increased ACSL6 protein levels in the prefrontal cortex and hippocampus of mice. These increases were suppressed by injection of the nicotinic receptor antagonist mecamylamine, suggesting that elevated expression ofACSL6 requires nicotinic receptor activation. These findings suggest that variations in theACSL6 gene may contribute to the quantity of cigarettes smoked. The independent associations of this locus with schizophrenia and with numCIG in non-schizophrenic subjects suggest that this locus may be a common liability to both conditions

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

Association of the CHRNA5-A3-B4 Gene Cluster With Heaviness of Smoking: A Meta-Analysis

Introduction: Variation in the CHRNA5-A3-B4 gene cluster is a promising candidate region for smoking behavior and has been linked to multiple smoking-related phenotypes (e.g., nicotine dependence) and diseases (e.g., lung cancer). Two single nucleotide polymorphisms (SNPs), rs16969968 in CHRNA5 and rs1051730 in CHRNA3, have generated particular interest. Methods: We evaluated the published evidence for association between rs16969968 (k = 27 samples) and rs1051730 (k = 44 samples) SNPs with heaviness of smoking using meta-analytic techniques. We explored which SNP provided a stronger genetic signal and investigated study-level characteristics (i.e., ancestry, disease state) to establish whether the strength of association differed across populations. We additionally tested for small study bias and explored the impact of year of publication. Results and Conclusions: Meta-analysis indicated compelling evidence of an association between the rs1051730/rs16966968 variants and daily cigarette consumption (fixed effects: B = 0.91, 95% CI = 0.77, 1.06, p < .001; random effects: B = 1.01, 95% CI = 0.81, 1.22, p < .001), equivalent to a per-allele effect of approximately 1 cigarette/day. SNP rs1051730 was found to provide a stronger signal than rs16966968 in stratified analyses (pdiff = .028), although this difference was only qualitatively observed in the subset of samples that provided data on both SNPs. While the functional relevance of rs1051730 is unknown, it may be a strong tagging SNP for functional haplotypes in this region

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

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

Homozygous Deletion of Six Olfactory Receptor Genes in a Subset of Individuals with Beta-Thalassemia

Progress in the functional studies of human olfactory receptors has been largely hampered by the lack of a reliable experimental model system. Although transgenic approaches in mice could characterize the function of individual olfactory receptors, the presence of over 300 functional genes in the human genome becomes a daunting task. Thus, the characterization of individuals with a genetic susceptibility to altered olfaction coupled with the absence of particular olfactory receptor genes will allow phenotype/genotype correlations and vindicate the function of specific olfactory receptors with their cognate ligands. We characterized a 118 kb β-globin deletion and found that its 3′ end breakpoint extends to the neighboring olfactory receptor region downstream of the β-globin gene cluster. This deletion encompasses six contiguous olfactory receptor genes (OR51V1, OR52Z1, OR51A1P, OR52A1, OR52A5, and OR52A4) all of which are expressed in the brain. Topology analysis of the encoded proteins from these olfactory receptor genes revealed that OR52Z1, OR52A1, OR52A5, and OR52A4 are predicted to be functional receptors as they display integral characteristics of G-proteins coupled receptors. Individuals homozygous for the 118 kb β-globin deletion are afflicted with β-thalassemia due to a homozygous deletion of the β-globin gene and have no alleles for the above mentioned olfactory receptors genes. This is the first example of a homozygous deletion of olfactory receptor genes in human. Although altered olfaction remains to be ascertained in these individuals, such a study can be carried out in β-thalassemia patients from Malaysia, Indonesia and the Philippines where this mutation is common. Furthermore, OR52A1 contains a γ-globin enhancer, which was previously shown to confer continuous expression of the fetal γ-globin genes. Thus, the hypothesis that β-thalassemia individuals, who are homozygous for the 118 kb deletion, may also have an exacerbation of their anemia due to the deletion of two copies of the γ-globin enhancer element is worthy of consideration