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.     

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.     

Genome-wide association study in Han Chinese identifies four new susceptibility loci for coronary artery disease

We performed a meta-analysis of 2 genome-wide association studies of coronary artery disease comprising 1,515 cases and 5,019 controls followed by replication studies in 15,460 cases and 11,472 controls, all of Chinese Han ancestry. We identify four new loci for coronary artery disease that reached the threshold of genome-wide significance (P < 5 × 10(-8)). These loci mapped in or near TTC32-WDR35, GUCY1A3, C6orf10-BTNL2 and ATP2B1. We also replicated four loci previously identified in European populations (in or near PHACTR1, TCF21, CDKN2A-CDKN2B and C12orf51). These findings provide new insights into pathways contributing to the susceptibility for coronary artery disease in the Chinese Han population.     

A genome-wide association study in Europeans and South Asians identifies five new loci for coronary artery disease

Genome-wide association studies have identified 11 common variants convincingly associated with coronary artery disease (CAD)1??, a modest number considering the apparent heritability of CAD?. All of these variants have been discovered in European populations. We report a meta-analysis of four large genome-wide association studies of CAD, with ～575,000 genotyped SNPs in a discovery dataset comprising 15,420 individuals with CAD (cases) (8,424 Europeans and 6,996 South Asians) and 15,062 controls. There was little evidence for ancestry-specific associations, supporting the use of combined analyses. Replication in an independent sample of 21,408 cases and 19,185 controls identified five loci newly associated with CAD (P < 5 × 10?? in the combined discovery and replication analysis): LIPA on 10q23, PDGFD on 11q22, ADAMTS7-MORF4L1 on 15q25, a gene rich locus on 7q22 and KIAA1462 on 10p11. The CAD-associated SNP in the PDGFD locus showed tissue-specific cis expression quantitative trait locus effects. These findings implicate new pathways for CAD susceptibility.     

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.     

Genome-wide strategies for detecting multiple loci that influence complex diseases

After nearly 10 years of intense academic and commercial research effort, large genome-wide association studies for common complex diseases are now imminent. Although these conditions involve a complex relationship between genotype and phenotype, including interactions between unlinked loci, the prevailing strategies for analysis of such studies focus on the locus-by-locus paradigm. Here we consider analytical methods that explicitly look for statistical interactions between loci. We show first that they are computationally feasible, even for studies of hundreds of thousands of loci, and second that even with a conservative correction for multiple testing, they can be more powerful than traditional analyses under a range of models for interlocus interactions. We also show that plausible variations across populations in allele frequencies among interacting loci can markedly affect the power to detect their marginal effects, which may account in part for the well-known difficulties in replicating association results. These results suggest that searching for interactions among genetic loci can be fruitfully incorporated into analysis strategies for genome-wide association studies.     

Two new susceptibility loci for Kawasaki disease identified through genome-wide association analysis

To find new candidate loci predisposing individuals to Kawasaki disease, an acute vasculitis that affects children, we conducted a genome-wide association study in 622 individuals with Kawasaki disease (cases) and 1,107 controls in a Han Chinese population residing in Taiwan, with replication in an independent Han Chinese sample of 261 cases and 550 controls. We report two new loci, one at BLK (encoding B-lymphoid tyrosine kinase) and one at CD40, that are associated with Kawasaki disease at genome-wide significance (P < 5 × 10(-8)). Our findings may lead to a better understanding of the role of immune activation and inflammation in Kawasaki disease pathogenesis.     

Segregation at three loci explains familial and population risk in Hirschsprung disease

Hirschsprung disease (HSCR), the most common hereditary cause of intestinal obstruction, shows considerable variation and complex inheritance. Coding sequence mutations in RET, GDNF, EDNRB, EDN3 and SOX10 lead to long-segment (L-HSCR) and syndromic HSCR but fail to explain the transmission of the much more common short-segment form (S-HSCR). We conducted a genome scan in families with S-HSCR and identified susceptibility loci at 3p21, 10q11 and 19q12 that seem to be necessary and sufficient to explain recurrence risk and population incidence. The gene at 10q11 is probably RET, supporting its crucial role in all forms of HSCR; however, coding sequence mutations are present in only 40% of linked families, suggesting the importance of noncoding variation. Here we show oligogenic inheritance of S-HSCR, the 3p21 and 19q12 loci as RET-dependent modifiers, and a parent-of-origin effect at RET. This study demonstrates by a complete genetic dissection why the inheritance pattern of S-HSCR is nonmendelian.     

Identification of two new loci at IL23R and RAB32 that influence susceptibility to leprosy

We performed a genome-wide association study with 706 individuals with leprosy and 5,581 control individuals and replicated the top 24 SNPs in three independent replication samples, including a total of 3,301 individuals with leprosy and 5,299 control individuals from China. Two loci not previously associated with the disease were identified with genome-wide significance: rs2275606 (combined P = 3.94 × 10(-14), OR = 1.30) on 6q24.3 and rs3762318 (combined P = 3.27 × 10(-11), OR = 0.69) on 1p31.3. These associations implicate IL23R and RAB32 as new susceptibility genes for leprosy. Furthermore, we identified evidence of interaction between the NOD2 and RIPK2 loci, which is consistent with the biological association of the proteins encoded by these genes (NOD2-RIPK2 complex) in activating the NF-κB pathway as a part of the host defense response to infection. Our findings have expanded the biological functions of IL23R by uncovering its involvement in infectious disease susceptibility and suggest a potential involvement of autophagocytosis in leprosy pathogenesis. The IL23R association supports previous observations of the marked overlap of susceptibility genes for leprosy and Crohn's disease, implying common pathogenesis mechanisms.     

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.     

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.     

ITPKC functional polymorphism associated with Kawasaki disease susceptibility and formation of coronary artery aneurysms

Kawasaki disease is a pediatric systemic vasculitis of unknown etiology for which a genetic influence is suspected. We identified a functional SNP (itpkc_3) in the inositol 1,4,5-trisphosphate 3-kinase C (ITPKC) gene on chromosome 19q13.2 that is significantly associated with Kawasaki disease susceptibility and also with an increased risk of coronary artery lesions in both Japanese and US children. Transfection experiments showed that the C allele of itpkc_3 reduces splicing efficiency of the ITPKC mRNA. ITPKC acts as a negative regulator of T-cell activation through the Ca2+/NFAT signaling pathway, and the C allele may contribute to immune hyper-reactivity in Kawasaki disease. This finding provides new insights into the mechanisms of immune activation in Kawasaki disease and emphasizes the importance of activated T cells in the pathogenesis of this vasculitis.     

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

Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture

Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10(-8)). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10(-4), Bonferroni corrected), of which six reached P < 5 × 10(-8), including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.     

Mapping autism risk loci using genetic linkage and chromosomal rearrangements

Autism spectrum disorders (ASDs) are common, heritable neurodevelopmental conditions. The genetic architecture of ASDs is complex, requiring large samples to overcome heterogeneity. Here we broaden coverage and sample size relative to other studies of ASDs by using Affymetrix 10K SNP arrays and 1,181 [corrected] families with at least two affected individuals, performing the largest linkage scan to date while also analyzing copy number variation in these families. Linkage and copy number variation analyses implicate chromosome 11p12-p13 and neurexins, respectively, among other candidate loci. Neurexins team with previously implicated neuroligins for glutamatergic synaptogenesis, highlighting glutamate-related genes as promising candidates for contributing to ASDs.     

Genome-wide association study identifies loci influencing concentrations of liver enzymes in plasma

Concentrations of liver enzymes in plasma are widely used as indicators of liver disease. We carried out a genome-wide association study in 61,089 individuals, identifying 42 loci associated with concentrations of liver enzymes in plasma, of which 32 are new associations (P = 10(-8) to P = 10(-190)). We used functional genomic approaches including metabonomic profiling and gene expression analyses to identify probable candidate genes at these regions. We identified 69 candidate genes, including genes involved in biliary transport (ATP8B1 and ABCB11), glucose, carbohydrate and lipid metabolism (FADS1, FADS2, GCKR, JMJD1C, HNF1A, MLXIPL, PNPLA3, PPP1R3B, SLC2A2 and TRIB1), glycoprotein biosynthesis and cell surface glycobiology (ABO, ASGR1, FUT2, GPLD1 and ST3GAL4), inflammation and immunity (CD276, CDH6, GCKR, HNF1A, HPR, ITGA1, RORA and STAT4) and glutathione metabolism (GSTT1, GSTT2 and GGT), as well as several genes of uncertain or unknown function (including ABHD12, EFHD1, EFNA1, EPHA2, MICAL3 and ZNF827). Our results provide new insight into genetic mechanisms and pathways influencing markers of liver function.     

Thirty new loci for age at menarche identified by a meta-analysis of genome-wide association studies

To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P = 5.4 × 10???) and 9q31.2 (P = 2.2 × 10?33), we identified 30 new menarche loci (all P < 5 × 10??) and found suggestive evidence for a further 10 loci (P < 1.9 × 10??). The new loci included four previously associated with body mass index (in or near FTO, SEC16B, TRA2B and TMEM18), three in or near other genes implicated in energy homeostasis (BSX, CRTC1 and MCHR2) and three in or near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and gene-set enrichment pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing.     

Genome-wide association identifies a susceptibility locus for coronary artery disease in the Chinese Han population

Coronary artery disease (CAD) causes more than 700,000 deaths each year in China. Previous genome-wide association studies (GWAS) in populations of European ancestry identified several genetic loci for CAD, but no such study has yet been reported in the Chinese population. Here we report a three-stage GWAS in the Chinese Han population. We identified a new association between rs6903956 in a putative gene denoted as C6orf105 on chromosome 6p24.1 and CAD (P = 5.00 × 10?3, stage 2 validation; P = 3.00 × 10?3, P = 1.19 × 10?? and P = 4.00 × 10?3 in three independent stage 3 replication populations; P = 4.87 × 10?12, odds ratio = 1.51 in the combined population). The minor risk allele A of rs6903956 is associated with decreased C6orf105 mRNA expression. We report the first GWAS for CAD in the Chinese Han population and identify a SNP, rs6903956, in C6orf105 associated with susceptibility to CAD in this population.