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.     

Genome-wide association study identifies new susceptibility loci for Crohn disease and implicates autophagy in disease pathogenesis

We present a genome-wide association study of ileal Crohn disease and two independent replication studies that identify several new regions of association to Crohn disease. Specifically, in addition to the previously established CARD15 and IL23R associations, we identified strong and significantly replicated associations (combined P < 10(-10)) with an intergenic region on 10q21.1 and a coding variant in ATG16L1, the latter of which was also recently reported by another group. We also report strong associations with independent replication to variation in the genomic regions encoding PHOX2B, NCF4 and a predicted gene on 16q24.1 (FAM92B). Finally, we demonstrate that ATG16L1 is expressed in intestinal epithelial cell lines and that functional knockdown of this gene abrogates autophagy of Salmonella typhimurium. Together, these findings suggest that autophagy and host cell responses to intracellular microbes are involved in the pathogenesis of Crohn disease.     

Genetic variation in the 5q31 cytokine gene cluster confers susceptibility to Crohn disease

Linkage disequilibrium (LD) mapping provides a powerful method for fine-structure localization of rare disease genes, but has not yet been widely applied to common disease. We sought to design a systematic approach for LD mapping and apply it to the localization of a gene (IBD5) conferring susceptibility to Crohn disease. The key issues are: (i) to detect a significant LD signal (ii) to rigorously bound the critical region and (iii) to identify the causal genetic variant within this region. We previously mapped the IBD5 locus to a large region spanning 18 cM of chromosome 5q31 (P<10(-4)). Using dense genetic maps of microsatellite markers and single-nucleotide polymorphisms (SNPs) across the entire region, we found strong evidence of LD. We bound the region to a common haplotype spanning 250 kb that shows strong association with the disease (P< 2 x 10(-7)) and contains the cytokine gene cluster. This finding provides overwhelming evidence that a specific common haplotype of the cytokine region in 5q31 confers susceptibility to Crohn disease. However, genetic evidence alone is not sufficient to identify the causal mutation within this region, as strong LD across the region results in multiple SNPs having equivalent genetic evidence-each consistent with the expected properties of the IBD5 locus. These results have important implications for Crohn disease in particular and LD mapping in general.     

Comprehensive analysis of the gene encoding filaggrin uncovers prevalent and rare mutations in ichthyosis vulgaris and atopic eczema

