Servicios personalizados
Español (pdf)
Articulo en XML
Referencias del artículo
Enviar articulo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
La drepanocitosis es una hemoglobinopatía de alta prevalencia mundial. La entidad está asociada con enfermedad cerebrovascular, hemólisis, crisis vasoclusiva (CVO), hipertensión pulmonar, entre otras.1,2) Junto al mecanismo clásico (falciformación-vasoclusión) se conoce que la hemólisis de los drepanocitos produce una importante disfunción endotelial por múltiples mecanismos que incluyen aumento de la hemoglobina libre, desregulación del metabolismo de la arginina, aumento del estrés oxidativo, disminución local de óxido nítrico, liberación de mediadores inflamatorios (interleukina 1, bradikinina, histamina y sustancia P, así como prostaglandina E2 e I2).3
Se conoce que la hiperproducción de especies reactivas de oxígeno (EROs) y de nitrógeno (ENOs) forma parte de su fisiopatología e inducen un estado de estrés oxidativo que se asocia con disfunción endotelial y daño tisular. Higuera y col evaluaron algunos biomarcadores de estrés oxidativo y su asociación con la gravedad clínica y encontraron un aumento significativo de sustancias reactivas al ácido tiobarbitúrico (TBARS) en los pacientes HbSS con fenotipo grave.4
El estado de CVO se diagnostica tardíamente debido a que los pacientes acuden a los servicios de urgencia, cuando el daño se hace evidente a través de las manifestaciones clínicas. Por lo que los especialistas carecen de una herramienta diagnóstica que les permita seguir la evolución de estos pacientes en períodos de estabilidad y detectar tempranamente el incremento del riesgo hacia la CVO. La aplicación en medicina de screening combinados enfocados a la predicción de diversos estados patológicos es un hecho en la actualidad.5,6,7
Entre los métodos utilizados con este fin, el análisis discriminante, es definido como “ una técnica estadística multivariante cuya finalidad es analizar si existen diferencias significativas entre grupos de objetos respecto a un conjunto de variables medidas sobre los mismos para, en el caso de que existan, explicar en qué sentido se dan y facilitar procedimientos de clasificación sistemática de nuevas observaciones de origen desconocido en uno de los grupos analizados”.8) El propósito básico del análisis discriminante es identificar el grupo de pertenencia de un sujeto a partir de la estimación de la relación entre una única variable dependiente no métrica (categórica) que suele ser dicotómica y un conjunto de variables independientes cuantitativas métricas, en esta forma general (ecuación 1):
El análisis implica obtener un valor teórico, es decir, una combinación lineal de dos o más variables independientes que mejor discriminen entre los grupos definidos a priori. La combinación lineal para el análisis discriminante, también conocida como función discriminante, se deriva de la ecuación (2)que adopta la siguiente forma, para el caso de dos grupos:
En la anemia drepanocítica, se desconoce la existencia de modelos matemáticos que establezcan el estado riesgo de estos pacientes hacia el estado de CVO. La experiencia obtenida de etapas previas de la investigación, con la caracterización de los estados asintomáticos y de CVO a partir de resultados del perfil lipídico y estado redox en una muestra de pacientes, justifica el estudio actual.
El objetivo del estudio consistió en establecer índices de riesgo en pacientes con drepanocitosis, a partir de la formulación de un modelo predictivo del estado vasoclusivo.
Métodos
Se realizó un estudio analítico transversal, en el Laboratorio de Ciencias Básicas Biomédicas (LCBB) de la Universidad de Ciencias Médicas en Santiago de Cuba, en pacientes con drepanocitosis en estados asintomático y de CVO, atendidos en el Centro Hematológico del Hospital Clínico Quirúrgico “Juan Bruno Zayas Alfonso” de Santiago de Cuba, durante el período de mayo-septiembre de 2021.
El criterio de inclusión consideró a un total de 76 pacientes con drepanocitosis; los que acudían a la consulta (estado asintomático) y al servicio de urgencias con sintomatología propia del estado vasoclusivo. Resultaron excluidos aquellos pacientes que padecían de neuropatías, cardiopatías y nefropatías, los cuales podrían contribuir a un sesgo en la investigación.
Se recogió, bajo la autorización del comité de ética y consejo científico del hospital, el consentimiento informado de los participantes, cumpliéndose estrictamente los principios éticos establecidos en las normas relacionadas con el manejo de información médica. De igual modo, se aplicaron los parámetros que determinan el reglamento ético interno del citado hospital. Los datos de identificación de los pacientes fueron anónimos, siguiendo los principios éticos de la investigación científica con seres humanos. Una entrevista estructurada incluyó sus datos generales, hábitos alimentarios-nutricionales, tóxicos y aspectos clínicos.
En una primera etapa del estudio se realizó una caracterización de la población de estudio atendiendo al perfil lipídico y el estado redox.9 En ese momento, mediante un análisis univariado resultaron variables de interés, en cuanto a la diferenciación de ambos grupos de estudio, los triacilglicéridos y HDL-colesterol, así como malonildialdehído (MDA) como indicador de daño oxidativo a lípidos y el glutatión reducido (GSH) y la catalasa (CAT), expresando respuesta antioxidante en los pacientes. A partir de los resultados obtenidos, se procedió en la presente etapa a la formulación del modelo predictivo e índices de riesgo del estado de CVO según se describe:
Métodos y procedimientos
Se consideraron las variables clínicas y de interés específico, a partir de entrevistas e historias clínicas de los pacientes. La selección estuvo basada en factores de riesgo de CVO y biomarcadores asociados a dicha entidad, a partir de una exhaustiva revisión bibliográfica. Los datos recolectados quedaron plasmados en una base de datos en el programa SPSS (V 21,0), que fue actualizada sistemáticamente.
En todos los contrastes de hipótesis realizados con técnicas estadísticas se aceptó la existencia de significación estadística para una confianza superior al 95 % (p< 0,05). El análisis estadístico de los datos fue realizado con el programa informático SPSS versión 21,0 (GNU General PublicLicense, Free Software Fundation).
Formulación del modelo predictivo
El procedimiento establecido en este trabajo se sustentó en la explicación detallada ofrecida por Torrado10 para su desarrollo por el programa SPSS según:
Selección de las variables dependientes e independientes :
En el análisis discriminante, la variable dependiente (presencia o no de CVO), es una variable categórica (dicotómica) que quedó explicada a través de las variables independientes (clínicas y bioquímicas), partiendo de los hallazgos reconocidos en la literatura para esta entidad.
Hipótesis del análisis :
Para asegurar la factibilidad del método de clasificación fueron inicialmente analizados los supuestos de normalidad multivariante, homocedasticidad multivariante y ausencia de multicolinealidad entre variables independientes.
Estimación de la función discriminante y valoración del ajuste global :
El método seleccionado para la estimación de la función fue stepwise (estimación “paso a paso”), el que permitió ir introduciendo las variables clasificadoras (desde la más significativa), hasta que se cumplió el criterio de detención. 10 Solo se incluyeron aquellas variables independientes que superaron ciertos niveles mínimos de poder explicativo según el estadígrafo Lambda de Wilks, asumiendo el ofrecido en el programa de un valor de F=3,84 (para entrar) y F= 2,71 (para salir).
Como medida de bondad de ajuste para explicar la capacidad explicativa de la función discriminante se usó la raíz cuadrada del “coeficiente eta cuadrado (ŋ2) o coeficiente de correlación canónica”. La evaluación de las diferencias entre grupos sobre las funciones discriminantes se determinó a partir de la comparación de los centroides (puntuación D discriminante media para los pacientes del grupo). La precisión en la predicción se determinó por medio de matrices de clasificación, que incluyó la proporción de aciertos y fallos, estableciendo el porciento de casos agrupados correctamente clasificados.
Elaboración de índices de riesgo de CVO
El conocimiento de los valores de media de las puntuaciones discriminantes arrojadas por los pacientes de cada grupo proporcionó la elaboración de los índices de riesgo que se proponen, partiendo de la siguiente metodología:
1) Se calcularon las 3 desviaciones estándar (SD) alrededor de la media (x̅; y̅) de cada grupo, para descartar la presencia de solapamiento entre ellos, que pudiera interferir en la estratificación de los índices.
2) Se realizó un ajuste de la metodología tradicional que ofrece el método, para este tipo de análisis en dos grupos, 11) se particionaron las distribuciones para tres categorías, calculando dos puntos de corte de manera que delimitaran tres zonas de riesgo equivalentes en amplitud.
3) Se establecieron tres zonas de riesgo de CVO: La primera zona (de menor valor) correspondería a la categoría de riesgo menor; seguida de la zona intermedia o de riesgo moderado y luego la de mayor valor para la categoría de riesgo mayor, siendo esta última la de mayor asociación con la ocurrencia de CVO.
Resultados
Un total de 76 pacientes con drepanocitosis fueron analizados mediante el uso de la estadística multivariada. El empleo del análisis discriminante proporcionó la formulación del modelo predictivo del estado de CVO, arribando a los siguientes resultados:
Modelo predictivo de CVO
El método “paso a paso” permitió incluir las variables independientes, una a una, según su capacidad discriminatoria. Una vez estimadas las funciones, se valoró el ajuste global del modelo según el valor que toman los estadígrafos que se muestran en Tabla 1. Como se observa, las variables que resultaron incluidas en la función acorde al valor Lambda de Wilks, fueron el MDA y el GSH.
Tabla 1 Variables introducidas en la función discriminante acorde al F exacto
| Pasos | Variables | Lambda de Wilks | |||||||
| Estadística | df1 | df2 | df3 | F Exacto | |||||
| Estadística | df1 | df2 | Sig. | ||||||
| 1 | MDA | 0,284 | 1 | 1 | 74,000 | 186,743 | 1 | 74,000 | 0,000 |
| 2 | GSH | 0,244 | 2 | 1 | 74,000 | 112,799 | 2 | 73,000 | 0,000 |
Como otra medida de evaluación de la bondad de ajuste de la función se utiliza el eta cuadrado (η2), coeficiente de determinación obtenido al realizar la regresión entre la variable dicotómica, que indica la pertenencia al grupo y las puntuaciones discriminantes (Tabla 2). El elevado autovalor (4,21) indica la proporción de varianza total explicada por la función discriminante, por lo que muestra cuán eficaz será el análisis para clasificar a los pacientes (alto poder discriminante de la función). El valor de la correlación canónica indicó la adecuada pertenencia de los mismos a los grupos conformados a priori al mostrar un 89,9%.
Tabla 2 Bondad de ajuste del modelo
| Función | Autovalor | % de varianza | % acumulado | Correlación canónica |
| Dimensión 1 | 4,216 | 100,0 | 100,0 | 0,899 |
Los centroides, como punto medio de la distribución de las puntuaciones individuales de los grupos, resumen la mejor clasificación. En la función discriminante resultante, los centroides de los grupos, mostraron valores distantes entre pacientes fuera de crisis (-1,418) y en crisis (2,894), que denota la excelente capacidad discriminante de la función.
Los resultados del proceso de clasificación se presentan de forma matricial (matriz de clasificación) (Tabla 3). Los elementos de la diagonal de la matriz representan el número de pacientes originalmente agrupados que resultaron correctamente clasificados (100 %), lo que brinda una evidencia de la robustez del modelo que se propone en función de la clasificación de los pacientes.
Tabla 3 Matriz de clasificacióna
| Clasificación | Predicción de grupos | Total | |||
| Paciente fuera de crisis | Paciente en crisis | ||||
| Original | n | Paciente fuera de crisis | 50 | 1 | 51 |
| Paciente en crisis | 2 | 23 | 25 | ||
| % | Paciente fuera de crisis | 98,0 | 2,0 | 100,0 | |
| Paciente en crisis | 8,0 | 92,0 | 100,0 | ||
a. 100% de los casos agrupados originalmente clasificados correctamente
El análisis discriminante identificó 2 de los 7 parámetros como discriminadores significativos del estado de CVO (p < 0,001): MDA y GSH. La función discriminante canónica obtenida fue la siguiente:
donde: GSH: glutatión reducido y MDA: malonildialdehído
Para brindar una clasificación estratificada en tres categorías de riesgo para esta condición (menor, moderado y mayor), se procedió al cálculo de índices de riesgo, con el siguiente resultado:
Índices de riesgo de CVO sustentados en el modelo predictivo
Los valores medios de cada población posibilitaron que se establecieran los índices de riesgo de CVO, para lo cual se calcularon dos puntos de corte (PC1 = 0,019 y PC2 = 1,456), a partir de su amplitud, quedando delimitada cada zona de riesgo de la siguiente manera:
D < 0,019 - Riesgo menor
D Є [0,019; 1,456] - Riesgo moderado
D > 1,456 - Riesgo mayor
Los índices expuestos ofrecen una ponderación del estado de riesgo de los pacientes, y brindan una aproximación al estado de salud real de estos. En la figura 1 se refleja la distribución de riesgos en los pacientes a partir de los índices obtenidos. Se muestra que a medida que aumenta el riesgo se incrementa el porcentaje de dicho evento, donde el 99 % de los casos de PE fueron clasificados con riesgo mayor, al mostrar un D>1,456. El resto de los casos clasificaron en la categoría de riesgo menor, asignando solo 1 paciente a la categoría de riesgo moderado.
Discusión
El modelo predictivo del estado de CVO en pacientes con drepanocitosis consideró como variables de interés al GSH y el MDA, al ser las que mantuvieron su significación estadística durante la clasificación en los grupos analizados. Estas variables, que ya habían mostrado su valor diagnóstico desde la estadística univariada, lo reafirman a través de este tipo de análisis de mayor complejidad.
La inclusión del MDA en la función, con su elevado coeficiente, revela su importante papel en la diferenciación entre grupos. Su influencia en el estado de CVO ha resultado un tema discutido en algunos estudios donde el proceso de peroxidación lipídica está íntimamente ligado a su ocurrencia, por activación de diferentes mecanismos moleculares involucrados.12,13La dislipidemia que acompaña a estos pacientes en el tránsito hacia la CVO constituye una importante fuente generadora de EO (sustrato oxidante) logrando superar, estas especies, la capacidad antioxidante y favoreciendo la peroxidación lipídica. A nivel endotelial existe un importante comprometimiento de la peroxidación lipídica, activada por diferentes mecanismos. La disminución de la biodisponibilidad de ON, con una generación ascendente de ERO como O2- y ONOO a este nivel, condiciona una reacción en cadena donde la peroxidación lipídica se incrementa en el tiempo y, junto a ella, las concentraciones de MDA.
En el caso particular del GSH, los estadígrafos obtenidos del análisis brindaron información satisfactoria sobre su elevado poder clasificatorio a través de la función obtenida. El valor negativo y elevado de su coeficiente se corresponde con la disminución de su concentración en el suero de pacientes afectados. Esta condición contribuye a reforzar la teoría de la depresión del sistema antioxidante para neutralizar a algunas especies reactivas de oxígeno, como el radical O2- y H2O2, con un gran comprometimiento en la fisiopatología de la drepanocitemia.
En el análisis, el estadígrafo Lambda de Wilks, que expresa la variabilidad total no debida a la diferencia entre grupos, permitió contrastar la hipótesis nula de que las medias multivariantes de los grupos son iguales, asegurando que la matriz de covarianza es una buena representación de las variables que contribuyen a la discriminación de los grupos. En Cuba, un estudio realizado por Coronados en el servicio de rehabilitación muestra un valor de 0,134 para este estadígrafo como evidencia del poder discriminante de su función, valor similar al obtenido en el estudio. 14) Igualmente, es superada la discriminación entre los grupos en el presente estudio a partir del valor de los centroides, donde el modelo propuesto muestra una apreciable distancia para los grupos definidos a priori.
Por su parte, el ajuste global del modelo quedó confirmado mediante el estadígrafo eta cuadrado(η2), coeficiente de determinación obtenido de la regresión entre la variable dicotómica y las puntuaciones discriminantes. Este mostró un valor elevado de vector propio, indicando cuán eficaz fue el análisis para clasificar a los pacientes y, por tanto, el alto poder discriminante de la función. Quedó confirmado este aspecto, además, por la matriz de clasificación, que reveló un 100% de acierto en los casos clasificados a priori. Además, el elevado valor de la correlación canónica indicó la adecuada pertenencia de los pacientes a los grupos conformados a priori. De esta manera quedó asegurado el ajuste global de la función resultante, dando una muestra del poder clasificatorio de las variables que la conforman.
Debido a la inexistencia de un modelo similar para esta condición, no fue posible establecer comparaciones en cuanto a las variables clasificatorias con la de otras investigaciones. La justificación de los hallazgos se sustentó en la revisión de la bibliografía especializada sobre el tema. Un estudio relacionado a la temática fue el de Malihe y col,15 quienes emplean el análisis discriminante longitudinal para obtener predictores de pre-eclampsia partir de perfiles de hemoglobina por cada trimestre de gestación, obteniendo un mejor resultado en el segundo trimestre.
El ser estos los biomarcadores de elección en la función resultante, denota su potencia estadística como discriminantes. De manera general fue muy importante la consistencia mostrada por el modelo que evidencia su valor diagnóstico.
Índices de riesgo de CVO
Resultó muy favorable para el estudio el resultado mostrado por el modelo en relación a la discriminación de los grupos, hecho que favoreció la estratificación de los estados de riesgo propuestos con un mayor nivel de precisión, posibilitando un diagnóstico evolutivo del riesgo de CVO en tres categorías. Los índices de riesgo de CVO calculados asignaron a cada paciente una clasificación, que en el estudio estuvo acorde con su estado real de salud. Considerando estos resultados, se muestra una excelente relación entre los casos observados y esperados en todos los niveles de riesgo, evidenciando el valor práctico del modelo que se propone.
El hecho de haber obtenido una medida del riesgo hacia la CVO que presentan los pacientes aquejados de anemia drepanocítica, constituye una posible herramienta diagnóstica de valor durante su seguimiento, acompañada de una adecuada valoración clínica. El empleo de la misma en los servicios de hematología contribuiría a la mejora de la atención médica y de la calidad de vida de estos pacientes.
En conclusión, se formuló un modelo predictivo del estado de CVO en la anemia drepanocítica, que incluye biomarcadores del estrés oxidativo, como predictores significativos. El modelo sustentó los índices predictivos propuestos, que determinan y estratifican el riesgo de la CVO en tres categorías. Queda justificado con estos resultados el valor de la evaluación de estos biomarcadores en el paciente drepanocítico, con vista a brindar un mejor seguimiento. Etapas posteriores de la investigación estarán dirigidas a la implementación en el servicio de hematología de ambos aportes, a partir del diseño de un sistema computarizado de apoyo al diagnóstico, que permita materializar estos resultados y minimizar el riesgo y frecuencia de estas complicaciones en los pacientes con el uso oportuno de la terapéutica adecuada.
