INTRODUCTION

Cocaine dependence is a devastating disorder with no FDA-approved pharmacological treatments available. Genetic studies estimate that 65–78% of the vulnerability risk for cocaine dependence is heritable (Kendler et al, 2000; Kendler and Prescott, 1998); however, identification of genetic risk factors remains difficult due to the complex mode of inheritance, clinical and genetic heterogeneity, and likely multiple genes involved, each contributing only a small effect to the overall risk. Dopaminergic brain systems have been documented to have a major role in drug reward (Dackis and O'Brien, 2005, 2001; Hyman et al, 2006), thus making genes involved in these circuits plausible candidates for susceptibility to substance use disorders (Lachman, 2006). The enzyme catechol-O-methyltransferase (COMT) is involved in the degradation of catecholamines, including dopamine (Axelrod and Tomchick, 1958). The functional COMT polymorphism Val158Met affects enzyme activity, with the Val allele resulting in higher enzyme activity relative to the Met allele (Aksoy et al, 1993; Boudikova et al, 1990; Lachman et al, 1996; Lotta et al, 1995; Spielman and Weinshilboum, 1981). Several studies suggest that the Val158Met polymorphism is involved in psychiatric phenotypes, including schizophrenia and bipolar disorder (Craddock et al, 2006; Tunbridge et al, 2006), and might further contribute to the comorbid substance abuse/dependence spectrum across psychiatric disorders. Although cocaine blocks the dopamine transporter primarily, leading to enhanced postsynaptic effects of dopamine signaling, COMT remains an important regulatory element in dopamine homeostasis. Emerging evidence suggests that COMT variation influences prefrontal cortex (PFC) dopamine regulation and might modulate aspects of cognition, emotions, and behavior (Egan et al, 2001; Tunbridge et al, 2006). PFC dysfunction might be an important component in cocaine dependence that contributes to loss of control and denial (Dackis and O'Brien, 2005). Individual differences in COMT activity might therefore influence vulnerability to cocaine dependence and other substance use disorders. In this study we tested the hypothesis that the functional Val158Met variation of the COMT gene increases susceptibility to cocaine dependence in individuals of African descent.

MATERIALS AND METHODS

DNA samples from cocaine-dependent individuals of African-American descent (n=330; 72% male subjects, mean age: 43) were collected during clinical studies of cocaine dependence at the University of Pennsylvania Treatment Research Center. Subjects were at least 18 years of age and were all assessed with the Structured Clinical Interview for DSM Disorders and urine drug screens were obtained. All patients had a clinical diagnosis of cocaine dependence as defined by DSM-IV. Family history was not obtained and ethnicity was determined by self-report. Exclusion criteria were all psychiatric axis I disorders except alcohol dependence/abuse and nicotine dependence. In addition to these exclusion criteria, participants were excluded if they had a history of a seizure disorder (except cocaine-induced seizures) or a severe medical illness, including a history of AIDS (but not merely of HIV+ status). Individuals currently being treated with psychotropic medications or with psychiatric symptoms, including psychosis, dementia, suicidal or homicidal ideation, mania or depression requiring antidepressant therapy were excluded from these studies. Control samples from persons of African descent (n=255; 29% male subjects, mean age: 40) were collected at the University of Pennsylvania, Thomas Jefferson University and through the National Institute of Mental Health Genetics Initiative (www.nimhgenetics.org). Control individuals were screened for history of substance use disorders or other psychiatric illness. Control subjects were not assessed with a urine drug screen and ethnicity determination was by self-report. Subjects with a history of substance dependence or a history of major psychiatric illness (schizophrenia and unipolar or bipolar illnesses) as defined by DSM-IV criteria were excluded from this study (Berrettini and Persico, 1996). All protocols were approved by the Institutional Review Boards at Thomas Jefferson University and the University of Pennsylvania and all subjects provided written informed consent before DNA sample collection.

Genotyping of three SNPs across the COMT gene was performed using the Applied Biosystems Inc. (Foster City, CA, USA) ‘Assays-on-demand’ (ABI) SNP genotyping assay as per manufacturers protocol (SNP1: rs737865; SNP2: rs4680 (Val158Met); SNP3: rs165599). SNPs were selected based on findings in previous studies (Shifman et al, 2002), including studies of nicotine dependence in African-American subjects (Berrettini et al, 2007) (Figure 1). Genotyping quality control was assured by genotyping 10% duplicates for cases and controls. Concordance rate of genotypes was >99.5%.

