Common variants in WFS1 confer risk of type 2 diabetes

We studied genes involved in pancreatic beta cell function and survival, identifying associations between SNPs in WFS1 and diabetes risk in UK populations that we replicated in an Ashkenazi population and in additional UK studies. In a pooled analysis comprising 9,533 cases and 11,389 controls, SNPs in WFS1 were strongly associated with diabetes risk. Rare mutations in WFS1 cause Wolfram syndrome; using a gene-centric approach, we show that variation in WFS1 also predisposes to common type 2 diabetes.     

Common variants on 8p12 and 1q24.2 confer risk of schizophrenia

Schizophrenia is a severe mental disorder affecting ～1% of the world population, with heritability of up to 80%. To identify new common genetic risk factors, we performed a genome-wide association study (GWAS) in the Han Chinese population. The discovery sample set consisted of 3,750 individuals with schizophrenia and 6,468 healthy controls (1,578 cases and 1,592 controls from northern Han Chinese, 1,238 cases and 2,856 controls from central Han Chinese, and 934 cases and 2,020 controls from the southern Han Chinese). We further analyzed the strongest association signals in an additional independent cohort of 4,383 cases and 4,539 controls from the Han Chinese population. Meta-analysis identified common SNPs that associated with schizophrenia with genome-wide significance on 8p12 (rs16887244, P = 1.27 × 10(-10)) and 1q24.2 (rs10489202, P = 9.50 × 10(-9)). Our findings provide new insights into the pathogenesis of schizophrenia.     

Common variants at CD40 and other loci confer risk of rheumatoid arthritis

To identify rheumatoid arthritis risk loci in European populations, we conducted a meta-analysis of two published genome-wide association (GWA) studies totaling 3,393 cases and 12,462 controls. We genotyped 31 top-ranked SNPs not previously associated with rheumatoid arthritis in an independent replication of 3,929 autoantibody-positive rheumatoid arthritis cases and 5,807 matched controls from eight separate collections. We identified a common variant at the CD40 gene locus (rs4810485, P = 0.0032 replication, P = 8.2 x 10(-9) overall, OR = 0.87). Along with other associations near TRAF1 (refs. 2,3) and TNFAIP3 (refs. 4,5), this implies a central role for the CD40 signaling pathway in rheumatoid arthritis pathogenesis. We also identified association at the CCL21 gene locus (rs2812378, P = 0.00097 replication, P = 2.8 x 10(-7) overall), a gene involved in lymphocyte trafficking. Finally, we identified evidence of association at four additional gene loci: MMEL1-TNFRSF14 (rs3890745, P = 0.0035 replication, P = 1.1 x 10(-7) overall), CDK6 (rs42041, P = 0.010 replication, P = 4.0 x 10(-6) overall), PRKCQ (rs4750316, P = 0.0078 replication, P = 4.4 x 10(-6) overall), and KIF5A-PIP4K2C (rs1678542, P = 0.0026 replication, P = 8.8 x 10(-8) overall).     

Common sequence variants on 20q11.22 confer melanoma susceptibility

We conducted a genome-wide association pooling study for cutaneous melanoma and performed validation in samples totaling 2,019 cases and 2,105 controls. Using pooling, we identified a new melanoma risk locus on chromosome 20 (rs910873 and rs1885120), with replication in two further samples (combined P < 1 x 10(-15)). The per allele odds ratio was 1.75 (1.53, 2.01), with evidence for stronger association in early-onset cases.     

Common variants near MC4R are associated with fat mass, weight and risk of obesity

To identify common variants influencing body mass index (BMI), we analyzed genome-wide association data from 16,876 individuals of European descent. After previously reported variants in FTO, the strongest association signal (rs17782313, P = 2.9 x 10(-6)) mapped 188 kb downstream of MC4R (melanocortin-4 receptor), mutations of which are the leading cause of monogenic severe childhood-onset obesity. We confirmed the BMI association in 60,352 adults (per-allele effect = 0.05 Z-score units; P = 2.8 x 10(-15)) and 5,988 children aged 7-11 (0.13 Z-score units; P = 1.5 x 10(-8)). In case-control analyses (n = 10,583), the odds for severe childhood obesity reached 1.30 (P = 8.0 x 10(-11)). Furthermore, we observed overtransmission of the risk allele to obese offspring in 660 families (P (pedigree disequilibrium test average; PDT-avg) = 2.4 x 10(-4)). The SNP location and patterns of phenotypic associations are consistent with effects mediated through altered MC4R function. Our findings establish that common variants near MC4R influence fat mass, weight and obesity risk at the population level and reinforce the need for large-scale data integration to identify variants influencing continuous biomedical traits.     

