Pooled genome wide association detects association upstream of FCRL3 with Graves' disease

Graves disease: latest understanding of pathogenesis and treatment options

Graves disease is the most common cause of hyperthyroidism in iodine-sufficient areas. The main responsible mechanism is related to autoantibodies that bind and activate the thyrotropin receptor (TSHR). Although Graves hyperthyroidism is relatively common, no causal treatment options are available. Established treatment modalities are antithyroid drugs, which reduce thyroid hormone synthesis, radioactive iodine and surgery. However, emerging drugs that target the main autoantigen (monoclonal antibodies, small molecules, peptides) or block the immune pathway have been recently tested in clinical trials. Graves disease can involve the thyroid exclusively or it can be associated with extrathyroidal manifestations, among which Graves orbitopathy is the most common. The presence of Graves orbitopathy can change the management of the disease. An established treatment for moderate-to-severe Graves orbitopathy is intravenous glucocorticoids. However, recent advances in understanding the pathogenesis of Graves orbitopathy have allowed the development of new target-based therapies by blocking pro-inflammatory cytokine receptors, lymphocytic infiltration or the insulin-like growth factor 1 receptor (IGF1R), with several clinical trials providing promising results. This article reviews the new discoveries in the pathogenesis of Graves hyperthyroidism and Graves orbitopathy that offer several important tools in disease management.

Uncovering the shared genetic components of thyroid disorders and reproductive health

OBJECTIVE

The aim of the study is to map the shared genetic component and relationships between thyroid and reproductive health traits to improve the understanding of the interplay between those domains.

DESIGN

A large-scale genetic analysis of thyroid traits (hyper- and hypothyroidism, and thyroid-stimulating hormone levels) was conducted in up to 743 088 individuals of European ancestry from various cohorts.

METHODS

We evaluated genetic associations using genome-wide association study (GWAS) meta-analysis, GWAS Catalog lookup, gene prioritization, mouse phenotype lookup, and genetic correlation analysis.

RESULTS

GWAS meta-analysis results for thyroid phenotypes showed that 50 lead variants out of 253 (including 5/52 of the novel hits) were linked to reproductive health in previous literature. Genetic correlation analyses revealed significant correlations between hypothyroidism and reproductive phenotypes. The results showed that 31.9% of thyroid-associated genes also had an impact on reproductive phenotypes, with the most affected functions being related to genitourinary tract issues.

CONCLUSIONS

The study discovers novel genetic loci linked to thyroid phenotypes and highlights the shared genetic determinants between thyroid function and reproductive health, providing evidence for the genetic pleiotropy and shared biological mechanisms between these traits in both sexes.

Estimation of genetic variation in vitiligo associated genes: Population genomics perspective

BACKGROUND

Vitiligo is an auto-immune progressive depigmentation disorder of the skin due to loss of melanocytes. Genetic risk is one of the important factors for development of vitiligo. Preponderance of vitiligo in certain ethnicities is known which can be analysed by understanding the distribution of allele frequencies across normal populations. Earlier GWAS identified 108 risk alleles for vitiligo in Europeans and East Asians. In this study, 64 of these risk alleles were used for analysing their enrichment and depletion across populations (1000 Genomes Project and IndiGen) with reference to 1000 Genomes dataset. Genetic risk scores were calculated and Fisher's exact test was performed to understand statistical significance of their variation in each population with respect to 1000 Genomes dataset as reference. In addition to SNPs reported in GWAS, significant variation in allele frequencies of 1079 vitiligo-related genes were also analysed. Two-tailed Chi-square test and Bonferroni's multiple adjustment values along with fixation index (≥ 0.5) and minimum allele frequency (≥ 0.05) were calculated and used to prioritise the variants based on pairwise comparison across populations.

RESULTS

Risk alleles rs1043101 and rs10768122 belong to 3 prime UTR of glutamate receptor gene SLC1A2 are found to be highly enriched in the South Asian population when compared with the 'global normal' population. Intron variant rs4766578 (ATXN2) was found to be deleted in SAS, EAS and AFR and enriched in EUR and AMR1. This risk allele is found to be under positive selection in SAS, AMR1 and EUR. From the ancillary vitiligo gene list, nonsynonymous variant rs16891982 was found to be enriched in the European and the Admixed American populations and depleted in all others. rs2279238 and rs11039155 belonging to the LXR-α gene involved in regulation of metalloproteinase 2 and 9 (melanocyte precursors) were found to be associated with vitiligo in the North Indian population (in earlier study).

