Identification of methylation quantitative trait loci (mQTLs) influencing promoter DNA methylation of alcohol dependence risk genes

Interaction of DNA methylation and sequence variants that are methylation quantitative trait loci (mQTLs) may influence susceptibility to diseases such as alcohol dependence (AD). We used genome-wide genotype data from 268 African Americans (AAs: 129 AD cases and 139 controls) and 143 European Americans (EAs: 129 AD cases and 14 controls) to identify mQTLs that were associated with promoter CpGs in 82 AD risk genes. 282 significant mQTL–CpG pairs (9.9 × 10^?100 ≤ P _nominal ≤ 7.7 × 10^?8) in AAs and 313 significant mQTL–CpG pairs (2.7 × 10^?53 ≤ P _nominal ≤ 9.9 × 10^?8) in EAs were identified [i.e., mQTL–CpG associations survived multiple-testing correction, q values (false discovery rate) ≤ 0.05]. The most significant mQTL was rs1800759, which was strongly associated with CpG cg12011299 in both AAs (P _nominal = 9.9 × 10^?100; q = 6.7 × 10^?91) and EAs (P _nominal = 2.7 × 10^?53; q = 1.4 × 10^?44). Rs1800759 (previously known to be associated to AD) and CpG cg12011299 (distance: 37 bp) are both located in alcohol dehydrogenase (ADH) 4 gene (ADH4) promoter region. In general, the strength of association between mQTLs and CpGs was inversely correlated with the distance between them. Association was also influenced by race and AD. Additionally, 48.3 % of the mQTLs identified in AAs and 65.6 % of the mQTLs identified in EAs were predicted to be expression QTLs. Three mQTLs (rs2173201, rs4147542, and rs4147541 in ADH1B-AHD1C gene cluster region) found in AAs were previously identified by our genome-wide association studies as being significantly associated with AD in AAs. Thus, DNA methylation, which can be influenced by sequence variants and is implicated in gene expression regulation, appears to at least partially underlie the association of genetic variation with AD.     

A Weighted False Discovery Rate Control Procedure Reveals Alleles at FOXA2 that Influence Fasting Glucose Levels

Association signals in GWAS are usually prioritized solely by p values. Here, we attempt to improve the power of GWAS by using a weighted false discovery rate control procedure to detect associations of low-frequency variants with effect sizes similar to or even larger than those of common variants. We used the Affymetrix Genome-Wide Human SNP Array 6.0 to test for association with fasting glucose levels in the Atherosclerosis Risk in Communities Study (ARIC) population. In addition to finding several previously identified sequence variations, we identified a low-frequency variant (rs1209523; minor allele frequency = 0.043) near FOXA2 that was associated with fasting glucose levels in European Americans (EAs) (n = 7428, p value = 1.3 × 10⁻⁵). The association between rs1209523 and glucose levels was also significant in African Americans (AAs) (n = 2029, p value = 6.7 × 10⁻³) of the ARIC and was confirmed by replication in both EAs and AAs of the Dallas Heart Study (n = 963 and 1571, respectively; p values = 5.3 × 10⁻³ and 5.8 × 10⁻⁴, respectively) and in EAs of the Cooper Center Longitudinal Study (n = 2862; p value = 1.6 × 10⁻²). A meta-analysis of these five populations yielded an estimated effect size of −1.31 mg/dl per minor allele (p value = 2.2 × 10⁻¹¹). This study reveals that there is a cache of less-frequent variants in GWAS arrays that can be identified via analytical approaches accounting for allele frequencies.     

A genome‐wide association study of coping behaviors suggests FBXO45 is associated with emotional expression

Individuals use coping behaviors to deal with unpleasant daily events. Such behaviors can moderate or mediate the pathway between psychosocial stress and health‐related outcomes. However, few studies have examined the associations between coping behaviors and genetic variants. We conducted a genome‐wide association study (GWAS) on coping behaviors in 14088 participants aged 35 to 69 years as part of the Japan Multi‐Institutional Collaborative Cohort Study. Five coping behaviors (emotional expression, emotional support seeking, positive reappraisal, problem solving and disengagement) were measured and analyzed. A GWAS analysis was performed using a mixed linear model adjusted for study area, age and sex. Variants with suggestive significance in the discovery phase (N = 6403) were further examined in the replication phase (N = 7685). We then combined variant‐level association evidence into gene‐level evidence using a gene‐based analysis. The results showed a significant genetic contribution to emotional expression and disengagement, with an estimation that the 19.5% and 6.6% variance in the liability‐scale was explained by common variants. In the discovery phase, 12 variants met suggestive significance (P < 1 × 10^?6) for association with the coping behaviors and perceived stress. However, none of these associations were confirmed in the replication stage. In gene‐based analysis, FBXO45, a gene with regulatory roles in synapse maturation, was significantly associated with emotional expression after multiple corrections (P < 3.1 × 10^?6). In conclusion, our results showed the existence of up to 20% genetic contribution to coping behaviors. Moreover, our gene‐based analysis using GWAS data suggests that genetic variations in FBXO45 are associated with emotional expression.     

A genome‐wide association study of bipolar disorder and comorbid migraine

Both migraine and bipolar affective disorder (BPAD) are complex phenotypes with significant genetic and nongenetic components. Epidemiological and clinical studies have showed a high degree of comorbidity between migraine and BPAD, and overlapping regions of linkage have been shown in numerous genome‐wide linkage studies. To identify susceptibility factors for the BPAD/migraine phenotype, we conducted a genome‐wide association study (GWAS) in 1001 cases with bipolar disorder collected through the NIMH Genetics Initiative for Bipolar Disorder and genotyped at 1 m single‐nucleotide polymorphisms (SNPs) as part of the Genetic Association Information Network (GAIN). We compared BPAD patients without any headache (n = 699) with BPAD patients with doctor diagnosed migraine (n = 56). The strongest evidence for association was found for several SNPs in a 317‐kb region encompassing the uncharacterized geneKIAA0564 {e.g. rs9566845 [OR = 4.98 (95% CI: 2.6–9.48), P = 7.7 × 10^?8] and rs9566867 (P = 8.2 × 10^?8)}. Although the level of signficance was significantly reduced when using the Fisher's exact test (as a result of the low count of cases with migraine), rs9566845 P = 1.4 × 10^?5 and rs9566867 P = 1.5 × 10^?5, this region remained the most prominent finding. Furthermore, marker rs9566845 was genotyped and found associated with migraine in an independent Norwegian sample of adult attention deficit hyperactivity disorder (ADHD) patients with and without comorbid migraine (n = 131 and n = 324, respectively), OR = 2.42 (1.18–4.97), P = 0.013. This is the first GWAS examining patients with bipolar disorder and comorbid migraine. These data suggest that genetic variants in the KIAA0564 gene region may predispose to migraine headaches in subgroups of patients with both BPAD and ADHD.     

Multiple cholinergic nicotinic receptor genes affect nicotine dependence risk in African and European Americans

Several independent studies show that the chromosome 15q25.1 region, which contains the CHRNA5–CHRNA3–CHRNB4 gene cluster, harbors variants strongly associated with nicotine dependence, other smoking behaviors, lung cancer and chronic obstructive pulmonary disease. We investigated whether variants in other cholinergic nicotinic receptor subunit (CHRN) genes affect the risk of nicotine dependence in a new sample of African Americans (AAs) (N = 710). We also analyzed this AA sample together with a European American (EA) sample (N = 2062, 1608 of which have been previously studied), allowing for differing effects in the two populations. Cases are current nicotine‐dependent smokers and controls are non‐dependent smokers. Variants in or near CHRND–CHRNG, CHRNA7 and CHRNA10 show modest association with nicotine dependence risk in the AA sample. In addition, CHRNA4, CHRNB3–CHRNA6 and CHRNB1 show association in at least one population. CHRNG and CHRNA4 harbor single nucleotide polymorphisms (SNPs) that have opposite directions of effect in the two populations. In each of the population samples, these loci substantially increase the trait variation explained, although no loci meet Bonferroni‐corrected significance in the AA sample alone. The trait variation explained by three key associated SNPs in CHRNA5–CHRNA3–CHRNB4 is 1.9% in EAs and also 1.9% in AAs; this increases to 4.5% in EAs and 7.3% in AAs when we add six variants representing associations at other CHRN genes. Multiple nicotinic receptor subunit genes outside chromosome 15q25 are likely to be important in the biological processes and development of nicotine dependence, and some of these risks may be shared across diverse populations.     

Extreme‐phenotype genome‐wide association study (XP‐GWAS): a method for identifying trait‐associated variants by sequencing pools of individuals selected from a diversity panel

Although approaches for performing genome‐wide association studies (GWAS) are well developed, conventional GWAS requires high‐density genotyping of large numbers of individuals from a diversity panel. Here we report a method for performing GWAS that does not require genotyping of large numbers of individuals. Instead XP‐GWAS (extreme‐phenotype GWAS) relies on genotyping pools of individuals from a diversity panel that have extreme phenotypes. This analysis measures allele frequencies in the extreme pools, enabling discovery of associations between genetic variants and traits of interest. This method was evaluated in maize (Zea mays) using the well‐characterized kernel row number trait, which was selected to enable comparisons between the results of XP‐GWAS and conventional GWAS. An exome‐sequencing strategy was used to focus sequencing resources on genes and their flanking regions. A total of 0.94 million variants were identified and served as evaluation markers; comparisons among pools showed that 145 of these variants were statistically associated with the kernel row number phenotype. These trait‐associated variants were significantly enriched in regions identified by conventional GWAS. XP‐GWAS was able to resolve several linked QTL and detect trait‐associated variants within a single gene under a QTL peak. XP‐GWAS is expected to be particularly valuable for detecting genes or alleles responsible for quantitative variation in species for which extensive genotyping resources are not available, such as wild progenitors of crops, orphan crops, and other poorly characterized species such as those of ecological interest.     

Racial differences in genome-wide methylation profiling and gene expression in breast tissues from healthy women

Breast cancer is more common in European Americans (EAs) than in African Americans (AAs) but mortality from breast cancer is higher among AAs. While there are racial differences in DNA methylation and gene expression in breast tumors, little is known whether such racial differences exist in breast tissues of healthy women. Genome-wide DNA methylation and gene expression profiling was performed in histologically normal breast tissues of healthy women. Linear regression models were used to identify differentially-methylated CpG sites (CpGs) between EAs (n = 61) and AAs (n = 22). Correlations for methylation and expression were assessed. Biological functions of the differentially-methylated genes were assigned using the Ingenuity Pathway Analysis. Among 485 differentially-methylated CpGs by race, 203 were hypermethylated in EAs, and 282 were hypermethylated in AAs. Promoter-related differentially-methylated CpGs were more frequently hypermethylated in EAs (52%) than AAs (27%) while gene body and intergenic CpGs were more frequently hypermethylated in AAs. The differentially-methylated CpGs were enriched for cancer-associated genes with roles in cell death and survival, cellular development, and cell-to-cell signaling. In a separate analysis for correlation in EAs and AAs, different patterns of correlation were found between EAs and AAs. The correlated genes showed different biological networks between EAs and AAs; networks were connected by Ubiquitin C. To our knowledge, this is the first comprehensive genome-wide study to identify differences in methylation and gene expression between EAs and AAs in breast tissues from healthy women. These findings may provide further insights regarding the contribution of epigenetic differences to racial disparities in breast cancer.     

Profile of common prostate cancer risk variants in an unscreened Romanian population

To find sequence variants affecting prostate cancer (PCA) susceptibility in an unscreened Romanian population we use a genome‐wide association study (GWAS). The study population included 990 unrelated pathologically confirmed PCA cases and 1034 male controls. DNA was genotyped using Illumina SNP arrays, and 24.295.558 variants were imputed using the 1000 Genomes data set. An association test was performed between the imputed markers and PCA. A systematic literature review for variants associated with PCA risk identified 115 unique variants that were tested in the Romanian sample set. Thirty of the previously reported SNPs replicated (P‐value < 0.05), with the strongest associations observed at: 8q24.21, 11q13.3, 6q25.3, 5p15.33, 22q13.2, 17q12 and 3q13.2. The replicated variants showing the most significant association in Romania are rs1016343 at 8q24.21 (P = 2.2 × 10^?4), rs7929962 at 11q13.3 (P = 2.7 × 10^?4) and rs9364554 at 6q25.2 (P = 4.7 × 10^?4). None of the variants tested in the Romanian GWAS reached genome‐wide significance (P‐value <5 × 10^?8) but 807 markers had P‐values <1 × 10^?4. Here, we report the results of the first GWAS of PCA performed in a Romanian population. Our study provides evidence that a substantial fraction of previously validated PCA variants associate with risk in this unscreened Romanian population.     

Genome‐wide association studies of alcohol dependence,DSM‐IV criterion count and individual criteria

Genome‐wide association studies (GWAS) of alcohol dependence (AD) have reliably identified variation within alcohol metabolizing genes (eg, ADH1B) but have inconsistently located other signals, which may be partially attributable to symptom heterogeneity underlying the disorder. We conducted GWAS of DSM‐IV AD (primary analysis), DSM‐IV AD criterion count (secondary analysis), and individual dependence criteria (tertiary analysis) among 7418 (1121 families) European American (EA) individuals from the Collaborative Study on the Genetics of Alcoholism (COGA). Trans‐ancestral meta‐analyses combined these results with data from 3175 (585 families) African‐American (AA) individuals from COGA. In the EA GWAS, three loci were genome‐wide significant: rs1229984 in ADH1B for AD criterion count (P = 4.16E?11) and Desire to cut drinking (P = 1.21E?11); rs188227250 (chromosome 8, Drinking more than intended, P = 6.72E?09); rs1912461 (chromosome 15, Time spent drinking, P = 1.77E?08). In the trans‐ancestral meta‐analysis, rs1229984 was associated with multiple phenotypes and two additional loci were genome‐wide significant: rs61826952 (chromosome 1, DSM‐IV AD, P = 8.42E?11); rs7597960 (chromosome 2, Time spent drinking, P = 1.22E?08). Associations with rs1229984 and rs18822750 were replicated in independent datasets. Polygenic risk scores derived from the EA GWAS of AD predicted AD in two EA datasets (P < .01; 0.61%‐1.82% of variance). Identified novel variants (ie, rs1912461, rs61826952) were associated with differential central evoked theta power (loss ? gain; P = .0037) and reward‐related ventral striatum reactivity (P = .008), respectively. This study suggests that studying individual criteria may unveil new insights into the genetic etiology of AD liability.     

Genome‐wide association study identifies loci associated with liability to alcohol and drug dependence that is associated with variability in reward‐related ventral striatum activity in African‐ and European‐Americans

Genetic influences on alcohol and drug dependence partially overlap, however, specific loci underlying this overlap remain unclear. We conducted a genome‐wide association study (GWAS) of a phenotype representing alcohol or illicit drug dependence (ANYDEP) among 7291 European‐Americans (EA; 2927 cases) and 3132 African‐Americans (AA: 1315 cases) participating in the family‐based Collaborative Study on the Genetics of Alcoholism. ANYDEP was heritable (h ² in EA = 0.60, AA = 0.37). The AA GWAS identified three regions with genome‐wide significant (GWS; P < 5E‐08) single nucleotide polymorphisms (SNPs) on chromosomes 3 (rs34066662, rs58801820) and 13 (rs75168521, rs78886294), and an insertion‐deletion on chromosome 5 (chr5:141988181). No polymorphisms reached GWS in the EA. One GWS region (chromosome 1: rs1890881) emerged from a trans‐ancestral meta‐analysis (EA + AA) of ANYDEP, and was attributable to alcohol dependence in both samples. Four genes (AA: CRKL, DZIP3, SBK3; EA: P2RX6) and four sets of genes were significantly enriched within biological pathways for hemostasis and signal transduction. GWS signals did not replicate in two independent samples but there was weak evidence for association between rs1890881 and alcohol intake in the UK Biobank. Among 118 AA and 481 EA individuals from the Duke Neurogenetics Study, rs75168521 and rs1890881 genotypes were associated with variability in reward‐related ventral striatum activation. This study identified novel loci for substance dependence and provides preliminary evidence that these variants are also associated with individual differences in neural reward reactivity. Gene discovery efforts in non‐European samples with distinct patterns of substance use may lead to the identification of novel ancestry‐specific genetic markers of risk.     

Genetic loci associated with skin pigmentation in African Americans and their effects on vitamin D deficiency

A recent genome-wide association study (GWAS) in African descent populations identified novel loci associated with skin pigmentation. However, how genomic variations affect skin pigmentation and how these skin pigmentation gene variants affect serum 25(OH) vitamin D variation has not been explored in African Americans (AAs). In order to further understand genetic factors that affect human skin pigmentation and serum 25(OH)D variation, we performed a GWAS for skin pigmentation with 395 AAs and a replication study with 681 AAs. Then, we tested if the identified variants are associated with serum 25(OH) D concentrations in a subset of AAs (n = 591). Skin pigmentation, Melanin Index (M-Index), was measured using a narrow-band reflectometer. Multiple regression analysis was performed to identify variants associated with M-Index and to assess their role in serum 25(OH)D variation adjusting for population stratification and relevant confounding variables. A variant near the SLC24A5 gene (rs2675345) showed the strongest signal of association with M-Index (P = 4.0 x 10⁻³⁰ in the pooled dataset). Variants in SLC24A5, SLC45A2 and OCA2 together account for a large proportion of skin pigmentation variance (11%). The effects of these variants on M-Index was modified by sex (P for interaction = 0.009). However, West African Ancestry (WAA) also accounts for a large proportion of M-Index variance (23%). M-Index also varies among AAs with high WAA and high Genetic Score calculated from top variants associated with M-Index, suggesting that other unknown genomic factors related to WAA are likely contributing to skin pigmentation variation. M-Index was not associated with serum 25(OH)D concentrations, but the Genetic Score was significantly associated with vitamin D deficiency (serum 25(OH)D levels less than 12 ng/mL) (OR, 1.30; 95% CI, 1.04–1.64). The findings support the hypothesis suggesting that skin pigmentation evolved responding to increased demand for subcutaneous vitamin D synthesis in high latitude environments.     

Common vitamin D pathway gene variants reveal contrasting effects on serum vitamin D levels in African Americans and European Americans

Vitamin D deficiency is more common among African Americans (AAs) than among European Americans (EAs), and epidemiologic evidence links vitamin D status to many health outcomes. Two genome-wide association studies (GWAS) in European populations identified vitamin D pathway gene single-nucleotide polymorphisms (SNPs) associated with serum vitamin D [25(OH)D] levels, but a few of these SNPs have been replicated in AAs. Here, we investigated the associations of 39 SNPs in vitamin D pathway genes, including 19 GWAS-identified SNPs, with serum 25(OH)D concentrations in 652 AAs and 405 EAs. Linear and logistic regression analyses were performed adjusting for relevant environmental and biological factors. The pattern of SNP associations was distinct between AAs and EAs. In AAs, six GWAS-identified SNPs in GC, CYP2R1, and DHCR7/NADSYN1 were replicated, while nine GWAS SNPs in GC and CYP2R1 were replicated in EAs. A CYP2R1 SNP, rs12794714, exhibited the strongest signal of association in AAs. In EAs, however, a different CYP2R1 SNP, rs1993116, was the most strongly associated. Our models, which take into account genetic and environmental variables, accounted for 20 and 28 % of the variance in serum vitamin D levels in AAs and EAs, respectively.     

Association of FTO Gene Variants With Adiposity in African‐American Adolescents

The prevalence of obesity continues to increase significantly, with the largest rise in the African‐American adolescents. Genetic contributions to obesity are being identified with the advent of genome‐wide association studies (GWAS). Specifically, variants of the fat mass and obesity associated (FTO) gene have been associated with obesity in populations of European descent. The studies in African Americans have been inconclusive. To further evaluate the association of the FTO gene and adiposity in African Americans, we genotyped 47 single‐nucleotide polymorphisms (SNPs), including seven SNPs previously reported to be significant in the literature in a cohort consisting of 561 non‐Hispanic white and 497 African‐American individuals. Analysis of our data showed 17 SNPs to be associated with BMI Z‐score (BMI‐Z) in our study population. The strongest association was found in the African Americans. The most significant SNP was rs8057044, which was associated with BMI‐Z in the African Americans (P = 0.00054). SNP rs9939609 was found to be significant in the non‐Hispanic white population (P = 0.028). Our data confirm the association between FTO and adiposity suggesting that FTO is a childhood obesity susceptibility gene. Our data also identify a novel SNP of the FTO gene (rs8057044) that is associated with measures of adiposity in the African‐American population.     

Genome‐wide association study of genetic factors related to confectionery intake: Potential roles of the ADIPOQ gene

Objective: The excessive consumption of confectionery might have adverse effects on human health. To screen genetic factors associated with confectionery‐intake frequency, a genome‐wide association study (GWAS) in Japan was conducted. Design and Methods: For the discovery phase (stage 1), we conducted a GWAS of 939 noncancer patients in a cancer hospital. Additive models were used to test associations between genotypes of approximately 500,000 single‐nucleotide polymorphisms (SNPs) and the confectionery‐intake score (based on intake frequency). We followed‐up association signals with P < 1 × 10^?5 and minor allele frequency >0.01 in stage 1 by genotyping the SNPs of 4,491 participants in a cross‐sectional study within a cohort (replication phase [stage 2]). Results: We identified 12 SNPs in stage 1 that were potentially related to confectionery intake. In stage 2, this association was replicated for one SNP (rs822396; P = 0.049 for stage 2 and 4.2 × 10^?5 for stage 1+2) in intron 1 of the ADIPOQ gene, which encodes the adipokine adiponectin. Conclusions: Given the biological plausibility and previous relevant findings, the association of an SNP in the ADIPOQ gene with a preference for confectionery is worthy of follow‐up and provides a good working hypothesis for experimental testing.     

Associations of Free Fatty Acids With Insulin Secretion and Action Among African‐American and European‐American Girls and Women

