The influence of regression models on genome-wide association studies of alcohol dependence: a comparison of binary and quantitative analyses.

INTRODUCTION: Genome-wide association studies (GWAS) of alcohol dependence syndrome (ADS) offer a platform to detect genetic risk loci. However, the majority of the ADS GWAS undertaken, to date, have utilized a case-control design and have failed to identify consistently replicable loci with the exception of protective variants within the alcohol metabolizing genes, notably ADH1B. The ADS phenotype shows considerable variability which means that the use of quantitative variables as a proxy for the severity of ADS has the potential to facilitate identification of risk loci by increasing statistical power. The current study aims to examine the influences of using binary and adjusted quantitative measures of ADS on GWAS outcomes and on calculated polygenic risk scores (PRS). METHODS: A GWAS was performed in 1251 healthy controls with no history of excess alcohol use and 739 patients with ADS classified using binary DMS-IV criteria. Two additional GWAS were undertaken using a quantitative score based on DSM-IV criteria, which were applied assuming both normal and non-normal distributions of the phenotypic variables. PRS analyses were performed utilizing the data from the binary and the quantitative trait analyses. RESULTS: No associations were identified at genome-wide significance in any of the individual GWAS; results were comparable in all three. The top associated single nucleotide polymorphism was located on the alcohol dehydrogenase gene cluster on chromosome 4, consistent with previous ADS GWAS. The quantitative trait analysis adjusted for the distribution of the criterion score and the associated PRS had the smallest standard errors and thus the greatest precision. CONCLUSION: Further exploitation of the use of qualitative trait analysis in GWAS in ADS is warranted.

Genetic variation in the miR-708 gene and its binding targets in bipolar disorder.

OBJECTIVE: rs12576775 was found to be associated with bipolar disorder (BD) in a genome-wide association study (GWAS). The GWAS signal implicates genes for the microRNAs miR-708 and miR-5579 and the first exon of the Odd Oz/ten-m homolog 4 gene (ODZ4). In the present study, miR-708, its surrounding region, and its targets were analyzed for potential BD-associated functional variants. METHODS: The miR-708 gene and surrounding regions were screened for variation using high-resolution melting (HRM) analysis in 1099 cases of BD, followed by genotyping of rare variants in an enlarged sample of 2078 subjects with BD, 1303 subjects with schizophrenia, and 1355 healthy controls. Whole-genome sequencing data from 99 subjects with BD were analyzed for variation in potential miR-708 binding sites. The minor allele frequencies (MAFs) of these variants were compared with those reported in reference individuals. RESULTS: Three variants detected by HRM were selected to be genotyped. rs754333774 was detected in three cases of BD, two cases of schizophrenia, and no controls. This variant is located 260 base pairs upstream from miR-708 and may play a role in controlling the expression of the miR. Four variants were identified in miR-708 targets binding sites. The MAFs of each of these variants were similar in BD and reference samples. CONCLUSIONS: We report a single recurrent variant located near the miR-708 gene that may have a role in BD and schizophrenia susceptibility. These findings await replication in independent cohorts, as do functional analyses of the potential consequences of this variant.

Analysis of ANK3 and CACNA1C variants identified in bipolar disorder whole genome sequence data.

OBJECTIVES: Genetic markers in the genes encoding ankyrin 3 (ANK3) and the α-calcium channel subunit (CACNA1C) are associated with bipolar disorder (BP). The associated variants in the CACNA1C gene are mainly within intron 3 of the gene. ANK3 BP-associated variants are in two distinct clusters at the ends of the gene, indicating disease allele heterogeneity. METHODS: In order to screen both coding and non-coding regions to identify potential aetiological variants, we used whole-genome sequencing in 99 BP cases. Variants with markedly different allele frequencies in the BP samples and the 1,000 genomes project European data were genotyped in 1,510 BP cases and 1,095 controls. RESULTS: We found that the CACNA1C intron 3 variant, rs79398153, potentially affecting an ENCyclopedia of DNA Elements (ENCODE)-defined region, showed an association with BP (p = 0.015). We also found the ANK3 BP-associated variant rs139972937, responsible for an asparagine to serine change (p = 0.042). However, a previous study had not found support for an association between rs139972937 and BP. The variants at ANK3 and CACNA1C previously known to be associated with BP were not in linkage disequilibrium with either of the two variants that we identified and these are therefore independent of the previous haplotypes implicated by genome-wide association. CONCLUSIONS: Sequencing in additional BP samples is needed to find the molecular pathology that explains the previous association findings. If changes similar to those we have found can be shown to have an effect on the expression and function of ANK3 and CACNA1C, they might help to explain the so-called 'missing heritability' of BP.