Discovery of common variants associated with low TSH levels and thyroid cancer risk

To search for sequence variants conferring risk of nonmedullary thyroid cancer, we focused our analysis on 22 SNPs with a P < 5 × 10(-8) in a genome-wide association study on levels of thyroid stimulating hormone (TSH) in 27,758 Icelanders. Of those, rs965513 has previously been shown to associate with thyroid cancer. The remaining 21 SNPs were genotyped in 561 Icelandic individuals with thyroid cancer (cases) and up to 40,013 controls. Variants suggestively associated with thyroid cancer (P < 0.05) were genotyped in an additional 595 non-Icelandic cases and 2,604 controls. After combining the results, three variants were shown to associate with thyroid cancer: rs966423 on 2q35 (OR = 1.34; P(combined) = 1.3 × 10(-9)), rs2439302 on 8p12 (OR = 1.36; P(combined) = 2.0 × 10(-9)) and rs116909374 on 14q13.3 (OR = 2.09; P(combined) = 4.6 × 10(-11)), a region previously reported to contain an uncorrelated variant conferring risk of thyroid cancer. A strong association (P = 9.1 × 10(-91)) was observed between rs2439302 on 8p12 and expression of NRG1, which encodes the signaling protein neuregulin 1, in blood.     

Mutations in BRIP1 confer high risk of ovarian cancer

Ovarian cancer causes more deaths than any other gynecologic malignancy in developed countries. Sixteen million sequence variants, identified through whole-genome sequencing of 457 Icelanders, were imputed to 41,675 Icelanders genotyped using SNP chips, as well as to their relatives. Sequence variants were tested for association with ovarian cancer (N of affected individuals = 656). We discovered a rare (0.41% allelic frequency) frameshift mutation, c.2040_2041insTT, in the BRIP1 (FANCJ) gene that confers an increase in ovarian cancer risk (odds ratio (OR) = 8.13, P = 2.8 × 10(-14)). The mutation was also associated with increased risk of cancer in general and reduced lifespan by 3.6 years. In a Spanish population, another frameshift mutation in BRIP1, c.1702_1703del, was seen in 2 out of 144 subjects with ovarian cancer and 1 out of 1,780 control subjects (P = 0.016). This allele was also associated with breast cancer (seen in 6/927 cases; P = 0.0079). Ovarian tumors from heterozygous carriers of the Icelandic mutation show loss of the wild-type allele, indicating that BRIP1 behaves like a classical tumor suppressor gene in ovarian cancer.     

A genome-wide association study identifies three new risk loci for Kawasaki disease

We performed a genome-wide association study (GWAS) of Kawasaki disease in Japanese subjects using data from 428 individuals with Kawasaki disease (cases) and 3,379 controls genotyped at 473,803 SNPs. We validated the association results in two independent replication panels totaling 754 cases and 947 controls. We observed significant associations in the FAM167A-BLK region at 8p22-23 (rs2254546, P = 8.2 × 10(-21)), in the human leukocyte antigen (HLA) region at 6p21.3 (rs2857151, P = 4.6 × 10(-11)) and in the CD40 region at 20q13 (rs4813003, P = 4.8 × 10(-8)). We also replicated the association of a functional SNP of FCGR2A (rs1801274, P = 1.6 × 10(-6)) identified in a recently reported GWAS of Kawasaki disease. Our findings provide new insights into the pathogenesis and pathophysiology of Kawasaki disease.     

A common variant at the TERT-CLPTM1L locus is associated with estrogen receptor-negative breast cancer

Estrogen receptor (ER)-negative breast cancer shows a higher incidence in women of African ancestry compared to women of European ancestry. In search of common risk alleles for ER-negative breast cancer, we combined genome-wide association study (GWAS) data from women of African ancestry (1,004 ER-negative cases and 2,745 controls) and European ancestry (1,718 ER-negative cases and 3,670 controls), with replication testing conducted in an additional 2,292 ER-negative cases and 16,901 controls of European ancestry. We identified a common risk variant for ER-negative breast cancer at the TERT-CLPTM1L locus on chromosome 5p15 (rs10069690: per-allele odds ratio (OR) = 1.18 per allele, P = 1.0 × 10(-10)). The variant was also significantly associated with triple-negative (ER-negative, progesterone receptor (PR)-negative and human epidermal growth factor-2 (HER2)-negative) breast cancer (OR = 1.25, P = 1.1 × 10(-9)), particularly in younger women (<50 years of age) (OR = 1.48, P = 1.9 × 10(-9)). Our results identify a genetic locus associated with estrogen receptor negative breast cancer subtypes in multiple populations.     

