Acessibilidade / Reportar erro

Cerebrovascular risk factors and their time-dependent effects on stroke survival in the EMMA cohort study

Abstract

To investigate the time-dependent effects of traditional risk factors on functional disability in all-cause mortality post-stroke, we evaluated data from a long-term stroke cohort. Baseline cerebrovascular risk factors (CVRF) and functionality at 1 and 6 months were evaluated in survivors from a prospective stroke cohort using the modified Rankin scale (m-RS), which classifies participants as improvement of disability, unchanged disability (at least moderate), and worsening disability. Cox regression models considering baseline risk factors, medication use, and functionality 6 months after stroke were fitted to identify their time-dependent effects up to 12 years of follow-up. Adjusted hazard ratios (HR) with 95% confidence intervals (CI) are presented. Among 632 survivors (median age 68, 54% male, 71% first-ever episode), age and functional disability (unchanged and worsening) 6 months after ischemic stroke had time-dependent effects on all-cause mortality risk up to 12 years of follow-up. The most impacting risk factors were unchanged (at least moderate) (HR, 2.99; 95%CI: 1.98-4.52) and worsening disability (HR, 2.85; 95%CI: 1.26-6.44), particularly in the first two years after a stroke event (Time 1: ≥6 mo to <2.5 y). Worsening disability also impacted mortality in the period from ≥2.5 to <7.5 years (Time 2) of follow-up (HR, 2.43 (95%CI: 1.03-5.73). Other baseline factors had a fixed high-risk effect on mortality during follow-up. Post-stroke and continuous medication use had a fixed protective effect on mortality. Functional disability was the main contributor with differential risks of mortality up to 12 years of follow-up.

Stroke epidemiology; Cerebrovascular risk factors; Stroke prevention


Introduction

Traditionally, the analysis of risk factors and their effects in longitudinal epidemiological studies is based on fixed baseline data from a single time measurement, regardless of the length of follow-up. However, the effect of a risk factor on long-term mortality may vary over time, resulting in a weakening or strengthening of associations. In other words, even for risk factors measured only once (baseline), their effects can be strong in the long term but not necessarily in the short term, and vice-versa (11. Dekker FW, de Mutsert R, van Dijk PC, Zoccali C, Jager KJ. Survival analysis: time-dependent effects and time-varying risk factors. Kidney Int 2008; 74: 994-997, doi: 10.1038/ki.2008.328.
https://doi.org/10.1038/ki.2008.328...
).

The proper analysis of effects over time should be driven by a clear research question. Two research questions, one for time-dependent effects (short-term versus long-term effects) of a fixed baseline risk factor and one for time-dependent risk factors (time-varying risk factors), can be analyzed with time-dependent Cox regression analysis. In this context, even sociodemographic determinants, especially low educational level, cerebrovascular risk factors, such as atrial fibrillation (AF), diabetes, hypertension, presence of multimorbidities, post-stroke depression, or disability, can have different impacts over time and might determine different mortality risks (22. Goulart AC, Olmos RD, Santos IS, Tunes G, Alencar AP, Thomas N, et al. The impact of atrial fibrillation and long-term oral anticoagulant use on all-cause and cardiovascular mortality: a 12-year evaluation of the prospective Brazilian Study of Stroke Mortality and Morbidity. Int J Stroke 2022; 17: 48-58, doi: 10.1177/1747493021995592.
https://doi.org/10.1177/1747493021995592...
-3. Tang Z, Zhou T, Luo Y, Xie C, Huo D, Tao L, et al. Risk factors for cerebrovascular disease mortality among the elderly in Beijing: a competing risk analysis. PLoS One 2014 ; 9: e87884, doi: 10.1371/journal.pone.0087884.
https://doi.org/10.1371/journal.pone.008...
4. Castro HHG, Alencar AP, Benseãor IM, Lotufo PA, Goulart AC. Multimorbidities Are associated to lower survival in ischaemic stroke: results from a Brazilian Stroke Cohort (EMMA Study). Cerebrovasc Dis 2017; 44: 232-239, doi: 10.1159/000479827.
https://doi.org/10.1159/000479827...
5. Goulart AC, Fernandes TG, Santos IS, Alencar AP, Benseãor IM, Lotufo PA. Predictors of long-term survival among first-ever ischemic and hemorrhagic stroke in a Brazilian stroke cohort. BMC Neurol 2013; 13: 51, doi: 10.1186/1471-2377-13-51.
https://doi.org/10.1186/1471-2377-13-51...
6. de Mello RF, Santos IS, Alencar AP, Benseãor IM, Lotufo PA, Goulart AC. Major depression as a predictor of poor long-term survival in a Brazilian Stroke Cohort (Study of Stroke Mortality and Morbidity in Adults) EMMA study. J Stroke Cerebrovasc Dis 2016; 25: 618-625, doi: 10.1016/j.jstrokecerebrovasdis.2015.11.021.
https://doi.org/10.1016/j.jstrokecerebro...
77. Hawkes MA, Blaginykh E, Ruff MW, Burrus T, Wijdicks EFM, Rabinstein AA. Long-term mortality, disability and stroke recurrence in patients with basilar artery occlusion. Eur J Neurol 2020; 27: 579-585, doi: 10.1111/ene.14126.
https://doi.org/10.1111/ene.14126...
). Despite the importance of analyzing the real effect of risk factors on mortality, particularly in long-term prospective cohorts, most previous studies (22. Goulart AC, Olmos RD, Santos IS, Tunes G, Alencar AP, Thomas N, et al. The impact of atrial fibrillation and long-term oral anticoagulant use on all-cause and cardiovascular mortality: a 12-year evaluation of the prospective Brazilian Study of Stroke Mortality and Morbidity. Int J Stroke 2022; 17: 48-58, doi: 10.1177/1747493021995592.
https://doi.org/10.1177/1747493021995592...
-3. Tang Z, Zhou T, Luo Y, Xie C, Huo D, Tao L, et al. Risk factors for cerebrovascular disease mortality among the elderly in Beijing: a competing risk analysis. PLoS One 2014 ; 9: e87884, doi: 10.1371/journal.pone.0087884.
https://doi.org/10.1371/journal.pone.008...
4. Castro HHG, Alencar AP, Benseãor IM, Lotufo PA, Goulart AC. Multimorbidities Are associated to lower survival in ischaemic stroke: results from a Brazilian Stroke Cohort (EMMA Study). Cerebrovasc Dis 2017; 44: 232-239, doi: 10.1159/000479827.
https://doi.org/10.1159/000479827...
5. Goulart AC, Fernandes TG, Santos IS, Alencar AP, Benseãor IM, Lotufo PA. Predictors of long-term survival among first-ever ischemic and hemorrhagic stroke in a Brazilian stroke cohort. BMC Neurol 2013; 13: 51, doi: 10.1186/1471-2377-13-51.
https://doi.org/10.1186/1471-2377-13-51...
6. de Mello RF, Santos IS, Alencar AP, Benseãor IM, Lotufo PA, Goulart AC. Major depression as a predictor of poor long-term survival in a Brazilian Stroke Cohort (Study of Stroke Mortality and Morbidity in Adults) EMMA study. J Stroke Cerebrovasc Dis 2016; 25: 618-625, doi: 10.1016/j.jstrokecerebrovasdis.2015.11.021.
https://doi.org/10.1016/j.jstrokecerebro...
77. Hawkes MA, Blaginykh E, Ruff MW, Burrus T, Wijdicks EFM, Rabinstein AA. Long-term mortality, disability and stroke recurrence in patients with basilar artery occlusion. Eur J Neurol 2020; 27: 579-585, doi: 10.1111/ene.14126.
https://doi.org/10.1111/ene.14126...
) demonstrated the effect of risk factors without considering time-stratified effects or time-dependent effects of fixed baseline risk factors (risk factors measured once, usually at baseline, considering effects analysis in different times during follow-up) or time-varying risk factors (risk factors measured multiple times during follow-up) (11. Dekker FW, de Mutsert R, van Dijk PC, Zoccali C, Jager KJ. Survival analysis: time-dependent effects and time-varying risk factors. Kidney Int 2008; 74: 994-997, doi: 10.1038/ki.2008.328.
https://doi.org/10.1038/ki.2008.328...
).

Therefore, we aimed to evaluate time-dependent effects of traditional risk factors (measured at baseline) and post-stroke disability (measured 1 and 6 months after stroke) in the short- and long-term all-cause mortality in a Brazilian stroke cohort, the Study of Stroke Mortality and Morbidity (EMMA). EMMA is an ongoing longitudinal study that investigates the impact of cardiovascular risk factors and the role of secondary prevention on long-term mortality among individuals afflicted by stroke living in a low-middle socioeconomic income (LMIC) community in the city of São Paulo, Brazil.

Material and Methods

Study design and population

Study subjects were participants of the EMMA cohort, a well-characterized prospective long-term stroke surveillance cohort that has been ongoing since 2006 (88. Goulart AC, Bustos IR, Abe IM, Pereira AC, Fedeli LM, Benseãor IM, et al. A stepwise approach to stroke surveillance in Brazil: the EMMA (Estudo de Mortalidade e Morbidade do Acidente Vascular Cerebral) study. Int J Stroke 2010; 5: 284-289, doi: 10.1111/j.1747-4949.2010.00441.x.
https://doi.org/10.1111/j.1747-4949.2010...
). The EMMA study was based on the Stepwise Approach to Stroke Surveillance (STEPS Stroke)-World Health Organization (WHO) (99. World Health Organization. WHO STEPS stroke manual: the WHO STEPwise approach to stroke surveillance/noncommunicable diseases and mental health. World Health Organization; 2006. [cited 12/06/2010]. Available from: World Health Organization http://www.who.int/chp/steps/stroke/manual/en/index.html.
http://www.who.int/chp/steps/stroke/manu...
,1010. Lotufo PA, Benseãor IM. Improving WHO STEPS Stroke in Brazil. Lancet Neurol 2007; 6: 387-388, doi: 10.1016/S1474-4422(07)70091-1.
https://doi.org/10.1016/S1474-4422(07)70...
). All patients older than 18 years with symptoms of acute stroke admitted to a single center, the Emergency Department (ED) of the HU-USP (University Hospital of the Universidade de São Paulo, Brazil), were invited to participate in the hospital phase of the EMMA study (88. Goulart AC, Bustos IR, Abe IM, Pereira AC, Fedeli LM, Benseãor IM, et al. A stepwise approach to stroke surveillance in Brazil: the EMMA (Estudo de Mortalidade e Morbidade do Acidente Vascular Cerebral) study. Int J Stroke 2010; 5: 284-289, doi: 10.1111/j.1747-4949.2010.00441.x.
https://doi.org/10.1111/j.1747-4949.2010...
). The HU-USP is a secondary community hospital located in a low-income area of approximately 500,000 inhabitants on the west side of São Paulo city. Of note, all potential candidates for reperfusion therapies are regularly transferred to our referral center (Stroke Unit Care, Department of Neurology, Hospital das Clínicas of the Universidade de São Paulo), located five miles from the HU-USP. Thus, these cases were not included in these analyses. Further information about the EMMA study has been published elsewhere (88. Goulart AC, Bustos IR, Abe IM, Pereira AC, Fedeli LM, Benseãor IM, et al. A stepwise approach to stroke surveillance in Brazil: the EMMA (Estudo de Mortalidade e Morbidade do Acidente Vascular Cerebral) study. Int J Stroke 2010; 5: 284-289, doi: 10.1111/j.1747-4949.2010.00441.x.
https://doi.org/10.1111/j.1747-4949.2010...
).

Between April 2006 and September 2014, 1,378 stroke cases (85.8% ischemic strokes [IS]) were enrolled in the EMMA cohort. For the present analyses, we excluded 195 hemorrhagic strokes, 304 participants with missing information on the modified Ranking Scale (mRS) at 1 and 6 months, 225 who died before the follow-up (6 months after stroke), and 22 without information on AF. Thus, the final sample comprised 632 participants.

Written informed consent was obtained from all EMMA participants or their legal representatives/caregivers (usually a close family member). The study was approved by the HU-USP Ethical Committee (# 593/05).

EMMA data collection

All data collection was performed by a trained interviewer (nurse) at hospital admission (personal interview) at 1 and 6 months (telephone interview) according to the WHO STEPS Stroke Manual instructions. Quality control was assured for cross-checking information, which was done by three medical coordinators of the EMMA study (88. Goulart AC, Bustos IR, Abe IM, Pereira AC, Fedeli LM, Benseãor IM, et al. A stepwise approach to stroke surveillance in Brazil: the EMMA (Estudo de Mortalidade e Morbidade do Acidente Vascular Cerebral) study. Int J Stroke 2010; 5: 284-289, doi: 10.1111/j.1747-4949.2010.00441.x.
https://doi.org/10.1111/j.1747-4949.2010...
,99. World Health Organization. WHO STEPS stroke manual: the WHO STEPwise approach to stroke surveillance/noncommunicable diseases and mental health. World Health Organization; 2006. [cited 12/06/2010]. Available from: World Health Organization http://www.who.int/chp/steps/stroke/manual/en/index.html.
http://www.who.int/chp/steps/stroke/manu...
). As an extension of WHO STEPS stroke, we evaluated EMMA participants yearly or until death for up to 12 years from baseline. All losses before 12 years were censored on the date they were last contacted or on which information about their vital status (alive) was confirmed by telephone contact or electronic hospital registers. Mortality data were confirmed by official death certificates in collaboration with the city of São Paulo's health statistics system (PRO-AIM, Program for Improvement of Mortality Information in the Municipality of São Paulo), State Health Offices (SEADE Foundation, São Paulo State Healthcare Data Analysis System), and the Brazilian Ministry of Health.

Stroke definition

Stroke was defined according to WHO criteria as “a focal (or at times global) neurological impairment of sudden onset lasting more than 24 h (or leading to death) of presumed vascular origin. This clinical definition has, therefore, four components: 1) a neurological impairment or deficit; 2) sudden onset; 3) lasting more than 24 h (or leading to death); and 4) of presumed vascular origin” (99. World Health Organization. WHO STEPS stroke manual: the WHO STEPwise approach to stroke surveillance/noncommunicable diseases and mental health. World Health Organization; 2006. [cited 12/06/2010]. Available from: World Health Organization http://www.who.int/chp/steps/stroke/manual/en/index.html.
http://www.who.int/chp/steps/stroke/manu...
). In addition to these four components, stroke diagnosis was validated by a medical doctor from the study and supported by no-contrast computed tomography (CT) performed within 24-48 of hospital admission for inclusion as a stroke case in the EMMA study (88. Goulart AC, Bustos IR, Abe IM, Pereira AC, Fedeli LM, Benseãor IM, et al. A stepwise approach to stroke surveillance in Brazil: the EMMA (Estudo de Mortalidade e Morbidade do Acidente Vascular Cerebral) study. Int J Stroke 2010; 5: 284-289, doi: 10.1111/j.1747-4949.2010.00441.x.
https://doi.org/10.1111/j.1747-4949.2010...
). Furthermore, stroke diagnosis was classified according to the International Classification of Diseases, 10th Edition (ICD-10: I60-I63.9). Here, we analyzed individuals with IS (code I63.X). Of note, history of stroke was based on information from the patient, caregiver, or hospital records.