Figure 1
figure 1

Diagram of the COMT gene and the three SNPs examined, with linkage disequilibrium D′ values. Values in bp refer to chromosome 22, taken from the March 2006 assembly of the human genome sequence at www.genome.ucsc.edu. Exons are indicated by shaded boxes. The unshaded box represents the 3′ UTR in exon 6. D′ values represent African-American sample results.

Statistical Analyses

Genotypes and allele frequencies were compared between groups using χ2 contingency analysis. A two-tailed type I error rate of 5% was chosen for the analysis. Linkage disequilibrium (LD) and haplotype frequencies were estimated using the COCAPHASE program (Dudbridge, 2003). The COCAPHASE program uses standard unconditional logistic regression analysis. Correction for multiple testing was performed using permutation correction by the COCAPHASE program. This approach corrects for multiple testing but takes into account the correlation between markers. It is thus less conservative than a Bonferroni correction, which is appropriate for independent tests such as unlinked markers. For the single-marker analyses, 10 000 permutations were carried out to estimate the significance of the best results, correcting for the three loci tested. Haplotype analysis was performed using a two- and three-marker window. Rare haplotypes were excluded from analysis as the EM algorithm does not accurately estimate haplotype frequencies <1% (Fallin and Schork, 2000). The most significant p-value was corrected by permutation analysis as described above.

Our sample size had reasonable power to detect a disease association at a p-value less than or equal to 0.05, assuming an odds ratio of 1.5 and a minor allele frequency of 30% (99% for a log-additive mode of inheritance, 92% for a dominant, and 56% for a recessive mode of inheritance). Power analysis was performed using the Quanto program (Gauderman, 2002).

RESULTS

None of the genotype counts deviated significantly from those expected from Hardy–Weinberg equilibrium for cases or controls. Genotype and allele frequencies differed significantly for the Val158Met polymorphism between cocaine-dependent individuals (f(Met)=35%) and normal controls (f(Met)=27%) (p=0.004; corrected p=0.014; OR 1.44; 95% CI 1.12–1.86). SNP1 and SNP3 did not show a statistical difference between cases and controls (Table 1). There was no significant gender effect for any of the tested SNPs. Allele frequencies were consistent with those reported in the literature for individuals of African descent (Berrettini et al, 2007). Haplotype analysis showed significant associations for two-marker analysis and a trend for the three-marker combination (Table 2); however, after correction for multiple testing only the rs737865–rs4680 haplotype remained statistically significant. The patient group carried the major allele of rs737865 and the Met158 allele more often (33%) compared to controls (26%) (p=0.005; corrected p=0.02; OR 1.44). Haplotype analysis results appear to be driven by the Val158Met SNP, supporting the hypothesis that the Val158Met polymorphism is a causative functional SNP.

Table 1 Genotype and Allele Frequencies of Variations in the COMT Gene
Table 2 Analysis of Common Haplotypes in the COMT Gene

DISCUSSION

In the present study, we show an association between the Val158Met polymorphism of the COMT gene and cocaine dependence in individuals of African descent. Furthermore, we identify a risk haplotype contributing to the susceptibility for cocaine dependence. To our knowledge, this is the first report of an association of the COMT Val158Met polymorphism in cocaine dependence in individuals of African descent; however, given the high comorbidity of polysubstance abuse in individuals using cocaine, our result might reflect an association of a broader phenotype of substance use disorders. In fact, several previous studies have implicated the Val158Met polymorphism in a variety of substance abuse disorders (Table 3) (Beuten et al, 2005; Enoch et al, 2006; Horowitz et al, 2000; Hosak et al, 2006; Li et al, 2004; Samochowiec et al, 2006; Sery et al, 2006; Tiihonen et al, 1999; Vandenbergh et al, 2000, 1997; Wang et al, 2001); however, others could not replicate results (Cevoli et al, 2006; Guo et al, 2007; Hallikainen et al, 2000; Kauhanen et al, 2000; Kweon et al, 2005) and there is no clear consensus on whether the Val allele or the Met allele increases risk (Table 3). This locus heterogeneity might indicate that the COMT Val158Met polymorphism is only one important variation in a cascade of regulatory systems, and depending on other genes or environmental factors, a higher or lower enzyme activity might predispose to substance use disorders. The COMT Val158Met polymorphism, rather then a single allele, might thus have an effect on pathways commonly shared between all substance use disorders. Such pathways include the reward system and cognitive functions influencing substance use behavior.

