
==== Front
Ther Adv Musculoskelet Dis
Ther Adv Musculoskelet Dis
TAB
sptab
Therapeutic Advances in Musculoskeletal Disease
1759-720X
1759-7218
SAGE Publications Sage UK: London, England

10.1177/1759720X241273039
10.1177_1759720X241273039
Original Research
Risk of motor vehicle accidents with profound injuries in patients with ankylosing spondylitis: a nationwide, population-based cohort study
https://orcid.org/0000-0002-2863-4350
Kao Chung-Mao Division of Allergy, Immunology, and Rheumatology, Department of Internal Medicine, Taichung Veterans General Hospital, Taichung, Taiwan
Division of Translational Medicine, Department of Medical Research, Taichung Veterans General Hospital, Taichung, Taiwan
Doctoral Program in Tissue Engineering and Regenerative Medicine, College of Medicine, National Chung Hsing University, Taichung, Taiwan
Department of Post-Baccalaureate Medicine, College of Medicine, National Chung Hsing University, Taichung, Taiwan
Investigation Writing – original draft Writing – review & editing
Ho Wei-Li Division of Allergy, Immunology and Rheumatology, Chiayi Branch, Taichung Veterans General Hospital, Chiayi, Taiwan
Conceptualization Data curation Funding acquisition Investigation Project administration
Chen Yi-Ming Division of Allergy, Immunology, and Rheumatology, Department of Internal Medicine, Taichung Veterans General Hospital, Taichung, Taiwan
Division of Translational Medicine, Department of Medical Research, Taichung Veterans General Hospital, Taichung, Taiwan
Department of Post-Baccalaureate Medicine, College of Medicine, National Chung Hsing University, Taichung, Taiwan
School of Medicine, National Yang Ming Chiao Tung University, Taipei, Taiwan
Institute of Biomedical Science and Rong Hsing Research Center for Translational Medicine, National Chung Hsing University, Taichung, Taiwan
Data curation Formal analysis Investigation Methodology
Hsieh Tsu-Yi Division of Allergy, Immunology, and Rheumatology, Department of Internal Medicine, Taichung Veterans General Hospital, Taichung, Taiwan
Department of Post-Baccalaureate Medicine, College of Medicine, National Chung Hsing University, Taichung, Taiwan
Department of Medical Education, Taichung Veterans General Hospital, Taichung, Taiwan
Resources Supervision Validation Visualization
Huang Wen-Nan Division of Allergy, Immunology, and Rheumatology, Department of Internal Medicine, Taichung Veterans General Hospital, Taichung, Taiwan
Department of Post-Baccalaureate Medicine, College of Medicine, National Chung Hsing University, Taichung, Taiwan
School of Medicine, National Yang Ming Chiao Tung University, Taipei, Taiwan
Supervision Validation Visualization
Chen Yi-Hsing Division of Allergy, Immunology, and Rheumatology, Department of Internal Medicine, Taichung Veterans General Hospital, Taichung, Taiwan
Department of Post-Baccalaureate Medicine, College of Medicine, National Chung Hsing University, Taichung, Taiwan
School of Medicine, National Yang Ming Chiao Tung University, Taipei, Taiwan
Supervision Validation Visualization
https://orcid.org/0000-0002-7304-4587
Chen Hsin-Hua Division of Allergy, Immunology, and Rheumatology, Department of Internal Medicine, Taichung Veterans General Hospital, No. 1650, Sec. 4, Taiwan Boulevard, Taichung 40705, Taiwan
Department of Post-Baccalaureate Medicine, College of Medicine, National Chung Hsing University, Taichung, Taiwan
School of Medicine, National Yang Ming Chiao Tung University, Taipei, Taiwan
Institute of Biomedical Science and Rong Hsing Research Center for Translational Medicine, National Chung Hsing University, Taichung, Taiwan
Department of Medical Education, Taichung Veterans General Hospital, Taichung, Taiwan
Department of Industrial Engineering and Enterprise Information, Tunghai University, Taichung, Taiwan
Institute of Public Health and Community Medicine Research Center, National Yang Ming Chiao Tung University, Taipei, Taiwan
Big Data Center, National Chung Hsing University, Taichung, Taiwan
Division of General Medicine, Department of Internal Medicine, Taichung Veterans General Hospital, Taichung, Taiwan
Conceptualization Formal analysis Funding acquisition Methodology Project administration Resources Software Supervision Writing – review & editing
shc5555@hotmail.com
20 9 2024
2024
16 1759720X2412730396 2 2024
10 7 2024
© The Author(s), 2024
2024
SAGE Publications Ltd unless otherwise noted. Manuscript content on this site is licensed under Creative Commons Licenses
https://creativecommons.org/licenses/by-nc/4.0/ This article is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License (https://creativecommons.org/licenses/by-nc/4.0/) which permits non-commercial use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/open-access-at-sage).
Background:

Patients with ankylosing spondylitis (AS) suffer from impaired physical activity and are prone to motor vehicle accidents (MVA), but definite instruction regarding the relationship between disease evolvement and MVA and potential risk factors is lacking.

Objectives:

To explore the risk factors and their impact on recorded MVA with profound injuries in AS patients with prescriptions.

Design:

Nationwide, population-based, matched retrospective cohort study.

Methods:

Using Taiwanese administrative healthcare databases, with available claims data from 2003 to 2013, we selected 30,911 newly diagnosed adult AS patients with concurrent prescriptions from 2006 to 2012 as AS patients, along with 309,110 non-AS individuals as the control group, matched in gender, age at index date and year of the index date. The risk of recorded MVA with profound injuries was compared between the two groups in terms of incidence rate ratio (IRR) and log-rank test p-value. Using Cox regression analysis, we studied associations between the risk and AS diagnosis.

Results:

The risk of recorded MVA with profound injuries in AS patients was significantly higher than in non-AS individuals, specifically 2 years after AS diagnosis (IRR, 2.00; 95% confidence interval (CI), 1.42–2.81). For patients with follow-up periods >2 years, the adjusted risk was positively associated with suburban residence (adjusted hazard ratio (aHR), 2.18; 95% CI, 1.55–3.06), rural residence (aHR, 1.89; 95% CI, 1.27–2.80), lower insured income (aHR, 1.35; 95% CI, 1.01–1.81) and recorded MVA with profound injuries before AS diagnosis (aHR, 6.16; 95% CI, 2.53–14.96). AS diagnosis (aHR, 1.81; 95% CI, 1.27–2.59) and frequency of ambulatory visits (aHR, 1.01; 95% CI, 1.004-–1.02) were specific associated factors for them compared with those with follow-up periods ⩽2 years.

Conclusion:

For adult AS patients in Taiwan, factors such as AS disease evolution and frequent ambulatory visits for disease control in the second year of the disease course may significantly increase the risk of recorded MVA with profound injuries beyond 2 years after AS diagnosis.

