Análisis de actividad de los nervios y conectividad mediante datos EEG intracraneal con Software SPM

Resumen

Presentamos dos protocolos analíticos que pueden utilizarse para analizar los datos de electroencefalografía intracraneal mediante el software de mapeo paramétrico estadístico (SPM): Análisis de mapeo paramétrico estadístico tiempo-frecuencia de actividad de los nervios y dinámica causal modelado de las respuestas inducidas para la conectividad intra e interregional.

Resumen

Medición de actividad de los nervios y la conectividad asociada a las funciones cognitivas en alta resolución espacial y temporal es un objetivo importante en Neurociencia cognitiva. Intracraneal electroencefalografía (EEG) directamente puede registrar la actividad eléctrica de los nervios y tiene la capacidad única para alcanzar esta meta. Tradicionalmente, con un promedio de análisis se ha aplicado para analizar los datos de EEG intracraneales; sin embargo, varias nuevas técnicas están disponibles para que la conectividad neuronal de la actividad y intra e interregional. Aquí, presentamos dos protocolos analíticos recientemente se aplicó para analizar los datos de EEG intracraneales mediante el software de mapeo paramétrico estadístico (SPM): Análisis de tiempo-frecuencia SPM para la actividad de los nervios y modelado causal dinámica de inducidas por conectividad intra e interregional. Divulgamos el análisis de datos de EEG intracraneales durante la observación de las caras como resultados representativos. Los resultados revelaron que el giro occipital inferior (IOG) demostró actividad gamma-banda en etapas muy tempranas (110 ms) en respuesta a caras y el GSI y amygdala demostrada rápida intra e interregional conectividad mediante varios tipos de oscilaciones. Estos protocolos analíticos tienen el potencial para identificar los mecanismos neuronales subyacentes a las funciones cognitivas con altos perfiles espaciales y temporales.

Introducción

Medición de actividad de los nervios y la conectividad asociada a las funciones cognitivas en altas resoluciones espaciales y temporales es uno de los objetivos principales de la Neurociencia cognitiva. Sin embargo, no es fácil lograr este objetivo. Un método popular usado para registrar la actividad neuronal es la proyección de imagen de resonancia magnética funcional (MRI). Aunque MRI funcional ofrece varias ventajas, tales como una alta resolución espacial a nivel milimétrico y grabación no invasivo, una desventaja clara de la RMN funcional es su baja resolución temporal. Además, MRI funcional mide señales dependientes de nivel del oxígeno de sangre, que sólo indirectamente reflejan la actividad eléctrica de los nervios. Métodos electrofisiológicos, como electroencefalografía (EEG) y magnetoencefalografía (MEG), tienen alta resolución temporal en el nivel de milisegundo. Sin embargo, tienen resolución espacial relativamente baja, ya que registran las señales eléctricas/magnéticas en el cuero cabelludo y debe resolver problemas inversos para representar la actividad cerebral.

EEG intracraneal puede grabar directamente la actividad eléctrica neuronal en alta temporal (milisegundo) y espacial (cm) resolución¹. Esta medida puede proporcionar valiosas oportunidades para entender la actividad de los nervios y la conectividad, aunque tiene limitaciones claras (por ejemplo, regiones medibles están restringidas a criterios clínicos). Varios estudios de EEG intracraneales han aplicado análisis promedio tradicional para representar la actividad de los nervios. Aunque con un promedio de análisis sensible puede detectar activación de banda de tiempo-bloqueada y baja frecuencia, no se detecta no-fase-bloqueado o activación de alta frecuencia (por ejemplo, la banda de gamma). Además, acoplamiento neuronal funcional no ha sido analizada en profundidad en la literatura en las grabaciones de EEG intracraneales. Varias nuevas técnicas se han desarrollado recientemente para representar neuronal actividad intra e interregional conectividad funcional grabaciones RM y EEG/MEG, que se puede aplicar para analizar los datos de EEG intracraneales.

