Este manuscrito describe un protocolo para la radiografía de neutrones y la tomografía computarizada de muestras biológicas utilizando una línea de luz CG-1D del reactor de isótopos de alto flujo (HFIR) para medir un implante metálico en un fémur de rata, un pulmón de ratón y un sistema de raíces / suelo de plantas herbáceas.
Los neutrones se han utilizado históricamente para una amplia gama de aplicaciones biológicas empleando técnicas como la dispersión de neutrones de ángulo pequeño, el eco de espín de neutrones, la difracción y la dispersión inelástica. A diferencia de las técnicas de dispersión de neutrones que obtienen información en el espacio recíproco, las imágenes de neutrones basadas en la atenuación miden una señal en el espacio real que se resuelve en el orden de decenas de micrómetros. El principio de la imagen de neutrones sigue la ley de Beer-Lambert y se basa en la medición de la atenuación de neutrones a granel a través de una muestra. Una mayor atenuación es exhibida por algunos elementos ligeros (en particular, hidrógeno), que son componentes principales de las muestras biológicas. Los agentes de contraste como el deuterio, el gadolinio o los compuestos de litio se pueden usar para mejorar el contraste de manera similar a como se hace en las imágenes médicas, incluidas técnicas como imágenes ópticas, imágenes por resonancia magnética, rayos X y tomografía por emisión de positrones. Para los sistemas biológicos, la radiografía de neutrones y la tomografía computarizada se han utilizado cada vez más para investigar la complejidad de la red de raíces de plantas subterráneas, su interacción con los suelos y la dinámica del flujo de agua in situ. Además, se han explorado los esfuerzos para comprender los detalles de contraste en muestras de animales, como tejidos blandos y huesos. Este manuscrito se centra en los avances en bioimagen de neutrones, como la preparación de muestras, instrumentación, estrategia de adquisición de datos y análisis de datos utilizando la línea de haz de imágenes de neutrones CG-1D del reactor de isótopos de alto flujo. Las capacidades antes mencionadas se ilustrarán utilizando una selección de ejemplos en fisiología vegetal (planta herbácea / sistema de raíces / suelo) y aplicaciones biomédicas (fémur de rata y pulmón de ratón).
El principio de la radiografía de neutrones (nR) se basa en la atenuación de los neutrones a través de la materia que atraviesan. A diferencia de los rayos X que son dispersados por la nube de electrones de un átomo, los neutrones pueden ser absorbidos o dispersados por su núcleo. Los neutrones son sensibles a los elementos ligeros, como el hidrógeno (H), y en consecuencia pueden utilizarse para aplicaciones biológicas radiográficas como tejidos animales 1,2,3,4,5,6,7 o humanos 8,9 y sistemas subterráneos de suelo/raíces 10,11,12,13,14 ,15. La imagen de neutrones es una técnica complementaria a la imagen de rayos X, que es capaz de detectar elementos pesados16,17,18. El nR basado en la atenuación se rige por los coeficientes de atenuación lineal de los materiales dentro de la muestra y por el espesor de la muestra, como se describe en la ley de Beer-Lambert, que establece que el haz transmitido es directamente proporcional a la cantidad de material y la longitud del camino a través del material. Por lo tanto, la transmitancia, T, se puede calcular como:
(1)
donde I 0 e I son, respectivamente, las intensidades del haz incidente y transmitido; μ y x son el coeficiente de atenuación lineal y el espesor de una muestra homogénea, respectivamente. El coeficiente de atenuación μ viene dado por:
(2)
donde σ es la sección transversal de atenuación de neutrones de la muestra (tanto dispersión como absorción), ρ es su densidad, NA es el número de Avogadro y M es su masa molar.
El contraste en la radiografía de muestras biológicas que utilizan neutrones de baja energía (es decir, energías inferiores a 0,5 eV) se debe principalmente a un cambio en la densidad de H (para un espesor de muestra fijo). Esto se debe a la probabilidad de interacción de un neutrón con el núcleo H, que es mayor que con otros núcleos presentes en muestras biológicas, y al hecho de que la densidad del átomo de H es primordial ya que es el átomo más abundante en muestras biológicas.
Desde sus primeras etapas, nR y la tomografía computarizada de neutrones (nCT) se han utilizado ampliamente para aplicaciones de materiales e ingeniería 19,20,21,22,23. Los primeros experimentos de demostración de sensibilidad neutrónica al H en muestras biológicas comenzaron a mediados de la década de 195024 con las mediciones de especímenes de plantas. El trabajo continuó durante la década de 1960 con, por ejemplo, la radiografía de un tórax humano25 o ratas26, en la que se exploró el uso de agentes de contraste, como el óxido de gadolinio (Gd2O3). Además, se planteó la hipótesis de que el contraste en el tejido tumoral humano frente al tejido normal se debía a un aumento local en el contenido de H. Durante estos ensayos iniciales, se concluyó que el aumento del flujo de neutrones y la resolución espacial mejoraría la calidad de nR y probablemente aumentaría su popularidad como técnica complementaria para aplicaciones industriales o biomédicas. Los estudios más recientes comprenden mediciones de nR y nCT realizadas en muestras de tejido canceroso1 y secciones de órganos animales 2,3,27 para aplicaciones biomédicas y forenses.
Ubicado en el Laboratorio Nacional de Oak Ridge, Oak Ridge, TN, el Reactor de Isótopos de Alto Flujo (HFIR) es una poderosa fuente de neutrones que produce neutrones por reacción de fisión. Estos neutrones tienen energías del orden de 2 MeV y son “enfriados” en la piscina del reactor por reacciones cinéticas con agua pesada para alcanzar energías del orden de 100-300 eV. La optimización de un experimento de neutrones, ya sea dispersión o imagen, comienza con la comprensión de la fuente de neutrones y las propiedades de la línea de luz, como la intensidad del haz, la distribución de energía y el efecto del fondo (neutrones rápidos, neutrones retardados, rayos gamma). En la sala de guía fría HFIR donde se encuentra la línea de luz de imágenes, los neutrones se “enfrían” aún más por interacciones cinéticas con un moderador H líquido. Luego se transportan en un sistema de guía curvo lejos de la línea de visión de la fuente, eliminando así los neutrones rápidos y la contaminación gamma. Como se ilustra en la Figura 1, la línea de haz de imágenes de neutronesCG-1D 28,29 se coloca en una guía fría, lo que implica que el rango de energía de neutrones varía de unos pocos meV a unas pocas decenas de eV (en este caso, la longitud de onda de neutrones utilizable correspondiente varía de 0.8 a 10 Å) con un flujo en el rango de 107 n / (cm2∙s) en la posición de la muestra. Un sistema motorizado de apertura/difusor define la geometría estenopeica del instrumento de imagen. Los neutrones viajan una distancia de 6,59 m en un tubo de vuelo lleno de helio (He) con ventanas de aluminio (Al) en cada extremo. Los tubos de vuelo se utilizan para transportar neutrones al tiempo que limitan la dispersión del aire de tal manera que la pérdida en la intensidad del haz sea mínima. Para las mediciones descritas en este manuscrito, el difusor está hecho de un nanopolvo de óxido de aluminio de 50 nm (Al2O3) de 1 mm de espesor encerrado en un recipiente de Al. El difusor reduce los artefactos del haz provenientes de la guía de neutrones (que se magnifican por la geometría estenopeica de una línea de luz de imagen), de lo contrario, las fluctuaciones de intensidad horizontales y verticales agudas son visibles en la radiografía y la normalización de los datos se vuelve un desafío. Para los experimentos ilustrados aquí, los neutrones se convierten en luz utilizando un fluoruro de litio-6 de 25 μm de espesor / fósforo de sulfuro de zinc (6LiF / ZnS: Ag).
La optimización de la colimación depende de la posición de la muestra al detector, la resolución espacial requerida y el tiempo de adquisición. Cuando la muestra se encuentra a pocos cm del centelleador, las colimaciones altas (L / D por encima de 800, donde L es la distancia desde la apertura del agujero de alfiler del diámetro, D y el detector) producen una mejor resolución espacial a costa del flujo de neutrones. La baja colimación (L/D por debajo de 800) es preferible para estudios dinámicos in situ cuando la resolución temporal prevalece sobre la resolución espacial. Para las mediciones descritas en este manuscrito, L/D y resolución espacial fueron aproximadamente 355 y 75 μm, respectivamente. La resolución temporal varió en función de la relación señal-ruido (SNR). La muestra se colocó lo más cerca posible del centelleador para reducir la distorsión geométrica, como el desenfoque. Las etapas de traslación y rotación están disponibles para colocar la muestra cerca de los detectores y realizar la tomografía computarizada (TC). CG-1D ofrece tres tipos de detectores: un dispositivo de carga acoplada (CCD) con 2048 píxeles x 2048 píxeles con un paso de píxeles de 13,5 μm, un detector científico complementario de semiconductores de óxido metálico (sCMOS) con 2560 píxeles x 2160 píxeles con un paso de píxeles de 6,5 μm, y un detector de placa de microcanal (MCP)30,31 con 512 píxeles x 512 píxeles con un tamaño de píxel de 55 μm. Los neutrones dispersos se absorben con goma de boro de ~ 5 mm de espesor para proteger el chip detector de ver neutrones. Esta absorción genera rayos gamma que pueden ser detenidos por el plomo (Pb) colocado entre la goma de boro y el detector. Cada detector está optimizado para un campo de visión (FOV) diferente, así como para resoluciones espaciales y temporales. Para las mediciones del fémur de rata y el pulmón del ratón, se utilizó el detector CCD por su gran capacidad FOV (~ 7 cm x 7 cm) y una resolución espacial razonable de aproximadamente 75 μm. La nCT del sistema radicular/suelo de la planta se realizó con el sCMOS, ya que el objetivo era adquirir nCT lo más rápido posible a costa del campo de visión (que se limitó a ~ 5 cm x 4,2 cm); Por lo tanto, la resolución espacial evidentemente sufrió. En estos detectores, los neutrones se convierten en luz o en una partícula alfa para fines de detección. La rotación de la muestra alrededor de su eje vertical y la adquisición de radiografías en ángulos de rotación consecutivos permite la adquisición de nCT. El modelo de renderizado volumétrico tridimensional de la muestra bajo investigación se obtiene utilizando el portátil interno de retroproyección filtrada (FBP) Jupyter basado en python iMARS3D, pyMBIR o un software comercial, todo lo cual se describe a continuación.
Finalmente, los neutrones que no han interactuado con la muestra o el detector se recogen en una posición de parada del haz aproximadamente 1 m aguas abajo del sistema detector para minimizar el ruido de fondo. El tope de viga CG-1D tiene 0,75 m de ancho, 0,5 m de alto y 35 mm de grosor y está hecho de B4C en epoxi. El tope del haz está reforzado con 10 mm de carbonato de litio enriquecido al 95% (6 Li2CO3) en un epoxi resistente al fuego donde el haz de neutrones golpea, con una cavidad revestida con 6Li, plomo (Pb) y acero diseñado para contener la alta tasa de rayos gamma secundarios. El tope de la viga está directamente unido a la pared de blindaje de acero de la línea de luz. Una fotografía de la línea de luz CG-1D se da en la Figura 2.
Se utilizaron tres programas de reconstrucción para reconstruir los tres datos experimentales en 3D, respectivamente. La reconstrucción de la muestra de pulmón de ratón se realizó utilizando Octopus32, un software de reconstrucción comercial que utiliza FBP. El software Octopus se encuentra en una PC servidor y se puede utilizar para reconstruir los datos recopilados en la línea de luz. Un software de reconstrucción, llamado iMARS3D, está disponible en CG-1D. Se basa en el código fuente abierto TomoPY33 con características añadidas como corrección automática de inclinación, filtros de post-procesamiento, etc. iMARS3D incluye el preprocesamiento de los datos (resta del fondo y el ruido), el recorte, el filtrado de medianas (para corregir los rayos gamma y los píxeles muertos), la corrección automatizada de la fluctuación de la intensidad del haz y la corrección de la inclinación de la muestra. Una vez que se crean los sinogramas, el procesamiento adicional de datos, como la eliminación de artefactos de anillo y el suavizado, son una opción. Los diferentes pasos de la reconstrucción se guardan en el servidor de análisis (y luego se mueven en la carpeta compartida de la propuesta), mientras que los segmentos 2D finales se almacenan inmediatamente en la carpeta compartida de la propuesta. El fémur de la rata fue reconstruido usando iMARS3D. La muestra de raíz / suelo de la planta se procesó previamente mediante el filtrado mediano de los datos utilizando TomoPY seguido de la corrección del eje de inclinación utilizando la biblioteca SciPy de Python. La reconstrucción se llevó a cabo utilizando un paquete de Python desarrollado internamente denominado pyMBIR (construido con kernels de la caja de herramientas ASTRA34) que implementa un conjunto de algoritmos tomográficos desde la línea de base FBP hasta técnicas avanzadas de reconstrucción iterativa basadas en modelos35 que pueden obtener reconstrucciones de alta calidad a partir de conjuntos de datos de neutrones extremadamente dispersos y ruidosos. Todos los volúmenes renderizados basados en las herramientas de reconstrucción mencionadas anteriormente se representan en contraste de atenuación. Toda la visualización se realizó utilizando el paquete de software comercial de visualización, segmentación y análisis de datos AMIRA36.
Este manuscrito tiene como objetivo demostrar el procedimiento de uso de imágenes de neutrones (nR y nCT) en la línea de luz HFIR CG-1D. Este estudio también ilustra las capacidades actuales de nR y nCT de vanguardia para muestras biológicas, específicamente un pulmón de ratón, un hueso de rata y sistemas de raíces / suelo de plantas. El pulmón de ratón fue elegido para ilustrar la complementariedad de los neutrones para medir el tejido pulmonar, mientras que los rayos X son principalmente sensibles a los huesos. La muestra ósea, un fémur de rata, tenía un implante de titanio (Ti), lo que ilustra el contraste entre el hueso y el metal, y la oportunidad de ver la interfaz hueso/metal (que es difícil de medir con rayos X ya que los metales los atenúan fuertemente4). Finalmente, el sistema de agua planta-raíz ilustra la capacidad tridimensional (3D) de nCT para medir los sistemas de raíces / suelo in situ. Además, muestra las ventajas / desventajas de usar nR para muestras biológicas. Evidentemente, este método puede usarse con seguridad para medir la dinámica del agua en un sistema planta-raíz, pero no puede considerarse como una técnica de imagen animal o humana viva debido a los riesgos asociados con la exposición a la radiación, limitando así los estudios a ratones (muertos) o mediciones similares a la patología en las que, por ejemplo, una muestra de tejido se reseca de un paciente (animal o humano) y se prepara por fijación antes de medirse en un haz de neutrones.
La radiografía de neutrones y la TC de muestras biológicas son técnicas de imagen prometedoras que complementan las imágenes de rayos X o las imágenes por resonancia magnética. Los pasos críticos para realizar un experimento de imágenes de neutrones de una muestra biológica están relacionados con su preparación y su contención en la línea de luz. La optimización de un experimento es impulsada por la pregunta científica a responder. Si la pregunta científica requiere una alta resolución espacial para observar un fenómeno, entonces se requieren largos tiempos de adquisición, y el inconveniente de nCT (con campo de visión de tamaño cm) es que lleva horas realizar un escaneo. Esto se debe principalmente a la diferencia en el flujo total de neutrones disponible en un reactor en comparación con una fuente de sincrotrón, donde las tomografías computarizadas de rayos X pueden tomar de segundos a minutos para un campo de visión de unos pocos mm2 . Aunque el método puede aplicarse a muestras de tejido ex vivo extraídas de animales, no puede extenderse in vivo a animales vivos o humanos debido al riesgo de exposición a la radiación (como los rayos gamma producidos por neutrones y las interacciones de neutrones con los átomos de la muestra). Sin embargo, es muy adecuado para la obtención de imágenes de las interacciones entre la raíz y el suelo de la planta (Figura 7), como la dinámica de absorción de agua.
La ventaja de usar nCT rápido para la dinámica de la planta es la sensibilidad al H en el agua y la ausencia de daño por radiación a la planta, a diferencia de la TC de rayos X. Además, se puede observar un contraste único del uso de neutrones en muestras de hueso / metal como un fémur de rata donde el metal es relativamente transparente en comparación con los tejidos circundantes (Figura 5), evitando potencialmente los artefactos metálicos inducidos por la TC39 de rayos X. Los tejidos animales, como el pulmón de ratón (Figura 6), muestran una detección impresionante de la estructura de los tejidos blandos porque los neutrones son sensibles al H, pero la resolución espacial es un factor limitante en estas mediciones. El contraste es proporcionado por los átomos de H presentes en muestras biológicas19,39.
Con los avances de técnicas novedosas como la interferometría de rejilla de neutrones y la mejora en la resolución espacial (recientemente se han reportado algunas micras42,43), las imágenes de neutrones pueden ofrecer nuevos mecanismos de contraste para los tejidos biológicos con una resolución espacial mejorada. La exploración de neutrones de mayor energía (para permitir las mediciones de muestras gruesas) también promete la capacidad de medir secciones más grandes de un tejido animal, como un ratón intacto, ofreciendo así nuevas posibilidades para la investigación biomédica.
The authors have nothing to disclose.
Parte de esta investigación utilizó recursos en el Reactor de Isótopos de Alto Flujo, operado por ORNL, y patrocinado por el Departamento de Energía de los Estados Unidos, Oficina de Ciencia, Instalaciones de Usuario, bajo el contrato DE-AC05-00OR22725 con UT-Battelle, LLC. Parte de esta investigación fue apoyada por ORNL a través del programa Eugene Wigner Distinguished Staff Fellowship. Esta investigación también fue patrocinada por la Oficina de Ciencia del DOE, Oficina de Investigación Biológica y Ambiental. Las muestras femorales de rata se obtuvieron de experimentos realizados en colaboración con el Dr. Rick Sumner en el Centro Médico de la Universidad de Rush con fondos obtenidos del NIH (R01AR066562) y del premio Smith and Nephew de la Fundación de Investigación y Educación Ortopédica. El equipo quiere agradecer a los equipos de apoyo HFIR que permiten el uso de las líneas de haz de dispersión de neutrones.
Aluminum containers | custom | Made from aluminum plates or tubing (alternate is quartz), plant and mouse sample | |
Aluminum foil | Fisher | 01-213-100 | Mouse lung sample containment |
Deionized water or deuterium oxide | Water or D2O can be used to enhance contrast, plant sample | ||
Ethanol | Fisher | 04-355-223 | Mouse lung sample |
Gauze sponges | CardinalHealth | Fully submerged in phosphate-buffered saline (PBS) and used to wrap samples, rat femur sample | |
Growth chamber | Conviron | A1000 | Any growth chamber or greenhouse with controlled conditions would work, plant sample |
Laboratory balance | Weighing plant system can be used to measure actual water content in the soils, plant sample | ||
Pure silica sand | US Silica Co. | Flint#13 | Pure SiO2 provides low neutron attenuation compared to soils, plant sample |
Sprague-Dawley Rats | Harlan | Order Code: 002-US | Rat femur sample |
Titanium Rod | Goodfellow | TI007905 | Rat femur sample |