ankylosing spondylitis
claims data
cohort study
health-related quality of life
motor vehicle accident
population-based study
Taichung Veterans General Hospital https://doi.org/10.13039/501100010101 TCVGH-VHCY1088605 Taichung Veterans General Hospital https://doi.org/10.13039/501100010101 TCVGH-VHCY1138604 cover-dateJanuary-December 2024
typesetterts1
==== Body
pmcIntroduction

Ankylosing spondylitis (AS), currently termed radiographic axial spondyloarthritis (axSpA), is a systemic inflammatory disease characterised by spondyloarthritis, sacroiliitis and enthesitis. Many patients with AS suffer from peripheral arthritis, hip or shoulder involvement,1,2 ankylosis and extra-articular manifestations like acute uveitis, which lead to impaired physical activity and a considerable healthcare burden. 3 With systemic involvement of AS, they also face multiple comorbidities, such as fatigue, 4 fibromyalgia, 5 sleep disturbance 6 and psychiatric disorders. 7 In Taiwan, AS featured a standardised prevalence rate of 0.24% in 2010 and an incidence rate of 0.42–0.50 cases per 1000 person-years. 8 AS was diagnosed by Taiwanese rheumatologists based on 1984 modified New York criteria 9 or 2009 Assessment of Spondyloarthritis International Society classification criteria for axSpA,10,11 and current management in Taiwan is based on the updated regional consensus recommendations. 12 With growing global emphasis on optimising quality standards in AS patients, 13 rheumatologists in Taiwan are aiming, with a holistic approach, to identify modifiable factors related to health-related quality of life.14,15

High rates and high fatality in motor vehicle accidents (MVA) are significant public issues in Taiwan.16,17 AS patients may be predisposed to a higher risk of MVA as the disease evolves. 18 This issue deserves more attention from rheumatologists making individualised long-term healthcare plans. Little research is done to elucidate the relationship between the diagnosis of AS and the risk of MVA, leaving an unmet need for definite instructions for AS patients. Without clear clarification, limited attention can be aroused to preventive intervention for injuries related to their functional impairments. In Taiwan, this issue can be addressed with the aid of population-based administrative health database research, which has facilitated many epidemiological studies. 19 In this study, we investigated the associated risk factors of recorded MVA with profound injuries and their impact on the risk in AS patients with ambulatory visits for pharmacological control and followed for up to 7 years.

Materials and methods

Study design and source of data

We conducted a nationwide, population-based, matched retrospective cohort study. Claims data were acquired from the Taiwan National Health Insurance Research Database (NHIRD) and Longitudinal Health Insurance Database 2010 (LHID 2010). The National Health Insurance (NHI) Programme is an obligatory, single-payer copayment scheme that has been implemented since March 1, 1995. It is currently the only public health insurance programme and has covered >99% of the entire Taiwanese population. NHIRD is released from a database of the system for research purposes, and it can be interlinked with the National Death Registry. Information collected includes claims data regarding registration, demographics, dates of outpatient and admission services, diagnoses and procedures with their International Classification of Diseases-Ninth Revision-Clinical Modification (ICD-9-CM) codes, medication prescriptions with their Anatomical Therapeutic Chemical Classification codes and medical expenditures covered by NHI. However, it does not include information like weight and height, personal history, alcohol and tobacco use, substances use without NHI reimbursement, laboratory measures and AS activity, damage and functional measures. As a subset of the health dataset originating from NHIRD, the LHID 2010 database randomly sampled 1 million beneficiaries from NHIRD participants in 2010. It contains the most updated claims data of sampled individuals since 1997. It is validated by Taiwanese researchers to have representative power of the nationwide population. 19 The available claims data in both databases started from January 1, 2003 to December 31, 2013.

Study subjects and identification of AS patients and matched non-AS individuals

In our study, newly diagnosed AS patients aged ⩾20 years with pharmacological control were identified via NHIRD to include those having ⩾3 rheumatology ambulatory visits with ICD-9-CM code of AS and a concurrent prescription of non-steroidal anti-inflammatory drugs, methotrexate or sulfasalazine for disease control from January 1, 2006 to December 31, 2012 (Supplemental Table S1). When AS was identified as ⩾2 visits with the corresponding ICD-9-CM code given by a rheumatologist, the reported positive predictive value was 81% and sensitivity was 67% from experience in the United States. 20 The selected patients were called the AS patient group. The other subset of patients with at least one ambulatory visit from January 1, 2006 to December 31, 2012 was identified via LHID 2010. After excluding those with a diagnosis of AS, the remaining patients were defined as non-AS individual groups that acted as the control group. In both groups, subjects with incomplete records or loss of follow-up by 2012 were excluded. Patients in both groups were each assigned an index date, which was the date of their first AS diagnosis for AS patients, and the date of first ambulatory visit for non-AS individuals. After matching for the year of index dates and pairing individuals across the two groups, we then determined if the occurrence of potential risk factors within a year before AS diagnosis had influenced the risk after AS diagnosis. We matched the AS group with the non-AS group at a ratio of 1:10 for gender, age at the index date and year of the index date. After matching, selected study subjects in both groups were included for final analysis. The process for the selection of study subjects for ultimate analysis is shown in Figure 1.

Figure 1. Flow diagram of enrolment, categorisation and matching for comparison of the study population.

AS, ankylosing spondylitis; LHID, Longitudinal Health Insurance Database; NHIRD, National Health Insurance Research Database; OPD, outpatient department.

Study outcomes

The primary investigated outcome was MVA with profound injuries that necessitated emergent medical exams or management, such as fractures or acquired spondylolisthesis, recorded as the chief diagnosis at emergency department visits, which was identified by ICD-9-CM codes (E810-E819) and concomitant injury codes (800-999) (Supplemental Table S1). The follow-up period of each patient was the interval between the index date and the first of the following censoring events: occurrence of recorded MVA with profound injuries, December 31, 2013, date of withdrawal from the NHI programme for any reason or death.

For subjects experiencing recorded MVA with profound injuries at the emergency department, the associated complications (major fractures and acquired spondylolisthesis) during and from episodes of MVA, and short-term adverse outcomes (surgeries and death) within 1 month after MVA were demonstrated, which were identified by at least one ambulatory or inpatient visit with corresponding ICD-9-CM codes and interlinked data from National Death Registry (Supplemental Tables S1 and S2).

Potential confounders for MVAs

We listed several covariates to be adjusted, including the AS diagnosis, gender, age at the index date, socioeconomic status by urbanisation level and insured monthly income, and MVA with profound injuries recorded at the emergency department before the index date. Some covariates that had occurred within 1 year before the index date were also adjusted, including history of spinal or joint surgeries (Supplemental Table S2), frequency of recorded outpatient department visits, comorbidities, extra-articular manifestations and use of medications including systemic corticosteroids and those predisposing patients to trance. Age at index date was classified into four groups by quartile distribution to observe the trend of the risk concerning age of diagnosis. The Taiwanese stratification model of urbanisation level was assessed based on population density, physician availability, the proportion of elders, agricultural workers and residents with college-level education in each subject’s residential district. 21 Selected comorbidities for adjustment included myocardial infarction, cerebrovascular diseases, chronic obstructive pulmonary disease, peptic ulcer disease, mild liver disease, diabetes mellitus, anxiety disorders, depressive disorders, bipolar disorders, psychotic disorders and alcohol-related disorders. The extra-articular manifestations, including uveitis, psoriasis and inflammatory bowel disease, were diagnosed specifically by corresponding specialists. The presence of comorbidities and extra-articular manifestations was identified by ⩾3 ambulatory visits or 1 inpatient visit with corresponding ICD-9-CM codes within a year before the index date. The use of medications within a year before the index date was identified through Anatomical Therapeutic Chemical classification codes. These medications included hypnotics (benzodiazepine-related drugs and benzodiazepines), antihistamines, antidepressants, antipsychotics and systemic corticosteroids. The above diseases, events, manifestations, medications and their corresponding codes are summarised in Supplemental Table S1.