Aquí, presentamos protocolos analíticos que recientemente hemos aplicado para analizar los datos de EEG intracraneales mediante el software de mapeo paramétrico estadístico (SPM). Primero, a revelar Cuándo y con que frecuencia, podrían activar las regiones cerebrales, realizamos tiempo-frecuencia SPM análisis². Este análisis descompone los dominios de tiempo y frecuencia simultáneamente usando una transformación wavelet continua y corrige adecuadamente la tasa de error family-wise (EFT) en los mapas de frecuencia de tiempo usando la teoría de campo al azar. En segundo lugar, para descubrir cómo se comunican las regiones del cerebro, aplicamos modelado causal dinámica (DCM) de respuestas inducido⁴. DCM permite la investigación de la conectividad efectiva (es decir, las influencias causales y direccionales entre las regiones de cerebro⁵). Aunque DCM fue propuesto originalmente como una herramienta para el análisis funcional MRI datos⁵, DCM de respuestas inducidas ha sido ampliado para analizar los espectros de energía varían con el tiempo de las señales electrofisiológicas⁴. Este análisis permite la representación de ambos intra e interregional neural conectividad. Varios estudios neurofisiológicos han sugerido que los cálculos intra regionales locales y comunicación interregional a larga distancia utilizan principalmente las oscilaciones de la theta y gamma banda, respectivamente, y sus interacciones (por ejemplo, entrainments) puede reflejarse por theta gamma Cruz-frecuencia acoplamiento^3,6,7,8. Este informe se centra en el protocolo de análisis de datos; para un resumen de información de fondo^9,10 y de protocolos de grabación¹¹ de EEG intracraneal, refiera por favor a la literatura.

Protocolo

Nuestro estudio fue aprobado por el Comité de ética institucional local.

1. básica información

Nota: Los protocolos analíticos pueden aplicarse a varios tipos de datos sin restricciones en cuanto a participantes específicos, electrodos, métodos de referencia o electrodo ubicaciones. En nuestro ejemplo, hemos probado a seis pacientes que sufren de epilepsia focal intratable farmacológicamente. Probamos a pacientes que no tenían los focos epilépticos en las regiones de interés.

1. Registrar datos de EEG intracraneales durante el experimento cognitivo en los electrodos de destino.
   1. Implante de electrodos de profundidad utilizando el método estereotáctico¹².
   2. Utilizar electrodos de platino subdurales (diámetro: 2,3 mm) y electrodos de platino de profundidad (diámetro: 0,8 mm) para medir simultáneamente la actividad cortical y subcortical, respectivamente.
   3. Coloque los electrodos de referencia en la superficie del cráneo de la región frontal dorsal de la línea media, con los contactos de los electrodos hacia el cráneo para evitar activación referencial¹².
   4. Ampliar datos, filtro en línea (pase de banda: 0.5 – 300 Hz) y muestra a 1000 Hz.
   5. Para registrar y quitar estadísticamente artefactos asociados con movimientos de los ojos, además de grabar electrooculograms. Seleccionar los electrodos de destino basados en intereses teóricos. Además, utilizar datos de tomografía computada y MRI individual para comprobar la ubicación del electrodo.
2. De la muestra y preprocesar datos EEG intracraneales (figura 1).
   Nota: Los protocolos analíticos pueden aplicarse a varios tipos de datos sin restricción de longitud de datos específico o métodos de preprocesamiento.
   1. Aquí, datos de la muestra durante ms de 3.000 (antes estímulo: 1.000 ms; posterior al estímulo: 2.000 ms) para cada ensayo.
   2. Porque los participantes aquí demostrados anormalmente actividad de alta amplitud en algunos ensayos, posiblemente relacionados con la epilepsia, excluyen estos ensayos aislados con umbrales predefinidos (> 800 μV y > 5 SD). Otros pasos de preprocesamiento, incluyendo inspección visual y análisis de componentes independientes, pueden ser necesarios dependiendo de los objetivos experimentales y condiciones.
