Genetic variants regulating immune cell levels in health and disease.

The complex network of specialized cells and molecules in the immune system has evolved to defend against pathogens, but inadvertent immune system attacks on "self" result in autoimmune disease. Both genetic regulation of immune cell levels and their relationships with autoimmunity are largely undetermined. Here, we report genetic contributions to quantitative levels of 95 cell types encompassing 272 immune traits, in a cohort of 1,629 individuals from four clustered Sardinian villages. We first estimated trait heritability, showing that it can be substantial, accounting for up to 87% of the variance (mean 41%). Next, by assessing ∼8.2 million variants that we identified and confirmed in an extended set of 2,870 individuals, 23 independent variants at 13 loci associated with at least one trait. Notably, variants at three loci (HLA, IL2RA, and SH2B3/ATXN2) overlap with known autoimmune disease associations. These results connect specific cellular phenotypes to specific genetic variants, helping to explicate their involvement in disease.

GWAS of genetic factors affecting white blood cell morphological parameters in Sardinians uncovers influence of chromosome 11 innate immunity gene cluster on eosinophil morphology.

Few genome-wide association studies (GWAS) analyzing genetic regulation of morphological traits of white blood cells have been reported. We carried out a GWAS of 12 morphological traits in 869 individuals from the general population of Sardinia, Italy. These traits, included measures of cell volume, conductivity and light scatter in four white-cell populations (eosinophils, lymphocytes, monocytes, neutrophils). This analysis yielded seven statistically significant signals, four of which were novel (four novel, PRG2, P2RX3, two of CDK6). Five signals were replicated in the independent INTERVAL cohort of 11 822 individuals. The most interesting signal with large effect size on eosinophil scatter (P-value = 8.33 x 10-32, beta = -1.651, se = 0.1351) falls within the innate immunity cluster on chromosome 11, and is located in the PRG2 gene. Computational analyses revealed that a rare, Sardinian-specific PRG2:p.Ser148Pro mutation modifies PRG2 amino acid contacts and protein dynamics in a manner that could possibly explain the changes observed in eosinophil morphology. Our discoveries shed light on genetics of morphological traits. For the first time, we describe such large effect size on eosinophils morphology that is relatively frequent in Sardinian population.

PRF1 mutation alters immune system activation, inflammation, and risk of autoimmunity.

BACKGROUND: Defective alleles within the PRF1 gene, encoding the pore-forming protein perforin, in combination with environmental factors, cause familial type 2 hemophagocytic lymphohistiocytosis (FHL2), a rare, severe autosomal recessive childhood disorder characterized by massive release of cytokines-cytokine storm. OBJECTIVE: The aim of this study was to determine the function of hypomorph PRF1:p.A91V g.72360387 G > A on multiple sclerosis (MS) and type 1 diabetes (T1D). METHODS: We cross-compare the association data for PRF1:p.A91V mutation derived from GWAS on adult MS and pediatric T1D in Sardinians. The novel association with T1D was replicated in metanalysis in 12,584 cases and 17,692 controls from Sardinia, the United Kingdom, and Scotland. To dissect this mutation function, we searched through the coincident association immunophenotypes in additional set of general population Sardinians. RESULTS: We report that PRF1:p.A91V, is associated with increase of lymphocyte levels, especially within the cytotoxic memory T-cells, at general population level with reduced interleukin 7 receptor expression on these cells. The minor allele increased risk of MS, in 2903 cases and 2880 controls from Sardinia p = 2.06 × 10-4, odds ratio OR = 1.29, replicating a previous finding, whereas it protects from T1D p = 1.04 × 10-5, OR = 0.82. CONCLUSION: Our results indicate opposing contributions of the cytotoxic T-cell compartment to MS and T1D pathogenesis.

Whole-Genome Sequencing Coupled to Imputation Discovers Genetic Signals for Anthropometric Traits.

