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.     

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.     

PTPRD (protein tyrosine phosphatase receptor type delta) is associated with restless legs syndrome

We identified association of restless legs syndrome (RLS) with PTPRD at 9p23-24 in 2,458 affected individuals and 4,749 controls from Germany, Austria, Czechia and Canada. Two independent SNPs in the 5' UTR of splice variants expressed predominantly in the central nervous system showed highly significant P values (rs4626664, P(nominal/lambda corrected) = 5.91 x 10(-10), odds ratio (OR) = 1.44; rs1975197, P(nominal/lambda corrected) = 5.81 x 10(-9), OR = 1.31). This work identifies PTPRD as the fourth genome-wide significant locus for RLS.     

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.     

SNPs in KCNQ1 are associated with susceptibility to type 2 diabetes in East Asian and European populations

We conducted a genome-wide association study using 207,097 SNP markers in Japanese individuals with type 2 diabetes and unrelated controls, and identified KCNQ1 (potassium voltage-gated channel, KQT-like subfamily, member 1) to be a strong candidate for conferring susceptibility to type 2 diabetes. We detected consistent association of a SNP in KCNQ1 (rs2283228) with the disease in several independent case-control studies (additive model P = 3.1 x 10(-12); OR = 1.26, 95% CI = 1.18-1.34). Several other SNPs in the same linkage disequilibrium (LD) block were strongly associated with type 2 diabetes (additive model: rs2237895, P = 7.3 x 10(-9); OR = 1.32, 95% CI = 1.20-1.45, rs2237897, P = 6.8 x 10(-13); OR = 1.41, 95% CI = 1.29-1.55). The association of these SNPs with type 2 diabetes was replicated in samples from Singaporean (additive model: rs2237895, P = 8.5 x 10(-3); OR = 1.14, rs2237897, P = 2.4 x 10(-4); OR = 1.22) and Danish populations (additive model: rs2237895, P = 3.7 x 10(-11); OR = 1.24, rs2237897, P = 1.2 x 10(-4); OR = 1.36).     

A genome-wide association study in Han Chinese identifies multiple susceptibility loci for IgA nephropathy

We performed a two-stage genome-wide association study of IgA nephropathy (IgAN) in Han Chinese, with 1,434 affected individuals (cases) and 4,270 controls in the discovery phase and follow-up of the top 61 SNPs in an additional 2,703 cases and 3,464 controls. We identified associations at 17p13 (rs3803800, P = 9.40 × 10(-11), OR = 1.21; rs4227, P = 4.31 × 10(-10), OR = 1.23) and 8p23 (rs2738048, P = 3.18 × 10(-14), OR = 0.79) that implicated the genes encoding tumor necrosis factor (TNFSF13) and α-defensin (DEFA) as susceptibility genes. In addition, we found multiple associations in the major histocompatibility complex (MHC) region (rs660895, P = 4.13 × 10(-20), OR = 1.34; rs1794275, P = 3.43 × 10(-13), OR = 1.30; rs2523946, P = 1.74 × 10(-11), OR = 1.21) and confirmed a previously reported association at 22q12 (rs12537, P = 1.17 × 10(-11), OR = 0.78). We also found that rs660895 was associated with clinical subtypes of IgAN (P = 0.003), proteinuria (P = 0.025) and IgA levels (P = 0.047). Our findings show that IgAN is associated with variants near genes involved in innate immunity and inflammation.     

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.     

ASIP and TYR pigmentation variants associate with cutaneous melanoma and basal cell carcinoma

Fair color increases risk of cutaneous melanoma (CM) and basal cell carcinoma (BCC). Recent genome-wide association studies have identified variants affecting hair, eye and skin pigmentation in Europeans. Here, we assess the effect of these variants on risk of CM and BCC in European populations comprising 2,121 individuals with CM, 2,163 individuals with BCC and over 40,000 controls. A haplotype near ASIP, known to affect a similar spectrum of pigmentation traits as MC1R variants, conferred significant risk of CM (odds ratio (OR) = 1.45, P = 1.2 x 10(-9)) and BCC (OR = 1.33, P = 1.2 x 10(-6)). The variant in TYR encoding the R402Q amino acid substitution, previously shown to affect eye color and tanning response, conferred risk of CM (OR = 1.21, P = 2.8 x 10(-7)) and BCC (OR = 1.14, P = 6.1 x 10(-4)). An eye color variant in TYRP1 was associated with risk of CM (OR = 1.15, P = 4.6 x 10(-4)). The association of all three variants is robust with respect to adjustment for the effect of pigmentation.     

Genome-wide association study identifies a susceptibility locus for schizophrenia in Han Chinese at 11p11.2

To identify susceptibility loci for schizophrenia, we performed a two-stage genome-wide association study (GWAS) of schizophrenia in the Han Chinese population (GWAS: 746 individuals with schizophrenia and 1,599 healthy controls; validation: 4,027 individuals with schizophrenia and 5,603 healthy controls). We identified two susceptibility loci for schizophrenia at 6p21-p22.1 (rs1233710 in an intron of ZKSCAN4, P(combined) = 4.76 × 10(-11), odds ratio (OR) = 0.79; rs1635 in an exon of NKAPL, P(combined) = 6.91 × 10(-12), OR = 0.78; rs2142731 in an intron of PGBD1, P(combined) = 5.14 × 10(-10), OR = 0.79) and 11p11.2 (rs11038167 near the 5' UTR of TSPAN18, P(combined) = 1.09 × 10(-11), OR = 1.29; rs11038172, P(combined) = 7.21 × 10(-10), OR = 1.25; rs835784, P(combined) = 2.73 × 10(-11), OR = 1.27). These results add to previous evidence of susceptibility loci for schizophrenia at 6p21-p22.1 in the Han Chinese population. We found that NKAPL and ZKSCAN4 were expressed in postnatal day 0 (P0) mouse brain. These findings may lead to new insights into the pathogenesis of schizophrenia.     

Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21

In a genome-wide association study to identify loci associated with colorectal cancer (CRC) risk, we genotyped 555,510 SNPs in 1,012 early-onset Scottish CRC cases and 1,012 controls (phase 1). In phase 2, we genotyped the 15,008 highest-ranked SNPs in 2,057 Scottish cases and 2,111 controls. We then genotyped the five highest-ranked SNPs from the joint phase 1 and 2 analysis in 14,500 cases and 13,294 controls from seven populations, and identified a previously unreported association, rs3802842 on 11q23 (OR = 1.1; P = 5.8 x 10(-10)), showing population differences in risk. We also replicated and fine-mapped associations at 8q24 (rs7014346; OR = 1.19; P = 8.6 x 10(-26)) and 18q21 (rs4939827; OR = 1.2; P = 7.8 x 10(-28)). Risk was greater for rectal than for colon cancer for rs3802842 (P < 0.008) and rs4939827 (P < 0.009). Carrying all six possible risk alleles yielded OR = 2.6 (95% CI = 1.75-3.89) for CRC. These findings extend our understanding of the role of common genetic variation in CRC etiology.     

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.     

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.     

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

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.     

A genome-wide association study identifies new susceptibility loci for non-cardia gastric cancer at 3q13.31 and 5p13.1

Gastric cancer, including the cardia and non-cardia types, is the second leading cause of cancer-related deaths worldwide. To identify genetic risk variants for non-cardia gastric cancer, we performed a genome-wide association study in 3,279 individuals (1,006 with non-cardia gastric cancer and 2,273 controls) of Chinese descent. We replicated significant associations in an additional 6,897 subjects (3,288 with non-cardia gastric cancer and 3,609 controls). We identified two new susceptibility loci for non-cardia gastric cancer at 5p13.1 (rs13361707 in the region including PTGER4 and PRKAA1; odds ratio (OR) = 1.41; P = 7.6 × 10(-29)) and 3q13.31 (rs9841504 in ZBTB20; OR = 0.76; P = 1.7 × 10(-9)). Imputation analyses also confirmed previously reported associations of rs2294008 and rs2976392 on 8q24, rs4072037 on 1q22 and rs13042395 on 20p13 with non-cardia gastric cancer susceptibility in the Han Chinese population.     

Loci on 20q13 and 21q22 are associated with pediatric-onset inflammatory bowel disease

Inflammatory bowel disease (IBD) is a common inflammatory disorder with complex etiology that involves both genetic and environmental triggers, including but not limited to defects in bacterial clearance, defective mucosal barrier and persistent dysregulation of the immune response to commensal intestinal bacteria. IBD is characterized by two distinct phenotypes: Crohn's disease (CD) and ulcerative colitis (UC). Previously reported GWA studies have identified genetic variation accounting for a small portion of the overall genetic susceptibility to CD and an even smaller contribution to UC pathogenesis. We hypothesized that stratification of IBD by age of onset might identify additional genes associated with IBD. To that end, we carried out a GWA analysis in a cohort of 1,011 individuals with pediatric-onset IBD and 4,250 matched controls. We identified and replicated significantly associated, previously unreported loci on chromosomes 20q13 (rs2315008[T] and rs4809330[A]; P = 6.30 x 10(-8) and 6.95 x 10(-8), respectively; odds ratio (OR) = 0.74 for both) and 21q22 (rs2836878[A]; P = 6.01 x 10(-8); OR = 0.73), located close to the TNFRSF6B and PSMG1 genes, respectively.     

Chymotrypsin C (CTRC) variants that diminish activity or secretion are associated with chronic pancreatitis

Chronic pancreatitis is a persistent inflammatory disease of the pancreas, in which the digestive protease trypsin has a fundamental pathogenetic role. Here we have analyzed the gene encoding the trypsin-degrading enzyme chymotrypsin C (CTRC) in German subjects with idiopathic or hereditary chronic pancreatitis. Two alterations in this gene, p.R254W and p.K247_R254del, were significantly overrepresented in the pancreatitis group, being present in 30 of 901 (3.3%) affected individuals but only 21 of 2,804 (0.7%) controls (odds ratio (OR) = 4.6; confidence interval (CI) = 2.6-8.0; P = 1.3 x 10(-7)). A replication study identified these two variants in 10 of 348 (2.9%) individuals with alcoholic chronic pancreatitis but only 3 of 432 (0.7%) subjects with alcoholic liver disease (OR = 4.2; CI = 1.2-15.5; P = 0.02). CTRC variants were also found in 10 of 71 (14.1%) Indian subjects with tropical pancreatitis but only 1 of 84 (1.2%) healthy controls (OR = 13.6; CI = 1.7-109.2; P = 0.0028). Functional analysis of the CTRC variants showed impaired activity and/or reduced secretion. The results indicate that loss-of-function alterations in CTRC predispose to pancreatitis by diminishing its protective trypsin-degrading activity.     

A genome-wide association study identifies two new susceptibility loci for lung adenocarcinoma in the Japanese population

Lung adenocarcinoma is the most common histological type of lung cancer, and its incidence is increasing worldwide. To identify genetic factors influencing risk of lung adenocarcinoma, we conducted a genome-wide association study and two validation studies in the Japanese population comprising a total of 6,029 individuals with lung adenocarcinoma (cases) and 13,535 controls. We confirmed two previously reported risk loci, 5p15.33 (rs2853677, P(combined) = 2.8 × 10(-40), odds ratio (OR) = 1.41) and 3q28 (rs10937405, P(combined) = 6.9 × 10(-17), OR = 1.25), and identified two new susceptibility loci, 17q24.3 (rs7216064, P(combined) = 7.4 × 10(-11), OR = 1.20) and 6p21.3 (rs3817963, P(combined) = 2.7 × 10(-10), OR = 1.18). These data provide further evidence supporting a role for genetic susceptibility in the development of lung adenocarcinoma.     

Common variants near TARDBP and EGR2 are associated with susceptibility to Ewing sarcoma

Ewing sarcoma, a pediatric tumor characterized by EWSR1-ETS fusions, is predominantly observed in populations of European ancestry. We performed a genome-wide association study (GWAS) of 401 French individuals with Ewing sarcoma, 684 unaffected French individuals and 3,668 unaffected individuals of European descent and living in the United States. We identified candidate risk loci at 1p36.22, 10q21 and 15q15. We replicated these loci in two independent sets of cases and controls. Joint analysis identified associations with rs9430161 (P = 1.4 × 10(-20); odds ratio (OR) = 2.2) located 25 kb upstream of TARDBP, rs224278 (P = 4.0 × 10(-17); OR = 1.7) located 5 kb upstream of EGR2 and, to a lesser extent, rs4924410 at 15q15 (P = 6.6 × 10(-9); OR = 1.5). The major risk haplotypes were less prevalent in Africans, suggesting that these loci could contribute to geographical differences in Ewing sarcoma incidence. TARDBP shares structural similarities with EWSR1 and FUS, which encode RNA binding proteins, and EGR2 is a target gene of EWSR1-ETS. Variants at these loci were associated with expression levels of TARDBP, ADO (encoding cysteamine dioxygenase) and EGR2.