Skip to main content

Hashimoto’s thyroiditis increases the risk of new-onset systemic lupus erythematosus: a nationwide population-based cohort study



Previous studies have shown systemic lupus erythematosus (SLE) patients had a significantly higher prevalence of thyroid diseases and hypothyroidism than matched controls, and some case reports showed SLE may occur after Hashimoto’s thyroiditis (HT).


This study aimed to investigate the subsequent risk of SLE in patients with HT.


In this retrospective cohort study done by the Taiwan National Health Insurance Research Database, the HT group (exposure group) and the non-HT group (comparator group) were propensity score matched at a ratio of 1:2 by demographic data, comorbidities, medications, and the index date. We used Cox proportional hazards models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs). Several sensitivity analyses were done for cross-validation of our findings.


We identified 15,512 HT patients and matched 31,024 individuals. The incidence rate ratio of SLE was 3.58 (95% CI, 2.43–5.28; p < 0.01). Several sensitivity analyses show adjusted hazard ratio (aHR) (CIs) of 4.35 (3.28–5.76), 4.39 (3.31–5.82), 5.11 (3.75–6.98), and 4.70 (3.46–6.38), consistent with the results of the main model.


Our study showed an increased risk of SLE in the HT group after adjustment for baseline characteristics, comorbidities, and medical confounders compared with the reference group.


Systemic lupus erythematosus (SLE) is a systematic autoimmune disease affecting almost every organ, in which the immune system attacks tissues and cells leading to inflammation and damage [1]. Previous studies have shown that both the loss of B cell self-tolerance due to genetic factor and Th1 lymphocyte response are necessary for the development of SLE [2,3,4]. However, some studies indicated that genetic susceptibility alone is not sufficient to account for all SLE patients, with only 24% concordance within monozygotic twins with SLE, which is much lower than previous estimate [5]. Another study revealed that environmental exposures are more related to SLE. In this study, a higher concordance of autoimmune disease, including SLE, was presented among monozygotic twins than among dizygotic twins, which was explained as environment factors that monozygotic twins usually share are more similar than dizygotic do [6], rather than just considered as a genetic similarity. The environment factors for SLE include but not limited to ultraviolet radiation, infection, and hormone [7]. Besides, having a history of other autoimmune diseases is also a risk factor for another one, according to some studies about multiple autoimmune syndrome (MAS) [8, 9].

Studies also revealed that T3 and T4 act as modulators in the immune system [10, 11]. Thyroid hormone not only stimulates T cells, B cells [12, 13], neutrophils, and macrophage chemotaxis [14, 15] but also enhances the generation of ROS and IL-18 [16]. Despite that abundant studies about T3 and T4 alternating the immune status, we now only noticed that SLE patients are more likely to develop autoimmune thyroiditis or even hypothyroidism in some cohort study [2, 17], and there is no definite evidence that make sure whether hypothyroidism will cause SLE so far in spite of two case reports about two young girls and two women respectively evolving SLE after being diagnosed as hypothyroidism and Hashimoto’s thyroiditis (HT) [18, 19].

Based on the concept mentioned above, we hypothesized that a history of HT increases the risk of subsequent SLE, and, as the epidemiological correlation between HT and SLE has remained unclear, we designed a retrospective cohort study to investigate this issue.

Materials and methods


Study design

In this national wide, retrospective cohort study with propensity score matched (PSM), the data was extracted from the National Health Insurance Research Database (NHIRD), which is a database constructed by 95% of residents in Taiwan since 1995 through the National health care insurance (NHI) system. Disease profiles are based on the International Classification of Diseases, 9th Revision, and Clinical Modification (ICD-9-CM) systems; the NHIRD provides information including demographics, outpatient visits, and hospitalizations with dates, prescriptions codes, diagnostic codes, laboratory tests and interventional procedure codes, and medical costs. De-identification was done for the protection of personal privacy.

We also extract the major illness registry data, also called the catastrophic illness registry, which is a certification for patients who were diagnosed with some severe, chronic, or fatal disease including cancer, diabetes mellitus, major injury, and SLE. The certification provides a discount on admission and medical charge. The Longitudinal Health Insurance Database 2000 (LHID 2000) was also used in the comparator group without HT of this study; LHID 2000 has collected a population of 1 million which was randomly sampled from the beneficiaries’ registration files within the year 2000. The representative of gender and age distribution has been statistically confirmed between the LHID 2000 and the origin NHIRD data.

Study population and propensity score match

We identified our population as patients who have had at least three outpatient visits or one hospital admission for autoimmune thyroiditis (Hashimoto’s thyroiditis) (ICD9 code: 245.2) between 2003 and 2012. The index date for the corresponding matching was defined as the date of the first diagnosis of autoimmune thyroiditis, outpatient, or admission.

We also constructed a comparator group without HT sampled from the LHID 2000 data, which includes patients ever visited outpatient departments between 2005 and 2012. The index date of the comparator group without HT was defined as the first visit to the outpatient department each year. Those who ever had at least one outpatient visit plus one hospitalization under the diagnosis of disorders of the thyroid gland (ICD9 codes 240–246: 240 for simple and unspecified goiter, 241 for nontoxic nodular goiter, 242 for thyrotoxicosis with or without goiter, 243 for congenital hypothyroidism, 244 for acquired hypothyroidism, 245 for thyroiditis, 246 for other disorders of the thyroid) between 1997 and 2013 were excluded.

