Genomewide linkage scan for cocaine dependence and related traits: significant linkages for a cocaine-related trait and cocaine-induced paranoia

Application of polygenic scores to a deeply phenotyped sample enriched for substance use disorders reveals extensive pleiotropy with psychiatric and somatic traits

Co-occurring psychiatric, medical, and substance use disorders (SUDs) are common, but the complex pathways leading to such comorbidities are poorly understood. A greater understanding of genetic influences on this phenomenon could inform precision medicine efforts. We used the Yale-Penn dataset, a cross-sectional sample enriched for individuals with SUDs, to examine pleiotropic effects of genetic liability for psychiatric and somatic traits. Participants completed an in-depth interview that provides information on demographics, environment, medical illnesses, and psychiatric and SUDs. Polygenic scores (PGS) for psychiatric disorders and somatic traits were calculated in European-ancestry (EUR; n = 5691) participants and, when discovery datasets were available, for African-ancestry (AFR; n = 4918) participants. Phenome-wide association studies (PheWAS) were then conducted. In AFR participants, the only PGS with significant associations was bipolar disorder (BD), all of which were with substance use phenotypes. In EUR participants, PGS for major depressive disorder (MDD), generalized anxiety disorder (GAD), post-traumatic stress disorder (PTSD), schizophrenia (SCZ), body mass index (BMI), coronary artery disease (CAD), and type 2 diabetes (T2D) all showed significant associations, the majority of which were with phenotypes in the substance use categories. For instance, PGSMDD was associated with over 200 phenotypes, 15 of which were depression-related (e.g., depression criterion count), 55 of which were other psychiatric phenotypes, and 126 of which were substance use phenotypes; and PGSBMI was associated with 138 phenotypes, 105 of which were substance related. Genetic liability for psychiatric and somatic traits is associated with numerous phenotypes across multiple categories, indicative of the broad genetic liability of these traits.

Application of polygenic scores to a deeply phenotyped sample enriched for substance use disorders reveals extensive pleiotropy with psychiatric and medical traits

Co-occurring psychiatric, medical, and substance use disorders (SUDs) are common, but the complex pathways leading to such comorbidities are poorly understood. A greater understanding of genetic influences on this phenomenon could inform precision medicine efforts. We used the Yale-Penn dataset, a cross-sectional sample enriched for individuals with SUDs, to examine pleiotropic effects of genetic liability for psychiatric and medical traits. Participants completed an in-depth interview that provides information on demographics, environment, medical illnesses, and psychiatric and SUDs. Polygenic scores (PGS) for psychiatric disorders and medical traits were calculated in European-ancestry (EUR; n=5,691) participants and, when discovery datasets were available, for African-ancestry (AFR; n=4,918) participants. Phenome-wide association studies (PheWAS) were then conducted. In AFR participants, the only PGS with significant associations was bipolar disorder (BD), all of which were with substance use phenotypes. In EUR participants, PGS for major depressive disorder (MDD), generalized anxiety disorder (GAD), post-traumatic stress disorder (PTSD), schizophrenia (SCZ), body mass index (BMI), coronary artery disease (CAD), and type 2 diabetes (T2D) all showed significant associations, the majority of which were with phenotypes in the substance use categories. For instance, PGS MDD was associated with over 200 phenotypes, 15 of which were depression-related (e.g., depression criterion count), 55 of which were other psychiatric phenotypes, and 126 of which were substance use phenotypes; and PGS BMI was associated with 138 phenotypes, 105 of which were substance related. Genetic liability for psychiatric and medical traits is associated with numerous phenotypes across multiple categories, indicative of the broad genetic liability of these traits.

Familial associations in adolescent substance use disorder: a population-based cohort study

BACKGROUND AND AIMS

Family history of substance use disorder (SUD) affects a child's risk of the disorder through both genetic and shared environmental factors. We aimed to estimate the association between parental or older sibling SUD history with the risk of adolescent SUD diagnosis.

DESIGN, SETTING AND PARTICIPANTS

We conducted a population-based cohort study using administrative health-care databases in the Province of Manitoba, Canada, which has a universal and publicly funded health-care system. We included all children born from 1984 to 2000 who have linkages to both parents and were followed until age 18 years. We used generalized estimating equation models to produce unadjusted and adjusted relative risk (RR) estimates of adolescent SUD risk. The study cohort included 134 389 children and 31 307 full sibling pairs; 51.3% were male and 35.4% first-born.

MEASUREMENTS

The exposure was SUD diagnosis in a mother or father in either hospitalization or outpatient physician visit records before the children's age of 13 years. The secondary exposure was an adolescent SUD diagnosis in an older full sibling. The outcome was SUD diagnosis during adolescence (13 and 18 years of age) identified in either hospitalization or physician visit records. Children demographics and characteristics associated with SUD diagnosis were included in the models.

FINDINGS

Of the 134 389 children, 9.5% had a mother with a history of SUD, 11.3% had a father and 1.3% had an older sibling with a history of SUD diagnosis; 2566 (1.9%) had an adolescent SUD diagnosis. An increased risk of adolescent SUD was observed with SUD history in mothers [adjusted RR (aRR) = 2.50; 95% confidence interval (CI) = 2.26, 2.79], fathers (aRR = 2.15; 95% CI = 1.95, 2.37), both parents (aRR = 3.74; 95% CI = 3.24, 4.31) and older sibling (aRR = 3.85; 95% CI = 2.53, 5.87).

CONCLUSIONS

A family history of substance use disorder in parents or older siblings appears to be associated with increased SUD risk in adolescents.

Reliability and validity of DSM-IV and DSM-5 methamphetamine use disorder diagnoses using the Chinese Version of the Semi-Structured Assessment for Drug Dependence and Alcoholism (SSADDA)

BACKGROUND

The Semi-structured Assessment for Drug Dependence and Alcoholism (SSADDA) was developed to assess substance-use disorders and other psychiatric traits. We translated the SSADDA into Chinese and evaluated its inter-rater reliability and concurrent validity in diagnosing DSM-IV methamphetamine (MA) dependence and DSM-5 MA-use disorder (MUD).

METHODS

The sample comprised 231 participants who were interviewed using the Chinese SSADDA and the Mini-International Neuropsychiatric Interview (Chinese MINI) for concurrent validation. Of the 231 participants, 191 were interviewed by two different interviewers two weeks apart. We evaluated the inter-rater reliability and concurrent validity of the diagnoses using percent agreement and Cohen's kappa coefficient (κ). Cohen's linear weighted kappa was used to assess the reliability of DSM-5 MUD severity.

RESULTS

It showed good inter-rater reliability and no significant differences among the DSM-5 MUD (κ = 0.71), DSM-IV MA abuse or dependence (κ = 0.72), and the DSM-IV diagnoses of MA dependence (κ = 0.66) and abuse (κ = 0.68) tested separately. The weighted kappa was 0.67 across the three DSM-5 MUD severity levels. The reliability of each individual diagnostic criterion for DSM-5 MUD ranged from fair to excellent (κ = 0.41-0.80), except for "repeated attempts to quit/control use" (κ = 0.38). The concurrent validity based on MINI-derived diagnoses ranged from good to excellent (κ = 0.65-0.78).

CONCLUSIONS

This study shows that the Chinese version of SSADDA has good reliability and validity among Chinese MA users.

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.

An update on the role of common genetic variation underlying substance use disorders

PURPOSE OF THE REVIEW

Sample size increases have resulted in novel and replicable loci for substance use disorders (SUDs). We summarize some of the latest insights into SUD genetics and discuss some next steps in addiction genetics.

RECENT FINDINGS

Genome-wide association studies have substantiated the role of previously known variants (e.g., rs1229984 in ADH1B for alcohol) and identified several novel loci for alcohol, tobacco, cannabis, opioid and cocaine use disorders. SUDs are genetically correlated with psychiatric outcomes, while liability to substance use is inconsistently associated with these outcomes and more closely associated with lifestyle factors. Specific variant associations appear to differ somewhat across populations, although similar genes and systems are implicated.

SUMMARY

The next decade of human genetic studies of addiction should focus on expanding to non-European populations, consider pleiotropy across SUD and with other psychiatric disorders, and leverage human and cross-species functional data to elucidate the biological mechanisms underlying SUDs.

Prospects for finding the mechanisms of sex differences in addiction with human and model organism genetic analysis

