Email updates

Keep up to date with the latest news and content from Behavioral and Brain Functions and BioMed Central.

Open Access Research

The μ-opioid receptor gene and smoking initiation and nicotine dependence

Lan Zhang12, Kenneth S Kendler1 and Xiangning Chen1*

Author Affiliations

1 Virginia Institute of Psychiatric and Behavioral Genetics, Department of Psychiatry, Virginia Commonwealth University, Richmond, VA 23298, USA

2 Psychiatry Department, West China Hospital, Sichuan University, Chengdu, China

For all author emails, please log on.

Behavioral and Brain Functions 2006, 2:28  doi:10.1186/1744-9081-2-28

The electronic version of this article is the complete one and can be found online at: http://www.behavioralandbrainfunctions.com/content/2/1/28


Received:3 July 2006
Accepted:4 August 2006
Published:4 August 2006

© 2006 Zhang et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The gene encoding the mu-opioid receptor (OPRM1) is reported to be associated with a range of substance dependence. Experiments in knockout mice indicate that the mu-opioid receptor may mediate reinforcing effects of nicotine. In humans, opioid antagonist naltrexone may reduce the reinforcing effects of tobacco smoking. Additionally, the OPRM1 gene is located in a region showing linkage to nicotine dependence. The OPRM1 is thus a plausible candidate gene for smoking behavior. To investigate whether OPRM1 contributes to the susceptibility of smoking initiation and nicotine dependence, we genotyped 11 SNPs in the gene for 688 Caucasian subjects of lifetime smokers and nonsmokers. Three SNPs showed nominal significance for smoking initiation and one reached significance for nicotine dependence. The global test for three-marker (rs9479757-rs2075572-rs10485057) haplotypes was significant for smoking initiation (p = 0.0022). The same three-marker haplotype test was marginal (p = 0.0514) for nicotine dependence. These results suggest that OPRM1 may be involved in smoking initiation and nicotine dependence.

Background

Tobacco use is a leading cause of preventable diseases in the world and causes nearly 5 million tobacco-related deaths annually [1]. According to the World Health Organization, 1.3 billion people are smokers and half of the smokers in the world today will die of smoking related diseases [2]. Nicotine is the primary factor responsible for the addictive behaviors. In last two decades, family, twin and adoption studies have implicated that genetic factors strongly influence the behaviors of tobacco use and nicotine dependence [3,4]. In the review by Sullivan and Kendler, they estimated that genetic factors accounted for 56% of the variance in the liability of smoking initiation (SI) and 67% of the variance for progression to nicotine dependence (ND) [4].

The μ-opioid receptor encoded by the OPRM1 gene is a key factor contributing to drug addiction [5]. The μ-opioid receptor is a major site of action for endogenous opioid peptides and exogenous opioid drugs. More interestingly, some non-opioid substances including nicotine that have other primary sites of action are likely to induce the release of endogenous opioid peptides, and subsequently activate the μ receptor [5]. For instance, experiments have shown that nicotine causes a release of endogenous opioids in the brains of rat and mouse [6]. Recent studies using inbred and knockout mice strongly suggest that the μ-opioid receptor mediates both the positive and negative reinforcing effects of nicotine. Nicotine-induced antinociception, reward effects and dependence are substantially attenuated in mice lacking the μ-opioid receptor [7,8]. This implies that some aspects of the rewarding valence of nicotine require μ opioid receptors. In humans, opioid receptor antagonist naloxone may reduce the relative reinforcing effects of nicotine, thus it has been used as a smoking cessation drug [9-11]. Lerman and colleagues show that a variant in OPRM1 may predict the treatment responses to clinical nicotine replacement therapy. Smokers who have the Asp40 variant in the OPRM1 gene are likely to have a favorable response to the treatment [12]. Interestingly, the OPRM1 gene is located at 6q24-25, about 4 million basepairs from a suggestive linkage peak of nicotine dependence [13]. All these suggest that the OPRM1 is an attractive candidate gene that may influence the smoking behavior in humans.

More than 100 SNPs in the human OPRM1 gene have been identified, some of which have been tested for associations with a range of substance dependence [14-19]. Most studies focus on analyzing the variants in coding regions or 5' UTRs. The Asn40Asp variant or rs1799971 is one of the most studied polymorphisms in the gene. However, the results are not consistent. Of the studies conducted for the gene, many do not have comprehensive coverage. In this study, we use haplotype-tagging approach to select SNPs and carry out association analyses for both individual markers and haplotypes. Our results indicate that a major haplotype may be involved in tobacco smoking.

