Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2012; 7(9): e45300.
Published online 2012 September 17. doi:  10.1371/journal.pone.0045300
PMCID: PMC3444488

Identification of POMC Exonic Variants Associated with Substance Dependence and Body Mass Index

Uwe Rudolph, Editor



Risk of substance dependence (SD) and obesity has been linked to the function of melanocortin peptides encoded by the proopiomelanocortin gene (POMC).

Methods and Results

POMC exons were Sanger sequenced in 280 African Americans (AAs) and 308 European Americans (EAs). Among them, 311 (167 AAs and 114 EAs) were affected with substance (alcohol, cocaine, opioid and/or marijuana) dependence and 277 (113 AAs and164 EAs) were screened controls. We identified 23 variants, including two common polymorphisms (rs10654394 and rs1042571) and 21 rare variants; 12 of which were novel. We used logistic regression to analyze the association between the two common variants and SD or body mass index (BMI), with sex, age, and ancestry proportion as covariates. The common variant rs1042571 in the 3′UTR was significantly associated with BMI in EAs (Overweight: Padj = 0.005; Obese: Padj = 0.018; Overweight+Obese: Padj = 0.002) but not in AAs. The common variant, rs10654394, was not associated with BMI and neither common variant was associated with SD in either population. To evaluate the association between the rare variants and SD or BMI, we collapsed rare variants and tested their prevalence using Fisher’s exact test. In AAs, rare variants were nominally associated with SD overall and with specific SD traits (SD: PFET,1df = 0.026; alcohol dependence: PFET,1df = 0.027; cocaine dependence: PFET,1df = 0.007; marijuana dependence: PFET,1df = 0.050) (the P-value from cocaine dependence analysis survived Bonferroni correction). There was no such effect in EAs. Although the frequency of the rare variants did not differ significantly between the normal-weight group and the overweight or obese group in either population, certain rare exonic variants occurred only in overweight or obese subjects without SD.


These findings suggest that POMC exonic variants may influence risk for both SD and elevated BMI, in a population-specific manner. However, common and rare variants in this gene may exert different effects on these two phenotypes.


Substance dependence (SD) and obesity are two prevalent health problems. They are highly heritable, but the underlying genetic mechanisms are, for the most part, not well understood. Since both substances of abuse and food have rewarding properties, it is possible that overlapping reward pathways in the brain are involved in SD and obesity development. Deficits in neural reward responses and alterations in reward homeostasis are thought to be a common mechanism for obesity and SD [1]. A number of studies have demonstrated an inverse relationship between BMI and SD [2], [3], [4], [5], [6]. Although other studies failed to support this finding [7], [8], [9], it is still of interest to identify genetic variation contributing to SD, obesity or both. Identification of variation in genes participating in the common reward pathway for SD and obesity could help in the design of hypothesis-guided therapies to treat both conditions.

Both obesity and SD can result from dysfunction of melanocortin peptides, which are components of the hypothalamic-pituitary-adrenal (HPA) axis. Melanocortin peptides are encoded by the proopiomelanocortin (POMC) gene (POMC) on chromosome 2p23 [10], [11]. They are derived from POMC after extensive tissue-specific cleavage and include as many as 10 functionally different peptides such as adrenocorticotropic hormone (ACTH), α-, β-, and γ-melanocyte-stimulating hormone (MSH), β- and γ-lipotropin, corticotropin-like intermediate peptide (CLIP) and β-endorphin (Figure 1). Melanocortin peptides play important roles in pain [12], energy homeostasis [13], melanocyte stimulation [14] and immune modulation [15]. Of these peptides, ACTH and β-endorphin have been implicated in craving for substances of abuse. Plasma ACTH and β-endorphin levels were significantly lower in heavy drinkers than in non-drinkers [16]. Further, a significant decrease in plasma levels of ACTH and β-endorphin was shown in abstinent alcoholics [17]. The decreased level of β-endorphin in alcoholics and the increased release of β-endorphin after alcohol consumption support the theory of β-endorphin deficiency in alcoholism [18]. Melanocortin peptides and their receptors are also involved in hormonal regulation of pigmentation, weight maintenance, adrenal function and exocrine gland secretion. Animal studies elucidated a dual role of α-MSH in regulating food intake and influencing hair pigmentation [18]. Augmented ACTH and β-lipotrophin secretion was shown in patients with obesity [19]. Thus, functional alterations in melanocortin peptides due to variation in POMC may predispose to SD and obesity, as well as other traits.

Figure 1
Schematic genomic structure of POMC and locations of identified variants.

Four published studies have provided evidence that POMC variants can regulate SD risk. Xuei et al [20] reported that two single nucleotide polymorphisms (SNPs) in POMC intron 1 were associated with opioid dependence (OD) in a family-based study. Racz et al [21] found that a two-SNP haplotype in POMC was associated with alcohol dependence (AD) in females. We conducted both family- and population-based studies (3,088 subjects) and found that variants in POMC promoter and intronic regions conferred vulnerability to multiple forms of SD [22]. Recently, a genome-wide association study revealed that variants in POMC were nominally associated with AD in African Americans (AAs) [23].

