Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM, PXK, KIAA1542 and other loci

Systemic lupus erythematosus (SLE) is a common systemic autoimmune disease with complex etiology but strong clustering in families (lambda(S) = approximately 30). We performed a genome-wide association scan using 317,501 SNPs in 720 women of European ancestry with SLE and in 2,337 controls, and we genotyped consistently associated SNPs in two additional independent sample sets totaling 1,846 affected women and 1,825 controls. Aside from the expected strong association between SLE and the HLA region on chromosome 6p21 and the previously confirmed non-HLA locus IRF5 on chromosome 7q32, we found evidence of association with replication (1.1 x 10(-7) < P(overall) < 1.6 x 10(-23); odds ratio = 0.82-1.62) in four regions: 16p11.2 (ITGAM), 11p15.5 (KIAA1542), 3p14.3 (PXK) and 1q25.1 (rs10798269). We also found evidence for association (P < 1 x 10(-5)) at FCGR2A, PTPN22 and STAT4, regions previously associated with SLE and other autoimmune diseases, as well as at > or =9 other loci (P < 2 x 10(-7)). Our results show that numerous genes, some with known immune-related functions, predispose to SLE.     

Meta-analysis of genome-wide association data and large-scale replication identifies additional susceptibility loci for type 2 diabetes

Genome-wide association (GWA) studies have identified multiple loci at which common variants modestly but reproducibly influence risk of type 2 diabetes (T2D). Established associations to common and rare variants explain only a small proportion of the heritability of T2D. As previously published analyses had limited power to identify variants with modest effects, we carried out meta-analysis of three T2D GWA scans comprising 10,128 individuals of European descent and approximately 2.2 million SNPs (directly genotyped and imputed), followed by replication testing in an independent sample with an effective sample size of up to 53,975. We detected at least six previously unknown loci with robust evidence for association, including the JAZF1 (P = 5.0 x 10(-14)), CDC123-CAMK1D (P = 1.2 x 10(-10)), TSPAN8-LGR5 (P = 1.1 x 10(-9)), THADA (P = 1.1 x 10(-9)), ADAMTS9 (P = 1.2 x 10(-8)) and NOTCH2 (P = 4.1 x 10(-8)) gene regions. Our results illustrate the value of large discovery and follow-up samples for gaining further insights into the inherited basis of T2D.     

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 common variant of HMGA2 is associated with adult and childhood height in the general population

Human height is a classic, highly heritable quantitative trait. To begin to identify genetic variants influencing height, we examined genome-wide association data from 4,921 individuals. Common variants in the HMGA2 oncogene, exemplified by rs1042725, were associated with height (P = 4 x 10(-8)). HMGA2 is also a strong biological candidate for height, as rare, severe mutations in this gene alter body size in mice and humans, so we tested rs1042725 in additional samples. We confirmed the association in 19,064 adults from four further studies (P = 3 x 10(-11), overall P = 4 x 10(-16), including the genome-wide association data). We also observed the association in children (P = 1 x 10(-6), N = 6,827) and a tall/short case-control study (P = 4 x 10(-6), N = 3,207). We estimate that rs1042725 explains approximately 0.3% of population variation in height (approximately 0.4 cm increased adult height per C allele). There are few examples of common genetic variants reproducibly associated with human quantitativetraits; these results represent, to our knowledge, the first consistently replicated association with adult and childhood height.     

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

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.     

A genome-wide association scan identifies the hepatic cholesterol transporter ABCG8 as a susceptibility factor for human gallstone disease

With an overall prevalence of 10-20%, gallstone disease (cholelithiasis) represents one of the most frequent and economically relevant health problems of industrialized countries. We performed an association scan of >500,000 SNPs in 280 individuals with gallstones and 360 controls. A follow-up study of the 235 most significant SNPs in 1,105 affected individuals and 873 controls replicated the disease association of SNP A-1791411 in ABCG8 (allelic P value P(CCA) = 4.1 x 10(-9)), which was subsequently attributed to coding variant rs11887534 (D19H). Additional replication was achieved in 728 German (P = 2.8 x 10(-7)) and 167 Chilean subjects (P = 0.02). The overall odds ratio for D19H carriership was 2.2 (95% confidence interval: 1.8-2.6, P = 1.4 x 10(-14)) in the full German sample. Association was stronger in subjects with cholesterol gallstones (odds ratio = 3.3), suggesting that His19 might be associated with a more efficient transport of cholesterol into the bile.     

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.     

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.     

Lung cancer susceptibility locus at 5p15.33

We carried out a genome-wide association study of lung cancer (3,259 cases and 4,159 controls), followed by replication in 2,899 cases and 5,573 controls. Two uncorrelated disease markers at 5p15.33, rs402710 and rs2736100 were detected by the genome-wide data (P = 2 x 10(-7) and P = 4 x 10(-6)) and replicated by the independent study series (P = 7 x 10(-5) and P = 0.016). The susceptibility region contains two genes, TERT and CLPTM1L, suggesting that one or both may have a role in lung cancer etiology.     

Genome-wide association study identifies five loci associated with susceptibility to pancreatic cancer in Chinese populations

Pancreatic cancer has the lowest survival rate among human cancers, and there are no effective markers for its screening and early diagnosis. To identify genetic susceptibility markers for this cancer, we carried out a genome-wide association study on 981 individuals with pancreatic cancer (cases) and 1,991 cancer-free controls of Chinese descent using 666,141 autosomal SNPs. Promising associations were replicated in an additional 2,603 pancreatic cancer cases and 2,877 controls recruited from 25 hospitals in 16 provinces or cities in China. We identified five new susceptibility loci at chromosomes 21q21.3, 5p13.1, 21q22.3, 22q13.32 and 10q26.11 (P = 2.24 × 10(-13) to P = 4.18 × 10(-10)) in addition to 13q22.1 previously reported in populations of European ancestry. These results advance our understanding of the development of pancreatic cancer and highlight potential targets for the prevention or treatment of this cancer.     

Variants in KCNQ1 are associated with susceptibility to type 2 diabetes mellitus

We carried out a multistage genome-wide association study of type 2 diabetes mellitus in Japanese individuals, with a total of 1,612 cases and 1,424 controls and 100,000 SNPs. The most significant association was obtained with SNPs in KCNQ1, and dense mapping within the gene revealed that rs2237892 in intron 15 showed the lowest Pvalue (6.7 x 10(-13), odds ratio (OR) = 1.49). The association of KCNQ1 with type 2 diabetes was replicated in populations of Korean, Chinese and European ancestry as well as in two independent Japanese populations, and meta-analysis with a total of 19,930 individuals (9,569 cases and 10,361 controls) yielded a P value of 1.7 x 10(-42) (OR = 1.40; 95% CI = 1.34-1.47) for rs2237892. Among control subjects, the risk allele of this polymorphism was associated with impairment of insulin secretion according to the homeostasis model assessment of beta-cell function or the corrected insulin response. Our data thus implicate KCNQ1 as a diabetes susceptibility gene in groups of different ancestries.     

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.     

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.     

Genome-wide association study identifies three new melanoma susceptibility loci

We report a genome-wide association study for melanoma that was conducted by the GenoMEL Consortium. Our discovery phase included 2,981 individuals with melanoma and 1,982 study-specific control individuals of European ancestry, as well as an additional 6,426 control subjects from French or British populations, all of whom were genotyped for 317,000 or 610,000 single-nucleotide polymorphisms (SNPs). Our analysis replicated previously known melanoma susceptibility loci. Seven new regions with at least one SNP with P < 10(-5) and further local imputed or genotyped support were selected for replication using two other genome-wide studies (from Australia and Texas, USA). Additional replication came from case-control series from the UK and The Netherlands. Variants at three of the seven loci replicated at P < 10(-3): an SNP in ATM (rs1801516, overall P = 3.4 × 10(-9)), an SNP in MX2 (rs45430, P = 2.9 × 10(-9)) and an SNP adjacent to CASP8 (rs13016963, P = 8.6 × 10(-10)). A fourth locus near CCND1 remains of potential interest, showing suggestive but inconclusive evidence of replication (rs1485993, overall P = 4.6 × 10(-7) under a fixed-effects model and P = 1.2 × 10(-3) under a random-effects model). These newly associated variants showed no association with nevus or pigmentation phenotypes in a large British case-control series.     

A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1

We performed a genome-wide association study of 19,779 nonsynonymous SNPs in 735 individuals with Crohn disease and 368 controls. A total of 7,159 of these SNPs were informative. We followed up on all 72 SNPs with P 0.4), these data suggest that the underlying biological process may be specific to Crohn disease.