Methods and materials

Subjects

The sample was drawn from two large population-based twin studies of the Virginia Twin Registry. The sampling and ascertainment procedures were described elsewhere [20,21]. In this study, we used a subset of twins of European ancestry and randomly selected one twin from each pair. All the subjects were unrelated. All individuals were assessed with basic smoking history, the Fagerstrom Tolerance Questionnaire (FTQ) [22] and withdrawal symptoms. The FTQ was an eight-item questionnaire (score range 0–11) widely used to evaluate the severity of nicotine dependence. In this study, the non-smokers (NS) were defined as those who never smoked a cigarette up to the time of the assessment. Regular smokers with low ND (Low-ND) were defined as those who smoked at least 5 cigarettes per week for five years and their FTQ scores were between 0–2 at their lifetime maximum tobacco consumption. Regular smokers with high ND (High-ND) were those who smoked for five years or more and had an FTQ score between 7–11. In order to evaluate separately the influence of OPRM1 on SI and ND, two measurements with overlapping but not identical genetic effects [21,23], we used this 3-group design. The selection of subjects was based on this design. To estimate the influence of OPRM1 on SI, we compared the allele and genotype frequencies between the NS and regular smokers (which included both Low- and High-ND subjects). For that of ND, we compared the frequencies between the Low- and High-ND subjects. Smokers with FTQ scores between 2 and 7 were not used in this dichotomized design. Of the 688 subjects used in this study, 244 were NS, 215 were Low-ND and 229 were High-ND.

The buccal epithelial cell samples were collected using standard cytology brushes (Fisher Scientific, Fair Lawn, New Jersey). DNA was isolated from the brushes as reported previously [24]. All individuals provided informed consent for participation in this study. Sixteen unlinked microsatellite markers were genotyped to assess the potential population stratification and no evidence of stratification was found [24].

SNPs selection and genotyping

Eleven SNPs were chosen from the dbSNP and Celera genomics database with the following criteria: (1) SNPs located in coding or regulatory regions of the gene or reported associated with drug dependence in literature; (2) haplotype-tagged SNPs suggested by HapMap database [25] or SNPbrowser™ Software 3.0 [26] using the default parameters; and (3) SNPs with minor-allele frequency ≥ 0.10 in the dbSNP database. SNP rs1799971 is a non-synonymous variant (Asn40Asp) located in exon 1. Other SNPs are located in intronic regions of OPRM1. Table 1 listed the characteristics of the selected SNPs.

Table 1. Marker characteristics

Genotyping was performed with the TaqMan genotyping method [27]. Briefly, the PCRs were conducted with 384-well microplates. To ensure the quality of genotyping, negative control samples were included in each plate. The PCRs were performed with 5 ng of genomic DNA, 0.25 μl of TaqMan assay mix (Applied Biosystems, Inc., Foster City, CA, USA) and 2.5 μl of TaqMan universal PCR master mix in a total reaction volume of 5 μl. After activating the polymerase and denaturizing DNA by heating at 95°C for 10 minutes, 40 cycles of 92°C for 15 seconds and 55°C for 1 minute were performed. After the reaction, the fluorescence intensities of reporter 1 and 2 (reporter 1:VIC, excitation = 520 ± 10 nm, emission = 550 ± 10 nm; reporter 2: FAM, excitation = 490 ± 10 nm, emission = 510 ± 10 nm) were measured by the Analyst fluorescence plate reader (LJL Biosytems, Sunnyvale, CA). Based on the ratio of fluorescence intensities, genotypes were scored by a Euclidean clustering algorithm developed in our laboratory [28].

Statistical analysis

To distinguish the effects of the OPRM1 gene on SI and ND, we made two dichotomous comparisons. For the effect on SI, we compared the genotype and allele frequencies of NS with that of regular smokers including both the Low-ND and High-ND groups. For the effect on ND, we compared the genotype and allele frequencies of the Low-ND group with that of the High-ND group. For individual SNP association analyses, genotype and allele frequencies in the three groups were compared by χ2 test or Fisher's exact test when the expected counts were less than 5 in any cell of the contingency table. Pairwise linkage disequilibrium (LD) was estimated for all subjects by the Haploview 3.2 software [29]. Multi-marker haplotype analyses were conducted with the COCAPHASE module of the UNPHASED program [30].

Results