Despite substantial evidence for sex differences in addiction epidemiology, addiction-relevant behaviors and associated neurobiological phenomena, the mechanisms and implications of these differences remain unknown. Genetic analysis in model organism is a potentially powerful and effective means of discovering the mechanisms that underlie sex differences in addiction. Human genetic studies are beginning to show precise risk variants that influence the mechanisms of addiction but typically lack sufficient power or neurobiological mechanistic access, particularly for the discovery of the mechanisms that underlie sex differences. Our thesis in this review is that genetic variation in model organisms are a promising approach that can complement these investigations to show the biological mechanisms that underlie sex differences in addiction.

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 >1 year or <6 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 in PTPRD (p_{AA + EA} = 2.24 × 10⁻⁶), rs36098404 in MYOM2 (p_EA = 2.24 × 10⁻⁶), and rs592026 in SNAP25-AS1 (p_EA = 6.53 × 10⁻⁶). Significant pathways identified in persons of European ancestry (EA) are related to vitamin D metabolism (p = 3.79 × 10⁻²) and fibroblast growth factor (FGF) signaling (p = 2.39 × 10⁻²). UK Biobank traits including smoking and drinking cessation and chronic back pain were significantly associated with opioid cessation using GWAS-derived polygenic risk scores. These results provide evidence for genetic influences on opioid cessation, suggest genetic overlap with other relevant traits, and may indicate potential novel therapeutic targets for OUD.

A genome-wide association study of cocaine use disorder accounting for phenotypic heterogeneity and gene–environment interaction

BACKGROUND

Phenotypic heterogeneity and complicated gene–environment interplay in etiology are among the primary factors that hinder the identification of genetic variants associated with cocaine use disorder.

METHODS

To detect novel genetic variants associated with cocaine use disorder, we derived disease traits with reduced phenotypic heterogeneity using cluster analysis of a study sample (n = 9965). We then used these traits in genome-wide association tests, performed separately for 2070 African Americans and 1570 European Americans, using a new mixed model that accounted for the moderating effects of 5 childhood environmental factors. We used an independent sample (918 African Americans, 1382 European Americans) for replication.

RESULTS

The cluster analysis yielded 5 cocaine use disorder subtypes, of which subtypes 4 (n = 3258) and 5 (n = 1916) comprised heavy cocaine users, had high heritability estimates (h2 = 0.66 and 0.64, respectively) and were used in association tests. Seven of the 13 identified genetic loci in the discovery phase were available in the replication sample. In African Americans, rs114492924 (discovery p = 1.23 × E−8), a single nucleotide polymorphism in LINC01411, was replicated in the replication sample (p = 3.63 × E−3). In a meta-analysis that combined the discovery and replication results, 3 loci in African Americans were significant genome-wide: rs10188036 in TRAK2 (p = 2.95 × E−8), del-1:15511771 in TMEM51 (p = 9.11 × E−10) and rs149843442 near LPHN2 (p = 3.50 × E−8).

LIMITATIONS

Lack of data prevented us from replicating 6 of the 13 identified loci.

CONCLUSION

Our results demonstrate the importance of considering phenotypic heterogeneity and gene–environment interplay in detecting genetic variations that contribute to cocaine use disorder, because new genetic loci have been identified using our novel analytic method.

Statistical Methods and Software for Substance Use and Dependence Genetic Research

BACKGROUND

Substantial substance use disorders and related health conditions emerged dur-ing the mid-20th century and continue to represent a remarkable 21st century global burden of disease. This burden is largely driven by the substance-dependence process, which is a complex process and is influenced by both genetic and environmental factors. During the past few decades, a great deal of pro-gress has been made in identifying genetic variants associated with Substance Use and Dependence (SUD) through linkage, candidate gene association, genome-wide association and sequencing studies.

METHODS

Various statistical methods and software have been employed in different types of SUD ge-netic studies, facilitating the identification of new SUD-related variants.

CONCLUSION

In this article, we review statistical methods and software that are currently available for SUD genetic studies, and discuss their strengths and limitations.

Genetic associations with suicide attempt severity and genetic overlap with major depression

In 2015, ~800,000 people died by suicide worldwide. For every death by suicide there are as many as 25 suicide attempts, which can result in serious injury even when not fatal. Despite this large impact on morbidity and mortality, the genetic influences on suicide attempt are poorly understood. We performed a genome-wide association study (GWAS) of severity of suicide attempts to investigate genetic influences. A discovery GWAS was performed in Yale-Penn sample cohorts of European Americans (EAs, n = 2,439) and African Americans (AAs, n = 3,881). We found one genome-wide significant (GWS) signal in EAs near the gene LDHB (rs1677091, p = 1.07 × 10^-8) and three GWS associations in AAs: ARNTL2 on chromosome 12 (rs683813, p = 2.07 × 10^-8), FAH on chromosome 15 (rs72740082, p = 2.36 × 10^-8), and on chromosome 18 (rs11876255, p = 4.61 × 10^-8) in the Yale-Penn discovery sample. We conducted a limited replication analysis in the completely independent Army-STARRS cohorts. rs1677091 replicated in Latinos (LAT, p = 6.52 × 10^-3). A variant in LD with FAH rs72740082 (rs72740088; r² = 0.68) was replicated in AAs (STARRS AA p = 5.23 × 10^-3; AA meta, 1.51 × 10^-9). When combined for a trans-population meta-analysis, the final sample size included n = 20,153 individuals. Finally, we found significant genetic overlap with major depressive disorder (MDD) using polygenic risk scores from a large GWAS (r² = 0.007, p = 6.42 × 10^-5). To our knowledge, this is the first GWAS of suicide attempt severity. We identified GWS associations near genes involved in anaerobic energy production (LDHB), circadian clock regulation (ARNTL2), and catabolism of tyrosine (FAH). These findings provide evidence of genetic risk factors for suicide attempt severity, providing new information regarding the molecular mechanisms involved.

Sex-specific linkage scans in opioid dependence

Sex influences risk for opioid dependence (OD). We hypothesized that sex might interact with genetic loci that influence the risk for OD. Therefore we performed an analysis to identify sex-specific genomic susceptibility regions for OD using linkage. Over 6,000 single nucleotide polymorphism (SNP) markers were genotyped for 1,758 African- and European-American (AA and EA) individuals from 739 families, ascertained via affected sib-pairs with OD and/or cocaine dependence. Autosomewide non-parametric linkage scans, stratified by sex and population, were performed. We identified one significant linkage region, segregating with OD in EA men, at 71.1 cM on chromosome 4 (LOD = 3.29; point-wise P = 0.00005; empirical autosome-wide P = 0.042), which significantly differed from the linkage signal at the same location in EA women (empirical P = 0.002). Three suggestive linkage signals were identified at 181.3 cM on chromosome 7 (LOD = 2.18), 104 cM on chromosome 11 (LOD = 1.85), and 60.9 cM on chromosome 16 (LOD = 1.93) in EA women. In AA men, four suggestive linkage signals were detected at 201.1 cM on chromosome 3 (LOD = 2.32), 152.9 cM on chromosome 6 (LOD = 1.86), 16.8 cM on chromosome 7 (LOD = 1.95), and 36.1 cM on chromosome 17 (LOD = 1.99). The significant region, mapping to 4q12-4q13.1, harbors several OD candidate genes with interconnected functionality, including VEGFR, CLOCK, PDCL2, NMU, NRSF, and IGFBP7. In conclusion, these results provide an evidence for the existence of sex-specific and population-specific differences in OD. Furthermore, these results provide positional information that will facilitate the use of targeted next-generation sequencing to search for genes that contribute to sex-specific differences in OD. © 2016 Wiley Periodicals, Inc.

Biochemical Diagnosis in Substance and Non-substance Addiction

An optimal biochemical marker for addiction would be some easily traced molecules in body specimens, which indicates indulgent addictive behaviors, or susceptibility to certain addictive stimuli. In this chapter, we discussed existing literature about possible biomarkers, and classified them into three categories: origin forms and metabolites of substances, markers from biochemical responses to certain addiction, and genetic and epigenetic biomarkers suggesting susceptibility to addiction. In every category, we examined studies concerning certain type of addiction one by one, with focuses mainly on opiates, psychostimulants, and pathological gambling. Several promising molecules were highlighted, including those of neurotrophic factors, inflammatory factors, and indicators of vascular injury, and genetic and epigenetic biomarkers such as serum miRNAs. DNA methylation signatures and signal nucleotide polymorphism of candidate gene underlying the addiction.

A Review of Genome-Wide Association Studies of Stimulant and Opioid Use Disorders