Deep sequence-based imputation can enhance the discovery power of genome-wide association studies by assessing previously unexplored variation across the common- and low-frequency spectra. We applied a hybrid whole-genome sequencing (WGS) and deep imputation approach to examine the broader allelic architecture of 12 anthropometric traits associated with height, body mass, and fat distribution in up to 267,616 individuals. We report 106 genome-wide significant signals that have not been previously identified, including 9 low-frequency variants pointing to functional candidates. Of the 106 signals, 6 are in genomic regions that have not been implicated with related traits before, 28 are independent signals at previously reported regions, and 72 represent previously reported signals for a different anthropometric trait. 71% of signals reside within genes and fine mapping resolves 23 signals to one or two likely causal variants. We confirm genetic overlap between human monogenic and polygenic anthropometric traits and find signal enrichment in cis expression QTLs in relevant tissues. Our results highlight the potential of WGS strategies to enhance biologically relevant discoveries across the frequency spectrum.

Overexpression of the Cytokine BAFF and Autoimmunity Risk.

BACKGROUND: Genomewide association studies of autoimmune diseases have mapped hundreds of susceptibility regions in the genome. However, only for a few association signals has the causal gene been identified, and for even fewer have the causal variant and underlying mechanism been defined. Coincident associations of DNA variants affecting both the risk of autoimmune disease and quantitative immune variables provide an informative route to explore disease mechanisms and drug-targetable pathways. METHODS: Using case-control samples from Sardinia, Italy, we performed a genomewide association study in multiple sclerosis followed by TNFSF13B locus-specific association testing in systemic lupus erythematosus (SLE). Extensive phenotyping of quantitative immune variables, sequence-based fine mapping, cross-population and cross-phenotype analyses, and gene-expression studies were used to identify the causal variant and elucidate its mechanism of action. Signatures of positive selection were also investigated. RESULTS: A variant in TNFSF13B, encoding the cytokine and drug target B-cell activating factor (BAFF), was associated with multiple sclerosis as well as SLE. The disease-risk allele was also associated with up-regulated humoral immunity through increased levels of soluble BAFF, B lymphocytes, and immunoglobulins. The causal variant was identified: an insertion-deletion variant, GCTGT→A (in which A is the risk allele), yielded a shorter transcript that escaped microRNA inhibition and increased production of soluble BAFF, which in turn up-regulated humoral immunity. Population genetic signatures indicated that this autoimmunity variant has been evolutionarily advantageous, most likely by augmenting resistance to malaria. CONCLUSIONS: A TNFSF13B variant was associated with multiple sclerosis and SLE, and its effects were clarified at the population, cellular, and molecular levels. (Funded by the Italian Foundation for Multiple Sclerosis and others.).

Cooperative translational control of polymorphic BAFF by NF90 and miR-15a.

Polymorphisms in untranslated regions (UTRs) of disease-associated mRNAs can alter protein production. We recently identified a genetic variant in the 3'UTR of the TNFSF13B gene, encoding the cytokine BAFF (B-cell-activating factor), that generates an alternative polyadenylation site yielding a shorter, more actively translated variant, BAFF-var mRNA. Accordingly, individuals bearing the TNFSF13B variant had higher circulating BAFF and elevated risk of developing autoimmune diseases. Here, we investigated the molecular mechanisms controlling the enhanced translation of BAFF-var mRNA. We identified nuclear factor 90 (NF90, also known as ILF3) as an RNA-binding protein that bound preferentially the wild-type (BAFF-WT mRNA) but not BAFF-var mRNA in human monocytic leukemia THP-1 cells. NF90 selectively suppressed BAFF translation by recruiting miR-15a to the 3'UTR of BAFF-WT mRNA. Our results uncover a paradigm whereby an autoimmunity-causing BAFF polymorphism prevents NF90-mediated recruitment of microRNAs to suppress BAFF translation, raising the levels of disease-associated BAFF.

Methods for association analysis and meta-analysis of rare variants in families.

Advances in exome sequencing and the development of exome genotyping arrays are enabling explorations of association between rare coding variants and complex traits. To ensure power for these rare variant analyses, a variety of association tests that group variants by gene or functional unit have been proposed. Here, we extend these tests to family-based studies. We develop family-based burden tests, variable frequency threshold tests and sequence kernel association tests. Through simulations, we compare the performance of different tests. We describe situations where family-based studies provide greater power than studies of unrelated individuals to detect rare variants associated with moderate to large changes in trait values. Broadly speaking, we find that when sample sizes are limited and only a modest fraction of all trait-associated variants can be identified, family samples are more powerful. Finally, we illustrate our approach by analyzing the relationship between coding variants and levels of high-density lipoprotein (HDL) cholesterol in 11,556 individuals from the HUNT and SardiNIA studies, demonstrating association for coding variants in the APOC3, CETP, LIPC, LIPG, and LPL genes and illustrating the value of family samples, meta-analysis, and gene-level tests. Our methods are implemented in freely available C++ code.

