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.     

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.     

NAD(P)H:quinone oxidoreductase 1 NQO1*2 genotype (P187S) is a strong prognostic and predictive factor in breast cancer

NQO1 guards against oxidative stress and carcinogenesis and stabilizes p53. We find that a homozygous common missense variant (NQO1(*)2, rs1800566(T), NM_000903.2:c.558C>T) that disables NQO1 strongly predicts poor survival among two independent series of women with breast cancer (P = 0.002, N = 1,005; P = 0.005, N = 1,162), an effect particularly evident after anthracycline-based adjuvant chemotherapy with epirubicin (P = 7.52 x 10(-6)) and in p53-aberrant tumors (P = 6.15 x 10(-5)). Survival after metastasis was reduced among NQO1(*)2 homozygotes, further implicating NQO1 deficiency in cancer progression and treatment resistance. Consistently, response to epirubicin was impaired in NQO1(*)2-homozygous breast carcinoma cells in vitro, reflecting both p53-linked and p53-independent roles of NQO1. We propose a model of defective anthracycline response in NQO1-deficient breast tumors, along with increased genomic instability promoted by elevated reactive oxygen species (ROS), and suggest that the NQO1 genotype is a prognostic and predictive marker for breast cancer.     

A germline variant in the TP53 polyadenylation signal confers cancer susceptibility

To identify new risk variants for cutaneous basal cell carcinoma, we performed a genome-wide association study of 16 million SNPs identified through whole-genome sequencing of 457 Icelanders. We imputed genotypes for 41,675 Illumina SNP chip-typed Icelanders and their relatives. In the discovery phase, the strongest signal came from rs78378222[C] (odds ratio (OR) = 2.36, P = 5.2 × 10(-17)), which has a frequency of 0.0192 in the Icelandic population. We then confirmed this association in non-Icelandic samples (OR = 1.75, P = 0.0060; overall OR = 2.16, P = 2.2 × 10(-20)). rs78378222 is in the 3' untranslated region of TP53 and changes the AATAAA polyadenylation signal to AATACA, resulting in impaired 3'-end processing of TP53 mRNA. Investigation of other tumor types identified associations of this SNP with prostate cancer (OR = 1.44, P = 2.4 × 10(-6)), glioma (OR = 2.35, P = 1.0 × 10(-5)) and colorectal adenoma (OR = 1.39, P = 1.6 × 10(-4)). However, we observed no effect for breast cancer, a common Li-Fraumeni syndrome tumor (OR = 1.06, P = 0.57, 95% confidence interval 0.88-1.27).     

A synonymous variant in IRGM alters a binding site for miR-196 and causes deregulation of IRGM-dependent xenophagy in Crohn's disease

Susceptibility to Crohn's disease, a complex inflammatory disease, is influenced by common variants at many loci. The common exonic synonymous SNP (c.313C>T) in IRGM, found in strong linkage disequilibrium with a deletion polymorphism, has been classified as non-causative because of the absence of an alteration in the IRGM protein sequence or splice sites. Here we show that a family of microRNAs (miRNAs), miR-196, is overexpressed in the inflammatory intestinal epithelia of individuals with Crohn's disease and downregulates the IRGM protective variant (c.313C) but not the risk-associated allele (c.313T). Subsequent loss of tight regulation of IRGM expression compromises control of intracellular replication of Crohn's disease-associated adherent invasive Escherichia coli by autophagy. These results suggest that the association of IRGM with Crohn's disease arises from a miRNA-based alteration in IRGM regulation that affects the efficacy of autophagy, thereby implicating a synonymous polymorphism as a likely causal variant.     

Genetic determinants of hair, eye and skin pigmentation in Europeans

Hair, skin and eye colors are highly heritable and visible traits in humans. We carried out a genome-wide association scan for variants associated with hair and eye pigmentation, skin sensitivity to sun and freckling among 2,986 Icelanders. We then tested the most closely associated SNPs from six regions--four not previously implicated in the normal variation of human pigmentation--and replicated their association in a second sample of 2,718 Icelanders and a sample of 1,214 Dutch. The SNPs from all six regions met the criteria for genome-wide significance. A variant in SLC24A4 is associated with eye and hair color, a variant near KITLG is associated with hair color, two coding variants in TYR are associated with eye color and freckles, and a variant on 6p25.3 is associated with freckles. The fifth region provided refinements to a previously reported association in OCA2, and the sixth encompasses previously described variants in MC1R.     

A recurrent mutation in MED12 leading to R961W causes Opitz-Kaveggia syndrome

Opitz-Kaveggia syndrome (also known as FG syndrome) is an X-linked disorder characterized by mental retardation, relative macrocephaly, hypotonia and constipation. We report here that the original family for whom the condition is named and five other families have a recurrent mutation (2881C>T, leading to R961W) in MED12 (also called TRAP230 or HOPA), a gene located at Xq13 that functions as a thyroid receptor-associated protein in the Mediator complex.     

Somatic mosaic IDH1 and IDH2 mutations are associated with enchondroma and spindle cell hemangioma in Ollier disease and Maffucci syndrome

Ollier disease and Maffucci syndrome are non-hereditary skeletal disorders characterized by multiple enchondromas (Ollier disease) combined with spindle cell hemangiomas (Maffucci syndrome). We report somatic heterozygous mutations in IDH1 (c.394C>T encoding an R132C substitution and c.395G>A encoding an R132H substitution) or IDH2 (c.516G>C encoding R172S) in 87% of enchondromas (benign cartilage tumors) and in 70% of spindle cell hemangiomas (benign vascular lesions). In total, 35 of 43 (81%) subjects with Ollier disease and 10 of 13 (77%) with Maffucci syndrome carried IDH1 (98%) or IDH2 (2%) mutations in their tumors. Fourteen of 16 subjects had identical mutations in separate lesions. Immunohistochemistry to detect mutant IDH1 R132H protein suggested intraneoplastic and somatic mosaicism. IDH1 mutations in cartilage tumors were associated with hypermethylation and downregulated expression of several genes. Mutations were also found in 40% of solitary central cartilaginous tumors and in four chondrosarcoma cell lines, which will enable functional studies to assess the role of IDH1 and IDH2 mutations in tumor formation.     

Two newly identified genetic determinants of pigmentation in Europeans

We present results from a genome-wide association study for variants associated with human pigmentation characteristics among 5,130 Icelanders, with follow-up analyses in 2,116 Icelanders and 1,214 Dutch individuals. Two coding variants in TPCN2 are associated with hair color, and a variant at the ASIP locus shows strong association with skin sensitivity to sun, freckling and red hair, phenotypic characteristics similar to those affected by well-known mutations in MC1R.     

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.     

Common variants on 1p36 and 1q42 are associated with cutaneous basal cell carcinoma but not with melanoma or pigmentation traits

To search for new sequence variants that confer risk of cutaneous basal cell carcinoma (BCC), we conducted a genome-wide SNP association study of 930 Icelanders with BCC and 33,117 controls. After analyzing 304,083 SNPs, we observed signals from loci at 1p36 and 1q42, and replicated these associations in additional sample sets from Iceland and Eastern Europe. Overall, the most significant signals were from rs7538876 on 1p36 (OR = 1.28, P = 4.4 x 10(-12)) and rs801114 on 1q42 (OR = 1.28, P = 5.9 x 10(-12)). The 1p36 locus contains the candidate genes PADI4, PADI6, RCC2 and ARHGEF10L, and the gene nearest to the 1q42 locus is the ras-homolog RHOU. Neither locus was associated with fair pigmentation traits that are known risk factors for BCC, and no risk was observed for melanoma. Approximately 1.6% of individuals of European ancestry are homozygous for both variants, and their estimated risk of BCC is 2.68 times that of noncarriers.     

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.     

Common variants in P2RY11 are associated with narcolepsy