Substance use disorders (SUD) are a major contributor to disability and disease burden worldwide. Risk for developing SUDs is influenced by variation in the genome. Identifying the genetic variants that influence SUD risk may help us to understand the biological mechanisms for the disorders and improve treatments. Genome-wide association studies (GWAS) have been successful in identifying many regions of the genome associated with common human disorders. Here, findings from recent GWAS of SUDs that involve illicit substances will be reviewed. Several GWAS have been reported, including studies on opioid and stimulant use disorder (cocaine and methamphetamine). Several of these GWAS report associations that are biologically interesting and statistically robust. Replication of the associations in independent samples and functional studies to understand the basis for the statistical associations will be important next steps.

Genomewide Association Study for Maximum Number of Alcoholic Drinks in European Americans and African Americans

BACKGROUND

We conducted a genomewide association study (GWAS) for maximum number of alcoholic drinks consumed in a 24-hour period ("MaxDrinks"), in 2 independent samples comprised of over 9,500 subjects, following up on our GWAS for alcohol dependence (AD) in European Americans (EAs) and African Americans (AAs).

METHODS

The samples included our GWAS samples (Yale-UPenn) recruited for studies of the genetics of drug or AD, and a publicly available sample: the Study of Addiction: Genetics and Environment (SAGE). Genomewide association analysis was performed for ~890,000 single nucleotide polymorphisms (SNPs) using linear association random effects models. EAs and AAs were separately analyzed.

RESULTS

The results confirmed significant associations of the well-known functional loci at ADH1B with MaxDrinks in EAs (rs1229984 Arg48His p = 5.96 × 10(-15) ) and AAs (rs2066702 Arg370Cys, p = 2.50 × 10(-10) ). The region of significant association on chromosome 4 was extended to LOC100507053 in AAs but not EAs. We also identified potentially novel significant common SNPs for MaxDrinks in EAs in the Yale-UPenn sample: rs1799876 at SERPINC1 on chromosome 1 (4.00 × 10(-8) ) and rs2309169 close to ANKRD36 on chromosome 2 (p = 5.58 × 10(-9) ). After adjusting for the peak SNP rs1229984 on ADH1B, rs1799876 was nearly significant (p = 1.99 × 10(-7) ) and rs2309169 remained highly significant (2.12 × 10(-9) ).

CONCLUSIONS

The results provide further support that ADH1B modulates alcohol consumption. Future replications of potential novel loci are warranted. This is the largest MaxDrinks GWAS to date, the first in AAs.

A review of pharmacogenetic studies of substance-related disorders

BACKGROUND

Substance-related disorders (SRDs) are a major cause of morbidity and mortality worldwide. Family, twin, and adoption studies have demonstrated the substantial heritability of SRDs. To determine the impact of genetic variation on risk for SRD and the response to treatment, researchers have conducted a number of secondary data analyses and quasi-experimental studies that target one or more candidate gene variants.

METHODS

This review examines studies in which candidate polymorphisms were examined as mediator variables to identify pharmacogenetic effects on subjective responses to drug administration or cues or outcomes of medication trials for SRDs. Efforts to use a meta-analytic approach to quantify these effects are premature because the number of available studies using similar methods and outcomes is limited, so the present review is qualitative.

RESULTS

Findings from these studies provide preliminary evidence of clinically relevant pharmacogenetic effects. However, independent replication of these findings has been sparse.

CONCLUSIONS

Although this growing body of literature has produced conflicting results, improved statistical controls may help to clarify the findings. Additionally, the use of empirically derived sub-phenotypes (i.e., which serve to differentiate distinct groups of affected individuals) may also help to identify genetic mediators of pharmacologic response in relation to SRDs. The identification of genetic mediators can inform clinical care both by identifying risk factors for SRDs and predicting adverse events and therapeutic outcomes associated with specific pharmacotherapies.

Genome-wide association study of copy number variations (CNVs) with opioid dependence

Single-nucleotide polymorphisms that have been associated with opioid dependence (OD) altogether account for only a small proportion of the known heritability. Most of the genetic risk factors are unknown. Some of the 'missing heritability' might be explained by copy number variations (CNVs) in the human genome. We used Illumina HumanOmni1 arrays to genotype 5152 African-American and European-American OD cases and screened controls and implemented combined CNV calling methods. After quality control measures were applied, a genome-wide association study (GWAS) of CNVs with OD was performed. For common CNVs, two deletions and one duplication were significantly associated with OD genome-wide (eg, P=2 × 10(-8) and OR (95% CI)=0.64 (0.54-0.74) for a chromosome 18q12.3 deletion). Several rare or unique CNVs showed suggestive or marginal significance with large effect sizes. This study is the first GWAS of OD using CNVs. Some identified CNVs harbor genes newly identified here to be of biological importance in addiction, whereas others affect genes previously known to contribute to substance dependence risk. Our findings augment our specific knowledge of the importance of genomic variation in addictive disorders, and provide an addiction CNV pool for further research. These findings require replication.

Examining the role of common genetic variants on alcohol, tobacco, cannabis and illicit drug dependence: genetics of vulnerability to drug dependence

BACKGROUND AND AIMS

Twin and family studies suggest that genetic influences are shared across substances of abuse. However, despite evidence of heritability, genome-wide association and candidate gene studies have indicated numerous markers of limited effects, suggesting that much of the heritability remains missing. We estimated (1) the aggregate effect of common single nucleotide polymorphisms (SNPs) on multiple indicators of comorbid drug problems that are typically employed across community and population-based samples, and (2) the genetic covariance across these measures.

PARTICIPANTS

A total of 2596 unrelated subjects from the Study of Addiction: Genetics and Environment provided information on alcohol, tobacco, cocaine, cannabis and other illicit substance dependence. Phenotypic measures included: (1) a factor score based on DSM-IV drug dependence diagnoses (DD), (2) a factor score based on problem use (PU; i.e. 1+ DSM-IV symptoms) and (3) dependence vulnerability (DV; a ratio of DSM-IV symptoms to the number of substances used).

FINDINGS

Univariate and bivariate genome-wide complex trait analyses of this selected sample indicated that common SNPs explained 25-36% of the variance across measures, with DD and DV having the largest effects [h(2) SNP (standard error) = 0.36 (0.13) and 0.33 (0.13), respectively; PU = 0.25 (0.13)]. Genetic effects were shared across the three phenotypic measures of comorbid drug problems [rDD-PU = 0.92 (0.08), rDD-DV = 0.97 (0.08) and rPU-DV = 0.96 (0.07)].

CONCLUSION

At least 20% of the variance in the generalized vulnerability to substance dependence is attributable to common single nucleotide polymorphisms. The additive effect of common single nucleotide polymorphisms is shared across important indicators of comorbid drug problems.

Genome-wide association study of opioid dependence: multiple associations mapped to calcium and potassium pathways

BACKGROUND

We report a genome-wide association study (GWAS) of two populations, African-American and European-American (AA, EA) for opioid dependence (OD) in three sets of subjects, to identify pathways, genes, and alleles important in OD risk.

METHODS

The design employed three phases (on the basis of separate sample collections). Phase 1 included our discovery GWAS dataset consisting of 5697 subjects (58% AA) diagnosed with opioid and/or other substance dependence and control subjects. Subjects were genotyped with the Illumina OmniQuad microarray, yielding 890,000 single nucleotide polymorphisms (SNPs) suitable for analysis. Additional genotypes were imputed with the 1000 Genomes reference panel. Top-ranked findings were further evaluated in Phase 2 by incorporating information from the publicly available Study of Addiction: Genetics and Environment dataset, with GWAS data from 4063 subjects (32% AA). In Phase 3, the most significant SNPs from Phase 2 were genotyped in 2549 independent subjects (32% AA). Analyses were performed with case-control and ordinal trait designs.

RESULTS

Most significant results emerged from the AA subgroup. Genome-wide-significant associations (p < 5.0 × 10(-8)) were observed with SNPs from multiple loci-KCNG2*rs62103177 was most significant after combining results from datasets in every phase of the study. The most compelling results were obtained with genes involved in potassium signaling pathways (e.g., KCNC1 and KCNG2). Pathway analysis also implicated genes involved in calcium signaling and long-term potentiation.

CONCLUSIONS

This is the first study to identify risk variants for OD with GWAS. Our results strongly implicate risk pathways and provide insights into novel therapeutic and prevention strategies and might biologically bridge OD and other non-substance dependence psychiatric traits where similar pathways have been implicated.

Genome-wide association study of cocaine dependence and related traits: FAM53B identified as a risk gene