Methods to enhance diagnostic accuracy and sensitivity analyses

We included patients with ⩾3 ambulatory visits with a diagnostic code of AS and a concurrent prescription given merely from rheumatologists. The inclusion of covariates within 1 year before the index date was set as only ⩾3 ambulatory visits or one inpatient visit with administrative codes. The included extra-articular manifestations were concurrently recorded by corresponding specialists (i.e. uveitis also recorded by ophthalmologists). These measures were aimed at enhancing diagnostic accuracy.

We compared baseline characteristics in matched study subjects having >2 years of follow-up (Supplemental Table S3). We also compared incidence rate ratio (IRR) with log-rank test p-value, and conducted Cox regressions using the time horizon of within and beyond 2 years after index date for sensitivity analysis.

Statistical analyses

Continuous variables were represented as mean ± standard deviation and categorical variables as the number (percentage) of patients. We compared continuous variables between the two groups using an independent t-test, and categorical variables using the Chi-square test or Fisher’s exact test. Incidence rate in both groups, IRR and results of log-rank tests were generated. Kaplan–Meier curves of cumulative incidence in matched study subjects were compared. The association between AS diagnosis and the risk as determined with Cox proportional hazard regressions to adjust for covariates was presented as an adjusted hazard ratio (aHR) with 95% confidence intervals (CIs). A probability (p) value of <0.05 was considered statistically significant. Statistical calculations were performed using Statistical Package for the Social Sciences, Windows Version 13.0 (SPSS Inc., Chicago, IL, USA).

Reporting standards

The reporting of this study conforms to the Strengthening the Reporting of Observational Studies in Epidemiology statement. 22

Results

Baseline characteristics of the study population

We included 30,911 patients with newly diagnosed adult AS with prescriptions from NHIRD in the AS group and 309,110 matched individuals from LHID 2010 in the non-AS group, among which no subject was excluded due to missing data or loss of follow-up (Figure 1), and baseline characteristics of study subjects of the two groups are shown in Table 1. The study population consisted predominantly of male subjects (collectively 62.9% in the two groups). The mean age of AS patients at diagnosis for the first time was 42.1 ± 16.8 years, represented as the age at index date. We observed no preponderance in the age of AS diagnosis. In comparison, a higher proportion of AS patients resided in urban districts at diagnosis, received more than the median monthly income, had undergone spinal or joint surgeries within 1 year before AS diagnosis (p < 0.001) and had experienced previous emergency department visits due to MVA with injuries (p < 0.001). AS patients also had significantly more frequent outpatient department visits before AS diagnosis (p < 0.001). Statistically greater proportions of AS patients were diagnosed with comorbidities and all extra-articular manifestations before AS diagnosis, with the exceptions of cerebrovascular diseases and psychotic disorders. Regarding medications, greater proportions of AS patients had been prescribed all specified medications before their AS diagnosis (all p < 0.001).

Table 1. Baseline characteristics in matched study subjects.

Baseline characteristic	Non-AS individuals (n = 309,110)	AS patients (n = 30,911)	p-Value	
Gender			1.000	
 Female, n (%)	114,560 (37.1)	11,456 (37.1)		
 Male, n (%)	194,550 (62.9)	19,455 (62.9)		
Age at index date a (years), mean ± SD	42.1 ± 16.8	42.1 ± 16.8	1.000	
Groups of age at index date (years)			1.000	
 Age < 30	83,740 (27.1)	8374 (27.1)		
 30 ⩽ age < 40	65,190 (21.1)	6519 (21.1)		
 40 ⩽ age < 55	88,150 (28.5)	8815 (28.5)		
 55 ⩽ age	72,030 (23.3)	7203 (23.3)		
Urbanisation, at the index date			0.007	
 Urban, n (%)	93,939 (30.4)	9645 (31.2)		
 Suburban, n (%)	150,495 (48.7)	14,793 (47.9)		
 Rural, n (%)	64,676 (20.9)	6473 (20.9)		
Lower insured income, b at index date, n (%)	160,918 (52.1)	13,363 (43.2)	<0.001	
Recorded spinal or joint surgeries, within 1 year before index date, n (%)	3334 (1.1)	772 (2.5)	<0.001	
MVA with profound injuries recorded before the index date, n (%)	928 (0.3)	133 (0.4)	<0.001	
Frequency of OPD visits, within 1 year before the index date (times), mean ± SD	11.8 ± 14.6	23.1 ± 20.4	<0.001	
Comorbidities, within 1 year before the index date	
 Myocardial infarction, n (%)	519 (0.2)	77 (0.2)	0.001	
 Cerebrovascular diseases, n (%)	5262 (1.7)	559 (1.8)	0.170	
 Chronic obstructive pulmonary disease, n (%)	6904 (2.2)	1178 (3.8)	<0.001	
 Peptic ulcer disease, n (%)	8996 (2.9)	2288 (7.4)	<0.001	
 Mild liver disease, n (%)	5183 (1.7)	817 (2.6)	<0.001	
 Diabetes mellitus, n (%)	13,379 (4.3)	1615 (5.2)	<0.001	
 Anxiety disorders, n (%)	7481 (2.4)	1629 (5.3)	<0.001	
 Depressive disorders, n (%)	3807 (1.2)	822 (2.7)	<0.001	
 Bipolar disorders or psychotic disorders, n (%)	2314 (0.7)	255 (0.8)	0.139	
 Bipolar disorders, n (%)	629 (0.2)	122 (0.4)	<0.001	
 Psychotic disorders, n (%)	1783 (0.6)	150 (0.5)	0.041	
 Alcohol-related disorders, n (%)	775 (0.3)	107 (0.3)	0.002	
Extra-articular manifestations, within 1 year before the index date	
 Uveitis, n (%)	134 (0.04)	968 (3.1)	<0.001	
 Psoriasis, n (%)	295 (0.1)	195 (0.6)	<0.001	
 Inflammatory bowel disease, n (%)	21 (0.01)	11 (0.04)	<0.001	
Medications, within 1 year before the index date	
 Hypnotics, n (%)	45,583 (14.7)	9029 (29.2)	<0.001	
  Benzodiazepine-related drugs, n (%)	14,845 (4.8)	2995 (9.7)	<0.001	
  Benzodiazepines, n (%)	39,354 (12.7)	8076 (26.1)	<0.001	
 Antihistamines, n (%)	168,534 (54.5)	22,260 (72.0)	<0.001	
 Antidepressants, n (%)	12,940 (4.2)	3129 (10.1)	<0.001	
 Antipsychotics, n (%)	18,398 (6.0)	3301 (10.7)	<0.001	
 Corticosteroids, n (%)	42,586 (13.8)	9458 (30.6)	<0.001	
a Age at first AS diagnosis for AS patients, and date of first ambulatory visit for non-AS individuals.

b Insured monthly income lower than median monthly income (⩽21,000 New Taiwan dollars/month).

AS, ankylosing spondylitis; MVA, motor vehicle accident; OPD, outpatient department; SD, standard deviation.

For subjects with longer than 2 years of follow-up, we found no predominant distribution of residence between groups, a greater proportion of AS patients with a history of cerebrovascular diseases and no difference regarding history of psychotic disorders and alcohol-related disorders. Otherwise remained consistent (Supplemental Table S3).

Incidence rates of MVA with profound injuries in AS patients and non-AS individuals

The incidence rate was 58 per 100,000 person-years for AS patients, compared with the rate of 38 per 100,000 person-years for non-AS individuals, and AS patients recorded significantly more frequent MVA with profound injuries than non-AS individuals (IRR, 1.53; 95% CI, 1.18–1.98; log-rank test p-value: 0.001) (Table 2). Despite no statistical difference in incidence rate between the two groups based on the period within 2 years after the index date (IRR, 1.13; 95% CI, 0.76–1.69; log-rank test p-value, 0.545), the difference became significant when the period of comparison extended beyond 2 years after the index date (IRR, 2.00; 95% CI, 1.42–2.81; log-rank test p-value: <0.001). Kaplan–Meier curves of cumulative incidence between the two groups showed consistent results (Figure 2).

Table 2. IRR of motor vehicle accidents with profound injuries between matched study subjects.

Variable	Total, n	Event, n (%)	Total follow-up person-time (year)	Average/median follow-up time (year)	Incidence rate (105 per year)	IRR (95% CI)	p-Value	Log-rank test p-value	Proportional hazards assumption	
Entire follow-up period	
 Non-AS individuals	309,110	467 (0.15)	1,238,853	4/4	38	1 (Reference)	0.001	0.001	0.019	
 AS patients	30,911	66 (0.21)	114,253	4/4	58	1.53 (1.18–1.98)				
Within 2 years after the index date	
 Non-AS individuals	309,110	248 (0.08)	566,252	2/2	44	1 (Reference)	0.535	0.545	0.007	
 AS patients	30,911	27 (0.09)	54,371	2/2	50	1.13 (0.76–1.69)				
Beyond 2 years after the index date (maximally 5 years)	
 Non-AS individuals	230,673	219 (0.09)	672,600	3/3	33	1 (Reference)	<0.001	<0.001	0.853	
 AS patients	22,969	39 (0.17)	59,882	3/3	65	2.00 (1.42–2.81)				
AS, ankylosing spondylitis; CI, confidence interval; IRR, incidence rate ratio.

Figure 2. Kaplan–Meier curves of cumulative incidence of motor vehicle accidents with profound injuries among matched study subjects. (a) During the entire follow-up period. (b) Period within 2 years after the index date. (c) Period beyond 2 years after the index date.

Outcomes of study subjects experiencing MVA with profound injuries

Among subjects experiencing recorded MVA with profound injuries at the emergency department, 1.5% of them underwent spinal or joint surgeries within 1 year before MVA (Table 3). During their recorded episodes, no difference was found in the proportions of subjects complicated with major fractures and acquired spondylolisthesis, nor were proportions of those facing surgeries and death within 1 month after MVA as the short-term adverse outcomes. Regarding longer-term influences, no difference in the proportions of subjects with recorded major fractures and acquired spondylolisthesis from their recorded episodes. The results were unchanged in the periods within and beyond 2 years after the index date.

Table 3. Outcomes of study subjects experiencing recorded MVAs with profound injuries at the emergency department.

Event	Entire follow-up period	Within 2 years after the index date	Beyond 2 years after the index date	
	Non-AS individuals (n = 467)	AS patients (n = 66)	p-Value	Non-AS individuals (n = 248)	AS patients (n = 27)	p-Value	Non-AS individuals (n = 219)	AS patients (n = 39)	p-Value	
Within 1 year before the MVA with profound injuries	
 Spinal or joint surgeries, n (%)	7 (1.5)	1 (1.5)	0.992	1 (0.4)	0 (0.0)	0.741	6 (2.7)	1 (2.6)	0.950	
MVA with profound injuries at the emergency department	
 Fractures, n (%) (ICD-9-CM code: 800-829)	66 (14.1)	7 (10.6)	0.435	34 (13.7)	5 (18.5)	0.496	32 (14.6)	2 (5.1)	0.107	
 Fracture of cervical vertebra, n (%) (ICD-9-CM code: 805.0, 805.1)	3 (0.6)	1 (1.5)	0.442	2 (0.8)	1 (3.7)	0.169	1 (0.5)	0 (0.0)	0.672	
 Acquired spondylolisthesis, n (%) (ICD-9-CM code: 738.4)	0 (0.0)	0 (0.0)	Not applicable	0 (0.0)	0 (0.0)	Not applicable	0 (0.0)	0 (0.0)	Not applicable	
Within 1 month after the MVA with profound injuries	
 Spinal or joint surgeries, n (%)	37 (7.9)	4 (6.1)	0.595	21 (8.5)	4 (14.8)	0.276	16 (7.3)	0 (0.0)	0.081	
 All surgeries, a n (%)	74 (15.8)	10 (15.2)	0.885	41 (16.5)	5 (18.5)	0.793	33 (15.1)	5 (12.8)	0.715	
 Death, b n (%)	4 (0.9)	0 (0.0)	0.450	3 (1.2)	0 (0.0)	0.566	1 (0.5)	0 (0.0)	0.672	
From MVA with profound injuries to December 31, 2012	
 Fractures, n (%)	116 (24.8)	15 (22.7)	0.709	65 (26.2)	10 (37)	0.230	51 (23.3)	5 (12.8)	0.144	
 Fracture of cervical vertebra, n (%)	3 (0.6)	1 (1.5)	0.442	2 (0.8)	1 (3.7)	0.169	1 (0.5)	0 (0.0)	0.672	
 Acquired spondylolisthesis, n (%)	3 (0.6)	2 (3.0)	0.060	1 (0.4)	0 (0.0)	0.741	2 (0.9)	2 (5.1)	0.050	
A p-value < 0.05 is considered statistically significant.

a All recorded surgeries or procedures with obtainable ICD-9-CM codes.

b Information on the death of study subjects was acquired from the Taiwan National Death Registry.

AS, ankylosing spondylitis; ICD-9-CM, International Classification of Diseases-Ninth Revision-Clinical Modification; MVA, motor vehicle accident.

Independent risk factors associated with MVA with profound injuries in AS patients

Results of Cox proportional hazard regression analyses during the entire follow-up period showed that the risk would be associated with AS diagnosis (aHR, 1.42; 95% CI, 1.08–1.87), male gender (aHR, 1.35; 95% CI, 1.11–1.64), suburban residence (aHR, 1.82; 95% CI, 1.45–2.28), rural residence (aHR, 1.84; 95% CI, 1.42–2.39), lower insured income (aHR, 1.44; 95% CI, 1.19–1.74) and MVA with profound injuries recorded at emergency department before index date (aHR, 9.40; 95% CI, 5.70–15.48) (Supplemental Table S4). Based on the period within 2 years after index date, it would be associated with male gender (aHR, 1.46; 95% CI, 1.11–1.94), suburban residence (aHR, 1.73; 95% CI, 1.27–2.36), rural residence (aHR, 1.80; 95% CI, 1.26–2.58), lower insured income (aHR, 1.44; 95% CI, 1.10–1.87) and recorded MVA with profound injuries before index date (aHR, 11.57; 95% CI, 6.31–21.22) (Supplemental Table S5). When analysed based on the period beyond 2 years after the index date, it would be associated with AS diagnosis (aHR, 1.81; 95% CI, 1.27–2.59), suburban residence (aHR, 2.18; 95% CI, 1.55–3.06), rural residence (aHR, 1.89; 95% CI, 1.27–2.80), lower insured income (aHR, 1.35; 95% CI, 1.01–1.81), recorded MVA with profound injuries before the index date (aHR, 6.16; 95% CI, 2.53–14.96) and frequency of ambulatory visits within the second year after the index date (aHR, 1.01; 95% CI, 1.004–1.02) (Table 4). The risk would be negatively associated with age at the index date for those initially diagnosed before the age of 30 as reference, which remained consistent within both time horizons of follow-up periods.

Table 4. Risk factors of MVA with profound injuries in AS patients beyond 2 years after the index date.

Variable	Univariable analysis	Multivariable analysis	
	HR (95% CI)	p-Value	Adjusted HR (95% CI)	p-Value	
Diagnosis of AS	1.99 (1.42–2.81)	<0.001	1.81 (1.27–2.59)	0.001	
Male gender	1.28 (0.99–1.67)	0.065	1.26 (0.95–1.66)	0.106	
Age at 2 years after the index date a	0.99 (0.98–0.99)	0.007			
Groups of age at 2 years after the index date	
 Age < 30	1 (Reference)		1 (Reference)		
 30 ⩽ age < 40	0.26 (0.17–0.41)	<0.001	0.28 (0.18–0.44)	<0.001	
 40 ⩽ age < 55	0.43 (0.31–0.59)	<0.001	0.44 (0.31–0.62)	<0.001	
 55 ⩽ age	0.62 (0.46–0.84)	0.002	0.49 (0.34–0.71)	<0.001	
Urbanisation, at 2 years after the index date	
 Urban	1 (Reference)		1 (Reference)		
 Suburban	2.24 (1.59–3.14)	<0.001	2.18 (1.55–3.06)	<0.001	
 Rural	2.09 (1.42–3.08)	<0.001	1.89 (1.27–2.80)	0.002	
Lower insured income, b at 2 years after the index date	1.72 (1.30–2.27)	<0.001	1.35 (1.01–1.81)	0.046	
Recorded spinal or joint surgeries, within the second year after the index date	0.87 (0.28–2.73)	0.817	0.58 (0.19–1.83)	0.355	
MVA with profound injuries recorded before the index date	7.44 (3.07–18.03)	<0.001	6.16 (2.53–14.96)	<0.001	
Frequency of OPD visits, within the second year after the index date	1.01 (1.01–1.02)	<0.001	1.01 (1.004–1.02)	0.004	
Comorbidities, from 1 year before the index date to 2 years after the index date (total 3 years)	
 Myocardial infarction	1.25 (0.18–8.91)	0.823	0.96 (0.13–6.90)	0.963	
 Cerebrovascular diseases	1.08 (0.54–2.19)	0.828	0.76 (0.37–1.60)	0.473	
 Chronic obstructive pulmonary disease	1.13 (0.68–1.87)	0.641	0.86 (0.50–1.46)	0.569	
 Peptic ulcer disease	1.32 (0.90–1.95)	0.155	1.01 (0.65–1.54)	0.982	
 Mild liver disease	0.87 (0.47–1.65)	0.677	0.72 (0.38–1.39)	0.330	
 Diabetes mellitus	1.50 (0.99–2.28)	0.059	1.39 (0.87–2.20)	0.166	
 Anxiety disorders	1.32 (0.83–2.11)	0.244	0.97 (0.56–1.67)	0.910	
 Depressive disorders	1.70 (0.93–3.12)	0.084	1.38 (0.66–2.86)	0.389	
 Bipolar disorders or psychotic disorders	0.35 (0.05–2.46)	0.288	0.21 (0.03–1.56)	0.127	
 Alcohol-related disorders	1.98 (0.63–6.17)	0.241	1.90 (0.59–6.16)	0.285	
Extra-articular manifestations, from 1 year before the index date to 2 years after the index date (total 3 years)	
 Uveitis	0.64 (0.09–4.52)	0.651	0.41 (0.06–3.00)	0.383	
Medications, within the second year after the index date	
 Hypnotics	1.40 (1.05–1.88)	0.024			
  Benzodiazepine-related drugs	1.45 (0.93–2.27)	0.102	1.12 (0.66–1.90)	0.672	
  Benzodiazepines	1.42 (1.04–1.92)	0.026	1.22 (0.84–1.77)	0.301	
 Antihistamines	1.23 (0.95–1.59)	0.120	1.00 (0.76–1.32)	0.990	
 Antidepressants	1.37 (0.85–2.20)	0.202	0.85 (0.46–1.57)	0.596	
 Antipsychotics	1.49 (0.99–2.23)	0.056	1.26 (0.80–1.97)	0.318	
 Corticosteroids	1.25 (0.92–1.70)	0.150	1.02 (0.73–1.41)	0.930	
A p-value < 0.05 is considered statistically significant.

a Age at first AS diagnosis for AS patients, and date of first ambulatory visit for non-AS individuals.

b Insured monthly income lower than median monthly income (⩽21,900 New Taiwan dollars/month).

AS, ankylosing spondylitis; CI, confidence interval; HR, hazard ratio; MVA, motor vehicle accident; OPD, outpatient department.

Discussion

For adult AS patients with pharmacological control in Taiwan, we found a higher risk of recorded MVA with profound injuries 2 years after diagnosis of AS compared with non-AS individuals. Factors including suburban and rural residence, lower insured income and recorded MVA with profound injuries at the emergency department before AS diagnosis might be associated with the risk in the following 7 years of AS course. Nonetheless, for those with follow-up periods >2 years, AS evolution and frequent ambulatory visits for disease control in the second year of the disease course were independent factors and considered to significantly increase the risk beyond 2 years after AS diagnosis. Diagnosis of AS at older ages, specifically >30 years, might be associated with a lower risk.

We observed a higher risk of recorded MVA with profound injuries with a temporal dependency, specifically in the period beyond 2 years after AS diagnosis, suggested by the significantly higher IRR for AS patients in this period and that AS diagnosis was an independent risk factor for those with follow-up periods >2 years. How AS disease progression might have affected the risk was rarely studied in the available epidemiological literature. Functional impairment in AS is largely associated with structural progression 23 and is correlated with radiographic progression in the spine or sacroiliac joints. 24 Prior studies revealed that the progression of radiographic sacroiliitis could be found as early as <2 years in patients with early axSpA25,26 and that structural changes in sacroiliac joints could be detected radiographically in an estimated 5%–10% of patients with non-radiographic axSpA within the first 2 years, increasing to 5%–40% within 10 years after diagnosis. 27 Based on the correlation between structural deterioration and mobility impairment, 28 our results did not contradict the abovementioned findings because we also presented the reduction in physical function connected with the progression of AS. Restriction in driving in AS patients was reported to be attributed to cervical spine involvement. 29 Concerning cervical spine involvement along the course of disease evolution, a Moroccan observational study reported that 19.6% of patients had radiological involvement after 5 years, and such involvement kept increasing in percentage thereafter, which indicated that cervical spine involvement causes physical disability mainly in later disease stages. For Taiwanese AS patients, a higher risk of developing obstructive sleep apnoea was reported, 30 and this was likely correlated with MVA. Of note, as recorded MVAs with profound injuries at the emergency department were investigated only, people experiencing MVAs that did not combine with injuries or those who did not seek medical attention could not be estimated from this study.

We also demonstrated the delayed diagnosis of axSpA as a major clinical challenge in Taiwan, with a mean diagnostic age of 42.1 years, and over half of the study subjects were diagnosed over the age of 40. Delayed diagnosis of AS is known to have worse functional outcomes. 31 To date, there are no available data on the definite duration of diagnostic delay in Taiwan, nor studies reporting clinical presentations in patients with diagnostic delay as the disease propagates. Such potential bias needs to be considered in the interpretation of our results.

Independent risk factors associated with MVA with profound injuries in both periods within and beyond 2 years after the index date included suburban and rural residence, lower insured income and recorded MVA with profound injuries before diagnosis. Factors in the period beyond 2 years after the index date additionally included AS diagnosis and frequent ambulatory visits for AS in the second year of the disease course, which was considered to prominently raise the risk thereafter. Morton et al. 32 found that diagnosis of axSpA is a factor associated with driving difficulty and that patients with uveitis showed greater difficulty with dynamic driving scenarios from their univariate regression analysis. Though our cohort had more study subjects, the development of uveitis was not found to be associated with the risk in our multivariable Cox regression models. This discrepancy is likely related to differences in study design and regression models, as well as reasons causing uveitis not recorded by ophthalmologists. Male patients are known to have more severe radiological progression, specifically in their early disease stages,33,34 and our results showed no contradiction. Patients residing in nonmetropolitan districts and receiving lower monthly income were shown with a higher risk. It coincided with Taiwanese real-world presentations 16 and these patients might tend to be older, in physical-intensive roles, and face presenteeism and reduced work productivity. 35 History of MVA with profound injuries before AS diagnosis was recognised as the most dominant risk factor throughout the follow-up period. Because prior mechanical stress or injury could trigger immune activation and cause the onset and progression of AS in susceptible patients, 36 disabling sequelae and disease progression might incur following MVA. In Taiwanese real-world practice, frequent ambulatory visits likely reflect disease conditions, inflammatory severity or complications necessitating close monitoring. The association between frequent ambulatory visits within the second year of disease follow-up and the risk in patients with follow-up periods >2 years suggested a persistently higher inflammatory activity or complications associated with the higher risk. Lastly, we found those with AS initially diagnosed at ages >30 years to be negatively associated with MVA with profound injuries. However, one study suggested both age-related and AS-specific changes might make elder AS patients more at risk. 37 This finding did not coincide with what prior studies demonstrated,16,17 which might partially reflect that elder AS patients had suffered from severe physical impairments that restricted them from motor vehicle use. Also, juvenile-onset patients are more likely to manifest with peripheral involvement. 38 Therefore, early administration of effective treatment is necessary. 39

We revealed no significant influence of selected psychiatric comorbidities and psychoactive agents on the risk though they were considered risk factors for MVA.40–42 This might reflect the prominent underdiagnosis of included psychiatric comorbidities in Taiwan. The impact of dosages of selected psychoactive agents and subjects’ adherence to medications were not sufficiently explored in this study. These might potentially bias our results.

MVA is a great concern for patients as it can cause permanent damage or disabilities. We revealed the 2 years of disease evolution before a significantly increased risk of MVA with profound injuries, which is shorter than the pooled mean diagnostic delay of 6.7 years. 43 Our results indicated that strategies aiming for timely recognition, effective treatment and comprehensive monitoring of early axSpA might be mandatory. For rheumatologists, definitive instructions on the risk of MVA accompanying disease evolution and corresponding prevention deserve a thorough education for those patients with the risk factors. Our results can provide a reference for Taiwanese administrative authorities to support the implementation of prompt policies.

This study has several limitations. First, some limitations regarding database research exist. Diagnostic accuracy of AS, covariates and outcomes in administrative databases could not be fully confirmed though validated in previous studies, and thus the misclassification bias could not be prevented, either. Potential confounders such as occupation, education level, smoking status, 44 alcohol use, 45 disabilities without obtainable diagnostic codes, alternative substances use, dosages of prescriptions, patients’ compliance to prescriptions, stressful events, 46 family history, laboratory recordings and disease activity could not be adjusted. Neither could issues regarding regional terrain and climate, process of patient involvement in MVA, types of involved vehicles and profiles of regional traffic networks and habits. However, we had enhanced diagnostic accuracy by some measures mentioned in our methods. We matched for gender, age at index date and year of index date to adjust for age of initial AS diagnosis and AS disease duration. We also explored the association between the frequency of ambulatory visits and MVA as the proxy of the association between persistently higher inflammatory activity and MVA. We had adjusted for potential confounders, including psychiatric disorders and alcoholism, as possible under accessible resources of administrative database research. Second, the causal relationship could not be confirmed due to retrospective use of claims data, so we stated association and associated factors in this study. Third, the time horizon of within 1 year before the index date was set to detect continuous exposure to potential confounders and the range of available claims data from 2003 to 2013, determination of the influence of potential confounders occurring longer before the index date was limited. Fourth, because the study population does not represent the entire AS population, but rather represents relatively more severe patients, we cannot extrapolate the results to AS patients with milder severity or without prescriptions. Fifth, the generalisability of our results to regions other than Taiwan was limited owing to regional geographical and administrative differences.

Despite these limitations, population-based administrative databases provided reliable data and minimised selection bias. Also, we conducted sensitivity analyses using different time horizons of disease courses for risk stratification. Following studies bearing prospective designs, longer follow-up periods or methods allowing definitive recording of real exposed dosage and duration of medications, use of tobacco and alcohol, comorbidities, and information on disease activity are warranted to confirm our findings.

Conclusion

This nationwide, population-based cohort study on adult Taiwanese AS patients with pharmacological control showed that factors such as AS disease evolution and frequent ambulatory visits for disease control in the second year of the disease course may significantly increase the risk of recorded MVA with profound injuries beyond 2 years after AS diagnosis.

Supplemental Material

sj-docx-1-tab-10.1177_1759720X241273039 – Supplemental material for Risk of motor vehicle accidents with profound injuries in patients with ankylosing spondylitis: a nationwide, population-based cohort study

Supplemental material, sj-docx-1-tab-10.1177_1759720X241273039 for Risk of motor vehicle accidents with profound injuries in patients with ankylosing spondylitis: a nationwide, population-based cohort study by Chung-Mao Kao, Wei-Li Ho, Yi-Ming Chen, Tsu-Yi Hsieh, Wen-Nan Huang, Yi-Hsing Chen and Hsin-Hua Chen in Therapeutic Advances in Musculoskeletal Disease

sj-docx-2-tab-10.1177_1759720X241273039 – Supplemental material for Risk of motor vehicle accidents with profound injuries in patients with ankylosing spondylitis: a nationwide, population-based cohort study

Supplemental material, sj-docx-2-tab-10.1177_1759720X241273039 for Risk of motor vehicle accidents with profound injuries in patients with ankylosing spondylitis: a nationwide, population-based cohort study by Chung-Mao Kao, Wei-Li Ho, Yi-Ming Chen, Tsu-Yi Hsieh, Wen-Nan Huang, Yi-Hsing Chen and Hsin-Hua Chen in Therapeutic Advances in Musculoskeletal Disease

We would like to thank Ms. Chia-Hui Yu and the Biostatistics Task Force of Taichung Veterans General Hospital for their assistance in performing the statistical analysis.

Declarations

ORCID iDs: Chung-Mao Kao https://orcid.org/0000-0002-2863-4350

Hsin-Hua Chen https://orcid.org/0000-0002-7304-4587

Supplemental material: Supplemental material for this article is available online.

Ethics approval and consent to participate: This study was approved by the Institutional Review Board(II) of Taichung Veterans General Hospital (IRB TCVGH number: CE17100B). Claims data of the databases were open to research, with all personal information encrypted. Since data were anonymised before analyses, informed consent from patients was allowed to be waived by Institutional Review Board(II) at 106-B-05 Board Meeting.

Consent for publication: Not applicable.

Author contributions: Chung-Mao Kao: Investigation; Writing – original draft; Writing – review & editing.

Wei-Li Ho: Conceptualization; Data curation; Funding acquisition; Investigation; Project administration.

Yi-Ming Chen: Data curation; Formal analysis; Investigation; Methodology.

Tsu-Yi Hsieh: Resources; Supervision; Validation; Visualization.

Wen-Nan Huang: Supervision; Validation; Visualization.

Yi-Hsing Chen: Supervision; Validation; Visualization.

Hsin-Hua Chen: Conceptualization; Formal analysis; Funding acquisition; Methodology; Project administration; Resources; Software; Supervision; Writing – review & editing.

Funding: The authors disclosed receipt of the following financial support for the research, authorship and/or publication of this article: This study was supported by Chiayi Branch, Taichung Veterans General Hospital, Taiwan, and by Taichung Veterans General Hospital, Taiwan (grant no. TCVGH-VHCY1138604, TCVGH-VHCY1088605).

The authors declare that there is no conflict of interest.

Availability of data and materials: The datasets used and/or analysed in the current study are available from the corresponding author upon reasonable request.
==== Refs
References

