Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Drug Alcohol Depend. Author manuscript; available in PMC 2009 January 11.
Published in final edited form as:
PMCID: PMC2266629

Linkage scan for quantitative traits identifies new regions of interest for substance dependence in the Collaborative Study on the Genetics of Alcoholism (COGA) Sample


Dependence on alcohol and illicit drugs frequently co-occur. Results from a number of twin studies suggest that heritable influences on alcohol dependence and drug dependence may substantially overlap. Using large, genetically informative pedigrees from the Collaborative Study on the Genetics of Alcoholism (COGA), we performed quantitative linkage analyses using a panel of 1717 SNPs. Genome-wide linkage analyses were conducted for quantitative measures of DSM-IV alcohol dependence criteria, cannabis dependence criteria and dependence criteria across any illicit drug (including cannabis) individually and in combination as an average score across alcohol and illicit drug dependence criteria. For alcohol dependence, LOD scores exceeding 2.0 were noted on chromosome 1 (2.0 at 213 cM), 2 (3.4 at 234 cM) and 10 (3.7 at 60 cM). For cannabis dependence, a maximum LOD of 1.9 was noted at 95 cM on chromosome 14. For any illicit drug dependence, LODs of 2.0 and 2.4 were observed on chromosome 10 (116 cM) and 13 (64 cM) respectively. Finally, the combined alcohol and/or drug dependence symptoms yielded LODs > 2.0 on chromosome 2 (3.2, 234 cM), 10 (2.4 and 2.6 at 60 cM and 116 cM) and 13 (2.1 at 64 cM). These regions may harbor genes that contribute to the biological basis of alcohol and drug dependence.

Keywords: Linkage, alcohol, cannabis, illicit drugs, dependence, COGA


Twin studies have demonstrated that heritable influences play a prominent role in the etiology of both alcohol and illicit drug dependence. In adults, genetic factors contribute to 40-60% of the total variance in risk for alcohol dependence (Dick et al., 2006; Heath et al., 1997; Kendler et al., 1994; Prescott et al., 1995; Tsuang et al., 2001), while for illegal drug dependence, estimates of heritability range from 35-75% (Kendler et al., 2003a; Lynskey et al., 2002; Tsuang et al., 1996; Tsuang et al., 2001; van den Bree et al., 1998). Results from multivariate twin analyses suggest that (i) the liability for dependence on various illicit drug classes is governed, in part, by common genetic factors (Kendler et al., 2003a; Tsuang et al., 1998), and (ii) a significant proportion of the genetic influences on alcohol and illicit drug dependence may be overlapping (Bierut et al., 1998; Kendler et al., 2003b; McGue et al., 2000; Tsuang et al., 2001).

Linkage studies of alcohol dependence have identified several chromosomal regions, including chromosomes 1, 2, 4, 7 and 11 (see Dick et al., 2006 for a review). Gelernter et al. (2005, 2006) have found evidence for linkage on chromosomes 2, 3, 9, 11, 12, 17 and 18 for various definitions of cocaine and opioid dependence. Recently, Hopfer and colleagues (2006) reported linkage on chromosomes 3 and 9 for cannabis dependence criteria in adolescents. Stallings et al. (2003) performed linkage analyses, in the same sample of adolescent probands and controls for average number of dependence criteria across alcohol, tobacco and illicit drugs, and reported elevations on chromosomes 3 and 9. Uhl and colleagues (2001) used 1,494 single nucleotide polymorphisms (SNPs) to examine genome-wide association for general substance dependence vulnerability in pooled samples from a case-control study. Evidence for association was reported for regions on chromosomes 3, 4, 9, 11 and 13 in both European-Americans and African-Americans (Uhl et al., 2001). Subsequently, these authors conducted genome-wide association analyses using 639,401 SNPs in pooled samples to identify a number of cell-adhesion genes (e.g. CNTN4, CNTN5, CNTN6: Contactin genes) associated with a general vulnerability to substance abuse/dependence (Liu et al., 2006).

A number of linkage studies of alcohol or substance dependence have used dichotomous phenotypes (i.e. sib pairs affected or unaffected for alcohol dependence). Such a dichotomous phenotype is limited in power and relies on diagnostic thresholds (e.g. a diagnosis of dependence upon endorsement of three or more criteria), which have repeatedly been demonstrated to be methodologically limited (Helzer et al., 2006a; Helzer et al., 2006b; Saha et al., 2006). Alternatively, ordinal measures (e.g. dependence criteria), that are adjusted for influential covariates (e.g. age, gender, ethnicity), may be more informative phenotypes, especially as they are not restricted by diagnostic thresholds (e.g. affected if and only if atleast 3 criteria endorsed). Unlike an analysis of affected sibling pairs, ordinal and continuous measures allow for incorporation of phenotypic data from all eligible family members, and thus better utilize the entire sample.

In the present study, we used data from the Collaborative Study on the Genetics of Alcoholism (COGA) (Begleiter et al., 1995; Reich et al., 1996,1998) to perform regression-based linkage analyses (Sham et al. 2002) for DSM-IV criteria of alcohol dependence, cannabis dependence and any illicit drug dependence (including cannabis). However, instead of a microsatellite marker panel, we selected a panel of 1,717 SNPs (Edenberg et al., 2005), which yielded a significant increase in information content across the genome.


2.1 Sample for GAW14

Genotyping of a panel of single-nucleotide polymorphisms (SNPs) was performed on a subset of COGA families determined to be most informative for linkage analyses (Edenberg et al., 2005). In general, the COGA high-risk families consist of multigenerational pedigrees ascertained using probands who met diagnostic criteria for both DSM-IIIR alcohol dependence (American Psychiatric Association, 1987) and for definite alcoholism specified by Feighner et al. (1972). Multiplex alcoholic families that were not bilineal and had at least two affected first-degree relatives in addition to the proband, along with informative second and third degree relatives comprise the original COGA high-risk genetic sample (Nurnberger et al., 2004). For SNP-typing, a sub-set of 143 families with at least 6 members with genotypic and interview data, a total of 1,614 individuals, was drawn from the high risk genetic sample. Of these 1,614 individuals, 1,364 were genotyped with the SNP panel, with the remaining 250 family members included to link members of the family in a genotypically informative manner (i.e. pedigree cohesion). Further details can be found in Edenberg et al. (2005).