Common variants on chromosome 5p12 confer susceptibility to estrogen receptor-positive breast cancer

We carried out a genome-wide association study of breast cancer predisposition with replication and refinement studies involving 6,145 cases and 33,016 controls and identified two SNPs (rs4415084 and rs10941679) on 5p12 that confer risk, preferentially for estrogen receptor (ER)-positive tumors (OR = 1.27, P = 2.5 x 10(-12) for rs10941679). The nearest gene, MRPS30, was previously implicated in apoptosis, ER-positive tumors and favorable prognosis. A recently reported signal in FGFR2 was also found to associate specifically with ER-positive breast cancer.     

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.     

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.     

Common variants in DGKK are strongly associated with risk of hypospadias

Hypospadias is a common congenital malformation of the male external genitalia. We performed a genome-wide association study using pooled DNA from 436 individuals with hypospadias (cases) and 494 controls of European descent and selected the highest ranked SNPs for individual genotyping in the discovery sample, an additional Dutch sample of 133 cases and their parents, and a Swedish series of 266 cases and 402 controls. Individual genotyping of two SNPs (rs1934179 and rs7063116) in DGKK, encoding diacylglycerol kinase κ, produced compelling evidence for association with hypospadias in the discovery sample (allele-specific odds ratio (OR) = 2.5, P = 2.5 × 10?11 and OR = 2.3, P = 2.9 × 10??, respectively) and in the Dutch (OR = 3.9, P = 2.4 × 10?? and OR = 3.8, P = 3.4 × 10??) and Swedish (OR = 2.5, P = 2.6 × 10?? and OR = 2.2, P = 2.7 × 10??) replication samples. Expression studies showed expression of DGKK in preputial tissue of cases and controls, which was lower in carriers of the risk allele of rs1934179 (P = 0.047). We propose DGKK as a major risk gene for hypospadias.     

Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants

We have genotyped 14,436 nonsynonymous SNPs (nsSNPs) and 897 major histocompatibility complex (MHC) tag SNPs from 1,000 independent cases of ankylosing spondylitis (AS), autoimmune thyroid disease (AITD), multiple sclerosis (MS) and breast cancer (BC). Comparing these data against a common control dataset derived from 1,500 randomly selected healthy British individuals, we report initial association and independent replication in a North American sample of two new loci related to ankylosing spondylitis, ARTS1 and IL23R, and confirmation of the previously reported association of AITD with TSHR and FCRL3. These findings, enabled in part by increased statistical power resulting from the expansion of the control reference group to include individuals from the other disease groups, highlight notable new possibilities for autoimmune regulation and suggest that IL23R may be a common susceptibility factor for the major 'seronegative' diseases.     

Common variants on chromosomes 2q35 and 16q12 confer susceptibility to estrogen receptor-positive breast cancer

Familial clustering studies indicate that breast cancer risk has a substantial genetic component. To identify new breast cancer risk variants, we genotyped approximately 300,000 SNPs in 1,600 Icelandic individuals with breast cancer and 11,563 controls using the Illumina Hap300 platform. We then tested selected SNPs in five replication sample sets. Overall, we studied 4,554 affected individuals and 17,577 controls. Two SNPs consistently associated with breast cancer: approximately 25% of individuals of European descent are homozygous for allele A of rs13387042 on chromosome 2q35 and have an estimated 1.44-fold greater risk than noncarriers, and for allele T of rs3803662 on 16q12, about 7% are homozygous and have a 1.64-fold greater risk. Risk from both alleles was confined to estrogen receptor-positive tumors. At present, no genes have been identified in the linkage disequilibrium block containing rs13387042. rs3803662 is near the 5' end of TNRC9 , a high mobility group chromatin-associated protein whose expression is implicated in breast cancer metastasis to bone.     

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)).     

Common variants in ZNF365 are associated with both mammographic density and breast cancer risk

High-percent mammographic density adjusted for age and body mass index is one of the strongest risk factors for breast cancer. We conducted a meta analysis of five genome-wide association studies of percent mammographic density and report an association with rs10995190 in ZNF365 (combined P = 9.6 × 10(-10)). Common variants in ZNF365 have also recently been associated with susceptibility to breast cancer.     

Common genetic variants at the CRAC1 (HMPS) locus on chromosome 15q13.3 influence colorectal cancer risk

We mapped a high-penetrance gene (CRAC1; also known as HMPS) associated with colorectal cancer (CRC) in the Ashkenazi population to a 0.6-Mb region on chromosome 15 containing SCG5 (also known as SGNE1), GREM1 and FMN1. We hypothesized that the CRAC1 locus harbored low-penetrance variants that increased CRC risk in the general population. In a large series of colorectal cancer cases and controls, SNPs near GREM1 and SCG5 were strongly associated with increased CRC risk (for rs4779584, P = 4.44 x 10(-14)).     

Two variants on chromosome 17 confer prostate cancer risk, and the one in TCF2 protects against type 2 diabetes

We performed a genome-wide association scan to search for sequence variants conferring risk of prostate cancer using 1,501 Icelandic men with prostate cancer and 11,290 controls. Follow-up studies involving three additional case-control groups replicated an association of two variants on chromosome 17 with the disease. These two variants, 33 Mb apart, fall within a region previously implicated by family-based linkage studies on prostate cancer. The risks conferred by these variants are moderate individually (allele odds ratio of about 1.20), but because they are common, their joint population attributable risk is substantial. One of the variants is in TCF2 (HNF1beta), a gene known to be mutated in individuals with maturity-onset diabetes of the young type 5. Results from eight case-control groups, including one West African and one Chinese, demonstrate that this variant confers protection against type 2 diabetes.     

Common and rare variants in multifactorial susceptibility to common diseases

Here, we give a historical overview of the search for genetic variants that influence the susceptibility of an individual to a chronic disease, from RA Fisher's seminal work to the current excitement of whole-genome association studies (WGAS). We then discuss the concepts behind the identification of common variants as disease causal factors and contrast them to the basic ideas that underlie the rare variant hypothesis. The identification of rare variants involves the careful selection of candidate genes to examine, the availability of highly efficient resequencing techniques and the appropriate assessment of the functional consequences of the implicated variant. We believe that this strategy can be successfully applied at present in order to unravel the contribution of rare variants to the multifactorial inheritance of common diseases, which could lead to the implementation of much needed preventative screening schemes.     

A nonsynonymous SNP in PRKCH (protein kinase C eta) increases the risk of cerebral infarction

Cerebral infarction is the most common type of stroke and often causes long-term disability. To investigate the genetic contribution to cerebral infarction, we conducted a case-control study using 52,608 gene-based tag SNPs selected from the JSNP database. Here we report that a nonsynonymous SNP in a member of protein kinase C (PKC) family, PRKCH, was significantly associated with lacunar infarction in two independent Japanese samples (P = 5.1 x 10(-7), crude odds ratio of 1.40). This SNP is likely to affect PKC activity. Furthermore, a 14-year follow-up cohort study in Hisayama (Fukuoka, Japan) supported involvement of this SNP in the development of cerebral infarction (P = 0.03, age- and sex-adjusted hazard ratio of 2.83). We also found that PKCeta was expressed mainly in vascular endothelial cells and foamy macrophages in human atherosclerotic lesions, and its expression increased as the lesion type progressed. Our results support a role for PRKCH in the pathogenesis of cerebral infarction.     

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.     

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.     

Common variants near MBNL1 and NKX2-5 are associated with infantile hypertrophic pyloric stenosis

Infantile hypertrophic pyloric stenosis (IHPS) is a severe condition characterized by hypertrophy of the pyloric sphincter muscle. We conducted a genome-wide association study (GWAS) on 1,001 surgery-confirmed cases and 2,401 controls from Denmark. The six most strongly associated loci were tested in a replication set of 796 cases and 876 controls. Three SNPs reached genome-wide significance. One of these SNPs, rs11712066 (odds ratio (OR) = 1.61; P = 1.5 × 10(-17)) at 3p25.1, is located 150 kb upstream of MBNL1, which encodes a factor that regulates splicing transitions occurring shortly after birth. The second SNP, rs573872 (OR = 1.41; P = 4.3 × 10(-12)), maps to an intergenic region at 3p25.2 approximately 1.3 Mb downstream of MBNL1. The third SNP, rs29784 (OR = 1.42; P = 1.5 × 10(-15)) at 5q35.2, is 64 kb downstream of NKX2-5, which is involved in development of cardiac muscle tissue and embryonic gut development.