POMC variation also influences traits such as obesity and hair pigmentation. Loss-of-function mutations [e.g., a homozygous mutation (C3804A) in exon 2 (now exon 3) and a compound heterozygote for two mutations (G7013T and 1 bp-deletion C7133Δ) in exon 3 (now exon 4)], resulting in ACTH and α-MSH deficiency, caused severe early-onset obesity, adrenal insufficiency and red hair pigmentation [24]. Linkage analysis confirmed the trend (maximum LOD at D2S2337 = 2.03) towards linkage between polymorphic markers around POMC and obesity [25]. Heterozygosity for any of the three non-synonymous mutations (G3834C or Ser7Thr and C3840T or Ser9Leu within POMC signal peptide, and C7406G or Arg236Gly within β-endorphin) and a 9-bp insertion/deletion polymorphism (−/AGCAGCGGC or rs10654394) were found in obese children [26]. The 9-bp-insertion allele was also associated with elevated serum leptin levels [27]. However, Echwald et al [28] did not find an association between POMC exonic variants (including the 9-bp insertion/deletion polymorphism) and early-onset obesity is a sample of 156 obese Caucasians and 380 healthy controls.

To date, POMC exons have not been sequenced in a large sample. Some unknown exonic variants may have not been identified and their association with diseases such as SD is waiting to be studied. In the present study, we sequenced all POMC exons in a relatively large case-control sample. The identified common and rare exonic variants were analyzed for the association with SD and/or BMI.


Study Subjects

Subjects were recruited from three sites in the United States: the University of Connecticut Health Center (Farmington, Connecticut), Yale University School of Medicine (APT Foundation, New Haven, Connecticut), and the Medical University of South Carolina (Charleston, South Carolina). The study protocol was approved by each local institutional review board (the Institutional Review Board of the University of Connecticut, the Yale University Human Research Protection Program, and the Institutional Review Board for Human.

Research of the Medical University of South Carolina), and written informed consent was obtained from each subject. Subjects were interviewed by trained interviewers using the Semi-Structured Assessment for Drug Dependence and Alcoholism (SSADDA) instrument [29], [30], which yielded diagnoses for lifetime substance dependence (SD) and other psychiatric traits based on the criteria of the Diagnostic and Statistical Manual of Mental Disorders, 4th edition (DSM-IV) [31].

A total of 588 subjects [280 African Americans (AAs) and 308 European Americans (EAs)] were selected for this study. Among them, 311 (167 AAs and 114 EAs) were affected with substance (alcohol, cocaine, opioid and/or marijuana) dependence (AD, CD, OD and/or MjD, respectively) and 277 (113 AAs and 164 EAs) were screened controls. Of the 167 AA cases, 131 (78.4%) had AD, 118 (70.7%) had CD, 24 (14.4%) had OD, and 68 (40.7%) had MjD. Of the 144 EA cases, 112 (98.2%) had AD, 85 (74.5%) had CD, 50 (43.9%) had OD, and 44 (38.6%) had MjD. Control subjects were not affected with AD, CD, OD and/or MjD. Major psychotic disorders such as schizophrenia and bipolar disorders were exclusion criteria for both groups. Among the 588 subjects, 51 AA cases, 40 AA controls, 15 EA cases and 41 EA controls were included in our previous study [22]. Additionally, self-reports for current height and weight were obtained, and were used to estimate BMI [i.e., weight in kg/(height in meters)2]. According to the World Health Organization (WHO) criteria, three subjects (BMI ± S.D.:  = 18.4±0.1) were classified as underweight (BMI <18.5), 195 subjects (BMI ± S.D.: 22.4±1.6) were classified as normal weight (BMI: 18.5–24.9), 212 subjects (BMI ± S.D.: 27.1±1.4) were classified as overweight (BMI: 25.0–29.9), and 178 subjects (BMI ± S.D.: 35.5±5.1) were classified as obese (BMI ≥30). The mean ± S.D. of BMI for AA cases was 29.0±6.5, which was not significantly different from that of AA controls (29.2±7.1) (T  = 0.248, P  = 0.805). The mean ± S.D. of BMI for EA cases was 28.3±5.3, which was significantly higher than that of EA controls (26.1. ±5.4) (T  = 3.613, P  = 0.0004). Detailed demographic and clinical information on the study samples is presented in Table 1.

Table 1
Characteristics of case and control subjects.