Sociodemographic and clinical data

At baseline (hospital admission due to stroke), we collected information about socio-demographic factors (age, education, marital status, and race as self-reported skin color: white, mixed, black, and yellow), stroke episode (first-ever or recurrent stroke), pre-existing cerebrovascular risk factors (CVRF) such as smoking, alcohol consumption, and co-morbidities such as hypertension, diabetes mellitus, dyslipidemia, heart failure, coronary heart disease, and atrial fibrillation (AF). Data about clinical co-morbidities (before index event) and cardiovascular medication were validated by two senior medical researchers, based on medical registries. We considered all clinical diagnoses based on medical history and/or the use of medications. Diagnosis of AF at baseline was defined by ECG tracings (before or at hospital admission). Data on medication use such as anti-hypertensive, lipid-lowering, anti-diabetic, anti-platelet, and anti-coagulants were collected at hospital admission (study baseline) and 1 month after hospital discharge due to stroke. The degree of functional disability was evaluated by the modified Rankin Scale (mRS) at 1 and 6 months after stroke. Stroke functionality by mRS was classified as improvement of disability, unchanged disability (at least moderate), or worsening disability at 6 months after stroke compared to 1 month after stroke. In brief, the m-RS ranges from 0 to 6, with higher scores indicating greater disability and 0-2 is generally considered a good outcome with individuals assuming complete functional independence and 6 corresponds to an individual who is deceased (1111. Bamford JM, Sandercock PA, Warlow CP and Slattery J. Interobserver agreement for the assessment of handicap in stroke patients. Stroke 1989; 20: 828, doi: 10.1161/01.STR.20.6.828.
https://doi.org/10.1161/01.STR.20.6.828...
,1212. Dennis M, Mead G, Doubal F, Graham C. Determining the modified Rankin Score after stroke by postal and telephone questionnaires. Stroke 2012; 43: 851-853, doi: 10.1161/STROKEAHA.111.639708.
https://doi.org/10.1161/STROKEAHA.111.63...
). For the present analysis, we considered only survivors 6 months after the index event, thus, we reported the m-RS evaluation from 0 to 5 points at 1 and 6 months.