The full COGA sample also consists of 984 individuals from the community recruited using a variety of methods including driver’s license records, medical/dental clinics, mailed questionnaires and advertisements, who were neither ascertained nor excluded for psychiatric or substance-related disorders. This community sample was interviewed using the same diagnostic instruments as the high-risk sample; however, no genotypic information is available on these individuals.

2.2 SNP Panel

We used a thinned panel of the original 4,596 SNPs typed by Illumina. Details surrounding the genotyping procedure may be found elsewhere (Edenberg et al., 2005). COGA pedigrees selected for SNP-typing were analyzed for Mendelian inconsistencies and genetic relatedness, including issues of non-paternity, by two independent groups during the GAW14 workshop (Wang et al., 2005; Hinrichs & Suarez, 2005). Final genotypic cleaning was performed using PREST (McPeek & Sun, 2000). Here, we focus on the procedure used to account for the substantial inter-marker linkage disequilibrium in the original Illumina panel and the consequent thinning procedure. Huang et al., (2004) have shown that in the absence of parental genotypes, linkage disequilibrium can produce spurious inflations in IBD estimates and inflate information content as a consequence. While a proportion of the COGA parental generations have been genotyped, we opted to utilize a thinning procedure that would both reduce the likelihood of IBD inflation and would retain the information content afforded by the full panel of SNPs. Thinning was performed by deleting SNPs with r2 of 0.1 or greater with any other SNP within 1 Mb. This produced our final thinned map of 1,717 SNPs and constitutes the panel used for the present analyses. Whenever possible, SNPs with the highest minor allele frequency were retained. The thinned map provided similar information content across the genome when compared to the full panel of SNPs (Hinrichs et al., 2005a; Hinrichs et al., 2005b). In general, information content was high (between 80-95%) across the genome. For larger chromosomes, such as chromosome 1, a modest dip (to between 50-65%) in information was noted in telomeric regions. Further, as shown by Hinrichs and colleagues (2005a), using chromosome 7 as an example, the dense and thinned maps afforded identical information, ranging from 75-95% (with lower information at telomeres), across this chromosome, which was significantly greater than information afforded by microsatellites (range 55-90%).

2.3 Phenotypes

Interview data on alcohol and illicit drug dependence criteria were collected using the Semi-Structured Assessment for the Genetics of Alcoholism (Bucholz et al., 1994; Hesselbrock et al., 1999). DSM-IV dependence criteria were coded for alcohol and individual illicit drugs including cannabis, cocaine, sedatives, stimulants and opiates (American Psychiatric Association, 1994). This included 7 criteria (tolerance, using larger quantities/longer periods of time, spending time acquiring/using, multiple attempts to quit/cut back, recurrent use despite physiological/emotional problems, giving up or cutting down important activities to use the substance and withdrawal for each substance, except cannabis). The following phenotypes were log-transformed (using “log (n+1)” transformation, where n is the raw variable with addition of “1” to account for zeroes) and used for linkage analyses:

  • (i) DSM-IV alcohol dependence criteria (range 0-7);
  • (ii) DSM-IV cannabis dependence criteria (range 0-6; withdrawal syndrome syndrome excluded to conform to DSM-IV);
  • (iii) Average of DSM-IV dependence criteria across all illicit drugs (cannabis, cocaine, sedatives, stimulants and opiates) (range 0-7);
  • (iv) Average of DSM-IV dependence criteria across alcohol and illicit drugs (range 0-7);

Criterion counts were missing in those who reported never having used the substance even once in their lifetime. Average scores were created by dividing the total number of criteria by the number of substances used. Stallings et al. (2003) have previously demonstrated that this assessment of dependence vulnerability is highly heritable (h2=48%) and also least likely to be influenced by environmental influences shared by siblings.

2.4 Sample Characteristics

A total of 1,214 Caucasians (N=112 families, average family size of 10.8) and 150 African-Americans (N=11 families, average family size of 13.6), were used for linkage analyses. The mean age of participants was 40.6 years (range 17-91years). Table 1 shows the distribution of dependence criteria in 984 individuals from the community sample and 2,773 individuals from the high risk sample: 13 and 46.4% of the community and high risk participants endorsed 3 or more alcohol dependence criteria respectively. Similarly 3.1 and 11.5% of participants endorsed 3 or more cannabis dependence criteria in community and high risk sample respectively. The mean number of alcohol dependence criteria was 0.9 and 2.7 in the community and high risk samples, with the comparable estimates for mean number of cannabis dependence criteria being 0.4 and 1.3. The mean number of DSM-IV illicit drug dependence criteria was 0.7 for the community sample and 3.9 for the high risk sample. Mean number of DSM-IV alcohol and illicit drug dependence criteria was 1.2 and 5.4 for the community and high risk samples respectively.

Table 1
Distribution of DSM-IV alcohol and drug dependence criteria in a community (N=984) and high-risk (N=2773, 13% African-American) sample.

2.5 Linkage Analysis

We used the regression-based quantitative trait linkage analysis available through MERLIN-REGRESS (Sham et al., 2002). While this method is, in general, robust to distributional assumptions of the trait, it is sensitive to ascertainment effects (phenotypic) and to differences in allele frequencies (genotypic) across populations. To account for these effects, we employed the following analytic steps:

  • (i) Using the community sample with neither ascertainment nor exclusion criteria, we determined regression parameters for significant confounding covariates, including gender, age (linear and/or quadratic effects, only significant parameter retained) and race (Caucasian or African-American).
  • (ii) These parameters (intercept and beta coefficients) for significant covariates were then applied to adjust the trait scores in the high risk families. Therefore, the resulting trait ytrait = xtrait − (Intercept + β1*covariate1 + β2*covariate2 + β3*covariate3), where xtrait is the substance dependence phenotype in the ascertained sample.
  • (iii) Families were identified as Caucasian or African-American (determined by self-report, confirmed using STRUCTURE (Kauwe et al., 2005)) and separate pedigrees for the Caucasian and African-American families. The race-specific pedigrees were combined using the pedmerge option available in the MERLIN package. If markers are named uniquely in each sub-population (say, wSNP1 in Caucasians and correspondingly, bSNP1 in African-Americans), pedmerge appends the pedigrees from both populations while retaining the population-specific allele frequencies.
  • (iv) MERLIN-REGRESS models for multipoint linkage analyses were specified using means and variances based on residuals from the regression analyses in the community sample from (i). Means, therefore, were near zero in all instances. MERLIN-REGRESS also requires an estimate of heritability, which we fixed at 0.50 for these analyses. This heritability estimate is well within the range reported in the literature.