Ethnic differences in insulin secretion and action between African Americans (AAs) and European Americans (EAs) may influence mobilization of free fatty acids (FFAs). We tested the hypotheses that FFA concentrations would be associated with measures of insulin secretion and action before and during a glucose challenge test. Subjects were 48 prepubertal girls, 60 premenopausal women, and 46 postmenopausal women. Fasting insulin (insulin₀), the acute insulin response to glucose (AIR_g), the insulin sensitivity index (S_I), basal and nadir FFA (FFA₀, FFA_nadir), and nadir time (TIME_nadir) were determined during an intravenous glucose tolerance test (IVGTT). Stepwise multiple linear regression (MLR) analysis was conducted to identify associations of FFA₀, FFA_nadir, and TIME_nadir with ethnicity, age group, insulin measures, indexes of body composition from dual‐energy X‐ray absorptiometry, and measures of fat distribution from computed tomography scan. In this population, insulin₀ and AIR_g were higher among AAs vs. EAs, whereas S_I was lower, independent of age group. MLR analyses indicated that FFA₀ was best predicted by lean tissue mass (LTM), leg fat mass, ethnicity (lower in AAs), S_I, and insulin₀. FFA_nadir was best predicted by FFA₀, age group, and intra‐abdominal adipose tissue (IAAT). TIME_nadir was best predicted by leg fat mass, AIR_g, and S_I. In conclusion, indexes of insulin secretion and action were associated with FFA dynamics in healthy girls and women. Lower FFA₀ among AAs was independent of insulin₀ and S_I. Whether lower FFA₀ is associated with substrate oxidation or risk for obesity remains to be determined.     

A genome‐wide search for quantitative trait loci affecting the cortical surface area and thickness of Heschl's gyrus

Heschl's gyrus (HG) is a core region of the auditory cortex whose morphology is highly variable across individuals. This variability has been linked to sound perception ability in both speech and music domains. Previous studies show that variations in morphological features of HG, such as cortical surface area and thickness, are heritable. To identify genetic variants that affect HG morphology, we conducted a genome‐wide association scan (GWAS) meta‐analysis in 3054 healthy individuals using HG surface area and thickness as quantitative traits. None of the single nucleotide polymorphisms (SNPs) showed association P values that would survive correction for multiple testing over the genome. The most significant association was found between right HG area and SNP rs72932726 close to gene DCBLD2 (3q12.1; P = 2.77 × 10^?7). This SNP was also associated with other regions involved in speech processing. The SNP rs333332 within gene KALRN (3q21.2; P = 2.27 × 10^?6) and rs143000161 near gene COBLL1 (2q24.3; P = 2.40 × 10^?6) were associated with the area and thickness of left HG, respectively. Both genes are involved in the development of the nervous system. The SNP rs7062395 close to the X‐linked deafness gene POU3F4 was associated with right HG thickness (Xq21.1; P = 2.38 × 10^?6). This is the first molecular genetic analysis of variability in HG morphology.     

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

Background

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: P _adj = 0.005; Obese: P _adj = 0.018; Overweight+Obese: P _adj = 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: P _FET,1df = 0.026; alcohol dependence: P _FET,1df = 0.027; cocaine dependence: P _FET,1df = 0.007; marijuana dependence: P _FET,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.

Conclusion

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.     

Genome-wide association study of african and European americans implicates multiple shared and ethnic specific Loci in sarcoidosis susceptibility

Sarcoidosis is a systemic inflammatory disease characterized by the formation of granulomas in affected organs. Genome-wide association studies (GWASs) of this disease have been conducted only in European population. We present the first sarcoidosis GWAS in African Americans (AAs, 818 cases and 1,088 related controls) followed by replication in independent sets of AAs (455 cases and 557 controls) and European Americans (EAs, 442 cases and 2,284 controls). We evaluated >6 million SNPs either genotyped using the Illumina Omni1-Quad array or imputed from the 1000 Genomes Project data. We identified a novel sarcoidosis-associated locus, NOTCH4, that reached genome-wide significance in the combined AA samples (rs715299, P _AA-meta = 6.51×10⁻¹⁰) and demonstrated the independence of this locus from others in the MHC region in the same sample. We replicated previous European GWAS associations within HLA-DRA, HLA-DRB5, HLA-DRB1, BTNL2, and ANXA11 in both our AA and EA datasets. We also confirmed significant associations to the previously reported HLA-C and HLA-B regions in the EA but not AA samples. We further identified suggestive associations with several other genes previously reported in lung or inflammatory diseases.