Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 02 October 2020
Sec. Autoimmune and Autoinflammatory Disorders

Increased Risk of Inflammatory Bowel Disease Among Patients With Ankylosing Spondylitis: A 13-Year Population-Based Cohort Study

  • 1Department of Rheumatology, The First Affiliated Hospital of Harbin Medical University, Harbin, China
  • 2Functional Neurosurgery Division, Neurological Institute, Taichung Veterans General Hospital, Taichung, Taiwan
  • 3Department of Rehabilitation, Jen-Teh Junior College of Medicine, Nursing and Management, Miaoli County, Taiwan
  • 4School of Health Policy and Management, Chung Shan Medical University, Taichung, Taiwan
  • 5Department of Medical Research, Chung Shan Medical University Hospital, Taichung, Taiwan
  • 6Beijing Tsinghua Changgung Hospital, School of Clinical Medicine, Tsinghua University, Beijing, China
  • 7Department of Rheumatology & Immunology, Chung Shan Medical University Hospital, Taichung, Taiwan
  • 8Institute of Medicine, Chung Shan Medical University, Taichung, Taiwan
  • 9Graduate Institute of Integrated Medicine, China Medical University, Taichung, Taiwan

Aim: Ankylosing spondylitis (AS) primarily affects the axial skeleton and extraarticular structures. Small-scaled studies have reported that the incidence and prevalence of inflammatory bowel disease (IBD) are higher in patients with AS than in the general population. This study determined the incidence of IBD in patients with AS using a large scale population-based cohort dataset.

Methods: This was a retrospective cohort study. Patient data were collected from the Taiwan National Health Insurance Research Database from 2000 to 2012. We enrolled 3,804 patients with AS and 7,608 non-AS patients. The endpoint was IBD diagnosis by using International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) coding 555 and 556 after at least three outpatient visits or one hospital admission, until the end of 2012. The Kaplan–Meier analysis was performed to discriminate the cumulative incidence of IBD and the log-rank test was used to test the significance. A Cox proportional hazard model was used to estimate the hazard ratio (HR) for IBD between the AS and non-AS groups.

Results: Among the population as a whole the Cox proportional hazard regression indicated that patients aged ≥65 years [adjusted HR (aHR): 2.48, 95% confidence interval (CI): 1.38–4.47] or with comorbidity of cancer (aHR: 3.51, 95% CI: 1.40–8.80) had a higher HR for IBD. Kaplan–Meier curves of cumulative incidence proportion of IBD indicated that patients with AS had a higher risk of IBD than the non-AS group in the subgroup aged <40 years (HR: 2.85, 95% CI: 1.51–5.40, p = 0.001).

Conclusions: Patients with AS aged <40 years had a higher IBD risk than did those without AS in Taiwan. Clinicians and patients should be aware of this association.

Highlights

Patients with AS aged <40 years had a higher IBD risk than patients without AS in Taiwan.

Introduction

Inflammatory bowel disease (IBD), which includes Crohn’s disease (CD), ulcerative colitis (UC), and indeterminate colitis (IBD-unclassified), is a group of chronic relapsing inflammatory disorder affecting the gastrointestinal (GI) tract caused by a complex interaction of genetic factors, environmental exposures, and immune response dysregulation. CD can trigger inflammation in any part of the GI tract and results in various GI and systemic symptoms, such as abdominal pain, diarrhea, fistulae, weight loss, and anorexia. UC is characterized by inflammation in the colon, which leads to abdominal pain and subsequent diarrhea. However, IBD also displays certain extraintestinal manifestations. For instance, studies have reported an association between IBD and arthritis (1, 2).

Ankylosing spondylitis (AS) is a chronic inflammatory disease that affects the axial skeleton and extraarticular structures, such as eyes, skin, and bowels. IBD incidence and prevalence are higher in patients with psoriasis, psoriatic arthritis, and AS than in the general population.(3) The estimated overall incidence of IBD in patients with AS is 5%–10%, with approximately 60% of patients with AS presenting subclinical ileal inflammation (4). However, Wang et al. reported that the occurrence of IBD in Han Chinese patients with AS was rare.(5) Therefore, we examined the IBD risk in Taiwanese patients with AS using the large Taiwan National Health Insurance Research Database (NHIRD) and compared the risk in patients with AS with that in matched comparison cohort patients without AS.

Materials and Methods

Data Source

In this population-based retrospective cohort study, data were obtained from the NHIRD which enrolls approximately 99% of the 23 million beneficiaries in Taiwan. This database includes all the registrations and claims data, including diagnoses, drug prescriptions, inpatient care, outpatient visits, emergency hospitalization, and diagnoses according to the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM). Data collected from 2000 to 2012 were extracted for 1 million subjects, sampled from the 23 million beneficiaries. The 1 million longitudinal health insurance database was generated by randomly allocation to represent 23 millions whole population, provide by the Authority of National Health Insurance. The sampled database was de-identified as the Longitudinal Health Insurance Database (LHID), and the study was approved by the Institutional Review Board of Chung Shan Medical University Hospital.

Study Group and Outcome Measurement

Patients diagnosed as having AS between January 1, 2000, and December 31, 2012, were selected from the LHID (n = 4007). Patients receiving their first diagnosis of AS (ICD-9-CM code 720.0) after at least three outpatient visits or one hospitalization, were considered to be newly diagnosed as having AS. The date of the first AS diagnosis was set as the index date. Patients diagnosed as having inflammatory bowel disease (ICD-9-CM 555 and 556) before the index date were excluded (n = 203). The non-AS group comprised 988084 individuals who were never diagnosed as having AS, randomly selected from 1999 to 2013. A ratio of 1:8 was used to match the 30,432 comparison cohort to patients with AS based on age and sex. The index date for the non-AS group was set as the index date of their AS group counterpart.

Follow-up was begun on the index date and ended on the IBD diagnosis date, December 31, 2012, or the date of withdrawal from the national insurance system. New cases of IBD were identified from the database by the presence of records of ICD-9-CM code 555 or 556 with at least three outpatient visits or one hospitalization.

Covariates and Matching

To reduce selection bias, the AS group was 1:2 matched with the non-AS group by propensity scores based on covariates and baseline characteristics, including age, sex, hypertension (ICD-9-CM 401–405), hyperlipidemia (ICD-9-CM 272.0–272.4), diabetes (ICD-9-CM 250), chronic obstructive pulmonary disease (COPD; ICD-9-CM 490–496), chronic kidney disease (ICD-9-CM 585), cardiovascular disease (ICD-9-CM 410–414), stroke (ICD-9-CM 430–438), cancer (ICD-9-CM 140–208), allergic rhinitis (ICD-9-CM 477), and obesity (ICD-9-CM 278). These comorbidities were defined within 1 year before the index date based on at least three outpatient visits or one hospitalization and were considered covariates in the multivariate analysis. The AS group comprised 3,804 patients and the non-AS group comprised 7,608 patients. Figure 1 displays the flowchart of the sample selection.

FIGURE 1
www.frontiersin.org

Figure 1 Flowchart of the sample selection.

Statistical Analysis

The chi-square or Fisher exact test were performed to compare categorical variables of the demographic data of AS and non-AS groups, and Student t tests were performed to compare continuous variables. A Kaplan–Meier analysis was performed to assess the cumulative incidence of IBD and a log-rank test was used to test the significance. A Cox proportional hazard model was used to estimate the hazard ratio (HR) for IBD between the AS and non-AS groups. Subgroup analyses were also performed by age and sex. The statistical software used was SPSS (version 18.0; SPSS Inc., Chicago, IL, USA). Statistical significance was defined as a p of <0.05.

Results

Study Population

We identified 4,007 patients with AS with at least three outpatient visits or one hospital admission for AS. We further excluded patients diagnosed as having IBD before to the index date. A total of 11,412 patients were recruited from the claims database, with 3,804 in the AS group and 7,608 in the non-AS group (Figure 1).

Baseline Characteristics of Patients With AS

Among the patients with AS, 38.4% were women and 61.6% were men. Before propensity score matching, the AS group were more likely to have comorbid hypertension (p < 0.001), hyperlipidemia (p = 0.001), diabetes (p = 0.025), COPD (p < 0.001), cardiovascular disease (p < 0.001), stroke (p = 0.023), and allergic rhinitis (p < 0.001) than the non-AS group. By contrast, after propensity score matching, there were no differences in baseline characteristics between the AS and non-AS group, including age, sex, comorbidity hypertension, hyperlipidemia, diabetes, COPD, Chronic kidney disease, cardiovascular disease, stroke, cancer, allergic rhinitis, and obesity (Table 1).

TABLE 1
www.frontiersin.org

Table 1 Demographic characteristics of Ankylosing Spondylitis and Non-Ankylosing Spondylitis.

HR for IBD

Over the 13 years of follow-up period, 61 patients without and 42 patients with AS were diagnosed as having IBD. This incidence density was 1.0 and 1.4 per 1,000 person-year in patients without and with AS, respectively (Table 2). The Cox proportional hazard regression indicated that patients aged ≥65 years [adjusted HR (aHR): 2.48, 95% confidence interval (CI): 1.38–4.47] and with comorbidity of cancer (aHR: 3.51, 95% CI: 1.40–8.80) had a higher HR for IBD. However, no difference was identified in the HR for IBD between the non-AS and AS groups (aHR: 1.37, 95% CI: 0.92–2.02).

TABLE 2
www.frontiersin.org

Table 2 Cox proportional hazard model.

Time-to-Event Analysis and Subgroup Analysis

Kaplan–Meier curves of the cumulative incidence proportion of IBD are displayed in Figure 2; no difference was determined in incidence risk between the AS and non-AS group (log-rank, p = 0.114). The subgroup analyses revealed that patients with AS aged <40 years had a higher HR (2.85, 95% CI: 1.51–5.40) than patients without AS (Table 3). Kaplan–Meier curves of cumulative incidence proportion of IBD (age < 40) are displayed in Figure 3 (p = 0.001).

FIGURE 2
www.frontiersin.org

Figure 2 Kaplan-Meier curve of cumulative incidence proportion of IBD in the AS and the comparison cohort population.

TABLE 3
www.frontiersin.org

Table 3 Subgroup analysis of Cox proportional hazard model.

FIGURE 3
www.frontiersin.org

Figure 3 Kaplan-Meier curve of cumulative incidence proportion of IBD in the AS and the non-AS group whose age <40.

Discussion

In this population-based cohort study, we determined that patients aged ≥65 years (aHR: 2.48, 95% CI: 1.38–4.47) and patients with cancer (aHR: 3.51, 95% CI: 1.40–8.80) had a higher HR for IBD among the population as a whole. Patients with AS had a higher risk of IBD than non-AS in the subgroup aged <40 years (HR: 2.85, 95% CI: 1.51–5.40, p = 0.001).

A study reported that patients with AS had a higher risk than the general population of developing hypertension, hyperlipidemia, diabetes, cancer, and other cardiovascular, pulmonary, neurological, and renal complications.(6, 7) Concurrently, our findings indicated that patients with AS in Taiwan were more likely to have comorbid hypertension, hyperlipidemia, diabetes, COPD, cardiovascular disease, stroke, and allergic rhinitis at baseline (Table 1). Our results showed that, patients with cancer had a higher risk of IBD. It was reported that, the underlying inflammation generated from IBD can promote initate and promote cancer development. Immune cells infiltrate in colon lead to reactive oxygen species and reactive nitrogen species release which cause damage of DNA then initiate cancer. Moreover, inflammatory celss produce growth factors that drive tumor progression.(8, 9)

The highest prevalence of IBD has been reported in Europe (UC 5.1 per 1,000 persons in the southeast of Norway; CD 3.2 per 1,000 persons in Hesse, Germany) and North America (UC 2.8 per 1000 persons in Olmsted County, USA; CD 3.2 per 1,000 persons in Nova Scotia, Canada).(10) We identified a lower incidence rate of IBD in Taiwan; the incidence rate of IBD was 1.0 per 1,000 person-years in patients without AS and 1.4 per 1,000 person-years in Taiwan patients with AS. This may be because of ethnic and environmental differences.

CD and UC reportedly have a bimodal age distribution, with a first peak at approximately 25 years of age and a second peak after 60 years of age, according to population-based studies and hospitalization statistics.(1115) Our findings were similar. The Cox proportional hazard regression indicated that patients aged ≥65 years had a higher HR for IBD (aHR: 3.14, 95% CI: 1.04–9.44, p < 0.05; Table 2). The Cox proportional hazard analysis also indicated that the HR for IBD was higher in patients with comorbid cancer (p < 0.05), with an adjusted HR of 3.51 (95% CI: 1.40–8.77; Table 2). This finding is consistent with reports that malignancies such as colorectal cancer, liver cancer, lymphoma, lung cancer, and biliary tract tumor were significantly related to IBD mortality.(1619)

Stolwijk et al. determined that the incidence and prevalence of IBD were higher in patients with AS than in the general population.(3, 20, 21) Approximately 10% of patients with AS have clinical IBD and up to 70% have subclinical bowel inflammation, demonstrated histologically.(22) Our results did not indicate a difference in HRs between AS and non-AS groups. However, the subgroup analysis of the Cox proportional hazard model revealed that patients with AS had a higher risk of IBD than the non-AS group in the subgroup aged <40 years (HR: 2.85, 95% CI: 1.51–5.40; Table 3). Kaplan–Meier curves of the cumulative incidence proportion of IBD indicated that patients with AS aged <40 years had a higher risk of IBD than the non-AS group (p = 0.001), as illustrated in Figure 3.

In this population-based retrospective cohort study, data were obtained from the NHIRD, which is a peer-reviewed long-term population-based database that covers approximately 99% of the 23 million inhabitants of Taiwan. This is the major strength of our study. However, there are some limitations. First, we used ICD-9-CM diagnostic codes, instead of medical records to confirm diagnoses of IBD and AS. This may have led to the inclusion of misdiagnosed or over-diagnosed cases. However, ICD coding is a catastrophic illness diagnosis and is reviewed by expert clinicians according to Taiwan insurance reimbursement policy, which improves the accuracy. Second, some other confounding factors, such as disease activity, smoking habits, alcohol consumption, and exercise habits, were not collected in this study. To minimize this bias, we stratified lifestyle-related diseases, such as COPD, diabetes, hyperlipidemia as a proxy and propensity score-matched participants to improve baseline comparability. Third, NSAIDs may affect the entire gastrointestinal system, and biologics such as secukinumab, golimumab, and etanercept may increase IBD risk.(2325) However, we did not collect the data on drug prescriptions. These confounding factors will be included in the study design of our future works.

In conclusion, this large population-based cohort study demonstrated that patients with AS aged <40 years had a higher IBD risk than patients without AS in Taiwan. Physicians managing patients with AS should consider this association when monitoring patients and prescribing medication.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

This observational human study protocol and informed consent waiver were reviewed and approved by Chung Shan mMedical University Hospital Institutional Review Board, No. CS15134.

Author Contributions

SW, JC-C, and ZZ wrote the manuscript. Y-HW analyzed the data. All authors contributed to the article and approved the submitted version.

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.

Acknowledgments

This study is based in part on data from the National Health Insurance Research Database provided by the Bureau of National Health Insurance, the Department of Health, and managed by the National Health Research Institutes. The authors had permission from the Department of Health and managed to use those data for analysis with full access to all of the data in the study. SW, JW, and ZZ take full responsibility for the integrity of the data and the accuracy of the data analysis.

References

1. Garber A, Regueiro M. Extraintestinal Manifestations of Inflammatory Bowel Disease: Epidemiology, Etiopathogenesis, and Management. Curr Gastroenterol Rep (2019) 21(7):31. doi: 10.1007/s11894-019-0698-1

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Vavricka SR, Schoepfer A, Scharl M, Lakatos PL, Navarini A, Rogler G. Extraintestinal Manifestations of Inflammatory Bowel Disease. Inflammat Bowel Dis (2015) 21(8):1982–92. doi: 10.1097/MIB.0000000000000392

CrossRef Full Text | Google Scholar

3. Stolwijk C, Essers I, van Tubergen A, Boonen A, Bazelier MT, De Bruin ML, et al. The epidemiology of extra-articular manifestations in ankylosing spondylitis: a population-based matched cohort study. Ann Rheumat Dis (2015) 74(7):1373–8. doi: 10.1136/annrheumdis-2014-205253

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Fragoulis GE, Liava C, Daoussis D, Akriviadis E, Garyfallos A, Dimitroulas T. Inflammatory bowel diseases and spondyloarthropathies: From pathogenesis to treatment. World J Gastroenterol (2019) 25(18):2162–76. doi: 10.3748/wjg.v25.i18.2162

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Wang CR, Weng CT, Lee CT, Huang KY, Hsu SM, Liu MF. Rare occurrence of inflammatory bowel disease in a cohort of Han Chinese ankylosing spondylitis patients- a single institute study. Sci Rep (2017) 7(1):13165. doi: 10.1038/s41598-017-13573-z

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Walsh JA, Song X, Kim G, Park Y. Evaluation of the comorbidity burden in patients with ankylosing spondylitis using a large US administrative claims data set. Clin Rheumatol (2018) 37(7):1869–78. doi: 10.1007/s10067-018-4086-2

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Tanidjaja I, Damodaran S. Influence of amino acids on thermal stability and heat-set gelation of bovine serum albumin. Food Chem (2020) 337:127670. doi: 10.1016/j.foodchem.2020.127670

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Hussain SP, Hofseth LJ, Harris CC. Radical causes of cancer. Nat Rev Cancer (2003) 3(4):276–85. doi: 10.1038/nrc1046

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Francescone R, Hou V, Grivennikov SI. Cytokines, IBD, and colitis-associated cancer. Inflammat Bowel Dis (2015) 21(2):409–18. doi: 10.1097/MIB.0000000000000236

CrossRef Full Text | Google Scholar

10. Ng SC, Shi HY, Hamidi N, Underwood FE, Tang W, Benchimol EI, et al. Worldwide incidence and prevalence of inflammatory bowel disease in the 21st century: a systematic review of population-based studies. Lancet (2018) 390(10114):2769–78. doi: 10.1016/S0140-6736(17)32448-0

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Butter M, Weiler S, Biedermann L, Scharl M, Rogler G, Bischoff-Ferrari HA, et al. Clinical manifestations, pathophysiology, treatment and outcome of inflammatory bowel diseases in older people. Maturitas (2018) 110:71–8. doi: 10.1016/j.maturitas.2018.01.015

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Sonnenberg A. Age distribution of IBD hospitalization. Inflammat Bowel Dis (2010) 16(3):452–7. doi: 10.1002/ibd.21058

CrossRef Full Text | Google Scholar

13. Kurata JH, Kantor-Fish S, Frankl H, Godby P, Vadheim CM. Crohn’s disease among ethnic groups in a large health maintenance organization. Gastroenterology (1992) 102(6):1940–8. doi: 10.1016/0016-5085(92)90317-R

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Smith GD, Watson R, Thompson DR. Older people and inflammatory bowel disease: a systematic review. J Clin Nurs (2008) 17(11c):400–6. doi: 10.1111/j.1365-2702.2008.02577.x

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Hruz P, Juillerat P, Kullak-Ublick GA, Schoepfer AM, Mantzaris GJ, Rogler G. Management of the Elderly Inflammatory Bowel Disease Patient. Digestion (2020) 2020:1–15. doi: 10.1159/000503099

CrossRef Full Text | Google Scholar

16. Kim HJ, Hann HJ, Hong SN, Kim KH, Ahn IM, Song JY, et al. Incidence and natural course of inflammatory bowel disease in Korea, 2006-2012: a nationwide population-based study. Inflammat Bowel Dis (2015) 21(3):623–30. doi: 10.1097/MIB.0000000000000313

CrossRef Full Text | Google Scholar

17. Jess T, Loftus EV Jr., Harmsen WS, Zinsmeister AR, Tremaine WJ, Melton LJ 3rd, et al. Survival and cause specific mortality in patients with inflammatory bowel disease: a long term outcome study in Olmsted County, Minnesota, 1940-2004. Gut (2006) 55(9):1248–54. doi: 10.1136/gut.2005.079350

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Jussila A, Virta LJ, Pukkala E, Farkkila MA. Mortality and causes of death in patients with inflammatory bowel disease: a nationwide register study in Finland. J Crohn’s Colitis (2014) 8(9):1088–96. doi: 10.1016/j.crohns.2014.02.015

CrossRef Full Text | Google Scholar

19. Lin WC, Weng MT, Tung CC, Chang YT, Leong YL, Wang YT, et al. Trends and risk factors of mortality analysis in patients with inflammatory bowel disease: a Taiwanese nationwide population-based study. J Trans Med (2019) 17(1):414. doi: 10.1186/s12967-019-02164-3

CrossRef Full Text | Google Scholar

20. De Keyser F, Elewaut D, De Vos M, De Vlam K, Cuvelier C, Mielants H, et al. Bowel inflammation and the spondyloarthropathies. Rheumat Dis Clinics North America (1998) 24(4):785–813, ix-x. doi: 10.1016/S0889-857X(05)70042-9

CrossRef Full Text | Google Scholar

21. De Vos M, Cuvelier C, Mielants H, Veys E, Barbier F, Elewaut A. Ileocolonoscopy in seronegative spondylarthropathy. Gastroenterology (1989) 96(2 Pt 1):339–44. doi: 10.1016/0016-5085(89)91557-6

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Mielants H, Veys EM, Cuvelier C, De Vos M, Botelberghe L. HLA-B27 related arthritis and bowel inflammation. Part 2. Ileocolonoscopy and bowel histology in patients with HLA-B27 related arthritis. J Rheumatol (1985) 12(2):294–8.

PubMed Abstract | Google Scholar

23. Uskudar Cansu D, Uskudar Teke H, Temel T, Erturk A, Kahraman O, Korkmaz C. Do Anti-TNF Agents Increase the Risk of Inflammatory Bowel Disease Evolution in Patients with Ankylosing Spondylitis? Real Life Data. J Natl Med Assoc (2019) 111(3):262–9. doi: 10.1016/j.jnma.2018.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Fobelo Lozano MJ, Serrano Gimenez R, Castro Fernandez M. Emergence of Inflammatory Bowel Disease During Treatment With Secukinumab. J Crohn’s Colitis (2018) 12(9):1131–3. doi: 10.1093/ecco-jcc/jjy063

CrossRef Full Text | Google Scholar

25. Bawany MZ, Rafiq E, Thotakura R, Lay R, Silverman AL, Nawras A. Golimumab may induce exacerbation of inflammatory bowel disease when it is used for the treatment of ankylosing spondylitis: a case report with a review of literature. Am J Ther (2014) 21(1):e26–27. doi: 10.1097/MJT.0b013e31825e6089

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: ankylosing spondylitis, inflammatory bowel disease, hazard ratio, musculoskeletal disorders, rheumatology, inflammatory

Citation: Wang S, Tsou H-K, Chiou J-Y, Wang Y-H, Zhang Z and Wei JC-C (2020) Increased Risk of Inflammatory Bowel Disease Among Patients With Ankylosing Spondylitis: A 13-Year Population-Based Cohort Study. Front. Immunol. 11:578732. doi: 10.3389/fimmu.2020.578732

Received: 01 July 2020; Accepted: 14 September 2020;
Published: 02 October 2020.

Edited by:

Antonio Vitale, University of Siena, Italy

Reviewed by:

Matthew Stoll, University of Alabama at Birmingham, United States
Malcolm Koo, Tzu Chi University of Science and Technology, Taiwan

Copyright © 2020 Wang, Tsou, Chiou, Wang, Zhang and Wei. 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: Zhiyi Zhang, emhhbmd6aGl5aTIwMTRAMTYzLmNvbQ==

These authors have contributed equally to this work

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.