We report a genome-wide association study (GWAS) for cocaine dependence (CD) in three sets of African- and European-American subjects (AAs and EAs, respectively) to identify pathways, genes and alleles important in CD risk. The discovery GWAS data set (n=5697 subjects) was genotyped using the Illumina OmniQuad microarray (8 90 000 analyzed single-nucleotide polymorphisms (SNPs)). Additional genotypes were imputed based on the 1000 Genomes reference panel. Top-ranked findings were evaluated by incorporating information from publicly available GWAS data from 4063 subjects. Then, the most significant GWAS SNPs were genotyped in 2549 independent subjects. We observed one genome-wide-significant (GWS) result: rs2629540 at the FAM53B ('family with sequence similarity 53, member B') locus. This was supported in both AAs and EAs; P-value (meta-analysis of all samples)=4.28 × 10(-8). The gene maps to the same chromosomal region as the maximum peak we observed in a previous linkage study. NCOR2 (nuclear receptor corepressor 2) SNP rs150954431 was associated with P=1.19 × 10(-9) in the EA discovery sample. SNP rs2456778, which maps to CDK1 ('cyclin-dependent kinase 1'), was associated with cocaine-induced paranoia in AAs in the discovery sample only (P=4.68 × 10(-8)). This is the first study to identify risk variants for CD using GWAS. Our results implicate novel risk loci and provide insights into potential therapeutic and prevention strategies.

Association of gamma-aminobutyric acid A receptor α2 gene (GABRA2) with alcohol use disorder

