GWAS of thyroid stimulating hormone highlights pleiotropic effects and inverse association with thyroid cancer.

Thyroid stimulating hormone (TSH) is critical for normal development and metabolism. To better understand the genetic contribution to TSH levels, we conduct a GWAS meta-analysis at 22.4 million genetic markers in up to 119,715 individuals and identify 74 genome-wide significant loci for TSH, of which 28 are previously unreported. Functional experiments show that the thyroglobulin protein-altering variants P118L and G67S impact thyroglobulin secretion. Phenome-wide association analysis in the UK Biobank demonstrates the pleiotropic effects of TSH-associated variants and a polygenic score for higher TSH levels is associated with a reduced risk of thyroid cancer in the UK Biobank and three other independent studies. Two-sample Mendelian randomization using TSH index variants as instrumental variables suggests a protective effect of higher TSH levels (indicating lower thyroid function) on risk of thyroid cancer and goiter. Our findings highlight the pleiotropic effects of TSH-associated variants on thyroid function and growth of malignant and benign thyroid tumors.

Genome sequencing elucidates Sardinian genetic architecture and augments association analyses for lipid and blood inflammatory markers.

We report ∼17.6 million genetic variants from whole-genome sequencing of 2,120 Sardinians; 22% are absent from previous sequencing-based compilations and are enriched for predicted functional consequences. Furthermore, ∼76,000 variants common in our sample (frequency >5%) are rare elsewhere (<0.5% in the 1000 Genomes Project). We assessed the impact of these variants on circulating lipid levels and five inflammatory biomarkers. We observe 14 signals, including 2 major new loci, for lipid levels and 19 signals, including 2 new loci, for inflammatory markers. The new associations would have been missed in analyses based on 1000 Genomes Project data, underlining the advantages of large-scale sequencing in this founder population.

Whole-genome sequence-based analysis of thyroid function.

Normal thyroid function is essential for health, but its genetic architecture remains poorly understood. Here, for the heritable thyroid traits thyrotropin (TSH) and free thyroxine (FT4), we analyse whole-genome sequence data from the UK10K project (N=2,287). Using additional whole-genome sequence and deeply imputed data sets, we report meta-analysis results for common variants (MAF≥1%) associated with TSH and FT4 (N=16,335). For TSH, we identify a novel variant in SYN2 (MAF=23.5%, P=6.15 × 10(-9)) and a new independent variant in PDE8B (MAF=10.4%, P=5.94 × 10(-14)). For FT4, we report a low-frequency variant near B4GALT6/SLC25A52 (MAF=3.2%, P=1.27 × 10(-9)) tagging a rare TTR variant (MAF=0.4%, P=2.14 × 10(-11)). All common variants explain ≥20% of the variance in TSH and FT4. Analysis of rare variants (MAF<1%) using sequence kernel association testing reveals a novel association with FT4 in NRG1. Our results demonstrate that increased coverage in whole-genome sequence association studies identifies novel variants associated with thyroid function.

The association of thyroid peroxidase antibody risk loci with susceptibility to and phenotype of Graves' disease.

BACKGROUND: Despite great progress, the genetic basis of Graves' disease (GD) remains poorly understood. Recently, a population-based genomewide association study (GWAS) identified five novel loci (ATXN2/SH2B3, MAGI3, BACH2, TPO and KALRN) as significantly associated with the presence of thyroid peroxidase autoantibodies (TPOAbs), whereas several other loci showed suggestive association. METHODS: In this study, we investigated 16 single nucleotide polymorphisms (SNPs) associated with TPOAbs for the association with susceptibility to and phenotype of GD in a cohort of 647 patients with GD and 769 controls from a Polish Caucasian population. RESULTS: SNPs within/near HCP5 (rs3094228, P = 1·6 × 10(-12) , OR = 1·88), MAGI3 (rs1230666, P = 1·9 × 10(-5) , OR = 1·51) and ATXN2/SH2B3 (rs653178, P = 0·0015, OR = 1·28) loci were significantly associated with susceptibility to GD. Allele frequencies differed significantly in subgroups of patients with GD stratified by age of GD onset for HCP5 (P = 0·0014, OR = 1·50) and showed a suggestive difference for MAGI3 (P = 0·0035, OR = 1·50) SNPs. Although rs11675434 located near TPO showed no association with GD susceptibility, it was significantly associated with the presence of clinically evident Graves' ophthalmopathy (GO, P = 5·2 × 10(-5) , OR = 1·64), and this effect was independent from smoking status, age of GD onset and gender. CONCLUSIONS: This is the first study showing an association of the ATXN2/SH2B3 locus with susceptibility to GD. Furthermore, we observed a novel significant association within the HLA region at a SNP located near HCP5 and confirmed the association of the MAGI3 locus with GD susceptibility. HCP5 and MAGI3 SNPs were further correlated with age of GD onset. Finally, we identified TPO as a new susceptibility locus for GO.

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.