Common variants at 11p13 are associated with susceptibility to tuberculosis

After imputation of data from the 1000 Genomes Project into a genome-wide dataset of Ghanaian individuals with tuberculosis and controls, we identified a resistance locus on chromosome 11p13 downstream of the WT1 gene (encoding Wilms tumor 1). The strongest signal was obtained at the rs2057178 SNP (P = 2.63 × 10(-9)). Replication in Gambian, Indonesian and Russian tuberculosis case-control study cohorts increased the significance level for the association with this SNP to P = 2.57 × 10(-11).     

A genome-wide association study of Hodgkin's lymphoma identifies new susceptibility loci at 2p16.1 (REL), 8q24.21 and 10p14 (GATA3)

To identify susceptibility loci for classical Hodgkin's lymphoma (cHL), we conducted a genome-wide association study of 589 individuals with cHL (cases) and 5,199 controls with validation in four independent samples totaling 2,057 cases and 3,416 controls. We identified three new susceptibility loci at 2p16.1 (rs1432295, REL, odds ratio (OR) = 1.22, combined P = 1.91 × 10(-8)), 8q24.21 (rs2019960, PVT1, OR = 1.33, combined P = 1.26 × 10(-13)) and 10p14 (rs501764, GATA3, OR = 1.25, combined P = 7.05 × 10(-8)). Furthermore, we confirmed the role of the major histocompatibility complex in disease etiology by revealing a strong human leukocyte antigen (HLA) association (rs6903608, OR = 1.70, combined P = 2.84 × 10(-50)). These data provide new insight into the pathogenesis of cHL.     

Common sequence variants on 2p15 and Xp11.22 confer susceptibility to prostate cancer

We conducted a genome-wide SNP association study on prostate cancer on over 23,000 Icelanders, followed by a replication study including over 15,500 individuals from Europe and the United States. Two newly identified variants were shown to be associated with prostate cancer: rs5945572 on Xp11.22 and rs721048 on 2p15 (odds ratios (OR) = 1.23 and 1.15; P = 3.9 x 10(-13) and 7.7 x 10(-9), respectively). The 2p15 variant shows a significantly stronger association with more aggressive, rather than less aggressive, forms of the disease.     

Identification of low-frequency variants associated with gout and serum uric acid levels

We tested 16 million SNPs, identified through whole-genome sequencing of 457 Icelanders, for association with gout and serum uric acid levels. Genotypes were imputed into 41,675 chip-genotyped Icelanders and their relatives, for effective sample sizes of 968 individuals with gout and 15,506 individuals for whom serum uric acid measurements were available. We identified a low-frequency missense variant (c.1580C>G) in ALDH16A1 associated with gout (OR = 3.12, P = 1.5 × 10(-16), at-risk allele frequency = 0.019) and serum uric acid levels (effect = 0.36 s.d., P = 4.5 × 10(-21)). We confirmed the association with gout by performing Sanger sequencing on 6,017 Icelanders. The association with gout was stronger in males relative to females. We also found a second variant on chromosome 1 associated with gout (OR = 1.92, P = 0.046, at-risk allele frequency = 0.986) and serum uric acid levels (effect = 0.48 s.d., P = 4.5 × 10(-16)). This variant is close to a common variant previously associated with serum uric acid levels. This work illustrates how whole-genome sequencing data allow the detection of associations between low-frequency variants and complex traits.     

A genome-wide association study in Chinese men identifies three risk loci for non-obstructive azoospermia

Non-obstructive azoospermia (NOA) is one of the most severe forms of male infertility. Its pathophysiology is largely unknown, and few genetic influences have been defined. To identify common variants contributing to NOA in Han Chinese men, we performed a three-stage genome-wide association study of 2,927 individuals with NOA and 5,734 controls. The combined analyses identified significant (P < 5.0 × 10(-8)) associations between NOA risk and common variants near PRMT6 (rs12097821 at 1p13.3: odds ratio (OR) = 1.25, P = 5.7 × 10(-10)), PEX10 (rs2477686 at 1p36.32: OR = 1.39, P = 5.7 × 10(-12)) and SOX5 (rs10842262 at 12p12.1: OR = 1.23, P = 2.3 × 10(-9)). These findings implicate genetic variants at 1p13.3, 1p36.32 and 12p12.1 in the etiology of NOA in Han Chinese men.     

Genome-wide association identifies three new susceptibility loci for Paget's disease of bone

Paget's disease of bone (PDB) is a common disorder characterized by focal abnormalities of bone remodeling. We previously identified variants at the CSF1, OPTN and TNFRSF11A loci as risk factors for PDB by genome-wide association study. Here we extended this study, identified three new loci and confirmed their association with PDB in 2,215 affected individuals (cases) and 4,370 controls from seven independent populations. The new associations were with rs5742915 within PML on 15q24 (odds ratio (OR) = 1.34, P = 1.6 × 10(-14)), rs10498635 within RIN3 on 14q32 (OR = 1.44, P = 2.55 × 10(-11)) and rs4294134 within NUP205 on 7q33 (OR = 1.45, P = 8.45 × 10(-10)). Our data also confirmed the association of TM7SF4 (rs2458413, OR = 1.40, P = 7.38 × 10(-17)) with PDB. These seven loci explained ～13% of the familial risk of PDB. These studies provide new insights into the genetic architecture and pathophysiology of PDB.     

Sequence variant on 8q24 confers susceptibility to urinary bladder cancer

We conducted a genome-wide SNP association study on 1,803 urinary bladder cancer (UBC) cases and 34,336 controls from Iceland and The Netherlands and follow up studies in seven additional case-control groups (2,165 cases and 3,800 controls). The strongest association was observed with allele T of rs9642880 on chromosome 8q24, 30 kb upstream of MYC (allele-specific odds ratio (OR) = 1.22; P = 9.34 x 10(-12)). Approximately 20% of individuals of European ancestry are homozygous for rs9642880[T], and their estimated risk of developing UBC is 1.49 times that of noncarriers. No association was observed between UBC and the four 8q24 variants previously associated with prostate, colorectal and breast cancers, nor did rs9642880 associate with any of these three cancers. A weaker signal, but nonetheless of genome-wide significance, was captured by rs710521[A] located near TP63 on chromosome 3q28 (allele-specific OR = 1.19; P = 1. 15 x 10(-7)).     

A genome-wide association study identifies two new lung cancer susceptibility loci at 13q12.12 and 22q12.2 in Han Chinese

Lung cancer is the leading cause of cancer-related deaths worldwide. To identify genetic factors that modify the risk of lung cancer in individuals of Chinese ancestry, we performed a genome-wide association scan in 5,408 subjects (2,331 individuals with lung cancer (cases) and 3,077 controls) followed by a two-stage validation among 12,722 subjects (6,313 cases and 6,409 controls). The combined analyses identified six well-replicated SNPs with independent effects and significant lung cancer associations (P < 5.0 × 10(-8)) located in TP63 (rs4488809 at 3q28, P = 7.2 × 10(-26)), TERT-CLPTM1L (rs465498 and rs2736100 at 5p15.33, P = 1.2 × 10(-20) and P = 1.0 × 10(-27), respectively), MIPEP-TNFRSF19 (rs753955 at 13q12.12, P = 1.5 × 10(-12)) and MTMR3-HORMAD2-LIF (rs17728461 and rs36600 at 22q12.2, P = 1.1 × 10(-11) and P = 6.2 × 10(-13), respectively). Two of these loci (13q12.12 and 22q12.2) were newly identified in the Chinese population. These results suggest that genetic variants in 3q28, 5p15.33, 13q12.12 and 22q12.2 may contribute to the susceptibility of lung cancer in Han Chinese.     

A genome-wide association study identifies two susceptibility loci for duodenal ulcer in the Japanese population