Genome-wide association study of susceptibility loci for breast cancer in Sardinian population.

BACKGROUND: Despite progress in identifying genes associated with breast cancer, many more risk loci exist. Genome-wide association analyses in genetically-homogeneous populations, such as that of Sardinia (Italy), could represent an additional approach to detect low penetrance alleles. METHODS: We performed a genome-wide association study comparing 1431 Sardinian patients with non-familial, BRCA1/2-mutation-negative breast cancer to 2171 healthy Sardinian blood donors. DNA was genotyped using GeneChip Human Mapping 500 K Arrays or Genome-Wide Human SNP Arrays 6.0. To increase genomic coverage, genotypes of additional SNPs were imputed using data from HapMap Phase II. After quality control filtering of genotype data, 1367 cases (9 men) and 1658 controls (1156 men) were analyzed on a total of 2,067,645 SNPs. RESULTS: Overall, 33 genomic regions (67 candidate SNPs) were associated with breast cancer risk at the p < 0(-6) level. Twenty of these regions contained defined genes, including one already associated with breast cancer risk: TOX3. With a lower threshold for preliminary significance to p < 10(-5), we identified 11 additional SNPs in FGFR2, a well-established breast cancer-associated gene. Ten candidate SNPs were selected, excluding those already associated with breast cancer, for technical validation as well as replication in 1668 samples from the same population. Only SNP rs345299, located in intron 1 of VAV3, remained suggestively associated (p-value, 1.16 x 10(-5)), but it did not associate with breast cancer risk in pooled data from two large, mixed-population cohorts. CONCLUSIONS: This study indicated the role of TOX3 and FGFR2 as breast cancer susceptibility genes in BRCA1/2-wild-type breast cancer patients from Sardinian population.

The burden of multiple sclerosis variants in continental Italians and Sardinians.

BACKGROUND: Recent studies identified > 100 non-HLA (human leukocyte antigen) multiple sclerosis (MS) susceptibility variants in Northern European populations, but their role in Southern Europeans is largely unexplored. OBJECTIVE: We aimed to investigate the cumulative impact of those variants in two Mediterranean populations: Continental Italians and Sardinians. METHODS: We calculated four weighted Genetic Risk Scores (wGRS), using up to 102 non-HLA MS risk variants and 5 HLA MS susceptibility markers in 1691 patients and 2194 controls from continental Italy; and 2861 patients and 3034 controls from Sardinia. We then assessed the differences between populations using Nagelkerke's R(2) and the area under the Receiver Operating Characteristic (ROC) curves. RESULTS: As expected, the genetic burden (mean wGRS value) was significantly higher in MS patients than in controls, in both populations. Of note, the burden was significantly higher in Sardinians. Conversely, the proportion of variability explained and the predictive power were significantly higher in continental Italians. Notably, within the Sardinian patients, we also observed a significantly higher burden of non-HLA variants in individuals who do not carry HLA risk alleles. CONCLUSIONS: The observed differences in MS genetic burden between the two Mediterranean populations highlight the need for more genetic studies in South Europeans, to further expand the knowledge of MS genetics.

A genome-wide association scan on the levels of markers of inflammation in Sardinians reveals associations that underpin its complex regulation.

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals-5 of which were identified only with the custom arrays-and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10(-29)); for ESR, at the HBB (rs4910472, p = 2.31×10(-11)) and UCN119B/SPPL3 (rs11829037, p = 8.91×10(-10)) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10(-13)) and in CADM3 (rs3026968, p = 7.63×10(-13)); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10(-21)), near DARC (rs3845624, p = 1.43×10(-10)), UNC119B/SPPL3 (rs11829037, p = 1.50×10(-14)), and ICOSLG/AIRE (rs113459440, p = 1.54×10(-08)) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process.

