A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other

In this report, we describe a simple correction for multiple testing of single-nucleotide polymorphisms (SNPs) in linkage disequilibrium (LD) with each other, on the basis of the spectral decomposition (SpD) of matrices of pairwise LD between SNPs. This method provides a useful alternative to more computationally intensive permutation tests. Additionally, output from SNPSpD includes eigenvalues, principal-component coefficients, and factor "loadings" after varimax rotation, enabling the selection of a subset of SNPs that optimize the information in a genomic region.

Multiple sclerosis susceptibility-associated SNPs do not influence disease severity measures in a cohort of Australian MS patients

Recent association studies in multiple sclerosis (MS) have identified and replicated several single nucleotide polymorphism (SNP) susceptibility loci including CLEC16A, IL2RA, IL7R, RPL5, CD58, CD40 and chromosome 12q13–14 in addition to the well established allele HLA-DR15. There is potential that these genetic susceptibility factors could also modulate MS disease severity, as demonstrated previously for the MS risk allele HLA-DR15. We investigated this hypothesis in a cohort of 1006 well characterised MS patients from South-Eastern Australia. We tested the MS-associated SNPs for association with five measures of disease severity incorporating disability, age of onset, cognition and brain atrophy. We observed trends towards association between the RPL5 risk SNP and time between first demyelinating event and relapse, and between the CD40 risk SNP and symbol digit test score. No associations were significant after correction for multiple testing. We found no evidence for the hypothesis that these new MS disease risk-associated SNPs influence disease severity.

Genome-wide linkage analysis of multiple measures of neuroticism of 2 large cohorts from Australia and the Netherlands

CONTEXT People meeting diagnostic criteria for anxiety or depressive disorders tend to score high on the personality scale of neuroticism. Studying this personality dimension can give insights into the etiology of these important psychiatric disorders. OBJECTIVES To undertake a comprehensive genome-wide linkage study of neuroticism using large study samples that have been measured multiple times and to compare the results between countries for replication and across time within countries for consistency. DESIGN Genome-wide linkage scan. SETTING Twin individuals and their family members from Australia and the Netherlands. PARTICIPANTS Nineteen thousand six hundred thirty-five sibling pairs completed self-report questionnaires for neuroticism up to 5 times over a period of up to 22 years. Five thousand sixty-nine sibling pairs were genotyped with microsatellite markers. METHODS Nonparametric linkage analyses were conducted in MERLIN-REGRESS for the mean neuroticism scores averaged across time. Additional analyses were conducted for the time-specific measures of neuroticism from each country to investigate consistency of linkage results. RESULTS Three chromosomal regions exceeded empirically derived thresholds for suggestive linkage using mean neuroticism scores: 10p 5 Kosambi cM (cM) (Dutch study sample), 14q 103 cM (Dutch study sample), and 18q 117 cM (combined Australian and Dutch study sample), but only 14q retained significance after correction for multiple testing. These regions all showed evidence for linkage in individual time-specific measures of neuroticism and 1 (18q) showed some evidence for replication between countries. Linkage intervals for these regions all overlap with regions identified in other studies of neuroticism or related traits and/or in studies of anxiety in mice. CONCLUSIONS Our results demonstrate the value of the availability of multiple measures over time and add to the optimism reported in recent reviews for replication of linkage regions for neuroticism. These regions are likely to harbor causal variants for neuroticism and its related psychiatric disorders and can inform prioritization of results from genome-wide association studies.

Brain-derived neurotrophic factor (BDNF) gene : no major impact on antidepressant treatment response

