Bayesian inference analyses of the polygenic architecture of rheumatoid arthritis

The genetic architectures of common, complex diseases are largely uncharacterized. We modeled the genetic architecture underlying genome-wide association study (GWAS) data for rheumatoid arthritis and developed a new method using polygenic risk-score analyses to infer the total liability-scale variance explained by associated GWAS SNPs. Using this method, we estimated that, together, thousands of SNPs from rheumatoid arthritis GWAS explain an additional 20% of disease risk (excluding known associated loci). We further tested this method on datasets for three additional diseases and obtained comparable estimates for celiac disease (43% excluding the major histocompatibility complex), myocardial infarction and coronary artery disease (48%) and type 2 diabetes (49%). Our results are consistent with simulated genetic models in which hundreds of associated loci harbor common causal variants and a smaller number of loci harbor multiple rare causal variants. These analyses suggest that GWAS will continue to be highly productive for the discovery of additional susceptibility loci for common diseases.     

Newly identified genetic risk variants for celiac disease related to the immune response

Our genome-wide association study of celiac disease previously identified risk variants in the IL2-IL21 region. To identify additional risk variants, we genotyped 1,020 of the most strongly associated non-HLA markers in an additional 1,643 cases and 3,406 controls. Through joint analysis including the genome-wide association study data (767 cases, 1,422 controls), we identified seven previously unknown risk regions (P < 5 x 10(-7)). Six regions harbor genes controlling immune responses, including CCR3, IL12A, IL18RAP, RGS1, SH2B3 (nsSNP rs3184504) and TAGAP. Whole-blood IL18RAP mRNA expression correlated with IL18RAP genotype. Type 1 diabetes and celiac disease share HLA-DQ, IL2-IL21, CCR3 and SH2B3 risk regions. Thus, this extensive genome-wide association follow-up study has identified additional celiac disease risk variants in relevant biological pathways.     

Newly identified loci that influence lipid concentrations and risk of coronary artery disease

To identify genetic variants influencing plasma lipid concentrations, we first used genotype imputation and meta-analysis to combine three genome-wide scans totaling 8,816 individuals and comprising 6,068 individuals specific to our study (1,874 individuals from the FUSION study of type 2 diabetes and 4,184 individuals from the SardiNIA study of aging-associated variables) and 2,758 individuals from the Diabetes Genetics Initiative, reported in a companion study in this issue. We subsequently examined promising signals in 11,569 additional individuals. Overall, we identify strongly associated variants in eleven loci previously implicated in lipid metabolism (ABCA1, the APOA5-APOA4-APOC3-APOA1 and APOE-APOC clusters, APOB, CETP, GCKR, LDLR, LPL, LIPC, LIPG and PCSK9) and also in several newly identified loci (near MVK-MMAB and GALNT2, with variants primarily associated with high-density lipoprotein (HDL) cholesterol; near SORT1, with variants primarily associated with low-density lipoprotein (LDL) cholesterol; near TRIB1, MLXIPL and ANGPTL3, with variants primarily associated with triglycerides; and a locus encompassing several genes near NCAN, with variants strongly associated with both triglycerides and LDL cholesterol). Notably, the 11 independent variants associated with increased LDL cholesterol concentrations in our study also showed increased frequency in a sample of coronary artery disease cases versus controls.     

Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease

We performed a meta-analysis of 14 genome-wide association studies of coronary artery disease (CAD) comprising 22,233 individuals with CAD (cases) and 64,762 controls of European descent followed by genotyping of top association signals in 56,682 additional individuals. This analysis identified 13 loci newly associated with CAD at P < 5 × 10?? and confirmed the association of 10 of 12 previously reported CAD loci. The 13 new loci showed risk allele frequencies ranging from 0.13 to 0.91 and were associated with a 6% to 17% increase in the risk of CAD per allele. Notably, only three of the new loci showed significant association with traditional CAD risk factors and the majority lie in gene regions not previously implicated in the pathogenesis of CAD. Finally, five of the new CAD risk loci appear to have pleiotropic effects, showing strong association with various other human diseases or traits.     

Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47

Genome-wide association studies and candidate gene studies in ulcerative colitis have identified 18 susceptibility loci. We conducted a meta-analysis of six ulcerative colitis genome-wide association study datasets, comprising 6,687 cases and 19,718 controls, and followed up the top association signals in 9,628 cases and 12,917 controls. We identified 29 additional risk loci (P < 5 × 10(-8)), increasing the number of ulcerative colitis-associated loci to 47. After annotating associated regions using GRAIL, expression quantitative trait loci data and correlations with non-synonymous SNPs, we identified many candidate genes that provide potentially important insights into disease pathogenesis, including IL1R2, IL8RA-IL8RB, IL7R, IL12B, DAP, PRDM1, JAK2, IRF5, GNA12 and LSP1. The total number of confirmed inflammatory bowel disease risk loci is now 99, including a minimum of 28 shared association signals between Crohn's disease and ulcerative colitis.     

Genome-wide association identifies three new susceptibility loci for Paget's disease of bone

Paget's disease of bone (PDB) is a common disorder characterized by focal abnormalities of bone remodeling. We previously identified variants at the CSF1, OPTN and TNFRSF11A loci as risk factors for PDB by genome-wide association study. Here we extended this study, identified three new loci and confirmed their association with PDB in 2,215 affected individuals (cases) and 4,370 controls from seven independent populations. The new associations were with rs5742915 within PML on 15q24 (odds ratio (OR) = 1.34, P = 1.6 × 10(-14)), rs10498635 within RIN3 on 14q32 (OR = 1.44, P = 2.55 × 10(-11)) and rs4294134 within NUP205 on 7q33 (OR = 1.45, P = 8.45 × 10(-10)). Our data also confirmed the association of TM7SF4 (rs2458413, OR = 1.40, P = 7.38 × 10(-17)) with PDB. These seven loci explained ～13% of the familial risk of PDB. These studies provide new insights into the genetic architecture and pathophysiology of PDB.     

Resequencing of positional candidates identifies low frequency IL23R coding variants protecting against inflammatory bowel disease

Genome-wide association studies (GWAS) have identified dozens of risk loci for many complex disorders, including Crohn's disease. However, common disease-associated SNPs explain at most ～20% of the genetic variance for Crohn's disease. Several factors may account for this unexplained heritability, including rare risk variants not adequately tagged thus far in GWAS. That rare susceptibility variants indeed contribute to variation in multifactorial phenotypes has been demonstrated for colorectal cancer, plasma high-density lipoprotein cholesterol levels, blood pressure, type 1 diabetes, hypertriglyceridemia and, in the case of Crohn's disease, for NOD2 (refs. 14,15). Here we describe the use of high-throughput resequencing of DNA pools to search for rare coding variants influencing susceptibility to Crohn's disease in 63 GWAS-identified positional candidate genes. We identify low frequency coding variants conferring protection against inflammatory bowel disease in IL23R, but we conclude that rare coding variants in positional candidates do not make a large contribution to inherited predisposition to Crohn's disease.     

Seven prostate cancer susceptibility loci identified by a multi-stage genome-wide association study

Prostate cancer (PrCa) is the most frequently diagnosed male cancer in developed countries. We conducted a multi-stage genome-wide association study for PrCa and previously reported the results of the first two stages, which identified 16 PrCa susceptibility loci. We report here the results of stage 3, in which we evaluated 1,536 SNPs in 4,574 individuals with prostate cancer (cases) and 4,164 controls. We followed up ten new association signals through genotyping in 51,311 samples in 30 studies from the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (PRACTICAL) consortium. In addition to replicating previously reported loci, we identified seven new prostate cancer susceptibility loci on chromosomes 2p11, 3q23, 3q26, 5p12, 6p21, 12q13 and Xq12 (P = 4.0 × 10(-8) to P = 2.7 × 10(-24)). We also identified a SNP in TERT more strongly associated with PrCa than that previously reported. More than 40 PrCa susceptibility loci, explaining ～25% of the familial risk in this disease, have now been identified.     

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohn's disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohn's disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10(-16), odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.     

Meta-analysis of genome-wide association studies of asthma in ethnically diverse North American populations

Asthma is a common disease with a complex risk architecture including both genetic and environmental factors. We performed a meta-analysis of North American genome-wide association studies of asthma in 5,416 individuals with asthma (cases) including individuals of European American, African American or African Caribbean, and Latino ancestry, with replication in an additional 12,649 individuals from the same ethnic groups. We identified five susceptibility loci. Four were at previously reported loci on 17q21, near IL1RL1, TSLP and IL33, but we report for the first time, to our knowledge, that these loci are associated with asthma risk in three ethnic groups. In addition, we identified a new asthma susceptibility locus at PYHIN1, with the association being specific to individuals of African descent (P = 3.9 × 10(-9)). These results suggest that some asthma susceptibility loci are robust to differences in ancestry when sufficiently large samples sizes are investigated, and that ancestry-specific associations also contribute to the complex genetic architecture of asthma.     

Genome-wide association defines more than 30 distinct susceptibility loci for Crohn's disease

Several risk factors for Crohn's disease have been identified in recent genome-wide association studies. To advance gene discovery further, we combined data from three studies on Crohn's disease (a total of 3,230 cases and 4,829 controls) and carried out replication in 3,664 independent cases with a mixture of population-based and family-based controls. The results strongly confirm 11 previously reported loci and provide genome-wide significant evidence for 21 additional loci, including the regions containing STAT3, JAK2, ICOSLG, CDKAL1 and ITLN1. The expanded molecular understanding of the basis of this disease offers promise for informed therapeutic development.     

A genome-wide association study identifies two new risk loci for Graves' disease

Graves' disease is a common autoimmune disorder characterized by thyroid stimulating hormone receptor autoantibodies (TRAb) and hyperthyroidism. To investigate the genetic architecture of Graves' disease, we conducted a genome-wide association study in 1,536 individuals with Graves' disease (cases) and 1,516 controls. We further evaluated a group of associated SNPs in a second set of 3,994 cases and 3,510 controls. We confirmed four previously reported loci (in the major histocompatibility complex, TSHR, CTLA4 and FCRL3) and identified two new susceptibility loci (the RNASET2-FGFR1OP-CCR6 region at 6q27 (P(combined) = 6.85 × 10(-10) for rs9355610) and an intergenic region at 4p14 (P(combined) = 1.08 × 10(-13) for rs6832151)). These newly associated SNPs were correlated with the expression levels of RNASET2 at 6q27, of CHRNA9 and of a previously uncharacterized gene at 4p14, respectively. Moreover, we identified strong associations of TSHR and major histocompatibility complex class II variants with persistently TRAb-positive Graves' disease.     

Genome-wide association study identifies three new susceptibility loci for esophageal squamous-cell carcinoma in Chinese populations

Esophageal squamous-cell carcinoma (ESCC) is one of the most prevalent cancers worldwide and occurs at a relatively high frequency in China. To identify genetic susceptibility loci for ESCC, we conducted a genome-wide association study on 2,031 individuals with ESCC (cases) and 2,044 controls of Chinese descent using 666,141 autosomal SNPs. We evaluated promising associations in an additional 6,276 cases and 6,165 controls of Chinese descent from different areas of China. We identified seven susceptibility loci on chromosomes 5q11, 6p21, 10q23, 12q24 and 21q22 (ranging from P = 7.48 × 10(-12) to P = 2.44 × 10(-31)); among these loci, 5q11, 6p21 and 21q22 were newly identified. Three variants in high linkage disequilibrium on 12q24 confer their risks to ESCC in a gene-lifestyle interaction manner, with more pronounced risk enhancement seen in tobacco and alcohol users. Furthermore, the identified variants had a cumulative association with ESCC risk (P(trend) = 7.92 × 10(-56)). These findings highlight the involvement of multiple genetic loci and gene-environment interaction in the development of esophageal cancer.     

Genome-wide association study identifies 12 new susceptibility loci for primary biliary cirrhosis