Table 3 Case–Control Association Studies of the COMT Val158Met Polymorphism in Substance Use Disorders

Dysregulation of COMT, which is the major enzyme involved in the degradation of dopamine in the frontal cortex (Karoum et al, 1994), might have primarily effects on cognitive processes involved in substance dependence and indirect downstream effects on the reward system, as suggested by recent neuroimaging data (Goldstein and Volkow, 2002; Volkow et al, 2002). Our results suggest an increased Met allele (low activity allele) frequency in cocaine users (35%) compared to normal controls (27%) (p=0.004; corrected p=0.014) in individuals of African descent. This finding is in line with the results of Hosak et al (2006) who showed higher novelty seeking scores in methamphetamine-dependent individuals with a Met allele (Hosak et al, 2006) and the observation that Met/Met carriers have decreased efficiency of PFC information processing in response amphetamine (Mattay et al, 2003). Individuals with the low-activity COMT allele may have longer-lasting and more effective dopamine release in the brain, in particular in the PFC affecting higher cortical functions involved in substance use behavior.

Although our study provides evidence for an association between the Val158Met polymorphism and cocaine dependence, it could be possible that other variations that are in LD with the Val158Met SNP might contribute to the observed association or that a haplotype confers risk rather then a single SNP. Consistent with this possibility are several studies that found haplotypes to be associated with disease rather than with a single polymorphism (Berrettini et al, 2007; Shifman et al, 2002). Recent reports show that COMT haplotypes code for differences in enzyme activity (Diatchenko et al, 2005) and haplotype-specific mRNA secondary structure has functional effects on COMT protein synthesis and enzyme activity (Nackley et al, 2006). Our haplotype analysis indicates a risk haplotype for rs737865 and Val158Met (p=0.005; corrected p=0.02; OR 1.44). Interestingly, both markers (rs737865 and Val158Met) are in strong LD and incorporate the SNPs in the functional haplotypes described by Nackley et al (2006). Additional studies are necessary to elucidate potential functional haplotypic variations. Furthermore, it might be necessary to investigate this chromosomal region in more depth, given recent evidence of copy-number variations in this region (Wong et al, 2007) and the close proximity of COMT to other interesting candidate genes, like ARVCF, a member of the catenin family, which is involved in cell–cell communication.

Even though we report a positive association between the COMT gene and cocaine dependence, it is possible that our finding might be a false-positive result due to population stratification. Case–control association studies of subjects with self-reported ancestries are not immune to population stratification (Freedman et al, 2004). Undetected differences in population structure can mimic the signal of association and lead to false-positive results or real effects that are missed (Pritchard and Donnelly, 2001). This is a particular concern for analyses of samples of African-American descent, as recent studies indicate larger genetic admixture than previously thought (Parra et al, 2001; Pfaff et al, 2001; Tian et al, 2006; Zhu et al, 2004). In fact, the Met allele is less frequent in individuals of African descent (Ameyaw et al, 2000; DeMille et al, 2002; McLeod et al, 1994), and haplotypes in COMT show marked differences across populations (Palmatier et al, 1999). Possible strategies to control for these stratification issues are the use of genomic controls (Bacanu et al, 2000; Devlin and Roeder, 1999) and/or the use of a family-based association design, a strategy that matches the genotype of an affected offspring with parental alleles not inherited by the offspring (Spielman and Ewens, 1996). In addition to these potential issues of genetic heterogeneity and population stratification, it is also important to consider limitations of clinical heterogeneity. All patients were diagnosed according to DSM-IV criteria; however, comorbid use of alcohol and nicotine might have differed between patients. Furthermore, the control subjects were assessed using semistructured interviews but did not undergo urine drug testing. Although drug testing is useful in establishing a diagnosis, it might not be useful for assessment of controls as it does not rule out past exposure or substance use. Unreported or minimized substance abuse in the control population is thus an important limitation that needs to be considered. Finally, spurious positive association findings remain a valid concern as recently shown in a statistical simulation study of the COMT gene (Sullivan, 2007). Thus, our results should be interpreted with caution and ultimately require careful replication and confirmation in an independent population of patients and controls.

In summary, we show that the Val158Met polymorphism in the COMT gene is associated with cocaine dependence. In addition we have identified a COMT risk haplotype for cocaine dependence. Our results require confirmation in other populations, and additional studies are required to elucidate the role of COMT in the pathophysiology of substance use disorders.