Alteración de la entropía en la precuña y la corteza cingulada posterior en la enfermedad de Alzheimer: estudio de resonancia magnética funcional en reposo

Alteration of Entropy in the Precuneus and Posterior Cingulate Cortex in Alzheimer’s Disease: A Resting-State Functional Magnetic Resonance Study

DOI 10.22430/22565337.2118 Tablas Figuras

Recibido: 09 agosto 2021
Aceptado: 23 noviembre 2021
Disponible: 16 noviembre 20211

Cómo citar / How to cite
A. C. Puche; J. F. Ochoa-Gómez; Y. A. Agudelo-Londoño; J. K. Rodas-Marín; C. A. Tobón-Quintero, “Alteración de la entropía en la precuña y la corteza cingulada posterior en la enfermedad de Alzheimer: estudio de resonancia magnética funcional en reposo”, TecnoLógicas, vol. 24, nro. 52, e2118, 2021. https://doi.org/10.22430/22565337.2118
Resumen

El cerebro humano ha sido descrito como un sistema complejo. Su estudio por medio de señales neurofisiológicas ha desvelado la presencia de interacciones lineales y no lineales. En este contexto, se han utilizado métricas de entropía para descubrir el comportamiento cerebral en presencia y ausencia de alteraciones neurológicas. El mapeo de la entropía es de gran interés para el estudio de enfermedades neurodegenerativas progresivas como la enfermedad de Alzheimer. El objetivo de este estudio fue caracterizar la dinámica de las oscilaciones cerebrales en dicha enfermedad por medio de la entropía y la amplitud de las oscilaciones de baja frecuencia a partir de señales Bold de la red por defecto y la red de control ejecutivo en pacientes con Alzheimer e individuos sanos, utilizando una base de datos extraída de la serie de estudios de imágenes de acceso abierto. Los resultados revelaron mayor poder discriminatorio de la entropía por permutaciones en comparación a la amplitud de fluctuación de baja frecuencia y la amplitud fraccional de fluctuaciones de baja frecuencia. Se obtuvo un incremento de la entropía por permutaciones en regiones de la red por defecto y la red de control ejecutivo en pacientes. La corteza cingulada posterior y la precuña manifestaron característica diferencial al evaluar la entropía por permutaciones en ambos grupos. No hubo hallazgos al correlacionar las métricas con las escalas clínicas. Los resultados demostraron que la entropía por permutaciones permite caracterizar la función cerebral en pacientes con Alzheimer, además revela información sobre las interacciones no lineales complementaria a las características obtenidas por medio del cálculo de la amplitud de las oscilaciones de baja frecuencia.

Palabras clave: Resonancia magnética funcional, enfermedad de Alzheimer, entropía por permutaciones, red por defecto, red de control ejecutivo.

Abstract

The human brain has been described as a complex system. Its study using neurophysiological signals has revealed the presence of linear and non-linear interactions. In this context, entropy metrics have been used to discover brain behavior in the presence and absence of neurological alterations. Entropy mapping is of great interest for the study of progressive neurodegenerative diseases such as Alzheimer’s Disease (AD). The objective of this study was to characterize the dynamics of brain oscillations in AD using entropy and the Amplitude of Low-Frequency Fluctuations (ALFF) of BOLD signals from the default network and the executive control network in patients with AD and healthy individuals. For this purpose, the data was extracted from the Open Access Series of Imaging Studies (OASIS). The results revealed greater discriminatory power in Permutation Entropy (PE) than in ALFF and fractional ALFF metrics. An increase in PE was obtained in regions of the default network and the executive control network in patients. The posterior cingulate cortex and the precuneus exhibited a differential characteristic when PE was evaluated in both groups. There were no findings when the metrics were correlated with clinical scales. The results showed that PE can be used to characterize the brain function in patients with AD and reveals information about non-linear interactions complementary to the characteristics obtained by calculating the ALFF.

Keywords: Functional magnetic resonance imaging, Alzheimer’s disease, Permutation entropy, Medical image processing, Default mode network, Executive network.

1. INTRODUCCIÓN

La enfermedad de Alzheimer (EA) es la causa más común de demencia [1]. Según proyecciones realizadas en el 2014 [2], el número de personas con EA en Colombia para el 2020 es de 260 000 con un costo aproximado de 1 billón de pesos [3]. La EA se caracteriza por una acumulación excesiva de placas de β-amiloide extracelular y la presencia de ovillos neurofibrilares de proteína tau hiperfosforilada intracelulares, que da lugar a una configuración atípica que se extiende por toda la corteza cerebral [4]. La EA ha sido descrita por tres estadios: la etapa preclínica, en el que inicia la acumulación atípica de β-amiloide sin manifestaciones clínicas; el deterioro cognitivo leve (DCL), que puede ser de tipo amnésico y no amnésico; y, por último, la etapa de demencia, en la que se identifican alteración de la memoria y falta de independencia del individuo para realizar tareas de la vida diaria. Esta última etapa, a su vez, tiene tres fases: leve, moderada y severa [5].

Distintas modalidades de imágenes han sido ampliamente usadas para el estudio de la EA. En particular, la resonancia magnética funcional es considerada una tecnología que podría tener relevancia clínica debido a su prometedor potencial para la identificación de alteraciones cerebrales causadas por la enfermedad [6]. El estudio de las oscilaciones cerebrales por medio de la resonancia magnética funcional en reposo (rs-fMRI) a través del análisis de la señal dependiente del nivel de oxigenación sanguínea (BOLD, por sus siglas en inglés: Blood-Oxygenation-Level-Dependent) ha demostrado alteraciones en la red por defecto (RD) y la red de control ejecutivo (RCE) en pacientes con EA [7]–[10].

