Desenvolvimento e Validação de Modelos Preditores de Mortalidade Cardíaca e Transplante na Terapia de Ressincronização Cardíaca
Desenvolvimento e Validação de Modelos Preditores de Mortalidade Cardíaca e Transplante na Terapia de Ressincronização Cardíaca
Autores:
ARTIGO ORIGINAL
Arquivos Brasileiros de Cardiologia
versão impressa ISSN 0066-782Xversão On-line ISSN 1678-4170
Arq. Bras. Cardiol. vol.105 no.4 São Paulo out. 2015 Epub 07-Ago-2015
https://doi.org/10.5935/abc.20150093
Resumo
Fundamento:
A terapia de ressincronização cardíaca pode demonstrar resultados desfavoráveis em 30-40% dos casos.
Objetivo:
Este estudo teve por objetivo desenvolver modelos preditores para o desfecho combinado de morte cardíaca ou transplante (Tx) em diferentes estágios da terapia de ressincronização cardíaca (TRC).
Métodos:
Estudo prospectivo observacional de 116 pacientes com idade média de 64,8 ± 11,1 anos, dos quais 68,1% estavam em classe funcional (CF) III e 31,9% em classe IV ambulatorial. Variáveis clínicas, eletrocardiográficas e ecocardiográficas foram avaliadas com regressão de Cox e curvas de Kaplan-Meier.
Resultados:
O desfecho mortalidade/Tx cardíacos foi de 16,3% durante seguimento de 34,0 ± 17,9 meses. No período pré-implante, disfunção ventricular direita (DVD), fração de ejeção < 25% e uso de altas doses de diuréticos (ADD) aumentaram o risco de morte e Tx cardíacos em 3,9, 4,8 e 5,9 vezes, respectivamente. No primeiro ano após TRC, DVD, AHDD e hospitalização por insuficiência cardíaca congestiva elevaram o risco de morte (razões de risco de 3,5, 5,3 e 12,5, respectivamente). No segundo ano após TRC, DVD e CF III/IV foram fatores de risco significativos de mortalidade no modelo multivariado de Cox. As acurácias dos modelos foram 84,6% na pré-implante, 93% no primeiro ano após TRC e 90,5% no segundo ano após TRC. Os modelos foram validados por bootstrapping.
Conclusão:
Desenvolvemos modelos preditores de morte e Tx cardíacos em diferentes estágios de TRC com base na análise de variáveis clínicas e ecocardiográficas simples e de fácil obtenção. Os modelos mostraram boa acurácia e ajuste, foram validados internamente e são úteis para a seleção, o monitoramento e a orientação de pacientes indicados para TRC.
Palavras-chave: Transplante de Coração / mortalidade; Insuficiência Cardíaca; Terapia de Ressincronização Cardíaca; Seguimentos; Marca-Passo Artificial
Abstract
Background:
30-40% of cardiac resynchronization therapy cases do not achieve favorable outcomes.
Objective:
This study aimed to develop predictive models for the combined endpoint of cardiac death and transplantation (Tx) at different stages of cardiac resynchronization therapy (CRT).
Methods:
Prospective observational study of 116 patients aged 64.8 ± 11.1 years, 68.1% of whom had functional class (FC) III and 31.9% had ambulatory class IV. Clinical, electrocardiographic and echocardiographic variables were assessed by using Cox regression and Kaplan-Meier curves.
Results:
The cardiac mortality/Tx rate was 16.3% during the follow-up period of 34.0 ± 17.9 months. Prior to implantation, right ventricular dysfunction (RVD), ejection fraction < 25% and use of high doses of diuretics (HDD) increased the risk of cardiac death and Tx by 3.9-, 4.8-, and 5.9-fold, respectively. In the first year after CRT, RVD, HDD and hospitalization due to congestive heart failure increased the risk of death at hazard ratios of 3.5, 5.3, and 12.5, respectively. In the second year after CRT, RVD and FC III/IV were significant risk factors of mortality in the multivariate Cox model. The accuracy rates of the models were 84.6% at preimplantation, 93% in the first year after CRT, and 90.5% in the second year after CRT. The models were validated by bootstrapping.
Conclusion:
We developed predictive models of cardiac death and Tx at different stages of CRT based on the analysis of simple and easily obtainable clinical and echocardiographic variables. The models showed good accuracy and adjustment, were validated internally, and are useful in the selection, monitoring and counseling of patients indicated for CRT.
Keywords: Heart Transplantation / mortality; Heart Failure / physiopathology; Cardiac Resynchronization Therapy; Follow-Up Studies; Pacemaker, Artificial
Introdução
As principais diretrizes internacionais recomendam fortemente (classe I) a terapia de ressincronização cardíaca (TRC) para pacientes com insuficiência cardíaca congestiva (ICC) e classe funcional (CF) da New York Heart Association (NYHA) II ou III ou classe IV ambulatorial (pacientes ambulatoriais em uso de medicação oral) quando apresentam distúrbios de condução intraventricular e fração de ejeção (FE) ≤ 35% durante tratamento médico ótimo1.
Entretanto, 30%-40% dos casos de TRC não têm desfechos favoráveis, indicando que os pacientes são submetidos a cirurgia com altos riscos e custos, mas sem benefícios clínicos, hemodinâmicos ou de sobrevida2. Portanto, índices ou escores multifatoriais precisam ser desenvolvidos para identificar de maneira mais precisa fatores preditores de sobrevida e os respondedores ao tratamento3, 4. Esses índices devem envolver variáveis relacionadas à redução de mortalidade, com alta sensibilidade e especificidade.
Este artigo teve por objetivo desenvolver modelos preditores para o desfecho combinado de morte e transplante (Tx) cardíacos em diferentes estágios da TRC.
Métodos
Este estudo prospectivo observacional avaliou 116 pacientes submetidos a implante de marca-passo "multi-sítio" consecutivamente em um hospital universitário terciário entre janeiro de 2008 e março de 2013 (Tabela 1). Esses pacientes apresentavam CF da NYHA III ou IV ambulatorial, FE ≤ 35%, QRS ≥ 120 ms (bloqueio de ramo esquerdo [BRE] e bloqueio de ramo direito [BRD] com hemibloqueio ou ritmo de marca-passo) e tratamento otimizado. Os critérios de exclusão foram comorbidades graves, indicação prévia de implante de marca-passo, hospitalização por insuficiência cardíaca CF IV (NYHA), valvopatia primária e dados incompletos.
Tabela 1 Características basais e comparação dos resultados de algumas variáveis durante os períodos avaliados
| Variáveis | Tempo 1 | Tempo 2 | Valor de p | Tempo 3 | Valor de p* |
|---|---|---|---|---|---|
| Pacientes | 116 | 114 | - | 92 | - |
| Idade (anos) | 64,8 ± 11,1 | - | - | - | - |
| Sexo (masculino) | 69,8% | - | - | - | - |
| IMC | 25,8 ± 4,1 | - | - | - | - |
| Betabloqueadores | 88,7% | 89,2% | - | 91,8% | |
| Inibidores da ECA | 97,4% | 96% | 95,9% | ||
| Furosemida ≥ 80mg/dia | 31,9% | 17% | - | - | - |
| Miocardiopatia dilatada | 59,4% | - | - | - | - |
| Miocardiopatia isquêmica | 29,3% | - | - | - | - |
| Doença de Chagas | 11,2% | - | - | - | - |
| Fibrilação atrial | 12% | - | - | - | - |
| TRC-D | 54,2% | - | - | - | - |
| BRE | 71,55% | - | - | - | - |
| BRD com hemibloqueio | 12% | ||||
| Marca-passo | 16,3% | ||||
| Veia posterolateral | 45,4% | ||||
| Veias anterolaterais | 52,5% | ||||
| Largura prévia do QRS (ms) | 160 | 140 | < 0,001 | - | - |
| Número de hospitalizações por ICC | 108 | 24 | < 0,001 | 16 | 0,79* |
| FE (mediana) | 29% | 33% | < 0,001 | 35% | 0,03* |
| DDVE (mm) | 70 | 66 | < 0,001 | 65 | 0,73* |
| PA sistólica (mmHg) | 115 | 119,6 | < 0,001 | 121,8 | 0,84* |
| PA diastólica (mmHg) | 70 | 80 | 0,07 | 70 | 0,34* |
| CF IIl (NYHA) | 68,1% | 8,7% | < 0,001 | 12% | 0,07* |
| CF IV (NYHA) | 31,9% | 6,1% | < 0,001 | 7,6% | 0,07* |
| DD | < 0,001 | - | 0,06* | ||
| DD Grau I | 34,6% | 59,2% | - | 63,2% | - |
| DD Grau II | 23,7% | 27,1% | - | 13,9% | - |
| DD Grau III | 29,7% | 8,7% | - | 16,4% | - |
| DD Grau IV | 11,8% | 4,8% | - | 5,0% | - |
| RM | - | - | 0,008 | - | 0,009* |
| Sem RM | 3,4% | 5,3% | - | 15,3% | |
| Leve RM | 50,4% | 66,0% | - | 56,0% | - |
| Moderada RM | 30,4% | 18,7% | - | 18,6% | - |
| Grave RM | 15,6% | 9,8% | - | 9,8% | - |
| Disfunção de VD | 20,9% | 17% | 0,62 | 12% | 0,5* |
| Creatinina (mg/dl) | 1,1 | 1,1 | - | 1,2 | - |
Tempo 1, pré-implante; tempo 2, 1 ano; tempo 3, 2 anos.
*Análise do tempo 3 em relação ao tempo 2; largura do QRS, fração de ejeção, diâmetro diastólico do ventrículo esquerdo e pressão arterial foram variáveis sem distribuição normal (valores de mediana); IMC: índice de massa corporal; ECA: enzima conversora da angiotensina; TRC-D: cardioversor-desfibrilador com marca-passo biventricular; BRE: bloqueio de ramo esquerdo; BRD: bloqueio de ramo direito; ICC: insuficiência cardíaca congestiva; FE: fração de ejeção; DDVE: diâmetro diastólico do ventrículo esquerdo; PA: pressão arterial; CF: classe funcional (NYHA); DD: disfunção diastólica; RM: regurgitação mitral; VD: ventrículo direito.
Dos 147 pacientes submetidos à implante durante o período do estudo, apenas 116 foram incluídos pelas seguintes razões: 4 tinham FE > 35%; 3 tinham bloqueio atrioventricular total; 2 tinham valvopatia cardíaca primária; 2 apresentaram infecção do marca-passo; 7 tinham dados incompletos; 4 apresentaram perda de captura no eletrodo do ventrículo esquerdo; 2 não completaram o seguimento; 1 apresentou comorbidade grave; 5 foram hospitalizados por ICC classe IV no momento da inclusão; e 1 morreu de infecção respiratória prematura.
Os eletrodos do ventrículo direito foram posicionados preferencialmente na região apical (84%). Os modelos St. Jude Medical, Biotronik, Medtronic e Guidant foram usados em 92, 12, 10 e 2 pacientes, respectivamente. Pacientes com indicação para implante concomitante de cardioversor-desfibrilador (grupo TRC-D) (54% dos 116 pacientes) também foram incluídos neste estudo. Tal indicação foi para prevenção primária em 47 pacientes e para prevenção secundária em 16 pacientes.
As avaliações foram realizadas no período pré-implante (primeira análise), um ano após a implante (segunda análise) e dois anos após a implante (terceira análise) conforme protocolo fixo. Foram analisadas 12 variáveis clínicas, 8 eletrocardiográficas e 7 ecocardiográficas. As variáveis clínicas foram: idade; sexo; índice de massa corporal; caquexia cardíaca; CF; etiologia da miocardiopatia; veia cardíaca em que o eletrodo foi posicionado no ventrículo esquerdo; nível sérico de creatinina; pressão arterial sistólica e diastólica; uso de alta dose de diuréticos de alça (≥ 80 mg/dia de furosemida); e hospitalização por insuficiência cardíaca. As variáveis eletrocardiográficas foram: fibrilação atrial; BRE ou BRD; marca-passo cardíaco prévio; bloqueio atrioventricular de 1º grau; duração do QRS; estreitamento do QRS após implante; onda R na derivação V1 em pacientes com BRE; e eixo QRS no plano frontal após implante. As variáveis ecocardiográficas foram: diâmetro diastólico e sistólico ventricular esquerdo (VE); FE calculada usando o método de Simpson; grau de disfunção diastólica (DD) de I a IV; grau de regurgitação mitral de I a III; disfunção ventricular direita (DVD); e dissincronismo.
Registrou-se eletrocardiograma de superfície de 12 derivações na velocidade de 25 mm/s e amplitude de 10 mm/mV. Considerou-se a maior duração do QRS medida em uma das derivações no plano frontal ou horizontal, que foi a derivação com o mais alto valor, que permitia melhor avaliação. Definiu-se mortalidade cardíaca como morte por ICC em estágio avançado ou morte súbita.
Parâmetros ecocardiográficos
As diretrizes ecocardiográficas para análise de vários parâmetros ecocardiográficos, assim como as diretrizes de dissincronismo para análise de tais parâmetros, foram seguidas5, 6. Três médicos experientes realizaram os exames ecocardiográficos, dos quais, 72% foram feitos pelo mesmo especialista. O sistema GE Vivid 7 Ultrasound System (GE Healthcare, Fairfield, CT, EUA) foi utilizado para os exames.
A função sistólica das câmaras cardíacas foi analisada usando-se o método de Simpson no modo bidimensional. Os diâmetros ventriculares foram obtidos na ecocardiografia modo M, de acordo com a diretriz padrão5. A função ventricular direita foi analisada qualitativamente, diferenciada em presença e ausência de qualquer grau de disfunção5.
A função diastólica foi analisada avaliando-se o fluxo mitral (em repouso e após a manobra de Valsalva), imagens de Doppler tecidual e velocidade de propagação de fluxo com modo M colorido. Usaram-se os resultados para classificar a DD em quatro graus (0, ausente; I, leve; II, moderada; III, acentuada ou com disfunção restritiva; e IV, grave ou com disfunção restritiva irreversível)7.
Avaliou-se o grau de regurgitação mitral como a percentagem de enchimento do átrio esquerdo usando-se ecocardiografia com Doppler colorido. A percentagem foi inferior a 20% no refluxo leve, entre 20% e 40% no moderado e valores acima desses indicaram refluxo grave5. Nesse contexto prático, o efeito Coanda foi interpretado como refluxo moderado quando restrito à parede lateral do átrio, e acentuado quando se estendeu pelo polo superior do átrio esquerdo.
Todos os pacientes assinaram o termo de consentimento livre e informado, tendo o Comitê de Ética do hospital aprovado o estudo, cujo protocolo obedeceu às normas éticas da declaração de Helsinki.
Análise estatística
As variáveis categóricas foram apresentadas como frequências e percentagens, enquanto as contínuas foram apresentadas como média e desvio padrão, ou mediana. As variáveis categóricas foram comparadas usando-se os testes de McNemar, Stuart-Maxwell ou qui-quadrado. O teste t de Student foi usado para comparar as variáveis contínuas com distribuição normal, sendo o teste Wilcoxon/Mann-Whitney U usado para as contínuas sem distribuição normal. As distribuições foram consideradas significativamente diferentes para p < 0,05.
A relação univariada das variáveis clínicas, eletrocardiográficas e ecocardiográficas com mortalidade ou Tx cardíacos foi avaliada usando-se a curva de sobrevida de Kaplan-Meier, o teste log-rank e a análise de regressão de Cox. Algumas variáveis contínuas foram avaliadas para determinar o ponto de corte.
Os modelos de regressão múltipla de Cox foram desenvolvidos nos seguintes tempos para avaliar a contribuição independente de cada uma das variáveis significativas no modelo univariado de Cox: pré-implante (tempo 1); primeiro ano após a TRC (tempo 2); e segundo ano após a TRC (tempo 3). As variáveis com p < 10% foram consideradas potenciais confundidoras. Cada variável foi incluída no modelo multivariado de acordo com a ordem descendente de risco, sendo excluída quando p ≥ 5%. Após a obtenção do modelo final, as variáveis previamente excluídas foram de novo incluídas no modelo e testadas individualmente, usando-se o mesmo critério.
Conduziram-se as análises de regressão logística usando-se a razão de risco8 como variável independente e morte/Tx cardíacos como variável dependente. Testou-se a acurácia dos modelos com a curva receiver-operating characteristic (ROC), juntamente com sensibilidade e especificidade. Modelos foram preparados dividindo-se os escores de risco em categorias de acordo com o número de variáveis presentes, sendo classificados como de risco baixo (classe A), médio (classe B) e alto (classe C).
As curvas de sobrevida de Kaplan-Meier foram elaboradas individualmente para variáveis independentes e classes de risco, sendo comparadas usando-se o teste log-rank.
Para os modelos propostos, todas as variáveis foram testadas quanto à observância à suposição de riscos proporcionais com o teste de Schoenfeld e a análise visual dos resíduos contra o tempo de mortes ou censura. Para cada modelo, analisou-se o efeito de cada observação nos parâmetros estimados. Para tal, após a deleção de uma observação, o modelo foi estimado de novo, sendo as novas estimativas comparadas às anteriores. Os valores não devem mudar muito ou o modelo é demasiado sensível para determinada observação.
Para a obtenção dos intervalos de confiança bootstrap, os dados originais foram experimentados 10 mil vezes para a obtenção de 10 mil pseudo-amostras de tamanho 60. Assim, para cada pseudo-amostra, estimaram-se razões de risco dos três modelos. Tais razões de risco estimadas foram ordenadas, relatando-se o intervalo de confiança de 95%.
Os dados foram analisados usando-se Stata/SE, versão 12.1 (StataCorp LP, College Station, TX, EUA) e o programa "R" (2014 - "R": A language and environment for statistical computing. R Foundation for Statistical Computing, Viena, Áustria).
Resultados
Durante o estudo, registraram-se 29 mortes, correspondendo a uma taxa de mortalidade total de 25% durante o seguimento de 34,09 ± 17,9 meses. Mortalidade/Tx cardíacos corresponderam a 16,3% (19 pacientes). Seis pacientes foram submetidos a Tx durante o período do estudo, cinco por ICC refratária e um por tempestade arrítmica recorrente. Três pacientes de Tx morreram prematuramente devido à gravidade da doença por ocasião do Tx. Não houve morte súbita no grupo TRC-D, mas três mortes súbitas ocorreram no grupo TRC-P (marca-passo sem desfibrilador). No grupo TRC-D, seis pacientes com taquicardia ventricular rápida ou fibrilação foram tratados com choque eficaz. A Tabela 1 mostra as características basais dos pacientes e a comparação dos resultados durante o período de avaliação.
Não houve significativa evidência estatística de que a suposição de riscos proporcionais tenha sido violada. Analisou-se o efeito de cada observação nos parâmetros estimados para cada modelo. Os dados obtidos não sugerem observações influentes. Intervalos de confiança bootstrap para o nível de significância de 95% foram obtidos e confirmaram o significado estatístico das razões de risco. Tais resultados não rejeitaram o ajuste do modelo com as variáveis propostas (Tabela 2).
Tabela 2 Intervalo de confiança bootstrap de 95% e teste formal para suposição de riscos proporcionais
| Modelo 1 | |||
|---|---|---|---|
| Covariável | IC | ||
| FE | (1.7142; 14.2053) | ||
| DVD | (1.8754; 16.6939) | ||
| ADD | (2.2563; 18.7021) | ||
| Modelo 2 | |||
| ICC | (2.1747; 11.4814) | ||
| DVD | (3.0642; 10.6684) | ||
| ADD | (4.0963; 18.3712) | ||
| Modelo 3 | |||
| CF | (3.5177; 37.5661) | ||
| DVD | (6.0592; 46.8405) | ||
| Modelo 1 | |||
| Covariável | ρ | χ2 | valor de p |
| FE | 0,073 | 0,08099 | 0,776 |
| DVD | 0,124 | 0,28512 | 0,593 |
| ADD | -0,012 | 0,00259 | 0,939 |
| Global | 0,33714 | 0,953 | |
| Modelo 2 | |||
| ICC | 0,3713 | 1,785 | 0,182 |
| DVD | 0,1089 | 0,223 | 0,637 |
| ADD | -0,0934 | 0,167 | 0,683 |
| Global | 1,905 | 0,592 | |
| Modelo 3 | |||
| CF | -0,110 | 0,118 | 0,732 |
| DVD | 0,125 | 0,162 | 0,687 |
| Global | 0,254 | 0,881 | |
IC: intervalo de confiança; FE: fração de ejeção; DVD: disfunção ventricular direita; ICC: insuficiência cardíaca congestiva; ADD: altas doses de diurético (furosemida ≥ 80 mg/dia); CF: classe funcional (NYHA) III/IV comparada com I/II.
Análise das variáveis no tempo 1 (pré-implante)
Das 27 variáveis analisadas durante o primeiro período do estudo (pré-implante), 13 foram significativas no modelo de regressão univariada de Cox. Nesse modelo, RVD, FE < 25% e o uso de altas doses de diuréticos associaram-se independentemente com aumento do desfecho mortalidade/Tx cardíacos, com razões de risco de 3,9, 4,8 e 5,9, respectivamente (Tabela 3).
Tabela 3 Análise pelo modelo de Cox relacionada com mortalidade/Tx cardíacos no tempo 1 (pré-implante)
| Variável | HR | IC 95% | p | HR | IC 95% | p |
|---|---|---|---|---|---|---|
| Univariada | Multivariada | |||||
| Hospitalização ≥ 1 | 9,23 | 1,23-69,21 | 0,031 | |||
| Disfunção de VD | 5,01 | 1,97-12,76 | 0,001 | 3.95 | 1.45-10.74 | 0.007 |
| CF III / IV | 4,87 | 1,85-12,83 | 0,001 | |||
| Doença de Chagas | 4,73 | 1,77-12,63 | 0,002 | |||
| FE < 25 % | 4,43 | 1,77-11,05 | 0,001 | 4,85 | 1,71-13,73 | 0,003 |
| Diurético ↑ | 3,89 | 1,56-9,72 | 0,004 | 5,97 | 2,15-16,53 | 0,001 |
| PAS < 100 mmHg | 3,38 | 1,35-8,46 | 0,009 | |||
| Creatinina > 1,1 | 2,85 | 1,06-7,67 | 0,038 | |||
| DDVE > 80 mm | 2,68 | 1,00-7,15 | 0,048 | |||
| PAD < 60 mmHg | 2,63 | 1,02-6,75 | 0,044 | |||
| Inibidores da ECA | 4,34 | 0,98-19,17 | 0,052 | |||
| RM grau II | 2,50 | 0,89-7,41 | 0,08 | |||
| RM grau III | 2,80 | 0,87-9,43 | 0,08 | |||
HR: hazard ratio no modelo de Cox; IC: intervalo de confiança; P: nível de significado estatístico; Diurético ↑: ≥ 80 mg de furosemida; PAS: pressão arterial sistólica; PAD: pressão arterial diastólica; CF III / IV: percentagem de classe funcional (CF) III sobre a CFC IV; Hospitalização ≥ 1: uma ou mais hospitalizações por insuficiência cardíaca congestiva; VD: ventrículo direito; FE: fração de ejeção; DDVE: diâmetro diastólico do ventrículo esquerdo; ECA: enzima conversora da angiotensina; RM: regurgitação mitral.
As variáveis significativas no modelo multivariado também foram significativas no modelo de Kaplan-Meier ao serem comparadas usando-se o teste log-rank. A análise do modelo usando a curva ROC mostrou uma área sob a curva (AUC) de 0,81, com sensibilidade de 61,1%, especificidade de 89,5% e acurácia de 84,6% (Figura 1).
Figura 1 Curva de sobrevida de Kaplan-Meier das variáveis com valor independente na análise multivariada de Cox, comparada usando o teste log-rank, com a construção da curva ROC, área sob a curva (AUC) de 0,81, sensibilidade de 61%, especificidade de 89% e acurácia de 84%. À direita, parte inferior, modelo de risco, correspondendo a baixo risco de mortalidade/Tx cardíacos a ausência das três variáveis, furosemida > 80 mg/dia, disfunção ventricular direita (VD) e fração de ejeção (FE) < 25%, ou a presença de uma delas.
A partir de combinações dessas variáveis, desenvolvemos um modelo com três classes como se segue: classe A (baixo risco para morte/Tx cardíacos), ausência de variáveis ou presença de apenas uma das variáveis significativas na análise multivariada, correspondendo a uma taxa de ausência de evento cardíaco em 30 meses de 93%. A combinação de duas (classe B) e três variáveis (classe C) resultou em uma taxa de ausência de evento cardíaco em 30 meses de 61% e 0%, respectivamente.
Análise das variáveis no tempo 2 (primeiro ano após TRC)
No tempo 2 (primeiro ano após TRC), 13 variáveis mostraram-se significativas no modelo de regressão univariada de Cox, onde RVD, uso de altas doses de diuréticos e hospitalização por ICC relacionaram-se de maneira independente com o aumento do desfecho mortalidade/Tx cardíacos, com razões de risco de 3,5, 5,3 e 12,5, respectivamente.
As variáveis significativas no modelo multivariado também foram significativas no modelo de Kaplan-Meier ao serem comparadas usando-se o teste log-rank. A análise do modelo usando a curva ROC mostrou uma AUC de 0,910, com sensibilidade de 76,4%, especificidade de 96,3% e acurácia de 93% (Figura 2).
Figura 2 Curva de sobrevida de Kaplan-Meier da variável hospitalização por insuficiência cardíaca congestiva (ICC), que, associada com disfunção ventricular direita e uso de altas doses de diuréticos, formaram o modelo preditivo de morte/Tx cardíacos no tempo 2 (primeiro ano). A ausência das três variáveis ou a presença de apenas uma (baixo risco) indica uma taxa livre de eventos em 30 meses de 98%. Na parte superior à direita, curva ROC com área sob a curva (AUC) de 0,91, sensibilidade de 76,4%, especificidade de 96,3% e acurácia de 93%.
A partir de combinações dessas variáveis, desenvolvemos um modelo com três classes (Tabela 4): classe A (baixo risco de morte/Tx cardíacos), ausência de variáveis ou presença de apenas uma das variáveis significativas na análise multivariada, correspondendo a uma taxa de ausência de evento cardíaco em 30 meses de 98%. A combinação de duas (classe B) e três variáveis (classe C) resultou em uma taxa de ausência de evento cardíaco em 30 meses de 65% e 0%, respectivamente (Figura 2).
Tabela 4 Tabela 4 - Escores preditores de mortalidade cardíaca e transplante na terapia de ressincronização cardíaca Escore no tempo 1 (pré-implante)
| Variável | Hazard | N | Escores | Classe | Risco |
|---|---|---|---|---|---|
| Nenhuma | 1,0 | 45 | 0 | A, | Baixo |
| DVD | 3,9 | 8 | 3 | Baixo | |
| FE | 4,8 | 14 | 4 | Baixo | |
| Diurético ↑ | 5,9 | 17 | 5 | Baixo | |
| DVD + FE | 19,1 | 5 | 7 | B | Intermediário |
| DVD + Diurético ↑ | 23,6 | 4 | 8 | B | Intermediário |
| FE + Diurético ↑ | 29,0 | 6 | 9 | B | Intermediário |
| DVD + FE + Diurético ↑ | 114,0 | 5 | 12 | C | Alto |
| Escore no tempo 2 (1° ano) | |||||
| Variável | Hazard | N | Escores | Classe | Risco |
| Nenhuma | 1,0 | 62 | 0 | A | Baixo |
| DVD | 3,5 | 7 | 2 | A | Baixo |
| Diurético ↑ | 5,3 | 12 | 3 | A | Baixo |
| Hospitalização | 12,5 | 3 | 5 | A | Baixo |
| DVD + Diurético ↑ | 18,7 | 2 | 6* | B | Intermediário |
| DVD + Hospitalização | 44,0 | 2 | 7 | B | Intermediário |
| Diurético ↑ + Hospitalização | 66,3 | 6 | 8 | B | Intermediário |
| DVD + Hospitalização + Diurético ↑ | 234,0 | 6 | 10 | C | Alto |
| Escore no tempo 3 ( 2° ano) | |||||
| Variável | Hazard | N | Escores | Classe | Risco |
| Nenhuma | 1,0 | 55 | 0 | A | Baixo |
| CF III/IV | 7,7 | 10 | 8 | B | Intermediário |
| DVD | 12,1 | 4 | 13 | B | Intermediário |
| DVD + CF III/IV | 94,5 | 5 | 21 | C | Alto |
DVD: disfunção ventricular direita; FE: fração de ejeção < 25%; diurético ↑: uso de furosemida ≥ 80 mg; CF: classe funcional (NYHA); hospitalização: ≥ 1 hospitalização por insuficiência cardíaca congestiva; classe A: risco baixo; classe B: risco intermediário; classe C: risco alto.
O risco foi usado como variável independente no modelo de regressão logística para a preparação do escore. O escore foi obtido pela hazard ratio da variável dividida pelo valor mais alto.
* adicionou-se uma unidade para manter a proporção de risco. N: número de pacientes.
Análise das variáveis no tempo 3 (segundo ano após TRC)
Hospitalizações por ICC, uso de altas doses de diuréticos, CF, DD, DVD, FE <30%, doença de Chagas e pressão arterial sistólica <110 mmHg mostraram-se significativas no modelo de regressão univariada de Cox no segundo ano após TRC.
No modelo multivariado de Cox, DVD e CF III/IV relacionaram-se de modo independente com o desfecho mortalidade/Tx cardíacos, com razões de risco de 7,7 e 12,0, respectivamente. As variáveis significativas no modelo multivariado também mostraram-se significativas em separado no modelo de Kaplan-Meier quando comparadas usando-se o teste log-rank (p < 0,001). A análise do modelo usando a curva ROC mostrou uma AUC de 0,789, com sensibilidade de 40%, especificidade de 98,4% e acurácia de 90,5% (Figura 3).
Figura 3 A: curva ROC do modelo no tempo 3 (segundo ano), com as variáveis disfunção ventricular direita (VD) e classe funcional (CF) III e IV comparada com I e II, área sob a curva (AUC) de 0,789, sensibilidade de 40%, especificidade de 98,4% e acurácia de 90,5%. B: a curva de Kaplan-Meier mostra que a ausência das variáveis disfunção VD e CF III / IV (classe A - baixo risco) indica uma taxa livre de eventos em 30 meses de 97.5%.
A partir de combinações dessas variáveis, desenvolvemos um modelo com três classes: classe A (baixo risco de morte/Tx cardíacos), ausência das duas variáveis significativas na análise multivariada, correspondendo a uma taxa de ausência de evento cardíaco em 30 meses de 97,5%. A combinação de duas (classe B) e três variáveis (classe C) resultou em uma taxa de ausência de evento cardíaco em 30 meses de 83,1% e 38,5%, respectivamente.
Discussão
No presente estudo, desenvolvemos três modelos preditores para risco de morte e Tx cardíacos em diferentes estágios da TRC. Até onde sabemos, este é o primeiro estudo a analisar sequencial e prospectivamente variáveis preditores na mesma população e em diferentes estágios de desenvolvimento (pré-implante, e no primeiro e segundo anos após TRC) e a desenvolver modelos de risco para morte e Tx cardíacos. Os modelos identificaram variáveis simples que, quando presentes, associaram-se a alto risco de morte e Tx cardíacos.
A mortalidade total foi de 25% (29/116) aos 34 ± 17 meses. No estudo CARE-HF9, a mortalidade foi 30% no grupo sem intervenção, comparada a 20% no grupo com TRC, durante o seguimento de 29,4 meses. No estudo COMPANION10, a taxa de mortalidade foi 21% (131/617) no grupo TRC, comparada a 25% (77/308) no grupo controle, durante seguimento de 24 meses. Portanto, nossa mortalidade total acha-se dentro da variação descrita nos grandes estudos. Analisamos os desfechos combinados de mortalidade e Tx cardíacos, visando a identificar variáveis mais específicas relacionadas aos resultados da TRC11.
Vários estudos avaliaram fatores preditores de resposta ou morte em diferentes populações e com diferentes critérios de resposta, sendo os resultados inconsistentes. Entretanto, várias publicações identificaram os seguintes fatores preditores de resposta: miocardiopatia dilatada12; largura do QRS13; estreitamento do QRS14; presença de dissincronismo15; sexo feminino16; tipo de bloqueio de ramo17; diâmetro do ventrículo esquerdo18; integral da velocidade aórtica no tempo13; e DD19.
Os pacientes com DVD (20,9% do grupo) apresentaram pior evolução em todos os tempos da análise. Entretanto, observamos que seis pacientes com boa evolução apresentaram regressão das alterações no ventrículo direito. O estudo de Praus e cols.20 mostrou que a regressão do ventrículo direito ocorreu mais tarde (15 meses), enquanto Leong e cols.21 identificaram o ventrículo direito como um preditor independente de mortalidade. Logo, pacientes com DVD não deveriam ser excluídos da indicação para TRC, embora representem um subgrupo de maior risco para morte cardíaca ou Tx cardíaco após TRC22. A importância do ventrículo direito na TRC foi demonstrada em outros estudos recentes, mas não na elaboração dos modelos de risco para diferentes estágios de evolução23, 24.
Treze pacientes (11,2%) tinham doença de Chagas, dos quais cinco tinham DVD. A miocardiopatia chagásica foi relacionada a aumento de mortalidade na curva de sobrevida, à semelhança de outro estudo que a relacionou a pior desfecho25. Na análise multivariada, a doença de Chagas não se manteve como uma variável independente, provavelmente porque 41% dos pacientes tinham DVD, uma variável significativa em todas as ocasiões da análise. Logo, a relevância da DVD não esteve exclusivamente relacionada à doença de Chagas, pois 19 pacientes tinham DVD devida a outras etiologias.
Uma FE < 25% pré-implante identificou um subgrupo de mais alto risco para morte cardíaca. Linde e cols.26, em uma subanálise do estudo REVERSE, mostraram que uma FE basal < 30%, comparada a valores entre 30% e 40%, relacionou-se positivamente com sobrevida. Entretanto, Kronborg e cols.27 relataram que FE basal < 22,5% determinou maior mortalidade após TRC.
A hospitalização por ICC mostrou-se uma variável independente quanto à predição de mortalidade/Tx cardíacos no primeiro ano após TRC. O estudo representa, até onde sabemos, a primeira vez em que aquela variável foi incluída como independente na análise do risco de mortalidade no primeiro ano após TRC e não como uma parte do desfecho combinado com morte. Hospitalização por ICC é um bem conhecido fator de risco para miocardiopatia, tendo vários estudos demonstrado redução na incidência desse evento após TRC9, 10.
Outra variável clínica facilmente obtida que mostrou significativo valor no período pré-implante e primeiro ano após TRC foi o uso de altas doses de diurético de alça (furosemida ≥ 80 mg/dia). Van Boven e cols.28relataram uma associação entre o não uso crônico de diuréticos e a resposta à TRC. Entretanto, Cleland e cols.29observaram que o uso de altas doses de diurético relacionou-se a pior prognóstico apenas na análise univariada. Acreditamos que a descrição dessa variável como tendo um valor independente para morte cardíaca em dois momentos da análise de TRC no nosso estudo seja original.
Recentemente elaborou-se uma regra clínica de predição para identificar pacientes com risco elevado de morte precoce após TRC30, incluindo as seguintes quatro variáveis independentes: diâmetro diastólico final do ventrículo esquerdo (DDVE) > 65 mm, morfologia de não-BRE, nível de creatinina > 1,5, e não uso de betabloqueadores. No nosso estudo, DDVE e nível de creatinina foram significativos apenas nas análises univariadas. Hospitalização por ICC, uso de altas doses de diurético e DVD, algumas variáveis independentes no nosso estudo, não foram incluídas no estudo anterior.
Conseguimos uma significativa melhora na especificidade dos modelos preditores de mortalidade ou resposta após TRC, alcançando 96% no primeiro ano após TRC e 98% no segundo ano após TRC, quando comparada com especificidade de 22%-70% dos modelos previamente descritos em relação às mortalidades total e cardíaca. Tais resultados estão de acordo com os desfechos-alvo da TRC no tratamento de pacientes com doenças graves, com altos custos e riscos no procedimento31.Os modelos usados neste estudo mostraram boa acurácia, variando de 84,6% a 93%, podendo ser usados em três diferentes estágios de TRC, outra contribuição original do nosso estudo. Nos níveis de significância usuais, o modelo foi validado internamente e não revelou falta de ajuste ou sensibilidade exagerada dos dados.
Acreditamos que este estudo contribua, promovendo o avanço da busca por melhores critérios de avaliação prognóstica, com a composição de índices multifatoriais simples e com a inclusão de variáveis usadas na prática clínica de fácil obtenção. Os modelos serão úteis para seleção, monitoramento e orientação de pacientes com indicação para TRC.
Limitações do estudo
Não foram realizadas análises da variabilidade intra- e interobservador de parâmetros ecocardiográficos e eletrocardiográficos. Os pacientes não foram submetidos a otimização do intervalo atrioventricular após a cirurgia. Os modelos criados não foram validados externamente, apenas internamente. Outra limitação deste estudo é o pequeno número de pacientes, o grande número de excluídos e o fato de ter sido conduzido em um único centro. A função do ventrículo direito foi analisada qualitativamente devido à ausência de correlação entre as medidas desse ventrículo e o prognóstico no início do estudo. Tais resultados devem ser considerados no contexto da população do estudo: 59,4% tinham miocardiopatia dilatada, 11,2% tinham miocardiopatia chagásica, 12% tinham BRD e 16,3% tinham implante prévia de marca-passo cardíaco. Maiores estudos prospectivos são necessários para validar as importantes variáveis relacionadas a morte e Tx cardíacos após TRC.
Conclusão
Desenvolvemos modelos preditores de morte ou Tx cardíacos em diferentes estágios da TRC com base na análise de variáveis clínicas e ecocardiográficas simples e de fácil obtenção. Os modelos mostram boa acurácia e ajuste, foram validados internamente e são úteis para a seleção, o monitoramento e a orientação de pacientes indicados para TRC.