- 1Academy of Zhongjing, Henan University of Chinese Medicine, Zhengzhou, China
- 2Laboratory of Traditional Chinese Medicine (TCM) Syndrome and Prescription Signaling, Academy of Zhongjing, Zhengzhou, China
- 3Henan Key Laboratory of Traditional Chinese Medicine (TCM) Syndrome and Prescription Signaling, Henan International Joint, Zhengzhou, China
- 4Department of Orthopedic, Peking University Third Hospital, Beijing, China
Background and objectives: Hashimoto’s thyroiditis (HT), a chronic autoimmune disorder impacting thyroid function, is a growing public health concern. The relationship between Treg cells and HT has been extensively studied, with Treg cells considered crucial in suppressing HT progression. However, these studies have mainly been observational, limiting our understanding of Treg cells’ impact on HT risk. Leveraging large datasets, we utilized Mendelian randomization (MR) analysis to examine the causal association between Treg cell biomarkers and HT, providing additional validation for these relationships.
Methods: Comprehensive two-sample Mendelian randomization analysis was performed to determine the causal association between Treg cells signatures and HT in this study. Based on publicly available genetic data, we explored causal associations between 165 Treg cells signatures and HT risk.
Results: The European cohort study has identified five Treg cell phenotypes that causally protect against HT risk. Resting Treg %CD4 (OR = 0.975, 95% CI = 0.954~0.998, P = 0.030); CD4 on resting Treg (OR = 0.938, 95% CI = 0.882~0.997, P = 0.041; CD28- CD8dim %CD8dim (OR = 0.983, 95% CI = 0.969~0.998, P = 0.030); CD25 on CD39+ resting Treg (OR = 0.926, 95% CI = 0.864~0.991, P = 0.026); 5) CD28 on activated & secreting Treg (OR = 0.969, 95% CI = 0.942~0.996, P = 0.025). The Asian cohort study has identified four Treg cell phenotypes negatively correlated with the risk of HT. CD25hi %T cell (OR = 0.635, 95% CI = 0.473~852, P = 0.002); CD4 Treg %CD4 (OR = 0.829, 95% CI = 0.687~1.000, P = 0.050); CD127-CD8br %T cell (OR = 0.463, 95% CI =0.311~0.687, P< 0.001); CD3 on resting Treg (OR = 0.786, 95% CI = 0.621~0.994, P = 0.044).
Conclusion: Our study has demonstrated the close connection between Treg cells and HT by genetic means, thus providing foundational basis for future research.
1 Introduction
Hashimoto’s Thyroiditis (HT) is a common chronic autoimmune disease affecting thyroid function (1, 2), named after the Japanese scientist Hashimoto, who first described it as a type of autoimmune thyroiditis typically resulting in hypothyroidism (3). The primary pathological changes include extensive enlargement of the thyroid gland, follicular hyperplasia, and the presence of plasma cells and lymphocytes (4). The incidence of HT has been rapidly increasing worldwide in recent years (4), with rates ranging between 27 and 448 per 100,000 (5). Consequently, HT has become a critical public health issue. And studies have shown that, HT is a known risk factor for thyroid cancer development (6). Therefore, Precise diagnosis and effective treatment of HT are essential for enhancing patients’ quality of life and reducing their risk of thyroid cancer (7).
Despite the ongoing debate regarding the etiology and pathogenesis of HT, there is a prevailing agreement that autoimmunity maintains a crucial role in its pathogenesis (2, 8). Recently, regulatory T cells (Treg) have emerged as crucial factors in HT. Treg cells, a CD4 + T cell population, exert immunosuppressive effects on excessive immune responses (9–11). Animal experiments (12–14) indicate that the lack of Tregs could result in the spontaneous advancement of autoimmune thyroiditis. Clinical evidence (10, 15) further substantiates that Treg frequency in circulation is reduced among HT patients. These findings suggest that a decrease in Treg numbers or dysfunction leads to impaired immune tolerance towards autoantigenic stimuli, contributing to the development of Hashimoto’s thyroiditis. Consequently, Treg immunomodulation has emerged as a crucial matter in recent HT immunology research. However, conflicting findings exist in other studies. Glick et al. reported no significant difference in CD4+CD25+ Tregs frequency in HT (16). Even one study found a significantly higher proportion of CD4+Foxp3+ Tregs in autoimmune thyroid disease (AITD) patients (17). The variation may be ascribed to the dual regulatory function of the immune system and the heterogeneity induced by the multiple Treg cell subtypes (18). Therefore, investigating the association between Treg cell subtypes and HT could potentially have a crucial role in precise immunotherapy for HT.
Mendelian randomization (MR) utilizes single nucleotide polymorphisms (SNPs) as instrumental variables (IVs) to infer potential causal relationships between exposure and outcomes (19, 20). follows the Mendelian gamete random allocation principle as well as the principle of free combination. This approach effectively mitigates the impact of confounding factors and reverse causation commonly encountered in conventional observational studies (21). MR is widely used in genetic epidemiology (22). With the large datasets currently available, MR analysis may be a good tool for exploring novel biomarkers in Treg cells that have a causal impact on the risk of Hashimoto’s thyroiditis. Therefore, this study employs a two-sample Mendelian randomization approach, utilizing large-scale GWAS databases to explore the causal relationship between immune cell phenotypes and HT risk, To our knowledge, there haven’t been previous studies using MR to assess the causal impact of Treg cell traits on the risk of HT.
2 Materials and methods
2.1 Study design
The current study applied a two-sample MR method to analyze causal relationships between 165 Treg traits and HT. Therefore, to ensure the validity of causal inference, the IVs relies on three fundamental assumptions: 1) Genetic variants have a direct connection with the exposure; 2) Genetic variants are not associated with potential confounding factors between exposure and outcome; and 3) Genetic variants do not influence the outcome through pathways unrelated to the exposure. All MR analyses were conducted using publicly available summary statistics, thus eliminating the need for additional ethical approval or informed consent.
2.2 Data source
The GWAS summary statistics for a diverse range of Treg traits were retrieved from the GWAS Catalog (23), the GWAS Catalog Accession Numbers are detailed in Supplementary Table 1. The original GWAS for immune traits utilized data from 3,757 European individuals, these datasets encompass an extensive collection of 165 Treg traits. The GWAS summary statistics for HT were sourced from the GWAS Catalog, encompassing datasets from both European (accession number GCST90018855) (24) and Asian populations (accession number GCST90018635). The initial GWAS analysis for the European cohort involved 395,640 individuals, including 15,654 cases and 379,986 controls. Likewise, the original GWAS analysis for the Asian cohort included 173,193 participants of East Asian descent, predominantly of Japanese origin, with 537 cases and 172,656 controls.
2.3 Instrumental variable selection
Genetic variants were utilized as IVs in this study. Significant SNPs that were associated with the Treg cell phenotype with genome-wide significance (P<1×10-5) (25) were selected based on recent research. To minimize bias due to linkage disequilibrium (LD), the following conditions were imposed: R2< 0.001 and distance = 10,000 kb (26). IVs with an F< 10 were excluded to ensure strong correlation of SNPs with HT risk (Supplementary Table 2). The formula for calculating F is: F = , where K is the sample size, R2 is the proportion of variance explained by SNPs in the exposure database (27). The formula for calculating R2 is: , where eaf is the effect allele frequency, β is the allele effect value, K is the sample size, and se is the Standard Error. The exposure and outcome datasets were merged, and incompatible alleles and palindromic SNPs were eliminated. The phenotypes associated with the remaining SNPs were identified utilizing the Human Genotype-Phenotype Association Database (http://www.phenoscanner.medschl.cam.ac.uk/) (28). The remaining SNPs were ultimately utilized as instrumental variables for exposure.
2.4 Statistical analysis
The univariable MR analysis was performed using the R packages “TwoSampleMR”. Our analysis primarily utilized the inverse variance weighting (IVW) method (29). We employed four additional methods to assess the reliability and consistency of our results. These methods include MR-Egger regression (30), the weighted median estimator (WME) (31), mode-based simple estimation (32), and mode-based weighted estimation. All causal estimates were converted to odds ratios (ORs) for the outcome which was a dichotomous phenotype.
Heterogeneity among the individual estimates of genetic variation was evaluated using The Cochran Q test. No significant heterogeneity was observed when P > 0.05. We used the MR-Egger intercept to assess the effect of horizontal pleiotropy (33). Sensitivity analysis was conducted using the “leave-one-out” method to assess the impact of individual SNPs on the causal relationship. The Mendelian Randomization Pleiotropy RESidual Sum and Outlier (MR-PRESSO) method were utilized as indicators to evaluate and correct the level of pleiotropy (34), the MR-PRESSO was conducted using the R package “MRPRESSO” (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1347695/fendo-15-1347695-HTML/image_m/fendo-15-1347695-g001.jpg)
Figure 1 The steps of Mendelian randomization (MR)analysis. SNP, single-nucleotide polymorphism; IVW, inverse variance weighting; MR-PRESSO, Mendelian Randomization Pleiotropy RESidual Sum and Outlier.
3 Results
3.1 The MR analysis of results from European cohorts
The results of the European cohorts study revealed nine Treg cell phenotypes significantly linked with HT risk. Figure 2 presents these findings. Of the identified Treg cell features, we found five that were protective against HT risk: 1) CD25 on CD39+ resting Treg (IVW OR = 0.926, 95% CI = 0.864~0.991, P = 0.026); 2) CD28- CD8dim %CD8dim (IVW OR = 0.983, 95% CI = 0.969~0.998, P = 0.030); 3)CD28 on activated & secreting Treg (IVW OR = 0.969, 95% CI = 0.942~0.996, P = 0.025); 4)CD4 on resting Treg (IVW OR = 0.938, 95% CI = 0.882~0.997, P = 0.041); 5) Resting Treg %CD4 (IVW OR = 0.975, 95% CI = 0.954~0.998, P = 0.030), as well as four other Treg cell subtypes (CD3 on CD39+ secreting Treg, CD3 on CD4 Treg, CD3 on CD45RA+ CD4+, CD39+ activated Treg %activated Treg). Reverse MR analysis revealed a negative causal relationship between HT and CD25 on CD39+ resting Treg (IVW OR = 0.891, 95% CI = 0.798~0.996, P = 0.042) (see Figure 3). The sensitivity test results are displayed in Table 1. Additionally, the leave-one-out analysis corroborated this finding (see Supplementary Figure 1). Causal effect estimates were depicted using scatter plots (see Supplementary Figure 2). Moreover, funnel plots, providing a visual appraisal, revealed no heterogeneity in the causal effect (see Supplementary Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1347695/fendo-15-1347695-HTML/image_m/fendo-15-1347695-g002.jpg)
Figure 2 Forest plots showing significantly causal relationships between 9 Treg cell traits and HT (European cohorts). nsnp, number of single-nucleotide polymorphism; OR, odds ratio; CI, confidence interval.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1347695/fendo-15-1347695-HTML/image_m/fendo-15-1347695-g003.jpg)
Figure 3 Forest plots showing reverse MR analysis results (European cohorts). nsnp, number of single-nucleotide polymorphism; OR, odds ratio; CI, confidence interval.
3.2 The MR analysis of results from East Asian cohorts
The study’s findings in the East Asian cohorts indicate that seven Treg cell phenotypes are significantly linked to HT risk. Figure 4 demonstrates these associations. Four Treg cell features were identified as beneficial in mitigating HT risk: 1) CD127-CD8br %T cell (OR = 0.463, 95% CI =0.311~0.687, P< 0.001); 2) CD25hi %T cell (OR = 0.635, 95% CI = 0.473~852, P = 0.002); 3) CD3 on resting Treg (OR = 0.786, 95% CI = 0.621~0.994, P = 0.044); 4) CD4 Treg %CD4 (OR = 0.829, 95% CI = 0.687~1.000, P = 0.050). And 3 other Treg cell subtypes (CD25 on CD39+activated Treg, CD39+CD8br %CD8br, CD39+CD8br %T cell).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1347695/fendo-15-1347695-HTML/image_m/fendo-15-1347695-g004.jpg)
Figure 4 Forest plots showing significantly causal relationships between 7 Treg cell traits and HT (East Asian cohorts). nsnp, number of single-nucleotide polymorphism; OR, odds ratio; CI, confidence interval.
Reverse MR analysis revealed no significant causal relationship between HT and Treg cell phenotype. The results of the sensitivity test are shown in Table 2. In the same vein, the leave-one-out test also supported this observation (Supplementary Figure 4). The scatter plots see Supplementary Figure 5, funnel plots see Supplementary Figure 6.
3.3 Cross-validation results comparing the European and East Asian cohorts
Mutual validation was performed between European groups and East Asian groups. The MR results showed a total of five results for both datasets that showed a similar propensity to protect as the preliminary analysis (IVW OR<1), but the significance did not reach statistical significance, as shown in Figure 5. Specifically included were CD25 on CD39+ resting Treg (East Asian OR = 0.759, and European OR = 0.926), CD28- CD8dim %CD8dim (East Asian OR = 0.991, and European OR = 0.983), CD28 on activated & secreting Treg (East Asian OR = 0.991, and European OR = 0.969), Resting Treg %CD4 (East Asian OR = 0.857, and European OR = 0.975), and CD4 Treg %CD4 (East Asian OR = 0.829, and European OR = 0.964). And 4 other Treg cell subtypes (CD4 on resting Treg, CD127-CD8br %T cell, CD25hi %T cell, and CD3 on resting Treg).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1347695/fendo-15-1347695-HTML/image_m/fendo-15-1347695-g005.jpg)
Figure 5 The MR analysis results of mutual validation was performed between European groups and East Asian groups. nsnp, number of SNP. CI, confidence interval.
In conclusion, through MR analysis, we identified five treg cell characteristics of interest that were protective against the risk of HT. They include Resting Treg %CD4, CD4 Treg %CD4, CD28 on activated & secreting Treg, CD28- CD8dim %CD8dim, and CD25 on CD39+ resting Treg. These features exhibit a negative correlation with the risk of HT, suggesting a potential protective role against HT.
4 Discussion
HT is an autoimmune disease in which Treg cells play a crucial role in preventing the onset and progression by suppressing the production of harmful autoantibodies targeting thyroglobulin (TG) and thyroid peroxidase (TPO) (35) by autoreactive B cells, thus preventing the destruction caused by these autoreactive B cells (36, 37). Furthermore, Treg cells protect the thyroid tissue from attack by the release of immunosuppressive molecules such as IL-10 and TGF- β (38). In addition to their role in immune homeostasis, Treg cells play indispensable roles in self-tolerance (39). Numerous animal experiments (12, 13, 40) have demonstrated that the absence of Tregs can trigger the spontaneous development of autoimmune thyroiditis. Meanwhile, clinical evidence (10) suggests a decrease in the frequency of circulating Tregs among individuals diagnosed with HT. However, there are also reports indicating that Treg cells may contribute to the progression of HT (16, 17). Such discrepancies may be attributed to different detection markers or subsets with distinct phenotypes. Considering the dissection of Tregs into subpopulations and limited data reported in autoimmune thyroiditis, in order to delve deeper into the potential Treg cell phenotypes associated with reduced risk of HT, we comprehensively examined the causal relationship between 165 Treg cell phenotypes and HT utilizing public genetic data and employing Mendelian randomization methods (41). MR employs genetic variants as instrumental variables, which are fixed at conception, to perform causal inferences regarding the effects of modifiable risk factors. This approach can effectively mitigate certain types of confounding (42). Notably, this study represented a pioneering analysis utilizing Mendelian randomization to scrutinize the causal association between multiple Treg cell traits and HT. In our study, five protective factors were ultimately identified.
The first trait is Resting Treg %CD4. Previous research has demonstrated that by utilizing FoxP3 and CD45RA, three distinct subpopulations can be differentiated based on both phenotypic and functional characteristics: the CD45RA + FoxP3lo resting Treg (rTreg), the CD45RA-FoxP3hi effector Treg (eTreg), and the CD45RA-FoxP3lo cells (Fraction III, FIII). Resting Treg cells have the ability to suppress autoimmune reactions (43), which prevents the immune system from launching attacks against its own tissues. Resting-phase Treg has been demonstrated to decrease in autoimmune disease, indicating a potential role for Resting Treg in mitigating autoimmune disease (44, 45). Research on resting Tregs in HT is limited; currently, only one clinical study has reported no significant differences in the proportion of resting Tregs between HT patients and healthy volunteers, and there is no direct research establishing a link between resting Treg %CD4 and HT. However, HT is considered an autoimmune disease, and our MR analysis results also indicate a negative correlation between the percentage of resting Treg %CD4 and the risk of HT, suggesting a potential protective role of resting Treg %CD4 in HT.
The second protective trait is CD4 Treg %CD4. CD4 is a membrane protein expressed extensively on the surface of T cells, which is indispensable for recognizing most antigens in the body. The CD4 co-receptor augments the T cell sensitivities to antigens by 30- to 100-fold (46), and its primary function is to bind with MHC-II (major histocompatibility complex class II) (47) molecules that are present on the surfaces of antigen-presenting cells, mainly associating with the β-chain binding. This association allows CD4+ T cells to closely associate with antigen-presenting cells, which initiates T cell activation and guides their immune response (48, 49). This process is crucial for regulating the immune system’s normal functionality and preventing excessive immune responses. Hence, the activation and execution of Treg cells’ biological functions are integrally linked to CD4 molecules. Experimental in vivo studies have indicated a significant reduction in the proportion of CD4CD25 Treg cells in patients with HT compared to healthy controls (HCs) (50). One meta-analysis showed that newly diagnosed HT patients have lower levels of Tregs compared to HCs (51). Another meta-analysis investigated the percentage of peripheral blood Tregs among CD4+ T cells in patients with AITDs. The results revealed that the proportions of Tregs among CD4+ T cells in untreated AITDs were significantly lower than those in HCs. However, unfortunately, no significant differences were observed in untreated HT patients (52). Although our findings suggest a negative correlation between CD4 Treg %CD4 and HT, further studies are warranted to establish and validate this relationship, considering the limited research on the association between CD4 Treg %CD4 and HT.
The third and fourth involve CD28 on activated & secreting Treg, and CD28- CD8dim %CD8dim. CD28 is a type I transmembrane glycoprotein receptor typically expressed on most CD4+ T cells (53). It serves as a co-stimulatory receptor expressed differently at specific stages of T cell differentiation and within specific T cell subgroups to regulate immune responses (54). T cell activation involves two signals. The first signal involves the T cell receptor (TCR) recognizing and binding specifically to the antigenic peptide-major histocompatibility complex (MHC) on antigen-presenting cells. The second signal results from the interaction between T cells and a range of co-stimulatory proteins expressed on the surface of antigen-presenting cells, generating co-stimulatory signals (55). It is worth noting that the most crucial interaction among these proteins is the binding of CD28 to CD80/CD86 (56). There is ongoing clinical development of a solitary CD28 agonist antibody, specifically TAB 08 (TGN 1412) (57, 58). However, the potential risks and benefits of this immunological agonist remain unclear and require further study (59). Additionally, CD28 superagonistic monoclonal antibodies (CD28SAB) have been found to preferentially activate and expand immuno-suppressive regulatory T (Treg) cells. In preclinical tests, CD28SABs have demonstrated the ability to alleviate the pathogenesis, progression, and clinical course of autoimmune diseases (60). Although our results revealed a negative association between CD28 on activated & secreting Treg and the risk of HT, which is consistent with previous findings, it is noteworthy that the CD28- CD8dim %CD8dim phenotype was also negatively associated with the risk of HT. Previous studies indicated that since CD28 was expressed not only on Treg but also on effector T cells, over-activation of effector T cells aggravated the progression of thyroid tissue injury or inflammation (61). Hyperreactivity of B7/CD28 signaling is closely linked with the development of autoimmune diseases (62). Cytotoxic T-Lymphocyte Antigen-4 (CTLA-4) is an inhibitory protein that competitively binds to B7 molecules (CD80 and CD86) in competition with CD28 (63), exerting immune co-suppressive responses (64), making it a highly regarded therapeutic target in the treatment of autoimmune diseases. It has been reported that using CTLA-4 fusion protein to block the interaction between CD28 and B7 can inhibit the autoimmune response in experimental rheumatoid arthritis (RA) (65). Additionally, inhibiting the CD80-CD28-mediated signaling pathway via specific RNAi and antibodies to block the interaction between CD28 and B7 reduced immune cell activation and reversed the damage in a pathological lupus nephritis model (66). Therefore, We believe that solely activating or inhibiting the CD28 molecule may not be ideal for improving HT.
The final parameter is CD25 on CD39+ resting Treg. The α-chain of the interleukin-2 receptor (IL-2R) is CD25, which is expressed mainly on immune-related cells’ surfaces, including activated T cells, B cells, and NK cells. It is another Treg surface marker that interacts with the ligand IL-2 and is involved in a variety of immune responses (67, 68). In 1995, Sakaguchi et al. identified a subset of CD4+ T lymphocytes expressing abundant CD25 in thymectomized mice aged 3 days and other autoimmune models, demonstrating their capacity to impede disease progression (69). Similarly, the elevated counts of CD4+ CD25 high HLA-DR+ cells in HT patients may suggest compensatory expansion of this Treg subpopulation to suppress the immune response (70). In an in vivo study (71), CD4+CD25+Foxp3+ Tregs were revealed to play a role in the pathogenesis of experimental autoimmune thyroiditis (EAT), a murine model of HT. CD25 expression is indispensable for Tregs’ differentiation and adaptive proliferation. Mice lacking IL2Rα or IL2Rβ chains manifest severe systemic autoimmunity, a condition ameliorated by transferring CD4 + CD25 + T cells from wild-type mice (72). Furthermore, Mice with deficiencies in IL-2 or the IL-2 receptor exhibit diminished Treg generation, precipitating autoimmune disorders (73). Our findings indicate a negative correlation between HT risk and CD25 on CD39+ resting Treg, consistent with previous research. However, the reverse MR analysis revealed a significant causal relationship between HT risk and CD25 on CD39+ resting Treg, which may be related to unknown confounding factors. Our MR analysis attempted to exclude confounding factors, but there may still be unidentified factors that could impact the results.
Our study has inherent limitations. First, we cannot entirely exclude the potential influence of other confounding factors on this association. Considering the possible presence of additional variables, further exploration of the precise role of these Treg cell phenotypes in the pathogenesis of HT is warranted based on our study findings. Second, despite leveraging GWAS data from European and Asian populations, the generalizability of our findings to other ethnic groups remains subject to constraints. Third, the absence of individual-specific data precluded a more nuanced, stratified analysis of the population, introducing potential inaccuracies into our speculative conclusions. Finally, Our study reported on the Treg cell phenotypes that may decrease the risk of HT in both European and Asian populations. We conducted mutual validation in two study cohorts and applied a more lenient FDR correction for p-values. Although this effectively reduces the occurrence of false-negative results, it inevitably leads to an increased probability of type I errors.
5 Conclusions
In summary, we have established a causal association between five distinct Treg cell phenotypes and HT through an extensive two-way Mendelian randomization analysis, elucidating the intricate interplay between the immune system and HT. Furthermore, our research has effectively mitigated the impact of inevitable confounding factors, reverse causality, and other variables. This may serve as a foundational basis for researchers to delve into Treg cells as potential immunotherapeutic targets for HT.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
JG: Data curation, Methodology, Software, Visualization, Writing – original draft. GS: Investigation, Methodology, Software, Validation, Visualization, Writing – review & editing. FS: Conceptualization, Formal Analysis, Funding acquisition, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The research was funded by the Henan science and technology research project (NO.222102310187), Key scientific research projects of universities of Henan Provincial Department of Education (No:21A360018), 2022 Provincial Science and Technology R&D Plan Joint Fund (cultivation of superior disciplines) (No: 222301420024).
Acknowledgments
Thanks to all authors for their contributions.
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/fendo.2024.1347695/full#supplementary-material
References
1. Klubo-Gwiezdzinska J, Wartofsky L. Hashimoto thyroiditis: an evidence-based guide to etiology, diagnosis and treatment. Polish Arch Internal Med. (2022) 132:16222. doi: 10.20452/pamw.16222
2. Weetman AP. An update on the pathogenesis of Hashimoto's thyroiditis. J Endocrinol Invest. (2021) 44:883–90. doi: 10.1007/s40618-020-01477-1
3. Chen P, Li C, Zhao S, Wang L, Liu L, Fan Q. Effect of large dosage of Prunella on Hashimoto's thyroiditis: A protocol of systematic review and meta-analysis of randomized clinical trials. Medicine. (2020) 99:e23391. doi: 10.1097/md.0000000000023391
4. Ralli M, Angeletti D, Fiore M, D'Aguanno V, Lambiase A, Artico M, et al. Hashimoto's thyroiditis: An update on pathogenic mechanisms, diagnostic protocols, therapeutic strategies, and potential Malignant transformation. Autoimmun Rev. (2020) 19:102649. doi: 10.1016/j.autrev.2020.102649
5. Li L, Ying YX, Liang J, Geng HF, Zhang QY, Zhang CR, et al. Urinary iodine and genetic predisposition to Hashimoto's thyroiditis in a Chinese Han population: A case-control study. Thyroid. (2020) 30:1820–30. doi: 10.1089/thy.2020.0094
6. Feldt-Rasmussen U. Hashimoto's thyroiditis as a risk factor for thyroid cancer. Curr Opin Endocrinol diabetes Obes. (2020) 27:364–71. doi: 10.1097/med.0000000000000570
7. Strand SH, Rivero-Gutiérrez B, Houlahan KE, Seoane JA, King LM, Risom T, et al. Molecular classification and biomarkers of clinical outcome in breast ductal carcinoma in situ: Analysis of TBCRC 038 and RAHBT cohorts. Cancer Cell. (2022) 40:1521–1536.e1527. doi: 10.1016/j.ccell.2022.10.021
8. Pyzik A, Grywalska E, Matyjaszek-Matuszek B, Roliński J. Immune disorders in Hashimoto's thyroiditis: what do we know so far? J Immunol Res. (2015) 2015:979167. doi: 10.1155/2015/979167
9. Shao S, Yu X, Shen L. Autoimmune thyroid diseases and Th17/Treg lymphocytes. Life Sci. (2018) 192:160–5. doi: 10.1016/j.lfs.2017.11.026
10. Xue H, Yu X, Ma L, Song S, Li Y, Zhang L, et al. The possible role of CD4+CD25(high)Foxp3+/CD4+IL-17A+ cell imbalance in the autoimmunity of patients with Hashimoto thyroiditis. Endocrine. (2015) 50:665–73. doi: 10.1007/s12020-015-0569-y
11. Zha B, Huang X, Lin J, Liu J, Hou Y, Wu G. Distribution of lymphocyte subpopulations in thyroid glands of human autoimmune thyroid disease. J Clin Lab Anal. (2014) 28:249–54. doi: 10.1002/jcla.21674
12. Saitoh O, Nagayama Y. Regulation of Graves' hyperthyroidism with naturally occurring CD4+CD25+ regulatory T cells in a mouse model. Endocrinology. (2006) 147:2417–22. doi: 10.1210/en.2005-1024
13. Verginis P, Li HS, Carayanniotis G. Tolerogenic semimature dendritic cells suppress experimental autoimmune thyroiditis by activation of thyroglobulin-specific CD4+CD25+ T cells. J Immunol (Baltimore Md. (2005) 1950) 174:7433–9. doi: 10.4049/jimmunol.174.11.7433
14. Wei WZ, Morris GP, Kong YC. Anti-tumor immunity and autoimmunity: a balancing act of regulatory T cells. Cancer Immunol Immunother CII. (2004) 53:73–8. doi: 10.1007/s00262-003-0444-1
15. Li C, Yuan J, Zhu YF, Yang XJ, Wang Q, Xu J, et al. Imbalance of Th17/Treg in different subtypes of autoimmune thyroid diseases. Cell Physiol Biochem. (2016) 40:245–52. doi: 10.1159/000452541
16. Glick AB, Wodzinski A, Fu P, Levine AD, Wald DN. Impairment of regulatory T-cell function in autoimmune thyroid disease. Thyroid. (2013) 23:871–8. doi: 10.1089/thy.2012.0514
17. Marazuela M, García-López MA, Figueroa-Vega N, de la Fuente H, Alvarado-Sánchez B, Monsiváis-Urenda A, et al. Regulatory T cells in human autoimmune thyroid disease. J Clin Endocrinol Metab. (2006) 91:3639–46. doi: 10.1210/jc.2005-2337
18. Zhang R, Miao J, Zhu P. Regulatory T cell heterogeneity and therapy in autoimmune diseases. Autoimmun Rev. (2021) 20:102715. doi: 10.1016/j.autrev.2020.102715
19. Emdin CA, Khera AV, Kathiresan S. Mendelian randomization. Jama. (2017) 318:1925–6. doi: 10.1001/jama.2017.17219
20. Birney E. Mendelian randomization. Cold Spring Harbor Perspect Med. (2022) 12:a041302. doi: 10.1101/cshperspect.a041302
21. Ference BA, Holmes MV, Smith GD. Using Mendelian randomization to improve the design of randomized trials. Cold Spring Harbor Perspect Med. (2021) 11:a040980. doi: 10.1101/cshperspect.a040980
22. Paternoster L, Tilling K, Davey Smith G. Genetic epidemiology and Mendelian randomization for informing disease therapeutics: Conceptual and methodological challenges. PloS Genet. (2017) 13:e1006944. doi: 10.1371/journal.pgen.1006944
23. Orrù V, Steri M, Sidore C, Marongiu M, Serra V, Olla S, et al. Complex genetic signatures in immune cells underlie autoimmunity and inform therapy. Nat Genet. (2020) 52:1036–45. doi: 10.1038/s41588-020-0684-4
24. Sakaue S, Kanai M, Tanigawa Y, Karjalainen J, Kurki M, Koshiba S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. (2021) 53:1415–24. doi: 10.1038/s41588-021-00931-x
25. Yu XH, Yang YQ, Cao RR, Bo L, Lei SF. The causal role of gut microbiota in development of osteoarthritis. Osteoarthritis Cartilage. (2021) 29:1741–50. doi: 10.1016/j.joca.2021.08.003
26. Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, et al. A global reference for human genetic variation. Nature. (2015) 526:68–74. doi: 10.1038/nature15393
27. Pierce BL, Ahsan H, Vanderweele TJ. Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants. Int J Epidemiol. (2011) 40:740–52. doi: 10.1093/ije/dyq151
28. Kamat MA, Blackshaw JA, Young R, Surendran P, Burgess S, Danesh J, et al. PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinf (Oxford England). (2019) 35:4851–3. doi: 10.1093/bioinformatics/btz469
29. Burgess S, Small DS, Thompson SG. A review of instrumental variable estimators for Mendelian randomization. Stat Methods Med Res. (2017) 26:2333–55. doi: 10.1177/0962280215597579
30. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. (2015) 44:512–25. doi: 10.1093/ije/dyv080
31. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. (2016) 40:304–14. doi: 10.1002/gepi.21965
32. Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. (2017) 46:1985–98. doi: 10.1093/ije/dyx102
33. Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. (2017) 32:377–89. doi: 10.1007/s10654-017-0255-x
34. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
35. Luty J, Ruckemann-Dziurdzińska K, Witkowski JM, Bryl E. Immunological aspects of autoimmune thyroid disease - Complex interplay between cells and cytokines. Cytokine. (2019) 116:128–33. doi: 10.1016/j.cyto.2019.01.003
36. Gotot J, Gottschalk C, Leopold S, Knolle PA, Yagita H, Kurts C, et al. Regulatory T cells use programmed death 1 ligands to directly suppress autoreactive B cells in vivo. Proc Natl Acad Sci United States America. (2012) 109:10468–73. doi: 10.1073/pnas.1201131109
37. Iikuni N, Lourenço EV, Hahn BH, La Cava A. Cutting edge: Regulatory T cells directly suppress B cells in systemic lupus erythematosus. J Immunol (Baltimore Md. (2009) 1950) 183:1518–22. doi: 10.4049/jimmunol.0901163
38. Mirlekar B. Tumor promoting roles of IL-10, TGF-β, IL-4, and IL-35: Its implications in cancer immunotherapy. SAGE Open Med. (2022) 10:20503121211069012. doi: 10.1177/20503121211069012
39. Mikami N, Kawakami R, Sakaguchi S. New Treg cell-based therapies of autoimmune diseases: towards antigen-specific immune suppression. Curr Opin Immunol. (2020) 67:36–41. doi: 10.1016/j.coi.2020.07.004
40. Wei WZ, Jacob JB, Zielinski JF, Flynn JC, Shim KD, Alsharabi G, et al. Concurrent induction of antitumor immunity and autoimmune thyroiditis in CD4+ CD25+ regulatory T cell-depleted mice. Cancer Res. (2005) 65:8471–8. doi: 10.1158/0008-5472.Can-05-0934
41. Skrivankova VW, Richmond RC, Woolf BAR, Yarmolinsky J, Davies NM, Swanson SA, et al. Strengthening the reporting of observational studies in epidemiology using Mendelian randomization: the STROBE-MR statement. Jama. (2021) 326:1614–21. doi: 10.1001/jama.2021.18236
42. Davies NM, Holmes MV, Davey Smith G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ (Clinical Res ed). (2018) 362:k601. doi: 10.1136/bmj.k601
43. Miyara M, Yoshioka Y, Kitoh A, Shima T, Wing K, Niwa A, et al. Functional delineation and differentiation dynamics of human CD4+ T cells expressing the FoxP3 transcription factor. Immunity. (2009) 30:899–911. doi: 10.1016/j.immuni.2009.03.019
44. Truffault F, Nazzal D, Verdier J, Gradolatto A, Fadel E, Roussin R, et al. Comparative analysis of thymic and blood Treg in Myasthenia gravis: thymic epithelial cells contribute to thymic immunoregulatory defects. Front Immunol. (2020) 11:782. doi: 10.3389/fimmu.2020.00782
45. Verma ND, Lam AD, Chiu C, Tran GT, Hall BM, Hodgkinson SJ. Multiple sclerosis patients have reduced resting and increased activated CD4(+)CD25(+)FOXP3(+)T regulatory cells. Sci Rep. (2021) 11:10476. doi: 10.1038/s41598-021-88448-5
46. Glatzová D, Cebecauer M. Dual role of CD4 in peripheral T lymphocytes. Front Immunol. (2019) 10:618. doi: 10.3389/fimmu.2019.00618
47. König R, Huang LY, Germain RN. MHC class II interaction with CD4 mediated by a region analogous to the MHC class I binding site for CD8. Nature. (1992) 356:796–8. doi: 10.1038/356796a0
48. Lu JG, Ji P, French SW. The major histocompatibility complex class II-CD4 immunologic synapse in alcoholic hepatitis and autoimmune liver pathology: the role of aberrant major histocompatibility complex class II in hepatocytes. Am J Pathol. (2020) 190:25–32. doi: 10.1016/j.ajpath.2019.09.019
49. Johnson DK, Magoffin W, Myers SJ, Finnell JG, Hancock JC, Orton TS, et al. CD4 inhibits helper T cell activation at lower affinity threshold for full-length T cell receptors than single chain signaling constructs. Front Immunol. (2020) 11:561889. doi: 10.3389/fimmu.2020.561889
50. Yongjun C, Nan Q, Yumeng S, Xiaowen J, Weibo W. Dioscin alleviates hashimoto's thyroiditis by regulating the SUMOylation of IRF4 to promote CD4(+)CD25(+)Foxp3(+) treg cell differentiation. Autoimmunity. (2021) 54:51–9. doi: 10.1080/08916934.2020.1855428
51. Chen A, Huang L, Zhang L. Helper T cell 17 and regulatory T cell levels in peripheral blood of newly diagnosed patients with autoimmune thyroid disease: A meta-analysis. Hormone Metab Res Hormon- und Stoffwechselforschung Hormones Metabolisme. (2023) 55:40–50. doi: 10.1055/a-1972-5787
52. Chen Z, Wang Y, Ding X, Zhang M, He M, Zhao Y, et al. The proportion of peripheral blood Tregs among the CD4+ T cells of autoimmune thyroid disease patients: a meta-analysis. Endocr J. (2020) 67:317–26. doi: 10.1507/endocrj.EJ19-0307
53. Xia S, Chen Q, Niu B. CD28: A new drug target for immune disease. Curr Drug Targets. (2020) 21:589–98. doi: 10.2174/1389450120666191114102830
54. Azuma M. Co-signal molecules in T-cell activation : historical overview and perspective. Adv Exp Med Biol. (2019) 1189:3–23. doi: 10.1007/978-981-32-9717-3_1
55. Chen L, Flies DB. Molecular mechanisms of T cell co-stimulation and co-inhibition. Nature reviews. Immunology. (2013) 13:227–42. doi: 10.1038/nri3405
56. Porciello N, Kunkl M, Tuosto L. CD28 between tolerance and autoimmunity: the side effects of animal models. F1000Research. (2018) 7:F1000 Faculty Rev-682. doi: 10.12688/f1000research.14046.1
57. Suntharalingam G, Perry MR, Ward S, Brett SJ, Castello-Cortes A, Brunner MD, et al. Cytokine storm in a phase 1 trial of the anti-CD28 monoclonal antibody TGN1412. New Engl J Med. (2006) 355:1018–28. doi: 10.1056/NEJMoa063842
58. Chenoweth AM, Esparon S, Wines BD, Schuurman J, Labrijn AF, Hogarth PM. Mutation of the TGN1412 anti-CD28 monoclonal antibody lower hinge confers specific FcγRIIb binding and retention of super-agonist activity. Immunol Cell Biol. (2023) 101:657–62. doi: 10.1111/imcb.12646
59. Panoskaltsis N, McCarthy NE, Stagg AJ, Mummery CJ, Husni M, Arebi N, et al. Immune reconstitution and clinical recovery following anti-CD28 antibody (TGN1412)-induced cytokine storm. Cancer Immunol Immunother CII. (2021) 70:1127–42. doi: 10.1007/s00262-020-02725-2
60. Lin CH, Hünig T. Efficient expansion of regulatory T cells in vitro and in vivo with a CD28 superagonist. Eur J Immunol. (2003) 33:626–38. doi: 10.1002/eji.200323570
61. Kyritsi EM, Kanaka-Gantenbein C. Autoimmune thyroid disease in specific genetic syndromes in childhood and adolescence. Front Endocrinol. (2020) 11:543. doi: 10.3389/fendo.2020.00543
62. Hossen MM, Ma Y, Yin Z, Xia Y, Du J, Huang JY, et al. Current understanding of CTLA-4: from mechanism to autoimmune diseases. Front Immunol. (2023) 14:1198365. doi: 10.3389/fimmu.2023.1198365
63. Liu W, Yang Z, Chen Y, Yang H, Wan X, Zhou X, et al. The association between CTLA-4, CD80/86, and CD28 gene polymorphisms and rheumatoid arthritis: an original study and meta-analysis. Front Med. (2021) 8:598076. doi: 10.3389/fmed.2021.598076
64. Van Coillie S, Wiernicki B, Xu J. Molecular and cellular functions of CTLA-4. Adv Exp Med Biol. (2020) 1248:7–32. doi: 10.1007/978-981-15-3266-5_2
65. Cutolo M, Sulli A, Paolino S, Pizzorni C. CTLA-4 blockade in the treatment of rheumatoid arthritis: an update. Expert Rev Clin Immunol. (2016) 12:417–25. doi: 10.1586/1744666x.2016.1133295
66. Huang L, Kong Y, Wang J, Sun J, Shi Q, Qiu YH. Reducing progression of experimental lupus nephritis via inhibition of the B7/CD28 signaling pathway. Mol Med Rep. (2015) 12:4187–95. doi: 10.3892/mmr.2015.3953
67. Hernandez R, LaPorte KM, Hsiung S, Santos Savio A, Malek TR. High-dose IL-2/CD25 fusion protein amplifies vaccine-induced CD4(+) and CD8(+) neoantigen-specific T cells to promote antitumor immunity. J Immunother Cancer. (2021) 9:e002865. doi: 10.1136/jitc-2021-002865
68. Li Y, Li X, Geng X, Zhao H. The IL-2A receptor pathway and its role in lymphocyte differentiation and function. Cytokine Growth Factor Rev. (2022) 67:66–79. doi: 10.1016/j.cytogfr.2022.06.004
69. Sakaguchi S, Sakaguchi N, Asano M, Itoh M, Toda M. Immunologic self-tolerance maintained by activated T cells expressing IL-2 receptor alpha-chains (CD25). Breakdown of a single mechanism of self-tolerance causes various autoimmune diseases. J Immunol (Baltimore Md. (1995) 1950) 155:1151–64. doi: 10.4049/jimmunol.155.3.1151
70. Fountoulakis S, Vartholomatos G, Kolaitis N, Frillingos S, Philippou G, Tsatsoulis A. HLA-DR expressing peripheral T regulatory cells in newly diagnosed patients with different forms of autoimmune thyroid disease. Thyroid. (2008) 18:1195–200. doi: 10.1089/thy.2008.0089
71. Morris GP, Brown NK, Kong YC. Naturally-existing CD4(+)CD25(+)Foxp3(+) regulatory T cells are required for tolerance to experimental autoimmune thyroiditis induced by either exogenous or endogenous autoantigen. J Autoimmun. (2009) 33:68–76. doi: 10.1016/j.jaut.2009.03.010
72. Malek TR, Yu A, Vincek V, Scibelli P, Kong L. CD4 regulatory T cells prevent lethal autoimmunity in IL-2Rbeta-deficient mice. Implications for the nonredundant function of IL-2. Immunity. (2002) 17:167–78. doi: 10.1016/s1074-7613(02)00367-9
Keywords: Hashimoto`s thyroiditis, Treg cells, immunity, causal inference, MR analysis
Citation: Guo J, Si G and Si F (2024) Treg cells as a protective factor for Hashimoto`s thyroiditis: a Mendelian randomization study. Front. Endocrinol. 15:1347695. doi: 10.3389/fendo.2024.1347695
Received: 01 December 2023; Accepted: 26 February 2024;
Published: 08 March 2024.
Edited by:
Artur Bossowski, Faculty of Medicine, Medical University of Bialystok, PolandReviewed by:
Ifigenia Kostoglou-Athanassiou, Asklepeion General Hospital, Athens, GreeceAnna Malgorzata Kucharska, Medical University of Warsaw, Poland
Copyright © 2024 Guo, Si and Si. 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: Fuchun Si, sifc2000@hotmail.com
†These authors have contributed equally to this work and share first authorship