Genetics, Epigenetics, Cellular Immunology, and Gut Microbiota: Emerging Links With Graves’ Disease
Fangyu Zhou1 
Xin Wang2 
Lingjun Wang2 
Xin Sun2 
Guiqin Tan1 
Wenwen Wei1 
Guangbing Zheng1 
Xiaomin Ma1 
Dan Tian1 
Hongsong Yu1*- 1Department of Immunology, Special Key Laboratory of Ocular Diseases of Guizhou Province, Zunyi Medical University, Zunyi, China
- 2School of Basic Medical Sciences, Special Key Laboratory of Gene Detection and Therapy of Guizhou Province, Zunyi Medical University, Zunyi, China
Graves’ disease (GD) is a well-known organ-specific autoimmune disease characterized by hyperthyroidism, goiter, and exophthalmos. The incidence of GD is approximately 2.0–3.0% in China and 0.5–2.0% in Western countries. Due to the complex pathogenesis and etiology of GD, current treatment methods have great side effects that seriously endanger human health. Therefore, it is particularly important to understand the pathogenesis of GD. Various studies have shown that genetics, epigenetics, cellular immunology, and gut microbiota are all involved in the development of GD. Genetically, CD25 gene and VDR gene polymorphisms are involved in the development of GD by increasing the ratio of Th17/Treg cells. Epigenetically, miR-23a-3p and lncRNA-MEG3 lead to Th17/Treg imbalance and participate in the progression of GD. Moreover, commensal microbe deletion can disrupt Th17/Treg balance and participate in the occurrence of GD. The imbalance of Th17/Treg cells induced by genetics, epigenetics, and gut microbiota plays a vital role in the pathogenesis of GD. Therefore, this article reviews the role of genetics, epigenetics, cellular immunology, and gut microbiota in the pathogenic mechanism of GD. This may lead to the development of novel therapeutic strategies and providing promising therapeutic targets.
Introduction
Graves’ disease (GD), also known as toxic diffuse goiter, is one of the most common autoimmune thyroid diseases (AITDs). It is the main cause of hyperthyroidism, and hyperthyroidism syndrome, and the main clinical characteristics include varying degrees of goiter and exophthalmos (Wemeau et al., 2018). GD is believed to represent the autoimmune process of the thyroid, in which irritant autoantibodies combined with the thyroid-stimulating hormone receptor (TSHR) and activate thyroid function, leading to hyperthyroidism (Prabhakar et al., 2003). GD tends to occur in adult females aged 30–50, and the risk in women is six times higher than in men, with an annual morbidity of 20–50 cases per 100,000 people (Smith and Hegedus, 2016). The immunological characteristics of GD are the presence of thyroid-stimulating hormone receptor antibodies (TRAbs), thyroid peroxidase antibodies (TPOAbs), thyroglobulin antibodies (TgAbs), and other autoantibodies in the serum, leading to hyperthyroidism and diffuse thyroid enlargement (McLachlan and Rapoport, 2014).
GD is a T cell-mediated organ-specific autoimmune disorder. The infiltration of T lymphocytes in GD directly destroys the thyroid. Moreover, it can also stimulate B cells to differentiate into plasma cells that secrete antibodies (Li et al., 2019). Study has found that an increased number of CD4+ T cells and plasma cells were observed in GD patients (Ben-Skowronek et al., 2013). Naïve CD4+ T cells can differentiate into T helper 17 (Th17) and regulatory T cells (Tregs) when stimulated by specific antigens (Borst et al., 2018). A number of studies have shown that Th17/Tregs balance is vital to the pathogenesis of GD. Besides GD, Th17/Tregs balance may be also involved in the etiology of other autoimmune diseases, including rheumatoid arthritis (RA), psoriasis, multiple sclerosis, inflammatory bowel disease (Lee, 2018). Therefore, immune system alterations play an important role in the etiology of GD. In addition, genetics, epigenetics, and gut microbiota are also vital to the pathogenesis of GD.
Currently, three treatments are widely applied in clinical practice: radioiodine (RAI) therapy, antithyroid drugs (ATDs), and thyroidectomy. However, each method has side effects (Hu et al., 2020). Hence, the development of effective and specific treatments that can etiologically treat the disease is essential. Exploring the pathogenesis of GD may help to identify drug targets for novel therapies.
This article reviews the pathogenesis of GD from the aspects of genetics, epigenetics, cellular immunology, and gut microbiota. The aim is to provide a comprehensive overview of the basic principles to inform further development and clinical study of new drugs targeting the pathogenesis of GD.
Cellular Immunological Mechanisms of GD
Th1/Th2 Balance in GD
Stimulation by interleukin-12 (IL-12), interferon-γ (IFN-γ), and IL-2 and expression of the transcription factor T-bet induce naïve CD4+ T cells to differentiate into Th1 cells. This subset mainly synthesizes IL-1, IL-2, IFN-γ, and transforming growth factor beta (TGF-β), mainly through interactions with macrophages and other T lymphocytes. The presence of IL-4 suppresses the differentiation of naïve CD4+ T cells into Th1 cells and promotes the production of Th2 lymphocytes. GATA3 transcription factors are also implicated in Th2 development; this subset mainly synthesizes IL-4, IL-5, IL-6, IL-10, and IL-13. Th2 cells mainly interact with B and plasma cells, resulting in an increased production of antibodies and thereby mediating the humoral immune response (Eshaghkhani et al., 2016; Ramos-Levi and Marazuela, 2016).
Compared with controls, T-bet and IFN-γ mRNA levels in peripheral blood mononuclear cells (PBMCs) were prominently upregulated in GD patients, whereas GATA3 and IL-4 mRNA expression levels were downregulated. In addition, a significant decrease in plasma IL-4 levels was observed in GD patients, while IFN-γ levels were higher in patients than in controls. These results suggest a Th1/Th2 imbalance in GD, which may be related to the pathogenesis of the disease (Eshaghkhani et al., 2016). Another study showed that the serum levels of IL-2 and IL-10 were elevated, while IFN-γ levels were lower in patients than in controls. IL-4 levels were not significantly different between patients and controls (Esfahanian et al., 2013). However, the levels of IFN-γ and IL-4 are discordant with the results of Eshaghkhani et al., and may be attributable to different sample types. In mouse models, IL-10 and IFN-γ expression levels were significantly increased (Ueki et al., 2011).
Taken together, the abnormal expression of Th1/Th2 cells and the abnormal secretion of related cytokines disrupts of Th1/Th2 balance. However, their expression in GD remains controversial, and additional studies on the role of Th1/Th2 balance in GD are warranted.
Th17/Treg Balance in GD
Th17 cells are a subset of CD4+ T cells, primarily characterized by the generation of IL-17A, IL-17F, IL-21, and IL-22. The interaction of IL-1β, IL-6, and IL-23, the expression of the transcription factor retinoid-related orphan receptor gamma t (RORγt), and the activation of the STAT3 intracellular pathway play vital roles in the differentiation of Th17 cells (Tan et al., 2019). Treg lymphocytes are a heterogeneous population of lymphocytes characterized by their immunosuppressive function (Shao et al., 2018). Tregs act as negative regulatory cells through the synthesis of IL-10 and TGF-β (Tan et al., 2019). In recent years, the close relationship between Tregs and Th17 cells in GD has been explored in several studies (Li et al., 2016; Qin et al., 2017).
In Th17 cells and Tregs in PBMCs isolated from the peripheral blood of GD patients, the percentage of CD4+IL17+ T cells was significantly increased in GD patients compared with controls. In addition, the percentage of CD4+FoxP3+ Tregs was significantly reduced in GD patients in comparison to controls (Zhang et al., 2019). Analysis of Th17/Tregs ratio in the peripheral blood from GD patients indicated markedly lower ratios of CD4+IL17+/CD4+CD25+CD127− and CD4+IL17+/CD4+CD25+CD127−FoxP3+ than in the controls. In untreated GD patients, a statistically significant positive correlation in CD4+IL17+/CD4+CD25+CD127−, CD4+IL17+/CD4+CD25+CD127−FoxP3+ T cells, and TRAb levels was observed. A positive association was also reported between the proportion of CD4+IL-17+ T cells and thyroid-stimulating antibody (TSAb) levels (Bossowski et al., 2016). In a mouse model of GD, the number of CD4+CD25+Foxp3+Treg cells decreased and the expression of IL-17 mRNA increased, but the CD4+IL-17+T cell subpopulation did not significantly differ from the controls (Yuan et al., 2017). Treg cells and the balance of Th17/Treg might be associated with the pathogenesis of GD, yet Th17 cells have not been shown to promote GD (Zhou et al., 2012). IL-17 and TGF-β were all upregulated in GD patients and mouse models (Yuan et al., 2017; Han et al., 2020).
In summary, inhibition of Treg cells disrupts the Th17/Treg balance, which ultimately leads to the occurrence and development of GD. In addition, Th17/Treg balance is vital to the pathogenesis of other autoimmune diseases (Noack and Miossec, 2014).
Tfh Cells in GD
Follicular helper T (Tfh) cells have been identified as a new subgroup of effector helper T cells, which are vital to regulating the development of antigen-specific B cell immunity. Tfh cell differentiation is strictly regulated by the specific transcription factor Bcl-6, which is characteristically expressed in Tfh cells (Chen J et al., 2015).
An increasing body of evidence indicates that Tfh cells are significantly decreased in GD, which suggests that Tfh cells play an essential role in the pathogenesis of GD (Chen J et al., 2015; Wang et al., 2019). The proportions of effectors circulating Tfh (cTfh) and cTfh cell subsets (Tfh1, Tfh2, and Tfh17) in the peripheral blood of GD patients have been investigated. The proportion of effector cTfh cells and the Tfh2 subset were found to be increased in GD patients. Furthermore, a positive correlation between circulating Tfh2 (or PD-1+ Tfh) cells and serum TPOAb levels in GD patients has been reported (Liu et al., 2018). Tfhs and relevant factors in GD thyroid tissues were upregulated and protein expression levels of Tfh-related factors were also higher in GD thyroid tissues than in normal tissues (Zhang et al., 2015). These results imply that Tfh participates in the pathogenesis of GD.
These studies have revealed the abnormal expression of Tfh cells and related cytokines in GD patients. It shows that Tfh cells play a crucial role in the pathogenic mechanisms of GD and are thus potential therapeutic targets for GD. However, identification of the exact roles of Tfh cells in the pathogenesis of GD requires further investigation.
Breg Cells in GD
Breg cells are immunosuppressive cells that support immune tolerance. Through the production of IL-10, IL-35, and TGF-β, Breg cells inhibit immunopathology and prohibit the expansion of pathogenic T cells and other pro-inflammatory lymphocytes (Rosser and Mauri, 2015).
The proportion of IL-10-producing B (B10) cells in the peripheral blood of GD patients was decreased and negatively correlated with TRAb levels. In addition, mouse disease models have generated similar results (Ji et al., 2020). These findings indicate that B10 cells suppress autoantibody production, and their abundance is decreased in GD patients. With the onset of GD, the proportion of CD19+CD1dhiCD5+ Breg cells in the PBMCs of Chinese patients is reduced (Qin et al., 2017). In addition, a conspicuous decrease in the number of CD19+CD24hiCD27+IL-10+ and CD19+IL-10+ B lymphocytes has been observed in untreated GD patients compared with controls (Rydzewska et al., 2018b).
In conclusion, Th1/Th2 balance, Th17/Treg balance, Tfh cells, and Breg cells contribute to the occurrence and development of GD. However, their precise roles in disease pathogenesis remain unclear. Additional studies on the cellular immunological mechanisms underlying GD to elucidate their interactions and provide new therapeutic strategies for GD are necessary.
The Genetic Pathogenesis of GD
The Association Between HLA-Related Gene Polymorphisms and GD Susceptibility
HLA is located on chromosome 6p21.3, has a total length about 4,000 kb, and contains a set of closely related genes, which are expressed on the cell surface as HLA-I and HLA-II molecules. These genes play important roles in antigen presentation to CD4+ T and CD8+ T cells (Xiaoheng et al., 2017). HLA-A, -B, -C, and -DRA are overexpressed in GD patients, with HLA-C showing the most significant upregulated expression (Yin et al., 2014).
HLA-I genes have been identified to be significantly associated with GD. In Iran, people with HLA-A*68 (case/control: 80/180, p = 0.004, OR = 4.23) and -B*08 (case/control: 80/180, p = 0.030, OR = 3.72) genes were susceptible to GD, while HLA-A*33 (case/control: 80/180, p = 0.011) appeared to play a protective role (Mehraji et al., 2017). A meta-analysis on the correlation between GD and HLA-B*46 in Asian populations revealed that HLA-B*46 was associated with an increased risk of GD (Li et al., 2013). In addition, the coexistence of HLA-B*46 (case/control: 73/159, p < 2.4 × 10−8, OR = 5.455) and -Cw*01 (case/control: 73/159, p < 0.00016, OR = 2.208) may be a genetic marker of early-onset AITD in Koreans. Moreover, the HLA-A*02 (case/control: 41/159, p < 0.014, OR = 1.905) was shown to confer susceptibility to GD, whereas -Cw*07 (case/control: 41/159, p < 0.001, OR = 0.144) played a protective role (Cho et al., 2011).
Among the HLA-II genes, the allelic frequency of HLA-DRB1*07 (case/control: 41/159, p < 0.015, OR = 0.128), -DQB1*0201 (case/control: 80/180, p = 0.040, OR = 0.50) and -DQA1*0201 (case/control: 80/180, p = 0.045, OR = 0.37) was found to be significantly lower among GD patients, whereas the frequency of HLA-DRB1*08 (case/control: 41/159, p < 5.6 × 10−5, OR = 3.436) was higher in GD patients than in controls (Cho et al., 2011; Mehraji et al., 2017). In a Romanian population, HLA-DRB1*03 (case/control: 77/445, p < 0.0001, OR = 3.29) and -DRB1*11 (case/control: 77/445, p = 0.045, OR = 1.70) were identified as the potential main susceptibility alleles of HLA-DRB1 related to GD. In addition, GD patients with HLA-DRB1*03/11 had a higher FT4/TT3 ratio and anti-TgAb levels, but HLA-DRB1*01 (case/control: 77/445, p = 0.007, OR = 0.20) and -DRB1*15 (case/control: 77/445, p = 0.038, OR = 0.42) were likely to have a protective function (Martin et al., 2014). Some HLA alleles were also associated with clinical features in GD patients, HLA-DRB1*1301 and -DQB1*0603 carriers more frequently develop larger goiters (Park et al., 2005).
In conclusion, HLA is the most diverse gene in the human genome, and its variation may be associated with the occurrence of GD or play a protective role. Hence, individuals carrying certain HLA alleles are susceptible to GD.
The Association Between Non-HLA Gene Polymorphisms and GD Susceptibility
Association of T Cell-Related Gene Polymorphisms with GD Susceptibility
CTLA4 Gene Polymorphisms in GD
The CTLA4 gene is located on human chromosome 2q33 and has a total length of approximately 8.44 kb. It contains four exons and encodes 233 amino acids, generating two CTLA4 protein isoforms, namely, full-length CTLA4 (flCTLA4) and soluble CTLA4 (sCTLA4) (Fang et al., 2015). FlCTLA4 is expressed on the surface of the cellular membrane of activated T cells to play an immunosuppressive role. SCTLA4 may interfere with the coordinated stimulation signal and inhibit T cell proliferation (Patel et al., 2016).
The correlation between CTLA4 and the incidence of GD has been confirmed in multiple populations (Xiaoheng et al., 2017). In the Chinese Han population, the synergistic interactions of the CTLA4 SNPs, including rs231775 (case/control: 260/248, p = 0.002, OR = 1.521), rs231779 and rs3087243 (case/control: 260/248, p = 0.007, OR = 1.615), were associated with a significant increase in the risk of GD (Chen et al., 2018). Similarly, CTLA4/rs3087243 is also significantly associated with GD susceptibility in Kashmiri (case/control: 135/150, p < 0.001, OR = 2.21) as well as Brazilian (case/control: 282/308, p < 0.0001, OR = 2.593) populations (Shehjar et al., 2020; Bufalo et al., 2021). In addition, a significant association between CTLA4/rs231775 (case/control: 135/150, p < 0.001, OR = 1.85) and GD development has been reported in a Kashmiri population (Shehjar et al., 2020). Therefore, CTLA4 polymorphisms are correlated to GD susceptibility. In patients with ophthalmopathy, the frequency of the rs231775 G allele was higher than in those without ophthalmopathy, suggesting that GD patients with CTLA4/rs231775 were more likely to develop ocular symptoms (Esteghamati et al., 2009).
PTPN22 Gene Polymorphisms in GD
PTPN22 is located on chromosome 1p13.3-13.1, consists of 16 exons, and encodes the 110-kDa lymphoid protein tyrosine phosphatase (LYP). LYP plays a negative regulatory role in T cell signaling (Xiaoheng et al., 2017).
With respect to the association between PTPN22 and GD susceptibility, the most studied SNP is rs2476601. The SNP rs2476601 is located in the N-terminal proline-rich motif and results in the substitution of arginine (Arg) with tryptophan (Trp) at codon 620 (Zhebrun et al., 2011). A significant association between PTPN22/rs2476601 (case/control: 171/200, p < 0.05, OR = 4.23) polymorphism and GD has been identified in a population from Northwest Russia and in Caucasians (Zhebrun et al., 2011; Luo et al., 2012). Moreover, female teenagers with GD in Poland have also exhibited a similar correlation, and allele A (case/control: 142/160, p = 0.009, OR = 2.13) was found to be an important risk factor (Rydzewska et al., 2018a). However, Bufalo et al. reported the opposite findings in Brazil (Bufalo et al., 2021). On the other hand, C allele of PTPN22/rs3789604 may be associated with an increased risk of liver damage in GD patients (Zhang et al., 2017).
FoxP3 Gene Polymorphisms in GD
The FoxP3 gene, which is located on the X chromosome, consists of 11 exons and encodes 431 amino acids. FoxP3 plays an essential role in the development and function of T cells, as well as differentiation into natural Tregs (Shehjar et al., 2018).
In a Caucasian population, the FoxP3 rs3761549 polymorphism was suggested to be involved in GD development in women (case/control: 109/75, p = 0.03, OR = 2.174) (Bossowski et al., 2014). In a Kashmiri population, the AC genotype (case/control: 135/150, p < 0.001, OR = 3.48) of rs3761548 and the TT genotype (case/control: 135/150, p < 0.001, OR = 5.62) of rs3761549 were considered risk factors of GD development (Shehjar et al., 2018). Similarly, in the Chinese Han population, allele A and genotype AA of rs3761548 were correlated with susceptibility to GD (Zheng et al., 2015).
CD25 Gene Polymorphisms in GD
CD25, also known as IL-2RA, is located on chromosome 10p15.1 and has a length of 60 kb. The expression of CD25 is a pivotal feature and phenotypic marker of Tregs (Brand et al., 2007).
It has been reported that CD25/rs2104286 (case/control: 650/1300; AA genotype: p = 8.772 × 10−6, OR = 1.636; A allele: p = 0.004, OR = 1.322) confers to GD susceptibility in the Chinese Han population (Du et al., 2021). The minor allele A (case/control: 1474/1609, p = 0.00017, OR = 1.43) of rs41295061 and the homozygous AA genotype (case/control: 1474/1609, p = 0.0053, OR = 1.54) of rs11594656 conferred susceptibility to GD in a Russian population (Chistiakov et al., 2011).
In summary, the expression in T cell-related genes may affect the differentiation and function of T cell subsets, and polymorphisms may lead to an imbalance of T cell subsets, thereby contributing to the occurrence of GD.
The Association Between B Cell Related Gene Polymorphisms and GD Susceptibility
IKZF3 Gene Polymorphisms in GD
The IKZF3 gene is found on chromosome 17q12-q21.1 in humans, contains nine exons, and is 104 kb in length. IKZF3 is a crucial transcription factor that inhibits the proliferation and differentiation of B cells (Li et al., 2018).
Li et al. were the first to report the relationship between IKZF3 polymorphisms and GD. They found that the association between the minor alleles of rs2941522 (case/control: 604/814, p = 0.02, OR = 1.21), rs907091 (case/control: 604/814, p = 0.006, OR = 1.25), rs1453559 (case/control: 604/814, p = 0.007, OR = 1.25), rs12150079 (case/control: 604/814, p = 0.006, OR = 1.29), and rs2872507 (case/control: 604/814, p = 0.004, OR = 1.27) and GD was significant (Li et al., 2018). More functional studies are needed to determine the contribution of IKZF3 polymorphisms to GD pathogenesis.
BAFF Gene Polymorphisms in GD
BAFF is located on chromosome 13q33, and plays an essential role in regulating the maturation, proliferation, differentiation and survival of B cells (Lin et al., 2016a).
Lin et al. found that serum BAFF levels in Chinese GD patients were higher than those of healthy controls. Further research showed that serum BAFF levels were significantly correlated with TSHRAb and anti-TPOAb levels only in women with active GD (Lin et al., 2016a). Lin et al. also studied SNPs in BAFF and showed that the G allele (case/control: 223/243, p = 0.009, OR = 0.70) of rs2893321 may decrease the risk of GD in women (Lin et al., 2016b). Lane et al. found that the frequency of another SNP (rs4000607) (case/control: 444/447, p = 0.019, OR = 1.80) was significantly different between GD patients and controls (Lane et al., 2019).
CD40 Gene Polymorphisms in GD
CD40 is located on chromosome 20q12-q13 and encodes a co-stimulatory receptor protein that plays an important role in B lymphocyte differentiation and antibody secretion. It is a member of the tumor necrosis factor receptor superfamily that has been associated with the pathogenesis of multiple autoimmune diseases (Wang D. et al., 2017).
With respect to the association between CD40 and GD susceptibility, rs1883832 is the most widely studied SNP. The CC genotype and the C allele with rs1883832 were associated with GD in the Japanese (case/control: 61/42; CC genotype: p = 0.041, OR = 2.438; C allele: p = 0.031, OR = 1.972) and Chinese populations (case/control: 196/122; CC genotype: p = 0.003, OR = 2.043; C allele: p = 0.008, OR = 1.57) (Inoue et al., 2012; Wang D. et al., 2017). Nevertheless, there was no significant association between CD40/rs1883832 and GD in the Brazilian and Pakistani populations (Mustafa et al., 2018; Bufalo et al., 2021). In addition, the G allele (case/control: 196/122, p = 0.001, OR = 5.472) of C64610G have been described as susceptibility factors for GD (Wang D. et al., 2017).
BACH2 Gene Polymorphisms in GD
BACH2 is located on chromosome 6q15, and encodes the BACH2 protein, which regulates the transformation of B cells into plasma cells. BACH2 is expressed in all PBMC isoforms, with the highest expression in B cells. A genome-wide association study (GWAS) via three stages verified that BACH2 is associated with GD in the Chinese Han population. BACH2/rs2474619 (case/control: 8882/9431, p = 3.28 × 10−8, OR = 1.13) is significantly associated with GD. However, there was no significant association between the rs2474619 genotypes and the level of BACH2 gene expression (Liu et al., 2014).
FAM167A-BLK Gene Polymorphisms in GD
FAM167A-BLK is located on chromosome 8p23.1. The function of the FAM167A protein, also known as c8orf13, is unclear. BLK encodes a Src kinase that functions as a B cell signal transducer, is mainly expressed in B cells, and may influence the proliferation and differentiation of B cells. Song et al. found that the GG genotype of rs2618431 (case/control: 624/797, p = 0.04, OR = 1.246) may contribute to susceptibility to GD (Song et al., 2018).
In summary, B lymphocytes are the main effectors of antigen presentation and autoantibody production. B cell-related gene polymorphisms affect B cells, which are involved in the occurrence of GD.
Association Between Thyroid Hormone Related Gene Polymorphisms and GD Susceptibility
TSHR Gene Polymorphisms in GD
TSHR is located on chromosome 14q31 and encodes TSHR (Rydzewska et al., 2018a), a 764 amino acid polypeptide synthesized by a G protein-coupled receptor. It undergoes post-translational cleavage, producing TSHR A and B chains. The extracellular A subunit, which is shed, leads to the production of auto-antigens, facilitating the activation of non-self-tolerant CD4+ T cells. This ultimately results in the generation of stimulating antibodies (Hesarghatta Shyamasunder and Abraham, 2017).
The frequencies of the GG genotype (case/control: 180/111, p = 0.000105, OR = 2.602) and the G allele (case/control: 180/111, p = 0.0008, OR = 1.853) of rs4411444, the AA genotype (case/control: 180/111, p = 0.0228, OR = 1.740) of rs2300519, and the AA genotype (case/control: 180/111, p = 0.0052, OR = 1.984) and the A allele (case/control: 180/111, p = 0.0158, OR = 1.548) of rs179247 were higher in GD patients than in controls. The frequencies of these genotypes and alleles, as well as the rs2300519 A allele (case/control: 62/48, p = 0.0243, OR = 1.967) and the rs4903961 GG genotype (case/control: 62/48, p = 0.0147, OR = 2.588) and G allele (case/control: 62/48, p = 0.0166, OR = 2.061), were higher in patients with intractable GD than in controls and patients in GD remission (Fujii et al., 2017; Rydzewska et al., 2018a). In addition, rs179247 allele A was found significantly more frequently in GD patients without oculopathy than in patients with oculopathy, indicating that this allele is associated with a lower risk of GD in patients with oculopathy (Jurecka-Lubieniecka et al., 2014).
TG Gene Polymorphisms in GD
TG, which is located on chromosome 8q24, encodes a 660-kDa glycoprotein that supports the generation of thyroid hormones (Xuan et al., 2019). A significant relationship between rs2069550 (case/control: 436/316, p = 0.01, OR = 1.49) of the TG gene and GD has been reported (Gu et al., 2010).
In the Chinese Han population, rs2294025 (case/control: 9757/10626, p = 1.52 × 10−9, OR = 1.16) and rs7005834 (case/control: 9757/1062, p = 1.62 × 10−7, OR = 1.16) are two independent loci associated with susceptibility to GD (Xuan et al., 2019). Also, the frequency of the TT genotype (case/control: 131/89, p = 0.0283, OR = 0.484) of rs2703013 was significantly lower in GD patients than in controls. The distribution of the rs2958692 T allele (case/control: 50/40, p = 0.0055, OR = 2.382) was significantly different between patients with intractable GD and those with GD in remission (Mizuma et al., 2017).
In summary, thyroid hormone-related gene polymorphisms lead to the abnormal expression of related proteins. In addition, the production of autoantibodies leads to thyroid hormone secretion disorders, which are involved in the pathogenesis of GD.
The Association Between Apoptosis Related Gene Polymorphisms and GD Susceptibility
Apoptosis as a regulatory mechanism can remove excess unnecessary cells such as autoreactive lymphocytes. A lack of this ability causes the proliferation of cells that react with their own antigen, in turn leading to autoimmune diseases, including GD (Bossowski et al., 2008).
VDR Gene Polymorphisms in GD
VDR is located at chromosome 12q12-12q14 and consists of eight coding exons and 3 alternative 5′-noncoding exons spanning more than 75 kb of DNA (Meng et al., 2015). VDR is a nuclear receptor that controls the transcription of regulatory genes that are related to calcium metabolism and immune responses (Inoue et al., 2014).
Allele A (case/control: 417/301, p = 0.041, OR = 1.278) of rs7975232 has been associated with GD in the Chinese Han population (Meng et al., 2015). Similarly, the frequency of the C allele of VDR/rs7975232 was higher in patients with AITD (case/control: 255/76, p = 0.037, OR = 1.514), particularly GD patients (case/control: 139/76, p = 0.0349, OR = 1.594), than in controls (Inoue et al., 2014). Moreover, the A allele (case/control: 650/1209, p = 2.62 × 10−2, OR = 1.20) and AA genotype (case/control: 650/1209, p = 3.45 × 10−4, OR = 1.87) of rs7975232 were associated with a significant increase in susceptibility to GD. In GD patients with ophthalmopathy, the frequency of the rs7975232 AA genotype was higher than in those without ophthalmopathy, suggesting that GD patients with VDR/rs7975232 are more likely to develop ocular symptoms (Zhou et al., 2021). The frequency of the TT genotype of rs731236 was associated with the higher expression of VDR (Inoue et al., 2014). A meta-analysis found that the TT genotype (case/control: 2380/2235, p = 0.025, OR = 1.42) of rs731236 was associated with an increased risk of GD, but rs1544410 and rs7975232 of VDR were not different between GD patients and controls (Veneti et al., 2019).
FAS Gene Polymorphisms in GD
The FAS gene is located on chromosome 10q23. The receptor-ligand combination FAS-FAS ligand (FASL) is a chief mediator of apoptosis in many main physiologic processes. FAS and FASL both exist as membrane-bound and soluble types. The level of soluble FAS in subjects with the AA genotype of FAS -670 A/G was markedly higher than in those carrying the GG genotype (Mahfoudh et al., 2007).
There were no significant differences in rs2234767, rs1800682, and rs763110 between controls and AITD patients. However, serum soluble FASL levels in GD patients were significantly higher than in controls (Inoue et al., 2016).
Bcl-2 Gene Polymorphisms in GD
Bcl-2 is an anti-apoptosis factor, and the upregulated expression of Bcl-2 facilitates cell survival. In the Bcl-2/rs1800477 polymorphism, the frequency of the G allele (case/control: 264/79, p = 0.011, OR = 2.630), which is related to stronger anti-apoptotic function than the A allele, was markedly higher in AITD patients than in controls. No association between the polymorphism of the rs2279115 site and GD susceptibility has been observed (Inoue et al., 2016).
TNFR1 and TNFR2 Gene Polymorphisms in GD
The TNFR1 and TNFR2 genes, which encode TNF-RI and TNF-RII have been mapped to chromosomes 12p13 and 1p36, respectively (Glossop et al., 2005). TNF-RI possesses a death domain which can transduce the signal for cell death. While no differences in the genotype and allele frequency of TNFR1 rs2234649 have been observed, the G allele (case/control: 160/87, p = 0.038, OR = 1.827) of TNFR2 rs1061622 is associated with an increased risk of GD (Inoue et al., 2016).
RNASET2 Gene Polymorphisms in GD
RNASET2 is located on chromosome 6q27 and is the only RNase T2 family member in humans. RNASET2 plays a role in the initiation of human dendritic cells for the Th2 polarization of CD4+ T cells (Chen X J et al., 2015). There was a significant association between the risk allele G (case/control: 701/938, p = 0.005, OR = 1.225) of rs9355610 and a decrease in RNASET2 mRNA levels in a Chinese population (Chen et al., 2015). This association was also identified in the previous GWAS of GD (case/control: 5530/5026, p = 6.85 × 10−10, OR = 1.19) (Chu et al., 2011).
Genetic factors also play an important role in other autoimmune diseases. HLA-DRB1/rs13192471 confers a genetic predisposition to RA in a northeast Indian population (Das et al., 2018). FoxP3, VDR, PTPN22, CTLA4, and CD40 gene polymorphisms were associated with susceptibility to autoimmune diseases including RA, systemic lupus erythematosus, alopecia areata, and psoriasis (Wu et al., 2016; Bin Huraib et al., 2018; Hashemi et al., 2018; Moravvej et al., 2018; Fouad et al., 2019). With the extensive development of GWAS, many GD susceptibility genes have been identified that partly explain the clinical manifestations and pathogenesis of GD (Supplementary Table S1). However, due to differences in region, ethnicity, sample size, and genotyping methods, studies on the correlation between gene polymorphisms and GD susceptibility have generated inconsistent results. Therefore, it is necessary to increase sample sizes to elucidate the relationships between gene polymorphisms and GD susceptibility which can also provide reference for the identification of the susceptibility genes of other autoimmune diseases.
Epigenetic Pathogenesis of GD
The Role of DNA Methylation in GD
DNA methylation is a typical epigenetic modification (Wang B. et al., 2017). The methylation reaction is catalyzed by a cluster of enzymes, including DNA methyltransferases (DNMTs) (Picascia et al., 2015). The DNMT family includes DNMT1, DNMT2, DNMT3A, DNMT3B, and DNMT3L (Nawrocki et al., 2017). Among these, three enzymes are known to have DNMT activity: DNMT3a and DNMT3b are responsible for de novo methylation, while DNMT1 is required for maintaining DNA methylation (Jeltsch and Jurkowska., 2014). Compared to controls, the relative mRNA level of DNMT1 was significantly reduced in GD patients (Cai et al., 2015). Global hypomethylation and downregulation of DNMT1 mRNA expression in CD3+ T and CD19+ B lymphocytes in patients with newly diagnosed GD have also been observed, while there were no significant differences in DNMT3a and DNMT3b expression between GD patients and controls (Guo et al., 2018).
A genome-wide analysis of DNA methylation identified that 365 and 3322 CpG sites in CD4+ and CD8+ T cells are respectively differentially methylated in GD patients in comparison with controls. In addition, hypermethylation in the first intron of the TSHR gene has been associated with GD (Limbach et al., 2016). However, these findings are discordant with the results of Guo et al. on global hypomethylation patterns in CD3+ T and CD19+ B lymphocytes with GD patients, and may be attributable to different types of cells. Another study found a negative association between global methylation in B cells and serum TPOAb (Guo et al., 2018).
A statistically significant difference in methylation status intercellular adhesion molecule 1 (ICAM1) between GD patients and controls was observed, and there was a significant association between decreased ICAM1 methylation and exophthalmos in GD patients (Shalaby et al., 2019). Another study showed that DNA methylation decreased and DNA hydroxymethylation increased at the promoter of ICAM1, and this change was correlated with ICAM1 mRNA expression (Liu et al., 2017). Morita et al. assessed the methylation status of six CpG sites in the TNFA promoter region and found that the methylation level of the -72 CpG was notably higher in GD patients than in controls (Morita et al., 2018). Kyrgios et al. reported that IL-2RA gene promoter methylation status significantly differed among GD patients and controls, and the level of IL-2RA promoter DNA methylation was negatively correlated to anti-TPO and anti-TSI (Kyrgios et al., 2020).
Based on the above information, it can be hypothesized that abnormal DNA methylation is involved in the pathogenesis of GD. Nevertheless, it is still necessary to explore the abnormal DNA methylation of specific genes to provide new targets for the treatment of GD.
The Role of Histone Modifications in GD
Chromatin is a dynamic structure that helps encapsulate the entire eukaryotic genome into the nucleus, as well as regulates DNA-related metabolic processes, including DNA transcription, recombination, repair, and replication. The basic unit of chromatin is composed of about 146 base pairs of DNA wrapped around a histone octamer (Li et al., 2007; Venkatesh and Workman, 2015; Martire and Banaszynski, 2020). Compared with GD patients carrying no chromosomal abnormalities, those with Turner- or Down-syndrome exhibited a prominent higher rate of transition from the previous Hashimoto’s hypothyroidism to Graves’ hyperthyroidism. This implies that chromosomal abnormalities and the coexistence of Hashimoto’s thyroiditis may increase the risk of developing GD. In addition, the dynamics of chromatin structure are strictly regulated by a variety of mechanisms including histone modification (Aversa et al., 2014). Histone modifications include acetylation, phosphorylation, methylation, ubiquitylation and SUMOylation (Bannister and Kouzarides, 2011). Mounting evidence has shown that abnormal histone modification profiles contribute to an imbalance in immune response, which in turn leads to the development of various autoimmune diseases (Araki and Mimura, 2017). In general, high levels of acetylation and trimethylation of the Lys-4 residue of histone 3 (H3K4me3) were detected in the promoter regions of transcriptionally active genes in individuals with GD (Stefan et al., 2011).
Genome-wide analysis of DNA methylation in CD4+ and CD8+ T cells revealed that hypermethylation was related to reduced levels of H3K4me3 and H3K27ac marks in several T cell signaling genes. However, these results were discordant to the findings of Guo et al. on global hypomethylation in CD3+ T and CD19+ B lymphocytes of GD patients, which may be attributed to the use of different types of cells. Regions with lower H3K4me3 or H3K27ac signals were related to hypermethylated CpGs and areas with strong H3K4me3 or H3K27ac signals overlapping with hypomethylated CpGs. In addition, gene ontology and pathway analysis have suggested reduced H3K4me3 and H3K27ac signals in genes involved in T cell activation (Limbach et al., 2016). In PBMCs, global histone H3K9 methylation was significantly downregulated in GD patients compared with controls. Global H3K4 methylation was decreased in GD patients, but this difference was not statistically significant (Yan et al., 2019).
Yan et al. reported that the mRNA expression of histone deacetylases HDAC1 and HDAC2 in the PBMCs of GD patients was significantly higher than that of controls. This indicates that histone acetylation modifications are abnormal in PBMCs of GD patients (Yan et al., 2015).
In summary, histone expression is associated with methylation status and histone modifications of the genes that may regulate its expression. It is obvious that histone modifications influence the occurrence of GD. However, studies on the role of histone modification in GD are limited, and more investigations are needed to elucidate its role in GD.
The Role of miRNAs in GD
MicroRNAs are a class of small non-coding RNAs that negatively control gene expression by directly combining with the 3′ untranslated region of their mRNA targets (Mehta and Baltimore, 2016). MicroRNAs play a crucial regulatory role in Toll-like receptor signaling, resulting in the activation of the NF-κB, IRF, and AP-1 transcription factors, which regulate the expression of pro-inflammatory cytokines (Wu et al., 2015).
The effects of miRNAs in GD are described below. MiR-181d was found to be upregulated in GD patients and miR-346 was significantly downregulated. Furthermore, miR-346 regulates CD4+CXCR5+ T cells by acting on Bcl-6, a positive regulator of Tfh cells (Chen J et al., 2015). The expression level of miR-146a in plasma was significantly reduced in AITD patients in comparison with controls. However, the expression level of miR-146a in PBMCs was increased. The reason of discordant expression levels in miR-146a in plasma and PBMCs might be that miR-146a inhibits the release of miR-146a from PBMC into plasma, and thus inducing the excessive inflammation in the thyroid of AITD patients. Thus miR-146a may be related to the occurrence of AITD (Otsu et al., 2017). MiR-155 enhances Treg and Th17 cells differentiation by targeting SOCS1 (Yao et al., 2012). However, another study showed that overexpression of miR-155 results in increased Th1 differentiation. This result is complementary to the observed bias toward Th2 differentiation in CD4+ T cells lacking miR-155, thereby contributing to the regulation of Th1/Th2 balance (Banerjee et al., 2010). The expression of miR-154, miR-376b, and miR-431 was suppressed in the PBMCs of initial GD patients (Liu et al., 2012). In comparison to GD patients in remission, circulating miR-23b-5p and miR-92a-39 levels were reduced in intractable GD patients, whereas let-7g-3p and miR-339-5p were increased (Hiratsuka et al., 2016). The serum levels of miR-16, miR-22, miR-375, and miR-451 were increased in GD patients compared with controls (Yamada et al., 2014). MiR-23a-3p levels were significantly downregulated in GD patients. Animal experiments have demonstrated that miR-23a-3p overexpression strengthens Tregs function in vivo (Zhang et al., 2019). Yao et al. discovered five miRNAs that were differentially expressed between GD patients and controls. There were four upregulated miRNAs (has-miR-122-5p, 16-1-3p, 221-3p, and 762) and one downregulated miRNA (has-miR-144-3p) (Yao et al., 2019). Qin et al. reported that miR-22 and miR-183 were overexpressed in GD patients, but miR-101, miR-197, and miR-660 levels were downregulated in the thyroid tissue of GD patients (Qin et al., 2015).
Martínez-Hernández et al. indicated that miRNAs are associated with autoimmune antibody levels in GD. They found that miR-let7d-5p was negatively correlated to TPOAb levels, while miR-21-5p and miR-96-5p were positively correlated with TPOAb, TgAb, and TRAb levels. MiR-142-3p and miR-301a-3p were only positively associated with TRAb levels, and miR-6503-3p was correlated with TPOAb levels (Martinez-Hernandez et al., 2018).
The above studies indicated that miRNAs participate in the regulation of gene and autoimmune antibody expression, ultimately contributing to the occurrence and development of GD and may hence serve as biomarkers for the diagnosis, treatment, and prognosis of GD.
The Role of lncRNAs in GD
Long non-coding RNAs (LncRNAs) are non-coding RNAs with lengths of >200 nt that are transcribed by RNA polymerase II. Most lncRNAs are transcribed from the 5′ end to 3′ end of genes, but lncRNAs do not encode proteins (Ponting et al., 2009).
LncRNAs have a wide range of effects and are involved in transcriptional regulation, post-transcriptional regulation, and participate in chromatin remodeling and protein transport (Quinn and Chang, 2016). LncRNAs can be used as miRNA sponges that control the expression of target genes in various physiopathological processes (Zhang et al., 2018). Qiu et al. found that lncRNA-MEG3 may serve as a competing endogenous RNA that regulates RORγt expression via miR-17, thus affecting Th17/Treg balance. This indicates that the lncRNA/miRNA axis plays an important role in the etiology of GD (Qiu et al., 2019). The level of the RNA transcript Heg in monocytes is negatively correlated to TRAb levels in untreated GD patients (Christensen et al., 2008). However, Christensen et al. found that a decrease in TRAb levels in treated GD patients was not correlated to changes in Heg RNA levels (Christensen et al., 2011). These results suggest that Heg is related to the pathogenesis of GD.
Studies on the role of epigenetic factors in the development of GD have improved our understanding of the pathogenic mechanisms of GD (Supplementary Table S2). In parallel, epigenetics also play a pivotal role in other autoimmune diseases (Long et al., 2016; Karagianni and Tzioufas, 2019). Understanding the role of epigenetics in GD allows the generation of new ideas on the pathogenesis of other autoimmune diseases. Therefore, further studies are needed to: 1) elucidate the relationship between epigenetic factors and GD, 2) explore the causes and mechanisms of epigenetic changes, and 3) identify effective molecular targets for the diagnosis, treatment and prognosis of GD.
The Gut Microbiota and GD
The human microbiome is an essential but largely ignored part of the human genome (Grice and Segre, 2012). The collection of bacteria, archaea, and eukarya colonizing the human gastrointestinal (GI) tract is called the “gut microbiota” which has coevolved with the host for thousands of years to form a complex and mutually beneficial relationship. The microbiota brings many benefits to the host, through a series of physiological functions. Microbial composition disruption is known as dysbiosis. With the development of increasingly sophisticated methods for analyzing and characterizing complex ecosystems, the role of microbiota in a variety of intestinal and extraintestinal diseases has become increasingly clear (Thursby and Juge, 2017).
An analysis of the relative richness and diversity has indicated that the diversity of intestinal bacteria in GD patients was lower than in controls. The microbiota in GD patients mainly included four phyla (Firmicutes, Bacteroidetes, Proteobacteria, and Actinobacteria). One study showed that in GD patients, the abundance of Bacteroidetes and Proteobacteria was higher and the level of Firmicutes was lower than in controls (Ishaq et al., 2018). Moreover, another study revealed that the relative abundance of Bacteroidetes and Actinobacteria increased and Firmicutes decreased in GD patients compared with healthy controls (Chang et al., 2021). In the Han population of Southwest China, GD patients showed a significantly higher abundance of Firmicutes, Proteobacteria, and Actinobacillus. Additionally, the proportion of Firmicutes was significantly higher in GD patients than in controls, while the proportion of Bacteroidetes was significantly lower in GD patients than in controls (Yang et al., 2019). However, another study on the Han population of Shanghai China found that the proportion of Firmicutes in GD patients was significantly reduced, while the proportion of Bacteroidetes was significantly increased (Jiang et al., 2021). The reason for these opposite results may be that two studies collected samples from different regions and had small sample sizes (case/control: the former 15/15; the latter 46/59).
In the genus level, the abundance of Lactobacillus is higher in GD in several studies, which suggests Lactobacillus may be important to the pathogenesis of GD. The abundance of the genera Oribacterium, Lactobacillus, Aggregatibacter, and Mogibacterium was significantly higher in GD patients than in controls (case/control: 15/15, p < 0.05) (Yang et al., 2019). Furthermore, GD patients have higher counts of Bacteroides and Lactobacillus (case/control: 45/59, p < 0.05) (Jiang et al., 2021). Chen et al. found that the relative abundance of Lactobacillus, Veillonella and Streptococcus were significantly higher in GD (case/control: 15/14, p < 0.05). Moreover, Synergistetes and Phascolarctobacterium were negatively correlated to TRAb, which suggests that Synergistetes and Phascolarctobacterium play a protective role in GD by protecting the thyroid gland (Chen et al., 2021). Changes in gut microbiota composition have also been reported in other autoimmune diseases and generally involve Proteobacteria, Actinobacteria, and Lactobacillus (Xu et al., 2019).
Some specific intestinal bacteria also play a crucial role in the development of GD. One study showed that the abundance of the genus Prevotella was significantly higher in GD patients (Ishaq et al., 2018). The incidence of AITD was significantly greater in patients with Helicobacter pylori infection (Arslan et al., 2015). Similarly Yersinia enterocolitica infection may play a role in the etiology of GD in Turkey (Corapcioglu et al., 2002).
The gut microbiota have functions to regulate Th17 cells (Duan and Kasper, 2011; Kim et al., 2016). Similarly, the presence of gut microbiota can increase the induction of Tregs, and the butyrate produced by gut microbiota enhances Tregs differentiation and maturation (Furusawa et al., 2013). Propionic acid is an important metabolite of Bacteroides fragilis (B. fragilis). Su et al. found that B. fragilis regulates Th17/Treg balance via the propionic acid-mediated pathway (Su et al., 2020). We can infer that a disruption of the gut microbiota result in abnormal secretion of Th17 and Tregs, which leads to the development of GD.
Gut microbiota disorders may induce the abnormal expression of Th17 and Treg cells, which in turn influence the occurrence and development of GD. Therefore, in-depth studies on gut microbiota disorders may provide new approaches for the treatment and prevention of this disease (Supplementary Table S3).
Future Perspectives and Treatment of GD
Recently, Chen et al. observed that Tregs dysregulation and the expansion of pathogenic T cell clones might be involved in the long-lasting phase of GD via upregulating chemotaxis or inflammation response, and proposed a treatment with ATDs combined therapies, especially those aimed at improving Tregs frequencies or targeting specific expanded pathogenic TCR clones (Chen et al., 2021).
ATDs combination therapy provides new ideas and methods for the treatment of GD. In particular, targeting the Th17/Tregs balance may be a new approach for treating GD in the future.
Conclusion
Overall, genetics, epigenetics, cellular immunology, and gut microbiota play crucial roles in the pathogenesis of GD (Figure 1). In parallel, these factors also influence the occurrence and development of other autoimmune diseases. Consequently, if we elucidated the pathogenesis of GD, then it will also provide new ideas in the study of the pathogenesis of other autoimmune diseases and facilitate in the development of methods for the treatment of other autoimmune diseases.
FIGURE 1. The interaction of genetics, epigenetics, cellular immunology, and gut microbiota involved in the development of GD.
GD is an organ-specific autoimmune disease, in which the T cell immune imbalance of patients not only participates in the occurrence and development of the disease, but also affects the function of the thyroid gland. However, no studies have explained why the changes of PBMC only affect the thyroid gland, and further studies should focus on this issue.
In this Review, we found that the similarities between GD patients and mouse models (Supplementary Figure S1). This discovery helps us to have a clearer comprehension of GD, thereby providing a direction for future research.
Although we have some understanding of the pathogenesis of GD through the Review, our comprehension of the etiology of GD is still limited. In the future, we need to further investigate the factors mentioned above in order to find potential targets for GD therapy.
Author Contributions
HY conceived the idea. FZ drafted the manuscript and searched the literatures. XW, GT, WW, GZ, XM, DT, LW, XS, and HY participated in the manuscript preparation and revised it critically for important intellectual content. Final approval of the version was submitted by all authors.
Funding
This work was supported by National Key R&D Program of China (grant number: 2018YFC1004300), National Natural Science Foundation Project of China (grant number: 82160154, 81670844), the Key Project of Guizhou Provincial Science and Technology Department (grant number: QKH-JC-2019-1464), the Excellent Talent Support Program of Guizhou Provincial Education Department (grant number: QJH-KY-2017-077), the Science and Technology Foundation of Guizhou Province (grant number: QKH-PTRC-2018-5772-042), the Science and Technology Project of Zunyi (grant number: ZSKH-HZ-2020-35), and the Program for Excellent Young Talents of Zunyi Medical University (grant number: 18-ZY-001).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.794912/full#supplementary-material
References
Araki, Y., and Mimura, T. (2017). The Histone Modification Code in the Pathogenesis of Autoimmune Diseases. Mediators Inflamm. 2017, 2608605. doi:10.1155/2017/2608605
Arslan, M. S., Ekiz, F., Deveci, M., Sahin, M., Topaloglu, O., Karbek, B., et al. (2015). The Relationship between Cytotoxin-Associated Gene A Positive Helicobacter pylori Infection and Autoimmune Thyroid Disease. Endocr. Res. 40, 211–214. doi:10.3109/07435800.2015.1015727
Aversa, T., Lombardo, F., Corrias, A., Salerno, M., De Luca, F., and Wasniewska, M. (2014). In Young Patients with Turner or Down Syndrome, Graves’ Disease Presentation Is Often Preceded by Hashimoto’s Thyroiditis. Thyroid 24, 744–747. doi:10.1089/thy.2013.0452
Banerjee, A., Schambach, F., DeJong, C. S., Hammond, S. M., and Reiner, S. L. (2010). Micro-RNA-155 Inhibits IFN-Gamma Signaling in CD4+ T Cells. Eur. J. Immunol. 40, 225–231. doi:10.1002/eji.200939381
Bannister, A. J., and Kouzarides, T. (2011). Regulation of Chromatin by Histone Modifications. Cell Res 21, 381–395. doi:10.1038/cr.2011.22
Ben-Skowronek, I., Szewczyk, L., Kulik-Rechberger, B., and Korobowicz, E. (2013). The Differences in T and B Cell Subsets in Thyroid of Children with Graves' Disease and Hashimoto's Thyroiditis. World J. Pediatr. 9, 245–250. doi:10.1007/s12519-013-0398-0
Bin Huraib, G., Al Harthi, F., Arfin, M., Rizvi, S., and Al-Asmari, A. (2018). The Protein Tyrosine Phosphatase Nonreceptor 22 (PTPN22) R620W Functional Polymorphism in Psoriasis. Clin. Med. Insights Arthritis Musculoskelet. Disord. 11, 1179544117751434. doi:10.1177/1179544117751434
Borst, J., Ahrends, T., Bąbała, N., Melief, C. J. M., and Kastenmüller, W. (2018). CD4+ T Cell Help in Cancer Immunology and Immunotherapy. Nat. Rev. Immunol. 18, 635–647. doi:10.1038/s41577-018-0044-0
Bossowski, A., Borysewicz-Sanczyk, H., Wawrusiewicz-Kurylonek, N., Zasim, A., Szalecki, M., Wikiera, B., et al. (2014). Analysis of Chosen Polymorphisms in FoxP3 Gene in Children and Adolescents with Autoimmune Thyroid Diseases. Autoimmunity 47, 395–400. doi:10.3109/08916934.2014.910767
Bossowski, A., Czarnocka, B., Bardadin, K., Stasiak-Barmuta, A., Urban, M., Dadan, J., et al. (2008). Identification of Apoptotic Proteins in Thyroid Gland from Patients with Graves' Disease and Hashimoto's Thyroiditis. Autoimmunity 41, 163–173. doi:10.1080/08916930701727749
Bossowski, A., Moniuszko, M., Idzkowska, E., Grubczak, K., Singh, P., Bossowska, A., et al. (2016). Decreased Proportions of CD4 + IL17+/CD4 + CD25 + CD127- and CD4 + IL17+/CD4 + CD25 + CD127 - FoxP3+ T Cells in Children with Autoimmune Thyroid Diseases. Autoimmunity 49, 320–328. doi:10.1080/08916934.2016.1183654
Brand, O. J., Lowe, C. E., Heward, J. M., Franklyn, J. A., Cooper, J. D., Todd, J. A., et al. (2007). Association of the Interleukin-2 Receptor Alpha (IL-2Ralpha)/CD25 Gene Region with Graves' Disease Using a Multilocus Test and Tag SNPs. Clin. Endocrinol. (Oxf) 66, 508–512. doi:10.1111/j.1365-2265.2007.02762.x
Bufalo, N. E., Dos Santos, R. B., Rocha, A. G., Teodoro, L., Romaldini, J. H., and Ward, L. S. (2021). Polymorphisms of the Genes CTLA4, PTPN22, CD40, and PPARG and Their Roles in Graves' Disease: Susceptibility and Clinical Features. Endocrine 71, 104–112. doi:10.1007/s12020-020-02337-x
Cai, T. T., Muhali, F. S., Song, R. H., Qin, Q., Wang, X., Shi, L. F., et al. (2015). Genome-wide DNA Methylation Analysis in Graves' Disease. Genomics 105, 204–210. doi:10.1016/j.ygeno.2015.01.001
Chang, S. C., Lin, S. F., Chen, S. T., Chang, P. Y., Yeh, Y. M., Lo, F. S., et al. (2021). Alterations of Gut Microbiota in Patients with Graves' Disease. Front Cel Infect Microbiol 11, 663131. doi:10.3389/fcimb.2021.663131
Chen, J., Tian, J., Tang, X., Rui, K., Ma, J., Mao, C., et al. (2015). MiR-346 Regulates CD4(+)CXCR5(+) T Cells in the Pathogenesis of Graves' Disease. Endocrine 49, 752–760. doi:10.1007/s12020-015-0546-5
Chen, J., Wang, W., Guo, Z., Huang, S., Lei, H., Zang, P., et al. (2021). Associations between Gut Microbiota and Thyroidal Function Status in Chinese Patients with Graves' Disease. J. Endocrinol. Invest. 44, 1913–1926. doi:10.1007/s40618-021-01507-6
Chen, X., Hu, Z., Liu, M., Li, H., Liang, C., Li, W., et al. (2018). Correlation between CTLA-4 and CD40 Gene Polymorphisms and Their Interaction in Graves' Disease in a Chinese Han Population. BMC Med. Genet. 19, 171. doi:10.1186/s12881-018-0665-y
Chen, X. J., Gong, X. H., Yan, N., Meng, S., Qin, Q., Jiang, Y. F., et al. (2015). RNASET2 Tag SNP but Not CCR6 Polymorphisms Is Associated with Autoimmune Thyroid Diseases in the Chinese Han Population. BMC Med. Genet. 16, 11. doi:10.1186/s12881-015-0150-9
Chen, Z., Liu, Y., Hu, S., Zhang, M., Shi, B., and Wang, Y. (2021). Decreased Treg Cell and TCR Expansion Are Involved in Long-Lasting Graves' Disease. Front. Endocrinol. (Lausanne) 12, 632492. doi:10.3389/fendo.2021.632492
Chistiakov, D. A., Chistiakova, E. I., Voronova, N. V., Turakulov, R. I., and Savost'anov, K. V. (2011). A Variant of the Il2ra/Cd25 Gene Predisposing to Graves' Disease Is Associated with Increased Levels of Soluble Interleukin-2 Receptor. Scand. J. Immunol. 74, 496–501. doi:10.1111/j.1365-3083.2011.02608.x
Cho, W. K., Jung, M. H., Choi, E. J., Choi, H. B., Kim, T. G., and Suh, B. K. (2011). Association of HLA Alleles with Autoimmune Thyroid Disease in Korean Children. Horm. Res. Paediatr. 76, 328–334. doi:10.1159/000331134
Christensen, N. J., Habekost, G., and Bratholm, P. (2008). A RNA Transcript (Heg) in Mononuclear Cells Is Negatively Correlated with CD14 mRNA and TSH Receptor Autoantibodies. Clin. Exp. Immunol. 154, 209–215. doi:10.1111/j.1365-2249.2008.03744.x
Christensen, N. J., Habekost, G., and Bratholm, P. (2011). Decrease in TSH Receptor Autoantibodies during Antithyroid Treatment: Relationship with a Long Noncoding Heg RNA and Cdk1 mRNA in Mononuclear Cells. ISRN Endocrinol. 2011, 287052. doi:10.5402/2011/287052
Chu, X., Pan, C. M., Zhao, S. X., Liang, J., Gao, G. Q., Zhang, X. M., et al. (2011). A Genome-wide Association Study Identifies Two New Risk Loci for Graves' Disease. Nat. Genet. 43, 897–901. doi:10.1038/ng.898
Corapcioglu, D., Tonyukuk, V., Kiyan, M., Yilmaz, A. E., Emral, R., Kamel, N., et al. (2002). Relationship between Thyroid Autoimmunity and Yersinia Enterocolitica Antibodies. Thyroid 12, 613–617. doi:10.1089/105072502320288483
Das, S., Baruah, C., Saikia, A. K., and Bose, S. (2018). Associative Role of HLA-DRB1 SNP Genotypes as Risk Factors for Susceptibility and Severity of Rheumatoid Arthritis: A North-east Indian Population-Based Study. Int. J. Immunogenet. 45, 1–7. doi:10.1111/iji.12347
Du, J., Wang, X., Tan, G., Wei, W., Zhou, F., Liang, Z., et al. (2021). Predisposition to Graves' Disease and Graves' Ophthalmopathy by Genetic Variants of IL2RA. J. Mol. Med. (Berl) 99, 1487–1495. doi:10.1007/s00109-021-02111-0
Duan, J., and Kasper, D. L. (2011). Regulation of T Cells by Gut Commensal Microbiota. Curr. Opin. Rheumatol. 23, 372–376. doi:10.1097/BOR.0b013e3283476d3e
Esfahanian, F., Naimi, E., Doroodgar, F., and Jadali, Z. (2013). Th1/Th2 Cytokines in Patients with Graves' Disease with or without Ophthalmopathy. Iran J. Allergy Asthma Immunol. 12, 168–175.
Eshaghkhani, Y., Sanati, M. H., Nakhjavani, M., Safari, R., Khajavi, A., Ataei, M., et al. (2016). Disturbed Th1 and Th2 Balance in Patients with Graves' Disease. Minerva Endocrinol. 41, 28–36.
Esteghamati, A., Khalilzadeh, O., Mobarra, Z., Anvari, M., Tahvildari, M., Amiri, H. M., et al. (2009). Association of CTLA-4 Gene Polymorphism with Graves' Disease and Ophthalmopathy in Iranian Patients. Eur. J. Intern. Med. 20, 424–428. doi:10.1016/j.ejim.2008.12.005
Fang, W., Zhang, Z., Zhang, J., Cai, Z., Zeng, H., Chen, M., et al. (2015). Association of the CTLA4 Gene CT60/rs3087243 Single-Nucleotide Polymorphisms with Graves' Disease. Biomed. Rep. 3, 691–696. doi:10.3892/br.2015.493
Fouad, H., Yahia, S., Elsaid, A., Hammad, A., Wahba, Y., El-Gilany, A. H., et al. (2019). Oxidative Stress and Vitamin D Receptor BsmI Gene Polymorphism in Egyptian Children with Systemic Lupus Erythematosus: a Single center Study. Lupus 28, 771–777. doi:10.1177/0961203319846380
Fujii, A., Inoue, N., Watanabe, M., Kawakami, C., Hidaka, Y., Hayashizaki, Y., et al. (2017). TSHR Gene Polymorphisms in the Enhancer Regions Are Most Strongly Associated with the Development of Graves' Disease, Especially Intractable Disease, and of Hashimoto's Disease. Thyroid 27, 111–119. doi:10.1089/thy.2016.0345
Furusawa, Y., Obata, Y., Fukuda, S., Endo, T. A., Nakato, G., Takahashi, D., et al. (2013). Commensal Microbe-Derived Butyrate Induces the Differentiation of Colonic Regulatory T Cells. Nature 504, 446–450. doi:10.1038/nature12721
Glossop, J. R., Dawes, P. T., Nixon, N. B., and Mattey, D. L. (2005). Polymorphism in the Tumour Necrosis Factor Receptor II Gene Is Associated with Circulating Levels of Soluble Tumour Necrosis Factor Receptors in Rheumatoid Arthritis. Arthritis Res. Ther. 7, R1227–R1234. doi:10.1186/ar1816
Grice, E. A., and Segre, J. A. (2012). The Human Microbiome: Our Second Genome. Annu. Rev. Genomics Hum. Genet. 13, 151–170. doi:10.1146/annurev-genom-090711-163814
Gu, L. Q., Zhu, W., Zhao, S. X., Zhao, L., Zhang, M. J., Cui, B., et al. (2010). Clinical Associations of the Genetic Variants of CTLA-4, Tg, TSHR, PTPN22, PTPN12 and FCRL3 in Patients with Graves' Disease. Clin. Endocrinol. (Oxf) 72, 248–255. doi:10.1111/j.1365-2265.2009.03617.x
Guo, Q., Wu, D., Yu, H., Bao, J., Peng, S., Shan, Z., et al. (2018). Alterations of Global DNA Methylation and DNA Methyltransferase Expression in T and B Lymphocytes from Patients with Newly Diagnosed Autoimmune Thyroid Diseases after Treatment: A Follow-Up Study. Thyroid 28, 377–385. doi:10.1089/thy.2017.0301
Han, H., Fu, X., Huang, J., Zhang, X., and Yu, J. (2020). PD-1/PD-L1 Affects Graves Progression through Lymphocytes on the Proliferation, Apoptosis and Inflammatory Cytokine Secretion of Thyroid Follicular Epithelial Cells. J. Toxicol. Sci. 45, 701–711. doi:10.2131/jts.45.701
Hashemi, V., Farrokhi, A. S., Tanomand, A., Babaloo, Z., Hojjat-Farsangi, M., Anvari, E., et al. (2018). Polymorphism of Foxp3 Gene Affects the Frequency of Regulatory T Cells and Disease Activity in Patients with Rheumatoid Arthritis in Iranian Population. Immunol. Lett. 204, 16–22. doi:10.1016/j.imlet.2018.10.001
Hesarghatta Shyamasunder, A., and Abraham, P. (2017). Measuring TSH Receptor Antibody to Influence Treatment Choices in Graves' Disease. Clin. Endocrinol. (Oxf) 86, 652–657. doi:10.1111/cen.13327
Hiratsuka, I., Yamada, H., Munetsuna, E., Hashimoto, S., and Itoh, M. (2016). Circulating MicroRNAs in Graves' Disease in Relation to Clinical Activity. Thyroid 26, 1431–1440. doi:10.1089/thy.2016.0062
Hu, R. T., Liu, D. S., and Li, B. (2020). Predictive Factors for Early Hypothyroidism Following the Radioactive Iodine Therapy in Graves' Disease Patients. BMC Endocr. Disord. 20, 76. doi:10.1186/s12902-020-00557-w
Inoue, N., Watanabe, M., Ishido, N., Katsumata, Y., Kagawa, T., Hidaka, Y., et al. (2014). The Functional Polymorphisms of VDR, GC and CYP2R1 Are Involved in the Pathogenesis of Autoimmune Thyroid Diseases. Clin. Exp. Immunol. 178, 262–269. doi:10.1111/cei.12420
Inoue, N., Watanabe, M., Ishido, N., Kodu, A., Maruoka, H., Katsumata, Y., et al. (2016). Involvement of Genes Encoding Apoptosis Regulatory Factors (FAS, FASL, TRAIL, BCL2, TNFR1 and TNFR2) in the Pathogenesis of Autoimmune Thyroid Diseases. Hum. Immunol. 77, 944–951. doi:10.1016/j.humimm.2016.07.232
Inoue, N., Watanabe, M., Yamada, H., Takemura, K., Hayashi, F., Yamakawa, N., et al. (2012). Associations between Autoimmune Thyroid Disease Prognosis and Functional Polymorphisms of Susceptibility Genes, CTLA4, PTPN22, CD40, FCRL3, and ZFAT, Previously Revealed in Genome-wide Association Studies. J. Clin. Immunol. 32, 1243–1252. doi:10.1007/s10875-012-9721-0
Ishaq, H. M., Mohammad, I. S., Shahzad, M., Ma, C., Raza, M. A., Wu, X., et al. (2018). Molecular Alteration Analysis of Human Gut Microbial Composition in Graves' Disease Patients. Int. J. Biol. Sci. 14, 1558–1570. doi:10.7150/ijbs.24151
Jeltsch, A., and Jurkowska, R. Z. (2014). New Concepts in DNA Methylation. Trends Biochem. Sci. 39, 310–318. doi:10.1016/j.tibs.2014.05.002
Ji, X., Wan, J., Chen, R., Wang, H., Huang, L., Wang, S., et al. (2020). Low Frequency of IL-10-producing B Cells and High Density of ILC2s Contribute to the Pathological Process in Graves' Disease, Which May Be Related to Elevated-TRAb Levels. Autoimmunity 53, 78–85. doi:10.1080/08916934.2019.1698553
Jiang, W., Yu, X., Kosik, R. O., Song, Y., Qiao, T., Tong, J., et al. (2021). Gut Microbiota May Play a Significant Role in the Pathogenesis of Graves' Disease. Thyroid 31, 810–820. doi:10.1089/thy.2020.0193
Jurecka-Lubieniecka, B., Ploski, R., Kula, D., Szymanski, K., Bednarczuk, T., Ambroziak, U., et al. (2014). Association between Polymorphisms in the TSHR Gene and Graves' Orbitopathy. PLoS One 9, e102653. doi:10.1371/journal.pone.0102653
Karagianni, P., and Tzioufas, A. G. (2019). Epigenetic Perspectives on Systemic Autoimmune Disease. J. Autoimmun. 104, 102315. doi:10.1016/j.jaut.2019.102315
Kim, D., Yoo, S. A., and Kim, W. U. (2016). Gut Microbiota in Autoimmunity: Potential for Clinical Applications. Arch. Pharm. Res. 39, 1565–1576. doi:10.1007/s12272-016-0796-7
Kyrgios, I., Fragou, A., Kotanidou, E. P., Mouzaki, K., Efraimidou, S., Tzimagiorgis, G., et al. (2020). DNA Methylation Analysis within the IL2RA Gene Promoter in Youth with Autoimmune Thyroid Disease. Eur. J. Clin. Invest. 50, e13199. doi:10.1111/eci.13199
Lane, L. C., Allinson, K. R., Campbell, K., Bhatnagar, I., Ingoe, L., Razvi, S., et al. (2019). Analysis of BAFF Gene Polymorphisms in UK Graves' Disease Patients. Clin. Endocrinol. (Oxf) 90, 170–174. doi:10.1111/cen.13872
Lee, G. R. (2018). The Balance of Th17 versus Treg Cells in Autoimmunity. Int. J. Mol. Sci. 19, 730. doi:10.3390/ijms19030730
Li, B., Carey, M., and Workman, J. L. (2007). The Role of Chromatin during Transcription. Cell 128, 707–719. doi:10.1016/j.cell.2007.01.015
Li, C., Yuan, J., Zhu, Y. F., Yang, X. J., Wang, Q., Xu, J., et al. (2016). Imbalance of Th17/Treg in Different Subtypes of Autoimmune Thyroid Diseases. Cell Physiol Biochem 40, 245–252. doi:10.1159/000452541
Li, L., Ding, X., Wang, X., Yao, Q., Shao, X., An, X., et al. (2018). Polymorphisms of IKZF3 Gene and Autoimmune Thyroid Diseases: Associated with Graves' Disease but Not with Hashimoto's Thyroiditis. Cel Physiol Biochem 45, 1787–1796. doi:10.1159/000487870
Li, Q., Wang, B., Mu, K., and Zhang, J. A. (2019). The Pathogenesis of Thyroid Autoimmune Diseases: New T Lymphocytes - Cytokines Circuits beyond the Th1-Th2 Paradigm. J. Cel Physiol 234, 2204–2216. doi:10.1002/jcp.27180
Li, Y., Yao, Y., Yang, M., Shi, L., Li, X., Yang, Y., et al. (2013). Association between HLA-B*46 Allele and Graves Disease in Asian Populations: a Meta-Analysis. Int. J. Med. Sci. 10, 164–170. doi:10.7150/ijms.5158
Limbach, M., Saare, M., Tserel, L., Kisand, K., Eglit, T., Sauer, S., et al. (2016). Epigenetic Profiling in CD4+ and CD8+ T Cells from Graves' Disease Patients Reveals Changes in Genes Associated with T Cell Receptor Signaling. J. Autoimmun. 67, 46–56. doi:10.1016/j.jaut.2015.09.006
Lin, J. D., Wang, Y. H., Fang, W. F., Hsiao, C. J., Chagnaadorj, A., Lin, Y. F., et al. (2016a). Serum BAFF and Thyroid Autoantibodies in Autoimmune Thyroid Disease. Clin. Chim. Acta 462, 96–102. doi:10.1016/j.cca.2016.09.004
Lin, J. D., Yang, S. F., Wang, Y. H., Fang, W. F., Lin, Y. C., Lin, Y. F., et al. (2016b). Analysis of Associations of Human BAFF Gene Polymorphisms with Autoimmune Thyroid Diseases. PLoS One 11, e0154436. doi:10.1371/journal.pone.0154436
Liu, R., Ma, X., Xu, L., Wang, D., Jiang, X., Zhu, W., et al. (2012). Differential microRNA Expression in Peripheral Blood Mononuclear Cells from Graves' Disease Patients. J. Clin. Endocrinol. Metab. 97, E968–E972. doi:10.1210/jc.2011-2982
Liu, T., Sun, J., Wang, Z., Yang, W., Zhang, H., Fan, C., et al. (2017). Changes in the DNA Methylation and Hydroxymethylation Status of the Intercellular Adhesion Molecule 1 Gene Promoter in Thyrocytes from Autoimmune Thyroiditis Patients. Thyroid 27, 838–845. doi:10.1089/thy.2016.0576
Liu, W., Wang, H. N., Gu, Z. H., Yang, S. Y., Ye, X. P., Pan, C. M., et al. (2014). Identification of BACH2 as a Susceptibility Gene for Graves' Disease in the Chinese Han Population Based on a Three-Stage Genome-wide Association Study. Hum. Genet. 133, 661–671. doi:10.1007/s00439-013-1404-2
Liu, Y., Yuan, X., Li, X., Cui, D., and Xie, J. (2018). Constitutive Changes in Circulating Follicular Helper T Cells and Their Subsets in Patients with Graves' Disease. J. Immunol. Res. 2018, 8972572. doi:10.1155/2018/8972572
Long, H., Yin, H., Wang, L., Gershwin, M. E., and Lu, Q. (2016). The Critical Role of Epigenetics in Systemic Lupus Erythematosus and Autoimmunity. J. Autoimmun. 74, 118–138. doi:10.1016/j.jaut.2016.06.020
Luo, L., Cai, B., Liu, F., Hu, X., and Wang, L. (2012). Association of Protein Tyrosine Phosphatase Nonreceptor 22 (PTPN22) C1858T Gene Polymorphism with Susceptibility to Autoimmune Thyroid Diseases: a Meta-Analysis. Endocr. J. 59, 439–445. doi:10.1507/endocrj.ej11-0381
Mahfoudh, W., Bel Hadj Jrad, B., Romdhane, A., and Chouchane, L. (2007). A Polymorphism in FAS Gene Promoter Correlated with Circulating Soluble FAS Levels. Int. J. Immunogenet. 34, 209–212. doi:10.1111/j.1744-313X.2007.00676.x
Martin, S., Dutescu, M. I., Sirbu, A., Barbu, C., Albu, A., Florea, S., et al. (2014). The Clinical Value of Human Leukocyte Antigen HLA-DRB1 Subtypes Associated to Graves' Disease in Romanian Population. Immunol. Invest. 43, 479–490. doi:10.3109/08820139.2014.886261
Martinez-Hernandez, R., Sampedro-Nunez, M., Serrano-Somavilla, A., Ramos-Levi, A. M., de la Fuente, H., Trivino, J. C., et al. (2018). A MicroRNA Signature for Evaluation of Risk and Severity of Autoimmune Thyroid Diseases. J. Clin. Endocrinol. Metab. 103, 1139–1150. doi:10.1210/jc.2017-02318
Martire, S., and Banaszynski, L. A. (2020). The Roles of Histone Variants in fine-tuning Chromatin Organization and Function. Nat. Rev. Mol. Cel Biol 21, 522–541. doi:10.1038/s41580-020-0262-8
McLachlan, S. M., and Rapoport, B. (2014). Breaking Tolerance to Thyroid Antigens: Changing Concepts in Thyroid Autoimmunity. Endocr. Rev. 35, 59–105. doi:10.1210/er.2013-1055
Mehraji, Z., Farazmand, A., Esteghamati, A., Noshad, S., Sadr, M., Amirzargar, S., et al. (2017). Association of Human Leukocyte Antigens Class I and II with Graves' Disease in Iranian Population. Iran J. Immunol. 14, 223–230.
Mehta, A., and Baltimore, D. (2016). MicroRNAs as Regulatory Elements in Immune System Logic. Nat. Rev. Immunol. 16, 279–294. doi:10.1038/nri.2016.40
Meng, S., He, S. T., Jiang, W. J., Xiao, L., Li, D. F., Xu, J., et al. (2015). Genetic Susceptibility to Autoimmune Thyroid Diseases in a Chinese Han Population: Role of Vitamin D Receptor Gene Polymorphisms. Ann. Endocrinol. (Paris) 76, 684–689. doi:10.1016/j.ando.2015.01.003
Mizuma, T., Watanabe, M., Inoue, N., Arakawa, Y., Tomari, S., Hidaka, Y., et al. (2017). Association of the Polymorphisms in the Gene Encoding Thyroglobulin with the Development and Prognosis of Autoimmune Thyroid Disease. Autoimmunity 50, 386–392. doi:10.1080/08916934.2017.1344971
Moravvej, H., Tabatabaei-Panah, P. S., Abgoon, R., Khaksar, L., Sokhandan, M., Tarshaei, S., et al. (2018). Genetic Variant Association of PTPN22, CTLA4, IL2RA, as Well as HLA Frequencies in Susceptibility to Alopecia Areata. Immunol. Invest. 47, 666–679. doi:10.1080/08820139.2018.1480032
Morita, E., Watanabe, M., Inoue, N., Hashimoto, H., Haga, E., Hidaka, Y., et al. (2018). Methylation Levels of the TNFA Gene Are Different between Graves' and Hashimoto's Diseases and Influenced by the TNFA Polymorphism. Autoimmunity 51, 118–125. doi:10.1080/08916934.2018.1448078
Mustafa, S., Fatima, H., Fatima, S., Khosa, T., Akbar, A., Shaikh, R. S., et al. (2018). Association of Single Nucleotide Polymorphism in CD28(C/T-I3 + 17) and CD40 (C/T-1) Genes with the Graves' Disease. J. Pak Med. Assoc. 68, 3–9.
Nawrocki, M. J., Majewski, D., Puszczewicz, M., and Jagodzinski, P. P. (2017). Decreased mRNA Expression Levels of DNA Methyltransferases Type 1 and 3A in Systemic Lupus Erythematosus. Rheumatol. Int. 37, 775–783. doi:10.1007/s00296-017-3711-8
Noack, M., and Miossec, P. (2014). Th17 and Regulatory T Cell Balance in Autoimmune and Inflammatory Diseases. Autoimmun. Rev. 13, 668–677. doi:10.1016/j.autrev.2013.12.004
Otsu, H., Watanabe, M., Inoue, N., Masutani, R., and Iwatani, Y. (2017). Intraindividual Variation of microRNA Expression Levels in Plasma and Peripheral Blood Mononuclear Cells and the Associations of These Levels with the Pathogenesis of Autoimmune Thyroid Diseases. Clin. Chem. Lab. Med. 55, 626–635. doi:10.1515/cclm-2016-0449
Park, M. H., Park, Y. J., Song, E. Y., Park, H., Kim, T. Y., et al. (2005). Association of HLA-DR and -DQ Genes with Graves Disease in Koreans. Hum. Immunol. 66, 741–747. doi:10.1016/j.humimm.2005.03.001
Patel, H., Mansuri, M. S., Singh, M., Begum, R., Shastri, M., and Misra, A. (2016). Association of Cytotoxic T-Lymphocyte Antigen 4 (CTLA4) and Thyroglobulin (TG) Genetic Variants with Autoimmune Hypothyroidism. PLoS One 11, e0149441. doi:10.1371/journal.pone.0149441
Picascia, A., Grimaldi, V., Pignalosa, O., De Pascale, M. R., Schiano, C., and Napoli, C. (2015). Epigenetic Control of Autoimmune Diseases: from Bench to Bedside. Clin. Immunol. 157, 1–15. doi:10.1016/j.clim.2014.12.013
Ponting, C. P., Oliver, P. L., and Reik, W. (2009). Evolution and Functions of Long Noncoding RNAs. Cell 136, 629–641. doi:10.1016/j.cell.2009.02.006
Prabhakar, B. S., Bahn, R. S., and Smith, T. J. (2003). Current Perspective on the Pathogenesis of Graves' Disease and Ophthalmopathy. Endocr. Rev. 24, 802–835. doi:10.1210/er.2002-0020
Qin, J., Zhou, J., Fan, C., Zhao, N., Liu, Y., Wang, S., et al. (2017). Increased Circulating Th17 but Decreased CD4(+)Foxp3(+) Treg and CD19(+)CD1d(hi)CD5(+) Breg Subsets in New-Onset Graves' Disease. Biomed. Res. Int. 2017, 8431838. doi:10.1155/2017/8431838
Qin, Q., Wang, X., Yan, N., Song, R. H., Cai, T. T., Zhang, W., et al. (2015). Aberrant Expression of miRNA and mRNAs in Lesioned Tissues of Graves' Disease. Cel Physiol Biochem 35, 1934–1942. doi:10.1159/000374002
Qiu, Y. Y., Wu, Y., Lin, M. J., Bian, T., Xiao, Y. L., and Qin, C. (2019). LncRNA-MEG3 Functions as a Competing Endogenous RNA to Regulate Treg/Th17 Balance in Patients with Asthma by Targeting microRNA-17/RORgammat. Biomed. Pharmacother. 111, 386–394. doi:10.1016/j.biopha.2018.12.080
Quinn, J. J., and Chang, H. Y. (2016). Unique Features of Long Non-coding RNA Biogenesis and Function. Nat. Rev. Genet. 17, 47–62. doi:10.1038/nrg.2015.10
Ramos-Levi, A. M., and Marazuela, M. (2016). Pathogenesis of Thyroid Autoimmune Disease: the Role of Cellular Mechanisms. Endocrinol. Nutr. 63, 421–429. doi:10.1016/j.endonu.2016.04.003
Rosser, E. C., and Mauri, C. (2015). Regulatory B Cells: Origin, Phenotype, and Function. Immunity 42, 607–612. doi:10.1016/j.immuni.2015.04.005
Rydzewska, M., Goralczyk, A., Goscik, J., Wawrusiewicz-Kurylonek, N., Bossowska, A., Kretowski, A., et al. (2018a). Analysis of Chosen Polymorphisms Rs2476601 a/G - PTPN22, Rs1990760 C/T - IFIH1, Rs179247 a/G - TSHR in Pathogenesis of Autoimmune Thyroid Diseases in Children. Autoimmunity 51, 183–190. doi:10.1080/08916934.2018.1486824
Rydzewska, M., Jaromin, M., Pasierowska, I. E., Stozek, K., and Bossowski, A. (2018b). Role of the T and B Lymphocytes in Pathogenesis of Autoimmune Thyroid Diseases. Thyroid Res. 11, 2. doi:10.1186/s13044-018-0046-9
Shalaby, S. M., Mackawy, A. M. H., Atef, D. M., Atef, R. M., and Saeed, J. (2019). Promoter Methylation and Expression of Intercellular Adhesion Molecule 1 Gene in Blood of Autoimmune Thyroiditis Patients. Mol. Biol. Rep. 46, 5345–5353. doi:10.1007/s11033-019-04990-6
Shao, S., Yu, X., and Shen, L. (2018). Autoimmune Thyroid Diseases and Th17/Treg Lymphocytes. Life Sci. 192, 160–165. doi:10.1016/j.lfs.2017.11.026
Shehjar, F., Afroze, D., Misgar, R. A., Malik, S. A., and Laway, B. A. (2018). Association of FoxP3 Promoter Polymorphisms with the Risk of Graves' Disease in Ethnic Kashmiri Population. Gene 672, 88–92. doi:10.1016/j.gene.2018.06.023
Shehjar, F., Dil, A., Misgar, R. A., Malik, S. A., and Laway, B. A. (2020). A Significant Association of the CTLA4 Gene Variants with the Risk of Autoimmune Graves' Disease in Ethnic Kashmiri Population. Cell Immunol 347, 103995. doi:10.1016/j.cellimm.2019.103995
Smith, T. J., and Hegedus, L. (2016). Graves' Disease. N. Engl. J. Med. 375, 1552–1565. doi:10.1056/NEJMra1510030
Song, R. H., Li, Q., Jia, X., Yao, Q. M., Wang, B., and Zhang, J. A. (2018). Polymorphisms of FAM167A-BLK Region Confer Risk of Autoimmune Thyroid Disease. DNA Cel Biol 37, 932–940. doi:10.1089/dna.2018.4344
Stefan, M., Jacobson, E. M., Huber, A. K., Greenberg, D. A., Li, C. W., Skrabanek, L., et al. (2011). Novel Variant of Thyroglobulin Promoter Triggers Thyroid Autoimmunity through an Epigenetic Interferon Alpha-Modulated Mechanism. J. Biol. Chem. 286, 31168–31179. doi:10.1074/jbc.M111.247510
Su, X., Yin, X., Liu, Y., Yan, X., Zhang, S., Wang, X., et al. (2020). Gut Dysbiosis Contributes to the Imbalance of Treg and Th17 Cells in Graves' Disease Patients by Propionic Acid. J. Clin. Endocrinol. Metab. 105. doi:10.1210/clinem/dgaa511
Tan, Y., Chen, W., Liu, C., Zheng, X., Guo, A., and Long, J. (2019). Effect of IL-21 on the Balance of Th17 Cells/Treg Cells in the Pathogenesis of Graves' Disease. Endocr. Res. 44, 138–147. doi:10.1080/07435800.2019.1600535
Thursby, E., and Juge, N. (2017). Introduction to the Human Gut Microbiota. Biochem. J. 474, 1823–1836. doi:10.1042/BCJ20160510
Ueki, I., Abiru, N., Kawagoe, K., and Nagayama, Y. (2011). Interleukin 10 Deficiency Attenuates Induction of Anti-TSH Receptor Antibodies and Hyperthyroidism in a Mouse Graves' Model. J. Endocrinol. 209, 353–357. doi:10.1530/JOE-11-0129
Veneti, S., Anagnostis, P., Adamidou, F., Artzouchaltzi, A. M., Boboridis, K., and Kita, M. (2019). Association between Vitamin D Receptor Gene Polymorphisms and Graves' Disease: a Systematic Review and Meta-Analysis. Endocrine 65, 244–251. doi:10.1007/s12020-019-01902-3
Venkatesh, S., and Workman, J. L. (2015). Histone Exchange, Chromatin Structure and the Regulation of Transcription. Nat. Rev. Mol. Cel Biol 16, 178–189. doi:10.1038/nrm3941
Wang, B., Jia, X., Yao, Q., Li, Q., He, W., Li, L., et al. (2019). CEP128 Is a Crucial Risk Locus for Autoimmune Thyroid Diseases. Mol. Cel Endocrinol 480, 97–106. doi:10.1016/j.mce.2018.10.017
Wang, B., Shao, X., Song, R., Xu, D., and Zhang, J. A. (2017a). The Emerging Role of Epigenetics in Autoimmune Thyroid Diseases. Front. Immunol. 8, 396. doi:10.3389/fimmu.2017.00396
Wang, D., Chen, J., Zhang, H., Zhang, F., Yang, L., and Mou, Y. (2017b). Role of Different CD40 Polymorphisms in Graves' Disease and Hashimoto's Thyroiditis. Immunol. Invest. 46, 544–551. doi:10.1080/08820139.2017.1319382
Wemeau, J. L., Klein, M., Sadoul, J. L., Briet, C., and Velayoudom-Cephise, F. L. (2018). Graves' Disease: Introduction, Epidemiology, Endogenous and Environmental Pathogenic Factors. Ann. Endocrinol. (Paris) 79, 599–607. doi:10.1016/j.ando.2018.09.002
Wu, C. J., Guo, J., Luo, H. C., Wei, C. D., Wang, C. F., Lan, Y., et al. (2016). Association of CD40 Polymorphisms and Haplotype with Risk of Systemic Lupus Erythematosus. Rheumatol. Int. 36, 45–52. doi:10.1007/s00296-015-3345-7
Wu, G. C., Pan, H. F., Leng, R. X., Wang, D. G., Li, X. P., Li, X. M., et al. (2015). Emerging Role of Long Noncoding RNAs in Autoimmune Diseases. Autoimmun. Rev. 14, 798–805. doi:10.1016/j.autrev.2015.05.004
Xiaoheng, C., Yizhou, M., Bei, H., Huilong, L., Xin, W., Rui, H., et al. (2017). General and Specific Genetic Polymorphism of Cytokines-Related Gene in AITD. Mediators Inflamm. 2017, 3916395. doi:10.1155/2017/3916395
Xu, H., Liu, M., Cao, J., Li, X., Fan, D., Xia, Y., et al. (2019). The Dynamic Interplay between the Gut Microbiota and Autoimmune Diseases. J. Immunol. Res. 2019, 7546047. doi:10.1155/2019/7546047
Xuan, M., Zhao, S. X., Yan, C. Y., Yang, J., Li, Y., Song, L. G., et al. (2019). Fine Mapping of Thyroglobulin Gene Identifies Two Independent Risk Loci for Graves' Disease in Chinese Han Population. Ann. Transl Med. 7, 434. doi:10.21037/atm.2019.08.115
Yamada, H., Itoh, M., Hiratsuka, I., and Hashimoto, S. (2014). Circulating microRNAs in Autoimmune Thyroid Diseases. Clin. Endocrinol. (Oxf) 81, 276–281. doi:10.1111/cen.12432
Yan, N., Mu, K., An, X. F., Li, L., Qin, Q., Song, R. H., et al. (2019). Aberrant Histone Methylation in Patients with Graves' Disease. Int. J. Endocrinol. 2019, 1454617. doi:10.1155/2019/1454617
Yan, N., Zhou, J. Z., Zhang, J. A., Cai, T., Zhang, W., Wang, Y., et al. (2015). Histone Hypoacetylation and Increased Histone Deacetylases in Peripheral Blood Mononuclear Cells from Patients with Graves' Disease. Mol. Cel Endocrinol 414, 143–147. doi:10.1016/j.mce.2015.05.037
Yang, M., Sun, B., Li, J., Yang, B., Xu, J., Zhou, X., et al. (2019). Alteration of the Intestinal flora May Participate in the Development of Graves' Disease: a Study Conducted Among the Han Population in Southwest China. Endocr. Connect. 8, 822–828. doi:10.1530/EC-19-0001
Yao, Q., Wang, X., He, W., Song, Z., Wang, B., Zhang, J., et al. (2019). Circulating microRNA-144-3p and miR-762 Are Novel Biomarkers of Graves' Disease. Endocrine 65, 102–109. doi:10.1007/s12020-019-01884-2
Yao, R., Ma, Y. L., Liang, W., Li, H. H., Ma, Z. J., Yu, X., et al. (2012). MicroRNA-155 Modulates Treg and Th17 Cells Differentiation and Th17 Cell Function by Targeting SOCS1. PLoS One 7, e46082. doi:10.1371/journal.pone.0046082
Yin, X., Sachidanandam, R., Morshed, S., Latif, R., Shi, R., and Davies, T. F. (2014). mRNA-Seq Reveals Novel Molecular Mechanisms and a Robust Fingerprint in Graves' Disease. J. Clin. Endocrinol. 99, E2076–E2083. doi:10.1210/jc.2014-1735
Yuan, Q., Zhao, Y., Zhu, X., and Liu, X. (2017). Low Regulatory T Cell and High IL-17 mRNA Expression in a Mouse Graves' Disease Model. J. Endocrinol. Invest. 40, 397–407. doi:10.1007/s40618-016-0575-9
Zhang, D., Qiu, X., Li, J., Zheng, S., Li, L., and Zhao, H. (2019). MiR-23a-3p-regulated Abnormal Acetylation of FOXP3 Induces Regulatory T Cell Function Defect in Graves' Disease. Biol. Chem. 400, 639–650. doi:10.1515/hsz-2018-0343
Zhang, J., Ren, M., Zeng, H., Guo, Y., Zhuang, Z., Feng, Z., et al. (2015). Elevated Follicular Helper T Cells and Expression of IL-21 in Thyroid Tissues Are Involved in the Pathogenesis of Graves' Disease. Immunol. Res. 62, 163–174. doi:10.1007/s12026-015-8647-z
Zhang, Q., Liu, S., Guan, Y., Chen, Q., Zhang, Q., and Min, X. (2017). RNASET2, GPR174, and PTPN22 Gene Polymorphisms Are Related to the Risk of Liver Damage Associated with the Hyperthyroidism in Patients with Graves' Disease. J. Clin. Lab. Anal. 32, e22258. doi:10.1002/jcla.22258
Zhang, X. Y., Tang, X. Y., Li, N., Zhao, L. M., Guo, Y. L., Li, X. S., et al. (2018). GAS5 Promotes Airway Smooth Muscle Cell Proliferation in Asthma via Controlling miR-10a/BDNF Signaling Pathway. Life Sci. 212, 93–101. doi:10.1016/j.lfs.2018.09.002
Zhebrun, D., Kudryashova, Y., Babenko, A., Maslyansky, A., Kunitskaya, N., Popcova, D., et al. (2011). Association of PTPN22 1858T/T Genotype with Type 1 Diabetes, Graves' Disease but Not with Rheumatoid Arthritis in Russian Population. Aging (Albany NY) 3, 368–373. doi:10.18632/aging.100305
Zheng, L., Wang, X., Xu, L., Wang, N., Cai, P., Liang, T., et al. (2015). Foxp3 Gene Polymorphisms and Haplotypes Associate with Susceptibility of Graves' Disease in Chinese Han Population. Int. Immunopharmacol 25, 425–431. doi:10.1016/j.intimp.2015.02.020
Zhou, F., Liang, Z., Wang, X., Tan, G., Wei, W., Zheng, G., et al. (2021). The VDR Gene Confers a Genetic Predisposition to Graves' Disease and Graves' Ophthalmopathy in the Southwest Chinese Han Population. Gene 793, 145750. doi:10.1016/j.gene.2021.145750
Keywords: graves’ disease, pathogenesis, genetics, epigenetics, cellular immunology, gut microbiota
Citation: Zhou F, Wang X, Wang L, Sun X, Tan G, Wei W, Zheng G, Ma X, Tian D and Yu H (2022) Genetics, Epigenetics, Cellular Immunology, and Gut Microbiota: Emerging Links With Graves’ Disease. Front. Cell Dev. Biol. 9:794912. doi: 10.3389/fcell.2021.794912
Received: 14 October 2021; Accepted: 09 December 2021;
Published: 04 January 2022.
Edited by:
Wesley H. Brooks, University of South Florida, United StatesReviewed by:
Marcos Edgar Herkenhoff, University of São Paulo, BrazilEguzkine Ochoa, University of Cambridge, United Kingdom
Copyright © 2022 Zhou, Wang, Wang, Sun, Tan, Wei, Zheng, Ma, Tian and Yu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hongsong Yu, yuhongsong@163.com