Introducción

Las diversas enfermedades cerebrales adquiridas representan un problema de salud pública de alto costo [¹]. Por ello existe un interés en la búsqueda de evidencias que describan las variaciones de la conectividad funcional entre diferentes regiones del cerebro [²], es decir, se trata de estimar la posible relación (a través de medidas de conectividad) entre procesos neurofisiológicos espacialmente distribuidos en redes o circuitos glioneuronales [³]. Estimar la conectividad funcional cerebral, en condiciones de salud o en presencia de padecimiento cerebral, puede apoyar la hipótesis de que las funciones del cerebro están basadas en interacciones entre grupos de neuronas localizadas en diferentes áreas del cerebro, es decir, los patrones de conectividad reflejan los potenciales de campo que generan oscilaciones y que pueden relacionarse con las respuestas sensoriales, motoras y cognitivas. Las enfermedades cerebrales adquiridas deben afectar el funcionamiento de los circuitos cerebrales alterando las conexiones funcionales en las redes y consecuentemente implica un reto en la investigación clínica [¹,²].

Un elevado nivel de conectividad entre la actividad cerebral generada en diferentes ensambles glioneuronales de diversas regiones cerebrales, posiblemente indica una relación funcional entre dichas áreas. Se ha propuesto que de éstos cambios de la conectividad funcional resultan los procesos de sensaciones, percepciones, integraciones semánticas, memoria, funciones ejecutivas y en general los procesos cognitivos. Por lo tanto, las posibles alteraciones de la conectividad funcional pueden afectar dichas funciones cerebrales que involucran áreas especializadas de integración de la información [⁴]. La pregunta que se plantea es: ¿Cómo poder observar cambios en la conectividad funcional y de ser posible evaluarlos computacionalmente?

Tratar de entender un sistema complejo como el cerebro, requiere estudiar su evolución, estructura anatómica y funcional, y el ambiente en el que funciona, por lo tanto se requiere analizar cómo esos múltiples factores o variables interactúan entre sí, generando patrones dinámicos de conectividad. Éstos patrones de conectividad pueden representarse a partir de interacciones de dos simples componentes: nodos y enlaces. Con ésta abstracción se pueden generar redes de conectividad y a su vez estudiar el comportamiento de los elementos de la red en función de los cambios de la dinámica cerebral, de manera global o en grupos locales. Diversas aportaciones [⁵- ¹⁰] han encontrado parámetros propios de una red compleja que caracteriza diversos patrones de conectividad cerebral. Es posible estudiar el comportamiento funcional del cerebro humano aplicando las características de la red de mundo-pequeño, modelando la red de interconexiones de gran escala, de redes locales y globales [⁶]. Su aplicación empieza a dar resultados importantes en el conocimiento del funcionamiento normal, por ejemplo, en el estudio del desarrollo del cerebro [¹¹], de los circuitos que intervienen en el estado de sueño [¹²], del movimiento de las manos o de los pies o en la imaginación del mismo [¹³-¹⁴], y alterado por lesiones cerebrales [⁸,¹⁵], por epilepsia [¹⁶- ¹⁸], en Parkinson [¹⁹], Alzheimer [⁷,²⁰-²¹] y esquizofrenia [²²-²³].

Una gráfica es una representación de una red de objetos, a partir de nodos (vértices) que se vinculan entre sí mediante conexiones (aristas). Dicha conexión puede ser de tipo binaria (conectado, no conectado) o pesada, de acuerdo a un valor estimado de intensidad de la conexión. El EEG puede analizarse a partir de una gráfica que tiene por nodos cada uno de los electrodos del montaje, aristas que conectan a todos los electrodos entre sí (gráfica completa) y con pesos derivados de una medida de estimación de conectividad, tratando de cuantificar alguna relación entre las señales de los electrodos.

En este trabajo se describen técnicas de procesamiento de señales aplicados a registros del electroencefalograma (EEG) para el estudio de la conectividad funcional cerebral desde el punto de vista de la teoría de gráficas. Las técnicas, en conjunto con el enfoque de la teoría de gráficas, brindan diversas alternativas para poder representar las conexiones funcionales de los nodos y enlaces de circuitos glioneuronales cuando el sujeto está en reposo con los ojos cerrados y los cambios que ocurren por estimulaciones sensoriales y sus respuestas.

