Identification of a Systemic Lupus Erythematosus Risk Locus Spanning ATG16L2, FCHSD2, and P2RY2 in Koreans.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a chronic autoimmune disorder whose etiology is incompletely understood, but likely involves environmental triggers in genetically susceptible individuals. Using an unbiased genome-wide association (GWA) scan and replication analysis, we sought to identify the genetic loci associated with SLE in a Korean population. METHODS: A total of 1,174 SLE cases and 4,246 population controls from Korea were genotyped and analyzed with a GWA scan to identify single-nucleotide polymorphisms (SNPs) significantly associated with SLE, after strict quality control measures were applied. For select variants, replication of SLE risk loci was tested in an independent data set of 1,416 SLE cases and 1,145 population controls from Korea and China. RESULTS: Eleven regions outside the HLA exceeded the genome-wide significance level (P = 5 × 10(-8) ). A novel SNP-SLE association was identified between FCHSD2 and P2RY2, peaking at rs11235667 (P = 1.03 × 10(-8) , odds ratio [OR] 0.59) on a 33-kb haplotype upstream of ATG16L2. In the independent replication data set, the SNP rs11235667 continued to show a significant association with SLE (replication meta-analysis P = 0.001, overall meta-analysis P = 6.67 × 10(-11) ; OR 0.63). Within the HLA region, the SNP-SLE association peaked in the class II region at rs116727542, with multiple independent effects observed in this region. Classic HLA allele imputation analysis identified HLA-DRB1*1501 and HLA-DQB1*0602, each highly correlated with one another, as most strongly associated with SLE. Ten previously established SLE risk loci were replicated: STAT1-STAT4, TNFSF4, TNFAIP3, IKZF1, HIP1, IRF5, BLK, WDFY4, ETS1, and IRAK1-MECP2. Of these loci, previously unreported, independent second risk effects of SNPs in TNFAIP3 and TNFSF4, as well as differences in the association with a putative causal variant in the WDFY4 region, were identified. CONCLUSION: Further studies are needed to identify true SLE risk effects in other loci suggestive of a significant association, and to identify the causal variants in the regions of ATG16L2, FCHSD2, and P2RY2.

MicroRNA-650 in a copy number-variable region regulates the production of interleukin 6 in human osteosarcoma cells.

Copy number variation is a well-known genetic variation. microRNAs (miRNAs/miRs) are non-coding RNAs that mediate gene expression by regulating target mRNAs. In the present study, copy number deletions encompassing miRNA coding regions were investigated to determine the association between the deletion of miRNA and its phenotypic effects. A total of 38 human miRNAs in copy number variants were identified and miR-650, which is functional in the human osteosarcoma MG-63 cell line, was selected. Overexpression of miR-650 decreased the expression of inhibitor of growth family member 4 (ING4) in the MG-63 cells and increased interleukin (IL)6 transcription, as well as IL6 secretion in IL1B-stimulated cells. Furthermore, miR-650 downregulated the amount of nuclear factor of κ light polypeptide gene enhancer in B cells inhibitor α and increased the transcriptional activity of nuclear factor (NF)κB. Downregulation of ING4 also increased the production of IL6, similar to miR-650 overexpression. Taken together, these data indicate that miR-650 plays a significant role in the production of IL6 by regulating ING4 expression and NFκB signaling in IL1B-stimulated MG-63 osteosarcoma cells.

Prediction of breast cancer survival using clinical and genetic markers by tumor subtypes.