Growing evidence supports the hypothesis that narcolepsy with cataplexy is an autoimmune disease. We here report genome-wide association analyses for narcolepsy with replication and fine mapping across three ethnic groups (3,406 individuals of European ancestry, 2,414 Asians and 302 African Americans). We identify a SNP in the 3' untranslated region of P2RY11, the purinergic receptor subtype P2Y?? gene, which is associated with narcolepsy (rs2305795, combined P = 6.1 × 10?1?, odds ratio = 1.28, 95% CI 1.19-1.39, n = 5689). The disease-associated allele is correlated with reduced expression of P2RY11 in CD8(+) T lymphocytes (339% reduced, P = 0.003) and natural killer (NK) cells (P = 0.031), but not in other peripheral blood mononuclear cell types. The low expression variant is also associated with reduced P2RY11-mediated resistance to ATP-induced cell death in T lymphocytes (P = 0.0007) and natural killer cells (P = 0.001). These results identify P2RY11 as an important regulator of immune-cell survival, with possible implications in narcolepsy and other autoimmune diseases.     

The same sequence variant on 9p21 associates with myocardial infarction, abdominal aortic aneurysm and intracranial aneurysm

Recently, two common sequence variants on 9p21, tagged by rs10757278-G and rs10811661-T, were reported to be associated with coronary artery disease (CAD) and type 2 diabetes (T2D), respectively. We proceeded to further investigate the contributions of these variants to arterial diseases and T2D. Here we report that rs10757278-G is associated with, in addition to CAD, abdominal aortic aneurysm (AAA; odds ratio (OR) = 1.31, P = 1.2 x 10(-12)) and intracranial aneurysm (OR = 1.29, P = 2.5 x 10(-6)), but not with T2D. This variant is the first to be described that affects the risk of AAA and intracranial aneurysm in many populations. The association of rs10811661-T to T2D replicates in our samples, but the variant does not associate with any of the five arterial diseases examined. These findings extend our insight into the role of the sequence variant tagged by rs10757278-G and show that it is not confined to atherosclerotic diseases.     

Large-scale genetic fine mapping and genotype-phenotype associations implicate polymorphism in the IL2RA region in type 1 diabetes

Genome-wide association studies are now identifying disease-associated chromosome regions. However, even after convincing replication, the localization of the causal variant(s) requires comprehensive resequencing, extensive genotyping and statistical analyses in large sample sets leading to targeted functional studies. Here, we have localized the type 1 diabetes (T1D) association in the interleukin 2 receptor alpha (IL2RA) gene region to two independent groups of SNPs, spanning overlapping regions of 14 and 40 kb, encompassing IL2RA intron 1 and the 5' regions of IL2RA and RBM17 (odds ratio = 2.04, 95% confidence interval = 1.70-2.45; P = 1.92 x 10(-28); control frequency = 0.635). Furthermore, we have associated IL2RA T1D susceptibility genotypes with lower circulating levels of the biomarker, soluble IL-2RA (P = 6.28 x 10(-28)), suggesting that an inherited lower immune responsiveness predisposes to T1D.     

Postzygotic HRAS and KRAS mutations cause nevus sebaceous and Schimmelpenning syndrome

Nevus sebaceous is a common congenital cutaneous malformation. Affected individuals may develop benign and malignant secondary tumors in the nevi during life. Schimmelpenning syndrome is characterized by the association of nevus sebaceous with extracutaneous abnormalities. We report that of 65 sebaceous nevi studied, 62 (95%) had mutations in the HRAS gene and 3 (5%) had mutations in the KRAS gene. The HRAS c.37G>C mutation, which results in a p.Gly13Arg substitution, was present in 91% of lesions. Nonlesional tissues from 18 individuals had a wild-type sequence, confirming genetic mosaicism. The HRAS c.37G>C mutation was also found in 8 of 8 associated secondary tumors. Mosaicism for HRAS c.37G>C and KRAS c.35G>A mutations was found in two individuals with Schimmelpenning syndrome. Functional analysis of HRAS c.37G>C mutant cells showed constitutive activation of the MAPK and PI3K-Akt signaling pathways. Our results indicate that nevus sebaceous and Schimmelpenning syndrome are caused by postzygotic HRAS and KRAS mutations. These mutations may predispose individuals to the development of secondary tumors in nevus sebaceous.