Skip to main content

SYSTEMATIC REVIEW article

Front. Endocrinol., 13 October 2022
Sec. Thyroid Endocrinology

Association between rheumatoid arthritis and thyroid dysfunction: A meta-analysis and systematic review

  • 1Department of Rheumatology and Immunology, The Second Affiliated Hospital of Fujian Medical University, Quanzhou, China
  • 2Centre of Neurological and Metabolic Research, The Second Affiliated Hospital of Fujian Medical University, Quanzhou, China
  • 3Group of Neuroendocrinology, Garvan Institute of Medical Research, Sydney, NSW, Australia

Objective: Rheumatoid arthritis (RA) is an autoimmune disorder. Multiple studies have investigated the risk of thyroid dysfunction in patients with RA but have reached conflicting conclusions. This systematic review aimed to determine whether patients with RA are at higher risk of thyroid dysfunction.

Methods: We comprehensively reviewed online literature databases, including PubMed, Scopus, Embase, and the Cochrane Library, from their respective inception dates to March 25, 2022. Studies that provided data on at least one case of thyroid dysfunction in RA patients and their controls were included. Based on these data, we calculated pooled odds ratios (ORs) and their corresponding 95% confidence intervals (CIs) for thyroid dysfunction in RA and non-RA patients.

Results: Twenty-nine studies met the inclusion criteria, involving a total of 35,708 patients with RA. The meta-analysis showed that, compared with non-RA patients, RA patients had an increased risk of developing thyroid dysfunction, particularly hypothyroidism (OR 2.25, 95% CI 1.78–2.84). Subgroup analysis suggested that study type and sample source of control group were the source of heterogeneity.

Conclusions: Patients with RA are at increased risk of developing thyroid dysfunction, especially hypothyroidism. Routine biochemical examination of thyroid function in RA patients should be strengthened. Larger prospective studies are needed to explore the causal relationship between RA and thyroid dysfunction, and to investigate the impact of thyroid dysfunction on RA disease activity, drug efficacy, and medication safety.

Systematic review registration: https://www.crd.york.ac.uk/prospero/, identifier CRD42022331142.

Introduction

Rheumatoid arthritis (RA) is an autoimmune disease characterized by systemic inflammation, persistent synovitis, and autoantibodies (particularly against citrullinate peptides and rheumatoid factors) (1). With repeated episodes of joint inflammation, the normal structure of the joint is destroyed, ultimately resulting in reduced mobility and increased disability (2). In addition, many RA patients also present with extra-articular clinical symptoms (3). Recently, a meta-analysis (4) determined that the global prevalence of RA was approximately 0.46% over the past 40 years. Moreover, RA has been found to be potentially associated with a variety of autoimmune diseases, including thyroid dysfunction (5). Thyroid dysfunction mainly includes hyperthyroidism, hypothyroidism, subclinical hyperthyroidism and subclinical hypothyroidism. Autoimmune thyroid disease (AITD), which mainly includes Hashimoto’s thyroiditis and Graves’ disease, is the most common cause of various thyroid dysfunction and is manifested by the production of antibodies against thyroid peroxidase, thyroglobulin or thyrotropin receptor autoantigens (6, 7). In addition to clinical manifestations, laboratory tests are more important in the definition of several thyroid diseases. Common test indicators include thyroid-stimulating hormone (TSH), triiodothyronine (T3), thyroxine (T4), free triiodothyronine (fT3), free thyroxine (fT4) (8). The range of reference values given by different experimental equipment and detection kits varies to a certain extent, but the trend of indicator change has been widely recognized (9, 10). Hypothyroidism and hyperthyroidism can be divided into overt and subclinical stages. Based on hormone levels, hyperthyroidism and hypothyroidism are defined as excessive and insufficient thyroid hormones, respectively (11). Most thyroid diseases and rheumatoid arthritis are chronic diseases, and unfortunately, early signs of thyroid dysfunction are not specific and are often overshadowed by clinical manifestations of other diseases in the absence of biochemical tests of thyroid function (12). Thyroid dysfunction can be detrimental to health and, if left undiagnosed or treated, can have serious adverse consequences and, in exceptional cases, can be fatal (12).

Findings from studies on the risk of thyroid dysfunction in RA patients are conflicting (13). For example, a study by McCoy et al. (14) concluded that RA patients did not have significantly increased risk of hypothyroidism as compared with non-RA patients (14). Graves’ disease is the most common cause of hyperthyroidism and has a high comorbidity rate with RA (15). Considering the accumulation of evidence, we conducted a systematic review and meta-analysis to assess whether patients with RA are at a higher risk of the four major thyroid dysfunctions. We also used subgroup analysis to investigate the influence of control group origin, study type, regions, and publication date on the results of the overall analysis, as well as to explore sources of heterogeneity.

Methods

Search strategy

A systematic review and a meta-analysis were conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA) guidelines (16). We searched PubMed, Scopus, Embase, and the Cochrane Library databases, for all eligible studies, regardless of language, from their respective inception dates to 25 March 2022. The search terms used were “Thyroid Dysfunction,” “Thyroid Diseases,” “Hypothyroidism,” “Hyperthyroidism,” and “Arthritis, Rheumatoid” (please see appendix eTable 1 for the complete search strategy). This study is registered on PROSPERO, number CRD42022331142.

Inclusion and exclusion criteria

A Population, Intervention, Comparator, Outcome and Study design (PICOS) scheme was used to clarify our research objectives. Our meta-analysis was designed to answer the question: does rheumatoid arthritis increase the risk of thyroid dysfunction? Representatives of the PICOS regimen as follows: patients with RA (P); comparison with non-RA, non-systemic lupus erythematosus (SLE), or healthy subjects (C); the diagnosis of thyroid dysfunction (O); observational studies (S). In our meta-analysis, patients with thyroid dysfunction mainly included hyperthyroidism, hypothyroidism, subclinical hypothyroidism, and subclinical hyperthyroidism.

Studies that up to the following standards were incorporated in this meta-analysis: Observational studies that included RA and non-RA controls provided the number of samples in both groups and provided the number or prevalence of thyroid dysfunction in each group. Patients who met the American College of Rheumatology diagnostic criteria for RA or the European League Against Rheumatism classification criteria, or who had at least one documented rheumatoid arthritis diagnosis in the medical records of hospitals and medical databases, were included in the study. In addition, to be included in this study, at least one case of hypothyroidism, subclinical hypothyroidism, hyperthyroidism, or subclinical hyperthyroidism had to be detected in RA patients and/or controls. To test whether the findings are common across countries, we did not restrict the language used to publish the articles.

We excluded reviews, editorials, case reports, comments, books, etc. Studies without available abstracts or complete manuscripts were also excluded. If there was duplication in the study population among studies, only the latest article was included. Given that a strong association between thyroid disease and systemic lupus erythematosus has been reported for more than 50 years (17), studies with SLE patients as controls were also excluded.

Data extraction and quality assessment

Two reviewers (Y.L. and H.M.) independently fetched the data from eligible papers and then combined the data. Disagreements were settled through discussion to consensus. We extracted the following basic information and outcome-related data from studies that up to our inclusion criteria: year of publication, the name of the first author, nationality of the author, average duration of RA, study design, the number of cases and controls, the mean age and sex ratio among RA and non-RA patients, methodological quality of the study, and number of people with thyroid dysfunction. We focused on the statistical results of the four types of thyroid dysfunction in the RA group and the control groups. We preferentially extracted the raw data in the study report, and if the study did not provide the raw data, the prevalence value provided by the study was used for calculation. We used the Newcastle-Ottawa Scale (NOS) to assess the quality of the included literature (18). The overall score is 9, with 7-9 being high quality, 4-6 being medium quality, and 0-3 being low quality.

Statistical analysis

The association strength between RA and four types of thyroid dysfunction was calculated, expressed by the odds ratio (OR), with the corresponding 95% confidence interval (CI), using the sample sizes of RA and control groups as well as the statistical results of four types of thyroid dysfunction. OR > 1 indicates that RA increases the risk of thyroid dysfunction, OR = 1 indicates that there is no correlation between RA and thyroid dysfunction, and OR < 1 indicates that RA patients have a reduced risk of thyroid dysfunction.

P < 0.05 was recognized statistically significant. We used the I2 statistic and χ2 test to assess the heterogeneity of the pooled ORs (P < 0.1 considered statistically significant). I2 ≤ 25% indicated that heterogeneity was insignificant, I2 < 50% indicated low heterogeneity, I2 between 50% and 75% indicated medium heterogeneity, and I2 > 75% indicated high heterogeneity (19).

When the pooled OR was calculated, if I2 > 50%, the random effect model was used for sensitivity analysis and subgroup analysis; if I2 < 50%, the fixed-effect model was adopted (20). Egger’s test, Funnel plot, and Begg’s test were used to evaluate whether publication bias existed. P < 0.05 indicated that the detected publication bias was statistically significant, and the impact of potential publication bias was further assessed using the Duval and Tweedie non-parametric modified trim-and-fill method (21). All statistical analyses were accomplished using Stata V.16.0 software.

Results

Literature search results and study characteristics

A preliminary search of four electronic databases, PubMed, Scopus, Embase, and the Cochrane Library, yielded 9,746 articles, including 1,549 PubMed articles, 5,378 Embase articles, 2,673 Scopus articles, and 146 Cochrane Library articles. After removing duplicates, 6,287 records remained. Two authors (Y.L. and H.M.) scanned the titles and abstracts of these articles as a preliminary screen and selected 156 articles that appeared eligible for inclusion. Further, the two authors then read the full text of these articles respectively. If there were disputes about inclusion of any of these papers, a third author (Z.C.) made the final decision. Finally, 29 studies (5, 11, 13, 14, 2246) that met all inclusion criteria were absorbed in the meta-analysis. Figure 1 shows the detailed screening steps.

FIGURE 1
www.frontiersin.org

Figure 1 Flow chart of search results.

Table 1 describes the general characteristics of the literature included in our meta-analysis. Twenty-nine studies met the inclusion criteria, involving 35,708 RA patients and 149,421 non-RA controls. Among them, 15 studies (5, 11, 13, 22, 24, 25, 27, 31, 33, 34) provided (3740, 45) prevalence values for hyperthyroidism, 12 provided (5, 11, 22, 24, 31, 33, 34, 36, 37, 42, 44, 45) prevalence values for subclinical hyperthyroidism, 26 provided (5, 11, 13, 14, 2231, 3341, 43, 45, 46) prevalence values for hypothyroidism, and 17 provided (5, 11, 14, 22, 24, 3134, 3639, 4245) prevalence values for subclinical hypothyroidism. The NOS checklist was used to evaluate the quality of studies, and 7 studies (24.13%) assessed as high quality, and 18 studies (62.07%) identified as medium quality.

TABLE 1
www.frontiersin.org

Table 1 Characteristic of basic information.

Synthesis of results

RA and hyperthyroidism

Fifteen studies (5, 11, 13, 22, 24, 25, 27, 31) recorded the prevalence of hyperthyroidism in RA and control groups (33, 34, 3740, 45). However, in 5 literatures (27, 31, 33, 34, 45), there is no information on hyperthyroidism in RA patients or controls.

We executed a statistical analysis on the remaining 10 studies (Figure 2A). The risk of developing hyperthyroidism in RA patients was observably higher than in control groups (OR 1.65, 95% CI 1.24–2.19). Among the studies reporting statistically significant differences, some extremes occurred, such as in the study by AI-Awadhi et al. (39) suggested that RA patients had a higher risk of hyperthyroidism, while the study of Mahagna et al. (13) found that RA patients had a relatively low risk of hyperthyroidism.

FIGURE 2
www.frontiersin.org

Figure 2 Forest plot for the risk of thyroid dysfunction among RA patients. (A) Forest plot for the risk of hyperthyroidism among RA patients, (B) Forest plot for the risk of subclinical hyperthyroidism among RA patients, (C) Forest plot for the risk of hypothyroidism among RA patients, (D) Forest plot for the risk of subclinical hypothyroidism among RA patients.

RA and subclinical hyperthyroidism

After excluding two studies (31, 44), in which subclinical hyperthyroidism was not found in any subjects, the prevalence of subclinical hyperthyroidism in RA patients and controls was calculated in a total of 10 studies (5, 11, 22, 24, 33, 34, 36, 37, 42, 45) (Figure 2B). The pooled OR was 2.13 (95%CI 1.25–3.63) (P < 0.05), indicating that RA patients had more than twice the risk of subclinical hyperthyroidism as non-RA subjects, which was statistically significant. Due to low heterogeneity (I2 = 0.0%, P = 0.743), a fixed-effect model was used for statistical analysis.

RA and hypothyroidism

In the study conducted by Acay et al. (33), neither group included any cases of hypothyroidism. We thus conducted a combined analysis of the remaining 25 studies (5, 11, 13, 14, 2228, 31, 3441, 43, 45, 46) (Figure 2C). Due to the high I2 value (I2 = 75.7%), a random-effects model was applied. The pooled OR was 2.25 (95%CI 1.78–2.84) (P > 0.05). Our statistical results indicated that RA patients had a higher risk of hypothyroidism than of the other three types of thyroid dysfunction.

RA and subclinical hypothyroidism

A pooled analysis of 17 studies (5, 11, 14, 22, 24, 3134, 3639, 4245) involving the association between RA and subclinical hypothyroidism was performed. The pooled OR was 2.18 (95% CI 1.32–3.61) (P < 0.05), indicating that the risk of subclinical hypothyroidism in RA patients was significantly higher than that in controls (Figure 2D), with an I2 value of 63.8% (P < 0.001); thus, we used a random effects model.

Subgroup analyses

In addition to the pooled analysis of the risk of developing subclinical hyperthyroidism in patients with RA, there was significant heterogeneity in the pooled analysis of the other three types of thyroid dysfunctions. To further explore possible sources of heterogeneity and to determine whether specific factors influence the risk intensity for thyroid dysfunction among RA patients, we performed subgroup analyses by control group origin, study type, geographic region, and publication date. As shown in Table 2, the subgroup analysis showed that the origin of controls was the source of heterogeneity in the risk of hyperthyroidism, and the study type was the source of heterogeneity in the risk of subclinical hypothyroidism. Region and publication data did not significantly contribute to heterogeneity.

TABLE 2
www.frontiersin.org

Table 2 Subgroup analysis of region, year of Publication, study design, and control group origin.

Sensitivity analysis

Sensitivity analysis was carried out by removing each of the studies consecutively to verify the reliability and stability of the pooled results, observe the effect of individual results on the overall analysis, and look for possible sources of heterogeneity.

Sensitivity analyses were performed for association of RA and four thyroid dysfunctions. If the pooled results suggested an increased risk of thyroid dysfunction even after exclusion of any study, the pooled analysis was considered robust. As shown in Figure 3, when sensitivity analysis was performed on the pooled results of studies related to subclinical hyperthyroidism, hypothyroidism, and subclinical hypothyroidism, the combined OR remained within the overall 95% CI range after the exclusion of any study, indicating that the statistical data were robust and reliable. However, in the sensitivity analysis of hyperthyroidism-related studies, the combined OR did not remain within the 95% CI range of the overall OR after excluding the studies of Mahagna et al. (13) or Huang et al. (25), and these two studies were considered to have an impact on the overall efficacy. Heterogeneity was significantly reduced after we excluded these two articles (P = 0.186, I2 = 30.3%). Therefore, the two studies were considered to represent a source of heterogeneity. Without the two studies, the combined OR of the remaining studies was 3.03 (95%CI 1.33–6.91, P = 0.008).

FIGURE 3
www.frontiersin.org

Figure 3 Sensitivity analysis based on excluding each study in turn. (A) Sensitivity analyses for the association of RA and hyperthyroidism. (B) Sensitivity analyses for the association of RA and subclinical hyperthyroidism. (C) Sensitivity analyses for the association of RA and hypothyroidism. (D) Sensitivity analyses for the association of RA and subclinical hypothyroidism.

Publication bias

To detect whether publication bias existed, we created funnel plots, and conducted Begg’s and Egger’s tests. Typically, precise studies or studies with large numbers of subjects appear at the top of a funnel plot, while less precise studies or studies with fewer subjects appear at the bottom (21). As shown in Figure 4, these plots showed that there was no significant bias in these studies except for the abnormal funnel plot associated with RA and hypothyroidism (Figure 4C). The detection of publication bias was also confirmed in Egger’s test (eTable 2).

FIGURE 4
www.frontiersin.org

Figure 4 Publications bias detected by funnel plot on the risk of thyroid dysfunction in RA patients. (A) Publications bias detected by funnel plot on the risk of hyperthyroidism in RA patients. (B) Publications bias detected by funnel plot on the risk of subclinical hyperthyroidism in RA patients. (C) Publications bias detected by funnel plot on the risk of hypothyroidism in RA patients. (D) Publications bias detected by funnel plot on the risk of subclinical hypothyroidism in RA patients.

We attempted to evaluate the impact of publication bias on the aggregate results by making use of the trim-and-fill method (21). After four iterations, nine potential studies were added, and the results showed that publication bias had little impact on the stability of the results (eFigure 1).

Discussion

This meta-analysis involved 35,708 RA patients and 149,421 control subjects. Our meta-analysis of 29 selected papers on the association between RA and thyroid dysfunction revealed a strong correlation between RA and thyroid dysfunction, particularly hypothyroidism. This yielded a pooled OR of 2.25 (95%CI 1.78–2.84) (Figure 2). In addition, RA was closely associated with subclinical hypothyroidism as well as subclinical hyperthyroidism, with ORs of 2.18 (95%CI 1.32–3.61) and 2.13 (95%CI 1.25–3.63), respectively, suggesting an increased risk of both thyroid dysfunctions more than doubled. The OR for hyperthyroidism among RA patients was 1.65 (95%CI, 1.24–2.19).

RA is an autoimmune inflammatory disease affecting symmetrical joints, usually characterized by persistent joint pain, tenderness, and joint destruction (47), while thyroid dysfunction is one of the most common chronic endocrine diseases (48, 49). Thyroid dysfunction mainly discussed in this meta-analysis included hyperthyroidism (low thyroid stimulating hormone [TSH], high free thyroxine [FT4]), hypothyroidism (high TSH, low FT4), subclinical hyperthyroidism (low TSH, normal FT4), subclinical hypothyroidism (high TSH, normal FT4) (49, 50). Because the clinical manifestations of thyroid disease are highly variable and usually nonspecific, the diagnosis of thyroid dysfunction relies primarily on laboratory biochemical tests (49). Although many studies have assessed the risk of thyroid dysfunction in patients with RA, their conclusions have been inconsistent, perhaps due to differences in the populations and geographical regions across studies (51). For example, Nazary et al. compared thyroid function in 400 newly diagnosed RA patients with matched controls and found a higher prevalence of thyroid dysfunction in RA patients (25.25% vs. 11.5%; P value = 0.00001). However, McCoy et al. (14) found that the incidence of hypothyroidism in RA patients was not significantly increased compared with non-RA patients (14). Unfortunately, many studies have not investigated these four types of thyroid dysfunction. Therefore, given the accumulation of evidence, we here performed a systematic review and meta-analysis to address whether patients with RA are at higher risk of developing the four major thyroid disorders.

We performed a meta-analysis of the 29 included articles. Using the four types of thyroid dysfunction as categorical variables, the risk of thyroid dysfunction in RA patients increased to varying degrees, but the comprehensive analysis results also showed high heterogeneity. Subgroup analysis showed that the origin of control groups was the source of heterogeneity in the risk of hyperthyroidism, and study type was the source of heterogeneity for that of subclinical hypothyroidism in the correlation analysis.

When the control group consisted of healthy individuals, RA patients had a higher risk of hyperthyroidism compared with other non-RA patients. This suggests the need for further research into the relationship between thyroid dysfunction and other diseases. Such as the use of chemotherapy drugs as a systemic treatment for breast cancer increases the risk of thyroid dysfunction (52). In turn, a meta-analysis two years ago suggested that hyperthyroidism was associated with prostate, thyroid and breast cancers, compared with those with normal thyroid function (53). A Mendelian randomization study last year also provided further evidence of a causal link between thyroid disease and breast cancer risk (54). Subgroup analyses revealed that geographic region and year of publication were not the main sources of heterogeneity. However, subgroup analysis found that only studies conducted in Asia showed statistically significant differences between the RA group and the control group in terms of the prevalence of hyperthyroidism and subclinical hypothyroidism. In contrast, a subgroup analysis of hyperthyroidism and subclinical hyperthyroidism showed that the difference in the prevalence of these disorders between the RA and control groups was only statistically significant in the last decade compared with earlier studies. We speculated that this finding might be due to variations in detection device sensitivity, diagnostic thresholds, and iodine intake (Table 2) (49). The specific reasons need to be explored in future studies. In the sensitivity analysis of the primary variable, the combined results after excluding studies showed that the risk of developing thyroid dysfunction in RA patients remained significant, indicating that our findings were robust. However, in the sensitivity analysis of hyperthyroidism, after excluding the studies of Mahagna et al. (13) and Huang et al. (25), the combined OR value fluctuated greatly. Moreover, after excluding these two articles, heterogeneity was significantly reduced (P = 0.186, I2 = 30.3%). Therefore, they may be one of the main sources of heterogeneity. The retrospective cohort study by Huang et al. focused on the risk of hearing loss in patients with RA, not thyroid dysfunction. The relevant information could only be extracted from their baseline data. In addition, Huang et al.’s study included only newly diagnosed RA patients and excluded patients with a history of hearing loss, which may have confounded the data on the prevalence of thyroid dysfunction. Mahagna et al. used a cross-sectional study with inherent limitations and low methodological quality scores. Both studies used large databases and had large sample sizes compared to the other studies. However, irrespective of whether these two studies were included, the combined OR values suggested an increased risk of hyperthyroidism in RA patients. To provide more convincing evidence, more prospective studies are needed to more accurately assess the association between RA and thyroid disease. Although the funnel plot and Egger’s test indicated that there was indeed publication bias in meta-analysis of RA association with hypothyroidism, our combined results were proven to be robust by the trim-and-fill method, and that reporting bias did not significantly affect the results.

The underlying mechanisms of thyroid dysfunction in patients with RA are unclear (51). Thyroid dysfunction often occurs in autoimmune thyroid disease (AITD), particularly in Hashimoto’s thyroiditis and Graves’s disease. Moreover, it is now generally accepted that RA and AITD are caused by multiple factors, such as genetic susceptibility and environmental factors (55, 56). We found that patients with RA had a higher risk of developing hypothyroidism than other types of thyroid dysfunction. We hypothesized that the immunogenicity related to hyperthyroidism was lower than that related to hypothyroidism. The vast majority of hypothyroidism is caused by autoimmune diseases, such as Hashimoto’s thyroiditis, whereas autoimmune diseases lead to hyperthyroidism in relatively few cases and are influenced by age and environmental factors (57). A study using Mendelian randomization found a bidirectional causal relationship between RA and Graves’ disease, and good treatment and management of one may benefit the other (58). It is reasonable to speculate that thyroid dysfunction shares a common pathogenesis with RA (59). Exploring the co-pathogenesis pathway of these two diseases and appropriate application of drugs that can act on the co-pathogenesis targets may be an ideal way to prevent and treat co-diseases in the future (58, 60). In addition, a study by Raterman et al. also emphasized the need for careful drug selection in patients with RA complicated by autoimmune thyroiditis. They found that, in these patients, treatment with rituximab and l-thyroxine resulted in a transition to hyperthyroidism several months later (61).

This systematic review and meta-analysis still have the following limitations: First, the quality of the included studies was uneven. According to the NOS scale, 86.21% of the included studies were of medium and high quality, and 13.79% were of low quality. In order to include articles that were as fully compliant as possible, we had not excluded meeting abstracts and letters that provided important data, but which were not always sufficient in their description of detail. Most of these low-quality articles were meeting abstracts without available full text. This also affects, to some extent, the exploration of whether some important variables, such as the duration of RA and general sex distribution, are the main sources of heterogeneity. We performed a subgroup analysis of studies involving subclinical hyperthyroidism based on the mean age of RA patients and found that patients > 50 years had a significantly increased risk of having subclinical hyperthyroidism as compared with RA patients < 50 years (eFigure 2). The relationship between age and the occurrence of other types of thyroid dysfunction needs to be analyzed with more complete data.

Second, all included studies were observational studies, and the results of meta-analysis (62). However, sensitivity analyses indicated that our overall results were robust. Furthermore, using the trim-and-fill method, we confirmed that the combined results were not significantly affected by publication bias. Furthermore, due to the inclusion of observational studies, it is difficult to establish a causal relationship between RA and thyroid dysfunction, which needs to be confirmed by more prospective studies in the future.

Third, the results of meta-analyses with high heterogeneity should be interpreted with caution. Significant heterogeneity still could not be explained by subgroup analyses. Factors such as sex ratio, age, course of RA disease and other factors that were not statistically analyzed due to insufficient data may significantly influence the degree of risk of thyroid dysfunction in RA patients. In fact, gender and age are decisive risk factors for thyroid autoimmunity (15).

Fourth, the clinical criteria for assessing RA and thyroid dysfunction are inconsistent and should be addressed in future studies. In addition, because autoimmune diseases often co-occur, we excluded patients with SLE as controls, but the controls in this study were not all healthy individuals. The association between thyroid disease and SLE has been reported for over 50 years. We believe that using SLE as a control would have confounded the pooled results (17).

Fifth, a small number of reports were not included due to inability to find abstracts or full texts, and the results may have been partially affected. In addition, this study only discussed the risk of thyroid dysfunction in RA patients, and did not consider the incidence of thyroid diseases, such as Hashimoto’s thyroiditis and Graves’s disease. There are many causes of the four types of thyroid dysfunction, both primary and secondary, and their classification is complex. The lack of restriction on etiology is also part of the limitation of our meta-analysis. However, thyroid dysfunction is most common in autoimmune thyroid diseases, especially Graves’ disease and Hashimoto’s thyroiditis. Future studies are needed to more accurately compare the relationship between different thyroid diseases and RA (50, 57).

Our meta-analysis confirmed that RA patients have a high incidence of thyroid dysfunction, and routine biochemical testing of thyroid function in RA patients is reasonable. This not only enables the potential thyroid diseases of RA patients to be detected and treated in a timely manner, but more importantly, the therapeutic drugs for RA may aggravate the thyroid diseases, and the detection of thyroid function is helpful to clinicians for rational drug use (5). Common drugs used to treat RA, such as glucocorticoids (63) and leflunomide (64), interfere with thyroid function, whereas tumor necrosis factor inhibitors improve thyroid function in RA patients with hypothyroidism (65). It is suggested that the effects of various drugs related to the treatment of RA on thyroid function and the influence of thyroid dysfunction on the early treatment effect of RA should be studied more comprehensively (66). Compared with RA patients with negative thyroid autoantibodies, RA patients with positive thyroid autoantibodies had younger onset age and higher disease activity (67). Therefore, thyroid function should be tested regularly in younger RA patients. The close association between thyroid dysfunction and RA suggests a common pathogenic mechanism and potential therapeutic targets, which deserve further exploration.

Conclusion

This meta-analysis revealed that RA patients have a higher prevalence of thyroid dysfunction, especially hypothyroidism. We conclude that it is reasonable to promote routine biochemical testing of thyroid function in patients with RA regardless of the presence or absence of obvious suspicious clinical manifestations, and vice versa. Our findings suggest that rheumatologists should screen RA patients for thyroid function and refer to endocrinologists to determine effective strategies for preventing and treating thyroid dysfunction. Furthermore, since RA may share a common pathogenesis with most AITDs that cause thyroid dysfunction, there may also be corresponding treatment options. More potential mechanisms and treatment options should be explored in future studies.

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 authors.

Author contributions

YJL, ZC, and SL contributed to the conception and design of the systemic review. YJL and H-BM drafted and finalized the manuscript. ZC, H-BM, and SL revised the manuscript and provided critical advice on the content of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the Science and Technology Bureau of Quanzhou [grant number 2020CT003, 2017Z009]; the Natural Science Foundation of Fujian Province [grant number 2020J01219].

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.2022.1015516/full#supplementary-material

Abbreviations

CI, confidence interval; NOS, Newcastle–Ottawa Scale; OR, odds ratio; RA, rheumatoid arthritis; SLE, systemic lupus erythematosus; AITD, autoimmune thyroid disease; TSH, thyroid-stimulating hormone; T3, triiodothyronine; T4, thyroxine; fT3, free triiodothyronine; fT4, free thyroxine.

References

1. Scott DL, Wolfe F, Huizinga TW. Rheumatoid arthritis. Lancet (2010) 376(9746):1094–108. doi: 10.1016/S0140-6736(10)60826-4

CrossRef Full Text | Google Scholar

2. Li Q, Laumonnier Y, Syrovets T, Simmet T. Yeast two-hybrid screening of proteins interacting with plasmin receptor subunit: C-terminal fragment of annexin A2. Acta Pharmacol Sin (2011) 32(11):1411–8. doi: 10.1038/aps.2011.121

CrossRef Full Text | Google Scholar

3. Ozcan M, Karakus MF, Gunduz OH, Tuncel U, Sahin H. Hearing loss and middle ear involvement in rheumatoid arthritis. Rheumatol Int (2002) 22(1):16–9. doi: 10.1007/s00296-002-0185-z

CrossRef Full Text | Google Scholar

4. Almutairi K, Nossent J, Preen D, Keen H, Inderjeeth C. The global prevalence of rheumatoid arthritis: a meta-analysis based on a systematic review. Rheumatol Int (2021) 41(5):863–77. doi: 10.1007/s00296-020-04731-0

CrossRef Full Text | Google Scholar

5. Nazary K, Hussain N, Ojo RO, Anwar S, Kadurei F, Hafizyar F, et al. Prevalence of thyroid dysfunction in newly diagnosed rheumatoid arthritis patients. Cureus. (2021) 13(9):e18204. doi: 10.7759/cureus.18204

CrossRef Full Text | Google Scholar

6. Collins J, Gough S. Autoimmunity in thyroid disease. Eur J Nucl Med Mol Imaging (2002) 29 Suppl 2:S417–24. doi: 10.1007/s00259-002-0848-8

CrossRef Full Text | Google Scholar

7. Conigliaro P, D'Antonio A, Pinto S, Chimenti MS, Triggianese P, Rotondi M, et al. Autoimmune thyroid disorders and rheumatoid arthritis: A bidirectional interplay. Autoimmun Rev (2020) 19(6):102529. doi: 10.1016/j.autrev.2020.102529

CrossRef Full Text | Google Scholar

8. Sapin R, Schlienger JL. [Thyroxine (T4) and tri-iodothyronine (T3) determinations: techniques and value in the assessment of thyroid function]. Ann Biol Clin (Paris) (2003) 61(4):411–20.

Google Scholar

9. Chaker L, Razvi S, Bensenor IM, Azizi F, Pearce EN, Peeters RP. Hypothyroidism. Nat Rev Dis Primers (2022) 8(1):31. doi: 10.1038/s41572-022-00364-8

CrossRef Full Text | Google Scholar

10. Torpy JM, Lynm C, Golub RM. JAMA patient page. hyperthyroidism. JAMA (2011) 306(3):330. doi: 10.1001/jama.306.3.330

CrossRef Full Text | Google Scholar

11. Li Q, Wang B, Mu K, Zhang J, Yang Y, Yao W, et al. Increased risk of thyroid dysfunction among patients with rheumatoid arthritis. Front Endocrinol (2019) 9:799. doi: 10.3389/fendo.2018.00799

CrossRef Full Text | Google Scholar

12. Chaker L, Bianco AC, Jonklaas J, Peeters RP. Hypothyroidism. Lancet (2017) 390(10101):1550–62. doi: 10.1016/S0140-6736(17)30703-1

CrossRef Full Text | Google Scholar

13. Mahagna H, Caplan A, Watad A, Bragazzi NL, Sharif K, Tiosano S, et al. Rheumatoid arthritis and thyroid dysfunction: A cross-sectional study and a review of the literature. Best Pract Res Clin Rheumatol (2018) 32(5):683–91. doi: 10.1016/j.berh.2019.01.021

CrossRef Full Text | Google Scholar

14. McCoy SS, Crowson CS, Gabriel SE, Matteson EL. Hypothyroidism as a risk factor for development of cardiovascular disease in patients with rheumatoid arthritis. J Rheumatol (2012) 39(5):954–8. doi: 10.3899/jrheum.111076

CrossRef Full Text | Google Scholar

15. Ferrari SM, Fallahi P, Ruffilli I, Elia G, Ragusa F, Benvenga S, et al. The association of other autoimmune diseases in patients with graves' disease (with or without ophthalmopathy): Review of the literature and report of a large series. Autoimmun Rev (2019) 18(3):287–92. doi: 10.1016/j.autrev.2018.10.001

CrossRef Full Text | Google Scholar

16. Page MJ, McKenzie JE, Bossuyt PM, Boutron I, Hoffmann TC, Mulrow CD, et al. The PRISMA 2020 statement: An updated guideline for reporting systematic reviews. PloS Med (2021) 18(3):e1003583. doi: 10.1371/journal.pmed.1003583

CrossRef Full Text | Google Scholar

17. Klionsky Y, Antonelli M. Thyroid disease in lupus: An updated review. ACR Open Rheumatol (2020) 2(2):74–8. doi: 10.1002/acr2.11105

CrossRef Full Text | Google Scholar

18. Stang A. Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses. Eur J Epidemiol (2010) 25(9):603–5. doi: 10.1007/s10654-010-9491-z

CrossRef Full Text | Google Scholar

19. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ (Clinical Res ed). (2003) 327(7414):557–60. doi: 10.1136/bmj.327.7414.557

CrossRef Full Text | Google Scholar

20. Qiao Y, Wang Z, Li Y, Han Y, Zhou Y, Cao X. Rheumatoid arthritis risk in periodontitis patients: A systematic review and meta-analysis. Joint Bone Spine (2020) 87(6):556–64. doi: 10.1016/j.jbspin.2020.04.024

CrossRef Full Text | Google Scholar

21. Duval S, Tweedie R. Trim and fill: A simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics (2000) 56(2):455–63. doi: 10.1111/j.0006-341X.2000.00455.x

CrossRef Full Text | Google Scholar

22. Deseatnicova E, Frunze V, Groppa L. Bone and thyroid gland pathology in rheumatoid arthritis. Osteoporosis Int (2022) 32(SUPPL 1):S396. doi: 10.1007/s00198-021-06125-9

CrossRef Full Text | Google Scholar

23. Yadav B, Roy A, Fatima S. A cross-sectional study of different rheumatic diseases and their respective comorbidities at a tertiary care hospital in India. Indian J Rheumatol (2019) 14(1):42–8.

Google Scholar

24. Saqre IM, El-Bahnasawy AS, Farag SEDM, Bazeed FB. Autoimmune thyroid disease in Egyptian patients with rheumatoid arthritis. Egyptian Rheumatol (2019) 41(3):167–71. doi: 10.1016/j.ejr.2018.08.002

CrossRef Full Text | Google Scholar

25. Huang CM, Chen HJ, Huang PH, Tsay GJ, Lan JL, Sung FC. Retrospective cohort study on risk of hearing loss in patients with rheumatoid arthritis using claims data. BMJ Open (2018) 8(1):e018134. doi: 10.1136/bmjopen-2017-018134

CrossRef Full Text | Google Scholar

26. Guy A, Sharif K, Bragazzi NL, Krosser A, Gilburd B, Zeruya E, et al. Low levels of renin and high aldosterone-to-renin ratio among rheumatoid patients and ankylosing spondylitis patients: A prospective study. Israel Med Assoc J (2018) 20(10):632–6.

Google Scholar

27. Figueroa-Sánchez M, Núñez-Atahualpa L, Hernández-Zúñiga M, Martín-Marquez BT, Martínez-García EA, Gómez-Bañuelos E, et al. Clinical and ultrasonographic evaluation of the thyroid gland in patients with rheumatoid arthritis. Gaceta Med Mexico (2018) 154(4):432–7. doi: 10.24875/GMM.18003556

CrossRef Full Text | Google Scholar

28. El Achek MA. Are co-morbidities during rheumatoid arthritis different from those in ankylosing spondyloarthritis? Ann Rheumatic Dis (2018) 77:1358. doi: 10.1136/annrheumdis-2018-eular.5788

CrossRef Full Text | Google Scholar

29. Posselt RT, Coelho VN, Pigozzo DC, Guerrer MI, Fagundes MC, Nisihara R, et al. Prevalence of thyroid autoantibodies in patients with systematic autoimmune rheumatic diseases. Cross-sectional study. Sao Paulo Med J (2017) 135(6):535–40. doi: 10.1590/1516-3180.2017.0089110617

CrossRef Full Text | Google Scholar

30. Tascilar K, Dell'Aniello S, Hudson M, Suissa S. Statins and risk of rheumatoid arthritis: A nested case-control study. Arthritis Rheumatol (2016) 68(11):2603–11. doi: 10.1002/art.39774

CrossRef Full Text | Google Scholar

31. El-saadany H, Abd Elkhalik M, Moustafa T, Abd El bar E. Thyroid dysfunction in systemic lupus erythematosus and rheumatoid arthritis: Its impact as a cardiovascular risk factor. Egyptian Rheumatol (2014) 36(2):71–8. doi: 10.1016/j.ejr.2013.12.003

CrossRef Full Text | Google Scholar

32. Gomez R, Gutierrez JM, Fernandez-Avila DG, Diaz MC, Aschner P. Autoinmune thyroid disease in patients with rheumatoid arthritis: Case-control study. Ann Rheumatic Dis (2013) 71:170. doi: 10.1136/annrheumdis-2012-eular.2016

CrossRef Full Text | Google Scholar

33. Acay A, Ulu MS, Ahsen A, Eroglu S, Ozuguz U, Yuksel S, et al. Assessment of thyroid disorders and autoimmunity in patients with rheumatic diseases. Endocrine Metab Immune Disord Drug Targets (2014) 14(3):182–6. doi: 10.2174/1871530314666140626113111

CrossRef Full Text | Google Scholar

34. Mousa AA, Ghonem M, Hegazy A, El-Baiomy AA, El-Diasty A. Thyroid function and auto-antibodies in Egyptian patients with systemic lupus erythematosus and rheumatoid arthritis. Trends Med Res (2012) 7(1):25–33. doi: 10.3923/tmr.2012.25.33

CrossRef Full Text | Google Scholar

35. Fatima F, Rao Uppuluri RK, Moots R, Goodson N. Prevalence of co-morbid conditions in patients with rheumatoid arthritis: Case control study. Int J rheumatic diseases. (2010) 13:66. doi: 10.1111/j.1756-185X.2010.01496.x

CrossRef Full Text | Google Scholar

36. Mobini M, Kashi Z, Ravanbakhsh N. Thyroid disorders in rheumatoid arthritis and osteoarthritis. Pakistan J Med Sci (2011) 27(3):595–8.

Google Scholar

37. Przygodzka M, Filipowicz-Sosnowska A. Prevalence of thyroid diseases and antithyroid antibodies in women with rheumatoid arthritis. Polskie Archiwum Medycyny Wewnetrznej (2009) 119(1-2):39–43. doi: 10.20452/pamw.600

CrossRef Full Text | Google Scholar

38. Haghighi A. Prevalence of thyroid disorders in systemic lupus erythematosus and rheumatoid arthritis. Iranian J Endocrinol Metab (2009) 11(1).

Google Scholar

39. Al-Awadhi AM, Olusi S, Hasan EA, Abdullah A. Frequency of abnormal thyroid function tests in Kuwaiti arabs with autoimmune diseases. Med Principles Practice. (2008) 17(1):61–5. doi: 10.1159/000109592

CrossRef Full Text | Google Scholar

40. Andonopoulos AP, Siambi V, Makri M, Christofidou M, Markou C, Vagenakis AG. Thyroid function and immune profile in rheumatoid arthritis. A Controlled study. Clin Rheumatol (1996) 15(6):599–603. doi: 10.1007/BF02238551

CrossRef Full Text | Google Scholar

41. Shiroky JB, Cohen M, Ballachey ML, Neville C. Thyroid dysfunction in rheumatoid arthritis: A controlled prospective survey. Ann Rheumatic Diseases (1993) 52(6):454–6. doi: 10.1136/ard.52.6.454

CrossRef Full Text | Google Scholar

42. Antonelli A, Delle Sedie A, Fallahi P, Ferrari SM, Maccheroni M, Ferrannini E, et al. High prevalence of thyroid autoimmunity and hypothyroidism in patients with psoriatic arthritis. J Rheumatol (2006) 33(10):2026–8.

Google Scholar

43. Al-Awadhi AM, Olusi SO, Al-Zaid NS, Prabha K, Al-Ali N, Al-Jarallah KH, et al. Prevalence of hypothyroidism among arabs with rheumatoid arthritis [2]. Rheumatology. (1999) 38(5):472–3. doi: 10.1093/rheumatology/38.5.472

CrossRef Full Text | Google Scholar

44. Innocencio RM, Romaldini JH, Ward LS. Thyroid autoantibodies in autoimmune diseases. Medicina. (2004) 64(3):227–30.

Google Scholar

45. Kostić I, Živančević-Simonović S, Bukilica M, Dimitrijević L. Thyroid function and antithyroid autoantibodies in patients with connective tissue diseases. Medicus. (2006) 7(2):61–4.

Google Scholar

46. Santos MJ, Vinagre F, Da Silva JC, Gil V, Fonseca JE. Cardiovascular risk profile in systemic lupus erythematosus and rheumatoid arthritis: A comparative study of female patients. Acta reumatol portuguesa. (2010) 35(3):325–32. doi: 10.1186/1471-2474-11-144

CrossRef Full Text | Google Scholar

47. Padyukov L. Genetics of rheumatoid arthritis. Semin Immunopathol (2022) 44(1):47–62. doi: 10.1007/s00281-022-00912-0

CrossRef Full Text | Google Scholar

48. Biondi B, Kahaly GJ, Robertson RP. Thyroid dysfunction and diabetes mellitus: Two closely associated disorders. Endocr Rev (2019) 40(3):789–824. doi: 10.1210/er.2018-00163

CrossRef Full Text | Google Scholar

49. Taylor PN, Albrecht D, Scholz A, Gutierrez-Buey G, Lazarus JH, Dayan CM, et al. Global epidemiology of hyperthyroidism and hypothyroidism. Nat Rev Endocrinol (2018) 14(5):301–16. doi: 10.1038/nrendo.2018.18

CrossRef Full Text | Google Scholar

50. Chaker L, Razvi S, Bensenor IM, Azizi F, Pearce EN, Peeters RP. Hypothyroidism. Nat Rev Dis Primers (2022) 8(1):30. doi: 10.1038/s41572-022-00357-7

CrossRef Full Text | Google Scholar

51. Yin R, Li L, Xu L, Sui W, Niu M, Xu R, et al. Association between depression and sleep quality in patients with systemic lupus erythematosus: a systematic review and meta-analysis. Sleep Breath (2022) 26(1):429–41. doi: 10.1007/s11325-021-02405-0

CrossRef Full Text | Google Scholar

52. Mortezaee K, Ahmadi A, Haghi-Aminjan H, Khanlarkhani N, Salehi E, Shabani Nashtaei M, et al. Thyroid function following breast cancer chemotherapy: A systematic review. J Cell Biochem (2019) 120(8):12101–7. doi: 10.1002/jcb.28771

CrossRef Full Text | Google Scholar

53. Tran TV, Kitahara CM, de Vathaire F, Boutron-Ruault MC, Journy N. Thyroid dysfunction and cancer incidence: A systematic review and meta-analysis. Endocr Relat Cancer. (2020) 27(4):245–59. doi: 10.1530/ERC-19-0417

CrossRef Full Text | Google Scholar

54. Yuan S, Kar S, Vithayathil M, Carter P, Mason AM, Burgess S, et al. Causal associations of thyroid function and dysfunction with overall, breast and thyroid cancer: A two-sample mendelian randomization study. Int J Cancer. (2020) 147(7):1895–903. doi: 10.1002/ijc.32988

CrossRef Full Text | Google Scholar

55. McInnes IB, Schett G. The pathogenesis of rheumatoid arthritis. N Engl J Med (2011) 365(23):2205–19. doi: 10.1056/NEJMra1004965

CrossRef Full Text | Google Scholar

56. Weetman AP. Diseases associated with thyroid autoimmunity: Explanations for the expanding spectrum. Clin Endocrinol (2011) 74(4):411–8. doi: 10.1111/j.1365-2265.2010.03855.x

CrossRef Full Text | Google Scholar

57. Cooper DS. Hyperthyroidism. Lancet. (2003) 362(9382):459–68. doi: 10.1016/S0140-6736(03)14073-1

CrossRef Full Text | Google Scholar

58. Wu D, Xian W, Hong S, Liu B, Xiao H, Li Y. Graves' disease and rheumatoid arthritis: A bidirectional mendelian randomization study. Front Endocrinol (2021) 12:702482. doi: 10.3389/fendo.2021.702482

CrossRef Full Text | Google Scholar

59. Punzi L, Betterle C. Chronic autoimmune thyroiditis and rheumatic manifestations. Joint Bone Spine (2004) 71(4):275–83. doi: 10.1016/j.jbspin.2003.06.005

CrossRef Full Text | Google Scholar

60. Guo Y, Li H, Chen X, Yang H, Guan H, He X, et al. Novel roles of chloroquine and hydroxychloroquine in graves' orbitopathy therapy by targeting orbital fibroblasts. J Clin Endocrinol Metab (2020) 105(6): 1906–17. doi: 10.1210/clinem/dgaa161

CrossRef Full Text | Google Scholar

61. Raterman HG, Simsek S, Lems WF, Meesters EW, Dijkmans BA, Nurmohamed MT. Rituximab and thyroid function. Arch Internal Med (2009) 169(11):1073–4. doi: 10.1001/archinternmed.2009.116

CrossRef Full Text | Google Scholar

62. Egger M, Schneider M, Davey Smith G. Spurious precision? meta-analysis of observational studies. BMJ (Clinical Res ed). (1998) 316(7125):140–4. doi: 10.1136/bmj.316.7125.140

CrossRef Full Text | Google Scholar

63. Re RN, Kourides IA, Ridgway EC, Weintraub BD, Maloof F. The effect of glucocorticoid administration on human pituitary secretion of thyrotropin and prolactin. J Clin Endocrinol Metab (1976) 43(2):338–46. doi: 10.1210/jcem-43-2-338

CrossRef Full Text | Google Scholar

64. Grzywa M, Kadziela M. [The conversion of hypothyroidism into hyperthyroidism during leflunomide with povidone iodine treatment of rheumatoid arthritis]. Polski merkuriusz lekarski Organ Polskiego Towarzystwa Lekarskiego. (2013) 34(204):348–50.

Google Scholar

65. Bliddal S, Borresen SW, Feldt-Rasmussen U. Thyroid autoimmunity and function after treatment with biological antirheumatic agents in rheumatoid arthritis. Front endocrinol (2017) 8:179. doi: 10.3389/fendo.2017.00179

CrossRef Full Text | Google Scholar

66. Waldenlind K, Delcoigne B, Saevarsdottir S, Askling J. Does autoimmune thyroid disease affect rheumatoid arthritis disease activity or response to methotrexate? RMD Open (2020) 6(2):e001282. doi: 10.1136/rmdopen-2020-001282

CrossRef Full Text | Google Scholar

67. Chen YL, Lin JZ, Mo YQ, Liang JJ, Li QH, Zhou CJ, et al. Joint damage is amplified in rheumatoid arthritis patients with positive thyroid autoantibodies. PeerJ. (2018) 6:e4216. doi: 10.7717/peerj.4216

CrossRef Full Text | Google Scholar

Keywords: rheumatoid arthritis, subclinical hyperthyroidism, subclinical hypothyroidism, thyroid dysfunction, immune system

Citation: Liu Y-j, Miao H-b, Lin S and Chen Z (2022) Association between rheumatoid arthritis and thyroid dysfunction: A meta-analysis and systematic review. Front. Endocrinol. 13:1015516. doi: 10.3389/fendo.2022.1015516

Received: 09 August 2022; Accepted: 23 September 2022;
Published: 13 October 2022.

Edited by:

Sijie Fang, Shanghai Jiao Tong University, China

Reviewed by:

Dilia Giuggioli, University of Modena and Reggio Emilia, Italy
Maryam Masoumi, Qom University of Medical Sciences, Iran

Copyright © 2022 Liu, Miao, Lin and Chen. 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: Shu Lin, shulin1956@126.com; Zhen Chen, zhenchen180@fjmu.edu.cn

ORCID: Shu Lin, orcid.org/0000-0002-4239-2028

Disclaimer: 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.