SciELO - Scientific Electronic Library Online

 
vol.16 número1Solución basada en herramientas de software libre para la implementación del teletrabajo online en empresas cubanasEstudio experimental con varios enfoques de asignación de demanda en el problema de localización de máxima cobertura capacitado índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

  • No hay articulos citadosCitado por SciELO

Links relacionados

  • No hay articulos similaresSimilares en SciELO

Compartir


Revista Cubana de Ciencias Informáticas

versión On-line ISSN 2227-1899

Rev cuba cienc informat vol.16 no.1 La Habana ene.-mar. 2022  Epub 01-Mar-2022

 

Artículo original

Cuantificación de la perfusión cerebral por SPECT utilizando co-registro con IRM y atlas anatómicos

Quantification of cerebral perfusion SPECT based on MRI co-registration and anatomic atlas

Karla Batista García-Ramó1  * 
http://orcid.org/0000-0003-1042-7369

Nancy Pavón Fuentes2 
http://orcid.org/0000-0003-0998-8601

Lilia Morales Chacón3 
http://orcid.org/0000-0003-0205-0733

Angel Aguila Ruiz4 
http://orcid.org/0000-0002-8166-5877

Nelson Rodríguez Mesa5 
http://orcid.org/0000-0002-6827-2964

1 Departamento de Medicina Nuclear, Centro Internacional de Restauración Neurológica. Ave 25 No. 15805 entre 158 y 160, La Habana, Cuba. karlabg@infomed.sld.cu

2 Departamento de Neuroinmunología, Centro Internacional de Restauración Neurológica. Ave 25 No. 15805 entre 158 y 160, La Habana, Cuba. nancy@neuro.ciren.cu

3 Departamento de Neurofisiología, Centro Internacional de Restauración Neurológica. Ave 25 No. 15805 entre 158 y 160, La Habana, Cuba. lily@neuro.ciren.cu

4 Departamento de Medicina Nuclear, Centro Internacional de Restauración Neurológica. Ave 25 No. 15805 entre 158 y 160, La Habana, Cuba. angelfa@infomed.sld.cu

5 Departamento de Medicina Nuclear, Centro de Investigaciones Médico-Quirúrgicas. Calle 216 y 11b, La Habana, Cuba. nelsonrdguez@infomed.sld.cu

RESUMEN

La cuantificación relativa del flujo sanguíneo cerebral con tomografía por emisión de fotón único presenta distintas limitaciones relacionadas con la pobre resolución espacial, la selección de las regiones de interés y la selección de la región de referencia. El objetivo del presente trabajo es proponer una metodología para la cuantificación relativa del flujo sanguíneo cerebral a partir del co-registro de la imagen funcional con la imagen estructural de resonancia magnética y utilizando atlas anatómicos para la definición de las regiones. La metodología incluye tres pasos principales integrados en un pipeline en python: pre-procesamiento espacial, corrección del volumen parcial, y cálculo de los índices de perfusión y asimetría. La metodología se implementó en 17 pacientes con epilepsia y 10 pacientes con enfermedad de Parkinson. En los pacientes con epilepsia, los resultados mostraron un incremento tanto en índice de perfusión como en el índice de asimetría en el estudio ictal en comparación con el interictal en las regiones correspondientes a la zona epileptógena. En los pacientes con enfermedad de Parkinson se identificaron regiones con diferencias significativas en ambos índices entre pacientes en distintos estadios. La metodología para la cuantificación del flujo sanguíneo cerebral a partir de la tomografía por emisión de fotón único, resonancia magnética y atlas anatómicos se ofrece como herramienta en la evaluación pre-quirúrgica de pacientes con epilepsia y para la evaluación del cuadro neurodegenerativo en la enfermedad de Parkinson.

Palabras clave: tomografía por emisión de fotón único; perfusión cerebral; cuantificación; epilepsia; enfermedad de Parkinson

ABSTRACT

