INTRODUCCIÓN
La enfermedad renal crónica (ERC) se ha convertido en la tercera causa de muerte de más rápido crecimiento y se estima que para el 2040 sea la quinta causa de años de vida potencialmente perdidos. A nivel mundial, alrededor de 850 millones de personas conviven con la enfermedad renal, de estos, 4 millones reciben terapias de reemplazo renal para una tasa bruta de 298,4 pacientes por millón de habitantes.1
En Cuba, durante el año 2021 la ERC fue la causa número 14 de mayor mortalidad. Se reportaron 1576 muertes, 398 más que el año 2020, para una tasa bruta de 14,1 y ajustada de 7,0.2
Los países desarrollados presentan adecuados sistemas de datos de los pacientes con enfermedad renal. Es conveniente extrapolar estos resultados a los países en vías de desarrollo y generar un cambio de paradigma que disminuya las inequidades.3
La mortalidad por ERC en Cuba ha aumentado de forma considerable con respecto a años anteriores por lo que constituye una necesidad científica la realización de investigaciones que contribuyan a identificar de manera precoz factores pronósticos. El envejecimiento poblacional y el aumento de enfermedades como la hipertensión arterial y la diabetes mellitus son desafíos importantes, en aras de obtener mejores resultados en la salud renal de los individuos.4
En la actualidad, se requieren nuevos modelos con diseño metodológico y relaciones de causalidad diferentes que investiguen la relación de variables y factores para mejorar el pronóstico de la disfunción renal. La implementación de técnicas de big data, ciencia de datos y otras herramientas de la inteligencia artificial podrían generar mejoras sobre las crecientes cifras de mortalidad.5
Realizar pronósticos certeros y de forma precoz es clave en la evolución de la historia natural de la ERC y en la estratificación del riesgo de los pacientes. En la actualidad el pronóstico de esta enfermedad cobra mayor importancia si se aplican nuevos avances metodológicos que aporten consistencia sobre el desarrollo de modelos predictivos de mortalidad más reproducibles y con enfoque uniforme.6
Existen numerosas escalas para la predicción de mortalidad en la ERC, pero a pesar de poseer gran sensibilidad y especificidad, en Cuba presentan limitaciones reales ya que algunas de las variables que reúnen no se pueden medir. Es por ello que es de gran importancia el desarrollo de nuevas herramientas pronósticas más homogéneas que se puedan utilizar para la toma de decisiones clínicas importantes.
La presente investigación se realizó con el objetivo de diseñar una escala predictiva de mortalidad en pacientes con enfermedad renal crónica.
MÉTODOS
Contexto y clasificación de la investigación
Se realizó un estudio observacional, analítico, longitudinal prospectivo en los pacientes con enfermedad renal crónica estadio 4 atendidos en la consulta de enfermedad renal crónica avanzada (ERCA) de la provincia Las Tunas durante el período del 1 de enero de 2022 al 31 de diciembre de 2022. Durante los primeros 6 meses del año se analizaron las variables de estudio para el diseño de la escala predictiva. En los próximos 6 meses, los pacientes fueron seguidos para identificar la ocurrencia o no de la variable dependiente mortalidad.
Población y muestra
La población estuvo constituida por 302 pacientes con enfermedad renal crónica estadio 4, atendidos durante el período declarado anteriormente. Para calcular el tamaño de la muestra se utilizó la fórmula para muestras finitas o conocidas.7
Donde "n" representó el total de pacientes con enfermedad renal crónica que se necesitaron incluir en la muestra, Z= 1,96 para un nivel de confianza de 95 %, con un máximo error permisible del 5 % (e= 0,05), el valor de p y q se planteó en 0,05 teniendo en cuenta que no existen trabajos actuales publicados acerca de la tasa de mortalidad de la enfermedad renal crónica en Cuba.8) Se necesitaron incluir 169 pacientes en la muestra. Se estudiaron 62 pacientes fallecidos y 117 pacientes vivos, con una relación 1:2 mediante un muestreo aleatorio simple.
Criterios de inclusión
Constituyeron criterios de inclusión los pacientes mayores de 18 años de edad con ERC estadio 4 atendidos en la consulta de ERCA.
Variables en la investigación
Para la descripción de la muestra se analizaron un total de 20 variables divididas en variables epidemiológicas, clínicas y exámenes complementarios.
Variable dependiente: mortalidad.
Variables independientes:
Sexo: (0) femenino, (1) masculino.
Edad: según los años cumplidos al inicio de la investigación. Para el análisis de regresión se consideraron 2 categorías, (0) igual o menor a 65 años, (1) mayor a 65 años.9)
Enfermedad cardiovascular: si el paciente presentó enfermedad coronaria, insuficiencia cardíaca o arritmia. Para el análisis de regresión se consideraron 2 categorías, (0) ausencia, (1) presencia.10)
Diabetes mellitus: si el paciente presentaba el antecedente de la enfermedad; (0) ausencia, (1) presencia.
Hipertensión arterial: si el paciente presentaba el antecedente de la enfermedad y los que sin tener el antecedente, presentaron cifras tensionales superiores a 140/90 mmHg en al menos 2 ocasiones; (0) ausencia, (1) presencia.
Tabaquismo: si el paciente en el momento del estudio fumaba; (0) no fumador, (1) fumador.11)
Dislipidemia: si el paciente tenía el antecedente o si presentó cifras de colesterol superior a 5,7 mmol/L o triglicéridos superior a 1,7 mmol/L. Para el análisis de regresión se consideraron 2 categorías; (0) ausencia, (1) presencia.12
Obesidad: si el índice de masa corporal era superior a 30 kg/m2. El cálculo se realizó según la fórmula de Quetelet;13 (0) ausencia, (1) presencia.
Desnutrición: si el índice de masa corporal era inferior a 18,5 kg/m2. El cálculo se realizó según la fórmula de Quetelet;13 (0) ausencia, (1) presencia.13)
En el caso de las variables de los exámenes complementarios, el valor que se tomó para la investigación fue la primera determinación realizada en la primera consulta. En cuanto a las variables cuantitativas, se incluyeron en el análisis de regresión, albúmina < 35 g/L, hemoglobina < 10 g/L y ácido úrico > 390 mmol/L. Estos puntos de corte se consideraron según la revisión de la literatura.14)
Procesamiento estadístico
Los datos obtenidos fueron almacenados y procesados en el software, SPSS versión 20 para Windows. Para determinar la distribución normal o no de las variables cuantitativas se aplicó el test Kolmogorov Smirnov. Aquellas variables cuantitativas con distribución normal se expresaron en media y desviación estándar y sus medias fueron comparadas a través del test estadístico t de Student. Las variables que no presentaron distribución normal se expresaron en mediana y rango intercuartil y se realizó la comparación con el test U de Mann Whitney.
Para determinar las diferencias entre los grupos establecidos según variables cualitativas se utilizó el test estadístico ji cuadrado basado en su hipótesis de homogeneidad.
Luego, se implementó un análisis multivariado de regresión logística binaria para determinar las variables que constituyeron factores pronósticos de mortalidad. El proceso de selección de las variables para la inclusión en la regresión logística se realizó teniendo en cuenta la significación estadística de los atributos en el análisis univariado y aquellas variables que la literatura plantea como factores pronósticos de mortalidad; así se evitó el sesgo de confusión. Luego, aquellas variables que presentaron problemas de colinealidad no se incluyeron en el modelo. Por cada variable independiente que se añadió al modelo se tuvo en cuenta 10 frecuencias de la variable dependiente mortalidad.15
Se trabajó con una confiabilidad del 95 %.
Diseño de la escala predictiva16
A partir de los resultados de la regresión logística, se tomaron como pesos o ponderaciones para cada variable, los odds ratios (OR) redondeados a una sola cifra, OR= Exp (B); de esta forma, la escala predictiva resultante quedó compuesta por un índice global en forma de una combinación lineal: I= x1w1+ x2w2+…+xkwk, donde I constituye la escala propuesta; k, el número de variables, xk es la variable predictiva; y wk, el peso escogido para dicha variable en la escala.
Se calculó el valor de la escala predictiva para cada paciente, y con el propósito de buscar estratos de riesgos diferentes, se procedió a determinar los percentiles 25 y 50 de la distribución empírica de los valores de la escala.16) Quedaron conformadas 3 categorías: riesgo bajo, riesgo mediano y riesgo alto.
Luego se evaluó la capacidad discriminatoria de la herramienta confeccionada mediante la estimación del área bajo la curva del estadístico C. Además, se evaluaron las curvas de supervivencia a los 6 meses según la puntuación de la escala obtenida por cada paciente, las medias de supervivencia se compararon mediante el test Log-Rank.
Aspectos éticos
Se cumplió con las normas éticas de la Declaración de Helsinki.17 Se utilizaron las historias clínicas, y otros documentos que pudieran ofrecer datos de relevancia para la realización de la investigación. Además, la investigación fue aprobada por el comité de ética de la provincia Las Tunas.
RESULTADOS
Como se muestra en la tabla 1 se analizaron 169 pacientes, la mediana de edad global fue de 60 años (RIQ: 19), la mediana de edad fue mayor en el grupo de los pacientes fallecidos (p= 0,000); predominó el sexo masculino (68,00 %). En cuanto a los antecedentes, el 50 % de los pacientes tenían hipertensión arterial y el 28,40 % diabetes mellitus. La presencia de enfermedad cardiovascular (p= 0,000) y la dislipidemia (p= 0,000) tuvieron diferencias con significación estadística.
En relación con las variables clínicas y exámenes complementarios según ambos grupos, no existieron diferencias significativas entre la tensión arterial sistólica, tensión arterial diastólica, la frecuencia cardíaca y el ácido úrico. Sin embargo, las cifras de hemoglobina (p= 0,000), albúmina (p= 0,000), urea (p= 0,009), creatinina (p= 0,002) y filtrado glomerular (p= 0,000), sí tuvieron diferencias estadísticas significativas. Por otra parte, las medianas de glucemia y colesterol presentaron significación estadística (p= 0,001) y (p= 0,000) respectivamente.
En la tabla 2 se muestran las variables pronósticas de mortalidad en la regresión logística binaria: edad > 65 años [B= 1,831; p= 0,000; Exp (B)= 2,118], enfermedad cardiovascular [B= 2,831; p= 0,000; Exp (B)= 2,930], albúmina < 35 g/L [B= 1,839; p= 0,000; Exp (B)= 2,166], dislipidemia [B= 1,326; p= 0,027; Exp (B)= 1,458], hemoglobina < 10 g/L [B= 1,600; p= 0,015; Exp (B)= 1,324] y ácido úrico > 390 mmol/L [B= 1,243; p= 0,001; Exp (B)= 1,420]. Entre las pruebas utilizadas para evaluar el rendimiento del modelo de predicción, se muestra buena calibración. El valor del R cuadrado de Nagelkerke igual a 0,803 traduce que el 80,30 % de la varianza de la variable dependiente es explicada por el modelo.
Tabla 1 - Distribución de los pacientes según características epidemiológicas, clínicas y exámenes complementarios
|
|
p | ||
---|---|---|---|---|
Sexo masculino | 115 (68,00 %) | 40 (64,50 %) | 75 (70,10 %) | 0,454 |
Edad a | 60 (19) | 65 (13) | 54 (21) | 0,000b |
Enfermedad cardiovascular | 41 (24,30 %) | 36 (58,10 %) | 5 (4,7 %) | 0,000b |
Diabetes mellitus | 48 (28,40 %) | 22 (35,50 %) | 26 (24,30 %) | 0,120 |
Hipertensión arterial | 85 (50,30 %) | 34 (54,80 %) | 51 (47,70 %) | 0,369 |
Tabaquismo | 34 (20,10 %) | 11 (17,70 %) | 23 (21,50 %) | 0,557 |
Dislipidemia | 34 (20,10 %) | 25 (40,30 %) | 9 (8,4 %) | 0,000b |
Obesidad | 17 (10,10 %) | 7 (11,30 %) | 10 (9,3 %) | 0,685 |
Desnutrición | 45 (26,60 %) | 19 (30,60 %) | 26 (24,30 %) | 0,368 |
Tensión arterial sistólica c | 123,46±14,121 | 122,42±10,664 | 124,07±15,798 | 0,421 |
Tensión arterial diastólica c | 74,11±12,587 | 71,77±11,842 | 75,47±12,859 | 0,066 |
Frecuencia cardiaca c | 70,99±10,484 | 69,58±10,462 | 71,81±10,458 | 0,183 |
Hemoglobina a | 10 (3,0) | 10 (2,0) | 11 (2,0) | 0,000b |
Albúmina a | 37 (12) | 29 (12) | 40 (7) | 0,000b |
Urea a | 16 (5,0) | 17 (4,0) | 13 (4,0) | 0,009b |
Creatinina a | 258 (23) | 296 (34) | 224 (19) | 0,002b |
Filtrado glomerular c | 22,905±1,9798 | 17,274±1,3810 | 26,274±2,1787 | 0,000b |
Ácido úrico c | 351,86±96,432 | 347,73±80,949 | 354,01±103,859 | 0,696 |
Glucemia a | 4 (1,0) | 4 (1,8) | 5 (2,0) | 0,001b |
Colesterol a | 4 (2,0) | 4 (1,0) | 4 (1,0) | 0,000b |
a Resultado expresado en mediana y rango intercuartil. b Estadístico significativo. c Resultado expresado como media ± DE.
Tabla 2 - Regresión logística para la predicción de variables pronósticas de mortalidad
B | ||||
---|---|---|---|---|
Edad > 65 años | 1,831 | 0,000 | 2,118 | 1,755-4,326 |
Enfermedad cardiovascular | 2,831 | 0,000 | 2,930 | 1,906-5,720 |
Albúmina < 35 g/L | 1,839 | 0,000 | 2,166 | 1,810-4,995 |
Dislipidemia | 1,326 | 0,027 | 1,458 | 1,147-3,652 |
Hemoglobina < 10 g/L | 1,400 | 0,015 | 1,324 | 1,087-3,008 |
Ácido úrico > 390 mmol/L | 1,243 | 0,001 | 1,420 | 1,130-3,745 |
Constante | 4,209 | 0,000 | 0,000 | - |
Se calculó el valor de la escala predictiva para cada paciente y con el propósito de buscar estratos de riesgos diferentes, se procedió a determinar los percentiles 25 y 50 de la distribución empírica de los valores de la escala, los cuales fueron 4 y 7.16) Luego quedaron conformadas 3 categorías: riesgo bajo, pacientes con valores menores a 4; riesgo mediano, pacientes con valores entre 4 y 7 y riesgo alto, pacientes con valores superiores a 7 (tabla 3).
Tabla 3 - Ponderaciones de las variables que conformaron la herramienta predictiva
Edad > 65 años | 2 |
Enfermedad cardiovascular | 3 |
Albúmina < 35 g/L | 2 |
Dislipidemia | 1 |
Hemoglobina < 10 g/L | 1 |
Ácido úrico > 390 mmol/L | 1 |
< 4 puntos: riesgo bajo
5-7 puntos: riesgo mediano
> 7 puntos: riesgo alto
El valor correspondiente al área bajo la curva del estadístico C obtenida fue superior a 0,70, por tanto, se puede considerar que la herramienta predictiva tiene valor predictivo significativo. El poder discriminatorio de la escala, en aquellos pacientes con puntuaciones superiores a 4 puntos para predecir mortalidad fue bueno, área bajo la curva COR 0,856 (IC 95 %: 0,783-0,929) (Fig. 1).
En la figura 2 se muestran las funciones de supervivencia de los pacientes a los 6 meses según los valores obtenidos de la escala predictiva. Existieron diferencias significativas al comparar a los pacientes según el valor de la escala estuviese por encima o no del punto de corte. La tasa de mortalidad fue mayor en aquellos pacientes que presentaron valores iguales o mayores a 4 puntos con respecto a aquellos con cifras menores (65,75 % vs. 14,58 %, p< 0,001).
Los pacientes con valores menores a 4 puntos presentaron media de supervivencia de 149,438 ± 7,296 días. En cambio, los que tenían valores iguales o superiores a 4 puntos presentaron media de supervivencia de 93,128 ± 8,545 días.
DISCUSIÓN
Los resultados de la presente investigación muestran que la escala predictiva propuesta es capaz de discriminar aquellos pacientes con ERC avanzada con mayor riesgo de mortalidad.
Schmidt y otros18) en un estudio de cohorte prospectivo que analizó 18 282 pacientes con ERC avanzada con el objetivo de desarrollar un modelo de predicción de mortalidad a los 12 meses y en el que se realizó regresión logística binaria demostraron que la edad y la albúmina tenían una relación lineal con la mortalidad.
Como se evidencia en el estudio anterior, el modelo de regresión logística binaria es seguro y efectivo en la predicción de mortalidad de pacientes con ERC. Los valores de los indicadores de eficacia estadística varían teniendo en cuenta el tipo y la calidad de las variables analizadas, por lo que un modelo de regresión logística basado en variables seleccionadas de contextos propios puede identificar nuevos factores pronósticos y optimizar los ya conocidos.
La edad avanzada tiene un rol esencial tanto en la progresión como en la mortalidad de los pacientes con ERC. En este estudio los pacientes fallecidos presentaron mediana de edad superiores a los vivos (65 años vs. 54 años).
Un estudio de cohortes transversales y prospectivas que incluyó pacientes de China y EE. UU. en el cual se incluyeron personas mayores de 65 años demostró que los pacientes con edad avanzada tenían peores niveles de biomarcadores relacionados con la enfermedad renal y la mortalidad fue mayor.19)
La enfermedad cardiovascular es la primera causa en la mortalidad de los pacientes con ERC. Los pacientes con disfunción renal crónica presentan calcificación vascular, ateroesclerosis avanzada y rigidez valvular. Estos fenómenos son capaces de favorecer el desarrollo de hipertrofia ventricular izquierda, insuficiencia cardíaca, arritmias y paro cardíaco.20)
La tendencia a la disminución de la albúmina sérica se asocia a los resultados adversos y a la mortalidad. Zeng y otros21) investigaron la asociación entre las concentraciones de albúmina sérica y el riesgo de paro cardíaco en pacientes con ERC. En los 4990 pacientes estudiados, por cada 1 g/dL de aumento de la albúmina se relacionó con una disminución del 68 % del riesgo de mortalidad cardiovascular.
En esta investigación se identificó como variable pronóstico, la hemoglobina menor a 10 g/dL. En el grupo fallecidos y de forma general las cifras de hemoglobina estuvieron disminuidas. La anemia en el paciente con ERC es de causa multifactorial pero el uso de agentes estimulantes de la eritropoyetina y suplementos de hierro inciden considerablemente en los resultados de la supervivencia.
Wittbrodt y otros22) caracterizaron 22 720 pacientes con ERC en un estudio longitudinal retrospectivo realizado en EE. UU. El 23,3 % de los pacientes presentó anemia, la prevalencia de esta fue mayor en los pacientes en estadio 4 y aquellos en terapias en reemplazo renal y se asoció a mayor tasa de hospitalizaciones, progresión de la enfermedad renal y mortalidad.
La dislipidemia y el ácido úrico > 390 mmol/L fueron otras variables que se asociaron de forma significativa a la mortalidad. En cuanto al ácido úrico, es bien reconocida la importancia de su análisis de forma longitudinal, ya que su elevación se asocia a mayor riesgo de diálisis y mortalidad. Además, presenta una forma de "U" en la supervivencia teniendo en cuenta que aquellos pacientes con valores elevados presentan riesgo incrementado de mortalidad. Las sales de urato ocasionan disfunción endotelial, activación del sistema renina angiotensina aldosterona y atrofia tubular renal.
La presente investigación presentó como limitación que se utilizó un tamaño de muestra relativamente pequeño relacionado con la cantidad real de pacientes con ERC avanzada en el contexto estudiado, aun así, los resultados obtenidos pueden ser considerados en el seguimiento de los pacientes con enfermedad renal crónica. En la actualidad, se implementa el modelo desarrollado para su aplicación práctica en las consultas de enfermedad renal crónica avanzada.
La escala predictiva compuesta por 6 variables contribuye a la estratificación del riesgo de mortalidad de los pacientes con ERC. La escala muestra buena capacidad discriminatoria y calibración. Las variables incluidas son de fácil determinación e interpretación por lo que se considera como un modelo de útil implementación en la toma de decisiones médicas en el ámbito clínico actual.
Se recomienda realizar investigaciones de validación interna y externa de la escala predictiva en cohortes con mayores muestras de estimación.