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

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.     

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.     

A functional SNP in CILP, encoding cartilage intermediate layer protein, is associated with susceptibility to lumbar disc disease

Lumbar disc disease (LDD) is caused by degeneration of intervertebral discs of the lumbar spine. One of the most common musculoskeletal disorders, LDD has strong genetic determinants. Using a case-control association study, we identified a functional SNP (1184T --> C, resulting in the amino acid substitution I395T) in CILP, which encodes the cartilage intermediate layer protein, that acts as a modulator of LDD susceptibility. CILP was expressed abundantly in intervertebral discs, and its expression increased as disc degeneration progressed. CILP colocalized with TGF-beta1 in clustering chondrocytes and their territorial matrices in intervertebral discs. CILP inhibited TGF-beta1-mediated induction of cartilage matrix genes through direct interaction with TGF-beta1 and inhibition of TGF-beta1 signaling. The susceptibility-associated 1184C allele showed increased binding and inhibition of TGF-beta1. Therefore, we conclude that the extracellular matrix protein CILP regulates TGF-beta signaling and that this regulation has a crucial role in the etiology and pathogenesis of LDD. Our study also adds to the list of connective tissue diseases that are associated with TGF-beta.     

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.     

A Mal functional variant is associated with protection against invasive pneumococcal disease, bacteremia, malaria and tuberculosis

Toll-like receptors (TLRs) and members of their signaling pathway are important in the initiation of the innate immune response to a wide variety of pathogens. The adaptor protein Mal (also known as TIRAP), encoded by TIRAP (MIM 606252), mediates downstream signaling of TLR2 and TLR4 (refs. 4-6). We report a case-control study of 6,106 individuals from the UK, Vietnam and several African countries with invasive pneumococcal disease, bacteremia, malaria and tuberculosis. We genotyped 33 SNPs, including rs8177374, which encodes a leucine substitution at Ser180 of Mal. We found that heterozygous carriage of this variant associated independently with all four infectious diseases in the different study populations. Combining the study groups, we found substantial support for a protective effect of S180L heterozygosity against these infectious diseases (N = 6,106; overall P = 9.6 x 10(-8)). We found that the Mal S180L variant attenuated TLR2 signal transduction.     

Variant between CPT1B and CHKB associated with susceptibility to narcolepsy

Narcolepsy (hypocretin deficiency), a sleep disorder characterized by sleepiness, cataplexy and rapid eye movement (REM) sleep abnormalities, is tightly associated with HLA-DRB1*1501 (M17378) and HLA-DQB1*0602 (M20432). Susceptibility genes other than those in the HLA region are also likely involved. We conducted a genome-wide association study using 500K SNP microarrays in 222 Japanese individuals with narcolepsy and 389 Japanese controls, with replication of top hits in 159 Japanese individuals with narcolepsy and 190 Japanese controls, followed by the testing of 424 Koreans, 785 individuals of European descent and 184 African Americans. rs5770917, a SNP located between CPT1B and CHKB, was associated with narcolepsy in Japanese (rs5770917[C], odds ratio (OR) = 1.79, combined P = 4.4 x 10(-7)) and other ancestry groups (OR = 1.40, P = 0.02). Real-time quantitative PCR assays in white blood cells indicated decreased CPT1B and CHKB expression in subjects with the C allele, suggesting that a genetic variant regulating CPT1B or CHKB expression is associated with narcolepsy. Either of these genes is a plausible candidate, as CPT1B regulates beta-oxidation, a pathway involved in regulating theta frequency during REM sleep, and CHKB is an enzyme involved in the metabolism of choline, a precursor of the REM- and wake-regulating neurotransmitter acetylcholine.     

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.     

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.     

Common variants at 11p13 are associated with susceptibility to tuberculosis

After imputation of data from the 1000 Genomes Project into a genome-wide dataset of Ghanaian individuals with tuberculosis and controls, we identified a resistance locus on chromosome 11p13 downstream of the WT1 gene (encoding Wilms tumor 1). The strongest signal was obtained at the rs2057178 SNP (P = 2.63 × 10(-9)). Replication in Gambian, Indonesian and Russian tuberculosis case-control study cohorts increased the significance level for the association with this SNP to P = 2.57 × 10(-11).     

TGFB2 mutations cause familial thoracic aortic aneurysms and dissections associated with mild systemic features of Marfan syndrome

A predisposition for thoracic aortic aneurysms leading to acute aortic dissections can be inherited in families in an autosomal dominant manner. Genome-wide linkage analysis of two large unrelated families with thoracic aortic disease followed by whole-exome sequencing of affected relatives identified causative mutations in TGFB2. These mutations-a frameshift mutation in exon 6 and a nonsense mutation in exon 4-segregated with disease with a combined logarithm of odds (LOD) score of 7.7. Sanger sequencing of 276 probands from families with inherited thoracic aortic disease identified 2 additional TGFB2 mutations. TGFB2 encodes transforming growth factor (TGF)-β2, and the mutations are predicted to cause haploinsufficiency for TGFB2; however, aortic tissue from cases paradoxically shows increased TGF-β2 expression and immunostaining. Thus, haploinsufficiency for TGFB2 predisposes to thoracic aortic disease, suggesting that the initial pathway driving disease is decreased cellular TGF-β2 levels leading to a secondary increase in TGF-β2 production in the diseased aorta.     

An aspartic acid repeat polymorphism in asporin inhibits chondrogenesis and increases susceptibility to osteoarthritis

Osteoarthritis is the most common form of human arthritis. We investigated the potential role of asporin, an extracellular matrix component expressed abundantly in the articular cartilage of individuals with osteoarthritis, in the pathogenesis of osteoarthritis. Here we report a significant association between a polymorphism in the aspartic acid (D) repeat of the gene encoding asporin (ASPN) and osteoarthritis. In two independent populations of individuals with knee osteoarthritis, the D14 allele of ASPN is over-represented relative to the common D13 allele, and its frequency increases with disease severity. The D14 allele is also over-represented in individuals with hip osteoarthritis. Asporin suppresses TGF-beta-mediated expression of the genes aggrecan (AGC1) and type II collagen (COL2A1) and reduced proteoglycan accumulation in an in vitro model of chondrogenesis. The effect on TGF-beta activity is allele-specific, with the D14 allele resulting in greater inhibition than other alleles. In vitro binding assays showed a direct interaction between asporin and TGF-beta. Taken together, these findings provide another functional link between extracellular matrix proteins, TGF-beta activity and disease, suggesting new therapeutic strategies for osteoarthritis.     

Common variants at TRAF3IP2 are associated with susceptibility to psoriatic arthritis and psoriasis

Psoriatic arthritis (PsA) is an inflammatory joint disease that is distinct from other chronic arthritides and which is frequently accompanied by psoriasis vulgaris (PsV) and seronegativity for rheumatoid factor. We conducted a genome-wide association study in 609 German individuals with PsA (cases) and 990 controls with replication in 6 European cohorts including a total of 5,488 individuals. We replicated PsA associations at HLA-C and IL12B and identified a new association at TRAF3IP2 (rs13190932, P = 8.56 × 10?1?). TRAF3IP2 was also associated with PsV in a German cohort including 2,040 individuals (rs13190932, P = 1.95 × 10?3). Sequencing of the exons of TRAF3IP2 identified a coding variant (p.Asp10Asn, rs33980500) as the most significantly associated SNP (P = 1.13 × 10?2?, odds ratio = 1.95). Functional assays showed reduced binding of this TRAF3IP2 variant to TRAF6, suggesting altered modulation of immunoregulatory signals through altered TRAF interactions as a new and shared pathway for PsA and PsV.