The brain-derived neurotrophic factor (BDNF) has been suggested to play a pivotal role in the aetiology of affective disorders. In order to further clarify the impact of BDNF gene variation on major depression as well as antidepressant treatment response, association of three BDNF polymorphisms [rs7103411, Val66Met (rs6265) and rs7124442] with major depression and antidepressant treatment response was investigated in an overall sample of 268 German patients with major depression and 424 healthy controls. False discovery rate (FDR) was applied to control for multiple testing. Additionally, ten markers in BDNF were tested for association with citalopram outcome in the STAR*D sample. While BDNF was not associated with major depression as a categorical diagnosis, the BDNF rs7124442 TT genotype was significantly related to worse treatment outcome over 6 wk in major depression (p=0.01) particularly in anxious depression (p=0.003) in the German sample. However, BDNF rs7103411 and rs6265 similarly predicted worse treatment response over 6 wk in clinical subtypes of depression such as melancholic depression only (rs7103411: TT

Association study of the calcitonin gene-related polypeptide-alpha (CALCA) and the receptor activity modifying 1 (RAMP1) genes with migraine

Migraine is a common neurovascular brain disorder characterised by recurrent attacks of severe headache that may be accompanied by various neurological symptoms. Migraine is thought to result from activation of the trigeminovascular system followed by vasodilation of pain-producing intracranial blood vessels and activation of second-order sensory neurons in the trigeminal nucleus caudalis. Calcitonin gene-related peptide (CGRP) is a mediator of neurogenic inflammation and the most powerful vasodilating neuropeptide, and has been implicated in migraine pathophysiology. Consequently, genes involved in CGRP synthesis or CGRP receptor genes may play a role in migraine and/or increase susceptibility. This study investigates whether variants in the gene that encodes CGRP, calcitonin-related polypeptide alpha (CALCA) or in the gene that encodes a component of its receptor, receptor activity modifying protein 1 (RAMP1), are associated with migraine pathogenesis and susceptibility. The single nucleotide polymorphisms (SNPs) rs3781719 and rs145837941 in the CALCA gene, and rs3754701 and rs7590387 at the RAMP1 locus, were analysed in an Australian Caucasian population of migraineurs and matched controls. Although we find no significant association of any of the SNPs tested with migraine overall, we detected a nominally significant association (p = 0.031) of the RAMP1 rs3754701 variant in male migraine subjects, although this is non-significant after Bonferroni correction for multiple testing.

Association analysis of chromosome 1 migraine candidate genes

Migraine with aura (MA) is a subtype of typical migraine. Migraine with aura (MA) also encompasses a rare severe subtype Familial Hemiplegic Migraine (FHM) with several known genetic loci. The type 2 FHM (FHM-2) susceptibility locus maps to chromosome 1q23 and mutations in the ATP1A2 gene at this site have recently been implicated. We have previously provided evidence of linkage of typical migraine (predominantly MA) to microsatellite markers on chromosome 1, in the 1q31 and 1q23 regions. In this study, we have undertaken a large genomic investigation involving candidate genes that lie within the chromosome 1q23 and 1q31 regions using an association analysis approach. Methods We have genotyped a large population of case-controls (243 unrelated Caucasian migraineurs versus 243 controls) examining a set of 5 single nucleotide polymorphisms (SNPs) and the Fas Ligand dinucleotide repeat marker, located within the chromosome 1q23 and 1q31 regions. Results Several genes have been studied including membrane protein (ATP 1 subtype A4 and FasL), cytoplasmic glycoprotein (CASQ 1) genes and potassium (KCN J9 and KCN J10) and calcium (CACNA1E) channel genes in 243 migraineurs (including 85% MA and 15% of migraine without aura (MO)) and 243 matched controls. After correction for multiple testing, chi-square results showed non-significant P values (P > 0.008) across all SNPs (and a CA repeat) tested in these different genes, however results with the KCN J10 marker gave interesting results (P = 0.02) that may be worth exploring further in other populations. Conclusion These results do not show a significant role for the tested candidate gene variants and also do not support the hypothesis that a common chromosome 1 defective gene influences both FHM and the more common forms of migraine.

The GSTM1 null genotype confers an increased risk for Solar Keratosis Development in an Australian Caucasian Population

Solar keratoses affect approximately 50% of Australian Caucasians aged over 40 y. Solar keratoses can undergo malignant transformation into squamous cell carcinoma followed by possible metastasis and are risk factors for basal cell carcinoma, melanoma, and squamous cell carcinoma. The glutathione-S-transferase genes play a part in detoxification of carcinogens and mutagens, including some produced by ultraviolet radiation. This study examined the role of glutathione-S-transferase M1, T1, P1, and Z1 gene polymorphisms in susceptibility to solar keratoses development. Using DNA samples from volunteers involved in the Nambour Skin Cancer Prevention Trial, allele and genotype frequencies were determined using polymerase chain reaction and restriction enzyme digestion. No significant differences were detected in glutathione-S-transferase P1 and glutathione-S-transferase Z1 allele or genotype frequencies; however, a significant association between glutathione-S-transferase M1 genotypes and solar keratoses development was detected (p=0.003) with null individuals having an approximate 2-fold increase in risk for solar keratoses development (odds ratio: 2.1; confidence interval: 1.3-3.5) and a significantly higher increase in risk in conjunction with high outdoor exposure (odds ratio: 3.4; confidence interval: 1.9-6.3). Also, a difference in glutathione-S-transferase T1 genotype frequencies was detected (p=0.039), although considering that multiple testing was undertaken, this was found not to be significant. Fair skin and inability to tan were found to be highly significant risk factors for solar keratoses development with odds ratios of 18.5 (confidence interval: 5.7-59.9) and 7.4 (confidence interval: 2.6-21.0), respectively. Overall, glutathione-S-transferase M1 conferred a significant increase in risk of solar keratoses development, particularly in the presence of high outdoor exposure and synergistically with known phenotypic risk factors of fair skin and inability to tan.

Environmental, personal, and genetic determinants of response to vitamin D supplementation in older adults

CONTEXT AND OBJECTIVE: Suboptimal vitamin D status can be corrected by vitamin D supplementation, but individual responses to supplementation vary. We aimed to examine genetic and nongenetic determinants of change in serum 25-hydroxyvitamin D (25(OH)D) after supplementation. DESIGN AND PARTICIPANTS: We used data from a pilot randomized controlled trial in which 644 adults aged 60 to 84 years were randomly assigned to monthly doses of placebo, 30 000 IU, or 60 000 IU vitamin D3 for 12 months. Baseline characteristics were obtained from a self-administered questionnaire. Eighty-eight single-nucleotide polymorphisms (SNPs) in 41 candidate genes were genotyped using Sequenom MassArray technology. Serum 25(OH)D levels before and after the intervention were measured using the Diasorin Liaison platform immunoassay. We used linear regression models to examine associations between genetic and nongenetic factors and change in serum 25(OH)D levels. RESULTS: Supplement dose and baseline 25(OH)D level explained 24% of the variability in response to supplementation. Body mass index, self-reported health status, and ambient UV radiation made a small additional contribution. SNPs in CYP2R1, IRF4, MC1R, CYP27B1, VDR, TYRP1, MCM6, and HERC2 were associated with change in 25(OH)D level, although only CYP2R1 was significant after adjustment for multiple testing. Models including SNPs explained a similar proportion of variability in response to supplementation as models that included personal and environmental factors. CONCLUSION: Stepwise regression analyses suggest that genetic variability may be associated with response to supplementation, perhaps suggesting that some people might need higher doses to reach optimal 25(OH)D levels or that there is variability in the physiologically normal level of 25(OH)D.

Dysbindin (DTNBP1) variants are associated with hallucinations in schizophrenia

BACKGROUND: Dystrobrevin binding protein 1 (DTNBP1) is a schizophrenia susceptibility gene involved with neurotransmission regulation (especially dopamine and glutamate) and neurodevelopment. The gene is known to be associated with cognitive deficit phenotypes within schizophrenia. In our previous studies, DTNBP1 was found associated not only with schizophrenia but with other psychiatric disorders including psychotic depression, post-traumatic stress disorder, nicotine dependence and opiate dependence. These findings suggest that DNTBP1 may be involved in pathways that lead to multiple psychiatric phenotypes. In this study, we explored the association between DTNBP1 SNPs (single nucleotide polymorphisms) and multiple psychiatric phenotypes included in the Diagnostic Interview of Psychosis (DIP). METHODS: Five DTNBP1 SNPs, rs17470454, rs1997679, rs4236167, rs9370822 and rs9370823, were genotyped in 235 schizophrenia subjects screened for various phenotypes in the domains of depression, mania, hallucinations, delusions, subjective thought disorder, behaviour and affect, and speech disorder. SNP-phenotype association was determined with ANOVA under general, dominant/recessive and over-dominance models. RESULTS: Post hoc tests determined that SNP rs1997679 was associated with visual hallucination; SNP rs4236167 was associated with general auditory hallucination as well as specific features including non-verbal, abusive and third-person form auditory hallucinations; and SNP rs9370822 was associated with visual and olfactory hallucinations. SNPs that survived correction for multiple testing were rs4236167 for third-person and abusive form auditory hallucinations; and rs9370822 for olfactory hallucinations. CONCLUSION: These data suggest that DTNBP1 is likely to play a role in development of auditory related, visual and olfactory hallucinations which is consistent with evidence of DTNBP1 activity in the auditory processing regions, in visual processing and in the regulation of glutamate and dopamine activity

No genetic support for a contribution of prostaglandins to the aetiology of androgenetic alopecia

MADAM, Androgenetic alopecia (AGA) is a common age-dependent trait, characterized by a progressive loss of hair from the scalp. The hair loss may commence during puberty and up to 80% of white men experience some degree of AGA during their lifetime.1 Research has established that two essential aetiological factors for AGA are a genetic predisposition and the presence of androgens (male sex hormones).1,2 A recent meta-analysis of genome-wide association studies (GWAS) has increased the number of identified loci associated with this trait at the molecular level to a total of eight.3 However, despite these successes, a large fraction of the genetic contribution remains to be identified. One way to identify further genetic loci is to combine the resource of GWAS datasets with knowledge about specific biological factors likely to be involved in the development of disease. The focused evaluation of a limited number of candidate genes in GWAS datasets avoids the necessity for extensive correction for multiple testing, which typically limits the power for detecting genetic loci at a genome-wide level.4 Because the presence of genetic association suggests that candidate genes are likely to operate early in the causative chain of events leading to the phenotype, this approach may also function to favour biological pathways for their importance in the development of AGA.

A genome-wide association study of sleep habits and insomnia

Several aspects of sleep behavior such as timing, duration and quality have been demonstrated to be heritable. To identify common variants that influence sleep traits in the population, we conducted a genome-wide association study of six sleep phenotypes assessed by questionnaire in a sample of 2,323 individuals from the Australian Twin Registry. Genotyping was performed on the Illumina 317, 370, and 610K arrays and the SNPs in common between platforms were used to impute non-genotyped SNPs. We tested for association with more than 2,000,000 common polymorphisms across the genome. While no SNPs reached the genome-wide significance threshold, we identified a number of associations in plausible candidate genes. Most notably, a group of SNPs in the third intron of the CACNA1C gene ranked as most significant in the analysis of sleep latency (P = 1.3 x 10(-)(6)). We attempted to replicate this association in an independent sample from the Chronogen Consortium (n = 2,034), but found no evidence of association (P = 0.73). We have identified several other suggestive associations that await replication in an independent sample. We did not replicate the results from previous genome-wide analyses of self-reported sleep phenotypes after correction for multiple testing.

Association mapping

Association mapping seeks to identify marker alleles present at significantly different frequencies in cases carrying a particular disease or trait compared with controls. Genome-wide association studies are increasingly replacing candidate gene-based association studies for complex diseases, where a number of loci are likely to contribute to disease risk and the effect size of each particular risk allele is typically modest or low. Good study design is essential to the success of an association study, and factors such as the heritability of the disease under investigation, the choice of controls, statistical power, multiple testing and whether the association can be replicated need to be considered before beginning. Likewise, thorough quality control of the genotype data needs to be undertaken prior to running any association analyses. Finally, it should be kept in mind that a significant genetic association is not proof positive that a particular genetic locus causes a disease, but rather an important first step in discovering the genetic variants underlying a complex disease.

Educational attainment: A genome wide association study in 9538 Australians

BACKGROUND Correlations between Educational Attainment (EA) and measures of cognitive performance are as high as 0.8. This makes EA an attractive alternative phenotype for studies wishing to map genes affecting cognition due to the ease of collecting EA data compared to other cognitive phenotypes such as IQ. METHODOLOGY In an Australian family sample of 9538 individuals we performed a genome-wide association scan (GWAS) using the imputed genotypes of approximately 2.4 million single nucleotide polymorphisms (SNP) for a 6-point scale measure of EA. Top hits were checked for replication in an independent sample of 968 individuals. A gene-based test of association was then applied to the GWAS results. Additionally we performed prediction analyses using the GWAS results from our discovery sample to assess the percentage of EA and full scale IQ variance explained by the predicted scores. RESULTS The best SNP fell short of having a genome-wide significant p-value (p = 9.77x10(-7)). In our independent replication sample six SNPs among the top 50 hits pruned for linkage disequilibrium (r(2)<0.8) had a p-value<0.05 but only one of these SNPs survived correction for multiple testing--rs7106258 (p = 9.7*10(-4)) located in an intergenic region of chromosome 11q14.1. The gene based test results were non-significant and our prediction analyses show that the predicted scores explained little variance in EA in our replication sample. CONCLUSION While we have identified a polymorphism chromosome 11q14.1 associated with EA, further replication is warranted. Overall, the absence of genome-wide significant p-values in our large discovery sample confirmed the high polygenic architecture of EA. Only the assembly of large samples or meta-analytic efforts will be able to assess the implication of common DNA polymorphisms in the etiology of EA.