La electroencefalografía es un método de registro de datos electrofisiológicos no invasivo, fácil de realizar, repetible y sin ninguna contraindicación médica. Medir aún indirectamente las variaciones de actividad eléctrica de ensambles neuronales puede ayudar a explorar el mecanismo de participación y cooperación entre las diferentes áreas del cerebro, y establecer una relación funcional a través de una estimación de medidas dinámicas de conectividad funcional temporal.

El estudio del EEG registrado bajo condiciones de vigilia con ojos cerrados, se considera el estado basal de actividad cerebral [¹]. En cuando menos cuatro ritmos o bandas espectrales altamente estudiados, al presentarse un estímulo se activan y desactivan circuitos neuronales [²⁴], modificándose las relaciones de conectividad funcional entre diversos ensambles neuronales, lo que se manifiesta como disminución de la potencia absoluta de esos ritmos, en particular del ritmo alfa que se ha descrito como desincronización o bloqueo del ritmo alfa. Sin embargo, al presentarse la misma estimulación, sin cambios importantes para el sujeto (estímulo no significativo), la respuesta, pasa de desincronización a sincronización, éste fenómeno se denomina habituación [¹,²⁴]. Se propone que estos procesos resultan de los cambios de conexiones funcionales en los diferentes circuitos glioneuronales.

La teoría de gráficas no había sido incorporada de manera relevante al estudio de la conectividad funcional cerebral antes de la concepción de un modelo de red de pequeño-mundo, desarrollado por Watts y Strogatz [²⁵], y que describe el famoso experimento social de Milgram [²⁶]. Ésta red tiene dos características importantes: muchas conexiones de corta distancia y algunas conexiones de larga distancia, entendiendo por distancia a la cuenta de nodos que hay que atravesar para ir de un nodo a cualquier otro. Éstas dos características de la red de pequeño-mundo, considerando que el comportamiento del cerebro humano se puede modelar a través de una red de interconexiones de gran escala, pueden reflejar información funcional y han servido para encontrar diversos resultados en base al estudio de la conectividad funcional a través de diversas técnicas de adquisición de datos neurofisiológicos y diversas acciones o enfermedades cerebrales como los mencionados anteriormente.

El coeficiente de agrupamiento C de un nodo indica la probabilidad de que sus nodos vecinos estén conectados entre sí, es decir, evalúa el grado de interconectividad entre un grupo de nodos vecinos alrededor de un nodo central [⁵]. Este coeficiente se considera una medida local de interconexión de la red, se le relaciona con la resistencia y adaptación a fallas de la red: la pérdida de nodos en un conjunto local afecta de menor manera si todos los vecinos están conectados entre sí, es decir, el flujo de información no se detendrá encontrando caminos locales alternos [⁶]. El promedio de C de todos los nodos conforma a C-, indicando si la red global tiene grupos locales altamente conectados o no.

La distancia mínima entre dos nodos es la menor cantidad de nodos intermedios necesarios para que el nodo i y el nodo j establezcan un enlace de comunicación. Se puede asegurar que, no importando que se tenga una gran cantidad de nodos, existirá un camino corto entre cualquier par de nodos [⁹]. El promedio de las distancias mínimas de todos los pares de nodos resulta en la distancia mínima característica L-, que es una medida global que indica la cantidad de pasos intermedios promedio con los que se comunican los nodos dentro de la red. Una red con una L- alta indica que la información sufre retardos para llegar de un nodo a otro mientras que una L- baja indica que los nodos pueden comunicarse en pocos pasos, por lo tanto la información sufre menos retardos en su propagación [⁶]. La eficiencia de una red E, es un parámetro derivado de L-, e indica la capacidad de tráfico de la red. Ésta medida al estar relacionada con L- nos indica que tan eficiente es la red en la transmisión de la información, es decir, una red altamente conectada con L- baja puede no ser eficiente al momento de transferir la información ya que necesita pasar por una gran cantidad de nodos antes de llegar a su destino, por el contrario una red con una L- mayor pero no altamente conectada conserva conexiones de baja distancia a pesar que dos nodos se encuentren lejanos, es posible que dicha red sea eficiente para transmitir la información.

Metodología

El análisis de la conectividad funcional a través de la teoría de gráficas brinda un marco para analizar dichas conexiones cerebrales como una red compleja a través de seis pasos cuyo flujo se muestra en la Figura 1.

Figura 1 Diagrama de flujo de la metodología propuesta para el análisis de conectividad con teoría de gráficas. 

Registros de EEG

Los registros de EEG utilizados fueron realizados en el Laboratorio de EEG de la Dirección de Investigación del Instituto Nacional de Rehabilitación (INR); dichos registros fueron obtenidos de 50 estudiantes de licenciatura en Terapia Física o en Comunicación Humana, con una edad promedio de 20.3 ± 2.15 años, mujeres con un nivel socioeconómico y cultural similar, sin afecciones o padecimientos neurológicos, de acuerdo a las historias clínicas y evaluaciones psicológicas de tamizaje, que fueron abiertamente invitados (los sujetos se estratificaron según se describirá en la sección Análisis de datos). Se les explicaron los objetivos y procedimientos del trabajo; los que aceptaron firmaron su consentimiento. El protocolo fue aprobado por los Comités de Ética y de Investigación del INR, de acuerdo con las normas internacionales.

El arreglo de 16 electrodos utilizado fue de tipo monopolar con una distribución 10/20 tomando como referencia el hueso mastoides de cada lado (amplificador v44, Natus, WI, EUA). Durante la adquisición se realizó un filtrado rechaza banda a 60 Hz para eliminar el artefacto provocado por la línea de alimentación. Para la obtención de los registros de EEG los sujetos se colocaron en posición supina durante todo el tiempo del registro, se pidió que permanecieran con ojos cerrados (2-3 minutos), alternaron de 6 a 10 repeticiones de ojos abiertos/cerrados (periodos de 10-15 segundos en cada condición). Estando con los ojos cerrados, recibieron un estímulo visual (fotoestimulación) consistente en 20 series de estímulos lumínicos (luz blanca, 182.4 lúmenes, NicLed Photic Stimulator) con una frecuencia de 5Hz durante 2s y con intervalos intermedios de 20-25 segundos. Otras etapas del registro, que no son de interés para este trabajo, fueron una etapa de asociación con un procedimiento similar al de fotoestimulación, pero con la indicación de presionar un interruptor (normalmente abierto) ante cada fotoestimulación percibida y finalmente un periodo de hiperventilación (3 minutos); en la Figura 2 se puede observar un esquema del protocolo utilizado.

Figura 2 Esquema del protocolo de adquisición de registros de EEG. 

Estimación de conectividad

La estimación de la conectividad se realizó en el dominio de la frecuencia w utilizando la coherencia espectral definida en la Ecuación 1.

(1)

donde S_XXXX (w) y S_YY (w) son las densidades de los autoespectros de las series de tiempo X y Y, éstas series representan la información de dos electrodos del EEG, S_XY (w) es la densidad del espectro cruzado de ambas señales. Para calcular las densidades de los espectros de potencia se usó el método de Welch [³⁶]. La función de coherencia se basa en que la información de los espectros de potencia de las señales, puede describir un grado de sincronización entre la actividad cerebral en el dominio de la frecuencia entre diferentes regiones cerebrales. Una alta coherencia indica que la relación entre las señales en dichas áreas puede ser aproximada por una transformación lineal; esto no implica que la dinámica neocortical sea lineal [²⁷]. Las bandas de frecuencia que se utilizaron para este trabajo fueron las bandas ampliamente estudiadas en el EEG: delta (1.5−4 Hz), theta (4.5−8 Hz), alfa (8.5−13 Hz) y beta (13.5 − 30 Hz).

Las estimaciones de coherencia se obtuvieron por pares de electrodos, obteniendo 120 posibles combinaciones de 16 electrodos del montaje, el valor de coherencia se representa como peso de conexión entre dos electrodos; al calcular todos los posibles pesos se puede representar esta información en una gráfica de conectividad funcional. La gráfica de conectividad tendrá N nodos (electrodos) y N(N − 1) conexiones posibles entre todos los electrodos.

De acuerdo a la literatura de las múltiples conexiones del lóbulo frontal y su participación en la regulación de funciones cognitivas se asume que la conectividad funcional relevante en el proceso de habituación se establece entre los electrodos Fronto-Occipitales (F-O), Fronto-Temporales (F-T), Temporo-Occipitales (T-O), Fronto- Centrales (F-C), Parieto-Occipital (P-O) y Centro-Occipitales (C-O) [³⁷]. Debido a que se ha propuesto que la actividad frontal, activada por el sistema retículo-talámico, actúa regulando las áreas sensoriales como la visual (occipital), de asociación (temporales) y motoras (centrales) [³⁸,³⁹]. Para cada registro se detectó la actividad eléctrica generada en los dos segundos previos (Pre) y dos segundos durante (Rph) en las 20 series de fotoestimulaciones y se procedió a obtener una medida de conectividad en cada serie de estimulación fijando la atención en dichos pares de electrodos, obteniendo el comportamiento promedio de la población de estudio para cada serie en Pre y Rph.

Parámetros de redes complejas

Las redes complejas pueden ser de dos tipos: pesadas y binarias. Las redes binarias son ampliamente estudiadas ya que sus valores de conectividad son 1 (conexión) y 0 (sin conexión). Para las redes pesadas se tienen dos opciones para la estimación de sus parámetros: utilizar parámetros binarios o sus versiones pesadas. En este trabajo se decidió estimar parámetros binarios, debido a que en la literatura los parámetros pesados C- y L- se comparan con una versión de C- y L- de modelos aleatorios.

Un patrón de conectividad (binario) en un instante de tiempo i y los valores de los parámetros de red calculados, al mismo instante de tiempo, dependen del valor de umbralización aplicado sobre los pesos de conectividad de la red, es decir, al utilizar otro valor de umbral se puede obtener una gráfica de conectividad y una distribución de conexiones diferentes (ver Figura 3). Elegir un valor de conectividad y observar el comportamiento de las redes generadas podría no reflejar los posibles cambios en la estructura de las redes, por ello observar a diversos valores los diferentes patrones de conectividad brinda la oportunidad de estudiar si existe un intervalo de valores para los cuales existan cambios en los patrones de conectividad [⁴⁰], y se reflejen en la dinámica de parámetros de redes complejas; aún y cuando es posible que la variación de umbrales provoque inestabilidades locales de las medidas de conectividad [⁴¹], la utilización posterior de estimadores promedio reduce el sesgo introducido por la arbitrariedad de la elección de umbral [⁴²-⁴³]. Es posible que ésta elección pueda brindar más información acerca de los procesos de interconexión en diferentes áreas del cerebro durante el proceso de habituación, es decir, probablemente la presencia de un valor de conectividad alto no necesariamente sea mejor para realizar una tarea o proceso, posiblemente un valor bajo o intermedio de coherencia sea requerido para la función deseada, por ello decidimos no elegir un valor de conectividad único, ¿se puede, al variar el umbral de conectividad, estudiar cómo se comportan los patrones de conectividad a través de parámetros de redes complejas? Los parámetros como el grado de nodo, C- y L-, ¿pueden brindar información al variar el umbral de conectividad, que pueda establecer un intervalo de umbral en el que se detectan cambios? Esto nos lleva a pensar en la variabilidad de la densidad de conexiones como un factor significante en la organización del cerebro, a través de la generación de redes en las cuales varía la densidad de conexión [⁴⁴]. Los parámetros serán estimados y comparados en su comportamiento a partir de un barrido en el intervalo de conectividad, usando 101 valores. La elección de estudiar en pasos de 0.01 del valor de la coherencia tiene interés observar el comportamiento de cada parámetro en pequeños intervalos, es decir, es posible que al elegir un paso mayor de estudio (0.1,0.3) no se puedan observar cambios en la dinámica que arroje cada parámetro estimado, esta es una propuesta exploratoria de cómo abordar el problema de la umbralización y la dinámica de los parámetros de redes en función de su variabilidad.

Figura 3 Gráficas de conectividad generadas a partir de cuatro valores de umbral diferentes. 

Los parámetros de red que se utilizan en este trabajo se describen a continuación. El grado de un nodo es el número de conexiones que tiene un nodo con los demás presentes en la gráfica (Ecuación 2), donde a(i,j) es el valor de la conexión entre los nodos i y j.

(2)

El grado promedio de una gráfica es el promedio de los grados de los nodos de la gráfica (d(G)). El coeficiente de agrupamiento del nodo cj evalúa el grado de interconectividad entre un grupo de nodos vecinos alrededor de un nodo central (Ecuación 3).

(3)

El promedio de los coeficientes de agrupamiento de una gráfica se conoce como coeficiente de agrupamiento característico (C (G)). El promedio de las distancias mínimas de todos los pares de nodos se conoce como distancia mínima característica y se calcula con la Ecuación 4.

(4)

Para concentrar la información de los parámetros de red se propone una representación gráfica en función de los fotoestímulos y del nivel de umbral al que se calculan los parámetros, mostrando el valor obtenido a través de una paleta de colores. En la Figura 4, se muestra la representación para el grado de cada nodo, el número de conexiones de cada electrodo se indica mediante una paleta de colores. El comportamiento en Rph, para cada fotoestímulo, se indica mediante el símbolo ▲ y la etapa Pre está la izquierda de dicho símbolo, el comportamiento de cada fotoestímulo se divide por una barra azul y una blanca. En la Figura 5, se muestra el grado promedio, en función de cada fotoestímulo, en etapas Pre y Rph (Rph se indica de la misma forma que en la Figura 4), por cada valor de umbral se obtiene una red y ésta a su vez tiene un grado promedio, ese valor se indica con una paleta de colores, los resultados de los parámetros C- y L- se presentan en una gráfica similar.

Figura 4 Representación propuesta para observar la variación del grado de cada nodo en función de los fotoestímulos, indicando el número de conexiones de cada electrodo mediante una paleta de colores. La etapa Rph se indica mediante el símbolo ▲ y la etapa Pre está la izquierda de dicho símbolo, el comportamiento de cada fotoestímulo se divide por una barra azul y una blanca. 

Figura 5 Representación propuesta para observar la variación del grado promedio en función de los fotoestímulos, indicando el valor del grado promedio para cada umbral mediante una paleta de colores. La etapa Rph y Pre se indica de la misma forma que la Figura 4. 

Resultados y discusión

Con el enfoque del estudio de la conectividad funcional a través de la teoría de gráficas, es posible observar los diversos patrones de conectividad temporal (cambios dinámicos en los circuitos cerebrales) que se generan en los sujetos en el estado de reposo con los ojos cerrados y sus modificaciones por fotoestimulación repetida en la misma condición sin ningún efecto importante para el sujeto. Como se mencionó antes la respuesta del sujeto (en potencia absoluta) va disminuyendo conforme se repite la estimulación, proceso que representa la habituación y que debe resultar de modificaciones de las conexiones funcionales de una red compleja de interconexiones, de manera global (toda la red) o de manera local (grupos regionales). Por lo tanto, en condiciones de Pre-estimulación deben tenerse unos niveles de conectividad funcional que se modifican por la activación de las vías sensoriales específicas e inespecíficas que cambia dichas conexiones funcionales. Esta activación provoca aumento de la actividad en otros circuitos relacionables con los procesos de atención a la fotoestimulación no esperada por el sujeto, de identificación de la estimulación y de asociación con cambios en su homeostasis. Procesos que le permiten calificar de significativa o no a la fotoestimulación. En este procedimiento esperamos que la califique como no significativa y disminuya su respuesta, es decir, se habitúe, disminuyendo las altas conectividades entre las diferentes áreas corticales.

Los procesos que se estudian pueden caracterizarse con la simple variación de la conectividad, sin embargo, las redes obtenidas al cuantificar la coherencia, son complejas y aún y cuando se tengan pocos nodos, tratar de interpretar la cantidad de conexiones resultantes (dependientes de la cantidad de nodos, en este caso n2-n2), resulta difícil de resumir y estudiar la información. Por esa razón este trabajo busca un siguiente paso, tratar de caracterizar el comportamiento de la población a través de los parámetros de redes durante el proceso de habituación.

Los resultados para G1 muestran, en el ritmo alfa, que los electrodos en los lóbulos frontales, parietales y occipitales en condición Pre presentan mayor número de conexiones que en la condición Rph en los primeros fotoestímulos (Figura 4). Los electrodos del lóbulo frontal presentan de 10 a 11 conexiones, los electrodos del lóbulo parietal de 9 a 10 conexiones y el lóbulo temporal alrededor de 8 conexiones. Los electrodos C3 y C4, presentan alrededor de 7 conexiones. El grado promedio, en los mismos fotoestímulos en la etapa Pre, oscila entre 9 y 11 conexiones promedio a valores superiores de 0.4 de coherencia (Figura 5). Ambos parámetros, grado de nodo y grado promedio, permiten observar el estado funcional inicial de los sujetos: mayor número de conexiones en los lóbulos frontales, parietales y en menor medida los occipitales y centrales. En los mismos fotoestímulos, la condición Rph, presenta una disminución de las conexiones en los electrodos de los lóbulos frontales (de 8 a 9), parietales (7 a 8), occipitales (de 7 a 8) y centrales (de 6 a 7), se puede observar que los sujetos al recibir los fotoestímulos iniciales presentan menores conexiones en dichos lóbulos, esto se observa en el comportamiento del grado promedio de las redes funcionales, se tienen en promedio de 6 a 7 conexiones, esto se puede interpretar como la respuesta inicial a la fotoestimulación: una disminución de conexiones a diferencia de la preparación del sujeto antes de recibir los estímulos.