CONCLUSION

The differential enrichment/depletion profile of the risk alleles provides insight into the underlying inter-population variations. This would provide clues towards prioritisation of SNPs associated with vitiligo thereby elucidating its preponderance in different ethnic groups.

Keratoconus International Consortium (KIC)- advancing keratoconus research

CLINICAL RELEVANCE

The Keratoconus International Consortium (KIC) will allow better understanding of keratoconus.

BACKGROUND

Keratoconus is a disorder characterised by corneal elevation and thinning, leading to reduced vision. The current gaps in understanding of this disease will be discussed and the need for a multi-pronged and multi-centre engagement to enhance our understanding of keratoconus will be highlighted.

DESIGN

KIC has been established to address the gaps in our understanding of keratoconus with the aim of collecting baseline as well as longitudinal data on several fields.

PARTICIPANTS

Keratoconus and control (no corneal condition) subjects from different sites globally will be recruited in the study.

METHODS

KIC collects data using an online, secure database, which enables standardised data collection at member sites. Data fields collected include medical history, clinical features, quality of life and economic burden questionnaires and possible genetic sample collection from patients of different ethnicities across different geographical locations.

RESULTS

There are currently 40 Australian and international clinics or hospital departments who have joined the KIC. Baseline data has so far been collected on 1130 keratoconus patients and indicates a median age of 29.70 years with 61% being male. A total of 15.3% report a positive family history of keratoconus and 57.7% self-report a history of frequent eye rubbing.

CONCLUSION

The strength of this consortium is its international, collaborative design and use of a common data collection tool. Inclusion and analyses of cross-sectional and longitudinal data will help answer many questions that remain in keratoconus, including factors affecting progression and treatment outcomes.

Imputation of the major histocompatibility complex region identifies major independent variants associated with bullous pemphigoid and dermatomyositis in Han Chinese

As autoimmune skin diseases, both bullous pemphigoid (BP) and dermatomyositis (DM) show significant associations with the major histocompatibility complex (MHC) region. In fact, the coexistence of BP and DM has been previously reported. Therefore, we hypothesized that there may be a potential genetic correlation between BP and DM. Based on data for 312 BP patients, 128 DM patients, and 6793 healthy control subjects, in the MHC region, we imputed single-nucleotide polymorphisms (SNP), insertions and deletions (INDEL), and copy number variations (CNV) using the 1KGP phase 3 dataset and amino acids (AA) and SNP using a Han-MHC reference database. An association study revealed the most significant SNP associated with BP, namely, rs580921 (p = 1.06E-08, odds ratio [OR] = 1.61, 95% confidence interval [CI] = 1.37-1.90), which is located in the C6orf10 gene, and the most significant classic human leukocyte antigen (HLA) allele associated with DM, namely, HLA-DPB1*1701 (p = 6.56E-10, OR = 3.61, 95% CI = 2.40-5.42). Further stepwise regression analyses with rs580921 identified a threonine at position 163 of the HLA-B gene as a new independent disease-associated AA, and HLA-DPB1*1701 indicated that no loci were significant. Three-dimensional ribbon models revealed that the HLA-B AA position 163 (p = 3.93E-07, OR = 1.64, 95% CI = 1.35-1.98) located in the α2 domain of the HLA-B molecule was involved in the process of specific antigen presentation. The calculations showed that there was no significant genetic correlation between BP and DM. Our study identified three significant loci in the MHC region, proving that the HLA region was significantly correlated with BP and DM separately. Our research highlights the key role of the MHC region in disease susceptibility.

Genetic association and differential expression of HLAComplexGroup lncRNAs in pemphigus

BACKGROUND

Pemphigus is a group of bullous diseases characterized by acantholysis and skin blisters. As for other autoimmune diseases, the strongest genetic associations found so far for pemphigus foliaceus (PF) and vulgaris (PV) are with alleles of HLA genes. However, apart from protein-coding genes, the MHC region includes a set of poorly explored long non-coding RNA (lncRNA) genes, the HLA complex group (HCG).

