Chris Wallace

Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes

The Wellcome Trust Case Control Consortium (WTCCC) primary genome-wide association (GWA) scan on seven diseases, including the multifactorial autoimmune disease type 1 diabetes (T1D), shows associations at P < 5 × 10−7 between T1D and six chromosome regions: 12q24, 12q13, 16p13, 18p11, 12p13 and 4q27. Here, we attempted to validate these and six other top findings in 4,000 individuals with T1D, 5,000 controls and 2,997 family trios independent of the WTCCC study. We confirmed unequivocally the associations of 12q24, 12q13, 16p13 and 18p11 (Pfollow-up ≤ 1.35 × 10−9; Poverall ≤ 1.15 × 10−14), leaving eight regions with small effects or false-positive associations. We also obtained evidence for chromosome 18q22 (Poverall = 1.38 × 10−8) from a GWA study of nonsynonymous SNPs. Several regions, including 18q22 and 18p11, showed association with autoimmune thyroid disease. This study increases the number of T1D loci with compelling evidence from six to at least ten.

Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease.

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.

Meta-Analysis of 28,141 Individuals Identifies Common Variants within Five New Loci That Influence Uric Acid Concentrations

Elevated serum uric acid levels cause gout and are a risk factor for cardiovascular disease and diabetes. To investigate the polygenetic basis of serum uric acid levels, we conducted a meta-analysis of genome-wide association scans from 14 studies totalling 28,141 participants of European descent, resulting in identification of 954 SNPs distributed across nine loci that exceeded the threshold of genome-wide significance, five of which are novel. Overall, the common variants associated with serum uric acid levels fall in the following nine regions: SLC2A9 (p = 5.2×10−201), ABCG2 (p = 3.1×10−26), SLC17A1 (p = 3.0×10−14), SLC22A11 (p = 6.7×10−14), SLC22A12 (p = 2.0×10−9), SLC16A9 (p = 1.1×10−8), GCKR (p = 1.4×10−9), LRRC16A (p = 8.5×10−9), and near PDZK1 (p = 2.7×10−9). Identified variants were analyzed for gender differences. We found that the minor allele for rs734553 in SLC2A9 has greater influence in lowering uric acid levels in women and the minor allele of rs2231142 in ABCG2 elevates uric acid levels more strongly in men compared to women. To further characterize the identified variants, we analyzed their association with a panel of metabolites. rs12356193 within SLC16A9 was associated with DL-carnitine (p = 4.0×10−26) and propionyl-L-carnitine (p = 5.0×10−8) concentrations, which in turn were associated with serum UA levels (p = 1.4×10−57 and p = 8.1×10−54, respectively), forming a triangle between SNP, metabolites, and UA levels. Taken together, these associations highlight additional pathways that are important in the regulation of serum uric acid levels and point toward novel potential targets for pharmacological intervention to prevent or treat hyperuricemia. In addition, these findings strongly support the hypothesis that transport proteins are key in regulating serum uric acid levels.

Genome-wide association study identifies genes for biomarkers of cardiovascular disease: serum urate and dyslipidemia.

Many common diseases are accompanied by disturbances in biochemical traits. Identifying the genetic determinants could provide novel insights into disease mechanisms and reveal avenues for developing new therapies. Here, we report a genome-wide association analysis for commonly measured serum and urine biochemical traits. As part of the WTCCC, 500,000 SNPs genome wide were genotyped in 1955 hypertensive individuals characterized for 25 serum and urine biochemical traits. For each trait, we assessed association with individual SNPs, adjusting for age, sex, and BMI. Lipid measurements were further examined in a meta-analysis of genome-wide data from a type 2 diabetes scan. The most promising associations were examined in two epidemiological cohorts. We discovered association between serum urate and SLC2A9, a glucose transporter (p = 2 x 10(-15)) and confirmed this in two independent cohorts, GRAPHIC study (p = 9 x 10(-15)) and TwinsUK (p = 8 x 10(-19)). The odds ratio for hyperuricaemia (defined as urate >0.4 mMol/l) is 1.89 (95% CI = 1.36-2.61) per copy of common allele. We also replicated many genes previously associated with serum lipids and found previously recognized association between LDL levels and SNPs close to genes encoding PSRC1 and CELSR2 (p = 1 x 10(-7)). The common allele was associated with a 6% increase in nonfasting serum LDL. This region showed increased association in the meta-analysis (p = 4 x 10(-14)). This finding provides a potential biological mechanism for the recent association of this same allele of the same SNP with increased risk of coronary disease.

Pervasive sharing of genetic effects in autoimmune disease.