2.6 Genomewide significance

To estimate levels for genomewide significance, 1000 Caucasian and African-American pedigrees were simulated using the gene-dropping algorithm available from MERLIN (Abecasis et al., 2002). Simulations assign random chromosomes to founders and segregate these alleles through the pedigrees while retaining the phenotypic data and segregation patterns. Pedigrees were combined after simulation using pedmerge and multipoint linkage analyses were performed. Maximum LOD scores across simulations were collected and the empirical p-value [p with macron] = r + 1/1001, where r was the number of times a simulated LOD score exceeded an observed LOD score value (North et al., 2003).


3.1 Multipoint Linkage Analyses

For DSM-IV alcohol dependence criteria, we observed linkage peaks on chromosomes 1 (LOD 2.0, 213 cM), 2 (LOD 3.4, 234 cM), 10 (LOD 3.7, 60 cM) and 13 (LOD 1.4, 64 cM) (see Table 2). For cannabis dependence criteria, a maximum LOD score of 1.9 was noted at 95 cM on chromosome 14, with a LOD of 1.4 on chromosome 13, coinciding with the LOD of similar magnitude for alcohol dependence at the same position. This elevation on chromosome 13 was considerably higher for any illicit drug dependence (LOD 2.4, 64 cM). In addition, for any illicit drug dependence, a LOD of 2.0 was noted at 116 cM on chromosome 10. When criteria of alcohol dependence and drug dependence were averaged, two regions of the genome showed significant linkage: on chromosome 2, a LOD of 3.2 at 234 cM and on chromosome 10, a LOD of 2.4 at 60 cM. Both linkage peaks overlapped with linkage peaks for DSM-IV alcohol dependence criteria (see Figure 2). Other evidence for linkage with the combined phenotype included chromosome 1 (LOD 1.4, 213 cM) and chromosome 13 (LOD 2.1, 64 cM). Empirical p-values based on a 1,000 multipoint simulations suggest that our linkage peaks are suggestive (p-values of 0.10 or less) for LOD scores exceeding 3.4, with p-value of 0.05 for LOD scores of 3.7 and greater.

Table 2
Multipoint LOD scores (exceeding 1.0) across DSM-IV dependence criteria of alcohol and illicit drugs.


Using an informative panel of single nucleotide polymorphisms (SNPs) that were typed in a subset of individuals from the Collaborative Study on the Genetics of Alcoholism (COGA), we conducted multipoint linkage analyses for DSM-IV criteria of alcohol and drug dependence. Linkage peaks on chromosome 2 and 10 were noted for alcohol dependence and for an average score of alcohol and drug dependence.

Prior linkage studies, especially those utilizing COGA data and a categorical phenotype (i.e. affected sibling pairs), have noted linkage peaks on chromosome 2 for a variety of substance-related phenotypes but these peaks have converged to 120-135 cM (Bierut et al., 2004; Dick et al., 2004; Foroud et al., 2000; Hesselbrock et al., 2004; Reich et al., 1998). Our findings, in contrast, occur at 230-260 cM, where Schuckit et al., (Schuckit et al., 2001) report linkage for SRE (Subjective Response to Ethanol) scores during the first five times alcohol was consumed (FIRST 5) and by Nurnberger et al. (2001) for comorbid alcoholism and depression. In addition, Straub et al. (1999) report linkage in this 2q region for nicotine dependence in their Christchurch sample. Gelernter et al., (2005, 2006) have also reported linkage for late-onset cocaine dependence and for DSM-IV opioid dependence at 221 cM on chromosome 2, but only in their African-American pedigrees. Zubenko et al., (2003a, b) find linkage in this region for Depression Spectrum Disorder, and recently, Kuo et al., (2006) report linkage for an alcohol withdrawal criteria factor score. While there is considerable phenotypic, and potentially genetic overlap, across the phenotypes assessed in our manuscript and those used in the above referenced studies, it is also possible that imprecise localization of linkage signals, or random chance, may have contributed to convergence of linkage signals.

Our highest LOD score for alcohol dependence was on chromosome 10. The elevations at 60 and 116 cM are in regions where several studies have reported linkage for correlated psychiatric traits. Li and colleagues (2006) report a LOD of 4.17 for a peak ranging from 60-100 cM for number of cigarettes smoked per day in a sample of African-American families ascertained for nicotine dependence. In contrast, Uhl et al., 2001, in their genomewide association study found evidence for association between a polymorphism (WIAF-3336) at 86 cM on chromosome 10 and substance dependence vulnerability in their European-American cases and controls alone. While the African-American families contributed to our LOD score on chromosome 10, re-doing the linkage analyses without these families did not greatly reduce the LOD score (LOD was reduced to 3.5), suggesting that our finding on chromosome 10 is not specific to African-American families. Vink et al. (2004) showed LOD scores of similar magnitude of both smoking initiation and quantity at 50-60 cM on chromosome 10. Strikingly, Zubenko and colleagues have reported linkage for mood disorders at 76 and 113 cM in their family study of recurrent unipolar major depressive disorder (Zubenko et al., 2003b). The latter peak at 116 cM also maps to within 10 cM of the region where Schuckit et al, (2005) report linkage for the Subjective High Assessment Scale (SHAS), a measure of alcohol response.

We did not, however, find significant evidence for linkage peaks on chromosome 3, 9 or 17, where Stallings et al., (2003) report linkage for average scores of alcohol and drug dependence, as well as for cannabis dependence, in their sample of adolescent probands from treatment samples, and matched controls. In fact, our highest LOD score for cannabis dependence criteria (LOD 1.92) occurred on chromosome 14. There may be several reasons for this including different ascertainment strategies between the two studies (i.e. COGA used alcohol dependent probands and dense families, while Stallings et al., used probands with substance use or conduct problems) and the higher information content afforded by using SNP panels instead of microsatellites.