Al continuar la fotoestimulación se disminuyen el número de conexiones en Pre y aumentan en Rph en los electrodos del arreglo, es decir, se observa un comportamiento inverso en las etapas Pre y Rph, de tal forma que al finalizar los fotoestímulos, las conexiones en etapa Rph superan ligeramente a la cantidad de conexiones de la etapa Pre. Dichos cambios ocurren en un intervalo de umbral de conectividad de 0.4 a 0.55 en el grado promedio, por debajo de un valor de 0.4, las conexiones promedio en todos los fotoestímulos en ambas etapas, es alrededor de 8 conexiones y por arriba de un valor de 0.55, se tiene un promedio menor a 2 conexiones.

En el caso de C-, en la Figura 6, se observa que en los primeros fotoestímulos (hasta el séptimo) se tiene un mayor valor en Rph a umbrales mayores que la etapa Pre. En los estímulos intermedios (8 al 13), el comportamiento es similar entre Pre y Rph, en los últimos fotoestímulos C- presenta un ligero aumento en ambas etapas. Podemos interpretar que en el estado inicial se genera un patrón de conectividad que considera mayor cantidad de grupos locales de interconexión. En el estado final ambas etapas presentan un aumento.

Figura 6 Representación de la variación de C- en función de los fotoestímulos, indicando el valor de C- para cada umbral mediante una paleta de colores. La etapa Rph y Pre se indica de la misma forma que la Figura 4. 

Para C- se observa un intervalo donde se presentan cambios en su comportamiento, entre 0.45 y 0.6. Por debajo de 0.4 se tiene un valor de C- máximo (uno), es decir, la red está altamente conectada de tal forma que la probabilidad de que los vecinos en torno a cualquier nodo estén altamente conectados es máxima. Por arriba de 0.65 la red presenta escaso número de grupos locales.

En la Figura 7, se observa el comportamiento de L-. En el que consideramos el estado inicial de conexión en condición Pre se tiene un valor aproximado de 2.5 de distancia mínima promedio, hasta umbrales cercanos a 0.8, esto continua hasta el estímulo 7, mostrando disminución gradual durante los fotoestímulos intermedios. En el caso de Rph, en los primeros fotoestímulos es menor el umbral (0.7) al que se obtienen valores de L-. En los siguientes fotoestímulos (8 a 13) se observa que en Rph se aumenta ligeramente el valor de L- a umbrales superiores. En los últimos fotoestímulos el comportamiento es inverso, Rph presenta a umbrales mayores una variabilidad a diferencia de Pre. El intervalo para el cual L- presenta alteraciones es entre 0.6 y 0.75 ya que por arriba de un umbral de 0.75 y por debajo de 0.6 no hay variaciones de distancias promedio.

Figura 7 Representación de la variación de L- en función de los fotoestímulos, indicando el valor de L- para cada umbral mediante una paleta de colores. La etapa Rph y Pre se indica de la misma forma que la Figura 4. 

De los resultados obtenidos de los parámetros de redes, nuestra percepción de la variación del comportamiento de los patrones de conectividad muestra que:

Un segundo análisis se realizó para tratar de observar el comportamiento de cada parámetro al promediar los datos obtenidos en la variación de umbral y a través de un modelo de regresión lineal obtener curvas de comparación entre ambas condiciones Pre y Rph. En la Figura 8, se muestran las curvas obtenidas de la regresión lineal, para el grado promedio (Figura 8 (arriba)), se muestra que en el estado inicial se tienen más de 4 conexiones en Pre y en Rph por debajo de 4 conexiones en promedio, conforme se presentan más fotoestímulos disminuye el grado de conexión de la red global, en Rph se tiene un comportamiento inverso, ya que aumenta gradualmente, hasta que en los últimos fotoestímulos se tiene similares niveles de conectividad. Este comportamiento también se presenta en C- (Figura 8 (abajo)), donde se inicia con un valor alrededor de 0.6 y disminuye conforme se presentan las fotoestimulaciones en Pre, caso contrario Rph, aumenta hasta llegar a niveles similares que Pre. En el caso de L- (Figura 9), no muestra una variación en función de los fotoestímulos, su comportamiento en Pre y Rph es estático prácticamente.

Figura 8 Curvas de regresión lineal del promedio de medidas de red vs número de fotoestimulación del grupo G1 en ritmo alfa: condición Pre (verde), condición Rph (azul). Grado promedio (arriba). C- (abajo). 

Figura 9 Curva de regresión lineal del promedio de medidas de red vs número de fotoestimulación del grupo G1 en ritmo alfa: condición Pre (verde), condición Rph (azul). L-. 

Con estas gráficas que representan la información promediada de los valores obtenidos de cada parámetro al umbralizar la conectividad, se puede observar un comportamiento común, se tienen más conexiones, y por tanto grupos locales de interconexión en un estado inicial en preparación de la fotoestimulación, que disminuye conforme se presentan los fotoestímulos.

Los cambios del ritmo alfa se relacionaron con procesos cognitivos desde su primera descripción por Berger en 1930 y esta relación ha sido ampliamente confirmada por múltiples investigaciones de las que Bafar hace una revisión y actualiza la hipótesis del funcionamiento del ritmo alfa [⁴⁵-⁴⁷]. Por otro lado, el ritmo beta se ha relacionado con procesos de sincronización de diversas redes para regular la conducta adaptativa correcta [⁴⁸]. En la actualidad se propone que cada rango de frecuencias se relaciona con funciones sensoriales, motoras y en particular cognitivas con base en conexiones funcionales que empezamos a estudiar con procesos de coherencia [⁴⁹].

Conclusiones

En este trabajo proponemos el uso de técnicas para el estudio de la conectividad funcional cerebral, obteniendo parámetros de redes que indiquen el comportamiento de los patrones de conectividad durante el fenómeno de la habituación, es decir, se pasa de un estado de desincronización a sincronización en el ritmo alfa dado que el sujeto se adapta o acostumbra a un estímulo [⁴⁵- ⁴⁷,⁴⁹]. La posibilidad de observar dichos cambios a nivel de conectividad, puede brindar un acercamiento del conocimiento de las estructuras globales de conectividad, cómo se relacionan las diferentes áreas del cerebro y describir su comportamiento a base de parámetros de redes complejas.

El estudio de la conectividad funcional cerebral es complejo, no se puede saber el valor real de la conectividad entre las señales de los electrodos a través de datos experimentales, ya que dependen de funciones de densidad de probabilidad de procesos estocásticos [²⁸], solo es posible la estimación a partir de diversos métodos, basados en diversas hipótesis sobre el modelo de conexión entre las señales analizadas. Adicionalmente, la interpretación de resultados puede ser limitada por la posible adición de un valor de bias influenciada por la elección de una referencia. La referencia situada en los lóbulos de las orejas, nariz o mastoide, puede añadir una variabilidad eléctrica y posteriormente una estimación de conectividad inflada. En la literatura se sugiere evitar el uso de una referencia común [²⁸,³¹-³³,⁵⁰], y una de las posibles soluciones es el uso de una referencia promediada o mediante un arreglo espacial, como el Laplaciano, sin embargo, no está claro cuál técnica es la mejor que puede representar las propiedades espectrales sin distorsión. En este trabajo los resultados obtenidos usando una referencia promediada, con un arreglo de 16 electrodos, pueden tener un valor de bias de conectividad, por ello la interpretación de los resultados es cuidadosa, al ser una primera aproximación en esta temática.

El uso de un arreglo de 16 electrodos, aunque es un número bajo de electrodos, se orienta a su aplicación en estudios clínicos, aumentar el número de electrodos, como se sugiere en la literatura (64 a 128 electrodos), hace difícil su aplicación en estudios clínicos, ciertamente se aleja de la aplicación de los desarrollo tecnológicos de los usos clínicos, sin embargo, se busca mejorar los procesos de diagnóstico y seguimiento de la evaluación de pacientes, por lo que en trabajos futuros se puede usar la información generada en casos clínicos de relevancia. Recalcamos que las interpretaciones a partir de los resultados de este trabajo, deben ser consideradas de manera exploratoria.

En reposo se analizan las conexiones funcionales, a través de la teoría de gráficas, establecidas entre el sistema retículo-talámico-paleó y neo-cortical que representan la condición de la red en cuatro frecuencias de oscilación. Los resultados de los parámetros grado de nodo y grado promedio, que cuantifican las conexiones que tiene cada electrodo y la red en promedio, respectivamente, muestran que se puede considerar un estado inicial de conectividad funcional de los sujetos sanos, la cantidad de conexiones en los lóbulos frontales, parietales, centrales y temporales son ligeramente superiores en la etapa Pre (estado inicial en preparación a la fotoestimulación) que en la etapa Rph. Lo anterior puede sugerir un incremento de la actividad oscilatoria simultánea en dichas regiones que podemos interpretar como indicativo de los procesos de identificación de la señal (fotoestimulación), de su percepción, significado y de memorización [⁵⁰,⁵²-55]. Además de presentarse una mayor proporción de grupos locales de interconexión, haciendo a la red inmune a la pérdida de conexiones entre nodos vecinos, la información podrá llegar a su destino.

Conforme se presentan las siguientes fotoestimulaciones, se observa un ligero aumento de conexiones en la etapa Rph y una disminución gradual en la etapa Pre, es decir, posiblemente se representa el fenómeno de habituación en la etapa Pre, que puede corresponder a las estimulaciones no significativas, en otras palabras, los sujetos al presentarse un estímulo ya conocido, se habitúan y el comportamiento de las conexiones cambia. Resulta interesante que al terminar la fotoestimulación los parámetros de redes tienen al mismo comportamiento en ambas etapas, pareciera que también puede considerarse un estado final de estimulación, ante las repeticiones la habituación se refleje con dichos cambios.

Variar el nivel de conectividad para conocer el comportamiento de los parámetros de la red binaria obtenida a partir de una umbralización, modifica la densidad y patrones de conexiones cerebrales. A partir de los resultados mostrados, se puede interpretar que para el proceso de habituación en sujetos sanos, estudiar niveles de conectividad por arriba de 0.6, puede no reflejar información que muestre la dinámica cerebral, debido a que, en todos los parámetros calculados los valores eran estáticos; esto mismo sucede en valores menores a 0.4 de conectividad, es decir, es posible que niveles altos y bajos de conectividad no sean una respuesta normal de un sujeto ante la fotoestimulación bajo el protocolo usado. Además, podemos considerar que sea adecuado estudiar un intervalo de conectividad entre 0.4 y 0.6, ya que es el intervalo en el que el grado promedio y C- muestran la dinámica en relación a la cantidad de conexiones de cada red. Para el parámetro L, a partir de nuestros resultados consideramos que el intervalo de estudio debe ser de 0.6 a 0.75.

Adicionalmente, los resultados mostrados en la Figura 8, respaldan el comportamiento mostrado en las figuras anteriores: es posible indicar un estado de conectividad inicial y final, el comportamiento (de los parámetros) en la etapa Pre tiende a disminuir conforme se presentan las fotoestimulaciones y en Rph tienden a aumentar hasta llegar a valores similares. En el segundo análisis el valor de L- (Figura 9), no cambia entre fotoestimulaciones, la representación del parámetro puede ser afectada por el promedio realizado, ya que pueden existir valores extremos que impiden observar de la variabilidad, a diferencia de que en el grado promedio y C- si muestran información al respecto.

Es importante mencionar que estas conclusiones se aplican al ritmo alfa ya que en los ritmos delta, theta y beta prácticamente no se observaron cambios importantes en el funcionamiento de los nodos de la red. En el grupo G3 aparece un comportamiento similar al del grupo G1 pero de menor intensidad en el ritmo alfa, indicando que este grupo presenta el mismo fenómeno al iniciar con mayores conexiones en la condición Pre que Rph. En los ritmos delta, theta y beta no muestran cambios de conectividad.

En este trabajo proponemos una forma de abordar el problema de la umbralización de manera exploratoria; realizar un barrido de conectividad para conocer si los parámetros de redes presentan cambios en su comportamiento en la etapa de Pre habituación y durante la habituación, en función de un intervalo de coherencia. De este modo, en un trabajo futuro se puede focalizar el análisis de los intervalos propuestos.

Aplicar la teoría de gráficas al estudio de la conectividad funcional cerebral permite evaluar las hipótesis de la funcionalidad del cerebro con base en redes o circuitos propuestas desde la década de los sesentas. Con estos procedimientos estamos estableciendo patrones característicos, en sujetos sanos que tomaremos como posibles biomarcadores en el apoyo de los procesos diagnósticos y de rehabilitación integral de pacientes con enfermedades cerebrales adquiridas.