A genome-wide association study identifies colorectal cancer susceptibility loci on chromosomes 10p14 and 8q23.3

To identify colorectal cancer (CRC) susceptibility alleles, we conducted a genome-wide association study. In phase 1, we genotyped 550,163 tagSNPs in 940 familial colorectal tumor cases (627 CRC, 313 high-risk adenoma) and 965 controls. In phase 2, we genotyped 42,708 selected SNPs in 2,873 CRC cases and 2,871 controls. In phase 3, we evaluated 11 SNPs showing association at P < 10(-4) in a joint analysis of phases 1 and 2 in 4,287 CRC cases and 3,743 controls. Two SNPs were taken forward to phase 4 genotyping (10,731 CRC cases and 10,961 controls from eight centers). In addition to the previously reported 8q24, 15q13 and 18q21 CRC risk loci, we identified two previously unreported associations: rs10795668, located at 10p14 (P = 2.5 x 10(-13) overall; P = 6.9 x 10(-12) replication), and rs16892766, at 8q23.3 (P = 3.3 x 10(-18) overall; P = 9.6 x 10(-17) replication), which tags a plausible causative gene, EIF3H. These data provide further evidence for the 'common-disease common-variant' model of CRC predisposition.     

A genome-wide association scan of tag SNPs identifies a susceptibility variant for colorectal cancer at 8q24.21

Much of the variation in inherited risk of colorectal cancer (CRC) is probably due to combinations of common low risk variants. We conducted a genome-wide association study of 550,000 tag SNPs in 930 familial colorectal tumor cases and 960 controls. The most strongly associated SNP (P = 1.72 x 10(-7), allelic test) was rs6983267 at 8q24.21. To validate this finding, we genotyped rs6983267 in three additional CRC case-control series (4,361 affected individuals and 3,752 controls; 1,901 affected individuals and 1,079 controls; 1,072 affected individuals and 415 controls) and replicated the association, providing P = 1.27 x 10(-14) (allelic test) overall, with odds ratios (ORs) of 1.27 (95% confidence interval (c.i.): 1.16-1.39) and 1.47 (95% c.i.: 1.34-1.62) for heterozygotes and rare homozygotes, respectively. Analyses based on 1,477 individuals with colorectal adenoma and 2,136 controls suggest that susceptibility to CRC is mediated through development of adenomas (OR = 1.21, 95% c.i.: 1.10-1.34; P = 6.89 x 10(-5)). These data show that common, low-penetrance susceptibility alleles predispose to colorectal neoplasia.     

Meta-analysis of three genome-wide association studies identifies susceptibility loci for colorectal cancer at 1q41, 3q26.2, 12q13.13 and 20q13.33