The relative quantification of the regional cerebral blood flow by single-photon emission computed tomography presents different limitations related to the poor spatial resolution, the selection of the regions of interest and the selection of the reference region. The objective of the present work is to propose a methodology for the relative quantification of cerebral blood flow from the co-registration of functional and structural images and using anatomical atlases for the definition of regions. The methodology includes three main steps integrated into a pipeline in python: spatial pre-processing, partial volume correction, and calculation of the perfusion and asymmetry indexes. The methodology was implemented in 17 patients with epilepsy and 10 patients with Parkinson's disease. In patients with epilepsy, the results showed an increase in both perfusion and asymmetry index in the ictal study compared to the interictal one in the regions corresponding to the epileptogenic zone. In Parkinson's disease patients, regions with significant differences in both perfusion and asymmetry index were identified between patients at different stages. The methodology for the quantification of cerebral blood flow from single-photon emission computed tomography, magnetic resonance image and atlases is offered as a tool for pre-surgical evaluation of patients with epilepsy and for the evaluation of the neurodegenerative condition in Parkinson's disease.

Keywords: single-photon emission computed tomography; cerebral blood flow; quantification; epilepsy; Parkinson's disease

Introducción

Algunos métodos han sido propuestos para la evaluación cuantitativa de los estudios por tomografía por emisión de fotón único (SPECT) (Matsuda 2021; Sánchez-Catasús et al. 2002). También en los programas de procesamiento de algunos equipos de SPECT se proponen herramientas básicas para la cuantificación. No obstante, los métodos propuestos presentan algunas limitaciones relacionadas con la pobre resolución espacial propia de la técnica, o con la selección de las regiones de interés, o la región de referencia utilizada para la cuantificación relativa, entre otros. El objetivo del presente trabajo es proponer una metodología para la cuantificación relativa de la perfusión cerebral a partir del co-registro de la SPECT con la imagen de resonancia magnética (IRM) y utilizando atlas anatómicos para la definición de las regiones de interés (ROIs). La metodología es implementada en dos patologías distintas: epilepsia farmacorresistente y enfermedad de Parkinson (EP).

Adquisición de las imágenes de resonancia magnética (IRM) y de tomografía por emisión de fotón único (SPECT)

Cada paciente fue estudiado en un equipo de resonancia magnética Magnetom Aera SIEMENS de 1.5 Tesla. El protocolo utilizado incluye imágenes sagitales T1-MPRAGE 3D. Para la adquisición de la SPECT, a cada paciente se le administró una dosis de 630 MBq y se utilizó una cámara gamma de doble cabezal MiE (ECAM, Siemens). Se empleó un colimador de tipo baja energía y alta resolución (LEHR), matrices de 128x128 píxeles, 128 proyecciones en órbita circular de 360° y con una ventana energética del 20%.

Metodología para la cuantificación de la perfusión cerebral

La metodología incluye tres pasos principales: pre-procesamiento espacial, corrección del volumen parcial, y cálculo de los índices de perfusión.

Pre-procesamiento de las imágenes

El pre-procesamiento espacial incluye el co-registro IRM-SPECT y definición de las ROIs:

Delimitación del encéfalo a partir de la IRM anatómica T1 eliminando las estructuras que rodean al cerebro utilizando la herramienta BET del programa FSL (Jenkinson et al. 2012).

Segmentación de la IRM en mapas de sustancia gris, sustancia blanca y líquido cefalorraquídeo utilizando la herramienta FAST (Zhang, Brady, and Smith 2001) del FSL.

Multiplicación del atlas co-registrado con la máscara obtenida de la IRM de cada paciente.

Co-registro de la IRM anatómica, las máscaras de sustancia gris y sustancia blanca y el atlas individual de cada paciente con la SPECT utilizando FLIRT (Greve and Fischl 2009) del FSL.