El uso de métricas de conectividad basadas en teoría de grafos y el análisis de componentes independientes ha permitido identificar un deterioro acelerado en sujetos con EA en regiones específicas del cerebro, como es el caso de las áreas involucradas en la red por defecto, mientras que en sujetos con envejecimiento normal se han identificado disrupciones con ocurrencia no uniformemente en todo el cerebro [11]. Por otro lado, el uso de las métricas espectrales tales como la amplitud de oscilaciones de baja frecuencia (ALFF, por sus siglas en inglés: Amplitude of Low Frequency Fluctuations) y su versión fraccional (fALFF, por sus siglas en inglés: fractional ALFF) [12] han mejorado el desempeño de algoritmos de aprendizaje de máquina para la discriminación de sujetos con DCL y EA de individuos cognitivamente normales [13].

A pesar de que la rs-fMRI ha arrojado información relacionada con la dinámica cerebral, aún no es una técnica usualmente utilizada en entornos clínicos debido a factores como la presencia de artefactos en las señales y la ausencia de métricas de referencia (gold standard) para la estimación de la función cerebral [6], [14]. Por lo tanto, algunos esfuerzos se han concentrado en realizar innovaciones metodológicas que permitan comprender la función cerebral a partir del análisis de las interacciones de naturaleza no lineal que caracterizan a las señales neurofisiológicas [15]. El uso de métodos basados en teoría de la información permite capturar interacciones no lineales que pueden ser beneficiosas para el entendimiento de la función cerebral [15], además estudios recientes resaltan que el uso de enfoques no lineales para el estudio de la función cerebral brinda información relevante para el estudio en cerebros sanos y patológicos [16].

Varios enfoques se han utilizado para estimar la complejidad cerebral por medio de métricas no lineales, dentro de ellas se encuentra una variedad de enfoques para el cálculo de la entropía [17]. Un enfoque menos usado es la entropía por permutaciones, un método simple, rápido y robusto para el análisis de series de tiempo caóticas [18]. Esta métrica ha sido utilizada previamente para el análisis de la señal BOLD en sujetos de la iniciativa de Neuroimágenes de la Enfermedad de Alzheimer (ADNI, por sus siglas en inglés) con DCL, EA y controles en donde se encontró una disminución de la entropía en pacientes con EA en comparación con sujetos con DCL y controles. El estudio reportó hallazgos significativos de correlación entre la métrica y las escalas clínicas Mini Mental y valoración clínica de la demencia. El análisis fue realizado con señales BOLD de 130 puntos de tiempo y una dimensión embebida de 4. Los parámetros utilizados son resaltados por los autores como una limitación del estudio [19]. Hasta donde tienen conocimiento los autores, es el único estudio realizado en EA utilizando la métrica de entropía por permutaciones, lo que evidencia que a la fecha existen pocos estudios que hayan utilizado métricas de complejidad para el análisis de datos de fMRI en EA [20].

En este estudio se utilizaron datos de rs-fMRI con 164 puntos de tiempo de la iniciativa de acceso abierto OASIS-3 (en inglés: Open Access Series of Imaging Studies: Longitudinal MRI Data in Non-demented and Demented Older Adults) [21] de sujetos rotulados como cognitivamente normales (CN) y con EA. Se cuantificó la función cerebral por medio de la entropía por permutaciones y se comparó con el ALFF y fALFF para la red por defecto y la red de control ejecutivo; adicionalmente, las métricas fueron correlacionadas con escalas clínicas. Lo anterior, con el objetivo de evaluar si el uso de la entropía por permutaciones brinda información adicional a la obtenida con métricas espectrales, sobre la dinámica de las oscilaciones cerebrales y brinda características diferenciales en pacientes con EA y sujetos sanos al ser evaluadas en las redes cerebrales en reposo que se ven mayormente afectadas en la EA.

2. METODOLOGÍA

2.1 Sujetos

Datos de resonancia magnética estructural y resonancia magnética funcional etiquetados como EA y CN fueron descargados de la base de datos OASIS-3 [21], un proyecto longitudinal en donde a varios individuos se les realizó más de una adquisición. Dos expertos analizaron la información clínica disponible en la plataforma de los sujetos etiquetados como EA y CN [21], teniendo en cuenta información como: edad, antecedentes personales, edad de inicio de los síntomas, cuadro clínico registrado en el tiempo, evolución de las alteraciones cognitivas, puntajes en pruebas neuropsicológicas y escalas de funcionalidad disponibles. A partir de la revisión se excluyeron sujetos etiquetados como EA que durante el seguimiento longitudinal tuvieron un diagnóstico posterior de CN, cuadro clínico leve que sugería un síndrome demencial incierto, evolución de los síntomas que aparentaban una demencia debido a otra condición o enfermedad neurodegenerativa, presentar síntomas depresivos que pudieran explicar las quejas cognitivas, tener la resonancia previa al inicio de los síntomas o tener antecedente de una enfermedad neurológica mayor distinta a EA que pudiera alterar los resultados (como antecedente de enfermedad cerebrovascular o traumatismo encefalocraneano).

Luego de realizar la revisión clínica, se preprocesaron datos de 117 sujetos, de los cuales 23 fueron descartados por exceso de movimiento (EA: 12, CN:11) y 1 sujeto del grupo CN fue descartado por tener valores atípicos en varias regiones de la red por defecto. Finalmente se incluyeron datos de 93 sujetos (EA: 36, CN:57) para los análisis posteriores. Si bien la base de datos no proporciona específicamente información sobre la gravedad de la enfermedad, esta se determinó por la suma de cajas de la escala de valoración clínica de la demencia (CDR-SB, por sus siglas en inglés: Dementia Rating Scale Sum of Boxes) y la escala Mini Mental (MMSE, por sus siglas en inglés: Mini-Mental State Examination). Dado que los individuos del grupo EA están caracterizados por un 4.5 < CDR-SB < 9 pertenecen a la categoría Alzheimer leve [22]. La información demográfica detallada puede observarse en la Tabla 1.