Genome-wide association studies (GWAS) have identified ten loci harboring common variants that influence risk of developing colorectal cancer (CRC). To enhance the power to identify additional CRC risk loci, we conducted a meta-analysis of three GWAS from the UK which included a total of 3,334 affected individuals (cases) and 4,628 controls followed by multiple validation analyses including a total of 18,095 cases and 20,197 controls. We identified associations at four new CRC risk loci: 1q41 (rs6691170, odds ratio (OR) = 1.06, P = 9.55 × 10?1? and rs6687758, OR = 1.09, P = 2.27 × 10??, 3q26.2 (rs10936599, OR = 0.93, P = 3.39 × 10??), 12q13.13 (rs11169552, OR = 0.92, P = 1.89 × 10?1? and rs7136702, OR = 1.06, P = 4.02 × 10??) and 20q13.33 (rs4925386, OR = 0.93, P = 1.89 × 10?1?). In addition to identifying new CRC risk loci, this analysis provides evidence that additional CRC-associated variants of similar effect size remain to be discovered.     

A genome-wide association study shows that common alleles of SMAD7 influence colorectal cancer risk

To identify risk variants for colorectal cancer (CRC), we conducted a genome-wide association study, genotyping 550,163 tag SNPs in 940 individuals with familial colorectal tumor (627 CRC, 313 advanced adenomas) and 965 controls. We evaluated selected SNPs in three replication sample sets (7,473 cases, 5,984 controls) and identified three SNPs in SMAD7 (involved in TGF-beta and Wnt signaling) associated with CRC. Across the four sample sets, the association between rs4939827 and CRC was highly statistically significant (P(trend) = 1.0 x 10(-12)).     

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.     

Genome-wide association scan identifies a colorectal cancer susceptibility locus on chromosome 8q24

Using a multistage genetic association approach comprising 7,480 affected individuals and 7,779 controls, we identified markers in chromosomal region 8q24 associated with colorectal cancer. In stage 1, we genotyped 99,632 SNPs in 1,257 affected individuals and 1,336 controls from Ontario. In stages 2-4, we performed serial replication studies using 4,024 affected individuals and 4,042 controls from Seattle, Newfoundland and Scotland. We identified one locus on chromosome 8q24 and another on 9p24 having combined odds ratios (OR) for stages 1-4 of 1.18 (trend; P = 1.41 x 10(-8)) and 1.14 (trend; P = 1.32 x 10(-5)), respectively. Additional analyses in 2,199 affected individuals and 2,401 controls from France and Europe supported the association at the 8q24 locus (OR = 1.16, trend; 95% confidence interval (c.i.): 1.07-1.26; P = 5.05 x 10(-4)). A summary across all seven studies at the 8q24 locus was highly significant (OR = 1.17, c.i.: 1.12-1.23; P = 3.16 x 10(-11)). This locus has also been implicated in prostate cancer.     

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

Common variants at 11q12, 10q26 and 3p11.2 are associated with prostate cancer susceptibility in Japanese

We have previously reported multiple loci associated with prostate cancer susceptibility in a Japanese population using a genome-wide association study (GWAS). To identify additional prostate cancer susceptibility loci, we genotyped nine SNPs that were nominally associated with prostate cancer (P < 1 × 10(-4)) in our previous GWAS in three independent studies of prostate cancer in Japanese men (2,557 individuals with prostate cancer (cases) and 3,003 controls). In a meta-analysis of our previous GWAS and the replication studies, which included a total of 7,141 prostate cancer cases and 11,804 controls from a single ancestry group, three new loci reached genome-wide significance on chromosomes 11q12 (rs1938781; P = 1.10 × 10(-10); FAM111A-FAM111B), 10q26 (rs2252004; P = 1.98 × 10(-8)) and 3p11.2 (rs2055109; P = 3.94 × 10(-8)). We also found suggestive evidence of association at a previously reported prostate cancer susceptibility locus at 2p11 (rs2028898; P = 1.08 × 10(-7)). The identification of three new susceptibility loci should provide additional insight into the pathogenesis of prostate cancer and emphasizes the importance of conducting GWAS in diverse populations.     

Genome-wide association study of prostate cancer identifies a second risk locus at 8q24

Recently, common variants on human chromosome 8q24 were found to be associated with prostate cancer risk. While conducting a genome-wide association study in the Cancer Genetic Markers of Susceptibility project with 550,000 SNPs in a nested case-control study (1,172 cases and 1,157 controls of European origin), we identified a new association at 8q24 with an independent effect on prostate cancer susceptibility. The most significant signal is 70 kb centromeric to the previously reported SNP, rs1447295, but shows little evidence of linkage disequilibrium with it. A combined analysis with four additional studies (total: 4,296 cases and 4,299 controls) confirms association with prostate cancer for rs6983267 in the centromeric locus (P = 9.42 x 10(-13); heterozygote odds ratio (OR): 1.26, 95% confidence interval (c.i.): 1.13-1.41; homozygote OR: 1.58, 95% c.i.: 1.40-1.78). Each SNP remained significant in a joint analysis after adjusting for the other (rs1447295 P = 1.41 x 10(-11); rs6983267 P = 6.62 x 10(-10)). These observations, combined with compelling evidence for a recombination hotspot between the two markers, indicate the presence of at least two independent loci within 8q24 that contribute to prostate cancer in men of European ancestry. We estimate that the population attributable risk of the new locus, marked by rs6983267, is higher than the locus marked by rs1447295 (21% versus 9%).     

Genome-wide association scan of tag SNPs identifies a susceptibility locus for lung cancer at 15q25.1

To identify risk variants for lung cancer, we conducted a multistage genome-wide association study. In the discovery phase, we analyzed 315,450 tagging SNPs in 1,154 current and former (ever) smoking cases of European ancestry and 1,137 frequency-matched, ever-smoking controls from Houston, Texas. For replication, we evaluated the ten SNPs most significantly associated with lung cancer in an additional 711 cases and 632 controls from Texas and 2,013 cases and 3,062 controls from the UK. Two SNPs, rs1051730 and rs8034191, mapping to a region of strong linkage disequilibrium within 15q25.1 containing PSMA4 and the nicotinic acetylcholine receptor subunit genes CHRNA3 and CHRNA5, were significantly associated with risk in both replication sets. Combined analysis yielded odds ratios of 1.32 (P < 1 x 10(-17)) for both SNPs. Haplotype analysis was consistent with there being a single risk variant in this region. We conclude that variation in a region of 15q25.1 containing nicotinic acetylcholine receptors genes contributes to lung cancer risk.     

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 three loci associated with susceptibility to uterine fibroids

Uterine fibroids are a common benign tumor of the female genital tract. We conducted a genome-wide association study in which 457,044 SNPs were analyzed in 1,607 individuals with clinically diagnosed uterine fibroids and 1,428 female controls. SNPs showing suggestive associations (P < 5 × 10(-5)) were further genotyped in 3,466 additional cases and 3,245 female controls. Three loci on chromosomes 10q24.33, 22q13.1 and 11p15.5 revealed genome-wide significant associations with uterine fibroids. The SNPs showing the most significant association in a combination analysis at each of these loci were rs7913069 (P = 8.65 × 10(-14), odds ratio (OR) = 1.47), rs12484776 (P = 2.79 × 10(-12), OR = 1.23) and rs2280543 (P = 3.82 × 10(-12), OR = 1.39), respectively. Subsequent fine mapping of these regions will be necessary to pinpoint the causal variants. Our findings should shed light on the pathogenesis of uterine fibroids.     

Common variation near CDKN1A, POLD3 and SHROOM2 influences colorectal cancer risk

We performed a meta-analysis of five genome-wide association studies to identify common variants influencing colorectal cancer (CRC) risk comprising 8,682 cases and 9,649 controls. Replication analysis was performed in case-control sets totaling 21,096 cases and 19,555 controls. We identified three new CRC risk loci at 6p21 (rs1321311, near CDKN1A; P = 1.14 × 10(-10)), 11q13.4 (rs3824999, intronic to POLD3; P = 3.65 × 10(-10)) and Xp22.2 (rs5934683, near SHROOM2; P = 7.30 × 10(-10)) This brings the number of independent loci associated with CRC risk to 20 and provides further insight into the genetic architecture of inherited susceptibility to CRC.     

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

A genome-wide association study identifies six susceptibility loci for chronic lymphocytic leukemia

We conducted a genome-wide association study of 299,983 tagging SNPs for chronic lymphocytic leukemia (CLL) and performed validation in two additional series totaling 1,529 cases and 3,115 controls. We identified six previously unreported CLL risk loci at 2q13 (rs17483466; P = 2.36 x 10(-10)), 2q37.1 (rs13397985, SP140; P = 5.40 x 10(-10)), 6p25.3 (rs872071, IRF4; P = 1.91 x 10(-20)), 11q24.1 (rs735665; P = 3.78 x 10(-12)), 15q23 (rs7176508; P = 4.54 x 10(-12)) and 19q13.32 (rs11083846, PRKD2; P = 3.96 x 10(-9)). These data provide the first evidence for the existence of common, low-penetrance susceptibility to a hematological malignancy and new insights into disease causation in CLL.     

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

Sequence variants in IL10, ARPC2 and multiple other loci contribute to ulcerative colitis susceptibility

Inflammatory bowel disease (IBD) typically manifests as either ulcerative colitis (UC) or Crohn's disease (CD). Systematic identification of susceptibility genes for IBD has thus far focused mainly on CD, and little is known about the genetic architecture of UC. Here we report a genome-wide association study with 440,794 SNPs genotyped in 1,167 individuals with UC and 777 healthy controls. Twenty of the most significantly associated SNPs were tested for replication in three independent European case-control panels comprising a total of 1,855 individuals with UC and 3,091 controls. Among the four consistently replicated markers, SNP rs3024505 immediately flanking the IL10 (interleukin 10) gene on chromosome 1q32.1 showed the most significant association in the combined verification samples (P = 1.35 x 10(-12); OR = 1.46 (1.31-1.62)). The other markers were located in ARPC2 and in the HLA-BTNL2 region. Association between rs3024505 and CD (1,848 cases, 1,804 controls) was weak (P = 0.013; OR = 1.17 (1.01-1.34)). IL10 is an immunosuppressive cytokine that has long been proposed to influence IBD pathophysiology. Our findings strongly suggest that defective IL10 function is central to the pathogenesis of the UC subtype of IBD.     

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.     

Genome-wide association study identifies two new susceptibility loci for atopic dermatitis in the Chinese Han population

Atopic dermatitis is a chronic, relapsing form of inflammatory skin disorder that is affected by genetic and environmental factors. We performed a genome-wide association study of atopic dermatitis in a Chinese Han population using 1,012 affected individuals (cases) and 1,362 controls followed by a replication study in an additional 3,624 cases and 12,197 controls of Chinese Han ethnicity, as well as 1,806 cases and 3,256 controls from Germany. We identified previously undescribed susceptibility loci at 5q22.1 (TMEM232 and SLC25A46, rs7701890, P(combined) = 3.15 × 10(-9), odds ratio (OR) = 1.24) and 20q13.33 (TNFRSF6B and ZGPAT, rs6010620, P(combined) = 3.0 × 10(-8), OR = 1.17) and replicated another previously reported locus at 1q21.3 (FLG, rs3126085, P(combined) = 5.90 × 10(-12), OR = 0.82) in the Chinese sample. The 20q13.33 locus also showed evidence for association in the German sample (rs6010620, P = 2.87 × 10(-5), OR = 1.25). Our study identifies new genetic susceptibility factors and suggests previously unidentified biological pathways in atopic dermatitis.