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.

Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Genome-wide association study of migraine implicates a common susceptibility variant on 8q22.1

Migraine is a common episodic neurological disorder, typically presenting with recurrent attacks of severe headache and autonomic dysfunction. Apart from rare monogenic subtypes, no genetic or molecular markers for migraine have been convincingly established. We identified the minor allele of rs1835740 on chromosome 8q22.1 to be associated with migraine (P = 5.38 x 10(-)(9), odds ratio = 1.23, 95% CI 1.150-1.324) in a genome-wide association study of 2,731 migraine cases ascertained from three European headache clinics and 10,747 population-matched controls. The association was replicated in 3,202 cases and 40,062 controls for an overall meta-analysis P value of 1.69 x 10(-)(1)(1) (odds ratio = 1.18, 95% CI 1.127-1.244). rs1835740 is located between MTDH (astrocyte elevated gene 1, also known as AEG-1) and PGCP (encoding plasma glutamate carboxypeptidase). In an expression quantitative trait study in lymphoblastoid cell lines, transcript levels of the MTDH were found to have a significant correlation to rs1835740 (P = 3.96 x 10(-)(5), permuted threshold for genome-wide significance 7.7 x 10(-)(5). To our knowledge, our data establish rs1835740 as the first genetic risk factor for migraine.

Susceptibility locus on chromosome 1q23-25 for a schizophrenia subtype resembling deficit schizophrenia identified by latent class analysis

Context: Identifying susceptibility genes for schizophrenia may be complicated by phenotypic heterogeneity, with some evidence suggesting that phenotypic heterogeneity reflects genetic heterogeneity. Objective: To evaluate the heritability and conduct genetic linkage analyses of empirically derived, clinically homogeneous schizophrenia subtypes. Design: Latent class and linkage analysis. Setting: Taiwanese field research centers. Participants: The latent class analysis included 1236 Han Chinese individuals with DSM-IV schizophrenia. These individuals were members of a large affected-sibling-pair sample of schizophrenia (606 ascertained families), original linkage analyses of which detected a maximum logarithm of odds (LOD) of 1.8 (z = 2.88) on chromosome 10q22.3. Main Outcome Measures: Multipoint exponential LOD scores by latent class assignment and parametric heterogeneity LOD scores. Results: Latent class analyses identified 4 classes, with 2 demonstrating familial aggregation. The first (LC2) described a group with severe negative symptoms, disorganization, and pronounced functional impairment, resembling “deficit schizophrenia.” The second (LC3) described a group with minimal functional impairment, mild or absent negative symptoms, and low disorganization. Using the negative/deficit subtype, we detected genome-wide significant linkage to 1q23-25 (LOD = 3.78, empiric genome-wide P = .01). This region was not detected using the DSM-IV schizophrenia diagnosis, but has been strongly implicated in schizophrenia pathogenesis by previous linkage and association studies.Variants in the 1q region may specifically increase risk for a negative/deficit schizophrenia subtype. Alternatively, these results may reflect increased familiality/heritability of the negative class, the presence of multiple 1q schizophrenia risk genes, or a pleiotropic 1q risk locus or loci, with stronger genotype-phenotype correlation with negative/deficit symptoms. Using the second familial latent class, we identified nominally significant linkage to the original 10q peak region. Conclusion: Genetic analyses of heritable, homogeneous phenotypes may improve the power of linkage and association studies of schizophrenia and thus have relevance to the design and analysis of genome-wide association studies.

Replication of the association of common rs9939609 variant of FTO with increased BMI in an Australian adult twin population but no evidence for gene by environment (G x E) interaction

OBJECTIVE: To further investigate a common variant (rs9939609) in the fat mass- and obesity-associated gene (FTO), which recent genome-wide association studies have shown to be associated with body mass index (BMI) and obesity. DESIGN: We examined the effect of this FTO variant on BMI in 3353 Australian adult male and female twins. RESULTS: The minor A allele of rs9939609 was associated with an increased BMI (P=0.0007). Each additional copy of the A allele was associated with a mean BMI increase of approximately 1.04 kg/m(2) (approximately 3.71 kg). Using variance components decomposition, we estimate that this single-nucleotide polymorphism accounts for approximately 3% of the genetic variance in BMI in our sample (approximately 2% of the total variance). By comparing intrapair variances of monozygotic twins of different genotypes we were able to perform a direct test of gene by environment (G x E) interaction in both sexes and gene by parity (G x P) interaction in women, but no evidence was found for either. CONCLUSIONS: In addition to supporting earlier findings that the rs9939609 variant in the FTO gene is associated with an increased BMI, our results indicate that the associated genetic effect does not interact with environment or parity.

A genome-wide linkage scan provides evidence for both new and previously reported loci influencing common migraine

Latent class analysis was performed on migraine symptom data collected in a Dutch population sample (N = 12,210, 59% female) in order to obtain empirical groupings of individuals suffering from symptoms of migraine headache. Based on these heritable groupings (h(2) = 0.49, 95% CI: 0.41-0.57) individuals were classified as affected (migrainous headache) or unaffected. Genome-wide linkage analysis was performed using genotype data from 105 families with at least 2 affected siblings. In addition to this primary phenotype, linkage analyses were performed for the individual migraine symptoms. Significance levels, corrected for the analysis of multiple traits, were determined empirically via a novel simulation approach. Suggestive linkage for migrainous headache was found on chromosomes 1 (LOD = 1.63; pointwise P = 0.0031), 13 (LOD = 1.63; P = 0.0031), and 20 (LOD = 1.85; P = 0.0018). Interestingly, the chromosome 1 peak was located close to the ATP1A2 gene, associated with familial hemiplegic migraine type 2 (FHM2). Individual symptom analysis produced a LOD score of 1.97 (P = 0.0013) on chromosome 5 (photo/phonophobia), a LOD score of 2.13 (P = 0.0009) on chromosome 10 (moderate/severe pain intensity) and a near significant LOD score of 3.31 (P = 0.00005) on chromosome 13 (pulsating headache). These peaks were all located near regions previously reported in migraine linkage studies. Our results provide important replication and support for the presence of migraine susceptibility genes within these regions, and further support the utility of an LCA-based phenotyping approach and analysis of individual symptoms in migraine genetic research. Additionally, our novel "2-step" analysis and simulation approach provides a powerful means to investigate linkage to individual trait components.

Susceptibility variants for male-pattern baldness on chromosome 20p11

We carried out a genome-wide association study in 296 individuals with male-pattern baldness (androgenetic alopecia) and 347 controls. We then investigated the 30 best SNPs in an independent replication sample and found highly significant association for five SNPs on chromosome 20p11 (rs2180439 combined P = 2.7 x 10(-15)). No interaction was detected with the X-chromosomal androgen receptor locus, suggesting that the 20p11 locus has a role in a yet-to-be-identified androgen-independent pathway.

Significant evidence of one or more susceptibility loci for endometriosis with near-Mendelian inheritance on chromosome 7p13-15

BACKGROUND: Endometriosis is a common disease with a heritable component. The collaborative International Endogene Study consists of two data sets (Oxford and Australia) comprising 1176 families with multiple affected. The aim was to investigate whether the apparent concentration of cases in a proportion of families could be explained by one or more rare variants with (near-)Mendelian autosomal inheritance. METHODS AND RESULTS: Linkage analyses (aimed at finding chromosomal regions harbouring disease-predisposing genes) were conducted in families with three or more affected (Oxford: n = 52; Australia: n = 196). In the Oxford data set, a non-parametric linkage score (Kong & Cox (K&C) Log of ODds (LOD)) of 3.52 was observed on chromosome 7p (genome-wide significance P = 0.011). A parametric MOD score (equal to maximum LOD maximized over 357 possible inheritance models) of 3.89 was found at 65.72 cM (D7S510) for a dominant model with reduced penetrance. After including the Australian data set, the non-parametric K&C LOD of the combined data set was 1.46 at 57.3 cM; the parametric analysis found an MOD score of 3.30 at D7S484 (empirical significance: P = 0.035) for a recessive model with high penetrance. Critical recombinant analysis narrowed the probable region of linkage down to overlapping 6.4 Mb and 11 Mb intervals containing 48 and 96 genes, respectively. CONCLUSIONS: This is the first report to suggest that there may be one or more high-penetrance susceptibility loci for endometriosis with (near-)Mendelian inheritance.

KRAS variation and risk of endometriosis

Endometriosis is a common gynaecological disease with symptoms of pelvic pain and infertility which affects 7-10% of women in their reproductive years. Activation of an oncogenic allele of Kirsten rat sarcoma viral oncogene homologue (KRAS) in the reproductive tract of mice resulted in the development of endometriosis. We hypothesized that variation in KRAS may influence risk of endometriosis in humans. Thirty tagSNPs spanning a region of 60.7 kb across the KRAS locus were genotyped using iPLEX chemistry on a MALDI-TOF MassARRAY platform in 959 endometriosis cases and 959 unrelated controls, and data were analysed for association with endometriosis. Genotypes were obtained for most individuals with a mean completion rate of 99.1%. We identified six haplotype blocks across the KRAS locus in our sample. There were no significant differences between cases and controls in the frequencies of individual single-nucleotide polymorphisms (SNPs) or haplotypes. We also developed a rapid method to screen for 11 common KRAS and BRAF mutations on the Sequenom MassARRAY system. The assay detected all mutations previously identified by direct sequencing in a panel of positive controls. No germline variants for KRAS or BRAF were detected. Our results demonstrate that any risk of endometriosis in women because of common variation in KRAS must be very small.

Association study of the dystrobrevin-binding gene with schizophrenia in Australian and Indian samples

Numerous studies have reported association between variants in the dystrobrevin binding protein 1 (dysbindin) gene (DTNBP1) and schizophrenia. However, the pattern of results is complex and to date, no specific risk marker or haplotype has been consistently identified. The number of single nucleotide polymorphisms (SNPs) tested in these studies has ranged from 5 to 20. We attempted to replicate previous findings by testing 16 SNPs in samples of 41 Australian pedigrees, 194 Australian cases and 180 controls, and 197 Indian pedigrees. No globally significant evidence for association was observed in any sample, despite power calculations indicating sufficient power to replicate several previous findings. Possible explanations for our results include sample differences in background linkage disequilibrium and/or risk allele effect size, the presence of multiple risk alleles upon different haplotypes, or the presence of a single risk allele upon multiple haplotypes. Some previous associations may also represent false positives. Examination of Caucasian HapMap phase II genotype data spanning the DTNBP1 region indicates upwards of 40 SNPs are required to satisfactorily assess all nonredundant variation within DTNBP1 and its potential regulatory regions for association with schizophrenia. More comprehensive studies in multiple samples will be required to determine whether specific DTNBP1 variants function as risk factors for schizophrenia.

Genomewide linkage study in 1,176 affected sister pair families identifies a significant susceptibility locus for endometriosis on chromosome 10q26

Endometriosis is a common gynecological disease that affects up to 10% of women in their reproductive years. It causes pelvic pain, severe dysmenorrhea, and subfertility. The disease is defined as the presence of tissue resembling endometrium in sites outside the uterus. Its cause remains uncertain despite >50 years of hypothesis-driven research, and thus the therapeutic options are limited. Disease predisposition is inherited as a complex genetic trait, which provides an alternative route to understanding the disease. We seek to identify susceptibility loci, using a positional-cloning approach that starts with linkage analysis to identify genomic regions likely to harbor these genes. We conducted a linkage study of 1,176 families (931 from an Australian group and 245 from a U.K. group), each with at least two members--mainly affected sister pairs--with surgically diagnosed disease. We have identified a region of significant linkage on chromosome 10q26 (maximum LOD score [MLS] of 3.09; genomewide P = .047) and another region of suggestive linkage on chromosome 20p13 (MLS = 2.09). Minor peaks (with MLS > 1.0) were found on chromosomes 2, 6, 7, 8, 12, 14, 15, and 17. This is the first report of linkage to a major locus for endometriosis. The findings will facilitate discovery of novel positional genetic variants that influence the risk of developing this debilitating disease. Greater understanding of the aberrant cellular and molecular mechanisms involved in the etiology and pathophysiology of endometriosis should lead to better diagnostic methods and targeted treatments.

Separate and interacting effects within the catechol-O-methyltransferase (COMT) are associated with schizophrenia

Several lines of evidence have implicated the catechol-O-methyltransferase (COMT) gene as a candidate for schizophrenia (SZ) susceptibility, not only because it encodes a key dopamine catabolic enzyme but also because it maps to the velocardiofacial syndrome region of chromosome 22q11 which has long been associated with SZ predisposition. The interest in COMT as a candidate SZ risk factor has led to numerous case-control and family-based studies, with the majority placing emphasis on examining a functional Val/Met polymorphism within this enzyme. Unfortunately, these studies have continually produced conflicting results. To assess the genetic contribution of other COMT variants to SZ susceptibility, we investigated three single-nucleotide polymorphisms (SNPs) (rs737865, rs4633, rs165599) in addition to the Val/Met variant (rs4680) in a highly selected sample of Australian Caucasian families containing 107 patients with SZ. The Val/Met and rs4633 variants showed nominally significant associations with SZ (P<0.05), although neither of the individual SNPs remained significant after adjusting for multiple testing (most significant P=0.1174). However, haplotype analyses showed strong evidence of an association; the most significant being the three-marker haplotype rs737865-rs4680-rs165599 (global P=0.0022), which spans more than 26 kb. Importantly, conditional analyses indicated the presence of two separate and interacting effects within this haplotype, irrespective of gender. In addition, our results indicate the Val/Met polymorphism is not disease-causing and is simply in strong linkage disequilibrium with a causative effect, which interacts with another as yet unidentified variant approximately 20 kb away. These results may help explain the inconsistent results reported on the Val/Met polymorphism and have important implications for future investigations into the role of COMT in SZ susceptibility.

Dominant negative ATM mutations in breast cancer families

BACKGROUND: The ATM gene encoding a putative protein kinase is mutated in ataxia-telangiectasia (A-T), an autosomal recessive disorder with a predisposition for cancer. Studies of A-T families suggest that female heterozygotes have an increased risk of breast cancer compared with noncarriers. However, neither linkage analyses nor mutation studies have provided supporting evidence for a role of ATM in breast cancer predisposition. Nevertheless, two recurrent ATM mutations, T7271G and IVS10-6T-->G, reportedly increase the risk of breast cancer. We examined these two ATM mutations in a population-based, case-control series of breast cancer families and multiple-case breast cancer families. METHODS: Five hundred twenty-five or 262 case patients with breast cancer and 381 or 68 control subjects, respectively, were genotyped for the T7271G and IVS10-6T-->G ATM mutations, as were index patients from 76 non-BRCA1/2 multiple-case breast cancer families. Linkage and penetrance were analyzed. ATM protein expression and kinase activity were analyzed in lymphoblastoid cell lines from mutation carriers. All statistical tests were two-sided. RESULTS: In case and control subjects unselected for family history of breast cancer, one case patient had the T7271G mutation, and none had the IVS10-6T-->G mutation. In three multiple-case families, one of these two mutations segregated with breast cancer. The estimated average penetrance of the mutations was 60% (95% confidence interval [CI] = 32% to 90%) to age 70 years, equivalent to a 15.7-fold (95% CI = 6.4-fold to 38.0-fold) increased relative risk compared with that of the general population. Expression and activity analyses of ATM in heterozygous cell lines indicated that both mutations are dominant negative. CONCLUSION: At least two ATM mutations are associated with a sufficiently high risk of breast cancer to be found in multiple-case breast cancer families. Full mutation analysis of the ATM gene in such families could help clarify the role of ATM in breast cancer susceptibility.

Further evidence for linkage of Gilles de la Tourette syndrome (GTS) susceptibility loci on chromosomes 2p11, 8q22 and 11q23-24 in South African Afrikaners

Utilizing DNA samples from 91 Afrikaner nuclear families with one or more affected children, five genomic regions on chromosomes 2p, 8q, 11q, 20q, and 21q that gave evidence for association with GTS in previous case-control association studies were investigated for linkage and association with GTS. Highly polymorphic markers with mean heterozygosity of 0.77 were typed and resulting genotypes evaluated using single marker transmission disequilibrium (TDT), single marker haplotype relative risk (HRR), and multi-marker "extended" TDT and HRR methods. Single marker TDT analysis showed evidence for linkage or association, with p-values near 0.05, for markers D2S139, GATA28F12, and D11S1377 on chromosomes 2p11, 8q22 and 11q23-24, respectively. Extended, two-locus TDT and HRR analysis provided further evidence for linkage or association on chromosome 2 with p-values of 0.007 and 0.025, and chromosome 8 with p-values of 0.059 and 0.013, respectively. These results provide important additional evidence for the location of GTS susceptibility loci.

A genomewide screen for autism susceptibility loci

We report the analysis of 335 microsatellite markers genotyped in 110 multiplex families with autism. All families include at least two "affected" siblings, at least one of whom has autism; the remaining affected sibs carry diagnoses of either Asperger syndrome or pervasive developmental disorder. Affected sib-pair analysis yielded multipoint maximum LOD scores (MLS) that reach the accepted threshold for suggestive linkage on chromosomes 5, X, and 19. Nominal evidence for linkage (point-wise P<.05) was obtained on chromosomes 2, 3, 4, 8, 10, 11, 12, 15, 16, 18, and 20, and secondary loci were found on chromosomes 5 and 19. Analysis of families sharing alleles at the putative X chromosomal linked locus and one or more other putative linked loci produced an MLS of 3.56 for the DXS470-D19S174 marker combination. In an effort to increase power to detect linkage, scan statistics were used to evaluate the significance of peak LOD scores based on statistical evidence at adjacent marker loci. This analysis yielded impressive evidence for linkage to autism and autism-spectrum disorders with significant genomewide P values <.05 for markers on chromosomes 5 and 8 and with suggestive linkage evidence for a marker on chromosome 19.