Multiethnic meta-analysis of genome-wide association studies in >100 000 subjects identifies 23 fibrinogen-associated Loci but no strong evidence of a causal association between circulating fibrinogen and cardiovascular disease.

BACKGROUND: Estimates of the heritability of plasma fibrinogen concentration, an established predictor of cardiovascular disease, range from 34% to 50%. Genetic variants so far identified by genome-wide association studies explain only a small proportion (<2%) of its variation. METHODS AND RESULTS: We conducted a meta-analysis of 28 genome-wide association studies including >90 000 subjects of European ancestry, the first genome-wide association meta-analysis of fibrinogen levels in 7 studies in blacks totaling 8289 samples, and a genome-wide association study in Hispanics totaling 1366 samples. Evaluation for association of single-nucleotide polymorphisms with clinical outcomes included a total of 40 695 cases and 85 582 controls for coronary artery disease, 4752 cases and 24 030 controls for stroke, and 3208 cases and 46 167 controls for venous thromboembolism. Overall, we identified 24 genome-wide significant (P<5×10(-8)) independent signals in 23 loci, including 15 novel associations, together accounting for 3.7% of plasma fibrinogen variation. Gene-set enrichment analysis highlighted key roles in fibrinogen regulation for the 3 structural fibrinogen genes and pathways related to inflammation, adipocytokines, and thyrotrophin-releasing hormone signaling. Whereas lead single-nucleotide polymorphisms in a few loci were significantly associated with coronary artery disease, the combined effect of all 24 fibrinogen-associated lead single-nucleotide polymorphisms was not significant for coronary artery disease, stroke, or venous thromboembolism. CONCLUSIONS: We identify 23 robustly associated fibrinogen loci, 15 of which are new. Clinical outcome analysis of these loci does not support a causal relationship between circulating levels of fibrinogen and coronary artery disease, stroke, or venous thromboembolism.

The molecular genetic architecture of self-employment.

Economic variables such as income, education, and occupation are known to affect mortality and morbidity, such as cardiovascular disease, and have also been shown to be partly heritable. However, very little is known about which genes influence economic variables, although these genes may have both a direct and an indirect effect on health. We report results from the first large-scale collaboration that studies the molecular genetic architecture of an economic variable-entrepreneurship-that was operationalized using self-employment, a widely-available proxy. Our results suggest that common SNPs when considered jointly explain about half of the narrow-sense heritability of self-employment estimated in twin data (σ(g)(2)/σ(P)(2) = 25%, h(2) = 55%). However, a meta-analysis of genome-wide association studies across sixteen studies comprising 50,627 participants did not identify genome-wide significant SNPs. 58 SNPs with p<10(-5) were tested in a replication sample (n = 3,271), but none replicated. Furthermore, a gene-based test shows that none of the genes that were previously suggested in the literature to influence entrepreneurship reveal significant associations. Finally, SNP-based genetic scores that use results from the meta-analysis capture less than 0.2% of the variance in self-employment in an independent sample (p≥0.039). Our results are consistent with a highly polygenic molecular genetic architecture of self-employment, with many genetic variants of small effect. Although self-employment is a multi-faceted, heavily environmentally influenced, and biologically distal trait, our results are similar to those for other genetically complex and biologically more proximate outcomes, such as height, intelligence, personality, and several diseases.

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.

Meta-analysis of genome-wide association studies in >80 000 subjects identifies multiple loci for C-reactive protein levels.

BACKGROUND: C-reactive protein (CRP) is a heritable marker of chronic inflammation that is strongly associated with cardiovascular disease. We sought to identify genetic variants that are associated with CRP levels. METHODS AND RESULTS: We performed a genome-wide association analysis of CRP in 66 185 participants from 15 population-based studies. We sought replication for the genome-wide significant and suggestive loci in a replication panel comprising 16 540 individuals from 10 independent studies. We found 18 genome-wide significant loci, and we provided evidence of replication for 8 of them. Our results confirm 7 previously known loci and introduce 11 novel loci that are implicated in pathways related to the metabolic syndrome (APOC1, HNF1A, LEPR, GCKR, HNF4A, and PTPN2) or the immune system (CRP, IL6R, NLRP3, IL1F10, and IRF1) or that reside in regions previously not known to play a role in chronic inflammation (PPP1R3B, SALL1, PABPC4, ASCL1, RORA, and BCL7B). We found a significant interaction of body mass index with LEPR (P<2.9×10(-6)). A weighted genetic risk score that was developed to summarize the effect of risk alleles was strongly associated with CRP levels and explained ≈5% of the trait variance; however, there was no evidence for these genetic variants explaining the association of CRP with coronary heart disease. CONCLUSIONS: We identified 18 loci that were associated with CRP levels. Our study highlights immune response and metabolic regulatory pathways involved in the regulation of chronic inflammation.

Identification of a common variant in the TFR2 gene implicated in the physiological regulation of serum iron levels.

The genetic determinants of variation in iron status are actively sought, but remain incompletely understood. Meta-analysis of two genome-wide association (GWA) studies and replication in three independent cohorts was performed to identify genetic loci associated in the general population with serum levels of iron and markers of iron status, including transferrin, ferritin, soluble transferrin receptor (sTfR) and sTfR-ferritin index. We identified and replicated a novel association of a common variant in the type-2 transferrin receptor (TFR2) gene with iron levels, with effect sizes highly consistent across samples. In addition, we identified and replicated an association between the HFE locus and ferritin and confirmed previously reported associations with the TF, TMPRSS6 and HFE genes. The five replicated variants were tested for association with expression levels of the corresponding genes in a publicly available data set of human liver samples, and nominally statistically significant expression differences by genotype were observed for all genes, although only rs3811647 in the TF gene survived the Bonferroni correction for multiple testing. In addition, we measured for the first time the effects of the common variant in TMPRSS6, rs4820268, on hepcidin mRNA in peripheral blood (n = 83 individuals) and on hepcidin levels in urine (n = 529) and observed an association in the same direction, though only borderline significant. These functional findings require confirmation in further studies with larger sample sizes, but they suggest that common variants in TMPRSS6 could modify the hepcidin-iron feedback loop in clinically unaffected individuals, thus making them more susceptible to imbalances of iron homeostasis.

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.

Common variants at ten loci modulate the QT interval duration in the QTSCD Study.

The QT interval, a measure of cardiac repolarization, predisposes to ventricular arrhythmias and sudden cardiac death (SCD) when prolonged or shortened. A common variant in NOS1AP is known to influence repolarization. We analyze genome-wide data from five population-based cohorts (ARIC, KORA, SardiNIA, GenNOVA and HNR) with a total of 15,842 individuals of European ancestry, to confirm the NOS1AP association and identify nine additional loci at P < 5 x 10(-8). Four loci map near the monogenic long-QT syndrome genes KCNQ1, KCNH2, SCN5A and KCNJ2. Two other loci include ATP1B1 and PLN, genes with established electrophysiological function, whereas three map to RNF207, near LITAF and within NDRG4-GINS3-SETD6-CNOT1, respectively, all of which have not previously been implicated in cardiac electrophysiology. These results, together with an accompanying paper from the QTGEN consortium, identify new candidate genes for ventricular arrhythmias and SCD.

Variants of the serotonin transporter gene and NEO-PI-R Neuroticism: No association in the BLSA and SardiNIA samples.

The polymorphism in the serotonin transporter gene promoter region (5-HTTLPR) is by far the most studied variant hypothesized to influence Neuroticism-related personality traits. The results of previous studies have been mixed and appear moderated by the personality questionnaire used. Studies that used the TCI to assess Harm Avoidance or the EPQ to assess Neuroticism have found no association with the 5-HTTLPR. However, studies that used the NEO-PI-R or related instruments (NEO-PI, NEO-FFI) to measure Neuroticism have found some evidence of association. This study examines the association of variants in the serotonin transporter gene in a sample from a genetically isolated population within Sardinia (Italy) that is several times larger than previous samples that used the NEO-PI-R (N = 3,913). The association was also tested in a sample (N = 548) from the Baltimore Longitudinal Study of Aging (BLSA), in which repeated NEO-PI-R assessments were obtained. In the SardiNIA sample, we found no significant association of the 5-HTTLPR genotypes with Neuroticism or its facets (Anxiety, Angry-Hostility, Depression, Self-Consciousness, Impulsiveness, and Vulnerability). In the BLSA sample, we found lower scores on Neuroticism traits for the heterozygous group, which is inconsistent with previous studies. We also examined eight SNPs in the SardiNIA (N = 3,972) and nine SNPs in the BLSA (N = 1,182) that map within or near the serotonin transporter gene (SLC6A4), and found no association. Along with other large studies that used different phenotypic measures and found no association, this study substantially increases the evidence against a link between 5-HTT variants and Neuroticism-related traits.

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.

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.

Common variants in the GDF5-UQCC region are associated with variation in human height.

Identifying genetic variants that influence human height will advance our understanding of skeletal growth and development. Several rare genetic variants have been convincingly and reproducibly associated with height in mendelian syndromes, and common variants in the transcription factor gene HMGA2 are associated with variation in height in the general population. Here we report genome-wide association analyses, using genotyped and imputed markers, of 6,669 individuals from Finland and Sardinia, and follow-up analyses in an additional 28,801 individuals. We show that common variants in the osteoarthritis-associated locus GDF5-UQCC contribute to variation in height with an estimated additive effect of 0.44 cm (overall P < 10(-15)). Our results indicate that there may be a link between the genetic basis of height and osteoarthritis, potentially mediated through alterations in bone growth and development.

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.

IRAK-M is involved in the pathogenesis of early-onset persistent asthma.

Asthma is a multifactorial disease influenced by genetic and environmental factors. In the past decade, several loci and >100 genes have been found to be associated with the disease in at least one population. Among these loci, region 12q13-24 has been implicated in asthma etiology in multiple populations, suggesting that it harbors one or more asthma susceptibility genes. We performed linkage and association analyses by transmission/disequilibrium test and case-control analysis in the candidate region 12q13-24, using the Sardinian founder population, in which limited heterogeneity of pathogenetic alleles for monogenic and complex disorders as well as of environmental conditions should facilitate the study of multifactorial traits. We analyzed our cohort, using a cutoff age of 13 years at asthma onset, and detected significant linkage to a portion of 12q13-24. We identified IRAK-M as the gene contributing to the linkage and showed that it is associated with early-onset persistent asthma. We defined protective and predisposing SNP haplotypes and replicated associations in an outbred Italian population. Sequence analysis in patients found mutations, including inactivating lesions, in the IRAK-M coding region. Immunohistochemistry of lung biopsies showed that IRAK-M is highly expressed in epithelial cells. We report that IRAK-M is involved in the pathogenesis of early-onset persistent asthma. IRAK-M, a negative regulator of the Toll-like receptor/IL-1R pathways, is a master regulator of NF- kappa B and inflammation. Our data suggest a mechanistic link between hyperactivation of the innate immune system and chronic airway inflammation and indicate IRAK-M as a potential target for therapeutic intervention against asthma.