In this study, we genotyped 11 SNPs in the OPRM1 gene. Nine of the 11 makers were in the Hardy-Weinberg Equilibrium (HWE). Two SNPs (rs524731 and rs10485057) showed significant deviation from the HWE, Table 1. When HWE was examined separately for the 3 subject groups, none of the 3 groups was significantly deviated from HWE for rs524731, the P values were 0.182, 0.616 and 0.079 for the NS, Low-ND and High-ND groups respectively. For rs10485057, the HWE deviation was caused largely by the High-ND group, the P values for the NS, Low-ND and High-ND were 0.069, 0.979 and 0.000 respectively. The average scoring rate for the 11 SNPs was 95.6% (94.2–98.7%).

Table 2 summarized the results of allelic and genotypic association analyses. Three SNPs (rs2075572, rs10485057 and rs10485058) showed significant differences (P = 0.036, 0.012 and 0.033, respectively) in genotype frequencies between the non-smokers and regular smokers. One SNP (rs2075572) showed significant difference in allele frequency between the non-smokers and regular smokers (P = 0.016). For ND, only one SNP (rs10485057) reached nominal significance (P = 0.0297). The non-synonymous marker, rs1799971, was not significant for SI or ND.

Table 2. Single marker association analyses

Pairwise LD estimates were listed in Table 3, with D' listed below the diagonal and r2 above the diagonal. These analyses indicated that there was only one LD block, covering markers 6–11. From markers 1 to 4, most pairwise LDs were low. However, several markers, i.e., 1, 3 and 5, or rs1799971, rs524731 and rs9479757, were in high LD with markers 6–11 individually. Based on this LD structure, we conducted multi-marker haplotype analyses using only the markers showing substantial pairwise LDs. The markers included in these analyses were markers 1, 3, 5, and 6 to 11. Of these analyses, several combinations reached or approached nominal significance. The most significant combination was combination 5-6-7 or rs9479757-rs2075527-rs10485057, with global P values of 0.0021 and 0.0514 for SI and ND respectively, Table 4. In this combination, there were two haplotypes contributing to elevated risks to tobacco smoking: Haplotype 1-2-1, or G-C-A, had a frequency of 0.363 in the regular smokers and that of 0.307 in the non-smokers (P = 0.0373), and haplotype 1-2-2, or G-C-G, was observed only in the regular smokers (frequency = 0.013, P = 0.0023). For ND, only the minor haplotype 1-2-2 was significant (the frequency in the High-ND group was 0.022, that in the Low-ND group was 0.002, P = 0.010). By comparing the patterns of association, it was clear that haplotype 1-2-2 associated with both SI and ND, while haplotype 1-2-1 was associated only with SI. The most abundant haplotype, 1-1-1 for combination 5-6-7 was significantly underrepresented in the regular smokers (P = 0.0167). Combination 1-3-5-6-7 produced similar results as that of combination 5-6-7. The P values for the global and individual haplotype tests were all comparable. These results indicated that the major risk haplotype defined by markers 5–7 extended to marker 1 in exon 1, the non-synonymous Asn40Asp polymorphism. This risk haplotype carried the A (Asn) allele.

Table 3. Pairwise LD of the typed SNPs*

Table 4. Multi-marker haplotype analyses

Discussion

The evidence from animal experiments, treatment response and linkage study suggests that OPRM1 may be a gene contributing to the liability of tobacco smoking and nicotine dependence. In this study, we conduct association analyses to test whether the OPRM1 gene is associated with SI and ND. We use a 3-group design to compare the genotypic and allelic distribution of 11 SNPs between the non-smokers and regular smokers and that between the Low- and High-ND smokers. We find that three SNPs show genotypic associations with SI and only one SNP is associated with ND. In haplotype analyses, a core region covered by 3-marker combination 5-6-7 (rs9479757-rs2075572-rs10485057) reaches global significance for smoking initiation. A major haplotype, 1-2-1 or G-C-A, is overrepresented in the regular smokers as compared with the non-smokers. A minor haplotype for the same marker combination, 1-2-2 or G-C-G, is significantly associated with both SI and ND. Since the frequency of the minor haplotype is relatively low (0.013–0.026), there is a chance that this haplotype may be an artifact of the haplotype reconstruction program. It is well known that haplotypes with low frequencies have much higher error rate. On the other hand, the association of the major haplotype 1-2-1 is likely true. In fact, this haplotype remains significant when we restrict the analyses to those haplotypes observed at least once in our data (frequencies were 0.370 and 0.311 for the smokers and non-smokers, and P = 0.04084). This haplotype seems to extend to the first exon, carrying the A allele (Asn) of the Asp40Asn polymorphism.

Several studies have shown that the μ-opioid receptor is involved in the reinforcing effects of drugs. The most extensively studied variation in OPRM1 is Asn40Asp (rs1799171), but the findings are not consistent [31,32]. Asn40Asp is an A/G transition at the 118th nucleotide of the coding sequence, causing an amino acid change at position 40 from asparagine (Asn) to aspartate (Asp). This change leads to the loss of a putative N-glycosylation site in the extracellular region. The Asp40 allele is associated with higher β-endorphin affinity, lower blood cortisol levels and higher aggressive threat scores as compared to the Asn40 allele. A recent meta-analysis including 28 distinct samples and over 8000 subjects concludes that this polymorphism does not increase the risk of substance dependence [33]. In our study, this polymorphism itself is not significant. However, we do notice that this polymorphism is in high LD with the core markers defining the associations with both SI and ND. Marker combination 1-3-5-6-7 produces very similar p values as that observed with the core markers (see Table 4, comparing combinations 1-3-5-6-7 and 5-6-7). These results suggest that the Asn40Asp mutation is not itself causative. The associations observed in some studies are likely a reflection of its high LD with other causative variant(s) not yet identified. In a recent study of alcohol dependence, this Asn40Asp polymorphism is not associated the phenotypes [19]. Similar to our study, while this polymorphism is partitioned in a different LD block, it does share substantial LD with a risk haplotype found in a separate LD block. Since several markers (rs1799971, rs609148 and rs648893) in the alcohol dependence study are also typed in our study, this allows us to compare the 3 marker haplotype directly. We find that the extended risk haplotype for SI in our sample is 1-1-1-2-1-2-2-2 or A-C-G-C-A-A-C-T for markers 1-3-5-6-7-8-9-10 or rs1799971-rs524731-rs9479757-rs2075572-rs10485057-rs9322447-rs609148-rs648893. Extracting markers 1-9-10 from this extended haplotype, we have 1-2-2 or A-C-T for rs1799971- rs609148-rs648893. This overlaps with one of the risk haplotypes, A-A-C-C-T for markers rs1799971-rs3823010-rs495491-rs609148-rs648893, identified in the alcoholism study. Given the high comorbidity between alcohol drinking and tobacco smoking, this finding is interesting. Because the 5' end of the gene is not in high LD for those markers we typed in our study, we cannot be certain whether the same risk haplotype underlies the risks for both alcoholism and tobacco smoking. Further studies with more markers at the 5' end of the gene are necessary to clarify this issue.

There are several factors contributing to false positives in case control studies, including population stratification, sample size and genotyping error. In our study, we have taken measurements to reduce these risks. We genotyped unlinked microsatellite markers to assess potential stratification and used reasonable sample sizes. We examined HWE for all SNPs. For those two SNPs showing departure from HWE, we checked each subject groups separately. Only one (rs10485057) of these two SNPs involves in the core risk haplotype. For this marker, HWE deviation is caused largely, if not exclusively, by the High-ND smokers. This case-specific HWE deviation, by itself, can be interpreted as an association [34,35]. Multiple testing is another source of false positives. In this study, none of the tests would survive Bonferroni correction.

Since many of these tests are correlated (many SNPs are in substantial LD and there is overlap of markers in multi-marker analyses), Bonferroni correction seems excessive in this case. Based on this rationale, we decide to use permutation tests to obtain empirical p values for the multi-marker haplotype tests. With 10,000 permutations, combinations 5-6-7, 5-6-7-8 and 1-3-5-6-7 remain significant for SI (p = 0.0090, 0.0114 and 0.0259 respectively). Based on these results, we may reasonably conclude that OPRM1 is associated with smoking initiation. Whether OPRM11 is associated with nicotine dependence is less clear. Further studies with independent samples are necessary to resolve these issues.

Acknowledgements

This study was supported by Virginia Tobacco Settlement Foundation through the Virginia Youth Tobacco Project to Virginia Commonwealth University. We thank Dr. Linda Corey for assistance with the ascertainment of twins from the Virginia Twin Registry, now part of the Mid-Atlantic Twin Registry (MATR), currently directed by Dr. Judy Silberg. The MATR has received support from the National Institutes of Health, the Carman Trust and the WM Keck, John Templeton and Robert Wood Johnson Foundations.