PURPOSE: To identify the genetic variants associated with breast cancer survival, a genome-wide association study (GWAS) was conducted of Korean breast cancer patients. METHODS: From the Seoul Breast Cancer Study (SEBCS), 3,226 patients with breast cancer (1,732 in the discovery and 1,494 in the replication set) were included in a two-stage GWAS on disease-free survival (DFS) by tumor subtypes based on hormone receptor (HR) and human epidermal growth factor receptor 2 (HER2). The associations of the re-classified combined prognostic markers through recursive partitioning analysis (RPA) of DFS for breast cancer were assessed with the Cox proportional hazard model. The prognostic predictive values of the clinical and genetic models were evaluated by Harrell's C. RESULTS: In the two-stage GWAS stratified by tumor subtypes, rs166870 and rs10825036 were consistently associated with DFS in the HR+ HER2- and HR- HER2- breast cancer subtypes, respectively (Prs166870 = 2.88 × 10(-7) and Prs10825036 = 3.54 × 10(-7) in the combined set). When patients were classified by the RPA in each subtype, genetic factors contributed significantly to differentiating the high risk group associated with DFS inbreast cancer, specifically the HR+ HER2- (P discovery=1.18 × 10(-8) and P replication = 2.08 × 10(-5)) and HR- HRE2- subtypes (P discovery = 2.35 × 10(-4) and P replication = 2.60 × 10(-2)). The inclusion of the SNPs tended to improve the performance of the prognostic models consisting of age, TNM stage and tumor subtypes based on ER, PR, and HER2 status. CONCLUSION: Combined prognostic markers that include clinical and genetic factors by tumor subtypes could improve the prediction of survival in breast cancer.

Identification of FAM13A gene associated with the ratio of FEV1 to FVC in Korean population by genome-wide association studies including gene-environment interactions.

Chronic obstructive pulmonary disease (COPD) is a complex, multifactorial disease. Although smoking is a main risk factor for obstructive impairment, not all smokers develop this critical disease. We conducted a genome-wide association study to identify the association between genetic variants and pulmonary function and also examined how these variants relate to lung impairment in accordance with smoking behaviors. Using two community-based cohorts, the Ansan cohort (n=4319) and the Ansung cohort (n=3674), in the Korean Genome Epidemiology Study, we analyzed the association between genetic variants (single-nucleotide polymorphisms and haplotypes) and the ratio of FEV1 to FVC (FEV1/FVC) using multivariate linear regression models. Similar analyses were conducted after stratification by smoking status. Four genome-wide significant signals in the FAM13A gene (the strongest signal at rs2609264, P=1.76 × 10(-7) in a combined set) were associated with FEV1/FVC. For the association with ratio, the effect size in the CTGA haplotype (risk haplotype) was -0.57% (s.e., 0.11; P=2.10 × 10(-7)) as compared with the TCAG haplotype (reference haplotype) in a combined set. There was also a significant interaction of FAM13A haplotypes with heavy smoking on FEV1/FVC (P for interaction=0.028). We confirmed the previously reported association of FAM13A in 4q22.1 with pulmonary function. The FAM13A haplotypes also interacted with heavy smoking to affect the risk of reduced pulmonary function.

Genome-wide copy number variation study reveals KCNIP1 as a modulator of insulin secretion.