Tabla 1. Información demográfica de los sujetos de estudio. Fuente: elaboración propia.
Parámetro EA Media ± Desviación estándar CN Media ± Desviación estándar Valor p
n 36 57 -
Sexo (femenino: masculino) 16:20 24:33 0.48
Edad 74.70 ± 5.89 70.97 ± 6.52 <0.01
Escolaridad (años) 15.33 ± 2.77 15.68 ± 2.70 0.72
MMSE 25.19 ± 3.14 28.91 ± 1.24 <0.01
CDR-SB 3.33 ± 1.45 0.00 ± 0.00 <0.01

2.2 Adquisición de las imágenes

Las imágenes de rs-fMRI fueron adquiridas con un escáner SIEMENS 3T (MAGNETOM Trio Tim, Siemens, Erlangen, Germany). Se adquirieron 164 volúmenes de 36 cortes cada uno y tamaño de voxel de 4x4x4 mm con los parámetros: tiempo eco (TE) = 27 ms, tiempo de repetición (TR) = 2200 ms, ángulo de giro (FA) = 90°, tamaño de la matriz = 64. Información adicional relacionada con la adquisición de las imágenes se puede encontrar en [23].

2.3 Procesamiento de las imágenes

El preprocesamiento de las imágenes fue realizado en la toolbox CONN [24]. Se hizo una inspección visual para descartar imágenes de baja calidad. Posteriormente se realizó realineación de los datos tomando el primer corte como referencia, corrección temporal y registro con la imagen T1 correspondiente. Luego se inspeccionó el reporte de corrección de movimiento, se excluyeron los sujetos con movimiento de cabeza >2° en rotación y >2 mm en traslación en cualquier dirección. Las imágenes de los 117 sujetos fueron preprocesadas, de los cuales 23 sujetos (EA:12, CN:11) fueron descartados por exceso de movimiento.

Las imágenes de los 94 sujetos restantes fueron normalizadas al espacio estándar definido por el Instituto Neurológico de Montreal (MNI, por sus siglas en inglés: Montreal National Institute), conservando el tamaño de voxel 2x2x2 mm. Se realizó eliminación de tendencias de movimiento, señales relacionadas con la sustancia blanca y líquido cefalorraquídeo. Se aplicó un filtro pasa-banda en el rango de frecuencia de 0.01 Hz a 0.1 Hz.

2.4 Extracción de la señal BOLD

La extracción de las señales se realizó a partir de las plantillas propuestas por [25] para la red por defecto y la red de control ejecutivo (Figura 1). Se utilizaron máscaras de cada región con número de voxeles mayor a 50, por lo tanto, se extrajeron 17 señales BOLD para la red por defecto y 11 señales BOLD para la red de control ejecutivo. El proceso de obtención de las señales fue realizado con la herramienta DPARSFA [26] (DPABI_V4.1_190725). En la Tabla 2 se dan detalles de los nombres de cada región involucrada.

Regiones de interés de la red por defecto y la red de control ejecutivo
Figura 1. Regiones de interés de la red por defecto y la red de control ejecutivo. Fuente: elaboración propia.
Tabla 2. Nombre de las regiones que conforman las redes según
la plantilla de referencia
. Fuente: elaboración propia.
Etiqueta Red por defecto Etiqueta Red de Control Ejecutivo
1 Corteza Prefrontal Medial, Corteza Cingulada Anterior y Corteza Orbitofrontal (B) 1 Giro Frontal Medio, Giro Frontal Superior (I)
2 Giro Angular (I) 2 Giro Frontal Inferior, Giro Orbitario Frontal (I)
3 Giro Frontal Superior (D) 3 Lóbulo Parietal Superior, Lóbulo Parietal Inferior, Precuña, Giro Angular (I)
4 Corteza Cingulada Posterior, Precuña (B) 4 Giro Temporal Inferior, Giro Temporal Medio (I)
5 Corteza Cingulada Anterior Dorsal (B) 5 Crus I (D)
6 Tálamo (B) 6 Giro Frontal Medio, Giro Frontal Superior (D)
7 Hipocampo (I) 7 Giro Frontal Medio (D)
8 Hipocampo (D) 8 Lóbulo Parietal Inferior, Giro Supramarginal, Giro Angular (D)
9 Corteza Retrospenial, Corteza Cingulada Posterior (I) 9 Giro Frontal Superior (D)
10 Giro Frontal Medio (I) 10 Crus I, Crus II, Lóbulo VI (D)
11 Giro Parahipocampal (I) 11 Núcleo Caudado (D)
12 Giro Occipital Medio (I) - -
13 Corteza Retrospenial, Corteza Cingulada Posterior (D) - -
14 Precuña (B) - -
15 Giro Frontal Superior, Giro Frontal Medio (D) - -
16 Giro Parahipocampal (D) - -
17 Giro Angular. Giro Occipital Medio (D) - -
* D: regiones ubicadas en el hemisferio derecho, I: regiones ubicadas en el hemisferio izquierdo y B: regiones localizadas en ambos hemisferios.

2.5 Análisis de la señal BOLD

Para cada región de la RD y la RCE se aplicaron las medidas espectrales: ALFF y fALFF y la métrica no lineal: entropía por permutaciones (PE, por sus siglas en inglés: Permutation Entropy). Todas las métricas fueron calculadas sobre las señales en la banda de frecuencia de 0.01 Hz a 0.1 Hz.

2.5.1 Amplitud de las fluctuaciones de baja frecuencia

La ALFF se define como la suma de amplitudes en una banda de frecuencia, mientras que la fALFF se define como la fracción de la suma de amplitudes en una banda de frecuencia. Para una serie de tiempo de longitud N definida en (1), el cálculo de las métricas ALFF y fALFF está representado por (2) y (3), respectivamente [27].

                     (1)

                    (2)

                    (3)

2.5.2 Entropía por permutaciones

La entropía por permutaciones es un enfoque propuesto por [18] para el análisis de series de tiempo complejas y caóticas, la cual tiene la ventaja de dar resultados significativos en presencia de ruido observacional y dinámico, además, es un método simple, de cálculo rápido y robusto. Para una serie de tiempo {𝑥𝑡} = 1,…,. , 𝜋 permutaciones de orden 𝑛, la PE se define por (4), donde 𝑝(𝜋) está definido por (5), donde el símbolo # hace referencia a “número” [18].

                    (4)

                    (5)

Para el cálculo de la EP se utilizó un retraso embebido (𝑙) de (2) y una dimensión embebida (𝑚) de (4) [28], donde 𝑚 y 𝑙 cumplieron la condición definida por (6) para evitar el submuestreo de las señales.

                    (6)

N = 164, corresponde a la longitud de la serie de tiempo. El parámetro 𝑙 no tiene mucha influencia en la entropía de la serie de tiempo, su valor se tomó con base en la literatura reportada por estudios de electroencefalografía [28]. ¡La medida fue estandarizada por el algoritmo natural de 𝑚!, donde PE alcanza su valor máximo, (7):

                    (7)

2.5.3 Análisis estadístico

Se removió el efecto de la edad de las métricas por medio de una regresión lineal. Para la comparación entre grupos se utilizó una prueba estadística no paramétrica de dos muestras, descrita por [29], con nivel de significancia p<0.05 y se realizó corrección por medio del método (FDR, por sus siglas en inglés: False Discovery Rate, p_fdr). Se analizaron los valores-p sin corrección (p_unc) para comparaciones múltiples por tratarse de un estudio exploratorio sobre la métrica PE y por ser un estudio con tamaño de muestra diferente entre los grupos (CN:57, AD:36). Se calculó el tamaño del efecto con la g de Hedges, utilizando la herramienta The Measures of Effect Size Toolbox [30]; los hallazgos se discuten principalmente con esta medida. Adicionalmente, se realizó un análisis de correlación utilizando el índice de correlación de Spearman entre las escalas de valoración clínica MMSE, CDR-SB y las métricas obtenidas en el análisis de la señal BOLD del grupo EA.

Las métricas se llevaron al rango de 0 a 1 para mejorar la visualización de los resultados por medio de diagramas de cajas.

3. RESULTADOS

3.1 Red por defecto

Se observó un incremento estadísticamente significativo en la entropía del grupo EA en comparación con el grupo CN en la región 10 que involucra a la corteza retrosplenial y la corteza cingulada posterior del hemisferio izquierdo, con tamaño del efecto mediano de valor -0.57. El incremento en la entropía también se pudo observar en la región 4, que involucra la corteza cingulada posterior y precuña, con tamaño del efecto mediano de valor -0.57 (Tabla 3 y Figura 2).

Tabla 3. Valores p y tamaño del efecto para las 17 regiones de la red
por defecto
. Fuente: elaboración propia
ROI ALFF fALFF PE
p1_unc p2_unc ES p1_unc p2_unc ES p1_unc p2_unc ES
1 .40 .60 -.01 .43 .57 .05 .21 .79 .15
2 .29 .71 .08 .09 .91 .16 .08 .92 .35
3 .47 .53 -.03 .72 .28 -.08 .54 .46 .01
4 .47 .53 -.04 .46 .54 -.07 .99 .01 -.57
5 .53 .47 -.05 .77 .23 -.16 .47 .53 .06
6 .88 .12 -.29 .70 .30 -.10 .24 .76 .08
7 .85 .15 -.32 .45 .55 .03 .08 .92 .37
8 .50 .50 -.12 .27 .73 .12 .18 .82 .29
9 .45 .55 -.02 .51 .49 -.01 1.0 .00* -.57
10 .40 .60 -.03 .65 .35 -.12 .59 .41 -.04
11 .86 .14 -.28 .77 .23 -.09 .84 .16 -.19
12 .26 .74 .09 .36 .64 .07 .16 .84 .26
13 .59 .41 -.13 .52 .48 -.04 .81 .19 -.07
14 .29 .71 .09 .80 .20 -.25 .87 .13 -.28
15 .44 .56 -.07 .53 .47 .00 .87 .13 -.23
16 .82 .18 -.25 .83 .17 -.19 .58 .42 .00
17 .37 .63 .02 .59 .41 -.12 .55 .45 .06

p1_unc: valor-p no corregido para la cola CN > AD, p2_unc: valor-p no corregido
para la cola CN < AD, ES: tamaño del efecto, I: izquierdo, D: derecho, *: región con
p_fdr<0.05 (p1_fdr:1.0, p2_fdr:0.04)
.

EP en la red por defecto, diferencia estadísticamente significativa *p≤0.05 en la región que involucra a la corteza retroesplenial izquierda y la corteza del cíngulo posterior, ** sin corrección en la región que involucra la corteza cingulada posterior y precuña
Figura 2. EP en la red por defecto, diferencia estadísticamente significativa *p≤0.05 en la región que involucra a la corteza retroesplenial izquierda y la corteza del cíngulo posterior, ** sin corrección en la región que involucra la corteza cingulada posterior y precuña. Fuente: elaboración propia.

