Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes

Metformin is the most commonly used pharmacological therapy for type 2 diabetes. We report a genome-wide association study for glycemic response to metformin in 1,024 Scottish individuals with type 2 diabetes with replication in two cohorts including 1,783 Scottish individuals and 1,113 individuals from the UK Prospective Diabetes Study. In a combined meta-analysis, we identified a SNP, rs11212617, associated with treatment success (n = 3,920, P = 2.9 × 10(-9), odds ratio = 1.35, 95% CI 1.22-1.49) at a locus containing ATM, the ataxia telangiectasia mutated gene. In a rat hepatoma cell line, inhibition of ATM with KU-55933 attenuated the phosphorylation and activation of AMP-activated protein kinase in response to metformin. We conclude that ATM, a gene known to be involved in DNA repair and cell cycle control, plays a role in the effect of metformin upstream of AMP-activated protein kinase, and variation in this gene alters glycemic response to metformin.     

A common haplotype of interferon regulatory factor 5 (IRF5) regulates splicing and expression and is associated with increased risk of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a complex autoimmune disease characterized by activation of the type I interferon (IFN) pathway. Here we convincingly replicate association of the IFN regulatory factor 5 (IRF5) rs2004640 T allele with SLE in four independent case-control cohorts (P = 4.4 x 10(-16)) and by family-based transmission disequilibrium test analysis (P = 0.0006). The rs2004640 T allele creates a 5' donor splice site in an alternate exon 1 of IRF5, allowing expression of several unique IRF5 isoforms. We also identify an independent cis-acting variant associated with elevated expression of IRF5 and linked to the exon 1B splice site. Haplotypes carrying the variant associated with elevated expression and lacking the exon 1B donor site do not confer risk of SLE. Thus, a common IRF5 haplotype driving elevated expression of multiple unique isoforms of IRF5 is an important genetic risk factor for SLE, establishing a causal role for type I IFN pathway genes in human autoimmunity.     

Common variants near MBNL1 and NKX2-5 are associated with infantile hypertrophic pyloric stenosis

Infantile hypertrophic pyloric stenosis (IHPS) is a severe condition characterized by hypertrophy of the pyloric sphincter muscle. We conducted a genome-wide association study (GWAS) on 1,001 surgery-confirmed cases and 2,401 controls from Denmark. The six most strongly associated loci were tested in a replication set of 796 cases and 876 controls. Three SNPs reached genome-wide significance. One of these SNPs, rs11712066 (odds ratio (OR) = 1.61; P = 1.5 × 10(-17)) at 3p25.1, is located 150 kb upstream of MBNL1, which encodes a factor that regulates splicing transitions occurring shortly after birth. The second SNP, rs573872 (OR = 1.41; P = 4.3 × 10(-12)), maps to an intergenic region at 3p25.2 approximately 1.3 Mb downstream of MBNL1. The third SNP, rs29784 (OR = 1.42; P = 1.5 × 10(-15)) at 5q35.2, is 64 kb downstream of NKX2-5, which is involved in development of cardiac muscle tissue and embryonic gut development.     

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.     

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 common polymorphism of the growth hormone receptor is associated with increased responsiveness to growth hormone

Growth hormone is used to increase height in short children who are not deficient in growth hormone, but its efficacy varies largely across individuals. The genetic factors responsible for this variation are entirely unknown. In two cohorts of short children treated with growth hormone, we found that an isoform of the growth hormone receptor gene that lacks exon 3 (d3-GHR) was associated with 1.7 to 2 times more growth acceleration induced by growth hormone than the full-length isoform (P < 0.0001). In transfection experiments, the transduction of growth hormone signaling through d3-GHR homo- or heterodimers was approximately 30% higher than through full-length GHR homodimers (P < 0.0001). One-half of Europeans are hetero- or homozygous with respect to the allele encoding the d3-GHR isoform, which is dominant over the full-length isoform. These observations suggest that the polymorphism in exon 3 of GHR is important in growth hormone pharmacogenetics.     

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

Multiple apical plasma membrane constituents are associated with susceptibility to meconium ileus in individuals with cystic fibrosis