It is also worth noting that some linkage peaks previously identified in COGA on chromosome 4 and 7 (Saccone et al., 2000; Saccone et al., 2005; Foroud et al., 2000; Reich et al., 1998; Nurnberger, Jr. et al., 2001) were not identified in the current study. Those peaks harbor genes such as GABRA2, ADH (on chromosome 4) and CHRM2 and hTAS2R16 (on chromosome 7), which have subsequently found to be associated with risk for alcoholism (Edenberg et al.,2004; Wang et al., 2004; Edenberg et al., 2006a; Hinrichs et al., 2006). We did not detect these peaks in our analyses. In other samples, however, quantitative indices of alcohol dependence similar to those used here have yielded linkage on chromosome 4 (see Prescott et al. 2006 for summary). For instance, in the Irish high-density sample of alcoholics, Prescott and colleagues (2006) report linkage (LOD 4.5) on chromosome 4. Whether the lack of these linkage findings in the present study is related to the phenotype used here or is a false negative is unknown. While the non-replication poses some interpretive challenges, it is not unusual in the extant linkage literature for substance-related phenotypes.

Before considering plausible candidates in the genomic regions encompassing our highest LOD scores, some study limitations need to be considered. First, there is significant overlap between the high risk families that contributed to the linkage signal for alcohol dependence criteria and the families for illicit drug dependence criteria. In our sample, nearly half of those with alcoholism also meet criteria for DSM-IV drug dependence and less than 25% of those with drug dependence do not report a lifetime history of alcohol dependence. Therefore, the overlap of linkage findings is likely due to the high comorbidity in this sample. Second, due to the relatively small proportion of African-American families, we did not have sufficient power to conduct linkage analyses independently in these families. Differences in allele frequencies, however, were maintained by creating pedigrees in each racial group and using pedmerge to combine them. Third, the community-based sample was ascertained from various sources: driver’s license registries, dental clinics and health maintenance organizations. While these families were not enriched for substance dependence or psychiatric disorders, they were selected to be large, therefore they may provide limited representation of the general population. Notwithstanding this limitation, using community-based estimates is highly recommended for regression-based linkage analyses. Furthermore, re-doing the analyses without using the community sample (as discussed below) did not significantly alter our findings. Lastly, diagnostic criteria for nicotine dependence were not included in the initial COGA assessments, and hence dependence vulnerability to tobacco was not included in these analyses.

Several subsidiary linkage analyses were conducted to validate the robustness of our LOD scores. First, analyses were performed using means and variances from the high-risk sample (i.e. without scoring) and without regressing out the influence of covariates. Second, analyses were replicated using varying estimates of heritability (0.45-0.55). Third, analyses were also conducted using DSM-IIIR dependence criteria for alcohol and illicit drugs. Fourth, linkage analyses were repeated in the full panel of 4,596 SNPs and were also conducted using the original microsatellite panel. Across all subsidiary analyses, LOD scores remained consistent, although modest magnitude changes were noted in some instances. In addition, LOD scores remained unchanged when particular large families were trimmed from the analysis.

Our strongest evidence for linkage was observed on 2q and 10p-q. On the q-arm of chromosome 2, putative candidate genes include HTR2B (5-Hydroxytryptamine receptor 2B) receptors. Recently, Lin et al. (2004) reported association between 3 coding SNPs (2 non-synonymous, resulting in double-mutants of the protein, and one synonymous) in HTR2B and substance abuse vulnerability. Other possible candidate genes include GPR55 (G protein-coupled receptor 55), which is expressed in human brain tissue and involved in extracellular-intracellular signal transduction, and HDAC4 (histone deacetylase 4) which participates in epigenetic modification of core histones and. Work by Zubenko and colleagues has also suggested a role of CREB1 on mood disorders; this lies 20 cM centromeric to our linkage peak (Zubenko et al., 2003a). These authors have also identified a sex-specific association between the 124bp repeat allele of D2S944 and recurrent, early-onset major depression and with anxious depression in an independent U.S. and Dutch sample (Beem et al., 2006; Philibert et al., 2003; Zubenko et al., 2002). Due to the high level of comorbidity between alcohol and drug dependence and major depression, and as indexed by a high LOD score of 3.49 for comorbid alcoholism and major depression in prior linkage analyses using COGA data (Nurnberger, Jr. et al., 2001), this region on chromosome 2 may be tapping into a region of susceptibility for both disorders.

The linkage peak on chromosome 10 spans from 42 to 65 cM. GAD2 (glutamic acid decarboxylase 2), a possible candidate in this region, has been putatively implicated, although not unequivocally, for its role in acute ethanol consumption and withdrawal, in animal models (Fehr et al., 2003). Two recent association studies in humans, however report negative findings for the role of GAD2 in the risk for alcoholism while another found limited evidence for the role of this gene in association with anxiety disorders, depression and neuroticism (Hettema et al., 2006; Lappalainen et al., 2007; Loh et al., 2006). A smaller linkage peak at 90-110 cM encompasses candidate genes such as HTR7 (5-Hydroxytryptamine receptor 2B) and VMAT2 (vesicular monoamine oxidase 2), as well as KCNMA1 (calcium activated channel). SNPs in VMAT2 and KCNMA1 have been shown to be associated with alcoholism and the endophenotype of subjective responses to ethanol (Lin et al., 2005; Schuckit et al., 2005). Also located at 83 cM is CTNNA3 (catenin alpha 3) which was recently identified in genomewide association studies of substance dependence vulnerability and nicotine dependence (Bierut et al., 2007; Liu et al., 2006).

Finally, our finding on chromosome 14 for cannabis dependence, while modest, is in a region of biological interest. Our linkage peak on chromosome 14 harbors candidate genes, such as GPR68 (G-protein coupled receptor 68) involved in cAMP regulation, CKB (creatine kinase, brain), encoding a cytoplasmic enzyme involved in energy homeostasis, as well as SERPINA1 and SERPINA2 (serine-peptidase inhibitor, clade A, members 1 and 2), which were previously identified in the genomewide association study for substance abuse vulnerability by Liu and colleagues (and also by Bierut et al.,61 for SERPINA1 and nicotine dependence).

In the past, linkage analyses in COGA have been extremely successful in identifying putative candidate genes which have led to positive association results (Dick et al., 2006; Edenberg et al., 2006b). We have now identified additional genomic regions that may harbor genetic loci that contribute to the etiology of alcohol and illicit drug dependence. Future efforts will target candidate genes in these regions for association analyses. This continued effort at understanding the biological basis of alcohol and drug use disorders is critical as, even today, the challenge posed by the considerable morbidity and mortality due to alcohol and drug problems is a substantial concern (Centers for Disease Control (CDC), 2004; Williams et al., 1988).

Multipoint linkage for criteria of DSM-IV alcohol and cannabis dependence, and for average number of illicit drug dependence and alcohol and illicit drug dependence criteria in the Collaborative Study of the Genetics of Alcoholism (COGA).


In memory of Henri Begleiter and Theodore Reich, Principal and Co-Principal Investigators of COGA since its inception; we are indebted to their leadership in the establishment and nurturing of COGA, and acknowledge with great admiration their seminal scientific contributions to the field.

Acknowledgement and Role of Funding Source: This national collaborative study is supported by the NIH Grant U10AA008401 from the National Institute on Alcohol Abuse and Alcoholism (NIAAA) and the National Institute on Drug Abuse (NIDA). A.A. is also funded by DA023668. L.J.B. is supported by DA13423, DA019963 & PO1 CA89392. SNP data were obtained through participation in a collaborative effort between COGA, organizers of the 14th Genetic Analysis Workshop (GAW), the Center for Inherited Disease Research (CIDR: CIDR is fully funded through a federal contract from the National Institutes of Health to the Johns Hopkins University, contract number N01-HG-65403), Illumina and Affymetrix. Funding sources were not involved in conducting these analyses.


Contributors: The Collaborative Study on the Genetics of Alcoholism (COGA), Co-Principal Investigators B. Porjesz, V. Hesselbrock, H. Edenberg, L. Bierut, includes nine different centers where data collection, analysis, and storage take place. The nine sites and Principal Investigators and Co-Investigators are: University of Connecticut (V. Hesselbrock); Indiana University (H.J. Edenberg, J. Nurnberger Jr., P.M. Conneally, T. Foroud); University of Iowa (S. Kuperman); SUNY Downstate (B. Porjesz); Washington University in St. Louis (L. Bierut, A. Goate, J. Rice); University of California at San Diego (M. Schuckit); Howard University (R. Taylor); Rutgers University (J. Tischfield); Southwest Foundation (L. Almasy). Zhaoxia Ren serves as the NIAAA Staff Collaborator.

All authors provided feedback and revisions and approved the manuscript for submission. In addition to data collection (see above), individual contributions to this paper, in alphabetical order, are as follows: Agrawal: Conducted analyses, responsible for first draft and incorporating feedback into all subsequent drafts; Bertelsen: Conducted analyses; Bierut: Phenotypic & Bioinformatics expertise; on-site supervision; Bucholz: Phenotypic expertise; Cloninger: Phenotypic and statistical expertise; Dick: Phenotypic and statistical expertise; Dunn: Conducted analyses; Edenberg: Phenotypic and bioinformatics expertise; Foroud: Phenotypic and bioinformatics expertise; Goate: Bioinformatics expertise; Grucza: Phenotypic expertise; Hesselbrock: Phenotypic expertise; Hinrichs: Conducted analyses, statistical/bioinformatics expertise; Kramer: Phenotypic expertise; Kuperman: Phenotypic expertise; Nurnberger: Phenotypic and statistical expertise; Porjesz: Phenotypic and statistical expertise; Saccone N: Bioinformatics expertise; Saccone S: Bioinformatics expertise; Schuckit: Phenotypic and statistical expertise; Wang: Bioinformatics expertise;

Conflict of Interest: None

Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

Reference List

  • Abecasis GR, Cherny SS, Cookson WO, Cardon LR. Merlin--rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002;30:97–101. [PubMed]
  • American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders. 3. Washington, DC: American Psychiatric Association; 1987.
  • American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders. 4. Washington, DC: American Psychiatric Association; 1994.
  • Beem AL, Geus EJ, Hottenga JJ, Sullivan PF, Willemsen G, Slagboom PE, Boomsma DI. Combined linkage and association analyses of the 124-bp allele of marker D2S2944 with anxiety, depression, neuroticism and major depression. Behav Genet. 2006;36:127–136. [PubMed]
  • Begleiter H, Reich T, Hesselbrock V, Porjesz B, Li TK, Schuckit M, Edenberg H, Rice JP. The collaborative study on the genetics of alcoholism. Alcohol Health Res World. 1995;19:228–236.
  • Bierut LJ, Dinwiddie SH, Begleiter H, Crowe RR, Hesselbrock V, Nurnberger JI, Jr, Porjesz B, Schuckit MA, Reich T. Familial transmission of substance dependence: alcohol, marijuana, cocaine, and habitual smoking: a report from the Collaborative Study on the Genetics of Alcoholism. Arch Gen Psych. 1998;55:982–988. [PubMed]
  • Bierut LJ, Madden PA, Breslau N, Johnson EO, Hatsukami D, Pomerleau OF, Swan GE, Rutter J, Bertelsen S, Fox L, Fugman D, Goate AM, Hinrichs AL, Konvicka K, Martin NG, Montgomery GW, Saccone NL, Saccone SF, Wang JC, Chase GA, Rice JP, Ballinger DG. Novel genes identified in a high-density genome wide association study for nicotine dependence. Hum Mol Genet. 2007;16:24–35. [PMC free article] [PubMed]
  • Bierut LJ, Rice JP, Goate A, Hinrichs AL, Saccone NL, Foroud T, Edenberg HJ, Cloninger CR, Begleiter H, Conneally PM, Crowe RR, Hesselbrock V, Li TK, Nurnberger JI, Jr, Porjesz B, Schuckit MA, Reich T. A genomic scan for habitual smoking in families of alcoholics: common and specific genetic factors in substance dependence. Am J Med Genet A. 2004;124:19–27. [PubMed]
  • Bucholz KK, Cadoret RJ, Cloninger RC, Dinwiddie SH, Hesselbrock V, Nurnberger JI, Reich T, Schmidt I, Schuckit MA. A New, Semi-Structured Psychiatric Interview For Use In Genetic Linkage Studies. J Stud Alcohol. 1994;55:149–158. [PubMed]
  • Centers for Disease Control CDC. Alcohol-Attributable Deaths and Years of Potential Life Lost --- United States, 2001. MMWR. 2004;53:870. [PubMed]
  • Dick DM, Bierut LJ. The genetics of alcohol dependence. Curr Psychiatry Rep. 2006;8:151–157. [PubMed]
  • Dick DM, Li TK, Edenberg HJ, Hesselbrock V, Kramer J, Kuperman S, Porjesz B, Bucholz K, Goate A, Nurnberger J, Foroud T. A genome-wide screen for genes influencing conduct disorder. Mol Psychiatry. 2004;9:81–86. [PubMed]
  • Edenberg HJ, Dick DM, Xuei X, Tian H, Almasy L, Bauer LO, Crowe RR, Goate A, Hesselbrock V, Jones K, Kwon J, Li TK, Nurnberger JI, Jr, O’Connor SJ, Reich T, Rice J, Schuckit MA, Porjesz B, Foroud T, Begleiter H. Variations in GABRA2, encoding the alpha 2 subunit of the GABAA receptor, are associated with alcohol dependence and with brain oscillations. Am J Hum Gen. 2004;74:705–714. [PubMed]
  • Edenberg HJ, Bierut LJ, Boyce P, Cao M, Cawley S, Chiles R, Doheny KF, Hansen M, Hinrichs T, Jones K, Kelleher M, Kennedy GC, Liu G, Marcus G, McBride C, Murray SS, Oliphant A, Pettengill J, Porjesz B, Pugh EW, Rice JP, Rubano T, Shannon S, Steeke R, Tischfield JA, Tsai YY, Zhang C, Begleiter H. Description of the data from the Collaborative Study on the Genetics of Alcoholism COGA and single-nucleotide polymorphism genotyping for Genetic Analysis Workshop 14. BMC Genet. 2005;6(Suppl 1):S2. [PMC free article] [PubMed]
  • Edenberg HJ, Xuei X, Chen HJ, Tian H, Wetherill LF, Dick DM, Almasy L, Bierut L, Bucholz KK, Goate A, Hesselbrock V, Kuperman S, Nurnberger J, Porjesz B, Rice J, Schuckit M, Tischfield J, Begleiter H, Foroud T. Association of alcohol dehydrogenase genes with alcohol dependence: a comprehensive analysis. Hum Mol Genet. 2006a;15:1539–1549. [PubMed]
  • Edenberg HJ, Foroud T. The genetics of alcoholism: identifying specific genes through family studies. Addict Biol. 2006b;11:386–396. [PubMed]
  • Fehr C, Rademacher BL, Buck KJ. Evaluation of the glutamate decarboxylase genes Gad1 and Gad2 as candidate genes for acute ethanol withdrawal severity in mice. Prog Neuropsychopharmacol Biol Psych. 2003;27:467–472. [PubMed]
  • Feighner JP, Robins E, Guze SB, Woodruff RA, Jr, Winokur G, Munoz R. Diagnostic criteria for use in psychiatric research. Arch Gen Psych. 1972;26:57–63. [PubMed]
  • Foroud T, Edenberg HJ, Goate A, Rice J, Flury L, Koller DL, Bierut LJ, Conneally PM, Nurnberger JI, Bucholz KK, Li TK, Hesselbrock V, Crowe R, Schuckit M, Porjesz B, Begleiter H, Reich T. Alcoholism susceptibility loci: confirmation studies in a replicate sample and further mapping. Alcohol Clin Exp Res. 2000;24:933–945. [PubMed]
  • Gelernter J, Panhuysen C, Weiss R, Brady K, Hesselbrock V, Rounsaville B, Poling J, Wilcox M, Farrer L, Kranzler HR. Genomewide linkage scan for cocaine dependence and related traits: significant linkages for a cocaine-related trait and cocaine-induced paranoia. Am J Med Genet B Neuropsychiatr Genet. 2005;136:45–52. [PubMed]
  • Gelernter J, Panhuysen C, Wilcox M, Hesselbrock V, Rounsaville B, Poling J, Weiss R, Sonne S, Zhao H, Farrer L, Kranzler HR. Genomewide linkage scan for opioid dependence and related traits. Am J Hum Genet. 2006;78:759–769. [PubMed]
  • Heath AC, Bucholz KK, Madden PA, Dinwiddie SH, Slutske WS, Bierut LJ, Statham DJ, Dunne MP, Whitfield JB, Martin NG. Genetic and environmental contributions to alcohol dependence risk in a national twin sample: consistency of findings in women and men. Psychol Med. 1997;27:1381–1396. [PubMed]
  • Helzer JE, Bucholz KK, Bierut LJ, Regier DA, Schuckit MA, Guth SE. Should DSM-V include dimensional diagnostic criteria for alcohol use disorders? Alcohol Clin Exp Res. 2006a;30:303–310. [PubMed]
  • Helzer JE, van den Brink W, Guth SE. Should there be both categorical and dimensional criteria for the substance use disorders in DSM-V? Addiction. 2006b;101:17–22. [PubMed]
  • Hesselbrock M, Easton C, Bucholz KK, Schuckit M, Hesselbrock V. A validity study of the SSAGA--a comparison with the SCAN. Addiction. 1999;94:1361–1370. [PubMed]
  • Hesselbrock V, Dick D, Hesselbrock M, Foroud T, Schuckit M, Edenberg H, Bucholz K, Kramer J, Reich T, Goate A, Bierut L, Rice JP, Nurnberger JI., Jr The search for genetic risk factors associated with suicidal behavior. Alcohol Clin Exp Res. 2004;28:70S–76S. [PubMed]
  • Hettema JM, An SS, Neale MC, Bukszar J, van den Oord EJ, Kendler KS, Chen X. Association between glutamic acid decarboxylase genes and anxiety disorders, major depression, and neuroticism. Mol Psychiatry. 2006;11:752–762. [PubMed]
  • Hinrichs AL, Bertelsen S, Bierut LJ, Dunn G, Jin CH, Kauwe JS, Suarez BK. Multipoint identity-by-descent computations for single-point polymorphism and microsatellite maps. BMC Genet. 2005;6(Suppl 1):S34. [PMC free article] [PubMed]
  • Hinrichs AL, Suarez BK. Genotyping errors, pedigree errors, and missing data. Genet Epidemiol. 2005;29(Suppl 1):S120–4. [PubMed]
  • Hopfer CJ, Lessem JM, Hartman CA, Stallings MC, Cherny SS, Corley RP, Hewitt JK, Krauter KS, Mikulich-Gilbertson SK, Rhee SH, Smolen A, Young SE, Crowley TJ. A genome-wide scan for loci influencing adolescent cannabis dependence symptoms: Evidence for linkage on chromosomes 3 and 9. Drug Alcohol Depend. 2006;89:34–41. [PMC free article] [PubMed]
  • Huang Q, Shete S, Amos CI. Ignoring linkage disequilibrium among tightly linked markers induces false-positive evidence of linkage for affected sib pair analysis. Am J Hum Genet. 2004;75:1106–1112. [PubMed]
  • Kauwe JS, Bertelsen S, Bierut LJ, Dunn G, Hinrichs AL, Jin CH, Suarez BK. The efficacy of short tandem repeat polymorphisms versus single-nucleotide polymorphisms for resolving population structure. BMC Genet. 2005;6(Suppl 1):S84. [PMC free article] [PubMed]
  • Kendler KS, Jacobson KC, Prescott CA, Neale MC. Specificity of genetic and environmental risk factors for use and abuse/dependence of cannabis, cocaine, hallucinogens, sedatives, stimulants, and opiates in male twins. Am J Psych. 2003a;160:687–695. [PubMed]
  • Kendler KS, Neale MC, Heath AC, Kessler RC, Eaves LJ. A twin-family study of alcoholism in women. Am J Psych. 1994;151:707–715. [PubMed]
  • Kendler KS, Prescott CA, Myers J, Neale MC. The structure of genetic and environmental risk factors for common psychiatric and substance use disorders in men and women. Arch Gen Psych. 2003b;60:929–937. [PubMed]
  • Kuo PH, Neale MC, Riley BP, Webb BT, Sullivan PF, Vittum J, Patterson DG, Thiselton DL, van den Oord EJ, Walsh D, Kendler KS, Prescott CA. Identification of susceptibility loci for alcohol-related traits in the Irish Affected Sib Pair Study of Alcohol Dependence. Alcohol Clin Exp Res. 2006;30:1807–1816. [PubMed]
  • Lander E, Kruglyak L. Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nat Genet. 1995;11 [PubMed]
  • Lappalainen J, Krupitsky E, Kranzler HR, Luo X, Remizov M, Pchelina S, Taraskina A, Zvartau E, Rasanen P, Makikyro T, Somberg LK, Krystal JH, Stein MB, Gelernter J. Mutation screen of the GAD2 gene and association study of alcoholism in three populations. Am J Med Genet B Neuropsychiatr Genet. 2007;144:183–192. [PubMed]
  • Li MD, Payne TJ, Ma JZ, Lou XY, Zhang D, Dupont RT, Crews KM, Somes G, Williams NJ, Elston RC. A genomewide search finds major susceptibility loci for nicotine dependence on chromosome 10 in African Americans. Am J Hum Genet. 2006;79:745–751. [PubMed]
  • Lin Z, Walther D, Yu XY, Drgon T, Uhl GR. The human serotonin receptor 2B: coding region polymorphisms and association with vulnerability to illegal drug abuse. Pharmacogenetics. 2004;14:805–811. [PubMed]
  • Lin Z, Walther D, Yu XY, Li S, Drgon T, Uhl GR. SLC18A2 promoter haplotypes and identification of a novel protective factor against alcoholism. Hum Mol Genet. 2005;14:1393–1404. [PubMed]
  • Liu QR, Drgon T, Johnson C, Walther D, Hess J, Uhl GR. Addiction molecular genetics: 639,401 SNP whole genome association identifies many “cell adhesion” genes. Am J Med Genet B Neuropsychiatr Genet. 2006;141:918–925. [PubMed]
  • Loh e, Lane HY, Chen CH, Chang PS, Ku LW, Wang KH, Cheng AT. Glutamate decarboxylase genes and alcoholism in Han Taiwanese men. Alcohol Clin Exp Res. 2006;30:1817–1823. [PubMed]
  • Lynskey MT, Heath AC, Nelson EC, Bucholz KK, Madden PA, Slutske WS, Statham DJ, Martin NG. Genetic and environmental contributions to cannabis dependence in a national young adult twin sample. Psychol Med. 2002;32:195–207. [PubMed]
  • McGue M, Elkins I, Iacono WG. Genetic and environmental influences on adolescent substance use and abuse. Am J Med Genet. 2000;96:671–677. [PubMed]
  • McPeek MS, Sun L. Statistical tests for detection of misspecified relationships by use of genome-screen data. Am J Hum Genet. 2000;66:1076–1094. [PubMed]
  • North BV, Curtis D, Sham PC. A note on calculation of empirical P values from Monte Carlo procedure. Am J Hum Genet. 2003;72:498–499. [PubMed]
  • Nurnberger JI, Jr, Foroud T, Flury L, Su J, Meyer ET, Hu K, Crowe R, Edenberg H, Goate A, Bierut L, Reich T, Schuckit M, Reich W. Evidence for a locus on chromosome 1 that influences vulnerability to alcoholism and affective disorder. Am J Psychiatry. 2001;158:718–724. [PubMed]
  • Nurnberger JI, Jr, Wiegand R, Bucholz K, O’Connor S, Meyer ET, Reich T, Rice J, Schuckit M, King L, Petti T, Bierut L, Hinrichs AL, Kuperman S, Hesselbrock V, Porjesz B. A family study of alcohol dependence: coaggregation of multiple disorders in relatives of alcohol-dependent probands. Arch Gen Psychiatry. 2004;61:1246–1256. [PubMed]
  • Philibert R, Caspers K, Langbehn D, Troughton EP, Yucuis R, Sandhu HK, Cadoret RJ. The association of the D2S2944 124 bp allele with recurrent early onset major depressive disorder in women. Am J Med Genet B Neuropsychiatr Genet. 2003;121:39–43. [PubMed]
  • Prescott CA, Kendler KS. Genetic and Environmental Influences on Alcohol and Tobacco Dependence among Women. NIAAA Res Monogr. 1995;30:59–87.
  • Prescott CA, Sullivan PF, Kuo PH, Webb BT, Vittum J, Patterson DG, Thiselton DL, Myers JM, Devitt M, Halberstadt LJ, Robinson VP, Neale MC, van den Oord EJ, Walsh D, Riley BP, Kendler KS. Genomewide linkage study in the Irish affected sib pair study of alcohol dependence: evidence for a susceptibility region for symptoms of alcohol dependence on chromosome 4. Mol Psychiatry. 2006;11:603–611. [PubMed]
  • Reich T. A genomic survey of alcohol dependence and related phenotypes: results from the Collaborative Study on the Genetics of Alcoholism COGA. Alcohol Clin Exp Res. 1996;20:133A–137A. [PubMed]
  • Reich T, Edenberg HJ, Goate A, Williams JT, Rice JP, Van EP, Foroud T, Hesselbrock V, Schuckit MA, Bucholz K, Porjesz B, Li TK, Conneally PM, Nurnberger JI, Jr, Tischfield JA, Crowe RR, Cloninger CR, Wu W, Shears S, Carr K, Crose C, Willig C, Begleiter H. Genome-wide search for genes affecting the risk for alcohol dependence. Am J Med Genet. 1998;81:207–215. [PubMed]
  • Saccone NL, Kwon JM, Corbett J, Goate A, Rochberg N, Edenberg HJ, Foroud T, Li TK, Begleiter H, Reich T, Rice JP. A genome screen of maximum number of drinks as an alcoholism phenotype. Am J Med Genet. 2000;96:632–637. [PubMed]
  • Saccone SF, Saccone NL, Neuman RJ, Rice JP. Genetic analysis of the maximum drinks phenotype. BMC Genet. 2005;6(Suppl 1):S124. [PMC free article] [PubMed]
  • Saha TD, Chou SP, Grant BF. Toward an alcohol use disorder continuum using item response theory: results from the National Epidemiologic Survey on Alcohol and Related Conditions. Psychol Med. 2006;36:931–941. [PubMed]
  • Schuckit MA, Edenberg HJ, Kalmijn J, Flury L, Smith TL, Reich T, Bierut L, Goate A, Foroud T. A genome-wide search for genes that relate to a low level of response to alcohol. Alcohol Clin Exp Res. 2001;25:323–329. [PubMed]
  • Schuckit MA, Wilhelmsen K, Smith TL, Feiler HS, Lind P, Lange LA, Kalmijn J. Autosomal linkage analysis for the level of response to alcohol. Alcohol Clin Exp Res. 2005;29:1976–1982. [PubMed]
  • Sham PC, Purcell S, Cherny SS, Abecasis GR. Powerful regression-based quantitative-trait linkage analysis of general pedigrees. Am J Hum Genet. 2002;71:238–253. [PubMed]
  • Stallings MC, Corley RP, Hewitt JK, Krauter KS, Lessem JM, Mikulich SK, Rhee SH, Smolen A, Young SE, Crowley TJ. A genome-wide search for quantitative trait loci influencing substance dependence vulnerability in adolescence. Drug Alcohol Depend. 2003;70:295–307. [PubMed]
  • Straub RE, Sullivan PF, Ma Y, Myakishev MV, Harris-Kerr C, Wormley B, Kadambi B, Sadek H, Silverman MA, Webb BT, Neale MC, Bulik CM, Joyce PR, Kendler KS. Susceptibility genes for nicotine dependence: a genome scan and followup in an independent sample suggest that regions on chromosomes 2, 4, 10, 16, 17 and 18 merit further study. Mol Psychiatry. 1999;4:129–144. [PubMed]
  • Tsuang MT, Bar JL, Harley RM, Lyons MJ. The Harvard Twin Study of Substance Abuse: what we have learned. Harv Rev Psychiatry. 2001;9:267–279. [PubMed]
  • Tsuang MT, Lyons MJ, Eisen SA, Goldberg J, True W, Lin N, Meyer JM, Toomey R, Faraone SV, Eaves L. Genetic influences on DSM-III-R drug abuse and dependence: a study of 3,372 twin pairs. Am J Med Genet. 1996;67:473–477. [PubMed]
  • Tsuang MT, Lyons MJ, Meyer JM, Doyle T, Eisen SA, Goldberg J, True W, Lin N, Toomey R, Eaves L. Co-occurrence of abuse of different drugs in men: the role of drug-specific and shared vulnerabilities. Arch Gen Psychiatry. 1998;55:967–972. [PubMed]
  • Uhl GR, Liu QR, Walther D, Hess J, Naiman D. Polysubstance abuse-vulnerability genes: genome scans for association, using 1,004 subjects and 1,494 single-nucleotide polymorphisms. Am J Hum Genet. 2001;69:1290–1300. [PubMed]
  • van den Bree MB, Johnson EO, Neale MC, Pickens RW. Genetic and environmental influences on drug use and abuse/dependence in male and female twins. Drug Alcohol Depend. 1998;52:231–241. [PubMed]
  • Vink JM, Beem AL, Posthuma D, Neale MC, Willemsen G, Kendler KS, Slagboom PE, Boomsma DI. Linkage analysis of smoking initiation and quantity in Dutch sibling pairs. Pharmacogenomics J. 2004;4:274–282. [PubMed]
  • Wang JC, Hinrichs AL, Stock H, Budde J, Allen R, Bertelsen S, Kwon JM, Wu W, Dick DM, Rice J, Jones K, Nurnberger JI, Jr, Tischfield J, Porjesz B, Edenberg HJ, Hesselbrock V, Crowe R, Schuckit M, Begleiter H, Reich T, Goate AM, Bierut LJ. Evidence of common and specific genetic effects: association of the muscarinic acetylcholine receptor M2 CHRM2 gene with alcohol dependence and major depressive syndrome. Hum Mol Genet. 2004;13:1903–1911. [PubMed]
  • Wang KS, Liu M, Paterson AD. Evaluating outlier loci and their effect on the identification of pedigree errors. BMC Genet. 2005;6(Suppl 1):S155. [PMC free article] [PubMed]
  • Williams GD, Grant BF, Stinson FS, Zobeck TS, Aitken SS, Noble J. Trends in alcohol-related morbidity and mortality. Public Health Rep. 1988;103:592–597. [PMC free article] [PubMed]
  • Zubenko GS, Hughes HB, III, Stiffler JS, Zubenko WN, Kaplan BB. D2S2944 identifies a likely susceptibility locus for recurrent, early-onset, major depression in women. Mol Psychiatry. 2002;7:460–467. [PubMed]
  • Zubenko GS, Hughes HB, III, Stiffler JS, Brechbiel A, Zubenko WN, Maher BS, Marazita ML. Sequence variations in CREB1 cosegregate with depressive disorders in women. Mol Psychiatry. 2003a;8:611–618. [PubMed]
  • Zubenko GS, Maher B, Hughes HB, III, Zubenko WN, Stiffler JS, Kaplan BB, Marazita ML. Genome-wide linkage survey for genetic loci that influence the development of depressive disorders in families with recurrent, early-onset, major depression. Am J Med Genet B Neuropsychiatr Genet. 2003b;123:1–18. [PubMed]