No se observaron diferencias entre grupos para las métricas ALFF y fALFF, el diagrama de cajas puede visualizarse en las Figuras S1 y S2 del apéndice. Los valores obtenidos al correlacionar las métricas con las escalas clínicas CDR-SB y MMSE para las regiones de la red por defecto se encuentran en las Tablas S1 y S2, respectivamente del apéndice.

3.2 Red de control ejecutivo

Se observó aumento del ALFF en el núcleo caudado derecho en el grupo EA, con tamaño del efecto mediano de valor -0.43. Por otro lado, se obtuvo un incremento en la entropía del grupo EA en comparación con el grupo CN en las siguientes regiones: Crus I derecho con tamaño del efecto -0.38, giro frontal medio derecho con tamaño del efecto -0.35, núcleo caudado derecho con tamaño del efecto -0.37 y en la región que involucra al Crus I, Crus II y lóbulo VI derecho con tamaño del efecto -0.38. Las diferencias descritas previamente presentaron significancia estadística sin corrección FDR (Tabla 4, Figura 3 y Figura 4).

No se observaron diferencias entre grupos para la métrica fALFF. El diagrama de cajas puede visualizarse en la Figura S3 del apéndice.

Tabla 4. Valores p y tamaño del efecto para las 11 regiones de la red
de control ejecutivo
. Fuente: elaboración propia
ROI ALFF fALFF PE
p1_unc p2_unc ES p1_unc p2_unc ES p1_unc p2_unc ES
1 .47 .53 .01 .88 .12 -.24 .34 .66 .09
2 .21 .79 .17 .41 .59 .05 .56 .44 -.04
3 .41 .59 .06 .20 .80 .18 .24 .76 .16
4 .85 .15 -.23 .82 .18 -.20 .78 .22 -.17
5 .71 .29 -.12 .58 .42 -.04 .97 .03 -.38
6 .41 .59 .05 .83 .17 -.21 .29 .71 .12
7 .12 .88 .25 .08 .92 .31 .95 .05 -.35
8 .52 .48 -.01 .40 .60 .06 .17 .83 .22
9 .39 .61 .06 .14 .86 .25 .73 .27 -.13
10 .41 .59 .05 .79 .21 -.17 .96 .04 -.38
11 .98 .02 -.43 .13 .87 .25 .96 .04 -.37

p1_unc: valor-p no corregido para la cola CN > AD, p2_unc: valor-p no corregido
para la cola CN < AD, ES: tamaño del efecto, I: izquierdo, D: derecho
.

 

ALFF: en la red de control ejecutivo se encontraron diferencias estadísticamente significativas ** sin corrección en el núcleo caudado derecho
Figura 3. ALFF: en la red de control ejecutivo se encontraron diferencias estadísticamente significativas ** sin corrección en el núcleo caudado derecho. Fuente: elaboración propia.
EP: en la red de control ejecutivo se encontraron diferencias estadísticamente significativas ** sin corrección en la región Crus I derecho, giro frontal medio derecho, la región que involucra al Crus I, Crus II y lóbulo VI derecho, y en el núcleo caudado derecho
Figura 4. EP: en la red de control ejecutivo se encontraron diferencias estadísticamente significativas ** sin corrección en la región Crus I derecho, giro frontal medio derecho, la región que involucra al Crus I, Crus II y lóbulo VI derecho, y en el núcleo caudado derecho. Fuente: elaboración propia.

Los valores obtenidos al correlacionar las métricas con las escalas clínicas CDR-SB y MMSE para las regiones de la red de control ejecutivo se encuentran en las Tablas S1 y S2, respectivamente del apéndice.

4. DISCUSIÓN

En este estudio se investigó el potencial de la métrica entropía por permutaciones para extraer información sobre la complejidad cerebral en regiones que conforman la red por defecto y la red de control ejecutivo, se comparó su desempeño con las métricas espectrales ALFF y fALFF en pacientes con EA y sujetos sanos. Se encontró un aumento de la EP en pacientes con EA en ambas redes, adicionalmente se encontró un aumento en el ALFF en pacientes con EA en la red de control ejecutivo.

Al analizar la red por defecto en pacientes, se halló un incremento de la EP en la región que involucra a la corteza cingulada posterior y la precuña. Por otro lado, el estudio sobre la red de control ejecutivo reveló un incremento de la ALFF en el núcleo caudado y un incremento de la EP en las regiones del hemisferio derecho: giro frontal medio, núcleo caudado y la región que comprende al Crus I, Crus II y el lóbulo VI. En línea con estos resultados se ha reportado que la desorganización en la conectividad funcional de esta región es de importancia para el desarrollo de demencia en la EA [31]. La corteza posterior del cíngulo es considerada una región clave de la red por defecto en la que se han reportado cambios en la conectividad cerebral en EA [32], con alto poder discriminatorio entre esta y sujetos controles [33]. Desde un punto de vista patológico, la precuña y la corteza cingulada posterior se han asociado con el desarrollo completo de la EA, en donde se ha reportado retención de proteína tau y astrogliosis, además de alteraciones significativas en la conectividad funcional [31]. Un estudio estructural de resonancia magnética reportó que la corteza posterior del cíngulo es una de las regiones más vulnerables involucrada en la patofisiología de la EA. Se identificó un deterioro acelerado de múltiples dominios cognitivos en individuos con una mayor tasa de atrofia del volumen de la sustancia gris de esta región [34]. En un metaanálisis de conectividad con rs-fMRI se ha reportado en sujetos con EA sin conectividad entre la red por defecto y el sistema límbico, dado principalmente por las áreas del precuneus y el cíngulo posterior [8]. Otro estudio que evaluaba la conectividad intrared por defecto mostró aumento de la conectividad entre el cíngulo posterior con la corteza parahipocampal y la precuña en personas con EA comparada con controles [35]. Sin embargo, se debe tener en cuenta que los hallazgos de conectividad en reposo en EA son contradictorios. Al parecer, la conectividad en la red por defecto se puede modificar según el estadio de la demencia, pudiendo haber una disminución en áreas posteriores parietales y temporales al inicio de los síntomas y un aumento en áreas anteriores frontales con el progreso de la enfermedad, no se ha esclarecido si es debido a una posible compensación o un cambio debido a la fisiopatología misma del depósito de proteínas anormales [36].