Para el análisis de pacientes con enfermedad de Parkinson, el atlas se construyó co-registrando el atlas AAL (Automated Anatomical Labeling) (Rolls et al. 2020) con un atlas de estructuras subcorticales (Pauli, Nili, and Tyszka 2018) para incluir el núcleo subtalámico (STH) y la segregación del pallidum (externo e interno) en el análisis. Para el análisis de los pacientes con epilepsia, se utilizó el atlas propuesto por Oishi y col. (Oishi et al. 2009) porque incluye el tronco encefálico, una región que ha sido relacionada con la muerte súbita en la epilepsia.

Corrección parcial del volumen (CVP)

Para la CVP de las imágenes de SPECT se utilizó un método que utiliza la información de las segmentaciones de la IRM y calcula los valores de actividad media del ROI utilizando el método de matriz de transferencia geométrica, y luego realiza una corrección voxel a voxel (Thomas et al. 2011).

Cálculo de los índices de perfusión

El índice de perfusión (IP) se calcula dividiendo la actividad media en una región (AROI) por la actividad media del resto de la materia gris (Aglobal) (ecuación 1). También se utilizó el índice de asimetría (IA), donde los valores positivos indican dominancia del hemisferio izquierdo y los valores negativos del derecho (ecuación 2).

(1)

(2)

I representa la actividad media de los voxeles de la región de interés en el hemisferio izquierdo y D representa la actividad media de la región de interés en el hemisferio derecho.Figura 1

Fig. 1 Esquema de la metodología de cuantificación del FSCr. CVP se refiere a la corrección del volumen parcial.  

Selección de sujetos

La metodología propuesta fue aplicada retrospectivamente a 17 pacientes con epilepsia farmacorresistente (Tabla 1). Estos pacientes pertenecen al Programa de Cirugía en la Epilepsia de Cuba dirigido por el Centro Internacional de Restauración Neurológica (CIREN, Cuba) (Morales Chacon et al. 2018). La metodología también se aplicó de forma retrospectiva en 10 pacientes con EP (Tabla 2) provenientes del servicio de Mediciana Nuclear del CIREN.

Tabla 1 Características clínicas y demográficas de los pacientes con epilepsia farmacorresistente 

Tabla 2 Características demográficas y clínicas de pacientes con EP 

Resultados y discusión

El presente trabajo propone una metodología para la cuantificación relativa de la perfusión cerebral por SPECT. La metodología se basa en cuantificación por regiones integrando información estructural y funcional y utilizando atlas anatómicos. El co-registro de la IRM con la SPECT y la CVP son muy útiles en aras de tratar con el problema de la baja resolución de la SPECT. Los límites de las ROIs son definidos a partir de un atlas anatómico que se seleccione a priori para el estudio y, como parte de la metodología, se propone un pre-procesamiento para ajustar dicho atlas a la anatomía individual del sujeto. Con la utilización del atlas se eliminan las diferencias que se daban al dibujarlas por distintos usuarios en la rutina clínica. Las regiones de referencia más utilizadas son el cerebelo, tronco encefálico, puente, lóbulo temporal, materia gris y materia blanca. La no estandarización y la heterogeneidad para la normalización de la intensidad tiene un efecto en la cuantificación de la SPECT y en la comparación de los resultados de distintos estudios. En la metodología propuesta se utiliza como región de referencia la materia gris pero la metodología no restringe la selección de la misma. Un aspecto determinante en la cuantificación es la calidad de las imágenes, tanto de la SPECT y la IRM, y de cada uno de los métodos empleados en el procesamiento de ambas. Un elemento crucial lo constituye el efecto de volumen parcial y los métodos de corrección de dicho efecto (CVP), al punto de que métodos distintos arrojan diferentes resultados (Greve et al. 2016). La calidad de la IRM también es determinante, dado que para la metodología que se propone, influye en la definición de las regiones de interés y en la corrección del efecto de volumen parcial. No obstante, evaluaciones adicionales serán requeridas para determinar la influencia de los distintos métodos de CVP y de segmentación en la metodología que se propone. La metodología se implementó en 17 pacientes con epilepsia farmacorresistente, y 10 pacientes con EP. Para cada uno de estos trastornos se utilizaron distintos atlas anatómicos. Como parte de la evaluación pre-quirúrgica, se aplicó una metodología (Batista García-Ramó et al. 2019) que combina la estimación de las fuentes generadoras del EEG ictal (ESI) y la sustracción del SPECT interictal del ictal co-registrado con la IRM (SISCOM). Los resultados sugieren que la zona epileptógena (ZE) está localizada en lóbulo occipital derecho. La metodología para cuantificar el FSCr fue aplicado al SPECT ictal e interictal (Figura 2). Los IP muestran hiperperfusión en el estudio ictal con respecto al interictal en las ROI del lóbulo occipital derecho (Figura 2B), que corresponden a la ZE. Esto también se refleja en el IA mostrando una mayor asimetría en el estudio ictal. En la Figura 2C, observamos una superposición del SISCOM (rojo), ESI (azul) e IP (verde), esta última con mejor definición anatómica de la ZE. Este primer resultado refleja la capacidad de la metodología como una herramienta adicional para la localización de la ZE.

Fig. 2 Paciente 4. A) y B) son los estudios inter ictal e ictal, respectivamente. De izquierda a derecha: Co-registro de la IRM con estudio SPECT, co-registro de IRM con mapa de los IP para cada ROI y tabla con métricas de cuantificación del FSC en la ZE. C) Rojo: ZE estimada a través de la metodología SISCOM; Azul: ZE estimada por la solución inversa del EEG ictal (ESI). El mapa de IP para la ZE se representa en verde.  

A partir del método de k-medias se obtuvieron las regiones caracterizadas por una hiperperfusión (Figura 3) e hipoperfusión (Figura 4) en la epilepsia del lóbulo temporal. En ambas figuras se aprecian diferencias entre los distintos estados, ictal e interictal. Este resultado contribuye al enfoque de la epilepsia como trastorno de redes neurales permitiendo identificar cuáles serían las regiones que conforman la red epileptógena. Para ello se requeriría de un grupo de control.

Fig. 3 Patrones de hiperperfusión en el estado ictal e interictal de pacientes con epilepsia del lóbulo temporal. 

Fig. 4 Patrones de hipoperfusión en el estado ictal e interictal de pacientes con epilepsia del lóbulo temporal. 

De los 10 pacientes con EP, 5 pacientes se encuentran en los estadios iniciales y 5 en un estadio avanzado con evidencia de deterioro cognitivo. A partir del método de agrupamiento se identificaron las regiones con hipoperfusión (Figura 5) e hiperperfusión (Figura 6). Un resultado interesante es la hipoperfusión en el lóbulo frontal en ambos hemisferios para un estadio avanzado de la enfermedad. Esto pudiera estar en correspondencia con la presencia de deterioro cognitivo en estos pacientes.

Fig. 5 Patrones de hipoperfusión de ambos hemisferios en el estado inicial y avanzado de la EP. 

Fig. 6 Patrones de hiperperfusión de ambos hemisferios en el estado inicial y avanzado de la EP. 

En la Tabla 3 se muestran las regiones con coeficientes de correlación de Spearman significativos entre los IP y las distintas variables clínicas.

Tabla 3 Correlación entre IP y variables clínicas 

HC: hemisferio contralateral. HI: hemisferio ipsilateral

Llaman la atención regiones que no tienen correlación significativa con la edad pero sí con los años de evolución, la UPDRS y el H&Y. Según la Tabla 3, existe una correlación negativa entre la evolución de la EP y la perfusión en el giro de Heschl y el cingulado, no asociadas al envejecimiento. Estudios previos han encontrado diferencias significativas en esta región entre pacientes con EP y sujetos controles en estudios por neuroimágenes tanto funcionales como estructurales (Deng et al. 2021). Sin embargo, la correlación es positiva para el núcleo pulvinar del tálamo y el cerebelo. Owens-Walton y col. (Owens-Walton et al. 2019) reportan un aumento de la conectividad funcional en los núcleos del tálamo en los pacientes con EP con respecto al grupo control. Estos resultados sugieren que la cuantificación del FSC por SPECT pudiera contribuir a la búsqueda de biomarcadores basados en imágenes para el estudio de la EP.