OBJECTIVES

To investigate if HCG lncRNA alleles are associated with pemphigus susceptibility.

METHODS AND RESULTS

We analyzed SNPs in 13 HCG lncRNA genes, both in PV (Germany: 241 patients; 1,188 controls) and endemic PF (Brazil: 227 patients; 194 controls), applying multivariate logistic regression. We found 55 associations with PV (p_corr < 0.01) and nine with endemic PF (p_corr < 0.05), the majority located in TSBP1-AS1 (which includes HCG23) and HCG27 lncRNA genes, independently of HLA alleles previously associated with pemphigus. The association of TSBP1-AS1 rs3129949*A allele was further replicated in sporadic PF (p = 0.027, OR = 0.054; 75 patients and 150 controls, all from Germany). Next, we evaluated the expression levels of TSBP1-AS1, TSBP1, HCG23, and HCG27 in blood mononuclear cells of Brazilian patients and controls. HCG27 was upregulated in endemic PF (p = 0.035, log₂ FC = 1.3), while TSBP1-AS1 was downregulated in PV (p = 0.029, log₂ FC = -1.29). The same expression patterns were also seen in cultured keratinocytes stimulated with IgG antibodies from patients and controls from Germany. TSBP1 mRNA levels were also decreased in endemic PF blood cells (p = 0.042, log₂ FC = -2.14). TSBP1-AS1 and HCG27 were also observed downregulated in CD19⁺ cells of endemic PF (p < 0.01, log₂ FC = -0.226 and -0.46 respectively).

CONCLUSIONS

HCG lncRNAs are associated with susceptibility to pemphigus, being TSBP1-AS1 and HCG27 also differentially expressed in distinct cell populations. These results suggest a role for HCG lncRNAs in pemphigus autoimmunity.

Predisposition to Graves' disease and Graves' ophthalmopathy by genetic variants of IL2RA

Previous studies have identified that Th17/Treg cells were involved in the occurrence and development of Graves' disease (GD). This study aimed at clarifying the association between GD susceptibility and nine single nucleotide polymorphisms (SNPs) of Th17/Treg cell-related genes, including IL2RA, miR27a, miR182, and FoxO1. A two-stage association study was performed in 650 GD patients and 1300 healthy controls. PCR-RFLP assays, real-time PCR, and ELISA were performed. In the first stage, association analysis has identified that IL2RA/rs3118470 TT genotype (Pc = 0.027, OR = 1.688) and IL2RA/rs2104286 AA genotype (Pc = 0.027, OR = 1.658) has significantly increased frequencies in patients with GD than control subjects. In the second stage, the result of rs2104286 was consistent with the first-stage results (AA genotype: Pc = 0.006, OR = 1.618). The combined data showed that IL2RA/rs2104286 AA genotype had increased frequencies in patients with GD (Pc = 8.772 × 10^-6, OR = 1.636). Stratification analysis also revealed that rs2104286 AA genotype was significantly associated with Graves' ophthalmopathy (GO) susceptibility (Pc = 9.150 × 10^-4, OR = 1.851). Functional studies showed that carriers of the rs2104286 AA genotype had lower IL2RA mRNA expression than AG genotype carriers (P = 0.021). Cytokine analyses revealed that the rs2104286 AA genotype individuals had lower IL-10 levels (P = 0.015) and increased IL-17 levels than AG genotype carriers (P = 1.467 × 10^-4). In conclusion, our findings suggested that IL2RA/rs2104286 was associated with GD and GO susceptibility in Southwest Chinese Han population, which may be involved in the occurrence of GD and GO by affecting the mRNA expression of IL2RA gene and the cytokine production. KEY MESSAGES: We identified that IL2RA/rs2104286 locus contributed to the predisposition of Graves' disease (GD) and Graves' ophthalmopathy (GO). Functional analyses suggested that IL2RA/rs2104286 may participate in the occurrence of GD and GO by affecting the mRNA expression of IL2RA and cytokine (IL-10 and IL-17) secretion. We found that IL2RA (rs3118470, rs7093069), miR27a/rs895819, miR182/rs76481776, and FoxO1 (rs2297626, rs17592236, rs9549241, rs12585277) loci polymorphisms were not associated with GD susceptibility.

Investigating the GWAS-Implicated Loci for Rheumatoid Arthritis in the Pakistani Population

Rheumatoid arthritis (RA) is a complex and multifactorial autoimmune disorder with the involvement of multiple genetic and environmental factors. Genome-wide association studies (GWAS) have identified more than 50 RA genetic loci in European populations. Given the anticipated overlap of RA-relevant genes and pathways across different ethnic groups, we sought to replicate 58 GWAS-implicated SNPs reported in Europeans in Pakistani subjects. 1,959 unrelated subjects comprising 1,222 RA cases and 737 controls were collected from three rheumatology facilities in Pakistan. Genotyping was performed using iPLEX or TaqMan® methods. A total of 50 SNPs were included in the final association analysis after excluding those that failed assay design/run or postrun QC analysis. Fourteen SNPs (LINC00824/rs1516971, PADI4/rs2240336, CEP57/rs4409785, CTLA4/rs3087243, STAT4/rs13426947, HLA-B/MICA/rs2596565, C5orf30/rs26232, CCL21/rs951005, GATA3/rs2275806, VPS37C/rs595158, HLA-DRB1/rs660895, EOMES/rs3806624, SPRED2/rs934734, and RUNX1/rs9979383) were replicated in our Pakistani sample at false discovery rate (FDR) of <0.20 with nominal p values ranging from 4.73E-06 to 3.48E-02. Our results indicate that several RA susceptibility loci are shared between Pakistani and European populations, supporting the role of common genes/pathways.

Immunological Drivers in Graves' Disease: NK Cells as a Master Switcher

Graves' disease (GD) is a common autoimmune cause of hyperthyroidism, which is eventually related to the generation of IgG antibodies stimulating the thyrotropin receptor. Clinical manifestations of the disease reflect hyperstimulation of the gland, causing thyrocyte hyperplasia (goiter) and excessive thyroid hormone synthesis (hyperthyroidism). The above clinical manifestations are preceded by still partially unraveled pathogenic actions governed by the induction of aberrant phenotype/functions of immune cells. In this review article we investigated the potential contribution of natural killer (NK) cells, based on literature analysis, to discuss the bidirectional interplay with thyroid hormones (TH) in GD progression. We analyzed cellular and molecular NK-cell associated mechanisms potentially impacting on GD, in a view of identification of the main NK-cell subset with highest immunoregulatory role.

Linkage Analysis of the Chromosome 5q31-33 Region Identifies JAKMIP2 as a Risk Factor for Graves' Disease in the Chinese Han Population

Background

This study aimed to investigate susceptibility to Graves’s disease and the association with the 5q32–33.1 region on chromosome 5 in a Chinese Han population.

Material/Methods

Eighty Chinese Han multiplex families included first-degree and second-degree relatives with Graves’ disease. Eight microsatellite markers on chromosome 5 at the 5q32–33.1 region underwent linkage analysis and the association between the regions D5S1480–D5S2014 were studied.

Results

The maximal heterogeneity logarithm of the odds (HLOD) score of D5S2090 was 4.29 (α=0.42) and of D5S2014 was 4.01 (α=0.34). A nonparametric linkage (NPL) score of 3.14 (P<0.001) was found for D5S2014. The D5S1480–D5S2014 region on chromosome 5 was associated with Graves’ disease, with eight haplotype domains. There were significant differences in the sixth and eighth haplotype domains between patients with Graves’ disease compared with normal individuals. Tagging single nucleotide polymorphisms (SNPs) of the sixth and eighth haplotype domains showed that individuals with SNP62 (rs12653715 G/C) who were GG homozygous had a significantly increased risk of Graves’ disease compared GC heterozygous or CC homozygous individuals. The transmission disequilibrium test (TDT) indicated that SNP62 (rs12653715) and SNP63 (rs12653081) loci in the Janus kinase and microtubule interacting protein 2 (JAKMIP2) gene showed dominant transmission from heterozygous parents to the affected offspring, and SNPs in the secretoglobin family 3A member 2 (SCGB3A2) gene showed no transmission disequilibrium. The haplotype JAKMIP2-1 was identified as being particularly significant.