We excluded the data of index date which were unmatched among the study and comparator group without HT (2003–2004), SLE (ICD9 code: 710.0) diagnosed before the index date and death during follow-up.

Overall, we extracted a total population of 17,978 cases with HT between the years 2005 and 2012 as our study group; and the comparator group without HT contained 783,345 patients without any disorders of the thyroid gland. To reduce the confounding bias, propensity score matching (PSM) was used, which was estimated by logistic regression modeling. Predictors involved index date, gender, and selected co-morbidities. The 1:2 matched comparator shares the same propensity score as the exposure group.

Outcomes and comorbidities

The primary outcome of the study was SLE occurrence, which was defined as patients who were diagnosed with SLE (ICD9 code: 710.0) and were identified as having “major illness” according to the NHI document for ensuring only correctly diagnosed patients were included. The follow-up started on the respective index date for different individuals until SLE was diagnosed or withdrawn from NHI due to any cause such as death, leaving, loss of data, or end of the study (December 2013), whichever occurred first. Relevant data of background variation including gender, age, urbanization, low income, length of hospital stays, times of outpatient department visits, medication control, and co-morbidities were also extracted and listed in Table 1.

Table 1 Baseline characteristics among the Hashimoto’s thyroiditis group and non-Hashimoto’s thyroiditis group

To minimize surveillance bias, patients who were diagnosed with SLE during the period of 2 years before the index date were excluded. Besides, due to the chronic and latent nature of SLE, patients who were diagnosed with SLE whose follow-up time was less than 3 months and 6 months were excluded, respectively, in 2 scenarios, which were conducted to increase the accuracy, and in the background variations, groups of 3 months before the index date and 3 and 6 months after the index date were also corrected with length of hospital stays and times of outpatient department visits (Table 1).

Comorbidities were captured by tracing all the ambulatory care and admission records in the NHI database within 1 previous year of the index date and have had at least three outpatient visits or one hospital admission. We analyzed autoimmune disorders including rheumatoid arthritis (RA, ICD9 code: 714.0), Sjögren’s syndrome (SS, ICD9 code: 710.2), systemic sclerosis (SSc, ICD9 code: 710.1), and vasculitis (ICD9 code: 433.0) that not rarely occur with SLE. Other common comorbidities such as hypertension (ICD9 codes: 401–405), diabetes mellitus (ICD9 code: 250), hyperlipidemia (ICD9 codes: 272.0–272.4), coronary artery disease (ICD9 codes: 410–414), osteoporosis (ICD9 code: 733), cerebral vascular accident (ICD9 codes: 430–438), chronic obstructive pulmonary disease (COPD)/asthma (ICD9 codes: 490–496), chronic kidney disease (CKD, ICD9 code: 585), chronic liver diseases (ICD9 codes: 571, 573), pancreatitis (ICD9 codes: 577.0, 577.1), affective psychosis (ICD9 code: 296), ankylosing spondylitis (ICD9 code: 720.0), inflammatory bowel disease (ICD9 codes: 555–556), HIV infection (ICD9 codes: 042–044, V08), autoimmune hemolytic anemia (AIHA) (ICD9 code: 283.0), and idiopathic thrombocytopenic purpura (ITP) (ICD9 code: 287.3) were also included in the study. Baseline treatment of HT including (1) no drug admiration, (2) anti-thyroid medication (carbimazole, propylthiouracil, methimazole)/eltroxin only, and (3) HCQ/corticosteroid+/− anti-thyroid medication/eltroxin was also analyzed, and all treatments were given within 6 months after diagnosis. Besides, hyperthyroidism (ICD9 code: 242) and hypothyroidism (ICD9 codes: 243, 244) diagnosed before the index date were separately analyzed by multivariable statistical analysis, which was listed in Tables 5, 6, and 7.

Statistical analysis

To compare and increase the similarities between our exposure group of HT and the comparator group without HT, the chi-square (χ2) tests and the two-tailed T test was used for the baseline demographic characteristics such as gender, age, urbanization, income level, admission duration, and comorbidities. Time-to-event analysis was conducted based on the index date defined as the fixed time point (January 2005) for every participation. All participants were followed up from their respective index date until the occurrence of SLE, until withdrawal, or until the end of 2013, whichever occurred first.

We also constructed a multivariable Cox proportional hazard model to estimate the hazard ratios (HRs) and 95% confidence intervals (CIs) for the SLE incidence. In the 1:20 age- and gender-matched population, 3 models were conducted. The first would be the model of HT alone, which also analyzed other thyroid disorders. For the second model, hyperthyroid and hypothyroid disorders were excluded. Model 3 contains HT with demographic variables, medical utilization, and comorbidities, and for the 1:2 PSM population, models A and B were constructed by controlling variables such as HT, hyperthyroid, and hypothyroid disorders. All the data and statistics were processed and analyzed by the Statistics Analysis System (SAS) software version 9.3 (SAS Institute, Inc., Cary, NC), and a p-value less than 0.05 was considered to indicate statistical significance.

Sensitivity analysis

To test the reliability of our study results, we established 4 sensitivity analysis scenarios including SLE medication treatment and the exclusion of autoimmune thyroiditis with other autoimmune diseases. The SLE treatment was identified as systemic corticosteroids or disease-modifying anti-rheumatic drugs (DMARDs) (including hydroxychloroquine (HCQ) or azathioprine) within 6 months after the first diagnosis of SLE. In scenarios 1–3, adjusted hazard ratio (aHR) was analyzed based on a different definition of SLE event. The main finding without medication treatment analysis would be scenario 1, systemic corticosteroids or DMARDs treatment were brought into scenario 2, and systemic corticosteroids were excluded in scenario 3. Scenario 4 modified the exclusion criteria and exclusion of the patients with rheumatic arthritis (RA), Sjögren’s syndrome (SS), systematic sclerosis (SSc), vasculitis, ankylosing spondylitis (AS), and inflammatory bowel disease (IBD) at baseline; hence, the autoimmune thyroiditis accompanied with other autoimmune diseases could be ruled out. The sensitivity analysis scenarios were listed in Tables 8, 9, and 10.


After exclusion and 1:20 age match, we identified 15,751 patients among 25,018 HT patients as our study group; 315,020 patients were extracted for the comparator group without HT. Furthermore, the 1:2 PSM filtered out 15,512 cases for the study group with 31,024 cases for the comparator group without HT. The baseline demographic characteristics, medical utilizations, and comorbidities of both groups were listed in Table 1. There was a significant higher proportion of most listed comorbidities (p < 0.001) except for CKD, IBD, and HIV infection in the HT group. As for the baseline medical treatment, over half of the cases (54.4%) in this study were under anti-thyroid medication only, such as carbimazole, propylthiouracil, methimazole, and eltroxin; about 26.3% of cases without medication control; and 19.3% patients were taking HCQ or corticosteroids. Also, Table 1 contained the baseline proportion of thyroid function disorders, which leads to our further multivariable Cox regression analysis.

Tables 2, 3, and 4 present the time-to-event analysis of the SLE incidence rate, including sensitivity analysis scenarios listed in Tables 8, 9, and 10, including data before and after PSM. Before PSM, the incidence rate ratio was similar in scenarios 1–3 (6.83, 95% CI: 5.35–8.72; 6.86, 95% CI: 5.37–8.77; 7.35 95% CI: 5.60–9.63, respectively) and slightly lower in scenario 4 (5.53, 95% CI: 4.16–7.35) which excluded other autoimmune disorders that might have a symptom of thyroiditis. However, after 1:2 PSM, the incidence rate ratio in scenario 4 became slightly higher than in scenarios 1–2 and still lower than in scenario 3. These results were also noted even though we set more barriers on following time, which are presented in Tables 3 and 4. Overall, patients with HT presented a significant increasing risk of SLE in all 4 scenarios (p, long rank p < 0.001). The cumulative probability of SLE incidence after PSM 1:2 was presented via Kaplan-Meier curves, and they were analyzed with the shortest follow-up period of 3 months and 6 months (Fig. 1a–c). In Fig. 2a–c, the factors of hyperthyroidism and hypothyroidism were discussed also. Their data were shown in Table 11.

Table 2 Incidence rate. No restriction for follow-up duration
Table 3 Incidence rate. Follow-up duration of samples 3 months
Table 4 Incidence rate. Follow-up duration of samples 6 months
Fig. 1
figure 1

a The cumulative probability of SLE in non-HT and HT patients after PSM 1:2. b The cumulative probability of SLE in non-HT and HT patients after PSM 1:2. Follow-up duration of samples 3 months. c The cumulative probability of SLE in non-HT and HT patients after PSM 1:2. Follow-up duration of samples 6 months

Fig. 2
figure 2

a The cumulative probability of SLE in non-HT and HT patients with hyperthyroidism or hypothyroidism after PSM 1:2. b The cumulative probability of SLE in non-HT and HT patients with hyperthyroidism or hypothyroidism after PSM 1:2. Follow-up duration of samples 3 months. c The cumulative probability of SLE in non-HT and HT patients with hyperthyroidism or hypothyroidism after PSM 1:2. Follow-up duration of samples 6 months

Tables 5, 6, and 7 show the results of univariable and multivariable Cox regression analyses. The adjusted hazard ratio in HT exposure alone (model 1) was 6.79 (95% CI: 5.32–8.66), which indicates that the increased risk of SLE in the HT exposure and other thyroid disorders was also analyzed, and shows that patients with either HT, hyperthyroidism, or hypothyroidism were all supposed to have an increased SLE incident risk. The aHR of model 2a, with demographic adjustment including sex, age, urbanization, low income, length of hospital stays at baseline, and times of outpatient department visits, was 5.83 (95% CI: 4.50–7.56) and showed increasing risk on long hospital stay patients. Model 3a shows that the aHR after adjustment of demographic variables, medical utilization, and comorbidities at baseline was 4.35 (95% CI: 3.28–5.76). Autoimmune diseases listed in the table including RA, SS, SSc, vasculitis, AIHA, and ITP as comorbidities also increased the incidence rate of SLE compared with the comparator group without HT. Models 2b and 3b included the variations of HT, hyperthyroidism, and hypothyroidism; these models still reveal similar results. HT, hyperthyroidism, and hypothyroidism were all supposed to increase the SLE incident risk, and the HT only group shows the highest aHR of 6.52 (4.55–9.34) in model 2b and second highest aHR of 4.53 (3.08–6.66) in model 3b; compared to HT combined with hyperthyroidism, those who have combined HT and hypothyroidism were more likely to develop SLE. Demographic variables were similar to the previous description. As for the comorbidities analysis, a high hazard ratio of SLE was also found in other autoimmune diseases such as RA, SS, SSc, and vasculitis.

Table 5 Cox proportional hazard regressions for estimation of adjusted HRs
Table 6 Cox proportional hazard regressions for estimation of adjusted HRs. Follow-up duration of samples3 months
Table 7 Cox proportional hazard regressions for the estimation of adjusted HRs. Follow-up duration of samples 6 months

After 1:2 PSM, the HR of the conditional Cox model with HT exposure alone (model A) was 3.54 (95% CI: 2.40–5.22). In model B, the group of HT and hypothyroidism only presented the highest aHR of 4.27 (95% CI: 2.67–6.83), followed by HT, no hyperthyroidism, and no hypothyroidism: 3.47 (95% CI: 2.12–5.69) and HT, hyperthyroidism, and hypothyroidism: 3.77 (95% CI: 1.76–8.05). The least risk was HT and hyperthyroidism only, aHR: 1.72 (95% CI: 0.68–4.35), with no significant difference compared with the result before PSM. The group of the 3 months barrier for following up time shares the same results, which are presented in Table 6, but the group of 6 months barrier presents the highest aHR of 4.33 (95% CI: 2.01–9.32), presented in Table 7.

We also conducted the sensitivity analysis in the estimation of the SLE risk for HT exposure in age- and sex-matched population. In the 4 scenarios, 2 different SLE treatment plan and the exclusion of other autoimmune diseases were listed in Tables 8, 9, and 10. Under the constructive of model 3, the aHR was 4.35 (95% CI: 3.28–5.76) in scenario 1, 4.39 (95% CI: 3.31–5.82) in scenario 2, 5.11 (95% CI: 3.75–6.98) in scenario 3, and 4.70 (95% CI: 3.46–6.38) in scenario 4. The sensitivity analysis was also performed on the groups of the 3 and 6 months barrier, and all 4 scenarios in the 3 groups have high enough aHRs to support the major result.

Table 8 Sensitivity analysis in the estimation of the SLE risk for Hashimoto’s thyroiditis exposure in the age-matched and sex-matched population
Table 9 Sensitivity analysis in the estimation of the SLE risk for Hashimoto’s thyroiditis exposure in the age-matched and sex-matched population. Follow-up duration of samples 3 months
Table 10 Sensitivity analysis in the estimation of the SLE risk for Hashimoto’s thyroiditis exposure in the age-matched and sex-matched population. Follow-up duration of samples 6 months


Although previous studies about thyroid and SLE showed that SLE patients are prone to develop hypothyroidism [2, 17], this study indeed told us that HT might also be associated with SLE (Table 11). In this population-based study in Taiwan, we found patients with a history of HT (aHR: 6.79, 95% CI: 5.32–8.66) or HT with hypothyroidism (aHR: 7.59, 95% CI: 5.34–10.81) were vulnerable to develop SLE compared to non-HT, no hyperthyroidism, and no hypothyroidism. Besides, hyperthyroidism was also a minor risk factor for SLE with a less ratio (aHR: 3.21, 95% CI: 1.43–7.21). More interestingly, in those HT patients who were combined with hyperthyroidism, the incidence of SLE decreased slightly but still higher than in the comparator group without HT. On the other hand, if HT patients once had hypothyroidism, then whether they had hyperthyroidism or not, the incidence of SLE is hardly affected compared to those with hypothyroidism only.

Table 11 Number at risk after PSM 1:2

The reason that patients with HT are prone to develop SLE needs to be clarified. In our opinion, first, impairment of regulatory T cells (Treg) might be a key. Impaired Treg might cause the loss of self-toleration and increases the risk of autoimmune disease, including SLE [20], and according to previous studies, the loss of Treg function was found in both HT and SLE [21,22,23]. Second, interleukin-17 (IL-17) and Th17 which are known to participate in inflammation [24] also play important roles in autoimmune diseases, including HT and SLE [25, 26]. A study pointed out that the more IL-17 is produced by Th17, the more thyroid function is lost in HT patients [27]. As SLE shares a similar pathogenesis [26], elevated Th17 and IL-17 in HT might stimulate the progression of SLE. Third, the common presence of antinuclear antibodies (ANA) in HT patients might be a crucial factor to induce other autoimmune diseases, including SLE. A study evaluating HT patients showed that 47% of HT patients were ANA positive, and 72% of them have other autoimmunity parameters besides anti-thyroid peroxidase (anti-TPO) or antithyroglobulin (anti-Tg), and/or have an autoimmune disease besides HT [28], indicating that it is possible for HT patients to come out with other autoimmune diseases. ANA is also a highly sensitive (98%) screening marker for SLE [29]. Thus, the common presence of ANA in patients with HT might be a crucial factor to induce the development of SLE. Last, studies have shown that phagocytosis was stimulated by physiological concentrations of thyroid hormone [30], and it was decreased after thyroid suppression in an animal model [31]. It means HT patients with hypothyroidism might lose the ability to clean up the autoimmune complex and develop SLE [32], which might explain why patients with hypothyroidism have a higher risk of SLE than hyperthyroidism in our study.

The finding of HT being associated with subsequent SLE is important in clinical practice. HT is not a difficult disease to treat, however, its complications are usually forgotten, for example, thyroid lymphoma [33]. And now SLE might also be another potential complication of HT according to our study. Once SLE has been developed, it will affect the whole body from sleep disturbance [34] to pulmonary complications [35] and cardiovascular disease [36], deteriorating the health condition. Moreover, since genetic factors cannot be modified easily, it seems important to investigate the other risk factors of SLE. Early identification and intervention of SLE will minimize the deterioration and damage to the body [37].

Previous studies have told us some environmental factors causing SLE. Besides ultraviolet radiation [38], bacteria and virus infections have also been proven that they are linked to SLE by affecting the immune system [39, 40], including nontyphoidal– Salmonella [41] and Varicella zoster virus [42]. Hormones have also been investigated for their correlation with SLE. For example, estrogen can rescue autoreactive B cells from apoptosis [43], explaining why females are more prone to autoimmune disease. Progesterone can regulate CD4+ T cells [44] to prevent pregnant women from producing anti-fetal antibodies, so low levels of progesterone are considered as a predisposing factor for SLE [45]. According to many studies, thyroid hormone also can affect the immune system [10,11,12,13,14,15,16]. Acquiring any kind of autoimmune disease is another risk to developing another one with unclear reasons, and when there is more than three autoimmune coexistence, it is called as “multiple autoimmune syndrome (MAS)” [8, 9], and HT is classified as one of MAS [46], which might be a risk factor of developing SLE.

Despite that many studies indicated the function of thyroid hormone in the immune system [10,11,12,13,14,15,16] and the hypothesis about MAS increasing the risk of SLE, there are only two case reports about two young girls and two women, respectively, implying HT might be a risk factor of SLE [18, 19]. This time, we conducted a large-scale retrospective cohort study, which is more persuasive than a case report, attempting to prove this thought.

In our study, we not only surveyed the relation between HT and SLE, but also consider other factors. We have described above the thought of MAS being a factor in developing SLE, and now, our data supported it again with an extremely higher hazard ratio of SLE with other autoimmune diseases, including RA (aHR: 14.57, CI: 9.51–22.33), SS (aHR: 20.96, CI: 13.92–31.56), SSc (aHR: 61.44, CI: 19.78–190.79), vasculitis (aHR: 81.34, CI: 36.38–181.88), AIHA (aHR: 417.97, CI: 104.30–1674.87), and ITP (aHR: 67.09, CI: 27.81–161.85). On the other hand, we also found the possibility of HT alone increasing the risk of SLE by excluding other common autoimmune diseases (aHR: 4.70, CI: 3.46–6.38), which means SLE following HT might not be attributed to other autoimmune diseases.

In addition, we considered whether the status of thyroid function would affect the risk of SLE. The results showed a higher risk of SLE among patients with hypothyroidism. Despite the absolute rate being low, the increased hazard ratio of SLE gave clinical physicians a hint for caring for these HT patients. This outcome might be explained by the ability of phagocytosis of the autoimmune complex which we have mentioned above [30,31,32]. However, the detailed mechanism between HT and SLE still needs to be clarified.

Our study was validated enough to be a presentative of the general population by using the NHIRD of Taiwan, which has multiple advantages including a great sample size covering over 99% of nationals of Taiwan, and long-term comprehensive follow-up to assess the risk of new-onset SLE in patients with HT [47]. In addition, we performed a sensitivity analysis to confirm the conclusion of this study by using four different scenarios. In other subgroups, urbanization, low income, length of hospital stays, and other factors were examined. Gender, age, and other factors were adjusted appropriately in PSM to minimize the selection bias.

Despite the advantages mentioned above, there are still some limitations in our study. First, the use of a “major illness registry” might lead to an underestimated incidence of SLE because few patients might not get this identification. Second, smoking status, a confounder factor of SLE [48], was unavailable in the NHIRD. Hence, we used chronic obstructive pulmonary disease as a surrogate variable for cigarette smoking because of its close correlation with cigarettes [49]. Third, our study cannot explain why HT patients with hyperthyroidism had a decreased incidence of SLE despite the fact that excessive thyroid hormone can also impair Treg cells [50]. It cannot prove the pathogenesis of SLE following HT directly, either. Thus, the mechanism still needs to be investigated through more experiments. Fourth, although we have adjusted many confounders, we still missed genetic factors which are associated with both HT and SLE due to the difficulty of collecting genetic data. Last, because people recorded in NHIRD are usually Taiwanese, this study might not be applicable to non-Asian ethnic groups.

In conclusion, this population-based study suggested an increased risk of SLE in the HT group after adjustment for baseline characteristics, comorbidities, and medical confounders compared with the reference group. It could provide hints for further research to clarify the pathogenesis between HT, hypothyroidism, hyperthyroidism, and SLE.

Availability of data and materials

Summarized individual data are available on request to the corresponding author. The data set used in this study is managed by the Taiwan Ministry of Health and Welfare and, thus, cannot be made available publicly. Researchers interested in accessing this data set can submit a formal application to the Ministry of Health and Welfare to request access (the postal address No. 488, Section 6, Zhongxiao E Rd, Nan-gang District, Taipei City 115, Taiwan; website: DOS/cp-2516-3591-113.html).



Systemic lupus erythematosus


Hashimoto’s thyroiditis


Disease-modifying anti-rheumatic drugs




Rheumatic arthritis


Sjögren’s syndrome


Systematic sclerosis


Ankylosing spondylitis


Inflammatory bowel disease


Chronic obstructive pulmonary disease


Chronic kidney disease


Autoimmune hemolytic anemia


Idiopathic thrombocytopenic purpura


Antinuclear antibody


Multiple autoimmune syndrome


  1. Fava A, Petri M. Systemic lupus erythematosus: diagnosis and clinical management. J Autoimmun. 2019;96:1–13.

    Article  PubMed  Google Scholar 

  2. Ferrari SM, Elia G, Virili C, Centanni M, Antonelli A, Fallahi P. Systemic lupus erythematosus and thyroid autoimmunity. Front Endocrinol (Lausanne). 2017;8:138.

    Article  PubMed  Google Scholar 

  3. Gravano DM, Hoyer KK. Promotion and prevention of autoimmune disease by CD8+ T cells. J Autoimmun. 2013;45:68–79.

    Article  CAS  PubMed  Google Scholar 

  4. Yurasov S, Wardemann H, Hammersen J, Tsuiji M, Meffre E, Pascual V, et al. Defective B cell tolerance checkpoints in systemic lupus erythematosus. J Exp Med. 2005;201(5):703–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Deapen D, Escalante A, Weinrib L, Horwitz D, Bachman B, Roy-Burman P, et al. A revised estimate of twin concordance in systemic lupus erythematosus. Arthritis Rheum. 1992;35(3):311–8.

    Article  CAS  PubMed  Google Scholar 

  6. Cooper GS, Miller FW, Pandey JP. The role of genetic factors in autoimmune disease: implications for environmental research. Environ Health Perspect. 1999;107(Suppl 5):693–700.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Cooper GS, Dooley MA, Treadwell EL, St Clair EW, Parks CG, Gilkeson GS. Hormonal, environmental, and infectious risk factors for developing systemic lupus erythematosus. Arthritis Rheum. 1998;41(10):1714–24.

    Article  CAS  PubMed  Google Scholar 

  8. Cojocaru M, Cojocaru IM, Silosi I. Multiple autoimmune syndrome. Maedica (Bucur). 2010;5(2):132–4.

    CAS  PubMed  Google Scholar 

  9. Sloka S. Observations on recent studies showing increased co-occurrence of autoimmune diseases. J Autoimmun. 2002;18(3):251–7.

    Article  PubMed  Google Scholar 

  10. De Vito P, Incerpi S, Pedersen JZ, Luly P, Davis FB, Davis PJ. Thyroid hormones as modulators of immune activities at the cellular level. Thyroid. 2011;21(8):879–90.

    Article  PubMed  Google Scholar 

  11. Rubingh J, van der Spek A, Fliers E, Boelen A. The role of thyroid hormone in the innate and adaptive immune response during infection. Compr Physiol. 2020;10(4):1277–87.

    Article  PubMed  Google Scholar 

  12. Robinson MV, Obut TA, Melnikova EV, Trufakin VA. Parameters of cellular and humoral immunity in experimental hyperthyroidism and its correction. Bull Exp Biol Med. 2014;156(4):473–5.

    Article  CAS  PubMed  Google Scholar 

  13. Robinson MV, Obut TA, Mel’nikova EV, Trufakin VA. Effects of thyroxin and mercazolyl on immunological parameters of blood lymphocytes and lymphoid organs. Bull Exp Biol Med. 2013;156(2):236–8.

    Article  CAS  PubMed  Google Scholar 

  14. Marino F, Guasti L, Cosentino M, De Piazza D, Simoni C, Piantanida E, et al. Thyroid hormone regulation of cell migration and oxidative metabolism in polymorphonuclear leukocytes: clinical evidence in thyroidectomized subjects on thyroxine replacement therapy. Life Sci. 2006;78(10):1071–7.

    Article  CAS  PubMed  Google Scholar 

  15. Ortega E, Forner MA, Garcia JJ, Rodriguez AB, Barriga C. Enhanced chemotaxis of macrophages by strenuous exercise in trained mice: thyroid hormones as possible mediators. Mol Cell Biochem. 1999;201(1-2):41–7.

    Article  CAS  PubMed  Google Scholar 

  16. Magsino CH Jr, Hamouda W, Ghanim H, Browne R, Aljada A, Dandona P. Effect of triiodothyronine on reactive oxygen species generation by leukocytes, indices of oxidative damage, and antioxidant reserve. Metabolism. 2000;49(6):799–803.

    Article  CAS  PubMed  Google Scholar 

  17. Watad A, Mahroum N, Whitby A, Gertel S, Comaneshter D, Cohen AD, et al. Hypothyroidism among SLE patients: case-control study. Autoimmun Rev. 2016;15(5):484–6.

    Article  CAS  PubMed  Google Scholar 

  18. Bakr A, Laimon W, El-Ziny MA, Hammad A, El-Hawary AK, Elsharkawy AA, et al. The emergence of systemic lupus erythematosus in hypothyroid patients: two case reports and mini review. Lupus. 2014;23(8):825–8.

    Article  CAS  PubMed  Google Scholar 

  19. Dhir R, Ahluwalia AI, Sridhar J, Mani H, Pruthi HS, Shah KM. Autoimmune thyroiditis perdating the presentation of systemic lupus erythematosus: two cases and a review of literature. Indian J Dermatol Venereol Leprol. 2002;68(5):292–4.

    PubMed  Google Scholar 

  20. Li W, Deng C, Yang H, Wang G. The regulatory T cell in active systemic lupus erythematosus patients: a systemic review and meta-analysis. Front Immunol. 2019;10:159.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Glick AB, Wodzinski A, Fu P, Levine AD, Wald DN. Impairment of regulatory T-cell function in autoimmune thyroid disease. Thyroid. 2013;23(7):871–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Scheinecker C, Bonelli M, Smolen JS. Pathogenetic aspects of systemic lupus erythematosus with an emphasis on regulatory T cells. J Autoimmun. 2010;35(3):269–75.

    Article  CAS  PubMed  Google Scholar 

  23. Hu Y, Zhang L, Chen H, Liu X, Zheng X, Shi H, et al. Analysis of regulatory T cell subsets and their expression of Helios and PD-1 in patients with hashimoto thyroiditis. Int J Endocrinol. 2019;2019:5368473.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Miossec P, Kolls JK. Targeting IL-17 and TH17 cells in chronic inflammation. Nat Rev Drug Discov. 2012;11(10):763–76.

    Article  CAS  PubMed  Google Scholar 

  25. Tabarkiewicz J, Pogoda K, Karczmarczyk A, Pozarowski P, Giannopoulos K. The role of IL-17 and Th17 lymphocytes in autoimmune diseases. Arch Immunol Ther Exp (Warsz). 2015;63(6):435–49.

    Article  CAS  PubMed  Google Scholar 

  26. Pernis AB. Th17 cells in rheumatoid arthritis and systemic lupus erythematosus. J Intern Med. 2009;265(6):644–52.

    Article  CAS  PubMed  Google Scholar 

  27. Li D, Cai W, Gu R, Zhang Y, Zhang H, Tang K, et al. Th17 cell plays a role in the pathogenesis of Hashimoto’s thyroiditis in patients. Clin Immunol. 2013;149(3):411–20.

    Article  CAS  PubMed  Google Scholar 

  28. Fiducia M, Lauretta R, Lunghi R, Kyanvash S, Pallotti S. Hashimoto’s thyroiditis and autoimmunity parameters: descriptive study. Minerva Med. 2007;98(2):95–9.

    CAS  PubMed  Google Scholar 

  29. Wichainun R, Kasitanon N, Wangkaew S, Hongsongkiat S, Sukitawut W, Louthrenoo W. Sensitivity and specificity of ANA and anti-dsDNA in the diagnosis of systemic lupus erythematosus: a comparison using control sera obtained from healthy individuals and patients with multiple medical problems. Asian Pac J Allergy Immunol. 2013;31(4):292–8.

    Article  PubMed  Google Scholar 

  30. Perrotta C, Buldorini M, Assi E, Cazzato D, De Palma C, Clementi E, et al. The thyroid hormone triiodothyronine controls macrophage maturation and functions: protective role during inflammation. Am J Pathol. 2014;184(1):230–47.

    Article  CAS  PubMed  Google Scholar 

  31. Hampton LMT, Finch MG, Martyniuk CJ, Venables BJ, Jeffries MKS. Developmental thyroid disruption causes long-term impacts on immune cell function and transcriptional responses to pathogen in a small fish model. Sci Rep. 2021;11(1):14496.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Ma C, Xia Y, Yang Q, Zhao Y. The contribution of macrophages to systemic lupus erythematosus. Clin Immunol. 2019;207:1–9.

    Article  CAS  PubMed  Google Scholar 

  33. Dündar HZ, Sarkut P, Kırdak T, Korun N. Primary thyroid lymphoma. Ulus Cerrahi Derg. 2016;32(1):75–7.

    PubMed  Google Scholar 

  34. Zhao Q, Deng N, Chen S, Cui Y, Du X, Gu Z. Systemic lupus erythematosus is associated with negatively variable impacts on domains of sleep disturbances: a systematic review and meta-analysis. Psychol Health Med. 2018;23(6):685–97.

    Article  PubMed  Google Scholar 

  35. Hannah JR, D’Cruz DP. Pulmonary complications of systemic lupus erythematosus. Semin Respir Crit Care Med. 2019;40(2):227–34.

    Article  PubMed  Google Scholar 

  36. Torres A, Askari AD, Malemud CJ. Cardiovascular disease complications in systemic lupus erythematosus. Biomark Med. 2009;3(3):239–52.

    Article  CAS  PubMed  Google Scholar 

  37. Samnaliev M, Barut V, Weir S, Langham J, Langham S, Wang X, et al. Health-care utilization and costs in adults with systemic lupus erythematosus in the United Kingdom: a real-world observational retrospective cohort analysis. Rheumatol Adv Pract. 2021;5(3):rkab071.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Bijl M, Kallenberg CG. Ultraviolet light and cutaneous lupus. Lupus. 2006;15(11):724–7.

    Article  CAS  PubMed  Google Scholar 

  39. Zandman-Goddard G, Shoenfeld Y. Infections and SLE. Autoimmunity. 2005;38(7):473–85.

    Article  CAS  PubMed  Google Scholar 

  40. Navarra SV, Leynes MS. Infections in systemic lupus erythematosus. Lupus. 2010;19(12):1419–24.

    Article  CAS  PubMed  Google Scholar 

  41. Tu TY, Yeh CY, Hung YM, Chang R, Chen HH, Wei JC. Association between a history of nontyphoidal Salmonella and the risk of systemic lupus erythematosus: a population-based, case-control study. Front Immunol. 2021;12:725996.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Sun F, Chen Y, Wu W, Guo L, Xu W, Chen J, et al. Varicella zoster virus infections increase the risk of disease flares in patients with SLE: a matched cohort study. Lupus Sci Med. 2019;6(1):e000339.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Cohen-Solal JF, Jeganathan V, Hill L, Kawabata D, Rodriguez-Pinto D, Grimaldi C, et al. Hormonal regulation of B-cell function and systemic lupus erythematosus. Lupus. 2008;17(6):528–32.

    Article  CAS  PubMed  Google Scholar 

  44. Hughes GC, Clark EA, Wong AH. The intracellular progesterone receptor regulates CD4+ T cells and T cell-dependent antibody responses. J Leukoc Biol. 2013;93(3):369–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Arnalich F, Benito-Urbina S, Gonzalez-Gancedo P, Iglesias E, de Miguel E, Gijon-Baños J. Inadequate production of progesterone in women with systemic lupus erythematosus. Br J Rheumatol. 1992;31(4):247–51.

    Article  CAS  PubMed  Google Scholar 

  46. Shahbaz A, Aziz K, Umair M, Sachmechi I. Prolonged duration of hashitoxicosis in a patient with Hashimoto’s thyroiditis: a case report and review of literature. Cureus. 2018;10(6):e2804.

    PubMed  PubMed Central  Google Scholar 

  47. Hsing AW, Ioannidis JP. Nationwide population science: lessons from the Taiwan National Health Insurance Research Database. JAMA Intern Med. 2015;175(9):1527–9.

    Article  PubMed  Google Scholar 

  48. Barbhaiya M, Tedeschi SK, Lu B, Malspeis S, Kreps D, Sparks JA, et al. Cigarette smoking and the risk of systemic lupus erythematosus, overall and by anti-double stranded DNA antibody subtype, in the Nurses’ Health Study cohorts. Ann Rheum Dis. 2018;77(2):196–202.

    Article  CAS  PubMed  Google Scholar 

  49. Chung WS, Kung PT, Chang HY, Tsai WC. Demographics and medical disorders associated with smoking: a population-based study. BMC Public Health. 2020;20(1):702.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Zhong Y, Lu TT, Liu XM, Liu BL, Hu Y, Liu S, et al. High levels of thyroid hormone impair regulatory T cell function via reduced PD-1 expression. J Clin Endocrinol Metab. 2021;106(9):2738–53.

    Article  PubMed  Google Scholar 

Download references


This study was supported by Taichung Veterans General Hospital.

Patient and Public Involvement patients and/or the public were not involved in the design, conduct, reporting, or dissemination of plans of this research.

Patient consent for publication is not required because the data we used comprised a de-identified secondary data set which had been released for research purposes and analyzed anonymously.


This study did not receive any funding support.

Author information

Authors and Affiliations



All authors were involved in the drafting of the article or revising it, and all authors approved the final version to be published. Study conception and design: H-HC and JC-CW. Acquisition of data: H-HC and JC-CW. Analysis and interpretation of the data: H-HC, H-CL, H-MC, and Y-MH. Writing (original draft preparation): H-CL and H-MC. Writing (review and editing): Y-MH and RC. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Hsin-Hua Chen or James Cheng-Chung Wei.

Ethics declarations

Ethics approval and consent to participate

The Institutional Review Board of Taichung Veterans General Hospital (number: CE17100B) approved this study. The requirement for informed consent was waived as personal information was anonymized before data analyses.

Consent for publication

In submitting an article to any of the journals published by BMC, I certify that:

I am authorized by my co-authors to enter into these arrangements.

I warrant, on behalf of myself and my co-authors, that:

The article is original, has not been formally published in any other peer-reviewed journal, is not under consideration by any other journal, and does not infringe any existing copyright or any other third-party rights.

I am/we are the sole authors of the article and have full authority to enter into this agreement and in granting rights to BMC are not in breach of any other obligation.

The article contains nothing that is unlawful, libelous, or which would, if published, constitute a breach of contract or of confidence or of commitment given to secrecy.

I/we have taken due care to ensure the integrity of the article. To my/our—and currently accepted scientific—knowledge, all statements contained in it purporting to be facts are true, and any formula or instruction contained in the article will not, if followed accurately, cause any injury, illness, or damage to the user.

I, and all co-authors, agree that the article, if editorially accepted for publication, shall be licensed under the Creative Commons Attribution License 4.0. In line with BMC’s Open Data Policy, data included in the article shall be made available under the Creative Commons 1.0 Public Domain Dedication waiver, unless otherwise stated. If the law requires that the article be published in the public domain, I/we will notify BMC at the time of submission, and in such cases, not only the data but also the article shall be released under the Creative Commons 1.0 Public Domain Dedication waiver. For the avoidance of doubt, it is stated that sections 1 and 2 of this license agreement shall apply and prevail regardless of whether the article is published under Creative Commons Attribution License 4.0 or the Creative Commons 1.0 Public Domain Dedication waiver.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lin, HC., Chang, HM., Hung, YM. et al. Hashimoto’s thyroiditis increases the risk of new-onset systemic lupus erythematosus: a nationwide population-based cohort study. Arthritis Res Ther 25, 20 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Hashimoto’s thyroiditis
  • Systemic lupus erythematosus
  • Retrospective cohort study