References

  1. Graul AI, Prous JR: Executive summary: nicotine addiction.

    Timely Top Med Cardiovasc Dis 2005, 9:E24. PubMed Abstract | Publisher Full Text OpenURL

  2. Perkins KA, Benowitz N, Henningfield J, Newhouse P, Pomerleau O, Swan G: Society for Research on Nicotine and Tobacco.

    Addiction 1996, 91:129-144. PubMed Abstract | Publisher Full Text OpenURL

  3. Batra V, Patkar AA, Berrettini WH, Weinstein SP, Leone FT: The genetic determinants of smoking.

    Chest 2003, 123:1730-1739. PubMed Abstract | Publisher Full Text OpenURL

  4. Sullivan PF, Kendler KS: The genetic epidemiology of smoking.

    Nicotine Tob Res 1999, 1(Suppl 2):S51-S57. PubMed Abstract | Publisher Full Text OpenURL

  5. Contet C, Kieffer BL, Befort K: Mu opioid receptor: a gateway to drug addiction.

    Curr Opin Neurobiol 2004, 14:370-378. PubMed Abstract | Publisher Full Text OpenURL

  6. Pomerleau OF: Endogenous opioids and smoking: a review of progress and problems.

    Psychoneuroendocrinology 1998, 23:115-130. PubMed Abstract | Publisher Full Text OpenURL

  7. Berrendero F, Kieffer BL, Maldonado R: Attenuation of nicotine-induced antinociception, rewarding effects, and dependence in mu-opioid receptor knock-out mice.

    J Neurosci 2002, 22:10935-10940. PubMed Abstract | Publisher Full Text OpenURL

  8. Kieffer BL, Gaveriaux-Ruff C: Exploring the opioid system by gene knockout.

    Prog Neurobiol 2002, 66:285-306. PubMed Abstract | Publisher Full Text OpenURL

  9. Byars JA, Frost-Pineda K, Jacobs WS, Gold MS: Naltrexone augments the effects of nicotine replacement therapy in female smokers.

    J Addict Dis 2005, 24:49-60. PubMed Abstract | Publisher Full Text OpenURL

  10. Lerman C, Patterson F, Berrettini W: Treating tobacco dependence: state of the science and new directions.

    J Clin Oncol 2005, 23:311-323. PubMed Abstract | Publisher Full Text OpenURL

  11. Covey LS, Glassman AH, Stetner F: Naltrexone effects on short-term and long-term smoking cessation.

    J Addict Dis 1999, 18:31-40. PubMed Abstract OpenURL

  12. Lerman C, Wileyto EP, Patterson F, Rukstalis M, Audrain-McGovern J, Restine S, Shields PG, Kaufmann V, Redden D, Benowitz N, et al.: The functional mu opioid receptor (OPRM1) Asn40Asp variant predicts short-term response to nicotine replacement therapy in a clinical trial.

    Pharmacogenomics J 2004. OpenURL

  13. Straub RE, Sullivan PF, Ma Y, Myakishev MV, Harris-Kerr C, Wormley B, Kadambi B, Sadek H, Silverman MA, Webb BT, et al.: 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 Abstract | Publisher Full Text OpenURL

  14. Crowley JJ, Oslin DW, Patkar AA, Gottheil E, DeMaria PA Jr, O'Brien CP, Berrettini WH, Grice DE: A genetic association study of the mu opioid receptor and severe opioid dependence.

    Psychiatr Genet 2003, 13:169-173. PubMed Abstract | Publisher Full Text OpenURL

  15. Luo X, Kranzler HR, Zhao H, Gelernter J: Haplotypes at the OPRM1 locus are associated with susceptibility to substance dependence in European-Americans.

    Am J Med Genet B Neuropsychiatr Genet 2003, 120:97-108. PubMed Abstract | Publisher Full Text OpenURL

  16. Schinka JA, Town T, Abdullah L, Crawford FC, Ordorica PI, Francis E, Hughes P, Graves AB, Mortimer JA, Mullan M: A functional polymorphism within the mu-opioid receptor gene and risk for abuse of alcohol and other substances.

    Mol Psychiatry 2002, 7:224-228. PubMed Abstract | Publisher Full Text OpenURL

  17. Gelernter J, Kranzler H, Cubells J: Genetics of two mu opioid receptor gene (OPRM1) exon I polymorphisms: population studies, and allele frequencies in alcohol- and drug-dependent subjects.

    Mol Psychiatry 1999, 4:476-483. PubMed Abstract | Publisher Full Text OpenURL

  18. Kranzler HR, Gelernter J, O'Malley S, Hernandez-Avila CA, Kaufman D: Association of alcohol or other drug dependence with alleles of the mu opioid receptor gene (OPRM1).

    Alcohol Clin Exp Res 1998, 22:1359-1362. PubMed Abstract | Publisher Full Text OpenURL

  19. Zhang H, Luo X, Kranzler HR, Lappalainen J, Yang BZ, Krupitsky E, Zvartau E, Gelernter J: Association between two {micro}-opioid receptor gene (OPRM1) haplotype blocks and drug or alcohol dependence.

    Hum Mol Genet 2006, 15:807-819. PubMed Abstract | Publisher Full Text OpenURL

  20. Prescott CA, Kendler KS: Genetic and environmental contributions to alcohol abuse and dependence in a population-based sample of male twins.

    Am J Psychiatry 1999, 156:34-40. PubMed Abstract | Publisher Full Text OpenURL

  21. Kendler KS, Neale MC, Sullivan P, Corey LA, Gardner CO, Prescott CA: A population-based twin study in women of smoking initiation and nicotine dependence.

    Psychol Med 1999, 29:299-308. PubMed Abstract | Publisher Full Text OpenURL

  22. Fagerstrom KO: Measuring degree of physical dependence to tobacco smoking with reference to individualization of treatment.

    Addict Behav 1978, 3:235-241. PubMed Abstract | Publisher Full Text OpenURL

  23. Maes HH, Sullivan PF, Bulik CM, Neale MC, Prescott CA, Eaves LJ, Kendler KS: A twin study of genetic and environmental influences n tobacco initiation, regular tobacco use and nicotine dependence.

    Psychol Med 2004, 34:1251-1261. PubMed Abstract | Publisher Full Text OpenURL

  24. Silverman MA, Neale MC, Sullivan PF, Harris-Kerr C, Wormley B, Sadek H, Ma Y, Kendler KS, Straub RE: Haplotypes of four novel single nucleotide polymorphisms in the nicotinic acetylcholine receptor beta2-subunit (CHRNB2) gene show no association with smoking initiation or nicotine dependence.

    Am J Med Genet 2000, 96:646-653. PubMed Abstract | Publisher Full Text OpenURL

  25. HapMap Database [Http://www.Hapmap.org] webcite

  26. SNPbrowser software [Http://www.allsnps.com] webcite

  27. Livak KJ: Allelic discrimination using fluorogenic probes and the 5' nuclease assay.

    Genet Anal 1999, 14:143-149. PubMed Abstract OpenURL

  28. van den Oord EJ, Jiang Y, Riley BP, Kendler KS, Chen X: FP-TDI SNP scoring by manual and statistical procedures: a study of error rates and types.

    Biotechniques 2003, 34:610-20. PubMed Abstract OpenURL

  29. Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps.

    Bioinformatics 2005, 21:263-265. PubMed Abstract | Publisher Full Text OpenURL

  30. Dudbridge F: Pedigree disequilibrium tests for multilocus haplotypes.

    Genet Epidemiol 2003, 25:115-121. PubMed Abstract | Publisher Full Text OpenURL

  31. Bond C, LaForge KS, Tian M, Melia D, Zhang S, Borg L, Gong J, Schluger J, Strong JA, Leal SM, et al.: Single-nucleotide polymorphism in the human mu opioid receptor gene alters beta-endorphin binding and activity: possible implications for opiate addiction.

    Proc Natl Acad Sci USA 1998, 95:9608-9613. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Ikeda K, Ide S, Han W, Hayashida M, Uhl GR, Sora I: How individual sensitivity to opiates can be predicted by gene analyses.

    Trends Pharmacol Sci 2005, 26:311-317. PubMed Abstract | Publisher Full Text OpenURL

  33. Arias A, Feinn R, Kranzler HR: Association of an Asn40Asp (A118G) polymorphism in the mu-opioid receptor gene with substance dependence: A meta-analysis.

    Drug Alcohol Depend 2005. OpenURL

  34. Hoh J, Wille A, Ott J: Trimming, weighting, and grouping SNPs in human case-control association studies.

    Genome Res 2001, 11:2115-2119. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Wittke-Thompson JK, Pluzhnikov A, Cox NJ: Rational inferences about departures from Hardy-Weinberg equilibrium.

    Am J Hum Genet 2005, 76:967-986. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL