Increased Cumulative Incidence of Dermatomyositis in Ulcerative Colitis: a Nationwide Cohort Study.

On a molecular level, two autoimmune diseases: ulcerative colitis (UC) and dermatomyositis share common genetic determinants. On a clinical level, case reports evidenced the co-occurrence of these two diseases. We therefore hypothesize that UC is potentially associated with increased cumulative incidence of dermatomyositis. The goals of this retrospective cohort study were to evaluate whether UC is associated with increased cumulative incidence of dermatomyositis independent of sex and age. For comparison, we also assessed the cumulative incidence of polymyositis in UC and control subjects. The study enrolled 3,133 UC subjects and 14,726 control subjects. The cumulative incidence of dermatomyositis was significantly higher in UC than that of control subjects (p = 0.026), but the cumulative incidence of polymyositis was comparable between UC and control subjects (p = 0.596). UC was independently associated with the increased incident dermatomyositis (hazard ratio: 6.19, 95% confidence interval = 1.77-21.59, p = 0.004) after adjusting for sex, age, and concomitant rheumatoid arthritis, systemic lupus erythematosus, and systemic sclerosis. Similar trends of increased dermatomyositis in UC were observed when patients were stratified based on sex and age. In conclusion, our findings suggest that UC is probably associated with increased cumulative incidence of dermatomyositis, independent of sex, age, and concomitant autoimmune diseases.

Both the incidence and prevalence of ulcerative colitis (UC) is increasing worldwide, affecting approximately 0.5% of the general population in the Western world1. Epidemiological data, genetic evidence and clinical features suggest that UC is a polygenic disease. Despite advances in medical management, morbidity remains high, with 30% of affected individuals undergoing colectomy2 due to failure of medical therapy, among other reasons. Though the precise etiology is not completely understood, the current hypothesis is a dysregulated mucosal immune response to commensal gut flora or microbiota in genetically susceptible individuals. Recent meta-analyses have identified several susceptibility loci associated with UC, including vitamin D receptor (VDR)3 and interferon regulatory factor 5 (IRF5)4. These two susceptibility genes are also associated with dermatomyositis (DM)56, another rarer autoimmune disease compared to UC that affects the musculoskeletal system and respiratory tract. Interestingly, there are several case reports describing the overlapping of these two autoimmune diseases78. The low occurrence of either UC or DM has made accumulating sufficient cases to analyze cumulative incidence and characteristics of each challenging. A high cumulative incidence of cases with both UC and DM co-occurring would suggest genetic overlapping of the two diseases. Studies have shown that, the association of UC with other autoimmune diseases is an important factor in the prognosis of patients diagnosed with UC9. Additionally, early recognition of patients with concurrent DM is associated with higher remission rates, shorter treatment period, and better outcomes for patients diagnosed with DM10. We thus performed a nationwide cohort study to evaluate the cumulative incidence of DM in UC patients and examined the association between these two diseases. Because there are data indicating that DM and polymyositis (PM) share common pathological and immunological characteristics1112, we also utilized this nationwide cohort study to estimate the cumulative incidence of PM in UC patients.

Methods

This nationwide cohort study was established using data from the National Health Insurance Research Database (NHIRD) and data from the Registry of Catastrophic Illness. Because these two databases were released with identification numbers encrypted, the researchers were blinded to the identity of individuals. Thus the analysis of the databases was exempt from the ethical approval by the Institutional Review Board of Kaohsiung Medical University Hospital (KMUHIRB-EXEMPT(I)-20150051).

Source of Sample Subjects

Taiwan initiated its National Health Insurance (NHI) program in 1995 to provide affordable health care for all citizens, with coverage rate more than 99.6%13. The NHI medical claims database, managed by Taiwan’s National Health Research Institutes (NHRI), has collected the data of all ambulatory care and in-patient claims from NHI program enrollees. The NHI database diagnosis coding mirrors the International Classification of Diseases, Ninth Revision (ICD-9), Clinical Modification diagnostic criteria. The NHRI offered NHIRD 1,000,000 random subjects for this study, which comprised registration data on 1,000,000 beneficiaries randomly sampled from the original NHIRD. This was a representative national sample, with comparable distribution in sex, age, or average insured payroll-related amount between the claimants in the sampled data and the original NHIRD14. Every record in the NHIRD contains information about the date of doctor visit, age, sex, and diagnosis of patients. We used a longitudinal cohort of 1,000,000 subjects from January 1, 1998, to December 31, 2011.

Another database offered by NHRI is the Registry of Catastrophic Illness. Every record in the Registry of Catastrophic Illness contains information about the date on which the patient received their first diagnosis of catastrophic illness, the ICD-9 of catastrophic illness, the sex and age of the patient. To avoid financial hardship for families coping with major illnesses, the NHI specifies 31 categories of catastrophic illness (including UC, DM, and PM) for which patients are waived from any copayment. To be registered in the Registry of Catastrophic Illness for UC, the attending physician of these patients, usually a gastroenterologist, must submit related information, including the clinical history, endoscopic, and/or histological findings. Every case must be validated by another expert gastroenterologist to be formally approved and registered. Application of catastrophic illness certificate for DM or PM requires a thorough clinical and laboratory survey that fulfills the criteria proposed by Bohan and Peter1516. After chart review, particularly reviewing clinical symptoms of muscle weakness, dermatological manifestations, electromyography results, muscular biopsy results, and levels of muscle enzymes (e.g. creatine kinase), lactate dehydrogenase, aspartate aminotransferase and alanine aminotransferase, only those with the diagnosis of DM or PM are issued a catastrophic illness certificate for the diseases.

To ensure coding accuracy, the NHI Bureau routinely implements cross-checks of each hospital claim with medical charts. Punitive measures were taken against coding infractions to avoid diagnosis up-coding. However, NHI’s reimbursement system ties a hospital’s reimbursement to its patient severity level. As a result, it is in the best interest of a hospital to accurately code diagnoses and care items. These databases have also been utilized for epidemiologic studies of various autoimmune diseases in numerous articles, including UC, DM, and PM17181920. Therefore, the high standard set by the databases allows researchers to analyze the epidemiological profile using the data.

Cohort Selection

Patients with newly diagnosed UC (ICD-9 code: 556)17 from the Registry of Catastrophic Illness for January 1, 1998, to December 31, 2011 were identified as the UC cohort. We only included patients older than 20 years with no diagnosis of DM (ICD-9: 7103) or PM (ICD-9: 7104) before or within 1 year after their diagnosis of UC. Previous studies171819 also utilized ICD-9 to identify patients with UC, DM, and PM in these databases. We further assigned the date on which they received a diagnosis of UC in the Registry of Catastrophic Illness as the entry date for UC patients. For every UC patient, we randomly selected at most five control subjects without UC, matched in terms of sex, age, and entry date. Control subjects could have other diseases except UC, such as rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), and so on. We recruited control subjects from NHIRD. The entry date for UC patients was the date on which the patients received their diagnosis of UC in the Registry of Catastrophic Illness. Similarly for control subjects, the entry date was a matched date on which control subjects used any medical services. Moreover, we ensured that none of our control subjects had ever received a diagnosis of DM or PM prior to or within one year of their entry date. We used this process to ensure that both cohorts shared the same baseline follow-up time. Patients who were followed less than 1 year were excluded from the analysis due to inadequate duration to assess long-term exposure status21.

Outcome Definition

Any participant with a DM (ICD-9: 7103) or PM (ICD-9: 7104) diagnosis was identified from the Registry of Catastrophic Illness1819. The endpoint of follow-up was: the date of last visit before December 31, 2011, or the diagnosis date of DM or PM in the Registry of Catastrophic Illness, whichever came first.

Assessment of Covariates

Because other autoimmune diseases including: RA, SLE, and systemic sclerosis (SSc)22, might be associated with DM or PM, we also ascertained the presence of RA (ICD-9: 7140), SLE (ICD-9: 7100), and SSc (ICD-9: 7101) from the Registry of Catastrophic Illness. Previous studies1819 also utilized ICD-9 to identify patients with RA, SLE, and SSc in this database.

Statistical Analysis

Since patients with DM or PM diagnosed before or within 1 year after the diagnosis of UC were excluded from the study, the first year after an UC diagnosis was not included in the follow-up time. The follow-up time was therefore defined as starting one year after UC diagnosis for all subjects. The follow-up period was calculated from a year after the cohort entry date to the end of follow-up.

We first evaluated the incidence rate of DM/PM and the corresponding 95% confidence interval (95% CI) in patients with UC and control subjects with person-years as the denominator under the Poisson assumption using Exact method23. We then used the Kaplan-Meier analysis to estimate the cumulative incidence of DM or PM after enrollment and differences in the full time-to-event distributions between different groups were compared by a log-rank test. After confirming the assumption of proportional hazards by Schoenfeld residuals trend tests, which examined the interactions between predictors and event time, we constructed Cox proportional hazard regression models with age, sex, RA, SLE, SSc, and UC adjusted simultaneously in the model. Statistical significance was inferred at a two-sided p < 0.05. All analyses were performed with the SPSS (v19.3) after mining the national medical records using the PERL (v5.8).

Results

Patient Characteristics

This study included a total of 3,133 patients with UC and 14,726 control subjects for the period spanning 1998 through 2011. The mean age was 50.42 ± 16.52 in male UC subjects and 49.70 ± 16.00 in male control subjects, 53.76 ± 15.38 in female UC subjects and 53.02 ± 14.82 in female control subjects (Table 1). The distribution by age was comparable between UC patients and control subjects in both male (p = 0.298) and female (p = 0.244) groups.

Table 1

The age distribution of patients with ulcerative colitis and control subjects in baseline.

	Males		Females
Gender	UC	Control subjects	UC	Control subjects
Number	(n = 1704)	(n = 7955)	(n = 1429)	(n = 6771)
Age(years; mean ± SD)	50.42 ± 16.52	49.70 ± 16.00	53.76 ± 15.38	53.02 ± 14.82
Age group	n (%)	n (%)	n (%)	n (%)
> = 20, <30	202 (11.85)	948 (11.92)	89 (6.23)	419 (6.19)
> = 30, <40	303 (17.78)	1438 (18.08)	196 (13.72)	937 (13.84)
> = 40, <50	351 (20.60)	1732 (21.77)	295 (20.64)	1471 (21.73)
> = 50, <60	320 (18.78)	1527 (19.20)	318 (22.25)	1593 (23.53)
> = 60, <70	265 (15.55)	1260 (15.84)	271 (18.96)	1279 (18.89)
> = 70, <80	205 (12.03)	850 (10.69)	203 (14.21)	884 (13.06)
> = 80, <90	55 (3.23)	194 (2.44)	54 (3.78)	182 (2.69)
> = 90	3 (0.18)	6 (0.08)	3 (0.21)	6 (0.09)
p	0.298		0.244

SD: standard deviation.

DM and PM Incidence Rate in the UC Cohort

A total of 12 cases of incidental DM were identified during the follow-up period, 5 in the UC group and 7 in the control subjects. The incidence rate of DM was 0.36 cases per 1000 person-years (95% CI = 0.12–0.85) in the UC cohort, compared to 0.11 cases per 1000 person-years (95% CI = 0.04–0.22) in the control group (Table 2). UC tended to have higher incidence rate of DM than control subjects. Because age modified risk of extraintestinal manifestations24 and disease severity25 in previous studies, we also stratified patients by age. The DM incidence rate of UC was consistently higher than that of the control subjects in respective age groups.

Table 2

Incidence rate of dermatomyositis and polymyositis in ulcerative colitis patients and control subjects stratified by age.

Age at baseline (years)	UC (n = 3133)			Control subjects (n = 14726)
	No.a	follow-upb	Incidence ratec	95% CI	No.	follow-up	Incidence rate	95% CI
DM
20–69	4	12101.58	0.33	0.09–0.85	6	58502.42	0.10	0.04–0.22
> = 70	1	1699.92	0.59	0.01–3.28	1	7006.42	0.14	0.00–0.79
Total	5	13801.50	0.36	0.12–0.85	7	65508.83	0.11	0.04–0.22
PM
Total	3	13822.83	0.22	0.04–0.63	10	65481.17	0.15	0.07–0.28

aNumbers of new diagnosed DM/PM.

bThe total follow-up period (person-years).

cIncidence rate per 1000 person-years.

A total of 13 cases of incidental PM were identified during the follow-up period, 3 in the UC group and 10 in the control subjects. The incidence rate of PM was 0.22 cases per 1000 person-years (95% CI = 0.04–0.63) in the UC cohort, close to that in the control group (0.15 cases per 1000 person-years; 95% CI = 0.07–0.28) (Table 2).

DM and PM Incidence Rate in the Male Cohort

To clarify the relationship of sex in regards to the incidence rate of DM and PM, we stratified the patients based on sex. The incidence rate of DM was 0.13 cases per 1000 person-years (95% CI = 0.00–0.74) in the male UC cohort, compared to 0.03 cases per 1000 person-years (95% CI = 0.00–0.16) in the male control group (Table 3). The overall incidence rate of DM in the male UC cohort was consistently higher than that in the male control group.

Table 3

Incidence rate of dermatomyositis and polymyositis in men.

Age at baseline (years)	UC (n = 1704)			Control subjects (n = 7955)
	No.a	follow-upb	Incidence ratec	95% CI	No.	follow-up	Incidence rate	95% CI
DM
Total	1	7506.17	0.13	0.00–0.74	1	33785.25	0.03	0.00–0.16
PM
Total	1	7507.83	0.13	0.00–0.74	2	33784.67	0.06	0.01–0.21

aNumbers of new diagnosed DM/PM.

bThe total follow-up period (person-years).

cIncidence rate per 1000 person-years.

For PM, the incidence rate of PM was 0.13 cases per 1000 person-years (95% CI = 0.00–0.74) in the male UC cohort, compared to 0.06 cases per 1000 person-years (95% CI = 0.01–0.21) in the male control group (Table 3).

DM and PM Incidence Rate in the Female Cohort

The incidence rate of DM was 0.64 cases per 1000 person-years (95% CI = 0.17–1.63) in the female UC cohort, compared to 0.19 cases per 1000 person-years (95% CI = 0.07–0.41) in the female control group (Table 4). The incidence rate of DM in the female UC cohort was also consistently higher than that in the control group. For PM, the incidence rate of PM was 0.32 cases per 1000 person-years (95% CI = 0.04–1.14) in the female UC cohort, compared to 0.25 cases per 1000 person-years (95% CI = 0.11–0.50) in the female control group (Table 4). The overall incidence rate of PM in the female UC cohort approached to that in the female control group.

Table 4

Incidence rate of dermatomyositis and polymyositis in women.

Age at baseline (years)	UC (n = 1429)			Control subjects (n = 6771)
	No.a	follow-upb	Incidence ratec	95% CI	No.	follow-up	Incidence rate	95% CI
DM
Total	4	6295.33	0.64	0.17–1.63	6	31723.58	0.19	0.07–0.41
PM
Total	2	6315.00	0.32	0.04–1.14	8	31696.50	0.25	0.11–0.50

aNumbers of new diagnosed DM/PM.

bThe total follow-up period (person-years).

cIncidence rate per 1000 person-years.

Cumulative Incidences of Dermatomyositis and Polymyositis

The cumulative incidence of DM in the UC cohort was statistically higher than that in the control group (p = 0.026; Fig. 1a) while the cumulative incidence of PM in the UC cohort was comparable to that in the control group (p = 0.596; Fig. 1b). We next performed a stratified analysis according to sex. A recurring trend of increased cumulative incidence of DM in UC patients across males and females was observed. In the male cohort, the cumulative incidence of DM in UC was higher than that in the control group, but the difference was not statistically significant (p = 0.245; Supplementary Fig. 1). With regard to PM, the cumulative incidence of PM in UC was comparable to that in the control group (p = 0.484; Supplementary Fig. 1). In the female cohort, the cumulative incidence of DM in UC was statistically higher than that in the control group (p = 0.044; Supplementary Fig. 2), while the cumulative incidence of PM in UC was comparable to that in the control group (p = 0.784; Supplementary Fig. 2).

Figure 1

Cumulative incidence of DM and PM.

(a) The cumulative incidence of DM was higher in UC patients than in control subjects (p = 0.026, estimated by the log-rank test). (b) The cumulative incidence of PM was comparable between UC patients and control subjects (p = 0.596, estimated by the log-rank test).

Hazard Ratio of Ulcerative colitis and Other Covariates in Multivariable Analyses

To determine if UC was an independent risk factor for DM, we used a multivariate Cox proportional hazard regression to calculate the hazard ratio (HR) of UC and other covariates. UC (HR: 6.19, 95% CI = 1.77–21.59, p = 0.004) in addition to RA (HR: 6.30, 95% CI = 1.71–23.22, p = 0.006), SLE (HR: 13.76, 95% CI = 2.96–64.10, p = 0.001), and SSc (HR: 6.45, 95% CI = 1.32–31.59, p = 0.021) was associated with a higher incident DM after adjusting for sex and age, RA, SLE, SSc, and UC (Table 5) concomitantly and was statistically significant.

Table 5

The hazard ratio of various autoimmune diseases for dermatomyositis in multivariable analyses.

Variables	HRa (95% CI)	p
UC	6.19 (1.77–21.59)	0.004
RA	6.30 (1.71–23.22)	0.006
SLE	13.76 (2.96–64.10)	0.001
SSc	6.45 (1.32–31.59)	0.021

aHazard ratio, adjusted for age, sex, UC, RA, SLE, and SSc.

Discussion

In this study, the cumulative incidence of DM appears to be significantly higher in UC patients when compared with the control group, and UC can potentially be a risk factor for the co-occurrence of DM independent of other autoimmune diseases. In contrast, cumulative incidence of PM was comparable between UC patients and the control group. To the best of our knowledge, this is the first nationwide study suggesting a positive correlation of UC with incident DM. The findings of the present study support the speculation that UC could be associated with increased cumulative incidence of DM, and it leads to several interesting points of discussion.

In this study, the presence of an existing diagnosis of UC as an independent predictor of a future diagnosis of DM when compared to the control group, were consistent across sexes (0.13 vs 0.03 in men, 0.64 vs 0.19 in women) and age (0.33 vs 0.10 in young patients, 0.59 vs 0.14 in old patients) (Table 2). In previous studies, there were reports of an association between UC and DM7826. Although the evidence is limited, these observations provide a rationale for the increased cumulative incidence of DM in UC patients as shown.

UC continues to be a major worldwide health problem, accounting for substantial costs to the health care system and society27. Several studies have identified several indicators that confer a poor prognosis associated with UC. Examples include: the presence of antinuclear antibodies in DM patients28 which are associated with uveitis, steroid dependency29, and paradoxical rheumatological manifestations resulting from treatment involving anti- tumour necrosis factor agents30. Additionally, the presence of mucocutaneous manifestations occurring in DM were also associated with primary sclerosing cholangitis in UC patients31. Moreover, the presence of any autoimmune disorders have also been associated with adverse outcomes of ileal pouch-anal anastomosis32, increased pancolitis3334, and overall increased clinical severity935 in UC patients. These observations give an excellent rationale for future studies to clarify whether UC patients with DM have a more complicated clinical course than those without.

Furthermore, patients with UC have increased complications, lower remission rate of other autoimmune diseases3637, increased recurrence rate of other autoimmune diseases38, and higher incidental colorectal cancer rate39 than those without. These findings support the need for further study of DM in patients with UC, its relationship, and its impact.

The mechanisms by which UC increases susceptibility to DM but not PM are yet to be completely elucidated. Past studies observed that DM and PM are characterized with different immunohistological features40, cytokine expressions41, autoantibodies42, clinical course, and prognostic factors43. DM and PM also showed discrepant associations with cancer subtypes in previous epidemiology studies44. Additionally, common genetic factors including IRF5 rs4728142 polymorphism and VDR rs2228570 polymorphism are correlated with both the diagnosis of UC34 and DM56. However. there is no compelling evidence indicating that these single nucleotide polymorphisms increase PM. Thus similar gene polymorphisms in both UC and DM might explain the increased cumulative incidence of DM but not PM occurring in UC patients.

Our study has several strengths. First, we used the NHI dataset, which is a representative national sample, minimizing selection bias resulting from non-response or loss to follow-up of study subjects. In addition, the diagnoses of UC, DM, and PM were confirmed by their inclusion in the Registry of Critical Illness. However, this study also had some limitations. First, administrative databases do not contain information on statin prescription and ultraviolet exposure, which are associated with an increased DM4546. However, previous studies suggest that statin use in UC patients is comparable to that in the general population47, and there is no association between ultraviolet exposure and UC incidence48. Therefore, the lack of information on these variables is unlikely to have introduced substantial bias. Next, this study did not examine effects of a family history of DM. To date, there have been no studies documenting associations between a family history of DM and the cumulative incidence of DM, although it is generally believed that a family history of DM implies a higher cumulative incidence of DM. However, it supported the link between UC and DM rather than refuted it even when confoundings by a family history of DM occurred.

Moreover, our database does not contain laboratory data and clinical information, such as disease phenotypes and serology reports. Thus we didn’t investigate the risk factors for the co-occurrence of UC and DM from this study. This would have been of interest to determine whether specific types of autoantibodies or disease phenotypes could be predictive of UC and DM co-occurrence. Future studies are clearly needed to address this question.

Other limitations include the relatively small number of DM cases, even in a nationwide cohort study, which might make the statistical models unstable. Regardless, this is the first nationwide study to assess the cumulative incidence of DM in UC patients and to suggest an effect of UC on cumulative incidence of DM independent of sex, age, and concomitant autoimmune diseases.

Another concern is UC patients might initiate earlier consultations with healthcare professionals for DM, because of their frequent clinic visits. However, previous studies showed little diagnosis delay in DM49. Furthermore, these scenarios may have affected the true cumulative incidence of DM to a certain degree, but they cannot bias the results considerably because of the short diagnosis delay in DM4950 in addition to easy accessibility and high coverage of universal health insurance in Taiwan13. Lastly, we cannot rule out the possibility of miscoding in an established database that relies on physician-reported diagnoses. However, such miscoding and misclassification are likely to affect both UC and control cohorts equally and even if misclassification did result in underestimation of the association, our results still support the association of UC with DM.

In summary, this is the first preliminary observation of increased cumulative incidence of DM in patients with UC in a nationwide cohort study. Our study raises the possibility that co-occurrence of common as well as rare autoimmune diseases could help elucidate the pathophysiologic background of rare autoimmune diseases that are not addressed extensively due to their low prevalence. In addition, genetic components (IRF5 rs4728142 and VDR rs2228570) common to UC and DM suggest that a common treatment may be effective against both these diseases. Current treatments for UC probably become novel treatment options for DM, for which it has been difficult to develop specific treatment due to its low frequency. For example, tacrolimus has been shown to exert therapeutic effects similar to biologics in UC51. There are also reports demonstrating its efficacy in DM52. Moreover, co-existing UC and DM might represent unique entities different from either process alone. Further research is necessary to confirm our findings, to clarify the clinical course of UC associated with DM, and to determine whether these patients require unique therapeutic approaches to achieve remission of either disease.

Additional Information

How to cite this article: Tseng, C.-C. et al. Increased Cumulative Incidence of Dermatomyositis in Ulcerative Colitis: a Nationwide Cohort Study. Sci. Rep. 6, 28175; doi: 10.1038/srep28175 (2016).