Copy number variations (CNVs) have emerged as another important genetic marker in addition to SNP for understanding etiology of complex diseases. In light of this, we performed a genome-wide CNV study to identify type 2 diabetes (T2D)-associated CNV using an array comparative genomic hybridization from 3180 subjects for T2D cases (n=863) and controls (n=2,317). Thus, five CNV regions having a p-value threshold ≤0.05 were identified and evaluated by validation with quantitative PCR and comparison with previously reported CNV regions in the Database of Genomic Variants. Furthermore, we performed a functional experiment to assess the biological significance of a gene encompassing a CNV region. The inhibition of KCNIP1 led to increased insulin secretion in a glucose-dependent manner, but had no effect on insulin gene transcription as well as cell apoptosis. Taken together, these data indicate that KCNIP1 from CNV study might function as a T2D-susceptibility gene whose dysregulation alters insulin production.

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density.

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10(-8)) level: 14q24.2 (rs227425, P-value 3.98 × 10(-13), SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10(-9), CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n = 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis.

Effect of glucose ingestion in plasma markers of inflammation and oxidative stress: analysis of 16 plasma markers from oral glucose tolerance test samples of normal and diabetic patients.

Sixteen plasma markers of inflammation and oxidative stress were measured during OGTT in 54 subjects. Leptin, RBP4, CRP, OPN, ANG, MDC, and MCSF concentrations significantly decreased during OGTT (P<0.05). IL6, IL8, and MCP3 concentrations significantly increased during OGTT (P<0.05). These results provide evidence that glucose ingestion affects systemic inflammation and oxidative stress.

Recapitulation of previous genome-wide association studies with two distinct pathophysiological entities of gastric cancer in the Korean population.

Gastric cancer (GC) is the most common malignancy. The incidence rates remain remarkably high in East Asians. Although genome-wide association studies in the Han Chinese and Japanese populations have so far yielded susceptibility loci for GC, these findings need to be validated in an independent ethnic group. To identify the potential heterogeneity by histological classified subtypes (intestinal and diffuse), we examined the previously reported associations in the Korean population. PRKAA1 at 5p13.1 was found to be more strongly associated with intestinal type (odds ratio, OR=1.39, 95% CI (confidence interval) =1.22-1.58, P=3.77 × 10(-7)) than diffuse type. In addition, PSCA at 8q23.3 was significantly replicated in diffuse type (OR=1.49, 95% CI=1.32-1.67, P=2.43 × 10(-11)) but far less significant in intestinal type. In conclusion, these findings could bring additional insights into the etiologic heterogeneity in gastric carcinogenesis mechanisms.

Meta-analysis identifies a MECOM gene as a novel predisposing factor of osteoporotic fracture.

BACKGROUND: Osteoporotic fracture (OF) as a clinical endpoint is a major complication of osteoporosis. To screen for OF susceptibility genes, we performed a genome-wide association study and carried out de novo replication analysis of an East Asian population. METHODS: Association was tested using a logistic regression analysis. A meta-analysis was performed on the combined results using effect size and standard errors estimated for each study. RESULTS: In a combined meta-analysis of a discovery cohort (288 cases and 1139 controls), three hospital based sets in replication stage I (462 cases and 1745 controls), and an independent ethnic group in replication stage II (369 cases and 560 for controls), we identified a new locus associated with OF (rs784288 in the MECOM gene) that showed genome-wide significance (p=3.59×10(-8); OR 1.39). RNA interference revealed that a MECOM knockdown suppresses osteoclastogenesis. CONCLUSIONS: Our findings provide new insights into the genetic architecture underlying OF in East Asians.

MicroRNAs in human lymphoblastoid cell lines.

Human lymphoblastoid cell lines (LCLs) are generated by EBV-mediated B-cell transformation to provide unlimited genomic resources for human genetics and immunological studies. The LCL is a good in vitro cell model for assessing population differences in the basal expression of genes and miRNAs as well as in cellular responses to various stimulators. Recently, the utility of LCLs was extended to pharmacogenomic studies to discover genetic factors underlying individual variations in response to chemicals and environmental stresses. Although LCLs represent generally lymphoid tissue-specific biological characteristics, genomic signatures of LCLs can distinguish patients with brain-related diseases and nonlymphoid tumors from normal controls. MicroRNA is known to be an epigenetic transcriptional regulator, and its expression is induced in abnormal conditions such as perturbagen-stimulated, virus-infected, or cancer cells. The epigenetic regulation of gene expression mediated by microRNA and DNA methylation is important for understanding the pathogenesis of cancers and complex diseases as well as discovering for therapeutic targets. For integrative genomic analyses, LCLs can be utilized to generate cellular phenotypes and various genomic data (e.g., SNP, CNV, transcriptome, methylome, etc.), which can be linked to clinical information of donors. Here, we discuss miRNA-mediated gene expression in LCLs and its application to disease genomics and global transcriptional regulatory machinery studies.

Molecular signatures in response to Isoliquiritigenin in lymphoblastoid cell lines.

Isoliquiritigenin (ISL) has been known to induce cell cycle arrest and apoptosis of various cancer cells. However, genetic factors regulating ISL effects remain unclear. The aim of this study was to identify the molecular signatures involved in ISL-induced cell death of EBV-transformed lymphoblastoid cell lines (LCLs) using microarray analyses. For gene expression and microRNA (miRNA) microarray experiments, each of 12 LCL strains was independently treated with ISL or DMSO as a vehicle control for a day prior to total RNA extraction. ISL treatment inhibited cell proliferation of LCLs in a dose-dependent manner. Microarray analysis showed that ISL-treated LCLs represented gene expression changes in cell cycle and p53 signaling pathway, having a potential as regulators in LCL survival and sensitivity to ISL-induced cytotoxicity. In addition, 36 miRNAs including five miRNAs with unknown functions were differentially expressed in ISL-treated LCLs. The integrative analysis of miRNA and gene expression profiles revealed 12 putative mRNA-miRNA functional pairs. Among them, miR-1207-5p and miR-575 were negatively correlated with p53 pathway- and cell cycle-associated genes, respectively. In conclusion, our study suggests that miRNAs play an important role in ISL-induced cytotoxicity in LCLs by targeting signaling pathways including p53 pathway and cell cycle.

Genome-wide association study in east Asians identifies novel susceptibility loci for breast cancer.

Genetic factors play an important role in the etiology of both sporadic and familial breast cancer. We aimed to discover novel genetic susceptibility loci for breast cancer. We conducted a four-stage genome-wide association study (GWAS) in 19,091 cases and 20,606 controls of East-Asian descent including Chinese, Korean, and Japanese women. After analyzing 690,947 SNPs in 2,918 cases and 2,324 controls, we evaluated 5,365 SNPs for replication in 3,972 cases and 3,852 controls. Ninety-four SNPs were further evaluated in 5,203 cases and 5,138 controls, and finally the top 22 SNPs were investigated in up to 17,423 additional subjects (7,489 cases and 9,934 controls). SNP rs9485372, near the TGF-ß activated kinase (TAB2) gene in chromosome 6q25.1, showed a consistent association with breast cancer risk across all four stages, with a P-value of 3.8×10(-12) in the combined analysis of all samples. Adjusted odds ratios (95% confidence intervals) were 0.89 (0.85-0.94) and 0.80 (0.75-0.86) for the A/G and A/A genotypes, respectively, compared with the genotype G/G. SNP rs9383951 (P = 1.9×10(-6) from the combined analysis of all samples), located in intron 5 of the ESR1 gene, and SNP rs7107217 (P = 4.6×10(-7)), located at 11q24.3, also showed a consistent association in each of the four stages. This study provides strong evidence for a novel breast cancer susceptibility locus represented by rs9485372, near the TAB2 gene (6q25.1), and identifies two possible susceptibility loci located in the ESR1 gene and 11q24.3, respectively.

Genetic susceptibility on CagA-interacting molecules and gene-environment interaction with phytoestrogens: a putative risk factor for gastric cancer.

OBJECTIVES: To evaluate whether genes that encode CagA-interacting molecules (SRC, PTPN11, CRK, CRKL, CSK, c-MET and GRB2) are associated with gastric cancer risk and whether an interaction between these genes and phytoestrogens modify gastric cancer risk. METHODS: In the discovery phase, 137 candidate SNPs in seven genes were analyzed in 76 incident gastric cancer cases and 322 matched controls from the Korean Multi-Center Cancer Cohort. Five significant SNPs in three genes (SRC, c-MET and CRK) were re-evaluated in 386 cases and 348 controls in the extension phase. Odds ratios (ORs) for gastric cancer risk were estimated adjusted for age, smoking, H. pylori seropositivity and CagA strain positivity. Summarized ORs in the total study population (462 cases and 670 controls) were presented using pooled- and meta-analysis. Plasma concentrations of phytoestrogens (genistein, daidzein, equol and enterolactone) were measured using the time-resolved fluoroimmunoassay. RESULTS: SRC rs6122566, rs6124914, c-MET rs41739, and CRK rs7208768 showed significant genetic effects for gastric cancer in both the pooled and meta-analysis without heterogeneity (pooled OR = 3.96 [95% CI 2.05-7.65], 1.24 [95% CI = 1.01-1.53], 1.19 [95% CI = 1.01-1.41], and 1.37 [95% CI = 1.15-1.62], respectively; meta OR = 4.59 [95% CI 2.74-7.70], 1.36 [95% CI = 1.09-1.70], 1.20 [95% CI = 1.00-1.44], and 1.32 [95% CI = 1.10-1.57], respectively). Risk allele of CRK rs7208768 had a significantly increased risk for gastric cancer at low phytoestrogen levels (p interaction<0.05). CONCLUSIONS: Our findings suggest that SRC, c-MET and CRK play a key role in gastric carcinogenesis by modulating CagA signal transductions and interaction between CRK gene and phytoestrogens modify gastric cancer risk.

A genome-wide association study identifies a breast cancer risk variant in ERBB4 at 2q34: results from the Seoul Breast Cancer Study.

INTRODUCTION: Although approximately 25 common genetic susceptibility loci have been identified to be independently associated with breast cancer risk through genome-wide association studies (GWAS), the genetic risk variants reported to date only explain a small fraction of the heritability of breast cancer. Furthermore, GWAS-identified loci were primarily identified in women of European descent. METHODS: To evaluate previously identified loci in Korean women and to identify additional novel breast cancer susceptibility variants, we conducted a three-stage GWAS that included 6,322 cases and 5,897 controls. RESULTS: In the validation study using Stage I of the 2,273 cases and 2,052 controls, seven GWAS-identified loci [5q11.2/MAP3K1 (rs889312 and rs16886165), 5p15.2/ROPN1L (rs1092913), 5q12/MRPS30 (rs7716600), 6q25.1/ESR1 (rs2046210 and rs3734802), 8q24.21 (rs1562430), 10q26.13/FGFR2 (rs10736303), and 16q12.1/TOX3 (rs4784227 and rs3803662)] were significantly associated with breast cancer risk in Korean women (Ptrend < 0.05). To identify additional genetic risk variants, we selected the most promising 17 SNPs in Stage I and replicated these SNPs in 2,052 cases and 2,169 controls (Stage II). Four SNPs were further evaluated in 1,997 cases and 1,676 controls (Stage III). SNP rs13393577 at chromosome 2q34, located in the Epidermal Growth Factor Receptor 4 (ERBB4) gene, showed a consistent association with breast cancer risk with combined odds ratios (95% CI) of 1.53 (1.37-1.70) (combined P for trend = 8.8 × 10-14). CONCLUSIONS: This study shows that seven breast cancer susceptibility loci, which were previously identified in European and/or Chinese populations, could be directly replicated in Korean women. Furthermore, this study provides strong evidence implicating rs13393577 at 2q34 as a new risk variant for breast cancer.

Recapitulation of the association of the Val66Met polymorphism of BDNF gene with BMI in Koreans.

Recent evidence suggests that brain-derived neurotrophic factor (BDNF) regulates food intake and the control of body weight. A common polymorphism in human BDNF, Val66Met (single-nucleotide polymorphism database (dbSNP) no. rs6265), impairs intracellular trafficking, resulting in the reduced secretion of BDNF. Several European studies have indicated that Val66Met is associated with BMI. In this study, we examined the association of the Val66Met polymorphism with BMI in Koreans (n = 20,270) from three independent epidemiological cohorts. All three studies observed a consistent association of this polymorphism with BMI, and their combined analysis demonstrated a robust correlation (ß = -0.17 ± 0.03 and P = 5.6 × 10(-8)). We also examined the effect of smoking on the link between Val66Met and BMI. The association of Val66Met with BMI was statistically significant only in the smoking group, reflecting a possible interaction between smoking and the BDNF polymorphism for BMI. Thus, we have confirmed BDNF as a genetic risk factor for BMI in an Asian population and hypothesize that the Val66Met mutation influences individual differences in BMI. In addition, smoking might interact with BDNF Val66Met to modulate BMI.

A network-based approach to prioritize results from genome-wide association studies.

Genome-wide association studies (GWAS) are a valuable approach to understanding the genetic basis of complex traits. One of the challenges of GWAS is the translation of genetic association results into biological hypotheses suitable for further investigation in the laboratory. To address this challenge, we introduce Network Interface Miner for Multigenic Interactions (NIMMI), a network-based method that combines GWAS data with human protein-protein interaction data (PPI). NIMMI builds biological networks weighted by connectivity, which is estimated by use of a modification of the Google PageRank algorithm. These weights are then combined with genetic association p-values derived from GWAS, producing what we call 'trait prioritized sub-networks.' As a proof of principle, NIMMI was tested on three GWAS datasets previously analyzed for height, a classical polygenic trait. Despite differences in sample size and ancestry, NIMMI captured 95% of the known height associated genes within the top 20% of ranked sub-networks, far better than what could be achieved by a single-locus approach. The top 2% of NIMMI height-prioritized sub-networks were significantly enriched for genes involved in transcription, signal transduction, transport, and gene expression, as well as nucleic acid, phosphate, protein, and zinc metabolism. All of these sub-networks were ranked near the top across all three height GWAS datasets we tested. We also tested NIMMI on a categorical phenotype, Crohn's disease. NIMMI prioritized sub-networks involved in B- and T-cell receptor, chemokine, interleukin, and other pathways consistent with the known autoimmune nature of Crohn's disease. NIMMI is a simple, user-friendly, open-source software tool that efficiently combines genetic association data with biological networks, translating GWAS findings into biological hypotheses.

Carriage of the V279F null allele within the gene encoding Lp-PLA2 is protective from coronary artery disease in South Korean males.

BACKGROUND: The Asia-specific PLA2G7 994G-T transversion leads to V279F substitution within the lipoprotein-associated phospholipase-A2 (Lp-PLA2) and to absence of enzyme activity in plasma. This variant offers a unique natural experiment to assess the role of Lp-PLA2 in the pathogenesis of coronary artery disease (CAD) in humans. Given conflicting results from mostly small studies, a large two-stage case-control study was warranted. METHODOLOGY/PRINCIPAL FINDINGS: PLA2G7 V279F genotypes were initially compared in 2890 male cases diagnosed with CAD before age 60 with 3128 male controls without CAD at age 50 and above and subsequently in a second independent male dataset of 877 CAD cases and 1230 controls. In the first dataset, the prevalence of the 279F null allele was 11.5% in cases and 12.8% in controls. After adjustment for age, body mass index, diabetes, smoking, glucose and lipid levels, the OR (95% CI) for CAD for this allele was 0.80 (0.66-0.97, p = 0.02). The results were very similar in the second dataset, despite lower power, with an allele frequency of 11.2% in cases and 12.5% in controls, leading to a combined OR of 0.80 (0.69-0.92), p = 0.002. The magnitude and direction of this genetic effect were fully consistent with large epidemiological studies on plasma Lp-PLA2 activity and CAD risk. CONCLUSIONS: Natural deficiency in Lp-PLA2 activity due to carriage of PLA2G7 279F allele protects from CAD in Korean men. These results provide evidence for a causal relationship between Lp-PLA2 and CAD, and support pharmacological inhibition of this enzyme as an innovative way to prevent CAD.

Genome-wide association studies identify genetic loci related to alcohol consumption in Korean men.

BACKGROUND: Genome-wide association (GWA) studies regarding the quantitative trait of alcohol consumption are limited. OBJECTIVE: The objective of the study was to explore genetic loci associated with the amount of alcohol consumed. DESIGN: We conducted a GWA study with discovery data on single nucleotide polymorphisms (SNPs) for 1721 Korean male drinkers aged 40-69 y who were included in an urban population-based cohort. Another sample that comprised 1113 male drinkers who were from an independent cohort enrolled in a rural area served as a resource for replication. At baseline (18 June 2001 through 29 January 2003), members of both cohorts provided information on average daily alcohol consumptions, and their DNA samples were collected for genotyping. RESULTS: We tested 315,914 SNPs of discovery data by using multivariate linear regression analysis adjusted for age and smoking, and 12 SNPs on chromosome 12q24 had genome-wide significant associations with alcohol consumption; adjusted P values by using Bonferroni correction were 1.6 × 10(-5) through 5.8 × 10(-46). We observed most SNPs in intronic regions and showed that the genes that harbor SNPs were C12orf51, CCDC63, MYL2, OAS3, CUX2, and RPH3A. In particular, signals in or near C12orf51, CCDC63, and MYL2 were successfully replicated in the test for 317,951 SNPs; rs2074356 in C12orf51 was in high linkage disequilibrium with SNPs in ALDH2, but other SNPs were not. CONCLUSIONS: In a GWA study, we identified loci and alleles highly associated with alcohol consumption. The findings suggest the need for further investigations on the genetic propensity for drinking excessive amounts of alcohol.

Associations of variants in CHRNA5/A3/B4 gene cluster with smoking behaviors in a Korean population.

Multiple genome-wide and targeted association studies reveal a significant association of variants in the CHRNA5-CHRNA3-CHRNB4 (CHRNA5/A3/B4) gene cluster on chromosome 15 with nicotine dependence. The subjects examined in most of these studies had a European origin. However, considering the distinct linkage disequilibrium patterns in European and other ethnic populations, it would be of tremendous interest to determine whether such associations could be replicated in populations of other ethnicities, such as Asians. In this study, we performed comprehensive association and interaction analyses for 32 single-nucleotide polymorphisms (SNPs) in CHRNA5/A3/B4 with smoking initiation (SI), smoking quantity (SQ), and smoking cessation (SC) in a Korean sample (N = 8,842). We found nominally significant associations of 7 SNPs with at least one smoking-related phenotype in the total sample (SI: P = 0.015 approximately 0.023; SQ: P = 0.008 approximately 0.028; SC: P = 0.018 approximately 0.047) and the male sample (SI: P = 0.001 approximately 0.023; SQ: P = 0.001 approximately 0.046; SC: P = 0.01). A spectrum of haplotypes formed by three consecutive SNPs located between rs16969948 in CHRNA5 and rs6495316 in the intergenic region downstream from the 5' end of CHRNB4 was associated with these three smoking-related phenotypes in both the total and the male sample. Notably, associations of these variants and haplotypes with SC appear to be much weaker than those with SI and SQ. In addition, we performed an interaction analysis of SNPs within the cluster using the generalized multifactor dimensionality reduction method and found a significant interaction of SNPs rs7163730 in LOC123688, rs6495308 in CHRNA3, and rs7166158, rs8043123, and rs11072793 in the intergenic region downstream from the 5' end of CHRNB4 to be influencing SI in the male sample. Considering that fewer than 5% of the female participants were smokers, we did not perform any analysis on female subjects specifically. Together, our detected associations of variants in the CHRNA5/A3/B4 cluster with SI, SQ, and SC in the Korean smoker samples provide strong evidence for the contribution of this cluster to the etiology of SI, ND, and SC in this Asian population.