Genome-wide association (GWA) studies have identified numerous, replicable, genetic associations between common single nucleotide polymorphisms (SNPs) and risk of common autoimmune and inflammatory (immune-mediated) diseases, some of which are shared between two diseases. Along with epidemiological and clinical evidence, this suggests that some genetic risk factors may be shared across diseases—as is the case with alleles in the Major Histocompatibility Locus. In this work we evaluate the extent of this sharing for 107 immune disease-risk SNPs in seven diseases: celiac disease, Crohns disease, multiple sclerosis, psoriasis, rheumatoid arthritis, systemic lupus erythematosus, and type 1 diabetes. We have developed a novel statistic for Cross Phenotype Meta-Analysis (CPMA) which detects association of a SNP to multiple, but not necessarily all, phenotypes. With it, we find evidence that 47/107 (44%) immune-mediated disease risk SNPs are associated to multiple—but not all—immune-mediated diseases (SNP-wise P CPMA<0.01). We also show that distinct groups of interacting proteins are encoded near SNPs which predispose to the same subsets of diseases; we propose these as the mechanistic basis of shared disease risk. We are thus able to leverage genetic data across diseases to construct biological hypotheses about the underlying mechanism of pathogenesis.

SLC2A9 Is a High-Capacity Urate Transporter in Humans

Background Serum uric acid levels in humans are influenced by diet, cellular breakdown, and renal elimination, and correlate with blood pressure, metabolic syndrome, diabetes, gout, and cardiovascular disease. Recent genome-wide association scans have found common genetic variants of SLC2A9 to be associated with increased serum urate level and gout. The SLC2A9 gene encodes a facilitative glucose transporter, and it has two splice variants that are highly expressed in the proximal nephron, a key site for urate handling in the kidney. We investigated whether SLC2A9 is a functional urate transporter that contributes to the longstanding association between urate and blood pressure in man. Methods and Findings We expressed both SLC2A9 splice variants in Xenopus laevis oocytes and found both isoforms mediate rapid urate fluxes at concentration ranges similar to physiological serum levels (200–500 μM). Because SLC2A9 is a known facilitative glucose transporter, we also tested whether glucose or fructose influenced urate transport. We found that urate is transported by SLC2A9 at rates 45- to 60-fold faster than glucose, and demonstrated that SLC2A9-mediated urate transport is facilitated by glucose and, to a lesser extent, fructose. In addition, transport is inhibited by the uricosuric benzbromarone in a dose-dependent manner (K i = 27 μM). Furthermore, we found urate uptake was at least 2-fold greater in human embryonic kidney (HEK) cells overexpressing SLC2A9 splice variants than nontransfected kidney cells. To confirm that our findings were due to SLC2A9, and not another urate transporter, we showed that urate transport was diminished by SLC2A9-targeted siRNA in a second mammalian cell line. In a cohort of men we showed that genetic variants of SLC2A9 are associated with reduced urinary urate clearance, which fits with common variation at SLC2A9 leading to increased serum urate. We found no evidence of association with hypertension (odds ratio 0.98, 95% confidence interval [CI] 0.9 to 1.05, p > 0.33) by meta-analysis of an SLC2A9 variant in six case–control studies including 11,897 participants. In a separate meta-analysis of four population studies including 11,629 participants we found no association of SLC2A9 with systolic (effect size −0.12 mm Hg, 95% CI −0.68 to 0.43, p = 0.664) or diastolic blood pressure (effect size −0.03 mm Hg, 95% CI −0.39 to 0.31, p = 0.82). Conclusions This study provides evidence that SLC2A9 splice variants act as high-capacity urate transporters and is one of the first functional characterisations of findings from genome-wide association scans. We did not find an association of the SLC2A9 gene with blood pressure in this study. Our findings suggest potential pathogenic mechanisms that could offer a new drug target for gout.

A trans-acting locus regulates an anti-viral expression network and type 1 diabetes risk

Combined analyses of gene networks and DNA sequence variation can provide new insights into the aetiology of common diseases that may not be apparent from genome-wide association studies alone. Recent advances in rat genomics are facilitating systems-genetics approaches. Here we report the use of integrated genome-wide approaches across seven rat tissues to identify gene networks and the loci underlying their regulation. We defined an interferon regulatory factor 7 (IRF7)-driven inflammatory network (IDIN) enriched for viral response genes, which represents a molecular biomarker for macrophages and which was regulated in multiple tissues by a locus on rat chromosome 15q25. We show that Epstein–Barr virus induced gene 2 (Ebi2, also known as Gpr183), which lies at this locus and controls B lymphocyte migration, is expressed in macrophages and regulates the IDIN. The human orthologous locus on chromosome 13q32 controlled the human equivalent of the IDIN, which was conserved in monocytes. IDIN genes were more likely to associate with susceptibility to type 1 diabetes (T1D)—a macrophage-associated autoimmune disease—than randomly selected immune response genes (P = 8.85 × 10−6). The human locus controlling the IDIN was associated with the risk of T1D at single nucleotide polymorphism rs9585056 (P = 7.0 × 10−10; odds ratio, 1.15), which was one of five single nucleotide polymorphisms in this region associated with EBI2 (GPR183) expression. These data implicate IRF7 network genes and their regulatory locus in the pathogenesis of T1D.