In addition to the HLA locus, six genetic risk factors for primary biliary cirrhosis (PBC) have been identified in recent genome-wide association studies (GWAS). To identify additional loci, we carried out a GWAS using 1,840 cases from the UK PBC Consortium and 5,163 UK population controls as part of the Wellcome Trust Case Control Consortium 3 (WTCCC3). We followed up 28 loci in an additional UK cohort of 620 PBC cases and 2,514 population controls. We identified 12 new susceptibility loci (at a genome-wide significance level of P < 5 × 10??) and replicated all previously associated loci. We identified three further new loci in a meta-analysis of data from our study and previously published GWAS results. New candidate genes include STAT4, DENND1B, CD80, IL7R, CXCR5, TNFRSF1A, CLEC16A and NFKB1. This study has considerably expanded our knowledge of the genetic architecture of PBC.     

Identification of common variants influencing risk of the tauopathy progressive supranuclear palsy

Progressive supranuclear palsy (PSP) is a movement disorder with prominent tau neuropathology. Brain diseases with abnormal tau deposits are called tauopathies, the most common of which is Alzheimer's disease. Environmental causes of tauopathies include repetitive head trauma associated with some sports. To identify common genetic variation contributing to risk for tauopathies, we carried out a genome-wide association study of 1,114 individuals with PSP (cases) and 3,247 controls (stage 1) followed by a second stage in which we genotyped 1,051 cases and 3,560 controls for the stage 1 SNPs that yielded P ≤ 10(-3). We found significant previously unidentified signals (P < 5 × 10(-8)) associated with PSP risk at STX6, EIF2AK3 and MOBP. We confirmed two independent variants in MAPT affecting risk for PSP, one of which influences MAPT brain expression. The genes implicated encode proteins for vesicle-membrane fusion at the Golgi-endosomal interface, for the endoplasmic reticulum unfolded protein response and for a myelin structural component.     

A comprehensive linkage analysis for myocardial infarction and its related risk factors

Coronary artery disease and myocardial infarction (MI) are leading causes of death in the western world. Numerous studies have shown that risk factors such as diabetes mellitus, arterial hypertension and hypercholesterolemia contribute to the development of the disease. Although each risk factor by itself is partly under genetic control, a positive family history is an independent predictor, which suggests that there are additional susceptibility genes. We have scanned the whole genome in 513 families to identify chromosomal regions linked to myocardial infarction and related risk factors that are known to be under genetic control. Here we show, by using variance component analysis and incorporating risk factors, that risk of myocardial infarction maps to a single region on chromosome 14 with a significant lod score of 3.9 (pointwise P=0.00015, genome-wide P<0.05), providing evidence of a principal MI locus. To characterize this locus we analyzed each risk factor by itself. Serum concentrations of lipoprotein (a) show linkage to both the apolipoprotein (a) locus (lod score 26.99) and a new locus on chromosome 1 (lod score 3.8). There is suggestive linkage for diabetes mellitus on chromosome 6 (lod score 2.96), for hypertension on chromosomes 1 and 6, for high-density and low-density lipoprotein cholesterol on chromosomes 1 and 17, and for triglyceride concentrations on chromosome 9. Although some of these risk factors overlap with previously identified loci, none overlaps with the newly identified susceptibility locus for myocardial infarction and coronary artery disease.     

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.     

Detection of large-scale variation in the human genome

We identified 255 loci across the human genome that contain genomic imbalances among unrelated individuals. Twenty-four variants are present in > 10% of the individuals that we examined. Half of these regions overlap with genes, and many coincide with segmental duplications or gaps in the human genome assembly. This previously unappreciated heterogeneity may underlie certain human phenotypic variation and susceptibility to disease and argues for a more dynamic human genome structure.     

Genome-wide meta-analyses of nonsyndromic cleft lip with or without cleft palate identify six new risk loci

We have conducted the first meta-analyses for nonsyndromic cleft lip with or without cleft palate (NSCL/P) using data from the two largest genome-wide association studies published to date. We confirmed associations with all previously identified loci and identified six additional susceptibility regions (1p36, 2p21, 3p11.1, 8q21.3, 13q31.1 and 15q22). Analysis of phenotypic variability identified the first specific genetic risk factor for NSCLP (nonsyndromic cleft lip plus palate) (rs8001641; P(NSCLP) = 6.51 × 10(-11); homozygote relative risk = 2.41, 95% confidence interval (CI) 1.84-3.16).