1 Burki V Gossec L Payet J , et al . Prevalence and characteristics of hip involvement in spondyloarthritis: a single-centre observational study of 275 patients. Clin Exp Rheumatol 2012; 30 (4 ): 481–486.22513241
2 Will R Kennedy G Elswood J , et al . Ankylosing spondylitis and the shoulder: commonly involved but infrequently disabling. J Rheumatol 2000; 27 (1 ): 177–182.10648036
3 Ara RM Packham JC Haywood KL . The direct healthcare costs associated with ankylosing spondylitis patients attending a UK secondary care rheumatology unit. Rheumatology (Oxford) 2008; 47 (1 ): 68–71.18077492
4 Bedaiwi M Sari I Thavaneswaran A , et al . Fatigue in ankylosing spondylitis and nonradiographic axial spondyloarthritis: analysis from a longitudinal observation cohort. J Rheumatol 2015; 42 (12 ): 2354–2360.26523020
5 Jones GT Mallawaarachchi B Shim J , et al . The prevalence of fibromyalgia in axial spondyloarthritis. Rheumatol Int 2020; 40 (10 ): 1581–1591.32556474
6 Leverment S Clarke E Wadeley A , et al . Prevalence and factors associated with disturbed sleep in patients with ankylosing spondylitis and non-radiographic axial spondyloarthritis: a systematic review. Rheumatol Int 2017; 37 (2 ): 257–271.27796520
7 Shen CC Hu LY Yang AC , et al . Risk of psychiatric disorders following ankylosing spondylitis: a nationwide population-based retrospective cohort study. J Rheumatol 2016; 43 (3 ): 625–631.26834219
8 Hsieh M-Y Kuo C-F . FRI0428 epidemiology of ankylosing spondylitis in Taiwan: a nationwide population study. Ann Rheum Dis 2016; 75 (Suppl. 2 ): 590–591.
9 van der Linden S Valkenburg HA Cats A . Evaluation of diagnostic criteria for ankylosing spondylitis. A proposal for modification of the New York criteria. Arthritis Rheum 1984; 27 (4 ): 361–368.6231933
10 Rudwaleit M Landewe R van der Heijde D , et al . The development of Assessment of SpondyloArthritis international Society classification criteria for axial spondyloarthritis (part I): classification of paper patients by expert opinion including uncertainty appraisal. Ann Rheum Dis 2009; 68 (6 ): 770–776.19297345
11 Rudwaleit M van der Heijde D Landewe R , et al . The development of Assessment of SpondyloArthritis international Society classification criteria for axial spondyloarthritis (part II): validation and final selection. Ann Rheum Dis 2009; 68 (6 ): 777–783.19297344
12 Wei JC Liu CH Tseng JC , et al . Taiwan Rheumatology Association consensus recommendations for the management of axial spondyloarthritis. Int J Rheum Dis 2020; 23 (1 ): 7–23.31777200
13 Packham J . Optimizing outcomes for ankylosing spondylitis and axial spondyloarthritis patients: a holistic approach to care. Rheumatology (Oxford) 2018; 57 (Suppl. 6 ): vi29–vi34.
14 Kiltz U van der Heijde D Boonen A , et al . Development of a health index in patients with ankylosing spondylitis (ASAS HI): final result of a global initiative based on the ICF guided by ASAS. Ann Rheum Dis 2015; 74 (5 ): 830–835.24399232
15 Chen HH Chen YM Lai KL , et al . Gender difference in ASAS HI among patients with ankylosing spondylitis. PLoS One 2020; 15 (7 ): e0235678.
16 Pai CW Lin HY Tsai SH , et al . Comparison of traffic-injury related hospitalisation between bicyclists and motorcyclists in Taiwan. PLoS One 2018; 13 (1 ): e0191221.
17 Chang YH Li CY Lu TH , et al . Risk of injury and mortality among driver victims involved in single-vehicle crashes in Taiwan: comparisons between vehicle types. Int J Environ Res Public Health 2020; 17 (13 ): 4687.32610689
18 Dagfinrud H Kjeken I Mowinckel P , et al . Impact of functional impairment in ankylosing spondylitis: impairment, activity limitation, and participation restrictions. J Rheumatol 2005; 32 (3 ): 516–523.15742446
19 Hsieh CY Su CC Shao SC , et al . Taiwan’s National Health Insurance Research Database: past and future. Clin Epidemiol 2019; 11 : 349–358.31118821
20 Curtis JR Harrold LR Asgari MM , et al . Diagnostic prevalence of ankylosing spondylitis using computerized health care data, 1996 to 2009: underrecognition in a US Health Care Setting. Perm J 2016; 20 (4 ): 15–151.
21 Kao CM Wang JS Ho WL , et al . Factors associated with the risk of major adverse cardiovascular events in patients with ankylosing spondylitis: a nationwide, population-based case–control study. Int J Environ Res Public Health 2022; 19 (7 ): 4098.35409780
22 von Elm E Altman DG Egger M , et al . The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Lancet 2007; 370 (9596 ): 1453–1457.18064739
23 Aouad K Ziade N Baraliakos X . Structural progression in axial spondyloarthritis. Joint Bone Spine 2020; 87 (2 ): 131–136.31067506
24 Poddubnyy D Listing J Haibel H , et al . Functional relevance of radiographic spinal progression in axial spondyloarthritis: results from the GErman SPondyloarthritis Inception Cohort. Rheumatology (Oxford) 2018; 57 (4 ): 703–711.29373733
25 Poddubnyy D Rudwaleit M Haibel H , et al . Rates and predictors of radiographic sacroiliitis progression over 2 years in patients with axial spondyloarthritis. Ann Rheum Dis 2011; 70 (8 ): 1369–1374.21622969
26 Dougados M Demattei C van den Berg R , et al . Rate and predisposing factors for sacroiliac joint radiographic progression after a two-year follow-up period in recent-onset spondyloarthritis. Arthritis Rheumatol 2016; 68 (8 ): 1904–1913.26990518
27 Schwartzman S Ruderman EM . A road map of the axial spondyloarthritis continuum. Mayo Clin Proc 2022; 97 (1 ): 134–145.34801248
28 Machado P Landewé R Braun J , et al . Both structural damage and inflammation of the spine contribute to impairment of spinal mobility in patients with ankylosing spondylitis. Ann Rheum Dis 2010; 69 (8 ): 1465–1470.20498215
29 Holden W Taylor S Stevens H , et al . Neck pain is a major clinical problem in ankylosing spondylitis, and impacts on driving and safety. Scand J Rheumatol 2005; 34 (2 ): 159–160.16095017
30 Tsao CH Huang JY Huang HH , et al . Ankylosing spondylitis is associated with risk of new-onset obstructive sleep apnea: a nationwide population-based cohort study. Front Med (Lausanne) 2019; 6 : 285.31867336
31 Yi E Ahuja A Rajput T , et al . Clinical, economic, and humanistic burden associated with delayed diagnosis of axial spondyloarthritis: a systematic review. Rheumatol Ther 2020; 7 (1 ): 65–87.31965538
32 Morton L Macfarlane GJ Jones G , et al . Driving difficulties in patients with axial spondyloarthritis: results from the Scotland Registry for Ankylosing Spondylitis. Arthritis Care Res (Hoboken) 2022; 74 (9 ): 1541–1549.33734612
33 Ramiro S van der Heijde D van Tubergen A , et al . Higher disease activity leads to more structural damage in the spine in ankylosing spondylitis: 12-year longitudinal data from the OASIS cohort. Ann Rheum Dis 2014; 73 (8 ): 1455–1461.24812292
34 Rudwaleit M Haibel H Baraliakos X , et al . The early disease stage in axial spondylarthritis: results from the German Spondyloarthritis Inception Cohort. Arthritis Rheum 2009; 60 (3 ): 717–727.19248087
35 Hollick RJ Stelfox K Dean LE , et al . Outcomes and treatment responses, including work productivity, among people with axial spondyloarthritis living in urban and rural areas: a mixed-methods study within a national register. Ann Rheum Dis 2020; 79 (8 ): 1055–1062.32522742
36 Ramiro S Landewe R van Tubergen A , et al . Lifestyle factors may modify the effect of disease activity on radiographic progression in patients with ankylosing spondylitis: a longitudinal analysis. RMD Open 2015; 1 (1 ): e000153.
37 Uslu S Gülle S Urak Ö , et al . Biological treatment in elderly and young patients with ankylosing spondylitis: TURKBIO real-life data results. Arch Rheumatol 2024; 39 (2 ): 232–241.38933720
38 Jadon DR Ramanan AV Sengupta R . Juvenile versus adult-onset ankylosing spondylitis—clinical, radiographic, and social outcomes. A systematic review. J Rheumatol 2013; 40 (11 ): 1797–1805.24037547
39 Lin YC Liang TH Chen WS , et al . Differences between juvenile-onset ankylosing spondylitis and adult-onset ankylosing spondylitis. J Chin Med Assoc 2009; 72 (11 ): 573–580.19948434
40 Keski-Filppula T Hakko H Raty E , et al . A short survival time after last psychiatric hospitalization in drivers with psychotic disorder killed in fatal motor vehicle accidents. Schizophr Res 2020; 216 : 235–242.31813802
41 Brubacher JR Chan H Erdelyi S , et al . Medications and risk of motor vehicle collision responsibility in British Columbia, Canada: a population-based case–control study. Lancet Public Health 2021; 6 (6 ): e374–e385.
42 Verster JC Volkerts ER . Antihistamines and driving ability: evidence from on-the-road driving studies during normal traffic. Ann Allergy Asthma Immunol 2004; 92 (3 ): 294–303; quiz 303–305, 355.15049392
43 Zhao SS Pittam B Harrison NL , et al . Diagnostic delay in axial spondyloarthritis: a systematic review and meta-analysis. Rheumatology (Oxford) 2021; 60 (4 ): 1620–1628.33428758
44 Pederson LL Koval J Vingilis E , et al . The relationship between motor vehicle collisions and cigarette smoking in Ontario: analysis of CAMH survey data from 2002 to 2016. Prev Med Rep 2019; 13 : 327–331.30792948
45 Hou CC Chen SC Tan LB , et al . Psychoactive substance use and the risk of motor vehicle crash injuries in southern Taiwan. Prev Sci 2012; 13 (1 ): 36–42.21845473
46 Holt PL . Stressful life events preceding road traffic accidents. Injury 1981; 13 (2 ): 111–115.7327723