The metabochip, a custom genotyping array for genetic studies of metabolic, cardiovascular, and anthropometric traits.

Genome-wide association studies have identified hundreds of loci for type 2 diabetes, coronary artery disease and myocardial infarction, as well as for related traits such as body mass index, glucose and insulin levels, lipid levels, and blood pressure. These studies also have pointed to thousands of loci with promising but not yet compelling association evidence. To establish association at additional loci and to characterize the genome-wide significant loci by fine-mapping, we designed the "Metabochip," a custom genotyping array that assays nearly 200,000 SNP markers. Here, we describe the Metabochip and its component SNP sets, evaluate its performance in capturing variation across the allele-frequency spectrum, describe solutions to methodological challenges commonly encountered in its analysis, and evaluate its performance as a platform for genotype imputation. The metabochip achieves dramatic cost efficiencies compared to designing single-trait follow-up reagents, and provides the opportunity to compare results across a range of related traits. The metabochip and similar custom genotyping arrays offer a powerful and cost-effective approach to follow-up large-scale genotyping and sequencing studies and advance our understanding of the genetic basis of complex human diseases and traits.

Fine mapping of five loci associated with low-density lipoprotein cholesterol detects variants that double the explained heritability.

Complex trait genome-wide association studies (GWAS) provide an efficient strategy for evaluating large numbers of common variants in large numbers of individuals and for identifying trait-associated variants. Nevertheless, GWAS often leave much of the trait heritability unexplained. We hypothesized that some of this unexplained heritability might be due to common and rare variants that reside in GWAS identified loci but lack appropriate proxies in modern genotyping arrays. To assess this hypothesis, we re-examined 7 genes (APOE, APOC1, APOC2, SORT1, LDLR, APOB, and PCSK9) in 5 loci associated with low-density lipoprotein cholesterol (LDL-C) in multiple GWAS. For each gene, we first catalogued genetic variation by re-sequencing 256 Sardinian individuals with extreme LDL-C values. Next, we genotyped variants identified by us and by the 1000 Genomes Project (totaling 3,277 SNPs) in 5,524 volunteers. We found that in one locus (PCSK9) the GWAS signal could be explained by a previously described low-frequency variant and that in three loci (PCSK9, APOE, and LDLR) there were additional variants independently associated with LDL-C, including a novel and rare LDLR variant that seems specific to Sardinians. Overall, this more detailed assessment of SNP variation in these loci increased estimates of the heritability of LDL-C accounted for by these genes from 3.1% to 6.5%. All association signals and the heritability estimates were successfully confirmed in a sample of ∼10,000 Finnish and Norwegian individuals. Our results thus suggest that focusing on variants accessible via GWAS can lead to clear underestimates of the trait heritability explained by a set of loci. Further, our results suggest that, as prelude to large-scale sequencing efforts, targeted re-sequencing efforts paired with large-scale genotyping will increase estimates of complex trait heritability explained by known loci.

The role of genetically predicted serum iron levels on neurodegenerative and cardiovascular traits.

Iron is an essential mineral that supports numerous biological functions. Studies have reported associations between iron dysregulation and certain cardiovascular and neurodegenerative diseases, but the direction of influence is not clear. Our goal was to use computational approaches to better understand the role of genetically predicted iron levels on disease risk. We meta-analyzed genome-wide association study summary statistics for serum iron levels from two cohorts and two previous meta-analyses. We then obtained summary statistics from 11 neurodegenerative, cerebrovascular, cardiovascular or lipid traits to assess global and regional genetic correlation between iron levels and these traits. We used two-sample Mendelian randomization (MR) to estimate causal effects. Sex-stratified analyses were also carried out to identify effects potentially differing by sex. Overall, we identified three significant global correlations between iron levels and (i) coronary heart disease, (ii) triglycerides, and (iii) high-density lipoprotein (HDL) cholesterol levels. A total of 194 genomic regions had significant (after correction for multiple testing) local correlations between iron levels and the 11 tested traits. MR analysis revealed two potential causal relationships, between genetically predicted iron levels and (i) total cholesterol or (ii) non-HDL cholesterol. Sex-stratified analyses suggested a potential protective effect of iron levels on Parkinson's disease risk in females, but not in males. Our results will contribute to a better understanding of the genetic basis underlying iron in cardiovascular and neurological health in aging, and to the eventual identification of new preventive interventions or therapeutic avenues for diseases which affect women and men worldwide.

