spinal cord injury ( sci ) impairs a patient 's motor , sensory , or autonomic system , and sci - related costs incurred by health care systems often burden society .
respiratory failure , cardiovascular dysfunction , thromboembolism , and autonomic dysreflexia are the common complications of sci .
furthermore , recent studies have reported a higher risk of developing insulin resistance , glucose intolerance , and lipid abnormalities in sci patients than in able - bodied people .
supporting a patient severely affected by sci can cost up to us$1 million in the first year . during the past decades
up to 250,000 people worldwide sustain a sci every year , and their life expectancy continues to increase .
peripheral arterial disease ( pad ) is one of the most lethal diseases but is frequently neglected . even without previous ischemic stroke or myocardial infarction ,
patients with pad are at an equal risk of death as are patients with cardiovascular disease ( cvd ) or previous coronary or cerebrovascular disease .
a meta - analysis suggested that pad affects more than 202 million people worldwide . resting pain , claudication , and atypical leg discomfort are the symptoms of pad ; however , up to 50% of patients with pad are asymptomatic .
the confirmed risk factors of pad include diabetes mellitus ( dm ) , smoking , hypertension ( htn ) , and dyslipidemia .
patients with sci are at a greater risk of cvd and a relevant higher risk of death than able - bodied persons .
recent studies have advocated cvd as the leading cause of mortality in patients with sci .
several studies have proposed subclinical atherosclerotic markers for patients with sci including the carotid intima - media thickness and ankle - brachial index .
the relationship between sci and pad requires further large - scale investigation because pad is an atherosclerotic process affecting the noncoronary arterial system .
our retrospective nationwide cohort study investigated the relationship between sci and pad by using data from a national health insurance database .
data used in this study were obtained from the taiwan national health insurance research database ( nhird ) .
the national health insurance ( nhi ) program is a government - run , single - payer insurance system that was established in 1995 in taiwan .
this nhi program covers over 99% of the 23.74 million residents of taiwan ( http://www.nhi.gov.tw ) .
the national health research institute ( nhri ) audits and releases the nhird data for use in health service studies . in accordance with the personal information protection act , all patient reimbursement data are deidentified and linked to a patient identification number before being released for academic research . for this study , we used a subset of the nhird containing health care data including files of inpatient claims and the registry of beneficiaries .
all clinical diagnoses were recorded according to the international classification of diseases , ninth revision , clinical modification [ icd-9-cm ] codes . in our national health insurance ( nhi )
program , all insurance claims are under scrutiny of medical reimbursement specialists and anonymous peer reviews .
the definition of sci and pad was based on icd-9-cm codes determined by clinical physicians after strict inspections in the reimbursement process based on laboratory , imaging , and pathological data
. moreover , there are even severe penalties for physicians if inappropriate icd-9-cm codes were documented in clinical records .
therefore , the diagnoses or definition of sci , and pad are accurate and reliable .
the study was exempted from full review by the research ethics committee of china medical university and hospital ( cmuh104-rec2 - 105 ) .
we selected adult patients with a first diagnosis of sci ( icd-9-cm codes 806 and 952 ) from 2000 to 2010 as the sci cohort .
cervical spine ( c - spine ) ( icd-9-cm codes 806.0 , 806.1 , 952.0 , 952.00 , 952.01 , 952.02 , 952.03 , 952.04 , 952.05 , 952.06 , 952.07 , 952.08 , and 952.09 ) ; 2 .
thoracic spine ( t - spine ) ( icd-9-cm codes 806.2 , 806.21 , 806.26 , 806.3 , 952.1 , 952.11 , and 952.16 ) ; 3 . lumbar , sacral , or coccygeal spine ( l - s - co spine ) ( icd-9-cm codes 806.4 , 806.5 , 806.6 , 806.7 , 806.8 , 806.9 , 952.2 , 952.3 , 952.4 , 952.8 , and 952.9 ) ; 4 .
multiple spine sci ( any 2 or more than 2 lesions in the c - spine , t - spine , and l - s - co spine ) .
the exclusion criteria were an age younger than 20 years , incomplete age or sex information , and a history of pad ( icd-9-cm codes 440.2 , 440.3 , 440.8 , 440.9 , 443 , 444.22 , 444.8 , 447.8 , and 447.9 ) at the baseline .
a non - sci cohort was randomly selected from the nhi beneficiaries aged 20 years and older and frequency - matched for age ( every 5 years ) , sex , index year , and comorbidities , including diabetes , hypertension , hyperlipidemia , chronic obstructive pulmonary disease ( copd ) , heart failure , obesity , cad , stroke , and asthma with the sci cohort at a 4:1 ratio , and was subjected to the same exclusion criteria .
the main outcome was based on the admission claims data of pad diagnoses during follow - up .
each patient was followed up from the index date to the date of pad occurrence , withdrawal from the nhi , or december 31 , 2011 .
baseline comorbidities including dm ( icd-9-cm codes 250 ) , htn ( icd-9-cm codes 401405 ) , hyperlipidemia ( icd-9-cm codes 272 ) , copd ( icd-9-cm codes 491 , 492 , and 496 ) , heart failure ( icd-9-cm codes 428 ) , obesity ( icd-9-cm code 278 ) , coronary artery disease ( cad ) ( icd-9-cm codes 410414 ) , stroke ( icd-9-cm codes 430438 ) , and asthma ( icd-9-cm code 493 ) were identified according to diagnoses in hospitalization records before the index date .
the chi - squared test was used to examine and compare the distributions of the categorical variables of the sci and the non - sci cohort .
the student 's t test was used to examine the mean ages and mean follow - up times of the cohorts .
we estimated the cumulative incidence by using the kaplan  meier method , and the log - rank test was used to compare the cumulative incidence curves of the sci and non - sci cohorts .
the overall sex- , age- , and comorbidity - specific incidence densities of pad were measured for each cohort .
univariable and multivariable cox proportional hazards regression models were used to examine the association between sci and the risk of pad , which is expressed as a hazard ratio ( hr ) with a 95% confidence interval ( ci ) .
we selected the comorbidities as variables in multivariable analysis because they were statistically significant in the univariable model .
we further tested the interaction between sex and sci , age and sci , and between comorbidity and sci by including a cross - product term in the model .
we used the scaled schoenfeld residuals for testing the proportional hazard model assumption . because the proportional hazard model assumption was violated ( p  =  0.005 ) , we stratified the follow - up duration to deal with the violation of proportional hazard assumption .
further analysis was performed to assess the variations in the association of pad at different levels of sci .
data management and descriptive analyses were performed using the sas 9.2 statistical package ( sas institute inc . ,
the chi - squared test was used to examine and compare the distributions of the categorical variables of the sci and the non - sci cohort .
the student 's t test was used to examine the mean ages and mean follow - up times of the cohorts .
we estimated the cumulative incidence by using the kaplan  meier method , and the log - rank test was used to compare the cumulative incidence curves of the sci and non - sci cohorts .
the overall sex- , age- , and comorbidity - specific incidence densities of pad were measured for each cohort .
univariable and multivariable cox proportional hazards regression models were used to examine the association between sci and the risk of pad , which is expressed as a hazard ratio ( hr ) with a 95% confidence interval ( ci ) .
we selected the comorbidities as variables in multivariable analysis because they were statistically significant in the univariable model .
we further tested the interaction between sex and sci , age and sci , and between comorbidity and sci by including a cross - product term in the model .
we used the scaled schoenfeld residuals for testing the proportional hazard model assumption . because the proportional hazard model assumption was violated ( p  =  0.005 ) , we stratified the follow - up duration to deal with the violation of proportional hazard assumption .
further analysis was performed to assess the variations in the association of pad at different levels of sci .
data management and descriptive analyses were performed using the sas 9.2 statistical package ( sas institute inc
this study included 42,673 with sci and 170,389 patients without sci . in our study , men constituted the majority ( 63.2% versus 36.8% ) , and almost 71.7% of the patients were more than 50 years old .
the mean age of the sci cohort was 52.4    18.2 years and that of the non - sci cohort was 52.1    18.4 years .
the sci cohort and non- sci cohort were well matched for comorbidities , except for obesity . during the mean follow - up of 5.66 years for the sci cohort and 6.03 years for the non - sci cohort ,
the cumulative incidence of pad was significantly higher for patients in the sci cohort than for those in the non - sci cohort ( log - rank test p  <  0.001 ) ( figure 1 ) .
the overall incidence of pad was 29% higher in the sci cohort than in the non - sci cohort ( 1.68 versus 1.30 per 1000 person - years ) with an adjusted hr of 1.37 ( 95% ci  =  1.221.53 ) ( table 1 ) .
the sex - specific adjusted hr of pad for the sci and non - sci cohorts was significant for women ( hr  =  1.29 , 95% ci  =  1.071.53 ) and men ( hr  =  1.42 , 95% ci  =  1.231.65 ) .
the pad incidence increased with age in both cohorts . however , the age - specific relative risk of pad was the greatest for the youngest age group among the sci and non - sci cohorts ( 50 years , adjusted hr  =  4.02 ; 95% ci  =  2.227.26 ) .
the corresponding adjusted hr decreased to 1.23 ( 95% ci , 1.081.42 ) for the oldest age group ( > 65 years ) , although this group had the highest incidence of pad .
the pad incidence rate was greater in patients with comorbidities than in patients without comorbidity in both the cohorts . a significantly higher risk of pad was observed in sci patients without comorbidities ( hr  =  1.87 ; 95% ci  =  1.502.33 ) than in the non - sci patients without comorbidities , respectively . in the first year of follow - up , the sci cohort had a higher risk of pad than the non - sci cohort ( adjusted hr  =  1.98 , 95% ci  =  1.532.57 ) . moreover , the risk of pad in the sci cohort was still significantly higher than that in the non - sci cohort after 5 years of follow - up ( adjusted hr  =  1.36 , 95% ci  =  1.121.65 ) .
meier curve for the cumulative incidence of pad in the sci and non - sci cohorts during the 12-year follow - up . incidence ratio and hr of pad
according to sci status stratified by demographic factors and comorbidities the results of the univariable and multivariable cox proportional hazards regression models for analyzing the risk of variables contributing to pad are summarized in table 2 .
the adjusted hr of pad was increased 1.42-fold for men relative to women ( 95% ci  =  1.291.73 ) and increased 1.06-fold ( 95%
the risk of pad was greater in patients with comorbidities , namely dm ( adjusted hr  =  3.11 , 95% ci  =  2.803.44 ) , htn ( adjusted hr  =  1.50 , 95% ci  =  1.341.68 ) , hyperlipidemia ( adjusted hr  =  1.49 , 95% ci  =  1.291.73 ) , heart failure ( adjusted hr  =  1.81 , 95% ci  =  1.542.12 ) , obesity ( adjusted hr  =  2.33 , 95% ci  =  1.045.21 ) , cad ( adjusted hr  =  1.16 , 95% ci  =  1.021.31 ) , and stroke ( adjusted hr  =  1.57 , 95% ci  =  1.351.82 ) .
hr of pad in association with sex , age , and comorbidities in univariable and multivariable cox regression models table 3 lists the relative risk and hr of pad associated with different levels of sci lesion . compared with the non - sci cohort ,
patients with l - s - co - spine injury were at a 56% ( adjusted hr  =  1.56 , 95% ci  =  1.331.84 ) higher risk of pad , respectively .
patients with multiple sci lesions had a 2.11-fold higher risk of pad than the non - sci cohort ( 95% ci  =  1.592.79 ) .
this study was conducted in a region with a high prevalence of pad . in our study , comorbidities were significantly more prevalent in the sci group than in the non - sci group .
after adjustment for sex , age , and the aforementioned comorbidities , the sci group exhibited a 1.37-fold higher risk of pad than the non - sci group .
the overall incidence of pad was 29% higher in the sci cohort than in the non - sci cohort . during the 12-year follow - up ,
a significantly higher risk of pad was noted within the first year . in our study , l - s - co - spine and multiple spine sci were significantly associated with an increased risk of pad . a historical prospective study among patients surviving at least 20 years with sci stated that the risk of all cvd increased with severity of sci . the cvd in this study defined by icd/9 codes 390 to 448 ;
furthermore , the number of patients with pad in our study might have been underestimated because pad is commonly underdiagnosed .
hence , the true impact of sci on pad might be stronger than that reported here .
first , the patients with sci in our study had a significantly higher incidence of comorbidities than did the non - sci patients ( table 4 ) .
moreover , htn , dm , and hyperlipidemia were observed to be the major risk factors for pad .
thus , patients with these comorbidities possibly exhibited a higher risk of pad than did patients without comorbidities .
second , previous research has reported a susceptibility of patients with sci to metabolic disorders including carbohydrate and dyslipidemia .
lee et al reported that a quarter of the sci population presented with metabolic syndromes and insulin resistance .
the impact of sci on the metabolism provides evidence of the association between sci and subsequent pad .
third , impaired control of vessels and altered vascular reactivity in patients with sci possibly contribute to atherosclerosis in the lower extremities .
fourth , a sedentary lifestyle and subsequent loss of lean body mass , and fluctuations in the blood pressure are possible mechanisms underlying pad observed in sci patients .
demographic characteristics and comorbidities in study patients according to sci status the strengths of this study are that it is the first cohort study focusing on pad in sci patients with uniform data collection and a sufficiently large sample size to enable meaningful analyses . nevertheless , some limitations must be noted .
first , the diagnoses recorded in the nhird are not validated for academic purposes . moreover , the nhird diagnoses were documented using icd-9-cm codes , and data on environmental risk factors influencing pad such as smoking habits and family history could not be obtained .
second , despite the controls and adjustment for interfering factors , we probably did not completely control for the confounding effects of preexisting comorbidities of pad , potentially leading to an inaccurate estimation of the relationship between sci and pad risk .
third , details of severity of pad could not be obtained ; thus , further survey in our study was limited .
fourth , we have tried to minimize the basic difference between the sci and non - sci cohorts ; however , the confounding variables of smoking , blood pressure , and cholesterol level were not available in our database .
thus , the effect of residual confounding could not be completely excluded . in conclusion , patients with sci are associated with an increased risk of pad , regardless of preexisting comorbidities .
our findings improve physicians awareness of the risk of pad in patients with sci and facilitate developing strategies to prevent , detect , and manage pad in patients with sci .
 abstractthe aim of this study was to elucidate the relationship between spinal cord injury ( sci ) and the risk of peripheral arterial disease ( pad ) in a cohort study with a large representative sample.the national health insurance database was used to select patients who were diagnosed from 2000 to 2010 . 
 patients with a history of pad were excluded . 
 the sci group comprised 42,673 patients diagnosed with sci , and we enrolled 170,389 matched controls ( non - sci group ) . 
 we used a cox proportional hazards regression model to analyze the adjusted risk of pad between the case and control patients.patients with sci exhibited a significantly higher risk ( hazard ratio [ hr ]  =  1.37 ; 95% confidence interval [ ci ]  =  1.221.53 ) of pad than patients without sci . 
 patients with diabetes were at the highest risk of developing pad ( adjusted hr  =  3.11 , 95% ci  =  2.803.44 ) . among patients without comorbidity , sci patients exhibited a significantly higher risk of pad than non - sci patients . 
 furthermore , lumbar , sacral , or coccygeal spine , and multiple spine sci were significantly associated with an increased risk of pad ( hr  =  1.56 , 95% ci  =  1.331.84 , hr  =  2.11 , 
 95% ci  =  1.592.79 , respectively).sci is associated with an increased risk of pad . 
 future studies should focus on modifying risk factors to reduce pad risk among patients with sci . 
INTRODUCTION
METHODS
SAMPLED PARTICIPANTS
OUTCOME AND COMORBIDITIES
Statistical Analysis
RESULTS
DISCUSSION