We attempted replication of the association finding in a sample of 2231 EA male twins (9.68% substance dependent) ascertained through the Virginia Twin Study of Psychiatric and Substance Use Disorders. We used generalized linear models (nesting family ID and zygosity) to account for dependence structure induced by twin data by modeling the impact of genotype on phenotype while accounting for degree of relationship (monozygotic/dizygotic).