Statistical analysis

Categorical variables were analyzed by the chi-squared test and reported as absolute and relative frequencies. As continuous variables did not have a normal distribution, they were analyzed by the Wilcoxon test and reported as median values with respective interquartile ranges (IQR). Sociodemographic and clinical data and use of medication for clinical chronic conditions (anti-hypertensive, anti-diabetic, lipid-lowering medication, anti-platelet, and anti-coagulants) at hospital admission due to stroke, as well as functional status by mRS (at 1 and 6 months) were analyzed by sex.

For the all-cause mortality outcome, we calculated hazard ratios (HR) with respective 95% confidence intervals (CIs) derived from multivariable models adjusted by potential predictors: age, sex, education attainment, marital status, hypertension, diabetes, dyslipidemia, smoking, recurrent stroke (at baseline), functional disability up to 6 months, and medication use according to index event (pre-stroke, post-stroke, and continuous use).

Kaplan-Meier survival curve (1313. Kaplan EL, Meier P. Nonparametric Estimation from incomplete observations. J Am Stat Assoc 1958; 53 : 457-481, doi: 10.1080/01621459.1958.10501452.
https://doi.org/10.1080/01621459.1958.10...
) and Cox regression models (1414. Cox DR. Regression model and life tables. J Royal Stat Soc 1972; 34: 187-220.) were fitted to the data. We first fitted usual Cox models that did not include time-dependent effects. Analysis of Schoenfeld residuals and hypothesis testing for proportional hazards indicated that hazards were not proportional for the variables post-stroke functionality (P-value=0.011) and age (P-value=0.035) (1515. Grambsch P, Therneau T. Proportional hazards tests and diagnostics based on weighted residuals. Biometrika 1994; 81: 515-526, doi: 10.1093/biomet/81.3.515.
https://doi.org/10.1093/biomet/81.3.515...
). We then included time-dependent effects for these two variables and time change points of effects were chosen based on the plot of Schoenfeld residuals. Since Schoenfeld residuals reveal the functional form of the time-dependent effect, if the tendency line of the plot is close to a horizontal line, then there is evidence that a usual Cox model is appropriate; however, the lines in plots were very distant from a horizontal straight line and they suggested that a reasonable approximation of this time-dependent effect would be a step-function with three breaks: Time 1: ≥6 months to <2.5 years, Time 2: ≥2.5 to <7.5 years, and Time 3: ≥7.5 years or longer (Supplementary Figure S1)

Statistical analyses were performed with the statistical software SPSS version 27.0 (IBM, USA) and the R software (R Core Team). For all analyses, P-values less than 0.05 were considered significant.

Results

In the 632 IS survivors (median age 68 years; 71% first-ever stroke, 54% male) investigated, we found high frequencies of CVRF such as hypertension (94.1%), diabetes (58.1%), dyslipidemia (61.1%), and smoking (current 30.7%, past 14.1%) and relatively low continuous use of medications for controlling these risk factors (49.8%) (Table 1). Most individuals in this sample (57.8%) had mild or no disability; however, 25% had no improvement in their functionality (at least moderate) since the acute event, regardless of sex (Figure 1).

Table 1
Baseline characteristics of the 632 participants from the Study of Stroke Mortality and Morbidity (EMMA) study.
Figure 1
Functionality by the modified Ranking scale at 6 months after stroke in the Study of Stroke Mortality and Morbidity (EMMA) cohort.

Among 275 deaths attested along the follow-up, 22.5, 32, and 45.5% were reported as cerebrovascular, cardiovascular, and non-cardiovascular basic causes, respectively. About 56.5% (357/632) of stroke participants survived up to the 12 years of follow-up. Mortality did not differ between sexes, with P log rank=0.96 (Figure 2). Of all risk factors, age and functional disability (unchanged and worsening) had time-dependent effects on mortality risk. However, the most impacting risk factor was functional disability, particularly in the first two years (Time 1: ≥6 months to <2.5 years) after the stroke event; HR=2.99 (95%CI: 1.98-4.52) for unchanged disability (at least moderate) and HR=2.85 (95%CI: 1.26-6.44) for worsening disability. Worsening disability also impacted mortality in Time 2 (≥2.5 to <7.5 years of follow-up, HR=2.43 (95%CI: 1.03-5.73)). After Time 3 (≥7.5 years or longer), we found no evidence of an effect of age and disability on long-term mortality. Sociodemographic factors such as male gender, illiteracy, as well as some cerebrovascular risk factors as AF (HR=1.32; 95%CI: 0.97-1.80), diabetes (HR=1.51, 95%CI: 1.16-1.97), and recurrent stroke (HR=1.73; 95%CI: 1.34-2.24) had fixed high-risk effects on mortality up to 12 years of follow-up. On the other hand, the use of medications post-stroke (HR=0.53; 95%CI: 0.34-0.84) and continuous medication use (HR=0.46; 95%CI: 0.31-0.69) had a fixed protective effect on mortality (Figures 3 and 4, Tables 2 and 3).

Figure 2
Kaplan Meyer survival curve for all-cause mortality by sex among ischemic stroke patients during the 12-year follow-up in the Study of Stroke Mortality and Morbidity (EMMA) cohort.
Figure 3
Fixed baseline risk factors without time-dependent effect on all-cause mortality in the Study of Stroke Mortality and Morbidity (EMMA) cohort.
Figure 4
Fixed baseline risk factors with time-dependent effect on all-cause mortality in the Study of Stroke Mortality and Morbidity (EMMA) cohort.
Table 2
Hazard ratios (HR) of all-cause mortality for fixed baseline risk factors without time-dependent effect along 12-y follow-up in the Study of Stroke Mortality and Morbidity (EMMA) cohort.
Table 3
Hazard ratios (HR) of all-cause mortality for fixed baseline risk factors with time-dependent effect along 12-y follow-up in the Study of Stroke Mortality and Morbidity (EMMA) cohort.

Discussion

This study aimed to investigate the time-dependent effect of main traditional risk factors on all-cause mortality during 12 years of follow-up. Our main findings revealed the different impacts of those risk factors on all-cause mortality over time. Male gender, illiteracy, previous history of stroke (recurrent stroke at baseline), and diabetes were associated with increased, but fixed hazard ratios of long-term mortality. On the other hand, the use of medication to control traditional CVRF (post-stroke use and continuous use) caused a reduction in the risk of mortality by 50% up to 12 years of follow-up. Age and post-stroke disability were the only risk factors with time-dependent effects on mortality mainly between 6 months and two and a half years after stroke. Compared to those without disability, individuals who maintained or worsened their disability 6 months after the index event had almost three times the risk of death in the first two years after stroke.

Previously published articles, including data reported from the EMMA cohort (22. Goulart AC, Olmos RD, Santos IS, Tunes G, Alencar AP, Thomas N, et al. The impact of atrial fibrillation and long-term oral anticoagulant use on all-cause and cardiovascular mortality: a 12-year evaluation of the prospective Brazilian Study of Stroke Mortality and Morbidity. Int J Stroke 2022; 17: 48-58, doi: 10.1177/1747493021995592.
https://doi.org/10.1177/1747493021995592...
,44. Castro HHG, Alencar AP, Benseãor IM, Lotufo PA, Goulart AC. Multimorbidities Are associated to lower survival in ischaemic stroke: results from a Brazilian Stroke Cohort (EMMA Study). Cerebrovasc Dis 2017; 44: 232-239, doi: 10.1159/000479827.
https://doi.org/10.1159/000479827...
-5. Goulart AC, Fernandes TG, Santos IS, Alencar AP, Benseãor IM, Lotufo PA. Predictors of long-term survival among first-ever ischemic and hemorrhagic stroke in a Brazilian stroke cohort. BMC Neurol 2013; 13: 51, doi: 10.1186/1471-2377-13-51.
https://doi.org/10.1186/1471-2377-13-51...
66. de Mello RF, Santos IS, Alencar AP, Benseãor IM, Lotufo PA, Goulart AC. Major depression as a predictor of poor long-term survival in a Brazilian Stroke Cohort (Study of Stroke Mortality and Morbidity in Adults) EMMA study. J Stroke Cerebrovasc Dis 2016; 25: 618-625, doi: 10.1016/j.jstrokecerebrovasdis.2015.11.021.
https://doi.org/10.1016/j.jstrokecerebro...
), showed increased risk of mortality among those subjects with low education, cardiovascular risk factors, and post-stroke depression. Even though those previous studies did not account for risk factors' time-dependent effects, our study supports some of these findings. Here, low education, diabetes, male sex, and recurrent stroke were found to have the same positive association with mortality during the entire follow-up period (12 years).

Regarding functional disability, studies have found a significant association between low level of education and high level of disability, as well as post-stroke disability and mortality (77. Hawkes MA, Blaginykh E, Ruff MW, Burrus T, Wijdicks EFM, Rabinstein AA. Long-term mortality, disability and stroke recurrence in patients with basilar artery occlusion. Eur J Neurol 2020; 27: 579-585, doi: 10.1111/ene.14126.
https://doi.org/10.1111/ene.14126...
,1616. Fernandes TG, Goulart AC, Santos-Junior WR, Alencar AP, Benseãor IM, Lotufo PA. Educational levels and the functional dependence of ischemic stroke survivors. Cad Saude Publica 2012 ; 28: 1581-1590, doi: 10.1590/S0102-311X2012000800016.
https://doi.org/10.1590/S0102-311X201200...
). Our results expand these findings by showing the influence of time on the effect of disability on mortality, regardless of educational attainment. We found a stronger association between functional disability and mortality shortly after stroke onset (from 6 months up to 2.5 years).

The dependent effect of age in our study was evaluated for a 1-year difference in age at the index event, and the increased HR in the risk of dying was mainly verified in the first two years after stroke (Time 1). In previous analyses from the EMMA cohort, where follow-up data was analyzed in years 1, 2, 3, and 4 after stroke, age was a strong predictor of mortality, with HR around 3 times among those aged 80 years or more (55. Goulart AC, Fernandes TG, Santos IS, Alencar AP, Benseãor IM, Lotufo PA. Predictors of long-term survival among first-ever ischemic and hemorrhagic stroke in a Brazilian stroke cohort. BMC Neurol 2013; 13: 51, doi: 10.1186/1471-2377-13-51.
https://doi.org/10.1186/1471-2377-13-51...
). In the latter method, time of follow-up was stratified, but the changes in the dependent effect of age on long-term mortality were not considered. In both analyses (previous and current) mortality was calculated by Cox regression models, but in the present analyses, the measurement of the time-dependent effect was more accurate in quantifying the influence of each risk factor over time than simply considering different follow-up periods.

Moreover, studies have found medication use for traditional CVRF to have an important protective effect against mortality in stroke patients (1717. Hassan Y, Al-Jabi SW, Abd Aziz N, Looi I, Sa'ed HZ. Effect of prestroke use of angiotensin-converting enzyme inhibitors alone versus combination with antiplatelets and statin on ischemic stroke outcome. Clin Neuropharmacol 2011; 34: 234-240, doi: 10.1097/WNF.0b013e3182348abe.
https://doi.org/10.1097/WNF.0b013e318234...
-18. van Dongen MM, Aarnio K, Martinez-Majander N, Pirinen J, Sinisalo J, Lehto M, et al. Use of antihypertensive medication after ischemic stroke in young adults and its association with long-term outcome. Ann Med 2019; 51: 68-77, doi: 10.1080/07853890.2018.1564358.
https://doi.org/10.1080/07853890.2018.15...
1919. van Dongen MM, Aarnio K, Martinez-Majander N, Pirinen J, Sinisalo J, Lehto M, et al. Use of statins after ischemic stroke in young adults and its association with long-term outcome. Stroke 2019; 50: 3385-3392, doi: 10.1161/STROKEAHA.119.026992.
https://doi.org/10.1161/STROKEAHA.119.02...
). This study also supports those findings. When analyzing medication use after stroke and continuous medication use (including pre-stroke), we found a strong inverse association with all-cause mortality after 12 years. Medication use had a fixed and constant effect during the entire follow-up period.

To the best of our knowledge, no prospective stroke cohorts have investigated the time-dependent effects of baseline risk factors on long-term mortality. In a different scenario, Bellera et al. evaluated the varying effects of risk factors on mortality among women with breast cancer, showing that traditional Cox regression might mask the effects of certain variables in the short-term versus long-term follow-up (2020. Bellera CA, MacGrogan G, Debled M, de Lara CT, Brouste V, Mathoulin-Pélissier S. Variables with time-varying effects and the Cox model: some statistical concepts illustrated with a prognostic factor study in breast cancer. BMC Med Res Methodol 2010; 10: 20, doi: 10.1186/1471-2288-10-20.
https://doi.org/10.1186/1471-2288-10-20...
). They found that a particular risk factor increased the risk of metastases early on, becoming protective afterward. Their findings make evident the importance of accounting for the time-dependent effects of baseline variables, especially as the cohort ages during follow-up years. Other previous studies with stroke populations used different methodology to evaluate the time-varying risk of some selected risk factors (2121. Huang HC, Chang CH, Hu CJ, Shyu ML, Chen CI, Huang CS, et al. Time-varying effects of psychological distress on the functional recovery of stroke patients. Arch Phys Med Rehabil 2017; 98: 722-729, doi: 10.1016/j.apmr.2016.09.120.
https://doi.org/10.1016/j.apmr.2016.09.1...
,2222. Liu KF, Lin HR, Lee TY, Lin KC. Time-varying risk factors associated with the progress of functional recovery and psychological distress in first-ever stroke patients. J Neurosci Nurs 2022; 54: 80-85, doi: 10.1097/JNN.0000000000000631.
https://doi.org/10.1097/JNN.000000000000...
). In the study of Huang et al. (2121. Huang HC, Chang CH, Hu CJ, Shyu ML, Chen CI, Huang CS, et al. Time-varying effects of psychological distress on the functional recovery of stroke patients. Arch Phys Med Rehabil 2017; 98: 722-729, doi: 10.1016/j.apmr.2016.09.120.
https://doi.org/10.1016/j.apmr.2016.09.1...
), the time-varying effects of psychological distress on the functional recovery of stroke patients were analyzed. Participants were interviewed five times after stroke onset (from 5 days up to 6 months). They reported a significant effect of psychological distress on functional recovery over time but did not analyze its effects on mortality.

In the study of Liu et al. (2222. Liu KF, Lin HR, Lee TY, Lin KC. Time-varying risk factors associated with the progress of functional recovery and psychological distress in first-ever stroke patients. J Neurosci Nurs 2022; 54: 80-85, doi: 10.1097/JNN.0000000000000631.
https://doi.org/10.1097/JNN.000000000000...
), risk factors associated with functional recovery and psychological distress in stroke patients were evaluated using 3 different measurements in time, and it was demonstrated that the National Institutes of Health Stroke Scale scores (NIHSS) and depression levels affected functional recovery. Also, depression was associated with anxiety, low social support levels, and low education. The main conclusions of this study revealed the importance of early poststroke rehabilitation with 50% of patients with full functional recovery after 6 months, but the effect of time-varying risk factors on mortality was not included in the analyses of Liu et al. (2222. Liu KF, Lin HR, Lee TY, Lin KC. Time-varying risk factors associated with the progress of functional recovery and psychological distress in first-ever stroke patients. J Neurosci Nurs 2022; 54: 80-85, doi: 10.1097/JNN.0000000000000631.
https://doi.org/10.1097/JNN.000000000000...
). In our study, although we evaluated mostly variables measured at baseline, we also found that functional disability should be addressed as early as possible due to its stronger association with mortality in the first 2 years.

Strengths and limitations

The current research corroborates the knowledge gap about the different effects that risk factors can have over time, even when measured only at baseline. Of note, the classical Cox model relies on the proportional hazard's assumption, implying that the factors have a constant impact on the hazard - or risk - over time (1414. Cox DR. Regression model and life tables. J Royal Stat Soc 1972; 34: 187-220.). The present study had a robust validation process of stroke cases, including a collection of data on all cerebrovascular risk factors and medications (before and during follow-up after an index event), entirely supervised by the medical research team.

This study had some limitations, most of them related to selection bias. The EMMA cohort is based on data from a single center that is a secondary community hospital without specialized treatment options (e.g., thrombolysis/thrombectomy). Although patients with IS who report any acute therapy are usually transferred to our reference (tertiary) hospital, we cannot rule out a selection bias from a highly selected segment of stroke patients excluded from the EMMA cohort that limits the external generalization of our data.

Although we did not have NIHSS or different measures of disability at hospital admission and even stroke recurrence data during follow-up, we measured mRS at two points during follow-up (1 and 6 months after stroke). This information about disability after stroke allowed us to calculate the evolution of functionality (unchanged, improved, and worsened) to better evaluate the time-dependent effect of this relevant factor associated with both short- and long-term mortality.