Variants associated with meconium ileus in cystic fibrosis were identified in 3,763 affected individuals by genome-wide association study (GWAS). Five SNPs at two loci near SLC6A14 at Xq23-24 (minimum P = 1.28 × 10(-12) at rs3788766) and SLC26A9 at 1q32.1 (minimum P = 9.88 × 10(-9) at rs4077468) accounted for ~5% of phenotypic variability and were replicated in an independent sample of affected individuals (n = 2,372; P = 0.001 and 0.0001, respectively). By incorporating the knowledge that disease-causing mutations in CFTR alter electrolyte and fluid flux across surface epithelium into a hypothesis-driven GWAS (GWAS-HD), we identified associations with the same SNPs in SLC6A14 and SLC26A9 and established evidence for the involvement of SNPs in a third solute carrier gene, SLC9A3. In addition, GWAS-HD provided evidence of association between meconium ileus and multiple genes encoding constituents of the apical plasma membrane where CFTR resides (P = 0.0002; testing of 155 apical membrane genes jointly and in replication, P = 0.022). These findings suggest that modulating activities of apical membrane constituents could complement current therapeutic paradigms for cystic fibrosis.     

Common variants in DGKK are strongly associated with risk of hypospadias

Hypospadias is a common congenital malformation of the male external genitalia. We performed a genome-wide association study using pooled DNA from 436 individuals with hypospadias (cases) and 494 controls of European descent and selected the highest ranked SNPs for individual genotyping in the discovery sample, an additional Dutch sample of 133 cases and their parents, and a Swedish series of 266 cases and 402 controls. Individual genotyping of two SNPs (rs1934179 and rs7063116) in DGKK, encoding diacylglycerol kinase κ, produced compelling evidence for association with hypospadias in the discovery sample (allele-specific odds ratio (OR) = 2.5, P = 2.5 × 10?11 and OR = 2.3, P = 2.9 × 10??, respectively) and in the Dutch (OR = 3.9, P = 2.4 × 10?? and OR = 3.8, P = 3.4 × 10??) and Swedish (OR = 2.5, P = 2.6 × 10?? and OR = 2.2, P = 2.7 × 10??) replication samples. Expression studies showed expression of DGKK in preputial tissue of cases and controls, which was lower in carriers of the risk allele of rs1934179 (P = 0.047). We propose DGKK as a major risk gene for hypospadias.     

Mutations in the gene encoding the 3'-5' DNA exonuclease TREX1 are associated with systemic lupus erythematosus

TREX1 acts in concert with the SET complex in granzyme A-mediated apoptosis, and mutations in TREX1 cause Aicardi-Goutières syndrome and familial chilblain lupus. Here, we report monoallelic frameshift or missense mutations and one 3' UTR variant of TREX1 present in 9/417 individuals with systemic lupus erythematosus but absent in 1,712 controls (P = 4.1 x 10(-7)). We demonstrate that two mutant TREX1 alleles alter subcellular targeting. Our findings implicate TREX1 in the pathogenesis of SLE.     

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.     

Common variants in DVWA on chromosome 3p24.3 are associated with susceptibility to knee osteoarthritis

Susceptibility to osteoarthritis, the most common human arthritis, is known to be influenced by genetic factors. Through a genome-wide association study using approximately 100,000 SNPs, we have identified a previously unknown gene on chromosome 3p24.3, DVWA, which is associated with susceptibility to knee osteoarthritis. Expressed specifically in cartilage, DVWA encodes a 276-amino-acid protein with two regions corresponding to the von Willebrand factor type A domain (VWA domain). Several DVWA SNPs are significantly associated with knee osteoarthritis in two independent Japanese case-control cohorts. This association was replicated in a Japanese population cohort and a Han Chinese case-control cohort (combined P = 7.3 x 10(-11)). DVWA protein binds to beta-tubulin, and the binding is influenced by two highly associated missense SNPs (rs11718863 and rs7639618) located in the VWA domain. The Tyr169-Cys260 isoform of DVWA, which is overrepresented in knee osteoarthritis, showed weaker interaction. Our findings reveal a new paradigm for study of osteoarthritis etiology and pathogenesis.