POMC spans about an 8 kb genomic region (Chromosome 2: 25,383,722–25,391,772) and contains four exons. It encodes a protein of 247 amino acids. We designed primers for polymerase chain reactions (PCRs) according to POMC DNA sequence obtained from the Ensembl database (ENST00000380794). PCR products were treated with ExoSAP-IT PCR clean-up reagents (USB Corporation, Cleveland, OH, USA). DNA sequencing was performed on Applied Biosystems 3730 capillary instruments (Applied Biosystems, Foster City, CA, USA) using reagents in the BigDye Direct Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA) at the Yale Keck core facility. DNA sequencing was first conducted in the forward direction using forward primers, and then DNA sequences potentially harboring variants were sequenced in the reverse direction using reverse primers for validation. PCR and sequencing primer sequences for each exon are presented in Supporting Information Table S1 and PCR conditions for amplifying the four exons are listed in Supporting Information Table S2. DNA sequencing raw data were first analyzed using the Sequence Scanner v1.1 (Applied Biosystems, Foster City, CA, USA). Sequencing chromatograms were then loaded into the program CodonCode Aligner v3.7.1 (CodonCode Corporation, Dedham, MA, USA) to determine homozygous and heterozygous calls.

Statistical Analysis

To verify the self-reported race, we used a Bayesian model-based clustering method implemented in the program STRUCTURE [32], [33]. We estimated the African and European ancestry proportions of individual subjects, based on the genotype data of 41 ancestry informative markers (AIMs), including 36 short tandem repeat markers and five SNPs, as described previously [34], [35]. Subjects with ancestry proportion scores ≥0.50 were grouped as African Americans (AAs), and those with ancestry proportion scores <0.50 were grouped as European Americans (EAs). These two distinct groups were highly concordant with self-reported AA and EA group membership. In this study, we analyzed the association of the identified POMC exonic variants with: (1) SD overall; (2) specific forms of SD (i.e., AD, CD, OD, or MjD); and (3) BMI. BMI was converted into four categorical traits [underweight (BMI <18.5), normal weight (BMI: 18.5–24.9), overweight (BMI: 25.0–29.9), and obese (BMI ≥30)] as described in Table 1. Three models were considered (Model 1: overweight group vs. normal weight group; Model 2: obese group vs. normal weight group; Model 3: overweight + obese group vs. normal weight group).

We first analyzed the association of the two identified POMC common variants with SD or BMI. Tests for Hardy-Weinberg equilibrium (HWE) were conducted in AA and EA control subjects. Allelic association analysis was performed using Pearson’s 2×2 contingency table Chi-square (χ2) tests. Odds ratios (ORs) and 95% confidence intervals (CIs) were estimated using χ2 tests. To adjust for multiple testing, permutation tests were performed 10,000 times to calculate empirical P values. Multivariate logistic regression analysis was used to calculate P values, with covariance adjustment for sex, age, and ancestry proportion (for binary SD traits, BMI was taken as an additional covariate; for BMI in three models, SD status was taken as an additional covariate). The above association analyses were implemented by the program PLINK v1.07 [36]. Haplotypes of the two identified common variants were constructed using program Haploview v4.2. [37], and haplotype-based association analysis was performed using PLINK. Statistical power analysis was conducted by program PS (Power and Sample size Calculations, version 3.0.43) [38].

We also analyzed the association of identified POMC rare variants with SD or BMI (in three models). Rare variants were defined as SNP markers with a minor allele frequency (MAF) less than 1%, consistent with previous studies [39], [40]. We found 20 variants with MAFs less than 1% and one variant with a MAF of 0.011. All of these 21 variants were included in the aggregate analysis using Fisher exact tests. To increase statistical power, we used a collapsing method in which alleles of all identified rare variants were summed as a single variable and then compared between cases and controls using the Fisher’s Exact Test (FET) implemented in the R package version R 2.13.1 ( The total number of sequenced chromosomes was counted based on the assumption that each participant had two copies of each chromosome. Because the GC contents of the templates were different, the four POMC exons were not amplified with the same success rate. As a result, variable numbers of samples sequenced for each rare variant caused different expected numbers of segregating sites, as shown in previous studies [41]. The harmonic mean approach was thus applied to adjust for the sample size as described by Xie et al [42]: An external file that holds a picture, illustration, etc.
Object name is pone.0045300.e001.jpg, where Ni is the sample size of the i-th variant and n is the number of variants. The number of each rare variant was adjusted as pi*N, where pi is the frequency of i-th variant.

Bioinformatics Analysis

The Transcription Element Search System (TESS, was used to query for the presence of transcriptional factors (TFs) that could potentially bind to the DNA sequence harboring variants in the 5′ untranslated region (5′ UTR). The program PolyPhen ( was used to predict the effect of missense variants on protein structure and function. It gives three predictions: benign, possibly damaging, and probably damaging. The multiple sequence alignment program ClusterW, which is incorporated in the BioEdit v7.1.3 software package (, was used to identify conserved protein sequences. To predict the function of variants in the 3′ UTR region of POMC, miRNAs putatively bound to the sequence containing 3′ UTR variants were identified by the program TargetScanHuman ( The minimum free energy (MFE) for hybridization of miRNAs to target mRNA sequences was predicted using the program RNAhybrid ( The PhyloP program, which is built into the UCSC genome browser based on multiple alignments of all 46 vertebrate species, was used to calculate the evolutionary conservation score (PhyloP score) of each variant site. The absolute values of the scores represent -log P-values under a null hypothesis of neutral evolution. Positive scores indicated that the sites were possibly conserved and negative values indicated that the sites were predicted to be fast-evolving.


Identification of POMC Exonic Variants

We identified 23 POMC exonic variants by direct sequencing, including two common variants and 12 newly discovered rare variants (Figure 1 and Table 2). The identified variants have been submitted to the GeneBank (BankIt1533561). One common variant was a 9-bp deletion/insertion polymorphism (c.560_561−/AGCAGCGGC or rs10654394) in the POMC coding region (exon 4) and another was a single nucleotide polymorphism (SNP) (c.1130 C>T or rs1042571) in POMC 3′ UTR. The frequency of the minor (or 9-bp insertion) allele of marker rs10654394 was 28.4% and 4.8% in AA and EA controls, respectively; the minor allele (T) frequency of SNP rs1042571 was 14.0% and 21.0% in AA and EA controls, respectively. The other 21 variants identified were rare [20 had minor allele frequency (MAF) <1% and one variant (c.61 A>G) had MAF of 0.011 in SD cases] (Table 2).

Table 2
Common and rare variants identified in POMC exons.

Association of Two POMC Common Variants with SD or BMI

Genotype distributions of the two common variants (rs10654394 and rs1042571) were consistent with HWE expectations in both AAs and EAs (P>0.01, data not shown). As shown in Table 3, the common variant rs1042571 in the 3′ UTR was significantly associated with BMI in EAs (Normal weight vs. Overweight: Pobs  = 0.003, Pemp  = 0.003, Padj = 0.005; Normal weight vs. Obese: Pobs  = 0.012, Pemp  = 0.013, Padj  = 0.018; Normal weight vs. Overweight + Obese: Pobs  = 0.001, Pemp  = 0.002, Padj = 0.002), but not in AAs. However, there was no association between the common variant rs10654394 and BMI in either AAs or EAs (Table 3). The Haploview program [37] was used to determine the phase of the two common SNPs. They were found to be in tight LD (AAs, D′ = 0.87; EAs, D′ = 1.00) but low correlation (AAs, R2 = 0.04; EAs, R2 = 0.02). Haplotype “9 bp Del-T”, consisting of major allele “9 bp Del” of rs10654394 and minor allele T of rs1042571, was significantly associated with BMI in EAs (Normal weight vs. Overweight: Pobs  = 0.008, Pemp  = 0.008, Padj = 0.011; Normal weight vs. Overweight + Obese: Pobs  = 0.016, Pemp  = 0.019, Padj = 0.009). Haplotype “ 9 bp Del-C”, consisting of major allele “9 bp Del” of rs10654394 and major allele C of rs1042571, was significantly associated with BMI in EAs (Normal weight vs. Overweight: Pobs  = 0.027, Pemp  = 0.033, Padj = 0.040; Normal weight vs. Overweight + Obese: Pobs  = 0.019, Pemp  = 0.022, Padj = 0.016) (Supporting Information Table S3). Nevertheless, neither the single maker nor haplotype association analyses revealed association of the two common variants with SD or any specific SD traits in either population (Supporting Information Tables S4 and S5).

Table 3
Allelic association of two common POMC variants and BMI.

Compound Effects of Rare POMC Exonic Variants on SD or BMI

Fisher’s exact tests with 21 rare variants collapsed as a single variable were performed and the results of data analysis are summarized in Table 4. In AAs, rare POMC exonic variants were significantly more frequent in SD cases overall or in cases with AD, CD or MjD specifically than in controls (SD: PFET, 1df  = 0.026; AD: PFET, 1df  = 0.027; CD: PFET, 1df  = 0.007; MjD: PFET, 1df  = 0.050). However, in EAs, the frequency of rare POMC exonic variants did not differ significantly between SD cases and controls (Table 4). Bonferroni correction was used to adjust the empirical P values obtained from the Fisher Exact test, resulting in a significance level Pcorrection <0.01 (0.05/5, i.e., correction for five times of comparisons in each sample set). Only the P value (PFET, 1df  = 0.007) obtained from the analysis of CD data in AAs survived Bonferroni correction. In addition, rare POMC exonic variants did not occur significantly more frequently in the overweight or obese group than in the normal weight group in either population (Table 4). Nevertheless, we observed that certain exonic variants (c.432 G>T, c.609 C>T, c.969 C>G and c.1187 A>G) appeared only in overweight or obese subjects unaffected with SD (Figure 2).

Table 4
Collapse association of POMC rare variants with substance dependence (SD) or body mass index (BMI).
Figure 2
Scatter plots of subjects with rare variants according to substance dependence (SD) and BMI.

Functional Prediction of POMC Exonic Variants

The predicted function of POMC exonic variants is presented in Table 5. The two rare variants (c.61 A>G and c.259 G>A) in the 5′ UTR of POMC were predicted to be located in transcription factor binding sites. Variant c.61 A>G was observed in nine heterozygous individuals. The variant allele G was predicted to be located in the binding site of transcription factor E2F, and the common allele A was potentially located in the binding site of transcription factors Sp1 and CAC-binding protein. Variant c.259 G>A was found in only one heterozygous carrier. The variant allele A was predicted to be located in the binding site of transcription factor GR, and the common allele G was potentially located in the binding site of transcription factor myogenin.

Table 5
Functional prediction of identified POMC rare variants.

In the coding region, two variants (c.343 G>A and c.432 G>T) were predicted to be nonsense mutations that potentially terminate POMC mRNA translation. Each of them was identified in only one heterozygous carrier. Variant c.867_868 6-bp (GGGCCC) caused a two-amino acid (Ala-Gly) insertion. This insertion variant is located in the β-LPH domain. Seven missense mutations (c.267 C>T, c.524 C>A, c.685 T>G, c.840 G>T, c.846 G>A, c.904 A>G, and c.969 C>G) were predicted to be damaging and one missense mutations (c.421 A>G) were predicted to be benign. Levels of evolutionary conservation of rare variants in coding regions were described at both the nucleotide and amino acid levels. Comparing the POMC coding sequence of humans with those of 45 other vertebrate species showed that most of the identified rare variants had a positive PhyloP score, indicating that they tended to be conservative rather than fast-evolving (Table 5). Additionally, alignment of POMC amino acid sequence across multiple species showed that the rare variants in coding regions were located in highly conserved regions (Supporting Information Figure S1).

In the 3′UTR of POMC, variant c.1095 T del was predicted to be located in the target site of six miRNAs. The variant deletion allele may cause an increased MFE for four miRNAs (hsa-mir-4728-5p, hsa-mir-4665-5p, hsa-mir-1275 and hsa-mir-625) but a decreased MFE for two miRNAs (hsa-mir-4723-5p and hsa-mir-4525). Variant c.1130 G>T was predicted to be located in the target site of two miRNAs (hsa-mir-3715 and hsa-mir-1909) and the variant allele T may result in an increased MFE for these two miRNAs. Variant c.1187 A>G was not predicted to be located in the target site of any miRNAs (Table 5).


Melanocortin peptides such as ACTH, MSH and β-endorphin are derived from the precursor molecule POMC, which is encoded by the POMC gene (POMC). Variation in POMC may contribute to SD [20], [21], [22], [23] and obesity [24], [26]. The present study focused on the identification of POMC exonic variants and the association of POMC exonic variants with SD (specifically, AD, CD, OD and/or MjD) and BMI. Our findings suggest that POMC exonic variants may influence risk for both SD and BMI. Nevertheless, POMC common and rare variants may exert different effects on these two phenotypes.

In this study, two common POMC exonic variants [the 9-bp insertion/deletion polymorphism (c.560_561 9-bp Ins or rs10654394) in exon 4 and c.1130 G>T (or rs1042571) in the 3′ UTR] were identified and their association with SD or BMI was analyzed. A positive association between common variant rs1042571 and BMI was observed in EAs but not AAs. In a previous study, the common variant rs10654394 (or the 9-bp insertion/deletion polymorphism) was associated with obesity in children [26]. However, in the present study, we saw no association between this variant and BMI in either population. Additionally, the two common variants were not associated with SD in either AAs or EAs. In fact, no published studies have demonstrated an association between these two common exonic variants and SD. To understand whether our sample had sufficient statistical power to detect the association between these two common variants and SD or BMI, a retrospective power analysis was conducted. Assuming the statistical power equals to 80% and the type I error equals to 0.05, the minimum effect size (or odds ratio) that is detectable would be 1.67 (in AAs) and 2.50 (in EAs) for variant rs10654394 and 2.00 (in AAs) and 1.68 (in EAs) for variant rs1042571. Considering that the sample size for this study was moderate, further studies with a larger sample are needed to validate the above findings.

Moreover, the two common variants are potentially functional. The 9-bp insertion allele of rs10654394 (c.560_561 9-bp Ins) leads to three extra amino acids (Ser-Ser-Gly) at the carboxyl terminus of γ-MSH in the conserved region of the 16 KD fragment (Figure 1). This variant may influence mRNA stability or post-translational cleavage of the POMC peptide. The variant rs1042571 (c.1130 G>T) in the 3′ UTR was predicted to be located in the binding site of two miRNAs (hsa-mir-3715 and hsa-mir-1909) (Table 5). The variant allele T potentially increases the minimal free energy (MFE) for hybridization of these two miRNAs (hsa-mir-3715: from −25.0 kcal/mol to −22.1 kcal/mol; hsa-mir-1909: from −26.7 kcal/mol to −24.9 kcal/mol) to the target sequence in POMC 3′ UTR, thus reducing the binding of miRNAs to POMC 3′ UTR and increasing POMC expression.

Additionally, we identified 21 rare variants in POMC exons, most of which were predicted to cause a functional or structural change in POMC as described in Table 5. When these exonic variants were collapsed into a single variable and the frequency of which was compared between SD cases and controls, a significantly higher frequency of these rare variants was observed in AA cases affected with SD in general or AD, CD or MjD specifically than in AA controls. Nevertheless, in EAs, these rare variants were not significantly more frequent in SD cases than controls (Table 4). Thus, rare POMC exonic variants may contribute to the etiology of SD in a population-specific pattern. Moreover, the association of these rare variants with overweight or obesity was analyzed using the same approach but no positive findings were obtained (Table 4). In this study, 21 rare variants were included in Fisher’s exact tests. Because they were collapsed into a single variable to test their cumulative influence on SD or BMI, corrections for multiple testing (e.g., Bonferroni correction) were not applied to adjust the significance P values.

Although POMC rare variants were not found to be associated with BMI using the collapsing method and the Fisher’s exact test, we cannot exclude the possibility that certain rare POMC exonic variants may influence overweight or obesity. As shown in Figure 2, eight rare variants (c.267 C>T, c.343 G>A, c.432 G>T, c.609 C>T, c.969 C>G, c.989 C>A, c.1095 T deletion, and c.1187 A>G) were identified only in controls, and each of them was identified in a single heterozygous carrier. Four of them (c.432 G>T, c.609 C>T, c.969 C>G and c.1187 A>G) were only present in control subjects with BMI >25.0, suggesting that these rare variants may increase risk of becoming overweight or obese specifically. In addition, rare POMC exonic variants may confer risk for SD and overweight or obesity through a shared mechanism. As shown in Figure 2, five rare variants (c.259 G>A, c.524 C>A, c.656 C>A, c.685 T>G, and c.867_868 6-bp Ins) were identified only in SD cases who had overweight or obesity. Thus, these rare POMC exonic variants may influence the vulnerability to SD and overweight or obesity via a common biological pathway.

As we know, exonic rare variants may have a larger impact on gene transcription or protein activity than intronic or intergenic variants, and in some cases, they may cause the disease directly. Hence, it would be important to explore the biological function of the rare variants or mutations in POMC exons. As shown in Table 5, POMC exonic variants may change the affinity of transcription factors for their binding sites (e.g., variants in the 5′ UTR), introduce premature stop codons (e.g., nonsense mutations in coding regions), alter amino acid sequences (e.g., missense mutations in coding regions), or change the activity of regulatory miRNAs (e.g., variants in the 3′ UTR).

The findings from this study should be seen in the context of three main limitations: (1) a moderate sample size, (2) the function of exonic variants was predicted only by bioinformatics analyses, and (3) the population stratification issue was not considered in rare variant data analysis. A larger sample would increase the potential to identify POMC rare variants and provide a greater statistical power to analyze the association of identified variants and diseases. This is particularly important for detection of small or moderate effects of genes involved in complex disorders. Because the etiological role of the rare variants identified in this study was only predicted using bioinformatics analysis, experimental studies in vivo or in vitro are needed to validate their biological function. Another challenge was to control the influence of potential population substructure bias on association analysis results. Matheson and McVean [40] pointed out that population structure could alter the conclusion of genetic association studies involving either common or rare variants. Principle component analysis is the usual method used to control for population structure in genome-wide association studies, in which the top components are considered as covariates. In this study, we used the ancestry coefficients obtained from running STRUCTURE [32], [33] on 41 ancestry informative markers to exclude participants whose genetic ancestry conflicted with their self-reported race and took the score as a continuous covariate to adjust for the genetic background noise for the two identified common variants. However, there is no accepted method to correct the effect of population stratification on rare variant analysis results per our current design. One possible way, as suggested by Matheson and McVean [40], is that family-based association tests, by inclusion of other family members of participants, might help minimize the influence of structured populations on association analysis results.

In conclusion, this study provides evidence that both common and rare variants in POMC could increase the risk for SD and/or obesity. It also suggests that certain POMC variants may influence the vulnerability to SD and overweight or obesity in a common or specific biological pathway. Functional studies are needed to elucidate the shared or specific molecular mechanisms by which POMC variation influences the susceptibility to SD and/or obesity.

Supporting Information

Figure S1

Conserved analysis of POMC protein sequence across multiple species.


Table S1

Primers for sequencing four POMC exons.


Table S2

PCR conditions for amplifying four POMC exons.


Table S3

Haplotype association of two common POMC variants (rs1654394-rs1042571) and BMI.


Table S4

Allelic association of two POMC common variants and substance dependence (SD) traits.


Table S5

Haplotype association of two POMC common variants (rs10654394-rs1042571) and substance dependence (SD).


Funding Statement

This study was supported by the National Institutes of Health (NIH) Grants K99/R00 DA022891 (HZ), R01 DA12690 (JG), R01 AA11330 (JG), R01 DA12849 (JG), R01 AA017535 (JG), R01 DA018432 (HRK), K24 AA013736 (HRK) and P50 AA012870 (JG); and grants from the Veterans Affairs (VA) Connecticut (CT) Mental Illness Research, Education and Clinical Center (MIRECC) to JG; and the Alcoholic Beverage Medical Research Foundation (ABMRF) to HZ. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Johnson PM, Kenny PJ (2010) Dopamine D2 receptors in addiction-like reward dysfunction and compulsive eating in obese rats. Nat Neurosci 13: 635–641 [PMC free article] [PubMed]
2. Kleiner KD, Gold MS, Frost-Pineda K, Lenz-Brunsman B, Perri MG, et al. (2004) Body mass index and alcohol use. J Addict Dis 23: 105–118 [PubMed]
3. Breslow RA, Smothers BA (2005) Drinking patterns and body mass index in never smokers: National Health Interview Survey, 1997–2001. Am J Epidemiol 161: 368–376 [PubMed]
4. Gearhardt AN, Corbin WR (2009) Body mass index and alcohol consumption: family history of alcoholism as a moderator. Psychol Addict Behav 23: 216–225 [PubMed]
5. Duncan AE, Grant JD, Bucholz KK, Madden PA, Heath AC (2009) Relationship between body mass index, alcohol use, and alcohol misuse in a young adult female twin sample. J Stud Alcohol Drugs 70: 458–466 [PubMed]
6. Bluml V, Kapusta N, Vyssoki B, Kogoj D, Walter H, et al. (2011) Relationship between substance use and body mass index in young males. Am J Addict 21: 72–77 [PubMed]
7. Arif AA, Rohrer JE (2005) Patterns of alcohol drinking and its association with obesity: data from the Third National Health and Nutrition Examination Survey, 1988–1994. BMC Public Health 5: 126. [PMC free article] [PubMed]
8. Pickering RP, Grant BF, Chou SP, Compton WM (2007) Are overweight, obesity, and extreme obesity associated with psychopathology? Results from the national epidemiologic survey on alcohol and related conditions. J Clin Psychiatry 68: 998–1009 [PubMed]
9. Barry D, Petry NM (2009) Associations between body mass index and substance use disorders differ by gender: results from the National Epidemiologic Survey on Alcohol and Related Conditions. Addict Behav 34: 51–60 [PMC free article] [PubMed]
10. Owerbach D, Rutter WJ, Roberts JL, Whitfeld P, Shine J, et al. (1981) The proopiocortin (adrenocorticotropin/beta-lipoprotein) gene is located on chromosome 2 in humans. Somatic Cell Genet 7: 359–369 [PubMed]
11. Zabel BU, Naylor SL, Sakaguchi AY, Bell GI, Shows TB (1983) High-resolution chromosomal localization of human genes for amylase, proopiomelanocortin, somatostatin, and a DNA fragment (D3S1) by in situ hybridization. Proc Natl Acad Sci U S A 80: 6932–6936 [PubMed]
12. Hartwig AC (1991) Peripheral beta-endorphin and pain modulation. Anesth Prog 38: 75–78 [PMC free article] [PubMed]
13. Ellacott KL, Cone RD (2004) The central melanocortin system and the integration of short- and long-term regulators of energy homeostasis. Recent Prog Horm Res 59: 395–408 [PubMed]
14. Millington GW (2006) Proopiomelanocortin (POMC): the cutaneous roles of its melanocortin products and receptors. Clin Exp Dermatol 31: 407–412 [PubMed]
15. Heijnen CJ, Zijlstra J, Kavelaars A, Croiset G, Ballieux RE (1987) Modulation of the immune response by POMC-derived peptides. I. Influence on proliferation of human lymphocytes. Brain Behav Immun 1: 284–291 [PubMed]
16. Gianoulakis C, Dai X, Brown T (2003) Effect of chronic alcohol consumption on the activity of the hypothalamic-pituitary-adrenal axis and pituitary beta-endorphin as a function of alcohol intake, age, and gender. Alcohol Clin Exp Res 27: 410–423 [PubMed]
17. Marchesi C, Chiodera P, Ampollini P, Volpi R, Coiro V (1997) Beta-endorphin, adrenocorticotropic hormone and cortisol secretion in abstinent alcoholics. Psychiatry Res 72: 187–194 [PubMed]
18. Zalewska-Kaszubska J, Czarnecka E (2005) Deficit in beta-endorphin peptide and tendency to alcohol abuse. Peptides 26: 701–705 [PubMed]
19. Weaver JU, Kopelman PG, McLoughlin L, Forsling ML, Grossman A (1993) Hyperactivity of the hypothalamo-pituitary-adrenal axis in obesity: a study of ACTH, AVP, beta-lipotrophin and cortisol responses to insulin-induced hypoglycaemia. Clin Endocrinol (Oxf) 39: 345–350 [PubMed]
20. Xuei X, Flury-Wetherill L, Bierut L, Dick D, Nurnberger J Jr, et al. (2007) The opioid system in alcohol and drug dependence: family-based association study. Am J Med Genet B Neuropsychiatr Genet 144B: 877–884 [PubMed]
21. Racz I, Schurmann B, Karpushova A, Reuter M, Cichon S, et al. (2008) The opioid peptides enkephalin and beta-endorphin in alcohol dependence. Biol Psychiatry 64: 989–997 [PMC free article] [PubMed]
22. Zhang H, Kranzler HR, Weiss RD, Luo X, Brady KT, et al. (2009) Pro-opiomelanocortin gene variation related to alcohol or drug dependence: evidence and replications across family- and population-based studies. Biol Psychiatry 66: 128–136 [PMC free article] [PubMed]
23. Kendler KS, Kalsi G, Holmans PA, Sanders AR, Aggen SH, et al. (2011) Genomewide association analysis of symptoms of alcohol dependence in the molecular genetics of schizophrenia (MGS2) control sample. Alcohol Clin Exp Res 35: 963–975 [PMC free article] [PubMed]
24. Krude H, Biebermann H, Luck W, Horn R, Brabant G, et al. (1998) Severe early-onset obesity, adrenal insufficiency and red hair pigmentation caused by POMC mutations in humans. Nat Genet 19: 155–157 [PubMed]
25. Delplanque J, Barat-Houari M, Dina C, Gallina P, Clement K, et al. (2000) Linkage and association studies between the proopiomelanocortin (POMC) gene and obesity in caucasian families. Diabetologia 43: 1554–1557 [PubMed]
26. Miraglia del Giudice E, Cirillo G, Santoro N, D'Urso L, Carbone MT, et al. (2001) Molecular screening of the proopiomelanocortin (POMC ) gene in Italian obese children: report of three new mutations. Int J Obes Relat Metab Disord 25: 61–67 [PubMed]
27. Rosmond R, Ukkola O, Bouchard C, Bjorntorp P (2002) Polymorphisms in exon 3 of the proopiomelanocortin gene in relation to serum leptin, salivary cortisol, and obesity in Swedish men. Metabolism 51: 642–644 [PubMed]
28. Echwald SM, Sorensen TI, Andersen T, Tybjaerg-Hansen A, Clausen JO, et al. (1999) Mutational analysis of the proopiomelanocortin gene in Caucasians with early onset obesity. Int J Obes Relat Metab Disord 23: 293–298 [PubMed]
29. Pierucci-Lagha A, Gelernter J, Chan G, Arias A, Cubells JF, et al. (2007) Reliability of DSM-IV diagnostic criteria using the semi-structured assessment for drug dependence and alcoholism (SSADDA). Drug Alcohol Depend 91: 85–90 [PMC free article] [PubMed]
30. Pierucci-Lagha A, Gelernter J, Feinn R, Cubells JF, Pearson D, et al. (2005) Diagnostic reliability of the Semi-structured Assessment for Drug Dependence and Alcoholism (SSADDA). Drug Alcohol Depend 80: 303–312 [PubMed]
31. Association AP (1994) Diagnostic and StatisticalManual of Mental Disorders (4th ed.). Washington, DC: Author.
32. Falush D, Stephens M, Pritchard JK (2003) Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164: 1567–1587 [PubMed]
33. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959 [PubMed]
34. Yang BZ, Zhao H, Kranzler HR, Gelernter J (2005) Practical population group assignment with selected informative markers: characteristics and properties of Bayesian clustering via STRUCTURE. Genet Epidemiol 28: 302–312 [PubMed]
35. Xie P, Kranzler HR, Poling J, Stein MB, Anton RF, et al. (2009) Interactive effect of stressful life events and the serotonin transporter 5-HTTLPR genotype on posttraumatic stress disorder diagnosis in 2 independent populations. Arch Gen Psychiatry 66: 1201–1209 [PMC free article] [PubMed]
36. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575 [PubMed]
37. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265 [PubMed]
38. Dupont WD, Plummer WD Jr (1990) Power and sample size calculations. A review and computer program. Control Clin Trials 11: 116–128 [PubMed]
39. Dickson SP, Wang K, Krantz I, Hakonarson H, Goldstein DB (2010) Rare variants create synthetic genome-wide associations. PLoS Biol 8: e1000294. [PMC free article] [PubMed]
40. Mathieson I, McVean G (2012) Differential confounding of rare and common variants in spatially structured populations. Nat Genet 44: 243–246 [PMC free article] [PubMed]
41. Tajima F (1989) DNA polymorphism in a subdivided population: the expected number of segregating sites in the two-subpopulation model. Genetics 123: 229–240 [PubMed]
42. Xie P, Kranzler HR, Krauthammer M, Cosgrove KP, Oslin D, et al. (2011) Rare nonsynonymous variants in alpha-4 nicotinic acetylcholine receptor gene protect against nicotine dependence. Biol Psychiatry 70: 528–536 [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science