Multipoint logarithm of odds (LOD) scores estimated from DeFries-Fulker regression analyses across the genome. A, Dependence vulnerability. B, Conduct disorder symptoms. C, Dependence vulnerability and conduct disorder symptoms.
Multipoint logarithm of odds (LOD) scores estimated from DeFries-Fulker regression analyses for chromosome 3.
Multipoint logarithm of odds (LOD) scores estimated from DeFries-Fulker regression analyses for chromosome 9.
Multipoint logarithm of odds (LOD) scores estimated from DeFries-Fulker regression analyses for chromosome 17.
Customize your JAMA Network experience by selecting one or more topics from the list below.
Stallings MC, Corley RP, Dennehey B, et al. A Genome-Wide Search for Quantitative Trait Loci That Influence Antisocial Drug Dependence in Adolescence. Arch Gen Psychiatry. 2005;62(9):1042–1051. doi:10.1001/archpsyc.62.9.1042
Among adolescents, externalizing problem behavior and substance use disorders are often comorbid. Familial influences, including shared genetic risk factors, may account for part of this comorbidity. Previously we reported 2 chromosomal regions (3q24-3q25 and 9q34) likely to contain genes that influence substance dependence vulnerability (DV) in adolescence.
To identify quantitative trait loci (QTLs) that influence externalizing problem behavior in adolescence and to determine whether any identified QTL overlap chromosomal regions that influence DV.
Regression-based QTL mapping procedures designed for selected sibling pair samples.
Patient probands were drawn from consecutive admissions to residential and outpatient (milieu-type) treatment facilities for substance abuse and deliquency operated by the University of Colorado; most of these patients were referred for treatment by juvenile justice or social service agencies.
A total of 249 proband-sibling pairs from 191 families were selected for the study. Patient probands were 13 to 19 years of age; siblings of the probands ranged in age from 12 to 25 years.
Main Outcome Measures
A community-based sample of 4493 adolescents and young adults was used to define clinically significant, heritable, age- and sex-normed indexes of DV, conduct disorder symptoms (CDS), and a composite index of antisocial substance dependence (DV + CDS). Siblings and parents were genotyped for 374 microsatellite markers distributed across the 22 autosomes (mean intermarker distance, 9.2 centimorgans).
For both DV and CDS, there was evidence of linkage to the same region on chromosome 9q34, as well as to 3q24-3q25 for DV, and a novel region on chromosome 17q12 for CDS. Our composite index (DV + CDS) yielded the strongest evidence for linkage (logarithm of odds = 2.65) to the chromosome 9q34 region.
These results provide the first evidence of a potential molecular genetic basis for the comorbidity between DV and antisocial behavior.
Among adolescents, conduct disorder (CD) and substance use disorders (SUDs) are frequently comorbid in both clinical1-7 and general population samples.8-11 A recent meta-analysis of 51 twin and adoption studies suggests that approximately 40% of the variation in adolescent and adult antisocial behavior can be accounted for by genetic factors,12 and large-scale twin studies have estimated the heritability of CD at 24% to 74%13 and 70%.14 Family, twin, and adoption studies suggest that genetic variation influences the risk of abuse and dependence on alcohol and nicotine,15-25 and studies of illicit drug abuse and dependence also suggest important genetic and familial contributions to risk for the development of illicit SUDs.26-35
Although the etiologic mechanisms underlying the co-occurrence of CD and SUDs remain unclear, support is increasing for the hypothesis that common genetic risks may explain part of their comorbidity.36-39 Moreover, this common risk may extend to a broader array of disinhibitory psychopathologic features in general.40-53
Because of the strong comorbidity and the evidence that genetic risks contribute to both CD and SUDs, we have hypothesized that common genetic influences may contribute to these disorders. Our analytical approach requires probands who are highly selected for our phenotypes of interest (ie, extreme scoring). We have specifically recruited a selected sample of adolescent patient probands and their families (including sibling pairs) from consecutive admissions to treatment programs that specialize in the treatment of adolescents with SUDs and comorbid antisocial behavioral problems. In a recent genome-wide search for quantitative trait loci (QTLs) underlying substance dependence vulnerability (DV),54 we reported regions on chromosome 3q24-3q25 and chromosome 9q34 that may confer risk for substance dependence. In the same adolescent sample we now extend this work, searching for chromosomal regions that may be involved in comorbid SUDs and CD. Specifically, the gene-mapping goals of the present study were to identify QTLs that influence adolescent CD and to determine whether there are overlapping chromosomal regions likely to contain genes that influence comorbid risk. Indication of shared chromosomal regions would provide the first evidence of a molecular genetic basis for the comorbidity between these common disorders.
Details regarding sample ascertainment and inclusion and exclusion criteria for the adolescent clinical and community samples have been described previously.54 Patient probands were drawn from consecutive admissions to residential and outpatient (milieu-type) treatment facilities for substance abuse and deliquency operated by the University of Colorado; most of these patients were referred for treatment by juvenile justice or social service agencies. The clinical sample included 249 proband-sibling pairs from 191 families. Probands were 13 to 19 years of age (mean ± SD age, 15.9 ± 1.3 years) at the time of assessment. The ethnicity distribution was 7.8% African American, 36.5% Hispanic, 52.1% white, and 3.6% other. The community sample included 4493 individuals 12 to 25 years of age (mean ± SD age, 15.9 ± 2.1 years) drawn from the Colorado Twin Registry,55 the Colorado Adoption Project,56 and the Colorado Adolescent Substance Abuse family study.57 The ethnicity distribution of the community sample was 81.6% white, 12.1% Hispanic, 2.3% African American, and 4% of mixed or unknown ethnicity. We have previously reported55 that substance use patterns and prevalence rates in the community sample are comparable to those reported in national studies.58,59
Participants were administered an extensive battery of cognitive, psychiatric, and sociodemographic assessments that included structured diagnostic interviews and self-report questionnaires in the treatment facilities (probands) or homes (siblings) as described previously.54 All participants who were at least 18 years of age gave informed written consent to participate. For juveniles, a parent or legal guardian gave written consent for participation, and assent was obtained from all juveniles. All research protocols and consent forms were approved by the institutional review boards of the University of Colorado, Denver (for proband assessment), and the University of Colorado, Boulder (for assessments of all other relatives), and all subjects were paid for participation.
Substance use data were obtained using the Composite International Diagnostic Interview-Substance Abuse Module (CIDI-SAM).54,60 Conduct disorder was assessed using the Diagnostic Interview Schedule for Children (DISC) for all patient probands. Conduct disorder for siblings younger than 19 years was assessed using the DISC, whereas siblings 19 years or older were assessed using the Diagnostic Interview Schedule (DIS). Participants recruited between 1993 and 1997 were assessed with DSM-III-R versions of the DISC61 and DIS.62 Participants recruited between 1998 and 2001 were assessed using DSM-IV versions of the DISC63 and DIS,64 modified for scoring both DSM-III-R and DSM-IV symptoms and diagnoses. Adjustments for potential instrument effects are described herein.
We used DSM criterion (symptom) counts to create 3 quantitative phenotypes to assess DV, CD symptoms (CDS), and a composite index of antisocial substance dependence (DV + CDS). The mean number of DSM-IV dependence symptoms (ie, total symptom count across all classes of substances divided by the number of substances used repeatedly) was used as our index of DV (see Stallings et al54 for details regarding the validity of this index), and the total number of DSM-III-R CD symptoms was used as our quantitative measure of CDS. (DSM-III-R criteria were used for CDS because DSM-IV versions of the DISC and DIS assessments were unavailable at the outset of the study; the version of the CIDI-SAM that was used provided DSM-IV symptom counts for substance dependence.)
Because of marked age trends on these indexes in our community sample, the measures were age and sex corrected using standard regression procedures and then standardized within sex groups to express scores in our selected (clinical) sibling sample in relation to the population means and variability for adolescents and young adults. In addition, CDS scores were obtained from the DISC for probands and younger siblings but from the DIS for older siblings, and 2 versions of the DISC and DIS assessments were used. To account for any instrument effects that may not have been corrected for by our age adjustments, interview type was dummy coded and used to obtain residual CDS scores. (Instrument type accounted for only 1.2% of the variance in our age- and sex-corrected CDS scores in the community sample.)
Standardized DV and CDS scores were summed to form a composite index of antisocial substance dependence (DV + CDS). This index was then subsequently restandardized so that the metric of the composite (sum) index was the same as for our individual measures of DV and CDS (ie, means of 0 and SDs of 1 in our community sample).
The methods used for preparation of genomic DNA and genotyping of the 374 microsatellite markers (ABI PRISM LMS2-MD10 panels; PE Biosystems, Foster City, Calif) have been described previously.54 For 47% of the families, genotypes were available from both parents; 49% included genotypes from 1 parent, and for 4% of the sibships, no parental genotypes were available. Sibships included 147 pairs, 35 trios, 7 sibships of 4, and 3 sibships of 5. Sex-averaged marker maps were obtained from the Marshfield Center for Medical Genetics database (http://research.marshfieldclinic.org/genetics), and allele frequencies were estimated from the full sample using Merlin.65Table 1 gives marker information for each chromosome and means across the 22 autosomes.
Validation of genotypes and family relationships was performed with the Graphical Relationship Representation program,66 Pedstats,67 and Merlin.65 Details are described in our previous communication.54 Genetic analysis indicated 28 individuals in whom paternity and/or full-sibling status was questioned. In addition, 1 proband did not have complete CIDI-SAM and DISC assessments, and 11 individuals were genotyped for less than 50% of the markers. These individuals were removed from the sample before linkage analysis, resulting in a final sample of 249 sibling pairs. Twelve markers were omitted for low call rates; the mean call rate for the markers used was approximately 85%. Estimated marker haplotypes from Merlin indicated no evidence for excessive recombination events, and few genotypes (0.07%) were flagged as highly unlikely. Genotypes flagged as unlikely were not concentrated among certain markers, individuals, or families, and analyses that removed unlikely genotypes yielded essentially identical results to full-sample analyses, so these genotypes were retained. Identity-by-descent (IBD) sharing for all sibpairs was estimated at 1-centimorgan (cM) intervals (and also at each marker for single-point analyses) across the genome (autosomes) using Merlin.
The regression-based sibpair linkage method that we used was proposed by Fulker et al68,69 as an extension of the methods originally developed by Fulker, DeFries, and colleagues70 for the analysis of selected twin data. The rationale of the DeFries-Fulker (DF) linkage analysis for selected sibling samples is that the siblings of selected probands should demonstrate differential regression toward the population mean on a quantitative trait commensurate with the proportion of marker alleles they share IBD with the proband, given that a QTL is linked to a particular chromosomal marker. The following regression model was fit to data from 249 proband-sibling pairs:
C = α + b1P + b 2(πi) = ε,
where C is the quantitative score of the proband’s sibling, P is the proband’s score, πi is the estimated mean proportion of alleles shared IBD for the sibling pair at location (i) in the genome, and α and ε are the regression constant and disturbance term, respectively. A t test of the significance of the b2 regression coefficient provides a test of linkage, which can be expressed as an approximate logarithm of odds (LOD) score (t2/[2 × natural log of 10]) . We used the DF regression–based QTL linkage analysis as operationalized in QMS2,71 a series of integrated SAS macros available at http://ibgwww.colorado.edu/~lessem/software/qms2.html.
The percentages of participants who met the criteria for DSM diagnoses of substance dependence (DSM-IV) and CD (DSM-III-R) are given in Table 2. Note the greater prevalence of substance dependence in the clinical probands and their siblings compared with community adolescents. There is substantial polysubstance use, with more than half of the clinical probands (56.7%) meeting criteria for 2 or more substance dependence diagnoses. Siblings of the clinical probands also show substantially higher rates of substance dependence and dependence on multiple substances than community adolescents. Rates of CD are substantially higher among our clinical probands and their siblings compared with community adolescents.
On average, our clinical probands score 2 to 3 SDs higher than community adolescents on our measures of DV, CDS, and DV + CDS. Siblings of the probands show regression to the population mean but score approximately half an SD higher than community adolescents on DV and nearly 1 SD higher for CDsx and DV + CDS. Means and SDs in the community sample are 0 and 1, respectively, because scores were standardized on that sample. The mean and SD for DV differ slightly from those values because recruitment of the community sample is ongoing and the community sample is now larger (N = 4493 vs 3676) than for our previous report.54 For consistency with that investigation, the standardization of DV was not recalculated for the larger sample.
Our measures of DV and CDsx correlate substantially (r = 0.49) in our community adolescent sample (Table 2). As expected, this correlation is attenuated in our highly selected clinical probands. Biometrical analyses of our community twin sample indicate that these measures are heritable (h2DV = 0.34; h2CDsx = 0.38), and the genetic correlation between them was estimated at Rg = 0.47.
Table 3 gives the 29 chromosomal markers (several contiguous) where single-point DF regression analyses resulted in a significant t value (P<.05). Reported P values are actual values obtained without adjustment for multiple testing. At 3 marker locations (3q24-3q25, 9q34, and 17q12), single-point LOD scores exceeded 2.0 (corresponding to t values ≥3.03). Single-point analysis includes only informative sibling pairs at a given marker location, so some caution in interpretation is advised given that the number of informative pairs varies for each analysis. Despite this limitation, an advantage of single-point results is that they do not depend on a correct marker map and can be less sensitive to genotyping errors.72
Figure 1 shows LOD score plots from multipoint DF regression analyses performed at 1-cM intervals across the genome (22 autosomes). The advantage of multipoint analyses compared with single-point evaluations is that information from all 249 sibling pairs is used for each linkage test evaluation.
Multipoint mapping results indicated 3 regions with LOD scores greater than 1.0, generally corroborating the single-point results. These findings provide support for the hypothesis that common genes may explain part of the comorbidity between DV and CDS. The strongest evidence for linkage (LOD = 2.65) was found for our composite index of antisocial substance dependence (DV + CDS), suggesting evidence for a QTL that influences both DV and CDS on chromosome 9q34 (near markers D9S1826 and D9S1838). Evidence for QTLs that primarily influence DV on chromosome 3q24-3q25 (near markers D3S1279 and D3S1614) and CDS on chromosome 17q12 (near marker D17S798) was also found. These 3 chromosomal regions are shown in greater detail in Figures 2, 3, and 4. Descriptive statistics from our DF regression analyses at these locations are given in Table 4.
Table 5 gives the sibling means plotted as a function of the estimated mean number of alleles shared IBD with the proband for our maximum LOD score locations on chromosomes 3q24-3q25, 9q34, and 17q12. To compute category means, we collapsed the estimated IBD distribution into 3 categories: 0% to 25%, more than 25% to 75%, and more than 75% alleles shared IBD at each location. At our 3q24-3q25 peak location, the sibling means for DV, CDS, and DV + CDS show approximately half an SD difference between sibpairs with an estimated IBD of less than 25% compared with sibling pairs with an estimated IBD of greater than 75%. For our 9q34 location, the mean difference among these groups for DV is of similar magnitude, whereas mean differences for CDS and DV + CDS are approximately 80% of an SD. As expected from our regression results, mean differences are lower in magnitude at our 17q12 location. The pattern of means at the 3q24-3q25 location suggests the possibility of a nonadditive QTL, whereas the pattern of means at the 9q34 location is consistent with an additive QTL.
Although the regression-based methods we used are robust to departures from normality, we performed several tests on the assumptions of our analyses. Plots of residuals by predicted values (at the 3 maximum LOD score locations) indicated no obvious heteroscedasticity problems. Consistent with expectations, approximately 5% of the absolute values of the standardized residuals exceeded 2.0 and approximately 1% exceeded 3.0 for each of these analyses. Sensitivity analyses were also conducted at these locations, systematically dropping each sibling pair from the analysis (with replacement) and recomputing regression statistics from the remaining pairs. Evaluation of regression influence statistics indicated no particularly influential cases or outliers. Nonetheless, to minimize any bias introduced by outliers, we identified the 5% most extreme individuals (2.5% in each direction), truncated their scores to the value of the most extreme cases in the remaining 95% of the sample,73 and then reconducted our regression analyses for each phenotype at the 3 maximum LOD score locations. The LOD scores for DV, CDS, and DV + CDS at our 9q34 peak location increased to 1.98, 1.84, and 2.74, respectively; LOD scores for the 3 phenotypes at our 3q24-3q25 peak location were 1.55, 0.77, and 1.57, respectively; and LOD scores at the 17q12 peak location were 0.00, 1.25, and 0.28, respectively. Comparisons between these results and those from the full sample given in Table 4 indicate that our findings are not simply artifacts that result from the undue influence of a few outliers.
We used permutation simulations to derive empirical P values for our obtained results at each of our peak locations. Empirical null distributions were simulated by randomly permuting IBD estimates with phenotypic scores for a given sibling pair (maintaining the family structures, phenotypic scores, and IBD estimates from the real data set). We then conducted our regression analyses on 5000 random permutations of our data to obtain empirical P values for the regression test statistics obtained at each of our maximum LOD score locations (ie, reported P values are the proportion of times of 5000 simulations that we observed t values greater than or equal to the maximum values we obtained at 3q24-3q25, 9q34, and 17q12). Those P values, given in Table 4, are not substantially different from those assuming a t distribution, again suggesting that our results are not unduly biased by violations of the assumptions of our analyses.
To explore these data further, we used a variant of the affected sibling pair method. We defined affected status as scoring above a z-score cutoff of 1.96 for each of our 3 phenotypic measures. Sibling pairs were considered concordant extreme-scoring pairs if they both scored above the respective z-score cutoff. When only the proband exceeded the threshold and the sibling did not, the pair was considered a discordant pair. A z score greater than 1.96 would include approximately 2.5% of the general population. By defining affection status this way, we are not suggesting that many of our sibling pairs are unaffected. This analysis was used to determine whether sibling pairs concordant for extreme scores on our quantitative phenotypes also showed increased allele sharing IBD at the locations suggested by our regression analyses.
Table 6 gives the estimated proportion of alleles shared IBD for siblings defined as concordant high scoring and discordant pairs for our peak locations and for the flanking markers proximal to each of the peak locations. At susceptibility loci, the mean proportion of alleles that are IBD is expected to be less than 50% for discordant pairs and greater than 50% in concordant pairs (the null expectation for siblings is 50%). In general, there is a trend for increased allele sharing among extreme-scoring sibling pairs for our 3q24-3q25 and 9q34 peak LOD score regions. As expected, this trend is less evident for our 17q12 location, where only our CDS phenotype yielded an LOD score greater than 1.0. A trend for decreased allele sharing among discordant sibling pairs was found for all phenotypes at our 9q34 region. This was not found for our 3q24-3q25 region, again suggesting the possibility of a nonadditive QTL in this region.
This study extends our initial findings of putative QTLs underlying DV in adolescence54 by performing additional genome scans for CDS and a composite phenotype reflecting antisocial substance dependence (DV + CDS) on the same adolescent sample. Our findings provide preliminary support for the hypothesis that common genes may influence CD and substance dependence vulnerability in adolescence. Single-point analyses indicated 29 markers with nominally significant (P<.05) test results for linkage for 1 or more of our measures. Multipoint mapping analyses corroborated 3 locations, with the strongest evidence for QTLs influencing CD and DV detected on chromosome 9q34 (near markers D9S1826 and D9S1838). A QTL more strongly associated with DV was detected on chromosome 3q24-3q25 (near markers D3S1279 and D3S1614), and a putative QTL unique to CDS was reported on chromosome 17q12 (near marker D17S798).
Our maximum LOD score in the region of 9q34 meets criteria for suggestive linkage (LOD = 2.65). The evidence for linkage to 3q24-3q25 and 17q12 is less strong but of sufficient magnitude to warrant further study of these regions. Although these locations do not meet the criteria of Lander and Kruglyak74 for significant linkage, few studies of comparable sample size will individually have sufficient power to identify QTLs at this level of significance. A strength of our design was the use of a sample of 249 selected sibling pairs, which provides roughly the same power as 1000 to 2000 unselected sibling pairs,69 using our analytical methods. Even so, the power to detect QTLs of low effect size (<15% heritability) is limited. Clearly, these findings require replication in independent adolescent samples. However, a series of checks on our current analyses indicate that our linkage results are robust to outliers and the phenotypic distributional properties of our sample.
Other investigators have suggested QTLs in the vicinity of ours for substance-related behaviors. Uhl et al,75 investigating QTLs for polysubstance abuse vulnerability in an adult sample, reported a significant association with a region on chromosome 3 within 10 cM of our peak location (173-175 cM). Long et al76 also reported suggestive evidence for linkage to 3 markers on chromosome 3, spanning the region 165.3 to 181.1 cM in a genome scan for genetic linkage to alcohol dependence. Bergen et al77 reported a similar region associated with smoking status on chromosome 3 (171.7 cM) and a region close to our peak location on chromosome 9 (161 cM). A recent genome scan for QTLs associated with CD78 did not find the same locations as we did for our CDS measure. The authors did report a QTL on chromosome 3 (134 cM) near a secondary peak (LOD = 0.80) evident from our genome scan for our DV measure (Figure 2). However, their adult sample was selected for high-density alcoholism, and CD symptoms were retrospectively reported. In addition, their sample was not highly selected for CD (approximately 13% would meet criteria for a diagnosis), so our samples are not directly comparable.
There are several limitations of our study. First, endophenotypes that uniquely represent specific genetic contributions to CDS and DV are unknown. Our phenotypes were constructed from standard psychiatric symptom clusters (DSM criteria). These clusters, although useful for guiding currently available treatments, generally do not consider differences in etiologies.
Second, because CD and SUDs develop over time, we used age- and sex-normed quantitative phenotypes to express CD and DV in our selected sibling sample in relation to the population means and variability for community adolescents. Although CD and polysubstance use are characteristic of our high-risk probands and their siblings, not all of our study participants have fully passed through the age of risk, and this may limit the validity of our phenotypes.
Finally, our clinical adolescent sample was selected for both problem substance use and CD; thus, our findings may not generalize to broader substance-abusing populations without comorbid antisocial behavior. However, joint selection for both of these phenotypes would not be expected to lead to the finding of a common QTL in the absence of pleiotropy or strong linkage disequilibrium.
To explore our results further, we conducted a search for candidate genes within ±5 megabases (Mb) of each of our peak locations79 (http://genome.ucsc.edu). These searches yielded approximately 100 genes within each region. We highlight herein some candidates of interest (a full listing of identified genes in these regions can be found on our Web site at http://ibgwww.colorado.edu/cadd/a_center/research). Potential candidates in the region of 3q24-3q25 include the phospholipase D1 gene (PLD1); a serine proteinase inhibitor gene (SERPINI1); claudin 11, an oligodendrocyte transmembrane protein gene (CLDN11); and neuroligin 1 (NLGN1). A notable candidate near our 9q34 region is the dopamine β-hydroxylase gene (DBH). However, when we added a SNP marker within the DBH gene (DBH-1021C→T80) to our chromosome 9 map, there was no evidence of linkage near this marker, suggesting that our putative QTL is distal to the DBH locus. Other potential candidates in the region include GRIN1 (NMDAR1), a gene that codes for a subunit of the N-methyl-D-aspartate receptor; a neuroblastoma protein gene (AMY); a chloride intracellular channel gene (CLIC3); and a calcium channel gene (CACNA1B). Candidates within ±5 Mb of our 17q12 peak location include the serotonin transporter gene (SLC6A4); a cyclin-dependent kinase 5 regulatory subunit (CDK5R1) that codes for a protein that is a neuron-specific activator of CDK5; and OMG, a gene that codes for an oligodendrocyte myelin glycoprotein. Although we restricted our systematic search to ±5 Mb, we note with interest that DARPP-32,81 a dopamine- and cyclic adenosine monophosphate–regulated phosphoprotein that appears to be essential for dopaminergic and serotonergic neurotransmission, is located 6.3 Mb from our chromosome 17 peak. Putative gene function was ascertained from the GeneCards82 database (http://bioinfo1.weizmann.ac.il/genecards/).
In summary, CD and SUDs are frequently comorbid in adolescence.1-3,6-11 Although the precise etiologic mechanisms that underlie the frequent co-occurrence of these problem behaviors remain unclear, the current investigation was motivated by increasing evidence that common genetic risks may explain part of the comorbidity.36,37,39-51 The goal of our current program of research is to extend this biometrical evidence to localize chromosomal regions that contain genes that influence both these disorders. Results of this study provide the first evidence of a potential molecular genetic basis for the comorbidity between drug dependence vulnerability and antisocial behavior in adolescents.
Correspondence: Michael C. Stallings, PhD, Institute for Behavioral Genetics, Campus Box 447, University of Colorado, Boulder, CO 80309-0444 (Michael.Stallings@Colorado.edu).
Submitted for Publication: December 12, 2003; final revision received March 25, 2005; accepted March 28, 2005.
Funding/Support: This research was supported in part by grants DA-05131, DA-11015, and DA-12845 from the National Institute on Drug Abuse, Bethesda, Md.
Acknowledgment: This article is dedicated to the late David W. Fulker. We acknowledge Stacey Cherny, PhD, and Marissa Ehringer, PhD, for their statistical and bioinformatics consultation and Kristen Ehlers, BA, for data management. Genotyping was performed by Jane Hutchinson, MS, Liping Liu, MS, and Chinatsu McGeary, BA. We especially thank the participants of the clinical, family, twin, and adoption studies of the Center for the Genetics and Treatment of Antisocial Drug Dependence.