Although most individuals had mild or no disability, which could influence the magnitude of our findings by underestimating the impact of functional disability on long-term mortality, we were able to report important data about higher levels of functional disability at 6 months related to higher mortality during follow-up. This data is particularly important for communities with low-middle socioeconomic status and limited resources for rehabilitation after stroke (2323. GBD 2016 Stroke Collaborators. Global, regional, and national burden of stroke, 1990-2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet Neurol 2019; 18: 439-458, doi: 10.1016/S1474-4422(19)30034-1.
https://doi.org/10.1016/S1474-4422(19)30...
).

In addition, most individuals who are disabled after stroke tend to have more difficulty of locomotion to receive proper care and achieve a good recovery. In the present analysis, we found about 8.3% of functional deterioration six months after stroke. We speculated that the lack of rehabilitation after hospital discharge could have contributed to this worsening. About 45.7% of participants reported some kind of rehabilitation (physiotherapy and speech therapy) after their stroke event, however we did not know the reasons for the low rehabilitation adherence in the EMMA cohort.

Finally, an observational study design does not allow us to rule out the possibility of immortal time bias in our analyses (2424. Suissa S. Immortal time bias in pharmaco-epidemiology. Am J Epidemiol 2008; 167: 492-499, doi: 10.1093/aje/kwm324.
https://doi.org/10.1093/aje/kwm324...
,2525. Lévesque LE, Hanley JA, Kezouh A, Suissa S. Problem of immortal time bias in cohort studies: example using statins for preventing progression of diabetes. BMJ 2010; 12: b5087, doi: 10.1136/bmj.b5087.
https://doi.org/10.1136/bmj.b5087...
).

Conclusions

Functional disability was the main contributor to short- and long-term post-stroke mortality. Male sex, low education attainment, and traditional baseline CVRF such as diabetes, AF, and recurrent stroke had a fixed high-risk effect on mortality up to 12 years after the acute event. Medication use, particularly the continuous use of medications for CVRF, more than halved the risk of death throughout the 12-year follow-up period. Thus, the present study reinforces tertiary prevention by early rehabilitation among those with post-stroke functional disability and secondary prevention by the use of medications for controlling risk factors and thus prevent recurrent strokes and reduce mortality risk after stroke.

Supplementary Material

Click here to view [pdf].

Acknowledgments

This study was funded by the Research Support Foundation of the State of São Paulo (FAPESP, grant No. 2018/05512-3) and by the NIHR Global Health Research Group on AF (grant No. 17/63/121).