Conclusions

JAKMIP2 gene polymorphism require further study as potential risk factors for Graves’ disease in the Chinese Han population.

Correlation of TSHR and CTLA-4 Single Nucleotide Polymorphisms with Graves Disease

This study was designed to explore the association between Graves disease (GD) and thyroid-stimulating hormone receptor (TSHR) and cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4) single nucleotide polymorphisms (SNPs). We studied a total of 1217 subjects from a Han population in northern Anhui province in China. Six SNPs within TSHR (rs179247, rs12101261, rs2284722, rs4903964, rs2300525, and rs17111394) and four SNPs within CTLA-4 (rs10197319, rs231726, rs231804, and rs1024161) were genotyped via a Taqman probe technique using a Fluidigm EP1 platform. The TSHR alleles rs179247-G, rs12101261-C, and rs4903964-G were negatively correlated with GD, whereas the rs2284722-A and rs17111394-C alleles were positively correlated with GD. Analyzing TSHR SNPs at rs179247, rs2284722, rs12101261, and rs4903964 yielded 8 different haplotypes. There were positive correlations between GD risk and the haplotypes AGTA and AATA (OR = 1.27, 95%CI = 1.07-1.50, P = 0.005; OR = 1.45, 95%CI = 1.21-1.75, P < 0.001, respectively). There were negative correlations between GD risk and the haplotype GGCG (OR = 0.56, 95%CI = 0.46-0.67, P < 0.001). With respect to haplotypes based on SNPs at the TSHR rs2300525 and rs17111394 loci, the CC haplotype was positively correlated with GD risk (OR = 1.32, 95%CI = 1.08-1.60, P = 0.006). Analyzing CTLA-4 SNPs at rs231804, rs1024161, and rs231726 yielded four haplotypes, of which AAA was positively correlated with GD risk (OR = 1.21, 95%CI = 1.02-1.43, P = 0.029). Polymorphisms at rs179247, rs12101261, rs2284722, rs4903964, and rs17111394 were associated with GD susceptibility. Haplotypes of both TSHR and CTLA-4 were additionally related to GD risk.

Identification of novel alleles associated with insulin resistance in childhood obesity using pooled-DNA genome-wide association study approach

Background:

Recently, we witnessed great progress in the discovery of genetic variants associated with obesity and type 2 diabetes (T2D), especially in adults. Much less is known regarding genetic variants associated with insulin resistance (IR). We hypothesized that novel IR genes could be efficiently detected in a population of obese children and adolescents who may not exhibit comorbidities and other confounding factors.

Objectives:

This study aimed to determine whether a genome-wide association study (GWAS), using a DNA-pooling approach, could identify novel genes associated with IR.

Subjects:

The pooled-DNA GWAS analysis included Slovenian obese children and adolescents with and without IR matched for body mass index, gender and age. A replication study was conducted in another independent cohort with or without IR.

Methods:

For the pooled-DNA GWAS, we used HumanOmni5-Quad SNP array (Illumina). Allele frequency distributions were compared with modified t-tests and χ²-tests and ranked using PLINK. Top single nucleotide polymorphisms (SNPs) were validated using individual genotyping by high-resolution melting analysis and TaqMan assay.

Results:

We identified five top-ranking SNPs from the pooled-DNA GWAS analysis within the ECE1, IL1R2, GNPDA1, HLA-J and PYGB loci. All except SNP rs9261108 (HLA-J locus) were confirmed in the validation phase using individual genotyping. The SNP rs2258617 within PYGB remained statistically significant for both recessive and additive models in both cohorts and in a merged analysis of both cohorts and present the strongest novel candidate gene for IR.

Conclusion:

We report for the first time a pooled-DNA GWAS approach to identify five novel SNPs or genes for IR in a paediatric population. The four loci confirmed in the second validation phase study warrant further studies, especially the strongest SNP rs2258617 within PYGB, and provide targets for further basic research of IR mechanisms and for the development of potential new IR and T2D therapies.