Through a genome-wide association analysis with a total of 7,035 individuals with duodenal ulcer and 25,323 controls from Japan, we identified two susceptibility loci at the PSCA gene (encoding prostate stem cell antigen) at 8q24 and at the ABO blood group locus at 9q34. The C allele of rs2294008 at PSCA was associated with increased risk of duodenal ulcer (odds ratio (OR) = 1.84; P = 3.92 × 10(-33)) in a recessive model but was associated with decreased risk of gastric cancer (OR = 0.79; P = 6.79 × 10(-12)), as reported previously. The T allele of rs2294008 encodes a translation initiation codon upstream of the reported site and changes protein localization from the cytoplasm to the cell surface. rs505922 at ABO was also associated with duodenal ulcer in a recessive model (OR = 1.32; P = 1.15 × 10(-10)). Our findings demonstrate a role for genetic variants in the pathogenesis of duodenal ulcer.     

Genome-wide association study identifies three new susceptibility loci for adult asthma in the Japanese population

Bronchial asthma is a common inflammatory disease caused by the interaction of genetic and environmental factors. Through a genome-wide association study and a replication study consisting of a total of 7,171 individuals with adult asthma (cases) and 27,912 controls in the Japanese population, we identified five loci associated with susceptibility to adult asthma. In addition to the major histocompatibility complex and TSLP-WDR36 loci previously reported, we identified three additional loci: a USP38-GAB1 locus on chromosome 4q31 (combined P = 1.87 × 10(-12)), a locus on chromosome 10p14 (P = 1.79 × 10(-15)) and a gene-rich region on chromosome 12q13 (P = 2.33 × 10(-13)). We observed the most significant association with adult asthma at rs404860 in the major histocompatiblity complex region (P = 4.07 × 10(-23)), which is close to rs2070600, a SNP previously reported for association with FEV(1)/FVC in genome-wide association studies for lung function. Our findings offer a better understanding of the genetic contribution to asthma susceptibility.     

Common 5p15.33 and 6p21.33 variants influence lung cancer risk

We conducted a genome-wide association (GWA) study of lung cancer comparing 511,919 SNP genotypes in 1,952 cases and 1,438 controls. The most significant association was attained at 15q25.1 (rs8042374; P = 7.75 x 10(-12)), confirming recent observations. Pooling data with two other GWA studies (5,095 cases, 5,200 controls) and with replication in an additional 2,484 cases and 3,036 controls, we identified two newly associated risk loci mapping to 6p21.33 (rs3117582, BAT3-MSH5; P(combined) = 4.97 x 10(-10)) and 5p15.33 (rs401681, CLPTM1L; P(combined) = 7.90 x 10(-9)).     

Association analyses identify multiple new lung cancer susceptibility loci and their interactions with smoking in the Chinese population

To find additional susceptibility loci for lung cancer, we tested promising associations from our previous genome-wide association study (GWAS) of lung cancer in the Chinese population in an extended validation sample size of 7,436 individuals with lung cancer (cases) and 7,483 controls. We found genome-wide significant (P < 5.0 × 10(-8)) evidence for three additional lung cancer susceptibility loci at 10p14 (rs1663689, close to GATA3, P = 2.84 × 10(-10)), 5q32 (rs2895680 in PPP2R2B-STK32A-DPYSL3, P = 6.60 × 10(-9)) and 20q13.2 (rs4809957 in CYP24A1, P = 1.20 × 10(-8)). We also found consistent associations for rs247008 at 5q31.1 (IL3-CSF2-P4HA2, P = 7.68 × 10(-8)) and rs9439519 at 1p36.32 (AJAP1-NPHP4, P = 3.65 × 10(-6)). Four of these loci showed evidence for interactions with smoking dose (P = 1.72 × 10(-10), P = 5.07 × 10(-3), P = 6.77 × 10(-3) and P = 4.49 × 10(-2) for rs2895680, rs4809957, rs247008 and rs9439519, respectively). These results advance our understanding of lung cancer susceptibility and highlight potential pathways that integrate genetic variants and smoking in the development of lung cancer.     

Common variation at 3p22.1 and 7p15.3 influences multiple myeloma risk

To identify risk variants for multiple myeloma, we conducted a genome-wide association study of 1,675 individuals with multiple myeloma and 5,903 control subjects. We identified risk loci for multiple myeloma at 3p22.1 (rs1052501 in ULK4; odds ratio (OR) = 1.32; P = 7.47 × 10(-9)) and 7p15.3 (rs4487645, OR = 1.38; P = 3.33 × 10(-15)). In addition, we observed a promising association at 2p23.3 (rs6746082, OR = 1.29; P = 1.22 × 10(-7)). Our study identifies new genomic regions associated with multiple myeloma risk that may lead to new etiological insights.     

Common variants at CDKAL1 and KLF9 are associated with body mass index in east Asian populations

Obesity is a disorder with a complex genetic etiology, and its epidemic is a worldwide problem. Although multiple genetic loci associated with body mass index, the most common measure of obesity, have been identified in European populations, few studies have focused on Asian populations. Here we report a genome-wide association study and replication studies with 62,245 east Asian subjects, which identified two new body mass index-associated loci in the CDKAL1 locus at 6p22 (rs2206734, P = 1.4 × 10(-11)) and the KLF9 locus at 9q21 (rs11142387, P = 1.3 × 10(-9)), as well as several previously reported loci (the SEC16B, BDNF, FTO, MC4R and GIPR loci, P < 5.0 × 10(-8)). We subsequently performed gene-gene interaction analyses and identified an interaction (P = 2.0 × 10(-8)) between a SNP in the KLF9 locus (rs11142387) and one in the MSTN (also known as GDF8) locus at 2q32 (rs13034723). These findings should provide useful insights into the etiology of obesity.     

Meta-analysis of genome-wide association studies identifies three new risk loci for atopic dermatitis

Atopic dermatitis (AD) is a commonly occurring chronic skin disease with high heritability. Apart from filaggrin (FLG), the genes influencing atopic dermatitis are largely unknown. We conducted a genome-wide association meta-analysis of 5,606 affected individuals and 20,565 controls from 16 population-based cohorts and then examined the ten most strongly associated new susceptibility loci in an additional 5,419 affected individuals and 19,833 controls from 14 studies. Three SNPs reached genome-wide significance in the discovery and replication cohorts combined, including rs479844 upstream of OVOL1 (odds ratio (OR) = 0.88, P = 1.1 × 10(-13)) and rs2164983 near ACTL9 (OR = 1.16, P = 7.1 × 10(-9)), both of which are near genes that have been implicated in epidermal proliferation and differentiation, as well as rs2897442 in KIF3A within the cytokine cluster at 5q31.1 (OR = 1.11, P = 3.8 × 10(-8)). We also replicated association with the FLG locus and with two recently identified association signals at 11q13.5 (rs7927894; P = 0.008) and 20q13.33 (rs6010620; P = 0.002). Our results underline the importance of both epidermal barrier function and immune dysregulation in atopic dermatitis pathogenesis.     

Common variants at MS4A4/MS4A6E, CD2AP, CD33 and EPHA1 are associated with late-onset Alzheimer's disease

The Alzheimer Disease Genetics Consortium (ADGC) performed a genome-wide association study of late-onset Alzheimer disease using a three-stage design consisting of a discovery stage (stage 1) and two replication stages (stages 2 and 3). Both joint analysis and meta-analysis approaches were used. We obtained genome-wide significant results at MS4A4A (rs4938933; stages 1 and 2, meta-analysis P (P(M)) = 1.7 × 10(-9), joint analysis P (P(J)) = 1.7 × 10(-9); stages 1, 2 and 3, P(M) = 8.2 × 10(-12)), CD2AP (rs9349407; stages 1, 2 and 3, P(M) = 8.6 × 10(-9)), EPHA1 (rs11767557; stages 1, 2 and 3, P(M) = 6.0 × 10(-10)) and CD33 (rs3865444; stages 1, 2 and 3, P(M) = 1.6 × 10(-9)). We also replicated previous associations at CR1 (rs6701713; P(M) = 4.6 × 10(-10), P(J) = 5.2 × 10(-11)), CLU (rs1532278; P(M) = 8.3 × 10(-8), P(J) = 1.9 × 10(-8)), BIN1 (rs7561528; P(M) = 4.0 × 10(-14), P(J) = 5.2 × 10(-14)) and PICALM (rs561655; P(M) = 7.0 × 10(-11), P(J) = 1.0 × 10(-10)), but not at EXOC3L2, to late-onset Alzheimer's disease susceptibility.