Contrario a los resultados de esta investigación, estudios previos de rs-fMRI [19], [37], [38] han reportado reducción de la complejidad en EA en comparación con sujetos sanos. Uno de esos estudios [19] halló disminución de la EP en el grupo EA al comparar con el grupo control en regiones de los lóbulos temporal, occipital y frontal (hemisferio derecho: giro temporal inferior, cuña, giro occipital medio y giro occipital superior; hemisferio izquierdo: giro frontal medio, giro frontal superior, giro cingulado anterior y cuña). Por otro lado, el mismo estudio reportó incremento de la EP en el grupo de deterioro cognitivo leve de etapa temprana al comparar con el grupo control en las regiones del hemisferio derecho: giro temporal inferior, giro frontal medio, cuña, giro occipital medio y giro occipital superior. Es importante resaltar que el estudio referenciado no partió del estudio de regiones que forman parte de las redes cerebrales en reposo que se ven mayormente afectadas en la enfermedad de Alzheimer; en su lugar, el análisis fue realizado sobre 8 regiones pertenecientes a 5 clústeres que mostraron diferencias significativas luego de obtener la métrica EP sobre todos los voxeles del cerebro. Los datos del estudio fueron tomados de la iniciativa ADNI. Por otro lado, el estudio de [37] utilizó la métrica entropía multiescala para cuantificar la complejidad de las señales BOLD en sujetos con deterioro cognitivo leve temprano, deterioro cognitivo leve tardío y sujetos controles tomados de la iniciativa ADNI. El estudio reportó disminución de la complejidad en la región del giro fusiforme izquierdo y en la corteza cingulada anterior rostral en el grupo de deterioro cognitivo leve temprano. El enfoque metodológico consistió en la generación de mapas de entropía multiescala para cada sujeto sobre todos los voxeles del cerebro. Con respecto al estudio de [38], estos utilizaron la métrica entropía multiescala para cuantificar la complejidad de las señales BOLD en sujetos con deterioro cognitivo leve y sujetos controles reclutados en la Clínica de Memoria del Departamento de Geriatría del Hospital Universitario Karolinska, en Huddinge, Suecia. El análisis se realizó sobre regiones de la red por defecto definida por [25], de las cuales reportó disminución de la entropía a nivel global y a nivel nodal reportó reducción en el hipocampo derecho. Los resultados se reportan sin corrección para comparaciones múltiples. Como limitaciones comunes de los estudios citados previamente se resalta que la iniciativa ADNI no proporcionó públicamente información sobre los factores de riesgo relacionados con la EA [19] y se indica que estudios previos han sugerido que el criterio diagnóstico de deterioro cognitivo leve de la iniciativa ADNI tiene alta tasa de falsos positivos [37].

El estudio de [39] evaluó el impacto del envejecimiento en la demencia desde un enfoque celular, en el que la EA de aparición tardía puede verse como una manifestación de producción de energía reducida, dando lugar a un aumento en la entropía. Los autores sugieren que la acumulación de modificaciones relacionadas con la edad, así como el agotamiento de las mitocondrias, podrían ser responsables de la escasa capacidad del cerebro para adaptar sus estructuras y funciones cerebrales. Por otro lado, la relación entre la entropía durante el envejecimiento saludable y la EA ha sido explorada recientemente por medio de rs-fMRI y la métrica de mapeo de entropía cerebral (BEN) [40] en grupos de individuos controles (n=54, edad: 65-95 años), controles con compromiso de la memoria (n=27, edad 65-95 años), deterioro cognitivo leve temprano (n=58, edad: 56-89 años), deterioro cognitivo leve tardío (n=38, 57-88 años) y sujetos con EA (n=34, 56-87 años). El estudio referenciado reportó el siguiente comportamiento de la BEN en los cuatro grupos de interés. La comparación se describe con respecto al grupo control: leve aumento en los individuos controles con compromiso de la memoria, incremento más notorio en el grupo de deterioro cognitivo leve temprano. Para la etapa de deterioro cognitivo leve tardío se reportó una disminución leve que posteriormente se convirtió en una reducción más notoria en el grupo de EA, la reducción en la BEN en las etapas tardías de la demencia es relacionadas por el autor con un fenómeno fallido de compensación [40]. Según [41] y [42], la reducción notoria de la BEN presente en las etapas tardías de la demencia da lugar al deterioro acelerado de funciones cerebrales que requieren niveles de entropía altos para ser flexibles funcionalmente. Los hallazgos de este estudio, descritos como demencia tipo Alzheimer en la sección 2.1, están en línea con el comportamiento descrito por el grupo de deterioro cognitivo leve de inicio temprano en el estudio realizado por [40], resalta este aspecto por ser el único estudio que describe el comportamiento de la entropía en 5 estadios cognitivos de interés en el estudio de la EA. Sin embargo, a través de la referencia [38] se evidencian contradicciones en los hallazgos de estudios realizados a partir de datos de rs-fMRI en EA con el objetivo de cuantificar la complejidad de las señales BOLD en EA, lo que sugiere que es necesario realizar estudios adicionales en el campo, que además incluyan seguimiento longitudinal y proporcionen información adicional sobre los criterios de selección de sujetos.

La ALFF es una métrica que permite identificar la actividad neuronal espontánea de regiones específicas del cerebro a través de la intensidad de la señal BOLD. Esta ha sido utilizada como métrica característica para mejorar el desempeño de algoritmos de aprendizaje de máquina para la discriminación de sujetos con DCL y EA de individuos cognitivamente normales [13]. La ALFF y las métricas de entropía no están relacionadas teóricamente, sin embargo, en este estudio se encontró un aumento en la ALFF y la PE en el núcleo caudado. El estudio de [43] realizado con una cohorte de sujetos sanos reportó correlación positiva de magnitud mediana a alta entre el BEN y la fALFF en la corteza orbitofrontal y la corteza temporal inferior posterior. Por otro lado, reportó correlaciones negativas de gran magnitud entre el BEN y la fALFF en la corteza visual, corteza temporal inferior anterior, red motora, precuña y la corteza parietal lateral, lo que sugiere que la fALFF y la BEN son mutuamente independientes. El estudio referenciado previamente sugiere que la BEN puede proveer información más completa de la función cerebral, especialmente en regiones donde no hay asociaciones entre la BEN y la fALFF, en línea con los hallazgos de nuestro estudio.

Dentro de las limitaciones de esta investigación, se encuentra el uso de grupos no pareados por sexo y edad. Es importante destacar que OASIS-3 es una base de datos que incluye sujetos con edades comprendidas entre los 42 años y los 95 años que fueron reclutados en el transcurso de 30 años. Algunos de los sujetos incluidos en este estudio fueron inicialmente etiquetados como cognitivamente normales que desarrollaron demencia por EA con el tiempo y otros ingresaron a la iniciativa OASIS-3 con diagnóstico de demencia tipo Alzheimer, para ambos casos la revisión realizada por los expertos clínicos permitió excluir sujetos con etiqueta EA que podían tener síntomas explicados por otras condiciones clínicas o no contar con un diagnóstico claro, a partir de la información proporcionada por la historia clínica disponible en la base de datos. Se debe tener en cuenta que en la base de datos no se especifica el estadio clínico de la demencia en la que se encuentran los sujetos, por lo que se tuvo que determinar a partir de mediciones de escalas que, si bien son fiables, no igualan el criterio clínico. Los enfoques de análisis reportados por estudios previos limitan la comparación de los hallazgos de este estudio, principalmente por la incertidumbre acerca del estadio de los sujetos clasificados como EA o deterioro cognitivo leve y tardío de la iniciativa ADNI, además del uso de métricas diferentes a la EP para cuantificar la complejidad de las señales BOLD. Además, se encontraron algunas limitaciones relacionadas con la presencia de artefactos en las imágenes y cambios en el diagnóstico durante el seguimiento longitudinal que disminuyó la muestra inicial. Aunque el tamaño de las señales es de mayor duración en comparación al estudio de PE realizado por [19], la realización de este análisis sobre señales con más de 164 volúmenes sería beneficioso para evaluar la entropía en el continuum de la EA debido a los 24 estados definidos por la dimensión embebida igual a 4.

Si bien las pruebas estadísticas dan una inferencia considerable de la relación entre dos o más variables, permitiendo una interpretabilidad dentro de los datos, su precisión predictiva no es su principal fortaleza, por lo que en investigaciones futuras se deberían considerar los efectos potenciales del aprendizaje de máquina para complementar los resultados presentados, donde se tendría la capacidad de generar un modelo a partir de la métrica no lineal PE y obtener predicciones reproducibles, con un análisis más robusto en la búsqueda de patrones discriminantes entre grupos, en este caso EA vs CN, que permitiría, además, obtener un mayor poder predictivo con la adición de mayor cantidad de muestras de las diferentes poblaciones y reducir la presencia de artefactos [44]. En consecuencia, de acuerdo a ciertos estudios, [16], [45], se ha sugerido el uso de métricas de conectividad funcional de fluctuaciones dinámicas no lineales como datos de base para el análisis, ya que proporciona una mejor caracterización patofisiológica de las redes neuronales en poblaciones con enfermedades neurodegenerativas, además de ser datos que pueden ser eficientemente manipulados por las redes neuronales artificiales, puesto que permiten el modelamiento de relaciones no lineales en campos de gran complejidad.

5. CONCLUSIONES

En el estudio actual se identificó incremento de la entropía en regiones de la red por defecto y la red de control ejecutivo en pacientes con EA. El hallazgo diferencial de mayor magnitud se manifestó en la corteza cingulada posterior y la precuña. La metodología implementada permitió demostrar la capacidad de la métrica entropía por permutaciones para capturar información sobre la función cerebral en la enfermedad de Alzheimer en comparación a las métricas espectrales. La métrica debe ser explorada en estudios longitudinales y en señales de mayor duración para evaluar su potencial para capturar patrones funcionales relacionados con la EA.

6. AGRADECIMIENTOS

Este trabajo fue apoyado por el Departamento Administrativo de Ciencias, Tecnología e Innovación (Colciencias), por el proyecto ‘Identificación de Biomarcadores Preclínicos en Enfermedad de Alzheimer a través de un Seguimiento Longitudinal de la Actividad Eléctrica Cerebral en Poblaciones con Riesgo Genético’, código 111577757635.

Los datos fueron proporcionados por OASIS-3: Principal Investigators: T. Benzinger, D. Marcus, J. Morris; NIH P50 AG00561, P30 NS09857781, P01 AG026276, P01 AG003991, R01 AG043434, UL1 TR000448, R01 EB009352.