We recently reported two common filaggrin (FLG) null mutations that cause ichthyosis vulgaris and predispose to eczema and secondary allergic diseases. We show here that these common European mutations are ancestral variants carried on conserved haplotypes. To facilitate comprehensive analysis of other populations, we report a strategy for full sequencing of this large, highly repetitive gene, and we describe 15 variants, including seven that are prevalent. All the variants are either nonsense or frameshift mutations that, in representative cases, resulted in loss of filaggrin production in the epidermis. In an Irish case-control study, the five most common European mutations showed a strong association with moderate-to-severe childhood eczema (chi2 test: P = 2.12 x 10(-51); Fisher's exact test: heterozygote odds ratio (OR) = 7.44 (95% confidence interval (c.i.) = 4.9-11.3), and homozygote OR = 151 (95% c.i. = 20-1,136)). We found three additional rare null mutations in this case series, suggesting that the genetic architecture of filaggrin-related atopic dermatitis consists of both prevalent and rare risk alleles.     

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.     

A genome-wide association study identifies new psoriasis susceptibility loci and an interaction between HLA-C and ERAP1

To identify new susceptibility loci for psoriasis, we undertook a genome-wide association study of 594,224 SNPs in 2,622 individuals with psoriasis and 5,667 controls. We identified associations at eight previously unreported genomic loci. Seven loci harbored genes with recognized immune functions (IL28RA, REL, IFIH1, ERAP1, TRAF3IP2, NFKBIA and TYK2). These associations were replicated in 9,079 European samples (six loci with a combined P < 5 × 10?? and two loci with a combined P < 5 × 10??). We also report compelling evidence for an interaction between the HLA-C and ERAP1 loci (combined P = 6.95 × 10??). ERAP1 plays an important role in MHC class I peptide processing. ERAP1 variants only influenced psoriasis susceptibility in individuals carrying the HLA-C risk allele. Our findings implicate pathways that integrate epidermal barrier dysfunction with innate and adaptive immune dysregulation in psoriasis pathogenesis.     

Genetics of gene expression in primary immune cells identifies cell type-specific master regulators and roles of HLA alleles

Trans-acting genetic variants have a substantial, albeit poorly characterized, role in the heritable determination of gene expression. Using paired purified primary monocytes and B cells, we identify new predominantly cell type-specific cis and trans expression quantitative trait loci (eQTLs), including multi-locus trans associations to LYZ and KLF4 in monocytes and B cells, respectively. Additionally, we observe a B cell-specific trans association of rs11171739 at 12q13.2, a known autoimmune disease locus, with IP6K2 (P = 5.8 × 10(-15)), PRIC285 (P = 3.0 × 10(-10)) and an upstream region of CDKN1A (P = 2 × 10(-52)), suggesting roles for cell cycle regulation and peroxisome proliferator-activated receptor γ (PPARγ) signaling in autoimmune pathogenesis. We also find that specific human leukocyte antigen (HLA) alleles form trans associations with the expression of AOAH and ARHGAP24 in monocytes but not in B cells. In summary, we show that mapping gene expression in defined primary cell populations identifies new cell type-specific trans-regulated networks and provides insights into the genetic basis of disease susceptibility.     

Common variants at ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer's disease

We sought to identify new susceptibility loci for Alzheimer's disease through a staged association study (GERAD+) and by testing suggestive loci reported by the Alzheimer's Disease Genetic Consortium (ADGC) in a companion paper. We undertook a combined analysis of four genome-wide association datasets (stage 1) and identified ten newly associated variants with P ≤ 1 × 10(-5). We tested these variants for association in an independent sample (stage 2). Three SNPs at two loci replicated and showed evidence for association in a further sample (stage 3). Meta-analyses of all data provided compelling evidence that ABCA7 (rs3764650, meta P = 4.5 × 10(-17); including ADGC data, meta P = 5.0 × 10(-21)) and the MS4A gene cluster (rs610932, meta P = 1.8 × 10(-14); including ADGC data, meta P = 1.2 × 10(-16)) are new Alzheimer's disease susceptibility loci. We also found independent evidence for association for three loci reported by the ADGC, which, when combined, showed genome-wide significance: CD2AP (GERAD+, P = 8.0 × 10(-4); including ADGC data, meta P = 8.6 × 10(-9)), CD33 (GERAD+, P = 2.2 × 10(-4); including ADGC data, meta P = 1.6 × 10(-9)) and EPHA1 (GERAD+, P = 3.4 × 10(-4); including ADGC data, meta P = 6.0 × 10(-10)).     

Two independent alleles at 6q23 associated with risk of rheumatoid arthritis

To identify susceptibility alleles associated with rheumatoid arthritis, we genotyped 397 individuals with rheumatoid arthritis for 116,204 SNPs and carried out an association analysis in comparison to publicly available genotype data for 1,211 related individuals from the Framingham Heart Study. After evaluating and adjusting for technical and population biases, we identified a SNP at 6q23 (rs10499194, approximately 150 kb from TNFAIP3 and OLIG3) that was reproducibly associated with rheumatoid arthritis both in the genome-wide association (GWA) scan and in 5,541 additional case-control samples (P = 10(-3), GWA scan; P < 10(-6), replication; P = 10(-9), combined). In a concurrent study, the Wellcome Trust Case Control Consortium (WTCCC) has reported strong association of rheumatoid arthritis susceptibility to a different SNP located 3.8 kb from rs10499194 (rs6920220; P = 5 x 10(-6) in WTCCC). We show that these two SNP associations are statistically independent, are each reproducible in the comparison of our data and WTCCC data, and define risk and protective haplotypes for rheumatoid arthritis at 6q23.     

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.     

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.     

A genome-wide association study identifies susceptibility loci for Wilms tumor

Wilms tumor is the most common renal malignancy of childhood. To identify common variants that confer susceptibility to Wilms tumor, we conducted a genome-wide association study in 757 individuals with Wilms tumor (cases) and 1,879 controls. We evaluated ten SNPs in regions significantly associated at P < 5 × 10(-5) in two independent replication series from the UK (769 cases and 2,814 controls) and the United States (719 cases and 1,037 controls). We identified clear significant associations at 2p24 (rs3755132, P = 1.03 × 10(-14); rs807624, P = 1.32 × 10(-14)) and 11q14 (rs790356, P = 4.25 × 10(-15)). Both regions contain genes that are plausibly related to Wilms tumorigenesis. We also identified candidate association signals at 5q14, 22q12 and Xp22.     

CFH haplotypes without the Y402H coding variant show strong association with susceptibility to age-related macular degeneration

In developed countries, age-related macular degeneration is a common cause of blindness in the elderly. A common polymorphism, encoding the sequence variation Y402H in complement factor H (CFH), has been strongly associated with disease susceptibility. Here, we examined 84 polymorphisms in and around CFH in 726 affected individuals (including 544 unrelated individuals) and 268 unrelated controls. In this sample, 20 of these polymorphisms showed stronger association with disease susceptibility than the Y402H variant. Further, no single polymorphism could account for the contribution of the CFH locus to disease susceptibility. Instead, multiple polymorphisms defined a set of four common haplotypes (of which two were associated with disease susceptibility and two seemed to be protective) and multiple rare haplotypes (associated with increased susceptibility in aggregate). Our results suggest that there are multiple disease susceptibility alleles in the region and that noncoding CFH variants play a role in disease susceptibility.     

Common variation in GPC5 is associated with acquired nephrotic syndrome

Severe proteinuria is a defining factor of nephrotic syndrome irrespective of the etiology. Investigation of congenital nephrotic syndrome has shown that dysfunction of glomerular epithelial cells (podocytes) plays a crucial role in this disease. Acquired nephrotic syndrome is also assumed to be associated with podocyte injury. Here we identify an association between variants in GPC5, encoding glypican-5, and acquired nephrotic syndrome through a genome-wide association study and replication analysis (P value under a recessive model (P(rec)) = 6.0 × 10(-11), odds ratio = 2.54). We show that GPC5 is expressed in podocytes and that the risk genotype is associated with higher expression. We further show that podocyte-specific knockdown and systemic short interfering RNA injection confers resistance to podocyte injury in mouse models of nephrosis. This study identifies GPC5 as a new susceptibility gene for nephrotic syndrome and implicates GPC5 as a promising therapeutic target for reducing podocyte vulnerability in glomerular disease.     

Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility

Ankylosing spondylitis is a common form of inflammatory arthritis predominantly affecting the spine and pelvis that occurs in approximately 5 out of 1,000 adults of European descent. Here we report the identification of three variants in the RUNX3, LTBR-TNFRSF1A and IL12B regions convincingly associated with ankylosing spondylitis (P < 5 × 10(-8) in the combined discovery and replication datasets) and a further four loci at PTGER4, TBKBP1, ANTXR2 and CARD9 that show strong association across all our datasets (P < 5 × 10(-6) overall, with support in each of the three datasets studied). We also show that polymorphisms of ERAP1, which encodes an endoplasmic reticulum aminopeptidase involved in peptide trimming before HLA class I presentation, only affect ankylosing spondylitis risk in HLA-B27-positive individuals. These findings provide strong evidence that HLA-B27 operates in ankylosing spondylitis through a mechanism involving aberrant processing of antigenic peptides.     

Genome-wide association analysis identifies 20 loci that influence adult height

Adult height is a model polygenic trait, but there has been limited success in identifying the genes underlying its normal variation. To identify genetic variants influencing adult human height, we used genome-wide association data from 13,665 individuals and genotyped 39 variants in an additional 16,482 samples. We identified 20 variants associated with adult height (P < 5 x 10(-7), with 10 reaching P < 1 x 10(-10)). Combined, the 20 SNPs explain approximately 3% of height variation, with a approximately 5 cm difference between the 6.2% of people with 17 or fewer 'tall' alleles compared to the 5.5% with 27 or more 'tall' alleles. The loci we identified implicate genes in Hedgehog signaling (IHH, HHIP, PTCH1), extracellular matrix (EFEMP1, ADAMTSL3, ACAN) and cancer (CDK6, HMGA2, DLEU7) pathways, and provide new insights into human growth and developmental processes. Finally, our results provide insights into the genetic architecture of a classic quantitative trait.     

Chromosome 6q25 is linked to susceptibility to leprosy in a Vietnamese population

Leprosy, a chronic infectious disease caused by Mycobacterium leprae, affects an estimated 700,000 persons each year. Clinically, leprosy can be categorized as paucibacillary or multibacillary disease. These clinical forms develop in persons that are intrinsically susceptible to leprosy per se, that is, leprosy independent of its specific clinical manifestation. We report here on a genome-wide search for loci controlling susceptibility to leprosy per se in a panel of 86 families including 205 siblings affected with leprosy from Southern Vietnam. Using model-free linkage analysis, we found significant evidence for a susceptibility gene on chromosome region 6q25 (maximum likelihood binomial (MLB) lod score 4.31; P = 5 x 10(-6)). We confirmed this by family-based association analysis in an independent panel of 208 Vietnamese leprosy simplex families. Of seven microsatellite markers underlying the linkage peak, alleles of two markers (D6S1035 and D6S305) showed strong evidence for association with leprosy (P = 6.7 x 10(-4) and P = 5.9 x 10(-5), respectively).     

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.