Genome-wide association study of vitamin B6, vitamin B12, folate, and homocysteine blood concentrations.

The B vitamins are components of one-carbon metabolism (OCM) that contribute to DNA synthesis and methylation. Homocysteine, a by-product of OCM, has been associated with coronary heart disease, stroke and neurological disease. To investigate genetic factors that affect circulating vitamin B6, vitamin B12, folate and homocysteine, a genome-wide association analysis was conducted in the InCHIANTI (N = 1175), SardiNIA (N = 1115), and BLSA (N = 640) studies. The top loci were replicated in an independent sample of 687 participants in the Progetto Nutrizione study. Polymorphisms in the ALPL gene (rs4654748, p = 8.30 x 10(-18)) were associated with vitamin B6 and FUT2 (rs602662, [corrected] p = 2.83 x 10(-20)) with vitamin B12 serum levels. The association of MTHFR, a gene consistently associated with homocysteine, was confirmed in this meta-analysis. The ALPL gene likely influences the catabolism of vitamin B6 while FUT2 interferes with absorption of vitamin B12. These findings highlight mechanisms that affect vitamin B6, vitamin B12 and homocysteine serum levels.

Common variants in the SLCO1B3 locus are associated with bilirubin levels and unconjugated hyperbilirubinemia.

Bilirubin, resulting largely from the turnover of hemoglobin, is found in the plasma in two main forms: unconjugated or conjugated with glucuronic acid. Unconjugated bilirubin is transported into hepatocytes. There, it is glucuronidated by UGT1A1 and secreted into the bile canaliculi. We report a genome wide association scan in 4300 Sardinian individuals for total serum bilirubin levels. In addition to the two known loci previously involved in the regulation of bilirubin levels, UGT1A1 (P = 6.2 x 10(-62)) and G6PD (P = 2.5 x 10(-8)), we observed a strong association on chromosome 12 within the SLCO1B3 gene (P = 3.9 x 10(-9)). Our findings were replicated in an independent sample of 1860 Sardinians and in 832 subjects from the Old Order Amish (combined P < 5 x 10(-14)). We also show that SLC01B3 variants contribute to idiopathic mild unconjugated hyperbilirubinemia. Thus, SLC01B3 appears to be involved in the regulation of serum bilirubin levels in healthy individuals and in some bilirubin-related disorders that are only partially explained by other known gene variants.

Phosphodiesterase 8B gene variants are associated with serum TSH levels and thyroid function.

Thyroid-stimulating hormone (TSH) controls thyroid growth and hormone secretion through binding to its G protein-coupled receptor (TSHR) and production of cyclic AMP (cAMP). Serum TSH is a sensitive indicator of thyroid function, and overt abnormalities in thyroid function lead to common endocrine disorders affecting approximately 10% of individuals over a life span. By genotyping 362,129 SNPs in 4,300 Sardinians, we identified a strong association (p = 1.3 x 10(-11)) between alleles of rs4704397 and circulating TSH levels; each additional copy of the minor A allele was associated with an increase of 0.13 muIU/ml in TSH. The single-nucleotide polymorphism (SNP) is located in intron 1 of PDE8B, encoding a high-affinity cAMP-specific phosphodiesterase. The association was replicated in 4,158 individuals, including additional Sardinians and two genetically distant cohorts from Tuscany and the Old Order Amish (overall p value = 1.9 x 10(-20)). In addition to association of TSH levels with SNPs in PDE8B, our genome scan provided evidence for association with PDE10A and several biologically interesting candidates in a focused analysis of 24 genes. In particular, we found evidence for association of TSH levels with SNPs in the THRB (rs1505287, p = 7.3 x 10(-5)), GNAQ (rs10512065, p = 2.0 x 10(-4)), TG (rs2252696, p = 2.2 x 10(-3)), POU1F1 (rs1976324, p = 3.9 x 10(-3)), PDE4D (rs27178, p = 8.3 x 10(-3)), and TSHR (rs4903957, p = 8.6 x 10(-3)) loci. Overall, the results suggest a primary effect of PDE8B variants on cAMP levels in the thyroid. This would affect production of T4 and T3 and feedback to alter TSH release by the pituitary. PDE8B may thus provide a candidate target for the treatment of thyroid dysfunction.