Parte de la financiación de los autores fue aportada por el fortalecimiento de programas de ciencias médicas y de la salud con talento joven e impacto regional del Ministerio de Ciencia, Tecnología e Innovación (MinCiencias) – RC 752 – 2018. Estos recursos fueron administrados por la Corporación Académica para el Estudio de Patologías Tropicales de la Universidad de Antioquia

CONFLICTO DE INTERÉS

El proyecto no presenta ningún conflicto de interés.

CONTRIBUCIÓN DE LOS AUTORES

Aura C. Puche: diseño del experimento, procesamiento de datos, análisis de resultados y preparación inicial del manuscrito. Revisión, edición y aprobación del documento final.

John Fredy Ochoa Gómez: diseño del experimento y análisis de resultados; Revisión, edición y aprobación del documento final.

Yésika Alexandra Agudelo L.: análisis de resultados y preparación inicial del manuscrito. Revisión, edición y aprobación del documento final.

Jan Karlo Rodas Marín (experto clínico): selección de sujetos y análisis de resultados. Revisión, edición y aprobación del documento final.

Carlos Tobón Quintero (experto clínico): selección de sujetos y análisis de resultados. Revisión, edición y aprobación del documento final

7. REFERENCIAS

APÉNDICE

ALFF en la red por defecto, sin hallazgos estadísticamente significativos.
Figura S1. ALFF en la red por defecto, sin hallazgos estadísticamente significativos. Fuente: elaboración propia.

fALFF en la red por defecto, sin hallazgos estadísticamente significativos
Figura S2. fALFF en la red por defecto, sin hallazgos estadísticamente significativos. Fuente: elaboración propia.

fALFF en la red de control ejecutivo, sin hallazgos estadísticamente significativos.
Figura S3. fALFF en la red de control ejecutivo, sin hallazgos estadísticamente
significativos
. Fuente: elaboración propia.
Tabla S1. Correlación de las métricas con la escala CDR-SB para las regiones
de la red por defecto
. Fuente: elaboración propia.
ROI ALFF Falff PE
r P r P r p
1 .01 .94 -.13 .45 -.14 .41
2 -.10 .55 -.27 .11 -.04 .83
3 .07 .66 -.12 .48 .04 .82
4 .06 .74 .05 .79 .16 .34
5 .10 .58 .13 .46 -.17 .31
6 .05 .76 -.05 .79 .02 .92
7 .05 .77 .02 .90 -.13 .45
8 .06 .72 -.01 .96 .24 .17
9 .07 .70 -.01 .95 .06 .75
10 .13 .44 -.01 .94 .01 .97
11 .09 .59 .25 .15 -.16 .34
12 -.06 .73 .01 .94 .12 .49
13 .11 .54 .08 .66 .03 .86
14 .10 .55 .03 .84 .15 .38
15 .04 .80 .06 .71 -.06 .73
16 .05 .77 -.03 .87 -.20 .24
17 .05 .76 .27 .12 .21 .21
Tabla S2. Correlación de las métricas con la escala MMSE para las regiones
de la red por defecto.
Fuente: elaboración propia.
ROI ALFF fALFF PE
r P r P r p
1 .11 .53 .12 .48 .07 .68
2 .07 .68 .03 .88 .12 .47
3 .17 .32 -.06 .74 .07 .67
4 .15 .38 -.05 .78 -.03 .87
5 .03 .85 -.07 .68 .11 .54
6 .13 .44 -.02 .90 .11 .54
7 -.14 .43 .08 .66 .13 .46
8 -.08 .65 -.12 .50 -.31 .07
9 .11 .54 -.06 .74 .07 .68
10 .09 .61 .11 .51 .32 .06
11 -.13 .46 .05 .75 .14 .41
12 .05 .79 .07 .68 .07 .66
13 .04 .82 .05 .78 .09 .62
14 .13 .44 .02 .93 .02 .92
15 .12 .47 -.05 .75 .07 .67
16 -.08 .65 .05 .79 .47 .00
17 .02 .93 .00 .98 -.13 .44
Tabla S3.Correlación de las métricas con la escala CDR-SB para las regiones
de la red de control ejecutivo
. Fuente :elaboración propia.
ROI ALFF fALFF PE
r P R p r p
1 .16 .35 -.16 .35 -.18 .30
2 .08 .62 -.02 .89 -.28 .10
3 .01 .98 -.05 .79 .13 .45
4 .02 .92 -.01 .95 -.13 .46
5 .16 .35 .07 .68 .07 .71
6 .07 .69 -.15 .38 -.09 .60
7 -.02 .91 .03 .85 .03 .84
8 .03 .86 .13 .44 .02 .90
9 .01 .93 .10 .55 -.03 .85
10 .21 .22 -.06 .73 .31 .06
11 .04 .83 .34 .04 .16 .36

Tabla S4. Correlación de las métricas con la escala MMSE para las regiones
de la red de control ejecutivo
. Fuente: Fuente: elaboración propia.
ROI ALFF fALFF PE
r p R p r p
1 .14 .43 .07 .67 .19 .27
2 .04 .82 .21 .22 .15 .37
3 .04 .83 -.04 .83 .17 .33
4 -.18 .28 -.05 .76 -.04 .80
5 -.10 .54 .02 .93 -.16 .36
6 .11 .52 -.05 .76 .07 .67
7 .09 .62 -.12 .49 .27 .10
8 .02 .91 -.20 .24 .02 .92
9 .14 .42 -.08 .65 .01 .97
10 -.05 .78 .11 .51 .11 .53
11 .03 .84 -.25 .14 .04 .80