References

  • 1
    Dekker FW, de Mutsert R, van Dijk PC, Zoccali C, Jager KJ. Survival analysis: time-dependent effects and time-varying risk factors. Kidney Int 2008; 74: 994-997, doi: 10.1038/ki.2008.328.
    » https://doi.org/10.1038/ki.2008.328
  • 2
    Goulart AC, Olmos RD, Santos IS, Tunes G, Alencar AP, Thomas N, et al. The impact of atrial fibrillation and long-term oral anticoagulant use on all-cause and cardiovascular mortality: a 12-year evaluation of the prospective Brazilian Study of Stroke Mortality and Morbidity. Int J Stroke 2022; 17: 48-58, doi: 10.1177/1747493021995592.
    » https://doi.org/10.1177/1747493021995592
  • 3
    Tang Z, Zhou T, Luo Y, Xie C, Huo D, Tao L, et al. Risk factors for cerebrovascular disease mortality among the elderly in Beijing: a competing risk analysis. PLoS One 2014 ; 9: e87884, doi: 10.1371/journal.pone.0087884.
    » https://doi.org/10.1371/journal.pone.0087884
  • 4
    Castro HHG, Alencar AP, Benseãor IM, Lotufo PA, Goulart AC. Multimorbidities Are associated to lower survival in ischaemic stroke: results from a Brazilian Stroke Cohort (EMMA Study). Cerebrovasc Dis 2017; 44: 232-239, doi: 10.1159/000479827.
    » https://doi.org/10.1159/000479827
  • 5
    Goulart AC, Fernandes TG, Santos IS, Alencar AP, Benseãor IM, Lotufo PA. Predictors of long-term survival among first-ever ischemic and hemorrhagic stroke in a Brazilian stroke cohort. BMC Neurol 2013; 13: 51, doi: 10.1186/1471-2377-13-51.
    » https://doi.org/10.1186/1471-2377-13-51
  • 6
    de Mello RF, Santos IS, Alencar AP, Benseãor IM, Lotufo PA, Goulart AC. Major depression as a predictor of poor long-term survival in a Brazilian Stroke Cohort (Study of Stroke Mortality and Morbidity in Adults) EMMA study. J Stroke Cerebrovasc Dis 2016; 25: 618-625, doi: 10.1016/j.jstrokecerebrovasdis.2015.11.021.
    » https://doi.org/10.1016/j.jstrokecerebrovasdis.2015.11.021
  • 7
    Hawkes MA, Blaginykh E, Ruff MW, Burrus T, Wijdicks EFM, Rabinstein AA. Long-term mortality, disability and stroke recurrence in patients with basilar artery occlusion. Eur J Neurol 2020; 27: 579-585, doi: 10.1111/ene.14126.
    » https://doi.org/10.1111/ene.14126
  • 8
    Goulart AC, Bustos IR, Abe IM, Pereira AC, Fedeli LM, Benseãor IM, et al. A stepwise approach to stroke surveillance in Brazil: the EMMA (Estudo de Mortalidade e Morbidade do Acidente Vascular Cerebral) study. Int J Stroke 2010; 5: 284-289, doi: 10.1111/j.1747-4949.2010.00441.x.
    » https://doi.org/10.1111/j.1747-4949.2010.00441.x
  • 9
    World Health Organization. WHO STEPS stroke manual: the WHO STEPwise approach to stroke surveillance/noncommunicable diseases and mental health. World Health Organization; 2006. [cited 12/06/2010]. Available from: World Health Organization http://www.who.int/chp/steps/stroke/manual/en/index.html
    » http://www.who.int/chp/steps/stroke/manual/en/index.html
  • 10
    Lotufo PA, Benseãor IM. Improving WHO STEPS Stroke in Brazil. Lancet Neurol 2007; 6: 387-388, doi: 10.1016/S1474-4422(07)70091-1.
    » https://doi.org/10.1016/S1474-4422(07)70091-1
  • 11
    Bamford JM, Sandercock PA, Warlow CP and Slattery J. Interobserver agreement for the assessment of handicap in stroke patients. Stroke 1989; 20: 828, doi: 10.1161/01.STR.20.6.828.
    » https://doi.org/10.1161/01.STR.20.6.828
  • 12
    Dennis M, Mead G, Doubal F, Graham C. Determining the modified Rankin Score after stroke by postal and telephone questionnaires. Stroke 2012; 43: 851-853, doi: 10.1161/STROKEAHA.111.639708.
    » https://doi.org/10.1161/STROKEAHA.111.639708
  • 13
    Kaplan EL, Meier P. Nonparametric Estimation from incomplete observations. J Am Stat Assoc 1958; 53 : 457-481, doi: 10.1080/01621459.1958.10501452.
    » https://doi.org/10.1080/01621459.1958.10501452
  • 14
    Cox DR. Regression model and life tables. J Royal Stat Soc 1972; 34: 187-220.
  • 15
    Grambsch P, Therneau T. Proportional hazards tests and diagnostics based on weighted residuals. Biometrika 1994; 81: 515-526, doi: 10.1093/biomet/81.3.515.
    » https://doi.org/10.1093/biomet/81.3.515
  • 16
    Fernandes TG, Goulart AC, Santos-Junior WR, Alencar AP, Benseãor IM, Lotufo PA. Educational levels and the functional dependence of ischemic stroke survivors. Cad Saude Publica 2012 ; 28: 1581-1590, doi: 10.1590/S0102-311X2012000800016.
    » https://doi.org/10.1590/S0102-311X2012000800016
  • 17
    Hassan Y, Al-Jabi SW, Abd Aziz N, Looi I, Sa'ed HZ. Effect of prestroke use of angiotensin-converting enzyme inhibitors alone versus combination with antiplatelets and statin on ischemic stroke outcome. Clin Neuropharmacol 2011; 34: 234-240, doi: 10.1097/WNF.0b013e3182348abe.
    » https://doi.org/10.1097/WNF.0b013e3182348abe
  • 18
    van Dongen MM, Aarnio K, Martinez-Majander N, Pirinen J, Sinisalo J, Lehto M, et al. Use of antihypertensive medication after ischemic stroke in young adults and its association with long-term outcome. Ann Med 2019; 51: 68-77, doi: 10.1080/07853890.2018.1564358.
    » https://doi.org/10.1080/07853890.2018.1564358
  • 19
    van Dongen MM, Aarnio K, Martinez-Majander N, Pirinen J, Sinisalo J, Lehto M, et al. Use of statins after ischemic stroke in young adults and its association with long-term outcome. Stroke 2019; 50: 3385-3392, doi: 10.1161/STROKEAHA.119.026992.
    » https://doi.org/10.1161/STROKEAHA.119.026992
  • 20
    Bellera CA, MacGrogan G, Debled M, de Lara CT, Brouste V, Mathoulin-Pélissier S. Variables with time-varying effects and the Cox model: some statistical concepts illustrated with a prognostic factor study in breast cancer. BMC Med Res Methodol 2010; 10: 20, doi: 10.1186/1471-2288-10-20.
    » https://doi.org/10.1186/1471-2288-10-20
  • 21
    Huang HC, Chang CH, Hu CJ, Shyu ML, Chen CI, Huang CS, et al. Time-varying effects of psychological distress on the functional recovery of stroke patients. Arch Phys Med Rehabil 2017; 98: 722-729, doi: 10.1016/j.apmr.2016.09.120.
    » https://doi.org/10.1016/j.apmr.2016.09.120
  • 22
    Liu KF, Lin HR, Lee TY, Lin KC. Time-varying risk factors associated with the progress of functional recovery and psychological distress in first-ever stroke patients. J Neurosci Nurs 2022; 54: 80-85, doi: 10.1097/JNN.0000000000000631.
    » https://doi.org/10.1097/JNN.0000000000000631
  • 23
    GBD 2016 Stroke Collaborators. Global, regional, and national burden of stroke, 1990-2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet Neurol 2019; 18: 439-458, doi: 10.1016/S1474-4422(19)30034-1.
    » https://doi.org/10.1016/S1474-4422(19)30034-1
  • 24
    Suissa S. Immortal time bias in pharmaco-epidemiology. Am J Epidemiol 2008; 167: 492-499, doi: 10.1093/aje/kwm324.
    » https://doi.org/10.1093/aje/kwm324
  • 25
    Lévesque LE, Hanley JA, Kezouh A, Suissa S. Problem of immortal time bias in cohort studies: example using statins for preventing progression of diabetes. BMJ 2010; 12: b5087, doi: 10.1136/bmj.b5087.
    » https://doi.org/10.1136/bmj.b5087

Publication Dates

  • Publication in this collection
    22 Sept 2023
  • Date of issue
    2023

History

  • Received
    09 June 2023
  • Accepted
    09 Aug 2023
Associação Brasileira de Divulgação Científica Av. Bandeirantes, 3900, 14049-900 Ribeirão Preto SP Brazil, Tel. / Fax: +55 16 3315-9120 - Ribeirão Preto - SP - Brazil
E-mail: bjournal@terra.com.br