Conclusiones

El presente trabajo propone una metodología para la cuantificación relativa de la perfusión cerebral por SPECT. La metodología se basa en cuantificación por regiones integrando información estructural y funcional y utilizando atlas anatómicos. La novedad del presente trabajo consiste en el enfoque multimodal de la metodología y en su capacidad para describir patrones de perfusión. En la bibliografía consultada no se encuentran estudios previos que describan los patrones de perfusión a partir de la SPECT ni en la epilepsia ni en la EP. La principal limitación se debe al pequeño número de pacientes estudiados. Este impide llegar a conclusiones sobre las patologías estudiadas pero permiten mostrar las bondades de la metodología y sirven como prueba de concepto para fundamentar estudios futuros con mayor número de pacientes.

Referencias

Batista García-Ramó, Karla, Carlos A. Sanchez Catasus, Lilia Morales Chacón, Angel Aguila Ruiz, Abel Sánchez Corneaux, Pedro Rojas López, and Jorge Bosh Bayard. 2019. “A Novel Noninvasive Approach Based on SPECT and EEG for the Location of the Epileptogenic Zone in Pharmacoresistant Non-Lesional Epilepsy.” Medicina 55 (8): 478. [ Links ]

Deng, Xia, Zheng Liu, Qin Kang, Lin Lu, Yu Zhu, and Renshi Xu. 2021. “Cortical Structural Connectivity Alterations and Potential Pathogenesis in Mid-Stage Sporadic Parkinson’s Disease.” Frontiers in Aging Neuroscience 13 (May): 650371-650371. https://doi.org/10.3389/fnagi.2021.650371. [ Links ]

Greve, Douglas N, and Bruce Fischl. 2009. “Accurate and Robust Brain Image Alignment Using Boundary-Based Registration.” NeuroImage 48 (1): 63-72. https://doi.org/10.1016/j.neuroimage.2009.06.060. [ Links ]

Greve, Douglas N., David H. Salat, Spencer L. Bowen, David Izquierdo-Garcia, Aaron P. Schultz, Ciprian Catana, J. Alex Becker, et al. 2016. “Different Partial Volume Correction Methods Lead to Different Conclusions: An 18F-FDG-PET Study of Aging.” NeuroImage 132 (May): 334-43. https://doi.org/10.1016/j.neuroimage.2016.02.042. [ Links ]

Jenkinson, Mark, Christian F. Beckmann, Timothy E. J. Behrens, Mark W. Woolrich, and Stephen M. Smith. 2012. “FSL.” NeuroImage 62 (2): 782-90. https://doi.org/10.1016/j.neuroimage.2011.09.015. [ Links ]

Matsuda, Hiroshi. 2021. “Quantification in Brain SPECT: Noninvasive Cerebral Blood Flow Measurements Using 99m Tc-Labeled Tracers.” In PET and SPECT in Neurology, Second, 53-72. Springer Nature Switzerland AG. https://doi.org/10.1007/978-3-030-53168-3. [ Links ]

Morales Chacon, Lilia Maria, Ivan Garcia Maeso, Margarita M. Baez Martin, Juan E. Bender Del Busto, Maria Eugenia Garcia Navarro, Nelson Quintanal Cordero, Barbara Estupinan Diaz, et al. 2018. “Long-Term Electroclinical and Employment Follow up in Temporal Lobe Epilepsy Surgery. A Cuban Comprehensive Epilepsy Surgery Program.” Behavioral Sciences (Basel, Switzerland) 8 (2). https://doi.org/10.3390/bs8020019. [ Links ]

Oishi, Kenichi, Andreia Faria, Hangyi Jiang, Xin Li, Kazi Akhter, Jiangyang Zhang, John T Hsu, et al. 2009. “Atlas-Based Whole Brain White Matter Analysis Using Large Deformation Diffeomorphic Metric Mapping: Application to Normal Elderly and Alzheimer’s Disease Participants.” NeuroImage 46 (2): 486-99. https://doi.org/10.1016/j.neuroimage.2009.01.002. [ Links ]

Owens-Walton, Conor, David Jakabek, Brian D Power, Mark Walterfang, Dennis Velakoulis, Danielle van Westen, Jeffrey C L Looi, Marnie Shaw, and Oskar Hansson. 2019. “Increased Functional Connectivity of Thalamic Subdivisions in Patients with Parkinson’s Disease.” PloS One 14 (9): e0222002-e0222002. https://doi.org/10.1371/journal.pone.0222002. [ Links ]

Pauli, Wolfgang M., Amanda N. Nili, and J. Michael Tyszka. 2018. “A High-Resolution Probabilistic in Vivo Atlas of Human Subcortical Brain Nuclei.” Scientific Data 5 (1): 180063. https://doi.org/10.1038/sdata.2018.63. [ Links ]

Rolls, Edmund T., Chu-Chung Huang, Ching-Po Lin, Jianfeng Feng, and Marc Joliot. 2020. “Automated Anatomical Labelling Atlas 3.” NeuroImage 206 (February): 116189. https://doi.org/10.1016/j.neuroimage.2019.116189. [ Links ]

Sánchez-Catasús, CA, O. Díaz Rizo, M. Rodríguez Castillo, R. Rodríguez Rojas, and A. Aguila Ruiz. 2002. “Mejoramiento de La Cuantificación Relativa Del Flujo Sanguíneo Cerebral Con La Tecnología Médico-Nuclear Disponible En El País.” Nucleus, no. 31: 9-18. [ Links ]

Thomas, Benjamin A., Kjell Erlandsson, Marc Modat, Lennart Thurfjell, Rik Vandenberghe, Sebastien Ourselin, and Brian F. Hutton. 2011. “The Importance of Appropriate Partial Volume Correction for PET Quantification in Alzheimer’s Disease.” European Journal of Nuclear Medicine and Molecular Imaging 38 (6): 1104-19. https://doi.org/10.1007/s00259-011-1745-9. [ Links ]

Zhang, Y., M. Brady, and S. Smith. 2001. “Segmentation of Brain MR Images through a Hidden Markov Random Field Model and the Expectation-Maximization Algorithm.” IEEE Transactions on Medical Imaging 20 (1): 45-57. https://doi.org/10.1109/42.906424. [ Links ]

Contribuciones de los autores

5Adquisición de fondos: -

9Recursos: -

Recibido: 15 de Diciembre de 2021; Aprobado: 17 de Febrero de 2022

*Autor para la correspondencia. (karlabg@infomed.sld.cu)

El autor autoriza la distribución y uso de su artículo.

Conceptualización: Karla Batista García-Ramó.

Curación de datos: Karla Batista García-Ramó, Ángel Águila Ruiz y Nelson Rodríguez Mesa

Análisis formal: Karla Batista García-Ramó, Nancy Pavón Fuentes

Investigación: Karla Batista García-Ramó, Nancy Pavón Fuentes y Lilia Morales Chacón

Metodología: Karla Batista García-Ramó.

Administración del proyecto: Karla Batista García-Ramó.

Software: Karla Batista García-Ramó.

Supervisión: Nancy Pavón Fuentes, Lilia Morales Chacón.

Validación: Karla Batista García-Ramó y Nelson Rodríguez Mesa

Visualización: Karla Batista García-Ramó.

Redacción - borrador original: Karla Batista García-Ramó.

Redacción - revisión y edición: Karla Batista García-Ramó.

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons