Genome-wide association study identifies a locus at 7p15.2 associated with endometriosis

Endometriosis is a common gynecological disease associated with pelvic pain and subfertility. We conducted a genome-wide association study (GWAS) in 3,194 individuals with surgically confirmed endometriosis (cases) and 7,060 controls from Australia and the UK. Polygenic predictive modeling showed significantly increased genetic loading among 1,364 cases with moderate to severe endometriosis. The strongest association signal was on 7p15.2 (rs12700667) for 'all' endometriosis (P = 2.6 x 10(-)(7), odds ratio (OR) = 1.22, 95% CI 1.13-1.32) and for moderate to severe disease (P = 1.5 x 10(-)(9), OR = 1.38, 95% CI 1.24-1.53). We replicated rs12700667 in an independent cohort from the United States of 2,392 self-reported, surgically confirmed endometriosis cases and 2,271 controls (P = 1.2 x 10(-)(3), OR = 1.17, 95% CI 1.06-1.28), resulting in a genome-wide significant P value of 1.4 x 10(-)(9) (OR = 1.20, 95% CI 1.13-1.27) for 'all' endometriosis in our combined datasets of 5,586 cases and 9,331 controls. rs12700667 is located in an intergenic region upstream of the plausible candidate genes NFE2L3 and HOXA10.

The genetic association between personality and major depression or bipolar disorder. A polygenic score analysis using genome-wide association data

The relationship between major depressive disorder (MDD) and bipolar disorder (BD) remains controversial. Previous research has reported differences and similarities in risk factors for MDD and BD, such as predisposing personality traits. For example, high neuroticism is related to both disorders, whereas openness to experience is specific for BD. This study examined the genetic association between personality and MDD and BD by applying polygenic scores for neuroticism, extraversion, openness to experience, agreeableness and conscientiousness to both disorders. Polygenic scores reflect the weighted sum of multiple single-nucleotide polymorphism alleles associated with the trait for an individual and were based on a meta-analysis of genome-wide association studies for personality traits including 13,835 subjects. Polygenic scores were tested for MDD in the combined Genetic Association Information Network (GAIN-MDD) and MDD2000+ samples (N=8921) and for BD in the combined Systematic Treatment Enhancement Program for Bipolar Disorder and Wellcome Trust Case-Control Consortium samples (N=6329) using logistic regression analyses. At the phenotypic level, personality dimensions were associated with MDD and BD. Polygenic neuroticism scores were significantly positively associated with MDD, whereas polygenic extraversion scores were significantly positively associated with BD. The explained variance of MDD and BD, approximately 0.1%, was highly comparable to the variance explained by the polygenic personality scores in the corresponding personality traits themselves (between 0.1 and 0.4%). This indicates that the proportions of variance explained in mood disorders are at the upper limit of what could have been expected. This study suggests shared genetic risk factors for neuroticism and MDD on the one hand and for extraversion and BD on the other.

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.

Genome-wide association study identifies a new melanoma susceptibility locus at 1q21.3

We performed a genome-wide association study of melanoma in a discovery cohort of 2,168 Australian individuals with melanoma and 4,387 control individuals. In this discovery phase, we confirm several previously characterized melanoma-associated loci at MC1R, ASIP and MTAP-CDKN2A. We selected variants at nine loci for replication in three independent case-control studies (Europe: 2,804 subjects with melanoma, 7,618 control subjects; United States 1: 1,804 subjects with melanoma, 1,026 control subjects; United States 2: 585 subjects with melanoma, 6,500 control subjects). The combined meta-analysis of all case-control studies identified a new susceptibility locus at 1q21.3 (rs7412746, P = 9.0 x 10(-11), OR in combined replication cohorts of 0.89 (95% CI 0.85-0.95)). We also show evidence suggesting that melanoma associates with 1q42.12 (rs3219090, P = 9.3 x 10(-8)). The associated variants at the 1q21.3 locus span a region with ten genes, and plausible candidate genes for melanoma susceptibility include ARNT and SETDB1. Variants at the 1q21.3 locus do not seem to be associated with human pigmentation or measures of nevus density.

Variation in BMPR1B, TGFRB1 and BMPR2 and control of dizygotic twinning

Genes in the TGF9 signaling pathway play important roles in the regulation of ovarian follicle growth and ovulation rate. Mutations in three genes in this pathway, growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and the bone morphogenetic protein receptor B 1 (BMPRB1), influence dizygotic (DZ) twinning rates in sheep. To date, only variants in GDF9 and BMP15, but not their receptors transforming growth factor ss receptor 1 (TGFBR1), bone morphogenetic protein receptor 2 (BMPR2) and BMPR1B, have been investigated with respect to their roles in human DZ twinning. We screened for rare and novel variants in TGFBR1, BMPR2 and BMPR1B in mothers of dizygotic twins (MODZT) from twin-dense families, and assessed association between genotyped and imputed variants and DZ twinning in another large sample of MODZT. Three novel variants were found: a deep intronic variant in BMPR2, and one intronic and one non-synonymous exonic variant in BMPRB1 which would result in the replacement of glutamine by glutamic acid at amino acid position 294 (p.Gln294Glu). None of these variants were predicted to have major impacts on gene function. However, the p.Gln294Glu variant changes the same amino acid as a sheep BMPR1B functional variant and may have functional consequences. Six BMPR1B variants were marginally associated with DZ twinning in the larger case-control sample, but these were no longer significant once multiple testing was taken into account. Our results suggest that variation in the TGF9 signaling pathway type II receptors has limited effects on DZ twinning rates in humans.

Meta-analysis of genome-wide association for migraine in six population-based European cohorts

Migraine is a common neurological disorder with a genetically complex background. This paper describes a meta-analysis of genome-wide association (GWA) studies on migraine, performed by the Dutch-Icelandic migraine genetics (DICE) consortium, which brings together six population-based European migraine cohorts with a total sample size of 10,980 individuals (2446 cases and 8534 controls). A total of 32 SNPs showed marginal evidence for association at a P-value<10(-5). The best result was obtained for SNP rs9908234, which had a P-value of 8.00 x 10(-8). This top SNP is located in the nerve growth factor receptor (NGFR) gene. However, this SNP did not replicate in three cohorts from the Netherlands and Australia. Of the other 31 SNPs, 18 SNPs were tested in two replication cohorts, but none replicated. In addition, we explored previously identified candidate genes in the meta-analysis data set. This revealed a modest gene-based significant association between migraine and the metadherin (MTDH) gene, previously identified in the first clinic-based GWA study (GWAS) for migraine (Bonferroni-corrected gene-based P-value=0.026). This finding is consistent with the involvement of the glutamate pathway in migraine. Additional research is necessary to further confirm the involvement of glutamate.

A quantitative-trait genome-wide association study of alcoholism risk in the community: findings and implications

BACKGROUND Given moderately strong genetic contributions to variation in alcoholism and heaviness of drinking (50% to 60% heritability) with high correlation of genetic influences, we have conducted a quantitative trait genome-wide association study (GWAS) for phenotypes related to alcohol use and dependence. METHODS Diagnostic interview and blood/buccal samples were obtained from sibships ascertained through the Australian Twin Registry. Genome-wide single nucleotide polymorphism (SNP) genotyping was performed with 8754 individuals (2062 alcohol-dependent cases) selected for informativeness for alcohol use disorder and associated quantitative traits. Family-based association tests were performed for alcohol dependence, dependence factor score, and heaviness of drinking factor score, with confirmatory case-population control comparisons using an unassessed population control series of 3393 Australians with genome-wide SNP data. RESULTS No findings reached genome-wide significance (p = 8.4 x 10(-8) for this study), with lowest p value for primary phenotypes of 1.2 x 10(-7). Convergent findings for quantitative consumption and diagnostic and quantitative dependence measures suggest possible roles for a transmembrane protein gene (TMEM108) and for ANKS1A. The major finding, however, was small effect sizes estimated for individual SNPs, suggesting that hundreds of genetic variants make modest contributions (1/4% of variance or less) to alcohol dependence risk. CONCLUSIONS We conclude that: - 1) meta-analyses of consumption data may contribute usefully to gene discovery; - 2) translation of human alcoholism GWAS results to drug discovery or clinically useful prediction of risk will be challenging, and; - 3) through accumulation across studies, GWAS data may become valuable for improved genetic risk differentiation in research in biological psychiatry (e.g., prospective high-risk or resilience studies).

Association between ORMDL3, IL1RL1 and a deletion on chromosome 17q21 with asthma risk in Australia

Genome-wide association studies followed by replication provide a powerful approach to map genetic risk factors for asthma. We sought to search for new variants associated with asthma and attempt to replicate the association with four loci reported previously (ORMDL3, PDE4D, DENND1B and IL1RL1). Genome-wide association analyses of individual single nucleotide polymorphisms (SNPs), rare copy number variants (CNVs) and overall CNV burden were carried out in 986 asthma cases and 1846 asthma-free controls from Australia. The most-associated locus in the SNP analysis was ORMDL3 (rs6503525, P = 4.8 x 10(-)(7)). Five other loci were associated with P < 10(-)(5), most notably the chemokine CXC motif ligand 14 (CXCL14) gene (rs31263, P = 7.8 x 10(-)(6)). We found no evidence for association with the specific risk variants reported recently for PDE4D, DENND1B and ILR1L1. However, a variant in IL1RL1 that is in low linkage disequilibrium with that reported previously was associated with asthma risk after accounting for all variants tested (rs10197862, gene wide P = 0.01). This association replicated convincingly in an independent cohort (P = 2.4 x 10(-)(4)). A 300-kb deletion on chromosome 17q21 was associated with asthma risk, but this did not reach experiment-wide significance. Asthma cases and controls had comparable CNV rates, length and number of genes affected by deletions or duplications. In conclusion, we confirm the association between asthma risk and variants in ORMDL3 and identify a novel risk variant in IL1RL1. Follow-up of the 17q21 deletion in larger cohorts is warranted.

GWAS of butyrylcholinesterase activity identifies four novel loci, independent effects within BCHE and secondary associations with metabolic risk factors

Serum butyrylcholinesterase (BCHE) activity is associated with obesity, blood pressure and biomarkers of cardiovascular and diabetes risk. We have conducted a genome-wide association scan to discover genetic variants affecting BCHE activity, and to clarify whether the associations between BCHE activity and cardiometabolic risk factors are caused by variation in BCHE or whether BCHE variation is secondary to the metabolic abnormalities. We measured serum BCHE in adolescents and adults from three cohorts of Australian twin and family studies. The genotypes from approximately 2.4 million single-nucleotide polymorphisms (SNPs) were available in 8791 participants with BCHE measurements. We detected significant associations with BCHE activity at three independent groups of SNPs at the BCHE locus (P = 5.8 x 10(-262), 7.8 x 10(-47), 2.9 x 10(-12)) and at four other loci: RNPEP (P = 9.4 x 10(-16)), RAPH1-ABI2 (P = 4.1 x 10(-18)), UGT1A1 (P = 4.0 x 10(-8)) and an intergenic region on chromosome 8 (P = 1.4 x 10(-8)). These loci affecting BCHE activity were not associated with metabolic risk factors. On the other hand, SNPs in genes previously associated with metabolic risk had effects on BCHE activity more often than can be explained by chance. In particular, SNPs within FTO and GCKR were associated with BCHE activity, but their effects were partly mediated by body mass index and triglycerides, respectively. We conclude that variation in BCHE activity is due to multiple variants across the spectrum from uncommon/large effect to common/small effect, and partly results from (rather than causes) metabolic abnormalities.

Common SNPs explain a large proportion of the heritability for human height

SNPs discovered by genome-wide association studies (GWASs) account for only a small fraction of the genetic variation of complex traits in human populations. Where is the remaining heritability? We estimated the proportion of variance for human height explained by 294,831 SNPs genotyped on 3,925 unrelated individuals using a linear model analysis, and validated the estimation method with simulations based on the observed genotype data. We show that 45% of variance can be explained by considering all SNPs simultaneously. Thus, most of the heritability is not missing but has not previously been detected because the individual effects are too small to pass stringent significance tests. We provide evidence that the remaining heritability is due to incomplete linkage disequilibrium between causal variants and genotyped SNPs, exacerbated by causal variants having lower minor allele frequency than the SNPs explored to date.

A genome-wide association study of Cloninger's temperament scales: implications for the evolutionary genetics of personality

Variation in personality traits is 30-60% attributed to genetic influences. Attempts to unravel these genetic influences at the molecular level have, so far, been inconclusive. We performed the first genome-wide association study of Cloninger's temperament scales in a sample of 5117 individuals, in order to identify common genetic variants underlying variation in personality. Participants' scores on Harm Avoidance, Novelty Seeking, Reward Dependence, and Persistence were tested for association with 1,252,387 genetic markers. We also performed gene-based association tests and biological pathway analyses. No genetic variants that significantly contribute to personality variation were identified, while our sample provides over 90% power to detect variants that explain only 1% of the trait variance. This indicates that individual common genetic variants of this size or greater do not contribute to personality trait variation, which has important implications regarding the genetic architecture of personality and the evolutionary mechanisms by which heritable variation is maintained.

Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and approximately 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 x 10(-)(8)), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation.

A genome wide linkage scan for dizygotic twinning in 525 families of mothers of dizygotic twins

BACKGROUND: The tendency to conceive dizygotic (DZ) twins is a complex trait influenced by genetic and environmental factors. To search for new candidate loci for twinning, we conducted a genome-wide linkage scan in 525 families using microsatellite and single nucleotide polymorphism marker panels. METHODS AND RESULTS: Non-parametric linkage analyses, including 523 families containing a total of 1115 mothers of DZ twins (MODZT) from Australia and New Zealand (ANZ) and The Netherlands (NL), produced four linkage peaks above the threshold for suggestive linkage, including a highly suggestive peak at the extreme telomeric end of chromosome 6 with an exponential logarithm of odds \[(exp)LOD] score of 2.813 (P = 0.0002). Since the DZ twinning rate increases steeply with maternal age independent of genetic effects, we also investigated linkage including only families where at least one MODZT gave birth to her first set of twins before the age of 30. These analyses produced a maximum expLOD score of 2.718 (P = 0.0002), largely due to linkage signal from the ANZ cohort, however, ordered subset analyses indicated this result is most likely a chance finding in the combined dataset. Linkage analyses were also performed for two large DZ twinning families from the USA, one of which produced a peak on chromosome 2 in the region of two potential candidate genes. Sequencing of FSHR and FIGLA, along with INHBB in MODZTs from two large NL families with family specific linkage peaks directly over this gene, revealed a potentially functional variant in the 5' untranslated region of FSHR that segregated with the DZ twinning phenotype in the Utah family. CONCLUSION: Our data provide further evidence for complex inheritance of familial DZ twinning.

A variant in LIN28B is associated with 2D:4D finger-length ratio, a putative retrospective biomarker of prenatal testosterone exposure

The ratio of the lengths of an individual's second to fourth digit (2D:4D) is commonly used as a noninvasive retrospective biomarker for prenatal androgen exposure. In order to identify the genetic determinants of 2D:4D, we applied a genome-wide association approach to 1507 11-year-old children from the Avon Longitudinal Study of Parents and Children (ALSPAC) in whom 2D:4D ratio had been measured, as well as a sample of 1382 12- to 16-year-olds from the Brisbane Adolescent Twin Study. A meta-analysis of the two scans identified a single variant in the LIN28B gene that was strongly associated with 2D:4D (rs314277: p = 4.1 x 10(-8)) and was subsequently independently replicated in an additional 3659 children from the ALSPAC cohort (p = 1.53 x 10(-6)). The minor allele of the rs314277 variant has previously been linked to increased height and delayed age at menarche, but in our study it was associated with increased 2D:4D in the direction opposite to that of previous reports on the correlation between 2D:4D and age at menarche. Our findings call into question the validity of 2D:4D as a simplistic retrospective biomarker for prenatal testosterone exposure.

A versatile gene-based test for genome-wide association studies

We have derived a versatile gene-based test for genome-wide association studies (GWAS). Our approach, called VEGAS (versatile gene-based association study), is applicable to all GWAS designs, including family-based GWAS, meta-analyses of GWAS on the basis of summary data, and DNA-pooling-based GWAS, where existing approaches based on permutation are not possible, as well as singleton data, where they are. The test incorporates information from a full set of markers (or a defined subset) within a gene and accounts for linkage disequilibrium between markers by using simulations from the multivariate normal distribution. We show that for an association study using singletons, our approach produces results equivalent to those obtained via permutation in a fraction of the computation time. We demonstrate proof-of-principle by using the gene-based test to replicate several genes known to be associated on the basis of results from a family-based GWAS for height in 11,536 individuals and a DNA-pooling-based GWAS for melanoma in approximately 1300 cases and controls. Our method has the potential to identify novel associated genes; provide a basis for selecting SNPs for replication; and be directly used in network (pathway) approaches that require per-gene association test statistics. We have implemented the approach in both an easy-to-use web interface, which only requires the uploading of markers with their association p-values, and a separate downloadable application.