Gamma-aminobutyric acid (GABA) is a major inhibitory neurotransmitter in mammalian brain. GABA receptor are involved in a number of complex disorders, including substance abuse. No variants of the commonly studied GABA receptor genes that have been associated with substance dependence have been determined to be functional or pathogenic. To reconcile the conflicting associations with substance dependence traits, we performed a meta-analysis of variants in the GABAA receptor genes (GABRB2, GABRA6, GABRA1, and GABRG2 on chromosome 5q and GABRA2 on chromosome 4p12) using genotype data from 4739 cases of alcohol, opioid, or methamphetamine dependence and 4924 controls. Then, we combined the data from candidate gene association studies in the literature with two alcohol dependence (AD) samples, including 1691 cases and 1712 controls from the Study of Addiction: Genetics and Environment (SAGE), and 2644 cases and 494 controls from our own study. Using a Bonferroni-corrected threshold of 0.007, we found strong associations between GABRA2 and AD (P=9 × 10(-6) and odds ratio (OR) 95% confidence interval (CI)=1.27 (1.15, 1.4) for rs567926, P=4 × 10(-5) and OR=1.21 (1.1, 1.32) for rs279858), and between GABRG2 and both dependence on alcohol and dependence on heroin (P=0.0005 and OR=1.22 (1.09, 1.37) for rs211014). Significant association was also observed between GABRA6 rs3219151 and AD. The GABRA2 rs279858 association was observed in the SAGE data sets with a combined P of 9 × 10(-6) (OR=1.17 (1.09, 1.26)). When all of these data sets, including our samples, were meta-analyzed, associations of both GABRA2 single-nucleotide polymorphisms remained (for rs567926, P=7 × 10(-5) (OR=1.18 (1.09, 1.29)) in all the studies, and P=8 × 10(-6) (OR=1.25 (1.13, 1.38)) in subjects of European ancestry and for rs279858, P=5 × 10(-6) (OR=1.18 (1.1, 1.26)) in subjects of European ancestry. Findings from this extensive meta-analysis of five GABAA receptor genes and substance abuse support their involvement (with the best evidence for GABRA2) in the pathogenesis of AD. Further replications with larger samples are warranted.

Comparing the utility of homogeneous subtypes of cocaine use and related behaviors with DSM-IV cocaine dependence as traits for genetic association analysis

Because DSM-IV cocaine dependence (CD) is heterogeneous, it is not an optimal phenotype to identify genetic variation contributing to risk for cocaine use and related behaviors (CRBs). We used a cluster analytic method to differentiate homogeneous, highly heritable subtypes of CRBs and to compare their utility with that of the DSM-IV CD as traits for genetic association analysis. Clinical features of CRBs and co-occurring disorders were obtained via a poly-diagnostic interview administered to 9,965 participants in genetic studies of substance dependence. A subsample of subjects (N = 3,443) were genotyped for 1,350 single nucleotide polymorphisms (SNPs) selected from 130 candidate genes related to addiction. Cluster analysis of clinical features of the sample yielded five subgroups, two of which were characterized by heavy cocaine use and high heritability: a heavy cocaine use, infrequent intravenous injection group and an early-onset, heavy cocaine use, high comorbidity group. The utility of these traits was compared with the CD diagnosis through association testing of 2,320 affected subjects and 480 cocaine-exposed controls. Analyses examined both single SNP (main) and SNP-SNP interaction (epistatic) effects, separately for African-Americans and European-Americans. The two derived subtypes showed more significant P values for 6 of 8 main effects and 7 of 8 epistatic effects. Variants in the CLOCK gene were significantly associated with the heavy cocaine use, infrequent intravenous injection group, but not with the DSM-IV diagnosis of CD. These results support the utility of subtypes based on CRBs to detect risk variants for cocaine addiction.

ALDH2 is associated to alcohol dependence and is the major genetic determinant of "daily maximum drinks" in a GWAS study of an isolated rural Chinese sample

Alcohol dependence (AD) is a moderately heritable phenotype with a small number of known risk genes mapped via linkage or candidate gene studies. We considered 313 males from among 595 members of documented, extended pedigrees in which AD segregates collected in Northern Hunan Province, China. A joint analysis of both males and females could not be performed as the difference in alcohol consumption variance was too large. Genome-wide association analyses were performed for approximately 300,000 single nucleotide polymorphisms (SNPs). Significant associations found in the ALDH2 region for AD (minimum P = 4.73 × 10(-8)) and two AD-related phenotypes: flushing response (minimum P = 4.75 × 10(-26)) and maximum drinks in a 24-hr period (minimum P = 1.54 × 10(-16)). Association of previous candidate SNP, rs10774610 in CCDC63, was confirmed but resulted from linkage disequilibrium with ALDH2. ALDH2 is strongly associated with flushing response, AD, and maximum drinks in males, with nonsynonymous SNP rs671 explaining 29.2%, 7.9%, and 22.9% of phenotypic variation, respectively, in this sample. When rs671 was considered as a candidate SNP in females, it explained 23.6% of the variation in flushing response, but alcohol consumption rates were too low among females-despite familial enrichment for AD-for an adequate test of association for either AD or maximum drinks. These results support a mediating effect of aldehyde dehydrogenase deficiency on alcohol consumption in males and a secondary, culturally mediated limitation on alcohol consumption by females that should be appropriately modeled in future studies of alcohol consumption in populations where this may be a factor.

Butyrylcholinesterase genetic variants: association with cocaine dependence and related phenotypes

Objective

The search for genetic vulnerability factors in cocaine dependence has focused on the role that neuroplasticity plays in addiction. However, like many other drugs, the ability of an individual to metabolize cocaine can also influence susceptibility to dependence. Butyrylcholinesterase (BChE) metabolizes cocaine, and genetic variants of the BChE gene (BCHE) alter its catalytic activity. Therefore, we hypothesize that cocaine users with polymorphisms in BCHE can show diverse addictive behaviors due to differences in effective plasma concentrations of cocaine. Those polymorphisms might also influence users to prefer one of the two main preparations (crack or powder cocaine), despite having equal access to both. The present work investigates polymorphisms in BCHE and if those genetic variants constitute risk factors for cocaine dependence and for crack cocaine use.

Methods

A total of 1,436 individuals (698 cocaine-dependent patients and 738 controls) were genotyped for three single nucleotide polymorphisms (SNPs) in BCHE: rs1803274, rs4263329, and rs4680662.

Results

For rs4263329, a nominal difference was found between cases and controls. For rs1803274 (the functional SNP), a statistically significant difference was found between patients who used crack cocaine exclusively and those who used only powder cocaine (P = 0.027; OR = 4.36; 95% CI = 1.18–16.04). Allele frequencies and genotypes related to other markers did not differ between cases and controls or between the two cocaine subgroups.

Conclusions

Our findings suggest that the AA genotype of rs1803274 is a risk factor for crack cocaine use, which is more addictive than powder cocaine use. Further studies are needed in order to confirm this preliminary result and clarify the role of BCHE and its variants in cocaine dependence.

Critical evaluation of transcription factor Atf2 as a candidate modulator of alcohol preference in mouse and human populations

In prior work, congenic strains carrying the DBA/2Igb (D2) region of chromosome 2 (Chr2) for alcohol preference were bred onto a C57BL/6Ibg (B6) background and as predicted were found to reduce voluntary consumption. Subsequently, interval-specific congenic recombinant strains (ISCRS) were generated and also tested. These ISCRS strains reduced the quantitative trait loci (QTL) interval to a comparatively small 3.4 Mb region. Here, we have exploited an integrative approach using both murine and human populations to critically evaluate candidate genes within this region. First, we used bioinformatics tools to search for genes relevant to alcohol preference within the QTL region. Second, we searched for single nucleotide polymorphisms (SNPs) within exons of every gene in this region. Third, we conducted follow-up microarray analyses to identify differentially expressed genes between the B6 and ISCRS strains in mice from each group. Fourth, we analyzed correlations between the expression level of candidate genes and phenotypes of alcohol preference in a large family of BXD recombinant inbred strains derived from B6 and D2. Finally, we evaluated SNP segregation in both BXD mouse strains and in humans who were heavy alcohol drinkers or non-drinkers. Among several potential candidate genes in this region, we identified activating transcription factor 2 (Atf2) as the most plausible gene that would influence alcohol preference. However, the candidacy of Atf2 was only weakly supported when we used a genetic network approach and by focused reanalysis of genome-wide association study data from European-American and African-American populations. Thus, we cannot conclude that Atf2 plays a role in the regulation of the QTL of mouse Chr2.

Introduction to deep sequencing and its application to drug addiction research with a focus on rare variants

Through linkage analysis, candidate gene approach, and genome-wide association studies (GWAS), many genetic susceptibility factors for substance dependence have been discovered such as the alcohol dehydrogenase gene (ALDH2) for alcohol dependence (AD) and nicotinic acetylcholine receptor (nAChR) subunit variants on chromosomes 8 and 15 for nicotine dependence (ND). However, these confirmed genetic factors contribute only a small portion of the heritability responsible for each addiction. Among many potential factors, rare variants in those identified and unidentified susceptibility genes are supposed to contribute greatly to the missing heritability. Several studies focusing on rare variants have been conducted by taking advantage of next-generation sequencing technologies, which revealed that some rare variants of nAChR subunits are associated with ND in both genetic and functional studies. However, these studies investigated variants for only a small number of genes and need to be expanded to broad regions/genes in a larger population. This review presents an update on recently developed methods for rare-variant identification and association analysis and on studies focused on rare-variant discovery and function related to addictions.

GABRA2 genotype, impulsivity, and body mass

BACKGROUND

The goal of this study was to test a hypothesis associating impulsivity with an elevated body mass index (BMI).

METHODS

To this end, we examined associations of BMI with putative genetic, neurophysiological, psychiatric, and psychological indicators of impulsivity in 78 women and 74 men formerly dependent on alcohol or drugs. A second analysis was designed to test the replicability of the genetic findings in an independent sample of 109 women and 111 men with a similar history of substance dependence.

RESULTS

The results of the first analysis showed that BMI was positively correlated with Total and Nonplanning Scale Scores on the Barratt Impulsiveness Scale and the number of childhood symptoms of Attention-Deficit/Hyperactivity Disorder in women. It was also positively correlated, in women, with a GABRA2 variant previously implicated as a risk factor for substance dependence and an objective electroencephalographic feature previously associated with GABRA2 and relapse risk. The second analysis confirmed that the correlation between BMI and the substance-dependence-associated GABRA2 genotype was reliable and sex-specific.

CONCLUSIONS

We conclude that an elevated BMI is associated with genetic, neurophysiological, psychiatric, and psychological indicators of impulsivity. The sex difference may be explained by greater opportunities to eat and overeat, a preference for higher calorie foods, a longer duration of alcohol/drug abstinence, or previous pregnancies in women.

Ordered subset linkage analysis based on admixture proportion identifies new linkage evidence for alcohol dependence in African-Americans

Genetic heterogeneity could reduce the power of linkage analysis to detect risk loci for complex traits such as alcohol dependence (AD). Previously, we performed a genomewide linkage analysis for AD in African-Americans (AAs) (Biol Psychiatry 65:111-115, 2009). The power of that linkage analysis could have been reduced by the presence of genetic heterogeneity owing to differences in admixture among AA families. We hypothesized that by examining a study sample whose genetic ancestry was more homogeneous, we could increase the power to detect linkage. To test this hypothesis, we performed ordered subset linkage analysis in 384 AA families using admixture proportion as a covariate to identify a more homogeneous subset of families and determine whether there is increased evidence for linkage with AD. Statistically significant increases in lod scores in subsets relative to the overall sample were identified on chromosomes 4 (P = 0.0001), 12 (P = 0.021), 15 (P = 0.026) and 22 (P = 0.0069). In a subset of 44 families with African ancestry proportions ranging from 0.858 to 0.996, we observed a genomewide significant linkage at 180 cM on chromosome 4 (lod = 4.24, pointwise P < 0.00001, empirical genomewide P = 0.008). A promising candidate gene located there, GLRA3, which encodes a subunit of the glycine neurotransmitter receptor. Our results demonstrate that admixture proportion can be used as a covariate to reduce genetic heterogeneity and enhance the detection of linkage for AD in an admixed population such as AAs. This approach could be applied to any linkage analysis for complex traits conducted in an admixed population.

Linkage analysis followed by association show NRG1 associated with cannabis dependence in African Americans

BACKGROUND

A genetic contribution to cannabis dependence (CaD) has been established but susceptibility genes for CaD remain largely unknown.

METHODS

We employed a multistage design to identify genetic variants underlying CaD. We first performed a genome-wide linkage scan for CaD in 384 African American (AA) and 354 European American families ascertained for genetic studies of cocaine and opioid dependence. We then conducted association analysis under the linkage peak, first using data from a genome-wide association study from the Study of Addiction: Genetics and Environment, followed by replication studies of prioritized single nucleotide polymorphisms (SNPs) in independent samples.

RESULTS

We identified the strongest linkage evidence with CaD (logarithm of odds = 2.9) on chromosome 8p21.1 in AAs. In the association analysis of the Study of Addiction: Genetics and Environment sample under the linkage peak, we identified one SNP (rs17664708) associated with CaD in both AAs (odds ratio [OR] = 2.93, p = .0022) and European Americans (OR = 1.38, p = .02). This SNP, located at NRG1, a susceptibility gene for schizophrenia, was prioritized for further study. We replicated the association of rs17664708 with CaD in an independent AAs sample (OR = 2.81, p = .0068). The joint analysis of the two AA samples demonstrated highly significant association between rs17664708 and CaD with adjustment for either global (p = .00044) or local ancestry (p = .00075).

CONCLUSIONS

Our study shows that NRG1 is probably a susceptibility gene for CaD, based on convergent evidence of linkage and replicated associations in two independent AA samples.

Improved methods to identify stable, highly heritable subtypes of opioid use and related behaviors

Although there is evidence that opioid dependence (OD) is heritable, efforts to identify genes contributing to risk for the disorder have been hampered by its complex etiology and variable clinical manifestations. Decomposition of a complex set of opioid users into homogeneous subgroups could enhance genetic analysis. We applied a series of data mining techniques, including multiple correspondence analysis, variable selection and cluster analysis, to 69 opioid-related measures from 5390 subjects aggregated from family-based and case-control genetic studies to identify homogeneous subtypes and estimate their heritability. Novel aspects of this work include our use of (1) heritability estimates of specific clinical features of OD to enhance the heritability of the subtypes and (2) a k-medoids clustering method in combination with hierarchical clustering to yield replicable clusters that are less sensitive to noise than previous methods. We identified five homogeneous groups, including two large groups comprised of 762 and 1353 heavy opioid users, with estimated heritability of 0.69 and 0.76, respectively. These methods represent a promising approach to the identification of highly heritable subtypes in complex, heterogeneous disorders.

Polymorphisms in AKR1C4 and HSD3B2 and differences in serum DHEAS and progesterone are associated with paranoid ideation during mania or hypomania in bipolar disorder

Paranoia is commonly a mood-incongruent psychotic symptom of mania which may be related to dopamine dysregulation. Progesterone and its metabolite allopregnanolone (ALLO) have been found in animals to antagonize the effects of dopamine. We therefore examined serum progesterone, its endogenous antagonist DHEAS and polymorphisms of the genes coding for certain steroidogenetic enzymes (AKR1C4, HSD3B2, and SRD5A1) in 64 males and 96 females with bipolar 1 or 2 disorder with or without paranoid ideation during mood elevation. Euthymic morning serum progesterone, DHEAS and cortisol concentrations were measured in males and in premenopausal women who were in follicular phase and not taking oral contraceptives. In women only, SNPs in AKR1C4 reduced the likelihood of having exhibited paranoid ideation by circa 60%. The haplotype of all 4 SNPs in the AKR1C4 gene reduced the risk of exhibiting paranoia by 80% (OR 0.19, 95% CI 0.06-0.61, p=0.05). A history of paranoid ideation was not, however, related to progesterone or DHEAS concentration. Serum DHEAS and progesterone concentrations were lower in men who had shown paranoid ideation during mania/hypomania compared with those who had not (F=7.30, p=0.006) however this was not coupled to polymorphisms in the selected genes. The ancestral G in rs4659174 in HSD3B2 was in men associated with a lower risk of paranoid ideation (likelihood ratio χ(2) 3.97, p=0.046, OR 0.31 (95% CI 0.10-0.96)) but did not correlate with hormone concentrations. Hence, gene variants in the steroidogenetic pathway and steroids concentration differences may be involved in the susceptibility to paranoia during mood elevation.

The genetics of substance dependence

A large segment of the population suffers from addiction to alcohol, smoking, or illicit drugs. Not only do substance abuse and addiction pose a threat to health, but the consequences of addiction also impose a social and economic burden on families, communities, and nations. Genome-wide linkage and association studies have been used for addiction research with varying degrees of success. The most well-established genetic factors associated with alcohol dependence are in the genes encoding alcohol dehydrogenase (ADH), which oxidizes alcohol to acetaldehyde, and aldehyde dehydrogenase (ALDH2), which oxidizes acetaldehyde to acetate. Recently emerging genetic studies have linked variants in the genes encoding the α3, α5, and β4 nicotinic acetylcholine receptor subunits to smoking risk. However, the influence of these well-established genetic variants accounts for only a small portion of the heritability of alcohol and nicotine addiction, and it is likely that there are both common and rare risk variants yet to be identified. Newly developed DNA sequencing technologies could potentially advance the detection of rare variants with a larger impact on addiction risk.

Utility of genetically modified mice for understanding the neurobiology of substance use disorders

Advances in our ability to modify the mouse genome have enhanced our understanding of the genetic and neurobiological mechanisms contributing to addiction-related behaviors underlying substance use and abuse. These experimentally induced manipulations permit greater spatial and temporal specificity for modification of gene expression within specific cellular populations and during select developmental time periods. In this review, we consider the current mouse genetic model systems that have been employed to understand aspects of addiction and highlight significant conceptual advances achieved related to substance use and abuse. The mouse models reviewed herein include conventional knock-out and knock-in, conditional knockout, transgenic, inducible transgenic, mice suitable for optogenetic control of discrete neuronal populations, and phenotype-selected mice. By establishing a reciprocal investigatory relationship between genetic findings in humans and genomic manipulations in mice, a far better understanding of the discrete neuromechanisms underlying addiction can be achieved, which is likely to provide a strong foundation for developing and validating novel therapeutics for the treatment of substance abuse disorders.

Autosomal linkage scan for loci predisposing to comorbid dependence on multiple substances

Multiple substance dependence (MSD) trait comorbidity is common, and MSD patients are often severely affected clinically. While shared genetic risks have been documented, so far there has been no published report using the linkage scan approach to survey risk loci for MSD as a phenotype. A total of 1,758 individuals in 739 families [384 African American (AA) and 355 European American (EA) families] ascertained via affected sib-pairs with cocaine or opioid or alcohol dependence were genotyped using an array-based linkage panel of single-nucleotide polymorphism markers. Fuzzy clustering analysis was conducted on individuals with alcohol, cannabis, cocaine, opioid, and nicotine dependence for AAs and EAs separately, and linkage scans were conducted for the output membership coefficients using Merlin-regression. In EAs, we observed an autosome-wide significant linkage signal on chromosome 4 (peak lod = 3.31 at 68.3 cM; empirical autosome-wide P = 0.038), and a suggestive linkage signal on chromosome 21 (peak lod = 2.37 at 19.4 cM). In AAs, four suggestive linkage peaks were observed: two peaks on chromosome 10 (lod = 2.66 at 96.7 cM and lod = 3.02 at 147.6 cM] and the other two on chromosomes 3 (lod = 2.81 at 145.5 cM) and 9 (lod = 1.93 at 146.8 cM). Three particularly promising candidate genes, GABRA4, GABRB1, and CLOCK, are located within or very close to the autosome-wide significant linkage region for EAs on chromosome 4. This is the first linkage evidence supporting existence of genetic loci influencing risk for several comorbid disorders simultaneously in two major US populations.

Association of COL25A1 with comorbid antisocial personality disorder and substance dependence

BACKGROUND

Antisocial personality disorder (ASPD) is a psychiatric disorder characterized by a long-term pattern of manipulating, exploiting, or violating the rights of others.

METHODS

Subjects ascertained for genetic studies of substance dependence (SD) and diagnosed with ASPD and comorbid SD were included in a two-stage genetic association study. In the discovery stage, 627 single nucleotide polymorphisms (SNPs) located in 179 candidate genes for addiction were analyzed in a case-control cohort and family-based cohort. The significant findings were replicated in an independent case-control cohort.

RESULTS

One SNP, rs13134663, in the collagen XXV alpha 1 gene (COL25A1) was significantly associated with ASPD in both African Americans and European Americans (smallest p values were .0002 and .0004, respectively). There was also evidence of association with the same SNP in independent samples of African American and European American cases and control subjects (p = .035 and .033, respectively). Analysis of the combined set of case-control subjects yielded an allelic p value of 9 × 10(-6) with odds ratio (95% confidence interval) of 1.3 (1.16, 1.47) (smallest p = 1 × 10(-7); Bonferroni threshold p = .00012).

CONCLUSIONS

The COL25A1 gene, located at chromosome 4q25, encodes the collagen-like Alzheimer amyloid plaque component precursor, a type II transmembrane protein specifically expressed in neurons; it co-localizes with amyloid β in senile plaques in Alzheimer disease brains. This SNP maps to the transcription factor binding site and is conserved in 17 vertebrates, including mice and rats. Our findings suggest that COL25A1 may be associated with ASPD, especially in the context of SD.

5-HTTLPR as a potential moderator of the effects of adverse childhood experiences on risk of antisocial personality disorder

INTRODUCTION

Antisocial personality disorder (ASPD) frequently co-occurs with substance dependence (SD). A functional polymorphism (5-HTTLPR) in the serotonin transporter gene has been widely studied as a risk factor for a variety of psychopathologic conditions including aggressive/violent behavior. Childhood abuse is an important predictor of ASPD. We examined 5-HTTLPR genotype and adverse childhood events (ACEs) as risk factors for ASPD in a SD sample.

MATERIALS AND METHODS

Study participants [602 European-Americans (EAs) and 779 African-Americans (AAs)] were interviewed to obtain lifetime diagnoses of ASPD and SD and information on ACEs. Triallelic genotypes for 5-HTTLPR were obtained using standard methods. We used logistic generalized estimating equations regression to examine ACEs and 5-HTTLPR genotype and their interaction as predictors of ASPD, separately by population group.

RESULTS

There were 203 (14.7%) participants diagnosed with ASPD. The frequency of the low-activity 5-HTTLPR S' allele did not differ by ASPD diagnosis, and there was no overall 5-HTTLPR×ACE interaction. However, among European-Americans, male sex (odds ratio=3.36; P<0.001) and ACE history (odds ratio=1.47; P=0.002) were significant predictors of ASPD. Among AAs, there was a significant interaction of sex×5-HTTLPR genotype×ACEs (χ=13.92, P<0.001). Among AA men, each additional ACE significantly increased the odds of ASPD irrespective of genotype, whereas among AA women, the effect of ACEs on ASPD was significant only among S' homozygotes. However, these results are limited by the small sample size in each subgroup, (particularly AA women with S'S' genotype; N=7) and require replication.

CONCLUSIONS

Childhood maltreatment contributes to the risk of ASPD, an effect for which there is preliminary evidence of moderation by 5-HTTLPR genotype in AA women.

Childhood adversity increases risk for nicotine dependence and interacts with α5 nicotinic acetylcholine receptor genotype specifically in males

The relative importance of specific genetic and environmental factors in regulating nicotine dependence (ND) risk, including the effects on specific forms of childhood adversity on smoking risk, have been understudied. Genome-wide association studies and rodent models have demonstrated that the α5 nicotinic acetylcholine receptor gene (CHRNA5) is important in regulating nicotine intake. Childhood adversity increases the methylation level of the CHRNA5 promoter region in European Americans (EAs), an effect that was observed only in males (Zhang et al, submitted for publication). In view of this potential sex difference in the effects of early life experience on smoking, we investigated the presence of a sex-specific gene-by-environment effect of this marker on ND risk. A nonsynonymous SNP in CHRNA5 previously associated to ND and several related traits, rs16969968, was genotyped in 2206 EAs (1301 men and 905 women). The main and interactive effects of childhood adversity and rs16969968 genotype on diagnosis of ND and ND defined by dichotomized Fagerstrom test for ND (FTND) scores were explored. Men and women were analyzed separately to test for sex differences. Childhood adversity significantly increased ND risk in both sexes, and the effect in women was twice than that in men. Significant interactive effects of childhood adversity and rs16969968 genotype were observed in men (ND: OR=1.80, 95% CI=1.18-2.73, P=0.0044; FTND: OR=1.79, 95% CI=1.11-2.88, P=0.012). No interaction was found in women. This study provides evidence of a sex-specific gene × environment effect of CHRNA5 and childhood adversity on the risk for ND.

Variation in OPRM1 and risk of suicidal behavior in drug-dependent individuals

Completed suicide and nonfatal suicide-related outcomes (SROs), such as suicidal ideation and attempts, are heritable. A recent genetic association study in a sample of suicide victims reported a protective effect of the G allele of Asn40Asp (rs1799971) on risk for completed suicide. We examined the association of three OPRM1 single nucleotide polymorphisms (SNPs) (rs1799971, rs609148, and rs648893) with SRO in 426 European Americans, using GEE logistic regression analysis to examine the association of a lifetime history of SRO. There was no allelic association with the SRO phenotypes. A larger sample may be needed to identify risk variants that convey SRO risk. OPRM1 may not be important in the risk of SRO. 

A CRHR1 haplotype moderates the effect of adverse childhood experiences on lifetime risk of major depressive episode in African-American women

Adverse childhood experiences (ACEs) increase the risk for adult depression and substance dependence, possibly mediated by the corticotropin-releasing hormone type 1 receptor (CRHR1). In some studies, a three-SNP "T-A-T" haplotype in CRHR1, which encodes CRHR1, exerted a protective moderating effect on risk of depression in adults with ACEs. Other studies have shown a main or moderating effect of SNPs in CRHR1 on alcohol consumption. We tested the moderating effects of the three-SNP haplotype on lifetime risk of a major depressive episode (MDE) and alcohol dependence (AD) in 1,211 European-Americans (EAs) and 1,869 African-Americans (AAs), most of whom had a lifetime substance use disorder. There were no significant main or interaction effects of the TAT haplotype on AD. There was a significant interaction of ACE by TAT on risk of depression only in AA women (P = 0.005); each copy of the TAT haplotype reduced the odds of MDE by almost 40% (OR = 0.63). In AA women without an ACE and two TAT haplotypes, the risk of MDE was increased (OR = 1.51 for each copy). Our findings in relation to the TAT haplotype of CRHR1 extend those obtained in other populations to a largely substance-dependent one. The complex structure of CRHR1 may help to explain why some variants in the gene moderate the effects of an ACE only on depression risk while others moderate the effect of an ACE only on AD risk.

Linkage analyses of stimulant dependence, craving, and heavy use in American Indians

Amphetamine-type substances are the second most widely used illicit drugs in the United States. There is evidence to suggest that stimulant use (cocaine and methamphetamine) has a heritable component, yet the areas of the genome underlying these use disorders are yet to be identified. This study's aims were to map loci linked to stimulant dependence, heavy use, and craving in an American Indian community at high risk for substance dependence. DSM diagnosis of stimulant dependence, as well as indices of stimulant "craving," and "heavy use," were obtained using the Semi-Structured Assessment for the Genetics of Alcoholism (SSAGA). Genotypes were determined for a panel of 791 microsatellite polymorphisms in 381 members of multiplex families using SOLAR. Stimulant dependence, stimulant "craving," and "heavy stimulant use," were all found to be heritable. Analyses of multipoint variance component LOD scores, failed to yield evidence of linkage for stimulant dependence. For the stimulant "craving" phenotype, linkage analysis revealed a locus that had a LOD score of 3.02 on chromosome 15q25.3-26.1 near the nicotinic receptor gene cluster. A LOD score of 2.05 was found at this same site for "heavy stimulant use." Additional loci with LOD scores above 2.00 were found for stimulant "craving" on chromosomes 12p13.33-13.32 and 18q22.3. These results corroborate the importance of "craving" as an important phenotype that is associated with regions on chromosome 12, 15, and 18, that have been highlighted in prior segregation studies in this and other populations for substance dependence-related phenotypes.

A genomewide linkage scan of cocaine dependence and major depressive episode in two populations

Cocaine dependence (CD) and major depressive episode (MDE) frequently co-occur with poorer treatment outcome and higher relapse risk. Shared genetic risk was affirmed; to date, there have been no reports of genomewide linkage scans (GWLSs) surveying the susceptibility regions for comorbid CD and MDE (CD-MDE). We aimed to identify chromosomal regions and candidate genes susceptible to CD, MDE, and CD-MDE in African Americans (AAs) and European Americans (EAs). A total of 1896 individuals were recruited from 384 AA and 355 EA families, each with at least a sibling-pair with CD and/or opioid dependence. Array-based genotyping of about 6000 single-nucleotide polymorphisms was completed for all individuals. Parametric and non-parametric genomewide linkage analyses were performed. We found a genomewide-significant linkage peak on chromosome 7 at 183.4 cM for non-parametric analysis of CD-MDE in AAs (lod=3.8, genomewide empirical p=0.016; point-wise p=0.00001). A nearly genomewide significant linkage was identified for CD-MDE in EAs on chromosome 5 at 14.3 cM (logarithm of odds (lod)=2.95, genomewide empirical p=0.055; point-wise p=0.00012). Parametric analysis corroborated the findings in these two regions and improved the support for the peak on chromosome 5 so that it reached genomewide significance (heterogeneity lod=3.28, genomewide empirical p=0.046; point-wise p=0.00053). This is the first GWLS for CD-MDE. The genomewide significant linkage regions on chromosomes 5 and 7 harbor four particularly promising candidate genes: SRD5A1, UBE3C, PTPRN2, and VIPR2. Replication of the linkage findings in other populations is warranted, as is a focused analysis of the genes located in the linkage regions implicated here.

GABRG1 and GABRA2 variation associated with alcohol dependence in African Americans

BACKGROUND

GABRG1 and GABRA2, genes that encode the γ1 and α2 subunits, respectively, of the GABA-A receptor, are located in a cluster on chromosome 4p. Association of alcohol dependence (AD) with markers located at the 3' region of GABRA2 has been replicated in several studies, but recent studies suggested the possibility that the signal may be attributable to the adjacent gene, GABRG1, located 90 kb distant in the 3' direction. Owing to strong linkage disequilibrium (LD) in European Americans (EAs), the origin, or origins, of the association signal is very difficult to discern, but our previous population-based study suggested that decreased LD across the GABRG1-GABRA2 region in African Americans (AAs) may be useful for fine mapping and resolution of the association signal in that population.

METHODS

To examine these associations in greater detail, we genotyped 13 single nucleotide polymorphisms (SNPs) spanning GABRG1 and GABRA2 in 380 AAs with AD and in 253 AA controls.

RESULTS

Although there was no association between any individual SNP and AD, a highly significant difference was shown between AD subjects and controls in the frequency of a 3-SNP GABRA2 haplotype (global p = 0.00029). A similar level of significance was obtained in 6-SNP haplotypes that combined tagging SNPs from both genes (global p = 0.00994). High statistical significance was also shown with a 6-SNP haplotype (T-G-C-G-T-A), p = 0.0033. The T-G-C-G-T-A haplotype contains the most significant GABRA2 3-SNP haplotype (p = 0.00019), G-T-A.

CONCLUSIONS

These findings reflect the interrelationship between these 2 genes and the likelihood that risk loci exist in each of them. Study of an AA population allowed evaluation of these associations at higher genomic resolution than is possible in a EA population, owing to the much lower LD across these loci in AAs.

Association between polymorphisms in catechol-O-methyltransferase (COMT) and cocaine-induced paranoia in European-American and African-American populations

Catechol-O-methyltransferase (genetic locus, COMT) is a major enzyme involved in catecholamine metabolism and has been associated with numerous psychiatric phenotypes. We studied COMT SNPs and haplotypes in cocaine-induced paranoia (CIP) in African-American (AA) and European-American (EA) populations. We genotyped 17 SNPs across the COMT locus in 319 AA pedigrees (848 individuals) and 302 EA pedigrees (707 individuals). Family-controlled association analyses were conducted using FBAT. We found SNP rs737865 to be nominally significantly associated in the AA family population (P = 0.05). In EAs, the best-known marker, rs4680 (Val158Met), was nominally significant in additive models (P = 0.03). SNP rs174696 also showed nominal significance in additive models (P = 0.02). We considered the three SNPs (rs737866-rs4680-rs174696) together in haplotype analysis in both family populations, using HBAT. The A-A-T haplotype was significantly associated with CIP in EAs (Z = 2.845; P = 0.0044, global P = 0.020). We then studied COMT SNPs in an additional 738 AA and 404 EA unrelated cocaine dependent individuals with and without paranoia. The A-A-T haplotype was significantly associated to CIP in the AA unrelated population (P = 0.0015). Two haplotypes, A-G-C and A-A-C, were significant in the EA unrelated population (P = 0.001 and 0.0003). We also identified rs4680 and three other SNPs, rs933271, rs5993883, and rs740603, as potentially functional variants, as predicted by a signature of positive selection in unrelated EAs and AAs. Based on our robust family-controlled and unrelated-affected analyses, we conclude that COMT is associated with CIP, possibly as a result of its role in the metabolism of dopamine and norepinephrine.

Agitated depression in substance dependence

BACKGROUND

Depression with psychomotor agitation (PMA; "agitated depression") is a putative psychiatric phenotype that appears to associate with some forms of substance dependence. However, it is unclear whether such relationships extend across different substances and independent (I-MDE) versus substance-induced (SI-MDE) subtypes of major depressive episodes.

METHOD

We examined whether lifetime depression with (vs. without) PMA was associated with lifetime substance dependence across individuals with lifetime: (1) I-MDE only (n=575); and (2) SI-MDE only (n=1683). Data were pooled from several family and genetic studies of substance dependence in which participants received identical structured interviews to diagnose DSM-IV mental disorders.

RESULTS

In I-MDE, PMA was significantly associated with alcohol, cocaine, opioid, other drug (hallucinogen, inhalant, speed-ball), and sedative dependence. After controlling for demographic and clinical co-factors, PMA's relationship to dependence on opioids, other drugs, and sedatives remained significant, but not its relationship to alcohol or cocaine. In SI-MDE, PMA was significantly associated with alcohol, cocaine, opioid, and other drug dependence. After adjusting for co-factors, associations remained significant for dependence on cocaine and opioids, but not alcohol or other drugs. Relationships between PMA and opioid dependence were stronger in I-MDE than SI-MDE. Depression subtype (I-MDE vs. SI-MDE) did not moderate relations between PMA and non-opioid forms of substance dependence.

CONCLUSIONS

Agitated depression associates with certain forms of substance dependence, particularly opioid dependence. MDE subtype did not alter most PMA-dependence associations, which suggests that the mechanisms underlying this comorbidity are complex and potentially bidirectional.

Empirically derived subtypes of opioid use and related behaviors

AIMS

To identify and validate homogeneous subtypes of opioid use and related behaviors.

DESIGN

Family-based and case-control genetic studies of opioid and/or cocaine dependence.

SETTINGS

Clinical and general community samples from Connecticut, Massachusetts, Pennsylvania and South Carolina.

PARTICIPANTS

A total of 4061 individuals (2003 individuals from 835 families and 2058 unrelated individuals) recruited to participate in genetic studies.

MEASUREMENTS

The computer-assisted Semi-Structured Assessment for Drug Dependence and Alcoholism (SSADDA) was used to assess participants' demographics, medical history, substance use behaviors and disorders and other psychiatric disorders.

FINDINGS

Five homogeneous subtypes were identified, which differed on opioid-related measures, demographics and prevalence rates of substance use and psychiatric disorders. Heritability estimates for the two most severely affected subtypes exceeded 0.60.

CONCLUSIONS

An empirical approach based on opioid use and related behaviors can yield homogeneous subtypes that could be of value in gene finding for opioid dependence.

Association of CHRNA4 polymorphisms with smoking behavior in two populations

CHRNA4, the gene that encodes the nicotinic acetylcholine receptor α(4) subunit, is a potential candidate gene for nicotine dependence (ND). However, studies of the association of CHNRA4 with smoking behavior have shown inconsistent results. Our meta-analysis of linkage studies of smoking behavior identified a genome-wide significant linkage of the phenotype maximum number of cigarettes smoked in a 24-hour period to a region (20q13.12-q13.32) harboring CHRNA4. This motivated us to examine the association of CHRNA4 with smoking behavior in two independent samples. In this study, we examined five single nucleotide polymorphisms (SNPs) within CHRNA4 and three smoking-related behaviors: one quantitative trait [cigarettes smoked per day (CPD)], and two binary traits [DSM-IV diagnosis of ND and dichotomized Fagerstrom test of ND (FTND)], in 1,249 unrelated European-Americans (EAs) and 1,790 unrelated African-Americans (AAs). Using the combined sample with sex, age, and race as covariates, the synonymous SNP rs1044394 was significantly associated with ND (P = 0.001) and FTND (P = 0.01). Rs2236196, which has a low correlation with rs1044394, was also significantly associated with CPD (P = 0.003). The pattern of association for these SNPs was similar in AAs and EAs. After correction for multiple testing, the association between rs1044394 and ND in the combined sample remained significant (P = 0.033). In summary, our study supports association between CHRNA4 common variation and ND in AA and EA samples. Additional studies will be necessary to evaluate the role of rare variants at CHRNA4 for ND.

Interaction of FKBP5 with childhood adversity on risk for post-traumatic stress disorder

FKBP5 regulates the cortisol-binding affinity and nuclear translocation of the glucocorticoid receptor. Polymorphisms at the FKBP5 locus have been associated with increased recurrence risk of depressive episodes and rapid response to antidepressant treatment. A recent study showed that FKBP5 genotypes moderated the risk of post-traumatic stress disorder (PTSD) symptoms associated with childhood maltreatment. One thousand one hundred forty-three European Americans (EAs) and 1284 African Americans (AAs) recruited for studies of the genetics of substance dependence were also screened for lifetime PTSD. Four single-nucleotide polymorphisms (SNPs) in FKBP5, rs3800373, rs9296158, rs1360780, and rs9470080, were genotyped on the complete sample. Logistic regression analyses were performed to explore the interactive effect of FKBP5 polymorphisms and childhood adversity on the risk for PTSD. After correction for multiple testing, childhood adversity significantly increased the risk for PTSD. FKBP5 genotypes were not associated with the development of the disorder. In AAs, one of the SNPs, rs9470080, moderated the risk of PTSD that was associated with childhood abuse. Without childhood adverse experiences, participants with the TT genotype of this SNP had the lowest risk for PTSD, whereas they had the highest risk for PTSD after childhood adversity exposure. In addition, in EAs, alcohol dependence was observed to interact with childhood adverse experiences, and also FKBP5 polymorphisms, to increase the risk for PTSD. This study provides further evidence of a gene x environment effect of FKBP5 and childhood abuse on the risk for PTSD in AAs. Further study is required in other populations.

Identification of population substructure among Jews using STR markers and dependence on reference populations included

BACKGROUND

Detecting population substructure is a critical issue for association studies of health behaviors and other traits. Whether inherent in the population or an artifact of marker choice, determining aspects of a population's genetic history as potential sources of substructure can aid in design of future genetic studies. Jewish populations, among which association studies are often conducted, have a known history of migrations. As a necessary step in understanding population structure to conduct valid association studies of health behaviors among Israeli Jews, we investigated genetic signatures of this history and quantified substructure to facilitate future investigations of these phenotypes in this population.

RESULTS

Using 32 autosomal STR markers and the program STRUCTURE, we differentiated between Ashkenazi (AJ, N = 135) and non-Ashkenazi (NAJ, N = 226) Jewish populations in the form of Northern and Southern geographic genetic components (AJ north 73%, south 23%, NAJ north 33%, south 60%). The ability to detect substructure within these closely related populations using a small STR panel was contingent on including additional samples representing major continental populations in the analyses.

CONCLUSIONS

Although clustering programs such as STRUCTURE are designed to assign proportions of ancestry to individuals without reference population information, when Jewish samples were analyzed in the absence of proxy parental populations, substructure within Jews was not detected. Generally, for samples with a given grandparental country of birth, STRUCTURE assignment values to Northern, Southern, African and Asian clusters agreed with mitochondrial DNA and Y-chromosomal data from previous studies as well as historical records of migration and intermarriage.