Inherited Variation in Vitamin D Genes Is Associated With Predisposition to Autoimmune Disease Type 1 Diabetes

OBJECTIVE Vitamin D deficiency (25-hydroxyvitamin D [25(OH)D] <50 nmol/L) is commonly reported in both children and adults worldwide, and growing evidence indicates that vitamin D deficiency is associated with many extraskeletal chronic disorders, including the autoimmune diseases type 1 diabetes and multiple sclerosis. RESEARCH DESIGN AND METHODS We measured 25(OH)D concentrations in 720 case and 2,610 control plasma samples and genotyped single nucleotide polymorphisms from seven vitamin D metabolism genes in 8,517 case, 10,438 control, and 1,933 family samples. We tested genetic variants influencing 25(OH)D metabolism for an association with both circulating 25(OH)D concentrations and disease status. RESULTS Type 1 diabetic patients have lower circulating levels of 25(OH)D than similarly aged subjects from the British population. Only 4.3 and 18.6% of type 1 diabetic patients reached optimal levels (≥75 nmol/L) of 25(OH)D for bone health in the winter and summer, respectively. We replicated the associations of four vitamin D metabolism genes (GC, DHCR7, CYP2R1, and CYP24A1) with 25(OH)D in control subjects. In addition to the previously reported association between type 1 diabetes and CYP27B1 (P = 1.4 × 10−4), we obtained consistent evidence of type 1 diabetes being associated with DHCR7 (P = 1.2 × 10−3) and CYP2R1 (P = 3.0 × 10−3). CONCLUSIONS Circulating levels of 25(OH)D in children and adolescents with type 1 diabetes vary seasonally and are under the same genetic control as in the general population but are much lower. Three key 25(OH)D metabolism genes show consistent evidence of association with type 1 diabetes risk, indicating a genetic etiological role for vitamin D deficiency in type 1 diabetes.

The imprinted DLK1-MEG3 gene region on chromosome 14q32.2 alters susceptibility to type 1 diabetes

Genome-wide association (GWA) studies to map common disease susceptibility loci have been hugely successful, with over 300 reproducibly associated loci reported to date. However, these studies have not yet provided convincing evidence for any susceptibility locus subject to parent-of-origin effects. Using imputation to extend existing GWA datasets, we have obtained robust evidence at rs941576 for paternally inherited risk of type 1 diabetes (T1D; ratio of allelic effects for paternal versus maternal transmissions = 0.75; 95% confidence interval (CI) = 0.71–0.79). This marker is in the imprinted region of chromosome 14q32.2, which contains the functional candidate gene DLK1. Our meta-analysis also provided support at genome-wide significance for a T1D locus at chromosome 19p13.2. The highest association was at marker rs2304256 (odds ratio (OR) = 0.86; 95%CI = 0.82–0.90) in the TYK2 gene, which has previously been associated with systemic lupus erythematosus and multiple sclerosis.

Bayesian test for colocalisation between pairs of genetic association studies using summary statistics.

Genetic association studies, in particular the genome-wide association study (GWAS) design, have provided a wealth of novel insights into the aetiology of a wide range of human diseases and traits, in particular cardiovascular diseases and lipid biomarkers. The next challenge consists of understanding the molecular basis of these associations. The integration of multiple association datasets, including gene expression datasets, can contribute to this goal. We have developed a novel statistical methodology to assess whether two association signals are consistent with a shared causal variant. An application is the integration of disease scans with expression quantitative trait locus (eQTL) studies, but any pair of GWAS datasets can be integrated in this framework. We demonstrate the value of the approach by re-analysing a gene expression dataset in 966 liver samples with a published meta-analysis of lipid traits including >100,000 individuals of European ancestry. Combining all lipid biomarkers, our re-analysis supported 26 out of 38 reported colocalisation results with eQTLs and identified 14 new colocalisation results, hence highlighting the value of a formal statistical test. In three cases of reported eQTL-lipid pairs (SYPL2, IFT172, TBKBP1) for which our analysis suggests that the eQTL pattern is not consistent with the lipid association, we identify alternative colocalisation results with SORT1, GCKR, and KPNB1, indicating that these genes are more likely to be causal in these genomic intervals. A key feature of the method is the ability to derive the output statistics from single SNP summary statistics, hence making it possible to perform systematic meta-analysis type comparisons across multiple GWAS datasets (implemented online at http://coloc.cs.ucl.ac.uk/coloc/). Our methodology provides information about candidate causal genes in associated intervals and has direct implications for the understanding of complex diseases as well as the design of drugs to target disease pathways.