Novel associations for hypothyroidism include known autoimmune risk loci

Hypothyroidism is the most common thyroid disorder, affecting about 5% of the general population. Here we present the first large genome-wide association study of hypothyroidism, in 2,564 cases and 24,448 controls from the customer base of 23andMe, Inc., a personal genetics company. We identify four genome-wide significant associations, two of which are well known to be involved with a large spectrum of autoimmune diseases: rs6679677 near _PTPN22_ and rs3184504 in _SH2B3_ (p-values 3.5e-13 and 3.0e-11, respectively). We also report associations with rs4915077 near _VAV3_ (p-value 8.3e-11), another gene involved in immune function, and rs965513 near _FOXE1_ (p-value 3.1e-14). Of these, the association with _PTPN22_ confirms a recent small candidate gene study, and _FOXE1_ was previously known to be associated with thyroid-stimulating hormone (TSH) levels. Although _SH2B3_ has been previously linked with a number of autoimmune diseases, this is the first report of its association with thyroid disease. The _VAV3_ association is novel. These results suggest heterogeneity in the genetic etiology of hypothyroidism, implicating genes involved in both autoimmune disorders and thyroid function. Using a genetic risk profile score based on the top association from each of the four genome-wide significant regions in our study, the relative risk between the highest and lowest deciles of genetic risk is 2.1

Specificity of the STAT4 Genetic Association for Severe Disease Manifestations of Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. A polymorphism in the STAT4 gene has recently been established as a risk factor for SLE, but the relationship with specific SLE subphenotypes has not been studied. We studied 137 SNPs in the STAT4 region genotyped in 4 independent SLE case series (total n = 1398) and 2560 healthy controls, along with clinical data for the cases. Using conditional testing, we confirmed the most significant STAT4 haplotype for SLE risk. We then studied a SNP marking this haplotype for association with specific SLE subphenotypes, including autoantibody production, nephritis, arthritis, mucocutaneous manifestations, and age at diagnosis. To prevent possible type-I errors from population stratification, we reanalyzed the data using a subset of subjects determined to be most homogeneous based on principal components analysis of genome-wide data. We confirmed that four SNPs in very high LD (r2 = 0.94 to 0.99) were most strongly associated with SLE, and there was no compelling evidence for additional SLE risk loci in the STAT4 region. SNP rs7574865 marking this haplotype had a minor allele frequency (MAF) = 31.1% in SLE cases compared with 22.5% in controls (OR = 1.56, p = 10−16). This SNP was more strongly associated with SLE characterized by double-stranded DNA autoantibodies (MAF = 35.1%, OR = 1.86, p<10−19), nephritis (MAF = 34.3%, OR = 1.80, p<10−11), and age at diagnosis<30 years (MAF = 33.8%, OR = 1.77, p<10−13). An association with severe nephritis was even more striking (MAF = 39.2%, OR = 2.35, p<10−4 in the homogeneous subset of subjects). In contrast, STAT4 was less strongly associated with oral ulcers, a manifestation associated with milder disease. We conclude that this common polymorphism of STAT4 contributes to the phenotypic heterogeneity of SLE, predisposing specifically to more severe disease

Proteins encoded in genomic regions associated with immune-mediated disease physically interact and suggest underlying biology

Genome-wide association studies have uncovered hundreds of DNA changes associated with complex disease. The ultimate promise of these studies is the understanding of disease biology; this goal, however, is not easily achieved because each disease has yielded numerous associations, each one pointing to a region of the genome, rather than a specific causal mutation. Presumably, the causal variants affect components of common molecular processes, and a first step in understanding the disease biology perturbed in patients is to identify connections among regions associated to disease. Since it has been reported in numerous Mendelian diseases that protein products of causal genes tend to physically bind each other, we chose to approach this problem using known protein–protein interactions to test whether any of the products of genes in five complex trait-associated loci bind each other. We applied several permutation methods and find robustly significant connectivity within four of the traits. In Crohn's disease and rheumatoid arthritis, we are able to show that these genes are co-expressed and that other proteins emerging in the network are enriched for association to disease. These findings suggest that, for the complex traits studied here, associated loci contain variants that affect common molecular processes, rather than distinct mechanisms specific to each association.Massachusetts Institute of Technology (MIT IDEA2 Program)Harvard University. Biological and Biomedical Sciences ProgramEunice Kennedy Shriver National Institute of Child Health and Human Development (U.S.) (NICHD RO1 grant HD055150-03)National Institute of Arthritis and Musculoskeletal and Skin Diseases (U.S.) (K08 NIH-NIAMS career development award (AR055688))National Institute of Diabetes and Digestive and Kidney Diseases (U.S.) (DK083756)National Institute of Diabetes and Digestive and Kidney Diseases (U.S.) (DK086502)Denmark. Forskningsradet for Sundhed og SygdomCenter for the Study of Inflammatory Bowel Diseas

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, Crohn's 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 PCPMA<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

On the Wegener granulomatosis associated region on chromosome 6p21.3

BACKGROUND: Wegener granulomatosis (WG) belongs to the heterogeneous group of systemic vasculitides. The multifactorial pathophysiology of WG is supposedly caused by yet unknown environmental influence(s) on the basis of genetic predisposition. The presence of anti-neutrophil cytoplasmic antibodies (ANCA) in the plasma of patients and genetic involvement of the human leukocyte antigen system reflect an autoimmune background of the disease. Strong associations were revealed with WG by markers located in the major histocompatibility complex class II (MHC II) region in the vicinity of human leukocyte antigen (HLA)-DPB1 and the retinoid X receptor B (RXRB) loci. In order to define the involvement of the 6p21.3 region in WG in more detail this previous population-based association study was expanded here to the respective 3.6 megabase encompassing this region on chromosome 6. The RXRB gene was analysed as well as a splice-site variation of the butyrophilin-like (BTNL2) gene which is also located within the respective region. The latter polymorphism has been evaluated here as it appears as a HLA independent susceptibility factor in another granulomatous disorder, sarcoidosis. METHODS: 150–180 German WG patients and a corresponding cohort of healthy controls (n = 100–261) were used in a two-step study. A panel of 94 microsatellites was designed for the initial step using a DNA pooling approach. Markers with significantly differing allele frequencies between patient and control pools were individually genotyped. The RXRB gene was analysed for single strand conformation polymorphisms (SSCP) and restriction fragment length polymorphisms (RFLP). The splice-site polymorphism in the BTNL2 gene was also investigated by RFLP analysis. RESULTS: A previously investigated microsatellite (#1.0.3.7, Santa Cruz genome browser (UCSC) May 2004 Freeze localisation: chr6:31257596-34999883), which was used as a positive control, remained associated throughout the whole two-step approach. Yet, no additional evidence for association of other microsatellite markers was found in the entire investigated region. Analysis of the RXRB gene located in the WG associated region revealed associations of two variations (rs10548957 p(allelic )= 0.02 and rs6531 p(allelic )= 5.20 × 10(-5), OR = 1.88). Several alleles of markers located between HLA-DPB1, SNP rs6531 and microsatellite 1.0.3.7 showed linkage disequilibrium with r(2 )values exceeding 0.10. Significant differences were not demonstrable for the sarcoidosis associated splice-site variation (rs2076530 p(allelic )= 0.80) in our WG cohort. CONCLUSION: Since a microsatellite flanking the RXRB gene and two intragenic polymorphisms are associated significantly with WG on chromosome 6p21.3, further investigations should be focussed on extensive fine-mapping in this region by densely mapping with additional markers such as SNPs. This strategy may reveal even deeper insights into the genetic contributions of the respective region for the pathogenesis of WG

Hypoxanthine-guanine phosophoribosyltransferase (HPRT) deficiency: Lesch-Nyhan syndrome

Deficiency of hypoxanthine-guanine phosphoribosyltransferase (HPRT) activity is an inborn error of purine metabolism associated with uric acid overproduction and a continuum spectrum of neurological manifestations depending on the degree of the enzymatic deficiency. The prevalence is estimated at 1/380,000 live births in Canada, and 1/235,000 live births in Spain. Uric acid overproduction is present inall HPRT-deficient patients and is associated with lithiasis and gout. Neurological manifestations include severe action dystonia, choreoathetosis, ballismus, cognitive and attention deficit, and self-injurious behaviour. The most severe forms are known as Lesch-Nyhan syndrome (patients are normal at birth and diagnosis can be accomplished when psychomotor delay becomes apparent). Partial HPRT-deficient patients present these symptoms with a different intensity, and in the least severe forms symptoms may be unapparent. Megaloblastic anaemia is also associated with the disease. Inheritance of HPRT deficiency is X-linked recessive, thus males are generally affected and heterozygous female are carriers (usually asymptomatic). Human HPRT is encoded by a single structural gene on the long arm of the X chromosome at Xq26. To date, more than 300 disease-associated mutations in the HPRT1 gene have been identified. The diagnosis is based on clinical and biochemical findings (hyperuricemia and hyperuricosuria associated with psychomotor delay), and enzymatic (HPRT activity determination in haemolysate, intact erythrocytes or fibroblasts) and molecular tests. Molecular diagnosis allows faster and more accurate carrier and prenatal diagnosis. Prenatal diagnosis can be performed with amniotic cells obtained by amniocentesis at about 15–18 weeks' gestation, or chorionic villus cells obtained at about 10–12 weeks' gestation. Uric acid overproduction can be managed by allopurinol treatment. Doses must be carefully adjusted to avoid xanthine lithiasis. The lack of precise understanding of the neurological dysfunction has precluded development of useful therapies. Spasticity, when present, and dystonia can be managed with benzodiazepines and gamma-aminobutyric acid inhibitors such as baclofen. Physical rehabilitation, including management of dysarthria and dysphagia, special devices to enable hand control, appropriate walking aids, and a programme of posture management to prevent deformities are recommended. Self-injurious behaviour must be managed by a combination of physical restraints, behavioural and pharmaceutical treatments

High-Density SNP Screening of the Major Histocompatibility Complex in Systemic Lupus Erythematosus Demonstrates Strong Evidence for Independent Susceptibility Regions

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99×10−16). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53×10−12), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80×10−13. Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE–associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation