Front. Endocrinol. Frontiers in Endocrinology Front. Endocrinol. 1664-2392 Frontiers Media S.A. 10.3389/fendo.2021.761077 Endocrinology Original Research Recurrence of Graves’ Disease: What Genetics of HLA and PTPN22 Can Tell Us Vejrazkova Daniela 1 * Vcelak Josef 1 Vaclavikova Eliska 1 Vankova Marketa 1 Zajickova Katerina 2 Vrbikova Jana 2 Duskova Michaela 3 Pacesova Petra 2 Novak Zdenek 2 Bendlova Bela 1 1 Department of Molecular Endocrinology, Institute of Endocrinology, Prague, Czechia 2 Department of Clinical Endocrinology, Institute of Endocrinology, Prague, Czechia 3 Department of Steroids and Proteohormones, Institute of Endocrinology, Prague, Czechia

Edited by: Yoshiyuki Ban, Teikyo University Chiba Medical Center, Japan

Reviewed by: Tim Cheetham, Newcastle University, United Kingdom; Filippo Biscarini, National Research Council (CNR), Italy

*Correspondence: Daniela Vejrazkova, dvejrazkova@endo.cz

This article was submitted to Thyroid Endocrinology, a section of the journal Frontiers in Endocrinology

 23 11 2021 2021 12 761077 20 08 2021 01 11 2021 Copyright © 2021 Vejrazkova, Vcelak, Vaclavikova, Vankova, Zajickova, Vrbikova, Duskova, Pacesova, Novak and Bendlova 2021 Vejrazkova, Vcelak, Vaclavikova, Vankova, Zajickova, Vrbikova, Duskova, Pacesova, Novak and Bendlova

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.

 Background

Approximately half of patients diagnosed with Graves’ disease (GD) relapse within two years of thyreostatic drug withdrawal. It is then necessary to decide whether to reintroduce conservative treatment that can have serious side effects, or to choose a radical approach. Familial forms of GD indicate a significant genetic component. Our aim was to evaluate the practical benefits of HLA and PTPN22 genetic testing for the assessment of disease recurrence risk in the Czech population.

 Methods

In 206 patients with GD, exon 2 in the HLA genes DRB1, DQA1, DQB1 and rs2476601 in the gene PTPN22 were sequenced.

 Results

The risk HLA haplotype DRB1*03-DQA1*05-DQB1*02 was more frequent in our GD patients than in the general European population. During long-term retrospective follow-up (many-year to lifelong perspective), 87 patients relapsed and 26 achieved remission lasting over 2 years indicating a 23% success rate for conservative treatment of the disease. In 93 people, the success of conservative treatment could not be evaluated (thyroidectomy immediately after the first attack or ongoing antithyroid therapy). Of the examined genes, the HLA-DQA1*05 variant reached statistical significance in terms of the ability to predict relapse (p=0.03). Combinations with either both other HLA risk genes forming the risk haplotype DRB1*03-DQA1*05-DQB1*02 or with the PTPN22 SNP did not improve the predictive value.

 Conclusion

the DQA1*05 variant may be a useful prognostic marker in patients with an unclear choice of treatment strategy.

 Graves’ disease HLA variants PTPN22 gene genetic predictors treatment Introduction

Graves’ disease (GD) is the most common cause of hyperthyroidism, affecting approximately 0.5% of men and 3% of women. In Europe, the first choice treatment is most often the administration of antithyroid drugs (1, 2). However, approximately half of patients relapse within two years of drug withdrawal. It is then necessary to decide whether to reintroduce thyreostatic therapy that may have serious side effects, or to choose a radical approach - total thyroidectomy (TTE) or radioiodine treatment.

A significant genetic component is evident from the familial occurrence of the disease. Studies on twin pairs have shown that the contribution of genetic factors can be as high as 70-80% (3). The autoimmune nature of GD has been associated with the human leukocyte antigen complex (HLA) as well as with the gene PTPN22 (protein tyrosine phosphatase, nonreceptor type 22) on chromosome 1p13.3-13.1 encoding for protein tyrosine phosphatase-22, a powerful inhibitor of T-cell activation. Within the HLA, DRB1*03, DQA1*05, and DQB1*02 allelic groups have proven to be promising predictors of the development and recurrence of GD in some studies (4). Outside the HLA region, the PTPN22 genetic variant rs2476601 has also been shown to be a potential Graves’ disease predictor and together with the HLA variants was included in the Graves’ Events After Therapy + (GREAT+) score (5). When individuals with a sensitive genetic background are exposed to certain environmental risk factors such as stress (6–8), smoking (9), iodine overdose (10), the postpartum period in women (11), microbiome-associated immunological changes (12) and possibly their combination, the production of autoantigens against the TSH receptor is triggered and the disease begins to develop or relapse. In addition to hyperthyroidism, extrathyroidal manifestations like Graves’ orbitopathy, thyroid dermatopathy, and rarely acropachy may be present. The long-term conservative treatment of persistent and recurrent hyperthyroidism entails considerable medical expenses and may have serious side effects in some cases, notably liver disorders and agranulocytosis.

In this retrospective study, we followed up on a pilot study published three years ago (13) with the aim to evaluate practical benefits of HLA and PTPN22 genetic testing for assessments of the disease recurrence risk in the Czech population. Taking genetics into account as part of long-term follow-up would make it easier for physicians and patients to consider the suitability and optimal timing of radical and definitive approaches to treatment.

 Materials and Methods Participants

We analyzed the HLA-DRB1, HLA-DQA1 and HLA-DQB1 allelic groups as well as the PTPN22 polymorphism rs2476601 (referred to also as 1858C/T) in 206 retrospectively observed Czech patients who had been diagnosed with Graves’ disease according to hormonal profile [low/suppressed TSH with simultaneously elevated freeT4, freeT3, presence of thyrotropin receptor antibodies (TRAK)] as well as sonographic examination of the thyroid gland. Orbit ultrasound had been performed in patients with present and active thyroid eye disease as a complementary exam while considering corticosteroids. Patients were recruited during the years 2017-2020 mainly from the outpatients of the Institute of Endocrinology in Prague. These patients are continuously monitored at the Institute, which gives us detailed long-term or even lifelong retrospective information on the status of their remission/relapse. An additional 16 patients were recruited from the outpatient Clinic of Endocrinology and Diabetology, M-centrum, Chotěboř, 3 patients were treated in the outpatient Clinic of Endocrinology in Říčany, and 3 patients were from the 3rd Medical Department, 1st Faculty of Medicine, Charles University and General Faculty Hospital in Prague. All participants gave their informed consent for inclusion before they participated in the study. The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of the Institute of Endocrinology (MH CR-RVO 00023761, date of approval 26. 6. 2017).

There were 171 women and 35 men in our cohort, and the mean age of patients at the time of diagnosis was 42.0 ± 14.54 years with median being 41 years and the age range spanning from 16 to 78 years. The mean duration of thyreostatic administration was 36.3 ± 42.47 months. Methimazol, or in case of intolerance propylthiouracyl, were prescribed. Systemic treatment of orbitopathy by peroral prednisone and/or by intravenous methylprednisolone pulses was indicated in 22.9% of the patients. Details regarding sex differences at the time of diagnosis, or in chronically ill patients at the time of the last relapse of GD, as well as retrospective specification of treatment are given in Table 1 .

Clinical characterization of the whole cohort of GD patients.

Characterization before treatment Women (N = 171) Men (N = 35) p-level
Age (years) 42.0 ± 14.44 43.9 ± 15.20 0.71*
Volume of goiter (ml) 20.9 ± 12.10 29.1 ± 17.31 0.01*
TSH (mIU/l) - normal range 0.27 to 4.2 0.018 ± 0.0977 0.011 ± 0.0153 0.62*
fT4 (pmol/l) - normal range 12 to 22 48.4 ± 23.09 46.9 ± 24.84 0.66*
fT3 (pmol/l) - normal range 3.1 to 6.8 18.5 ± 12.70 21.2 ± 13.86 0.39*
Way of treatment
Thyreostatic administration (months) 37.8 ± 44.79 27.9 ± 25.54 0.28*
Systemic treatment of orbitopathy (%) 23.9 18.2 0.47†
Thyroid surgery (%) 36.8 36.7 0.99†
RI administration (%) 10.43 13.8 0.58†

Data are given as mean ± SD or as a percentage; p-level according to *Mann-Whitney test or according to †Chi-square test.

Statistically significant p-levels are in bold.

 Genetic Analysis

Next generation sequencing performed on a MiSeq (Illumina) was used for haplotyping. In the three HLA genes (DRB1, DQA1 and DQB1), exon 2 was amplified, since this part of the HLA molecule determines its antigenic properties. Amplification of DNA sequences was performed using certified kits (CE 0197, CE IVD) of allele-specific pre-coated PCR strips PROTRANS S4 HLA-DRB1 * (cat. No. 34 04), PROTRANS S3 HLA-DQB1 * (cat. No. 33 09) and PROTRANS S2 HLA-DQA1 * (cat. No. 32 08) from the manufacturer PentaGen, Protrans, Germany. A public IPD-IMGT/HLA (International ImMunoGeneTics Information System/HLA database) alignment tool (https://www.ebi.ac.uk/ipd/imgt/hla/align.html) (14) was used to identify the allelic groups after sequencing. For the DRB1 and DQB1 HLA allelic groups, the Allele Frequencies website (http://www.allelefrequencies.net/) (15) allowed a comparison of the obtained allele frequencies of GD patients with the Czech bone marrow donors (n=5099).

The polymorphism rs2476601; c.1858C> T p. (Arg620Trp) was analyzed by sequencing exon 14 of the PTPN22 gene (protein tyrosine phosphatase, nonreceptor type 22) by next generation sequencing on a MiSeq sequencer (Illumina) using the Nextera XT kit (Illumina) to prepare libraries.

 Biochemical Analysis

Serum TSH, fT4 and fT3 concentrations were measured using the ECLIA method (Cobas Integra, Roche Diagnostics). For detection of autoantibodies to the TSH receptor, a fully automated electrochemiluminescence immunoassay Elecsys Anti-TSHR (TRAK) was used (Elecsys 2010, Cobas e601, Roche). The precision of this method is less than 6%, sensitivity 97% and specificity 99,5%.

 Statistical Analysis

Statistical analyses were performed using the NCSS/Pass 2004 software (NCSS, LLC, Kaysville, Utah, USA). Data are presented as means ± SD or as a percentage. The Chi-square test was used to compare the distribution of allelic groups in the particular cohorts. Odds ratios, relative risk and 95% confidence intervals were calculated in MedCalc Software. Differences in anamnestic data between men and women were tested by the non-parametric Mann-Whitney test. All tests were two-tailed (both positive and negative differences were considered). The p-value<0.05 threshold was used to suggest statistical significant differences.

 Results Evaluation of the Success of Conservative Treatment

As regards the distinction between successfully and unsuccessfully treated patients, we followed these criteria: patients remaining in remission at least 2 years since discontinuation of thyreostatic treatment and with no relapse throughout life were considered successfully treated. On the contrary, those who had relapsed at least once during monitoring at our institution (many-year to lifelong perspective) as well as patients requiring thyreostatics for more than 5 years were considered unsuccessfully treated. Of the 206 patients examined, 87 met the criteria for unsuccessful treatment, 26 patients met the criteria for successful treatment, indicating a 23% success rate for conservative treatment of the disease in long term follow-up. In 93 people, the success of treatment could not be evaluated because they either underwent TTE immediately after the first attack (so it was not possible to determine whether the disease would relapse), they are still on antithyroid therapy (for a period of less than 5 years, so they cannot yet be classified as unsuccessfully treated), or their remission has not yet reached two years (so they cannot yet be classified as successfully treated).

Comparison of patients with successful and unsuccessful conservative treatment is given in Table 2 . Significant difference was observed in fT3 levels. Unsuccessfully treated patients showed twice as high concentrations compared to patients in remission. For fT4 levels, such differences between the two groups were not apparent.

Comparison of patients with successful and unsuccessful conservative treatment.

Characterization before treatment Successful treatment (N = 26) Unsuccessful treatment (N = 87) p-level
Age (years) 41.0 ± 13.05 44.9 ± 14.23 0.17*
Volume of goiter (ml) 17.8 ± 7.82 22.0 ± 13.43 0.27*
TSH (mIU/l) - normal range 0.27 to 4.2 0.011 ± 0.0160 0.013 ± 0.0198 0.80*
fT4 (pmol/l) - normal range 12 to 22 43.6 ± 19.12 46.6 ± 24.27 1.00*
fT3 (pmol/l) - normal range 3.1 to 6.8 10.4 ± 4.78 21.5 ± 15.37 0.03*
Way of treatment
Thyreostatic administration (months) 25.1 ± 17.02 56.1 ± 55.11 0.002*
Systemic treatment of orbitopathy (%) 9.1 24.1 12.4†
Thyroid surgery (%) 0 31.7 –
RI administration (%) 0 12.0 –

Data are given as mean ± SD or as a percentage; p-level according to *Mann-Whitney test or according to †Chi-square test.

Statistically significant p-levels are in bold.

 Frequencies of the HLA Allelic Groups

The risk haplotype DRB1*03-DQA1*05-DQB1*02 was present in 35% of the entire cohort of 206 GD patients, 5 patients were homozygotes in this haplotype, 67 were heterozygotes. Haplotype frequency in the whole cohort was 18.7%. This highly exceeds the frequencies in non-GD populations as given in the Allele Frequencies website. For the individual HLA genes and the proportion of allelic groups in the whole cohort of GD patients, see Table 3 .

Frequencies of the HLA allelic groups in the whole cohort of GD patients, N = 206.

HLA genes Allelic groups of individual HLA genes
*01 *02 *03 *04 *05 *06 *07 *08 *09 *10 *11 *12 *13 *14 *15 *16
DRB1 9 – 21 10 – – 6 4 1 0.5 16 1.5 12 2 13 4
DQA1 37 6 11 3 42 1 – – – – – – – – – –
DQB1 – 26 34 3 16 21 – – – – – – – – – –

Data are given as %, frequencies of the risk variants of the particular genes are in bold.

The asterix represent an established form of HLA genes nomenclature in accordance with the current HLA nomenclature (http://hla.alleles.org/nomenclature/nomenc_updates.html).

For the frequencies of the HLA allelic groups analyzed separately in the unsuccessfully and successfully treated cohort as well as in the patients whose success of treatment could not be evaluated, see Tables 4 – 6 , respectively.

Frequencies of the HLA allelic groups in unsuccessfully treated GD patients, N = 87.

HLA genes Allelic groups of individual HLA genes
*01 *02 *03 *04 *05 *06 *07 *08 *09 *10 *11 *12 *13 *14 *15 *16
DRB1 10 – 23 9 – – 3 4 0.5 0.5 18 2 15 3 10 2
DQA1 35 3 10 3 48 1 – – – – – – – – – –
DQB1 – 25 37 3 15 20 – – – – – – – – – –

Data are given as %, frequencies of the risk variants of the particular genes are in bold.

The asterix represent an established form of HLA genes nomenclature in accordance with the current HLA nomenclature (http://hla.alleles.org/nomenclature/nomenc_updates.html).

Frequencies of the HLA allelic groups in successfully treated GD patients, N = 26.

HLA genes Allelic groups of individual HLA genes
*01 *02 *03 *04 *05 *06 *07 *08 *09 *10 *11 *12 *13 *14 *15 *16
DRB1 10 – 18 10 – – 10 4 0 2 10 0 10 2 18 6
DQA1 46 10 10 4 30 0 – – – – – – – – – –
DQB1 – 26 22 6 20 26 – – – – – – – – – –

Data are given as %, frequencies of the risk variants of the particular genes are in bold.

The asterix represent an established form of HLA genes nomenclature in accordance with the current HLA nomenclature (http://hla.alleles.org/nomenclature/nomenc_updates.html).

Frequencies of the HLA allelic groups in GD patients whose success of treatment could not be evaluated, N = 93.

HLA genes Allelic groups of individual HLA genes
*01 *02 *03 *04 *05 *06 *07 *08 *09 *10 *11 *12 *13 *14 *15 *16
DRB1 9 – 19.5 11 – – 8 3 1 0 15.5 1.5 10.5 1 14 6
DQA1 36 8.5 12.5 2.5 40 0.5 – – – – – – – – – –
DQB1 – 26.5 35 2.5 15.5 20.5 – – – – – – – – – –

Data are given as %, frequencies of the risk variants of the particular genes are in bold.

The asterix represent an established form of HLA genes nomenclature in accordance with the current HLA nomenclature (http://hla.alleles.org/nomenclature/nomenc_updates.html).

Evaluated in the overall cohort of 206 participants, the allele DRB1*03 was more frequent in GD patients (21%) compared to the Czech bone marrow donors (12%) we used as a reference, see details in Methods. On the contrary, the DRB1*07 allele was less frequent in our GD group (6%) in comparison with bone marrow donors (14%), indicating its protective effect in terms of disease development. It is noteworthy that, in the unsuccessfully treated subgroup, the DRB1*07 allele frequency reached only 3%, while in the successfully treated subgroup it reached 10% (p=0.05), indicating a protective effect of the *07 variant even in terms of disease relapse. The risk variant DQA1*05 was present in 67.0% of our GD patients with an allele frequency of 42%, 17.5% of the patients (36 individuals) were homozygotes in this allelic group. The proportion of these homozygotes reached 18.1% in unsuccessfully treated patients, whereas it was only 12.0% in successfully treated patients (p=0.035). Unfortunately, for the HLA-DQA1 gene, a comparison with the Czech bone marrow donors was not available; however, according to the Allele Frequencies website, the allele frequency observed in our study in Czech GD patients was markedly higher than in the other listed European populations, see the Discussion for details. The frequency of the risk variant DQB1*02 observed in our group of patients (26%) was comparable to the frequency in Czech bone marrow donors (23%). No significant differences in the distribution of the particular risk allelic groups DRB1*03, DQA1*05 and DQB1*02 or of the risk haplotype DRB1*03-DQA1*05-DQB1*02 were observed between men and women.

Complete allele, phenotype and haplotype frequencies of our cohort of 206 GD patients were published in public repository Figshare.com (16) with DOI for public link: 10.6084/m9.figshare.16446954. The raw sequencing data are available in the public repository NCBI (17), accession ID to cite for the Sequence Read Archive (SRA) data: PRJNA776790, public link to access the SRA records: https://www.ncbi.nlm.nih.gov/sra/PRJNA776790.

 Success of Conservative Treatment in Relation to Genotype

An evaluation of the three HLA genes and the PTPN22 SNP in terms of their association with treatment success showed that only the HLA-DQA1*05 variant reached statistical significance (p=0.03), as detailed in Table 7 :

Distribution of the studied genetic variants in the group of GD patients depending on the success of the treatment, N = 113.

DRB1*03 presence Successful conservative treatment Statistics
no yes Chi2-test=0.59
no 53 18 power=0.12
yes 34 8 p=0.44
DQB1*02 presence Successful conservative treatment Statistics
no yes Chi2-test=0.0003
no 47 14 power=0.05
yes 40 12 p=0.99
DQA1*05 presence Successful conservative treatment Statistics
no yes Chi2-test=4.69
no 21 12 power=0.6
yes 66 14 p=0.03
DRB1*03-DQB1*02-DQA1*05 haplotype Successful conservative treatment Statistics
no yes Chi2-test=0.44
no 54 18 power=0.10
yes 33 8 p=0.51
PTPN22 variant Successful conservative treatment Statistics
no yes Chi2-test=1.62
no 59 21 power=0.25
yes 28 5 p=0.20

Statistically significant p-levels are in bold.

The asterix represent an established form of HLA genes nomenclature in accordance with the current HLA nomenclature (http://hla.alleles.org/nomenclature/nomenc_updates.html).

There were 14 successfully treated patients among the 80 DQA1*05 carriers (17.5%), while there were 66 DQA1*05 carriers among the 87 unsuccessfully treated patients (75.9%). The combination of the DQA1*05 variant with both other risk HLA genetic variants forming the risk haplotype DRB1*03-DQA1*05-DQB1*02 and with the PTPN22 SNP did not improve the ability to predict relapse. The odds ratio (OR) for relapse (or more precisely, for unsuccessful conservative treatment, that is, for relapse of the disease or for prolonged inconclusive treatment lasting more than 5 years, which we observed in several cases) calculated as the ratio of the odds of the unsuccessful treatment occurring in the group of risk allele DQA1*05 carriers to the odds of it occurring in the group of risk allele DQA1*05 non-carriers was 2.7 (95%CI: 1.08-6.72); p=0.05. The relative risk was 1.3 (95%CI: 0.98-1.71); p=0.07.

The DRB1*03 variant was not associated with treatment outcome (p=0.44), nor was the allelic group DQB1*02 (p=0.73), see Table 7 .

As regards rs2476601 (1858C/T) in the PTPN22 gene, the minor allele was present in 27% of patients, with an allelic frequency of 14.1%. The allelic frequency in unsuccessfully treated GD patients was 16.7%, in successfully treated patients 12.0% (Chi2-test=0.45; p=0.50). In patients with impossible evaluation of the treatment success, the allelic frequency was 12.4%. The OR calculated for relapse was 1.99, (95%CI: 0.68-5.83); p=0.30, and the relative risk was 1.15 (95%CI: 0.95-1.40); p=0.16. Although this polymorphism causing the substitution of arginine to tryptophan at codon 620 was found to be associated with GD in many European studies, its association with disease relapse was not apparent in our data.

 Discussion

Studies in the past two decades have clearly shown that multiple factors are involved in the development and recurrence of GD. The interaction of susceptible genes with variable environmental factors, mediated through very complex endogenous communication including dynamic epigenetic modulation and gene expression changes, may lead to impaired immunological tolerance and the outbreak of the disease (18–20). In addition to HLA variants, many non-HLA genes have been confirmed to be involved in the etiology of GD; some are unique to the recurrence (21, 22) of GD, while others are common to both the disease onset and recurrence, or pose a risk for the development of an even wider range of autoimmune diseases (23).

According to our data, the risk HLA haplotype DRB1*03-DQA1*05-DQB1*02 was more frequent in GD patients than in the general European population, as given in the Allele Frequencies website, with the two exceptions of Sardinia and Slovenia, that are represented by small groups (n ≤ 140). For the relatively small Czech control population sample in this website (n=180), a frequency of 9.2% is given for the haplotype DRB1*03-DQA1*05:01-DQB1*02:01. The haplotype frequency observed in our cohort of patients (18.7%) corresponds with a metaanalysis on GD patients by Gough (4). The risk allelic group DRB1*03 was more frequent in the GD group compared to the Czech bone marrow donors, which is in agreement with previous findings in other populations (24). On the other hand, the DRB1*07 allele was less frequent in our GD group in comparison with bone marrow donors and its frequency was clearly the lowest in unsuccessfully treated patients. The DRB1*07 allele has been reported to be protective for GD in UK Caucasians (25); however, due to its low incidence, a larger cohort of patients would be needed to demonstrate its significant protective role in the Czech population.

The allele frequency of the risk variant DQA1*05 was markedly higher in our study (42%) than in other European populations, ranging from 16.6% in France to 30.3% in Belgium according to the Allele Frequencies website. The DQA1*05 variant showed the best, though borderline, predictive ability in terms of disease relapse, with an OR=2.7. Although it cannot be considered a reliable recurrence predictor, its evaluation may be useful in patients with a complex clinical picture and unclear treatment strategy.

The tyrosine phosphatase-22 protein encoded by the PTPN22 gene inhibits T-cell activation. The substitution of arginine to tryptophan at codon 620 (rs2476601) disrupts an interaction motif in the protein (26). This SNP has been associated with GD in many studies. These associations have shown high reproducibility in Caucasian populations, with an OR of 1.9 in British Caucasians (27), 1.7 in a Polish population (28), and even 4.2 in a Russian population (29); however, in Asian and African populations the minor allele was sporadic or absent (30). The concept of our study does not allow the calculation of OR in this sense, because we worked only with a group of patients and have no comparison with the Czech healthy controls for this SNP. However, the OR calculated for relapse did not show statistical significance. A comparison with some reference European populations suggests that there is higher minor allele frequency in our GD patients (14.1%). For example, the minor allele frequency in controls was reported as 10.4% in the United Kingdom (31), 10.0% in Germany (32) and 9% in Italy (33). However, the northernmost and easternmost countries of the continent show minor allele frequencies in the general population similar to our GD cohort (15% in Finland and 14.1% in the Ukraine (34). Either way, our study did not demonstrate the predictive potential of the PTPN22 variant in terms of disease recurrence.

It is worth noting the observation of significantly higher fT3 values - but not fT4 values - in patients whose conservative treatment has not been successful for a long time. It is consistent with the conclusions of the original (35) and review articles (36, 37) on this topic, which describe a higher fT3 to fT4 ratio in patients who failed to respond to antithyroid drug therapy. This association may be mediated through variability in deiodinase type 2 gene, which affects activity of type 2 deiodinase regulating the transformation of T4 into T3 (38).

The absence of a control group of healthy Czech volunteers, ideally without history of autoimmune disorders, is one limitation for an interpretation of the data presented here. Furthermore, the size of the cohort suitable for assessing the predictive ability in terms of disease relapse should be larger, considering that only 113 people could be included in the statistical evaluation. On the other hand, an advantage of our study is the availability of medical history data and retrospective course of the disease of almost all patients in the sample, as they are persons receiving long-term and often life-long care in our institution. This approach indicates a 23% success rate for conservative treatment. Available medical history will also allow for a re-evaluation after several years, when it will be interesting to verify how many participants now belonging to the group of successfully treated patients in remission will move to the group of those who relapsed, and how these potential changes will affect the statistical results.

In conclusion, we analyzed the HLA DRB1, DQA1 and DQB1 allelic groups and rs2476601 in the PTPN22 gene in 206 patients diagnosed with GD. According to our data, the proportions of the risk variant HLA-DQA1*05 as well as the risk HLA haplotype DRB1*03-DQA1*05-DQB1*02 were significantly higher in GD patients than in the general European population as given in the Allele Frequencies website. Although the tested genes individually or in combination cannot be considered reliable relapse predictors in the Czech population, the HLA-DQA1*05 allelic group showed a statistically significant association with relapse or long-term unsuccessful conservative treatment. Therefore, this variant may be a useful prognostic marker in GD patients with a difficult interpretation of their clinical picture, when treatment strategies remain unclear.

 Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI SRA Bio Project, accession no: PRJNA776790.

 Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of the Institute of Endocrinology. The patients/participants provided their written informed consent to participate in this study.

 Author Contributions

Conceptualization and design of the work, formal analysis, project administration, original draft preparation and writing, DV. Methodology of genetic analyzes, JVc and EV. Review & editing, statistics, MV. Data curation, critical revision of the work, KZ, JVr, MD, PP, and ZN. Supervision, BB. All authors contributed to the article and approved the submitted version.

 Funding

The study was supported by Ministry of Health of the Czech Republic, grant RVO 00023761.

 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.

 Acknowledgments

The authors would like to thank Mgr. Nela Vondrková and to all physicians of the Institute of Endocrinology for their excellent assistance. Thanks also belong to MUDr. Mgr. Sylvie Špitálníková, Ph.D. from the Clinic of Endocrinology and Diabetology in Chotěboř, to doc. Zdeňka Límanová, CSc. from the 1st Faculty of Medicine, Charles University and General Faculty Hospital in Prague, and to MUDr. Magdalena Matějková Běhanová, Ph.D. from the Clinic of Endocrinology in Říčany.

 References Ross DS Burch HB Cooper DS Greenlee MC Laurberg P Maia AL . 2016 American Thyroid Association Guidelines for Diagnosis and Management of Hyperthyroidism and Other Causes of Thyrotoxicosis. Thyroid (2016) 26:1343–421. doi: 10.1089/thy.2016.0229 Kahaly GJ Bartalena L Hegedüs L Leenhardt L Poppe K Pearce SH . 2018 European Thyroid Association Guideline for the Management of Graves’ Hyperthyroidism. Eur Thyroid J (2018) 7:167–86. doi: 10.1159/000490384 Brix TH Kyvik KO Christensen K Hegedüs L . Evidence for a Major Role of Heredity in Graves’ Disease: A Population-Based Study of Two Danish Twin Cohorts. J Clin Endocrinol Metab (2001) 86:930–4. doi: 10.1210/jcem.86.2.7242 Gough SC . The Genetics of Graves’ Disease. Endocrinol Metab Clin North Am (2000) 29:255–66. doi: 10.1016/s0889-8529(05)70130-4 Vos XG Endert E Zwinderman AH Tijssen JG Wiersinga WM . Predicting the Risk of Recurrence Before the Start of Antithyroid Drug Therapy in Patients With Graves’ Hyperthyroidism. J Clin Endocrinol Metab (2016) 101:1381–9. doi: 10.1210/jc.2015-3644 Effraimidis G Tijssen JG Brosschot JF Wiersinga WM . Involvement of Stress in the Pathogenesis of Autoimmune Thyroid Disease: A Prospective Study. Psychoneuroendocrinology (2012) 37:1191–8. doi: 10.1016/j.psyneuen.2011.12.009 Vita R Lapa D Trimarchi F Benvenga S . Stress Triggers the Onset and the Recurrences of Hyperthyroidism in Patients With Graves’ Disease. Endocrine (2015) 48:254–63. doi: 10.1007/s12020-014-0289-8 Falgarone G Heshmati HM Cohen R Reach G . Mechanisms in Endocrinology. Role of Emotional Stress in the Pathophysiology of Graves’ Disease. Eur J Endocrinol (2012) 168:R13–8. doi: 10.1530/EJE-12-0539 Effraimidis G Wiersinga WM . Mechanisms in Endocrinology: Autoimmune Thyroid Disease: Old and New Players. Eur J Endocrinol (2014) 170:R241–52. doi: 10.1530/EJE-14-0047 Pedersen IB Knudsen N Carlé A Vejbjerg P Jørgensen T Perrild H . A Cautious Iodization Programme Bringing Iodine Intake to a Low Recommended Level is Associated With an Increase in the Prevalence of Thyroid Autoantibodies in the Population. Clin Endocrinol (Oxf) (2011) 75:120–6. doi: 10.1111/j.1365-2265.2011.04008.x Alexander EK Pearce EN Brent GA Brown RS Chen H Dosiou C . 2017 Guidelines of the American Thyroid Association for the Diagnosis and Management of Thyroid Disease During Pregnancy and the Postpartum. Thyroid (2017) 27:315–89. doi: 10.1089/thy.2016.0457 Astl J Šterzl I . Activation of Helicobacter Pylori Causes Either Autoimmune Thyroid Diseases or Carcinogenesis in the Digestive Tract. Physiol Res (2015) 64:S291–301. doi: 10.33549/physiolres.933118 Vejrazkova D Vcelak J Vaclavikova E Vankova M Zajickova K Duskova M . Genetic Predictors of the Development and Recurrence of Graves’ Disease. Physiol Res (2018) 67:S431–9. doi: 10.33549/physiolres.934018 The IPD-IMGT/HLA Alignment Tool. Available at: https://www.ebi.ac.uk/ipd/imgt/hla/align.html (Accessed on 30. 11. 2020). The Allele Frequencies Webside. Available at: http://www.allelefrequencies.net/ (Accessed on 30. 11. 2020). Figshare Data Repository. Available at: https://figshare.com (Accessed on 16. 9. 2021). NCBI Sequence Read Archive (SRA) Data Repository. Available at: https://www.ncbi.nlm.nih.gov/sra (Accessed on 27. 10. 2021). Wang B Shao X Song R Xu D Zhang JA . The Emerging Role of Epigenetics in Autoimmune Thyroid Diseases. Front Immunol (2017) 8:396. doi: 10.3389/fimmu.2017.00396 Imgenberg-Kreuz J Carlsson Almlöf J Leonard D Alexsson A Nordmark G Eloranta ML . DNA Methylation Mapping Identifies Gene Regulatory Effects in Patients With Systemic Lupus Erythematosus. Ann Rheum Dis (2018) 77:736–43. doi: 10.1136/annrheumdis-2017-212379 Wiersinga WM . Graves’ Disease: Can It Be Cured? Endocrinol Metab (Seoul) (2019) 34:29–38. doi: 10.3803/EnM.2019.34.1.29 Liu J Fu J Duan Y Wang G . Predictive Value of Gene Polymorphisms on Recurrence After the Withdrawal of Antithyroid Drugs in Patients With Graves’ Disease. Front Endocrinol (Lausanne) (2017) 8:258. doi: 10.3389/fendo.2017.00258 Azizi F Amouzegar A Tohidi M Hedayati M Khalili D Cheraghi L . Increased Remission Rates After Long-Term Methimazole Therapy in Patients With Graves’ Disease: Results of a Randomized Clinical Trial. Thyroid (2019) 29:1192–200. doi: 10.1089/thy.2019.0180 Tomer Y . Mechanisms of Autoimmune Thyroid Diseases: From Genetics to Epigenetics. Annu Rev Pathol (2014) 9:147–56. doi: 10.1146/annurev-pathol-012513-104713 Tomer Y Davies TF . Searching for the Autoimmune Thyroid Disease Susceptibility Genes: From Gene Mapping to Gene Function. Endocr Rev (2003) 24:694–717. doi: 10.1210/er.2002-0030 Simmonds MJ Howson JM Heward JM Cordell HJ Foxall H Carr-Smith J . Regression Mapping of Association Between the Human Leukocyte Antigen Region and Graves Disease. Am J Hum Genet (2005) 76:157–63. doi: 10.1086/426947 Stanford SM Bottini N . PTPN22: The Archetypal non-HLA Autoimmunity Gene. Nat Rev Rheumatol (2014) 10:602–11. doi: 10.1038/nrrheum.2014.109 Velaga MR Wilson V Jennings CE Owen CJ Herington S Donaldson PT . The Codon 620 Tryptophan Allele of the Lymphoid Tyrosine Phosphatase (LYP) Gene is a Major Determinant of Graves’ Disease. J Clin Endocrinol Metab (2004) 89:5862–5. doi: 10.1210/jc.2004-1108 Skórka A Bednarczuk T Bar-Andziak E Nauman J Ploski R . Lymphoid Tyrosine Phosphatase (PTPN22/LYP) Variant and Graves’ Disease in a Polish Population: Association and Gene Dose-Dependent Correlation With Age of Onset. Clin Endocrinol (Oxf) (2005) 62:679–82. doi: 10.1111/j.1365-2265.2005.02279.x Zhebrun D Kudryashova Y Babenko A Maslyansky A Kunitskaya N Popcova D . Association of PTPN22 1858t/T Genotype With Type 1 Diabetes, Graves’ Disease But Not With Rheumatoid Arthritis in Russian Population. Aging (Albany NY) (2011) 3:368–73. doi: 10.18632/aging.100305 Nabi G Akhter N Wahid M Bhatia K Mandal RK Dar SA . Meta-Analysis Reveals PTPN22 1858c/T Polymorphism Confers Susceptibility to Rheumatoid Arthritis in Caucasian But Not in Asian Population. Autoimmunity (2016) 49:197–210. doi: 10.3109/08916934.2015.1134514 Smyth D Cooper JD Collins JE Heward JM Franklyn JA Howson JM . Replication of an Association Between the Lymphoid Tyrosine Phosphatase Locus (LYP/PTPN22) With Type 1 Diabetes, and Evidence for its Role as a General Autoimmunity Locus. Diabetes (2004) 53:3020–3. doi: 10.2337/diabetes.53.11.3020 Burkhardt H Hüffmeier U Spriewald B Böhm B Rau R Kallert S . Association Between Protein Tyrosine Phosphatase 22 Variant R620W in Conjunction With the HLA-DRB1 Shared Epitope and Humoral Autoimmunity to an Immunodominant Epitope of Cartilage-Specific Type II Collagen in Early Rheumatoid Arthritis. Arthritis Rheum (2006) 54:82–9. doi: 10.1002/art.21498 Totaro MC Tolusso B Napolioni V Faustini F Canestri S Mannocci A . PTPN22 1858c>T Polymorphism Distribution in Europe and Association With Rheumatoid Arthritis: Case-Control Study and Meta-Analysis. PloS One (2011) 6:e24292. doi: 10.1371/journal.pone.0024292 Burn GL Svensson L Sanchez-Blanco C Saini M Cope AP . Why is PTPN22 a Good Candidate Susceptibility Gene for Autoimmune Disease? FEBS Lett (2011) 585:3689–98. doi: 10.1016/j.febslet.2011.04.032 Miao J Zhao YJ Wang S Jiang XH Zhao ZF Gu LQ . Prognostic Factors in the Relapse of Graves Disease. Zhonghua Nei Ke Za Zhi (2008) 47:185–8. Liu J Fu J Xu Y Wang G . Antithyroid Drug Therapy for Graves’ Disease and Implications for Recurrence. Int J Endocrinol (2017) 2017:3813540. doi: 10.1155/2017/3813540 Shi H Sheng R Hu Y Liu X Jiang L Wang Z . Risk Factors for the Relapse of Graves’ Disease Treated With Antithyroid Drugs: A Systematic Review and Meta-Analysis. Clin Ther (2020) 42:662–675.e4. doi: 10.1016/j.clinthera.2020.01.022 Babenko A Popkova D Freylihman O Solncev V Kostareva A Grineva E . Thr92Ala Polymorphism of Human Type 2 Deiodinase Gene (Hd2) Affects the Development of Graves’ Disease, Treatment Efficiency, and Rate of Remission. Clin Dev Immunol (2012) 2012:340542. doi: 10.1155/2012/340542