Neuroticism, depressive symptoms, and serum BDNF.

OBJECTIVE: Animal models and clinical studies suggest that brain-derived neurotrophic factor (BDNF) is involved in the pathophysiology of depression. We test whether serum and plasma levels of BDNF are associated with trait neuroticism and its facets and with state measures of depressive symptoms. METHODS: In a community-based cohort (N = 2099), we measured serum and plasma BDNF concentrations and administered the Revised NEO Personality Inventory and the Center for Epidemiological Studies Depression Scale. Covariates included age, sex, cigarette smoking, obesity, and antidepressant use. RESULTS: Serum BDNF concentrations were inversely related to neuroticism (r = -0.074, p < .001), in particular the depression facet (r = -0.08, p < .001). Lower BDNF concentrations were also associated with severe depressive symptoms (Center for Epidemiological Studies Depression Scale ≥ 28; odds ratio = 0.906; 95% confidence interval = 0.851-0.965). The association of serum BDNF with neuroticism was independent of depressive symptoms, indicating that serum BDNF might represent a biological correlate of neuroticism and not just of transient depressive states. Plasma BDNF was not associated with measures of depression. CONCLUSIONS: Our study suggests that lower serum BDNF is associated with both a dispositional vulnerability to depression and acute depressive states in the general population.

Genome-wide association study shows BCL11A associated with persistent fetal hemoglobin and amelioration of the phenotype of beta-thalassemia.

beta-Thalassemia and sickle cell disease both display a great deal of phenotypic heterogeneity, despite being generally thought of as simple Mendelian diseases. The reasons for this are not well understood, although the level of fetal hemoglobin (HbF) is one well characterized ameliorating factor in both of these conditions. To better understand the genetic basis of this heterogeneity, we carried out genome-wide scans with 362,129 common SNPs on 4,305 Sardinians to look for genetic linkage and association with HbF levels, as well as other red blood cell-related traits. Among major variants affecting HbF levels, SNP rs11886868 in the BCL11A gene was strongly associated with this trait (P < 10(-35)). The C allele frequency was significantly higher in Sardinian individuals with elevated HbF levels, detected by screening for beta-thalassemia, and patients with attenuated forms of beta-thalassemia vs. those with thalassemia major. We also show that the same BCL11A variant is strongly associated with HbF levels in a large cohort of sickle cell patients. These results indicate that BCL11A variants, by modulating HbF levels, act as an important ameliorating factor of the beta-thalassemia phenotype, and it is likely they could help ameliorate other hemoglobin disorders. We expect our findings will help to characterize the molecular mechanisms of fetal globin regulation and could eventually contribute to the development of new therapeutic approaches for beta-thalassemia and sickle cell anemia.

The GLUT9 gene is associated with serum uric acid levels in Sardinia and Chianti cohorts.

High serum uric acid levels elevate pro-inflammatory-state gout crystal arthropathy and place individuals at high risk for cardiovascular morbidity and mortality. Genome-wide scans in the genetically isolated Sardinian population identified variants associated with serum uric acid levels as a quantitative trait. They mapped within GLUT9, a Chromosome 4 glucose transporter gene predominantly expressed in liver and kidney. SNP rs6855911 showed the strongest association (p = 1.84 x 10(-16)), along with eight others (p = 7.75 x 10(-16) to 6.05 x 10(-11)). Individuals homozygous for the rare allele of rs6855911 (minor allele frequency = 0.26) had 0.6 mg/dl less uric acid than those homozygous for the common allele; the results were replicated in an unrelated cohort from Tuscany. Our results suggest that polymorphisms in GLUT9 could affect glucose metabolism and uric acid synthesis and/or renal reabsorption, influencing serum uric acid levels over a wide range of values.