3. Convertir el formato nativo del sistema de EEG en un formato basado en MATLAB SPM (cf. SPM12 Manual de 12.1 y 12.2).
   Nota: La mayoría de formatos de datos EEG puede ser importada directamente en software SPM seleccionar conversión en la interfaz del Editor de lote de SPM y especificando todos los parámetros de entrada requeridos. Otra forma posible es usar un script de ejemplo "spm_eeg_convert_arbitrary_data.m" en el subdirectorio de hombre/example_scripts del directorio de programa SPM. Esta secuencia de comandos proporciona una forma conveniente para convertir un archivo ASCII o estera que se puede exportar por muchos sistemas de EEG con formato SPM.

2. tiempo-frecuencia SPM análisis

1. Configurar SPM12 (http://www.fil.ion.ucl.ac.uk/spm/) y utiliza el M/EEG menú analítica¹³ (figura 2).
2. Realizar el tiempo – frecuencia SPM análisis² seleccionando "análisis tiempo-frecuencia" en el menú SPM de preprocesados datos EEG intracraneales de cada ensayo mediante descomposición wavelet continua con wavelets de Morlet basadas en parámetros predefinidos ( Figura 3).
   Nota: Transformaciones Wavelet revelan la evolución temporal de las componentes espectrales por convolución tamaño datos de EEG intracraneales con wavelets de múltiples frecuencias¹⁴.
   1. Aquí, llevar a cabo utilizando ciclo de siete Morlet wavelets para la época toda la descomposición wavelet (-1,000 – 2.000 ms) y rango de frecuencia de 4 – 300 Hz.
   2. Determinar la wavelet madre y número de ciclos basada en un anterior estudio¹⁵. Tenga en cuenta que el número de ciclos en la wavelet controla el tiempo – frecuencia resoluciones y se recomienda ser mayor de 5 para asegurar la estabilidad de valoración¹³.
   3. Determinar rangos de tiempo y frecuencia basadas en el interés de investigación.
3. Los mapas de tiempo – frecuencia resultantes de cultivos automáticamente seleccionando "Cultivo" en el menú SPM para eliminar los efectos de borde. Aquí, los mapas de tiempo – frecuencia de los cultivos en -200 – 500 ms.
4. Realizar la transformación de datos (opcional) y corrección de línea base mediante la selección de "reescalar tiempo-frecuencia" en el menú SPM para los mapas de tiempo – frecuencia visualizar mejor los cambios de energía relacionados con eventos y mejorar la normalidad de los datos.
   Nota: Aquí, los mapas de tiempo – frecuencia fueron transformados registro y línea de base (-200-0 ms)-corregido.
5. Convertir los mapas de tiempo – frecuencia imágenes bidimensionales (2D) seleccionando "Convert2Images" en el menú SPM.
   1. Liso con un núcleo Gaussiano con un valor de (FWHM) predefinido lleno-ancho máximo de la mitad para compensar la variabilidad entre sujeta y se ajustan a los supuestos de la teoría de campo al azar en la inferencia estadística^2,13 .
      Nota: Aquí, los mapas de tiempo – frecuencia se suavizan con un kernel gaussiano de FWHM de 96 ms en el dominio del tiempo y 12 Hz en el dominio de la frecuencia basado en un estudio anterior².
6. Entrar las imágenes 2D en el modelo lineal general seleccionando especificar 1 "nivel" en el menú SPM.
7. Estimar el modelo lineal general seleccionando "Estimación de modelo" en el menú SPM.
8. Realizar inferencias estadísticas sobre los datos de tiempo – frecuencia SPM {T} basados en la teoría del campo al azar² seleccionando "Resultados" en el menú SPM. Detectar activados mucho tiempo – frecuencia racimos con umbrales predefinidos (posiblemente corregidos para comparaciones múltiples).
   Nota: Aquí, el umbral de medida de p < 0.05, que es EFT-corregido para comparaciones múltiples, con un límite de altura de p < 0.001 (sin corregir) fue utilizado.

3. DCM de respuestas inducidas

1. Configurar SPM12 (http://www.fil.ion.ucl.ac.uk/spm/) y utiliza el M/EEG menú analítica¹³ (figura 4).
   1. Iniciar análisis de DCM haciendo clic en botón "DCM" en el menú SPM. Activar DCM inducido respuestas seleccionando "IND" en el cuadro de lista. Importar los datos de EEG intracraneales preprocesados haciendo clic en "nuevos datos" en el DCM para menú M/EEG.
2. Especificar la ventana de tiempo de interés, las condiciones de interés, contrastes para las condiciones seleccionadas (Esto define las entradas de modulación utilizadas más adelante en la especificación de red), ventana de frecuencia de interés y el número de los ciclos de la wavelet en el DCM de menú () M/EEG Figura 5).
   1. Utilizar cinco-ciclo Morlet wavelets (4 – 100 Hz en pasos de 1 Hz) y la ventana de tiempo a ms de 1, 500.
   2. Determinar el ciclo de wavelet con arreglo a la configuración por defecto. Tenga en cuenta que la recomendación de software es el valor mayor de 5 para asegurar la estabilidad de valoración¹³. Se determinaron los rangos de tiempo – frecuencia basada en nuestro interés de investigación. Tenga en cuenta que una ventana de tiempo con un ± adicional 512 ms automáticamente se utilizó durante el cómputo para eliminar los efectos de borde.
3. Basado en el marco de DCM^4,5, definen el (1) conducción de entradas, que representan las entradas sensoriales de Estados neuronales; (2) las conexiones intrínsecas, que incorporan la conectividad de línea de base entre los Estados neuronales y conexiones por cuenta propia; y (3) efectos moduladores en las conexiones intrínsecas vía experimental manipulaciones para modelos nulos e hipotéticas. También definir el tipo de modulación lineal (en frecuencia) o no lineales (entre frecuencias).
   1. Especificar conexiones (lineales y no lineales) intrínsecas, conducción y entradas de modulación en el DCM para menú M/EEG.
   2. Modificar la configuración por defecto de algunos parámetros relacionados (p. ej., tiempo de inicio del estímulo previo y duración) si es necesario. Estimar los modelos seleccionando "DCM de invertir" en el DCM para menú M/EEG. Después de seleccione Guardar resultados como img guardar modulador frecuencia acoplamiento imágenes de parámetro.
4. Realizar un análisis efectos aleatorios Bayesiano modelo selección (BMS)¹⁷ haciendo clic en "BMS" en el DCM para menú M/EEG para identificar el modelo de red óptimo. Utilizar las probabilidades del modelo esperado o probabilidades de excedencia como medidas de evaluación.
5. Hacer inferencias con respecto a los patrones de la Cruz-frecuencia de las conexiones del moduladores usando los parámetros del modelo ganador mediante el menú SPM (véase paso 2).
   1. Suavizar las imágenes del parámetro de acoplamiento modulador seleccionando "Convert2Images" en el menú SPM.
   2. Realizar análisis de modelo lineal general seleccionando "Especificar nivel de 2 º" en el menú SPM.
   3. Calcular los valores SPM {T} 2D seleccionando "Resultados" en el menú SPM.
      Nota: Aquí, el FWHM se fijó en 8HZ, basado en un estudio anterior⁴. Valores significativos exploratorily se identificaron usando un umbral de altura de p < 0.05 (sin corregir).

Resultados

Utilizando el protocolo aquí presentado, se analizaron datos de EEG intracraneales en respuesta a caras^18,19. Registraron datos de seis pacientes durante la observación pasiva de caras, mosaicos y casas en posición vertical y había invertida orientaciones. Los contrastes de comparación vertical mosaicos de caras verticales y caras vertical versus casas vertical revelaron el efecto de la cara (es decir,

Discusión

Los protocolos analíticos para los datos de EEG intracraneales mediante el software SPM introducido aquí tienen varias ventajas en comparación con MRI funcional. En primer lugar, los protocolos pueden representar activación neural en una alta resolución temporal. Por lo tanto, los resultados indican que si se implementan los correlatos cognitivos de la activación neuronal en etapas iniciales o finales de procesamiento. En nuestro ejemplo, el efecto de la cara fue identificado durante las primeras etapas (es dec...

Materiales

Name	Company	Catalog Number	Comments
	none